• <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ì)探討
    两个人视频免费观看高清| 三级毛片av免费| 免费在线观看亚洲国产| 波多野结衣巨乳人妻| 在线观看www视频免费| 巨乳人妻的诱惑在线观看| 国产一区二区三区在线臀色熟女| 欧美色欧美亚洲另类二区| 精品久久久久久久久久免费视频| 身体一侧抽搐| 亚洲熟妇熟女久久| 99在线视频只有这里精品首页| 日本三级黄在线观看| 亚洲最大成人中文| 村上凉子中文字幕在线| 久久中文字幕一级| 可以免费在线观看a视频的电影网站| 精品久久久久久久人妻蜜臀av| 欧美亚洲日本最大视频资源| 亚洲精品美女久久久久99蜜臀| 亚洲av熟女| 99热6这里只有精品| 精品国内亚洲2022精品成人| 女同久久另类99精品国产91| 国产私拍福利视频在线观看| 午夜精品在线福利| 中文字幕av电影在线播放| 国产99白浆流出| 久久99热这里只有精品18| 国产成人精品无人区| 首页视频小说图片口味搜索| www日本在线高清视频| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久男人| x7x7x7水蜜桃| av电影中文网址| 97碰自拍视频| 少妇裸体淫交视频免费看高清 | 国产日本99.免费观看| 热re99久久国产66热| 一进一出抽搐gif免费好疼| 色综合亚洲欧美另类图片| 天堂√8在线中文| 亚洲男人天堂网一区| 在线免费观看的www视频| 国产精品野战在线观看| 高潮久久久久久久久久久不卡| 美国免费a级毛片| 成人av一区二区三区在线看| av中文乱码字幕在线| 在线永久观看黄色视频| 黄色视频不卡| 最近最新中文字幕大全免费视频| 国产色视频综合| 欧美一级毛片孕妇| 色尼玛亚洲综合影院| 色播在线永久视频| 麻豆成人av在线观看| 亚洲中文字幕日韩| 亚洲午夜精品一区,二区,三区| 一边摸一边做爽爽视频免费| 天天躁夜夜躁狠狠躁躁| 在线视频色国产色| 国产欧美日韩一区二区三| 亚洲无线在线观看| 精品久久久久久久末码| 18禁黄网站禁片午夜丰满| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美日韩高清在线视频| 无限看片的www在线观看| 久久香蕉精品热| 村上凉子中文字幕在线| 51午夜福利影视在线观看| 国产成人欧美在线观看| 亚洲av成人一区二区三| 99国产精品一区二区三区| 日本熟妇午夜| 国产精品乱码一区二三区的特点| 亚洲一区中文字幕在线| 麻豆久久精品国产亚洲av| 国产三级在线视频| 国产精品久久久人人做人人爽| 欧美av亚洲av综合av国产av| 男女午夜视频在线观看| 男女之事视频高清在线观看| 国产高清videossex| 女同久久另类99精品国产91| 久久久水蜜桃国产精品网| av福利片在线| 看片在线看免费视频| 日本免费一区二区三区高清不卡| 香蕉丝袜av| 欧美丝袜亚洲另类 | 亚洲欧美精品综合久久99| 亚洲一区中文字幕在线| 深夜精品福利| 男人操女人黄网站| 午夜久久久在线观看| 美女 人体艺术 gogo| 人人妻人人澡人人看| 99久久精品国产亚洲精品| 亚洲av片天天在线观看| 999精品在线视频| 国产一卡二卡三卡精品| 国产精品av久久久久免费| 男人的好看免费观看在线视频 | x7x7x7水蜜桃| 久久人妻av系列| 日韩精品中文字幕看吧| 2021天堂中文幕一二区在线观 | 韩国av一区二区三区四区| 不卡一级毛片| 日韩欧美一区视频在线观看| 一二三四社区在线视频社区8| 亚洲七黄色美女视频| 高清毛片免费观看视频网站| 日本一区二区免费在线视频| 老司机午夜福利在线观看视频| 亚洲av中文字字幕乱码综合 | 午夜福利成人在线免费观看| 精品国产亚洲在线| 免费看美女性在线毛片视频| 日韩欧美国产在线观看| 好男人电影高清在线观看| 国产精品av久久久久免费| 特大巨黑吊av在线直播 | 别揉我奶头~嗯~啊~动态视频| 国产1区2区3区精品| 日本五十路高清| 狠狠狠狠99中文字幕| 午夜成年电影在线免费观看| 午夜福利高清视频| 亚洲av熟女| 日韩欧美国产在线观看| 18禁观看日本| 男女之事视频高清在线观看| 精品国产乱子伦一区二区三区| 好男人电影高清在线观看| 国产亚洲精品久久久久5区| www.精华液| 欧美三级亚洲精品| 亚洲性夜色夜夜综合| 午夜老司机福利片| 中文字幕高清在线视频| 亚洲精品久久成人aⅴ小说| 久久性视频一级片| 一本大道久久a久久精品| 日韩欧美一区视频在线观看| 久久久精品国产亚洲av高清涩受| 国产黄片美女视频| 久久久国产精品麻豆| 精品久久蜜臀av无| 亚洲中文字幕一区二区三区有码在线看 | 岛国在线观看网站| 啦啦啦 在线观看视频| 免费在线观看亚洲国产| 成人国产综合亚洲| 国产亚洲欧美在线一区二区| 亚洲国产精品久久男人天堂| 此物有八面人人有两片| 欧美黑人巨大hd| 怎么达到女性高潮| 精品国产一区二区三区四区第35| 在线观看免费午夜福利视频| 黄色片一级片一级黄色片| 成人永久免费在线观看视频| 免费在线观看完整版高清| 国产国语露脸激情在线看| 欧美绝顶高潮抽搐喷水| 国产精品久久久人人做人人爽| 一级作爱视频免费观看| 777久久人妻少妇嫩草av网站| 天堂√8在线中文| 国产精品亚洲av一区麻豆| 久热爱精品视频在线9| 宅男免费午夜| 欧美黄色淫秽网站| 国产一级毛片七仙女欲春2 | 国产一级毛片七仙女欲春2 | 一区二区三区高清视频在线| 免费看美女性在线毛片视频| 精品国内亚洲2022精品成人| 亚洲第一av免费看| 国产极品粉嫩免费观看在线| 国产精品98久久久久久宅男小说| 亚洲人成77777在线视频| 亚洲国产中文字幕在线视频| 亚洲成人久久性| 少妇裸体淫交视频免费看高清 | 热re99久久国产66热| 国内揄拍国产精品人妻在线 | 亚洲中文字幕日韩| 欧美成人性av电影在线观看| 丁香六月欧美| 少妇被粗大的猛进出69影院| av片东京热男人的天堂| 欧美亚洲日本最大视频资源| 51午夜福利影视在线观看| 91老司机精品| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 婷婷亚洲欧美| 亚洲九九香蕉| 欧美最黄视频在线播放免费| 久久99热这里只有精品18| 热99re8久久精品国产| 激情在线观看视频在线高清| 黄色视频不卡| svipshipincom国产片| 18禁观看日本| 成人手机av| av电影中文网址| 欧美黄色片欧美黄色片| 午夜成年电影在线免费观看| 色在线成人网| 91国产中文字幕| 中文字幕精品免费在线观看视频| 女生性感内裤真人,穿戴方法视频| 午夜福利成人在线免费观看| 2021天堂中文幕一二区在线观 | 免费看a级黄色片| 嫁个100分男人电影在线观看| 亚洲一区二区三区不卡视频| 一二三四社区在线视频社区8| 欧美成人免费av一区二区三区| 日韩欧美 国产精品| 亚洲七黄色美女视频| 黄色a级毛片大全视频| 天天躁狠狠躁夜夜躁狠狠躁| 麻豆国产av国片精品| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 国产成人欧美在线观看| 亚洲 国产 在线| 欧美zozozo另类| 可以在线观看毛片的网站| 在线观看一区二区三区| x7x7x7水蜜桃| 悠悠久久av| 日本a在线网址| 黄色a级毛片大全视频| 国产精品久久久av美女十八| 亚洲成人久久爱视频| 午夜激情av网站| 久久久久久大精品| 男女那种视频在线观看| 亚洲精品久久成人aⅴ小说| 日本一区二区免费在线视频| 精品国产超薄肉色丝袜足j| 中文字幕另类日韩欧美亚洲嫩草| 人人澡人人妻人| 搡老岳熟女国产| 欧美人与性动交α欧美精品济南到| 99国产精品一区二区三区| 亚洲精品美女久久久久99蜜臀| 精品欧美一区二区三区在线| 亚洲全国av大片| 18美女黄网站色大片免费观看| 午夜日韩欧美国产| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 在线观看一区二区三区| 欧美中文日本在线观看视频| 久久久久久九九精品二区国产 | 高清在线国产一区| 巨乳人妻的诱惑在线观看| 精品熟女少妇八av免费久了| 午夜两性在线视频| 一级片免费观看大全| 9191精品国产免费久久| 久久精品影院6| 午夜久久久在线观看| 男女午夜视频在线观看| 99热6这里只有精品| 精品一区二区三区视频在线观看免费| 99久久无色码亚洲精品果冻| 亚洲va日本ⅴa欧美va伊人久久| 亚洲男人天堂网一区| 桃红色精品国产亚洲av| 亚洲精品中文字幕一二三四区| 午夜亚洲福利在线播放| 中出人妻视频一区二区| 久久久久国产精品人妻aⅴ院| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91| 欧美另类亚洲清纯唯美| 国产又爽黄色视频| 在线av久久热| 午夜视频精品福利| 男人的好看免费观看在线视频 | 999久久久精品免费观看国产| 日日干狠狠操夜夜爽| 伦理电影免费视频| 国产真人三级小视频在线观看| 久久久久久久久免费视频了| 大型黄色视频在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 亚洲一区二区三区不卡视频| 成人亚洲精品一区在线观看| 免费在线观看日本一区| 亚洲熟妇中文字幕五十中出| 国内精品久久久久精免费| 少妇被粗大的猛进出69影院| 免费在线观看亚洲国产| 亚洲专区中文字幕在线| 在线观看日韩欧美| 亚洲人成网站在线播放欧美日韩| 国产黄色小视频在线观看| 午夜日韩欧美国产| 亚洲狠狠婷婷综合久久图片| 中文字幕av电影在线播放| 亚洲精品色激情综合| 国产片内射在线| 欧美在线黄色| 女人被狂操c到高潮| 韩国精品一区二区三区| 精品一区二区三区四区五区乱码| 999精品在线视频| 黄片大片在线免费观看| 亚洲精品一区av在线观看| 久久九九热精品免费| 国产精品98久久久久久宅男小说| av超薄肉色丝袜交足视频| 日本熟妇午夜| 亚洲aⅴ乱码一区二区在线播放 | 2021天堂中文幕一二区在线观 | 亚洲精品久久成人aⅴ小说| 老司机午夜福利在线观看视频| av中文乱码字幕在线| 国产麻豆成人av免费视频| 成人av一区二区三区在线看| 久99久视频精品免费| 身体一侧抽搐| 99精品久久久久人妻精品| 亚洲 欧美 日韩 在线 免费| 欧美乱妇无乱码| 亚洲最大成人中文| 69av精品久久久久久| 老司机深夜福利视频在线观看| 叶爱在线成人免费视频播放| 熟妇人妻久久中文字幕3abv| 日韩视频一区二区在线观看| 国产熟女xx| 国产精品自产拍在线观看55亚洲| netflix在线观看网站| av有码第一页| 99久久99久久久精品蜜桃| 久久精品影院6| 夜夜躁狠狠躁天天躁| 亚洲av美国av| 91老司机精品| 日本精品一区二区三区蜜桃| 啦啦啦韩国在线观看视频| 免费高清视频大片| 脱女人内裤的视频| 女性生殖器流出的白浆| 窝窝影院91人妻| 国产午夜精品久久久久久| 免费在线观看日本一区| 亚洲中文av在线| 99精品在免费线老司机午夜| 亚洲精品国产精品久久久不卡| 国产在线精品亚洲第一网站| 亚洲成av片中文字幕在线观看| 久久久久久久久免费视频了| 久久久水蜜桃国产精品网| 精品国内亚洲2022精品成人| 色综合站精品国产| 99久久99久久久精品蜜桃| 亚洲av电影在线进入| 午夜日韩欧美国产| 香蕉丝袜av| 欧美绝顶高潮抽搐喷水| 色播亚洲综合网| 精品福利观看| 日韩大码丰满熟妇| 熟女电影av网| 制服诱惑二区| 成人国产一区最新在线观看| xxxwww97欧美| 亚洲国产毛片av蜜桃av| 可以免费在线观看a视频的电影网站| 欧美av亚洲av综合av国产av| 欧美日韩亚洲国产一区二区在线观看| 国产成人系列免费观看| 特大巨黑吊av在线直播 | 脱女人内裤的视频| 90打野战视频偷拍视频| 国产熟女xx| 少妇裸体淫交视频免费看高清 | 成人特级黄色片久久久久久久| 一本久久中文字幕| 午夜视频精品福利| 在线看三级毛片| 村上凉子中文字幕在线| 校园春色视频在线观看| 亚洲精品久久成人aⅴ小说| 老司机福利观看| 欧美黑人精品巨大| 成熟少妇高潮喷水视频| 欧美黑人精品巨大| 国产aⅴ精品一区二区三区波| 88av欧美| 中文字幕av电影在线播放| 国产一区二区激情短视频| 一区二区三区国产精品乱码| 90打野战视频偷拍视频| 一本综合久久免费| 久久久精品欧美日韩精品| 中文字幕高清在线视频| 午夜亚洲福利在线播放| 国产精品99久久99久久久不卡| 又大又爽又粗| 999精品在线视频| 一进一出抽搐动态| 亚洲国产精品999在线| 久久草成人影院| 中文字幕人成人乱码亚洲影| 亚洲国产中文字幕在线视频| 久久精品国产亚洲av香蕉五月| 法律面前人人平等表现在哪些方面| 国产蜜桃级精品一区二区三区| 一本精品99久久精品77| 亚洲成av片中文字幕在线观看| 精品国产乱子伦一区二区三区| 久久精品国产亚洲av香蕉五月| 国产成人啪精品午夜网站| 国产一卡二卡三卡精品| 亚洲精品av麻豆狂野| 久久久精品欧美日韩精品| 中文字幕精品免费在线观看视频| 久久亚洲真实| 欧美 亚洲 国产 日韩一| 99久久无色码亚洲精品果冻| 精品国产亚洲在线| 中文字幕精品亚洲无线码一区 | 午夜精品在线福利| 99国产极品粉嫩在线观看| 久久精品国产99精品国产亚洲性色| 国产精品香港三级国产av潘金莲| 久久午夜综合久久蜜桃| 国产黄a三级三级三级人| 在线国产一区二区在线| 99国产精品一区二区蜜桃av| 精品国产亚洲在线| 亚洲人成伊人成综合网2020| 美女午夜性视频免费| 欧美国产日韩亚洲一区| 天堂动漫精品| 亚洲国产日韩欧美精品在线观看 | 亚洲国产高清在线一区二区三 | 午夜两性在线视频| 黄色片一级片一级黄色片| 亚洲国产毛片av蜜桃av| 欧美大码av| 国产黄a三级三级三级人| 好看av亚洲va欧美ⅴa在| 色播亚洲综合网| 精品国产超薄肉色丝袜足j| 国产极品粉嫩免费观看在线| 一级a爱片免费观看的视频| 丰满人妻熟妇乱又伦精品不卡| 日本精品一区二区三区蜜桃| 91九色精品人成在线观看| 日本a在线网址| 国产精品亚洲一级av第二区| bbb黄色大片| 久久久国产欧美日韩av| 欧美国产日韩亚洲一区| 久久久水蜜桃国产精品网| 美国免费a级毛片| a级毛片在线看网站| 久久婷婷成人综合色麻豆| 国产成人一区二区三区免费视频网站| 啦啦啦免费观看视频1| 国产aⅴ精品一区二区三区波| 国产一区二区激情短视频| 老司机深夜福利视频在线观看| 男女视频在线观看网站免费 | 国产午夜精品久久久久久| 熟妇人妻久久中文字幕3abv| 高清毛片免费观看视频网站| 欧美又色又爽又黄视频| 国产久久久一区二区三区| 亚洲av成人av| 久久精品亚洲精品国产色婷小说| 这个男人来自地球电影免费观看| 伊人久久大香线蕉亚洲五| 欧美成人免费av一区二区三区| 黄片小视频在线播放| 国产成人精品久久二区二区免费| 中文字幕另类日韩欧美亚洲嫩草| 中文亚洲av片在线观看爽| 非洲黑人性xxxx精品又粗又长| 国产男靠女视频免费网站| 国产久久久一区二区三区| 变态另类丝袜制服| 日韩国内少妇激情av| 一区二区三区精品91| 成人国产一区最新在线观看| 欧美中文综合在线视频| 久久久水蜜桃国产精品网| 欧美绝顶高潮抽搐喷水| 成人手机av| 国产三级在线视频| 最新美女视频免费是黄的| or卡值多少钱| 麻豆成人av在线观看| 午夜老司机福利片| 午夜福利一区二区在线看| 麻豆一二三区av精品| 国产精品久久久久久人妻精品电影| 久久久久久免费高清国产稀缺| 免费在线观看成人毛片| 欧美午夜高清在线| 欧美日本视频| 午夜激情av网站| 别揉我奶头~嗯~啊~动态视频| 国产在线观看jvid| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久5区| 一区二区三区激情视频| 最新在线观看一区二区三区| 欧美激情 高清一区二区三区| 国产亚洲精品av在线| 中国美女看黄片| 午夜福利在线在线| 正在播放国产对白刺激| or卡值多少钱| 亚洲人成网站在线播放欧美日韩| 两人在一起打扑克的视频| 久久天堂一区二区三区四区| 亚洲av成人一区二区三| 这个男人来自地球电影免费观看| 女人高潮潮喷娇喘18禁视频| 成人18禁高潮啪啪吃奶动态图| 精品久久久久久久毛片微露脸| 一级黄色大片毛片| 亚洲人成伊人成综合网2020| 在线观看一区二区三区| 真人做人爱边吃奶动态| 精品久久久久久久末码| av天堂在线播放| 在线播放国产精品三级| 日韩视频一区二区在线观看| 亚洲三区欧美一区| www国产在线视频色| 亚洲av成人av| 激情在线观看视频在线高清| bbb黄色大片| 国产蜜桃级精品一区二区三区| 亚洲av片天天在线观看| 久久久久九九精品影院| 亚洲美女黄片视频| 怎么达到女性高潮| 啪啪无遮挡十八禁网站| 韩国精品一区二区三区| 欧美大码av| 欧美日韩乱码在线| 在线观看www视频免费| 精品久久蜜臀av无| 国产午夜精品久久久久久| 国产黄片美女视频| 国产日本99.免费观看| 岛国视频午夜一区免费看| 天天一区二区日本电影三级| 欧美日韩福利视频一区二区| 男女做爰动态图高潮gif福利片| 午夜日韩欧美国产| 午夜福利视频1000在线观看| 国产又爽黄色视频| 真人做人爱边吃奶动态| 18禁裸乳无遮挡免费网站照片 | 国产成人一区二区三区免费视频网站| 国产欧美日韩一区二区三| 国产成+人综合+亚洲专区| 久9热在线精品视频| av欧美777| 精品高清国产在线一区| 又大又爽又粗| 国产成人系列免费观看| 日日摸夜夜添夜夜添小说| 在线观看午夜福利视频| 村上凉子中文字幕在线| 精品国内亚洲2022精品成人| av超薄肉色丝袜交足视频| 男女午夜视频在线观看| 中文字幕久久专区| 亚洲精华国产精华精| 国产野战对白在线观看| 欧美最黄视频在线播放免费| 久久天躁狠狠躁夜夜2o2o| 亚洲自拍偷在线| 可以在线观看毛片的网站| 在线观看午夜福利视频| 97人妻精品一区二区三区麻豆 | 在线免费观看的www视频| 色婷婷久久久亚洲欧美| 99热6这里只有精品| 亚洲成av片中文字幕在线观看| 成年版毛片免费区| 国产视频内射| 天堂√8在线中文| 在线观看66精品国产| 免费看日本二区| 欧美日韩中文字幕国产精品一区二区三区| 女同久久另类99精品国产91| 精品欧美一区二区三区在线| 久久久久久大精品| 国产一级毛片七仙女欲春2 | av有码第一页| 成人特级黄色片久久久久久久| 欧美日韩亚洲国产一区二区在线观看|