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

    近年來厄爾尼諾期間北赤道流輸運的年際變化

    2016-07-06 10:30:13吳國麗王慶業(yè)翟方國胡敦欣中國科學院海洋研究所海洋環(huán)流與波動重點實驗室山東青島26607中國科學院大學北京0009中國石油大學華東理學院山東青島266580中國海洋大學海洋環(huán)境學院山東青島26600
    海洋科學 2016年2期
    關鍵詞:厄爾尼諾氣旋赤道

    吳國麗, 王慶業(yè), 翟方國, 胡敦欣(. 中國科學院海洋研究所 海洋環(huán)流與波動重點實驗室, 山東 青島 26607; 2. 中國科學院大學, 北京0009; 3. 中國石油大學(華東) 理學院, 山東 青島266580; . 中國海洋大學 海洋環(huán)境學院, 山東 青島26600)

    ?

    近年來厄爾尼諾期間北赤道流輸運的年際變化

    吳國麗1, 2, 3, 王慶業(yè)1, 翟方國4, 胡敦欣1
    (1. 中國科學院海洋研究所 海洋環(huán)流與波動重點實驗室, 山東 青島 266071; 2. 中國科學院大學, 北京100049; 3. 中國石油大學(華東) 理學院, 山東 青島266580; 4. 中國海洋大學 海洋環(huán)境學院, 山東 青島266100)

    為了研究近年來厄爾尼諾期間北赤道流輸運的年際變化, 本文利用海洋客觀分析數(shù)據(jù)MOAA GPV(Grid Point Value of the Monthly Objective Analysis)以及P-vector方法計算了北太平洋絕對地轉流,探討了2001~2013年期間厄爾尼諾與北赤道流輸運之間的關系。在此期間發(fā)生的4次厄爾尼諾事件中,北赤道流輸運在2002~2003、2006~2007、2009~2010年的厄爾尼諾成熟期都出現(xiàn)了明顯的增強, 但是在2004~2005年的厄爾尼諾成熟期并沒有明顯的增強。進一步分析發(fā)現(xiàn), 在2002~2003年、2006~2007年、2009~2010年的厄爾尼諾成熟期, 10°N以南的熱帶西北太平洋區(qū)域出現(xiàn)了負的海面高度異常和氣旋式環(huán)流異常, 這主要是由熱帶環(huán)流區(qū)域出現(xiàn)的西風異常和正的Ekman抽吸通過Rossby波西傳到熱帶西太平洋區(qū)域所致; 但是在2004~2005年厄爾尼諾成熟期, 海面溫度異常的分布明顯不同, 西風異常和正的Ekman抽吸異常明顯北移, 導致負的海面高度異常和氣旋式環(huán)流異常出現(xiàn)在了10°N以北的西北太平洋區(qū)域, 使得北赤道流輸運在2004~2005年的厄爾尼諾成熟期沒有明顯的增強。

    北赤道流輸運; 厄爾尼諾; 年際變異

    [Foundation: “Global Change and Air-Sea Interaction” Project, No. GASI-03-01-01-05; National Basic Research Program of China, No. 2012CB417401, No. 2013CB956202]

    在熱帶北太平洋, 北赤道流(North Equatorial Current, NEC)自東向西流動, 到達菲律賓海岸附近分叉形成向北流動的黑潮(Kuroshio Current, KC)和向南流動的棉蘭老流(Mindanao Current, MC), 組成了NEC-MC-KC(NMK)海流系統(tǒng)[1-3]。NEC作為NMK海流系統(tǒng)的核心, 它對西太平洋暖池的熱量收支[4-5]、全球的熱鹽環(huán)流[6]等都具有非常重要的作用, 幾十年以來, 受到世界各國物理海洋學家們廣泛的關注和研究。前人通過對觀測數(shù)據(jù)、海洋同化數(shù)據(jù)的分析以及數(shù)值模式的模擬等重點研究了NEC的分叉和輸運等問題, 其中對NEC輸運在時間變異上的研究主要集中在季節(jié)和年際尺度上[4, 7-15]。已有的研究結果表明, NEC輸運的年際變化與厄爾尼諾-南方濤動(ENSO)事件密切相關, 在ENSO暖事件(厄爾尼諾)時期, NEC輸運增強, 在ENSO冷事件(拉尼娜)時期, NEC輸運減弱[4, 13-16]。但是, Qiu等[8]以及Zhai等[14]分別利用沿137°E觀測數(shù)據(jù)和衛(wèi)星高度計數(shù)據(jù)的研究結果指出, 并不是所有的NEC輸運的極值都與ENSO事件有關。因此, NEC輸運與ENSO之間存在著復雜的關系。近年來, 主要基于Argo浮標資料構建的三維客觀分析網(wǎng)格數(shù)據(jù)MOAA GPV(Grid Point Value of the Monthly Objective Analysis)被廣泛地用來研究全球變化問題[17]以及一些復雜的區(qū)域問題[18-19]。因此, 本文將基于利用MOAA GPV數(shù)據(jù)計算的北太平洋絕對地轉流來檢驗2001~2013年厄爾尼諾與NEC輸運之間的關系, 并討論其動力學過程。

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

    1.1 數(shù)據(jù)

    本文主要用到的MOAA GPV數(shù)據(jù)提供了2001 年1月~2013年12月月平均的溫度和鹽度資料[20], 其水平分辨率為1°×1°, 垂直方向從10 dbar到2 000 dbar分為25層。MOAA GPV數(shù)據(jù)是主要利用Argo浮標資料, 同時結合TRITON(Triangle Trans-OceanBuoy-Netword)浮標資料以及現(xiàn)有的CTD(conductivity-tem-perature-depth)資料通過客觀分析方法所構建。關于構建過程中所涉及的數(shù)據(jù)處理過程、網(wǎng)格化技術和誤差分析等方面的詳細描述可參考文獻[20]。

    本文用到的衛(wèi)星高度計數(shù)據(jù)(由網(wǎng)站http: //www. aviso.oceanobs.com/下載)包括衛(wèi)星觀測的海面高度以及由地轉關系得到的海表地轉流。其原始數(shù)據(jù)為周數(shù)據(jù), 為了分析變量的低頻變化, 周數(shù)據(jù)被平均到了月數(shù)據(jù)。另外, 本文用到的月平均的風場數(shù)據(jù)來自ECMWF(European Centre for Medium-Range Weather Forecasts)海洋分析/再分析系統(tǒng)3(ECMWF ORA-S3)[21]。

    1.2 P-vector方法

    本文利用P-vector方法計算北太平洋絕對地轉流。P-vector方法是在地轉平衡和Boussinesq近似下基于位勢密度和位勢渦度守恒的一種計算絕對地轉流的方法[22-23]。其中地轉流的方向由等位勢面與等密度面交叉的方向決定, 地轉流的方向也被稱為P-vector, 地轉流的大小由任意兩層之間的熱成風關系來計算。以往的研究表明P-vector方法能夠較好描述大洋和邊緣海中海洋環(huán)流的基本特征[22, 24-25]。

    由于在上混合層海洋運動不能精確的遵守位勢密度和位勢渦度守恒, 因此, 本文參照Zhang等[26]和Yuan等[27]的做法, 僅僅對大洋中層的地轉流使用P-vector方法重構, 即800~2 000 dbar之間的絕對地轉流使用P-vector方法計算, 然后用800 dbar層的絕對地轉流作為參考面, 通過動力計算方法獲得800 dbar以上的絕對地轉流。已有的研究結果表明,如此計算的絕對地轉流能夠很好地描述位于北太平洋的NEC等海洋環(huán)流的結構特征及動力特征[26-27]。

    1.3 2001~2013年的厄爾尼諾事件

    根據(jù)CPC(Climate Prediction Center, http: //www. cpc.ncep.noaa.gov/data/indices)Ni?o-3.4指數(shù), 本文定義Ni?o-3.4指數(shù)≥0.5°C并持續(xù)5個月的事件為一個厄爾尼諾事件。2001年1月~2013年12月一共發(fā)生了4次厄爾尼諾事件(分別記為2002~2003、2004~ 2005、2006~2007和2009~2010年厄爾尼諾事件), 其中在2004~2005年厄爾尼諾期間, 海面溫度(Sea Surface Temperature, SST)正異常的極值發(fā)生在2004年的9月, 其它厄爾尼諾期間的SST正異常的極值分別發(fā)生在2002年的11月、2006年的12月和2009年的12月(圖1a)。

    圖1 Ni?o-3.4 SST異常(a)和NEC分別基于MOAA GPV數(shù)據(jù)(b)以及衛(wèi)星高度計數(shù)據(jù)(c)沿130°~180°E輸運的時間序列Fig. 1 Ni?o-3.4 SST anomaly (a) and the NEC transport at 130°~180° E from MOAA GPV (b) and altimeter products (c), respectively

    1.4 月平均數(shù)據(jù)年際信號的提取

    本文利用3階Butterworth濾波器對月平均變量進行13個月的低通濾波去除高頻信號, 得到月平均數(shù)據(jù)的年際信號。為了獲取月平均數(shù)據(jù)的年際異常信號, 在對月平均變量進行13個月的低通濾波后,再減去多年平均。

    2 主要結果

    2.1 NEC輸運的年際變化

    根據(jù)MOAA GPV數(shù)據(jù)中137°E斷面上的年平均緯向速度圖(圖2), NEC沿137°E的主體位于8°~ 18°N、26.7σθ等密度面以上, 在其它的經(jīng)度斷面上也有類似的分布特征。因此, 本文定義NEC沿某一經(jīng)度的輸運為8°~18°N、26.7σθ等密度面以上的西向流速的積分, 其中西向流速來自利用MOAA GPV數(shù)據(jù)通過P-vector方法計算的北太平洋絕對地轉流。圖1a是Ni?o-3.4 SST異常的時間序列圖, 其中黃色帶代表4次厄爾尼諾事件, 圖1b是NEC分別沿130°~ 180°E輸運的時間序列圖, 其中圖1a、圖1b中的紅色橫直線對應的是4次厄爾尼諾期間SST正異常最大的月份。由圖1b可以看出, 在2002~2003、2006~ 2007、2009~2010年厄爾尼諾期間SST正異常最大的月份, NEC輸運都出現(xiàn)了明顯的增強, 但是, 在2004~2005年厄爾尼諾期間SST正異常最大的月份, NEC輸運并沒有出現(xiàn)明顯的增強。為進一步增加這一現(xiàn)象的可信性, 本文利用衛(wèi)星高度計從2001年1月~ 2013年8月的數(shù)據(jù)資料計算了NEC沿130°~180°E的海表輸運(圖1c)。由圖1c可以看出, 厄爾尼諾期間SST正異常與NEC輸運之間的對應關系和利用MOAA GPV數(shù)據(jù)得到的結果非常吻合。因此, 以上描述的NEC輸運與4次厄爾尼諾之間的對應關系是可信的。

    圖2 MOAA GPV數(shù)據(jù)中沿137°E斷面的年平均緯向西向速度(cm/s)Fig. 2 Annual mean zonal westward velocities (cm/s) across 137°E from the MOAA GPV

    為了更清楚地研究厄爾尼諾和NEC輸運之間的關系, 圖3給出Ni?o-3.4 SST異常和基于MOAA GPV數(shù)據(jù)計算的140°E斷面上NEC輸運異常的時間序列, 為方便比較, 這里取NEC的西向輸運為正值。由圖可見, 自2001年以來的4次厄爾尼諾事件中, NEC輸運在2002~2003、2006~2007、2009~2010年厄爾尼諾成熟期都出現(xiàn)了西向的異常極值, 但是, 在2004~2005年厄爾尼諾成熟期, NEC輸運并沒有出現(xiàn)西向的異常極值, 而是出現(xiàn)了東向異常。另外, 沿140°E斷面, NEC輸運異常超前Ni?o-3.4 SST異常大約2個月(相關系數(shù)為0.59, 高于95%置信水平), 這與Zhai等[15]得出的NEC沿138°~156°E的輸運異常超前Ni?o-3.4指數(shù)1~2個月的結果一致。

    圖3 Ni?o-3.4 SST異常和NEC沿140°E輸運異常的時間序列Fig. 3 Ni?o-3.4 SST and NEC transport anomalies at 140°E

    2.2 流場的演變

    由于NEC輸運的年際變異與北太平洋海洋環(huán)流的變化有關, 因此利用合成分析法研究2002~2003、2006~2007、2009~2010年以及2004~2005年厄爾尼諾期間熱帶北太平洋環(huán)流異常的演變過程。圖4是2002~2003、2006~2007、2009~2010年厄爾尼諾期間熱帶北太平洋上200 m平均的流異常和海面高度異常的演變合成圖, 其中T=0表示每次厄爾尼諾期間Ni?o-3.4 SST正異常最大的月份, 負值表示T=0之前的月份, 正值表示T=0之后的月份(以下類似)。由圖可見, 從T= –6個月開始, 在10°N以南、赤道以北的熱帶西太平洋區(qū)域, 出現(xiàn)了負的海面高度異常和氣旋式環(huán)流異常。隨著時間的演變, 負的海面高度異常和氣旋式環(huán)流異常逐漸加強, 在T=0個月的時候達到極值, 之后逐漸的減弱。在西北太平洋, 氣旋式環(huán)流異常將會加強NEC的輸運。

    為了進一步研究熱帶西北太平洋環(huán)流變化的垂直分布, 圖5給出的是2002~2003、2006~2007、2009~2010年厄爾尼諾期間140°E斷面上緯向流異常的合成圖。由圖可以看出, T= –8個月的時候, 在大約12°N的南北兩側分別出現(xiàn)了西/東向流異常。隨著時間的演變, 西向流異常逐漸地加強并向北移動, 到T=2個月的時候, 在10°N以南完全變成了東向流異常, 而西向流異常的核心移動到大約18°N附近。因此, NEC輸運在2002~2003、2006~2007、2009~2010年厄爾尼諾的成熟期達到極大值, 然后逐漸減弱。

    圖6是2004~2005年厄爾尼諾期間熱帶北太平洋上200 m平均的流異常和海面高度異常圖。由圖可見, 在厄爾尼諾成熟期之前, 從T= –8個月到T= –2個月, 在10°N以南的熱帶西太平洋區(qū)域出現(xiàn)了比較強的負的海面高度異常和氣旋式的環(huán)流異常,由此導致在2004~2005年厄爾尼諾成熟期之前出現(xiàn)了NEC輸運增強的現(xiàn)象。但是, 從T=0個月開始, 負的海面高度異常和氣旋式的環(huán)流異常逐漸擴展到10°N以北的區(qū)域。

    由2004~2005年厄爾尼諾期間沿140°E斷面的緯向流異常圖(圖7)可以看出, 在T= –8個月的時候, 10°N的南北兩側分別出現(xiàn)了較強的東/西向流異常,這些緯向流異常構成了圖6中氣旋式環(huán)流異常的一部分, 并且隨著時間的演變逐漸減弱。從T= –2個月開始, 西向流異常的核心移動到15°N以北, 而東向流異常逐漸移動到10°N以北, 并逐漸加強。因此, 在2004~2005年厄爾尼諾的成熟期并沒有出現(xiàn)NEC輸運的加強, 反而出現(xiàn)了NEC輸運的東向異常。

    2.3 風場的影響

    已有的研究結果表明, 在熱帶西北太平洋, 海洋環(huán)流的低頻變異主要由海面風強迫所導致[4, 8, 11-14]。因此, 圖8給出的是2002~2003、2006~2007、2009~2010年厄爾尼諾期間風應力異常和相應的Ekman抽吸速度異常的合成圖。由圖可見, 從T= –6個月到T= –2個月, 在赤道以北、20°N以南的熱帶北太平洋的中西部都出現(xiàn)了比較強的西風異常和正的Ekman抽吸速度異常, 之后逐漸減弱。正的Ekman抽吸速度異常將會引起西向傳播的上升Rossby波, 導致在熱帶西北太平洋出現(xiàn)負的海面高度異常和氣旋式環(huán)流異常。

    圖9是2004~2005年厄爾尼諾期間風應力異常和相應的Ekman抽吸速度異常的時間演變圖。由圖可以看出, 從T= –8個月到T= –6個月, 在15°N以南的熱帶北太平洋的中西部也出現(xiàn)了較強的正的Ekman抽吸速度異常。但是, 隨著時間的演變, 正的Ekman抽吸速度異常逐漸的向北移動至大約20°N附近。因此, 在2004~2005年厄爾尼諾成熟期, 氣旋式環(huán)流異常的中心比2002~2003、2006~2007、2009~ 2010年厄爾尼諾成熟期時偏北。

    圖5 2002~2003、2006~2007、2009~2010年厄爾尼諾期間沿140°E斷面的緯向流速度異常的合成Fig. 5 Composite anomalies of the zonal velocity at 140° E during the 2002~2003, 2006~2007, 2009~2010 El Ni?o events

    圖6 2004~2005年厄爾尼諾期間上200 m平均的流異常和海面高度異常Fig. 6 Anomalies of the upper 200 m mean flow and sea surface height during the 2004~2005 El Ni?o event

    圖7 2004~2005年厄爾尼諾期間沿140°E斷面的緯向流速度異常Fig. 7 Anomalies of the zonal velocity at 140° E during the 2004~2005 El Ni?o event

    圖9 2004~2005年厄爾尼諾期間風應力異常和相應的Ekman抽吸速度異常Fig. 9 Anomalies of the wind stress vector and associated Ekman pumping velocity during the 2004~2005 El Ni?o event

    在熱帶北太平洋, 海面風場和Ekman抽吸與海面溫度分布有關, 因此, 圖10給出了2002~ 2003、2006~2007、2009~2010和2004~2005年厄爾尼諾成熟期(厄爾尼諾期間SST正異常最大的月份以及之前和之后的一個月)平均的SST異常圖。由圖可見, 在2002~2003、2006~2007、2009~2010年厄爾尼諾成熟期, 太平洋中東部的赤道區(qū)域出現(xiàn)了較強的SST正異常; 而在2004~2005年厄爾尼諾成熟期, 太平洋中部的赤道區(qū)域雖然也出現(xiàn)了SST正異常, 但強度較弱, 大范圍的正的SST異常出現(xiàn)在北太平洋區(qū)域, 特別是在20°N以北出現(xiàn)了較強的SST正異常。因此2004~2005年厄爾尼諾成熟期的Ekman抽吸速度正異常相對于其它3次厄爾尼諾成熟期偏北。

    另外值得關注的是, Wu[28]指出ENSO對西北太平洋的影響與北太平洋年代際濤動(PDO)的位相有關, 在PDO的暖位相時, ENSO和季風的年際關系比較弱。2004~2005的厄爾尼諾事件恰好發(fā)生在新世紀的一次暖的PDO時期, 而暖的PDO位相會影響東亞夏季風, 進而使得ENSO對風場的調(diào)制作用變?nèi)? 導致2004~2005年厄爾尼諾時期的風場相對于其它3次厄爾尼諾時期偏北。

    圖10 2002~2003、2006~2007、2009~2010以及2004~2005年厄爾尼諾成熟期平均的SST異常Fig. 10 Anomalies of the mean SST in mature phases during the 2002~2003, 2006~2007, 2009~2010, and 2004~2005 El Ni?o events

    3 結論

    已有的研究結果普遍認為NEC輸運的年際變化與ENSO事件密切相關, 在ENSO暖事件時期, NEC輸運增強[4, 13-16], 但是, 也有研究結果表明, 并不是所有的NEC輸運的極值都和ENSO事件有關[8, 14],因此, NEC輸運與ENSO之間存在著復雜的關系。本文基于利用海洋客觀分析數(shù)據(jù)MOAA GPV通過P-vector方法計算出的北太平洋絕對地轉流, 研究了2001~2013年期間厄爾尼諾與NEC輸運之間的關系。結果顯示, 2001年以來的4次厄爾尼諾事件中, NEC輸運在2002~2003、2006~2007、2009~2010年的厄爾尼諾成熟期都出現(xiàn)了明顯的增強, 但是在2004~ 2005年厄爾尼諾成熟期并沒有明顯的增強。通過分析熱帶北太平洋的環(huán)流異常和海面高度異常發(fā)現(xiàn), 在2002~2003、2006~2007、2009~2010年的厄爾尼諾成熟期, 10°N以南的熱帶西北太平洋區(qū)域出現(xiàn)了較強的負的海面高度異常和氣旋式環(huán)流異常, 而在2004~2005年厄爾尼諾成熟期, 負的海面高度異常和氣旋式環(huán)流異常出現(xiàn)在了10°N以北的西北太平洋區(qū)域。進一步分析北太平洋的風場發(fā)現(xiàn), 在2002~ 2003、2006~2007、2009~2010年的厄爾尼諾成熟期之前, 赤道以北、20°N以南的熱帶北太平洋的中西部出現(xiàn)了比較強的西風異常和正的Ekman抽吸速度異常, 之后逐漸的減弱。正的Ekman抽吸速度異常將引起西向傳播的上升Rossby波, 導致在熱帶西北太平洋出現(xiàn)負的海面高度異常和氣旋式環(huán)流異常,而氣旋式的環(huán)流異常有利于加強北赤道流的輸運。在2004~2005年厄爾尼諾的成熟期之前, 15°N以南的熱帶北太平洋中西部也出現(xiàn)了較強的正的Ekman抽吸速度異常, 但是之后, 正的Ekman抽吸速度異常明顯北移。因此, 在2004~2005年厄爾尼諾成熟期,氣旋式環(huán)流異常的中心比2002~2003、2006~2007、2009~2010年厄爾尼諾成熟期時偏北。最后, 由熱帶太平洋平均的SST異常分布圖發(fā)現(xiàn), 在2004~2005年厄爾尼諾成熟期, 北太平洋出現(xiàn)了大范圍的正的SST異常, 特別是在20°N以北的區(qū)域強度較大, 而在赤道區(qū)域強度較弱, 因此, 2004~2005年厄爾尼諾時期的Ekman抽吸速度正異常相對于其他三次厄爾尼諾時期偏北。另外, 暖的PDO位相也使得在2004~2005年厄爾尼諾時期, ENSO對風場的調(diào)制能力變?nèi)?。在未來的工作? 我們將繼續(xù)關注PDO和NEC輸運在年際尺度上的關系。

    [1] Toole J, Millard R, Wang Z, et al. Observations of the Pacific North Equatorial bifurcation at the Philippinecoast[J]. J Phys Oceanogr, 1990, 20(2): 307-318.

    [2] Hu D X, Cui M C. The western boundary current of the Pacific and its role in the climate[J]. Chin J Oceanol Limnol, 1991, 9(1): 1-14.

    [3] Yaremchuk M, Qu T D. Seasonal Variability of the Large-Scale Currents near the Coast of the Philippines[J]. J Phys Oceanogr, 2004, 34(4): 844-855.

    [4] Qiu B, Lukas R. Seasonal and interannual variability of the North Equatorial Current, the Mindanao Current and the Kuroshio along the Pacific western boundary[J]. J Geophys Res, 1996, 101(C5): 12315-12330.

    [5] Qu T D, Meyers G, Godfrey J S, et al. Upper ocean dynamics and its role in maintaining the annual mean western Pacific warm pool in a global GCM[J]. Int J Climatol, 1997, 17: 711-724.

    [6] Gordon A. Interocean exchanges of thermocline water[J]. J Geophys Res, 1986, 91: 5037-5046.

    [7] White W, Hasunuma K. Interannual variability in the baroclinic gyre structure of the western North Pacificfrom 1954–1974[J]. J Mar Res, 1980, 38: 651-672.

    [8] Qiu B, Joyce T. Interannual variability in the mid- and low-latitude western North Pacific[J]. J Phys Oceanogr, 1992, 22: 1062-1079.

    [9] Gu Y H. Zonal volume transports at 137oE in the northwestern Pacific[J]. Chin J Oceanol Limnol, 1996, 27(1): 79-85.

    [10] Wang F, Chang P, Hu D X, et al. Circulation in the western tropical Pacific Ocean and its seasonalvariation[J]. Chin Sci Bull, 2002, 47(7): 591-595.

    [11] Qu T D, Mitsudera H, Yamagata T. On the western boundary currents in the Philippine Sea[J]. J Geophys Res, 1998, 103(C4): 7537-7548.

    [12] Kashino Y, Espa?a N, Syamsudin F, et al. Observations of the North Equatorial Current, Mindanao Current, and Kuroshio Current System during the 2006/07 El Ni?o and 2007/08 La Ni?a[J]. J Oceanogr, 2009, 65(3): 325-333.

    [13] Kim Y, Qu T D, Jensen T, et al. Seasonal and interannual variations of the North Equatorial Current bifurcation in a high-resolution OGCM[J]. J Geophys Res, 2004, 109(C03040): 1-19.

    [14] Zhai F G, Hu D X. Interannual variability of transport and bifurcation of the North Equatorial Current in the tropical North Pacific Ocean[J]. Chin J Oceanol Limnol, 2012, 30(1): 177-185.

    [15] Zhai F G, Hu D X. Revisit the interannual variability of the North Equatorial Current transport with ECMWF ORA-S3[J]. J Geophys Res, 2013, 118: 1349-1366.

    [16] Wang Q Y, Hu D X. Bifurcation of the North Equatorial Current derived from altimetry in the Pacific Ocean[J]. J Hydrodyn, 2006, 18B: 620-626.

    [17] Hosoda S, Suga T, Shikama N, et al. Global surface layer salinity change detected by Argo and its implication for hydrological cycle intensification[J]. J Oceanogr, 2009, 65: 579-586.

    [18] Kashino Y, Ishida A, Hosoda S. Observed ocean variability in the Mindanao Dome region[J]. J Phys Oceanogr, 2011, 41: 287-302.

    [19] Kobashi F, Xie S. Interannual variability of the North Pacific Subtropical Countercurrent: Role of local ocean–atmosphere interaction[J]. J Oceanogr, 2012, 68: 113-126.

    [20] Hosoda S, Ohira T, Nakamura T. A monthly mean dataset of global oceanic temperature and salinity derived from Argo float observations[J]. Jam Rep Res Dev, 2008, 8: 47-59.

    [21] Balmaseda M, Vidard A, Anderson D. The ECMWF Ocean Analysis System: ORA-S3[J]. Mon Wea Rev, 2008, 136: 3018-3034.

    [22] Chu P C. P-vector method for determining absolute velocity from hydrographic data[J]. Marine Tech Soc J, 1995, 29(2): 3-14.

    [23] Chu P C. P-vector spiral and determination of absolute velocities[J]. J Oceanogr, 2000, 56: 591-599.

    [24] Chu P C, Fan C, Lozano C, et al. An airborne expandable bathythermograph (AXBT) survey of the South China Sea, May 1995[J]. J Geophys Res, 1998, 103(10): 21637-21652.

    [25] Chu P C, Lan J, Fan C. Japan Sea circulation and thermohaline structure Part I, Climatology[J]. J Phys Oceanogr, 2001, 31: 244-271.

    [26] Zhang Z C, Yuan D L, Chu, P C. Meridional Volume Transport in the Tropical Northwest Pacific Ocean Based on Argo Profiles[J]. Chin J Oceanol Limnol, 2013, 31(3): 656-664.

    [27] Yuan D L, Zhang Z C, Chu P C, et al. GeostrophicCirculation in the Tropical North Pacific Ocean Based on Argo Profiles[J]. J Phys Oceanogr, 2014, 44(2): 558-575.

    [28] Wu C R. Interannual modulation of the Pacific Decadal Oscillation(PDO) on the low-latitude western North Pacific[J]. P Oceanogr, 2013, 110: 49-58.

    (本文編輯: 李曉燕)

    Interannual variations of the North Equatorial Current transport in recent years during El Ni?o events

    WU Guo-li1, 2, 3, WANG Qing-ye1, ZHAI Fang-guo4, HU Dun-xin1
    (1. Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China; 3. College of Science, China University of Petroleum, Qingdao 266580, China; 4. College of Physical and Environmental Oceanography, Ocean University of China, Qingdao 266100, China)

    Jan., 30, 2015

    the North Equatorial Current transport; El Ni?o; interannual variation

    In this study, we calculated the absolute geostrophic currents in the North Pacific Ocean using the P-vector method and ocean Grid Point Value of the Monthly Objective Analysis using Argo data (MOAA GPV). We then investigated the relationship between the North Equatorial Current (NEC) transport and El Ni?o events occurring between 2001~2013. During the four El Ni?o events since 2001, the NEC transport significantly increased during the mature phases of the 2002~2003、2006~2007 and 2009~2010 El Ni?o events, but did not significantly increase in the mature phase of the 2004~2005 El Ni?o event. We performed further analyses and found that negative sea surface height and cyclonic gyre anomalies took place south of 10° N in the tropical North Pacific Ocean during the mature phases of the 2002~2003、2006~2007 and 2009~2010 El Ni?o events. These were mainly induced by the westerly wind anomaly and positive Ekman pumping in the tropical gyre area through westward-propagating Rossby waves. However, during the mature phase of the 2004~2005 El Ni?o event, the distribution of sea surface temperature anomalies was obviously different. The westerly wind anomaly and positive Ekman pumping moved significantly northward, which led to negative sea surface height and cyclonic gyre anomalies taking place north of 10° N in the northwest Pacific Ocean. As such, the NEC transport did not significantly increase during the mature phase of the 2004~2005 El Ni?o.

    P731.27

    A

    1000-3096(2016)02-0118-10

    10.11759/hykx20150130002

    2015-01-30;

    2015-03-27

    “全球變化與海氣相互作用”任務專項(GASI-03-01-01-05);國家基礎研究項目(2012CB417401, 2013CB956202)

    吳國麗(1981-), 女, 山東菏澤人, 講師, 博士研究生, 主要從事物理海洋學研究, 電話: 18254230560, E-mail: wuguoli81@163.com

    猜你喜歡
    厄爾尼諾氣旋赤道
    “雪龍2”號過赤道,救下一船人!
    軍事文摘(2024年2期)2024-01-10 01:58:44
    2003年3月北大西洋上兩個爆發(fā)性氣旋的“吞并”過程及發(fā)展機制分析?
    氨還原劑噴槍中保護氣旋流噴射作用的研究
    能源工程(2021年1期)2021-04-13 02:05:56
    距赤道越遠越喜愛黃色
    厄爾尼諾現(xiàn)象橫行中國
    最熱的地方不在赤道
    厄爾尼諾現(xiàn)象及其影響
    北太平洋上一個爆發(fā)性氣旋族的結構分析?
    2014年3月大西洋上一個爆發(fā)性氣旋的研究
    三步圖解厄爾尼諾與拉尼娜
    地理教學(2015年17期)2016-01-06 12:00:28
    亚洲九九香蕉| 一边摸一边做爽爽视频免费| 精品视频人人做人人爽| 午夜福利在线观看吧| 欧美日韩成人在线一区二区| 精品欧美一区二区三区在线| 91在线观看av| 婷婷精品国产亚洲av在线 | 国产单亲对白刺激| 不卡av一区二区三区| 亚洲专区中文字幕在线| 女性生殖器流出的白浆| 制服诱惑二区| 窝窝影院91人妻| 曰老女人黄片| 午夜福利一区二区在线看| а√天堂www在线а√下载 | 一区福利在线观看| a级毛片在线看网站| 自拍欧美九色日韩亚洲蝌蚪91| 精品少妇久久久久久888优播| 国产精品国产高清国产av | 老司机福利观看| 一本综合久久免费| 免费少妇av软件| 日韩视频一区二区在线观看| 久99久视频精品免费| 欧美久久黑人一区二区| 少妇粗大呻吟视频| 露出奶头的视频| 日本一区二区免费在线视频| 王馨瑶露胸无遮挡在线观看| 三上悠亚av全集在线观看| 十八禁网站免费在线| 久久性视频一级片| 免费观看精品视频网站| 亚洲专区中文字幕在线| 国产有黄有色有爽视频| 可以免费在线观看a视频的电影网站| 久久 成人 亚洲| 婷婷成人精品国产| 美女高潮到喷水免费观看| 午夜精品在线福利| 人妻丰满熟妇av一区二区三区 | 老司机影院毛片| 国产国语露脸激情在线看| 免费观看精品视频网站| 人人妻,人人澡人人爽秒播| 久久香蕉国产精品| 人人妻人人澡人人看| 一级毛片女人18水好多| av线在线观看网站| 精品免费久久久久久久清纯 | 日韩欧美一区视频在线观看| a在线观看视频网站| 亚洲精品国产色婷婷电影| 午夜激情av网站| 国产不卡一卡二| 最近最新中文字幕大全免费视频| 免费在线观看亚洲国产| 看黄色毛片网站| 欧美+亚洲+日韩+国产| 女人久久www免费人成看片| 在线观看日韩欧美| 男男h啪啪无遮挡| 18禁观看日本| 90打野战视频偷拍视频| 狂野欧美激情性xxxx| 高清欧美精品videossex| tube8黄色片| 天堂中文最新版在线下载| 亚洲欧美激情综合另类| 亚洲午夜精品一区,二区,三区| 日韩成人在线观看一区二区三区| 香蕉丝袜av| 18禁裸乳无遮挡动漫免费视频| 在线国产一区二区在线| 国产成人系列免费观看| 中亚洲国语对白在线视频| 高清在线国产一区| 999久久久精品免费观看国产| 精品国产一区二区三区四区第35| 中文字幕另类日韩欧美亚洲嫩草| 精品欧美一区二区三区在线| 天堂俺去俺来也www色官网| 一区二区日韩欧美中文字幕| 午夜福利在线观看吧| 黄色丝袜av网址大全| 精品少妇一区二区三区视频日本电影| 国产av一区二区精品久久| 亚洲av欧美aⅴ国产| 高清在线国产一区| 18禁观看日本| 国产成人精品久久二区二区91| 少妇的丰满在线观看| 久久精品亚洲熟妇少妇任你| 1024香蕉在线观看| 两性夫妻黄色片| 看免费av毛片| 亚洲精品在线观看二区| 精品国产乱码久久久久久男人| 美女 人体艺术 gogo| 在线国产一区二区在线| 日韩 欧美 亚洲 中文字幕| 日韩人妻精品一区2区三区| 亚洲精品国产色婷婷电影| 建设人人有责人人尽责人人享有的| 两个人免费观看高清视频| 又大又爽又粗| 亚洲av成人一区二区三| 精品一品国产午夜福利视频| 国产精华一区二区三区| 日日爽夜夜爽网站| 久久国产精品男人的天堂亚洲| 性色av乱码一区二区三区2| 人妻久久中文字幕网| 欧美日韩瑟瑟在线播放| 午夜福利,免费看| 中文欧美无线码| 国产日韩一区二区三区精品不卡| 激情视频va一区二区三区| 一二三四社区在线视频社区8| 大香蕉久久成人网| 怎么达到女性高潮| 新久久久久国产一级毛片| 久久人妻熟女aⅴ| 美女扒开内裤让男人捅视频| 一边摸一边做爽爽视频免费| 久久久久精品人妻al黑| 日本撒尿小便嘘嘘汇集6| 欧美国产精品va在线观看不卡| 女人高潮潮喷娇喘18禁视频| а√天堂www在线а√下载 | 久久久国产精品麻豆| 啦啦啦免费观看视频1| 亚洲九九香蕉| 国产aⅴ精品一区二区三区波| 黄网站色视频无遮挡免费观看| 国产又色又爽无遮挡免费看| 国产精品久久久久久人妻精品电影| 18禁裸乳无遮挡动漫免费视频| 丝袜美足系列| 亚洲精品久久成人aⅴ小说| 视频在线观看一区二区三区| 精品久久久久久久久久免费视频 | 自线自在国产av| 亚洲av电影在线进入| 色婷婷av一区二区三区视频| 黄色成人免费大全| 久久久国产精品麻豆| 精品欧美一区二区三区在线| 建设人人有责人人尽责人人享有的| 一进一出好大好爽视频| 欧美日韩福利视频一区二区| 国产成人免费无遮挡视频| 久久精品国产a三级三级三级| 国内毛片毛片毛片毛片毛片| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 最新的欧美精品一区二区| 一级毛片女人18水好多| avwww免费| 久久性视频一级片| 热99国产精品久久久久久7| 国产蜜桃级精品一区二区三区 | 午夜视频精品福利| 老鸭窝网址在线观看| 免费在线观看日本一区| 搡老岳熟女国产| 18在线观看网站| 精品高清国产在线一区| 欧美午夜高清在线| 女警被强在线播放| 欧美激情久久久久久爽电影 | tube8黄色片| 国产黄色免费在线视频| 欧美乱码精品一区二区三区| 99精品欧美一区二区三区四区| 青草久久国产| 免费在线观看视频国产中文字幕亚洲| 亚洲视频免费观看视频| 夜夜夜夜夜久久久久| 在线观看66精品国产| 啦啦啦 在线观看视频| 精品欧美一区二区三区在线| 天天躁夜夜躁狠狠躁躁| 国产成人影院久久av| 看黄色毛片网站| 午夜福利乱码中文字幕| 美女 人体艺术 gogo| 老司机在亚洲福利影院| 午夜福利在线免费观看网站| 中文字幕最新亚洲高清| 日韩成人在线观看一区二区三区| 久久九九热精品免费| 咕卡用的链子| 亚洲成人国产一区在线观看| 亚洲人成伊人成综合网2020| 国产亚洲一区二区精品| 黄频高清免费视频| 久久人妻av系列| 日韩欧美国产一区二区入口| а√天堂www在线а√下载 | 欧美激情高清一区二区三区| 久久青草综合色| 波多野结衣av一区二区av| 91国产中文字幕| 欧美乱码精品一区二区三区| 捣出白浆h1v1| 高清在线国产一区| 日本欧美视频一区| 精品国产一区二区久久| 丰满饥渴人妻一区二区三| 极品教师在线免费播放| 一级毛片女人18水好多| av视频免费观看在线观看| 国产欧美日韩综合在线一区二区| 9热在线视频观看99| 大香蕉久久网| 91大片在线观看| 国产野战对白在线观看| 在线观看日韩欧美| 久久久精品区二区三区| 亚洲黑人精品在线| 国产在视频线精品| 岛国毛片在线播放| 最新美女视频免费是黄的| 精品亚洲成国产av| 国产国语露脸激情在线看| 国产精品免费视频内射| 男女免费视频国产| 亚洲色图综合在线观看| 亚洲精品一二三| 国产成人免费观看mmmm| 男人舔女人的私密视频| 久久午夜亚洲精品久久| 在线观看66精品国产| 国产av精品麻豆| 国产日韩一区二区三区精品不卡| 国产xxxxx性猛交| 69精品国产乱码久久久| 在线观看66精品国产| 俄罗斯特黄特色一大片| 啦啦啦在线免费观看视频4| 久久久精品免费免费高清| 日韩一卡2卡3卡4卡2021年| 中文亚洲av片在线观看爽 | 变态另类成人亚洲欧美熟女 | 18禁国产床啪视频网站| 亚洲第一欧美日韩一区二区三区| 少妇粗大呻吟视频| 欧美精品一区二区免费开放| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲 | 欧美激情极品国产一区二区三区| 在线播放国产精品三级| 国产淫语在线视频| 国产精品亚洲av一区麻豆| 不卡一级毛片| 国产精品 国内视频| 天天躁日日躁夜夜躁夜夜| 人妻一区二区av| 久久久久视频综合| 国产野战对白在线观看| 久久这里只有精品19| 丰满饥渴人妻一区二区三| 亚洲av熟女| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久| 动漫黄色视频在线观看| 天天操日日干夜夜撸| 看黄色毛片网站| 99精品欧美一区二区三区四区| 亚洲一区高清亚洲精品| 久久久久久人人人人人| 国产色视频综合| 国产免费男女视频| 露出奶头的视频| 黄色怎么调成土黄色| 久热这里只有精品99| 在线av久久热| 999久久久国产精品视频| 十八禁人妻一区二区| 欧美激情久久久久久爽电影 | 成年动漫av网址| 精品无人区乱码1区二区| 丰满饥渴人妻一区二区三| 国产精品九九99| 亚洲国产中文字幕在线视频| 搡老岳熟女国产| 成人黄色视频免费在线看| 在线十欧美十亚洲十日本专区| 99国产精品免费福利视频| 一级片'在线观看视频| 久久精品国产a三级三级三级| 一边摸一边抽搐一进一小说 | 一级,二级,三级黄色视频| 天天躁日日躁夜夜躁夜夜| 亚洲精品粉嫩美女一区| 免费在线观看日本一区| 天堂√8在线中文| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品| svipshipincom国产片| 国产野战对白在线观看| x7x7x7水蜜桃| 在线观看www视频免费| 最近最新中文字幕大全免费视频| netflix在线观看网站| 啦啦啦 在线观看视频| 欧美另类亚洲清纯唯美| 亚洲av成人不卡在线观看播放网| 叶爱在线成人免费视频播放| 王馨瑶露胸无遮挡在线观看| 90打野战视频偷拍视频| 在线国产一区二区在线| 国产男女超爽视频在线观看| 看免费av毛片| 亚洲精品国产精品久久久不卡| 美女扒开内裤让男人捅视频| 一区福利在线观看| 多毛熟女@视频| 亚洲国产毛片av蜜桃av| 久久精品国产99精品国产亚洲性色 | 又大又爽又粗| 又紧又爽又黄一区二区| 这个男人来自地球电影免费观看| 精品福利观看| 亚洲一码二码三码区别大吗| 色婷婷久久久亚洲欧美| 亚洲av熟女| av不卡在线播放| 99精国产麻豆久久婷婷| 午夜91福利影院| 大香蕉久久网| 男男h啪啪无遮挡| 国产欧美亚洲国产| xxxhd国产人妻xxx| 大香蕉久久网| 午夜日韩欧美国产| 视频区欧美日本亚洲| 国精品久久久久久国模美| 国产欧美日韩一区二区三| 人妻久久中文字幕网| 极品少妇高潮喷水抽搐| 日本wwww免费看| 欧美+亚洲+日韩+国产| 在线观看免费午夜福利视频| 国产乱人伦免费视频| 深夜精品福利| 国产精品美女特级片免费视频播放器 | 精品少妇久久久久久888优播| 免费av中文字幕在线| 18禁观看日本| 最新的欧美精品一区二区| 欧美最黄视频在线播放免费 | 三级毛片av免费| 91九色精品人成在线观看| 国产深夜福利视频在线观看| 日韩熟女老妇一区二区性免费视频| 无人区码免费观看不卡| 欧美 亚洲 国产 日韩一| 少妇被粗大的猛进出69影院| 久久九九热精品免费| 欧美乱码精品一区二区三区| 80岁老熟妇乱子伦牲交| 两性夫妻黄色片| 好看av亚洲va欧美ⅴa在| 国产高清激情床上av| 男人操女人黄网站| 在线国产一区二区在线| 一区福利在线观看| 在线观看日韩欧美| 午夜激情av网站| 国产精品自产拍在线观看55亚洲 | 成年动漫av网址| 亚洲人成伊人成综合网2020| 久久香蕉精品热| 欧美在线一区亚洲| 91在线观看av| 欧美日韩中文字幕国产精品一区二区三区 | 在线十欧美十亚洲十日本专区| 天天影视国产精品| 中文字幕人妻熟女乱码| 手机成人av网站| 久久精品国产亚洲av香蕉五月 | 一区二区三区激情视频| 亚洲九九香蕉| 国产精品一区二区在线不卡| 999精品在线视频| 男女之事视频高清在线观看| 国产精品一区二区在线不卡| 国产亚洲精品第一综合不卡| 黄片大片在线免费观看| 亚洲午夜理论影院| 国产精品免费一区二区三区在线 | 久久久久久久久免费视频了| 亚洲欧美色中文字幕在线| 18禁裸乳无遮挡免费网站照片 | 在线国产一区二区在线| 少妇被粗大的猛进出69影院| 国产成人免费观看mmmm| 操美女的视频在线观看| 一本大道久久a久久精品| 看片在线看免费视频| 午夜影院日韩av| 久久精品亚洲av国产电影网| 97人妻天天添夜夜摸| 女性生殖器流出的白浆| www.精华液| 成人国产一区最新在线观看| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲一级av第二区| 中文欧美无线码| 老熟妇仑乱视频hdxx| 人人妻人人澡人人看| 国产日韩欧美亚洲二区| 亚洲综合色网址| 国产av一区二区精品久久| 精品第一国产精品| 国产午夜精品久久久久久| 亚洲av成人一区二区三| 亚洲一码二码三码区别大吗| 黄网站色视频无遮挡免费观看| 免费女性裸体啪啪无遮挡网站| 成人永久免费在线观看视频| 精品国产乱码久久久久久男人| 国产在视频线精品| 日韩人妻精品一区2区三区| 亚洲熟女毛片儿| 69精品国产乱码久久久| 国产日韩欧美亚洲二区| 国产成人啪精品午夜网站| 欧美老熟妇乱子伦牲交| 亚洲 欧美一区二区三区| 日韩三级视频一区二区三区| 久久国产乱子伦精品免费另类| 国产熟女午夜一区二区三区| 国产深夜福利视频在线观看| 久久草成人影院| 成年动漫av网址| 久久人妻av系列| 十八禁网站免费在线| 丝袜美足系列| 亚洲aⅴ乱码一区二区在线播放 | 欧美性长视频在线观看| 久久天躁狠狠躁夜夜2o2o| 妹子高潮喷水视频| 国产精品偷伦视频观看了| 国产视频一区二区在线看| 成年人免费黄色播放视频| 最新美女视频免费是黄的| 精品国产亚洲在线| 新久久久久国产一级毛片| 久9热在线精品视频| 女同久久另类99精品国产91| 一级片免费观看大全| 手机成人av网站| 久久久久久久精品吃奶| 欧美乱色亚洲激情| 国产精品av久久久久免费| 热re99久久国产66热| 99国产精品免费福利视频| 久久久久精品国产欧美久久久| 12—13女人毛片做爰片一| 国产精品秋霞免费鲁丝片| 波多野结衣av一区二区av| 成人18禁在线播放| 国产一区二区激情短视频| 亚洲一区二区三区不卡视频| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 久久久精品国产亚洲av高清涩受| 亚洲成人免费av在线播放| √禁漫天堂资源中文www| 亚洲久久久国产精品| 看黄色毛片网站| 久久久久久久精品吃奶| 国产一区二区三区视频了| 日本vs欧美在线观看视频| avwww免费| 亚洲欧美日韩高清在线视频| 中文字幕人妻丝袜制服| 亚洲五月色婷婷综合| 欧美日韩一级在线毛片| 国产真人三级小视频在线观看| 在线播放国产精品三级| 久久久久久久午夜电影 | 亚洲熟女毛片儿| 无人区码免费观看不卡| 午夜福利乱码中文字幕| av国产精品久久久久影院| 国产有黄有色有爽视频| 久久香蕉精品热| 久久这里只有精品19| 久久亚洲精品不卡| 国产99久久九九免费精品| 老司机影院毛片| 欧美精品亚洲一区二区| 一进一出抽搐动态| 午夜福利欧美成人| 脱女人内裤的视频| 香蕉丝袜av| 国产亚洲精品久久久久久毛片 | 亚洲精品国产区一区二| 首页视频小说图片口味搜索| 国产一卡二卡三卡精品| 韩国av一区二区三区四区| 午夜福利在线观看吧| 大型黄色视频在线免费观看| 法律面前人人平等表现在哪些方面| 美女高潮喷水抽搐中文字幕| 亚洲专区字幕在线| 十八禁人妻一区二区| 欧美+亚洲+日韩+国产| www.精华液| 亚洲中文av在线| 曰老女人黄片| 18禁黄网站禁片午夜丰满| 欧美在线黄色| 黄色怎么调成土黄色| 久久人妻熟女aⅴ| 妹子高潮喷水视频| 亚洲第一欧美日韩一区二区三区| 国产精品国产av在线观看| 在线观看免费高清a一片| 啪啪无遮挡十八禁网站| 身体一侧抽搐| 中文字幕色久视频| 国产精品二区激情视频| 国产免费男女视频| 国产欧美日韩一区二区精品| 巨乳人妻的诱惑在线观看| 亚洲在线自拍视频| 性色av乱码一区二区三区2| 999久久久精品免费观看国产| 国产成人av教育| 一区福利在线观看| 一级毛片高清免费大全| 亚洲精品自拍成人| 新久久久久国产一级毛片| 国产精品一区二区在线不卡| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区精品视频观看| 后天国语完整版免费观看| 天天影视国产精品| 日日夜夜操网爽| 啪啪无遮挡十八禁网站| 日韩制服丝袜自拍偷拍| 丰满人妻熟妇乱又伦精品不卡| 国产午夜精品久久久久久| 老司机在亚洲福利影院| 国产精品美女特级片免费视频播放器 | 色婷婷av一区二区三区视频| 亚洲精品一二三| 在线观看舔阴道视频| 大香蕉久久成人网| 一级a爱片免费观看的视频| 夜夜躁狠狠躁天天躁| 天堂中文最新版在线下载| 悠悠久久av| 大片电影免费在线观看免费| 成人免费观看视频高清| 亚洲精品成人av观看孕妇| 国产97色在线日韩免费| 一二三四社区在线视频社区8| 久久精品亚洲精品国产色婷小说| 十八禁人妻一区二区| 亚洲国产毛片av蜜桃av| 精品国内亚洲2022精品成人 | 久久青草综合色| 亚洲国产精品一区二区三区在线| 亚洲精品国产区一区二| 精品欧美一区二区三区在线| 美女高潮到喷水免费观看| 国产乱人伦免费视频| aaaaa片日本免费| 日韩欧美三级三区| 超碰成人久久| 久久香蕉激情| 久久亚洲精品不卡| 搡老熟女国产l中国老女人| 成年人黄色毛片网站| 亚洲国产看品久久| 色综合婷婷激情| 亚洲五月婷婷丁香| 亚洲国产看品久久| 亚洲av片天天在线观看| 麻豆av在线久日| tube8黄色片| 国产激情欧美一区二区| av超薄肉色丝袜交足视频| 日本五十路高清| 精品少妇久久久久久888优播| 久久精品国产亚洲av高清一级| 黄色 视频免费看| 黑人猛操日本美女一级片| 成在线人永久免费视频| 国产精品久久久久久精品古装| 国产精品亚洲av一区麻豆| 岛国毛片在线播放| 国产视频一区二区在线看| 可以免费在线观看a视频的电影网站| 91精品三级在线观看| 欧美精品亚洲一区二区| 女人精品久久久久毛片| 在线天堂中文资源库| 亚洲专区国产一区二区| 少妇的丰满在线观看| 日韩熟女老妇一区二区性免费视频| 变态另类成人亚洲欧美熟女 | 波多野结衣一区麻豆| 一区二区日韩欧美中文字幕| 久久精品亚洲精品国产色婷小说|