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

    熱帶大氣對單一型赤道非對稱熱源的響應(yīng)

    2014-08-04 08:51:56邢楠1李建平1李耀錕3
    大氣科學(xué) 2014年6期
    關(guān)鍵詞:海溫氣旋赤道

    邢楠1, 2, 3 李建平1, 3 李耀錕3

    ?

    熱帶大氣對單一型赤道非對稱熱源的響應(yīng)

    邢楠李建平李耀錕

    1中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室,北京100029;2中國科學(xué)院大學(xué),北京100049;3北京師范大學(xué)全球變化與地球系統(tǒng)科學(xué)研究院,北京100875

    本文采用Gill模式得到了熱帶大氣對單一型赤道非對稱熱源響應(yīng)的理論解析通解,從理論上完善了單一型赤道非對稱熱源激發(fā)的赤道非對稱的大氣響應(yīng)結(jié)果。同時(shí)在單一型赤道非對稱熱源的位置、強(qiáng)度及范圍變化對大氣響應(yīng)的影響方面做了詳細(xì)的研究。當(dāng)熱源中心位置北移,北半球氣旋強(qiáng)度增加、位置北移,同時(shí)赤道輻合氣流減弱而越赤道氣流增強(qiáng);當(dāng)熱源強(qiáng)度增強(qiáng)(減弱),熱源激發(fā)的大氣響應(yīng)整體增強(qiáng)(減弱),但大氣分布型不發(fā)生變化;當(dāng)熱源范圍不斷增大(減小)時(shí),北半球氣旋強(qiáng)度增強(qiáng)(減弱)、位置西移(東移)、范圍增大(減?。?,同時(shí)越赤道氣流增強(qiáng)(減弱)。將上述結(jié)論應(yīng)用于分析孟加拉灣地區(qū)海溫對夏季風(fēng)爆發(fā)影響的研究,指出當(dāng)孟加拉灣地區(qū)經(jīng)向最大暖海溫位于赤道附近時(shí),其兩側(cè)表現(xiàn)為Rossby波響應(yīng)的Gill型氣旋環(huán)流,而海溫暖軸北移后,其南側(cè)激發(fā)出有利于季風(fēng)爆發(fā)的越赤道氣流,這是Rossby重力混合波對熱源響應(yīng)的結(jié)果。這是上述理論結(jié)果的一個(gè)很好例證,同時(shí)也為孟加拉灣夏季風(fēng)的爆發(fā)給出一種動(dòng)力學(xué)解釋。

    單一型赤道非對稱熱源 熱源強(qiáng)迫 孟加拉灣夏季風(fēng)爆發(fā) 海溫

    1 引言

    熱帶大氣除自身的動(dòng)力學(xué)約束外,也受非絕熱加熱的直接驅(qū)動(dòng)(Lorenz,1955;Webster,1972)。作為大氣重要外強(qiáng)迫源的海洋,其海表面溫度(SST)與大尺度環(huán)流系統(tǒng)如Hadley及Walker環(huán)流變化緊密相關(guān)(Bjerknes,1966,1969;Feng et al., 2011,2013;Feng and Li,2013;李建平等,2011)。所以,熱源通過對大氣環(huán)流的強(qiáng)迫影響局地及更大范圍地區(qū)的天氣氣候(Latif et al., 1998;Saji et al., 1999)。

    一些理論研究與數(shù)值模擬工作解釋及驗(yàn)證 了熱帶熱源對熱帶及熱帶外大氣環(huán)流的重要影響。Webster(1972)采用包含緯向基流的線性方 程,得到了赤道對稱熱源的解析解,而且用赤道Kelvin波響應(yīng)解釋了赤道對稱熱源東側(cè)激發(fā)出的東風(fēng)異常。Gill(1980)采用簡單的淺水波方程,在Webster(1972)理論解的基礎(chǔ)上,用東傳的Kelvin波及西傳的Rossby波解釋赤道對稱熱源兩側(cè)的緯向非對稱的東、西風(fēng)異常及赤道外南、北半球?qū)ΨQ的氣旋環(huán)流異常的大氣響應(yīng),該赤道對稱熱源的大氣響應(yīng)被稱為Gill響應(yīng)。Gill響應(yīng)作為基礎(chǔ)理論被學(xué)者用于解釋實(shí)際大氣中的環(huán)流特征(Guan et al., 2003;劉琳和于衛(wèi)東,2006;薛洪斌等,2006)。同時(shí),他將理想的赤道對稱熱源與反對稱熱源解析解簡單疊加得到赤道非對稱熱源的大氣響應(yīng)特征。對于較為復(fù)雜的熱源形式,很多學(xué)者也采用了數(shù)值模式模擬大氣對熱源強(qiáng)迫的響應(yīng)(Webster, 1972;Hoskins and Karoly,1981;Sashegyi and Geisler,1987;Jin and Hoskins,1995;Bretherton and Sobel,2003),得到的共同結(jié)論是赤道對稱熱源激發(fā)出赤道對稱的大氣響應(yīng),而赤道非對稱熱源激發(fā)出非對稱的大氣響應(yīng)。綜上研究,Gill模式理論解(Gill,1980)很好地解釋了熱帶大尺度環(huán)流系統(tǒng)對海溫強(qiáng)迫的響應(yīng),對于一般熱源具有普適性。但Gill(1980)文中關(guān)于赤道非對稱熱源由赤道對稱及反對稱兩個(gè)熱源疊加而成,這兩項(xiàng)恰好是Weber函數(shù)展開的前兩項(xiàng)。那么,如果采用一般形式的非對稱熱源,其解析形式如何,目前還尚未有人研究。由于模式、外強(qiáng)迫熱源不統(tǒng)一,非對稱熱源的大氣響應(yīng)中的越赤道氣流大小、氣旋位置、走向及強(qiáng)度不盡相同。并且赤道熱源的加熱尺度及緯度位置對激發(fā)的大氣環(huán)流結(jié)構(gòu)有一定的影響(Webster,1981;巢紀(jì)平和王彰貴,1991),那么系統(tǒng)分析不同經(jīng)向位置、強(qiáng)度及范圍的赤道非對稱熱源激發(fā)的大氣響應(yīng)對于深入理解熱源強(qiáng)迫的大氣環(huán)流特征尤為重要,但目前這方面的研究還比較少。

    實(shí)際的大氣和海洋中大都可以找到類似的單一型赤道非對稱熱源,如局地強(qiáng)對流活動(dòng)的凝結(jié)潛熱釋放(Ren and Huang,2002)和局地海表升溫等效應(yīng)可近似視為自高值中心向四周漸進(jìn)減弱的單一型分布。赤道非對稱的熱源對低空越赤道氣流存在一定影響,局地?zé)嵩磳υ匠嗟罋饬鞯挠绊懜语@著,并且其強(qiáng)弱也受熱源經(jīng)向梯度的影響(唐東昇和謝立安,1986;吳秋英和朱乾根,1990;劉秦玉和王韶霞,2000;黃亮等,2012)。Jiang and Li (2011)研究發(fā)現(xiàn)局地經(jīng)向SST的最大值位置(海溫暖軸)的北移提前于季風(fēng)爆發(fā),這種局地單一型熱源位置的北移對季風(fēng)爆發(fā)的影響也得到了浮標(biāo)觀測數(shù)據(jù)的證實(shí)(Yu et al., 2012)。因此,研究這種單一型赤道非對稱熱源對大氣環(huán)流及天氣可能造成的影響具有實(shí)際意義,也為進(jìn)一步研究復(fù)合型赤道非對稱熱源對熱帶大氣的強(qiáng)迫提供一定的參考。

    為了完善單一型赤道非對稱熱源在Gill模式中的理論解析解,本研究采用能表征外強(qiáng)迫源在靜止大氣中的影響的理論模式(Gill,1980),對以下科學(xué)問題進(jìn)行分析:當(dāng)單一型熱源表現(xiàn)為赤道非對稱時(shí),能否在Gill模式中得到相應(yīng)的解析解?并且熱源的強(qiáng)度、中心位置及范圍變化時(shí),大氣的響應(yīng) 是如何變化的?由于局地加熱與季風(fēng)爆發(fā)密切相關(guān),那么利用帶有外部熱源Gill模式的解析解能否為孟加拉灣夏季風(fēng)建立的重要過程提供可能的動(dòng)力學(xué)解釋?

    2 赤道非對稱熱源、理論模式和資料

    2.1 單一型赤道非對稱熱源

    由于實(shí)際大氣中更多的熱源是自高值中心向四周漸進(jìn)減弱的單一型熱源分布(Jiang and Li,2011;Yu et al., 2012),所以基于實(shí)際大氣,本文的赤道非對稱熱源只有一個(gè)熱源中心。我們采用與Gill(1980)文中相似的熱源函數(shù),通過改變方向熱源的中心位置,得到赤道非對稱的熱源,其公式如下:

    其中,代表熱源強(qiáng)度,為赤道到熱源中心的經(jīng)向距離,>0、<0分別代表熱源中心位于南、北半球,2為熱源的緯向?qū)挾取?/p>

    連續(xù)的函數(shù)()可以展開為Weber函數(shù)(拋物線柱函數(shù))

    其中

    通過公式(2)展開單一型赤道非對稱熱源公式(1)可得

    其中

    , (4)

    (5)

    2.2 Gill模型中對應(yīng)單一型非對稱熱源的解析解

    利用無量綱的正壓淺水方程,考慮方程的定常解并引入耗散,方程簡化如下(Matsuno,1966;Gill,1980):

    , (7)

    , (8)

    (10)

    (12)

    接下來將熱源展開項(xiàng)逐一帶入公式(10)、(11)及(12),求其對應(yīng)的方程解析解。

    當(dāng)=0時(shí),該解代表東傳的阻尼Kelvin波響應(yīng),

    (13)

    各變量的解的形式如下:

    當(dāng)=1時(shí),該解代表西傳的Rossby波解,

    (15)

    各變量的解的形式如下:

    當(dāng)=0時(shí),該解代表羅斯貝重力混合波響應(yīng),

    當(dāng)=2時(shí),該解代表Rossby波響應(yīng),

    (18)

    各變量的解的形式如下:

    , (20)

    所以,對于>1的每個(gè)熱源展開項(xiàng),方程滿足=+1,只有一種解的形式,其通解如下:

    其中

    各變量的解的形式如下:

    所以,單一型非對稱熱源的解是各熱源展開項(xiàng)帶入方程得到的解析解之和。取衰減系數(shù),下面我們將給出解析解的分布圖,以分析大氣對單一型非對稱熱源的響應(yīng)特征。

    2.3 資料

    本文使用的資料包括:(1)1982~2010年NCEP/ NCAR月平均再分析資料(Kalnay et al., 1996);(2)1982~2010年NOAA第二版的高分辨逐月海表面溫度(OISST)資料(Reynolds et al., 2002)。

    3 單一型赤道非對稱熱源的分布及其大氣響應(yīng)特征

    3.1 單一型赤道非對稱冷、熱源的分布

    取冷(熱)源強(qiáng)度A為-1(1)、范圍2為4以及中心位置=1、=-1,由公式(1)得到四組理想單一型赤道非對稱冷、熱源(圖1)。冷、熱源分布在20°S~20°N(1個(gè)單位約為10°),熱源中心分別位于10°S(圖1a)、10°N(圖1b)以及冷源中心分別位于10°S(圖1c)、10°N(圖1d),冷、熱源強(qiáng)度向東西及南北兩側(cè)遞減。

    圖1 由公式(1)得到的單一型赤道非對稱熱、冷源的分布。(a)熱源中心位于d=1;(b)同(a)但為d=-1;(c)同(a)但為冷源;(d)同(b)但為冷源。這里d為赤道到熱源中心的經(jīng)向距離,d>0、d<0分別表示熱源位于南、北半球。實(shí)線表示熱源,虛線表示冷源,等值線間隔為0.2

    圖1所示的單一型赤道非對稱冷、熱源可以展開為多冷、熱源分布型疊加的形式。以北半球單一型赤道非對稱熱源(圖1b)為例,其展開的熱源分布如圖2所示。由圖可見,單一型赤道非對稱熱源可以分解為赤道對稱型(圖2a、c、e)與反對稱型熱源(圖2b、d),經(jīng)向熱源中心數(shù)隨著展開項(xiàng)數(shù)的增加而增加,同時(shí)熱源強(qiáng)度減小。熱源展開項(xiàng)的前五項(xiàng)的合成(圖2f)與公式(1)得到的單一型赤道非對稱熱源的絕對誤差小于5%(圖略),所以以下將單一型赤道非對稱熱源對應(yīng)的解析解近似為前五項(xiàng)展開熱源對應(yīng)的解析解之和。

    圖2 圖1b中的單一型赤道非對稱熱源的前五個(gè)展開熱源的(a–e)分布及其(f)合成,這里n代表拋物柱面函數(shù)展開的階數(shù)

    圖3 圖1中南、北半球的單一型赤道非對稱熱、冷源分別激發(fā)的低層水平風(fēng)(箭頭)、垂直速度(填色)以及氣壓場(粗等值線)。(a)圖1a中的熱源得到的大氣變量解析解;(b)同(a)但為圖1b中的熱源;(c)同(a)但為圖1c中的熱源;(d)同(a)但為圖1d中的熱源。粗實(shí)線、虛線分別代表正、負(fù)氣壓值,等值線間隔是0.3

    3.2 大氣對單一型赤道非對稱冷、熱源強(qiáng)迫的響應(yīng)

    圖3是單一型赤道非對稱熱源通過Gill模式得到的風(fēng)、垂直速度及氣壓解析解的分布。由圖3a、圖3b可見,當(dāng)熱源中心南(北)移,垂直上升中心也出現(xiàn)南(北)移,上升中心與熱源中心保持一致,同時(shí)較弱的下沉氣流位于上升中心西側(cè);赤道附近,熱源東、西兩側(cè)分別是與Rossby波及Kelvin波響應(yīng)相對的東、西風(fēng)向熱源中心輻合,熱源中心西南(北)側(cè)出現(xiàn)Rossby波響應(yīng)的氣旋環(huán)流;與赤道對稱熱源激發(fā)的大氣響應(yīng)相比,赤道非對稱熱源激發(fā)的南(北)半球的低壓強(qiáng)度大于北(南)半球的,所以南北半球出現(xiàn)了經(jīng)向氣壓梯度力,在Rossby重力混合波的作用下,北(南)半球的氣旋環(huán)流消失,取而代之的是越赤道氣流;赤道外風(fēng)場滿足準(zhǔn)地轉(zhuǎn)平衡,氣旋中心與低壓中心一致。單一型赤道非對稱熱源的赤道非對稱大氣響應(yīng)與先前的結(jié)果一致,但本研究中理論解得到的氣旋環(huán)流基本呈東西走向,所以在氣旋環(huán)流形狀、位置等細(xì)節(jié)上的結(jié)果與之前存在一定差別(Webster,1972;Bretherton and Sobel,2003)。當(dāng)冷源中心南(北)移后,下沉氣流也隨之南(北)移,較弱的上升氣流位于下沉中心西側(cè);赤道附近,冷源東、西兩側(cè)分別是與Rossby波及Kelvin波響應(yīng)相對的西、東風(fēng),冷源中心西南(北)側(cè)出現(xiàn)Rossby波響應(yīng)的反氣旋環(huán)流,赤道非對稱冷源所激發(fā)的Rossby重力混合波響應(yīng)也引起向北(南)半球的越赤道氣流;風(fēng)場同樣滿足準(zhǔn)地轉(zhuǎn)平衡,所以反氣旋環(huán)流中心與高壓中心對應(yīng)(圖3c和d)。總之,赤道非對稱的冷、熱源激發(fā)赤道非對稱的大氣響應(yīng),而冷、熱源激發(fā)的大氣響應(yīng)相反。

    實(shí)際大氣中,由于外強(qiáng)迫源位置、強(qiáng)度、范圍各異,同時(shí)由Gill模式理論解可見,赤道對稱及非對稱外強(qiáng)迫源所激發(fā)的大氣分布存在顯著差別(Gill,1980),所以不同分布型的外強(qiáng)迫源的大氣響應(yīng)也是本文研究的重點(diǎn)問題。由圖3可見,中心位置相同的南、北半球熱源激發(fā)的大氣響應(yīng)呈對稱分布,并且中心位置相同的同一半球冷、熱源激發(fā)的大氣響應(yīng)正好相反。所以,在后面的研究中我們以北半球熱源為例,通過一系列的試驗(yàn)設(shè)計(jì)來考察熱源的中心位置、強(qiáng)度、范圍變化引起的大氣響應(yīng)的變化。

    4 大氣對不同熱源分布型的響應(yīng)

    4.1 大氣對赤道強(qiáng)迫熱源的敏感性

    本文給出三組試驗(yàn)設(shè)計(jì),各試驗(yàn)中理想熱源的影響因子如表1所示。試驗(yàn)1中,我們設(shè)定強(qiáng)度及緯向?qū)挾炔蛔兊耐鈴?qiáng)迫熱源,而熱源中心位置分別設(shè)定在=0、=-0.5、=-1及=-1.5,其所在的緯度分別為0°(試驗(yàn)1A)、5°N(試驗(yàn)1B)、10°N(試驗(yàn)1C)以及15°N(試驗(yàn)1D),用來考察熱源中心位置改變引起的大氣響應(yīng)的變化;試驗(yàn)2中,我們設(shè)定緯向?qū)挾燃爸行奈恢貌蛔兊耐鈴?qiáng)迫熱源,而熱源強(qiáng)度分別設(shè)定為(試驗(yàn)2A)、=1(試驗(yàn)2B)和=1.5(試驗(yàn)2C),用來研究熱源強(qiáng)度對大氣響應(yīng)的影響;試驗(yàn)3中,我們設(shè)定強(qiáng)度及中心位置不變的外強(qiáng)迫熱源,而熱源緯向范圍分別設(shè)定為2=2(試驗(yàn)3A),3=4(試驗(yàn)3B)以及2=6(試驗(yàn)3C),以研究熱源范圍對大氣響應(yīng)的影響。

    表1 試驗(yàn)中理想熱源的影響因子的無量綱值

    4.1.1 熱源中心經(jīng)向位置變化

    為了分析赤道非對稱熱源與赤道對稱熱源強(qiáng)迫的垂直速度、水平風(fēng)及氣壓的差異,圖4給出了兩類熱源得到的解析解的差值。圖4a為中心在=-0.5(5°N)與=0(0°)的熱源強(qiáng)迫的大氣響應(yīng)的差值,上升氣流中心位于北半球,其西北側(cè)為氣旋環(huán)流,對應(yīng)一個(gè)低壓中心;較弱的下沉中心位于南半球,其西南側(cè)為反氣旋環(huán)流,對應(yīng)著弱的高壓中心;經(jīng)向的熱源梯度形成了向北的越赤道氣流。中心位置分別取=-1(10°N)及=-15(15°N)的熱源激發(fā)的大氣響應(yīng)與赤道對稱熱源的大氣響應(yīng)的差值結(jié)果表明,隨著熱源中心的北移,北半球氣旋環(huán)流增強(qiáng)且北移,南半球的反氣旋環(huán)流增強(qiáng)且北移,赤道上緯向氣流向東西兩側(cè)輻散,同時(shí)越赤道氣流也隨著經(jīng)向熱源梯度的增強(qiáng)而增強(qiáng)(圖4b和c)。所以隨著熱源北移,熱源中心西北側(cè)的氣旋環(huán)流不斷北移且強(qiáng)度逐漸增強(qiáng),南半球的氣旋環(huán)流逐漸減弱甚至消失,同時(shí)增強(qiáng)的熱源經(jīng)向梯度加強(qiáng)了越赤道氣流,而向熱源輻合的赤道緯向氣流減弱。本文選取最大越赤道氣流、最低氣壓、赤道最大西風(fēng)及最大東風(fēng)為指標(biāo),定量分析這些大氣指標(biāo)隨熱源北移的演變情況(表2)。隨熱源中心北移,最大越赤道氣流以及最低氣壓增強(qiáng),而赤道最大西風(fēng)及東風(fēng)減弱。

    表2 試驗(yàn)1中最大越赤道氣流、最低氣壓、最大赤道西風(fēng)及最大赤道東風(fēng)值

    所以,當(dāng)熱源中心北移,北半球低壓強(qiáng)度大于南半球,熱源的經(jīng)向梯度形成向北的越赤道氣流。并且隨著熱源北移,熱源西北側(cè)的低壓中心北移且增強(qiáng),越赤道氣流也隨著熱源經(jīng)向梯度的增強(qiáng)而不斷增強(qiáng),而赤道上水平氣流減弱。

    圖4 試驗(yàn)(a)1B、(b)1C、(c)1D分別得到的低層垂直速度、水平風(fēng)以及氣壓與試驗(yàn)1A的差值。圖中陰影區(qū)的細(xì)實(shí)線代表上升氣流,細(xì)虛線代表下沉氣流,等值間隔為0.15;粗實(shí)線及虛線分別代表正氣壓及負(fù)氣壓值,等值線間隔是0.3

    4.1.2 熱源強(qiáng)度變化

    圖5為試驗(yàn)2A和2C減去試驗(yàn)2B得到的垂直速度、水平風(fēng)及氣壓的差值分布圖。強(qiáng)度的熱源減去的熱源的大氣響應(yīng)特征為:下沉氣流中心位于北半球,其西北側(cè)分布著一個(gè)與高壓中心對應(yīng)的反氣旋環(huán)流,赤道上緯向氣流向東西兩側(cè)輻散,同時(shí)存在向南的越赤道氣流(圖5a)。而強(qiáng)度=1.5與=1的熱源激發(fā)的大氣響應(yīng)的差值(圖5b)與圖5a的剛好相反,北半球有一個(gè)上升氣流中心,其西北側(cè)分布一個(gè)與低壓中心對應(yīng)的氣旋環(huán)流,赤道上緯向風(fēng)向氣流上升中心輻合,同時(shí)存在向北的越赤道氣流。由于本組試驗(yàn)設(shè)計(jì)的熱源經(jīng)向位置、緯向?qū)挾炔蛔?,所以?qiáng)度較小的熱源減去強(qiáng)度較大的熱源得到的是負(fù)值,大氣響應(yīng)表現(xiàn)為冷源激發(fā)的大氣響應(yīng)特征(圖5a),反之大氣響應(yīng)則為熱源強(qiáng)迫的大氣特征(圖5b)。通過定量分析熱源強(qiáng)度與大氣響應(yīng)變量的關(guān)系,越赤道氣流最大值、最低氣壓、赤道最大西風(fēng)及最大東風(fēng)與熱源強(qiáng)度呈正比(表3)。所以,熱源強(qiáng)度越強(qiáng),大氣響應(yīng)越強(qiáng)。

    總之,當(dāng)赤道非對稱熱源增強(qiáng)或減弱時(shí),激發(fā)的大氣環(huán)流分布型不變,但大氣響應(yīng)的強(qiáng)度隨熱源強(qiáng)度的增強(qiáng)而增強(qiáng),隨熱源強(qiáng)度的減弱而減弱。

    表3 同表2,但為試驗(yàn)2的結(jié)果

    4.1.3 熱源緯向?qū)挾茸兓?/p>

    圖6為試驗(yàn)3A和3C減去試驗(yàn)3B得到的垂直速度、水平風(fēng)及氣壓的差值分布圖。由圖6a可見,北半球分布著兩個(gè)下沉中心,并且下沉中心的西北側(cè)各有一個(gè)與高壓中心相對應(yīng)的反氣旋,赤道上緯向氣流向東西兩側(cè)輻散,同時(shí)存在向南的越赤道氣流。由經(jīng)向范圍與的熱源強(qiáng)迫的大氣響應(yīng)的差值可見,北半球?yàn)閮蓚€(gè)上升中心,并且其西北側(cè)各分布一個(gè)與低壓中心相對應(yīng)的氣旋環(huán)流,赤道上緯向氣流向上升中心輻合,并且存在向北的越赤道氣流(圖6b)。對比圖6的兩個(gè)試驗(yàn)差值分布可見,不同范圍的熱源激發(fā)的大氣響應(yīng)的差值分布都存在兩個(gè)大氣響應(yīng)中心;相對于圖6a,圖6b熱源中心兩側(cè)的氣旋環(huán)流距離較遠(yuǎn),同時(shí)赤道西風(fēng)覆蓋范圍偏東;從大氣響應(yīng)差值的強(qiáng)度來看,圖6a中的氣旋及垂直速度的強(qiáng)度大于圖6b中的反氣旋及垂直速度的強(qiáng)度。這些結(jié)果表明不同緯向?qū)挾鹊臒嵩醇ぐl(fā)的氣旋環(huán)流中心位置及強(qiáng)度不同,熱源范圍越大,熱源激發(fā)的氣旋環(huán)流范圍越廣、位置越偏離熱源中心,同時(shí)大氣響應(yīng)增幅隨緯向?qū)挾鹊脑黾佣兙彙>唧w到大氣變量,隨熱源范圍其變化如下:低壓中心、越赤道氣流、赤道西風(fēng)及東風(fēng)最大值隨熱源范圍的增加而增強(qiáng)、隨熱源寬度的減小而減弱;并且對大氣變量的定量分析也表明,隨著熱源范圍的增加,大氣變量的增幅變緩(表4)。

    圖 6 同圖4,但為(a)試驗(yàn)3A減去3B以及(b)3C減去3B

    表4 同表2,但為試驗(yàn)3的結(jié)果

    所以,當(dāng)熱源范圍減小,北半球氣旋強(qiáng)度減弱、位置東移,向熱源輻合的赤道上的東、西風(fēng)及越赤道氣流減弱;熱源范圍增大時(shí),北半球氣旋強(qiáng)度增強(qiáng)、范圍增大,中心位置偏西,同時(shí)向熱源輻合的越赤道氣流以及赤道上的東、西風(fēng)增強(qiáng)。

    4.2 孟加拉灣地區(qū)海溫對季風(fēng)爆發(fā)的影響

    上面的理論分析表明,不同類型的局地加熱對大氣響應(yīng)起著非常重要的作用。大范圍熱源的季節(jié)性移動(dòng)與夏季風(fēng)的爆發(fā)密切相關(guān),而海溫分布與熱源的季節(jié)移動(dòng)有密切聯(lián)系。所以海溫通常被用來分析和解釋熱帶海氣相互作用。作為亞洲夏季風(fēng)爆發(fā)最早的孟加拉灣地區(qū),已有的研究表明局地?zé)嵩捶植寂c季風(fēng)爆發(fā)之間存在一定的相關(guān)關(guān)系(Liu et al., 2002; Jiang and Li, 2011; Yu et al., 2012)。Jiang and Li(2011)以及Yu et al.(2012)用再分析資料及浮標(biāo)觀測數(shù)據(jù)研究指出海溫暖軸的北跳會激發(fā)其南側(cè)的西南風(fēng),繼而改變局地風(fēng)向、引發(fā)季風(fēng)的爆發(fā)。那么,可否從理論結(jié)果出發(fā)給予孟加拉灣季風(fēng)爆發(fā)一個(gè)很好的動(dòng)力學(xué)解釋呢?

    為了進(jìn)一步分析不同海溫分布下的實(shí)際大氣分布特征,同時(shí)與理想熱源分布下的大氣響應(yīng)特征比對,圖7給出了孟加拉灣地區(qū)海溫與風(fēng)場的分布及與海溫分布相似的理想熱源激發(fā)的大氣響應(yīng)分布。3月份海溫暖軸位于赤道附近,海溫暖軸兩側(cè)表現(xiàn)為氣旋環(huán)流,赤道地區(qū)表現(xiàn)為向印度洋東部輻合的特征(圖7a)。觀測結(jié)果與赤道對稱熱源激發(fā)的大氣響應(yīng)結(jié)果一致(圖7b),關(guān)于赤道對稱的南北半球氣旋環(huán)流是赤道對稱熱源激發(fā)出的Rossby波西傳的結(jié)果(Gill,1980;陳哲和李崇銀,2006)。當(dāng)海溫暖軸北移至孟加拉灣中部,暖軸南側(cè)為西南風(fēng),經(jīng)向海溫梯度形成向北的越赤道氣流,同時(shí)海溫暖軸的北側(cè)表現(xiàn)為氣旋切變特征(圖7b)。孟加拉灣地區(qū)海氣相互作用的大氣環(huán)流特征與理想的赤道非對稱熱源激發(fā)的大氣響應(yīng)一致,即在Rossby重力混合波的作用下,暖中心南側(cè)為南風(fēng)且出現(xiàn)大量越赤道氣流(圖7d)。所以,觀測結(jié)果驗(yàn)證了理想熱源激發(fā)的大氣響應(yīng),雖然實(shí)際大氣分布更復(fù)雜,但我們可以看到分布型相似的實(shí)際赤道熱源對應(yīng)的熱帶大氣與理想熱源得到的大氣響應(yīng)的分布型幾乎一致。因此,在孟加拉灣地區(qū)氣旋環(huán)流及越赤道氣流的形成過程中,海洋加熱的強(qiáng)迫將起著重要的作用。

    圖7 孟加拉灣地區(qū)(a)3 月、(b)4月海溫(單位:°C)與風(fēng)場(單位:m/s)分布以及分別與(c)3月、(d)4月海溫分布型一致的理想熱源及其大氣響應(yīng)分布圖。(a)和(b)中黃色實(shí)線代表經(jīng)向海溫最大值,(c)和(d)中紅色實(shí)線代表熱源,等值線間隔是0.3

    5 總結(jié)與討論

    本文采用解析解得到了單一型赤道非對稱熱源的結(jié)果,并分析了不同熱源分布型下大氣的響應(yīng),從動(dòng)力學(xué)角度闡述了非對稱加熱對季風(fēng)爆發(fā)的重要作用。主要結(jié)論如下:

    (1)我們將單一型赤道非對稱熱源做Weber函數(shù)展開,然后將熱源展開項(xiàng)逐一代入Gill模式進(jìn)行求解。對于>1的熱源展開項(xiàng),解析方程存在唯一解,所以我們得到了赤道非對稱熱源的解析通解。該通解從理論上完善了單一型赤道非對稱熱源激發(fā)的赤道非對稱的大氣響應(yīng)結(jié)果。

    (2)熱源的強(qiáng)度、位置及范圍變化對結(jié)果都具有一定的影響。當(dāng)熱源強(qiáng)度增強(qiáng)(減弱),熱源激發(fā)的大氣響應(yīng)整體增強(qiáng)(減弱),但大氣分布型不發(fā)生變化;當(dāng)熱源中心位置不斷北移,北半球氣旋強(qiáng)度增加、位置北移,同時(shí)赤道輻合氣流減弱而越赤道氣流增強(qiáng);當(dāng)熱源范圍不斷增大時(shí),北半球氣旋強(qiáng)度增強(qiáng)、位置西移、范圍增大,同時(shí)越赤道氣流增強(qiáng)。

    (3)大范圍熱源的季節(jié)性移動(dòng)與夏季風(fēng)的爆發(fā)密切相關(guān),而熱帶海溫分布與熱源的季節(jié)移動(dòng)有密切聯(lián)系。本文以孟加拉灣地區(qū)海溫對夏季風(fēng)爆發(fā)的影響作為個(gè)例進(jìn)行研究。孟加拉灣地區(qū)3月份最暖海溫位于赤道附近,4月份最暖海溫北移至10°N,所對應(yīng)的大氣響應(yīng)如下:當(dāng)最大海溫位于赤道附近時(shí),其兩側(cè)表現(xiàn)為Rossby波響應(yīng)的氣旋環(huán)流,而當(dāng)海溫暖軸北移后,其南側(cè)激發(fā)出大范圍越赤道氣流,這是大氣對Rossby重力混合波響應(yīng)的結(jié)果。觀測結(jié)果與分布型相似的理想熱源的大氣響應(yīng)一致,同時(shí)也證明了海洋加熱的強(qiáng)迫對大氣環(huán)流分布型起著重要的作用。

    本文利用Gill模式分析了赤道非對稱熱源的大氣響應(yīng),但是仍有一些問題有待進(jìn)一步的研究。由于熱源分布的多樣性及復(fù)雜性,如兩類厄爾尼諾事件、拉尼娜事件中海溫分布型以及印度洋偶極子型等(Wallace et al., 1998; Saji et al., 1999; Ashok et al., 2007),這些復(fù)合型熱源的大氣響應(yīng)特征有待于進(jìn)一步的分析。

    致謝 感謝兩位審稿專家和編輯老師的建議,感謝課題組同學(xué)提供的幫助和有益的討論。

    Ashok K, Behera S K, Rao S A, et al. 2007. El Ni?o Modoki and its possible teleconnection [J]. J. Geophys. Res., 112, C11007, doi: 10.1029/2006JC003798.

    Bjerknes J. 1966. A possible response of the atmospheric Hadley circulation to equatorial anomalies of ocean temperature [J]. Tellus, 18 (4): 820–829.

    Bjerknes J. 1969. Atmospheric teleconnections from the equatorial Pacific [J]. Mon. Wea. Rev., 97 (3): 163–172.

    Bretherton C S, Sobel A H. 2003. The Gill model and the weak temperature gradient approximation [J]. J. Atmos. Sci., 60 (2): 451–460.

    巢紀(jì)平, 王彰貴. 1991. 海洋加熱尺度對熱帶大氣垂直環(huán)流圈結(jié)構(gòu)的影響 [J]. 南京氣象學(xué)院學(xué)報(bào), 14 (1): 10–17. Chao Jiping, Wang Zhanggui. 1991. Influence of the ocean heating scale on the structure of the vertical circulation cell in the tropical atmosphere [J]. Journal of Nanjing Institute of Meteorology (in Chinese), 14 (1): 10–17.

    陳哲, 李崇銀. 2006. 亞洲夏季風(fēng)爆發(fā)與熱源強(qiáng)迫下的熱帶大氣Rossby波 [J]. 大氣科學(xué), 30 (6): 1227–1235. Chen Zhe, Li Chongyin. 2006. Asian summer monsoon onset and tropical atmosphere Rossby wave forced by the heating source [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 30 (6): 1227–1235.

    Feng J, Li J P. 2013. Contrasting impacts of two types of ENSO on the boreal spring Hadley circulation [J]. J. Climate, 26: 4773–4789.

    Feng J, Li J P, Xie F. 2013. Long-term variation of the principal mode of boreal spring Hadley circulation linked to SST over the Indo-Pacific warm pool [J]. J. Climate, 26: 532–544.

    Feng R, Li J P, Wang J C. 2011. Regime change of the boreal summer Hadley circulation and its connection with the tropical SST [J]. J. Climate, 24: 3867–3877.

    Gill A E. 1975. Models of equatorial currents in numerical models of ocean circulation [C]// Proceedings of a Symposium in Durham, N. H. Washington: National Academy of Sciences, 181–203.

    Gill A E. 1980. Some simple solutions for heat-induced tropical circulation [J]. Quart. J. Roy. Meteor. Soc., 106 (449): 447–462.

    Guan Z Y, Ashok K, Yamagata T. 2003. Summertime response of the tropical atmosphere to the Indian Ocean dipole sea surface temperature anomalies [J]. J. Meteor. Soc. Japan, 81 (3): 533–561.

    Hoskins B J, Karoly D J. 1981. The steady linear response of a spherical atmosphere to thermal and orographic forcing [J]. J. Atmos. Sci., 38 (6): 1179–1196.

    黃亮, 譚桂容, 劉向文. 2012. 熱力強(qiáng)迫對澳洲北部越赤道氣流強(qiáng)弱變化的影響 [J]. 大氣科學(xué)學(xué)報(bào), 35 (3): 311–321. Huang Liang, Tan Guirong, Liu Xiangwen. 2012. The impacts of thermal forcing on the cross-equatorial flow over North Australia [J]. Transactions of Atmospheric Sciences (in Chinese), 35 (3): 311–321.

    Jiang X W, Li J P. 2011. Influence of the annual cycle of sea surface temperature on the monsoon onset [J]. J. Geophys. Res., 116 (D10), doi: 10.1029/2010jd015236.

    Jin F F, Hoskins B J. 1995. The direct response to tropical heating in a baroclinic atmosphere [J]. J. Atmos. Sci., 52 (3): 307–319.

    Kalnay E, Kanamitsu M, Kistler R, et al. 1996. The NCEP/NCAR 40-year reanalysis project [J]. Bull. Amer. Meteor. Soc., 77 (3): 437–471.

    Latif M, Anderson D, Barnett T, et al. 1998. A review of the predictability and prediction of ENSO [J]. J. Geophys. Res., 103 (C7): 14375–14393.

    李建平, 吳國雄, 胡敦欣, 等. 2011. 亞印太交匯區(qū)海氣相互作用及其對我國短期氣候的影響(上卷) [M]. 北京: 氣象出版社, 1–516. Li Jianping, Wu Guoxiong, Hu Dunxin, et al. 2011. Ocean?Atmosphere Interaction over the Joining Area of Asia and Indian?Pacific Ocean and Its Impact on the Short-term Climate Variation in China (Volume 1) (in Chinese) [M]. Beijing: China Meteorological Press, 1–516.

    劉琳, 于衛(wèi)東. 2006. ENSO循環(huán)之中的海氣相互響應(yīng)——對次表層海溫和大氣旋度的資料分析 [J]. 地球物理學(xué)報(bào), 49 (1): 45–51. Liu Lin, Yu Weidong. 2006. Analysis of the characteristic time scale during ENSO [J]. Chinese Journal of Geophysics (in Chinese), 49 (1): 45–51.

    劉秦玉, 王韶霞. 2000. 印度尼西亞群島鄰近海域SST對越赤道氣流的影響 [J]. 大氣科學(xué), 24 (5): 634–640. Liu Qinyu, Wang Shaoxia. 2000. Influence of the sea surface temperature in the vicinity of the Indonesia Archipelago on the cross equator flow [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 24 (5): 634–640.

    Liu Y M, Chan J C L, Mao J Y, et al. 2002. The role of Bay of Bengal convection in the onset of the 1998 South China Sea summer monsoon [J]. Mon. Wea. Rev., 130 (11): 2731–2744.

    Lorenz E N. 1955. Available potential energy and the maintenance of the general circulation [J]. Tellus, 7 (2): 157–167.

    Matsuno T. 1966. Quasi-geostrophic motions in the equatorial area [J]. J. Meteor. Soc. Japan, 44 (1): 25–43.

    Ren B H, Huang R H. 2002. 10-25-Day intraseasonal variations of convection and circulation associated with thermal state of the western Pacific warm pool during boreal summer [J]. Advances in Atmospheric Sciences, 19 (2): 321–336.

    Reynolds R W, Rayner N A, Smith T M, et al. 2002. An improved in situ and satellite SST analysis for climate [J]. J. Climate, 15 (13): 1609– 1625.

    Saji N, Goswami B N, Vinayachandran P, et al. 1999. A dipole mode in the tropical Indian Ocean [J]. Nature, 401 (6751): 360–363.

    Sashegyi K D, Geisler J E. 1987. A linear model study of cross-equatorial flow forced by summer monsoon heat sources [J]. J. Atmos. Sci., 44 (13): 1706?1722.

    唐東昇, 謝立安. 1986. 加熱場對低空越赤道氣流影響的數(shù)值試驗(yàn)研究 [J]. 熱帶氣象, 2 (4): 311–319. Tang Dongsheng, Xie Li’an. 1986. Numerical experiments of the effects of the heating field on the low-level cross-equatorial flows [J]. Journal of Tropical Meteorology (in Chinese), 2 (4): 311–319.

    Wallace J M, Rasmusson E M, Mitchell T P, et al. 1998. On the structure and evolution of ENSO-related climate variability in the tropical Pacific: Lessons from TOGA [J]. J. Geophys. Res., 103 (C7): 14241–14259.

    Webster P J. 1972. Response of the tropical atmosphere to local, steady forcing [J]. Mon. Wea. Rev., 100 (7): 518–541.

    Webster P J. 1981. Mechanisms determining the atmospheric response to sea surface temperature anomalies [J]. J. Atmos. Sci., 38 (3): 554–571.

    吳秋英, 朱乾根. 1990. 熱帶大氣環(huán)流對低緯太平洋SST暖異常的響應(yīng) [J]. 南京氣象學(xué)院學(xué)報(bào), 13 (1): 11–22. Wu Qiuying, Zhu Qian’gen. 1990. Response of the tropical atmospheric circulation to the low-latitude Pacific SST warm anomaly [J]. Journal of Nanjing Institute of Meteorology (in Chinese), 13 (1): 11–22.

    薛洪斌, 張銘, 王業(yè)桂. 2006. 熱帶太平洋SST異常對IAP-9 LAGCM年際變率影響的模擬 [J]. 氣象科學(xué), 26 (1): 58–65. Xue Hongbin, Zhang Ming, Wang Yegui. 2006. The simulation of impact on climate interannual variability due to SST variation in the tropics [J]. Scientia Meteorologica Sinica (in Chinese), 26 (1): 58–65.

    Yu W D, Shi J W, Liu L, et al. 2012. The onset of the monsoon over the Bay of Bengal: The observed common features for 2008–2011 [J]. Atmos. Oceanic Sci. Lett., 5 (4): 314–318.

    邢楠, 李建平, 李耀錕. 2014. 熱帶大氣對單一型赤道非對稱熱源的響應(yīng)[J]. 大氣科學(xué), 38 (6): 1147?1158, doi:10.3878/j.issn.1006-9895.1401.13275. Xing Nan, Li Jianping, Li Yaokun. 2014. Response of the tropical atmosphere to isolated equatorially asymmetric heating [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 38 (6): 1147?1158.

    Response of the Tropical Atmosphere to Isolated Equatorially Asymmetric Heating

    XING Nan, LI Jianping, and LI Yaokun

    1,,1000292,1000493,,100875

    An analysis of the response of the tropical atmosphere to isolated equatorially asymmetric heating is performed using the Gill model. General solutions of the Gill model are calculated, and the equatorially asymmetric atmospheric response is examined in detail within the general solutions. Heating strength, location, and extent all have some influence on the atmospheric response. When the center of heating moves northward, the cyclone in the Northern Hemisphere is strengthened and moves northward; at the same time, equatorial convergence weakens and cross-equatorial flow is enhanced. When heating strength increases (decreases), the atmospheric response forced by the heating is enhanced (weakened), however, the atmospheric pattern does not vary. When the heating extent increases, the cyclonic strength and extent are enhanced, the cyclone is displaced westward, and cross-equatorial flow is enhanced. The above theoretical conclusions are applied to study of the influence of SST on the onset of monsoons in the Bay of Bengal. When the meridional warmest SST axis (WSSTA) is located at the equator, the Gill pattern appears on either side of the equator; however, when the WSSTA moves northward, cyclonic circulation disappears and cross-equatorial flow is induced south of the WSSTA. This effect represents the response of mixed Rossby-gravity waves to external heat forcing. This work shows the utility of the above analytic approach, and also explains the dynamics responsible for summer monsoon onset in the Bay of Bengal.

    Isolated equatorially asymmetric heating, Heat forcing, Summer monsoon onset in the Bay of Bengal, Sea surface temperature

    1006?9895(2014)06?1147?12

    P433

    A

    10.3878/j.issn.1006-9895.1401.13275

    2013?10?08,2014?01?24 收修定稿

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973計(jì)劃)項(xiàng)目2010CB950400,國家自然科學(xué)基金重點(diǎn)項(xiàng)目41030961,國家海洋局項(xiàng)目

    邢楠,女,1988年出生,博士研究生,主要從事季風(fēng)與海氣相互作用研究。E-mail: nxing0923@163.com

    李建平,E-mail: ljp@lasg.iap.ac.cn

    猜你喜歡
    海溫氣旋赤道
    “雪龍2”號過赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    2003年3月北大西洋上兩個(gè)爆發(fā)性氣旋的“吞并”過程及發(fā)展機(jī)制分析?
    基于深度學(xué)習(xí)的海溫觀測數(shù)據(jù)質(zhì)量控制應(yīng)用研究
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    氨還原劑噴槍中保護(hù)氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關(guān)系
    距赤道越遠(yuǎn)越喜愛黃色
    最熱的地方不在赤道
    北太平洋上一個(gè)爆發(fā)性氣旋族的結(jié)構(gòu)分析?
    2014年3月大西洋上一個(gè)爆發(fā)性氣旋的研究
    看免费av毛片| 国产成人精品久久二区二区91| 日本成人三级电影网站| 中文字幕精品亚洲无线码一区 | 国内久久婷婷六月综合欲色啪| 日韩欧美一区二区三区在线观看| 日韩中文字幕欧美一区二区| 亚洲第一欧美日韩一区二区三区| 国产一区二区三区在线臀色熟女| 日本免费a在线| 精品国产美女av久久久久小说| 亚洲国产欧美日韩在线播放| 每晚都被弄得嗷嗷叫到高潮| 少妇熟女aⅴ在线视频| 色老头精品视频在线观看| 亚洲专区字幕在线| 黄网站色视频无遮挡免费观看| 女同久久另类99精品国产91| 欧美 亚洲 国产 日韩一| 久久国产精品人妻蜜桃| 少妇裸体淫交视频免费看高清 | av电影中文网址| 法律面前人人平等表现在哪些方面| 一本精品99久久精品77| 精品久久蜜臀av无| 欧美成狂野欧美在线观看| 日本 欧美在线| 日韩欧美免费精品| 色精品久久人妻99蜜桃| 久久中文字幕一级| 搡老岳熟女国产| 国产成人系列免费观看| 国内揄拍国产精品人妻在线 | 久久九九热精品免费| 一本一本综合久久| 日本撒尿小便嘘嘘汇集6| 久久久久久亚洲精品国产蜜桃av| 欧美成狂野欧美在线观看| 久99久视频精品免费| 日本一区二区免费在线视频| 亚洲av电影在线进入| 一个人观看的视频www高清免费观看 | 自线自在国产av| 久久香蕉精品热| 黄色丝袜av网址大全| 伦理电影免费视频| 久久精品国产综合久久久| 在线十欧美十亚洲十日本专区| 成年版毛片免费区| 欧美乱码精品一区二区三区| 在线国产一区二区在线| 亚洲国产精品合色在线| 亚洲精品一卡2卡三卡4卡5卡| 精品国产超薄肉色丝袜足j| 精品电影一区二区在线| 色在线成人网| 国产精品二区激情视频| 欧美性猛交黑人性爽| 亚洲av熟女| 成年人黄色毛片网站| 亚洲成a人片在线一区二区| 欧美一级毛片孕妇| 亚洲 欧美一区二区三区| 久久精品亚洲精品国产色婷小说| 欧美av亚洲av综合av国产av| 久久中文字幕一级| 国产精品 欧美亚洲| 国产乱人伦免费视频| 草草在线视频免费看| 男女视频在线观看网站免费 | 好看av亚洲va欧美ⅴa在| 神马国产精品三级电影在线观看 | 99久久精品国产亚洲精品| 午夜日韩欧美国产| 巨乳人妻的诱惑在线观看| 免费女性裸体啪啪无遮挡网站| 最新在线观看一区二区三区| 精品卡一卡二卡四卡免费| 午夜免费成人在线视频| 人成视频在线观看免费观看| 成人18禁高潮啪啪吃奶动态图| 又紧又爽又黄一区二区| 天堂√8在线中文| 又大又爽又粗| 成人精品一区二区免费| 久久人人精品亚洲av| 俄罗斯特黄特色一大片| www.精华液| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| 欧美黄色淫秽网站| 白带黄色成豆腐渣| 亚洲精品美女久久久久99蜜臀| 后天国语完整版免费观看| 看免费av毛片| 亚洲国产精品sss在线观看| 性色av乱码一区二区三区2| 久久精品国产99精品国产亚洲性色| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆 | 精品卡一卡二卡四卡免费| 免费看a级黄色片| 国产黄色小视频在线观看| 特大巨黑吊av在线直播 | 亚洲色图 男人天堂 中文字幕| 欧美一级a爱片免费观看看 | 亚洲精品美女久久av网站| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 少妇粗大呻吟视频| 天天添夜夜摸| 午夜成年电影在线免费观看| 国产亚洲精品久久久久久毛片| 18禁国产床啪视频网站| 精品熟女少妇八av免费久了| 女性被躁到高潮视频| 久久精品夜夜夜夜夜久久蜜豆 | 国产伦人伦偷精品视频| 国产精品永久免费网站| 天天躁夜夜躁狠狠躁躁| 色播亚洲综合网| 国产高清videossex| 国产午夜福利久久久久久| 19禁男女啪啪无遮挡网站| 久久这里只有精品19| 午夜日韩欧美国产| 国产精品久久久久久精品电影 | 嫩草影视91久久| 久久中文字幕人妻熟女| 日韩欧美国产在线观看| 亚洲av成人不卡在线观看播放网| 人成视频在线观看免费观看| 久久 成人 亚洲| 2021天堂中文幕一二区在线观 | 不卡一级毛片| 欧美日韩乱码在线| 黄色丝袜av网址大全| 亚洲精品中文字幕在线视频| 国产精品日韩av在线免费观看| 国产精品久久视频播放| 成熟少妇高潮喷水视频| 成人永久免费在线观看视频| 亚洲精品国产区一区二| 国产高清视频在线播放一区| 久久99热这里只有精品18| 亚洲九九香蕉| 99riav亚洲国产免费| 狂野欧美激情性xxxx| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 午夜影院日韩av| 久久久水蜜桃国产精品网| 亚洲黑人精品在线| 亚洲免费av在线视频| 我的亚洲天堂| 亚洲精品国产区一区二| 免费在线观看完整版高清| 在线观看免费午夜福利视频| 免费看日本二区| 亚洲成人久久爱视频| 日日爽夜夜爽网站| 日本 av在线| tocl精华| 亚洲va日本ⅴa欧美va伊人久久| 国产av一区二区精品久久| 99国产精品一区二区蜜桃av| 成年免费大片在线观看| 人妻丰满熟妇av一区二区三区| 欧美不卡视频在线免费观看 | 久久久久久九九精品二区国产 | 亚洲,欧美精品.| 国产成人精品久久二区二区91| 美女 人体艺术 gogo| 久久久久久久久免费视频了| 无限看片的www在线观看| 自线自在国产av| 日韩三级视频一区二区三区| 久久天躁狠狠躁夜夜2o2o| www国产在线视频色| 99国产精品99久久久久| 麻豆久久精品国产亚洲av| 日韩欧美一区视频在线观看| 俄罗斯特黄特色一大片| 99久久国产精品久久久| 国产亚洲精品一区二区www| 在线观看一区二区三区| 午夜精品在线福利| 婷婷亚洲欧美| 国产欧美日韩精品亚洲av| 精品欧美国产一区二区三| 久久精品国产综合久久久| 草草在线视频免费看| 可以在线观看的亚洲视频| 亚洲最大成人中文| 久久久久久国产a免费观看| 成人亚洲精品av一区二区| 亚洲一区二区三区色噜噜| 亚洲九九香蕉| 91大片在线观看| av中文乱码字幕在线| 久久天躁狠狠躁夜夜2o2o| 黄色毛片三级朝国网站| 巨乳人妻的诱惑在线观看| 一本一本综合久久| 久久久久国产精品人妻aⅴ院| 亚洲精品色激情综合| 女人被狂操c到高潮| 熟女电影av网| 狠狠狠狠99中文字幕| 91老司机精品| 久久狼人影院| 日韩一卡2卡3卡4卡2021年| 欧美精品亚洲一区二区| 久久中文字幕一级| 国产激情偷乱视频一区二区| 他把我摸到了高潮在线观看| 夜夜夜夜夜久久久久| 国产亚洲精品久久久久久毛片| 欧美性猛交╳xxx乱大交人| 国产高清videossex| 日韩精品青青久久久久久| 一个人观看的视频www高清免费观看 | 免费高清在线观看日韩| 50天的宝宝边吃奶边哭怎么回事| 国产在线观看jvid| 18禁黄网站禁片午夜丰满| 欧美日韩精品网址| 嫩草影院精品99| 国产人伦9x9x在线观看| 午夜福利欧美成人| 可以在线观看的亚洲视频| 婷婷精品国产亚洲av| 露出奶头的视频| 亚洲欧洲精品一区二区精品久久久| 亚洲真实伦在线观看| 桃红色精品国产亚洲av| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 精品一区二区三区av网在线观看| 国产在线观看jvid| 亚洲成a人片在线一区二区| 亚洲三区欧美一区| 人成视频在线观看免费观看| 国产亚洲精品av在线| 1024香蕉在线观看| 99热这里只有精品一区 | 亚洲av片天天在线观看| 亚洲成人国产一区在线观看| 少妇的丰满在线观看| 欧美亚洲日本最大视频资源| 丝袜人妻中文字幕| 亚洲自偷自拍图片 自拍| 天堂√8在线中文| 国产久久久一区二区三区| 亚洲真实伦在线观看| 成年人黄色毛片网站| www日本在线高清视频| 亚洲熟妇熟女久久| 淫妇啪啪啪对白视频| 国产又爽黄色视频| 热99re8久久精品国产| 国产精品亚洲一级av第二区| 日韩av在线大香蕉| 欧美精品啪啪一区二区三区| 香蕉丝袜av| 高潮久久久久久久久久久不卡| 波多野结衣av一区二区av| av福利片在线| 久久婷婷成人综合色麻豆| 欧美在线一区亚洲| 超碰成人久久| 黑丝袜美女国产一区| 高潮久久久久久久久久久不卡| 久久亚洲精品不卡| 91麻豆av在线| 亚洲自拍偷在线| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 亚洲国产欧美一区二区综合| 天天一区二区日本电影三级| 免费观看人在逋| 一级毛片精品| 久久精品影院6| 伦理电影免费视频| 黄色女人牲交| svipshipincom国产片| 国产亚洲精品第一综合不卡| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 一级片免费观看大全| 无限看片的www在线观看| 成人免费观看视频高清| АⅤ资源中文在线天堂| 最新美女视频免费是黄的| 一级a爱片免费观看的视频| 亚洲成人久久爱视频| 动漫黄色视频在线观看| 一本久久中文字幕| 久久国产精品影院| 亚洲五月天丁香| 超碰成人久久| 搡老妇女老女人老熟妇| 亚洲专区字幕在线| 免费观看人在逋| 国产激情偷乱视频一区二区| 18禁国产床啪视频网站| 精品熟女少妇八av免费久了| 老司机深夜福利视频在线观看| 女人高潮潮喷娇喘18禁视频| 日本一区二区免费在线视频| 国内少妇人妻偷人精品xxx网站 | 69av精品久久久久久| 看黄色毛片网站| 国产精品久久久人人做人人爽| 欧美黄色淫秽网站| 久久国产乱子伦精品免费另类| 热99re8久久精品国产| 一边摸一边做爽爽视频免费| 亚洲精品国产精品久久久不卡| 最新在线观看一区二区三区| 精品人妻1区二区| 一级作爱视频免费观看| 精品一区二区三区视频在线观看免费| 国产精华一区二区三区| 日韩中文字幕欧美一区二区| 男人的好看免费观看在线视频 | 国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 久久婷婷人人爽人人干人人爱| 中出人妻视频一区二区| 亚洲第一av免费看| 欧美日韩亚洲国产一区二区在线观看| 国产不卡一卡二| 成年女人毛片免费观看观看9| 99国产极品粉嫩在线观看| 日韩有码中文字幕| 男人舔女人的私密视频| 久久这里只有精品19| 久久久久久九九精品二区国产 | 日日摸夜夜添夜夜添小说| 男女床上黄色一级片免费看| xxx96com| 国产欧美日韩精品亚洲av| 日本一本二区三区精品| 亚洲av美国av| 免费观看精品视频网站| 日韩三级视频一区二区三区| 亚洲天堂国产精品一区在线| 亚洲,欧美精品.| 日韩欧美国产在线观看| 99热只有精品国产| 淫妇啪啪啪对白视频| 香蕉av资源在线| 日本在线视频免费播放| 香蕉av资源在线| 久久久久精品国产欧美久久久| 日本 av在线| 精品久久久久久成人av| 国产精品影院久久| 国产日本99.免费观看| 神马国产精品三级电影在线观看 | 久久国产亚洲av麻豆专区| 亚洲精品美女久久久久99蜜臀| 国产成人av激情在线播放| 亚洲男人天堂网一区| 最近最新中文字幕大全电影3 | 亚洲欧美精品综合久久99| 波多野结衣av一区二区av| 在线av久久热| 欧美在线黄色| 日本三级黄在线观看| av在线播放免费不卡| 成熟少妇高潮喷水视频| 久久99热这里只有精品18| 高清在线国产一区| 深夜精品福利| 最新美女视频免费是黄的| 欧美激情 高清一区二区三区| 正在播放国产对白刺激| 国产精品二区激情视频| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 国产99久久九九免费精品| 欧美激情极品国产一区二区三区| 亚洲中文日韩欧美视频| 黑丝袜美女国产一区| 久久久久免费精品人妻一区二区 | 黑人欧美特级aaaaaa片| 亚洲国产精品999在线| 男女视频在线观看网站免费 | 99久久久亚洲精品蜜臀av| 国产一区二区激情短视频| 女人被狂操c到高潮| 日本黄色视频三级网站网址| 久久人妻福利社区极品人妻图片| 色综合欧美亚洲国产小说| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 精品久久久久久久人妻蜜臀av| 久久婷婷人人爽人人干人人爱| 午夜激情福利司机影院| 亚洲三区欧美一区| 亚洲av第一区精品v没综合| 可以在线观看毛片的网站| 久久婷婷人人爽人人干人人爱| 欧美色欧美亚洲另类二区| 成人特级黄色片久久久久久久| 日日干狠狠操夜夜爽| 国产激情偷乱视频一区二区| 久久中文看片网| 色av中文字幕| 亚洲av成人不卡在线观看播放网| 亚洲人成77777在线视频| 国产成人欧美在线观看| 日韩精品青青久久久久久| 又大又爽又粗| 国产乱人伦免费视频| 18美女黄网站色大片免费观看| 侵犯人妻中文字幕一二三四区| 男人的好看免费观看在线视频 | 变态另类丝袜制服| 91九色精品人成在线观看| 午夜亚洲福利在线播放| 亚洲专区国产一区二区| 久久九九热精品免费| 久久久久久九九精品二区国产 | 青草久久国产| 好看av亚洲va欧美ⅴa在| 18禁美女被吸乳视频| 无限看片的www在线观看| cao死你这个sao货| 在线观看www视频免费| 在线观看免费午夜福利视频| 精品久久久久久,| 国产激情久久老熟女| 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 91成年电影在线观看| 欧美绝顶高潮抽搐喷水| 亚洲国产欧美一区二区综合| 不卡一级毛片| 两人在一起打扑克的视频| 日韩 欧美 亚洲 中文字幕| 伦理电影免费视频| 国产三级黄色录像| 天天一区二区日本电影三级| 俺也久久电影网| 国产爱豆传媒在线观看 | 久久精品亚洲精品国产色婷小说| 国产精品 国内视频| 亚洲精品在线美女| 女生性感内裤真人,穿戴方法视频| 国产1区2区3区精品| 中文字幕另类日韩欧美亚洲嫩草| 婷婷精品国产亚洲av| 午夜福利视频1000在线观看| 国产爱豆传媒在线观看 | 日本熟妇午夜| www.www免费av| 国产成人欧美| 国产97色在线日韩免费| 久热这里只有精品99| 人人澡人人妻人| 国产一区二区在线av高清观看| 国产熟女xx| 亚洲性夜色夜夜综合| 又黄又爽又免费观看的视频| 大香蕉久久成人网| 此物有八面人人有两片| 美女扒开内裤让男人捅视频| 免费高清视频大片| 韩国av一区二区三区四区| 日韩av在线大香蕉| 老司机靠b影院| 国产亚洲欧美98| 亚洲精品久久成人aⅴ小说| 淫妇啪啪啪对白视频| 国产精品自产拍在线观看55亚洲| 亚洲午夜精品一区,二区,三区| 禁无遮挡网站| 男女午夜视频在线观看| 天堂影院成人在线观看| 视频区欧美日本亚洲| 国产精品乱码一区二三区的特点| 桃红色精品国产亚洲av| 国产av在哪里看| 日韩欧美在线二视频| 亚洲熟妇熟女久久| 亚洲最大成人中文| 亚洲色图 男人天堂 中文字幕| 人人妻人人澡欧美一区二区| 伦理电影免费视频| 国产精品爽爽va在线观看网站 | 免费看美女性在线毛片视频| 一二三四在线观看免费中文在| 国产精品久久视频播放| 女同久久另类99精品国产91| 国产亚洲精品久久久久5区| 日韩av在线大香蕉| 脱女人内裤的视频| 美女高潮喷水抽搐中文字幕| av片东京热男人的天堂| 黑人巨大精品欧美一区二区mp4| 成在线人永久免费视频| 人妻丰满熟妇av一区二区三区| 国产高清激情床上av| 亚洲欧美精品综合久久99| 搡老岳熟女国产| 欧美黄色淫秽网站| 2021天堂中文幕一二区在线观 | 黑人巨大精品欧美一区二区mp4| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 欧美成人性av电影在线观看| 搞女人的毛片| 日本免费a在线| 一本大道久久a久久精品| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 午夜福利免费观看在线| 久久久久久久久中文| 免费高清在线观看日韩| 亚洲午夜理论影院| 国产精品精品国产色婷婷| 国产亚洲精品一区二区www| 国产精品av久久久久免费| 91在线观看av| 黑丝袜美女国产一区| 亚洲一区二区三区不卡视频| 久久中文字幕一级| 777久久人妻少妇嫩草av网站| 男人的好看免费观看在线视频 | 精品久久久久久,| 精品第一国产精品| 午夜亚洲福利在线播放| 身体一侧抽搐| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区在线av高清观看| 色老头精品视频在线观看| 99精品久久久久人妻精品| 精品国产一区二区三区四区第35| 国产精品电影一区二区三区| 人成视频在线观看免费观看| 大型黄色视频在线免费观看| 亚洲第一电影网av| 亚洲狠狠婷婷综合久久图片| 国产亚洲av高清不卡| 桃色一区二区三区在线观看| 18禁国产床啪视频网站| 成人永久免费在线观看视频| 制服诱惑二区| 久久国产精品男人的天堂亚洲| 啦啦啦韩国在线观看视频| 69av精品久久久久久| 亚洲全国av大片| 欧美中文综合在线视频| or卡值多少钱| 午夜两性在线视频| 免费看日本二区| 在线天堂中文资源库| 欧美激情久久久久久爽电影| 一进一出抽搐动态| 欧美午夜高清在线| 搡老熟女国产l中国老女人| 亚洲五月婷婷丁香| 巨乳人妻的诱惑在线观看| 国产激情欧美一区二区| 男女视频在线观看网站免费 | 欧美日韩亚洲综合一区二区三区_| 在线观看www视频免费| 欧美成人免费av一区二区三区| 国产av不卡久久| 午夜激情av网站| 精品无人区乱码1区二区| 亚洲,欧美精品.| 日本熟妇午夜| bbb黄色大片| 国产私拍福利视频在线观看| 老司机深夜福利视频在线观看| 一级作爱视频免费观看| 午夜福利一区二区在线看| 香蕉久久夜色| 最近在线观看免费完整版| 日韩欧美一区视频在线观看| 在线观看舔阴道视频| 成人手机av| 18禁黄网站禁片免费观看直播| 久久 成人 亚洲| 又大又爽又粗| 亚洲第一av免费看| 99国产精品99久久久久| 嫁个100分男人电影在线观看| 国产又爽黄色视频| 热re99久久国产66热| 国产亚洲av嫩草精品影院| 久久狼人影院| 黑人操中国人逼视频| 真人一进一出gif抽搐免费| 一区二区日韩欧美中文字幕| 美国免费a级毛片| 国产真实乱freesex| 久久精品aⅴ一区二区三区四区| 韩国av一区二区三区四区| 好男人在线观看高清免费视频 | 国产1区2区3区精品| 桃红色精品国产亚洲av| 国产不卡一卡二| 欧美日本视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产蜜桃级精品一区二区三区| 国产免费男女视频| 操出白浆在线播放| 国产精品免费一区二区三区在线| www日本在线高清视频| 亚洲成人久久爱视频|