亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        低空急流與賀蘭山東麓暴雨過程的相關(guān)性研究*

        2023-11-18 12:06:00陳豫英張一星
        氣象 2023年10期
        關(guān)鍵詞:雨強(qiáng)賀蘭山急流

        蘇 洋 陳豫英 楊 侃 張一星 楊 銀

        1 中國氣象局旱區(qū)特色農(nóng)業(yè)氣象災(zāi)害監(jiān)測預(yù)警與風(fēng)險(xiǎn)管理重點(diǎn)實(shí)驗(yàn)室,銀川 750002 2 寧夏氣象防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室,銀川 750002 3 寧夏氣象臺(tái),銀川 750002 4 寧夏石嘴山市氣象局,石嘴山 753000

        提 要:利用2006—2021年逐時(shí)降水、常規(guī)氣象探測、銀川CA雷達(dá)和ERA5高分辨率再分析資料,研究了低空急流與賀蘭山東麓暴雨過程的時(shí)間和空間的相關(guān)性,并初步探討了低空急流影響暴雨過程發(fā)生發(fā)展的可能機(jī)制。結(jié)果表明:影響賀蘭山東麓暴雨過程的低空急流有三個(gè)關(guān)鍵區(qū),分別為河套南部、寧夏東南部和山西西南部,對應(yīng)700 hPa南風(fēng)急流、775 hPa偏南急流和850 hPa東南急流;寧夏東南部作為三支低空急流匯合后繼續(xù)北上西進(jìn)的關(guān)鍵區(qū)域,對賀蘭山東麓暴雨過程的發(fā)生發(fā)展有著更為重要的影響。依據(jù)低空急流核所在高度,將影響賀蘭山東麓暴雨過程的低空急流分為七類,其中三層急流型出現(xiàn)頻率最高,占總過程的54.5%,其次是700和775 hPa急流同時(shí)出現(xiàn)的過程,占36.5%。暴雨過程與低空急流在時(shí)間上存在一致性:700、775、850 hPa急流建立較暴雨開始平均提前了18、10、7 h,700、775 hPa急流最大風(fēng)速較暴雨過程最大雨強(qiáng)平均提前了54、18 min,而850 hPa急流最大風(fēng)速較暴雨最大雨強(qiáng)平均滯后了12 min;850 hPa的1級(jí)急流與775 hPa的2級(jí)急流頻率分別對20~40、40~60 mm·h-1的短時(shí)暴雨頻率指示性更強(qiáng),而河套南部關(guān)鍵區(qū)的700 hPa平均風(fēng)速對暴雨過程的最大雨強(qiáng)量級(jí)指示性更強(qiáng)。暴雨過程與低空急流在空間上也存在一致性:隨著低空急流建立、加強(qiáng)北抬或西進(jìn)、減弱東退或南壓,賀蘭山東麓暴雨開始、增強(qiáng)、減弱;暴雨落區(qū)位于急流軸的左前方。低空急流北上西進(jìn)與賀蘭山地形結(jié)合,在東坡山前觸發(fā)多個(gè)對流單體、合并加強(qiáng)形成移動(dòng)緩慢、發(fā)展強(qiáng)盛、組織化程度高、列車效應(yīng)明顯的帶狀線性回波,易在賀蘭山區(qū)形成局地性強(qiáng)對流暴雨。

        引 言

        低空急流與暴雨的關(guān)系非常密切,因而得到了廣泛關(guān)注(陶詩言,1980;Higgins et al,1997;Chen et al,2005;Cook and Vizy,2010;鄭婧等,2020),相關(guān)研究可以有效提高暴雨預(yù)報(bào)的準(zhǔn)確率(Huang and Luo,2017)。低空急流對暴雨的影響主要體現(xiàn)在四個(gè)方面:(1)暴雨的日變化特征與低空急流的發(fā)展密切相關(guān)(王東阡和張耀存,2012;徐娟和陳勇明,2013;周靜等,2017),我國東部暖區(qū)暴雨和梅雨鋒暴雨過程在清晨出現(xiàn)的降水峰值與夜間西南低空急流加強(qiáng)有關(guān)(Chen et al,2010;Liu,2012;Jiang et al,2017);(2)低空急流為暴雨區(qū)輸送水汽和不穩(wěn)定能量(Tollerud et al,2008;顧清源等,2009;吳海英等,2010;許彬等,2019;陳芳麗等,2021;蘇愛芳等,2022),其左前方慣性穩(wěn)定度大值區(qū)的存在有利于低層濕熱能量的積累(張文龍和周軍,2003;付煒等,2020);(3)低空急流具有的強(qiáng)烈的不穩(wěn)定性使得急流軸上的風(fēng)速出現(xiàn)強(qiáng)風(fēng)速脈動(dòng)傳播現(xiàn)象,從而觸發(fā)中尺度低渦及暴雨雨團(tuán)發(fā)生(孫淑清和翟國慶,1980;張京英等,2005;趙嫻婷等,2020;胡雅君等,2020),低空急流在地形抬升的作用下造成強(qiáng)烈的垂直上升運(yùn)動(dòng),激發(fā)對流云團(tuán)強(qiáng)烈發(fā)展,也能夠觸發(fā)暴雨過程(Zhang and Meng,2019;黃小彥等,2020;汪小康等,2022;黃美金等,2022);(4)低空急流與暴雨是相互促進(jìn)的,低空急流的存在有利于暴雨的發(fā)生,而暴雨的發(fā)生又促進(jìn)了急流的形成和維持(賽瀚和苗峻峰,2012;劉鴻波等,2014)。孫繼松(2005)、Saulo et al(2007)分別證明了邊界層低空急流和局地暴雨、低空急流和中尺度對流及高空急流存在明顯的正反饋?zhàn)饔谩?/p>

        地處西北地區(qū)東部的賀蘭山東麓地區(qū)易發(fā)生局地性強(qiáng)、突發(fā)性強(qiáng)的暴雨過程,部分學(xué)者圍繞賀蘭山東麓地區(qū)暴雨的低空急流活動(dòng)也開展了一些研究。陳豫英等(2018;2021)、Chen et al(2021)對比研究了低空急流對不同類型賀蘭山東麓極端暴雨的影響,進(jìn)一步指出賀蘭山東麓暴雨分布與低空偏東南急流夜間增強(qiáng)并配合地形作用在賀蘭山東坡山前觸發(fā)或增強(qiáng)暴雨中小尺度系統(tǒng)造成地形處降水增幅密切相關(guān)。楊侃等(2020)研究也表明賀蘭山東側(cè)迎風(fēng)坡阻擋低空急流,強(qiáng)迫暖濕空氣抬升,促使低渦切變強(qiáng)烈發(fā)展,從而影響賀蘭山地區(qū)降水強(qiáng)度、落區(qū)及中心位置。

        雖然賀蘭山東麓暴雨過程相關(guān)的低空急流特征研究取得了一定的成果,但研究僅限于個(gè)別極端暴雨過程,尚未開展系統(tǒng)性研究,影響賀蘭山東麓暴雨過程的低空急流活動(dòng)是否存在關(guān)鍵區(qū)?賀蘭山東麓暴雨過程與低空急流的日變化特征是否存在一致性?低空急流如何影響賀蘭山東麓暴雨過程的降水時(shí)段、落區(qū)和強(qiáng)度?討論這些問題不僅對了解西北干旱地區(qū)復(fù)雜地形條件下的低空急流特征有一定意義,也有助于更準(zhǔn)確地預(yù)報(bào)暴雨。本文利用2006—2021年的逐時(shí)降水、常規(guī)氣象探測和ERA5高分辨率再分析資料,研究低空急流與賀蘭山東麓暴雨過程的時(shí)間和空間的相關(guān)性,以期為賀蘭山東麓暴雨預(yù)報(bào)提供參考。

        1 資料和方法

        1.1 資 料

        2006—2021年賀蘭山東麓地區(qū)(37.8°~39.4°N、105.7°~107.0°E)512個(gè)站逐時(shí)降水、常規(guī)氣象探測、銀川CA雷達(dá)、ERA5高分辨率再分析資料。其中,ERA5資料時(shí)間分辨率為1 h、水平分辨率為0.25°×0.25°、垂直分辨率為23 層(對流層低層垂直分辨率加密,1000~750 hPa為25 hPa/層,750~600 hPa為50 hPa/層,600~100 hPa為100 hPa/層)。文中時(shí)間均為北京時(shí)。

        1.2 暴雨過程標(biāo)準(zhǔn)

        參照Chen et al(2021)研究,定義賀蘭山東麓暴雨過程標(biāo)準(zhǔn):至少一個(gè)測站任意24 h累計(jì)雨量≥100 mm為一次暴雨過程,否則至少相鄰5個(gè)站≥25 mm,其中至少相鄰3個(gè)站≥50 mm;短時(shí)暴雨:1 h 雨量≥20 mm。

        1.3 低空急流標(biāo)準(zhǔn)

        由于低空急流出現(xiàn)的高度、范圍、風(fēng)速及垂直風(fēng)切變均有一定的差異,迄今為止,對低空急流的定義尚未形成統(tǒng)一的標(biāo)準(zhǔn)(劉鴻波等,2014)。參照Zhang and Meng(2019)對華南暖區(qū)暴雨過程的低空急流的選取標(biāo)準(zhǔn),本文對影響賀蘭山東麓暴雨過程的低空急流選取標(biāo)準(zhǔn)定義為:在30°~40°N、100°~120°E范圍,600 hPa以下存在風(fēng)速≥10 m·s-1的氣流帶(長度超過100 km),低層風(fēng)速最大值超出該風(fēng)速出現(xiàn)高度上方風(fēng)速最小值至少2 m·s-1,且暴雨區(qū)與低空急流軸的距離不超過4個(gè)經(jīng)距或緯距。

        考慮到銀川平原地區(qū)海拔高度在1102~1122 m(大約是875 hPa),選取ERA5資料的875、850、825、800、775、750、700、650和600 hPa共9個(gè)高度層進(jìn)行低空急流的統(tǒng)計(jì)分析。定義急流最大頻率:賀蘭山東麓暴雨過程期間,某個(gè)高度層、某個(gè)格點(diǎn)上,急流出現(xiàn)的最多次數(shù)。

        依據(jù)上述標(biāo)準(zhǔn),選取了2006—2021年賀蘭山東麓23次暴雨過程,其中,有22次過程伴隨低空急流出現(xiàn)(表1),占總暴雨過程的95.7%。

        1.4 對比ERA5與探空數(shù)據(jù)低層風(fēng)場

        參考楊程等(2018)對浙江探空數(shù)據(jù)和ERA-Interim 數(shù)據(jù)評(píng)估的方法,選取西北地區(qū)東部的6個(gè)探空站(銀川、延安、崆峒、武都、漢中和安康;圖1),對比探空和ERA5的低層風(fēng)場(700、850 hPa的風(fēng)向和風(fēng)速)。

        注:紅色方框?yàn)橘R蘭山東麓地區(qū)。圖1 西北地區(qū)東部的探空站(紅色圓點(diǎn))分布和地形高度(填色)Fig.1 Distribution of sounding stations (red dot) and terrain height (colored) in the eastern of Northwest China

        從表2中可以看出,兩類數(shù)據(jù)的相關(guān)性高,相關(guān)系數(shù)普遍在0.92以上(通過了α=0.01的顯著性水平檢驗(yàn))。風(fēng)向的絕對誤差在7.31°~10.94°,均方根誤差在9.66°~14.17°;風(fēng)速的絕對誤差在0.83~1.07 m·s-1,均方根誤差在1.26~1.49 m·s-1。從平均絕對偏差來看,銀川站、崆峒站和漢中站的風(fēng)向?yàn)樨?fù)值,即相對于探空風(fēng)向,ERA5的風(fēng)向平均偏右,相反地,對于延安站、武都站和安康站而言,ERA5的風(fēng)向平均偏左;相對于探空風(fēng)速,ERA5的風(fēng)速普遍略偏大。風(fēng)向的相對誤差均低于8%,除武都站外,其他5個(gè)站的風(fēng)速相對誤差均低于16%,這可能與武都站地形有關(guān)。因此,ERA5低層風(fēng)場數(shù)據(jù)總體上是可信的,可用于研究暴雨過程與低空急流的相關(guān)性。

        表2 22次賀蘭山東麓暴雨過程期間37個(gè)時(shí)次的ERA5再分析數(shù)據(jù)與西北地區(qū)東部的6個(gè)探空站數(shù)據(jù)的評(píng)估結(jié)果Table 2 Evaluation result of the ERA5 data and the observed data at 37 times by 6 sounding stations in the eastern of Northwest China during the 22 rainstorm processes in the eastern foot of Helan Mountains

        2 影響賀蘭山東麓暴雨過程的低空急流活動(dòng)關(guān)鍵區(qū)和關(guān)鍵高度

        圖2a為22次賀蘭山東麓暴雨過程的875~600 hPa急流最大頻率分布;因各高度急流發(fā)生次數(shù)差異較大,選取急流次數(shù)最多的高度層(700 hPa,793次)為標(biāo)準(zhǔn),統(tǒng)計(jì)這793個(gè)時(shí)次內(nèi)每個(gè)格點(diǎn)的急流最大頻率,即研究同一個(gè)較長時(shí)間段內(nèi)的急流最大頻率分布。圖2b為對應(yīng)圖2a中急流最大頻率出現(xiàn)的高度分布。由圖可見,影響賀蘭山東麓暴雨過程的低空急流分布大值區(qū)有三個(gè):A區(qū)位于河套南部、B區(qū)位于寧夏東南部、C區(qū)位于山西西南部,其中,A區(qū)范圍最大、頻率最高(中心達(dá)到793次),B區(qū)次之(中心為627次),C區(qū)范圍最小、頻率也最低(中心為510次);A、B和C區(qū)對應(yīng)的平均高度分別為700、775和850 hPa。從各區(qū)最大頻率中心水平風(fēng)速的垂直廓線(圖2c)上也可以看出:A區(qū)低空急流出現(xiàn)的高度最高且強(qiáng)度最強(qiáng),高度為700 hPa,平均風(fēng)速達(dá)到12.9 m·s-1;B區(qū)低空急流出現(xiàn)在775 hPa,急流強(qiáng)度較強(qiáng)(11.0 m·s-1);C區(qū)低空急流高度為850 hPa,強(qiáng)度最弱,最大風(fēng)速為9.8 m·s-1。

        注:圖a,圖b中黑色方框表示急流關(guān)鍵區(qū),A、B、C分別表示三個(gè)關(guān)鍵區(qū)的最大頻率中心,+表示最大頻率中心位置。圖2 22次賀蘭山東麓暴雨過程875~600 hPa的(a)低空急流最大頻率及其(b)對應(yīng)高度的分布,(c)關(guān)鍵區(qū)最大頻率中心水平風(fēng)速的垂直廓線Fig.2 Distribution of (a) maximum frequency and (b) the height corresponding to maximum frequency of low-level jet at 875-600 hPa during the 22 rainstorm processes in the eastern foot of Helan Mountains, and (c) horizontal wind speed vertical profiles of maximum frequency center in the key areas

        從圖3可以看出,700 hPa急流為南風(fēng)急流,775 hPa為偏南急流,850 hPa為東南急流,且急流核與頻率分布大值區(qū)(圖2a)有很好的對應(yīng)關(guān)系,分別位于A、B和C三個(gè)關(guān)鍵區(qū)。775 hPa急流區(qū)距離賀蘭山東麓地區(qū)最近,對該地區(qū)暴雨有直接影響,而700 hPa急流軸最前端恰好處于B區(qū)775 hPa急流核所在位置,850 hPa急流出口區(qū)前方恰好對應(yīng)775 hPa急流入口區(qū);結(jié)合圖1的地形分布可知,850 hPa急流經(jīng)地形抬升可造成B區(qū)775 hPa急流加強(qiáng)。由此可見:B區(qū)是三支低空急流匯合后繼續(xù)北上西進(jìn)的中轉(zhuǎn)站,對賀蘭山東麓暴雨過程的發(fā)生發(fā)展有著極其重要的作用。

        注:藍(lán)色圓點(diǎn)表示合成平均風(fēng)速通過了α=0.01的顯著性水平檢驗(yàn),紅色方框?yàn)橘R蘭山東麓地區(qū)。圖3 22次賀蘭山東麓暴雨過程合成平均的(a)700 hPa,(b)775 hPa和(c)850 hPa低空急流風(fēng)場(填色和風(fēng)矢)Fig.3 Composite average wind field of low-level jet at (a) 700 hPa, (b) 775 hPa and (c) 850 hPa during the 22 rainstorm processes in the eastern foot of Helan Mountains (colored and wind vector)

        3 低空急流分類

        參考Zhang and Meng(2019)、賽瀚和苗峻峰(2012)研究,結(jié)合寧夏地形,考慮B區(qū)的海拔高度為1400~1600 m(圖1),775 hPa急流所在高度低于1 km,依據(jù)低空急流最大風(fēng)速軸所在高度對其進(jìn)行分類,影響賀蘭山東麓暴雨過程的775 hPa、850 hPa急流均為邊界層急流,700 hPa急流為低空急流。在此基礎(chǔ)上,進(jìn)一步將影響賀蘭山東麓暴雨過程的低空急流過程分為七類,定義方法及其出現(xiàn)頻率統(tǒng)計(jì)見表3。

        表3 不同類型低空急流出現(xiàn)的頻率Table 3 Frequency of different types of low-level jets

        2006—2021年22次伴隨低空急流出現(xiàn)的賀蘭山東麓暴雨過程中有20次過程(約為90.9%)同時(shí)伴有低空急流和邊界層急流出現(xiàn),其中,三層低空急流型占比為54.5%,雙層低空急流Ⅰ型占比為36.5%,低空急流型和邊界層急流Ⅰ型過程各有1次(占4.5%),其他類型過程均沒有發(fā)生(表3)。這樣的分布表明:850 hPa有邊界層急流出現(xiàn),775 hPa 邊界層急流也一定會(huì)出現(xiàn),可見850 hPa急流對775 hPa急流的發(fā)生發(fā)展有著重要的影響。

        4 低空急流與賀蘭山東麓暴雨過程的相關(guān)性

        4.1 不同高度低空急流與暴雨過程的時(shí)間相關(guān)性

        從急流的發(fā)生頻率來看(圖4a):700、775、850 hPa三個(gè)高度上的急流總頻率分布呈單峰變化,峰值出現(xiàn)在00時(shí),13—14時(shí)頻率最低,其中700 hPa急流出現(xiàn)次數(shù)最多(793時(shí)次),775 hPa急流次之(597時(shí)次),850 hPa急流最少(212時(shí)次),700、775、850 hPa急流頻率的峰值分別出現(xiàn)在22—00時(shí)、23—00時(shí)、00時(shí),急流頻率最低值分別出現(xiàn)在13時(shí)和18時(shí)、12時(shí)和14時(shí)、12—14時(shí);分時(shí)段來看,上午(08—13時(shí))、下午(14—19時(shí))、前半夜(20—01時(shí))和后半夜(02—07時(shí))急流總頻率分別為350、348、464、440時(shí)次,各高度急流夜間發(fā)生頻率均較白天明顯偏多,700、775、850 hPa急流分別偏多57、93、56時(shí)次。這樣的急流日變化表明,急流的夜間增多特征顯著,前半夜急流發(fā)生最頻繁,下午急流出現(xiàn)最少,其中775 hPa夜間增多特征最為明顯。進(jìn)一步統(tǒng)計(jì)表明,22次暴雨過程的平均時(shí)長為18 h,700、775、850 hPa急流的平均持續(xù)時(shí)間分別為36、27、10 h,可見700 hPa、775 hPa急流的生命史普遍較暴雨過程更長,850 hPa 急流持續(xù)時(shí)間最短。

        注:圖c,圖d中,0表示低空急流與暴雨過程開始時(shí)間同步,-1、1分別表示低空急流建立較暴雨過程開始早、晚1 h,以此類推。圖4 22次賀蘭山東麓暴雨過程期間不同高度低空急流與暴雨過程的時(shí)間相關(guān)性(a)不同高度的低空急流頻率和低空急流總頻率及短時(shí)暴雨總頻率分布,(b)低空急流建立時(shí)間的頻率分布,(c)低空急流建立相對暴雨過程開始時(shí)間的超前/滯后關(guān)系,(d)低空急流最大風(fēng)速相對暴雨最大雨強(qiáng)時(shí)刻的超前/滯后關(guān)系Fig.4 Temporal correlation between low-level jets at different heights and rainstorm processes during the 22 rainstorm processes in the eastern foot of Helan Mountains: (a) distributions of frequency of low-level jets at different heights, total frequency of low-level jet and total frequency of short-time rainstorm, (b) distribution of frequency onset time of low-level jet, (c) lead-lag relationship between the onset of low-level jet and the start time of rainstorm process, (d) lead-lag relationship between the maximum wind speed of low-level jet and the maximum rainfall intensity of rainstorm

        從圖4a短時(shí)暴雨總頻率分布可以看出,短時(shí)暴雨總頻率呈雙峰結(jié)構(gòu),主峰值出現(xiàn)在22時(shí),次峰值出現(xiàn)在16時(shí)。19—22時(shí)短時(shí)暴雨頻率急劇增加,與此時(shí)低空急流增多有密切關(guān)系。下文討論急流建立、最大風(fēng)速時(shí)間與暴雨開始、最大雨強(qiáng)時(shí)間的對應(yīng)關(guān)系。

        從急流建立時(shí)間的頻率分布(圖4b)可以看出:低空急流普遍在中午到前半夜建立,其中,700 hPa和775 hPa 急流均在19時(shí)建立頻率最高,850 hPa急流在20時(shí)建立頻率最高。圖4c為低空急流建立相對暴雨過程開始的超前/滯后時(shí)間關(guān)系,可以看到:80%的低空急流在暴雨過程發(fā)生前建立,進(jìn)一步統(tǒng)計(jì)可知700、775、850 hPa急流建立較暴雨開始分別平均提前18、10、7 h。值得注意的是,只有10%的700 hPa急流和20%的775 hPa急流建立時(shí)間滯后于暴雨過程的開始時(shí)間,而接近一半的850 hPa急流(41.7%)在暴雨過程開始后才建立,因此700 hPa 和775 hPa急流的建立時(shí)間相對暴雨過程的開始時(shí)間有更好的指示作用。

        圖4d為低空急流最大風(fēng)速時(shí)間和暴雨過程最大雨強(qiáng)時(shí)間的超前/滯后關(guān)系??梢钥吹?40.7%的低空急流最大風(fēng)速較暴雨最大雨強(qiáng)發(fā)生時(shí)間提前,48%的急流最大風(fēng)速滯后于最大雨強(qiáng)時(shí)刻。結(jié)合圖4a,4b急流分布,表明低空急流與暴雨過程之間存在相互反饋,即急流增強(qiáng)時(shí),不僅促進(jìn)水汽輸送,還使得垂直切變增強(qiáng),超地轉(zhuǎn)現(xiàn)象明顯,不穩(wěn)定性增大,有利于對流發(fā)展,從而造成降水強(qiáng)度增大(陶詩言,1980);同時(shí),上升運(yùn)動(dòng)和凝結(jié)潛熱釋放所誘發(fā)的中尺度環(huán)流系統(tǒng)對低空急流的最大風(fēng)速起到了顯著的加強(qiáng)作用(斯公望,1994;Qian et al,2004)??傮w上,急流最大風(fēng)速超前最大雨強(qiáng)1 h到滯后1 h的這3個(gè)時(shí)次內(nèi)出現(xiàn)的頻率達(dá)到63%,這與王東阡和張耀存(2012)對中國東部西南低空急流日變化特征的研究結(jié)論一致,進(jìn)一步統(tǒng)計(jì)表明:700 hPa、775 hPa 急流最大風(fēng)速較暴雨最大雨強(qiáng)分別平均提前了54 min、18 min,而850 hPa急流最大風(fēng)速較暴雨最大雨強(qiáng)平均滯后12 min。

        上述分析表明,賀蘭山東麓暴雨的夜發(fā)性與低空急流的夜間增多密切相關(guān),且兩者存在相互反饋,其中700 hPa和775 hPa急流的建立、最大風(fēng)速時(shí)間對暴雨過程的開始、最大雨強(qiáng)時(shí)間的指示性更明顯。

        4.2 不同強(qiáng)度低空急流和短時(shí)暴雨的時(shí)間相關(guān)性

        根據(jù)最大風(fēng)速將低空急流強(qiáng)度分為3個(gè)等級(jí):10~14、14~18、18 m·s-1以上分別為1級(jí)、2級(jí)、3級(jí)急流。統(tǒng)計(jì)表明,2級(jí)急流出現(xiàn)頻率最高(654時(shí)次),1級(jí)急流次之(623時(shí)次),3級(jí)急流出現(xiàn)頻率最低(325時(shí)次)??梢?影響賀蘭山東麓暴雨過程的低空急流風(fēng)速普遍為10~18 m·s-1。從圖5可以看出,不同強(qiáng)度的低空急流頻率分布的日變化特征差異明顯,1級(jí)急流普遍在白天出現(xiàn)頻率高,夜間降低,而較強(qiáng)等級(jí)的2級(jí)、3級(jí)急流在夜間顯著增多,其中775 hPa和850 hPa的2級(jí)、3級(jí)急流的高頻時(shí)段集中在22時(shí)至01時(shí),而700 hPa的較強(qiáng)急流從23時(shí)持續(xù)到06時(shí)。這表明低空急流普遍在夜間強(qiáng)度較白天更強(qiáng),并且在前半夜迅速增強(qiáng),后半夜逐漸減弱。不同強(qiáng)度的短時(shí)暴雨頻率呈類似的分布,短時(shí)暴雨雨強(qiáng)以20~40 mm·h-1為主,20~40 mm·h-1短時(shí)暴雨頻率在16時(shí)達(dá)到最高,而40~60 mm·h-1和60 mm·h-1以上的短時(shí)暴雨均在22時(shí)出現(xiàn)次數(shù)最多,即短時(shí)暴雨在夜間強(qiáng)度更強(qiáng)。由此可見,賀蘭山東麓暴雨和低空急流的強(qiáng)度在時(shí)間上呈一致性變化,夜間隨著低空急流增強(qiáng),暴雨量級(jí)增大,而21時(shí)至01時(shí)20~40 mm·h-1短時(shí)暴雨維持較高頻率,也與這一時(shí)段內(nèi)低空急流強(qiáng)度明顯增強(qiáng)有關(guān)。但是,不同強(qiáng)度的低空急流和短時(shí)暴雨的頻率相關(guān)程度不同,其中,20~40 mm·h-1短時(shí)暴雨與850 hPa 的1級(jí)急流顯著正相關(guān)(相關(guān)系數(shù)為0.44),40~60 mm·h-1短時(shí)暴雨與775 hPa的2級(jí)急流相關(guān)系數(shù)為0.43,均通過了α=0.01的顯著性水平檢驗(yàn)。

        圖5 22次賀蘭山東麓暴雨過程期間不同強(qiáng)度的低空急流和短時(shí)暴雨頻率(a)700 hPa,(b)775 hPa,(c)850 hPaFig.5 Frequency of low-level jet and short-time rainstorm with different intensities during the 22 rainstorm processes in the eastern foot of Helan Mountains(a) 700 hPa, (b) 775 hPa, (c) 850 hPa

        圖6給出了暴雨過程的最大雨強(qiáng)與對應(yīng)時(shí)刻的三個(gè)關(guān)鍵區(qū)內(nèi)最大風(fēng)速和平均風(fēng)速分布,可以發(fā)現(xiàn),A區(qū)700 hPa急流和B區(qū)775 hPa急流的最大風(fēng)速和平均風(fēng)速變化趨勢一致性極高,且二者與最大雨強(qiáng)的分布也比較一致,其中最大雨強(qiáng)與A區(qū)700 hPa 急流的平均風(fēng)速相關(guān)程度最高,相關(guān)系數(shù)為0.22,通過了α=0.01的顯著性水平檢驗(yàn)。

        圖6 22次賀蘭山東麓暴雨過程期間最大雨強(qiáng)與對應(yīng)時(shí)刻急流關(guān)鍵區(qū)的(a)最大風(fēng)速和(b)平均風(fēng)速分布Fig.6 Distribution of the (a) maximum wind speed and (b) average wind speed of low-level jet in the key area corresponding to the maximum rainfall intensities during the 22 rainstorm processes in the eastern foot of Helan Mountains

        上述分析進(jìn)一步表明,賀蘭山東麓暴雨的夜間增強(qiáng)與低空急流夜間增強(qiáng)密切相關(guān),其中850 hPa的1級(jí)急流、775 hPa的2級(jí)急流分別對20~40 mm·h-1、40~60 mm·h-1短時(shí)暴雨的頻率指示性更強(qiáng),而700 hPa平均風(fēng)速對暴雨過程的最大雨強(qiáng)量級(jí)指示性更強(qiáng)。

        4.3 低空急流與暴雨過程的空間一致性

        賀蘭山東麓暴雨過程有量級(jí)大、強(qiáng)度強(qiáng)、降水效率高、對流性強(qiáng)、局地性和地形性顯著等特征(Chen et al,2021)。除個(gè)別持續(xù)性暴雨過程外,大部分對流性暴雨過程均存在雨強(qiáng)大于20 mm·h-1的站數(shù)明顯突增的時(shí)刻。借鑒羅輝等(2020)以20 mm·h-1降水站數(shù)的突增為時(shí)間節(jié)點(diǎn)對四川盆地暖區(qū)暴雨演變時(shí)段的劃分標(biāo)準(zhǔn),定義突增前為暴雨初始階段,其后為發(fā)展階段,突減時(shí)刻視為暴雨過程減弱階段。類似地,根據(jù)低空急流強(qiáng)度的突增和突減來定義低空急流的建立、發(fā)展和減弱階段。

        圖7給出了22次暴雨過程合成平均的不同階段降水量分布和低空急流的風(fēng)場分布。急流建立階段(圖7d),各關(guān)鍵區(qū)均出現(xiàn)了風(fēng)速小于11 m·s-1的弱急流。低空急流將充沛水汽和不穩(wěn)定能量輸向暴雨區(qū)(圖8a,8d,8g,8j),700 hPa和775 hPa上水汽輻合中心位于大于10 g·hPa-1·cm-1·s-1的水汽輸送帶前部(寧夏中部地區(qū)),暴雨區(qū)水汽輻散,只有850 hPa上位于水汽輻合區(qū)。大于500 J·kg-1的對流有效位能(CAPE)和大于31℃的K指數(shù)分布顯示暴雨區(qū)有較強(qiáng)的對流不穩(wěn)定能量,同時(shí),在對流層低層都有上升運(yùn)動(dòng)。相應(yīng)的該時(shí)段,賀蘭山東麓地區(qū)普遍出現(xiàn)5 mm以下的降水(圖7a)。

        急流發(fā)展階段(圖7e),700 hPa南風(fēng)急流(急流核13 m·s-1)前端北抬到775 hPa偏南急流核37.4°N 附近(最大風(fēng)速為12 m·s-1),775 hPa偏南急流前端也進(jìn)一步北抬至暴雨區(qū)東部,850 hPa東南急流核風(fēng)速也增大至12 m·s-1,前端北抬西進(jìn)至775 hPa急流入口區(qū)。隨著低空急流加強(qiáng)北抬,大于10 g·hPa-1·cm-1·s-1的水汽通量區(qū)也進(jìn)一步北抬,700 hPa的水汽通量輻合中心位于775 hPa 偏南急流關(guān)鍵區(qū)(B區(qū)),775 hPa水汽輻合中心也位于B區(qū)西部,對應(yīng)775 hPa急流核和700 hPa 南風(fēng)急流前端所在地,更加有利于低空急流將水汽向暴雨區(qū)輸送,且暴雨區(qū)在對流層低層都位于水汽輻合區(qū),其中850 hPa更是位于水汽輻合大值區(qū);隨著降水增強(qiáng),賀蘭山東麓地區(qū)的CAPE降低至300 J·kg-1,K指數(shù)變化不大,表明對流不穩(wěn)定能量釋放對應(yīng)暴雨增強(qiáng)階段,而賀蘭山東麓地區(qū)在低層的上升運(yùn)動(dòng)較暴雨開始階段明顯增強(qiáng)(圖8b,8e,8h,8k)??梢?低空急流的加強(qiáng)北抬或西進(jìn)進(jìn)一步增強(qiáng)了暴雨發(fā)生發(fā)展的環(huán)境條件。從圖7b 也可以看到,大于10 mm的降水區(qū)域已經(jīng)遍布整個(gè)賀蘭山東麓地區(qū),尤其是賀蘭山區(qū),賀蘭山沿山平羅到永寧段出現(xiàn)大于25 mm的強(qiáng)降水區(qū),有2個(gè)站超過30 mm,可見隨著急流加強(qiáng)北抬或西進(jìn),暴雨范圍擴(kuò)大、量級(jí)增強(qiáng)。

        急流減弱階段(圖7f),700 hPa轉(zhuǎn)為西南急流(急流核11 m·s-1)并南壓東退至河套東南部(前端位于36°N附近),775 hPa偏南急流和850 hPa東南急流基本東退出境。700 hPa和775 hPa水汽輸送明顯減弱東移(圖8c,8f,8i,8l),775 hPa和850 hPa 水汽輻合中心東移至內(nèi)蒙古鄂托克旗附近,與850 hPa東南急流出口區(qū)對應(yīng),除了700 hPa上暴雨區(qū)呈水汽弱輻合,暴雨區(qū)在775 hPa已轉(zhuǎn)為水汽輻散,暴雨區(qū)北部在850 hPa也轉(zhuǎn)為水汽輻散;同時(shí),暴雨區(qū)CAPE降至200 J·kg-1,K指數(shù)降低到30℃,垂直上升運(yùn)動(dòng)也減弱,對應(yīng)暴雨區(qū)降水量普遍低于10 mm(圖7c),暴雨過程趨于減弱結(jié)束。

        綜上所述,低空急流和賀蘭山東麓暴雨過程在空間分布上存在一致性變化,即隨著低空急流建立、加強(qiáng)北抬或西進(jìn)、減弱東退或南壓,賀蘭山東麓暴雨開始、增強(qiáng)、減弱,暴雨落區(qū)位于急流軸的左前方。低空急流為暴雨區(qū)帶來充沛水汽和不穩(wěn)定能量,有利于上升運(yùn)動(dòng)增強(qiáng),關(guān)鍵區(qū)B則是水汽和熱力不穩(wěn)定疊加區(qū)。當(dāng)700 hPa急流前端疊加在775 hPa急流核上方,850 hPa急流經(jīng)地形抬升后使775 hPa急流風(fēng)速增大,共同作用造成B區(qū)垂直風(fēng)切變增大,更加有利于急流將水汽和不穩(wěn)定能量進(jìn)一步向暴雨區(qū)輸送。

        5 低空急流影響賀蘭山東麓暴雨過程的發(fā)生發(fā)展機(jī)制初步探討

        2018年7月22日夜間賀蘭山東麓出現(xiàn)了特大暴雨過程,最大累計(jì)降水量和最大雨強(qiáng)均出現(xiàn)在賀蘭山滑雪場,分別為277.6 mm和74.1 mm·h-1(22時(shí))(表1),過程累計(jì)降水量位列寧夏歷史第一。本節(jié)以這次暴雨過程為例初步探討低空急流影響賀蘭山東麓暴雨過程發(fā)生發(fā)展的可能機(jī)制。

        22日14時(shí)(圖略),關(guān)鍵區(qū)A區(qū)700 hPa出現(xiàn)大于12 m·s-1的偏南急流,急流北端伸至寧夏中南部同心至海原一線。17時(shí)(圖9a),775 hPa急流和850 hPa急流開始建立,700 hPa急流核增強(qiáng)到14 m·s-1。19時(shí)急流發(fā)展北抬或西進(jìn)(圖9b,9e),700、775和850 hPa急流核分別增強(qiáng)到15、14和12 m·s-1,775 hPa急流軸前端北抬西進(jìn)至賀蘭山東坡。此時(shí),暴雨區(qū)附近在對流層中低層600 hPa以下存在東南西北向的次級(jí)環(huán)流,賀蘭山東側(cè)暴雨區(qū)位于上升運(yùn)動(dòng)大值區(qū),中心為-0.9 Pa·s-1(位于775 hPa附近),賀蘭山西側(cè)的下沉運(yùn)動(dòng)中心(0.6 Pa·s-1)也位于775 hPa附近,暴雨中心恰好處于上升運(yùn)動(dòng)中心偏西的梯度帶上,其上空700 hPa 附近有70%以上的相對濕度(圖9e)??梢?隨著低空急流北抬西進(jìn),降水落區(qū)從賀蘭山東麓的平原地區(qū)移至賀蘭山東坡山區(qū),降水強(qiáng)度顯著增大,20時(shí)有4個(gè)站出現(xiàn)10 mm以上降水,最大小時(shí)雨強(qiáng)為48.1 mm·h-1(圖9g,9h)。

        注:圖a中,黑色圓點(diǎn)表示暴雨中心滑雪場站位置,AB直線為剖線;圖e,圖f中,因資料分辨率低無法顯示賀蘭山地形,賀蘭山東西寬約30 km,故以白色三角形表示賀蘭山位置,白色矩形區(qū)域?yàn)橘R蘭山東麓地區(qū)。圖9 2018年7月22日暴雨過程期間的(a~c)急流演變(填色:700 hPa急流,紅色箭矢:775 hPa急流,綠色箭矢:850 hPa急流,單位:m·s-1),(e,f)沿圖9a中AB剖線的相對濕度(填色)、垂直速度(紅線,單位:Pa·s-1)、低空急流(綠線,單位:m·s-1)和風(fēng)場(風(fēng)矢,單位:m·s-1)的垂直剖面,(g~i)小時(shí)降水量分布(填色),及(d)22日17時(shí)至23日03時(shí)暴雨中心滑雪場站的相對濕度(填色)、假相當(dāng)位溫(黑線,單位:K)、垂直速度(紅線,單位:Pa·s-1)和風(fēng)場(風(fēng)羽)時(shí)間-高度剖面(a)17時(shí),(b,e)19時(shí),(c,f)21時(shí),(g)17—18時(shí),(h)19—20時(shí),(i)21—22時(shí)Fig.9 (a-c) Evolution of low-level jet (colored: 700 hPa jet, red vector: 775 hPa jet, green vector: 850 hPa jet, unit: m·s-1), (e, f) vertical cross-sections along line AB line of relative humidity (colored), vertical velocity (red contour, unit: Pa·s-1), low-level jet (green contour, unit: m·s-1) and wind vector (unit: m·s-1), and (g-i) distribution of hourly precipitation (colored) on 22 July 2018 respectively at (a) 17:00 BT, (b, e) 19:00 BT, (c, f) 21:00 BT, (g) 17:00-18:00 BT, (h) 19:00-20:00 BT, (i) 21:00-22:00 BT, and (d) time-height cross-section of relative humidity (colored), pseudo-equivalent potential temperature (black contour, unit: K), vertical velocity (red contour, unit: Pa·s-1) and wind barb at Huaxuechang Station in the rainstorm center from 17:00 BT 22 to 03:00 BT 23 July 2018

        21時(shí)(圖9c),低空急流進(jìn)一步增強(qiáng)并繼續(xù)北抬西進(jìn),700 hPa和775 hPa急流核均增強(qiáng)到17 m·s-1,三支急流的前端位置重合,均落在賀蘭山東側(cè)山區(qū),暴雨中心滑雪場站位于三支低空急流軸前端的風(fēng)速輻合大值區(qū);與19時(shí)相比,沿著急流前進(jìn)方向的次級(jí)環(huán)流明顯加強(qiáng),其中,850 hPa急流經(jīng)關(guān)鍵區(qū)B區(qū)地形抬升后,風(fēng)速明顯加大,急流核位于775 hPa 高度上,而后在B區(qū)的西北方向背風(fēng)坡氣流下沉(下沉速度中心為0.9 Pa·s-1,高度在775 hPa 附近),下沉支急流的前部落在暴雨區(qū),急流繼續(xù)北上西進(jìn),遇到賀蘭山地形阻擋,地形抬升和動(dòng)力抬升作用疊加,東坡山前的上升運(yùn)動(dòng)增強(qiáng),上升速度中心風(fēng)速達(dá)到最大,為-1.5 Pa·s-1(775 hPa高度附近),賀蘭山西側(cè)的下沉運(yùn)動(dòng)中心強(qiáng)度也增強(qiáng)到0.9 Pa·s-1,此時(shí),整個(gè)暴雨區(qū)在對流層中低層均為強(qiáng)烈的上升運(yùn)動(dòng)區(qū),暴雨中心仍處于上升運(yùn)動(dòng)中心偏西的梯度帶上(圖9f);同時(shí),顯著增強(qiáng)的低空急流還為暴雨區(qū)帶來了充沛的水汽,整個(gè)對流層中低層相對濕度明顯增大,700 hPa以下相對濕度大于80%,500 hPa以上相對濕度低于70%,呈上干下濕的分布(圖9d,9f)。對應(yīng)21—22時(shí)有3個(gè)站出現(xiàn)60 mm以上的強(qiáng)降水,最大小時(shí)雨強(qiáng)為74.1 mm·h-1(圖9i)。

        另外,從暴雨中心滑雪場站的時(shí)間-高度剖面(圖9d)可見,暴雨期間(22日19時(shí)至23日03時(shí)),對流層中低層600 hPa以下,(?θse/?p)<0,假相當(dāng)位溫θse的水平梯度≈0 K,700 hPa以下θse>350 K, 850 hPa以下θse>370 K,表明暴雨中心上空在對流層中下層高溫、高濕、高能,且存在強(qiáng)盛的對流不穩(wěn)定。

        從銀川雷達(dá)徑向速度圖(圖10a,10d)上也可以看到,從19:21至21:37,急流核強(qiáng)度由21.5 m·s-1增強(qiáng)到23.5 m·s-1,急流方向幾乎與山體垂直,隨著急流加強(qiáng)發(fā)展并進(jìn)一步北抬西進(jìn),暴雨中心滑雪場附近的氣旋式輻合也顯著增強(qiáng),旋轉(zhuǎn)速度由12 m·s-1增強(qiáng)到14.5 m·s-1(俞小鼎等,2006)。穩(wěn)定維持并加強(qiáng)發(fā)展的低空急流遇到賀蘭山東坡地形抬升后不斷觸發(fā)對流單體反復(fù)生消形成強(qiáng)烈的帶狀線性回波沿山體傳播(圖10b,10e),向北緩慢移動(dòng),形成列車效應(yīng)。并且60 dBz以上的強(qiáng)回波中心在經(jīng)過暴雨中心滑雪場站附近時(shí)維持了5個(gè)體掃時(shí)間(約30 min),對流活動(dòng)發(fā)展旺盛且穩(wěn)定少動(dòng),造成了暴雨中心滑雪場74.1 mm·h-1的強(qiáng)降水。反射率因子的剖面圖上可見強(qiáng)回波質(zhì)心有所下降(圖10c,10f),也表明了強(qiáng)降水效率增大。

        注:圖a,圖d中,黑色圓圈表示氣旋式輻合位置, 紅色圓圈表示急流核位置, 黑色箭頭指示急流方向;圖b,圖e中,五角星表示滑雪場測站位置, 白色直線表示剖線位置。圖10 2018年7月22日(a~c)19:21,(d~f)21:37銀川CA雷達(dá)的(a)1.5°和(d)2.4°仰角(a,d)徑向速度,(b,e)2.4°仰角反射率因子,及(c,f)沿白色剖線的反射率因子垂直剖面Fig.10 (a, d) Radial velocity of Yinchuan CA Radar at (a) 1.5° and (d) 2.4° elevations, (b, e) reflectivity factor at 2.4° elevation and (c, f) cross-section of reflectivity factor along the white line at (a-c) 19:21 BT and (d-f) 21:37 BT 22 July 2018

        可見,急流為暴雨區(qū)帶來了充沛的水汽,暴雨區(qū)呈上干下濕的不穩(wěn)定層結(jié),低層持續(xù)的高溫高濕高能為中尺度系統(tǒng)的觸發(fā)和維持提供了重要的熱動(dòng)力條件,急流遇到賀蘭山東坡地形抬升后不斷觸發(fā)對流單體反復(fù)生消形成強(qiáng)烈的帶狀線性回波沿山體傳播,經(jīng)過暴雨中心時(shí)穩(wěn)定少動(dòng),造成了局地強(qiáng)對流暴雨。

        6 結(jié)論和討論

        影響賀蘭山東麓暴雨過程的低空急流有三個(gè)關(guān)鍵區(qū),分別為河套南部、寧夏東南部和山西西南部,對應(yīng)700 hPa南風(fēng)急流、775 hPa偏南急流和850 hPa 東南急流。寧夏東南部的關(guān)鍵區(qū)作為三支低空急流匯合后繼續(xù)北上西進(jìn)的中轉(zhuǎn)站,對賀蘭山東麓暴雨過程的發(fā)生發(fā)展有著極其重要的影響。

        賀蘭山東麓暴雨過程同時(shí)受三層急流的影響較大,占比多達(dá)54.5%,其次是700和775 hPa急流同時(shí)出現(xiàn)的過程,占比為36.5%。

        賀蘭山東麓暴雨過程與低空急流在時(shí)間上存在一致性。暴雨的夜發(fā)性與夜間低空急流增強(qiáng)增多密切相關(guān),其中,700、775、850 hPa急流建立較暴雨開始平均提前18、10、7 h,700、775 hPa急流最大風(fēng)速較暴雨最大雨強(qiáng)分別平均提前54、18 min,而850 hPa 急流最大風(fēng)速較暴雨最大雨強(qiáng)平均滯后12 min。850 hPa的1級(jí)急流、775 hPa的2級(jí)急流頻率分別對20~40、40~60 mm·h-1的短時(shí)暴雨頻率指示性更強(qiáng),而河套南部關(guān)鍵區(qū)的700 hPa平均風(fēng)速對暴雨過程的最大雨強(qiáng)量級(jí)指示性更強(qiáng)。

        賀蘭山東麓暴雨過程和低空急流在空間上也存在一致性。隨著低空急流建立、加強(qiáng)北抬或西進(jìn)、減弱東退或南壓,賀蘭山東麓暴雨開始、增強(qiáng)、減弱,暴雨落區(qū)位于急流軸的左前方。低空急流北上西進(jìn)與賀蘭山地形結(jié)合,在東坡山前觸發(fā)多個(gè)對流單體形成帶狀線性回波沿山體傳播是造成賀蘭山區(qū)局地強(qiáng)對流暴雨的主要原因。

        本文系統(tǒng)地分析了低空急流與賀蘭山東麓暴雨過程的相關(guān)性,找到了影響暴雨過程的急流關(guān)鍵區(qū)和關(guān)鍵高度,揭示了低空急流與暴雨過程在時(shí)間和空間上的一致性變化,建立了急流相對暴雨開始時(shí)間的定量指標(biāo),初步探討了急流與地形結(jié)合影響暴雨中尺度系統(tǒng)發(fā)生發(fā)展的環(huán)境場和雷達(dá)特征,對暴雨預(yù)報(bào)預(yù)警有很好的指示意義,但由于再分析資料分辨率不夠高,無法反映實(shí)際的賀蘭山地形分布,對低空急流在地形抬升作用下觸發(fā)暴雨中尺度系統(tǒng)發(fā)展演變的分析不夠細(xì)致,尚需更深入的研究。下一步將選取典型暴雨過程進(jìn)行數(shù)值模擬,借助高分辨率數(shù)值模擬結(jié)果,深入研究地形作用下低空急流對暴雨中尺度系統(tǒng)觸發(fā)、發(fā)展和維持的影響機(jī)制。

        猜你喜歡
        雨強(qiáng)賀蘭山急流
        賀蘭山的鳥兒們
        幼兒100(2024年11期)2024-03-27 08:32:36
        賀蘭山巖畫:刻在石頭上的遠(yuǎn)古呼喚
        智海急流
        智海急流
        智海急流
        寧夏賀蘭山森林生態(tài)系統(tǒng)經(jīng)營技術(shù)
        智海急流
        靈魂只能獨(dú)行之六
        金秋(2019年24期)2019-06-17 11:31:00
        托里降雨時(shí)間分布特征分析
        雨強(qiáng)及坡度對黃土區(qū)草地坡面水流流速的影響
        人民黃河(2018年4期)2018-09-10 15:14:42
        337p日本欧洲亚洲大胆| 黑森林福利视频导航| 亚洲男人av天堂午夜在| 久久久久久av无码免费看大片| 久久亚洲国产欧洲精品一| 国产精品久久久一本精品| 国产一区二区高清不卡在线| 国产精品成人黄色大片| 日日麻批免费高清视频| 18禁免费无码无遮挡不卡网站| 久久久无码精品亚洲日韩按摩 | 国产一级二级三级在线观看av| 激情综合色五月丁香六月欧美| 精品无码中文字幕在线| 妺妺窝人体色www聚色窝韩国| 中文字幕无码免费久久9一区9| 丰满人妻一区二区三区免费| 国产亚洲午夜精品久久久| 国产免费又爽又色又粗视频| 亚洲av永久无码国产精品久久 | 日本黄色一区二区三区| 久久午夜av一区二区三区| 极品粉嫩嫩模大尺度无码视频| 无码gogo大胆啪啪艺术| 午夜理论片yy44880影院| 国产suv精品一区二区| 中国免费av网| 开心五月激情五月天天五月五月天 | 好看的日韩精品视频在线 | 亚洲成av在线免费不卡| 自拍成人免费在线视频| 亚洲乱码av中文一区二区| 亚洲国产日韩精品一区二区三区| 欧美大肥婆大肥bbbbb| 久久久久亚洲av无码专区| 最新国产美女一区二区三区| 国产影院一区二区在线| 伊人久久大香线蕉av波多野结衣 | 236宅宅理论片免费 | 久久道精品一区二区三区| 天堂av国产一区二区熟女人妻|