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

    北半球大氣對(duì)春季北極海冰異常響應(yīng)的數(shù)值模擬

    2011-12-15 07:06:10張若楠武炳義
    大氣科學(xué) 2011年5期
    關(guān)鍵詞:斜壓北極海位勢(shì)

    張若楠 武炳義

    中國(guó)氣象科學(xué)研究院,北京 100081

    北半球大氣對(duì)春季北極海冰異常響應(yīng)的數(shù)值模擬

    張若楠 武炳義

    中國(guó)氣象科學(xué)研究院,北京 100081

    在CAM 3.0模式中,通過(guò)設(shè)計(jì)一系列數(shù)值模擬試驗(yàn)來(lái)研究北半球大氣環(huán)流對(duì)北極海冰密集度 (Sea Ice Concentration,簡(jiǎn)稱(chēng)SIC)強(qiáng)迫的響應(yīng)。將海冰密集度 EOF第二模態(tài) (EOF2)的時(shí)間系數(shù)回歸到原場(chǎng),并把這個(gè)空間異常場(chǎng)逐月的加到SIC春季氣候態(tài)上,作為敏感性試驗(yàn)的外強(qiáng)迫。試驗(yàn)結(jié)果表明:在季節(jié)尺度上,大氣環(huán)流由春季逐漸向夏季演變,春季SIC異常與春、夏季北半球大氣環(huán)流、地表溫度以及降水異常有顯著關(guān)系,尤其是對(duì)春季中國(guó)東部降水負(fù)異常,夏季東北部、東南沿海降水正異常有很大貢獻(xiàn)。在逐日和逐周時(shí)間尺度瞬態(tài)演變上,局地高、低層大氣環(huán)流對(duì)海冰強(qiáng)迫的響應(yīng)呈現(xiàn)斜壓結(jié)構(gòu),兩周之后大氣由斜壓結(jié)構(gòu)逐漸轉(zhuǎn)變?yōu)檎龎航Y(jié)構(gòu),強(qiáng)迫響應(yīng)可以通過(guò)遙相關(guān)波列傳播到北極以外地區(qū),到了第6周大氣的響應(yīng)基本達(dá)到正壓平衡狀態(tài)。這種遙相關(guān)的產(chǎn)生與極區(qū)熱通量異常所激發(fā)的波活動(dòng)通量有關(guān):北極SIC異常首先通過(guò)直接熱力強(qiáng)迫過(guò)程來(lái)改變表層熱通量空間分布,這種表面熱通量異常通過(guò)與大氣環(huán)流的相互作用可以激發(fā)出大尺度羅斯貝波。在初始斜壓區(qū)的低層有波動(dòng)能量向上傳播,并在對(duì)流層中高層向南傳播,最后通過(guò)直接熱力強(qiáng)迫和大氣內(nèi)部動(dòng)力學(xué)相互作用引發(fā)的遙相關(guān)過(guò)程將能量頻散到東亞地區(qū),進(jìn)而影響東亞地區(qū)的天氣和氣候。

    海冰 海冰密集度 夏季降水 敏感性試驗(yàn) 控制試驗(yàn) 瞬變過(guò)程 波活動(dòng)通量

    1 引言

    海冰是氣候系統(tǒng)五大圈層的重要組成部分,它的高反照率和隔絕海—?dú)鉄崃亢蛣?dòng)量交換作用對(duì)區(qū)域乃至全球氣候變化均有重要的調(diào)節(jié)作用 (Walsh,1983;彭公炳等,1992)。在全球變暖過(guò)程中,氣溫變化幅度最大的是北半球的中高緯地區(qū),這與北極海冰減少對(duì)氣溫增幅的放大作用有密切的關(guān)系(Screen and Simmonds,2010)。海冰自身變化必然對(duì)氣候變化產(chǎn)生重要的反饋?zhàn)饔?。首?海冰的存在顯著地改變了海洋和大氣之間的熱量、動(dòng)量和物質(zhì)交換。其次,海冰的高反照率減少了極地對(duì)太陽(yáng)輻射的吸收,使得極區(qū)成為全球氣候系統(tǒng)的冷源以及冷空氣的源地。再者,海冰的凍結(jié)和融化可以改變海洋上層溫鹽垂直結(jié)構(gòu),進(jìn)而潛在影響海洋表層的垂直混合和溫鹽環(huán)流。海冰的這些局地效應(yīng),可以通過(guò)大氣環(huán)流和海洋的熱力、動(dòng)力過(guò)程對(duì)遙遠(yuǎn)地區(qū)的氣候變化產(chǎn)生影響。影響我國(guó)夏季降水變化的因素主要分為大氣自身的動(dòng)力、熱力過(guò)程和外部強(qiáng)迫因子 (張慶云等,2007)。而夏季降水分布的不均衡,不僅是由大尺度季風(fēng)環(huán)流和熱帶海溫的異常變化所引發(fā) (Wu et al.,2009c;Zhou et al.,2009),來(lái)自高緯度地區(qū)大氣環(huán)流和外部強(qiáng)迫異常也會(huì)對(duì)東亞地區(qū)的降水變化產(chǎn)生影響 (武炳義等,1999;Wu et al.,2008)。

    北極海冰對(duì)北半球,尤其中國(guó)氣候變率影響的可能機(jī)理已得到廣泛關(guān)注。北極海冰異常對(duì)大氣環(huán)流和氣候的影響是通過(guò)激發(fā)全球大氣異常遙相關(guān)型來(lái)實(shí)現(xiàn)的 (謝倩和黃士松,1990;楊修群和黃士松,1993;黃士松等,1995)。Plumb and A lan(1985)認(rèn)為準(zhǔn)定常行星波的三維傳播是由地形效應(yīng)、非絕熱加熱及其與瞬變渦旋的相互作用三種機(jī)制來(lái)實(shí)現(xiàn),并沿一定的波導(dǎo)傳播。

    北極海冰變化對(duì)北半球氣候有顯著的影響。A lexander et al.(2004)利用CCM 3.6研究了冬季大氣在副極地海域局地響應(yīng)明顯,高度一般低于700 hPa[這與Wu et al.(2004)觀測(cè)分析結(jié)果是一致的],當(dāng)海冰范圍異常減少時(shí),有巨大的向上熱通量異常、近表層加熱、降水增加以及異常偏低的海平面氣壓;當(dāng)海冰異常偏多時(shí)則相反;對(duì)大尺度大氣環(huán)流的反饋?zhàn)饔迷谔窖蠛^(qū)主要表現(xiàn)為正反饋,而在大西洋地區(qū)主要是負(fù)反饋。Deser et al.(2004)研究了模式大氣對(duì)北大西洋海溫和海冰異常的直接和間接響應(yīng):間接響應(yīng)主導(dǎo)了半球尺度的位勢(shì)高度異常,類(lèi)似于準(zhǔn)正壓的NAO和北半球環(huán)狀模;而直接響應(yīng)僅局限于強(qiáng)迫源區(qū)附近,在垂直方向呈現(xiàn)斜壓結(jié)構(gòu)。Deser et al.(2007)研究了冬季大氣響應(yīng)的瞬態(tài)演變的過(guò)程:初始響應(yīng)在強(qiáng)迫區(qū)域附近表現(xiàn)為斜壓結(jié)構(gòu),由低層異常非絕熱加熱強(qiáng)迫所維持,最大振幅持續(xù)2~3周后漸變?yōu)闇?zhǔn)正壓結(jié)構(gòu),這種調(diào)整在第2~2.5個(gè)月時(shí)達(dá)到平衡,表現(xiàn)為類(lèi)似于北大西洋濤動(dòng)—北半球環(huán)狀模 (NAO或NAM)的特征。

    北極海冰對(duì)中國(guó)氣候的影響也表現(xiàn)在很多方面。王嘉和郭裕福 (2004)研究全球大氣環(huán)流模式A GCM R15L 9,海冰的季節(jié)和跨季節(jié)潛在可預(yù)報(bào)性結(jié)果的局地特征比較明顯,北半球某些區(qū)域海冰在面積發(fā)生較大異常的時(shí)候,可能對(duì)同期或 (及)后期環(huán)流有著比較重要的影響。Wu et al.(2004)研究表明冬季巴倫支海和喀拉海海冰變化通過(guò)影響大氣遙相關(guān)型和冷空氣活動(dòng)進(jìn)而影響東亞冬季風(fēng)。巴倫支海6月海冰偏多時(shí),造成局地500 hPa位勢(shì)高度場(chǎng)異常偏低,貝加爾湖地區(qū)和鄂霍次克海地區(qū)異常升高,這種歐亞波列結(jié)構(gòu)會(huì)造成中國(guó)東南部地區(qū)降水偏多,北部偏少 (Wang and Guo,2004)。Wu et al.(2009a,2009b)研究了北冰洋、格陵蘭海春季海冰密集度 (Sea Ice Concentration,簡(jiǎn)稱(chēng)SIC)減少 (增多),歐亞大陸邊緣海區(qū)海冰SIC增多 (減少),中國(guó)東北、黃河長(zhǎng)江之間的中部地區(qū)夏季降水會(huì)異常增多 (減少),南部降水異常減少 (增多),并提出了兩種可能的聯(lián)系機(jī)理:(1)春季北極海冰和歐亞大陸積雪一致變化,及其對(duì)后期夏季歐亞大陸緯向遙相關(guān)波列的影響可以解釋春季北極海冰與中國(guó)夏季降水異??臻g分布的關(guān)系;(2)春季海冰密集度的減少可能導(dǎo)致夏季北極偶極子的異常,對(duì)中國(guó)東北夏季降水有重要影響 (Wu et al.,2008)。

    因?yàn)榇髿鈱?duì)海冰具有優(yōu)勢(shì)影響作用,所以直接由海冰觀測(cè)分析揭示對(duì)大氣反饋?zhàn)饔煤芾щy。目前大多數(shù)研究關(guān)注于冬季海冰變化對(duì)同期大氣環(huán)流的影響,關(guān)注春季海冰的不多;數(shù)值模擬也多關(guān)注于北極、大西洋地區(qū)和太平洋的波列結(jié)構(gòu)、ENSO型,具體關(guān)注東亞地區(qū)氣候狀況的不多;大氣的初始響應(yīng)及其動(dòng)力機(jī)制的研究方面也存在不足,需要進(jìn)一步得到研究、驗(yàn)證和解釋。因此,本文運(yùn)用美國(guó)大氣環(huán)流模式 (Community A tmosphere Model,CAM 3.0)來(lái)模擬研究春季北極海冰變化在不同時(shí)間尺度上對(duì)東亞春、夏季氣候變化的可能影響,從歐亞波列遙相關(guān)和波活動(dòng)通量角度上揭示了其影響機(jī)理和可預(yù)報(bào)性。關(guān)注季節(jié)不同、研究地區(qū)不同、使用模式不同和試驗(yàn)方法不同,使得本文有別于其他的海冰研究。

    2 資料和方法

    2.1 資料

    自1978年以來(lái),被動(dòng)式微波輻射計(jì) (PMR)衛(wèi)星資料的使用,為整個(gè)北極地區(qū)提供了最綜合和全面的海冰 (SIC)信息。由此,形成了很多從PMR數(shù)據(jù)資料中獲取SIC數(shù)據(jù)的算法,比如美國(guó)宇航局 (NASA)的 TEAM(1991)和BOOTSTRAP(1997)算法等。之前,很多學(xué)者經(jīng)過(guò)對(duì)比發(fā)現(xiàn)BOOTSTRAP算法更為可靠一些 (Comiso et al.,1997;Singarayer and Bamber,2003)。本文采用的海冰資料為經(jīng)由BOO TSTRAP算法而獲取的SIC數(shù)據(jù),具體如下:

    月平均北極海冰密集度 1.0°(緯度)×1.0°(經(jīng)度)格點(diǎn)資料,取自英國(guó)大氣數(shù)據(jù)中心 (BADC,http:∥badc.nerc.ac.uk/data/hadisst/[2009-07-06])的 Had ISST ICE v.1.1資料集 (Rayner et al.,2000),時(shí)間從1960年1月至 2008年 12月共 49年。海冰密集度值為0~100,0表示單位網(wǎng)格點(diǎn)當(dāng)月月平均無(wú)冰覆蓋,100表示網(wǎng)格點(diǎn)范圍當(dāng)月月平均全部為冰覆蓋,1表示網(wǎng)格點(diǎn)范圍有1%為冰覆蓋,其它值依此類(lèi)推。

    位勢(shì)高度場(chǎng)資料:全球500 hPa位勢(shì)高度場(chǎng)2.5°×2.5°格點(diǎn)資料,取自美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心(NCEP,下載時(shí)間:2010年7月21日),時(shí)間長(zhǎng)度為1950年~2010年。

    降水資料:中國(guó)180站站點(diǎn)降水資料,取自國(guó)家氣象信息中心,時(shí)間長(zhǎng)度為1951~2010年。

    2.2 方法

    北極SIC和北半球氣候變化顯著相關(guān),為了進(jìn)一步研究北極海冰對(duì)東亞地區(qū)大氣環(huán)流影響的物理過(guò)程及影響機(jī)理,我們進(jìn)行數(shù)值模式試驗(yàn)來(lái)探討和研究問(wèn)題。本文采用的數(shù)值模式是NCAR第六代大氣環(huán)流模式CAM 3(Collins et al.,2006),它對(duì)東亞氣候和大型天氣系統(tǒng)的模擬能力已經(jīng)得到檢驗(yàn),如氣候態(tài)和季節(jié)循環(huán)的東亞季風(fēng)系統(tǒng):西藏高壓、東太平洋副熱帶高壓等 (Chen et al.,2010)。CAM 3歐拉譜模式選用 T85分辨率為1.4°×1.4°,在水平方向上,緯向和經(jīng)向?yàn)榫鶆蚍植嫉?56×128個(gè)高斯格點(diǎn),在垂直方向上采用σ-p混合坐標(biāo)系,共分為 26層,坐標(biāo)層頂在 2.917 hPa,其中從83.1425 hPa向上為純氣壓坐標(biāo),中間層為混合σ-p坐標(biāo),底層采用純?chǔ)易鴺?biāo)。CAM 3選用非耦合的DOM海洋模式,SST選用季節(jié)平均的氣候態(tài)資料,模式中指定北半球冰厚為2m,南半球冰厚為0.5 m。

    2.2.1 控制試驗(yàn)

    將CAM 3模式積分14年,其中北極SIC強(qiáng)迫取模式中的氣候態(tài)資料。將模擬試驗(yàn)中后10年大氣環(huán)流場(chǎng)等要素場(chǎng)和環(huán)流場(chǎng)的平均作為控制試驗(yàn)的結(jié)果。

    2.2.2 敏感性試驗(yàn)

    對(duì)觀測(cè)的1968~2005年50.5°N以北的春季北極SIC資料作EOF分析,第一模態(tài) (13%)反映了春季海冰的總體減少趨勢(shì),第二模態(tài) (10%)反映了春季海冰變率的主要特征。由于數(shù)值模擬試驗(yàn)是為了驗(yàn)證春季海冰與夏季氣候,尤其是降水的關(guān)系,首先,夏季降水 EOF的第一模態(tài)的解釋方差也為10%,反映了其變化的主要特征;其次,春季海冰和后期夏季降水SVD第一模態(tài)的海冰空間分布型與海冰EOF第二模態(tài) (EOF2)的空間分布型基本一致,體現(xiàn)了EOF2這種分布型與夏季降水相關(guān)強(qiáng)烈,最能體現(xiàn)其變化特征,因此,EOF2完全可以體現(xiàn)春季海冰變率的特征。取第二模態(tài)(EOF2)的時(shí)間系數(shù) (圖1b)對(duì)SIC原場(chǎng)進(jìn)行一元線(xiàn)性回歸,將回歸的空間異常場(chǎng) (圖1a)逐月的加到SIC春季氣候態(tài)上,作為敏感性試驗(yàn)的 SIC強(qiáng)迫。當(dāng)EOF2的時(shí)間系數(shù)為正時(shí),定義該試驗(yàn)為正位相試驗(yàn),相應(yīng)的強(qiáng)迫為正位相強(qiáng)迫;反之,時(shí)間系數(shù)為負(fù)時(shí),定義為負(fù)位相試驗(yàn),相應(yīng)的強(qiáng)迫為負(fù)位相強(qiáng)迫。

    由圖1a可以看出,在春季正位相時(shí),負(fù)的SIC異常區(qū)位于格陵蘭海、白令海東部、北極海盆以及巴芬灣部分地區(qū),正的SIC異常區(qū)位于北歐亞沿海,尤其是巴倫支海以及鄂霍次克海。負(fù)位相時(shí)SIC異常區(qū)的空間分布型則正好相反。

    圖1 春季北極SIC的 (a)正位相異??臻g分布型和 (b)EOF第二模態(tài)的時(shí)間序列Fig.1 (a)The spatial distribution of positive phase of sp ring A rctic SIC;(b)time series of the second EOFmode of sp ring A rctic SIC

    在正位相試驗(yàn)中,將模式第5~14年分別由3月1日積分到8月31日,初始場(chǎng)為每年的2月22日和2月28日,可得到20個(gè)試驗(yàn)個(gè)例;負(fù)位相試驗(yàn)亦然。因此,在相同的初始場(chǎng)下,分別可以得到20個(gè)正位相和20個(gè)負(fù)位相敏感性試驗(yàn)的結(jié)果。

    為了獲取大氣對(duì)指定春季海冰SIC強(qiáng)迫的季節(jié)響應(yīng),我們用正位相試驗(yàn)減去負(fù)位相試驗(yàn)作為敏感性試驗(yàn)的結(jié)果進(jìn)行分析。統(tǒng)計(jì)性檢驗(yàn)采用雙邊t檢驗(yàn)方法。

    3 試驗(yàn)結(jié)果分析

    3.1 模式大氣環(huán)流對(duì)北極海冰密集度變化的響應(yīng)

    季節(jié)尺度模式大氣環(huán)流響應(yīng)定義為,將季節(jié)平均的大氣環(huán)流響應(yīng)的正位相試驗(yàn)結(jié)果與負(fù)位相試驗(yàn)結(jié)果作差,可以得到20個(gè)個(gè)例的試驗(yàn)結(jié)果。目的在于既剔除了模式初始條件對(duì)模擬結(jié)果的影響,又使得大氣異常響應(yīng)加倍,異常形勢(shì)更加顯著。春季為3~5月,夏季為6~8月。以下為各個(gè)大氣環(huán)流場(chǎng)和要素場(chǎng)的響應(yīng)型。

    圖2 500 hPa位勢(shì)高度場(chǎng)異常響應(yīng)的EOF第一模態(tài) (單位:m):(a)春季 (16%);(b)夏季 (28%)Fig.2 (a)The spatial distribution of thefirst EOFmode of differences in sp ring 500-hPa height betw een positive and negative phases in sensitivity tests(units:m);(b)same as in(a),but for summer.The outer circle is 10°N

    3.1.1 位勢(shì)高度場(chǎng)響應(yīng)

    對(duì)20個(gè)個(gè)例的模式結(jié)果的500 hPa位勢(shì)高度場(chǎng) (Z500)作 EOF分析,第一模態(tài)的方差貢獻(xiàn)在春、夏季分別為16%、28%。結(jié)果發(fā)現(xiàn)500 hPa位勢(shì)高度場(chǎng)響應(yīng)的第一模態(tài),在春季表現(xiàn)為緯向上的波列結(jié)構(gòu),兩個(gè)正中心分別位于西歐和日本海,負(fù)中心位于西歐沿海和中西伯利亞地區(qū);在夏季同樣表現(xiàn)為波列結(jié)構(gòu),但振幅相對(duì)減弱,正中心位于烏拉爾山和鄂霍次克海,負(fù)中心位于東歐和中國(guó)東北部 (圖2)。從 EOF結(jié)果還可以看出,高緯度擾動(dòng)中心振幅在春季要大于夏季,這可能與海冰外強(qiáng)迫加在極地和高緯有關(guān),這時(shí)候大氣環(huán)流的活動(dòng)不只存在著大氣內(nèi)部調(diào)整,還會(huì)一直存在外強(qiáng)迫的直接熱力響應(yīng)。

    對(duì)Z500作EOF分析后,依據(jù)第一模態(tài)的標(biāo)準(zhǔn)化時(shí)間系數(shù)作合成分析,來(lái)分別選取正、負(fù)位相9個(gè)個(gè)例,并將二者均值相減,發(fā)現(xiàn)其結(jié)果與20個(gè)個(gè)例集合平均的結(jié)果大致趨勢(shì)相同,但是信號(hào)更強(qiáng)(如圖3)。以后試驗(yàn)也是如此。

    春季的500 hPa場(chǎng)與 EOF第一模態(tài)類(lèi)似表現(xiàn)為-+-的緯向波列結(jié)構(gòu),且由春季逐漸東移向夏季演變,位于歐亞大陸的正負(fù)中心均振幅增強(qiáng)范圍增大,位于西歐沿海的負(fù)中心到了夏季逐漸消退,負(fù)中心變?yōu)橛杀睒O延伸至歐洲北部,中國(guó)東北地區(qū)被高度負(fù)異常占據(jù),烏拉爾山和鄂霍次克海地區(qū)被正異常占據(jù),可看出歐亞波列由西向東逐漸傳播。將模式模擬結(jié)果與觀測(cè)結(jié)果作對(duì)比,發(fā)現(xiàn)二者在歐亞地區(qū)結(jié)果非常相似,波列型式也非常一致,因此可以驗(yàn)證數(shù)值模擬結(jié)果具有可靠性,春季北極海冰的自然變率對(duì)大氣環(huán)流的調(diào)整影響非常顯著。

    根據(jù)天氣學(xué)原理,結(jié)合表層熱通量場(chǎng) (圖4b)、溫度場(chǎng) (圖5b)分析,中國(guó)東北地區(qū)在春夏季都被一個(gè)冷性渦旋 (圖6a、b)占據(jù),這將非常有利于東北地區(qū)降水的發(fā)生。Tao and Chen(1987)、張慶云和陶詩(shī)言 (1998)研究表明,烏拉爾山和鄂霍次克海地區(qū)在夏季被雙阻結(jié)構(gòu)占據(jù),亞洲中高緯度及東亞?wèn)|部地區(qū)的距平場(chǎng)易分別形成+-+的距平波列,東亞地區(qū)距平波列的這種雙阻形勢(shì),及其相互作用常常造成東亞夏季特別是梅雨期降水偏多 (圖7b)。

    3.1.2 海平面氣壓響應(yīng)

    通常海平面氣壓場(chǎng)的變化是與冷熱源的分布緊密相連的,會(huì)與表層溫度產(chǎn)生相一致的變化。春季,北極大部分地區(qū)被負(fù)的海平面氣壓 (SLP)異常所占據(jù),歐亞大陸為正的SLP異常,與500 hPa位勢(shì)高度場(chǎng)相一致,在垂直方向上表現(xiàn)出正壓的特征;溫度場(chǎng)上,在北極地區(qū)表現(xiàn)為正的溫度異常,對(duì)應(yīng)該區(qū)偏少的海冰。夏季,北極地區(qū)的負(fù)異常范圍集中在北極海盆,振幅增強(qiáng);中高緯的正異常在鄂霍次克海附近,負(fù)異常在中國(guó)東部地區(qū),與垂直方向上的500 hPa位勢(shì)高度場(chǎng)同樣呈現(xiàn)了正壓結(jié)構(gòu)(圖8),由于夏季不存在海冰外強(qiáng)迫,因此極區(qū)的溫度異常很弱。

    圖3 模擬結(jié)果 (a、b)和觀測(cè)結(jié)果 (c、d)合成的500 hPa位勢(shì)高度場(chǎng) (等值線(xiàn))響應(yīng) (單位:m):(a、c)春季;(b、d)夏季。淺、深陰影:雙邊t檢驗(yàn)置信度超過(guò)90%、95%的區(qū)域 (下同)Fig.3 (a,b)Simulated and(c,d)observed geopotential height responses(units:m)at 500 hPa(contour):(a,c)sp ring;(b,d)summer.Light and dark shadings denote values significant at the 90%and 95%confidence levels based on a two-sided Student's t test(the same below)

    3.1.3 表層熱通量

    北極SIC的變化首先導(dǎo)致了極區(qū)表層熱通量的異常變化,并通過(guò)熱動(dòng)力學(xué)過(guò)程影響位于其上的大氣,最后通過(guò)波動(dòng)過(guò)程和大氣環(huán)流將能量異常傳播至北極以外的地區(qū)。在北極SIC偏少的區(qū)域,對(duì)應(yīng)于熱通量由海洋向大氣的凈向上輸送,SIC偏多區(qū)域反之。春季,在海冰的異常減少區(qū),如北冰洋中心區(qū)域、格陵蘭島東部有凈向上的熱通量,并且通過(guò)遙相關(guān)影響,在波列傳播過(guò)程中東太平洋和大西洋也產(chǎn)生了向上的熱通量;在海冰增多區(qū),如巴倫支海、巴芬灣和白令海有凈向下的熱通量,西太平洋、冰島地區(qū)也有凈向下的熱通量。夏季,在格陵蘭海、歐洲地區(qū)和中國(guó)東部地區(qū)有凈向上的熱通量,在大氣內(nèi)部動(dòng)力學(xué)作用下,可能導(dǎo)致中國(guó)東北地區(qū)降水的發(fā)生 (圖4)。通過(guò)研究還可以發(fā)現(xiàn),陸地上的熱通量異常很小,而海洋上熱通量變化卻很大,尤其是中低緯地區(qū),這可能與夏季短波輻射由南向北的遞減以及熱帶海洋的潛熱能釋放有關(guān)。凈熱通量異常是由感熱、潛熱和長(zhǎng)短波輻射通量共同作用造成的,短波、感熱和潛熱通量異常在春季量級(jí)相當(dāng),長(zhǎng)波通量異常量級(jí)較小約為前者的一半;夏季的熱通量均有增強(qiáng),這時(shí)候短波、潛熱通量占據(jù)主要地位,這與夏季太陽(yáng)輻射的增強(qiáng)和降水的增多有直接關(guān)系。

    圖4 合成的表層熱通量響應(yīng) (單位:W/m2):(a)春季;(b)夏季。等值線(xiàn)間隔為5(零線(xiàn)已省略);最外圈表示30°N;陰影同圖3Fig.4 Composite surface heat flux(units:W/m2):(a)Sp ring;(b)summer.The contour interval is 5 W/m2,and the zero contour isomitted;outer circle is 30°N

    圖5 同圖4,但為地表溫度場(chǎng) (單位:℃)。等值線(xiàn)間隔:0.5Fig.5 A s in Fig.4,but for surface air temperature(units:℃).Contour interval is 0.5

    3.1.4 地表溫度場(chǎng)響應(yīng)和垂直方向溫度響應(yīng)

    海冰異常區(qū)在垂直方向上存在熱通量異常,而表層溫度 (Ts)與表層熱通量在大部分地區(qū)的分布通常是一致的,其一致性表明大氣環(huán)流驅(qū)動(dòng)溫度平流可以對(duì)冷熱源以外遙遠(yuǎn)的地區(qū)產(chǎn)生一定程度的熱力響應(yīng) (Wu et al.,2009b)。

    在春季,表層溫度響應(yīng)在北極海盆、格陵蘭海以及巴芬灣為正的異常,對(duì)應(yīng)SIC偏少區(qū),甚至正異常還傳播至歐亞大陸大部分地區(qū),包括中國(guó)東部,并且與降水有相似的振幅。相反地,北歐亞沿海、拉布拉多海、鄂霍次克海表層溫度響應(yīng)為負(fù)異常,對(duì)應(yīng)SIC偏多區(qū)。需要特別指出,最大正值區(qū)位于格陵蘭海東部>1.5℃,最大負(fù)值區(qū)位于巴倫支海和鄂霍次克海<-1.0℃,分別對(duì)應(yīng)SIC異常的極值區(qū) (圖5a)。到了夏季,表層溫度異常向南傳播,極區(qū)異常不明顯,亞歐大陸和北美大陸響應(yīng)范圍異常增大,垂直方剖面上緯向?yàn)樵鰪?qiáng)的波列結(jié)構(gòu),與500 hPa位勢(shì)高度場(chǎng)變化完全一致。在110°E~130°E東北地區(qū)內(nèi)可以看到冷性渦旋結(jié)構(gòu)的存在,高層負(fù)的溫度異常大于低層,在低層加熱作用下大氣重力能減弱,垂直運(yùn)動(dòng)發(fā)展,易于形成對(duì)流性降水 (圖7b)。

    圖6 同圖4,但為850 hPa的矢量風(fēng)場(chǎng) (單位:m/s)Fig.6 As in Fig.4,but for wind vecto rat 850 hPa

    3.1.5 850 hPa風(fēng)場(chǎng)響應(yīng)

    850 hPa風(fēng)場(chǎng) (V850)在春季表現(xiàn)為中國(guó)東部顯著的北風(fēng),在中國(guó)東南部形成反氣旋性環(huán)流 (圖6)。東部地區(qū)處于弱氣旋性渦旋的后方,低層輻合不明顯,雖然偏北氣流將北部的干冷空氣帶至南方,但無(wú)暖濕空氣的匯合,使得水汽條件不充足,不利于降水的發(fā)生,因此東部地區(qū)以降水負(fù)異常為主 (圖7a)。在夏季,印度季風(fēng)、南海季風(fēng)和副熱帶季風(fēng)全面爆發(fā),東亞副熱帶季風(fēng)異常水汽輸送也隨之增強(qiáng),渤海灣地區(qū)氣旋性渦旋結(jié)構(gòu)明顯,中國(guó)南部至西北太平洋地區(qū)為顯著的西南氣流,北方地區(qū)為偏北氣流,渤海灣地區(qū)上升運(yùn)動(dòng)強(qiáng)烈,這與500 hPa位勢(shì)高度場(chǎng)是動(dòng)力一致的,動(dòng)力和水汽條件滿(mǎn)足后非常有利于東部沿海地區(qū)降水的發(fā)生。結(jié)合高度場(chǎng)分析蒙古地區(qū)因水汽條件不明顯,華南地區(qū)因動(dòng)力機(jī)制不足都降水偏少。

    3.1.6 降水響應(yīng)

    大氣環(huán)流場(chǎng)的異常分布形勢(shì)互相配置,導(dǎo)致了東亞地區(qū)降水的異常分布。在降水場(chǎng)上,春季,中國(guó)東部地區(qū)整體降水偏少為-80 mm,顯著中心分別位于東北、長(zhǎng)江黃河的下游區(qū)域和華南地區(qū),這是由于東部無(wú)水汽輸送及動(dòng)力機(jī)制很弱。夏季,與春季差異較大,中國(guó)東北部地區(qū)以及長(zhǎng)江和黃河之間的中部地區(qū)被降水正異常區(qū)所覆蓋,尤其是東北地區(qū)降水很顯著,最大異常可以達(dá)到180 mm以上,而蒙古、華南地區(qū)降水偏少 (圖7b),這與渤海地區(qū)氣旋結(jié)構(gòu)、海洋水汽輸送以及副熱帶季風(fēng)的加強(qiáng)有很重要的關(guān)系。與觀測(cè)合成的結(jié)果相對(duì)比,可以發(fā)現(xiàn)二者在東部降水的空間分布型上很相似,即東北和長(zhǎng)江、黃河之間地區(qū)降水顯著偏多,華南地區(qū)降水顯著偏少,其中模擬的東北降水較之觀測(cè)更偏東,這與模擬的Z500東北低渦位置偏東是一致的。

    可見(jiàn),在春季極區(qū)SIC存在負(fù)異常,通過(guò)熱力和動(dòng)力機(jī)制,表層熱通量發(fā)生改變,不僅會(huì)對(duì)局地的SLP、高度場(chǎng)、溫度場(chǎng)造成影響,并且可以通過(guò)大氣環(huán)流和遙相關(guān)過(guò)程,尤其是在500 hPa高度場(chǎng)上的歐亞波列結(jié)構(gòu),對(duì)高層大氣和遙遠(yuǎn)地區(qū)造成影響,將能量由極區(qū)傳送至東亞地區(qū),從而影響東亞氣候。

    圖8 同圖4,但為合成的SLP場(chǎng)響應(yīng) (單位:hPa):(a)春季;(b)夏季Fig.8 As in Fig.4,but for sea level pressure(hPa):(a)Sp ring;(b)summer

    3.2 模式大氣環(huán)流對(duì)北極SIC響應(yīng)的演變過(guò)程

    以上研究是從季節(jié)平均的角度來(lái)討論的海冰對(duì)模式大氣的影響,現(xiàn)在進(jìn)一步分析模式大氣對(duì)于強(qiáng)迫的響應(yīng)是如何產(chǎn)生并逐步演變的,演變過(guò)程從時(shí)間尺度上可分為逐日和逐周響應(yīng)。試驗(yàn)方法設(shè)計(jì)上稍稍不同于季節(jié)平均響應(yīng),每個(gè)個(gè)例的模式初始場(chǎng)設(shè)定為3月1日,并積分至8月31日,模擬結(jié)果為9個(gè)個(gè)例的合成分析。

    大氣響應(yīng)可以分為直接熱力響應(yīng)和間接熱力響應(yīng)兩部分,分別對(duì)應(yīng)于海冰外強(qiáng)迫和大氣內(nèi)部變率的間接強(qiáng)迫。1000 hPa和300 hPa位勢(shì)高度場(chǎng)的初始響應(yīng)是海冰直接熱力強(qiáng)迫的結(jié)果,如圖9所示為響應(yīng)前12天的逐日分布:時(shí)間尺度為3月2日~3月13日。位勢(shì)高度場(chǎng)初始響應(yīng)的高低層配置顯示了很強(qiáng)的反位相結(jié)構(gòu),在垂直方向上為斜壓狀態(tài):第一周,高低層異常響應(yīng)之間存在著一定的空間負(fù)相關(guān),第二天的相關(guān)系數(shù)達(dá)到了-0.22,異常響應(yīng)很弱,都維持在SIC強(qiáng)迫附近地區(qū),格陵蘭海的北部在低層開(kāi)始為異常低壓,后來(lái)逐漸被異常高壓所取代,而在高層始終維持著異常高壓;從3月9日開(kāi)始高度異常區(qū)范圍增大、南擴(kuò)加強(qiáng),振幅逐漸增強(qiáng)到40 m,在緯向上形成波列結(jié)構(gòu)緩慢向東傳播,圖10給出響應(yīng)的均方根振幅也顯示了由0異常非線(xiàn)性增大的過(guò)程,但反位相結(jié)構(gòu)開(kāi)始削弱,高低層空間相關(guān)系數(shù)由負(fù)到正;到3月13日,響應(yīng)幅度達(dá)到了半球尺度,除了喀拉海和白令海還存在較強(qiáng)的斜壓性以外,其他地區(qū)高低層配置都開(kāi)始緩慢向準(zhǔn)正壓型轉(zhuǎn)變,空間相關(guān)逐漸增強(qiáng),且呈現(xiàn)了與高度場(chǎng)響應(yīng)相同的變化趨勢(shì)??梢?jiàn),局地的初始斜壓響應(yīng)逐步向半球尺度準(zhǔn)正壓響應(yīng)轉(zhuǎn)變,這個(gè)過(guò)程發(fā)展較快,與北極SIC強(qiáng)迫范圍較小、強(qiáng)度較弱、強(qiáng)迫時(shí)間較短以及斜壓響應(yīng)的區(qū)域局限性有關(guān),轉(zhuǎn)變機(jī)制在下文會(huì)進(jìn)行探討。通過(guò)斜壓轉(zhuǎn)化,高緯地區(qū)存在能量由低層向高層、由高緯向低緯逐漸傳播,使得中高緯、高低層能量配置逐漸達(dá)到平衡。

    圖9 3月2日~3月13日1000 hPa和300 hPa逐日的位勢(shì)高度場(chǎng)響應(yīng) (單位:m)Fig.9 1000-and 300-hPa geopotential height responses from 2 Mar to 13 Mar in composite analysis

    圖10 Z1000和Z300逐日分布的空間相關(guān)場(chǎng) (點(diǎn)虛線(xiàn))及各自異常響應(yīng)的振幅 (虛、實(shí)線(xiàn),單位:m)Fig.10 Spatial correlation coefficients betw een daily geopotential height responses at 1000 and 300 hPa and their anomalous response amplitudes

    為了探討大氣對(duì)SIC強(qiáng)迫初始斜壓響應(yīng)的持續(xù)性以及波列結(jié)構(gòu)的形成過(guò)程,從3月起采用逐周時(shí)間尺度進(jìn)行分析。位勢(shì)高度場(chǎng)響應(yīng)隨時(shí)間的變化形勢(shì)如圖10(圖中振幅是40°N以北地區(qū)高度場(chǎng)響應(yīng)的均方根作區(qū)域平均的結(jié)果),分別給出了40°N以北1000 hPa異常位勢(shì)高度場(chǎng) (Z1000)和300 hPa異常位勢(shì)高度場(chǎng) (Z300)的逐日空間相關(guān)系數(shù)場(chǎng)和區(qū)域平均的振幅場(chǎng),在計(jì)算之前對(duì)異常場(chǎng)作加權(quán)和3點(diǎn)平滑。

    結(jié)合空間場(chǎng)圖11(見(jiàn)文后彩圖),從模式試驗(yàn)結(jié)果可以得到,局地的斜壓結(jié)構(gòu)具有一定的持續(xù)性:在第1~2周時(shí),低層的正負(fù)強(qiáng)迫附近分別形成偶極子結(jié)構(gòu)的正負(fù)異常中心,中高層正好與之是反位相,這種斜壓性在冰島地區(qū)、格陵蘭海、白令海、巴倫支海以及喀拉海可以持續(xù)到第3周,之后減弱,而到第6周大氣在垂直方向上完全轉(zhuǎn)為正壓狀態(tài),高度場(chǎng)振幅也達(dá)到最大的階段。相對(duì)比相關(guān)系數(shù)圖 (圖10),40°N以北整個(gè)中高緯地區(qū)的斜壓性較弱,在第一周就由最大的負(fù)相關(guān)-0.22減少至0并逐漸轉(zhuǎn)正,第5周之后維持很高的正相關(guān),并呈準(zhǔn)兩周和月周期變化。高低層空間場(chǎng)的振幅變化與相關(guān)場(chǎng)一致,大氣異常響應(yīng)由開(kāi)始的0異常到第4~6周達(dá)到高峰,之后緩慢下降達(dá)到平衡階段,并呈現(xiàn)出以月為周期的振蕩型。還可以看出Z300異常響應(yīng)振幅大約是Z1000異常響應(yīng)的1.5~2倍。

    歐亞遙相關(guān)型的形成和維持是能量得以從極地東傳到東亞地區(qū)的原因,其機(jī)制將于下一節(jié)中介紹。隨著初始斜壓場(chǎng)向正壓的演變,歐亞大陸中高緯整層大氣都有波列結(jié)構(gòu)的形成,并緩慢向東傳播。由于整層的正壓性,從500 hPa異常高度場(chǎng)來(lái)研究整層波列結(jié)構(gòu)的演變,其中以中心位于格陵蘭海、西北歐、中西伯利亞的負(fù)正負(fù)歐亞波列,逐漸向中心位于西北歐、中西伯利亞、東亞的歐亞波列的東移演變?yōu)橹?(圖11),通過(guò)這種遙相關(guān)的方式大氣環(huán)流將北極、歐洲與中國(guó)地區(qū)的氣候聯(lián)系起來(lái),從而造成對(duì)中國(guó)氣候的影響。因此,通過(guò)研究大氣的瞬態(tài)演變過(guò)程,既可以看到垂直方向上高低層大氣由斜壓向正壓結(jié)構(gòu)的轉(zhuǎn)化,又可以看到在水平方向上大氣波列結(jié)構(gòu)響應(yīng)由北極至東亞大陸的逐步演變,揭示了北極海冰是如何通過(guò)大氣環(huán)流來(lái)對(duì)中國(guó)氣候造成顯著的影響。

    3.3 動(dòng)力機(jī)制研究

    由大氣對(duì)北極SIC變化的瞬態(tài)響應(yīng)可以清楚地看到北半球大氣響應(yīng)從斜壓到正壓的緩慢演變過(guò)程,但是響應(yīng)的動(dòng)力機(jī)制卻沒(méi)有給出,下面嘗試從波活動(dòng)通量的角度來(lái)解釋海冰異常變化所產(chǎn)生能量在垂直和水平方向產(chǎn)生和傳播的過(guò)程,用大氣內(nèi)部的動(dòng)力學(xué)過(guò)程來(lái)解釋能量的維持和轉(zhuǎn)化。由于北極海冰的存在顯著地改變了海洋和大氣之間的熱量、動(dòng)量和物質(zhì)交換,海冰的高反照率也大大減少了極地對(duì)太陽(yáng)輻射的吸收,所以當(dāng)存在SIC的異常分布時(shí),首先改變的是表層熱通量的分布,然后位勢(shì)高度場(chǎng)發(fā)生相應(yīng)的改變,由于地形、非絕熱加熱及其與渦旋的相互作用是產(chǎn)生波活動(dòng)通量的機(jī)制,這種熱力和高度異常可以激發(fā)出波活動(dòng)通量向高層和水平方向傳播,并在位勢(shì)高度場(chǎng)上產(chǎn)生一致的準(zhǔn)靜止羅斯貝波結(jié)構(gòu)向東傳播影響東亞地區(qū) (Honda et al.,2009)。

    圖12 (a-c)沿10°E的位勢(shì)高度場(chǎng)異常 (等值線(xiàn),單位:m)和波活動(dòng)通量異常 (箭頭,垂直方向×100,單位:m2/s2)緯度—高度剖面圖;(d-f)1000 hPa熱通量異常 (陰影)和垂直波活動(dòng)通量異常 (等值線(xiàn),間隔0.1,零等值線(xiàn)已省略;單位:m2/s2);(g-i)300 hPa位勢(shì)高度場(chǎng)異常 (陰影,單位:m)和水平波活動(dòng)通量異常 (箭頭,單位:m2/s2)。(a、d、g)3月2~5日平均;(b、e、h)3月6~9日平均;(c、f、i)3月10~13日平均Fig.12 (a-c)Latitude-height cross section of geopotential height responses(contour,units:m)and wave activity flux responses(arrow with vertical componentmultip lied by 100,units:m2/s2)along 10°E;(d-f)heat flux responses(shaded)and vertical wave activity flux responses(contour with interval 0.1 m2/s2)at 1000 hPa;(g-i)geopotential height responses(shaded)and horizontal wave activity flux responses(arrow)at 300 hPa.(a,d,g)Averaged over 2-5 Mar;(b,e,h)averaged over 6-9 Mar;(c,f,i)averaged over 10-13 Mar

    極地大氣的初始響應(yīng)在垂直方向上具有斜壓結(jié)構(gòu),這是海冰直接熱力強(qiáng)迫的結(jié)果,由于75°N~90°N地區(qū)SIC異常偏少,所以在響應(yīng)的前4天,位勢(shì)高度場(chǎng)上表現(xiàn)為大氣表層負(fù)異常,925 hPa以上正異常,這種初始斜壓是由 Ekman抽吸機(jī)制所引起的 (Wu et al.,2004)。表層熱通量在格陵蘭附近表現(xiàn)為正異常 (圖未給出),由低層激發(fā)出向上的垂直波活動(dòng)通量;而在巴倫支海和拉布拉多海附近熱通量表現(xiàn)為負(fù)異常,垂直方向波活動(dòng)通量?jī)粝蛳聻樨?fù)值,可見(jiàn)最大熱通量異常和最大垂直波活動(dòng)通量在空間上是一致存在的 (圖12d);水平方向的波活動(dòng)通量局限于格陵蘭島的高度場(chǎng)異常梯度附近(圖12g)。在響應(yīng)的第5~8天,在垂直方向上,位勢(shì)高度場(chǎng)正異常振幅增強(qiáng),低層有凈向下的熱通量和波活動(dòng)通量,這是由于70°N處海冰異常較弱,在大氣環(huán)流調(diào)整下,高度場(chǎng)在低層由負(fù)異常轉(zhuǎn)變?yōu)檎惓?溫度異常也相應(yīng)由正轉(zhuǎn)負(fù),低層空氣不再是熱強(qiáng)迫源,因此不能產(chǎn)生向上的波活動(dòng)通量。方向向南向下的波活動(dòng)通量很強(qiáng),與位勢(shì)高度梯度方向一致 (圖12b);水平方向波活動(dòng)通量表現(xiàn)為在負(fù)位勢(shì)梯度區(qū)由大西洋向東傳播的波活動(dòng)通量 (圖12h)。在響應(yīng)的第9~12天,位勢(shì)高度場(chǎng)、熱通量和波活動(dòng)通量均增強(qiáng),垂直向上傳播的波活動(dòng)通量在70°N附近上升,在高層加強(qiáng)南傳,并在中緯度地區(qū)下沉 (圖12c),45°N和55°N附近下沉的形成是由于其處于中層高壓后部,北風(fēng)異常很強(qiáng),形成向下的波活動(dòng)通量,這與 Plumb and A lan(1985)在大西洋波活動(dòng)通量向東南方向傳播的結(jié)論一致;水平方向波活動(dòng)通量由歐洲向東亞地區(qū)東傳增強(qiáng),通量最大值位于歐洲的高度正異常區(qū),這是與阿爾卑斯山的地形強(qiáng)迫作用緊密相連的 (圖12i)。由此可見(jiàn),我們可認(rèn)為低層異常非絕熱加熱是波活動(dòng)通量上傳和東傳的強(qiáng)迫機(jī)制,這種非絕熱加熱是由異常SIC直接熱力強(qiáng)迫形成的,還可以推斷出熱通量就是大氣響應(yīng)的冷熱源,波活動(dòng)通量的上傳和下傳正好位于這對(duì)冷熱中心上方。因此,歐亞波列的遙相關(guān)是北極海冰影響東亞地區(qū)的主要因素,而由直接熱力強(qiáng)迫產(chǎn)生的波活動(dòng)通量正是遙相關(guān)型形成的直接動(dòng)力機(jī)制,由此能量可以通過(guò)三維準(zhǔn)定常行星波由極地向北半球和東亞地區(qū)傳播。

    因此,大氣遙相關(guān)型的產(chǎn)生可以解釋為大氣對(duì)海冰直接外強(qiáng)迫的響應(yīng),通過(guò)準(zhǔn)定常羅斯貝波向北半球傳播,但是其維持機(jī)制卻與大氣基本態(tài)結(jié)構(gòu)也有很重要的關(guān)系,因此也可以稱(chēng)作間接強(qiáng)迫機(jī)制,它通過(guò)大氣內(nèi)部的動(dòng)力學(xué)過(guò)程吸取能量并激發(fā)新的強(qiáng)迫產(chǎn)生,再激發(fā)新的大氣異常產(chǎn)生。從春季5月平均來(lái)看,大氣環(huán)流場(chǎng)已經(jīng)很穩(wěn)定,達(dá)到了平衡狀態(tài),位勢(shì)高度場(chǎng)在垂直方向演變?yōu)闇?zhǔn)正壓結(jié)構(gòu),這時(shí)候可以看到明顯的波活動(dòng)通量由極區(qū)向東傳播,在西部渦旋中心的東南側(cè)存在波活動(dòng)通量的最大輻散,這是與阿爾卑斯山地形也有一定關(guān)系的,而在東部渦旋的西側(cè)波活動(dòng)通量有最大輻合,其輻散和輻合過(guò)程都是由渦度源匯決定的 [圖13(見(jiàn)文后彩圖)]。波活動(dòng)通量的傳播過(guò)程正是正壓結(jié)構(gòu)的羅斯貝波頻散能量的過(guò)程,其中間接強(qiáng)迫的過(guò)程是非常復(fù)雜的,當(dāng)其傳至東亞季風(fēng)區(qū)時(shí)可以激發(fā)斜壓結(jié)構(gòu)的產(chǎn)生,釋放對(duì)流性潛熱;也可以從西風(fēng)帶等大氣環(huán)流基本態(tài)通過(guò)正壓能量轉(zhuǎn)換來(lái)吸取能量,維持大氣擾動(dòng)和遙相關(guān)型的存在 (楊修群和黃士松,1993)。另外,外強(qiáng)迫除了海冰的非絕熱加熱外,地形效應(yīng)和渦旋結(jié)構(gòu)也可以對(duì)波活動(dòng)通量的形成和傳播造成很大影響。大氣環(huán)流遙相關(guān)型是極冰外強(qiáng)迫作用所激發(fā)的波活動(dòng)通量與大氣內(nèi)部動(dòng)力學(xué)過(guò)程相互作用的結(jié)果,極冰的變化由此可以強(qiáng)烈地影響東亞地區(qū)的氣候。

    4 結(jié)論和討論

    通過(guò)對(duì)北極SIC資料的統(tǒng)計(jì)分析,發(fā)現(xiàn)春季北極海冰存在著一定的自然變率,這體現(xiàn)在 EOF的第二模態(tài),這種變率將對(duì)氣候變化產(chǎn)生很重要的影響。我們?cè)O(shè)計(jì)了一系列數(shù)值模擬試驗(yàn)來(lái)研究CAM 3中大氣對(duì)北極SIC強(qiáng)迫的響應(yīng)。將SIC第二模態(tài) (EOF2)的時(shí)間系數(shù)回歸到原場(chǎng),并把這個(gè)空間異常場(chǎng)逐月加到SIC春季氣候態(tài)上,作為敏感性試驗(yàn)的外強(qiáng)迫來(lái)研究大氣響應(yīng)過(guò)程。

    試驗(yàn)結(jié)果:在季節(jié)尺度上,大氣環(huán)流由春向夏演變,春季的SIC與 Z500、SLP、Ts、V850異常場(chǎng),通過(guò)熱動(dòng)力過(guò)程相互配置起來(lái),對(duì)春季中國(guó)東部負(fù)降水異常、夏季東北部正降水異常有很大貢獻(xiàn)。在逐日和逐周時(shí)間尺度上,體現(xiàn)了大氣響應(yīng)的瞬態(tài)演變過(guò)程,高低層位勢(shì)高度場(chǎng)的初始響應(yīng)在垂直方向上表現(xiàn)出很明顯的反位相結(jié)構(gòu),并局限于強(qiáng)迫區(qū)附近,隨后異常響應(yīng)范圍增大、振幅增強(qiáng),持續(xù)到第3周后逐漸轉(zhuǎn)為準(zhǔn)正壓結(jié)構(gòu),到了第6周完全達(dá)到平衡階段,振幅最大并呈周期性振蕩。

    用波活動(dòng)通量和大氣內(nèi)部動(dòng)力學(xué)過(guò)程來(lái)解釋遙相關(guān)的傳播和維持,分別作為直接強(qiáng)迫機(jī)制和間接強(qiáng)迫機(jī)制。直接熱力強(qiáng)迫的結(jié)果改變的是表層熱通量,并使位勢(shì)高度場(chǎng)產(chǎn)生異常,二者共同作用激發(fā)準(zhǔn)靜止羅斯貝波在三維方向的傳播,波活動(dòng)通量在低層初始斜壓區(qū)上傳,在高層產(chǎn)生一致的波列結(jié)構(gòu)并緩慢向東向南傳,最后在中低緯度下沉。間接強(qiáng)迫是大氣內(nèi)部動(dòng)力學(xué)過(guò)程的通過(guò)正壓轉(zhuǎn)換和從大氣基本態(tài)吸取能量等方式維持波列的傳播。通過(guò)兩種機(jī)制的綜合效應(yīng),歐亞遙相關(guān)型將能量傳遞到東亞地區(qū),從而影響當(dāng)?shù)氐臍夂颉?/p>

    運(yùn)用不同的大氣環(huán)流模式、不同的外強(qiáng)迫因子會(huì)得到不同的大氣響應(yīng)結(jié)果。因此,本文結(jié)果不僅非常具有模式依賴(lài)性,而且也有強(qiáng)迫因子的片面性,在以后的研究中可以嘗試用多種模式、多種外強(qiáng)迫因子來(lái)做模擬試驗(yàn),以得到更加科學(xué)的結(jié)論。除了SIC還存在多種因素的影響,導(dǎo)致了模式結(jié)果和觀測(cè)結(jié)果的不一致性,本文的機(jī)理研究只局限于波列遙相關(guān)、波活動(dòng)通量的熱動(dòng)力作用,大氣內(nèi)部動(dòng)力學(xué)機(jī)制的研究沒(méi)有得到完善,其他的機(jī)理研究也還需要進(jìn)一步的探索。

    References)

    Alexander M A,Bhatt U S,Walsh J E,et al.2004.The atmospheric response to realistic A rctic sea ice anomalies in an AGCM during winter[J].J.Climate,17:890-905.

    Chen H M,Zhou T J,Neale R B,et al.2010.Perfo rmance of the New NCAR CAM 3.5 in East A sian summer monsoon simulations:Sensitivity to modifications of the convection scheme[J].J.Climate,23:3657-3675.

    Collins WD,Bitz CM,Blackmon M L,et al.2006.The Community Climate System Model version 3(CCSM 3)[J].J.Climate,19:2122-2143.

    Comiso J C,Cavalieri D J,Parkinson C L,et al.1997.Passivemicrowave algorithms for sea ice concentration:A comparison of two techniques[J].Remote Sens.Environ.,60:357-384.

    Deser C,Magnusdottir G,Saravanan R.2004.The effectsof North A tlantic SST and sea ice anomalies on the winter circulation in CCM 3.Part II:Direct and indirect components of the response[J].J.Climate,17:877-889.

    Deser C,Tomas RA,Peng Shiling.2007.The transient atmospheric circulation response to North A tlantic SST and sea ice anomalies[J].J.Climate,20:4751-4767.

    Honda M,Inoue J,Yamane S.2009.Influence oflow A rctic sea-ice minima on anomalously cold Eurasian winters[J].Geophys.Res.Lett.,36.

    黃士松,楊修群,蔣全榮.1995.極地海冰變化對(duì)氣候的影響 [J].氣象科學(xué),15(4):46-56. Huang S S,Yang X Q,Jiang Q R.1995.The effectsof the polar sea ice on climate[J].Scientia Meteo rologica Sinica(in Chinese),15(4):46-56.

    彭公炳,李倩,錢(qián)步東.1992.氣候與冰雪覆蓋[M].北京:氣象出版社,66-84.Peng GB,Li Q,Qian B D.1992.Climate and Ice-Snow Cover(in Chinese)[M].Beijing:China Meteorological Press,66-84.

    Plumb R A.1985.On the three-dimensional propagation of stationary waves[J].J.A tmos.Sci.,42:217-229.

    Rayner N A,Parker D E,Frich P.2000.SST and sea-ice fields for ERA 40[C].Proc.Second WCRP Int.Conf.on Reanalysis,18-21.

    ScreenJ A,Simmonds I.2010.The central role of diminishing sea ice in recent A rctic temperature amplification[J].Nature,464:1334-1337,doi:10.1038/nature09051.

    Singarayer J S,Bamber J L.2003.EOF analysis of three recordsof sea-ice concentration spanning the last 30 years[J].Geophys.Res.Lett.,30:1251-1254.

    Tao S Y,Chen L X.1987.A review of recent research on the East Asia summer monsoon in China[M]//Chang C P,Krishnamurti T N.Monsoon Meteo rology.U K:Oxfo rd University Press,60-92.

    王嘉,郭裕福.2004.北半球海冰強(qiáng)迫作用下大氣可預(yù)報(bào)性研究[J].應(yīng)用氣象學(xué)報(bào),15(3):291-305. Wang J,Guo Y F.2004.Assessing potential predictability associated with bo real sea ice with an AGCM[J].Quarterly Journal of App lied Meteo rology(in Chinese),15(3):291-305.

    Wang J,Guo Y F.2004.Possible impacts of Barents Sea Ice on the Eurasian atmospheric circulation and the rainfall of East China in the beginning of summer[J].Adv.A tmos.Sci,21:662-674.

    Walsh J E.1983.The role of Sea Ice in climate variability:Theo ries and evidence[J].A tmosphere-Ocean,21:229-242.

    武炳義,黃榮輝,高登義.1999.冬季北極喀拉海、巴倫支海海冰面積變化對(duì)東亞冬季風(fēng)的影響 [J].大氣科學(xué),23(3):267-275.Wu B Y,Huang R H,Gao D Y.1999.The impact of variation of sea-ice extent in the Kara Sea and the Barents Seas in winter on the winter monsoon over East A sia[J].Chinese J.A tmos.Sci.(in Chinese),23(3):267-275.

    Wu B Y,Wang J,Walsh J.2004.Possible feedback of winter sea ice in the Greenland and the Barents Sea on the local atmosphere[J].Mon.Wea.Rev.,132:1868-1876.

    Wu B Y,Zhang R H,Rosanne D'A rrigo.2008.A rctic Dipole A-nomaly and Summer Rainfall in Northeast China[J].Chinese Science Bulletin,53:2222-2229.

    Wu B Y,Zhang R H,Wang B.2009a.On the association between sp ring A rctic sea ice concentration and Chinese summer rainfall[J]. Geophys. Res. Lett.,36,L 09501,doi: 10.1029/2009GL 037299.

    Wu B Y,Zhang R H,Wang B.2009b.On the association between sp ring A rctic sea ice concentration and Chinese summer rainfall:A further study[J].Adv.A tmos.Sci.,26:666-678.

    Wu B Y,Zhou T J,Li T.2009c.Seasonally evolving dominant interannual variabilitymodesof East Asian climate[J].J.Climate,22:2992-3005.

    謝倩,黃士松.1990.冬季赤道中東太平洋海溫和北極海冰異常對(duì)大氣環(huán)流影響的研究 [J].氣象科學(xué),10:325-338. Xie Q,Huang S S.1990.A study of the effects of anomalies of the cen-tral-eastern equatorial Pacific sea surtace temperature and the A rctic sea ice cover on the atmospheric general circulation during the northern winter[J].Scientia Meteorologica Sinica(in Chinese),10:325-338.

    楊修群,黃士松.1993.外強(qiáng)迫引起的夏季大氣環(huán)流異常及其機(jī)制討論[J].大氣科學(xué),17(6):697-702. Yang X Q,Huang S S.1993.On the mechanism of summertime atmospheric anomalies induced by external forcing[J].Chinese J.A tmos.Sci.(Scientia A tmospherica Sinica)(in Chinese),17(6):697-702.

    張慶云,陶詩(shī)言.1998.亞洲中高緯度環(huán)流對(duì)東亞夏季降水的影響[J].氣象學(xué)報(bào),6(2). Zhang Q Y,Tao S Y.1998.Influence of Asian mid high latitude circulation on East Asian summer rainfall[J].Acta Meteor.Sin.(in Chinese),6(2).

    張慶云,呂俊梅,楊蓮梅,等.2007.夏季中國(guó)降水型的年代際變化與大氣內(nèi)部動(dòng)力過(guò)程及外強(qiáng)迫因子關(guān)系[J].大氣科學(xué),24(6):1290-1300.Zhang Q Y,LüJ M,Yang L M.2007.The interdecadal variation of precipitation pattern over China during Summer and its relationship with the atmospheric internal dynamic processes and extra-forcing factors[J].Chinese J.A tmos.Sci.(in Chinese),24(6):1290-1300.

    Zhou T J,Gong D Y,Li J.2009.Detecting and understanding the multi-decadal variability of the East Asian summer monsoon:Recent progressand stateof affairs[J].Meteorologische Zeitschrift,18(4):455-467.

    圖11 合成的500 hPa逐周的位勢(shì)高度場(chǎng)響應(yīng) (3月2日到8月30日)Fig.11 Weekly averages of 500-hPa geopotential height responses during the integration(2 Mar-30 Aug)

    圖13 (a)5月平均沿60°N的位勢(shì)高度場(chǎng)異常 (等值線(xiàn),單位:m)和波活動(dòng)通量異常 (箭頭,垂直分量×100,單位:m2/s2)經(jīng)度—高度剖面圖;(b)同圖12d,但為5月平均;(c)同圖12g,但為5月平均Fig.13 Longitude-height cross section of geopotential height responses(contour,units:m)and wave activity flux responses(arrow with vertical component multiplied by 100,units:m2/s2)along 60°N;(b)same as Fig.12d,but for May average;(c)same as Fig.12d,but for May average

    The Northern Hem isphere Atmospheric Response to Spring Arctic Sea Ice Anomalies in CAM 3.0 Model

    ZHANG Ruonan and WU Bingyi

    ChineseAcademyofMeteorologicalSciences,Beijing100081

    The influence of A rctic Sea Ice Concentration(SIC)anomalieson the atmospheric general circulation during sp ring and summer is investigated with version 3.0 of the NCAR Community A tmosphere Model(CAM 3.0).Regress the second EOFmodeof ICE to observed ICE,and then added to climatic ICE as the external force of test.The results show that on the intraseasonal time scale,the atmospheric circulation anomalies evolved from sp ring to summer,and the sp ring A rctic SIC anomalies were thermally and dynamically consistent with atmospheric circulation,surface temperature,and rainfall anomalies in sp ring and summer,especially contributing to less sp ring rainfall in eastern China and more summer rainfall in northeastern and central China.While on short time scales,the initial adjustment of the atmospheric circulation is characterized by relationship of outof phase in geopotential height anom-alies in the lower and upper troposphere locally,and after two weeks,the out of phase turned progressively to mo re barotropic with the responses propagating to remote areas,at last,the quasi-equilibrium stage of adjustment is reached in the sixth week.The remote responses are regarded as a stationary Rossby wave generated thermally and dynamically through an anomalous turbulent heat fluxes and atmospheric circulation internal varieties.A t first,the surface heat fluxes were changed by anomalies in the A rctic SIC,and then a large scale stationary Rossby wave was triggered through the interaction with atmospheric circulation.The lower tropospheric response is baroclinic and thus favo rs upward emanation of wave activity flux in the negative height anomaly area;w hile,in upper levels,the energy is dispersed to East Asia through teleconnection,the internal varieties keeping the energy,and then affects the climate in this area.

    sea ice,sea ice concentration,summer rainfall,sensitivity test,control test,transient process,wave activity flux

    1006-9895(2011)05-0847-16

    P461

    A

    張若楠,武炳義.2011.北半球大氣對(duì)春季北極海冰異常響應(yīng)的數(shù)值模擬 [J].大氣科學(xué),35(5):847-862. Zhang Ruonan,Wu Bingyi.2011.The Northern Hemisphere atmospheric response to sp ring A rctic sea ice anomalies in CAM 3.0model[J].Chinese Journalof A tmospheric Sciences(in Chinese),35(5):846-862.

    2010-09-26,2011-03-29收修定稿

    公益性行業(yè) (氣象)科研專(zhuān)項(xiàng) GYHY200906017、GYHY200706005,國(guó)家自然科學(xué)基金資助項(xiàng)目40875052、40921003,中國(guó)氣象科學(xué)研究院資助項(xiàng)目“氣候系統(tǒng)模式”

    張若楠,女,1986年出生,碩士研究生,研究方向:中高緯外強(qiáng)迫對(duì)大氣環(huán)流的影響。E-mail:xiaoruo5201314@126.com

    武炳義,E-mail:w by@cams.cma.gov.cn

    猜你喜歡
    斜壓北極海位勢(shì)
    含Hardy位勢(shì)的非線(xiàn)性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類(lèi)帶強(qiáng)制位勢(shì)的p-Laplace特征值問(wèn)題
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    含變號(hào)位勢(shì)的ρ-Kirchhoff型方程組無(wú)窮多個(gè)高能量解的存在性
    含位勢(shì)的非線(xiàn)性雙調(diào)和方程解的存在性
    Rossby波的線(xiàn)性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    南、北極海冰的長(zhǎng)期變化趨勢(shì)及其與大氣環(huán)流的聯(lián)系
    中部型El Nino與北極海冰變化的聯(lián)系
    累積海冰密集度及其在認(rèn)識(shí)北極海冰快速變化的作用
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    三级毛片av免费| 老司机在亚洲福利影院| 亚洲人成电影免费在线| 日本精品一区二区三区蜜桃| 欧美成人午夜精品| 日韩欧美在线乱码| 老司机午夜福利在线观看视频| 欧美色视频一区免费| 18禁美女被吸乳视频| 国产精品av视频在线免费观看| www.自偷自拍.com| 欧美日韩瑟瑟在线播放| 天堂影院成人在线观看| 99久久精品国产亚洲精品| 制服丝袜大香蕉在线| 午夜a级毛片| 三级国产精品欧美在线观看 | 我的老师免费观看完整版| 波多野结衣巨乳人妻| 人妻久久中文字幕网| 久久九九热精品免费| 日日干狠狠操夜夜爽| 88av欧美| 男女下面进入的视频免费午夜| 一区二区三区高清视频在线| 又紧又爽又黄一区二区| 日日夜夜操网爽| 久久草成人影院| www.www免费av| 露出奶头的视频| 日韩精品免费视频一区二区三区| 三级国产精品欧美在线观看 | 正在播放国产对白刺激| 欧美一级毛片孕妇| 一级作爱视频免费观看| 最近视频中文字幕2019在线8| 日韩 欧美 亚洲 中文字幕| 欧美日韩黄片免| 欧美日韩瑟瑟在线播放| 日韩欧美在线乱码| 首页视频小说图片口味搜索| 久久精品成人免费网站| 天天添夜夜摸| 国产成年人精品一区二区| 国产爱豆传媒在线观看 | 国产不卡一卡二| 999精品在线视频| 好看av亚洲va欧美ⅴa在| 亚洲人成网站高清观看| 免费在线观看亚洲国产| 制服诱惑二区| 欧美日韩亚洲综合一区二区三区_| 18禁国产床啪视频网站| 国产av在哪里看| 99久久久亚洲精品蜜臀av| 亚洲av五月六月丁香网| 搡老熟女国产l中国老女人| 正在播放国产对白刺激| 校园春色视频在线观看| 在线观看日韩欧美| 亚洲精华国产精华精| 叶爱在线成人免费视频播放| av国产免费在线观看| 午夜福利在线在线| 久久国产精品影院| 男男h啪啪无遮挡| 高潮久久久久久久久久久不卡| 中出人妻视频一区二区| 久99久视频精品免费| 免费看十八禁软件| 精品免费久久久久久久清纯| 露出奶头的视频| 国产精品av视频在线免费观看| 亚洲自拍偷在线| 午夜福利在线在线| 国产一级毛片七仙女欲春2| 特级一级黄色大片| 男男h啪啪无遮挡| 夜夜躁狠狠躁天天躁| 99热这里只有精品一区 | 午夜福利高清视频| 亚洲午夜精品一区,二区,三区| 日日摸夜夜添夜夜添小说| 精品久久久久久久久久免费视频| 淫妇啪啪啪对白视频| 日韩中文字幕欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 久久人妻av系列| 99精品久久久久人妻精品| 天堂动漫精品| 欧美日韩亚洲综合一区二区三区_| 天堂√8在线中文| 五月玫瑰六月丁香| 日本精品一区二区三区蜜桃| 变态另类丝袜制服| 国产精品爽爽va在线观看网站| 他把我摸到了高潮在线观看| 色噜噜av男人的天堂激情| 人人妻人人澡欧美一区二区| 国产亚洲精品第一综合不卡| 少妇人妻一区二区三区视频| 久久久国产欧美日韩av| 啦啦啦免费观看视频1| 国产精品久久久人人做人人爽| av欧美777| 少妇被粗大的猛进出69影院| www.自偷自拍.com| 12—13女人毛片做爰片一| 国产午夜精品论理片| 啦啦啦韩国在线观看视频| 国产成人一区二区三区免费视频网站| 亚洲精品久久国产高清桃花| 在线观看美女被高潮喷水网站 | 91麻豆av在线| 免费在线观看成人毛片| 欧美中文日本在线观看视频| 色综合欧美亚洲国产小说| 精品久久蜜臀av无| 亚洲欧美日韩高清专用| 老鸭窝网址在线观看| 啪啪无遮挡十八禁网站| 亚洲性夜色夜夜综合| 国产午夜精品久久久久久| 日韩欧美精品v在线| 久久久久久九九精品二区国产 | 亚洲熟妇熟女久久| 99热只有精品国产| 国产精品亚洲美女久久久| 国产一级毛片七仙女欲春2| 国产精华一区二区三区| 97超级碰碰碰精品色视频在线观看| 久热爱精品视频在线9| 天堂av国产一区二区熟女人妻 | 777久久人妻少妇嫩草av网站| 激情在线观看视频在线高清| 麻豆一二三区av精品| 日韩大尺度精品在线看网址| 国产99久久九九免费精品| 日本一本二区三区精品| 九色国产91popny在线| 精品久久久久久久久久久久久| 人妻久久中文字幕网| 最新在线观看一区二区三区| 老司机午夜福利在线观看视频| 香蕉国产在线看| 黄片小视频在线播放| 美女黄网站色视频| 国产精品,欧美在线| 999久久久国产精品视频| 久热爱精品视频在线9| 真人做人爱边吃奶动态| 99精品在免费线老司机午夜| 午夜精品一区二区三区免费看| 国产aⅴ精品一区二区三区波| 国产久久久一区二区三区| 国产成人欧美在线观看| 非洲黑人性xxxx精品又粗又长| 久久性视频一级片| 成人午夜高清在线视频| 国产成人精品无人区| 性色av乱码一区二区三区2| 久久精品国产亚洲av高清一级| 亚洲精品久久国产高清桃花| 国产熟女午夜一区二区三区| 久9热在线精品视频| 国内精品久久久久久久电影| 欧美丝袜亚洲另类 | 全区人妻精品视频| 激情在线观看视频在线高清| 欧美乱色亚洲激情| 精品福利观看| 国产午夜福利久久久久久| 变态另类成人亚洲欧美熟女| 熟女电影av网| 欧美日韩黄片免| 久久精品人妻少妇| 天天躁狠狠躁夜夜躁狠狠躁| www日本黄色视频网| 在线免费观看的www视频| 亚洲aⅴ乱码一区二区在线播放 | 99久久综合精品五月天人人| 老司机福利观看| 久久久久国产一级毛片高清牌| 老司机午夜福利在线观看视频| 亚洲第一欧美日韩一区二区三区| 狂野欧美激情性xxxx| 亚洲av电影不卡..在线观看| 在线播放国产精品三级| 观看免费一级毛片| 亚洲国产欧美网| 国内精品一区二区在线观看| 国产亚洲精品久久久久久毛片| 18禁观看日本| 欧美在线黄色| 丰满的人妻完整版| 久久久久久大精品| 欧美zozozo另类| 国产欧美日韩一区二区精品| 国产区一区二久久| 蜜桃久久精品国产亚洲av| 一本精品99久久精品77| 亚洲 欧美一区二区三区| 天堂av国产一区二区熟女人妻 | 日本黄大片高清| 精品人妻1区二区| 国产野战对白在线观看| 国产精品乱码一区二三区的特点| 日本三级黄在线观看| 午夜精品一区二区三区免费看| 久久 成人 亚洲| 手机成人av网站| 国产亚洲av高清不卡| 老熟妇仑乱视频hdxx| 一级a爱片免费观看的视频| 麻豆国产97在线/欧美 | 免费在线观看黄色视频的| bbb黄色大片| 午夜久久久久精精品| 免费在线观看视频国产中文字幕亚洲| 不卡av一区二区三区| 免费在线观看影片大全网站| 天天躁狠狠躁夜夜躁狠狠躁| 日韩精品免费视频一区二区三区| a级毛片a级免费在线| 久久天躁狠狠躁夜夜2o2o| 久99久视频精品免费| 欧美国产日韩亚洲一区| 波多野结衣巨乳人妻| 亚洲精品一卡2卡三卡4卡5卡| 国产一级毛片七仙女欲春2| www.自偷自拍.com| 99re在线观看精品视频| 91成年电影在线观看| 丝袜人妻中文字幕| 可以免费在线观看a视频的电影网站| www国产在线视频色| 色尼玛亚洲综合影院| 午夜日韩欧美国产| 999久久久精品免费观看国产| ponron亚洲| 久久婷婷人人爽人人干人人爱| 波多野结衣高清作品| 精品国产乱码久久久久久男人| 男女视频在线观看网站免费 | 白带黄色成豆腐渣| 97碰自拍视频| 国产精品久久久人人做人人爽| 香蕉丝袜av| 欧美人与性动交α欧美精品济南到| 校园春色视频在线观看| 可以在线观看毛片的网站| 国产精品免费视频内射| 日韩精品青青久久久久久| 欧美黑人欧美精品刺激| 18禁裸乳无遮挡免费网站照片| 99久久综合精品五月天人人| 日日摸夜夜添夜夜添小说| 在线a可以看的网站| 精华霜和精华液先用哪个| 亚洲精品美女久久av网站| 一级片免费观看大全| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 免费在线观看影片大全网站| 日韩欧美三级三区| 老鸭窝网址在线观看| 日韩有码中文字幕| 欧美久久黑人一区二区| 欧美精品啪啪一区二区三区| 久久精品人妻少妇| 深夜精品福利| 丝袜人妻中文字幕| 日本在线视频免费播放| 欧美性长视频在线观看| 两个人的视频大全免费| 国产精品一及| 99国产精品一区二区蜜桃av| 久久久久久亚洲精品国产蜜桃av| 观看免费一级毛片| 国产乱人伦免费视频| 搡老岳熟女国产| 制服丝袜大香蕉在线| 亚洲人成伊人成综合网2020| 淫妇啪啪啪对白视频| 午夜激情av网站| 搡老岳熟女国产| 日韩欧美免费精品| 日韩有码中文字幕| 成人永久免费在线观看视频| 村上凉子中文字幕在线| 中亚洲国语对白在线视频| 黄色a级毛片大全视频| 国产亚洲精品av在线| 在线观看www视频免费| 日韩欧美一区二区三区在线观看| 亚洲av熟女| av片东京热男人的天堂| 亚洲av中文字字幕乱码综合| 午夜两性在线视频| 村上凉子中文字幕在线| 舔av片在线| 亚洲人成伊人成综合网2020| 国产99久久九九免费精品| 色哟哟哟哟哟哟| 国产成年人精品一区二区| 国产精华一区二区三区| 性欧美人与动物交配| 99国产精品99久久久久| 国产爱豆传媒在线观看 | x7x7x7水蜜桃| 欧美乱码精品一区二区三区| 日本五十路高清| 久久香蕉国产精品| 国产欧美日韩一区二区精品| 欧美国产日韩亚洲一区| 国语自产精品视频在线第100页| 一a级毛片在线观看| 亚洲欧美一区二区三区黑人| 国产精品亚洲av一区麻豆| 欧美绝顶高潮抽搐喷水| 美女午夜性视频免费| 午夜福利免费观看在线| 欧美成狂野欧美在线观看| 亚洲av第一区精品v没综合| 亚洲成人久久爱视频| 国产av麻豆久久久久久久| 亚洲一区二区三区色噜噜| 黄片小视频在线播放| 在线观看午夜福利视频| 男女那种视频在线观看| 一边摸一边做爽爽视频免费| 久久亚洲精品不卡| 亚洲欧美一区二区三区黑人| 免费看日本二区| netflix在线观看网站| 婷婷精品国产亚洲av在线| 国产av不卡久久| 日本一区二区免费在线视频| 精品国内亚洲2022精品成人| 亚洲熟女毛片儿| 男女做爰动态图高潮gif福利片| 精品免费久久久久久久清纯| 一卡2卡三卡四卡精品乱码亚洲| 午夜激情av网站| 色噜噜av男人的天堂激情| 一个人观看的视频www高清免费观看 | 亚洲avbb在线观看| 无限看片的www在线观看| avwww免费| 毛片女人毛片| 欧美一级a爱片免费观看看 | 国内毛片毛片毛片毛片毛片| 国产精品美女特级片免费视频播放器 | 18美女黄网站色大片免费观看| 国产又黄又爽又无遮挡在线| 国产视频一区二区在线看| 国产蜜桃级精品一区二区三区| 免费看十八禁软件| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 日本 欧美在线| 啦啦啦观看免费观看视频高清| 18禁观看日本| 亚洲无线在线观看| 欧美成狂野欧美在线观看| 中文在线观看免费www的网站 | 久久九九热精品免费| 超碰成人久久| 国产一区二区在线av高清观看| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 黄色视频不卡| 美女黄网站色视频| 97超级碰碰碰精品色视频在线观看| 黄色 视频免费看| av有码第一页| 久久久久久久久免费视频了| 久久99热这里只有精品18| 久久久精品国产亚洲av高清涩受| 黄色成人免费大全| 久久久精品欧美日韩精品| 淫妇啪啪啪对白视频| 国产精品av视频在线免费观看| 国产69精品久久久久777片 | 亚洲精品中文字幕在线视频| a在线观看视频网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av第一区精品v没综合| 91九色精品人成在线观看| 久久精品综合一区二区三区| 伦理电影免费视频| 少妇人妻一区二区三区视频| 好男人在线观看高清免费视频| 国产精品一区二区精品视频观看| 日韩欧美国产一区二区入口| 国内精品一区二区在线观看| 国产在线观看jvid| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| 香蕉av资源在线| 黄色片一级片一级黄色片| 成人高潮视频无遮挡免费网站| 天堂√8在线中文| 亚洲欧美精品综合久久99| 一区二区三区高清视频在线| 亚洲aⅴ乱码一区二区在线播放 | 男人舔奶头视频| 色精品久久人妻99蜜桃| 麻豆一二三区av精品| 亚洲国产看品久久| 欧美成狂野欧美在线观看| 国产精品乱码一区二三区的特点| 久久中文看片网| 国产日本99.免费观看| 国产视频一区二区在线看| 免费看美女性在线毛片视频| 国产高清videossex| 成人国产一区最新在线观看| 性色av乱码一区二区三区2| 精品国产乱子伦一区二区三区| 欧美av亚洲av综合av国产av| 欧美zozozo另类| 99热这里只有是精品50| 午夜免费成人在线视频| 18禁裸乳无遮挡免费网站照片| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 狂野欧美白嫩少妇大欣赏| 亚洲午夜精品一区,二区,三区| 亚洲av成人av| www.精华液| 最好的美女福利视频网| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 啦啦啦免费观看视频1| 国产熟女午夜一区二区三区| 亚洲专区字幕在线| 一二三四社区在线视频社区8| 成年女人毛片免费观看观看9| 在线观看66精品国产| 国产av不卡久久| 色av中文字幕| 久久久久久九九精品二区国产 | 91大片在线观看| 波多野结衣高清无吗| 毛片女人毛片| 欧美乱色亚洲激情| 久久精品影院6| 日韩欧美在线二视频| 91字幕亚洲| 中国美女看黄片| 又黄又粗又硬又大视频| 国产精品免费一区二区三区在线| 久久99热这里只有精品18| 19禁男女啪啪无遮挡网站| aaaaa片日本免费| 国产99白浆流出| 一个人观看的视频www高清免费观看 | 亚洲美女视频黄频| 日本一区二区免费在线视频| www日本在线高清视频| 男女床上黄色一级片免费看| 熟女少妇亚洲综合色aaa.| 国产一区二区三区视频了| 国产伦人伦偷精品视频| 日日爽夜夜爽网站| 精品久久久久久,| 给我免费播放毛片高清在线观看| 不卡av一区二区三区| 90打野战视频偷拍视频| 成人18禁高潮啪啪吃奶动态图| 亚洲国产精品sss在线观看| 99精品久久久久人妻精品| 久久久久国产精品人妻aⅴ院| 国产精品爽爽va在线观看网站| 久久国产精品人妻蜜桃| 97人妻精品一区二区三区麻豆| 久9热在线精品视频| 婷婷精品国产亚洲av| a级毛片a级免费在线| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 精品无人区乱码1区二区| 999精品在线视频| 午夜福利在线观看吧| 国产免费男女视频| 国产欧美日韩精品亚洲av| 国产成人啪精品午夜网站| 国产亚洲精品第一综合不卡| av福利片在线| 亚洲 欧美 日韩 在线 免费| 国产成人精品无人区| 亚洲va日本ⅴa欧美va伊人久久| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 色噜噜av男人的天堂激情| 美女 人体艺术 gogo| 黄色 视频免费看| 91国产中文字幕| 亚洲专区国产一区二区| 亚洲aⅴ乱码一区二区在线播放 | 欧美一区二区国产精品久久精品 | 欧美一区二区精品小视频在线| 俺也久久电影网| 黑人巨大精品欧美一区二区mp4| 精品免费久久久久久久清纯| 色综合婷婷激情| 国产一区二区激情短视频| 免费观看精品视频网站| 国产精品,欧美在线| 日本在线视频免费播放| 国产又色又爽无遮挡免费看| 狂野欧美激情性xxxx| 人人妻,人人澡人人爽秒播| 亚洲成av人片免费观看| 一二三四在线观看免费中文在| 亚洲黑人精品在线| 亚洲一区中文字幕在线| 中文字幕熟女人妻在线| 免费一级毛片在线播放高清视频| 免费高清视频大片| 久久国产精品人妻蜜桃| 国产精品一及| 日本三级黄在线观看| 婷婷六月久久综合丁香| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久久久亚洲av鲁大| 操出白浆在线播放| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 国产蜜桃级精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 又粗又爽又猛毛片免费看| 最近在线观看免费完整版| 日本黄色视频三级网站网址| 青草久久国产| 成人特级黄色片久久久久久久| 免费看十八禁软件| 亚洲欧美精品综合一区二区三区| cao死你这个sao货| 欧美精品啪啪一区二区三区| 亚洲欧美一区二区三区黑人| 国产免费男女视频| 少妇裸体淫交视频免费看高清 | 欧美av亚洲av综合av国产av| 国产精品亚洲av一区麻豆| 亚洲真实伦在线观看| 1024手机看黄色片| 一级毛片女人18水好多| 一区二区三区国产精品乱码| 久久 成人 亚洲| 黄色毛片三级朝国网站| 亚洲第一欧美日韩一区二区三区| 国产成人欧美在线观看| 亚洲 欧美一区二区三区| 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 99热只有精品国产| 久久久精品国产亚洲av高清涩受| 久久久久免费精品人妻一区二区| 久久天堂一区二区三区四区| 国产精品久久电影中文字幕| 男女下面进入的视频免费午夜| 美女扒开内裤让男人捅视频| 日本 av在线| √禁漫天堂资源中文www| 国产精品久久久人人做人人爽| 伦理电影免费视频| 美女黄网站色视频| 夜夜躁狠狠躁天天躁| 国产欧美日韩精品亚洲av| 亚洲精品av麻豆狂野| 亚洲精品国产精品久久久不卡| 黄色毛片三级朝国网站| 亚洲狠狠婷婷综合久久图片| 国产精品亚洲美女久久久| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 91国产中文字幕| 黄色 视频免费看| 国产麻豆成人av免费视频| 免费看a级黄色片| 欧美成人免费av一区二区三区| av欧美777| 日韩欧美国产在线观看| 男女床上黄色一级片免费看| 精品国内亚洲2022精品成人| 亚洲中文av在线| 国产亚洲精品综合一区在线观看 | 丁香六月欧美| 90打野战视频偷拍视频| 欧美乱码精品一区二区三区| 日日干狠狠操夜夜爽| 成年版毛片免费区| 黄片大片在线免费观看| 久久久久久人人人人人| 久久久久性生活片| 日韩成人在线观看一区二区三区| 亚洲欧美精品综合久久99| 在线永久观看黄色视频| 久久久精品大字幕| 精品久久久久久久毛片微露脸| 欧美性长视频在线观看| a级毛片在线看网站| 色哟哟哟哟哟哟| 非洲黑人性xxxx精品又粗又长| 一级a爱片免费观看的视频| www.精华液| 中文字幕人妻丝袜一区二区| 一区二区三区国产精品乱码| 一夜夜www| 小说图片视频综合网站| 久久精品aⅴ一区二区三区四区| 国产片内射在线|