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

    一次西南印度洋熱帶氣旋過程發(fā)展機制的數(shù)值研究*

    2017-12-14 05:34:10郭敬天曹永正劉桂艷
    海洋與湖沼 2017年4期
    關(guān)鍵詞:潛熱氣旋印度洋

    王 彬 郭敬天① 曹永正 張 薇 劉桂艷

    (1.國家海洋局北海預報中心 青島 266100; 2.山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室 青島 266100)

    一次西南印度洋熱帶氣旋過程發(fā)展機制的數(shù)值研究*

    王 彬1,2郭敬天1,2①曹永正1張 薇1劉桂艷2

    (1.國家海洋局北海預報中心 青島 266100; 2.山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室 青島 266100)

    本文利用盡可能多的觀測資料和WRF-3.4.1模式(Weather Research and Forecasting Model)對2012年1月19日至28日發(fā)生在西南印度洋上空的1次強熱帶氣旋進行了研究, 并分析其時空結(jié)構(gòu)和發(fā)展機制。對該熱帶氣旋的移動路徑、強度及內(nèi)部結(jié)構(gòu)的數(shù)值模擬結(jié)果與實際符合較好。分析表明, 感熱和凝結(jié)潛熱貫串于熱帶氣旋發(fā)展的整個過程, 其中感熱對氣旋發(fā)展的影響較弱, 凝結(jié)潛熱是氣旋發(fā)展的主要能量來源, CISK機制可解釋該熱帶氣旋的發(fā)展過程。

    西南印度洋; 熱帶氣旋; WRF模式; 發(fā)展機制

    熱帶氣旋是發(fā)生在較低緯度海洋上并具有暖心結(jié)構(gòu)的大型渦旋, 是一種強災(zāi)害性天氣系統(tǒng), 它引發(fā)的大風和大浪天氣影響范圍廣, 嚴重威脅海上船舶的航行和作業(yè)安全。探討不同海區(qū)上空發(fā)生的熱帶氣旋的特點和發(fā)展機制, 能加強對熱帶氣旋的認識, 為海洋災(zāi)害性天氣預報提供參考。

    目前熱帶氣旋研究多集中在北半球, 經(jīng)過半個多世紀的探索, 學者們已經(jīng)在熱帶氣旋生成所需的氣候背景場、擾動源以及生成和發(fā)展機制等方面有了一定了解(張慶紅等, 2008)。我國作為世界上受熱帶氣旋影響最嚴重的國家之一, 過去10年在熱帶氣旋運動突變、結(jié)構(gòu)和強度突變、臺風暴雨的突然增幅、路徑預報方法研究等方面都取得了新的進展(陳聯(lián)壽等, 2001)。相比而言, 受觀測資料缺乏、地理環(huán)境復雜等條件的限制, 針對南大洋尤其是西南印度洋海域熱帶氣旋的研究較少。從2012年開始, 我國將對西南印度洋開展連續(xù)10年的礦產(chǎn)勘探工作, 多種調(diào)查設(shè)備均對海況有較高的要求, 因此對西南印度洋天氣系統(tǒng)特別是熱帶氣旋發(fā)展過程開展研究十分重要。

    熱帶氣旋的發(fā)展因擾動能量來源不同, 相應(yīng)地存在著兩種不同假設(shè): 第二類條件不穩(wěn)定理論和海-氣相互作用理論。Charney等(1964)提出了第二類條件不穩(wěn)定理論CISK(Conditional Instability of Second Kind)機制, 該理論認為積云對流釋放的凝結(jié)潛熱是熱帶氣旋發(fā)展的能量來源。Emanuel(1986)提出了海-氣相互作用理論 WISHE(Wind-Induced Surface Heat Exchange), 它把熱帶氣旋的發(fā)展過程總結(jié)為一個初始擾動轉(zhuǎn)變?yōu)楹Q蟊砻骒释亢惋L場之間正反饋的過程。應(yīng)該指出, 兩種理論解釋熱帶氣旋的生成和發(fā)展都有其局限性, 前者沒有清楚說明臺風中心不在對流云區(qū)而是在云外晴空區(qū), 后者忽略了各云塊間下沉時的絕熱增溫使臺風中心增暖的作用, 前人對此也有很多的探討和研究(沈如金等, 1982; 劉裕祿等,2009)。

    近年來隨著高分辨率氣象資料的出現(xiàn), 熱帶氣旋的數(shù)值模擬工作也有了較大的發(fā)展。袁金南等(2003)采用廣州有限區(qū)域數(shù)值預報模式, 研究了對流凝結(jié)潛熱對 0214號臺風“黃蜂”登陸的影響, 結(jié)果表明, 模式中有無對流凝結(jié)潛熱加熱對臺風移動路徑和降水等均有明顯的影響。曹鈺等(2012)利用WRF(Weather Research and Forecasting Model)模式較成功地模擬了2009年第8號臺風“莫拉克”登陸的過程, 通過考慮對流凝結(jié)潛熱加熱, 探討了對非均勻飽和大氣中非地轉(zhuǎn)濕Q矢量的改進。陳聯(lián)壽(2010)認為, 以全球模式為背景嵌套有限區(qū)中尺度模式, 并結(jié)合先進的三維變分或四維變分同化系統(tǒng)來開展臺風預報和研究已占主導地位。

    南半球目前關(guān)于熱帶氣旋發(fā)展機制的研究相對較少, 主要的研究內(nèi)容集中在熱帶氣旋的活動路徑、季節(jié)變化等。楊亞新(2005)根據(jù)世界氣象組織(WMO)和全球熱帶氣旋業(yè)務(wù)網(wǎng)站收集的1968—1990年資料,闡述了全球熱帶氣旋的情況, 其中西南印度洋平均每年發(fā)生10.4個熱帶氣旋。冷梅等(2002)利用聯(lián)合臺風警報中心(JTWC)提供的 2000—2002年氣候資料,統(tǒng)計了南半球熱帶氣旋的活動, 結(jié)果表明南印度洋熱帶氣旋發(fā)生的頻率和強度大, 生命期長, 初始位置的平均緯度大約在15°S, 基本上向西偏南移動。

    基于上述西南印度洋熱帶氣旋的研究現(xiàn)狀, 本文選取了2012年1月的一次西南印度洋熱帶氣旋過程, 首先利用衛(wèi)星和再分析資料進行分析, 然后使用WRF模式模擬氣旋的發(fā)展過程, 最后通過分析感熱通量和凝結(jié)潛熱, 探討影響氣旋發(fā)展的能量來源和機制。

    1 資料分析

    1.1 統(tǒng)計分析

    利用FNL(Final Operational Global Analysis)再分析資料, 對2003—2012年10年期間影響西南印度洋的熱帶氣旋進行統(tǒng)計后發(fā)現(xiàn)(表 1): 每年影響西南印度洋的熱帶氣旋在 5—25個之間, 平均每年13.8個;10年間熱帶氣旋影響的頻率呈上升趨勢, 從2003年到 2006年的 5—12個上升至 2007年到 2012年的13—25個; 熱帶氣旋影響較多的月份為4—10月, 目前大洋科考任務(wù)多安排在 11月至翌年 4月間, 也是為了盡可能規(guī)避熱帶氣旋的影響。該區(qū)域熱帶氣旋的發(fā)展比較旺盛, 中心氣壓都低于 1000hPa, 個別氣旋甚至達到930hPa左右(如本文研究的熱帶氣旋中心最低氣壓為 936hPa), 它們主要在馬達加斯加島附近海域生成, 然后直接南下, 對作業(yè)海區(qū)的影響較強。

    表1 2003—2012年影響西南印度洋的熱帶氣旋數(shù)量Tab.1 The amount of tropical cyclone that influenced the southwest Indian Ocean in 2003—2012

    1.2 研究個例

    表1表明, 熱帶氣旋對西南印度洋的影響在增加,本文選取2012年1月19日至28日發(fā)生在西南印度洋海域的一次強熱帶氣旋過程進行研究。利用MODIS(Moderate Resolution Imaging Spectoradiometer)衛(wèi)星可見光云圖, 可以清晰地看到本次熱帶氣旋的發(fā)展過程(圖1)。氣旋于19日06UTC (Universal Time Coordinated)在馬達加斯加島西側(cè)的莫桑比克海峽生成, 其中心位于(17.3°S, 40.7°E)附近。在 20 日的云圖上(圖 1a), 螺旋云帶中的氣旋式彎曲已非常明顯, 但其中心無云眼區(qū)還沒有形成, 至23日00UTC(圖1b),衛(wèi)星云圖上已出現(xiàn)明顯的臺風眼結(jié)構(gòu), 表明氣旋的發(fā)展進入成熟階段, 此階段一直持續(xù)到 25日00UTC(圖 1c), 其眼區(qū)結(jié)構(gòu)明顯, 并逐漸向南移動。此后氣旋進入消亡階段, 27日00UTC(圖1d)眼區(qū)已不能清楚分辨, 28日00UTC, 螺旋狀云系的結(jié)構(gòu)開始變得松散。

    圖1 西南印度洋海域2012年1月20—27日MODIS紅外衛(wèi)星云圖Fig.1 The MODIS satellite infrared images of the southwest Indian Ocean, 20 to 27, January, 2012

    利用 FNL再分析數(shù)據(jù), 重點考察熱帶氣旋在成熟階段的空間結(jié)構(gòu)。在24日12UTC的海平面氣壓場上(圖 2a), 熱帶氣旋的中心位于(21.5°S, 39°E)附近,中心氣壓值約為 984hPa。為了考察熱帶氣旋的垂直結(jié)構(gòu), 沿通過氣旋中心的A-B線(其位置見圖2a)做垂直剖面分析。圖2b為沿A-B線的位溫和垂直速度分布圖, 可以看到垂直速度在21°S附近的600—200hPa之間出現(xiàn)小值區(qū), 而且兩側(cè)伴有強烈的垂直上升運動, 表明熱帶氣旋的中心眼區(qū)位于該位置。等位溫線的分布在熱帶氣旋中心附近出現(xiàn)了一個“漏斗”狀結(jié)構(gòu), 這種位溫分布表明氣旋具有“暖心”結(jié)構(gòu)。

    圖2 2012年1月24日12UTC空間結(jié)構(gòu)分析Fig.2 Structured analysis at 12UTC on 24 January 2012

    2 數(shù)值試驗

    2.1 數(shù)值模式

    本文選用 WRF模式(3.4.1版本)作為研究工具,它是由美國國家大氣研究中心(NCAR)等機構(gòu)聯(lián)合開發(fā)的新一代中尺度非靜力預報模式, 經(jīng)過不斷地改進和更新, 目前該模式內(nèi)部參數(shù)化方案豐富, 對物理過程的刻畫也較為細致。

    2.2 試驗設(shè)置

    模擬區(qū)域設(shè)置如圖 3所示, 模擬時間為 12年 1月22日00UTC至1月28日00UTC, 模式的具體參數(shù)設(shè)置見表2。

    圖3 WRF模擬的計算區(qū)域Fig.3 Simulation domain in WRF

    試驗采用 FNL再分析數(shù)據(jù)為初始場提供背景數(shù)據(jù)與時變側(cè)邊界。Fu(2001)指出高質(zhì)量的初始場對氣旋的準確預報有至關(guān)重要的作用, 因此本文采用了高山紅等(2010a, b)設(shè)計和發(fā)展的循環(huán)三維變分同化手段, 同化了初始時刻之前12h的觀測數(shù)據(jù)。同化的觀測數(shù)據(jù)包括固定站點的探空(12h/次)和地面站(3h/次)等 GTS(Global Telecommunication System)數(shù)據(jù)、衛(wèi)星反演的溫度廓線 AIRS(Atmospheric Infrared Sounder)、大氣可降水量 SSMI(Special Sensor Microwave Imager)以及BUFR格式的船舶、飛機報等觀測數(shù)據(jù)。

    表2 WRF模式設(shè)置Tab.2 Specifications of WRF model

    在利用循環(huán) 3DVAR同化數(shù)據(jù)之前, 需要給定背景誤差協(xié)方差, 由于WRF模式提供的CV3背景誤差的模擬效果不理想, 所以本文選用了 NMC(National Meteorological Center)方法生成的 CV5背景誤差(Parrish and Derber, 1992)。利用FNL數(shù)據(jù), 以初始時刻為中心前后各7天, 針對本次過程作每天2次(起點分別為00UTC和12UTC)的24h后報, 然后利用后報結(jié)果統(tǒng)計出CV5背景誤差。

    3 結(jié)果與分析

    3.1 模擬結(jié)果初步分析

    首先對模擬的結(jié)果進行驗證, 在圖 4a中可以看出模擬和觀測的路徑比較接近, 說明模式較好地模擬出了氣旋的移動。圖4b為熱帶氣旋中心氣壓變化圖, 紅線為 WRF模擬結(jié)果, 黑線為實測資料。對比來看, WRF模擬的中心氣壓在氣旋發(fā)展前期較實測偏高, 后期略有偏低, 但曲線的走向基本一致, 都是在24日12UTC達到最小值。

    下面分析模擬的熱帶氣旋在發(fā)展最成熟時期(24日12UTC)的結(jié)構(gòu)。通過對比模擬的海平面氣壓場(圖4c)和 FNL再分析資料的結(jié)果(圖 2a), 可以發(fā)現(xiàn)兩者的形勢十分接近, 但模擬得到的中心最低氣壓(928hPa)更接近實測值(936hPa)。圖 4d是沿圖 4c中C-D 線(39.5°E)的位溫和垂直速度的剖面圖, 從中看到等位溫線在(20.5°S, 39.5°E)附近有明顯的向下彎曲,表明模擬結(jié)果很好地反映了熱帶氣旋的“暖心”結(jié)構(gòu);同時氣旋中心的兩側(cè)從900hPa到200hPa的垂直速度中心, 說明存在較強的對流運動, 模擬的氣旋中心北側(cè)垂直速度要遠高于南側(cè), 這與再分析資料的結(jié)果(圖2b)有差異。

    圖4 模擬的熱帶氣旋路徑、強度和結(jié)構(gòu)分析Fig.4 Simulated tropical cyclone tracks, central pressure and structure

    從以上分析來看, 本次模擬較為合理地反映了氣旋的移動路徑、強度和內(nèi)部結(jié)構(gòu), 使用該模擬結(jié)果進一步分析是可信的。

    3.2 氣旋發(fā)展機制

    從天氣形勢來看, 熱帶氣旋在發(fā)生發(fā)展時沒有外來能量如冷空氣或者低空急流等的注入, 結(jié)合上文提到的CISK和WHISE理論, 通過計算熱帶氣旋發(fā)展過程中的海面感熱通量和凝結(jié)潛熱, 分析影響熱帶氣旋發(fā)展的能量來源, 以探討其發(fā)展的機理。

    洋面上感熱通量計算公式如下(丁一匯, 1993):Va是海表面10m高度上的風向量,Tw是海表溫度,Ta是海表面 10m 高度上的氣溫, 此時Fh單位為W/m2。這里給出了氣旋在23日12UTC至25日12UTC的感熱通量變化(圖5)。根據(jù)圖4b中觀測的氣旋中心氣壓強度, 該熱帶氣旋在23—25日經(jīng)歷了一個發(fā)展、成熟然后逐漸衰弱的過程。但從圖5中發(fā)現(xiàn), 感熱通量的最大值(60W/m2)出現(xiàn)在23日12UTC, 此后持續(xù)減小, 至25日最大感熱通量已減小為30W/m2。分析原因, 是熱帶氣旋在發(fā)展增強的階段, 海面風雖然繼續(xù)增大, 但因海表氣溫差減小值大于風速增大值, 使感熱通量減小, 感熱作用持續(xù)減弱。綜上所述, 感熱作用雖然貫串熱帶氣旋發(fā)展、成熟、減弱整個過程,但其最大值出現(xiàn)的時次與氣旋最大強度還不一致,并不是熱帶氣旋發(fā)展和增強的主要原因。

    凝結(jié)潛熱的計算需要大氣條件是不穩(wěn)定的, 水汽是輻合的。這里采用郭曉嵐(Kuo, 1965, 1974)所提出的積云對流參數(shù)化方法計算凝結(jié)潛熱Hc(m2/s3):

    其中I為單位空氣柱總的水汽流入率:

    圖5 2012年1月23日12UTC至25日12UTC熱帶氣旋感熱通量(W/m2)變化Fig.5 Sensible heat flux (W/m2) of tropical cyclone from 12UTC on 23 to 12UTC on 25 January 2012

    式中Pt為云頂高度,Pb為云底高度,T為云外環(huán)境溫度,TS為云中溫度。

    這里計算了包含整個氣旋(38°—42°E, 20°—26°S)區(qū)域大氣氣柱內(nèi)的凝結(jié)潛熱釋放。從潛熱隨時間的演變圖(圖 6)中可以看到, 熱帶氣旋在發(fā)展過程中伴隨著大量的對流凝結(jié)潛熱釋放, 并且在25日00UTC達到最大值, 為 1938m2/s3, 這比觀測到的氣旋最大強度的時間(24日12UTC)滯后了12小時, 但凝結(jié)潛熱釋放與氣旋強度變化的趨勢(圖4b)是一致的。25日以后當熱帶氣旋逐漸減弱時,Hc值略有減小, 但仍然保持在 1756m2/s3。氣旋在 27日 00UTC時Hc減小到1212m2/s3。以上說明: 對流凝結(jié)潛熱, 其大小與熱帶氣旋強度呈正相關(guān)關(guān)系, 在熱帶氣旋發(fā)展成熟階段,對流凝結(jié)潛熱釋放顯著增大, 在逐漸衰弱階段, 凝結(jié)潛熱仍然保持較大值。由此推論, 凝結(jié)潛熱是熱帶氣旋維持和發(fā)展的主要能量來源。

    圖6 2012年22日00UTC—27日00UTC氣柱內(nèi)總的對流凝結(jié)潛熱(m2/s3)的變化Fig.6 Distribution of total convective condensation latent heat (m2/s3) from 00UTC on 22 to 00UTC on 27 January 2012

    為進一步驗證上述分析, 給出了沿氣旋中心的凝結(jié)潛熱和水汽的垂直剖面圖(圖7)。24日12UTC(圖7a), 熱帶氣旋凝結(jié)潛熱最大值為10 m2/s3, 此時氣旋中心兩側(cè)水汽有明顯的向上輻合, 對流發(fā)展旺盛。而在26日12UTC(圖7b), 熱帶氣旋凝結(jié)潛熱最大值為7m2/s3, 此時氣旋內(nèi)部水汽分布逐漸均勻, 對流強度也明顯減弱, 這進一步說明了凝結(jié)潛熱釋放對氣旋發(fā)展作用巨大, 它不僅是影響垂直上升運動的關(guān)鍵因子, 還影響水汽的分布。

    綜上分析, 可以用CISK理論解釋本次熱帶氣旋發(fā)展的機制。19日, 在莫桑比克海峽上空有一弱小的熱帶低壓生成, 通過高溫洋面的摩擦作用, 大量熱帶潮濕空氣輻合到氣旋中心, 并產(chǎn)生上升運動, 暖空氣在上空凝結(jié)釋放出大量潛熱, 使低壓中心上空的大氣溫度升高, 形成暖心。同時高層等壓面抬高形成輻散流出, 使地面氣壓降低、環(huán)流加強。環(huán)流變強進一步使低空暖濕空氣向中心摩擦輻合加強, 更多的水汽向中心集中, 對流更旺盛, 凝結(jié)潛熱釋放更多, 中心變得更暖, 地面低壓更低。如此循環(huán), 促進氣旋的發(fā)展。25日以后, 氣旋移至中緯度洋面, 凝結(jié)潛熱提供的能量不足以維持氣旋的發(fā)展, 氣旋逐漸減弱。

    圖7 經(jīng)過氣旋中心的凝結(jié)潛熱Hc(m2/s3, 藍線)和水汽(g/kg, 黑線)剖面圖Fig.7 The vertical structure of condensation latent heat(m2/s3, blue line)and water vapor(g/kg, black line) along the center of tropical cyclone

    4 結(jié)論與展望

    本文利用觀測分析和數(shù)值模擬的手段, 對 2012年 1月的一次西南印度洋熱帶氣旋的發(fā)展過程和機制進行了研究, 主要結(jié)論如下:

    (1) 利用再分析資料對西南印度洋熱帶氣旋進行統(tǒng)計發(fā)現(xiàn), 2003—2012年間熱帶氣旋影響的頻率呈上升趨勢, 從2003到2006年的5—12個上升至2007到2012年的13—25個;

    (2) 該熱帶氣旋成熟時期在衛(wèi)星云圖上出現(xiàn)了螺旋云帶和臺風眼結(jié)構(gòu)。使用WRF模式對該氣旋的模擬是成功的, 模擬的路徑、強度和空間結(jié)構(gòu), 都合理地再現(xiàn)了氣旋的特征, 為進一步研究其發(fā)展機制提供了必要條件;

    (3) 在熱帶氣旋發(fā)展過程中, 感熱和凝結(jié)潛熱作用始終貫串于整個過程, 其中感熱對氣旋發(fā)展的作用比較微弱;

    (4) 凝結(jié)潛熱的變化與氣旋強度呈正相關(guān), 是氣旋發(fā)展的主要能量來源, CISK機制可以解釋本次氣旋發(fā)展的過程。

    針對此次個例, 今后工作展望如下:

    (1) 本文重點探究了熱帶氣旋形成以后的發(fā)展機制, 未來可以研究它前期的形成過程;

    (2) 莫桑比克海峽的存在會對熱帶氣旋的移動路徑產(chǎn)生影響, 可以考慮設(shè)計地形的敏感性試驗進行分析;

    (3) 此次個例是一次典型的影響西南印度洋的熱帶氣旋過程, 但缺乏普適性, 今后可以采用統(tǒng)計學的方法, 對影響西南印度洋的熱帶氣旋的特征進行總結(jié)。

    致謝 中國海洋大學高山紅教授對本文寫作以及循環(huán)三維變分和繪圖方面提出了寶貴意見, 在此表示感謝。

    丁一匯, 1993.天氣動力學中的診斷分析方法.北京: 中國科技大學研究生院, 118—120

    劉裕祿, 方祥生, 金飛勝等, 2009.臺風鳳凰形成發(fā)展過程中對流凝結(jié)潛熱和感熱的作用.氣象, 35(12): 51—57

    楊亞新, 2005.全球熱帶氣旋概述.世界海運, 28(2): 3—6

    冷 梅, 王 輝, 2002.船舶接收北大西洋熱帶氣旋信息的途徑.見: 中國航海學會海洋船舶防臺學術(shù)會議論文集.大連: 中國航海學會

    沈如金, 張寶嚴, 1982.凝結(jié)潛熱加熱對臺風降水分布的影響.大氣科學, 6(3): 249—257

    張慶紅, 郭春蕊, 2008.熱帶氣旋生成機制的研究進展.海洋學報, 30(4): 1—11

    陳聯(lián)壽, 2010.熱帶氣象災(zāi)害及其研究進展.氣象, 36(7): 101—110

    陳聯(lián)壽, 孟智勇, 2001.我國熱帶氣旋研究十年進展.大氣科學, 25(3): 420—432

    袁金南, 萬齊林, 2003.島嶼地形和對流凝結(jié)潛熱對登陸臺風“黃蜂”影響的數(shù)值研究.熱帶氣象學報, 19(S1): 81—87

    高山紅, 齊伊玲, 張守寶等, 2010.利用循環(huán)3DVAR改進黃海海霧數(shù)值模擬初始場Ⅰ: WRF數(shù)值試驗.中國海洋大學學報, 40(10): 1—9

    高山紅, 張守寶, 齊伊玲等, 2010.利用循環(huán)3DVAR改進黃海海霧數(shù)值模擬初始場Ⅱ: RAMS數(shù)值試驗.中國海洋大學學報, 40(11): 1—10, 18

    曹 鈺, 苗春生, 岳彩軍等, 2012.引入對流凝結(jié)潛熱作用對非均勻飽和大氣中非地轉(zhuǎn)濕 Q矢量的改進研究.高原氣象, 31(1): 76—86

    Charney J G, Eliassen A, 1964.On the growth of the hurricane depression.J Atmos Sci, 21(1): 68—75

    Chen F, Mitchell K, Schaake Jet al, 1996.Modeling of land surface evaporation by four schemes and comparison with FIFE observations.J Geophys Res, 101(D3): 7251—7266

    Dudhia J, 1996.A multi-layer soil temperature model for MM 5.In: Sixth Annual PSU/NCAR Mesoscale Model Users’Workshop.Boulder Colorado: National Center for Atmospheric Research, 49—51

    Emanuel K A, 1986.An air-sea interaction theory for tropical cyclones.PartΙ: steady-state maintenance.J Atmos Sci,43(6): 585—605

    Fu G, 2001.Polar Lows: Intense Cyclones in Winter.Beijing:China Meteorological Press, 218

    Hong S Y, Dudhia J, Chen S H, 2004.A revised approach to ice microphysical processes for the bulk parameterization of clouds and precipitation.Mon Wea Rev, 132(1): 103—120

    Hong S Y, Noh Y, Dudhia J, 2006.A new vertical diffusion package with an explicit treatment of entrainment processes.Mon Wea Rev, 134(9): 2318—2341

    Kain J S, Fritsch J M, 1990.A one-dimensional entraining/detraining plume model and its application in convective parameterization.J Atmos Sci, 47(23): 2784—2802

    Kuo H L, 1965.On formation and intensification of tropical cyclones through latent heat release by cumulus convection.J Atmos Sci, 22(1): 40—63

    Kuo H L, 1974.Further studies of the Parameterization of the influence of cumulus convection on large-scale flow.J Atmos Sci, 31(5): 1232—1240

    Mlawer E J, Taubman S J, Brown P Det al, 1997.Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave.J Geophys Res,102(D14): 16663—16682

    Parrish D F, Derber J C, 1992.The national meteorological center’s spectral statistical-interpolation analysis system.Mon Wea Rev, 120(8): 1747—1763

    NUMERICAL ANALYSES ON DEVELOPMENT MECHANISM OF A TYPICAL TROPICAL CYCLONE OVER SOUTHWESTERN INDIAN OCEAN

    WANG Bin1,2, GUO Jing-Tian1,2, CAO Yong-Zheng1, ZHANG Wei1, LIU Gui-Yan2
    (1.North China Sea Marine Forecasting Center of State Oceanic Administration,Qingdao266100,China; 2.Shandong Provincial Key Laboratory of Marine Ecological Environment and Disaster Prevention and Mitigation,Qingdao266100,China)

    The structure, evolution, and physical mechanism of a strong tropical cyclone that occurred over the Southwestern Indian Ocean from 19 to 28 January 2012 were investigated by using all types of observational data and WRF-3.4.1 (Weather Research Forecasting Model) model.The model reproduced the tropical cyclone reasonably including the moving path, strength and inner structure.The analysis based upon the simulation result indicated that the sensible heat flux and latent heat of condensation run through the whole process of tropical cyclone development.The simulating result shows that the latent heat of condensation is the main source of energy for development of the cyclone, and the role of sensible heat flux is relatively weak.The CISK (Conditional Instability of Second Kind) mechanism can explain the development of this tropical cyclone well.

    Southwestern Indian Ocean; tropical cyclone; WRF; development mechanism

    P435

    10.11693/hyhz20161100259

    * 國際海域資源調(diào)查與開發(fā)“十二五”項目, DY125-22-QY-26號; 海洋公益性行業(yè)科研專項項目, 201205010號; 國家自然科學基金項目, 41306028號。王彬, 助理工程師, E-mail: sgwangbin@126.com

    ① 通訊作者: 郭敬天, 教授級高工, E-mail: guojingtian@bhfj.gov.cn

    2016-11-25,收修改稿日期:2017-01-18

    猜你喜歡
    潛熱氣旋印度洋
    2003年3月北大西洋上兩個爆發(fā)性氣旋的“吞并”過程及發(fā)展機制分析?
    與南亞高壓相聯(lián)的歐亞大陸-印度洋經(jīng)向環(huán)流
    氨還原劑噴槍中保護氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    Effect of moxibustion combined with acupoint application on enteral nutrition tolerance in patients with severe acute pancreatitis
    工業(yè)革命時期蒸汽動力的應(yīng)用與熱力學理論的關(guān)系
    去印度洋
    滇池(2018年1期)2018-01-17 21:34:30
    北太平洋上一個爆發(fā)性氣旋族的結(jié)構(gòu)分析?
    2014年3月大西洋上一個爆發(fā)性氣旋的研究
    青藏高原東部夏季降水凝結(jié)潛熱變化特征分析
    堿回收爐空氣加熱器冷凝水系統(tǒng)
    中國造紙(2015年1期)2015-03-05 05:07:12
    人人澡人人妻人| 亚洲,欧美精品.| 人人妻人人添人人爽欧美一区卜| 成人黄色视频免费在线看| 90打野战视频偷拍视频| 国产一区有黄有色的免费视频| 欧美中文综合在线视频| 亚洲精品美女久久久久99蜜臀 | 久久这里只有精品19| 久久99一区二区三区| 久久久欧美国产精品| 老熟女久久久| 国产亚洲av高清不卡| netflix在线观看网站| 亚洲欧美一区二区三区久久| 免费不卡黄色视频| 亚洲一码二码三码区别大吗| 久久亚洲国产成人精品v| 天天躁狠狠躁夜夜躁狠狠躁| 女的被弄到高潮叫床怎么办| 亚洲四区av| 黑人欧美特级aaaaaa片| 97人妻天天添夜夜摸| 亚洲精品av麻豆狂野| 久久久久视频综合| 国产精品无大码| 久久天躁狠狠躁夜夜2o2o | 亚洲,欧美精品.| 在线观看一区二区三区激情| 老司机在亚洲福利影院| 少妇猛男粗大的猛烈进出视频| √禁漫天堂资源中文www| 日韩精品有码人妻一区| 婷婷色麻豆天堂久久| 天堂中文最新版在线下载| 久久久久精品人妻al黑| 欧美精品高潮呻吟av久久| 亚洲精品国产区一区二| 午夜91福利影院| 青春草亚洲视频在线观看| 操美女的视频在线观看| 满18在线观看网站| 亚洲第一青青草原| 国产日韩欧美亚洲二区| 一二三四在线观看免费中文在| 日韩成人av中文字幕在线观看| 亚洲免费av在线视频| 日韩精品有码人妻一区| 欧美日韩成人在线一区二区| 满18在线观看网站| 国产免费现黄频在线看| 国产av一区二区精品久久| 精品亚洲乱码少妇综合久久| 天堂俺去俺来也www色官网| 精品久久久精品久久久| 9色porny在线观看| 国产精品人妻久久久影院| 免费高清在线观看视频在线观看| 另类亚洲欧美激情| 国产欧美日韩一区二区三区在线| 亚洲自偷自拍图片 自拍| 日本猛色少妇xxxxx猛交久久| 久久精品亚洲av国产电影网| 90打野战视频偷拍视频| 国产野战对白在线观看| 97精品久久久久久久久久精品| 欧美久久黑人一区二区| 精品国产乱码久久久久久小说| 日韩,欧美,国产一区二区三区| 免费黄频网站在线观看国产| 国产老妇伦熟女老妇高清| 热re99久久国产66热| 国产片内射在线| 国产97色在线日韩免费| 电影成人av| 丝袜美腿诱惑在线| 国产一区二区激情短视频 | 国产淫语在线视频| 欧美日韩av久久| av在线播放精品| 欧美少妇被猛烈插入视频| 少妇被粗大猛烈的视频| 色精品久久人妻99蜜桃| 欧美日韩亚洲综合一区二区三区_| 亚洲,欧美精品.| 极品少妇高潮喷水抽搐| 黄色视频在线播放观看不卡| 黄片播放在线免费| 亚洲免费av在线视频| 日韩一区二区视频免费看| 亚洲欧洲精品一区二区精品久久久 | 国产精品久久久久久精品电影小说| 国产亚洲最大av| av福利片在线| av.在线天堂| 两个人免费观看高清视频| 久久久久人妻精品一区果冻| 黄片无遮挡物在线观看| 久久精品亚洲熟妇少妇任你| 亚洲精品美女久久久久99蜜臀 | 国产精品一区二区精品视频观看| 免费黄频网站在线观看国产| 亚洲国产精品成人久久小说| 黄色毛片三级朝国网站| 欧美少妇被猛烈插入视频| 亚洲成国产人片在线观看| 伊人亚洲综合成人网| 国产一区二区三区综合在线观看| 高清在线视频一区二区三区| 国产精品一二三区在线看| 热99国产精品久久久久久7| 亚洲国产精品一区三区| 亚洲av福利一区| 亚洲一卡2卡3卡4卡5卡精品中文| 18在线观看网站| 国产成人精品无人区| 夫妻性生交免费视频一级片| 久久婷婷青草| 亚洲少妇的诱惑av| 韩国高清视频一区二区三区| 在线 av 中文字幕| 欧美人与性动交α欧美软件| 久久99热这里只频精品6学生| 久久 成人 亚洲| 青春草亚洲视频在线观看| av天堂久久9| av天堂久久9| 综合色丁香网| 啦啦啦在线观看免费高清www| av福利片在线| 日韩精品有码人妻一区| 亚洲av成人精品一二三区| 免费观看a级毛片全部| 国产黄色免费在线视频| 一级,二级,三级黄色视频| 男的添女的下面高潮视频| 叶爱在线成人免费视频播放| 国产精品av久久久久免费| 中文字幕人妻丝袜制服| 国产免费福利视频在线观看| 制服丝袜香蕉在线| 天堂中文最新版在线下载| 亚洲一码二码三码区别大吗| 国产成人精品久久久久久| 免费女性裸体啪啪无遮挡网站| 中文字幕另类日韩欧美亚洲嫩草| 国产av精品麻豆| 中文字幕人妻丝袜制服| 亚洲人成77777在线视频| 国产在线免费精品| 国产精品一区二区在线观看99| 中文字幕人妻丝袜制服| 一区二区三区精品91| 国产精品一国产av| 少妇被粗大猛烈的视频| 69精品国产乱码久久久| 99香蕉大伊视频| 美女主播在线视频| 日本av免费视频播放| 一区在线观看完整版| bbb黄色大片| 最新在线观看一区二区三区 | 欧美黄色片欧美黄色片| 亚洲精品美女久久久久99蜜臀 | 80岁老熟妇乱子伦牲交| 亚洲人成网站在线观看播放| 欧美激情 高清一区二区三区| 日韩精品免费视频一区二区三区| 一级爰片在线观看| 国产又爽黄色视频| 午夜免费观看性视频| 亚洲精品久久成人aⅴ小说| 日韩精品免费视频一区二区三区| 亚洲精品久久成人aⅴ小说| 亚洲第一av免费看| 久久久国产欧美日韩av| 波多野结衣av一区二区av| 精品亚洲成国产av| 亚洲国产精品国产精品| 久久久国产欧美日韩av| 在线观看免费视频网站a站| 伊人久久国产一区二区| 国产无遮挡羞羞视频在线观看| 国产日韩欧美视频二区| 日韩中文字幕欧美一区二区 | 免费av中文字幕在线| 中国三级夫妇交换| 国产一区二区 视频在线| 精品久久蜜臀av无| av网站在线播放免费| 一本色道久久久久久精品综合| 操美女的视频在线观看| 成人国语在线视频| 亚洲美女视频黄频| 香蕉国产在线看| 哪个播放器可以免费观看大片| 日韩 亚洲 欧美在线| 好男人视频免费观看在线| 咕卡用的链子| 亚洲欧美清纯卡通| 亚洲五月色婷婷综合| 超碰97精品在线观看| 日韩一区二区三区影片| 一边摸一边做爽爽视频免费| 天堂俺去俺来也www色官网| 999精品在线视频| 午夜影院在线不卡| 久久久精品免费免费高清| 亚洲男人天堂网一区| 久久久久精品人妻al黑| 免费高清在线观看日韩| 中文字幕人妻丝袜制服| 国产极品天堂在线| 大香蕉久久成人网| 日本一区二区免费在线视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲一级一片aⅴ在线观看| 欧美日韩亚洲国产一区二区在线观看 | 成人亚洲欧美一区二区av| 国产成人精品久久二区二区91 | 国产亚洲欧美精品永久| 人人妻人人添人人爽欧美一区卜| 中文字幕高清在线视频| 水蜜桃什么品种好| 亚洲精品国产av蜜桃| 欧美在线黄色| 国产一卡二卡三卡精品 | 午夜日本视频在线| 又大又黄又爽视频免费| 熟女少妇亚洲综合色aaa.| 久久久久人妻精品一区果冻| 色播在线永久视频| 免费黄色在线免费观看| 黄色 视频免费看| 国产日韩欧美亚洲二区| 波多野结衣av一区二区av| 亚洲精华国产精华液的使用体验| 成年人免费黄色播放视频| 90打野战视频偷拍视频| 纵有疾风起免费观看全集完整版| 国产精品无大码| 中文字幕色久视频| 观看美女的网站| 又黄又粗又硬又大视频| 国产精品免费视频内射| 亚洲精品在线美女| 国产在线一区二区三区精| 国产一区有黄有色的免费视频| 18禁国产床啪视频网站| 夫妻性生交免费视频一级片| 交换朋友夫妻互换小说| 少妇人妻 视频| 国产精品三级大全| 亚洲 欧美一区二区三区| 2021少妇久久久久久久久久久| 日韩中文字幕欧美一区二区 | 亚洲一区中文字幕在线| 激情视频va一区二区三区| videos熟女内射| 午夜免费观看性视频| 久久久久精品人妻al黑| 国产高清不卡午夜福利| 十八禁人妻一区二区| 少妇人妻 视频| 日本猛色少妇xxxxx猛交久久| 国精品久久久久久国模美| a 毛片基地| 在线观看www视频免费| 丁香六月天网| 久久精品aⅴ一区二区三区四区| kizo精华| 国产熟女欧美一区二区| 在线观看免费日韩欧美大片| 99久久精品国产亚洲精品| 国产人伦9x9x在线观看| av国产精品久久久久影院| 精品卡一卡二卡四卡免费| 国产精品久久久久久精品电影小说| 三上悠亚av全集在线观看| 赤兔流量卡办理| 精品国产露脸久久av麻豆| 亚洲精品自拍成人| 一区二区av电影网| 久久午夜综合久久蜜桃| 王馨瑶露胸无遮挡在线观看| 国产日韩欧美视频二区| 纯流量卡能插随身wifi吗| 亚洲国产看品久久| 日本黄色日本黄色录像| 国产熟女午夜一区二区三区| 日本猛色少妇xxxxx猛交久久| 啦啦啦中文免费视频观看日本| 欧美日韩精品网址| 亚洲少妇的诱惑av| 中文字幕人妻丝袜一区二区 | 久久久久人妻精品一区果冻| 亚洲成人手机| 国产亚洲av高清不卡| 亚洲欧美色中文字幕在线| 亚洲熟女精品中文字幕| 久久久久久久久免费视频了| 中文字幕制服av| 人体艺术视频欧美日本| 熟女少妇亚洲综合色aaa.| 国产精品国产三级国产专区5o| 熟女av电影| 美女主播在线视频| 毛片一级片免费看久久久久| 丁香六月天网| 永久免费av网站大全| 人人妻人人澡人人爽人人夜夜| 51午夜福利影视在线观看| 日韩av免费高清视频| 久久久精品区二区三区| 深夜精品福利| 精品国产一区二区久久| 女人爽到高潮嗷嗷叫在线视频| 日韩一区二区视频免费看| 欧美日韩一级在线毛片| 国产老妇伦熟女老妇高清| 精品国产露脸久久av麻豆| 91精品伊人久久大香线蕉| 国产av一区二区精品久久| 在线观看三级黄色| 欧美97在线视频| 精品一区二区免费观看| 高清在线视频一区二区三区| 日日撸夜夜添| 午夜福利视频在线观看免费| 亚洲av男天堂| 高清在线视频一区二区三区| av片东京热男人的天堂| 精品亚洲乱码少妇综合久久| 视频区图区小说| 悠悠久久av| 中文欧美无线码| 久久人人爽人人片av| 满18在线观看网站| 久久久久网色| 国产精品久久久久久人妻精品电影 | 在线观看一区二区三区激情| av片东京热男人的天堂| 亚洲美女搞黄在线观看| 国产麻豆69| 飞空精品影院首页| 日韩伦理黄色片| 亚洲成人一二三区av| 只有这里有精品99| 一边摸一边做爽爽视频免费| 午夜日本视频在线| 一区二区三区激情视频| 免费人妻精品一区二区三区视频| 激情视频va一区二区三区| 国产极品天堂在线| 日韩 亚洲 欧美在线| 中文天堂在线官网| 欧美日韩av久久| 交换朋友夫妻互换小说| 2021少妇久久久久久久久久久| 色网站视频免费| 少妇猛男粗大的猛烈进出视频| 成年人免费黄色播放视频| 久久久久久久精品精品| 亚洲成人国产一区在线观看 | 在线免费观看不下载黄p国产| 51午夜福利影视在线观看| 女人爽到高潮嗷嗷叫在线视频| 色94色欧美一区二区| 黄色怎么调成土黄色| 蜜桃在线观看..| 夫妻午夜视频| 人妻 亚洲 视频| 国产精品欧美亚洲77777| 巨乳人妻的诱惑在线观看| 久久久久精品久久久久真实原创| 中文字幕亚洲精品专区| 香蕉国产在线看| 一二三四在线观看免费中文在| 亚洲一码二码三码区别大吗| 9热在线视频观看99| 久久人妻熟女aⅴ| 日本欧美国产在线视频| 国产成人精品久久二区二区91 | 久久婷婷青草| 国产人伦9x9x在线观看| 中文字幕人妻丝袜一区二区 | av.在线天堂| 2021少妇久久久久久久久久久| 国产伦人伦偷精品视频| 水蜜桃什么品种好| 只有这里有精品99| 飞空精品影院首页| 日韩不卡一区二区三区视频在线| 老熟女久久久| 十分钟在线观看高清视频www| 波多野结衣av一区二区av| 精品久久蜜臀av无| 欧美变态另类bdsm刘玥| 91老司机精品| 久久久久久久久久久免费av| 99九九在线精品视频| 夫妻性生交免费视频一级片| 亚洲四区av| 纵有疾风起免费观看全集完整版| 亚洲精品第二区| 国产黄频视频在线观看| 性色av一级| 日本午夜av视频| 成年av动漫网址| 99久久精品国产亚洲精品| 日本91视频免费播放| 精品久久蜜臀av无| 一区二区日韩欧美中文字幕| 国产片内射在线| 国产午夜精品一二区理论片| 黄色怎么调成土黄色| 国产一区二区激情短视频 | 久久久久久久久免费视频了| 捣出白浆h1v1| 欧美日韩国产mv在线观看视频| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| 国产精品国产三级专区第一集| 建设人人有责人人尽责人人享有的| 国产精品一国产av| 狠狠精品人妻久久久久久综合| 亚洲国产av新网站| 天天躁狠狠躁夜夜躁狠狠躁| svipshipincom国产片| 国产成人午夜福利电影在线观看| 国产成人av激情在线播放| 不卡视频在线观看欧美| 啦啦啦视频在线资源免费观看| 不卡视频在线观看欧美| 大片免费播放器 马上看| 欧美日韩综合久久久久久| 热99国产精品久久久久久7| 99国产综合亚洲精品| 免费在线观看视频国产中文字幕亚洲 | 亚洲图色成人| 各种免费的搞黄视频| 日韩一区二区视频免费看| 在线观看国产h片| 亚洲精华国产精华液的使用体验| 一区二区三区四区激情视频| 女性生殖器流出的白浆| 国产在视频线精品| 母亲3免费完整高清在线观看| 亚洲熟女毛片儿| av免费观看日本| 国产爽快片一区二区三区| 最新的欧美精品一区二区| 最近的中文字幕免费完整| 久久性视频一级片| 自线自在国产av| 中文字幕色久视频| 波多野结衣av一区二区av| 欧美国产精品一级二级三级| 国产精品二区激情视频| 色播在线永久视频| 欧美精品人与动牲交sv欧美| 亚洲精品在线美女| 日日撸夜夜添| 午夜免费观看性视频| 精品第一国产精品| 天天添夜夜摸| 9热在线视频观看99| 中国国产av一级| 欧美激情极品国产一区二区三区| 人妻人人澡人人爽人人| 亚洲精品视频女| 久久久久久久大尺度免费视频| 国产一区二区在线观看av| 亚洲欧美色中文字幕在线| 亚洲一级一片aⅴ在线观看| 亚洲精品久久成人aⅴ小说| 欧美日韩视频高清一区二区三区二| av有码第一页| 亚洲一区中文字幕在线| 亚洲熟女精品中文字幕| 免费黄网站久久成人精品| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美一区二区三区久久| 亚洲精品美女久久av网站| 欧美久久黑人一区二区| 一本一本久久a久久精品综合妖精| 久久久久久久久免费视频了| 欧美日本中文国产一区发布| 看免费av毛片| 欧美老熟妇乱子伦牲交| www.精华液| 国产又爽黄色视频| av一本久久久久| 宅男免费午夜| 青青草视频在线视频观看| 国产成人一区二区在线| 99热网站在线观看| 午夜影院在线不卡| 午夜激情久久久久久久| 国产免费又黄又爽又色| 日韩大片免费观看网站| 十八禁高潮呻吟视频| 国产成人免费观看mmmm| 亚洲少妇的诱惑av| av片东京热男人的天堂| 午夜福利视频在线观看免费| 国产精品亚洲av一区麻豆 | av片东京热男人的天堂| 97在线人人人人妻| 日韩一卡2卡3卡4卡2021年| 啦啦啦中文免费视频观看日本| 丝袜在线中文字幕| 青草久久国产| 丝袜脚勾引网站| 日韩一本色道免费dvd| 国产精品久久久久久精品电影小说| 午夜福利,免费看| 搡老乐熟女国产| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 熟女av电影| 久久综合国产亚洲精品| 国产探花极品一区二区| 色视频在线一区二区三区| 91精品三级在线观看| xxxhd国产人妻xxx| 91精品三级在线观看| 青青草视频在线视频观看| 深夜精品福利| 国产精品女同一区二区软件| 欧美老熟妇乱子伦牲交| 久久这里只有精品19| 亚洲av综合色区一区| 啦啦啦 在线观看视频| 黄片播放在线免费| av线在线观看网站| 女人高潮潮喷娇喘18禁视频| 成人免费观看视频高清| 青草久久国产| 麻豆精品久久久久久蜜桃| 亚洲精品久久久久久婷婷小说| 日韩一区二区三区影片| 免费观看av网站的网址| 午夜福利影视在线免费观看| 亚洲第一av免费看| 久久久国产一区二区| 久久国产精品男人的天堂亚洲| 视频在线观看一区二区三区| 久热爱精品视频在线9| 大香蕉久久网| 最近中文字幕高清免费大全6| 午夜影院在线不卡| 亚洲综合色网址| 中文字幕制服av| 免费在线观看黄色视频的| 一本久久精品| 免费观看人在逋| 国产精品女同一区二区软件| 亚洲,欧美精品.| xxx大片免费视频| 欧美av亚洲av综合av国产av | 十八禁人妻一区二区| 丝袜美足系列| 嫩草影视91久久| 久久久亚洲精品成人影院| 啦啦啦视频在线资源免费观看| 精品酒店卫生间| 国产精品三级大全| 国产淫语在线视频| 亚洲欧洲国产日韩| 精品国产超薄肉色丝袜足j| 黄频高清免费视频| 欧美久久黑人一区二区| 欧美日韩亚洲高清精品| 人成视频在线观看免费观看| 欧美国产精品一级二级三级| kizo精华| 国产男人的电影天堂91| 美女中出高潮动态图| 国产国语露脸激情在线看| 亚洲国产欧美日韩在线播放| 男女免费视频国产| 久久久久精品性色| 性高湖久久久久久久久免费观看| 国产免费视频播放在线视频| 久热这里只有精品99| 中国国产av一级| 国产精品久久久av美女十八| 国产成人免费观看mmmm| 亚洲av综合色区一区| 亚洲一卡2卡3卡4卡5卡精品中文| 在现免费观看毛片| 国产福利在线免费观看视频| 午夜福利视频精品| 国产成人精品在线电影| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 一边亲一边摸免费视频| 两个人免费观看高清视频| 欧美亚洲 丝袜 人妻 在线| 国产1区2区3区精品| 日本欧美视频一区| 免费人妻精品一区二区三区视频| 国产精品国产三级国产专区5o| av有码第一页| 热re99久久国产66热| 蜜桃在线观看..| 成人漫画全彩无遮挡| 国产又爽黄色视频| 午夜久久久在线观看| 精品亚洲乱码少妇综合久久| 一级毛片 在线播放| 欧美成人午夜精品| 夫妻性生交免费视频一级片| 女人爽到高潮嗷嗷叫在线视频|