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

    1521 號(hào)臺(tái)風(fēng)“杜鵑”生成階段中尺度對(duì)流系統(tǒng)和降水演變特征

    2021-04-16 07:37:06王科陳光華呂欣宇李興良
    大氣科學(xué) 2021年1期
    關(guān)鍵詞:層云渦度杜鵑

    王科 陳光華 呂欣宇 ,3 李興良

    1 中國(guó)科學(xué)院大氣物理研究所云降水物理與強(qiáng)風(fēng)暴實(shí)驗(yàn)室,北京100029

    2 天津海洋中心氣象臺(tái),天津 300074

    3 云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院,昆明 650500

    4 國(guó)家氣象中心(中央氣象臺(tái)),北京 100081

    1 引言

    熱帶氣旋(Tropical Cyclone,簡(jiǎn)稱TC)是發(fā)生在熱帶或亞熱帶洋面上的氣旋性環(huán)流,強(qiáng)烈的TC 往往伴隨著狂風(fēng)、暴雨、風(fēng)暴潮等天氣現(xiàn)象,造成嚴(yán)重的財(cái)產(chǎn)損失和人員傷亡,是氣象災(zāi)害中最具破壞性的天氣現(xiàn)象之一。我國(guó)由于幅員遼闊,所以沿海地區(qū)飽受臺(tái)風(fēng)的肆虐,是世界上受TC 影響最嚴(yán)重的國(guó)家之一。隨著沿海地區(qū)社會(huì)經(jīng)濟(jì)的不斷發(fā)展,TC 造成的人民生命安全和財(cái)產(chǎn)經(jīng)濟(jì)損失日益嚴(yán)重。例如,2017 年8 月23 日登陸廣東珠海的第13 號(hào)強(qiáng)臺(tái)風(fēng)“天鴿”造成我國(guó)6 省247.8 萬人受災(zāi),直接經(jīng)濟(jì)損失290.3 億元。因此,準(zhǔn)確地預(yù)報(bào)西北太平洋(WNP)TC 的生成、移動(dòng)路徑、強(qiáng)度變化和登陸地點(diǎn)具有重要的科學(xué)價(jià)值以及實(shí)際意義?,F(xiàn)有關(guān)于TC 的預(yù)報(bào)業(yè)務(wù)主要依賴于各種數(shù)值模式(馮文等,2016),其中TC 生成預(yù)報(bào)使用最多的是全球模式,隨著全球模式分辨率不斷提升,TC 的生成預(yù)報(bào)的準(zhǔn)確度也在不斷增強(qiáng),但是由于目前研究對(duì)于TC 生成機(jī)制認(rèn)識(shí)的局限,所以TC生成的準(zhǔn)確預(yù)報(bào)一直是數(shù)值預(yù)報(bào)中的難點(diǎn)。

    學(xué)界普遍認(rèn)為TC 生成的機(jī)制是多尺度過程共同作用的結(jié)果,天氣尺度、次天氣尺度和中尺度過程都可能有利于TC 的生成(Holland,1995)。Gray(1968,1998)提出的有利于TC 形成的大尺度背景條件包括:高的海表溫度(SST)、條件不穩(wěn)定、對(duì)流層中層的高相對(duì)濕度、對(duì)流層低層的氣旋性絕對(duì)渦度、對(duì)流層高層的反氣旋性相對(duì)渦度以及弱的垂直風(fēng)切變(Cheung,2004)。西北太平洋(WNP)是全球TC 形生成最頻繁的地區(qū),夏季W(wǎng)NP 的季風(fēng)槽區(qū)通常具備上述大部分背景條件,因此每年超過30%的TC 在WNP 生成(李憲之,1983)。TC 的生成除了具備適當(dāng)?shù)臒崃W(xué)條件,還需要有較為強(qiáng)烈的外部擾動(dòng)源。外部擾動(dòng)使對(duì)流層低層大規(guī)模高溫、高濕大氣中的不穩(wěn)定能量釋放,低層空氣輻合上升形成氣旋。TC 生成的外部擾動(dòng)源主要有:東風(fēng)波、季節(jié)內(nèi)振蕩、熱帶對(duì)流層上部槽(TUTT)、熱帶云團(tuán)、冷空氣等(張慶紅和郭春蕊,2008)。其中熱帶地區(qū)的波動(dòng)會(huì)在不同的時(shí)間尺度上調(diào)節(jié)全球的TC 活動(dòng),目前許多研究探討了熱帶波動(dòng)與TC 頻率、位置及強(qiáng)度的關(guān)系(Dickinson and Molinari, 2002; Frank and Roundy, 2006)。

    在研究TC 生成的觸發(fā)機(jī)制的早期階段,前人提出了第二類條件不穩(wěn)定(CISK)機(jī)制(Charney and Eliassen, 1964)和 風(fēng) 驅(qū) 動(dòng) 的 海 氣 熱 量 交 換(WISHE)理論(Emanuel,1986),但這些理論適用于TC 基本環(huán)流結(jié)構(gòu)建立之后的階段,其中有些TC 的暖心結(jié)構(gòu)已經(jīng)建立。在TC 生成的過程中,從無組織的低壓擾動(dòng)發(fā)展到熱帶低壓階段的時(shí)間非常短,其中具體的物理過程還不是十分明確,而且洋面上空觀測(cè)數(shù)據(jù)缺乏,所以這是臺(tái)風(fēng)領(lǐng)域研究的難點(diǎn)(Lee et al., 2008)。在TC 生成的觸發(fā)機(jī)制方面目前研究的重點(diǎn)是中小尺度過程。大量觀測(cè)資料和研究結(jié)果表明,季風(fēng)槽或熱帶波動(dòng)中的TC 形成過程中往往伴隨著中尺度對(duì)流系統(tǒng)(MCS)的發(fā)生和發(fā)展(Simpson et al., 1997; Ritchie and Holland,1997; Cheung and Elsberry, 2004)。關(guān) 于MCS 在TC 生成中的作用主要有兩種理論。第一種是“topdown”理論,即在層狀降雨區(qū)中有MCS 生成,MCS 云砧釋放的潛熱和融化層冷卻可以生成中尺度對(duì)流渦旋(MCV),雨滴的蒸發(fā)增加了低層的相對(duì)濕度,產(chǎn)生下沉運(yùn)動(dòng)并使蒸發(fā)層變低,同時(shí)下沉運(yùn)動(dòng)將MCV 向下平流,增加了海平面的風(fēng)速和渦度,增強(qiáng)了潛熱通量,導(dǎo)致對(duì)流的生成和暖心結(jié)構(gòu) 的 建 立(Bister and Emanuel, 1997)。同 時(shí)MCV 通過自身的動(dòng)力學(xué)條件配合有利的環(huán)境場(chǎng),可以激發(fā)出新的對(duì)流,使MCS 可以維持更長(zhǎng)時(shí)間(Trier and Davis, 2002)。第二種理論是“bottomup”理論,這種理論提出在低層渦度增強(qiáng)的情況下會(huì)爆發(fā)深對(duì)流(Zehr,1992),10~20 km 的渦旋熱塔(VHT)在該過程中起著重要作用(Hendricks et al., 2004; Tory et al., 2006a, 2006b)。VHT 的 存在周期很短,通過VHT 不斷地生成及消亡持續(xù)影響周圍大氣,可以使兩個(gè)或多個(gè)MCV 合并形成TC(Ritchie and Holland, 1997)。Tory and Frank(2010)指出這兩個(gè)理論并不是互相排斥的,只是關(guān)注TC 生成過程中的時(shí)間節(jié)點(diǎn)不同。“top-down”理論沒有考慮TC 生成過程中對(duì)流渦旋的加強(qiáng)作用,但他們并不否認(rèn)這一現(xiàn)象的存在,只是認(rèn)為對(duì)流渦旋的加強(qiáng)只有在TC 生成過程的后期才會(huì)出現(xiàn)?!皌op-down”和“bottom-up”理論本質(zhì)上分別是與層云和對(duì)流云的發(fā)展在TC 生成中的作用有關(guān),所以在討論這兩種理論共同作用時(shí),主要是分析兩種不同降水類型的演變規(guī)律。

    因?yàn)榕cTC 生成有關(guān)的MCS 演變主要發(fā)生在遠(yuǎn)海地區(qū),資料的缺乏使得大多數(shù)研究只能借助高分辨率數(shù)值模式。本文基于前人的已有研究,通過數(shù)值模式對(duì)2015 年第21 號(hào)超強(qiáng)臺(tái)風(fēng)“杜鵑”生成過程進(jìn)行時(shí)間和空間高分辨率的模擬,分析有利于“杜鵑”生成的天氣尺度環(huán)境背景場(chǎng);通過對(duì)MCS 的識(shí)別,重點(diǎn)研究在“杜鵑”生成過程中MCS 的演變特征,其中包括MCS 生成個(gè)數(shù),覆蓋面積,以及在MCS 中層云和對(duì)流云面積比例和強(qiáng)度變化等,揭示有利于“杜鵑”初始生成的天氣尺度系統(tǒng)和MCS 的演變特征。

    2 資料與方法

    2.1 資料來源

    本文使用的資料如下:(1)云頂亮溫?cái)?shù)據(jù)來源于葵花8 靜止衛(wèi)星(Himawari-8),空間分辨率為0.05°;(2)近表面降水率與云類數(shù)據(jù)來自GPM(Global Precipitation Measurement)極軌衛(wèi)星,空間分辨率為5 km;(3)“杜鵑”生成位置(熱帶低壓時(shí)刻)數(shù)據(jù)來源于美國(guó)關(guān)島聯(lián)合臺(tái)風(fēng)警報(bào)中心(JTWC)的最佳路徑集資料;(4)風(fēng)場(chǎng)和模式背景場(chǎng)數(shù)據(jù)來源于美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NECP)/美國(guó)國(guó)家大氣研究中心(NCAR)6 h一次的FNL 全球分析資料,空間分辨率為1°。

    2.2 模式介紹

    為了研究“杜鵑”生成過程中MCS 的演變特征,本文利用WRF-ARW3.6.1 模式對(duì)“杜鵑”的生成過程進(jìn)行高時(shí)空分辨率模擬。模擬選擇WSM6 的微物理方案(Hong and Lim, 2006),積云參數(shù)化方案為Tiedtke 方案(Zhang et al., 2011),YSU 行星邊界層方案(Shin et al., 2012),RRTMG長(zhǎng) 波 輻 射 方 案 和 短 波 輻 射 方 案(Iacono et al.,2008)。模 擬 采 用 三 重 嵌 套 網(wǎng) 格(27/9/3 km),格點(diǎn)數(shù)分別為340×340、784×940、700×901(圖1)。垂直方向σ 層個(gè)數(shù)為45 層,模式層頂定為50 hPa。10 km 以下的網(wǎng)格區(qū)域沒有使用積云參數(shù)化方案。模擬時(shí)間從2015 年9 月16 日18:00(協(xié)調(diào)世界時(shí),下同),比JTWC 最佳路徑集資料記錄的熱帶低壓生成時(shí)間(20 日18:00)提前了96 h,至9 月20 日18:00 結(jié)束。模擬的初始條件和邊界條件均來源于FNL 資料。

    3 結(jié)果分析

    3.1 “杜鵑”概況與生成期間背景場(chǎng)分析

    圖1 WRF 模式三重嵌套的模擬區(qū)域。黑色實(shí)線代表“杜鵑”移動(dòng)路徑(JTWC,2015 年9 月20 日18:00 至29 日18:00),紅色實(shí)線表示“杜鵑”生成前的移動(dòng)路徑(2015 年9 月17 日18:00 至20 日18:00),圓點(diǎn)代表每日的00:00。圖中00Z18 表示9 月18日00:00Fig.1 The triple nested grid domain configuration of the WRF model.The black solid line is the track of Dujuan (JTWC, 1800 UTC 17 to 1800 UTC September 29, 2015), red solid line is the track before Dujuan formation (1800 UTC 17 to 1800 UTC September 20, 2015).The dots represent 0000 UTC of each day 0000 UTC, 00Z18 means 0000 UTC on September 18, 2015

    2015 年第21 號(hào)超強(qiáng)臺(tái)風(fēng)“杜鵑”于2015 年9 月20 日18:00 在西北太平洋洋面上(15.1°N,145.3°E)達(dá)到季風(fēng)低壓(MD)強(qiáng)度[最大風(fēng)速為20 kt(1 kt=1.852 km h?1)],并被JTWC 正式編號(hào)。隨后“杜鵑”繼續(xù)向西北方向移動(dòng)(圖1 黑色實(shí)線),于23 日12:00 增強(qiáng)為熱帶風(fēng)暴(TS,最大風(fēng)速為45 kt),9 月25 日00:00 達(dá)到臺(tái)風(fēng)(TY)強(qiáng)度(最大風(fēng)速為65 kt),27 日18:00 升級(jí)為超強(qiáng)臺(tái)風(fēng)(ST,最大風(fēng)速為130 kt),登陸后強(qiáng)度減弱,JTWC 于29 日18:00 停止對(duì)其記錄(最大風(fēng)速為30 kt)?!岸霹N”于9 月28 日登陸日本八重山群島,最大瞬時(shí)風(fēng)速達(dá)到81.1 m s?1,突破1994年臺(tái)風(fēng)“德格”70.2 m s?1的記錄。登陸臺(tái)灣期間,受其影響,全臺(tái)灣超過208 萬戶停電,個(gè)別地區(qū)發(fā)生洪澇災(zāi)害,風(fēng)暴直接造成2 人死亡,324 人受傷。

    由于JTWC 最佳路徑集資料大多是從熱帶氣旋強(qiáng)度達(dá)到熱帶低壓時(shí)刻開始記錄的,而達(dá)到熱帶低壓強(qiáng)度的氣旋往往已經(jīng)具有一定的組織化結(jié)構(gòu),所以本文從“杜鵑”生成時(shí)刻(2015 年9 月20 日18:00)向前追蹤三天至9 月17 日18:00,來探究有利于“杜鵑”生成的大尺度環(huán)境條件。向前追蹤的方法為:將850 hPa 的風(fēng)場(chǎng)進(jìn)行3~8 天的濾波(Peng et al., 2012;Fu et al., 2012),求取相應(yīng)的流函數(shù),定義流函數(shù)的最低點(diǎn)為擾動(dòng)中心。

    有利于“杜鵑”生成的大尺度背景場(chǎng)主要與熱帶對(duì)流層上部槽和季風(fēng)渦旋有關(guān)。熱帶對(duì)流層上部槽是暖季形成于北太平洋中部及北大西洋中部熱帶地區(qū)對(duì)流層上部的低壓槽,一般在200 hPa 最為明顯,又稱洋中槽,簡(jiǎn)稱TUTT。國(guó)內(nèi)外的一些研究發(fā)現(xiàn),TUTT 附近的中尺度對(duì)流活動(dòng)會(huì)通過加熱加濕大氣來促進(jìn)TC 生成(Sadler,1978;Kelley and Mock, 1982)。另一方面,TUTT 南部的強(qiáng)西風(fēng)切變及其他與TUTT 東側(cè)相關(guān)的環(huán)境條件通常不利于熱帶云團(tuán)的發(fā)展和TC 的形成,通常使WNP 上TC 的生成位置無法向東延伸(Gray,1968)。TUTT的東進(jìn)和西退主要影響在WNP 東部生成的TC,西伸的TUTT 會(huì)抑制在WNP 東側(cè)的TC 生成活動(dòng)(Wang and Wu, 2016)。

    WNP 夏季季風(fēng)環(huán)流在對(duì)流層低層通常表現(xiàn)為季風(fēng)槽結(jié)構(gòu)(Holland,1995),但有時(shí)季風(fēng)槽會(huì)被一個(gè)直徑約為2500 km 的大型季風(fēng)渦旋(MG)取 代(Lander,1994),MG 會(huì) 對(duì)WNP 上TC 的生成位置和移動(dòng)路徑產(chǎn)生重要的影響(Ritchie and Holland, 1999)。據(jù)Wu et al.(2013)的統(tǒng)計(jì)結(jié)果,在2000~2005 年5~10 月期間共出現(xiàn)37 個(gè)季風(fēng)渦旋,其中31 個(gè)MG 伴隨著42 個(gè)TC 的生成,占熱帶氣旋生成總數(shù)的19.8%。10 天以上850 hPa 低通濾波風(fēng)場(chǎng)可以較好地顯示MG 的結(jié)構(gòu)特征(Wu et al.,2013),圖4 顯示MG 的位置約在10~20°N,130~150°E,“杜鵑”在17 日18:00(約6°N,158°E)向西或西北偏西方向移動(dòng),到20 日18:00(約14°N,145°E)發(fā)展成為熱帶低壓。以往的研究發(fā)現(xiàn),MG 環(huán)流東南側(cè)存在強(qiáng)烈的氣流輻合區(qū),有利于東部西傳的天氣尺度波動(dòng)加強(qiáng)發(fā)展成TC(Webster and Chang, 1988;Holland,1995)。在“杜鵑”生成的前72 h(圖4a),擾動(dòng)中心位于季風(fēng)渦旋東南側(cè)的渦度大值區(qū),越赤道氣流的西南風(fēng)和偏東風(fēng)在MG 的東部邊緣匯合,形成強(qiáng)的東西輻合區(qū),為TC 的發(fā)展提供有利環(huán)境;在“杜鵑”生成前的48 h(圖4b),擾動(dòng)中心向西北移動(dòng),由于此時(shí)擾動(dòng)中心正好位于東西風(fēng)輻合區(qū)的交匯處,“杜鵑”中心位置附近的渦度有較為明顯的增強(qiáng);在“杜鵑”生成前的24 h(圖4c),擾動(dòng)中心位置附近仍為渦度大值區(qū),但是渦度值小于前48 h(圖4b);在“杜鵑”生成時(shí)刻(圖4d),渦旋中心移動(dòng)到MG 中心的東側(cè)。

    圖2 2015 年9 月(a)19 日18:00、(b)20 日18:00 FNL 在分析資料的200 hPa 無濾波風(fēng)場(chǎng)(矢量,單位:m s?1)和500 hPa 位勢(shì)高度場(chǎng)(陰影,單位:gpm),黑色圓點(diǎn)表示該時(shí)刻“杜鵑”處于擾動(dòng)狀態(tài)的中心位置,粗黑線代表洋中槽TUTT 槽線Fig.2 The unfiltered 200-hPa wind field (vector, units: m s?1) and 500-hPa potential height field (shaded, units: gpm) at (a) 1800 UTC 19 September and (b) 1800 UTC 20 September, 2015 from FNL (Final Operational Global Analysis) data. Black dots are the central positions of Dujuan in the stage of tropical disturbance. Black lines represent the trough lines of TUTT (tropical upper tropospheric trough)

    3.2 模式結(jié)果驗(yàn)證

    在TUTT 和MG 有利的大尺度環(huán)境場(chǎng)共同作用下,“杜鵑”逐漸從弱的擾動(dòng)發(fā)展成有組織的熱帶低壓。較粗時(shí)空分辨率的分析資料無法更好地體現(xiàn)“杜鵑”組織化的過程,本文將繼續(xù)用數(shù)值模式來探究“杜鵑”生成過程中MCS 的演變特征。為了驗(yàn)證模式結(jié)果的可靠性,本節(jié)將模式最外層網(wǎng)格的模擬結(jié)果插值成1°×1°后,利用3~8 天850 hPa風(fēng)場(chǎng)濾波的追蹤方法確定擾動(dòng)中心,與用FNL 資料追蹤的結(jié)果進(jìn)行對(duì)比驗(yàn)證(圖5)。由圖5 可見,模擬結(jié)果與FNL 資料中的“杜鵑”路徑從17 日18:00 至20 日18:00 均是向西北方向移動(dòng),趨勢(shì)一致,局部的路徑重合,說明模式模擬的路徑結(jié)果較好。

    除了對(duì)比路徑結(jié)果,本文再選擇三個(gè)時(shí)次的觀測(cè)資料與模式模擬結(jié)果數(shù)據(jù)進(jìn)行對(duì)比。除了葵花8衛(wèi)星的亮溫產(chǎn)品,我們還選擇GPM 衛(wèi)星的近表面降水率產(chǎn)品作為參考。在這三個(gè)時(shí)次中,GPM 衛(wèi)星的掃描帶均位于距離擾動(dòng)中心半徑500 km 的圓域范圍內(nèi)。對(duì)比觀測(cè)資料(圖6a、b、c)和模式結(jié)果(圖6d、e、f)可以發(fā)現(xiàn):圖6a 中GPM 觀測(cè)的降水大值區(qū)位于(4°N~6°N,152°E~154°E),圖6d 中模擬的降水在相應(yīng)的區(qū)域也模擬出相同量級(jí)的降水。雖然圖6d 中擾動(dòng)中心偏北的云區(qū)覆蓋范圍比觀測(cè)要小,但是模擬的云頂亮溫很好地反映出了實(shí)際擾動(dòng)中的對(duì)流特征。從圖6b 與圖6e 的對(duì)比可見,擾動(dòng)中心附近的降水大值區(qū)位置以及低亮溫的區(qū)域都有很好的對(duì)應(yīng)。圖6c 中GPM 軌道掃到的降水格點(diǎn)較少,大范圍云區(qū)位于擾動(dòng)西側(cè),圖6f 中大范圍的云區(qū)也同樣位于擾動(dòng)中心西側(cè),只是位置較圖6c 略偏東。模式結(jié)果較好地反映了實(shí)際的對(duì)流活動(dòng),在余下各小節(jié),本文將利用模式結(jié)果對(duì)“杜鵑”整個(gè)環(huán)流(距離擾動(dòng)中心半徑500 km的圓域范圍)內(nèi)的MCS 以及MCS 中的層云降水和對(duì)流降水演變特征進(jìn)行研究。

    圖3 2015 年9 月20 日12:00(a)200 hPa、(b)850 hPa 的FNL 無濾波風(fēng)場(chǎng)(矢量,單位:m s?1)和垂直風(fēng)速切變(陰影,單位:m s?1),黑色圓點(diǎn)表示“杜鵑”的中心位置Fig.3 The unfiltered wind fields (vector, units: m s?1) at 1200 UTC 20 September on (a) 200 hPa and (b) 850 hPa from FNL data, which are overlaid by the vertical wind shear (shaded, units: m s?1) between 200-hPa and 850-hPa levels. Black dots are the center locations of Dujuan

    3.3 MCS 的演變特征

    “杜鵑”從擾動(dòng)階段發(fā)展成為有組織的熱帶低壓必定伴隨著MCS 的生成、合并和發(fā)展過程。Houze(1993)將MCS 定義為至少在某一個(gè)方向能產(chǎn)生100 km 持續(xù)降水的積雨云系統(tǒng),這個(gè)定義聯(lián)系了降水的范圍大小和強(qiáng)度。在衛(wèi)星圖像中定義MCS 有很多種方法:利用日本靜止氣象衛(wèi)星(GMS)數(shù)據(jù),Ritchie and Holland(1999)將云頂亮溫小于214 K(?56.15°C),面積大于4×104km2,偏心率大于0.5 的深對(duì)流系統(tǒng)定義為MCS;利用偏振修正亮溫(PCT,利用TRMM 85 GHz 微波產(chǎn)品計(jì)算得到)數(shù)據(jù),Nesbitt et al.(2000)對(duì)PCT≤250 K 的衛(wèi)星格點(diǎn)數(shù)據(jù)的連續(xù)性進(jìn)行判斷,如果其連續(xù)面積大于等于2000 km2,并且其中PCT≤225 K的連續(xù)面積大于等于185 km2,則這個(gè)降水系統(tǒng)被定義為MCS?;谇叭说难芯浚疚睦靡韵屡袚?jù)在模式中來定義MCS:在距離擾動(dòng)中心500 km范圍內(nèi)(1)格點(diǎn)降水率大于等于1.5 mm h?1;(2)連續(xù)降水格點(diǎn)面積大于2000 km2;(3)對(duì)流降水(降雨率大于或等于4.5 mm h?1)面積大于185 km2;(4)任意方向上有100 km 以上(MCS 長(zhǎng)度)的連續(xù)降水。

    圖4 2015 年9 月(a)17 日18:00、(b)18 日18:00、(c)19 日18:00 和(d)20 日18:00 的850 hPa 10 天以上低通濾波風(fēng)場(chǎng)(矢量,單位:m s?1)和850 hPa 渦度場(chǎng)(陰影,單位:10?5 s?1),黑色圓點(diǎn)為“杜鵑”生成前的位置Fig.4 The 10-day-low-pass filtered 850-hPa wind (vectors, units: m s?1) and 850-hPa vorticity fields (shaded, units: 10?5 s?1) at (a) 1800 UTC 17,(b) 1800 UTC 18, (c) 1800 UTC 19, and (d) 1800 UTC 20 September 2015. Black dots are the center positions before Dujuan formation

    圖5 模式模擬結(jié)果(藍(lán)線)和FNL 分析資料追蹤(紅線)的“杜鵑”路徑,較大的圓點(diǎn)表示“杜鵑”中心每日00:00 所在位置Fig.5 The tracks of Dujuan traced by model results (blue solid line)and FNL data (red solid line). The larger dots indicate the center positions of Dujuan at 0000 UTC of each day

    根據(jù)Lee et al.(2008)提出的方法,利用生成前48 小時(shí)的925 hPa 和850 hPa 風(fēng)場(chǎng)數(shù)據(jù),以擾動(dòng)中心為中心、將5°×5°的區(qū)域分為四個(gè)象限。隨后判斷每個(gè)象限中平均緯向風(fēng)和經(jīng)向風(fēng)的大小及方向,風(fēng)速>5 m s?1定義為強(qiáng)風(fēng),<2 m s?1定義為弱風(fēng),并根據(jù)這個(gè)標(biāo)準(zhǔn)將影響TC 生成的天氣形勢(shì)分為六種:東風(fēng)波型(EW)、東北風(fēng)主導(dǎo)型(NE)、東北風(fēng)與西南風(fēng)共存型(NE-SW)、西南季風(fēng)主導(dǎo)型(SW)、季風(fēng)輻合型(MC)和季風(fēng)切變型(MS)。按照此標(biāo)準(zhǔn)計(jì)算“杜鵑”生成前的風(fēng)場(chǎng)數(shù)據(jù),得到東北象限緯向風(fēng)為?8.0 m s?1,西南象限緯向風(fēng)為9.5 m s?1,確定影響“杜鵑”生成的天氣形勢(shì)為MS 型。據(jù)Lee et al.(2008)的統(tǒng)計(jì)結(jié)果,MS 型中的對(duì)流活動(dòng)傾向于在擾動(dòng)中心的西側(cè)形成,這與“杜鵑”中MCS 的水平空間分布一致,面積最大的MCS 通常位于擾動(dòng)中心的西側(cè)附近(圖7),遠(yuǎn)離擾動(dòng)中心的MCS 面積較小且分布零散。在“杜鵑”的生成過程中,MCS 每小時(shí)的水平空間分布都有較大的不同,說明MCS 的生成、合并、消亡周期較短。MCS 在擾動(dòng)的各個(gè)象限均有生成,消亡的MCS 通常面積較小。合并過程多發(fā)生在位于擾動(dòng)中心西側(cè)最大的MCS 附近,并且MCS 的合并通常伴隨著面積較大的MCS 從內(nèi)核區(qū)域的外側(cè)向內(nèi)側(cè)收縮的過程。同時(shí)MS 型的天氣形勢(shì)更有利于熱帶擾動(dòng)在發(fā)展過程中生成多個(gè)MCS 以及多個(gè)MCS 同時(shí)存在,MCS 內(nèi)部的層云降水區(qū)容易滿足MCV 生成的熱力學(xué)條件,如果MCS 中生成兩個(gè)或多個(gè)MCV,并且這些MCV 的距離足夠近,通過MCV 之間的相互作用或合并,可以形成更強(qiáng)的渦旋(Lee et al., 2008)。

    圖6 2015 年9 月(a)18 日09:00、(b)18 日21:00 和(c)20 日09:00 用FNL 資料追蹤的“杜鵑”中心位置(紅色圓點(diǎn)),葵花8 衛(wèi)星的亮溫產(chǎn)品(黑白陰影,單位:°C)和GPM 衛(wèi)星的近表面降水率產(chǎn)品(彩色陰影,單位:mm h?1)。(d?f)與(a?c)相同,但是為數(shù)值模擬結(jié)果。白色實(shí)線表示GPM 衛(wèi)星的軌道范圍。紅色實(shí)線為距離“杜鵑”中心半徑500 km 的范圍Fig.6 The center positions (red dots) of Dujuan traced by FNL data, the infrared brightness temperature (shaded by black and white, units: °C) from Himawari-8 satellite, and the near-surface precipitation rate (colorfully shaded, units: mm h?1) from GPM satellite at (a) 0900 UTC 18, (b) 2100 UTC 18, and (c) 0900 UTC 20 September 2015. (d), (e) and (f) are same as (a), (b) and (c) respectively, but for the numerical model results. White solid lines are the orbital boundaries of GPM swath, and the red circles indicate the radius of 500 km from Dujuan center

    除了MCS 的水平空間分布特征,本節(jié)利用MCS 中的相對(duì)渦度和動(dòng)能繼續(xù)分析MCS 垂直結(jié)構(gòu)隨時(shí)間的演變(圖8)。正相對(duì)渦度的大值區(qū)一般分布在對(duì)流層中低層,高層通常為負(fù)渦度。18 日12:00 至18 日18:00,即“杜鵑”生成的早期階段,渦度有從對(duì)流層中層向下發(fā)展的特征,逐漸增強(qiáng)的低層渦度可以增加潛熱通量,促使擾動(dòng)中對(duì)流活動(dòng)的爆發(fā)和暖心結(jié)構(gòu)形成(Bister and Emanuel, 1997)。19 日06:00~12:00,900~600 hPa 的渦度有較為明顯的增強(qiáng),隨后渦度振蕩增強(qiáng),大值區(qū)出現(xiàn)在19日15:00 至20 日03:00 以及“杜鵑”生成前2 h 至生成時(shí)刻(20 日16:00~18:00)的對(duì)流層中低層附近。

    17 日18:00 至19 日04:00 1000~500 hPa 的MCS 動(dòng)能較小,沒有大值區(qū)域分布。19 日08:00對(duì)流層低層(900 hPa 附近)出現(xiàn)了動(dòng)能大值中心,與渦度場(chǎng)的大值區(qū)域出現(xiàn)的時(shí)間和高度匹配,中低層渦度增強(qiáng),此時(shí)對(duì)流層高層的動(dòng)能依舊較小。在20 日04:00~08:00 和16:00~18:00,MCS 低層均出現(xiàn)動(dòng)能的極大值區(qū)域,高層的動(dòng)能也有較明顯的增加,與MCS 渦度場(chǎng)的配合較好,低層渦度增強(qiáng)伴隨著動(dòng)能增加的過程,有利于深對(duì)流的爆發(fā),導(dǎo)致對(duì)流云降雨率增加,本文將在3.4 節(jié)討論MCS 內(nèi)低層渦度增強(qiáng)、動(dòng)能增加的時(shí)刻與對(duì)流云降雨率變化的對(duì)應(yīng)情況。

    圖7 2015 年9 月(a?p)18 日10:00 至19 日01:00 距離“杜鵑”中心500 km 范圍內(nèi)的MCS 的空間分布。每個(gè)時(shí)次中不同的MCS 用不同顏色表示,但不同時(shí)次相同顏色的MCS 并不一定代表同一個(gè)MCS。黑色圓點(diǎn)為“杜鵑”中心位置Fig.7 The spatial distribution of MCSs (Mesoscale Convective System) from (a?p) 1000 UTC 18 to 0100 UTC 19 September 2015 within 500 km of Dujuan center. Different MCS are represented by different color in each figure. But MCSs with a same color in different figures don't necessarily point to a same MCS. Black dot indicates the center location of Dujuan

    本文利用10 m 高度方位角平均切向風(fēng)最大值(Vmax)來代表“杜鵑”的強(qiáng)度(圖9a),并對(duì)每小時(shí)MCS 的特征進(jìn)行統(tǒng)計(jì),其中包括MCS 的總面積、MCS 的個(gè)數(shù)和MCS 的平均降雨率(圖9b?d)。最大方位角平均切向風(fēng)速顯示,“杜鵑”的強(qiáng)度在9 月17 日18:00 至19 日00:00 增加幅度較為平緩,于19 日00:00~19 日15:00 迅速增加,之后Vmax呈振蕩增加的趨勢(shì)。雖然20 日18:00 強(qiáng)度達(dá)到了16.6 m s?1,強(qiáng)度偏強(qiáng),但是結(jié)合衛(wèi)星觀測(cè)的亮溫、降水以及路徑,模擬還是很好地再現(xiàn)了“杜鵑”的生成過程。

    圖8 2015 年9 月17 日18:00 至20 日18:00 所有MCS 區(qū)域平均的(a)相對(duì)渦度(單位:10?5 s?1)和(b)動(dòng)能(單位:J)的高度—時(shí)間演變Fig.8 The areal averaged height–time evolution of (a) relative vorticity (units: 10?5 s?1) and (b) kinetic energy (units: J) of all MCSs from 1800 UTC 17 to 1800 UTC 20 September 2015

    圖9b 表示MCS 的總面積隨時(shí)間的演變,MCS 的總面積于9 月18 日06:00 至19 日00:00 增加,19 日00:00~12:00 減小,于19 日00:00 左右達(dá)到最大值。對(duì)比圖9a 和b 可以發(fā)現(xiàn),在19 日00:00~12:00,MCS 面積的減小伴隨著Vmax的迅速增加。由于MCS 面積減小的過程通常伴隨著MCS 的合并,以及原先較為分散的MCS 變得更加緊湊并且向擾動(dòng)中心收縮(圖7),因此MCS 面積的收縮有利于熱帶擾動(dòng)的增強(qiáng)。

    逐小時(shí)MCS 的個(gè)數(shù)隨時(shí)間存在著顯著的震蕩(圖9c),說明MCS 的發(fā)生、合并和消亡周期較短。此外MCS 日個(gè)數(shù)分別為109(第一天:17 日19:00 至18 日18:00)、124(第二天:18 日19:00至19 日18:00)和127(第三天:19 日19:00 至20 日18:00),第一天至第二天MCS 個(gè)數(shù)的增長(zhǎng)幅度明顯大于第二天至第三天的增長(zhǎng)幅度,并且從第一天至第二天有一個(gè)明顯的渦度向下發(fā)展過程(圖8),這些過程均發(fā)生在“杜鵑”生成的早期階段。根據(jù)“top-down”理論,在TC 的生成過程中,層云降雨區(qū)的熱力條件可以在MCS 中生成MCV,MCV 通過下沉運(yùn)動(dòng)向下平流,增強(qiáng)了對(duì)流層中低層的渦度,同時(shí)通過MCV 自身的動(dòng)力學(xué)條件配合有利的環(huán)境場(chǎng),可以激發(fā)出新的對(duì)流,這些對(duì)流又可以組織發(fā)展形成新的MCS(Trier and Davis, 2002)。圖9a 中擾動(dòng)的強(qiáng)度在第一天有小幅度的增加,第二天的增加趨勢(shì)更為明顯,說明在“杜鵑”生成的早期階段,MCS 個(gè)數(shù)的增加,即新生成更多的MCS 有利于“杜鵑”的初始發(fā)展。

    圖9 2015 年9 月17 日18:00 至20 日18:00 逐小時(shí)(a)“杜鵑”中心附近10 m 高度上方位角平均切向風(fēng)最大值(Vmax,單位:m s?1)、(b)MCS 總面積(單位:km2)、(c)MCS 個(gè)數(shù)和(d)MCS 的平均降雨率(單位:mm h?1)的演變Fig.9 The temporal evolution of (a) the maximum of azimuthal averaged tangential wind speed at 10 m height (Vmax, units: m s?1) near the center of Dujuan, (b) the total areal coverage of MCSs (units: km2),(c) the total number of MCSs, and (d) the averaged precipitation rate of all MCSs (units: mm h?1) from 1800 UTC 17 to 1800 UTC 20 September 2015

    本文使用MCS 的平均降雨率(圖9d)來代表MCS 的強(qiáng)度,MCS 強(qiáng)度隨時(shí)間的變化(圖9d)與面積隨時(shí)間的變化(圖9b)在“杜鵑”生成中后期大致呈反位相關(guān)系。MCS 的總面積在19 日00:00 左右達(dá)到最大值,而平均降雨率為最小值;MCS 的降雨率在19 日20:00 左右達(dá)到最強(qiáng)時(shí),MCS 的總面積卻為局部極小值;20 日14:00 左右,MCS 的強(qiáng)度與面積也呈相反的變化趨勢(shì)。聯(lián)系圖9c 和d,MCS 的平均降雨率與個(gè)數(shù)同樣呈反位相的演變特征,例如:MCS 的個(gè)數(shù)在“杜鵑”生成的后期階段(19 日15:00)達(dá)到局部極大值,之后MCS 的個(gè)數(shù)減小,平均降雨率開始增加,在19 日19:00 平均降雨率達(dá)到最大值,此時(shí)MCS 僅有2 個(gè),這與MCS 自身的消亡和多個(gè)MCS 之間的合并有關(guān)。在“杜鵑”生成的前期,擾動(dòng)內(nèi)核區(qū)有很多MCS 生成,這些MCS 會(huì)在較大的范圍內(nèi)形成中低壓,而中低壓則有利于氣旋性風(fēng)速的增加和海平面氣壓的減小。在“杜鵑”生成后期,強(qiáng)的對(duì)流爆發(fā),在對(duì)流層低層產(chǎn)生大的渦度并向上輸送,同時(shí)中低層動(dòng)能都有較明顯的增加過程(圖8),這一過程伴隨著海表熱通量的增加、低層輻合和氣旋性渦度增加等現(xiàn)象,有利于MCS 的合并加強(qiáng)(Kieu and Zhang, 2008)。

    3.4 MCS 中層云和對(duì)流云的演變特征

    熱帶大氣中的降水主要分為兩類:層云降水和對(duì)流云降水。這兩種類型的降水在動(dòng)力和熱力特征上具有明顯的不同:層云降水的特征是高低層輻散、中層輻合;對(duì)流云降水則為低層輻合、高層輻散。相應(yīng)的在層云降水中,融化層以下為下沉氣流和非絕熱冷卻,以上為上升氣流和非絕熱增溫;在對(duì)流云降水中,強(qiáng)的上升氣流和水汽的凝結(jié)加熱可以分布于整個(gè)對(duì)流層低層(Houze,1997)。

    因?yàn)閷釉平邓蛯?duì)流云降水的動(dòng)力和熱力特征差異十分明顯,本文利用兩類降水非絕熱加熱廓線的特征作為判據(jù)來區(qū)分兩類降水:2 km、3 km、4 km 的非絕熱加熱率小于0,6 km 的非絕熱加熱率大于0 的格點(diǎn)定義為層云降水;2 km、3 km、4 km、6 km 的非絕熱加熱率均大于0 的格點(diǎn)定義為對(duì)流云降水。

    為了驗(yàn)證使用上述分類方法得到的結(jié)果與GPM 衛(wèi)星觀測(cè)數(shù)據(jù)是否吻合,本文將分類結(jié)果與圖6a,b 中降水格點(diǎn)數(shù)較多的軌道進(jìn)行對(duì)比。圖10為GPM 衛(wèi)星產(chǎn)品與依靠WRF 模擬結(jié)果分類得到的層云降水和對(duì)流云降水占總降水格點(diǎn)的百分比。衛(wèi)星產(chǎn)品中層云降水所占比例分別為77.1%(軌道號(hào)8840)和80.7%(軌道號(hào)8833),略低于模式對(duì)應(yīng)時(shí)刻的結(jié)果83.6%和84.8%;而GPM 產(chǎn)品中對(duì)流云降水比例則略高于模式結(jié)果。Fritz et al.(2016)將TC 生成前降水分為四種類型:層云降水、淺對(duì)流降水、中層對(duì)流降水和深對(duì)流降水,它們占總降水格點(diǎn)的百分比在TC 生成前3 天到生成后1 天變化不大,大約分別為80%、7%、10%和3%,本模擬所得到的層云和對(duì)流云降水比例與前人研究基本吻合,表明本研究所采用的分類方法能夠合理的再現(xiàn)層云和對(duì)流云的分布特征。在所有降水格點(diǎn)中,層云降水所占總降水格點(diǎn)的比例遠(yuǎn)大于對(duì)流云降水所占的比例,高的層云比例可以代表擾動(dòng)較好的組織化程度。

    圖10 距離“杜鵑”中心500 km 范圍內(nèi),對(duì)流降水和層云降水格點(diǎn)分別占這兩類降水總格點(diǎn)數(shù)的百分比。8833 和8840 代表GPM衛(wèi)星的軌道號(hào)(圖6a 和b),Con 和Str 分別表示對(duì)流降水和層云降水。紅色柱狀圖為GPM 2ADPR 產(chǎn)品,藍(lán)色柱狀圖為利用WRF模擬結(jié)果在相應(yīng)時(shí)刻GPM 軌道范圍內(nèi)的分類結(jié)果Fig.10 The percentage of the grid number of convective/stratiform precipitation to the total grid number of these two types of precipitation within 500 km of Dujuan center. 8833 and 8840 indicate the GPM swath number as shown in Fig.6a and Fig.6b, respectively. Con and Str indicate convective and stratiform precipitation, respectively. The red histograms represent the statistics from the GPM 2ADPR product,and the blue histograms are the classification results from WRF outputs within the same areal coverage of GPM swath at the same time

    為了研究不同降水類型對(duì)“杜鵑”生成的貢獻(xiàn),圖11 為所有MCS 內(nèi)層云和對(duì)流云降水所占比例及降雨率的逐小時(shí)變化。9 月17 日18:00 至19 日15:00 層云(藍(lán)色實(shí)線)和對(duì)流云降水(藍(lán)色虛線)所占比例保持相對(duì)穩(wěn)定,之后呈振蕩變化。通過圖11 和圖9a 可以看出,19 日15:00~18:00 層云降水所占比例減小,對(duì)流云降水所占比例增加,同時(shí)Vmax也減小;19 日18:00 至20 日03:00 層云降水所占比例先增加再減小,與Vmax呈同位相變化,對(duì)流云降水與層云降水所占比例的變化相反,表明層云降水百分比值的增加在“杜鵑”的增強(qiáng)過程中扮演更為重要的角色。通過圖11 和圖9b 可以發(fā)現(xiàn),19 日15:00 之后,層云降水所占比例隨時(shí)間的變化與MCS 面積大小的演變規(guī)律相同,這同時(shí)也表明在“杜鵑”增強(qiáng)至熱帶低壓時(shí)刻的前一天,MCS的面積大小主要取決于層云面積的大小。

    圖11 所有MCS 中層云降水(藍(lán)色實(shí)線)和對(duì)流降水(紅色實(shí)線)的格點(diǎn)數(shù)占這兩類降水總格點(diǎn)數(shù)的百分比,以及層云降水(藍(lán)色虛線)和對(duì)流降水(紅色虛線)的降雨率(單位:mm h?1)隨時(shí)間的演變。Fig.11 The temporal evolution of the percentage of the grid number of stratiform (blue solid line) and convective (red solid line) precipitation to the total grid number of the precipitating pixels of these two precipitation types and the precipitation rate (units: mm h?1) of stratiform (blue dashed line)and convective (red dashed line) precipitation in all MCS

    本文利用降雨率來代表層云降水(紅色實(shí)線)和對(duì)流降水(紅色虛線)的強(qiáng)度(圖11),層云降雨率的變化幅度比對(duì)流云降雨率的變化幅度更小。通過聯(lián)系兩類降水的強(qiáng)度(圖11)和MCS 的強(qiáng)度演變特征(圖9d)可以發(fā)現(xiàn),對(duì)流云降雨率的變化與MCS 強(qiáng)度的變化一致。在19 日06:00~12:00“杜鵑”增強(qiáng)的過程中,低層的渦度和動(dòng)能均為先增強(qiáng)再減弱(圖8),與對(duì)流降雨率在同一時(shí)間段的變化一致。Zehr(1992)指出低層渦度增強(qiáng)誘發(fā)深對(duì)流活動(dòng),深對(duì)流通過CISK 和WISHE機(jī)制正反饋?zhàn)饔眠M(jìn)一步加強(qiáng)渦旋(Montgomery et al., 2006);占較大百分比的層云則可以不斷加濕大 氣,促 進(jìn)MCV 的 合 并(Bister and Emanuel,1997)。因此層云降水的覆蓋面積決定了MCS 的面積,對(duì)流降水的降水強(qiáng)度決定了MCS 的強(qiáng)度,層云與對(duì)流云的共同作用促進(jìn)了MCS 的發(fā)展和“杜鵑”的生成。

    4 結(jié)論

    TC 生成初期是從熱帶擾動(dòng)向熱帶低壓轉(zhuǎn)變的過程,這一過程伴隨著MCS 的生成、組織化、合并和發(fā)展過程,目前對(duì)于在TC 生成過程中MCS的演變特征及其所扮演的角色還不是很了解。本文通過FNL 分析資料和高分辨率的數(shù)值模式從“杜鵑”生成時(shí)刻前3 天開始追蹤和模擬,探討了“杜鵑”生成過程中MCS 的發(fā)展演變及MCS 中的層云降水和對(duì)流降水特征,結(jié)果表明:

    (1)在臺(tái)風(fēng)“杜鵑”生成的大尺度環(huán)境中,TUTT 較平均位置偏東,減弱的垂直風(fēng)切有利于“杜鵑”未來的生成發(fā)展,同時(shí)“杜鵑”生成過程中的移動(dòng)路徑是從季風(fēng)渦旋的東南側(cè)邊緣向西北方向移動(dòng),季風(fēng)渦旋的東南側(cè)的輻合區(qū)有助于“杜鵑”渦度的增大。

    (2)在“杜鵑”生成過程中,影響“杜鵑”生成的天氣形勢(shì)表現(xiàn)為季風(fēng)切變型(MS),面積最大的MCS 通常位于擾動(dòng)中心的西側(cè),遠(yuǎn)離擾動(dòng)中心的MCS 面積較小且分布零散。在“杜鵑”生成的前期,新生成更多的MCS 有利于“杜鵑”的增強(qiáng),同時(shí)MCS 中層的相對(duì)渦度有一個(gè)明顯的向下發(fā)展過程,說明在這個(gè)階段“杜鵑”經(jīng)歷了“top-down”發(fā)展;在“杜鵑”生成的中后期階段,對(duì)流爆發(fā),MCS 發(fā)生合并導(dǎo)致MCS 個(gè)數(shù)減少,并且在合并過程中,MCS 變得更加緊湊,面積最大的MCS 逐漸向擾動(dòng)中心收縮。在這一階段低層渦度和動(dòng)能有突然增強(qiáng)的現(xiàn)象,表明擾動(dòng)中有深對(duì)流爆發(fā),深對(duì)流通過CISK 和WISHE 機(jī)制的正反饋?zhàn)饔瞄g歇發(fā)展,符合“bottom-up”理論提出的觀點(diǎn)。

    (3)通過研究不同類型的降水在“杜鵑”生成過程中的貢獻(xiàn)作用可以發(fā)現(xiàn),MCS 中層云降水的覆蓋面積較大,是對(duì)流云降水的4 倍以上,更能體現(xiàn)擾動(dòng)和降水的組織化程度;對(duì)流云降水的降雨率比層云的降雨率更大,小時(shí)變化幅度更明顯。層云降水占總降水格點(diǎn)的百分比增加與“杜鵑”增強(qiáng)的過程關(guān)系密切,對(duì)流降水降雨率的增加有利于MCS 的增強(qiáng),層云降水和對(duì)流降水的共同作用促進(jìn)了“杜鵑”的生成。

    本文雖然探討了臺(tái)風(fēng)“杜鵑”生成過程中MCS 相關(guān)物理特征的演變規(guī)律,但是在TC 發(fā)生發(fā)展的過程中,對(duì)MCS 的生成及消亡的研究還不夠細(xì)致,并且這個(gè)過程實(shí)際上存在復(fù)雜的多尺度的相互作用。由于海上觀測(cè)資料的稀缺,初步分析得到的規(guī)律特征仍然需要更精細(xì)的觀測(cè)資料加以驗(yàn)證。

    猜你喜歡
    層云渦度杜鵑
    杜鵑紅
    心聲歌刊(2021年3期)2021-08-05 07:43:52
    中秋遣懷(新韻)
    杜鵑
    南海相對(duì)渦度的時(shí)空變化特征分析*
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    百里杜鵑百里歌
    民族音樂(2018年5期)2018-11-17 08:20:00
    百里杜鵑
    有壓迫感的層云
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    国产成人aa在线观看| 日日啪夜夜爽| 不卡视频在线观看欧美| 亚洲精品,欧美精品| 一级毛片 在线播放| 久久综合国产亚洲精品| 久久久久国产精品人妻一区二区| 联通29元200g的流量卡| 在线天堂最新版资源| 亚洲成人中文字幕在线播放| 有码 亚洲区| 亚洲精品国产av成人精品| 国产精品一二三区在线看| 成年女人在线观看亚洲视频 | 久久国内精品自在自线图片| 尤物成人国产欧美一区二区三区| 一个人看的www免费观看视频| 成人黄色视频免费在线看| 久久精品综合一区二区三区| 国产黄色视频一区二区在线观看| 另类亚洲欧美激情| 91午夜精品亚洲一区二区三区| 亚洲国产欧美在线一区| 精品一区二区三区视频在线| 午夜免费观看性视频| 永久网站在线| 日韩av免费高清视频| 高清欧美精品videossex| 丝瓜视频免费看黄片| 国产欧美亚洲国产| 一区二区av电影网| av卡一久久| 狠狠精品人妻久久久久久综合| 欧美人与善性xxx| 日日啪夜夜撸| 综合色丁香网| 亚洲成人一二三区av| 欧美老熟妇乱子伦牲交| 深夜a级毛片| 久久影院123| 少妇 在线观看| 成人欧美大片| 国产永久视频网站| 成人漫画全彩无遮挡| 狂野欧美激情性bbbbbb| 欧美日本视频| 日本午夜av视频| 色视频在线一区二区三区| 天堂俺去俺来也www色官网| 国产成人午夜福利电影在线观看| 国产精品无大码| 青春草亚洲视频在线观看| 51国产日韩欧美| 国产欧美日韩精品一区二区| 国产精品国产三级专区第一集| 一级毛片电影观看| 青春草国产在线视频| 另类亚洲欧美激情| 亚洲精品成人久久久久久| 91在线精品国自产拍蜜月| 干丝袜人妻中文字幕| 国产精品99久久久久久久久| 又大又黄又爽视频免费| 国语对白做爰xxxⅹ性视频网站| 久久精品国产亚洲av天美| 成人免费观看视频高清| 久久久久精品久久久久真实原创| 国模一区二区三区四区视频| 久久精品熟女亚洲av麻豆精品| 我的老师免费观看完整版| 精品午夜福利在线看| 国产精品熟女久久久久浪| 亚洲av成人精品一区久久| 欧美三级亚洲精品| 久久久国产一区二区| 老司机影院成人| 欧美97在线视频| 交换朋友夫妻互换小说| 99九九线精品视频在线观看视频| 中文字幕制服av| 大香蕉久久网| 免费不卡的大黄色大毛片视频在线观看| 国模一区二区三区四区视频| 久久这里有精品视频免费| 22中文网久久字幕| 午夜亚洲福利在线播放| 九九爱精品视频在线观看| 一个人看的www免费观看视频| 亚洲综合色惰| 国内精品宾馆在线| 国产精品国产av在线观看| 亚洲av中文av极速乱| 18禁在线播放成人免费| 亚洲美女搞黄在线观看| 三级国产精品欧美在线观看| 亚洲欧洲日产国产| 香蕉精品网在线| 看免费成人av毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美人与善性xxx| 亚洲精品成人久久久久久| 一个人观看的视频www高清免费观看| 欧美高清性xxxxhd video| 久久女婷五月综合色啪小说 | 欧美 日韩 精品 国产| 十八禁网站网址无遮挡 | 午夜福利在线在线| 一二三四中文在线观看免费高清| 亚洲国产精品国产精品| 永久网站在线| 日韩成人av中文字幕在线观看| 极品教师在线视频| 色综合色国产| 少妇人妻久久综合中文| 日韩亚洲欧美综合| 国产毛片在线视频| 在线观看一区二区三区| 国产片特级美女逼逼视频| 亚洲av中文字字幕乱码综合| 久久久午夜欧美精品| 深爱激情五月婷婷| 精品亚洲乱码少妇综合久久| 制服丝袜香蕉在线| 久久精品综合一区二区三区| 欧美成人午夜免费资源| 日韩一本色道免费dvd| 久久影院123| 亚洲av不卡在线观看| 婷婷色综合大香蕉| 国产精品爽爽va在线观看网站| av女优亚洲男人天堂| 晚上一个人看的免费电影| 大陆偷拍与自拍| 亚洲av免费在线观看| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av涩爱| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 色综合色国产| 夫妻午夜视频| 别揉我奶头 嗯啊视频| 高清午夜精品一区二区三区| 亚洲国产欧美人成| 亚洲精品国产av蜜桃| 99精国产麻豆久久婷婷| 少妇的逼好多水| 免费在线观看成人毛片| av天堂中文字幕网| 国产一区二区亚洲精品在线观看| 欧美高清性xxxxhd video| 欧美xxⅹ黑人| 久久综合国产亚洲精品| 亚洲av成人精品一区久久| 亚洲精品日韩在线中文字幕| av在线天堂中文字幕| 亚洲成人一二三区av| 免费av不卡在线播放| 五月玫瑰六月丁香| 男人狂女人下面高潮的视频| 亚洲电影在线观看av| 亚洲国产欧美人成| 精品亚洲乱码少妇综合久久| 免费看a级黄色片| 欧美zozozo另类| 国产亚洲av片在线观看秒播厂| 国产女主播在线喷水免费视频网站| 亚洲最大成人av| 国产精品蜜桃在线观看| 成人毛片60女人毛片免费| 国产高清国产精品国产三级 | 少妇高潮的动态图| 久久久久久久午夜电影| 高清毛片免费看| 日本-黄色视频高清免费观看| 国产亚洲午夜精品一区二区久久 | 亚洲精品乱久久久久久| 91精品国产九色| 国产老妇女一区| 大片电影免费在线观看免费| 国产视频内射| 精品国产一区二区三区久久久樱花 | 国产黄色免费在线视频| 日韩欧美精品免费久久| 伦精品一区二区三区| 精品国产乱码久久久久久小说| 成人无遮挡网站| 人人妻人人看人人澡| av在线观看视频网站免费| 国产成人精品婷婷| 美女国产视频在线观看| 久久久精品免费免费高清| 高清视频免费观看一区二区| av福利片在线观看| 99久久精品热视频| 91久久精品国产一区二区成人| 久久人人爽人人片av| 大又大粗又爽又黄少妇毛片口| 夫妻性生交免费视频一级片| 中文天堂在线官网| 爱豆传媒免费全集在线观看| 久久午夜福利片| 久久久久久久大尺度免费视频| 国产精品国产av在线观看| 天天躁夜夜躁狠狠久久av| 国产成人精品婷婷| 插阴视频在线观看视频| 亚洲国产成人一精品久久久| 中国美白少妇内射xxxbb| 深夜a级毛片| 午夜老司机福利剧场| 可以在线观看毛片的网站| 国产极品天堂在线| 国产黄片视频在线免费观看| 在线观看一区二区三区激情| 日韩中字成人| 亚洲精品日本国产第一区| 欧美人与善性xxx| av又黄又爽大尺度在线免费看| 日本黄色片子视频| 国产乱人视频| 禁无遮挡网站| 国产精品.久久久| 精品久久国产蜜桃| 国产一区有黄有色的免费视频| 一区二区三区精品91| 精品久久久精品久久久| 伊人久久精品亚洲午夜| 自拍欧美九色日韩亚洲蝌蚪91 | videossex国产| 中国美白少妇内射xxxbb| 亚州av有码| 国产精品久久久久久精品古装| av.在线天堂| 国产综合精华液| 国产精品麻豆人妻色哟哟久久| 国产午夜精品久久久久久一区二区三区| 一本久久精品| 国产毛片a区久久久久| 精品久久久久久久久av| 国产综合懂色| 美女主播在线视频| 国产乱来视频区| 精品熟女少妇av免费看| 欧美性感艳星| 美女cb高潮喷水在线观看| 国产爱豆传媒在线观看| 午夜日本视频在线| 国产精品成人在线| 国产高清有码在线观看视频| 亚洲精品成人av观看孕妇| 超碰97精品在线观看| 久久精品国产亚洲av天美| 人人妻人人爽人人添夜夜欢视频 | 日本wwww免费看| 男插女下体视频免费在线播放| 国产一区有黄有色的免费视频| 制服丝袜香蕉在线| 高清av免费在线| 国产精品一区www在线观看| 美女被艹到高潮喷水动态| 纵有疾风起免费观看全集完整版| 五月玫瑰六月丁香| 亚洲国产色片| 欧美区成人在线视频| 免费看a级黄色片| 日本av手机在线免费观看| 免费黄频网站在线观看国产| 少妇人妻 视频| 性色av一级| 1000部很黄的大片| 一区二区三区乱码不卡18| 一个人看视频在线观看www免费| 国产亚洲91精品色在线| 亚洲图色成人| 久久精品久久久久久噜噜老黄| 最近2019中文字幕mv第一页| 日韩欧美 国产精品| 亚洲国产最新在线播放| 成年版毛片免费区| 日本爱情动作片www.在线观看| 午夜视频国产福利| 男女边吃奶边做爰视频| 夜夜爽夜夜爽视频| 丰满人妻一区二区三区视频av| 晚上一个人看的免费电影| 日日啪夜夜撸| 亚洲精品成人av观看孕妇| 免费观看在线日韩| 亚洲精品影视一区二区三区av| 久久亚洲国产成人精品v| 国产国拍精品亚洲av在线观看| 高清毛片免费看| 少妇人妻一区二区三区视频| 麻豆精品久久久久久蜜桃| 中文天堂在线官网| 久久久久精品性色| av在线蜜桃| 卡戴珊不雅视频在线播放| 丝袜喷水一区| 久久精品人妻少妇| 99久久中文字幕三级久久日本| 观看免费一级毛片| 看十八女毛片水多多多| 在线观看三级黄色| 啦啦啦啦在线视频资源| 看十八女毛片水多多多| 97在线视频观看| 99久久人妻综合| 肉色欧美久久久久久久蜜桃 | 亚洲美女搞黄在线观看| 男人舔奶头视频| 性色avwww在线观看| 伦精品一区二区三区| 日本猛色少妇xxxxx猛交久久| 在线亚洲精品国产二区图片欧美 | 国产精品女同一区二区软件| 亚洲国产精品国产精品| 亚洲av国产av综合av卡| 女人十人毛片免费观看3o分钟| 黄色欧美视频在线观看| 日韩精品有码人妻一区| 51国产日韩欧美| 99视频精品全部免费 在线| 国产视频首页在线观看| 直男gayav资源| 边亲边吃奶的免费视频| 一个人观看的视频www高清免费观看| 久久久久久久亚洲中文字幕| 国产乱人视频| 亚洲高清免费不卡视频| a级一级毛片免费在线观看| 日韩伦理黄色片| 男人和女人高潮做爰伦理| 在线精品无人区一区二区三 | 久久久午夜欧美精品| 国产黄片视频在线免费观看| 国产成人一区二区在线| 久久久精品94久久精品| 日韩欧美一区视频在线观看 | 久久精品久久久久久噜噜老黄| 男人狂女人下面高潮的视频| 亚洲成人中文字幕在线播放| 免费看日本二区| 国产精品成人在线| 成人鲁丝片一二三区免费| 一区二区三区精品91| 人妻系列 视频| 中文字幕av成人在线电影| 亚洲欧美中文字幕日韩二区| 精品人妻偷拍中文字幕| 在线 av 中文字幕| 午夜福利网站1000一区二区三区| 男人爽女人下面视频在线观看| 国产免费福利视频在线观看| 老司机影院成人| 寂寞人妻少妇视频99o| 最近手机中文字幕大全| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 亚洲精品亚洲一区二区| 性色av一级| 黄片无遮挡物在线观看| 久久久久久久大尺度免费视频| 日韩中字成人| 久久精品国产a三级三级三级| 成人国产麻豆网| 国产爽快片一区二区三区| 精品国产乱码久久久久久小说| 成人亚洲精品av一区二区| 久久国产乱子免费精品| 亚洲久久久久久中文字幕| 成人国产麻豆网| 午夜精品国产一区二区电影 | 最近2019中文字幕mv第一页| 2022亚洲国产成人精品| 日本黄色片子视频| 看非洲黑人一级黄片| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| 精品一区二区免费观看| 免费黄频网站在线观看国产| 2018国产大陆天天弄谢| 亚洲av中文字字幕乱码综合| 搡老乐熟女国产| 国产亚洲av片在线观看秒播厂| 久久久久久久亚洲中文字幕| 午夜精品国产一区二区电影 | 91精品一卡2卡3卡4卡| 成人毛片60女人毛片免费| 免费播放大片免费观看视频在线观看| 一本一本综合久久| 身体一侧抽搐| 少妇人妻精品综合一区二区| 免费看a级黄色片| 女人被狂操c到高潮| 亚洲成人中文字幕在线播放| 成年av动漫网址| 精品人妻视频免费看| 免费看日本二区| av免费观看日本| 成年av动漫网址| 搡女人真爽免费视频火全软件| 国产欧美日韩一区二区三区在线 | 亚洲三级黄色毛片| 男女边吃奶边做爰视频| 国产伦在线观看视频一区| 在线观看美女被高潮喷水网站| 亚洲性久久影院| 最近中文字幕高清免费大全6| 岛国毛片在线播放| av国产久精品久网站免费入址| 久久久久久久久久久免费av| 精品久久久精品久久久| 男的添女的下面高潮视频| 水蜜桃什么品种好| 国产一区二区三区av在线| 成人亚洲精品av一区二区| 亚洲国产av新网站| 国产色婷婷99| 舔av片在线| 伦精品一区二区三区| 欧美激情在线99| 亚洲成色77777| 丰满少妇做爰视频| 中文天堂在线官网| 欧美3d第一页| 亚洲伊人久久精品综合| 成人无遮挡网站| 精品少妇黑人巨大在线播放| 在线观看免费高清a一片| 波野结衣二区三区在线| 国产免费一区二区三区四区乱码| 国产精品久久久久久久电影| 亚洲成人久久爱视频| 一区二区三区精品91| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 免费黄频网站在线观看国产| 欧美另类一区| 亚洲最大成人中文| 丰满少妇做爰视频| 日韩人妻高清精品专区| 女人被狂操c到高潮| 少妇人妻 视频| videos熟女内射| 亚洲av福利一区| 欧美高清成人免费视频www| 国国产精品蜜臀av免费| 十八禁网站网址无遮挡 | 69人妻影院| 网址你懂的国产日韩在线| 国语对白做爰xxxⅹ性视频网站| 边亲边吃奶的免费视频| 成人二区视频| 高清日韩中文字幕在线| 成人欧美大片| 最近中文字幕高清免费大全6| 一个人看的www免费观看视频| 2021少妇久久久久久久久久久| 日韩一本色道免费dvd| 免费看日本二区| 国产精品久久久久久精品电影| 大又大粗又爽又黄少妇毛片口| 人人妻人人看人人澡| 久久影院123| 看免费成人av毛片| 2021少妇久久久久久久久久久| 在线观看av片永久免费下载| 亚洲精品中文字幕在线视频 | 色网站视频免费| 中文乱码字字幕精品一区二区三区| 男女边摸边吃奶| 韩国高清视频一区二区三区| 精品国产乱码久久久久久小说| av在线app专区| 色哟哟·www| 日韩国内少妇激情av| 久久午夜福利片| 亚洲av欧美aⅴ国产| 国产精品久久久久久精品古装| 草草在线视频免费看| 老女人水多毛片| 国产精品三级大全| 国产成人精品福利久久| 日本欧美国产在线视频| 中文欧美无线码| 97在线人人人人妻| 国产精品久久久久久精品电影| 免费黄网站久久成人精品| 五月开心婷婷网| 国产精品.久久久| 深爱激情五月婷婷| 欧美97在线视频| 黑人高潮一二区| 亚洲不卡免费看| 成人黄色视频免费在线看| 国产成人aa在线观看| 成年女人看的毛片在线观看| 一级黄片播放器| 成人国产麻豆网| 夜夜爽夜夜爽视频| 最近中文字幕高清免费大全6| 九九在线视频观看精品| 亚洲成人久久爱视频| 国产一区二区在线观看日韩| 久久久久久久亚洲中文字幕| 久久影院123| 人人妻人人澡人人爽人人夜夜| 老司机影院成人| 少妇人妻一区二区三区视频| 在线观看人妻少妇| av在线天堂中文字幕| 老司机影院毛片| 51国产日韩欧美| 精品国产三级普通话版| 精品久久久噜噜| 黄色怎么调成土黄色| 成年女人看的毛片在线观看| 成人鲁丝片一二三区免费| 久久久久久久久大av| kizo精华| 亚洲综合精品二区| 国产永久视频网站| 国产爽快片一区二区三区| 九色成人免费人妻av| 十八禁网站网址无遮挡 | 精品久久久久久久久av| 日本熟妇午夜| 亚洲四区av| 亚洲精品日韩av片在线观看| 午夜精品一区二区三区免费看| 在线观看美女被高潮喷水网站| 最近中文字幕2019免费版| 人妻系列 视频| 国产精品一及| kizo精华| 亚洲精品国产av成人精品| 一区二区三区精品91| 又爽又黄a免费视频| 我要看日韩黄色一级片| 在线a可以看的网站| 黄色视频在线播放观看不卡| 亚洲美女视频黄频| 人妻 亚洲 视频| 国精品久久久久久国模美| 日本黄大片高清| 女的被弄到高潮叫床怎么办| 国产成人精品婷婷| 国产成人精品福利久久| 国产黄色视频一区二区在线观看| 亚洲欧美成人综合另类久久久| 可以在线观看毛片的网站| 在线观看免费高清a一片| 久久久久久九九精品二区国产| 免费黄网站久久成人精品| 久久久久久久午夜电影| 久久久久久久久久人人人人人人| 亚洲人与动物交配视频| 99热网站在线观看| 伦精品一区二区三区| 好男人在线观看高清免费视频| 少妇被粗大猛烈的视频| 男人添女人高潮全过程视频| 久久韩国三级中文字幕| 亚洲精品久久久久久婷婷小说| videos熟女内射| 一级二级三级毛片免费看| 国产熟女欧美一区二区| 蜜桃亚洲精品一区二区三区| 一级片'在线观看视频| 免费大片黄手机在线观看| 亚洲欧美一区二区三区国产| 99热网站在线观看| 91狼人影院| 日本欧美国产在线视频| 男女下面进入的视频免费午夜| 午夜免费男女啪啪视频观看| 丰满少妇做爰视频| 国产 一区精品| 午夜免费观看性视频| 六月丁香七月| 免费高清在线观看视频在线观看| 在线观看人妻少妇| 91在线精品国自产拍蜜月| 亚洲av欧美aⅴ国产| 久久精品国产鲁丝片午夜精品| 久久人人爽av亚洲精品天堂 | 中文精品一卡2卡3卡4更新| 亚洲最大成人中文| 乱系列少妇在线播放| 久久久久久久久大av| 十八禁网站网址无遮挡 | 免费观看的影片在线观看| 国产精品福利在线免费观看| 国产探花在线观看一区二区| 性插视频无遮挡在线免费观看| 久久99蜜桃精品久久| 国产日韩欧美亚洲二区| 国产黄片视频在线免费观看| 久久久亚洲精品成人影院| 久久久a久久爽久久v久久| 欧美激情在线99| 免费高清在线观看视频在线观看| 夜夜爽夜夜爽视频| videossex国产| 97人妻精品一区二区三区麻豆| 成人午夜精彩视频在线观看| 97热精品久久久久久| 免费观看a级毛片全部| 久久这里有精品视频免费| 一本久久精品| 在线观看一区二区三区| 中国美白少妇内射xxxbb| 国产精品久久久久久精品电影小说 |