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

    2021年夏季河南極端降雨水汽來源分析與西太平洋副熱帶高壓的可能作用*

    2023-01-17 05:54:14孫其明胡石建趙義軍
    海洋與湖沼 2023年1期
    關(guān)鍵詞:副熱帶西太平洋水汽

    孫其明 胡石建 趙義軍

    2021年夏季河南極端降雨水汽來源分析與西太平洋副熱帶高壓的可能作用*

    孫其明1, 2, 3胡石建2, 3, 4①趙義軍1

    (1. 山東科技大學(xué) 山東青島 266590; 2. 中國科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室 中國科學(xué)院海洋研究所 山東青島 266071; 3. 中國科學(xué)院大學(xué) 北京 100049; 4. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室 山東青島 266237)

    為研究西太平洋海域在河南“7·20”極端降雨事件中的作用, 利用海表面溫度、海平面氣壓、位勢高度和風(fēng)場等再分析數(shù)據(jù)與降水觀測數(shù)據(jù), 分析了2021年夏季河南極端降雨事件的水汽輸送來源及其與西太平洋副熱帶高壓多尺度變化特征的關(guān)系。數(shù)據(jù)顯示, 2021年夏季河南極端降雨期間, 西太平洋副熱帶高壓偏強(qiáng)偏北, 河南地區(qū)以東至西北太平洋中低空出現(xiàn)顯著的東風(fēng)異常和西向水汽輸運(yùn)。結(jié)果表明, 2021年夏季河南暴雨區(qū)中高空相對渦度、位勢高度和垂直上升運(yùn)動(dòng)存在顯著異常, 從西北太平洋海域通過中低空異常東風(fēng)輸送的水汽是主要來源。經(jīng)驗(yàn)?zāi)B(tài)分解結(jié)果顯示, 副熱帶高壓強(qiáng)度指數(shù)的第二模態(tài)(intrinsic mode functions of west Pacific subtropical high, IMFWPSH)與河南地區(qū)降雨第二模態(tài)存在統(tǒng)計(jì)顯著關(guān)系并在此次極端降雨事件中扮演重要角色。診斷分析發(fā)現(xiàn), 在IMFWPSH正位相期間, 副熱帶西太平洋存在異常增強(qiáng)的上升運(yùn)動(dòng)并在中低空產(chǎn)生異常東風(fēng), 河南地區(qū)上升氣流異常增強(qiáng), 因此, IMFWPSH可為河南地區(qū)降雨提供重要條件。

    河南; 極端降雨; 西太平洋副熱帶高壓; 臺(tái)風(fēng)

    極端降雨事件對社會(huì)經(jīng)濟(jì)、生產(chǎn)和生活具有極其重要的影響, 其形成過程和機(jī)制是一個(gè)被廣泛關(guān)注的研究熱點(diǎn)。在全球變暖的大背景下, 極端降雨等極端天氣事件的發(fā)生更加頻繁、影響更加顯著(Zhai, 2005; IPCC, 2019), 對極端降雨過程和物理機(jī)制的研究具有重要意義。

    我國極端降雨事件的發(fā)生受印太海溫變化、熱帶氣旋、副熱帶高壓移動(dòng)和局地陸面過程等多種復(fù)雜因素的影響。研究表明, 北印度洋和赤道中東太平洋海溫異常對我國極端降雨事件的發(fā)生頻次有影響, 當(dāng)前者冬季海溫異常偏暖時(shí), 我國華南地區(qū)夏季極端降雨次數(shù)會(huì)偏少, 而在長江流域及其以北地區(qū)極端降雨次數(shù)會(huì)偏多(Yang, 2009; 范蘇丹, 2010)。北半球冬季印度洋海溫由于自身長期記憶, 會(huì)向太平洋激發(fā)出斜壓開爾文波, 同時(shí)在南海局部海氣相互作用下向西太平洋發(fā)展, 從而對東亞大氣環(huán)流和降雨產(chǎn)生影響(He, 2015)。Studholme等(2022)指出, 由于氣候變化的影響, 颶風(fēng)和臺(tái)風(fēng)正在不斷擴(kuò)大其緯度活動(dòng)范圍, 對我國降雨將產(chǎn)生重要影響。中國東部夏季降雨異常與大尺度海溫和環(huán)流異常有關(guān)(Dong, 2016), 中緯度中太平洋夏季海溫異??梢鹬袊鴸|部地區(qū)同期降雨異常(Yu, 2007), 熱帶印度洋和太平洋海溫與副高存在相關(guān)關(guān)系, 熱帶印太海溫異常會(huì)通過調(diào)節(jié)副高的位置、面積和強(qiáng)度變動(dòng)、進(jìn)而影響中國中東部地區(qū)降雨的變化(Cai, 2010; 錢代麗, 2019; Chen, 2020)。黃櫻等(2011)利用小波分析方法, 指出除了海溫外, 南亞高壓的位置變化也會(huì)對淮河流域極端降雨產(chǎn)生重要影響。Chen等(2014a, 2014b)分析指出中國中東部持續(xù)性極端降雨事件具有雙阻塞高壓型和單阻塞高壓型兩種特征環(huán)流類型, 但在對流層上部, 南亞高壓東伸能夠促進(jìn)這兩種環(huán)流模式輻散。此外, 地形也會(huì)為極端降雨的發(fā)生起到“催化劑”的作用(Flesch, 2012; Ntwali, 2016), 比如華北中部夏季強(qiáng)降雨存在兩種類型——持續(xù)山地降雨和山地平原降雨(Sun, 2018), 太行山和伏牛山系的嵩山形成的喇叭口地形被認(rèn)為是河南“7·20”極端降雨事件產(chǎn)生的一個(gè)重要因素(冉令坤等, 2021; 張霞等, 2021; 梁旭東等, 2022)。

    近年來, 我國長江中下游地區(qū)和華北中原地區(qū)的夏季降雨異常偏高(柯丹等, 2014; Luo, 2016; Bai, 2021; Yu, 2021), 極端降雨發(fā)生的概率增加(Wei, 2010)。2020年6~7月, 長江中下游流域發(fā)生超強(qiáng)梅雨(劉蕓蕓等, 2020; 郭媛媛等, 2021), 產(chǎn)生的嚴(yán)重的洪澇導(dǎo)致150億美元以上的經(jīng)濟(jì)損失和279人死亡(翟建青等, 2021)。研究表明, 2020年超強(qiáng)梅雨的主要原因是熱帶印度洋海溫異常增暖導(dǎo)致西太平洋反氣旋加強(qiáng), 西太平洋副熱帶高壓偏強(qiáng)、西伸、脊線偏北, 南亞高壓偏強(qiáng)東伸, 異常反氣旋西北側(cè)引起西南風(fēng)異常將南海和孟加拉灣大量的水汽輸送到長江中下游地區(qū)(Fang, 2021; 郭媛媛等, 2021; Cai, 2022), 印度洋海溫異常變暖與長江中下游地區(qū)極端降雨密切相關(guān)(Tang, 2021), 梅雨期降雨具有非常顯著的年際變化(Xia, 2021)。Zheng等(2021)還指出, 2020年的極端梅雨事件與印度洋、太平洋和大西洋三大洋都有關(guān)系, 但大西洋的海溫異常起到主要作用。持續(xù)極端降雨和短時(shí)極端降雨期間, 水汽散度的變化與季節(jié)內(nèi)變化和天氣風(fēng)場異常有關(guān), 濕度場變化作用較小(Li, 2015)。

    河南省位于中國的中部、華北平原的黃河中下游地區(qū)。歷史上, 1975年8月河南就發(fā)生過特大暴雨事件(簡稱“75·8”大暴雨事件)。丁一匯(2015)總結(jié)了河南“75·8”大暴雨事件的發(fā)生機(jī)制, 指出源自東海和黃海且挾帶大量水汽的熱低壓系統(tǒng)與南下的強(qiáng)冷空氣相遇, 輻合作用劇烈、垂直上升運(yùn)動(dòng)加強(qiáng), 且受“喇叭口”地形的抬升作用, 導(dǎo)致在駐馬店一帶引發(fā)歷史罕見的強(qiáng)對流性暴雨(丁一匯等, 1978)。2016年7月河南豫北連續(xù)發(fā)生“7·09”和“7·19”兩次特大暴雨事件, 水汽均主要來自熱帶和副熱帶的暖濕氣流, “7·09”事件的整層可降水量更具有極端性, “7·19”事件垂向速度異常更顯著(王君, 2019), 栗晗等(2018)指出豫北“7·19”暴雨事件在太行山地形和低渦氣旋大尺度環(huán)流背景下, 中尺度地面輻合線的發(fā)展和維持是造成局部極端降雨的主要原因。

    2021年7月20日至21日, 河南極端降雨事件再次發(fā)生。此次事件的災(zāi)害中心是鄭州, 具有持久性、累積性、范圍廣和極端性等特點(diǎn)(蘇愛芳等, 2021), 7月20日每小時(shí)最大降水量達(dá)到201.9 mm, 鄭州國家站日累計(jì)降水量高達(dá)624.1 mm (陳剛等, 2022; 梁旭東等, 2022; 齊道日娜等, 2022; 楊浩等, 2022)。據(jù)新華社(2022)報(bào)道, 此次河南省極端降雨共造成1 478.6萬人受災(zāi), 398人死亡失蹤, 造成直接經(jīng)濟(jì)損失1 200.6億元(http://www.gov.cn/xinwen/2022-01/21/ content_5669723.htm)。

    河南“7·20”極端降雨事件受南亞高壓東擴(kuò)、大陸高壓、副高西伸北抬偏強(qiáng)和黃淮地區(qū)低渦或切變線以及臺(tái)風(fēng)“煙花”和“查帕卡”遠(yuǎn)距離水汽輸送等多尺度天氣系統(tǒng)共同影響(冉令坤等, 2021; 張霞等, 2021; 王軍等, 2022; 祝傳棟等, 2022)。西南季風(fēng)將南海的水汽輸送至西北太平洋, 通過臺(tái)風(fēng)“煙花”作用向北輸送,西北太平洋的水汽在臺(tái)風(fēng)“煙花”北部的異常偏東氣流和維持在日本上空副高的異常東南氣流共同作用下, 持續(xù)不斷地向河南地區(qū)輸送水汽(梁旭東等, 2022)。研究表明從華南到河南存在一個(gè)“西南-東北”水汽輸送帶(孫躍等, 2021), 布和朝魯?shù)?2022)通過水汽通量分析和LAGRANTO (Lagrangian Analysis Tool)模式軌跡分析發(fā)現(xiàn)在“7·20”極端降雨事件期間河南南側(cè)的經(jīng)向水汽輸送急劇增強(qiáng), 但Nie等(2022)運(yùn)用拉格朗日方法指出此次極端降雨事件中水汽主要來源是中國南方和西北太平洋地區(qū), 其貢獻(xiàn)率分別為38.1%和30.0%, 而南海貢獻(xiàn)率僅為4.2%, 而且西北太平洋海域的水汽來源可以追溯到180°E。Cheng等(2022)使用全球大氣環(huán)流的水分收支方程和三模式分解方法, 發(fā)現(xiàn)與經(jīng)向環(huán)流相比, 緯向環(huán)流在該異常降水事件中占主導(dǎo)作用。

    綜上可知, 2021年河南“7·20”極端降雨事件的成因與歷史上發(fā)生的極端事件似乎有所不同, 目前對“7·20”極端降雨事件水汽來源的結(jié)論尚有分歧。盡管前人研究已經(jīng)證明西太平洋副熱帶高壓對我國降雨具有重要影響, 但是副熱帶高壓與河南降雨之間是否以及存在何種聯(lián)系也不太清楚。因此, 本文在前人的大量深入研究基礎(chǔ)上, 分析2021年7月河南極端降雨的海洋大氣特征和水汽來源, 探索副熱帶高壓和河南降雨之間的關(guān)系和機(jī)制。第1節(jié)將介紹所用的數(shù)據(jù)和方法, 第2節(jié)將展示主要分析結(jié)果, 包括河南降雨的基本特征、“7.20”事件的水汽來源和西太平洋副熱帶高壓的作用等。第3節(jié)將簡要總結(jié)結(jié)論并對本文存在的問題進(jìn)行討論。

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

    1.1 數(shù)據(jù)

    本文使用的降雨數(shù)據(jù)包括: (1) 美國國家海洋大氣管理局氣候預(yù)測中心(National Oceanic and AtmosphericAdministration/Climate Prediction Center, NOAA/CPC)逐日降水資料, 分辨率為0.5°×0.5°; (2) 美國國家海洋大氣管理局 (Precipitation Reconstruction over Land, PREC/L)降水資料, 分辨率是1°×1°, 時(shí)間范圍為1979~2021年; (3) 中國2 479個(gè)臺(tái)站的逐日降水資料, 從中選取河南省121個(gè)臺(tái)站的逐日降水觀測資料, 時(shí)間范圍為1979~2019年(朱雪妍等, 2018)。

    大氣場資料使用的是國家環(huán)境預(yù)測中心(National Center for Environmental Prediction, NCEP)/國家大氣研究中心(National Center for Atmospheric Research, NCAR)再分析數(shù)據(jù)集Ⅰ, 包括海平面氣壓(sea level pressure, SLP)、位勢高度(geopotential height, GPH)、風(fēng)場、比濕和地表氣壓, 空間分辨率是2.5°×2.5°, 所選用的時(shí)間段都是1979~2021年。

    海表面溫度(sea surface temperature, SST)數(shù)據(jù)使用的是美國國家海洋大氣局(NOAA)提供的逐月延長重建海表面溫度(extended reconstructed sea surface temperature, ERSST) V5版本數(shù)據(jù), 分辨率是2°×2°, 時(shí)間段是1979~2021年。海表面蒸發(fā)數(shù)據(jù)采用歐洲中期天氣預(yù)報(bào)中心(European Centre for Medium-Range Weather Forecasts, ECMWF)的ERA5 (ECMWF Reanalysis v5)數(shù)據(jù), 分辨率是0.25°×0.25°, 時(shí)間段為1979~2021年。

    此外, 為了研究西太平洋副熱帶高壓的作用, 本文還采用了國家氣候中心發(fā)布的1979~2021年逐月西太副高指數(shù)(劉蕓蕓等, 2012)。

    1.2 方法

    利用風(fēng)場、比濕和地表氣壓數(shù)據(jù), 可計(jì)算整層水汽通量和整層水汽通量散度(Sun, 2011; Nie, 2022; 布和朝魯?shù)? 2022), 其中, 整層水汽通量H為

    式中,指風(fēng)速矢量,指比濕,指氣壓,t指上邊界氣壓,s指地表氣壓,為重力加速度, 取9.8 m/s2, 單位為kg/(m?s)。

    整層水汽通量可以寫成緯向u和經(jīng)向v分量的形式:

    式中,指緯向風(fēng),指經(jīng)向風(fēng)。

    由于水汽通量和風(fēng)速都屬于向量, 且方向相同, 所以水汽通量散度可由水平散度得出:

    垂直積分水平水汽通量散度可以表示:

    為研究西太副高和河南降雨的變化特征, 我們對西太副高強(qiáng)度指數(shù)和河南降雨時(shí)間序列開展了集合經(jīng)驗(yàn)?zāi)B(tài)分解(ensemble empirical mode decomposition, EEMD) (Wu等,2009), 并對EEMD分解得到的本征模函數(shù)(intrinsic mode functions, IMF)進(jìn)行分析, 研究對河南降雨影響最主要的IMF模態(tài)。

    2 結(jié)果

    2.1 河南降雨的基本特征

    圖1a和圖1b分別給出了2021年7月18~21日降水分布和2021年7月平均的降水異常(相對于1979~2020年氣候態(tài)平均), 結(jié)果可見, 華北地區(qū)和華東部分地區(qū)在2021年7月具有很強(qiáng)的降水正異常, 其中7月18~21日的降雨主要集中在河南鄭州地區(qū)。

    為便于分析, 本文將河南地區(qū)(110.5°~116.5°E, 31.5°~36.5°N)的平均降水定義為河南降水指數(shù)。如圖1c所示, 2021年7月河南降水指數(shù)是1979年以來的最高紀(jì)錄。去除長期線性趨勢后, 2021年7月的河南降水指數(shù)仍是1979年以來的最高值(圖1d)。

    根據(jù)圖1c和1d所示, 河南7月降水最強(qiáng)的三年分別是2007年、2020年和2021年。在1979~2019年期間氣候態(tài)的河南7月累積降雨量約為150 mm, 但這三年的7月累積降雨量均超過了250 mm, 在2021年達(dá)到了歷年最高的300 mm以上, 其中在7月20日之前的一周多時(shí)間內(nèi), 集中降雨量很大、短期內(nèi)超過了100多毫米(圖2)。

    圖1 2021年7月河南降水異常

    注: a: 基于CPC GLOBAL PRCP V1.0降水?dāng)?shù)據(jù)2021年7月18~21日河南的降水分布, 圖中紅色三角符號(hào)代表鄭州位置; b: 基于PREC/L數(shù)據(jù)計(jì)算的2021年7月河南降水異常, 圖中紅色方框表示河南區(qū)域(110.5°~116.5°E, 31.5°~36.5°N); c: 基于站點(diǎn)觀測數(shù)據(jù)和PREC/L數(shù)據(jù)計(jì)算的河南降水指數(shù)的標(biāo)準(zhǔn)化時(shí)間序列(柱狀圖), 黑實(shí)線為降水指數(shù)時(shí)間序列的線性趨勢; d: 去除線性趨勢后的河南降水指數(shù)的時(shí)間序列

    圖2 7月1~31日河南逐日累加降水量

    注: 基于CPC數(shù)據(jù)和河南省121個(gè)氣象站地面降水觀測數(shù)據(jù)計(jì)算的7月1~31日累加降水量

    2.2 河南“7·20”極端降雨的水汽來源

    圖3展示了2021年7月以河南為中心的亞印太區(qū)域在不同高度的大氣環(huán)流異常和水汽輸運(yùn)異常。在500和850 hPa等壓面上, 副熱帶西北太平洋存在兩個(gè)異常反氣旋(圖3a、3b), 同時(shí)在我國臺(tái)灣以東海域存在一個(gè)大范圍的異常氣旋(圖3a、3b)。在這些氣旋和反氣旋的作用下, 在30°~40°N存在著顯著的東風(fēng)異常和西向水汽輸運(yùn)異常, 且延伸至河南地區(qū)(圖3d), 河南地區(qū)呈現(xiàn)上升氣流。結(jié)合前人研究表明, 在此期間副熱帶西北太平洋的異常反氣旋, 可能與副高偏強(qiáng)、偏北且維持在日本附近海域有關(guān)系(王旭棟, 2021),而我國臺(tái)灣以東海域的異常氣旋主要是由臺(tái)風(fēng)“煙花”引起的。在500 hPa等壓面上, 異常氣旋的緯向影響范圍相較于850 hPa等壓面更大、更強(qiáng)且向北擴(kuò)展(圖3a、3b)。

    圖3f給出了500 hPa等壓面上氣候態(tài)和2021年7月的西太平洋副熱帶高壓范圍和脊線。通過和氣候態(tài)對比可知, 2021年7月西太副高比往年偏強(qiáng)、位置偏東, 脊線變化不大, 但西太副高的北邊界偏北, 其中2021年7月20~21日平均西太副高范圍偏強(qiáng)、偏北, 維持在30°N附近西太平洋海域(圖3f)。在200 hPa等壓面上, 河南地區(qū)上空為不太顯著的位勢高度負(fù)異常, 而且2021年7月南亞高壓范圍和脊線與氣候態(tài)進(jìn)行比較, 可以看出南亞高壓偏強(qiáng)東伸, 脊線在東海地區(qū)開始偏南, 呈現(xiàn)出階段性演變特征, 其中2021年7月20~21日平均南亞高壓偏強(qiáng)東伸, 北邊界北移(圖3c)。同時(shí), 在500 hPa等壓面上, 臺(tái)風(fēng)“查帕卡”和“煙花”雙臺(tái)風(fēng)的影響下導(dǎo)致顯著的垂向速度負(fù)異常, 而在日本西部和南部海域垂向速度正異常, 在河南地區(qū)出現(xiàn)顯著的垂向速度負(fù)異常, 即存在顯著的垂直上升運(yùn)動(dòng), 為極端降雨形成營造了一個(gè)重要天氣尺度條件(圖3e)。大氣環(huán)流的上述變化特征與水汽輸運(yùn)的空間特征基本一致, 即水汽從副熱帶西太平洋輸運(yùn)至河南等華東地區(qū)(圖3d)。因此, 從圖3可見, 在2021年7月, 西北太平洋的異常氣旋和反氣旋相關(guān)的東風(fēng)異常和源自副熱帶西北太平洋的強(qiáng)烈水汽輸送為產(chǎn)生極端降雨建立了大尺度條件。

    圖3 2021年7月的大氣環(huán)流異常

    注: a: 2021年7月的850 hPa位勢高度異常和風(fēng)場異常(矢量); b: 500 hPa等壓面上, 位勢高度和風(fēng)場異常(矢量); c: 200 hPa等壓面上, 位勢高度和風(fēng)場異常(矢量), 實(shí)線(12520線)和虛線分別是南亞高壓的范圍和脊線, 其中紅(藍(lán))線表示2021年7月(氣候態(tài)7月)南亞高壓的范圍及脊線, 紫線表示2021年7月20~21日平均南亞高壓的范圍; d: 2021年7月整層水汽通量異常(矢量)及其整層水汽通量散度的異常(負(fù)值代表水汽的輻合); e: 2021年7月20日的700~850 hPa垂向平均的垂向速度異常(負(fù)值代表向上運(yùn)動(dòng)); f: 500 hPa等壓面的西太副高范圍(黑色粗實(shí)線, 5880線)及脊線(虛線), 其中紅粗(藍(lán)粗)線表示2021年7月(氣候態(tài)7月)西太副高范圍及脊線, 紫粗線表示2021年7月20~21日平均西太副高范圍, 細(xì)線表示位勢高度

    冷暖空氣相遇是產(chǎn)生局部極端降雨的重要因素(劉蕓蕓等, 2020), 前人研究表明, 在2021年7月, 副熱帶高壓和大陸高壓異常偏強(qiáng), “兩高對峙”阻礙了中高緯度南下的冷空氣(梁旭東等, 2022)。由圖4可見, 在500 hPa位勢高度上, 7月20~21日期間河南受冷空氣的影響不大, 直到7月底, 河南才重新被強(qiáng)冷空氣控制(圖4)。

    前人研究表明, 西太平洋副熱帶高壓的相對強(qiáng)度與相對渦度有很好的相關(guān)性(聶俊等, 2021)。圖5給出了2021年7月的相對渦度異常。在850 hPa等壓面上, 日本附近地區(qū)出現(xiàn)相對渦度負(fù)異常, 與圖3a中的異常反氣旋對應(yīng), 而在我國臺(tái)灣以東出現(xiàn)相對渦度正異常, 正好對應(yīng)著圖3a中的異常氣旋, 而河南地區(qū)位于正渦度異常的邊緣處(圖5a)。在500 hPa和200 hPa等壓面上, 相對渦度的正異常和負(fù)異常位置分別與圖3b和圖3c中異常氣旋和異常反氣旋相一致(圖5b和圖5c)。在河南地區(qū)緯度帶(27.5°~35°N)上, 河南地區(qū)主要受相對渦度正異??刂? 而在同緯度西太平洋(大約150°E以東)則主要受相對渦度負(fù)異常影響(圖5d), 這與前文圖3所見的垂向速度異常相吻合, 表明在此期間河南地區(qū)中低空的異常垂向上升氣流進(jìn)一步促進(jìn)了極端降雨的發(fā)生。

    圖4 100°~120°E平均的500 hPa位勢高度在2021年7月的逐日演變(填色)

    注: 位勢高度5 840 m等值線(黑粗線)可指示北方冷空氣所能到達(dá)的緯度位置

    圖5 2021年7月的相對渦度異常

    注: a~c分別指2021年7月850, 500和200 hPa的相對渦度(relative vorticity, RV)異常; d: 2021年7月27.5°~35°N緯度帶上經(jīng)向平均的相對渦度異常

    通過圖3可知水汽主要受異常東風(fēng)和垂直運(yùn)動(dòng)的影響, 來自于河南以東方向的西太平洋海域。為了進(jìn)一步研究河南極端暴雨水汽輸送的具體來源, 我們在圖6中展示了30°~35°N經(jīng)向平均的風(fēng)場(,)和水汽通量隨等壓面和經(jīng)度的分布和演變。氣候態(tài)平均場(1979~2020年)顯示在該緯度帶上西太平洋中低空主要盛行西風(fēng), 其中在河南地區(qū)雖然垂向速度下沉明顯, 但主要為西風(fēng)(圖6a)。但是, 與氣候態(tài)大氣環(huán)流相比, 2021年7月的(,)異常截然不同。對比圖6a和6b可見, 2021年7月在河南至西太平洋的大范圍內(nèi)由氣候態(tài)的盛行西風(fēng)轉(zhuǎn)為東風(fēng)異常, 同時(shí)在125°~145°E、160°E~170°W的中西太平洋海域出現(xiàn)顯著的下沉氣流, 而在河南地區(qū)出現(xiàn)顯著的異常上升氣流。如圖6c所示, 在河南以東的中西太平洋大范圍內(nèi)存在顯著的西向水汽輸送。綜合起來看, 這些西向水汽輸運(yùn)從中西太平洋一直延伸到大約115°E, 并與河南地區(qū)異常強(qiáng)的上升氣流共同作用下, 對2021年7月河南地區(qū)的極端降雨的發(fā)生產(chǎn)生重要影響, 這與前人提出的水汽來源有顯著差異(王新偉等, 2018)。

    圖6 經(jīng)向平均(30°~35°N)的(u, w)和水汽通量異常

    注: a: 30°~35°N經(jīng)向平均的7月氣候態(tài)平均(1979~2020年)垂向速度(填色, 正值表示上升運(yùn)動(dòng))和緯向風(fēng)場(箭頭); b: 30°~35°N經(jīng)向平均的2021年7月垂向速度異常和緯向風(fēng)異常場; c: 30°~35°N經(jīng)向平均的2021年7月緯向水汽通量異常(填色); 黑色三角表示鄭州位置

    2.3 西太平洋副熱帶高壓的作用

    從上節(jié)中對“7·20”極端降雨的水汽來源和過程分析可知, 西太平洋副熱帶高壓可能具有非常重要的作用。比如, Yang等(2022)指出西太平洋副熱帶高壓的異常加強(qiáng)會(huì)增加熱帶到東亞的水汽輸送, 導(dǎo)致夏季極端降雨事件的發(fā)生, 未來與西太平洋副熱帶高壓異常相關(guān)的極端氣候會(huì)變得更加頻繁和嚴(yán)重。當(dāng)夏季副熱帶高壓主體北部SLP正異常, 而主體南部SLP負(fù)異常, 對應(yīng)的副高位置就會(huì)偏北(何敏等, 1991)。圖7a展示了2021年7月SLP異常和500 hPa風(fēng)場異常??梢钥闯? 在西太平洋附近, 由高緯度到低緯度, 海平面氣壓異常分布呈現(xiàn)“-+-”結(jié)構(gòu), 中緯度的副熱帶西太平洋海平面氣壓異常升高, 而在我國臺(tái)灣東部和日本南部的西太平洋地區(qū)海平面氣壓異常降低, 表明在此期間副高位置偏北, 而在異常高壓和異常低壓之間的緯度上, 亦即河南地區(qū)緯度上, 產(chǎn)生顯著的東風(fēng)異常(圖7a)。

    圖7 2021年7月的海平面氣壓(SLP)和500 hPa風(fēng)場異常以及海表面溫度(SST)和蒸發(fā)異常

    注: a: 2021年7月SLP (填色)和500 hPa風(fēng)場(箭頭)異常; b: 2021年7月SST (等值線, 單位: °C)和蒸發(fā)(填色, 正值表示蒸發(fā), 負(fù)值表示冷凝)異常

    前人研究表明, 北太平洋海溫異常與西太平洋副熱帶高壓強(qiáng)度聯(lián)系緊密, 對華北地區(qū)強(qiáng)降水的發(fā)生有重要影響(Li, 2017); 西北太平洋夏季海溫異常偏高可能會(huì)通過東亞上空東亞-太平洋遙相關(guān)型異常大氣環(huán)流影響中國東部的降水異常(Zhang, 2013)。2021年7月熱帶印度洋和熱帶太平洋的海表面溫度異常不太明顯, 但是在以大約40°N為核心的中緯度北太平洋地區(qū)海表面溫度顯著異常偏高(圖7b)、副高強(qiáng)度增強(qiáng)。2021年7月蒸發(fā)異常主要位于西太平洋緯向25°~35°N和菲律賓東部海域(120°~ 140°E), 其中在緯向帶(25°~35°N)西太平洋海域, SST異常偏高, 蒸發(fā)異常偏高, 有助于增強(qiáng)西太平洋作為水汽來源的能力。因此, 西太平洋副熱帶高壓在2021年7月河南極端降雨中可能起到非常重要的作用。

    為了進(jìn)一步研究副熱帶高壓和河南地區(qū)降雨之間的關(guān)系, 我們在圖8a中給出1979~2021年的河南逐月降水和副熱帶高壓強(qiáng)度指數(shù)的標(biāo)準(zhǔn)化時(shí)間序列。對比發(fā)現(xiàn), 河南逐月降水和副熱帶高壓強(qiáng)度指數(shù)之間存在統(tǒng)計(jì)顯著的相關(guān)性, 但其相關(guān)系數(shù)僅為0.267 (副熱帶高壓的其他指數(shù), 如面積、脊線和西伸脊點(diǎn)等,與河南降水的相關(guān)系數(shù)不如強(qiáng)度高), 在2020年和2021年, 西太平洋副熱帶高壓強(qiáng)度為正異常(圖8a)。

    西太平洋副熱帶高壓自身存在多個(gè)時(shí)間尺度的變化(如吳國雄等, 2003)。我們對河南降水和副高強(qiáng)度時(shí)間序列進(jìn)行了EEMD模態(tài)分解, 并分別對比了河南降水和副高強(qiáng)度的各個(gè)模態(tài)之間的關(guān)系。結(jié)果發(fā)現(xiàn)第一模態(tài)IMF1相對高頻, 且相關(guān)性較低, 而河南降水和副高強(qiáng)度的第二模態(tài)IMF2存在顯著的統(tǒng)計(jì)相關(guān)性, 各正負(fù)位相可較好地一一對應(yīng)起來, 二者相關(guān)系數(shù)達(dá)到0.5左右, 其中2021年河南大暴雨期間的副高強(qiáng)度IMF2模態(tài)為顯著正異常, 與河南降水IMF2在2021年的顯著正異常吻合(圖8b)。從圖9可以看出, IMF1為高頻變化項(xiàng), 周期大致為顯著的0.34 a, 而IMF2為較低頻變化項(xiàng), 其主周期為顯著的0.99 a和1.43 a。

    圖8 1979~2021年逐月河南降雨數(shù)據(jù)與西太副高(west Pacific subtropical high, WPSH)強(qiáng)度指數(shù)的時(shí)間序列及其本征模函數(shù)第二模態(tài)(the second mode of intrinsic mode functions, IMF2)的時(shí)間序列

    注: 河南降雨逐月數(shù)據(jù)和WPSH強(qiáng)度指數(shù)均進(jìn)行了13個(gè)月的滑動(dòng)平均和標(biāo)準(zhǔn)化處理, 藍(lán)線指降雨, 紅線指WPSH強(qiáng)度

    圖9 副高強(qiáng)度第一模態(tài)和第二模態(tài)的功率譜

    為研究副高強(qiáng)度指數(shù)IMF2模態(tài)對河南降雨的影響, 我們對副高強(qiáng)度指數(shù)IMF2大于其1倍標(biāo)準(zhǔn)差的年份(7月)的SST、SLP、蒸發(fā)、風(fēng)場和垂向速度場進(jìn)行了合成分析。結(jié)果顯示, 在副高強(qiáng)度指數(shù)IMF2正位相期間, 副熱帶西北太平洋海表面高度異常增強(qiáng)、SST顯著升高、海面蒸發(fā)增強(qiáng)、下沉氣流運(yùn)動(dòng)增強(qiáng), 而河南地區(qū)的垂向上升運(yùn)動(dòng)增強(qiáng), 從副熱帶西太平洋到華東地區(qū)的中低空出現(xiàn)強(qiáng)烈的東風(fēng)異常(圖10)。副高強(qiáng)度指數(shù)IMF2正位相期間的大氣環(huán)流形態(tài)是有助于為河南地區(qū)降雨提供水汽來源和垂向上升氣流條件的, 也與前文分析的2021年7月河南極端降雨期間的形態(tài)非常吻合, 說明副高強(qiáng)度指數(shù)的第二模態(tài)對河南降雨具有重要的意義。此外, Li等(2020)研究表明, 西北太平洋(大約20°~50°N, 150°~180°E)的海表溫度冷異常和熱帶太平洋的海表溫度暖異常有助于導(dǎo)致副熱帶高壓向西擴(kuò)展。從圖10a展示的SST異常場可見, 在IMF2期間西北太平洋海溫偏冷而熱帶太平洋偏暖, 這或許有助于副熱帶高壓西伸, 對河南等東部地區(qū)降雨產(chǎn)生影響。

    圖10 合成的IMF2正位相的相關(guān)物理場分布

    注: 取WPSH強(qiáng)度IMF2序列中超過標(biāo)準(zhǔn)差57.61的月份, 對相關(guān)物理量進(jìn)行合成分析. a: SST (等值線)和蒸發(fā)(填色, 正值表示蒸發(fā), 負(fù)值表示冷凝); b: 垂向速度(填色, 正值表示上升運(yùn)動(dòng))和緯向風(fēng)(箭頭); c: SLP (填色)和850~700 hPa平均風(fēng)場(矢量). 紅色實(shí)心三角表示鄭州位置

    總地來說, 副高強(qiáng)度指數(shù)的第二模態(tài)在河南降雨變化中有著非常重要的意義, 副高強(qiáng)度第二模態(tài)在2021年的極端高值可能對2021年河南極端降雨事件產(chǎn)生了重要影響。

    3 結(jié)論和討論

    極端降雨事件直接關(guān)系到人民的生命安全和財(cái)產(chǎn)安全, 2021年7月河南極端降雨給河南地區(qū)帶來了嚴(yán)重的旱澇災(zāi)害, 造成了不可估計(jì)的損失。本文利用站點(diǎn)數(shù)據(jù)、再分析數(shù)據(jù)等數(shù)據(jù)和EEMD方法, 研究了河南地區(qū)2021年7月份極端降雨的基本特征、水汽來源以及西太平洋副熱帶高壓的作用, 特別地分析了副熱帶高壓強(qiáng)度變化主要模態(tài)與河南降雨的關(guān)系。主要結(jié)論如下:

    (1) 副熱帶西太平洋是2021年7月河南地區(qū)極端降雨的主要水汽來源, 期間河南地區(qū)強(qiáng)烈的異常上升運(yùn)動(dòng)具有重要作用。在此次事件期間, 西太平洋副熱帶高壓顯著偏強(qiáng)、偏北, 在與此相關(guān)的異常反氣旋和臺(tái)風(fēng)“煙花”相關(guān)氣旋共同影響下, 副熱帶西太平洋到河南地區(qū)的大范圍內(nèi)產(chǎn)生顯著的東風(fēng)異常, 導(dǎo)致大量的水汽通過中低空從西太平洋海域輸送到河南地區(qū), 并在河南地區(qū)的異常強(qiáng)烈的低空上升氣流作用下, 形成極端降雨(圖11)。

    (2) 河南地區(qū)降水時(shí)間序列與副熱帶高壓強(qiáng)度指數(shù)之間存在統(tǒng)計(jì)顯著的相關(guān)關(guān)系, 其中副熱帶高壓的第二模態(tài)IMF2與河南地區(qū)降水時(shí)間序列的第二模態(tài)存在非常顯著的關(guān)系。在副熱帶高壓的第二模態(tài)IMF2正位相期間, 副高強(qiáng)度增強(qiáng)、SST異常偏高、蒸發(fā)異常偏強(qiáng), 中低空大氣出現(xiàn)顯著的東風(fēng)異常, 從而對河南地區(qū)降雨產(chǎn)生重要影響。在2021年7月, 副熱帶高壓的第二模態(tài)IMF2位于極端峰值, 相關(guān)海洋和大氣環(huán)流過程更加顯著, 對此期間河南地區(qū)極端降雨可能具有非常重要的作用。

    圖11 2021年7月河南極端降雨水汽來源示意圖

    注: 紅色箭頭表示異常上升氣流; 黑色箭頭表示異常矢量風(fēng); 藍(lán)色箭頭代表水汽輸送路徑; “+”表示水汽上升運(yùn)動(dòng)增強(qiáng)

    本項(xiàng)研究尚存在許多不足之處。一方面, 本文對于關(guān)于河南極端降雨的形成原因缺乏深入細(xì)致的分析。大氣環(huán)流背景場的變化和水汽來源的變化僅僅是極端降雨的必要條件而非充分條件, 只有對除了大尺度環(huán)流背景外的其他因素——比如中小尺度天氣動(dòng)力過程等, 進(jìn)行了充分的研究后才能對此次極端事件的形成原因有一個(gè)完整的理解。另一方面, 盡管本文試圖從副熱帶高壓多尺度變化模態(tài)方面研究其對河南降雨的影響, 揭示了副高強(qiáng)度第二模態(tài)可能的重要作用, 但是副高為何存在這種模態(tài)、以及其影響河南地區(qū)降雨的具體演變過程, 本文尚未研究, 則需要未來開展相關(guān)的工作。此外, 以往研究表明印度洋的暖異常在中國降雨過程中也扮演著非常重要的作用(Zhou, 2021), 印度洋的熱力變化是否以及如何在2021年7月河南極端降雨事件中發(fā)揮作用亦不清楚, 需要進(jìn)一步研究。

    致謝 感謝國家海洋環(huán)境預(yù)報(bào)中心馮立成研究員、宋春陽博士為我們提供了降雨的站點(diǎn)觀測數(shù)據(jù)和指導(dǎo), 感謝中國氣象科學(xué)研究院劉伯奇研究員、復(fù)旦大學(xué)周震強(qiáng)博士為本文提出的富有啟發(fā)性的意見和建議, 非常感謝兩位匿名審稿人提出的極具建設(shè)性的重要意見和建議。

    丁一匯, 2015. 論河南“75.8”特大暴雨的研究: 回顧與評(píng)述[J]. 氣象學(xué)報(bào), 73(3): 411-424, doi: 10.11676/qxxb2015.067.

    丁一匯, 蔡則怡, 李吉順, 1978. 1975年8月上旬河南特大暴雨的研究[J]. 大氣科學(xué), 2(4): 276-289, doi: 10.3878/j.issn. 1006-9895.1978.04.02.

    王旭棟, 2021. 夏季西北太平洋異常反氣旋的季節(jié)內(nèi)至年際尺度變化特征與機(jī)理研究[D]. 南京: 南京信息工程大學(xué): 47-60, doi: 10.27248/d.cnki.gnjqc.2021.000017.

    王軍, 吳迪, 王超杰, 等, 2022. “21·7”河南極端降水的遠(yuǎn)距離臺(tái)風(fēng)作用分析[J]. 氣象與環(huán)境科學(xué), 45(2): 75-85, doi: 10. 16765/j.cnki.1673-7148.2022.02.008.

    王君, 2019. 豫北兩次特大暴雨事件的物理量極端性和中尺度特征[J]. 干旱氣象, 37(3): 419-429, doi: 10.11755/j.issn. 1006-7639(2019)-03-0419.

    王新偉, 平凡, 姬興杰, 等, 2018. 近40a河南省夏季降水變化特征分析[J]. 暴雨災(zāi)害, 37(3): 219-229, doi: 10.3969/j.issn. 1004-9045.2018.03.004.

    布和朝魯, 諸葛安然, 謝作威, 等, 2022. 2021年“7.20”河南暴雨水汽輸送特征及其關(guān)鍵天氣尺度系統(tǒng)[J]. 大氣科學(xué), 46(3): 725-744, doi: 10.3878/j.issn.1006-9895.2202.21226.

    冉令坤, 李舒文, 周玉淑, 等, 2021. 2021年河南“7.20”極端暴雨動(dòng)、熱力和水汽特征觀測分析[J]. 大氣科學(xué), 45(6): 1366-1383, doi: 10.3878/j.issn.1006-9895.2109.21160.

    朱雪妍, 黃安寧, 2018. 不同分辨率BCC_CSM氣候系統(tǒng)模式對中國夏季極端降水事件模擬能力評(píng)估[C] // 第35屆中國氣象學(xué)會(huì)年會(huì)S25研究生論壇. 合肥: 中國氣象學(xué)會(huì): 256-263.

    劉蕓蕓, 丁一匯, 2020. 2020年超強(qiáng)梅雨特征及其成因分析[J]. 氣象, 46(11): 1393-1404, doi: 10.7519/j.issn.1000-0526. 2020.11.001.

    劉蕓蕓, 李維京, 艾秀, 等, 2012. 月尺度西太平洋副熱帶高壓指數(shù)的重建與應(yīng)用[J]. 應(yīng)用氣象學(xué)報(bào), 23(4): 414-423.

    齊道日娜, 何立富, 王秀明, 等, 2022. “7·20”河南極端暴雨精細(xì)觀測及熱動(dòng)力成因[J]. 應(yīng)用氣象學(xué)報(bào), 33(1): 1-15, doi: 10.11898/1001-7313.20220101.

    孫躍, 肖輝, 楊慧玲, 等, 2021. 基于遙感數(shù)據(jù)光流場的2021年鄭州“7·20”特大暴雨動(dòng)力條件和水凝物輸送特征分析[J]. 大氣科學(xué), 45(6): 1384-1399, doi: 10.3878/j.issn.1006- 9895.2109.21155.

    蘇愛芳, 呂曉娜, 崔麗曼, 等, 2021. 鄭州“7.20”極端暴雨天氣的基本觀測分析[J]. 暴雨災(zāi)害, 40(5): 445-454, doi: 10. 3969/j.issn.1004-9045.2021.05.001.

    楊浩, 周文, 汪小康, 等, 2022. “21·7”河南特大暴雨降水特征及極端性分析[J]. 氣象, 48(5): 571-579, doi: 10.7519/j.issn. 1000-0526.2021.111201.

    吳國雄, 丑紀(jì)范, 劉屹岷, 等, 2003. 副熱帶高壓研究進(jìn)展及展望[J]. 大氣科學(xué), 27(4): 503-517.

    何敏, 董守玉, 劉引鴿, 1991. 夏季500hPa西太平洋副高與海平面氣壓的關(guān)系[J]. 氣象, 17(2): 12-17, doi: 10.7519/j.issn. 1000-0526.1991.02.003.

    張霞, 楊慧, 王新敏, 等, 2021. “21·7”河南極端強(qiáng)降水特征及環(huán)流異常性分析[J]. 大氣科學(xué)學(xué)報(bào), 44(5): 672-687, doi: 10.13878/j.cnki.dqkxxb.20210907001.

    陳剛, 趙坤, 呂迎輝, 等, 2022. 河南“21·7”特大暴雨過程微物理特征變化分析[J]. 中國科學(xué): 地球科學(xué), 52(10): 1887-1904, zdoi: 10.1360/N072022-0029.

    范蘇丹, 2010. 中國極端降水變化的季節(jié)性差異及其與海溫異常的可能聯(lián)系[D]. 南京: 南京信息工程大學(xué): 26-30, doi: 10.7666/d.y1694863.

    柯丹, 管兆勇, 2014. 華中地區(qū)夏季區(qū)域性極端日降水事件變化特征及環(huán)流異常[J]. 氣象學(xué)報(bào), 72(3): 478-493, doi: 10.11676/qxxb2014.037.

    祝傳棟, 陳正洪, 李亞飛, 等, 2022. 平流層高位渦侵入對河南“21·7”極端暴雨的影響[J]. 氣象與環(huán)境科學(xué), 45(2): 27-37, doi: 10.16765/j.cnki.1673-7148.2022.02.004.

    聶俊, 劉鵬, 趙燦, 2021. 西北太平洋副熱帶高壓的不同指數(shù)及其與中國東部夏季降水關(guān)系的探討[J]. 大氣科學(xué), 45(4): 833-850, doi: 10.3878/j.issn.1006-9895.2009.20160.

    栗晗, 王新敏, 張霞, 等, 2018. 河南“7·19”豫北罕見特大暴雨降水特征及極端性分析[J]. 氣象, 44(9): 1136-1147, doi: 10.7519/j.issn.1000-0526.2018.09.002.

    錢代麗, 2019. 印度洋-太平洋熱帶區(qū)域異常熱力強(qiáng)迫對夏季西太平洋副熱帶高壓變化的影響機(jī)理研究[D]. 南京: 南京信息工程大學(xué): 13-24.

    郭媛媛, 張芮潔, 溫之平, 等, 2021. 熱帶西印度洋年代際增暖對2020年超強(qiáng)梅雨的可能貢獻(xiàn)[J]. 中國科學(xué): 地球科學(xué), 51(10): 1624-1636.

    黃櫻, 錢永甫, 黃丹青, 等, 2011. 全球變暖背景下淮河流域極端降水的演變規(guī)律[C] // 第28屆中國氣象學(xué)會(huì)年會(huì)—— S3天氣預(yù)報(bào)災(zāi)害天氣研究與預(yù)報(bào). 廈門: 中國氣象學(xué)會(huì): 1239-1250.

    梁旭東, 夏茹娣, 寶興華, 等, 2022. 2021年7月河南極端暴雨過程概況及多尺度特征初探[J]. 科學(xué)通報(bào), 67(10): 997-1011, doi: 10.1360/TB-2021-0827.

    新華社, 2022. 河南鄭州“7·20”特大暴雨災(zāi)害調(diào)查報(bào)告公布[J]. 中國防汛抗旱, 32(2): 5.

    翟建青, 代潭龍, 王國復(fù), 2021. 2020年全球氣候特征及重大天氣氣候事件[J]. 氣象, 47(4): 471-477, doi: 10.7519/j. issn.000-0526.2021.04.008.

    BAI S Y, GAO J X, XUE Y,, 2021. Spatial and temporal characteristics of rainfall anomalies in 1961-2010 in the Yangtze River Basin, China [J]. Atmosphere, 12(8): 960, doi: 10.3390/atmos12080960.

    CAI Y N, CHEN Z S, DU Y, 2022. The role of Indian ocean warming on extreme rainfall in central China during early summer 2020: without El Ni?o influence [J]. Climate Dynamics, 59: 951-960, doi: 10.21203/rs.3.rs-748847/v1.

    CAI Y J, TAN L X, CHENG H,, 2010. The variation of summer monsoon precipitation in central China since the last deglaciation [J]. Earth and Planetary Science Letters, 291(1/2/3/4): 21-31, doi: 10.1016/j.epsl.2009.12.039.

    CHEN Y, ZHAI P M, 2014a. Two types of typical circulation pattern for persistent extreme precipitation in Central-Eastern China [J]. Quarterly Journal of the Royal Meteorological Society, 140(682): 1467-1478, doi: 中0.1002/qj.2231.

    CHEN Y, ZHAI P M, 2014b. Precursor circulation features for persistent extreme precipitation in central-eastern China [J]. Weather and Forecasting, 29(2): 226-240, doi: 10.1175/ WAF-D-13-00065.1.

    CHEN Z H, ZHANG J, 2020. The characteristics of late summer extreme precipitation in northern China and associated large‐scale circulations [J]. International Journal of Climatology, 40(12): 5170-5187, doi: 10.1002/joc.6512.

    CHENG J B, ZHAO Y H, ZHI R,, 2022. Analysis of the July 2021 extreme precipitation in Henan using the novel moisture budget equation [J]. Theoretical and Applied Climatology, 149(1): 15-24, doi: 10.1007/s00704-022- 04022-7.

    DONG X, 2016. Influences of the Pacific Decadal Oscillation on the East Asian Summer Monsoon in non-ENSO years [J]. Atmospheric Science Letters, 17(1): 115-120, doi: 10.1002/ asl.634.

    FANG C X, LIU Y, CAI Q F,, 2021. Why does extreme rainfall occur in Central China during the summer of 2020 after a weak El Ni?o? [J]. Advances in Atmospheric Sciences, 38(12): 2067-2081, doi: 10.1007/s00376-021-1009-y.

    FLESCH T K, REUTER G W, 2012. WRF model simulation of two Alberta flooding events and the impact of topography [J]. Journal of Hydrometeorology, 13(2): 695-708, doi: 10.1175/jhm-d-11-035.1.

    HE S P, 2015. Potential connection between the Australian summer monsoon circulation and summer precipitation over Central China [J]. Atmospheric and Oceanic Science Letters, 8(3): 120-126, doi: 10.3878/AOSL20140091.

    IPCC, 2019. Summary for policymakers [C] // P?RTNER H O, ROBERTS D C, MASSON-DELMOTTE VIPCC Special Report on the Ocean and Cryosphere in a Changing Climate. Cambridge: Cambridge University Press.

    LI H X, CHEN H P, WANG H J, 2017. Influence of North Pacific SST on heavy precipitation events in Autumn over North China [J]. Atmospheric and Oceanic Science Letters, 10(1): 21-28, doi: 10.1080/16742834.2017.1237256.

    LI H, XU F H, LIN Y L, 2020. The impact of SST on the zonal variability of the western Pacific subtropical high in boreal summer [J]. Journal of Geophysical Research: Atmospheres, 125(11): e2019JD031720, doi: 10.1029/2019JD031720.

    LI R C Y, ZHOU W, 2015. Multiscale control of summertime persistent heavy precipitation events over South China in association with synoptic, intraseasonal, and low-frequency background [J]. Climate Dynamics, 45(3/4): 1043-1057, doi: 10.1007/s00382-014-2347-6.

    LUO Y L, WU M W, REN F M,, 2016. Synoptic situations of extreme hourly precipitation over China [J]. Journal of Climate, 29(24): 8703-8719, doi: 10.1175/JCLI-D-16-0057.1.

    NIE Y B, SUN J Q, 2022. Moisture sources and transport for extreme precipitation over Henan in July 2021 [J]. Geophysical Research Letters, 49(4): e2021GL097446, doi: 10.1029/2021GL097446.

    NTWALI D, OGWANG B A, ONGOMA V, 2016. The impacts of topography on spatial and temporal rainfall distribution over Rwanda based on WRF Model [J]. Atmospheric and Climate Sciences, 6(2): 145-157, doi: 10.4236/acs.2016.62013.

    STUDHOLME J, FEDOROV A, GULEV S K,, 2022. Poleward expansion of tropical cyclone latitudes in warming climates [J]. Nature Geoscience, 15(1): 14-28, doi: 10.1038/s41561-021-00859-1.

    SUN W, LI J, YU R C,, 2018. Circulation structures leading to propagating and non-propagating heavy summer rainfall in central North China [J]. Climate Dynamics, 51(9/10): 3447-3465, doi: 10.1007/s00382-018-4090-x.

    SUN B, ZHU Y L, WANG H J, 2011. The recent interdecadal and interannual variation of water vapor transport over Eastern China [J]. Advances in Atmospheric Sciences, 28(5): 1039-1048, doi: 10.1007/s00376-0100093-1.

    TANG S L, LUO J J, HE J Y,, 2021. Toward understanding the extreme floods over Yangtze River Valley in June–July 2020: role of Tropical Oceans [J]. Advances in Atmospheric Sciences, 38(12): 2023-2039, doi: 10.1007/s00376-021- 1036-8.

    WEI F Y, ZHANG T, 2010. Oscillation characteristics of summer precipitation in the Huaihe River valley and relevant climate background [J]. Science China Earth Sciences, 53(2): 301-316, doi: 10.1007/s11430-009-0151-7.

    WU Z H, HUANG N E, 2009. Ensemble empirical mode decomposition: a noise-assisted data analysis method [J]. Advances in Adaptive Data Analysis, 1(1): 1-41, doi: 10.1142/S1793536909000047.

    XIA Y C, HUANG Q, YAO S X,, 2021. Multiscale causes of persistent heavy rainfall in the Meiyu Period over the middle and lower reaches of the Yangtze River [J]. Frontiers in Earth Science, 9: 700878, doi: 10.3389/feart.2021.700878.

    YANG J H, 2009. Relationship of Pacific SSTA to summer extreme precipitation events over Eastern China [C] // 2009 3rd International Conference on Bioinformatics and Biomedical Engineering. Beijing: IEEE: 1-6, doi: 10.1109/ ICBBE.2009.5163299.

    YANG K, CAI W J, HUANG G,, 2022. Increased variability of the western Pacific subtropical high under greenhouse warming [J]. Proceedings of the National Academy of Sciences of the United States of America, 119(23): e2120335119, doi: 10.1073/pnas.2120335119.

    YU Z T, SUN Z B, NI D H,, 2007. Effects of Pacific SSTA on summer precipitation over eastern China, part II: numerical simulations [J]. Journal of Tropical Meteorology, 13(2): 113-119.

    YU R, ZHAI P M, 2021. Changes in summer persistent precipitation over the middle-lower reaches of the Yangtze River and associated atmospheric circulation patterns [J]. Journal of Meteorological Research, 35(3): 393-401, doi: 10.1007/s13351-021-0186-z.

    ZHAI P M, ZHANG X B, WAN H,, 2005. Trends in total precipitation and frequency of daily precipitation extremes over China [J]. Journal of Climate, 18(7): 1096-1108, doi: 10.1175/JCLI-3318.1.

    ZHENG J Y, WANG C Z, 2021. Influences of three oceans on record-breaking rainfall over the Yangtze River Valley in June 2020 [J]. Science China Earth Sciences, 64(10): 1607-1618.

    ZHANG R H, WU B Y, HAN J P,, 2013. Effects on summer monsoon and rainfall change over China duo to Eurasian snow cover and ocean thermal conditions [M]. In: Singh B R. Climate Change-Realities, Impacts over Ice Cap, Sea Level and Risks. Rijeka: InTech: 227-250, doi: 10.5772/54831.

    ZHOU Z Q, XIE S P, ZHANG R H, 2021. Historic Yangtze flooding of 2020 tied to extreme Indian Ocean conditions [J]. Proceedings of the National Academy of Sciences of the United States of America, 118(12): e2022255118, doi: 10.1073/pnas.2022255118.

    MOISTURE SOURCE OF EXTREME RAINFALL IN HENAN IN SUMMER 2021 AND THE POSSIBLE ROLE OF THE WESTERN PACIFIC SUBTROPICAL HIGH

    SUN Qi-Ming1, 2, 3, HU Shi-Jian2, 3, 4, ZHAO Yi-Jun1

    (1. Shandong University of Science and Technology, Qingdao 266590, China; 2. CAS Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Qingdao 266071, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China; 4. Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China)

    In order to study the role of the Western Pacific Ocean in the “7.20” extreme rainfall event in Henan Province, water vapor transportation, origination, and relationship with multiscale variation of the western Pacific subtropical high during extreme rainfall events in Henan Province, central China in summer 2021 were analyzed using reanalysis data of sea surface temperature, sea surface pressure, potential height, and wind field, and observational data of precipitation. Results show that during the extreme rainfall in Henan in summer 2021, the western Pacific subtropical high was relatively stronger and more northerly, and significant east wind anomalies and westward water vapor transport occurred in the middle and low altitude from the east of Henan to the northwest Pacific Ocean. Significant anomalies in relative vorticity, potential height, and vertical upward movement at upper level existed in the summer of 2021 in Henan rainstorm area, and the water vapor transported from the northwest Pacific Ocean via abnormal easterly winds at the middle and lower levels was the main water source. In addition, analysis with an Empirical Mode Decomposition method revealed a statistically significant relationship between the second mode of subtropical high intensity index (IMF2WPSH) and the second mode of rainfall in Henan Province, and the IMF2WPSHplays an important role in this extreme rainfall event. During the positive phase of the IMF2WPSH, there was an abnormally enhanced upward movement in the subtropical western Pacific and an abnormal east wind in the low-middle altitude, by which updraft was abnormally enhanced. The IMF2WPSHcould be used to predict extreme rainfall event in Henan and the neighboring regions.

    Henan; extreme rainfall; western Pacific subtropical high; typhoon

    P466

    10.11693/hyhz20220500121

    *中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)項(xiàng)目, XDB42010403號(hào); 國家自然科學(xué)基金項(xiàng)目, 42022040號(hào); 山東省自然科學(xué)基金資助項(xiàng)目, ZR2020JQ18號(hào); 中國科學(xué)院創(chuàng)新交叉團(tuán)隊(duì)項(xiàng)目, JCTD-2020-12號(hào)。孫其明, 碩士研究生, E-mail: sunqiming@qdio.ac.cn

    胡石建, 研究員, 博士生導(dǎo)師, E-mail: sjhu@qdio.ac.cn

    2022-05-07,

    2022-09-09

    猜你喜歡
    副熱帶西太平洋水汽
    青藏高原上空平流層水汽的時(shí)空演變特征
    2016年1月寒潮與東亞副熱帶急流關(guān)系的研究
    1979~2011年間平流層溫度及平流層水汽的演變趨勢
    西藏南部暴雪與西太平洋副熱帶高壓的關(guān)系
    西藏科技(2016年5期)2016-09-26 12:16:42
    深圳“5·11”特大暴雨過程的水汽輸送特征分析
    副熱帶東南太平洋海溫對東北夏季降水的影響及可能機(jī)制
    地基GPS觀測在水汽監(jiān)測中的應(yīng)用
    副熱帶高壓
    亚洲五月婷婷丁香| 12—13女人毛片做爰片一| 成年动漫av网址| 人人妻人人澡人人爽人人夜夜| 亚洲五月婷婷丁香| 亚洲精品久久午夜乱码| 久久久精品区二区三区| 精品一品国产午夜福利视频| 欧美精品人与动牲交sv欧美| 校园春色视频在线观看| 一边摸一边抽搐一进一出视频| 波多野结衣av一区二区av| videosex国产| 99久久国产精品久久久| 亚洲午夜理论影院| 涩涩av久久男人的天堂| 91av网站免费观看| 日本撒尿小便嘘嘘汇集6| 国产1区2区3区精品| 在线观看www视频免费| 欧美午夜高清在线| 香蕉国产在线看| 国产人伦9x9x在线观看| 一区二区三区激情视频| 91麻豆精品激情在线观看国产 | 日本一区二区免费在线视频| 91老司机精品| 久久久久久人人人人人| 国产午夜精品久久久久久| 成年人黄色毛片网站| 伦理电影免费视频| 天堂俺去俺来也www色官网| 一进一出好大好爽视频| 成年人黄色毛片网站| 两个人看的免费小视频| 欧美 日韩 精品 国产| 亚洲av美国av| 99精品久久久久人妻精品| 午夜免费鲁丝| 三上悠亚av全集在线观看| 啦啦啦 在线观看视频| 香蕉国产在线看| 淫妇啪啪啪对白视频| 成人特级黄色片久久久久久久| 久久天堂一区二区三区四区| 欧美亚洲 丝袜 人妻 在线| 高清av免费在线| 老司机亚洲免费影院| 国产精品久久久av美女十八| 亚洲精华国产精华精| 日韩欧美在线二视频 | 欧美成狂野欧美在线观看| 亚洲精品国产色婷婷电影| 成年版毛片免费区| 日韩欧美三级三区| 亚洲熟女精品中文字幕| 亚洲av成人一区二区三| 国产又爽黄色视频| 国产成人影院久久av| 高清在线国产一区| 亚洲精品国产区一区二| 日韩熟女老妇一区二区性免费视频| 久99久视频精品免费| 久久青草综合色| 黄色 视频免费看| 精品人妻熟女毛片av久久网站| 国产精品偷伦视频观看了| 国产aⅴ精品一区二区三区波| 麻豆乱淫一区二区| av有码第一页| 久久精品国产99精品国产亚洲性色 | 99在线人妻在线中文字幕 | 欧美性长视频在线观看| 又紧又爽又黄一区二区| 每晚都被弄得嗷嗷叫到高潮| 无遮挡黄片免费观看| 国产一卡二卡三卡精品| 国产精品一区二区在线观看99| 精品免费久久久久久久清纯 | 亚洲成人免费电影在线观看| 国产亚洲精品第一综合不卡| 人妻 亚洲 视频| 老汉色av国产亚洲站长工具| 国产有黄有色有爽视频| 国产精品久久视频播放| 精品一区二区三区四区五区乱码| 欧美精品av麻豆av| 搡老熟女国产l中国老女人| 亚洲专区国产一区二区| 极品少妇高潮喷水抽搐| 丝袜美腿诱惑在线| 亚洲欧洲精品一区二区精品久久久| 又紧又爽又黄一区二区| 视频区欧美日本亚洲| 很黄的视频免费| 日韩免费高清中文字幕av| 欧美av亚洲av综合av国产av| 亚洲国产精品一区二区三区在线| 国产片内射在线| 亚洲成人免费电影在线观看| 免费高清在线观看日韩| 狂野欧美激情性xxxx| 国产精品乱码一区二三区的特点 | 99久久人妻综合| 俄罗斯特黄特色一大片| 国内久久婷婷六月综合欲色啪| 亚洲欧美精品综合一区二区三区| 亚洲av成人av| 可以免费在线观看a视频的电影网站| 在线十欧美十亚洲十日本专区| 在线观看免费视频网站a站| 亚洲五月天丁香| 久久影院123| 黄色 视频免费看| 一个人免费在线观看的高清视频| 成人18禁在线播放| 亚洲av成人不卡在线观看播放网| 免费人成视频x8x8入口观看| 99精品在免费线老司机午夜| 国产一区有黄有色的免费视频| 免费人成视频x8x8入口观看| 国精品久久久久久国模美| 亚洲少妇的诱惑av| 嫁个100分男人电影在线观看| 亚洲色图综合在线观看| 精品人妻熟女毛片av久久网站| 免费在线观看黄色视频的| 亚洲国产欧美网| 国产一区有黄有色的免费视频| 69av精品久久久久久| 国产高清videossex| 成人特级黄色片久久久久久久| 岛国毛片在线播放| 大片电影免费在线观看免费| 999精品在线视频| 久久国产精品大桥未久av| 99久久国产精品久久久| 亚洲专区国产一区二区| 午夜福利在线免费观看网站| 欧美一级毛片孕妇| 国产成人欧美| 久久精品国产亚洲av高清一级| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品中文字幕在线视频| 高清黄色对白视频在线免费看| 国产高清国产精品国产三级| 老熟女久久久| 久久精品成人免费网站| 亚洲 欧美一区二区三区| 国产在线一区二区三区精| 久久99一区二区三区| 正在播放国产对白刺激| 十八禁高潮呻吟视频| 亚洲av欧美aⅴ国产| 国产在线精品亚洲第一网站| 国产在线观看jvid| 亚洲熟女精品中文字幕| 搡老熟女国产l中国老女人| 午夜福利乱码中文字幕| 99久久国产精品久久久| 亚洲五月色婷婷综合| 黄色视频不卡| 国产片内射在线| 999久久久国产精品视频| 日韩欧美免费精品| 国产aⅴ精品一区二区三区波| 欧美国产精品va在线观看不卡| 午夜福利欧美成人| av超薄肉色丝袜交足视频| 欧美黄色淫秽网站| 动漫黄色视频在线观看| 亚洲精品中文字幕一二三四区| 搡老乐熟女国产| 不卡一级毛片| 日本a在线网址| 国产高清videossex| 天天躁夜夜躁狠狠躁躁| 这个男人来自地球电影免费观看| 亚洲欧美一区二区三区久久| 99国产精品免费福利视频| 国产极品粉嫩免费观看在线| 男女高潮啪啪啪动态图| 国产成人免费无遮挡视频| 久久99一区二区三区| 熟女少妇亚洲综合色aaa.| 久久青草综合色| 亚洲精品国产区一区二| a级片在线免费高清观看视频| 亚洲国产欧美一区二区综合| 欧洲精品卡2卡3卡4卡5卡区| 丝袜美足系列| 女人高潮潮喷娇喘18禁视频| 免费观看人在逋| 亚洲综合色网址| 午夜激情av网站| 午夜免费观看网址| 久久精品亚洲av国产电影网| 18禁黄网站禁片午夜丰满| 久99久视频精品免费| 91精品国产国语对白视频| 国产精品影院久久| 亚洲中文av在线| 大香蕉久久成人网| 黑人猛操日本美女一级片| 91在线观看av| 成人特级黄色片久久久久久久| 91精品三级在线观看| 精品熟女少妇八av免费久了| 黄色丝袜av网址大全| 女人爽到高潮嗷嗷叫在线视频| 午夜视频精品福利| 天堂√8在线中文| 嫁个100分男人电影在线观看| 男女高潮啪啪啪动态图| 97人妻天天添夜夜摸| 国产又色又爽无遮挡免费看| 精品国产乱码久久久久久男人| 久久久精品国产亚洲av高清涩受| 亚洲精品av麻豆狂野| 超碰97精品在线观看| 日本精品一区二区三区蜜桃| 桃红色精品国产亚洲av| 国内久久婷婷六月综合欲色啪| 国产精品99久久99久久久不卡| 亚洲国产毛片av蜜桃av| 高清黄色对白视频在线免费看| 中出人妻视频一区二区| 激情在线观看视频在线高清 | 免费看十八禁软件| 久久人人爽av亚洲精品天堂| 黄网站色视频无遮挡免费观看| 99re6热这里在线精品视频| 99久久国产精品久久久| 精品国产乱子伦一区二区三区| 婷婷成人精品国产| 91麻豆精品激情在线观看国产 | 这个男人来自地球电影免费观看| 亚洲精品国产一区二区精华液| 欧美激情久久久久久爽电影 | 亚洲 欧美一区二区三区| 精品国产亚洲在线| 天堂动漫精品| 欧美色视频一区免费| 精品国产乱子伦一区二区三区| 女人被躁到高潮嗷嗷叫费观| 久久久水蜜桃国产精品网| 久久香蕉激情| 少妇猛男粗大的猛烈进出视频| 老熟妇仑乱视频hdxx| 国产精品.久久久| 女人精品久久久久毛片| 色老头精品视频在线观看| 午夜福利一区二区在线看| 精品国产一区二区久久| 建设人人有责人人尽责人人享有的| 国内毛片毛片毛片毛片毛片| 男男h啪啪无遮挡| 国产精品免费一区二区三区在线 | 每晚都被弄得嗷嗷叫到高潮| 高清毛片免费观看视频网站 | 亚洲,欧美精品.| 性少妇av在线| 大型黄色视频在线免费观看| 久9热在线精品视频| 97人妻天天添夜夜摸| 天堂√8在线中文| 国产精品欧美亚洲77777| 人人妻,人人澡人人爽秒播| 一区二区三区精品91| 欧美国产精品va在线观看不卡| 日韩熟女老妇一区二区性免费视频| 91大片在线观看| 欧美av亚洲av综合av国产av| 女人精品久久久久毛片| 十八禁高潮呻吟视频| 欧美色视频一区免费| 在线观看66精品国产| 亚洲中文av在线| 丰满人妻熟妇乱又伦精品不卡| 久久国产乱子伦精品免费另类| 波多野结衣一区麻豆| 在线观看66精品国产| 午夜亚洲福利在线播放| 久久国产亚洲av麻豆专区| 日韩欧美三级三区| 最近最新中文字幕大全电影3 | 国产97色在线日韩免费| av中文乱码字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 婷婷成人精品国产| 如日韩欧美国产精品一区二区三区| 97人妻天天添夜夜摸| 天堂中文最新版在线下载| 一区二区三区国产精品乱码| 日韩三级视频一区二区三区| 怎么达到女性高潮| www.999成人在线观看| 国产片内射在线| 国产不卡一卡二| 最新在线观看一区二区三区| 亚洲全国av大片| 欧美成狂野欧美在线观看| 一级a爱视频在线免费观看| 欧美日韩乱码在线| 精品国产国语对白av| 成年人免费黄色播放视频| 国产高清国产精品国产三级| 国产主播在线观看一区二区| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美亚洲国产| 欧美最黄视频在线播放免费 | 国产三级黄色录像| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 高潮久久久久久久久久久不卡| 一进一出好大好爽视频| 中文字幕最新亚洲高清| 热99久久久久精品小说推荐| 黄色丝袜av网址大全| 国产精品香港三级国产av潘金莲| 大码成人一级视频| 欧美精品高潮呻吟av久久| 亚洲成a人片在线一区二区| 自线自在国产av| 99精品久久久久人妻精品| 妹子高潮喷水视频| 两性夫妻黄色片| 国产视频一区二区在线看| 一夜夜www| 亚洲成人手机| 亚洲av电影在线进入| av国产精品久久久久影院| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 免费观看精品视频网站| 捣出白浆h1v1| 精品国内亚洲2022精品成人 | 高清在线国产一区| 久久久精品免费免费高清| 大片电影免费在线观看免费| 亚洲精品久久成人aⅴ小说| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| 十分钟在线观看高清视频www| 久久精品国产a三级三级三级| 亚洲人成电影免费在线| 日韩制服丝袜自拍偷拍| 很黄的视频免费| 丝瓜视频免费看黄片| 无人区码免费观看不卡| 日韩欧美国产一区二区入口| 人人妻人人澡人人看| 少妇粗大呻吟视频| bbb黄色大片| 亚洲 国产 在线| 欧美日韩福利视频一区二区| 亚洲全国av大片| 女警被强在线播放| av中文乱码字幕在线| aaaaa片日本免费| 国产在线一区二区三区精| 国产精品秋霞免费鲁丝片| 自线自在国产av| 日本wwww免费看| 日韩制服丝袜自拍偷拍| 久久久久久亚洲精品国产蜜桃av| 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 在线永久观看黄色视频| 精品第一国产精品| 天堂动漫精品| 美女福利国产在线| 天天添夜夜摸| 超碰成人久久| 亚洲黑人精品在线| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| 欧美乱色亚洲激情| 亚洲熟妇中文字幕五十中出 | 老司机午夜十八禁免费视频| 手机成人av网站| 99在线人妻在线中文字幕 | 啦啦啦免费观看视频1| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| 精品久久久久久久毛片微露脸| 中文字幕最新亚洲高清| 巨乳人妻的诱惑在线观看| 久99久视频精品免费| 国产亚洲精品久久久久久毛片 | 搡老熟女国产l中国老女人| 国产精品免费一区二区三区在线 | www日本在线高清视频| 亚洲av成人av| 国产精品欧美亚洲77777| videos熟女内射| 日韩有码中文字幕| av免费在线观看网站| 久久人妻熟女aⅴ| 日韩免费高清中文字幕av| 久久人妻熟女aⅴ| 欧美激情久久久久久爽电影 | 一夜夜www| 一边摸一边抽搐一进一出视频| 亚洲国产欧美一区二区综合| 在线观看www视频免费| 国产欧美日韩一区二区精品| 久久人人97超碰香蕉20202| 日本撒尿小便嘘嘘汇集6| 午夜福利视频在线观看免费| 久久久国产一区二区| 精品一品国产午夜福利视频| 91大片在线观看| 国产欧美日韩精品亚洲av| 日日夜夜操网爽| 国产精品一区二区在线观看99| 亚洲一区二区三区不卡视频| 国产精品久久久久久人妻精品电影| 少妇 在线观看| 十八禁人妻一区二区| 999精品在线视频| 99久久国产精品久久久| 亚洲欧美色中文字幕在线| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看 | 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| 美女高潮到喷水免费观看| 亚洲欧美一区二区三区久久| 亚洲一码二码三码区别大吗| 麻豆成人av在线观看| 欧美日韩视频精品一区| 国产在线一区二区三区精| 老司机午夜福利在线观看视频| 欧美黑人精品巨大| 两个人看的免费小视频| 露出奶头的视频| 久99久视频精品免费| 99久久人妻综合| 国产成人欧美在线观看 | 国产精品电影一区二区三区 | 老熟妇乱子伦视频在线观看| 在线观看免费日韩欧美大片| 亚洲av成人不卡在线观看播放网| 日韩人妻精品一区2区三区| 777久久人妻少妇嫩草av网站| 大香蕉久久成人网| 老熟女久久久| 亚洲精品国产一区二区精华液| a级片在线免费高清观看视频| 国产激情欧美一区二区| 丝袜美腿诱惑在线| 一级黄色大片毛片| 精品熟女少妇八av免费久了| 国产一区有黄有色的免费视频| 黑人猛操日本美女一级片| 99re6热这里在线精品视频| 国产99久久九九免费精品| 91成年电影在线观看| 在线十欧美十亚洲十日本专区| 免费看十八禁软件| 91字幕亚洲| 久久精品国产99精品国产亚洲性色 | 亚洲情色 制服丝袜| 成人免费观看视频高清| 热99国产精品久久久久久7| 最近最新免费中文字幕在线| 一进一出抽搐gif免费好疼 | 一边摸一边抽搐一进一出视频| 高清毛片免费观看视频网站 | 亚洲片人在线观看| 50天的宝宝边吃奶边哭怎么回事| 日韩视频一区二区在线观看| 激情在线观看视频在线高清 | 精品亚洲成国产av| 国产精品永久免费网站| 欧美日本中文国产一区发布| 无限看片的www在线观看| 视频在线观看一区二区三区| 久久久久久久国产电影| 精品午夜福利视频在线观看一区| 亚洲男人天堂网一区| 亚洲午夜理论影院| 色精品久久人妻99蜜桃| 亚洲全国av大片| 人妻久久中文字幕网| 亚洲第一青青草原| 黑人操中国人逼视频| 久久亚洲真实| 久久久久久久国产电影| 成人三级做爰电影| 久久久久精品人妻al黑| 不卡av一区二区三区| 大香蕉久久成人网| videos熟女内射| 亚洲精品国产色婷婷电影| 国产精品.久久久| 亚洲一区中文字幕在线| 成人黄色视频免费在线看| 亚洲成人免费av在线播放| 在线播放国产精品三级| 五月开心婷婷网| 色婷婷久久久亚洲欧美| 波多野结衣一区麻豆| tube8黄色片| 久久国产精品人妻蜜桃| 叶爱在线成人免费视频播放| 女同久久另类99精品国产91| 三级毛片av免费| 午夜福利乱码中文字幕| 真人做人爱边吃奶动态| 欧美精品人与动牲交sv欧美| 999久久久国产精品视频| 久久婷婷成人综合色麻豆| 亚洲成a人片在线一区二区| 另类亚洲欧美激情| av国产精品久久久久影院| 男女免费视频国产| 少妇粗大呻吟视频| 免费久久久久久久精品成人欧美视频| 日韩三级视频一区二区三区| a级毛片黄视频| 久久久国产成人免费| 巨乳人妻的诱惑在线观看| 高潮久久久久久久久久久不卡| 狠狠婷婷综合久久久久久88av| 国产精品99久久99久久久不卡| 99热国产这里只有精品6| 啦啦啦免费观看视频1| 免费观看精品视频网站| 69av精品久久久久久| 国产高清视频在线播放一区| 免费久久久久久久精品成人欧美视频| av有码第一页| 亚洲专区国产一区二区| 老司机午夜十八禁免费视频| 国产91精品成人一区二区三区| 精品卡一卡二卡四卡免费| 久久精品成人免费网站| 国产日韩欧美亚洲二区| 久久午夜综合久久蜜桃| 黄色视频,在线免费观看| 中文亚洲av片在线观看爽 | 91麻豆精品激情在线观看国产 | 国产精品久久久av美女十八| 精品国产乱子伦一区二区三区| 亚洲七黄色美女视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产中文字幕在线视频| 成人亚洲精品一区在线观看| 亚洲欧美精品综合一区二区三区| 免费在线观看亚洲国产| 宅男免费午夜| 欧美 亚洲 国产 日韩一| 成人黄色视频免费在线看| 精品一区二区三卡| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一小说 | 欧美精品av麻豆av| 香蕉丝袜av| 不卡一级毛片| 国产精品免费大片| 国产男女超爽视频在线观看| 国产精品亚洲一级av第二区| av片东京热男人的天堂| 一本一本久久a久久精品综合妖精| av电影中文网址| tocl精华| 精品一品国产午夜福利视频| 精品卡一卡二卡四卡免费| 一区二区三区精品91| 免费少妇av软件| 啦啦啦在线免费观看视频4| 最近最新免费中文字幕在线| 免费看a级黄色片| 9色porny在线观看| 99国产精品一区二区蜜桃av | 欧美日韩亚洲国产一区二区在线观看 | 亚洲七黄色美女视频| 亚洲精品一卡2卡三卡4卡5卡| 一边摸一边抽搐一进一小说 | 国产成人精品久久二区二区91| 国产精品一区二区精品视频观看| 国产黄色免费在线视频| 免费在线观看影片大全网站| 很黄的视频免费| 国产又色又爽无遮挡免费看| 成人18禁高潮啪啪吃奶动态图| 欧美日本中文国产一区发布| 午夜91福利影院| 人成视频在线观看免费观看| 国产亚洲一区二区精品| 精品国产一区二区久久| 亚洲国产看品久久| 亚洲情色 制服丝袜| 国产成人精品在线电影| 757午夜福利合集在线观看| 日韩大码丰满熟妇| 亚洲全国av大片| 婷婷丁香在线五月| 久久精品国产清高在天天线| 日日夜夜操网爽| 亚洲av第一区精品v没综合| 麻豆乱淫一区二区| 美女 人体艺术 gogo| 免费少妇av软件| 精品第一国产精品| 视频区欧美日本亚洲| 757午夜福利合集在线观看| 大香蕉久久成人网| 老熟女久久久| 黄色a级毛片大全视频| 久久久久久久久久久久大奶| 男女午夜视频在线观看|