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

    歐亞大陸不同生態(tài)區(qū)植被生長對降水響應的季節(jié)變化規(guī)律

    2018-12-20 06:42:50孫立群李晴嵐
    生態(tài)學報 2018年22期
    關(guān)鍵詞:歐亞大陸生態(tài)區(qū)降水

    孫立群,李晴嵐,陳 驥,牛 俊

    1 中國科學院深圳先進技術(shù)研究院,深圳 518055 2 香港大學土木工程系,香港 9990773 中國農(nóng)業(yè)大學中國農(nóng)業(yè)水問題研究中心,北京 100083

    歐亞大陸土地廣袤,涵蓋了我國“一帶一路”戰(zhàn)略規(guī)劃的大部分國家。歐亞大陸人口眾多,不同國家和地區(qū)的地理條件、經(jīng)濟發(fā)展差異極大[1]。類似于中國的“胡煥庸線”難題,水資源對歐亞大陸很多國家,尤其是中亞干旱地區(qū)國家的生態(tài)安全和經(jīng)濟發(fā)展具有重要的意義[2]。了解歐亞大陸植被生長對降水響應的時空格局是“一帶一路”戰(zhàn)略順利實施的科學與生態(tài)基礎(chǔ)。在氣候變化的背景下,如何因地制宜地應對及適應氣候變化,從而促進當?shù)厣鐣c經(jīng)濟的可持續(xù)性發(fā)展,是歐亞大陸每一個國家和地區(qū)的人們所需要面對和亟待解決的問題[3]。

    植被生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的主體,在不同地理環(huán)境和氣象條件下對氣候變化有著不同響應機制。例如,在中國的內(nèi)蒙古地區(qū),植被生長受水分影響多一些[4]。而在中國東北地區(qū),則對溫度變化更為敏感[5]。植被生長是一個極為復雜的過程,在排除霜凍、病蟲害、水利灌溉和施肥等影響后,植被生長主要受到溫度、水分條件、日照輻射3個自然條件的脅迫[6]。選取合適的方法,分析陸地植被生長在不同地區(qū)對降水變化的響應關(guān)系是本文研究的重點。

    經(jīng)過數(shù)十年遙感數(shù)據(jù)的積累,在大陸尺度研究每個像元上植被生長與氣候變化的關(guān)系已經(jīng)成為可能。Peng等[7]使用歸一化植被指數(shù)NDVI(Normalized Difference Vegetation Index)和全球氣象格點數(shù)據(jù),采用偏相關(guān)方法(Partial Correlation Coefficient)單獨計算了日間和夜間增溫對北半球高緯度地區(qū)植被生態(tài)系統(tǒng)的影響,發(fā)現(xiàn)濕潤地區(qū)和干旱地區(qū)的植被對晝夜增溫的響應是不一致的。Shen等[8]使用SPOT(Satellite Pour l′Observation de la Terre)NDVI分析了青藏高原上降水對春季植被物候的影響,指出干旱地區(qū)的植被為了最大限度地利用降水來促進自身的生長,會比濕潤地區(qū)的植被表現(xiàn)得對降水更為敏感。因此,本文將以生態(tài)區(qū)為視角,使用NDVI植被指數(shù)表征植被生長狀態(tài),分析歐亞大陸植被生長對降水的響應關(guān)系在不同生態(tài)區(qū)的季節(jié)差異和該響應關(guān)系的變化趨勢,從而分析得到歐亞大陸植被生長對降水響應的時空格局,為歐亞大陸不同區(qū)域應對和適應氣候變化提供決策依據(jù)。

    1 數(shù)據(jù)和方法

    本文以歐亞大陸為研究區(qū)域,以陸域植被生態(tài)系統(tǒng)為研究對象,利用GIMMS(Global Inventory Modeling and Mapping Studies)長序列(1982—2015年)第3代植被指數(shù)NDVI3g數(shù)據(jù)來表征陸地生態(tài)系統(tǒng)的植被生長狀況。NDVI3g采用多顆AVHRR(Advanced Very High Resolution Radiometer)衛(wèi)星數(shù)據(jù)計算得到,數(shù)據(jù)的空間分辨率為0.0833度,每月的前15天產(chǎn)生一個a值,余下的天數(shù)的數(shù)據(jù)產(chǎn)生第二個b值,在數(shù)據(jù)發(fā)布時已經(jīng)做了初步的大氣校正[9]。為了最大限度去除云的干擾,在計算月度數(shù)據(jù)時,本研究將取每個月a、b兩值的最大值作為當月的NDVI值[10- 11]。相關(guān)研究區(qū)域的氣象數(shù)據(jù)來自英國東安哥拉大學氣候研究中心發(fā)布的空間分辨率為0.5度的全球網(wǎng)格數(shù)據(jù)集Climatic Research Unit Time-series data sets第4版本(CRUTS V4)[12]。

    根據(jù)歐亞大陸的地理位置,參考已有研究成果對季節(jié)的劃分方案[13],本文設(shè)置4至5月為春季,6至8月為夏季,9至10月為秋季,從而將月值植被數(shù)據(jù)和氣象數(shù)據(jù)整合生成春、夏、秋3個季節(jié)的均值。對季節(jié)平均NDVI值小于0.1的像素設(shè)置為無植被,將每個季節(jié)的NDVI和氣象數(shù)據(jù)做去趨勢化處理[13],再將氣象數(shù)據(jù)從0.5度重新采樣至與NDVI數(shù)據(jù)保持一致的空間分辨率(1/12度),參與其后的相關(guān)性計算[7]。最后使用世界野生動物基金會(World Wildlife Fund,WWF)發(fā)布的全球陸地生態(tài)區(qū)劃圖[14],對植被-降水關(guān)系在每個生態(tài)區(qū)上進行分異研究。WWF發(fā)布的生態(tài)區(qū)劃地圖是目前學界使用最為廣泛的區(qū)劃方案[15],在本文研究的歐亞大陸范圍內(nèi),共有328個生態(tài)區(qū)(eco-region)。

    為了單獨研究區(qū)域內(nèi)植被生長對降水響應的時空變化規(guī)律,本文通過計算偏相關(guān)系數(shù)(Partial Correlation Coefficient)[16],排除溫度和日照輻射的共同作用[6]。例如,在計算某像素點上的植被生長對降水的響應關(guān)系,可通過計算二階偏相關(guān)系數(shù),單獨計算植被指數(shù)(NDVI)與降水(Precipitation)的偏相關(guān)系數(shù)(RNDVI-Prep),并排除平均氣溫和日照輻射的共同作用,其計算公式為:

    (1)

    式(1)中,4個計算變量分別為:1=植被指數(shù)(NDVI);2=降水;3=平均溫度;4 =日照輻射(使用云量數(shù)據(jù)替代)。式(1)中的R12,3是排除平均溫度共同作用之后計算得到的NDVI與降水的偏相關(guān)系數(shù),該相關(guān)系數(shù)的計算公式如式(2)所示:

    (2)

    同理,參考公式(2),可計算得到R14,3和R24,3。其中R14,3是排除平均溫度共同作用之后計算得到的NDVI與日照輻射的偏相關(guān)系數(shù);R24,3是排除平均溫度共同作用之后計算得到的降水與日照輻射的偏相關(guān)系數(shù)。公式(2)中,R12為NDVI與降水的相關(guān)系數(shù),R13為NDVI與溫度的相關(guān)系數(shù),R23為降水與溫度的相關(guān)系數(shù)。

    本文首先分析1982—2015整個時間段上NDVI與降水的偏相關(guān)系數(shù)(RNDVI-Prep)的分布規(guī)律,并排除平均氣溫和日照輻射的共同作用。然后以17年為一個時間窗口,把1982—2015整個時間段共34年平均分兩個時間段。第一個17年為1982—1998;最后一個17年為1999—2015。然后使用時間滑動窗口,每次向后滑動一年,共產(chǎn)生18個時間窗口。在每個時間窗口i內(nèi)計算歐亞大陸陸域植被在不同季節(jié)對降水的響應關(guān)系,即排除溫度和日照輻射共同作用后的NDVI與降水的偏相關(guān)關(guān)系RNDVI-Prep(i)。將18個時間窗口的偏相關(guān)系數(shù)采用最小二乘法做趨勢分析[17],進而得到RNDVI-Prep(i)的斜率,即為TrendR。

    同時,為了表征每個像素上植被對降水的依賴程度和以及該依賴程度的變化趨勢,我們將像素分為四類。在每個像素(x,y)上,將18個滑動窗口上的偏相關(guān)系數(shù)RNDVI-Prep(i)做一次篩選,將偏相關(guān)系數(shù)為正,且置信度大于90%(即偏相關(guān)系數(shù)的絕對值|R|≥0.437)的第i個時間窗口選出來,并將T(i)標記為1,不滿足這個置信度的則標記為0,公式如下所示:

    (3)

    計算每個像素上ξpositive(x,y)的值,如果ξpositive(x,y)≥10,則表明該像素上不同時間滑動窗口偏相關(guān)系數(shù)RNDVI-Prep(i)為正且顯著相關(guān)(P<0.1)的時間段超過9個,即為總時間窗口的一半以上。這表明該像素上植被對降水的依賴程度較高。如果再疊加上植被與降水的偏相關(guān)關(guān)系RNDVI-Prep(i)的趨勢TrendR,則可以較為精確地描述每個地區(qū)的植被與降水的關(guān)系。例如,某個像素的ξpositive(x,y) ≥10,而且該像素上偏相關(guān)系數(shù)的趨勢TrendR為正,則表明像素上對降水的依賴程度較高,且在1982至2015年間有增強的趨勢。

    相反,在每個像素(x,y)上將18個滑動窗口上的相關(guān)系數(shù)為負且置信度大于90%的第i個時間窗口選出來,并將T(i)標記為1,不滿足這個置信度的則標記為0,公式如下所示:

    (4)

    計算每個像素上ξnegative(x,y)的值,如果ξnegative(x,y) ≥10,則表明該像素上不同時間滑動窗口RNDVI-Prep(i)為負且偏相關(guān)系數(shù)顯著相關(guān)(P<0.1)的時間段為總時間窗口數(shù)的一半以上。這表明該像素上植被對降水的依賴程度較低,且過多降水會導致溫度和日照輻射相應的減少,越不利于植被生長。

    2 水文氣象條件

    根據(jù)Nemani等[6]提出的植被生產(chǎn)力與光照、溫度、水分因子這3個自然變量的脅迫關(guān)系模型,本文將該模型由年尺度分解到每個季節(jié)。如圖1所示,春季受水分脅迫(紅色區(qū)域)影響最大的研究區(qū)域西起歐洲伊比利亞半島南部,貫穿地中海北岸、阿拉伯半島、伊朗高原、南亞次大陸,直至孟加拉灣地區(qū)。春季受日照輻射脅迫(綠色區(qū)域)的地區(qū)主要集中在歐洲、中國華南地區(qū)和東南亞地區(qū)。高緯度地區(qū)在春季主要受溫度脅迫(藍色區(qū)域)。在春季,中國北方地區(qū)、蒙古高原、中亞地區(qū)的植被同時受水分和溫度的脅迫(粉色區(qū)域)。

    圖1 植被生長所受自然脅迫條件的季節(jié)分布Fig.1 Distribution of seasonal potential climatic constraints to plant growth紅色代表水分脅迫區(qū)域(Water,W),藍色代表溫度(Temperature,T)脅迫區(qū)域,綠色代表日照輻射(Radiation,R)脅迫區(qū)域

    如圖1所示,在夏季,由于氣溫升高,俄羅斯境內(nèi)大部分地區(qū)的植被從春季時受溫度脅迫轉(zhuǎn)變?yōu)槭苋照蛰椛涿{迫(綠色區(qū)域),只有俄羅斯遠東少部分地區(qū)和中國的青藏高原部分地區(qū)依然保持受溫度脅迫。植被受水分條件脅迫的區(qū)域(紅色區(qū)域)主要集中在伊比利亞半島、亞平寧半島、巴爾干半島、阿拉伯半島和中亞等地區(qū)。其余地區(qū),包括中國東部和南部地區(qū)、南亞及東南亞地區(qū),其植被大多是受日照輻射脅迫。圖1顯示,相對于夏季,在秋季受水分脅迫的地區(qū)有向南移動的趨勢,在歐洲和中亞地區(qū)有所減少,但是在印度西北部地區(qū)呈增加趨勢。歐洲大部分地區(qū)和中國南部、日本、南亞及中南半島、東南亞地區(qū)的植被主要受日照輻射脅迫。相對于夏季,中國西北部和蒙古國地區(qū)的植被在秋季同時受水分和溫度脅迫的區(qū)域(粉色區(qū)域)有所擴展。

    通過上述3個自然脅迫條件的分季節(jié)分析可以發(fā)現(xiàn),在歐亞大陸植被生長對降水的響應有明顯的時空差異,其研究必須在時間尺度和空間尺度同時展開。

    3 季節(jié)變化

    通過計算1982—2015年整個時間段上NDVI與降水的偏相關(guān)系數(shù)RNDVI-Prep(排除了溫度與日照輻射的共同作用),得到歐亞大陸每個像素上春、夏、秋3個季節(jié)植被生長對降水的響應關(guān)系(圖2)。在春季,植被生長與降水呈顯著正相關(guān)關(guān)系(P<0.01)的地區(qū)與WWF發(fā)布的部分生態(tài)區(qū)邊界有較高的一致性,這些區(qū)域主要集中在東歐大草原(Pontic steppe,PS)、中亞地區(qū)的干旱草原與荒漠地區(qū)(例如:Central Asian northern desert,CAND),中國東北和蒙古國東部的草原區(qū)(例如:Mongolian-Manchurian grassland, MG)。與此相對,植被生長與降水的關(guān)系為顯著負相關(guān)(P<0.01)的區(qū)域所處的緯度相對正相關(guān)地區(qū)的緯度普遍較高,主要集中北歐瑞典、芬蘭以及俄羅斯境內(nèi)60°N以北的高緯度地區(qū)。這些地區(qū)的負相關(guān)可能與春季的降雪影響有關(guān),加之該地區(qū)在春季主要受到溫度脅迫,過多的降水會使得春季的溫度升溫來的緩慢,從而影響這些地區(qū)的植被生長。

    夏季由于溫度的增長,RNDVI-Prep呈顯著正相關(guān)(P<0.01)的區(qū)域明顯增加。這些區(qū)域廣泛地分布在35°—60°N之間(圖2中虛線所示)的歐亞大陸干旱、半干旱地區(qū)。它們西起東歐大草原,貫穿中亞草原和干旱地區(qū),東至我國東北和蒙古國地區(qū)的草原生態(tài)區(qū)。另外巴基斯坦和印度西部的干旱森林和荒漠地區(qū),以及印度南部疏林生態(tài)區(qū)的植被也表現(xiàn)出與降水呈顯著正相關(guān)關(guān)系。需要指出的是,在我國華北平原和中歐部分地區(qū),農(nóng)業(yè)區(qū)域占比很高,這些地區(qū)的植被生長與降水也呈現(xiàn)很高的正相關(guān)關(guān)系。與此相對,植被生長與降水呈負相關(guān)(P<0.01)的區(qū)域?qū)Ρ却杭久黠@減少,這些區(qū)域主要集中在北歐及俄羅斯境內(nèi)的部分高緯度地區(qū)。這些地區(qū)的植被在夏季雖然主要受到日照輻射脅迫,但是夏季的降水還是會導致區(qū)域性的降溫,因此可能對植被生長造成不利的影響。值得注意的是,在阿富汗南部和巴基斯坦西南部以及阿拉伯半島的荒漠地區(qū),也有部分區(qū)域表現(xiàn)出植被生長與降水的負相關(guān)關(guān)系,這可能主要是由于本文選取的最小NDVI閾值為0.1,并沒有很好地將該地區(qū)的非真NDVI像素剔除所致(圖2)。

    圖2 植被生長與降水的偏相關(guān)系數(shù)分布圖Fig.2 Distribution of partial correlation coefficient between vegetation growth and precipitationR=±0.447, R=±0.408, R=±0.349, R=±0.295, R=±0.232對應的置信度分別為99%, 98%, 95%, 90% 和80%;CAND,Central Asian northern desert;PS, Pontic steppe;MG,Mongolian-Manchurian grassland;YGT,Yamal-Gydan tundra;WST,West Siberian Taiga;EST East Siberian Taiga

    相對于春季和夏季,植被與降水的相關(guān)關(guān)系在秋季則更為集中。顯著正相關(guān)(P<0.01)的區(qū)域主要集中環(huán)黑海的中東歐地區(qū),印度西北部的干旱地區(qū),中國北方干旱的黃土高原和鄂爾多斯高原,以及印度尼西亞蘇門答臘和加里曼丹島南部的雨林地區(qū)。與此相對,植被生長與降水呈顯著負相關(guān)的地區(qū)主要分布在北歐斯堪的納維亞半島中東部的泰加森林(Taiga),并且與該半島南端的混合林生態(tài)區(qū)有著明顯的差別,可明顯看出WWF生態(tài)區(qū)劃邊界南北兩側(cè)的差異。在俄羅斯,西西伯利亞泰加森林(West Siberian Taiga,WST)和東西伯利亞泰加森林(East Siberian Taiga,EST)對降水的響應也有著明顯差異,西西伯利亞地區(qū)植被對降水的關(guān)系以負相關(guān)居多,東西伯利亞地區(qū)植被對降水的關(guān)系以正相關(guān)居多。值得關(guān)注的是,在俄羅斯極地地區(qū)的亞馬爾-格達苔原(Yamal-Gydan tundra,YGT),植被生長與降水表現(xiàn)出極強的正相關(guān)關(guān)系,與其他高緯度生態(tài)區(qū)的響應關(guān)系的表現(xiàn)極為不同。通過分析該生態(tài)區(qū)的植被與降水的時間序列可知,這主要是由于該地區(qū)近十多年來秋季降水顯著減少,并且植被指數(shù)NDVI在2013—2015也相對于往年減少較明顯所致。

    4 趨勢分析

    通過計算1982—2015年間的18個時間窗口的偏相關(guān)系數(shù)的斜率,進而得到每個像素上植被指數(shù)與降水的響應關(guān)系的趨勢(圖3),即為TrendR。如果一個像素上NDVI與降水的相關(guān)性的趨勢為正值,則表明該像素上的植被對降水的依賴在1982—2015年之間越來越高。反之,則表明該像素上的植被對降水的依賴在這34年之間在減小。如圖3所示,植被與降水的偏相關(guān)關(guān)系RNDVI-Prep(i)的趨勢TrendR在春、夏、秋3個季節(jié)上都表現(xiàn)出顯著的空間分異特征。

    圖3 1982—2015年間植被生長與降水的偏相關(guān)系數(shù)的趨勢分布圖Fig.3 Trend of the partial correlation coefficient between vegetation growth and precipitation from 1982 to 2015

    Piao等[18]通過計算植被生產(chǎn)力與溫度的偏相關(guān)系數(shù)的趨勢發(fā)現(xiàn)北半球植被與溫度的相關(guān)關(guān)系在減弱。同理,針對RNDVI-Prep(i)的趨勢TrendR分析雖然能表征植被對降水相關(guān)系數(shù)的變化趨勢,但通過判斷TrendR的正負關(guān)系和置信度來判斷植被生長對降水的響應關(guān)系是有一定的局限性的。這主要是因為18個滑動窗口得到如果是非顯著相關(guān)的RNDVI-Prep(i),但通過趨勢計算也可能得到一個顯著增強的TrendR。例如,當某上像素上植被與降水的偏相關(guān)關(guān)系由顯著負相關(guān)慢慢轉(zhuǎn)變?yōu)榉秋@著的正相關(guān),如果只依靠TrendR來描述該地區(qū)植被與降水的關(guān)系,便會表述成該地區(qū)植被對降水的依賴程度在增加,這顯然是不客觀、不準確的。因此,本文結(jié)合TrendR和本文數(shù)據(jù)方法部分介紹的ξ(x,y)這兩個要素來共同表征植被對降水的依賴程度,從而將研究區(qū)域的像素分成4類。

    表1中,第一類像素(在圖4中標記為紅色)表示在1982至2015年間,在18個滑動時間窗口上,RNDVI-Prep(i)呈顯著正相關(guān)(P<0.1)的時間窗口超過一半(即ξ(x,y)≥10),因此本文設(shè)定該像素的植被在大部分時間段對降水形成顯著依賴,而且由于TrendR為正,這說明該像素上植被生長對降水的依賴程度在增加。

    表1 植被生長對降水的依賴程度分類

    圖4 植被生長對降水的依賴程度Fig.4 Dependency of the vegetation growth to precipitationBMF,Baltic mixed forests;MG,Mongolian-Manchurian grassland;SRT,Scandinavian and Russian taiga;WST,West Siberian taiga;PS,Pontic steppe;KS,Kazakh steppe;TD,Thar Desert;SDF,South Deccan Plateau dry deciduous forests

    同理,第二類像素(橙色)表征的是植被在大部分時間段對降水形成顯著依賴,但是TrendR為負,這說明該像素上植被生長對降水的依賴程度在減小。依此類推,第三類像素(綠色)表征的是植被在大部分時間段對降水形成顯著相斥,但是TrendR為正,這說明該像素上植被生長對降水的響應關(guān)系在由相斥轉(zhuǎn)為逐步依賴。第四類像素(藍色)表征的是植被在大部分時間段對降水是顯著相斥,而且TrendR為負,這說明該像素上植被生長對降水表現(xiàn)出相斥的程度在繼續(xù)加大。

    通過對春、夏、秋3個季節(jié)植被對降水依賴程度的分析(圖4),可以看出在春季植被生長對降水表現(xiàn)出相斥的區(qū)域零星地分布在高緯度地區(qū)(圖4),包括英國南部和北歐斯堪的納維亞半島部分生態(tài)區(qū),以及俄羅斯的西伯利亞部分區(qū)域。與此相對,植被對降水表現(xiàn)出依賴的地區(qū)主要集中在35°N至60°N之間的地區(qū)。在歐洲表現(xiàn)最為明顯的生態(tài)區(qū)是波蘭西北部的波羅的海混合林(Baltic mixed forests,BMF)。另外,在中亞地區(qū)的干旱草原和荒漠地區(qū),中國東北和蒙古國東部的草原區(qū)(MG),植被不僅對降水在大部分時間段表現(xiàn)出顯著依賴(橙色像素),而且部分地區(qū)的依賴程度還在不斷的加強(紅色像素)。

    相對于春季,在夏季(圖4),植被對降水表現(xiàn)出依賴的地區(qū)(紅色和橙色像素)緯度有明顯抬升,這些地區(qū)主要集中在東歐大草原(Pontic steppe,PS)和中亞地區(qū)的哈薩克草原(Kazakh steppe,KS)以及中國東北和蒙古國東部的草原區(qū)(MG)。另外,在巴基斯坦和印度西北部地區(qū)的塔爾沙漠(Thar Desert,TD)和南部的南德干高原干旱落葉森林(South Deccan Plateau dry deciduous forests,SDF),以及印度尼西亞的雨林等地區(qū),植被對降水都表現(xiàn)出依賴,但是依賴程度在減弱(橙色像素)。在夏季植被對降水表現(xiàn)出相斥的地區(qū)(綠色和藍色像素)主要集中在60°N以北更高緯度生態(tài)區(qū),如斯堪的納維亞和俄羅斯泰加林(Scandinavian and Russian taiga, SRT)和西西伯利亞泰加林區(qū)(West Siberian taiga, WST)。

    由圖4所示,在秋季,植被對降水表現(xiàn)出依賴的地區(qū)相對于春季和夏季則表現(xiàn)的較為集中。在歐洲,這些地區(qū)主要集中在東歐大草原(PS)和環(huán)黑海的眾多混合林生態(tài)區(qū)。在亞洲,植被對降水表現(xiàn)出依賴的地區(qū)主要集中在巴基斯坦和印度西北部的干旱荒漠地區(qū)(例如,塔爾沙漠,TD)和印度尼西亞蘇門答臘島和婆羅洲南部的雨林地區(qū)。中國的黃土高原混合林和鄂爾多斯高原草地生態(tài)區(qū)的植被對降水也有很強的依賴關(guān)系,且依賴程度在增強。與此相對,在歐洲地區(qū),植被對降水表現(xiàn)出相斥的區(qū)域(綠色和藍色像素)主要集中在法國東部的西歐闊葉林區(qū)和瑞典中部的泰加林區(qū)地區(qū)。在亞洲,植被對降水表現(xiàn)出相斥的地區(qū)主要集中在中亞的干旱荒漠地區(qū)以及俄羅斯的西西伯利亞泰加林區(qū)。

    5 討論

    本文使用了目前最新的第3代植被指數(shù)NDVI3g數(shù)據(jù)來表征歐亞大陸陸域植被生長狀態(tài)。與之前文獻所采用的方法不同[7, 13, 18],本研究采用最新的全球空間分辨率為0.5度的氣象格點數(shù)據(jù)(CRUTS V4)通過最簡單的臨近采樣法重新采樣至與植被數(shù)據(jù)精度(1/12度)保持一致,參與其后的相關(guān)系數(shù)計算。這樣的處理方法的有兩個優(yōu)點:一方面將對植被數(shù)據(jù)的細節(jié)進行最大程度的保留;另一方面,考慮到氣象要素中,除了降水以外,氣溫和日照輻射在一定的范圍內(nèi)都有很好的空間一致性。因此把0.5度的氣象數(shù)據(jù)采樣成1/12度能最大程度上保留植被指數(shù)信息的細節(jié),得到更精細的分析結(jié)果。限于篇幅,在未來可考慮采用ANUSPLIN等專業(yè)氣象插值軟件,對比不同采樣方法和插值工具對研究結(jié)果的影響。

    另外,本文通過改進Nemani等[6]提出的氣象條件對植被生長的脅迫模型,將該模型分解到春、夏、秋3個季節(jié),分析了歐亞大陸不同季節(jié)上氣象條件對植被生長脅迫的基本分布情況。進一步,本文以像素為基本研究對象,通過偏相關(guān)關(guān)系計算,分析了歐亞大陸每個像素上植被指數(shù)NDVI對降水在不同季節(jié)的響應關(guān)系,并同時排除了溫度和日照輻射的共同作用。研究表明這些地區(qū)植被生長對水資源的響應在不同季節(jié)、不同地點表現(xiàn)出明顯的時間和空間分異特征,并且其空間差異的邊界與現(xiàn)有的生態(tài)區(qū)邊界在某些地區(qū)吻合極好,有可能成為未來生產(chǎn)更精細化的生態(tài)區(qū)劃圖的一個重要參考。

    通過統(tǒng)計1982—2015整個時間段上植被NDVI與降水的響應關(guān)系,我們發(fā)現(xiàn)在春季(圖2),RNDVI-Prep為顯著正相關(guān)的地區(qū)基本限制在35°—60°N之間的東歐大草原、中亞地區(qū)的干旱草原與荒漠地區(qū)以及中國東北和蒙古國東部的草原區(qū)。而在夏季(圖2),由于溫度升高,蒸發(fā)加大,植被與降水呈正相關(guān)的地區(qū)相對于春季明顯增加。而且這些區(qū)域西起東歐大草原,東至我國東北和蒙古國地區(qū)的草原生態(tài)區(qū),貫穿了整個歐亞大陸的干旱和半干旱地區(qū)。不僅如此,甚至在60°N以北的高緯度泰加林生態(tài)區(qū),很多地區(qū)的植被也表現(xiàn)出與降水的正相關(guān)關(guān)系。相對于春季和夏季,植被與降水的相關(guān)關(guān)系在秋季(圖2)則更為集中,主要集中環(huán)黑海的中東歐地區(qū),南亞次大陸西北部的干旱地區(qū),以及印度尼西亞的雨林地區(qū)。在中國,北方干旱的黃土高原混合林生態(tài)區(qū)和鄂爾多斯高原草甸構(gòu)成一條東北至西南走向的植被與降水呈正相關(guān)關(guān)系的生物帶。

    為了更好地了解NDVI與降水的響應關(guān)系在不同時間段上的變化,本文使用了時間滑動窗口來分析該響應關(guān)系的在每個季節(jié)上的變化趨勢。我們發(fā)現(xiàn)在春、夏、秋3個季節(jié),植被生長與降水響應關(guān)系在滑動時間窗口上的變化趨勢表現(xiàn)出極大的空間異質(zhì)特性(圖3)。在不同地點不同季節(jié),其植被對降水的依賴程度不同,依賴程度隨時間變化的趨勢也不同??紤]到不同時間窗口對偏相關(guān)系數(shù)計算結(jié)果的影響,本文同時使用偏相關(guān)關(guān)系的趨勢TrendR和RNDVI-Prep(i)呈顯著相關(guān)的所有時間窗口的個數(shù)ξ(x,y)來共同表征植被對降水的依賴程度,將植被對降水敏感的區(qū)域分為四類(圖4)。35°N以北地區(qū),在春秋兩季,植被生長對降水呈依賴關(guān)系地區(qū)的緯度要低于那些植被對降水呈相斥關(guān)系地區(qū)的緯度。通過對比圖4,可以明顯發(fā)現(xiàn)在夏季受增溫影響,蒸發(fā)加大,植被對降水呈依賴關(guān)系的區(qū)域比春季要多,而且有向高緯度地區(qū)擴張的趨勢。35°N以南地區(qū),植被對降水呈依賴關(guān)系的地區(qū)主要集中在夏、秋兩季的南亞次大陸西北部的干旱、半干旱地區(qū)及東南亞的雨林地區(qū)。

    本研究目前只考慮了氣象因素對植被生長的影響,并未考慮諸如人工施肥與灌溉、病蟲害、霜凍、火災等影響。例如,本文在圖2中發(fā)現(xiàn)俄羅斯極地地區(qū)的亞馬爾-格達苔原(Yamal-Gydan tundra,YGT)的植被和降水的響應關(guān)系在秋季與極地其他地區(qū)明顯不同,但是這是否也是和該地區(qū)馴鹿的過度養(yǎng)殖[19]有關(guān)還有待進行進一步研究。

    6 結(jié)論

    采用1982—2015年GIMMS NDVI3g數(shù)據(jù)和氣象格點數(shù)據(jù)對歐亞大陸不同生態(tài)區(qū)上植被對降水響應的季節(jié)規(guī)律進行了分析,得到如下結(jié)論:

    (1)在春節(jié),植被對降水的響應關(guān)系(RNDVI-Prep)為顯著正相關(guān)的地區(qū)廣泛地分布在東歐和中亞地區(qū)的草原及荒漠地帶,以及中國東北和蒙古國東部的草原地區(qū);在夏季,RNDVI-Prep為顯著正相關(guān)的地區(qū)廣泛地分布在35°N—60°N之間的歐亞大陸干旱、半干旱地區(qū)以及南亞次大陸西北和南部地區(qū);在秋季,RNDVI-Prep為顯著正相關(guān)的地區(qū)分布則相對于春夏兩季顯得更為集中,主要分布在環(huán)黑海的中東歐地區(qū),印度西北部的干旱地區(qū),中國北方干旱的黃土高原和鄂爾多斯高原,以及印度尼西亞蘇門答臘和加里曼丹島南部的雨林地區(qū)。

    (2)在歐亞大陸,夏季植被呈現(xiàn)對降水依賴的地區(qū)所處緯度要高于春、秋兩季。中亞干旱、半干旱生態(tài)區(qū)的植被在春、夏兩季對降水依賴較強;南亞次大陸西北部地區(qū)的干旱、半干旱及東南亞的雨林地區(qū)的植被生長在夏、秋兩季對降水依賴較強。

    (3)以生態(tài)區(qū)為視角,了解歐亞大陸的植被對降水響應在不同季節(jié)的變化規(guī)律可以為不同區(qū)域應對和適應氣候變化提供決策依據(jù)。

    猜你喜歡
    歐亞大陸生態(tài)區(qū)降水
    黑龍江省玉米生長季自然降水與有效降水對比分析
    黑龍江氣象(2021年2期)2021-11-05 07:07:00
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經(jīng)向環(huán)流
    新時代中俄全面合作與歐亞大陸經(jīng)濟空間再拓展
    俄羅斯研究(2020年3期)2020-08-11 09:09:28
    如何在幾分鐘之內(nèi)往返歐亞大陸
    少兒科技(2019年10期)2019-09-10 16:07:27
    考古發(fā)現(xiàn)的舌頭琴與歐亞大陸東部的古代游牧文化
    為什么南極降水很少卻有很厚的冰層?
    家教世界(2018年16期)2018-06-20 02:22:00
    大都市生態(tài)區(qū)綜合規(guī)劃思路研究
    江西建材(2018年4期)2018-04-10 12:36:58
    降水現(xiàn)象儀模擬軟件設(shè)計與實現(xiàn)
    打造社區(qū)O2O生態(tài)區(qū) 彩生活住宅新模式問世
    IT時代周刊(2015年8期)2015-11-11 05:50:49
    世界
    生物進化(2015年4期)2015-02-14 02:12:08
    国产综合懂色| 国产男女超爽视频在线观看| 久久99精品国语久久久| 97在线视频观看| 最后的刺客免费高清国语| 欧美丝袜亚洲另类| 91精品一卡2卡3卡4卡| 国产成人精品久久久久久| 少妇人妻 视频| 午夜亚洲福利在线播放| 女人被狂操c到高潮| 人妻 亚洲 视频| 亚洲第一区二区三区不卡| 日韩精品有码人妻一区| 久久韩国三级中文字幕| 国产91av在线免费观看| 欧美+日韩+精品| 神马国产精品三级电影在线观看| 九色成人免费人妻av| 一级毛片我不卡| 久久久久久久国产电影| 欧美激情久久久久久爽电影| 国产黄频视频在线观看| 超碰av人人做人人爽久久| 亚洲精品第二区| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 少妇人妻 视频| 亚洲,一卡二卡三卡| 欧美丝袜亚洲另类| 天天一区二区日本电影三级| 精品人妻视频免费看| 国产精品爽爽va在线观看网站| 久久久久网色| 国产欧美日韩精品一区二区| 男女啪啪激烈高潮av片| 男女国产视频网站| 麻豆精品久久久久久蜜桃| 日韩欧美精品免费久久| 精品午夜福利在线看| 大话2 男鬼变身卡| 亚洲国产色片| av网站免费在线观看视频| 大又大粗又爽又黄少妇毛片口| 久久久精品免费免费高清| 日本一本二区三区精品| 亚洲图色成人| 联通29元200g的流量卡| 久久精品国产亚洲av涩爱| 成年版毛片免费区| 国产毛片a区久久久久| 国产男女超爽视频在线观看| 三级男女做爰猛烈吃奶摸视频| 人人妻人人澡人人爽人人夜夜| 国产一区有黄有色的免费视频| av在线app专区| 80岁老熟妇乱子伦牲交| 亚洲成色77777| 成人毛片a级毛片在线播放| 精品国产三级普通话版| 久久久精品94久久精品| 国语对白做爰xxxⅹ性视频网站| 简卡轻食公司| 国产乱人偷精品视频| 日韩成人伦理影院| 亚洲欧美精品自产自拍| 3wmmmm亚洲av在线观看| 亚洲精品国产成人久久av| av国产精品久久久久影院| 午夜福利视频精品| 午夜激情福利司机影院| 国产综合懂色| 久久久精品免费免费高清| 好男人在线观看高清免费视频| 国产黄a三级三级三级人| 午夜免费鲁丝| 美女cb高潮喷水在线观看| 能在线免费看毛片的网站| 简卡轻食公司| 狠狠精品人妻久久久久久综合| 中文字幕免费在线视频6| 可以在线观看毛片的网站| 人人妻人人澡人人爽人人夜夜| 日韩欧美精品v在线| 欧美xxⅹ黑人| av播播在线观看一区| 97在线视频观看| 好男人视频免费观看在线| 女人十人毛片免费观看3o分钟| 国产av国产精品国产| 欧美zozozo另类| 免费电影在线观看免费观看| 黄片无遮挡物在线观看| freevideosex欧美| 日产精品乱码卡一卡2卡三| 亚洲av国产av综合av卡| 男女边摸边吃奶| 日韩一本色道免费dvd| 日韩av不卡免费在线播放| 久久热精品热| 亚洲精品久久久久久婷婷小说| 日本午夜av视频| 久久人人爽av亚洲精品天堂 | 老师上课跳d突然被开到最大视频| 女人久久www免费人成看片| 26uuu在线亚洲综合色| 国产免费视频播放在线视频| 亚洲av成人精品一区久久| 日本免费在线观看一区| 精品久久久久久久久av| 国产精品久久久久久精品古装| 欧美日本视频| 精品国产三级普通话版| 日韩亚洲欧美综合| 欧美变态另类bdsm刘玥| 亚洲精品日韩av片在线观看| 国语对白做爰xxxⅹ性视频网站| tube8黄色片| 亚洲国产高清在线一区二区三| 国产精品一及| 中文乱码字字幕精品一区二区三区| 亚洲三级黄色毛片| 免费高清在线观看视频在线观看| 亚洲欧美日韩卡通动漫| 尤物成人国产欧美一区二区三区| 成人毛片60女人毛片免费| 国产国拍精品亚洲av在线观看| 亚洲国产色片| 久久精品国产亚洲av涩爱| 久久精品国产自在天天线| 18禁在线无遮挡免费观看视频| 亚洲图色成人| 日韩人妻高清精品专区| 69人妻影院| 亚洲精品国产av蜜桃| 国产精品.久久久| 亚洲综合色惰| 亚洲精品成人av观看孕妇| 中国美白少妇内射xxxbb| 精华霜和精华液先用哪个| 国产高清有码在线观看视频| 丝袜喷水一区| 2022亚洲国产成人精品| 一二三四中文在线观看免费高清| 国产综合懂色| 免费高清在线观看视频在线观看| 一区二区三区免费毛片| 看黄色毛片网站| 日本与韩国留学比较| 色视频在线一区二区三区| 中文字幕免费在线视频6| 美女脱内裤让男人舔精品视频| 亚洲精品色激情综合| 亚洲国产色片| 国产精品一区二区三区四区免费观看| 国产成年人精品一区二区| 男人爽女人下面视频在线观看| 成人午夜精彩视频在线观看| 美女主播在线视频| 在线 av 中文字幕| 十八禁网站网址无遮挡 | 亚洲av欧美aⅴ国产| 国产亚洲最大av| 免费看日本二区| 国产白丝娇喘喷水9色精品| 免费电影在线观看免费观看| 国产亚洲一区二区精品| 最近中文字幕2019免费版| xxx大片免费视频| 国产 精品1| 日日摸夜夜添夜夜添av毛片| 日日撸夜夜添| 日韩一区二区三区影片| 国产男女超爽视频在线观看| 嫩草影院精品99| 亚洲人成网站在线播| av天堂中文字幕网| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品| 1000部很黄的大片| 一二三四中文在线观看免费高清| a级一级毛片免费在线观看| 国产精品偷伦视频观看了| 亚洲aⅴ乱码一区二区在线播放| 中文在线观看免费www的网站| 天天躁夜夜躁狠狠久久av| 久久人人爽av亚洲精品天堂 | 久久女婷五月综合色啪小说 | 99热这里只有精品一区| 欧美精品人与动牲交sv欧美| 国产免费视频播放在线视频| 赤兔流量卡办理| 精品人妻视频免费看| 最后的刺客免费高清国语| 成年女人在线观看亚洲视频 | 高清在线视频一区二区三区| 国产精品福利在线免费观看| 好男人在线观看高清免费视频| 亚洲国产精品999| 国产精品国产三级专区第一集| av在线播放精品| 老师上课跳d突然被开到最大视频| 日韩人妻高清精品专区| 一级毛片我不卡| 久热这里只有精品99| 欧美日韩在线观看h| 国产免费一区二区三区四区乱码| 少妇的逼水好多| 少妇被粗大猛烈的视频| 在线免费观看不下载黄p国产| 人妻制服诱惑在线中文字幕| 亚洲精品国产成人久久av| 午夜精品国产一区二区电影 | 亚洲国产精品成人综合色| 久久久久久久亚洲中文字幕| 精华霜和精华液先用哪个| 亚洲av免费高清在线观看| 熟女av电影| 亚洲成人中文字幕在线播放| 久久久久精品性色| 国产精品一区二区三区四区免费观看| 噜噜噜噜噜久久久久久91| 久热久热在线精品观看| av网站免费在线观看视频| 亚洲,一卡二卡三卡| 国产日韩欧美在线精品| 日韩中字成人| 99热6这里只有精品| .国产精品久久| 国产一区二区亚洲精品在线观看| 亚洲精品一区蜜桃| av免费观看日本| 亚洲精品国产av成人精品| 国产91av在线免费观看| 麻豆久久精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 午夜爱爱视频在线播放| 三级男女做爰猛烈吃奶摸视频| 综合色av麻豆| 夜夜爽夜夜爽视频| 亚洲欧美成人综合另类久久久| 狂野欧美白嫩少妇大欣赏| 国产精品一区二区在线观看99| 国产男女内射视频| 国产乱人偷精品视频| 亚洲真实伦在线观看| 97热精品久久久久久| 国产成人午夜福利电影在线观看| av又黄又爽大尺度在线免费看| 女人被狂操c到高潮| 制服丝袜香蕉在线| 欧美成人一区二区免费高清观看| 搞女人的毛片| 国产亚洲午夜精品一区二区久久 | 国产成人精品婷婷| 国产综合懂色| 菩萨蛮人人尽说江南好唐韦庄| 99热全是精品| 身体一侧抽搐| 午夜老司机福利剧场| 内射极品少妇av片p| 又粗又硬又长又爽又黄的视频| 中文字幕制服av| 男人舔奶头视频| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 青青草视频在线视频观看| 成人一区二区视频在线观看| 国产精品一及| 国产精品蜜桃在线观看| 亚洲真实伦在线观看| 一本一本综合久久| 在线天堂最新版资源| 久久久久久久精品精品| 日韩在线高清观看一区二区三区| 久久精品综合一区二区三区| 免费少妇av软件| 国产伦精品一区二区三区视频9| 久久人人爽人人片av| 寂寞人妻少妇视频99o| 欧美zozozo另类| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 亚洲精品一区蜜桃| 少妇人妻 视频| 天堂俺去俺来也www色官网| 午夜福利在线观看免费完整高清在| 国精品久久久久久国模美| 欧美日本视频| 日韩av免费高清视频| 精品亚洲乱码少妇综合久久| 精品国产乱码久久久久久小说| 久久久久国产精品人妻一区二区| 国内精品美女久久久久久| 王馨瑶露胸无遮挡在线观看| 网址你懂的国产日韩在线| 中国国产av一级| 夜夜爽夜夜爽视频| 在线观看三级黄色| 国产日韩欧美在线精品| 亚州av有码| 内地一区二区视频在线| 最新中文字幕久久久久| 国产v大片淫在线免费观看| freevideosex欧美| 99热6这里只有精品| 1000部很黄的大片| 国产老妇伦熟女老妇高清| 亚洲成色77777| 91在线精品国自产拍蜜月| 波多野结衣巨乳人妻| 亚洲精品乱久久久久久| 中文字幕人妻熟人妻熟丝袜美| 日本午夜av视频| 麻豆精品久久久久久蜜桃| 22中文网久久字幕| av女优亚洲男人天堂| 波野结衣二区三区在线| 色播亚洲综合网| 色5月婷婷丁香| 国产亚洲午夜精品一区二区久久 | 嫩草影院精品99| 亚洲丝袜综合中文字幕| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 国产男人的电影天堂91| 一级爰片在线观看| 大码成人一级视频| 欧美一级a爱片免费观看看| av免费在线看不卡| 观看免费一级毛片| 在线观看一区二区三区| 亚洲成人久久爱视频| 亚洲精华国产精华液的使用体验| 中文字幕av成人在线电影| 卡戴珊不雅视频在线播放| 久久精品熟女亚洲av麻豆精品| 亚洲天堂av无毛| videossex国产| 日韩在线高清观看一区二区三区| 亚洲国产精品专区欧美| 丝袜脚勾引网站| 日本熟妇午夜| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 狂野欧美白嫩少妇大欣赏| av网站免费在线观看视频| 好男人在线观看高清免费视频| 欧美三级亚洲精品| 久久久色成人| 午夜福利高清视频| 丝袜脚勾引网站| 国产女主播在线喷水免费视频网站| 能在线免费看毛片的网站| 人人妻人人澡人人爽人人夜夜| 91久久精品国产一区二区成人| 国产一区二区三区综合在线观看 | 男人爽女人下面视频在线观看| 天天一区二区日本电影三级| 免费看a级黄色片| 国产精品国产三级国产专区5o| 国产黄片美女视频| 免费观看性生交大片5| 熟妇人妻不卡中文字幕| 欧美xxⅹ黑人| 久久人人爽av亚洲精品天堂 | 中文字幕人妻熟人妻熟丝袜美| 欧美日韩在线观看h| 成年人午夜在线观看视频| 亚洲最大成人中文| 国产 一区精品| 亚洲人与动物交配视频| 欧美日韩一区二区视频在线观看视频在线 | 日韩欧美精品v在线| 一本一本综合久久| 黄色日韩在线| 亚洲欧美日韩卡通动漫| 日日摸夜夜添夜夜添av毛片| 国国产精品蜜臀av免费| 丝袜喷水一区| 国产精品av视频在线免费观看| 国产69精品久久久久777片| 1000部很黄的大片| 少妇猛男粗大的猛烈进出视频 | 男女国产视频网站| 欧美 日韩 精品 国产| 黄色一级大片看看| 日本猛色少妇xxxxx猛交久久| 女的被弄到高潮叫床怎么办| 国产永久视频网站| 一个人观看的视频www高清免费观看| 亚洲av中文字字幕乱码综合| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载 | 别揉我奶头 嗯啊视频| 国产亚洲5aaaaa淫片| 男插女下体视频免费在线播放| 久久久久久久久久成人| 国产av不卡久久| 国产乱人偷精品视频| 亚洲av国产av综合av卡| 人妻系列 视频| 欧美亚洲 丝袜 人妻 在线| 国产乱来视频区| 欧美亚洲 丝袜 人妻 在线| 欧美丝袜亚洲另类| 亚洲精品视频女| 在线免费观看不下载黄p国产| 亚洲真实伦在线观看| 亚洲成人一二三区av| 少妇被粗大猛烈的视频| videos熟女内射| 日韩大片免费观看网站| 国产色婷婷99| 欧美少妇被猛烈插入视频| 国产精品不卡视频一区二区| 黄色怎么调成土黄色| 伦理电影大哥的女人| 成人无遮挡网站| 在线天堂最新版资源| 好男人在线观看高清免费视频| 久久热精品热| 婷婷色av中文字幕| 精品久久久久久电影网| 精品一区二区三卡| 国产精品爽爽va在线观看网站| 麻豆成人av视频| 精品少妇久久久久久888优播| 亚洲最大成人中文| 一区二区三区乱码不卡18| 亚洲国产高清在线一区二区三| 国产淫片久久久久久久久| 国产爱豆传媒在线观看| 婷婷色av中文字幕| 精品久久久久久电影网| 联通29元200g的流量卡| av一本久久久久| 特级一级黄色大片| 毛片一级片免费看久久久久| 欧美xxⅹ黑人| 久久久久性生活片| 最新中文字幕久久久久| 久久久久久久亚洲中文字幕| 91久久精品国产一区二区三区| 中文资源天堂在线| 欧美zozozo另类| 精品国产露脸久久av麻豆| 男人爽女人下面视频在线观看| 黄色怎么调成土黄色| 人妻制服诱惑在线中文字幕| 18禁裸乳无遮挡免费网站照片| 国产午夜福利久久久久久| 亚洲国产精品成人久久小说| 蜜臀久久99精品久久宅男| 晚上一个人看的免费电影| 日本wwww免费看| 国产 一区 欧美 日韩| 精品久久久久久久久av| 内射极品少妇av片p| 18禁在线无遮挡免费观看视频| 寂寞人妻少妇视频99o| 在线免费观看不下载黄p国产| 日本色播在线视频| 中文字幕人妻熟人妻熟丝袜美| 一边亲一边摸免费视频| 在线精品无人区一区二区三 | av线在线观看网站| 大码成人一级视频| 国产伦理片在线播放av一区| 丝袜美腿在线中文| 免费观看av网站的网址| 久久久久性生活片| 亚洲怡红院男人天堂| 国产色婷婷99| 青青草视频在线视频观看| 久久女婷五月综合色啪小说 | 大陆偷拍与自拍| 在现免费观看毛片| 91aial.com中文字幕在线观看| 直男gayav资源| 国产av不卡久久| 人人妻人人爽人人添夜夜欢视频 | 亚洲av成人精品一二三区| 中文在线观看免费www的网站| 国产一区二区亚洲精品在线观看| 成人高潮视频无遮挡免费网站| 日韩免费高清中文字幕av| 亚洲激情五月婷婷啪啪| 91精品国产九色| 日韩大片免费观看网站| 91精品一卡2卡3卡4卡| 久久6这里有精品| 国产午夜福利久久久久久| 麻豆乱淫一区二区| 久久99热这里只有精品18| 精品视频人人做人人爽| 啦啦啦中文免费视频观看日本| 99精国产麻豆久久婷婷| 一级黄片播放器| 日韩中字成人| 青春草视频在线免费观看| 我要看日韩黄色一级片| 成年版毛片免费区| 丝袜脚勾引网站| 精华霜和精华液先用哪个| 国产精品麻豆人妻色哟哟久久| 一级片'在线观看视频| 亚洲无线观看免费| videossex国产| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲av天美| 国产黄a三级三级三级人| 伊人久久精品亚洲午夜| 高清欧美精品videossex| 只有这里有精品99| 久久久久久久久久久丰满| 国产黄a三级三级三级人| 天天一区二区日本电影三级| 男人和女人高潮做爰伦理| 国产一区二区三区av在线| 真实男女啪啪啪动态图| 免费大片18禁| 天堂中文最新版在线下载 | 亚洲av欧美aⅴ国产| 欧美 日韩 精品 国产| 在线观看三级黄色| 国产探花极品一区二区| 春色校园在线视频观看| 黄片无遮挡物在线观看| 国产一区有黄有色的免费视频| 欧美一区二区亚洲| 亚洲电影在线观看av| 久久人人爽人人片av| 亚洲伊人久久精品综合| 亚洲精华国产精华液的使用体验| 成人无遮挡网站| 亚洲精品久久久久久婷婷小说| 日日啪夜夜爽| 欧美性感艳星| 天堂中文最新版在线下载 | 日韩av不卡免费在线播放| 免费黄网站久久成人精品| 春色校园在线视频观看| 国产午夜精品久久久久久一区二区三区| 伊人久久国产一区二区| 听说在线观看完整版免费高清| 白带黄色成豆腐渣| 亚洲va在线va天堂va国产| 欧美xxⅹ黑人| 狂野欧美白嫩少妇大欣赏| 成人美女网站在线观看视频| 日韩人妻高清精品专区| 午夜激情福利司机影院| 亚洲精品日韩av片在线观看| 精品99又大又爽又粗少妇毛片| 欧美日韩在线观看h| 性色av一级| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 久久99蜜桃精品久久| 欧美 日韩 精品 国产| 99热国产这里只有精品6| 久久久久久久久久久丰满| 高清日韩中文字幕在线| 国产色爽女视频免费观看| 午夜激情久久久久久久| 欧美变态另类bdsm刘玥| 超碰av人人做人人爽久久| 日韩国内少妇激情av| 嫩草影院入口| 你懂的网址亚洲精品在线观看| 老司机影院成人| 久久久午夜欧美精品| 国产精品久久久久久精品古装| 日韩不卡一区二区三区视频在线| 国产 一区 欧美 日韩| 精品久久久久久久久av| 日日摸夜夜添夜夜添av毛片| 亚洲aⅴ乱码一区二区在线播放| 国产熟女欧美一区二区| 麻豆成人午夜福利视频| 91午夜精品亚洲一区二区三区| 国产精品伦人一区二区| 久久99热这里只有精品18| 国产免费又黄又爽又色| 色哟哟·www| 黄色配什么色好看| 狂野欧美白嫩少妇大欣赏| 亚洲av.av天堂| 一个人看视频在线观看www免费| 下体分泌物呈黄色| www.色视频.com| 在线精品无人区一区二区三 | 香蕉精品网在线| 国国产精品蜜臀av免费| 九九在线视频观看精品| 欧美丝袜亚洲另类| 亚洲自偷自拍三级| 伦理电影大哥的女人| 日韩一区二区三区影片| 在线观看三级黄色| 国产黄片美女视频| 国产亚洲一区二区精品| 久久人人爽人人爽人人片va| 黄色日韩在线| 免费看av在线观看网站| 国产精品三级大全| 老司机影院毛片| 国产精品久久久久久精品古装| 一级a做视频免费观看| 九九在线视频观看精品| 欧美成人精品欧美一级黄| 久久ye,这里只有精品|