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

    近35年青藏高原植被帶變化對(duì)氣候變化及人類活動(dòng)的響應(yīng)

    2019-05-28 06:36:00魏彥強(qiáng)蘆海燕王金牛王旭峰
    草業(yè)科學(xué) 2019年4期
    關(guān)鍵詞:青藏高原海拔高原

    魏彥強(qiáng),蘆海燕,王金牛,孫 建,王旭峰

    (1.中國(guó)科學(xué)院西北生態(tài)環(huán)境資源研究院 甘肅省遙感重點(diǎn)實(shí)驗(yàn)室,甘肅 蘭州 730000;2.蘭州大學(xué) 經(jīng)濟(jì)學(xué)院,甘肅 蘭州 730000;3.蘭州財(cái)經(jīng)大學(xué) 會(huì)計(jì)學(xué)院,甘肅 蘭州 730101;4.中國(guó)科學(xué)院成都生物研究所,四川 成都 610041;5.中國(guó)科學(xué)院地理科學(xué)與資源研究所,北京 100101)

    近些年來(lái),當(dāng)代全球暖溫化問(wèn)題(contemporary global warming),尤其是在高海拔、高緯度等氣候變化敏感區(qū)的環(huán)境變化已引起了全球的廣泛關(guān)注[1-2]。在高緯度和高海拔地區(qū)的積雪、冰川、凍土和植被等對(duì)氣候變化的響應(yīng)顯得尤為明顯,是較為典型和敏感的氣候指示因子[3-4]。目前的研究結(jié)果表明,高海拔地區(qū)氣候暖化較低海拔地區(qū)明顯[5-7],全面的暖化對(duì)寒冷的高海拔地區(qū)并不意味著都是負(fù)面的影響,例如增溫使得嚴(yán)寒地區(qū)開(kāi)始改善植物的生長(zhǎng)環(huán)境,有利于植物向高海拔地區(qū)遷移和動(dòng)物及人類的進(jìn)入,是全球變暖中物種的避難所[8-10]。因此,對(duì)高海拔地區(qū)氣候變化的研究在當(dāng)代氣候變化研究中具有代表性和典型性。

    氣候制約著植被的地理分布,植被是區(qū)域氣候特征的反映和指示,二者之間存在密不可分的聯(lián)系。植被帶往往被認(rèn)為是氣候帶的典型指示物。在大尺度的植被特征研究中,植被帶的概念與因海拔、經(jīng)度、緯度變化而產(chǎn)生的水熱梯度差異相結(jié)合,以指示植被類型的劃分[11-15]。由于植物的分布受控于氣候模式,最低溫度對(duì)限制物種的分布范圍和向南北極高緯度地區(qū)及高海拔地區(qū)的極化擴(kuò)展起著決定性的作用[16-17]。可以預(yù)測(cè),全球全面暖化將改變物種分布的適宜環(huán)境和生存范圍的極限,因而使得物種在暖化過(guò)程中得到重新分布[16-17]。當(dāng)前在高海拔地區(qū)關(guān)于全球變暖引起的分布界限變化研究主要集中于以下幾個(gè)方面:1)雪線/零物質(zhì)平衡線 (snowline / equilibrium line altitude, ELA)隨冰川退縮的變動(dòng)[18-19];2)樹(shù)線/林線(timberline)和生態(tài)交錯(cuò)帶(ecotone)隨全球變暖向兩極高緯度地區(qū)及高海拔地區(qū)的極化移動(dòng)[20-22];3)氣候帶(climate belt)或物種(species)向兩極高緯度地區(qū)及高海拔地區(qū)的遷移[8, 16-17];4)土地利用變化 (land-use and land-cover change, LUCC)[23-25]。Colwell等[16]和 Jump 等[17]研究認(rèn)為,由于在垂直方向上氣溫遞減率的存在以及在南北方向上溫度隨緯度的遞減,在全球增溫情形下,樹(shù)線的分布范圍將向高緯度高海拔地區(qū)擴(kuò)展。植被的重分布及其他物種的遷移將隨之發(fā)生[21, 24, 26]。與此同時(shí),暖溫化使得物種的適宜生存條件下限隨之改變,人口的遷移及以游牧為主的牲畜也隨之?dāng)U展其活動(dòng)范圍,表現(xiàn)為向兩極及高海拔地區(qū)的遷入,人類活動(dòng)增強(qiáng)[16, 27-28]。因此高山區(qū)和高緯度地區(qū)生態(tài)系統(tǒng)對(duì)全球暖化的響應(yīng)往往表現(xiàn)為植被帶的分布范圍及其面積的變化,以及植被生長(zhǎng)狀況的變化。而目前在該領(lǐng)域長(zhǎng)時(shí)間序列的監(jiān)測(cè)和研究,歸一化植被指數(shù) (normalized difference vegetation index, NDVI)被認(rèn)為是一個(gè)很好的監(jiān)測(cè)指標(biāo)[9, 29]。

    青藏高原 (Qinghai-Tibetan Plateau, QTP)因其對(duì)氣候變化的高敏感性和在氣候變化中其氣候因子的指示意義較同緯度低海拔地區(qū)更具代表性和典型性而倍受關(guān)注[30-31]。根據(jù)最近的研究結(jié)果,青藏高原地區(qū)在過(guò)去的幾十年中經(jīng)歷了一個(gè)加速的暖化過(guò)程,其降水變化較為平緩,而潛在蒸散能力則有下降趨勢(shì)[32-34]。根據(jù)最新對(duì)青藏高原冰川及凍土退化和環(huán)境效應(yīng)的調(diào)查顯示,從1950s開(kāi)始青藏高原暖溫化趨勢(shì)明顯,其凍土和山岳冰川遭受了大規(guī)模的退縮和面積縮減[35-37]。大量事實(shí)已表明,因高原的全面暖化,其雪線/零物質(zhì)平衡線及樹(shù)線的位置因冰川退縮而得到了抬升,植被覆蓋發(fā)生了明顯的變化[19, 22]。氣候條件的變化已經(jīng)引起了其生態(tài)環(huán)境和生態(tài)系統(tǒng)的變化[38-39],而人為因素如人口增加、基礎(chǔ)設(shè)施工程建設(shè)、城鎮(zhèn)化發(fā)展、過(guò)度放牧等的干擾已對(duì)生態(tài)系統(tǒng)產(chǎn)生了明顯的影響,并在一些區(qū)域超過(guò)了氣候變化影響本身[27-28, 40]。另外,大量的氣候模擬情景和預(yù)測(cè)結(jié)果表明,青藏高原在21世紀(jì)將經(jīng)歷持續(xù)的暖化過(guò)程,并且會(huì)進(jìn)一步緩解其嚴(yán)寒而干冷的氣候環(huán)境[5, 41-42]。目前已有很多關(guān)于青藏高原植被變化及其對(duì)氣候變化響應(yīng)的研究[38, 42-44],但通過(guò)林線及生態(tài)交錯(cuò)帶位置變化而研究植被帶變動(dòng),以及除了考慮氣候因素之外,綜合考慮人為活動(dòng)影響的研究還相對(duì)較少[27-28]。且在人為影響因素的解釋上,指標(biāo)的選取具有較強(qiáng)的任意性,概念化、定性化的描述較多而定量化的數(shù)值模擬和分析較為少見(jiàn)[27-28, 35],往往與氣候影響難以剝離,缺乏客觀性評(píng)價(jià)[43-45]。近些年來(lái),隨著人口規(guī)模、基礎(chǔ)設(shè)施建設(shè)及城鎮(zhèn)化的發(fā)展,越來(lái)越多的研究表明,人類活動(dòng)增強(qiáng)是一個(gè)不可忽視的重要影響方面,在區(qū)域環(huán)境變化中往往超過(guò) 氣 候 變 化 影 響 本 身[24, 27-28, 44, 46], 而 目 前 在 青 藏 高原地區(qū)對(duì)該方面的研究相對(duì)較少,相關(guān)研究較為薄弱。因此,本研究從NDVI的變化出發(fā),分析植被帶變化趨勢(shì)和氣候變化、人類活動(dòng)之間的關(guān)系,進(jìn)而剖析植被帶變化對(duì)氣候變化及人類活動(dòng)的響應(yīng)。

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

    1.1 研究區(qū)域

    本研究選擇氣候變化較為劇烈且垂直地帶性較為典型的青藏高原地區(qū)作為研究區(qū)。該區(qū)域海拔普遍較高,且擁有相對(duì)較為寬闊的海拔范圍、垂直帶譜較為發(fā)育,其對(duì)氣候變化響應(yīng)較為敏感。依據(jù)Zhang等[47]界定的包括其南部的喜馬拉雅山脈、北部的昆侖山及祁連山、西至帕米爾高原及喀喇昆侖山、東至成都平原的范圍目前在學(xué)界普遍比較接受。本研究以此為參考,南部以中國(guó)國(guó)界線為界,其北部及東部選擇海拔在3 500 m以上的地區(qū)為其邊緣。本研究中青藏高原總面積約2.57 ×106km2,平均海拔 4 390 m 左右,其經(jīng)度在 70°-105° E,緯度為 25°-40° N(圖 1)。

    1.2 數(shù)據(jù)

    1.2.1 NDVI數(shù)據(jù)

    由于植被地上生物量及凈初級(jí)生產(chǎn)力與歸一化植被指數(shù)(NDVI)之間有很高的相關(guān)性,借用遙感數(shù)據(jù)對(duì)地上植被狀況的分析通常以NDVI來(lái)代替植被覆蓋狀況[9, 48-50]。本研究選擇美國(guó)航空航天局(NASA)全球監(jiān)測(cè)與模型研究組GIMMS (Global Inventory Monitoring and Modelling Studies)提 供的NOAA AVHRR-NDVI-3g V1.0 數(shù)據(jù)集 (https://ecocast.arc.nasa.gov/data/pub/gimms/3g.v1/),其空間分辨率為8 km × 8 km,屬半月合成數(shù)據(jù),時(shí)間跨度為 1981 年7月1日-2015年12月31日,共計(jì)69個(gè)NetCDF-4半月合成圖像文件。該數(shù)據(jù)集的NDVI已經(jīng)過(guò)幾何精糾正、輻射校正、大氣校正等預(yù)處理,且已采用最大值合成法(MVC)以減少云、大氣、太陽(yáng)高度角的影響。

    1.2.2 氣候及其他數(shù)據(jù)

    本研究的氣象臺(tái)站數(shù)據(jù)由中國(guó)氣象局國(guó)家氣候中心提供,時(shí)間范圍為1961-2015年。其中位于高原上的站點(diǎn)91個(gè),因站點(diǎn)變動(dòng),記錄缺失等原因,其中有完整月平均氣溫和月降水記錄的站點(diǎn)87個(gè)。由于高原上的站點(diǎn)分布較少,尤其是廣大的中西部地區(qū)長(zhǎng)時(shí)間序列的站點(diǎn)較少,本研究另外選取了環(huán)境數(shù)據(jù)分析中心(CEDA)的長(zhǎng)序列再分析數(shù)據(jù) CRUTS 4.01 (http://catalogue.ceda.ac.uk/intro)。該數(shù)據(jù)集包括了每天的平均溫度、降水等常用指標(biāo),數(shù)據(jù)包括1900年1月1日到2016年12月31日的逐日氣象要素,空間分辨率為 0.5° × 0.5°。人口、牲畜數(shù)量等經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)均來(lái)自各個(gè)省市的統(tǒng)計(jì)年鑒。由于該類數(shù)據(jù)均以行政單元(縣)來(lái)統(tǒng)計(jì),但青藏高原是一個(gè)自然地理邊界,在高原上的部分以其所占總縣域面積的比例來(lái)進(jìn)行分割。另外DEM 數(shù)據(jù)由 NASA提供,空間分辨率為90 m × 90 m。植被覆蓋/土地利用類型數(shù)據(jù)為基于GLC2000的中國(guó)科學(xué)院1∶100萬(wàn)土地利用數(shù)據(jù)集(圖1A);生態(tài)區(qū)劃數(shù)據(jù)采用Zheng[51]擬定的中國(guó)地理生態(tài)區(qū)域系統(tǒng)的框架方案,將青藏高原分為11個(gè)自然區(qū)(圖1B,表1),包含了3個(gè)溫度帶和4個(gè)干濕區(qū)。

    圖1 青藏高原地區(qū)根據(jù)GLC2000分類的植被覆蓋/土地利用類型(A),多年平均NDVI以及依照鄭度等[51]擬定的生態(tài)氣候分區(qū)(B)Figure 1 Vegetation types/land covers in the Tibetan Plateau according to the GLC2000 categories (A),the mean NDVI in recent years and eco-climatic zones in the TP according to [51] (B)

    1.3 數(shù)據(jù)處理

    首先對(duì)GIMMS-NDVI的半月合成數(shù)據(jù)通過(guò)國(guó)際通用的最大值合成法 MVC (maximum value composites)合成月NDVI,該方法可以進(jìn)一步消除云、大氣、太陽(yáng)高度角等部分干擾。并利用12個(gè)月月數(shù)據(jù)合成4個(gè)季節(jié)NDVI數(shù)據(jù)和年數(shù)據(jù)。合成公式為:

    式(1)中,NDVIi為半月NDVI數(shù)據(jù)或已經(jīng)合成的月值數(shù)據(jù),NDVIm為合成后的NDVI季節(jié)數(shù)據(jù)或年數(shù)據(jù)??紤]到青藏高原地區(qū)植被覆蓋普遍較低,且大多為溫帶草原及高寒草甸,本研究取NDVI的閥值為0.05[9, 52],即小于此值的地區(qū)被認(rèn)為是無(wú)植被帶覆蓋。為清楚地了解NDVI在空間上的變化趨勢(shì),逐柵格的趨勢(shì)分析能直觀地反映其在空間上的變化趨勢(shì)。一般地,可以通過(guò)計(jì)算植被綠度變化率 (greenness rate of change, GRC)來(lái)表征植被帶變化趨勢(shì)[10, 52]。而常用的則是采用容易實(shí)現(xiàn)的最小二乘法擬合年均NDVI的變化斜率[10]。但最小二乘法線性回歸的缺點(diǎn)則是樣本量較小時(shí),極易受到異常觀測(cè)點(diǎn)的影響。本研究采用了Theil-Sen線性趨勢(shì)分析法。該方法是一種非參數(shù)的估計(jì)算法,能有效消除長(zhǎng)時(shí)間序列分析中的異常值的影響[53-54],是對(duì)最小二乘線性回歸方法的改進(jìn)。其計(jì)算公式為:

    式(2)中:slope為植被綠度變化率,變量i,j為年的序號(hào),即i,j= 1, 2, …, 35。n為研究的時(shí)間序列長(zhǎng)度,n= 35;NDVIi為第i年的年最大 NDVI值。其中slope> 0說(shuō)明NDVI在這35年中的變化趨勢(shì)是增加的,反之則是減少的,其值的大小反映出植被帶變化的劇烈程度。以上線性趨勢(shì)分析均在PyCharm平臺(tái)上利用Python 3.7.1計(jì)算完成。

    2 結(jié)果

    2.1 植被帶的年際變化

    時(shí)間序列的變化分析可以從植被帶演化更替的角度來(lái)分析整個(gè)青藏高原整體的狀況。為了綜合分析植被帶隨時(shí)間的變化,本研究計(jì)算了整個(gè)青藏高原NDVI平均值的歷年變化情況,并按季節(jié)平均值統(tǒng)計(jì)其在過(guò)去35年的變化。從統(tǒng)計(jì)結(jié)果來(lái)看,NDVI在以月為尺度的變化中呈現(xiàn)出倒“U”字形的周期性變化特點(diǎn),即年內(nèi)的變化在冬春兩季較低,平均值約為0.16,而夏秋兩季達(dá)到年內(nèi)植被生長(zhǎng)的最好水平,平均值夏季為0.29,秋季為0.26左右。年內(nèi)的變化顯現(xiàn)出地上生物量與四季的更替及氣溫的變化較為一致的特點(diǎn)。從總體變動(dòng)趨勢(shì)上看,最近35年的NDVI年際變化較小,且波動(dòng)不大,總體上處于平穩(wěn)的波動(dòng)態(tài)勢(shì)。

    表1 青藏高原生態(tài)地理系統(tǒng)分區(qū)Table 1 Eco-climate zones in the Qinghai-Tibetan Plateau

    由于青藏高原整體上全年氣溫普遍較低,冬季漫長(zhǎng)而嚴(yán)寒,生長(zhǎng)季(5-9月)相對(duì)短暫,因此植被帶的發(fā)育主要集中在以夏季為主的生長(zhǎng)季。從生長(zhǎng)季的NDVI變化趨勢(shì)來(lái)看,過(guò)去35年生長(zhǎng)季的NDVI呈現(xiàn)出一定的波動(dòng)性,但總體趨勢(shì)平緩,最大值出現(xiàn)在1988年,為0.296;最小值出現(xiàn)在1995年,為0.267;多年生長(zhǎng)季平均值為0.282左右,有輕微的上升趨勢(shì)。以上結(jié)果表明,一方面,較低的NDVI反映了青藏高原總體上植被覆蓋較差,整個(gè)柴達(dá)木盆地、藏北高原等地區(qū)植被分布稀疏,幾近于裸地。由于生長(zhǎng)季節(jié)較短,地上生物量普遍相對(duì)較小,使得區(qū)域平均值較低。另一方面,青藏高原內(nèi)部植被類型多樣、植被帶較多、生態(tài)氣候區(qū)種類較多,氣候的垂直地帶性分布較為明顯,復(fù)雜的地形和氣候類型狀況使得區(qū)域間的植被發(fā)育狀況差異明顯。有些因增溫而使得植被帶生長(zhǎng)條件改善的地區(qū),很可能因一些地區(qū)因生長(zhǎng)環(huán)境的惡化以致植被覆蓋減少,但在整個(gè)高原求其平均值時(shí)而被相互抵消掉,使得總體上植被帶變化趨勢(shì)波動(dòng)性不大。因此,以青藏高原為整體來(lái)研究其平均的植被帶生長(zhǎng)狀況的變化顯得意義甚微。為此,在空間上分區(qū)域和植被帶類型的生長(zhǎng)狀況變化研究顯得頗為必要。

    2.2 植被帶的空間變化

    2.2.1 總體變化

    逐柵格的趨勢(shì)分析能直觀地反映植被帶在空間上的變化趨勢(shì),能較好地反映出每個(gè)柵格在時(shí)間序列上的變化差異,因此是在空間上研究植被帶對(duì)氣候及人為因素響應(yīng)差異的有效途徑。利用公式(2)可以計(jì)算出最近35年植被帶的空間變化趨勢(shì)(圖2)。

    圖2 青藏高原1981-2015年間植被帶生長(zhǎng)狀況變化情況Figure 2 Vegetation trends during 1981-2015 in the Tibetan Plateau

    從整個(gè)區(qū)域來(lái)看,大部分地區(qū)的植被覆蓋趨于好轉(zhuǎn)。在高原西南部及中部、柴達(dá)木盆地周圍的廣大地區(qū)植被都出現(xiàn)了增長(zhǎng)情況。部分地區(qū)的年增長(zhǎng)率達(dá)到了0.002 5及以上。尤其是在H1C1、H1C2兩個(gè)高原亞寒帶高寒草甸草原區(qū)及H2C2藏南高山灌叢草原區(qū)的增長(zhǎng)尤為明顯。而處于惡化趨勢(shì)的植被帶分布于山南林芝、云南北部、川西北及青海湖附近的低海拔地區(qū),高原的東部邊緣是這一趨勢(shì)較為集中的區(qū)域。另外H1B1果洛那曲高原亞寒帶灌叢草甸區(qū)的植被覆蓋退化也十分明顯(圖2)。

    2.2.2 不同分區(qū)間的差異

    由于植被帶的分布與海拔、植物種類、氣候類型等緊密相關(guān),本研究依據(jù)海拔帶、植被覆蓋種類和生態(tài)氣候區(qū)統(tǒng)計(jì)在各個(gè)類型區(qū)中的植被帶生長(zhǎng)狀況變化情況(圖3)。

    從生態(tài)氣候區(qū)的內(nèi)部統(tǒng)計(jì)結(jié)果來(lái)看(圖3A1),在整個(gè)青藏高原的11個(gè)生態(tài)氣候區(qū)中,植被帶生長(zhǎng)狀況平均水平處于較明顯改善的有6個(gè)區(qū),其中有H1C1、H1C2和H1B1三個(gè)高原亞寒帶氣候區(qū)以及H2C2、H2D1、H2D2三個(gè)高原溫帶氣候區(qū)。而植被帶平均生長(zhǎng)狀況趨于退化的有4個(gè)區(qū),均為植被生長(zhǎng)條件相對(duì)較好的亞熱帶闊葉林區(qū)(H3)和2個(gè)溫度較高的高原溫帶針葉林(H2AB1和H2C1)及高原溫帶荒漠區(qū)(H2D3)。從統(tǒng)計(jì)結(jié)果來(lái)看,位于高原亞寒帶的H1D1昆侖山高寒荒漠區(qū)變化不明顯。由于趨于改善的生態(tài)區(qū)集中于高原亞寒帶和高原溫帶,而退化的區(qū)域處于平均溫度等相對(duì)較高的亞熱帶及高原溫帶地區(qū),因此推測(cè),全球暖化對(duì)改善高海拔地區(qū)嚴(yán)寒條件下的植被帶生長(zhǎng)狀況有積極作用,溫度升高有可能是樹(shù)線和林線在這一地區(qū)升高、植被趨于改善的主導(dǎo)因素。而在氣候條件相對(duì)較好的亞熱帶、高原溫帶地區(qū)的植被退化極有可能是升溫后蒸散發(fā)加劇或人為因素影響的結(jié)果。但具體原因還待進(jìn)一步詳細(xì)分析。

    為了進(jìn)一步的深入分析植被帶生長(zhǎng)狀況在空間上的變化差異,以植被覆蓋類型來(lái)統(tǒng)計(jì)各個(gè)類型區(qū)內(nèi)的植被帶變化情況顯得十分必要。植被帶的分布與植被覆蓋/土地利用變化之間有著直接的聯(lián)系,植被覆蓋/土地利用變化的種類界定了植被帶的基本類型[16-17, 55]。本研究以 GLC 2000 的用地分類系統(tǒng)為參考,對(duì)其6個(gè)大類的25個(gè)小類進(jìn)行重新合并調(diào)整,將水體和冰川等進(jìn)行歸并,最后形成12個(gè)土地覆蓋類型。以此類型為基礎(chǔ),分別對(duì)各個(gè)區(qū)內(nèi)的植被帶生長(zhǎng)狀況變化率進(jìn)行了統(tǒng)計(jì)(圖3B1)。

    從統(tǒng)計(jì)的結(jié)果可以看出,在12個(gè)大的植被覆蓋/土地利用類型中,有4個(gè)顯著增長(zhǎng)的類型區(qū),即L31高覆蓋度(> 50%)草地、L32中覆蓋度(20%~50%)草地、L4水體及L33低覆蓋度(5%~20%)草地。其中平均增長(zhǎng)速度最快的為L(zhǎng)33低覆蓋度(5%~20%)草地。從植被帶的平均變化速度來(lái)看,保持基本不變的為L(zhǎng)6沙漠/戈壁,這與實(shí)際的情況相吻合。而在其他所剩的幾個(gè)植被帶減小的區(qū)域中,L21森林、L22高覆蓋度灌叢及L23疏林地等的減少較為明顯,是3種減少最為劇烈的類型區(qū)。從植被覆蓋在空間上的分布來(lái)看,增長(zhǎng)的區(qū)域L31、L32、L33基本上是位于高海拔地區(qū)的草地,由于水體中在分類時(shí)已包含了冰川和積雪,其覆蓋度的增加推測(cè)應(yīng)當(dāng)是氣溫增加后雪線升高、植被帶的上線苔原線隨之升高的結(jié)果;而高海拔地區(qū)整體上植被的好轉(zhuǎn)則反映了對(duì)氣候轉(zhuǎn)暖后植被帶生存狀況改善、植被長(zhǎng)勢(shì)趨好及覆蓋范圍擴(kuò)大的結(jié)果。在減少的幾個(gè)區(qū)域中,L21森林、L22高覆蓋度灌叢及L23疏林地均分布于高原邊緣谷地及南部海拔相對(duì)較低、溫度、降水等生長(zhǎng)環(huán)境相對(duì)較好的地區(qū)。

    在以上幾個(gè)分析結(jié)果中,無(wú)論是基于生態(tài)氣候分區(qū)的分析還是從用地類型分區(qū)分析,均顯示出植被帶的空間變化和海拔之間的密切聯(lián)系,在不同的海拔上其變化趨勢(shì)表現(xiàn)出明顯的差異性。由于本研究的植被帶位置移動(dòng)假設(shè)是以植被帶的垂直帶譜變化為依據(jù),據(jù)此,按海拔分區(qū)來(lái)分析植被帶的空間變化顯得很有必要。本研究以500 m等間隔做海拔分區(qū),分別劃分出各個(gè)范圍的海拔帶,即 0-500 m、500-1 000 m 和 1 000-1 500 m等不同的分區(qū),分別對(duì)每個(gè)區(qū)內(nèi)植被帶的變化率進(jìn)行了統(tǒng)計(jì)(圖3C1)。從計(jì)算的結(jié)果來(lái)看,在劃分的15個(gè)海拔區(qū)中,有5個(gè)均表現(xiàn)出明顯的增長(zhǎng)趨勢(shì),除去最高海拔的一個(gè)分區(qū),幾個(gè)植被好轉(zhuǎn)的地區(qū)均集中分布于4 500-6 500 m的高海拔地區(qū),并且占整個(gè)高原絕大部分的面積(圖3C2)。而4 000-4 500 m似乎是植被帶變化的一個(gè)臨界位置,該范圍內(nèi)沒(méi)有統(tǒng)計(jì)意義上的明顯變化,而在該海拔以下各個(gè)植被帶均表現(xiàn)出了較為明顯的退化趨勢(shì)。從平均減少速率來(lái)看,1 000-2 500 m 之間的 3 個(gè)較低海拔區(qū)減少較快,處于明顯的退化水平。

    從以上結(jié)果可以看出,植被帶生長(zhǎng)狀況的變化與生長(zhǎng)季節(jié)、植被帶類型、生態(tài)氣候區(qū)和海拔等有關(guān)。以青藏高原為整體的時(shí)間序列分析并不能體現(xiàn)出植被帶的變化差異,而在考慮到按生態(tài)氣候區(qū)、植被覆蓋類型和海拔后,其變化的差異性過(guò)程才得到了明顯體現(xiàn)。以分區(qū)為視角的空間差異分析體現(xiàn)出青藏高原植被帶在最近35年的變化中與生態(tài)氣候區(qū)、植被覆蓋類型及海拔之間的密切關(guān)系。

    3 討論

    3.1 植被帶變化對(duì)氣候變化的響應(yīng)

    一般認(rèn)為,植被類型、生長(zhǎng)狀況及其分布受氣候因子的直接影響,在特定的氣候模式制約下在地理分布上形成特定的植被帶和生態(tài)系統(tǒng)類型[16, 55]。因而植被帶類型及其生長(zhǎng)狀態(tài)的變化是對(duì)一定氣候模式和生長(zhǎng)環(huán)境變化的直接反應(yīng),氣候條件的改變使得植被的種類和分布范圍隨之發(fā)生相應(yīng)的變化以適應(yīng)新的氣候特點(diǎn)[26, 55-56]。在眾多的氣候因子中,光照、溫度、降水、蒸散發(fā)等被認(rèn)為是植被生長(zhǎng)的最重要且起決定性作用的氣候因子,而通常以日照、溫度和降水的影響最為直接[57]??紤]到數(shù)據(jù)的可獲得性,本研究在分析植被帶與氣候變化的相互關(guān)系時(shí),選取了與之相對(duì)應(yīng)的年平均溫度和年總降水作為主要?dú)夂騾⒖家蜃印?/p>

    從臺(tái)站的氣象記錄來(lái)看(圖4),整個(gè)青藏高原的年平均溫度(87個(gè)站點(diǎn)的平均值)呈現(xiàn)出明顯的增長(zhǎng)趨勢(shì),在Theil-Sen法線性擬合的結(jié)果中,年溫度增長(zhǎng)率為 0.035 3 ℃·a-1(R2= 0.688,P< 0.001)。表明整個(gè)青藏高原最近幾十年經(jīng)歷了較為快速的增溫暖化過(guò)程,這一暖化過(guò)程與NDVI緩慢增長(zhǎng)的趨勢(shì)基本一致,可以認(rèn)為是對(duì)氣候暖化的正向響應(yīng)。

    與溫度相比較,年總降水量在總體上呈現(xiàn)出略微的增長(zhǎng),但增長(zhǎng)幅度較小,為每年0.831 mm(R2=0.041,P< 0.32),線性增長(zhǎng)趨勢(shì)不太明顯,呈現(xiàn)出較大的年際變化差異。因此,從整個(gè)高原總體的NDVI及溫度和降水變化的趨勢(shì)來(lái)看,NDVI的增長(zhǎng)趨勢(shì)與溫度有較好的一致性,但與降水增加的關(guān)系不太明顯。

    從同化的數(shù)據(jù)來(lái)看(圖5A、B),氣溫的增長(zhǎng)率在H1C2、H1C1、H2C2、H2AB1及H2DI等藏北高原、高原南部河谷及柴達(dá)木盆地等地區(qū)相對(duì)較高,而在H2C1、H1B1高原的東北部及果洛那曲亞寒帶高原草甸區(qū)的增長(zhǎng)率相對(duì)較小。而降水在整個(gè)高原的增加則主要集中在西藏和青海兩省 (區(qū))的中部地區(qū),尤其以H1C1、H2C2和H2C2等區(qū)的增加為主。而其他生態(tài)氣候區(qū)內(nèi)則基本上呈現(xiàn)出減少趨勢(shì)。以溫度和降水的變化區(qū)分氣候變化的類型,即暖濕、暖干、冷濕和冷干4種氣候變化模式。其中以暖濕為主的地區(qū)主要集中于高原中部及南部的H1C1、H2C2、H2C2等高海拔地區(qū)。在其他幾個(gè)暖干化為主的區(qū)域中NDVI變化出現(xiàn)了不同的變化態(tài)勢(shì),其中H2D1柴達(dá)木盆地周邊地區(qū)有較為顯著的增長(zhǎng),但在H2C1及H1B1兩區(qū)中總體上有不同程度的減少。因此,總體上可以判斷,氣候的暖濕化有利于植被的生長(zhǎng),植被帶在暖濕氣候條件下出現(xiàn)了明顯的恢復(fù)和擴(kuò)展;而暖干化則在帶來(lái)部分地區(qū)植被好轉(zhuǎn)的同時(shí),退化較為明顯。值得一提的是,在氣候條件較好的高原南部林芝地區(qū)及滇西北的NDVI減小趨勢(shì)明顯。雖然年均溫度上升較快,但近些年降水減少明顯,尤其是近些年的夏季干旱極有可能是這一地區(qū)植被退化的主要原因[58-59]。

    以上分析可以總結(jié)為,最近35年整個(gè)青藏高原的暖溫化對(duì)植被狀況的改善具有明顯的促進(jìn)作用,尤其是在高原中部及西南部的H1C1、H1C2、H2C2三個(gè)高海拔區(qū)暖濕化促進(jìn)作用明顯,但在H2AB1區(qū)中尤其是川西北地區(qū),無(wú)論是年均溫還是降水都不明顯。由此可以推斷,暖化的過(guò)程使得在高海拔地區(qū)的植被帶生長(zhǎng)環(huán)境改善,有利于植被的生長(zhǎng)和植被帶向高海拔,高緯度地區(qū)擴(kuò)展。但在H2AB1區(qū)中,尤其是川西北地區(qū),無(wú)論是年均溫還是降水都與NDVI變化的關(guān)系并不明顯,則說(shuō)明這些地區(qū)的植被帶變化原因中,氣候要素變化并不是其最主要原因,可能是人為因素影響等其他因素作用的結(jié)果。

    3.2 植被帶變化對(duì)人類活動(dòng)的響應(yīng)

    圖5 2010s年平均氣溫T(A)和平均降水量P(B)與1980s年平均氣溫T(A)和平均降水量P(B)的差值Figure 5 The mean temperature (A2) and mean precipitation increasing maps (B2) were generated from the differences between the two periods (Mean T/P2010-2016 - Mean T/P1981-1990)

    由于植被覆蓋狀況受氣候條件和人為干擾兩個(gè)方面的影響,尤其是在生態(tài)環(huán)境比較脆弱的青藏高原高寒草甸及草原生態(tài)區(qū),其生長(zhǎng)季較短,植被地上生物量較少,破壞后恢復(fù)能力弱,很容易受氣候及人為因素的影響而使植被帶退化[27, 55, 60]。除了氣候因素之外,人為因素的影響其實(shí)一直為人 們 所 關(guān) 注[24, 27-28, 43]。 由 于 青 藏 高 原 地 區(qū) 是 我 國(guó) 四大牧區(qū)之一,畜牧業(yè)在該區(qū)域中的比重歷來(lái)較高[28]。近些年來(lái),隨著氣候轉(zhuǎn)暖及西部大開(kāi)發(fā)等國(guó)家發(fā)展政策的實(shí)施,青藏鐵路的開(kāi)通等使得這一地區(qū)的遷移及流動(dòng)人口逐漸增多,人口規(guī)模持續(xù)地增加[35],旅游業(yè)等第三產(chǎn)業(yè)的發(fā)展使得該地區(qū)的城鎮(zhèn)化發(fā)展速度加快[44, 46]。由于人口的增加及城鎮(zhèn)的發(fā)展,伴隨人口增長(zhǎng)的是對(duì)奶類、肉類、毛皮制品等畜產(chǎn)品的需求增加,并且牧民的收入與牲畜數(shù)量直接相關(guān),收入需求的刺激使得畜牧業(yè)超規(guī)模發(fā)展,人類活動(dòng)對(duì)高寒草原草甸生態(tài)系統(tǒng)的影響在持續(xù)加強(qiáng)[28, 43-44]。近些年來(lái)在牲畜數(shù)量持續(xù)增長(zhǎng)的同時(shí),畜牧業(yè)的發(fā)展對(duì)高寒草甸草地生態(tài)系統(tǒng)的壓力進(jìn)一步加大,過(guò)度放牧、鼠蟲害等使草地退化的報(bào)道越來(lái)越多[35, 43-44]。為了分析人類活動(dòng)對(duì)青藏高原地區(qū)植被帶變化的影響,本研究選取了人為因素影響中最重要的兩個(gè)指標(biāo),即反映人為因素規(guī)模和城鎮(zhèn)規(guī)模的人口規(guī)模,以及代表畜牧業(yè)規(guī)模的牲畜數(shù)量,并以二者與NDVI變化間的相互關(guān)系進(jìn)行分析,進(jìn)而探討人類活動(dòng)對(duì)植被帶的影響。

    從統(tǒng)計(jì)的結(jié)果來(lái)看,近些年來(lái),青藏高原地區(qū)無(wú)論是人口規(guī)模還是牲畜規(guī)模都得到了較快的增長(zhǎng)(圖6)。人口數(shù)量從1971年的813萬(wàn)人很快增長(zhǎng)為2010年的1 553.5萬(wàn)人,平均年增長(zhǎng)率達(dá)到18.51萬(wàn)人·年-1,總?cè)丝谠鲩L(zhǎng)了近兩倍。而以羊單位換算的整個(gè)青藏高原牲畜的數(shù)量則由1971年的1.002 8 ×108羊單位增長(zhǎng)到 2010 年的 1.343 6 × 108羊單位,平均增長(zhǎng)速度為85.2萬(wàn)羊單位·年-1,總規(guī)模增長(zhǎng)了1.34倍。由于青藏高原地區(qū)以畜牧業(yè)為主的第一產(chǎn)業(yè)占有較大的比重,從第一產(chǎn)業(yè)的發(fā)展及整個(gè)經(jīng)濟(jì)的增長(zhǎng)規(guī)模來(lái)看,青藏高原地區(qū)近些年的經(jīng)濟(jì)增長(zhǎng)較為明顯,尤其是第二、三產(chǎn)業(yè)增長(zhǎng)速度一度處于加速階段(圖7)。另外近些年來(lái),以礦產(chǎn)開(kāi)采及能源開(kāi)發(fā)為主的第二產(chǎn)業(yè)以及以旅游和服務(wù)業(yè)為主的第三產(chǎn)業(yè)發(fā)展勢(shì)頭較快,其比重越來(lái)越高,并逐漸占據(jù)了經(jīng)濟(jì)發(fā)展的主體地位。另外人口的快速發(fā)展也進(jìn)一步地推動(dòng)了青藏高原城鎮(zhèn)數(shù)量和規(guī)模的發(fā)展,近些年城鎮(zhèn)化的發(fā)展速度較快[46]。

    由于青藏高原地區(qū)的植被生長(zhǎng)季節(jié)較短,高海拔高寒及缺氧的生態(tài)環(huán)境使得生態(tài)系統(tǒng)較為脆弱,極易受人為活動(dòng)的影響和干擾。另外與其他低海拔氣候較適宜地區(qū)相比,青藏高原地區(qū)無(wú)論是從植被的生物總量、活動(dòng)狀況、污染的自然降解水平還是破壞后的自我恢復(fù)能力等都比其他暖濕氣候條件下的生態(tài)系統(tǒng)差,脆弱的生態(tài)環(huán)境使得這一地區(qū)對(duì)人為因素的干擾尤其敏感[43, 60]。

    圖6 青藏高原地區(qū)1971-2010年人口規(guī)模及以羊單位換算的牲畜數(shù)量的增長(zhǎng)情況Figure 6 Population and sheep-unit livestock numbers for the Qinghai-Tibetan Plateau during 1971-2010

    圖7 1978-2010年青海、西藏地區(qū)三次產(chǎn)業(yè)的規(guī)模(以當(dāng)年價(jià)格計(jì)算)及城鎮(zhèn)化率Figure 7 The sizes of the three industries and urbanization rate of the Qinghai Province and Tibet

    圖8 青藏高原地區(qū)最近35年的人口(A)及牲畜(B)密度線性增長(zhǎng)率Figure 8 Mean population (A) and mean livestock increase rates (B) over the past 40 years

    為了從空間上探討人為因素對(duì)植被帶的影響,本研究分析了青藏高原人口的分布及近40年的人口增長(zhǎng)速度(圖8),以及以羊單位換算的整個(gè)青藏高原的牲畜增長(zhǎng)速度(圖8)。受高海拔及高寒缺氧等惡劣的自然條件的限制,人口基本上分布于海拔高度相對(duì)較低、植被狀況較好的低山河谷地區(qū),并在高原邊緣接近低海拔的地區(qū)有相對(duì)較為集中的分布。從生態(tài)氣候分區(qū)上來(lái)看,H2C1、H1B1、H2AB1及H3等地區(qū)的人口密度較高。而在高原廣闊的中西部地區(qū)人口分布則較為稀疏。從人口的年平均增長(zhǎng)率也可以看出,高原邊緣地區(qū)、藏南及中部地區(qū)的人口增長(zhǎng)較快。牲畜的空間分布與人口的分布相似,主要分布于海拔高度相對(duì)較低、氣候環(huán)境相對(duì)較好的低海拔河谷地區(qū),H2C1、H1B1、H2AB1及H2C2等地區(qū)是分布相對(duì)較為集中的地區(qū)。與人口分布不同的是,牲畜的增長(zhǎng)速度近些年在川西北及高原的西部高海拔地區(qū)增長(zhǎng)較快,而在高原中部的H1C1、H1B1兩個(gè)地區(qū)有減少的態(tài)勢(shì)。由于這兩個(gè)地區(qū)是三江源自然保護(hù)區(qū),近些年來(lái)的退牧還草等生態(tài)移民及植被保護(hù)政策可能是其主要原因。

    以上無(wú)論是人口規(guī)模的分布還是牲畜規(guī)模的分布,如果對(duì)比分析最近幾十年的植被帶生長(zhǎng)狀況在空間上的變化(圖2),則可以發(fā)現(xiàn),在H2C1、H1B1、H2AB1等人口及牲畜密度較高的地區(qū),雖然植被的生長(zhǎng)環(huán)境相較高原北部及西部高海拔地區(qū)為好,但與其植被改善相反,這些人口及牲畜密度較高的地區(qū)其植被卻普遍出現(xiàn)了統(tǒng)計(jì)意義上的明顯下降態(tài)勢(shì)。由于人為活動(dòng)對(duì)植被帶的影響除了人為的對(duì)礦山開(kāi)采、森林的砍伐及城鎮(zhèn)、基礎(chǔ)設(shè)施建設(shè)對(duì)土地的侵占之外,最為直接的則是以發(fā)展畜牧業(yè)為主的牲畜放牧對(duì)植被帶的影響[27, 43-44]。為了檢驗(yàn)人為因素對(duì)植被帶的壓力和影響,本研究按縣統(tǒng)計(jì)了最近35年來(lái)青藏高原人口及牲畜的密度及其與NDVI變化率之間的關(guān)系。由于密度的數(shù)量級(jí)差異較大,本研究在取密度的自然對(duì)數(shù)后與NDVI增長(zhǎng)率之間進(jìn)行了Theil-Sen線性回歸分析 (圖 9A,9B)。

    從在青藏高原上的207個(gè)縣區(qū)的統(tǒng)計(jì)分析結(jié)果來(lái)看,人口、牲畜密度對(duì)數(shù)和NDVI變化率之間均出現(xiàn)了較為一致的負(fù)相關(guān)性。其中NDVI變化率NDVItrend和人口密度自然對(duì)數(shù)ln(Population Density)之 間 的 回 歸 方 程 為y= -0.000 2x+ 0.000 2(R2=0.166,P< 0.001),NDVI變化率與牲畜密度自然對(duì)數(shù)ln(Livestock Density)之間的回歸關(guān)系為y= -0.000 3x+0.009(R2= 0.142,P< 0.001),二者的回歸方程檢驗(yàn)均較為顯著?;貧w結(jié)果表明,隨著人口、牲畜密度的增加,NDVI增長(zhǎng)率出現(xiàn)了明顯的遞減,反映出隨著人口、牲畜對(duì)植被帶的壓力加大后,植被帶出現(xiàn)了退化現(xiàn)象。這與在NDVI變化生態(tài)氣候分區(qū)統(tǒng)計(jì)中的結(jié)果基本一致。

    圖9 以及1981-2015年以縣為統(tǒng)計(jì)單位人口密度(取自然對(duì)數(shù))與NDVI增長(zhǎng)率之間的回歸關(guān)系(A)及牲畜密度(取自然對(duì)數(shù))與NDVI增長(zhǎng)率之間的回歸關(guān)系(B)Figure 9 The correlations between the NDVI trend and ln(population density) (A)/ln(livestock density) (B)during 1981-2015 in each county

    由于人口及牲畜分布較為集中的地區(qū)為H2C1、H2AB1及H3等海拔相對(duì)較低、氣候環(huán)境相對(duì)較好的地區(qū),但這些地區(qū)的植被生長(zhǎng)狀況卻普遍出現(xiàn)了退化。在分析人口/牲畜密度與NDVI變化率之間的關(guān)系后可知,氣候的暖化使這些地區(qū)嚴(yán)寒的氣候條件緩解,變得更適宜于人類的游牧和定居,人口的遷入和畜牧業(yè)的快速發(fā)展卻對(duì)這一地區(qū)的植被帶帶來(lái)了壓力和干擾、產(chǎn)生了負(fù)面的影響,干擾加大從而使植被帶總體上處于退化狀態(tài)。在這些地區(qū)植被帶生長(zhǎng)狀況惡化的原因中,人為因素的影響已經(jīng)超過(guò)了平穩(wěn)而緩慢的氣候變化影響本身,已成為其退化的主要原因。

    4 結(jié)論

    本研究以對(duì)氣候變化和人為影響較為敏感的青藏高原地區(qū)植被帶變化趨勢(shì)為視角,采用GIMMSNDVI來(lái)反映植被帶的生長(zhǎng)狀況,并就其時(shí)空間變化與氣候因子、人為影響之間的響應(yīng)關(guān)系進(jìn)行了分析,得出以下結(jié)論:

    ⅰ)最近35年來(lái)整個(gè)青藏高原的植被帶生長(zhǎng)狀況有緩慢的增長(zhǎng)趨勢(shì),這與近些年來(lái)的氣候全面暖化和降水的輕微增加有正向的相關(guān)性。

    ⅱ)在高原中部及西南部等高海拔地區(qū),暖濕化的氣候變化趨勢(shì)是其植被好轉(zhuǎn)的直接原因。即隨著氣候暖化,高海拔地區(qū)嚴(yán)寒的植物生長(zhǎng)環(huán)境得到改善,林線上升,植被帶向高海拔、高緯度地區(qū)擴(kuò)展,NDVI相應(yīng)地增加。

    ⅲ)在高原東北部及其東部邊緣、東南部等海拔高度相對(duì)較低、氣候環(huán)境相對(duì)較好的地區(qū)植被普遍出現(xiàn)惡化。與緩慢的氣候暖化相比,人口及牲畜數(shù)量的增加等人為影響增強(qiáng)是其主要原因。即隨著氣候的暖化,一些地區(qū)因嚴(yán)寒的氣候環(huán)境改善而開(kāi)始適宜于游牧及定居,人口遷入,人為活動(dòng)的影響加劇,當(dāng)超過(guò)了植被帶的承載能力時(shí),植被帶開(kāi)始退化,表現(xiàn)為NDVI相應(yīng)地減少。

    致謝:感謝兩位匿名審稿人對(duì)本文提出的修改建議。中國(guó)氣象局國(guó)家氣象信息中心(NMIC)及英國(guó)環(huán)境數(shù)據(jù)分析中心(CEDA)提供了氣象資料,在此一并致謝。

    猜你喜歡
    青藏高原海拔高原
    青藏高原上的“含羞花”
    高海拔的甘南
    高原往事
    迸射
    高原往事
    高原往事
    為了讓青藏高原的天更藍(lán)、水更綠、草原更美
    南美三大高海拔球場(chǎng)
    海拔4000米的父愛(ài)
    世界上海拔最高、最雄偉的山脈
    河北遙感(2014年4期)2014-07-10 13:54:58
    一级片'在线观看视频| 妹子高潮喷水视频| 黄色 视频免费看| 亚洲精品在线观看二区| 他把我摸到了高潮在线观看 | bbb黄色大片| 黑丝袜美女国产一区| 国产精品偷伦视频观看了| 免费观看a级毛片全部| 精品国产一区二区久久| 国产成人av激情在线播放| 成人精品一区二区免费| 黄色视频在线播放观看不卡| 精品国产超薄肉色丝袜足j| 免费少妇av软件| 亚洲视频免费观看视频| 99国产综合亚洲精品| 久久中文字幕人妻熟女| 宅男免费午夜| 少妇裸体淫交视频免费看高清 | 男人操女人黄网站| av网站在线播放免费| 久久久久久久久免费视频了| 国产精品久久久久成人av| 亚洲五月色婷婷综合| 18禁观看日本| 热99久久久久精品小说推荐| 欧美日韩亚洲综合一区二区三区_| 日韩熟女老妇一区二区性免费视频| 18禁美女被吸乳视频| 十八禁人妻一区二区| 中文字幕人妻熟女乱码| 多毛熟女@视频| 国产av国产精品国产| 欧美国产精品va在线观看不卡| 国产亚洲一区二区精品| av欧美777| 99re6热这里在线精品视频| 欧美黑人精品巨大| 99国产精品免费福利视频| 80岁老熟妇乱子伦牲交| 亚洲精品粉嫩美女一区| 久久精品国产亚洲av高清一级| 天天操日日干夜夜撸| 青草久久国产| 大码成人一级视频| 午夜免费鲁丝| 老熟妇乱子伦视频在线观看| 乱人伦中国视频| 两个人免费观看高清视频| 国产一区二区激情短视频| 叶爱在线成人免费视频播放| 热99re8久久精品国产| 亚洲国产中文字幕在线视频| 757午夜福利合集在线观看| 精品少妇久久久久久888优播| 亚洲中文日韩欧美视频| 国产亚洲av高清不卡| www.自偷自拍.com| 99国产精品免费福利视频| 最近最新中文字幕大全免费视频| 夜夜爽天天搞| 女人精品久久久久毛片| 高潮久久久久久久久久久不卡| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美激情在线| 精品一品国产午夜福利视频| 午夜福利在线免费观看网站| 成年人黄色毛片网站| 国产极品粉嫩免费观看在线| 黑人猛操日本美女一级片| 成人精品一区二区免费| 欧美 日韩 精品 国产| 一二三四社区在线视频社区8| 亚洲av第一区精品v没综合| 青草久久国产| 可以免费在线观看a视频的电影网站| 久久ye,这里只有精品| 真人做人爱边吃奶动态| 亚洲熟女毛片儿| 在线观看人妻少妇| 成人免费观看视频高清| 亚洲色图 男人天堂 中文字幕| 色综合欧美亚洲国产小说| 大码成人一级视频| 在线av久久热| 汤姆久久久久久久影院中文字幕| 欧美一级毛片孕妇| 欧美乱码精品一区二区三区| 12—13女人毛片做爰片一| 精品一区二区三卡| 十八禁网站免费在线| 一进一出好大好爽视频| 黄色毛片三级朝国网站| 中文字幕最新亚洲高清| 亚洲人成电影观看| 久久国产亚洲av麻豆专区| 69精品国产乱码久久久| 精品久久久久久久毛片微露脸| 免费观看人在逋| 美女高潮到喷水免费观看| 久久av网站| 精品人妻熟女毛片av久久网站| 久久亚洲精品不卡| 午夜福利乱码中文字幕| 久久精品aⅴ一区二区三区四区| 国产精品久久久人人做人人爽| 1024香蕉在线观看| 亚洲成av片中文字幕在线观看| 亚洲av电影在线进入| 国产伦人伦偷精品视频| 天天操日日干夜夜撸| 亚洲av日韩在线播放| 满18在线观看网站| 母亲3免费完整高清在线观看| 99re6热这里在线精品视频| 一区福利在线观看| 日韩 欧美 亚洲 中文字幕| 少妇被粗大的猛进出69影院| 黄片小视频在线播放| 99re在线观看精品视频| 国产精品国产av在线观看| 日本欧美视频一区| 亚洲欧美色中文字幕在线| 999精品在线视频| 成人特级黄色片久久久久久久 | 精品少妇黑人巨大在线播放| 日韩 欧美 亚洲 中文字幕| 欧美精品高潮呻吟av久久| 午夜福利视频在线观看免费| 亚洲 国产 在线| 午夜精品国产一区二区电影| 久久人人97超碰香蕉20202| 在线亚洲精品国产二区图片欧美| 老汉色∧v一级毛片| 国产高清视频在线播放一区| 国产深夜福利视频在线观看| 另类精品久久| 91字幕亚洲| 精品人妻1区二区| 午夜福利视频在线观看免费| 久久国产精品大桥未久av| 视频在线观看一区二区三区| 国产人伦9x9x在线观看| 色婷婷久久久亚洲欧美| 在线观看一区二区三区激情| 精品国产乱码久久久久久男人| 在线观看免费午夜福利视频| 国产精品一区二区在线不卡| 一级黄色大片毛片| 极品少妇高潮喷水抽搐| 欧美激情高清一区二区三区| 操美女的视频在线观看| 国产成人av激情在线播放| 乱人伦中国视频| 五月开心婷婷网| 久9热在线精品视频| 91大片在线观看| 欧美性长视频在线观看| 天堂俺去俺来也www色官网| 精品一区二区三卡| 在线观看免费午夜福利视频| 老司机深夜福利视频在线观看| 亚洲精品国产色婷婷电影| 国精品久久久久久国模美| 自线自在国产av| av有码第一页| 欧美一级毛片孕妇| 午夜精品国产一区二区电影| 久久精品国产亚洲av高清一级| 国产亚洲欧美精品永久| 一区二区三区激情视频| 国产在线视频一区二区| 99久久国产精品久久久| 无限看片的www在线观看| tube8黄色片| 纵有疾风起免费观看全集完整版| 中文字幕最新亚洲高清| 老汉色∧v一级毛片| 国产三级黄色录像| 日本wwww免费看| 中文字幕高清在线视频| 12—13女人毛片做爰片一| 日韩熟女老妇一区二区性免费视频| 亚洲国产av新网站| 岛国在线观看网站| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 无限看片的www在线观看| 日本欧美视频一区| 免费人妻精品一区二区三区视频| 成人精品一区二区免费| 国产不卡一卡二| 老鸭窝网址在线观看| 成人影院久久| 精品人妻1区二区| 国产成人精品无人区| 午夜久久久在线观看| 欧美性长视频在线观看| 久久久久精品人妻al黑| 国产区一区二久久| 王馨瑶露胸无遮挡在线观看| 国产真人三级小视频在线观看| 啪啪无遮挡十八禁网站| 高潮久久久久久久久久久不卡| 999久久久国产精品视频| 精品福利永久在线观看| 黄色怎么调成土黄色| 91字幕亚洲| 天堂中文最新版在线下载| www.熟女人妻精品国产| 国产精品久久久久久精品古装| 在线天堂中文资源库| 亚洲成a人片在线一区二区| 欧美精品啪啪一区二区三区| 亚洲第一青青草原| 777米奇影视久久| 亚洲五月色婷婷综合| 19禁男女啪啪无遮挡网站| 精品一区二区三区四区五区乱码| 高清欧美精品videossex| 欧美精品高潮呻吟av久久| 午夜福利免费观看在线| 2018国产大陆天天弄谢| 亚洲精品在线美女| 亚洲精品国产精品久久久不卡| 老司机午夜十八禁免费视频| 国产精品国产高清国产av | 啦啦啦免费观看视频1| 亚洲av日韩在线播放| 男女高潮啪啪啪动态图| 建设人人有责人人尽责人人享有的| 一本一本久久a久久精品综合妖精| 国产黄频视频在线观看| 久久 成人 亚洲| 国产日韩欧美亚洲二区| 色精品久久人妻99蜜桃| 黄色毛片三级朝国网站| 国产又色又爽无遮挡免费看| 久久久久久久精品吃奶| 日韩有码中文字幕| 亚洲国产av影院在线观看| 一个人免费看片子| 亚洲av成人不卡在线观看播放网| 欧美+亚洲+日韩+国产| 亚洲男人天堂网一区| 国产亚洲欧美在线一区二区| 69av精品久久久久久 | 国产亚洲欧美精品永久| 99精品在免费线老司机午夜| 国产精品久久久久久精品电影小说| 自线自在国产av| 男女床上黄色一级片免费看| 免费女性裸体啪啪无遮挡网站| 女人精品久久久久毛片| 中文字幕精品免费在线观看视频| 国产精品久久久人人做人人爽| 午夜精品久久久久久毛片777| 少妇猛男粗大的猛烈进出视频| 国产xxxxx性猛交| 超色免费av| 91精品国产国语对白视频| 久久久久久人人人人人| 肉色欧美久久久久久久蜜桃| 91老司机精品| 精品人妻1区二区| 天天操日日干夜夜撸| 丰满迷人的少妇在线观看| 91字幕亚洲| 狠狠狠狠99中文字幕| 欧美精品av麻豆av| 亚洲全国av大片| 黄频高清免费视频| 色视频在线一区二区三区| 精品亚洲成国产av| 交换朋友夫妻互换小说| 日日爽夜夜爽网站| 久久久国产欧美日韩av| 久久天躁狠狠躁夜夜2o2o| 脱女人内裤的视频| 一区福利在线观看| 亚洲人成伊人成综合网2020| 97人妻天天添夜夜摸| 亚洲久久久国产精品| 久久这里只有精品19| 一级a爱视频在线免费观看| 久久久水蜜桃国产精品网| 91麻豆av在线| 成人av一区二区三区在线看| 成人亚洲精品一区在线观看| 视频区图区小说| 一二三四在线观看免费中文在| 精品少妇黑人巨大在线播放| 黑人猛操日本美女一级片| 亚洲av成人不卡在线观看播放网| 久久久久国内视频| 精品少妇内射三级| 亚洲第一青青草原| 日韩大片免费观看网站| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 国产精品 欧美亚洲| 三上悠亚av全集在线观看| 大香蕉久久网| 久久免费观看电影| 99re在线观看精品视频| 大陆偷拍与自拍| 国产在线视频一区二区| 丝袜美腿诱惑在线| 国产成人欧美在线观看 | 波多野结衣一区麻豆| 一级毛片精品| 水蜜桃什么品种好| 精品人妻在线不人妻| 日韩视频一区二区在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 日本av手机在线免费观看| 大型黄色视频在线免费观看| 日韩熟女老妇一区二区性免费视频| 国产极品粉嫩免费观看在线| 精品国产乱子伦一区二区三区| 午夜视频精品福利| 中文字幕人妻丝袜一区二区| tocl精华| 久久99热这里只频精品6学生| 国产精品久久久久久精品古装| 他把我摸到了高潮在线观看 | 午夜老司机福利片| 最新美女视频免费是黄的| 亚洲欧美激情在线| 波多野结衣av一区二区av| a级片在线免费高清观看视频| 男女之事视频高清在线观看| 女人精品久久久久毛片| 淫妇啪啪啪对白视频| 亚洲国产欧美网| 蜜桃国产av成人99| 国产野战对白在线观看| 日韩熟女老妇一区二区性免费视频| 免费不卡黄色视频| 久久精品亚洲精品国产色婷小说| 十八禁高潮呻吟视频| 99国产综合亚洲精品| 精品午夜福利视频在线观看一区 | 国产亚洲一区二区精品| 亚洲av成人不卡在线观看播放网| 亚洲精品国产精品久久久不卡| 青草久久国产| 夫妻午夜视频| 高清欧美精品videossex| 国产亚洲精品第一综合不卡| 国产精品九九99| 久久99一区二区三区| 性少妇av在线| 午夜免费鲁丝| 丝袜在线中文字幕| 免费黄频网站在线观看国产| 亚洲男人天堂网一区| 日韩欧美免费精品| 99久久精品国产亚洲精品| 色老头精品视频在线观看| 欧美老熟妇乱子伦牲交| 菩萨蛮人人尽说江南好唐韦庄| 精品国产国语对白av| 亚洲国产成人一精品久久久| 国产精品九九99| 色婷婷av一区二区三区视频| av电影中文网址| 国产精品影院久久| 久久中文字幕一级| 国产精品亚洲av一区麻豆| 日日夜夜操网爽| 久久人人97超碰香蕉20202| 国产一区二区激情短视频| 少妇的丰满在线观看| 一本大道久久a久久精品| 少妇精品久久久久久久| 免费观看av网站的网址| 亚洲国产中文字幕在线视频| 成人特级黄色片久久久久久久 | 国产伦理片在线播放av一区| 一区在线观看完整版| 在线观看舔阴道视频| 在线观看一区二区三区激情| 中国美女看黄片| 啦啦啦在线免费观看视频4| 极品少妇高潮喷水抽搐| 亚洲黑人精品在线| 精品少妇久久久久久888优播| 久久精品国产a三级三级三级| 最新的欧美精品一区二区| 老司机深夜福利视频在线观看| 午夜福利免费观看在线| 国产精品一区二区在线不卡| 国产精品久久久人人做人人爽| 中文字幕色久视频| 日韩免费av在线播放| 久久精品亚洲av国产电影网| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久人妻精品电影 | 国产成人啪精品午夜网站| 精品人妻1区二区| 99精品久久久久人妻精品| 国产精品久久久久久精品电影小说| 大陆偷拍与自拍| 怎么达到女性高潮| 午夜福利,免费看| 欧美国产精品一级二级三级| 波多野结衣av一区二区av| 一区二区av电影网| 老熟妇仑乱视频hdxx| 日韩大码丰满熟妇| 又紧又爽又黄一区二区| 一级,二级,三级黄色视频| 757午夜福利合集在线观看| 精品国产一区二区三区久久久樱花| 窝窝影院91人妻| 黄片小视频在线播放| 大片电影免费在线观看免费| 丁香六月天网| 亚洲人成伊人成综合网2020| 欧美 日韩 精品 国产| 男女之事视频高清在线观看| 成人永久免费在线观看视频 | 色播在线永久视频| 啦啦啦 在线观看视频| 精品第一国产精品| 日韩成人在线观看一区二区三区| 久久久久视频综合| 国产日韩欧美亚洲二区| 99热网站在线观看| 国产一区二区三区视频了| 汤姆久久久久久久影院中文字幕| 国产精品免费视频内射| 51午夜福利影视在线观看| 久久人妻熟女aⅴ| 黄片播放在线免费| 天堂动漫精品| 国产成人av激情在线播放| 1024香蕉在线观看| 国产一区二区 视频在线| 日韩欧美一区视频在线观看| 国产成人啪精品午夜网站| 国产真人三级小视频在线观看| 首页视频小说图片口味搜索| 丝袜喷水一区| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 亚洲熟女精品中文字幕| 国产日韩欧美视频二区| 桃红色精品国产亚洲av| 亚洲中文字幕日韩| 高清黄色对白视频在线免费看| 久久精品国产综合久久久| 国产三级黄色录像| 黄色毛片三级朝国网站| 成人免费观看视频高清| 国产成人免费无遮挡视频| 岛国毛片在线播放| 无限看片的www在线观看| 国产男靠女视频免费网站| 久久性视频一级片| 成人国产一区最新在线观看| 久久精品国产99精品国产亚洲性色 | 国产欧美日韩一区二区三| 免费在线观看日本一区| 成在线人永久免费视频| 成年版毛片免费区| 波多野结衣一区麻豆| 日韩制服丝袜自拍偷拍| 大片电影免费在线观看免费| 男女床上黄色一级片免费看| 咕卡用的链子| 亚洲精品自拍成人| 91字幕亚洲| av又黄又爽大尺度在线免费看| 午夜福利乱码中文字幕| 啦啦啦中文免费视频观看日本| 久久中文字幕一级| 亚洲欧美日韩另类电影网站| 黄色片一级片一级黄色片| 丝袜喷水一区| 精品第一国产精品| 大陆偷拍与自拍| 午夜激情av网站| 青草久久国产| 91成人精品电影| 国产成人啪精品午夜网站| 99国产精品免费福利视频| 91九色精品人成在线观看| 欧美精品人与动牲交sv欧美| 欧美日韩黄片免| 亚洲国产av新网站| 80岁老熟妇乱子伦牲交| 国产成人免费无遮挡视频| 久久精品国产亚洲av高清一级| 侵犯人妻中文字幕一二三四区| 一区在线观看完整版| 免费黄频网站在线观看国产| 最新美女视频免费是黄的| 91成年电影在线观看| 下体分泌物呈黄色| 亚洲人成电影免费在线| 丁香欧美五月| 午夜免费鲁丝| 51午夜福利影视在线观看| 一区二区三区激情视频| 桃花免费在线播放| 久久中文看片网| 一区二区三区国产精品乱码| 母亲3免费完整高清在线观看| 免费在线观看影片大全网站| 777米奇影视久久| 99精品欧美一区二区三区四区| 一区二区三区乱码不卡18| 久久亚洲真实| 大型黄色视频在线免费观看| 法律面前人人平等表现在哪些方面| 成年人免费黄色播放视频| 超碰成人久久| 男人操女人黄网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久影院123| 一个人免费在线观看的高清视频| 久热爱精品视频在线9| av福利片在线| 在线观看免费视频网站a站| 啦啦啦在线免费观看视频4| 国产成人精品无人区| 亚洲国产欧美日韩在线播放| 色综合欧美亚洲国产小说| 高清视频免费观看一区二区| 免费在线观看视频国产中文字幕亚洲| a在线观看视频网站| 国产男靠女视频免费网站| 亚洲欧洲精品一区二区精品久久久| 亚洲中文av在线| 精品福利观看| 欧美日韩一级在线毛片| 精品视频人人做人人爽| 国产免费视频播放在线视频| 亚洲国产欧美一区二区综合| 丝瓜视频免费看黄片| 在线播放国产精品三级| 国产免费现黄频在线看| 水蜜桃什么品种好| 国产麻豆69| 另类亚洲欧美激情| netflix在线观看网站| 汤姆久久久久久久影院中文字幕| videos熟女内射| 亚洲情色 制服丝袜| 国产精品av久久久久免费| 欧美成人午夜精品| 在线av久久热| 国产精品二区激情视频| www日本在线高清视频| 99精品久久久久人妻精品| 在线观看免费视频网站a站| 在线观看免费高清a一片| av网站免费在线观看视频| 嫩草影视91久久| 欧美国产精品一级二级三级| 久久国产精品大桥未久av| 成人特级黄色片久久久久久久 | 亚洲 国产 在线| 男女高潮啪啪啪动态图| 久久中文字幕人妻熟女| 欧美日韩黄片免| 一级片'在线观看视频| 女人久久www免费人成看片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产av影院在线观看| 日韩欧美一区二区三区在线观看 | 侵犯人妻中文字幕一二三四区| 99久久精品国产亚洲精品| 午夜福利,免费看| 欧美精品啪啪一区二区三区| 国产极品粉嫩免费观看在线| 在线永久观看黄色视频| 丰满少妇做爰视频| 一区二区三区精品91| tube8黄色片| 日日爽夜夜爽网站| 亚洲av国产av综合av卡| 少妇精品久久久久久久| 窝窝影院91人妻| 1024视频免费在线观看| 无遮挡黄片免费观看| 丰满迷人的少妇在线观看| 免费观看人在逋| 国产高清视频在线播放一区| 久久久国产一区二区| 一个人免费在线观看的高清视频| 日韩中文字幕欧美一区二区| 成人三级做爰电影| 国产成人精品在线电影| 亚洲成人免费电影在线观看| 欧美成人免费av一区二区三区 | 高清黄色对白视频在线免费看| 91字幕亚洲| 久久热在线av| 欧美在线一区亚洲| 十八禁高潮呻吟视频| 精品少妇久久久久久888优播| 国产精品 国内视频| 新久久久久国产一级毛片| 精品国内亚洲2022精品成人 | 三级毛片av免费| 亚洲男人天堂网一区| 老司机在亚洲福利影院| 极品人妻少妇av视频| 一二三四社区在线视频社区8| 亚洲欧美精品综合一区二区三区| 成人国语在线视频|