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

    北大西洋濤動位相轉(zhuǎn)換的數(shù)值研究

    2012-12-15 03:00:54李忠燕羅德海
    大氣科學(xué) 2012年2期
    關(guān)鍵詞:正位西風(fēng)流場

    李忠燕 羅德海

    1 中國海洋大學(xué)海洋環(huán)境學(xué)院,青島 266100

    2 中國科學(xué)院大氣物理研究所東亞區(qū)域氣候—環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100029

    北大西洋濤動位相轉(zhuǎn)換的數(shù)值研究

    李忠燕1羅德海2

    1 中國海洋大學(xué)海洋環(huán)境學(xué)院,青島 266100

    2 中國科學(xué)院大氣物理研究所東亞區(qū)域氣候—環(huán)境重點(diǎn)實(shí)驗(yàn)室,北京 100029

    借助一個非線性準(zhǔn)地轉(zhuǎn)正壓模型,模擬了基本西風(fēng)氣流的強(qiáng)弱、不同空間分布的天氣尺度擾動渦對北大西洋濤動 (NAO)位相轉(zhuǎn)換的作用,以及大尺度雙波地形對其的影響。通過一系列的對比試驗(yàn)發(fā)現(xiàn),適當(dāng)弱的基本西風(fēng)氣流有利于下游系統(tǒng)西退使NAO由正 (負(fù))位相轉(zhuǎn)換為負(fù) (正)位相;振幅強(qiáng)、活動范圍大的天氣尺度擾動波是NAO由負(fù)位相轉(zhuǎn)換為正位相的有利條件,而加強(qiáng)并向北移的天氣尺度擾動波是NAO由正位相轉(zhuǎn)換為負(fù)位相的有利條件;疊加上大尺度雙波地形后,NAO由正(負(fù))位相轉(zhuǎn)換為負(fù) (正)位相的個例數(shù)減少。

    北大西洋濤動 (NAO) 位相轉(zhuǎn)換 數(shù)值研究 天氣尺度波 西風(fēng)氣流

    1 引言

    北大西洋濤動 (North Atlantic Oscillation,簡稱NAO)是指北大西洋上冰島低壓與亞速爾高壓的氣壓變化之間存在一種南北向 “蹺蹺板”的現(xiàn)象,即當(dāng)亞速爾地區(qū)氣壓偏高時,冰島地區(qū)的氣壓偏低,反之當(dāng)冰島地區(qū)氣壓偏高時,亞速爾地區(qū)的氣壓偏低。在全球氣候系統(tǒng)中,NAO占有重要的位置。它不僅直接影響北大西洋及附近地區(qū)的氣候,甚至對整個北半球的氣候都有重要影響。沈?qū)W順和木本昌秀 (2007)研究發(fā)現(xiàn)春季歐亞大陸中高緯地表氣溫的年際變化與前期冬季大西洋濤動有著顯著的正相關(guān)。Wu et al.(2009)通過觀測和數(shù)值試驗(yàn)指出春季異常NAO對東亞夏季風(fēng)有重要影響,且NAO對近30年東亞夏季風(fēng)與ENSO之間的關(guān)系的加強(qiáng)有重要作用 (Wu et al.,2011)。通常當(dāng)NAO強(qiáng)時,格陵蘭島以西及南歐偏冷,歐洲北部和美國東部偏暖,歐洲降水偏多,而地中海地區(qū)則降水偏少,反之亦然 (Hurrell,1995a)。

    許多研究表明,對于NAO不同位相,天氣尺度波的強(qiáng)迫是一個重要的驅(qū)動機(jī)制 (Hurrell,1995b;Feldstein and Lee,1998;Feldstein,2003;Benedict et al.,2004;Vallis et al.,2004;Luo et al.,2007a,2007b)。急流區(qū)是大氣的強(qiáng)斜壓區(qū),而在準(zhǔn)定常急流下游稍北處存在千公里尺度天氣尺度瞬變擾動渦 (Blackmon et al.,1984;Hoskins and Valdes,1990),將這種渦動的帶通瞬變方差局地最大值區(qū)稱為風(fēng)暴軸。在北半球冬季渦動的振幅最大,并且這種天氣尺度的渦動帶在東太平洋和北大西洋有兩個峰值,這兩個峰值帶即通常所說的太平洋和大西洋風(fēng)暴軸。Chang and Fu(2002)利用EOF方法對51年的300hPa經(jīng)向風(fēng)場進(jìn)行診斷分析,指出冬季北大西洋和北太平洋風(fēng)暴軸在20世紀(jì)80年代中期存在由偏弱轉(zhuǎn)為偏強(qiáng)的年代際位相轉(zhuǎn)換。胡躍文和楊小怡 (2007)研究表明,NAO模態(tài)的北大西洋中心在1975年后表現(xiàn)出不同程度的東移。Luo et al.(2010b)對不同時期250hPa冬季EKE(eddy kinetic energy)空間形態(tài)進(jìn)行診斷分析,指出,1979~1990年和1998~2007年兩個時間段的渦動動能相比1958~1977年具有年代際東移的特征。

    不同位相的NAO與北大西洋急流和風(fēng)暴路徑的強(qiáng)度和位置以及經(jīng)向和緯向上水汽和熱量的輸送等相聯(lián)系,最終導(dǎo)致溫度和降水場的變化。Archambault et al.(2010)研究表明,當(dāng)NAO由正位相轉(zhuǎn)換為負(fù)位相時,美國東北地區(qū)降水顯著增加,而當(dāng)NAO由負(fù)位相轉(zhuǎn)換為正位相時,美國東北地區(qū)降水減少,所以了解NAO位相轉(zhuǎn)換機(jī)理,對預(yù)測災(zāi)害天氣有重要的指導(dǎo)作用。目前,有關(guān)NAO位相轉(zhuǎn)換的理論問題還沒有得到解決?;疚黠L(fēng)氣流、天氣尺度擾動渦以及大尺度雙波地形這些因素對NAO的流場結(jié)構(gòu)造成怎樣的影響?NAO的位相轉(zhuǎn)換是否與這些因素有關(guān)?本文運(yùn)用一個非線性準(zhǔn)地轉(zhuǎn)正壓模式數(shù)值 (Luo et al.,2010a),模擬了基本西風(fēng)氣流強(qiáng)弱、不同強(qiáng)度和不同分布的天氣尺度擾動波、大尺度雙波地形對NAO位相轉(zhuǎn)換的影響。

    2 數(shù)值模式

    2.1 模式方程

    Luo et al.(2010a)提出了一個非線性準(zhǔn)地轉(zhuǎn)正壓模型,它既能考慮基本西風(fēng)氣流對不同位相的NAO初始場的空間分布和傳播的影響,又能考慮天氣尺度瞬變擾動以及大尺度雙波地形對行星尺度流場的強(qiáng)迫作用。NAO就其形態(tài)來講通常表現(xiàn)為兩種:即NAO的正位相和NAO的負(fù)位相。本文運(yùn)用這一正壓模式研究在NAO不同位相的初始流場條件下,基本西風(fēng)氣流、天氣尺度瞬變擾動等因素對NAO位相轉(zhuǎn)換的影響。為了探討在NAO兩種不同位相的初始流場條件下,大尺度雙波地形對NAO位相轉(zhuǎn)換的強(qiáng)迫作用,本文也進(jìn)行了以下兩種模擬:(1)一個NAO正 (負(fù))初始場在不同強(qiáng)度的基本西風(fēng)氣流和不同空間分布的天氣尺度擾動波的強(qiáng)迫下經(jīng)歷一個NAO由正 (負(fù))位相轉(zhuǎn)換為NAO負(fù) (正)位相的演變過程;(2)在上述的初始流場上疊加一個大尺度雙波地形,初始場在基本西風(fēng)氣流、天氣尺度擾動波和地形的共同強(qiáng)迫下經(jīng)歷位相轉(zhuǎn)換的演變過程。從而揭示大尺度雙波地形對NAO位相轉(zhuǎn)換的可能影響。

    在β平面通道假設(shè)下,考慮地形的受渦動強(qiáng)迫的無量綱的正壓渦度方程為

    其中,Ψ是大氣總的流函數(shù),h是大尺度地形,F(xiàn)=(L/Rd)2,L和Rd分別是水平長度的特征尺度和羅斯貝變形半徑,β=β0L2/U0,是緯度為φ0的無量綱科里奧利參數(shù)的經(jīng)向梯度,F(xiàn)s是擾動渦源。在經(jīng)向方向滿足以下南北邊界條件:

    2.2 初始場

    將初始流場分解為三個部分,即基本流、行星尺度波流場和天氣尺度波流場:

    其中,Ψ是大氣總的流函數(shù);是由基本流和行星尺度波流場疊加而成的流場,即擾動的背景流場;是天氣尺度波流場 (或瞬變渦度流場);U(y)=u0是基本西風(fēng)氣流;B0=0.4是初始流場行星尺度波的無量綱振幅;k=2k0(k0=1/[6.371cos(φ0)]),m=±2π/(Ly/2),m為正 (負(fù))時代表 NAO 正(負(fù))位相,Ly分別為β平面通道的寬度;γ=0.5;初始天氣尺度波流場

    2.3 擾動場 (強(qiáng)迫場)

    在數(shù)值模式中擾動場由兩個部分組成,一是由位于背景流場上游的造波器產(chǎn)生的天氣尺度波擾動場,二是大尺度雙波地形。其中,天氣尺度波擾動場是位于背景流場上游的造波器不斷產(chǎn)生,在基本流的作用下向下游移動,形成一渦列,用以模擬風(fēng)暴軸。其控制表達(dá)式為:

    其中,b0是無量綱擾動渦源振幅;2)是波數(shù)為9和11的天氣尺度波的角頻率;ms=m/2;c0是控制擾動渦活動范圍的量,c0越小,擾動渦活動范圍越大;xs表示擾動渦源的中心位置,即位于背景流上游的位置,xs越小,擾動渦活動越偏東,這里令xs=2(表示2000km);α=±1,當(dāng)α為負(fù)時代表NAO負(fù)位相。圖1是不同振幅、不同活動范圍的天氣尺度波擾動源分布圖。

    加上大尺度雙波地形后,NAO不同位相的初始流場將在基本西風(fēng)氣流、天氣尺度擾動波和地形的共同強(qiáng)迫下發(fā)生位相轉(zhuǎn)換過程。本文所采用的大尺度雙波地形如圖2所示,其控制表達(dá)式為:

    其中,h0是無量綱地形的振幅,k=2k0(k0=1/[6.371cos(φ0)]),ms=m/2,m=±2π/(Ly/2),k和m分別表示x和y方向的波數(shù)。

    差分計(jì)算中空間差分使用的是中央差分格式,時間用蛙跳格式,使用超松馳迭代法解泊松方程。在初始值和強(qiáng)迫源條件下積分,得到NAO不同位相初始流場在強(qiáng)迫源的作用下發(fā)生位相轉(zhuǎn)換的演變情況。

    表1 無地形強(qiáng)迫條件下不同強(qiáng)度的基本西風(fēng)氣流和天氣尺度波條件下NAO由負(fù)位相轉(zhuǎn)為正位相的個例數(shù)Table 1 Without topographic forcing,the number of the cases in which the NAO transforms from negative phase to positive phase under the conditions of different intensity of basic westerly flow and synoptic-scale waves

    3 NAO負(fù)位相轉(zhuǎn)正位相的數(shù)值模擬

    3.1 無地形條件下的數(shù)值模擬

    設(shè)定h0=0,即不考慮大尺度雙波地形的影響。為了討論在NAO負(fù)位相初始條件下,來自上游的不同強(qiáng)度天氣尺度擾動波和不同強(qiáng)度基本西風(fēng)氣流對NAO由負(fù)位相轉(zhuǎn)為正位相的影響,本文進(jìn)行了多次數(shù)值模擬試驗(yàn)。按照如圖1的五種天氣尺度波分布,再疊加不同強(qiáng)度的基本西風(fēng)氣流背景流進(jìn)行數(shù)值模擬,統(tǒng)計(jì)出NAO由負(fù)位相轉(zhuǎn)換成正位相的個例數(shù),統(tǒng)計(jì)結(jié)果如表1所示。其中括號內(nèi)的數(shù)和符號“*”表示位于下游的系統(tǒng)西退至大西洋使NAO的位相轉(zhuǎn)換的情況,其余的表示從上游或局地發(fā)展使位相發(fā)生轉(zhuǎn)換的個例數(shù)。從表1可以看出,弱基本西風(fēng)、強(qiáng)(振幅大,活動范圍廣)的天氣尺度擾動波是NAO由負(fù)位相轉(zhuǎn)換為正位相的有利條件。

    圖1 天氣尺度波擾動源水平分布圖。橫坐標(biāo) (向東):x方向的無量綱長度,縱坐標(biāo) (向北):y方向的無量綱長度;等值線:無量綱流函數(shù) (等值線間距0.2),實(shí)線:異常反氣旋;虛線:異常氣旋Fig.1 Horizontal distributions of the wavemaker Fs.xcoordinate and ycoordinate represent nondimensional length;the contour represents stream function with contour interval(CI)of 0.2;the solid(dashed)contour denotes abnormal anticyclone(cyclone)

    圖3給出的是其中一例下游系統(tǒng)西退使NAO由負(fù)位相轉(zhuǎn)換為正位相總流場的演變過程,從圖中可以看出,初始時刻 (第0天)為NAO負(fù)位相的初始場,在第1天時,上游的強(qiáng)迫源不斷激發(fā)出天氣尺度擾動波,與行星尺度波相互作用,第4天流函數(shù)的振幅出現(xiàn)閉合高壓中心,到第8天環(huán)流增強(qiáng)至最大,出現(xiàn)多渦結(jié)構(gòu),到第12天位于下游的弱NAO正位相事件開始逐漸西退,到第16天時已經(jīng)西退至大西洋使原中心的NAO負(fù)位相發(fā)生轉(zhuǎn)換,接著轉(zhuǎn)換的NAO正位相發(fā)展再逐漸消亡。

    圖4給出的是對應(yīng)時間的行星尺度波流場的演變過程,從中也可以明顯地看出下游的大尺度北負(fù)南正的距平中心西退至大西洋使NAO由負(fù)位相轉(zhuǎn)為正位相的演變過程。

    圖3 模擬的一次NAO由負(fù)位相轉(zhuǎn)換為正位相的演變過程 (無量綱總流場)。等值線間距:0.2Fig.3 Simulated NAO transition process from negative phase to positive phase(total stream function field).CI=0.2

    3.2 大尺度雙波地形條件下的數(shù)值模擬

    設(shè)定h0=0.2,即考慮大尺度雙波地形的影響。大尺度雙波地形提供了兩槽兩脊的形勢。在地形的強(qiáng)迫下,類似第一部分,通過調(diào)整基本西風(fēng)氣流大小,天氣尺度擾動波的振幅、活動范圍進(jìn)行數(shù)值模擬,統(tǒng)計(jì)出NAO由負(fù)位相轉(zhuǎn)為正位相的個例數(shù),統(tǒng)計(jì)結(jié)果如表2所示。括號內(nèi)的數(shù)和符號 “*”表示位于下游的系統(tǒng)西退至大西洋使NAO的位相轉(zhuǎn)換的情況,其余的表示從上游東移或局地發(fā)展使位相發(fā)生轉(zhuǎn)換的個例數(shù)。從表2可以看出,疊加大尺度雙波地形強(qiáng)迫后,弱基本西風(fēng)、強(qiáng) (振幅大,活動范圍廣)的天氣尺度擾動波仍是NAO由負(fù)位相轉(zhuǎn)為正位相的有利條件,但相比于無地形強(qiáng)迫條件下的數(shù)值模擬,發(fā)現(xiàn)大尺度雙波地形的存在使位相轉(zhuǎn)換的個例數(shù)減少了,并且大尺度雙波地形的存在造成流場的南北不對稱形勢,高壓得到加強(qiáng),低壓被弱化,這也使得閉合型流場的時間提前,閉合高壓的強(qiáng)度略微增強(qiáng),并且地形的存在破壞了天氣尺度波的南北結(jié)構(gòu),雖也分裂成2支,但在阻塞區(qū)域北支明顯強(qiáng)于南支,行星尺度波雖出現(xiàn)弱偶極子結(jié)構(gòu),由于高壓強(qiáng)度明顯強(qiáng)于低壓、活動范圍大于低壓區(qū),以致弱偶極子減弱形成Ω結(jié)構(gòu)。刁一娜等 (2004)通過動力學(xué)研究得出地形強(qiáng)迫下會產(chǎn)生“Ω”型阻塞。

    圖4 模擬的一次NAO由負(fù)位相轉(zhuǎn)換為正位相的演變過程 (無量綱行星尺度流場)(等值線間距:0.1)Fig.4 Simulated NAO transition process from negative phase to positive phase(planetary-scale stream function field,CI=0.1)

    表2 同表1,但為地形強(qiáng)迫條件下Table 2 As in Table 1,but with topographic forcing

    表3 無地形強(qiáng)迫條件下不同強(qiáng)度的基本西風(fēng)氣流和天氣尺度波條件下NAO由正位相轉(zhuǎn)換為負(fù)位相的個例數(shù)Table 3 Without topographic forcing,the number of the cases in which the NAO transforms from positive phase to negative phase under the conditions of different intensity of basic westerly flow and synoptic-scale waves

    4 NAO正位相轉(zhuǎn)負(fù)位相的數(shù)值模擬

    4.1 無地形條件的數(shù)值模擬

    設(shè)定h0=0,即不考慮大尺度雙波地形的影響。需要指出的是,當(dāng)初始流場為NAO正位相時,在原有模式中調(diào)節(jié)基本西風(fēng)氣流或擾動源的強(qiáng)度,NAO正位相很難轉(zhuǎn)化為NAO負(fù)位相。根據(jù)歷年觀測資料分析發(fā)現(xiàn),風(fēng)暴軸瞬變擾動渦列在上世紀(jì)末前后有變動,其活動中心東移并略向北移,強(qiáng)度增強(qiáng) (Chang and Fu,2002,Luo et al.2010b),根據(jù)這個觀測結(jié)果,我們調(diào)節(jié)擾動源的分布,即在模式中將擾動源公式變?yōu)椋?/p>

    其中,令km=8,調(diào)節(jié)b0、b1、α的大小,得到不同強(qiáng)度的天氣尺度波空間分布圖 (如圖5所示)。

    根據(jù)圖5中擾動源的分布調(diào)節(jié)模式中b0、b1、α的大小,并同時調(diào)節(jié)基本流u0的大小,使其從0.5遞增到0.9,統(tǒng)計(jì)出下游阻塞西退至大西洋,使NAO由正位相轉(zhuǎn)換為NAO負(fù)位相的情況,根據(jù)西退、局地發(fā)展、東移三種不同情況進(jìn)行統(tǒng)計(jì),結(jié)果如表3所示。

    從表3可以看出,適當(dāng)弱的基本西風(fēng)氣流有利于歐洲大陸上的阻塞西退至大西洋,使NAO由正位相轉(zhuǎn)換為負(fù)位相;而當(dāng)基本西風(fēng)急流強(qiáng)時,會出現(xiàn)局地發(fā)展或從上游東移而來的阻塞使NAO由正位相轉(zhuǎn)換為負(fù)位相;強(qiáng)度加強(qiáng),活動中心東移且北移的天氣尺度擾動波是NAO由正位相轉(zhuǎn)換為負(fù)位相的有利條件。

    圖6給出的是其中一例下游阻塞西退至大西洋使NAO由正位相轉(zhuǎn)換為負(fù)位相總流場的演變過程,從圖中可以看出,在第27天中心流函數(shù)的振幅最大,到第29天位于下游的阻塞開始逐漸西退,到第33天位于下游的阻塞已經(jīng)西退至大西洋,然后開始不斷發(fā)展,到35天達(dá)到最強(qiáng),隨后又不斷減弱,到第41天,又出現(xiàn)西退使位相發(fā)生轉(zhuǎn)換。

    表4 同表3,但為地形強(qiáng)迫條件下Table 4 As in Table 3,but with topographic forcing

    圖5 不同強(qiáng)度的天氣尺度波擾動分布圖 (等值線間距:0.4,實(shí)線表示異常反氣旋,虛線表示異常氣旋)Fig.5 Horizontal distributions of the wave maker.CI=0.4;the solid(dashed)contours denote anticyclones(cyclones)

    4.2 大尺度雙波地形條件下的數(shù)值模擬

    設(shè)定h0=0.2,考慮大尺度雙波地形的影響。類似上一部分,通過調(diào)整基本西風(fēng)氣流大小以及天氣尺度擾動波的振幅,根據(jù)西退、局地發(fā)展、東移三種不同情況統(tǒng)計(jì)出NAO由正位相轉(zhuǎn)為負(fù)位相的個例數(shù),計(jì)結(jié)果如表4所示。從表4中可以得到類似的結(jié)論,適當(dāng)弱的基本西風(fēng)氣流有利于歐洲大陸上的阻塞西退至大西洋,使NAO由正位相轉(zhuǎn)換為負(fù)位相;而當(dāng)基本西風(fēng)急流強(qiáng)時,會出現(xiàn)局地發(fā)展的阻塞或從上游東移而來的阻塞使中心的NAO由正位相轉(zhuǎn)換為負(fù)位相;強(qiáng)度加強(qiáng),活動中心東移且北移的天氣尺度擾動波是NAO由正位相轉(zhuǎn)換為負(fù)位相的有利條件;但相比于無地形強(qiáng)迫條件下的數(shù)值模擬,發(fā)現(xiàn)大尺度雙波地形的存在使下游阻塞西退到大西洋使NAO位相轉(zhuǎn)換的個例數(shù)略有減少。

    圖6 模擬的一次NAO由正位相轉(zhuǎn)換為負(fù)位相的演變過程 (無量綱總流場)(等值線間距:0.4)Fig.6 Simulated NAO transition process from positive phase to negative phase(total stream function field,CI=0.4)

    5 結(jié)論與討論

    本文借助一個β通道正壓數(shù)值模型來模擬了基本西風(fēng)氣流和天氣尺度渦對NAO位相轉(zhuǎn)換的影響,并討論了大尺度雙波地形對其的影響。得到了主要結(jié)論如下:

    (1)弱的基本西風(fēng)氣流有利于下游系統(tǒng)西退至大西洋,使NAO由負(fù) (正)位相轉(zhuǎn)換為NAO正(負(fù))位相。對這一結(jié)論的理解,可采用有限振幅的羅斯貝波相速度公式 (Luo,2000;Luo et al.,2011):

    其中,Cp代表波動相速度,u0代表基本西風(fēng)基本流,k,m分別為緯向和經(jīng)向波數(shù),M0代表波動振幅大小。取δ>0以及F≈1(Luo,2000)。從式(6)中可以得到,u0越小,導(dǎo)致Cp越容易出現(xiàn)負(fù)值,即下游系統(tǒng)越容易西退至中心,而原中心的模態(tài)將不斷減弱,最終使中心的位相發(fā)生位相轉(zhuǎn)換。

    (2)振幅強(qiáng)、活動范圍大的天氣尺度擾動波是NAO由負(fù)位相轉(zhuǎn)換為正位相的有利條件;而天氣尺度擾動波的加強(qiáng)并北移時有利于NAO由正位相轉(zhuǎn)換為負(fù)位相。這是因?yàn)?,?dāng)天氣尺度波越強(qiáng),就會產(chǎn)生較大振幅的波動,即式 (6)中M0的值越大,就會導(dǎo)致Cp越容易出現(xiàn)負(fù)值;且當(dāng)天氣尺度擾動波加強(qiáng)北移時,有利于北歐高緯度地區(qū)1波加強(qiáng)西退,最終使中心的位相發(fā)生位相轉(zhuǎn)換。

    (3)大尺度雙波地形的存在使NAO由正 (負(fù))位相轉(zhuǎn)換為負(fù)(正)個例數(shù)略有減少。這是因?yàn)榇蟪叨入p波地形的存在使中高緯度的基本西風(fēng)氣流增強(qiáng) (Luo et al.,2007a),從而使得式 (6)中Cp越容易出現(xiàn)正值,不利于下游系統(tǒng)西退。

    需要指出的是,除了基本流、天氣尺度渦以及地形會對NAO位相轉(zhuǎn)換的有影響,在大氣中海溫、急流以及大氣斜壓性等多種要素對NAO位相轉(zhuǎn)換也有著重要的影響。因此,綜合其它因素對NAO位相轉(zhuǎn)換又有著怎樣的影響,這些仍需要進(jìn)一步的研究。

    致謝兩位審稿人對本文的修改提出了許多寶貴的意見,使作者受益匪淺。在此,作者對兩位審稿人表示由衷的感謝。

    (References)

    Archambault H M,Keyser D,Bosart L F.2010.Relationships between large-scale regime transitions and major cool-season precipitation events in the Northeastern United States[J].Mon.Wea.Rev.,138:3454-3473.

    Benedict J J,Lee S,F(xiàn)eldstein S B.2004.Synoptic view of the North Atlantic Oscillation[J].J.Atmos.Sci.,61:121-144.

    Blackmon M L,Lee Y H,Wallace J M.1984.Horizontal structure of 500mb height fluctuations with long,intermediate and short time scales[J].J.Atmos.Sci.,41:961-979.

    Chang E K M,F(xiàn)u Y F.2002.Interdecadal variations in Northern Hemisphere winter storm track intensity [J].J.Climate,15:642-658.

    刁一娜,李建平,羅德海.2004.阻塞流場與瞬變渦動相互作用的動力學(xué)研究 [J].大氣科學(xué),28:901-923. Diao Yina,Li Jianping,Luo Dehai.2004.A dynamic study of the interaction between transient eddies and blocking[J].Chinese Journal of Atmospheric Sciences(in Chinese),28:901-923.

    Feldstein S B.2003.The dynamics of NAO teleconnection pattern growth and decay[J].Quart.J.Roy.Meteor.Soc.,129:901-924.

    Feldstein S B,Lee S.1998.Is the atmospheric zonal index driven by an eddy feedback?[J].J.Atmos.Sci.,55:3077-3086.

    Hoskins B J,Valdes P J.1990.On the existence of storm-tracks[J].J.Atmos.Sci.,47:1854-1864.

    胡躍文,楊小怡.2007.北極濤動與北大西洋濤動的低頻變化特征[J].氣象科學(xué),27(3):216-322. Hu Yuewen,Yang Xiaoyi.2007.Low-frequency variability of AO and NAO [J].Scientia Meteorologica Sinica(in Chinese),27(3):216-322.

    Hurrell J W.1995a.Decadal trends in the North Atlantic Oscillation:regional temperatures and precipitation[J].Science,269:676-679.

    Hurrell J W.1995b.Transient eddy forcing of the rotational flow during northern winter[J].J.Atmos.Sci.,52:2286-2301.

    Luo D H.2000.Planetary-scale baroclinic envelope Rossby solitons in a two-layer model and their interaction with synoptic-scale eddies[J].Dyn.Atmos.Oceans,32:27-74.

    Luo D H,Diao Y N,F(xiàn)eldstein S B.2011.The variability of the Atlantic storm track and the North Atlantic Oscillation:A link between intraseasonal and interannual variability [J].J.Atmos.Sci.,68:577-601.

    Luo D,Gong T,Diao Y.2007a.Dynamics of eddy-driven low-frequency dipole modes.Part III:Meridional displacement of westerly jet anomalies during two phases of NAO [J].J.Atmos.Sci.,64:3232-3248.

    Luo D H,Lupo A R,Wan H.2007b.Dynamics of eddy-driven lowfrequency dipole modes.Part I:A simple model of North Atlantic Oscillation[J].J.Atmos.Sci.,64:3-38.

    Luo D H,Zhong L H,Ren R C,et al.2010a.Spatial pattern and zonal shift of the North Atlantic Oscillation.Part II:Numerical experiments[J].J.Atmos.Sci.,67:2827-2853.

    Luo D H,Zhu Z H,Ren R C,et al.2010b.Spatial pattern and zonal shift of the North Atlantic Oscillation.Part I:A dynamical interpretation[J].J.Atmos.Sci.,67:2805-2826.

    沈?qū)W順,木本昌秀.2007.春季歐亞大陸地表氣溫變化特征的氣候意義 [J].大氣科學(xué),31(1):19-27. Shen Xueshun,Kimoto M.2007.Studies of the interannual variability of springtime Eurasian surface air temperature[J].Chinese Journal of Atmospheric Sciences(in Chinese),31(1):19-27.

    Vallis G K,Gerber E P,Kushner P J,et al.2004.A mechanism and simple dynamical model of the North Atlantic Oscillation and annular modes[J].J.Atmos.Sci.,61:264-280.

    Wu Z W,Li J P,Jiang Z H,et al.2011.Possible effects of the North Atlantic Oscillation on the strengthening relationship between the East Asian summer monsoon and ENSO [J].International Journal of Climatology,31.

    Wu Z W,Wang B,Li J P,et al.2009.An empirical seasonal prediction model of the East Asian summer monsoon using ENSO and NAO [J].J.Geophys.Res.,114:D18120,doi:10.1029/2009JD011733.

    A Numerical Study of the Phase Transition of the North Atlantic Oscillation

    LI Zhongyan1and LUO Dehai2

    1CollegeofPhysicalandEnvironmentalOceanography,OceanUniversityofChina,Qingdao266100
    2KeyLaboratoryofRegionalClimate-EnvironmentforEastAsia,InstituteofAtmosphericPhysics,ChineseAcademyofSciences,Beijing100029

    Using a nonlinear quasi-geostrophic barotropic model,the paper investigates numerically how the spatial pattern of synoptic scale waves and the westerly wind strength affect the phase transition of the North Atlantic Oscillation(NAO),and the role played by the large-scale topography.It is found that the weak basic westerly wind is favorable for the westward shift of the downstream large-scale system,favoring the transition of the NAO from positive(negative)phase into negative(positive)phase.It is noted that when the synoptic scale waves are stronger and concentrated in a wider region,the NAO event can transform more easily from the negative phase into positive phase;when the synoptic scale waves are stronger and shift northward,the NAO can transform more easily from the positive into negative phase.In the presence of a large scale topography,the number of the cases in which the NAO transforms from positive(negative)phase to negative(positive)phase is reduced.

    North Atlantic Oscillation(NAO),phase transition;numerical study,synoptic scale waves,westerly wind strength

    1006-9895(2012)02-0374-11

    P434

    A

    10.3878/j.issn.1006-9895.2011.11055

    李忠燕,羅德海.2012.北大西洋濤動位相轉(zhuǎn)換的數(shù)值研究 [J].大氣科學(xué),36(2):374-384,

    10.3878/j.issn.1006-9895.2011.11055.Li Zhongyan,Luo Dehai.2012.A numerical study of the phase transition of the North Atlantic Oscillation[J].Chinese Journal of Atmospheric Sciences(in Chinese),36(2):374-384.

    2011-03-14,2011-06-28收修定稿

    國家自然科學(xué)基金資助項(xiàng)目41075042,中國科學(xué)院 “百人計(jì)劃”資助項(xiàng)目y163011

    李忠燕,女,1986年出生,碩士,主要從事北大西洋濤動研究。E-mail:lizhongyan@ouc.edu.cn

    猜你喜歡
    正位西風(fēng)流場
    滌綸正位/換位紡復(fù)合結(jié)構(gòu)紗的試制及性能評價
    古道西風(fēng)(水彩)
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計(jì)算
    關(guān)于鐵路貨車制動管系漏泄原因分析及對策建議
    秋漸濃,分付西風(fēng)此夜涼
    文苑(2019年20期)2019-11-16 08:52:46
    浣溪沙 誰念西風(fēng)獨(dú)自涼
    飛魔幻A(2019年3期)2019-06-12 09:22:35
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    金融委的意義在于為金融“正位”
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計(jì)分析
    需要給合作學(xué)習(xí)“正位”
    亚洲av美国av| 日韩电影二区| 精品一品国产午夜福利视频| 午夜福利乱码中文字幕| 日韩欧美一区二区三区在线观看 | 成年动漫av网址| 国产91精品成人一区二区三区 | 成人免费观看视频高清| 成年人免费黄色播放视频| 丝袜美腿诱惑在线| 天堂中文最新版在线下载| 最新的欧美精品一区二区| 久久香蕉激情| 国产成人啪精品午夜网站| 久久ye,这里只有精品| 日日爽夜夜爽网站| 超碰成人久久| 日韩一区二区三区影片| 久久久久久人人人人人| 午夜福利视频在线观看免费| 国产野战对白在线观看| 搡老岳熟女国产| 国产不卡av网站在线观看| 在线精品无人区一区二区三| 国产av精品麻豆| 在线看a的网站| 免费高清在线观看视频在线观看| 女警被强在线播放| 亚洲成国产人片在线观看| www.熟女人妻精品国产| 国产99久久九九免费精品| 国产成+人综合+亚洲专区| 精品少妇一区二区三区视频日本电影| 国产成人欧美| 日本欧美视频一区| 国产伦理片在线播放av一区| 国产亚洲欧美在线一区二区| 日本a在线网址| 日韩人妻精品一区2区三区| 如日韩欧美国产精品一区二区三区| 在线av久久热| 久久国产精品人妻蜜桃| 午夜91福利影院| 欧美日韩精品网址| 12—13女人毛片做爰片一| 欧美+亚洲+日韩+国产| 国产精品香港三级国产av潘金莲| 1024香蕉在线观看| www.999成人在线观看| 成年美女黄网站色视频大全免费| 成人手机av| 两人在一起打扑克的视频| 男女国产视频网站| 国产1区2区3区精品| 国产精品99久久99久久久不卡| 国产高清视频在线播放一区 | 色综合欧美亚洲国产小说| 大码成人一级视频| 国产在线视频一区二区| 99国产极品粉嫩在线观看| av天堂在线播放| 国产视频一区二区在线看| 欧美亚洲日本最大视频资源| 午夜福利视频在线观看免费| 黄片小视频在线播放| 国产有黄有色有爽视频| 中国美女看黄片| 国产日韩一区二区三区精品不卡| 爱豆传媒免费全集在线观看| 中国国产av一级| 欧美日韩国产mv在线观看视频| 黄网站色视频无遮挡免费观看| 天天躁夜夜躁狠狠躁躁| 国产免费av片在线观看野外av| av有码第一页| cao死你这个sao货| 午夜精品久久久久久毛片777| 色婷婷av一区二区三区视频| 精品一区在线观看国产| 男男h啪啪无遮挡| 下体分泌物呈黄色| 老熟妇仑乱视频hdxx| 国产在线观看jvid| 99精品欧美一区二区三区四区| 天天添夜夜摸| 老鸭窝网址在线观看| 视频区图区小说| 丰满人妻熟妇乱又伦精品不卡| 精品高清国产在线一区| 欧美精品亚洲一区二区| 亚洲国产欧美日韩在线播放| 日韩欧美国产一区二区入口| 蜜桃在线观看..| 狠狠婷婷综合久久久久久88av| 亚洲av欧美aⅴ国产| 日韩有码中文字幕| www.熟女人妻精品国产| av电影中文网址| 啦啦啦啦在线视频资源| 欧美在线黄色| 亚洲中文av在线| 色精品久久人妻99蜜桃| 亚洲免费av在线视频| 黄片播放在线免费| 一本大道久久a久久精品| 天堂俺去俺来也www色官网| 国产精品久久久av美女十八| 高清视频免费观看一区二区| 日韩欧美免费精品| av福利片在线| 亚洲欧洲精品一区二区精品久久久| 一本久久精品| 国产亚洲一区二区精品| 国产成人系列免费观看| 日韩欧美一区二区三区在线观看 | 脱女人内裤的视频| 美女国产高潮福利片在线看| 少妇粗大呻吟视频| 免费少妇av软件| 国产欧美日韩精品亚洲av| 欧美在线一区亚洲| 在线观看免费午夜福利视频| 精品一区在线观看国产| 一本久久精品| 亚洲精品中文字幕在线视频| av视频免费观看在线观看| 岛国毛片在线播放| 亚洲精品中文字幕一二三四区 | 国产成人欧美在线观看 | 啪啪无遮挡十八禁网站| 人人妻人人澡人人看| √禁漫天堂资源中文www| 91九色精品人成在线观看| 成人av一区二区三区在线看 | 中国国产av一级| 老司机深夜福利视频在线观看 | 99香蕉大伊视频| 人妻久久中文字幕网| 欧美乱码精品一区二区三区| 欧美日韩黄片免| 国产在线观看jvid| 成在线人永久免费视频| 成人黄色视频免费在线看| 亚洲国产中文字幕在线视频| 国产精品久久久久成人av| 纯流量卡能插随身wifi吗| 国产精品久久久久成人av| 在线观看www视频免费| 日韩电影二区| 欧美日韩亚洲综合一区二区三区_| 美女主播在线视频| www.熟女人妻精品国产| 亚洲国产精品成人久久小说| 新久久久久国产一级毛片| 国产在线免费精品| 国产黄色免费在线视频| 日韩大片免费观看网站| 99国产精品99久久久久| 12—13女人毛片做爰片一| 免费在线观看完整版高清| 一级,二级,三级黄色视频| 国产亚洲欧美精品永久| 亚洲欧美清纯卡通| 精品人妻1区二区| 人人妻人人澡人人爽人人夜夜| 免费观看a级毛片全部| 搡老乐熟女国产| 午夜激情av网站| 男女无遮挡免费网站观看| 国产伦人伦偷精品视频| 国产熟女午夜一区二区三区| 久久热在线av| 久久久久国产精品人妻一区二区| 在线观看一区二区三区激情| 大香蕉久久网| tocl精华| 搡老乐熟女国产| 天堂俺去俺来也www色官网| 女性被躁到高潮视频| 国产无遮挡羞羞视频在线观看| 亚洲精品一二三| 亚洲国产欧美一区二区综合| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 999久久久精品免费观看国产| 久久久欧美国产精品| cao死你这个sao货| 亚洲国产欧美在线一区| 国产免费视频播放在线视频| 国产精品一二三区在线看| 亚洲成人国产一区在线观看| 久久99一区二区三区| 成人国语在线视频| 国产成人免费无遮挡视频| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人欧美精品刺激| 两性夫妻黄色片| 亚洲欧洲精品一区二区精品久久久| 老司机影院成人| 美女中出高潮动态图| 精品国产乱码久久久久久男人| 国产精品秋霞免费鲁丝片| 国产成人a∨麻豆精品| www.自偷自拍.com| av有码第一页| 99久久精品国产亚洲精品| 国产免费视频播放在线视频| 亚洲国产成人一精品久久久| 黄色视频不卡| 成年人黄色毛片网站| 中文字幕最新亚洲高清| 99热全是精品| 精品人妻在线不人妻| 免费观看人在逋| 成人国产av品久久久| 伦理电影免费视频| 亚洲精品国产av成人精品| 亚洲伊人色综图| 91av网站免费观看| 中文字幕色久视频| 欧美精品啪啪一区二区三区 | 欧美精品一区二区大全| 欧美一级毛片孕妇| 国产精品熟女久久久久浪| 美女中出高潮动态图| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| 美女国产高潮福利片在线看| 久久久欧美国产精品| 久久天躁狠狠躁夜夜2o2o| 丰满人妻熟妇乱又伦精品不卡| 97精品久久久久久久久久精品| 久久精品国产亚洲av香蕉五月 | 国产又爽黄色视频| 亚洲一码二码三码区别大吗| 国产精品免费大片| 老司机亚洲免费影院| 国产真人三级小视频在线观看| 一级,二级,三级黄色视频| 久久久精品94久久精品| 久久国产精品影院| kizo精华| 中文欧美无线码| 一级毛片女人18水好多| 黄色视频在线播放观看不卡| 亚洲国产看品久久| 伊人亚洲综合成人网| 制服人妻中文乱码| 妹子高潮喷水视频| 中文字幕人妻丝袜一区二区| 后天国语完整版免费观看| 日本黄色日本黄色录像| 又大又爽又粗| 亚洲精品国产色婷婷电影| 欧美日韩精品网址| 国产91精品成人一区二区三区 | 欧美黄色淫秽网站| 五月天丁香电影| av欧美777| a级毛片在线看网站| 久久国产精品人妻蜜桃| 国产一级毛片在线| 日韩熟女老妇一区二区性免费视频| 日韩一卡2卡3卡4卡2021年| 国产亚洲精品一区二区www | 热re99久久国产66热| 99久久人妻综合| 少妇裸体淫交视频免费看高清 | 超色免费av| 亚洲欧美日韩高清在线视频 | 日韩精品免费视频一区二区三区| 午夜福利在线观看吧| 两个人看的免费小视频| 黄色 视频免费看| av有码第一页| 女人精品久久久久毛片| 99国产精品一区二区三区| 少妇被粗大的猛进出69影院| 十八禁人妻一区二区| 久久 成人 亚洲| 午夜福利视频在线观看免费| 欧美在线一区亚洲| 欧美中文综合在线视频| 亚洲三区欧美一区| 19禁男女啪啪无遮挡网站| 亚洲人成电影观看| 高潮久久久久久久久久久不卡| 欧美日韩黄片免| 久久国产精品大桥未久av| 大片电影免费在线观看免费| 国产亚洲av高清不卡| 亚洲久久久国产精品| 黄色片一级片一级黄色片| 丝袜美足系列| 美女午夜性视频免费| 亚洲精品日韩在线中文字幕| 一区二区三区四区激情视频| 丰满饥渴人妻一区二区三| 少妇的丰满在线观看| 国产高清videossex| 亚洲欧美激情在线| av天堂在线播放| 国产精品av久久久久免费| 国产成人欧美在线观看 | 日韩一卡2卡3卡4卡2021年| 欧美日韩精品网址| 国产在线免费精品| 亚洲精品乱久久久久久| 丝袜脚勾引网站| 人妻一区二区av| 女人精品久久久久毛片| 黄色视频在线播放观看不卡| 一本一本久久a久久精品综合妖精| 在线观看免费午夜福利视频| 中文字幕精品免费在线观看视频| a级片在线免费高清观看视频| 久久精品久久久久久噜噜老黄| 中文字幕精品免费在线观看视频| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 国产一区二区 视频在线| 韩国精品一区二区三区| 国产成人精品无人区| 一级毛片女人18水好多| tube8黄色片| 国产精品偷伦视频观看了| 午夜日韩欧美国产| 在线观看一区二区三区激情| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 美女中出高潮动态图| 男人操女人黄网站| 99国产精品一区二区三区| 国产精品国产av在线观看| 一区在线观看完整版| 又黄又粗又硬又大视频| 国产一卡二卡三卡精品| 国产精品秋霞免费鲁丝片| 水蜜桃什么品种好| 十八禁网站免费在线| 亚洲av男天堂| 亚洲熟女毛片儿| 欧美日韩中文字幕国产精品一区二区三区 | 午夜成年电影在线免费观看| 不卡一级毛片| 国产精品av久久久久免费| 中文字幕最新亚洲高清| 99久久人妻综合| 丝袜人妻中文字幕| 久久人人爽av亚洲精品天堂| 美女视频免费永久观看网站| 国产成人系列免费观看| 亚洲美女黄色视频免费看| 又黄又粗又硬又大视频| 亚洲av美国av| 精品久久蜜臀av无| 动漫黄色视频在线观看| 最近中文字幕2019免费版| 国产亚洲欧美精品永久| 久久久久久久大尺度免费视频| 欧美激情极品国产一区二区三区| 欧美另类亚洲清纯唯美| 一级a爱视频在线免费观看| 成年人免费黄色播放视频| 777米奇影视久久| 在线亚洲精品国产二区图片欧美| 欧美激情高清一区二区三区| 国产成人啪精品午夜网站| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲综合一区二区三区_| 欧美黄色淫秽网站| 亚洲精品在线美女| 日本a在线网址| 久久青草综合色| 女人高潮潮喷娇喘18禁视频| 久久人妻熟女aⅴ| √禁漫天堂资源中文www| 自线自在国产av| 色精品久久人妻99蜜桃| 日日爽夜夜爽网站| 他把我摸到了高潮在线观看 | 狂野欧美激情性xxxx| 国产视频一区二区在线看| 中文字幕人妻丝袜制服| 亚洲av电影在线观看一区二区三区| 伊人久久大香线蕉亚洲五| 国产免费福利视频在线观看| 国产av国产精品国产| 国产精品国产av在线观看| 少妇人妻久久综合中文| 国产精品国产三级国产专区5o| 欧美激情极品国产一区二区三区| 精品欧美一区二区三区在线| 老熟妇仑乱视频hdxx| 一级,二级,三级黄色视频| 搡老熟女国产l中国老女人| kizo精华| 大码成人一级视频| 老汉色∧v一级毛片| 成年美女黄网站色视频大全免费| av在线播放精品| 午夜福利乱码中文字幕| 人人妻人人澡人人看| 建设人人有责人人尽责人人享有的| 欧美xxⅹ黑人| 啦啦啦中文免费视频观看日本| 丁香六月天网| 最黄视频免费看| 美女国产高潮福利片在线看| 国产一区二区三区在线臀色熟女 | 丝袜脚勾引网站| 视频区欧美日本亚洲| 亚洲欧洲精品一区二区精品久久久| av在线老鸭窝| 久久香蕉激情| 国产精品久久久久久精品电影小说| 欧美日本中文国产一区发布| 国产男人的电影天堂91| 自线自在国产av| 亚洲一卡2卡3卡4卡5卡精品中文| 丰满饥渴人妻一区二区三| cao死你这个sao货| 人妻人人澡人人爽人人| 深夜精品福利| 色婷婷av一区二区三区视频| 色婷婷久久久亚洲欧美| 高清视频免费观看一区二区| 精品一区在线观看国产| 国产精品久久久av美女十八| 两人在一起打扑克的视频| 国产欧美日韩一区二区三区在线| 满18在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 国产精品.久久久| 久久午夜综合久久蜜桃| 久久狼人影院| 一本—道久久a久久精品蜜桃钙片| 性色av一级| 男女下面插进去视频免费观看| 高清欧美精品videossex| 五月开心婷婷网| 高清在线国产一区| 久久狼人影院| 首页视频小说图片口味搜索| 黑人巨大精品欧美一区二区蜜桃| 热99国产精品久久久久久7| 女人久久www免费人成看片| 91老司机精品| 午夜福利一区二区在线看| 日本撒尿小便嘘嘘汇集6| 老熟妇仑乱视频hdxx| 制服诱惑二区| 午夜激情av网站| 最近最新中文字幕大全免费视频| 亚洲欧美清纯卡通| 欧美日韩中文字幕国产精品一区二区三区 | 美女视频免费永久观看网站| 欧美亚洲 丝袜 人妻 在线| 老司机亚洲免费影院| 大陆偷拍与自拍| 亚洲色图综合在线观看| 成年女人毛片免费观看观看9 | 黄色视频不卡| 精品国产乱子伦一区二区三区 | 国产高清videossex| avwww免费| 成年女人毛片免费观看观看9 | av福利片在线| 国产成人欧美| 成人手机av| 动漫黄色视频在线观看| 超色免费av| 国产精品国产三级国产专区5o| 精品免费久久久久久久清纯 | 青草久久国产| 国产精品免费视频内射| 亚洲国产精品一区二区三区在线| 咕卡用的链子| 91成年电影在线观看| 在线观看免费高清a一片| 久久久久精品人妻al黑| 久久人人爽av亚洲精品天堂| 欧美黄色片欧美黄色片| 制服人妻中文乱码| 亚洲国产精品999| 丰满迷人的少妇在线观看| 久久久欧美国产精品| 日韩中文字幕视频在线看片| 欧美精品高潮呻吟av久久| 两人在一起打扑克的视频| 黑人巨大精品欧美一区二区蜜桃| 午夜老司机福利片| 亚洲精品在线美女| 日本av免费视频播放| 国产一级毛片在线| 久9热在线精品视频| 黄色 视频免费看| 中文字幕人妻熟女乱码| 婷婷成人精品国产| 两性夫妻黄色片| 丝袜在线中文字幕| 操出白浆在线播放| 国产深夜福利视频在线观看| 免费在线观看黄色视频的| 久久人人97超碰香蕉20202| 别揉我奶头~嗯~啊~动态视频 | 后天国语完整版免费观看| 国产精品一二三区在线看| 欧美精品一区二区免费开放| 欧美亚洲日本最大视频资源| 欧美日韩一级在线毛片| 成人av一区二区三区在线看 | 国产欧美日韩一区二区三区在线| 免费在线观看黄色视频的| 日韩制服丝袜自拍偷拍| 热99re8久久精品国产| 新久久久久国产一级毛片| 成人国语在线视频| 巨乳人妻的诱惑在线观看| 人妻久久中文字幕网| 久久亚洲精品不卡| 多毛熟女@视频| 18禁国产床啪视频网站| 我的亚洲天堂| 19禁男女啪啪无遮挡网站| 国产成人精品久久二区二区免费| 俄罗斯特黄特色一大片| 欧美国产精品va在线观看不卡| 51午夜福利影视在线观看| 欧美精品啪啪一区二区三区 | 淫妇啪啪啪对白视频 | 夜夜夜夜夜久久久久| 国产日韩欧美在线精品| √禁漫天堂资源中文www| 黄片小视频在线播放| 午夜免费鲁丝| 视频区欧美日本亚洲| 丝袜美足系列| 美女高潮到喷水免费观看| 亚洲国产欧美网| 久久人妻福利社区极品人妻图片| 成年人午夜在线观看视频| 中国美女看黄片| 999精品在线视频| 日韩制服骚丝袜av| 精品亚洲成国产av| 国产成人啪精品午夜网站| 久久毛片免费看一区二区三区| 中文字幕最新亚洲高清| 久久av网站| 国产高清国产精品国产三级| 国产伦理片在线播放av一区| 亚洲精品av麻豆狂野| 久久人妻熟女aⅴ| 又紧又爽又黄一区二区| av福利片在线| 久久久国产欧美日韩av| 国产激情久久老熟女| 亚洲国产精品999| 免费一级毛片在线播放高清视频 | 久久av网站| 黄色怎么调成土黄色| 国产一区二区三区在线臀色熟女 | 午夜福利在线观看吧| 另类精品久久| 黄色片一级片一级黄色片| 一个人免费看片子| 丰满人妻熟妇乱又伦精品不卡| 亚洲久久久国产精品| netflix在线观看网站| 亚洲一码二码三码区别大吗| 午夜两性在线视频| 国产欧美日韩综合在线一区二区| 久久精品国产a三级三级三级| 国产精品亚洲av一区麻豆| 美女午夜性视频免费| 久久影院123| 岛国毛片在线播放| 美国免费a级毛片| 午夜福利在线观看吧| 12—13女人毛片做爰片一| 成人国产av品久久久| 成年动漫av网址| 久久午夜综合久久蜜桃| 美女高潮到喷水免费观看| 国产成+人综合+亚洲专区| 一级,二级,三级黄色视频| 亚洲国产成人一精品久久久| 狂野欧美激情性xxxx| 国产91精品成人一区二区三区 | 老司机在亚洲福利影院| 色精品久久人妻99蜜桃| 一级黄色大片毛片| 中文字幕人妻熟女乱码| 国产精品久久久av美女十八| 久久ye,这里只有精品| 亚洲中文av在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲av欧美aⅴ国产| 丝袜在线中文字幕| 亚洲精品乱久久久久久| 欧美老熟妇乱子伦牲交| 国产成人免费无遮挡视频| 免费观看人在逋| 久久国产精品影院| 国产熟女午夜一区二区三区| 丝袜喷水一区| 免费av中文字幕在线| 亚洲精品粉嫩美女一区| 中文字幕精品免费在线观看视频| 久久天躁狠狠躁夜夜2o2o| 动漫黄色视频在线观看| 一区二区av电影网| 亚洲精品久久成人aⅴ小说| 国产av国产精品国产|