• <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)鍵詞:趨勢(shì)區(qū)域研究

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

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

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

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

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

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

    京津風(fēng)沙源治理區(qū);歸一化植被指數(shù);植被恢復(fù);變化趨勢(shì);可持續(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ù)來(lái)源[3].歸一化植被指數(shù)(normalized difference vegetation index,NDVI)具有良好的植被信息表達(dá)能力和數(shù)據(jù)提取過(guò)程中較強(qiáng)的抗干擾能力,在植被變化研究中應(yīng)用廣泛[4].目前常用的NDVI數(shù)據(jù)包括SPOT VGT/NDVI、MODIS/NDVI和GIMMS/NDVI等[5].其中,GIMMS/NDVI數(shù)據(jù)時(shí)間序列最長(zhǎng),覆蓋范圍廣,因此在長(zhǎng)時(shí)期、大尺度的植被宏觀動(dòng)態(tài)監(jiān)測(cè)中應(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ù),采用線性趨勢(shì)法、Sen+Mann-Kendall方法和Hurst指數(shù)對(duì)不同時(shí)空尺度上的植被變化趨勢(shì)進(jìn)行了探討分析[1-2].

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

    1 資料來(lái)源和研究方法

    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個(gè)縣(旗、市、區(qū)),總面積45.8萬(wàn)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來(lái),本區(qū)氣候呈現(xiàn)暖干化趨勢(shì),干旱區(qū)和干旱亞濕潤(rùn)區(qū)的界限逐步向東推進(jìn)[16].京津風(fēng)沙源治理工程于2000年試點(diǎn),2002年全面實(shí)施,到2010年為止,工程區(qū)共完成退耕還林及造林6.0萬(wàn)km2,草地治理8.67萬(wàn)km2,生態(tài)環(huán)境得到顯著改善.按照國(guó)家規(guī)劃,工程治理區(qū)分為北部干旱草原治理區(qū)、農(nóng)牧交錯(cuò)帶沙化土地治理區(qū)、燕山丘陵山地水源保護(hù)區(qū)和渾善達(dá)克沙地治理區(qū).

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

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

    1.3研究方法

    1.3.1趨勢(shì)分析

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

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

    1.3.2可持續(xù)性分析

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

    其中,1982—2000年間的p(x,y)值與該時(shí)間段的s(x,y)值相等,其他時(shí)間段的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,那么該時(shí)段的p(x,y)=0.以此類推.

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

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

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

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

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

    2 結(jié)果與分析

    2.1NDVI時(shí)間變化特征

    通過(guò)分別統(tǒng)計(jì)京津風(fēng)沙源治理區(qū)及4大分區(qū)2001—2011年的逐年NDVI值,制作NDVI年際變化圖,分析其植被的時(shí)間變化趨勢(shì),結(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ū)植被整體上呈增加趨勢(shì),平均增速為每10 a增加0.02.4大分區(qū)在2001—2011年間植被也都呈增加趨勢(shì),但增速存在區(qū)域差異.其中,燕山丘陵山地水源保護(hù)區(qū)的植被增加速度最快,平均增速為每10 a增加0.05;農(nóng)牧交錯(cuò)帶沙化土地治理區(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ù)性相加,計(jì)算得到NDVI變化趨勢(shì)的持續(xù)性疊加圖,以此來(lái)反映研究區(qū)2001—2011年間植被變化趨勢(shì)持續(xù)性的時(shí)空特征,結(jié)果如圖3所示.

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

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

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

    2.3NDVI變化趨勢(shì)的未來(lái)持續(xù)時(shí)間特征

    2.3.1整個(gè)研究區(qū)NDVI變化趨勢(shì)的未來(lái)持續(xù)時(shí)間

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

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

    2.3.24大治理區(qū)NDVI變化趨勢(shì)的未來(lái)持續(xù)時(shí)間

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

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

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

    3.2結(jié)論

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

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

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

    [1]楊林山,李常斌,王帥兵,等.氣候變化和人類活動(dòng)對(duì)洮河流域植被動(dòng)態(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é)報(bào),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)境時(shí)空變化的遙感監(jiān)測(cè)分析[J].地球信息科學(xué)學(xué)報(bào),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]呂洋,董國(guó)濤,楊勝天,等.雅魯藏布江流域NDVI時(shí)空變化及其與降水和高程的關(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的中國(guó)東部植被指數(shù)對(duì)比分析[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ù)時(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]范德芹,朱文泉,潘耀忠,等.青藏高原小嵩草高寒草甸返青期遙感識(shí)別方法篩選[J].遙感學(xué)報(bào),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ū)植被覆蓋時(shí)空演變特征[J].生態(tài)學(xué)報(bào),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ū)植被的時(shí)空變化及其對(duì)氣候因子的響應(yīng)[J].中國(guó)沙漠,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

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

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

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

    猜你喜歡
    趨勢(shì)區(qū)域研究
    FMS與YBT相關(guān)性的實(shí)證研究
    遼代千人邑研究述論
    趨勢(shì)
    視錯(cuò)覺在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    初秋唇妝趨勢(shì)
    Coco薇(2017年9期)2017-09-07 21:23:49
    SPINEXPO?2017春夏流行趨勢(shì)
    關(guān)于四色猜想
    分區(qū)域
    基于嚴(yán)重區(qū)域的多PCC點(diǎn)暫降頻次估計(jì)
    国产精品久久久久久精品电影| 精品第一国产精品| 久久久精品欧美日韩精品| 成在线人永久免费视频| 岛国在线免费视频观看| 男女下面进入的视频免费午夜| 久久中文看片网| 亚洲国产欧洲综合997久久,| 制服丝袜大香蕉在线| 一级片免费观看大全| a在线观看视频网站| 免费在线观看黄色视频的| 真人一进一出gif抽搐免费| 免费在线观看黄色视频的| 亚洲avbb在线观看| 欧美一级a爱片免费观看看 | 午夜福利18| 午夜福利18| 久久香蕉精品热| 国产精品98久久久久久宅男小说| 成人永久免费在线观看视频| 国产激情偷乱视频一区二区| 窝窝影院91人妻| 色综合婷婷激情| 搞女人的毛片| 亚洲精品久久成人aⅴ小说| 国产伦人伦偷精品视频| 99热6这里只有精品| 91麻豆av在线| 成人国产综合亚洲| 99riav亚洲国产免费| 一个人观看的视频www高清免费观看 | 国产精品av久久久久免费| 女生性感内裤真人,穿戴方法视频| 99热只有精品国产| 五月伊人婷婷丁香| 久久久国产欧美日韩av| 欧美日本亚洲视频在线播放| av国产免费在线观看| 免费无遮挡裸体视频| 亚洲精品中文字幕在线视频| 国产精品久久久久久精品电影| 色播亚洲综合网| 老司机福利观看| 国产精品乱码一区二三区的特点| 久久精品国产亚洲av香蕉五月| 99久久无色码亚洲精品果冻| 久久国产精品人妻蜜桃| 99久久久亚洲精品蜜臀av| 波多野结衣高清作品| 午夜精品在线福利| 久久午夜综合久久蜜桃| 欧美精品亚洲一区二区| 一a级毛片在线观看| 成在线人永久免费视频| 青草久久国产| 三级国产精品欧美在线观看 | 国产成人精品无人区| 日本a在线网址| 欧美性猛交╳xxx乱大交人| 久久久久久久精品吃奶| 色精品久久人妻99蜜桃| 国内毛片毛片毛片毛片毛片| 亚洲狠狠婷婷综合久久图片| av视频在线观看入口| 久久人妻福利社区极品人妻图片| 性色av乱码一区二区三区2| 国产黄色小视频在线观看| 久久精品成人免费网站| 美女 人体艺术 gogo| 亚洲五月婷婷丁香| 又紧又爽又黄一区二区| 欧美日韩黄片免| 亚洲专区中文字幕在线| 99久久久亚洲精品蜜臀av| 久久久久久免费高清国产稀缺| 一进一出抽搐gif免费好疼| 丁香欧美五月| av有码第一页| 亚洲第一欧美日韩一区二区三区| 亚洲av五月六月丁香网| 精品乱码久久久久久99久播| 热99re8久久精品国产| 在线观看舔阴道视频| 人人妻,人人澡人人爽秒播| 一个人免费在线观看的高清视频| 99在线视频只有这里精品首页| 少妇的丰满在线观看| e午夜精品久久久久久久| 女人高潮潮喷娇喘18禁视频| 久久久久久久久久黄片| 男女视频在线观看网站免费 | 手机成人av网站| 亚洲精品在线美女| 一级a爱片免费观看的视频| 哪里可以看免费的av片| 一个人观看的视频www高清免费观看 | 精品乱码久久久久久99久播| 一级片免费观看大全| 首页视频小说图片口味搜索| 看免费av毛片| 国产亚洲av高清不卡| 天天躁夜夜躁狠狠躁躁| 久久久久亚洲av毛片大全| av中文乱码字幕在线| av在线播放免费不卡| 国产一区二区在线av高清观看| 99在线人妻在线中文字幕| 久久九九热精品免费| 免费观看精品视频网站| 嫩草影视91久久| 久久中文字幕一级| 日本三级黄在线观看| 日韩中文字幕欧美一区二区| 亚洲狠狠婷婷综合久久图片| 国产久久久一区二区三区| 国产精品久久久久久久电影 | 色综合亚洲欧美另类图片| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| av欧美777| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 亚洲国产精品999在线| 欧美在线一区亚洲| 天堂影院成人在线观看| 成人国产一区最新在线观看| 亚洲黑人精品在线| 男女那种视频在线观看| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 窝窝影院91人妻| 村上凉子中文字幕在线| 一本久久中文字幕| 91九色精品人成在线观看| av国产免费在线观看| 岛国在线免费视频观看| 毛片女人毛片| 国产精品影院久久| 老熟妇仑乱视频hdxx| 巨乳人妻的诱惑在线观看| 国产精品自产拍在线观看55亚洲| 国产熟女午夜一区二区三区| 最新美女视频免费是黄的| 在线a可以看的网站| 国内揄拍国产精品人妻在线| 丝袜美腿诱惑在线| 美女 人体艺术 gogo| 看片在线看免费视频| 国产精品亚洲一级av第二区| 久久久久久人人人人人| 99精品欧美一区二区三区四区| 在线观看午夜福利视频| 国产激情久久老熟女| 波多野结衣巨乳人妻| 欧美乱色亚洲激情| 亚洲色图av天堂| 国产精品影院久久| 琪琪午夜伦伦电影理论片6080| 麻豆成人午夜福利视频| 日韩欧美国产在线观看| 色老头精品视频在线观看| 午夜免费成人在线视频| 精品国产美女av久久久久小说| a级毛片在线看网站| 国产av一区在线观看免费| 欧美不卡视频在线免费观看 | 最近最新免费中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产看品久久| 男女之事视频高清在线观看| 桃色一区二区三区在线观看| 中文字幕久久专区| 手机成人av网站| 999久久久国产精品视频| 国产欧美日韩一区二区三| av有码第一页| 一个人观看的视频www高清免费观看 | 婷婷精品国产亚洲av在线| 99在线人妻在线中文字幕| 一a级毛片在线观看| 99久久久亚洲精品蜜臀av| 亚洲成人中文字幕在线播放| 九色成人免费人妻av| 最近最新中文字幕大全免费视频| 好男人在线观看高清免费视频| 国产精品亚洲av一区麻豆| 国产伦在线观看视频一区| 国内少妇人妻偷人精品xxx网站 | av福利片在线| 一级黄色大片毛片| 少妇裸体淫交视频免费看高清 | 国产不卡一卡二| 成年人黄色毛片网站| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美一区二区三区在线观看| 两人在一起打扑克的视频| 熟女电影av网| 欧美zozozo另类| 观看免费一级毛片| 国产又色又爽无遮挡免费看| 国产日本99.免费观看| 小说图片视频综合网站| 男女之事视频高清在线观看| 成人精品一区二区免费| 亚洲成人国产一区在线观看| 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 1024视频免费在线观看| 国产免费男女视频| 琪琪午夜伦伦电影理论片6080| 国产成人啪精品午夜网站| 日本免费一区二区三区高清不卡| 欧美黑人巨大hd| 99国产精品99久久久久| 国产91精品成人一区二区三区| aaaaa片日本免费| 熟女电影av网| 99久久无色码亚洲精品果冻| 熟妇人妻久久中文字幕3abv| 色综合站精品国产| 99热这里只有是精品50| 久久人人精品亚洲av| 午夜福利视频1000在线观看| 搞女人的毛片| 亚洲人成77777在线视频| 日韩欧美国产在线观看| 色在线成人网| 一级作爱视频免费观看| 亚洲人成网站在线播放欧美日韩| 叶爱在线成人免费视频播放| 国产精品野战在线观看| 黄片小视频在线播放| 狂野欧美激情性xxxx| 999精品在线视频| cao死你这个sao货| 国产成+人综合+亚洲专区| 久久久久久久久久黄片| 国产精品永久免费网站| 成人三级黄色视频| 欧美色欧美亚洲另类二区| 国产99白浆流出| 久久中文字幕人妻熟女| 精品国产乱子伦一区二区三区| 99re在线观看精品视频| 成年免费大片在线观看| 欧美日韩黄片免| 免费人成视频x8x8入口观看| 搞女人的毛片| 国产亚洲欧美在线一区二区| 日韩大尺度精品在线看网址| 国产私拍福利视频在线观看| 女人被狂操c到高潮| 日韩欧美 国产精品| 国产精品 国内视频| 亚洲成av人片在线播放无| 深夜精品福利| 精品久久久久久久久久久久久| 亚洲乱码一区二区免费版| 18禁裸乳无遮挡免费网站照片| 久久精品亚洲精品国产色婷小说| 嫩草影院精品99| 国产精品一区二区免费欧美| 一进一出好大好爽视频| 精品欧美一区二区三区在线| 18禁观看日本| 精品熟女少妇八av免费久了| 不卡一级毛片| 男插女下体视频免费在线播放| 又黄又爽又免费观看的视频| 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添小说| 香蕉av资源在线| 久久草成人影院| ponron亚洲| 桃红色精品国产亚洲av| 小说图片视频综合网站| 两个人视频免费观看高清| 中文字幕av在线有码专区| 天堂动漫精品| 99热这里只有是精品50| 岛国视频午夜一区免费看| 亚洲男人的天堂狠狠| 欧美又色又爽又黄视频| 一级作爱视频免费观看| 免费在线观看影片大全网站| www.www免费av| 岛国在线观看网站| 日日干狠狠操夜夜爽| 色综合婷婷激情| 老汉色av国产亚洲站长工具| 黄色a级毛片大全视频| 少妇被粗大的猛进出69影院| 国产精品 欧美亚洲| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类 | 免费观看人在逋| 精品午夜福利视频在线观看一区| 精品久久久久久,| 亚洲人与动物交配视频| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 少妇的丰满在线观看| 又粗又爽又猛毛片免费看| 我要搜黄色片| 日韩大尺度精品在线看网址| 天堂√8在线中文| 欧美大码av| 啦啦啦免费观看视频1| 一本一本综合久久| 在线播放国产精品三级| 午夜精品久久久久久毛片777| 日韩欧美国产在线观看| 一进一出抽搐gif免费好疼| 哪里可以看免费的av片| 很黄的视频免费| 亚洲欧美日韩东京热| 亚洲片人在线观看| 丝袜人妻中文字幕| 亚洲午夜理论影院| 草草在线视频免费看| 一边摸一边抽搐一进一小说| 国产精品永久免费网站| 麻豆久久精品国产亚洲av| 国产一区二区三区在线臀色熟女| 欧美3d第一页| 亚洲国产精品合色在线| 免费av毛片视频| 1024手机看黄色片| 麻豆一二三区av精品| 欧美日本亚洲视频在线播放| 叶爱在线成人免费视频播放| 成人三级做爰电影| 欧美成人午夜精品| 亚洲av美国av| 欧美丝袜亚洲另类 | 亚洲人成伊人成综合网2020| 99久久99久久久精品蜜桃| 国内揄拍国产精品人妻在线| 后天国语完整版免费观看| 好看av亚洲va欧美ⅴa在| 午夜免费观看网址| 色综合亚洲欧美另类图片| 中文字幕人妻丝袜一区二区| 国产真实乱freesex| 亚洲一区二区三区不卡视频| 欧美在线黄色| 嫁个100分男人电影在线观看| 国产av麻豆久久久久久久| 精品人妻1区二区| 免费观看人在逋| 亚洲精品国产一区二区精华液| 欧美三级亚洲精品| 91成年电影在线观看| 精品久久久久久成人av| 不卡一级毛片| 国产精品爽爽va在线观看网站| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 欧美性猛交╳xxx乱大交人| 亚洲中文日韩欧美视频| 叶爱在线成人免费视频播放| 欧美性长视频在线观看| 一个人观看的视频www高清免费观看 | 亚洲男人的天堂狠狠| 国产不卡一卡二| 色av中文字幕| 国产精品久久久人人做人人爽| 国产久久久一区二区三区| 亚洲片人在线观看| 中文在线观看免费www的网站 | 久久久久久久久中文| 亚洲无线在线观看| 久久精品亚洲精品国产色婷小说| 亚洲av熟女| 757午夜福利合集在线观看| 欧美zozozo另类| 亚洲午夜精品一区,二区,三区| 一夜夜www| 欧美日韩瑟瑟在线播放| 国产精品久久视频播放| 妹子高潮喷水视频| 成人av在线播放网站| 日本免费一区二区三区高清不卡| xxxwww97欧美| 亚洲国产精品999在线| 精品日产1卡2卡| 一a级毛片在线观看| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 欧美中文综合在线视频| 黄片小视频在线播放| 精品久久久久久久久久免费视频| 国产三级黄色录像| 91九色精品人成在线观看| 我要搜黄色片| 国产黄片美女视频| 在线观看免费日韩欧美大片| 一进一出抽搐gif免费好疼| 一本一本综合久久| 午夜激情av网站| 久久久久亚洲av毛片大全| 中文字幕人成人乱码亚洲影| 久久性视频一级片| 国产免费男女视频| 天堂√8在线中文| 中文在线观看免费www的网站 | 国产精品美女特级片免费视频播放器 | 在线国产一区二区在线| 日韩欧美三级三区| 亚洲欧美精品综合久久99| 999久久久精品免费观看国产| 成人午夜高清在线视频| 一本综合久久免费| 久久久久久久久久黄片| 国内精品久久久久久久电影| 一级黄色大片毛片| 日本a在线网址| 免费一级毛片在线播放高清视频| 国产一区在线观看成人免费| 欧美色视频一区免费| 两性午夜刺激爽爽歪歪视频在线观看 | 国模一区二区三区四区视频 | 免费无遮挡裸体视频| 亚洲狠狠婷婷综合久久图片| 一本精品99久久精品77| 女人高潮潮喷娇喘18禁视频| 欧美色视频一区免费| 国产精品99久久99久久久不卡| 在线观看美女被高潮喷水网站 | 男男h啪啪无遮挡| 99精品欧美一区二区三区四区| 精品国内亚洲2022精品成人| 丝袜美腿诱惑在线| 国产精品98久久久久久宅男小说| 成年人黄色毛片网站| 男人舔女人的私密视频| 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 国产探花在线观看一区二区| 久久中文字幕一级| 亚洲午夜精品一区,二区,三区| 日韩大码丰满熟妇| 91字幕亚洲| 国产单亲对白刺激| 久久天堂一区二区三区四区| 国产精品1区2区在线观看.| 叶爱在线成人免费视频播放| 国产麻豆成人av免费视频| 久久精品亚洲精品国产色婷小说| 五月伊人婷婷丁香| 久久精品91无色码中文字幕| 久久午夜综合久久蜜桃| 黑人巨大精品欧美一区二区mp4| 18禁黄网站禁片免费观看直播| 18禁裸乳无遮挡免费网站照片| 成人精品一区二区免费| 欧美+亚洲+日韩+国产| 两性夫妻黄色片| 亚洲国产高清在线一区二区三| 黄片大片在线免费观看| 欧美三级亚洲精品| 后天国语完整版免费观看| 亚洲成人免费电影在线观看| 久久精品成人免费网站| 国产亚洲av高清不卡| 在线观看免费日韩欧美大片| 国产1区2区3区精品| 精品久久久久久久毛片微露脸| 波多野结衣巨乳人妻| 亚洲无线在线观看| 欧美乱码精品一区二区三区| 一个人免费在线观看的高清视频| 亚洲国产日韩欧美精品在线观看 | 久久久久久久久免费视频了| 国产私拍福利视频在线观看| 制服人妻中文乱码| 日本精品一区二区三区蜜桃| 国内久久婷婷六月综合欲色啪| 亚洲激情在线av| 最新美女视频免费是黄的| 国产精品,欧美在线| 两个人看的免费小视频| 精品乱码久久久久久99久播| 日韩大尺度精品在线看网址| 亚洲精品美女久久av网站| 午夜精品一区二区三区免费看| 变态另类丝袜制服| 中文字幕av在线有码专区| 天天添夜夜摸| 国产亚洲精品av在线| 国产黄色小视频在线观看| 欧美乱色亚洲激情| 国产成人精品久久二区二区91| 又黄又爽又免费观看的视频| 中亚洲国语对白在线视频| 国产黄片美女视频| 久久精品91无色码中文字幕| 熟女电影av网| 两性夫妻黄色片| 亚洲 欧美 日韩 在线 免费| 国产精华一区二区三区| 国产午夜福利久久久久久| a在线观看视频网站| 色哟哟哟哟哟哟| 成人av在线播放网站| 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 日本五十路高清| 最新在线观看一区二区三区| 久久人人精品亚洲av| 亚洲第一欧美日韩一区二区三区| 亚洲国产看品久久| 日本三级黄在线观看| 男人舔女人下体高潮全视频| 一级作爱视频免费观看| 久久婷婷成人综合色麻豆| 久久九九热精品免费| ponron亚洲| 成熟少妇高潮喷水视频| 国产午夜精品久久久久久| 少妇熟女aⅴ在线视频| 成人手机av| 午夜影院日韩av| 久久人妻av系列| 欧美日韩亚洲综合一区二区三区_| 一边摸一边做爽爽视频免费| 亚洲成a人片在线一区二区| 久久精品国产清高在天天线| 黑人操中国人逼视频| 免费看十八禁软件| 亚洲精品一区av在线观看| 制服人妻中文乱码| 香蕉丝袜av| 精品一区二区三区视频在线观看免费| 国产69精品久久久久777片 | 99国产综合亚洲精品| 亚洲欧美精品综合一区二区三区| 非洲黑人性xxxx精品又粗又长| 久久这里只有精品中国| 日日夜夜操网爽| 亚洲精品中文字幕一二三四区| 人妻夜夜爽99麻豆av| 久久草成人影院| 国产av不卡久久| 亚洲最大成人中文| 在线观看免费午夜福利视频| 一区福利在线观看| 久久 成人 亚洲| 搡老妇女老女人老熟妇| 少妇粗大呻吟视频| 亚洲av电影不卡..在线观看| 麻豆国产97在线/欧美 | 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 首页视频小说图片口味搜索| 亚洲av成人av| 老司机在亚洲福利影院| 久久这里只有精品19| 国产精品一及| 国产av麻豆久久久久久久| 日本 欧美在线| 国产三级中文精品| 亚洲天堂国产精品一区在线| 免费观看人在逋| 老司机靠b影院| 成人亚洲精品av一区二区| 琪琪午夜伦伦电影理论片6080| 在线a可以看的网站| 两人在一起打扑克的视频| 欧美黑人欧美精品刺激| 亚洲一区中文字幕在线| 99久久国产精品久久久| 丰满的人妻完整版| 亚洲专区字幕在线| 麻豆成人午夜福利视频| 国产av一区二区精品久久| 色综合婷婷激情| 亚洲五月婷婷丁香| 丁香欧美五月| 国产激情欧美一区二区| 国产精品 国内视频| 男女午夜视频在线观看| 啦啦啦观看免费观看视频高清| 国产精品久久久久久人妻精品电影| 亚洲av成人不卡在线观看播放网| 九色国产91popny在线| 成人午夜高清在线视频| 日韩欧美国产一区二区入口| 男女那种视频在线观看| 欧美在线一区亚洲| av天堂在线播放| 国内精品久久久久精免费| 免费搜索国产男女视频| 亚洲精品一区av在线观看| 久久久精品大字幕| 国产精品久久久久久人妻精品电影| av天堂在线播放| 怎么达到女性高潮| 最近最新免费中文字幕在线| 亚洲成人精品中文字幕电影| 国产精品野战在线观看| 午夜免费激情av| 国产又色又爽无遮挡免费看| 91在线观看av| 99在线视频只有这里精品首页| 99久久精品热视频| 老司机深夜福利视频在线观看| 久99久视频精品免费| 久久天堂一区二区三区四区| 在线十欧美十亚洲十日本专区| 久久久精品国产亚洲av高清涩受| 99精品欧美一区二区三区四区| 国产乱人伦免费视频|