• <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
    日本wwww免费看| 免费高清在线观看视频在线观看| 丁香六月天网| 少妇精品久久久久久久| 久久精品国产a三级三级三级| 自线自在国产av| 欧美成人午夜免费资源| 欧美av亚洲av综合av国产av | 蜜桃国产av成人99| 亚洲国产精品国产精品| videos熟女内射| 成年动漫av网址| 亚洲国产精品999| 黄色一级大片看看| 欧美激情高清一区二区三区 | 国产成人精品久久久久久| 99国产综合亚洲精品| 成人免费观看视频高清| 亚洲欧美清纯卡通| 久久 成人 亚洲| av国产久精品久网站免费入址| 美国免费a级毛片| 欧美人与性动交α欧美软件| 亚洲精品成人av观看孕妇| 香蕉精品网在线| 精品人妻一区二区三区麻豆| 日本91视频免费播放| 久久久久久久大尺度免费视频| www.自偷自拍.com| 巨乳人妻的诱惑在线观看| 久热这里只有精品99| a级片在线免费高清观看视频| 国产在线视频一区二区| 久久久久久久久久人人人人人人| av.在线天堂| 1024视频免费在线观看| 精品一区二区免费观看| 青春草亚洲视频在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲综合精品二区| 国产 精品1| 中文字幕制服av| 天天影视国产精品| 国产毛片在线视频| 亚洲av电影在线进入| 日本免费在线观看一区| 大香蕉久久成人网| 久久这里只有精品19| 欧美日韩视频精品一区| av网站在线播放免费| 综合色丁香网| 国产精品嫩草影院av在线观看| 精品国产超薄肉色丝袜足j| 亚洲国产成人一精品久久久| tube8黄色片| 欧美黄色片欧美黄色片| 久久女婷五月综合色啪小说| 男女啪啪激烈高潮av片| 精品久久久久久电影网| 日韩av免费高清视频| 久久人人爽av亚洲精品天堂| 国产av国产精品国产| 美国免费a级毛片| 麻豆精品久久久久久蜜桃| www.自偷自拍.com| 国产午夜精品一二区理论片| 极品少妇高潮喷水抽搐| 18禁国产床啪视频网站| 久久久久精品人妻al黑| 狠狠婷婷综合久久久久久88av| 亚洲精品aⅴ在线观看| 18在线观看网站| 免费高清在线观看视频在线观看| 亚洲国产av影院在线观看| 多毛熟女@视频| 在线观看国产h片| 黄色怎么调成土黄色| 在现免费观看毛片| 久久久精品区二区三区| 国产一区有黄有色的免费视频| 高清视频免费观看一区二区| 女性生殖器流出的白浆| 电影成人av| 女人精品久久久久毛片| 午夜激情久久久久久久| 精品第一国产精品| av国产久精品久网站免费入址| 丝袜脚勾引网站| 成人二区视频| 天天影视国产精品| videosex国产| 在线看a的网站| 亚洲熟女精品中文字幕| 久久久久久伊人网av| 我的亚洲天堂| 中文字幕精品免费在线观看视频| 少妇 在线观看| 国产淫语在线视频| 成人免费观看视频高清| 永久免费av网站大全| 99热国产这里只有精品6| 中文精品一卡2卡3卡4更新| 有码 亚洲区| 一二三四中文在线观看免费高清| 久久 成人 亚洲| a 毛片基地| 飞空精品影院首页| 少妇被粗大的猛进出69影院| 亚洲久久久国产精品| 日本免费在线观看一区| 两性夫妻黄色片| 国产精品秋霞免费鲁丝片| 国产人伦9x9x在线观看 | www.av在线官网国产| 精品福利永久在线观看| 青草久久国产| 如日韩欧美国产精品一区二区三区| 性少妇av在线| 男女国产视频网站| 91精品伊人久久大香线蕉| 777久久人妻少妇嫩草av网站| 18禁动态无遮挡网站| 欧美97在线视频| 老汉色av国产亚洲站长工具| 寂寞人妻少妇视频99o| 男人操女人黄网站| 欧美人与善性xxx| av免费在线看不卡| 久久青草综合色| 日韩 亚洲 欧美在线| 乱人伦中国视频| 亚洲精华国产精华液的使用体验| 伦理电影大哥的女人| 免费女性裸体啪啪无遮挡网站| 成人黄色视频免费在线看| 精品一区二区三区四区五区乱码 | 亚洲精品日本国产第一区| 国产精品久久久久久精品古装| 国产xxxxx性猛交| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄频高清免费视频| 久久精品亚洲av国产电影网| 一本色道久久久久久精品综合| 亚洲欧洲日产国产| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区| 国产乱来视频区| 天天躁狠狠躁夜夜躁狠狠躁| 久久婷婷青草| 99国产综合亚洲精品| 免费观看a级毛片全部| 国产精品人妻久久久影院| 久久99热这里只频精品6学生| 韩国av在线不卡| 老女人水多毛片| videosex国产| 香蕉国产在线看| 一级片'在线观看视频| 免费观看性生交大片5| 少妇熟女欧美另类| 一区二区日韩欧美中文字幕| 人人妻人人澡人人爽人人夜夜| 少妇被粗大猛烈的视频| 亚洲精品国产色婷婷电影| 国产精品久久久久成人av| 久久人人97超碰香蕉20202| 久久久欧美国产精品| 男女啪啪激烈高潮av片| 一级毛片电影观看| 国产精品99久久99久久久不卡 | 午夜福利一区二区在线看| 亚洲国产精品一区二区三区在线| 亚洲国产av新网站| 大片免费播放器 马上看| 女人被躁到高潮嗷嗷叫费观| 宅男免费午夜| 国产黄色免费在线视频| 亚洲内射少妇av| 老司机亚洲免费影院| 国产在线一区二区三区精| 成年人免费黄色播放视频| 日本av免费视频播放| www.熟女人妻精品国产| 人人妻人人添人人爽欧美一区卜| 韩国高清视频一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 如日韩欧美国产精品一区二区三区| 精品一区在线观看国产| 9热在线视频观看99| 国产乱人偷精品视频| 熟女电影av网| 精品午夜福利在线看| 国精品久久久久久国模美| 亚洲欧美一区二区三区国产| 青春草视频在线免费观看| 成人毛片60女人毛片免费| 中文精品一卡2卡3卡4更新| 一本久久精品| 伊人久久国产一区二区| 黑丝袜美女国产一区| 国产极品粉嫩免费观看在线| 91精品伊人久久大香线蕉| 天天躁夜夜躁狠狠躁躁| 久久久久久久久久久免费av| 国产极品粉嫩免费观看在线| 天美传媒精品一区二区| 久久久久久久大尺度免费视频| 亚洲精品美女久久av网站| 中文欧美无线码| 精品少妇久久久久久888优播| 久久精品久久久久久久性| 老汉色av国产亚洲站长工具| 美女主播在线视频| 亚洲精品久久午夜乱码| 伊人久久大香线蕉亚洲五| 色网站视频免费| 亚洲国产欧美在线一区| 99久久精品国产国产毛片| av在线老鸭窝| 免费观看在线日韩| 日韩大片免费观看网站| 2021少妇久久久久久久久久久| 少妇精品久久久久久久| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 人妻少妇偷人精品九色| 亚洲精品久久成人aⅴ小说| 亚洲精品aⅴ在线观看| 日本色播在线视频| 观看美女的网站| 国产欧美日韩一区二区三区在线| 亚洲中文av在线| av免费在线看不卡| 国产精品久久久久成人av| 男女无遮挡免费网站观看| 成人毛片a级毛片在线播放| 亚洲综合色网址| 久久99热这里只频精品6学生| 纯流量卡能插随身wifi吗| 91精品三级在线观看| 亚洲国产av新网站| 日韩免费高清中文字幕av| 美女中出高潮动态图| 秋霞在线观看毛片| 亚洲第一av免费看| 高清黄色对白视频在线免费看| 久久久国产一区二区| 免费黄色在线免费观看| 久久久久精品人妻al黑| 人人妻人人澡人人看| 久久国内精品自在自线图片| 久久精品熟女亚洲av麻豆精品| a级毛片在线看网站| 人人澡人人妻人| 美女视频免费永久观看网站| 新久久久久国产一级毛片| 午夜免费男女啪啪视频观看| 亚洲欧美精品自产自拍| 欧美日韩视频精品一区| 亚洲国产最新在线播放| tube8黄色片| 男的添女的下面高潮视频| 多毛熟女@视频| 三级国产精品片| 国产探花极品一区二区| 日本色播在线视频| 日韩免费高清中文字幕av| 少妇的逼水好多| 欧美bdsm另类| 美国免费a级毛片| 日韩视频在线欧美| 亚洲人成电影观看| 久久精品久久精品一区二区三区| 久久精品国产a三级三级三级| 老鸭窝网址在线观看| 搡老乐熟女国产| 天堂俺去俺来也www色官网| 99re6热这里在线精品视频| 美女主播在线视频| 国产成人免费无遮挡视频| 日日啪夜夜爽| 高清黄色对白视频在线免费看| 一级毛片电影观看| 一本久久精品| 男女无遮挡免费网站观看| 91午夜精品亚洲一区二区三区| 在线观看国产h片| 亚洲av综合色区一区| 欧美日本中文国产一区发布| 少妇被粗大猛烈的视频| 午夜日韩欧美国产| 久久久a久久爽久久v久久| 国产成人精品无人区| 自拍欧美九色日韩亚洲蝌蚪91| 99久国产av精品国产电影| 午夜福利乱码中文字幕| 亚洲av国产av综合av卡| 国产综合精华液| 女的被弄到高潮叫床怎么办| videosex国产| 777米奇影视久久| 国产亚洲一区二区精品| 中文字幕人妻丝袜一区二区 | 亚洲国产av新网站| 高清欧美精品videossex| 精品一区二区三卡| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂| 久久精品国产亚洲av天美| 免费女性裸体啪啪无遮挡网站| 亚洲内射少妇av| 老汉色∧v一级毛片| 国产成人精品在线电影| 国产精品亚洲av一区麻豆 | 在线观看www视频免费| 亚洲成色77777| 国产亚洲欧美精品永久| 日韩中文字幕欧美一区二区 | 精品久久久精品久久久| 国产av一区二区精品久久| 一区二区三区精品91| 少妇猛男粗大的猛烈进出视频| 国产成人精品一,二区| 国产熟女欧美一区二区| 欧美少妇被猛烈插入视频| 中文字幕亚洲精品专区| 成年美女黄网站色视频大全免费| 亚洲情色 制服丝袜| 高清在线视频一区二区三区| 欧美精品一区二区大全| av在线播放精品| 天天影视国产精品| 美女大奶头黄色视频| 91在线精品国自产拍蜜月| 丝瓜视频免费看黄片| 最近最新中文字幕大全免费视频 | 久久精品人人爽人人爽视色| √禁漫天堂资源中文www| 热re99久久精品国产66热6| 国产高清不卡午夜福利| 一边亲一边摸免费视频| 国产成人欧美| 欧美日韩综合久久久久久| 国产精品一区二区在线不卡| 999精品在线视频| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩一区二区三区在线| 成人毛片60女人毛片免费| 久久影院123| 999久久久国产精品视频| 中文天堂在线官网| av国产精品久久久久影院| 80岁老熟妇乱子伦牲交| 最黄视频免费看| 美国免费a级毛片| 亚洲国产精品999| 18禁裸乳无遮挡动漫免费视频| 一二三四在线观看免费中文在| 久久久久久久久久久免费av| 国产亚洲午夜精品一区二区久久| 欧美 日韩 精品 国产| 亚洲一码二码三码区别大吗| 亚洲情色 制服丝袜| 老女人水多毛片| 久久精品亚洲av国产电影网| 丰满少妇做爰视频| 亚洲人成网站在线观看播放| 一区在线观看完整版| 热99国产精品久久久久久7| 国产在线免费精品| 亚洲国产欧美在线一区| 久久久久久免费高清国产稀缺| 日本-黄色视频高清免费观看| 美女xxoo啪啪120秒动态图| 久久狼人影院| 26uuu在线亚洲综合色| 曰老女人黄片| 国产在线免费精品| 夫妻午夜视频| 纯流量卡能插随身wifi吗| 日韩av在线免费看完整版不卡| 亚洲国产欧美网| 国产成人午夜福利电影在线观看| 久久久精品国产亚洲av高清涩受| 亚洲久久久国产精品| 成人国产麻豆网| 秋霞在线观看毛片| 亚洲成人av在线免费| 在线观看国产h片| 女的被弄到高潮叫床怎么办| 亚洲av欧美aⅴ国产| 色播在线永久视频| 中文字幕人妻丝袜一区二区 | av一本久久久久| 国产日韩欧美亚洲二区| 亚洲精品日韩在线中文字幕| 99久久综合免费| 国产又色又爽无遮挡免| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲一区二区精品| 狠狠婷婷综合久久久久久88av| 久久久久人妻精品一区果冻| 巨乳人妻的诱惑在线观看| 人人妻人人澡人人看| 亚洲精品在线美女| 成人毛片a级毛片在线播放| 99九九在线精品视频| 老汉色av国产亚洲站长工具| 国产成人精品婷婷| av线在线观看网站| 久久精品亚洲av国产电影网| 美国免费a级毛片| 九草在线视频观看| 2018国产大陆天天弄谢| 欧美 日韩 精品 国产| 汤姆久久久久久久影院中文字幕| 天天躁日日躁夜夜躁夜夜| 欧美老熟妇乱子伦牲交| 中文精品一卡2卡3卡4更新| 亚洲国产成人一精品久久久| 1024视频免费在线观看| 亚洲av免费高清在线观看| 国产精品.久久久| 美女xxoo啪啪120秒动态图| 国产av一区二区精品久久| 最近最新中文字幕免费大全7| 欧美亚洲日本最大视频资源| 男女啪啪激烈高潮av片| 午夜福利影视在线免费观看| 99国产精品免费福利视频| 欧美另类一区| 少妇精品久久久久久久| 一区二区av电影网| 韩国高清视频一区二区三区| 亚洲av日韩在线播放| 天天躁日日躁夜夜躁夜夜| 国产免费又黄又爽又色| 99久久中文字幕三级久久日本| 女的被弄到高潮叫床怎么办| 超色免费av| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 国产av国产精品国产| 尾随美女入室| 亚洲欧洲精品一区二区精品久久久 | 国产乱来视频区| 美女大奶头黄色视频| www日本在线高清视频| 你懂的网址亚洲精品在线观看| 日日摸夜夜添夜夜爱| 欧美激情 高清一区二区三区| 欧美少妇被猛烈插入视频| 美女午夜性视频免费| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 亚洲天堂av无毛| 亚洲三区欧美一区| 男女啪啪激烈高潮av片| 久久精品国产鲁丝片午夜精品| 精品一区二区三卡| 26uuu在线亚洲综合色| 国产一级毛片在线| 这个男人来自地球电影免费观看 | 日韩一区二区三区影片| 美女大奶头黄色视频| 一区二区三区乱码不卡18| 国产精品国产av在线观看| 色哟哟·www| 国产亚洲精品第一综合不卡| 又大又黄又爽视频免费| 美女视频免费永久观看网站| 亚洲精品久久久久久婷婷小说| 人妻少妇偷人精品九色| 久久99一区二区三区| 国产精品久久久久成人av| 最近中文字幕2019免费版| 免费观看性生交大片5| 亚洲成av片中文字幕在线观看 | 一边摸一边做爽爽视频免费| 制服丝袜香蕉在线| 嫩草影院入口| 美女高潮到喷水免费观看| 制服人妻中文乱码| 久久久欧美国产精品| 黑人猛操日本美女一级片| 精品视频人人做人人爽| 亚洲av福利一区| 精品亚洲成a人片在线观看| 亚洲av成人精品一二三区| 三上悠亚av全集在线观看| 日本色播在线视频| 亚洲成人一二三区av| 最近最新中文字幕大全免费视频 | 在线观看免费日韩欧美大片| 久久人人爽人人片av| a级毛片在线看网站| 日韩 亚洲 欧美在线| 男女下面插进去视频免费观看| 免费日韩欧美在线观看| 黄频高清免费视频| 亚洲四区av| 亚洲成人av在线免费| 大片免费播放器 马上看| 看十八女毛片水多多多| 国产1区2区3区精品| 日韩,欧美,国产一区二区三区| 午夜免费鲁丝| 黄片小视频在线播放| 久久综合国产亚洲精品| 国产不卡av网站在线观看| 97在线视频观看| 少妇的逼水好多| 国产男人的电影天堂91| 亚洲国产看品久久| 国产精品人妻久久久影院| 日韩精品免费视频一区二区三区| 亚洲欧洲日产国产| 国产精品二区激情视频| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 成人国语在线视频| 亚洲国产看品久久| 精品午夜福利在线看| 天天影视国产精品| 一级片'在线观看视频| h视频一区二区三区| 男女高潮啪啪啪动态图| 亚洲国产精品999| 久久久久久久久免费视频了| 人妻一区二区av| 我的亚洲天堂| 老汉色av国产亚洲站长工具| 啦啦啦中文免费视频观看日本| 老汉色av国产亚洲站长工具| av在线播放精品| 国产精品.久久久| 在线观看免费高清a一片| 波野结衣二区三区在线| 丰满少妇做爰视频| 免费观看av网站的网址| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 亚洲欧美中文字幕日韩二区| 日韩人妻精品一区2区三区| 97人妻天天添夜夜摸| 亚洲一区二区三区欧美精品| 久久久精品区二区三区| 婷婷色综合大香蕉| 黑人巨大精品欧美一区二区蜜桃| 日韩不卡一区二区三区视频在线| av在线app专区| 国产成人精品婷婷| 亚洲精品久久成人aⅴ小说| 在线观看国产h片| 亚洲精品美女久久av网站| 91精品伊人久久大香线蕉| 熟女电影av网| 久久久久网色| 麻豆av在线久日| 妹子高潮喷水视频| 又大又黄又爽视频免费| 午夜福利视频在线观看免费| 最近最新中文字幕免费大全7| 国产xxxxx性猛交| 国产亚洲最大av| 欧美人与善性xxx| 91精品国产国语对白视频| 久久国产精品大桥未久av| 国产精品熟女久久久久浪| 老鸭窝网址在线观看| 亚洲精品第二区| 两性夫妻黄色片| 自线自在国产av| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 亚洲内射少妇av| 欧美激情 高清一区二区三区| 久久青草综合色| 亚洲成色77777| 亚洲av免费高清在线观看| 日韩av不卡免费在线播放| 亚洲av综合色区一区| 国产麻豆69| 精品一区二区三卡| 中国国产av一级| av有码第一页| 不卡av一区二区三区| 老鸭窝网址在线观看| 国产一级毛片在线| 精品久久久久久电影网| 亚洲国产看品久久| 国产一区二区 视频在线| 国产在线免费精品| 亚洲精品国产av成人精品| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区国产| 国产精品国产三级专区第一集| 亚洲国产成人一精品久久久| 性少妇av在线| 午夜福利乱码中文字幕| 黄片播放在线免费| 精品酒店卫生间| 高清欧美精品videossex| 日韩人妻精品一区2区三区| 妹子高潮喷水视频| 精品国产一区二区久久| 欧美97在线视频| 成人漫画全彩无遮挡| 性高湖久久久久久久久免费观看| 美国免费a级毛片| 亚洲av中文av极速乱| 久久久久久久大尺度免费视频| 一级毛片 在线播放| 两个人免费观看高清视频|