• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于抽樣調(diào)查資料估算區(qū)域土壤侵蝕量

    2017-04-18 07:09:02王振穎何遠(yuǎn)梅
    中國(guó)水土保持科學(xué) 2017年1期
    關(guān)鍵詞:坡長(zhǎng)模數(shù)土壤侵蝕

    王振穎,何遠(yuǎn)梅

    (1.遼寧省水利水電科學(xué)院,110003,沈陽(yáng);2.北京林業(yè)大學(xué)水土保持學(xué)院,100083,北京)

    基于抽樣調(diào)查資料估算區(qū)域土壤侵蝕量

    王振穎1,何遠(yuǎn)梅2

    (1.遼寧省水利水電科學(xué)院,110003,沈陽(yáng);2.北京林業(yè)大學(xué)水土保持學(xué)院,100083,北京)

    區(qū)域土壤侵蝕量估算是土壤侵蝕調(diào)查的重點(diǎn)和難點(diǎn),為了利用抽樣調(diào)查數(shù)據(jù)定量估算區(qū)域土壤侵蝕狀況,以2011年第1次全國(guó)水利普查中的水力侵蝕抽樣調(diào)查資料為基礎(chǔ),利用中國(guó)土壤流失方程(CSLE),通過(guò)地理信息系統(tǒng)估算2011年遼寧省的降雨侵蝕力、土壤可蝕性、地形、水土保持措施等侵蝕因子及侵蝕量,在此基礎(chǔ)上對(duì)全省的水力侵蝕強(qiáng)度進(jìn)行分級(jí)和面積統(tǒng)計(jì),并與抽樣調(diào)查結(jié)果進(jìn)行對(duì)比。結(jié)果顯示:全省各土壤侵蝕因子的估算結(jié)果與抽樣調(diào)查結(jié)果較為接近;土壤水蝕以微度和輕度侵蝕為主,每年全省約85%的土地處于微度侵蝕(<1 000 t/km2),而強(qiáng)度侵蝕(>5 000 t/km2)面積僅占2%;抽樣調(diào)查得到的獨(dú)立工礦用地、果園、旱地和其他土地的侵蝕模數(shù)較估算得出的侵蝕模數(shù)大,而在草地、其他林地等土地利用類型上抽樣調(diào)查得到的侵蝕模數(shù)卻較估算得出的侵蝕模數(shù)小,主要原因可能是省域和抽樣單元計(jì)算地形因子的DEM數(shù)據(jù)精度不同??傮w上,該研究為如何利用抽樣調(diào)查數(shù)據(jù)估算全國(guó)土壤侵蝕狀況提供參考。

    抽樣調(diào)查資料;中國(guó)土壤流失方程(CSLE);GIS應(yīng)用;區(qū)域土壤侵蝕;遼寧省

    在較長(zhǎng)時(shí)段和較大空間尺度上的區(qū)域土壤侵蝕研究多依賴于土壤侵蝕試驗(yàn)與GIS技術(shù),其研究成果能較好地為水土保持宏觀決策提供支持,且與全球變化研究密切相關(guān)[1-4]。我國(guó)已經(jīng)開展了3次全國(guó)性的土壤侵蝕調(diào)查[5]:20世紀(jì)50年代,采用人工調(diào)查法,根據(jù)植被蓋度與坡度,在地形圖上進(jìn)行侵蝕強(qiáng)度的等級(jí)劃分,該法耗時(shí)耗力,且無(wú)法定量評(píng)價(jià)土壤侵蝕量;20世紀(jì)80年代和90年代末期,在水利部組織下進(jìn)行的2次全國(guó)土壤侵蝕遙感調(diào)查與公告,是基于遙感影像資料和GIS技術(shù)的大尺度調(diào)查,主要以植被覆蓋度和坡度為依據(jù)確定土壤侵蝕類型和強(qiáng)度[6],而并未考慮降雨、土壤可蝕性、坡長(zhǎng)、水土保持措施等重要因子。長(zhǎng)期以來(lái),土壤侵蝕動(dòng)態(tài)監(jiān)測(cè)、定量評(píng)價(jià)等方面的研究大多是在坡面和小流域等較小尺度上進(jìn)行的,對(duì)于區(qū)域土壤侵蝕的定量評(píng)價(jià)雖有研究[7-10],但尚未提出比較全面而實(shí)用的尺度轉(zhuǎn)換方法。在區(qū)域土壤侵蝕調(diào)查方法方面,美國(guó)農(nóng)業(yè)部提出的土壤侵蝕抽樣調(diào)查方法是最具有特色的,該法基于統(tǒng)計(jì)學(xué)原理在全國(guó)布設(shè)樣點(diǎn),利用USLE計(jì)算土壤流失量,再統(tǒng)計(jì)得到全國(guó)各類土壤侵蝕面積及土壤侵蝕量[11]。

    楊勤科等[12]對(duì)此種抽樣普查方法進(jìn)行了初步討論,建議采用我國(guó)土壤侵蝕最新模型,并結(jié)合計(jì)算機(jī)模擬實(shí)現(xiàn)土壤侵蝕強(qiáng)度的定量估算,同時(shí)指出調(diào)查內(nèi)容應(yīng)包括主要的侵蝕影響因子;李智廣等[13]闡述了我國(guó)水力侵蝕抽樣調(diào)查方法;張巖等[14]為了測(cè)試抽樣方法和土壤侵蝕模型在土壤侵蝕普查中的適用性,以陜西吳起縣為試點(diǎn),對(duì)均勻抽樣方法和基于遙感數(shù)據(jù)的水蝕分級(jí)分類方法進(jìn)行比較。這些研究為第1次全國(guó)水利普查提供了理論依據(jù)。在第1次全國(guó)水利普查水土保持情況普查中,雖然采用抽樣調(diào)查法獲取抽樣調(diào)查單元的實(shí)地?cái)?shù)據(jù),但如何快速準(zhǔn)確呈現(xiàn)區(qū)域土壤水力侵蝕的空間分布狀況與水力侵蝕量仍是亟待解決的問題。本文基于GIS技術(shù),利用抽樣調(diào)查資料和中國(guó)土壤流失方程(CSLE),估算遼寧省的降雨侵蝕力、土壤可蝕性、地形、水土保持措施等因子及土壤侵蝕量,在此基礎(chǔ)上對(duì)全省土壤侵蝕強(qiáng)度進(jìn)行分級(jí)統(tǒng)計(jì),并將計(jì)算結(jié)果與抽樣調(diào)查結(jié)果進(jìn)行對(duì)比,該研究可為探索抽樣調(diào)查數(shù)據(jù)估算區(qū)域土壤侵蝕狀況的方法提供參考。

    1 研究區(qū)概況

    遼寧省位于我國(guó)東北地區(qū)南部,含14個(gè)地級(jí)市,56個(gè)縣,地理坐標(biāo)為E 118°50'~125°47',N 38° 43'~43°29'。總面積30.61萬(wàn)km2,其中陸地面積14.59萬(wàn)km2,海域面積15.02萬(wàn)km2。遼東為海拔800 m左右的山地,遼西為海拔500 m以下的丘陵,山地和丘陵約占全省總面積的60%。遼寧省中部為海拔200 m以下的平原,約占全省總面積的33%,水面及其他約占7%。遼寧省屬溫帶大陸性季風(fēng)氣候,從東南向西北逐漸由濕潤(rùn)氣候向半濕潤(rùn)、半干旱氣候過(guò)渡。年平均氣溫在5.2~11.7℃之間,年平均降水量為500~1 000 mm,且由東向西逐漸減少。遼寧省的土壤地帶性分布明顯,從南至北呈現(xiàn)暖溫落葉闊葉林棕壤—中溫帶森林草甸黑土—草甸草原黑鈣土的緯度地帶性,從東至西呈現(xiàn)遼東山地丘陵為棕壤—下遼河平原為草甸土—遼西山地丘陵為褐土的經(jīng)度地帶性。地帶性森林植被主要為紅松(Pinus koraiensis)、沙松(Abies holophyllaMaxim)、闊葉混交林和暖溫帶落葉闊葉林。

    2 材料和方法

    2.1 數(shù)據(jù)來(lái)源

    抽樣調(diào)查數(shù)據(jù)來(lái)源于2011年第1次全國(guó)水利普查中水土保持情況普查的水力侵蝕抽樣調(diào)查。遼寧省水蝕調(diào)查抽樣單元共1 287個(gè),其中1 235個(gè)調(diào)查單元是按照1%的密度布點(diǎn),另外52個(gè)點(diǎn)(分布在4個(gè)縣)按照0.25%的密度布點(diǎn)。基礎(chǔ)地理與專題地圖數(shù)據(jù)主要包括DEM、行政區(qū)劃圖、流域分界圖、土地利用圖。降水量數(shù)據(jù)為遼寧56個(gè)縣1981—2010年逐日降水量≥12 mm的資料。

    2.2 土壤流失量估算方法

    本研究所采用的土壤流失方程CSLE[15]計(jì)算公式如下:

    式中:A為多年平均土壤流失量,t/hm;R為降雨侵蝕力因子,MJ·mm/(hm2·h);K為土壤可蝕性因子, t·hm2·h/(hm2·MJ·mm);S為坡度因子;L為坡長(zhǎng)因子;B為水土保持生物措施因子;E為水土保持工程措施因子;T為水土保持耕作措施因子。

    2.2.1 降雨侵蝕力采用章文波等[16]提出的方法,利用日降雨量直接估算逐年半月降雨侵蝕力的簡(jiǎn)易計(jì)算模型計(jì)算降雨侵蝕力,得出各雨量點(diǎn)的多年平均降雨侵蝕力值后在ArcGIS下用Kriging內(nèi)插法進(jìn)行插值,從而得到全省的降雨侵蝕力空間分布圖。

    2.2.2 土壤可蝕性因子在一定區(qū)域內(nèi)土壤類型與土地利用方式密切相關(guān),所以土壤可蝕性與土地利用也密切相關(guān)。本研究基于1∶10萬(wàn)土地利用圖和抽樣調(diào)查結(jié)果推算全省的土壤可蝕性,步驟如下:1)計(jì)算縣級(jí)各土地利用條件下抽樣調(diào)查單元的平均土壤可蝕性因子;2)使用空間關(guān)聯(lián)方法對(duì)各縣不同土地利用斑塊賦予相應(yīng)的土壤可蝕性值,從而得到全省的土壤可蝕性因子空間分布圖。

    2.2.3 坡度和坡長(zhǎng)因子土壤侵蝕普查成果中已經(jīng)基于1∶1萬(wàn)地形圖計(jì)算了各抽樣單元中各地塊的坡度和坡長(zhǎng)因子,但本研究中只有1∶5萬(wàn)地形圖生成的DEM,在進(jìn)行坡度和坡長(zhǎng)因子估算時(shí):1)根據(jù)1∶5萬(wàn)地形圖生成的DEM,計(jì)算并提取與抽樣單元相對(duì)應(yīng)地塊的S和L值;2)比較基于2種比例尺(1∶1萬(wàn)和1∶5萬(wàn))地形圖計(jì)算的S和L值,發(fā)現(xiàn)雖然兩者獲取的地形因子差異較大,但在特定區(qū)域仍具有較好的線性相關(guān)性;3)對(duì)于坡長(zhǎng)因子和坡度因子>4(坡度≥25°)的部分,分別以地級(jí)市為單位建立2種土地利用類型(大農(nóng)業(yè)用地、居民點(diǎn)及建設(shè)用地)基于1∶1萬(wàn)地形圖與基于1∶5萬(wàn)地形圖計(jì)算的地形因子的線性回歸方程,提取回歸方程的常數(shù)項(xiàng)和斜率作為校正參數(shù)圖層;4)把參數(shù)圖層和基于1∶5萬(wàn)地形圖計(jì)算的全省S和L圖層進(jìn)行運(yùn)算,得到校正后的坡度和坡長(zhǎng)因子圖。

    2.2.4 生物措施因子、工程措施因子和耕作措施因子本研究根據(jù)全省1∶10萬(wàn)土地利用現(xiàn)狀圖和抽樣調(diào)查結(jié)果推算B、E、T值。具體方法如下:1)以縣為單位統(tǒng)計(jì)抽樣調(diào)查單元中各類土地利用類型的B、E、T均值;2)將B、E、T值聯(lián)結(jié)到土地利用類型圖的屬性表中;3)將土地利用類型圖進(jìn)行數(shù)據(jù)格式轉(zhuǎn)換,分別以B、E、T值作為字段,轉(zhuǎn)化成柵格大小為10 m×10 m的grid文件。

    3 結(jié)果與分析

    3.1 土壤侵蝕影響因子估算

    3.1.1 降雨侵蝕力估算利用遼寧省抽樣單元與全省降雨侵蝕力的統(tǒng)計(jì)數(shù)據(jù),繪制遼寧省抽樣單元與全省降雨侵蝕力頻率分布對(duì)比圖(圖1)。

    圖1 遼寧省抽樣單元與全省降雨侵蝕力頻率分布對(duì)比圖Fig.1 Frequency distribution of rainfall erosivity in the whole Liaoning Province and the sampling plots

    從圖2可以看出:遼寧省的降雨侵蝕力主要分布在1 500~3 500 MJ·mm/(hm2·h)范圍內(nèi),并且1 500~3 500 MJ·mm/(hm2·h)的侵蝕力在抽樣調(diào)查和全省估算結(jié)果中所占比例分別為84.026%和75.551%。在1 500~3 500 MJ·mm/(hm2·h)的侵蝕力范圍中又以2 500~3 000 MJ·mm/(hm2·h)的比例為最大,占抽樣調(diào)查和全省估算結(jié)果的比例分別為36.038%和27.237%。通過(guò)對(duì)比可以看出2次估算的降雨侵蝕力差異不大。

    3.1.2 土壤可蝕性因子估算由抽樣單元與全省各土地利用類型K均值的比較結(jié)果(圖2)可知,抽樣調(diào)查和估算值的差異范圍為0~0.001 5,相對(duì)于抽樣單元的相對(duì)誤差均小于10%。全省估算值與抽樣調(diào)查的K均值誤差較大的是居民點(diǎn)、商服用地、獨(dú)立工礦用地及草地,但是這4種土地利用類型所占面積比例很小??梢姳狙芯抗浪愕娜「魍恋乩妙愋偷腒均值與抽樣單元K均值非常接近。

    3.1.3 地形因子估算根據(jù)1∶1萬(wàn)地形圖獲取抽樣單元中的地形因子,對(duì)基于1∶5萬(wàn)地形圖獲取的地形因子進(jìn)行線性校正。從圖3可以看出,經(jīng)過(guò)校正后坡度因子大于6的頻率大幅減少,校正后頻率分布情況與抽樣單元的坡度因子分布情況更為接近。

    抽樣單元根據(jù)1∶1萬(wàn)地形圖計(jì)算的坡長(zhǎng)因子都在2.6以內(nèi),基于1∶5萬(wàn)地形圖獲取的坡長(zhǎng)因子也都小于3。從圖3可以看出抽樣單元的坡長(zhǎng)因子接近于正態(tài)分布,主要集中在0.5到2,占85%左右。全省坡長(zhǎng)因子經(jīng)過(guò)調(diào)整,大于2的減少,1~2之間的增加,0~0.5之間的減少,與抽樣單元的坡長(zhǎng)因子分布情況趨于一致,但在0~1之間的調(diào)整結(jié)果不理想。

    圖2 遼寧省抽樣單元與全省在不同土地利用類型下平均土壤可蝕性K值對(duì)比圖Fig.2 Values of soil erodibility K under different land use in the whole Liaoning Province and the sampling plots

    圖3 遼寧省坡度和坡長(zhǎng)因子校正前后與抽樣單元坡度坡長(zhǎng)因子頻率分布對(duì)比圖Fig.3 Frequency distribution of gradient and slope length factor in the whole Liaoning Province and the sampling plots before and after revising

    3.1.4 生物措施因子B、工程措施因子E、耕作措施因子T的估算根據(jù)1 287個(gè)抽樣單元的調(diào)查結(jié)果,將遼寧省18類土地利用類型合并為8大類,分別為:耕地(水田、水澆地、旱地)、園地(果園、其他園地)、林地(有林地、灌木林地、其他林地)、草地(天然牧草地、人工牧草地、其他草地)、商服住宅公共用地及工礦用地(城鎮(zhèn)居民點(diǎn)、農(nóng)村居民點(diǎn)、獨(dú)立工礦用地、商服及公共用地)、交通運(yùn)輸用地、水域及水利設(shè)施用地、其他土地。遼寧省總面積是1 530.19萬(wàn)hm2,抽樣調(diào)查單元的總面積是9.794萬(wàn)hm2。從圖4可知,抽樣單元中各土地利用類型的面積比例與全省各土地利用類型的面積比例差異不大,因此抽樣調(diào)查的水土保持措施因子能很好地反映各土地利用類型的水土保持措施因子。

    圖4 遼寧省抽樣單元與全省不同土地利用類型所占比例對(duì)比圖Fig.4 Proportions of different land use types in the whole Liaoning Province and the sampling plots

    由圖5可知,抽樣調(diào)查和估算的全省水田、水澆地、旱地、獨(dú)立工礦用地的B均值都在0.9以上,其他土地的B均值在0.6左右,林地、草地和果園等的B均值都在0.2以下。僅通過(guò)生物措施因子來(lái)看,耕地和獨(dú)立工礦用地的水土保持效益最差。除其他土地外,估算的全省各土地利用類型的B均值與抽樣調(diào)查單元的B均值雖有差異,但相對(duì)誤差均小于10%,說(shuō)明全省生物措施因子估算方法合理。水域及水利設(shè)施用地、其他土地和水澆地等的E均值接近于1,說(shuō)明這些土地利用類型少有水土保持工程措施,而交通用地、工礦用地及商服用地等的E均值都小于0.3。耕地及人工牧草地的T值較小,均<0.5,其余土地利用類型的T值接近或等于1。抽樣單元的E值和T值也與全省估算值較為一致,說(shuō)明工程和耕作措施因子估算方法也是可行的。

    3.2 土壤侵蝕強(qiáng)度比較

    3.2.1 各級(jí)土壤侵蝕強(qiáng)度面積比較由遼寧省抽樣單元和省域土壤侵蝕強(qiáng)度匯總表(表1)可知,遼寧省土壤水蝕以微度和輕度為主,抽樣調(diào)查結(jié)果和估算結(jié)果都表明處于微度侵蝕的地區(qū)約占85%,強(qiáng)度侵蝕以上的地區(qū)所占面積不到2%。

    圖5 遼寧省抽樣單元與全省各土地利用的B值、E值和T值對(duì)比圖Fig.5 Values of biological measure factor B,engineering measure factor E,and tillage measure factor T under different land use in the whole Liaoning Province and the sampling plots

    3.2.2 不同土地利用類型土壤侵蝕模數(shù)比較對(duì)比遼寧省抽樣調(diào)查和估算的不同土地利用條件下土壤水蝕模數(shù)(表2)可知,土壤侵蝕模數(shù)最大的是其他土地,但其他土地在全省的面積比例很小,可忽略,其次是旱地和草地,各種林地的侵蝕量居中,居民點(diǎn)、交通用地和水域、水田等土壤侵蝕較小。值得注意的是工礦用地的平均土壤侵蝕強(qiáng)度雖然小,但遼寧省作為工礦業(yè)重省,其工礦用地的總侵蝕量不容忽視。

    從表2可以看出,抽樣調(diào)查得到土壤侵蝕模數(shù)和估算的土壤侵蝕模數(shù)之間存在一定差異,主要表現(xiàn)為抽樣調(diào)查得到的獨(dú)立工礦用地、果園、旱地、其他土地的侵蝕模數(shù)與估算得到的侵蝕模數(shù)相比較大,而在草地、其他林地等土地利用類型上抽樣調(diào)查得到的侵蝕模數(shù)卻較小,主要原因可能是估算時(shí)采用的坡度、坡長(zhǎng)因子的地形數(shù)據(jù)精度與抽樣調(diào)查時(shí)不同。

    表1 2011年遼寧省抽樣單元和省域土壤侵蝕強(qiáng)度比較Tab.1 Soil erosion intensities in the whole Liaoning Province and in the sampling plots in 2011

    表2 2011年遼寧省不同土地利用土壤水蝕模數(shù)比較Tab.2 Soil erosion modulus under different land use in the Liaoning Province in 2011t/km2

    3.3 侵蝕估算結(jié)果差異分析

    抽樣調(diào)查結(jié)果和全省估算結(jié)果都說(shuō)明遼寧省土壤水蝕以微度和輕度為主(表1),但不同土地利用類型的平均土壤侵蝕模數(shù)存在較大差異(表2),其原因主要有2方面:1)抽樣調(diào)查單元在實(shí)地勾繪地塊邊界的基礎(chǔ)上,較全面地記錄了土地利用、水土保持措施及土壤方面的信息,對(duì)于區(qū)域具有很好的代表性,但是2011年遼寧省抽樣調(diào)查單元的面積只占全省面積的0.64%,存在差異是正常的;2)省域估算和抽樣單元計(jì)算坡度、坡長(zhǎng)因子的地形數(shù)據(jù)精度不同。

    4 討論

    本文根據(jù)土壤侵蝕抽樣調(diào)查成果推求了全省范圍內(nèi)的土壤侵蝕模數(shù),其關(guān)鍵在于土壤侵蝕影響因子的估算。降雨侵蝕力因子的估算使用與抽樣調(diào)查完全相同的數(shù)據(jù),因此降雨侵蝕力因子的精度不存在問題。土壤可蝕性因子以及生物、工程和耕作措施因子是基于抽樣調(diào)查成果和土地利用類型圖估算得到的,由于土地利用類型圖中缺少工程措施和耕作措施因子信息,只能根據(jù)抽樣調(diào)查單元中各種土地利用條件下的水土保持措施因子均值推算區(qū)域水土保持措施因子,所以對(duì)估算精度有一定影響。

    計(jì)算省域土壤侵蝕模數(shù)最大的問題是地形因子的估算,前人的研究和本研究的分析都表明基于1∶1萬(wàn)地形圖和1∶5萬(wàn)地形圖計(jì)算的坡度和坡長(zhǎng)因子差異較大。雖然本研究以地級(jí)市為單元對(duì)1∶5萬(wàn)地形圖獲取的坡度和坡長(zhǎng)因子進(jìn)行了校正,但地形因子的計(jì)算精度不高,會(huì)直接影響省域土壤侵蝕模數(shù)的估算精度。

    通過(guò)對(duì)比全省和抽樣調(diào)查單元的土壤侵蝕因子值及各級(jí)土壤侵蝕強(qiáng)度面積,發(fā)現(xiàn)本研究方法合理,且估算精度較高,若能獲取更大比例尺的土壤類型圖和地形圖,省域土壤侵蝕的估算精度可以得到進(jìn)一步提升。

    5 結(jié)論

    1)除坡度與坡長(zhǎng)因子外,遼寧全省土壤侵蝕因子的估算結(jié)果與抽樣調(diào)查結(jié)果較為接近。

    2)土壤侵蝕強(qiáng)度估算結(jié)果和抽樣調(diào)查結(jié)果均表明:遼寧省土壤水蝕以微度侵蝕和輕度侵蝕為主,全省約85%的土地處于微度侵蝕,而強(qiáng)度侵蝕以上面積僅占2%。

    3)抽樣調(diào)查得到的獨(dú)立工礦用地、果園、旱地、其他土地的侵蝕模數(shù)比估算得出的侵蝕模數(shù)大,而在草地、其他林地等土地利用類型上抽樣調(diào)查得到的侵蝕模數(shù)卻比估算得出的侵蝕模數(shù)小,主要原因可能是省域和抽樣單元計(jì)算地形因子的DEM數(shù)據(jù)精度不同。

    4)基于抽樣調(diào)查資料和侵蝕模型估算的區(qū)域土壤侵蝕模數(shù)和強(qiáng)度具有很高的可信度,此法為實(shí)現(xiàn)定期和持續(xù)的區(qū)域土壤侵蝕估算提供參考。

    [1] 李銳,楊勤科.區(qū)域水土流失快速調(diào)查與管理信息系統(tǒng)研究[M].鄭州:黃河水利出版社,2000:1. LI Rui,YANG Qinke.The quick survey of regional water and soil loss and management information system[M]. Zhengzhou:The Yellow River Water Conservancy Press, 2000:1.

    [2] 冷疏影,馮仁國(guó),李銳,等.土壤侵蝕與水土保持科學(xué)重點(diǎn)研究領(lǐng)域與問題[J].水土保持學(xué)報(bào),2004,18 (1):1. LENG Suying,FENG Renguo,LI Rui,et al.Key research of soil erosion and conservation in China[J].Journal of Soil and Water Conservation,2004,18(1):1.

    [3] 王禮先,張有實(shí),李銳,等.關(guān)于我國(guó)水土保持科學(xué)技術(shù)的重點(diǎn)研究領(lǐng)域[J].中國(guó)水土保持科學(xué),2005,3 (1):1. WANG Lixian,ZHANG Youshi,LI Rui,et al.On key research domain of science and technology for soil and water conservation in China[J].Science of Soil and Water Conser1ation,2005,3(1):1.

    [4] 楊勤科,李銳,曹明明.區(qū)域土壤侵蝕定量研究的國(guó)內(nèi)外進(jìn)展[J].地球科學(xué)進(jìn)展,2006,21(8):849. YANG Qinke,LI Rui,CAO Mingming.Advances of quantitative assessment on regional soil erosion[J].Advances in Earth Science,2006,21(8):849.

    [5] 郭索彥,李智廣.我國(guó)水土保持監(jiān)測(cè)的發(fā)展歷程與成就[J].中國(guó)水土保持科學(xué),2009,7(5):19. GUO Suoyan,LI Zhiguang.Development and achievements of soil and water conservation monitoring in China [J].Science of Soil and Water Conser1ation,2009,7 (5):19.

    [6] 土壤侵蝕分類分級(jí)標(biāo)準(zhǔn):SL190—2007[S].北京:中國(guó)水利水電出版社,1997:4. Standards for classification and gradation of soil erosion: SL 190-2007[S].Beijing:China Water Power Press, 1997:4.

    [7] KIRKBY M J,ABRAHART R,MCMAHON M D,et al. MEDALUS soil erosion models for global change[J].Geomorphology,1998,24(1):35.

    [8] KIRKBY M J,IMESON A C,BERGKAMP G J J,et al. Scaling up processes and models from the field plot to the watershed and regional areas[J].Journal of Soil and Water Conservation,1996,51(5):391.

    [9] ZHANG X Y,DRAKE N A,WAINWRIGHT J W.Scaling land surface parameters for global-scale soil erosion estimation[J].Water Resources Research,2002,38(9): 1180.

    [10]KING D,FOX D M,DAROUSSIN J,et al.Upscaling a simple erosion model from small areas to a large region [J].Nutrient Cycling in Agroecosystems,1998,50(1): 143.

    [11]NUSSER S M,GOEBEL J J.The national resources inventory:a long term resource monitoring programme[J].Environmental and Ecological Statistics,1997,4(3):182.

    [12]楊勤科,李銳,劉詠梅.區(qū)域土壤侵蝕普查方法的初步討論[J].中國(guó)水土保持科學(xué),2008,6(3):1. YANG Qinke,LI Rui,LIU Yongmei.Discussions on methods of national regional soil erosion survey[J].Science of Soil and Water Conservation,2008,6(3):1.

    [13]李智廣,符素華,劉寶元.我國(guó)水力侵蝕抽樣調(diào)查方法[J].中國(guó)水土保持科學(xué),2012,10(1):77. LI Zhiguang,FU Suhua,LIU Baoyuan.Sampling program of water erosion inventory in the first national water resource survey[J].Science of Soil and Water Conservation,2012,10(1):77.

    [14]張巖,劉憲春,李智廣,等.利用侵蝕模型普查黃土高原土壤侵蝕狀況[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(10): 165. ZHANG Yan,LIU Xianchun,LI Zhiguang,et al.Surveying soil erosion condition in the Loess Plateau using soil erosion model[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(10):165.

    [15]LIU Baoyuan,ZHANG Keli,XIE Yun.An empirical soil loss equation[C]∥Proceedings 12th International Soil Conservation Organization Conference.Beijing,China: Tsinghua University Press,2002,3:21.

    [16]章文波,付金生.不同類型雨量資料估算降雨侵蝕力[J].資源科學(xué),2003,25(1):36. ZHANG Wenbo,FU Jinsheng.Rainfall erosivity estimation under different rainfall amount[J].Resources Science, 2003,25(1):36.

    Estimation of regional soil erosion based on sampling survey

    WANG Zhenying1,HE Yuanmei2
    (1.Liaoning Academy of Water Resources and Hydropower,110003,Shenyang,China; 2.College of Soil and Water Conservation,Beijing Forestry University,100083,Beijing,China)

    [Background]The estimation of regional soil erosion is an important but difficult topic for soil erosion research.During the past years,artificial survey method had been used to make national soil erosion survey,classifying the soil erosion intensity in the topographic map according to vegetation coverage and slope gradient was time-consuming.The national soil erosion investigation based on remote sensing image and GIS technology in large scale had also be conducted,the soil erosion type and intensity were evaluated only according to vegetation coverage and slope,but without rainfall,soil erodibility and other important factors.Thus,how to present the regional spatial distribution of soil erosion quickly and accurately is an urgent issue.[Methods]In order to study method of estimating the regional soil erosion by survey data,Liaoning Province was taken as the study area,soil erosion factors(rainfall erosivity,soil erodibility,topography,soil and water conservation measures)and soil erosion modulus in 2011 were estimated using Chinese Soil Loss Equation(CSLE)based on GIS technology and sampling data from the first nationwide water resources survey in 2011.The soil erosion intensity of the whole Liaoning Province were graded,and compared with the statistical results in sampling plots.[Results]The results showed that the estimated soil erosion factors were relatively consistent with the sampling results apart from slope length and gradient.And the difference of scale accuracy of the topographic data mainly contributed to the difference of slope length and gradient between the whole Liaoning Province and sampling plots.Thedata from sampling plots and the whole province both showed the similar results that the intensity of soil erosion in Liaoning Province were mostly slight and light.Every year the intensity less than 1 000 t/km2accounted for about 85%of the area,while the one greater than 5 000 t/km2accounted for only 2%of the area.The soil erosion modulus of other land was the largest(1 777.65 t/km2),folllowed by the dryland(842.20 t/km2)and grassland(794.87 t/km2),the urban residence(11.69 t/km2),paddy field and water areas(17.14 t/km2)were the smallest.Compared to the estimated results of the whole province,the erosion modulus of sampling plots in the independent industrial and mining area,other garden and other land were larger with average values of 3 974.89 t/km2,1 656.95 t/km2and 807.74 t/ km2respectively,while the erosion modulus of sampling plots in grassland and other woodland were smaller with average values of 551.48 t/km2and 197.55 t/km2respectively,and these differences might resulted from the different resolutions of DEM data between the whole province and sampling plots[Conclusions]The research suggests that it is highly credible to estimate regional soil erosion modulus relying on sampling survey data and soil erosion model.And it can provide a reference for the estimation of regional soil erosion using sampling survey data.

    sampling data;CSLE;ArcGIS;regional soil erosion;Liaoning Province

    S157.1

    :A

    :2096-2673(2017)01-0015-07

    10.16843/j.sswc.2017.01.003

    2015- 10- 28

    2016- 11- 17

    王振穎(1966—),男,博士,教授級(jí)研究員。主要研究方向:3S技術(shù)應(yīng)用。E-mail:173667330@qq.com

    猜你喜歡
    坡長(zhǎng)模數(shù)土壤侵蝕
    基于單片機(jī)和模數(shù)化設(shè)計(jì)的低壓側(cè)電壓監(jiān)視與保護(hù)裝置
    能源工程(2021年2期)2021-07-21 08:40:02
    地邊截水地物對(duì)黑土區(qū)小流域坡長(zhǎng)因子計(jì)算的影響
    模數(shù)化設(shè)計(jì)方法在景觀鋪裝設(shè)計(jì)中的應(yīng)用
    綠色科技(2020年11期)2020-08-01 02:23:58
    鄉(xiāng)村聚落土壤侵蝕環(huán)境與水土流失研究綜述
    坡長(zhǎng)對(duì)貴州喀斯特區(qū)黃壤坡耕地土壤侵蝕的影響
    基于LID模式的城區(qū)排澇模數(shù)探析
    一種新型的RSA密碼體制模數(shù)分解算法
    海壇島土壤侵蝕問題研究
    大別山區(qū)土壤侵蝕動(dòng)態(tài)變化及趨勢(shì)預(yù)測(cè)
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟(jì)損失估算
    亚洲av日韩精品久久久久久密| 国产精品一区二区精品视频观看| 国产精品野战在线观看| 9色porny在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 非洲黑人性xxxx精品又粗又长| 在线观看免费午夜福利视频| 两性夫妻黄色片| 波多野结衣av一区二区av| 男女床上黄色一级片免费看| 亚洲一码二码三码区别大吗| 大码成人一级视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲专区中文字幕在线| 国产精品二区激情视频| 亚洲精品久久国产高清桃花| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品合色在线| 日韩 欧美 亚洲 中文字幕| 久久中文字幕一级| 欧美黄色淫秽网站| 一本大道久久a久久精品| 国产欧美日韩一区二区三区在线| 88av欧美| 制服诱惑二区| 中文字幕久久专区| 国产视频一区二区在线看| 国产xxxxx性猛交| 久久精品人人爽人人爽视色| 成年人黄色毛片网站| www.熟女人妻精品国产| 女同久久另类99精品国产91| 18禁国产床啪视频网站| 一本久久中文字幕| 中文字幕最新亚洲高清| 亚洲人成网站在线播放欧美日韩| avwww免费| 免费观看人在逋| 久99久视频精品免费| 午夜久久久久精精品| 色婷婷久久久亚洲欧美| 啦啦啦免费观看视频1| 精品久久久久久久毛片微露脸| 真人一进一出gif抽搐免费| 在线播放国产精品三级| 午夜激情av网站| 国产真人三级小视频在线观看| 久久久精品欧美日韩精品| 亚洲色图 男人天堂 中文字幕| 色尼玛亚洲综合影院| 无限看片的www在线观看| 午夜免费成人在线视频| 久久久国产精品麻豆| 在线观看免费午夜福利视频| 在线观看舔阴道视频| av电影中文网址| 悠悠久久av| 欧美黄色片欧美黄色片| 这个男人来自地球电影免费观看| 久久影院123| 九色亚洲精品在线播放| 亚洲第一av免费看| 久久亚洲真实| 久久午夜综合久久蜜桃| 亚洲电影在线观看av| 精品国产美女av久久久久小说| av网站免费在线观看视频| 精品欧美一区二区三区在线| 在线观看66精品国产| 国产免费av片在线观看野外av| 99国产综合亚洲精品| 老司机午夜十八禁免费视频| 国产免费男女视频| 一a级毛片在线观看| 免费看美女性在线毛片视频| 成人18禁在线播放| 老汉色∧v一级毛片| 日本精品一区二区三区蜜桃| 久久午夜亚洲精品久久| 91精品国产国语对白视频| 一区二区三区高清视频在线| 变态另类丝袜制服| 性少妇av在线| 亚洲国产精品合色在线| 久久亚洲精品不卡| 视频在线观看一区二区三区| xxx96com| 身体一侧抽搐| 亚洲熟妇熟女久久| 美女 人体艺术 gogo| 国产成人免费无遮挡视频| 女人被躁到高潮嗷嗷叫费观| 亚洲无线在线观看| 欧美中文日本在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 91成年电影在线观看| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 久久香蕉精品热| 1024香蕉在线观看| 别揉我奶头~嗯~啊~动态视频| 我的亚洲天堂| 久久婷婷成人综合色麻豆| 午夜福利,免费看| 国产色视频综合| 正在播放国产对白刺激| 国产精品综合久久久久久久免费 | 国产成人啪精品午夜网站| 一边摸一边抽搐一进一出视频| 国产成人系列免费观看| 日本 欧美在线| 欧美国产日韩亚洲一区| 国产午夜福利久久久久久| 欧美黄色片欧美黄色片| 18禁黄网站禁片午夜丰满| 精品国产乱子伦一区二区三区| 精品少妇一区二区三区视频日本电影| 国产精品久久久久久精品电影 | 亚洲第一欧美日韩一区二区三区| 一级a爱视频在线免费观看| 女人高潮潮喷娇喘18禁视频| 男女之事视频高清在线观看| 中文字幕色久视频| 久热爱精品视频在线9| 日韩精品免费视频一区二区三区| 母亲3免费完整高清在线观看| 亚洲av电影不卡..在线观看| 精品一区二区三区视频在线观看免费| 在线天堂中文资源库| 后天国语完整版免费观看| 少妇的丰满在线观看| 午夜免费鲁丝| 亚洲国产欧美日韩在线播放| 国产一区二区三区在线臀色熟女| av片东京热男人的天堂| 亚洲在线自拍视频| 女人被躁到高潮嗷嗷叫费观| 99精品欧美一区二区三区四区| 叶爱在线成人免费视频播放| 免费久久久久久久精品成人欧美视频| 亚洲三区欧美一区| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 在线国产一区二区在线| 99re在线观看精品视频| 视频在线观看一区二区三区| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 99在线人妻在线中文字幕| 久9热在线精品视频| 欧美最黄视频在线播放免费| 国产野战对白在线观看| 麻豆av在线久日| 亚洲国产精品久久男人天堂| 国产成人欧美| 国产精品精品国产色婷婷| 首页视频小说图片口味搜索| 变态另类成人亚洲欧美熟女 | 波多野结衣巨乳人妻| 亚洲全国av大片| 国产亚洲欧美在线一区二区| 久久久久国产一级毛片高清牌| 国产精品国产高清国产av| 一区二区日韩欧美中文字幕| 欧美成狂野欧美在线观看| 脱女人内裤的视频| 91成人精品电影| 亚洲欧美一区二区三区黑人| 中文字幕最新亚洲高清| 精品人妻1区二区| 亚洲九九香蕉| xxx96com| 免费看a级黄色片| 国产伦人伦偷精品视频| 久久久久精品国产欧美久久久| 999久久久精品免费观看国产| 啦啦啦 在线观看视频| 在线观看免费午夜福利视频| 美女国产高潮福利片在线看| 岛国视频午夜一区免费看| 女警被强在线播放| 色婷婷久久久亚洲欧美| 亚洲成av片中文字幕在线观看| 在线观看免费午夜福利视频| 欧美精品亚洲一区二区| 老汉色av国产亚洲站长工具| 色综合欧美亚洲国产小说| 亚洲 欧美一区二区三区| 嫩草影视91久久| 亚洲欧美一区二区三区黑人| 精品国内亚洲2022精品成人| 亚洲第一av免费看| 亚洲色图 男人天堂 中文字幕| 亚洲av成人一区二区三| 少妇被粗大的猛进出69影院| 日本一区二区免费在线视频| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 伦理电影免费视频| 久久久国产欧美日韩av| 男人操女人黄网站| 亚洲av五月六月丁香网| 女生性感内裤真人,穿戴方法视频| 操美女的视频在线观看| 一二三四在线观看免费中文在| 国产黄a三级三级三级人| 女人被狂操c到高潮| 国产乱人伦免费视频| 亚洲精品美女久久久久99蜜臀| 精品无人区乱码1区二区| 亚洲va日本ⅴa欧美va伊人久久| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 99国产精品一区二区三区| 久久亚洲精品不卡| e午夜精品久久久久久久| 少妇熟女aⅴ在线视频| 久久天堂一区二区三区四区| 97超级碰碰碰精品色视频在线观看| 欧美乱色亚洲激情| 亚洲男人天堂网一区| 国产麻豆69| 日韩精品中文字幕看吧| 国产三级黄色录像| 高清在线国产一区| 在线视频色国产色| 国产极品粉嫩免费观看在线| 亚洲人成网站在线播放欧美日韩| а√天堂www在线а√下载| 黄色a级毛片大全视频| 在线视频色国产色| 如日韩欧美国产精品一区二区三区| 91在线观看av| 国产精品久久久久久精品电影 | 国内精品久久久久精免费| 天天躁夜夜躁狠狠躁躁| www.www免费av| 中文字幕av电影在线播放| 国产成人系列免费观看| 在线观看免费午夜福利视频| 午夜两性在线视频| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| avwww免费| 久99久视频精品免费| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 丰满人妻熟妇乱又伦精品不卡| 午夜精品国产一区二区电影| 男人的好看免费观看在线视频 | 中文字幕av电影在线播放| 国产不卡一卡二| 欧美日韩福利视频一区二区| 少妇 在线观看| 国产三级黄色录像| 久久人人爽av亚洲精品天堂| 久久精品aⅴ一区二区三区四区| 久久香蕉激情| 亚洲中文字幕日韩| 中文字幕另类日韩欧美亚洲嫩草| 国产aⅴ精品一区二区三区波| 神马国产精品三级电影在线观看 | 久久香蕉国产精品| 亚洲欧洲精品一区二区精品久久久| 久久久久久人人人人人| 色播亚洲综合网| 免费观看精品视频网站| 亚洲片人在线观看| 国内精品久久久久久久电影| 亚洲三区欧美一区| 日本免费a在线| 久久影院123| 国产亚洲精品第一综合不卡| 制服诱惑二区| 手机成人av网站| 在线天堂中文资源库| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 在线观看午夜福利视频| 无限看片的www在线观看| 国产私拍福利视频在线观看| 国产成人欧美在线观看| 国产一级毛片七仙女欲春2 | 免费高清在线观看日韩| 欧美久久黑人一区二区| 日韩 欧美 亚洲 中文字幕| 99香蕉大伊视频| 亚洲av五月六月丁香网| 国产精品自产拍在线观看55亚洲| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区| www国产在线视频色| 久久人人97超碰香蕉20202| 日韩高清综合在线| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| 波多野结衣巨乳人妻| 中文字幕人妻丝袜一区二区| 级片在线观看| 在线观看一区二区三区| 99久久综合精品五月天人人| 午夜影院日韩av| 禁无遮挡网站| 国产成人av激情在线播放| 亚洲精品国产一区二区精华液| 男女之事视频高清在线观看| 国产亚洲av嫩草精品影院| 脱女人内裤的视频| 久久精品aⅴ一区二区三区四区| 亚洲人成伊人成综合网2020| 亚洲最大成人中文| 91精品三级在线观看| 亚洲免费av在线视频| 国产亚洲av高清不卡| 久久精品人人爽人人爽视色| 伦理电影免费视频| 久久久久久人人人人人| 国产伦人伦偷精品视频| 亚洲第一电影网av| 97碰自拍视频| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 欧美+亚洲+日韩+国产| 亚洲七黄色美女视频| 国产aⅴ精品一区二区三区波| 婷婷精品国产亚洲av在线| a级毛片在线看网站| 精品国产乱子伦一区二区三区| 亚洲专区中文字幕在线| 国产av精品麻豆| 亚洲成av人片免费观看| 亚洲av片天天在线观看| 亚洲国产精品999在线| 亚洲黑人精品在线| 最近最新中文字幕大全电影3 | 女人被狂操c到高潮| 精品国产一区二区久久| 50天的宝宝边吃奶边哭怎么回事| 久久精品亚洲熟妇少妇任你| 手机成人av网站| 久久精品国产亚洲av高清一级| 9色porny在线观看| 波多野结衣av一区二区av| 97超级碰碰碰精品色视频在线观看| 午夜两性在线视频| 男女床上黄色一级片免费看| 色尼玛亚洲综合影院| 一进一出抽搐gif免费好疼| 亚洲黑人精品在线| 国产精品久久电影中文字幕| 成在线人永久免费视频| 午夜福利在线观看吧| 久久久精品国产亚洲av高清涩受| 精品一区二区三区av网在线观看| 老司机在亚洲福利影院| 神马国产精品三级电影在线观看 | 91成人精品电影| 成人亚洲精品一区在线观看| 午夜精品在线福利| 精品久久久久久久人妻蜜臀av | 男人舔女人下体高潮全视频| 制服丝袜大香蕉在线| 天堂影院成人在线观看| 免费无遮挡裸体视频| 色在线成人网| 色综合站精品国产| 亚洲情色 制服丝袜| 97碰自拍视频| 好男人在线观看高清免费视频 | 亚洲一区二区三区色噜噜| 91av网站免费观看| 在线十欧美十亚洲十日本专区| 国产三级黄色录像| 成人精品一区二区免费| 午夜福利18| 国产私拍福利视频在线观看| 亚洲熟女毛片儿| 日本欧美视频一区| 日本免费a在线| 热99re8久久精品国产| 久久国产精品人妻蜜桃| 午夜激情av网站| 国产精品久久久av美女十八| 天天添夜夜摸| tocl精华| 亚洲一区高清亚洲精品| 99国产精品一区二区三区| 欧美av亚洲av综合av国产av| 国产一区二区三区在线臀色熟女| 一区在线观看完整版| 久久午夜综合久久蜜桃| 久久精品亚洲熟妇少妇任你| 午夜福利高清视频| 国产蜜桃级精品一区二区三区| 国产欧美日韩一区二区精品| 欧美国产精品va在线观看不卡| 18禁美女被吸乳视频| 亚洲情色 制服丝袜| 精品国产一区二区三区四区第35| 色av中文字幕| 日韩免费av在线播放| av网站免费在线观看视频| 亚洲av五月六月丁香网| 亚洲无线在线观看| 亚洲第一欧美日韩一区二区三区| 此物有八面人人有两片| 十八禁人妻一区二区| 精品高清国产在线一区| 制服丝袜大香蕉在线| 国产人伦9x9x在线观看| 久久精品91蜜桃| 一个人免费在线观看的高清视频| 亚洲国产欧美网| 91成人精品电影| 国产精品 国内视频| 男男h啪啪无遮挡| 国产精品精品国产色婷婷| 成人18禁高潮啪啪吃奶动态图| 亚洲久久久国产精品| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费 | 丁香欧美五月| 人妻久久中文字幕网| 男人舔女人的私密视频| 国产av又大| 91字幕亚洲| 亚洲人成电影观看| 看黄色毛片网站| 国产精品精品国产色婷婷| 久久国产精品人妻蜜桃| 久热这里只有精品99| 别揉我奶头~嗯~啊~动态视频| 韩国精品一区二区三区| 午夜免费激情av| 久久人人爽av亚洲精品天堂| 视频区欧美日本亚洲| 高清在线国产一区| 欧美日韩精品网址| a级毛片在线看网站| 色综合婷婷激情| 性色av乱码一区二区三区2| 人妻丰满熟妇av一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 成年女人毛片免费观看观看9| 成人av一区二区三区在线看| av超薄肉色丝袜交足视频| 一边摸一边做爽爽视频免费| 一夜夜www| 怎么达到女性高潮| 精品久久久久久成人av| 欧美一区二区精品小视频在线| 啪啪无遮挡十八禁网站| 99久久99久久久精品蜜桃| videosex国产| 亚洲欧美日韩高清在线视频| 97碰自拍视频| 国产成人影院久久av| 咕卡用的链子| 一级片免费观看大全| 亚洲国产日韩欧美精品在线观看 | 欧美日韩亚洲国产一区二区在线观看| 日韩免费av在线播放| 国产精品日韩av在线免费观看 | 在线观看免费视频网站a站| 午夜精品久久久久久毛片777| 亚洲欧美激情在线| 夜夜爽天天搞| 身体一侧抽搐| 中国美女看黄片| 人人澡人人妻人| 久久久久久大精品| 亚洲在线自拍视频| 欧美一区二区精品小视频在线| 日本撒尿小便嘘嘘汇集6| 久久婷婷成人综合色麻豆| 99久久精品国产亚洲精品| 悠悠久久av| 久久久久国内视频| 97碰自拍视频| 欧美日韩福利视频一区二区| 国产精品久久久av美女十八| 精品一区二区三区四区五区乱码| 久久婷婷成人综合色麻豆| 狠狠狠狠99中文字幕| 真人做人爱边吃奶动态| 亚洲av日韩精品久久久久久密| 少妇熟女aⅴ在线视频| 日韩欧美免费精品| 国产精品98久久久久久宅男小说| 搡老妇女老女人老熟妇| 动漫黄色视频在线观看| a级毛片在线看网站| 成人永久免费在线观看视频| 18禁裸乳无遮挡免费网站照片 | 日韩欧美在线二视频| av天堂在线播放| 欧美在线一区亚洲| 日本黄色视频三级网站网址| 国产一区二区在线av高清观看| 国产亚洲精品一区二区www| 午夜福利在线观看吧| 国语自产精品视频在线第100页| 满18在线观看网站| 九色国产91popny在线| ponron亚洲| 给我免费播放毛片高清在线观看| 神马国产精品三级电影在线观看 | 久久婷婷人人爽人人干人人爱 | 搡老岳熟女国产| 国产人伦9x9x在线观看| 欧美在线黄色| 欧美中文综合在线视频| 国产精品九九99| 国产成人免费无遮挡视频| 久久久久久久精品吃奶| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 9191精品国产免费久久| 黑丝袜美女国产一区| 午夜激情av网站| 亚洲av第一区精品v没综合| 香蕉久久夜色| 美女高潮喷水抽搐中文字幕| 免费在线观看亚洲国产| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频 | 欧美不卡视频在线免费观看 | 国产高清视频在线播放一区| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av在线| 国产精品秋霞免费鲁丝片| 久久国产精品影院| 亚洲午夜理论影院| 亚洲av五月六月丁香网| 俄罗斯特黄特色一大片| 国产片内射在线| 亚洲激情在线av| 97超级碰碰碰精品色视频在线观看| 九色亚洲精品在线播放| 18美女黄网站色大片免费观看| 天天一区二区日本电影三级 | x7x7x7水蜜桃| 在线观看免费视频日本深夜| 如日韩欧美国产精品一区二区三区| 欧美日本视频| av视频免费观看在线观看| 亚洲欧洲精品一区二区精品久久久| 国产三级在线视频| 久99久视频精品免费| 欧美精品啪啪一区二区三区| 免费高清视频大片| 国产精品久久电影中文字幕| 国产精品影院久久| 亚洲精华国产精华精| 中文字幕精品免费在线观看视频| 99热只有精品国产| 9色porny在线观看| av视频免费观看在线观看| 亚洲天堂国产精品一区在线| 韩国精品一区二区三区| 久久久久国内视频| 身体一侧抽搐| 淫秽高清视频在线观看| 日韩大尺度精品在线看网址 | 国产精品 国内视频| 男男h啪啪无遮挡| 国产一区二区激情短视频| 久久久久久久久免费视频了| 国产欧美日韩精品亚洲av| 久久精品国产99精品国产亚洲性色 | 午夜老司机福利片| 一区二区三区激情视频| 亚洲五月天丁香| 少妇 在线观看| 国产精品久久视频播放| 国产免费男女视频| 51午夜福利影视在线观看| 国产精品电影一区二区三区| or卡值多少钱| 国产精品一区二区在线不卡| 久久人人精品亚洲av| 久久这里只有精品19| 神马国产精品三级电影在线观看 | 成人三级做爰电影| 国产高清视频在线播放一区| 日本一区二区免费在线视频| 少妇 在线观看| 90打野战视频偷拍视频| 香蕉久久夜色| 国产成人精品在线电影| www日本在线高清视频| 欧美日韩亚洲综合一区二区三区_| 无限看片的www在线观看| 99国产精品免费福利视频| 日本一区二区免费在线视频| 免费看a级黄色片| 日韩有码中文字幕| 韩国精品一区二区三区| 一区二区三区精品91| 国产精品1区2区在线观看.| 一区二区三区精品91| 成人国产综合亚洲| 一级,二级,三级黄色视频| 亚洲专区字幕在线| 久久精品91无色码中文字幕| 色哟哟哟哟哟哟| 99riav亚洲国产免费| 久久香蕉精品热| 久99久视频精品免费| 精品国产亚洲在线| 国产欧美日韩一区二区精品| 精品一区二区三区视频在线观看免费| 青草久久国产| 国产一区二区三区视频了| 免费看美女性在线毛片视频|