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

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

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

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

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

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

    引 言

    低空急流與暴雨的關(guān)系非常密切,因而得到了廣泛關(guān)注(陶詩言,1980;Higgins et al,1997;Chen et al,2005;Cook and Vizy,2010;鄭婧等,2020),相關(guān)研究可以有效提高暴雨預(yù)報的準(zhǔn)確率(Huang and Luo,2017)。低空急流對暴雨的影響主要體現(xiàn)在四個方面:(1)暴雨的日變化特征與低空急流的發(fā)展密切相關(guān)(王東阡和張耀存,2012;徐娟和陳勇明,2013;周靜等,2017),我國東部暖區(qū)暴雨和梅雨鋒暴雨過程在清晨出現(xiàn)的降水峰值與夜間西南低空急流加強有關(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)低空急流具有的強烈的不穩(wěn)定性使得急流軸上的風(fēng)速出現(xiàn)強風(fēng)速脈動傳播現(xiàn)象,從而觸發(fā)中尺度低渦及暴雨雨團發(fā)生(孫淑清和翟國慶,1980;張京英等,2005;趙嫻婷等,2020;胡雅君等,2020),低空急流在地形抬升的作用下造成強烈的垂直上升運動,激發(fā)對流云團強烈發(fā)展,也能夠觸發(fā)暴雨過程(Zhang and Meng,2019;黃小彥等,2020;汪小康等,2022;黃美金等,2022);(4)低空急流與暴雨是相互促進的,低空急流的存在有利于暴雨的發(fā)生,而暴雨的發(fā)生又促進了急流的形成和維持(賽瀚和苗峻峰,2012;劉鴻波等,2014)。孫繼松(2005)、Saulo et al(2007)分別證明了邊界層低空急流和局地暴雨、低空急流和中尺度對流及高空急流存在明顯的正反饋作用。

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

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

    1 資料和方法

    1.1 資 料

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

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

    參照Chen et al(2021)研究,定義賀蘭山東麓暴雨過程標(biāo)準(zhǔn):至少一個測站任意24 h累計雨量≥100 mm為一次暴雨過程,否則至少相鄰5個站≥25 mm,其中至少相鄰3個站≥50 mm;短時暴雨: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個經(jīng)距或緯距。

    考慮到銀川平原地區(qū)海拔高度在1102~1122 m(大約是875 hPa),選取ERA5資料的875、850、825、800、775、750、700、650和600 hPa共9個高度層進行低空急流的統(tǒng)計分析。定義急流最大頻率:賀蘭山東麓暴雨過程期間,某個高度層、某個格點上,急流出現(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ù)評估的方法,選取西北地區(qū)東部的6個探空站(銀川、延安、崆峒、武都、漢中和安康;圖1),對比探空和ERA5的低層風(fēng)場(700、850 hPa的風(fēng)向和風(fēng)速)。

    注:紅色方框為賀蘭山東麓地區(qū)。圖1 西北地區(qū)東部的探空站(紅色圓點)分布和地形高度(填色)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的顯著性水平檢驗)。風(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)向為負值,即相對于探空風(fēng)向,ERA5的風(fēng)向平均偏右,相反地,對于延安站、武都站和安康站而言,ERA5的風(fēng)向平均偏左;相對于探空風(fēng)速,ERA5的風(fēng)速普遍略偏大。風(fēng)向的相對誤差均低于8%,除武都站外,其他5個站的風(fēng)速相對誤差均低于16%,這可能與武都站地形有關(guān)。因此,ERA5低層風(fēng)場數(shù)據(jù)總體上是可信的,可用于研究暴雨過程與低空急流的相關(guān)性。

    表2 22次賀蘭山東麓暴雨過程期間37個時次的ERA5再分析數(shù)據(jù)與西北地區(qū)東部的6個探空站數(shù)據(jù)的評估結(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 影響賀蘭山東麓暴雨過程的低空急流活動關(guān)鍵區(qū)和關(guān)鍵高度

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

    注:圖a,圖b中黑色方框表示急流關(guān)鍵區(qū),A、B、C分別表示三個關(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三個關(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急流加強。由此可見:B區(qū)是三支低空急流匯合后繼續(xù)北上西進的中轉(zhuǎn)站,對賀蘭山東麓暴雨過程的發(fā)生發(fā)展有著極其重要的作用。

    注:藍色圓點表示合成平均風(fēng)速通過了α=0.01的顯著性水平檢驗,紅色方框為賀蘭山東麓地區(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)速軸所在高度對其進行分類,影響賀蘭山東麓暴雨過程的775 hPa、850 hPa急流均為邊界層急流,700 hPa急流為低空急流。在此基礎(chǔ)上,進一步將影響賀蘭山東麓暴雨過程的低空急流過程分為七類,定義方法及其出現(xiàn)頻率統(tǒng)計見表3。

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

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

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

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

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

    注:圖c,圖d中,0表示低空急流與暴雨過程開始時間同步,-1、1分別表示低空急流建立較暴雨過程開始早、晚1 h,以此類推。圖4 22次賀蘭山東麓暴雨過程期間不同高度低空急流與暴雨過程的時間相關(guān)性(a)不同高度的低空急流頻率和低空急流總頻率及短時暴雨總頻率分布,(b)低空急流建立時間的頻率分布,(c)低空急流建立相對暴雨過程開始時間的超前/滯后關(guān)系,(d)低空急流最大風(fēng)速相對暴雨最大雨強時刻的超前/滯后關(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短時暴雨總頻率分布可以看出,短時暴雨總頻率呈雙峰結(jié)構(gòu),主峰值出現(xiàn)在22時,次峰值出現(xiàn)在16時。19—22時短時暴雨頻率急劇增加,與此時低空急流增多有密切關(guān)系。下文討論急流建立、最大風(fēng)速時間與暴雨開始、最大雨強時間的對應(yīng)關(guān)系。

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

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

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

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

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

    圖5 22次賀蘭山東麓暴雨過程期間不同強度的低空急流和短時暴雨頻率(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給出了暴雨過程的最大雨強與對應(yīng)時刻的三個關(guān)鍵區(qū)內(nèi)最大風(fēng)速和平均風(fēng)速分布,可以發(fā)現(xiàn),A區(qū)700 hPa急流和B區(qū)775 hPa急流的最大風(fēng)速和平均風(fēng)速變化趨勢一致性極高,且二者與最大雨強的分布也比較一致,其中最大雨強與A區(qū)700 hPa 急流的平均風(fēng)速相關(guān)程度最高,相關(guān)系數(shù)為0.22,通過了α=0.01的顯著性水平檢驗。

    圖6 22次賀蘭山東麓暴雨過程期間最大雨強與對應(yīng)時刻急流關(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

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

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

    賀蘭山東麓暴雨過程有量級大、強度強、降水效率高、對流性強、局地性和地形性顯著等特征(Chen et al,2021)。除個別持續(xù)性暴雨過程外,大部分對流性暴雨過程均存在雨強大于20 mm·h-1的站數(shù)明顯突增的時刻。借鑒羅輝等(2020)以20 mm·h-1降水站數(shù)的突增為時間節(jié)點對四川盆地暖區(qū)暴雨演變時段的劃分標(biāo)準(zhǔn),定義突增前為暴雨初始階段,其后為發(fā)展階段,突減時刻視為暴雨過程減弱階段。類似地,根據(jù)低空急流強度的突增和突減來定義低空急流的建立、發(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ū)有較強的對流不穩(wěn)定能量,同時,在對流層低層都有上升運動。相應(yīng)的該時段,賀蘭山東麓地區(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偏南急流前端也進一步北抬至暴雨區(qū)東部,850 hPa東南急流核風(fēng)速也增大至12 m·s-1,前端北抬西進至775 hPa急流入口區(qū)。隨著低空急流加強北抬,大于10 g·hPa-1·cm-1·s-1的水汽通量區(qū)也進一步北抬,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ū);隨著降水增強,賀蘭山東麓地區(qū)的CAPE降低至300 J·kg-1,K指數(shù)變化不大,表明對流不穩(wěn)定能量釋放對應(yīng)暴雨增強階段,而賀蘭山東麓地區(qū)在低層的上升運動較暴雨開始階段明顯增強(圖8b,8e,8h,8k)??梢?低空急流的加強北抬或西進進一步增強了暴雨發(fā)生發(fā)展的環(huán)境條件。從圖7b 也可以看到,大于10 mm的降水區(qū)域已經(jīng)遍布整個賀蘭山東麓地區(qū),尤其是賀蘭山區(qū),賀蘭山沿山平羅到永寧段出現(xiàn)大于25 mm的強降水區(qū),有2個站超過30 mm,可見隨著急流加強北抬或西進,暴雨范圍擴大、量級增強。

    急流減弱階段(圖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)為水汽輻散;同時,暴雨區(qū)CAPE降至200 J·kg-1,K指數(shù)降低到30℃,垂直上升運動也減弱,對應(yīng)暴雨區(qū)降水量普遍低于10 mm(圖7c),暴雨過程趨于減弱結(jié)束。

    綜上所述,低空急流和賀蘭山東麓暴雨過程在空間分布上存在一致性變化,即隨著低空急流建立、加強北抬或西進、減弱東退或南壓,賀蘭山東麓暴雨開始、增強、減弱,暴雨落區(qū)位于急流軸的左前方。低空急流為暴雨區(qū)帶來充沛水汽和不穩(wěn)定能量,有利于上升運動增強,關(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)定能量進一步向暴雨區(qū)輸送。

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

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

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

    注:圖a中,黑色圓點表示暴雨中心滑雪場站位置,AB直線為剖線;圖e,圖f中,因資料分辨率低無法顯示賀蘭山地形,賀蘭山東西寬約30 km,故以白色三角形表示賀蘭山位置,白色矩形區(qū)域為賀蘭山東麓地區(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)小時降水量分布(填色),及(d)22日17時至23日03時暴雨中心滑雪場站的相對濕度(填色)、假相當(dāng)位溫(黑線,單位:K)、垂直速度(紅線,單位:Pa·s-1)和風(fēng)場(風(fēng)羽)時間-高度剖面(a)17時,(b,e)19時,(c,f)21時,(g)17—18時,(h)19—20時,(i)21—22時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時(圖9c),低空急流進一步增強并繼續(xù)北抬西進,700 hPa和775 hPa急流核均增強到17 m·s-1,三支急流的前端位置重合,均落在賀蘭山東側(cè)山區(qū),暴雨中心滑雪場站位于三支低空急流軸前端的風(fēng)速輻合大值區(qū);與19時相比,沿著急流前進方向的次級環(huán)流明顯加強,其中,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ù)北上西進,遇到賀蘭山地形阻擋,地形抬升和動力抬升作用疊加,東坡山前的上升運動增強,上升速度中心風(fēng)速達到最大,為-1.5 Pa·s-1(775 hPa高度附近),賀蘭山西側(cè)的下沉運動中心強度也增強到0.9 Pa·s-1,此時,整個暴雨區(qū)在對流層中低層均為強烈的上升運動區(qū),暴雨中心仍處于上升運動中心偏西的梯度帶上(圖9f);同時,顯著增強的低空急流還為暴雨區(qū)帶來了充沛的水汽,整個對流層中低層相對濕度明顯增大,700 hPa以下相對濕度大于80%,500 hPa以上相對濕度低于70%,呈上干下濕的分布(圖9d,9f)。對應(yīng)21—22時有3個站出現(xiàn)60 mm以上的強降水,最大小時雨強為74.1 mm·h-1(圖9i)。

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

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

    注:圖a,圖d中,黑色圓圈表示氣旋式輻合位置, 紅色圓圈表示急流核位置, 黑色箭頭指示急流方向;圖b,圖e中,五角星表示滑雪場測站位置, 白色直線表示剖線位置。圖10 2018年7月22日(a~c)19:21,(d~f)21:37銀川CA雷達的(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ā)和維持提供了重要的熱動力條件,急流遇到賀蘭山東坡地形抬升后不斷觸發(fā)對流單體反復(fù)生消形成強烈的帶狀線性回波沿山體傳播,經(jīng)過暴雨中心時穩(wěn)定少動,造成了局地強對流暴雨。

    6 結(jié)論和討論

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

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

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

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

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

    猜你喜歡
    雨強賀蘭山急流
    賀蘭山的鳥兒們
    幼兒100(2024年11期)2024-03-27 08:32:36
    賀蘭山巖畫:刻在石頭上的遠古呼喚
    華人時刊(2023年7期)2023-05-17 09:04:30
    智海急流
    智海急流
    智海急流
    寧夏賀蘭山森林生態(tài)系統(tǒng)經(jīng)營技術(shù)
    智海急流
    靈魂只能獨行之六
    金秋(2019年24期)2019-06-17 11:31:00
    托里降雨時間分布特征分析
    雨強及坡度對黃土區(qū)草地坡面水流流速的影響
    人民黃河(2018年4期)2018-09-10 15:14:42
    亚洲一卡2卡3卡4卡5卡精品中文| 热re99久久精品国产66热6| 亚洲成人免费av在线播放| 欧美日本中文国产一区发布| 麻豆国产av国片精品| 国产精品久久电影中文字幕 | 久久久久久久国产电影| 大香蕉久久网| tocl精华| 日日夜夜操网爽| 一级黄色大片毛片| 亚洲第一青青草原| 男人的好看免费观看在线视频 | 中文字幕精品免费在线观看视频| 人人妻,人人澡人人爽秒播| 亚洲少妇的诱惑av| 99在线人妻在线中文字幕 | 欧美乱码精品一区二区三区| 国产男靠女视频免费网站| 亚洲精品粉嫩美女一区| 视频区图区小说| 黄色怎么调成土黄色| 一级黄色大片毛片| 99精品欧美一区二区三区四区| 亚洲国产欧美一区二区综合| 伦理电影免费视频| 激情视频va一区二区三区| 视频区欧美日本亚洲| 天天添夜夜摸| 18禁美女被吸乳视频| 久久精品91无色码中文字幕| 后天国语完整版免费观看| 18禁裸乳无遮挡动漫免费视频| 成人18禁在线播放| 国产亚洲欧美精品永久| 国产三级黄色录像| 午夜福利免费观看在线| av网站免费在线观看视频| 男人操女人黄网站| 成人三级做爰电影| 亚洲精品国产色婷婷电影| 欧美久久黑人一区二区| 成人精品一区二区免费| av不卡在线播放| 精品国产亚洲在线| 久久青草综合色| 欧美日韩国产mv在线观看视频| 亚洲午夜精品一区,二区,三区| 老汉色∧v一级毛片| 亚洲美女黄片视频| 丰满的人妻完整版| 欧美日韩中文字幕国产精品一区二区三区 | 久热这里只有精品99| 一夜夜www| 亚洲成国产人片在线观看| 日韩免费高清中文字幕av| 99国产极品粉嫩在线观看| 人人妻人人澡人人爽人人夜夜| 视频区图区小说| 99国产精品99久久久久| 91字幕亚洲| 国产精品.久久久| av福利片在线| 亚洲专区字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 电影成人av| 久久久久精品国产欧美久久久| 九色亚洲精品在线播放| 亚洲久久久国产精品| 大香蕉久久成人网| 国产成人免费无遮挡视频| 999久久久国产精品视频| 日本撒尿小便嘘嘘汇集6| 搡老乐熟女国产| 在线观看免费视频日本深夜| 日韩欧美三级三区| 如日韩欧美国产精品一区二区三区| 91老司机精品| 男人舔女人的私密视频| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索| 97人妻天天添夜夜摸| 精品国产美女av久久久久小说| 午夜福利在线观看吧| 男女高潮啪啪啪动态图| 丝袜在线中文字幕| 露出奶头的视频| 日本vs欧美在线观看视频| 无限看片的www在线观看| 91成人精品电影| 中文字幕另类日韩欧美亚洲嫩草| 欧洲精品卡2卡3卡4卡5卡区| av免费在线观看网站| 另类亚洲欧美激情| 国产av又大| 亚洲色图综合在线观看| 精品亚洲成国产av| 一级片免费观看大全| 国产高清videossex| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 精品久久久久久电影网| 欧美国产精品va在线观看不卡| 一a级毛片在线观看| 亚洲av欧美aⅴ国产| 国内久久婷婷六月综合欲色啪| 国产视频一区二区在线看| 欧美激情 高清一区二区三区| 欧美日韩成人在线一区二区| 亚洲欧美一区二区三区久久| 成年动漫av网址| 久久精品成人免费网站| 丁香六月欧美| 午夜福利视频在线观看免费| 一边摸一边抽搐一进一小说 | 亚洲九九香蕉| 国产精品久久电影中文字幕 | 欧美在线一区亚洲| 亚洲精品国产色婷婷电影| 身体一侧抽搐| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩一区二区三区精品不卡| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 久久影院123| 欧美 日韩 精品 国产| 女性生殖器流出的白浆| 欧美另类亚洲清纯唯美| 热re99久久精品国产66热6| 国产精品.久久久| 18禁裸乳无遮挡动漫免费视频| 国产乱人伦免费视频| 两个人免费观看高清视频| 高清毛片免费观看视频网站 | 韩国av一区二区三区四区| 美女午夜性视频免费| 又大又爽又粗| 91在线观看av| 国产精品久久久人人做人人爽| 国产色视频综合| 久久 成人 亚洲| 极品少妇高潮喷水抽搐| 男人的好看免费观看在线视频 | 午夜免费成人在线视频| 久久久久久久久久久久大奶| 国产1区2区3区精品| 成人国产一区最新在线观看| 色精品久久人妻99蜜桃| 国产激情欧美一区二区| √禁漫天堂资源中文www| 国产免费男女视频| 国产高清激情床上av| 大香蕉久久成人网| 亚洲中文日韩欧美视频| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 天堂√8在线中文| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 超碰成人久久| 99国产极品粉嫩在线观看| 国产一区二区激情短视频| 在线观看免费日韩欧美大片| 高清av免费在线| 久久精品亚洲精品国产色婷小说| 韩国精品一区二区三区| 99re6热这里在线精品视频| 亚洲精品自拍成人| 高清av免费在线| av不卡在线播放| 日本五十路高清| 成人亚洲精品一区在线观看| 欧美乱妇无乱码| 免费久久久久久久精品成人欧美视频| 香蕉久久夜色| 韩国精品一区二区三区| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 成人18禁高潮啪啪吃奶动态图| 老司机在亚洲福利影院| 免费少妇av软件| 黄色视频,在线免费观看| aaaaa片日本免费| 黄色怎么调成土黄色| 国产一区二区三区视频了| 国产91精品成人一区二区三区| 亚洲欧美一区二区三区黑人| 精品一区二区三区四区五区乱码| 黑人欧美特级aaaaaa片| 国产亚洲av高清不卡| 精品一品国产午夜福利视频| videosex国产| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲色图综合在线观看| 天堂俺去俺来也www色官网| 国产亚洲精品第一综合不卡| 亚洲精品美女久久av网站| 国产欧美日韩一区二区精品| 久久国产亚洲av麻豆专区| 国产99白浆流出| 99re6热这里在线精品视频| 欧美精品一区二区免费开放| 国产片内射在线| 一二三四在线观看免费中文在| 热99久久久久精品小说推荐| 老鸭窝网址在线观看| 国产精品亚洲av一区麻豆| 在线国产一区二区在线| 伊人久久大香线蕉亚洲五| 日韩熟女老妇一区二区性免费视频| 香蕉国产在线看| 搡老乐熟女国产| 久久国产精品影院| 国产日韩欧美亚洲二区| 国产一区二区三区在线臀色熟女 | 精品久久久久久电影网| 精品一区二区三区av网在线观看| 热re99久久精品国产66热6| 99国产精品免费福利视频| 国产精品偷伦视频观看了| 正在播放国产对白刺激| 亚洲精品成人av观看孕妇| 精品国产美女av久久久久小说| a在线观看视频网站| 精品第一国产精品| 999久久久精品免费观看国产| 一二三四社区在线视频社区8| 别揉我奶头~嗯~啊~动态视频| xxxhd国产人妻xxx| 欧美+亚洲+日韩+国产| 国产精品一区二区在线观看99| x7x7x7水蜜桃| 热re99久久国产66热| 国产亚洲精品第一综合不卡| а√天堂www在线а√下载 | 免费少妇av软件| 在线看a的网站| 一级片'在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 无人区码免费观看不卡| 亚洲 欧美一区二区三区| 天天躁夜夜躁狠狠躁躁| 久久久久久久午夜电影 | 久久久久视频综合| 欧美乱色亚洲激情| 亚洲片人在线观看| 亚洲aⅴ乱码一区二区在线播放 | 99久久国产精品久久久| 丁香六月欧美| av中文乱码字幕在线| 欧美激情高清一区二区三区| 高清av免费在线| 亚洲精品成人av观看孕妇| 精品人妻熟女毛片av久久网站| 国产精品美女特级片免费视频播放器 | 亚洲全国av大片| 成人亚洲精品一区在线观看| aaaaa片日本免费| 欧美日韩国产mv在线观看视频| 国产成人系列免费观看| 国产精品 国内视频| 久久影院123| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 欧美丝袜亚洲另类 | 99久久综合精品五月天人人| 曰老女人黄片| 亚洲精品久久成人aⅴ小说| 亚洲精华国产精华精| 中文字幕最新亚洲高清| 黄色 视频免费看| 色94色欧美一区二区| 亚洲成a人片在线一区二区| 日韩免费高清中文字幕av| 国产精品香港三级国产av潘金莲| 一边摸一边抽搐一进一小说 | 日本精品一区二区三区蜜桃| 精品熟女少妇八av免费久了| 精品高清国产在线一区| 精品人妻在线不人妻| 女人被狂操c到高潮| 女人被躁到高潮嗷嗷叫费观| 欧美乱妇无乱码| 精品福利观看| 男女下面插进去视频免费观看| 精品国产国语对白av| 宅男免费午夜| 色尼玛亚洲综合影院| 免费黄频网站在线观看国产| 免费在线观看黄色视频的| av免费在线观看网站| 亚洲性夜色夜夜综合| 亚洲 欧美一区二区三区| 热99久久久久精品小说推荐| 高清在线国产一区| 国产av一区二区精品久久| 久久久久久亚洲精品国产蜜桃av| 国产深夜福利视频在线观看| 久9热在线精品视频| 亚洲专区国产一区二区| 中亚洲国语对白在线视频| 美国免费a级毛片| 国产午夜精品久久久久久| 成熟少妇高潮喷水视频| 国产又色又爽无遮挡免费看| 欧美成人免费av一区二区三区 | 欧美日韩福利视频一区二区| 精品国产亚洲在线| 国产精华一区二区三区| 亚洲中文字幕日韩| 在线观看午夜福利视频| 女同久久另类99精品国产91| 满18在线观看网站| 一边摸一边抽搐一进一出视频| 午夜福利在线免费观看网站| 亚洲av成人一区二区三| 天堂√8在线中文| 久久热在线av| 久热这里只有精品99| 精品久久蜜臀av无| 欧美在线黄色| 熟女少妇亚洲综合色aaa.| 亚洲成人免费av在线播放| 久久这里只有精品19| 午夜精品国产一区二区电影| 一进一出抽搐动态| 飞空精品影院首页| 麻豆国产av国片精品| 激情视频va一区二区三区| 免费看十八禁软件| 国产精品亚洲一级av第二区| 午夜福利在线观看吧| 成年版毛片免费区| 亚洲精华国产精华精| 曰老女人黄片| a级毛片在线看网站| 9热在线视频观看99| 天堂中文最新版在线下载| 天天躁日日躁夜夜躁夜夜| 日韩熟女老妇一区二区性免费视频| 人妻丰满熟妇av一区二区三区 | svipshipincom国产片| av在线播放免费不卡| 99久久综合精品五月天人人| 日韩 欧美 亚洲 中文字幕| 日本vs欧美在线观看视频| 国产深夜福利视频在线观看| 色婷婷av一区二区三区视频| 波多野结衣av一区二区av| 亚洲欧美激情综合另类| 极品少妇高潮喷水抽搐| 国产单亲对白刺激| 亚洲第一av免费看| 美女扒开内裤让男人捅视频| 久久精品aⅴ一区二区三区四区| 欧美一级毛片孕妇| 午夜福利乱码中文字幕| 夫妻午夜视频| 涩涩av久久男人的天堂| 久久精品亚洲av国产电影网| 亚洲欧美一区二区三区黑人| 99久久精品国产亚洲精品| 久久香蕉精品热| videosex国产| www.精华液| 自线自在国产av| 免费看十八禁软件| 久久精品亚洲熟妇少妇任你| 亚洲av成人av| 老司机午夜十八禁免费视频| av免费在线观看网站| 亚洲成av片中文字幕在线观看| 国产精品美女特级片免费视频播放器 | 国产国语露脸激情在线看| 欧美国产精品va在线观看不卡| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 在线观看www视频免费| 国产伦人伦偷精品视频| 91老司机精品| 国产一区在线观看成人免费| av福利片在线| 黄色成人免费大全| 日韩欧美免费精品| 国产成人av激情在线播放| 亚洲欧美色中文字幕在线| 国产国语露脸激情在线看| 久久香蕉精品热| 久久99一区二区三区| 日日夜夜操网爽| 深夜精品福利| 久久国产亚洲av麻豆专区| 成人亚洲精品一区在线观看| 中文字幕精品免费在线观看视频| 天天操日日干夜夜撸| 亚洲色图av天堂| 亚洲一区二区三区欧美精品| 丰满饥渴人妻一区二区三| 大香蕉久久成人网| 中文字幕制服av| 热99re8久久精品国产| av天堂在线播放| 夜夜夜夜夜久久久久| 欧美精品高潮呻吟av久久| 91麻豆av在线| 热99久久久久精品小说推荐| 两个人看的免费小视频| 美女高潮到喷水免费观看| 精品卡一卡二卡四卡免费| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品久久成人aⅴ小说| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩精品亚洲av| 99国产精品一区二区蜜桃av | 精品人妻1区二区| 一级,二级,三级黄色视频| 少妇猛男粗大的猛烈进出视频| 天天躁狠狠躁夜夜躁狠狠躁| 99精品久久久久人妻精品| 在线观看免费视频日本深夜| 久久精品国产亚洲av高清一级| 少妇猛男粗大的猛烈进出视频| 欧美老熟妇乱子伦牲交| 夜夜夜夜夜久久久久| 亚洲精品成人av观看孕妇| 亚洲欧美一区二区三区黑人| 搡老乐熟女国产| av天堂在线播放| 别揉我奶头~嗯~啊~动态视频| 欧美成人午夜精品| 久久精品成人免费网站| 久久国产精品影院| 日韩人妻精品一区2区三区| 色综合婷婷激情| 天堂动漫精品| 午夜91福利影院| 老司机福利观看| 午夜91福利影院| 亚洲av成人一区二区三| 啦啦啦 在线观看视频| av一本久久久久| 精品一品国产午夜福利视频| 成熟少妇高潮喷水视频| 国产日韩欧美亚洲二区| 久久性视频一级片| 老司机深夜福利视频在线观看| 深夜精品福利| 在线播放国产精品三级| 无遮挡黄片免费观看| 男人操女人黄网站| 女同久久另类99精品国产91| 中文字幕人妻丝袜制服| 国产成人精品久久二区二区免费| 亚洲精品一二三| 国产成人欧美在线观看 | 亚洲一码二码三码区别大吗| 99国产精品99久久久久| 欧美黄色片欧美黄色片| 国产精品欧美亚洲77777| 亚洲第一av免费看| e午夜精品久久久久久久| 日韩免费av在线播放| 嫁个100分男人电影在线观看| 亚洲国产欧美网| 在线天堂中文资源库| 久久午夜综合久久蜜桃| 99热只有精品国产| 午夜日韩欧美国产| 国产成人精品无人区| 国产亚洲一区二区精品| 黄色片一级片一级黄色片| 天天影视国产精品| 两个人看的免费小视频| 大陆偷拍与自拍| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 99久久人妻综合| 大片电影免费在线观看免费| 一二三四社区在线视频社区8| 欧美日韩国产mv在线观看视频| 久久久久久久久久久久大奶| 午夜视频精品福利| 欧美一级毛片孕妇| 亚洲成av片中文字幕在线观看| 久久精品国产综合久久久| 亚洲黑人精品在线| 国产高清视频在线播放一区| 精品欧美一区二区三区在线| 深夜精品福利| 国产成+人综合+亚洲专区| 亚洲国产欧美一区二区综合| 国产欧美日韩一区二区三区在线| 中文字幕色久视频| 黄片小视频在线播放| 99久久99久久久精品蜜桃| 久久久久久久久久久久大奶| 日韩免费av在线播放| 777米奇影视久久| 深夜精品福利| 欧美人与性动交α欧美软件| 国产精华一区二区三区| 国产亚洲精品第一综合不卡| 国产精品一区二区在线不卡| 黑丝袜美女国产一区| 97人妻天天添夜夜摸| 欧美日韩福利视频一区二区| 国产淫语在线视频| 脱女人内裤的视频| 亚洲国产精品一区二区三区在线| 丝袜在线中文字幕| 人妻久久中文字幕网| 一本一本久久a久久精品综合妖精| 咕卡用的链子| 99riav亚洲国产免费| 99国产极品粉嫩在线观看| 黑人猛操日本美女一级片| 女性生殖器流出的白浆| 97人妻天天添夜夜摸| 国产精品美女特级片免费视频播放器 | 无限看片的www在线观看| 嫁个100分男人电影在线观看| 99国产极品粉嫩在线观看| 在线播放国产精品三级| 老司机影院毛片| 又紧又爽又黄一区二区| aaaaa片日本免费| 热99久久久久精品小说推荐| netflix在线观看网站| 视频在线观看一区二区三区| 婷婷精品国产亚洲av在线 | 久久精品人人爽人人爽视色| 一二三四在线观看免费中文在| 无人区码免费观看不卡| 少妇被粗大的猛进出69影院| 精品国产国语对白av| 一边摸一边抽搐一进一小说 | 男女高潮啪啪啪动态图| 国产一区二区三区在线臀色熟女 | 精品乱码久久久久久99久播| 国产在视频线精品| 久久精品亚洲av国产电影网| 亚洲欧美色中文字幕在线| 高潮久久久久久久久久久不卡| 香蕉丝袜av| 欧美日韩黄片免| 精品欧美一区二区三区在线| 在线天堂中文资源库| 欧美不卡视频在线免费观看 | av电影中文网址| 久久九九热精品免费| 国产精品久久久av美女十八| 国产99久久九九免费精品| 国产成人精品无人区| 午夜精品国产一区二区电影| 丰满人妻熟妇乱又伦精品不卡| 国产精品一区二区在线观看99| 亚洲精品成人av观看孕妇| 两人在一起打扑克的视频| 宅男免费午夜| 不卡av一区二区三区| 欧美日韩一级在线毛片| 日韩欧美一区视频在线观看| 久久精品国产综合久久久| 一区二区三区精品91| 在线观看一区二区三区激情| 免费观看精品视频网站| 久久亚洲真实| 日韩制服丝袜自拍偷拍| 欧美大码av| 黄频高清免费视频| 亚洲精品国产区一区二| 免费久久久久久久精品成人欧美视频| 久久精品人人爽人人爽视色| 亚洲欧美激情在线| 久久精品国产99精品国产亚洲性色 | 色综合婷婷激情| 超色免费av| 巨乳人妻的诱惑在线观看| av国产精品久久久久影院| 亚洲黑人精品在线| 午夜免费成人在线视频| 淫妇啪啪啪对白视频| 亚洲七黄色美女视频| 男人舔女人的私密视频| 国产aⅴ精品一区二区三区波| 免费女性裸体啪啪无遮挡网站| 午夜视频精品福利| 久久久久精品国产欧美久久久| 99国产综合亚洲精品| 久久婷婷成人综合色麻豆| √禁漫天堂资源中文www| 俄罗斯特黄特色一大片| 国产成人欧美| 99久久综合精品五月天人人| 日韩熟女老妇一区二区性免费视频| 亚洲精品成人av观看孕妇| 久久草成人影院| 国产欧美亚洲国产| 黑人巨大精品欧美一区二区mp4| 国产成人欧美在线观看 | 一边摸一边抽搐一进一出视频| 日本wwww免费看| 国产精品二区激情视频| 久久久久国产精品人妻aⅴ院 | 80岁老熟妇乱子伦牲交| 久久精品亚洲精品国产色婷小说| 99精品久久久久人妻精品| 国产不卡av网站在线观看| 国产精品乱码一区二三区的特点 | 两性夫妻黄色片| 久久精品aⅴ一区二区三区四区| 日韩免费高清中文字幕av| 日韩有码中文字幕| 国产一卡二卡三卡精品| 超碰成人久久|