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

    京津風(fēng)沙源治理區(qū)植被變化的可持續(xù)性分析

    2016-11-01 06:52:03楊艷麗孫艷玲王中良
    關(guān)鍵詞:趨勢區(qū)域研究

    楊艷麗,孫艷玲,王中良,康 杰,

    (天津師范大學(xué)a.城市與環(huán)境科學(xué)學(xué)院,b.天津市水資源與水環(huán)境重點實驗室,天津 300387)

    京津風(fēng)沙源治理區(qū)植被變化的可持續(xù)性分析

    楊艷麗a,孫艷玲a,王中良b,康杰a,b

    (天津師范大學(xué)a.城市與環(huán)境科學(xué)學(xué)院,b.天津市水資源與水環(huán)境重點實驗室,天津 300387)

    利用京津風(fēng)沙源治理區(qū)1982—2011年的GIMMS/NDVI3g數(shù)據(jù)集,對治理區(qū)植被變化的可持續(xù)性進(jìn)行定量分析.結(jié)果表明:(1)歸一化植被指數(shù)(NDVI)時間變化方面,2001—2011年間研究區(qū)植被整體上呈增加趨勢,4大分區(qū)植被增速存在區(qū)域差異;(2)NDVI變化的可持續(xù)性方面,2001—2011年間整個研究區(qū)植被變化呈現(xiàn)出顯著的長期持續(xù)性特征,且存在明顯的時空差異,自2000年治理工程實施以來,治理區(qū)植被明顯恢復(fù),一期工程成效顯著,但局部地區(qū)植被仍在持續(xù)減少,因此有必要進(jìn)行二期工程以鞏固和提高一期工程的成果,進(jìn)一步改善環(huán)境;(3)NDVI變化趨勢的未來持續(xù)時間方面,整個治理區(qū)植被的未來變化狀況不容樂觀,4大分區(qū)中農(nóng)牧交錯帶沙化土地治理區(qū)植被的未來發(fā)展趨勢優(yōu)于其他3個治理區(qū),工程治理效果也相對較好.

    京津風(fēng)沙源治理區(qū);歸一化植被指數(shù);植被恢復(fù);變化趨勢;可持續(xù)性

    植被作為生態(tài)系統(tǒng)的重要組成部分,是連接自然界中各要素的紐帶和表征區(qū)域環(huán)境演變的敏感指示器[1],在調(diào)節(jié)區(qū)域氣候、維持地表能量平衡方面扮演著重要角色.植被變化作為生態(tài)環(huán)境變化的直接結(jié)果,能夠反映地表生境變化和區(qū)域環(huán)境總體狀況,在全球變化研究領(lǐng)域受到高度關(guān)注[2].遙感數(shù)據(jù)為植被變化研究提供了重要的數(shù)據(jù)來源[3].歸一化植被指數(shù)(normalized difference vegetation index,NDVI)具有良好的植被信息表達(dá)能力和數(shù)據(jù)提取過程中較強的抗干擾能力,在植被變化研究中應(yīng)用廣泛[4].目前常用的NDVI數(shù)據(jù)包括SPOT VGT/NDVI、MODIS/NDVI和GIMMS/NDVI等[5].其中,GIMMS/NDVI數(shù)據(jù)時間序列最長,覆蓋范圍廣,因此在長時期、大尺度的植被宏觀動態(tài)監(jiān)測中應(yīng)用較廣[5-6].與傳統(tǒng)的GIMMS/NDVI數(shù)據(jù)相比,GIMMS/NDVI3g數(shù)據(jù)精度更高,可以更加準(zhǔn)確地反映區(qū)域植被變化特征及生態(tài)環(huán)境狀況[7-9].目前很多研究都基于NDVI序列數(shù)據(jù),采用線性趨勢法、Sen+Mann-Kendall方法和Hurst指數(shù)對不同時空尺度上的植被變化趨勢進(jìn)行了探討分析[1-2].

    京津風(fēng)沙源治理區(qū)整體上屬于干旱與半干旱區(qū),大部分處在我國的農(nóng)牧交錯帶上,生態(tài)環(huán)境比較脆弱.京津風(fēng)沙源治理工程自2000年試點以來,已受到國家和社會各界的廣泛關(guān)注.近年來,不斷有學(xué)者對治理區(qū)的植被變化狀況進(jìn)行調(diào)查與評估研究,如嚴(yán)恩萍等[10]和裴亮等[11]分別利用2000—2012年和2000—2008年的MODIS/NDVI時序數(shù)據(jù)對京津風(fēng)沙源治理區(qū)植被覆蓋的時空變化特征及影響因子進(jìn)行了分析,得出研究區(qū)植被覆蓋總體呈上升趨勢、植被變化整體呈中強持續(xù)性特征的結(jié)論.但這些研究主要是對植被變化狀況及未來發(fā)展趨勢的定性描述,難以準(zhǔn)確反映出植被在2000年以后的可持續(xù)性變化特征,也無法定量預(yù)測出植被變化趨勢在未來持續(xù)的時間長度.由Lanfredi等[12]提出的計算植被變化可持續(xù)性的方法,不僅可以計算出植被在研究時間段內(nèi)的可持續(xù)性狀況,還能定量預(yù)測出植被變化趨勢未來持續(xù)的時間長度,目前已被很多學(xué)者用來分析植被變化的可持續(xù)性特征[13-15].因此,本文利用京津風(fēng)沙源治理區(qū)1982—2011年的GIMMS/NDVI3g數(shù)據(jù),采用Lanfredi等[12]的方法定量分析該地區(qū)植被變化的可持續(xù)性狀況,揭示研究區(qū)植被的可持續(xù)性變化特征和未來演變趨勢,評價工程區(qū)生態(tài)治理效果,為區(qū)域植被恢復(fù)工作與生態(tài)建設(shè)提供一定的決策支持.

    1 資料來源和研究方法

    1.1京津風(fēng)沙源治理區(qū)概況

    研究區(qū)(109°30′~119°20′E、38°50′~46°40′N)西起內(nèi)蒙古達(dá)茂旗,東至河北省平泉縣,南起山西省代縣,北至內(nèi)蒙古東烏珠穆沁旗,涉及北京、天津、河北、山西及內(nèi)蒙古5?。ㄗ灾螀^(qū)、直轄市)的75個縣(旗、市、區(qū)),總面積45.8萬km2,如圖1所示.

    圖1 京津風(fēng)沙源治理區(qū)位置示意圖Fig.1 Location of Beijing-Tianjin Sandstorm Source Control Region of China

    研究區(qū)由平原、山地、高原3大地貌類型組成.區(qū)內(nèi)氣候復(fù)雜,區(qū)域差異性明顯,氣候干旱、熱量偏低、多風(fēng)是本區(qū)氣候的典型特征.近30 a來,本區(qū)氣候呈現(xiàn)暖干化趨勢,干旱區(qū)和干旱亞濕潤區(qū)的界限逐步向東推進(jìn)[16].京津風(fēng)沙源治理工程于2000年試點,2002年全面實施,到2010年為止,工程區(qū)共完成退耕還林及造林6.0萬km2,草地治理8.67萬km2,生態(tài)環(huán)境得到顯著改善.按照國家規(guī)劃,工程治理區(qū)分為北部干旱草原治理區(qū)、農(nóng)牧交錯帶沙化土地治理區(qū)、燕山丘陵山地水源保護(hù)區(qū)和渾善達(dá)克沙地治理區(qū).

    1.2數(shù)據(jù)來源與處理

    采用京津風(fēng)沙源治理區(qū)1982—2011年的GIMMS/NDVI3g數(shù)據(jù),這是由美國國家航天航空局(NASA)全球監(jiān)測與模型研究組(Global Inventory Modeling and MappingStudies,GIMMS)發(fā)布的第3代GIMMS時序NDVI數(shù)據(jù),空間分辨率為0.082 5°,時間分辨率為15d.該數(shù)據(jù)集是經(jīng)過天頂角、氣溶膠、云層覆蓋等因素的影響校正后發(fā)布的相對標(biāo)準(zhǔn)的數(shù)據(jù),較其他的NDVI數(shù)據(jù)精度更高、誤差更小,適用于植被覆蓋變化的長期監(jiān)測[7-9].進(jìn)行數(shù)據(jù)處理時,首先對GIMMS/NDVI3g數(shù)據(jù)源進(jìn)行數(shù)據(jù)格式轉(zhuǎn)換,然后借助ArcGIS10.2軟件利用研究區(qū)邊界對其進(jìn)行裁剪,再采用最大值合成法(MVC)分別對每年每半月的NDVI值進(jìn)行最大化處理,獲取每年的NDVI最大值,代表當(dāng)年的NDVI值.

    1.3研究方法

    1.3.1趨勢分析

    利用一元線性回歸方法,逐像元分別計算1982—2000年、1982—2001年、…、1982—2011年間的NDVI變化趨勢,即Slope值,來反映研究區(qū)各時段內(nèi)的植被變化趨勢[4].其計算公式為:

    式中:Slope為像元NDVI的斜率值,若Slope>0,表示區(qū)域植被呈增加趨勢,且值愈大表明植被增加趨勢愈明顯,反之則呈減少趨勢;若值為0,表示植被沒有變化;n為研究時段的長度;i代表年份;為平均年份;NDVIi為第i年的NDVI值;為NDVI多年平均值.

    1.3.2可持續(xù)性分析

    Lanfredi等[12]的方法計算的前提是要選取一個參考時間段,也稱為初始時間段.首先,根據(jù)研究區(qū)情況,選取1982—2000年作為研究的初始時間段,用來分析2001—2011年間研究區(qū)植被變化的可持續(xù)性特征;之后采用一元線性回歸方法逐像元分別計算1982—2000、1982—2001、1982—2002、1982—2003、1982—2004、…、1982—2011年的NDVI變化趨勢,用s(x,y,t)表示,其中(x,y)代表研究區(qū)內(nèi)每個像元的坐標(biāo)位置,t代表時間節(jié)點,分別為2000、2001、…、2011.如果植被呈增加趨勢,則s(x,y,t)=1,否則s(x,y,t)= -1.在此基礎(chǔ)上,逐像元依次計算1982—2000、1982—2001、1982—2002、1982—2003、1982—2004、…、1982—2011年的植被變化的可持續(xù)性.用p(x,y,t)表示,其中(x,y)代表研究區(qū)內(nèi)每個像元的坐標(biāo)位置,t代表時間節(jié)點,分別為2000、2001、…、2011.其計算公式為:

    其中,1982—2000年間的p(x,y)值與該時間段的s(x,y)值相等,其他時間段的p(x,y)值根據(jù)公式(2)確定.以1982—2001年為例,假定某像元上1982—2000年間s(x,y)=1,則1982—2000年間的p(x,y)=1.若1982—2001年間s(x,y)=1,那么1982—2001年間的p(x,y)=1;若1982—2001年s(x,y)=-1,那么該時段的p(x,y)=0.以此類推.

    將1982—2001、1982—2002、1982—2003、…、1982—2011年11個時間段的p(x,y)相加,得到研究區(qū)NDVI變化的可持續(xù)性累加圖,以此來分析整個研究區(qū)在2001—2011年間植被變化的可持續(xù)性狀況.其計算公式為:

    接著分別計算研究區(qū)各時段內(nèi)NDVI增加(減少)趨勢的持續(xù)性概率q(t),即1982—2001、1982—2002、1982—2003、…、1982—2011年的NDVI增加(減少)趨勢的柵格數(shù)與1982—2000年NDVI增加(減少)趨勢的柵格數(shù)之比.

    再分別計算各時段內(nèi)NDVI增加(減少)趨勢未來持續(xù)的時間,計算公式為:

    式中:τ代表NDVI增加(減少)趨勢未來持續(xù)的時間長度;1982—2001、1982—2002、1982—2003、…、1982—2011對應(yīng)的t值分別為1、2、3、…、11.τ值越大,說明這種趨勢持續(xù)時間愈長,可持續(xù)性愈強.

    最后,將各時段內(nèi)NDVI增加(減少)趨勢的τ值相加求平均值,得到NDVI增加(減少)趨勢未來持續(xù)的平均時間.

    2 結(jié)果與分析

    2.1NDVI時間變化特征

    通過分別統(tǒng)計京津風(fēng)沙源治理區(qū)及4大分區(qū)2001—2011年的逐年NDVI值,制作NDVI年際變化圖,分析其植被的時間變化趨勢,結(jié)果如圖2所示.

    圖2 2001—2011年研究區(qū)NDVI年際變化Fig.2 Inter-annual variation of NDVI for study areas from 2001 to 2011

    由圖2可以看出,2001—2011年間研究區(qū)植被整體上呈增加趨勢,平均增速為每10 a增加0.02.4大分區(qū)在2001—2011年間植被也都呈增加趨勢,但增速存在區(qū)域差異.其中,燕山丘陵山地水源保護(hù)區(qū)的植被增加速度最快,平均增速為每10 a增加0.05;農(nóng)牧交錯帶沙化土地治理區(qū)和北部干旱草原治理區(qū)次之,平均增速分別為每10 a增加0.02和0.01;渾善達(dá)克沙地治理區(qū)的植被平均增速最小,每10 a增加0.003.

    2.2NDVI變化的可持續(xù)性特征

    利用京津風(fēng)沙源治理區(qū) 1982—2011年的GIMMS/NDVI3g數(shù)據(jù)集,采用可持續(xù)性分析方法將1982—2001、1982—2002、1982—2003、…、1982—2011年的持續(xù)性相加,計算得到NDVI變化趨勢的持續(xù)性疊加圖,以此來反映研究區(qū)2001—2011年間植被變化趨勢持續(xù)性的時空特征,結(jié)果如圖3所示.

    圖3 2001—2011年研究區(qū)NDVI變化趨勢的可持續(xù)性疊加圖Fig.3 Cumulative persistence map of NDVI change trends for study areas from 2001 to 2011

    從圖3可以看出,在整個研究區(qū)內(nèi),植被變化可持續(xù)性累加值為11的區(qū)域分布面積最大,占研究區(qū)總面積的34.91%,集中分布在農(nóng)牧交錯帶沙化土地治理區(qū)的幾乎全部區(qū)域、渾善達(dá)克沙地治理區(qū)和燕山丘陵山地水源保護(hù)區(qū)的南部、北部干旱草原治理區(qū)的西北部和西南部,這些區(qū)域在2001—2011年間植被呈增加趨勢,表明其植被恢復(fù)效果顯著.累加值為-11的區(qū)域分布面積次之,占12.17%,主要分布在北部干旱草原治理區(qū)東北部的東烏珠穆沁旗和西烏珠穆沁旗、燕山丘陵山地水源保護(hù)區(qū)的北部和西北部、渾善達(dá)克沙地治理區(qū)北部和中西部部分區(qū)域,這些區(qū)域在2000—2011年間植被呈減少趨勢,表明其植被在持續(xù)退化,盡管在這些區(qū)域?qū)嵤┝酥脖换謴?fù)治理措施,但仍未收到良好效果,有關(guān)部門應(yīng)予以足夠重視.累加值為0的區(qū)域占9.89%,零星分布在北部干旱草原治理區(qū)的東部和西部部分區(qū)域、渾善達(dá)克沙地治理區(qū)的北部和西部小部分地區(qū)、燕山丘陵山地水源保護(hù)區(qū)北部的局部地區(qū),這些區(qū)域在2001—2011年間植被變化趨勢的持續(xù)性無法確定.在所有正值中,累加值為7的區(qū)域所占面積比例最小,為1.19%,呈散點狀分布在燕山丘陵山地水源保護(hù)區(qū)、渾善達(dá)克沙地治理區(qū)的小部分區(qū)域,這些區(qū)域在2001—2007年間植被呈增加趨勢,從2008年開始這一趨勢發(fā)生改變,說明2008年在該區(qū)域?qū)嵤┑纳鷳B(tài)治理工程效果欠佳,也可能是由于這一年該區(qū)域的一些不合理人類活動干擾了植被生長.在所有負(fù)值中,值為-6的區(qū)域分布面積最少,占0.06%,只出現(xiàn)在北部干旱草原治理區(qū)的包頭市市轄區(qū)內(nèi)和燕山丘陵山地水源保護(hù)區(qū)的豐寧滿族自治縣極小部分區(qū)域,這些區(qū)域植被減少趨勢一直持續(xù)了6 a,直到2007年這一趨勢才發(fā)生改變,表明2007年該地區(qū)植被生長狀況發(fā)生好轉(zhuǎn),生態(tài)治理工作效益凸顯.應(yīng)該指出的是,增加(減少)趨勢持續(xù)性并不表示在2001—2011年間植被一定是逐年增加(減少)的,而是說明植被變化的可持續(xù)性較好,仍保持著最初的增加(減少)趨勢.

    綜上可見,京津風(fēng)沙源治理區(qū)在2001—2011年間植被增加趨勢持續(xù)了11 a的區(qū)域分布最廣,表明自2000年治理工程實施以來,研究區(qū)植被明顯恢復(fù),生態(tài)環(huán)境得到改善.但是減少趨勢持續(xù)了11 a的區(qū)域占研究區(qū)總面積的12.17%,說明雖然一期治理工程已取得顯著成效,但是治理區(qū)仍存在一些問題,局部地區(qū)植被仍在持續(xù)減少.因此,有必要在二期工程實施中加強植被恢復(fù)工作來鞏固和提升一期工程建設(shè)效果,進(jìn)一步改善環(huán)境.此外,值為11和-11的區(qū)域占整個研究區(qū)的47.09%,表明整個研究區(qū)植被變化呈現(xiàn)出顯著的長期持續(xù)性特征且存在明顯的時空差異,因此生態(tài)治理工作應(yīng)結(jié)合治理區(qū)內(nèi)部差異,分階段、分地域進(jìn)行,積極關(guān)注植被持續(xù)減少區(qū)域,及時采取措施,提高研究區(qū)整體生態(tài)效益.

    2.3NDVI變化趨勢的未來持續(xù)時間特征

    2.3.1整個研究區(qū)NDVI變化趨勢的未來持續(xù)時間

    植被變化趨勢的持續(xù)性概率圖可以反映植被增加(減少)趨勢持續(xù)性的概率和時間之間的對應(yīng)關(guān)系,以及植被變化趨勢未來持續(xù)的時間長度.本研究分別計算了研究區(qū)各時段內(nèi)NDVI增加(減少)趨勢持續(xù)的時間,將其相加求平均值得到整個研究區(qū)NDVI增加(減少)趨勢未來持續(xù)的平均時間,以此來評估植被變化趨勢的未來持續(xù)時間特征,結(jié)果如圖4所示.由圖4可以看出,整個研究區(qū)植被增加趨勢持續(xù)性概率在每個時間段都小于減少趨勢持續(xù)性概率,而且下降幅度也大于減少趨勢持續(xù)性概率.對于生長狀態(tài)良好的植被,其增加趨勢要比減少趨勢持續(xù)時間更長.整個研究區(qū)植被減少趨勢未來持續(xù)的平均時間τ-為29 a,增加趨勢未來持續(xù)的平均時間τ+為12 a,增加趨勢比減少趨勢持續(xù)時間更短,說明從整體上來看研究區(qū)未來植被變化狀況不容樂觀,系統(tǒng)的穩(wěn)定性仍較差.因此,京津風(fēng)沙源區(qū)治理工作還有待進(jìn)一步加強實施,以延長植被增加趨勢的持續(xù)時間,增強系統(tǒng)的持續(xù)性和穩(wěn)定性,保證區(qū)域生態(tài)安全.

    圖4 研究區(qū)植被變化趨勢的持續(xù)性概率Fig.4 Persistence probability of vegetation change trends for the whole study area

    2.3.24大治理區(qū)NDVI變化趨勢的未來持續(xù)時間

    分別分析京津風(fēng)沙源治理區(qū)內(nèi)4大分區(qū)的NDVI變化趨勢未來持續(xù)時間,結(jié)果如圖5所示.

    從持續(xù)性概率來看,4大治理分區(qū)中,北部干旱草原治理區(qū)、燕山丘陵山地水源保護(hù)區(qū)、渾善達(dá)克沙地治理區(qū)這3個治理區(qū)的植被減少趨勢持續(xù)性概率都大于同時段內(nèi)的增加趨勢持續(xù)性概率;而農(nóng)牧交錯帶沙化土地治理區(qū)從2002年開始,增加趨勢的持續(xù)性概率都大于同時段的減少趨勢持續(xù)性概率,這與京津風(fēng)沙源一期治理工程全面實施時間正好吻合,說明該治理工程對農(nóng)牧交錯帶沙化土地治理區(qū)的植被恢復(fù)起到了立竿見影的效果.從植被變化趨勢未來持續(xù)時間來看,除農(nóng)牧交錯帶沙化土地治理區(qū)外,其余3個治理區(qū)的植被減少趨勢未來持續(xù)時間都大于增加趨勢未來持續(xù)時間.增加趨勢未來持續(xù)時間最長的是農(nóng)牧交錯帶沙化土地治理區(qū),達(dá)到36 a.而減少趨勢持續(xù)時間最長的是渾善達(dá)克沙地治理區(qū),達(dá)到40 a.說明4大治理區(qū)中,農(nóng)牧交錯帶沙化土地治理區(qū)的植被恢復(fù)效果最好,工程治理效果最優(yōu).而其余3個治理區(qū)可能受城市化、人類活動等因素的干擾,工程治理效果和植被恢復(fù)狀況相對較差,積極恢復(fù)治理這些區(qū)域的植被是維持整個京津風(fēng)沙源治理區(qū)系統(tǒng)持久性和穩(wěn)定性的關(guān)鍵.因此,應(yīng)從整體和系統(tǒng)的角度出發(fā),重點加強北部干旱草原治理區(qū)、燕山丘陵山地水源保護(hù)區(qū)、渾善達(dá)克沙地治理區(qū)這3個治理區(qū)植被的改善工作,以從整體上提高研究區(qū)的生態(tài)治理效果,促進(jìn)區(qū)域整體的可持續(xù)發(fā)展.

    圖5 4大治理區(qū)植被變化趨勢的持續(xù)性概率Fig.5 Persistence probability of vegetation change trends for four subregions

    3 討論與結(jié)論

    3.1討論

    本文基于京津風(fēng)沙源治理區(qū)1982—2011年的GIMMS/NDVI3g數(shù)據(jù)集,采用植被變化可持續(xù)性分析方法,定量分析了該區(qū)植被變化的可持續(xù)性特征及未來持續(xù)時間.京津風(fēng)沙源治理區(qū)2001—2011年植被呈增加趨勢,表明治理工程正向作用顯著,但研究區(qū)內(nèi)4大分區(qū)植被變化的可持續(xù)性存在差異.植被變化是氣候變化和人類活動綜合作用的結(jié)果[1].研究區(qū)氣候干旱,因此降水是影響其植被變化的主要氣候因子,尤其是以干旱為主的北部干旱草原治理區(qū).已有研究發(fā)現(xiàn),近30 a來,京津風(fēng)沙源治理區(qū)氣溫顯著上升,降水和濕潤指數(shù)總體下降,位于本區(qū)東部的部分區(qū)域已經(jīng)由亞濕潤干旱區(qū)變?yōu)楦珊祬^(qū)[16].再加上不合理的人類活動和土地利用方式的轉(zhuǎn)變,如渾善達(dá)克沙地治理區(qū)東部林地變草地和耕地、北部干旱草原治理區(qū)的過度放牧和工礦建設(shè)等導(dǎo)致的草地退化、燕山丘陵山地水源保護(hù)區(qū)草地開墾和居民地擴張侵占耕地[3]等.因此,北部干旱草原治理區(qū)、燕山丘陵山地水源保護(hù)區(qū)和渾善達(dá)克沙地治理區(qū)的植被增加趨勢的未來持續(xù)時間較短,工程治理效果不顯著;而農(nóng)牧交錯帶沙化土地治理區(qū)是京津風(fēng)沙源工程的重點治理區(qū)域[10],在營造農(nóng)田防護(hù)林、變放牧為圈養(yǎng)、退耕還林等工程建設(shè)的推動下,植被增加趨勢的未來持續(xù)時間較長,系統(tǒng)穩(wěn)定性較好.這些結(jié)果可為研究區(qū)制定生態(tài)保護(hù)政策提供有效參考.另一方面,由于該工程涉及不同省市,工程治理措施和側(cè)重點因地而異,植被變化也相應(yīng)地存在一定地域差異.因此應(yīng)加強小尺度的植被調(diào)查研究,為京津風(fēng)沙源治理工程提供實用的實施對策和建議,以保證該治理工程實現(xiàn)生態(tài)效益和社會效益的最大化,促進(jìn)區(qū)域整體的可持續(xù)發(fā)展.

    3.2結(jié)論

    (1)從NDVI時間變化特征來看,2001—2011年間京津風(fēng)沙源治理區(qū)植被整體上呈增加趨勢,平均增速為每10 a增加0.02.4大分區(qū)植被也呈增加趨勢,但增速存在區(qū)域差異.

    (2)從NDVI變化的可持續(xù)性特征來看,研究區(qū)植被變化呈現(xiàn)出顯著的長期持續(xù)性特征,且存在明顯的時空差異.2001—2011年間,植被增加趨勢持續(xù)了11 a的區(qū)域占總面積的34.91%,減少趨勢持續(xù)了11 a的區(qū)域占12.17%,表明自2000年京津風(fēng)沙源治理工程實施以來,治理區(qū)植被明顯恢復(fù),一期工程已取得顯著成效.但是治理區(qū)仍存在一些問題,局部地區(qū)植被仍在持續(xù)減少.

    (3)從NDVI變化趨勢的未來持續(xù)時間來看,整個治理區(qū)植被減少趨勢持續(xù)的平均時間為29 a,增加趨勢持續(xù)的平均時間為12 a,說明從整體上來看治理區(qū)植被未來變化狀況不容樂觀.4大分區(qū)中,除農(nóng)牧交錯帶沙化土地治理區(qū)外,其余3個治理區(qū)的減少趨勢持續(xù)時間都長于增加趨勢持續(xù)時間,說明農(nóng)牧交錯帶沙化土地治理區(qū)植被的未來發(fā)展趨勢優(yōu)于其他3個治理區(qū),工程治理效果也相對較好.

    [1]楊林山,李常斌,王帥兵,等.氣候變化和人類活動對洮河流域植被動態(tài)的影響研究[J].資源科學(xué),2014,36(9):1941-1948.YANG L S,LI C B,WANG S B,et al.Impacts of climate change and human activities on vegetation dynamics in two geo-ecologic regions in the Taohe river basin,northwest China[J].Resources Science,2014,36(9):1941-1948(in Chinese).

    [2]溫曉金,劉焱序,楊新軍.恢復(fù)力視角下生態(tài)型城市植被恢復(fù)空間分異及其影響因素[J].生態(tài)學(xué)報,2015,35(13):4377-4389.WEN X J,LIU Y X,YANG X J.A resilience-based analysis on the spatial heterogeneity of vegetation restoration and its affecting factors in the construction of eco-cities:a case study of Shangluo,Shanxi[J].Acta Ecologica Sinica,2015,35(13):4377-4389(in Chinese).

    [3]劉洛,徐新良,段建南,等.京津風(fēng)沙源區(qū)生態(tài)環(huán)境時空變化的遙感監(jiān)測分析[J].地球信息科學(xué)學(xué)報,2011,13(6):819-824.LIU L,XU X L,DUAN J N,et al.The spatial-temporal changes monitoring of ecological environment in source area of the sand and dust endangering Beijing-Tianjin by remote sensing[J].Journal of Geo-Information Science,2011,13(6):819-824(in Chinese).

    [4]呂洋,董國濤,楊勝天,等.雅魯藏布江流域NDVI時空變化及其與降水和高程的關(guān)系[J].資源科學(xué),2014,36(3):603-611.Lü Y,DONG G T,YANG S T,et al.Spatio-temporal variation in NDVI in the Yarlung Zangbo river basin and its relationship with precipitation and elevation[J].Resources Science,2014,36(3):603-611(in Chinese).

    [5]侯美亭,趙海燕,王箏,等.基于GIMMS、VGT和MODIS的中國東部植被指數(shù)對比分析[J].遙感技術(shù)與應(yīng)用,2013,28(2):290-299.HOU M T,ZHAO H Y,WANG Z,et al.Comparison of GIMMS,VGT and MODIS vegetation index datasets in eastern China[J].Remote Sensing Technology and Application,2013,28(2):290-299(in Chinese).

    [6]孫艷玲,郭鵬.1982—2006年華北植被指數(shù)時空變化特征[J].干旱區(qū)研究,2012,29(2):187-193.SUNY L,GUO P.Spatio-temporal variation of vegetation coverage index in North China during the period from 1982 to 2006[J].Arid Zone Research,2012,29(2):187-193(in Chinese).

    [7]XU G,ZHANG H F,CHEN B Z,et al.Changes in vegetation growth dynamics and relations with climate over China's landmass from 1982 to 2011[J].Remote Sensing,2014,6(4):3263-3283.

    [8]ZHU Z,BI J,PAN Y,et al.Global data sets of vegetation leaf area index(LAI)3g and fraction of photosynthetically active radiation(FPAR)3g derived from global inventory modeling and mapping studies(GIMMS)normalized difference vegetation index(NDVI3g)for the period 1981 to 2011[J].Remote Sens,2013,5(2):927-948.

    [9]范德芹,朱文泉,潘耀忠,等.青藏高原小嵩草高寒草甸返青期遙感識別方法篩選[J].遙感學(xué)報,2014,18(5):1117-1127.FAN D Q,ZHU W Q,PAN Y Z,et al.Identifying an optimal method for estimating green-up data of Kobresia pygmaea alpine meadow in Qinghai-Tibetan Plateau[J].Journal of Remote Sensing,2014,18(5):1117-1127(in Chinese).

    [10]嚴(yán)恩萍,林輝,黨永鋒,等.2000—2012年京津風(fēng)沙源治理區(qū)植被覆蓋時空演變特征[J].生態(tài)學(xué)報,2014,34(17):5007-5020.YAN E P,LIN H,DANG Y F,et al.The spatiotemporal changes of vegetation cover in Beijing-Tianjin sandstorm source control region during 2000—2012[J].Acta Ecologica Sinica,2014,34(17):5007-5020(in Chinese).

    [11]裴亮,黃森旺,陳麗萍.京津風(fēng)沙源區(qū)植被的時空變化及其對氣候因子的響應(yīng)[J].中國沙漠,2013,33(5):1593-1597.PEI L,HUANG S W,CHEN L P.Vegetation spatio-temporal changes and the relationship with climate factors in Beijing-Tianjin sand source region[J].Journal of Desert Research,2013,33(5):1593-1597(in Chinese).

    [12]LANFREDI M,SIMONIELLO T,MACCHIATO M.Temporal persistence in vegetation cover changes observed from satellite:development of an estimation procedure in the test site of the Mediterranean Italy[J].Remote Sens Environ,2004,93(4):565-576.

    [13]SIMONIELLO T,LANFREDI M,LIBERTI M,et al.Estimation of vegetation cover resilience from satellite time series[J].Hydrol Earth Syst Sci,2008,12(4):1053-1064.

    [14]COPPOLA R,CUOMO V,D'EMILIO M,et al.Terrestrial vegetation cover activity as a problem of fluctuating surfaces[J].Int J Mod Phys B,2009,23(28):5444-5452.

    [15]HARRIS A,CARR A S,DASH J.Remote sensing of vegetation cover dynamics and resilience across southern Africa[J].International Journal of Applied Earth Observation and Geoinformation,2014,28(5):131-139.

    [16]孫斌,高志海,王紅巖,等.近30年京津風(fēng)沙源區(qū)氣候干濕變化分析[J].干旱區(qū)資源與環(huán)境,2014,28(11):164-170.SUN B,GAO Z H,WANG H Y,et al.Dry/wet variation of Beijing-Tianjin dust and sandstorm source region during 1981—2010[J].Journal of Arid Land Resources and Environment,2014,28(11):164-170(in Chinese).

    (責(zé)任編校紀(jì)翠榮)

    Persistence in vegetation cover changes in Beijing-Tianjin Sandstorm Source Control Region

    YANG Yanlia,SUN Yanlinga,WANG Zhongliangb,KANG Jiea,b
    (a.College of Urban and Environment Science,b.Tianjin Key Laboratory of Water Resources and Environment,Tianjin Normal University,Tianjin 300387,China)

    Based on the GIMMS/NDVI3g data set of Beijing-Tianjin Sandstorm Source Control Region from 1982 to 2011,the persistence in vegetation cover changes in this region was quantitatively analyzed.The results showed that:(1)As for the temporal variation of normalized difference vegetation index(NDVI),the vegetation showed an increasing trend on the whole in Beijing-Tianjin Sandstorm Source Control Region from 2000 to 2011.And there were regional differences in the vegetation growth rate among four regions;(2)As for the persistence in NDVI changes,the vegetation changes in the whole study area showed significant long-time persistence characteristics,and there were obvious differences in time and space for the persistence from 2000 to 2011.Since the Beijing-Tianjin Sandstorm Source Control Project was implemented from 2000,the vegetation in the study area had significantly restored and the first phase project had achieved remarkable results.But the vegetation in some areas still continued decreasing.So it is necessary to adopt the second phase project to consolidate and enhance the construction effect of the first phase project,and further improve the environment;(3)As for the future duration of NDVI change trends,the vegetation future change situation was not optimistic on the whole.In the four sandstorm source control regions,the vegetation future change trends and the project treatment effect in the pastoral transitional zone desertified land control region were better than those in the other three sandstorm source control regions.

    Beijing-Tianjin Sandstorm Source Control Region;normalized difference vegetation index;vegetation restoration;change trend;persistence

    Q146

    A

    1671-1114(2016)02-0047-07

    2015-11-09

    國家自然科學(xué)基金資助項目(41001022,31270510);國家科技支撐計劃資助項目(2012BAC07B02);教育部新世紀(jì)優(yōu)秀人才支持計劃資助項目(NCET-10-0954);天津師范大學(xué)市級重點實驗室開發(fā)研究基金資助項目(YF11700102).

    楊艷麗(1992—),女,碩士研究生.

    孫艷玲(1977—),女,副教授,主要從事資源環(huán)境遙感與全球變化方面的研究.

    猜你喜歡
    趨勢區(qū)域研究
    FMS與YBT相關(guān)性的實證研究
    遼代千人邑研究述論
    趨勢
    視錯覺在平面設(shè)計中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    初秋唇妝趨勢
    Coco薇(2017年9期)2017-09-07 21:23:49
    SPINEXPO?2017春夏流行趨勢
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點暫降頻次估計
    電測與儀表(2015年5期)2015-04-09 11:30:52
    日韩欧美 国产精品| 在线观看三级黄色| 边亲边吃奶的免费视频| 亚洲欧美精品自产自拍| 国产爽快片一区二区三区| 亚洲综合精品二区| 欧美性感艳星| 人妻制服诱惑在线中文字幕| 午夜免费男女啪啪视频观看| 亚洲图色成人| 欧美日韩av久久| 国产日韩一区二区三区精品不卡 | 欧美区成人在线视频| 亚洲人成网站在线观看播放| 亚洲国产精品国产精品| 日韩一区二区三区影片| 久久久久精品性色| 色视频在线一区二区三区| 黑人猛操日本美女一级片| 日韩三级伦理在线观看| 国产男人的电影天堂91| 国产免费一区二区三区四区乱码| 另类精品久久| 免费少妇av软件| 黄色配什么色好看| av在线观看视频网站免费| 少妇人妻 视频| 97在线视频观看| av国产久精品久网站免费入址| 天堂俺去俺来也www色官网| 成人二区视频| 亚洲精品一二三| 精品熟女少妇av免费看| www.色视频.com| 一级毛片久久久久久久久女| 免费在线观看成人毛片| 午夜免费男女啪啪视频观看| 亚洲精品亚洲一区二区| 亚洲国产欧美在线一区| 亚洲精品国产av成人精品| 晚上一个人看的免费电影| 赤兔流量卡办理| 欧美区成人在线视频| 看非洲黑人一级黄片| 伦理电影免费视频| 丝袜脚勾引网站| 青春草国产在线视频| 久久久久久久亚洲中文字幕| 永久网站在线| 免费不卡的大黄色大毛片视频在线观看| 免费观看无遮挡的男女| 国产亚洲一区二区精品| 在线播放无遮挡| 国产女主播在线喷水免费视频网站| 亚洲怡红院男人天堂| 青青草视频在线视频观看| 欧美日韩av久久| 精品人妻一区二区三区麻豆| av福利片在线观看| 亚洲怡红院男人天堂| 亚洲色图综合在线观看| 中文字幕制服av| 成人午夜精彩视频在线观看| 午夜日本视频在线| 精品一区在线观看国产| 97精品久久久久久久久久精品| 成人美女网站在线观看视频| 午夜福利在线观看免费完整高清在| 久久99一区二区三区| 99久久精品热视频| 久久午夜福利片| 国产精品蜜桃在线观看| 免费大片黄手机在线观看| 亚洲一区二区三区欧美精品| 97在线视频观看| 国产精品久久久久久av不卡| 免费看日本二区| 亚洲欧美成人精品一区二区| 一区二区三区精品91| 精品久久国产蜜桃| 国产片特级美女逼逼视频| 伊人久久国产一区二区| 国产av精品麻豆| 久久精品国产自在天天线| 偷拍熟女少妇极品色| 欧美日韩国产mv在线观看视频| 国产精品久久久久久久久免| 国产成人精品福利久久| 草草在线视频免费看| 校园人妻丝袜中文字幕| 免费看日本二区| 亚洲性久久影院| 亚洲精品乱码久久久久久按摩| 日韩中文字幕视频在线看片| 日韩亚洲欧美综合| 久久久a久久爽久久v久久| h日本视频在线播放| 国产熟女午夜一区二区三区 | 一区二区三区乱码不卡18| 丝袜在线中文字幕| 伊人久久精品亚洲午夜| 国产极品天堂在线| 综合色丁香网| 亚洲精品乱码久久久v下载方式| 99re6热这里在线精品视频| 国产精品国产三级专区第一集| 日韩大片免费观看网站| 精品国产国语对白av| 欧美成人午夜免费资源| 波野结衣二区三区在线| 亚洲国产精品专区欧美| 嫩草影院入口| 在线观看一区二区三区激情| 女性生殖器流出的白浆| 少妇人妻久久综合中文| 美女内射精品一级片tv| 男人狂女人下面高潮的视频| 在线精品无人区一区二区三| av网站免费在线观看视频| 中文精品一卡2卡3卡4更新| 最近手机中文字幕大全| 国产免费又黄又爽又色| 伊人亚洲综合成人网| 中文精品一卡2卡3卡4更新| 在现免费观看毛片| 秋霞伦理黄片| 免费高清在线观看视频在线观看| 国产亚洲最大av| 51国产日韩欧美| 日韩电影二区| 亚洲精品亚洲一区二区| 18禁在线播放成人免费| 特大巨黑吊av在线直播| 亚洲成人av在线免费| 国产一区二区在线观看av| 18禁在线无遮挡免费观看视频| 久久久久久人妻| 少妇高潮的动态图| 亚洲欧洲国产日韩| 欧美 亚洲 国产 日韩一| 国产91av在线免费观看| 欧美三级亚洲精品| 黑人高潮一二区| 亚洲va在线va天堂va国产| 极品少妇高潮喷水抽搐| 岛国毛片在线播放| 99热网站在线观看| 男人爽女人下面视频在线观看| 日韩视频在线欧美| 在线观看免费日韩欧美大片 | 成人无遮挡网站| 三上悠亚av全集在线观看 | 我要看日韩黄色一级片| 青春草国产在线视频| 久久精品国产亚洲av涩爱| 成人黄色视频免费在线看| 搡老乐熟女国产| 99热这里只有是精品在线观看| 国产探花极品一区二区| 亚洲国产av新网站| 在线观看免费高清a一片| 香蕉精品网在线| 国产一区二区三区综合在线观看 | 国产成人精品婷婷| 男男h啪啪无遮挡| 欧美成人午夜免费资源| 男女免费视频国产| 男女边吃奶边做爰视频| 久热久热在线精品观看| 最近中文字幕2019免费版| 国产精品国产三级国产av玫瑰| 99久久精品国产国产毛片| 一个人免费看片子| 日韩av不卡免费在线播放| 欧美性感艳星| 精品亚洲成a人片在线观看| 成年美女黄网站色视频大全免费 | 精品国产乱码久久久久久小说| 日日爽夜夜爽网站| 亚洲精品国产av蜜桃| 亚洲欧美成人精品一区二区| 国产精品秋霞免费鲁丝片| 成人无遮挡网站| tube8黄色片| 亚洲欧美一区二区三区黑人 | 国产一级毛片在线| 夫妻性生交免费视频一级片| 亚洲内射少妇av| 丝袜脚勾引网站| 免费看不卡的av| 国产精品成人在线| 性色avwww在线观看| videos熟女内射| 国产精品一区二区性色av| 男女国产视频网站| 午夜精品国产一区二区电影| 久久99一区二区三区| 成人特级av手机在线观看| 麻豆精品久久久久久蜜桃| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 欧美日韩国产mv在线观看视频| 日韩电影二区| 成人黄色视频免费在线看| videossex国产| 美女中出高潮动态图| 汤姆久久久久久久影院中文字幕| 校园人妻丝袜中文字幕| 中文字幕久久专区| 久久99蜜桃精品久久| 国产精品久久久久久av不卡| 免费人成在线观看视频色| 汤姆久久久久久久影院中文字幕| 美女内射精品一级片tv| 伊人久久精品亚洲午夜| 国产亚洲欧美精品永久| 精品亚洲乱码少妇综合久久| 男的添女的下面高潮视频| 国产高清国产精品国产三级| 特大巨黑吊av在线直播| 亚洲精品aⅴ在线观看| a级毛片免费高清观看在线播放| 免费看av在线观看网站| 永久网站在线| 人人澡人人妻人| 国产91av在线免费观看| 成人毛片a级毛片在线播放| 人人澡人人妻人| 日日啪夜夜爽| 成人毛片60女人毛片免费| 搡老乐熟女国产| 国产91av在线免费观看| 高清午夜精品一区二区三区| 狂野欧美激情性xxxx在线观看| 亚洲av综合色区一区| 国产在线免费精品| av又黄又爽大尺度在线免费看| 51国产日韩欧美| 国产伦在线观看视频一区| 看非洲黑人一级黄片| 人妻人人澡人人爽人人| 在线观看免费视频网站a站| 91aial.com中文字幕在线观看| 久久97久久精品| 久久久久久久亚洲中文字幕| 最新的欧美精品一区二区| 日本猛色少妇xxxxx猛交久久| 国产视频首页在线观看| 91精品国产国语对白视频| 亚洲av免费高清在线观看| 免费黄色在线免费观看| 亚洲图色成人| 成年女人在线观看亚洲视频| h日本视频在线播放| 少妇丰满av| 嫩草影院新地址| 黄色视频在线播放观看不卡| 精品久久久噜噜| 高清黄色对白视频在线免费看 | 一区二区三区四区激情视频| 狂野欧美白嫩少妇大欣赏| 观看免费一级毛片| 性高湖久久久久久久久免费观看| 国产 一区精品| 欧美精品国产亚洲| 成人无遮挡网站| 亚洲不卡免费看| 看免费成人av毛片| 亚洲一区二区三区欧美精品| 有码 亚洲区| 国内精品宾馆在线| av线在线观看网站| 大陆偷拍与自拍| 国国产精品蜜臀av免费| 久久精品国产鲁丝片午夜精品| 在线观看av片永久免费下载| 国精品久久久久久国模美| videos熟女内射| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 国产精品欧美亚洲77777| 啦啦啦啦在线视频资源| 亚洲av日韩在线播放| 综合色丁香网| av视频免费观看在线观看| 国产极品粉嫩免费观看在线 | 国语对白做爰xxxⅹ性视频网站| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 99久久精品国产国产毛片| 亚洲久久久国产精品| 色婷婷av一区二区三区视频| 欧美xxxx性猛交bbbb| 久久久久国产网址| 亚洲av男天堂| 春色校园在线视频观看| 伦精品一区二区三区| 亚洲欧美精品专区久久| 久久青草综合色| 乱码一卡2卡4卡精品| 欧美亚洲 丝袜 人妻 在线| av卡一久久| 少妇被粗大的猛进出69影院 | 久久精品久久久久久久性| av在线老鸭窝| 成人亚洲精品一区在线观看| 亚洲国产最新在线播放| 国产69精品久久久久777片| 一个人免费看片子| 肉色欧美久久久久久久蜜桃| av又黄又爽大尺度在线免费看| 午夜免费鲁丝| 韩国高清视频一区二区三区| 80岁老熟妇乱子伦牲交| 在线看a的网站| 草草在线视频免费看| 少妇裸体淫交视频免费看高清| 亚洲av综合色区一区| 国产一区二区三区综合在线观看 | 欧美日韩视频高清一区二区三区二| 男人狂女人下面高潮的视频| 九草在线视频观看| 国产精品秋霞免费鲁丝片| 嫩草影院入口| 老女人水多毛片| 国产免费福利视频在线观看| √禁漫天堂资源中文www| 麻豆乱淫一区二区| 久久人妻熟女aⅴ| av在线播放精品| 亚洲人成网站在线观看播放| 日韩大片免费观看网站| 亚洲国产精品国产精品| 日韩av在线免费看完整版不卡| 中文字幕av电影在线播放| av视频免费观看在线观看| 亚州av有码| 亚洲人成网站在线播| 亚洲精品乱码久久久久久按摩| 欧美 亚洲 国产 日韩一| 不卡视频在线观看欧美| 国产亚洲最大av| 女的被弄到高潮叫床怎么办| 黄片无遮挡物在线观看| 国产老妇伦熟女老妇高清| 色94色欧美一区二区| 色视频在线一区二区三区| 人妻 亚洲 视频| 日本vs欧美在线观看视频 | 亚洲精品日韩在线中文字幕| 春色校园在线视频观看| 久久久欧美国产精品| 日韩中文字幕视频在线看片| 亚洲不卡免费看| 国产精品一区二区性色av| 日日摸夜夜添夜夜爱| 欧美高清成人免费视频www| 爱豆传媒免费全集在线观看| 国产爽快片一区二区三区| 国产精品福利在线免费观看| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 高清不卡的av网站| 伦精品一区二区三区| 激情五月婷婷亚洲| 人妻 亚洲 视频| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产成人一精品久久久| 如日韩欧美国产精品一区二区三区 | 最黄视频免费看| 少妇的逼好多水| 老司机亚洲免费影院| av在线app专区| 精品久久久久久久久亚洲| 国产亚洲最大av| 另类精品久久| 亚洲第一区二区三区不卡| 少妇的逼好多水| 亚洲成人一二三区av| www.av在线官网国产| 丰满人妻一区二区三区视频av| 久久99热6这里只有精品| 全区人妻精品视频| 久久午夜综合久久蜜桃| 狠狠精品人妻久久久久久综合| 大话2 男鬼变身卡| 成人漫画全彩无遮挡| 黑人巨大精品欧美一区二区蜜桃 | 熟女av电影| 国产亚洲5aaaaa淫片| 热re99久久国产66热| 精品亚洲成a人片在线观看| 亚洲电影在线观看av| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av蜜桃| 永久免费av网站大全| 一区二区三区四区激情视频| 成年女人在线观看亚洲视频| 国产欧美日韩综合在线一区二区 | 国产一区有黄有色的免费视频| 最新的欧美精品一区二区| 午夜福利视频精品| 精品亚洲成国产av| 亚洲欧美一区二区三区黑人 | 亚洲,一卡二卡三卡| 亚洲av不卡在线观看| 亚洲经典国产精华液单| 国产成人免费观看mmmm| 国语对白做爰xxxⅹ性视频网站| 国产欧美日韩一区二区三区在线 | 国产高清三级在线| av福利片在线观看| 久久久久视频综合| 国产 精品1| 精品少妇黑人巨大在线播放| 九九久久精品国产亚洲av麻豆| 我要看黄色一级片免费的| 少妇的逼水好多| 亚洲美女视频黄频| 色婷婷久久久亚洲欧美| 蜜桃在线观看..| 亚洲欧美成人综合另类久久久| av不卡在线播放| 尾随美女入室| 精品一区二区免费观看| 久久人人爽人人片av| 免费av不卡在线播放| 夫妻午夜视频| 国产美女午夜福利| 亚洲激情五月婷婷啪啪| 日本黄色片子视频| 精品国产乱码久久久久久小说| 狠狠精品人妻久久久久久综合| 久久99蜜桃精品久久| av在线播放精品| 免费观看性生交大片5| 午夜91福利影院| 免费黄频网站在线观看国产| 日韩熟女老妇一区二区性免费视频| 亚洲熟女精品中文字幕| 在线观看免费高清a一片| 欧美97在线视频| 深夜a级毛片| 一本大道久久a久久精品| 国产成人午夜福利电影在线观看| 亚洲欧美成人综合另类久久久| 免费黄色在线免费观看| 国产精品一区二区三区四区免费观看| av又黄又爽大尺度在线免费看| 日本黄色片子视频| 日本-黄色视频高清免费观看| 成人漫画全彩无遮挡| 亚洲av二区三区四区| 免费观看性生交大片5| .国产精品久久| 美女cb高潮喷水在线观看| 老司机影院毛片| 99久久精品国产国产毛片| 在线观看av片永久免费下载| 日本黄色片子视频| 精品久久久久久久久亚洲| 国产日韩欧美亚洲二区| 国产国拍精品亚洲av在线观看| 久久午夜福利片| 日韩,欧美,国产一区二区三区| 午夜免费鲁丝| 尾随美女入室| 欧美性感艳星| h日本视频在线播放| 午夜91福利影院| 国产白丝娇喘喷水9色精品| 国产综合精华液| 国产白丝娇喘喷水9色精品| www.色视频.com| 丁香六月天网| 噜噜噜噜噜久久久久久91| 精品酒店卫生间| 2021少妇久久久久久久久久久| 大陆偷拍与自拍| 91精品国产国语对白视频| 一级毛片aaaaaa免费看小| 亚洲,欧美,日韩| 免费久久久久久久精品成人欧美视频 | 黄色怎么调成土黄色| 国产淫语在线视频| 亚洲怡红院男人天堂| 三上悠亚av全集在线观看 | 午夜免费观看性视频| 青春草亚洲视频在线观看| 国产成人精品福利久久| 色婷婷久久久亚洲欧美| 国产精品成人在线| 亚洲美女搞黄在线观看| 国产精品秋霞免费鲁丝片| 久久久久久久国产电影| 哪个播放器可以免费观看大片| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 五月天丁香电影| 男的添女的下面高潮视频| 乱系列少妇在线播放| 日韩av不卡免费在线播放| 欧美日韩av久久| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 美女大奶头黄色视频| 成年人午夜在线观看视频| 高清午夜精品一区二区三区| 久久韩国三级中文字幕| 国产成人精品福利久久| 国内精品宾馆在线| 久久国产乱子免费精品| 丁香六月天网| 国产精品女同一区二区软件| 精品国产一区二区三区久久久樱花| 久久精品夜色国产| 国产精品秋霞免费鲁丝片| 秋霞伦理黄片| 亚洲欧美日韩卡通动漫| 国产成人精品久久久久久| 国产视频内射| 日本vs欧美在线观看视频 | 最近中文字幕高清免费大全6| 欧美三级亚洲精品| 建设人人有责人人尽责人人享有的| 亚洲美女视频黄频| 午夜视频国产福利| 精品卡一卡二卡四卡免费| 欧美精品人与动牲交sv欧美| 国产精品蜜桃在线观看| 亚洲内射少妇av| 日本黄色日本黄色录像| 国产淫片久久久久久久久| 国产真实伦视频高清在线观看| 久久久久久久亚洲中文字幕| 日本av手机在线免费观看| 欧美三级亚洲精品| 男人狂女人下面高潮的视频| 亚洲精品乱久久久久久| 国产成人一区二区在线| 一级毛片我不卡| 免费人妻精品一区二区三区视频| 精品人妻一区二区三区麻豆| 久久午夜综合久久蜜桃| 欧美日韩一区二区视频在线观看视频在线| 另类精品久久| 又粗又硬又长又爽又黄的视频| 欧美日本中文国产一区发布| 简卡轻食公司| 免费播放大片免费观看视频在线观看| 国产在线男女| 亚洲婷婷狠狠爱综合网| 蜜臀久久99精品久久宅男| 国产乱人偷精品视频| 亚洲精品国产色婷婷电影| 国产精品一区二区在线观看99| 下体分泌物呈黄色| 国产熟女午夜一区二区三区 | 高清黄色对白视频在线免费看 | 日韩三级伦理在线观看| 亚洲欧美日韩东京热| 简卡轻食公司| 欧美 日韩 精品 国产| 大码成人一级视频| 精品少妇黑人巨大在线播放| 一本一本综合久久| 99热国产这里只有精品6| 在线观看免费视频网站a站| 六月丁香七月| 久久免费观看电影| 国产精品蜜桃在线观看| 少妇的逼好多水| av有码第一页| 春色校园在线视频观看| 国产精品国产三级专区第一集| 久久精品久久久久久久性| 精品人妻偷拍中文字幕| 视频中文字幕在线观看| 在线精品无人区一区二区三| 国产高清有码在线观看视频| 色94色欧美一区二区| 国模一区二区三区四区视频| 精品国产露脸久久av麻豆| 99久国产av精品国产电影| 国产成人精品无人区| 少妇 在线观看| 国产精品福利在线免费观看| 中文字幕制服av| 日韩不卡一区二区三区视频在线| 免费黄频网站在线观看国产| 又粗又硬又长又爽又黄的视频| 18+在线观看网站| 国产精品女同一区二区软件| 久久久久久久久久久免费av| 2021少妇久久久久久久久久久| 男女边摸边吃奶| 国产无遮挡羞羞视频在线观看| 久久国产乱子免费精品| 亚洲国产精品国产精品| 大码成人一级视频| 另类亚洲欧美激情| 校园人妻丝袜中文字幕| 久久久久久久久久人人人人人人| 国产精品熟女久久久久浪| 91午夜精品亚洲一区二区三区| 男人狂女人下面高潮的视频| 国产成人91sexporn| 中文字幕久久专区| 国产高清国产精品国产三级| 精品少妇黑人巨大在线播放| 亚洲成人手机| 人妻夜夜爽99麻豆av| 3wmmmm亚洲av在线观看| 91精品国产九色| 亚洲欧美成人综合另类久久久|