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

    颮線組織化過程對環(huán)境垂直風(fēng)切變和水汽的響應(yīng)

    2016-07-27 09:32:47張建軍王詠青鐘瑋1湖南省氣候中心長沙410118南京信息工程大學(xué)大氣科學(xué)學(xué)院太平洋臺風(fēng)研究中心南京100443解放軍理工大學(xué)氣象海洋學(xué)院南京111014南京大學(xué)中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室南京10093
    大氣科學(xué) 2016年4期
    關(guān)鍵詞:數(shù)值模擬

    張建軍王詠青鐘瑋1湖南省氣候中心,長沙410118南京信息工程大學(xué)大氣科學(xué)學(xué)院太平洋臺風(fēng)研究中心,南京100443解放軍理工大學(xué)氣象海洋學(xué)院,南京111014南京大學(xué)中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室,南京10093

    ?

    颮線組織化過程對環(huán)境垂直風(fēng)切變和水汽的響應(yīng)

    張建軍1, 2王詠青2鐘瑋3, 4
    1湖南省氣候中心,長沙410118
    2南京信息工程大學(xué)大氣科學(xué)學(xué)院太平洋臺風(fēng)研究中心,南京210044
    3解放軍理工大學(xué)氣象海洋學(xué)院,南京211101
    4南京大學(xué)中尺度災(zāi)害性天氣教育部重點(diǎn)實(shí)驗(yàn)室,南京210093

    摘 要利用ARPS模式對颮線發(fā)生發(fā)展過程進(jìn)行二維理想數(shù)值試驗(yàn),討論了低層環(huán)境垂直風(fēng)切變和水汽條件變化時(shí),颮線內(nèi)部物理因子配置變化及其與系統(tǒng)強(qiáng)度演變的聯(lián)系。研究表明,颮線發(fā)展過程中出現(xiàn)的動(dòng)量、熱量和水汽的再分配過程,造成系統(tǒng)內(nèi)垂直環(huán)流結(jié)構(gòu)和擾動(dòng)溫濕場分布發(fā)生變化,從而影響系統(tǒng)內(nèi)部深對流的組織化過程和颮線強(qiáng)度的發(fā)展。基于低層環(huán)境垂直風(fēng)切變和水汽兩個(gè)要素的敏感性試驗(yàn)研究表明,低層環(huán)境垂直風(fēng)切變增大(減?。r(shí),颮線移速減慢(加快),冷池前沿激發(fā)的新對流與中高層的垂直運(yùn)動(dòng)相互貫通(分離),颮線系統(tǒng)強(qiáng)度隨之增強(qiáng)(減弱)。此外,當(dāng)?shù)蛯铀黾樱p少)時(shí),會導(dǎo)致輸送到中層的水汽增加(減少),中層凝結(jié)潛熱釋放增多(減少),該層垂直運(yùn)動(dòng)增強(qiáng)(減弱);同時(shí),颮線系統(tǒng)區(qū)域環(huán)境釋放的對流有效位能(CAPE)增大(減小),新生對流的強(qiáng)度增強(qiáng)(減弱)。低層水汽條件通過水汽輸送和能量釋放,改變冷池前沿新對流與中高層垂直環(huán)流的組織化結(jié)構(gòu),從而影響颮線強(qiáng)度。

    關(guān)鍵詞颮線 數(shù)值模擬 冷池 環(huán)境垂直風(fēng)切變 水汽1

    張建軍, 王詠青, 鐘瑋. 2016. 颮線組織化過程對環(huán)境垂直風(fēng)切變和水汽的響應(yīng) [J]. 大氣科學(xué), 40 (4): 689–702. Zhang Jianjun, Wang Yongqing, Zhong Wei. 2016. Impact of vertical wind shear and moisture on the organization of squall lines [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 40 (4): 689–702, doi:10.3878/j.issn.1006-9895.1505.14337.

    1 引言

    颮線是有組織的雷暴群帶狀排列形成的中尺度對流系統(tǒng),往往伴隨有局地暴雨、短時(shí)大風(fēng)、冰雹等強(qiáng)災(zāi)害性的強(qiáng)對流天氣,維持時(shí)間有時(shí)能夠達(dá)到或超過十幾小時(shí)。對于這類有較長生命史的災(zāi)害性天氣現(xiàn)象,探討其在不同環(huán)境條件下的發(fā)展變化過程對防災(zāi)減災(zāi)具有重要作用。

    早在20世紀(jì)50年代,Newton(1950)指出水平風(fēng)速存在垂直切變的條件下,低層輻合高層輻散的配置有利于雷暴的維持和發(fā)展。后續(xù)研究表明垂直于颮線的環(huán)境垂直風(fēng)切變是颮線發(fā)展演變的關(guān)鍵因子(Fujita, 1955)。統(tǒng)計(jì)資料也顯示大多數(shù)颮線的環(huán)境垂直風(fēng)切變與颮線正交(Bluestein and Jain, 1985; Bluestein et al., 1987; Wyss and Emanuel, 1988)。近三十年來,颮線的強(qiáng)對流運(yùn)動(dòng)被認(rèn)為是冷池與低層(0~2.5 km)環(huán)境垂直風(fēng)切變相互作用的結(jié)果更是得到了廣泛認(rèn)可(Thorpe et al., 1982; Rotunno et al., 1988;Weisman and Rotunno, 2004; Takemi, 2006, 2007, 2014; Coniglio et al., 2012;陳明軒和王迎春,2012)。其中Rotunno et al. (1988) 提出的理論被稱之為RKW理論,RKW理論認(rèn)為低層環(huán)境垂直風(fēng)切變與近地面冷池強(qiáng)度相當(dāng)時(shí),颮線強(qiáng)度最強(qiáng),兩者相差較大時(shí)則颮線減弱。

    除了低層環(huán)境垂直風(fēng)切變對颮線的發(fā)展維持至關(guān)重要外,其他環(huán)境因子的作用也不容忽視。Takemi(2006)指出在較濕的環(huán)境下,更強(qiáng)的切變更有利于強(qiáng)對流的發(fā)展,且保持環(huán)境中水汽總量一致,低層更濕的環(huán)境所生成的颮線強(qiáng)度會更強(qiáng)。Takemi(2007)還指出在對流不穩(wěn)定層中,靜力穩(wěn)定度在決定颮線強(qiáng)度方面起著重要作用;在相同靜力穩(wěn)定度環(huán)境下,對流有效位能(CAPE)及可降水量能有效的表示颮線強(qiáng)度。Parker(2008)通過不斷減小近地面層的溫度使得颮線所處環(huán)境的溫度要素不斷發(fā)生變化,并探討了這種條件下颮線的發(fā)展變化過程。

    以往的理想數(shù)值試驗(yàn)及個(gè)例分析,研究了環(huán)境垂直風(fēng)切變、水汽等環(huán)境要素在特定配置條件下,颮線觸發(fā)和組織化的物理過程,并通過冷池強(qiáng)度與環(huán)境垂直風(fēng)切變的比值(簡稱為冷切比)來解釋颮線發(fā)展變化的原因,并未從颮線內(nèi)部環(huán)流結(jié)構(gòu)的變化出發(fā)來探討颮線強(qiáng)度變化的機(jī)理。實(shí)際上在颮線發(fā)展演變過程中,系統(tǒng)內(nèi)部出現(xiàn)了動(dòng)量、熱量和水汽的輸送、集中等再分配過程,并依此對颮線組織化程度產(chǎn)生影響。本文通過理想數(shù)值試驗(yàn),研究適宜環(huán)境配置條件下,颮線觸發(fā)和發(fā)展階段系統(tǒng)內(nèi)部動(dòng)量、熱量和水汽的再分配過程,及其對颮線系統(tǒng)強(qiáng)度變化的影響,并在此基礎(chǔ)上討論當(dāng)?shù)蛯迎h(huán)境垂直風(fēng)切變和水汽等環(huán)境要素場發(fā)生變化時(shí),系統(tǒng)內(nèi)部各物理因子的響應(yīng)過程及其對颮線系統(tǒng)強(qiáng)度變化的影響。文章第二部分基于RKW理論最優(yōu)環(huán)境條件配置模擬了颮線的發(fā)生發(fā)展過程,深入分析了系統(tǒng)內(nèi)部動(dòng)量、熱量和水汽的再分配過程及其與颮線系統(tǒng)強(qiáng)度變化之間的聯(lián)系;第三、四部分則分別介紹颮線發(fā)展過程中其強(qiáng)度對低層環(huán)境垂直風(fēng)切變和低層水汽的響應(yīng)結(jié)果,并通過分析颮線的內(nèi)部環(huán)流結(jié)構(gòu)探討了颮線強(qiáng)度變化的原因;最后對本文進(jìn)行了歸納總結(jié)。

    2 颮線觸發(fā)及發(fā)展過程的數(shù)值試驗(yàn)

    觀測研究表明,颮線具有明顯的準(zhǔn)二維特征,即沿颮線方向上各個(gè)風(fēng)暴單體的發(fā)展演變特征比較相似(Thorpe et al., 1982; Fovell and Ogura, 1988; Houze et al., 1989; Parker, 1998; Parker and Johnson, 2004),因此本文采用Xue et al.(2000)開發(fā)的非靜力風(fēng)暴尺度高級區(qū)域預(yù)報(bào)系統(tǒng)(ARPS)模式(5.3.4版本),設(shè)計(jì)了二維理想數(shù)值試驗(yàn),再現(xiàn)颮線的觸發(fā)和發(fā)展過程,以及颮線發(fā)展過程中系統(tǒng)內(nèi)部物理量再分配的物理圖像。

    2.1 試驗(yàn)設(shè)計(jì)

    假定理想二維剖面為x–z平面,水平(x)方向區(qū)域?yàn)?60 km,分辨率為600 m;垂直(z)高度為18 km,為盡可能節(jié)省計(jì)算時(shí)間并模擬出精細(xì)的颮線低層垂直結(jié)構(gòu),垂直方向上進(jìn)行了網(wǎng)格拉伸,即在z=3 km以下格距固定為100 m,3~15 km高度上則采用雙曲正切函數(shù)拉伸,最大格距為574 m,15 km以上格距都為574 m。這樣的分辨率已足以模擬出颮線的結(jié)構(gòu)及其組織化過程(Weisman and Rotunno, 2004)。模式頂層及底層為剛性邊界,水平方向上為開放邊界。為了減小模式頂層向下反射的影響,14 km高度以上開啟了瑞利阻尼選項(xiàng)(Takemi, 2006)。

    試驗(yàn)中溫濕背景場參考了Weisman and Klemp (1982)所設(shè)計(jì)的理想探空曲線(圖1),其計(jì)算公式如下:

    其中,對流層頂高度 ztr=12 km,對流層頂處的溫度Ttr=213 K,對流層頂處的位溫θtr=343 K;θ0表示地面位溫,恒定為300 K;g為重力加速度;cp為干空氣比定壓熱容;q(z)為基態(tài)位溫;r為相對濕度?;旌蠈拥母叨仍O(shè)置為1.2 km,在混合層中相對濕度最大值為0.95,依據(jù)RKW理想試驗(yàn)的最優(yōu)配置,我們將對照試驗(yàn)的最底層比濕設(shè)置為 qv=14 g kg?1,垂直切變在垂直方向上的范圍為0~2.5 km,水平風(fēng)速自地面至 2.5 km高度的變化值(Us)為17 m s?1,即地面風(fēng)速為?17 m s?1,線性增強(qiáng)至2.5 km處,以上風(fēng)速都為0 m s?1。

    對照試驗(yàn)初始對流啟動(dòng)采用的是熱泡激發(fā)方案,即給定一個(gè)中心位于模式底層x=90 km處(將此處設(shè)為零點(diǎn)),水平半徑為 10 km,垂直半徑為2.5 km的熱泡,初始擾動(dòng)位溫在中心點(diǎn)處設(shè)定為2 K,并按余弦函數(shù)向半徑范圍內(nèi)遞減,范圍以外則都為零。模擬采用了冰微物理過程(Lin et al., 1983)。此外,試驗(yàn)沒有考慮鋒、急流等颮線環(huán)境,也忽略了科氏力的影響。模式資料每 10分鐘輸出一次。

    2.2 颮線觸發(fā)試驗(yàn)結(jié)果

    理論研究(Thorpe et al., 1982; Rotunno et al., 1988)和敏感性試驗(yàn)(圖略)均表明:對照試驗(yàn)的這種環(huán)境配置條件,能夠使得對流觸發(fā)過程所能達(dá)到的最旺盛階段維持時(shí)間最長,強(qiáng)度最大。圖2給出了系統(tǒng)觸發(fā)過程不同發(fā)展階段的內(nèi)部環(huán)流結(jié)構(gòu)、擾動(dòng)位溫和擾動(dòng)比濕的垂直剖面圖。積分1.5 h時(shí)(圖 2a),颮線處于初步形成階段,在熱泡的中心點(diǎn)(x=0 km)處存在伸展高度接近6 km的垂直上升運(yùn)動(dòng)帶,對流層高層(8 km)處出現(xiàn)了以x=0 km為中心向正負(fù)兩個(gè)方向的輻散氣流;由于上升氣流不斷地將低層水汽輸送到中層,使得上升帶中擾動(dòng)比濕為正值,同時(shí)凝結(jié)作用產(chǎn)生的大量潛熱也使得上升帶3~5 km高度上出現(xiàn)明顯的擾動(dòng)增溫;上升氣流左側(cè)存在的明顯下沉運(yùn)動(dòng),一方面使得凝結(jié)的水汽下落形成降水,另一方面將高層較為干冷的空氣向下輸送,并在低層環(huán)境風(fēng)影響下向水平負(fù)方向擴(kuò)散,造成熱泡中心左側(cè)大片2 km高度以下的區(qū)域出現(xiàn)了負(fù)的擾動(dòng)比濕,同時(shí)1 km以下出現(xiàn)大面積擾動(dòng)降溫,說明對流啟動(dòng)后,在降水和干冷空氣下傳堆積的影響下,冷池開始形成。

    颮線進(jìn)一步發(fā)展后系統(tǒng)內(nèi)部結(jié)構(gòu)出現(xiàn)明顯變化(圖2b):首先由于近地面位溫以及比濕的擾動(dòng)降幅增大,冷池強(qiáng)度增強(qiáng),地面冷池前沿移動(dòng)至x=6 km處,且其后部水平風(fēng)速明顯減弱,表明冷池增強(qiáng)后干冷空氣質(zhì)量堆積產(chǎn)生的外流逐步削弱了近地面層基本環(huán)境氣流的影響;其次雖然冷池前沿的右移使得地面強(qiáng)迫上升的啟動(dòng)位置隨之右移,但對流層中層上升運(yùn)動(dòng)和潛熱釋放的大值中心與地面冷池前沿的水平位置相差不大,因此對流層中層上升運(yùn)動(dòng)區(qū)的局地增溫,以及高空輻散外流的抽氣和通風(fēng)效應(yīng),與地面冷池強(qiáng)迫垂直貫通,有利于垂直環(huán)流的進(jìn)一步增強(qiáng),同時(shí)這種垂直方向貫通的傾斜對流使得颮線系統(tǒng)上升運(yùn)動(dòng)和下沉運(yùn)動(dòng)分離,也有利于颮線系統(tǒng)的組織化和其后部降水增強(qiáng)。在颮線維持階段(圖 2c、2d),由于累積降水的不斷增加使得冷池強(qiáng)度進(jìn)一步加強(qiáng),在冷空氣堆質(zhì)量輻散的強(qiáng)大作用下,2 km高度以下擾動(dòng)位溫負(fù)值覆蓋區(qū)域水平風(fēng)速開始出現(xiàn)逆基本氣流的分布特征,因此地面冷池前沿及其引起的地面強(qiáng)迫上升的起始位置開始快速向右推進(jìn)。在此影響下,系統(tǒng)的環(huán)流結(jié)構(gòu)出現(xiàn)明顯的傾斜狀態(tài),對流層中層上升運(yùn)動(dòng)以及凝結(jié)和高空輻散造成擾動(dòng)增溫的大值中心均明顯落后于低層,這種垂直分布導(dǎo)致中高層凝結(jié)和輻散引起的對流運(yùn)動(dòng)與地面冷池造成的強(qiáng)迫抬升難以形成有效的連通。到積分6.0~6.5 h時(shí)(圖2e、2f),這一現(xiàn)象尤其明顯,因此圖2e、2f中雖然地面降溫、地面上升速度和中高層擾動(dòng)增溫相比前一階段(圖2c、2d)在量值上均增大,系統(tǒng)強(qiáng)度卻開始減弱。

    以上分析表明,颮線的發(fā)展伴隨著系統(tǒng)內(nèi)部垂直環(huán)流以及冷池強(qiáng)度的變化,冷池發(fā)展到一定強(qiáng)度后,不僅對颮線內(nèi)部結(jié)構(gòu)產(chǎn)生影響,還能夠通過與環(huán)境垂直風(fēng)切變等環(huán)境因素共同作用,對颮線的組織化和維持具有重要影響(RKW理論)。為了定量描述颮線不同階段冷池的強(qiáng)度、結(jié)構(gòu)和演變過程,本文對颮線發(fā)展過程中的冷池強(qiáng)度進(jìn)行定量計(jì)算。RKW理論給出的單點(diǎn)冷池強(qiáng)度(C)定義公式為

    其中,H為冷池的高度,B則表示浮力,具體定義如下:

    由圖2分析可知,冷池的形成和發(fā)展均出現(xiàn)在地面上升運(yùn)動(dòng)后部的大片區(qū)域中,在垂直方向上也存在伸展,單點(diǎn)強(qiáng)度很難準(zhǔn)確描述冷池的發(fā)展過程,需要對其垂直高度和水平覆蓋范圍進(jìn)行界定。

    圖1 三種比濕條件下的Skew-T圖。其中綠色點(diǎn)線、黑色實(shí)線、紅色虛線分別為12 g kg?1、14 g kg?1、16 g kg?1條件下的露點(diǎn)曲線,紅色實(shí)線為溫度曲線Fig. 1 Diagram of Skew-T under three specific humidity used in the numerical simulations. The green dotted line (qv= 12 g kg?1), black solid line (qv= 14 g kg?1), and red dashed line (qv= 16g kg?1) are dew point curves, the red solid line is temperature curve

    v在實(shí)際計(jì)算中,冷池的垂直高度H邊界取擾動(dòng)位溫自底層向上遞增至?1 K時(shí)的高度(Takemi, 2006);考慮到冷池的快速水平擴(kuò)散及其對颮線系統(tǒng)的影響,采用Takemi(2007)和Parker(2008)的方法,選取地面冷池前沿線至其后部20 km范圍作為冷池強(qiáng)度計(jì)算的水平區(qū)域。綜上討論,本文在計(jì)算過程中,取界定區(qū)域內(nèi)冷池強(qiáng)度的平均值作為系統(tǒng)的冷池強(qiáng)度。需要說明的是,本文設(shè)定的冷池前沿線為擾動(dòng)位溫梯度最大值處。

    圖2 對照試驗(yàn)不同發(fā)展階段的垂直剖面圖。陰影為擾動(dòng)位溫(單位:K),紅色等值線為擾動(dòng)比濕(間隔:1.0 g kg?1,未繪制零線),箭頭表示風(fēng)場(單位:m s?1)。橫坐標(biāo)為水平距離,零點(diǎn)為熱泡中心點(diǎn),縱坐標(biāo)為垂直距離Fig. 2 Vertical cross sections of control simulation in different stages of development. Color scales represent perturbation potential temperature (units: K), red contours represent perturbation specific humidity (interval: 1.0 g kg?1; zero line not plotted), and arrows represent wind field (units: m s?1). The x-axis is horizontal distance, the point of zero represent the center of thermal bubble, and the vertical axis is vertical distance

    為更具體的說明颮線強(qiáng)度與冷池強(qiáng)度之間的聯(lián)系,圖3給出了積分9 h內(nèi)5 km高度上的垂直運(yùn)動(dòng)、地面降溫以及冷切比。由于Us始終為17 m s?1,因此冷切比與冷池強(qiáng)度等價(jià)。由圖分析可知,積分初始0.7 h內(nèi),整個(gè)模擬區(qū)域在5 km高度上垂直速度均小于4 m s?1,說明此時(shí)對流強(qiáng)度較弱,地面擾動(dòng)位溫大于?1 K,冷池強(qiáng)度基本為零。積分0.7 h后,垂直運(yùn)動(dòng)增強(qiáng)并開始出現(xiàn)超過4 m s?1的上升運(yùn)動(dòng)大值中心,其位置基本位于x=0 km處,且大值中心在時(shí)間—空間的二維平面上呈現(xiàn)連續(xù)分布的狀態(tài),說明該區(qū)域開始出現(xiàn)強(qiáng)度較強(qiáng)且伸展高度較高的深厚對流,并且這類深厚對流在該區(qū)域能夠穩(wěn)定地激發(fā)和維持,組織化的颮線系統(tǒng)開始形成;同時(shí)地面也出現(xiàn)小于?1 K的擾動(dòng)位溫,冷池前沿的位置與系統(tǒng)上升運(yùn)動(dòng)大值中心位置一致,計(jì)算得到冷池強(qiáng)度在半小時(shí)內(nèi)快速增強(qiáng),冷切比達(dá)到0.65左右后穩(wěn)定持續(xù)至2.5 h。隨后,冷池強(qiáng)度進(jìn)入了近2個(gè)小時(shí)的持續(xù)快速增長期,到積分4.5 h時(shí)冷切比超過 1.5,此時(shí)上升運(yùn)動(dòng)大值區(qū)不僅在時(shí)間和強(qiáng)度上保持連續(xù)發(fā)展的態(tài)勢,水平范圍也明顯擴(kuò)大;地面擾動(dòng)位溫減小進(jìn)一步加劇,冷池前沿線開始出現(xiàn)較為明顯的向前移動(dòng),5 km高度上升運(yùn)動(dòng)中心稍落后于冷池前沿,但兩者水平距離較小。隨后到積分6 h時(shí)間段內(nèi),雖然冷切比仍呈現(xiàn)較為平緩的增大趨勢,但是上升運(yùn)動(dòng)大值中心的分布呈現(xiàn)強(qiáng)度減弱,結(jié)構(gòu)不連續(xù)的特征,說明颮線強(qiáng)度開始減弱。同時(shí)冷池前沿隨時(shí)間快速向x正方向推進(jìn),且逐漸遠(yuǎn)離上升運(yùn)動(dòng)中心;到積分6 h后,冷切比維持在 2.2左右,地面出現(xiàn)超過?13 K的降溫中心,但5 km高度上的上升運(yùn)動(dòng)結(jié)構(gòu)零散,強(qiáng)度迅速減弱直至消失。Weisman和Rotunno(2004)指出當(dāng)冷切比介于1.0~1.5時(shí),最適宜颮線的發(fā)展,從上述分析來看本文結(jié)果與其一致。

    由此可見,在RKW理論最優(yōu)配置條件下,颮線系統(tǒng)經(jīng)歷了從對流啟動(dòng)到深對流組織化,最后減弱消散的演變過程。在這一過程中,冷池強(qiáng)度的變化與颮線強(qiáng)度有很好的對應(yīng)關(guān)系,這與RKW理論結(jié)果一致,與此同時(shí)系統(tǒng)內(nèi)部位溫、比濕以及風(fēng)場發(fā)生了顯著的變化,即動(dòng)量、熱量以及水汽出現(xiàn)區(qū)域性的再分配,與之相對應(yīng)的則是系統(tǒng)內(nèi)部環(huán)流結(jié)構(gòu)發(fā)生改變,并依此對颮線的組織化程度產(chǎn)生影響。接下來,我們將在颮線發(fā)展過程中改變低層環(huán)境垂直風(fēng)切變或水汽,探討這種情形下系統(tǒng)內(nèi)部各物理因子的響應(yīng)過程及其對颮線系統(tǒng)強(qiáng)度變化的影響。

    3 低層環(huán)境垂直風(fēng)切變的變化對颮線強(qiáng)度的影響

    從對照試驗(yàn)可以看出,對流發(fā)展及其伴隨的降水過程會使得冷池強(qiáng)度持續(xù)增強(qiáng),與此同時(shí),整體環(huán)境風(fēng)的垂直切變變化卻并不明顯,即冷切比不斷增大。因此,依據(jù)RKW理論及Weisman and Rotunno (2004)結(jié)果,冷切比介于1~1.5時(shí)颮線強(qiáng)度最強(qiáng),隨著冷切比的持續(xù)增加,颮線強(qiáng)度開始減弱。而環(huán)境垂直風(fēng)切變是影響冷切比的關(guān)鍵因子,因此本節(jié)將低層環(huán)境垂直風(fēng)切變作為敏感性因子,探討在切變改變的條件下颮線強(qiáng)度如何變化,并從颮線系統(tǒng)的內(nèi)部環(huán)流結(jié)構(gòu)出發(fā)解釋這一強(qiáng)度變化原因。

    對照試驗(yàn)結(jié)果表明2.5~4.5 h這一階段是颮線發(fā)展最強(qiáng)盛的階段,隨后垂直運(yùn)動(dòng)減弱,因此我們選取了對照試驗(yàn)積分3 h這一積分時(shí)刻的颮線狀態(tài)作為敏感性試驗(yàn)的初始對流,之后颮線所處環(huán)境相較對照試驗(yàn)發(fā)生變化(指垂直風(fēng)切變或水汽),其他配置則與對照試驗(yàn)相同。環(huán)境垂直風(fēng)切變只考慮低層切變一種類型(0~2.5 km),但設(shè)計(jì)了三種不同大小的Us值,即14 m s?1,17 m s?1,20 m s?1。各試驗(yàn)命名見表1。

    表1 各試驗(yàn)命名。其中,第一個(gè)字符L表示低層切變,第二個(gè)字符L、M、H分別表示低層比濕較低(qv=12 g kg?1)、適中(qv=14 g kg?1)、較高(qv=16g kg?1)Table 1 Details of the experiments, in which the first character L represents the low level shear, the second characters L, M, and H represent the lower (qv=12 g kg?1), moderate (qv=14 g kg?1), and higher (qv=16 g kg?1) specific humidity, respectively

    3.1 颮線強(qiáng)度對低層環(huán)境垂直風(fēng)切變的響應(yīng)

    圖3 對照試驗(yàn)z=5 km高度上垂直速度、地面擾動(dòng)位溫的時(shí)空分布以及冷池強(qiáng)度隨時(shí)間變化。左圖中,暖色陰影表示垂直速度(只給出速度大于4 m s?1的值),冷色陰影為地面擾動(dòng)位溫分布(只給出小于?1 K的值),黑色實(shí)線表示冷池前沿線位置。右圖中,紅色實(shí)線表示冷池強(qiáng)度隨時(shí)間變化Fig. 3 Depiction of the control simulation: (left) temporal–spatial diagram for vertical velocity (warm colors; values less than ?4 m s?1not plotted) at 5 km and the surface perturbation potential temperature (cold colors; values greater than ?1 K not plotted), in which the black solid line represents the leading edge of the surface cold pool; (right) the strength of the cold pool

    圖4 低層切變條件下各試驗(yàn)的(a)最大垂直速度(wmax)的時(shí)間(3~7 h)平均以及(b)3~7 h時(shí)間內(nèi)大于1 m s?1的垂直速度所覆蓋區(qū)域(w-area)的平均值隨高度的分布Fig. 4 Temporally averaged (3–7 h) vertical profiles of the (a) maximum updraft speed (wmax) and (b) fractional area of the updraft speed (w-area) more than 1 m s?1

    颮線強(qiáng)度用最大垂直速度以及超過1 m s?1的垂直速度所覆蓋區(qū)域來表示(Takemi, 2014)。圖4給出了最大垂直速度和垂直上升運(yùn)動(dòng)所覆蓋區(qū)域在3~7 h內(nèi)時(shí)間平均(時(shí)間t包含了對照試驗(yàn)的前3小時(shí)積分時(shí)間)的垂直分布。總的來看,最大垂直速度在1 km高度以上開始顯著增大,而在9 km以上則明顯減小,峰值集中在5~8 km。在混合層比濕相同的情況下,Us值愈大則峰值愈大,如LM-20試驗(yàn)的峰值明顯大于 LM-17試驗(yàn),LM-17試驗(yàn)又明顯大于LM-14試驗(yàn)。說明增大低層環(huán)境垂直風(fēng)切變有利于颮線強(qiáng)度增強(qiáng),反之亦然。各試驗(yàn)的垂直上升運(yùn)動(dòng)所覆蓋區(qū)域在中低層無明顯差異,高層卻明顯不同。值得注意的是,各試驗(yàn)高層垂直上升運(yùn)動(dòng)覆蓋區(qū)域的大小順序與其中層最大垂直速度的大小順序相同,即用這兩種方式表示颮線的強(qiáng)度具有一致性。

    為了反應(yīng)颮線強(qiáng)度的變化趨勢,計(jì)算了各小時(shí)內(nèi)最大垂直速度均值及其標(biāo)準(zhǔn)偏差(圖5)。同樣在水汽適中的條件下,增大低層環(huán)境垂直風(fēng)切變(LM-20)所對應(yīng)的各時(shí)間段的均值較對照試驗(yàn)或減小低層環(huán)境垂直風(fēng)切變(LM-14)更大,這與上述的結(jié)論是一致的。另外,從均值隨時(shí)間的變化來看,LM-17試驗(yàn)與LM-20試驗(yàn)的均值先增強(qiáng)后減弱,而LM-14試驗(yàn)則單調(diào)遞減。即減小低層環(huán)境垂直風(fēng)切變會明顯減弱颮線的強(qiáng)度。

    圖5 低層切變條件下不同試驗(yàn)各小時(shí)內(nèi)最大垂直速度均值(不同記號表示不同時(shí)間段內(nèi)的均值)及其標(biāo)準(zhǔn)偏差(虛線)。橫坐標(biāo)為各個(gè)試驗(yàn)名稱(C表示對照試驗(yàn)),縱坐標(biāo)則為最大垂直速度均值大?。▎挝唬簃 s?1)Fig. 5 The means (symbols) and standard deviations (dashed lines) of the maximum updraft speed during integrating 2–3 h (filled circle), 3–4 h (asterisks), 4–5 h (circles), 5–6 h (squares), and 6–7 h (triangles) for all experiments. The x-axis is the name of experiments (the character C represents the control experiment). The vertical axis is the mean of maximum vertical

    3.2 颮線強(qiáng)度對環(huán)境垂直風(fēng)切變的響應(yīng)分析

    圖6、圖7分別為減?。↙M-14)和增加(LM-20)低層環(huán)境垂直風(fēng)切變試驗(yàn)不同發(fā)展時(shí)刻的擾動(dòng)位溫、擾動(dòng)比濕、風(fēng)場的垂直剖面圖。減小低層切變后,無論是對流高度還是擾動(dòng)比濕大小(負(fù)值區(qū),可大致反應(yīng)降水強(qiáng)度)相較對照試驗(yàn)(圖 2)都明顯減弱,這與3.1節(jié)中結(jié)論一致。從移動(dòng)速度來看,t=1 h時(shí)對照試驗(yàn)的冷池前沿位于x=12 km附近,而LM-14試驗(yàn)已達(dá)到x=18 km處,即環(huán)境垂直風(fēng)切變的減小使得颮線移動(dòng)速度加快。從圖8a中LM-14,LM-17,LM-20三個(gè)試驗(yàn)的低層平均水平速度(水平方向上為冷池前沿線至其后側(cè)20 km,垂直方向?yàn)?~1 km共10層的水平速度平均)也可以看出,環(huán)境垂直風(fēng)切變愈大,平均水平風(fēng)速愈小。Wakimoto(1982)以及Mueller and Carbon(1987)的觀測分析表明對于大多數(shù)對流風(fēng)暴來說其傳播速度由風(fēng)暴產(chǎn)生的冷池的移動(dòng)速度決定,因?yàn)樾碌娘L(fēng)暴單體極易被冷池前沿的陣風(fēng)鋒觸發(fā)從而維持整個(gè)風(fēng)暴系統(tǒng)例如颮線的發(fā)展維持。在圖 6a中上升運(yùn)動(dòng)存在兩個(gè)中心,一個(gè)位于冷池前沿(x=18 km),一個(gè)位于x=12 km處。即在開始階段,颮線移速的加快使得冷池前沿陣風(fēng)鋒觸發(fā)的對流與中層凝結(jié)潛熱釋放激發(fā)的上升運(yùn)動(dòng)出現(xiàn)分離,t=5.0 h時(shí)這種分離不再明顯,而隨著颮線的發(fā)展,這一分離再次顯現(xiàn)(圖6c、6d),且在t=6.5 h時(shí)出現(xiàn)明顯分離。增加低層切變后(圖7),颮線移速減小,上升運(yùn)動(dòng)并沒有出現(xiàn)如對照試驗(yàn)所出現(xiàn)的不連續(xù)現(xiàn)象,且無論是從擾動(dòng)比濕(負(fù)值區(qū))還是對流高度來看,都較對照試驗(yàn)明顯增強(qiáng)。依據(jù)第二部分中冷池的計(jì)算公式,給出了這三組試驗(yàn)的冷切比隨時(shí)間的變化曲線(圖 8b),可以發(fā)現(xiàn)低層切變越小,冷切比越大,即在切變發(fā)生改變的情況下,RKW理論定量的解釋了颮線強(qiáng)度變化的原因。而通過分析環(huán)流結(jié)構(gòu)的變化,則具體的展現(xiàn)了颮線強(qiáng)度變化的內(nèi)部機(jī)理。

    以上分析表明,颮線前側(cè)濕空氣被抬升至中高層,形成颮線系統(tǒng)的上升氣流,濕空氣抬升凝結(jié)釋放潛熱使得氣流進(jìn)一步上升,濕空氣凝結(jié)形成雨滴降至地面,降水同時(shí)拖拽空氣下沉;雨水的蒸發(fā)使得近地面溫度降低,形成中尺度高壓區(qū)域。這兩者共同作用下在冷池區(qū)域形成輻散氣流,向x正方向的輻散氣流(陣風(fēng)鋒)與入流形成輻合從而觸發(fā)新的對流產(chǎn)生,并使得颮線向前推進(jìn)。低層環(huán)境垂直風(fēng)切變的減小使得輻合中心向右側(cè)推進(jìn)速度加快,從而使得冷池前沿陣風(fēng)鋒觸發(fā)的新的對流與中層凝結(jié)潛熱釋放激發(fā)的上升運(yùn)動(dòng)分離,颮線強(qiáng)度減弱,反之亦然。

    圖6 同圖2,但為LM-14試驗(yàn)Fig. 6 As in Fig. 2, but for experiment LM-14

    4 低層水汽變化對颮線強(qiáng)度的影響

    Meng et al.(2013)通過統(tǒng)計(jì)中國東部的颮線發(fā)現(xiàn),在同樣的不穩(wěn)定條件下不同程度的比濕影響著颮線的發(fā)生程度,即比濕的大小會影響颮線的發(fā)生發(fā)展。因此本章將低層比濕作為敏感性因子,探討在比濕改變的條件下颮線強(qiáng)度如何變化,并從颮線系統(tǒng)的內(nèi)部環(huán)流結(jié)構(gòu)出發(fā)解釋這一強(qiáng)度變化原因。

    與垂直風(fēng)切變的敏感性試驗(yàn)類似,在這組敏感性試驗(yàn)中將敏感性因子改為低層比濕,而垂直風(fēng)切變則與對照試驗(yàn)一致。試驗(yàn)設(shè)計(jì)了三種不同的混合層比濕:12 g kg?1、14 g kg?1、16 g kg?1,代表低層的水汽含量較低(L)、適中(M)、及較高(H)。不同的混合層比濕對應(yīng)的環(huán)境CAPE也不相同,分別為1144 J kg?1、2219 J kg?1、3500 J kg?1。就颮線的發(fā)展演變來看,CAPE的變化范圍很大,一般而言,CAPE在1000 J kg?1左右時(shí),對形成一般的颮線系統(tǒng)比較有利,而當(dāng)其超過2200 J kg?1時(shí),則有可能形成較強(qiáng)的颮線系統(tǒng)(Bluestein et al., 1987;Parker and Johnson, 2000;Weckwerth, 2000)。各試驗(yàn)命名見表1。

    4.1 颮線強(qiáng)度對低層水汽的響應(yīng)

    圖4也可反應(yīng)颮線強(qiáng)度對低層水汽的響應(yīng)。從圖中可以看出,在Us相同的情況下,低層比濕愈大則峰值愈大,如LH-17試驗(yàn)的峰值明顯大于LM-17試驗(yàn),而LM-17試驗(yàn)又明顯大于LL-17試驗(yàn)。圖4b所呈現(xiàn)的圖像反應(yīng)出了同樣的結(jié)果。從颮線強(qiáng)度的變化趨勢圖(圖 5)也同樣可以看出,隨著低層比濕的增加,各時(shí)間段的均值會隨之增大。因此以得出這樣的結(jié)論:低層水汽供應(yīng)越充分,越有利于颮線強(qiáng)度的增強(qiáng)。孫建華等(2014)對實(shí)際個(gè)例的水汽敏感性試驗(yàn)也印證了這一結(jié)論。

    圖7 同圖2,但為LM-20試驗(yàn)Fig. 7 As in Fig. 2, but for experiment LM-20

    圖8 LM-14(紅線)、LM-17(黑線)、LM-20(藍(lán)線)三個(gè)試驗(yàn)(a)平均水平速度和(b)冷池強(qiáng)度與低層環(huán)境垂直風(fēng)切變的比值(C/Us)隨時(shí)間變化Fig. 8 Time series of (a) average horizontal velocity and (b) the ratio between the strength of cold pool and low-level VWS (C/Us) for experiments LM-14 (red line), LM-17 (black line), and LM-20 (blue line)

    4.2 颮線強(qiáng)度對低層水汽的響應(yīng)分析

    圖9 同圖2,但為LL-17試驗(yàn)Fig. 9 As in Fig. 2, but for experiment LL-17

    圖9、圖10分別為減少(LL-17)和增加(LH-17)低層比濕試驗(yàn)不同發(fā)展時(shí)刻的擾動(dòng)位溫、擾動(dòng)比濕、風(fēng)場的垂直剖面圖。與對照試驗(yàn)(圖2)相比較可以發(fā)現(xiàn),三者在颮線推進(jìn)速度上存在差異,但并不明顯,即改變低層比濕后,颮線強(qiáng)度變化的原因與改變環(huán)境垂直風(fēng)切變的有所不同。而在垂直上升運(yùn)動(dòng)、擾動(dòng)位溫強(qiáng)度(尤其是正值),擾動(dòng)比濕強(qiáng)度上,三者差異明顯。LL-17試驗(yàn)最弱,LH-17最強(qiáng),LM-17居中。即低層比濕愈大,颮線強(qiáng)度愈強(qiáng),這與4.1的結(jié)論一致。從擾動(dòng)比濕(正值)來看,LH-17試驗(yàn)較LL-17試驗(yàn)擾動(dòng)幅度更大,即在低層增加比濕的條件下,輸送到中層的水汽更多,相應(yīng)的凝結(jié)潛熱釋放也更多(從擾動(dòng)位溫正值可看出);比濕的增加使得颮線降水強(qiáng)度也相應(yīng)增加,這從冷池前沿線后方的擾動(dòng)比濕負(fù)值大小可以看出。降水愈強(qiáng),說明空氣中的水汽凝結(jié)成液態(tài)水愈多,比濕變化也就越大。減小低層比濕的情況下(圖9),中層釋放的凝結(jié)潛熱相對較少,從而使得颮線發(fā)展后期中層不再出現(xiàn)明顯的上升運(yùn)動(dòng)或上升運(yùn)動(dòng)較弱(圖9c、9d),而增加低層比濕則與此相反;低層比濕的減小使得冷池前沿陣風(fēng)鋒觸發(fā)的對流強(qiáng)度減小,這與CAPE的釋放大小有關(guān)。

    為了反映各試驗(yàn)在颮線移至?xí)r環(huán)境CAPE的釋放大小,我們對x=0 km至x=30 km范圍內(nèi)的CAPE進(jìn)行了區(qū)域平均,并給出了其隨時(shí)間的變化曲線(圖11)。從圖中可以看出,所有試驗(yàn)的平均CAPE隨時(shí)間基本呈現(xiàn)出單調(diào)遞減的趨勢,并最終都趨近于零,這是颮線過境后,環(huán)境CAPE釋放的結(jié)果。而在不同的條件下,平均CAPE的初始值及變化曲線存在顯著差異,在低層切變相同的情況下,低層比濕越大,則平均CAPE的變化幅度越大,即CAPE釋放值越大。對比圖中LH-17與LL-17試驗(yàn)的變化曲線可以發(fā)現(xiàn),前者較后者所釋放的CAPE更大,這決定了各自在冷池前沿新激發(fā)的對流強(qiáng)度LL-17試驗(yàn)要顯著弱于LH-17試驗(yàn)。從圖中還可以看出,低層比濕相同的情況下,切變越大,則平均CAPE從初始值遞減至零所需時(shí)間越長,即颮線移動(dòng)的速度越慢。而在低層切變一致時(shí),無論低層比濕條件如何,遞減至零所需時(shí)間基本相同。即颮線的移動(dòng)速度主要與低層切變有關(guān),低層比濕則影響不大。

    以上分析表明,在保持低層環(huán)境垂直風(fēng)切變不變的情況下,改變低層比濕條件,颮線的移動(dòng)速度未發(fā)生明顯變化,而輸送至颮線系統(tǒng)中高層的水汽卻存在不同,這使得水汽凝結(jié)釋放潛熱量不同,進(jìn)而使得中高層的垂直上升運(yùn)動(dòng)強(qiáng)度出現(xiàn)差異;此外,低層比濕的改變也使得環(huán)境CAPE發(fā)生變化,從而影響新激發(fā)的對流強(qiáng)度。即減?。ㄔ黾樱┑蛯颖葷褚环矫鏁沟美涑厍把氐年囷L(fēng)鋒觸發(fā)的新的對流強(qiáng)度減弱(增強(qiáng)),另一方面輸送到中層的水汽減少(增加),相應(yīng)的凝結(jié)潛熱釋放減小(增大),進(jìn)而使得中層的上升運(yùn)動(dòng)減弱(增強(qiáng)),兩處上升運(yùn)動(dòng)的減小及分離(增強(qiáng)并相互貫通)使得颮線強(qiáng)度減弱(增強(qiáng))。

    圖10 同圖2,但為LH-17試驗(yàn)Fig. 10 As in Fig. 2, but for experiment LH-20

    圖11 各試驗(yàn)區(qū)域平均(x=0~30 km)的對流有效位能(CAPE)值隨時(shí)間變化Fig. 11 Time series of regional average (x = 0–30 km) convective available potential energy (CAPE) for all experiments

    圖12 環(huán)境參數(shù)影響下颮線內(nèi)部結(jié)構(gòu)調(diào)整及其對颮線強(qiáng)度影響的示意圖Fig. 12 Schematic diagram of squall line internal restructuring and variation in intensity for the effects of environmental parameters

    5 結(jié)論與討論

    本文利用ARPS模式對颮線發(fā)生發(fā)展過程進(jìn)行了二維理想數(shù)值試驗(yàn),并從颮線系統(tǒng)內(nèi)部動(dòng)量、熱量、水汽等物理因子配置的變化過程入手,分析了颮線系統(tǒng)組織化對流的形成、發(fā)展和消亡機(jī)制,討論了當(dāng)?shù)蛯迎h(huán)境垂直風(fēng)切變和水汽等環(huán)境要素場發(fā)生變化時(shí),系統(tǒng)內(nèi)部各物理因子的響應(yīng)過程及其對颮線系統(tǒng)強(qiáng)度變化的影響。

    研究表明,在深對流演變過程中,颮線系統(tǒng)內(nèi)部出現(xiàn)動(dòng)量、熱量以及水汽的區(qū)域性再分配,導(dǎo)致系統(tǒng)內(nèi)環(huán)流結(jié)構(gòu)、冷池強(qiáng)度發(fā)生改變,并依此對颮線的組織化程度產(chǎn)生影響。通過颮線強(qiáng)度對低層環(huán)境垂直風(fēng)切變和水汽的敏感性試驗(yàn)發(fā)現(xiàn),環(huán)境因素的變化同樣會使颮線系統(tǒng)內(nèi)部出現(xiàn)動(dòng)量、熱量、水汽的再分配過程,并依此對颮線強(qiáng)度產(chǎn)生影響。改變低層環(huán)境垂直風(fēng)切變時(shí),颮線系統(tǒng)內(nèi)部動(dòng)量、熱量以及水汽的再分配使得颮線移速發(fā)生變化,從而影響颮線強(qiáng)度:當(dāng)增大(減?。┑蛯迎h(huán)境垂直風(fēng)切變時(shí),颮線后側(cè)低層的水平速度減?。ㄔ黾樱?,從而使得颮線移速減慢(加快),新激發(fā)的對流與中高層的垂直運(yùn)動(dòng)相互貫通(分離),颮線強(qiáng)度得以增強(qiáng)(減弱);改變低層水汽時(shí),颮線系統(tǒng)內(nèi)部動(dòng)量、熱量以及水汽的再分配則主要影響中高層的垂直運(yùn)動(dòng)以及新觸發(fā)的對流,從而使得颮線強(qiáng)度發(fā)生變化:增加(減少)低層水汽時(shí),輸送到中高層的水汽增加(減少),水汽凝結(jié)釋放潛熱得以增加(減少),這使得中高層的垂直運(yùn)動(dòng)增強(qiáng)(減弱);同時(shí),增加(減少)低層水汽使得颮線移至?xí)r環(huán)境所釋放的 CAPE增大(減?。涑厍把厮ぐl(fā)的對流增強(qiáng)(減弱)。兩者作用下新激發(fā)的對流與中高層的垂直運(yùn)動(dòng)得以相互貫通(分離),颮線強(qiáng)度增強(qiáng)(減弱)。颮線強(qiáng)度對低層切變及水汽的響應(yīng)可由示意圖圖12來表示。

    本文對颮線發(fā)展過程中其強(qiáng)度對低層環(huán)境垂直風(fēng)切變和水汽的響應(yīng)及其原因進(jìn)行了詳細(xì)的研究,但文章涉及的數(shù)值試驗(yàn)皆為兩維結(jié)構(gòu),這與實(shí)際颮線的三維結(jié)構(gòu)存在差異,因此后續(xù)的研究所有試驗(yàn)將擴(kuò)展到三維。此外,理想試驗(yàn)只代表理想環(huán)境配置條件下颮線的發(fā)展變化過程,下一步,將利用實(shí)際颮線個(gè)例的觀測和數(shù)值模擬結(jié)果進(jìn)行診斷分析。

    參考文獻(xiàn)(References)

    Bluestein H B, Jain M H. 1985. Formation of mesoscale lines of precipitation: Severe squall lines in Oklahoma during the spring [J]. J. Atmos. Sci., 42 (16): 1711–1732.

    Bluestein H B, Marx G T, Jain M H. 1987. Formation of mesoscale lines of precipitation: Non-Severe squall lines in Oklahoma during the spring [J]. Mon. Wea. Rev., 115 (11): 2719–2727.

    陳明軒, 王迎春. 2012. 低層垂直風(fēng)切變和冷池相互作用影響華北地區(qū)一次颮線過程發(fā)展維持的數(shù)值模擬 [J]. 氣象學(xué)報(bào), 70 (3): 371–386. Chen Mingxuan, Wang Yingchun. 2012. Numerical simulation study of interactional effects of the low-level vertical wind shear with the cold pool on a squall line evolution in North China [J]. Acta Meteorologica Sinica (in Chinese), 70 (3): 371–386.

    Coniglio M C, Corfidi S F, Kain J S. 2012. Views on applying RKW theory: An illustration using the 8 May 2009 derecho-producing convective system [J]. Mon. Wea. Rev., 140 (3): 1023–1043.

    Fovell R G, Ogura Y. 1988. Numerical simulation of a mid-latitude squall line in two dimensions [J]. J. Atmos. Sci., 45 (24): 3846–3879.

    Fujita T. 1955. Results of detailed synoptic studies of squall lines [J]. Tellus, 7 (4): 405–436.

    Houze Jr R A, Biggerstaff M I, Rutledge S A, et al. 1989. Interpretation of Doppler weather radar displays of midlatitude mesoscale convective systems [J]. Bull. Amer. Meteor. Soc., 70 (6): 608–619.

    Lin Y L, Farley R D, Orville H D. 1983. Bulk parameterization of the snow field in a cloud model [J]. J. Climate Appl. Meteor., 22 (6): 1065–1092.

    Meng Z Y, Yan D C, Zhang Y J. 2013. General features of squall lines in East China [J]. Mon. Wea. Rev., 141 (5): 1629–1647.

    Mueller C K, Carbone R E. 1987. Dynamics of a thunderstorm outflow [J].J Atmos. Sci., 44 (15): 1879–1898.

    Newton C W. 1950. Structure and mechanism of the prefrontal squall line [J]. J. Meteor., 7 (3): 210–222.

    Parker D J. 1998. The dependence of cold-pool depth on source conditions [J]. Mon. Wea. Rev., 126 (2): 516–520.

    Parker M D. 2008. Response of simulated squall lines to low-level cooling [J]. J. Atmos. Sci., 65 (4): 1323–1341.

    Parker M D, Johnson R H. 2000. Organizational modes of midlatitude mesoscale convective systems [J]. Mon. Wea. Rev., 128 (10): 3413–3436.

    Parker M D, Johnson R H. 2004. Structures and dynamics of quasi-2D mesoscale convective systems [J]. J. Atmos. Sci., 61 (5): 545–567.

    Rotunno R, Klemp J B, Weisman M L. 1988. A theory for strong, long-lived squall lines [J]. J. Atmos. Sci., 45 (3): 463–485.

    孫建華, 鄭淋淋, 趙思雄. 2014. 水汽含量對颮線組織結(jié)構(gòu)和強(qiáng)度影響的數(shù)值試驗(yàn) [J]. 大氣科學(xué), 38 (4): 742–755. Sun Jianhua, Zheng Linlin, Zhao Sixiong. 2014. Impact of moisture on the organizational mode and intensity of squall lines determined through numerical experiments [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 38 (4): 742?755.

    Takemi T. 2006. Impacts of moisture profile on the evolution and organization of midlatitude squall lines under various shear conditions [J]. Atmospheric Research, 82 (1–2): 37–54.

    Takemi T. 2007. A sensitivity of squall-line intensity to environmental static stability under various shear and moisture conditions [J]. Atmospheric Research, 84 (4): 374–389.

    Takemi T. 2014. Convection and precipitation under various stability and shear conditions: Squall lines in tropical versus midlatitude environment [J]. Atmospheric Research, 142: 111–123.

    Thorpe A J, Miller M J, Moncrieff M W. 1982. Two-dimensional convection in non-constant shear: A model of mid-latitude squall lines [J]. Quart. J. Roy. Meteor. Soc., 108 (458): 739–762.

    Wakimoto R M. 1982. The life cycle of thunderstorm gust fronts as viewed with Doppler radar and rawinsonde data [J]. Mon. Wea. Rev., 110 (8): 1060–1082.

    Weckwerth T M. 2000. The effect of small-scale moisture variability on thunderstorm initiation [J]. Mon. Wea. Rev., 128 (12): 4017–4030.

    Weisman M L, Klemp J B. 1982. The dependence of numerically simulated convective storms on vertical wind shear and buoyancy [J]. Mon. Wea. Rev., 110 (6): 504–520.

    Weisman M L, Rotunno R. 2004. “A theory for strong long-lived squall lines” revisited [J]. J. Atmos. Sci., 61 (4): 361–382.

    Wyss J, Emanuel K A. 1988. The pre-storm environment of midlatitude prefrontal squall lines [J]. Mon. Wea. Rev., 116 (3): 790–794.

    Xue M, Droegemeier K K, Wong V. 2000. The Advanced Regional Prediction System (ARPS)—A multi-scale nonhydrostatic atmospheric simulation and prediction model. Part I: Model dynamics and verification [J]. Meteor. Atmos. Phys., 75 (3–4): 161–193.

    資助項(xiàng)目 國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目2013CB430103,國家自然科學(xué)基金項(xiàng)目41275002、41230421,江蘇省高校自然科學(xué)研究重大項(xiàng)目14KJA170005,江蘇省“333高層次人才培養(yǎng)工程”項(xiàng)目,災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室開放課題(2014LASW-B08)

    Funded by National Key Basic Research and Development Project of China (973 Program) (Grant 2013CB430103), National Natural Science Foundation of China (Grants 41275002 and 41230421), The Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant 14KJA170005), “333 High-Level Talent Development Program” in Jiangsu Province, National Key Laboratory of Disaster Weather (Grant 2014LASW-B08)

    文章編號1006-9895(2016)04-0689-14 中圖分類號 P446

    文獻(xiàn)標(biāo)識碼A

    doi:10.3878/j.issn.1006-9895.1505.14337

    收稿日期2014-12-12;網(wǎng)絡(luò)預(yù)出版日期 2015-07-07

    作者簡介張建軍,男,1989年出生,碩士研究生,主要從事中尺度研究。E-mail: zhangjj_nuist@126.com

    通訊作者鐘瑋,E-mail: wzhong_vivian@126.com

    Impact of Vertical Wind Shear and Moisture on the Organization of Squall Lines

    ZHANG Jianjun1, 2, WANG Yongqing2, and ZHONG Wei3, 4
    1 Hunan Climate Center, Changsha 410118
    2 Pacific Typhoon Research Center, School of Atmospheric Sciences, Nanjing University of Information Science &amp; Technology, Nanjing 210044
    3 Institute of Meteorology and Ocean, PLA University of Science and Technology, Nanjing 211101
    4 Key Laboratory of Mesoscale Severe Weather (Nanjing University), Ministry of Education, Nanjing 210093

    AbstractThe redistribution of physical factors and its impacts on the intensity of squall lines under the influence oflow-level Vertical Wind Shear (VWS) and moisture content are examined through two-dimensional idealized simulations with the ARPS model (the University of Oklahoma’s Advanced Research Prediction System). It shows that the redistribution of momentum, heat and moisture during the evolution of squall lines leads to the change of inner vertical circulation and the configuration of perturbation temperature and humidity, which affects the organization of deep convection and the intensity of the system. The results of sensitivity tests of low-level VWS and moisture content show that increasing (decreasing) the low-level VWS decelerates (accelerates) the propagation of the squall line, and makes the connection (separation) between the mid-level upward current and the new forced updrafts at the front edge of the cold pool, which corresponds to the intensification (weakening) of the squall line. On the other hand, increasing (decreasing) the low-level moisture content results in an increase (decrease) of moisture delivery from the low to middle level, which enhances (weakens) the mid-level latent heating and upward movement. Energy analysis indicates that the low-level moisture change influences the release of Convective Available Potential Energy (CAPE), and the intensity of the new convection. The combined effects of latent heating and CAPE released from low-level moisture change also affect the squall line intensity through exerting an influence on the organization of the upper-level upward currents and the new forced updrafts at the front edge of the cold pool.

    KeywordsSquall line, Numerical simulation, Cold pool, Vertical Wind Shear, Moisture

    猜你喜歡
    數(shù)值模擬
    基于AMI的雙色注射成型模擬分析
    錐齒輪精密冷擺輾成形在“材料成型數(shù)值模擬”課程教學(xué)中的應(yīng)用
    基于氣象信息及風(fēng)場信息的風(fēng)機(jī)輪轂處風(fēng)速預(yù)測
    鉆孔灌注樁樁底沉渣對樁體承載特性影響的模擬分析
    西南地區(qū)氣象資料測試、預(yù)處理和加工研究報(bào)告
    科技資訊(2016年18期)2016-11-15 08:01:18
    張家灣煤礦巷道無支護(hù)條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護(hù)參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機(jī)翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機(jī)的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    欧美日本中文国产一区发布| 美女大奶头黄色视频| 一二三四在线观看免费中文在| 97人妻天天添夜夜摸| 亚洲精品国产av成人精品| 国产精品99久久99久久久不卡 | 日韩免费高清中文字幕av| 成年人午夜在线观看视频| 国产精品久久久久久精品电影小说| 日韩中文字幕欧美一区二区 | 少妇被粗大猛烈的视频| 香蕉丝袜av| 成年美女黄网站色视频大全免费| 1024视频免费在线观看| 亚洲精品,欧美精品| 一区在线观看完整版| 日本午夜av视频| 老鸭窝网址在线观看| 午夜福利在线免费观看网站| 色网站视频免费| 老司机影院毛片| 看十八女毛片水多多多| 日本-黄色视频高清免费观看| a级毛片在线看网站| 成人黄色视频免费在线看| 婷婷色综合大香蕉| 99热网站在线观看| 午夜久久久在线观看| 久久人人爽av亚洲精品天堂| 天堂8中文在线网| 精品卡一卡二卡四卡免费| 亚洲四区av| 亚洲三区欧美一区| 久久精品国产鲁丝片午夜精品| 热99久久久久精品小说推荐| 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 亚洲欧洲国产日韩| 国产极品天堂在线| 亚洲经典国产精华液单| 大香蕉久久网| 香蕉国产在线看| 免费av中文字幕在线| 老司机亚洲免费影院| 色视频在线一区二区三区| av天堂久久9| 男人爽女人下面视频在线观看| 久久人人97超碰香蕉20202| 国产精品国产av在线观看| 欧美日韩一区二区视频在线观看视频在线| 99国产综合亚洲精品| 亚洲成色77777| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲色图综合在线观看| 国产乱来视频区| 国产成人精品一,二区| 高清视频免费观看一区二区| 亚洲精品aⅴ在线观看| 精品国产一区二区久久| 国产精品女同一区二区软件| 亚洲综合精品二区| a 毛片基地| 狠狠精品人妻久久久久久综合| 日韩电影二区| 精品酒店卫生间| 伊人亚洲综合成人网| 免费在线观看黄色视频的| 国产亚洲最大av| 一级毛片我不卡| 如何舔出高潮| 一个人免费看片子| 亚洲第一青青草原| 亚洲,欧美,日韩| 大陆偷拍与自拍| 波多野结衣一区麻豆| freevideosex欧美| 精品一区二区三区四区五区乱码 | 国产老妇伦熟女老妇高清| 婷婷色av中文字幕| 国产成人精品一,二区| 成人午夜精彩视频在线观看| 精品酒店卫生间| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 午夜福利影视在线免费观看| 少妇人妻精品综合一区二区| 电影成人av| 国产精品无大码| 一本色道久久久久久精品综合| 在线观看免费高清a一片| 在线 av 中文字幕| 女人高潮潮喷娇喘18禁视频| 久久99蜜桃精品久久| 久久精品国产自在天天线| 久久ye,这里只有精品| 九色亚洲精品在线播放| 午夜福利视频在线观看免费| 久久久久久久久久久免费av| 日韩一本色道免费dvd| 国产麻豆69| 亚洲精品国产色婷婷电影| 人妻人人澡人人爽人人| 丝袜在线中文字幕| 91成人精品电影| 少妇被粗大的猛进出69影院| 亚洲中文av在线| 亚洲精品aⅴ在线观看| 一区二区av电影网| 777米奇影视久久| 久久综合国产亚洲精品| 我要看黄色一级片免费的| 久久精品久久久久久久性| 国产精品久久久av美女十八| 在线精品无人区一区二区三| 高清黄色对白视频在线免费看| 久久精品国产自在天天线| av一本久久久久| 欧美 亚洲 国产 日韩一| 亚洲欧美中文字幕日韩二区| 熟女av电影| 中文欧美无线码| 啦啦啦视频在线资源免费观看| 久久精品国产a三级三级三级| 欧美激情高清一区二区三区 | 啦啦啦视频在线资源免费观看| 黑人巨大精品欧美一区二区蜜桃| 边亲边吃奶的免费视频| xxxhd国产人妻xxx| 一个人免费看片子| 一区二区三区精品91| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| 国产有黄有色有爽视频| 国产亚洲最大av| 亚洲成人一二三区av| videosex国产| 欧美国产精品一级二级三级| 精品国产超薄肉色丝袜足j| 久久鲁丝午夜福利片| 亚洲伊人久久精品综合| 久久久久国产精品人妻一区二区| 国产成人精品在线电影| 国产亚洲最大av| 如何舔出高潮| 熟女电影av网| 人人澡人人妻人| 99久久人妻综合| 亚洲欧美精品自产自拍| 观看美女的网站| 国产深夜福利视频在线观看| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 午夜精品国产一区二区电影| 高清黄色对白视频在线免费看| 青春草亚洲视频在线观看| 久久久久久人人人人人| 欧美精品一区二区大全| 男男h啪啪无遮挡| av在线老鸭窝| 成人黄色视频免费在线看| 欧美日韩国产mv在线观看视频| 男人添女人高潮全过程视频| av免费观看日本| 天天躁日日躁夜夜躁夜夜| 午夜福利乱码中文字幕| 精品一区二区三区四区五区乱码 | 国产亚洲午夜精品一区二区久久| 午夜精品国产一区二区电影| a级毛片在线看网站| 久久99热这里只频精品6学生| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片| 免费在线观看视频国产中文字幕亚洲 | 三上悠亚av全集在线观看| 一级毛片黄色毛片免费观看视频| 午夜免费鲁丝| 国产日韩欧美在线精品| 一区二区三区精品91| videossex国产| 成年女人在线观看亚洲视频| 欧美日韩一级在线毛片| 欧美精品亚洲一区二区| 亚洲人成77777在线视频| 亚洲精品国产一区二区精华液| 免费看av在线观看网站| 边亲边吃奶的免费视频| 9191精品国产免费久久| 亚洲国产精品999| 国产成人免费观看mmmm| 亚洲av综合色区一区| 国产av精品麻豆| 97在线视频观看| 女的被弄到高潮叫床怎么办| 91午夜精品亚洲一区二区三区| 亚洲 欧美一区二区三区| 老女人水多毛片| 只有这里有精品99| 国产片特级美女逼逼视频| 午夜福利,免费看| 99热网站在线观看| 飞空精品影院首页| 午夜日本视频在线| av电影中文网址| 少妇被粗大的猛进出69影院| 久久久久视频综合| 另类精品久久| 91午夜精品亚洲一区二区三区| av卡一久久| 久久久精品免费免费高清| 99热国产这里只有精品6| 国产精品嫩草影院av在线观看| 亚洲久久久国产精品| 日韩大片免费观看网站| 亚洲国产欧美日韩在线播放| 午夜日本视频在线| 国产精品熟女久久久久浪| 亚洲成av片中文字幕在线观看 | 自线自在国产av| 边亲边吃奶的免费视频| 一本久久精品| 黄色一级大片看看| 国产精品三级大全| 日本91视频免费播放| 91aial.com中文字幕在线观看| 亚洲欧美清纯卡通| videosex国产| 国产成人av激情在线播放| 久久久久久久亚洲中文字幕| 国产国语露脸激情在线看| 国产精品99久久99久久久不卡 | 一区在线观看完整版| 精品一区二区三区四区五区乱码 | 欧美日韩综合久久久久久| 亚洲精品乱久久久久久| 热re99久久国产66热| 亚洲伊人色综图| a 毛片基地| 在线精品无人区一区二区三| av免费在线看不卡| 欧美国产精品va在线观看不卡| 日韩电影二区| 一区二区三区精品91| 久久国产亚洲av麻豆专区| 久久久久久人妻| 母亲3免费完整高清在线观看 | www日本在线高清视频| 一二三四在线观看免费中文在| 国产日韩欧美亚洲二区| 免费观看在线日韩| 两性夫妻黄色片| 成人18禁高潮啪啪吃奶动态图| 欧美成人精品欧美一级黄| 亚洲视频免费观看视频| 久久99蜜桃精品久久| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说| 日本欧美视频一区| 汤姆久久久久久久影院中文字幕| 国产高清国产精品国产三级| 青春草国产在线视频| 岛国毛片在线播放| av网站在线播放免费| 亚洲成人一二三区av| 亚洲精品一二三| 91aial.com中文字幕在线观看| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 男女下面插进去视频免费观看| 成人影院久久| 妹子高潮喷水视频| 国产白丝娇喘喷水9色精品| 汤姆久久久久久久影院中文字幕| 久久韩国三级中文字幕| 韩国av在线不卡| 久久久久久久大尺度免费视频| 中文字幕人妻丝袜制服| 国产片内射在线| 婷婷色av中文字幕| 午夜福利网站1000一区二区三区| 美女xxoo啪啪120秒动态图| 伊人久久大香线蕉亚洲五| 中文字幕另类日韩欧美亚洲嫩草| 久久热在线av| 韩国精品一区二区三区| 免费观看在线日韩| 精品少妇久久久久久888优播| 午夜福利视频在线观看免费| 青草久久国产| 日本欧美国产在线视频| 欧美人与善性xxx| 精品一区二区三区四区五区乱码 | 看十八女毛片水多多多| 下体分泌物呈黄色| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 免费不卡的大黄色大毛片视频在线观看| 2021少妇久久久久久久久久久| 日韩欧美一区视频在线观看| 成人毛片a级毛片在线播放| 高清不卡的av网站| 精品一区二区三卡| 一区二区三区精品91| 午夜日韩欧美国产| 中文字幕精品免费在线观看视频| 日韩欧美一区视频在线观看| 制服诱惑二区| 2021少妇久久久久久久久久久| 岛国毛片在线播放| 国产精品 欧美亚洲| 午夜老司机福利剧场| 超碰97精品在线观看| 久久久欧美国产精品| 麻豆av在线久日| 亚洲,一卡二卡三卡| 久久久国产欧美日韩av| 国产av码专区亚洲av| 国产一区亚洲一区在线观看| 麻豆精品久久久久久蜜桃| 国产有黄有色有爽视频| 9色porny在线观看| 熟妇人妻不卡中文字幕| 一级片免费观看大全| 欧美日韩视频高清一区二区三区二| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费| 亚洲内射少妇av| 韩国av在线不卡| 久久久久国产精品人妻一区二区| 欧美日韩亚洲高清精品| 免费女性裸体啪啪无遮挡网站| 91久久精品国产一区二区三区| 久久久久久人妻| 最新的欧美精品一区二区| 色播在线永久视频| 成年动漫av网址| 精品国产一区二区久久| 日韩一区二区三区影片| 交换朋友夫妻互换小说| 男男h啪啪无遮挡| 尾随美女入室| 国产精品免费大片| 亚洲av综合色区一区| 日韩,欧美,国产一区二区三区| 久久久欧美国产精品| 久久精品久久精品一区二区三区| 欧美bdsm另类| tube8黄色片| 国产探花极品一区二区| 在线观看免费日韩欧美大片| 国产毛片在线视频| 国产女主播在线喷水免费视频网站| 亚洲欧美成人综合另类久久久| 亚洲av电影在线观看一区二区三区| 亚洲五月色婷婷综合| 午夜福利一区二区在线看| 在线天堂最新版资源| 中文精品一卡2卡3卡4更新| 777久久人妻少妇嫩草av网站| 国产日韩欧美在线精品| 看十八女毛片水多多多| 国产熟女午夜一区二区三区| 国产精品秋霞免费鲁丝片| 啦啦啦视频在线资源免费观看| 伊人久久大香线蕉亚洲五| 国产高清不卡午夜福利| 免费播放大片免费观看视频在线观看| 久久久久精品久久久久真实原创| 免费久久久久久久精品成人欧美视频| 日本欧美国产在线视频| 天天影视国产精品| 寂寞人妻少妇视频99o| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| av不卡在线播放| 久久久久久久国产电影| 亚洲国产精品999| 成人亚洲欧美一区二区av| 纵有疾风起免费观看全集完整版| 日韩人妻精品一区2区三区| 99香蕉大伊视频| av一本久久久久| 一二三四在线观看免费中文在| 国产成人精品一,二区| 夫妻性生交免费视频一级片| 少妇人妻久久综合中文| 亚洲国产最新在线播放| 熟妇人妻不卡中文字幕| 成人亚洲欧美一区二区av| 国产精品人妻久久久影院| 亚洲国产色片| 秋霞伦理黄片| 亚洲精品日韩在线中文字幕| 精品第一国产精品| 熟妇人妻不卡中文字幕| 成人毛片a级毛片在线播放| 成人国语在线视频| 麻豆精品久久久久久蜜桃| 视频在线观看一区二区三区| 一边摸一边做爽爽视频免费| 成人手机av| 天堂俺去俺来也www色官网| 亚洲av国产av综合av卡| 国产精品一二三区在线看| 两性夫妻黄色片| 国产精品久久久久久精品古装| 久久99蜜桃精品久久| 黄网站色视频无遮挡免费观看| 国产极品天堂在线| 韩国av在线不卡| 国产极品天堂在线| 久热久热在线精品观看| 最新中文字幕久久久久| 伊人久久国产一区二区| 18禁裸乳无遮挡动漫免费视频| 免费黄频网站在线观看国产| 777米奇影视久久| 久久国产精品男人的天堂亚洲| 一本色道久久久久久精品综合| 在线观看人妻少妇| 少妇的丰满在线观看| 亚洲国产毛片av蜜桃av| 欧美日韩亚洲国产一区二区在线观看 | 波野结衣二区三区在线| 亚洲国产最新在线播放| 国产成人精品久久二区二区91 | 麻豆精品久久久久久蜜桃| 两性夫妻黄色片| 女人被躁到高潮嗷嗷叫费观| 精品一区二区免费观看| 少妇精品久久久久久久| 熟女少妇亚洲综合色aaa.| 一区福利在线观看| 制服人妻中文乱码| 国产毛片在线视频| 91精品伊人久久大香线蕉| 波野结衣二区三区在线| 老司机影院成人| 人成视频在线观看免费观看| 天堂俺去俺来也www色官网| 婷婷色综合大香蕉| 乱人伦中国视频| 女人精品久久久久毛片| 欧美97在线视频| www.精华液| 熟女av电影| 久久久久久久精品精品| av在线播放精品| 在线精品无人区一区二区三| 大片免费播放器 马上看| 国产男女超爽视频在线观看| 欧美人与性动交α欧美软件| 久久久久久久久久久久大奶| 建设人人有责人人尽责人人享有的| 街头女战士在线观看网站| 美女视频免费永久观看网站| 一区二区av电影网| 欧美精品高潮呻吟av久久| 久热这里只有精品99| 多毛熟女@视频| 18禁国产床啪视频网站| 亚洲av电影在线观看一区二区三区| 日本黄色日本黄色录像| 午夜福利视频精品| 香蕉丝袜av| 亚洲av免费高清在线观看| 伊人亚洲综合成人网| 欧美成人午夜精品| 亚洲欧美成人综合另类久久久| 波多野结衣av一区二区av| 自拍欧美九色日韩亚洲蝌蚪91| 久久久欧美国产精品| 最近中文字幕高清免费大全6| 女人久久www免费人成看片| 久久精品国产综合久久久| 狠狠婷婷综合久久久久久88av| 曰老女人黄片| 亚洲精品久久久久久婷婷小说| 国产无遮挡羞羞视频在线观看| 99热国产这里只有精品6| 午夜福利视频在线观看免费| √禁漫天堂资源中文www| 18禁国产床啪视频网站| 亚洲成人手机| 高清黄色对白视频在线免费看| 日韩一区二区三区影片| 国产精品秋霞免费鲁丝片| 欧美在线黄色| 亚洲情色 制服丝袜| av有码第一页| 99久国产av精品国产电影| 日韩av在线免费看完整版不卡| 婷婷色av中文字幕| 性少妇av在线| 久久这里有精品视频免费| 国产精品女同一区二区软件| 亚洲av电影在线进入| www日本在线高清视频| 婷婷色麻豆天堂久久| 18禁观看日本| 亚洲精品国产一区二区精华液| 国产伦理片在线播放av一区| 毛片一级片免费看久久久久| 午夜福利网站1000一区二区三区| 999精品在线视频| 精品一区二区三区四区五区乱码 | 在线 av 中文字幕| 国产黄频视频在线观看| 91成人精品电影| 中文字幕最新亚洲高清| 女人久久www免费人成看片| 老司机影院成人| 欧美精品一区二区免费开放| 国产乱人偷精品视频| 日韩av免费高清视频| 老汉色av国产亚洲站长工具| 春色校园在线视频观看| 婷婷色麻豆天堂久久| 老熟女久久久| 国产视频首页在线观看| 女性生殖器流出的白浆| 男人爽女人下面视频在线观看| 大香蕉久久成人网| 亚洲精品一区蜜桃| 国产精品 欧美亚洲| 亚洲国产日韩一区二区| 亚洲一区中文字幕在线| 如日韩欧美国产精品一区二区三区| 伊人久久国产一区二区| 国产视频首页在线观看| 国产男女超爽视频在线观看| 在线免费观看不下载黄p国产| 精品国产一区二区久久| 日韩一区二区视频免费看| 国产免费视频播放在线视频| 满18在线观看网站| 国产精品久久久久久精品古装| 亚洲成国产人片在线观看| 国产精品国产av在线观看| 色播在线永久视频| 久久免费观看电影| 天天躁夜夜躁狠狠躁躁| 欧美国产精品一级二级三级| 久久女婷五月综合色啪小说| 精品少妇一区二区三区视频日本电影 | 日本免费在线观看一区| 麻豆精品久久久久久蜜桃| 一级毛片 在线播放| 久久久欧美国产精品| av国产久精品久网站免费入址| 婷婷色麻豆天堂久久| 97在线人人人人妻| 久久ye,这里只有精品| 午夜福利在线免费观看网站| 狂野欧美激情性bbbbbb| 视频区图区小说| 秋霞伦理黄片| 精品少妇内射三级| a级毛片在线看网站| 国产爽快片一区二区三区| 黑人猛操日本美女一级片| 国产成人精品福利久久| 曰老女人黄片| 久久这里只有精品19| 国产不卡av网站在线观看| 精品久久久精品久久久| 91aial.com中文字幕在线观看| 巨乳人妻的诱惑在线观看| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| 婷婷色综合www| 这个男人来自地球电影免费观看 | 叶爱在线成人免费视频播放| 国产精品国产三级国产专区5o| 搡女人真爽免费视频火全软件| 国产色婷婷99| 男女边吃奶边做爰视频| 天天影视国产精品| 亚洲av福利一区| 丝袜喷水一区| 哪个播放器可以免费观看大片| 午夜精品国产一区二区电影| 少妇被粗大猛烈的视频| 免费看av在线观看网站| 男的添女的下面高潮视频| 色婷婷久久久亚洲欧美| 看十八女毛片水多多多| 啦啦啦中文免费视频观看日本| 在线天堂中文资源库| 久久韩国三级中文字幕| 90打野战视频偷拍视频| 国产成人免费无遮挡视频| 国产成人精品在线电影| av卡一久久| 日本免费在线观看一区| 黄色视频在线播放观看不卡| 免费高清在线观看视频在线观看| 国产成人欧美| 色婷婷av一区二区三区视频| 在线观看一区二区三区激情| 国产亚洲午夜精品一区二区久久| 亚洲,欧美精品.| 一二三四中文在线观看免费高清| 精品国产一区二区三区久久久樱花| 五月天丁香电影| 欧美亚洲 丝袜 人妻 在线| 国产成人欧美| 亚洲一级一片aⅴ在线观看| 成年动漫av网址| 午夜免费鲁丝| 午夜福利在线免费观看网站| a级毛片黄视频| 成人手机av| 在线观看人妻少妇| 久久久久人妻精品一区果冻| 亚洲av电影在线进入|