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

    POM模式在日本南部黑潮路徑變異研究中的應(yīng)用

    2016-07-06 10:30:25鄒廣安中國(guó)科學(xué)院海洋研究所山東青島266071中國(guó)科學(xué)院大學(xué)北京100049河南大學(xué)河南開封475004
    海洋科學(xué) 2016年2期
    關(guān)鍵詞:斜壓黑潮渦度

    鄒廣安(1. 中國(guó)科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 河南大學(xué), 河南 開封475004)

    ?

    POM模式在日本南部黑潮路徑變異研究中的應(yīng)用

    鄒廣安1, 2, 3
    (1. 中國(guó)科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國(guó)科學(xué)院大學(xué), 北京 100049; 3. 河南大學(xué), 河南 開封475004)

    日本南部黑潮路徑變異對(duì)北太平洋地區(qū)的氣候和環(huán)境具有顯著的影響, 對(duì)黑潮路徑變異的研究具有重要的意義。本文利用POM (Princeton Ocean Model)數(shù)值模式模擬了日本南部黑潮的路徑變異情況, 分析了黑潮大彎曲路徑形成的可能機(jī)制。研究結(jié)果表明, 當(dāng)黑潮處于非大彎曲路徑時(shí), 相對(duì)位勢(shì)渦度的平均值呈現(xiàn)遞減趨勢(shì), 說明日本南部低位勢(shì)渦度水在不斷積累, 這樣會(huì)使得四國(guó)再循環(huán)流的強(qiáng)度增強(qiáng), 迫使黑潮保持平直路徑, 同時(shí), 近岸黑潮垂直流速剪切增大, 斜壓不穩(wěn)定性的作用也逐漸增大;當(dāng)黑潮從非大彎曲路徑向大彎曲路徑過渡時(shí), 再循環(huán)流強(qiáng)度的減弱會(huì)導(dǎo)致黑潮的流速剪切減小。根據(jù)海表高度異常場(chǎng)以及海洋上層流場(chǎng)信息發(fā)現(xiàn), 近岸黑潮附近的氣旋渦會(huì)隨著再循環(huán)流區(qū)域反氣旋渦的東側(cè)向南運(yùn)動(dòng), 最終導(dǎo)致黑潮大彎曲的發(fā)生。分析渦流的能量, 結(jié)果顯示, 黑潮大彎曲路徑的形成與斜壓不穩(wěn)定性密切相關(guān)。

    黑潮路徑變異; POM (Princeton Ocean Model); 斜壓不穩(wěn)定性

    [Foundation: National Natural Science Foundation of China, No. 41230420]

    黑潮(Kuroshio Current)是北太平洋副熱帶環(huán)流系統(tǒng)中非常著名的一支西邊界流, 又稱日本暖流,是世界上第二大暖流, 僅次于墨西哥灣暖流。黑潮因其水色深藍(lán), 遠(yuǎn)看似黑色而得名。一般認(rèn)為, 黑潮的起源主要是由北赤道暖流在菲律賓群島東部向北分叉形成的, 向南的一支稱為棉蘭老流。黑潮主軸沿著巴士海峽的東側(cè)北上, 經(jīng)臺(tái)灣東部進(jìn)入中國(guó)東海,從吐噶喇海峽北部流出東海, 進(jìn)入日本南部的北太平洋海域。黑潮具有高溫、高鹽、流速快、流量大等特點(diǎn), 其從低緯向中高緯輸送了大量的熱量, 對(duì)流經(jīng)區(qū)域乃至全球的氣候起著重要的作用。研究發(fā)現(xiàn), 黑潮流經(jīng)日本南部海域時(shí)會(huì)出現(xiàn)兩種典型的路徑形態(tài)[1]: 大彎曲路徑(LM)和非大彎曲路徑(NLM)。這兩種路徑形態(tài)和位置的變化是否發(fā)生大彎曲, 對(duì)北太平洋地區(qū)的氣候、漁業(yè)、航海安全以及海洋資源等都具有顯著的影響[2-3]。此外, 日本南部的黑潮路徑變異對(duì)我國(guó)東部沿岸以及長(zhǎng)江流域的降水也有較大的影響[4-5]。因此, 弄清黑潮路徑變異的原因和機(jī)制, 對(duì)黑潮路徑變異進(jìn)行預(yù)測(cè)具有重要的意義。

    前人對(duì)于黑潮大彎曲路徑形成的原因給出了多種可能的機(jī)制解釋。Kawabe[6]利用簡(jiǎn)單的動(dòng)力學(xué)模式研究發(fā)現(xiàn), 日本九州島東南部的黑潮小彎曲對(duì)于下游黑潮路徑有明顯的影響。Akitomo等[7]采用兩層海洋模式進(jìn)行研究, 認(rèn)為日本九州島東南部渦的增強(qiáng)有利于黑潮從平直路徑向大彎曲路徑轉(zhuǎn)換。Qiu等[8]利用兩層原始方程模式進(jìn)行研究, 認(rèn)為日本南部黑潮路徑變異取決于黑潮內(nèi)部的自持震蕩系統(tǒng)。Usui 等[9]利用高分辨率的大洋環(huán)流(OGCM)模式研究了2004年發(fā)生的典型黑潮大彎曲事件, 認(rèn)為日本九州島南部的黑潮觸發(fā)小彎曲是黑潮大彎曲形成的重要原因。Miyazawa等[10]研究認(rèn)為2004年黑潮大彎曲是由于臺(tái)灣以東強(qiáng)的反氣旋渦引起九州島南部形成觸發(fā)小彎曲, 進(jìn)而形成黑潮大彎曲。Tsujino等[11]認(rèn)為日本四國(guó)再循環(huán)流區(qū)域的渦的相互作用是黑潮大彎曲路徑形成的主要原因。

    目前, 由于海洋觀測(cè)資料有限, 人們對(duì)于黑潮大彎曲的形成原因和機(jī)制沒有統(tǒng)一的認(rèn)識(shí), 高分辨率的大洋環(huán)流模式(OGCM)為研究黑潮路徑變異現(xiàn)象提供了很大的便利和可能性。目前, 采用POM(Princeton Ocean Model)模式對(duì)黑潮的研究工作相對(duì)較少。Kagimoto[12]利用POM模式模擬了黑潮的季節(jié)變化。Miyazawa等[13]利用POM模式對(duì)黑潮路徑變異進(jìn)行了集合預(yù)報(bào)研究。Waseda等[14]采用修訂的POM模式研究了渦的相互作用對(duì)黑潮的影響。但是, 這些工作并沒有利用POM模式將黑潮的自持系統(tǒng)和渦的相互作用結(jié)合起來, 也沒有從渦流能量分析角度來考察黑潮的路徑變異情況。本文采用POM模式模擬了日本南部黑潮的路徑變異過程, 進(jìn)一步分析和討論了黑潮大彎曲形成的可能機(jī)制和原因。

    1 模式設(shè)置和模擬結(jié)果

    POM模式是由美國(guó)普林斯頓大學(xué)Alan Blumberg和George L. Mellor于1977年建立起來的一個(gè)基于三維斜壓原始方程海洋數(shù)值模式[15], 該模式經(jīng)過多年的不斷發(fā)展和完善, 目前已經(jīng)成為國(guó)內(nèi)外廣泛使用的海洋數(shù)值模式之一。

    模式的模擬區(qū)域設(shè)置為(30°S~60°N, 110°E~ 80°W), 覆蓋了整個(gè)北太平洋區(qū)域。模式的水平分辨率為0.2°×0.2°, 垂直方向分16層, 外模式的時(shí)間步長(zhǎng)為20 s, 內(nèi)模式的時(shí)間步長(zhǎng)為90 s。模式采用的是ETOPO5的分辨率為1/12°×1/12°的地形數(shù)據(jù), 模式的溫鹽資料采用的是Levitus WOA97分辨率為1/4°×1/4°的月平均氣候態(tài)溫鹽場(chǎng), 模式的風(fēng)應(yīng)力強(qiáng)迫場(chǎng)為NCEP/NCAR再分析的月平均風(fēng)應(yīng)力資料(真實(shí)時(shí)間的月平均風(fēng)場(chǎng)資料)。

    完成了以上的模式基本設(shè)置后, 首先將模式積分20 a, 每10 d輸出一次結(jié)果, 我們將后10 a的模式輸出結(jié)果作為分析數(shù)據(jù), 用來分析和檢驗(yàn)POM模式對(duì)黑潮路徑變異的模擬能力。為了考察黑潮路徑隨時(shí)間的變化, Qiu等[8]采用136°~140°E黑潮流軸離岸的平均距離定義黑潮路徑指數(shù)。類似于該定義, 我們也定義了一個(gè)指數(shù)來表征黑潮路徑所處的狀態(tài)。該黑潮路徑指數(shù)定義為: 黑潮流軸在136°~141°E最南端點(diǎn)的緯度, 黑潮流軸用海表高度(SSH)30 cm等值線表示。由定義可知, 黑潮路徑指數(shù)值較小時(shí), 表明黑潮流軸與日本南岸的距離越遠(yuǎn), 黑潮處于大彎曲路徑, 反之, 指數(shù)值較大時(shí), 黑潮處于非大彎曲路徑。圖1給出了POM模式模擬的黑潮路徑指數(shù)的時(shí)間序列, 可以看出, 黑潮的路徑變異主要呈現(xiàn)年際變化特征, 這10 a間共發(fā)生了3次大彎曲事件, 每次大彎曲事件發(fā)生的強(qiáng)度不同, 并且黑潮在大彎曲路徑和非大彎曲路徑的維持時(shí)間也有所相同。與觀測(cè)資料對(duì)比, POM模式模擬的黑潮大彎曲路徑的維持時(shí)間相對(duì)較短, 這與其他學(xué)者用POM模式模擬的黑潮路徑情況類似[10, 13], 這也是目前許多OGCM模式都存在的模擬缺陷[9, 16-17]。

    圖1 黑潮路徑指數(shù)的時(shí)間序列Fig. 1 Time series of the Kuroshio path index

    為了顯示POM模式模擬的黑潮路徑的變異情況, 圖2給出了模式第17~18年的一次黑潮大彎曲路徑的發(fā)生、發(fā)展和消亡過程。由圖2可知, 模式在第17年的2月時(shí), 黑潮處于典型的非大彎曲路徑, 在第17年的3月有個(gè)小彎曲形成(圖2b), 小彎曲逐漸發(fā)展, 到第17年的4月時(shí)形成了典型的大彎曲路徑(圖2c), 大彎曲維持了一段時(shí)間, 到第17年的8月黑潮的流軸開始形成一個(gè)渦旋(圖2e), 隨著渦旋的脫落并且西傳, 黑潮在第17年的12月又回到了典型的非大彎曲路徑(圖2h)。這一過程與Qiu等[8]采用2.5層淺水模式模擬的黑潮路徑變異過程也比較類似。該模式模擬的黑潮大彎曲的發(fā)生過程與觀測(cè)到的黑潮大彎曲發(fā)生比較一致[9]。這說明POM模式也能夠捕捉到黑潮路徑變異的基本特征, 能夠模擬黑潮路徑變異的基本動(dòng)力學(xué)過程。盡管模擬的黑潮大彎曲的維持時(shí)間相對(duì)較短, 但是, 由于本文主要討論的是黑潮大彎曲路徑形成的機(jī)制, 因此, 該模式對(duì)黑潮大彎曲路徑維持時(shí)間的模擬不足不會(huì)給我們的研究造成較大影響。

    2 黑潮大彎曲路徑形成的機(jī)制分析

    本文對(duì)模擬的三個(gè)大彎曲事件的形成過程進(jìn)行了分析, 三次大彎曲的形成過程和原因基本一致。因此, 主要以第17~18年黑潮大彎曲事件(圖2)為例給出詳細(xì)的分析和描述過程。

    根據(jù)Qiu等[8]的定義, 將相對(duì)渦度定義如公式(1):

    圖2 模式第17~18年發(fā)生的黑潮大彎曲事件的海表高度場(chǎng)序列Fig. 2 Sequence of SSH field including the Kuroshio large meander event in model year 17-18

    其中, u和v分別表示緯向和經(jīng)向流速, 相對(duì)渦度?表示的是四國(guó)再循環(huán)流區(qū)域(25°~35°N, 132°~140°E)內(nèi)的平均值, 緯向和徑向流速是由海洋上層500 m的流場(chǎng)平均值來計(jì)算的。

    為了表征黑潮再循環(huán)流的強(qiáng)度S, 根據(jù)Qiu等的定義[18], 定義了公式(2):

    其中, t表示時(shí)間, A表示四國(guó)再循環(huán)流區(qū)域(25°~ 35°N, 132°~140°E)內(nèi)SSH>0.2 m的區(qū)域, h表示海表面高度。

    此外, 我們將近岸黑潮的垂直流速剪切κ定義為公式(3):

    通過上述定義公式(1)~公式(3)進(jìn)行計(jì)算, 圖3給出了黑潮大彎曲及大彎曲事件發(fā)生前半年的相對(duì)渦度平均、四國(guó)再循環(huán)流的強(qiáng)度和近岸黑潮垂直流速剪切的時(shí)間序列。由圖可見, 當(dāng)黑潮處于非大彎曲路徑時(shí), 相對(duì)渦度的平均值出現(xiàn)了明顯的遞減趨勢(shì)(圖3a), 也就是說, 低位勢(shì)渦度水在不斷積累, 這樣會(huì)引起黑潮再循環(huán)流的強(qiáng)度增強(qiáng)(圖3b), 迫使黑潮保持平直路徑沿著日本南岸向東流, 此外, 四國(guó)再循環(huán)流的增強(qiáng)會(huì)導(dǎo)致黑潮的近岸流速剪切增大(圖3c)。當(dāng)黑潮從非大彎曲路徑向大彎曲路徑過渡時(shí),日本南部再循環(huán)流的強(qiáng)度減小, 黑潮的流速剪切也隨著減小, 由于正壓或斜壓不穩(wěn)定性的作用使得能量向渦動(dòng)能轉(zhuǎn)化, 最終形成大彎曲事件, 這與Qiu等用簡(jiǎn)單的兩層海洋模式得出的結(jié)論[8]一致。

    圖3 相對(duì)渦度平均(a)、黑潮再循環(huán)流強(qiáng)度(b)和近岸黑潮流速剪切(c)的時(shí)間序列Fig. 3 Time series of the average relative vorticity (a), the strength of the Shikoku recirculation gyre (b), and the vertical shear of the alongshore Kuroshio (c)

    為了考察渦漩對(duì)黑潮的影響, 圖4給出了黑潮大彎曲發(fā)生過程中海表高度異常(SSHA)場(chǎng)的信號(hào)情況(用背景場(chǎng)的值減去氣候態(tài)平均場(chǎng)的值), 可以看出, 在黑潮大彎曲形成過程中, 日本九州島東南部的冷異常信號(hào)向東傳播(圖4a, 圖4b), 并且冷異常在四國(guó)再循環(huán)流區(qū)域的強(qiáng)度逐漸增大(圖4c), 由于黑潮對(duì)日本南岸的流速剪切減小, 冷異常信號(hào)會(huì)向南運(yùn)動(dòng), 最終在黑潮流軸的彎曲處附近形成一個(gè)大的氣旋渦(圖4d), 這個(gè)氣旋渦對(duì)黑潮大彎曲的形成具有重要的影響。

    圖4 海表高度異常以及黑潮流軸(黑色實(shí)線)Fig. 4 SSH anomaly and the Kuroshio axes (black solid line)

    此外, 我們給出了黑潮大彎曲發(fā)生前一個(gè)月的海洋上層500 m平均水平流速場(chǎng)和海表高度場(chǎng), 如圖5所示, 可以看出, 在黑潮再循環(huán)流區(qū)域內(nèi)有一個(gè)反氣旋渦([135°E, 32°N]附近), 而在黑潮流軸北部([137.5°E, 33.5°N]附近)有一個(gè)氣旋渦存在, 當(dāng)黑潮的近岸流速剪切減小時(shí), 再循環(huán)流區(qū)域反氣旋渦的東側(cè)向南運(yùn)動(dòng), 同時(shí)會(huì)攜帶著氣旋渦向南運(yùn)動(dòng), 最終導(dǎo)致黑潮大彎曲的形成。

    采用渦流的能量分析來考察正壓或者斜壓不穩(wěn)定性對(duì)黑潮大彎曲的影響, 根據(jù)Tsujino等[11]提出的將背景流場(chǎng)分成平均場(chǎng)(前一個(gè)時(shí)間段的月平均值)和渦度場(chǎng)(當(dāng)前時(shí)間段的月平均值與前一時(shí)間段月平均值的差), 將平均動(dòng)能(MKE)向渦動(dòng)能(EKE)的正壓(barotropic, BT)轉(zhuǎn)化率RBT定義如下:

    圖5 海表高度場(chǎng)(單位: m)和海洋上層500m的平均流速場(chǎng)(單位: m/s)Fig. 5 SSH field (unit: m) and mean velocity field (unit: m/s) at the upper 500 m

    其中, ρ0表示Boussinesq近似下的海水密度,()表示平均場(chǎng)中的緯向(經(jīng)向)流速, u′( v′ )表示渦度場(chǎng)中的緯向(經(jīng)向)擾動(dòng)流速, 平均流速和擾動(dòng)流速都采用的是上層300 m的平均值, 該定義表明了正壓不穩(wěn)定性的作用。

    將平均位能(MPE)向渦位能(EPE)的斜壓(baroclinic, BC)轉(zhuǎn)化率RBC定義為:

    其中, g表示重力加速度, u′=(u′ , v′ ,0)表示擾動(dòng)流速場(chǎng),和δρ′分別平均場(chǎng)和渦度場(chǎng)密度,表示背景密度場(chǎng)在密度z處的垂向變化率, ▽表示水平梯度算子。

    另外, 考慮到渦浮力通量項(xiàng)δρ′u′的旋轉(zhuǎn)分量在平均場(chǎng)和渦度場(chǎng)的轉(zhuǎn)換中作用不大[11], 因此, 采用斜壓轉(zhuǎn)化率方程中渦浮力通量項(xiàng)的散度(div)分量來計(jì)算動(dòng)力斜壓轉(zhuǎn)換率RdynBC, 考慮動(dòng)力過程對(duì)能量轉(zhuǎn)化的影響:

    上述的渦能量轉(zhuǎn)化率公式的詳細(xì)推導(dǎo)以及具體計(jì)算過程可參閱文獻(xiàn)[11]。

    圖6 300 m深度黑潮大彎曲發(fā)生過程中的能量轉(zhuǎn)化率(單位:  ×10-3 kg /(m· s3))以及海表面高度(單位: m)Fig. 6 Energy conversion rates (unit:  ×10-3 kg /(m· s3)) at 300 m

    通過三個(gè)能量轉(zhuǎn)化公式的計(jì)算, 圖6給出了黑潮大彎曲發(fā)生過程中上層300 m處的正壓能量轉(zhuǎn)化率(第一行)、斜壓能量轉(zhuǎn)化率(第二行)和動(dòng)力斜壓能量轉(zhuǎn)化率(第三行)情況, 垂直深度之所以選取300 m,該深度的渦能量轉(zhuǎn)化率比較大, 標(biāo)準(zhǔn)差也較大。從圖中可以發(fā)現(xiàn), 渦度能量分布的主要區(qū)域位于黑潮流軸及其流軸附近, 這可能是由于黑潮流軸的流速較大的原因。此外, 斜壓和動(dòng)力斜壓能量轉(zhuǎn)化率的信息要比正壓能量轉(zhuǎn)化率的信息量大, 表明黑潮大彎曲發(fā)生過程中斜壓不穩(wěn)定性起著主要的作用。同時(shí), 我們計(jì)算了黑潮大彎曲事件發(fā)生及其發(fā)生前半年時(shí)間的能量轉(zhuǎn)化率在日本四國(guó)再循環(huán)流區(qū)域(25°~35°N, 132°~140°E)內(nèi)絕對(duì)值的平均值, 如圖7所示, 我們發(fā)現(xiàn)斜壓轉(zhuǎn)化率的值要大于正壓轉(zhuǎn)化率的值, 并且, 隨著近岸流速剪切的增大(圖3c), 斜壓不穩(wěn)定的值也在增大。也就是說,斜壓不穩(wěn)定在黑潮大彎曲發(fā)生過程中的作用要比正壓不穩(wěn)定性的作用大, 是引發(fā)黑潮大彎曲事件的重要因素。

    圖7 3種能量轉(zhuǎn)化率在區(qū)域(25°~35°N, 132°~140°E)內(nèi)絕對(duì)值的平均值的時(shí)間序列(單位: ×10–8 kg /(m· s3))Fig. 7 Time series of the mean absolute value of three conversion rates over the region (25°–35°N, 132°–140°E) (unit: × 10–8 kg /(m· s3)).

    3 結(jié)論

    本文基于前人對(duì)日本南部黑潮大彎曲現(xiàn)象的研究結(jié)果, 利用POM模式成功模擬了黑潮路徑變異情況, 并且分析了黑潮大彎曲形成的可能原因和機(jī)制。結(jié)果表明, 黑潮的自身動(dòng)力過程以及南部再循環(huán)區(qū)域渦的作用能夠促使黑潮發(fā)生大彎曲事件。日本南部低位渦水的不斷積累, 引起黑潮再循環(huán)流的強(qiáng)度增強(qiáng), 使得黑潮保持平直路徑, 增強(qiáng)的再循環(huán)流會(huì)導(dǎo)致黑潮的近岸流速剪切增大, 從而, 斜壓不穩(wěn)定性的作用也在逐漸增大; 當(dāng)黑潮從非大彎曲路徑向大彎曲路徑開始過渡時(shí), 日本南部再循環(huán)流的強(qiáng)度減小, 黑潮的流速剪切也隨著減小, 黑潮流軸附近的上層氣旋渦會(huì)隨著再循環(huán)流區(qū)域反氣旋渦的東側(cè)向南運(yùn)動(dòng), 導(dǎo)致黑潮大彎曲的形成。利用渦度能量分析理論進(jìn)一步驗(yàn)證了黑潮大彎曲的發(fā)生與斜壓不穩(wěn)定性密切相關(guān)。鑒于該模式能夠很好地模擬黑潮大彎曲事件, 接下來將采用該模式來進(jìn)行數(shù)值預(yù)報(bào)研究, 考慮目標(biāo)觀測(cè)效果對(duì)黑潮大彎曲路徑預(yù)報(bào)的影響。

    致謝: 本文的研究工作得到了中國(guó)科學(xué)院海洋研究所穆穆院士的悉心指導(dǎo), 在此表示衷心的感謝!

    [1] Taft BA. Characteristics of the flow of the Kuroshio south of Japan[C] //Stommel H, Yoshida K.Kuroshio: Its Physical Aspects.Tokyo: University of Tokyo Press, 1972: 165-216.

    [2] 馮士筰, 李鳳岐, 李少菁. 海洋科學(xué)導(dǎo)論[M]. 北京:高等教育出版社, 1999, 170-171.

    Feng Shizuo, Li Fengqi, Li Shaojing. An Introduction to Marine Science[M]. Beijing: Higher Education Press, 1999, 170-171.

    [3] Xu H M, Tokinaga H, Xie S P. Atmospheric effects of the Kuroshio large meander during 2004-05[J]. J Climate, 2010, 23: 4704-4715.

    [4] 于振娟. 東海黑潮流軸的變化及日本以南黑潮大彎曲同青島降水量的關(guān)系[J]. 海洋科學(xué), 1998, 4: 6-11.

    Yu Zhenjuan. Relationship of variation of axial position of the Kuroshio in the East China sea and its meanders south of Japan with the precipitation in Qingdao[J]. Marine Sciences, 1998, 4: 6-11.

    [5] 侍茂崇. 物理海洋學(xué)[M]. 濟(jì)南: 山東教育出版社, 2004: 157.

    Shi Maochong. Physical oceanography[M]. Jinan: Shandong Education Press, 2004: 157.

    [6] Kawabe M. Model study of flow conditions causing the large meander of the Kuroshio in relation with the large meander[J]. J Phys Oceanogr, 1996, 26: 2449-2461.

    [7] Akitomo K, Kurogi M. Path transition of the Kuroshio due to mesoscale eddies: A two-layer, wind-driven experiment [J]. J Oceanogr, 2001, 57(6): 735-741.

    [8] Qiu B, Miao W F. Kuroshio path variations south of Japan: Bimodality as a self-sustained internal oscillation[J]. J Phys Oceanogr, 2000, 30: 2124-2137.

    [9] Usui N, Tsujino H, Nakano H, et al. Formation process of the Kuroshio large meander in 2004 [J].J Geophys Res, 2008, 113: C08047.

    [10] MiyazawaY, Kagimoto T, Guo X, et al. The Kuroshio large meander formation in 2004 analyzed by an eddy-resolving ocean forecast system[J]. J Geophys Res, 2008, 113: C10015.

    [11] Tsujino H, Usui N, Nakano H. Dynamics of Kuroshio path variations in a high-resolution general circulationmodel[J]. J Geophys Res, 2006, 111: C11001.

    [12] Kagimoto T, YamagataT. Seasonal Transport Variations of the Kuroshio: An OGCM Simulation[J]. J Phys Oceanogr, 1997, 27(3): 403-418.

    [13] Miyazawa Y, Yamane S, Guo X, et al. Ensemble forecast of the Kuroshio meandering[J]. J Geophys Res, 2005, 110: C10026.

    [14] Waseda T, Mitsudera H, Taguchi B, et al. On the eddy-Kuroshio interaction: Meander formation processes[J]. J Geophys Res, 2003, 108: C73220.

    [15] Blumberg A F, Mellor G L. Diagnostic and prognostic numerical circulation studies of the South Atlantic Bight[J]. J Geophys Res, 1983, 88: 4579-4592.

    [16] Usui N, Tsujino H, Fujii Y, et al. Generation of a trigger meander for the 2004 Kuroshio large meander[J]. J Geophys Res, 2008, 113: C01012.

    [17] Masumoto Y. Generation of small meanders of the Kuroshio south of Kyushu in a high-resolution ocean general circulation model[J]. J Oceanogr, 2004, 60: 313-320.

    [18] QiuB, Chen S M. Variability of the Kuroshio Extension jet, recirculation gyre and mesoscale eddies on decadal time scales[J]. J Phys Oceanogr, 2005, 35: 2090-2103.

    (本文編輯: 李曉燕)

    Application of the POM model to study the path variations of Kuroshio currents in southern Japan seas

    ZOU Guang-an1, 2, 3
    (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. Henan University, Kaifeng 475004, China)

    Mar., 31, 2015

    Kuroshio path variations; POM model; Baroclinic instability

    Path variations of the Kuroshio currents in the southern Japan seas have a significant impact on the climate and environment. In this study, we used the Princeton Ocean Model (POM) to simulate the Kuroshio path variations south of Japan, and analyzed the formation mechanism of the Kuroshio large meander events. The results show a decreasing trend of the average relative vorticity when the Kuroshio is taken as the non-large path, indicating an accumulation of the low potential vorticity (PV) south of Japan which intensified the strength of the Shikoku recirculation gyre, thereby forcing the Kuroshio into its straight path. Meanwhile, the vertical velocity shear of the alongshore Kuroshio was increasing, and the impact of the baroclinic instability was also gradually increasing. During the transition period from the Kuroshio non-large meander to the large meander, a sudden release of velocity shear corresponded well to the weakening of the Shikoku recirculation gyre. According to the sea surface height (SSH) anomaly fields and the mean velocity field in the upper ocean, we found that the cyclonic eddy around the offshore Kuroshio moved southward, carried by a southward flow in the eastern part of the anticyclonic eddy in the region of the Shikoku recirculation gyre. This led to the eventual development of the Kuroshio large meander. The results of our eddy energetic analysis indicate that the formation of the Kuroshio large meander is closely associated with baroclinic instability.

    P731.21

    A

    1000-3096(2016)02-0151-08

    10.11759/hykx20150127001

    2015-01-27;

    2015-03-31

    國(guó)家自然科學(xué)基金(41230420)

    鄒廣安(1985-), 男, 陜西安康人, 博士研究生, 主要從事西邊界流的數(shù)值模擬和預(yù)報(bào)研究, 電話: 15736863632, E-mail: zouguangan00@163.com

    猜你喜歡
    斜壓黑潮渦度
    譜松弛動(dòng)力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    南海相對(duì)渦度的時(shí)空變化特征分析*
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計(jì)數(shù)據(jù)的黑潮變異特征
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    渤海夏季第一斜壓羅斯貝變形半徑的計(jì)算與分析
    東海冷渦對(duì)黑潮熱輸運(yùn)年代際躍變的響應(yīng)
    免费观看在线日韩| 一区二区三区免费毛片| 人人妻人人澡欧美一区二区| 久久久久久久久久成人| 久久久久久久久久久丰满| 简卡轻食公司| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 亚洲不卡免费看| 亚洲高清免费不卡视频| 久久精品综合一区二区三区| 99视频精品全部免费 在线| 欧美最黄视频在线播放免费| 亚洲国产精品成人久久小说 | 97超视频在线观看视频| 亚洲欧美精品自产自拍| 亚洲无线观看免费| 国产激情偷乱视频一区二区| 麻豆精品久久久久久蜜桃| 日本av手机在线免费观看| 久久这里只有精品中国| 欧美精品一区二区大全| 国产免费男女视频| 高清在线视频一区二区三区 | 免费av毛片视频| ponron亚洲| 国产伦在线观看视频一区| 又爽又黄a免费视频| 极品教师在线视频| 亚洲欧美成人综合另类久久久 | 看免费成人av毛片| 观看免费一级毛片| 日韩亚洲欧美综合| 18禁在线播放成人免费| 日日撸夜夜添| 热99在线观看视频| 成人特级黄色片久久久久久久| 国产亚洲av片在线观看秒播厂 | 亚洲av电影不卡..在线观看| 在现免费观看毛片| 色综合站精品国产| 丰满人妻一区二区三区视频av| 亚洲精品亚洲一区二区| 成人特级黄色片久久久久久久| 亚洲第一区二区三区不卡| 日本-黄色视频高清免费观看| 青春草视频在线免费观看| 啦啦啦啦在线视频资源| 亚洲va在线va天堂va国产| 97在线视频观看| www.色视频.com| 亚洲中文字幕一区二区三区有码在线看| 特级一级黄色大片| 一个人看视频在线观看www免费| 麻豆av噜噜一区二区三区| 欧美成人a在线观看| 18禁在线播放成人免费| 少妇裸体淫交视频免费看高清| 中国国产av一级| 最近中文字幕高清免费大全6| a级毛片a级免费在线| 久久99热这里只有精品18| 国产精品综合久久久久久久免费| 大香蕉久久网| 亚洲丝袜综合中文字幕| 欧美一级a爱片免费观看看| 久久久久久久亚洲中文字幕| 精品国内亚洲2022精品成人| 亚洲第一电影网av| 最近2019中文字幕mv第一页| 91麻豆精品激情在线观看国产| 国产乱人视频| 国产乱人视频| 美女脱内裤让男人舔精品视频 | 国产精品乱码一区二三区的特点| 精品国内亚洲2022精品成人| 欧美极品一区二区三区四区| 给我免费播放毛片高清在线观看| 日日干狠狠操夜夜爽| АⅤ资源中文在线天堂| 婷婷六月久久综合丁香| 只有这里有精品99| 成年女人看的毛片在线观看| 天天一区二区日本电影三级| 国产伦精品一区二区三区四那| 色尼玛亚洲综合影院| 日本撒尿小便嘘嘘汇集6| 国产精华一区二区三区| 欧美激情在线99| 狂野欧美激情性xxxx在线观看| 国产精品,欧美在线| 精华霜和精华液先用哪个| 亚洲精品日韩在线中文字幕 | 99久久中文字幕三级久久日本| 又粗又爽又猛毛片免费看| 寂寞人妻少妇视频99o| 你懂的网址亚洲精品在线观看 | 日本黄色片子视频| 亚洲国产精品成人久久小说 | 国产探花在线观看一区二区| 午夜亚洲福利在线播放| 久久精品国产清高在天天线| 悠悠久久av| 国产高清不卡午夜福利| av福利片在线观看| 天堂中文最新版在线下载 | 桃色一区二区三区在线观看| 人妻少妇偷人精品九色| 久久九九热精品免费| 禁无遮挡网站| 国产麻豆成人av免费视频| 国产黄色视频一区二区在线观看 | 亚洲精品日韩av片在线观看| 久久精品人妻少妇| 国产精品一区二区性色av| 美女脱内裤让男人舔精品视频 | 91aial.com中文字幕在线观看| 91久久精品电影网| 一级黄色大片毛片| 黄色一级大片看看| 欧美激情国产日韩精品一区| 久久久久久久久中文| 成人一区二区视频在线观看| 亚洲在线自拍视频| av在线老鸭窝| 中文字幕人妻熟人妻熟丝袜美| 精品99又大又爽又粗少妇毛片| 欧美丝袜亚洲另类| 精品人妻熟女av久视频| 国产亚洲av嫩草精品影院| videossex国产| 少妇猛男粗大的猛烈进出视频 | 精品久久久久久久末码| 舔av片在线| 国产精品日韩av在线免费观看| 中文精品一卡2卡3卡4更新| 91麻豆精品激情在线观看国产| 亚洲国产欧美人成| 色播亚洲综合网| 免费观看人在逋| 日日干狠狠操夜夜爽| 中国美女看黄片| 18禁裸乳无遮挡免费网站照片| 男人和女人高潮做爰伦理| 久久九九热精品免费| 国产老妇伦熟女老妇高清| 少妇人妻精品综合一区二区 | 国产成人午夜福利电影在线观看| 欧美日本视频| 乱码一卡2卡4卡精品| 亚洲图色成人| 18+在线观看网站| 夫妻性生交免费视频一级片| 又粗又爽又猛毛片免费看| 日韩成人伦理影院| 国产精品不卡视频一区二区| 国产激情偷乱视频一区二区| 久久精品国产亚洲av涩爱 | 久久精品夜夜夜夜夜久久蜜豆| 五月伊人婷婷丁香| 日韩中字成人| avwww免费| 男女做爰动态图高潮gif福利片| 男女做爰动态图高潮gif福利片| 精品国产三级普通话版| 亚洲欧美成人综合另类久久久 | 18禁黄网站禁片免费观看直播| 深夜a级毛片| eeuss影院久久| 精品国产三级普通话版| 精品一区二区免费观看| 国产私拍福利视频在线观看| 听说在线观看完整版免费高清| 午夜亚洲福利在线播放| 婷婷色综合大香蕉| 久久久精品94久久精品| 国产成人a∨麻豆精品| 国产av在哪里看| 最好的美女福利视频网| 91aial.com中文字幕在线观看| 男女做爰动态图高潮gif福利片| АⅤ资源中文在线天堂| 日日干狠狠操夜夜爽| av女优亚洲男人天堂| 午夜免费激情av| 麻豆成人午夜福利视频| 欧美3d第一页| 精品无人区乱码1区二区| 在线播放国产精品三级| 老女人水多毛片| 国产伦理片在线播放av一区 | 欧美激情国产日韩精品一区| 色5月婷婷丁香| av在线亚洲专区| 国产精品国产高清国产av| 欧美精品一区二区大全| 九九久久精品国产亚洲av麻豆| 免费不卡的大黄色大毛片视频在线观看 | 国产国拍精品亚洲av在线观看| 国产男人的电影天堂91| 欧美日韩乱码在线| 一本久久精品| 女人被狂操c到高潮| 亚洲七黄色美女视频| 欧美日韩乱码在线| 亚洲欧美日韩卡通动漫| 欧美xxxx黑人xx丫x性爽| 国产精品嫩草影院av在线观看| 成人av在线播放网站| 麻豆精品久久久久久蜜桃| 女的被弄到高潮叫床怎么办| 精品久久久久久久久久久久久| 亚洲人成网站高清观看| 国产精品电影一区二区三区| 国产在线男女| 看免费成人av毛片| 欧美日韩一区二区视频在线观看视频在线 | 69av精品久久久久久| 一区二区三区免费毛片| 国产精品一区www在线观看| 夜夜夜夜夜久久久久| 精品99又大又爽又粗少妇毛片| 久久久久久久亚洲中文字幕| 日韩欧美在线乱码| 亚洲电影在线观看av| 国产不卡一卡二| 亚洲真实伦在线观看| 国产一区二区在线观看日韩| 长腿黑丝高跟| 舔av片在线| 亚洲成人久久爱视频| 国产精品麻豆人妻色哟哟久久 | 亚洲真实伦在线观看| 国产一区二区在线观看日韩| 免费不卡的大黄色大毛片视频在线观看 | 最近手机中文字幕大全| 精品一区二区三区视频在线| 中文字幕人妻熟人妻熟丝袜美| 日韩国内少妇激情av| 我的老师免费观看完整版| 日本与韩国留学比较| 亚洲乱码一区二区免费版| 久久久久久久久久黄片| 最近视频中文字幕2019在线8| 性色avwww在线观看| 欧美xxxx性猛交bbbb| 精品一区二区三区视频在线| 欧美成人免费av一区二区三区| 亚洲欧美日韩高清在线视频| 一级黄色大片毛片| 欧美一区二区精品小视频在线| 亚洲乱码一区二区免费版| 免费av毛片视频| 最近视频中文字幕2019在线8| videossex国产| 韩国av在线不卡| 老师上课跳d突然被开到最大视频| 日韩欧美在线乱码| 一级黄色大片毛片| 亚洲精品久久久久久婷婷小说 | 国产午夜精品论理片| 亚洲中文字幕日韩| 村上凉子中文字幕在线| 日韩欧美三级三区| 欧美日韩乱码在线| 亚洲欧美中文字幕日韩二区| 欧美极品一区二区三区四区| 久久久a久久爽久久v久久| 看免费成人av毛片| 国产精品福利在线免费观看| 亚洲欧美成人综合另类久久久 | 自拍偷自拍亚洲精品老妇| 99热网站在线观看| 一区二区三区四区激情视频 | 婷婷精品国产亚洲av| 寂寞人妻少妇视频99o| 看非洲黑人一级黄片| 少妇裸体淫交视频免费看高清| 日韩 亚洲 欧美在线| 国产精品电影一区二区三区| 12—13女人毛片做爰片一| 黄片wwwwww| 乱人视频在线观看| 国产黄色视频一区二区在线观看 | 欧美另类亚洲清纯唯美| 色尼玛亚洲综合影院| 看片在线看免费视频| 欧美最黄视频在线播放免费| 亚洲,欧美,日韩| 大型黄色视频在线免费观看| 国产男人的电影天堂91| 能在线免费看毛片的网站| 免费看日本二区| a级一级毛片免费在线观看| 亚洲成人av在线免费| 久久国产乱子免费精品| 久久久久久国产a免费观看| 中国美女看黄片| 国产一级毛片在线| 麻豆一二三区av精品| 亚洲国产高清在线一区二区三| 黄片wwwwww| av在线天堂中文字幕| 99在线视频只有这里精品首页| 少妇高潮的动态图| 国产精品1区2区在线观看.| 精华霜和精华液先用哪个| av在线播放精品| 99久国产av精品| 成人永久免费在线观看视频| 免费av不卡在线播放| 亚洲成人中文字幕在线播放| av在线天堂中文字幕| 人妻系列 视频| 日韩成人伦理影院| 国产av在哪里看| 校园人妻丝袜中文字幕| 欧美3d第一页| 国产精品不卡视频一区二区| 国产乱人视频| 欧美性猛交╳xxx乱大交人| 国产人妻一区二区三区在| 午夜精品一区二区三区免费看| 成人毛片60女人毛片免费| 草草在线视频免费看| 亚洲国产最新在线播放| 涩涩av久久男人的天堂| 永久网站在线| 国产国拍精品亚洲av在线观看| 久久久久精品性色| 成人综合一区亚洲| 亚洲第一区二区三区不卡| 一区二区av电影网| 一二三四中文在线观看免费高清| 国产69精品久久久久777片| 91精品一卡2卡3卡4卡| 久久免费观看电影| 亚洲精品乱久久久久久| 精品国产国语对白av| 熟妇人妻不卡中文字幕| 久久久国产精品麻豆| 色婷婷av一区二区三区视频| 校园人妻丝袜中文字幕| 麻豆精品久久久久久蜜桃| 欧美变态另类bdsm刘玥| 两个人的视频大全免费| 爱豆传媒免费全集在线观看| 午夜老司机福利剧场| 久久精品国产鲁丝片午夜精品| 性高湖久久久久久久久免费观看| 日本午夜av视频| 青青草视频在线视频观看| 青春草国产在线视频| 97超碰精品成人国产| 亚洲激情五月婷婷啪啪| 色哟哟·www| a 毛片基地| 日韩精品免费视频一区二区三区 | 国产69精品久久久久777片| 熟女av电影| 国产精品偷伦视频观看了| 搡老乐熟女国产| 少妇人妻精品综合一区二区| 国产不卡av网站在线观看| 日韩人妻高清精品专区| 亚洲精品国产色婷婷电影| 啦啦啦啦在线视频资源| 日韩熟女老妇一区二区性免费视频| 久久久久久伊人网av| 插阴视频在线观看视频| 精品99又大又爽又粗少妇毛片| 五月开心婷婷网| 草草在线视频免费看| 99九九在线精品视频| 一级毛片aaaaaa免费看小| 超色免费av| 综合色丁香网| 欧美日韩亚洲高清精品| 伊人久久精品亚洲午夜| 午夜av观看不卡| 91精品伊人久久大香线蕉| 少妇精品久久久久久久| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 久久精品国产鲁丝片午夜精品| 国产av码专区亚洲av| 亚洲综合色网址| 建设人人有责人人尽责人人享有的| 欧美激情极品国产一区二区三区 | 边亲边吃奶的免费视频| 欧美另类一区| 精品99又大又爽又粗少妇毛片| 国产免费福利视频在线观看| av天堂久久9| 纯流量卡能插随身wifi吗| 性高湖久久久久久久久免费观看| 国产极品粉嫩免费观看在线 | 国产一级毛片在线| 岛国毛片在线播放| 极品人妻少妇av视频| 春色校园在线视频观看| 黄色欧美视频在线观看| 欧美精品一区二区大全| av一本久久久久| 婷婷色av中文字幕| 性高湖久久久久久久久免费观看| 色婷婷av一区二区三区视频| 亚洲四区av| 午夜视频国产福利| 日本与韩国留学比较| av国产精品久久久久影院| 国产日韩欧美亚洲二区| 国产免费福利视频在线观看| 中文字幕最新亚洲高清| 亚洲av在线观看美女高潮| 中文字幕人妻丝袜制服| 久久久久久久精品精品| tube8黄色片| 欧美精品一区二区免费开放| 性色avwww在线观看| 中国美白少妇内射xxxbb| 成人国产麻豆网| 亚洲av免费高清在线观看| 久久久精品免费免费高清| 国产 一区精品| 中文字幕精品免费在线观看视频 | 男人添女人高潮全过程视频| 观看av在线不卡| 国产 精品1| 少妇高潮的动态图| 我的女老师完整版在线观看| 日韩亚洲欧美综合| 午夜福利,免费看| 2018国产大陆天天弄谢| 国产日韩欧美视频二区| 女性生殖器流出的白浆| 亚州av有码| 国产精品久久久久久精品古装| 免费高清在线观看日韩| av一本久久久久| 91精品国产九色| av网站免费在线观看视频| 一二三四中文在线观看免费高清| 大又大粗又爽又黄少妇毛片口| 午夜av观看不卡| 亚洲国产最新在线播放| 国产午夜精品一二区理论片| 精品国产乱码久久久久久小说| 大香蕉97超碰在线| 亚洲精品美女久久av网站| 日韩av免费高清视频| 在线亚洲精品国产二区图片欧美 | av在线播放精品| 国产日韩欧美在线精品| 免费高清在线观看视频在线观看| 99国产综合亚洲精品| 3wmmmm亚洲av在线观看| 中文字幕精品免费在线观看视频 | 亚洲精品一区蜜桃| 91久久精品电影网| 亚洲精品久久午夜乱码| 午夜福利视频在线观看免费| 2022亚洲国产成人精品| 好男人视频免费观看在线| 日韩一区二区视频免费看| 成人黄色视频免费在线看| 国产精品久久久久成人av| 日日摸夜夜添夜夜添av毛片| 极品人妻少妇av视频| 狂野欧美白嫩少妇大欣赏| 五月开心婷婷网| 一边亲一边摸免费视频| 晚上一个人看的免费电影| 欧美精品一区二区大全| 日日摸夜夜添夜夜爱| h视频一区二区三区| 国产一区二区在线观看av| 国产亚洲一区二区精品| 制服丝袜香蕉在线| 一二三四中文在线观看免费高清| 亚洲人成网站在线观看播放| 最黄视频免费看| 狂野欧美激情性bbbbbb| 国产精品99久久久久久久久| 国产欧美另类精品又又久久亚洲欧美| 成人二区视频| 久久热精品热| 夜夜骑夜夜射夜夜干| 欧美+日韩+精品| 国产国拍精品亚洲av在线观看| 欧美成人精品欧美一级黄| 亚洲av日韩在线播放| 黄色怎么调成土黄色| 久久久精品免费免费高清| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 午夜av观看不卡| 亚洲五月色婷婷综合| 永久网站在线| 一区二区三区四区激情视频| h视频一区二区三区| 亚洲国产精品专区欧美| 亚洲av不卡在线观看| 亚洲精品自拍成人| 一区二区三区乱码不卡18| 伊人久久精品亚洲午夜| 欧美国产精品一级二级三级| 丝袜喷水一区| 精品久久久久久电影网| 欧美 日韩 精品 国产| 日产精品乱码卡一卡2卡三| 国产片特级美女逼逼视频| 久久女婷五月综合色啪小说| 一级二级三级毛片免费看| 久久久久久久亚洲中文字幕| 哪个播放器可以免费观看大片| 欧美日韩av久久| 亚洲成人av在线免费| 色视频在线一区二区三区| 午夜免费观看性视频| 久久久久久久久大av| 飞空精品影院首页| 亚洲精品自拍成人| 免费观看的影片在线观看| 少妇高潮的动态图| 久久影院123| 日韩精品免费视频一区二区三区 | a级片在线免费高清观看视频| 亚洲av福利一区| 亚洲av在线观看美女高潮| 制服诱惑二区| 亚洲精品第二区| 国产女主播在线喷水免费视频网站| 王馨瑶露胸无遮挡在线观看| 一个人看视频在线观看www免费| √禁漫天堂资源中文www| 久久久a久久爽久久v久久| 久久久久精品久久久久真实原创| 国产视频内射| 国产精品 国内视频| 日本欧美国产在线视频| 黑丝袜美女国产一区| 卡戴珊不雅视频在线播放| 高清视频免费观看一区二区| 国产精品三级大全| 亚洲欧美色中文字幕在线| 少妇人妻久久综合中文| 亚洲美女搞黄在线观看| 日韩中字成人| 女的被弄到高潮叫床怎么办| 久久精品久久久久久噜噜老黄| 色婷婷久久久亚洲欧美| 成人毛片60女人毛片免费| 免费日韩欧美在线观看| av.在线天堂| 亚洲婷婷狠狠爱综合网| 99re6热这里在线精品视频| 国产高清国产精品国产三级| h视频一区二区三区| 国产精品人妻久久久影院| 免费黄频网站在线观看国产| 精品国产一区二区久久| 久久久久久人妻| 26uuu在线亚洲综合色| 在线亚洲精品国产二区图片欧美 | 精品一区二区三区视频在线| 高清午夜精品一区二区三区| 亚洲高清免费不卡视频| 女人久久www免费人成看片| 久久韩国三级中文字幕| 看十八女毛片水多多多| 日日爽夜夜爽网站| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 99re6热这里在线精品视频| 在线观看美女被高潮喷水网站| 热99久久久久精品小说推荐| 精品一品国产午夜福利视频| 七月丁香在线播放| 国产精品麻豆人妻色哟哟久久| 日韩中文字幕视频在线看片| 夜夜爽夜夜爽视频| 免费久久久久久久精品成人欧美视频 | 日本vs欧美在线观看视频| av在线app专区| av又黄又爽大尺度在线免费看| 在线观看一区二区三区激情| 亚洲内射少妇av| 熟女av电影| 午夜福利视频精品| 制服丝袜香蕉在线| 少妇熟女欧美另类| 最黄视频免费看| 亚洲美女搞黄在线观看| 99热网站在线观看| 亚洲第一av免费看| 一级a做视频免费观看| 亚洲国产精品999| 日韩成人av中文字幕在线观看| 亚洲成色77777| 九色亚洲精品在线播放| 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 爱豆传媒免费全集在线观看| 美女中出高潮动态图| 丰满乱子伦码专区| 亚洲性久久影院| 久久国产亚洲av麻豆专区| 成人亚洲欧美一区二区av| 夜夜骑夜夜射夜夜干| 男的添女的下面高潮视频| 麻豆乱淫一区二区| 国产午夜精品久久久久久一区二区三区| 久久久精品区二区三区| 亚洲色图 男人天堂 中文字幕 |