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

    2021 年河南“7.20”極端暴雨動、熱力和水汽特征觀測分析

    2021-12-14 07:18:44冉令坤李舒文周玉淑楊帥馬淑萍周括申冬冬焦寶峰李娜
    大氣科學 2021年6期
    關(guān)鍵詞:云團中尺度低層

    冉令坤 李舒文 周玉淑 楊帥 馬淑萍 周括 申冬冬 焦寶峰 李娜

    1 中國科學院大氣物理研究所,北京 100029

    2 中國科學院大學,北京 100049

    1 引言

    河南地處我國中部,是我國暴雨災害較為嚴重的地區(qū)。河南省地勢西高東低,北、西、南三面被太行山、伏牛山、桐柏山和大別山環(huán)繞,西南部為南陽盆地,由于地形抬升作用,河南暴雨多分布在太行山迎風坡和山區(qū),具有強度大、持續(xù)時間長、突發(fā)性顯著等特點,經(jīng)常引發(fā)洪澇、泥石流和滑坡等地質(zhì)災害,嚴重威脅人民生命財產(chǎn)安全。1975年8 月河南南部及淮河上游發(fā)生罕見特大暴雨(“75·8”暴雨),最大過程雨量達到1631 mm,造成水庫大壩垮塌,數(shù)十萬人員傷亡(陶詩言,1980);2007 年三門峽市“7·29”大暴雨引發(fā)洪澇和山體滑坡,造成89 人死亡和失蹤,轉(zhuǎn)移安置5.7 萬人,造成經(jīng)濟損失14.1 億元(任軼, 2013);2016 年7 月9 日河南北部新鄉(xiāng)地區(qū)發(fā)生特大暴雨,暴雨中心7 h 累計雨量達到420.1 mm,共造成32.90 萬人受災,直接經(jīng)濟損失2.30 億元(徐姝等,2019; 劉冠華和??〕? 2020)。有研究表明(張遠, 2014),河南暴雨日和大暴雨日數(shù)均呈弱增加的趨勢。

    很多學者對導致河南暴雨的多種天氣尺度系統(tǒng)相互作用、中尺度系統(tǒng)活動以及地形效應開展深入研究。河南暴雨主要受到江淮氣旋、江淮倒槽、中低層低槽、鋒面、低渦和切變線等多尺度天氣系統(tǒng)的影響(梁鈺等, 2020)。丁一匯等(1978)研究指出天氣尺度系統(tǒng)為“75·8”暴雨提供有利條件,中尺度系統(tǒng)和地形對觸發(fā)暴雨起著重要作用;暴雨反過來又影響天氣尺度溫壓場和風場分布, 改變后的氣象要素場使暴雨發(fā)展更強烈。蘇愛芳等(2016)指出低渦切變線、西南急流和南下弱冷空氣的共同作用是產(chǎn)生2010 年“8·13”黃淮北部暴雨的主要原因。胡燕平等(2009)分析認為,2008 年7 月22 日河南黃淮地區(qū)暴雨是高低空急流耦合區(qū)內(nèi)西南渦沿切變線移出,弱冷空氣侵入暖倒槽觸發(fā)不穩(wěn)定能量釋放造成的。李博等(2009)分析了2007年6 月19~20 日入梅期首場強降雨,指出該降雨首先產(chǎn)生在湖北、河南南部地區(qū),隨后迅速北跳至山東南部,500 hPa 橫槽對首場強降水的維持起重要作用,低層渦旋系統(tǒng)沿著鋒生帶在槽前氣流引導下向東北方向移動,致使雨帶北跳。李宏宇等(2006)分析2002 年4 月4~5 日河南省冷鋒降水過程,指出云系存在明顯的“播種供給”云降水機制,鋒區(qū)附近擬降水效率達40%~60%,在積層混合區(qū)擬降水效率可高達70%以上。

    中尺度系統(tǒng)是造成暴雨的直接影響系統(tǒng),人們對中尺度系統(tǒng)結(jié)構(gòu)、發(fā)生發(fā)展動力學問題開展深入研究和探索(矯梅燕等, 2006)。栗晗等(2018)對2016 年“7·19”暴雨過程進行了分析,結(jié)果表明,中尺度輻合線發(fā)展維持、多個地面中尺度氣旋移動造成的列車效應是導致局地特大暴雨的主要原因。牛淑貞等(2001)分析指出2000 年7 月5~6日鄭州特大暴雨具有典型的中尺度擾動及強對流特征。覃丹宇(2010)對中尺度云團形態(tài)進行分析,利用改進Shuman-Shapiro 濾波方法有效地分離出中尺度擾動,結(jié)果表明,近圓形中尺度云團主要由低層中尺度氣旋組織形成,而寬帶狀云團主要由低層中尺度切變線組織產(chǎn)生。孫建華等(2004)對2002 年中國暴雨試驗期間南陽盆地一次中尺度對流系統(tǒng)進行分析和模擬,結(jié)果表明,由于對流層中層正渦度中心的強迫和潛熱釋放產(chǎn)生了氣旋,低層的暖平流可能是低渦東移發(fā)展的原因之一。

    河南省西高東低的地形地貌對暴雨觸發(fā)和增幅發(fā)揮重要作用。蔡則怡和宇如聰(1997)利用ETA 坐標有限區(qū)域數(shù)值模式對1975 年8 月河南特大暴雨進行數(shù)值敏感性試驗,結(jié)果表明,暴雨強度對山脈相對高度、陡峭程度和山脈相對位置等非常敏感,地形稍有改變,暴雨強度將減少3/4 以上。馬月枝等(2017)對2016 年7 月9 日新鄉(xiāng)暖區(qū)特大暴雨成因進行綜合分析,結(jié)果表明,太行山地形抬升、中低層低渦和切變線使新鄉(xiāng)上空出現(xiàn)深厚垂直上升運動是暴雨形成的動力機制。研究指出,業(yè)務(wù)數(shù)值模式對“7·21”極端暴雨預報失敗的主要原因可能是預報的低渦和地面倒槽強度和位置有較大誤差(孟智勇等, 2019)。

    2021 年7 月17~22 日河南遭遇極端暴雨。7月20 日00 時至21 日00 時(協(xié)調(diào)世界時,下同),鄭州市內(nèi)24 小時降雨量達到624.1 mm,超過2019 年全年降水量(509.5 mm)。鄭州、新鄉(xiāng)、開封、周口、洛陽等10 個國家級氣象觀測站日雨量突破有氣象記錄以來歷史極值。由于暴雨突發(fā)性、極端性和不確定性,再加上業(yè)務(wù)數(shù)值模式對極端降水的模擬預報能力有限,增加了預報難度(王寧等, 2014; 莊瀟然等, 2017; 孟智勇等, 2019)。此次暴雨發(fā)生在多尺度系統(tǒng)相互作用背景下,受臺風、副熱帶高壓(副高)和西北槽共同影響,高、低空急流耦合,中尺度渦旋長期穩(wěn)定徘徊,在鄭州等地引發(fā)持續(xù)強降水,具有一定的典型性和代表性,值得進行分析總結(jié)。為此本文利用多源數(shù)據(jù),從降水特點、環(huán)境條件和中尺度系統(tǒng)演變等方面對此次河南極端暴雨進行綜合分析,以期厘清動、熱力學特點,歸納科學問題,促進對此次極端降水事件的認識理解,為未來此類極端天氣的分析和預報提供參考。

    2 資料說明與降水實況

    2.1 資料說明

    本文采用的資料包括中國氣象局MICAPS 系統(tǒng)自動站降水觀測數(shù)據(jù)、FY-4A 衛(wèi)星云頂亮溫(TBB)資料和歐洲中期天氣預報中心ERA5 數(shù)據(jù),其中MICAPS 自動站降水數(shù)據(jù)覆蓋全國范圍,時間間隔1 小時。TBB 亮溫數(shù)據(jù)是FY-4A 衛(wèi)星成像儀全圓盤4KML1 數(shù)據(jù),由第12 通道掃描成像輻射計0 級源包數(shù)據(jù)經(jīng)過質(zhì)量檢驗、地理定位、輻射定標處理后得到的產(chǎn)品數(shù)據(jù),水平分辨率為4 km。ERA5 數(shù)據(jù)是歐洲中期天氣預報中心第五代大氣再分析數(shù)據(jù),起始時間為 1979 年,并持續(xù)更新至實時5 天以內(nèi),水平分辨率為0.25°、垂直方向37 層,逐小時間隔(https://confluence.ecmwf.int/display/CKB/ERA5%3A+data+documentation [2020-07-15])(劉鴻波等, 2021)。

    2.2 降水實況

    2021 年7 月17~22 日河南中北部出現(xiàn)暴雨、大暴雨和局地特大暴雨,鄭州、鶴壁和新鄉(xiāng)等地局地累計降雨量超過900 mm,強降水中心位于鄭州,最強降水時段出現(xiàn)在20 日。此次過程最早發(fā)生在2021 年7 月17 日00 時河南東北部,隨后加強西移。如圖1 所示,19 日強降水中心移至鄭州附近,呈現(xiàn)中尺度雨團狀,不斷增強,長時間滯留,20日產(chǎn)生創(chuàng)紀錄日降水量;21 日降水中心位于河南東北部,降水落區(qū)成南北向帶狀,向河北西部伸展;22 日河南降水逐漸減弱,20 時降水結(jié)束。

    圖1 2021 年7 月(a)18 日,(b)19 日,(c)20 日和(d)21 日24 小時累計降水量分布(單位:mm),棕色線為省界,黑點代表鄭州位置Fig. 1 Distribution of 24-hour cumulative precipitation on (a) July 18, (b) July 19, (c) July 20 and (d) July 21, 2021 (units: mm). The brown line represents the provincial boundary, and the black dot represents the location of Zhengzhou

    此次過程屬極端降水事件,具有持續(xù)時間長、小時雨強和累計雨量大、范圍廣、時段集中和極端性突出等特征。河南省境內(nèi)1923 個觀測站降水量大于100 mm,606 個測站大于250 mm;鄭州、嵩山、偃師、新密、伊川和登封等10 個站點突破建站以來日降水量歷史極值。20 日00 時至21 日00時最大降水量出現(xiàn)在鄭州市尖崗,達到696.9 mm,超過鄭州全年平均降水總量;最大小時雨強出現(xiàn)在鄭州本站201.9 mm(20 日08~09 時),突破了中國大陸有氣象記錄以來小時雨強極值。20 日22 時至21 日22 時安陽、鶴壁、新鄉(xiāng)和焦作等地出現(xiàn)局地特大暴雨,其中鶴壁市淇濱區(qū)科創(chuàng)中心最大降水量為675.5 mm,新鄉(xiāng)市牧野區(qū)最大小時雨強達到149.9 mm。

    3 天氣形勢分析

    此次極端暴雨發(fā)生在特殊天氣形勢下,由多尺度天氣系統(tǒng)共同作用造成。西太平洋副高和大陸高壓穩(wěn)定維持在日本海和我國西北地區(qū),受其阻擋,低壓系統(tǒng)在黃淮和華北地區(qū)徘徊。深厚的東南風低空急流和穩(wěn)定的低渦切變線在太行山區(qū)和嵩山等地形強迫作用下引起輻合抬升,降水系統(tǒng)穩(wěn)定少動,在河南中北部造成長時間降水。

    3.1 大尺度環(huán)流

    在200 hPa 等壓面上(圖2a),19 日12 時我國中部逐步建立“兩槽一脊”流型,20 日00 時高壓脊頂部位于河北與河南交界處,東西兩側(cè)為南北向低槽,西低槽大約位于105°E,東低槽大約位于125°E;受高壓脊影響,河南地區(qū)存在西南風向西北風轉(zhuǎn)換的反氣旋環(huán)流,位于高空急流區(qū)的南側(cè),氣流輻散非常顯著。

    在500 hPa 等壓面上(圖2b),20 日00 時短波槽影響我國西部地區(qū)。副高南部的臺風“煙花”大約位于(23°N,132°E)。臺風“查帕卡”在廣東陽江登陸,黃淮氣旋位于陜西中部,河南位于黃淮氣旋的東部。副高和臺風“煙花”之間的東南氣流向河南地區(qū)輸送水汽,與西南氣流匯聚在河南鄭州,為暴雨提供動力和水汽輻合條件。

    在700 hPa 等壓面上(圖2c),20 日00 時黃淮氣旋位于(33°N~36°N,112°E~115°E),鄭州位于其東北側(cè),西南氣流(來自“查帕卡”)和東南氣流(來自副高西側(cè))輻合明顯,鄭州附近東南風速最大為12 m s-1。

    受副高和臺風“煙花”影響,20 日00 時 925 hPa等壓面上從東海至河南鄭州存在強勁的東南氣流(圖2d),把洋面水汽輸送到鄭州;同時“查帕卡”偏南氣流穿越水汽高值區(qū),也向鄭州輸送部分水汽,這兩股輸送水汽的氣流受太行山阻擋,在山前迎風坡(鄭州西部)堆積輻合抬升,促進降水系統(tǒng)發(fā)生發(fā)展。

    圖2 2021 年7 月20 日00 時(協(xié)調(diào)世界時,下同)(a)200 hPa 位勢高度(等值線,單位:m)和風速(填色,單位:m s-1)以及風矢量(風羽,單位:m s-1),(b)500 hPa 位勢高度(等值線;單位:m)和風速(填色,單位:m s-1),(c)700 hPa 位勢高度(等值線,單位:m)和風速(填色,單位:m s-1)以及風矢量(風羽,單位:m s-1)和(d)925 hPa 位勢高度(等值線,單位:m)和水汽比濕(填色,單位:g kg-1)以及風矢量(風羽,單位:m s-1)。紅點代表鄭州位置Fig. 2 (a) 200 hPa geopotential height (isoline, units: m), wind speed (shaded, units: m s-1) and wind vector (wind barb, units: m s-1), (b) 500 hPa geopotential height (isoline, units: m) and wind speed (shaded, units: m s-1), (c) 700 hPa geopotential height (isoline, units: m) and wind speed (shaded,units: m s-1) and wind vector (wind barb, units: m s-1) and (d) 925 hPa geopotential height (isoline, units: m) and water vapor specific humidity(shaded, units: g kg-1) and wind vector (wind barb, units: m s-1) at 0000 UTC 20 July 2021. The red dot represents the location of Zhengzhou

    綜上所述,此次暴雨大尺度環(huán)流特點顯著,高空兩槽一脊流型創(chuàng)造高層強烈輻散的動力條件;中層短波槽、西部大陸高壓、東部副高和臺風“煙花”、南部臺風“查帕卡”是主要的天氣尺度影響系統(tǒng),最重要的是緩慢北移的黃淮氣旋,它是此次暴雨的直接參與者。低層西南氣流與副高和臺風“煙花”之間的東南氣流長時間控制河南地區(qū),供應充沛水汽。這種高、低空系統(tǒng)配置提供有利的動、熱力和水汽條件,促使河南暴雨長時間維持,產(chǎn)生極端降水量。

    3.2 水汽條件

    利用ERA5 再分析數(shù)據(jù)計算了整層大氣可降水量(簡稱:PW),如圖3 所示,可降水量高值區(qū)呈南北向帶狀,從廣東和廣西向北伸展,經(jīng)過湖南、江西和湖北,伸展到河南北部,高值中心位于臺風“帕查卡”和河南中部,臺風倒槽和黃淮氣旋的西南氣流穿過這條高濕帶向河南輸送水汽。另一條相對較弱的可降水量高值區(qū)從東海經(jīng)浙江、江蘇和安徽,伸展到河南,副高和臺風“煙花”引導東南氣流穿越該濕區(qū),向河南輸送水汽。這兩條水汽輸送帶匯合在河南北部,在鄭州附近形成水汽通量輻合中心,為暴雨供應水汽。

    圖3 2021 年7 月20 日(a)00 時和(b)12 時大氣可降水量(填色,單位:mm)與850 hPa 風矢量(箭頭,單位:m s-1),(c)00 時和(d)12 時1000~500 hPa 水汽通量的垂直積分(箭頭,單位:10 kg m-1 s-1)及其散度的垂直積分(填色,單位:10-4 kg m-2 s-1)的水平分布。紅點代表鄭州位置Fig. 3 Atmospheric precipitable water (shaded, units: mm) and 850 hPa wind vector (arrow, units: m s-1) at (a) 0000 UTC and (b) 1200 UTC 20 July 2021. 1000-500 hPa vertical integral of water vapor flux (arrow, units: 10 kg m-1 s-1) and its divergence (shaded, unit: 10-4 kg m-2 s-1) at (c) 0000 UTC and (d) 1200 UTC 20 July 2021. The red dot represents the location of Zhengzhou

    在垂直剖面內(nèi)(圖4),20 日00 時在鄭州西側(cè)迎風坡從近地面到500 hPa 都為水汽通量輻合區(qū),強輻合中心位于950 hPa 附近,主要由邊界層東風急流輸送引起。鄭州近地面東風分量比南風分量更強勁,受西側(cè)地形阻擋,邊界層水汽通量輻合抬升,形成水汽通量輻合強中心。由于水汽濕層深厚,水汽輻合層從近地面伸展到550 hPa,為降水持續(xù)發(fā)展源源不斷供應水汽。從流場結(jié)構(gòu)來看,鄭州東側(cè)低層為東風入流,沿鄭州西側(cè)山地迎風坡爬升,在650 hPa 以下從后部流出,低層垂直運動主要由地形強迫抬升產(chǎn)生;在中高層,垂直運動發(fā)展旺盛,一支上升氣流在300 hPa 向后流出,另一支氣流在400~200 hPa 向前流出。這種流場結(jié)構(gòu)與典型強對流的流場結(jié)構(gòu)類似。

    圖4 2021 年7 月20 日(a)00 時和(b)08 時沿34.43°N 水汽通量散度(填色,單位:10-7 g cm-2 hPa-1 s-1)和風矢量(箭頭,m s-1)的緯向—垂直分布,其中紅三角代表鄭州位置Fig. 4 Zonal-vertical distribution of moisture flux divergence (shaded, units: 10-7 g cm-2 hPa-1 s-1) and wind vector (arrow, units: m s-1) along 34.43°N at (a) 0000 UTC and (b) 0800 UTC on July 20, 2021, where the red triangle represents the location of Zhengzhou

    一般情況下,環(huán)境相對濕度越高、雨滴蒸發(fā)率越小,降水效率越高。20 日08 時強降水開始前(圖5a),鄭州整層相對濕度已經(jīng)飽和,并且部分時段達到過飽和(相對濕度大于100%),在較低溫度下,更容易凝結(jié),暖云發(fā)展更旺盛。20 日08~09 時隨著強降水的發(fā)生,大量水汽被消耗,低層和高層的相對濕度迅速降低,降到90%以下,但在邊界層急流持續(xù)輸送下,水汽得到快速補充,20 日 10 時相對濕度很快達到并超過飽和。如圖5b所示,20 日08 時鄭州上空暖云(雨水+云水)層非常厚,從近地面垂直伸展到400 hPa。研究表明,暖云層越厚,越有利于高降水效率的產(chǎn)生。

    極端暴雨通常具有較高的降水效率。本文采用水汽保留率和水汽消耗率來分析水汽轉(zhuǎn)化效率,這兩個量與降水效率有一定的關(guān)聯(lián)。水汽保留率可表示為當前時刻的可降水量和前一時刻的可降水量與水汽通量散度時間積分之差的比值,即,

    如圖5c 和5d 所示,18 日01 時至19 日06 時,鄭州保持較高的水汽保留率和較低水汽消耗率。在19 日06 時至21 日08 時強降水時段,水汽保留率明顯下降,水汽消耗率上升到較高水平,平均接近8%,個別時刻達到20%;21 日08 時至22 日00 時隨著降水減弱結(jié)束,水汽消耗率也逐漸降到較低水平。水汽消耗率與降水效率有一定關(guān)系。水汽消耗主要來自水汽凝結(jié)和凝華;如果降水效率高,那么云水轉(zhuǎn)化和雪融化引起的雨水生成量會增大,需要水汽凝結(jié)補充云水,消耗大量的水汽,導致較高的水汽消耗率。

    圖5 2021 年7 月19 日00 時至21 日00 時鄭州(34.43°N,113.39°E)(a)相對濕度、(b)液態(tài)水混合比含量(單位:10-4 kg kg-1)的高度—時間演變和2021 年7 月18 日01 時至22 日00 時(c)水汽保留率、(d)水汽消耗率的時間演變Fig. 5 Height-time evolution of (a) relative humidity, (b) mixing ratio content of liquid water (units: 10-4 kg kg-1) in Zhengzhou (34.43°N, 113.39°E)from 0000 UTC 19 July to 0000 UTC 21 July 2021, and time evolution of (c) water vapor reserved rate, (d) water vapor consumption rate from 0100 UTC 18 July to 0000 UTC 22 July 2021

    上述分析表明,西南氣流和東南氣流輸送水汽到河南,在地形阻擋作用下輻合抬升。鄭州低層具有大范圍深厚濕層,可降水量異常偏高,其變化與降水演變有一定關(guān)聯(lián)。流場具有低層東風流入爬坡、高層后部流出和前部流出的典型強對流結(jié)構(gòu)特征。另外,大范圍高濕環(huán)境、深厚暖云層、水汽過飽和以及較高的水汽消耗率是鄭州高降水效率的有利條件。

    3.3 不穩(wěn)定條件分析

    楊博雷等(2016)研究結(jié)果表明,“7·21”暴雨前期對流層中高層高濕位渦的冷空氣擴散南下,冷空氣到達華北地區(qū)上空時,在有利等熵面的引導下從穩(wěn)定層結(jié)向不穩(wěn)定層結(jié)快速下滑,產(chǎn)生了劇烈的正渦度個別變化,使得低渦得到發(fā)展加強。此次暴雨也存在干冷空氣侵入,影響低層大氣的穩(wěn)定性。在等壓坐標系中位渦可寫為

    其中,u和v為緯向和經(jīng)向風速分量, θ為位溫,f0為科氏力參數(shù)。在位渦垂直剖面內(nèi)(圖6a),20日00 時鄭州附近位渦從100 hPa 高層位渦源區(qū)下傳到低層,代表干冷空氣侵入,有利于形成上干冷下暖濕的不穩(wěn)定層結(jié);同時也有利于把高層渦度傳到低層,加強低層渦旋的發(fā)展。從相當位溫垂直分布來看(圖6b),鄭州上空是弱層結(jié)不穩(wěn)定,其南側(cè)低層大氣是強層結(jié)不穩(wěn)定。鄭州925~500 hPa 相當位溫等值線幾乎垂直伸展,表明大氣是近乎中性或者弱不穩(wěn)定層結(jié)。由于飽和濕絕熱過程中氣塊沿著等相當位溫線運動,因此中性層結(jié)非常有利于垂直運動發(fā)展(許煥斌和丁正平, 1997; 李俊等, 2012)。

    圖6 2021 年7 月20 日00 時沿113.39°E(a)位渦(單位:PUV)和(b)相當位溫(等值線,單位:K)、層結(jié)穩(wěn)定度?θse/?p(填色,單位:10-4 K Pa-1)的經(jīng)向—垂直分布。紅三角形代表鄭州位置Fig. 6 Meridional-vertical distribution of (a) potential vorticity (units: PUV), (b) equivalent potential temperature (isoline; units: K) and stratification stability ?θse/?p (shaded, units: 10-4 K Pa-1) along 113.39°E at 0000 UTC 20 July 2021. The red triangle represents the location of Zhengzhou

    在時間演變上(圖7),鄭州850 hPa 和600 hPa高度保持弱層結(jié)不穩(wěn)定,邊界層內(nèi)穩(wěn)定與強層結(jié)不穩(wěn)定交替出現(xiàn)。7 月18 日00 時至22 日00 時鄭州經(jīng)歷四次強層結(jié)不穩(wěn)定與低層弱層結(jié)不穩(wěn)定的轉(zhuǎn)化。在強降水發(fā)生前期,19 日06 時大氣不穩(wěn)定度明顯增強;隨著強降水發(fā)展,19 日16 時不穩(wěn)定能量逐漸被消耗,大氣趨于弱不穩(wěn)定;當弱不穩(wěn)定度減小到一定程度,降水減弱;20 日10 時不穩(wěn)定度又開始增強,為后續(xù)強降水積蓄能量。

    層結(jié)不穩(wěn)定變化可以用層結(jié)穩(wěn)定度傾向方程來解釋(劉璐等, 2015)。該方程可以寫為

    在方程(4)中,作為層結(jié)穩(wěn)定度個別變化的強迫項,位勢散度包括垂直切變項和散度項(周圍等,2018)。劉璐等(2015)研究表明臨近暴雨發(fā)生時刻及暴雨初期,大氣低層主要以對流不穩(wěn)定為主。計算結(jié)果表明(圖7),位勢散度在低層存在五個負高值區(qū),與層結(jié)穩(wěn)定度減弱時段相對應,主要作用是增強層結(jié)穩(wěn)定度;其中散度項的貢獻是削弱不穩(wěn)定層結(jié),促進穩(wěn)定度增強,這是因為地形阻擋抬升使得近地面始終存在氣流輻合((?u/?x+?v/?y)<0),同時低層大氣又是層結(jié)不穩(wěn)定的(?θse/?p>0),因此散度項(?u/?x+?v/?y)?θse/?p<0,促使穩(wěn)定度增強。垂直切變項[-(?u/?p)(?θse/?x)-(?v/?p)(?θse/?y)]的高值區(qū)在低層,正負相間,由水平風垂直切變和大氣的斜壓性共同作用造成。除此之外,層結(jié)穩(wěn)定度的水平平流項[-u(?2θse/?x?p)-v(?2θse/?y?p)]在低層900 hPa 以下促使穩(wěn)定度增強,在900~850 hPa 高度附近增強層結(jié)不穩(wěn)定,而垂直平流項[-ω(?2θse/?2p)]則幾乎與水平平流項反位相相關(guān),二者在對流層中低層大部分相互抵消。

    圖7 2021 年7 月18 日00 時至22 日00 時鄭州(34.43°N,113.39°E)(a)相當位溫(等值線,單位:K)以及層結(jié)穩(wěn)定度(填色,單位:10-4 K Pa-1)、(b)層結(jié)穩(wěn)定度局地變化(單位:10-8 K s-1 Pa-1)、(c)水平平流項(單位:10-8 K s-1 Pa-1)、(d)垂直平流項(單位:10-8 K s-1 Pa-1)、(e)垂直風切變強迫項(單位:10-8 K s-1 Pa-1)和(f)散度強迫項(單位:10-8 K s-1 Pa-1)的高度—時間演變Fig. 7 Height-time evolution from 0000 UTC 18 July to 0000 UTC 22 July 2021, over Zhengzhou (34.43°N, 113.39°E): (a) Equivalent potential temperature (isoline, units: K) and stratification stability (shaded, units: 10-4 K Pa-1); (b) local changes of stratification stability (units: 10-8 K s-1 Pa-1);(c) horizontal advection term (units: 10-8 K s-1 Pa-1); (d) vertical advection term (units: 10-8 K s-1 Pa-1); (e) vertical wind shear forcing term (units:10-8 K s-1 Pa-1); (f) divergence forcing term (units: 10-8 K s-1 Pa-1)

    3.4 動力抬升條件

    20 日00 時至21 日00 時鄭州上空950~350 hPa為氣旋環(huán)流(圖8),近地面900 hPa 以下氣流持續(xù)輻合。20 日00 時垂直運動加強,從地面到200 hPa均為上升運動區(qū),強上升中心位于550 hPa;20 日08 時由于降水粒子下沉拖曳作用,鄭州上空800 hPa以上垂直運動以下沉為主,伴有輻散,10 時轉(zhuǎn)為上升運動。

    圖8 2021 年7 月19 日00 時至21 日00 時鄭州(34.43°N,113.39°E)(a)渦度(單位:單位:10-4 s-1),(b)散度(單位:10-4 s-1)和(c)垂直速度(單位:Pa s-1)的高度—時間演變Fig. 8 Height-time evolution of (a) vorticity (units: 10-4 s-1), (b) divergence (units: 10-4 s-1), and (c) vertical velocity (units: Pa s-1) in Zhengzhou(34.43°N, 113.39°E) from 0000 UTC 19 July to 0000 UTC 21 July 2021

    如圖9 所示,20 日08 時500 hPa 在黃淮氣旋內(nèi)部存在多個中尺度渦旋(圖9a),其中鄭州西側(cè)的中尺度渦旋較強。700~850 hPa 湖北至河南中北部存在東北—西南向的切變線,其北端有中尺度渦旋生成發(fā)展,該渦旋非常深厚,從950 hPa 垂直伸展到350 hPa,并且高值區(qū)位于800~540 hPa。在整個強降水時段,切變線和中尺度渦旋在鄭州地區(qū)穩(wěn)定少動,其西南氣流與副高西側(cè)的東南氣流形成輻合區(qū)。925 hPa 鄭州地區(qū)北側(cè)為氣旋性環(huán)流,風速約為12~16 m s-1,達到邊界層急流標準。同時,東北—西南向的太行山與東西向的嵩山形成喇叭口地形(圖9d),氣旋性急流灌入喇叭口地形,遇到地形阻擋,在鄭州西側(cè)山前形成風速輻合(圖9f);輻合層深厚,從950 hPa 伸展到500 hPa,輻合中心位于迎風坡950 hPa,輻散中心位于350 hPa。這種深厚的中尺度渦旋和輻合層穩(wěn)定存在,有利于垂直運動維持,促使暴雨長時間發(fā)展。

    圖9 2021 年7 月20 日00 時(a)200 hPa、(b)500 hPa、(c)700 hPa、(d)925 hPa 風羽和渦度(填色,單位:10-4 s-1)水平分布以及(e)沿34.43°N 渦度(單位:10-4 s-1)和(f)散度(單位:10-4 s-1)的緯向—垂直分布。箭頭代表風矢量(單位:m s-1) ,紅點和紅三角代表鄭州位置,藍色實線代表地形高度(單位:m)Fig. 9 Horizontal distribution of wind barbs and vorticity (shaded, units: 10-4 s-1) at 0000 UTC July 20 2021: (a) 200 hPa; (b) 500 hPa; (c) 700 hPa;(d) 925 hPa. Zonal-vertical distribution along 34.43°N: (e) vorticity (units: 10-4 s-1); (f) divergence (units: 10-4 s-1). The arrow represents the wind vector (units: m s-1), the red dot and triangle represent the location of Zhengzhou, blue lines indicate the terrain height (units: m)

    K=(u2+v2)/2

    水平動能 [] 分析表明(圖10a),20 日00 時鄭州上空動能高值區(qū)垂直伸展,高、低空動能中心分別位于875 hPa 和175 hPa,分別對應著低空急流和高空急流,二者耦合產(chǎn)生低層輻合和高層輻散的動力抽吸效應,加強上升運動。利用動能傾向方程可分析低空急流的產(chǎn)生機制,在等壓坐標系中,水平動能方程可寫為

    圖10 2021 年7 月20 日00 時沿34.43°N(a)K(等值線,單位:m2 s2)和 ?K / ?t(填色,單位:10-3 m2 s)、(b)-u·?φ/?x(單位:10-3 m2 s)、(c)-v·?φ/?y(單位:m2 s)的緯向—垂直分布以及925 hPa(d) φ(單位:m2 s-2)、(e) u(單位:m s-1)和(f)?φ/?x(單位:m s-2)的水平分布(填色)。藍色線為地形高度(單位:m)等值線,紅點和紅三角代表鄭州位置Fig. 10 Latitude-vertical distribution of (a) K (isoline, units: m2 s2) and ?K/?t (shaded, units: 10-3 m2 s), (b)-u·?φ/?x (units: 10-3 m2 s),(c)-v·?φ/?y (units: 10-3 m2 s) along 34.43°N and horizontal distribution (shaded) of (d) φ(units: m2 s-2), (e) u (units: m s-1), and (f) ?φ/?x (unit:m s-2) 925 hPa at 0000 UTC 20 July 2021. Blue lines indicate the terrain height (units: m), and the red dot and triangle represent the location of Zhengzhou

    鄭州西部為屬于伏牛山系的嵩山(海拔350~1512 米),與河南西北部太行山余脈(海拔高度1000~1600 m)構(gòu)成喇叭口地形。如圖11 所示,20 日00 時近地面東北氣流和偏東氣流灌入喇叭口地形,在中尺度渦旋(中心位于34.75°N)引導下,氣流在嵩山北側(cè)和東側(cè)爬坡,地形阻擋強迫產(chǎn)生輻合上升運動。20 日09 時中尺度渦旋北移至35°N,引導偏東氣流灌入喇叭口,在底部爬升,強勁東風從嵩山的東側(cè)和北側(cè)爬坡,風速達到8~12 m s-1,使得山前大量水汽和動能堆積,激發(fā)和加強對流發(fā)展,促使降水明顯增幅。

    圖11 2021 年7 月20 日(a)00 時和(b)09 時地形追隨坐標系模式面第一層的風矢量(單位:m s-1)。填色代表地形高度(單位:km),紅點代表鄭州位置Fig. 11 Wind vector (units: m s-1) of the first layer on the model plane of the terrain-following coordinates at (a) 0000 UTC and (b) 0900 UTC 20 July 2021. The shaded represents the terrain height (units: km), and the red dot represents the location of Zhengzhou

    綜上所述,降水開始前及其期間,高、低空急流耦合加強抽吸效應,中低層切變線和黃淮氣旋內(nèi)部中尺度渦旋發(fā)展,偏東風低空急流遭遇地形阻擋產(chǎn)生強烈持續(xù)的輻合,這些動力條件促使河南鄭州上空長時間存在大范圍深厚上升運動,地形抬升作用在鄭州暴雨中發(fā)揮重要作用。

    4 中尺度對流系統(tǒng)的演變特征

    此次暴雨有明顯的中尺度對流系統(tǒng)活動,為了揭示中尺度系統(tǒng)特征,利用FY-4A 衛(wèi)星紅外云圖和廣義濕位渦對中尺度對流系統(tǒng)的發(fā)展演變進行分析。

    4.1 云團特征

    在衛(wèi)星云圖上中尺度對流系統(tǒng)定義為水平尺度在20~200 km 之間、生命史3 h 或以上、云頂亮溫(TBB)小于或等于-32°C 的云團(趙思雄等,2007)。FY-4A 衛(wèi)星TBB 分析表明(圖12),河南省中部和東北部始終存在中尺度對流系統(tǒng),并經(jīng)歷三次發(fā)展過程,7 月19 日00~20 時、19 日21時至20 日16 時和20 日17 時至21 日12 時。20日00 時西南氣流和東南氣流匯聚在河南中部和東北部,形成強大的結(jié)構(gòu)密實的中尺度云團,呈現(xiàn)獨立團狀,TBB 中心值降至-56°C,穩(wěn)定少動,降水開始明顯增大。20 日09 時湖北和安徽有小云團生成,在西南氣流引導下,向東北方向移動,不斷地匯入河南中尺度云團,使其不斷加強,中心TBB降到-60°C,長時間停滯在河南中北部,造成20日08~09 時鄭州1 h 降水量達到201.9 mm。16 時云團緩慢南移,同時有較強的小云團從西部移入合并,中心TBB 進一步下降到-64°C。21 日00 時云團恢復到孤立結(jié)構(gòu)密實階段,中心移到河南東南部,鄭州脫離其影響。中尺度云團的三次發(fā)展過程類似,都是在經(jīng)歷長達十幾小時緩慢發(fā)展,合并西南方向移入的小尺度云團后,形成尺度較大的孤立、結(jié)構(gòu)密實的中尺度云團。此次暴雨是中尺度云團三次迭代發(fā)展的結(jié)果,云團合并有利于降水增幅,20 日08~09 時鄭州強降水主要發(fā)生在中尺度和小尺度云團合并期間。

    圖12 2021 年7 月20 日(a)00 時、(b)09 時、(c)16 時和(d)7 月21 日00 時 FY-4A 云頂亮溫(小于-32 °C,單位:°C)分布。紅點代表鄭州位置Fig. 12 Distribution of FY-4A cloud-top brightness temperature (less than -32°C, units: °C) at (a) 0000 UTC, (b) 0900 UTC, and (c) 0000 UTC 20 July 2021 and at (d) 0000 UTC on 21 July 2021. The red dot represents the location of Zhengzhou

    這些中尺度系統(tǒng)在動、熱力學上如何表達刻畫呢?廣義濕位渦是利用廣義位溫定義的位渦,可以綜合描述中尺度系統(tǒng)的動、熱力學特點,其表達式為

    其中,φ(T)為凍結(jié)概率函數(shù),qv為水汽比濕,qs為飽和水汽比濕,k=9為經(jīng)驗指數(shù)(Wang and Huang, 2018)。計算結(jié)果表明(圖13),20 日00 時鄭州附近為廣義位溫高值區(qū),說明存在較強的潛熱釋放,與降水區(qū)相對應。在垂直結(jié)構(gòu)上,鄭州上空廣義位溫高值柱顯著,從300 hPa 垂直下探到近地面,由于包含凝結(jié)和凍結(jié)過程的潛熱釋放,因此廣義位溫能夠更好地表征降水系統(tǒng)高濕和高溫的熱力狀況,與降水落區(qū)和移動關(guān)系密切。

    圖13 2021 年7 月20 日00 時(a)925 hPa 廣義位溫(單位:K)水平分布和(b)沿113.39°E 廣義位溫的經(jīng)向—垂直分布。紅點和紅三角代表鄭州的位置Fig. 13 Generalized potential temperature (units: K) at 0000 UTC 20 July 2021: (a) Horizontal distribution at 925 hPa; (b) meridional-vertical distribution along 113.39°E. The red dot and triangle represent the location of Zhengzhou

    如圖14a 和14b 所示,20 日00 時和09 時廣義濕位渦在降水區(qū)上空表現(xiàn)為異常高值區(qū),從100 hPa 位渦源區(qū)向下延伸到近地面層,在925 hPa形成廣義濕位渦異常中心,主要與邊界層急流垂直風切變、中尺度渦旋和濕斜壓性等密切相關(guān)。高層位渦源與低層廣義濕位渦異常的性質(zhì)不同。由于高層水汽含量很小,廣義位溫退化為位溫,因此100 hPa以上的位渦源主要是位溫垂直梯度和反氣旋環(huán)流的耦合作用造成的。在位渦下傳過程中性質(zhì)發(fā)生改變,由于降水區(qū)上空水汽相變釋放潛熱,廣義位溫變?yōu)橄喈斘粶?,等值線密集,造成降水區(qū)低層廣義濕位渦異常。如圖14c-f 所示,廣義濕位渦的異常區(qū)覆蓋觀測降水落區(qū),二者對應關(guān)系良好,表明廣義濕位渦對降水區(qū)有一定的指示意義,可以為降水趨勢預報提供參考,這主要是因為廣義濕位渦可以表征降水區(qū)上空中尺度系統(tǒng)的垂直風切變、渦度以及大氣濕斜壓性和層結(jié)不穩(wěn)定等動、熱力垂直結(jié)構(gòu)特點。

    圖14 2021 年7 月20 日(a)00 時和(b)09 時沿113.39°E 廣義位渦(單位:PVU)經(jīng)向—垂直分布以及(c)00 時、(b)11 時、(e)13 時和(f)22 時 廣義位渦絕對值的垂直積分(填色,單位:PVU)和觀測1 小時降水量(等值線,單位:mm)的分布。紅點和紅三角代表鄭州的位置Fig. 14 Meridional-vertical distribution of generalized potential vorticity (units: PVU) along 113.39°E on 20 July 2021: (a) 0000 UTC and (b) 0900 UTC. Vertical integration of absolute values of generalized potential vorticity (shaded, units: PVU) and distribution of observed 1-hour precipitation(isoline, units: mm) on 20 July 2021: (c) 0000 UTC; (b) 1100 UTC; (e) 1300 UTC; (f) 2200 UTC. The red dot and triangle represent the location of Zhengzhou

    綜上所述,此次暴雨合并小尺度云團后形成孤立、結(jié)構(gòu)密實的中尺度云團。由于降水區(qū)上空較強的潛熱釋放,廣義位溫高值柱顯著。廣義濕位渦在降水區(qū)上空為異常高值,主要與邊界層急流垂直風切變、中尺度渦旋和濕斜壓性等密切相關(guān)。由于廣義濕位渦可以表征降水區(qū)上空動、熱力垂直結(jié)構(gòu)特點,對降水區(qū)有一定的指示意義,可以為降水趨勢預報提供參考。

    5 結(jié)論

    利用MICAPS 觀測數(shù)據(jù)、ERA5 再分析數(shù)據(jù)和FY-4A 衛(wèi)星云頂亮溫數(shù)據(jù)對2021 年7 月20 日河南極端暴雨的動、熱力和水汽特征進行診斷分析,初步得到以下結(jié)論:

    (1) 大尺度環(huán)流特點顯著,高層兩槽一脊環(huán)流型、中層短波槽、西部大陸高壓、東部西伸北抬的副高和西移的臺風“煙花”以及南部剛登陸的臺風“查帕卡”是主要天氣影響系統(tǒng),中部黃淮氣旋是直接參與者,其西南氣流與副高和“煙花”之間的東南氣流輸送水汽到河南,在地形阻擋作用下輻合抬升。

    (2) 降水區(qū)上空濕層非常深厚、可降水量異常偏高,存在水汽過飽和、暖云層深厚、水汽消耗率高的特點,為高降水效率創(chuàng)造有利條件。邊界層內(nèi)弱層結(jié)不穩(wěn)定與強層結(jié)不穩(wěn)定交替出現(xiàn),位勢散度主要貢獻是增強層結(jié)穩(wěn)定度,低空垂直風切變和氣流輻合輻散對層結(jié)穩(wěn)定度變化有重要影響。

    (3)黃淮氣旋內(nèi)部中尺度渦旋發(fā)展,偏東風低空急流遭遇地形阻擋產(chǎn)生強烈持續(xù)的輻合,這些動力條件促使河南鄭州上空長時間維持大范圍深厚上升運動。

    (4)孤立、結(jié)構(gòu)密實的中尺度云團合并小尺度云團對暴雨發(fā)展有重要影響。降水區(qū)上空廣義濕位渦異常,由于廣義濕位渦可以表征中尺度系統(tǒng)的垂直風切變、渦度以及濕斜壓性和層結(jié)不穩(wěn)定等動、熱力垂直結(jié)構(gòu)特點,因而廣義濕位渦異常對降水區(qū)有一定的指示意義。

    本文著重分析此次暴雨的水汽、不穩(wěn)定和動力條件,但還有一些重要科學問題需要進一步研究,包括:此次暴雨中尺度對流系統(tǒng)長時間維持的機制,多尺度天氣系統(tǒng)相互作用的貢獻,地形和邊界層急流作用以及極端雨強中關(guān)鍵云微物理過程和產(chǎn)生高降水效率的機理等等,這些科學問題需要進一步通過數(shù)值敏感性試驗開展深入研究。

    猜你喜歡
    云團中尺度低層
    巴蜀少年齊上一堂云團課
    少先隊活動(2021年3期)2021-12-04 13:08:26
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運動研究
    海洋通報(2020年5期)2021-01-14 09:26:52
    基于深度學習的中尺度渦檢測技術(shù)及其在聲場中的應用
    2016年7月四川持續(xù)性強降水的中尺度濾波分析
    關(guān)于低層房屋建筑工程造價的要點及控制措施探討探索
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    面陣探測下的污染云團紅外光譜仿真?
    物理學報(2017年11期)2017-08-09 07:34:04
    住八樓以上的人,早亡風險低
    益壽寶典(2017年34期)2017-02-26 08:27:20
    葡萄牙現(xiàn)明亮橘色云團似握著火球的拳頭
    潤·文摘(2016年4期)2016-07-13 04:13:03
    低層高密度住宅設(shè)計探討
    99热6这里只有精品| 午夜福利成人在线免费观看| 久久精品国产99精品国产亚洲性色| 特大巨黑吊av在线直播| 欧美变态另类bdsm刘玥| 久久久久国产网址| 日韩av不卡免费在线播放| 97超视频在线观看视频| 美女内射精品一级片tv| 日本三级黄在线观看| 免费一级毛片在线播放高清视频| 黄色一级大片看看| 国产精品福利在线免费观看| 99在线人妻在线中文字幕| av视频在线观看入口| 好男人在线观看高清免费视频| 99久久人妻综合| 亚洲av中文字字幕乱码综合| 国产免费一级a男人的天堂| 亚洲美女视频黄频| 六月丁香七月| 亚洲图色成人| 亚洲色图av天堂| 国产三级在线视频| 一级毛片久久久久久久久女| 国产黄色视频一区二区在线观看 | 寂寞人妻少妇视频99o| 成人国产麻豆网| 欧美一区二区亚洲| 一级黄片播放器| 亚洲在线观看片| 夫妻性生交免费视频一级片| www日本黄色视频网| a级毛片免费高清观看在线播放| 一区二区三区免费毛片| 99久久成人亚洲精品观看| www.av在线官网国产| 国产精品无大码| 国产精品爽爽va在线观看网站| 国产精品嫩草影院av在线观看| 亚洲久久久久久中文字幕| 国产伦理片在线播放av一区| freevideosex欧美| 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| 日韩成人伦理影院| 波多野结衣一区麻豆| 多毛熟女@视频| 我要看黄色一级片免费的| 久久精品国产鲁丝片午夜精品| 女人久久www免费人成看片| 91aial.com中文字幕在线观看| 色哟哟·www| 春色校园在线视频观看| 国产又爽黄色视频| 国产淫语在线视频| 国产女主播在线喷水免费视频网站| 精品国产乱码久久久久久小说| 看免费av毛片| 久久久久久久国产电影| 波野结衣二区三区在线| 免费在线观看黄色视频的| av在线播放精品| 国产精品久久久久久久久免| 欧美日本中文国产一区发布| 欧美精品一区二区免费开放| 天堂8中文在线网| 麻豆精品久久久久久蜜桃| 国产精品人妻久久久久久| 久久免费观看电影| 亚洲情色 制服丝袜| 天堂俺去俺来也www色官网| 18禁裸乳无遮挡动漫免费视频| 久久青草综合色| 国产成人欧美| 老司机影院毛片| 菩萨蛮人人尽说江南好唐韦庄| 卡戴珊不雅视频在线播放| 久热久热在线精品观看| 国产探花极品一区二区| 成人国产麻豆网| 久久午夜福利片| 成年av动漫网址| 看免费成人av毛片| 26uuu在线亚洲综合色| 天堂8中文在线网| 久久人人爽人人片av| 国产又爽黄色视频| 亚洲av中文av极速乱| 国产亚洲午夜精品一区二区久久| 午夜激情av网站| 午夜视频国产福利| 欧美精品一区二区免费开放| 人人妻人人添人人爽欧美一区卜| 亚洲成人一二三区av| 国产午夜精品一二区理论片| 欧美精品一区二区大全| 制服人妻中文乱码| 97精品久久久久久久久久精品| 在线观看免费视频网站a站| 欧美xxⅹ黑人| 中国国产av一级| 国产日韩欧美视频二区| 人妻一区二区av| 91国产中文字幕| 黑人欧美特级aaaaaa片| 国产精品蜜桃在线观看| 老司机影院毛片| 一本色道久久久久久精品综合| 久久国内精品自在自线图片| 精品国产一区二区三区四区第35| 中文天堂在线官网| 日本av手机在线免费观看| 国产成人欧美| 免费女性裸体啪啪无遮挡网站| av免费在线看不卡| a 毛片基地| 国产不卡av网站在线观看| 最新的欧美精品一区二区| 黄片播放在线免费| 三级国产精品片| 亚洲国产精品国产精品| 午夜免费观看性视频| 亚洲久久久国产精品| 国产熟女欧美一区二区| 男女啪啪激烈高潮av片| 菩萨蛮人人尽说江南好唐韦庄| 自拍欧美九色日韩亚洲蝌蚪91| 女人精品久久久久毛片| 亚洲一级一片aⅴ在线观看| 在线观看美女被高潮喷水网站| 卡戴珊不雅视频在线播放| 少妇熟女欧美另类| 日本91视频免费播放| 99国产精品免费福利视频| 国产有黄有色有爽视频| 两性夫妻黄色片 | 日本-黄色视频高清免费观看| 免费高清在线观看视频在线观看| 欧美日韩视频高清一区二区三区二| 性色av一级| 不卡视频在线观看欧美| 菩萨蛮人人尽说江南好唐韦庄| 精品国产国语对白av| 久久精品国产鲁丝片午夜精品| 你懂的网址亚洲精品在线观看| 日韩av免费高清视频| 日韩 亚洲 欧美在线| 免费在线观看完整版高清| 少妇被粗大猛烈的视频| 国产老妇伦熟女老妇高清| 少妇被粗大的猛进出69影院 | 国产精品久久久久久av不卡| 乱码一卡2卡4卡精品| 成年美女黄网站色视频大全免费| 一边亲一边摸免费视频| 亚洲精品456在线播放app| 久久精品国产亚洲av涩爱| 亚洲色图 男人天堂 中文字幕 | 免费人成在线观看视频色| 永久网站在线| 校园人妻丝袜中文字幕| 久久久国产欧美日韩av| 欧美日韩成人在线一区二区| 久久久久久久精品精品| 久久久久久久亚洲中文字幕| 国产极品天堂在线| av电影中文网址| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 久久人人爽av亚洲精品天堂| 国产黄色视频一区二区在线观看| 丝袜脚勾引网站| 看非洲黑人一级黄片| 亚洲av中文av极速乱| 日韩中文字幕视频在线看片| 一区二区三区四区激情视频| 最后的刺客免费高清国语| 国产成人欧美| 午夜免费男女啪啪视频观看| 男女啪啪激烈高潮av片| 日本vs欧美在线观看视频| 久久ye,这里只有精品| 97人妻天天添夜夜摸| 人妻一区二区av| 国产国语露脸激情在线看| 一本久久精品| 午夜福利,免费看| 极品人妻少妇av视频| 精品一区二区三区视频在线| 亚洲国产精品999| 日韩伦理黄色片| 欧美国产精品va在线观看不卡| 午夜免费男女啪啪视频观看| 久久这里只有精品19| 欧美 亚洲 国产 日韩一| 欧美变态另类bdsm刘玥| 精品国产乱码久久久久久小说| 国产精品熟女久久久久浪| 18禁国产床啪视频网站| 亚洲精品乱久久久久久| 天堂俺去俺来也www色官网| 国产免费现黄频在线看| 国产精品女同一区二区软件| av在线app专区| 欧美精品一区二区大全| 久久青草综合色| 日本欧美视频一区| 国产精品一国产av| 99九九在线精品视频| 日韩视频在线欧美| 中文字幕人妻熟女乱码| 最近的中文字幕免费完整| 黄色一级大片看看| 91久久精品国产一区二区三区| 欧美成人精品欧美一级黄| 大陆偷拍与自拍| 国产视频首页在线观看| tube8黄色片| 国产极品天堂在线| 中文乱码字字幕精品一区二区三区| 深夜精品福利| 爱豆传媒免费全集在线观看| 国产精品蜜桃在线观看| 丝袜人妻中文字幕| 中文字幕精品免费在线观看视频 | 一区二区三区四区激情视频| 久热这里只有精品99| 国产国语露脸激情在线看| 国产精品一区www在线观看| 欧美激情极品国产一区二区三区 | 国产精品不卡视频一区二区| 精品一区二区免费观看| 国产亚洲一区二区精品| 国产免费一级a男人的天堂| 亚洲精品av麻豆狂野| 一区二区日韩欧美中文字幕 | 久久99热6这里只有精品| 夫妻性生交免费视频一级片| 久久人人爽av亚洲精品天堂| 成年动漫av网址| 又大又黄又爽视频免费| 在线观看www视频免费| 肉色欧美久久久久久久蜜桃| 日本av手机在线免费观看| 日韩电影二区| 国产麻豆69| 麻豆精品久久久久久蜜桃| 日本免费在线观看一区| 18禁在线无遮挡免费观看视频| 97超碰精品成人国产| 亚洲成国产人片在线观看| 国产精品蜜桃在线观看| 中文欧美无线码| 一区二区三区乱码不卡18| 国产成人一区二区在线| 大片电影免费在线观看免费| 内地一区二区视频在线| 久久99蜜桃精品久久| 欧美人与性动交α欧美软件 | 又黄又粗又硬又大视频| 国产精品.久久久| 亚洲精品国产色婷婷电影| 91成人精品电影| 观看av在线不卡| 曰老女人黄片| 99视频精品全部免费 在线| 在线观看美女被高潮喷水网站| 国产xxxxx性猛交| 成人国产麻豆网| 青春草国产在线视频| 永久网站在线| 又黄又粗又硬又大视频| 欧美日韩综合久久久久久| 青春草亚洲视频在线观看| 国产精品麻豆人妻色哟哟久久| 国产精品一国产av| 韩国高清视频一区二区三区| av网站免费在线观看视频| 亚洲婷婷狠狠爱综合网| 国产成人精品福利久久| 观看美女的网站| 成人午夜精彩视频在线观看| 日韩大片免费观看网站| av黄色大香蕉| 欧美人与性动交α欧美精品济南到 | 伦理电影大哥的女人| 国产精品免费大片| 男人添女人高潮全过程视频| 日日撸夜夜添| 中文字幕制服av| 丝袜人妻中文字幕| 亚洲高清免费不卡视频| 在线观看一区二区三区激情| 日本黄大片高清| 国产精品国产三级国产av玫瑰| 99久久综合免费| 久久久久久久亚洲中文字幕| 国产永久视频网站| 久久国产亚洲av麻豆专区| 亚洲国产成人一精品久久久| 成人亚洲欧美一区二区av| 一本—道久久a久久精品蜜桃钙片| 曰老女人黄片| 尾随美女入室| 搡女人真爽免费视频火全软件| 在线 av 中文字幕| a级片在线免费高清观看视频| 日本av免费视频播放| 嫩草影院入口| 久久婷婷青草| 免费人妻精品一区二区三区视频| 久久 成人 亚洲| 夜夜爽夜夜爽视频| 婷婷色麻豆天堂久久| 亚洲成人手机| 成人二区视频| 国产熟女欧美一区二区| 久久久精品94久久精品| 国产精品蜜桃在线观看| 亚洲欧美精品自产自拍| 久久热在线av| 国产片内射在线| 99热6这里只有精品| 永久网站在线| 亚洲av免费高清在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲精品久久午夜乱码| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品秋霞免费鲁丝片| 久热这里只有精品99| 久久精品国产自在天天线| 妹子高潮喷水视频| 如日韩欧美国产精品一区二区三区| 亚洲精品第二区| 22中文网久久字幕| 亚洲四区av| 熟女电影av网| 欧美人与性动交α欧美软件 | 99久国产av精品国产电影| 黑丝袜美女国产一区| 国产免费视频播放在线视频| 女性生殖器流出的白浆| 美女福利国产在线| 女性生殖器流出的白浆| 99久久综合免费| 五月开心婷婷网| 久久97久久精品| 捣出白浆h1v1| av国产精品久久久久影院| 亚洲性久久影院| 国产日韩一区二区三区精品不卡| 久热这里只有精品99| 一边亲一边摸免费视频| 日韩欧美一区视频在线观看| 色哟哟·www| 日韩电影二区| 欧美少妇被猛烈插入视频| 久久这里有精品视频免费| 成人国语在线视频| 夜夜爽夜夜爽视频| 男女边吃奶边做爰视频| 精品一品国产午夜福利视频| 欧美少妇被猛烈插入视频| 久久这里有精品视频免费| 九色成人免费人妻av| 男女无遮挡免费网站观看| 在线观看三级黄色| h视频一区二区三区| 国产男女内射视频| 久久这里有精品视频免费| 国产男女内射视频| 久久久欧美国产精品| 伦理电影免费视频| 蜜桃在线观看..| www.熟女人妻精品国产 | 亚洲 欧美一区二区三区| 亚洲在久久综合| 亚洲精品久久午夜乱码| 精品酒店卫生间| 亚洲精品国产色婷婷电影| 久久毛片免费看一区二区三区| 大码成人一级视频| 精品一区二区三卡| www.熟女人妻精品国产 | 亚洲欧美日韩卡通动漫| 国产探花极品一区二区| 在线观看免费高清a一片| 免费大片黄手机在线观看| 久久97久久精品| 精品酒店卫生间| 卡戴珊不雅视频在线播放| 国产精品.久久久| 国产精品久久久久久av不卡| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 亚洲成人av在线免费| 亚洲欧美日韩卡通动漫| 成人亚洲欧美一区二区av| 午夜91福利影院| www.av在线官网国产| 国产1区2区3区精品| 国产一区二区激情短视频 | 天堂俺去俺来也www色官网| 大香蕉久久成人网| 久久精品久久久久久久性| 成人黄色视频免费在线看| 欧美日韩亚洲高清精品| 日韩精品有码人妻一区| 高清毛片免费看| 两性夫妻黄色片 | 精品一区二区免费观看| 国产成人av激情在线播放| 成年av动漫网址| 一级片免费观看大全| 午夜免费鲁丝| 免费人成在线观看视频色| 日产精品乱码卡一卡2卡三| av电影中文网址| 免费看光身美女| 精品一区在线观看国产| tube8黄色片| 亚洲国产av新网站| 亚洲精品中文字幕在线视频| 美女大奶头黄色视频| 亚洲婷婷狠狠爱综合网| 日本黄色日本黄色录像| 国产片内射在线| 日韩视频在线欧美| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色av中文字幕| 男女啪啪激烈高潮av片| 中文字幕精品免费在线观看视频 | 国产精品女同一区二区软件| 男女无遮挡免费网站观看| 亚洲综合色惰| 久久久久久久大尺度免费视频| 亚洲精品乱码久久久久久按摩| 麻豆精品久久久久久蜜桃| 日本色播在线视频| 国产乱人偷精品视频| 精品午夜福利在线看| 考比视频在线观看| 日本爱情动作片www.在线观看| 亚洲成色77777| 大香蕉97超碰在线| 婷婷色麻豆天堂久久| 亚洲精品456在线播放app| 草草在线视频免费看| 久久久久久久亚洲中文字幕| 国产亚洲最大av| 91精品三级在线观看| 嫩草影院入口| 国产女主播在线喷水免费视频网站| 一边摸一边做爽爽视频免费| 最新中文字幕久久久久| 国产精品久久久av美女十八| av一本久久久久| 捣出白浆h1v1| 又黄又粗又硬又大视频| 这个男人来自地球电影免费观看 | 欧美国产精品va在线观看不卡| 黄色 视频免费看| 亚洲熟女精品中文字幕| 伊人久久国产一区二区| 国产精品三级大全| 亚洲精品美女久久av网站| 日韩欧美精品免费久久| 久久久亚洲精品成人影院| 亚洲第一av免费看| 搡老乐熟女国产| av有码第一页| 大香蕉97超碰在线| 大香蕉久久成人网| 精品亚洲乱码少妇综合久久| 精品久久国产蜜桃| 精品少妇黑人巨大在线播放| 一级毛片我不卡| 亚洲av欧美aⅴ国产| 免费看av在线观看网站| 亚洲av男天堂| 久热这里只有精品99| 午夜免费鲁丝| av.在线天堂| 亚洲五月色婷婷综合| 欧美bdsm另类| 国产又爽黄色视频| 综合色丁香网| 欧美日韩一区二区视频在线观看视频在线| 少妇 在线观看| 免费在线观看黄色视频的| 在线观看一区二区三区激情| 成人18禁高潮啪啪吃奶动态图| 国产精品.久久久| 老司机影院成人| 黄片无遮挡物在线观看| 乱码一卡2卡4卡精品| 男女高潮啪啪啪动态图| 精品酒店卫生间| 国产精品嫩草影院av在线观看| 亚洲激情五月婷婷啪啪| 九九爱精品视频在线观看| 国产亚洲欧美精品永久| 欧美+日韩+精品| 精品久久久精品久久久| 丝袜人妻中文字幕| 伊人久久国产一区二区| av一本久久久久| 成人免费观看视频高清| 久久久久精品性色| 欧美bdsm另类| 精品一品国产午夜福利视频| 中文字幕av电影在线播放| 亚洲国产精品一区二区三区在线| 亚洲精品中文字幕在线视频| 爱豆传媒免费全集在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久韩国三级中文字幕| 高清黄色对白视频在线免费看| 亚洲成人一二三区av| 日韩大片免费观看网站| 亚洲,欧美精品.| 午夜老司机福利剧场| 亚洲av中文av极速乱| 高清毛片免费看| 国产精品欧美亚洲77777| 亚洲成av片中文字幕在线观看 | 国产精品99久久99久久久不卡 | 男女啪啪激烈高潮av片| 亚洲美女搞黄在线观看| 黑人猛操日本美女一级片| 性高湖久久久久久久久免费观看| 国产 精品1| 免费在线观看黄色视频的| 国产淫语在线视频| 国产成人免费无遮挡视频| 交换朋友夫妻互换小说| 男人舔女人的私密视频| 免费av中文字幕在线| 久久久国产精品麻豆| 国产精品99久久99久久久不卡 | 亚洲av国产av综合av卡| 乱人伦中国视频| 欧美精品av麻豆av| 国产成人精品在线电影| 少妇熟女欧美另类| 我的女老师完整版在线观看| 国产毛片在线视频| 又黄又粗又硬又大视频| 中文欧美无线码| 亚洲欧美色中文字幕在线| 人人妻人人添人人爽欧美一区卜| 亚洲av日韩在线播放| 欧美日韩亚洲高清精品| 国产亚洲午夜精品一区二区久久| videosex国产| 少妇 在线观看| 色婷婷av一区二区三区视频| kizo精华| 午夜老司机福利剧场| 亚洲av中文av极速乱| 蜜桃国产av成人99| 亚洲欧洲国产日韩| av在线app专区| 18禁裸乳无遮挡动漫免费视频| 青青草视频在线视频观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品少妇黑人巨大在线播放| 91精品国产国语对白视频| 在线天堂最新版资源| 熟女人妻精品中文字幕| 国产精品偷伦视频观看了| 亚洲欧美一区二区三区国产| 久久99精品国语久久久| 中文字幕av电影在线播放| 久久久久久久久久人人人人人人| av不卡在线播放| 国产精品久久久久久久电影| 中文乱码字字幕精品一区二区三区| 男女啪啪激烈高潮av片| 国产日韩欧美视频二区| 1024视频免费在线观看| 欧美3d第一页| 日本wwww免费看| 国产xxxxx性猛交| 成人影院久久| 精品少妇黑人巨大在线播放| 大片免费播放器 马上看| 国产成人精品久久久久久| 久久精品夜色国产| 欧美日韩视频精品一区| 在线亚洲精品国产二区图片欧美| 亚洲av.av天堂| 在线观看免费视频网站a站| 黑人欧美特级aaaaaa片| 久久久国产精品麻豆| 老司机影院成人| 欧美日韩成人在线一区二区| 国产爽快片一区二区三区| 国产精品人妻久久久影院| 国产又色又爽无遮挡免| 黄色怎么调成土黄色| 久久久久久久久久人人人人人人| 成年人午夜在线观看视频| 大香蕉久久成人网| 日韩制服丝袜自拍偷拍| 又粗又硬又长又爽又黄的视频| a级片在线免费高清观看视频| 啦啦啦视频在线资源免费观看| 大陆偷拍与自拍| 亚洲精品日本国产第一区| 国产国语露脸激情在线看| 亚洲欧美成人精品一区二区| 亚洲av免费高清在线观看| 日韩不卡一区二区三区视频在线|