• <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| 婷婷色综合大香蕉| 国产综合精华液| 国产老妇伦熟女老妇高清| 亚洲天堂国产精品一区在线| 老司机影院毛片| 少妇被粗大猛烈的视频| 女人十人毛片免费观看3o分钟| 久久久久国产精品人妻一区二区| 日韩欧美一区视频在线观看 | 免费黄色在线免费观看| 欧美潮喷喷水| 亚洲精品乱码久久久久久按摩| 亚洲av福利一区| 亚洲综合精品二区| 日本一本二区三区精品| 国产女主播在线喷水免费视频网站| 成人国产av品久久久| 美女高潮的动态| a级一级毛片免费在线观看| 日本-黄色视频高清免费观看| 少妇人妻 视频| 啦啦啦中文免费视频观看日本| 啦啦啦啦在线视频资源| 日本熟妇午夜| 亚洲天堂av无毛| 欧美性感艳星| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 免费看a级黄色片| 国产亚洲最大av| 欧美xxxx性猛交bbbb| 久久精品熟女亚洲av麻豆精品| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| 狂野欧美激情性bbbbbb| av在线亚洲专区| 国产黄色免费在线视频| av一本久久久久| 51国产日韩欧美| 国产伦精品一区二区三区四那| 两个人的视频大全免费| 综合色av麻豆| 国产一级毛片在线| 国产白丝娇喘喷水9色精品| 精品人妻视频免费看| 国产精品女同一区二区软件| 亚洲av不卡在线观看| 十八禁网站网址无遮挡 | 美女脱内裤让男人舔精品视频| 精品久久久久久电影网| 久久久色成人| 少妇的逼水好多| 日韩av在线免费看完整版不卡| 精品久久国产蜜桃| 麻豆精品久久久久久蜜桃| kizo精华| 少妇人妻一区二区三区视频| 午夜亚洲福利在线播放| 亚洲美女搞黄在线观看| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 免费大片18禁| 国产精品久久久久久精品电影小说 | 高清在线视频一区二区三区| 亚洲色图av天堂| 日韩精品有码人妻一区| 国产在视频线精品| 午夜免费观看性视频| 国产永久视频网站| 午夜免费鲁丝| 久久精品国产亚洲av天美| 亚洲成色77777| 久久久久久伊人网av| 亚洲精品国产av成人精品| 国产 一区 欧美 日韩| 亚洲电影在线观看av| 亚洲av免费在线观看| 男人和女人高潮做爰伦理| 成人综合一区亚洲| 欧美亚洲 丝袜 人妻 在线| 18禁在线播放成人免费| 亚洲精品国产av蜜桃| 网址你懂的国产日韩在线| 久久99热这里只有精品18| 精品久久久精品久久久| 日本一本二区三区精品| 中文天堂在线官网| 综合色丁香网| 国产乱人偷精品视频| 亚洲电影在线观看av| 又粗又硬又长又爽又黄的视频| 下体分泌物呈黄色| 欧美激情国产日韩精品一区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av成人精品| 欧美97在线视频| 交换朋友夫妻互换小说| 少妇裸体淫交视频免费看高清| 国产精品一二三区在线看| 伊人久久精品亚洲午夜| 一本色道久久久久久精品综合| 99久国产av精品国产电影| 欧美潮喷喷水| 日韩中字成人| 国产亚洲精品久久久com| 夜夜爽夜夜爽视频| 亚洲激情五月婷婷啪啪| av在线播放精品| av线在线观看网站| 久久久亚洲精品成人影院| 欧美性猛交╳xxx乱大交人| 国产一级毛片在线| 国产成人精品久久久久久| 国产成人福利小说| 中文字幕免费在线视频6| 老司机影院毛片| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 精品酒店卫生间| 久久久久久伊人网av| 男人和女人高潮做爰伦理| 小蜜桃在线观看免费完整版高清| 自拍偷自拍亚洲精品老妇| 一区二区三区四区激情视频| av在线老鸭窝| 国产精品无大码| 男男h啪啪无遮挡| tube8黄色片| 日韩人妻高清精品专区| 一个人观看的视频www高清免费观看| 涩涩av久久男人的天堂| 国产女主播在线喷水免费视频网站| 国产高清国产精品国产三级 | 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 日韩亚洲欧美综合| 日日摸夜夜添夜夜添av毛片| 亚洲成人中文字幕在线播放| 一级毛片我不卡| 成人综合一区亚洲| 欧美bdsm另类| 熟女人妻精品中文字幕| 亚洲色图av天堂| 亚洲av福利一区| 青春草国产在线视频| 99久久九九国产精品国产免费| 国产精品一及| 成年免费大片在线观看| 亚洲经典国产精华液单| 舔av片在线| 午夜精品国产一区二区电影 | 女人十人毛片免费观看3o分钟| 国产综合懂色| av播播在线观看一区| 大片免费播放器 马上看| 欧美一级a爱片免费观看看| 日产精品乱码卡一卡2卡三| 少妇猛男粗大的猛烈进出视频 | 国产高清国产精品国产三级 | 五月开心婷婷网| 18禁裸乳无遮挡动漫免费视频 | 亚洲av成人精品一二三区| 老司机影院成人| 熟女av电影| 又大又黄又爽视频免费| 99热6这里只有精品| 成人高潮视频无遮挡免费网站| 中国三级夫妇交换| 久久精品国产亚洲网站| 久久久久久久久久人人人人人人| 国产精品女同一区二区软件| 男人添女人高潮全过程视频| 久久99精品国语久久久| 97热精品久久久久久| 男女无遮挡免费网站观看| 国产精品偷伦视频观看了| 欧美日韩国产mv在线观看视频 | 中文字幕av成人在线电影| 国产av国产精品国产| 亚洲av成人精品一区久久| 欧美亚洲 丝袜 人妻 在线| videossex国产| 欧美潮喷喷水| 欧美日韩精品成人综合77777| 日本欧美国产在线视频| 成人亚洲精品一区在线观看 | 日韩av免费高清视频| 色视频在线一区二区三区| 亚洲精品aⅴ在线观看| 亚洲精品久久久久久婷婷小说| 国产精品秋霞免费鲁丝片| 精品少妇久久久久久888优播| 亚洲av在线观看美女高潮| 亚洲成人av在线免费| 大码成人一级视频| 日韩 亚洲 欧美在线| 色5月婷婷丁香| 哪个播放器可以免费观看大片| 国产成人免费无遮挡视频| 欧美日韩国产mv在线观看视频 | 六月丁香七月| 乱码一卡2卡4卡精品| 女人久久www免费人成看片| 色吧在线观看| 精品人妻熟女av久视频| 国产成人一区二区在线| 久久久久精品久久久久真实原创| 色播亚洲综合网| 丝瓜视频免费看黄片| 少妇高潮的动态图| 高清日韩中文字幕在线| 亚洲欧美中文字幕日韩二区| 色视频www国产| 亚洲人成网站在线播| 亚洲电影在线观看av| 国产国拍精品亚洲av在线观看| 少妇熟女欧美另类| 成人高潮视频无遮挡免费网站| 成人国产麻豆网| 大香蕉97超碰在线| 一级毛片aaaaaa免费看小| 日日啪夜夜撸| 日本爱情动作片www.在线观看| 97在线人人人人妻| av在线老鸭窝| 免费高清在线观看视频在线观看| av在线播放精品| 国产91av在线免费观看| 大话2 男鬼变身卡| 亚洲av电影在线观看一区二区三区 | 久久久久久久亚洲中文字幕| 日本-黄色视频高清免费观看| 2018国产大陆天天弄谢| 欧美xxxx性猛交bbbb| 高清视频免费观看一区二区| 精品一区二区免费观看| 亚洲国产色片| 中文资源天堂在线| 日日啪夜夜撸| 日韩成人av中文字幕在线观看| 麻豆乱淫一区二区| 国产成人精品福利久久| 黑人高潮一二区| 狂野欧美白嫩少妇大欣赏| 国产成人freesex在线| 亚洲av男天堂| 下体分泌物呈黄色| 天天一区二区日本电影三级| 最近最新中文字幕免费大全7| 丝袜美腿在线中文| 狠狠精品人妻久久久久久综合| 精品国产露脸久久av麻豆| 久久久精品94久久精品| 成人鲁丝片一二三区免费| 色综合色国产| 3wmmmm亚洲av在线观看| 久久人人爽av亚洲精品天堂 | 欧美一区二区亚洲| 国产亚洲一区二区精品| 国产欧美亚洲国产| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区| 国产毛片在线视频| www.色视频.com| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 欧美一区二区亚洲| 在线观看一区二区三区激情| 日韩欧美精品v在线| 亚洲最大成人手机在线| 尾随美女入室| 七月丁香在线播放| 天堂俺去俺来也www色官网| 肉色欧美久久久久久久蜜桃 | 熟女av电影| 亚洲精品乱久久久久久| 国产精品久久久久久久久免| 欧美xxⅹ黑人| 欧美激情在线99| 久久这里有精品视频免费| 久久99热这里只有精品18| 三级国产精品片| 日本-黄色视频高清免费观看| 嫩草影院新地址| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄| 国产黄色视频一区二区在线观看| 免费av不卡在线播放| 久久精品综合一区二区三区| 丝瓜视频免费看黄片| 亚洲av中文av极速乱| 久久久久九九精品影院| 天天一区二区日本电影三级| 国产av不卡久久| 男人狂女人下面高潮的视频| 国产精品久久久久久久电影| 亚洲av电影在线观看一区二区三区 | 国产伦在线观看视频一区| 久久久久久伊人网av| 秋霞在线观看毛片| 久久女婷五月综合色啪小说 | a级毛色黄片| av福利片在线观看| 亚洲国产欧美人成| 久久精品国产亚洲av天美| 精品人妻偷拍中文字幕| 女人十人毛片免费观看3o分钟| 成人美女网站在线观看视频| 日本午夜av视频| 国产精品熟女久久久久浪| 国产伦在线观看视频一区| 爱豆传媒免费全集在线观看| 日韩 亚洲 欧美在线| 久久这里有精品视频免费| 少妇人妻一区二区三区视频| 五月天丁香电影| 99视频精品全部免费 在线| 毛片一级片免费看久久久久| 99久久九九国产精品国产免费| 午夜亚洲福利在线播放| 伊人久久精品亚洲午夜| 久久97久久精品| 欧美激情国产日韩精品一区| 国产黄片视频在线免费观看| 日韩国内少妇激情av| 欧美一区二区亚洲| 亚洲在久久综合| 久久人人爽av亚洲精品天堂 | 在线看a的网站| 一区二区三区乱码不卡18| 成人亚洲精品一区在线观看 | 久久久久九九精品影院| 肉色欧美久久久久久久蜜桃 | 久久精品人妻少妇| 国产片特级美女逼逼视频| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| 99视频精品全部免费 在线| 精品久久久久久久末码| 大陆偷拍与自拍| 亚洲aⅴ乱码一区二区在线播放| 熟女电影av网| 大香蕉97超碰在线| 国产免费一级a男人的天堂| 男插女下体视频免费在线播放| 大陆偷拍与自拍| 国产免费一区二区三区四区乱码| 精品国产露脸久久av麻豆| 亚洲国产最新在线播放| 九九在线视频观看精品| 禁无遮挡网站| 人妻夜夜爽99麻豆av| 亚洲最大成人手机在线| 欧美精品一区二区大全| 少妇人妻精品综合一区二区| 亚洲精品,欧美精品| 久久久久性生活片| 欧美日韩视频精品一区| 成年av动漫网址| 最新中文字幕久久久久| 亚洲国产精品国产精品| 高清日韩中文字幕在线| 成人美女网站在线观看视频| 国产黄a三级三级三级人| 男人添女人高潮全过程视频| 久久久精品94久久精品| 在线播放无遮挡| 国产亚洲最大av| 一级毛片aaaaaa免费看小| 少妇人妻久久综合中文| 嘟嘟电影网在线观看| 中文资源天堂在线| 亚洲精华国产精华液的使用体验| 欧美日韩在线观看h| 亚洲高清免费不卡视频| 久热这里只有精品99| 五月伊人婷婷丁香| 丝袜脚勾引网站| 九草在线视频观看| 男女边吃奶边做爰视频| 日本三级黄在线观看| 久久99蜜桃精品久久| 国产亚洲精品久久久com| 久久久久国产网址| 色哟哟·www| 别揉我奶头 嗯啊视频| 亚洲av中文av极速乱| 男插女下体视频免费在线播放| 内射极品少妇av片p| 麻豆国产97在线/欧美| 免费黄网站久久成人精品| 日韩av不卡免费在线播放| 日韩欧美精品免费久久| 国产探花极品一区二区| 97超视频在线观看视频| 晚上一个人看的免费电影| 2022亚洲国产成人精品| 在线免费十八禁| 国产男女超爽视频在线观看| 偷拍熟女少妇极品色| 成年女人在线观看亚洲视频 | 好男人在线观看高清免费视频| 欧美最新免费一区二区三区| 久久99热这里只频精品6学生| 国产精品国产av在线观看| 精品午夜福利在线看| 精品亚洲乱码少妇综合久久| 欧美日韩视频精品一区| a级毛色黄片| 只有这里有精品99| 丝袜脚勾引网站| 国产精品久久久久久久电影| 国产精品嫩草影院av在线观看| 纵有疾风起免费观看全集完整版| 三级国产精品片| 97人妻精品一区二区三区麻豆| 蜜桃久久精品国产亚洲av| 日韩精品有码人妻一区| 亚洲天堂国产精品一区在线| 18禁裸乳无遮挡动漫免费视频 | 欧美性感艳星| 三级经典国产精品| 日韩亚洲欧美综合| 亚洲真实伦在线观看| 国产精品久久久久久精品古装| 少妇丰满av| 乱码一卡2卡4卡精品| 又粗又硬又长又爽又黄的视频| 波多野结衣巨乳人妻| 久久人人爽av亚洲精品天堂 | 五月伊人婷婷丁香| 在线观看人妻少妇| 一边亲一边摸免费视频| 伦精品一区二区三区| 日日摸夜夜添夜夜爱| 欧美区成人在线视频| 亚洲精品日本国产第一区| 国产黄色免费在线视频| 国产91av在线免费观看| 大码成人一级视频| 日韩成人伦理影院| 性插视频无遮挡在线免费观看| 全区人妻精品视频| 国产黄色视频一区二区在线观看| 亚洲精品日韩av片在线观看| 极品少妇高潮喷水抽搐| 国产一区二区亚洲精品在线观看| 伊人久久精品亚洲午夜| 日韩免费高清中文字幕av| 日韩欧美 国产精品| 午夜激情福利司机影院| 日本午夜av视频| av在线app专区| 高清av免费在线| 亚洲国产精品999| 久久韩国三级中文字幕| 亚洲av日韩在线播放| 亚洲精品成人av观看孕妇| 美女主播在线视频| 久久久久久久亚洲中文字幕| 我要看日韩黄色一级片| 黄色日韩在线| 黄色视频在线播放观看不卡| 五月伊人婷婷丁香| 99热这里只有是精品50| 我要看日韩黄色一级片| av福利片在线观看| 日韩 亚洲 欧美在线| 只有这里有精品99| 欧美极品一区二区三区四区| 大码成人一级视频| av在线亚洲专区| 亚洲精品影视一区二区三区av| 日本-黄色视频高清免费观看| 国产中年淑女户外野战色| 国产 一区精品| 尤物成人国产欧美一区二区三区| 久久久久久国产a免费观看| 伊人久久精品亚洲午夜| 大香蕉97超碰在线| av播播在线观看一区| 高清在线视频一区二区三区| 久久精品综合一区二区三区| 五月开心婷婷网| 老司机影院成人| 国产精品成人在线| 高清在线视频一区二区三区| 91久久精品国产一区二区三区| 国产毛片a区久久久久| 国内精品宾馆在线| 国产精品无大码| 极品少妇高潮喷水抽搐| 噜噜噜噜噜久久久久久91| 久久久久久久久久久免费av| 男人爽女人下面视频在线观看| 国产一区二区三区av在线| 涩涩av久久男人的天堂| 噜噜噜噜噜久久久久久91| 亚洲成色77777| 国产 精品1| 看非洲黑人一级黄片| 国产一区有黄有色的免费视频| 狂野欧美激情性xxxx在线观看| 国产白丝娇喘喷水9色精品| 亚洲av免费高清在线观看| 国产午夜精品久久久久久一区二区三区| 久久久久久久久久久丰满| 国国产精品蜜臀av免费| h日本视频在线播放| 人妻少妇偷人精品九色| 美女主播在线视频| 日韩电影二区| 亚洲aⅴ乱码一区二区在线播放| 我的老师免费观看完整版| 又爽又黄无遮挡网站| 高清日韩中文字幕在线| 中国美白少妇内射xxxbb| 日韩欧美精品免费久久| 只有这里有精品99| 久久久a久久爽久久v久久| 久久影院123| 九九在线视频观看精品| 99热这里只有是精品在线观看| 亚洲伊人久久精品综合| 国产成人午夜福利电影在线观看| videossex国产| 2022亚洲国产成人精品| 乱码一卡2卡4卡精品| 一级二级三级毛片免费看| 插逼视频在线观看| 啦啦啦啦在线视频资源| 国产精品人妻久久久影院| 精品少妇黑人巨大在线播放| 男插女下体视频免费在线播放| 亚洲欧美精品专区久久| 国产日韩欧美在线精品| 日韩成人av中文字幕在线观看| 国产午夜福利久久久久久| 日本黄色片子视频| 亚洲av免费在线观看| 久久精品国产亚洲av涩爱| 日韩亚洲欧美综合| 欧美bdsm另类| 波多野结衣巨乳人妻| 久久午夜福利片| 蜜臀久久99精品久久宅男| 国产综合懂色| 校园人妻丝袜中文字幕| 身体一侧抽搐| 久久久久久久精品精品| 纵有疾风起免费观看全集完整版| 亚洲人成网站在线播| 久久久久久久午夜电影| 人人妻人人澡人人爽人人夜夜| 久久人人爽人人片av| 一级毛片 在线播放| 欧美+日韩+精品| 老司机影院成人| 26uuu在线亚洲综合色| 成人午夜精彩视频在线观看| av国产久精品久网站免费入址| 久久久久久九九精品二区国产| 国产久久久一区二区三区| 噜噜噜噜噜久久久久久91| 欧美日韩一区二区视频在线观看视频在线 | 天堂网av新在线| 国产精品精品国产色婷婷| 成年女人看的毛片在线观看| 欧美日韩亚洲高清精品| 一级毛片电影观看| 插逼视频在线观看| 亚洲精品456在线播放app| 国产成人福利小说| 国产高清不卡午夜福利| 免费看av在线观看网站| 日本猛色少妇xxxxx猛交久久| 波野结衣二区三区在线| 久久精品久久久久久久性| 日韩三级伦理在线观看| 国产毛片a区久久久久| 色婷婷久久久亚洲欧美| 免费看a级黄色片| 亚洲精品乱码久久久v下载方式| 777米奇影视久久| 国产又色又爽无遮挡免| av国产精品久久久久影院| 亚洲在久久综合| 国产免费一区二区三区四区乱码| 高清午夜精品一区二区三区| 国产成人aa在线观看| 五月伊人婷婷丁香| 国产亚洲av片在线观看秒播厂| 国产成人福利小说| 搞女人的毛片| 成人无遮挡网站| 亚洲精品aⅴ在线观看| 国产亚洲一区二区精品| 亚洲国产av新网站| 成年女人在线观看亚洲视频 | 午夜免费男女啪啪视频观看| 国精品久久久久久国模美| 全区人妻精品视频| 秋霞在线观看毛片| 日韩伦理黄色片| 又黄又爽又刺激的免费视频.| 久久ye,这里只有精品| 成年版毛片免费区| 老司机影院成人| 国产精品无大码| av在线app专区| 噜噜噜噜噜久久久久久91| 看十八女毛片水多多多|