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

    2000—2014年黃土高原植被葉面積指數(shù)時空變化特征

    2017-04-18 07:09:08王志慧姚文藝湯秋鴻楊二孔祥兵王玲玲肖培青
    中國水土保持科學 2017年1期
    關鍵詞:黃土高原覆蓋度變化率

    王志慧,姚文藝,湯秋鴻,楊二,孔祥兵,王玲玲,肖培青

    (1.黃河水利科學研究院水利部黃土高原水土流失過程與控制重點實驗室,450003,鄭州; 2.中國科學院地理科學與資源研究所,100101,北京)

    2000—2014年黃土高原植被葉面積指數(shù)時空變化特征

    王志慧1,2,姚文藝1?,湯秋鴻2,楊二1,孔祥兵1,王玲玲1,肖培青1

    (1.黃河水利科學研究院水利部黃土高原水土流失過程與控制重點實驗室,450003,鄭州; 2.中國科學院地理科學與資源研究所,100101,北京)

    研究植被葉面積指數(shù)(LAI)時空變化特征,對植被的水土保持效具有重要意義。利用MOD15A2H遙感產(chǎn)品,基于Mann-Kendall趨勢檢驗與Sen斜率分析方法,提取區(qū)域尺度與像素尺度上的植被LAI變化特征,并基于不同子流域、坡度、坡向及植被覆蓋類型,對植被LAI的變化特征進行分析。基于MOD44B遙感產(chǎn)品,利用線性回歸和偏相關系數(shù),分析植被LAI的變化原因。結果表明:1)黃土高原2000—2014年,植被LAI呈顯著增加趨勢,其年絕對變化幅度為0.042,年相對變化程度為2.71%。2)空間上,在黃土高原58.6%的區(qū)域,LAI呈現(xiàn)顯著增加趨勢,僅有0.9%的區(qū)域LAI呈現(xiàn)顯著減少趨勢。植被LAI劇烈增加,主要發(fā)生在河口—龍門區(qū)間,包括皇甫川、窟野河、無定河和延河。植被在15°~35°的坡度上,LAI變化程度最劇烈,其變化在各坡向上沒有顯著差異,農(nóng)田和草地的LAI變化程度最劇烈。3)與植被總覆蓋度相比,植被垂直維結構與黃土高原植被LAI的變化更為相關,其中樹木覆蓋度的增加,是植被垂直維結構變化的重要原因之一。

    MODIS產(chǎn)品;黃土高原;植被葉面積指數(shù);時空變化分析

    黃土高原是我國乃至世界上水土流失最嚴重的地區(qū)。黃河上中游地區(qū)人類活動對自然植被的破壞是造成該生態(tài)問題的主要原因,尤其是對坡地植被的破壞,大大增加水土流失和洪澇災害的風險[1-4]。為了解決這一生態(tài)環(huán)境問題,中國政府于1999年啟動“退耕還林還草”重點生態(tài)修復工程[5]。目前大量研究已表明:近10年內,人類活動對黃土高原下墊面影響巨大,植被覆蓋明顯增加[6-8];同時,對黃土高原的產(chǎn)流產(chǎn)沙也產(chǎn)生了顯著的影響[9-10]。但是,也有相關研究指出,該生態(tài)恢復工程的效益被夸大,且存在一定的負影響[11-13]。因此,對近十幾年黃土高原的植被變化進行時空變化分析,有助于全面了解該區(qū)域植被變化面積與變化程度,并對科學評價流域生態(tài)環(huán)境和制定有效的水土保持措施,均具有重要意義。

    大量研究表明,植被的結構特征是影響植被水土保持效益的重要因素[14-15]。葉面積指數(shù)(Leaf Area Index,LAI)是最重要的植被結構參數(shù)之一,也是陸地生態(tài)、水文模型模擬中不可或缺的植被參數(shù)。LAI是指單位土地面積上,植物葉片總面積占土地面積的倍數(shù),該參數(shù)能夠同時反映植被水平與垂直結構特征。目前,遙感技術已經(jīng)成為監(jiān)測區(qū)域尺度植被動態(tài)變化的最有效工具。NDVI數(shù)據(jù)是進行區(qū)域植被監(jiān)測的常用數(shù)據(jù)源,如AVHRR[1617]、 SPOT[18]和MODIS[7,19];但NDVI不能直接代表特定植被結構參數(shù),且研究中所選遙感數(shù)據(jù)空間尺度較大,難以反映地表植被的真實變化情況。目前,針對黃土高原植被LAI的時空變化分析研究還相對較少,筆者首次利用NASA最新發(fā)布的高分辨率LAI產(chǎn)品,從區(qū)域到像素尺度,對黃土高原植被LAI的時空變化特征進行提取,從不同角度,全面了解2000—2014年黃土高原植被LAI的變化情況,同時對植被LAI的變化原因進行分析。

    1 研究區(qū)概況

    黃河是中國第二長河,是世界第五長河流,也是世界上含沙量最多的河流。黃土高原(E 100.9°~114.5°,N 33.7°~41.3°)是世界上水土流失最嚴重和生態(tài)環(huán)境最脆弱的地區(qū)之一,也是黃河泥沙主要來源區(qū),黃土高原海拔一般在1 000~1 300 m,地貌起伏不平,坡陡溝深,千溝萬壑,溝谷面積占40%~50%。黃土高原主要產(chǎn)流產(chǎn)沙的8個子流域分別為皇甫川、窟野河、無定河、延河、北洛河、涇河、渭河和汾河,其地理位置見圖1。

    2 研究方法

    2.1 數(shù)據(jù)獲取與處理

    圖1 黃土高原與8個主要子流域地理位置Fig.1 Geographic location of the Loess Plateau and eight sub-basins

    2.1.1 遙感數(shù)據(jù)使用2015年NASA最新發(fā)布的MODIS LAI產(chǎn)品(MOD15A2H),對2000—2014年黃土高原植被LAI進行時空變化分析研究。其產(chǎn)品是基于地表反射率數(shù)據(jù)和土地覆蓋類型數(shù)據(jù),綜合利用經(jīng)驗模型和三維輻射傳輸模型查找表法,對植被真實LAI進行反演得到[20],其空間分辨率500 m,時間分辨率8 d。其他數(shù)據(jù)包括MOD44B產(chǎn)品、ASTER-DEM和MCD12Q1產(chǎn)品。

    2.1.2 遙感數(shù)據(jù)處理利用MODIS數(shù)據(jù)處理軟件MRT(MODIS Reprojection Tools),對MODIS產(chǎn)品進行數(shù)據(jù)預處理。MOD44B產(chǎn)品表示年內植被最大覆蓋度,將250 m分辨率的植被覆蓋度,上推到500 m分辨率,尺度上推方法采用聚合平均法。對于DEM數(shù)據(jù)處理,采用基于30 m分辨率的ASTER-DEM,計算得到局地坡度和坡向。

    對于MOD15A2H產(chǎn)品,由于在實際成像過程中,受到云與氣溶膠的影響,導致遙感影像中存在不可避免的噪聲信號,從而不能客觀的反映地表真實LAI情況。針對這一問題,筆者采用Chen Jin等[21]提出的迭代濾波方法,對LAI數(shù)據(jù)進行濾波處理,某一典型植樹造林區(qū)域像素的LAI時間序列重構前后對比情況見圖2。

    從圖2可知,由于受到云的影響,2003和2008年夏季LAI數(shù)據(jù)上下波動劇烈。這顯然不符合植被的生長規(guī)律,而重構算法對受云影響的LAI數(shù)據(jù)進行平滑,使得重構數(shù)據(jù)曲線更符合植被物候變化規(guī)律,從而可以大大提高數(shù)據(jù)分析精度。

    按照上述方法對MODIS產(chǎn)品處理后,得到2000—2014年時間分辨率為8 d的LAI數(shù)據(jù)序列。為了反映植被的年內最佳生長狀況,對每年的LAI數(shù)據(jù)序列中最大值進行提取,然后構建2000—2014年LAI年最大值數(shù)據(jù)序列。

    2.2 LAI時空變化特征

    不同區(qū)域尺度植被LAI變化發(fā)生的空間位置、變化范圍和變化強度都是植被變化監(jiān)測研究的重要內容。

    圖2 典型植樹造林區(qū)域2000—2014年LAI數(shù)據(jù)序列重構前后對比Fig.2 Comparison between original and reconstructed LAI time-series during 2000—2014 in a typical plantation region

    2.2.1 變化面積比例采用Mann-Kendall(M-K)非參數(shù)檢驗方法[22],對黃河流域植被LAI時間序列變化趨勢進行統(tǒng)計分析。若時間序列M-K檢驗顯著性水平α<0.05,則認為該趨勢變化為顯著性變化;若顯著性水平>0.05,則認為該趨勢沒有統(tǒng)計學意義,即表示沒有發(fā)生顯著性變化。

    LAI變化面積比例(D)是指某個區(qū)域范圍內, LAI顯著變化(P<0.05)區(qū)域占該區(qū)域范圍總面積的比重,該指標可反映LAI在某個區(qū)域內的變化劇烈程度。計算公式如下:

    式中:A為LAI顯著變化的區(qū)域面積,km2;S為某區(qū)域范圍的總面積,km2。

    2.2.2 絕對與相對變化率采用Kendall傾斜率β表示每年LAI絕對變化率,該指標可反映LAI的年絕對變化幅度,β>0表示上升趨勢,β<0表示下降趨勢,其變化趨勢顯著性水平參考M-K檢驗結果。計算公式如下:

    式中:xj和xi分別為第j年和第i年的觀測數(shù)值,j>i;n為年總數(shù)。

    R為每年LAI相對變化率,該指標可反映LAI的年相對變化程度,其計算公式如下:

    式中:β為LAI絕對變化率;xi為第i年的觀測數(shù)值; n為年總數(shù)。由于該值利用年均值對β進行標準化,從而使得不同空間像素的趨勢變化程度具有可比性[23]。

    2.2.3 區(qū)域/流域尺度LAI時空變化分析為了反映黃土高原以及其內部不同流域間LAI的變化差異,對黃土高原以及其內部流域面積較大,且產(chǎn)流產(chǎn)沙量較大的8個子流域的植被LAI進行時空變化分析。采用絕對變化率和相對變化率2個指標,對區(qū)域/流域尺度LAI變化特征進行評估。

    2.2.4 像素尺度LAI時空變化分析為了反映植被LAI變化特征在黃土高原、不同子流域、不同坡度和坡向以及不同植被覆蓋類型的差異性,采用變化面積比和相對變化率2個指標,對不同區(qū)域范圍內的植被LAI變化特征差異進行評估。其中:子流域為所選8個典型流域;坡度分類執(zhí)行水利部頒布的SL 190—2007《土壤侵蝕分類分級標準》[24],將坡度劃分為6個等級:<5°,5°~8°, 8°~15°,15°~25°,25°~35°,>35°;坡向分類按照45°夾角,由北開始,順時針方向旋轉分為9個坡向:平地(-1°)、北坡(N,0~22.5°,337.5°~360°)、東北坡(NE,22.5°~67.5°)、東(E,67.5°~112.5°)、東南坡(SE,112.5°~157.5°)、南坡(S, 157.5°~202.5°)、西南坡(SW,202.5°~247.5°)、西坡(W,247.5°~292.5°)、西北坡(NW,292.5°~337.5°)[25];植被覆蓋類型采用2001年的MOD12Q1土地覆蓋產(chǎn)品,提取其中的農(nóng)田、草地、灌木以及森林類型進行分析。

    2.3 LAI時空變化原因

    退耕還林還草政策,對黃土高原植被LAI的改善起到至關重要的作用,但不同區(qū)域植被LAI變化原因的空間分布規(guī)律并沒有被探討。由于在物理意義上,LAI同時反映植被水平覆蓋度與垂直維結構,為了分析不同維植被結構特征對植被LAI的影響,首先將植被垂直維結構特征定義為垂直維LAI,指單位植被覆蓋面積上,植物葉片總面積占植被覆蓋面積的倍數(shù)。其具體公式如下:

    式中:LAI為植被葉面積指數(shù);VFC為植被總覆蓋度;LAIver為垂直維LAI;Sleaf、Sveg和Sland分別為植物葉片總面積、植被覆蓋面積和土地面積。

    利用MOD44B產(chǎn)品,提取樹木覆蓋度與草本植被覆蓋度,二者相加即為植被總覆蓋度,并依據(jù)公式(4),計算得到垂直維LAI。首先將區(qū)域尺度LAI與4種植被參數(shù)(植被總覆蓋度、垂直維LAI、樹木覆蓋度、草本植被覆蓋度)分別進行線性回歸分析。為了更進一步分析植被類型變化對LAI的影響,根據(jù)公式(5)和(6),分別計算像素尺度LAI與樹木覆蓋度和草本植被覆蓋度的偏相關系數(shù),并對偏相關系數(shù)進行顯著性檢驗。通過偏相關分析空間分布圖和LAI趨勢變化空間分布圖的疊加分析,同時確定每個像素LAI的變化趨勢和與LAI變化顯著相關的變量,從而得到植被LAI顯著變化區(qū)域的變化原因空間分布圖。

    式中:r12,3為偏相關系數(shù);r12、r13、r23分別為不同要素之間的Pearson相關系數(shù);n為樣本數(shù),q為階數(shù);統(tǒng)計量t服從n-q-2個自由度的t分布。若t檢驗顯著性水平α<0.05,則認為2個要素之間存在顯著的相關性。

    3 結果與分析

    3.1 區(qū)域/流域尺度LAI時空變化特征

    2000—2014年,整個黃土高原以及8個典型子流域的植被LAI年最大值變化趨勢和變化特征統(tǒng)計結果見圖3和表1,整個黃土高原及其子流域的植被LAI均呈顯著上升趨勢(P<0.01)。整個黃土高原地區(qū),植被LAI的年絕對變化幅度為0.042,年相對變化程度為2.71%。在8個子流域中,汾河流域的LAI年絕對變化幅度最大,達到0.081,從大到小排列順序依次為汾河>延河>北洛河>涇河>渭河>無定河>皇甫川>窟野河。延河流域的LAI年相對變化程度最大,達到5.37%,從大到小排列順序依次為延河>皇甫川>無定河>汾河>窟野河>涇河>北洛河>渭河。

    表1 2000—2014年黃土高原與8個子流域植被結構參數(shù)平均值及其變化特征統(tǒng)計結果Tab.1 Statistical result of variation characteristics for the Loess Plateau and eight different sub-basins every year from 2000 to 2014

    圖3 2000—2014年黃土高原以及8個子流域的年最大LAI變化趨勢Fig.3 Change trend of annual maximum LAI for the Loess Plateau and 8 different sub-basins during 2000-2014

    圖4 像素尺度LAI絕對變化率與相對變化率空間分布圖Fig.4 Spatial distribution map of absolute change rate and relative change rate for LAI at pixel scale

    3.2 像素尺度LAI時空變化特征

    3.2.1 LAI變化總體情況利用公式(2)和(3),基于像素尺度計算每個像素的絕對變化率和相對變化率,計算結果見圖4。結果表明:絕對變化率呈現(xiàn)出由西北向東南方向遞增的趨勢,這與黃土高原氣候濕潤條件的空間梯度分布趨勢較為一致,說明在濕潤地區(qū),植被由于能夠獲取充足的水分,導致其生長較快, LAI絕對變化幅度也較大;相反,干旱地區(qū)植被由于缺水,導致其生長較慢,LAI絕對變化幅度也較小。因此,LAI絕對變化率在一定程度上反映了地區(qū)濕潤條件。反觀LAI相對變化率,消除了地區(qū)濕潤條件對植被變化的影響,其空間分布并沒有隨著地區(qū)濕潤條件的變化而變化,可真實表達LAI變化程度,且不同區(qū)域的LAI變化程度可直接進行比較。

    利用公式(1)對整個黃土高原的變化面積比進行統(tǒng)計。結果表明:58.6%的區(qū)域,LAI呈現(xiàn)顯著增加趨勢,增加程度最大的區(qū)域,主要集中在河口-龍門區(qū)間;僅有0.9%的區(qū)域,LAI呈現(xiàn)顯著減少趨勢,主要分布在渭河、伊洛河和湟水流域零星區(qū)域;其余40.5%的區(qū)域,LAI保持穩(wěn)定狀態(tài)。LAI變化趨勢空間分布見圖4。

    3.2.2 LAI變化的流域差異LAI變化的流域差異如圖5所示。結果表明:根據(jù)變化面積比和相對變化率的排序比較,延河流域LAI變化劇烈程度最大,渭河流域LAI變化劇烈程度最小。除渭河流域外,其余流域內變化面積比均超過69%。

    圖5 8個子流域植被LAI變化面積比與相對變化率統(tǒng)計結果Fig.5 Statistics of changed area percentage and relative change rate of LAI in 8 sub-basins

    3.2.3 LAI變化的坡度差異將LAI重采樣到30 m尺度,與30 m尺度坡度數(shù)據(jù)進行疊加分析。LAI變化的坡度差異見圖6。結果表明:15°~25°和25°~35°2個坡度的植被LAI變化程度最為劇烈,變化面積比均超過60%,年相對變化率均超過3%;>35°區(qū)域的植被LAI變化程度最小,可能是由于坡度較大,導致水分較易流失,植被生長較為緩慢。

    圖6 黃土高原不同坡度區(qū)域植被LAI變化面積比與相對變化率統(tǒng)計結果Fig.6 Statistics of changed area percentage and relative change rate of LAI for different regions with different slopes in the Loess Plateau

    3.2.4 LAI變化的坡向差異將LAI重采樣到30 m尺度,與30 m尺度坡度數(shù)據(jù)進行疊加分析。LAI變化的坡向差異見圖7。結果表明,不同坡向植被LAI的變化面積比和相對變化率基本一致。說明黃土高原LAI的變化,在坡向上沒有顯著差異,在一定程度上反映出植被變化并沒有明顯的坡向規(guī)律。

    3.2.5 LAI變化的植被覆蓋差異LAI變化的植被覆蓋差異見圖8。結果表明,不同植被類型的LAI變化劇烈程度從大到小排序為農(nóng)田>草地>灌木>森林。這也反映出退耕還林還草政策的實施效果,使得原來的農(nóng)田和草地區(qū)域植被LAI變化劇烈。

    3.3 LAI時空變化原因分析

    區(qū)域尺度LAI與4種植被參數(shù)的線性回歸分析結果見圖9。從圖中可知,對于整個黃土高原植被總覆蓋度與植被LAI相關性較低(R2=0.21),垂直維LAI與植被LAI之間存在顯著的正相關性(R2=0.94),樹木覆蓋度與植被LAI也存在顯著正相關性(R2=0.82),說明黃土高原植被LAI的變化與植被垂直維結構的變化更為相關,而樹木覆蓋度的持續(xù)增加,則是植被垂直維結構變化的重要原因之一。

    圖7 黃土高原不同坡向區(qū)域植被LAI變化面積比與相對變化率統(tǒng)計結果Fig.7 Statistics of changed area percentage and relative change rate of LAI for different regions with different aspects in the Loess Plateau

    圖8 黃土高原不同植被類型LAI變化面積比與相對變化率統(tǒng)計結果Fig.8 Statistics of changed area percentage and relative change rate of LAI for different vegetation types in the Loess Plateau

    通過偏相關分析可知,某地區(qū)的植被LAI變化主要是由木本植被或草本植被的改變所致。將像素尺度LAI與樹木覆蓋度和草本植被覆蓋度偏相關分析空間分布圖,與LAI趨勢變化空間分布圖進行疊加分析見圖10。結果表明,在LAI顯著增加區(qū)域,42%的區(qū)域是由樹木覆蓋度的增加所致,4.3%的區(qū)域是由草本植被覆蓋度的增加所致,7.5%的區(qū)域是由樹木和草本植被覆蓋度共同增加所致。在LAI顯著減少區(qū)域中,9.5%的區(qū)域是由樹木覆蓋度的減少所致,16.8%的區(qū)域是由草本植被覆蓋度減少所致,4.9%的區(qū)域是由樹木和草本覆蓋度共同減少所致。

    4 結論

    1)黃土高原植被LAI的年絕對變化幅度為0.042,年相對變化程度為2.71%。在8個子流域中,汾河流域的LAI年絕對變化幅度最大,達到0.081,從大到小排列順序為汾河>延河>北洛河>涇河>渭河>無定河>皇甫川>窟野河。延河流域的LAI年相對變化程度最大,達到5.37%,從大到小排列順序為延河>皇甫川>無定河>汾河>窟野河>涇河>北洛河>渭河。

    2)在空間上,黃土高原58.6%的區(qū)域LAI呈現(xiàn)顯著增加趨勢,僅有0.9%的區(qū)域LAI呈現(xiàn)顯著減少趨勢,其余40.5%的區(qū)域LAI保持穩(wěn)定狀態(tài)。延河流域LAI的變化面積比最大,達到97%,而渭河流域LAI的變化面積比最小,僅為47.5%。15°~25°和25°~35°2個坡度的植被LAI變化程度最為劇烈,變化面積比均超過60%,年相對變化率均超過3%。黃土高原LAI的變化在坡向上沒有顯著差異。不同植被類型的LAI變化劇烈程度,從大到小排序為農(nóng)田>草地>灌木>森林。

    3)對于整個黃土高原與植被總覆蓋度相比,植被垂直維結構與植被LAI的變化更為相關;其中,樹木覆蓋度的持續(xù)增加是植被垂直維結構變化的重要原因之一。在LAI顯著增加區(qū)域中, 42%的區(qū)域是由樹木覆蓋度的增加所致,4.3%的區(qū)域是由草本植被覆蓋度的增加所致,7.5%的區(qū)域是由樹木和草本植被覆蓋度共同增加所致。在LAI顯著減少區(qū)域中,9.5%的區(qū)域是由樹木覆蓋度的減少所致,16.8%的區(qū)域是由草本植被覆蓋度減少所致,4.9%的區(qū)域是由樹木和草本覆蓋度共同減少所致。

    5 討論

    圖9 4種植被覆蓋度與平均植被LAI線性回歸分析Fig.9 Linear regression analysis of average LAI and coverage by 4 types of vegetation over the Loess Plateau

    圖10 像素尺度植被LAI的偏相關分析與LAI顯著變化趨勢疊加分析結果Fig.10 Overlay analysis of result derived from partial relation analysis and significant change trend of LAI at pixel scale

    在已有研究報道中,均使用低分辨率NDVI或植被覆蓋度數(shù)據(jù)對黃土高原植被動態(tài)變化進行分析,但NDVI在物理意義上不能直接代表植被結構參數(shù)。筆者選用500 m空間分辨率的LAI數(shù)據(jù)進行分析,與MODIS 1km[26]、LTDR 5km[7]、GIMMS 8km[27]NDVI數(shù)據(jù)的分析結果相比,能夠更為精細的反映黃土高原植被垂直結構的時空變化特征。另外,與已有成果相比,本研究利用相對變化率指標,對黃土高原植被變化劇烈程度的空間異質性規(guī)律進行分析和總結。在此需要指出的是,MODIS LAI數(shù)據(jù)是全球尺度產(chǎn)品,對于黃土高原這一特定區(qū)域必然存在絕對誤差,但LAI年際相對變化趨勢仍具有一定的可信度。今后可利用多尺度和多源遙感數(shù)據(jù),對黃土高原植被LAI變化進行多尺度分析研究。

    [1] CHEN Liding,WEI Wei,FU Bojie,et al.Soil and water conservation on the Loess Plateau in China:review and perspective[J].Progress in Physical Geography,2007, 31(4):389.

    [2] XIN Zhongbao,YU Xinxiao,LU Xixi.Factors controlling sediment yield in China's Loess Plateau[J].Earth Sur-face Processes&Landforms,2011,36(6):816.

    [3] FU Bojie,MENG Q H,Qiu Y,et al.Effects of land use on soil erosion and nitrogen loss in the hilly area of the Loess Plateau,China[J].Land Degradation&Development,2003,15(1):87.

    [4] YANG Xiaohui,ZHANG Ke,JIA Baoquan,et al.Desertification assessment in China:an overview[J].Journal of Arid Environments,2005,63(2):517.

    [5] LYU Yihe,FU Bojie,FENG Xiaoming.A policy-driven large scale ecological restoration:quantifying ecosystem services changes in the Loess Plateau of China[J].Plos One,2012,7(2):e31782.

    [6] LIU Yu,FU Bojie,LYU Yihe,et al.Hydrological responses and soil erosion potential of abandoned cropland in the Loess Plateau,China[J].Geomorphology,2012, 138(1):404.

    [7] SUN Wenyi,SONG Xiaoyan,MU Xingmin,et al.Spatiotemporal vegetation cover variations associated with climate change and ecological restoration in the Loess Plateau[J].Agricultural&Forest Meteorology,2015,209: 87.

    [8] SUN Wenyi,SHAO Quanqin,LIU Jiyuan.Soil erosion and its response to the changes of precipitation and vegetation cover on the Loess Plateau[J].Journal of Geographical Sciences,2013,23(6):1091.

    [9] 劉曉燕,劉昌明,楊勝天,等.基于遙感的黃土高原林草植被變化對河川徑流的影響分析[J].地理學報, 2014,69(11):1595. LIU Xiaoyan,LIU Changmin,YANG Shengtian,et al. Influences of shrubs-herbs-arbor vegetation coverage on the runoff based on the remote sensing data in Loess Plateau[J].Acta Geographica Sinica,2014,69(11):1595.

    [10]劉曉燕,楊勝天,王富貴,等.黃土高原現(xiàn)狀梯田和林草植被的減沙作用分析[J].水利學報,2014,45 (11):1293. LIU Xiaoyan,YANG Shengtian,WANG Fugue,et al. Sediment reduction effect of level terrace in the hilly-gully region in the Loess Plateau[J].Journal of Hydraulic Engineering,2014,45(11):1293.

    [11]JIAO Juying,ZHANG Zhenguo,BAI Wenjuan,et al. Assessing the ecological success of restoration by afforestation on the Chinese Loess Plateau[J].Restoration Ecology,2012,20(2):240.

    [12]ZHANG Qiang,XU Chongyu,YANG Tao.Variability of water resource in the Yellow River basin of past 50 Years,China[J].Water Resources Management,2009, 23(6):1157.

    [13]WANG Yunqiang,SHAO Mingan,LIU Zhipeng.Vertical distribution and influencing factors of soil water content within 21-m profile on the Chinese Loess Plateau[J]. Geoderma,2013,193(2):300.

    [14]孫佳佳,于東升,史學正,等.植被葉面積指數(shù)與覆蓋度定量表征紅壤區(qū)土壤侵蝕關系的對比研究[J].土壤學報,2010,47(6):1060. SUN Jiajia,YU Dongshen,SHI Xuezheng,et al.Comparison of between LAI and VFC in relationship with soil erosion in the red soil hilly region of South China[J]. Acta Pedologica Sinica.2010,47(6):1060.

    [15]馮強,趙文武.USLE/RUSLE中植被覆蓋與管理因子研究進展[J].生態(tài)學報,2014,34(16):4461. FENG Qiang,ZHAO Wenwu.The study on cover-management factor in USLE and RUSLE:a review[J].Acta Ecologica Sinica,2014,34(16):4461.

    [16]WANG Yanfang,SHEN Yanjun,CHEN Yaning,et al. Vegetation dynamics and their response to hydroclimatic factors in the Tarim River Basin,China[J].Ecohydrology,2013,6(6):927.

    [17]楊勝天,劉昌明,孫睿.近20年來黃河流域植被覆蓋變化分析[J].地理學報,2002,57(6):679. YANG Shengtian,LIU Changming,SUN Rui.The vegetation cover over last 20 years in Yellow River Basin[J]. Acta Geographica Sinica,2002,57(6):679.

    [18]CAO Lei,XU Jianhua,CHEN Yaning,et al.Understanding the dynamic coupling between vegetation cover and climatic factors in a semiarid region—a case study of Inner Mongolia,China[J].Ecohydrology,2013,6(6): 917.

    [19]WANG Jing,WANG Kelin,ZHANG Mingyang,et al. Impacts of climate change and human activities on vegetation cover in hilly southern China[J].Ecological Engineering,2015,81:451.

    [20]MYNENI R,KNYAZIKHINY,PARKT,etal. MOD15A2H MODIS/Terra Leaf Area Index/FPAR 8-Day L4 Global 500m SIN Grid V006[R/OL].http:∥doi. org/10.5067/MODIS/MOD15A2H.006.

    [21]CHEN Jin,JONSSON P,TAMURA M,et al.A simple method for reconstructing a high-quality NDVI time-series data set based on the Savitzky-Golay filter[J].Remote Sensing of Environment,2004,91(3):332.

    [22]李占玲,徐宗學.黑河上游山區(qū)徑流變化特征分析[J].干旱區(qū)資源與環(huán)境,2012,26(9):51. LI Zhanlin,XU Zhongxue.Characteristics of runoff time series for the upper reaches of Heihe River basin[J]. Journal of Arid Land Resources and Environment,2012, 26(9):51.

    [23]TANG Qiuhong,OKIL Taikan,KANAEL Shinjiro,et al.A spatial analysis of hydro-climatic and vegetation condition trends in the Yellow River basin[J].Hydrological Process,2008,22(3):451.

    [24]QUAN Bin,R?MKENS M J M,LI Rui,et al.Effect of land use and land cover change on soil erosion and the spatio-temporal variation in Liupan Mountain Region, southern Ningxia,China[J].Frontiers of Environmental Science&Engineering,2011,5(4):564

    [25]羅婭,楊勝天,劉曉燕,等.黃河河口鎮(zhèn)—潼關區(qū)間1998—2010年土地利用變化特征[J].地理學報, 2014,69(1):42. LUO Ya,YANG Shengtian,LIU Xiaoyan,et al.Land use change in the reach from Hekouzhen to Tongguan of the Yellow River during 1998-2010[J].Acta Geographica Sinica,2014,69(1):42.

    [26]劉英,岳輝.基于MODIS的2000—2014年黃土高原植被覆蓋度估算及其時空變化分析[J].環(huán)境與可持續(xù)發(fā)展,2016,41(3):181 LIU Ying,YUE Hui.Estimation and temporal variation analysis of vegetation coverage in the Loess Plateau in 2000-2014 based on MODIS[J].Environment and Sustainable Development,2016,41(3):181.

    [27]木熱提江·阿不拉,張曉萍,陳利利,等.基于GIMMS NDVI的黃土高原地區(qū)荒漠化時空特征分析[J].中國水土保持科學,2015,13(3):24. MURAT Abla,ZHANG Xiaoping,CHEN Lili,et al. Spatial-temporal distribution of desertification on the Loess Plateau using the GIMMS NDVI data[J].Science of Soil and Water Conservation,2015,13(3):24.

    Spatio-temporal change analysis of leaf area index in the Loess Plateau during 2000-2014

    WANG Zhihui1,2,YAO Wenyi1,TANG Qiuhong2,YANG Er1,KONG Xiangbing1,WANG Lingling1,XIAO Peiqing1
    (1.Yellow River Institute of Hydraulic Research,Yellow River Conservancy Commission,450003,Zhengzhou,China; 2.Institute of Geographic Sciences and Natural Resources Research,Chinese Academy of Sciences,100101,Beijing,China)

    [Background]Vegetation structural characteristic is an important factor determining soil and water conservation effect of vegetation,and leaf area index(LAI)is one of the most important vegetation structural characteristics.Remote sensing technology is the most effective tool to monitor LAI dynamics at regional scale.[Methods]In this paper,the spatio-temporal variation characteristics of LAI,including changed area percentage,absolute change rate and relative change rate,at regional scale and pixel scale were investigated respectively using Mann-Kendall trend analysis and Sen-slope analysis based on MOD15A2H product.And variation characteristics of LAI for different sub-basins,slopes,aspects and vegetation types were analyzed and compared.In addition,the reason of variation of LAI was exploredusing linear regression method and partial relation analysis based on MOD44B product.[Results]1)The change trend of LAI during 2000-2014 in the Loess Plateau increased significantly,and the absolute change rate and relative change rate of LAI were 0.042/Year and 2.71%/Year respectively.2)For the spatial distribution of variation characteristics of LAI,areas with a significant increase trend of LAI accounted for 58.6%in the Loess Plateau,and areas with a significant decrease trend of LAI accounted for 0.9%in the Loess Plateau.The variation extent of LAI in the regions between Hekou and Longmen, including Huangfuchuan,Kuyehe,Wudinghe and Yanhe were tremendous.The change rate and area percentage of changed regions for LAI in the area where slope was between 15°and 25°were the highest. There was no significant difference for LAI in different aspects.The change rate and area percentage of changed regions for LAI in cropland and grassland were the highest.3)In the Loess Plateau,the increase of tree coverage was the main reason for increase of vegetation LAI,and there were no significant relationship between LAI and vegetation fractional coverage.[Conclusions]The change trend of LAI during 2000-2014 in the Loess Plateau is significantly positive.Compared with absolute change rate, relative change rate is more able to illustrate the variation extent of LAI.According to statistical analysis, compared with total vegetation fractional coverage,LAI is more related to the vegetation vertical structure in the Loess Plateau,and the increase of tree coverage caused by“Grain to Green”project carried out across the Loess Plateau could be one of the important reasons for improvement of vegetation vertical structure.

    MODIS products;the Loess Plateau;LAI;spatio-temporal change analysis

    S157.1

    :A

    :2096-2673(2017)01-0071-10

    10.16843/j.sswc.2017.01.009

    2016- 05- 11

    2016- 10- 15

    項目名稱:黃河水利科學研究院科技發(fā)展基金項目“基于遙感數(shù)據(jù)的黃丘區(qū)林草覆蓋變化動態(tài)監(jiān)測方法研究”(201602);河南省科技攻關計劃項目“黃丘區(qū)降雨-侵蝕-植被響應臨界及模擬”(162101510004);國家自然科學基金“黃土丘陵區(qū)植被作用下產(chǎn)流機制及侵蝕動力響應”(41571276)

    王志慧(1985—),男,博士,工程師。主要研究方向:植被生態(tài)遙感與土壤侵蝕模擬與監(jiān)測。E-mail:wzh8588 @aliyun.com.cn

    ?通信作者簡介:姚文藝(1957—),男,博士,教授級高工。主要研究方向:土土壤侵蝕與水土保持、河流泥沙。E-mail:wyyao @yrihr.com.cn

    猜你喜歡
    黃土高原覆蓋度變化率
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
    基于NDVI的晉州市植被覆蓋信息提取
    基于電流變化率的交流濾波器失諧元件在線辨識方法
    湖南電力(2021年4期)2021-11-05 06:44:42
    例談中考題中的變化率問題
    低覆蓋度CO分子在Ni(110)面的吸附研究
    選舉 沸騰了黃土高原(下)
    公民與法治(2016年3期)2016-05-17 04:09:00
    選舉沸騰了黃土高原(上)
    公民與法治(2016年1期)2016-05-17 04:07:56
    利用基波相量變化率的快速選相方法
    川滇地區(qū)地殼應變能密度變化率與強震復發(fā)間隔的數(shù)值模擬
    地震地質(2015年3期)2015-12-25 03:29:53
    灑向黃土高原的愛
    中國火炬(2015年7期)2015-07-31 17:39:57
    噜噜噜噜噜久久久久久91| 国产精品美女特级片免费视频播放器| 狂野欧美白嫩少妇大欣赏| 精品亚洲乱码少妇综合久久| 欧美激情久久久久久爽电影| 91av网一区二区| 有码 亚洲区| 免费观看的影片在线观看| 国产乱人偷精品视频| 欧美激情久久久久久爽电影| 欧美bdsm另类| 国产成人a∨麻豆精品| 麻豆av噜噜一区二区三区| 亚洲三级黄色毛片| 精品国产一区二区三区久久久樱花 | 天堂中文最新版在线下载 | 国产 一区精品| 日韩欧美精品v在线| 亚洲激情五月婷婷啪啪| 别揉我奶头 嗯啊视频| 好男人在线观看高清免费视频| 亚洲精品久久午夜乱码| 国精品久久久久久国模美| 国产在视频线精品| 国产精品国产三级国产专区5o| 中文字幕久久专区| 你懂的网址亚洲精品在线观看| 日韩欧美三级三区| 亚洲精品日本国产第一区| 亚洲在线观看片| 欧美区成人在线视频| 亚洲欧美成人综合另类久久久| 两个人视频免费观看高清| 亚洲成人一二三区av| 久久久久久久久中文| 精品少妇黑人巨大在线播放| .国产精品久久| 亚洲精品久久久久久婷婷小说| 热99在线观看视频| 91久久精品国产一区二区三区| 插阴视频在线观看视频| 男人狂女人下面高潮的视频| 国产在视频线在精品| 日韩电影二区| 中文资源天堂在线| 亚洲精品视频女| 国产男人的电影天堂91| 春色校园在线视频观看| 久久精品熟女亚洲av麻豆精品 | 亚洲精品乱久久久久久| 亚洲va在线va天堂va国产| 色尼玛亚洲综合影院| 看十八女毛片水多多多| 99久久精品热视频| 大片免费播放器 马上看| 亚洲国产成人一精品久久久| 亚洲丝袜综合中文字幕| 熟妇人妻不卡中文字幕| 久久久久久国产a免费观看| or卡值多少钱| 69人妻影院| 日本一二三区视频观看| 韩国av在线不卡| 亚洲人成网站在线播| 成人欧美大片| 99久久九九国产精品国产免费| 永久免费av网站大全| 看非洲黑人一级黄片| 国产成年人精品一区二区| 青青草视频在线视频观看| 亚洲国产最新在线播放| 亚洲一区高清亚洲精品| 天堂网av新在线| 麻豆久久精品国产亚洲av| 九色成人免费人妻av| 国产精品久久久久久精品电影小说 | 男女国产视频网站| 国产精品蜜桃在线观看| 精品国产一区二区三区久久久樱花 | a级一级毛片免费在线观看| 亚洲精品aⅴ在线观看| 九色成人免费人妻av| 国产麻豆成人av免费视频| 免费看不卡的av| 亚洲精品久久久久久婷婷小说| 欧美一区二区亚洲| 亚洲人成网站在线播| 夫妻性生交免费视频一级片| 日本猛色少妇xxxxx猛交久久| 久99久视频精品免费| 少妇的逼水好多| 久久热精品热| 免费观看性生交大片5| 精品99又大又爽又粗少妇毛片| 国产精品.久久久| 日日撸夜夜添| 麻豆成人午夜福利视频| 久久精品国产亚洲av涩爱| 亚洲国产精品sss在线观看| 亚洲精品久久午夜乱码| 99久国产av精品| 网址你懂的国产日韩在线| 国产精品久久久久久精品电影| 免费观看的影片在线观看| 免费人成在线观看视频色| 午夜日本视频在线| 一个人观看的视频www高清免费观看| 国产久久久一区二区三区| 久久久久久久久久人人人人人人| 国产成年人精品一区二区| 亚洲欧美中文字幕日韩二区| 日韩,欧美,国产一区二区三区| 精品久久久久久久久亚洲| 嫩草影院入口| 欧美潮喷喷水| 日韩中字成人| 午夜福利视频1000在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲伊人久久精品综合| 精品人妻熟女av久视频| 国内精品一区二区在线观看| 最近最新中文字幕免费大全7| 亚洲国产色片| 人妻制服诱惑在线中文字幕| 国产亚洲91精品色在线| 国产精品国产三级国产专区5o| 久久鲁丝午夜福利片| 久久这里只有精品中国| 日本欧美国产在线视频| 午夜视频国产福利| 亚洲怡红院男人天堂| 春色校园在线视频观看| 国产成人一区二区在线| 久久99热这里只有精品18| 午夜免费激情av| 99久久精品一区二区三区| 亚洲精品亚洲一区二区| 亚洲国产精品专区欧美| 国产一区二区在线观看日韩| 狂野欧美激情性xxxx在线观看| 国产黄片视频在线免费观看| 久久99热这里只频精品6学生| 久久久久久九九精品二区国产| 亚洲av国产av综合av卡| 婷婷六月久久综合丁香| 国产高清国产精品国产三级 | 成人毛片60女人毛片免费| 久久99热这里只有精品18| 啦啦啦中文免费视频观看日本| 超碰97精品在线观看| 丝瓜视频免费看黄片| 久久精品国产亚洲av涩爱| 一个人看的www免费观看视频| 国内少妇人妻偷人精品xxx网站| 中国国产av一级| 成人亚洲精品av一区二区| av.在线天堂| av一本久久久久| 色视频www国产| 国产黄片美女视频| 国产精品精品国产色婷婷| 国产精品久久久久久久久免| 欧美精品一区二区大全| 日日摸夜夜添夜夜爱| 最近最新中文字幕大全电影3| 天天一区二区日本电影三级| 男女边摸边吃奶| 天堂√8在线中文| 欧美丝袜亚洲另类| .国产精品久久| 精品久久久久久久久av| 午夜精品国产一区二区电影 | 欧美日韩在线观看h| 亚洲高清免费不卡视频| 欧美高清成人免费视频www| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 国内精品一区二区在线观看| 2018国产大陆天天弄谢| 女的被弄到高潮叫床怎么办| 亚洲精品乱码久久久久久按摩| 欧美zozozo另类| 精品久久久噜噜| 神马国产精品三级电影在线观看| 又粗又硬又长又爽又黄的视频| 国国产精品蜜臀av免费| 22中文网久久字幕| av播播在线观看一区| 国国产精品蜜臀av免费| 国产精品久久久久久精品电影小说 | 日本-黄色视频高清免费观看| 久久久久久久久大av| 97热精品久久久久久| 熟妇人妻不卡中文字幕| 伦理电影大哥的女人| 久久精品久久久久久噜噜老黄| 日日干狠狠操夜夜爽| 极品教师在线视频| 特级一级黄色大片| 高清午夜精品一区二区三区| 国产成人一区二区在线| 边亲边吃奶的免费视频| 中国国产av一级| 性色avwww在线观看| 极品少妇高潮喷水抽搐| 亚洲色图av天堂| 国产成人91sexporn| 成年人午夜在线观看视频 | 我的老师免费观看完整版| 日韩强制内射视频| 午夜精品国产一区二区电影 | 网址你懂的国产日韩在线| av又黄又爽大尺度在线免费看| 国产亚洲精品久久久com| 欧美一级a爱片免费观看看| .国产精品久久| 精华霜和精华液先用哪个| 午夜福利视频精品| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 五月天丁香电影| 18+在线观看网站| 日韩精品有码人妻一区| 亚洲高清免费不卡视频| 99热这里只有是精品50| 黄色日韩在线| 最近最新中文字幕免费大全7| 国产麻豆成人av免费视频| 尾随美女入室| 日韩欧美精品v在线| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 天堂av国产一区二区熟女人妻| 亚洲精品自拍成人| 国产亚洲最大av| a级毛色黄片| 欧美人与善性xxx| 3wmmmm亚洲av在线观看| 91精品伊人久久大香线蕉| 老师上课跳d突然被开到最大视频| 在线天堂最新版资源| 国产精品人妻久久久影院| 中文乱码字字幕精品一区二区三区 | 精品国产一区二区三区久久久樱花 | 亚洲美女搞黄在线观看| av网站免费在线观看视频 | 日韩国内少妇激情av| 国产亚洲91精品色在线| 欧美成人午夜免费资源| 婷婷色麻豆天堂久久| 毛片女人毛片| 人人妻人人澡人人爽人人夜夜 | 美女被艹到高潮喷水动态| 精品国产三级普通话版| 欧美成人精品欧美一级黄| 天天躁夜夜躁狠狠久久av| 丰满少妇做爰视频| 中文资源天堂在线| 91久久精品国产一区二区三区| 在线观看人妻少妇| 亚洲欧美成人综合另类久久久| 99热这里只有是精品在线观看| 日本黄大片高清| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 99久久中文字幕三级久久日本| 欧美一区二区亚洲| 热99在线观看视频| 国产黄色小视频在线观看| 韩国av在线不卡| 人妻制服诱惑在线中文字幕| 少妇裸体淫交视频免费看高清| 国产精品麻豆人妻色哟哟久久 | 一级av片app| 国产永久视频网站| 老司机影院成人| 国产白丝娇喘喷水9色精品| 嫩草影院精品99| 成人综合一区亚洲| 91aial.com中文字幕在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人freesex在线| 91久久精品国产一区二区三区| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄| 国产精品国产三级国产专区5o| 1000部很黄的大片| 午夜福利视频1000在线观看| 国产有黄有色有爽视频| 一区二区三区四区激情视频| 日韩,欧美,国产一区二区三区| 嫩草影院精品99| 成年女人在线观看亚洲视频 | 午夜精品在线福利| 精品熟女少妇av免费看| 一本久久精品| 精品国产露脸久久av麻豆 | 中文资源天堂在线| 有码 亚洲区| 国产色爽女视频免费观看| 久久这里只有精品中国| 九九爱精品视频在线观看| 久久国产乱子免费精品| 国产一级毛片在线| 久久国产乱子免费精品| 久久草成人影院| 只有这里有精品99| 人妻系列 视频| 国产永久视频网站| 51国产日韩欧美| 人妻夜夜爽99麻豆av| 最近手机中文字幕大全| 在线天堂最新版资源| 欧美成人a在线观看| 免费大片黄手机在线观看| 国产精品女同一区二区软件| 看黄色毛片网站| 精品久久国产蜜桃| 久久久久久久亚洲中文字幕| 亚洲国产欧美在线一区| 欧美一级a爱片免费观看看| 一级毛片 在线播放| 国产v大片淫在线免费观看| 亚洲色图av天堂| 日韩欧美 国产精品| 亚洲av日韩在线播放| 淫秽高清视频在线观看| 99热这里只有是精品在线观看| 久久久久免费精品人妻一区二区| 日韩精品有码人妻一区| 老司机影院成人| 99re6热这里在线精品视频| 亚洲国产av新网站| 亚洲精品成人av观看孕妇| 网址你懂的国产日韩在线| 麻豆国产97在线/欧美| 亚洲精品中文字幕在线视频 | av福利片在线观看| 久久久久久国产a免费观看| 欧美极品一区二区三区四区| 国产成人精品久久久久久| 在线观看美女被高潮喷水网站| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 免费在线观看成人毛片| 蜜臀久久99精品久久宅男| 久久精品夜色国产| av在线观看视频网站免费| 在线免费观看不下载黄p国产| 日韩一本色道免费dvd| 国产淫片久久久久久久久| 亚洲欧美一区二区三区黑人 | 免费观看a级毛片全部| 色视频www国产| 亚洲精品,欧美精品| 免费在线观看成人毛片| 日韩欧美精品v在线| 啦啦啦啦在线视频资源| 精品人妻视频免费看| 婷婷色av中文字幕| 国产一区二区在线观看日韩| 亚洲图色成人| 国产 亚洲一区二区三区 | 九九爱精品视频在线观看| 男女下面进入的视频免费午夜| 久久久久久久大尺度免费视频| 亚洲欧美日韩无卡精品| 国产精品蜜桃在线观看| 国产伦精品一区二区三区视频9| 亚洲精品国产av蜜桃| 在线观看一区二区三区| 久久精品国产鲁丝片午夜精品| www.色视频.com| av线在线观看网站| 色吧在线观看| 亚洲三级黄色毛片| 色综合站精品国产| 久久久久免费精品人妻一区二区| 久久99热这里只频精品6学生| 成人高潮视频无遮挡免费网站| 久久6这里有精品| 黄色配什么色好看| 成人无遮挡网站| 最近中文字幕2019免费版| 成人无遮挡网站| 日本免费在线观看一区| 成年人午夜在线观看视频 | 亚洲欧美清纯卡通| 永久免费av网站大全| 成人性生交大片免费视频hd| 天天躁夜夜躁狠狠久久av| 在线免费观看不下载黄p国产| 国产精品国产三级专区第一集| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区| 婷婷色av中文字幕| 草草在线视频免费看| 国产在视频线在精品| 国产欧美另类精品又又久久亚洲欧美| 久久这里有精品视频免费| 91久久精品国产一区二区三区| 日韩成人伦理影院| 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| 欧美日韩国产mv在线观看视频 | 亚洲av二区三区四区| 能在线免费看毛片的网站| 亚洲第一区二区三区不卡| 在线 av 中文字幕| 免费高清在线观看视频在线观看| 如何舔出高潮| 日日撸夜夜添| 国产色婷婷99| 噜噜噜噜噜久久久久久91| 欧美人与善性xxx| 欧美zozozo另类| 大陆偷拍与自拍| 午夜福利视频1000在线观看| 中文字幕av成人在线电影| 久久99热6这里只有精品| 亚洲成人中文字幕在线播放| 精品人妻一区二区三区麻豆| 国产国拍精品亚洲av在线观看| 少妇的逼水好多| 在线观看美女被高潮喷水网站| 伦理电影大哥的女人| 午夜老司机福利剧场| 国产精品一区二区性色av| 水蜜桃什么品种好| 欧美日韩在线观看h| 成年av动漫网址| 美女黄网站色视频| 小蜜桃在线观看免费完整版高清| 少妇裸体淫交视频免费看高清| 欧美不卡视频在线免费观看| 国产亚洲精品久久久com| 日韩人妻高清精品专区| 亚洲一区高清亚洲精品| 菩萨蛮人人尽说江南好唐韦庄| 亚洲人成网站高清观看| 国产欧美另类精品又又久久亚洲欧美| 亚洲美女搞黄在线观看| 日本-黄色视频高清免费观看| 午夜福利在线观看免费完整高清在| 青春草亚洲视频在线观看| 99热6这里只有精品| 亚洲精品,欧美精品| 蜜桃亚洲精品一区二区三区| 国产精品久久久久久精品电影| 成人综合一区亚洲| 深爱激情五月婷婷| 成年版毛片免费区| 99久国产av精品国产电影| ponron亚洲| 成人鲁丝片一二三区免费| 久久久久精品久久久久真实原创| 国产成人a区在线观看| 色综合站精品国产| 黄色一级大片看看| 大香蕉97超碰在线| 国产黄色免费在线视频| 成年av动漫网址| 国产一区亚洲一区在线观看| 久久亚洲国产成人精品v| 91久久精品国产一区二区成人| 高清午夜精品一区二区三区| 国产中年淑女户外野战色| 欧美日韩精品成人综合77777| 午夜久久久久精精品| 日本猛色少妇xxxxx猛交久久| 麻豆成人av视频| 18禁在线无遮挡免费观看视频| av在线蜜桃| 国产亚洲一区二区精品| 国产精品久久久久久精品电影| 一级毛片aaaaaa免费看小| 欧美性感艳星| 黄色一级大片看看| 好男人在线观看高清免费视频| 亚洲真实伦在线观看| 91精品国产九色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产在线一区二区三区精| 又粗又硬又长又爽又黄的视频| 色5月婷婷丁香| 亚洲精品影视一区二区三区av| 1000部很黄的大片| 啦啦啦啦在线视频资源| 亚洲欧美精品专区久久| 久久99热6这里只有精品| 只有这里有精品99| 日韩制服骚丝袜av| 少妇人妻精品综合一区二区| 纵有疾风起免费观看全集完整版 | 国产亚洲av嫩草精品影院| 国产精品一区二区在线观看99 | 十八禁国产超污无遮挡网站| 国产伦精品一区二区三区视频9| 成年av动漫网址| 免费播放大片免费观看视频在线观看| 日本黄色片子视频| 日韩av不卡免费在线播放| 能在线免费观看的黄片| 一级毛片黄色毛片免费观看视频| 一级毛片我不卡| 欧美区成人在线视频| av在线播放精品| 大香蕉久久网| 在线观看免费高清a一片| 国内精品一区二区在线观看| 免费观看性生交大片5| 天天一区二区日本电影三级| 国产精品女同一区二区软件| 免费av不卡在线播放| 精品人妻熟女av久视频| 久久久久精品久久久久真实原创| 美女被艹到高潮喷水动态| 国产av国产精品国产| 麻豆久久精品国产亚洲av| 久久草成人影院| xxx大片免费视频| 搞女人的毛片| 亚洲av免费高清在线观看| 深夜a级毛片| 国产日韩欧美在线精品| 久久99热这里只频精品6学生| 亚洲欧美精品自产自拍| 99久国产av精品| 一级片'在线观看视频| 亚洲最大成人中文| 国产精品熟女久久久久浪| 中文欧美无线码| 真实男女啪啪啪动态图| 午夜亚洲福利在线播放| 麻豆久久精品国产亚洲av| 国产精品蜜桃在线观看| 床上黄色一级片| 白带黄色成豆腐渣| 最近2019中文字幕mv第一页| 久久99热这里只有精品18| 日韩欧美一区视频在线观看 | 一区二区三区高清视频在线| 国产91av在线免费观看| 久久久精品94久久精品| 色网站视频免费| 三级国产精品片| 亚洲一级一片aⅴ在线观看| 欧美+日韩+精品| 日韩亚洲欧美综合| 美女主播在线视频| 欧美潮喷喷水| 国产精品久久视频播放| 免费高清在线观看视频在线观看| av在线老鸭窝| 啦啦啦中文免费视频观看日本| 91久久精品电影网| 岛国毛片在线播放| 日韩 亚洲 欧美在线| 国产一区二区三区av在线| 黄色欧美视频在线观看| 一二三四中文在线观看免费高清| 天天躁夜夜躁狠狠久久av| 直男gayav资源| 欧美性猛交╳xxx乱大交人| 精品国产露脸久久av麻豆 | 亚洲精品乱码久久久久久按摩| 久久精品国产亚洲av天美| 乱人视频在线观看| 永久网站在线| 99热这里只有精品一区| 成人特级av手机在线观看| 成人亚洲欧美一区二区av| 啦啦啦啦在线视频资源| 2018国产大陆天天弄谢| 人人妻人人澡人人爽人人夜夜 | 国产爱豆传媒在线观看| 欧美高清性xxxxhd video| 亚洲最大成人av| 夜夜看夜夜爽夜夜摸| 欧美成人精品欧美一级黄| 欧美另类一区| 精品国产露脸久久av麻豆 | 亚洲精品久久久久久婷婷小说| 免费黄频网站在线观看国产| 日本免费a在线| 国产女主播在线喷水免费视频网站 | 国产爱豆传媒在线观看| 又爽又黄a免费视频| 午夜亚洲福利在线播放| 午夜福利高清视频| 直男gayav资源| 免费大片18禁| 欧美精品一区二区大全| 国产乱人偷精品视频| 日韩强制内射视频| 国产精品.久久久| 免费观看的影片在线观看| 午夜福利在线在线| 亚洲精品成人av观看孕妇| 亚洲欧美成人精品一区二区| 国产黄片视频在线免费观看| 天堂中文最新版在线下载 | 国产精品精品国产色婷婷| 大陆偷拍与自拍| 免费在线观看成人毛片| 国内揄拍国产精品人妻在线| 一级毛片aaaaaa免费看小| 天美传媒精品一区二区| 欧美区成人在线视频| 狂野欧美白嫩少妇大欣赏| 日韩中字成人| 亚洲在线观看片| 一级二级三级毛片免费看| 在线a可以看的网站|