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

    渦度收支與潛熱釋放對西南低渦形成的作用

    2014-09-11 01:06:30陳鵬徐海明林永輝
    大氣科學(xué)學(xué)報(bào) 2014年5期
    關(guān)鍵詞:平流渦度潛熱

    陳鵬,徐海明,林永輝

    (1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院,江蘇 南京 210044;2.中國氣象科學(xué)研究院 災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,北京 100081;3.重慶市氣象臺,重慶 401147)

    渦度收支與潛熱釋放對西南低渦形成的作用

    陳鵬1,2,3,徐海明1,林永輝2

    (1.南京信息工程大學(xué) 大氣科學(xué)學(xué)院,江蘇 南京 210044;2.中國氣象科學(xué)研究院 災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室,北京 100081;3.重慶市氣象臺,重慶 401147)

    利用中尺度WRF模式對2008年6月30日—7月1日生成于川東南地區(qū)的一個(gè)西南低渦的發(fā)生發(fā)展過程進(jìn)行了數(shù)值模擬研究。模擬結(jié)果顯示低渦首先出現(xiàn)在850 hPa上,幾個(gè)小時(shí)后700 hPa上才有低渦生成,850 hPa低渦的形成與西南低空急流有著密切的聯(lián)系。通過ω方程的診斷分析表明,渦度的水平平流項(xiàng)和輻散項(xiàng)對850 hPa低渦的形成起主要作用,而潛熱釋放對850 hPa低渦的形成作用不大;潛熱加熱是700 hPa氣流不斷輻合從而形成低渦的主要因子。干敏感性試驗(yàn)研究進(jìn)一步證實(shí)了潛熱釋放對850 hPa低渦的影響不明顯,但是會導(dǎo)致700 hPa上氣旋性的切變加強(qiáng)輻合從而形成低渦。

    西南低渦;低空急流;潛熱釋放;暴雨

    0 引言

    西南低渦是在我國青藏高原地形影響下以及在一定的環(huán)流背景形勢下產(chǎn)生的中間尺度天氣系統(tǒng)(陶詩言,1980)。吳國雄和劉還珠(1999)指出西南低渦的產(chǎn)生是地形作用的結(jié)果。Kuo et al.(1986)認(rèn)為500 hPa上東移短波槽的強(qiáng)迫動(dòng)力作用會產(chǎn)生西南低渦。沈如桂等(1983)研究發(fā)現(xiàn)西南低渦的形成是季風(fēng)槽的氣旋性渦度輸送和來自孟加拉灣的熱帶季風(fēng)的水汽在低渦區(qū)輻合所致。鄒波和陳忠明(2000)發(fā)現(xiàn)大氣邊界層頂?shù)姆瞧胶鈩?dòng)力強(qiáng)迫及500 hPa 正渦度平流的動(dòng)力作用對低渦的形成和發(fā)展有重要貢獻(xiàn)。高守亭(1987)采用定常二層模式討論較小地形及高、低層流場配置對西南低渦形成的動(dòng)力作用,指出西南低渦的形成與盆地、河谷以及其上空氣流分層有關(guān),上、下為西風(fēng)分層時(shí),低層的淺薄暖濕西風(fēng)有利于西南低渦的形成;在上、下為東、西風(fēng)分層時(shí),上層淺薄東風(fēng)有利于西南低渦的形成。王賽西(1992)從角動(dòng)量輸送方面指出低渦源地正角動(dòng)量的大量增加為西南渦的形成提供了必須的動(dòng)力,對西南渦的生成具有一定的促進(jìn)作用。陳忠明(2007)應(yīng)用動(dòng)力分析方法發(fā)現(xiàn),在大氣運(yùn)動(dòng)處于不平衡狀態(tài)下,慣性重力內(nèi)波的發(fā)展可能是西南低渦發(fā)展的一種物理機(jī)制。另外當(dāng)氣塊沿等θ面絕熱下滑而誘發(fā)的垂直渦度快速增長遠(yuǎn)遠(yuǎn)超過了下沉氣柱由于密度壓縮誘發(fā)的渦度增長值時(shí),也會導(dǎo)致了西南低渦的形成(吳國雄和劉還珠,1999)。

    實(shí)際上,大氣熱力作用(尤其是潛熱作用)在西南低渦發(fā)展中的重要地位也一直為大家所認(rèn)同(劉紅武和李國平,2008)。李國平等(1991)利用熱成風(fēng)適應(yīng)理論表明由于地面感熱加熱與暖平流作用在西南低渦源地形成較大的非熱成風(fēng)渦度,在一定的層結(jié)和尺度條件下,其適應(yīng)調(diào)整在低層形成暖性熱低壓(即暖性西南低渦)。盧敬華(1988)通過研究也認(rèn)為正的非熱成風(fēng)渦度是低渦生成的一個(gè)重要指標(biāo)。俞樟孝等(1982)研究指出低層暖平流的作用會破壞西南地區(qū)的平衡,產(chǎn)生出非熱成風(fēng)渦度,通過非熱成風(fēng)渦度引起的調(diào)整變化,以及地形作用,從而低層產(chǎn)生一個(gè)暖性低渦。對于凝結(jié)潛熱釋放的作用,一般認(rèn)為是西南渦得以維持和發(fā)展的主要機(jī)制(濮梅娟等,1989;Wang et al.,1993;Chen et al.,1997)。趙平和孫淑清(1991)經(jīng)過試驗(yàn)發(fā)現(xiàn)潛熱通過加強(qiáng)西南低渦上空高層輻散和低層輻合,使該低渦發(fā)展。郭大梅和丁治英(2005)認(rèn)為潛熱釋放與降水之間的正反饋有利于氣旋的加深發(fā)展。

    以上的研究由于所選擇的個(gè)例不同,得到的形成機(jī)理也不相同,這充分說明了西南低渦形成發(fā)展的多樣性和復(fù)雜性,對其形成機(jī)理無疑值得進(jìn)一步研究和探討。本文對2008年6月30日—7月1日生成于川東南地區(qū)的一個(gè)西南低渦的發(fā)生發(fā)展過程進(jìn)行了數(shù)值模擬研究。

    1 數(shù)值模擬

    利用WRF模式對這個(gè)西南低渦進(jìn)行數(shù)值模擬,以2008年6月30日00時(shí)(世界時(shí),下同)為初始時(shí)刻,積分到7月2日00時(shí)。模式使用粗細(xì)嵌套網(wǎng)格方案,粗細(xì)網(wǎng)格中水平分辨率分別為45 km和15 km。模擬格點(diǎn)數(shù)分別為100×112 、142×148,中心位于105°E、35°N。網(wǎng)格和次網(wǎng)格尺度降水采用Kain-Fritsch積云對流參數(shù)化方案,行星邊界層物理過程采用Yonsei University方案,側(cè)邊界條件采用NCEP 每6 h一次的1°×1°的格點(diǎn)資料。

    對比2008年6月30日12時(shí)850 hPa流場實(shí)況(圖1a)和模擬(圖1b),可以看出在川東南地區(qū)模擬流場和實(shí)況流場都有氣旋式渦旋,渦旋中心的位置大體吻合,整個(gè)流場的形勢也基本上一致,模擬效果較好。所以本文利用WRF模式的高分辨率輸出資料對此次低渦的發(fā)生、發(fā)展過程進(jìn)行研究分析。

    2 西南低渦的發(fā)生發(fā)展

    2.1 850 hPa西南低渦

    從模擬的850 hPa流場可以看到氣旋性閉合環(huán)流初見于2008年6月30日12時(shí),位于(105.6°E,28.1°N)附近。6月30日14時(shí),川東南地區(qū)低渦形成,低渦中心略向北移動(dòng)。16時(shí)低渦中心移至(105.6°E,28.3°N)附近,強(qiáng)度逐漸加強(qiáng),低渦的尺度也進(jìn)一步擴(kuò)大。隨后低渦強(qiáng)度繼續(xù)增大,其中心仍向北移動(dòng)(圖略)。

    圖1 2008年6月30日12時(shí)850 hPa流場(陰影為地形高度大于1 500 m) a.實(shí)況;b.模擬Fig.1 The 850 hPa stream fields at 1200 UTC 30 June 2008(shaded area:terrain height>1 500 m) a.observation;b.simulation

    圖2 2008年6月30日全風(fēng)速(a,c;單位:m/s)和渦度場(b,d;單位:10-4 s-1)沿105.6°E的緯度—高度剖面 a,b.12時(shí);c,d.14時(shí)Fig.2 Latitude-height vertical cross sections of (a,c)wind speed (units:m/s) and (b,d)vorticity fields (units:10-4 s-1) along 105.6°E on 30 June 2008 a,b.1200 UTC;c,d.1400 UTC

    低空急流對低渦的發(fā)展有著重要的作用(何編等,2012)。由圖2a可以看到,初始低渦形成時(shí)(6月30日12時(shí)),850 hPa上存在強(qiáng)度為16 m/s的低空急流,在急流的北面有一個(gè)渦度值為5×10-4s-1的輻合中心,位于(105.6°E,28.1°N)(圖2b),這和氣旋性閉合環(huán)流的中心吻合。6月30日14時(shí)(圖2c),低層低空急流強(qiáng)度增大為18 m/s,川東南地區(qū)低層的渦度輻合中心的強(qiáng)度也在增大,強(qiáng)度為6×10-4s-1(圖2d),此時(shí)氣旋性閉合環(huán)流由于強(qiáng)的輻合作用已經(jīng)發(fā)展成為中尺度低渦。6月30日16時(shí)低渦中心向北移動(dòng)的同時(shí),急流軸向北移動(dòng),中心強(qiáng)度增強(qiáng),相應(yīng)的低渦中心的渦度值增大(圖略)。可以看出,低空急流核前方為風(fēng)速梯度最大的地方,也就是氣流輻合最大的地方。在低空急流軸的北側(cè),有氣旋性風(fēng)速切變產(chǎn)生正渦度,為低空的輻合起到了加速的作用,這有利于低渦的形成(王歡和倪允琪,2006)。西南低空急流的發(fā)展和維持,是暴雨低渦形成和持續(xù)發(fā)展的重要?jiǎng)恿W(xué)條件(程麟生等,1995)。

    下面通過渦度收支方程(朱愛軍和潘益農(nóng),2007)分析方程中各項(xiàng)對低渦形成的作用,渦度收支方程的表達(dá)式如下:

    圖3 2008年6月30日渦度收支方程各項(xiàng)區(qū)域(104.9~105.5°E,28.9~29.5°N)平均的垂直分布(單位:10-8 s-2) a.12時(shí);b.14時(shí);c.16時(shí)Fig.3 The vorticity budget equation at (a)1200 UTC,(b)1400 UTC,(c)1600 UTC 30 June 2008:The regional average (28.9—29.5°N,104.9—105.5°E) vertical distribution of horizontal advection(long dotted line with ○)、vertical transportation(short dotted line with solid ○)、horizontal convergence and divergence(dotted line with □)、twist(attunements line with solid □) and the total vorticity balance(solid line with *)(units:10-8 s-2)

    從圖3可以看出渦度水平平流項(xiàng)(A項(xiàng))是渦旋變化的最大貢獻(xiàn)項(xiàng)。水平輻合項(xiàng)(C項(xiàng))對中低層渦度也為正貢獻(xiàn)、對高層則為負(fù)貢獻(xiàn),但是水平輻合項(xiàng)的作用相對水平平流項(xiàng)的作用要弱一些。在低渦形成過程中,垂直輸送項(xiàng)在低層對正渦度有正貢獻(xiàn),對于高層正渦度的貢獻(xiàn)基本很小。

    6月30日12時(shí)總的渦度收支分布是:低層900~700 hPa之間為正渦度收支,主要是渦度水平平流項(xiàng)的作用,垂直輸送項(xiàng)對正渦度的收支也有一定貢獻(xiàn),此時(shí)850 hPa出現(xiàn)氣旋性輻合。700~500 hPa之間,主要為負(fù)渦度收支(圖3a)。6月30日14時(shí),850 hPa正渦度收支增大為2.6×10-8s-1,水平平流項(xiàng)的值相比12時(shí)也增大0.5個(gè)單位(圖3b)。6月30日16時(shí),低層水平平流項(xiàng)(A項(xiàng))和水平輻合輻散項(xiàng)(C項(xiàng))對總渦度收支的正貢獻(xiàn)作用是最主要的。相比形成初期,中層的水平輻散加大,低層850 hPa正渦度收支的最大值繼續(xù)增大,低渦的強(qiáng)度進(jìn)一步增強(qiáng)(圖3c)。

    由渦度收支的各項(xiàng)可以知道,水平渦度平流項(xiàng)和輻散項(xiàng)均起到最主要的作用。對低層850 hPa低渦的形成發(fā)展,水平平流項(xiàng)起直接作用;水平平流最高層為負(fù),中低層為正,產(chǎn)生上升運(yùn)動(dòng),從而引起中低層輻合,對低渦發(fā)展起作用;垂直輸送項(xiàng)對中低層總渦度始終是正貢獻(xiàn),對低層低渦的發(fā)展有間接作用;而扭轉(zhuǎn)項(xiàng)主要是對總渦度發(fā)展的貢獻(xiàn)不大。

    850 hPa渦度圖(圖略)可以看到6月30日12時(shí)在川東南地區(qū)存在一條正的渦度帶,并且有一渦度中心。在正渦度帶以西有一負(fù)渦度帶。兩小時(shí)后正渦度帶仍然穩(wěn)定在川東南地區(qū),負(fù)渦度帶略有東移向正渦度帶靠近,正渦度帶的低渦中心向北緩慢移動(dòng),中心強(qiáng)度也有增強(qiáng)。6月30日15時(shí)正渦度低渦中心繼續(xù)北移,中心強(qiáng)度繼續(xù)增強(qiáng),相應(yīng)的850 hPa的低渦中心也在向北移動(dòng)。同時(shí)15時(shí)水平輻合輻散項(xiàng)作用增強(qiáng),這說明水平輻合使得移動(dòng)的低渦系統(tǒng)加強(qiáng),而低渦中心不斷北移主要是因?yàn)樗狡搅鬏斔偷淖饔谩?/p>

    模擬區(qū)域平均的視熱源Q1和視水汽匯Q2的時(shí)間—高度演變(圖略)顯示,在低層Q1和Q2都為正值,6月30日14時(shí)最大的加熱中心位于850 hPa,Q2的強(qiáng)度略強(qiáng)于Q1。這表明低層水汽輻合上升會釋放凝結(jié)潛熱,這將會導(dǎo)致低層位勢高度降低,在降水中心附近會造成額外的氣壓梯度,從而加強(qiáng)流向輻合區(qū)的中尺度低空急流(喬林等,2009),低空急流的增強(qiáng)進(jìn)而引起低層輻合的加強(qiáng),低渦進(jìn)一步加強(qiáng)。

    通過以上的分析可以看出,氣流繞過高原后通過水平平流的輸送作用在川東南地區(qū)匯合,而且低空急流的發(fā)展正好對應(yīng)著川東南地區(qū)850 hPa中尺度低渦的發(fā)生,所以此次個(gè)例中低層850 hPa低渦的形成與西南低空急流有著密切的關(guān)系,而潛熱釋放對850 hPa低渦的后續(xù)發(fā)展和維持有著較為重要的作用。

    2.2 700 hPa西南低渦

    700 hPa上閉合環(huán)流的出現(xiàn)晚于850 hPa的渦旋,6月30日21時(shí)700 hPa才出現(xiàn)低渦(圖4a)。30日22時(shí)低渦強(qiáng)度增強(qiáng),中西北移(圖4b)。6月30日18時(shí)、19時(shí)700 hPa上川東南地區(qū)位勢高度逐漸減弱,在6月30日19時(shí)形成低值中心,中心大致在(105.1°E,29.1°N),風(fēng)場也在不斷的向低值中心輻合(圖略)。在低值中心形成之前川東南地區(qū)就已經(jīng)有對流降水的發(fā)生,但降水強(qiáng)度在不斷減小,降水一直持續(xù)到6月30日22時(shí),降水中心和700 hPa低渦中心沒有重合,位置相對低渦中心偏南。

    假相當(dāng)位溫θse、相對濕度、風(fēng)矢量的高度—緯度剖面顯示在6月30日18時(shí)一直到21時(shí)在800~600 hPa間始終存在著暖濕空氣帶(圖5),這和降水過程中潛熱釋放加熱有一定的關(guān)系。由于暖濕帶的存在,700 hPa上先是有弱低渦生成,此后氣體濕度的增大,低渦中心的輻合強(qiáng)度也增大。相比之下,在800 hPa以下為干冷空氣帶,這主要是由于空氣在下降時(shí)的凝結(jié)冷卻所致,從而在川東南地區(qū)形成一個(gè)冷池,而且冷卻作用會加強(qiáng)空氣的下沉運(yùn)動(dòng)以及中層的輻合(Lafoe and Moncrieff,1989),所以700 hPa的潛熱加熱和800 hPa以下的凝結(jié)冷卻對700 hPa的低渦形成都有影響。

    圖4 2008年6月30日700 hPa流場(陰影為地形高度大于3 000 m) a.21時(shí);b.22時(shí)Fig.4 The 700 hPa stream fields at (a)2100 UTC and (b) 2200 UTC 30 June 2008(shaded area:terrain height>3 000 m)

    圖5 2008年6月30日θse(實(shí)線;單位:K)、相對濕度(虛線;單位:%)、風(fēng)矢量(單位:10-1 m/s)沿105.2°E的緯度—高度剖面 a.18時(shí);b.19時(shí)Fig.5 The latitude-height vertical cross sections of θse (solid line;units:K),relative humidity (dotted line;units:%),wind vector(units:10-1 m/s) along 105.2°E at (a)1800 UTC and (b)1900 UTC 30 June 2008

    圖5中還可以看到在低渦中心上部,上升運(yùn)動(dòng)引起輻合增強(qiáng),而在低渦中心下部的下沉運(yùn)動(dòng)是較弱的,這使得潛熱加熱對低渦產(chǎn)生的作用強(qiáng)于低層凝結(jié)冷卻對低渦產(chǎn)生的作用,潛熱釋放應(yīng)該對氣流的上升以及低渦的生成和加強(qiáng)起到很重要的作用。另一方面低渦產(chǎn)生后也通過輻合上升運(yùn)動(dòng)向高層帶去更多暖濕空氣從而加強(qiáng)潛熱加熱作用,潛熱加熱再作用低渦的加強(qiáng),形成反饋?zhàn)饔谩?/p>

    圖6 2008年6月30日18時(shí)沿105.2°E溫度擾動(dòng)(a,單位:℃;等值線間隔為0.5 ℃)、擾動(dòng)θse(b,單位:K;等值線間隔為0.5 K)、擾動(dòng)位渦(c,單位:PVU;等值線間隔為0.5 PVU)垂直剖面Fig.6 The vertical cross sections of (a)temperature perturbation (units:℃;contour lines spacing of 0.5 ℃),(b)θse perturbation(units:K;contour lines spacing of 0.5 K),(c)PV perturbation vertical section(units:PVU;contour lines spacing of 0.5 PVU) at 1800 UTC 30 June 2008

    Hoskins and Berrisford(1988)指出位渦是研究中緯度現(xiàn)象的一種簡單而又簡明的方法。丁治英等(2010)發(fā)現(xiàn)低層的正位渦通過上升運(yùn)動(dòng)向上傳遞,導(dǎo)致了高層位渦正異常,高層位渦的正異常又可導(dǎo)致低層的氣旋性渦度進(jìn)一步加大并使降水加大。Zhang and Fritsch(1988)通過集合預(yù)報(bào)的分析發(fā)現(xiàn)對流層中上層潛熱釋放的產(chǎn)生以及位渦值的增大會增強(qiáng)中層渦旋。而且一些研究結(jié)果表明潛熱釋放使得位渦再分配,從而使得負(fù)位渦異常區(qū)和與之相關(guān)的反氣旋環(huán)流在最大加熱層上部形成。相反,正位渦異常區(qū)與氣旋環(huán)流在最大加熱層下部發(fā)生發(fā)展(Thorpe,1985;Haynes and McIntyre,1987;Raymond and Jiang,1990)。

    沿29.1°N的經(jīng)向風(fēng)和位渦的經(jīng)度—高度剖面(圖略)顯示,在700 hPa低渦中心的位渦值由6月30日18時(shí)的1 PVU(1 PVU=10-6m2·K·s-1·kg-1)增大到30日20時(shí)的2.5 PVU,并且低渦中心的渦度值也由1.5×10-4s-1變?yōu)?.0×10-4s-1,垂直上升運(yùn)動(dòng)也在增強(qiáng)。這表明在潛熱釋放過程,正位渦異常大值對應(yīng)暖濕加熱區(qū)域,對低渦的形成和發(fā)展很著重要的影響。

    此外也可通過擾動(dòng)量的變化來分析潛熱釋放對低層低渦形成的影響。6月30日18時(shí),從擾動(dòng)量的垂直分布圖6a中可以看到,700 hPa為相對暖的區(qū)域,是θse高值區(qū),這主要是暖濕氣流上升過程中引起的潛熱釋放的結(jié)果;700 hPa以下為相對冷區(qū),是θse低值區(qū),這是凝結(jié)冷卻使得低層溫度降低的結(jié)果。位渦垂直剖面顯示,700 hPa為位渦大值區(qū)域,6月30日19時(shí)(圖7),川東南地區(qū)700 hPa有暖區(qū)中心形成,此時(shí)低渦也輻合生成,在低渦中心也有位渦大值中心生成。從渦度廓線也可以看到700 hPa上渦度的值也在略微的隨時(shí)間增加,這說明潛熱釋放引起了對流層低層位渦的增加,導(dǎo)致了渦度的增加。而且氣流在川東南地區(qū)低層的風(fēng)速隨時(shí)間也在不斷的增大,低層的輻合最大值由1.5×10-4s-1變?yōu)?.0×10-4s-1。低層的強(qiáng)的輻合導(dǎo)致強(qiáng)的上升運(yùn)動(dòng),暖濕氣流在上升過程中釋放潛熱導(dǎo)致了位渦的增加,對應(yīng)的正渦度也在增加,流場的氣旋性彎曲也顯著增強(qiáng),從而逐漸形成氣旋性的渦旋(圖7b)。所以潛熱釋放是導(dǎo)致低層700 hPa渦度增加從而不斷輻合形成低渦的主要因子。

    圖7 2008年6月30日19時(shí)沿105.2°E溫度擾動(dòng)(a,單位:℃;等值線間隔為0.5 ℃)、擾動(dòng)θse(b,單位:K;等值線間隔為0.5 K)、擾動(dòng)位渦(c,單位:PVU;等值線間隔為0.5 PVU)垂直剖面Fig.7 The vertical cross sections of (a)temperature perturbation (units:℃;contour lines spacing of 0.5 ℃),(b)θse perturbation(units:K;contour line spacing of 0.5 K),(c)PV perturbation vertical section(units:PVU;contour line spacing of 0.5 PVU) at 1900 UTC 30 June 2008

    圖8 2008年6月30日18時(shí)區(qū)域平均(104.9~105.5°E,28.9~29.5°N)凝結(jié)加熱率(實(shí)線;單位:K·d-1)、散度(空心○虛線;單位:10-4 s-1)、渦度(點(diǎn)劃線;單位:10-4 s-1)和垂直速度(實(shí)心○虛線;單位:m/s)的垂直廓線及ω方程1、2、3三項(xiàng)垂直運(yùn)動(dòng)分量(b)Fig.8 (a)The regional average(28.9—29.5°N,104.9—105.5°E) vertical profile of condensation heating rate(solid line;units:K·d-1),the divergence(dashed lines;units:10-4 s-1),the vorticity(attunements line;units:10-4 s-1) and vertical velocity(solid line;units:m/s) and (b)vertical movement of the three components are calculated respectively in ω equation at 1800 UTC 30 June 2008

    再利用潛熱加熱率公式(Kuo,1974),對低渦發(fā)生發(fā)展進(jìn)行診斷分析。6月30日18時(shí)(圖8a)區(qū)域平均的凝結(jié)加熱率、散度、渦度和垂直速度的垂直廓線顯示,最強(qiáng)的上升運(yùn)動(dòng)處于中層500 hPa,中層處于輻合狀態(tài),凝結(jié)潛熱加熱率最大值位于800 hPa。利用考慮了非絕熱加熱的ω方程,進(jìn)一步診斷各項(xiàng)物理因子對垂直運(yùn)動(dòng)的影響。其中:等號右邊第一項(xiàng)為渦度平流垂直變化項(xiàng);第二項(xiàng)為溫度平流拉普拉斯項(xiàng);第三項(xiàng)為非絕熱加熱項(xiàng),分別稱為1、2、3項(xiàng),在對各項(xiàng)計(jì)算過程中采用的都是一種近似的計(jì)算。圖8b為6月30日18時(shí)各項(xiàng)對垂直速度的影響,可以看出在500 hPa沒有明顯的溫度平流,但是有明顯的正渦度平流。在對流層低層,潛熱加熱對垂直運(yùn)動(dòng)的貢獻(xiàn)較為重要,尤其是低層700 hPa潛熱加熱造成的上升運(yùn)動(dòng)超過了低層暖平流的作用。此時(shí)在800 hPa以下有冷平流的作用,造成下沉運(yùn)動(dòng)。

    圖9 2008年6月30日19時(shí)區(qū)域平均(104.9~105.5°E,28.9~29.5°N)凝結(jié)加熱率(實(shí)線;單位:K·d-1)、散度(空心○虛線;單位:10-4 s-1)、渦度(點(diǎn)劃線;單位:10-4 s-1)和垂直速度(實(shí)心○虛線;單位:m/s)的垂直廓線及ω方程1、2、3三項(xiàng)垂直運(yùn)動(dòng)分量(b)Fig.9 (a)The regional average(28.9—29.5°N,104.9—105.5°E) vertical profile of condensation heating rate(solid line;units:K·d-1),the divergence(dashed lines;units:10-4 s-1),the vorticity(attunements line;units:10-4 s-1) and vertical velocity(solid line;units:m/s) and (b)vertical movement of the three component are calculated respectively in ω equation at 1900 UTC 30 June 2008

    圖9a可見,6月30日19時(shí)渦度和散度廓線表明最強(qiáng)的輻合仍然位于對流層低層,800~600 hPa之間,輻合最大值為2×10-4s-1,相比前一個(gè)小時(shí)變化不大。對流層中層以上400 hPa存在一個(gè)較大的輻散層,上升速度中心最大值為1.5 m/s,高度在500 hPa左右。散度的結(jié)構(gòu)和渦度結(jié)構(gòu)近于反位相,這和朱愛軍和潘益農(nóng)(2007)研究是一致的。對流層低層有在兩個(gè)正渦度中心,一個(gè)在850 hPa,另一個(gè)在700 hPa。此時(shí)低層700 hPa凝結(jié)加熱率比18時(shí)有所增大,但是仍然為負(fù)。ω方程診斷的各強(qiáng)迫項(xiàng)對應(yīng)的垂直速度表明(圖9b),凝結(jié)加熱項(xiàng)強(qiáng)迫上升運(yùn)動(dòng)在對流層低層700 hPa上的作用是最重要的,最大值達(dá)到2 m/s,相比前一個(gè)小時(shí)有所增大。這表明降水導(dǎo)致的凝結(jié)加熱對700 hPa低渦輻合起很大的作用。此時(shí)渦度平流垂直變化對上升運(yùn)動(dòng)的作用為負(fù)貢獻(xiàn),但是溫度平流在低層700 hPa上表現(xiàn)為正貢獻(xiàn)。而且800 hPa以下仍然有冷平流導(dǎo)致一定的下沉運(yùn)動(dòng),在對流層中層由于湍流混合作用較強(qiáng),溫度平流較弱,對垂直運(yùn)動(dòng)基本沒有貢獻(xiàn)(喬林等,2009)。

    6月30日20時(shí),850 hPa的渦度值減小,而700 hPa低渦中心的渦度值略有增大,低層潛熱加熱率繼續(xù)增大,約為1 K/h,同時(shí)潛熱加熱對低層輻合和上升運(yùn)動(dòng)的影響進(jìn)一步增強(qiáng)(圖略),在高層溫度平流對上升運(yùn)動(dòng)的影響很小,但是渦度平流對垂直上升運(yùn)動(dòng)有著很強(qiáng)的貢獻(xiàn)。

    通過上面的分析可以表明,隨著低渦的形成和成熟,低渦的下部出現(xiàn)與下沉運(yùn)動(dòng)相伴的輻散結(jié)構(gòu)。這種結(jié)構(gòu)與Maddox(1983)及Menard and Fritsch(1989)提出的成熟時(shí)的中尺度對流輻合體相類似,即對流層高層與低層都為輻散外流,而對流層中層出現(xiàn)氣旋性結(jié)構(gòu)。

    3 數(shù)值敏感試驗(yàn)

    通過干敏感性試驗(yàn)進(jìn)一步研究潛熱釋放對700 hPa低渦的影響,它與全因子的不同之處就是模擬當(dāng)中沒有降水凝結(jié)潛熱釋放的作用,因此它與全因子試驗(yàn)的差異,可以認(rèn)為是由凝結(jié)潛熱過程所致。

    從6月30日全因子模擬850 hPa流場和干試驗(yàn)?zāi)M850 hPa流場對比(圖略)可以看到,去除了潛熱釋放后在川東南地區(qū)仍有中尺度的低渦生成,低渦中心和尺度基本上和全因子模擬的結(jié)果相近,這表明潛熱釋放對850 hPa上低渦的影響可以忽略。而從700 hPa流場可以看到,6月30日21時(shí)在四川盆地中部形成一個(gè)低渦(圖4a)。30日22時(shí)該中尺度渦旋逐漸加強(qiáng),低渦中心略向北移動(dòng)(圖4c)。如果將模式中的潛熱釋放關(guān)掉后,可以看到在6月30日21時(shí)、22時(shí)川東南附近的地區(qū)沒有中尺度渦旋形成(圖略),這說明在700 hPa上降水發(fā)生過程中潛熱的釋放會導(dǎo)致氣旋性的切變加強(qiáng)輻合從而形成低渦。

    圖10 2008年6月30日18時(shí)(a,b)、19(c,d)區(qū)域平均(104.9~105.5°E,28.9~29.5°N)的位勢高度(點(diǎn)線;單位:dagpm)、溫度(實(shí)線;單位℃)的垂直廓線(a,c)和低渦中心(105.2°E,29.1°N)的渦度(點(diǎn)線;單位:10-4 s-1)和散度(實(shí)線;單位:10-4 s-1)(b,d)的垂直廓線Fig.10 The regional average (28.9—29.5°N,104.9—105.5°E) (a,c)height (point line;units:dagpm),temperature (solid line;units:℃),vertical profile and (b,d)low vortex center (29.1°N,105.2°E) vorticity (point line;line;units:10-4 s-1),divergence(solid line;units:10-4 s-1) at (a,b)1800 UTC and (c,d)1900 UTC 30 June 2008

    為了更進(jìn)一步通過數(shù)值試驗(yàn)來反映潛熱釋放對700 hPa低渦的作用,利用全因子模擬和去掉潛熱釋放后模擬的各個(gè)變量的差值進(jìn)行分析。首先全因子模擬和去掉潛熱釋放后模擬的溫度場和高度場之差的區(qū)域平均的垂直分布(圖10a、10c)中可以看出,潛熱釋放使得800~400 hPa之間的大氣都減壓,而400 hPa以上是大氣加壓,并且高層的加壓比低層的減壓大。從溫度垂直廓線中看以看到凝結(jié)潛熱加熱了800~400 hPa之間較深厚的大氣層,最大加熱中心位于400 hPa附近。從兩個(gè)實(shí)驗(yàn)的渦度和散度的差值顯示,在6月30日18時(shí)低渦形成初期,凝結(jié)潛熱對暴雨區(qū)正渦度的貢獻(xiàn)較弱,700 hPa正渦度差值中心值為1.0×10-4s-1(圖10b)。6月30日19時(shí)低渦形成,700 hPa正渦度差值明顯增大,中心強(qiáng)度達(dá)到2.6×10-4s-1(圖10d),700 hPa以下有很明顯的低層正渦度輻合場,700 hPa以上為負(fù)渦度場。通過上面的分析可知,降水過程中凝結(jié)潛熱的釋放加熱了大氣層,由熱成風(fēng)的關(guān)系,暴雨區(qū)北側(cè)高空西北氣流加強(qiáng),空氣質(zhì)量大量流出,即高空輻散加強(qiáng),同時(shí)低層出現(xiàn)減壓和輻合(張曉芳和陸漢城,2006),這使得700 hPa上氣旋性切邊加強(qiáng)輻合,從而轉(zhuǎn)變成氣旋性渦旋。所以潛熱加熱直接增暖了中低層大氣,增強(qiáng)上升運(yùn)動(dòng),促使700 hPa低渦的形成。這和趙平和孫淑清(1991)潛熱加熱不影響西南低渦的形成的研究結(jié)果不同。

    通過上面的分析,可以認(rèn)為降水過程中由于水汽充足,低層強(qiáng)輻合引起上升運(yùn)動(dòng)的同時(shí),引起水汽的凝結(jié)。由于凝結(jié)潛熱的釋放,增強(qiáng)了川東南地區(qū)低層輻合、高層輻散以及中低層正渦度的發(fā)展,所以水汽凝結(jié)、釋放潛熱是700 hPa低渦形成的主要原因。

    4 結(jié)論與討論

    本文運(yùn)用WRF模式模擬了2008年6月30日到7月1日生成于川東南地區(qū)的一個(gè)西南低渦的發(fā)生發(fā)展過程。數(shù)值模擬診斷分析表明:

    1)閉合性的氣旋中心首先在850 hPa上生成,850 hPa低渦的形成與西南低空急流有著密切的關(guān)系,潛熱釋放對850 hPa低渦的后續(xù)發(fā)展和維持有著較為重要的作用。

    2)渦度收支的各項(xiàng)表明水平渦度平流項(xiàng)和輻散項(xiàng)均起到最主要的作用。對低層850 hPa低渦的形成發(fā)展,水平平流項(xiàng)起直接作用;垂直平流項(xiàng)對中低層總渦度始終是正貢獻(xiàn),對低層低渦的發(fā)展有間接作用;而扭轉(zhuǎn)項(xiàng)主要是對總渦度發(fā)展的貢獻(xiàn)不大。

    3)通過凝結(jié)潛熱率公式和非絕熱加熱的ω方程表明,潛熱加熱對700 hPa流場的輻合有著重要的影響,能夠引起該層渦旋不斷輻合從而形成低渦。

    4)干敏感試驗(yàn)表明潛熱釋放對850 hPa低渦的生成可以忽略,但在700 hPa上會導(dǎo)致氣旋性的切變加強(qiáng)輻合從而形成低渦,潛熱釋放對700 hPa低渦的生成起主要是作用。

    通過本例的研究分析,對西南渦形成過程有了一些初淺的認(rèn)識,但是對西南地區(qū)來講,由于受到青藏高原大地形的影響,西南渦的形成原因本身就是一個(gè)難點(diǎn)。因此還應(yīng)該針對多個(gè)個(gè)例開展深入分析,進(jìn)一步深入研究西南渦形成的機(jī)理和特點(diǎn)。

    陳忠明.2007.暴雨激發(fā)和維持的正、斜壓強(qiáng)迫機(jī)制的理論研究[J].大氣科學(xué),31(2):292-297.

    程麟生,彭新東,馬艷.1995.“91.7”江淮暴雨低渦發(fā)展結(jié)構(gòu)和演變的中尺度數(shù)值模擬[J].高原氣象,14(3):270-280.

    丁治英,羅靜,沈新勇.2010.2008年6月20—21日一次β中尺度切變線、低渦降水機(jī)制研究[J].大氣科學(xué)學(xué)報(bào),33(6):657-666.

    高守亭.1987.流場配置及地形對西南低渦形成的動(dòng)力作用[J].大氣科學(xué),11(3):263-271.

    郭大梅,丁治英.2005.潛熱和慣性不穩(wěn)定與氣旋爆發(fā)性發(fā)展關(guān)系的數(shù)值研究[J].大氣科學(xué)學(xué)報(bào),28(6):794-800.

    何編,孫照渤,李忠賢.2012.一次華南持續(xù)性暴雨的動(dòng)力診斷分析和數(shù)值模擬[J].大氣科學(xué)學(xué)報(bào),35(4):466-476.

    孔期,Ghulam R,趙思雄.2005.一次引發(fā)南亞大暴雨的季風(fēng)低壓結(jié)構(gòu)、渦度與水汽收支分析[J].氣候與環(huán)境研究,10(3):526-543.

    李國平,萬軍,盧敬華.1991.暖性西南低渦生成的一種可能機(jī)制[J].應(yīng)用氣象學(xué)報(bào),2(1):91-99.

    劉紅武,李國平.2008.近三十年西南低渦研究的回顧與展望[J].高原山地氣象研究,28(2):68-72.

    盧敬華.1988.利用熱成風(fēng)適應(yīng)原理對暖性西南低渦生成機(jī)制的再分析[J].高原氣象,7(4):345-356.

    濮梅娟,劉富明,沈如金.1989.一次夏季西南低渦形成機(jī)理的數(shù)值試驗(yàn)[J].高原氣象,8(4):322-330.

    喬林,陳濤,路秀娟.2009.黔西南一次中尺度暴雨的數(shù)值模擬診斷研究[J].大氣科學(xué),33(3):350-358.

    沈如桂,林新彬,夏志強(qiáng),等.1983.印度季風(fēng)槽的活動(dòng)對我國西南低渦形成及發(fā)展的作用[J].中山大學(xué)學(xué)報(bào):自然科學(xué)版,(2):65-72.

    陶詩言.1980.中國之暴雨[M].北京:科學(xué)出版社.

    王歡,倪允琪.2006.2003年淮河汛期一次中尺度強(qiáng)暴雨過程的診斷分析和數(shù)值模擬研究[J].氣象學(xué)報(bào),64(6):735-742.

    王賽西.1992.西南低渦形成的氣候特征與角動(dòng)量輸送的關(guān)系[J].高原氣象,11(2):144-151.

    吳國雄,劉還珠.1999.全型垂直渦度傾向方程和傾斜渦度發(fā)展[J].氣象學(xué)報(bào),57(1):2-16.

    俞樟孝,翟國慶,楊厥正.1982.大氣中的不平衡狀態(tài)與西南低渦發(fā)生、發(fā)展的關(guān)系[J].杭州大學(xué)學(xué)報(bào):自然科學(xué)版,9(3):340-348.

    張曉芳,陸漢城.2006.一次梅雨鋒暴雨過程的數(shù)值模擬及機(jī)理分析[J].熱帶氣象學(xué)報(bào),22(6):626-631.

    趙平,孫淑清.1991.一次西南低渦形成過程的數(shù)值試驗(yàn)和診斷(一)——地形動(dòng)力作用和潛熱作用對西南低渦影響的數(shù)值試驗(yàn)對比分析[J].大氣科學(xué),15(6):46-52.

    朱愛軍,潘益農(nóng).2007.中國東部地區(qū)一個(gè)中尺度對流渦旋的渦度收支分析[J].南京大學(xué)學(xué)報(bào):自然科學(xué)版,43(7):260-269.

    鄒波,陳忠明.2000.一次西南低渦發(fā)生發(fā)展的中尺度診斷[J].高原氣象,19(2):142-149.

    Chen Y L,Chen X A,Chen S,et al.1997.A numerical study of the low-level jet during TAMEX IOP5[J].Mon Wea Rev,125:2583-2604.

    Haynes P H,McIntyre M E.1987.On the evolution of vorticity and potential vorticity in the presence of diabatic heating and frictional or other forces[J].J Atmos Sci,44:828-841.

    Hoskins P,Berrisford.1988.A potential vorticity perspective of the storm of 15—16 October 1987[J].Weather,43:122-129.

    Kuo H L.1974.Further studies of the parameterization of the influence of cumulus convection on large-scale flow[J].J Atmos Sci,31:1232-1240.

    Kuo Y H,Cheng L,Anthes R A.1986.Mesoscale analyses of Sichuan flood catastrophe 11—15 July 1981[J].Mon Wea Rev,114:1984-2003.

    Lafoe J P,Moncrieff M W.1989.A numerical investigation of the organization and interaction of the convective and stratiform regions of tropical squall lines[J].J Atmos Sci,46:521-544.

    Maddox R A.1983.Large scale meteorological conditions associated with midlatitude,mesoscale convective complexes[J].Mon Wea Rev,111:1475-1493.

    Menard R D,Fritsch J M.1989.A mesoscale convective complex generated inertially stable warm core vortex[J].Mon Wea Rev,117:1237-1261.

    Raymond D J,Jiang H.1990.A theory for long-lived mesoscale convective systems[J].J Atmos Sci,47:3067-3077.

    Thorpe A J.1985.Diagnosis of balanced vortex structure using potential vorticity[J].J Atmos Sci,42:397-406.

    Wang W,Kuo Y H,Warner T T.1993.A diabatically driven mesoscale vortex in the lee of the Tibetan Plateau[J].Mon Wea Rev,121:2542-2561.

    Wu G X,Chen S J.1985.The effect of mechanical forcing on the formation of a mesoscale vortex[J].Quart J Roy Meteor Soc,111:1049-1070.

    Zhang D L.1992.The formation of a cooling induced mesovortex in the trailing stratiform region of a midlatitude squall Line[J].Mon Wea Rev,120:2763-2785.

    Zhang D L,Fritsch M J.1988.A numerical simulation of a convectively generated inertially stable,warm-core extratropical mesovortex over land.Part Ⅰ:Structure and evolution[J].Mon Wea Rev,116:2660-2687.

    (責(zé)任編輯:劉菲)

    Theinfluenceoflow-leveljetandlatentheatreleaseontheformationofasouthwestvortex

    CHEN Peng1,2,3,XU Hai-ming1,LIN Yong-hui2

    (1.School of Atmospheric Science,NUIST,Nanjing 210044,China;2.State Key Laboratory of Severe Weather,Chinese Academy of Meteorological Sciences,Beijing 100081,China;3.Chongqing Meteorological Observatory,Chongqing 401147,China)

    Using mesoscale simulation WRF model,the formation and development of a southwest vortex occurred over Sichuan province on June 30 and July 1,2008 were investigated.The results indicated that the southwest vortex primarily emerged at 850 hPa and then formed at 700 hPa after several hours.The formation of the southwest vortex at 850 hPa was directly connected with low-level jet.The diagnostic analysis ofωequation showed that the horizontal advection of vorticity and divergence mainly influenced the formation of the southwest vortex at 850 hPa while latent heat release was not important to its genesis.The results also indicated that the latent heating was the key factor for the formation of the southwest vortex at 700 hPa.The dry sensitive experiment further confirmed that latent heat release had little influence on the southwest vortex at 850 hPa,but it caused strong convergence of cyclonic shear circulation and finally formed the vortex at 700 hPa.

    southwest vortex;low-level jet;latent heat release;heavy rain

    2012-10-16;改回日期2012-12-17

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目(2012CB417202);國家自然科學(xué)基金資助項(xiàng)目(40745026;40875035);公益性行業(yè)(氣象)科研專項(xiàng)(GYHYQX2007-6-20);中國氣象科學(xué)研究院基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2007Y002);中國氣象科學(xué)研究院災(zāi)害天氣國家重點(diǎn)實(shí)驗(yàn)室自主研究專項(xiàng)(2008LASWZ103)

    陳鵬,碩士,研究方向?yàn)橹行〕叨葦?shù)值模擬及天氣分析,chenpeng8495@126.com.

    10.13878/j.cnki.dqkxxb.20121016001.

    1674-7097(2014)05-0575-10

    P433

    A

    10.13878/j.cnki.dqkxxb.20121016001

    陳鵬,徐海明,林永輝.2014.渦度收支與潛熱釋放對西南低渦形成的作用[J].大氣科學(xué)學(xué)報(bào),37(5):575-584.

    Chen Peng,Xu Hai-ming,Lin Yong-hui.2014.The influence of low-level jet and latent heat release on the formation of a southwest vortex[J].Trans Atmos Sci,37(5):575-584.(in Chinese)

    猜你喜歡
    平流渦度潛熱
    2007—2020年威海機(jī)場4—9月平流低云的變化特征及影響因子分析
    南海相對渦度的時(shí)空變化特征分析*
    Effect of moxibustion combined with acupoint application on enteral nutrition tolerance in patients with severe acute pancreatitis
    斜壓渦度的變化與臺風(fēng)暴雨的關(guān)系研究
    基于衛(wèi)星和Argo觀測的阿拉伯海中北部海表鹽度季節(jié)和年際變化
    工業(yè)革命時(shí)期蒸汽動(dòng)力的應(yīng)用與熱力學(xué)理論的關(guān)系
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    荊州市一次局地濃霧天氣特征分析
    浦東機(jī)場一次低云低能見度天氣氣象服務(wù)總結(jié)
    青藏高原東部夏季降水凝結(jié)潛熱變化特征分析
    国产av一区二区精品久久| 国产精品久久电影中文字幕| e午夜精品久久久久久久| av片东京热男人的天堂| x7x7x7水蜜桃| 成人精品一区二区免费| 精品少妇一区二区三区视频日本电影| 999精品在线视频| 每晚都被弄得嗷嗷叫到高潮| 嫩草影视91久久| 久久精品国产清高在天天线| 久久人妻福利社区极品人妻图片| 国产精华一区二区三区| 国产精品亚洲av一区麻豆| 国产亚洲欧美精品永久| 国产精品野战在线观看| 国语自产精品视频在线第100页| 国产亚洲精品久久久久久毛片| 大型黄色视频在线免费观看| 色综合亚洲欧美另类图片| 变态另类丝袜制服| 免费搜索国产男女视频| 操出白浆在线播放| 国产精品 欧美亚洲| 久久人妻av系列| 色综合婷婷激情| 国产精品国产高清国产av| 久热爱精品视频在线9| 亚洲精品中文字幕在线视频| 亚洲第一av免费看| 日韩三级视频一区二区三区| 亚洲精品色激情综合| av在线天堂中文字幕| 色精品久久人妻99蜜桃| 桃色一区二区三区在线观看| 国产精品电影一区二区三区| 亚洲五月色婷婷综合| 日韩精品青青久久久久久| 在线观看一区二区三区| 琪琪午夜伦伦电影理论片6080| 一夜夜www| 91国产中文字幕| 桃色一区二区三区在线观看| 中文字幕最新亚洲高清| 欧美一级毛片孕妇| 亚洲精品中文字幕在线视频| 黑人欧美特级aaaaaa片| 99热这里只有精品一区 | 国产av在哪里看| 一区二区三区高清视频在线| 国产91精品成人一区二区三区| 男女床上黄色一级片免费看| 亚洲精品一区av在线观看| 国产精品一区二区三区四区久久 | 亚洲色图 男人天堂 中文字幕| 中文字幕精品亚洲无线码一区 | 欧美精品啪啪一区二区三区| 色哟哟哟哟哟哟| 日本五十路高清| 亚洲精品久久成人aⅴ小说| 天堂√8在线中文| 国产黄色小视频在线观看| 国产成人影院久久av| 午夜激情福利司机影院| 一本精品99久久精品77| 少妇熟女aⅴ在线视频| 久久久久久免费高清国产稀缺| 国产成人精品久久二区二区免费| 欧美成人一区二区免费高清观看 | av福利片在线| 在线播放国产精品三级| 色播在线永久视频| 午夜福利18| 亚洲国产欧洲综合997久久, | 一级a爱视频在线免费观看| 91字幕亚洲| 国产av一区二区精品久久| 国内毛片毛片毛片毛片毛片| 黄色片一级片一级黄色片| 精品久久久久久久久久久久久 | 午夜福利一区二区在线看| 天天添夜夜摸| 午夜激情福利司机影院| 日日爽夜夜爽网站| 国产视频一区二区在线看| 无遮挡黄片免费观看| 久久香蕉国产精品| 嫁个100分男人电影在线观看| 制服人妻中文乱码| 精品不卡国产一区二区三区| 中文字幕人成人乱码亚洲影| 久久人妻福利社区极品人妻图片| 成年免费大片在线观看| 国产久久久一区二区三区| 亚洲中文日韩欧美视频| 免费在线观看影片大全网站| 天堂动漫精品| 久久九九热精品免费| 91麻豆av在线| 欧美 亚洲 国产 日韩一| 久久精品人妻少妇| 99在线人妻在线中文字幕| 亚洲精华国产精华精| 成人国语在线视频| 黑人操中国人逼视频| 麻豆成人av在线观看| 一本大道久久a久久精品| 欧美+亚洲+日韩+国产| 精品人妻1区二区| 国产免费av片在线观看野外av| 一区二区三区激情视频| 亚洲午夜精品一区,二区,三区| 亚洲一区二区三区不卡视频| 久久国产精品影院| 国产爱豆传媒在线观看 | 国产精华一区二区三区| 美女 人体艺术 gogo| 国产精品永久免费网站| 丰满人妻熟妇乱又伦精品不卡| 久久伊人香网站| 亚洲第一电影网av| 日韩大尺度精品在线看网址| 美国免费a级毛片| 波多野结衣巨乳人妻| 亚洲熟女毛片儿| 黄色 视频免费看| 亚洲成人国产一区在线观看| 久久亚洲真实| 国产片内射在线| 免费无遮挡裸体视频| 色在线成人网| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 日韩国内少妇激情av| 久久久久亚洲av毛片大全| 亚洲av美国av| 免费搜索国产男女视频| 久久午夜综合久久蜜桃| 精品国产乱子伦一区二区三区| 欧美一级毛片孕妇| 老汉色av国产亚洲站长工具| 午夜福利一区二区在线看| 99久久无色码亚洲精品果冻| 久久久精品国产亚洲av高清涩受| 99精品在免费线老司机午夜| 亚洲一区二区三区色噜噜| 少妇被粗大的猛进出69影院| 久久精品91无色码中文字幕| 丰满人妻熟妇乱又伦精品不卡| 久久久久国产精品人妻aⅴ院| 日韩欧美国产一区二区入口| a级毛片a级免费在线| 成人特级黄色片久久久久久久| 成人特级黄色片久久久久久久| 午夜视频精品福利| 亚洲色图 男人天堂 中文字幕| 天堂影院成人在线观看| 中文字幕精品免费在线观看视频| 人人妻,人人澡人人爽秒播| 黑丝袜美女国产一区| 国内精品久久久久久久电影| 久久青草综合色| 身体一侧抽搐| 一级a爱片免费观看的视频| АⅤ资源中文在线天堂| 精品久久久久久久毛片微露脸| av在线天堂中文字幕| 丰满的人妻完整版| www日本在线高清视频| 女性被躁到高潮视频| 麻豆一二三区av精品| 99热这里只有精品一区 | 国产乱人伦免费视频| 欧美国产日韩亚洲一区| 亚洲一区高清亚洲精品| 久久久久久亚洲精品国产蜜桃av| 午夜激情福利司机影院| 亚洲人成伊人成综合网2020| bbb黄色大片| 啪啪无遮挡十八禁网站| cao死你这个sao货| 久久久国产精品麻豆| av天堂在线播放| 成年版毛片免费区| 丁香六月欧美| 日韩一卡2卡3卡4卡2021年| 国产欧美日韩精品亚洲av| av中文乱码字幕在线| 亚洲专区中文字幕在线| 欧美激情久久久久久爽电影| 国产精品日韩av在线免费观看| 亚洲av成人一区二区三| 性欧美人与动物交配| 精品熟女少妇八av免费久了| 最近最新免费中文字幕在线| 国产国语露脸激情在线看| 女人被狂操c到高潮| 欧美日韩亚洲国产一区二区在线观看| 91麻豆av在线| 欧美成人免费av一区二区三区| 非洲黑人性xxxx精品又粗又长| 十八禁人妻一区二区| 亚洲性夜色夜夜综合| 国产一区二区激情短视频| 久久久久久免费高清国产稀缺| 国产伦在线观看视频一区| 国内久久婷婷六月综合欲色啪| 亚洲一区高清亚洲精品| 免费在线观看影片大全网站| 亚洲av成人一区二区三| 久久久久精品国产欧美久久久| 亚洲欧洲精品一区二区精品久久久| 18美女黄网站色大片免费观看| 久久人妻福利社区极品人妻图片| 久久久久久久久久黄片| 亚洲中文字幕一区二区三区有码在线看 | 成人精品一区二区免费| 久久久久九九精品影院| 嫁个100分男人电影在线观看| av超薄肉色丝袜交足视频| 桃红色精品国产亚洲av| 丝袜人妻中文字幕| 制服丝袜大香蕉在线| 国产精品久久久久久人妻精品电影| 国产av又大| 十分钟在线观看高清视频www| 色av中文字幕| 麻豆一二三区av精品| 日日爽夜夜爽网站| 我的亚洲天堂| 免费在线观看视频国产中文字幕亚洲| 国产精品乱码一区二三区的特点| 亚洲精品在线美女| 日韩中文字幕欧美一区二区| 欧美日韩黄片免| 好男人电影高清在线观看| 午夜影院日韩av| 日韩三级视频一区二区三区| 久久香蕉激情| 成人精品一区二区免费| 国产男靠女视频免费网站| 国产亚洲欧美在线一区二区| 国内精品久久久久久久电影| 久久久久免费精品人妻一区二区 | 久久久久久久久久黄片| 国产色视频综合| 国内精品久久久久精免费| 免费在线观看影片大全网站| 国产亚洲欧美在线一区二区| 久久这里只有精品19| 久久久久久久久中文| 日韩精品中文字幕看吧| 美女免费视频网站| 久久精品成人免费网站| 级片在线观看| 久久香蕉精品热| 亚洲成av片中文字幕在线观看| 亚洲成av人片免费观看| 可以免费在线观看a视频的电影网站| 少妇熟女aⅴ在线视频| 在线观看一区二区三区| 在线看三级毛片| 999久久久精品免费观看国产| 国产亚洲精品久久久久5区| 午夜免费鲁丝| 少妇被粗大的猛进出69影院| 色综合欧美亚洲国产小说| 日韩精品青青久久久久久| 日韩高清综合在线| 亚洲色图av天堂| 天天添夜夜摸| 久久青草综合色| 一区二区日韩欧美中文字幕| 在线av久久热| 大型黄色视频在线免费观看| 两个人视频免费观看高清| 欧美日韩亚洲国产一区二区在线观看| 欧美亚洲日本最大视频资源| 欧美激情极品国产一区二区三区| 精品无人区乱码1区二区| 精品第一国产精品| 精品高清国产在线一区| 国内少妇人妻偷人精品xxx网站 | 日本熟妇午夜| 国产真人三级小视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 一卡2卡三卡四卡精品乱码亚洲| 亚洲一区中文字幕在线| 亚洲成人国产一区在线观看| 国产欧美日韩一区二区精品| 成人精品一区二区免费| 一a级毛片在线观看| 99riav亚洲国产免费| 法律面前人人平等表现在哪些方面| 国语自产精品视频在线第100页| 欧美黄色片欧美黄色片| 国产视频内射| 日本熟妇午夜| 日韩欧美三级三区| 女警被强在线播放| 搞女人的毛片| 日本精品一区二区三区蜜桃| 国产成人一区二区三区免费视频网站| 亚洲国产精品合色在线| 国产免费av片在线观看野外av| 好男人在线观看高清免费视频 | 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 很黄的视频免费| 亚洲av成人一区二区三| 99精品欧美一区二区三区四区| 99热只有精品国产| 国产亚洲欧美精品永久| 国产精品 欧美亚洲| 精品国产国语对白av| 不卡av一区二区三区| 成年版毛片免费区| 欧美+亚洲+日韩+国产| 丝袜美腿诱惑在线| www.精华液| 中文亚洲av片在线观看爽| 亚洲成人国产一区在线观看| 国产极品粉嫩免费观看在线| 不卡av一区二区三区| 国产免费av片在线观看野外av| 麻豆国产av国片精品| 国产熟女午夜一区二区三区| 国产单亲对白刺激| www日本在线高清视频| 国产精品综合久久久久久久免费| 99精品欧美一区二区三区四区| 99riav亚洲国产免费| 色av中文字幕| netflix在线观看网站| 成人精品一区二区免费| 国内少妇人妻偷人精品xxx网站 | 91在线观看av| 日韩免费av在线播放| 久久人人精品亚洲av| 欧美激情 高清一区二区三区| 国产亚洲欧美在线一区二区| 免费女性裸体啪啪无遮挡网站| 国产单亲对白刺激| 一区二区三区激情视频| 19禁男女啪啪无遮挡网站| av视频在线观看入口| 久99久视频精品免费| 午夜福利成人在线免费观看| 精品久久久久久久人妻蜜臀av| 99久久国产精品久久久| 亚洲免费av在线视频| 黄片播放在线免费| 国语自产精品视频在线第100页| 欧美三级亚洲精品| 国产精品久久久人人做人人爽| 亚洲精品国产精品久久久不卡| 久久人妻福利社区极品人妻图片| 国产成人av教育| 亚洲欧美日韩高清在线视频| 久久 成人 亚洲| 日韩欧美 国产精品| 免费在线观看影片大全网站| 国产av在哪里看| cao死你这个sao货| 中文字幕av电影在线播放| 欧美av亚洲av综合av国产av| 在线观看午夜福利视频| 成人永久免费在线观看视频| 天堂动漫精品| 久久精品国产99精品国产亚洲性色| 亚洲av成人不卡在线观看播放网| 日日干狠狠操夜夜爽| 精品久久久久久,| 丁香欧美五月| 中文字幕人妻丝袜一区二区| 99久久无色码亚洲精品果冻| 桃色一区二区三区在线观看| 欧美中文综合在线视频| 亚洲在线自拍视频| 精品久久久久久久久久免费视频| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 国产欧美日韩精品亚洲av| 亚洲人成伊人成综合网2020| 日韩国内少妇激情av| 成人国产一区最新在线观看| 中国美女看黄片| 日本一本二区三区精品| 国产主播在线观看一区二区| 搡老妇女老女人老熟妇| 丝袜美腿诱惑在线| 丝袜美腿诱惑在线| 国产av又大| 男女那种视频在线观看| 美女午夜性视频免费| 国产真实乱freesex| 无限看片的www在线观看| 韩国精品一区二区三区| 亚洲男人天堂网一区| 国产成人欧美在线观看| 欧美国产精品va在线观看不卡| 国产不卡一卡二| 在线天堂中文资源库| 丝袜在线中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品精品国产色婷婷| 国产成人精品无人区| 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 国产高清激情床上av| 99久久无色码亚洲精品果冻| 成人国语在线视频| 久久久久国产精品人妻aⅴ院| 国产免费av片在线观看野外av| 在线观看午夜福利视频| 久久午夜综合久久蜜桃| 亚洲人成网站在线播放欧美日韩| 国内揄拍国产精品人妻在线 | www.999成人在线观看| 久久久久九九精品影院| 亚洲人成网站高清观看| 亚洲熟女毛片儿| 亚洲av美国av| 国产极品粉嫩免费观看在线| 亚洲国产欧美网| 国产男靠女视频免费网站| 免费观看人在逋| 熟女少妇亚洲综合色aaa.| 国产亚洲av高清不卡| 十八禁网站免费在线| 午夜激情福利司机影院| 欧美国产精品va在线观看不卡| 日韩 欧美 亚洲 中文字幕| 丁香欧美五月| 国产精品自产拍在线观看55亚洲| 成人亚洲精品av一区二区| 亚洲 欧美 日韩 在线 免费| 女生性感内裤真人,穿戴方法视频| 婷婷丁香在线五月| 国产乱人伦免费视频| 99精品欧美一区二区三区四区| 十分钟在线观看高清视频www| 日韩国内少妇激情av| 亚洲成人国产一区在线观看| 夜夜夜夜夜久久久久| 黄色成人免费大全| videosex国产| 色av中文字幕| 黄色 视频免费看| 黄色视频,在线免费观看| 级片在线观看| 国产精品影院久久| 亚洲av熟女| 国产成人av激情在线播放| 色在线成人网| 91国产中文字幕| 欧美绝顶高潮抽搐喷水| 国产成年人精品一区二区| 欧美黄色淫秽网站| 香蕉丝袜av| 国产av一区在线观看免费| 国产亚洲av高清不卡| 亚洲性夜色夜夜综合| 精品久久蜜臀av无| 999久久久精品免费观看国产| 男女视频在线观看网站免费 | 日韩欧美在线二视频| 亚洲精品国产区一区二| 狠狠狠狠99中文字幕| 真人做人爱边吃奶动态| 日韩欧美一区二区三区在线观看| 黄色片一级片一级黄色片| 亚洲专区国产一区二区| 99在线人妻在线中文字幕| 999精品在线视频| 亚洲精品国产精品久久久不卡| 精品午夜福利视频在线观看一区| 国产精品综合久久久久久久免费| 18禁美女被吸乳视频| 丁香六月欧美| 丰满的人妻完整版| 美女 人体艺术 gogo| 在线观看免费视频日本深夜| 美女高潮到喷水免费观看| 国产人伦9x9x在线观看| 欧美日本亚洲视频在线播放| 中亚洲国语对白在线视频| 中文字幕高清在线视频| 在线永久观看黄色视频| 中文字幕精品亚洲无线码一区 | 国产真人三级小视频在线观看| 国产真实乱freesex| 免费无遮挡裸体视频| 亚洲精品国产精品久久久不卡| 色尼玛亚洲综合影院| 91成人精品电影| 久久香蕉激情| 变态另类成人亚洲欧美熟女| 免费高清在线观看日韩| 高潮久久久久久久久久久不卡| 国内揄拍国产精品人妻在线 | 亚洲欧美日韩高清在线视频| 美女午夜性视频免费| 别揉我奶头~嗯~啊~动态视频| 国产精品免费视频内射| 搡老妇女老女人老熟妇| 啦啦啦韩国在线观看视频| 国产不卡一卡二| 国产日本99.免费观看| 久久久久久免费高清国产稀缺| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人| 亚洲精品久久成人aⅴ小说| 午夜激情福利司机影院| 999久久久精品免费观看国产| 男女之事视频高清在线观看| 精品久久久久久久久久免费视频| 丁香欧美五月| 女警被强在线播放| 激情在线观看视频在线高清| 久久久久久大精品| 精品第一国产精品| 无遮挡黄片免费观看| tocl精华| 少妇 在线观看| 国产精品精品国产色婷婷| 久久久久久久久免费视频了| 国产成人系列免费观看| 欧美乱妇无乱码| 757午夜福利合集在线观看| 国产黄a三级三级三级人| 性欧美人与动物交配| 亚洲精品在线美女| 亚洲精品国产精品久久久不卡| 色尼玛亚洲综合影院| 久久久国产欧美日韩av| 欧美人与性动交α欧美精品济南到| 国产黄片美女视频| av在线天堂中文字幕| 久久热在线av| 国产在线观看jvid| 国产极品粉嫩免费观看在线| 午夜两性在线视频| 两个人视频免费观看高清| 韩国av一区二区三区四区| 国产精品一区二区三区四区久久 | 国产激情欧美一区二区| 日本 av在线| 久久精品国产99精品国产亚洲性色| 大型黄色视频在线免费观看| 麻豆av在线久日| 欧美精品啪啪一区二区三区| 亚洲中文av在线| 精品第一国产精品| 一区二区日韩欧美中文字幕| 国产亚洲欧美98| 成人av一区二区三区在线看| 久久性视频一级片| 麻豆成人av在线观看| 亚洲av熟女| 欧美乱码精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 日韩一卡2卡3卡4卡2021年| 午夜福利高清视频| 午夜精品在线福利| 国语自产精品视频在线第100页| 久久中文字幕一级| xxx96com| 国产不卡一卡二| 夜夜夜夜夜久久久久| 久久久久久九九精品二区国产 | 老鸭窝网址在线观看| 亚洲专区字幕在线| 久久久久久国产a免费观看| 观看免费一级毛片| 免费av毛片视频| 成人18禁在线播放| 国产精品香港三级国产av潘金莲| 99久久精品国产亚洲精品| 男人舔女人的私密视频| 久久午夜亚洲精品久久| 亚洲精品在线美女| 一本久久中文字幕| 久久香蕉精品热| 中亚洲国语对白在线视频| 国产黄片美女视频| 大型av网站在线播放| 天天一区二区日本电影三级| 丝袜美腿诱惑在线| 色播在线永久视频| 91成年电影在线观看| 国产激情欧美一区二区| 黄色视频不卡| 久久久久久久久免费视频了| 18禁国产床啪视频网站| svipshipincom国产片| 久久亚洲精品不卡| 亚洲国产中文字幕在线视频| 亚洲精品久久成人aⅴ小说| 人妻久久中文字幕网| 午夜影院日韩av| 欧美黄色淫秽网站| 国产三级在线视频| 十分钟在线观看高清视频www| 男女午夜视频在线观看| 国产成人一区二区三区免费视频网站| 97碰自拍视频| 人人妻人人澡人人看| 国产爱豆传媒在线观看 | 大型av网站在线播放| 欧美绝顶高潮抽搐喷水| 18禁美女被吸乳视频| 色播在线永久视频| 国产三级黄色录像| 国产1区2区3区精品|