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

    退耕還林還草工程實(shí)施對(duì)洛河流域土壤侵蝕的影響

    2016-10-26 03:24:36崔曉臨郝寧燕付立萌李貝貝
    水土保持研究 2016年5期
    關(guān)鍵詞:洛河土壤侵蝕降雨

    崔曉臨,雷 剛,王 濤,郝寧燕,付立萌,李貝貝

    (西安科技大學(xué) 測繪科學(xué)與技術(shù)學(xué)院,西安 710054)

    ?

    退耕還林還草工程實(shí)施對(duì)洛河流域土壤侵蝕的影響

    崔曉臨,雷 剛,王 濤,郝寧燕,付立萌,李貝貝

    (西安科技大學(xué) 測繪科學(xué)與技術(shù)學(xué)院,西安 710054)

    退耕還林還草工程是中國實(shí)施的重要生態(tài)環(huán)境建設(shè)與保護(hù)工程,對(duì)區(qū)域植被覆蓋及土壤侵蝕產(chǎn)生重要影響。以洛河流域(陜北黃土高原部分)為研究對(duì)象,利用流域通用土壤侵蝕方程(RUSLE),結(jié)合流域降雨、土壤類型、DEM、植被覆蓋等數(shù)據(jù),定量分析了2000—2010年退耕還林還草工程實(shí)施對(duì)流域土壤侵蝕的影響。結(jié)果表明:(1)洛河流域2000—2010年耕地面積減少,林地、草地面積增加,土地利用變化主要發(fā)生在2000—2005年;(2)洛河流域2000—2010年土地利用變化導(dǎo)致植被NDVI平均值增大,耕地變化區(qū)域植被NDVI值增加幅度高于耕地未變化區(qū)域,表明耕地變化區(qū)域植被NDVI增加對(duì)耕地區(qū)域總體植被NDVI值增加貢獻(xiàn)較大;(3)降雨侵蝕力和退耕還林還草工程實(shí)施對(duì)土壤侵蝕具有明顯的影響。受降雨侵蝕力增大影響,2000—2010年洛河流域土壤侵蝕呈增加趨勢;不考慮降雨侵蝕力變化情況下,洛河流域土壤侵蝕呈減少趨勢,反映出退耕還林還草工程實(shí)施對(duì)土壤侵蝕的減緩作用。

    土壤侵蝕; 植被覆蓋; RUSLE; 洛河流域

    退耕還林還草工程是中國規(guī)模較大的生態(tài)環(huán)境建設(shè)與保護(hù)工程,旨在改善區(qū)域生態(tài)環(huán)境,促進(jìn)農(nóng)村經(jīng)濟(jì)結(jié)構(gòu)調(diào)整及提高農(nóng)民生活水平[1],1999年開始在全國部分縣區(qū)開展退耕工作。陜北黃土高原是中國水土流失嚴(yán)重、生態(tài)環(huán)境脆弱區(qū)域,退耕還林還草工程實(shí)施,對(duì)區(qū)域植被覆蓋以及土壤侵蝕產(chǎn)生了重要影響。盡管退耕政策已實(shí)施多年,但土壤侵蝕狀況仍較嚴(yán)重,是黃河及其支流泥沙的主要來源地[2],因此陜北黃土高原土壤侵蝕研究具有重要理論與現(xiàn)實(shí)意義。

    土壤侵蝕是指土壤及其母質(zhì)在水力、風(fēng)力、凍融、重力等外力作用下,被破壞、剝蝕、搬運(yùn)和沉積的過程[3]。長期的土壤侵蝕會(huì)造成河流泥沙淤積、土地生產(chǎn)力下降,還會(huì)加劇洪澇和干旱等自然災(zāi)害[4],因此,土壤侵蝕已成為生態(tài)環(huán)境建設(shè)與保護(hù)不可忽視的問題。Wischmeier等[5]提出了通用土壤流失方程,在此基礎(chǔ)上一些學(xué)者提出了ANSWERS,CREAMS和AGNPS等土壤侵蝕模型[6-8]。20世紀(jì)80年代以后,國內(nèi)研究人員利用美國通用土壤流失方程(USLE)結(jié)合GIS(Geographical Information System)和RS(Remote Sensing)技術(shù),根據(jù)研究區(qū)的實(shí)際情況算出因子參數(shù),從而得到土壤侵蝕經(jīng)驗(yàn)?zāi)P?。目前來看,盡管已經(jīng)建立基于中國國情的CSLE(Chinese Soil Loss Equation)模型[9],但RUSLE模型仍是土壤侵蝕研究當(dāng)中應(yīng)用較為廣泛的模型[10]。

    洛河和延河均是陜北黃土高原地區(qū)重要的河流,但延河流域研究明顯多于洛河流域,如謝紅霞等[11]利用RUSLE模型分析了1997年和2000年延河流域土壤侵蝕狀況,認(rèn)為退耕還林還草工程實(shí)施引起的植被覆蓋改善以及降雨侵蝕力的增加對(duì)流域土壤侵蝕均產(chǎn)生了重要影響,但降雨侵蝕力增加導(dǎo)致的土壤侵蝕量明顯高于植被覆蓋改善減少的土壤侵蝕量。汪邦穩(wěn)等[12]利用1986年、1997年和2000年數(shù)據(jù)分析了延河流域退耕前后土壤侵蝕變化??梢?,退耕還林還草工程對(duì)延河流域土壤影響受到更多的關(guān)注,而洛河流域相對(duì)較少。

    鑒于謝紅霞[11]、汪邦穩(wěn)[12]等使用數(shù)據(jù)僅限于2000年前,而第一輪退耕還林還草工程在2008年已經(jīng)結(jié)束,且于2014年底又重啟了第二輪,有必要對(duì)退耕還林還草工程長期實(shí)施對(duì)土壤侵蝕影響進(jìn)行分析,可為新一輪工程實(shí)施提供理論支持與政策建議?;诖耍疚囊匝芯枯^少的洛河流域(陜北黃土高原部分)為對(duì)象,利用應(yīng)用較為成熟、廣泛的RUSLE模型,結(jié)合流域降雨數(shù)據(jù)、DEM、植被覆蓋、土壤類型等數(shù)據(jù),分析2000—2010年退耕還林還草工程實(shí)施對(duì)土壤侵蝕的影響。

    1 試驗(yàn)材料與方法

    1.1研究區(qū)概況

    研究對(duì)象主要為洛河流域位于陜北黃土高原的部分,地處中國西北部陜西省境內(nèi)(107°33′—109°51′E,35°19′—37°19′N),包括榆林地區(qū)的定邊縣和靖邊縣部分區(qū)域,延安地區(qū)的吳起縣、志丹縣、甘泉縣、富縣、黃陵縣和洛川縣,以及黃龍縣部分區(qū)域。流域內(nèi)海拔為688~1 900 m,西北部地勢較高,中部洛河河谷區(qū)域地勢較低。地貌屬黃土丘陵溝壑地貌,植被覆蓋南部區(qū)域高于北部,且南部區(qū)域以林地為主,而北部區(qū)域以草地為主。流域所在區(qū)域?qū)贉貛Т箨懶约撅L(fēng)氣候,年平均氣溫為14.1℃,年均降水量為533 mm,降水主要集中在夏秋季節(jié),同時(shí)也是土壤侵蝕的高發(fā)時(shí)期。夏秋季的集中降水、地表植被覆蓋狀況以及地勢、土壤類型等決定了洛河流域也是土壤侵蝕較為嚴(yán)重的區(qū)域。

    1.2研究方法

    1.2.1數(shù)據(jù)來源研究使用數(shù)據(jù)包括:(1)陜北黃土高原地區(qū)30 m分辨率DEM數(shù)據(jù),來源于中國科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國際科學(xué)數(shù)據(jù)鏡像網(wǎng)站(http:∥www.gscloud.cn)。DEM的用途有二:一是利用DEM提取洛河流域范圍;二是利用DEM生成坡度,用于計(jì)算土壤侵蝕因子。(2)1991—2010年7個(gè)站點(diǎn)(榆林、橫山、綏德、延安、定邊、吳起、洛川)逐日降水?dāng)?shù)據(jù),結(jié)合洛河流域范圍,文中選擇定邊、吳起、延安和洛川4個(gè)站點(diǎn),用于計(jì)算2000年、2005年和2010年流域降雨侵蝕力。(3)1∶50萬陜西省土壤圖,利用洛河流域邊界,提取流域土壤類型空間分布圖。(4)250 m分辨率MODIS NDVI 13Q1-Level 3數(shù)據(jù)(2000年、2005年和2010年),收集自http:∥ladsweb.nasa.gov。(5)陜西省1∶10萬土地利用類型矢量圖(2000年、2005年和2010年),來源于“全國生態(tài)環(huán)境十年變化(2000—2010年)遙感調(diào)查與評(píng)估”項(xiàng)目。

    1.2.2土壤侵蝕計(jì)算土壤侵蝕計(jì)算采用RUSLE模型,表達(dá)式如下[13]:

    A=R·K·L·S·C·P

    (1)

    式中:A為年均土壤侵蝕量;R為降雨侵蝕力因子;K為土壤可侵蝕性因子;L為坡長因子;S為坡度因子;C為地表覆蓋因子;P為水土保護(hù)措施因子。

    (1)降雨侵蝕力因子(R)。R因子反映的是區(qū)域內(nèi)降雨對(duì)地表土壤的侵蝕能力。采用章文波等[14-15]提出的計(jì)算方法:

    (2)

    α=21.586β-7.1891

    (3)

    (4)

    式中:Ri為第i年的降雨侵蝕力[(MJ·mm)/(hm2·h·a)];k為一年中形成侵蝕性降水的天數(shù)(d);Pi為第i天的侵蝕性降雨量(mm);降雨≥12 mm/d是形成侵蝕性降雨的條件,小于該值則記為0 mm[16];α和β為系數(shù),由Pd12和Py12計(jì)算得到,其中Pd12為降雨≥12 mm/d的多年平均日降雨量(mm),Py12為降雨≥12 mm/d的年平均降水量(mm)。

    (2)土壤侵蝕因子(K)。K因子反映土壤抵抗降雨侵蝕的能力,通過EPIC模型計(jì)算[17]:

    (5)

    式中:M,F(xiàn),N為土壤中砂粒、粉沙粒和黏粒質(zhì)量分?jǐn)?shù)(%);T為土壤有機(jī)碳質(zhì)量分?jǐn)?shù)(%);δ=1-M/100。計(jì)算出的結(jié)果乘以0.1317轉(zhuǎn)化為國際制單位。

    (3)地形因子(LS)。LS因子包括坡長(L)因子和坡度(S)因子,基于DEM并參考Renard[18]、McCool[19]和Liu[20]等的計(jì)算方法得到,計(jì)算方法如下:

    (6)

    (7)

    (8)

    (9)

    式中:λ為坡長;α為坡長指數(shù);θ為坡度(%)。

    (4)植被覆蓋因子(C)。C因子反映地表覆蓋狀況對(duì)土壤侵蝕的影響,主要基于MODIS NDVI數(shù)據(jù)計(jì)算得到[21],計(jì)算方法如下:

    AVHRR NDVI=0.18MODIS NDVI+0.131

    (10)

    (11)

    式中:α和β取值分別為2,1;C值通過AVHRR NDVI計(jì)算得到,需要先將MODIS NDVI轉(zhuǎn)化為AVHRR NDVI再計(jì)算。

    (5)水土保持措施因子(P)。P因子計(jì)算通過以下公式[22]:

    P=0.2+0.03θ

    (12)

    式中:θ為坡度(%)。

    2 結(jié)果與分析

    2.1土地利用變化

    洛河流域2000年、2005年和2010年土地利用結(jié)構(gòu)均以林地、草地和耕地為主,三者占流域比例在95%以上,其他用地(水域、建設(shè)用地、未利用地)面積比例低于5%。受退耕還林還草工程實(shí)施影響,洛河流域2000—2010年林地、草地呈增加趨勢,二者所占面積比例分別由2000年的44.48%和33.70%增加到2010年的45.15%和36.86%,分別增加0.66%和3.16%;耕地面積呈減少趨勢,面積比例由2000年20.97%減少到2010年的16.87%,減少4.09%(圖1A),可見,2000—2010年洛河流域耕地面積減少最大,而草地面積增加最多。

    圖1 洛河流域2000—2010年土地利用結(jié)構(gòu)及各時(shí)段變化量

    變化量上,草地、耕地是洛河流域面積變化量最大的土地利用類型,其次為林地。2000—2010年,洛河流域草地和林地分別由2000年的6 044.21 km2,7 977.26 km2增加到2010年的6 623.81 km2,8 113.07 km2,分別增加了579.60 km2,135.81 km2;耕地面積由2000年的3 759.67 km2減少到2010年的3 032.19 km2,減少727.48 km2(圖1B)。變化階段上,2000—2005年是洛河流域土地利用劇烈變化時(shí)期,2005—2010年土地利用變化較緩。具體表現(xiàn)為,2000—2005年林地、草地增加面積分別占2000—2010年總增加面積的92.72%和84.29%,耕地減少面積占93.11%。以上分析表明,2000—2010年洛河流域土地利用以林地、草地和耕地為主,受退耕還林還草工程實(shí)施影響,流域內(nèi)草地面積增加最大,而耕地面積減少最多,并且土地利用劇烈變化主要發(fā)生在2000—2005年。

    2.2植被NDVI變化

    退耕還林還草工程實(shí)施引起的洛河流域土地利用變化,必然影響流域內(nèi)植被覆蓋變化,植被NDVI是植被覆蓋的重要反映。基于2000—2010年洛河流域耕地?cái)?shù)據(jù),獲取2000年耕地區(qū)域、2000—2010年耕地未發(fā)生變化區(qū)域和2000—2010年耕地變化區(qū)域,結(jié)合2000年、2005年和2010年流域MODIS NDVI數(shù)據(jù),得到不同耕地區(qū)域植被NDVI變化情況,見表1。

    由表1可知,洛河流域2000—2010年植被NDVI均值較高,反映出流域內(nèi)植被覆蓋水平較高。對(duì)于2000年為耕地的區(qū)域,總體上植被NDVI值由2000年的0.501 7增加至2010年的0.625 9,增加24.76%。植被NDVI值的增加得益于退耕還林還草工程實(shí)施引起的耕地向林地、草地的轉(zhuǎn)變。結(jié)合2000—2010年耕地區(qū)域是否發(fā)生變化,即2000年是耕地,2010年是否轉(zhuǎn)變?yōu)槠渌恋乩妙愋?,劃分為耕地未變化區(qū)域和耕地變化區(qū)域??芍?,2000年耕地區(qū)域的植被NDVI變化,受到該兩部分區(qū)域的影響。對(duì)于2000—2010年耕地未變化區(qū)域,植被NDVI值由2000年的0.511 6增加至2010年的0.627 5,增加22.65%;而2000—2010年耕地變化區(qū)域,植被NDVI值由2000年的0.461 3增加至2010年的0.619 5,增加34.39%。對(duì)比可知,2000—2010年耕地未變化區(qū)域植被NDVI值明顯高于同期耕地變化區(qū)域,但耕地變化區(qū)域植被NDVI值增加對(duì)耕地區(qū)域植被NDVI值增加貢獻(xiàn)較大。其原因可能為,退耕還林還草工程實(shí)施,主要針對(duì)坡耕地及生產(chǎn)力較低的耕地,該部分耕地植被NDVI較??;同時(shí),由于耕地轉(zhuǎn)變?yōu)榱值亍⒉莸爻跗?,植被覆蓋度較小,植被NDVI值較低,通過較長時(shí)間的生長發(fā)育,植被覆蓋率升高,植被NDVI值逐漸增加。

    表1 不同耕地區(qū)域植被NDVI均值變化

    2.3土壤侵蝕變化

    退耕還林還草工程實(shí)施,促使洛河流域耕地向林地、草地的轉(zhuǎn)變,從而引起植被NDVI的變化,植被覆蓋得到一定程度的改善。對(duì)于RUSLE模型估算流域土壤侵蝕,模型中除降雨侵蝕力(R)、植被覆蓋和管理因子(C)每年均有變化外,其他因子各年份的變動(dòng)很小,故不同年份土壤侵蝕的估算值,主要受到R因子和C因子的影響。結(jié)合前文可知,退耕還林還草工程實(shí)施影響到C因子的變化,而R因子作為自然因素,不受控于人為作用。通過計(jì)算R因子可知,2000年、2005年和2010年洛河流域降雨侵蝕力平均值分別為1 236,1 326,1 550 (MJ·mm)/(hm2·h·a),呈增加趨勢。

    利用RUSLE模型計(jì)算得到洛河流域2000年、2005年和2010年年土壤侵蝕量,參照國家標(biāo)準(zhǔn)[23],將土壤侵蝕程度劃分為6個(gè)等級(jí),分別是1級(jí)[微度侵蝕,<1 000 t/(km2·a)]、2級(jí)[輕度侵蝕,1 000~2 500 t/(km2·a)]、3級(jí)[中度侵蝕,2 500~5 000 t/(km2·a)],4級(jí)[強(qiáng)烈侵蝕,5 000~8 000 t/(km2·a)]、5級(jí)[極強(qiáng)烈侵蝕,8 000~15 000 t/(km2·a)]、6級(jí)[劇烈侵蝕,≥15 000 t/(km2·a)]。統(tǒng)計(jì)各等級(jí)所占面積比例及面積比例在研究時(shí)段內(nèi)的變化情況,得到圖2A和圖2B。鑒于降雨侵蝕力(R)具有增加趨勢,為了反映退耕還林還草工程實(shí)施的影響,在計(jì)算土壤侵蝕時(shí),利用2000年降雨侵蝕力重新計(jì)算2005年和2010年土壤侵蝕量,統(tǒng)計(jì)各等級(jí)面積比例及面積比例在研究時(shí)段的變化情況,得到圖2C和圖2D。

    注:A和B利用2000年、2005年和2010年降雨侵蝕力計(jì)算得到;C和D利用2000年降雨侵蝕力重新計(jì)算得到。

    圖2洛河流域2000-2010年不同土壤侵蝕等級(jí)面積比例及變化

    圖2A和圖2C表明,洛河流域2000—2010年土壤侵蝕以微度、輕度和中度侵蝕為主,其中微度侵蝕面積占80%以上,強(qiáng)烈、極強(qiáng)烈和劇烈侵蝕面積比例均較小。2000年、2005年和2010年降雨侵蝕力計(jì)算結(jié)果表明,洛河流域2000—2010年除微度侵蝕面積呈減少趨勢外,其他各土壤侵蝕等級(jí)均呈增加趨勢,同時(shí)可以看出,2005—2010年是土壤侵蝕減少和增加的主要時(shí)期(圖2A和圖2B)。利用2000年降雨侵蝕力重新計(jì)算2005年和2010年土壤侵蝕,結(jié)果表明,與圖2A和圖2B結(jié)果相反,除微度侵蝕呈增加趨勢外,其他各土壤侵蝕等級(jí)均呈減少趨勢,且2000—2005年是土壤侵蝕增加和減少的主要時(shí)期。

    圖2B反映出受到降雨侵蝕力增加及植被覆蓋改善影響,洛河流域2000—2010年土壤侵蝕呈增加狀況;圖2D反映出,當(dāng)降雨侵蝕力維持在2000年水平時(shí),洛河流域2000—2010年土壤侵蝕受植被覆蓋改善的影響,表現(xiàn)為土壤侵蝕的減少。綜合二者可以認(rèn)為,降雨侵蝕力和植被覆蓋對(duì)土壤侵蝕具有顯著的影響,洛河流域2000—2010年降雨侵蝕力對(duì)土壤侵蝕的增加作用明顯高于植被覆蓋的減少作用,總體呈現(xiàn)為土壤侵蝕增加趨勢。

    2.4結(jié)果驗(yàn)證

    查閱黃河水資源和泥沙公報(bào)(2000—2010年)中洛河流域狀頭水文站2000—2010年年徑流量和年輸沙量數(shù)據(jù)(http:∥www.yellowriver.gov.cn),繪制得到圖3。洛河流域2000—2010年年徑流量呈波動(dòng)變化,但總體變化維持在6億m3左右,但年輸沙量呈明顯的減少趨勢(圖3)。

    圖3 洛河流域狀頭水文站2000-2010年年徑流量和年輸沙量變化

    在年徑流量變動(dòng)不大情況下,年輸沙量呈顯著的減少趨勢,反映出洛河流域2000—2010年土壤侵蝕呈下降趨勢,顯然文中RUSLE模型計(jì)算的土壤侵蝕結(jié)果與水文站實(shí)測結(jié)果存在差異。原因可能為RUSLE模型計(jì)算得到土壤侵蝕結(jié)果是空間靜態(tài)狀況,與實(shí)際河流的輸沙量之間存在不一致性,二者之間需要通過坡面徑流相聯(lián)系,計(jì)算結(jié)果表明土壤侵蝕呈增加趨勢,而實(shí)測結(jié)果為減少趨勢,可能與退耕還林還草工程實(shí)施引起的地表植被覆蓋改善,減緩坡面徑流、涵蓄水源能力增加有關(guān)。

    3 結(jié) 論

    (1)洛河流域2000—2010年耕地面積減少,共減少727.48 km2,減少4.09%;林地、草地面積增加,分別增加579.60,135.81 km2,增加0.66%和3.16%;洛河流域土地利用變化主要發(fā)生在2000—2005年期間。

    (2)洛河流域2000—2010年土地利用變化提高了植被平均NDVI值,流域內(nèi)耕地區(qū)域植被NDVI值增加了24.76%,其中,耕地變化區(qū)域植被NDVI值增加34.29%,高于耕地未變化區(qū)域的22.65%。表明耕地變化區(qū)域植被NDVI增加對(duì)耕地區(qū)域總體植被NDVI值增加貢獻(xiàn)較大。

    (3)降雨侵蝕力和退耕還林還草工程實(shí)施對(duì)土壤侵蝕具有明顯的影響。受降雨侵蝕力增大影響,2000—2010年洛河流域土壤侵蝕呈增加趨勢;不考慮降雨侵蝕力變化情況下,洛河流域土壤侵蝕呈減少趨勢,反映出退耕還林還草工程實(shí)施對(duì)土壤侵蝕的減緩作用。

    致謝:非常感謝西北大學(xué)城市與環(huán)境學(xué)院楊勤科教授、評(píng)審專家及編輯部老師提出的寶貴建議,感謝澳大利亞聯(lián)邦科學(xué)與工業(yè)研究組織(CSIRO Land and Water)傅國斌研究員對(duì)文章英文摘要的修改。特此致謝!

    [1]劉會(huì)靜,姜志德,王繼軍.黃土高原退耕區(qū)農(nóng)業(yè)后續(xù)產(chǎn)業(yè)發(fā)展影響因素的多層線性分析[J].經(jīng)濟(jì)地理,2014,34(2):125-130.

    [2]甘枝茂,桑廣書.關(guān)于黃土高原退耕還林(草)問題[J].干旱區(qū)資源與環(huán)境,2002,16(1):62-65.

    [3]韓富偉,張柏,宋開山,等.長春市土壤侵蝕潛在危險(xiǎn)度分級(jí)及侵蝕背景的空間分析[J].水土保持學(xué)報(bào),2007,21(1):39-43.

    [4]何耀幫,趙永蘭,田國行,等.基于RS和GIS的伊河流域土壤侵蝕研究[J].中國農(nóng)學(xué)通報(bào),2012,28(20):237-242.

    [5]Wischmeier W H,Smith D D.Predicting rainfall erosion losses from cropland east of the Rocky Mountains : A guide to conservation planning[M]∥Agriculture Handbook No.537.Washington D C: USDA,1978.

    [6]Beasley D B,Huggins L F,Monke E J.Answers: A model for watershed planning[J].Transactions of the ASAE,1980,23(4):938-944.

    [7]Young R A,Onstad C A,Bosch D D,et al.AGNPS: A nonpoint-source pollution model for evaluating agricultural watersheds[J].Journal of Soil and Water Conservation,1989,44(2):168-173.

    [8]Knisel W G.CREAMS: A field-scale model for chemicals,runoff and erosion from agricultural management systems[R].Washington D C: USDA,2010.

    [9]Liu Baoyuan,Zhang Keli,Xie Yun.An empirical soil loss equation[C]∥Proceedings 12th International Soil Conservation Organization Conference.Beijing: Tsinghua University Press,2002.

    [10]Fu G,Chen S,McCool D K.Modeling the impacts of no-till practice on soil erosion and sediment yield with RUSLE,SEDD,and ArcView GIS[J].Soil and Tillage Research,2006,85(1):38-49.

    [11]謝紅霞,李銳,楊勤科,等.退耕還林(草)和降雨變化對(duì)延河流域土壤侵蝕的影響[J].中國農(nóng)業(yè)科學(xué),2009,42(2):569-576.

    [12]汪邦穩(wěn),楊勤科,劉志紅,等.延河流域退耕前后土壤侵蝕強(qiáng)度的變化[J].中國水土保持科學(xué),2007,5(4):27-33.

    [13]Yoder D,Lown J.The future of RUSLE: inside the new revised universal soil loss equation[J].Journal of Soil and Water Conservation,1995,50(5):484-489.

    [14]章文波,謝云,劉寶元.利用日雨量計(jì)算降雨侵蝕力的方法研究[J].地理科學(xué),2002,22(6):705-711.

    [15]章文波,付金生.不同類型雨量資料估算降雨侵蝕力[J].資源科學(xué),2003,25(1):35-41.

    [16]謝云,劉寶元.侵蝕性降雨標(biāo)準(zhǔn)研究[J].水土保持學(xué)報(bào),2000,14(4):6-11.

    [17]Sharpley A N,Williams J R.EPIC-erosion/productivity impact calculator: 1.Model documentation[J].Technical Bulletin-United States Department of Agriculture,1990,4(4):206-207.

    [18]Renard K G,Foster G R,Weesies G A,et al.Predicting rainfall erosion by water: a guide to conservation planning with the Revised Universal Soil Loess Equation (RUSLE)[M]∥Agricultural Handbook N0.730.Washington D C: USDA,1997.

    [19]McCool D K,Brown L C,Foster G R,et al.Revised slope steepness factor for the Universal Soil Loss Equation[J].Transactions of the ASAE,1987,30(5):1387-1396.

    [20]Liu B Y,Nearing M A,Risse L M.Slope gradient effects on soil loss for steep slopes[J].Transactions of the ASAE,1994,37(6):1835-1840.

    [21]李天宏,鄭麗娜.基于RUSLE模型的延河流域2001—2010年土壤侵蝕動(dòng)態(tài)變化[J].自然資源學(xué)報(bào),2012,27(7):1164-1175.

    [22]Lufafa A,Tenywa M M,Isabirye M,et al.Prediction of soil erosion in a Lake Victoria basin catchment using a GIS-based Universal Soil Loss model[J].Agricultural Systems,2003,76(3):883-894.

    [23]中華人民共和國水利部.土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)SL190—2007[S].北京:中國水利水電出版社,2007.

    Impacts of Grain for Green Project on Soil Erosion in Luohe River Basin of Northern Shaanxi Province,China

    CUI Xiaolin,LEI Gang,WANG Tao,HAO Ningyan,FU Limeng,LI Beibei

    (College of Geomatics,Xi′an University of Science and Technology,Xi′an 710054,China)

    Grain for Green Project (GGP)is an important ecological and environmental construction and protection engineering project in China,and it has great influence on vegetation cover and soil erosion.The Revised Universal Soil Loss Equation (RUSLE),together with precipitation,soil type,DEM and vegetation cover,was used to investigate the impacts of GGP on soil erosion from 2000 to 2010 in Luohe River Basin,northern Shaanxi Province,China.Results showed that: (1)farmland decreased while forestland and grassland increased due to GGP,and land use change mainly occurred from 2000 to 2005; (2)the increase of vegetation NDVI value induced because of land use change from 2000 to 2010 in Luohe River Basin,and the increase rate of vegetation NDVI in farmland changing area was higher than that in the unchanged area; (3)both rainfall erosivity and GGP implementation have the significant impact on soil erosion.The increase of soil erosion in Luohe River Basin during 2000—2010 mainly resulted from the increase of rainfall erosivity.If rainfall erosivity keeps constant,then soil erosion in the basin would decrease because of the effects of GGP on soil erosion.

    soil erosion; vegetation cover; RUSLE; Luohe River Basin

    2015-09-15

    2015-10-05

    國家自然科學(xué)基金資助項(xiàng)目(41271103);黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(A314021402-1616);西安科技大學(xué)培育基金資助項(xiàng)目(2014007);西安科技大學(xué)博士啟動(dòng)金項(xiàng)目(2015QDJ061)

    崔曉臨(1965—),女,陜西寶雞人,博士,副教授,主要從事土地利用變化與環(huán)境效應(yīng)研究。E-mail:xiaolin_cui@xust.edu.cn

    王濤(1984—),男,河南湯陰人,博士,講師,主要從事區(qū)域環(huán)境變化研究。E-mail:wht432@163.com

    S157.1

    A

    1005-3409(2016)05-0068-06

    猜你喜歡
    洛河土壤侵蝕降雨
    鄉(xiāng)村聚落土壤侵蝕環(huán)境與水土流失研究綜述
    洛 河
    洛河
    滄州市2016年“7.19~7.22”與“8.24~8.25”降雨對(duì)比研究
    海壇島土壤侵蝕問題研究
    紅黏土降雨入滲的定量分析
    彬長礦區(qū)白堊系洛河組不同層段水力聯(lián)系的綜合分析
    地球(2016年7期)2016-04-14 22:00:20
    大別山區(qū)土壤侵蝕動(dòng)態(tài)變化及趨勢預(yù)測
    科洛河的傳說
    南方降雨不斷主因厄爾尼諾
    叶爱在线成人免费视频播放| 亚洲av日韩精品久久久久久密| 国产高清视频在线播放一区 | 新久久久久国产一级毛片| 性色av一级| xxxhd国产人妻xxx| 日日夜夜操网爽| 免费在线观看日本一区| 亚洲专区字幕在线| 无遮挡黄片免费观看| 久久毛片免费看一区二区三区| 精品免费久久久久久久清纯 | 极品少妇高潮喷水抽搐| av在线老鸭窝| 国产日韩欧美视频二区| 在线观看www视频免费| 777久久人妻少妇嫩草av网站| 老熟妇乱子伦视频在线观看 | 国产日韩欧美亚洲二区| 亚洲专区中文字幕在线| 国产成人欧美| 国产亚洲精品第一综合不卡| 久久人人97超碰香蕉20202| 久久久久久久大尺度免费视频| 精品国产超薄肉色丝袜足j| 99热全是精品| 欧美97在线视频| 黑丝袜美女国产一区| 欧美成狂野欧美在线观看| 18禁黄网站禁片午夜丰满| 国产成人欧美| 国产精品 国内视频| 国产亚洲精品久久久久5区| 天天添夜夜摸| 91av网站免费观看| 丝袜美腿诱惑在线| 亚洲精品中文字幕一二三四区 | 最近最新免费中文字幕在线| 欧美日本中文国产一区发布| 亚洲国产毛片av蜜桃av| 99国产极品粉嫩在线观看| 亚洲精品中文字幕在线视频| 欧美 日韩 精品 国产| 涩涩av久久男人的天堂| 男女午夜视频在线观看| 久久久国产欧美日韩av| 在线观看免费视频网站a站| h视频一区二区三区| 另类亚洲欧美激情| 国产欧美日韩一区二区三 | 精品福利永久在线观看| 亚洲国产精品一区二区三区在线| 久久久久国内视频| 叶爱在线成人免费视频播放| 国产精品一区二区精品视频观看| 性色av一级| 精品亚洲成国产av| 男人添女人高潮全过程视频| 美国免费a级毛片| 久久午夜综合久久蜜桃| 窝窝影院91人妻| 少妇猛男粗大的猛烈进出视频| 久久精品亚洲熟妇少妇任你| 亚洲三区欧美一区| 久久精品国产a三级三级三级| 国产91精品成人一区二区三区 | 美女中出高潮动态图| 在线精品无人区一区二区三| 日本黄色日本黄色录像| 性高湖久久久久久久久免费观看| 国产淫语在线视频| 久热这里只有精品99| 久久精品aⅴ一区二区三区四区| 宅男免费午夜| 亚洲av电影在线进入| 久久中文字幕一级| 伊人亚洲综合成人网| 香蕉国产在线看| 一本色道久久久久久精品综合| 欧美一级毛片孕妇| 久久国产精品人妻蜜桃| 日韩免费高清中文字幕av| 久久久久国产精品人妻一区二区| 亚洲精品成人av观看孕妇| 久久毛片免费看一区二区三区| 性色av乱码一区二区三区2| 精品国产一区二区久久| 国产精品熟女久久久久浪| videos熟女内射| 性色av乱码一区二区三区2| 99久久精品国产亚洲精品| 一级毛片精品| 亚洲第一av免费看| 免费女性裸体啪啪无遮挡网站| 国产97色在线日韩免费| 亚洲国产精品成人久久小说| 无遮挡黄片免费观看| 久久99一区二区三区| 亚洲精品乱久久久久久| 一区二区三区乱码不卡18| 欧美av亚洲av综合av国产av| 夫妻午夜视频| 亚洲欧美日韩高清在线视频 | 国产av又大| 成人三级做爰电影| 真人做人爱边吃奶动态| 美国免费a级毛片| 国产精品一区二区在线观看99| 飞空精品影院首页| 黄色视频在线播放观看不卡| 午夜日韩欧美国产| 一级毛片女人18水好多| 91精品国产国语对白视频| 9热在线视频观看99| 免费在线观看日本一区| 国产亚洲欧美精品永久| 国产日韩一区二区三区精品不卡| 在线观看免费高清a一片| 下体分泌物呈黄色| av天堂在线播放| 国产精品国产av在线观看| 久久久久久久久久久久大奶| 国产日韩一区二区三区精品不卡| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 午夜久久久在线观看| 97人妻天天添夜夜摸| 一区二区日韩欧美中文字幕| 免费高清在线观看视频在线观看| 国产精品一区二区在线不卡| 每晚都被弄得嗷嗷叫到高潮| 日韩三级视频一区二区三区| 国产成人欧美在线观看 | av福利片在线| 色播在线永久视频| 亚洲精品第二区| av视频免费观看在线观看| 亚洲自偷自拍图片 自拍| av天堂在线播放| 人成视频在线观看免费观看| 狂野欧美激情性xxxx| 国产亚洲精品一区二区www | 777久久人妻少妇嫩草av网站| 午夜影院在线不卡| 免费看十八禁软件| 免费av中文字幕在线| 999精品在线视频| 国产一区二区三区综合在线观看| 视频区图区小说| 免费日韩欧美在线观看| 18禁国产床啪视频网站| 欧美日韩视频精品一区| 日日摸夜夜添夜夜添小说| 午夜福利视频在线观看免费| 亚洲黑人精品在线| 黄色怎么调成土黄色| 蜜桃国产av成人99| 老汉色av国产亚洲站长工具| av网站在线播放免费| 无限看片的www在线观看| 午夜福利乱码中文字幕| 久久性视频一级片| 一本大道久久a久久精品| 久久久国产欧美日韩av| 欧美成狂野欧美在线观看| 老熟妇仑乱视频hdxx| 91精品三级在线观看| 免费不卡黄色视频| 在线天堂中文资源库| 最近最新中文字幕大全免费视频| 正在播放国产对白刺激| 国产亚洲精品第一综合不卡| 亚洲激情五月婷婷啪啪| 91精品国产国语对白视频| 黑人巨大精品欧美一区二区蜜桃| 久久久国产精品麻豆| av又黄又爽大尺度在线免费看| 免费日韩欧美在线观看| 国产欧美亚洲国产| av片东京热男人的天堂| 久久久久国产精品人妻一区二区| 国产在线视频一区二区| 成年女人毛片免费观看观看9 | 亚洲,欧美精品.| 一区二区三区精品91| 亚洲中文日韩欧美视频| 美女午夜性视频免费| 免费人妻精品一区二区三区视频| 午夜福利视频精品| 男人舔女人的私密视频| 欧美xxⅹ黑人| 亚洲欧洲精品一区二区精品久久久| 嫩草影视91久久| 日韩欧美免费精品| 亚洲 欧美一区二区三区| 纯流量卡能插随身wifi吗| 狠狠婷婷综合久久久久久88av| 亚洲国产欧美在线一区| 国产成人欧美| 啪啪无遮挡十八禁网站| 国产在视频线精品| 久久中文看片网| 亚洲伊人久久精品综合| 欧美日韩精品网址| 国产精品av久久久久免费| 日本猛色少妇xxxxx猛交久久| 不卡一级毛片| 肉色欧美久久久久久久蜜桃| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 亚洲欧美一区二区三区黑人| 国产成人a∨麻豆精品| av在线老鸭窝| 久久精品久久久久久噜噜老黄| 中文字幕精品免费在线观看视频| 国产在线视频一区二区| 狂野欧美激情性xxxx| 亚洲专区国产一区二区| av线在线观看网站| 欧美老熟妇乱子伦牲交| 成年美女黄网站色视频大全免费| 热99re8久久精品国产| 国产成人av教育| 午夜久久久在线观看| 免费日韩欧美在线观看| 极品少妇高潮喷水抽搐| 亚洲精品国产精品久久久不卡| 欧美日韩视频精品一区| 国产高清国产精品国产三级| 一本大道久久a久久精品| 国产一级毛片在线| 宅男免费午夜| 美女大奶头黄色视频| 亚洲国产av新网站| 免费一级毛片在线播放高清视频 | 在线观看一区二区三区激情| 色精品久久人妻99蜜桃| 国产成+人综合+亚洲专区| 久久精品国产亚洲av香蕉五月 | 国产在线免费精品| 一本色道久久久久久精品综合| 亚洲精品中文字幕在线视频| 人成视频在线观看免费观看| 国产精品久久久久成人av| 国产av一区二区精品久久| 国产精品秋霞免费鲁丝片| 久久久久网色| 国产免费现黄频在线看| 欧美少妇被猛烈插入视频| 国产亚洲av片在线观看秒播厂| 麻豆av在线久日| av电影中文网址| 久热这里只有精品99| 99热网站在线观看| 99热全是精品| 性色av乱码一区二区三区2| 久久热在线av| 男人舔女人的私密视频| 国产一区二区三区av在线| 精品少妇黑人巨大在线播放| 色综合欧美亚洲国产小说| 久久久久久久大尺度免费视频| 在线看a的网站| 亚洲精品av麻豆狂野| 欧美在线黄色| 国产一区二区三区av在线| 一本大道久久a久久精品| 国产一卡二卡三卡精品| 大片电影免费在线观看免费| av在线app专区| 1024视频免费在线观看| 亚洲va日本ⅴa欧美va伊人久久 | 免费在线观看黄色视频的| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 久久精品久久久久久噜噜老黄| 亚洲专区字幕在线| 老汉色av国产亚洲站长工具| 久久久精品94久久精品| 免费黄频网站在线观看国产| 一区二区av电影网| 亚洲欧美激情在线| 一级片免费观看大全| 国产成人精品久久二区二区免费| 久久久水蜜桃国产精品网| 国产男人的电影天堂91| 亚洲少妇的诱惑av| 美女主播在线视频| 国产成人a∨麻豆精品| 狠狠精品人妻久久久久久综合| 日本精品一区二区三区蜜桃| 69精品国产乱码久久久| 国产精品一区二区免费欧美 | 男女免费视频国产| 亚洲伊人久久精品综合| 亚洲专区字幕在线| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| 欧美另类亚洲清纯唯美| 欧美成人午夜精品| 午夜福利,免费看| 国产麻豆69| 亚洲国产欧美网| 波多野结衣av一区二区av| 国产黄色免费在线视频| 91老司机精品| 在线亚洲精品国产二区图片欧美| 欧美日韩视频精品一区| 侵犯人妻中文字幕一二三四区| 在线看a的网站| 日日摸夜夜添夜夜添小说| 高清视频免费观看一区二区| 午夜成年电影在线免费观看| 波多野结衣av一区二区av| 久9热在线精品视频| 黑人巨大精品欧美一区二区mp4| 色94色欧美一区二区| 50天的宝宝边吃奶边哭怎么回事| 一区二区日韩欧美中文字幕| 国产成人啪精品午夜网站| www.自偷自拍.com| 久久久国产成人免费| 女性被躁到高潮视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文日韩欧美视频| 亚洲精品粉嫩美女一区| 久久久久国产一级毛片高清牌| 午夜激情av网站| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区久久| 国产精品九九99| 我的亚洲天堂| 久久久欧美国产精品| 免费久久久久久久精品成人欧美视频| 国产亚洲av片在线观看秒播厂| 亚洲人成电影观看| 777久久人妻少妇嫩草av网站| 国产伦人伦偷精品视频| 在线观看舔阴道视频| 亚洲中文av在线| 国产人伦9x9x在线观看| 蜜桃在线观看..| 狂野欧美激情性xxxx| 人人妻,人人澡人人爽秒播| 国产av又大| www日本在线高清视频| 大陆偷拍与自拍| 9色porny在线观看| 宅男免费午夜| svipshipincom国产片| 久久人人97超碰香蕉20202| 亚洲一卡2卡3卡4卡5卡精品中文| av免费在线观看网站| 国产精品av久久久久免费| 久久青草综合色| 亚洲精品美女久久久久99蜜臀| 熟女少妇亚洲综合色aaa.| 国产日韩欧美在线精品| 精品一区二区三区四区五区乱码| 日韩欧美国产一区二区入口| 后天国语完整版免费观看| 婷婷色av中文字幕| 十八禁人妻一区二区| 一级黄色大片毛片| 一级,二级,三级黄色视频| 亚洲精品一卡2卡三卡4卡5卡 | 波多野结衣一区麻豆| 秋霞在线观看毛片| 女人久久www免费人成看片| 亚洲国产av影院在线观看| 亚洲自偷自拍图片 自拍| 80岁老熟妇乱子伦牲交| 国产精品亚洲av一区麻豆| 亚洲伊人色综图| 久9热在线精品视频| 国产精品免费视频内射| 国产免费福利视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 免费高清在线观看日韩| 99国产精品99久久久久| 人人妻人人澡人人爽人人夜夜| 我的亚洲天堂| 亚洲国产日韩一区二区| 啦啦啦中文免费视频观看日本| 欧美激情极品国产一区二区三区| 人人妻人人添人人爽欧美一区卜| 欧美日本中文国产一区发布| 看免费av毛片| 亚洲av日韩精品久久久久久密| 午夜福利一区二区在线看| 久久久久国产一级毛片高清牌| 欧美在线黄色| 亚洲国产欧美在线一区| 国产精品久久久人人做人人爽| 女人精品久久久久毛片| 国产精品久久久av美女十八| 亚洲精品国产一区二区精华液| 欧美精品高潮呻吟av久久| 可以免费在线观看a视频的电影网站| av线在线观看网站| 两人在一起打扑克的视频| 12—13女人毛片做爰片一| 久久久国产成人免费| 精品国产超薄肉色丝袜足j| 波多野结衣一区麻豆| 黄网站色视频无遮挡免费观看| 欧美日韩成人在线一区二区| 欧美成人午夜精品| 别揉我奶头~嗯~啊~动态视频 | 免费在线观看黄色视频的| 亚洲欧美激情在线| 国产av又大| 国产精品免费视频内射| 叶爱在线成人免费视频播放| 欧美日韩一级在线毛片| 99九九在线精品视频| 欧美激情 高清一区二区三区| 视频区图区小说| 久久人人97超碰香蕉20202| av福利片在线| 日本91视频免费播放| 成人国产av品久久久| a 毛片基地| 成人免费观看视频高清| 真人做人爱边吃奶动态| 久久人妻福利社区极品人妻图片| 人妻人人澡人人爽人人| 91九色精品人成在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 大陆偷拍与自拍| netflix在线观看网站| 日韩欧美免费精品| 亚洲av片天天在线观看| 国产福利在线免费观看视频| 亚洲精品成人av观看孕妇| 亚洲一区二区三区欧美精品| 亚洲精品第二区| 亚洲激情五月婷婷啪啪| 欧美 日韩 精品 国产| 乱人伦中国视频| 亚洲色图 男人天堂 中文字幕| 丝袜脚勾引网站| 亚洲精品日韩在线中文字幕| 久久人人爽人人片av| 成年人黄色毛片网站| 欧美日韩黄片免| 人人妻人人澡人人看| 欧美 日韩 精品 国产| 菩萨蛮人人尽说江南好唐韦庄| 黑人巨大精品欧美一区二区蜜桃| 亚洲综合色网址| 欧美在线黄色| 黄片播放在线免费| 国产成人精品在线电影| 一本综合久久免费| 亚洲国产欧美一区二区综合| 侵犯人妻中文字幕一二三四区| 亚洲天堂av无毛| 黄色 视频免费看| 久久九九热精品免费| 久久人人97超碰香蕉20202| 一进一出抽搐动态| 免费少妇av软件| 亚洲精品国产一区二区精华液| 亚洲国产精品999| 天堂8中文在线网| 国产精品一区二区在线观看99| 99久久99久久久精品蜜桃| 亚洲精品粉嫩美女一区| 啦啦啦视频在线资源免费观看| 人人妻人人澡人人爽人人夜夜| 99国产精品99久久久久| 精品国产乱码久久久久久男人| 午夜福利视频精品| 欧美激情高清一区二区三区| 一二三四在线观看免费中文在| 久久精品成人免费网站| 三上悠亚av全集在线观看| 国产精品秋霞免费鲁丝片| 欧美97在线视频| 91麻豆av在线| 在线精品无人区一区二区三| 女人爽到高潮嗷嗷叫在线视频| 欧美xxⅹ黑人| 动漫黄色视频在线观看| 99香蕉大伊视频| videos熟女内射| 国产男人的电影天堂91| 99国产精品免费福利视频| 中文字幕最新亚洲高清| 国产在线视频一区二区| 我要看黄色一级片免费的| 精品免费久久久久久久清纯 | 欧美日韩亚洲国产一区二区在线观看 | 大香蕉久久成人网| 男人爽女人下面视频在线观看| 视频在线观看一区二区三区| 国产精品一二三区在线看| 美女大奶头黄色视频| 午夜两性在线视频| 日本av免费视频播放| 久久久久久久国产电影| cao死你这个sao货| 国产日韩欧美视频二区| 最黄视频免费看| 国产xxxxx性猛交| 在线av久久热| 人妻一区二区av| 精品国产国语对白av| 国产免费福利视频在线观看| 婷婷色av中文字幕| 精品久久久久久电影网| 国产av一区二区精品久久| 欧美少妇被猛烈插入视频| 女人久久www免费人成看片| 亚洲专区字幕在线| 男女边摸边吃奶| 久久久久久久大尺度免费视频| 两人在一起打扑克的视频| 他把我摸到了高潮在线观看 | 国产日韩欧美视频二区| 国产欧美日韩一区二区精品| 天天躁狠狠躁夜夜躁狠狠躁| 夫妻午夜视频| 午夜福利在线免费观看网站| 高清欧美精品videossex| 成人手机av| 丰满饥渴人妻一区二区三| 国产成人av教育| 久久久精品免费免费高清| 9色porny在线观看| 性色av乱码一区二区三区2| 老司机影院成人| 久久久国产成人免费| 国产成人系列免费观看| 男女高潮啪啪啪动态图| 久久精品人人爽人人爽视色| 久久精品亚洲熟妇少妇任你| 国产熟女午夜一区二区三区| 啦啦啦免费观看视频1| 久久久精品94久久精品| 高清视频免费观看一区二区| 热99国产精品久久久久久7| 成人黄色视频免费在线看| 成人av一区二区三区在线看 | 欧美av亚洲av综合av国产av| 黄色a级毛片大全视频| 亚洲久久久国产精品| 老司机影院成人| 亚洲人成电影免费在线| 午夜影院在线不卡| 啦啦啦在线免费观看视频4| 国产91精品成人一区二区三区 | 国产区一区二久久| 国产亚洲一区二区精品| 欧美成狂野欧美在线观看| av天堂在线播放| 两个人看的免费小视频| 90打野战视频偷拍视频| 一边摸一边做爽爽视频免费| 后天国语完整版免费观看| 国产亚洲精品第一综合不卡| 国产成人av激情在线播放| 91成年电影在线观看| 侵犯人妻中文字幕一二三四区| 国产av又大| 久久精品熟女亚洲av麻豆精品| 一级毛片女人18水好多| 性色av乱码一区二区三区2| 12—13女人毛片做爰片一| 国产97色在线日韩免费| 一二三四社区在线视频社区8| 女性被躁到高潮视频| 如日韩欧美国产精品一区二区三区| 欧美av亚洲av综合av国产av| 中文字幕最新亚洲高清| 久久人妻熟女aⅴ| 秋霞在线观看毛片| 成人手机av| 男人添女人高潮全过程视频| 精品第一国产精品| 高清黄色对白视频在线免费看| 中国美女看黄片| 啦啦啦啦在线视频资源| 久久精品久久久久久噜噜老黄| videos熟女内射| 亚洲欧美精品综合一区二区三区| 午夜福利在线观看吧| 十八禁高潮呻吟视频| 两个人看的免费小视频| 国产男女超爽视频在线观看| 久久精品人人爽人人爽视色| 久热爱精品视频在线9| 色婷婷av一区二区三区视频| 亚洲国产精品999| 亚洲国产精品999| 午夜日韩欧美国产| 香蕉国产在线看| 亚洲少妇的诱惑av| 精品国内亚洲2022精品成人 | 国产三级黄色录像| 啦啦啦 在线观看视频| 女性被躁到高潮视频| 精品一区二区三区四区五区乱码| 久久精品国产亚洲av香蕉五月 | 18禁观看日本| 99热网站在线观看| 日本黄色日本黄色录像| 黄色视频在线播放观看不卡| 亚洲精品自拍成人| 性色av一级| 久久 成人 亚洲| 久久精品人人爽人人爽视色| 热re99久久国产66热| 国产精品久久久久久人妻精品电影 |