999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

“21?7”河南極端強降水特征及環流異常性分析

2021-12-11 11:05:27張霞楊慧王新敏申琳王迪栗晗
大氣科學學報 2021年5期

張霞 楊慧 王新敏 申琳 王迪 栗晗

摘要 基于國家自動站及區域自動站降水觀測資料、歐洲中心大氣再分析資料(ERA5)分析了河南“21·7”過程的降水特征及環流和物理量的異常性,并對比1981年以來鄭州和鶴壁50 mm以上強降水過程的物理量特征。結果表明:1)“21·7”強降水過程在累計降水量、強降水覆蓋范圍、日雨量和小時雨強等方面均表現出顯著極端性,過程累計降水量超過400 mm的站點集中分布在太行山東麓沿山地區和伏牛山東側迎風坡一側,與地形關系十分密切。2)南亞高壓增強東伸,副熱帶高壓異常偏強偏北,低緯度地區活躍的低值系統等大氣環流異常,導致了水汽穩定持久向河南輸送,太行山和伏牛山沿山一帶水汽輻合偏離氣候態最強超過-10σ,表現出顯著極端性。3)“21·7”過程中,動力條件的異常特征十分顯著,200 hPa的輻散中心分別位于伏牛山和太行山東麓沿山一帶,相較歷史氣候態偏離度達到2σ~5σ;伏牛山沿山一帶850 hPa渦度偏離氣候態程度較太行山東麓一帶更大,達6σ;而700 hPa上升運動則是太行山東麓一帶極端性更強,標準差達-3σ~-5σ。4)與1981年以來同區域暴雨過程相比,“21·7”過程中,850 hPa渦度和700 hPa垂直速度的標準差為歷次過程最大(最?。┗虼未螅ù涡。┱?,對暴雨極端性有指示意義,地形附近歷次暴雨過程物理量統計顯示,伏牛山和太行山東麓的850 hPa輻合及700 hPa垂直速度平均偏離氣候態超過3σ(-3σ),且偏離程度與日雨量呈正相關。

關鍵詞極端強降水;環流異常;標準差;地形抬升

近年來,學者們對極端強降水的研究聚焦于極端強降水的氣候變化特征(高濤和謝立安,2014;孫繼松等,2015)、大氣環流異常(張玲和智協飛,2010;陳海山和陳健康,2017;張夢珂和金大超,2019)、主要影響系統的演變機制和環境場特征(趙洋洋等,2013;雷蕾等,2017;曾智琳等,2018)、中尺度對流及微物理特征(徐珺等,2018;劉晶等,2019;陳豫英等,2021)等方面,取得了較多成果。對中國1951—2010年極端降水事件的研究表明,全球變暖背景下極端降水事件的頻率和強度具有升高趨勢(高濤和謝立安,2014)。王倩等(2019)認為異常偏強、偏北的副熱帶高壓以及增強的南亞高壓與東亞地區極端降水事件直接相關。丁一匯等(2020)研究認為東亞夏季風水汽輸送的強度、影響范圍和持續性在極端強降水過程中起著關鍵的作用,環流的異常及其相伴隨的季風水汽輸送帶的北推加強,是造成北方持續性大暴雨的重要原因。孫軍等(2012)診斷了北京2012年“7·21”極端強降水過程中物理量偏離氣候平均值的程度,發現大氣可降水量、850 hPa水汽通量、850 hPa全風速偏離氣候態平均程度均達3σ以上,而850 hPa水汽通量更是達6σ。王婧羽等(2014)研究發現2012年北京“7·21”特大暴雨的產生與異常的水汽輸送密切相關,經向水汽輸送在此次暴雨過程中起主要作用。肖遞祥等(2017)對四川盆地的23次極端強降水天氣過程診斷發現,環境場具有低層高比濕、整層相對濕度大、暖云層厚、垂直風切變小等特征,850 hPa比濕和假相當位溫具有顯著正距平。王叢梅等(2017)研究認為太行山地形通過增強輻合上升運動、增大垂直風切變使雷暴下山加強,是小時雨量超過50 mm的極端短時強降水(俞小鼎,2013)天氣產生的重要原因。

河南省地處中緯度南北氣候過渡帶,西部為連綿的丘陵山地,東部為廣闊的豫東平原,地形復雜,暴雨災害頻發,降雨量的季節和區域分布極不均勻,暴雨具有強度大、突發性強、次生災害嚴重等特點。每年的6—8月隨著副熱帶高壓季節性北抬和南落,河南處于暴雨頻發時段。低渦切變線和臺風是造成河南極端強降水事件的主要影響天氣系統(張霞等,2020a),新中國成立以來,受臺風影響,河南曾先后出現了“58·7”、“75·8”、“82·8”、“96·8”和“18·8”等極端強降水事件(“58·7”暴雨研究組,1987;丁一匯,2015;趙培娟等,2019),其中“75·8”特大洪水災害世界矚目(陶詩言,1980;李澤椿等,2015)。受低渦切變線影響,2003年、2007年和2020年的7月,淮河上游河南段出現持續性強降雨天氣,致使王家壩水位超警并泄洪(張霞,2020b)。

河南省極端強降水的研究成果,集中在對極端強降水個例的多尺度天氣系統相互作用、復雜地形對強降水的影響及物理量的極端性特征等方面。河南“75·8”特大暴雨成因分析表明,行星尺度環流的調整、天氣尺度系統的相互作用、中小尺度系統的劇烈活動和有利地形作用這四方面有利條件的配合,構成了這次極端特大暴雨過程(“75·8”暴雨會戰組,1977a,1977b)。對2016年“7·9”和“7·19”豫北兩次極端強降水研究表明(馬月枝等,2017),太行山地形的動力抬升、山前中尺度地形輻合線的發展和維持是導致局地特大暴雨的主要原因;遠距離臺風登陸促成了一條伸向內陸的暖濕氣流輸送帶,并在豫北太行山東麓迎風坡喇叭口處匯聚,形成長達6 h以上超強水汽輻合傾斜上升運動,為“7·9”極端強降水的發生提供了充足的水汽條件(司福意等,2021)。極端強降水的物理量極端性特征方面,張霞等通過對1981年以來河南省13例極端強降水過程的動力、熱力因子診斷,發現大氣環流異常導致環境場的熱力、動力等物理參數異常是造成暴雨極端性的重要原因(張霞等,2020a)。冀翠華和李姝霞(2021)對比分析了2000—2016年開封22例暴雨以上降水的動力和水汽條件,發現日雨量最大的2016年“7·19”暴雨過程中,其700 hPa和850 hPa垂直速度和整層可降水量均為當地2000年以來暴雨過程之最。楊舒楠和端義宏(2020)發現臺風“溫比亞”極端強降水過程中存在極端低層輻合和高空輻散,同時假相當位溫、整層可降水量和水汽通量散度均具有顯著極端性。

2021年7月17—22日,河南省出現歷史罕見的極端強降水(簡稱“21·7”),最大累計降雨量高達1 122.6 mm,鄭州國家站最大1 h降雨量達201.9 mm。數值模式和業務預報對于此次降水中心預報有偏差,對于降水強度和極端性估計不足。那么,本次發生在我國黃淮地區的極端強降水具有哪些主要特征?其大氣環流有何異常?環境場的水汽、熱動力參數偏離氣候態達到什么程度?這些問題是本文分析和研究的重點。

1 資料與方法

1.1 資料說明

所用資料包括2021年7月17—22日河南省2 611個區域自動站逐小時降水觀測資料、河南省119個國家自動站建站以來日降水和小時降水資料(來自河南省信息中心);1951—2021年全國國家氣象站1 h降雨量資料(來自中國氣象局氣象信息中心),歐洲中期天氣預報中心第五代大氣再分析資料(ERA-5,https://cds.climate.copernicus.eu/),垂直方向1 000~100 hPa 共27層,水平分辨率0.25°×0.25°,時間間隔1 h。

使用的地形數據來自航天飛機雷達地形探測任務(Shuttle Radar Topography Mission,簡稱SRTM-90 m)數字高程模型資料(Digital Elevation Mode,簡稱DEM)。

1.2 方法介紹

1.2.1 環境場標準差倍數計算

使用環境場各物理量相對于歷史同期氣候平均值的標準差倍數表征其偏離氣候態的程度,由于不同時期環境場物理量值差異較大,為統一衡量其異常性,去除各物理量的季、月差異,文中的歷史同期氣候平均值均采用了1981—2010年30 a間暴雨日當天及其前后兩天共5 d的數據參與計算。

用于計算環境場標準差倍數公式如下:

其中:k為環境場的標準差倍數;X為計算時刻的環境參數值;為環境參數歷史同期氣候平均值;σ為計算日當天及其前后兩天共5 d的氣候態標準差。

標準差的計算采用公式(2)取得:

其中:xi為1981—2010年計算日及其前后兩天共5 d逐日物理量值;意義同式(1);n為該組數的個數。

1.2.2 副高指數計算

文中副高指數計算方法如下:

副高脊線位置:取110°~140°E范圍內副高體脊線與經線交點緯度的平均值作為副高脊線位置。

副高強度指數:對平均位勢高度大于5 880 gpm網格點的位勢高度值與5 870 gpm之差進行累計,此累計值為副高強度指數。

副高指數的多年平均指采用上述方法計算的1981—2010年30 a逐日的副高脊線位置和副高強度指數氣候平均值。

2 “21·7”強降水特征分析

7月17—22日,河南省出現了歷史罕見的特大暴雨,強降雨中心位于鄭州、鶴壁、新鄉、焦作和安陽等地,該五地市的累計平均降水量408~623 mm。最強時段為19—21日,19日夜里至20日強降雨中心位于鄭州,21日北移至豫北的鶴壁、新鄉、焦作和安陽等地。暴雨過程具有持續時間長、累計雨量大、強降雨范圍廣、短時雨量極強等特征。全省1/6國家站日降雨量突破歷史極值,最大累計降雨量高達1 122.6 mm,鄭州最大小時降雨量突破我國內陸氣象數據小時降雨量歷史極值(表1)。

極端強降雨及次生災害導致鄭州、鶴壁、新鄉、安陽等城市發生嚴重內澇,一些河流、水庫出現超警水位和保證水位,賈魯河、衛河部分河段出現漫堤潰堤,致使一些農田和村莊被淹,部分鐵路停運、航班取消,造成了嚴重的災害損失。

2.1 強降水階段性演變特征

依據降水演變特征將“21·7”極端強降水過程分為四個階段:第一階段:17日08時—19日08時,豫北沿山一帶對流性降水階段(圖1a)。該階段全省以陣性降水為主,暴雨點較分散,17日11—20時,大于等于20 mm/h的短時強降水主要出現在豫東南信陽和豫西南的南陽一帶。18日16時起,豫北太行山沿山一帶不斷有對流新生、發展合并,致使該地區持續出現強降水,安陽、鶴壁、新鄉、焦作累計降水量達50 mm以上。該階段全省有98個氣象站降水量超過100 mm,最大為新鄉輝縣的萬仙山站達235.7 mm。全省平均降水量為28 mm,鶴壁平均降雨量最大為68 mm。

第二階段:19日08時—21日08時,豫中強降水集中階段。隨著低空東風急流加強,19日起河南降水加強。19日08時至21日08時,河南省淮河以北大部地區出現暴雨,特大暴雨主要集中在以省會鄭州為中心的豫中一帶。鄭州地區連續兩天持續特大暴雨,累計平均降水量達474 mm,最大降水量出現在鄭州新密縣的白寨站860.8 mm(圖1b),超過新密站多年平均降水量(新密國家氣象站年均降水量666 mm)近200 mm。鄭州全市有180個氣象站降水量超過250 mm,500 mm以上的氣象站達104個,是鄭州有氣象記錄以來特大暴雨影響范圍最廣的一次降水。該階段是“21·7”過程中強降水覆蓋范圍最廣時段,全省2 611個氣象站的平均降水量達到了141 mm,近一半站點(1 217個氣象站,占比46.6%)降水量超過100 mm。

第三階段:強降水中心移至豫北階段(21日08時—22日08時)。21日,豫中強降水趨于減弱,特大暴雨落區北移到豫北安陽、新鄉、焦作、鶴壁四地區,降水范圍較第二階段明顯縮小,全省平均降水量51 mm。250 mm以上的特大暴雨有139站,均集中在豫北太行山沿山一帶,最大降水量為鶴壁市科創中心站達777.5 mm(圖1c),超過鶴壁市歷年年均降水量接近200 mm(鶴壁國家站年均降水量625 mm),1 h雨強最大為149.9 mm/h(新鄉市牧野鄉)。

第四階段:降水明顯減弱階段(22日08時—23日08時)。22日,豫北降水明顯減弱,分散性暴雨和大暴雨分布于豫北太行山東側沿山一帶和伏牛山區。全省平均降水量9 mm,全省共103站出現暴雨,其中24站出現大暴雨。

2.2 強降水的極端性特征

本次強降水過程在累計降水量、強降水覆蓋范圍、日雨量和小時雨強等方面均表現出顯著極端性。

17日08時至23日08時,河南省2 611個氣象站中,有1 644站(占63%)降水量超過100 mm,其中487個氣象站累計降水量超過400 mm,集中在河南省中北部。最大累計降水量出現在鶴壁市科創中心站1 122.6 mm,是當地年均降水量的近2倍(鶴壁市年均降水量625 mm)(圖2a),鄭州、輝縣等11個國家氣象站過程雨量超過本站年降水量的歷年平均值,國家站中以鄭州站過程雨量最大,達820.9 mm,遠遠超過該站641 mm的多年平均年降水量,是7月平均月降水量的5.5倍(圖2b)。

統計顯示,全省119個國家氣象站中,鄭州、新密等19站(16%)日降水量突破建站以來歷史極值;鄭州、輝縣等32站(27%)突破建站以來最大連續3 d降水量歷史極值(圖2c),其中鄭州國家站最大日降水量624.1 mm,達建站以來最大日降水量3.3倍(189.4 mm,出現在1972年7月2日),20日16—17時鄭州站1 h雨強達201.9 mm(圖2d),刷新了該站建站以來記錄,居1951年以來我國氣象數據中1 h降雨量第一位(不含港澳臺地區),遠超第二名江蘇如皋,當地1 h降雨量為161.9 mm(出現在2019年7月17日,表1)。

與河南省歷史出現的幾次強降水過程對比(表2),“21·7”暴雨除過程累計雨量(6 d累計最大1 122.6 mm)略小于“75.8”暴雨(5 d累計1 631 mm)外,日最大降水量、全省平均降水量及日雨量破歷史極值站數均超過“75·8”、“63·8”等極端強降水過程。

2.3 地形與強降水關系

豫北大暴雨與太行山地形關系密切,栗晗等(2018)對2016年“7·19”特大暴雨過程研究發現,過程累計雨量超過250 mm的站點分布在太行山東麓臨近地區,大于500 mm的強降水中心站點海拔均大于300 m。“21·7”過程中,強降水分布有類似特征。如圖3a,累計降水量在400~800 mm(黑圓點)的站點集中分布在太行山東麓臨近地區和伏牛山東側迎風坡一側,極個別站點分布在平原。其中大于等于800 mm(紅圓點)的兩個強降水中心位于豫北新鄉、鶴壁兩地區西部緊臨太行山東麓沿山一帶和鄭州西部伏牛山東側一帶,與地形關系密切。過鄭州(圖2b)極值中心同一緯度地區地形和降水分布可看出,伏牛山東側迎風坡降水量隨海拔高度升高而增強,600 mm以上的強降水多分布于海拔300~400 m高度。過科創中心(圖2c) 同一緯度地區地形和降水分布具有相似特征,600 mm以上降水大多出現在太行山東麓迎風坡一側具有一定地形高度處,且以300~900 m 高度最為集中。兩地均有山前迎風坡一側降水量遠遠超過其東側平原地區的分布特征。

3 大氣環流及物理量異常度分析

“21·7”河南極端強降水是在南亞高壓東伸、副熱帶高壓異常偏北、低緯的孟加拉灣低壓、南海臺風“查帕卡”和西太平洋臺風“煙花”共同作用下,暖濕氣流自海上穩定持續輸送至河南,并在太行山、伏牛山前輻合抬升,多系統共同作用下形成了這次持續性罕見的極端強降水。

3.1 大氣環流異常

200 hPa高空圖上,17日08時(圖略),南亞高壓脊線位于30°N,東脊點伸至105°E,日本海低壓發展后西移至朝鮮半島附近,我國華南經長江中下游伸向黃淮一帶為大片暖區,在暖平流作用下,湖北經河南伸向河北南部的高壓脊發展,與南亞高壓之間形成一低壓槽,且低槽隨高脊發展而加深。過程平均200 hPa位勢高度、風場和散度場顯示(圖4a),河南省中西部、北部上空處在高空槽前分流區中,散度達4×10-5/s,為明顯輻散區,高空輻散對低層系統發展和上升運動增強十分有利。

“21·7”過程平均500 hPa高度場上,副熱帶高壓位于我國東北地區東部到朝鮮半島附近,貝加爾湖至東北的中高緯地區為顯著正距平區,副高位置較常年同期異常偏北(圖4b),副高脊線逐日變化顯示,自7月中旬以來,副高呈持續北抬的變化趨勢,17—22日,副高脊線位于36~42°N,較常年偏北超過10個緯距(圖4c),7月以來,副高強度指數維持在6 000以上波動,較常年同期明顯偏強(圖4d);而在10~30°N的中低緯地區,則為一東西帶狀的顯著負距平區,該區域內活躍著今年的第6號臺風“煙花”和第7號臺風“查帕卡”以及孟加拉灣低壓,活躍的熱帶低值系統為河南暴雨的發生提供了源源不斷的水汽供應。

3.2 環境場物理量異常特征

3.2.1 水汽輸送異常特征

北方持續性強降雨的發生和強弱主要決定于源自熱帶的季風水汽輸送帶(丁一匯等,2020)。由于異常偏北偏強的副熱帶高壓穩定維持,其南側的臺風“煙花”受偏東氣流引導和7號臺風“查帕卡”的影響緩慢西行,向臺灣靠近,在上述環流共同影響下,“21·7”過程的水汽輸送和水汽輻合異常加強,持久維持影響河南。圖5a為19日08時—22日08時強降水集中時段1 000~300 hPa整層水汽通量積分,可以看到來自印度洋的強西南季風水汽輸送帶,輸送水汽流入臺風“查帕卡”和“煙花”,使兩臺風長時間維持,相互作用,移動緩慢,臺風“煙花”與其北側副熱帶高壓相互作用使副高外圍的偏東風加強,沿副高外圍將東南風水汽輸送至河南,在太行山和伏牛山前形成強烈輻合。

20日是鄭州地區降雨最強時段,20日08時850 hPa水汽通量顯示(圖5b),河南中部為水汽通量大值區,中心位于鄭州附近,達24 g/(cm·hPa·s),隨著低層東南風增強,該水汽大值中心略西移,穩定維持在太行山南麓與伏牛山余脈之間的喇叭口地形處,最大值達28 g/(cm·hPa·s),該時段內,水汽通量散度輻合區一直維持在鄭州及其以西伏牛山前(圖略)。

20日23時起,200 hPa高空槽東北移,低層切變線隨之緩慢東移,水汽輸送大值帶東北移至鄭州以東和豫北一帶,21日08時,水汽通量大值中心位于太行山前迎風坡一側,中心值達28 g/(cm·hPa·s)(圖5c),水汽輻合也北抬至太行山東側沿山一帶(圖略),并在此維持至22日5時,此時段強降水中心維持在豫北。

選取19日08時至21日08時以鄭州為中心的豫中強降水區域(112.4°~114.1°E,34.2°~35.2°N)和21日08時至22日08時包含科創中心在內的豫北強降水區域(113.6°~114.3°E,35.3°~36.4°N),計算了700 hPa、850 hPa和925 hPa三層的水汽通量和水汽通量散度逐小時區域最大標準差倍數(圖6),用以分析“21·7”過程水汽輸送和水汽輻合的異常程度。

19日20時起,豫中的水汽輸送三層標準差倍數均維持在1.5σ以上,尤其是20日白天,水汽輸送和輻合均異常偏強,水汽通量的標準差維持在3σ以上(圖6a),而水汽通量散度標準差自19日08時起均過-3σ,20日期間更是維持在-6σ以下(圖6b),異常偏強且持久維持的水汽輸送和水汽輻合是鄭州產生極端暴雨的關鍵因素。

豫北強降水集中時段,水汽輸送和水汽輻合同樣具有顯著偏強特征,與豫中不同的是,對流層中層700 hPa水汽通量的標準差遠大于對流層下層850 hPa和925 hPa(圖6c),而水汽輻合則反之,以對流層下層更強,偏離氣候態程度更大(圖6d),因此,豫北極端強降水與發展深厚的濕層和對流層下層水汽輻合異常偏強密切相關。

3.2.2 動力異常特征

分別計算了豫中和豫北兩個強降水集中時段的高層散度、低層渦度和垂直速度等有利降水物理量及其標準差倍數分布(計算區域選取同圖6),可以看到,19日08時至21日08時豫中強降水階段,河南上空200 hPa為明顯輻散區,輻散中心位于伏牛山東側與太行山南麓形成的喇叭口地形區域,相較歷史氣候態達到了2σ~3σ(圖7a1),異常的高空輻散有利于對流層低層產生強烈輻合,上升運動增強,850 hPa渦度場上全省均表現為輻合區,鄭州以西的伏牛山東側沿山一帶為強輻合中心,其標準差最大達到了6σ(圖7b1),850 hPa垂直速度場上,河南省中西部有強上升運動,上升運動中心的標準差倍數達-3σ以上(圖7c1)。

豫北強降水集中時段(21日08時至22日08時),河南上空200 hPa上同樣為輻散區,輻散中心位于太行山東麓與豫北交界的河北省南部,標準差倍數達3σ~5σ(圖7a2),豫北一帶散度為10×10-5/s,與氣候態偏離度不高,但對流層低層850 hPa的輻合中心位于豫北至河北南部,沿太行山呈南北帶狀分布,標準差超過2σ(圖7b2),此處的上升運動異常度更大,中心標準差達-3σ~-5σ(圖7c2)。綜上分析,豫中和豫北特大暴雨的發生均與對流層高層持續維持的強輻散、低層維持強輻合密切相關,增強的上升運動,將“煙花”外圍和副高南側輸送至豫中和豫北的水汽抬升至更高層次,形成深厚濕層并長時間輻合,利于強降水長時間維持從而形成特大暴雨。

3.3 地形與水汽和動力條件的關系

對兩個降水集中時段的水汽和動力條件分析表明,以鄭州為中心的豫中和科創中心為中心的豫北特大暴雨的產生,其水汽輻合及動力抬升條件異常偏離氣候態的區域分布在地形迎風坡或喇叭口區域,強降水分布與地形關系分析也表明,“21·7”過程累計降水量超過400 mm的站點集中分布在太行山東麓臨近地區和伏牛山東側迎風坡一側,與地形關系十分密切。

3.3.1 “21·7”過程中地形與水汽的關系

過鄭州(34.7°N)站的水汽通量和水汽通量散度剖面顯示,“煙花”臺風外圍及副高南側的水汽沿東南氣流向內陸輸送,伏牛山地形阻擋了對流層下層水汽而使其在地形迎風坡一側形成輻合,大于12 g/(cm·hPa·s)的水汽通量在地形迎風坡處向上伸至700 hPa,濕層深厚,925 hPa和850 hPa上,地形起伏處有多個水汽輻合中心存在(圖8a1);強降水集中時段豫中地形區域(112.4°~114.1°E,34.2°~35.2°N)和其東側同緯度平原區域(114.2~115.4°E,34.2~35.2°N)逐小時水汽條件變化對比顯示(圖8 b1、8 c1),地形區域的水汽輻合強度和伸展高度均明顯高于平原區域,水汽通量大值區伸展的高度和中心值具有類似特征。

豫北過科創中心(35.7°N)的水汽通量和水汽通量散度剖面圖顯示,水汽輸送的大值中心位于太行山前地形高度較低地段,而水汽輻合位置則位于太行山迎風坡地形高度稍高區域,中心在850~925 hPa之間(圖8a2),與稍高地形區降水強于較低地形區域,有地形區域降水強于平原區的特征相一致。豫北有地形區與平原區域相比,有地形區域的水汽輻合增強、輻合伸展更高的特征較豫中更明顯。

豫中和豫北的逐小時水汽輸送和水汽輻合變化顯示,兩個強降水集中時段,山前迎風坡一帶持續維持強而深厚的水汽輸送和水汽輻合,利于強降水持續而達特大暴雨。

3.3.2 “21·7”過程中地形與動力條件的關系

過豫中強降水中心的垂直速度剖面顯示,上升運動中心分布于伏牛山地形前迎風坡一側,最強中心位于600 hPa附近,向上伸展至300 hPa(圖9a1),逐小時垂直速度變化顯示(圖9b1),強降水集中時段,伏牛山前對流層一直維持強上升運動,尤以20日08—20時期間更強,中心更是超過了50×10-1Pa/s,而其東側同緯度平原區域,上升運動雖維持,但強度較地形區域明顯偏弱(圖9c1),因此,無論小時雨強、累計雨量均以地形區域更強。過豫北強降水中心的垂直速度剖面顯示,上升運動中心分布于太行山東麓迎風坡一側,強中心位于700 hPa附近,向上伸展至500 hPa(圖9a2),逐小時垂直速度變化特征與豫中相似,強降水集中時段太行山迎風坡一側維持強上升運動,且有地形區域強于平原區域的特征(圖9b2、9c2)。與豫中地形區上升運動對比,豫北地形區域的上升運動中心更強,但伸展高度較伏牛山前低。

綜上分析,伏牛山和太行山地形對沿東南氣流輸送至豫中和豫北的水汽起阻擋作用,使水汽在地形迎風坡一側匯聚;加強的東南風遇地形抬升,山前上升運動增強,使得對流層下層水汽被上升氣流抬升至對流層中層,形成深厚濕層,水汽輻合和動力上升條件均強于平原地區,有利于地形處降水增強。

3.3.3 1981年以來地形附近暴雨以上降水過程的物理量特征

為了客觀地對比“21·7”特大暴雨過程中鄭州和科創中心兩站與歷史上發生的強降水過程之間的水汽和動力特征差異,也為了發現地形附近區域強降水過程中水汽和動力條件的共性特征,本文采用1981—2021年6—8月期間的國家自動站歷史觀測雨量資料和ERA5再分析數據集,選取鄭州站和科創中心最臨近的鶴壁和淇縣站(兩站只要有一站達標準,同時達標準時取雨量最大者)50 mm以上的強降水過程,計算了鄭州為中心的豫中強降水區域(112.4°~114.1°E,34.2°~35.2°N)和包含科創中心的豫北強降水區(113.6°~114.3°E,35.3°~36.4°N)700 hPa和850 hPa比濕、整層可降水量、200 hPa散度、850 hPa渦度和700 hPa垂直速度等物理量各個暴雨日的平均標準差最大(垂直速度選最?。┍稊等鐖D10所示。

1981年以來,鄭州站6—8月共發生50 mm以上強降水日數39 d,“21·7”過程中,20日降水量達624.1 mm,排第一位,該日850 hPa的比濕標準差為2.3σ,接近所有過程的平均值,與同區域日雨量不足100 mm的過程相比,極端性并不顯著;而700 hPa的比濕標準差為4.1σ,整層水汽含量標準差達7.1σ,二者均遠大于同區域暴雨以上降水的平均值,表現出明顯極端性(圖10a),尤其是整層水汽含量,在所有過程中排在第一位。結合圖6b,由于地形阻擋和抬升,“21·7”過程中水汽輻合異常偏強且輻合長時間持續(圖6b),持久且充沛的水汽輸送和輻合利于強降水長時間維持。值得關注的是這次過程的動力條件,200 hPa散度的標準差偏離氣候態達4.5σ,高于暴雨過程的平均值,850 hPa渦度和700 hPa的垂直速度的標準差分別達7.8σ和-6.8σ,居所有過程的第一、二位,動力條件的極端性更加顯著(圖10b)。

以鶴壁為代表的豫北,1981年以來共發生45 d暴雨日(日雨量大于等于50 mm),“21·7”過程中,豫北共產生了3個暴雨日,以21日雨量最大436.1 mm,居1981年以來首位,其水汽條件特征與鄭州站不盡相同,700 hPa比濕和整層可降水量的標準差分別達4.4σ和2.9σ,均明顯高于同區域暴雨過程的平均值,但并非所有過程最大者,而850 hPa比濕標準差為1.8σ,較平均值偏低(圖10c)。從水汽條件上,并不能區分“21·7”過程與同區域其他暴雨過程。 200 hPa散度標準差達4.2σ,略高于45例暴雨平均值,850 hPa渦度標準差達7.9σ,僅略低于2016年“7.19”暴雨過程(8.5σ),居第二位,而700 hPa垂直速度達-12.6σ,居同區域所有暴雨個例首位,動力條件的極端性較水汽條件更加顯著,對極端暴雨更有指示意義。

1981年以來地形區域的暴雨過程物理量平均標準差顯示,影響豫北和豫中極端暴雨的關鍵因子具有共性特征,即動力條件較氣候態偏離程度較水汽條件更大,100 mm以上的大暴雨過程,多數個例的850 hPa渦度和700 hPa垂直速度強于平均值,這一特征說明地形區域異常增強的動力條件更有利于暴雨增強。多個例物理量異常性分析也同時顯示出強降水影響因子的復雜性,并非最極端的降水過程,其所有物理量偏離氣候態均達最大,但異常偏強的降水過程中,總有一些物理量的表現是最極端的,這一特征應在預報業務中引起足夠警惕,多物理量綜合偏離氣候態程度或可更好將極端暴雨與一般暴雨區分開來(張霞等,2020a)。

4 結論和討論

本文對“21·7”河南省持續強降水過程的降水特征及環流和物理量的異常性進行詳細分析,并與1981年以來豫中和豫北相同強降水區域物理量做對比,得到如下結論:

1)“21·7”強降水過程在累計降水量、強降水覆蓋范圍、日雨量和小時雨強等方面均表現出顯著極端性。日最大降水量、全省平均降水量及日雨量破歷史極值站數超過“75·8”、“63·8”等新中國成立以來河南的極端強降水過程。過程累計降水量超過400 mm的站點集中分布在太行山東麓臨近地區和伏牛山東側迎風坡一側,與地形關系十分密切。

2)南亞高壓增強東伸,副熱帶高壓異常偏強偏北,低緯度地區活躍著臺風“查帕卡”、“煙花”和孟灣低壓,來自印度洋的西南季風水汽輸送使“煙花”長時間維持,其北側與副高之間東風加強向河南穩定持久輸送水汽,在太行山和伏牛山前輻合,導致了豫中和豫北兩個特大暴雨中心的產生。

3)“21·7”強降水過程中,豫中和豫北對流層下層水汽通量散度均表現出輻合持續偏強、異常偏離氣候態的極端性特征;而水汽通量波動性大,極端特征不如水汽通量散度顯著。伏牛山和太行山東麓沿山一帶的動力條件偏離氣候態程度大,200 hPa散度的標準差達2σ~5σ;豫中伏牛山沿山一帶850 hPa渦度偏離氣候態程度較太行山東麓一帶更大,達6σ,而700 hPa上升運動則是豫北太行山東麓一帶更極端,標準差達-3σ~-5σ。

4)1981年以來豫中和豫北同區域暴雨以上量級降水的物理量對比分析顯示,“21·7”過程中,850 hPa渦度和700 hPa的垂直速度標準差為歷次過程最大(最?。┗虼未螅ù涡。┱撸瑯O端性顯著,可明顯區分同區域其他暴雨過程;豫中整層水汽含量居所有暴雨過程首位,對暴雨極端性有指示意義,而豫北水汽條件與其他過程相比排位并不靠前。所選地形附近歷次暴雨過程的850 hPa渦度和700 hPa的垂直速度平均標準差偏離氣候態超過±3σ,且大暴雨以上降水過程中,多數個例的850 hPa渦度和700 hPa垂直速度強于平均值,表明地形區域動力作用的增強與暴雨增幅有明顯正相關。

本文僅對“21·7”強降水過程特征、大氣環流和環境物量偏離氣候態的程度做了較詳細分析,研究所得可為極端暴雨預報業務借鑒。然而,鄭州超過200 mm/h的極端雨強、接近年降水量的日雨量與城市發展及熱動力效應有何關聯?地形和復雜下墊面對暴雨的增幅和影響機制?極端強降水的可預報性及預報技術等等諸多科學問題值得進一步研究。

參考文獻(References)

陳海山,陳健康,2017.東亞夏季風指數的分類及物理特征分析[J].大氣科學學報,40(3):299-309. Chen H S,Chen J K,2017.Classification of East Asian summer monsoon indicesand their basic physical features[J].Trans Atmos Sci,40(3):299-309.doi:10.13878/j.cnki.dqkxxb.20160323001.(in Chinese).

陳豫英,蘇洋,楊銀,等,2021.賀蘭山東麓極端暴雨的中尺度特征[J].高原氣象,40(1):47-60. Chen Y Y,Su Y,Yang Y,et al.,2021.The mesoscale characteristics of extreme rainstorm in the eastern region of Helan mountain[J].Plateau Meteor,40(1):47-60.doi:10.7522/j.issn.1000-0534.2020.00012.(in Chinese).

丁一匯,2015.論河南“75·8”特大暴雨的研究:回顧與評述[J].氣象學報,73(3):411-424. Ding Y H,2015.On the study of the unprecedented heavy rainfall in Henan Province during 4-8 August 1975[J].Acta Meteorol Sin,73(3):411-424.doi:10.11676/qxxb2015.067.(in Chinese).

丁一匯,柳艷菊,宋亞芳,2020.東亞夏季風水汽輸送帶及其對中國大暴雨與洪澇災害的影響[J].水科學進展,31(5):629-643. Ding Y H,Liu Y J,Song Y F,2020.East Asian summer monsoon moisture transport belt and its impact on heavy rainfalls and floods in China[J].Adv Water Sci,31(5):629-643.doi:10.14042/j.cnki.32.1309.2020.05.001.(in Chinese).

高濤,謝立安,2014.近50年來中國極端降水趨勢與物理成因研究綜述[J].地球科學進展,29(5):577-589. Gao T,Xie L A,2014.Study on progress of the trends and physical causes of extreme precipitation in China during the last 50 years[J].Adv Earth Sci,29(5):577-589.doi:10.11867/j.issn.1001-8166.2014.05.0577.(in Chinese).

冀翠華,李姝霞,2021.開封“7·19”大暴雨天氣過程的極端性分析[J].氣象與環境科學,44(4):53-62. Ji C H,Li S X,2021.Analysis on extremity of the “7·19”severe rainstorm process in Kaifeng[J].Meteor Environ Sci,44(4):53-62.doi:10.16765/j.cnki.1673-7148.2021.04.007.(in Chinese).

雷蕾,孫繼松,何娜,等,2017.“7·20”華北特大暴雨過程中低渦發展演變機制研究[J].氣象學報,75(5):685-699. Lei L,Sun J S,He N,et al.,2017.A study on the mechanism for the vortex system evolution and development during the torrential rain event in North China on 20 July 2016[J].Acta Meteorol Sin,75(5):685-699.doi:10.11676/qxxb2017.054.(in Chinese).

栗晗,王新敏,張霞,等,2018.河南 “7·19” 豫北罕見特大暴雨降水特征及極端性分析[J].氣象,44(9):1136-1147. Li H,Wang X M,Zhang X,et al.,2018.Analysis on extremity and characteristics of the 19 July 2016 severe torrential rain in the north of Henan Province[J].Meteor Mon,44(9):1136-1147.doi:10.7519/j.issn.1000-0526.2018.09.002.(in Chinese).

李澤椿,諶蕓,張芳華,等,2015.由河南“75·8”特大暴雨引發的思考[J].氣象與環境科學,38(3):1-12. Li Z C,Chen Y,Zhang F H,et al.,2015.Consideration by “75·8” extreme heavy rainfall event in Henan[J].Meteorological and Environmental Sciences,38(3):1-12(in Chinese).

劉晶,周雅蔓,楊蓮梅,等,2019.伊犁河谷"極端暴雨過程不穩定性及其觸發機制研究[J].大氣科學,43(6):1204-1218. Liu J,Zhou Y M,Yang L M,et al.,2019.The instability and its trigger mechanism of extreme precipitation event in the Yili river valley on 31 July 2016[J].Chin J Atmos Sci,43(6):1204-1218.doi:10.3878/j.issn.1006-9895.1901.18155.(in Chinese).

馬月枝,張霞,胡燕平,2017.2016年7月9日新鄉暖區特大暴雨成因分析[J].暴雨災害,36(6):557-565. Ma Y Z,Zhang X,Hu Y P,2017.Cause analysis of a warm-sector excessive heavy rainfall event in Xinxiang on 9 July 2016[J].Torrential Rain Disasters,36(6):557-565.doi:10.3969/j.issn.1004-9045.2017.06.009.(in Chinese).

司福意,張一平,趙海青,等,2021.豫北“7·9”特大暴雨的形成與極端性探析[J].氣象與環境科學,44(1):65-73. Si F Y,Zhang Y P,Zhao H Q,et al.,2021.Analysis on formation and extremity of the 9 July extreme torrential rain in the north of Henan Province[J].Meteor Environ Sci,44(1):65-73.doi:10.16765/j.cnki.1673-7148.2021.01.008.(in Chinese).

孫繼松,雷蕾,于波,等,2015.近10年北京地區極端暴雨事件的基本特征[J].氣象學報,(4):609-623. Sun J S,Lei L,Yu B,et al.,2015.The fundamental features of the extreme severe rain events in the recent 10 years in the Beijing area[J].Acta Meteorol Sin,(4):609-623.doi:10.11676/qxxb2015.044.(in Chinese).

孫軍,諶蕓,楊舒楠,等,2012.北京721特大暴雨極端性分析及思考(二)極端性降水成因初探及思考[J].氣象,38(10):1267-1277. Sun J,Chen Y,Yang S N,et al.,2012.Analysis and thinking on the extremes of the 21 July 2012 torrential rain in Beijing part Ⅱ:preliminary causation analysis and thinking[J].Meteor Mon,38(10):1267-1277.doi:10.7519/j.issn.1000-0526.2012.10.013.(in Chinese).

陶詩言,1980.中國之暴雨[M].北京:科學出版社:147-161. Tao S Y,1980.Torrential rain in China[M].Beijing:Science Press:147-161.(in Chinese).

王叢梅,俞小鼎,李芷霞,等,2017.太行山地形影響下的極端短時強降水分析[J].氣象,43(4):425-433. Wang C M,Yu X D,Li Z X,et al.,2017.Investigation of extreme flash-rain events on the impact of Taihang mountain[J].Meteor Mon,43(4):425-433.doi:10.7519/j.issn.1000-0526.2017.04.004.(in Chinese).

王婧羽,崔春光,王曉芳,等,2014.2012年7月21日北京特大暴雨過程的水汽輸送特征[J].氣象,40(2):133-145. Wang J Y,Cui C G,Wang X F,et al.,2014.Analysis on water vapor transport and budget of the severe torrential rain over Beijing region on 21 July 2012[J].Meteor Mon,40(2):133-145.doi:10.7519/j.issn.1000-0526.2014.02.001.(in Chinese).

王倩,翟盤茂,余榮,2019.2018年7月北半球極端天氣氣候事件及環流特征分析[J].大氣科學學報,42(1):28-35. Wang Q,Zhai P M,Yu R,2019.Analysis of extreme weather and climate events and circulation characteristics in the Northern Hemisphere in July 2018[J].Trans Atmos Sci,42(1):28-35.doi:10.13878/j.cnki.dqkxxb.20181129001.(in Chinese).

肖遞祥,楊康權,俞小鼎,等,2017.四川盆地極端暴雨過程基本特征分析[J].氣象,43(10):1165-1175. Xiao D X,Yang K Q,Yu X D,et al.,2017.Characteristics analyses of extreme rainstorm events in Sichuan basin[J].Meteor Mon,43(10):1165-1175.doi:10.7519/j.issn.1000-0526.2017.10.001.(in Chinese).

徐珺,畢寶貴,諶蕓,等,2018.“5·7”廣州局地突發特大暴雨中尺度特征及成因分析[J].氣象學報,76(4):511-524. Xu J,Bi B G,Chen Y,et al.,2018.Mesoscale characteristics and mechanism analysis of the unexpected local torrential rain in Guangzhou on 7 May 2017[J].Acta Meteorol Sin,76(4):511-524.doi:10.11676/qxxb2018.016.(in Chinese).

楊舒楠,端義宏,2020.臺風溫比亞(1818)降水及環境場極端性分析[J].應用氣象學報,31(3):290-302. Yang S N,Duan Y H,2020.Extremity analysis on the precipitation and environmental field of typhoon rumbia in 2018[J].J Appl Meteor Sci,31(3):290-302.doi:10.11898/1001-7313.20200304.(in Chinese).

俞小鼎,2013.短時強降水臨近預報的思路與方法[J].暴雨災害,32(3):202-209. Yu X D,2013.Nowcasting thinking and method of flash heavy rain[J].Torrential Rain Disasters,32(3):202-209.doi:10.3969/j.issn.1004-9045.2013.03.002.(in Chinese).

曾智琳,諶蕓,王東海,2020.2018年8月華南超歷史極值降水事件的觀測分析與機理研究[J].大氣科學,44(4):695-715. Zeng Z L,Chen Y,Wang D H,2020.Observation and mechanism analysis for a record-breaking heavy rainfall event over Southern China in August 2018[J].Chin J Atmos Sci,44(4):695-715.doi:10.3878/j.issn.1006-9895.1906.18265.(in Chinese).

張玲,智協飛,2010.南亞高壓和西太副高位置與中國盛夏降水異常[J].氣象科學,30(4):438-444. Zhang L,Zhi X F,2010.South Asian high and the subtropical western Pacific high and its relation to the mid-summer precipitation anomalies over China[J].Sci Meteor Sin,30(4):438-444.doi:10.3969/j.issn.1009-0827.2010.04.002.(in Chinese).

張夢珂,金大超,2019.2015年夏季東亞和南亞降水異常與熱帶太平洋海溫異常的聯系及機理[J].大氣科學學報,42(5):745-754. Zhang M K,Jin D C,2019.Relationship of East and South Asia summer precipitation anomalies with tropical Pacific SSTA in 2015 and its mechanisms[J].Trans Atmos Sci,42(5):745-754.doi:10.13878/j.cnki.dqkxxb.20190322001.(in Chinese).

張霞,王新敏,栗晗,等,2020a.基于環境參數的極端暴雨指數構建及其應用[J].氣象,46(7):898-912. Zhang X,Wang X M,Li H,et al.,2020a.Construction and application of extreme rainstorm index based on environmental parameters[J].Meteor Mon,46(7):898-912.doi:10.7519/j.issn.1000-0526.2020.07.003.(in Chinese).

張霞,趙培娟,張淵萌,等,2020b.淮河上游致洪暴雨對王家壩水位影響分析及估測[J].氣象與環境科學,43(4):1-10. Zhang X,Zhao P J,Zhang Y M,et al.,2020b.Analysis and estimation of the impact of flood-causing rainstorm on the Wangjiaba dam water level in the upper reaches of Huaihe river[J].Meteor Environ Sci,43(4):1-10.doi:10.16765/j.cnki.1673-7148.2020.04.001.(in Chinese).

趙培娟,邵宇翔,張霞,2019.相似路徑臺風“摩羯”“溫比亞”登陸后環境場對比分析[J].氣象與環境科學,42(3):17-28. Zhao P J,Shao Y X,Zhang X,2019.Comparative analysis of the after landed environmental field of the similar track typhoons yagi and rumbia[J].Meteor Environ Sci,42(3):17-28.doi:10.16765/j.cnki.1673-7148.2019.03.003.(in Chinese).

趙洋洋,張慶紅,杜宇,等,2013.北京“特大暴雨環流形勢極端性客觀分析[J].氣象學報,71(5):817-824. Zhao Y Y,Zhang Q H,Du Y,et al.,2013.Objective analysis of the extreme of circulation patterns during the 21 July 2012 torrential rain event in Beijing[J].Acta Meteorol Sin,71(5):817-824.doi:10.11676/qxxb2013.066.(in Chinese).

“75·8”暴雨會戰組,1977a.河南“75·8”特大暴雨成因的初步分析(一)[J].氣象,3(7):3-5. “75·8”Rainstorm Battle Group,1977a.Preliminary analysis on the formation of “75·8” torrential rain in Henan(Ⅰ)[J].Meteor Mon,3(7):3-5(in Chinese).

“75·8”暴雨會戰組,1977b.河南“75·8”特大暴雨成因的初步分析(二)[J].氣象,3(7):6-8. “75·8” Rainstorm Battle Group,1977b.Preliminary analysis on the formation of “75·8” torrential rain in Henan(Ⅱ)[J].Meteor Mon,3(7):6-8(in Chinese).

“58·7”暴雨研究組,1987.黃河中游“58·7”大暴雨成因的天氣學分析[J].大氣科學,11(1):100-107. Research Group of “58·7” Heavy Rainfall,1987.A synoptic study of the “58·7” persistent rainstorm over the middle Huanghe (yellow) River[J].Chine J Atmos Sci,11(1):100-107.

Analysis on characteristic and abnormality of atmospheric circulations of the July 2021 extreme precipitation in Henan

ZHANG Xia1,2,YANG Hui2,WANG Xinmin1,2,SHEN Lin3,WANG Di2,LI Han2

1Henan Key Laboratory of Agrometeorological Support and Applied Technique,CMA,Zhengzhou 450003,China;

2Henan Meteorological Observatory,Zhengzhou 450003,China;

3Yangtze Ecology and Environmental Co.,Ltd.,Wuhan 430062,China

The severe torrential rain attacked Henan province during July 17—22 in 2021 (abbreviation “21·7”) whenever the process total rainfall amount reaches 1 122.6 mm and the hourly rain intensity at Zhengzhou national observation is up to 201.9 mm/h which broke through the meteorological record of hourly rain intensity extreme value in Chinese inland since 1951.The precipitation features and abnormality of atmospheric circulations and the environmental physical quantities were analyzed with the comparative analysis of physical quantities in the process of heavy rain over 50 mm in Zhengzhou and Hebi since 1981 based on the national and regional automatic meteorological observation data and the high resolution European central atmospheric reanalysis (ERA5).As shown by the results,the “21·7” severe rainfall process presents significant extremes in terms of cumulative precipitation,heavy rainfall coverage,daily rainfall,and hourly rainfall intensity.The eastern foot of Taihang Mountain and the windward slope of eastern Funiu Mountain are the only two locations with accumulated precipitation over 400 mm,which are closely related to topography.Anomalies in atmospheric circulation such as the enhancement and eastward extension of the South Asian high,the northward anomaly of the subtropical high,and the low-value system in low-latitude regions are responsible for the stable and lasting water vapor transport to Henan province.Water vapor convergence along Taihang and Funiu Mountain was below -10σ,showing significant extremes.Dynamic conditions are also apparent in the “21·7” process.With a deviation of 2σ—5σ from the historical climate,the divergence centers are located along the eastern foothills of Taihang Mountain and Funiu Mountain in 200 hPa.The deviation of 850 hPa vorticity from the climate state along Funiu Mountain is higher than that on the eastern foot of Taihang Mountain while the deviation of 700 hPa vertical velocity from the climate is more extreme on the eastern foot.For the same region since 1981,the standard deviation of 850 hPa vorticity and 700 hPa vertical velocity is the largest (least) or the secondary maximum (minimum) in the process of “21·7”,which indicates the extreme rain.In the eastern foot of Taihang Mountain and Funiu Mountain,the deviations of 850 hPa convergence and 700 hPa vertical velocity are consistently greater than 3σ(-3σ) on average,positively correlated with daily rainfall.

extreme precipitation;circulation anomaly;standard deviation;terrain elevation

doi:10.13878/j.cnki.dqkxxb.20210907001

(責任編輯:劉菲)

主站蜘蛛池模板: 美女内射视频WWW网站午夜 | 无码免费视频| 老色鬼欧美精品| 亚洲人成影院在线观看| 色吊丝av中文字幕| 亚洲欧美成aⅴ人在线观看| 亚洲区一区| 国模私拍一区二区三区| 无码专区国产精品一区| 欧美一区精品| 99热国产这里只有精品无卡顿" | 国产精品成人免费综合| 色哟哟色院91精品网站| 亚洲成年人片| 亚洲高清无码久久久| 自慰网址在线观看| 亚洲午夜福利在线| 一区二区理伦视频| 日韩麻豆小视频| 亚洲日韩日本中文在线| 色综合a怡红院怡红院首页| 国产精品久久久久无码网站| 中文字幕不卡免费高清视频| 波多野一区| 日韩毛片免费视频| 成人在线第一页| 亚洲男人在线天堂| 国产在线小视频| www精品久久| 亚洲第一视频网| 九九线精品视频在线观看| 亚洲αv毛片| 五月天福利视频| 国产欧美中文字幕| 久久久久久久久久国产精品| 国产在线小视频| 国产乱视频网站| 香蕉国产精品视频| 久久国产乱子伦视频无卡顿| 亚洲九九视频| 国外欧美一区另类中文字幕| 国产高清在线观看| 伊人成色综合网| 九色91在线视频| 日韩国产 在线| 日韩免费毛片视频| 国内精品久久久久鸭| 久久久波多野结衣av一区二区| 欧美一区二区三区国产精品| 蜜臀AV在线播放| 亚卅精品无码久久毛片乌克兰| 熟女视频91| 亚洲区第一页| 视频一区视频二区中文精品| 国产视频自拍一区| 欧美第一页在线| 欧洲极品无码一区二区三区| 91色老久久精品偷偷蜜臀| 亚洲日韩精品欧美中文字幕| 久久精品国产免费观看频道| 国产一级做美女做受视频| 高清精品美女在线播放| 拍国产真实乱人偷精品| 午夜在线不卡| 日韩欧美在线观看| 国模极品一区二区三区| 极品尤物av美乳在线观看| 国产第四页| 国产本道久久一区二区三区| 国产区精品高清在线观看| 欧美精品一区二区三区中文字幕| 久久这里只有精品2| 日韩欧美中文| 久久香蕉国产线看观| 国产乱人免费视频| 91九色最新地址| 免费A级毛片无码免费视频| 欧美不卡视频在线观看| 五月婷婷亚洲综合| a国产精品| 亚洲自偷自拍另类小说| 91福利免费视频|