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

    輻射對(duì)高原渦形成和發(fā)展影響的模擬研究

    2022-06-01 07:15:34鄧中仁葛旭陽姚秀萍陳明誠
    大氣科學(xué) 2022年3期
    關(guān)鍵詞:渦的渦度低層

    鄧中仁 葛旭陽 ,2 姚秀萍 陳明誠

    1 南京信息工程大學(xué)氣象災(zāi)害教育部重點(diǎn)實(shí)驗(yàn)室/氣候與環(huán)境變化國際合作聯(lián)合實(shí)驗(yàn)室, 南京 210044

    2 安陽國家氣候觀象臺(tái), 河南安陽 455000

    3 中國氣象局氣象干部培訓(xùn)學(xué)院, 北京 100081

    1 引言

    青藏高原對(duì)整個(gè)東亞地區(qū)的大氣環(huán)流具有重要的動(dòng)力及熱力作用(Ye, 1950; 葉篤正等, 1957)。高原低渦作為北半球夏季青藏高原上特殊的低壓系統(tǒng),其水平尺度為400~800 km,垂直尺度可達(dá)2~3 km(葉篤正和高由禧, 1979),是直接影響高原地區(qū)降水的主要系統(tǒng)。在一定的環(huán)流背景下,部分高原渦可以東移出高原地區(qū),從而導(dǎo)致我國長江中下游地區(qū)的強(qiáng)對(duì)流天氣過程(楊克明等, 2001; 何光碧等, 2010; 郁淑華和高文良, 2010),因此研究高原渦發(fā)生及發(fā)展的機(jī)制具有重要意義。

    國內(nèi)外已有大量研究針對(duì)高原渦發(fā)生、發(fā)展機(jī)制進(jìn)行探討(Dell’ Osso and Chen, 1986; 李國平和劉紅武, 2006; Zheng et al., 2013)。一些研究認(rèn)為在高原渦生成階段,凝結(jié)潛熱起到了至關(guān)重要的作用(Wang, 1987; Li et al., 2011);另一些研究則指出地表感熱對(duì)于高原渦的生成更為關(guān)鍵,而凝結(jié)潛熱對(duì)于高原渦的發(fā)展及維持具有顯著的影響(Shen et al., 1986a; 李黎等, 2019)。Li et al.(2014, 2018)研究表明,在白天高原渦發(fā)生頻率較低,而夜間其發(fā)生頻率迅速增加。這意味著高原渦的發(fā)生頻率具有明顯的日變化特征。輻射過程調(diào)制對(duì)流活動(dòng)日變化的物理機(jī)制相對(duì)清楚。Webster and Stephens(1980)指出,夜間云頂存在長波輻射冷卻,導(dǎo)致溫度遞減率加大,從而促使了對(duì)流的發(fā)展;而在白天,云頂受到短波輻射加熱的影響,溫度遞減率減小,靜力穩(wěn)定度增加抑制對(duì)流過程。作為熱帶地區(qū)典型的災(zāi)害性天氣系統(tǒng),熱帶氣旋也具有明顯的日變化特征(Kossin, 2002; Dunion et al., 2014)。大量研究表明輻射日變化對(duì)熱帶氣旋的生成有明顯的影響(Melhauser and Zhang, 2014; Cao et al., 2014;Tang and Zhang, 2016), Melhauser and Zhang(2014)指出,夜間由于高層輻射冷卻大于低層,使得靜力穩(wěn)定度減小,從而促進(jìn)了低渦的發(fā)展。Tang and Zhang(2016)指出,太陽短波輻射有利于對(duì)流層高層增溫,加強(qiáng)大氣靜力穩(wěn)定度,從而抑制對(duì)流活動(dòng)發(fā)展,由此可見,輻射日變化對(duì)熱帶氣旋的生成有明顯的影響。由于高原渦與熱帶氣旋具有一定的相似性,如:兩者都具有眼的結(jié)構(gòu),并且均為暖心結(jié)構(gòu),但因?yàn)闊釒庑膹?qiáng)度、所在區(qū)域與高原渦不一樣,所以,兩者暖心的強(qiáng)度和所處高度也會(huì)有所不同(羅四維和楊洋, 1992; 羅四維等,1993)。因此,本文借鑒前人研究熱帶氣旋的方法(Cao et al., 2014; Tang and Zhang, 2016),通過個(gè)例模擬來探究輻射日變化對(duì)高原渦發(fā)展的影響機(jī)制。

    2 數(shù)值模式

    2.1 模式及數(shù)值試驗(yàn)設(shè)計(jì)

    本文采用中尺度數(shù)值模式WRF-ARW(3.9 版本),模 擬 時(shí) 間 為2013 年6 月3 日00:00 至7 日00:00(協(xié)調(diào)世界時(shí),下同)。模式的初始條件及邊界條件采用逐6 小時(shí),水平分辨率為0.5°×0.5°的ERA-Interim 再分析資料。模式區(qū)域采用雙層嵌套,中心為(32.5°N,102.5°E),網(wǎng)格點(diǎn)分別為295×176 和298×151,水平分辨率為27 km 和9 km,模式垂直層次為55 層。采用的物理參數(shù)化方案:Eta(Ferrier)云微物理方案、RRTM 長波輻射方案、Dudhia 短波輻射方案、MM5 Monin-Obukhov近地面層方案、Noah 陸面過程方案以及Betts-Miller-Janjic 積云參數(shù)化方案。

    在Dudhia 短波輻射方案(Dudhia, 1989)中,太陽短波輻射的強(qiáng)度隨太陽與地球表面之間的夾角而變化。例如在中午時(shí)分天頂角最大,即對(duì)應(yīng)的短波輻射強(qiáng)度最大。反之,在午夜短波輻射強(qiáng)度則為零。因此通過改變模式中的當(dāng)?shù)貢r(shí)間,可以調(diào)整太陽高度角,進(jìn)而改變太陽短波輻射強(qiáng)度。本文設(shè)計(jì)以下三組試驗(yàn)(表1):控制試驗(yàn)(CTL)具有完整的短波輻射日循環(huán);在All_day 試驗(yàn)中,將模式局地時(shí)間固定在正午12:00,以此得到最強(qiáng)的短波輻射強(qiáng)度;在All_night 試驗(yàn)中,將整個(gè)模式積分過程中局地時(shí)間設(shè)置為午夜00:00,以此去掉太陽短波輻射的影響。

    表1 試驗(yàn)的描述Table 1 Description of the experiments

    2.2 高原渦個(gè)例

    圖1 給出了本次高原渦的觀測路徑以及強(qiáng)度,其中強(qiáng)度以500 hPa 渦旋中心最低位勢高度表示。此次高原低渦于6 月4 日12:00 生成在高原東部曲麻萊地區(qū),5 日06:00 東移出高原,此后一路向東,于8 日06:00 入海,并逐漸減弱消亡。本文重點(diǎn)討論,高原渦在高原上的時(shí)段。為檢驗(yàn)控制試驗(yàn)的模擬效果,圖2 給出了2013 年6 月4 日06:00、12:00基于ERA-Interim 再分析資料給出的此次高原渦個(gè)例及CTL 試驗(yàn)?zāi)M的高原渦情況。由圖2a、b 看出,在4 日06:00,低層500 hPa 有豎切變線的存在。CTL 試驗(yàn)中高原渦生成前的環(huán)流場和高原渦生成的位置與再分析資料中環(huán)流場形勢較為一致,說明CTL 試驗(yàn)較好的再現(xiàn)了高原渦生成過程。因此,本文將利用CTL 試驗(yàn)結(jié)果對(duì)該高原渦演變過程進(jìn)行診斷分析。

    圖1 2013 年6 月4 日12:00 至10 日00:00 高原渦觀測路徑。紅色數(shù)字表示時(shí)間,括號(hào)內(nèi)數(shù)值為500 hPa 最低位勢高度(單位: dagpm)Fig. 1 Observation track of the TPV (Tibetan Plateau vortex) from 1200 UTC 4 to 0000 UTC 10 June 2013. Red numbers indicate the time (UTC),the lowest geopotential height (units: dagpm) at 500 hPa is indicated in the brackets

    圖2 2013 年6 月4 日(a、c)06:00、(b、d)12:00(a、b)基于再分析資料以及(c、d)CTL 試驗(yàn)?zāi)M的500 hPa 等壓面上的流場(矢量,單位:m s-1)分布,等值線為地形高度在3000 m 以上的區(qū)域,三角形為高原渦所處的位置Fig. 2 Distributions of wind field (vectors, units: m s-1) from (a, b) reanalysis data, (c, d) CTL (control run) experiment at 500 hPa at 0600 UTC (left column), 1200 UTC (right column) June 4, 2013. Contours outline the Tibetan Plateau with an altitude higher than 3000 m; the triangle shows the location of the TPV

    3 模擬結(jié)果分析

    3.1 高原渦的演變特征

    3.1.1 高原渦發(fā)展情況

    前人研究指出,小尺度對(duì)流單體對(duì)熱帶氣旋的生成、加強(qiáng)起著重要的作用(Hendricks et al., 2004;Montgomery et al., 2006; 郭春蕊和張慶紅, 2012)??紤]到高原渦與熱帶氣旋之間存在諸多的相似性,本文首先探討小尺度對(duì)流單體的發(fā)展、組織情況,這對(duì)了解高原渦生成就有著十分重要的意義。本文利用小尺度系統(tǒng)的垂直渦度來表征強(qiáng)對(duì)流單體的活動(dòng)情況,需要說明的是:為了區(qū)分不同尺度系統(tǒng)之間的相互作用,首先進(jìn)行了空間濾波處理(Cao et al., 2014; 徐夢婷等, 2016)。在本文中將波長大于100 km 以上的系統(tǒng)當(dāng)作背景環(huán)境場;相應(yīng)的,波長在100 km 以下的分量則作為小尺度系統(tǒng)。通過改變不同濾波波長閾值(如100~300 km),所得結(jié)果基本接近。圖3 揭示了2013 年6 月4 日09:00、15:00 和21:00,三個(gè)時(shí)次500 hPa 大尺度背景環(huán)流和中小尺度相對(duì)渦度隨時(shí)間演變特征,其中綠色方框?yàn)橄挛挠?jì)算平均值所用到的渦旋中心區(qū)域(300 km×300 km)。

    從動(dòng)力學(xué)角度來看,根據(jù)“渦度隔離”效應(yīng)(渦旋β 效應(yīng)),即當(dāng)環(huán)境場存在渦度梯度時(shí),小尺度正(負(fù))渦度擾動(dòng)趨向(背向)于環(huán)境場最大渦度中心,而隔離速率與環(huán)境場渦度梯度成正比(徐夢婷等, 2016)。從圖3 中可以看到4 日09:00,CTL 試驗(yàn)中前期環(huán)境場相對(duì)渦度較小,小尺度對(duì)流單體分布散亂,沒有明顯的向環(huán)境場渦度中心聚集特征。對(duì)比圖3d 不難發(fā)現(xiàn),All_night 試驗(yàn)中的環(huán)境場渦度較CTL 試驗(yàn)要大,且存在明顯的徑向梯度,即該組試驗(yàn)中“渦度隔離”效應(yīng)使得小尺度對(duì)流渦旋更早的向大尺度低渦中心聚集,為高原渦的形成提供了更好的初始條件,這與熱帶氣旋中積云熱塔的作用類似(Montgomery et al., 2006; Ge et al., 2015)。在4 日15:00,CTL 試驗(yàn)開始呈現(xiàn)出較為明顯的徑向梯度,小尺度渦旋也隨之開始向低渦中心聚集。此時(shí)All_night 試驗(yàn)中,隨著小尺度對(duì)流系統(tǒng)向內(nèi)聚集,有機(jī)會(huì)聚集成一定組織化結(jié)構(gòu),這使得All_night 試驗(yàn)中高原渦強(qiáng)度發(fā)展的更為強(qiáng)盛(圖3e)。到了4 日21:00,CTL 試驗(yàn)中小尺度對(duì)流系統(tǒng)在低渦中心附近形成了一個(gè)組織化的系統(tǒng),促使高原渦進(jìn)一步的加強(qiáng)。相較之下,在All_day試驗(yàn)中,環(huán)境場相對(duì)渦度較同期的CTL 和All_night 試驗(yàn)都要小,導(dǎo)致“渦度隔離”效應(yīng)不顯著,致使小尺度對(duì)流系統(tǒng)分布散亂,沒有明顯的向低渦中心聚集,主要集中在渦區(qū)外圍,這與高原渦發(fā)展緩慢相符。值得關(guān)注的是,到4 日21:00,All_night試驗(yàn)中環(huán)境場相對(duì)渦度減弱,小尺度對(duì)流系統(tǒng)的組織化衰減,高原渦較早地進(jìn)入消亡階段。關(guān)于All_night 試驗(yàn)中高原渦快速衰亡的原因?qū)⒃谙挛脑敿?xì)討論。

    圖3 2013 年6 月4 日09:00(左列)、15:00(中間列)和21:00(右列)模擬的500 hPa 等壓面上的環(huán)境場大尺度相對(duì)渦度(陰影,單位:10-5 s-1)、小尺度系統(tǒng)相對(duì)渦度(等值線,單位:10-5 s-1)和風(fēng)場(矢量,單位:m s-1)分布:(a-c)CTL 試驗(yàn);(d-f)All_night 試驗(yàn);(g-i)All_day 試驗(yàn)。綠色方框?yàn)闇u旋中心區(qū)域(300 km×300 km)Fig. 3 Distributions of large-scale relative vorticity (shaded, units: 10-5 s-1), small-scale system vorticity (contours, units: 10-5 s-1) and wind field(vectors, units: m s-1) from (a-c) CTL, (d-f) All_night experiment, and (g-i) All_day experiment at 500 hPa at 0900 UTC (left column), 1500 UTC(middle column), and 2100 UTC (right column) 4 June 2013. The green box indicates the central area of the TPV(300 km×300 km)

    上述分析表明,高原渦發(fā)展快慢與其所處區(qū)域?qū)α骰顒?dòng)有關(guān)。為進(jìn)一步比較三組試驗(yàn)中高原渦內(nèi)區(qū)對(duì)流的活動(dòng)情況,圖4 統(tǒng)計(jì)了4 日12:00~18:00時(shí)間段內(nèi),低渦內(nèi)區(qū)不同高度下具有特定大小的垂直速度的格點(diǎn)數(shù)占區(qū)域內(nèi)總格點(diǎn)數(shù)的比例CFAD(Contoured Frequency by Altitude Diagram;Yuter and Houze, 1995)。統(tǒng)計(jì)范圍取低渦中心附近300 km×300 km 區(qū)域。這里橫坐標(biāo)表征垂直速度的大小,陰影為特定垂直速度所占的比例。由圖4 可見,三組試驗(yàn)中較小的垂直速度(-0.2 m s-1<w<0.2 m s-1)占據(jù)了大多數(shù)網(wǎng)格點(diǎn),大于0.4 m s-1的強(qiáng)上升運(yùn)動(dòng)所占據(jù)的格點(diǎn)數(shù)很小,此現(xiàn)象表明在低渦附近僅有一小部分強(qiáng)對(duì)流活動(dòng)。這與前人研究熱帶氣旋的情況一致(Steranka et al., 1986; Hendricks et al.,2004)。而對(duì)比CTL 和All_night 試驗(yàn),可以清楚地看到,由于All_night 試驗(yàn)中高原渦發(fā)展較快,因此,在這一時(shí)間段內(nèi)的All_night 試驗(yàn)中,大于0.4 m s-1的垂直速度所占比例更大,這意味著對(duì)流活動(dòng)更強(qiáng)。而All_day 試驗(yàn)中,該比例則最小,說明對(duì)流活動(dòng)較弱,低渦發(fā)展最慢。

    圖4 2013 年6 月4 日高原渦中心附近(300 km×300 km)12:00~18:00 時(shí)間平均的CFAD(Contoured Frequency by Altitude Diagram,陰影)分布:(a)CTL 試驗(yàn);(b)All_night 試驗(yàn);(c)All_day 試驗(yàn)Fig. 4 Averaged CFAD ( Contoured Frequency by Altitude Diagram; shaded) during 1200 UTC-1800 UTC 4 June 2013, calculated over a box area with a radius of 300 km around the vortex center: (a) CTL; (b) All_night experiment; (c) All_day experiment

    為比較高原渦的發(fā)展過程中結(jié)構(gòu)差異,圖5 給出了三組試驗(yàn)中低渦內(nèi)區(qū)相對(duì)渦度隨高度變化的演變過程。在CTL 試驗(yàn)中,高原渦在白天發(fā)展較為緩慢,最大渦度中心出現(xiàn)在450~400 hPa 之間,而隨著時(shí)間的推移,到了夜間4 日15:00 以后,低渦強(qiáng)度逐漸加強(qiáng),同時(shí)最大渦度中心開始向低層延伸。在All_day 試驗(yàn)中,正渦度強(qiáng)度較小,且最大渦度中心到4 日18:00 后才開始緩慢向低層延伸,因此高原渦未得到明顯的發(fā)展。相較于CTL 試驗(yàn),All_night 試驗(yàn)中低渦發(fā)展速度的更快,最大渦度中心更早的向低層延伸。4 日06:00 低渦迅速加強(qiáng),到了4 日12:00~15:00 低渦達(dá)到成熟階段而后逐漸減弱。

    圖5 2013 年6 月4 日00:00~23:00 低渦中心附近區(qū)域平均(300 km×300 km)相對(duì)渦度(陰影,單位: 10-5 s-1)的垂直分布:(a)CTL 試驗(yàn);(b)All_night 試驗(yàn);(c)All_day 試驗(yàn)Fig. 5 Time evolution of the vertical profile of vorticity (shaded, units: 10-5 s-1) during 0000 UTC-2300 UTC 4 June 2013, averaged over a box area with a radius of 300 km around the vortex center: (a) CTL experiment; (b) All_night experiment; (c) All_day experiment

    綜上所述,高原渦正渦度中心開始出現(xiàn)在450~400 hPa(中層)之間,隨著低渦的發(fā)展,正渦度中心開始向下延伸,這類似于熱帶氣旋生成中Top-down 理 論(Bister and Emanuel,1997)。隨著中尺度渦旋的發(fā)展并向下延伸時(shí),有利于低層氣旋性環(huán)流發(fā)展。

    Zhang and Fritsch(1987)指出,在中尺度對(duì)流系統(tǒng)中,層狀云往往對(duì)應(yīng)著低層冷卻下沉,中層凝結(jié)加熱上升。根據(jù)連續(xù)方程,即存在中層輻合,低層輻散;而對(duì)流云呈現(xiàn)為整層均為凝結(jié)加熱特征,動(dòng)力場上對(duì)應(yīng)著低層(高層)輻合(輻散)。圖6、圖7 給出了4 日00:00~23:00 三組試驗(yàn)中凝結(jié)潛熱及散度垂直分布特征。由圖6 可見,4 日00:00~12:00,CTL 試驗(yàn)中呈現(xiàn)出低層冷卻、中層加熱的情況;與此同時(shí),最大輻合中心位于中層,這與層狀云中加熱、散度廓線相似。到了4 日15:00 以后,中層加熱強(qiáng)度增大,并且低層冷卻逐漸轉(zhuǎn)為凝結(jié)加熱;同時(shí)最強(qiáng)輻合中心開始出現(xiàn)在中低層,即表現(xiàn)為低層輻合、高層輻散的情形。這說明高原渦經(jīng)歷了從層狀云向強(qiáng)對(duì)流云特征轉(zhuǎn)化。All_day 試驗(yàn)中4 日00:00~23:00,低層都有凝結(jié)冷卻,而中層凝結(jié)加熱強(qiáng)度小,從散度場分布來看低層輻合始終較CTL 和All_night 試驗(yàn)弱,表現(xiàn)為層狀云特征,這也意味著強(qiáng)對(duì)流發(fā)展緩慢。在4 日06:00~15:00期間,All_night 試驗(yàn)中,低層和中層均為凝結(jié)加熱,最大輻合中心位于低層并且與中層輻散中心配合,對(duì)應(yīng)著強(qiáng)對(duì)流的發(fā)展;到了4 日15:00 以后,凝結(jié)加熱作用減弱,低層輻合強(qiáng)度減小,高原渦進(jìn)入衰減階段。

    圖6 同圖5,但為凝結(jié)潛熱(陰影,單位:K h-1)Fig. 6 Same as Fig. 5, but for latent heating (shaded, units: K h-1)

    圖7 同圖5,但為散度垂直分布(陰影,單位:10-5 s-1)Fig. 7 Same as Fig. 5, but for the divergence (shaded, units: 10-5 s-1)

    3.1.2 Okubo-Weiss(OW)指數(shù)

    以上分析表明,對(duì)應(yīng)不同的輻射過程,三組試驗(yàn)中高原渦發(fā)展存在明顯的差異。為進(jìn)一步揭示差異,圖8 給出了三組試驗(yàn)中伴隨高原渦發(fā)展過程中Okubo-Weiss(OW)指數(shù)的演變特征。這里OW=ζ2-Sn2-S2s, 其中,相對(duì)渦度ζ=?v/?x-?u/?y,Ss=?u/?x-?v/?y,Sn=?v/?x+?u/?y。前人研究(Isern-Fontanet et al., 2004; 管靚等, 2016)指出,OW 指數(shù)的大小可以反映出氣旋的強(qiáng)度。當(dāng)OW>0 時(shí),即旋轉(zhuǎn)效應(yīng)大于形變效應(yīng),有利于對(duì)流的發(fā)展。反之,負(fù)的OW 區(qū)域表明形變效應(yīng)占主導(dǎo)作用,則不利于對(duì)流發(fā)展及維持。由圖8 可以看到,在4日09:00,All_night 試驗(yàn)中OW 指數(shù)的形態(tài)、位置與渦度中心的形態(tài)和位置比較吻合,且數(shù)值比CTL 和All_day 試驗(yàn)要大,表明該時(shí)刻系統(tǒng)旋轉(zhuǎn)性較強(qiáng),有利于氣旋性環(huán)流的產(chǎn)生。類似地,4 日15:00,CTL 試驗(yàn)中,OW 指數(shù)大值中心逐漸與相對(duì)渦度中心重合,旋轉(zhuǎn)性加強(qiáng)意味著高原渦的發(fā)展。相較之下,OW 指數(shù)數(shù)值較同時(shí)刻All_night 試驗(yàn)要小很多,這與高原渦強(qiáng)度比All_night 試驗(yàn)中要弱的情況一致。而在All_day 試驗(yàn)中,OW 指數(shù)的值最小,甚至出現(xiàn)大范圍的負(fù)值區(qū),即形變效應(yīng)大于旋轉(zhuǎn)效應(yīng)。該時(shí)刻低渦四周風(fēng)場以變形作用為主,有利于切變線的生成和維持,不利于氣旋性環(huán)流的形成。在4 日21 時(shí),All_night 試驗(yàn)中OW 指數(shù)數(shù)值較4 日15:00 有明顯的減小,即高原渦逐漸衰減。此后,CTL 試驗(yàn)中OW 指數(shù)中心數(shù)值明顯大于All_night 試驗(yàn)。而All_day 試驗(yàn)中OW 指數(shù)大值中心的范圍與強(qiáng)度都明顯偏弱,表明高原渦發(fā)展速度較慢。這些演變特征與圖3 中相對(duì)渦度演變過程一致。

    圖8 2013 年6 月4 日09:00(左列)、15:00(中間列)和21:00(右列)500 hPa 等壓面上的環(huán)境場大尺度相對(duì)渦度(陰影,單位:10-5 s-1)、OW 指數(shù)(等值線,單位:10-9 s-2)分布:(a)CTL 試驗(yàn);(b)All_night 試驗(yàn);(c)All_day 試驗(yàn)Fig. 8 Distributions of environmental vorticity (shaded, units: 10-5 s-1), OW index (contour, units: 10-9 s-2) at 500 hPa at 0900 UTC (left column),1500 UTC (middle column) and 2100 UTC (right column) 4 June 2013: (a) CTL experiment; (b) All_night experiment; (c) All_day experiment

    3.2 動(dòng)力條件

    3.2.1 非平衡項(xiàng)分析

    上述模擬結(jié)果比較清楚的表明了三組試驗(yàn)的差異,接下來本文將針對(duì)熱、動(dòng)力學(xué)過程展開討論。對(duì)于中尺度高原渦系統(tǒng)而言,散度場有著至關(guān)重要的影響(馬林等, 2003; 李山山, 2017)。為此,本文利用散度方程來討論三組試驗(yàn)中高原渦生成過程的差異。

    圖9 2013 年6 月4 日09:00(左列)、15:00(中間列)和21:00(右列)500 hPa 等壓面上非平衡項(xiàng)(陰影,單位:10-9 s-2)、風(fēng)場(矢量,單位:m s-1)分布:(a)CTL 試驗(yàn);(b)All_night 試驗(yàn);(c)All_day 試驗(yàn)Fig. 9 Distributions of non-balance terms (shaded, units: 10-9 s-2 ), wind field (vectors, units: m s-1) at 500 hPa at 0900 UTC (left column), 1500 UTC(middle column) and 2100 UTC (right column) 4 June 2013: (a) CTL experiment; (b) All_night experiment; (c) All_day experiment

    上述分析表明,在高原渦成熟階段,其中心存在明顯的非平衡項(xiàng)正值中心,而負(fù)值中心位于低渦外側(cè)。這可能源于高原渦具有熱帶氣旋類似結(jié)構(gòu)導(dǎo)致的結(jié)果。為驗(yàn)證這種可能性,本文繪制了三組試驗(yàn)中溫度距平、垂直運(yùn)動(dòng)和切向風(fēng)的軸對(duì)稱分量的半徑—高度的垂直剖面圖(圖10 和圖11)。從圖10、圖11 可以看到,當(dāng)?shù)蜏u中心對(duì)應(yīng)非平衡項(xiàng)的正值中心時(shí),低渦中心存在下沉運(yùn)動(dòng),同時(shí)低渦中心外圍有較強(qiáng)的上升運(yùn)動(dòng)。一般而言,強(qiáng)的下沉運(yùn)動(dòng)利于絕熱增溫。對(duì)比圖10 可以看出,由于CTL 和All_night 試驗(yàn)中高原渦中心都對(duì)應(yīng)有下沉運(yùn)動(dòng),而下沉絕熱增溫使得兩組試驗(yàn)中暖心強(qiáng)度均比All_day 試驗(yàn)要強(qiáng),這也表明高原渦中心的下沉運(yùn)動(dòng)對(duì)其暖心結(jié)構(gòu)的形成起到了一定的作用(圖11)。CTL 試驗(yàn)中,切向風(fēng)最大風(fēng)速半徑為70 km,且最大風(fēng)速達(dá)到8 m s-1。同時(shí)低層為氣旋性環(huán)流,到了200 hPa 轉(zhuǎn)為反氣旋,即切向風(fēng)隨高度遞減。為滿足熱成風(fēng)平衡關(guān)系,應(yīng)具有暖心結(jié)構(gòu)。以上分析表明高原渦具有與熱帶氣旋相近的動(dòng)力和熱力特征。此外,在All_night 試驗(yàn),最大風(fēng)速半徑更?。?5 km)、最大風(fēng)速為7 m s-1。而在All_day 試驗(yàn)中,最大風(fēng)速半徑則明顯更大(180 km)、最大風(fēng)速為4 m s-1。這也清楚的反映出三者強(qiáng)度差異。

    圖10 2013 年6 月4 日21:00(a)CTL、(b)All_night 和(c)All_day 試驗(yàn)中垂直運(yùn)動(dòng)(陰影,單位:10-2 m s-1)和切向風(fēng)速(等值線,單位:m s-1)的軸對(duì)稱分量的半徑—高度的垂直剖面Fig. 10 Radius-height cross sections of the axisymmetric component of vertical velocity (shaded, units: 10-2 m s-1) and tangential wind speed(contour, units: m s-1) in (a) CTL, (b) All_night, and (c) All_day experiments at 2100 UTC 4 June 2013

    圖11 同圖10,但為溫度距平(單位:K)Fig. 11 Same as Fig 10, but for the temperature anomaly (units: K)

    3.2.2 正壓動(dòng)能轉(zhuǎn)換

    高原渦形成過程涉及多尺度環(huán)流的相互作用。為了探究不同尺度分量對(duì)高原渦生成的動(dòng)力學(xué)影響,本文采用正壓動(dòng)能方程(BKE)對(duì)擾動(dòng)動(dòng)能進(jìn)行診斷分析。

    圖12 2013年6月4日(a-c)CTL、(d-f)All_night和(g-i)All_day試驗(yàn)的500hPa 等壓面上公式(2)中擾動(dòng)動(dòng)能傾向(左列)、緯向風(fēng)輻合-(中間列)和徑向風(fēng)輻合(右列)00:00~12:00時(shí)間平均結(jié)果分布(陰影,單位:10-5 m2s-3),三角形為高原渦所處的位置Fig.12Horizontalpattern of kineticenergy tendenc′(leftcolumn), latitudinal wind convergence - (middle column)andradialwind convergence?(right column )inEq.2 on the 500 hPa (shaded,units:10-5m2 s-3),the triangle reflects thelocation of theTPVin(a-c)CTL,(df) All_night, and (g-i) All_day. Time averaged from 0000 UTC-1200 UTC 4 June 2013

    3.3 熱力學(xué)條件

    3.3.1 水汽條件

    非絕熱加熱過程是影響高原渦生成的重要因素(Wang et al., 1987; 羅四維和楊洋, 1992; 羅四維等,1993; 宋雯雯等, 2012; 田珊儒等, 2015)。因此,本文比較了三組試驗(yàn)中相對(duì)濕度場的演變特征(圖13)。All_night 試驗(yàn)中去掉了短波輻射,在長波輻射冷卻的作用下環(huán)境場溫度降低,使得相對(duì)濕度增加,有利于降水的產(chǎn)生。鑒于降水效率是表征降水的一個(gè)重要指標(biāo)(Li and Gao, 2012; Zhou et al., 2014),為此本文計(jì)算了降水效率(PE)=(降水量/大氣可降水量)。對(duì)比圖13 和圖14 可知,前期CTL 試驗(yàn)中在短波輻射的影響下,白天大氣溫度升高導(dǎo)致相對(duì)濕度較低(圖13),從而不利于降水的產(chǎn)生。而到了夜間4 日15:00 以后,隨著相對(duì)濕度的增加,降水效率迅速增強(qiáng),為高原渦的發(fā)展提供了大量的潛熱能量。而在All_night 試驗(yàn)中,強(qiáng)降水出現(xiàn)最早(主要發(fā)生在4 日03:00~15:00),且降水效率最高達(dá)到0.38,大量凝結(jié)潛熱釋放為低渦的快速發(fā)展提供了條件。在4 日15:00 以后降水強(qiáng)度快速減弱,高原渦趨于衰亡,這可能與水汽供應(yīng)有關(guān)。在高原這樣特殊地形情況下,水汽大量消耗的同時(shí)往往得不到充分補(bǔ)充,因此水汽條件快速減弱進(jìn)而導(dǎo)致對(duì)流消亡。在整個(gè)All_day 試驗(yàn)中,相對(duì)濕度及降水強(qiáng)度都保持較低水平,潛熱釋放最弱,高原渦發(fā)展緩慢。由此可見,水汽條件在一定程度上可以決定高原渦的生命史。從觀測事實(shí)中可以發(fā)現(xiàn),大多數(shù)高原渦的發(fā)生、發(fā)展和消亡都位于高原上,較少的高原渦會(huì)東移出高原,這也可能與水汽的輸送和補(bǔ)充等條件有關(guān),值得進(jìn)一步討論。

    圖13 2013 年6 月4 日00:00~23:00(a)CTL、(b)All_night 和(c)All_day 試驗(yàn)中低渦中心附近區(qū)域平均(300 km×300 km)相對(duì)濕度(陰影)垂直分布隨時(shí)間的變化Fig. 13 Time evolution of the vertical profile of relative humidity (shaded), averaged over a box area centered at the TPV during 0000 UTC-2300 UTC 4 June 2013: (a) CTL experiment; (b) All_night experiment; (c) All_day experiment

    圖14 2013 年6 月4 日00:00~23:00(a)CTL、(b)All_night 和(c)All_day 試驗(yàn)中低渦中心附近區(qū)域平均(300 km×300 km)的降水效率隨時(shí)間的變化Fig. 14 Time evolution of PE (Precipitation Efficiency) in CTL,All_night, and All_day experiments during 0000 UTC-2300 UTC 4 June 2013. Averaged over a box area with a radius of 300 km around the vortex center

    3.3.2 穩(wěn)定度

    在對(duì)流系統(tǒng)中,云頂上部大量吸收短波輻射時(shí),可能會(huì)增加其靜力穩(wěn)定度,這將抑制垂直運(yùn)動(dòng)的強(qiáng)度。反之,長波輻射冷卻則會(huì)減弱靜力穩(wěn)定度,從而有利于對(duì)流的發(fā)展(Hobgood, 1986; Xu and Randall, 1995; Tang and Zhang, 2016)。圖15 給出了與輻射過程相關(guān)聯(lián)的凈加熱廓線(短波輻射和長波輻射之和,其中所選區(qū)域如圖3 綠色方框所示),以及夜間與控制試驗(yàn)和白天與控制試驗(yàn)之間低渦中心附近溫度場差值隨時(shí)間的演變特征。從圖15a 可以看到,CTL 試驗(yàn)中低層輻射冷卻較弱,最大冷卻出現(xiàn)在200 hPa 高度上。在All_night 試驗(yàn)中從低層到高層都有輻射冷卻,最大冷卻出現(xiàn)在250 hPa附近。而在All_day 試驗(yàn)中由輻射導(dǎo)致的加熱一直延伸到了200 hPa 附近,最大加熱在300 hPa 附近。從圖15b 和c 也可以看到,夜間試驗(yàn)相較于控制試驗(yàn),對(duì)流層中高層有明顯的冷異常,有利于對(duì)流的產(chǎn)生。相反,白天試驗(yàn)中對(duì)流層中高層增暖較大,使得大氣趨于穩(wěn)定狀態(tài),因此不利于對(duì)流的觸發(fā)。這結(jié)論與Melhauser and Zhang(2014)相一致。圖16 進(jìn)一步繪制了三組試驗(yàn)中?θ/?p的垂直廓線,以此反映靜力穩(wěn)定度程度。在白天4 日09:00,在CTL和All_day 試驗(yàn)中存在短波輻射的影響,從低層到高層均為靜力穩(wěn)定的。而在All_night 試驗(yàn)中,長波輻射冷卻降低了靜力穩(wěn)定度,在低層出現(xiàn)了靜力不穩(wěn)定的現(xiàn)象,這意味著對(duì)流的發(fā)展較快。到了夜間4 日15:00,CTL 試驗(yàn)中低層也呈現(xiàn)出靜力不穩(wěn)定的現(xiàn)象,但不穩(wěn)定的強(qiáng)度明顯小于All_night 試驗(yàn)。在4 日21:00,All_night 試驗(yàn)中低層不穩(wěn)定明顯減弱,對(duì)流強(qiáng)度減弱,這與水汽條件變差有關(guān),而該時(shí)段CTL 試驗(yàn)中低層不穩(wěn)定性進(jìn)一步加大,對(duì)流加強(qiáng)。類推可知,All_day 試驗(yàn)中,持續(xù)存在的短波輻射增加了大氣靜力穩(wěn)定度,這極大地抑制對(duì)流活動(dòng),最終導(dǎo)致低渦緩慢發(fā)展。

    圖15 2013 年6 月4 日00:00~23:00(a)三組試驗(yàn)中低渦中心附近區(qū)域(300 km×300 km),輻射導(dǎo)致的平均加熱廓線(單位:K h-1),低渦中心附近區(qū)域(300 km×300 km)平均溫度垂直分布隨時(shí)間變化的差異,(b)All_night 減CTL 試驗(yàn)(c)All_day 減 CTL 試驗(yàn)(等值線,單位:K, 實(shí)線為正值,虛線為負(fù)值)Fig. 15 (a) Vertical profiles of radiative heating (units: K h-1) in the three experiments, which are averaged over a box area centered at TPV during 0000-2300 UTC 4 June 2013. Time evolution of the vertical profile of temperature contrast (units: K, solid lines indicate positive value, dash lines indicate negative value) which are averaged over the box area centered at the TPV for (b) All_night minus CTL and (c) All_day minus CTL

    圖16 2013 年6 月4 日09:00、15:00 和21:00(a)CTL、(b)All_night 和(c)All_day 試驗(yàn)中低渦中心附近(300 km×300 km)區(qū)域平均的? θ/?p垂直廓線(單位:K hPa-1)分布Fig. 16 Vertical profiles of ? θ/?p(units: K hPa-1) at 0900 UTC, 1500 UTC, and 2100 UTC 4 June 2013, which are averaged over the box area centered at the TPV: (a) CTL experiment; (b) All_night experiment; (c) All_day experiment

    4 結(jié)論與討論

    為探究短波輻射日變化對(duì)高原渦發(fā)生的影響,本文利用WRF-ARW 模式對(duì)2013 年6 月4 日一次高原渦個(gè)例進(jìn)行模擬研究。通過設(shè)置不同太陽短波輻射強(qiáng)度的三組試驗(yàn),發(fā)現(xiàn)輻射日變化對(duì)高原低渦生成具有重要的影響。主要結(jié)論如下:

    (1)初始生成階段,高原渦在夜間發(fā)展的速率比白天更快,這與強(qiáng)對(duì)流系統(tǒng)在夜間更為活躍的觀測事實(shí)一致。從熱力學(xué)角度來看,在夜間,云頂長波輻射冷卻加強(qiáng)了對(duì)流層溫度遞減率,從而減弱了大氣靜力穩(wěn)定度,有利于對(duì)流層低層出現(xiàn)位勢不穩(wěn)定。同時(shí),由于大氣溫度的降低使得夜間相對(duì)濕度增大,進(jìn)而促使對(duì)流的發(fā)展,與此同時(shí),凝結(jié)潛熱的釋放會(huì)促進(jìn)和維持低渦的強(qiáng)度。反之,白天云頂吸收太陽短波輻射,使得對(duì)流層高層增溫大于地面,加強(qiáng)了大氣靜力穩(wěn)定度,從而抑制對(duì)流活動(dòng)發(fā)展。以上結(jié)論與輻射影響強(qiáng)對(duì)流日變化演變特征相符。

    (2)從數(shù)值試驗(yàn)結(jié)果表明,高原渦的發(fā)展與熱帶氣旋具有一定的相似性。作為中尺度系統(tǒng)高原渦,本文中該系統(tǒng)初始最大渦度出現(xiàn)在中層。從動(dòng)力學(xué)角度來看,其發(fā)展過程類似于熱帶氣旋“topdown”過程。而“渦度隔離”效應(yīng)對(duì)初始低渦的加強(qiáng)具有重要作用,小尺度強(qiáng)對(duì)流單體向中尺度低渦中心聚集,為高原渦的形成提供了擾動(dòng)渦度,這與熱帶氣旋生成中“積云熱塔”的作用類似。BKE 診斷結(jié)果表明,夜間低層輻合更為強(qiáng)盛,有利于上升運(yùn)動(dòng)的加強(qiáng)并誘發(fā)高原渦的生成。此外,從非平衡角度來看,在高原渦環(huán)流中心區(qū)域存在大值區(qū),這對(duì)應(yīng)中心的下沉運(yùn)動(dòng);而低渦四周為明顯的負(fù)值區(qū),對(duì)應(yīng)低渦周圍存在上升運(yùn)動(dòng)??偠灾?,高原渦具有暖心結(jié)構(gòu)特征,即低渦中心存在下沉運(yùn)動(dòng),而下沉絕熱增溫會(huì)加強(qiáng)其暖心,使得低渦發(fā)展。

    需要指出的是,前人研究指出云輻射—對(duì)流相互作用可以影響對(duì)流活動(dòng)(葛旭陽等, 2018; Ruppert et al., 2020),而本文僅針對(duì)一次高原渦個(gè)例進(jìn)行了模擬,此外云輻射—對(duì)流相互作用還沒有涉及,同時(shí)高原渦生成的背景環(huán)境場存在多樣性,這意味著不同背景環(huán)流場(如,切變線、輻合帶等)下低渦發(fā)展的機(jī)制可能存在差異。因此,下一步需要進(jìn)行更多的實(shí)際個(gè)例模擬,以得到較為完善的結(jié)論。

    猜你喜歡
    渦的渦度低層
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動(dòng)研究
    南海相對(duì)渦度的時(shí)空變化特征分析*
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    關(guān)于低層房屋建筑工程造價(jià)的要點(diǎn)及控制措施探討探索
    溝槽對(duì)湍流邊界層中展向渦影響的實(shí)驗(yàn)研究
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    開縫圓柱縫隙傾斜角對(duì)脫落渦的影響
    住八樓以上的人,早亡風(fēng)險(xiǎn)低
    益壽寶典(2017年34期)2017-02-26 08:27:20
    三角柱體對(duì)機(jī)翼根部馬蹄渦的影響
    低層高密度住宅設(shè)計(jì)探討
    精品久久久精品久久久| 窝窝影院91人妻| 麻豆一二三区av精品| 一二三四在线观看免费中文在| 欧美av亚洲av综合av国产av| 亚洲狠狠婷婷综合久久图片| 亚洲五月婷婷丁香| 精品第一国产精品| 动漫黄色视频在线观看| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看| 久久久国产成人精品二区| 男男h啪啪无遮挡| 久久婷婷成人综合色麻豆| 1024香蕉在线观看| 18禁裸乳无遮挡免费网站照片 | 久久午夜综合久久蜜桃| 一区二区三区激情视频| 久久中文看片网| 日韩av在线大香蕉| 午夜久久久在线观看| 国产免费男女视频| 亚洲欧美日韩无卡精品| 久久中文字幕一级| 成人国语在线视频| 在线十欧美十亚洲十日本专区| 亚洲伊人色综图| 国产精品一区二区三区四区久久 | 亚洲专区字幕在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品国产色婷婷电影| 美女扒开内裤让男人捅视频| 91麻豆av在线| 亚洲欧美精品综合久久99| www.精华液| 久久午夜综合久久蜜桃| 久久青草综合色| 精品卡一卡二卡四卡免费| 国产熟女xx| 无限看片的www在线观看| 岛国在线观看网站| 夜夜看夜夜爽夜夜摸| 色在线成人网| 两个人视频免费观看高清| 色婷婷久久久亚洲欧美| 国产三级在线视频| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 又黄又爽又免费观看的视频| 日本 av在线| 日本a在线网址| 亚洲成a人片在线一区二区| 国产亚洲精品av在线| 老司机靠b影院| 亚洲成av片中文字幕在线观看| 国产精华一区二区三区| 一级毛片精品| 欧美av亚洲av综合av国产av| 午夜福利在线观看吧| 亚洲国产中文字幕在线视频| 中文亚洲av片在线观看爽| 欧美激情 高清一区二区三区| 母亲3免费完整高清在线观看| 老司机福利观看| 欧美乱色亚洲激情| 久久婷婷人人爽人人干人人爱 | 亚洲成人久久性| 国产欧美日韩精品亚洲av| 如日韩欧美国产精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 一级黄色大片毛片| 久久中文看片网| 亚洲成av人片免费观看| 国产成人免费无遮挡视频| 成人av一区二区三区在线看| 国产精品自产拍在线观看55亚洲| 亚洲第一欧美日韩一区二区三区| 一级a爱视频在线免费观看| 两个人免费观看高清视频| 老鸭窝网址在线观看| 亚洲男人的天堂狠狠| 曰老女人黄片| 美女高潮喷水抽搐中文字幕| 18禁裸乳无遮挡免费网站照片 | 国产蜜桃级精品一区二区三区| 成人国语在线视频| 国产麻豆成人av免费视频| av视频在线观看入口| 日日干狠狠操夜夜爽| 最近最新中文字幕大全免费视频| cao死你这个sao货| 最近最新中文字幕大全免费视频| 桃红色精品国产亚洲av| 日本一区二区免费在线视频| 国产97色在线日韩免费| 国产高清videossex| 一级毛片精品| 国产精品乱码一区二三区的特点 | videosex国产| 精品无人区乱码1区二区| 成人永久免费在线观看视频| 久久狼人影院| 电影成人av| 国产高清videossex| 999久久久精品免费观看国产| 亚洲av成人av| 日韩一卡2卡3卡4卡2021年| 女人爽到高潮嗷嗷叫在线视频| 日本欧美视频一区| 1024视频免费在线观看| 国产精品野战在线观看| 在线视频色国产色| www.www免费av| 亚洲在线自拍视频| 亚洲伊人色综图| 男女之事视频高清在线观看| 欧美大码av| 999精品在线视频| 亚洲无线在线观看| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 免费高清视频大片| 免费在线观看亚洲国产| 欧美不卡视频在线免费观看 | 97人妻精品一区二区三区麻豆 | 免费在线观看黄色视频的| 久久久久久久久中文| 国产三级在线视频| 国产成人影院久久av| 欧美一级a爱片免费观看看 | 黑人巨大精品欧美一区二区蜜桃| 亚洲黑人精品在线| 精品第一国产精品| 波多野结衣av一区二区av| 亚洲国产精品成人综合色| 国产高清videossex| 精品一区二区三区视频在线观看免费| 日韩视频一区二区在线观看| 在线av久久热| 操美女的视频在线观看| 国产免费av片在线观看野外av| 亚洲全国av大片| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 午夜两性在线视频| 精品高清国产在线一区| 日本黄色视频三级网站网址| av在线播放免费不卡| 香蕉丝袜av| bbb黄色大片| 夜夜躁狠狠躁天天躁| av免费在线观看网站| 巨乳人妻的诱惑在线观看| 黑丝袜美女国产一区| 69精品国产乱码久久久| 久久人妻av系列| 欧美日韩精品网址| 搡老岳熟女国产| 亚洲狠狠婷婷综合久久图片| 欧美日韩福利视频一区二区| 9色porny在线观看| 亚洲人成网站在线播放欧美日韩| av视频免费观看在线观看| 免费女性裸体啪啪无遮挡网站| av在线播放免费不卡| 国产欧美日韩精品亚洲av| 91av网站免费观看| 亚洲人成电影免费在线| √禁漫天堂资源中文www| 国产99久久九九免费精品| 少妇 在线观看| 亚洲国产精品999在线| 长腿黑丝高跟| 1024视频免费在线观看| 成年女人毛片免费观看观看9| 热re99久久国产66热| 亚洲电影在线观看av| 国产xxxxx性猛交| 俄罗斯特黄特色一大片| 18禁黄网站禁片午夜丰满| 免费久久久久久久精品成人欧美视频| 成人欧美大片| 在线观看免费视频日本深夜| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久精品电影 | 人人妻人人爽人人添夜夜欢视频| 啪啪无遮挡十八禁网站| 欧美一级毛片孕妇| 好看av亚洲va欧美ⅴa在| 天天躁夜夜躁狠狠躁躁| 91精品三级在线观看| 母亲3免费完整高清在线观看| 精品欧美一区二区三区在线| 亚洲欧美精品综合久久99| 嫁个100分男人电影在线观看| 中文字幕最新亚洲高清| 亚洲电影在线观看av| 亚洲熟妇中文字幕五十中出| 黄片大片在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 美女国产高潮福利片在线看| avwww免费| АⅤ资源中文在线天堂| 午夜免费激情av| 亚洲一区高清亚洲精品| 亚洲黑人精品在线| 免费搜索国产男女视频| 国产av又大| 黄片播放在线免费| 亚洲成人免费电影在线观看| 午夜免费观看网址| 18美女黄网站色大片免费观看| 久久伊人香网站| 国产乱人伦免费视频| 国产蜜桃级精品一区二区三区| av视频在线观看入口| 亚洲成国产人片在线观看| 免费在线观看视频国产中文字幕亚洲| av天堂在线播放| 女性生殖器流出的白浆| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 亚洲三区欧美一区| 国产一卡二卡三卡精品| 亚洲人成电影观看| 九色国产91popny在线| ponron亚洲| а√天堂www在线а√下载| 久久人人爽av亚洲精品天堂| 精品人妻1区二区| 日本a在线网址| 欧美另类亚洲清纯唯美| 老熟妇乱子伦视频在线观看| 亚洲精品久久成人aⅴ小说| 两个人视频免费观看高清| 国产激情久久老熟女| 九色亚洲精品在线播放| 一区福利在线观看| 又大又爽又粗| 最近最新中文字幕大全电影3 | 免费观看精品视频网站| 成人av一区二区三区在线看| 成人免费观看视频高清| 人人妻人人澡人人看| 9191精品国产免费久久| 波多野结衣高清无吗| 国产一区二区三区在线臀色熟女| 性少妇av在线| 亚洲午夜理论影院| 亚洲欧美日韩无卡精品| 国产成人精品久久二区二区91| 亚洲avbb在线观看| 又紧又爽又黄一区二区| 美女高潮喷水抽搐中文字幕| 操出白浆在线播放| 亚洲午夜理论影院| 国产xxxxx性猛交| 久久天躁狠狠躁夜夜2o2o| 国产av一区在线观看免费| 亚洲精品在线美女| АⅤ资源中文在线天堂| 美女 人体艺术 gogo| 精品久久久精品久久久| 国产精品香港三级国产av潘金莲| 免费人成视频x8x8入口观看| 欧美中文综合在线视频| 又紧又爽又黄一区二区| 免费高清视频大片| 操美女的视频在线观看| 欧美 亚洲 国产 日韩一| 久久人人精品亚洲av| 久久人妻熟女aⅴ| svipshipincom国产片| 亚洲中文字幕日韩| 在线天堂中文资源库| 非洲黑人性xxxx精品又粗又长| av电影中文网址| 午夜福利欧美成人| 久久久久九九精品影院| 国产激情欧美一区二区| 人人妻,人人澡人人爽秒播| 亚洲精品中文字幕一二三四区| svipshipincom国产片| 999久久久精品免费观看国产| 亚洲中文av在线| 欧美黄色淫秽网站| 少妇 在线观看| 国产一区二区三区综合在线观看| 极品人妻少妇av视频| 国产免费男女视频| av视频免费观看在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲精品国产精品久久久不卡| 18禁裸乳无遮挡免费网站照片 | 久久精品国产亚洲av高清一级| 久久久国产成人精品二区| 丁香六月欧美| 亚洲美女黄片视频| 老鸭窝网址在线观看| 黄色a级毛片大全视频| 自线自在国产av| 国产成人精品久久二区二区91| 伦理电影免费视频| 男人舔女人的私密视频| 丝袜美腿诱惑在线| 国产精品永久免费网站| 12—13女人毛片做爰片一| 国内毛片毛片毛片毛片毛片| 国产97色在线日韩免费| 亚洲专区国产一区二区| 亚洲精品国产精品久久久不卡| 久久久久国产一级毛片高清牌| 精品国产亚洲在线| 免费观看人在逋| 天堂影院成人在线观看| 村上凉子中文字幕在线| 最近最新免费中文字幕在线| 日日摸夜夜添夜夜添小说| 欧美激情 高清一区二区三区| 日本三级黄在线观看| 欧美av亚洲av综合av国产av| 欧美中文日本在线观看视频| 国产成人系列免费观看| 夜夜夜夜夜久久久久| 亚洲一区二区三区色噜噜| 国产精品 国内视频| 在线十欧美十亚洲十日本专区| 极品教师在线免费播放| 亚洲精品av麻豆狂野| 亚洲熟妇中文字幕五十中出| 怎么达到女性高潮| 一二三四社区在线视频社区8| 国产av一区在线观看免费| 天天躁狠狠躁夜夜躁狠狠躁| 午夜精品在线福利| 久久香蕉激情| www.www免费av| 免费久久久久久久精品成人欧美视频| 亚洲av片天天在线观看| 日本欧美视频一区| 欧美成人一区二区免费高清观看 | 91大片在线观看| 此物有八面人人有两片| 母亲3免费完整高清在线观看| 日韩大尺度精品在线看网址 | av在线天堂中文字幕| 老汉色av国产亚洲站长工具| 丝袜人妻中文字幕| 久久 成人 亚洲| avwww免费| 可以在线观看的亚洲视频| 久久精品亚洲熟妇少妇任你| 国产精品免费一区二区三区在线| 制服诱惑二区| 少妇粗大呻吟视频| 热99re8久久精品国产| 天天躁夜夜躁狠狠躁躁| aaaaa片日本免费| 天天躁夜夜躁狠狠躁躁| 日韩一卡2卡3卡4卡2021年| 国产精品久久久久久亚洲av鲁大| 国内精品久久久久久久电影| av天堂在线播放| aaaaa片日本免费| 人人妻,人人澡人人爽秒播| 国产麻豆69| 男人舔女人的私密视频| 99久久精品国产亚洲精品| 午夜免费成人在线视频| 国产区一区二久久| 国产一区在线观看成人免费| 美女扒开内裤让男人捅视频| 一进一出好大好爽视频| 国产成人av激情在线播放| www.999成人在线观看| 久久久久久久久免费视频了| 啦啦啦 在线观看视频| 啦啦啦观看免费观看视频高清 | 99久久99久久久精品蜜桃| av有码第一页| 好看av亚洲va欧美ⅴa在| 香蕉丝袜av| 99国产精品一区二区三区| 久久精品国产亚洲av高清一级| 亚洲最大成人中文| 久久天躁狠狠躁夜夜2o2o| 亚洲五月色婷婷综合| 黄片大片在线免费观看| 国产午夜福利久久久久久| tocl精华| 嫩草影视91久久| aaaaa片日本免费| 亚洲美女黄片视频| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片 | 精品人妻在线不人妻| 国产av一区在线观看免费| 宅男免费午夜| 日日夜夜操网爽| 久久久久久大精品| 好看av亚洲va欧美ⅴa在| 黄频高清免费视频| av免费在线观看网站| a在线观看视频网站| 亚洲欧美激情在线| 国产欧美日韩一区二区精品| 亚洲全国av大片| 免费在线观看亚洲国产| 啦啦啦免费观看视频1| 日韩欧美在线二视频| 在线观看一区二区三区| 精品欧美一区二区三区在线| 男女午夜视频在线观看| av福利片在线| 在线观看舔阴道视频| 天堂√8在线中文| √禁漫天堂资源中文www| 一进一出抽搐gif免费好疼| 久久久久久久久久久久大奶| 男男h啪啪无遮挡| 久久狼人影院| 精品一区二区三区四区五区乱码| 亚洲 欧美 日韩 在线 免费| 欧美日韩亚洲综合一区二区三区_| 50天的宝宝边吃奶边哭怎么回事| av欧美777| 色综合欧美亚洲国产小说| 婷婷精品国产亚洲av在线| 精品免费久久久久久久清纯| 老熟妇乱子伦视频在线观看| 青草久久国产| 国产精品影院久久| 成熟少妇高潮喷水视频| 亚洲,欧美精品.| 久久青草综合色| 国产一区二区三区在线臀色熟女| 757午夜福利合集在线观看| 中亚洲国语对白在线视频| 久久香蕉精品热| 一级作爱视频免费观看| 亚洲伊人色综图| 精品午夜福利视频在线观看一区| 国产一卡二卡三卡精品| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 午夜激情av网站| 成人手机av| 欧美+亚洲+日韩+国产| 少妇的丰满在线观看| 在线观看一区二区三区| 91老司机精品| 久久久国产精品麻豆| 免费一级毛片在线播放高清视频 | 正在播放国产对白刺激| 久久久水蜜桃国产精品网| 国产在线精品亚洲第一网站| videosex国产| 岛国在线观看网站| 两性夫妻黄色片| 黑人巨大精品欧美一区二区蜜桃| 成人亚洲精品av一区二区| 宅男免费午夜| 亚洲av片天天在线观看| 夜夜爽天天搞| 99久久久亚洲精品蜜臀av| av片东京热男人的天堂| 久久精品国产亚洲av香蕉五月| 国内精品久久久久精免费| 久久久精品欧美日韩精品| 不卡av一区二区三区| 丝袜人妻中文字幕| 麻豆成人av在线观看| 亚洲精华国产精华精| 又黄又爽又免费观看的视频| 亚洲七黄色美女视频| 国产成+人综合+亚洲专区| 亚洲人成77777在线视频| 久久 成人 亚洲| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 午夜福利高清视频| 男人操女人黄网站| 黄网站色视频无遮挡免费观看| 桃红色精品国产亚洲av| 国产精品1区2区在线观看.| 两个人视频免费观看高清| 亚洲自偷自拍图片 自拍| 中文字幕另类日韩欧美亚洲嫩草| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频| 91大片在线观看| 久久中文看片网| 欧美久久黑人一区二区| 美女高潮喷水抽搐中文字幕| 制服人妻中文乱码| 在线永久观看黄色视频| 两个人看的免费小视频| 精品乱码久久久久久99久播| 女人被狂操c到高潮| 亚洲人成伊人成综合网2020| 欧美精品亚洲一区二区| 免费人成视频x8x8入口观看| 国产激情欧美一区二区| videosex国产| 香蕉久久夜色| 国产野战对白在线观看| 国产麻豆69| 亚洲美女黄片视频| 大码成人一级视频| av在线播放免费不卡| 黄色丝袜av网址大全| 天堂动漫精品| 中文字幕最新亚洲高清| 黄片小视频在线播放| 精品国产乱子伦一区二区三区| 日韩三级视频一区二区三区| 真人做人爱边吃奶动态| 国内毛片毛片毛片毛片毛片| 国产精品野战在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产主播在线观看一区二区| 欧美一级毛片孕妇| 精品国产美女av久久久久小说| 女人被躁到高潮嗷嗷叫费观| 在线永久观看黄色视频| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| 夜夜躁狠狠躁天天躁| 国产精品日韩av在线免费观看 | 露出奶头的视频| 亚洲精品av麻豆狂野| 亚洲av美国av| 色综合婷婷激情| 国产一区二区在线av高清观看| 亚洲精品av麻豆狂野| 俄罗斯特黄特色一大片| 亚洲午夜精品一区,二区,三区| 又大又爽又粗| ponron亚洲| 成人国语在线视频| 久久狼人影院| 级片在线观看| 国产欧美日韩精品亚洲av| 波多野结衣巨乳人妻| 亚洲九九香蕉| 9热在线视频观看99| 1024视频免费在线观看| 精品一区二区三区视频在线观看免费| 久久性视频一级片| 纯流量卡能插随身wifi吗| 两人在一起打扑克的视频| 亚洲国产精品sss在线观看| www.精华液| 欧美黑人欧美精品刺激| 久久 成人 亚洲| 成在线人永久免费视频| 99久久99久久久精品蜜桃| 久久久精品欧美日韩精品| 久久久久久久久中文| 国产精品一区二区三区四区久久 | 一个人观看的视频www高清免费观看 | 亚洲欧美激情在线| 欧美成人一区二区免费高清观看 | 国产乱人伦免费视频| 18美女黄网站色大片免费观看| 激情视频va一区二区三区| 精品熟女少妇八av免费久了| 国产精品影院久久| 国产麻豆69| 精品久久久久久,| 成人亚洲精品av一区二区| 午夜福利一区二区在线看| 色尼玛亚洲综合影院| 1024香蕉在线观看| 窝窝影院91人妻| 国产精品日韩av在线免费观看 | 久久中文字幕一级| 老熟妇乱子伦视频在线观看| www.www免费av| ponron亚洲| 性欧美人与动物交配| 男女做爰动态图高潮gif福利片 | 99在线人妻在线中文字幕| 精品福利观看| 日韩一卡2卡3卡4卡2021年| 亚洲男人天堂网一区| 欧美性长视频在线观看| 国产国语露脸激情在线看| 久久天堂一区二区三区四区| 搞女人的毛片| 99精品久久久久人妻精品| 他把我摸到了高潮在线观看| 亚洲自偷自拍图片 自拍| 国产精品美女特级片免费视频播放器 | 午夜精品在线福利| 午夜a级毛片| 91九色精品人成在线观看| 一区在线观看完整版| 国产免费av片在线观看野外av| 欧美乱色亚洲激情| 中文字幕人妻熟女乱码| 在线观看免费视频网站a站| 韩国av一区二区三区四区| 91国产中文字幕| 亚洲国产看品久久| 9191精品国产免费久久| 精品欧美国产一区二区三| 欧美黑人精品巨大| 精品高清国产在线一区| 黄片大片在线免费观看| 午夜福利欧美成人| 91在线观看av| 麻豆国产av国片精品| 精品人妻在线不人妻| 熟妇人妻久久中文字幕3abv|