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

    厄爾尼諾次年夏季中國地區(qū)云的變化特征及其與降水的關(guān)系*

    2022-11-02 07:34:30彭艷玉郜倩倩
    氣象學(xué)報(bào) 2022年5期

    彭艷玉 劉 煜 郜倩倩

    災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,中國氣象局大氣化學(xué)重點(diǎn)實(shí)驗(yàn)室,中國氣象科學(xué)研究院,北京,100081

    1 引言

    厄爾尼諾事件是指赤道中、東太平洋海表溫度大范圍持續(xù)異常偏高的熱帶海-氣相互作用的一種氣候現(xiàn)象,太平洋東部海溫異常升高并向西擴(kuò)展稱為東部型厄爾尼諾,赤道中太平洋海溫異常升高并向東擴(kuò)展稱為中部型厄爾尼諾(任福民等,2012)。中部型厄爾尼諾對南美、北美西海岸、甚至日本和新西蘭氣候的影響可能會與東部型厄爾尼諾的影響 完 全 相 反(Larkin,et al,2005;Ashok,et al,2007)。中國位于東亞季風(fēng)區(qū)、太平洋西海岸,鄰近厄爾尼諾發(fā)生區(qū),氣候必然會受到厄爾尼諾現(xiàn)象的影響。降水是氣候最基本的組成因素之一,許多研究表明:東部型厄爾尼諾事件次年夏季中國形成南、北兩條異常雨帶,江淮地區(qū)少雨(Huang,et al,1989;陳文,2002;Feng,et al,2011;王欽等,2012;吳萍等,2017);中部型厄爾尼諾事件次年夏季中國形成一條異常雨帶,長江以南少雨(Feng,et al,2011;王欽等,2012;吳萍等,2017)。不同類型厄爾尼諾事件的大氣環(huán)流特征以及氣候影響存在明 顯 差 異(Yuan,et al,2012a,2012b;袁 媛 等,2012;Chen,et al,2014;李麗平等,2015;汪婉婷等,2018)。

    云是氣候系統(tǒng)中最復(fù)雜的因素之一,在氣候系統(tǒng)的能量收支和水分循環(huán)過程中起著重要作用,云相關(guān)物理量的任何變化都可能對全球氣候產(chǎn)生重大影響。中國地區(qū)多年平均總云量分布的總體形勢是南多北少、東多西少??傇屏?、高云量和低云量呈減少趨勢,中云量稍有增加(Kaiser,1998;丁守國等,2004,2005;吳澗等,2011)。但是,段皎等(2011)的結(jié)果則顯示中國大部分地區(qū)總云量沒有顯著的變化趨勢,在華南地區(qū)和渤海灣以及西北部分地區(qū)有增多趨勢,青藏高原中部有減少趨勢;中國大部分地區(qū)高云量沒有明顯的變化趨勢;華南和東北地區(qū)以及青藏高原西部中云量有增多趨勢。東亞季風(fēng)區(qū)低云發(fā)生頻率最大,夏季深對流云頻率增加明顯(劉建軍等,2017)。

    在云和降水的關(guān)系方面,已有較多研究成果。中國年平均總云量與年降水量存在良好的線性正相關(guān)關(guān)系(曾昭美等,1993)。對流云多產(chǎn)生陣性降水,時(shí)間短強(qiáng)度大,層狀云可產(chǎn)生連續(xù)性降水,雨強(qiáng)相對較小但持續(xù)時(shí)間長(胡亮等,2010;杜振彩等,2011);對流云降水主要來自深對流云和卷層云;層云降水主要來自高層云和層積云(李昀英等,2015)。西南地區(qū)夏季雨層云、高層云和深對流云與雨量和雨日有著較好的一致關(guān)系(李躍清等,2014);西北地區(qū)云水路徑值較大的層狀云類的云量多寡與降水多寡一致,積狀云類和層積云類云量多少與降水不存在一定的關(guān)系(陳勇航等,2005)。也有一些學(xué)者利用風(fēng)云氣象衛(wèi)星相關(guān)數(shù)據(jù)針對中國部分地區(qū)探討了云特征參數(shù)與降水量的相關(guān),認(rèn)為有些云結(jié)構(gòu)特征參數(shù)對降水有著很好的指示作用(周毓荃等,2011;張中波等,2017;王磊等,2019)。

    中國云量年際振蕩的主要特征與ENSO相同,ENSO年份中國云量少,降水也偏少;全國普遍存在的ENSO尺度振蕩主頻率有由南向北、由東向西推遲的傾向(曾昭美等,1993);ENSO通過影響高空大氣環(huán)流狀況,進(jìn)而影響中國各地云量變化而最終影響各地的降水量(盧愛剛等,2009)。李月洪(1987)在其研究中提到:厄爾尼諾發(fā)生年冬季西太平洋赤道地區(qū)150°E以東云量增加,以西云量減少。云的產(chǎn)生條件和海溫有密切的關(guān)系,在研究海洋和大氣相互作用時(shí),探究云與海洋的關(guān)系是十分重要的。ENSO包括不同類型的厄爾尼諾事件和不同類型的拉尼娜事件,且事件的影響包括事件發(fā)生當(dāng)年和事件發(fā)生次年,不同季節(jié)影響也有很大的不同,以上關(guān)于海溫異常對云影響的相關(guān)研究均過于簡單、籠統(tǒng)。眾多已有研究成果表明,厄爾尼諾事件對中國次年夏季風(fēng)進(jìn)程和夏季降水有著極其重要的影響,而云的形成是產(chǎn)生降水的必要條件,但將厄爾尼諾事件與云的變化以及云和降水關(guān)系聯(lián)系起來的研究成果相對缺乏。如今全球變暖仍在持續(xù),近百年中國變暖趨勢高于全球平均水平(嚴(yán)中偉等,2020),全球變暖更穩(wěn)定地體現(xiàn)在海洋(Cheng,et al,2019),在此背景下將厄爾尼諾、降水和云三者聯(lián)系起來分析中國氣候?qū)Σ煌愋投驙柲嶂Z事件的響應(yīng)規(guī)律和敏感性,對于了解中國各地氣候特征與海溫異常的關(guān)系、全面認(rèn)識中國氣候變化、改進(jìn)天氣預(yù)報(bào)和短期氣候預(yù)測、更好地應(yīng)對全球變化等都具有重要意義。因此,本研究期望能夠?qū)⒍驙柲嶂Z和云、降水聯(lián)系起來研究,以探討兩類厄爾尼諾事件次年夏季中國云相關(guān)參數(shù)的變化特征及其與降水變化的聯(lián)系。

    2 資料和方法

    2.1 資料簡介

    地面常規(guī)觀測資料的觀測時(shí)間序列長、局地代表性強(qiáng),但測站分布極不均勻;衛(wèi)星遙感資料具有覆蓋范圍廣、時(shí)間分辨率高、能夠反演云的多種物理量等特點(diǎn),但由于對云識別技術(shù)還不完善,對云的分層描述和低云的準(zhǔn)確識別還不理想(王旻燕等,2009)。并且,地基與天基的云觀測視角不同,反映云的特性也不同,如:深對流云,地基觀測根據(jù)云底高度,其為低云,衛(wèi)星觀測根據(jù)云頂高度,其為高云??紤]到每種數(shù)據(jù)都有其局限性和視角的不同,文中主要用到中國地面臺站常規(guī)觀測的降水和云量數(shù)據(jù)、氣象信息綜合分析處理系統(tǒng)(MICAPS)降水和云相關(guān)數(shù)據(jù)、國際衛(wèi)星云氣候?qū)W計(jì)劃(ISCCP)的D2數(shù)據(jù)、中分辨率成像光譜儀(MODIS)數(shù)據(jù)(MOD08_M3)和熱帶測雨衛(wèi)星(TRMM)3A25數(shù)據(jù)等數(shù)據(jù)集的夏季(6—8月)數(shù)據(jù)來做分析,以盡可能提高所得結(jié)果的可靠性。時(shí)間長、包含個(gè)例多的站點(diǎn)數(shù)據(jù)用于分析厄爾尼諾次年夏季降水和云量的變化特征以及二者的關(guān)系;含較多氣象觀測要素的MICAPS數(shù)據(jù)用于補(bǔ)充低云狀出現(xiàn)頻率的變化;含多種云相關(guān)變量的ISCCP和MODIS數(shù)據(jù)用于分析云頂氣壓、溫度和光學(xué)厚度等特性與降水異常的關(guān)系;含對流降水和層云降水的TRMM衛(wèi)星數(shù)據(jù)用于分析厄爾尼諾事件對不同類型降水的影響。用到的5種資料簡介如下:

    (1)中國國家級地面氣象站1961—2010年逐日降水和云量數(shù)據(jù)。數(shù)據(jù)集由中國氣象局國家氣象信息中心提供,原始數(shù)據(jù)經(jīng)過嚴(yán)格的質(zhì)量檢測與控制。(2)MICAPS(Meteorological Information Combine Analysis and Process System)2000—2013年第一類地面全要素填圖數(shù)據(jù)。該數(shù)據(jù)集是氣象觀測臺站每3 h的原始觀測數(shù)據(jù),包含云狀等數(shù)據(jù)集(1)中沒有的信息。本研究主要用到降水、云量和低云狀等要素,用于數(shù)據(jù)集(1)所得結(jié)果的補(bǔ)充。(3)ISCCP(International Satellite Cloud Climatology Project)1984—2009年的D2數(shù)據(jù)集。D2數(shù)據(jù)集的空間分辨率為2.5°×2.5°,共有130個(gè)與云相關(guān)的物理量,含3 h和日平均數(shù)據(jù),本研究主要用到日平均的云量、云頂氣壓、云頂溫度、光學(xué)厚度和云水路徑等氣象要素。(4)MODIS(Moderate-Resolution Imaging Spectroradiometer)2000—2017年的MOD08_M3數(shù)據(jù)集。MOD08_M3是MODIS三級全球1°×1°格點(diǎn)的月平均大氣標(biāo)準(zhǔn)產(chǎn)品,是全球氣溶膠、水汽和云月合成產(chǎn)品。(5)TRMM(Tropical Rainfall Measuring Mission satellite)1998—2014年 的3A25降水 數(shù) 據(jù)。3A25產(chǎn)品是TRMM星載降水雷達(dá)(PR)格點(diǎn)化月平均降水資料3級產(chǎn)品。分辨率是0.5°×0.5°,數(shù)據(jù)集經(jīng)過嚴(yán)格的質(zhì)量控制。

    2.2 方法

    根據(jù)2017年5月中國氣象局頒布的氣象行業(yè)標(biāo)準(zhǔn)《厄爾尼諾/拉尼娜事件判別方法》(QB/T33666-2017)將各數(shù)據(jù)相應(yīng)時(shí)間長度范圍內(nèi)的厄爾尼諾事件進(jìn)行分類。以東部型(EP)和中部型(CP)來區(qū)分不同分布型的厄爾尼諾事件。本研究將厄爾尼諾發(fā)生的后一年稱為厄爾尼諾次年。各類數(shù)據(jù)時(shí)段內(nèi)對應(yīng)的兩類事件年份和數(shù)據(jù)使用詳細(xì)情況說明見表1。

    表1 各數(shù)據(jù)集的選用說明Table 1 Instruction for the selection of datasets

    以厄爾尼諾次年夏季氣象要素合成值相對于夏季氣候平均態(tài)的距平百分率來表示厄爾尼諾事件對夏季中國地區(qū)降水和云等各種氣象要素的影響。用到的研究方法(黃嘉佑,2004)主要有合成分析、差異分析(距平:數(shù)據(jù)相對于氣候平均的差值,距平百分率:距平相對于氣候平均的百分比)、相關(guān)分析(線性相關(guān))、回歸分析(線性回歸),方法較簡單常見。差異分析的顯著性檢驗(yàn)用t檢驗(yàn)。

    3 結(jié)果分析

    3.1 厄爾尼諾對云量等多種云相關(guān)物理量的影響

    地面臺站數(shù)據(jù)經(jīng)過缺測率和遷址條件控制后一共選取1678個(gè)符合條件的站點(diǎn)(圖1),這些站點(diǎn)大部分位于中國的季風(fēng)氣候區(qū),正是本研究分析的主要區(qū)域。厄爾尼諾事件次年夏季總降水異常(圖2)表明,東部型事件次年夏季中國地區(qū)形成南北兩條異常雨帶(圖2a),長江流域、內(nèi)蒙古、華北和東北西部地區(qū)多雨;江淮地區(qū)少雨。中部型事件次年夏季中國地區(qū)形成一條異常雨帶(圖2b),華北南部、東北南部、江淮地區(qū)以及青藏高原等地區(qū)多雨;內(nèi)蒙古、黑龍江和長江以南地區(qū)少雨。該分析結(jié)果與吳萍等(2017)的研究結(jié)果基本一致,在西北地區(qū)和雨帶范圍上存在一些差異,可能是由于選取的氣候背景場和站點(diǎn)不一致所引起的。白天和夜間降水(圖略)與總降水的距平百分率分布形勢基本相同。

    圖1 常規(guī)站點(diǎn)數(shù)據(jù)進(jìn)行條件控制后選取的站點(diǎn)分布情況Fig.1 Distribution of stations selected after conditional control of conventional station dataset

    圖2 厄爾尼諾(a.東部型,b.中部型)事件次年夏季累計(jì)降水相對于氣候平均值的距平百分率(單位:%,站點(diǎn)數(shù)據(jù);通過置信度為95%顯著性檢驗(yàn)的站點(diǎn)以*標(biāo)注)Fig.2 Percentage anomalies(unit:%,station dataset)of cumulative precipitation in the El Ni?o(a.EPEl Ni?o,b.CP El Ni?o)decaying summers relative to the climatic average(stations with values passing significance test at the 95%confidence level are marked with*)

    在分析云的數(shù)據(jù)之前,對資料進(jìn)行了相互對比和驗(yàn)證,資料對比所用的時(shí)段和氣象要素選為所對比資料共有的時(shí)間長度(6—8月)和氣象要素。常規(guī)站點(diǎn)、ISCCP資料、MODIS的數(shù)據(jù)總云量分布形勢是一致的(圖略),高值中心均位于西南地區(qū),由南向北遞減,低值中心位于西北地區(qū);在黃河以南、華北和東北大部分地區(qū),三類數(shù)據(jù)總云量差異較小,基本在6%以內(nèi);在內(nèi)蒙古和西北地區(qū)三類數(shù)據(jù)的差異較明顯,ISCCP總云量比站點(diǎn)數(shù)據(jù)多6%—12%,MODIS總云量比站點(diǎn)數(shù)據(jù)少6%—12%。站點(diǎn)數(shù)據(jù)和MICAPS數(shù)據(jù)的降水和云量分布形勢基本吻合(圖略),在長江以南地區(qū)降水MICAPS數(shù)據(jù)略小于站點(diǎn)數(shù)據(jù);總云量MICAPS數(shù)據(jù)整體略小于站點(diǎn)數(shù)據(jù),而低云量MICAPS數(shù)據(jù)整體略大于站點(diǎn)數(shù)據(jù)。總體而言,各種數(shù)據(jù)之間的偏差主要在西部地區(qū),且數(shù)值偏差較小,比如總云量的偏差基本在1成云量以內(nèi),數(shù)據(jù)質(zhì)量好,具有較高的使用價(jià)值,可用于云的變化特征及其與降水關(guān)系的研究。

    3.1.1 站點(diǎn)資料:厄爾尼諾事件對總云量和低云量的影響

    厄爾尼諾事件次年夏季日平均云量合成值相對于氣候平均值的異常(圖3)表明,東部型事件對應(yīng)的白天總云量(圖3a)在大部分地區(qū)呈增多趨勢,增多大值區(qū)在華北北部和內(nèi)蒙古中部,增加幅度超過5%;在江漢、江淮和江南地區(qū)是減少的,減少大值區(qū)位于江淮地區(qū),減小幅度超過5%。白天低云量(圖3c)在黃河流域及其以北和西北地區(qū)東部有所增多,黃河中游地區(qū)增多最明顯,幅度超過20%;在西南、江淮和長江以南地區(qū)是減少的,減少幅度基本不超過8%。夜間低云量(圖3e)與白天基本一致,但在東北地區(qū)夜間低云量是減少的。中部型厄爾尼諾事件對應(yīng)的白天總云量(圖3b)在東北、華北、長江—黃河之間的區(qū)域呈增多趨勢,增多極大值區(qū)在華北東部沿海地區(qū),增幅在5%左右;在內(nèi)蒙古、江南和華南東部地區(qū)是減少的,減少幅度在4%以內(nèi)。低云量的異常分布(圖3d、f)大致同總云量,但低云量異常增多大值區(qū)位于華北南部,且幅度比總云量略大,接近10%。

    圖3 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季日平均云量(a、b.白天總云量,c、d.白天低云量,e、f.夜間低云量)的距平百分率(單位:%,站點(diǎn)數(shù)據(jù))Fig.3 Percentage anomalies(unit:%,station dataset) of daily average cloud amount in the El Ni?o decaying summers (a,c,e.EPEl Ni?o,b,d,f.CPEl Ni?o;a,b.total cloud amount at daytime,c,d.low cloud amount at daytime,e,f.low cloud amount at nighttime)

    3.1.2 ISCCP資料:厄爾尼諾事件對多種云量的影響

    厄爾尼諾事件次年站點(diǎn)數(shù)據(jù)的夏季降水的異常分布(氣候平均態(tài)是ISCCP數(shù)據(jù)對應(yīng)的時(shí)段1984—2009年)中,兩種類型事件對應(yīng)的降水異常分布形勢(圖略)與圖2基本一致,不過東部型厄爾尼諾事件對應(yīng)的南、北雨帶分界不明顯,北部異常雨區(qū)東邊界線西移,華北東部地區(qū)降水變化幅度不大,降水異常大值區(qū)在黃河中游流域。

    圖4是ISCCP數(shù)據(jù)中厄爾尼諾事件次年夏季云量的異常。由于衛(wèi)星觀測的方向是由上向下,故而其對低云量的觀測結(jié)果不甚理想(低云量和中云量圖略)。東部型厄爾尼諾對應(yīng)的總云量(圖4a)在西北、內(nèi)蒙古和西南地區(qū)是增多的;高云量(圖4c)在內(nèi)蒙古—華北以西和長江中上游流域及其以北是增多的;卷層云云量(圖4e)增減分界線較高云量更偏西;除了東北和華南南部部分地區(qū),中國東部地區(qū)深對流云云量都是增多的(圖4g)。中部型厄爾尼諾對應(yīng)的云量異常分布與降水(圖略)較為一致,但云量增多的范圍比降水增加范圍更大??傇屏浚▓D4b)在全國大部分地區(qū)都是增多的;高云量(圖4d)在長江以北地區(qū)都是增多的;卷層云云量(圖4f)在青海、四川和108°E以東是增多的;深對流云云量(圖4h)在中國西部和東部的長江—黃河之間是增多的。

    圖4 厄爾尼諾(a、c、e、g.東部型,b、d、f、h.中部型)事件次年夏季云量(a、b.總云量,c、d.高云量,e、f.卷層云云量,g、h.深對流云云量)的距平百分率(單位:%,ISCCP數(shù)據(jù))Fig.4 Percentage anomalies (unit:%,ISCCPdataset)of cloud amount in the El Ni?o decaying summers(a,c,e,g.EP El Ni?o,b,d,f,h.CPEl Ni?o;a,b.total cloud amount,c,d.high cloud amount,e,f.cirrostratus cloud amount,g,h.deep convective cloud amount)

    3.1.3 MODIS數(shù)據(jù):厄爾尼諾事件對多種云相關(guān)物理量的影響

    厄爾尼諾事件次年站點(diǎn)數(shù)據(jù)的夏季降水的異常情況(氣候平均態(tài)是MODIS數(shù)據(jù)對應(yīng)的時(shí)段2000—2017年)中,兩種類型事件對應(yīng)的降水異常分布形勢(圖略)與圖2存在差異,東部型厄爾尼諾事件仍然對應(yīng)南、北兩條雨帶(圖略),但在長江—黃河中下游流域南北雨帶連在一起,黃淮、江漢和江淮地區(qū)的降水是增加的,北部異常雨區(qū)擴(kuò)展到整個(gè)新疆和西藏西部,南部雨帶降水增加的幅度也更大;中部型厄爾尼諾對應(yīng)的雨帶分布(圖略)差異不大,只是東北雨區(qū)擴(kuò)展到整個(gè)東北地區(qū)。由于MODIS數(shù)據(jù)長度只有17 a,東部型和中部型事件個(gè)例分別為2個(gè)和3個(gè),降水分布的差異來源于個(gè)例數(shù)量少。

    從MODIS數(shù)據(jù)厄爾尼諾事件次年夏季云量的異常分布(圖5)可見,東部型厄爾尼諾對應(yīng)的總云量(圖5a)、高云量(圖5c)和卷云量(圖5e)也分為南、北兩條異常增多云帶,但在長江—黃河中下游流域之間云量是減少的,南異常云帶比南異常雨帶位置偏南,基本位于長江以南。中部型厄爾尼諾對應(yīng)的總云量(圖5b)、高云量(圖5d)和卷云量(圖5f)增多區(qū)分布在長江與黃河之間、西北、西南和東北地區(qū),幾乎與降水異常增加形勢(圖略)一致,但云南除外。圖6是MODIS數(shù)據(jù)厄爾尼諾事件次年夏季與總云相關(guān)的物理量的異常分布。東部型厄爾尼諾對應(yīng)的總云云頂氣壓(圖6a)和云頂溫度(圖6c)的異常分布形勢與降水異常相反,在降水增加的區(qū)域上云頂氣壓和云頂溫度是下降的;光學(xué)厚度(圖6e)與降水距平百分率增加區(qū)雖然存在較大差異,但在降水異常增加大值區(qū)光學(xué)厚度是增大的。中部型厄爾尼諾對應(yīng)的總云云頂氣壓(圖6b)和云頂溫度(圖6d)與降水的距平百分率分布形勢都呈大致相反趨勢,即:云頂氣壓減小和云頂溫度下降的區(qū)域?qū)?yīng)降水增加,云頂氣壓增大和云頂溫度升高的區(qū)域?qū)?yīng)降水減少;光學(xué)厚度(圖6f)和降水異常分布差異較小且正大值區(qū)基本重合。夏季平均總云不同相態(tài)云水路徑的異常分布形勢(圖略)基本與光學(xué)厚度一致。也就是說在降水異常增加區(qū)域,不僅總云云量增加,光學(xué)厚度增大,云水含量升高,而且云頂氣壓和云頂溫度下降,對流活動(dòng)增強(qiáng),云頂高度增高,云層增厚。

    圖5 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季云量(a、b.總云量,c、d.高云量,e、f.卷云云量)的距平百分率(單位:%,MODIS數(shù)據(jù))Fig.5 Percentage anomalies(unit:%,MODIS dataset)of cloud amount in the El Ni?o decaying summers(a,c,e.EP El Ni?o,b,d,f.CPEl Ni?o;a,b.total cloud amount,c,d.high cloud amount,e,f.cirrus cloud amount)

    圖6 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季與總云相關(guān)的物理量(a、b.云頂氣壓,c、d. 云頂溫度,e、f.光學(xué)厚度)距平百分率(單位:%,MODIS數(shù)據(jù))Fig.6 Percentage anomalies(unit:%,MODIS dataset) of physical quantities of total cloud in the El Ni?o decaying summers(a,c,e.EPEl Ni?o,b,d,f.CPEl Ni?o;a,b.pressure in cloud top,c,d.temperature in cloud top,e,f.optical thickness)

    3.1.4 MICAPS數(shù)據(jù):厄爾尼諾事件對云量和低云狀出現(xiàn)頻率的影響

    由于數(shù)據(jù)長度只有14 a,東部型和中部型事件個(gè)例分別為1個(gè)和3個(gè),所得結(jié)果可能不具有統(tǒng)計(jì)上的顯著意義和廣泛代表性,但其包含的低云狀等量可用于對常規(guī)站點(diǎn)數(shù)據(jù)進(jìn)行一些補(bǔ)充分析。

    從MICAPS數(shù)據(jù)中厄爾尼諾事件次年夏季白天降水和云量的異常分布(圖7)可見,東部型厄爾尼諾對應(yīng)的降水(圖7a)在全國范圍內(nèi)基本都是大幅度增加的,只黑龍江和新疆部分區(qū)域減少;總云量(圖7c)在東北、內(nèi)蒙古東部、新疆、西藏、四川西部、華南西部和江南東部是減少的,在其余地區(qū)是增多的;低云量的分布形勢(圖7e)同總云量,但云量減少區(qū)域的范圍更大。東部型厄爾尼諾總云量和低云量的異常分布與降水相似,但異常增多區(qū)域略小一些,云量的增減幅度都在20%以內(nèi)。中部型厄爾尼諾降水(圖7b)在東北、華北南部、黃淮西部、江淮、西南北部、青海—西藏大部分區(qū)域是增多的,其余地區(qū)是減少的;總云量(圖7d)在內(nèi)蒙古中西部、華北、江淮和江南部分地區(qū)是減少的,其余地區(qū)是增多的;低云量的分布形勢(圖7f)大致同總云量,但云量減少區(qū)較小。與降水相比,中部型厄爾尼諾的云量減少區(qū)域面積小且幅度基本不超過10%,增多區(qū)域面積大,部分地區(qū)幅度超過25%。東部型厄爾尼諾對應(yīng)的降水和云量異常與站點(diǎn)數(shù)據(jù)差異較大,可能是因?yàn)镸ICAPS數(shù)據(jù)的東部型事件只有一個(gè)個(gè)例,不具有代表性;中部型厄爾尼諾對應(yīng)的降水和云量異常與站點(diǎn)數(shù)據(jù)較為一致。夜間降水和云量的異常分布形勢基本同白天(圖略)。

    圖7 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季白天降水(a、b)和云量(c、d.總云量,e、f.低云量)距平百分率(單位:%,MICAPS數(shù)據(jù))Fig.7 Percentage anomalies(unit:%,MICAPSdataset)of precipitation and cloud amount in the El Ni?o decaying summers at daytime(a,c,e.EPEl Ni?o,b,d,f. CPEl Ni?o;a,b.precipitation,c,d.total cloud amount,e,f.low cloud amount)

    MICAPS數(shù)據(jù)將低云狀分為9種,表2給出其編碼和伴隨的天氣現(xiàn)象(中國氣象局,2003)。參考高翠翠等(2015)關(guān)于云出現(xiàn)頻率的統(tǒng)計(jì)量:各站點(diǎn)各類云的發(fā)生頻率等于各站點(diǎn)各類云出現(xiàn)的次數(shù)與各站點(diǎn)的觀測次數(shù)的比值,以百分?jǐn)?shù)(%)表示。由于低云狀種類較多,細(xì)分之后會出現(xiàn)氣候平均值為0的情況,無法計(jì)算距平百分率,故而低云狀的出現(xiàn)頻率異常用距平數(shù)據(jù)進(jìn)行分析。接下來對厄爾尼諾事件次年中國夏季9種低云的出現(xiàn)頻率進(jìn)行差異分析。

    由表2可知,淡積云常出現(xiàn)在晴天,其他8種低云都有可能產(chǎn)生降水。圖8所示為厄爾尼諾事件次年夏季平均白天部分低云(層積云(圖8a、b),碎雨云(圖8c、d),鬃狀/砧狀積雨云(圖8e、f))出現(xiàn)頻率相對于氣候平均態(tài)的距平。兩類厄爾尼諾事件次年夏季,中國地區(qū)白天低云量異常增多區(qū)域上層積云、碎雨云和積雨云等3類低云的出現(xiàn)頻率是增加的,頻率增加幅度基本在3.5%以內(nèi)。東部型厄爾尼諾事件中,華北西部和西北東部地區(qū)主要是層積云和鬃狀/砧狀積雨云的出現(xiàn)頻率增加;陜西、江漢和黃淮地區(qū)主要是層積云和碎雨云的出現(xiàn)頻率增加;西南東部和北部、華南東部、江南和江淮地區(qū)則是鬃狀/砧狀積雨云的出現(xiàn)頻率增加。中部型厄爾尼諾事件中,東北地區(qū)主要是層積云和鬃狀/砧狀積雨云的出現(xiàn)頻率增加;黃河上游流域?qū)臃e云的出現(xiàn)頻率增加;西南北部和長江—黃河之間地區(qū)是層積云和碎雨云的出現(xiàn)頻率增加。夜間低云狀出現(xiàn)頻率的距平情況與白天基本一致(圖略)。結(jié)合圖7a、b可知,兩類厄爾尼諾事件異常雨區(qū)的形成與層積云、雨層云和積雨云等低云狀的出現(xiàn)頻率增加密切相關(guān)。

    圖8 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季白天低云(a、b.普通層積云,c、d.碎雨云,e、f.鬃狀/砧狀積雨云)出現(xiàn)頻率的距平(單位:%,MICAPS數(shù)據(jù))Fig.8 Anomalies(unit:%,MICAPS dataset)of frequency of low cloud forms in the El Ni?o decaying summers at daytime(a,c,e.EPEl Ni?o,b,d,f.CPEl Ni?o;a,b.stratocumulus,c,d.fractonimbus,e,f.cumulonimbus capillatus and incus)

    表2 MICAPS數(shù)據(jù)低云狀的類別及其伴隨天氣現(xiàn)象Table2 Low cloud forms and accompanied weather phenomena in the MICAPSdataset

    3.2 厄爾尼諾對降水和云關(guān)系的影響

    曾昭美等(1993)在其研究中提到,中國年平均總云量與年降水量有良好的線性正相關(guān)關(guān)系(相關(guān)系數(shù)0.624,達(dá)到0.001顯著性水平),本研究對厄爾尼諾次年夏季降水與云多種相關(guān)物理量的異常進(jìn)行了線性相關(guān)分析。

    3.2.1 站點(diǎn)數(shù)據(jù)的降水與云量距平百分率的相關(guān)分析

    從厄爾尼諾事件次年夏季降水與云量異常的相關(guān)分析(圖9)可見,東部型厄爾尼諾次年夏季,除內(nèi)蒙古和新疆的小部分區(qū)域外,降水距平百分率與白天平均總云量(圖9a)、白天平均低云量(圖9c)和夜間平均低云量(圖9e)距平百分率基本呈現(xiàn)正相關(guān),中國東部地區(qū)較多站點(diǎn)都通過了顯著性水平為0.05的相關(guān)顯著性檢驗(yàn)。中部型厄爾尼諾對應(yīng)的情況(圖9b、d、f)更復(fù)雜一些,在內(nèi)蒙古、華北和西南等地區(qū)的小部分區(qū)域存在負(fù)相關(guān),但總體也呈正相關(guān),正相關(guān)區(qū)域中也有不少站點(diǎn)通過了相關(guān)顯著性檢驗(yàn)。

    圖9 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季降水與云量距平百分率的相關(guān)系數(shù)(站點(diǎn)數(shù)據(jù);a、b.白天降水與總云量,c、d.白天降水與低云量,e、f.夜間降水與低云量)Fig.9 Correlation between percentage anomalies of precipitation and cloud amount in the El Ni?o decaying summers(station dataset;a,c,e.EPEl Ni?o,b,d,f.CPEl Ni?o;a,b. precipitation and total cloud amount at daytime,c,d. precipitation and low cloud amount at daytime,e,f.precipitation and low cloud amount at nighttime)

    3.2.2 站點(diǎn)數(shù)據(jù)的降水與云量距平百分率的線性回歸分析

    將通過降水與云量正相關(guān)顯著性檢驗(yàn)的站點(diǎn)挑選出來,對這些站點(diǎn)的云量與降水距平百分率進(jìn)行線性回歸分析,結(jié)果如圖10所示,東部型厄爾尼諾對應(yīng)的云量與降水的線性關(guān)系非常明顯,且都通過了較高水平的顯著性檢驗(yàn);白天總云量與降水異常的增加比例(圖10a)為1:2.82,低云量(圖10c為白天,圖10e為夜間)與降水異常的增加比例為1:0.82—0.91。中部型厄爾尼諾對應(yīng)的情況則較為復(fù)雜,云量與降水的線性關(guān)系不如東部型厄爾尼諾明顯,但也都通過了顯著性檢驗(yàn);總云量與降水異常的增加比例(圖10b)為1:3.00,低云量(圖10d為白天,圖10f為夜間)與降水異常的增加比例為1:1.11—1.17。云量異常與降水異常有著很好的相關(guān)關(guān)系,低云量異常與降水異常的線性回歸關(guān)系更為明顯(顯著性水平高2—3個(gè)量級):總云量異常增加1%,降水異常增加約為3%;低云量異常增加1%,降水異常隨之增加約為1%。

    圖10 厄爾尼諾(a、c、e.東部型,b、d、f.中部型)事件次年夏季降水與云量距平百分率的線性回歸(站點(diǎn)數(shù)據(jù),公式為回歸方程,Num表示擬合站點(diǎn)數(shù),p為回歸方程F顯著性檢驗(yàn)的p值;a、b.白天降水與總云量,c、d. 白天降水與低云量,e、f.夜間降水與低云量)Fig.10 Linear regression between percentage anomalies of precipitation and cloud amount in the El Ni?o decaying summers(station dataset,the formula on the figure is the regression equation,"Num"represents the number of fitting stations,and "p"is the p-value of F test of regression equation;a,c,e.EP El Ni?o,b,d,f.CP El Ni?o;a,b.precipitation and total cloud amount at daytime,c,d.precipitation and low cloud amount at daytime,e,f.precipitation and low cloud amount at nighttime)

    3.2.3 站點(diǎn)數(shù)據(jù)的降水與ISCCP云相關(guān)物理量距平百分率的相關(guān)分析

    圖11是厄爾尼諾事件次年夏季降水與ISCCP云量距平百分率的相關(guān)系數(shù)分布。東部型厄爾尼諾對應(yīng)的總云量(圖11a)和高云量(圖11c)與降水距平百分率基本呈正相關(guān)關(guān)系,高云中的卷層云(圖11e)和深對流云云量(圖11g)與降水異常有很好的正相關(guān)關(guān)系。中部型厄爾尼諾對應(yīng)的總云量(圖11b)與降水異常的相關(guān)較差,高云量(圖11d)以及高云中的卷層云(圖11f)和深對流云(圖11h)云量與降水距平百分率有很好的正相關(guān)。由圖11可知,高云量和深對流云云量與降水異常的正相關(guān)比較強(qiáng)。

    圖11 厄爾尼諾(a、c、e、g.東部型,b、d、f、h.中部型)事件次年夏季站點(diǎn)降水與ISCCP云量距平百分率的相關(guān)系數(shù)(a、b.總降水與總云量,c、d.總降水與高云量,e、f.總降水與卷層云云量,g、h.總降水與深對流云云量;通過0.05水平顯著性檢驗(yàn)的區(qū)域打點(diǎn)標(biāo)注)Fig.11 Correlation between percentage anomalies of precipitation of station dataset and cloud amount in the ISCCPdataset in the El Ni?o decaying summers(a,c,e,g.EPEl Ni?o,b,d,f,h.CPEl Ni?o;a,b.total precipitation and total cloud amount,c,d.total precipitation and high cloud amount,e,f.total precipitation and cirrostratus cloud amount,g,h.total precipitation and deep convective cloud amount;dots denote passing the test of 0.05 significance level)

    圖12是厄爾尼諾次年夏季ISCCP數(shù)據(jù)總云的其他相關(guān)物理量與站點(diǎn)降水距平百分率的相關(guān)系數(shù)分布。在大部分地區(qū),兩類事件對應(yīng)的總云光學(xué)厚度和云水路徑與降水異常呈正相關(guān),云頂氣壓、云頂溫度與降水異常的相關(guān)關(guān)系(圖略)較為復(fù)雜,正、負(fù)相關(guān)同時(shí)存在。對于卷層云和深對流云(圖略),云頂氣壓、云頂溫度、光學(xué)厚度和云水路徑與降水距平百分率的相關(guān)較差,正、負(fù)相關(guān)同時(shí)存在。

    圖12 厄爾尼諾(a、c.東部型,b、d.中部型)事件次年夏季站點(diǎn)降水與ISCCP總云相關(guān)物理量距平百分率的相關(guān)系數(shù)(a、b.總降水與總云光學(xué)厚度,c、d.總降水與云水路徑)Fig.12 Correlation between percentage anomalies of precipitation from station dataset and total cloud related physical quantities from the ISCCPdataset in the El Ni?o decaying summer (a,c.EPEl Ni?o,b,d.CPEl Ni?o;total precipitation and total cloud (a,b)optical thickness and (c,d)cloud water path)

    3.3 厄爾尼諾對不同類型降水的影響

    對于厄爾尼諾事件次年夏季站點(diǎn)降水率與TRMM地表降水率的分布形勢進(jìn)行對比(圖略)。在空間分布上,TRMM衛(wèi)星與站點(diǎn)降水分布形勢一致,降水率從西北向東南遞增;在降水幅度上,TRMM衛(wèi)星在全域存在高估降水量的現(xiàn)象,與楊東等(2019)關(guān)于洞庭湖流域TRMM降水?dāng)?shù)據(jù)精度的相關(guān)研究結(jié)論一致。利用TRMM衛(wèi)星的降水資料研究厄爾尼諾事件對不同類型降水的影響,給出厄爾尼諾事件次年夏季近地表對流降水和層云降水總和的距平百分率以及總降水異常中二者所占的比例(圖13)。東部型厄爾尼諾事件次年夏季,淮河流域、長江流域和珠江流域總降水都是增多的(圖13a),云南、江南東部和華南南部等地降水減少;降水增加極大值超過20%,大值區(qū)位于江漢北部和西南東部地區(qū)。圖13c、e是東部型厄爾尼諾事件對流降水和層云降水異常對總降水異常的貢獻(xiàn),在整個(gè)南方區(qū)域,對流降水異常的貢獻(xiàn)都是正值,表明對流降水異常與總降水異常分布趨勢一致,且貢獻(xiàn)率基本在80%以上;除江南西部等小部分地區(qū)是負(fù)值外,層云降水異常的貢獻(xiàn)率也基本都是正值,大部分地區(qū)貢獻(xiàn)率在30%以內(nèi)。中部型厄爾尼諾事件對應(yīng)的降水異常增加區(qū)主要在長江以北地區(qū)(圖13b),降水大值區(qū)位于華北南部和西北地區(qū)南部,增加幅度基本不超過20%。圖13d、f是中部型厄爾尼諾事件對流降水和層云降水異常對總降水異常的貢獻(xiàn),整個(gè)區(qū)域內(nèi)對流降水對總降水異常的貢獻(xiàn)也都是正的,貢獻(xiàn)率也基本在80%以上;層云降水異常對總降水的貢獻(xiàn)也幾乎均為正值,表現(xiàn)為與總降水變化一致,在西南東部和華南東部等地區(qū)占比幅度超過30%。對比圖13a、b和相同年份站點(diǎn)數(shù)據(jù)兩類厄爾尼諾事件的降水距平(圖略)可知,TRMM降水異常分布與站點(diǎn)降水大致相同。綜上可知,東部型厄爾尼諾事件次年夏季黃河以南地區(qū)對流降水和層云降水均有增加,中部型厄爾尼諾事件次年夏季長江以北地區(qū)兩類降水同步增多、長江以南地區(qū)兩類降水同步減少,兩類厄爾尼諾事件中對流降水增加對總降水變化的貢獻(xiàn)率在大部分地區(qū)都達(dá)到了80%。孫齊穎等(2017)利用CMORPH 3 h衛(wèi)星降水?dāng)?shù)據(jù)分析了北半球夏季對流性降水與厄爾尼諾事件(1個(gè)東部型和3個(gè)中部型合成)的關(guān)系,也得出了厄爾尼諾事件次年中國陸地(尤其是黃河以南)對流性降水明顯偏多、非對流性降水則表現(xiàn)為長江—黃河之間偏多和江南偏少的結(jié)論。

    圖13 厄爾尼諾 (a、c、e.東部型,b、d、f.中部型)事件次年夏季近地表總降水的距平百分率(a、b)以及對流降水 (c、d) 和層云降水(e、f)異常占總降水異常的百分?jǐn)?shù)(單位:%,TRMM數(shù)據(jù))Fig.13 Percentage anomalies(unit:%,TRMM dataset)of total precipitation (a,b)and the ratios accounted for by convective(c,d)and stratus(e,f)precipitation anomalies in the El Ni?o decaying summers(a,c,e.EPEl Ni?o,b,d,f.CPEl Ni?o)

    4 結(jié)論與討論

    文中主要探討了不同類型厄爾尼諾事件次年夏季中國地區(qū)云(云量、云頂氣壓、云頂溫度、云水路徑、光學(xué)厚度和低云狀出現(xiàn)頻率)的變化特征及其與降水變化的聯(lián)系(相關(guān)關(guān)系和線性回歸關(guān)系)。研究中使用到多種數(shù)據(jù),由于站點(diǎn)數(shù)據(jù)時(shí)間長度較長、缺測率較低且經(jīng)過嚴(yán)格的質(zhì)量控制,以站點(diǎn)資料為參考進(jìn)行比較,各種數(shù)據(jù)之間的偏差主要在西部地區(qū),且數(shù)值偏差較小。對比表明這些數(shù)據(jù)質(zhì)量高,所得結(jié)果具有較高的可信度。

    (1)東部型厄爾尼諾次年夏季中國地區(qū)形成南、北兩條異常雨帶,中部型厄爾尼諾次年夏季中國地區(qū)形成一條異常雨帶,這個(gè)結(jié)果與已有研究結(jié)果基本一致。厄爾尼諾次年夏季的低云量和總云量的異常分布形勢大致與降水異常分布形勢相同,降水異常增加的地區(qū)對應(yīng)的云量也是增多的。云量與降水異常存在很好的正相關(guān)。在通過正相關(guān)關(guān)系顯著性檢驗(yàn)的站點(diǎn)中,總云量和低云量與降水異常分別存在1:3和1:1的線性關(guān)系,且都通過了0.05水平的顯著性檢驗(yàn)。也就是說,厄爾尼諾次年夏季,中國東部地區(qū)總云量異常增多1%,降水異常相應(yīng)增加約3%;低云量異常增多1%,降水異常相應(yīng)增加約1%。而在異常雨帶所在地區(qū)積雨云、雨層云和層積云的出現(xiàn)頻率有所增加。

    (2)兩類厄爾尼諾事件對應(yīng)的高云量和深對流云云量異常的分布與降水異常形勢較為一致,總云云頂氣壓和云頂溫度與降水的距平百分率分布形勢大致相反,光學(xué)厚度和云水路徑與降水異常分布形勢相似且正大值區(qū)基本重合。厄爾尼諾事件的發(fā)生使得次年夏季中國季風(fēng)區(qū)對流活動(dòng)增強(qiáng),云水含量和對流云云量增加,云層增厚,云頂向上發(fā)展,有利于降水的發(fā)生。

    (3)黃河以南地區(qū)兩類厄爾尼諾事件對應(yīng)雨帶的形成是對流降水和層云降水同步異常增加的結(jié)果,但以對流降水增加為主,在大部分地區(qū)其對總降水變化的貢獻(xiàn)率都達(dá)到了80%。中部型事件中長江以南地區(qū)干旱少雨也主要是對流降水減少造成的,不過部分地區(qū)層云降水減少的貢獻(xiàn)率超過30%。

    在厄爾尼諾事件次年夏季,站點(diǎn)數(shù)據(jù)的云量與降水異常的線性回歸分析得到總云量和低云量與降水增加的比例分別接近1:3和1:1的結(jié)論,說明低云對降水增加具有較大的影響;MICAPS數(shù)據(jù)對站點(diǎn)數(shù)據(jù)補(bǔ)充分析表明,低云中積雨云、雨層云和層積云的出現(xiàn)頻率增加,對降水影響大;ISCCP觀測數(shù)據(jù)的相關(guān)分析則表明高云和深對流云云量與降水異常存在顯著的正相關(guān);TRMM衛(wèi)星數(shù)據(jù)相關(guān)分析則表明,黃河以南地區(qū)對流降水異常對總降水變化的貢獻(xiàn)基本高達(dá)80%。ISCCP衛(wèi)星數(shù)據(jù)根據(jù)云頂高度分高中低云,站點(diǎn)數(shù)據(jù)則根據(jù)云底高度分高中低云,通過對比可知ISCCP數(shù)據(jù)中高云的深對流云和站點(diǎn)數(shù)據(jù)低云中深厚的積雨云與TRMM數(shù)據(jù)中對流降水的變化是一致的。由此可知,厄爾尼諾事件次年夏季,中國東部地區(qū)低云對降水具有良好的指示作用;總降水的異常主要來自對流降水的變化;對流降水異常增加與深對流云的異常增加有關(guān)。

    東部型厄爾尼諾事件次年夏季,東亞-太平洋遙相關(guān)型建立,副高西側(cè)的強(qiáng)西南氣流將來自太平洋的大量水汽持續(xù)輸送至中國中東部地區(qū);同時(shí),中高緯度異常雙阻塞型環(huán)流建立,將北冰洋的水汽帶到中國北部地區(qū)(吳萍等,2017)。兩條水汽通道導(dǎo)致中國除江淮地區(qū)外,各地云量的異常增多,由此形成南北兩條異常雨帶。中部型厄爾尼諾事件次年夏季,菲律賓異常反氣旋強(qiáng),影響范圍廣,夏季風(fēng)將來自西太平洋的水汽集中輸送到江淮地區(qū),中國北方地區(qū)水汽不足而長江以南地區(qū)被副熱帶高壓控制,故而在長江—黃河流域云量異常增多,形成了一條異常雨帶(吳萍等,2017)。賈子康等(2020)的研究也指出,夏季華南型和江淮型高溫分別與東部型和中部型厄爾尼諾衰減密切相關(guān),中部型事件中副熱帶高壓偏強(qiáng)西伸北抬,控制范圍更偏北,故而雨帶北移。

    東部型事件對應(yīng)的降水異常與云相關(guān)物理量間的關(guān)系呈現(xiàn)出一定的規(guī)律性,結(jié)果較為理想;中部型事件對應(yīng)的降水異常與云相關(guān)物理量間的關(guān)系則無規(guī)律可循,也表明雖然中部型都是中等和弱厄爾尼諾事件,但其對氣候的影響更大、更復(fù)雜和具有更大的不確定性。全球變暖背景下東太平洋厄爾尼諾變率增大(Cai,et al,2018);20世紀(jì)90年代以來,中部型厄爾尼諾事件發(fā)生頻率顯著增加,可能與全球變暖,尤其是赤道太平洋海溫的年代際升高有關(guān)(Yeh,et al,2009)。Tang等(2021)認(rèn)為模式中的系統(tǒng)性偏差導(dǎo)致全球變暖背景下的極端厄爾尼諾頻率變化被顯著高估,基于訂正后的氣候態(tài)變化,極端厄爾尼諾頻率將幾乎不變。全球變暖背景下,同樣振幅的ENSO能導(dǎo)致更大的對流層水汽異常,進(jìn)而造成更大的全球大氣環(huán)流、氣溫和降水異常,揭示了全球變暖下ENSO的氣候影響會顯著增強(qiáng),未來ENSO可能會造成更大的氣候?yàn)?zāi)害(Hu,et al,2021)。厄爾尼諾事件是影響東亞夏季風(fēng)最顯著的自然因子,溫室效應(yīng)至今仍不斷增強(qiáng),全球變暖是21世紀(jì)最重要的環(huán)境問題;全球變暖背景下,厄爾尼諾事件發(fā)生更加頻繁、復(fù)雜,不確定性更大,對氣候產(chǎn)生的影響也具有更大的不確定性。

    目前全球變暖仍在持續(xù),極端天氣氣候事件發(fā)生的頻率增加、強(qiáng)度增強(qiáng)。近30年,厄爾尼諾事件的強(qiáng)度有增大的趨勢,且其復(fù)雜性、多樣性和不確定性也在增大。氣象觀測歷史記錄中,厄爾尼諾事件的樣本為數(shù)不多,本研究將事件分為東部型和中部型進(jìn)行分析,樣本數(shù)量有限,故而所得結(jié)果在統(tǒng)計(jì)上較難通過顯著性檢驗(yàn)。計(jì)算能力快速增強(qiáng)和科學(xué)技術(shù)的不斷進(jìn)步推動(dòng)了氣候模式的發(fā)展和成熟,如今氣候模式被廣泛應(yīng)用于模擬當(dāng)前氣候和預(yù)測未來氣候的變化。今后可利用模式進(jìn)行數(shù)值模擬試驗(yàn),獲得更多厄爾尼諾事件的樣本,以期在厄爾尼諾事件影響云和降水方面進(jìn)行更細(xì)致深入的研究工作,為更好地應(yīng)對氣候變化和防災(zāi)減災(zāi)工作提供科學(xué)決策依據(jù)。

    日韩av在线大香蕉| 亚洲婷婷狠狠爱综合网| av在线观看视频网站免费| 国产白丝娇喘喷水9色精品| 看黄色毛片网站| 久久久国产成人免费| 久久久久久久午夜电影| 久久亚洲精品不卡| 欧美日本亚洲视频在线播放| 91av网一区二区| 久久久久久九九精品二区国产| av黄色大香蕉| 国模一区二区三区四区视频| av免费在线看不卡| 精品久久久久久电影网 | 一级毛片aaaaaa免费看小| 99久久精品国产国产毛片| 国产国拍精品亚洲av在线观看| 边亲边吃奶的免费视频| 国产精品人妻久久久久久| 精品无人区乱码1区二区| 久久综合国产亚洲精品| 久久久久精品久久久久真实原创| 一级黄色大片毛片| 大又大粗又爽又黄少妇毛片口| 免费看av在线观看网站| 99久久无色码亚洲精品果冻| 日韩强制内射视频| 国产欧美另类精品又又久久亚洲欧美| 九九久久精品国产亚洲av麻豆| 不卡视频在线观看欧美| 亚洲国产精品成人久久小说| 日韩成人av中文字幕在线观看| 18禁在线播放成人免费| 国产精品爽爽va在线观看网站| 床上黄色一级片| 久久99蜜桃精品久久| 赤兔流量卡办理| 国内精品宾馆在线| 成年女人永久免费观看视频| 久久99精品国语久久久| 哪个播放器可以免费观看大片| 深夜a级毛片| 亚洲美女视频黄频| 亚洲欧洲日产国产| 国产亚洲最大av| 欧美bdsm另类| 日本五十路高清| 只有这里有精品99| 国产女主播在线喷水免费视频网站 | 全区人妻精品视频| 久久久久久久久久久免费av| 99在线视频只有这里精品首页| 午夜福利成人在线免费观看| 少妇人妻一区二区三区视频| 亚洲欧美成人综合另类久久久 | 久热久热在线精品观看| 成年av动漫网址| 亚洲av熟女| 久久99蜜桃精品久久| 男的添女的下面高潮视频| 成人亚洲欧美一区二区av| 中国美白少妇内射xxxbb| 午夜精品一区二区三区免费看| 别揉我奶头 嗯啊视频| 性插视频无遮挡在线免费观看| 色综合站精品国产| 国产私拍福利视频在线观看| 波多野结衣高清无吗| 69人妻影院| 久久精品久久久久久噜噜老黄 | 一夜夜www| 天堂√8在线中文| 搡女人真爽免费视频火全软件| 在线观看av片永久免费下载| 亚洲av二区三区四区| 久久久久久久久中文| 乱码一卡2卡4卡精品| 一级爰片在线观看| 一区二区三区乱码不卡18| 激情 狠狠 欧美| 亚洲熟妇中文字幕五十中出| 久久久久精品久久久久真实原创| av播播在线观看一区| 精品国产露脸久久av麻豆 | 简卡轻食公司| 免费黄色在线免费观看| 桃色一区二区三区在线观看| 亚洲美女视频黄频| 亚洲欧美日韩无卡精品| 免费看a级黄色片| 亚洲精品成人久久久久久| 免费观看a级毛片全部| 成人国产麻豆网| 婷婷六月久久综合丁香| 国产极品天堂在线| 国产精品福利在线免费观看| 五月伊人婷婷丁香| 少妇人妻一区二区三区视频| 日韩欧美三级三区| 日本av手机在线免费观看| av免费在线看不卡| 高清午夜精品一区二区三区| 亚洲人成网站在线观看播放| 大又大粗又爽又黄少妇毛片口| 国产成人午夜福利电影在线观看| 狂野欧美白嫩少妇大欣赏| 日韩av在线大香蕉| 97在线视频观看| 国产伦精品一区二区三区视频9| 身体一侧抽搐| 亚洲精品亚洲一区二区| 国产免费福利视频在线观看| 最近中文字幕高清免费大全6| 日韩精品青青久久久久久| 天堂网av新在线| 少妇人妻一区二区三区视频| 国产精品综合久久久久久久免费| 中文字幕av成人在线电影| 99视频精品全部免费 在线| 久久热精品热| 大话2 男鬼变身卡| 日韩强制内射视频| 搡女人真爽免费视频火全软件| 亚洲国产高清在线一区二区三| 大香蕉久久网| 成年女人看的毛片在线观看| 大香蕉久久网| 亚洲中文字幕一区二区三区有码在线看| 国产成人freesex在线| 天堂影院成人在线观看| 中文字幕亚洲精品专区| 亚洲av电影不卡..在线观看| 简卡轻食公司| videossex国产| 国产黄片视频在线免费观看| 欧美激情在线99| 日韩欧美精品免费久久| 69av精品久久久久久| 22中文网久久字幕| 大又大粗又爽又黄少妇毛片口| 舔av片在线| 亚洲综合色惰| 九草在线视频观看| 99久久九九国产精品国产免费| 久久久国产成人精品二区| 一级毛片电影观看 | 天堂网av新在线| 婷婷色av中文字幕| 一区二区三区乱码不卡18| 亚洲精品,欧美精品| 别揉我奶头 嗯啊视频| 在线播放国产精品三级| 国产 一区 欧美 日韩| 直男gayav资源| 国语对白做爰xxxⅹ性视频网站| 亚洲国产欧美人成| 一区二区三区高清视频在线| 91精品一卡2卡3卡4卡| 在线免费观看不下载黄p国产| 久久精品91蜜桃| 国产成人一区二区在线| 久久精品国产亚洲av天美| 天天躁日日操中文字幕| 丰满乱子伦码专区| 国产精品av视频在线免费观看| 欧美最新免费一区二区三区| 啦啦啦韩国在线观看视频| 嫩草影院新地址| 午夜亚洲福利在线播放| 国产精品久久久久久av不卡| 国产人妻一区二区三区在| 亚洲av一区综合| 激情 狠狠 欧美| 哪个播放器可以免费观看大片| 在线a可以看的网站| 尾随美女入室| 99久久九九国产精品国产免费| 人人妻人人澡人人爽人人夜夜 | 欧美zozozo另类| 午夜老司机福利剧场| 中文字幕精品亚洲无线码一区| 久久精品久久久久久久性| 男人狂女人下面高潮的视频| 极品教师在线视频| 国产成人a∨麻豆精品| av又黄又爽大尺度在线免费看 | 免费观看a级毛片全部| 久久久亚洲精品成人影院| 日韩成人伦理影院| 亚洲在线自拍视频| 亚洲国产成人一精品久久久| 免费黄色在线免费观看| 丝袜美腿在线中文| 少妇的逼好多水| 国产成人精品久久久久久| 久久精品国产亚洲av涩爱| 亚洲人成网站在线播| 亚洲国产色片| 69av精品久久久久久| 99国产精品一区二区蜜桃av| 欧美一区二区国产精品久久精品| 在线观看一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲第一区二区三区不卡| 不卡视频在线观看欧美| 级片在线观看| 亚洲成人久久爱视频| 国产老妇伦熟女老妇高清| 水蜜桃什么品种好| 综合色av麻豆| 日韩中字成人| 精品熟女少妇av免费看| 亚洲精品亚洲一区二区| 午夜a级毛片| 美女国产视频在线观看| 高清av免费在线| 美女黄网站色视频| 一级爰片在线观看| 亚洲av电影不卡..在线观看| videossex国产| 亚洲婷婷狠狠爱综合网| 99久国产av精品国产电影| 中国国产av一级| av又黄又爽大尺度在线免费看 | 日韩视频在线欧美| 亚洲天堂国产精品一区在线| 中文字幕久久专区| 午夜久久久久精精品| 三级毛片av免费| 精品欧美国产一区二区三| 亚洲四区av| 一级二级三级毛片免费看| av福利片在线观看| 搡老妇女老女人老熟妇| 日韩欧美三级三区| 久久久久久大精品| 中文字幕精品亚洲无线码一区| 在线免费观看的www视频| 亚洲无线观看免费| 欧美丝袜亚洲另类| 国产精品久久视频播放| 99久久精品热视频| 亚洲真实伦在线观看| 永久网站在线| 我的老师免费观看完整版| 国产成人精品久久久久久| 高清视频免费观看一区二区 | 国产极品天堂在线| 日韩av不卡免费在线播放| 美女脱内裤让男人舔精品视频| 中文亚洲av片在线观看爽| 亚洲性久久影院| 国产一区有黄有色的免费视频 | 中文字幕精品亚洲无线码一区| 日韩亚洲欧美综合| 床上黄色一级片| 亚洲国产欧洲综合997久久,| 精品午夜福利在线看| 村上凉子中文字幕在线| 国内精品美女久久久久久| 国产亚洲精品av在线| 欧美日韩国产亚洲二区| 午夜免费男女啪啪视频观看| 高清视频免费观看一区二区 | 亚洲欧美精品专区久久| 精品一区二区三区人妻视频| 人妻夜夜爽99麻豆av| 看非洲黑人一级黄片| 男女国产视频网站| 国产美女午夜福利| 美女内射精品一级片tv| 亚洲国产色片| 久久久a久久爽久久v久久| 国产午夜福利久久久久久| 在线免费观看的www视频| 免费观看的影片在线观看| 91精品国产九色| 成人二区视频| 久久人妻av系列| 国产日韩欧美在线精品| 亚洲国产成人一精品久久久| 色视频www国产| 日韩 亚洲 欧美在线| 乱系列少妇在线播放| 国产精品一区www在线观看| 人妻少妇偷人精品九色| 最近的中文字幕免费完整| 日韩欧美精品免费久久| 国内精品宾馆在线| 非洲黑人性xxxx精品又粗又长| 简卡轻食公司| 在线免费观看不下载黄p国产| 国产老妇伦熟女老妇高清| 国产视频首页在线观看| 乱人视频在线观看| 欧美一区二区亚洲| 能在线免费观看的黄片| 国产真实乱freesex| 精品无人区乱码1区二区| 成人二区视频| 老师上课跳d突然被开到最大视频| 一级毛片我不卡| 亚洲国产最新在线播放| 久久精品久久久久久久性| 久久99热6这里只有精品| 日韩制服骚丝袜av| 纵有疾风起免费观看全集完整版 | 国产极品天堂在线| 亚洲av成人精品一二三区| 偷拍熟女少妇极品色| 久久午夜福利片| 午夜a级毛片| www.色视频.com| 欧美日韩精品成人综合77777| 久久久久久久久大av| 久久久久久久午夜电影| 观看美女的网站| 51国产日韩欧美| 观看免费一级毛片| 欧美激情国产日韩精品一区| 日本免费在线观看一区| 欧美日韩精品成人综合77777| 长腿黑丝高跟| 国产精华一区二区三区| 色尼玛亚洲综合影院| a级一级毛片免费在线观看| 色综合色国产| 三级国产精品欧美在线观看| 一边亲一边摸免费视频| 日本猛色少妇xxxxx猛交久久| 精品一区二区三区视频在线| 乱人视频在线观看| 亚洲精品日韩av片在线观看| 久久久久久九九精品二区国产| 精品一区二区三区人妻视频| 永久免费av网站大全| 黄片无遮挡物在线观看| 国产av码专区亚洲av| 插逼视频在线观看| 男人舔奶头视频| 日韩欧美精品v在线| 亚洲伊人久久精品综合 | 国产精品麻豆人妻色哟哟久久 | 精品一区二区三区视频在线| 97热精品久久久久久| 美女xxoo啪啪120秒动态图| 亚洲av电影在线观看一区二区三区 | 国产一级毛片在线| 精品99又大又爽又粗少妇毛片| 国产单亲对白刺激| 久久久久久久国产电影| 秋霞伦理黄片| 日本免费一区二区三区高清不卡| 婷婷六月久久综合丁香| 色播亚洲综合网| 女的被弄到高潮叫床怎么办| 国产在线男女| 国产伦在线观看视频一区| 天天一区二区日本电影三级| kizo精华| 全区人妻精品视频| 成人毛片60女人毛片免费| 少妇人妻一区二区三区视频| 国产精品福利在线免费观看| 3wmmmm亚洲av在线观看| 淫秽高清视频在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品无大码| 欧美日韩精品成人综合77777| 在线观看一区二区三区| 精品久久久久久电影网 | 国产 一区精品| 五月玫瑰六月丁香| 亚洲最大成人av| 午夜福利在线观看吧| 久久精品国产亚洲网站| 国产在视频线在精品| a级一级毛片免费在线观看| 乱码一卡2卡4卡精品| 两性午夜刺激爽爽歪歪视频在线观看| 色吧在线观看| 乱人视频在线观看| 人妻少妇偷人精品九色| 国产精品一区www在线观看| 国产91av在线免费观看| 午夜久久久久精精品| 久久久久久久久久成人| 成人三级黄色视频| 狂野欧美激情性xxxx在线观看| 国产精华一区二区三区| 国产 一区 欧美 日韩| 久久久久久久久久久免费av| 亚洲乱码一区二区免费版| 久久精品夜夜夜夜夜久久蜜豆| 99热网站在线观看| 亚洲aⅴ乱码一区二区在线播放| 麻豆乱淫一区二区| 久久久欧美国产精品| 免费观看在线日韩| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 免费看av在线观看网站| 国产探花极品一区二区| 国产高潮美女av| 99热网站在线观看| 午夜免费激情av| 日韩视频在线欧美| 黄色日韩在线| 国产乱人视频| 亚洲伊人久久精品综合 | 国产午夜精品一二区理论片| 国语自产精品视频在线第100页| 一区二区三区高清视频在线| 亚洲精品亚洲一区二区| 久久久久久久久大av| 最后的刺客免费高清国语| 国产久久久一区二区三区| 一个人看视频在线观看www免费| 亚洲怡红院男人天堂| 男女视频在线观看网站免费| 欧美zozozo另类| 久久精品国产99精品国产亚洲性色| 午夜免费男女啪啪视频观看| av在线蜜桃| 久久婷婷人人爽人人干人人爱| 综合色丁香网| 亚洲四区av| 蜜桃久久精品国产亚洲av| 国产探花在线观看一区二区| 观看美女的网站| 免费搜索国产男女视频| 91精品国产九色| 久久久久久久久久成人| 国内精品美女久久久久久| av免费在线看不卡| 免费不卡的大黄色大毛片视频在线观看 | 2021天堂中文幕一二区在线观| 国产亚洲午夜精品一区二区久久 | 2022亚洲国产成人精品| 国产亚洲午夜精品一区二区久久 | 最近视频中文字幕2019在线8| 啦啦啦观看免费观看视频高清| 成人二区视频| 久久久精品大字幕| 天堂中文最新版在线下载 | 亚洲在久久综合| 看免费成人av毛片| 日本与韩国留学比较| 如何舔出高潮| 日韩一区二区视频免费看| 九草在线视频观看| 永久免费av网站大全| 国产黄a三级三级三级人| 日韩高清综合在线| 成人鲁丝片一二三区免费| 亚洲av中文字字幕乱码综合| 国产av码专区亚洲av| 哪个播放器可以免费观看大片| 国产精品三级大全| av国产久精品久网站免费入址| 国产v大片淫在线免费观看| 18+在线观看网站| 日日啪夜夜撸| 综合色av麻豆| 国产精品一二三区在线看| 18+在线观看网站| 亚洲精品日韩在线中文字幕| 久久精品综合一区二区三区| 久99久视频精品免费| 99久久精品国产国产毛片| 亚洲综合精品二区| 久久精品熟女亚洲av麻豆精品 | 免费观看a级毛片全部| 少妇猛男粗大的猛烈进出视频 | 国产黄色小视频在线观看| 午夜亚洲福利在线播放| 午夜福利在线观看吧| 亚洲人成网站在线播| h日本视频在线播放| 亚洲av.av天堂| videossex国产| 老司机福利观看| 久久精品国产亚洲av涩爱| 久久久久久伊人网av| 国产高清有码在线观看视频| av线在线观看网站| 内射极品少妇av片p| 久久精品国产自在天天线| 久久午夜福利片| 国产在视频线在精品| 日韩av不卡免费在线播放| 最后的刺客免费高清国语| 人妻制服诱惑在线中文字幕| 我要看日韩黄色一级片| 午夜免费男女啪啪视频观看| 国产国拍精品亚洲av在线观看| 国产av一区在线观看免费| 亚洲成人中文字幕在线播放| 人妻少妇偷人精品九色| 九九在线视频观看精品| 天美传媒精品一区二区| 久久精品熟女亚洲av麻豆精品 | 伦精品一区二区三区| 床上黄色一级片| 国产成人精品一,二区| 国产成人91sexporn| 色尼玛亚洲综合影院| 国产精品无大码| 丝袜美腿在线中文| 亚洲美女搞黄在线观看| 久久精品91蜜桃| 国产一区二区在线观看日韩| 亚洲欧美日韩卡通动漫| 亚洲图色成人| 我的老师免费观看完整版| 一级毛片电影观看 | 久久综合国产亚洲精品| av线在线观看网站| 精品福利永久在线观看| 久久久精品94久久精品| 欧美亚洲日本最大视频资源| 国产精品国产三级专区第一集| av国产久精品久网站免费入址| 免费在线观看完整版高清| 亚洲国产欧美在线一区| 亚洲美女搞黄在线观看| 久久99精品国语久久久| 欧美人与性动交α欧美软件 | 一区二区av电影网| 麻豆精品久久久久久蜜桃| 亚洲欧美精品自产自拍| 少妇猛男粗大的猛烈进出视频| 宅男免费午夜| 日韩不卡一区二区三区视频在线| 99热6这里只有精品| 日韩精品免费视频一区二区三区 | 菩萨蛮人人尽说江南好唐韦庄| 秋霞在线观看毛片| 久久99一区二区三区| 9191精品国产免费久久| 如何舔出高潮| 在线观看免费日韩欧美大片| 另类精品久久| 18禁国产床啪视频网站| 婷婷色综合大香蕉| 久久精品国产综合久久久 | 国产亚洲精品第一综合不卡 | 欧美bdsm另类| 91精品伊人久久大香线蕉| 日本欧美视频一区| 美女大奶头黄色视频| 国产亚洲最大av| 国产乱人偷精品视频| 国产精品国产三级专区第一集| 精品卡一卡二卡四卡免费| av视频免费观看在线观看| 国产高清不卡午夜福利| 在线观看三级黄色| a级毛色黄片| 日韩,欧美,国产一区二区三区| 制服诱惑二区| 咕卡用的链子| a级片在线免费高清观看视频| 最黄视频免费看| 熟女人妻精品中文字幕| 岛国毛片在线播放| 一级毛片电影观看| 精品亚洲乱码少妇综合久久| 精品人妻偷拍中文字幕| 一区二区三区四区激情视频| 国产精品国产av在线观看| 制服人妻中文乱码| 亚洲国产毛片av蜜桃av| 中文字幕人妻熟女乱码| 97精品久久久久久久久久精品| 最近中文字幕2019免费版| 中文乱码字字幕精品一区二区三区| 免费看光身美女| 国产精品成人在线| 丰满少妇做爰视频| 国产精品一区二区在线观看99| 黑人高潮一二区| 欧美丝袜亚洲另类| 制服丝袜香蕉在线| 日日爽夜夜爽网站| kizo精华| 成人黄色视频免费在线看| 看免费av毛片| 新久久久久国产一级毛片| 日本黄大片高清| 久久精品久久精品一区二区三区| tube8黄色片| 我的女老师完整版在线观看| 97精品久久久久久久久久精品| 90打野战视频偷拍视频| 一级毛片电影观看| 精品亚洲乱码少妇综合久久| 国产极品天堂在线| 欧美人与性动交α欧美软件 | 最新的欧美精品一区二区| 夜夜骑夜夜射夜夜干| 日本黄色日本黄色录像| 爱豆传媒免费全集在线观看| 老司机影院成人| 久久av网站| 少妇高潮的动态图| 成人影院久久| 久热这里只有精品99| 国产极品粉嫩免费观看在线| 国产成人aa在线观看| www.av在线官网国产| 97在线视频观看| 日日啪夜夜爽| 只有这里有精品99|