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

    2003~2018年干旱對中國森林綠度異常的影響

    2022-01-21 02:11:44趙洪飛他志杰
    中國環(huán)境科學(xué) 2022年1期
    關(guān)鍵詞:綠度異?,F(xiàn)象針葉林

    侯 鑫,趙 杰,趙洪飛,他志杰,岳 超,3*

    2003~2018年干旱對中國森林綠度異常的影響

    侯 鑫1,趙 杰1,趙洪飛1,他志杰2,岳 超1,3*

    (1.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院,陜西 楊凌 712100;2.西安外國語大學(xué)旅游學(xué)院·人文地理研究所,陜西 西安 710128;3.中國科學(xué)院水利部水土保持研究所,陜西 楊凌 712100)

    基于衛(wèi)星觀測和氣象再分析數(shù)據(jù)提取2003~2018年間中國森林綠度異?,F(xiàn)象(月尺度LAI長期趨勢的累積偏差),并分析干旱對森林綠度異常影響的時(shí)空動態(tài)變化,探究不同森林類型的干旱敏感性.結(jié)果表明:①干旱導(dǎo)致的綠度異?,F(xiàn)象頻率分布呈現(xiàn)出南高北低的空間格局,而高強(qiáng)度異?,F(xiàn)象主要分布于中國的東北與西南地區(qū).②2003~2018年間綠度異常嚴(yán)重程度呈現(xiàn)顯著增加趨勢(即綠度異常強(qiáng)度指數(shù)下降),變化率為:-0.06/a(<0.05),而綠度異常面積比率下降,變化率為-0.0049/a.③中國森林對干旱脅迫較為敏感的區(qū)域主要分布于云南北部與大興安嶺北部,這是由區(qū)域的地質(zhì)和氣候特點(diǎn)以及植被類型共同作用導(dǎo)致的.④干旱主導(dǎo)區(qū)不同森林類型的干旱敏感性依次為:闊葉林>針葉林>混交林,綠度異常與SPEI回歸斜率依次為1.8>1.3>1.2.⑤針葉林遭受的干旱強(qiáng)度最高,2003~2018年間統(tǒng)計(jì)的SPEI中位數(shù)為-1.65,綠度異?,F(xiàn)象也最為嚴(yán)重,多年綠度異常指數(shù)中位數(shù)為-1.81,說明干旱強(qiáng)度而非敏感性主導(dǎo)了森林綠度異常的現(xiàn)象.

    中國;綠度異常;干旱;森林;敏感性

    干旱是影響森林綠度異常主要的脅迫因子之一[1],其通過影響植物光合作用與呼吸作用對森林的生長帶來負(fù)面影響[2-3].持續(xù)性嚴(yán)重干旱甚至?xí)斐缮执竺娣e死亡[4].

    不同森林類型對干旱的敏感性存在差異.Liu等[5]基于不同生物群落中極端低降水頻率與植被綠度異常頻率的回歸斜率定義敏感性,發(fā)現(xiàn)闊葉林的干旱敏感性高于針葉林.Li等[6]通過樹木年輪和植物性狀數(shù)據(jù),表明裸子植物具有更高的干旱抵抗力.與之相反的是,Huang等[7]的研究表明闊葉林在干旱期間保持原有增強(qiáng)型植被指數(shù)(EVI)水平的能力比針葉林更強(qiáng).以上結(jié)果的差異說明有必要在全國范圍內(nèi)開展不同森林類型的干旱敏感性研究.

    中國是干旱災(zāi)害發(fā)生頻率最高且影響最嚴(yán)重的國家之一,近年我國南方地區(qū)干旱頻次和強(qiáng)度明顯增加[8-9],干旱已經(jīng)造成了我國森林的大面積死亡[10].國內(nèi)有關(guān)干旱脅迫對森林影響的研究多位于干旱與半干旱區(qū)域[11-13].由于干旱半干旱區(qū)的氣候條件限制了森林的最大蓋度[14],在大尺度研究中可被遙感產(chǎn)品識別的森林多位于濕潤區(qū).目前星載遙感技術(shù)已經(jīng)被廣泛應(yīng)用于研究森林?jǐn)_動[15-17],因此本研究采用一種新的方法通過2003~2018年間中分辨率成像光譜儀(MODIS)葉面積指數(shù)刻畫干旱影響下中國森林綠度異常的時(shí)空動態(tài),量化干旱對中國森林的影響,并更加深入的探究不同森林類型干旱敏感性差異,提高對于干旱引起的森林生長異常甚至死亡的風(fēng)險(xiǎn)認(rèn)識;為森林管理者對森林實(shí)施區(qū)域化管理,制定合適的策略以提高森林生態(tài)系統(tǒng)韌性提供理論參考.

    1 研究方法與數(shù)據(jù)來源

    1.1 數(shù)據(jù)來源與預(yù)處理

    葉面積指數(shù)(LAI)是核算生態(tài)系統(tǒng)光合作用、呼吸作用等過程的關(guān)鍵變量[18].作為全球植被的基本屬性,LAI已經(jīng)被全球氣候變化研究界列為必需氣候變量[19].McDowell等[20]的研究表明干旱時(shí)高溫缺水的環(huán)境條件使得飽和水汽壓差增加,致使樹木葉面積減少以適應(yīng)環(huán)境.Liu等[21]的研究也指出干旱導(dǎo)致了樹木莖葉生長的減少.因此本研究選取LAI作為反應(yīng)森林冠層狀況的指標(biāo),用于定量描述植被綠度異?,F(xiàn)象.Steltzer等[22]的研究發(fā)現(xiàn)不同植被類型的光合特性會影響歸一化植被指數(shù)(NDVI)對葉面積的估算.本文數(shù)據(jù)選用美國國家航空航天局(NASA)發(fā)布的MODIS的葉面積指數(shù)遙感產(chǎn)品MCD15A2H(https://ladsweb.modaps.eosdis.nasa.gov/).產(chǎn)品算法選用美國宇航局的Terra和Aqua衛(wèi)星在8d時(shí)段中的最佳影像數(shù)據(jù).時(shí)間跨度為2003年1月~2018年12月,空間分辨率為500m,時(shí)間間隔為8d,共計(jì)736期影像.通過求均值將LAI的空間分辨率聚合到0.1°,保留原時(shí)間分辨率不變.

    土地覆蓋類型數(shù)據(jù)使用第6版MODIS土地覆蓋類型(MCD12Q1)數(shù)據(jù)產(chǎn)品.該產(chǎn)品提供2003~ 2018年空間分辨率為500m的年尺度土地覆蓋類型數(shù)據(jù).該產(chǎn)品給出了基于6種不同土地覆蓋分類方案的數(shù)據(jù),本研究采用國際地球生物圈計(jì)劃(IGBP)分類方案.選定土地覆蓋類型中的1~5(常綠針葉林、常綠闊葉林、落葉針葉林、落葉闊葉林、混交林)作為森林.通過眾數(shù)法將空間分辨率轉(zhuǎn)為0.1°,獲得0.1°的中國森林類型分布.

    氣象數(shù)據(jù)采用中國區(qū)域地面氣象要素驅(qū)動數(shù)據(jù)集(CMFD),該數(shù)據(jù)是專門為研究中國陸地表層過程而開發(fā)的高時(shí)空分辨率近地表氣象數(shù)據(jù)集[23].本研究采用的溫度和降水?dāng)?shù)據(jù)時(shí)間跨度為1979~2018年,空間分辨率為0.1°,時(shí)間分辨率為1d.

    利用CMFD氣象數(shù)據(jù)集中的降水和溫度估算了1979~2018年的標(biāo)準(zhǔn)化降水蒸散發(fā)指數(shù)(SPEI),選取其中的2003~2018年用作研究.SPEI是利用降水和潛在蒸散發(fā)(PET)之間的月差值來量化水分盈余或虧缺的程度,其中月尺度PET按照Thornthwaite[24]的方法估算獲得,然后將月差值最終歸一化為對數(shù)邏輯概率分布,得到SPEI序列[25].本研究中,SPEI的時(shí)間尺度選取為12個(gè)月,代表了數(shù)據(jù)所在月前1a截止至該月的累計(jì)水量平衡.長時(shí)間尺度(12個(gè)月)的SPEI表征更多的是水文干旱,能夠反映深層土壤水分狀態(tài)變化[26].森林具有較深的根系,可以汲取深層土壤水,使其更傾向于對長時(shí)間尺度的干旱產(chǎn)生響應(yīng),因此采用12月時(shí)間尺度的SPEI是適當(dāng)?shù)腫7].

    本研究將不同尺度的數(shù)據(jù)通過GIS實(shí)現(xiàn)數(shù)據(jù)的投影轉(zhuǎn)換與重采樣,使得數(shù)據(jù)空間分辨率統(tǒng)一為0.1°,投影統(tǒng)一為WGS84.通過將土地覆蓋類型數(shù)據(jù)與LAI數(shù)據(jù)疊加,提取中國森林覆蓋區(qū)域的LAI數(shù)據(jù)集.

    1.2 研究方法

    1.2.1 森林綠度異常數(shù)據(jù)的提取 綠度異常指的是在極端氣候影響下,植被葉面積指數(shù)對于其長期生長趨勢的負(fù)向偏離.由于人類土地利用和森林管理、氣候變化、二氧化碳施肥以及氮沉降等諸多因素的作用,2000~2017年間中國森林LAI總體呈現(xiàn)增長趨勢[27].因此本文采用去除線性趨勢法,對LAI數(shù)據(jù)的每個(gè)步長(8d)去趨勢.該方法除了去除LAI的時(shí)間增長趨勢以外,也一并去除了森林生長的季節(jié)周期[28-29],便于觀察年際間植被LAI的異常表現(xiàn).對去趨勢后LAI進(jìn)一步統(tǒng)計(jì)獲得最終綠度異常指數(shù) (GAI),便于量化分析中國森林綠度異常的時(shí)空動態(tài).

    圖1 中國某0.1°像元2005年去除線性趨勢的值

    1.2.2 干旱影響下綠度異常時(shí)空分異特征分析 森林的綠度異常現(xiàn)象可能由多種極端氣候?qū)е?為獲取僅發(fā)生于干旱情況下的綠度異常,采用與Schwalm類似的方法,將SPEI<-1定義為像元出現(xiàn)中度或重度干旱[30],結(jié)合SPEI與綠度異常指數(shù),提取干旱情況下的森林綠度異?,F(xiàn)象.在此基礎(chǔ)上,統(tǒng)計(jì)了2003~2018年間綠度異常的頻率和強(qiáng)度,探究干旱導(dǎo)致綠度異常的頻率和強(qiáng)度的空間分布特征

    式中:表示多年間平均綠度異常強(qiáng)度,表示2003~ 2018年間綠度異?,F(xiàn)象出現(xiàn)的年份次數(shù)(即頻率).

    為了探究干旱情況下發(fā)生森林綠度異常的面積、強(qiáng)度的時(shí)間變化趨勢,累加獲取全國干旱情況下逐年綠度異常面積比率(綠度異常面積與森林總面積的比值),并求取每年全國綠度異常強(qiáng)度的均值,采用線性回歸分析(最小二乘法擬合)獲取隨時(shí)間的變化趨勢.

    2 結(jié)果與分析

    2.1 我國森林類型空間分布特征

    如圖2所示,我國森林空間分布差異化明顯,針葉林主要分布于大興安嶺北部以及四川省西南部.混交林主要分布我國中部一帶,在我國云南東北部也存在混交林聚集分布.闊葉林分布最為廣泛,北至我國黑龍江、吉林、遼寧一帶,中至在我國陜西省南部,同時(shí)我國貴州、廣西、廣東、福建、江西幾省主要森林類型均為闊葉林.

    圖2 中國森林類型的空間分布

    基于國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號為GS(2020)4632號的標(biāo)準(zhǔn)地圖制作,底圖無修改,下同

    2.2 干旱影響下我國森林綠度異常時(shí)空分異特征

    森林綠度異?,F(xiàn)象可能由多種極端氣候事件(干旱、臺風(fēng))驅(qū)動,研究僅保留干旱情況下(SPEI<-1)出現(xiàn)的森林綠度異常.2003~2018年間,累計(jì)出現(xiàn)綠度異?,F(xiàn)象的森林面積約有606萬km2,說明干旱是致使樹木生長異常的重要因素.本文結(jié)果成功刻畫了受到廣泛關(guān)注與報(bào)道的我國森林受干旱影響的事件(圖3),包括:2004年我國南方發(fā)生大范圍綠度異常現(xiàn)象,涉及廣東、浙江、福建和江西省[32].全國森林受干旱影響(即出現(xiàn)綠度異常)總面積約63萬km2,平均綠度異常指數(shù)為0.98,森林出現(xiàn)較大面積綠度異?,F(xiàn)象但程度較輕.2010年我國云南發(fā)生了嚴(yán)重的綠度異?,F(xiàn)象[33].綠度異??偯娣e約47萬km2,像元平均綠度異常指數(shù)為-0.72,即干旱程度低于通常嚴(yán)重程度約一個(gè)標(biāo)準(zhǔn)差,受災(zāi)森林LAI顯著降低.2012, 2013年云南地區(qū)又連續(xù)發(fā)生嚴(yán)重的綠度異?,F(xiàn)象[33-34].2018年我國東南部(廣東、福建和湖南)出現(xiàn)綠度異?,F(xiàn)象[35].總體綠度異常面積約47萬km2,平均綠度異常指數(shù)為-0.25.東南地區(qū)出現(xiàn)綠度異常面積較大但強(qiáng)度較低.這些綠度異常事件的成功識別表明本研究中所用綠度異常的刻畫方法具有較強(qiáng)的可靠性.

    圖3 干旱條件下森林綠度異常分布

    2003~2018年間我國干旱情況下綠度異常強(qiáng)度較高(<0)的區(qū)域主要分布在藏東南、云南省大部、四川南部、貴州廣西和湖南三省交界、華北地區(qū)西部以及大興安嶺一帶(圖4a).綠度異常強(qiáng)度較高的森林面積占全部有林區(qū)的31.05%.綠度異常強(qiáng)度總體上呈現(xiàn)出東北與西南較高,其他地區(qū)較低的空間格局.2003~2018年16a間,綠度異常頻率較高(出現(xiàn)次數(shù)>4)的區(qū)域多位于云南西部與東部.此外,西藏東部、四川南部、湖北與重慶交界以及浙江與安徽交界處也有少量綠度異?,F(xiàn)象聚集(圖4b).綠度異常頻率較高的像元面積占總體有林區(qū)的45.71%.空間分布呈現(xiàn)出南高北低的格局并且高頻像元集中于我國西南一帶.這表明我國南方尤其是西南地區(qū)森林受干旱影響明顯較為嚴(yán)重,綠度異常頻率與強(qiáng)度均較高.干旱情況下極端綠度異常(I<-2)的像元2003~2018年間出現(xiàn)頻率均低于3次,干旱頻次7次以上的像元未出現(xiàn)過極端綠度異常現(xiàn)象,這意味著高頻率與高強(qiáng)度綠度異常事件往往不在同一區(qū)域同時(shí)發(fā)生.

    圖4 干旱條件下森林綠度異常強(qiáng)度與頻率分布

    圖5 干旱條件下森林綠度異常的時(shí)間變化趨勢

    2.3 綠度異常受干旱影響的空間格局

    森林綠度異常指數(shù)與SPEI的回歸擬合優(yōu)度(<0.05)空間差異化明顯.SPEI對于云南省大部、大興安嶺北部以及廣東、福建和浙江南部的綠度異常解釋程度較高(圖6a),表明上述地區(qū)森林冠層葉面積指數(shù)與水分條件的耦合程度更高.回歸斜率則表明綠度異常隨SPEI的變化,即綠度異常對干旱的敏感性.其中云南省北部與大興安嶺北部斜率為正值,表明隨著干旱程度加重綠度異常加重.廣東、福建和浙江南部斜率為負(fù)值,則說明隨著干旱的出現(xiàn)綠度異?,F(xiàn)象反而減輕(圖6b),干旱不是森林綠度異常的主要驅(qū)動因素.結(jié)合兩者判斷,云南省北部和大興安嶺北部和中部的森林綠度異常受到干旱影響較為嚴(yán)重.

    2.4 不同森林類型綠度異常及對干旱敏感性的差異

    圖7 2003~2018年SPEI分布與干旱導(dǎo)致的綠度異常指數(shù)分布統(tǒng)計(jì)

    在分析不同森林類型的干旱敏感性時(shí)將研究區(qū)限制在干旱主導(dǎo)區(qū),即上述綠度異常受干旱影響的空間格局中回歸分析的2>0.4,且回歸斜率為正值的區(qū)域,以排除其他極端氣候事件的干擾.按照不同森林類型劃分干旱主導(dǎo)區(qū)的綠度異常數(shù)據(jù),并通過線性回歸斜率的大小來表征不同森林類型對干旱的敏感性.圖7a中闊葉林的回歸斜率為1.8,針葉林的回歸斜率為1.3,混交林的回歸斜率為1.2.對應(yīng)不同森林類型的敏感性強(qiáng)弱依次為:闊葉林>針葉林>混交林.這表明處于相同的干旱強(qiáng)度時(shí),針葉林的綠度異常現(xiàn)象要比闊葉林輕微,混交林綠度異常表現(xiàn)最小.與干旱敏感性的結(jié)果相反,闊葉林的多年綠度異常指數(shù)中位數(shù)為-1.21,針葉林的多年綠度異常指數(shù)中位數(shù)為-1.81(圖7c).這說明2003~2018年間干旱敏感性較低的針葉林實(shí)際產(chǎn)生的綠度異?,F(xiàn)象比闊葉林更加嚴(yán)重.圖7b的結(jié)果部分解釋了這一現(xiàn)象:闊葉林多年間SPEI值的中位數(shù)為-1.45,針葉林多年間SPEI值的中位數(shù)為-1.65.表明2003~2018年間,針葉林區(qū)域的干旱強(qiáng)度相對闊葉林區(qū)域更高.意味著可能干旱強(qiáng)度主導(dǎo)了森林綠度異常表現(xiàn),而非干旱敏感性.

    3 討論

    為了探究干旱對于森林綠度異常的影響,結(jié)合SPEI與綠度異常指數(shù)獲得我國2003~2018年干旱影響下綠度異常的空間分布.結(jié)果表明干旱影響下綠度異常現(xiàn)象的頻率與強(qiáng)度具有明顯的空間異質(zhì)性.極端綠度異?,F(xiàn)象往往發(fā)生于干旱頻率相對較低的區(qū)域.這可能是由于這些區(qū)域的植被缺少對于極端水分虧缺的應(yīng)對機(jī)制,這種現(xiàn)象同樣有可能發(fā)生于干旱適應(yīng)性強(qiáng)的裸子植物[36].對綠度異常時(shí)間變化趨勢的研究表明,綠度異常的強(qiáng)度顯著增強(qiáng).這可能是由于我國南方極端干旱事件明顯增多,西南地區(qū)特大干旱事件頻發(fā),干旱中心的南移使得中度到重度干旱事件的頻率與強(qiáng)度增加[37].綠度異常面積的下降趨勢可能是由于面積變化存在周期性波動,王強(qiáng)等[17]統(tǒng)計(jì)的農(nóng)作物旱災(zāi)受災(zāi)面積在1951~ 2015年存在總體增大趨勢,在這期間的2000~2015年受災(zāi)面積卻表現(xiàn)出顯著下降趨勢.本文森林干旱影響下綠度異常面積可能與之相似,選取的時(shí)間序列處于大周期的波動下降階段.

    綠度異常的干旱主導(dǎo)區(qū)主要分布于云南省北部和大興安嶺北部和中部.從云南地質(zhì)和氣候特點(diǎn)角度考慮,云南地區(qū)干季持續(xù)時(shí)間長,從11月到次年4月,長達(dá)6個(gè)月[38],干旱事件持續(xù)時(shí)間長,強(qiáng)度大[9].該地區(qū)的巖溶地貌地表水易下滲成為地下水,使得土壤水分難以保持,容易導(dǎo)致土壤水分虧缺[39].從森林類型角度考慮,我國南方多為被子植物,相較于裸子植物對干旱的抵抗能力更弱[6,40].這些條件共同造成了云南北部干旱主導(dǎo)區(qū)的集中現(xiàn)象.大興安嶺中部和北部地區(qū)干季為11月到次年3月,占年降水量的約10%~12%[41].落葉松是大興安嶺森林生態(tài)系統(tǒng)的優(yōu)勢樹種.大興安嶺干旱主導(dǎo)區(qū)單一的生態(tài)系統(tǒng)格局導(dǎo)致物種多樣性較低,減弱了對于干旱的抵抗力[42].雖然研究普遍認(rèn)為裸子植物本身相較于被子植物具有更高的抵抗力,但森林在頻繁或長期的干旱下,仍有可能產(chǎn)生綠度異?,F(xiàn)象[10,6,36,43-44].另一方面,大興安嶺區(qū)域的綠度異常可能來自于火災(zāi)這種通常與干旱并發(fā)的災(zāi)害.火災(zāi)是興安落葉松生長的主要干擾之一.由于干旱導(dǎo)致的氣溫升高、枯落物堆積以及盛行強(qiáng)風(fēng),使得春季很容易發(fā)生森林大火[45].

    我國東南部沿海廣東、福建和浙江南部的SPEI與綠度異常同樣存在高相關(guān)性,但SPEI與綠度異常呈現(xiàn)負(fù)相關(guān).這可能是由于該區(qū)域輻射為森林生長的主要限制因子[10].干旱情況下降水減少,隨著蒸散發(fā)降低,更透明的非降雨云與云量降低使得太陽輻射強(qiáng)度增加,促進(jìn)了植被生長[46]. Morimoto等[47]的研究表明降水是預(yù)測伴隨暴雨的臺風(fēng)災(zāi)害對森林干擾的重要變量,結(jié)合該區(qū)域的綠度異常與水熱條件的高2表現(xiàn)(圖6a)以及該區(qū)域臺風(fēng)災(zāi)害頻發(fā)的特點(diǎn)[48],伴隨暴雨的臺風(fēng)災(zāi)害可能是該區(qū)域森林綠度異常的主要驅(qū)動因素.

    研究通過SPEI與綠度異常數(shù)據(jù)的回歸斜率來表征干旱敏感性,以保證在同一干旱嚴(yán)重程度下展開分析[5].結(jié)果表明針葉林對干旱的敏感性更低,這可能源于針葉林木質(zhì)部水力特征使其在干旱情況下不易產(chǎn)生空化或者栓塞[6].除此之外,其能通過調(diào)節(jié)氣孔開閉控制水分流失[49].該結(jié)果與Huang等[7]的差異可能因?yàn)槠鋵Ω珊得舾行缘亩x未考慮森林是否處于同一干旱強(qiáng)度下.不同森林類型2003~ 2018年間的總體綠度異常表明,針葉林在干旱期間保持原有LAI水平的能力弱于闊葉林,這與兩者的敏感性表現(xiàn)相反.這說明植被LAI對于干旱事件的響應(yīng)可能受到干旱強(qiáng)度、環(huán)境和植被敏感性等因素的復(fù)合影響[7].圖7b也印證了多年間針葉林遭受的干旱程度更高.在更高強(qiáng)度的干旱下,抵抗力更強(qiáng)的針葉林表現(xiàn)出比闊葉林更嚴(yán)重的異?,F(xiàn)象.

    本研究還存在一些不足之處.首先,限于目前氣象數(shù)據(jù)的空間分辨率較低以及遙感數(shù)據(jù)的精度限制無法在大尺度上開展更為精細(xì)的相關(guān)研究.其次,本研究僅通過SPEI所表征的水量平衡來定義干旱,具有一定局限性.因?yàn)榫G度異常可能是由任何與干旱同時(shí)發(fā)生的次生災(zāi)害(如火災(zāi)、蟲災(zāi))造成的[50].更為復(fù)雜的干旱影響機(jī)制需要進(jìn)一步細(xì)化干旱情況下綠度異常的驅(qū)動因素.因此,今后有必要結(jié)合其他森林干擾分布數(shù)據(jù)厘清干旱對森林綠度異常的影響及其驅(qū)動機(jī)制.

    4 結(jié)論

    4.1 干旱情況下中國東北與西南地區(qū)綠度異常強(qiáng)度較高,綠度異常頻率空間分布呈現(xiàn)出南高北低的格局.

    4.2 干旱導(dǎo)致的綠度異常嚴(yán)重程度呈顯著上升趨勢,變化率為-0.06/a,綠度異常面積比率呈下降趨勢,變化率為-0.0049/a.

    4.3 中國森林綠度異常的干旱主導(dǎo)區(qū)集中于云南省北部和大興安嶺北部.

    4.4 干旱主導(dǎo)區(qū)不同森林類型對于干旱的敏感性為:闊葉林>針葉林>混交林,3種森林類型綠度異常與SPEI回歸斜率依次為1.8>1.3>1.2.2003~2018年間干旱主導(dǎo)區(qū)總體干旱強(qiáng)度為:針葉林>闊葉林>混交林,統(tǒng)計(jì)的SPEI中位數(shù)依次為:-1.65<-1.45<-1.42.針葉林總體干旱強(qiáng)度最高,闊葉林與混交林相當(dāng).多年間由干旱導(dǎo)致的森林綠度異?,F(xiàn)象強(qiáng)度為:針葉林>闊葉林>混交林,綠度異常指數(shù)中位數(shù)依次為:-1.21> -1.34>-1.81.針葉林由干旱導(dǎo)致的綠度異常現(xiàn)象最為嚴(yán)重,闊葉林和混交林表現(xiàn)相當(dāng).受干旱強(qiáng)度影響,總體綠度異常表現(xiàn)與敏感性并不完全一致.

    [1] Seidl R, Thom D, Kautz M, et al. Forest disturbances under climate change [J]. Nature Climate Change, 2017,7(6):395-402.

    [2] Griffin J J, Ranney T G, Pharr D M. Heat and drought influence photosynthesis, water relations, and soluble carbohydrates of two ecotypes of redbud (Cercis canadensis) [J]. Journal of the American Society for Horticultural Science, 2004,129(4):497-502.

    [3] Ciais P, Reichstein M, Viovy N, et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003 [J]. Nature, 2005,437(7058):529-533.

    [4] Young D J, Stevens J T, Earles J M, et al. Long‐term climate and competition explain forest mortality patterns under extreme drought [J]. Ecology letters, 2017,20(1):78-86.

    [5] Liu G, Liu H, Yin Y. Global patterns of NDVI-indicated vegetation extremes and their sensitivity to climate extremes [J]. Environmental Research Letters, 2013,8(2):279-288.

    [6] Li X, Piao S, Wang K, et al. Temporal trade-off between gymnosperm resistance and resilience increases forest sensitivity to extreme drought [J]. Nature Ecology & Evolution, 2020,4(8):1075-1083.

    [7] Huang K, Xia J. High ecosystem stability of evergreen broadleaf forests under severe droughts [J]. Global Change Biology, 2019,25(10): 3494-3503.

    [8] 任國玉,封國林,嚴(yán)中偉.中國極端氣候變化觀測研究回顧與展望 [J]. 氣候與環(huán)境研究, 2010,15(4):337-353.

    Ren G Y, Feng G L, Yan Z W, et al. Progresses in observation studies of climate extremes and changes in mainland China [J]. Climatic and Environmental Research, 2010,15(4):337-353.

    [9] 張 強(qiáng),姚玉璧,李耀輝,等.中國干旱事件成因和變化規(guī)律的研究進(jìn)展與展望 [J]. 氣象學(xué)報(bào), 2020,78(3):500-521.

    Zhang Q, Yao Y B, Li Y H, et al. Progress and prospect on the study of causes and variation regularity of droughts in China [J]. Acta Meteorologica Sinica, 2020,78(3):500-521.

    [10] Allen C D, Macalady A K, Chenchouni H, et al. A global overview of drought and heat-induced tree mortality reveals emerging climate change risks for forests [J]. Forest Ecology and Management, 2010, 259(4):660-684.

    [11] 何志斌,杜 軍,陳龍飛,等.干旱區(qū)山地森林生態(tài)水文研究進(jìn)展 [J]. 地球科學(xué)進(jìn)展, 2016,31(10):1078-1089.

    He Z B, Du J, Chen L F, et al. Review on montane forest eco- hydrology in arid area [J]. Advances in Earth Science, 2016,31(10): 1078-1089.

    [12] 趙鵬武,管立娟,劉兵兵,等.我國半干旱區(qū)東段森林動態(tài)研究現(xiàn)狀及展望 [J]. 世界林業(yè)研究, 2021,34(2):74-79.

    Zhao P W, Guan L J, Liu B B, et al. Current research and prospect of forest dynamics in eastern section of Semi-arid Area in China [J]. World Forestry Research, 2021,34(2):74-79.

    [13] 郭永強(qiáng),王乃江,褚曉升,等.基于Google Earth Engine分析黃土高原植被覆蓋變化及原因 [J]. 中國環(huán)境科學(xué), 2019,39(11):4804-4811.

    Guo Y Q, Wang N J, Chu X S, et al. Analyzing vegetation coverage changes and its reasons on the Loess Plateau based on Google Earth Engine. China Environmental Science, 2019,39(11):4804-4811.

    [14] Gunin P D, Vostokova E A, Dorofeyuk N I, et al. Vegetation dynamics of Mongolia [M]. Springer Science & Business Media, 2013:79-129.

    [15] Chen J, Sun L. Using MODIS EVI to detect vegetation damage caused by the 2008 ice and snow storms in south China [J]. Journal of Geophysical Research: Biogeosciences, 2010,115(G3):1-12.

    [16] Hansen M C, Potapov P V, Moore R, et al. High-resolution global maps of 21st-century forest cover change [J]. Science, 2013,342 (6160):850-853.

    [17] 王 強(qiáng),張 勃,戴聲佩,等.三北防護(hù)林工程區(qū)植被覆蓋變化與影響因子分析 [J]. 中國環(huán)境科學(xué), 2012,32(7):1302-1308.

    Wang Q, Zhang B, Dai S P, et al. Analysis of the vegetation cover chang and its relationship with factors in the Three-North Shelter Forest Program [J]. China Environmental Science, 2012,32(7):1302- 1308.

    [18] Fang H, Baret F, Plummer S, et al. An overview of global leaf area index (LAI): Methods, products, validation, and applications [J]. Reviews of Geophysics, 2019,57(3):739-799.

    [19] Mason P, Zillman J, Simmons A, et al. Implementation plan for the global observing system for climate in support of the UNFCCC (2010Update) [Z]. GCOS Rep. 2009:28-29.

    [20] Mcdowell N G, Allen C D. Darcy's law predicts widespread forest mortality under climate warming [J]. Nature Climate Change, 2015, 5(7):669-672.

    [21] Liu J F, Arend M, Yang W J, et al. Effects of drought on leaf carbon source and growth of European beech are modulated by soil type [J]. Scientific Reports, 2017,7(1):1-9.

    [22] Steltzer H, Welker J M. Modeling the effect of photosynthetic vegetation properties on the NDVI–LAI relationship [J]. Ecology, 2006,87(11):2765-2772.

    [23] Yang K, He J, Tang W J, et al. On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan Plateau [J]. Agricultural and Forest Meteorology, 2010,150(1): 38-46.

    [24] Thornthwaite C W. An approach toward a rational classification of climate [J]. Geographical Review, 1948,38(1):55-94.

    [25] Vicente-Serrano S M, Beguería S, López-Moreno J I. A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index [J]. Journal of Climate, 2010, 23(7):1696-1718.

    [26] Wang H, Vicente-Serrano S M, Tao F, et al. Monitoring winter wheat drought threat in Northern China using multiple climate-based drought indices and soil moisture during 2000~2013 [J]. Agricultural and Forest Meteorology, 2016,228:1-12.

    [27] Chen C, Park T, Wang X, et al. China and India lead in greening of the world through land-use management [J]. Nature Sustainability, 2019,2(2):122-129.

    [28] Zhang Y, Liu C, Tang Y, et al. Trends in pan evaporation and reference and actual evapotranspiration across the Tibetan Plateau [J]. Journal of Geophysical Research: Atmospheres, 2007,112:D11210.

    [29] Cuo L, Zhang Y, Gao Y, et al. The impacts of climate change and land cover/use transition on the hydrology in the upper Yellow River Basin, China [J]. Journal of Hydrology, 2013,502:37-52.

    [30] Schwalm C R, Anderegg W R, Michalak A M, et al. Global patterns of drought recovery [J]. Nature, 2017,548(7666):202-205.

    [31] 趙 杰,劉雪佳,杜自強(qiáng),等.晝夜增溫速率的不對稱性對新疆地區(qū)植被動態(tài)的影響 [J]. 中國環(huán)境科學(xué), 2017,37(6):2316-2321.

    Zhao J, Liu X J, Du Z Q, et al. Effects of the asymmetric diurnal- warming on vegetation dynamics in Xinjiang [J]. China Environmental Science, 2017,37(6):2316-2321.

    [32] 王維國,饒曉琴,康志明,等.2004年我國重大災(zāi)害性天氣和大尺度環(huán)流特征 [J]. 氣象, 2005,31(5):26-31.

    Wang W G, Rao X Q, Kang Z M, et al. Analysis of the calamitous weather and large-scale circulation in 2004 [J]. Meteorological Monthly, 2005,31(5):26-31.

    [33] 黃榮輝,劉 永,王 林,等.2009年秋至2010年春我國西南地區(qū)嚴(yán)重干旱的成因分析 [J]. 大氣科學(xué), 2012,36(3):443-457.

    Huang R H, Liu Y, Wang L, et al. Analyses of the causes of severe drought occurring in Southwest China from the fall of 2009 to the Spring of 2010 [J]. Chinese Journal of Atmospheric Sciences, 2012, 36(3):443-457.

    [34] 段海霞,王素萍,馮建英.2013年全國干旱狀況及其影響與成因 [J]. 干旱氣象, 2014,32(2):310-316.

    Duan H X, Wang S P, Feng J Y. Drought events and its influence in 2013 in China [J]. Journal of Arid Meteorology, 2014,32(2):310-316.

    [35] 張 宇,王芝蘭,沙 莎,等.2018年夏季全國干旱狀況及其成因 [J]. 干旱氣象, 2018,36(5):884-892.

    Zhang Y, Wang Z L, Sha S, et al. Drought events and its causes in Summer of 2018 in China [J]. Journal of Arid Meteorology, 2018, 36(5):884-892.

    [36] Zhang X, Li X, Manzanedo R D, et al. High risk of growth cessation of planted larch under extreme drought [J]. Environmental Research Letters, 2021,16(1):014040.

    [37] 韓蘭英,張 強(qiáng),賈建英,等.氣候變暖背景下中國干旱強(qiáng)度、頻次和持續(xù)時(shí)間及其南北差異性 [J]. 中國沙漠, 2019,39(5):1-10.

    Han L Y, Zhang Q, Jia J Y, et al. Drought severity, frequency, duration and regional differences in China [J]. Journal of Desert Research, 2019,39(5):1-10.

    [38] 黃中艷.1961~2007年云南干季干濕氣候變化研究 [J]. 氣候變化研究進(jìn)展, 2010,6(2):113-118.

    Huang Z Y. Dry-wet climate change of the dry season in Yunnan During 1961~2007 [J]. Climate Change Research, 2010,6(2):113-118.

    [39] Liu M, Xu X, Wang D, et al. Karst catchments exhibited higher degradation stress from climate change than the non-karst catchments in southwest China: An ecohydrological perspective [J]. Journal of Hydrology, 2016,535:173-180.

    [40] Lü L, Cai H, Yang Y, et al. Geographic patterns and environmental determinants of gymnosperm species diversity in China [J]. Biodiversity Science, 2018,26(11):1133-1146.

    [41] Jiang Y, Zhang J, Han S, et al. Radial growth response of larix gmelinii to climate along a latitudinal gradient in the Greater Khingan Mountains, Northeastern China [J]. Forests, 2016,7(12):295.

    [42] Grossiord C, Granier A, Ratcliffe S, et al. Tree diversity does not always improve resistance of forest ecosystems to drought [J]. Proceedings of the National Academy of Sciences, 2014,111(41): 14812-14815.

    [43] Trugman A T, Detto M, Bartlett M K, et al. Tree carbon allocation explains forest drought-kill and recovery patterns [J]. Ecology Letters, 2018,21(10):1552-1560.

    [44] Li X, Li Y, Chen A, et al. The impact of the 2009/2010 drought on vegetation growth and terrestrial carbon balance in Southwest China [J]. Agricultural and Forest Meteorology, 2019,269:239-248.

    [45] Fan Q, Wang C, Zhang D, et al. Environmental influences on forest fire regime in the Greater Hinggan Mountains, Northeast China [J]. Forests, 2017,8(10):372.

    [46] Josefsson W, Landelius T. Effect of clouds on UV irradiance: As estimated from cloud amount, cloud type, precipitation, global radiation and sunshine duration [J]. Journal of Geophysical Research: Atmospheres, 2000,105(D4):4927-4935.

    [47] Morimoto J, Aiba M, Furukawa F, et al. Risk assessment of forest disturbance by typhoons with heavy precipitation in northern Japan [J]. Forest Ecology and Management, 2021,479:118521.

    [48] 殷 潔,戴爾阜,吳紹洪.中國臺風(fēng)災(zāi)害綜合風(fēng)險(xiǎn)評估與區(qū)劃 [J]. 地理科學(xué), 2013,33(11):1370-1376.

    Yin J, Dai E F, Wu S H. Integrated risk assessment and zoning of typhoon disasters in China [J]. Scientia Geographica Sinica, 2013, 33(11):1370-1376.

    [49] Brodribb T J, Mcadam S A M, Jordan G J, et al. Conifer species adapt to low-rainfall climates by following one of two divergent pathways [J]. Proceedings of the National Academy of Sciences, 2014,111(40): 14489-14493.

    [50] 馬明衛(wèi),韓宇平,嚴(yán)登華,等.特大干旱事件災(zāi)害孕育機(jī)理及影響研究進(jìn)展 [J]. 水資源保護(hù), 2020,36(5):11-21.

    Ma M W, Han Y P, Yan D H, et al. Research progress on the mechanism and influence of extreme drought-induced disasters [J]. Water Resources Protection, 2020,36(5):11-21.

    Drought effects on forest greenness anomalies in China from 2003 to 2018.

    HOU Xin1, ZHAO Jie1, ZHAO Hong-fei1, TA Zhi-jie2, YUE Chao1,3*

    (1.College of Natural Resources and Environment, North West Agriculture and Forestry University, Yangling 712100, China;2.School of Tourism & Research Institute of Human Geography, Xi'an International Studies University, Xi'an 710128, China;3.Institute of Soil and Water Conservation, CAS & MWR, Yangling 712100, China)., 2022,42(1):336~344

    This study was to evaluate the effects of drought on forest greenness anomalies in China. Satellite observations of forest leaf area index (LAI) and climate reanalysis data were used to examine forest canopy damage induced by drought during the period from 2003 to 2018. We assigned the forest Greenness Anomaly index (GAI) as a proxy for the forest canopy damage, defined as cumulative monthly negative LAI anomalies deviated from the long-term monthly LAI trend, then followed by standardization using the standard deviation. Using this index, we analyzed the spatiotemporal dynamics of drought-induced forest damage and the drought sensitivity of different forest types. Our results showed that: (1) The drought-induced forest canopy damage was more widespread in South China than in North China, while high-intensity drought damage mainly occurred in the Northeast and Southwest of China; (2) From 2003 to 2018, the severity of drought-induced canopy damage increased significantly with the GAI change rate of -0.06/a (a lower index means a higher severity), while the forest area ratio subjected to the drought damage decreased with the change rate of -0.0049/a; (3) Forests with a high sensitivity of drought damage mainly occurred in the north of Yunnan and the north part of the Greater Khingan Mountains in China, due to the common effects of geological, climate and vegetation type characteristics of the region; (4) Broadleaf forest was most sensitive to drought, followed by coniferous forest, while mixed forest was the least sensitive forest type and the corresponding regression slopes of GAI and SPEI for these three forest types were 1.8, 1.3 and 1.2, respectively; and (5) from 2003 to 2018, coniferous forest suffered the highest drought intensity and forest damage, with the median of SPEI as -1.65. As a result forest damage was more dominated by the drought intensity than by drought sensitivity.

    China;greenness anomaly;drought;forest;sensitivity

    X173

    A

    1000-6923(2022)01-0336-09

    侯 鑫(1995-),男,河南新鄉(xiāng)人,西北農(nóng)林科技大學(xué)碩士研究生,主要從事氣候變化與植被響應(yīng)方向研究.

    2021-06-08

    陜西省自然科學(xué)基礎(chǔ)研究計(jì)劃(2021JQ-771)

    * 責(zé)任作者, 研究員, chaoyue@ms.iswc.ac.cn

    猜你喜歡
    綠度異?,F(xiàn)象針葉林
    寧夏草地綠度時(shí)空變化圖譜及驅(qū)動因素分析
    基于GEE的黃土高原植被綠度線推移變化研究
    金沙江上游與其他流域典型植被碳含量差異性探究
    全球森林綠度異常對干旱事件的響應(yīng)
    大興安嶺主要針葉樹種移植容器苗的培育技術(shù)與造林效果的研究
    施肥對油松人工林針葉養(yǎng)分變化的影響
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    淺談唐卡中綠度母的造型特點(diǎn)
    變壓器運(yùn)行中的異?,F(xiàn)象與故障處理
    電子制作(2017年7期)2017-06-05 09:36:14
    前向散射能見度儀的常見異常現(xiàn)象處理及日常維護(hù)
    亚洲欧美精品综合一区二区三区| 岛国视频午夜一区免费看| 精品国产美女av久久久久小说| 午夜福利免费观看在线| 99久久国产精品久久久| 日韩大尺度精品在线看网址| 国产精品久久久久久久电影 | 熟女电影av网| 国产高清激情床上av| av有码第一页| 欧美乱妇无乱码| 性欧美人与动物交配| 精品国产乱子伦一区二区三区| 国语自产精品视频在线第100页| 亚洲乱码一区二区免费版| 中文字幕人妻丝袜一区二区| 午夜免费激情av| 午夜福利成人在线免费观看| 亚洲一区中文字幕在线| av天堂在线播放| 亚洲精品粉嫩美女一区| 国产精品美女特级片免费视频播放器 | 国产精品永久免费网站| 婷婷亚洲欧美| 巨乳人妻的诱惑在线观看| 亚洲av日韩精品久久久久久密| 午夜精品在线福利| 中文字幕最新亚洲高清| 日韩中文字幕欧美一区二区| 男女床上黄色一级片免费看| 成人手机av| 色尼玛亚洲综合影院| 国产精品亚洲美女久久久| 日韩精品青青久久久久久| 波多野结衣高清作品| 精品久久久久久久末码| 天堂影院成人在线观看| 在线观看66精品国产| 国产一区在线观看成人免费| 日日干狠狠操夜夜爽| 久久久久性生活片| 黄色毛片三级朝国网站| 97碰自拍视频| 舔av片在线| 精品熟女少妇八av免费久了| 亚洲 欧美一区二区三区| 成人特级黄色片久久久久久久| avwww免费| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色| netflix在线观看网站| 日韩欧美 国产精品| 一a级毛片在线观看| 亚洲欧美一区二区三区黑人| 日本三级黄在线观看| 又黄又爽又免费观看的视频| 亚洲美女视频黄频| 久久精品国产亚洲av高清一级| 又爽又黄无遮挡网站| 波多野结衣巨乳人妻| 黄色 视频免费看| 久久久久久久久久黄片| 午夜免费激情av| 色尼玛亚洲综合影院| 欧美精品亚洲一区二区| 日韩欧美三级三区| 91国产中文字幕| 制服人妻中文乱码| 色哟哟哟哟哟哟| 美女高潮喷水抽搐中文字幕| 国产亚洲av嫩草精品影院| 91老司机精品| 嫩草影视91久久| 亚洲熟妇熟女久久| 日韩有码中文字幕| 国产熟女xx| 最好的美女福利视频网| 免费观看人在逋| 久久草成人影院| 97碰自拍视频| 国产av一区二区精品久久| 69av精品久久久久久| 国产v大片淫在线免费观看| 又黄又粗又硬又大视频| avwww免费| 啪啪无遮挡十八禁网站| 两个人看的免费小视频| 国产高清视频在线播放一区| 91成年电影在线观看| 操出白浆在线播放| 亚洲va日本ⅴa欧美va伊人久久| 欧美av亚洲av综合av国产av| 一级片免费观看大全| 欧美另类亚洲清纯唯美| 久久午夜综合久久蜜桃| 国产成人啪精品午夜网站| 欧美中文日本在线观看视频| 最近最新中文字幕大全电影3| 在线国产一区二区在线| www.精华液| 国产久久久一区二区三区| 蜜桃久久精品国产亚洲av| 亚洲av中文字字幕乱码综合| 成人国产综合亚洲| 18禁观看日本| 久久99热这里只有精品18| 人人妻,人人澡人人爽秒播| 国产真人三级小视频在线观看| 欧美一区二区国产精品久久精品 | 欧美人与性动交α欧美精品济南到| 首页视频小说图片口味搜索| 一级毛片女人18水好多| 午夜日韩欧美国产| 天天一区二区日本电影三级| 老汉色∧v一级毛片| 午夜福利高清视频| 欧美av亚洲av综合av国产av| 精品乱码久久久久久99久播| 久久精品综合一区二区三区| 搡老妇女老女人老熟妇| 一级作爱视频免费观看| 非洲黑人性xxxx精品又粗又长| 精品久久久久久成人av| 国产探花在线观看一区二区| 国产乱人伦免费视频| 成人国产一区最新在线观看| www日本黄色视频网| 91麻豆精品激情在线观看国产| 免费看美女性在线毛片视频| 国模一区二区三区四区视频 | 欧美一区二区国产精品久久精品 | 免费观看人在逋| 两个人视频免费观看高清| 一边摸一边抽搐一进一小说| 禁无遮挡网站| 在线观看免费日韩欧美大片| 性欧美人与动物交配| 欧美高清成人免费视频www| 亚洲 欧美 日韩 在线 免费| 精品国产美女av久久久久小说| 国产一区二区激情短视频| 极品教师在线免费播放| 99在线视频只有这里精品首页| 蜜桃久久精品国产亚洲av| 可以在线观看毛片的网站| 亚洲国产精品久久男人天堂| 久久精品国产综合久久久| 久久国产精品影院| 男女午夜视频在线观看| 99久久精品热视频| 亚洲国产欧洲综合997久久,| АⅤ资源中文在线天堂| 久热爱精品视频在线9| 午夜成年电影在线免费观看| 亚洲欧美精品综合久久99| 在线看三级毛片| 亚洲国产精品成人综合色| 国内揄拍国产精品人妻在线| 欧美+亚洲+日韩+国产| 国产午夜福利久久久久久| 国产精品综合久久久久久久免费| 黑人操中国人逼视频| 国产av一区在线观看免费| 18禁黄网站禁片午夜丰满| 在线看三级毛片| 免费看十八禁软件| 国产黄a三级三级三级人| 国产av不卡久久| 香蕉av资源在线| 美女免费视频网站| 亚洲精品一区av在线观看| 国产精品乱码一区二三区的特点| 国产99久久九九免费精品| √禁漫天堂资源中文www| 美女 人体艺术 gogo| 欧美日韩亚洲国产一区二区在线观看| 成年免费大片在线观看| 国产高清videossex| 两性夫妻黄色片| 欧美zozozo另类| 欧美日韩中文字幕国产精品一区二区三区| 久久天堂一区二区三区四区| 波多野结衣高清无吗| 亚洲人成网站在线播放欧美日韩| 女同久久另类99精品国产91| 国产亚洲av高清不卡| 国产三级在线视频| 精华霜和精华液先用哪个| 他把我摸到了高潮在线观看| 露出奶头的视频| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 国产成人系列免费观看| 日韩大码丰满熟妇| 国产精品一区二区免费欧美| 每晚都被弄得嗷嗷叫到高潮| 久久久国产精品麻豆| 在线观看日韩欧美| 窝窝影院91人妻| 一进一出抽搐gif免费好疼| 日韩欧美在线二视频| 欧美精品啪啪一区二区三区| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 久久久久久久久免费视频了| 亚洲成人久久性| 在线免费观看的www视频| 51午夜福利影视在线观看| 精品久久蜜臀av无| 国产精品98久久久久久宅男小说| 国产主播在线观看一区二区| 成年人黄色毛片网站| 久久久久性生活片| 国产乱人伦免费视频| 国产精品,欧美在线| 好男人在线观看高清免费视频| 久久久久久久久免费视频了| 日本免费一区二区三区高清不卡| av超薄肉色丝袜交足视频| 制服人妻中文乱码| www.自偷自拍.com| 夜夜爽天天搞| 脱女人内裤的视频| 亚洲乱码一区二区免费版| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻 | 琪琪午夜伦伦电影理论片6080| 狂野欧美白嫩少妇大欣赏| 国产av麻豆久久久久久久| 午夜福利在线在线| 90打野战视频偷拍视频| 久久久久国产一级毛片高清牌| 激情在线观看视频在线高清| avwww免费| 午夜激情福利司机影院| 亚洲自偷自拍图片 自拍| 亚洲国产看品久久| 麻豆av在线久日| 九色国产91popny在线| 国产野战对白在线观看| 亚洲国产精品成人综合色| 黑人欧美特级aaaaaa片| 久久久久久九九精品二区国产 | 色av中文字幕| 国产亚洲精品久久久久久毛片| 少妇被粗大的猛进出69影院| 好男人在线观看高清免费视频| 在线国产一区二区在线| 久久精品aⅴ一区二区三区四区| 久久性视频一级片| 超碰成人久久| 国产精品一区二区三区四区久久| 黄色片一级片一级黄色片| 欧美日韩国产亚洲二区| 国产91精品成人一区二区三区| 国产欧美日韩一区二区精品| 国产精品国产高清国产av| 久久99热这里只有精品18| 麻豆国产av国片精品| 欧美精品亚洲一区二区| 久久久久久国产a免费观看| 给我免费播放毛片高清在线观看| 中文字幕av在线有码专区| 丁香欧美五月| 别揉我奶头~嗯~啊~动态视频| 人妻久久中文字幕网| 久久久久久免费高清国产稀缺| 床上黄色一级片| 中出人妻视频一区二区| 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 中出人妻视频一区二区| 久久精品91蜜桃| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区| 日本五十路高清| 一级片免费观看大全| 一级毛片女人18水好多| 日韩 欧美 亚洲 中文字幕| 99热6这里只有精品| 亚洲欧美日韩东京热| 亚洲成人国产一区在线观看| 国产日本99.免费观看| 成人国产综合亚洲| 91成年电影在线观看| 国产成人av教育| 中文亚洲av片在线观看爽| 亚洲中文av在线| 熟女少妇亚洲综合色aaa.| 国产亚洲精品久久久久久毛片| 一个人免费在线观看的高清视频| 国产日本99.免费观看| 小说图片视频综合网站| 国产一区二区在线av高清观看| 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| 18禁美女被吸乳视频| 啦啦啦韩国在线观看视频| 啦啦啦免费观看视频1| 桃红色精品国产亚洲av| 国产亚洲精品av在线| 国产精品香港三级国产av潘金莲| 美女 人体艺术 gogo| 日韩中文字幕欧美一区二区| 国产黄色小视频在线观看| 91大片在线观看| 欧美又色又爽又黄视频| 青草久久国产| 欧美大码av| 亚洲片人在线观看| 国产精品美女特级片免费视频播放器 | 狂野欧美激情性xxxx| 熟女电影av网| 大型av网站在线播放| 精品久久久久久久久久免费视频| 久久国产乱子伦精品免费另类| 午夜激情福利司机影院| 亚洲精品在线美女| 黄频高清免费视频| 丁香六月欧美| 美女黄网站色视频| 国产av一区在线观看免费| 亚洲av成人不卡在线观看播放网| 国产成人影院久久av| 亚洲国产日韩欧美精品在线观看 | tocl精华| 99国产精品一区二区三区| 国产91精品成人一区二区三区| 久久久久性生活片| 两个人的视频大全免费| 免费搜索国产男女视频| 在线十欧美十亚洲十日本专区| 日本精品一区二区三区蜜桃| 欧美极品一区二区三区四区| 日韩精品青青久久久久久| 精品熟女少妇八av免费久了| 精品少妇一区二区三区视频日本电影| 日日干狠狠操夜夜爽| 国内精品久久久久精免费| 黄片小视频在线播放| 首页视频小说图片口味搜索| 久久久久久九九精品二区国产 | 免费无遮挡裸体视频| 琪琪午夜伦伦电影理论片6080| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美精品v在线| 国产成人aa在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产又色又爽无遮挡免费看| 夜夜夜夜夜久久久久| 精品国产亚洲在线| 久久精品国产亚洲av高清一级| 日韩大尺度精品在线看网址| 色噜噜av男人的天堂激情| 免费看日本二区| 成人av一区二区三区在线看| 亚洲精品久久国产高清桃花| 在线十欧美十亚洲十日本专区| 欧美日韩乱码在线| 精品国内亚洲2022精品成人| 精品国产美女av久久久久小说| 最新在线观看一区二区三区| 日本黄大片高清| 日日干狠狠操夜夜爽| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美三级三区| bbb黄色大片| 欧美乱码精品一区二区三区| 十八禁网站免费在线| 欧美日韩一级在线毛片| 婷婷亚洲欧美| 丝袜美腿诱惑在线| 亚洲一区二区三区不卡视频| 变态另类丝袜制服| 制服人妻中文乱码| 国产精品影院久久| 精品一区二区三区四区五区乱码| 亚洲专区字幕在线| 黄片小视频在线播放| 婷婷丁香在线五月| 欧美又色又爽又黄视频| 久久久国产成人精品二区| 黄色 视频免费看| 亚洲精华国产精华精| 搞女人的毛片| videosex国产| 中文字幕av在线有码专区| 国产91精品成人一区二区三区| 99久久综合精品五月天人人| 又粗又爽又猛毛片免费看| 在线观看免费午夜福利视频| 精品国产亚洲在线| 久久久久九九精品影院| 国产伦一二天堂av在线观看| 日韩欧美免费精品| 男男h啪啪无遮挡| 天堂动漫精品| 国产激情偷乱视频一区二区| 免费在线观看完整版高清| 欧洲精品卡2卡3卡4卡5卡区| 日韩 欧美 亚洲 中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲乱码一区二区免费版| 精品熟女少妇八av免费久了| 麻豆久久精品国产亚洲av| 一个人免费在线观看电影 | 99久久综合精品五月天人人| 日本黄大片高清| 欧美人与性动交α欧美精品济南到| 国产亚洲av高清不卡| 亚洲国产精品成人综合色| 精品午夜福利视频在线观看一区| 亚洲精品在线观看二区| 午夜老司机福利片| 久久久久久人人人人人| 热99re8久久精品国产| 悠悠久久av| 国产精华一区二区三区| 香蕉久久夜色| 一个人免费在线观看电影 | 日韩精品中文字幕看吧| 免费看美女性在线毛片视频| 91九色精品人成在线观看| 性色av乱码一区二区三区2| 一区二区三区国产精品乱码| 丰满人妻一区二区三区视频av | 精品国内亚洲2022精品成人| 午夜免费观看网址| 欧美久久黑人一区二区| 香蕉av资源在线| 午夜成年电影在线免费观看| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 亚洲av成人不卡在线观看播放网| 国产成人av教育| 午夜影院日韩av| 免费在线观看亚洲国产| 欧美+亚洲+日韩+国产| 成人av在线播放网站| 18禁观看日本| 亚洲成人久久爱视频| 免费在线观看影片大全网站| 久久中文字幕一级| 亚洲午夜精品一区,二区,三区| 日本a在线网址| 国产亚洲精品综合一区在线观看 | 国产精品永久免费网站| 午夜福利在线在线| 99久久综合精品五月天人人| 欧美在线黄色| 欧美日韩亚洲国产一区二区在线观看| 精华霜和精华液先用哪个| 欧美国产日韩亚洲一区| 色噜噜av男人的天堂激情| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久亚洲av鲁大| 中出人妻视频一区二区| 丁香欧美五月| 成人三级做爰电影| 亚洲片人在线观看| 法律面前人人平等表现在哪些方面| 国产精品久久久av美女十八| 免费看日本二区| 日本黄色视频三级网站网址| 日本一区二区免费在线视频| 久久精品人妻少妇| 无人区码免费观看不卡| 99久久久亚洲精品蜜臀av| 亚洲成a人片在线一区二区| av天堂在线播放| xxxwww97欧美| 日本 av在线| 色av中文字幕| 国产精品av视频在线免费观看| aaaaa片日本免费| 国产精品香港三级国产av潘金莲| 亚洲精品国产一区二区精华液| 色噜噜av男人的天堂激情| 黄色丝袜av网址大全| 深夜精品福利| 怎么达到女性高潮| 淫秽高清视频在线观看| 看片在线看免费视频| 舔av片在线| 在线国产一区二区在线| 视频区欧美日本亚洲| 午夜精品一区二区三区免费看| 国产又色又爽无遮挡免费看| 成人三级做爰电影| 精华霜和精华液先用哪个| www国产在线视频色| 日本黄色视频三级网站网址| 黄色片一级片一级黄色片| 精品电影一区二区在线| av福利片在线| 一区二区三区高清视频在线| 97超级碰碰碰精品色视频在线观看| 男人舔女人的私密视频| 日本一区二区免费在线视频| 国产精品综合久久久久久久免费| 女人被狂操c到高潮| 亚洲美女黄片视频| 成年免费大片在线观看| 亚洲乱码一区二区免费版| 99国产精品一区二区蜜桃av| 亚洲av五月六月丁香网| 免费搜索国产男女视频| 国产精品1区2区在线观看.| 亚洲色图 男人天堂 中文字幕| 夜夜看夜夜爽夜夜摸| 一边摸一边抽搐一进一小说| 欧美成人午夜精品| 母亲3免费完整高清在线观看| 久久久久久久久中文| 日韩欧美一区二区三区在线观看| 亚洲一区高清亚洲精品| 亚洲美女黄片视频| 又爽又黄无遮挡网站| 国产成人系列免费观看| 精品国产亚洲在线| 久久久国产成人免费| 亚洲美女黄片视频| 无人区码免费观看不卡| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 一区福利在线观看| www国产在线视频色| 亚洲成人精品中文字幕电影| 精品欧美一区二区三区在线| 少妇粗大呻吟视频| 免费看日本二区| 美女大奶头视频| 婷婷丁香在线五月| 中文字幕人成人乱码亚洲影| 青草久久国产| 一区二区三区高清视频在线| 无人区码免费观看不卡| 又粗又爽又猛毛片免费看| 亚洲欧美精品综合久久99| 国产一区在线观看成人免费| 欧美极品一区二区三区四区| 他把我摸到了高潮在线观看| 波多野结衣巨乳人妻| or卡值多少钱| 亚洲一区二区三区不卡视频| 久久中文字幕一级| 欧美成人午夜精品| tocl精华| 亚洲全国av大片| 操出白浆在线播放| 午夜福利欧美成人| 法律面前人人平等表现在哪些方面| 黄色女人牲交| 99精品在免费线老司机午夜| 欧美性猛交黑人性爽| 人人妻人人看人人澡| 久久久精品欧美日韩精品| 黄色成人免费大全| 桃红色精品国产亚洲av| 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 欧美av亚洲av综合av国产av| 日本三级黄在线观看| 精品久久久久久,| 欧美久久黑人一区二区| 久久精品综合一区二区三区| 久久精品亚洲精品国产色婷小说| 男插女下体视频免费在线播放| 国产一区二区三区视频了| 国产成人欧美在线观看| 国产精品影院久久| 久久午夜综合久久蜜桃| 九色成人免费人妻av| 久久精品国产综合久久久| 最近视频中文字幕2019在线8| 一本大道久久a久久精品| 岛国在线免费视频观看| 曰老女人黄片| x7x7x7水蜜桃| 男人舔女人下体高潮全视频| 久久 成人 亚洲| 欧美黑人巨大hd| 婷婷精品国产亚洲av在线| 91大片在线观看| 久久香蕉国产精品| 麻豆成人av在线观看| 亚洲成av人片免费观看| 全区人妻精品视频| 亚洲一区高清亚洲精品| 日韩欧美一区二区三区在线观看| 高潮久久久久久久久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 久久天堂一区二区三区四区| 亚洲成人久久爱视频| www.自偷自拍.com| 最近最新中文字幕大全电影3| videosex国产| 色综合婷婷激情| 最近最新中文字幕大全免费视频| 在线观看66精品国产| 麻豆国产av国片精品| 99久久精品热视频| 一级毛片精品| 叶爱在线成人免费视频播放| 悠悠久久av| 少妇的丰满在线观看| 天堂动漫精品| 好看av亚洲va欧美ⅴa在| 亚洲国产欧美一区二区综合| 神马国产精品三级电影在线观看 | 日韩国内少妇激情av| 麻豆av在线久日| 国产亚洲精品av在线| 丝袜美腿诱惑在线| 少妇粗大呻吟视频|