• <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
    成人手机av| 欧美日韩精品网址| 色在线成人网| 国产一区二区三区在线臀色熟女 | 久久中文字幕一级| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 午夜免费成人在线视频| 国产在线视频一区二区| tocl精华| 亚洲精品av麻豆狂野| 欧美成人午夜精品| 人人妻人人爽人人添夜夜欢视频| 成年人黄色毛片网站| 亚洲精品国产色婷婷电影| 欧美精品啪啪一区二区三区| 久久精品熟女亚洲av麻豆精品| 99国产精品免费福利视频| 成年人免费黄色播放视频| 十八禁网站免费在线| 99精品欧美一区二区三区四区| 欧美日本中文国产一区发布| 久久九九热精品免费| 久久久久精品国产欧美久久久| 757午夜福利合集在线观看| 亚洲专区中文字幕在线| 午夜福利在线免费观看网站| 99香蕉大伊视频| 丰满迷人的少妇在线观看| 亚洲专区国产一区二区| 美女高潮喷水抽搐中文字幕| 免费观看a级毛片全部| 最近最新中文字幕大全电影3 | 在线观看www视频免费| 久久天堂一区二区三区四区| 青草久久国产| 日日爽夜夜爽网站| 亚洲精品在线观看二区| 国产免费现黄频在线看| 两人在一起打扑克的视频| 午夜福利视频精品| 黑人巨大精品欧美一区二区蜜桃| 国产精品av久久久久免费| 国产麻豆69| 精品高清国产在线一区| 日本五十路高清| 国产男女超爽视频在线观看| 亚洲成国产人片在线观看| 99九九在线精品视频| 日韩欧美国产一区二区入口| 国产一区二区三区在线臀色熟女 | 制服诱惑二区| 在线观看免费高清a一片| 国产高清视频在线播放一区| 成年人免费黄色播放视频| 亚洲一卡2卡3卡4卡5卡精品中文| 王馨瑶露胸无遮挡在线观看| 亚洲精品在线观看二区| 久久99一区二区三区| 可以免费在线观看a视频的电影网站| 一本大道久久a久久精品| 亚洲少妇的诱惑av| 亚洲一码二码三码区别大吗| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 日本欧美视频一区| 久久99一区二区三区| 2018国产大陆天天弄谢| 建设人人有责人人尽责人人享有的| 少妇被粗大的猛进出69影院| 久久性视频一级片| 亚洲 欧美一区二区三区| 亚洲国产欧美日韩在线播放| 日本黄色日本黄色录像| 99热网站在线观看| 777久久人妻少妇嫩草av网站| 亚洲国产毛片av蜜桃av| 欧美黑人欧美精品刺激| 韩国精品一区二区三区| av天堂久久9| av不卡在线播放| 妹子高潮喷水视频| 人妻 亚洲 视频| 国产一区二区三区在线臀色熟女 | 国产区一区二久久| 两个人看的免费小视频| 咕卡用的链子| 午夜精品国产一区二区电影| 人妻 亚洲 视频| 亚洲精品中文字幕一二三四区 | 日韩欧美国产一区二区入口| 天堂动漫精品| a级毛片在线看网站| 国产单亲对白刺激| 国产成人精品久久二区二区免费| 黑人操中国人逼视频| 超碰97精品在线观看| 另类亚洲欧美激情| 亚洲国产欧美一区二区综合| 日韩中文字幕视频在线看片| 精品少妇内射三级| 国产高清视频在线播放一区| 国产伦理片在线播放av一区| 免费在线观看黄色视频的| 久久影院123| 蜜桃国产av成人99| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 亚洲一区中文字幕在线| 日本精品一区二区三区蜜桃| 两人在一起打扑克的视频| 欧美国产精品va在线观看不卡| 午夜福利欧美成人| 国产色视频综合| 亚洲 欧美一区二区三区| 人人妻,人人澡人人爽秒播| 国产精品av久久久久免费| 老司机福利观看| 下体分泌物呈黄色| 欧美老熟妇乱子伦牲交| 极品人妻少妇av视频| 免费在线观看视频国产中文字幕亚洲| 欧美黄色片欧美黄色片| 性少妇av在线| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 国产欧美日韩综合在线一区二区| 国产精品熟女久久久久浪| 一区二区日韩欧美中文字幕| 一区在线观看完整版| 少妇猛男粗大的猛烈进出视频| 免费看a级黄色片| 51午夜福利影视在线观看| 国产精品 欧美亚洲| 精品一区二区三区视频在线观看免费 | 18禁黄网站禁片午夜丰满| 国产欧美日韩精品亚洲av| 欧美老熟妇乱子伦牲交| 亚洲精品粉嫩美女一区| 日本a在线网址| 国产精品偷伦视频观看了| 成人国语在线视频| 丁香六月欧美| 国产精品亚洲一级av第二区| 亚洲国产毛片av蜜桃av| 啦啦啦 在线观看视频| 亚洲黑人精品在线| 精品久久久久久电影网| 亚洲精品国产色婷婷电影| 夜夜骑夜夜射夜夜干| 国产成人欧美| 国产一卡二卡三卡精品| 在线看a的网站| 午夜激情久久久久久久| 无限看片的www在线观看| 亚洲少妇的诱惑av| 国产色视频综合| 午夜91福利影院| 国产主播在线观看一区二区| 亚洲国产av新网站| 亚洲精品成人av观看孕妇| 国产一区二区激情短视频| 久久精品国产99精品国产亚洲性色 | 美女福利国产在线| 国产免费av片在线观看野外av| 中文欧美无线码| 精品国产国语对白av| 欧美黑人精品巨大| 久9热在线精品视频| 亚洲中文日韩欧美视频| 热99久久久久精品小说推荐| 国产成人免费观看mmmm| 亚洲精品中文字幕在线视频| 欧美精品一区二区大全| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 欧美日韩成人在线一区二区| 亚洲自偷自拍图片 自拍| 国产精品国产av在线观看| 欧美日韩福利视频一区二区| 欧美激情久久久久久爽电影 | 热re99久久国产66热| av在线播放免费不卡| 亚洲熟女精品中文字幕| a级毛片黄视频| 亚洲成av片中文字幕在线观看| 免费不卡黄色视频| 国产视频一区二区在线看| 色94色欧美一区二区| 人妻久久中文字幕网| 妹子高潮喷水视频| 夜夜爽天天搞| 咕卡用的链子| 欧美日本中文国产一区发布| 久久久精品94久久精品| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 搡老熟女国产l中国老女人| 另类亚洲欧美激情| 亚洲第一欧美日韩一区二区三区 | 欧美在线一区亚洲| 成人精品一区二区免费| 国产精品国产av在线观看| 国产麻豆69| 精品少妇一区二区三区视频日本电影| 精品免费久久久久久久清纯 | 女人被躁到高潮嗷嗷叫费观| 色精品久久人妻99蜜桃| 午夜日韩欧美国产| 国产精品久久久人人做人人爽| 一进一出抽搐动态| 亚洲成人免费电影在线观看| 国产精品一区二区在线观看99| 午夜久久久在线观看| 久久人人97超碰香蕉20202| av超薄肉色丝袜交足视频| 男男h啪啪无遮挡| 久久久精品免费免费高清| 黄色片一级片一级黄色片| 国产精品1区2区在线观看. | 亚洲av第一区精品v没综合| 欧美精品av麻豆av| 免费看a级黄色片| 亚洲五月婷婷丁香| 九色亚洲精品在线播放| 久久人妻福利社区极品人妻图片| 日日爽夜夜爽网站| 久久天堂一区二区三区四区| 欧美另类亚洲清纯唯美| 一个人免费看片子| 热99re8久久精品国产| av免费在线观看网站| 老司机影院毛片| 久久影院123| 丝袜喷水一区| 国产成人影院久久av| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲专区中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 久久精品91无色码中文字幕| 国产又爽黄色视频| 母亲3免费完整高清在线观看| 久久精品国产99精品国产亚洲性色 | 男女边摸边吃奶| 熟女少妇亚洲综合色aaa.| av不卡在线播放| 丁香欧美五月| 黄频高清免费视频| 婷婷丁香在线五月| 国产精品久久久久成人av| 国产精品麻豆人妻色哟哟久久| 一夜夜www| 国产亚洲欧美精品永久| 香蕉久久夜色| 午夜福利视频在线观看免费| 叶爱在线成人免费视频播放| 人妻久久中文字幕网| 亚洲三区欧美一区| 女人久久www免费人成看片| 91麻豆av在线| 黄色视频在线播放观看不卡| 国产av一区二区精品久久| 青草久久国产| 亚洲欧美一区二区三区久久| 国产又爽黄色视频| 亚洲黑人精品在线| 精品第一国产精品| 91字幕亚洲| 亚洲av日韩在线播放| 最新在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 丁香六月天网| 久久99一区二区三区| 国产av国产精品国产| 12—13女人毛片做爰片一| 十八禁网站免费在线| 亚洲精品乱久久久久久| 亚洲成人免费电影在线观看| 黄色视频,在线免费观看| 久久影院123| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品成人在线| 又黄又粗又硬又大视频| 亚洲久久久国产精品| 日韩中文字幕视频在线看片| 亚洲五月色婷婷综合| av在线播放免费不卡| 一区二区av电影网| 欧美乱码精品一区二区三区| 蜜桃在线观看..| 无人区码免费观看不卡 | 国产精品一区二区在线观看99| 一夜夜www| 午夜激情久久久久久久| 国产老妇伦熟女老妇高清| 精品一品国产午夜福利视频| 久久久久国内视频| 国产精品一区二区在线不卡| 激情视频va一区二区三区| 19禁男女啪啪无遮挡网站| 国产一区二区在线观看av| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av香蕉五月 | 国产熟女午夜一区二区三区| 99国产精品99久久久久| 丁香六月天网| 18在线观看网站| 精品少妇内射三级| 亚洲av日韩精品久久久久久密| 91老司机精品| 亚洲视频免费观看视频| 国产精品一区二区在线不卡| 宅男免费午夜| 搡老乐熟女国产| 亚洲精品久久成人aⅴ小说| 狂野欧美激情性xxxx| 人人妻,人人澡人人爽秒播| 中文字幕人妻熟女乱码| 三上悠亚av全集在线观看| 制服诱惑二区| 成人av一区二区三区在线看| 最黄视频免费看| 久久久久久免费高清国产稀缺| 久久精品成人免费网站| 亚洲欧美精品综合一区二区三区| 久久免费观看电影| 国产精品1区2区在线观看. | 国产熟女午夜一区二区三区| 亚洲精品在线美女| 制服诱惑二区| 欧美日韩亚洲高清精品| av福利片在线| 欧美+亚洲+日韩+国产| 国产av又大| 亚洲中文字幕日韩| 国产精品成人在线| 蜜桃国产av成人99| 麻豆国产av国片精品| 亚洲欧美激情在线| 国产精品成人在线| 美女视频免费永久观看网站| 国产欧美日韩一区二区精品| 在线观看www视频免费| 亚洲人成电影观看| 少妇 在线观看| 精品福利永久在线观看| 一级a爱视频在线免费观看| 亚洲av日韩在线播放| av不卡在线播放| 免费观看av网站的网址| 欧美国产精品va在线观看不卡| 国产亚洲精品久久久久5区| 脱女人内裤的视频| 岛国毛片在线播放| 在线观看免费高清a一片| 亚洲中文av在线| 丁香六月欧美| 亚洲人成77777在线视频| 亚洲成人手机| 大陆偷拍与自拍| 99久久99久久久精品蜜桃| 色94色欧美一区二区| 一级毛片女人18水好多| 久久亚洲精品不卡| 国产成人免费观看mmmm| 国产欧美日韩一区二区三区在线| 一进一出抽搐动态| 十八禁高潮呻吟视频| 精品一区二区三区四区五区乱码| 涩涩av久久男人的天堂| 亚洲人成伊人成综合网2020| 女人高潮潮喷娇喘18禁视频| 中文字幕色久视频| 国产又色又爽无遮挡免费看| 亚洲国产中文字幕在线视频| 亚洲视频免费观看视频| 人人妻人人添人人爽欧美一区卜| 岛国毛片在线播放| 午夜91福利影院| 俄罗斯特黄特色一大片| 色婷婷久久久亚洲欧美| 中亚洲国语对白在线视频| 国产精品98久久久久久宅男小说| 亚洲av国产av综合av卡| 精品国产超薄肉色丝袜足j| 在线观看免费日韩欧美大片| 少妇粗大呻吟视频| 三级毛片av免费| 国产精品影院久久| 欧美黄色淫秽网站| 国产淫语在线视频| 男女高潮啪啪啪动态图| 日本精品一区二区三区蜜桃| 国产免费视频播放在线视频| 丝袜美腿诱惑在线| 岛国在线观看网站| 视频在线观看一区二区三区| 亚洲精品在线美女| 人人妻人人澡人人看| 午夜免费成人在线视频| 精品欧美一区二区三区在线| 午夜福利一区二区在线看| 伦理电影免费视频| 老汉色av国产亚洲站长工具| 精品少妇内射三级| 一夜夜www| 欧美日韩av久久| 中国美女看黄片| 一区二区三区激情视频| 一区在线观看完整版| 18禁观看日本| 最新在线观看一区二区三区| 成年动漫av网址| 国产精品欧美亚洲77777| 热99久久久久精品小说推荐| xxxhd国产人妻xxx| 大陆偷拍与自拍| 欧美性长视频在线观看| 亚洲精品自拍成人| 老汉色∧v一级毛片| av福利片在线| 亚洲成人免费av在线播放| 欧美一级毛片孕妇| 日韩欧美一区二区三区在线观看 | 啦啦啦在线免费观看视频4| 国产精品.久久久| 久久精品国产综合久久久| 91av网站免费观看| av免费在线观看网站| 大香蕉久久成人网| 欧美+亚洲+日韩+国产| 十分钟在线观看高清视频www| 精品国产一区二区三区四区第35| 搡老乐熟女国产| 亚洲精品美女久久久久99蜜臀| 美女福利国产在线| 亚洲av日韩在线播放| 少妇 在线观看| 香蕉国产在线看| 好男人电影高清在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲精品久久成人aⅴ小说| 成人18禁高潮啪啪吃奶动态图| 99在线人妻在线中文字幕 | 极品少妇高潮喷水抽搐| 国产精品久久电影中文字幕 | 免费高清在线观看日韩| 国产精品国产高清国产av | av天堂久久9| 国产欧美亚洲国产| 夜夜骑夜夜射夜夜干| 亚洲国产欧美一区二区综合| 欧美亚洲日本最大视频资源| bbb黄色大片| 亚洲美女黄片视频| 高清在线国产一区| 亚洲精品美女久久av网站| 亚洲av成人一区二区三| 高清黄色对白视频在线免费看| 精品亚洲乱码少妇综合久久| 一区二区av电影网| 精品一区二区三卡| 青青草视频在线视频观看| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频| 免费看十八禁软件| 狠狠精品人妻久久久久久综合| 首页视频小说图片口味搜索| 热99久久久久精品小说推荐| 在线观看免费高清a一片| 一夜夜www| 麻豆乱淫一区二区| 无人区码免费观看不卡 | 天堂8中文在线网| videosex国产| 国产精品久久久久久精品电影小说| 国产黄频视频在线观看| 十八禁高潮呻吟视频| 亚洲欧洲日产国产| 久久久久精品人妻al黑| 欧美 亚洲 国产 日韩一| 亚洲av国产av综合av卡| 日韩制服丝袜自拍偷拍| 国产男女内射视频| 99在线人妻在线中文字幕 | 天天影视国产精品| 最新美女视频免费是黄的| 国产亚洲欧美在线一区二区| 岛国在线观看网站| 后天国语完整版免费观看| videos熟女内射| www.熟女人妻精品国产| 久久精品人人爽人人爽视色| 在线 av 中文字幕| 欧美日韩视频精品一区| 欧美av亚洲av综合av国产av| 黑人欧美特级aaaaaa片| 国产有黄有色有爽视频| 亚洲三区欧美一区| 亚洲免费av在线视频| 国产免费av片在线观看野外av| 窝窝影院91人妻| 免费观看av网站的网址| 在线观看www视频免费| 色婷婷av一区二区三区视频| 他把我摸到了高潮在线观看 | cao死你这个sao货| 亚洲精品国产区一区二| 久久久久久亚洲精品国产蜜桃av| 国产91精品成人一区二区三区 | 亚洲国产欧美一区二区综合| 狠狠精品人妻久久久久久综合| 成人三级做爰电影| 一二三四社区在线视频社区8| 高清欧美精品videossex| 超色免费av| 亚洲少妇的诱惑av| 制服人妻中文乱码| 成在线人永久免费视频| 亚洲五月色婷婷综合| 国产精品秋霞免费鲁丝片| 久热这里只有精品99| 亚洲精品av麻豆狂野| 久久久久久久精品吃奶| 亚洲国产欧美日韩在线播放| 亚洲综合色网址| 国产一区二区三区综合在线观看| www日本在线高清视频| 亚洲熟女精品中文字幕| 日韩视频在线欧美| 国产99久久九九免费精品| 亚洲国产精品一区二区三区在线| 悠悠久久av| 人妻久久中文字幕网| 久久久久精品人妻al黑| 国产精品国产高清国产av | 成人精品一区二区免费| 中亚洲国语对白在线视频| 日韩制服丝袜自拍偷拍| 亚洲欧洲精品一区二区精品久久久| 国产成人精品在线电影| 国产亚洲av高清不卡| 老司机午夜十八禁免费视频| 亚洲精品一二三| 一本综合久久免费| 制服诱惑二区| 国产不卡一卡二| 国产99久久九九免费精品| 亚洲一码二码三码区别大吗| 亚洲国产成人一精品久久久| 国产精品免费大片| 久久久国产精品麻豆| 一级黄色大片毛片| 国产97色在线日韩免费| 人人妻人人爽人人添夜夜欢视频| 欧美国产精品va在线观看不卡| 欧美人与性动交α欧美精品济南到| 日韩 欧美 亚洲 中文字幕| 国精品久久久久久国模美| 久久久久久久久久久久大奶| 99久久人妻综合| 制服人妻中文乱码| 免费女性裸体啪啪无遮挡网站| 国产无遮挡羞羞视频在线观看| 少妇的丰满在线观看| 99re6热这里在线精品视频| 在线 av 中文字幕| 国产国语露脸激情在线看| 国产精品久久久人人做人人爽| 欧美日韩国产mv在线观看视频| 超色免费av| 天天影视国产精品| 俄罗斯特黄特色一大片| 在线 av 中文字幕| 日本av免费视频播放| 91成人精品电影| 国产精品欧美亚洲77777| 高清视频免费观看一区二区| 中文字幕高清在线视频| 久久免费观看电影| 国产野战对白在线观看| 淫妇啪啪啪对白视频| 99国产精品一区二区蜜桃av | 精品久久久久久久毛片微露脸| 一级黄色大片毛片| 18禁裸乳无遮挡动漫免费视频| 99在线人妻在线中文字幕 | 国产不卡av网站在线观看| 夜夜骑夜夜射夜夜干| 午夜精品久久久久久毛片777| 亚洲精品久久成人aⅴ小说| netflix在线观看网站| 欧美成狂野欧美在线观看| 日本欧美视频一区| 亚洲中文日韩欧美视频| 久久性视频一级片| 国产精品香港三级国产av潘金莲| 在线观看免费高清a一片| 欧美日韩成人在线一区二区| 欧美在线黄色| 成人精品一区二区免费| 亚洲五月婷婷丁香| 国产又色又爽无遮挡免费看| 亚洲av成人一区二区三| 久久午夜亚洲精品久久| 久久久精品免费免费高清| 日韩免费高清中文字幕av| 精品熟女少妇八av免费久了| 热99久久久久精品小说推荐| 99国产精品一区二区三区| 99久久国产精品久久久| 亚洲男人天堂网一区| 欧美日韩精品网址|