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

    植被對(duì)干旱趨勢(shì)的影響

    2016-03-23 00:34:59劉永強(qiáng)
    大氣科學(xué) 2016年1期
    關(guān)鍵詞:干旱氣候變化植被

    劉永強(qiáng)

    Center for Forest Disturbance Science, USDA Forest Service, 320 Green St., Athens, Georgia 30602, USA

    ?

    植被對(duì)干旱趨勢(shì)的影響

    劉永強(qiáng)

    Center for Forest Disturbance Science, USDA Forest Service, 320 Green St., Athens, Georgia 30602, USA

    摘 要?dú)v史干旱事件的觀測(cè)和數(shù)值研究表明,植被可通過地—?dú)馑?、能量和其他通量交換影響和反饋干旱。本研究旨在了解氣候變化情形下植被對(duì)干旱趨勢(shì)的影響和機(jī)制。應(yīng)用美國(guó)大陸七個(gè)動(dòng)力氣候降尺度區(qū)域氣候變化情景,計(jì)算和分析了現(xiàn)在和未來的干旱指數(shù)、空間分布和季節(jié)變化。通過比較同一氣候區(qū)兩種植被類型區(qū)域干旱強(qiáng)度和頻率理解植被的影響。集成分析結(jié)果表明,未來美國(guó)干旱很可能增加,其中大平原中部所有季節(jié)都很顯著,而東南和西南地區(qū)夏秋更為顯著。植被對(duì)干旱趨勢(shì)的影響和氣候區(qū)有關(guān)。在溫暖和潮濕/干燥氣候區(qū),林地(草地)未來干旱強(qiáng)度和頻率的增幅大于對(duì)應(yīng)的農(nóng)田(荒漠)區(qū)域,因此植被可以放大未來干旱的風(fēng)險(xiǎn)。相反,在寒冷和潮濕氣候區(qū),林地(草地)區(qū)域未來干旱強(qiáng)度和頻率增幅較小,表明植被放大未來干旱的作用可能只在某些氣候情形下出現(xiàn)。這種植被對(duì)未來干旱影響的復(fù)雜性和對(duì)氣候區(qū)的依賴性對(duì)氣候模式提供可靠的干旱模擬和預(yù)測(cè)及森林管理部門制定適應(yīng)和減緩氣候變化的策略提出了新的挑戰(zhàn)。

    關(guān)鍵詞植被 干旱 氣候變化 干旱指數(shù) 美國(guó)氣候區(qū)

    資助項(xiàng)目 美國(guó)自然科學(xué)基金會(huì)、農(nóng)業(yè)部和能源部聯(lián)合項(xiàng)目:地球系統(tǒng)模式預(yù)報(bào)年代際區(qū)域氣候NIFA-2013-35100-20516

    Funded by US NSF, USDA and DOE Joint Program on Decadal and Regional Climate Prediction using Earth System Models (EaSM) NIFA-2013-35100-20516

    Impacts of Vegetation on Drought Trends

    LIU Yongqiang

    Center for Forest Disturbance Science, USDA Forest Service, 320 Green St., Athens, Georgia 30602, USA

    Abstract Vegetation conditions interact with and feed back to atmospheric anomalies such as droughts through water, energy and other exchanges on the land surface. Various mechanisms have been proposed to understand the related physical processes based on observational and numerical analyses of historical droughts. This study aims to understand the complex impacts and mechanisms of vegetation on future droughts under changing climate. Future trends in droughts were projected by comparing a drought index for present and future periods using seven dynamically downscaled regional climate change scenarios for the continental United States. The impacts of vegetation were examined by comparing changes in drought intensity and frequency between two types of land cover in each of six climate regimes. The ensemble results indicate that future droughts would increase in both intensity and frequency, mainly in the central Great Plains during all seasons, and in the Southeast and Southwest U.S. during warm seasons. The roles of vegetation depend on climate regime. Future droughts would increase more significantly on forest lands (grasslands) than the corresponding farmlands (drylands) in warm and moist/dry climate regimes. Vegetation, therefore, can amplify future droughts. In contrast, they would be less significant in cold and moist regimes, suggesting that the capacity of vegetation in amplifying future droughts would exist only in certain climate circumstances. The complex impacts and regional dependence raise new challenges for climate models to provide reliable drought simulation and projection, and for forest management todevelop adaptation and mitigation strategies for climate change.

    Keywords Vegetation, Drought, Climate change, Drought index, U.S. climate regime

    前言

    恩師葉篤正先生于2013年10月溘然長(zhǎng)逝,離開了為之奮斗了七十余載的科學(xué)事業(yè)。作者于上世紀(jì)80年代后期師從葉先生進(jìn)行博士生的學(xué)習(xí)和研究,有幸親耳聆聽了葉先生正直做人和嚴(yán)謹(jǐn)做學(xué)問的教誨、親身感受了對(duì)科學(xué)的摯愛和追求、親眼目睹了對(duì)全球變化研究?jī)A注的巨大心血和作出的卓越貢獻(xiàn)。

    作者感受最深、受益最大的是葉先生的視野、敏銳性和開拓精神。上世紀(jì)70年代,科學(xué)家在探索與人類活動(dòng)息息相關(guān)的重大氣象和環(huán)境問題(如北非70年代初開始的持續(xù)干旱,南美森林砍伐的氣候效應(yīng),我國(guó)北方干旱和沙漠化,全球溫室氣體效應(yīng)等)時(shí)發(fā)現(xiàn),僅僅研究大氣過程已無(wú)法全面、深入地理解其物理機(jī)制和成因。由此提出了“氣候系統(tǒng)”的概念,研究大氣與地球表明其他組成部分之間的相互作用和反饋。葉先生的求學(xué)和早期主要研究領(lǐng)域是大氣動(dòng)力學(xué),但開闊的視野、對(duì)科學(xué)問題的敏銳性和開拓新的研究領(lǐng)域的勇氣,使他成為探索氣候系統(tǒng)的先驅(qū)者。葉先生與GFDL科學(xué)家合作,很快在地—?dú)庀嗷プ饔妙I(lǐng)域做出了開創(chuàng)性工作(Yeh et al., 1984)。此后,葉先生繼續(xù)傾其精力和智慧致力于中國(guó)的全球變化研究。

    作者跟從葉先生進(jìn)行地—?dú)庀嗷プ饔玫难芯?,博士論文?duì)地—?dú)庀到y(tǒng)異常的時(shí)間尺度進(jìn)行了理論分析和數(shù)值模擬,并探討了陸地過程在減緩溫室氣體效應(yīng)中可能的作用。此后二十多年一直繼續(xù)這一領(lǐng)域的研究。本文探討植被對(duì)未來干旱趨勢(shì)的影響,以此紀(jì)念葉先生百年華誕和輝煌人生。

    1 引言

    干旱是一段時(shí)間(季節(jié)或更長(zhǎng))降水明顯低于平均水平,由此導(dǎo)致水資源短缺,并影響植被生長(zhǎng)和人們?nèi)粘I畹默F(xiàn)象(NWS,2006)。干旱可分為氣象、水文和農(nóng)業(yè)干旱。降水偏少是氣象干旱的標(biāo)志,而土壤水分不足是水文和農(nóng)業(yè)干旱的基本特征之一。干旱對(duì)不同時(shí)空尺度、不同地區(qū)的生態(tài)系統(tǒng)結(jié)構(gòu)和過程有著極為不利的影響(例如,Wu and Chen,2013;Vicente-Serrano et al.,2013)。植物主要通過蒸騰過程獲得土壤中的養(yǎng)分;土壤缺水導(dǎo)致植物根系受損,吸水能力下降,同時(shí)氣孔阻力增加,因而蒸騰減少。這意味著光合作用減弱,植物生長(zhǎng)減慢或停止。Zhang et al.(2014)發(fā)現(xiàn),從1986到2006年,氣候變化導(dǎo)致北京干旱增加,樹木死亡率明顯增加。

    植被是氣候系統(tǒng)的組成部分之一。它可以通過其物理和生理過程影響地表能量、水分、動(dòng)量和生物化學(xué)物質(zhì)的交換,由此影響大氣狀況(Bonan, 2008)。與裸地相比,植被覆蓋的地表面反射率較小,因而凈太陽(yáng)輻射較大,更多的能量可用于顯熱和潛熱。植被對(duì)全球輻射能量的影響與溫室氣體對(duì)氣候變化的影響同等重要(Betts,2000)。植被覆蓋的地表面粗糙度較大,湍流更強(qiáng),空氣和水汽混合動(dòng)力阻力較小(Liu et al., 2007),因而蒸騰較大,向大氣輸送的水汽可導(dǎo)致區(qū)域性大氣濕度的增加,改變當(dāng)?shù)亟邓退倪^程(Ellison et al., 2012)。植被覆蓋度高的地表面葉面積較大,根系更深。發(fā)達(dá)的根系能增加植被適應(yīng)干旱,穩(wěn)定水分脅迫條件下水分儲(chǔ)存的能力。這種功能是生態(tài)系統(tǒng)適應(yīng)氣候變化的重要機(jī)制(Jones et al., 2012)。大量觀測(cè)和模擬研究揭示了植被通過地—?dú)饨粨Q過程對(duì)氣候和水文過程的重要影響(符淙斌和袁慧玲,2001;Zheng et al., 2002;Liu et al., 2006;高學(xué)杰等,2007;Liu et al., 2008; 曾紅玲等,2010;Chen et al., 2012; Ma et al., 2013)。

    干旱—植被相互作用研究的一個(gè)新的重要議題是溫室效應(yīng)導(dǎo)致的氣候變化對(duì)干旱和植被的影響。氣候模型預(yù)測(cè)結(jié)果表明,溫室效應(yīng)將導(dǎo)致本世紀(jì)氣候發(fā)生顯著的變化(IPCC,2013)。例如,北美洲西部可能更為干燥(Cook et al.,2014),干旱發(fā)展可能更快更強(qiáng)(Trenberth et al.,2014)。中國(guó)21世紀(jì)與降水有關(guān)的事件更趨于極端化,極端降水強(qiáng)度可能增強(qiáng)(Jiang et al., 2011),干旱可能加重(Wang and Chen, 2014)。未來極端氣候的發(fā)生將對(duì)生態(tài)系統(tǒng)可能產(chǎn)生巨大的影響(吳紹洪等,2007)。開放灌叢和荒漠草原是受影響最為嚴(yán)重的類型,落葉闊葉林、有林草地和常綠針葉林也會(huì)受到影響。干旱和熱浪已導(dǎo)致世界各地森林枯死,預(yù)計(jì)這種情況會(huì)隨著氣候變化而增加(Anderegg et al.,2014)。

    氣候變化也會(huì)改變干旱—植被相互作用的另一個(gè)過程,也就是植被對(duì)氣候的影響 (Falloon et al.,2012)。盡管預(yù)測(cè)未來全球干旱趨勢(shì)的不確定性很大,一些研究表明一些特定區(qū)域的預(yù)測(cè)有較大的一致性。未來植被稀少的干燥地區(qū)的干旱會(huì)比其他地區(qū)增加的可能性更大。Seager et al.(2007)預(yù)測(cè),北美干燥的西南部會(huì)強(qiáng)烈干旱化,而且這種變干趨勢(shì)正在發(fā)生,今后幾十年美國(guó)西南部將再現(xiàn)沙塵暴肆虐的情形。相對(duì)于干旱對(duì)生態(tài)系統(tǒng)的影響,我們對(duì)植被在未來干旱變化中的作用的認(rèn)識(shí)非常有限,這限制了預(yù)測(cè)氣候變化情形下干旱趨勢(shì)的能力。

    干旱是一種極端天氣氣候事件。與其他一些極端事件一樣,當(dāng)前和氣候變化情形下未來干旱的預(yù)測(cè)都具有極大的挑戰(zhàn)性。多種因素(主要是海洋表面和大氣的異常)會(huì)導(dǎo)致干旱的形成和發(fā)展,而包括植被在內(nèi)的陸面狀況也是一個(gè)影響干旱強(qiáng)度和持續(xù)時(shí)間,以及空間和季節(jié)性變化的因素。因此,了解植物的作用將有助于減少預(yù)測(cè)未來干旱趨勢(shì)的不確定性。此外,減緩氣候變化的可能影響是科學(xué)研究、管理和決策部門共同關(guān)心的問題和職責(zé)。森林具有通過光合作用減少大氣中的碳的能力,從而減緩氣候變化影響。此外,全球增暖會(huì)增加植被蒸騰消耗的熱量,這會(huì)一定程度上抵消全球增暖(Liu, 2010)。因此,研究氣候變化情形下植被對(duì)干旱的影響,對(duì)于發(fā)展減緩氣候變化影響的策略有著重要的意義。

    觀測(cè)、動(dòng)力和數(shù)值分析揭示了植被對(duì)干旱的復(fù)雜影響。多數(shù)研究表明植被—干旱之間的負(fù)反饋機(jī)制,即植被增加(減少),干旱減弱(加強(qiáng))。丁一匯等(2005)發(fā)現(xiàn),內(nèi)蒙古地區(qū)植被覆蓋減少會(huì)導(dǎo)致中國(guó)北方大部分地區(qū)降水減少,加劇了華北、西北地區(qū)的干旱。長(zhǎng)時(shí)間干旱將導(dǎo)致植被狀況變化,由此形成的植被—降水反饋是干旱增強(qiáng)和持續(xù)的原因之一。Charney (1975), Charney et al. (1977) 認(rèn)為,過度放牧和砍伐森林增加地表反照率,減少地表吸收的太陽(yáng)幅射和感熱,這將減少對(duì)流降水,植被覆蓋也因而進(jìn)一步減少。這一反饋機(jī)制被用來解釋上世紀(jì)70年代初到80年代中期北非薩赫勒地區(qū)的持續(xù)干旱。植被—降水相互作用的研究多偏重氣象干旱。由于與植物根系的直接聯(lián)系,作為水文和農(nóng)業(yè)干旱特征之一的土壤水分不足常常是干旱影響生態(tài)系統(tǒng)研究的一個(gè)核心。Zeng et al.(1999)認(rèn)為,植被覆蓋減少導(dǎo)致較小蒸騰和土壤水分向大氣的傳輸。這將減少降水,植被覆蓋也相應(yīng)變得更小。由此提出了一種植被影響土壤濕度、蒸騰蒸散、和降水的反饋機(jī)制。

    然而,最近的一些研究表明,植被也可能加強(qiáng)干旱。Teuling et al.(2013)分析了兩次歐洲干旱事件的觀測(cè)結(jié)果,發(fā)現(xiàn)干旱期間減少的云量導(dǎo)致陸地表面吸收更多的太陽(yáng)輻射,因此溫度較高。這會(huì)增加氣候濕潤(rùn)、植被覆蓋較大地區(qū)的蒸騰,增強(qiáng)土壤水分不足的程度。土壤貯水量也因此變得更小,生態(tài)干旱變得更加嚴(yán)重。Meng et al.(2014)模擬了澳大利亞東南部2002年的干旱,發(fā)現(xiàn)土壤水分、植被與降水變化之間存在緩慢的反饋過程。土壤水分由于干旱而變小,由此引起的蒸發(fā)減少不利于對(duì)流和降水,短期內(nèi)導(dǎo)致土壤濕度進(jìn)一步減少。但如果土壤水分持續(xù)不足,植被將減少,蒸發(fā)量會(huì)減少,從而減少土壤水分下降的速度,減緩干旱。

    上述研究表明,植被可能減少干旱,也可能加強(qiáng)干旱。那么,這種不確定性對(duì)未來氣候變化情形下干旱趨勢(shì)有什么影響?本文通過分析美國(guó)干旱強(qiáng)度和頻率的變化趨勢(shì)與植被的關(guān)系對(duì)這一問題進(jìn)行探索。具體研究?jī)?nèi)容包括預(yù)測(cè)未來干旱趨勢(shì)及其地理分布和季節(jié)變化,比較高、低植被覆蓋度區(qū)域未來干旱強(qiáng)度和頻率,分析植被影響的機(jī)制和出現(xiàn)的氣候條件,及探討對(duì)適應(yīng)和減緩氣候變化的意義。

    2 方法

    2.1 干旱指數(shù)

    本文采用Keetch-Byram干旱指數(shù)(KBDI)(Keetch and Byram, 1968)表示干旱強(qiáng)度。KBDI是土壤回到飽和狀態(tài)所需的水分,由簡(jiǎn)化土壤水分平衡方程獲得:

    Qt和Qt?1是當(dāng)前和前一天的KBDI,單位為10?2英寸(1英寸等于2.54 cm),ΔQ (=Qt-Qt?1) 是土壤水分短缺(KBDI)在 Δt時(shí)段(這里為日)的變化。Qmax為土壤飽和時(shí)所需水分(假定為8英寸或0.2 m)。KBDI分別隨蒸發(fā)QE和降水Pt增加和減少。蒸發(fā)由三個(gè)因子決定。一是當(dāng)日最高溫度Tt(K);溫度越高,潛在蒸發(fā)越大。二是多年平均降水Pave,近似表示植被狀況(假設(shè)高降水區(qū)域能支持更多的植被);植被越多,蒸騰蒸發(fā)也越大。三是前一天缺水程度;程度越小,蒸發(fā)越大。ci(i=0, …, 5) 為經(jīng)驗(yàn)常數(shù)。此外,還假設(shè)溫度低于283 K (10°C)時(shí)KBDI不再變化。如果前一天無(wú)降水,僅當(dāng)日降水量大于0.20英寸(0.5 cm)的部分用于計(jì)算。

    KBDI對(duì)氣候變化的響應(yīng)(溫度和降水變化,假設(shè)植被不變)為:

    干旱增幅(ΔQt)由下列因子決定:溫度Tt和其變化ΔTt(正比)、降水變化ΔtP(反比)、當(dāng)前平均降水(植被覆蓋度)Pave(正比)、當(dāng)前土壤水分短缺或干旱程度Qt?1(反比)。Pave越大,公式(4)分母中e指數(shù)項(xiàng)越小,Δt前分?jǐn)?shù)項(xiàng)越大,因此ΔQt與Pave成正比。這意味著,在其他氣象條件相同的情況下,植被覆蓋度越大,蒸騰蒸發(fā)能力越強(qiáng),土壤水分損失越快,干旱發(fā)展越快。

    KBDI最大值為800。干旱程度分為四級(jí):低(KBDI = 0~200)、中度(200~400)、高(400~600)和極端(600~800)。 KBDI的經(jīng)驗(yàn)系數(shù)根據(jù)美國(guó)東南部觀測(cè)資料獲得。一些研究發(fā)現(xiàn)通常情況下KBDI可用于美國(guó)及世界其他地區(qū) (Xanthopoulos et al.,2006)。Liu et al. (2010a, 2013) 應(yīng)用KBDI研究了氣候變化對(duì)全球和美國(guó)森林干旱和火災(zāi)的影響。

    目前有多種干旱指數(shù)可供選擇,每一種指數(shù)側(cè)重一種或多種干旱特征和影響(Heim,2002;Zargar et al.,2011)。其中帕爾默干旱指數(shù)(PDSI)和標(biāo)準(zhǔn)化降水指數(shù)(SPI)應(yīng)用較廣。PDSI(Palmer,1965)已被用來研究各種尺度上的干旱(Dai, 2011;Zhang et al.,2013)。PDSI和KBDI都是基于水分平衡方程,溫度和降水為主要因子。兩者都依賴于前期的狀況,因而能較好地反映長(zhǎng)時(shí)間缺水的累積效應(yīng)。PDSI更為復(fù)雜全面,主要是使用月值,而KBDI多使用日值。KBDI和PDSI都始于上世紀(jì)六十年代。由于缺乏獲得土壤濕度狀況的手段,它們應(yīng)用經(jīng)驗(yàn)關(guān)系(如公式2)根據(jù)氣象狀況估算蒸發(fā)。目前氣候模式中包括模擬地—?dú)庀嗷プ饔玫年懙剡^程模式,能直接模擬蒸發(fā),為評(píng)估干旱提供了新的途徑。盡管如此,許多研究(包括本文)仍采用干旱指數(shù)中的公式計(jì)算蒸發(fā)。原因之一是有些陸地過程模式在模擬極端天氣情形下土壤濕度的變化有一定的不確定性。例如,Liu et al.(2010b)比較了應(yīng)用一個(gè)區(qū)域氣候模式模擬的氣象因子和土壤濕度分別計(jì)算1988年北美大旱期間的KBDI,發(fā)現(xiàn)前者更能反映干旱的強(qiáng)度。

    2.2 區(qū)域氣候變化

    研究區(qū)域是美國(guó)大陸。氣候變化資料由北美區(qū)域氣候變化評(píng)估項(xiàng)目(NARCCAP)提供(Mearns et al., 2012)。NARCCAP產(chǎn)生高分辨率氣候變化預(yù)測(cè)結(jié)果,用以評(píng)估未來區(qū)域氣候預(yù)測(cè)的不確定性,并為研究氣候變化的影響提供氣候資料。資料為1971~2000年和2041~2070兩個(gè)時(shí)期,多數(shù)氣象要素輸出頻率為三小時(shí),空間分辯率為50公里。未來氣候預(yù)測(cè)降尺度只包含SRES A2排放情景。

    區(qū)域氣候變化情景由動(dòng)力降尺度獲得。大氣環(huán)流模型(GCM)包括全球耦合氣候模式(CGCM3)、地球物理流體動(dòng)力學(xué)實(shí)驗(yàn)室氣候模型(GFDL)、哈德利中心氣候模式(HadCM3)和氣候系統(tǒng)模式(CCSM)。區(qū)域氣候模型(RCM)包括加拿大區(qū)域氣候模式(CRCM)、區(qū)域譜模式(RSM)、高分辨率區(qū)域模型(HRM)、中尺度氣象模式(MM5I)區(qū)域氣候模式3(RCM3)和天氣研究和預(yù)報(bào)模式(WRFG)。NARCCAP總共提供10 個(gè)GCM與RCM組合產(chǎn)生的氣候變化情景,即CGCM3與CRCM、RCM3和WRFG組合,GFDL與HRM3 和RCM3組合,HadCM3與HRM3和MM5I組合,CCSM與CRCM、MM5I和WRFG組合。本研究使用前七種組合。由CCSM產(chǎn)生的三種組合缺每個(gè)時(shí)期最后一年的資料,因此未用于計(jì)算。值得一提的是,NARCCAP所用的GCM和RCM均未模擬氣候變化情形下植被類型、結(jié)構(gòu)、分布等特征的變化。資料可從http://www.narccap.ucar.edu[2015-06-21]獲得。

    2.3 分析植被對(duì)干旱的影響

    應(yīng)用NARCCAP氣候變化情形分別計(jì)算了美國(guó)大陸現(xiàn)在和未來兩個(gè)時(shí)期逐日KBDI。每個(gè)格點(diǎn)KBDI計(jì)算步驟如下:(i)前五年每一年初值為零,計(jì)算逐日KBDI。(ii)獲得每年最后一天KBDI的平均值。(iii)以此值作為初值,計(jì)算30年逐日KBDI。除了KBDI,還計(jì)算了強(qiáng)干旱(KBDI>400)和極端干旱(KBDI>600)日出現(xiàn)的百分率。通過比較每個(gè)氣候區(qū)兩種植被類型(高、低植被覆蓋度)之間KBDI,高干旱及極端干旱風(fēng)險(xiǎn)的出現(xiàn)頻率的差異,分析植被對(duì)氣候變化情形下干旱趨勢(shì)的影響。

    氣候區(qū)根據(jù)Koppen-Geiger氣候分類(Kottek et al., 2006)劃分。每一氣候區(qū)由主要?dú)夂蝾愋汀⒔邓穸龋┖蜏囟忍卣鳑Q定。美國(guó)大陸主要?dú)夂蝾愋桶ǜ稍?、溫暖、寒冷;降水(濕度)類型包括濕?rùn)、干夏、荒漠、沙漠;溫度類型包括熱夏、暖夏、冷夏、熱干、冷干。美國(guó)大陸主要分為六個(gè)氣候區(qū)(圖1a):東南(SE)溫暖、濕潤(rùn)、夏季炎熱區(qū),東北(NE)寒冷、潮濕、夏季溫暖或炎熱區(qū),西部?jī)?nèi)地(W)干旱和高寒草原區(qū),西南(SW)干旱草原和沙漠以及夏季干熱或干冷區(qū),落基山脈(RM)區(qū)(與NE氣候類似),太平洋沿岸(PC)溫暖區(qū)(夏季干燥冬季濕潤(rùn))。

    美國(guó)東南地區(qū)主要土地覆蓋類型為落葉林、針葉林、混交林與濕地,農(nóng)業(yè)分布在密西西比河流域、德克薩斯東部、以及大西洋沿岸的一些地區(qū)。東北地區(qū)主要為落葉林,農(nóng)業(yè)區(qū)主要在中西部地區(qū)。太平洋沿海地區(qū)以針葉林和草為主,而西部山區(qū)和荒漠地區(qū)主要為針葉林、干草和沙漠。

    不同區(qū)域氣候模式中的土地類型方案可能有所不同。這里采用RCM3中擴(kuò)充的BATS土地類型方案(在原BATS方案中增加了與城市和城鎮(zhèn)有關(guān)的兩種類型)(Dickinson et al.,1993;Elguindi et al., 2004)。每個(gè)氣候區(qū)選擇兩種土地覆蓋類型代表高、低植被覆蓋度區(qū)域進(jìn)行比較(表1)。東南和東北氣候區(qū)為落葉林和農(nóng)田。這兩個(gè)氣候區(qū)也包含較多的林地/城鎮(zhèn),由于其許多特征與落葉林相近,特別是難以區(qū)分林地與城鎮(zhèn)面積的比例,所以未選其代表低植被覆蓋度區(qū)域。西部和西南氣候區(qū)為草地和荒漠草原,落基山和太平洋沿岸氣候區(qū)為針葉林和荒漠草原。為敘述方便,以下通稱落葉林和針葉林為林地,簡(jiǎn)稱荒漠草原為荒漠。

    表1 主要植被類型的格點(diǎn)數(shù)。括號(hào)中m和l為進(jìn)行比較所選擇高、低植被覆蓋度類型Table 1 Grid point numbers of major land cover types. The letters ‘m’ and ‘l’ in parentheses indicate selected types as more and less vegetated lands for comparison

    3 結(jié)果

    3.1 氣候特征和變化

    冬季最高氣溫(圖1)從美國(guó)東南地區(qū)約290 K降至東北地區(qū)西北部的260 K。落基山地區(qū)為265~270 K,西部地區(qū)為270~275 K,均較同緯度的其他地區(qū)寒冷,而太平洋沿岸地區(qū)為275~285 K,較同緯度其他地區(qū)溫暖。夏天地理分布與冬季類似,但氣溫變得炎熱,東南和西南地區(qū)達(dá)到300~310 K。春季和秋季介于冬夏之間,為280~300 K(落基山區(qū)氣溫較低)。

    未來全美國(guó)所有季節(jié)的最高氣溫度都有所增加(圖2)。冬季東北地區(qū)、春季西南地區(qū)增溫達(dá)到2.5 K。夏季全國(guó)大部分地區(qū)增溫3~5 K,東南地區(qū)中部和西部地區(qū)最大。秋季增溫地理分布與夏季類似,但幅度較小,僅在西部地區(qū)東南部達(dá)到2~3 K。

    冬季降水(圖3)主要分布在東南地區(qū)和太平洋沿岸地區(qū)北部,降水量約4~6 mm d?1。落基山區(qū)也有較大降水。春季降水區(qū)東南地區(qū)變大,而太平洋沿岸地區(qū)變小。夏季降水主要集中在東北地區(qū),大西洋沿岸降水量達(dá)4~4.5 mm d?1,但東南地區(qū)降水減少,太平洋沿岸地區(qū)降水區(qū)幾乎消失。秋天降水分布和冬季相似,但降雨量較小。

    未來美國(guó)冬天降水量(圖4)普遍增加,東南和東北地區(qū)大西洋沿岸增幅達(dá)到0.4 mm d?1。春季西南地區(qū)、西部地區(qū)南部、和太平洋沿岸降水開始減少,而夏天除東北地區(qū)西北部和西部地區(qū)東北部,全國(guó)范圍內(nèi)降水減少,東南地區(qū)西部和西部地區(qū)東南部減幅達(dá)到0.6 mm d?1。秋天全國(guó)范圍降水又開始增加,但除西部地區(qū)的一些地方,增幅較小。

    同一氣候區(qū)不同植被區(qū)域的溫度和降水存在很大差別。在美國(guó)東南、東北、落基山和太平洋沿岸氣候區(qū),與農(nóng)田(荒漠)相比,林地(草地)通常溫度較低,降水較大(圖5、6)。唯一例外是東北氣候區(qū)冬季,盡管林地(草地)降水較大,溫度卻較高。與上述氣候區(qū)相反,美國(guó)西部和西南氣候區(qū)林地(草地)溫度與對(duì)應(yīng)的農(nóng)田(荒漠)溫度相近或更大。降水量差別取決于季節(jié)。春夏季林地(草地)降水較大。

    圖1 當(dāng)前最高氣溫(單位:K)(7個(gè)區(qū)域氣候變化情景平均)。(a)、(b)、(c)和(d)分別為冬、春、夏和秋季。(a)中SE、NE、W、SW、RM 和PC分別為東南、東北、西部、西南、落基山、太平洋沿岸氣候區(qū),白線為近似區(qū)界Fig. 1 Present maximum temperature (units: K) averaged over seven regional climate change scenarios in (a) winter, (b) spring, (c) summer, and (d) fall. SE, NE, W, SW, RM and PC in (a) represent the southeast, northeast, west, southwest, Rocky Mountains, and Pacific coast climate regimes, respectively. White lines are region boundaries

    圖2 未來最高氣溫變化(單位:K)(7個(gè)區(qū)域氣候變化情景平均)。其他同圖1Fig. 2 As in Fig. 1 but for future changes in maximum temperature (units: K)

    圖3 當(dāng)前降水量(mm d?1)(7個(gè)區(qū)域氣候變化情景平均)。其他同圖1Fig. 3 As in Fig. 1 but for present precipitation (units: mm d?1)

    圖4 未來降水量變化(單位:mm d?1)(7個(gè)區(qū)域氣候變化情景平均)。其他同圖1Fig. 4 As in Fig. 1 but for future changes in seasonal precipitation (units: mm d?1)

    圖5 不同植被類型區(qū)域最高氣溫(單位:K)。(a)、(b)、(c)和(d)分別為冬、春、夏和秋季。T和ΔT分別為當(dāng)前最高氣溫和未來變化。SE、NE、W、SW、RM和PC分別為東南、東北、西部、西南、落基山和太平洋沿岸氣候區(qū)。每一氣候區(qū)左側(cè)表示林地(草地),右側(cè)表示對(duì)應(yīng)的農(nóng)田(荒漠)。七個(gè)區(qū)域氣候變化情景平均,其最小和最大值區(qū)間由垂直線表示Fig. 5 Maximum temperature (units: K) with different land cover types in (a) winter, (b) spring, (c) summer and (d) fall. T and ΔT represent present conditions and future changes, respectively. SE, NE, W, SW, RM and PC represent the southeast, northeast, west, southwest, Rocky Mountains, and Pacific coast climate regimes, respectively. Left and right bars for each region represent forest land (rangeland) and the corresponding farmland (dryland), respectively. Averaged over seven climate change scenarios, with ranges indicated by vertical lines

    3.2 干旱特征和趨勢(shì)

    冬季太平洋沿岸和西部地區(qū)西部KBDI超過500(圖7),德克薩斯為350~400。KBDI從西向東減小,東北地區(qū)和大部分東南地區(qū)低于100。春季的地理分布不變,但數(shù)值略有減小。夏季和秋季東南和東北地區(qū)KBDI顯著增加,其中東南地區(qū)夏季KBDI達(dá)到400,秋天達(dá)到500。

    未來冬季KBDI普遍增加,增幅達(dá)到100以上(圖8),但德克薩斯和加利福尼亞州南部略有減少。春天KBDI變化相似,但數(shù)值較小。此后,全國(guó)范圍KBDI急劇增加,夏季東南地區(qū),秋季大西洋中部沿岸和西部地區(qū)東部增幅達(dá)150。夏秋季KBDI的顯著增加和這些地區(qū)強(qiáng)烈增暖和變干一致。 強(qiáng)干旱和極端干旱日出現(xiàn)頻率(圖未顯示)的地理分布、季節(jié)變化、及未來趨勢(shì)與KBDI的特征類似。

    3.3 植被的影響

    與農(nóng)田(荒漠)區(qū)域相比,對(duì)應(yīng)的林地(草地)所有氣候區(qū)所有季節(jié)的KBDI都較?。▓D9)。這一差別對(duì)每個(gè)氣候變化情景都一樣。在氣候與植被均為動(dòng)態(tài)系統(tǒng)的模式中,這一關(guān)系可能表明植被會(huì)降低干旱風(fēng)險(xiǎn),也可能表明干旱風(fēng)險(xiǎn)低的區(qū)域植被覆蓋度高。由于獲得NARCCAP氣候變化的區(qū)域氣候模式中這兩個(gè)系統(tǒng)只有氣候是動(dòng)態(tài)的,上述關(guān)系應(yīng)代表第一種可能性,即表明植被會(huì)降低干旱風(fēng)險(xiǎn)。

    圖6 不同植被類型區(qū)域降水(mm d?1)。P和Δ P分別為當(dāng)前降水量和未來變化。其他同圖5Fig. 6 As in Fig. 5 but for precipitation (units: mm d?1). P and Δ P represent present conditions and future changes, respectively

    但是植被對(duì)未來干旱增強(qiáng)的影響卻因氣候區(qū)而異。六個(gè)氣候區(qū)可分為三種情形:(i)東南、西部和西南氣候區(qū):與農(nóng)田(荒漠)地相比,所有季節(jié)對(duì)應(yīng)林地(草地)KBDI增幅較大。唯一例外是東南區(qū)春天。(ii)東北和落基山氣候區(qū):所有季節(jié)林地(草地)KBDI增幅較小。(iii)太平洋沿岸氣候區(qū):該區(qū)屬于地中海氣候,夏秋溫暖但少雨,冬春低溫但多雨,這和美國(guó)其他區(qū)域溫度和降水的對(duì)應(yīng)關(guān)系(即雨季和暖季、干季和冷季同時(shí)出現(xiàn))相反。該區(qū)植被對(duì)未來干旱增強(qiáng)的影響兼有上述兩種情形,主要由溫度特征決定:夏秋季和東南、西部、和西南氣候區(qū)一樣(情形i),林地(草地)KBDI增幅較大;冬春季和東北、落基山氣候區(qū)一樣(情形ii),林地(草地)KBDI增幅較小。

    植被對(duì)未來高干旱和極端干旱風(fēng)險(xiǎn)日出現(xiàn)頻率增加會(huì)產(chǎn)生同樣的影響。 圖10為高干旱的情形。因此,在溫暖和濕潤(rùn)/干燥氣候區(qū)(Koppen-Geiger氣候分類中主要?dú)夂蝾愋虲和B),林地(草地)干旱強(qiáng)度和頻率增幅較大,而在寒冷和濕潤(rùn)氣候區(qū)(主要?dú)夂蝾愋虳),林地(草地)干旱強(qiáng)度和頻率增幅較小。

    3.4 物理機(jī)制

    3.4.1 當(dāng)前干旱狀況

    如公式(1)和(2)所示,KBDI與溫度和降水直接相關(guān);植被會(huì)影響溫度以及降水,從而間接影響KBDI。過去很多研究結(jié)果表明,植被覆蓋的土壤通常溫度較低,降水較大,這和現(xiàn)在美國(guó)東南、東北、落基山和太平洋沿岸氣候區(qū)植被狀況與溫度和降水的關(guān)系一致(圖5、6)。因此,這些氣候區(qū)林地(草地)的KBDI較小。

    圖7 當(dāng)前KBDI(Keetch-Byram干旱指數(shù))(7個(gè)區(qū)域氣候變化情景平均)。其他同圖1Fig. 7 As in Fig. 1 but for present KBDI (Keetch-Byram Drought Index)

    圖8 未來KBDI變化(7個(gè)區(qū)域氣候變化情景平均)。其他同圖1Fig. 8 As in Fig. 1 but for future changes in seasonal KBDI

    圖9 不同植被類型區(qū)域KBDI。KBDI和Δ KBDI分別為當(dāng)前干旱指數(shù)和未來變化。其他同圖5Fig. 9 As in Fig. 5 but for KBDI with different land cover types. KBDI and Δ KBDI represent present conditions and future changes, respectively

    美國(guó)西部和西南氣候區(qū)林地(草地)的溫度與對(duì)應(yīng)的農(nóng)田(荒漠)區(qū)域相近或更大。春夏季降水林地(草地)較大??紤]到這兩個(gè)氣候區(qū)水分為地—?dú)饨粨Q更為關(guān)健的因素,可以假定較大降水比較高氣溫起著更大的作用,因此KBDI較低。冬秋季降水林地(草地)較小或差別不大。較高氣溫應(yīng)導(dǎo)致較大KBDI。這和計(jì)算結(jié)果不一致。一種可能的解釋是,KBDI取決于其前期值;由于夏季KBDI較小,因而秋冬天的KBDI相應(yīng)地較小。

    3.4.2 未來干旱強(qiáng)度變化

    如公式(3)和(4)所示,氣溫和降水量的未來變化也是未來KBDI變化的主要因素。與農(nóng)田(荒漠)相比,對(duì)應(yīng)的林地(草地)區(qū)域溫度、降水和KBDI的變化可能有以下情形:(i)增暖和變干更明顯(后者可以是降水減少更多或降水增加更?。?;預(yù)期KBDI增幅更明顯(三種變化特征分別表示為+ΔT,-ΔP和+ΔQ)。(ii)增暖和變干較弱(后者可以是降水減少較小或降水增加更大);預(yù)期KBDI增幅較?。ǎ,+ΔP和-ΔQ)。(iii)增暖更明顯但變干較弱(+ΔT和+ΔP),或增暖較弱但變干更明顯(ΔT和-ΔP); 難以預(yù)期KBDI增幅較大還是較?。ā捆)。

    東南、西部和西南氣候區(qū)(所有季節(jié))和太平洋沿岸氣候區(qū)(夏秋季)溫度、降水,或KBDI變化在不同區(qū)域和季節(jié)之間很不一致。例如,ΔT符號(hào)為正、負(fù)和幾乎為零的次數(shù)為4、 4和6(表2)。其結(jié)果是,預(yù)期和計(jì)算的KBDI變化同號(hào)的次數(shù)與總次數(shù)的比例僅為3/14,故情形(iii)出現(xiàn)最多。因此無(wú)法通過未來氣溫和降水的變化來判斷植被對(duì)未來干旱增強(qiáng)的影響。這意味著不得不尋找其他機(jī)制。一個(gè)可能的機(jī)制為,如公式(3)和(4)所示,ΔT對(duì)KBDI增加的重要性與當(dāng)前干旱水平(KBDI)成反比。由于KBDI明顯較小,這些氣候區(qū)林地(草地)區(qū)域由氣候變暖導(dǎo)致的干旱增幅度更明顯。

    相反,東北、落基山(所有季節(jié))和太平洋沿岸氣候區(qū)(冬春季)溫度、降水,或KBDI變化不同區(qū)域和季節(jié)之間較為一致。例如,ΔT符號(hào)為正、負(fù)和幾乎為零的次數(shù)為1、8和1。其結(jié)果是,預(yù)期的和計(jì)算的KBDI變化同號(hào)的比例較高,達(dá)到6/10,故情形(ii)出現(xiàn)最多。四次不同符號(hào)情形有三次出現(xiàn)在東北氣候區(qū)。該區(qū)KBDI很小,因此增溫作用很重要。ΔT都為負(fù),由此導(dǎo)致林地(草地)干旱增加較弱。因此,這些氣候區(qū)植被通過影響未來氣溫和降水變化(特別是氣溫)間接影響干旱變化是一個(gè)主要物理機(jī)制。

    表2 比較不同植被類型之間未來溫度(ΔT)、降水(ΔP)和KBDI(Keetch-Byram干旱指數(shù))(ΔQ)變化?!埃薄ⅰ埃焙汀啊碧?hào)分別表示林地(草地)比對(duì)應(yīng)的農(nóng)田(荒漠)區(qū)域變幅更大、更小和相近。預(yù)計(jì)的ΔQ符號(hào)根據(jù)ΔT和ΔP的符號(hào)推斷,x表示推斷與計(jì)算的ΔQ同號(hào)Table 2 Comparisons of changes in temperature (ΔT), precipitation (ΔP), and KBDI (Keetch-Byram Drought Index) (ΔQ) between different land cover types. The ‘+’,‘-’ and ‘~’ symbols indicate the magnitude of change on more vegetated lands is larger, smaller, and almost same, respectively. The expected sign of ΔQ is obtained based on the signs of ΔT and ΔP, with an ‘x’ indicating the same expected and calculated signs

    4 討論

    (1)未來干旱趨勢(shì)

    一些研究對(duì)氣候變化情形下未來干旱趨勢(shì)進(jìn)行了預(yù)測(cè)(Dai,2011)。 Cook et al.(2015)根據(jù)CMIP5氣候預(yù)測(cè)結(jié)果分析了PDSI和土壤水分,發(fā)現(xiàn)美國(guó)中部大平原和西南地區(qū)夏季干旱將急劇增加。本文分析的KBDI也發(fā)現(xiàn)類似結(jié)果。所不同的是KBDI夏季最大增加出現(xiàn)在美國(guó)東南部,而中部大平原干旱增加也出現(xiàn)在其他季節(jié)。Liu et al.(2013)使用一個(gè)NARCCAP氣候變化情景(HadCM3-HRM3)計(jì)算了未來KBDI變化。七個(gè)NARCCAP氣候變化情景的集成分析結(jié)果顯示了類似的干旱增加趨勢(shì)及其地理分布和季節(jié)變化,但冬季和春季干旱增加主要出現(xiàn)在大平原的中部而不是整個(gè)區(qū)域;此外,西南地區(qū)干旱增加的幅度較小。

    (2)植被影響干旱的機(jī)理

    如引言所述,已有研究表明了植被對(duì)干旱影響的復(fù)雜性,植被可能減少干旱(Charney,1975;Zeng et al.,1999),也可能加強(qiáng)干旱(Teuling et al.,2013;Meng et al.,2014)。本研究為未來氣候變化情形下這一特征的存在提供了證據(jù),并揭示其與氣候區(qū)的關(guān)系和機(jī)制。在溫暖和潮濕/干燥的氣候區(qū),當(dāng)前植被覆蓋度較高的地區(qū)干旱風(fēng)險(xiǎn)較小,但未來氣候變化情形下干旱強(qiáng)度和頻率增幅更大。然而,寒冷和潮濕氣候區(qū),植被導(dǎo)致相對(duì)較小的增暖,干旱強(qiáng)度和頻率增加的幅度相應(yīng)地也較小。這些結(jié)果有助于增加對(duì)植被與干旱相互作用的理解,改進(jìn)當(dāng)前和未來干旱模擬和預(yù)測(cè)的能力。這些特征根據(jù)在美國(guó)的區(qū)域氣候變化模擬結(jié)果獲得,是否在世界其它地區(qū)也有類似的特征有待進(jìn)一步的研究。

    值得一提的是,氣候變化情形下植被本身會(huì)改變。例如,Grimm et al.(2013)表明,由于氣候變化,2100年物種種類將發(fā)生廣泛的變化,由此可能會(huì)改變5%~20%美國(guó)陸地的生物群落成分。這種影響可能會(huì)改變公式(2)中降水和植被之間的關(guān)系,也會(huì)改變植被對(duì)干旱趨勢(shì)的影響。本文所應(yīng)用的區(qū)域氣候變化情形未包括氣候變化情形下植被的變化。因此無(wú)法研究未來植被變化對(duì)干旱趨勢(shì)的影響變化。這是今后研究需關(guān)注的一個(gè)問題。

    另一個(gè)需關(guān)注的問題是如何證明本文分析的不同干旱趨勢(shì)是植被的影響,而不是不同氣候特性干旱自身的演變特性。本文研究植被影響采用的方法是比較同一氣候區(qū)兩種不同植被類型干旱趨勢(shì)的差別。但由于氣候區(qū)范圍較大,地形也可能不一樣,兩種植被區(qū)域大氣的狀態(tài)和過程本身會(huì)有所差別,并對(duì)干旱趨勢(shì)產(chǎn)生一定影響。解決這一問題的途徑之一是應(yīng)用區(qū)域氣候模式進(jìn)行數(shù)值實(shí)驗(yàn),將兩種植被區(qū)域變成同一種植被(例如,在東南氣候區(qū)將農(nóng)田變成林地,反之亦然)。如果兩個(gè)區(qū)域的干旱趨勢(shì)差別不復(fù)存在,則可證明不同干旱趨勢(shì)主要是植被的影響,否則是其他因素的影響。此實(shí)驗(yàn)超出本文分析的內(nèi)容,但是今后研究值得嘗試的課題。

    圖10 不同植被類型區(qū)域高干旱風(fēng)險(xiǎn)日出現(xiàn)頻率。HD和ΔHD分別為當(dāng)前高干旱風(fēng)險(xiǎn)日出現(xiàn)頻率和未來變化。其他同圖5Fig. 10 As in Fig. 5 but for the percentage of high drought (HD) risk days with different land cover types. HD and ΔHD represent present conditions and future changes, respectively

    (3)減緩氣候變化

    由于森林在減少大氣二氧化碳中重要作用,廣泛認(rèn)為森林恢復(fù)和人工造林是減緩氣候變化的一個(gè)十分有效的手段。本研究結(jié)果表明,森林可能對(duì)未來干旱增加有著相反的作用。一方面,在寒冷和潮濕氣候區(qū),增加森林將減小未來干旱增加的幅度。植被的這種作用進(jìn)一步加強(qiáng)了森林對(duì)減緩氣候變化的意義。但另一方面,在溫暖和潮濕/干燥氣候區(qū),更多的植被可能增加未來的干旱,加大氣候變化的影響。其中一個(gè)影響是未來林火趨勢(shì)。植被越多,干旱越嚴(yán)重,發(fā)生森林火災(zāi)的可能行也越大。這也許是森林管理部門制定適應(yīng)和減緩氣候變化計(jì)劃和策略時(shí)必須考慮的一個(gè)問題。

    (4)KBDI

    KBDI代表的干旱與氣象條件之間的關(guān)系較為簡(jiǎn)單,而后者受植被狀況的影響。因此,比較容易通過考察溫度和降水的變化認(rèn)識(shí)植被對(duì)干旱的影響。但是,KBDI有一定的局限性。首先,KBDI量度絕對(duì)水分不足;不同氣候區(qū)域同一土壤缺水量可能意味著不同的干旱程度。例如,KBDI=400在美國(guó)東北地區(qū)較為少見,但在西南地區(qū)很常見,因此實(shí)際上干旱僅在東北出現(xiàn)。其次,KBDI只考慮氣溫和降水的貢獻(xiàn)。植被也可以改變風(fēng)及其他氣象條件 (Dai et al., 2004)。NARCCAP氣候變化情景預(yù)測(cè)未來美國(guó)中緯度地區(qū)風(fēng)速會(huì)減小(Liu et al.,2013)。這將影響未來地表水分和能量交換,從而影響到植被對(duì)干旱的影響。應(yīng)用KBDI無(wú)法考察這種影響。此外,KBDI使用的降水和植被之間的關(guān)系在有些地區(qū)不一定有效(Liu et al.,2013)。

    (5)氣候變化情景

    NARCCAP氣候變化情景由CMIP3全球氣候變化預(yù)測(cè)產(chǎn)生,并只考慮了A2排放情景。新的CMIP5全球氣候變化預(yù)測(cè)采用了新的排放情景,氣候模式也有所更新和改進(jìn)。CMIP5干旱預(yù)測(cè)結(jié)果與CMIP3有所不同。例如,CMIP5的分析表明全球干旱趨勢(shì)預(yù)測(cè)的不確定性有所增加,預(yù)測(cè)的中緯度地區(qū)干旱趨勢(shì)有所減小 (IPCC,2013)。應(yīng)用CMIP5降尺度氣候變化情景是加深對(duì)植被影響未來干旱趨勢(shì)認(rèn)識(shí)的一個(gè)有效途徑。一個(gè)所面臨的問題是,動(dòng)力氣候降尺度需大量的計(jì)算資源,因此將需要很長(zhǎng)一段時(shí)間制定并實(shí)施像NARCCAP一樣的大規(guī)模動(dòng)力降尺度研究項(xiàng)目。

    5 結(jié)論

    本研究考察了氣候變化情形下植被對(duì)干旱趨勢(shì)的影響和機(jī)制。結(jié)果表明,在美國(guó)大陸溫暖和潮濕/干燥氣候區(qū),高植被覆蓋度地區(qū)干旱增加更為明顯,但在寒冷和潮濕氣候區(qū),高植被覆蓋度區(qū)域干旱增加較小。由此可以得出結(jié)論,植被對(duì)美國(guó)干旱趨勢(shì)的影響與氣候區(qū)有關(guān)。在特定的氣候區(qū),植被可以放大全球變暖期間未來干旱的風(fēng)險(xiǎn)。

    致謝 作者感謝審稿人極有價(jià)值和建設(shè)性的問題和建議。此項(xiàng)研究是由美國(guó)自然科學(xué)基金會(huì)、農(nóng)業(yè)部和能源部聯(lián)合項(xiàng)目:地球系統(tǒng)模式預(yù)報(bào)年代際區(qū)域氣候 (NIFA-2013-35100-20516)。區(qū)域氣候變化情景數(shù)據(jù)由北美區(qū)域氣候變化評(píng)估計(jì)劃(NARCCAP)提供。

    參考文獻(xiàn)(References)

    Anderegg W R L, Anderegg L D L, Berry J A, et al. 2014. Loss of whole?tree hydraulic conductance during severe drought and multi?year forest die?off [J]. Oecologia, 175 (1): 11–23, doi:10.1007/s00442-013- 2875-5.

    Betts R A. 2000. Offset of the potential carbon sink from boreal forestation by decreases in surface albedo [J]. Nature, 408: 187–190, doi:10.1038/35041545. Bonan G B. 2008. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests [J]. Science, 320: 1444–1449, doi:10.1126/ science.1155121.

    Charney J, Quirk W J, Chow S -H, et al. 1977. A comparative study of the effects of albedo change on drought in semi-arid regions [J]. J. Atmos. Sci., 34: 1366–1385, doi:10.1175/1520-0469(1977)034<1366: ACSOTE> 2.0.CO; 2.

    Charney J G. 1975. Dynamics of deserts and drought in the Sahel [J]. Q. J. R. Meteor. Soc., 101 (428): 193–202, doi:10.1002/qj.49710142802.

    Chen G S, Notaro M, Liu Z Y, et al. 2012. Simulated local and remote biophysical effects of afforestation over the southeast United States in boreal summer [J]. J. Climate, 25: 4511–4522, doi:10.1175/JCLI-D-11-00317.1.

    Cook B I, Ault T R, Smerdon J E. 2015. Unprecedented 21st century drought risk in the American Southwest and Central Plains [J]. Sci. Adv., 1(1): e1400082, doi:10.1126/sciadv.1400082.

    Cook B I, Smerdon J E, Seager R, et al. 2014. Global warming and 21st century drying [J]. Clim. Dyn., 43 (9–10): 2607–2627. doi:10.1007/ s00382-014-2075-y.

    Dai A. 2011. Drought under global warming: A review [J]. Wiley Interdisciplinary Reviews: Climate Change, 2: 45–65, doi:10. 1002/wcc.81.

    Dai Y J, Dickinson R R, Wang Y P. 2004. A two-big-leaf model for canopy temperature, photosynthesis, and stomatal conductance [J]. J. Climate, 17 (12): 2281–2299, doi:10.1175/1520-0442(2004)017<2281:ATMFCT>2.0.CO;2.

    Dickinson R E, Henderson-Sellers A, Kennedy P J. 1993. Biosphereatmosphere Transfer Scheme (BATS) Version 1e as Coupled to the NCAR Community Climate Model [R]. NCAR Technical Note NCAR/ TN-387+STR, doi:10.5065/D67W6959.

    丁一匯 李巧萍 董文杰. 2005. 植被變化對(duì)中國(guó)區(qū)域氣候影響的數(shù)值模擬研究 [J]. 氣象學(xué)報(bào), 63 (5): 613–621. Ding Y H, Li Q P, Dong W J. 2005. A numerical simulation study of the impacts of vegetation changes on regional climate in China [J]. Acta Meteor. Sinica, 63 (5): 613–621, doi:10.11676/qxxb2005.060

    Elguindi N, Bi X, Giorgi F, et al. 2004. RegCM Version 3 User’s guide [R]. PWCG Abdus Salam ICTP, Trieste, Italy, 48pp.

    Ellison D, Futter M N, Bishop K. 2012. On the forest cover–water yield debate: From demand- to supply-side thinking [J]. Global Change Biol., 18: 806–820, doi:10.1111/j.1365-2486.2011.02589.x.

    Falloon P D, Dankers R, Betts R A, et al. 2012. Role of vegetation change in future climate under the A1B scenario and a climate stabilisation scenario, using the HadCM3C Earth system model [J]. Biogeosciences, 9: 4739–4756, doi:10.5194/bg-9-4739-2012.

    符淙斌, 袁慧玲. 2001. 恢復(fù)自然植被對(duì)東亞夏季氣候和環(huán)境影響的一個(gè)虛擬試驗(yàn) [J]. 科學(xué)通報(bào), 46 (8): 691-695. Fu C B, Yuan H L. 2001. An virtual numerical experiment to understand the impacts of recovering natural vegetation on the summer climate and environmental conditions in East Asia [J]. Chin. Sci. Bull., 46: 1199–1203, doi:10.1007/ BF02900602.

    高學(xué)杰, 張冬峰, 陳仲新,等. 2007. 中國(guó)當(dāng)代土地利用對(duì)區(qū)域氣候影響的數(shù)值模擬 [J]. 中國(guó)科學(xué)(D輯): 地球科學(xué), 37 (3): 397–404. Gao X J, Zhang D F, Chen Z X, et al. 2007. Land use effects on climate in China as simulated by a regional climate model [J]. Science in China (Series D): Earth Sciences, 50 (4): 620–628, doi:10.1007/s11430-007-2060-y.

    Grimm N B, Chapin Ⅲ F S, Bierwagen B, et al. 2013. The impacts of climate change on ecosystem structure and function [J]. Front. Ecol. Environ., 11 (9): 474–482, doi:10.1890/120282.

    Heim R R. 2002. A review of twentieth-century drought indices used in the United States [J]. Bull. Am. Meteor. Soc., 83: 1149–1166.

    IPCC. 2013. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [M]// Stocker T F, Qin D, Plattner G -K, et al. (eds.). Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press, 1535pp, doi: 10.1017/CBO9781107415324.

    Jiang Z H, Song J, Li L, et al. 2011. Extreme climate events in China: IPCC-AR4 model evaluation and projection [J]. Climatic Change, 110 (1–2): 385–401doi:10.1007/s10584-011-0090-0.

    Jones J A, Creed I F, Hatcher K L, et al. 2012. Ecosystem processes and human influences regulate streamflow response to climate change at long-term ecological research sites [J]. BioScience, 62 (4): 390–404. doi:10.1525/bio.2012.62.4.10.

    Keetch J J, Byram G M. 1968. A Drought Index for Forest Fire Control [R]. USDA Forest Service, Southeast Forest Experiment Station Research Paper SE-38. Asheville, NC, 35pp.

    Kottek M, Grieser J, Beck C, et al. 2006. World Map of the K?ppen-Geiger climate classification updated [J]. Meteor. Z., 15 (3): 259–263, doi:10.1127/0941-2948/2006/0130.

    Liu S M, Lu L, Mao D, et al. 2007. Evaluating parameterizations of aerodynamic resistance to heat transfer using field measurements [J]. Hydrol. Earth Syst. Sci., 11: 769–783, doi:10.5194/hess-11-769-2007.

    Liu Y Q. 2010. A numerical study on hydrological impacts of forest restoration in the southern United States [J]. Ecohydrology, 4 (2): 299–314, doi:10.1002/eco.178.

    Liu Y Q, Stanturf J, Lu H Q. 2008. Modeling the potential of the northern China forest shelterbelt in improving hydroclimate conditions [J]. J. Am. Water Resour. Assoc., 44 (5): 1176–1192, doi:10.1111/j.1752-1688.2008. 00240.x.

    Liu Y Q, Stanturf J, Goodrick S. 2010a. Trends in global wildfire potential in a changing climate [J]. Forest Ecol. Manag., 259 (4): 378–1127, 10.1016/j.foreco.2009.09.002.

    Liu Y Q, Stanturf J, Goodrick S. 2010b. Wildfire potential evaluation during a drought event with a regional climate model and NDVI [J]. Ecolog. Inform., 5 (5): 418–428, doi:10.1016/j.ecoinf.2010.04.001.

    Liu Y Q, Goodrick S L, Stanturf J A. 2013. Future U.S. wildfire potential trends projected using a dynamically downscaled climate change scenario [J]. Forest Ecol. Manag., 294: 120–135, doi:10.1016/j.foreco.2012.06.049.

    Liu Z Y, Notaro M, Kutzbach J, et al. 2006. Assessing global vegetationclimate feedbacks from the observation [J]. J. Climate, 19 (5): 787–814, doi:10.1175/JCLI3658.1.

    Ma D, Notaro M, Liu Z Y, et al. 2013. Simulated impacts of afforestation in East China monsoon region as modulated by ocean variability [J]. Climate Dyn., 41: 2439–2450, doi:10.1007/s00382-012-1592-9.

    Mearns L O, Arritt R, Biner S, et al. 2012. The North American regional climate change assessment program: Overview of phase I results [J]. Bull. Amer. Met. Soc., 93: 1337–1362, doi:10.1175/BAMS-D-11-00223.1.

    Meng X H, Evans J P, McCabe M F. 2014. The impact of observed vegetation changes on land–atmosphere feedbacks during drought [J]. J. Hydrometeor, 15: 759–776, doi:10.1175/JHM-D-13-0130.1.

    NWS (U. S. National Weather Service). 2006 [2015-6-21]. Drought [OL].. http://www.nws.noaa.gov/om/brochures/climate/Drought.pdf.

    Palmer W C. 1965. Meteorological drought [R]. Research Paper 45, U. S. Dept. of Commerce, 58pp. Seager R, Ting M F, Held I, et al. 2007. Model projections of an imminent transition to a more arid climate in southwestern North America [J]. Science, 316 (5828): 1181–1184, doi:10.1126/science.1139601.

    Teuling A J, van Loon A F, Seneviratne S I, et al. 2013. Evapotranspiration amplifies European summer drought [J]. Geophys. Res. Lett., 40: 2071–2075, doi:10.1002/grl.50495.

    Trenberth K E, Dai A, van der Schrier G, et al. 2014. Global warming and changes in drought [J]. Nature Climate Change, 4: 17–22, doi:10.1038/ NCLIMATE2067.

    Vicente-Serrano S M, Gouveia C, Camarero J J, et al. 2013. Response of vegetation to drought time-scales across global land biomes [J]. Proc. Natl Acad. Sci. USA, 110 (1): 52–57, doi:10.1073/pnas.1207068110.

    Wang L, Chen W. 2014. A CMIP5 multimodel projection of future temperature, precipitation, and climatological drought in China [J]. Int. J. Climatol., 34 (6): 2059–2078, doi:10.1002/joc.3822.

    Wu C Y, Chen J M. 2013. Diverse responses of vegetation production to interannual summer drought in North America [J]. Int. J. Appl. Earth Observ. Geoinform., 21: 1–6, doi:10.1016/j.jag.2012.08.001.

    吳紹洪, 戴爾阜, 黃玫, 等. 2007. 21世紀(jì)未來氣候變化情景(B2)下我國(guó)生態(tài)系統(tǒng)的脆弱性研究 [J]. 科學(xué)通報(bào), 52 (7): 811–817. Wu S H, Dai E F, Huang M, et al. 2007. Ecosystem vulnerability of China under B2 climate scenario in the 21st century [J]. Chin. Sci. Bull., 52 (10): 1379–1386, doi:10.1007/s11434-007-0197-x.

    Xanthopoulos G, Maheras G, Gouma V, et al. 2006. Is the Keetch–Byram drought index (KBDI) directly related to plant water stress? [J] Forest Ecol. Manag., 234 (Suppl. 1): S27.

    Yeh T -C, Wetherald R T, Manabe S. 1984. The effect of soil moisture on the short-term climate and hydrology change—A numerical experiment [J]. Mon. Wea. Rev., 112: 474–490, doi:10.1175/1520-0493(1984)112 <0474: TEOSMO>2.0.CO;2.

    Zargar A, Sadiq R, Naser B, et al. 2011. A review of drought indices [J]. Environ. Rev., 19 (NA): 333–349, doi:10.1139/a11-013.

    Zhang B Q, Wu P Y, Zhao X N, et al. 2013. Spatiotemporal analysis of climate variability (1971–2010) in spring and summer on the Loess Plateau, China [J]. Hydrological Processes, 28 (4): 1689–1702, doi:10.1002/hyp.9724.

    Zhang X Q, Lei Y C, Pang Y, et al. 2014. Tree mortality in response to climate change induced drought across Beijing, China [J]. Climatic Change, 124 (1–2): 179–190, doi:10.1007/s10584-014-1089-0.

    曾紅玲, 季勁鈞, 吳國(guó)雄. 2010. 全球植被分布對(duì)氣候影響的數(shù)值試驗(yàn)[J]. 大氣科學(xué), 34 (1): 1–11. Zeng H L, Ji J J, Wu G X. 2010. Numerical experiment of the influence of global vegetation distribution on climate [J]. Chinese J. Atmos. Sci., 34 (1): 1–11.

    Zeng N, Neelin J D, Lau K -M, et al. 1999. Enhancement of interdecadal climate variability in the Sahel by vegetation interaction [J]. Science, 286: 1537–1540, doi:10.1126/science.286.5444.1537.

    Zheng Y Q, Yu G, Qian Y F, et al. 2002. Simulations of regional climatic effects of vegetation change in China [J]. Q. J. R. Meteor. Soc., 128 (584): 2089–2114, doi:10.1256/003590002320603557.

    劉永強(qiáng). 2016. 植被對(duì)干旱趨勢(shì)的影響 [J]. 大氣科學(xué), 40 (1): 142?156. Liu Yongqiang. 2016. Impacts of vegetation on drought trends [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 40 (1): 142?156, doi:10.3878/j.issn.1006-9895.1508.15146.

    作者簡(jiǎn)介劉永強(qiáng),男,1959年出生, 博士, 研究員, 主要從事氣候與生態(tài)系統(tǒng)相互作用研究。E-mail: yliu@fs.fed.us

    收稿日期2015-03-08;網(wǎng)絡(luò)預(yù)出版日期 2015-08-06

    doi:10.3878/j.issn.1006-9895.1508.15146

    文章編號(hào)1006-9895(2016)01-0142-15

    中圖分類號(hào)P461

    文獻(xiàn)標(biāo)識(shí)碼A

    猜你喜歡
    干旱氣候變化植被
    基于植被復(fù)綠技術(shù)的孔植試驗(yàn)及應(yīng)用
    《應(yīng)對(duì)氣候變化報(bào)告(2022)》發(fā)布
    氣候變化與環(huán)保法官
    氣候變化:法官的作用
    綠色植被在溯溪旅游中的應(yīng)用
    基于距平的白城地區(qū)干旱時(shí)間分布特征分析
    臨夏地區(qū)干旱特征及干濕氣候區(qū)劃
    夏季高溫干旱時(shí)節(jié)高山蔬菜種植管理策略
    小說里的世界 虛擬里的真實(shí)
    基于原生植被的長(zhǎng)山群島植被退化分析
    日韩免费高清中文字幕av| 国产又色又爽无遮挡免| 亚洲丝袜综合中文字幕| 久久av网站| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 亚洲第一区二区三区不卡| 九九在线视频观看精品| 美女脱内裤让男人舔精品视频| 日韩欧美 国产精品| 日韩一本色道免费dvd| 国产熟女午夜一区二区三区 | 亚洲av免费高清在线观看| 午夜老司机福利剧场| 一本大道久久a久久精品| 亚洲精品国产av成人精品| 乱系列少妇在线播放| 日日爽夜夜爽网站| 我要看日韩黄色一级片| 人妻人人澡人人爽人人| 中文乱码字字幕精品一区二区三区| 久久人人爽人人片av| 中文字幕亚洲精品专区| 日本91视频免费播放| 久久狼人影院| 国产在线免费精品| freevideosex欧美| 国产成人91sexporn| 九色成人免费人妻av| 成人影院久久| av国产久精品久网站免费入址| 久久国产亚洲av麻豆专区| 老司机亚洲免费影院| 韩国av在线不卡| 乱人伦中国视频| 午夜影院在线不卡| 久久亚洲国产成人精品v| 欧美+日韩+精品| 三级经典国产精品| freevideosex欧美| 久久人人爽人人爽人人片va| 久久97久久精品| 一本一本综合久久| 亚洲不卡免费看| h视频一区二区三区| 看免费成人av毛片| 亚洲av日韩在线播放| 中文欧美无线码| 九九久久精品国产亚洲av麻豆| 91久久精品国产一区二区三区| 一级毛片久久久久久久久女| 人妻人人澡人人爽人人| 超碰97精品在线观看| 亚洲av在线观看美女高潮| 精品视频人人做人人爽| 插逼视频在线观看| 少妇人妻 视频| 日韩av不卡免费在线播放| 久久精品国产自在天天线| 国产精品人妻久久久影院| 91久久精品国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 亚洲精品国产av成人精品| 欧美3d第一页| 免费人成在线观看视频色| 国产成人精品福利久久| 亚洲激情五月婷婷啪啪| 色吧在线观看| 午夜免费鲁丝| 国产成人aa在线观看| 国产亚洲精品久久久com| 青春草国产在线视频| 精品亚洲乱码少妇综合久久| 欧美人与善性xxx| 免费大片黄手机在线观看| 国产免费一级a男人的天堂| 欧美日韩视频精品一区| 国产亚洲最大av| 少妇猛男粗大的猛烈进出视频| 国产淫片久久久久久久久| 在线亚洲精品国产二区图片欧美 | av不卡在线播放| 精品久久久精品久久久| 啦啦啦啦在线视频资源| 嘟嘟电影网在线观看| av专区在线播放| www.av在线官网国产| 亚洲真实伦在线观看| 交换朋友夫妻互换小说| 草草在线视频免费看| 男女国产视频网站| 国产精品久久久久成人av| 国产中年淑女户外野战色| 人妻 亚洲 视频| 女性生殖器流出的白浆| 一个人看视频在线观看www免费| 国产午夜精品久久久久久一区二区三区| 偷拍熟女少妇极品色| 香蕉精品网在线| 欧美激情极品国产一区二区三区 | 欧美日韩一区二区视频在线观看视频在线| 欧美精品国产亚洲| 亚洲精品日本国产第一区| 久久久久久久久久久丰满| 久久久久久久久久人人人人人人| 一本色道久久久久久精品综合| 亚洲,一卡二卡三卡| 女性被躁到高潮视频| 涩涩av久久男人的天堂| 少妇猛男粗大的猛烈进出视频| 91精品国产国语对白视频| 蜜臀久久99精品久久宅男| 欧美97在线视频| 久久久久久久大尺度免费视频| 久久精品国产a三级三级三级| 久久狼人影院| 少妇被粗大的猛进出69影院 | 久久久久久久久久久丰满| 中文字幕制服av| 男人狂女人下面高潮的视频| 成年av动漫网址| 伊人久久国产一区二区| 一个人免费看片子| 国产又色又爽无遮挡免| 91精品一卡2卡3卡4卡| 欧美高清成人免费视频www| 成人综合一区亚洲| 日韩一区二区三区影片| 一本一本综合久久| 性色avwww在线观看| 曰老女人黄片| 一本一本综合久久| 毛片一级片免费看久久久久| 久久热精品热| 亚洲欧美一区二区三区国产| 亚洲精品国产av成人精品| 亚洲成人手机| 免费观看无遮挡的男女| 91精品国产国语对白视频| 人人妻人人添人人爽欧美一区卜| av天堂久久9| 偷拍熟女少妇极品色| 在线观看美女被高潮喷水网站| 国产精品三级大全| 亚洲国产欧美在线一区| 国产精品秋霞免费鲁丝片| 婷婷色综合大香蕉| 一级毛片 在线播放| av一本久久久久| 青春草亚洲视频在线观看| 人人澡人人妻人| 亚洲国产精品一区三区| 国产无遮挡羞羞视频在线观看| 久久狼人影院| 欧美精品一区二区大全| 爱豆传媒免费全集在线观看| 成人无遮挡网站| 热re99久久精品国产66热6| 夜夜骑夜夜射夜夜干| 中文字幕免费在线视频6| 久久久精品免费免费高清| 成年美女黄网站色视频大全免费 | 久久久欧美国产精品| 亚洲第一区二区三区不卡| 国产成人免费观看mmmm| 国产欧美另类精品又又久久亚洲欧美| 最近手机中文字幕大全| 精品亚洲成a人片在线观看| 女人久久www免费人成看片| 99热国产这里只有精品6| kizo精华| 日本与韩国留学比较| 黑人猛操日本美女一级片| 91精品国产国语对白视频| 极品教师在线视频| 人人妻人人看人人澡| 天天躁夜夜躁狠狠久久av| 国产亚洲欧美精品永久| 99久久精品国产国产毛片| 亚洲精品一区蜜桃| 免费看不卡的av| 久久热精品热| 久久久久人妻精品一区果冻| 性色av一级| 午夜视频国产福利| 最新的欧美精品一区二区| 2018国产大陆天天弄谢| 日本欧美国产在线视频| 少妇人妻精品综合一区二区| 久久久精品免费免费高清| 纵有疾风起免费观看全集完整版| 777米奇影视久久| 婷婷色综合www| 好男人视频免费观看在线| 婷婷色av中文字幕| 久久久久久久亚洲中文字幕| 午夜日本视频在线| 18禁裸乳无遮挡动漫免费视频| 中文字幕av电影在线播放| 亚洲av男天堂| 一级二级三级毛片免费看| 极品人妻少妇av视频| 国产黄片视频在线免费观看| 久久久久久久大尺度免费视频| 欧美激情极品国产一区二区三区 | 久久久久久久精品精品| 国产一区有黄有色的免费视频| 中文字幕人妻丝袜制服| 日本免费在线观看一区| 国产一区亚洲一区在线观看| 嫩草影院入口| 乱人伦中国视频| 大香蕉久久网| 最黄视频免费看| 日韩伦理黄色片| av卡一久久| 久久精品国产a三级三级三级| 国产高清不卡午夜福利| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| 乱码一卡2卡4卡精品| 久久久久久久久久人人人人人人| 日本黄大片高清| 亚洲国产欧美在线一区| 韩国av在线不卡| 亚洲,欧美,日韩| 久久99精品国语久久久| 久久99热这里只频精品6学生| 美女视频免费永久观看网站| 精品酒店卫生间| 国内揄拍国产精品人妻在线| 夜夜骑夜夜射夜夜干| 乱人伦中国视频| 国产精品久久久久久久电影| √禁漫天堂资源中文www| 麻豆成人午夜福利视频| 国产一区二区三区综合在线观看 | 亚洲怡红院男人天堂| 亚洲一区二区三区欧美精品| 一级片'在线观看视频| 女性被躁到高潮视频| 少妇的逼水好多| 欧美亚洲 丝袜 人妻 在线| 亚洲欧洲精品一区二区精品久久久 | 日韩欧美一区视频在线观看 | 亚洲欧洲精品一区二区精品久久久 | 黄色视频在线播放观看不卡| 丝袜喷水一区| 大香蕉久久网| 国产精品偷伦视频观看了| 色哟哟·www| 免费黄色在线免费观看| 欧美+日韩+精品| 3wmmmm亚洲av在线观看| 王馨瑶露胸无遮挡在线观看| 久久久国产欧美日韩av| 国产女主播在线喷水免费视频网站| 少妇的逼好多水| 免费在线观看成人毛片| 国产伦理片在线播放av一区| 国产 一区精品| 日日摸夜夜添夜夜爱| 免费黄频网站在线观看国产| 欧美三级亚洲精品| 国产亚洲午夜精品一区二区久久| 日韩一区二区三区影片| 五月开心婷婷网| 女人精品久久久久毛片| 精品一区二区三区视频在线| 十八禁网站网址无遮挡 | 国产亚洲av片在线观看秒播厂| 高清黄色对白视频在线免费看 | 亚洲丝袜综合中文字幕| 自拍偷自拍亚洲精品老妇| 91久久精品国产一区二区三区| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 国产极品粉嫩免费观看在线 | 国产日韩一区二区三区精品不卡 | 亚洲av免费高清在线观看| 午夜视频国产福利| 国产真实伦视频高清在线观看| 性色avwww在线观看| 日日摸夜夜添夜夜添av毛片| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 夜夜爽夜夜爽视频| 国产永久视频网站| 人人妻人人添人人爽欧美一区卜| 国产在线一区二区三区精| 日韩欧美 国产精品| 18禁在线播放成人免费| 纯流量卡能插随身wifi吗| 精品国产露脸久久av麻豆| 成人18禁高潮啪啪吃奶动态图 | 五月天丁香电影| 国产伦精品一区二区三区四那| 日本av免费视频播放| 国产精品一区www在线观看| 高清午夜精品一区二区三区| 日韩欧美精品免费久久| h日本视频在线播放| 天堂8中文在线网| 国产欧美日韩精品一区二区| 亚洲怡红院男人天堂| 美女视频免费永久观看网站| 如日韩欧美国产精品一区二区三区 | 91精品国产国语对白视频| 人妻一区二区av| 日本vs欧美在线观看视频 | 99热国产这里只有精品6| 啦啦啦啦在线视频资源| 国产精品无大码| 99热全是精品| 亚洲精品自拍成人| 水蜜桃什么品种好| 亚洲精品亚洲一区二区| a级毛片免费高清观看在线播放| 在线观看www视频免费| 五月开心婷婷网| av又黄又爽大尺度在线免费看| 丰满少妇做爰视频| 偷拍熟女少妇极品色| 欧美 日韩 精品 国产| 一个人免费看片子| 熟女电影av网| 久久99蜜桃精品久久| 午夜91福利影院| 少妇的逼水好多| 亚洲综合色惰| 欧美日韩av久久| 亚洲第一av免费看| 在线观看免费日韩欧美大片 | 国产精品免费大片| 国产精品人妻久久久影院| 亚洲精品中文字幕在线视频 | 国产高清有码在线观看视频| 亚洲丝袜综合中文字幕| 哪个播放器可以免费观看大片| 欧美日韩国产mv在线观看视频| 午夜免费观看性视频| 亚洲美女黄色视频免费看| 91精品国产国语对白视频| 中文字幕久久专区| 哪个播放器可以免费观看大片| 少妇裸体淫交视频免费看高清| 久久精品国产亚洲网站| 日韩在线高清观看一区二区三区| 午夜av观看不卡| av国产精品久久久久影院| av在线app专区| 午夜福利网站1000一区二区三区| 久久午夜综合久久蜜桃| 日韩伦理黄色片| av国产精品久久久久影院| 涩涩av久久男人的天堂| 国产女主播在线喷水免费视频网站| 亚洲情色 制服丝袜| 日本黄色片子视频| 国产在线免费精品| 黑人高潮一二区| 在现免费观看毛片| 内射极品少妇av片p| 日本wwww免费看| 老司机影院成人| 中文资源天堂在线| 99久久精品热视频| av卡一久久| 久久久国产欧美日韩av| 精品人妻偷拍中文字幕| 成人18禁高潮啪啪吃奶动态图 | 国产成人精品无人区| 最近最新中文字幕免费大全7| 男人和女人高潮做爰伦理| 少妇被粗大的猛进出69影院 | 激情五月婷婷亚洲| 乱系列少妇在线播放| 久久av网站| 亚洲美女视频黄频| 97精品久久久久久久久久精品| 九草在线视频观看| 夫妻午夜视频| 亚洲美女搞黄在线观看| 亚洲精品亚洲一区二区| 人妻人人澡人人爽人人| 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91 | 99精国产麻豆久久婷婷| 一级毛片久久久久久久久女| 国产色婷婷99| 最黄视频免费看| 97精品久久久久久久久久精品| 在线天堂最新版资源| 国产精品女同一区二区软件| 又黄又爽又刺激的免费视频.| 久久99蜜桃精品久久| 卡戴珊不雅视频在线播放| av视频免费观看在线观看| 18禁动态无遮挡网站| 在线观看人妻少妇| 国产免费一级a男人的天堂| 一级二级三级毛片免费看| 精品久久久噜噜| 26uuu在线亚洲综合色| 欧美xxxx性猛交bbbb| 三级经典国产精品| 我要看日韩黄色一级片| 成人美女网站在线观看视频| 最新中文字幕久久久久| 免费黄色在线免费观看| 免费观看在线日韩| 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| 欧美精品人与动牲交sv欧美| 亚洲精品国产av成人精品| www.av在线官网国产| 搡老乐熟女国产| 一级爰片在线观看| 欧美bdsm另类| 国产乱来视频区| 男女免费视频国产| videossex国产| 成年人免费黄色播放视频 | 亚洲国产精品成人久久小说| 只有这里有精品99| 国产欧美另类精品又又久久亚洲欧美| 丁香六月天网| 国产黄色免费在线视频| 国产美女午夜福利| 亚洲精品456在线播放app| 国产一级毛片在线| 国产成人一区二区在线| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久av不卡| 男女无遮挡免费网站观看| 日韩制服骚丝袜av| 精品人妻熟女av久视频| 亚洲av二区三区四区| av在线播放精品| 欧美+日韩+精品| 国产日韩一区二区三区精品不卡 | 人妻少妇偷人精品九色| 插逼视频在线观看| 毛片一级片免费看久久久久| 久久精品国产亚洲av涩爱| 久久久国产一区二区| 少妇 在线观看| 我的老师免费观看完整版| 久久综合国产亚洲精品| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | 国产爽快片一区二区三区| 日日撸夜夜添| 亚洲人成网站在线观看播放| 日日摸夜夜添夜夜爱| 欧美日韩一区二区视频在线观看视频在线| 最近的中文字幕免费完整| 国产精品成人在线| 高清在线视频一区二区三区| 亚洲欧洲日产国产| 精品一品国产午夜福利视频| 国产国拍精品亚洲av在线观看| 中文字幕久久专区| 国产精品女同一区二区软件| av福利片在线观看| 热99国产精品久久久久久7| 久久人妻熟女aⅴ| 欧美xxxx性猛交bbbb| 免费观看a级毛片全部| 国产深夜福利视频在线观看| 精品久久久久久久久av| 日日摸夜夜添夜夜添av毛片| 久久午夜福利片| 国产精品一区二区性色av| a级毛色黄片| 精品国产一区二区久久| 在线观看免费视频网站a站| 色视频在线一区二区三区| 黑丝袜美女国产一区| 乱系列少妇在线播放| 成人黄色视频免费在线看| 中文在线观看免费www的网站| 多毛熟女@视频| 亚洲在久久综合| 久久午夜福利片| 美女大奶头黄色视频| 亚洲欧美日韩另类电影网站| 一区二区av电影网| a级一级毛片免费在线观看| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 国产精品国产三级国产专区5o| 亚洲欧洲国产日韩| 美女中出高潮动态图| 午夜免费鲁丝| 日本-黄色视频高清免费观看| a 毛片基地| 最黄视频免费看| 久久久国产欧美日韩av| 亚洲自偷自拍三级| 99久久精品一区二区三区| av在线播放精品| 国产免费福利视频在线观看| 在线观看www视频免费| 日日撸夜夜添| 尾随美女入室| videossex国产| 亚洲精品乱码久久久久久按摩| 久久久久久人妻| 免费久久久久久久精品成人欧美视频 | 国产精品久久久久久av不卡| 最近中文字幕2019免费版| 国产中年淑女户外野战色| 国产伦精品一区二区三区视频9| h视频一区二区三区| 亚洲欧美日韩卡通动漫| 中文天堂在线官网| 国产91av在线免费观看| 国产精品欧美亚洲77777| 精品酒店卫生间| a 毛片基地| 精品久久久久久久久亚洲| 亚洲av成人精品一区久久| 欧美日韩视频高清一区二区三区二| 国产熟女欧美一区二区| 久久精品久久久久久久性| 美女国产视频在线观看| 国产一区有黄有色的免费视频| 一级毛片久久久久久久久女| 天堂中文最新版在线下载| 国产精品久久久久久精品古装| 日韩三级伦理在线观看| 伊人久久精品亚洲午夜| 麻豆成人午夜福利视频| 国内揄拍国产精品人妻在线| av视频免费观看在线观看| 久久6这里有精品| 人体艺术视频欧美日本| 欧美3d第一页| 欧美精品国产亚洲| 岛国毛片在线播放| 热99国产精品久久久久久7| 交换朋友夫妻互换小说| 一级黄片播放器| 乱系列少妇在线播放| 成人免费观看视频高清| 国产成人a∨麻豆精品| 亚洲成人av在线免费| 午夜91福利影院| 天堂中文最新版在线下载| 亚洲精品日本国产第一区| 久久久亚洲精品成人影院| 99热全是精品| 日本av免费视频播放| 国产精品蜜桃在线观看| 日韩伦理黄色片| 又爽又黄a免费视频| 最近最新中文字幕免费大全7| 老女人水多毛片| 在线观看三级黄色| 自线自在国产av| 久久女婷五月综合色啪小说| 男男h啪啪无遮挡| 精品熟女少妇av免费看| 亚洲精品第二区| 一区在线观看完整版| 国产69精品久久久久777片| 国产成人一区二区在线| 好男人视频免费观看在线| 国产女主播在线喷水免费视频网站| 国产日韩欧美视频二区| 深夜a级毛片| 高清毛片免费看| 99久久中文字幕三级久久日本| 欧美日韩亚洲高清精品| 日韩中文字幕视频在线看片| 男人添女人高潮全过程视频| 免费大片18禁| 国产高清不卡午夜福利| 人妻少妇偷人精品九色| 深夜a级毛片| 极品少妇高潮喷水抽搐| 多毛熟女@视频| 久久久久国产网址| 国产在视频线精品| 我的老师免费观看完整版| 丰满饥渴人妻一区二区三| 国产一区二区在线观看av| 久热这里只有精品99| 91在线精品国自产拍蜜月| 黑丝袜美女国产一区| 日韩 亚洲 欧美在线| 如日韩欧美国产精品一区二区三区 | 国产成人91sexporn| 国产伦精品一区二区三区四那| 嫩草影院新地址| 亚洲电影在线观看av| a级毛色黄片| 国产精品无大码| 如何舔出高潮| 久久久久久久久大av| 少妇的逼水好多| 嫩草影院入口| 久久久亚洲精品成人影院| 午夜精品国产一区二区电影| 国产亚洲91精品色在线| 内射极品少妇av片p| 高清欧美精品videossex| 国产精品福利在线免费观看| 极品少妇高潮喷水抽搐| 91久久精品国产一区二区三区| 少妇高潮的动态图| 亚洲精品久久午夜乱码|