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

    黑潮末端與黑潮延伸體的流場異常模態(tài)分析*

    2016-12-12 02:38:50姁,徐建,盧艷,張
    海洋科學(xué)進(jìn)展 2016年1期
    關(guān)鍵詞:黑潮中尺度年際

    盧 姁,徐 建,盧 艷,張 銘

    (1.中國人民解放軍 61741部隊(duì),北京 100094;2.解放軍理工大學(xué) 氣象海洋學(xué)院 大氣環(huán)流與短期氣候預(yù)測實(shí)驗(yàn)室,江蘇 南京 211101;3.北京市計(jì)量檢測科學(xué)研究院,北京 100029)

    ?

    黑潮末端與黑潮延伸體的流場異常模態(tài)分析*

    盧 姁1,2,徐 建1,2,盧 艷2,3,張 銘2

    (1.中國人民解放軍 61741部隊(duì),北京 100094;2.解放軍理工大學(xué) 氣象海洋學(xué)院 大氣環(huán)流與短期氣候預(yù)測實(shí)驗(yàn)室,江蘇 南京 211101;3.北京市計(jì)量檢測科學(xué)研究院,北京 100029)

    對北太平洋黑潮末端及黑潮延伸體的流場異常采用復(fù)EOF分析,結(jié)果發(fā)現(xiàn):該海域流場明顯異常在全年四季都主要發(fā)生在黑潮末端和黑潮延伸體上,其異常的主要形式在前2個(gè)模態(tài)均表現(xiàn)為一條由中尺度渦旋組成的渦旋帶,帶上通常有兩兩旋轉(zhuǎn)方向相反的渦旋偶。各季第一模態(tài)方差貢獻(xiàn)表現(xiàn)的季節(jié)變化特征與日本以南的黑潮流量相一致,秋季最大,夏季次之。ENSO的機(jī)制是副熱帶太平洋經(jīng)圈環(huán)流與熱帶太平洋緯向環(huán)流系統(tǒng)之間的相互作用,其在各季第一模態(tài)的年際變化上均有所反映。冬季PDO在年代際變化上表現(xiàn)明顯,其他三季則NPGO表現(xiàn)明顯。該海域流場異常的性質(zhì)是海洋Rossby波的異常,而赤道太平洋流場異常的性質(zhì)是赤道俘獲波的異常,但兩者均與風(fēng)應(yīng)力異常強(qiáng)迫有關(guān)。

    流場異常;黑潮末端;黑潮延伸體;復(fù)EOF分析

    黑潮是位于北太平洋西側(cè)的一支西邊界流,具有高溫、高鹽、流幅窄、流速快、流量大等特征。日本列島以南的黑潮通常稱之為黑潮末端。黑潮延伸體是指太平洋西邊界流在140°E,35°N附近與日本海岸分離后進(jìn)入北太平洋的水體,其呈現(xiàn)2個(gè)大彎曲形態(tài),體現(xiàn)了“蛇曲”的現(xiàn)象,且各年彎曲程度和彎曲位置均有較明顯的季節(jié)變化。黑潮延伸體的主軸南北兩側(cè)都有渦旋分布,這些渦旋各年也存在季節(jié)變化,在春季較弱,夏季最強(qiáng),秋末冬初較強(qiáng)。在離開日本海岸后,黑潮延伸體的流量增加,不穩(wěn)定性增強(qiáng),因此該區(qū)域是海洋中尺度渦旋最活躍的區(qū)域之一[1-2]。分析表明,黑潮末端和黑潮延伸體所在的海域是北太平洋流場異常最明顯的海域[3]。

    衛(wèi)星高度計(jì)資料目前已得到廣泛應(yīng)用,黑潮末端和黑潮延伸體所在海域也成為研究的熱點(diǎn)海域;學(xué)者利用該資料分析了該海域的SSH和SSHA,以了解流和中尺度渦旋及其異常的結(jié)構(gòu)和演變,探討了其季節(jié)和年際變化,并得到了很多有意義的結(jié)果[4-6]。在準(zhǔn)地轉(zhuǎn)的框架下,海表流場及其異常的結(jié)構(gòu)雖可由SSH和SSHA推得,然而畢竟不是直接的流動圖像,很難用于表層之下的上層海洋;此外,目前衛(wèi)星高度計(jì)資料的時(shí)段還不長,不能討論年代際變化。目前,在該海域,除利用衛(wèi)星高度計(jì)資料進(jìn)行研究外,利用SST資料在這方面所做的研究也很多[7-8],但SST受輻射、感熱等熱力因子的影響,其與流場的關(guān)系較復(fù)雜。次表層海溫與流場的關(guān)系較SST要單純,因其排除了熱力因子,故目前很多研究已轉(zhuǎn)到利用次表層的海溫資料上來[9-10]。在該海域采用年代較長且較可靠的海洋上層流場資料,對黑潮末端和黑潮延伸體的流場和中尺度渦旋直接開展研究,則可彌補(bǔ)衛(wèi)星高度計(jì)資料的不足,并特別有助于開展年代際變化的研究,但有關(guān)的研究尚不多見。

    年代際氣候變率已成為全球氣候中的熱點(diǎn)問題。太平洋海氣系統(tǒng)既有最強(qiáng)的年際信號,也存在明顯的年代際變化。研究表明,太平洋年際變化優(yōu)勢區(qū)位于熱帶太平洋,而年代際變化優(yōu)勢區(qū)位于中緯北太平洋,包含了黑潮末端和黑潮延伸體海域[11]。Mantua等把這種太平洋年代際振蕩現(xiàn)象稱為PDO(Pacific Decadal Oscillation)[12]。PDO可以用北太平洋SSTA的EOF第一模態(tài)來表征。目前關(guān)于PDO已有大量研究成果[13-14]。然而PDO還不能解釋東北太平洋中鹽度、營養(yǎng)物質(zhì)、含氯量及魚儲量等的年代際變化。為此,2008年Di Lorenzo定義了一個(gè)新的氣候模態(tài)NPGO (North Pacific Gyre Oscillation)[15],NPGO是東北太平洋SSHA的EOF分解第二模態(tài),因其與SSTA的EOF分解第二模態(tài)有很好相關(guān)[16],故也可用SSTA的EOF第二模態(tài)來表征。PDO和NPGO是北太平洋2個(gè)主要的氣候模態(tài)。

    因此,本文采用較可靠的海洋流場資料,對黑潮末端和黑潮延伸體海域的流場做了復(fù)EOF診斷,分析了其第一、二模態(tài)的空間場和時(shí)間系數(shù),以便探討黑潮末端、黑潮延伸體和中尺度渦旋流場異常的結(jié)構(gòu),給出其年際特別是年代際變化,并對各模態(tài)與ENSO、PDO和NPGO的關(guān)系做了討論。

    1 資料和診斷方法

    本文所用的資料為美國UMD的SODA 2.2.4逐月全球海洋同化分析資料[17](http://soda.tamu.edu/assim/SODA_2.2.4/),該資料提供了1871-2008年深度在5~5 375 m共40層上的各月平均洋流,網(wǎng)格距為0.5°×0.5°。此前我們曾用該資料對1951-2008年共58 a的赤道外北太平洋上層洋流的異常做過分析[3];分析中所取的海洋深度為:112,96,82,70,57,46,35,25,15,5 m,因?yàn)榇?0層上資料的可靠性較好。分析發(fā)現(xiàn),北太平洋流場異常明顯的海域位于(27°15′~39°45′N, 134°15′~154°45′E),該海域中包含了黑潮末端與黑潮延伸體以及中尺度渦旋,且將該海域稱之為關(guān)鍵區(qū)[3]。本文研究的范圍就取此關(guān)鍵區(qū),而時(shí)間段也取1951-2008年。

    為揭示此關(guān)鍵區(qū)內(nèi)的全年上層洋流的異常規(guī)律,本文對各月該洋流異常做了復(fù)EOF分析,其原理可參見文獻(xiàn)[18-19],本文不再贅述。具體操作步驟:首先對上述10個(gè)深度上的各年各月平均流場按月求其58 a的平均值,將各年的各月平均流場相應(yīng)減去該平均值,則可得各年各月這10層上的偏差流場;然后將各月這10層上的偏差流場作為一個(gè)整體進(jìn)行復(fù)EOF分析,這樣各層有相同的時(shí)間系數(shù),這也是該復(fù)EOF分析方法的特色和優(yōu)點(diǎn)。此時(shí)各月各EOF模態(tài)的空間場和時(shí)間系數(shù)場都是復(fù)數(shù);各月空間場的模表示各月各模態(tài)流場異常的流速大小,輻角則表示其流向;而時(shí)間系數(shù)表示這58 a中各月流速、流向隨年份的變化。

    本文主要討論區(qū)域深度為112和15 m的2層。前者是本文所取的最深層次,在該深度海溫變化已很難受到輻射、感熱等熱力因子影響;其主要由海水的垂直運(yùn)動來決定,上升處有降溫,下沉處有增溫,此即為海溫的動力變化。在以上10層中,因前者處于次表層水團(tuán)中,以下為敘述方便,就直接稱其為次表層(次表層水是在大洋表層之下,以躍層為界形成的水團(tuán),其厚度通常為200~300 m);而后者可稱之為近表層。

    2 診斷結(jié)果分析

    復(fù)EOF分析的結(jié)果表明,所得到的前2個(gè)模態(tài)都通過了顯著性檢驗(yàn)[20]。各月第一、二模態(tài)的方差貢獻(xiàn)見表1,前者的各月平均為21.4%,是最重要的模態(tài),后者為13.5%,仍具重要性。本文以下僅給出第一、二模態(tài)的結(jié)果并進(jìn)行討論。

    由表1還可見,第一模態(tài)在秋季(9-11月)的方差貢獻(xiàn)最大,夏季(6-8月)次之,冬季(12月至翌年2月)最??;第二模態(tài)在夏季最大,秋季次之,冬季最小。夏季和秋季的流場異常要比冬季和春季顯著(參見表1中累計(jì)方差一欄)。研究表明[21],日本以南(137°E斷面)的黑潮流量表現(xiàn)出秋季最大,夏季次之的季節(jié)特點(diǎn),且1—7月的流量較年平均流量要偏小,而8—12月的則偏大。這與第一模態(tài)的方差貢獻(xiàn)相一致。以下分別給出1,4,7,10月份的復(fù)EOF分析結(jié)果,并分別代表冬、春、夏、秋四季。

    表1 北太平洋各月各模態(tài)的方差貢獻(xiàn)及其累計(jì)方差貢獻(xiàn)Table 1 Variance and accumulated variance contributions of the two modes in the North Pacific in each calendar month

    2.1 第一模態(tài)

    在北太平洋上層,從診斷得到的各季第一模態(tài)空間場的總體分布可知:在上述10層上,明顯的偏差流場出現(xiàn)在黑潮末端和黑潮延伸體上,表現(xiàn)為一條由中尺度渦旋組成的渦旋帶,這些渦旋的尺度(直徑)約為200~400 km,并通常兩兩構(gòu)成旋轉(zhuǎn)方向相反的渦旋偶形式;距本州島東海岸最近的那些渦旋其強(qiáng)度最大,且它們的軸心呈東北西南走向,大致平行于該海岸線,在黑潮延伸體上的渦旋帶則呈緯向分布(圖1)。隨著季節(jié)變化,上述渦旋帶中渦旋的形態(tài)也有些變化,春、夏季這些渦旋的強(qiáng)度更大些。在這10層上該空間場的差異很小,這表明海洋上層的流場異常具有正壓性,該處海水可看作是不可壓縮的,為此該模態(tài)次表層的圖略。

    圖1 近表層第一模態(tài)的空間場(m·s-1)

    從各季第一模態(tài)時(shí)間系數(shù)春、夏、秋季的輻角集中在0°和±180°附近;冬季輻角分布的離散度略大,但仍集中在-50°和130°附近。這表明四季輻角的分布均大體對應(yīng)于2個(gè)狀態(tài),這2個(gè)狀態(tài)也可用數(shù)學(xué)符號“+”、“-“來標(biāo)注。第一模態(tài)時(shí)間系數(shù)各季的模值各年也有所不同。

    對于輻角分布存在2個(gè)狀態(tài)的情況,可將輻角和模的這2個(gè)時(shí)間系數(shù)序列綜合為一個(gè)時(shí)間序列。具體的做法是:將春、夏、秋三季每年時(shí)間系數(shù)的模值乘以該年輻角的余弦,而冬季則每年將時(shí)間系數(shù)的模乘以該年輻角加上50°的余弦。如此各年四季的余弦值均接近±1,這樣對冬、春、夏、秋四季,各年均能得到一個(gè)新的實(shí)數(shù)序列,稱之為實(shí)時(shí)間系數(shù)序列,該序列能綜合反映輻角和模的時(shí)間演變。

    此方法僅適用于時(shí)間系數(shù)的輻角分布在兩個(gè)狀態(tài)的情況。此時(shí)對于春、夏、秋三季,因其第一模態(tài)輻角分布在0°和±180°,故有結(jié)論:該模態(tài)各年的流向就約等于此模態(tài)空間場的流向,而其強(qiáng)度則由實(shí)時(shí)間系數(shù)序列值的絕對值來決定,若該序列值為負(fù),則流向要反向。對冬季,因第一模態(tài)輻角分布在-50°和130°,故首先要對此模態(tài)空間場的流向做-50°也即順時(shí)針50°的旋轉(zhuǎn),然后將旋轉(zhuǎn)流向后的場作為新的空間場,對該新空間場,則有以上相同的結(jié)論。由于旋轉(zhuǎn)角度(50°)不算太大,故該新空間場中,上述渦旋帶中渦旋的數(shù)目不變,但各渦旋的散度絕對值減小而渦度絕對值增大,新空間場上中尺度渦旋表現(xiàn)得更明顯。對于時(shí)間系數(shù)的輻角分布在3個(gè)及以上狀態(tài)的情況,則不能得到實(shí)時(shí)間系數(shù)序列,此時(shí)必須要用輻角和模的時(shí)間系數(shù)序列來討論問題。

    圖2為各季第一模態(tài)的實(shí)時(shí)間系數(shù)序列圖。由圖2可知,該序列有明顯的年際變化和年代際變化。為了更準(zhǔn)確地反映這種變化,現(xiàn)對各季第一模態(tài)實(shí)時(shí)間系數(shù)序列做了小波分析。圖3為其小波全譜,從該圖可見:冬季有約6 a的年際變化,12~13 a和18~19 a的年代際變化,并以最后者為最顯著;春、夏、秋三季都有約7 a的年際變化和約15 a的年代際變化,均以后者最顯著。總之,各季年代際變化都較年際變化要顯著。

    圖2 第一模態(tài)實(shí)時(shí)間系數(shù)圖

    2.2 第二模態(tài)

    各季第二模態(tài)的空間場分布見圖4,其總體仍表現(xiàn)為上述渦旋帶上的渦旋,但渦旋分布形式與第一模態(tài)有些不同,兩者渦旋有所錯(cuò)位(可對照圖1與圖4)。第二模態(tài)與第一模態(tài)一樣,各層上空間場的差異也很小,流場異常的正壓性明顯,因此次表層的圖也略。

    圖3 第一模態(tài)實(shí)時(shí)間系數(shù)序列的小波全譜

    圖4 近表層第二模態(tài)的空間場(m·s-1)

    與第一模態(tài)相似,各季第二模態(tài)時(shí)間系數(shù)的輻角分布也有2個(gè)狀態(tài),不過其在這2個(gè)狀態(tài)上分布的離散度要較第一模態(tài)大些。各季第二模態(tài)時(shí)間系數(shù)的模值大小各年也都有明顯差別。采用與第一模態(tài)同樣的做法,也能得到各季第二模態(tài)的實(shí)時(shí)間系數(shù)序列(圖5);對該序列也做了小波分析,并可見第二模態(tài)各季的小波全譜總體呈單峰分布,年際變化均不很顯著,而年代際變化都十分顯著;冬季有約20 a,春、夏、秋三季均有約15 a的年代際變化周期(圖6)。

    圖5 第二模態(tài)實(shí)時(shí)間系數(shù)圖

    圖6 第二模態(tài)實(shí)時(shí)間系數(shù)序列的小波全譜

    3 有關(guān)中尺度渦旋的討論

    綜上所述,黑潮末端和黑潮延伸體海域是北太平洋流場異常最明顯的海域。在該海域即關(guān)鍵區(qū)中海洋上層58 a的各月份平均流場上可見,各月份各層次的平均流場差異都很小,因此這里僅給出1,7月(冬、夏季)近表層的平均流場圖(圖7)。圖7中最主要的流系有3個(gè):呈“蛇曲”(波動)狀的黑潮延伸體,其波動振幅向大洋中部衰減,“蛇曲”上的兩個(gè)主脊分別位于143°30′E和147°36′E,兩脊之間的主槽則位于145°30′E。位于日本本州島南面的黑潮末端。黑潮末端南側(cè)中心在(30°18′N,135°18′E)的橢圓狀反氣旋渦旋,即回流渦旋。

    圖7 冬、夏季多年平均流場(m·s-1)圖

    黑潮末端和黑潮延伸體均呈現(xiàn)急流形式。在多年平均上,這3個(gè)系統(tǒng)都很穩(wěn)定,冬夏變化很小。除以上回流渦旋外,冬、夏季在黑潮延伸體急流以南還有較弱的3個(gè)中尺度反氣旋渦旋,中心分別位于(33°N, 142°36′E)、(33°N,148°E)、(33°N,151°36′E)附近(圖7)。因該圖是多年平均圖,給出了中尺度渦旋分布的背景狀況,各年各月的具體情況則會有所不同,這由流場異常各模態(tài)的空間場和時(shí)間系數(shù)來決定。

    從各季第一、二模態(tài)的空間場可見(圖1,4),流場異常都表現(xiàn)為黑潮末端和黑潮延伸體上的渦旋帶處,在該帶之外流場異常小,在該帶上則多表現(xiàn)為旋轉(zhuǎn)方向相反的中尺度渦旋偶;這表明中尺度渦旋運(yùn)動是流場異常的最主要形態(tài)。比較平均流場(圖7)與異常流場的第一、二模態(tài)(圖1,4)可見,兩者形式完全不同,前者以急流形式為主,后者則以中尺度渦旋為主,這表明各月關(guān)鍵區(qū)中3個(gè)流系的變化(異常),主要由各模態(tài)中的中尺度渦旋來決定。從動能分布來看,平均場動能以急流的動能為主,而偏差場動能則以渦旋動能(Eddy Kinetic Energy, EKE)為主,也即動能異常主要由偏差場的中尺度渦旋動能來體現(xiàn)。

    實(shí)時(shí)觀測表明,在水深200 m以上,有多個(gè)中尺度渦旋分布在黑潮末端和黑潮延續(xù)體的急流兩側(cè),并也構(gòu)成帶狀;在此急流與渦旋并存,這是黑潮末端和黑潮延續(xù)體上的顯著特點(diǎn)[23-24]??傮w而言,這些渦旋在春季較弱,夏季最強(qiáng),秋末冬初較強(qiáng)。需注意的是,以上特點(diǎn)看來與四季第一、二模態(tài)的空間場有些類似,但在該空間場上出現(xiàn)的那些渦旋是復(fù)EOF分析的結(jié)果,反映了該處的流場異常;故與觀測的結(jié)果不能混淆。不過觀測場的多年平均就是氣候平均場(背景場),而各模態(tài)空間場則是該背景場的偏差(異常),故兩者分布有些類似是很自然的。因觀測和復(fù)EOF得到的中尺度渦旋(兩者合稱為前者)都分布在黑潮末端和黑潮延伸體(這兩者合稱為后者)急流兩側(cè),故前者的渦旋動能應(yīng)與后者密切有關(guān);至于前者產(chǎn)生的原因,有學(xué)者認(rèn)為與后者上急流的正、斜壓不穩(wěn)定有關(guān),為背景流的失穩(wěn)所導(dǎo)致[25-26],此外風(fēng)應(yīng)力的強(qiáng)迫也是前者產(chǎn)生的原因之一[27];本文復(fù)EOF的診斷結(jié)果表明,各模態(tài)流場異常的正壓性明顯;且后者呈急流形式(圖7),滿足正壓不穩(wěn)定的必要條件,故而支持正壓不穩(wěn)定的說法,當(dāng)然也不排除上述其他原因;至于這些原因中何種是最主要的則有待于更深入的研究。

    第一、二模態(tài)中主要表現(xiàn)為中尺度渦旋形式的流場異常,其仍是準(zhǔn)地轉(zhuǎn)的,這是由于其水平尺度L為100 km,其流速尺度U為1 m/s,取地轉(zhuǎn)參數(shù)f=10-4s-1,則其Rossby數(shù)R0=U/(fL)=0.1≤1的緣故;還可知,其出現(xiàn)在一條渦旋帶中,而此帶處于氣候平均的黑潮末端和黑潮延伸體的急流上(將圖1,4與圖7做比較);如此則其主要體現(xiàn)了該急流的異常??烧J(rèn)為該急流是由海底地形導(dǎo)致的海洋Rossby波[28](這與大氣中氣流過山形成的大氣Rossby波相類似),故其反映了該海洋Rossby波的異常。

    4 各模態(tài)的年際和年代際變化

    4.1 季節(jié)與年際變化

    各模態(tài)的季節(jié)變化比較明顯,這體現(xiàn)在各模態(tài)四季的空間場上渦旋帶中的渦旋分布都不相同。以第一模態(tài)冬、夏(1,7月)為例,兩者空間場上最大的不同在于冬季在本州島以東海域(33°30′N以北,144°30′E以西)有2個(gè)旋轉(zhuǎn)方向相反的中尺度強(qiáng)渦旋,組成渦旋偶,在本州島以南海域有1個(gè)中尺度強(qiáng)渦旋;而夏季在該島以東的該海域只有1個(gè)中尺度強(qiáng)渦旋,在該島以南則有2個(gè)旋轉(zhuǎn)方向相反的中尺度強(qiáng)渦旋,并組成渦旋偶。造成這種狀況的原因是在冬、夏的風(fēng)向異常兩者不同,前者主要為偏西風(fēng)異常,后者主要為偏南風(fēng)異常。

    在年際變化上,黑潮延伸區(qū)的海平面異常和海表面溫度異常的低頻分量與Nino3指數(shù)具有較高的相關(guān)性[29],而各季第一模態(tài)大致都有6~7 a的明顯年際變化,其中夏、秋兩季更要顯著(圖3c,3d),這應(yīng)是ENSO在年際變化上的反映。以上結(jié)果表明黑潮末端和黑潮延伸體流場的年際變化參與了ENSO循環(huán),而黑潮與熱帶太平洋緯向環(huán)流系統(tǒng)間的相互作用是溝通ENSO循環(huán)的內(nèi)在環(huán)節(jié)[30]。第二模態(tài)實(shí)時(shí)間系數(shù)小波全譜圖中年際變化則不如第一模態(tài)明顯,這說明ENSO對該模態(tài)影響較第一模態(tài)要小。

    4.2 年代際變化

    第一模態(tài)冬季有約20 a的最強(qiáng)年代際變化周期。其與北太平洋年代際濤動(PDO)的周期相同,應(yīng)為PDO在流場異常上的反映[31],而冬季阿留申低壓也最強(qiáng)大。約15 a的年代際變化在各季都很明顯,除冬季外,在小波全譜圖上該周期成為主峰,特別在夏秋季表現(xiàn)則更突出。該約15 a的周期則與NPGO約13 a的周期很接近,應(yīng)為NPGO在流場異常上的反映,并與副熱帶太平洋環(huán)流有關(guān)[32],而該環(huán)流在夏季強(qiáng)度最強(qiáng),位置也最偏北。

    在年代際變化上,各季第二模態(tài)也有約20 a的周期,也應(yīng)為PDO在流場異常上的反映。除冬季外,其他三季最明顯的年代際變化同樣約為15 a,為NPGO在流場異常上的反映。因NPGO與SSHA與SSTA的第二模態(tài)有關(guān)[15-16],對于氣候尺度,表層流場異常與SSHA有準(zhǔn)地轉(zhuǎn)關(guān)系,在一定條件下與SSTA也有診斷關(guān)系[7,32]。既然SSHA與SSTA的第二模態(tài)能反映NPGO,則除冬季外流場異常的第二模態(tài)反映了NPGO也很自然。冬季因阿留申低壓十分強(qiáng)大,且位置偏南,故NPGO為PDO所掩蓋,在小波全譜圖上約15 a的周期則成為次高峰。

    以上對第一、二模態(tài)的分析表明,黑潮末端和黑潮延伸體海域(關(guān)鍵區(qū))的年代際變化與北太平洋流(North Pacific Current, NPC)以及該流以北的阿拉斯加環(huán)流(Alaskan Gyre)和以南的副熱帶環(huán)流(Subtropical Gyre)的異常有關(guān);該兩個(gè)環(huán)流的異常又分別與大氣中的阿留申低壓異常和太平洋副熱帶高壓異常有關(guān),而它們之間的海氣相互作用應(yīng)是造成PDO和NPGO的原因。本文的關(guān)鍵區(qū)則抓住了這些流動的關(guān)鍵部分。綜上,冬季第一模態(tài)為PDO的主模態(tài),第二模態(tài)為PDO的次模態(tài),其他季節(jié)第一模態(tài)為NPGO的主模態(tài),第二模態(tài)為NPGO的次模態(tài)。

    前面已發(fā)現(xiàn)日本以南的黑潮流量與第一模態(tài)方差貢獻(xiàn)的季節(jié)變化特征相一致,根據(jù)Sverdrup等[33]、Munk[34]和Stommel[35]提出的有關(guān)均質(zhì)海洋的風(fēng)生海洋環(huán)流理論,作為北太平洋副熱帶環(huán)流中的西邊界流黑潮,其流量由北太平洋副熱帶海面風(fēng)應(yīng)力旋度所決定。研究表明[13,36],年代際振蕩的時(shí)間尺度取決于副熱帶海洋環(huán)流的調(diào)整時(shí)間,而后者的緩慢調(diào)整涉及兩個(gè)物理過程,一個(gè)是大氣風(fēng)應(yīng)力旋度強(qiáng)迫產(chǎn)生的海洋Rossby波向西傳播,在若干年后到達(dá)西邊界,從而影響副熱帶環(huán)流西邊界流的向極輸送以及向黑潮末端和黑潮延伸體的熱量輸送;另一個(gè)是副熱帶環(huán)流對溫度的平流輸送過程[37]。

    5 北太平洋與熱帶太平洋洋流異常的比較

    文獻(xiàn)[38]中我們研究了熱帶太平洋上層洋流異常(前者),這里將前者與北太平洋黑潮末端和黑潮延伸體海域(關(guān)鍵區(qū))的上層洋流異常(后者)作一比較;從明顯異常的緯向范圍看,前者異常分布在整個(gè)赤道太平洋海域,范圍較后者要大。從異常的性質(zhì)看,前者是海洋赤道波動(含Kelvin波)的異常,其表現(xiàn)為赤道俘獲波的形式;而后者上面已論及,有R0≤1,是準(zhǔn)地轉(zhuǎn)的,并是海洋Rossby波的異常,該異常表現(xiàn)為中尺度渦旋的形式。

    該兩者雖有以上不同之處,但也有共同的地方,即這兩者異常則均與風(fēng)應(yīng)力異常密切有關(guān)。有研究表明,中緯度風(fēng)應(yīng)力旋度強(qiáng)迫的海洋Rossby波西傳到太平洋西邊界可能和向赤道傳播的海岸Kelvin波耦合,從而使中緯度信號與熱帶信號相關(guān)聯(lián)[39-40],這就將本文對北太平洋關(guān)鍵區(qū)的研究與以前對熱帶太平洋的研究聯(lián)系到一起了,同時(shí)這也說明了本文第一模態(tài)的年際變化受ENSO影響大的原因。

    6 結(jié) 語

    本文采用復(fù)EOF分析方法,對全年四季北太平洋黑潮末端和黑潮延伸體海域(關(guān)鍵區(qū))的海洋上層流場異常做了統(tǒng)計(jì)動力診斷,得到了以下主要結(jié)論:

    1)各季該海域上層流場的明顯異常主要發(fā)生在黑潮末端和黑潮延伸體上即關(guān)鍵區(qū)中,其異常的第一、二模態(tài)均表現(xiàn)為一條由中尺度渦旋組成的渦旋帶,帶上通常有兩兩旋轉(zhuǎn)方向相反的渦旋偶,這是該流場異常的主要形式。

    2)第一模態(tài)在秋季的方差貢獻(xiàn)最大,夏季次之,冬季最??;第二模態(tài)則在夏季最大,秋季次之,冬季也最??;各季第一模態(tài)方差貢獻(xiàn)表現(xiàn)的季節(jié)變化特征與日本以南的黑潮流量相一致。

    3)在第一模態(tài)的年際變化上,ENSO在各季均有所反映,其機(jī)制是副熱帶太平洋經(jīng)圈環(huán)流與熱帶太平洋緯向環(huán)流系統(tǒng)之間的相互作用;而第二模態(tài)ENSO的影響不大。

    4)在年代際變化上,冬季PDO表現(xiàn)明顯,其他三季NPGO表現(xiàn)明顯;冬季第一模態(tài)為PDO的主模態(tài),第二模態(tài)為PDO的次模態(tài),其他三季第一模態(tài)為NPGO的主模態(tài),第二模態(tài)為NPGO的次模態(tài)。

    5)以上兩模態(tài)異常的性質(zhì)是海洋Rossby波的異常,這與赤道太平洋是赤道俘獲波(含Kelvin波)的異常不同,但兩者則均與風(fēng)應(yīng)力異常強(qiáng)迫有關(guān)。

    最后要說明的是,本文主要工作是進(jìn)行診斷分析,對造成異常的動力學(xué)原因涉及不夠,這是本文局限所在,也是今后要做的工作。

    [1] QIU B, CHEN S, HACKER P, et al. The Kuroshio extension northern recirculation gyre: Profiling float measurements and forcing mechanism[J].Journal of Physical Oceanography,2008,38(8):1764-1779.

    [2] SCHARFFENBERG M G, STAMMER D. Annual variations of geostrophic currents and eddy kinetic energy inferred from TOPEX/Poseidon and Jason-1 tandem mission data[J].Ocean Sciences Meeting, Orlando,2008:3-7.

    [3] LU K C, LU X, ZHANG M. Analysis for upper abnormal current in North Pacific without equatorial area[J].Marine Science Bullent, 2011, 30(1): 29-36. 路凱程,盧姁,張銘.赤道外北太平洋上層洋流異常分析[J].海洋通報(bào),2011,30(1):741-748.

    [4] DBCHELTON D B, SCHLER M G, SAMELSON R M. Global observations of nonlinear mesoscale eddies[J].Progress in Oceanography,2011,91(2):167-216.

    [5] TAGUCHI B, QIU B, NONAKA M, et al. Decadal variability of the Kuroshio extension: Mesoscale eddies and recirculations[J].Ocean Dynamics,2010,60(3):673-691.

    [6] ZHANG Y C, ZHANG L F, WANG Y G. Interannual sea level variability in the North Pacific Ocean and its mechanisms[J].Journal of Geophysics, 2010, 53(2): 247-255. 張永垂,張立鳳,王業(yè)桂. 北太平洋海表面高度的年際變化及其機(jī)制[J].地球物理學(xué)報(bào),2010,53(2):247-255.

    [7] Lü Q P, LU K C, ZHANG M. NPGO mode of the upper sea temperature anomalies in the North Pacific during winter[J].Climatic and Environmental Research, 2013, 18(2): 210-220. 呂慶平,路凱程,張銘.北太平洋冬季上層海溫異常的NPGO模態(tài)[J].氣候與環(huán)境研究,2013,18(2):210-220.

    [8] WANG L, LI T, ZHOU T J. Interseasonal SST variability and air-sea interaction over the kuroshio extension region during boreal summer[J].Journal of Climate,2012,25(5):1619-1634.

    [9] WU H Y, LI C Y. The preliminary reseach of equatorial Pacific-Indian Ocean temperature anomaly mode and subsurface ocean temperature anomaly[J].Acta Oceanologica Sinica,2009,31(2):24-33. 吳海燕,李崇銀.赤道太平洋-印度洋海溫異常綜合模與次表層海溫異常[J].海洋學(xué)報(bào),2009,31(2):24-33.

    [10] XU K, ZHU C W. Estimation of the total dust column and dry eeposition flux over the Yellow Sea, China based on shipboard sun photometer measurements:case study[J].Atmospheric and Oceanic Science Letters,2010,3(2):106-110.

    [11] GIESE B S, CARTON J A. Interanuual and decadal variability in the tropical and midlatitude Pacific Ocean[J].Journal of Climate,1999,12(12):3402-3418.

    [12] MANTUA N J, HARE S R, ZHANG Y, et al. A Pacific interdecadal climate oscillation with impacts on salmon production[J].Bulletion of the American Meteorological Society,1997,78(6):1069-1079.

    [13] YANG X Q, ZHU Y M,XIE Q, et al. Advances in studies of Pacific Decadal Oscillation[J].Chinese Journal of Atmospherice, 2004, 28(6): 979-992. 楊修群,朱益民,謝倩,等. 太平洋年代際振蕩的研究進(jìn)展[J].大氣科學(xué),2004,28(6):979-992.

    [14] ALEXENDER M A. Extratropical air-sea interaction, SST variability and the Pacific Decadal Oscillation[M].Climate dynamics: why does climate vary? Sun D, Bryan F. Washington D C: Am. Geophys Uninon,2010:123-148.

    [15] LORENZO D E. North Pacific Gyre Oscillation links ocean climate and ecosystem change[J].Geophys Research Letters, 2008,35(8):1-6.

    [16] BOND N A, OVERL J E, SPILLANE M, et al. Recent shifts in the state of the North Pacific[J].Geophysical Research Letters,2003,30(23):2183,doi:10.1029/2003GL018597.

    [17] CARTON J A, GRIESE B S. A reanalysis of ocean climate using simple ocean data assimilation(SODA)[J].Monthly Weather Review,2008,136(8):2999-3017.

    [18] ZENG Q C. The theory of atmospheric infrared remote sense[M].Beijing: China Scientific Press,1974:160-166.曾慶存. 大氣紅外遙感原理[M].北京:科學(xué)出版社,1974:160-166.

    [19] HUANG J Y. Statistic analysis and forecast methods in meteorology[M].Beijing: China Meteorological Press, 2000: 135-139. 黃嘉佑. 氣象統(tǒng)計(jì)分析與預(yù)報(bào)方法[M].北京:氣象出版社,2000:135-139.

    [20] NORTH G R, BELLl T L, CAHALAN R F, et al. Sampling errors in the estimation of empirical orthogonal function[J].Monthly Weather Review,1982,110(7):699-706.

    [21] CAI R S, ZHANG Q L, QI Q H. Character of transport variations at the source and adjacent area of Kuroshio[J].Journal of Oceanography in Taiwan Strait, 2009, 28(3): 299-307.蔡榕碩,張啟龍,齊慶華.源地黑潮及其上下游流量的變化特征[J].臺灣海峽,2009,28(3):299-307.

    [22] LIU T, CHEN C C, MILO L. Accurate standard cell characterization and statistical timing analysis using multivariate adaptive regression splines[C]∥IEEE.Proc. Int. Symp. Quality Electronic Systems, 2015:272-279.

    [23] WANG D X, FANG G H, WANG Q, et al. The tropical Pacific Ocean circulation variability and the interaction of air and sea[M]. Beijing:Ocean Press,2009:310-313. 王東曉,方國洪,王啟,等. 熱帶太平洋環(huán)流變異與海氣相互作用[M].北京:海洋出版社,2009:310-313.

    [24] WATERMAN S, HOGG N G, JAYNE S R. Eddy-mean flow interaction in the Kuroshio extension region[J].Journal of Physical Oceanography,2011,41(6):1182-1208.

    [25] STAURT P B, FRANK O B. A comparison of mesoscale eddy heat fluxes from observations and a high-resolution ocean model simulation of the Kuroshio extension[J].Journal of Physical Oceanography,2013,43:2563-2570.

    [26] SYUART P B, RANDOLPH D W, KATHLEAN A D. Divergent eddy heat fluxes in the Kuroshio extension at 144°-148° E. Part I: Mean Structure[J].Journal of Physical Oceanography,2013,43:1533-1550.

    [27] Lü Q P, ZHANG W F, ZHANG M. Study of two-layer ocean model forced by the climate wind anomalies[J].Climatic and Environmental Research, 2013, 18(1): 124-134. 呂慶平,張維鋒,張銘.兩層海洋對風(fēng)場氣候異常響應(yīng)的解析解及其討論[J].氣候與環(huán)境研究,2013,18(1):124-134.

    [28] MIZUON K, WHITE W B. Annual and interannual variability in the Kuroshio current system[J].Journal of Physical Oceanography,1983,13(10):1847-1867.

    [29] KRAMER W, DIJKSTRA H A. Optimal localized observations for advancing beyond the ENSO predictability barrier[J].Nonlinear Processes in Geophysics,2013,20(2):221-230.

    [30] WANG D X, XIE Q, LIU Y, et al. The research progress of the Pacific decadal variability[J].Journal of Tropical Oceanoraphy,2003,22(1):76-83. 王東曉,謝強(qiáng),劉赟,等.太平洋年代際海洋變率研究進(jìn)展[J].熱帶海洋學(xué)報(bào),2003,22(1):76-83.

    [31] QIU B, CHEN S. Eddy-mean flow interaction in the decadally modulating Kuroshio Extension system[J].Deep Sea Research,2010,57(13-14):1098-1110.

    [32] ZHANG L F, Lü Q P. The relationship between the npgo mode and npo mode in the north pacific during winter[EB/OL]. [2014-08-10].http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGQX201111004045.htm.張立鳳,呂慶平.冬季北太平洋NPGO模態(tài)與大氣NPO模態(tài)的關(guān)系[EB/OL]. [2014-08-10]. http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGQX201111004045.htm.模態(tài)與大氣NPO模態(tài)的關(guān)系[EB/OL].[2014-08-10].http://cpfd.cnki.com.cn/Article/CPFDTOTAL-ZGQX201111004045.htm.

    [33] SVERDRUP H U, JOHNSON M W, FLEMING R H. The oceans: Their physics,chemistry and general biology[M].New Jersey:Prentice Hall,Englewood Ciliffs,1942.

    [34] MUNK W H. On the wind-dricen ocean circulation[J].Journal of Meteorological,1950,7(3):79-93.

    [35] STOMMEL H. The westward intensification of wind-driven ocean currents[J].Trans of the American Geophys Union,1948,29(2):202-206.

    [36] JASON C, FURTADO, EMANUELE DL, et al. North Pacific decadal variability and climate change in the IPCC AR4 models[J].Journal of Climate,2011,24(2):3049-3067.

    [37] ROBERTSONA W. Interdecadal variability over the North Pacific in a multi-century climate simulation[J].Climate Dynamic,1996,12:227-241.

    [38] LU X, ZHANG D L. A statistical diagnosis of upper-ocean currents in tropical Pacific Ocean in May[J].Journal of Tropical Oceanography,2009,28(2):22-30.盧姁,張東凌.熱帶太平洋5月份上層洋流的動力統(tǒng)計(jì)診斷[J].熱帶海洋學(xué)報(bào),2009,28(2):22-30.

    [39] JIN F F, KINOTO M, WANG X C. A model of decadal ocean-atmosphere interaction in the North Pacific basin[J].Geophysical Research Letters,2001,28(8):1531-1534.

    [40] CAPOTONDI A, ALEXANDER M A. Rossby waves in the tropical North Pacific and their role in decadal thermocline variability[J].Journal of Physical Oceanography, 2001, 31(12): 3496-3515.

    Received: September 22, 2014

    Analysis of Abnormal Current in Kuroshio End and Kuroshio Extension

    LU Xu1,2, XU Jian1,2, LU Yan2,3,ZHANG Ming2

    (1.Unit61741,P.L.A,Beijing 100094,China;2.LaboratoryofAtmosphericCirculationandShort-rangeClimateForecastMeteorologicalCollege,P.L.A.UniversityofScienceandTechnology, Nanjing 211101, China; 3.BeijingInstituteofMetrology,Beijing 100029,China)

    Anomalous circulations of four seasons in the upper layer in Kuroshio end and Kuroshio extension are studied using the complex empirical orthogonal function (CEOF) method. The results show that the anomalous circulation in the North Pacific mainly occurs in Kuroshio end and Kuroshio extension in all the seasons, with major anomalies represented by the first and the second modes with a vortex belt consisting of mesoscale eddies, which usually shows pairs of vortices rotating in opposite direction. Seasonality of variance contribution the first mode in each season is consistent with the strength of Kuroshio to the south of Japan, which is greatest in fall and secondary in summer. ENSO has influence on interannual variations of the first mode in all seasons through interactions between subtropical meridional cell and tropical zonal cell in the Pacific. In interdecadal variations, PDO's influence is obvious in winter while NPGO's is obvious in the other three seasons. The anomalous circulations represented by the two modes may be resulted from oceanic Rossby waves, which are different from those trapped waves in equatorial Pacific, but both have a close relation with anomalous wind stresses.

    anomalous circulation; Kuroshio end; Kuroshio extension; CEOF

    2014-09-22

    國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目——西北太平洋海洋多尺度變化過程、機(jī)理及可預(yù)報(bào)性 (2013CB956203)

    盧 姁(1982-),女,江蘇南通人,工程師,博士,主要從事氣候方面研究.E-mail: xlu2006918@163.com

    (陳 靖 編輯)

    P732.6

    A

    1671-6647(2016)01-0010-12

    10.3969/j.issn.1671-6647.2016.01.002

    猜你喜歡
    黑潮中尺度年際
    譜松弛動力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動研究
    基于深度學(xué)習(xí)的中尺度渦檢測技術(shù)及其在聲場中的應(yīng)用
    2016年7月四川持續(xù)性強(qiáng)降水的中尺度濾波分析
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計(jì)數(shù)據(jù)的黑潮變異特征
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    伊人久久大香线蕉亚洲五| 国产成人一区二区在线| 99热国产这里只有精品6| 成年人免费黄色播放视频| 女的被弄到高潮叫床怎么办| 在线观看国产h片| 十八禁人妻一区二区| 五月开心婷婷网| h视频一区二区三区| 久久精品国产a三级三级三级| a级毛片黄视频| 婷婷色综合www| 最近中文字幕高清免费大全6| 国产麻豆69| 免费高清在线观看视频在线观看| 亚洲欧美精品自产自拍| 看免费成人av毛片| kizo精华| 亚洲,欧美精品.| 亚洲成av片中文字幕在线观看| 午夜影院在线不卡| 久久国产精品男人的天堂亚洲| 亚洲精品一区蜜桃| 一边摸一边做爽爽视频免费| 精品少妇内射三级| 欧美日韩成人在线一区二区| 国产成人精品福利久久| 亚洲精品久久成人aⅴ小说| 国产成人系列免费观看| 国产无遮挡羞羞视频在线观看| www.精华液| 如日韩欧美国产精品一区二区三区| 欧美日本中文国产一区发布| 亚洲精品一区蜜桃| 夜夜骑夜夜射夜夜干| av又黄又爽大尺度在线免费看| 人成视频在线观看免费观看| 亚洲视频免费观看视频| 麻豆精品久久久久久蜜桃| 日本av手机在线免费观看| 欧美人与性动交α欧美精品济南到| 国产精品.久久久| 国产亚洲精品第一综合不卡| 亚洲欧美激情在线| 精品少妇一区二区三区视频日本电影 | 久久精品亚洲熟妇少妇任你| 国产成人a∨麻豆精品| 国产视频首页在线观看| 日本wwww免费看| 这个男人来自地球电影免费观看 | 丰满乱子伦码专区| 又大又黄又爽视频免费| 1024视频免费在线观看| 永久免费av网站大全| 亚洲,欧美,日韩| 久热爱精品视频在线9| 亚洲国产av影院在线观看| 久久久久久久大尺度免费视频| av福利片在线| 自线自在国产av| 国产成人91sexporn| 日日啪夜夜爽| av又黄又爽大尺度在线免费看| 久久狼人影院| 超碰成人久久| 黄频高清免费视频| 久久人人97超碰香蕉20202| av一本久久久久| 国产成人精品在线电影| 精品国产超薄肉色丝袜足j| 午夜久久久在线观看| 国产成人精品久久久久久| 狂野欧美激情性bbbbbb| 亚洲天堂av无毛| 免费黄频网站在线观看国产| 成年人午夜在线观看视频| 亚洲欧美一区二区三区国产| www.自偷自拍.com| 免费看不卡的av| 免费在线观看完整版高清| 美女中出高潮动态图| 久久久国产欧美日韩av| 超碰97精品在线观看| 午夜福利影视在线免费观看| 天天躁日日躁夜夜躁夜夜| 性少妇av在线| 国产精品 国内视频| 精品视频人人做人人爽| 最近的中文字幕免费完整| 丁香六月欧美| 欧美中文综合在线视频| 免费观看a级毛片全部| 夫妻性生交免费视频一级片| 亚洲男人天堂网一区| 新久久久久国产一级毛片| 日韩 亚洲 欧美在线| xxxhd国产人妻xxx| 久久久国产一区二区| 日日撸夜夜添| 99香蕉大伊视频| 你懂的网址亚洲精品在线观看| 久久久欧美国产精品| 狠狠精品人妻久久久久久综合| 一区二区三区乱码不卡18| 国产精品免费大片| 自拍欧美九色日韩亚洲蝌蚪91| 成人手机av| 考比视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 黄频高清免费视频| 国产欧美亚洲国产| 日韩视频在线欧美| 国产精品久久久久久久久免| 中文字幕精品免费在线观看视频| 亚洲av电影在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片| 婷婷色麻豆天堂久久| 19禁男女啪啪无遮挡网站| 成年av动漫网址| 日韩欧美精品免费久久| av电影中文网址| 久久ye,这里只有精品| xxx大片免费视频| 99精品久久久久人妻精品| 亚洲一级一片aⅴ在线观看| 免费在线观看黄色视频的| 满18在线观看网站| 纯流量卡能插随身wifi吗| 美女扒开内裤让男人捅视频| 1024香蕉在线观看| 99久久综合免费| 亚洲精品久久久久久婷婷小说| 亚洲伊人色综图| 中文乱码字字幕精品一区二区三区| 免费久久久久久久精品成人欧美视频| 国产女主播在线喷水免费视频网站| 精品国产一区二区三区四区第35| 亚洲国产欧美日韩在线播放| 午夜福利视频精品| 成年女人毛片免费观看观看9 | 欧美激情极品国产一区二区三区| 欧美黄色片欧美黄色片| kizo精华| 精品免费久久久久久久清纯 | 久久热在线av| 国产亚洲午夜精品一区二区久久| 高清av免费在线| xxxhd国产人妻xxx| 丝袜脚勾引网站| 国产亚洲一区二区精品| 精品亚洲成a人片在线观看| 免费少妇av软件| 少妇精品久久久久久久| 飞空精品影院首页| 国产xxxxx性猛交| 亚洲,欧美精品.| 久久女婷五月综合色啪小说| 天天操日日干夜夜撸| 天天躁夜夜躁狠狠久久av| av在线观看视频网站免费| 视频区图区小说| 夫妻午夜视频| 日日撸夜夜添| 菩萨蛮人人尽说江南好唐韦庄| 精品国产乱码久久久久久小说| 欧美人与性动交α欧美精品济南到| 欧美亚洲 丝袜 人妻 在线| 超色免费av| 又大又爽又粗| 亚洲av成人精品一二三区| 日本vs欧美在线观看视频| 又粗又硬又长又爽又黄的视频| 国产一区亚洲一区在线观看| 欧美中文综合在线视频| 多毛熟女@视频| 综合色丁香网| 色视频在线一区二区三区| 国产成人免费无遮挡视频| 不卡av一区二区三区| 中文字幕制服av| 一级爰片在线观看| 欧美黑人欧美精品刺激| 叶爱在线成人免费视频播放| 麻豆av在线久日| 日韩制服骚丝袜av| 国产成人91sexporn| 视频区图区小说| 成人国语在线视频| www.精华液| 91精品国产国语对白视频| 2021少妇久久久久久久久久久| 啦啦啦啦在线视频资源| 国产xxxxx性猛交| 久久鲁丝午夜福利片| 中文字幕人妻丝袜制服| 桃花免费在线播放| 国产免费视频播放在线视频| 亚洲国产欧美日韩在线播放| 老汉色av国产亚洲站长工具| 麻豆av在线久日| 中文精品一卡2卡3卡4更新| 黄色毛片三级朝国网站| 青春草国产在线视频| 久久狼人影院| 国产有黄有色有爽视频| 欧美日韩成人在线一区二区| 亚洲欧洲国产日韩| 亚洲欧洲精品一区二区精品久久久 | 在线观看免费视频网站a站| 久久精品久久久久久久性| 水蜜桃什么品种好| 亚洲欧美日韩另类电影网站| 看免费av毛片| 又粗又硬又长又爽又黄的视频| 一边亲一边摸免费视频| 黑人猛操日本美女一级片| 欧美另类一区| 国产亚洲午夜精品一区二区久久| 亚洲精品国产av成人精品| av线在线观看网站| 18禁国产床啪视频网站| 国产毛片在线视频| 亚洲精品美女久久久久99蜜臀 | 晚上一个人看的免费电影| 丁香六月天网| 在线精品无人区一区二区三| 男女下面插进去视频免费观看| 男人操女人黄网站| 十分钟在线观看高清视频www| 亚洲欧洲精品一区二区精品久久久 | 日韩av在线免费看完整版不卡| 亚洲国产精品999| 亚洲成人av在线免费| 一级,二级,三级黄色视频| 菩萨蛮人人尽说江南好唐韦庄| 大码成人一级视频| 一二三四中文在线观看免费高清| 好男人视频免费观看在线| 亚洲精品国产区一区二| 丁香六月天网| 欧美变态另类bdsm刘玥| 国产免费现黄频在线看| 欧美在线一区亚洲| 国产成人91sexporn| 亚洲伊人久久精品综合| 日本vs欧美在线观看视频| 久久人人爽av亚洲精品天堂| www日本在线高清视频| 纯流量卡能插随身wifi吗| 欧美变态另类bdsm刘玥| 丝袜美腿诱惑在线| 国产欧美日韩综合在线一区二区| 国产精品 欧美亚洲| av在线播放精品| 亚洲精品日本国产第一区| 久久久久久人人人人人| 少妇精品久久久久久久| 中文字幕人妻丝袜一区二区 | 国产精品亚洲av一区麻豆 | 人妻 亚洲 视频| 亚洲久久久国产精品| 天天躁夜夜躁狠狠久久av| 爱豆传媒免费全集在线观看| 天美传媒精品一区二区| 国产在视频线精品| 欧美精品av麻豆av| 69精品国产乱码久久久| 又粗又硬又长又爽又黄的视频| 欧美日韩国产mv在线观看视频| 水蜜桃什么品种好| 亚洲第一av免费看| 久久亚洲国产成人精品v| 制服丝袜香蕉在线| 亚洲精品久久久久久婷婷小说| 精品一区二区免费观看| av在线播放精品| 在线免费观看不下载黄p国产| 欧美精品高潮呻吟av久久| 啦啦啦 在线观看视频| 大香蕉久久成人网| 国产人伦9x9x在线观看| 99九九在线精品视频| 97人妻天天添夜夜摸| 乱人伦中国视频| 久久久欧美国产精品| 欧美xxⅹ黑人| 制服诱惑二区| 国产成人欧美在线观看 | 在线亚洲精品国产二区图片欧美| tube8黄色片| 免费观看av网站的网址| 少妇人妻久久综合中文| 老司机深夜福利视频在线观看 | 日韩中文字幕视频在线看片| 超碰成人久久| 秋霞伦理黄片| 国产精品偷伦视频观看了| 欧美变态另类bdsm刘玥| 欧美精品人与动牲交sv欧美| 午夜激情久久久久久久| 一边摸一边做爽爽视频免费| 高清视频免费观看一区二区| 91国产中文字幕| 一区二区三区四区激情视频| 97在线人人人人妻| 精品国产露脸久久av麻豆| 成人手机av| av网站在线播放免费| 久久精品久久精品一区二区三区| 9热在线视频观看99| 女人高潮潮喷娇喘18禁视频| 国产av精品麻豆| 久久性视频一级片| 十八禁人妻一区二区| 黑丝袜美女国产一区| 19禁男女啪啪无遮挡网站| 宅男免费午夜| av视频免费观看在线观看| 国产xxxxx性猛交| 久久狼人影院| 欧美av亚洲av综合av国产av | 精品福利永久在线观看| 中文字幕色久视频| 老司机亚洲免费影院| 久久精品国产a三级三级三级| 亚洲久久久国产精品| 最近最新中文字幕大全免费视频 | 人体艺术视频欧美日本| 美女脱内裤让男人舔精品视频| 国产欧美亚洲国产| 欧美黑人精品巨大| 18禁动态无遮挡网站| 欧美成人精品欧美一级黄| 韩国高清视频一区二区三区| 久久久久久人人人人人| 国产爽快片一区二区三区| 一本大道久久a久久精品| 欧美人与性动交α欧美软件| 777久久人妻少妇嫩草av网站| 亚洲国产精品成人久久小说| 国产乱人偷精品视频| 人人妻人人爽人人添夜夜欢视频| 国产男女超爽视频在线观看| 电影成人av| 午夜老司机福利片| 国产视频首页在线观看| 国产精品欧美亚洲77777| 七月丁香在线播放| √禁漫天堂资源中文www| 哪个播放器可以免费观看大片| 亚洲精品久久久久久婷婷小说| 天天躁夜夜躁狠狠久久av| 精品免费久久久久久久清纯 | 欧美精品人与动牲交sv欧美| 国产精品 国内视频| 99热国产这里只有精品6| 久久久久精品国产欧美久久久 | 咕卡用的链子| 久久久久国产精品人妻一区二区| 少妇精品久久久久久久| 久热这里只有精品99| 蜜桃在线观看..| 国产成人精品福利久久| 女人久久www免费人成看片| 久久狼人影院| 久久婷婷青草| 亚洲男人天堂网一区| 国产成人一区二区在线| 亚洲成人av在线免费| 亚洲精华国产精华液的使用体验| 尾随美女入室| 狠狠婷婷综合久久久久久88av| 91aial.com中文字幕在线观看| 日韩中文字幕视频在线看片| 午夜日韩欧美国产| 悠悠久久av| 另类亚洲欧美激情| 国产av码专区亚洲av| 天堂中文最新版在线下载| 91国产中文字幕| 777米奇影视久久| 欧美中文综合在线视频| 精品一区二区三区av网在线观看 | 日本欧美视频一区| 人体艺术视频欧美日本| 日韩欧美一区视频在线观看| 欧美日韩亚洲综合一区二区三区_| 国产精品秋霞免费鲁丝片| 日本av手机在线免费观看| 建设人人有责人人尽责人人享有的| 亚洲欧美色中文字幕在线| av网站免费在线观看视频| 人人澡人人妻人| 免费黄网站久久成人精品| 国产精品二区激情视频| 亚洲一码二码三码区别大吗| 日韩大码丰满熟妇| 啦啦啦在线观看免费高清www| av天堂久久9| 国产片特级美女逼逼视频| 久久久久久久大尺度免费视频| a级片在线免费高清观看视频| 免费观看av网站的网址| 91国产中文字幕| 精品午夜福利在线看| 亚洲综合精品二区| 悠悠久久av| 80岁老熟妇乱子伦牲交| 母亲3免费完整高清在线观看| √禁漫天堂资源中文www| 日韩欧美精品免费久久| 乱人伦中国视频| 成人亚洲精品一区在线观看| 国产 一区精品| 日本av手机在线免费观看| 九色亚洲精品在线播放| 久久婷婷青草| 男女免费视频国产| 一本—道久久a久久精品蜜桃钙片| 好男人视频免费观看在线| 国产精品久久久久久精品古装| 精品一品国产午夜福利视频| 精品免费久久久久久久清纯 | 亚洲一卡2卡3卡4卡5卡精品中文| 欧美久久黑人一区二区| 婷婷成人精品国产| 欧美人与善性xxx| 男人操女人黄网站| 男女高潮啪啪啪动态图| 一边摸一边抽搐一进一出视频| 久久精品国产a三级三级三级| 女性生殖器流出的白浆| 成人黄色视频免费在线看| 亚洲欧洲精品一区二区精品久久久 | 久久久亚洲精品成人影院| 亚洲欧洲日产国产| 免费黄色在线免费观看| 国产午夜精品一二区理论片| 麻豆av在线久日| 久久热在线av| 国产精品久久久久久久久免| 亚洲第一青青草原| 91国产中文字幕| 久久久久久久久免费视频了| 蜜桃国产av成人99| 久久婷婷青草| 欧美亚洲 丝袜 人妻 在线| 久久ye,这里只有精品| 一边摸一边做爽爽视频免费| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 亚洲欧美色中文字幕在线| 国产精品久久久av美女十八| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站| 欧美激情 高清一区二区三区| 成人国产av品久久久| 免费在线观看黄色视频的| 久久久久网色| 国产在线视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成人一二三区av| 制服丝袜香蕉在线| 国产99久久九九免费精品| 另类亚洲欧美激情| 七月丁香在线播放| 成年av动漫网址| 国产精品国产三级专区第一集| 午夜91福利影院| 久久久欧美国产精品| 一级毛片 在线播放| 婷婷成人精品国产| 制服丝袜香蕉在线| 久久97久久精品| 日韩一区二区视频免费看| 午夜老司机福利片| 日韩一区二区三区影片| 满18在线观看网站| 国语对白做爰xxxⅹ性视频网站| 黄色视频不卡| 国产男人的电影天堂91| 新久久久久国产一级毛片| 日本一区二区免费在线视频| 老司机在亚洲福利影院| 国产av精品麻豆| xxxhd国产人妻xxx| 亚洲精品一二三| 高清av免费在线| 精品人妻一区二区三区麻豆| 午夜老司机福利片| 天堂俺去俺来也www色官网| 一级,二级,三级黄色视频| 一级毛片我不卡| 精品酒店卫生间| 男女高潮啪啪啪动态图| 亚洲国产毛片av蜜桃av| 操出白浆在线播放| 三上悠亚av全集在线观看| 久久性视频一级片| 久久 成人 亚洲| 不卡视频在线观看欧美| 丰满乱子伦码专区| 久久精品熟女亚洲av麻豆精品| 性少妇av在线| 乱人伦中国视频| 亚洲av成人精品一二三区| 亚洲欧美精品综合一区二区三区| 日韩,欧美,国产一区二区三区| 一边亲一边摸免费视频| 波多野结衣av一区二区av| 国产黄频视频在线观看| 老司机在亚洲福利影院| 国产精品一区二区在线观看99| 日本猛色少妇xxxxx猛交久久| 悠悠久久av| 亚洲精品视频女| 2021少妇久久久久久久久久久| 黄色一级大片看看| 国产精品女同一区二区软件| 久热这里只有精品99| 精品一区二区三区四区五区乱码 | 色婷婷av一区二区三区视频| 国产精品蜜桃在线观看| 免费人妻精品一区二区三区视频| 亚洲欧美精品综合一区二区三区| 国产福利在线免费观看视频| 99久久99久久久精品蜜桃| 国产在线视频一区二区| svipshipincom国产片| 日本91视频免费播放| 午夜老司机福利片| 老司机靠b影院| 一区二区三区精品91| 无限看片的www在线观看| 十八禁网站网址无遮挡| 国产精品欧美亚洲77777| 日日啪夜夜爽| 国产黄频视频在线观看| 成年人午夜在线观看视频| 五月开心婷婷网| 9色porny在线观看| 国产精品国产av在线观看| 国产高清国产精品国产三级| 国产一卡二卡三卡精品 | 99久久综合免费| 九草在线视频观看| 丰满少妇做爰视频| 国产福利在线免费观看视频| 国产黄色视频一区二区在线观看| 亚洲激情五月婷婷啪啪| 又黄又粗又硬又大视频| 黑丝袜美女国产一区| 午夜福利视频精品| 成人免费观看视频高清| 欧美日韩视频高清一区二区三区二| 精品人妻在线不人妻| 国产精品 欧美亚洲| 精品国产一区二区三区四区第35| 一级毛片 在线播放| 97在线人人人人妻| e午夜精品久久久久久久| 国产野战对白在线观看| 性少妇av在线| 成年av动漫网址| 亚洲欧洲国产日韩| 国产成人欧美| 国产精品无大码| 超碰97精品在线观看| 大香蕉久久成人网| av国产久精品久网站免费入址| 国产男女超爽视频在线观看| 国产在视频线精品| 久久99一区二区三区| 老司机在亚洲福利影院| 母亲3免费完整高清在线观看| 一级a爱视频在线免费观看| 国产精品国产三级专区第一集| 日韩中文字幕欧美一区二区 | 男女免费视频国产| av在线播放精品| 亚洲成色77777| 久久精品国产综合久久久| 国产高清不卡午夜福利| 欧美人与性动交α欧美软件| 满18在线观看网站| 777久久人妻少妇嫩草av网站| 免费女性裸体啪啪无遮挡网站| 99九九在线精品视频| 国产精品一区二区在线观看99| 日日摸夜夜添夜夜爱| 亚洲国产日韩一区二区| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 制服人妻中文乱码| 亚洲精品一二三| 成人漫画全彩无遮挡| 亚洲视频免费观看视频| www.熟女人妻精品国产| 一区二区日韩欧美中文字幕| 日本猛色少妇xxxxx猛交久久| 国产xxxxx性猛交| 成年女人毛片免费观看观看9 | 亚洲精品美女久久久久99蜜臀 | 久久热在线av| 十八禁网站网址无遮挡| 日本av手机在线免费观看| 成年人免费黄色播放视频| 日韩成人av中文字幕在线观看| 老汉色av国产亚洲站长工具| 国精品久久久久久国模美| 成人亚洲精品一区在线观看| 一级毛片我不卡| 欧美在线一区亚洲| 国产日韩欧美视频二区| 亚洲欧美成人综合另类久久久|