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

    甘肅祁連山自然保護(hù)區(qū)草地時(shí)空變化及其對(duì)氣候的響應(yīng)

    2019-10-22 02:28:08宋偉宏王莉娜張金龍
    草業(yè)科學(xué) 2019年9期
    關(guān)鍵詞:祁連山草甸自然保護(hù)區(qū)

    宋偉宏,王莉娜,張金龍

    (1.甘肅省環(huán)境監(jiān)測(cè)中心站,甘肅 蘭州 730020;2.中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院/甘肅省遙感重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000)

    祁連山是由森林、草地交錯(cuò)分布組成的復(fù)合生態(tài)系統(tǒng)[1-3]。草地作為祁連山生態(tài)系統(tǒng)的重要組成部分,在維持區(qū)域生態(tài)平衡,阻止騰格里、巴丹吉林和庫(kù)姆塔格沙漠南侵,保障黃河和內(nèi)陸河徑流補(bǔ)給等方面具有不可忽視的作用[4-5]。同時(shí),草地也是祁連山地區(qū)畜牧業(yè)發(fā)展的重要生產(chǎn)生活資料,與牧民的家庭經(jīng)濟(jì)收益及區(qū)域經(jīng)濟(jì)發(fā)展密不可分[6-7]。及時(shí)掌握祁連山草地植被生長(zhǎng)狀況、演變趨勢(shì)及其與氣候的關(guān)系,對(duì)區(qū)域復(fù)合生態(tài)系統(tǒng)功能的科學(xué)評(píng)估具有參考價(jià)值,對(duì)指導(dǎo)地方生態(tài)環(huán)境系統(tǒng)修復(fù)治理具有重要借鑒意義。

    遙感技術(shù)以其速度快、范圍廣、連續(xù)性好和成本低等優(yōu)點(diǎn)被廣泛應(yīng)用在草地變化的監(jiān)測(cè)中[8-11]。Akiyama等[12]曾提出遙感監(jiān)測(cè)結(jié)合地理信息系統(tǒng)在草地監(jiān)測(cè)方面很有前景,而且隨著衛(wèi)星傳感器的發(fā)展和地理信息系統(tǒng)的推廣,精準(zhǔn)、實(shí)時(shí)的草地監(jiān)測(cè)也將會(huì)成為可能。Zhao等[13]基于遙感技術(shù),實(shí)現(xiàn)了對(duì)陰山北麓草地生態(tài)功能區(qū)植被覆蓋率的動(dòng)態(tài)監(jiān)測(cè)。馮雙雙[8]基于Landsat TM和Landsat 8 OLI TIRS影像建立了草地植被退化遙感監(jiān)測(cè)模型,實(shí)現(xiàn)對(duì)壩上草原草地退化的動(dòng)態(tài)監(jiān)測(cè)。Eckert等[14]研究發(fā)現(xiàn)MODIS NDVI時(shí)間序列分析對(duì)于探測(cè)植被變化區(qū)域和識(shí)別土地退化和再生是適宜的。葛靜等[15]基于UAV和MODIS遙感數(shù)據(jù)對(duì)黃河源東部地區(qū)的高寒草地蓋度動(dòng)態(tài)進(jìn)行了監(jiān)測(cè)。李肖娟[16]基于遙感數(shù)據(jù)對(duì)祁連山NPP時(shí)空變化特征和趨勢(shì)做了分析。Shao等[17]基于多期Landsat衛(wèi)星圖像對(duì)青藏高原阿壩縣草地保護(hù)工程前后的草地變化進(jìn)行了監(jiān)測(cè),從中得到的草地退化和再生的變化規(guī)律有效地輔助青藏高原環(huán)境恢復(fù)措施的制定。氣候變化對(duì)植被的影響已經(jīng)被國(guó)內(nèi)外眾多學(xué)者關(guān)注并展開研究,且主要體現(xiàn)在植被NDVI與氣溫、降水的關(guān)系和響應(yīng)特征上[18-21]。Gu等[22]基于MOD13Q1 NDVI數(shù)據(jù)對(duì)紅河流域的植被覆蓋與氣候因子的響應(yīng)進(jìn)行研究。Duo等[23]對(duì)華北平原過(guò)去33年的植被覆蓋度進(jìn)行了監(jiān)測(cè)。有研究表明,降水在影響西北地區(qū)植被變化的自然因素中起著最主要作用[24]。張禹舜[25]基于MODIS遙感數(shù)據(jù)研究了近11年祁連山年際、年內(nèi)、生長(zhǎng)季及不同季節(jié)植被NPP變化及其與氣候因子的響應(yīng)。Guo等[26]通過(guò)重構(gòu)SPOT_Vegetation NDVI數(shù)據(jù)集實(shí)現(xiàn)對(duì)中國(guó)北方錫林郭勒草原生長(zhǎng)氣的遙感監(jiān)測(cè)并分析了與氣象因子的關(guān)系。王巧玲等[1]對(duì)祁連山北麓草原植被NDVI與氣候的關(guān)系研究發(fā)現(xiàn),天祝縣和肅北縣的NDVI分別與氣溫和年降水量顯著正相關(guān)。

    本研究以祁連山自然保護(hù)區(qū)為研究區(qū),在前人研究的基礎(chǔ)上結(jié)合多期土地覆被數(shù)據(jù)、植被指數(shù)和氣象資料,分析祁連山自然保護(hù)區(qū)內(nèi)不同草原類型NDVI空間變化特征及趨勢(shì),并從季節(jié)尺度分析草地演變對(duì)氣溫和降水的響應(yīng)特征,以期進(jìn)一步認(rèn)識(shí)該區(qū)域氣候變化對(duì)草地演變的影響程度。

    1 材料與方法

    1.1 研究區(qū)概括

    甘肅祁連山自然保護(hù)區(qū)位于甘肅省境內(nèi)祁連山自然保護(hù)區(qū)中、東段,地跨武威、金昌、張掖3市 8 個(gè)縣區(qū),地理位置 97°23′34″-103°45′49″ E,36°29′57″-39°43′39″ N,北至祁連山森林分布下線與河西走廊相鄰,南沿祁連山主脈與青海省接壤,西至肅南縣界,東至天??h界,平均海拔4 000-4 500 m。由于特殊的外營(yíng)力作用及復(fù)雜的地質(zhì)發(fā)展期,使得區(qū)域內(nèi)的地形與氣候形成了東西南北向差異明顯的地貌形態(tài),植被和土壤的水平與垂直地帶性規(guī)律表現(xiàn)明顯[27],自上而下分為5個(gè)生物氣候帶,分別為山地荒漠草原氣候帶、山地草原氣候帶、山地森林草原氣候帶、亞高山灌叢草甸氣候帶和高山亞冰雪稀疏植被氣候帶。在涵養(yǎng)水源、維護(hù)區(qū)域生態(tài)平衡方面具有重要地位的祁連山區(qū)通過(guò)其內(nèi)部的森林、草地等生態(tài)系統(tǒng)構(gòu)成的復(fù)合生態(tài)系統(tǒng)調(diào)蓄區(qū)域降水、冰川和積雪融水。

    區(qū)域內(nèi)全年降水量主要集中在5-9月,年均降水量257.2~389.9 mm,年均降水量隨海拔高度升高增多,在達(dá)到最大降水高度后,又隨海拔高度升高而逐漸減少。祁連山區(qū)大部分地區(qū)年平均氣溫1~4 ℃,林區(qū)極端最高氣溫34 ℃,極端最低氣溫-30.2 ℃?!?10 ℃活動(dòng)積溫在海拔2 500 m以下地區(qū)為1 500~2 500 ℃·d,海拔3 000 m以上地區(qū)為200~500 ℃·d[28]。研究區(qū)具體地理位置如圖1所示。

    1.2 數(shù)據(jù)源與分析方法

    1.2.1 數(shù)據(jù)源及處理

    MODIS NDVI數(shù)據(jù)來(lái)源于美國(guó)地質(zhì)勘探局(https://lpdaac.usgs.gov)的MOD13Q1,該產(chǎn)品是一個(gè)采用Sinusoidal投影方式的3級(jí)陸地植被數(shù)據(jù),其空間分辨率為250 m,時(shí)間分辨率為16 d。本研究選用2000-2018年間共414景影像,基于MRT工具對(duì)其進(jìn)行拼接、波段提取及投影轉(zhuǎn)換,并利用祁連山自然保護(hù)區(qū)矢量邊界對(duì)其裁剪,再通過(guò)ArcGIS軟件提取2000-2018年16 d合成的NDVI數(shù)據(jù)。

    土地覆被數(shù)據(jù):基于國(guó)家土地覆被分類標(biāo)準(zhǔn),將祁連山自然保護(hù)區(qū)土地覆被分為耕地、林地、草地、水域、城鄉(xiāng)工礦居民用地和未利用地六大類。Landsat影像(表1)來(lái)自于美國(guó)地質(zhì)勘探局網(wǎng)站,運(yùn)用ENVI 5.1對(duì)影像數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正等預(yù)處理,通過(guò)ArcGIS10.3中人機(jī)交互的方式解譯了祁連山自然保護(hù)區(qū)2000年、2005年、2010年和2015年4期土地覆被數(shù)據(jù)。經(jīng)驗(yàn)證,kappa系數(shù)達(dá)到86.93%,解譯精度較高。

    草原類型數(shù)據(jù):來(lái)源于中國(guó)科學(xué)院1∶100萬(wàn)全國(guó)植被類型圖,基于1∶100萬(wàn)中國(guó)植被圖數(shù)據(jù)將甘肅祁連山草地植被分為溫帶叢生禾草草原、嵩草雜類草高寒草甸、禾草苔草高寒草原和溫帶叢生矮禾草矮半灌木荒漠草原四類,文中依次簡(jiǎn)稱典型草原、高寒草甸、高寒草原和荒漠草原。

    氣候數(shù)據(jù):來(lái)源于中國(guó)氣象數(shù)據(jù)中心逐日氣溫與降水量數(shù)據(jù)(圖1),通過(guò)ANUSPLIN軟件對(duì)其進(jìn)行空間插值。ANUSPLIN是基于局部薄盤樣條光滑函數(shù)開發(fā)的較為成熟的氣象要素插值軟件,可以同時(shí)批量處理多個(gè)空間表面,且在長(zhǎng)時(shí)間序列氣象數(shù)據(jù)的處理中有很大優(yōu)勢(shì),將最終插值結(jié)果同植被NDVI做響應(yīng)分析。

    1.2.2 趨勢(shì)分析法

    采用一元線性方程分析研究區(qū)生長(zhǎng)季NDVI對(duì)應(yīng)像元的時(shí)間變化趨勢(shì),從而反映祁連山自然保護(hù)區(qū)整體空間變化規(guī)律及不同草原區(qū)域格局的變化,其公式表示如下:

    圖1 祁連山自然保護(hù)區(qū)地理位置及數(shù)字高程模型Figure 1 Location and digital elevation model of the Qilian Mountain National Nature Reserve

    表1 祁連山自然保護(hù)區(qū)2000、2005、2010、2015年Landsat影像信息Table 1 Landsat image information for the Qilian Mountain National Nature Reserve for 2000, 2005, 2010, and 2015

    式中:SLOPE代表斜率,其值大于0,說(shuō)明變化趨勢(shì)呈現(xiàn)改善的趨勢(shì),相反則是退化的趨勢(shì);n為年數(shù);i為年序號(hào),NDVIi表示第i年的NDVI最大。參照前人的研究[29],將SLOPE劃分為7個(gè)等級(jí),分別為嚴(yán)重退化(≤-0.009 1)、中度退化(-0.009 0~-0.004 6)、輕微退化(-0.004 5~-0.001 0)、基本不變(-0.000 9~0.000 9)、輕微改善(0.001 0~0.004 5)、中度改善(0.004 6~0.009 0)、明顯改善(≥ 0.009 1)。

    1.2.3 氣象插值方法

    ANUSPLIN是基于薄盤樣條理論的專業(yè)插值軟件,包括普通薄盤和局部薄盤樣條函數(shù)。由于局部薄盤光滑樣條法除普通的樣條自變量外,還允許引入線性協(xié)變量子模型,諸如溫度和海拔之間的相關(guān)關(guān)系[30],因此本研究采用該方法對(duì)所需氣象站點(diǎn)數(shù)據(jù)(氣溫與降水量)進(jìn)行插值。其原理表述如下:

    式中:Zi為i點(diǎn)的因變量;χi為d維插值自變量向量,f(χi)是χi的 未知光滑函數(shù);yi是p維的獨(dú)立協(xié)變量矩陣,b為yi的 未知p維系數(shù)矩陣;ei是期望為0、方差為wiσ2的獨(dú)立變量,wi是i的相對(duì)誤差方差,σ2為常數(shù),是所有數(shù)據(jù)點(diǎn)上的誤差方差。需注意:當(dāng)p=0時(shí),即缺少bTyi,模型可簡(jiǎn)化為普通薄盤光滑樣條,而若f(χi)為0時(shí),即沒(méi)有獨(dú)立樣條變量,ANUSPLIN不可用。本研究以經(jīng)度、緯度作為獨(dú)立變量,高程為獨(dú)立協(xié)變量,對(duì)研究區(qū)氣象點(diǎn)數(shù)據(jù)進(jìn)行空間插值。

    1.2.4 相關(guān)關(guān)系分析

    本研究從柵格尺度上對(duì)不同草原植被指數(shù)NDVI與氣候因素包括年降水量、年平均氣溫進(jìn)行相關(guān)性運(yùn)算,從而在空間上表達(dá)NDVI與氣候因素的相關(guān)程度的分布特征。計(jì)算兩個(gè)時(shí)間序列xt和yt的相關(guān)系數(shù):

    式中:σx,σy表示x和y的均方差;和y的均值;Cxy為x和y的協(xié)方差;rxy表示x和y的相關(guān)系數(shù);n為時(shí)間序列長(zhǎng)度。

    相關(guān)性系數(shù)表達(dá)了兩個(gè)變量之間的相關(guān)性程度,其值的范圍[-1,1],絕對(duì)值越大,說(shuō)明相關(guān)性越好。本研究的相關(guān)性研究劃分標(biāo)準(zhǔn)如表2所列。

    表2 相關(guān)系數(shù)與相關(guān)程度分類Table 2 R values and correlation classification

    2 結(jié)果與分析

    2.1 祁連山自然保護(hù)區(qū)2000-2015年草地面積時(shí)空動(dòng)態(tài)變化分析

    分析2000年、2005年、2010年和2015年祁連山自然保護(hù)區(qū)草地的面積變化發(fā)現(xiàn)(表3),草地是祁連山自然保護(hù)區(qū)內(nèi)重要的土地覆被類型,其面積占比在47%左右上下波動(dòng)??傮w來(lái)說(shuō),祁連山自然保護(hù)區(qū)內(nèi)草地的面積變化主要表現(xiàn)在前期(2000-2005年)減少,后期2010年和2015年又相對(duì)增加;2000-2015年其面積變化總體上呈現(xiàn)為減少,減少了1.98 km2。

    表3 祁連山自然保護(hù)區(qū)2000-2015年草地面積變化Table 3 Grassland area change in the Qilian Mountain National Nature Reserve from 2000 to 2015

    在退耕還林還草等一系列生態(tài)保護(hù)與修復(fù)工程實(shí)施的背景下,草地在生長(zhǎng)過(guò)程中受到地表水熱組合和人類活動(dòng)的雙重影響,所以祁連山自然保護(hù)區(qū)內(nèi)草地面積的變化與氣候變化和人類活動(dòng)關(guān)系緊密。進(jìn)一步分析2000年與2015年草地與其他土地覆被類型間的時(shí)空轉(zhuǎn)移及變化特征(表4、圖2),結(jié)果表明:草地地類轉(zhuǎn)移最明顯的圖譜單元為“32”和“23”,即草地轉(zhuǎn)為林地和林地轉(zhuǎn)為草地,轉(zhuǎn)移面積分別為799.01和786.96 km2,轉(zhuǎn)移區(qū)域主要分布在祁連山自然保護(hù)區(qū)的中部與東南部;其次為“63”和“36”,即未利用地轉(zhuǎn)為草地和草地轉(zhuǎn)為未利用地,轉(zhuǎn)移面積分別為500.3和469.69 km2,轉(zhuǎn)移區(qū)域主要集中在祁連山自然保護(hù)區(qū)的北部與西南邊緣;再次為“31”和“13”,即草地轉(zhuǎn)為耕地和耕地轉(zhuǎn)為草地,轉(zhuǎn)移面積分別為94.86和88.17 km2,主要分布于自然保護(hù)區(qū)的中部邊緣和東南部邊緣區(qū)域;再次為“34”和“43”,即草地轉(zhuǎn)為水域和水域轉(zhuǎn)為草地,轉(zhuǎn)移面積分別為34.86和22.71 km2,其空間分布差異顯著,主要集中分布在保護(hù)區(qū)的北部、中部、南部;最后為“35”和“53”,轉(zhuǎn)移面積占比相對(duì)較小,草地轉(zhuǎn)出為農(nóng)村居民用地的面積為5.57 km2,城鄉(xiāng)工礦居民用地轉(zhuǎn)入草地的面積為4.12 km2,零星分布在研究區(qū)的東南邊緣。

    表4 祁連山自然保護(hù)區(qū)2000和2015年草地與其他土地覆被類型間的轉(zhuǎn)移矩陣Table 4 Transitions between grassland and other land-use types in the Qilian Mountain National Nature Reserve between 2000 and 2015 km2

    圖2 草地地類轉(zhuǎn)移圖譜單元空間差異(a)與面積變化(b)Figure 2 Grassland class transition map (a) and area change (b)

    綜上可得,草地轉(zhuǎn)為林地明顯,這說(shuō)明區(qū)域內(nèi)退耕還林還草等一系列生態(tài)保護(hù)與修復(fù)工程的實(shí)施以及生態(tài)經(jīng)濟(jì)林的建設(shè)使得祁連山自然保護(hù)區(qū)生態(tài)功能區(qū)的生態(tài)功能進(jìn)一步加強(qiáng),有助于區(qū)域水源涵養(yǎng)功能的發(fā)揮。草地轉(zhuǎn)為未利用地和城鄉(xiāng)工礦居民用地的面積大于轉(zhuǎn)入面積,數(shù)據(jù)量?jī)H為保護(hù)區(qū)草地面積的0.4%,隨著近年生態(tài)環(huán)境整改的進(jìn)一步落實(shí),數(shù)據(jù)會(huì)進(jìn)一步減少。

    2.2 不同草原類型的遙感監(jiān)測(cè)

    2.2.1 季節(jié)尺度下不同草原類型的時(shí)空變化趨勢(shì)

    計(jì)算分析2000-2018年典型草原、高寒草甸、高寒草原和荒漠草原4類草地春季、夏季、秋季和冬季NDVI的變化趨勢(shì)(圖3-圖6)。

    2000-2018年春季,4種類型草地長(zhǎng)勢(shì)變化差異不明顯,以無(wú)明顯變化為主(圖3、表5)。整個(gè)自然保護(hù)區(qū)輕微改善的草地類型中,典型草原的面積最大,占保護(hù)區(qū)草地總面積的18.54%;其次是高寒草甸,面積占保護(hù)區(qū)草地總面積的7.58%。此外,高寒草甸的退化面積占自然保護(hù)區(qū)草地總面積的8.64%,在4個(gè)草地類型中退化最為顯著,主要分布在自然保護(hù)區(qū)的中部和東部。

    2000-2018年夏季,4種類型草地均以輕微改善和中度改善為主(圖4、表6)。相比其他草地類型,高寒草甸輕微改善的面積最大,占到保護(hù)區(qū)草地總面積的20.06%。中度改善的草地中,典型草原的面積最大,占保護(hù)區(qū)草地總面積的23.44%,其次是高寒草甸,面積占保護(hù)區(qū)草地總面積的13.52%,高寒草原中度改善的面積占保護(hù)區(qū)草地總面積的8.13%,中度改善區(qū)域重要分布在保護(hù)區(qū)的中部和西北部。

    2000-2018年秋季,4種類型草地均以輕微改善為主(圖5、表7)。其中,典型草原輕微改善的面積最大,占保護(hù)區(qū)草地總面積的31.51%;其次是高寒草甸,面積占比為28.07%。

    圖3 2000-2018年祁連山自然保護(hù)區(qū)春季不同草原類型植被NDVI的空間變化趨勢(shì)Figure 3 Spatial variation trends in the NDVI for different grassland types in the Qilian Mountain National Nature Reserve in spring of 2000-2018

    圖4 2000-2018年祁連山自然保護(hù)區(qū)夏季不同草原類型植被NDVI的空間變化趨勢(shì)Figure 4 Spatial variation trends in the NDVI of different grassland types in the Qilian Mountain National Nature Reserve in summer of 2000-2018

    2000-2018年冬季,4種類型草地同樣以輕微改善為主(圖6、表8)。相比其他草地類型,典型草原輕度改善的面積最大,占保護(hù)區(qū)草地總面積的27.28%;其次是高寒草甸,面積占比為24.76%。結(jié)合空間分布來(lái)看,高寒草甸在張掖市東南部靠近青海省門源縣存在少量草地輕微退化,典型草原在張掖市東南部靠近金昌市永昌縣同樣分布了一些退化草地,其他區(qū)域草地普遍輕微改善。

    綜上,2000-2018年,祁連山自然保護(hù)區(qū)不同類型草地的長(zhǎng)勢(shì)變化在季節(jié)尺度上均存在較明顯的差異,同一季節(jié),4種草地類型草地的主要變化趨勢(shì)相似。整個(gè)保護(hù)區(qū)除了春季,其他季節(jié)草地均以改善為主。草地在春季的退化面積要明顯大于其他季節(jié),草地在夏季的改善程度要明顯大于其他季節(jié),秋季改善程度大于冬季。

    2.2.2 不同草原類型年際時(shí)空變化趨勢(shì)

    2000-2018年,祁連山自然保護(hù)區(qū)典型草原、高寒草甸、高寒草原、荒漠草原的NDVI總體呈現(xiàn)明顯的增加趨勢(shì),增加速率分別為0.005 3·(19a)-1、0.003·(19a)-1、0.005 3·(19a)-1和 0.006 2·(19a)-1(圖7)。各個(gè)草原類型的NDVI值均在2001年達(dá)到最小,而2018年達(dá)到最大值,其中,典型草原、高寒草原、高寒草甸和荒漠草原的NDVI值范圍分別為0.38~0.53、0.35~0.50、0.39~0.47 和 0.26~0.41,可見高寒草甸的NDVI值比較穩(wěn)定,而其他3種類型的NDVI值波動(dòng)較大。

    2.2.3 不同草原類型變化對(duì)氣候的時(shí)空響應(yīng)特征

    本研究從柵格尺度上對(duì)不同草原植被指數(shù)NDVI與氣候因素包括年降水量、年平均氣溫進(jìn)行相關(guān)性運(yùn)算,得到NDVI與氣候因素相關(guān)程度的空間分布特征及對(duì)應(yīng)的面積統(tǒng)計(jì)結(jié)果(圖8、表9)。

    圖5 2000-2018年祁連山自然保護(hù)區(qū)秋季不同草原類型植被NDVI的空間變化趨勢(shì)Figure 5 Spatial variation trends in the NDVI for different grassland types in the Qilian Mountain National Nature Reserve in autumn of 2000-2018

    從典型草原植被NDVI與年降水的相關(guān)性空間分布來(lái)看,整體上NDVI與年降水量相關(guān)性較氣溫高。NDVI與年降水的相關(guān)性表達(dá)以輕度和中度相關(guān)為主,相關(guān)性表達(dá)為顯著相關(guān)的區(qū)域主要分布在自然保護(hù)區(qū)的西部,占整個(gè)區(qū)域總面積的23.06%;NDVI與氣溫的相關(guān)性表達(dá)以輕度相關(guān)為主,相關(guān)性表達(dá)為中度相關(guān)、顯著相關(guān)的區(qū)域主要分布在張掖市東部與金昌市接壤的區(qū)域,其中,顯著相關(guān)的區(qū)域僅占整個(gè)區(qū)域總面積的3.62%。從高寒草甸植被NDVI與年降水的相關(guān)性空間分布來(lái)看,整體上NDVI與年降水的相關(guān)性表達(dá)以輕度相關(guān)為主,顯著相關(guān)的區(qū)域主要集中在保護(hù)區(qū)的西南部邊緣地區(qū),占區(qū)域總面積的17.02%;NDVI與年平均氣溫的相關(guān)性較低,整體上以輕度相關(guān)為主,中度相關(guān)的區(qū)域在自然保護(hù)區(qū)中部零星分布,顯著相關(guān)區(qū)域僅占區(qū)域總面積的3.59%。從高寒草原植被NDVI與年降水的相關(guān)性空間分布來(lái)看,整體上NDVI與年降水的輕度、中度和顯著相關(guān)的面積差異較小,顯著相關(guān)的面積占區(qū)域總面積的35.33%,在4類草地類型中面積占比最大。相比年降水,年平均氣溫與NDVI的相關(guān)性較低,整體上以輕度相關(guān)為主,輕度相關(guān)的比重占區(qū)域總面積的82.02%?;哪菰脖籒DVI與年降水的相關(guān)性表達(dá)以中度相關(guān)為主,顯著相關(guān)占區(qū)域總面積的20.33%。與其他3類草地類型一樣,年平均氣溫與NDVI的相關(guān)性水平要低于與年降水的相關(guān)性水平,顯著相關(guān)的區(qū)域面積占10.16%。

    總體而言,甘肅祁連山自然保護(hù)區(qū)內(nèi)不同草原類型植被NDVI與年降水、年平均氣溫的相關(guān)性時(shí)空響應(yīng)特征具有一定的差異,但整體來(lái)看,草地植被NDVI與年降水的相關(guān)性水平均高于年平均氣溫,與年平均氣溫的相關(guān)性水平主要表達(dá)為輕度相關(guān)為主。說(shuō)明區(qū)域年降水比年平均氣溫對(duì)草原植被NDVI的影響程度更大。

    分析祁連山自然保護(hù)區(qū)不同草原類型季節(jié)NDVI與當(dāng)季、上一季氣溫相關(guān)關(guān)系表明(表10),4類草原類型的季節(jié)NDVI與當(dāng)季氣溫的相關(guān)性均不顯著。其中,高寒草甸植被NDVI與氣溫的相關(guān)性最大,相關(guān)性系數(shù)為0.311,其余草原類型與植被NDVI的相關(guān)系數(shù)均小于0.3。不同草原類型植被NDVI與氣溫的相關(guān)性除高寒草甸外其余3個(gè)類型草原以夏季最高,而秋季最差。4類草原植被NDVI與上季度氣溫的相關(guān)性大多呈不顯著特點(diǎn),秋季各草原類型NDVI與氣溫的相關(guān)性最高,滯后性最為明顯。其中,僅有典型草原秋季NDVI與夏季氣溫相關(guān)關(guān)系達(dá)顯著水平。

    圖6 2000-2018年祁連山自然保護(hù)區(qū)冬季不同草原類型植被NDVI的空間變化趨勢(shì)Figure 6 Spatial variation trends in the NDVI for different grassland types in the Qilian Mountain National Nature Reserve in winter of 2000-2018

    表5 祁連山自然保護(hù)區(qū)不同草原類型春季變化趨勢(shì)的面積統(tǒng)計(jì)Table 5 Spatial trends for different grassland types in the Qilian Mountain National Nature Reserve in spring 2000-2018

    表6 祁連山自然保護(hù)區(qū)不同草原類型夏季變化趨勢(shì)的面積統(tǒng)計(jì)Table 6 Spatial trends of different grassland types in the Qilian Mountain National Nature Reserve in summer of 2000-2018

    表7 祁連山自然保護(hù)區(qū)不同草原類型秋季變化趨勢(shì)的面積統(tǒng)計(jì)Table 7 Spatial trends for different grassland types in the Qilian MountainNational Nature Reserve in autumn of 2000-2018

    表8 祁連山自然保護(hù)區(qū)不同草原類型冬季變化趨勢(shì)的面積統(tǒng)計(jì)Table 8 Spatial trends for different grassland types in the Qilian Mountain National Nature Reserve in winter of 2000-2018

    分析祁連山自然保護(hù)區(qū)不同草原類型季節(jié)NDVI與當(dāng)季、上一季降水的相關(guān)關(guān)系表明(表11),不同草原類型植被NDVI與降水量的相關(guān)性均高于氣溫,各草原類型與降水量的相關(guān)性大于0.3,其中,荒漠草原植被NDVI與降水量的相關(guān)系數(shù)為0.612,達(dá)到0.01顯著水平。在4類草原類型中,夏季與秋季的植被NDVI與降水量的相關(guān)性明顯高于春季與冬季。除夏季外,4類草原類型植被NDVI與上季降水量的相關(guān)系數(shù)基本上大于0.3,其中,典型草原春季NDVI與冬季降水量的相關(guān)系數(shù)為0.553(P<0.05);荒漠草原春季NDVI與冬季降水量的相關(guān)系數(shù)為0.626 (P<0.01);荒漠草原秋季NDVI與夏季降水量顯著相關(guān),荒漠草原冬季NDVI與冬季降水量顯著負(fù)相關(guān)。春季不同草原植被NDVI與上季度降水量的滯后性最為明顯,其次是冬季,而夏季滯后性最弱。

    圖7 2000-2018年祁連山自然保護(hù)區(qū)不同草原類型NDVI的變化趨勢(shì)Figure 7 Variation in the NDVI for different grassland types in the Qilian Mountain National Nature Reserve from 2000 to 2018

    3 討論與結(jié)論

    本研究在分析2000-2015年祁連山自然保護(hù)區(qū)草地面積變化的基礎(chǔ)上,基于MOD13Q1 NDVI產(chǎn)品,結(jié)合草原類型,從季節(jié)、年際兩個(gè)尺度分別分析了祁連山自然保護(hù)區(qū)草地植被NDVI的時(shí)空變化趨勢(shì),并討論了不同草原類型NDVI對(duì)氣候變化響應(yīng)特征。結(jié)果表明:

    1) 2000-2015年草地面積變化中,草地作為祁連山自然保護(hù)區(qū)內(nèi)重要的土地覆被類型,其變化主要體現(xiàn)為2000-2005年減少,于2010年和2015年又相對(duì)增加。草地轉(zhuǎn)出以轉(zhuǎn)為林地最明顯,轉(zhuǎn)出為耕地和水域較為明顯,草地轉(zhuǎn)出為林地、耕地和水域的面積分別小于各自轉(zhuǎn)變?yōu)椴莸氐拿娣e;草地轉(zhuǎn)出為城鄉(xiāng)工礦未利用地和未利用地的面積大于其轉(zhuǎn)入為草地的面積。

    2) 2000-2018年祁連山自然保護(hù)區(qū)除春季外,其他季節(jié)均以輕微改善為主,夏季植被NDVI的改善狀況最佳。其中,春季高寒草甸的退化面積占自然保護(hù)區(qū)草地總面積的8.64%,而典型草原輕度改善草地面積占草地總面積的18.54%;夏季輕微改善的草地中高寒草甸的面積最大,占到保護(hù)區(qū)草地總面積的20.06%,中度改善的草地中典型草原的面積最大,占保護(hù)區(qū)草地總面積的23.44%;秋季保護(hù)區(qū)草地改善程度整體優(yōu)于冬季。

    3)祁連山自然保護(hù)區(qū)不同草原植被NDVI均呈現(xiàn)增加的趨勢(shì),且以荒漠草原的增加速率最大0.006 2·(19a)-1,其次是典型草原與高寒草原,均為0.005 3·(19a)-1,而高寒草甸的增加速率最小,僅為 0.003·(19a)-1。戴聲佩等[31]基于 SPOT VGT-NDVI數(shù)據(jù)對(duì)祁連山不同草地植被NDVI變化趨勢(shì)進(jìn)行研究,發(fā)現(xiàn)典型草原和平原草地NDVI增加速率高于高寒草甸,雖本研究采用的數(shù)據(jù)和方法與其不盡相同,但得到的結(jié)果基本一致。

    圖8 2000-2018年不同草原類型NDVI與氣溫(左)、降水量(右)的相關(guān)性Figure 8 Correlations between NDVI and temperature (left panels), and precipitation (right)for different grassland types from 2000 to 2018

    表9 不同草原類型NDVI與氣溫、降水量的面積比例Table 9 Percentage area of different grassland types in the NDVI with temperature and precipitation under different degrees of correlation %

    表10 不同草原類型季節(jié)NDVI與氣溫的相關(guān)系數(shù)Table 10 Correlation coefficients between NDVI and temperature for different grassland types

    表11 不同草原類型季節(jié)NDVI與降水量的相關(guān)系數(shù)Table 11 Correlation coefficients between NDVI and precipitation for different grassland types

    4)祁連山自然保護(hù)區(qū)不同草原類型植被NDVI與年降水、年平均氣溫的相關(guān)性呈現(xiàn)出一定的差異,但整體上與年降水的相關(guān)性水平均要高于年平均氣溫,與年平均氣溫的相關(guān)性水平主要表達(dá)為輕度相關(guān)且鮮有顯著相關(guān)。說(shuō)明年降水比年平均氣溫對(duì)草原植被NDVI的影響程度更大。張禹舜[25]研究祁連山凈初級(jí)生產(chǎn)力NPP變化對(duì)氣候的響應(yīng)中也同樣指出草甸、草原和荒漠植被對(duì)降水的響應(yīng)程度要明顯大于氣溫。

    5)祁連山自然保護(hù)區(qū)降水量變化對(duì)不同草原植被NDVI的影響高于氣溫,且與草原植被NDVI的滯后性較為明顯,其中以春季的滯后性最強(qiáng),其次是冬季,而夏季滯后性最弱。黃德青等[32]的研究表明祁連山自然保護(hù)區(qū)生物量與降水量正相關(guān),而與氣溫未達(dá)到顯著水平,這與本研究降水量是影響不同草原植被NDVI的主要因素的結(jié)論吻合。

    猜你喜歡
    祁連山草甸自然保護(hù)區(qū)
    高山草甸
    圖志
    發(fā)展(2024年1期)2024-04-25 03:45:30
    The world's narrowest river
    山地草甸
    七姊妹山國(guó)家級(jí)自然保護(hù)區(qū)
    祁連山下
    武功山山地草甸的成因調(diào)查
    活力(2019年21期)2019-04-01 12:17:10
    保護(hù)自然保護(hù)區(qū)刻不容緩
    祁連山草原:如夢(mèng)如幻近高天
    淺淡黃崗梁自然保護(hù)區(qū)的資源保護(hù)與開發(fā)利用
    精品久久久精品久久久| 王馨瑶露胸无遮挡在线观看| 国产三级黄色录像| 亚洲全国av大片| 久久国产精品影院| 日日夜夜操网爽| 日本一区二区免费在线视频| 操出白浆在线播放| 久久精品国产99精品国产亚洲性色 | 精品久久久精品久久久| 色综合婷婷激情| 热re99久久国产66热| 熟女少妇亚洲综合色aaa.| xxxhd国产人妻xxx| 欧美日韩精品网址| 日韩免费高清中文字幕av| 一级毛片女人18水好多| 制服人妻中文乱码| 又大又爽又粗| 亚洲五月婷婷丁香| 久久这里只有精品19| √禁漫天堂资源中文www| 欧美另类亚洲清纯唯美| 真人做人爱边吃奶动态| 欧美精品啪啪一区二区三区| 国产视频一区二区在线看| 9191精品国产免费久久| 亚洲一区中文字幕在线| 国产精品久久久久久精品电影小说| 在线观看www视频免费| 亚洲美女黄片视频| 1024香蕉在线观看| 日韩三级视频一区二区三区| 国产亚洲欧美在线一区二区| 18禁美女被吸乳视频| 亚洲av第一区精品v没综合| 黄片大片在线免费观看| 欧美人与性动交α欧美精品济南到| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美在线一区| 国产片内射在线| 脱女人内裤的视频| 在线十欧美十亚洲十日本专区| 亚洲三区欧美一区| 精品亚洲成国产av| 亚洲,欧美精品.| 久久人妻福利社区极品人妻图片| 免费人妻精品一区二区三区视频| bbb黄色大片| 视频区图区小说| videosex国产| 久久精品亚洲av国产电影网| 亚洲人成伊人成综合网2020| 黄色a级毛片大全视频| 色婷婷久久久亚洲欧美| 不卡一级毛片| 精品国产一区二区三区四区第35| av有码第一页| cao死你这个sao货| 一区二区av电影网| 久久国产精品人妻蜜桃| 免费日韩欧美在线观看| 韩国精品一区二区三区| 欧美黑人欧美精品刺激| 男女之事视频高清在线观看| 99香蕉大伊视频| 亚洲综合色网址| 色播在线永久视频| 亚洲专区国产一区二区| 日韩欧美国产一区二区入口| 久久久欧美国产精品| 最近最新中文字幕大全免费视频| 在线观看免费日韩欧美大片| 蜜桃在线观看..| 国产精品一区二区在线观看99| 色精品久久人妻99蜜桃| 老司机午夜福利在线观看视频 | 欧美人与性动交α欧美软件| 国产一区二区三区视频了| 午夜免费成人在线视频| 日韩欧美三级三区| 欧美激情高清一区二区三区| 五月开心婷婷网| 国产在线免费精品| 法律面前人人平等表现在哪些方面| 国产精品成人在线| 交换朋友夫妻互换小说| 深夜精品福利| 国产真人三级小视频在线观看| 久久久久久人人人人人| 欧美精品亚洲一区二区| 女人爽到高潮嗷嗷叫在线视频| 男男h啪啪无遮挡| 一级毛片电影观看| 国产亚洲午夜精品一区二区久久| 大片电影免费在线观看免费| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 九色亚洲精品在线播放| 巨乳人妻的诱惑在线观看| 亚洲av美国av| 久久久久久久久久久久大奶| 免费不卡黄色视频| 久久ye,这里只有精品| 一本久久精品| www.自偷自拍.com| 日本黄色日本黄色录像| 久久久久久久精品吃奶| 亚洲国产av新网站| 国产在线精品亚洲第一网站| 无限看片的www在线观看| 免费在线观看影片大全网站| 蜜桃在线观看..| 日韩人妻精品一区2区三区| 欧美成狂野欧美在线观看| 老司机午夜福利在线观看视频 | 色婷婷av一区二区三区视频| 欧美久久黑人一区二区| 少妇 在线观看| 丁香六月天网| 在线天堂中文资源库| 男女之事视频高清在线观看| 人人妻人人澡人人看| 中文字幕最新亚洲高清| 黄色视频在线播放观看不卡| 日韩精品免费视频一区二区三区| 18禁美女被吸乳视频| 精品一区二区三区四区五区乱码| 精品少妇久久久久久888优播| 国产精品亚洲一级av第二区| 欧美在线黄色| 男女床上黄色一级片免费看| 国产成人精品无人区| 中文字幕色久视频| 黄色毛片三级朝国网站| 亚洲欧美日韩另类电影网站| 老熟女久久久| 日韩 欧美 亚洲 中文字幕| 人妻 亚洲 视频| 国产精品av久久久久免费| 中文字幕另类日韩欧美亚洲嫩草| 国产一卡二卡三卡精品| 色在线成人网| av线在线观看网站| www.自偷自拍.com| 一进一出好大好爽视频| 叶爱在线成人免费视频播放| 国产在线免费精品| 欧美精品高潮呻吟av久久| 黄色视频在线播放观看不卡| 菩萨蛮人人尽说江南好唐韦庄| 蜜桃在线观看..| 999久久久精品免费观看国产| 国产熟女午夜一区二区三区| 五月开心婷婷网| 自线自在国产av| 窝窝影院91人妻| 交换朋友夫妻互换小说| 大型黄色视频在线免费观看| 精品一区二区三区av网在线观看 | 乱人伦中国视频| 久久人妻av系列| 91大片在线观看| 夫妻午夜视频| 国产成人免费观看mmmm| 国产精品成人在线| 一区二区三区精品91| 国产精品麻豆人妻色哟哟久久| 一边摸一边做爽爽视频免费| 日日摸夜夜添夜夜添小说| 高潮久久久久久久久久久不卡| 亚洲欧洲精品一区二区精品久久久| 丰满迷人的少妇在线观看| av福利片在线| 黄片大片在线免费观看| 亚洲一码二码三码区别大吗| 视频区欧美日本亚洲| 91字幕亚洲| 波多野结衣av一区二区av| 少妇精品久久久久久久| 老司机午夜十八禁免费视频| 啪啪无遮挡十八禁网站| 中国美女看黄片| 老鸭窝网址在线观看| 国产免费福利视频在线观看| 国产一区二区三区综合在线观看| 欧美老熟妇乱子伦牲交| 成人免费观看视频高清| 在线天堂中文资源库| 国产一区二区激情短视频| 午夜福利欧美成人| 一个人免费看片子| 国产高清视频在线播放一区| 亚洲中文字幕日韩| 国产精品1区2区在线观看. | 啦啦啦 在线观看视频| 国产精品久久久av美女十八| 嫩草影视91久久| 丝袜在线中文字幕| 丝袜在线中文字幕| 国产成人欧美| 国产高清视频在线播放一区| 亚洲美女黄片视频| 精品亚洲乱码少妇综合久久| 亚洲全国av大片| 老熟女久久久| 女人高潮潮喷娇喘18禁视频| 丝袜在线中文字幕| 国产精品影院久久| 波多野结衣av一区二区av| 1024香蕉在线观看| 99精品欧美一区二区三区四区| 精品人妻熟女毛片av久久网站| 久久精品人人爽人人爽视色| 精品少妇一区二区三区视频日本电影| 精品国产一区二区三区久久久樱花| 91成人精品电影| 免费在线观看视频国产中文字幕亚洲| av片东京热男人的天堂| 色精品久久人妻99蜜桃| 国产欧美日韩精品亚洲av| 亚洲少妇的诱惑av| 女人爽到高潮嗷嗷叫在线视频| 久久久久精品国产欧美久久久| 99国产精品99久久久久| 99国产精品一区二区三区| 久久久精品区二区三区| 韩国精品一区二区三区| av欧美777| 18禁国产床啪视频网站| 欧美黑人精品巨大| kizo精华| 露出奶头的视频| 精品熟女少妇八av免费久了| 老司机午夜福利在线观看视频 | 亚洲精品自拍成人| 男女边摸边吃奶| 午夜福利视频在线观看免费| 精品福利观看| 久久久精品区二区三区| 搡老岳熟女国产| 欧美激情极品国产一区二区三区| 久久久久久久精品吃奶| 国产成人精品无人区| 黄片大片在线免费观看| 精品国产一区二区久久| 亚洲国产av影院在线观看| 国产单亲对白刺激| 正在播放国产对白刺激| 中文字幕av电影在线播放| 亚洲自偷自拍图片 自拍| 一本大道久久a久久精品| 麻豆成人av在线观看| 嫩草影视91久久| av在线播放免费不卡| 日韩欧美三级三区| 在线天堂中文资源库| 麻豆av在线久日| 久9热在线精品视频| 一本—道久久a久久精品蜜桃钙片| 免费在线观看黄色视频的| 高清av免费在线| 黄色怎么调成土黄色| 亚洲 欧美一区二区三区| 制服人妻中文乱码| 一区二区三区乱码不卡18| 老司机亚洲免费影院| 亚洲全国av大片| 精品视频人人做人人爽| 国产成人系列免费观看| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av高清一级| 欧美日韩中文字幕国产精品一区二区三区 | 欧美老熟妇乱子伦牲交| 涩涩av久久男人的天堂| 蜜桃在线观看..| av超薄肉色丝袜交足视频| 国产精品秋霞免费鲁丝片| 自拍欧美九色日韩亚洲蝌蚪91| 黄色丝袜av网址大全| 午夜91福利影院| 日韩中文字幕欧美一区二区| 精品国产一区二区三区四区第35| 亚洲成人国产一区在线观看| 性色av乱码一区二区三区2| 男女边摸边吃奶| 精品欧美一区二区三区在线| 国产在线精品亚洲第一网站| 超碰成人久久| 国产视频一区二区在线看| 亚洲av成人不卡在线观看播放网| 欧美激情 高清一区二区三区| 少妇裸体淫交视频免费看高清 | 在线亚洲精品国产二区图片欧美| 日韩成人在线观看一区二区三区| 久久精品aⅴ一区二区三区四区| 久久久水蜜桃国产精品网| 久久精品成人免费网站| 一进一出抽搐动态| 夫妻午夜视频| 91麻豆av在线| 黄片播放在线免费| 欧美日韩国产mv在线观看视频| 亚洲五月色婷婷综合| 女人被躁到高潮嗷嗷叫费观| 啦啦啦免费观看视频1| 亚洲七黄色美女视频| 99香蕉大伊视频| 欧美久久黑人一区二区| 一本大道久久a久久精品| 欧美日韩精品网址| 女人高潮潮喷娇喘18禁视频| a级片在线免费高清观看视频| www.999成人在线观看| 肉色欧美久久久久久久蜜桃| 国产麻豆69| 国产精品久久电影中文字幕 | 国产精品一区二区在线观看99| 男女边摸边吃奶| 久久久国产成人免费| 99久久国产精品久久久| 国产麻豆69| 久久久久久久久免费视频了| 日韩一区二区三区影片| avwww免费| 19禁男女啪啪无遮挡网站| 夜夜爽天天搞| 在线av久久热| 9色porny在线观看| 午夜91福利影院| 12—13女人毛片做爰片一| 在线观看免费视频网站a站| 黄频高清免费视频| 三级毛片av免费| 啦啦啦在线免费观看视频4| 免费日韩欧美在线观看| 国产亚洲精品第一综合不卡| 亚洲第一青青草原| 少妇裸体淫交视频免费看高清 | 97在线人人人人妻| 深夜精品福利| 90打野战视频偷拍视频| 国产黄频视频在线观看| 777米奇影视久久| 亚洲情色 制服丝袜| 成人永久免费在线观看视频 | 女人久久www免费人成看片| 桃红色精品国产亚洲av| av又黄又爽大尺度在线免费看| 女性被躁到高潮视频| 纯流量卡能插随身wifi吗| 亚洲熟妇熟女久久| 69精品国产乱码久久久| 午夜久久久在线观看| 99九九在线精品视频| 国产精品亚洲一级av第二区| 国产不卡一卡二| cao死你这个sao货| 欧美日韩亚洲国产一区二区在线观看 | 久久久久网色| 成人特级黄色片久久久久久久 | 在线观看一区二区三区激情| 欧美黑人欧美精品刺激| 别揉我奶头~嗯~啊~动态视频| 国产精品电影一区二区三区 | 欧美乱码精品一区二区三区| 亚洲精品国产区一区二| 亚洲男人天堂网一区| 变态另类成人亚洲欧美熟女 | 制服人妻中文乱码| 水蜜桃什么品种好| 久热这里只有精品99| 欧美精品人与动牲交sv欧美| 色综合欧美亚洲国产小说| 老司机亚洲免费影院| 久久久久久久大尺度免费视频| 久久毛片免费看一区二区三区| 亚洲一码二码三码区别大吗| 一区二区av电影网| 精品高清国产在线一区| 久久精品国产a三级三级三级| 国产成人精品在线电影| 国产精品美女特级片免费视频播放器 | bbb黄色大片| 淫妇啪啪啪对白视频| 在线观看免费视频日本深夜| 国产在线视频一区二区| 亚洲成a人片在线一区二区| 亚洲av片天天在线观看| 无遮挡黄片免费观看| 人妻一区二区av| 国产成人影院久久av| 精品一区二区三区四区五区乱码| 啦啦啦在线免费观看视频4| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 国产黄色免费在线视频| 黄色a级毛片大全视频| 日韩免费av在线播放| 日本a在线网址| 多毛熟女@视频| 精品一区二区三区视频在线观看免费 | 久久久久网色| 美女午夜性视频免费| 国产精品久久久久久人妻精品电影 | 一进一出抽搐动态| 国产精品一区二区精品视频观看| 国产淫语在线视频| 国产精品一区二区在线观看99| 日韩一卡2卡3卡4卡2021年| 多毛熟女@视频| 亚洲欧美精品综合一区二区三区| 久久久精品免费免费高清| 国产伦人伦偷精品视频| 一区二区av电影网| 少妇猛男粗大的猛烈进出视频| 国产日韩欧美亚洲二区| 51午夜福利影视在线观看| 窝窝影院91人妻| 午夜福利乱码中文字幕| 免费观看a级毛片全部| 久久国产精品影院| 满18在线观看网站| 91九色精品人成在线观看| 两个人免费观看高清视频| 亚洲人成伊人成综合网2020| 老熟妇乱子伦视频在线观看| 高清视频免费观看一区二区| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 亚洲熟妇熟女久久| 欧美另类亚洲清纯唯美| av免费在线观看网站| 他把我摸到了高潮在线观看 | xxxhd国产人妻xxx| 黑人巨大精品欧美一区二区蜜桃| 久久人妻熟女aⅴ| 一级毛片电影观看| 好男人电影高清在线观看| 一边摸一边抽搐一进一小说 | 久久九九热精品免费| 欧美另类亚洲清纯唯美| 激情在线观看视频在线高清 | 在线观看66精品国产| 少妇被粗大的猛进出69影院| 黄色视频,在线免费观看| 午夜久久久在线观看| 悠悠久久av| 国产成人精品久久二区二区91| 久久 成人 亚洲| 亚洲视频免费观看视频| 中文字幕人妻熟女乱码| 午夜福利视频精品| 伊人久久大香线蕉亚洲五| 午夜福利在线免费观看网站| 亚洲成a人片在线一区二区| 亚洲成人手机| 新久久久久国产一级毛片| 久久午夜亚洲精品久久| 香蕉久久夜色| 捣出白浆h1v1| 日韩中文字幕视频在线看片| 精品国产超薄肉色丝袜足j| 俄罗斯特黄特色一大片| 两人在一起打扑克的视频| 国产99久久九九免费精品| 色在线成人网| 亚洲一码二码三码区别大吗| 日韩大片免费观看网站| 中文字幕高清在线视频| 男女午夜视频在线观看| 亚洲av成人一区二区三| 在线观看免费视频网站a站| 精品国产乱码久久久久久小说| 欧美国产精品va在线观看不卡| 亚洲国产看品久久| 免费看十八禁软件| 脱女人内裤的视频| 午夜福利视频在线观看免费| 一级片免费观看大全| 黄色怎么调成土黄色| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人| 一级毛片电影观看| 亚洲美女黄片视频| 丝袜人妻中文字幕| 一本综合久久免费| 宅男免费午夜| 国产精品秋霞免费鲁丝片| 丝袜美足系列| 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 又黄又粗又硬又大视频| www.999成人在线观看| 国产精品亚洲av一区麻豆| 一区二区三区激情视频| 色在线成人网| 欧美日韩一级在线毛片| 亚洲欧美精品综合一区二区三区| 亚洲人成77777在线视频| 国产精品一区二区在线不卡| 黑人巨大精品欧美一区二区蜜桃| 久久免费观看电影| 国产激情久久老熟女| 少妇粗大呻吟视频| 亚洲欧美一区二区三区黑人| 久久久精品94久久精品| 亚洲伊人色综图| 亚洲国产欧美一区二区综合| 两个人看的免费小视频| 日韩视频在线欧美| 最新在线观看一区二区三区| 久久久久久久国产电影| 丝袜人妻中文字幕| 女性被躁到高潮视频| 国产欧美亚洲国产| 久久久久精品国产欧美久久久| 大型黄色视频在线免费观看| 日韩大码丰满熟妇| 无限看片的www在线观看| kizo精华| 国产单亲对白刺激| 久久久久久久大尺度免费视频| 电影成人av| 色综合婷婷激情| 捣出白浆h1v1| 亚洲伊人久久精品综合| 亚洲欧美激情在线| 久久久久久久精品吃奶| 中文字幕精品免费在线观看视频| kizo精华| 亚洲国产欧美在线一区| 日韩免费av在线播放| 欧美亚洲 丝袜 人妻 在线| av不卡在线播放| 一二三四社区在线视频社区8| av又黄又爽大尺度在线免费看| 美女视频免费永久观看网站| 满18在线观看网站| 窝窝影院91人妻| 99国产精品免费福利视频| 亚洲avbb在线观看| 一个人免费看片子| 欧美黄色淫秽网站| 午夜免费成人在线视频| 国产精品偷伦视频观看了| 777久久人妻少妇嫩草av网站| 国产av精品麻豆| 国产成人免费无遮挡视频| 久久av网站| 国产成人啪精品午夜网站| 国产男靠女视频免费网站| 国产熟女午夜一区二区三区| 中文亚洲av片在线观看爽 | 一区二区日韩欧美中文字幕| 国产单亲对白刺激| 精品人妻1区二区| 国产成人一区二区三区免费视频网站| 另类精品久久| 老鸭窝网址在线观看| 成年人午夜在线观看视频| 国产精品av久久久久免费| 亚洲av片天天在线观看| www.999成人在线观看| 午夜成年电影在线免费观看| 国产在线一区二区三区精| 别揉我奶头~嗯~啊~动态视频| 飞空精品影院首页| 成人黄色视频免费在线看| 欧美成人午夜精品| 久久精品亚洲av国产电影网| 免费在线观看影片大全网站| 在线观看免费日韩欧美大片| 欧美成人午夜精品| 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址 | 国产在线视频一区二区| 怎么达到女性高潮| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 另类亚洲欧美激情| 亚洲国产中文字幕在线视频| 日韩视频在线欧美| av片东京热男人的天堂| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 国产亚洲精品第一综合不卡| 精品午夜福利视频在线观看一区 | 一本一本久久a久久精品综合妖精| 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 黄色丝袜av网址大全| 热99久久久久精品小说推荐| 中文字幕最新亚洲高清| 久久中文看片网| 国产日韩一区二区三区精品不卡| 精品人妻在线不人妻| 91麻豆av在线| 午夜福利欧美成人| 国产精品 欧美亚洲| 国产高清激情床上av| 欧美中文综合在线视频| 国产一卡二卡三卡精品| 色视频在线一区二区三区| 国产精品二区激情视频| 免费在线观看日本一区| 国产欧美亚洲国产| 在线观看免费高清a一片| 色94色欧美一区二区| 欧美乱码精品一区二区三区| 成在线人永久免费视频| 欧美乱码精品一区二区三区| 欧美日韩亚洲综合一区二区三区_| 女人爽到高潮嗷嗷叫在线视频| 99精国产麻豆久久婷婷|