• <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麻豆狂野| 精品国产一区二区三区久久久樱花| 欧美亚洲日本最大视频资源| 黑人欧美特级aaaaaa片| 这个男人来自地球电影免费观看 | 日韩大片免费观看网站| 久久久国产精品麻豆| 欧美3d第一页| 夜夜骑夜夜射夜夜干| 九色成人免费人妻av| 欧美精品一区二区免费开放| 免费高清在线观看日韩| tube8黄色片| 下体分泌物呈黄色| 制服人妻中文乱码| 在线观看免费日韩欧美大片 | 国产精品人妻久久久影院| 91久久精品电影网| 久久久国产欧美日韩av| 人人澡人人妻人| 成人黄色视频免费在线看| 精品少妇黑人巨大在线播放| 最近2019中文字幕mv第一页| 一级毛片电影观看| av播播在线观看一区| 国产精品一区二区三区四区免费观看| 久久久久久久国产电影| 高清不卡的av网站| 亚洲怡红院男人天堂| 国产白丝娇喘喷水9色精品| av在线app专区| 王馨瑶露胸无遮挡在线观看| 久久久久久伊人网av| 2022亚洲国产成人精品| 久久久久久久久久久免费av| 精品一区二区三卡| 蜜桃国产av成人99| 一区二区三区免费毛片| 午夜激情久久久久久久| 国产精品国产三级专区第一集| 十八禁网站网址无遮挡| 久久久久久久久久久丰满| 一本一本综合久久| av天堂久久9| 最黄视频免费看| 九九在线视频观看精品| 欧美一级a爱片免费观看看| 亚洲欧美一区二区三区国产| 国产高清国产精品国产三级| 国产熟女午夜一区二区三区 | 人人澡人人妻人| 日韩大片免费观看网站| 国产午夜精品久久久久久一区二区三区| 国产高清不卡午夜福利| 两个人免费观看高清视频| 亚洲美女搞黄在线观看| 国产精品一区www在线观看| 91久久精品国产一区二区成人| 国产欧美另类精品又又久久亚洲欧美| 久久久午夜欧美精品| 久久国产精品大桥未久av| 久久毛片免费看一区二区三区| 久久久久久久大尺度免费视频| 久久久久网色| 日本爱情动作片www.在线观看| 成人国语在线视频| 久久99蜜桃精品久久| 国产精品久久久久久久电影| 国产伦理片在线播放av一区| 国产高清三级在线| 晚上一个人看的免费电影| 欧美日韩av久久| 久久午夜福利片| 成人18禁高潮啪啪吃奶动态图 | 久久久久人妻精品一区果冻| 色婷婷av一区二区三区视频| 日产精品乱码卡一卡2卡三| 日韩熟女老妇一区二区性免费视频| 久久狼人影院| 中文字幕免费在线视频6| 人妻少妇偷人精品九色| 日日摸夜夜添夜夜爱| av又黄又爽大尺度在线免费看| 亚洲精品,欧美精品| 色网站视频免费| 18+在线观看网站| 精品99又大又爽又粗少妇毛片| 精品久久久噜噜| 三级国产精品片| 精品少妇内射三级| 国产极品粉嫩免费观看在线 | 日韩av不卡免费在线播放| 日韩成人伦理影院| 老熟女久久久| 日韩 亚洲 欧美在线| 在线观看免费高清a一片| 精品久久国产蜜桃| 69精品国产乱码久久久| 满18在线观看网站| 精品人妻熟女毛片av久久网站| 国产日韩欧美亚洲二区| 免费黄色在线免费观看| 大片免费播放器 马上看| 观看美女的网站| 久久综合国产亚洲精品| 中国三级夫妇交换| 国产精品秋霞免费鲁丝片| 蜜臀久久99精品久久宅男| 亚洲成人av在线免费| 亚洲,欧美,日韩| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美日韩另类电影网站| 午夜老司机福利剧场| 天堂8中文在线网| 黑人高潮一二区| 欧美bdsm另类| 免费高清在线观看视频在线观看| 大陆偷拍与自拍| 热99久久久久精品小说推荐| 一个人免费看片子| 嘟嘟电影网在线观看| 成年人午夜在线观看视频| 国产成人a∨麻豆精品| 黑人猛操日本美女一级片| 色婷婷av一区二区三区视频| 91aial.com中文字幕在线观看| 纵有疾风起免费观看全集完整版| 国产探花极品一区二区| 人妻 亚洲 视频| 新久久久久国产一级毛片| 在线亚洲精品国产二区图片欧美 | 国产av精品麻豆| 成人免费观看视频高清| 成人免费观看视频高清| 国产欧美亚洲国产| 国产探花极品一区二区| 人妻一区二区av| 少妇的逼水好多| 秋霞伦理黄片| 日本欧美视频一区| 少妇精品久久久久久久| 国产深夜福利视频在线观看| 丝袜在线中文字幕| 妹子高潮喷水视频| 日本黄色片子视频| 久久精品久久久久久久性| 99热全是精品| 亚洲五月色婷婷综合| 最近2019中文字幕mv第一页| 国产精品一区www在线观看| 大香蕉97超碰在线| 一本大道久久a久久精品| 亚洲成人手机| 嫩草影院入口| 久久韩国三级中文字幕| 啦啦啦中文免费视频观看日本| 综合色丁香网| 国产高清有码在线观看视频| 久久ye,这里只有精品| 一级毛片电影观看| 国产精品欧美亚洲77777| 大片电影免费在线观看免费| 爱豆传媒免费全集在线观看| 女人精品久久久久毛片| 亚洲综合色惰| 我的女老师完整版在线观看| 日日撸夜夜添| 菩萨蛮人人尽说江南好唐韦庄| 国产一区二区在线观看日韩| 男人操女人黄网站| 男人操女人黄网站| 91aial.com中文字幕在线观看| 婷婷色麻豆天堂久久| 这个男人来自地球电影免费观看 | 99视频精品全部免费 在线| 欧美精品亚洲一区二区| 国产片特级美女逼逼视频| 热99国产精品久久久久久7| 2021少妇久久久久久久久久久| 亚洲欧美清纯卡通| 99久国产av精品国产电影| 国产片特级美女逼逼视频| a 毛片基地| 99热全是精品| 热re99久久精品国产66热6| 嘟嘟电影网在线观看| 纯流量卡能插随身wifi吗| 三上悠亚av全集在线观看| 国产乱人偷精品视频| 国产日韩欧美亚洲二区| 午夜影院在线不卡| 秋霞伦理黄片| 国产成人一区二区在线| 免费看光身美女| 波野结衣二区三区在线| 日本-黄色视频高清免费观看| 国产免费福利视频在线观看| 我的老师免费观看完整版| 人妻夜夜爽99麻豆av| 男女啪啪激烈高潮av片| 超色免费av| 婷婷色综合大香蕉| 色婷婷av一区二区三区视频| 一本色道久久久久久精品综合| 亚洲国产欧美在线一区| 多毛熟女@视频| 成人18禁高潮啪啪吃奶动态图 | 午夜91福利影院| 亚洲精品国产av成人精品| 日韩 亚洲 欧美在线| 麻豆精品久久久久久蜜桃| 如何舔出高潮| 久久精品国产亚洲av天美| 国产高清三级在线| 赤兔流量卡办理| 不卡视频在线观看欧美| 国产高清不卡午夜福利| 亚洲国产精品一区二区三区在线| 午夜日本视频在线| 日韩伦理黄色片| 成人国语在线视频| 搡老乐熟女国产| 91久久精品国产一区二区成人| 夜夜骑夜夜射夜夜干| 日本av手机在线免费观看| 黄色怎么调成土黄色| 日日摸夜夜添夜夜爱| 又粗又硬又长又爽又黄的视频| 精品久久蜜臀av无| 国产精品偷伦视频观看了| 人人妻人人澡人人看| xxx大片免费视频| av天堂久久9| 热re99久久国产66热| 日本欧美国产在线视频| 少妇熟女欧美另类| 国产深夜福利视频在线观看| 精品一区二区三卡| 99re6热这里在线精品视频| 久久久久久久久久人人人人人人| 边亲边吃奶的免费视频| 午夜av观看不卡| 色婷婷av一区二区三区视频| 日韩av免费高清视频| 人人妻人人添人人爽欧美一区卜| 丰满少妇做爰视频| 午夜激情福利司机影院| 另类精品久久| 少妇精品久久久久久久| 日韩成人伦理影院| 欧美日韩综合久久久久久| 日韩人妻高清精品专区| av福利片在线| 各种免费的搞黄视频| av一本久久久久| 国产成人免费无遮挡视频| 国产老妇伦熟女老妇高清| 国产69精品久久久久777片| 久久久久久久久久久丰满| 国产av一区二区精品久久| 大香蕉久久网| 母亲3免费完整高清在线观看 | 亚洲av福利一区| 国产精品一区二区在线不卡| 毛片一级片免费看久久久久| 亚洲精品日本国产第一区| 国产亚洲一区二区精品| 精品人妻一区二区三区麻豆| 日韩成人伦理影院| 美女脱内裤让男人舔精品视频| 亚洲一区二区三区欧美精品| 日韩不卡一区二区三区视频在线| 少妇的逼好多水| 免费观看无遮挡的男女| 丝袜脚勾引网站| 蜜桃久久精品国产亚洲av| 精品一品国产午夜福利视频| 国产精品.久久久| 亚洲av欧美aⅴ国产| 永久网站在线| 我的老师免费观看完整版| 国产精品一国产av| 色吧在线观看| 久久影院123| 能在线免费看毛片的网站| 在线免费观看不下载黄p国产| 91精品国产九色| 一级爰片在线观看| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 亚洲一区二区三区欧美精品| 日韩不卡一区二区三区视频在线| 一边亲一边摸免费视频| 最近中文字幕高清免费大全6| 亚洲国产精品成人久久小说| 国产深夜福利视频在线观看| 国产高清不卡午夜福利| 高清av免费在线| 亚洲精品国产av成人精品| 人妻系列 视频| 亚洲第一av免费看| 夫妻性生交免费视频一级片| 免费高清在线观看视频在线观看| 一本色道久久久久久精品综合| 校园人妻丝袜中文字幕| 老司机亚洲免费影院| 中文字幕制服av| 免费看不卡的av| 精品熟女少妇av免费看| 3wmmmm亚洲av在线观看| 制服人妻中文乱码| 在现免费观看毛片| 久久久久久人妻| 国产毛片在线视频| 欧美精品一区二区大全| 亚洲av.av天堂| 一区二区三区四区激情视频| 99热6这里只有精品| 又黄又爽又刺激的免费视频.| 国产欧美亚洲国产| 日韩av在线免费看完整版不卡| av免费在线看不卡| 中文字幕av电影在线播放| 亚洲精品成人av观看孕妇| 免费少妇av软件| 亚洲人成77777在线视频| 成人黄色视频免费在线看| 男女高潮啪啪啪动态图| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 国产伦理片在线播放av一区| 欧美精品一区二区大全| 18禁裸乳无遮挡动漫免费视频| 91久久精品国产一区二区成人| 性色av一级| 久久这里有精品视频免费| 大陆偷拍与自拍| 在线观看免费日韩欧美大片 | 黄色欧美视频在线观看| 国产白丝娇喘喷水9色精品| 国产黄频视频在线观看| 高清av免费在线| 亚洲国产最新在线播放| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| .国产精品久久| 99热国产这里只有精品6| av国产精品久久久久影院| 在线观看www视频免费| 欧美老熟妇乱子伦牲交| 亚洲精品视频女| 国产精品欧美亚洲77777| 在现免费观看毛片| 国产色爽女视频免费观看| 亚洲av二区三区四区| 亚洲色图 男人天堂 中文字幕 | 在线精品无人区一区二区三| 只有这里有精品99| 国产精品女同一区二区软件| 97精品久久久久久久久久精品| 久久久久久久亚洲中文字幕| 99视频精品全部免费 在线| 26uuu在线亚洲综合色| 黄色怎么调成土黄色| 99热全是精品| av视频免费观看在线观看| av网站免费在线观看视频| 精品久久久久久久久av| 女性被躁到高潮视频| 国产色婷婷99| 精品国产一区二区三区久久久樱花| 日韩伦理黄色片| 免费观看无遮挡的男女| 精品99又大又爽又粗少妇毛片| 欧美激情国产日韩精品一区| 一区二区三区精品91| 久久ye,这里只有精品| 在线免费观看不下载黄p国产| 秋霞伦理黄片| a级毛片在线看网站| 伦理电影免费视频| 麻豆乱淫一区二区| 精品一品国产午夜福利视频| 国产永久视频网站| 男女高潮啪啪啪动态图| 狂野欧美激情性bbbbbb| 亚洲人成网站在线观看播放| 国产黄色免费在线视频| 制服人妻中文乱码| 国产精品一二三区在线看| 国产精品不卡视频一区二区| 蜜臀久久99精品久久宅男| 国产精品一区二区在线不卡| 一边摸一边做爽爽视频免费| 尾随美女入室| 少妇丰满av| 美女内射精品一级片tv| 久久久久久久久久久丰满| 99国产综合亚洲精品| 18禁裸乳无遮挡动漫免费视频| 欧美 亚洲 国产 日韩一| 黄色欧美视频在线观看| 99精国产麻豆久久婷婷| 亚洲av男天堂| 国产亚洲精品久久久com| 国产一区有黄有色的免费视频| 狠狠婷婷综合久久久久久88av| 97超视频在线观看视频| 99久久精品国产国产毛片| 丝袜喷水一区| 男女免费视频国产| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 青春草视频在线免费观看| 人妻夜夜爽99麻豆av| 国产男人的电影天堂91| 水蜜桃什么品种好| 国产伦理片在线播放av一区| 国产 一区精品| 精品少妇内射三级| 精品久久久噜噜| 伊人亚洲综合成人网| av天堂久久9| 婷婷成人精品国产| 国产免费现黄频在线看| 一本—道久久a久久精品蜜桃钙片| 一本一本综合久久| 国产av一区二区精品久久| 国产精品偷伦视频观看了| 啦啦啦中文免费视频观看日本| 精品99又大又爽又粗少妇毛片| 久久久久久久久久人人人人人人| 亚洲精品乱久久久久久| 亚洲国产精品成人久久小说| 80岁老熟妇乱子伦牲交| 免费人妻精品一区二区三区视频| 在线观看美女被高潮喷水网站| 亚洲精品中文字幕在线视频| 少妇高潮的动态图| 一本色道久久久久久精品综合| 免费日韩欧美在线观看| 日韩一区二区视频免费看| av在线观看视频网站免费| 91在线精品国自产拍蜜月| 欧美xxxx性猛交bbbb| 日本欧美国产在线视频| 久久久精品免费免费高清| 久久免费观看电影| 人人澡人人妻人| 国产乱人偷精品视频| 亚洲av日韩在线播放| 99九九线精品视频在线观看视频| 自线自在国产av| 满18在线观看网站| 国产高清国产精品国产三级| 18禁动态无遮挡网站| 欧美三级亚洲精品| 久久影院123| 国产在线一区二区三区精| 两个人的视频大全免费| 久久综合国产亚洲精品| 久久久久久人妻| 三级国产精品欧美在线观看| 在线观看国产h片| 亚洲精品乱码久久久v下载方式| 日本黄大片高清| 在线观看美女被高潮喷水网站| 乱人伦中国视频| 22中文网久久字幕| 校园人妻丝袜中文字幕| 乱人伦中国视频| 纵有疾风起免费观看全集完整版| 久久青草综合色| 丝袜喷水一区| 免费不卡的大黄色大毛片视频在线观看| 卡戴珊不雅视频在线播放| 日韩视频在线欧美| 观看av在线不卡| 超碰97精品在线观看| 欧美日韩一区二区视频在线观看视频在线| 18禁在线播放成人免费| 丁香六月天网| 免费观看a级毛片全部| 下体分泌物呈黄色| 欧美最新免费一区二区三区| 亚洲国产精品一区三区| 亚洲第一av免费看| 我的女老师完整版在线观看| 在线观看人妻少妇| 成人国产麻豆网| 如何舔出高潮| 精品熟女少妇av免费看| 国产成人av激情在线播放 | 考比视频在线观看| 少妇人妻 视频| 狂野欧美激情性bbbbbb| 国产白丝娇喘喷水9色精品| 日韩av免费高清视频| 欧美成人精品欧美一级黄| 中文天堂在线官网| 美女福利国产在线| 久热这里只有精品99| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻精品综合一区二区| 亚洲人成77777在线视频| 中文字幕最新亚洲高清| 国产亚洲精品第一综合不卡 | 成人亚洲欧美一区二区av| 丰满乱子伦码专区| 亚洲国产精品国产精品| 亚洲精品久久久久久婷婷小说| 午夜老司机福利剧场| 午夜福利,免费看| 亚州av有码| 飞空精品影院首页| 欧美精品一区二区大全| 卡戴珊不雅视频在线播放| 亚洲人成网站在线观看播放| 国产伦精品一区二区三区视频9| 久久精品夜色国产| 国产欧美日韩一区二区三区在线 | 美女内射精品一级片tv| 日韩在线高清观看一区二区三区| 一本色道久久久久久精品综合| 亚洲精品中文字幕在线视频| 满18在线观看网站| 热re99久久国产66热| 在线播放无遮挡| 精品午夜福利在线看| 亚洲人成网站在线播| 国产精品偷伦视频观看了| 丝袜脚勾引网站| videosex国产| 97在线视频观看| 午夜福利视频在线观看免费| 男人操女人黄网站| 日韩 亚洲 欧美在线| 欧美日韩综合久久久久久| 色视频在线一区二区三区| 两个人免费观看高清视频| 欧美日韩成人在线一区二区| 亚洲精品久久午夜乱码| 80岁老熟妇乱子伦牲交| 日本vs欧美在线观看视频| 婷婷色av中文字幕| 看十八女毛片水多多多| 久久久久久伊人网av| 亚洲精品456在线播放app| 看十八女毛片水多多多| 美女国产视频在线观看| 晚上一个人看的免费电影| 91久久精品电影网| 美女cb高潮喷水在线观看| 五月开心婷婷网| 在线观看美女被高潮喷水网站| 秋霞伦理黄片| 久久久久视频综合| 蜜桃久久精品国产亚洲av| 日本av免费视频播放| 亚洲精品第二区| 亚洲av日韩在线播放| 久久午夜福利片| 久久鲁丝午夜福利片| 一级二级三级毛片免费看| 久久婷婷青草| 亚洲精品色激情综合| 亚洲不卡免费看| 午夜av观看不卡| a级毛色黄片| 亚洲精品亚洲一区二区| 人妻少妇偷人精品九色| 亚洲欧美一区二区三区黑人 | 久久国产精品男人的天堂亚洲 | 大香蕉久久成人网| 午夜激情av网站| 午夜福利网站1000一区二区三区| 国语对白做爰xxxⅹ性视频网站| 人妻制服诱惑在线中文字幕| 人人澡人人妻人| 日韩成人伦理影院| 搡老乐熟女国产| 国产精品人妻久久久久久| 午夜久久久在线观看| 国产日韩欧美在线精品| 女人精品久久久久毛片| 久久久久久久久久久丰满| 国产极品粉嫩免费观看在线 | 久久久久久久亚洲中文字幕| 亚洲国产精品专区欧美| 久久精品国产亚洲av天美| 久久精品国产自在天天线| 亚洲精品中文字幕在线视频| 这个男人来自地球电影免费观看 | 日本猛色少妇xxxxx猛交久久| 亚洲第一av免费看| 18禁裸乳无遮挡动漫免费视频| 亚洲不卡免费看| 久久久久精品久久久久真实原创| 欧美 日韩 精品 国产| 制服诱惑二区| 亚洲av综合色区一区| 成人黄色视频免费在线看| 丝袜在线中文字幕| 制服诱惑二区| 如何舔出高潮| 欧美bdsm另类| 中文字幕免费在线视频6| 最黄视频免费看| 高清毛片免费看| 日韩精品有码人妻一区| 狠狠精品人妻久久久久久综合| 观看av在线不卡| 午夜精品国产一区二区电影| 少妇被粗大猛烈的视频|