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

    南半球風(fēng)暴路徑的加強(qiáng)與南北擴(kuò)張?

    2018-01-18 05:48:00張榮望
    關(guān)鍵詞:南半球急流西風(fēng)

    張榮望, 黃 菲,2? ?

    (1. 中國(guó)海洋大學(xué)物理海洋教育部重點(diǎn)實(shí)驗(yàn)室 青島海洋科學(xué)與技術(shù)國(guó)家實(shí)驗(yàn)室,山東 青島 266100; 2. 寧波大學(xué)寧波市非線性海洋和大氣災(zāi)害系統(tǒng)協(xié)同創(chuàng)新中心,浙江 寧波 315000)

    風(fēng)暴路徑(Storm track)是中緯度天氣尺度瞬變擾動(dòng)的最強(qiáng)信號(hào)之一,通常由位勢(shì)高度場(chǎng)周期不大于一周的帶通濾波方差來(lái)表征[1]。風(fēng)暴路徑在北半球和南半球均有存在,反映的是大尺度斜壓不穩(wěn)定波的活動(dòng),對(duì)天氣和氣候系統(tǒng)有著重要影響[2]。不同于北半球風(fēng)暴路徑,南半球風(fēng)暴路徑的下墊面為廣闊的南大洋,在空間分布、時(shí)間演變等方面有許多區(qū)別于北半球風(fēng)暴路徑的特征。例如,南半球風(fēng)暴路徑的緯向分布更加均勻?qū)ΨQ[3],具有顯著的半年周期循環(huán)信號(hào)[4],夏季可能比冬季更強(qiáng)[5]等等。近年來(lái),南半球的氣候變化及其對(duì)全球變化的響應(yīng)和影響受到越來(lái)越多的關(guān)注[6-7]。在南半球氣候變異因子中已經(jīng)得到廣泛共識(shí)的就是南半球環(huán)狀模[8](Southern Annular Mode, SAM)正處于正位相時(shí)期[9]。SAM的正位相意味著南極繞極渦旋增強(qiáng),同時(shí)繞極西風(fēng)急流也會(huì)增強(qiáng)。事實(shí)上,近年來(lái)南半球西風(fēng)帶的增強(qiáng)和向極偏移已經(jīng)得到越來(lái)越多的證據(jù)支持[10-13]。在這樣的背景下,南半球的風(fēng)暴路徑也出現(xiàn)了對(duì)背景流場(chǎng)變化的響應(yīng)。例如,Simmonds和Keay[14]指出南半球氣旋活動(dòng)越來(lái)越多地生成在相對(duì)較高的緯度。Frederiksen and Frederiksen[15]系統(tǒng)地分析了南半球冬季風(fēng)暴路徑的年代際變化。Luo等[16]根據(jù)理論計(jì)算模型確認(rèn)了SAM與風(fēng)暴路徑之間的聯(lián)系。Yin[17]根據(jù)IPCC AR4多模式的結(jié)果,發(fā)現(xiàn)在未來(lái)情景下風(fēng)暴活動(dòng)出現(xiàn)了持續(xù)的向極偏移。

    以往對(duì)南半球風(fēng)暴路徑強(qiáng)度和位置變化的研究都是較為定性的描述,并且缺乏對(duì)某一層的全面深入分析,所研究的時(shí)間跨度也普遍偏短。本文利用1956—2015年共60年的數(shù)據(jù),在進(jìn)一步確認(rèn)南半球風(fēng)暴路徑具有的基本特征的基礎(chǔ)上,通過(guò)定義一種能夠反映風(fēng)暴活動(dòng)帶變化的指數(shù),定量地分析了南半球300 hPa風(fēng)暴活動(dòng)在過(guò)去60年間的趨勢(shì)變化及其可能原因。

    1 數(shù)據(jù)與方法

    本文主要用到了美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心和美國(guó)國(guó)家大氣研究中心的NCEP/NCAR Reanalysis 1(NCEP1)資料集,包括位勢(shì)高度場(chǎng)、風(fēng)場(chǎng)和溫度場(chǎng)數(shù)據(jù),水平分辨率為2.5°×2.5°,垂直方向按等壓面共分為17層(數(shù)據(jù)下載地址:https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html)。此外,本文還用到了英國(guó)氣象局哈得萊中心的HadISST海表面溫度(Sea Surface Temperature, SST)數(shù)據(jù),水平分辨率為1°×1°(數(shù)據(jù)下載地址:http://www.metoffice.gov.uk/hadobs/hadsst2/data/download.html)。本文所選取的時(shí)間范圍為1956—2015年。根據(jù)所研究問(wèn)題的時(shí)間尺度,以上變量除了位勢(shì)高度場(chǎng)選用日平均數(shù)據(jù)以外,其他變量均選用月平均數(shù)據(jù)。

    為了得到風(fēng)暴路徑數(shù)據(jù),首先對(duì)逐日的位勢(shì)高度場(chǎng)數(shù)據(jù)進(jìn)行了2~7 d 周期的帶通濾波,然后逐月計(jì)算濾波數(shù)據(jù)在該月的方差(均方差),該方差(均方差)值一般就作為該月風(fēng)暴活動(dòng)的強(qiáng)度值。本文使用均方差來(lái)表征風(fēng)暴活動(dòng)的強(qiáng)度,這樣可以使以此定義的風(fēng)暴路徑強(qiáng)度和位勢(shì)高度的單位保持一致,均為位勢(shì)米(gpm)。本文用到的主要分析方法包括趨勢(shì)分析與回歸分析。

    2 南半球風(fēng)暴路徑的變化趨勢(shì)

    在分析南半球風(fēng)暴路徑的變化趨勢(shì)之前,本文首先給出了1956—2015年間氣候態(tài)的南半球風(fēng)暴路徑的三維結(jié)構(gòu)特征,然后從風(fēng)暴路徑的強(qiáng)度和范圍兩個(gè)方面闡述了其過(guò)去60年來(lái)的變化趨勢(shì)。

    2.1 三維結(jié)構(gòu)特征

    從圖1可以清楚地看到南半球風(fēng)暴路徑垂向、緯向和經(jīng)向的三維結(jié)構(gòu)特征。垂向上,300 hPa附近的風(fēng)暴路徑信號(hào)最強(qiáng)(見(jiàn)圖1(a),(b))。經(jīng)向上,風(fēng)暴路徑呈單峰型變化(見(jiàn)圖1(a)),中緯度地區(qū)強(qiáng)度最大,極大值出現(xiàn)在50°S附近,然后向赤道和極地兩個(gè)方向逐漸衰減,赤道地區(qū)最弱。從圖1(a)還可以看出,中緯度風(fēng)暴路徑的強(qiáng)信號(hào)可以一直從對(duì)流層低層延伸到平流層,表明中緯度地區(qū)的風(fēng)暴活動(dòng)是一個(gè)非常深厚的系統(tǒng),最強(qiáng)中心從300 hPa隨高度向極傾斜。根據(jù)圖1(a)中反映的風(fēng)暴路徑在中緯度最強(qiáng)的特征,作者選取40°S~70°S緯度范圍進(jìn)行平均,得到了風(fēng)暴路徑沿緯向的分布??梢钥闯?,風(fēng)暴路徑隨經(jīng)度的變化不如隨緯度的變化那么劇烈,緯向分布相對(duì)比較均勻(見(jiàn)圖1(b))。在300 hPa附近可以看到風(fēng)暴路徑在南半球三個(gè)大洋上存在明顯的活動(dòng)中心,其中南印度洋的風(fēng)暴活動(dòng)中心強(qiáng)度最大,系統(tǒng)更加深厚,并且和南大西洋的風(fēng)暴活動(dòng)中心連接;相比之下,南太平洋的風(fēng)暴活動(dòng)最弱,系統(tǒng)更加淺薄。

    ((a)中每個(gè)緯度上的數(shù)據(jù)為整個(gè)緯圈平均的結(jié)果;(b)中每個(gè)經(jīng)度上的數(shù)據(jù)為40°S~70°S平均的結(jié)果。(a) Data at every latitude is the mean result of the whole circle; (b) It is the mean result of 40°S~70°S.)

    圖1 南半球風(fēng)暴路徑的高度-緯度剖面(a)和高度-經(jīng)度剖面的氣候態(tài)(b)分布圖
    Fig. 1 Climatology of height-latitude profile (a) and height-longitude profile(b) of the Southern Hemisphere storm track

    通過(guò)以上分析,可以明確南半球風(fēng)暴路徑在300 hPa信號(hào)最強(qiáng)的特點(diǎn),為此,本文選取出300 hPa風(fēng)暴路徑數(shù)據(jù),進(jìn)一步分析了南半球風(fēng)暴路徑的空間分布特征。從圖2可以看出,南半球風(fēng)暴路徑呈現(xiàn)繞極的帶狀分布,中緯度地區(qū)的強(qiáng)度最大,軸線大約在50°S附近(見(jiàn)圖2綠色虛線)。和圖1(b)的結(jié)果類似,在300 hPa高度上,也能發(fā)現(xiàn)在南印度洋上存在著風(fēng)暴活動(dòng)最強(qiáng)的中心,這個(gè)中心向西與南大西洋東部的風(fēng)暴活動(dòng)中心連接在一起,向東一直延伸到新西蘭南部。這個(gè)橫貫整個(gè)東半球的強(qiáng)風(fēng)暴活動(dòng)帶,與西風(fēng)急流的位置有非常好的對(duì)應(yīng)關(guān)系(見(jiàn)圖2黑色實(shí)線)。特別是南印度洋和南大西洋上的強(qiáng)風(fēng)暴活動(dòng)中心,幾乎都處于西風(fēng)急流所控制的區(qū)域,且風(fēng)暴路徑最大中心位于西風(fēng)急流軸的向極一側(cè)。相比而言,南太平洋的風(fēng)暴路徑活動(dòng)和西風(fēng)急流沒(méi)有明顯的位置上的對(duì)應(yīng)關(guān)系,整個(gè)中緯度地區(qū)的西風(fēng)急流都比較弱,而最強(qiáng)的急流軸出現(xiàn)在靠近赤道的30°S附近。這種西風(fēng)急流軸與風(fēng)暴活動(dòng)中心位置的不匹配,可能是造成前面提到的南太平洋風(fēng)暴活動(dòng)偏弱的原因之一。關(guān)于西風(fēng)急流與風(fēng)暴路徑之間的關(guān)系,將在下一節(jié)詳細(xì)討論。

    (綠線表示50°S緯線所在的位置;灰色弧線表示30°S和60°S緯線所在的位置;紅線表示風(fēng)暴路徑強(qiáng)度等于50 gpm的等值線,在后文中將用到。Green circle line denotes the latitude line of 50°S; Gray curves denote the latitude of 30°S and 60°S; Respectively; Red line denote the storm track intensity of 50 gmp, which will be used in the following text.)

    圖2 300 hPa風(fēng)暴路徑的氣候態(tài)(填色)和西風(fēng)帶(黑色等值線)分布圖

    Fig. 2 Climatology of storm track (shaded) and westerly (black contour) at 300 hPa

    2.2 強(qiáng)度變化趨勢(shì)

    從上述氣候態(tài)的角度已經(jīng)看出,風(fēng)暴活動(dòng)的空間分布與西風(fēng)急流有很好的對(duì)應(yīng)關(guān)系。因此,在南半球西風(fēng)急流出現(xiàn)增強(qiáng)的趨勢(shì)下,有必要考察風(fēng)暴活動(dòng)是否也出現(xiàn)了某種趨勢(shì)變化。圖3給出了風(fēng)暴路徑強(qiáng)度在4個(gè)季節(jié)里隨時(shí)間和緯度變化的異常值??傮w而言,無(wú)論哪個(gè)季節(jié),風(fēng)暴路徑強(qiáng)度都出現(xiàn)了由負(fù)異常到正異常的年代際轉(zhuǎn)變。其中,2000年以來(lái)的正異常值變得更加顯著,表明過(guò)去60年來(lái)南半球風(fēng)暴活動(dòng)出現(xiàn)了增強(qiáng)的趨勢(shì),并在近10多年變得更加明顯。分季節(jié)來(lái)看,風(fēng)暴活動(dòng)在南半球夏季(DJF)的增強(qiáng)幅度略強(qiáng)于其他三個(gè)季節(jié),并且集中在60°S附近(見(jiàn)圖3(a));秋季(MAM)相對(duì)于DJF,增強(qiáng)的位置集中在60°S以北(見(jiàn)圖3(b));冬季(JJA)相對(duì)于DJF,增強(qiáng)的位置集中在60°S以南(見(jiàn)圖3(c));春季(SON)的增強(qiáng)趨勢(shì)則在30°S以南的區(qū)域都分布得相對(duì)較為均勻(見(jiàn)圖3(d))。四個(gè)季節(jié)組合在一起可以刻畫(huà)出南半球風(fēng)暴活動(dòng)帶一致增強(qiáng)的趨勢(shì),而這種趨勢(shì)不僅僅存在于風(fēng)暴活動(dòng)最活躍的中緯度,在高緯度地區(qū)甚至直到極區(qū),也能發(fā)現(xiàn)風(fēng)暴活動(dòng)增強(qiáng)的信號(hào),且高緯度極區(qū)的風(fēng)暴活動(dòng)信號(hào)早在1990s初甚至更早就開(kāi)始出現(xiàn)并向北傳播。這種信號(hào)反映的是風(fēng)暴活動(dòng)的另一種變化趨勢(shì),即位置的向南偏移,其具體特征將在下文中進(jìn)一步分析。

    2.3 范圍變化趨勢(shì)

    為了定量地描述南半球風(fēng)暴活動(dòng)帶范圍的變化,需要確定一個(gè)客觀標(biāo)準(zhǔn)來(lái)表征其位置的變化。本文選取50 gpm等值線作為確定風(fēng)暴活動(dòng)邊界的參考值,定義了風(fēng)暴路徑邊界指數(shù)(Storm Track Boundary Index, STBI)。選取50 gpm的理由是:首先,如果選取的參考值過(guò)小,則北邊的界線就太靠近熱帶,南邊的界線太靠近極區(qū),失去了考察風(fēng)暴活動(dòng)帶南北移動(dòng)的意義;如果選取的值過(guò)大,得到的風(fēng)暴活動(dòng)范圍太窄,甚至?xí)斐赡承┰路菀蛱窖笊系娘L(fēng)暴活動(dòng)較弱而使得南北兩個(gè)環(huán)狀的等值線無(wú)法閉合,無(wú)法考察風(fēng)暴活動(dòng)帶的整體移動(dòng)。其次,兩條50 gpm等值線的多年平均位置分別位于40°S和65°S附近(見(jiàn)圖2紅色等值線),這一范圍基本包含了整個(gè)風(fēng)暴活動(dòng)帶變化最強(qiáng)的范圍。有趣的是,對(duì)研究南半球氣候具有重要意義的SAM指數(shù)就是根據(jù)40°S和65°S兩條緯度帶上的海平面氣壓來(lái)定義的,說(shuō)明這兩條緯度帶上的物理量對(duì)南半球氣候研究具有重要的參考意義。最后,經(jīng)過(guò)敏感性試驗(yàn),發(fā)現(xiàn)選取45~55 gpm之間的數(shù),對(duì)結(jié)果均沒(méi)有顯著影響,得到的變化特征都是相似的。

    綜上所述,本文以50 gpm等值線確定風(fēng)暴活動(dòng)帶的南北邊界,考察其活動(dòng)范圍的變化特征。需要指出的是,這個(gè)參考值并不唯一,可能因采用的數(shù)據(jù)集、濾波方法或選取的時(shí)間段不同存在差異。

    以下是計(jì)算STBI的具體步驟:

    (1) 首先檢驗(yàn)所有月份經(jīng)向上的風(fēng)暴路徑強(qiáng)度最大值,發(fā)現(xiàn)其最大值均超過(guò)50 gpm??紤]到風(fēng)暴路徑在經(jīng)向上是單峰型的變化,說(shuō)明按照50 gpm的標(biāo)準(zhǔn),一定能在經(jīng)向上找到一南一北兩個(gè)點(diǎn)的數(shù)值等于或最接近50 gmp,稱之為風(fēng)暴路徑的南界和北界;

    (2) 由于NCEP1數(shù)據(jù)經(jīng)向的間隔為比較粗的2.5°,不利于更加精確地判斷風(fēng)暴活動(dòng)帶的邊界??紤]到風(fēng)暴活動(dòng)在經(jīng)向上呈單峰型變化,適合將其做簡(jiǎn)單的線性插值以提高分辨率,因此本文將每條經(jīng)線上的數(shù)據(jù)都線性插值成0.1°的間隔,這樣更便于判斷出50 gpm等值線所在的緯度;

    (3) 逐月判斷每條經(jīng)線上數(shù)值最接近50 gpm的點(diǎn)所對(duì)應(yīng)的2個(gè)緯度值,靠近極地一側(cè)的為風(fēng)暴活動(dòng)帶南界,靠近赤道一側(cè)的為北界;

    (4) 計(jì)算所有經(jīng)線上得到的南界和北界的平均值,得到相應(yīng)的STBI時(shí)間序列。

    在過(guò)去的研究中,雖然南半球風(fēng)暴路徑向極偏移的趨勢(shì)已經(jīng)得到確認(rèn)[18],但從本文的STBI指數(shù)可以看出,南半球風(fēng)暴路徑在過(guò)去60年間并非只是單側(cè)地向極偏移,而南北兩側(cè)都在擴(kuò)張。圖4展示了根據(jù)以上步驟計(jì)算得到的風(fēng)暴路徑南界和北界指數(shù)在1956—2015年期間的距平時(shí)間序列,可以看出,北界存在向赤道偏移的趨勢(shì)(見(jiàn)圖4(a)),南界存在向極地偏移的趨勢(shì)(見(jiàn)圖4(b))。兩者具體的變化趨勢(shì)分別為每10年0.34°和0.84°,表明過(guò)去60年間南半球風(fēng)暴活動(dòng)帶正在變得越來(lái)越寬闊。同時(shí)由于南界向極偏移趨勢(shì)是北界向赤道偏移趨勢(shì)的約2.5倍,所以南半球風(fēng)暴路徑向極擴(kuò)張的速度更快,整體的位置變得更加靠近極地。

    圖 3 300 hPa風(fēng)暴路徑在(a)DJF、(b)MAM、(c)JJA和(d)SON 四個(gè)季節(jié)異常值合成的緯度-時(shí)間剖面圖Fig. 3 Latitude-time profiles of the anomalies of storm track for (a) DJF, (b) MAM, (c) JJA and (d) SON at 300 hPa

    除了表現(xiàn)出趨勢(shì)變化以外,南半球風(fēng)暴路徑的南界和北界在過(guò)去60年間似乎還存在若干次比較明顯的突變。圖5給出的對(duì)南界和北界指數(shù)滑動(dòng)T檢驗(yàn)的結(jié)果證實(shí)了這一觀點(diǎn)。1956年以來(lái),發(fā)生突變最明顯的是南界,一共有三次突變超過(guò)99%顯著性檢驗(yàn)(見(jiàn)圖5(b))。第一次突變發(fā)生在1970s初,表現(xiàn)為南界偏北的幅度出現(xiàn)較為明顯的減弱。到了1980s初期,南界出現(xiàn)第二突變,表現(xiàn)為從偏北的位置移動(dòng)到偏南的位置。隨后,南界的位置逐漸向南偏移,大約在1990s前后由此前較為一致的正異常轉(zhuǎn)為此后較為一致的負(fù)異常。到2000年以后,南界又出現(xiàn)第三次突變,表現(xiàn)出一致偏南的特征。相對(duì)于南界,北界發(fā)生過(guò)兩次突變(見(jiàn)圖5(a))。第一次突變時(shí)間大約在1960s初期,表現(xiàn)為北界的位置由偏南轉(zhuǎn)為偏北。但這次突變之后,北界一直維持著在偏南和偏北之間弱幅振蕩的狀態(tài),直到1990s初期,這種南北振蕩開(kāi)始轉(zhuǎn)變?yōu)槌掷m(xù)的偏北變化,發(fā)生第二次突變。

    需要注意的是,自南界最后一次突變以來(lái),南界的偏南異常值比之前顯著加強(qiáng),這說(shuō)明進(jìn)入21世紀(jì)以來(lái),南半球風(fēng)暴路徑向極偏移的幅度有增加的趨勢(shì)。南偏事件正變得越來(lái)越頻繁,在2012年2月份甚至發(fā)生了一次南界向極異常偏移超過(guò)15個(gè)緯度的極端事件,這勢(shì)必對(duì)南半球中高緯地區(qū)的天氣和氣候變化產(chǎn)生一定影響。例如,已有研究指出,南半球風(fēng)暴路徑的向極偏移,帶來(lái)了南半球中高緯地區(qū)更加頻繁的氣旋 活動(dòng)和降水發(fā)生[2]。因此,近年來(lái)南半球風(fēng)暴路徑增強(qiáng)并向極偏移的趨勢(shì),應(yīng)值得特別關(guān)注。

    ((a)北界指數(shù);(b)南界指數(shù);正值代表向北的異常;負(fù)值代表向南的異常;圖中"Trend"后面的數(shù)字代表STBI的變化趨勢(shì)。(a) Northern STBI; (b) Southern STBI; Positive value means northward anomaly; Negative value means southward anomaly; Values behind the string "Trend" denote the variation trend of STBI.)

    圖4 風(fēng)暴路徑邊界指數(shù)距平值時(shí)間序列
    Fig. 4 Time series of STBI anomalies

    ((a)北界指數(shù);(b)南界指數(shù)的突變檢驗(yàn);紅色和藍(lán)色柱狀圖代表STBI距平值;黑線為T(mén)檢驗(yàn)的檢驗(yàn)值;灰色虛線為99%置信度的臨界值。T-test of (a) Northern STBI; (b) southern STBI; Red and blue bars denote the STBI anomalies; Black curves denote the test values of T-test; Gray dashed lines mean the 99% confidence level.)

    圖5 風(fēng)暴路徑邊界指數(shù)的突變檢驗(yàn)結(jié)果
    Fig. 5 Results of the T-test of STBI

    3 南半球風(fēng)暴路徑趨勢(shì)變化的可能原因

    鑒于南半球風(fēng)暴活動(dòng)是具有顯著三維結(jié)構(gòu)的斜壓系統(tǒng),本文從高層的背景西風(fēng)急流變化和低層的SST變化兩個(gè)方面,分析了南半球風(fēng)暴路徑出現(xiàn)的趨勢(shì)變化的可能原因。

    3.1 西風(fēng)急流的可能影響

    由圖2已經(jīng)知道,南半球風(fēng)暴活動(dòng)中心和西風(fēng)急流中心的平均位置十分吻合,表明兩者之間可能具有比較密切的聯(lián)系。事實(shí)上,強(qiáng)盛的繞極西風(fēng)急流是南半球大氣運(yùn)動(dòng)的顯著特征之一,在這樣的背景下天氣尺度瞬變渦旋活動(dòng)也更趨于活躍。因此,南半球風(fēng)暴路徑的強(qiáng)度和位置變化可能就是對(duì)作為背景流的西風(fēng)急流變化的一種響應(yīng)。

    圖6(a)給出的是1956—2015年期間西風(fēng)急流和風(fēng)暴路徑的變化趨勢(shì)圖,可以看出風(fēng)暴路徑和西風(fēng)急流趨勢(shì)變化明顯的區(qū)域具有很好的符合關(guān)系。對(duì)于風(fēng)暴路徑而言,增強(qiáng)趨勢(shì)幾乎都集中在其平均軸線位置50°S以南的區(qū)域,三個(gè)大洋上的增強(qiáng)都比較明顯,最顯著的區(qū)域是南印度洋和南大西洋。風(fēng)暴活動(dòng)帶增強(qiáng)的區(qū)域,幾乎全部包含在西風(fēng)急流增強(qiáng)的區(qū)域(見(jiàn)圖6(a)中50°S附近的0等值線)。另一方面,南印度洋和南太平洋風(fēng)暴路徑在靠近赤道一側(cè)都出現(xiàn)了一定程度的減弱,與之對(duì)應(yīng)的則是同位置的西風(fēng)急流也出現(xiàn)了減弱。由此可見(jiàn),南半球風(fēng)暴路徑無(wú)論增強(qiáng)還是減弱,都與同位置西風(fēng)急流的變化緊密聯(lián)系。從本文的分析中可以定量地看到,風(fēng)暴路徑軸線向極一側(cè)即50°S~65°S之間,西風(fēng)急流的增強(qiáng)趨勢(shì)超過(guò)每10年0.5 m·s-1,且?guī)缀趵@極一周,最大增強(qiáng)趨勢(shì)出現(xiàn)在南印度洋,超過(guò)每10年1 m·s-1。風(fēng)暴路徑增強(qiáng)最顯著的區(qū)域和西風(fēng)急流是幾乎一致的,最大增強(qiáng)趨勢(shì)出現(xiàn)在南印度洋和南大西洋,超過(guò)每10年3 gmp。

    ((a)填色為風(fēng)暴路徑,藍(lán)色等值線為西風(fēng)(單位:m·s-1·decade-1);(b)填色為SST,(a)和(b)中黑線表示風(fēng)暴路徑南界和北界平均位置;綠線表示風(fēng)暴路徑軸線位置。(a) Storm trcak (shaded) and westerlies (bule contour, unit: m·s-1·decade-1); (b) SST (shaded). Black lines in (a) and (b) denote the mean position of northern and southern boundaries; Green line denotes the axis of storm track.)

    圖6 1956—2015年期間風(fēng)暴路徑、西風(fēng)和SST的變化趨勢(shì)
    Fig. 6 Trends of storm track, westerlies and SST in the period of 1956—2015

    為了分析風(fēng)暴路徑位置的變化與西風(fēng)急流變化之間的聯(lián)系,本文將南北界指數(shù)距平值回歸到緯向風(fēng)距平場(chǎng)。從圖7(a)和圖7(c)中可以看出,風(fēng)暴路徑南界和北界的移動(dòng)顯著地受到西風(fēng)變化的影響。由于北界指數(shù)的變化趨勢(shì)相對(duì)較弱,而南界指數(shù)和西風(fēng)的變化趨勢(shì)較強(qiáng),因此北界指數(shù)對(duì)西風(fēng)的回歸系數(shù)普遍要低于南界。就回歸場(chǎng)總體的分布而言,風(fēng)暴路徑對(duì)其軸線50°S以南和以北西風(fēng)急流的變化十分敏感,尤其是風(fēng)暴路徑南界。結(jié)合西風(fēng)急流的實(shí)際變化(見(jiàn)圖6(a)藍(lán)色等值線),風(fēng)暴路徑軸線以南西風(fēng)的增強(qiáng)趨勢(shì)造成了風(fēng)暴路徑向極地偏移,以北西風(fēng)的減弱趨勢(shì)造成了風(fēng)暴路徑的向赤道偏移,最終表現(xiàn)為風(fēng)暴活動(dòng)帶的南北擴(kuò)張。值得注意的是,西風(fēng)對(duì)風(fēng)暴路徑北界的影響呈現(xiàn)三極型特征(見(jiàn)圖7(a)),尤其在南太平洋區(qū)域表現(xiàn)得更為明顯。當(dāng)40°S~50°S之間的西風(fēng)增強(qiáng)(減弱)、30°S以北和65°S以南的西風(fēng)減弱(增強(qiáng))時(shí),風(fēng)暴路徑的北界出現(xiàn)向赤道(極地)的變化。由此可見(jiàn),風(fēng)暴路徑的北界變化較為復(fù)雜,同時(shí)受副熱帶、中緯度和副極地地區(qū)西風(fēng)急流的影響,這與Nakamura and Shimpo[19]的結(jié)論是一致的。需要特別指出的是,西風(fēng)急流盡管能造成風(fēng)暴路徑邊界南北的移動(dòng),但是對(duì)風(fēng)暴路徑軸線(50°S)的位置幾乎沒(méi)有影響,這說(shuō)明風(fēng)暴路徑軸線的位置相對(duì)穩(wěn)定,不隨風(fēng)暴活動(dòng)帶南北的移動(dòng)而發(fā)生明顯改變。

    (北界指數(shù)回歸(a)西風(fēng)和(b)SST,南界指數(shù)回歸(c)西風(fēng)和(d)SST,所有數(shù)據(jù)均為距平值;(a)~(d)中黑線表示風(fēng)暴路徑南界和北界平均位置;綠線表示風(fēng)暴路徑軸線位置;黑點(diǎn)表示過(guò)99%置信度檢驗(yàn)的區(qū)域。Regression of (a) westerlies and (b) SST with northern STBI, (c) and (d) are the same but for southern STBI. All the data are processed as anomaly values; Black lines in (a)~(d) denote the mean position of northern; Southern boundaries and green line denotes the axis of storm track; Black dots denote areas exceeding the 99% confidence level.)

    圖7 風(fēng)暴路徑邊界指數(shù)與西風(fēng)和SST經(jīng)向梯度的線性回歸圖
    Fig. 7 Linear regression of westerlies and SST with STBI

    概括而言,無(wú)論是從局地還是非局地的角度,無(wú)論是強(qiáng)度還是位置的變化,南半球風(fēng)暴路徑的變化都與西風(fēng)急流的變化緊密聯(lián)系在一起。

    3.2 SST的可能影響

    前面已經(jīng)提到,南半球風(fēng)暴路徑的一個(gè)顯著特點(diǎn)就是下墊面幾乎全部為海洋,海洋的加熱和冷卻對(duì)風(fēng)暴路徑的作用不容忽視。已有研究表明,在南半球開(kāi)闊的大洋上,氣旋活動(dòng)的生成位置與SST經(jīng)向梯度最大的位置密切相關(guān)[14, 20]。Inatsu and Hoskins[21]用數(shù)值試驗(yàn)進(jìn)一步證實(shí),南半球低層的風(fēng)暴路徑受SST變化的影響。下面將分析SST的變化對(duì)300 hPa風(fēng)暴路徑變化的可能影響。

    從圖6(b)中可以看出,三個(gè)大洋上的SST以50°S附近為界,以北的區(qū)域以增暖趨勢(shì)為主,以南的區(qū)域以變冷趨勢(shì)為主。對(duì)比圖6(a)和6(b)可以看出,風(fēng)暴路徑和SST的變化趨勢(shì)有較好的空間對(duì)應(yīng)關(guān)系:SST變冷的區(qū)域?qū)?yīng)的是風(fēng)暴路徑的增強(qiáng)趨勢(shì),SST增暖的區(qū)域?qū)?yīng)的是風(fēng)暴路徑變化較小甚至減弱趨勢(shì)。SST北邊暖南邊冷的變化趨勢(shì),增大了SST經(jīng)向梯度,有利于風(fēng)暴路徑這樣深厚的斜壓系統(tǒng)的發(fā)展。將風(fēng)暴路徑南北界指數(shù)距平值回歸到SST距平場(chǎng)。結(jié)合SST的變化趨勢(shì)可以看出,SST增暖的區(qū)域(見(jiàn)圖6(b)正值區(qū)域)幾乎都對(duì)應(yīng)著風(fēng)暴路徑北界的向赤道偏移(見(jiàn)圖7(b)正值區(qū)域)和南界的向極地偏移(見(jiàn)圖7(d)負(fù)值區(qū)域),而在中高緯SST變冷的區(qū)域,風(fēng)暴路徑北界和南界的變化趨勢(shì)正好相反。這表明,SST的增暖有利于風(fēng)暴路徑變得更寬,而SST變冷則會(huì)讓風(fēng)暴路徑出現(xiàn)收縮的趨勢(shì)。由于SST的增暖趨勢(shì)更強(qiáng),范圍更廣,因此風(fēng)暴路徑最終體現(xiàn)的是擴(kuò)張的趨勢(shì)。

    4 結(jié)論與討論

    本文通過(guò)1956—2015年的位勢(shì)高度場(chǎng)2~7 d帶通濾波數(shù)據(jù),首先確認(rèn)了南半球風(fēng)暴路徑的三維結(jié)構(gòu)特征。緯向上,風(fēng)暴路徑呈單峰型變化,中緯度區(qū)域最強(qiáng),最大值大約在50°S附近,并且該軸線能自下而上一直延伸到平流層;經(jīng)向上,風(fēng)暴路徑在三個(gè)大洋上存在強(qiáng)活動(dòng)中心,其中南印度洋最強(qiáng),南太平洋最弱;垂向上,300 hPa附近風(fēng)暴路徑的信號(hào)最強(qiáng),平流層和對(duì)流層低層相對(duì)較弱。

    在此基礎(chǔ)上,本文選取300 hPa的風(fēng)暴路徑數(shù)據(jù),重點(diǎn)分析了該層南半球風(fēng)暴活動(dòng)在過(guò)去60年間的變化趨勢(shì)。結(jié)果表明,南半球風(fēng)暴路徑具有顯著的增強(qiáng)趨勢(shì),增強(qiáng)最明顯的區(qū)域位于50°S以南,最大中心在南印度洋和南大西洋上,趨勢(shì)可超過(guò)每10年3 gpm。通過(guò)構(gòu)造風(fēng)暴路徑邊界指數(shù)發(fā)現(xiàn),南半球風(fēng)暴路徑不僅存在向極地的偏移,還存在向赤道的偏移,表明風(fēng)暴活動(dòng)帶的范圍正在變得越來(lái)越寬。但是,南界的偏移速度更快,大約是北界的2.5倍。

    分析表明,南半球風(fēng)暴路徑與背景場(chǎng)西風(fēng)急流和下墊面SST的變化存在密切關(guān)系。中高緯地區(qū)的西風(fēng)急流增強(qiáng),使得斜壓不穩(wěn)定波從背景流獲取到更多的擾動(dòng)能量,風(fēng)暴路徑因此強(qiáng)度變得更大,南界更加靠近極地。強(qiáng)盛的繞極西風(fēng)還加劇了下墊面海水的抽吸和混合,表現(xiàn)為中高緯地區(qū)SST的變冷趨勢(shì)。在中低緯地區(qū),西風(fēng)是減弱的趨勢(shì),缺少能量來(lái)源的風(fēng)暴活動(dòng)變化很小甚至出現(xiàn)了減弱的趨勢(shì),而西風(fēng)減弱也意味著風(fēng)暴路徑會(huì)有一定程度的向赤道偏移。另外,中低緯地區(qū)的SST具有明顯的增暖趨勢(shì),這種加熱作用有利于風(fēng)暴路徑的“受熱膨脹”,使得風(fēng)暴活動(dòng)范圍變寬。與此同時(shí),SST南冷北暖的變化趨勢(shì)也使得SST的經(jīng)向梯度增加,有利于斜壓系統(tǒng)的維持和發(fā)展,對(duì)風(fēng)暴活動(dòng)的發(fā)展也有促進(jìn)作用。

    綜上所述,西風(fēng)急流南強(qiáng)北弱和SST南冷北暖的變化趨勢(shì)是南半球風(fēng)暴路徑增強(qiáng)和南北擴(kuò)張的重要原因,其中以西風(fēng)急流的影響最為顯著。受西風(fēng)急流明顯向極偏移趨勢(shì)的影響,南半球風(fēng)暴路徑也主要表現(xiàn)為向極地的偏移。

    需要特別注意的是,進(jìn)入21世紀(jì)以來(lái),南半球風(fēng)暴活動(dòng)帶的南界向極地偏移的幅度有增大的趨勢(shì),極端南偏的事件變得更為頻繁,可以預(yù)見(jiàn)這將對(duì)南半球中高緯地區(qū)的極端天氣和氣候變化產(chǎn)生較大影響,應(yīng)值得關(guān)注。

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

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

    [3] Trenberth K E. Storm tracks in the Southern Hemisphere [J]. J Atmos Sci, 1991, 48(19): 2159-2178.

    [4] Meehl G A, Hurrell J W, Loon H van. A modulation of the mechanism of the semiannual oscillation in the Southern Hemisphere [J]. Tellus, 1998, 50(4): 442-450.

    [5] Hoskins B J, Hodges K I. A new perspective on Southern Hemisphere storm tracks [J]. J Climate, 2005, 18(20): 4108-4129.

    [6] Kushner P J, Held I M, Delworth T L. Southern Hemisphere atmospheric circulation response to global warming [J]. J Climate, 2001, 14(10): 2238-2249.

    [7] Gillett N P, Thompson D W J. Simulation of recent Southern Hemisphere climate change [J]. Science, 2003, 302(5643): 273-275.

    [8] Limpasuvan V, Hartmann D L. Eddies and the annular modes of climate variability [J]. Geophys Res Lett, 1999, 26(20): 3133-3136.

    [9] Hall A, Visbeck M . Synchronous variability in the Southern Hemisphere atmosphere, sea ice, and ocean resulting from the annular mode [J]. J Climate, 2002, 15(21): 3043-3057.

    [10] Thompson D W J, Solomon S. Interpretation of recent Southern Hemisphere climate change [J]. Science, 2002, 296(5569): 895-899.

    [11] Marshall G J. Trends in the Southern Annular Mode from observations and reanalyses[J]. J Climate, 2003, 16(16): 4134-4143.

    [12] Arblaster J M, Meehl G A. Contributions of external forcings to Southern Annular Mode trends [J]. J Climate, 2006, 19(12): 2896-2905.

    [13] Spence P, Grifies S M, England M H, et al. Rapid subsurface warming and circulation changes of Antarctic coastal waters by poleward shifting winds [J]. Geophys Res Lett, 2014, 41(13): 4601-4610.

    [14] Simmonds I, Keay K. Mean Southern Hemisphere extratropical cyclone behavior in the 40-year NCEP-NCAR reanalysis [J]. J Climate, 2000, 13(5): 873-885.

    [15] Frederiksen J, Frederiksen C. Interdecadal changes in southern hemisphere winter storm track modes [J]. Tellus, 2007, 59(5): 599-617.

    [16] Luo D, Gong T, Diao Y, et al. Storm tracks and annular modes [J]. Geophys Res Lett, 2007, 34(17): L17701.

    [17] Yin J H. A consistent poleward shift of the storm tracks in simulations of 21st century climate [J]. Geophys Res Lett, 2005, 32(18): L18701.

    [18] Son S W, Tandon N F, Polvani L M, et al. Ozong hole and Southern Hemisphere climate change [J]. Geophys Res Lett, 2009, 36(15): L15705.

    [19] Nakamura H, Shimpo A. Seasonal variations in the Southern Hemisphere storm tracks and jet streams in a reanalysis dataset [J]. J Climate, 2004, 17(9): 1828-1844.

    [20] Sinclair M R. A climatology of cyclogenesis for the Southern Hemisphere [J]. Mon Wea Rev, 1995, 123(6): 1601-1619.

    [21] Inatsu M, Hoskins B J. The zonal asymmetry of the Southern Hemisphere winter storm track [J]. J Climate, 2004, 53(17): 4882-4892.

    猜你喜歡
    南半球急流西風(fēng)
    智海急流
    智海急流
    古道西風(fēng)(水彩)
    悉尼中國(guó)城:南半球華人生活縮影
    智海急流
    智海急流
    澳大利亞:南半球的動(dòng)物王國(guó)
    南半球有個(gè)奇異國(guó)
    秋漸濃,分付西風(fēng)此夜涼
    文苑(2019年20期)2019-11-16 08:52:46
    浣溪沙 誰(shuí)念西風(fēng)獨(dú)自涼
    飛魔幻A(2019年3期)2019-06-12 09:22:35
    成人精品一区二区免费| 夜夜看夜夜爽夜夜摸| 999久久久国产精品视频| 男人舔女人下体高潮全视频| 欧美激情高清一区二区三区| 别揉我奶头~嗯~啊~动态视频| www.精华液| 老司机深夜福利视频在线观看| 色哟哟哟哟哟哟| 免费av毛片视频| 国产精品久久久久久人妻精品电影| 日韩一卡2卡3卡4卡2021年| 久久青草综合色| 亚洲欧美精品综合一区二区三区| av在线播放免费不卡| 香蕉丝袜av| x7x7x7水蜜桃| 韩国av一区二区三区四区| 亚洲精品在线美女| 亚洲熟妇中文字幕五十中出| 精品熟女少妇八av免费久了| 麻豆一二三区av精品| 国产成人欧美| 视频区欧美日本亚洲| 99香蕉大伊视频| 亚洲欧美日韩高清在线视频| 12—13女人毛片做爰片一| 日韩欧美一区二区三区在线观看| 日本 欧美在线| 国产亚洲欧美精品永久| 午夜激情av网站| 午夜福利欧美成人| 中亚洲国语对白在线视频| 校园春色视频在线观看| 黄色成人免费大全| 在线免费观看的www视频| 亚洲精品久久国产高清桃花| av天堂久久9| 一二三四在线观看免费中文在| 国产精品香港三级国产av潘金莲| 狠狠狠狠99中文字幕| 中国美女看黄片| а√天堂www在线а√下载| 日韩成人在线观看一区二区三区| 最好的美女福利视频网| 巨乳人妻的诱惑在线观看| 他把我摸到了高潮在线观看| 叶爱在线成人免费视频播放| 亚洲国产欧美一区二区综合| 如日韩欧美国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 欧美中文日本在线观看视频| 一本大道久久a久久精品| 午夜影院日韩av| 国产精品av久久久久免费| 精品国产一区二区久久| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 搡老妇女老女人老熟妇| √禁漫天堂资源中文www| 国产乱人伦免费视频| 欧美在线黄色| 午夜免费激情av| 国产日韩一区二区三区精品不卡| 国产欧美日韩一区二区三区在线| 叶爱在线成人免费视频播放| 日本三级黄在线观看| 老司机在亚洲福利影院| 亚洲成av片中文字幕在线观看| av超薄肉色丝袜交足视频| 亚洲狠狠婷婷综合久久图片| 一本大道久久a久久精品| 身体一侧抽搐| 国产精品,欧美在线| 免费人成视频x8x8入口观看| 日本vs欧美在线观看视频| 国产av一区在线观看免费| 国产激情欧美一区二区| 亚洲午夜精品一区,二区,三区| 波多野结衣高清无吗| 日韩欧美一区视频在线观看| 性色av乱码一区二区三区2| 大香蕉久久成人网| 校园春色视频在线观看| 一本大道久久a久久精品| 亚洲全国av大片| 亚洲第一av免费看| 黄片小视频在线播放| 黄频高清免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人啪精品午夜网站| 国内久久婷婷六月综合欲色啪| 搡老熟女国产l中国老女人| 少妇粗大呻吟视频| 久久人人精品亚洲av| 黄片小视频在线播放| 国产精品日韩av在线免费观看 | 99久久99久久久精品蜜桃| 久久欧美精品欧美久久欧美| 叶爱在线成人免费视频播放| 精品欧美一区二区三区在线| 女同久久另类99精品国产91| 亚洲av日韩精品久久久久久密| av片东京热男人的天堂| 亚洲欧美精品综合一区二区三区| 黄色视频不卡| 亚洲欧美精品综合久久99| 精品福利观看| 国产午夜精品久久久久久| 18禁黄网站禁片午夜丰满| 亚洲午夜理论影院| av在线播放免费不卡| 最近最新免费中文字幕在线| 9色porny在线观看| 亚洲精品一区av在线观看| 男女之事视频高清在线观看| 国产成人av激情在线播放| 少妇粗大呻吟视频| а√天堂www在线а√下载| 国内久久婷婷六月综合欲色啪| 亚洲国产高清在线一区二区三 | 97人妻精品一区二区三区麻豆 | 性色av乱码一区二区三区2| 久久久精品国产亚洲av高清涩受| 国产精品国产高清国产av| 欧美黑人精品巨大| 青草久久国产| 无限看片的www在线观看| av欧美777| 亚洲五月婷婷丁香| 欧美乱码精品一区二区三区| 亚洲第一av免费看| 好男人在线观看高清免费视频 | 亚洲精品在线观看二区| 日本a在线网址| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区视频了| 亚洲av日韩精品久久久久久密| 首页视频小说图片口味搜索| 日本黄色视频三级网站网址| 99国产综合亚洲精品| 自拍欧美九色日韩亚洲蝌蚪91| 九色亚洲精品在线播放| 久久久久久久午夜电影| 国产亚洲欧美精品永久| 国产欧美日韩一区二区三区在线| 久久中文字幕一级| 亚洲av电影不卡..在线观看| 18禁观看日本| av网站免费在线观看视频| 亚洲情色 制服丝袜| 日韩免费av在线播放| 国产精品av久久久久免费| 欧洲精品卡2卡3卡4卡5卡区| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 国产激情久久老熟女| 国产视频一区二区在线看| 女性生殖器流出的白浆| 欧洲精品卡2卡3卡4卡5卡区| 久久人人97超碰香蕉20202| 欧美激情高清一区二区三区| 黄色a级毛片大全视频| av中文乱码字幕在线| 99国产精品一区二区蜜桃av| 国产精品一区二区免费欧美| 久久久国产欧美日韩av| 久9热在线精品视频| 91九色精品人成在线观看| 男女做爰动态图高潮gif福利片 | 亚洲五月色婷婷综合| 国产精品一区二区三区四区久久 | 黄色丝袜av网址大全| 欧美亚洲日本最大视频资源| 最近最新中文字幕大全免费视频| 国产精品野战在线观看| 国产精品久久久人人做人人爽| 一区在线观看完整版| 国产伦一二天堂av在线观看| 岛国视频午夜一区免费看| 欧美色欧美亚洲另类二区 | 在线观看免费午夜福利视频| 国产乱人伦免费视频| 精品乱码久久久久久99久播| 一夜夜www| 久久天躁狠狠躁夜夜2o2o| 日韩 欧美 亚洲 中文字幕| videosex国产| 免费在线观看黄色视频的| 久久人人爽av亚洲精品天堂| 黑人操中国人逼视频| 日本黄色视频三级网站网址| 日本免费a在线| 日韩欧美国产在线观看| 美女免费视频网站| 人人妻人人澡人人看| 成人国产综合亚洲| 国产精品自产拍在线观看55亚洲| 成人18禁在线播放| 欧美日本亚洲视频在线播放| 欧美人与性动交α欧美精品济南到| 亚洲专区字幕在线| 亚洲精品国产精品久久久不卡| 99久久精品国产亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线av高清观看| 精品午夜福利视频在线观看一区| 亚洲专区中文字幕在线| 亚洲五月婷婷丁香| 成年版毛片免费区| 少妇裸体淫交视频免费看高清 | 国产精品av久久久久免费| 亚洲欧洲精品一区二区精品久久久| 亚洲 国产 在线| 日本 欧美在线| 色综合婷婷激情| 国产蜜桃级精品一区二区三区| bbb黄色大片| 一级a爱视频在线免费观看| 亚洲中文av在线| 精品久久久久久成人av| 91字幕亚洲| 久久久久国内视频| 久久精品亚洲熟妇少妇任你| 露出奶头的视频| 在线播放国产精品三级| 在线视频色国产色| 亚洲欧美激情综合另类| 不卡一级毛片| 久久久久国内视频| 一级毛片精品| 99国产极品粉嫩在线观看| 中亚洲国语对白在线视频| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 搡老妇女老女人老熟妇| 精品卡一卡二卡四卡免费| 免费av毛片视频| 亚洲一区高清亚洲精品| 亚洲精品在线美女| 免费看十八禁软件| 亚洲av五月六月丁香网| 18禁美女被吸乳视频| 欧美一级a爱片免费观看看 | 亚洲第一av免费看| aaaaa片日本免费| 女同久久另类99精品国产91| 夜夜爽天天搞| 香蕉久久夜色| 久久久久久国产a免费观看| 丰满人妻熟妇乱又伦精品不卡| 免费看a级黄色片| 久久午夜亚洲精品久久| 曰老女人黄片| 51午夜福利影视在线观看| 亚洲欧美日韩高清在线视频| 欧美午夜高清在线| 中出人妻视频一区二区| 午夜福利18| 香蕉国产在线看| 又黄又爽又免费观看的视频| 日韩欧美免费精品| 丝袜美足系列| 亚洲久久久国产精品| 免费久久久久久久精品成人欧美视频| 97人妻精品一区二区三区麻豆 | 日韩视频一区二区在线观看| 怎么达到女性高潮| 日韩欧美三级三区| 午夜成年电影在线免费观看| 人人妻人人澡欧美一区二区 | 波多野结衣高清无吗| www.www免费av| 在线十欧美十亚洲十日本专区| 91精品三级在线观看| 精品欧美一区二区三区在线| 国产三级在线视频| 极品人妻少妇av视频| 亚洲精品中文字幕在线视频| 最近最新免费中文字幕在线| 精品福利观看| 国产精品久久久久久人妻精品电影| 欧美日本视频| or卡值多少钱| 精品国产乱子伦一区二区三区| 午夜福利一区二区在线看| 欧美中文日本在线观看视频| 亚洲美女黄片视频| 很黄的视频免费| 每晚都被弄得嗷嗷叫到高潮| 后天国语完整版免费观看| 久久国产精品影院| 中国美女看黄片| www日本在线高清视频| 国产一区二区三区在线臀色熟女| 成人欧美大片| 99精品欧美一区二区三区四区| 欧洲精品卡2卡3卡4卡5卡区| avwww免费| 操美女的视频在线观看| 国产av在哪里看| 麻豆国产av国片精品| 久久国产精品影院| 最近最新免费中文字幕在线| 国产精品98久久久久久宅男小说| 男女下面插进去视频免费观看| 韩国av一区二区三区四区| 女生性感内裤真人,穿戴方法视频| 久久热在线av| 在线免费观看的www视频| 多毛熟女@视频| 免费少妇av软件| 国产区一区二久久| 国产野战对白在线观看| 亚洲免费av在线视频| 亚洲 国产 在线| 69精品国产乱码久久久| 12—13女人毛片做爰片一| 91国产中文字幕| 一级黄色大片毛片| 久久精品91无色码中文字幕| 国产精品一区二区精品视频观看| 亚洲国产精品成人综合色| 精品久久蜜臀av无| 天堂影院成人在线观看| 国产私拍福利视频在线观看| 日韩欧美国产在线观看| 他把我摸到了高潮在线观看| 制服人妻中文乱码| 国产高清激情床上av| 一卡2卡三卡四卡精品乱码亚洲| 在线视频色国产色| 91麻豆av在线| 午夜a级毛片| 可以免费在线观看a视频的电影网站| 此物有八面人人有两片| 精品无人区乱码1区二区| 国产精品综合久久久久久久免费 | 一边摸一边抽搐一进一出视频| 久久久久久久午夜电影| 亚洲av片天天在线观看| 亚洲九九香蕉| 天天一区二区日本电影三级 | 老司机在亚洲福利影院| 免费看美女性在线毛片视频| 欧美绝顶高潮抽搐喷水| 亚洲av第一区精品v没综合| 一本大道久久a久久精品| 成在线人永久免费视频| 国产熟女xx| 亚洲性夜色夜夜综合| 真人做人爱边吃奶动态| 日韩视频一区二区在线观看| 久久精品亚洲精品国产色婷小说| 久久精品国产亚洲av高清一级| 中文字幕精品免费在线观看视频| 最近最新中文字幕大全免费视频| 亚洲色图av天堂| 高清黄色对白视频在线免费看| 欧美 亚洲 国产 日韩一| 俄罗斯特黄特色一大片| 国产精品乱码一区二三区的特点 | 国产精品99久久99久久久不卡| 国产伦人伦偷精品视频| 欧美激情久久久久久爽电影 | 动漫黄色视频在线观看| 久久久国产成人免费| 成人av一区二区三区在线看| 无人区码免费观看不卡| 精品人妻1区二区| 99久久国产精品久久久| 91成人精品电影| 女人被狂操c到高潮| 欧美av亚洲av综合av国产av| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美视频一区| 熟女少妇亚洲综合色aaa.| 搡老岳熟女国产| 久久久久国产一级毛片高清牌| 欧美在线黄色| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱 | 亚洲熟妇熟女久久| 一级毛片精品| 极品教师在线免费播放| 久久久久国产精品人妻aⅴ院| 亚洲aⅴ乱码一区二区在线播放 | 欧美乱色亚洲激情| 亚洲精品中文字幕一二三四区| 色老头精品视频在线观看| 免费人成视频x8x8入口观看| 久热爱精品视频在线9| 1024视频免费在线观看| 妹子高潮喷水视频| 久99久视频精品免费| 自线自在国产av| 午夜老司机福利片| 18禁观看日本| 青草久久国产| 一本大道久久a久久精品| 男女做爰动态图高潮gif福利片 | www日本在线高清视频| 少妇裸体淫交视频免费看高清 | 亚洲欧美精品综合久久99| 成人三级黄色视频| 欧美日韩亚洲国产一区二区在线观看| 美女午夜性视频免费| 校园春色视频在线观看| 91麻豆精品激情在线观看国产| 九色国产91popny在线| 亚洲成av片中文字幕在线观看| 国产99久久九九免费精品| 琪琪午夜伦伦电影理论片6080| 亚洲午夜理论影院| 满18在线观看网站| 在线观看66精品国产| 青草久久国产| 自线自在国产av| 国产亚洲欧美98| 亚洲国产精品999在线| 熟女少妇亚洲综合色aaa.| 少妇被粗大的猛进出69影院| 日韩大尺度精品在线看网址 | 在线十欧美十亚洲十日本专区| cao死你这个sao货| 久久精品91蜜桃| 午夜精品久久久久久毛片777| 伊人久久大香线蕉亚洲五| 国产亚洲欧美98| 中文字幕人成人乱码亚洲影| 国产日韩一区二区三区精品不卡| 成在线人永久免费视频| 欧美色欧美亚洲另类二区 | 久久午夜综合久久蜜桃| 99在线人妻在线中文字幕| 一级毛片女人18水好多| 亚洲av电影在线进入| 欧美性长视频在线观看| 久久亚洲真实| 999久久久精品免费观看国产| 18禁裸乳无遮挡免费网站照片 | 欧美+亚洲+日韩+国产| 久久久久久亚洲精品国产蜜桃av| 在线观看免费视频网站a站| 久久国产精品人妻蜜桃| 淫秽高清视频在线观看| 国产精品亚洲美女久久久| 色播亚洲综合网| 一边摸一边做爽爽视频免费| 女警被强在线播放| 日本欧美视频一区| 精品第一国产精品| or卡值多少钱| 国产乱人伦免费视频| 午夜福利视频1000在线观看 | 校园春色视频在线观看| 88av欧美| 午夜日韩欧美国产| 国产高清有码在线观看视频 | 9色porny在线观看| 亚洲狠狠婷婷综合久久图片| 9热在线视频观看99| a在线观看视频网站| 12—13女人毛片做爰片一| 亚洲成a人片在线一区二区| 免费少妇av软件| 欧美大码av| 国产精品日韩av在线免费观看 | 成人手机av| 日韩免费av在线播放| 久久精品国产99精品国产亚洲性色 | 午夜福利欧美成人| 黄色a级毛片大全视频| 嫩草影院精品99| 老熟妇乱子伦视频在线观看| 99在线视频只有这里精品首页| 最好的美女福利视频网| 国产精品永久免费网站| 亚洲精品国产精品久久久不卡| 熟妇人妻久久中文字幕3abv| 日韩高清综合在线| 亚洲av成人不卡在线观看播放网| www.精华液| 一二三四社区在线视频社区8| 午夜日韩欧美国产| av超薄肉色丝袜交足视频| 悠悠久久av| 欧美中文日本在线观看视频| 国产精品,欧美在线| 丝袜美足系列| 中亚洲国语对白在线视频| 女人被躁到高潮嗷嗷叫费观| 精品久久蜜臀av无| 99精品欧美一区二区三区四区| netflix在线观看网站| 亚洲国产欧美一区二区综合| 午夜精品国产一区二区电影| 成人亚洲精品av一区二区| 深夜精品福利| 一本综合久久免费| 亚洲五月婷婷丁香| 欧美日本亚洲视频在线播放| 欧美大码av| 亚洲av电影在线进入| 国产精品99久久99久久久不卡| 99精品久久久久人妻精品| 97人妻精品一区二区三区麻豆 | 咕卡用的链子| 1024视频免费在线观看| 日韩欧美国产一区二区入口| 国产激情久久老熟女| 变态另类丝袜制服| 欧美激情久久久久久爽电影 | 国产精品久久久av美女十八| 国产成人av激情在线播放| 欧美激情高清一区二区三区| 动漫黄色视频在线观看| 夜夜夜夜夜久久久久| 午夜免费观看网址| 久久久久久久久中文| 亚洲国产欧美一区二区综合| 欧美+亚洲+日韩+国产| 欧美中文综合在线视频| 成人免费观看视频高清| 国产精品一区二区精品视频观看| 三级毛片av免费| 男女做爰动态图高潮gif福利片 | 长腿黑丝高跟| 日韩国内少妇激情av| 午夜精品在线福利| 亚洲第一av免费看| 日日夜夜操网爽| 日韩av在线大香蕉| 国产精品综合久久久久久久免费 | 国产三级在线视频| 精品第一国产精品| 性色av乱码一区二区三区2| 999久久久精品免费观看国产| 性欧美人与动物交配| 欧美黄色片欧美黄色片| 亚洲自偷自拍图片 自拍| 两性夫妻黄色片| 国产在线观看jvid| 国产精品一区二区在线不卡| 久久人妻福利社区极品人妻图片| 亚洲 欧美 日韩 在线 免费| 亚洲一区高清亚洲精品| 99国产精品免费福利视频| 亚洲最大成人中文| 99在线视频只有这里精品首页| 精品国产一区二区三区四区第35| 免费在线观看日本一区| 人成视频在线观看免费观看| 自线自在国产av| 正在播放国产对白刺激| 丝袜美腿诱惑在线| 国产一区二区三区视频了| 精品午夜福利视频在线观看一区| 19禁男女啪啪无遮挡网站| 午夜激情av网站| x7x7x7水蜜桃| 午夜福利,免费看| 一本大道久久a久久精品| 精品欧美国产一区二区三| 97人妻精品一区二区三区麻豆 | 少妇熟女aⅴ在线视频| 如日韩欧美国产精品一区二区三区| 久久国产精品男人的天堂亚洲| 成人欧美大片| 18美女黄网站色大片免费观看| 国产男靠女视频免费网站| 精品第一国产精品| 国产高清有码在线观看视频 | 国产亚洲精品av在线| 桃红色精品国产亚洲av| 亚洲第一av免费看| 国产1区2区3区精品| 亚洲精品久久国产高清桃花| 无遮挡黄片免费观看| 亚洲 欧美一区二区三区| 日本 欧美在线| 国产av在哪里看| 精品一区二区三区av网在线观看| 51午夜福利影视在线观看| 亚洲国产欧美日韩在线播放| 国产精品 欧美亚洲| 叶爱在线成人免费视频播放| 妹子高潮喷水视频| 欧美另类亚洲清纯唯美| 免费高清在线观看日韩| 999久久久精品免费观看国产| 精品久久久久久成人av| 嫩草影院精品99| 亚洲精品国产区一区二| 亚洲国产精品久久男人天堂| av在线天堂中文字幕| 欧美亚洲日本最大视频资源| 亚洲国产精品久久男人天堂| 丝袜美腿诱惑在线| 女人被躁到高潮嗷嗷叫费观| 亚洲第一电影网av| 丝袜美腿诱惑在线| 久久久久久人人人人人| 国产精品永久免费网站| 午夜亚洲福利在线播放| 国产一区二区三区视频了| 色综合欧美亚洲国产小说| 啦啦啦 在线观看视频| 在线观看免费视频日本深夜| 91麻豆av在线| 久久久精品欧美日韩精品| 国产一区二区三区视频了| 免费久久久久久久精品成人欧美视频| 大陆偷拍与自拍| 精品国产超薄肉色丝袜足j|