• <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
    观看av在线不卡| 欧美精品人与动牲交sv欧美| 18禁国产床啪视频网站| 老熟女久久久| 女人精品久久久久毛片| av视频免费观看在线观看| 国产精品嫩草影院av在线观看| 黄片无遮挡物在线观看| 在线观看一区二区三区激情| 少妇猛男粗大的猛烈进出视频| 欧美性感艳星| 久久人人爽人人爽人人片va| 国产国拍精品亚洲av在线观看| 一级片'在线观看视频| 亚洲欧洲国产日韩| 国产色婷婷99| 晚上一个人看的免费电影| 日本猛色少妇xxxxx猛交久久| 一级毛片 在线播放| 99久久中文字幕三级久久日本| 欧美精品一区二区免费开放| 婷婷色综合大香蕉| 香蕉国产在线看| 国产永久视频网站| 97超碰精品成人国产| 97在线人人人人妻| 观看美女的网站| 免费大片18禁| 久久ye,这里只有精品| 日本色播在线视频| av在线观看视频网站免费| 国产有黄有色有爽视频| 久久狼人影院| a级毛片黄视频| 中文乱码字字幕精品一区二区三区| 亚洲精品日本国产第一区| 久久精品国产综合久久久 | 91精品伊人久久大香线蕉| 校园人妻丝袜中文字幕| 看免费av毛片| 啦啦啦中文免费视频观看日本| 极品人妻少妇av视频| av天堂久久9| 丁香六月天网| 精品一区二区三卡| 在线观看免费视频网站a站| 国产又色又爽无遮挡免| 久久毛片免费看一区二区三区| 国产日韩欧美在线精品| 老女人水多毛片| 涩涩av久久男人的天堂| 午夜影院在线不卡| 韩国高清视频一区二区三区| 麻豆乱淫一区二区| 国产欧美日韩一区二区三区在线| 久热久热在线精品观看| 成年av动漫网址| 侵犯人妻中文字幕一二三四区| 日韩 亚洲 欧美在线| 婷婷成人精品国产| 国产在线一区二区三区精| 国产欧美日韩一区二区三区在线| 亚洲精品456在线播放app| 女人被躁到高潮嗷嗷叫费观| 在线免费观看不下载黄p国产| 国产亚洲一区二区精品| 999精品在线视频| 久久久久久人妻| 久久久久久久大尺度免费视频| 久久精品国产亚洲av天美| 男女午夜视频在线观看 | freevideosex欧美| 男的添女的下面高潮视频| 男女下面插进去视频免费观看 | 一级毛片电影观看| 国产高清国产精品国产三级| 国产成人91sexporn| 亚洲伊人久久精品综合| 国产男人的电影天堂91| 男女边吃奶边做爰视频| 色哟哟·www| 精品第一国产精品| 国产日韩一区二区三区精品不卡| 中文字幕最新亚洲高清| 91精品国产国语对白视频| 少妇猛男粗大的猛烈进出视频| av卡一久久| 我要看黄色一级片免费的| 亚洲欧美一区二区三区国产| 免费日韩欧美在线观看| 亚洲精品第二区| 超碰97精品在线观看| 国产欧美日韩一区二区三区在线| 亚洲国产成人一精品久久久| 丰满饥渴人妻一区二区三| 飞空精品影院首页| 丁香六月天网| 国产成人一区二区在线| 成人亚洲精品一区在线观看| 九九在线视频观看精品| 狠狠婷婷综合久久久久久88av| 成人亚洲欧美一区二区av| 夫妻性生交免费视频一级片| 日韩av不卡免费在线播放| 成人毛片60女人毛片免费| 深夜精品福利| 日本黄色日本黄色录像| 久久婷婷青草| www日本在线高清视频| 波野结衣二区三区在线| 精品福利永久在线观看| 成人黄色视频免费在线看| 久久久久人妻精品一区果冻| 在现免费观看毛片| 自线自在国产av| 久久久久久久久久久免费av| 亚洲av中文av极速乱| 交换朋友夫妻互换小说| 交换朋友夫妻互换小说| 69精品国产乱码久久久| 精品久久蜜臀av无| 亚洲欧美日韩卡通动漫| 哪个播放器可以免费观看大片| 热99国产精品久久久久久7| av国产精品久久久久影院| 一本久久精品| 国产精品久久久久久久久免| 在线观看人妻少妇| 色5月婷婷丁香| 中文字幕另类日韩欧美亚洲嫩草| 成人免费观看视频高清| 日韩熟女老妇一区二区性免费视频| 女人久久www免费人成看片| 肉色欧美久久久久久久蜜桃| 日韩av不卡免费在线播放| 欧美人与性动交α欧美精品济南到 | 美女主播在线视频| 亚洲精品色激情综合| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产成人一精品久久久| 男女午夜视频在线观看 | 欧美精品国产亚洲| 久久99热这里只频精品6学生| 国内精品宾馆在线| 人妻系列 视频| 精品一区二区三卡| 日韩,欧美,国产一区二区三区| 永久免费av网站大全| 91aial.com中文字幕在线观看| 亚洲国产欧美在线一区| 午夜免费鲁丝| 婷婷色综合www| 久久人人爽人人片av| 日本91视频免费播放| av视频免费观看在线观看| 亚洲成色77777| 免费高清在线观看视频在线观看| 国产一区亚洲一区在线观看| 观看美女的网站| 久久久久精品性色| 国产乱人偷精品视频| 国产不卡av网站在线观看| 高清黄色对白视频在线免费看| 日本午夜av视频| www.熟女人妻精品国产 | 免费观看av网站的网址| 街头女战士在线观看网站| 少妇的逼好多水| 人妻人人澡人人爽人人| 欧美 日韩 精品 国产| 赤兔流量卡办理| 日本免费在线观看一区| 国产精品偷伦视频观看了| 男女国产视频网站| 男人爽女人下面视频在线观看| 涩涩av久久男人的天堂| 亚洲av电影在线进入| 男女啪啪激烈高潮av片| 草草在线视频免费看| 制服人妻中文乱码| 成年人午夜在线观看视频| 高清在线视频一区二区三区| 欧美日本中文国产一区发布| 欧美亚洲 丝袜 人妻 在线| 啦啦啦在线观看免费高清www| 国产精品久久久久久av不卡| 老司机影院成人| 天天影视国产精品| 亚洲经典国产精华液单| av免费观看日本| 亚洲av电影在线观看一区二区三区| 老司机影院成人| 国产色爽女视频免费观看| 99视频精品全部免费 在线| 视频区图区小说| 午夜免费鲁丝| 最近2019中文字幕mv第一页| 女人被躁到高潮嗷嗷叫费观| 在线观看免费视频网站a站| 男女免费视频国产| 亚洲国产精品专区欧美| 亚洲成人手机| 亚洲av综合色区一区| 亚洲精华国产精华液的使用体验| 男的添女的下面高潮视频| 日本wwww免费看| 久久 成人 亚洲| 久久久久久人妻| 看免费成人av毛片| 亚洲天堂av无毛| 爱豆传媒免费全集在线观看| 视频中文字幕在线观看| 免费黄色在线免费观看| 精品人妻在线不人妻| 亚洲欧美精品自产自拍| 精品久久久久久电影网| 免费大片18禁| 国产日韩欧美视频二区| 国产不卡av网站在线观看| 成年美女黄网站色视频大全免费| 亚洲欧美成人综合另类久久久| 国产成人av激情在线播放| 人妻少妇偷人精品九色| 新久久久久国产一级毛片| av有码第一页| 精品熟女少妇av免费看| 在线看a的网站| 久久精品国产自在天天线| 日本色播在线视频| 日日爽夜夜爽网站| 成人二区视频| 欧美人与善性xxx| 男女啪啪激烈高潮av片| 免费久久久久久久精品成人欧美视频 | 色哟哟·www| 在线 av 中文字幕| 极品少妇高潮喷水抽搐| 在线精品无人区一区二区三| 男女边摸边吃奶| 色哟哟·www| 国产 精品1| 欧美精品高潮呻吟av久久| 久久精品熟女亚洲av麻豆精品| 国产av码专区亚洲av| 国产精品秋霞免费鲁丝片| 国产成人a∨麻豆精品| 亚洲精品中文字幕在线视频| 免费看不卡的av| 18在线观看网站| 国产精品人妻久久久影院| 一区二区三区乱码不卡18| 老熟女久久久| 亚洲成人av在线免费| 久久久久久伊人网av| 国产在视频线精品| 国产亚洲av片在线观看秒播厂| 熟女电影av网| 尾随美女入室| 国产精品 国内视频| 欧美另类一区| 亚洲精品国产av成人精品| 国产一区有黄有色的免费视频| videosex国产| 七月丁香在线播放| 2021少妇久久久久久久久久久| 国产色爽女视频免费观看| 精品国产一区二区三区久久久樱花| 国产成人91sexporn| 精品久久久精品久久久| 纯流量卡能插随身wifi吗| 国产探花极品一区二区| 22中文网久久字幕| 久久久国产一区二区| 国产黄色视频一区二区在线观看| 成年动漫av网址| 国产女主播在线喷水免费视频网站| 国产亚洲欧美精品永久| 亚洲精品国产av成人精品| 国产精品久久久av美女十八| 日本黄大片高清| 一本—道久久a久久精品蜜桃钙片| 丰满少妇做爰视频| 少妇精品久久久久久久| 久久99一区二区三区| 亚洲内射少妇av| 亚洲欧美日韩另类电影网站| 国产av一区二区精品久久| 亚洲精品乱码久久久久久按摩| 国产精品三级大全| 秋霞在线观看毛片| 宅男免费午夜| 婷婷色综合大香蕉| 亚洲国产最新在线播放| 亚洲国产毛片av蜜桃av| 亚洲av欧美aⅴ国产| 久久久久人妻精品一区果冻| 欧美亚洲日本最大视频资源| 国产成人精品在线电影| 成人国产av品久久久| 国产综合精华液| 少妇熟女欧美另类| 午夜福利视频在线观看免费| 日本黄大片高清| 欧美国产精品一级二级三级| 超碰97精品在线观看| 人人妻人人澡人人爽人人夜夜| 国产成人免费观看mmmm| 免费日韩欧美在线观看| 高清视频免费观看一区二区| 91aial.com中文字幕在线观看| 卡戴珊不雅视频在线播放| 国产色爽女视频免费观看| 久久99蜜桃精品久久| 亚洲精品av麻豆狂野| 九九爱精品视频在线观看| 日本黄大片高清| 亚洲国产精品一区三区| 国产一级毛片在线| 国产成人精品在线电影| 成人国产av品久久久| 一个人免费看片子| 91午夜精品亚洲一区二区三区| 男女边吃奶边做爰视频| 欧美丝袜亚洲另类| 亚洲,欧美精品.| 日韩精品有码人妻一区| 我要看黄色一级片免费的| 午夜免费男女啪啪视频观看| 国产极品天堂在线| 黄色毛片三级朝国网站| 精品人妻一区二区三区麻豆| 大话2 男鬼变身卡| 欧美国产精品一级二级三级| 国产亚洲一区二区精品| 久久女婷五月综合色啪小说| 中国三级夫妇交换| 久久久久国产精品人妻一区二区| 蜜桃在线观看..| 久久精品国产自在天天线| 春色校园在线视频观看| 国产日韩欧美在线精品| 国产永久视频网站| 久久青草综合色| 亚洲少妇的诱惑av| 亚洲激情五月婷婷啪啪| 免费观看无遮挡的男女| 国产日韩欧美亚洲二区| 纵有疾风起免费观看全集完整版| 精品少妇久久久久久888优播| 精品亚洲成a人片在线观看| 男女啪啪激烈高潮av片| 91在线精品国自产拍蜜月| 一个人免费看片子| 汤姆久久久久久久影院中文字幕| 亚洲 欧美一区二区三区| 国产精品久久久av美女十八| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 精品亚洲成a人片在线观看| 国产精品国产三级专区第一集| 91精品国产国语对白视频| 夜夜爽夜夜爽视频| 国产永久视频网站| 亚洲 欧美一区二区三区| 熟妇人妻不卡中文字幕| 欧美精品一区二区大全| 久久人人爽av亚洲精品天堂| 极品少妇高潮喷水抽搐| 满18在线观看网站| 国产av国产精品国产| 国产免费又黄又爽又色| 国产精品一国产av| 欧美国产精品va在线观看不卡| 国产成人一区二区在线| 亚洲内射少妇av| 男女午夜视频在线观看 | 少妇熟女欧美另类| 精品熟女少妇av免费看| 在线观看美女被高潮喷水网站| 久久久久久久国产电影| 国产永久视频网站| 亚洲成人手机| 亚洲丝袜综合中文字幕| 日本午夜av视频| 99久国产av精品国产电影| 蜜桃在线观看..| 国产成人91sexporn| 亚洲精品成人av观看孕妇| 国产精品一区二区在线观看99| 黑人猛操日本美女一级片| 老熟女久久久| 超色免费av| 国产精品久久久久成人av| 中文字幕制服av| 男女下面插进去视频免费观看 | 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 男人爽女人下面视频在线观看| 国产一区二区激情短视频 | 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 自线自在国产av| 满18在线观看网站| 天堂8中文在线网| 亚洲国产毛片av蜜桃av| 香蕉国产在线看| 99久国产av精品国产电影| 久久久久精品人妻al黑| 18禁裸乳无遮挡动漫免费视频| 另类精品久久| 免费观看av网站的网址| 久久综合国产亚洲精品| 国产爽快片一区二区三区| 边亲边吃奶的免费视频| 欧美另类一区| 新久久久久国产一级毛片| 日本与韩国留学比较| 日本欧美国产在线视频| 国产成人一区二区在线| 在线观看三级黄色| 亚洲精品中文字幕在线视频| 国产精品久久久久久久电影| 天堂8中文在线网| 一级爰片在线观看| 美女主播在线视频| 国产精品99久久99久久久不卡 | 一级毛片我不卡| 久久久久国产网址| 精品久久国产蜜桃| videos熟女内射| av播播在线观看一区| 我要看黄色一级片免费的| 捣出白浆h1v1| 97在线视频观看| 精品亚洲乱码少妇综合久久| 91aial.com中文字幕在线观看| 国产精品嫩草影院av在线观看| 亚洲中文av在线| 97在线视频观看| 一级爰片在线观看| 少妇人妻 视频| 韩国av在线不卡| 寂寞人妻少妇视频99o| 亚洲精品久久午夜乱码| 考比视频在线观看| 七月丁香在线播放| 亚洲精品国产av蜜桃| 国产极品天堂在线| 亚洲国产色片| 久久久国产精品麻豆| 国产高清国产精品国产三级| 国产色爽女视频免费观看| 五月伊人婷婷丁香| 国产精品久久久久久精品电影小说| 大香蕉久久网| 在线观看国产h片| 亚洲三级黄色毛片| 亚洲综合色惰| 成年av动漫网址| 亚洲天堂av无毛| 秋霞在线观看毛片| 交换朋友夫妻互换小说| 性色avwww在线观看| 国产在线免费精品| 久久综合国产亚洲精品| 精品一品国产午夜福利视频| 草草在线视频免费看| 国产成人欧美| 成人毛片a级毛片在线播放| 亚洲欧美清纯卡通| 国产白丝娇喘喷水9色精品| 免费大片18禁| 赤兔流量卡办理| 高清av免费在线| 99re6热这里在线精品视频| 国产男女超爽视频在线观看| 国产综合精华液| 精品视频人人做人人爽| 最黄视频免费看| 久久精品夜色国产| 中文字幕人妻丝袜制服| 日日摸夜夜添夜夜爱| 少妇被粗大的猛进出69影院 | 精品久久久精品久久久| 国产探花极品一区二区| 又粗又硬又长又爽又黄的视频| 男女午夜视频在线观看 | 狠狠精品人妻久久久久久综合| 欧美日韩成人在线一区二区| 成年人午夜在线观看视频| 欧美日韩一区二区视频在线观看视频在线| 91精品国产国语对白视频| 捣出白浆h1v1| 在线观看国产h片| 精品人妻在线不人妻| 精品久久蜜臀av无| 国产高清三级在线| 日韩制服骚丝袜av| 国产深夜福利视频在线观看| 91国产中文字幕| a 毛片基地| 亚洲一码二码三码区别大吗| 久久午夜福利片| 日本黄色日本黄色录像| 丰满迷人的少妇在线观看| 亚洲,欧美精品.| 久久久国产一区二区| 久久精品国产a三级三级三级| 中文字幕人妻熟女乱码| 欧美成人午夜精品| 亚洲精品久久久久久婷婷小说| 下体分泌物呈黄色| 成年美女黄网站色视频大全免费| av国产久精品久网站免费入址| 亚洲在久久综合| 欧美性感艳星| a级片在线免费高清观看视频| 国国产精品蜜臀av免费| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡动漫免费视频| 久久久久久人妻| 亚洲国产av影院在线观看| 自线自在国产av| 亚洲欧美成人综合另类久久久| 欧美 日韩 精品 国产| 久久99精品国语久久久| 秋霞在线观看毛片| av在线观看视频网站免费| a级毛色黄片| 2022亚洲国产成人精品| 高清欧美精品videossex| 色网站视频免费| 人妻人人澡人人爽人人| 国产在线视频一区二区| 人人妻人人添人人爽欧美一区卜| kizo精华| 美女xxoo啪啪120秒动态图| 热99久久久久精品小说推荐| 成人影院久久| 宅男免费午夜| 狂野欧美激情性xxxx在线观看| 蜜桃在线观看..| 国产欧美另类精品又又久久亚洲欧美| 免费看av在线观看网站| 狠狠婷婷综合久久久久久88av| 亚洲情色 制服丝袜| 伊人久久国产一区二区| 成人毛片a级毛片在线播放| 一级毛片电影观看| 丝袜在线中文字幕| 亚洲av电影在线进入| 99久久中文字幕三级久久日本| 亚洲人与动物交配视频| 欧美激情极品国产一区二区三区 | 乱码一卡2卡4卡精品| 日本wwww免费看| 香蕉精品网在线| 成人黄色视频免费在线看| 99热这里只有是精品在线观看| av播播在线观看一区| 丝袜在线中文字幕| 欧美bdsm另类| 久久精品国产亚洲av天美| av免费在线看不卡| 天天影视国产精品| 亚洲内射少妇av| 亚洲精华国产精华液的使用体验| 伦理电影免费视频| a级片在线免费高清观看视频| 国产高清国产精品国产三级| 国产毛片在线视频| 欧美性感艳星| 免费在线观看完整版高清| 午夜免费鲁丝| 国产一区二区三区综合在线观看 | 99re6热这里在线精品视频| 国产1区2区3区精品| 人体艺术视频欧美日本| 女性被躁到高潮视频| 男女边摸边吃奶| 欧美成人午夜精品| 成人黄色视频免费在线看| 男人舔女人的私密视频| 99热网站在线观看| 国产免费又黄又爽又色| 国产精品欧美亚洲77777| av视频免费观看在线观看| 成人漫画全彩无遮挡| 国产欧美亚洲国产| 9191精品国产免费久久| 中文字幕另类日韩欧美亚洲嫩草| 国产精品嫩草影院av在线观看| 美女大奶头黄色视频| 黄网站色视频无遮挡免费观看| 男人操女人黄网站| 久久精品国产亚洲av天美| 美女国产视频在线观看| 国产又爽黄色视频| 欧美少妇被猛烈插入视频| 中国美白少妇内射xxxbb| 最近中文字幕高清免费大全6| 波野结衣二区三区在线| 国内精品宾馆在线| 亚洲精品视频女| 飞空精品影院首页| 草草在线视频免费看| 亚洲第一av免费看| av一本久久久久| 亚洲精品美女久久久久99蜜臀 | 99热这里只有是精品在线观看| 看非洲黑人一级黄片| av在线app专区| 午夜久久久在线观看| 亚洲在久久综合| 国语对白做爰xxxⅹ性视频网站|