• <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码专区亚洲av| 国产成人精品在线电影| 99国产综合亚洲精品| 国产精品av久久久久免费| 欧美精品一区二区大全| 熟女少妇亚洲综合色aaa.| 亚洲免费av在线视频| 桃花免费在线播放| av网站在线播放免费| 两个人看的免费小视频| 女人爽到高潮嗷嗷叫在线视频| 欧美日韩精品网址| 秋霞伦理黄片| 国产高清不卡午夜福利| 亚洲一卡2卡3卡4卡5卡精品中文| 久久97久久精品| 又黄又粗又硬又大视频| 免费少妇av软件| 成人亚洲欧美一区二区av| 一进一出抽搐动态| 国产精品电影一区二区三区| 精品国产超薄肉色丝袜足j| 99re在线观看精品视频| 国产一区二区在线av高清观看| 欧美激情极品国产一区二区三区| 香蕉丝袜av| 亚洲精品中文字幕在线视频| 麻豆av在线久日| 午夜福利在线观看吧| 亚洲黑人精品在线| 可以在线观看的亚洲视频| 无人区码免费观看不卡| 一卡2卡三卡四卡精品乱码亚洲| 91麻豆av在线| 婷婷丁香在线五月| 亚洲国产精品久久男人天堂| 亚洲第一电影网av| 欧美一级毛片孕妇| 亚洲第一欧美日韩一区二区三区| 午夜老司机福利片| 可以免费在线观看a视频的电影网站| 欧美中文日本在线观看视频| 99香蕉大伊视频| 女警被强在线播放| aaaaa片日本免费| 可以在线观看毛片的网站| 精品国产超薄肉色丝袜足j| 黄片小视频在线播放| 久久精品成人免费网站| 黑人欧美特级aaaaaa片| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站 | 亚洲第一青青草原| 久久午夜亚洲精品久久| 久久 成人 亚洲| 黑人欧美特级aaaaaa片| 国产精品1区2区在线观看.| 99re在线观看精品视频| 99riav亚洲国产免费| 欧美国产日韩亚洲一区| 国产精品98久久久久久宅男小说| 欧美午夜高清在线| 后天国语完整版免费观看| 人人澡人人妻人| 亚洲中文字幕日韩| 欧美黄色片欧美黄色片| 久久香蕉激情| av电影中文网址| 自线自在国产av| 亚洲无线在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 久久狼人影院| www.www免费av| 级片在线观看| 黄网站色视频无遮挡免费观看| 他把我摸到了高潮在线观看| 夜夜爽天天搞| 日韩欧美国产在线观看| 精品国产超薄肉色丝袜足j| 国产又爽黄色视频| 少妇被粗大的猛进出69影院| 亚洲avbb在线观看| 免费久久久久久久精品成人欧美视频| 亚洲天堂国产精品一区在线| 国产激情欧美一区二区| 亚洲一区中文字幕在线| 午夜视频精品福利| 一区二区日韩欧美中文字幕| 久久久精品欧美日韩精品| 丝袜人妻中文字幕| 夜夜夜夜夜久久久久| 欧美一级a爱片免费观看看 | 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久久久电影| 亚洲在线自拍视频| 一边摸一边抽搐一进一出视频| 乱人伦中国视频| 久久久国产欧美日韩av| 人妻丰满熟妇av一区二区三区| 国产精品一区二区免费欧美| 国产精品亚洲美女久久久| 少妇 在线观看| 国产真人三级小视频在线观看| 日韩欧美一区二区三区在线观看| 国产一区二区在线av高清观看| 久久人妻熟女aⅴ| 一区在线观看完整版| √禁漫天堂资源中文www| 看免费av毛片| 亚洲国产精品合色在线| 国产精品综合久久久久久久免费 | 国产精品99久久99久久久不卡| 国产午夜福利久久久久久| 亚洲精品中文字幕在线视频| 女人高潮潮喷娇喘18禁视频| avwww免费| 亚洲av成人一区二区三| 国产午夜福利久久久久久| 亚洲专区中文字幕在线| 欧美av亚洲av综合av国产av| 他把我摸到了高潮在线观看| 亚洲中文av在线| 91成年电影在线观看| 搡老妇女老女人老熟妇| av超薄肉色丝袜交足视频| 国产单亲对白刺激| 午夜久久久在线观看| www日本在线高清视频| 亚洲人成网站在线播放欧美日韩| 免费av毛片视频| 91成年电影在线观看| 国产国语露脸激情在线看| 一级毛片精品| 757午夜福利合集在线观看| 亚洲欧美精品综合一区二区三区| 宅男免费午夜| 99精品久久久久人妻精品| 国产片内射在线| 妹子高潮喷水视频| 久久久国产精品麻豆| 如日韩欧美国产精品一区二区三区| 国产激情欧美一区二区| 久久国产乱子伦精品免费另类| 淫秽高清视频在线观看| 午夜免费激情av| 一级毛片高清免费大全| 亚洲全国av大片| 天堂√8在线中文| 极品人妻少妇av视频| tocl精华| 国产区一区二久久| 午夜成年电影在线免费观看| 美女 人体艺术 gogo| 成人18禁高潮啪啪吃奶动态图| 人妻丰满熟妇av一区二区三区| 黑人欧美特级aaaaaa片| 久久国产亚洲av麻豆专区| 又黄又粗又硬又大视频| 亚洲三区欧美一区| 女人被躁到高潮嗷嗷叫费观| 亚洲伊人色综图| 午夜影院日韩av| 亚洲 国产 在线| 精品人妻在线不人妻| 黄色视频不卡| 亚洲av第一区精品v没综合| 在线观看免费日韩欧美大片| 无人区码免费观看不卡| 午夜成年电影在线免费观看| 欧美 亚洲 国产 日韩一| 激情视频va一区二区三区| av视频免费观看在线观看| 免费在线观看完整版高清| 国产av精品麻豆| 欧美乱妇无乱码| 日本在线视频免费播放| 日韩欧美国产在线观看| 中文字幕av电影在线播放| 老司机在亚洲福利影院| 999久久久国产精品视频| 久久国产乱子伦精品免费另类| 亚洲国产精品sss在线观看| 国产三级黄色录像| 亚洲国产精品成人综合色| 中文字幕人妻熟女乱码| 亚洲国产中文字幕在线视频| 国产日韩一区二区三区精品不卡| 精品卡一卡二卡四卡免费| 亚洲色图av天堂| 老熟妇乱子伦视频在线观看| 精品国产亚洲在线| 亚洲精品一区av在线观看| 午夜a级毛片| 99国产精品一区二区蜜桃av| 欧美乱妇无乱码| 欧美中文日本在线观看视频| 亚洲天堂国产精品一区在线| 久热爱精品视频在线9| 女人高潮潮喷娇喘18禁视频| 男女床上黄色一级片免费看| 免费观看人在逋| 午夜精品久久久久久毛片777| 在线av久久热| 香蕉久久夜色| 久久久久久人人人人人| 国产一区在线观看成人免费| 在线观看免费日韩欧美大片| 午夜日韩欧美国产| 少妇被粗大的猛进出69影院| 中文字幕人妻熟女乱码| 两性午夜刺激爽爽歪歪视频在线观看 | 国产亚洲精品一区二区www| 久久国产乱子伦精品免费另类| 满18在线观看网站| 久久久久久久久免费视频了| 老司机靠b影院| 亚洲精品中文字幕一二三四区| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 国产亚洲欧美在线一区二区| 人妻久久中文字幕网| 在线观看66精品国产| 亚洲av成人av| 国产成人欧美在线观看| 动漫黄色视频在线观看| 亚洲精品美女久久av网站| 久9热在线精品视频| 欧美最黄视频在线播放免费| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 色尼玛亚洲综合影院| 久久久久久免费高清国产稀缺| 国产三级黄色录像| 国产欧美日韩综合在线一区二区| 欧美精品亚洲一区二区| 美女高潮喷水抽搐中文字幕| 免费在线观看亚洲国产| 757午夜福利合集在线观看| 两个人看的免费小视频| 女性被躁到高潮视频| 一区二区日韩欧美中文字幕| 精品欧美国产一区二区三| 搡老熟女国产l中国老女人| 18禁国产床啪视频网站| 亚洲人成伊人成综合网2020| 欧美一区二区精品小视频在线| 满18在线观看网站| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产精品999在线| 国产欧美日韩一区二区精品| 亚洲自偷自拍图片 自拍| 欧美黑人精品巨大| 欧美人与性动交α欧美精品济南到| netflix在线观看网站| av在线播放免费不卡| 一进一出抽搐动态| 一级黄色大片毛片| 国产精品九九99| 黄色丝袜av网址大全| ponron亚洲| 亚洲欧美一区二区三区黑人| 免费一级毛片在线播放高清视频 | 69精品国产乱码久久久| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久久久99蜜臀| 99久久精品国产亚洲精品| 成人国产综合亚洲| 国产一区二区激情短视频| 亚洲国产精品久久男人天堂| 人人妻,人人澡人人爽秒播| 黄片大片在线免费观看| 欧美乱色亚洲激情| 成在线人永久免费视频| 精品第一国产精品| 亚洲视频免费观看视频| 亚洲人成电影免费在线| 亚洲国产欧美日韩在线播放| 人成视频在线观看免费观看| 日日夜夜操网爽| 国产真人三级小视频在线观看| 女警被强在线播放| 曰老女人黄片| 十八禁网站免费在线| 成人国语在线视频| 色哟哟哟哟哟哟| 亚洲av成人av| 国产精品亚洲av一区麻豆| 国产乱人伦免费视频| 免费高清视频大片| 中出人妻视频一区二区| 久久久久久人人人人人| 亚洲精华国产精华精| 亚洲国产欧美网| 精品电影一区二区在线| 黄色视频不卡| av网站免费在线观看视频| 国产精品美女特级片免费视频播放器 | 大陆偷拍与自拍| 男女下面插进去视频免费观看| 久9热在线精品视频| 人人妻人人澡人人看| 国产精品 欧美亚洲| 欧美乱色亚洲激情| 99国产极品粉嫩在线观看| 久久中文看片网| 久久精品成人免费网站| 女人精品久久久久毛片| 精品一区二区三区av网在线观看| 老司机在亚洲福利影院| 在线观看免费视频网站a站| 性色av乱码一区二区三区2| 99精品久久久久人妻精品| 欧美久久黑人一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产1区2区3区精品| 国产精品亚洲一级av第二区| 女同久久另类99精品国产91| 999久久久精品免费观看国产| 欧美激情极品国产一区二区三区| 狠狠狠狠99中文字幕| 日本一区二区免费在线视频| 久久精品亚洲精品国产色婷小说| 97碰自拍视频| 午夜福利免费观看在线| 欧美黄色淫秽网站| 手机成人av网站| 黄色女人牲交| 国产精品久久视频播放| 99久久综合精品五月天人人| 国产亚洲精品综合一区在线观看 | 中国美女看黄片| 法律面前人人平等表现在哪些方面| 色播在线永久视频| www.自偷自拍.com| 中文字幕人妻丝袜一区二区| 母亲3免费完整高清在线观看| 亚洲七黄色美女视频| 91国产中文字幕| 久久伊人香网站| 一级毛片女人18水好多| 欧美成人午夜精品| 宅男免费午夜| 久久中文字幕人妻熟女| 18禁黄网站禁片午夜丰满| 丝袜美足系列| 99国产极品粉嫩在线观看| 99国产综合亚洲精品| 可以免费在线观看a视频的电影网站| 又黄又粗又硬又大视频| 97人妻精品一区二区三区麻豆 | 欧美日韩亚洲国产一区二区在线观看| 热re99久久国产66热| 一区二区日韩欧美中文字幕| 两个人视频免费观看高清| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 国产不卡一卡二| 国产精华一区二区三区| 可以在线观看毛片的网站| 午夜免费成人在线视频| 国产成人精品无人区| 久久久久精品国产欧美久久久| 波多野结衣高清无吗| 91av网站免费观看| 美女扒开内裤让男人捅视频| netflix在线观看网站| 国产成人精品久久二区二区免费| 久久婷婷成人综合色麻豆| 男女之事视频高清在线观看| 50天的宝宝边吃奶边哭怎么回事| 侵犯人妻中文字幕一二三四区| 一夜夜www| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看影片大全网站| 亚洲狠狠婷婷综合久久图片| 无人区码免费观看不卡| 老司机深夜福利视频在线观看| 久久婷婷人人爽人人干人人爱 | 国产aⅴ精品一区二区三区波| 一区二区三区高清视频在线| 丝袜人妻中文字幕| 黄色丝袜av网址大全| 搡老熟女国产l中国老女人| 午夜福利一区二区在线看| 久久久久国内视频| 精品第一国产精品| 成人国产一区最新在线观看| 欧美成人一区二区免费高清观看 | e午夜精品久久久久久久| 国产免费男女视频| 亚洲精品美女久久久久99蜜臀| 操美女的视频在线观看| 9热在线视频观看99| 99热只有精品国产| 女人精品久久久久毛片| 757午夜福利合集在线观看| 亚洲人成77777在线视频| 九色国产91popny在线| 91精品三级在线观看| 伦理电影免费视频| 波多野结衣高清无吗| 午夜老司机福利片| 成人国产一区最新在线观看| 亚洲成a人片在线一区二区| 日韩免费av在线播放| 国产精品 国内视频| 国产精品免费视频内射| 老鸭窝网址在线观看| 精品国产乱子伦一区二区三区| 久久天堂一区二区三区四区| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品sss在线观看| 亚洲人成电影免费在线| 国产亚洲精品第一综合不卡| 在线观看免费午夜福利视频| 欧美日本亚洲视频在线播放| 国产麻豆69| 夜夜看夜夜爽夜夜摸| 国产av又大| 精品不卡国产一区二区三区| 亚洲精品国产区一区二| 又黄又粗又硬又大视频| 香蕉丝袜av| 老熟妇乱子伦视频在线观看| 亚洲精品久久国产高清桃花| 色婷婷久久久亚洲欧美| 午夜免费激情av| 搡老熟女国产l中国老女人| 亚洲激情在线av| 91老司机精品| 黄片小视频在线播放| 可以在线观看毛片的网站| 国产不卡一卡二| 好男人电影高清在线观看| 大陆偷拍与自拍| 色综合欧美亚洲国产小说| 欧美一区二区精品小视频在线| 精品免费久久久久久久清纯| 老司机福利观看| 国产一级毛片七仙女欲春2 | 中文字幕色久视频| 一区二区三区激情视频| 国产蜜桃级精品一区二区三区| 好男人在线观看高清免费视频 | 久99久视频精品免费| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 国产国语露脸激情在线看| 国产高清视频在线播放一区| 国产欧美日韩一区二区精品| 12—13女人毛片做爰片一| 免费在线观看亚洲国产| 好男人电影高清在线观看| 欧美色视频一区免费| 一边摸一边做爽爽视频免费| 99riav亚洲国产免费| 9热在线视频观看99| 男人的好看免费观看在线视频 | 欧美久久黑人一区二区| 亚洲激情在线av| 欧洲精品卡2卡3卡4卡5卡区| 精品一品国产午夜福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人国产一区在线观看| 精品第一国产精品| 国产精品免费一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 男女下面进入的视频免费午夜 | 久久性视频一级片| 一级a爱视频在线免费观看| 国产精品久久久久久人妻精品电影| 熟妇人妻久久中文字幕3abv| 十八禁人妻一区二区| 亚洲熟女毛片儿| www.www免费av| 国产精品一区二区精品视频观看| 久久九九热精品免费| 亚洲精品在线美女| 亚洲色图 男人天堂 中文字幕| 日本黄色视频三级网站网址| 亚洲色图 男人天堂 中文字幕| 涩涩av久久男人的天堂| 久久久国产成人精品二区| 大型黄色视频在线免费观看| 免费在线观看亚洲国产| 亚洲美女黄片视频| 99久久综合精品五月天人人| 午夜福利欧美成人| 成人18禁在线播放| 亚洲成人精品中文字幕电影| 91国产中文字幕| 亚洲精品在线美女| 美女 人体艺术 gogo| 女人高潮潮喷娇喘18禁视频| 性欧美人与动物交配| 成人av一区二区三区在线看| 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| www国产在线视频色| 欧美av亚洲av综合av国产av| 国产在线精品亚洲第一网站| 日韩三级视频一区二区三区| 精品欧美一区二区三区在线| 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 成人av一区二区三区在线看| 日本一区二区免费在线视频| 国产男靠女视频免费网站| 久久久久国内视频| 两性午夜刺激爽爽歪歪视频在线观看 | 丁香欧美五月| 后天国语完整版免费观看| 岛国在线观看网站| 亚洲中文av在线| 国产亚洲欧美精品永久| 日本免费一区二区三区高清不卡 | 一级毛片高清免费大全| 国产一区在线观看成人免费| 最近最新中文字幕大全免费视频| 色综合欧美亚洲国产小说| 精品国产亚洲在线| av网站免费在线观看视频| 搡老熟女国产l中国老女人| 国产成年人精品一区二区| 啦啦啦观看免费观看视频高清 | avwww免费| 男人操女人黄网站| 一进一出抽搐gif免费好疼| 欧美乱妇无乱码| 神马国产精品三级电影在线观看 | 久久精品人人爽人人爽视色| 国产精品野战在线观看| 久久久久久久午夜电影| 久热爱精品视频在线9| 亚洲成av人片免费观看| 日本欧美视频一区| 大码成人一级视频| 丰满的人妻完整版| 中文字幕高清在线视频| 18禁美女被吸乳视频| 身体一侧抽搐| 操出白浆在线播放| 自线自在国产av| 国产亚洲欧美98| 午夜免费观看网址| 人妻丰满熟妇av一区二区三区| 狠狠狠狠99中文字幕| 侵犯人妻中文字幕一二三四区| 老司机午夜十八禁免费视频| 精品国产一区二区三区四区第35| 精品无人区乱码1区二区| 一进一出好大好爽视频| 九色国产91popny在线| 国产亚洲精品久久久久5区| 欧美不卡视频在线免费观看 | 免费少妇av软件| 50天的宝宝边吃奶边哭怎么回事| 久久久久久久午夜电影| 精品少妇一区二区三区视频日本电影| 真人一进一出gif抽搐免费| 亚洲国产看品久久| 黑人欧美特级aaaaaa片| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| 国产精品精品国产色婷婷| 动漫黄色视频在线观看| 精品国产一区二区久久| 亚洲成人国产一区在线观看| 国产1区2区3区精品| 色综合站精品国产| 久久精品亚洲熟妇少妇任你| 日韩av在线大香蕉| 此物有八面人人有两片| 人妻久久中文字幕网| videosex国产| 超碰成人久久| 亚洲人成电影免费在线| 在线视频色国产色| 日日爽夜夜爽网站| 亚洲av第一区精品v没综合| 老熟妇仑乱视频hdxx| 国产成人欧美在线观看| 久久香蕉国产精品| 亚洲av片天天在线观看| 久久九九热精品免费| 午夜免费成人在线视频| 亚洲人成伊人成综合网2020| 一级a爱片免费观看的视频| 国产1区2区3区精品| 狠狠狠狠99中文字幕| 韩国精品一区二区三区| 亚洲av电影不卡..在线观看| 国产精品爽爽va在线观看网站 | 精品国产国语对白av| 夜夜爽天天搞| 亚洲国产高清在线一区二区三 | av在线天堂中文字幕| 亚洲精品国产区一区二| 亚洲精品久久国产高清桃花| 黄色 视频免费看| 国产av又大| 国产精品一区二区免费欧美| 免费在线观看影片大全网站| 大型av网站在线播放| 亚洲一区二区三区不卡视频| 国产成人精品久久二区二区91| 欧美最黄视频在线播放免费| 国产主播在线观看一区二区| 级片在线观看| 国产高清激情床上av| 精品午夜福利视频在线观看一区|