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

    北太平洋東部風暴軸的時空演變特征

    2013-09-22 05:38:38朱偉軍袁凱陳懿妮
    大氣科學 2013年1期
    關鍵詞:位勢風暴太平洋

    朱偉軍 袁凱, 陳懿妮,

    1 南京信息工程大學氣象災害省部共建教育部重點實驗室,南京 210044

    2 湖北省武漢市氣象局,武漢 430040

    3 浙江省氣象局,杭州 310002

    1 引言

    北半球天氣尺度(2.5~6 d)帶通瞬變擾動方差的極大值集中分布在兩個緯向拉長的區(qū)域,位置與地面氣旋/反氣旋活動路徑一致,這就是北半球中緯度太平洋和大西洋上空著名的兩大風暴軸。由于風暴軸對應著強烈的熱量輸送、動量輸送和水汽輸送,因而對北半球大氣環(huán)流的維持和全球天氣氣候異常具有十分重要的作用(Held et al., 1989)。

    風暴軸最早由 Blackmon(1976)發(fā)現(xiàn),隨后大量的觀測分析研究(如Lau,1978,1979,1988;Trenberth,1981;Blackmon, et al., 1984a,1984b;Hoskins and Valdes,1990;朱偉軍和孫照渤,2000),揭示了風暴軸的許多氣候平均特征:(1)冬季,北半球天氣尺度的渦動活動被組織成兩大風暴軸,分別位于中緯度的太平洋和大西洋上;(2)若以 500 hPa天氣尺度渦動的位勢高度方差極值中心來表征風暴軸的強度和位置,則就冬季平均而言,在太平洋區(qū)域,最強中心位于(42.5°N,175°E)附近,強度達40 dagpm2左右,而在大西洋區(qū)域,風暴軸更強,中心強度可達56 dagpm2左右,位于(45°N,50°W)附近;(3)在風暴軸區(qū)域,天氣尺度渦動的位勢高度方差和擾動動能的最大值、渦動西風動量通量的最強輻合以及渦動對位渦的最大經向輸送,都出現(xiàn)在對流層上部300 hPa層附近,而渦動向極和向上的熱量通量最大值,則位于對流層低層850 hPa層附近。

    此外,傅剛等(2009)發(fā)現(xiàn),在北太平洋和北大西洋的對流層低層,還分別存在“副熱帶風暴軸”和“副極地風暴軸”這種雙風暴軸現(xiàn)象;風暴軸不僅北半球有,在南半球也存在,而且在南半球,由于中緯度海陸差異較小,風暴軸有其獨特的特點(Trenberth,1984;Nakamura and Shimpo,2004)。

    多年來,更多的關注點是風暴軸的變率問題,現(xiàn)已逐漸成為研究熱點。人們從月、季、年、年代等不同尺度上,考察了風暴軸的這種短期氣候異常(如 Lau,1988;鄧興秀和孫照渤,1994;胡增臻和黃榮輝,1997;朱偉軍和孫照渤,2000;Norris,2000;Chang,2001;Chang and Fu,2002,2003;Nakamura et al.,2002;Compo and Sardeshmukh,2004;丁葉風等,2006;任雪娟和張耀存,2007;韓博等,2007;Wettstein and Wallace,2010;朱偉軍和李瑩,2010),主要研究成果可以簡要概括如下:(1)北半球風暴軸存在顯著的月際和季節(jié)變化,一般而言,風暴軸在冬季月份最強,位置偏南,在夏季月份最弱,位置偏北,而春、秋兩季為過渡時期;(2)風暴軸還存在顯著的年際變化和年代際變化,如北太平洋風暴軸在20世紀80年代中期,存在由偏弱轉為偏強的年代際位相轉換;(3)EOF等分析結果顯示:在月以上各種時間尺度上,北太平洋風暴軸的這種短期氣候異常,具有相似的空間變化特征,主要有兩種變化模態(tài):第一模態(tài)是風暴軸在其氣候平均位置增強或減弱的主體一致變化型,第二模態(tài)是風暴軸特別是其中東部在氣候平均位置南北兩側振蕩的經向異常型。

    但是,以上研究基本上把整個北太平洋區(qū)域作為研究對象,而李瑩等(2010)最近研究發(fā)現(xiàn),從緯向結構上看,北太平洋風暴軸不僅存在“多中心”現(xiàn)象(傅剛等,2009),而且最強中心還經常出現(xiàn)在160°W以東地區(qū)。目前,對這種在北太平洋東部出現(xiàn)風暴軸最強中心的現(xiàn)象本身及其原因,認識都還非常有限,而風暴軸東部的異常,不僅可以直接引起附近阿留申低壓和阻塞形勢的異常(陸日宇,2001;Huang et al.,2002;刁一娜等,2004;Nakamura and Wallace,1990),而且對整個北半球大氣環(huán)流和我國天氣氣候的異常,具有重要影響。因此,本文將重點探討北太平洋東部地區(qū)風暴軸的時空演變究竟有何獨特特征。

    2 資料與方法

    文中所用資料取自歐洲中期天氣預報中心(ECMWF)的再分析資料集(ERA-40)(http://www.ecmwf.int/products/data/archive/descriptions/e4/index.html [2008-09-04]):逐日的位勢高度場(z)、風場(u、v、w)、溫度場(T)資料,水平分辨率為2.5°×2.5°,覆蓋時段為1957年9月~2002年8月。

    本文首先采用 31點對稱數(shù)字濾波器(李瑩和朱偉軍,2009),從逐日原始資料直接濾波出2.5~6 d的瞬變渦動,然后每月為一段,對每一段各自計算其方差,得到每月的月平均帶通濾波方差(以下簡稱濾波方差)。其后所用的經驗正交函數(shù)(EOF)分解、回歸分析以及顯著性檢驗等方法,原理可參見Storch and Zwiers(1999),不在此贅述。

    此外,本文冬、春、夏、秋四季分別包含12~2月、3~5月、6~8月和9~11月。

    3 多中心現(xiàn)象與東部型風暴軸的統(tǒng)計特征

    傅剛等(2009)研究北太平風暴軸三維結構時,發(fā)現(xiàn)無論冬夏季,位勢高度方差、擾動動能、向極的熱量輸送在各個層次都存在“多中心”現(xiàn)象,且這種“多中心”現(xiàn)象在夏季更加明顯,但其多中心僅為圖中直觀觀察,沒有進行客觀統(tǒng)計。下面,我們以 500 hPa位勢高度濾波方差場作為代表,對1957年12月~2001年11月期間44年528個月北太平洋區(qū)域(30°N~60°N,120°E~120°W)月平均風暴軸的多中心數(shù)目以及最強中心位置進行客觀統(tǒng)計。其中,每月月平均風暴軸的多中心數(shù)目和最強中心位置的定義為:該區(qū)域每月濾波方差最大值中心所在的經緯度即為該月風暴軸最強中心的位置;而該區(qū)域每月濾波方差次大值達到當月最強中心值的三分之一以上,且有獨立的閉合中心,則記為該區(qū)域該月風暴軸存在一個次大值中心,多中心數(shù)目即為次大值中心累計數(shù)目加上1(最強中心)。

    3.1 多中心現(xiàn)象的統(tǒng)計結果

    表1給出了1957年12月~2001年11月期間44年528個月北太平洋區(qū)域(30°N~60°N,120°E~120°W)月平均風暴軸的多中心數(shù)目的統(tǒng)計結果。由表1可以看出,在總數(shù)528個月中,有500個月出現(xiàn)兩個及兩個以上的中心,其中出現(xiàn) 2~3個中心的月數(shù)有312個,這表明,“多中心”現(xiàn)象確實普遍存在,概率高達94.7%,而且最典型的分布形勢就是存在 2~3個中心。從季節(jié)來看,春季是最容易出現(xiàn)“多中心”現(xiàn)象的季節(jié),秋季和冬季相對較少,而4個及4個以上的“多中心”現(xiàn)象則更容易出現(xiàn)在夏季。進一步分析表明,每月的多中心數(shù)還存在著強烈的年際變化,但除個別月份外,其長期變化趨勢并不顯著(圖略)。

    表1 1957年12月~2001年11月北太平洋區(qū)域(30°N~60°N,120°E~120°W)500 hPa位勢高度濾波方差多中心數(shù)目逐月統(tǒng)計結果Table 1 The statistics result of numbers for maxima of 500-hPa geopotential height variance over the North Pacific Ocean (30°N–60°N,120°E–120°W) through Dec 1957 to Nov 2001

    3.2 東部型風暴軸的統(tǒng)計結果

    圖1 1957年12月~2001年11月北太平洋500 hPa位勢高度濾波方差最強中心的逐月位置圖(圖中不同符號代表出現(xiàn)的次數(shù))Fig.1 The monthly position of maximum value of 500-hPa geopotential height variance over the North Pacific Ocean through Dec 1957 to Nov 2001(Different marks represent occurrence times)

    圖1為1957年12月~2001年11月北太平洋500 hPa位勢高度濾波方差最強中心的逐月位置。根據(jù)風暴軸最強中心的落點位置,我們還從東西方向上,將北太平洋地區(qū)分為西部(120°E~160°E)、中部(160°E~160°W)和東部(160°W~120°W)三個區(qū)域,分別統(tǒng)計了風暴軸最強中心落在這三個區(qū)域的月數(shù)(表2)。如圖1和表2顯示,雖然,風暴軸最強中心的氣候平均位置位于日界線以西附近(朱偉軍和李瑩,2010),但在總數(shù)528個月中,有 63個月位于西部(120°E~160°E),295個月位于中部(160°E~160°W),即只有358個月(約 68%)位于160°W以西,其余170個月(約32%)則位于160°W以東,東西振蕩十分顯著,這一比率與李瑩等(2010)對 44年冬季平均的統(tǒng)計結果基本一致;同時,最強中心偏西時,平均位置偏南,最強中心偏東時,平均位置偏北,南北振蕩也十分明顯。

    表2 1957年12月~2001年11月北太平洋區(qū)域(30°N~60°N,120°E~120°W)500 hPa 位勢高度濾波方差最強中心位置逐月統(tǒng)計結果Table 2 The statistics result of position for maximum value of 500-hPa geopotential height variance over the North Pacific Ocean (30°N–60°N,120°E–120°W) through Dec 1957 to Nov 2001

    本文把最強中心出現(xiàn)在 160°W 以東地區(qū)時的整個風暴軸分布結構定義為東部型風暴軸,此時整個太平洋風暴軸(120°E~120°W)呈東強西弱分布。由表2顯示,從月份上來看,7月相對最容易出現(xiàn)東部型風暴軸(21/44),1月和2月最難(10/44);從季節(jié)上來看,夏季相對最容易出現(xiàn)東部型風暴軸(52/132),冬季最難(31/132)。進一步考察發(fā)現(xiàn)(圖略),當風暴軸最強中心位于160°W以東時,由于西邊的強度一般也不弱,甚至相當,因而多出現(xiàn)“多中心”現(xiàn)象;同時,東部型風暴軸出現(xiàn)的頻次及其強度還存在著顯著的年際變化。

    4 風暴軸的垂直結構

    根據(jù)上節(jié)中表 2的統(tǒng)計結果,在 1957年 12月~2001年11月的44年里,冬季和夏季總共分別有31和52個月出現(xiàn)東部型風暴軸。本節(jié)我們將分別對這31和52個月求平均,用以代表冬、夏兩季的氣候平均狀況,并同樣從東西方向上,將北太平洋地區(qū)分為西部(120°E~160°E)、中部(160°E~160°W)和東部(160°W~120°W)三個區(qū)域,重點考察東部型風暴軸所對應的各瞬變量在冬、夏兩季和不同區(qū)域平均時的緯向—垂直結構。圖 2~6分別給出了冬夏兩季以上三個區(qū)域平均的東部型風暴軸所對應的天氣尺度位勢高度方差、擾動動能動量輸送向極熱量輸送和向上熱量輸送的緯向—垂直分布。值得注意的是,與以往研究的不同之處在于(如Lau,1978,1979;Trenberth,1981;Blackmon, et al.,1984a,1984b;Hoskins and Valdes,1990),這里我們單獨挑出了東部型風暴軸來求季節(jié)平均,當然為了對比,我們也計算了不單獨挑出東部型風暴軸來求季節(jié)平均(132個月)的結果,但限于篇幅和不重復以往研究之慮,在此從略。

    通過對比發(fā)現(xiàn),是否單獨挑出東部型風暴軸來求季節(jié)平均的緯向—垂直結構特征是有些異同的。其中,相似點為:天氣尺度渦動的位勢高度方差(圖2)、擾動動能(圖 3)的最大值以及渦動西風動量通量的最強輻合(圖4)都出現(xiàn)在對流層上部300~200 hPa層附近,而渦動向極(圖5)和向上(圖6)的熱量通量最大值則位于對流層低層850~700 hPa層附近;并且,在北太平洋中、西部(圖 2~6中的a、d和b、e),上述量位置偏南,而在北太平洋東部(圖2~6中的c、f),上述量位置偏北;在冬季(圖 2~6中的a-c),上述量位置偏南,而在夏季(圖2~6中的d-f),上述量位置偏北。

    不同點為:無論冬夏,東部型風暴軸所揭示的天氣尺度渦動的位勢高度方差(圖 2)和擾動動能(圖 3)的最大值在北太平洋東部地區(qū)最強,而渦動西風動量通量的最強輻合(圖 4)以及渦動向極(圖 5)和向上(圖 6)的熱量通量最大值則在北太平洋中部地區(qū)最強,但不單獨挑出東部型風暴軸來求季節(jié)平均的結果顯示,這些量的最大值都是一致對應在北太平洋中部地區(qū)最強。

    研究表明(如 Lau, 1978, 1979; Hoskins and Valdes, 1990; Chang et al., 2002),斜壓性的強弱對風暴軸的維持和發(fā)展起到了至關重要的作用,如果不單獨挑出東部型風暴軸來求季節(jié)平均的話,那么在北太平洋中部地區(qū),由于斜壓性最強,因而與斜壓性密切聯(lián)系的渦動向極和向上熱量通量最強,進而風暴軸的強度也最強;而在北太平洋160°W以東地區(qū),由于斜壓性較弱,一般認為是風暴軸衰亡的地區(qū),因而渦動向極和向上熱量通量也較弱,進而強度也較弱(圖略)?,F(xiàn)在,東部型風暴軸的結果顯示,在北太平洋160°W以東地區(qū),風暴軸的強度(天氣尺度位勢高度方差和擾動動能)可以最強,但與斜壓性密切聯(lián)系的渦動向極和向上熱量通量的最大值卻并不是最強,因此,是何原因可以導致該地區(qū)出現(xiàn)風暴軸最強中心或東部型風暴軸,機理有待進一步診斷研究。

    5 風暴軸的EOF分析特征

    圖2 1957~2001年氣候平均的北太平洋西部(a、d)、中部(b、e)以及東部(c、f)區(qū)域平均的天氣尺度位勢高度濾波方差的緯度—高度剖面圖:(a–c)冬季;(d–f)夏季。等值線間隔:8.0 dagpm2Fig.2 Pressure–latitude cross sections of climatologically and regionally mean synoptic geopotential height variance over (a, d) western,(b, e) central, and (c,f) eastern North Pacific storm track during 1957–2001: (a–c) Winter; (d–f) summer.Contour interval: 8.0 dagpm2

    圖3 同圖2,但為擾動動能,等值線間隔:8.0 m2 s–2Fig.3 Same as Fig.2, but for disturbed kinetic energy at contour interval of 8.0 m2 s–2

    圖4 同圖2,但為動量輸送,等值線間隔:2.0 m2 s–2Fig.4 Same as Fig.2, but for momentum transport at contour interval of 2.0 m2 s–2

    目前,利用經驗正交函數(shù)(EOF)分解方法來研究風暴軸之變化的研究雖然已有很多(如Lau,1988;朱偉軍和孫照渤,2000;Wettstein and Wallace,2010;朱偉軍和李瑩,2010),但基本上都是針對整體的風暴軸來進行的。而由第3節(jié)的討論結果顯示,北太平洋風暴軸在水平結構上是有差異的,特別是在其東部區(qū)域,風暴軸有1/3概率可以出現(xiàn)最強中心,因此,能否可以通過對比不同區(qū)域的EOF結果來揭示一些新的事實正是本節(jié)討論的目的。

    圖5 同圖2,但為向極熱量輸送,等值線間隔:1.0(a–c)和0.5(d–f) K m s–1Fig.5 Same as Fig.2, but for meridional heat transport at contour intervals of (a–c) 1.0 and (d– f) 0.5 K m s–1

    圖6 同圖2,但為向上熱量輸送,等值線間隔:0.02 (a–c) 和0.01 (d–f) K Pa s–1Fig.6 Same as Fig.2, but for vertical heat transport at contour intervals of (a–c) 0.02 and (d–f) 0.01 K Pa s–1

    下面,我們以500 hPa位勢高度濾波方差場作為代表來對北太平洋風暴軸進行 EOF分析??紤]到風暴軸的東西擺動和季節(jié)南北跳躍以及對比需要,本文 EOF分解的區(qū)域為:北太平洋地區(qū)西部(30°N~60°N,120°E~160°E)、中部(30°N~60°N,160°E~160°W)、東部(30°N~60°N,160°W~120°W)和整體(30°N~60°N,120°E~120°W)四個區(qū)域;時間為:1957年12月~2001年11月44年的冬季(132個月)和夏季(132個月),此時沒有單獨強調東部型風暴軸;1957年12月~2001年11月44年里出現(xiàn)東部型風暴軸的冬季(31個月)和夏季(52個月),此時僅單獨強調了東部型風暴軸。

    表3和表4分別給出了冬季和夏季EOF分解的前三個模態(tài)的解釋方差,因為前兩個模態(tài)已解釋了大部分的方差,所以后文的討論僅限于前兩個模態(tài)。圖 7~9(圖 10~12)則給出了冬季(夏季)EOF結果的前兩個模態(tài)分布及其對應的時間系數(shù)以及回歸結果,圖中陰影區(qū)分別對應了各自時間段氣候平均的風暴軸位置。

    表3 冬季北太平洋各區(qū)域風暴軸EOF展開的前三個模態(tài)的解釋方差(其中括號內外的數(shù)據(jù)分別對應于有無單獨強調東部型風暴軸的結果)Table 3 Explained variance of the first three EOF modes over different areas of the North Pacific storm track during winter (Values in/out of the brackets indicate the situation with/without only considering of the eastern type of storm tracks)

    5.1 冬季狀況

    下面先來考察冬季狀況(圖 7、8)。在沒有單獨強調東部型風暴軸的情況下,EOF結果的第一模態(tài)反映出北太平洋風暴軸在西部(圖7a1)、中部(圖7a2)和整體(圖7a4)三個區(qū)域的空間變化分布幾乎一致,即在各自的氣候平均位置處一致增強或減弱;而在東部區(qū)域,第一模態(tài)則表現(xiàn)為風暴軸在氣候平均位置南北相反的偶極子變化型(圖7a3)。與此相反,EOF結果的第二模態(tài)雖也反映出風暴軸在西部(圖 8a1)、中部(圖 8a2)和整體(圖 8a4)三個區(qū)域的空間變化分布幾乎一致,但卻是在氣候平均位置南北的反位相變化;而在東部區(qū)域,第二模態(tài)則表現(xiàn)為風暴軸在氣候平均位置處(略偏南)一致增強或減弱(圖8a3)。

    表4 夏季北太平洋各區(qū)域風暴軸EOF展開的前三個模態(tài)的解釋方差(其中括號內外的數(shù)據(jù)分別對應于有無單獨強調東部型風暴軸的結果)Table 4 Explained variance of the first three EOF modes over different areas of the North Pacific storm track during summer (Values in/out of the brackets indicate the situation with/without only considering of the eastern type of storm tracks)

    而在僅單獨強調東部型風暴軸的情況下,與沒有單獨強調東部型風暴軸的情況相比,就北太平洋西部、中部和整體三個區(qū)域來講,兩者EOF分析的結果非常相似,即第一模態(tài)反映了風暴軸在各自氣候平均位置附近一致增強或減弱的變化(圖 7b1、b2、b4),第二模態(tài)反映了風暴軸在氣候平均位置南北的反位相變化(圖8b1、b2、b4),只是單獨考慮東部型風暴軸與否,風暴軸的氣候平均位置也有相應的變化而已;而對北太平洋東部區(qū)域而言,兩者EOF分析的結果有所不同,在僅單獨強調東部型風暴軸的情況下,第一模態(tài)(圖7b3)和第二模(圖8b3)都表現(xiàn)為風暴軸在氣候平均位置南北相反的偶極子變化型,但此時第一模態(tài)似乎更突出了風暴軸在氣候平均位置北面的變化幅度,第二模態(tài)則更突出了風暴軸在氣候平均位置南面的變化幅度。

    圖9為冬季北太平洋東部區(qū)域EOF結果的前兩個模態(tài)所對應的時間系數(shù)分別回歸到500 hPa位勢高度濾波方差場的結果。通過對比圖7、圖8和圖9可以發(fā)現(xiàn),在冬季,無論單獨強調東部型風暴軸與否,北太平洋東部區(qū)域EOF結果第一模態(tài)時間系數(shù)的回歸分布(圖 9a、c)都與相應北太平洋整體區(qū)域EOF結果的第二模態(tài)分布(圖8a4、b4)非常相似,而北太平洋東部區(qū)域EOF結果第二模態(tài)時間系數(shù)的回歸分布(圖 9b、d)則都與相應北太平洋整體區(qū)域 EOF結果的第一模態(tài)分布非常相似(圖7a4、b4)。

    5.2 夏季狀況

    圖7 冬季北太平洋 (a1、b1) 西部、(a2、b2) 中部、(a3、b3) 東部和 (a4、b4) 整體區(qū)域500 hPa位勢高度濾波方差EOF展開結果的第一模態(tài)及其(a5,b5)時間系數(shù):(a1、a2、a3、a4、a5)無僅單獨強調東部型風暴軸;(b1、b2、b3、b4、b5)僅單獨強調東部型風暴軸。陰影區(qū)為各區(qū)域相應時段氣候平均的風暴軸位置,淺色(深色)陰影表示風暴軸強度>18(24)dagpm2Fig.7 The first EOF modes of 500-hPa geopotential height variance over (a1, b1) western, (a2, b2) central, (a3, b3) eastern, and (a4, b4) entire storm track during winter, and (a5, b5) their time coefficients: (a1, a2, a3, a4, a5) Without only considering of the eastern type of storm tracks; (b1, b2, b3, b4, b5) only considering of the eastern type of storm tracks.Shaded area refers to as the climatologically mean position of the relative storm track, with light (dark) shaded areas showing the track intensity larger than 18 (24) dagpm2

    圖8 同圖7,但為冬季第二模態(tài)及其時間系數(shù)Fig.8 Same as Fig.7, but for the second modes and their time coefficients during winter

    圖9 冬季北太平洋東部區(qū)域EOF結果的第一模態(tài)(a、c)和第二模態(tài)(b、d)時間系數(shù)回歸的500 hPa位勢高度濾波方差的空間分布。(a、b)無僅單獨強調東部型風暴軸;(c、d)僅單獨強調東部型風暴軸。等值線間隔:1.0 dagpm2;淺色(深色)陰影表示通過95%(99%)的顯著性檢驗Fig.9 Regressions of the time coefficients of the first (a, c) and second (b, d) EOF modes for the eastern Pacific area on the synoptic 500-hPa geopotential height variance during winter (contour interval: 1.0 dagpm2): (a, b) Without only considering of the eastern type of storm tracks; (c, d) only considering of the eastern type of storm tracks.Light (dark) shading shows correlation values with confidence level larger than 95% (99%)

    下面再來考察夏季狀況(圖10、11)。在沒有單獨強調東部型風暴軸的情況下,與冬季情況不同的是,夏季EOF結果的第一模態(tài)雖也反映出北太平洋風暴軸在西部(圖10a1)、中部(圖10a2)和整體(圖10a4)三個區(qū)域的空間變化分布幾乎一致,但此時反映的是風暴軸在氣候平均位置偏西北和偏東南方向的反位相變化,并且似乎更突出了風暴軸在氣候平均位置偏西北方向的變化幅度;而在東部區(qū)域,第一模態(tài)則表現(xiàn)為風暴軸在氣候平均位置處(略偏南)一致增強或減弱(圖 10a3)。與此相反,EOF結果的第二模態(tài)雖也反映出風暴軸在西部(圖11a1)、中部(圖11a2)和整體(圖11a4)三個區(qū)域的空間變化分布幾乎一致,即在氣候平均位置偏西北和偏東南方向的反位相變化,但此時卻似乎更突出了風暴軸在氣候平均位置偏東南方向的變化幅度;而在東部區(qū)域,第二模態(tài)則表現(xiàn)為風暴軸在氣候平均位置南北的反位相變化(略呈西北—東南向)(圖11a3)。

    而在僅單獨強調東部型風暴軸的情況下,夏季北太平洋西部、中部、東部和整體四個區(qū)域的EOF分析結果非常一致,即第一模態(tài)反映了風暴軸在各自氣候平均位置處一致增強或減弱的變化(圖10b1、b2、b3、b4),第二模態(tài)則反映了風暴軸在氣候平均位置南北的反位相變化(圖11b1、b2、b3、b4)。

    圖10 同圖7,但為夏季第一模態(tài)及其時間系數(shù),圖中淺色(深色)陰影表示風暴軸強度>6(10)dagpm2Fig.10 Same as Fig.7, but for the first modes and their time coefficients during summer with light (dark) shaded areas showing the track intensity larger than 6 (10) dagpm2

    圖11 同圖10,但為夏季第二模態(tài)及其時間系數(shù)Fig.11 Same as Fig.10, but for the second modes and their time coefficients during summer

    圖12為夏季北太平洋東部區(qū)域EOF結果的前兩個模態(tài)所對應的時間系數(shù)分別回歸到500 hPa位勢高度濾波方差場的結果。通過對比圖 10、圖 11和圖 12可以發(fā)現(xiàn),與冬季情況相似,夏季在沒有單獨強調東部型風暴軸的情況下,北太平洋東部區(qū)域EOF結果第一模態(tài)時間系數(shù)的回歸分布(圖12a)也與相應北太平洋整體區(qū)域 EOF結果的第二模態(tài)分布(圖 11a4)非常相似,而北太平洋東部區(qū)域EOF結果第二模態(tài)時間系數(shù)的回歸分布(圖12b)則與相應北太平洋整體區(qū)域 EOF結果的第一模態(tài)分布非常相似(圖 10a4)。而在單獨強調東部型風暴軸的情況下,夏季不再出現(xiàn)模態(tài)交叉對應,即北太平洋東部區(qū)域 EOF結果第一模態(tài)時間系數(shù)的回歸分布(圖12c)與相應北太平洋整體區(qū)域EOF結果的第一模態(tài)分布(圖10b4)非常相似,而北太平洋東部區(qū)域 EOF結果第二模態(tài)時間系數(shù)的回歸分布(圖12d)則與相應北太平洋整體區(qū)域EOF結果的第二模態(tài)分布非常相似(圖11b4)。

    以上結果表明,在不同季節(jié)、不同區(qū)域以及是否單獨考慮東部型風暴軸的情況下,北太平洋風暴軸的EOF分析結果雖然表現(xiàn)出了一定的差異,但都反映出在北太平洋東部區(qū)域風暴軸的變化特征有一些獨特特點,如在該區(qū)域風暴軸的主要變化模態(tài)并不一定時時與其他區(qū)域的主要變化模態(tài)一一對應。因此,是什么原因導致北太平洋東部區(qū)域風暴軸的獨特變化?是內部因素(如Lau,1978,1979;Hoskins and Valdes,1990;Chang et al.,2002;等)還是外強迫因素(如朱偉軍和孫照渤,2000;朱偉軍和李瑩,2010;等)?目前尚不十分清楚,值得在以后的研究中深入探討。

    圖12 同圖9,但為夏季。等值線間隔:(a、b)0.2 dagpm2,(c、d)0.5 dagpm2Fig.12 Same as Fig.9, but for summer with contour intervals of (a, b) 0.2 and (c, d) 0.5 dagpm2

    6 結論

    本文利用ECMWF逐日再分析資料(ERA-40),以500 hPa位勢高度濾波方差為代表,對1957年12月~2001年11月期間44年528個月北太平洋區(qū)域(30°N~60°N,120°E~120°W)月平均風暴軸的多中心數(shù)目和最強中心位置進行了客觀統(tǒng)計,在此基礎上,對北太平洋區(qū)域進行了分區(qū),通過對比各區(qū)域風暴軸的時間演變和結構變化,重點揭示了北太平洋東部地區(qū)風暴軸的時空演變特征。主要結論如下:

    (1)逐月來看,北太平洋風暴軸“多中心”現(xiàn)象普遍存在,概率高達94.7%,最典型的分布呈2~3個中心分布。從季節(jié)上來看,春季是“多中心”現(xiàn)象最容易出現(xiàn)的季節(jié),秋季和冬季相對較少,而4個及以上的“多中心”現(xiàn)象則更容易出現(xiàn)在夏季。

    (2)若把最強中心出現(xiàn)在160°W以東地區(qū)的北太平洋風暴軸定義為東部型風暴軸,那么從月份上來看,7月相對最容易出現(xiàn)東部型風暴軸,1月和2月最難;從季節(jié)上來看,夏季相對最容易出現(xiàn)東部型風暴軸,冬季最難;總的來看,出現(xiàn)東部型風暴軸的頻數(shù)大約占總頻數(shù)的三分之一。

    (3)按照本文的分區(qū),從垂直結構上看,東部型風暴軸的結果顯示,在北太平洋 160°W 以東地區(qū),風暴軸的強度(天氣尺度位勢高度方差和擾動動能)可以最強,但與斜壓性密切聯(lián)系的渦動向極和向上熱量通量的最大值卻并不是最強。

    (4)EOF分析的結果表明,在不同季節(jié)、不同區(qū)域以及是否單獨考慮東部型風暴軸的情況下,風暴軸的變化雖然表現(xiàn)出了一定的差異,但都反映出在北太平洋東部區(qū)域風暴軸的變化特征有其獨特特點。具體來講,冬季:在北太平洋風暴軸東部區(qū)域,在沒有單獨強調東部型風暴軸的情況下,第一模態(tài)表現(xiàn)為風暴軸在氣候平均位置南北相反的偶極子變化型,第二模態(tài)表現(xiàn)為風暴軸在氣候平均位置處(略偏南)一致增強或減弱,而在僅單獨強調東部型風暴軸的情況下,第一模態(tài)和第二模態(tài)都表現(xiàn)為風暴軸在氣候平均位置南北相反的偶極子變化型,但此時第一模態(tài)似乎更突出了在氣候平均位置北面的變化幅度,第二模態(tài)則更突出了在氣候平均位置南面的變化幅度;與此相反,無論有無單獨強調東部型風暴軸,此時在北太平洋風暴軸西部、中部和整體三個區(qū)域,第一模態(tài)都一致表現(xiàn)為風暴軸在各自氣候平均位置處一致增強或減弱,第二模態(tài)都表現(xiàn)為風暴軸在各自氣候平均位置南北的反位相變化。夏季:在北太平洋風暴軸東部區(qū)域,無論有無單獨強調東部型風暴軸,第一模態(tài)都表現(xiàn)為風暴軸在氣候平均位置處(略偏南)一致增強或減弱,第二模態(tài)都表現(xiàn)為風暴軸在氣候平均位置南北的反位相變化(略呈西北—東南向);而此時在北太平洋風暴軸西部、中部和整體三個區(qū)域,在僅單獨強調東部型風暴軸的情況下,其EOF分析結果與東部區(qū)域一致,而在沒有單獨強調東部型風暴軸的情況下,第一模態(tài)和第二模態(tài)都反映的是風暴軸在氣候平均位置偏西北和偏東南方向的反位相變化,但此時第一模態(tài)似乎更突出了在氣候平均位置偏西北方向的變化幅度,而第二模態(tài)則似乎更突出了在氣候平均位置偏東南方向的變化幅度。

    值得指出的是,本文雖揭示了北太平洋東部風暴軸許多獨特的時空演變特征,但對其成因和機理并沒有太多涉及,這有待于下一步進行深入研究。

    References)

    Blackmon M L.1976.A climatological spectral study of the 500 mb geopotential height of the Northern Hemisphere [J].J.Atmos.Sci., 33 (8):1607–1623.

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

    Blackmon M L, Lee Y H, Wallace J M, et al.1984b.Time variation of 500 mb height fluctuations with long, intermediate and short time scales as deduced from lag-correlation statistics [J].J.Atmos.Sci., 41 (6):981–991.

    Chang E K M.2001.GCM and observational diagnoses of the seasonal and interannual variations of the Pacific storm track during the cool season [J].J.Atmos.Sci., 58 (13): 1784–1800.

    Chang E K M, Fu Y F.2002.Interdecadal variations in Northern Hemisphere winter storm track intensity [J].J.Climate, 15 (6): 642–658.

    Chang E K M, Fu Y F.2003.Using mean flow change as a proxy to infer interdecadal storm track variability [J].J.Climate, 16 (13): 2178–2196.

    Chang E K M, Lee S, Swanson K L.2002.Storm track dynamics [J].J.Climate, 15 (16): 2163–2183.

    Compo G P, Sardeshmukh P D.2004.Storm track predictability on seasonal and decadal scales [J].J.Climate, 17 (19): 3701–3720.

    鄧興秀, 孫照渤.1994.北半球風暴軸的時間演變特征 [J].南京氣象學院學報, 17 (2): 165–170. Deng Xingxiu, Sun Zhaobo.1994.Characteristics of temporal evolution of northern storm tracks [J].Journal of Nanjing Institute of Meteorology (in Chinese), 17 (2): 165–170.

    刁一娜, 李建平, 羅德海.2004.阻塞流場與瞬變渦動相互作用的動力學研究 [J].大氣科學, 28 (6): 901–924. 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 (6): 901–924.

    丁葉風, 任雪娟, 韓博.2006.北太平洋風暴軸的氣候特征及其變化的初步研究 [J].氣象科學, 26 (3): 237–242. Ding Yefeng, Ren Xuejuan,Han Bo.2006.Primiary analysis of the climatic characteristics and variability of the North Pacific storm track [J].Scientia Meteorologica Sinica (in Chinese), 26 (3): 237–242.

    傅剛, 畢瑋, 郭敬天.2009.北太平洋風暴軸的三維空間結構 [J].氣象學報, 67 (2): 189–200. Fu Gang, Bi Wei, Guo Jingtian.2009.Threedimensional structure of storm track over the North Pacific [J].Acta Meteorologica Sinica (in Chinese), 67 (2): 189–200.

    韓博, 任雪娟, 楊修群.2007.北太平洋風暴軸異常變化特征及其與緯向風場的關系研究 [J].氣象科學, 27 (3): 237–245. Han Bo, Ren Xuejuan, Yang Xiuqun.2007.Analysis of the North Pacific storm track anomaly and the relationship with zonal wind [J].Scientia Meteorologica Sinica (in Chinese), 27 (3): 237–245.

    Held I M, Lyons S W, Nigam S.1989.Transients and the extratropical response to EL Ni?o [J].J.Atmos.Sci., 46: 163–174.

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

    胡增臻, 黃榮輝.1997.冬季熱帶西太平洋對流活動異常的年際變化及其對北太平洋風暴軸的影響 [J].大氣科學, 21 (5): 513–522. Hu Zengzhen, Huang Ronghui.1997.The interannual variation of the convective activity in the tropical West Pacific in winter and its effect on the storm track in the North Pacific [J].Chinese Journal of Atmospheric Sciences (Scientia Atmospherica Sinica) (in Chinese), 21(5): 513–522.

    Huang F, Zhou F X, Qian X D.2002.Interannual and decadal variability of the North Pacific blocking and its relationship to SST, teleconnection and storm tracks [J].Adv.Atmos.Sci., 19 (5): 807–820.

    Lau N C.1978.On the three-dimensional structure of the observed transient eddy statistics of the Northern Hemisphere wintertime circulation [J].J.Atmos.Sci., 35 (10): 1900–1923.

    Lau N C.1979.The structure and energetics of transient disturbances in the Northern Hemisphere wintertime circulation [J].J.Atmos.Sci., 36 (6):982–995.

    Lau N C.1988.Variability of the observed midlatitude storm tracks in relation to low-frequency changes in the circulation pattern [J].J.Atmos.Sci., 45 (19): 2718–2743.

    李瑩, 朱偉軍.2009.不同數(shù)字濾波方法在風暴軸研究中的性能比較 [J].大氣科學學報, 32 (4): 565–573. Li Ying, Zhu Weijun.2009.Performance of different digital filters in storm track study [J].Transactions of Atmospheric Sciences (in Chinese), 32 (4): 565–573.

    李瑩, 朱偉軍, 魏建蘇.2010.冬季北太平洋風暴軸指數(shù)的評估及其改進 [J].大氣科學, 34 (5): 1001–1010. Li Ying, Zhu Weijun, Wei Jiansu.2010.Reappraisal and improvement of winter storm track indices in the North Pacific [J].Chinese Journal of Atmospheric Sciences (in Chinese), 34 (5): 1001–1010.

    陸日宇.2001.夏季東北亞阻塞形勢維持時的天氣尺度波 [J].大氣科學,25 (3): 289–302. Lu Riyu.2001.Eddies during the blocking maintenance over the northeastern Asia in summer [J].Chinese Journal of Atmospheric Sciences (Scientia Atmospherica Sinica) (in Chinese), 25(3): 289–302.

    Nakamura H, Wallace J M.1990.Observed changes in baroclinic wave activity during the life cycles of low-frequency circulation anomalies [J].J.Atmos.Sci., 47 (9): 1100–1116.

    Nakamura H, Izumi T, Sampe T.2002.Interannual and decadal modulations recently observed in the Pacific storm track activity and East Asian winter monsoon [J].J.Climate, 15 (14): 1855–1874.

    Nakamura H, Shompo A.2004.Seasonal variations in the Southern Hemisphere storm tracks and jet streams as revealed in a reanalysis dataset [J].J.Climate, 17: 1828–1844.

    Norris J R.2000.Interannual and interdecadal variability in the storm track,cloudiness, and sea surface temperature over the summertime North Pacific [J].J.Climate, 13 (2): 422–430.

    任雪娟, 張耀存.2007.冬季 200 hPa 西太平洋急流異常與海表加熱和大氣瞬變擾動的關系探討 [J].氣象學報, 65 (4): 550–560. Ren Xuejuan, Zhang Yaocun.2007.Association of winter western Pacific jet stream anomalies at 200 hPa with ocean surface heating and atmospheric transient eddies [J].Acta Meteorologica Sinica (in Chinese), 65 (4):550–560.

    Storch H V, Zwiers F W.1999.Statistical Analysis in Climate Research [M].Cambridge: Cambridge University Press, 484pp.

    Trenberth K E.1981.Observed Southern Hemisphere eddy statistics at 500 mb: Frequency and spatial dependence [J].J.Atmos.Sci., 38 (12):2585–2605.

    Trenberth K E.1984.Interannual variability of the Southern Hemisphere circulation: Representativeness of the year of the global weather experiment [J].Mon.Wea.Rev., 43 (1): 108–123.

    Wettstein J J, Wallace J M.2010.Observed patterns of month-to-month storm-track variability and their relationship to the background flow [J].J.Atmos.Sci., 67 (5): 1420–1437.

    朱偉軍, 孫照渤.2000.冬季北太平洋風暴軸的年際變化及其與500 hPa高度以及熱帶和北太平洋海溫的聯(lián)系 [J].氣象學報, 58 (3): 309–320.Zhu Weijun, Sun Zhaobo.2000.Interannual variability of northern winter Pacific storm track and its association with 500 hPa height and tropical and northern Pacific sea surface temperature [J].Acta Meteorologica Sinica (in Chinese), 58 (3): 309–320.

    朱偉軍, 李瑩.2010.冬季北太平洋風暴軸的年代際變化特征及其可能影響機制 [J].氣象學報, 68 (4): 477–486. Zhu Weijun, Li Ying.2010.Inter-decadal variation characteristics of winter North Pacific storm tracks and its possible influencing mechanism [J].Acta Meteorologica Sinica (in Chinese), 68 (4): 477–486.

    猜你喜歡
    位勢風暴太平洋
    含Hardy位勢的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類帶強制位勢的p-Laplace特征值問題
    愛管閑事的“太平洋警察”
    決勝太平洋
    跨越太平洋的愛戀
    含變號位勢的ρ-Kirchhoff型方程組無窮多個高能量解的存在性
    含位勢的非線性雙調和方程解的存在性
    《風暴中奔跑的樹》
    散文詩(2017年12期)2017-08-23 16:07:44
    玩轉腦風暴
    可怕的大風暴
    搖籃(2016年13期)2016-02-28 20:53:17
    美女免费视频网站| 在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 国产精品98久久久久久宅男小说| 人人妻人人澡欧美一区二区| 成人三级黄色视频| 欧美av亚洲av综合av国产av| 男人舔女人下体高潮全视频| 制服丝袜大香蕉在线| 嫁个100分男人电影在线观看| 夜夜夜夜夜久久久久| 亚洲欧洲精品一区二区精品久久久| 亚洲av中文字字幕乱码综合| 欧美黄色淫秽网站| 久久久国产成人精品二区| 午夜精品在线福利| 两性夫妻黄色片| 成年女人毛片免费观看观看9| 国产激情欧美一区二区| 亚洲熟女毛片儿| 亚洲国产精品999在线| 亚洲精品美女久久久久99蜜臀| 99国产精品一区二区蜜桃av| 国产单亲对白刺激| 搡老熟女国产l中国老女人| 亚洲精品在线美女| 国产私拍福利视频在线观看| www.999成人在线观看| 亚洲aⅴ乱码一区二区在线播放 | 欧美成人午夜精品| 亚洲国产看品久久| 搡老妇女老女人老熟妇| 成年免费大片在线观看| 麻豆一二三区av精品| 久久国产精品影院| 亚洲美女视频黄频| 日韩av在线大香蕉| 91麻豆av在线| 欧美绝顶高潮抽搐喷水| 国产野战对白在线观看| 最近视频中文字幕2019在线8| 成年人黄色毛片网站| 亚洲男人的天堂狠狠| 熟妇人妻久久中文字幕3abv| 日韩成人在线观看一区二区三区| 国产伦人伦偷精品视频| 精品久久蜜臀av无| 成人高潮视频无遮挡免费网站| 欧美又色又爽又黄视频| 国产精品精品国产色婷婷| 一级作爱视频免费观看| 看免费av毛片| ponron亚洲| 免费电影在线观看免费观看| 丝袜人妻中文字幕| 正在播放国产对白刺激| 一级黄色大片毛片| 亚洲成人国产一区在线观看| 一边摸一边抽搐一进一小说| 日本一二三区视频观看| a在线观看视频网站| 国产亚洲av嫩草精品影院| 亚洲国产精品成人综合色| 我要搜黄色片| a在线观看视频网站| 亚洲男人的天堂狠狠| 亚洲av成人一区二区三| 色在线成人网| а√天堂www在线а√下载| 桃色一区二区三区在线观看| 99国产精品一区二区三区| 国产精品影院久久| 制服丝袜大香蕉在线| 夜夜夜夜夜久久久久| 精品久久久久久久人妻蜜臀av| 麻豆久久精品国产亚洲av| 久久99热这里只有精品18| cao死你这个sao货| 国产人伦9x9x在线观看| 午夜两性在线视频| 在线免费观看的www视频| 99久久无色码亚洲精品果冻| 丰满人妻熟妇乱又伦精品不卡| 欧美性长视频在线观看| 国产av一区二区精品久久| av视频在线观看入口| 欧美中文日本在线观看视频| 欧美精品啪啪一区二区三区| 在线观看一区二区三区| 麻豆av在线久日| 久久伊人香网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久国产欧美日韩av| 最近最新中文字幕大全免费视频| 亚洲avbb在线观看| 在线观看免费视频日本深夜| 麻豆国产97在线/欧美 | 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 亚洲人成伊人成综合网2020| 亚洲一区中文字幕在线| 精品乱码久久久久久99久播| 女警被强在线播放| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲| 久久久久九九精品影院| 成人精品一区二区免费| 国产一区在线观看成人免费| 一本久久中文字幕| 久久久久久久久中文| 欧美又色又爽又黄视频| 啪啪无遮挡十八禁网站| 亚洲欧美日韩高清在线视频| 久久久久亚洲av毛片大全| 黄色视频,在线免费观看| 夜夜夜夜夜久久久久| 99国产精品一区二区三区| 午夜a级毛片| 1024视频免费在线观看| 国产激情欧美一区二区| av国产免费在线观看| 大型av网站在线播放| www.精华液| 小说图片视频综合网站| 亚洲中文av在线| 亚洲成a人片在线一区二区| 天天一区二区日本电影三级| 国产成年人精品一区二区| 最近在线观看免费完整版| 啦啦啦韩国在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 男女做爰动态图高潮gif福利片| 97人妻精品一区二区三区麻豆| 久久久久久免费高清国产稀缺| 国产亚洲精品久久久久久毛片| 国内少妇人妻偷人精品xxx网站 | 国产精品久久久久久久电影 | 香蕉国产在线看| 九色国产91popny在线| 亚洲av五月六月丁香网| 变态另类丝袜制服| 欧美一级毛片孕妇| netflix在线观看网站| 亚洲国产高清在线一区二区三| 中文字幕人成人乱码亚洲影| 搡老妇女老女人老熟妇| 夜夜夜夜夜久久久久| 久久热在线av| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 欧美性猛交黑人性爽| 一本一本综合久久| 国产在线精品亚洲第一网站| 亚洲成人免费电影在线观看| 好男人电影高清在线观看| 久久人妻福利社区极品人妻图片| 久久久久久亚洲精品国产蜜桃av| 国产成+人综合+亚洲专区| 国产高清视频在线播放一区| АⅤ资源中文在线天堂| 亚洲色图 男人天堂 中文字幕| 一区二区三区国产精品乱码| 动漫黄色视频在线观看| 久久天堂一区二区三区四区| 日本一二三区视频观看| 黄色a级毛片大全视频| 少妇熟女aⅴ在线视频| 久久天堂一区二区三区四区| 成人手机av| www.熟女人妻精品国产| 国产精品久久久久久亚洲av鲁大| 免费高清视频大片| 欧美另类亚洲清纯唯美| 99久久久亚洲精品蜜臀av| 久久精品影院6| 亚洲成av人片免费观看| 9191精品国产免费久久| 亚洲成人精品中文字幕电影| 亚洲成av人片免费观看| 一级毛片女人18水好多| 国产不卡一卡二| 国内精品久久久久精免费| 夜夜看夜夜爽夜夜摸| 淫秽高清视频在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲成人久久性| 久久久久久久久久黄片| 日韩高清综合在线| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 亚洲精品美女久久av网站| 中出人妻视频一区二区| 欧美日韩精品网址| 国产一区二区三区在线臀色熟女| aaaaa片日本免费| 精品福利观看| 亚洲人成77777在线视频| 国产真实乱freesex| 国产成+人综合+亚洲专区| 一级a爱片免费观看的视频| 久久中文字幕人妻熟女| 精品久久久久久久毛片微露脸| 99国产精品一区二区三区| 国产精品久久久久久久电影 | 黑人操中国人逼视频| 丰满人妻熟妇乱又伦精品不卡| 久久热在线av| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产日韩欧美精品在线观看 | 成人精品一区二区免费| 欧美一级毛片孕妇| 91老司机精品| 最新美女视频免费是黄的| aaaaa片日本免费| 午夜激情福利司机影院| 最新在线观看一区二区三区| 久久久久国产一级毛片高清牌| 亚洲人成电影免费在线| x7x7x7水蜜桃| 成人一区二区视频在线观看| 亚洲一区二区三区不卡视频| 国产黄片美女视频| 美女大奶头视频| 亚洲国产精品999在线| 18禁美女被吸乳视频| 欧美在线一区亚洲| 精品乱码久久久久久99久播| 亚洲性夜色夜夜综合| 免费在线观看亚洲国产| 国产人伦9x9x在线观看| 一本精品99久久精品77| 在线观看66精品国产| 久久精品综合一区二区三区| 99久久精品国产亚洲精品| 国产伦一二天堂av在线观看| 黄色女人牲交| av有码第一页| 欧美大码av| 国产一区二区三区视频了| 亚洲av熟女| 婷婷丁香在线五月| 欧美日本视频| 18禁国产床啪视频网站| a在线观看视频网站| 久久久久精品国产欧美久久久| 国内少妇人妻偷人精品xxx网站 | 久久中文字幕一级| 麻豆成人午夜福利视频| 50天的宝宝边吃奶边哭怎么回事| 一级a爱片免费观看的视频| 啪啪无遮挡十八禁网站| 久久热在线av| 亚洲真实伦在线观看| 欧美日韩亚洲综合一区二区三区_| 日韩高清综合在线| 欧美一级a爱片免费观看看 | 啪啪无遮挡十八禁网站| 色av中文字幕| 天天一区二区日本电影三级| 美女黄网站色视频| aaaaa片日本免费| 亚洲va日本ⅴa欧美va伊人久久| 日韩免费av在线播放| 国产蜜桃级精品一区二区三区| 夜夜夜夜夜久久久久| 特级一级黄色大片| 国产精品一区二区三区四区久久| 在线看三级毛片| 999精品在线视频| 亚洲精品美女久久av网站| 一进一出抽搐动态| 国产亚洲精品一区二区www| 色在线成人网| 嫁个100分男人电影在线观看| 少妇熟女aⅴ在线视频| 久久亚洲精品不卡| 日韩大码丰满熟妇| 真人做人爱边吃奶动态| 色在线成人网| 久久精品成人免费网站| 熟妇人妻久久中文字幕3abv| 99久久精品热视频| 操出白浆在线播放| 日本熟妇午夜| 视频区欧美日本亚洲| 男女下面进入的视频免费午夜| 一进一出好大好爽视频| 可以在线观看的亚洲视频| 欧美日本视频| 最近视频中文字幕2019在线8| 久久久久久亚洲精品国产蜜桃av| 久久国产精品影院| 精品久久久久久久毛片微露脸| 50天的宝宝边吃奶边哭怎么回事| 人人妻,人人澡人人爽秒播| 一个人免费在线观看电影 | 精品午夜福利视频在线观看一区| 国产精品一区二区精品视频观看| 深夜精品福利| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲电影在线观看av| 亚洲美女黄片视频| 免费无遮挡裸体视频| 香蕉av资源在线| 嫩草影院精品99| 老司机福利观看| 国产精品免费视频内射| 国产成年人精品一区二区| www.熟女人妻精品国产| 999久久久精品免费观看国产| АⅤ资源中文在线天堂| 成人av在线播放网站| 日韩欧美在线二视频| 免费在线观看成人毛片| 久久久久国产一级毛片高清牌| 黄色毛片三级朝国网站| 国产乱人伦免费视频| 啦啦啦观看免费观看视频高清| 在线十欧美十亚洲十日本专区| 欧美成人午夜精品| 无遮挡黄片免费观看| 一二三四社区在线视频社区8| 亚洲成av人片免费观看| 成年女人毛片免费观看观看9| 日韩国内少妇激情av| 婷婷精品国产亚洲av在线| 免费高清视频大片| 99国产精品一区二区蜜桃av| 久久午夜综合久久蜜桃| 无限看片的www在线观看| 欧美在线黄色| 国模一区二区三区四区视频 | 国产精品 欧美亚洲| 熟女电影av网| 欧美乱妇无乱码| 国产亚洲精品av在线| 美女大奶头视频| 久久久久久亚洲精品国产蜜桃av| 日本黄大片高清| 日韩高清综合在线| 757午夜福利合集在线观看| 欧美一区二区国产精品久久精品 | 国产av麻豆久久久久久久| 又粗又爽又猛毛片免费看| 黄片大片在线免费观看| 91在线观看av| 他把我摸到了高潮在线观看| 这个男人来自地球电影免费观看| 日本成人三级电影网站| 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 精品人妻1区二区| 18禁黄网站禁片午夜丰满| 欧美精品啪啪一区二区三区| 一区福利在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产高清在线一区二区三| 国产成人影院久久av| 黑人操中国人逼视频| 搡老岳熟女国产| 69av精品久久久久久| 天堂影院成人在线观看| 欧美日韩精品网址| 丝袜人妻中文字幕| 国产亚洲精品一区二区www| 国产精品免费一区二区三区在线| 午夜免费激情av| 午夜两性在线视频| 亚洲人成77777在线视频| 午夜福利成人在线免费观看| 999久久久国产精品视频| 欧美日韩中文字幕国产精品一区二区三区| 日韩欧美在线乱码| 欧美午夜高清在线| 欧美色欧美亚洲另类二区| 日本成人三级电影网站| 大型av网站在线播放| 国产精品国产高清国产av| 哪里可以看免费的av片| 成年女人毛片免费观看观看9| 黄色 视频免费看| 国产一区二区激情短视频| 国产1区2区3区精品| 啦啦啦免费观看视频1| 在线观看免费日韩欧美大片| 亚洲专区字幕在线| 啪啪无遮挡十八禁网站| 啦啦啦韩国在线观看视频| 国产一区在线观看成人免费| 好男人在线观看高清免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 人人妻,人人澡人人爽秒播| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 哪里可以看免费的av片| 丰满人妻熟妇乱又伦精品不卡| 日韩免费av在线播放| 性色av乱码一区二区三区2| 91老司机精品| 成人三级黄色视频| 成人三级做爰电影| 99久久久亚洲精品蜜臀av| 狠狠狠狠99中文字幕| 国产在线精品亚洲第一网站| 长腿黑丝高跟| 国产精品久久电影中文字幕| 男人舔奶头视频| 在线观看免费日韩欧美大片| 999久久久国产精品视频| 国产成人精品无人区| 亚洲av中文字字幕乱码综合| 成人三级黄色视频| 午夜福利在线观看吧| 精品午夜福利视频在线观看一区| 丰满人妻熟妇乱又伦精品不卡| 午夜精品在线福利| 成人永久免费在线观看视频| 国产三级中文精品| 久久性视频一级片| 天天一区二区日本电影三级| 亚洲精品美女久久av网站| 国产1区2区3区精品| 国产真人三级小视频在线观看| 久久香蕉激情| 婷婷亚洲欧美| 99国产精品一区二区蜜桃av| 国产片内射在线| 老司机靠b影院| 99国产精品一区二区蜜桃av| 国产1区2区3区精品| 国产av不卡久久| 婷婷丁香在线五月| 此物有八面人人有两片| 国产成人精品久久二区二区91| 国产精品国产高清国产av| 国产在线精品亚洲第一网站| 五月玫瑰六月丁香| svipshipincom国产片| 精品国产美女av久久久久小说| 亚洲一区高清亚洲精品| 一区二区三区高清视频在线| 久久精品国产亚洲av高清一级| 精品日产1卡2卡| av在线播放免费不卡| 久久久精品国产亚洲av高清涩受| 夜夜看夜夜爽夜夜摸| 大型av网站在线播放| 啦啦啦观看免费观看视频高清| 久久精品aⅴ一区二区三区四区| 女人被狂操c到高潮| www.自偷自拍.com| 久久久久九九精品影院| 久久久久久人人人人人| 18美女黄网站色大片免费观看| 亚洲国产日韩欧美精品在线观看 | 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品久久久久久毛片| 久久精品国产亚洲av高清一级| netflix在线观看网站| 亚洲av成人av| xxx96com| 国产一区在线观看成人免费| 中文在线观看免费www的网站 | 成人av在线播放网站| 国产精品99久久99久久久不卡| 一本久久中文字幕| 在线a可以看的网站| 亚洲av熟女| 小说图片视频综合网站| 成人午夜高清在线视频| 久久香蕉精品热| 久久中文字幕一级| 国产精品98久久久久久宅男小说| 给我免费播放毛片高清在线观看| 老司机在亚洲福利影院| 国产精品av久久久久免费| 窝窝影院91人妻| 天天添夜夜摸| 亚洲午夜理论影院| 免费观看人在逋| 怎么达到女性高潮| 久久久精品欧美日韩精品| 国产1区2区3区精品| 亚洲欧美一区二区三区黑人| 亚洲av成人一区二区三| 亚洲无线在线观看| 夜夜躁狠狠躁天天躁| 国产成年人精品一区二区| 99热这里只有精品一区 | 亚洲午夜理论影院| 日韩av在线大香蕉| 中文字幕久久专区| 91字幕亚洲| 日本 欧美在线| 波多野结衣巨乳人妻| 国产av不卡久久| 欧美绝顶高潮抽搐喷水| 国产精华一区二区三区| 亚洲欧美精品综合一区二区三区| 国产欧美日韩一区二区精品| 全区人妻精品视频| 非洲黑人性xxxx精品又粗又长| 亚洲中文字幕一区二区三区有码在线看 | 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 99热只有精品国产| 久久人妻av系列| 嫩草影院精品99| 看免费av毛片| 日日干狠狠操夜夜爽| 国产成人av激情在线播放| 极品教师在线免费播放| 男男h啪啪无遮挡| 亚洲国产精品成人综合色| 久久久久免费精品人妻一区二区| 久久精品成人免费网站| 男男h啪啪无遮挡| а√天堂www在线а√下载| 韩国av一区二区三区四区| 1024视频免费在线观看| 精品福利观看| 一进一出抽搐动态| 黄色 视频免费看| 久久久久国产一级毛片高清牌| 一进一出抽搐gif免费好疼| 精品午夜福利视频在线观看一区| 美女午夜性视频免费| 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 精品久久久久久久人妻蜜臀av| 精品乱码久久久久久99久播| 国产午夜精品论理片| 精品久久久久久久毛片微露脸| 国产精品免费视频内射| 免费搜索国产男女视频| 亚洲 国产 在线| 国产精华一区二区三区| 少妇人妻一区二区三区视频| 亚洲欧美激情综合另类| 国产私拍福利视频在线观看| netflix在线观看网站| 91av网站免费观看| 不卡一级毛片| 桃红色精品国产亚洲av| 日韩欧美精品v在线| 一级作爱视频免费观看| 一区福利在线观看| 人妻久久中文字幕网| 19禁男女啪啪无遮挡网站| 中文字幕久久专区| cao死你这个sao货| 欧美不卡视频在线免费观看 | 日本一本二区三区精品| 香蕉丝袜av| 天天一区二区日本电影三级| 一级作爱视频免费观看| 午夜视频精品福利| 久久伊人香网站| 丁香欧美五月| 窝窝影院91人妻| 色在线成人网| 啦啦啦观看免费观看视频高清| 色综合亚洲欧美另类图片| 久久精品国产综合久久久| 日韩成人在线观看一区二区三区| 日韩大尺度精品在线看网址| 国产精品免费一区二区三区在线| 欧美绝顶高潮抽搐喷水| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 嫁个100分男人电影在线观看| 久久人妻av系列| 精品免费久久久久久久清纯| 色综合婷婷激情| 黄片大片在线免费观看| 搡老岳熟女国产| 精品欧美一区二区三区在线| 欧美黄色片欧美黄色片| 制服人妻中文乱码| 亚洲欧美一区二区三区黑人| www日本在线高清视频| 久久久国产成人免费| 999久久久精品免费观看国产| 麻豆成人午夜福利视频| 人人妻人人澡欧美一区二区| www.精华液| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 成人av在线播放网站| 日韩欧美 国产精品| 亚洲国产中文字幕在线视频| 美女午夜性视频免费| 50天的宝宝边吃奶边哭怎么回事| 亚洲黑人精品在线| 国内精品一区二区在线观看| 在线观看免费午夜福利视频| 人人妻人人澡欧美一区二区| 国产探花在线观看一区二区| 欧美黑人欧美精品刺激| 岛国在线观看网站| 亚洲色图 男人天堂 中文字幕| 国产精品影院久久| 久久久精品大字幕| 神马国产精品三级电影在线观看 | 亚洲美女黄片视频| 天堂影院成人在线观看| 日本a在线网址| 国产激情久久老熟女| 俄罗斯特黄特色一大片| 国产99白浆流出| 女警被强在线播放| av天堂在线播放| 午夜免费激情av| 啦啦啦免费观看视频1| 成人18禁高潮啪啪吃奶动态图|