• <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
    我的女老师完整版在线观看| 婷婷亚洲欧美| 日本五十路高清| 久久久久久久久大av| www.色视频.com| 99久久久亚洲精品蜜臀av| 亚洲中文字幕日韩| 成人av一区二区三区在线看| 伊人久久精品亚洲午夜| 亚洲av不卡在线观看| 中出人妻视频一区二区| 十八禁网站免费在线| 国内揄拍国产精品人妻在线| 亚洲欧美日韩无卡精品| 久久国产乱子免费精品| 久久婷婷人人爽人人干人人爱| 十八禁网站免费在线| 听说在线观看完整版免费高清| 丝袜美腿在线中文| 精品欧美国产一区二区三| 黄色丝袜av网址大全| 亚洲最大成人中文| 亚洲片人在线观看| 成人无遮挡网站| 欧洲精品卡2卡3卡4卡5卡区| 99久久九九国产精品国产免费| 久久国产乱子免费精品| 天天躁日日操中文字幕| 国产在视频线在精品| 国产白丝娇喘喷水9色精品| 亚洲乱码一区二区免费版| 欧美性猛交黑人性爽| 国产淫片久久久久久久久 | 99精品在免费线老司机午夜| 国产精品综合久久久久久久免费| 综合色av麻豆| 日韩有码中文字幕| 精品国内亚洲2022精品成人| 国产高清视频在线播放一区| 嫩草影院精品99| x7x7x7水蜜桃| 欧美激情国产日韩精品一区| 18禁黄网站禁片免费观看直播| 黄色视频,在线免费观看| 一级黄片播放器| av专区在线播放| 国语自产精品视频在线第100页| 丁香欧美五月| 国产精品久久视频播放| 亚洲在线自拍视频| 一区二区三区高清视频在线| 国产精品一及| 午夜免费激情av| 级片在线观看| 一区二区三区高清视频在线| a级毛片免费高清观看在线播放| av女优亚洲男人天堂| 国产探花在线观看一区二区| 国产乱人视频| 国产一区二区亚洲精品在线观看| 能在线免费观看的黄片| 18+在线观看网站| 亚洲精品在线观看二区| 免费观看人在逋| aaaaa片日本免费| 在线播放无遮挡| av天堂在线播放| 我要看日韩黄色一级片| 免费在线观看成人毛片| 欧美绝顶高潮抽搐喷水| 黄色女人牲交| 久久久色成人| 99久久99久久久精品蜜桃| 中文字幕av成人在线电影| 少妇的逼好多水| av女优亚洲男人天堂| 99热精品在线国产| 国内少妇人妻偷人精品xxx网站| 麻豆国产97在线/欧美| 真人一进一出gif抽搐免费| 淫妇啪啪啪对白视频| 精品一区二区三区人妻视频| 欧美性猛交黑人性爽| 日韩有码中文字幕| 韩国av一区二区三区四区| 一区二区三区高清视频在线| 国产精品免费一区二区三区在线| 久久国产乱子免费精品| 欧美3d第一页| 欧美色欧美亚洲另类二区| 极品教师在线视频| 国产av一区在线观看免费| 亚洲在线观看片| 国内精品久久久久精免费| 国产精品一区二区三区四区免费观看 | 国产成+人综合+亚洲专区| 日本免费一区二区三区高清不卡| 免费大片18禁| 国产精品伦人一区二区| 国产精品1区2区在线观看.| 亚洲精品成人久久久久久| 亚洲,欧美精品.| 久久久色成人| 欧美日韩国产亚洲二区| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 又黄又爽又刺激的免费视频.| 欧美三级亚洲精品| a级毛片a级免费在线| 亚洲成a人片在线一区二区| 日日摸夜夜添夜夜添av毛片 | 成人特级av手机在线观看| 狂野欧美白嫩少妇大欣赏| 精品久久久久久久人妻蜜臀av| 91午夜精品亚洲一区二区三区 | 1024手机看黄色片| 久久久久免费精品人妻一区二区| 3wmmmm亚洲av在线观看| 久久午夜福利片| 天堂影院成人在线观看| 国产爱豆传媒在线观看| 俄罗斯特黄特色一大片| 老司机福利观看| 欧美不卡视频在线免费观看| 精品久久久久久久久av| 成年女人毛片免费观看观看9| 免费看光身美女| 久久99热这里只有精品18| 午夜老司机福利剧场| 国产成人a区在线观看| 国产精品综合久久久久久久免费| 亚洲性夜色夜夜综合| 国产淫片久久久久久久久 | 国产91精品成人一区二区三区| 青草久久国产| 欧美高清成人免费视频www| 国产av在哪里看| 搡老妇女老女人老熟妇| 国产综合懂色| 久久热精品热| 免费人成视频x8x8入口观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品免费久久久久久久清纯| 观看免费一级毛片| 欧美又色又爽又黄视频| 国产精品久久视频播放| 亚洲国产日韩欧美精品在线观看| 搡女人真爽免费视频火全软件 | 我的女老师完整版在线观看| 1024手机看黄色片| 国产精品99久久久久久久久| 色吧在线观看| 99热这里只有是精品在线观看 | 美女cb高潮喷水在线观看| 99热精品在线国产| 国产高清激情床上av| 久久久久九九精品影院| 欧美日韩国产亚洲二区| 在线免费观看的www视频| 1024手机看黄色片| 伦理电影大哥的女人| 黄片小视频在线播放| 国产人妻一区二区三区在| 亚洲,欧美,日韩| 男人舔女人下体高潮全视频| 亚洲人成电影免费在线| 亚洲国产精品合色在线| 99在线人妻在线中文字幕| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 国产欧美日韩精品亚洲av| 在线看三级毛片| 国产毛片a区久久久久| 高清日韩中文字幕在线| 国产成人aa在线观看| 内射极品少妇av片p| 亚洲美女视频黄频| 1024手机看黄色片| 欧美绝顶高潮抽搐喷水| 亚洲综合色惰| 性欧美人与动物交配| 网址你懂的国产日韩在线| 波多野结衣高清作品| 首页视频小说图片口味搜索| 国产午夜福利久久久久久| 淫秽高清视频在线观看| 亚洲久久久久久中文字幕| 亚洲精品在线美女| 午夜亚洲福利在线播放| 十八禁国产超污无遮挡网站| 久久国产乱子伦精品免费另类| 色精品久久人妻99蜜桃| 好看av亚洲va欧美ⅴa在| 一区福利在线观看| 一进一出抽搐动态| 免费无遮挡裸体视频| 十八禁网站免费在线| 国产精品伦人一区二区| 三级国产精品欧美在线观看| 最近最新免费中文字幕在线| 亚洲一区二区三区不卡视频| 波多野结衣高清作品| 麻豆av噜噜一区二区三区| 小说图片视频综合网站| 午夜亚洲福利在线播放| 精品国内亚洲2022精品成人| 国产伦精品一区二区三区视频9| 午夜精品久久久久久毛片777| 欧美乱色亚洲激情| АⅤ资源中文在线天堂| 在线播放国产精品三级| 1000部很黄的大片| 午夜福利在线观看吧| 国产亚洲精品综合一区在线观看| 国产精品一及| 色av中文字幕| 神马国产精品三级电影在线观看| 国产亚洲av嫩草精品影院| 搞女人的毛片| 女生性感内裤真人,穿戴方法视频| 国产精品久久久久久久久免 | 国产精品影院久久| АⅤ资源中文在线天堂| 欧美日韩瑟瑟在线播放| 97碰自拍视频| 国产成年人精品一区二区| 日韩av在线大香蕉| 久久草成人影院| 久久久久久久久久成人| 久久国产乱子免费精品| 免费观看人在逋| 在线看三级毛片| 亚洲国产日韩欧美精品在线观看| 最好的美女福利视频网| 国产三级黄色录像| 无人区码免费观看不卡| 亚洲中文日韩欧美视频| 91字幕亚洲| or卡值多少钱| 蜜桃亚洲精品一区二区三区| 在线观看舔阴道视频| 又粗又爽又猛毛片免费看| 色综合欧美亚洲国产小说| 国产在线精品亚洲第一网站| 亚洲色图av天堂| 精品国产三级普通话版| 91久久精品国产一区二区成人| 直男gayav资源| 国产午夜福利久久久久久| 天堂√8在线中文| 午夜精品久久久久久毛片777| 成人精品一区二区免费| 十八禁网站免费在线| 91麻豆av在线| 日韩人妻高清精品专区| 免费观看的影片在线观看| 欧美绝顶高潮抽搐喷水| 舔av片在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美成人a在线观看| 成人精品一区二区免费| 国产一区二区在线av高清观看| 国产日本99.免费观看| 中文资源天堂在线| 精品一区二区三区人妻视频| 国产一区二区激情短视频| 亚洲一区高清亚洲精品| 成人毛片a级毛片在线播放| 97热精品久久久久久| 欧美日本视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| bbb黄色大片| av天堂中文字幕网| 性色av乱码一区二区三区2| 国内精品久久久久精免费| 国内揄拍国产精品人妻在线| 精品日产1卡2卡| 亚洲精品亚洲一区二区| 国产欧美日韩一区二区精品| 精品一区二区三区视频在线观看免费| 亚洲精品影视一区二区三区av| 亚洲一区二区三区不卡视频| 十八禁网站免费在线| 黄色丝袜av网址大全| 午夜影院日韩av| 亚洲国产色片| 亚洲成人久久性| 免费搜索国产男女视频| 91九色精品人成在线观看| 国产黄a三级三级三级人| 国内久久婷婷六月综合欲色啪| 村上凉子中文字幕在线| 成人av在线播放网站| 变态另类成人亚洲欧美熟女| 51午夜福利影视在线观看| 亚洲中文日韩欧美视频| 国产精品三级大全| 99精品在免费线老司机午夜| 亚洲18禁久久av| 别揉我奶头 嗯啊视频| 亚洲精品在线观看二区| 真人一进一出gif抽搐免费| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 久久久久久大精品| 免费一级毛片在线播放高清视频| АⅤ资源中文在线天堂| 国内精品久久久久精免费| 高清日韩中文字幕在线| 黄色丝袜av网址大全| 18禁黄网站禁片免费观看直播| 国产视频内射| 此物有八面人人有两片| 国产成人aa在线观看| 一区二区三区激情视频| 国内毛片毛片毛片毛片毛片| 国产精品久久电影中文字幕| 欧美乱色亚洲激情| 日韩大尺度精品在线看网址| 午夜影院日韩av| 天天躁日日操中文字幕| 欧美最黄视频在线播放免费| 丁香欧美五月| 午夜福利视频1000在线观看| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 欧美精品国产亚洲| 可以在线观看的亚洲视频| 老司机深夜福利视频在线观看| 少妇人妻一区二区三区视频| 久久精品综合一区二区三区| 久久性视频一级片| 一级av片app| 亚洲三级黄色毛片| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 99精品在免费线老司机午夜| 好男人电影高清在线观看| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 欧美一区二区精品小视频在线| 亚洲av第一区精品v没综合| 丰满人妻熟妇乱又伦精品不卡| 美女大奶头视频| 麻豆av噜噜一区二区三区| av天堂在线播放| 一级黄片播放器| 日韩欧美国产一区二区入口| 色视频www国产| 亚洲精品在线观看二区| 麻豆av噜噜一区二区三区| 国产免费男女视频| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 亚洲av日韩精品久久久久久密| 成人欧美大片| 国产一区二区激情短视频| 窝窝影院91人妻| 变态另类丝袜制服| 看黄色毛片网站| 成年女人永久免费观看视频| 中文字幕av在线有码专区| 成年女人永久免费观看视频| 精品午夜福利视频在线观看一区| 两个人的视频大全免费| 日本黄大片高清| 午夜福利18| 日韩高清综合在线| 窝窝影院91人妻| 亚洲乱码一区二区免费版| 亚洲五月婷婷丁香| 欧美日本视频| 99久久精品热视频| 欧美xxxx性猛交bbbb| 免费av不卡在线播放| 国产精品女同一区二区软件 | 99精品久久久久人妻精品| 午夜激情福利司机影院| 91九色精品人成在线观看| 亚洲国产欧美人成| 亚洲av成人精品一区久久| 搡女人真爽免费视频火全软件 | 国产一区二区三区视频了| 90打野战视频偷拍视频| 级片在线观看| 欧美+亚洲+日韩+国产| 青草久久国产| 两性午夜刺激爽爽歪歪视频在线观看| 高清日韩中文字幕在线| 免费黄网站久久成人精品 | 非洲黑人性xxxx精品又粗又长| 999久久久精品免费观看国产| 久久亚洲真实| 日本熟妇午夜| or卡值多少钱| 窝窝影院91人妻| 我的女老师完整版在线观看| 亚洲国产精品久久男人天堂| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 人妻夜夜爽99麻豆av| 97热精品久久久久久| 久久久久久久久久成人| 哪里可以看免费的av片| 在线播放国产精品三级| 亚洲国产精品sss在线观看| 一进一出抽搐动态| 成年女人看的毛片在线观看| 丝袜美腿在线中文| 欧美日韩国产亚洲二区| 欧美一区二区国产精品久久精品| 天堂网av新在线| 亚洲av免费在线观看| 欧美国产日韩亚洲一区| 少妇丰满av| 免费看日本二区| av中文乱码字幕在线| 18禁黄网站禁片午夜丰满| 不卡一级毛片| 嫩草影院入口| 在线十欧美十亚洲十日本专区| 国产国拍精品亚洲av在线观看| 男人狂女人下面高潮的视频| 一进一出抽搐gif免费好疼| 久久久久国内视频| 欧美+日韩+精品| 日本精品一区二区三区蜜桃| 色综合欧美亚洲国产小说| av在线观看视频网站免费| 精品一区二区三区视频在线| 国产精品国产高清国产av| 综合色av麻豆| 窝窝影院91人妻| 我的女老师完整版在线观看| 性欧美人与动物交配| 久久久久久九九精品二区国产| 日本一二三区视频观看| 成年版毛片免费区| 久久久久久大精品| 十八禁国产超污无遮挡网站| 亚洲国产精品sss在线观看| 亚洲欧美激情综合另类| 一进一出抽搐gif免费好疼| 欧美色视频一区免费| 午夜两性在线视频| 狠狠狠狠99中文字幕| 亚洲,欧美,日韩| 国产精品野战在线观看| 激情在线观看视频在线高清| 国产成人影院久久av| 看十八女毛片水多多多| 禁无遮挡网站| 男人舔奶头视频| 麻豆av噜噜一区二区三区| 亚洲成人中文字幕在线播放| 国产综合懂色| 有码 亚洲区| 一本综合久久免费| 老鸭窝网址在线观看| 小说图片视频综合网站| 日韩有码中文字幕| 日韩成人在线观看一区二区三区| 亚洲第一电影网av| 少妇的逼好多水| 国产伦一二天堂av在线观看| 欧美成人性av电影在线观看| 亚洲第一区二区三区不卡| 国产成人aa在线观看| 亚洲国产精品999在线| 日韩欧美在线二视频| 国产高清视频在线观看网站| 免费av观看视频| 狠狠狠狠99中文字幕| 久久国产精品影院| 简卡轻食公司| 美女大奶头视频| 国产成人啪精品午夜网站| 在线免费观看的www视频| 国产伦精品一区二区三区四那| 观看美女的网站| 波野结衣二区三区在线| 午夜精品在线福利| av在线观看视频网站免费| 国产成人啪精品午夜网站| 亚洲精品色激情综合| 中文字幕高清在线视频| 一夜夜www| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 午夜福利欧美成人| 国产成人a区在线观看| 亚洲性夜色夜夜综合| 成年版毛片免费区| 成人性生交大片免费视频hd| 少妇人妻一区二区三区视频| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 网址你懂的国产日韩在线| 精品欧美国产一区二区三| 国产午夜精品久久久久久一区二区三区 | 波野结衣二区三区在线| 99精品久久久久人妻精品| 亚洲精品日韩av片在线观看| av黄色大香蕉| 亚洲最大成人av| 激情在线观看视频在线高清| 成人av在线播放网站| 无遮挡黄片免费观看| 午夜久久久久精精品| 床上黄色一级片| 国产探花极品一区二区| 禁无遮挡网站| 亚洲欧美日韩无卡精品| 丰满人妻熟妇乱又伦精品不卡| 久久久久久久午夜电影| 99久久九九国产精品国产免费| 99在线视频只有这里精品首页| av中文乱码字幕在线| 又爽又黄a免费视频| 国产黄a三级三级三级人| 赤兔流量卡办理| 丰满乱子伦码专区| 中文字幕高清在线视频| 别揉我奶头 嗯啊视频| 午夜福利高清视频| 天天躁日日操中文字幕| 国产精品亚洲美女久久久| 男女床上黄色一级片免费看| 亚洲第一电影网av| 脱女人内裤的视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲乱码一区二区免费版| 麻豆成人午夜福利视频| 小说图片视频综合网站| 日本a在线网址| 成人永久免费在线观看视频| 51午夜福利影视在线观看| 内射极品少妇av片p| 亚洲av第一区精品v没综合| 网址你懂的国产日韩在线| 欧美成人a在线观看| 变态另类成人亚洲欧美熟女| 嫁个100分男人电影在线观看| 一区二区三区激情视频| 性色avwww在线观看| 亚洲电影在线观看av| 精品人妻熟女av久视频| www日本黄色视频网| 亚洲在线自拍视频| а√天堂www在线а√下载| 不卡一级毛片| 亚洲av电影不卡..在线观看| 九色成人免费人妻av| 看十八女毛片水多多多| 国产精品美女特级片免费视频播放器| 国内揄拍国产精品人妻在线| 日日摸夜夜添夜夜添小说| 91久久精品国产一区二区成人| 亚洲 欧美 日韩 在线 免费| 亚洲欧美精品综合久久99| АⅤ资源中文在线天堂| 看十八女毛片水多多多| bbb黄色大片| 香蕉av资源在线| 国产精品人妻久久久久久| 桃红色精品国产亚洲av| 亚洲真实伦在线观看| 99热这里只有是精品50| 精品免费久久久久久久清纯| 午夜福利高清视频| 日韩欧美在线二视频| 最近视频中文字幕2019在线8| 国产视频一区二区在线看| 亚洲一区高清亚洲精品| 亚洲精品色激情综合| 国产毛片a区久久久久| 熟女人妻精品中文字幕| 香蕉av资源在线| 欧美日韩国产亚洲二区| 99在线视频只有这里精品首页| 中文字幕人成人乱码亚洲影| 中文字幕av在线有码专区| 亚洲精品久久国产高清桃花| 亚洲自拍偷在线| 国产欧美日韩精品一区二区| 欧美黑人巨大hd| 婷婷精品国产亚洲av在线| 最近最新中文字幕大全电影3| 嫩草影院新地址| 12—13女人毛片做爰片一| www.熟女人妻精品国产| 亚洲最大成人av| 欧美日韩黄片免| 国产v大片淫在线免费观看| 琪琪午夜伦伦电影理论片6080| 变态另类成人亚洲欧美熟女| 欧美中文日本在线观看视频| 日本精品一区二区三区蜜桃| 精品福利观看| 97超级碰碰碰精品色视频在线观看| 亚洲自偷自拍三级| 婷婷亚洲欧美| 色视频www国产| 色哟哟·www| 麻豆成人av在线观看| 首页视频小说图片口味搜索| 久久精品影院6| 变态另类成人亚洲欧美熟女| 欧美丝袜亚洲另类 | 亚洲内射少妇av| 国产精品98久久久久久宅男小说| 欧美黄色淫秽网站|