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

    東北冷渦背景下中尺度低渦的演變成因及其對(duì)MCS的影響分析

    2021-10-21 02:29:12喬娜錢鵬周勍張晨昕吳昕悅
    暴雨災(zāi)害 2021年5期
    關(guān)鍵詞:冷渦緯向散度

    喬娜,錢鵬,周勍,張晨昕,吳昕悅

    (江蘇省鎮(zhèn)江市氣象局,鎮(zhèn)江 212000)

    引言

    暴雨是中國(guó)主要的災(zāi)害性天氣之一,國(guó)內(nèi)外學(xué)者早在20世紀(jì)80年代就已經(jīng)開始關(guān)注引發(fā)強(qiáng)降水的中小尺度系統(tǒng)。研究表明,中尺度對(duì)流系統(tǒng)(MCS)在同一個(gè)地區(qū)重復(fù)出現(xiàn)往往是產(chǎn)生持續(xù)性暴雨的直接原因(諶偉等,2017),MCS的發(fā)生發(fā)展不僅受到天氣尺度系統(tǒng)的影響,同時(shí)受到中尺度系統(tǒng)的制約(Orlanski,1975;貝耐芳等,2002),因此在研究強(qiáng)對(duì)流的同時(shí),加強(qiáng)對(duì)不同系統(tǒng)之間相互影響的分析十分必要。

    東北冷渦是指我國(guó)東北地區(qū)具有一定強(qiáng)度,能維持一段時(shí)間,且有深厚冷空氣的高空冷性的氣旋性渦旋。氣象研究中常用的定義為500 hPa 天氣圖上,115°—45°E、35°—60°N 范圍內(nèi),具有一定強(qiáng)度(至少有一根閉合等高線),能維持3 d 或以上,且有冷中心或者冷槽配合的高空氣旋性渦旋(鄭秀雅,1992;Zhao et al.,2007;Fu and Sun,2012)。東北冷渦的發(fā)生、發(fā)展、滯留、衰亡過(guò)程均對(duì)大氣環(huán)流起著至關(guān)重要的作用,在其演變過(guò)程中,它與高空急流、鋒面等系統(tǒng)相配合形成溫壓不對(duì)稱結(jié)構(gòu),進(jìn)而為強(qiáng)對(duì)流發(fā)生提供良好的環(huán)境場(chǎng),并且它經(jīng)常與中尺度系統(tǒng)相互作用,為強(qiáng)對(duì)流的發(fā)展提供有利的環(huán)流背景(白人海等,1998;孫力等,2000;張立祥等,2008,2009;王宇欣等,2014)。大尺度流場(chǎng)為強(qiáng)對(duì)流的發(fā)生提供了有利的環(huán)流背景,而中尺度系統(tǒng)同樣也是影響強(qiáng)對(duì)流發(fā)生發(fā)展的關(guān)鍵,研究表明,中尺度低渦的活動(dòng)會(huì)直接影響強(qiáng)對(duì)流的演變及發(fā)展,進(jìn)而引發(fā)暴雨的出現(xiàn)(東高紅等,2013;何光碧等,2014;高守亭等,2019)。程麟生等(2001)利用MM5 模式對(duì)“98.7”特大暴雨過(guò)程進(jìn)行模擬指出,低層環(huán)境場(chǎng)中存在一個(gè)明顯的中尺度低渦,低渦區(qū)域強(qiáng)烈的輻合作用是這次大暴雨出現(xiàn)的主要原因。不少學(xué)者指出對(duì)流層中低層正垂直渦度擾動(dòng)有利于中尺度低渦的生成,而低渦輻合作用又為對(duì)流增強(qiáng)提供了有利的環(huán)境場(chǎng),進(jìn)而引發(fā)暴雨出現(xiàn)(徐亞梅等,2002;王智等,2003;孫建華等,2004;王芬等,2015)。張慶紅等(2000)對(duì)一次MCS 過(guò)程進(jìn)行數(shù)值模擬指出,當(dāng)中尺度對(duì)流發(fā)生時(shí),潛熱釋放會(huì)引發(fā)中尺度低壓出現(xiàn),該低壓會(huì)直接影響低空急流的強(qiáng)度,進(jìn)而影響對(duì)流的發(fā)展。王海東等(2008)指出低渦東移會(huì)引發(fā)大暴雨的生成。

    以上研究成果表明,關(guān)于東北冷渦和中尺度低渦分別從環(huán)流場(chǎng)上影響強(qiáng)對(duì)流的分析較多,但是對(duì)于不同系統(tǒng)相互作用及其如何影響強(qiáng)對(duì)流發(fā)生發(fā)展的研究較少,因此本文利用WRF數(shù)值模式對(duì)東北冷渦背景下的一次中尺度低渦伴隨MCS 的過(guò)程進(jìn)行模擬分析與診斷,深入了解東北冷渦影響下中尺度低渦發(fā)生發(fā)展的原因及其對(duì)MCS的物理影響機(jī)制,以期提高低渦影響下強(qiáng)對(duì)流預(yù)報(bào)的準(zhǔn)確率。

    1 資料和方法

    1.1 選取資料

    使用的資料包括:美國(guó)國(guó)家大氣研究中心(NCAR)和美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP)提供的FNL (Final Operational Global Analysis)全球分析資料,時(shí)間間隔為每6 h 一次,空間分辨率為1°×1°,屬于GRIB2 碼全球資料;全國(guó)常規(guī)地面觀測(cè)和探空資料;SWAN系統(tǒng)處理后的多普勒雷達(dá)資料,包括濱州、滄州等11 個(gè)雷達(dá)站的觀測(cè)資料,時(shí)間間隔為每6 min一次,空間分辨率為0.01°×0.01°;中尺度數(shù)值模式WRFV3.6(ARW)輸出的高分辨率模擬資料。

    1.2 模式方案設(shè)計(jì)

    本文采用非靜力中尺度數(shù)值模式WRFV3.6 進(jìn)行數(shù)值模擬研究。模擬時(shí)間從2015年8月3日00時(shí)—3日18時(shí)(世界時(shí),下同)共18 h,初始場(chǎng)及側(cè)邊界場(chǎng)采用NCEP 1°×1°逐6 h的FNL全球分析資料。模擬采用三重雙向嵌套,中心經(jīng)緯度為(117.037°E、36.316°N),三層嵌套網(wǎng)格數(shù)分別是349×349,454×496,598×595,水平分辨率分別是18.0 km、6.0 km、2.0 km,垂直方向共分為35 層,第一、二層每1 h 輸出一次模擬結(jié)果,第三層每10 min 輸出一次模擬結(jié)果。采用的參數(shù)化方案包括:Thompson aerosol-aware 微物理方案(胡向軍等,2008),RRTM 長(zhǎng)波輻射方案(Mlawer et al.,1997),Dud?hia 短波輻射方案(Lacis,1974),K-F 積云對(duì)流方案(第三層不使用) (Kain et al.,1990,1993),Monin-Obukhov近地面層方案(章國(guó)材,2004),YSU 邊界層方案(Hong et al.,2006),Noah 陸面過(guò)程方案(趙向軍,2017)等,文中主要分析第三區(qū)域D03的模擬結(jié)果。

    2 天氣過(guò)程實(shí)況

    2.1 環(huán)流背景

    從天氣形勢(shì)來(lái)看,2015 年8 月3 日06 時(shí)(世界時(shí),下同) (圖1a),500 hPa 呈現(xiàn)兩脊一槽的形勢(shì),高壓脊分別位于貝加爾湖和鄂霍次克海西部,槽線位于內(nèi)蒙古高原-黃土高原-四川盆地一帶。東北冷渦位于黑龍江以西,中心在(120°E、53°N)附近,閉合環(huán)流尺度達(dá)1 500 km 以上,強(qiáng)度在556 dagpm 以下,是深厚的冷性低壓系統(tǒng)。低緯西北太平洋副熱帶高壓(以下簡(jiǎn)稱副高)脊線西伸至110°E 附近并穩(wěn)定維持,副高西北側(cè)的西南暖濕氣流為強(qiáng)對(duì)流天氣發(fā)生提供良好水汽條件。200 hPa 高空急流位于東北冷渦中心位置以南,中心最大風(fēng)速達(dá)50 m·s-1以上,強(qiáng)度較強(qiáng)。中尺度低渦和對(duì)流發(fā)生在500 hPa 東北冷渦槽前的偏西氣流中,200 hPa高空急流入口區(qū)右側(cè),高層對(duì)應(yīng)強(qiáng)輻散區(qū)。850 hPa(圖1c)由于對(duì)流區(qū)和高空急流右后方強(qiáng)烈的輻散加強(qiáng)低空急流左前方輻合,從而導(dǎo)致中尺度低渦(400 km)生成,隨后10 時(shí)(圖略)低渦輻合作用進(jìn)一步促進(jìn)對(duì)流單體組織成帶狀中尺度對(duì)流系統(tǒng)(MCS)。

    12 時(shí)(圖1b),500 hPa 東北冷渦緩慢東移,中心位置向東偏移,槽線移至東北平原—華北平原一帶,200 hPa 高空急流中心也略微東移,強(qiáng)度有小幅度減弱,中尺度低渦和MCS 仍位于東北冷渦槽前、高空急流的右后方。850 hPa 低空急流向東北方向移動(dòng)(圖1d),風(fēng)速增強(qiáng),低渦中心略微東北移,伴隨143 dagpm等高線范圍的擴(kuò)大而強(qiáng)度增加,高度延伸至800 hPa,MCS處于低空急流與中尺度低渦之間輻合區(qū)內(nèi)。

    圖1 2015年8月3日06時(shí)(a)、12時(shí)(b)的200 hPa全風(fēng)速場(chǎng)(陰影,單位:m·s-1)和500 hPa高度場(chǎng)(黑色實(shí)線,單位:dagpm)、溫度場(chǎng)(紅色實(shí)線,單位:℃)和風(fēng)場(chǎng)(風(fēng)向標(biāo),單位:m·s-1)與06時(shí)(c)、12時(shí)(d)的850 hPa全風(fēng)速場(chǎng)(陰影,單位:m·s-1)、高度場(chǎng)(藍(lán)色實(shí)線,單位:dagpm)和流場(chǎng)(黑色流線,單位:m·s-1)(紅色三角表示強(qiáng)對(duì)流區(qū),棕色實(shí)線表示500 hPa槽線)Fig.1 Wind speed field(shadow,unit:m·s-1)at 200 hPa and geopotential height field(black solid line,unit:dagpm),temperature field(red solid line,unit:℃)and wind field(wind vane,unit:m·s-1)at 500 hPa at(a)06∶00 UTC and(b)12∶00 UTC on 3August 2015,and wind speed field(shadow,unit:m·s-1),geopotential height field(blue solid line,unit:dagpm)and flow field(black line,unit:m·s-1)at 850 hPa at(c)06∶00 UTC and(d)12∶00 UTC on 3 August 2015(The red triangle represents severe convection zone and the solid brown line represents the trough-line at 500 hPa).

    2.2 中尺度低渦和MCS實(shí)況分析

    低渦是低壓在流場(chǎng)上的表現(xiàn),雷達(dá)回波分布可以體現(xiàn)MCS的演變特征。通過(guò)上節(jié)分析可知,東北冷渦為中尺度低渦和對(duì)流的發(fā)生提供了良好的天氣尺度環(huán)流背景,同時(shí)中尺度低渦為MCS的組織生成有重要貢獻(xiàn)。由圖2 可知,06 時(shí)(圖2a)東北冷渦與高低空急流等系統(tǒng)共同作用使得對(duì)流單體在山東西部、河南北部生成,同時(shí)由于對(duì)流區(qū)上方對(duì)應(yīng)假相當(dāng)位溫的密集帶,鋒生使得對(duì)流單體快速組織發(fā)展;高空急流右后方的輻散加強(qiáng)了低空急流左前方氣流的輻合,同時(shí)對(duì)流區(qū)輻合作用推進(jìn)了偏北風(fēng)與西南風(fēng)交匯,進(jìn)而導(dǎo)致低層中尺度低渦在此處生成。12時(shí)(圖2b)中尺度低渦范圍擴(kuò)大,強(qiáng)度增加,它與低空急流間的切變輻合加強(qiáng),風(fēng)場(chǎng)的輻合作用使得散亂的對(duì)流單體組織成帶狀中尺度對(duì)流系統(tǒng),即MCS生成發(fā)展。

    圖2 2015年8月3日06時(shí)(a)、12時(shí)(b)的850 hPa實(shí)況雷達(dá)回波反射率(陰影,單位:dBz)、流場(chǎng)(黑色流線,單位:m·s-1)和假相當(dāng)位溫分布(紅色實(shí)線,單位:K,黑色實(shí)心點(diǎn)代表滄州站)Fig.2 Observed radar reflectivity(shadow,unit:dBz),flow field(black line,unit:m·s-1)and distribution of pseudo-potential temperature(red solid line,unit:K)at 850 hPa at(a)06∶00 UTC and(b)12∶00 UTC on 3August 2015(The black solid point represents Cangzhou station).

    由于14 時(shí)左右為中尺度低渦和MCS 發(fā)展強(qiáng)盛的階段,因此圖3 給出了該時(shí)刻滄州站上空的雷達(dá)探測(cè)資料,由徑向風(fēng)速(圖3a)可以看出,山東西北部存在完整的中尺度低渦形態(tài),此時(shí)低渦正強(qiáng)盛發(fā)展。同時(shí),由滄州站上空雷達(dá)反射率(圖3b)可見,緊鄰低渦東南部的帶狀對(duì)流系統(tǒng)清晰,可見低渦與MCS的發(fā)展存在密切聯(lián)系。

    圖3 2015年8月3日14時(shí)滄州站仰角1.5°的徑向風(fēng)(a,單位:m·s-1,黃圈表示中尺度低渦)和仰角0.5°的雷達(dá)反射率(b,單位:dBz)Fig.3 (a)Radial wind(unit:m·s-1)at elevation angle 1.5°and(b)radar reflectance(unit:dBz)at elevation angle 0.5°from Cangzhou station at 14∶00 UTC 3 August 2015(The yellow circle represents the mesoscale vortex).

    總體而言,東北冷渦等天氣尺度系統(tǒng)為中尺度低渦和強(qiáng)對(duì)流的發(fā)生提供了有利的環(huán)流背景,同時(shí)中尺度低渦的發(fā)展又是促進(jìn)MCS生成的重要因素,因此此次過(guò)程中尺度低渦的演變成因和它對(duì)MCS 發(fā)展的影響將是本文討論的主要問(wèn)題。

    3 模擬結(jié)果檢驗(yàn)

    圖4 為本次模擬得到的環(huán)流場(chǎng),從中可見,06 時(shí)(圖4a、c)500 hPa呈現(xiàn)兩脊一槽形態(tài),東北冷渦位于黑龍江西部,低緯副高脊線明顯西伸,200 hPa高空急流位于東北冷渦南部,850 hPa存在明顯的低空急流,高低層系統(tǒng)共同作用促使對(duì)流和中尺度低渦生成,與實(shí)況類似。12 時(shí)(圖4b、d) 各層系統(tǒng)東移,強(qiáng)度增加,MCS 處于低空急流和中尺度低渦之間的切變輻合區(qū)內(nèi),與實(shí)況一致。

    圖4 2015年8月3日06時(shí)(a)、12時(shí)(b)模擬的200 hPa全風(fēng)速場(chǎng)(陰影,單位:m·s-1)和500 hPa高度場(chǎng)(黑色實(shí)線,單位:dagpm)、溫度場(chǎng)(紅色實(shí)線,單位:℃)和風(fēng)場(chǎng)(風(fēng)向標(biāo),單位:m·s-1)與06時(shí)(c)、12時(shí)(d)模擬的850 hPa全風(fēng)速場(chǎng)(陰影,單位:m·s-1)、高度場(chǎng)(藍(lán)色實(shí)線,單位:dagpm)和流場(chǎng)(黑色流線,單位:m·s-1)Fig.4 Simulated wind speed field(shadow,unit:m·s-1)at 200 hPa and geopotential height field(black solid line,unit:dagpm),temperature field(red solid line,unit:℃)and wind field(wind vane,unit:m·s-1)at 500 hPa at(a)06∶00 UTC and(b)12∶00 UTC on 3 August 2015,and simulated wind speed field(shadow,unit:m·s-1),geopotential height field(blue solid line,unit:dagpm)and flow field(black line,unit:m·s-1)at 850 hPa at(c)06∶00 UTC and(d)12∶00 UTC on 3August 2015.

    通過(guò)模擬的中尺度低渦(圖5)可以看出:一個(gè)中尺度低渦于06 時(shí)(圖5a)在113°E、38°N 附近生成,此處正是對(duì)流發(fā)生區(qū)域,隨后低渦中心東移,強(qiáng)度增加,范圍擴(kuò)大;12 時(shí)(圖5b)低渦范圍擴(kuò)大到400 km,散亂的對(duì)流單體已組織成東北-西南向的帶狀中尺度對(duì)流系統(tǒng)(MCS);14 時(shí)(圖5c)低渦和MCS 向東北方向移動(dòng),強(qiáng)度和范圍不斷增加;15時(shí)30分(圖5d)低渦繼續(xù)向東北方向移動(dòng),并逐漸減弱消失,MCS 也逐漸東移減弱。通過(guò)與實(shí)況(圖2、3)對(duì)比可知,本次模擬較好地再現(xiàn)了中尺度低渦發(fā)生的位置、時(shí)間和演變過(guò)程,以及它與強(qiáng)對(duì)流的對(duì)應(yīng)關(guān)系,因此可以利用本次模擬結(jié)果進(jìn)行進(jìn)一步的中尺度分析。

    圖5 2015年8月3日06時(shí)(a)、12時(shí)(b)、14時(shí)(c)、15時(shí)30分(d)的850 hPa模擬雷達(dá)回波反射率(陰影,單位:dBz)和流場(chǎng)(黑色流線,單位:m·s-1,紫色實(shí)心圓代表滄州站)Fig.5 Simulated radar reflectivity(shadow,unit:dBz)and flow field(black line,unit:m·s-1)at 850 hPa at(a)06∶00 UTC,(b)12∶00 UTC,(c)14∶00 UTC,and(d)15∶30 UTC on 3 August 2015(The purple solid point represents Cangzhou station).

    4 中尺度低渦和MCS的渦度與能量分析

    由上述分析可知,在東北冷渦的背景下,中尺度低渦的發(fā)生發(fā)展和MCS的演變關(guān)系密切,因此下面將從渦度方程出發(fā)探究中尺度低渦發(fā)生發(fā)展的原因,并利用能量轉(zhuǎn)換方程分析低渦是如何影響MCS發(fā)展的。

    4.1 中尺度低渦渦度收支分析

    為了研究中尺度低渦的發(fā)生發(fā)展機(jī)制,下面利用渦度方程(式(1))對(duì)低渦區(qū)域進(jìn)行渦度收支分析。p 坐標(biāo)系下的渦度方程(不考慮摩擦項(xiàng))(朱乾根等,2007)

    其中,A為相對(duì)渦度的局地變化項(xiàng),B為水平平流項(xiàng),C為垂直輸送項(xiàng),D為水平散度項(xiàng),E為傾側(cè)項(xiàng),u、v為水平風(fēng)速,ω為垂直速度,ζ為相對(duì)渦度,f為地轉(zhuǎn)渦度,β=

    4.1.1 局地相對(duì)渦度變化

    為了更好的觀察中尺度低渦發(fā)生發(fā)展過(guò)程中垂直方向上渦度的變化,選取低渦初生到消亡階段每半小時(shí)進(jìn)行時(shí)間平均,低渦區(qū)域5°×4°的范圍進(jìn)行區(qū)域平均,進(jìn)而對(duì)相對(duì)渦度進(jìn)行計(jì)算分析。從圖6可見,在中尺度低渦演變過(guò)程中,500 hPa以下始終以正渦度為主,最大正渦度集中在700—800 hPa,500 hPa 以上以負(fù)渦度為主;低渦最高延伸至800 hPa,06—15時(shí)低渦不斷加強(qiáng),800 hPa以下相對(duì)渦度不斷增大,850 hPa相對(duì)渦度最大增至8×10-5s-1以上,可以推測(cè)在此期間有中層正渦度持續(xù)的垂直下傳;15 時(shí)之后,低渦開始減弱,800 hPa以下相對(duì)渦度有所減小,850 hPa相對(duì)渦度減小明顯。

    圖6 區(qū)域(5°×4°)時(shí)間(每0.5 h)平均的中尺度低渦相對(duì)渦度的垂直廓線(單位:10-5s-1)(不同顏色的線條表示不同時(shí)次,橫虛線表示850 hPa高度,橫實(shí)線表示950 hPa高度)Fig.6 The vertical profile of relative vorticity of the mesoscale vortex with region(5°×4°)and time(0.5 h)mean(unit:10-5 s-1)(Different color lines indicate different times,and the horizontal dotted line represents 850 hPa,and the horizontal solid line represents 950 hPa).

    4.1.2 渦度方程各項(xiàng)貢獻(xiàn)

    為了探究在中尺度低渦發(fā)生發(fā)展過(guò)程中渦度變化的機(jī)制,根據(jù)渦度方程,分別選取了中尺度低渦發(fā)生(06 時(shí))、強(qiáng)盛(14 時(shí)30 分)、減弱(15 時(shí))階段的三個(gè)時(shí)次,計(jì)算垂直速度、散度和渦度方程中各項(xiàng)的變化來(lái)分析渦度變化的原因。

    在中尺度低渦形成時(shí)的8月3日06時(shí),由圖7a可知,低渦區(qū)域整層的垂直速度為正值,最大值在600 hPa附近,600 hPa以下為水平輻合,600 hPa以上為水平輻散,850 hPa的垂直渦度為3×10-5s-1;由圖7d可知,低渦區(qū)域中低層正的垂直渦度主要受到垂直輸送項(xiàng)和水平散度項(xiàng)的影響,其中700 hPa 以下以水平散度項(xiàng)的正貢獻(xiàn)為主,700—500 hPa 以垂直輸送項(xiàng)的正貢獻(xiàn)為主,中高層負(fù)的垂直渦度主要受到水平平流項(xiàng)和傾側(cè)項(xiàng)影響。結(jié)合圖7a、d可知,低層水平輻合引起垂直渦度的增加,垂直速度的峰值和垂直輸送項(xiàng)的峰值均集中在600 hPa 附近,這表明低層由于水平散度的輻合作用增加的正渦度,將通過(guò)垂直平流向上輸送,從而有利于增強(qiáng)氣旋性渦度的局地變化。

    在中尺度低渦強(qiáng)盛時(shí)的8月3日14時(shí)30分,由圖7b 可知,低渦區(qū)域整層的垂直速度為正值,峰值有所增加,且位于300 hPa 附近,300 hPa 以下為弱輻合,300 hPa 以上為強(qiáng)輻散,這可能與200 hPa 高空急流的東移增強(qiáng)有關(guān),850 hPa 的垂直渦度增大至8×10-5s-1;由圖7e可知,低渦區(qū)域中低層正渦度主要受到垂直輸送項(xiàng)和水平散度項(xiàng)的積極作用,其中800 hPa 以下以水平散度項(xiàng)的積極作用為主,且強(qiáng)度比06 時(shí)有所增加,800—400 hPa以垂直輸送項(xiàng)的積極作用為主,高層負(fù)的垂直渦度主要受到傾側(cè)項(xiàng)影響。結(jié)合圖7b、e 可知,低層水平散度項(xiàng)的正作用加強(qiáng),使得低層低渦區(qū)域有正的相對(duì)渦度局地變化,正渦度加大,由于800 hPa以上垂直速度和垂直輸送項(xiàng)均為正值且有所加強(qiáng),因此將低層正渦度不斷向中高層輸送,使得800—400 hPa相對(duì)渦度局地變化加強(qiáng),氣旋性渦度增加。

    在中尺度低渦減弱時(shí)的8月3日15時(shí),由圖7c可知,低渦區(qū)域的上升運(yùn)動(dòng)有所減弱,峰值仍位于300 hPa附近,300 hPa 以下有弱輻合,300 hPa 以上有弱輻散,這與200 hPa 高空急流的減弱相關(guān)聯(lián),850 hPa 的垂直渦度減小至6×10-5s-1;由圖7f 可知,800 hPa 以下水平散度項(xiàng)的正值減小,800—400 hPa 垂直輸送項(xiàng)的正值減小,850 hPa以下垂直輸送項(xiàng)變?yōu)樨?fù)值,中高層負(fù)垂直渦度仍主要受到傾側(cè)項(xiàng)作用。結(jié)合圖7c、f 可知,800 hPa以下水平散度項(xiàng)的正貢獻(xiàn)減小,垂直輸送項(xiàng)變?yōu)樨?fù)貢獻(xiàn),相對(duì)渦度的局地變化由正值變?yōu)樨?fù)值,使得低層氣旋性渦度開始減弱,并且中高層垂直運(yùn)動(dòng)和垂直輸送項(xiàng)有所減弱,中尺度低渦開始減弱。

    圖7 2015年8月3日06時(shí)(a)、14時(shí)30分(b)、15時(shí)(c)的中尺度低渦區(qū)域平均(5°×4°)和時(shí)間平均(每0.5 h)的垂直速度(W,單位:10-2m·s-1)、渦度(Vor,單位:10-5s-1)和散度(Div,單位:10-5s-1)的垂直廓線和06時(shí)(d)、14時(shí)30分(e)、15時(shí)(f)的渦度方程各項(xiàng)(單位:10-9s-2)的垂直廓線(A為相對(duì)渦度的局地變化項(xiàng),B為水平平流項(xiàng),C為垂直輸送項(xiàng),D為水平散度項(xiàng),E為傾側(cè)項(xiàng))Fig.7 Vertical profiles on the mesoscale vortex region(5°×4°)and time(0.5 h)mean of vertical velocity(W,unit:10-2 m·s-1),vertical vorticity(Vor,unit:10-5 s-1)and divergence(Div,unit:10-5 s-1)at(a)06∶00 UTC,(b)14∶30 UTC,and(c)15∶00 UTC on 3 August 2015,and vertical profiles of terms of the vorticity equation(unit:10-9 s-2)at(d)06∶00 UTC,(e)14∶30 UTC,and(f)15∶00 UTC on 3 August 2015(A represents local variation term of relative vorticity,and B represents horizontal advection term,and Crepresents vertical transport term,and D represents horizontal divergence term,and E represents heeling term).

    4.1.3 渦度方程各項(xiàng)變化分析

    為進(jìn)一步探究渦度方程各項(xiàng)變化,作各項(xiàng)的時(shí)間-高度剖面如圖8。06—15時(shí)為中尺度低渦發(fā)生發(fā)展的階段,水平散度項(xiàng)和垂直輸送項(xiàng)正貢獻(xiàn)最大,水平平流項(xiàng)和傾側(cè)項(xiàng)負(fù)貢獻(xiàn)最大,800 hPa以下水平散度項(xiàng)正貢獻(xiàn)更明顯,即輻合作用強(qiáng),水平平流項(xiàng)負(fù)貢獻(xiàn)更明顯,800—400 hPa 垂直輸送項(xiàng)正貢獻(xiàn)更明顯,即向上輸送強(qiáng),傾側(cè)項(xiàng)負(fù)貢獻(xiàn)更明顯。水平散度項(xiàng)在800 hPa以下對(duì)局地渦度變化有明顯的積極作用,強(qiáng)輻合引起正的渦度收支,使得低層氣旋性渦度增大,中層有較弱的消極作用,引起負(fù)的渦度收支,使得中層反氣旋性渦度略微增大,這種低層氣旋性渦度,中層反氣旋性渦度的配置有利于渦旋和上升運(yùn)動(dòng)的發(fā)展。垂直輸送項(xiàng)在800—400 hPa 始終表現(xiàn)為很強(qiáng)的積極作用,使得局地相對(duì)渦度穩(wěn)健增長(zhǎng),并將低層正渦度向上輸送,有利于渦旋的發(fā)展。水平平流項(xiàng)在800 hPa以下對(duì)局地相對(duì)渦度為明顯的負(fù)貢獻(xiàn),有正渦度向周圍流失,800 hPa以上存在正貢獻(xiàn),周圍有正渦度向渦旋中心輸送,這種渦度平流下負(fù)上正的垂直分布形態(tài),有利于上升運(yùn)動(dòng)的發(fā)展,從而進(jìn)一步加強(qiáng)低層輻合。傾側(cè)項(xiàng)在800—400 hPa表現(xiàn)為負(fù)貢獻(xiàn),多有垂直渦度向水平渦度的轉(zhuǎn)換,產(chǎn)生負(fù)渦度收支,不利于氣旋性渦旋向上伸展。

    圖8 中尺度低渦區(qū)域平均(5°×4°)和時(shí)間平均(0.5 h)的渦度方程各項(xiàng)中水平平流項(xiàng)(a)、垂直輸送項(xiàng)(b)、水平散度項(xiàng)(c)、傾側(cè)項(xiàng)(d)和相對(duì)渦度的局地變化項(xiàng)(e)的時(shí)間-高度剖面圖(單位:10-9s-2)Fig.8 Time-height cross-sections for(a)horizontal advection term,(b)vertical transport term,(c)horizontal divergence term,(d)heeling term and(e)local variation term of relative vorticity(unit:10-9 s-2)of the vorticity equation on the average mesoscale vortex.

    15 時(shí)之后中尺度低渦開始減弱,800 hPa 以下水平散度項(xiàng)逐漸減弱,由正值轉(zhuǎn)為負(fù)值,垂直輸送項(xiàng)負(fù)值區(qū)有所伸展,水平平流項(xiàng)仍為負(fù)值,使得相對(duì)渦度的局地變化由正值變?yōu)樨?fù)值,正的渦度收支有所減弱,負(fù)的渦度收支逐漸出現(xiàn),從而引起低層的氣旋性渦度環(huán)流減弱;800—400 hPa垂直輸送項(xiàng)有所減弱,使得正的渦度收支減弱,傾側(cè)項(xiàng)有所減弱,但仍表現(xiàn)為負(fù)值,引起局地相對(duì)渦度負(fù)變化。

    綜合上述分析,中尺度低渦形成和發(fā)展時(shí)期,低層的水平輻合和低層向高層正渦度的垂直輸送,使得正的垂直渦度增加,氣旋性環(huán)流加強(qiáng),低渦發(fā)展;中尺度低渦減弱時(shí)期,低層輻合減弱、垂直輸送變?yōu)樨?fù)值、垂直輸送項(xiàng)有所減小,使得正渦度制造減弱,負(fù)渦度制造增加,氣旋性環(huán)流減弱,低渦消亡。

    4.2 強(qiáng)對(duì)流帶能量收支分析

    上面利用渦度方程分析了中尺度低渦發(fā)生發(fā)展的原因,而伴隨著低渦的演變,MCS生成并發(fā)展,那么中尺度低渦是如何影響強(qiáng)對(duì)流的,下面將利用Jiang等(1995)的能量收支方法,詳細(xì)探討強(qiáng)對(duì)流帶與中尺度低渦間的能量轉(zhuǎn)換。

    將相關(guān)變量分為緯向平均場(chǎng)和擾動(dòng)場(chǎng),進(jìn)而將大氣動(dòng)能和有效位能分成緯向平均部分和擾動(dòng)部分來(lái)進(jìn)行研究,其中有效位能可用來(lái)衡量在絕熱過(guò)程中位能轉(zhuǎn)化為動(dòng)能的最大量。本節(jié)主要分析式(2)和(3)中涉及的三種轉(zhuǎn)化過(guò)程和一種生成過(guò)程,其中忽略S(KE)的作用,相應(yīng)轉(zhuǎn)化過(guò)程和生成過(guò)程的計(jì)算方法如式(4)—(8)所示。

    式(2)和(3)分別表示擾動(dòng)動(dòng)能和擾動(dòng)有效位能隨時(shí)間的變化,其中AE表示擾動(dòng)有效位能,AZ表示緯向平均有效位能,KE表示擾動(dòng)動(dòng)能,KZ表示緯向平均動(dòng)能,G(AE)表示擾動(dòng)有效位能在非絕熱加熱過(guò)程中的生成率,S(KE)表示非保守外力對(duì)能量的變化率。C(X→Y)表示能量從X 形式轉(zhuǎn)化到Y(jié) 形式。在以上方程中,方括號(hào)表示物理量的緯向平均,物理量右上角的撇號(hào)表示相對(duì)于該平均值的偏差,其中u'、v'、w'分別表示x、y、z方向上的擾動(dòng)速度,θ'表示擾動(dòng)位溫,H'表示絕熱加熱率的擾動(dòng)。

    圖9是根據(jù)上式計(jì)算出的模式第三層D03(111°—123°E、31°—41°N)區(qū)域時(shí)間平均的各種能量轉(zhuǎn)換隨時(shí)間的變化,具體算法如下:計(jì)算各變量緯向平均值、擾動(dòng)值和擾動(dòng)乘積的緯向平均值等;依據(jù)式(4)—(8)計(jì)算各轉(zhuǎn)換項(xiàng);計(jì)算各轉(zhuǎn)換項(xiàng)的垂直積分;計(jì)算各轉(zhuǎn)換項(xiàng)垂直積分后的區(qū)域時(shí)間平均;進(jìn)而分析擾動(dòng)動(dòng)能和擾動(dòng)有效位能的變化。其中由于中尺度低渦范圍很大,基本覆蓋D03區(qū)域,因此緯向平均能量可以代表中尺度低渦的能量,而擾動(dòng)能量可以表示強(qiáng)對(duì)流的能量。由圖可知,06—12時(shí),擾動(dòng)有效位能的生成項(xiàng)G(AE)和緯向平均有效位能向擾動(dòng)有效位能的轉(zhuǎn)化項(xiàng)C(AZ→AE)均為正值,且G(AE)在10時(shí)達(dá)到最大,最大值為20×10-4W,C(AZ→AE)在12時(shí)達(dá)到最大,最大值約為34×10-4W,擾動(dòng)有效位能向擾動(dòng)動(dòng)能的轉(zhuǎn)化項(xiàng)C(AE→KE)較小,說(shuō)明在此期間強(qiáng)對(duì)流獲得并積累了較高的擾動(dòng)有效位能,而非絕熱加熱項(xiàng)和緯向平均有效位能向擾動(dòng)有效位能的轉(zhuǎn)化項(xiàng)均是擾動(dòng)有效位能增加的主要來(lái)源。同時(shí)緯向平均動(dòng)能向擾動(dòng)動(dòng)能的轉(zhuǎn)化項(xiàng)C(KZ→KE)也為正值,且不斷增長(zhǎng),說(shuō)明此段時(shí)間中尺度低渦不斷向強(qiáng)對(duì)流區(qū)域輸送動(dòng)能,以助其發(fā)生發(fā)展。12—15時(shí),擾動(dòng)有效位能的生成項(xiàng)G(AE)和緯向平均有效位能向擾動(dòng)有效位能的轉(zhuǎn)化項(xiàng)C(AZ→AE)呈下降趨勢(shì),緯向平均動(dòng)能向擾動(dòng)動(dòng)能的轉(zhuǎn)化項(xiàng)C(KZ→KE)也開始減弱,但仍表現(xiàn)為正值,同時(shí),擾動(dòng)有效位能向擾動(dòng)動(dòng)能的轉(zhuǎn)化項(xiàng)C(AE→KE)呈現(xiàn)出明顯的上升趨勢(shì),說(shuō)明在此期間,擾動(dòng)有效位能轉(zhuǎn)化為擾動(dòng)動(dòng)能明顯,擾動(dòng)動(dòng)能的增加促進(jìn)MCS加強(qiáng)。15時(shí)之后中尺度低渦和強(qiáng)對(duì)流帶開始減弱,不做具體分析。整個(gè)過(guò)程中尺度低渦和MCS之間的主要能量轉(zhuǎn)化如圖10所示,這樣的能量轉(zhuǎn)化有利于中尺度對(duì)流系統(tǒng)的發(fā)展增強(qiáng)。

    圖9 D03區(qū)域時(shí)間平均(每1 h)的垂直積分的不同形式能量轉(zhuǎn)化隨時(shí)間的變化(單位:10-4W)(黑色實(shí)線表示C(AE→KE),紅色長(zhǎng)虛線表示C(AZ→AE),深紫色短虛線表示C(KZ→KE),青色點(diǎn)劃線表示G(AE))Fig.9 The energy transformation chart of vertical integral different forms of D03 region and time(1 h)mean over time(unit:10-4W),and the solid black line represents C(AE→KE),and the long dashed red line represents C(AZ→AE),and the short dashed dark purple line represents C(KZ→KE),and the dotted blue line represents G(AE)

    圖10 能量轉(zhuǎn)化示意圖Fig.10 Diagram of energy transformation.

    依據(jù)上述分析,進(jìn)一步分析風(fēng)場(chǎng)對(duì)于強(qiáng)對(duì)流帶的影響。圖11 給出了D03 區(qū)域時(shí)間平均的垂直積分的緯向平均風(fēng)場(chǎng)和擾動(dòng)風(fēng)場(chǎng)隨時(shí)間的變化,緯向平均風(fēng)場(chǎng)代表中尺度低渦,擾動(dòng)風(fēng)場(chǎng)代表強(qiáng)對(duì)流。分析可知,11時(shí)之前,緯向平均風(fēng)場(chǎng)主要以西風(fēng)為主,伴有較弱的南風(fēng),垂直上升運(yùn)動(dòng)較強(qiáng),且不斷加大,擾動(dòng)風(fēng)場(chǎng)主要以西風(fēng)為主,伴有較弱的北風(fēng)分量,存在一定的下沉運(yùn)動(dòng),由此可知,中尺度低渦輸送的南方暖空氣與對(duì)流區(qū)域的北方冷空氣交匯,產(chǎn)生垂直運(yùn)動(dòng),進(jìn)而促進(jìn)能量轉(zhuǎn)化,強(qiáng)對(duì)流發(fā)生。11—15 時(shí),緯向平均風(fēng)場(chǎng)西風(fēng)減弱,南風(fēng)加強(qiáng),12時(shí)之后主要以西南風(fēng)為主,垂直運(yùn)動(dòng)減弱,但仍為正值,擾動(dòng)風(fēng)場(chǎng)西風(fēng)先增強(qiáng)后減弱,北風(fēng)擾動(dòng)逐漸變?yōu)槟巷L(fēng)擾動(dòng),14 時(shí)之后主要以西南風(fēng)為主,下沉運(yùn)動(dòng)減弱,出現(xiàn)上升運(yùn)動(dòng),由此可知,由于中尺度低渦帶來(lái)的強(qiáng)烈西南風(fēng)作用,使得強(qiáng)對(duì)流帶處的北風(fēng)逐漸消失,隨后主要受到西南氣流的控制,水汽條件較好,且擾動(dòng)下沉轉(zhuǎn)為上升運(yùn)動(dòng),進(jìn)而促進(jìn)強(qiáng)對(duì)流發(fā)展。

    圖11 D03區(qū)域時(shí)間平均(每1 h)的垂直積分的緯向平均風(fēng)場(chǎng)(a,單位:m·s-1)和擾動(dòng)風(fēng)場(chǎng)(b,單位:10-2m·s-1)隨時(shí)間的變化Fig.11 Vertical integral charts of(a)zonal mean wind field(unit:m·s-1)and(b)disturbed wind field(unit:10-2 m·s-1)of D03 region and time(1 h)mean over time.

    5 結(jié)論

    由于中尺度低渦和MCS常引發(fā)暴雨等災(zāi)害,本文對(duì)2015年8月3日東北冷渦背景下的一次中尺度低渦伴隨MCS過(guò)程進(jìn)行數(shù)值模擬與分析,研究了中尺度低渦演變成因和它對(duì)強(qiáng)對(duì)流發(fā)生發(fā)展的影響,得到以下結(jié)論:

    (1)此次過(guò)程發(fā)生于200 hPa高空急流右后方輻散區(qū),500 hPa 東北冷渦南側(cè),槽前偏西氣流中,850 hPa低空急流西北側(cè),高低層系統(tǒng)相配合為中尺度低渦和強(qiáng)對(duì)流的生成提供了有利的環(huán)流背景,同時(shí),中尺度低渦的輻合作用促使散亂的對(duì)流單體組織成帶狀中尺度對(duì)流系統(tǒng),即低渦的生成有利于MCS的組織形成。

    (2)利用渦度方程分析中尺度低渦發(fā)生發(fā)展的原因,低層正渦度的變化主要受到水平散度項(xiàng)和垂直輸送項(xiàng)的影響,其中水平散度項(xiàng)作用更明顯,而中層正渦度的變化主要受到垂直輸送項(xiàng)的影響,高層負(fù)渦度的變化主要受到水平平流項(xiàng)和傾側(cè)項(xiàng)的影響,其中傾側(cè)項(xiàng)作用更顯著。中尺度低渦形成和發(fā)展時(shí)期,低層的水平輻合和低層向上正渦度的垂直輸送,使得正渦度增加,低渦發(fā)展;中尺度低渦減弱時(shí)期,低層輻合減弱、垂直輸送變?yōu)樨?fù)值,使得正渦度制造減弱,負(fù)渦度制造增加,低渦消亡。

    (3)利用能量收支方程分析中尺度低渦和MCS之間的能量轉(zhuǎn)換,緯向平均場(chǎng)代表中尺度低渦,擾動(dòng)場(chǎng)代表強(qiáng)對(duì)流,可知在MCS 發(fā)展過(guò)程中,潛熱釋放作用和緯向平均有效位能向擾動(dòng)有效位能的轉(zhuǎn)化是擾動(dòng)有效位能增加的主要來(lái)源,而風(fēng)場(chǎng)變化所引起的垂直上升下沉運(yùn)動(dòng)是擾動(dòng)有效位能轉(zhuǎn)化為擾動(dòng)動(dòng)能的主要機(jī)制,擾動(dòng)動(dòng)能增加,強(qiáng)對(duì)流發(fā)展。中尺度低渦不僅為強(qiáng)對(duì)流的發(fā)展提供能量,低渦附近風(fēng)場(chǎng)的變化也是影響強(qiáng)對(duì)流演變的重要原因。

    猜你喜歡
    冷渦緯向散度
    4—9 月東北大陸型冷渦消亡方式及其統(tǒng)計(jì)特征*
    帶勢(shì)加權(quán)散度形式的Grushin型退化橢圓算子的Dirichlet特征值的上下界
    紗線強(qiáng)力對(duì)純棉平紋面料強(qiáng)力的影響
    利用掩星溫度數(shù)據(jù)推算大氣月平均緯向風(fēng)場(chǎng)
    2000-2017年?yáng)|北冷渦活動(dòng)氣候特征及分析
    黑龍江氣象(2020年1期)2020-05-20 05:57:06
    具有部分BMO系數(shù)的非散度型拋物方程的Lorentz估計(jì)
    一次冷渦減弱階段產(chǎn)生的颮線過(guò)程分析
    溫度對(duì)絲綢面料粘襯熱縮率的影響
    絲綢(2018年10期)2018-10-15 09:54:16
    H型群上一類散度形算子的特征值估計(jì)
    柞蠶絲面料在粘襯過(guò)程中的熱縮率變化分析
    中文字幕免费在线视频6| 自线自在国产av| 秋霞伦理黄片| 国产成人精品久久久久久| 最新中文字幕久久久久| 黄色毛片三级朝国网站| 亚洲美女搞黄在线观看| .国产精品久久| 亚洲色图 男人天堂 中文字幕 | 国产高清国产精品国产三级| 国产成人aa在线观看| 日本欧美视频一区| 成人毛片60女人毛片免费| 成人国语在线视频| 好男人视频免费观看在线| 香蕉精品网在线| 欧美激情 高清一区二区三区| 交换朋友夫妻互换小说| 又粗又硬又长又爽又黄的视频| 啦啦啦视频在线资源免费观看| 免费看av在线观看网站| 日韩成人伦理影院| 免费高清在线观看视频在线观看| 午夜免费鲁丝| a 毛片基地| 亚洲av成人精品一区久久| 日韩熟女老妇一区二区性免费视频| 日本午夜av视频| 黑人欧美特级aaaaaa片| 国产视频内射| 成年人免费黄色播放视频| 成人手机av| a级毛片免费高清观看在线播放| 免费日韩欧美在线观看| 在线观看一区二区三区激情| 夜夜看夜夜爽夜夜摸| 精品一区二区免费观看| 国产精品熟女久久久久浪| 老女人水多毛片| 亚洲av二区三区四区| 国产极品粉嫩免费观看在线 | 亚洲精华国产精华液的使用体验| 美女大奶头黄色视频| 777米奇影视久久| 大香蕉97超碰在线| 午夜免费观看性视频| 男的添女的下面高潮视频| 大片电影免费在线观看免费| 最近的中文字幕免费完整| 日韩一区二区视频免费看| 全区人妻精品视频| 亚洲av福利一区| 99久国产av精品国产电影| 欧美精品高潮呻吟av久久| 高清午夜精品一区二区三区| 一级a做视频免费观看| 国产老妇伦熟女老妇高清| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 伦理电影免费视频| 日本黄色片子视频| 女人久久www免费人成看片| 街头女战士在线观看网站| 在线观看国产h片| 亚洲精品第二区| 少妇的逼水好多| 午夜日本视频在线| 人人妻人人添人人爽欧美一区卜| 免费黄网站久久成人精品| 五月开心婷婷网| 看非洲黑人一级黄片| 91在线精品国自产拍蜜月| 欧美亚洲日本最大视频资源| 大陆偷拍与自拍| 久久久久久久久大av| 亚洲欧美清纯卡通| 亚洲人成网站在线观看播放| 国产精品国产三级专区第一集| 中文精品一卡2卡3卡4更新| 极品人妻少妇av视频| 黄色一级大片看看| 欧美日韩综合久久久久久| 国产有黄有色有爽视频| 亚洲第一av免费看| 色网站视频免费| 精品一区二区免费观看| 精品少妇内射三级| 国产精品久久久久久久久免| a级片在线免费高清观看视频| 99视频精品全部免费 在线| 搡女人真爽免费视频火全软件| 国产精品秋霞免费鲁丝片| 亚洲国产av影院在线观看| 亚洲国产av影院在线观看| 一本—道久久a久久精品蜜桃钙片| 国产精品女同一区二区软件| 一本一本综合久久| 99热这里只有精品一区| 亚洲丝袜综合中文字幕| 免费av不卡在线播放| 久久99蜜桃精品久久| 在线观看三级黄色| 精品少妇内射三级| 亚洲av综合色区一区| 男人添女人高潮全过程视频| .国产精品久久| 欧美亚洲 丝袜 人妻 在线| 黑人高潮一二区| 一二三四中文在线观看免费高清| 亚洲精品成人av观看孕妇| 97在线人人人人妻| 国产成人精品一,二区| 日本黄色日本黄色录像| 国产精品 国内视频| av免费在线看不卡| 日日爽夜夜爽网站| 建设人人有责人人尽责人人享有的| 国产高清三级在线| 女人久久www免费人成看片| 一区二区三区精品91| 国产国拍精品亚洲av在线观看| 久热这里只有精品99| a级片在线免费高清观看视频| 中文字幕人妻丝袜制服| 久久精品久久久久久噜噜老黄| 制服人妻中文乱码| 午夜激情av网站| 国产成人精品一,二区| 欧美日本中文国产一区发布| 久久久午夜欧美精品| 亚洲综合色网址| 最后的刺客免费高清国语| 如何舔出高潮| 久久精品人人爽人人爽视色| 国产极品天堂在线| 亚州av有码| 亚洲内射少妇av| 亚洲精品av麻豆狂野| 亚洲精品成人av观看孕妇| 黄色欧美视频在线观看| 成人亚洲欧美一区二区av| 国产 一区精品| 成人18禁高潮啪啪吃奶动态图 | 日本欧美国产在线视频| 亚洲欧美清纯卡通| 国产一级毛片在线| 在线观看www视频免费| 免费黄频网站在线观看国产| 热99久久久久精品小说推荐| 国产在线免费精品| 欧美精品一区二区大全| 成人黄色视频免费在线看| 中国三级夫妇交换| 亚洲av电影在线观看一区二区三区| 久久久久网色| 日日啪夜夜爽| 大香蕉久久网| 在线观看免费日韩欧美大片 | 亚洲丝袜综合中文字幕| 99视频精品全部免费 在线| 国产极品粉嫩免费观看在线 | 久久久久久人妻| www.av在线官网国产| 国产精品一国产av| 在线观看三级黄色| 欧美精品亚洲一区二区| 成年人免费黄色播放视频| 一本—道久久a久久精品蜜桃钙片| 九色亚洲精品在线播放| 人体艺术视频欧美日本| 天堂8中文在线网| 精品人妻一区二区三区麻豆| 久久久午夜欧美精品| 午夜福利视频精品| 亚洲精华国产精华液的使用体验| 亚洲综合色惰| 美女视频免费永久观看网站| 亚洲欧美精品自产自拍| 久久99热这里只频精品6学生| 日本91视频免费播放| 精品人妻偷拍中文字幕| 最后的刺客免费高清国语| 欧美变态另类bdsm刘玥| 最新中文字幕久久久久| 麻豆精品久久久久久蜜桃| .国产精品久久| 97在线视频观看| 特大巨黑吊av在线直播| 国产免费福利视频在线观看| 国产成人a∨麻豆精品| 日本91视频免费播放| freevideosex欧美| 高清欧美精品videossex| 亚洲欧美成人精品一区二区| 色吧在线观看| 久热久热在线精品观看| 男的添女的下面高潮视频| 国精品久久久久久国模美| 精品久久久噜噜| 亚洲美女黄色视频免费看| 国语对白做爰xxxⅹ性视频网站| 欧美变态另类bdsm刘玥| 黑人高潮一二区| 欧美丝袜亚洲另类| 亚洲欧美精品自产自拍| av卡一久久| 考比视频在线观看| 国产一区亚洲一区在线观看| 内地一区二区视频在线| 久久久亚洲精品成人影院| 亚洲精品av麻豆狂野| 久久精品熟女亚洲av麻豆精品| 最黄视频免费看| 五月天丁香电影| 九草在线视频观看| 欧美性感艳星| av专区在线播放| 国产极品粉嫩免费观看在线 | 免费观看在线日韩| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 亚洲av.av天堂| 一级黄片播放器| 精品国产露脸久久av麻豆| 丰满迷人的少妇在线观看| 精品亚洲成国产av| 美女主播在线视频| 亚洲五月色婷婷综合| 波野结衣二区三区在线| 午夜免费鲁丝| 国产精品人妻久久久久久| 精品一区二区免费观看| 中文字幕制服av| 母亲3免费完整高清在线观看 | 色吧在线观看| 内地一区二区视频在线| 亚洲精品日韩av片在线观看| 韩国av在线不卡| 一级毛片黄色毛片免费观看视频| 久久国产亚洲av麻豆专区| 精品一品国产午夜福利视频| 国产欧美另类精品又又久久亚洲欧美| 中文字幕免费在线视频6| 亚洲精品乱码久久久久久按摩| 大片免费播放器 马上看| 亚洲欧美成人精品一区二区| 日本猛色少妇xxxxx猛交久久| 欧美成人午夜免费资源| 国产精品无大码| 女性被躁到高潮视频| 少妇高潮的动态图| 免费大片18禁| 日韩精品免费视频一区二区三区 | 午夜久久久在线观看| 九九爱精品视频在线观看| 亚洲精品乱码久久久久久按摩| 午夜老司机福利剧场| 人妻一区二区av| 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 日韩精品有码人妻一区| 久热久热在线精品观看| 亚洲美女视频黄频| 中文字幕亚洲精品专区| av播播在线观看一区| 精品一区二区三区视频在线| 免费观看av网站的网址| 黄色怎么调成土黄色| 久久久精品免费免费高清| 一级毛片 在线播放| 欧美 日韩 精品 国产| 男女啪啪激烈高潮av片| 中文乱码字字幕精品一区二区三区| 性高湖久久久久久久久免费观看| av在线播放精品| 国产又色又爽无遮挡免| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 成人黄色视频免费在线看| 亚洲人与动物交配视频| 精品人妻偷拍中文字幕| 免费观看的影片在线观看| 精品酒店卫生间| 中文字幕最新亚洲高清| 考比视频在线观看| 97在线视频观看| 国产一区有黄有色的免费视频| 26uuu在线亚洲综合色| 韩国av在线不卡| 天天操日日干夜夜撸| 男女高潮啪啪啪动态图| 女性被躁到高潮视频| 一级黄片播放器| 人妻少妇偷人精品九色| 午夜免费观看性视频| 欧美日韩成人在线一区二区| 国产高清有码在线观看视频| av在线app专区| 中文天堂在线官网| 制服人妻中文乱码| 精品99又大又爽又粗少妇毛片| 一级毛片黄色毛片免费观看视频| 97在线视频观看| 在线亚洲精品国产二区图片欧美 | 卡戴珊不雅视频在线播放| videos熟女内射| 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 国产亚洲av片在线观看秒播厂| 亚洲av欧美aⅴ国产| 好男人视频免费观看在线| 国产一区二区三区综合在线观看 | 啦啦啦啦在线视频资源| 中文欧美无线码| 久久久久国产精品人妻一区二区| 91午夜精品亚洲一区二区三区| 日韩亚洲欧美综合| 26uuu在线亚洲综合色| 国产无遮挡羞羞视频在线观看| 美女大奶头黄色视频| 在现免费观看毛片| 我的老师免费观看完整版| 亚洲精品国产av成人精品| 国产成人精品在线电影| 99久久综合免费| 国产在线免费精品| 91精品伊人久久大香线蕉| 色5月婷婷丁香| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 国产乱来视频区| 狠狠精品人妻久久久久久综合| 国产永久视频网站| 午夜老司机福利剧场| 在线观看一区二区三区激情| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| 亚洲精品一二三| av又黄又爽大尺度在线免费看| 国产69精品久久久久777片| 国产高清不卡午夜福利| 中文字幕最新亚洲高清| 亚洲国产成人一精品久久久| 亚洲av欧美aⅴ国产| 亚洲欧美日韩卡通动漫| 丁香六月天网| 69精品国产乱码久久久| av卡一久久| 午夜视频国产福利| 久久99精品国语久久久| 啦啦啦在线观看免费高清www| 中文字幕亚洲精品专区| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 午夜免费观看性视频| 女人精品久久久久毛片| 最近的中文字幕免费完整| 免费av不卡在线播放| 亚洲欧洲精品一区二区精品久久久 | 夫妻性生交免费视频一级片| √禁漫天堂资源中文www| 极品少妇高潮喷水抽搐| 欧美日韩在线观看h| 日韩免费高清中文字幕av| 天天影视国产精品| 只有这里有精品99| 亚洲成人手机| 国产午夜精品一二区理论片| 亚洲av男天堂| 丰满饥渴人妻一区二区三| 欧美日韩一区二区视频在线观看视频在线| 一级片'在线观看视频| 晚上一个人看的免费电影| 亚洲精品自拍成人| 国精品久久久久久国模美| 亚洲精华国产精华液的使用体验| 99热网站在线观看| 美女视频免费永久观看网站| 亚洲国产精品一区二区三区在线| 9色porny在线观看| 久久精品国产鲁丝片午夜精品| 制服丝袜香蕉在线| 久久ye,这里只有精品| 亚洲不卡免费看| 男女边摸边吃奶| 啦啦啦中文免费视频观看日本| 亚洲少妇的诱惑av| 久久久亚洲精品成人影院| 国产免费视频播放在线视频| 日本欧美视频一区| 日本av免费视频播放| 日韩欧美精品免费久久| 精品一品国产午夜福利视频| 看非洲黑人一级黄片| 在线 av 中文字幕| 一个人看视频在线观看www免费| 日本vs欧美在线观看视频| 熟女人妻精品中文字幕| 十八禁网站网址无遮挡| 日韩中文字幕视频在线看片| 日韩一区二区三区影片| 亚洲成色77777| 欧美另类一区| 18禁动态无遮挡网站| 中文字幕人妻丝袜制服| 免费大片18禁| 国产一区有黄有色的免费视频| 中文欧美无线码| 精品久久久精品久久久| 亚洲国产精品999| 日韩精品有码人妻一区| 久久久久国产精品人妻一区二区| 少妇高潮的动态图| 精品久久久久久久久av| av线在线观看网站| 欧美激情国产日韩精品一区| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| 超色免费av| 成人无遮挡网站| 国产成人午夜福利电影在线观看| 各种免费的搞黄视频| 中文字幕制服av| 欧美日韩成人在线一区二区| a级毛片黄视频| .国产精品久久| 久久午夜综合久久蜜桃| av在线老鸭窝| 国产午夜精品久久久久久一区二区三区| 久久综合国产亚洲精品| 涩涩av久久男人的天堂| 国产精品偷伦视频观看了| 高清欧美精品videossex| 看免费成人av毛片| 国产精品欧美亚洲77777| 国产乱人偷精品视频| 中文欧美无线码| 热99国产精品久久久久久7| 99国产精品免费福利视频| 2018国产大陆天天弄谢| 亚洲国产成人一精品久久久| 九九在线视频观看精品| 亚洲伊人久久精品综合| 99热全是精品| 久久人妻熟女aⅴ| 久久韩国三级中文字幕| 91精品一卡2卡3卡4卡| 五月玫瑰六月丁香| 晚上一个人看的免费电影| 日韩人妻高清精品专区| 中文字幕精品免费在线观看视频 | 能在线免费看毛片的网站| 亚洲欧美成人精品一区二区| 日日撸夜夜添| 高清在线视频一区二区三区| 国产精品秋霞免费鲁丝片| 久久精品久久久久久久性| 日本黄大片高清| 2021少妇久久久久久久久久久| 在线播放无遮挡| 黑人猛操日本美女一级片| 久久这里有精品视频免费| 久久婷婷青草| 我的老师免费观看完整版| 人人澡人人妻人| 在线观看一区二区三区激情| 精品久久久久久久久亚洲| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 亚洲国产av新网站| 狂野欧美激情性bbbbbb| 久久久久久久久久久免费av| 夫妻午夜视频| 欧美97在线视频| 黑人高潮一二区| 精品人妻一区二区三区麻豆| 国产一区二区三区av在线| 免费日韩欧美在线观看| 夫妻性生交免费视频一级片| 卡戴珊不雅视频在线播放| 欧美激情国产日韩精品一区| 在线观看国产h片| h视频一区二区三区| 中国三级夫妇交换| 免费黄网站久久成人精品| 制服诱惑二区| 嘟嘟电影网在线观看| 欧美精品一区二区大全| 插阴视频在线观看视频| 久久久久国产网址| 搡老乐熟女国产| 少妇精品久久久久久久| 日韩制服骚丝袜av| 97精品久久久久久久久久精品| 九九在线视频观看精品| 久久韩国三级中文字幕| 在线精品无人区一区二区三| 男女国产视频网站| 久久这里有精品视频免费| 精品国产一区二区久久| 欧美丝袜亚洲另类| 午夜福利网站1000一区二区三区| 中文字幕精品免费在线观看视频 | videos熟女内射| 观看av在线不卡| 国产成人aa在线观看| 色网站视频免费| 亚洲成人手机| 老司机影院毛片| 亚洲综合精品二区| 在线天堂最新版资源| 爱豆传媒免费全集在线观看| 91久久精品国产一区二区成人| 欧美最新免费一区二区三区| 精品卡一卡二卡四卡免费| 黄色一级大片看看| 91精品国产国语对白视频| 精品酒店卫生间| 亚洲成人手机| 久久国产精品男人的天堂亚洲 | 桃花免费在线播放| 麻豆精品久久久久久蜜桃| .国产精品久久| 99九九线精品视频在线观看视频| 国产又色又爽无遮挡免| 制服人妻中文乱码| 亚洲性久久影院| 中文字幕av电影在线播放| 欧美日本中文国产一区发布| 亚洲精品av麻豆狂野| 黄色一级大片看看| 美女cb高潮喷水在线观看| 亚洲av中文av极速乱| 纯流量卡能插随身wifi吗| 中文欧美无线码| 久久久久久人妻| 十分钟在线观看高清视频www| 欧美精品高潮呻吟av久久| 97在线视频观看| 国产伦理片在线播放av一区| 久久99热这里只频精品6学生| 成人亚洲精品一区在线观看| 丰满少妇做爰视频| 亚洲国产精品一区二区三区在线| 国产精品久久久久久精品电影小说| 亚洲一区二区三区欧美精品| 女性被躁到高潮视频| 熟女电影av网| 老司机亚洲免费影院| 亚洲av日韩在线播放| 亚洲四区av| 高清黄色对白视频在线免费看| 精品一区二区三区视频在线| 中文天堂在线官网| 日韩亚洲欧美综合| 91午夜精品亚洲一区二区三区| 少妇熟女欧美另类| freevideosex欧美| 国产乱来视频区| 欧美日韩在线观看h| 纯流量卡能插随身wifi吗| 成人黄色视频免费在线看| 久热久热在线精品观看| 国产精品一国产av| 自线自在国产av| 国产成人免费观看mmmm| 免费少妇av软件| 精品久久国产蜜桃| 99久久精品一区二区三区| 免费久久久久久久精品成人欧美视频 | 建设人人有责人人尽责人人享有的| 久久久久久久国产电影| 18在线观看网站| 一边摸一边做爽爽视频免费| 亚洲av中文av极速乱| 少妇熟女欧美另类| 日本与韩国留学比较| 另类精品久久| 婷婷色综合大香蕉| 99久久人妻综合| 亚州av有码| 精品人妻在线不人妻| 久久ye,这里只有精品| 91精品一卡2卡3卡4卡| 9色porny在线观看| 岛国毛片在线播放| 卡戴珊不雅视频在线播放| 久久久久久久久久人人人人人人| 插阴视频在线观看视频| 永久免费av网站大全| 免费黄色在线免费观看| 欧美亚洲 丝袜 人妻 在线| 欧美3d第一页| 一区在线观看完整版| 欧美 日韩 精品 国产| 又大又黄又爽视频免费| 免费日韩欧美在线观看| 欧美精品人与动牲交sv欧美| 建设人人有责人人尽责人人享有的| 国产黄色免费在线视频| 婷婷成人精品国产| 人人妻人人澡人人看| 十八禁高潮呻吟视频| 一级毛片 在线播放| 成人综合一区亚洲| 久久午夜综合久久蜜桃| 久久久久国产网址| 成人无遮挡网站| 人妻夜夜爽99麻豆av| 欧美精品亚洲一区二区| 欧美日韩视频精品一区| 97在线人人人人妻| 男女边吃奶边做爰视频| 国产在线视频一区二区| 最后的刺客免费高清国语|