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

    呂宋海峽附近海域海表面高度季節(jié)內(nèi)變化

    2020-04-17 07:06:10王慶業(yè)
    海洋科學(xué) 2020年3期
    關(guān)鍵詞:斜壓呂宋海峽

    袁 欣, 王慶業(yè)

    呂宋海峽附近海域海表面高度季節(jié)內(nèi)變化

    袁 欣1, 2, 3, 王慶業(yè)1, 2

    (1. 中國(guó)科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國(guó)科學(xué)院 海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071; 3. 中國(guó)科學(xué)院大學(xué), 北京 100049)

    利用1993~2017年海表面高度異常數(shù)據(jù)集, 分析研究了西北太平洋季節(jié)內(nèi)變化(20~120 d)的整體分布特征, 結(jié)果表明空間上季節(jié)內(nèi)信號(hào)在20°N附近海域(16°~24°N)最強(qiáng), 時(shí)間上在6~8月達(dá)到一年中的最大值。在呂宋海峽東側(cè)(123.875°E, 20.125°N)季節(jié)內(nèi)信號(hào)周期(70 d)和傳播速度(10.7~12.7cm/s)均大于呂宋海峽西側(cè)(119.625°E, 20.125°N)(60 d, 6.5~7.8cm/s)。在大洋內(nèi)部(123°~140°E, 18°~24°N)存在準(zhǔn)90 d的周期信號(hào), 傳播速度約10.3cm/s。傳播路徑受黑潮的影響發(fā)生改變, 由沿緯度西傳轉(zhuǎn)向向西北方向傳播。第一斜壓Rossby波理論對(duì)海表面高度季節(jié)內(nèi)變化的周期和傳播速度具有很好的解釋性。

    西北太平洋; 呂宋海峽; 海表面高度; 季節(jié)內(nèi)變化

    西北太平洋具有復(fù)雜的海洋環(huán)流系統(tǒng)[1-4]。西向的北赤道流遇菲律賓東岸一分為二[5-6], 向北是典型的西邊界流的黑潮, 繼而構(gòu)成太平洋亞熱帶環(huán)流; 向南是棉蘭老海流, 與北赤道逆流一起構(gòu)成太平洋熱帶環(huán)流系統(tǒng)。西北太平洋海域環(huán)流在全球和區(qū)域海洋與氣候變異中具有重要作用[7]。

    近20年來, 西北太平洋海表面高度季節(jié)內(nèi)變化引起國(guó)內(nèi)外學(xué)者關(guān)注。劉秦玉等[8]利用衛(wèi)星高度計(jì)數(shù)據(jù)和渦分辨率的海洋環(huán)流模式數(shù)據(jù), 研究表明西北太平洋在5°N、10°N、20°N三個(gè)緯度帶海面高度分別存在30、60、90 d的準(zhǔn)周期振蕩, 后續(xù)胡瑞金[9]利用小波分析方法同樣也表明在15°N處和20°N存在準(zhǔn)60、90 d的振蕩, 并且均呈緯向分布[10-11]。

    關(guān)于海平面高度季節(jié)內(nèi)變化機(jī)制, 不同學(xué)者在不同海域得出的結(jié)論不同。如孟加拉灣海平面季節(jié)內(nèi)變化是由渦旋導(dǎo)致[12-13], 南海區(qū)域海表面高度季節(jié)內(nèi)變化是由海洋不穩(wěn)定和季節(jié)內(nèi)風(fēng)應(yīng)力強(qiáng)迫引起[14-15],印尼貫穿流區(qū)和北赤道流區(qū)上層海洋準(zhǔn)60 d振蕩的主要機(jī)制是斜壓不穩(wěn)定[10]。另外, Polito等[16]分析了全球海洋不同周期的海表面高度異常信號(hào)的相速度、周期、波長(zhǎng)和振幅等波動(dòng)要素, 結(jié)果表明Rossby波的相速度與線性自由Rossby波理論值吻合較好, 以及Lin等[17]和翟萍[18]將緯向分布的不同周期的海表面高度變化解釋為季節(jié)內(nèi)Rossby波。

    對(duì)于海表面高度變化周期呈緯向分布特征在此前眾多學(xué)者研究工作中也都有提及, 但同緯度的海表面高度季節(jié)內(nèi)變化是否在存在空間/時(shí)間差異?如果存在, 其物理現(xiàn)象又該如何解釋?本文利用25 a (1993~2017年)的高度計(jì)數(shù)據(jù)資料研究20°N附近緯度帶上, 呂宋海峽東西兩側(cè)海表面高度的季節(jié)內(nèi)變化差異, 并給出可能的解釋。

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

    1.1 數(shù)據(jù)

    本文所用的觀測(cè)數(shù)據(jù)是CLS (Collecte Localization Satellite)中心和CNES(Centre National d’Etudes Spatiales)中心的海平面異常場(chǎng)(Sea level anomaly, SLA) (https://climate.copernicus.eu/)。該數(shù)據(jù)融合了TOPEX/Poseidon、Jason-1、OSTM/Jason-2和Jason-3多種衛(wèi)星數(shù)據(jù)資料。選取數(shù)據(jù)時(shí)間范圍1993年1月1日~ 2017年12月31日, 空間范圍0~25°N, 115°~155°E, 時(shí)間分辨率是1 d, 空間分辨率為0.25°×0.25°的格點(diǎn)數(shù)據(jù)。

    此外, 還使用了Chelton等[19]根據(jù)美國(guó)國(guó)家海洋中心(National Oceanographic Data Center, NODC)提供的全球氣候平均態(tài)的1°×1°的溫鹽資料求出的第一斜壓Rossby波變形半徑和重力波速數(shù)據(jù)資料(http: // www-po.coas.oregonstate.edu/research/po/research/ros-sby_radius/)。

    1.2 方法及公式

    在提取季節(jié)內(nèi)信號(hào)分析時(shí)用到Lanczos濾波器[20], 時(shí)滯相關(guān)分析和經(jīng)驗(yàn)正交分解(Empirical Orthogonal Function, EOF)三種常規(guī)統(tǒng)計(jì)方法, 在周期分析用到功率譜分析(Power Spectral Density, PSD)。

    根據(jù)Rossby波的頻散關(guān)系, 其公式為:

    式中,是頻率,和分別是緯向和經(jīng)向方向上的波數(shù),1是第一斜壓Rossby波半徑,是科氏參數(shù)的經(jīng)向梯度。

    2 結(jié)果與分析

    2.1 季節(jié)內(nèi)信號(hào)處理

    本文中季節(jié)內(nèi)變化是指周期在20~120 d內(nèi)的信號(hào), 圖1是海表面高度季節(jié)內(nèi)變化(20~120 d帶通濾波后的海表面高度異常場(chǎng), SLA’)的標(biāo)準(zhǔn)差分布。季節(jié)內(nèi)變化在16°~24°N緯狀帶海域最強(qiáng)(標(biāo)準(zhǔn)差大于5 cm), 從呂宋海峽向東逐漸減弱, 核心區(qū)域位于臺(tái)灣東南海域, 最大值約9.3 cm(圖1)。

    圖1 20~120 d帶通濾波的海表面高度異常的標(biāo)準(zhǔn)差分布

    圖2是不同季節(jié)海表面高度季節(jié)內(nèi)變化的標(biāo)準(zhǔn)差。對(duì)比發(fā)現(xiàn)在6~8月季節(jié)內(nèi)海表面變化強(qiáng)度達(dá)到最大值(圖2c), 說明海表面高度的季節(jié)內(nèi)變化強(qiáng)度存在明顯的季節(jié)差異, 這可能與背景流場(chǎng)的季節(jié)變化導(dǎo)致的海洋不穩(wěn)定有關(guān)[21]。

    圖2 不同季節(jié)20~120 d帶通濾波的海表面高度異常的標(biāo)準(zhǔn)差分布

    注: a: 12~2月; b: 3~5月; c: 6~8月; d: 9~11月; A: 呂宋海峽西側(cè)季節(jié)內(nèi)海表面高度標(biāo)準(zhǔn)差最大值點(diǎn); B: 呂宋海峽東側(cè)季節(jié)內(nèi)海表面高度標(biāo)準(zhǔn)差最大值點(diǎn)

    在呂宋海峽西側(cè)季節(jié)內(nèi)海表面高度標(biāo)準(zhǔn)差相對(duì)較大(5 cm), 根據(jù)呂宋海峽東西兩側(cè)季節(jié)內(nèi)海表面高度標(biāo)準(zhǔn)差最大值的分布(圖2c), 選取A、B兩個(gè)點(diǎn)單獨(dú)研究, 其坐標(biāo)分別是119.625°E, 20.125°N和123.875°E, 20.125°N。圖3是A、B兩點(diǎn)處季節(jié)內(nèi)海表面高度變化的不同時(shí)間滯后相關(guān)性的分布圖, 經(jīng)過Radon換算計(jì)算A點(diǎn)(圖3a)的傳播速度約6.5~ 7.8 cm/s, B點(diǎn)(圖3b)的傳播速度較大, 約10.7~ 12.7 cm/s。為了進(jìn)一步驗(yàn)證呂宋海峽兩側(cè)季節(jié)內(nèi)信號(hào)周期的不同, 本文對(duì)A、B兩站點(diǎn)進(jìn)行了功率譜分析(圖4)。在A點(diǎn)(圖4a)和B(圖4b)點(diǎn)的季節(jié)內(nèi)海表面高度變化周期分別是62天和63、73天, 可以看成準(zhǔn)60 d和準(zhǔn)70 d周期。可見呂宋海峽東側(cè)(B站點(diǎn))海流季節(jié)內(nèi)信號(hào)周期大于西側(cè)(A站點(diǎn))周期。為揭示不同經(jīng)度上季節(jié)內(nèi)信號(hào)的分布我們給出20.125°N功率譜(圖5)??梢钥闯鲈?23°E東西兩側(cè)季節(jié)內(nèi)海表面高度的周期呈現(xiàn)出不同的分布狀態(tài)。123°E以西周期徘徊于62 d左右, 這也和之前眾多學(xué)者得出呂宋海峽附近存在的準(zhǔn)60 d周期一致[22]; 在123°E以東海域周期從75 d增長(zhǎng)到90 d左右, 也驗(yàn)證了上文的結(jié)果(圖4), 表明在20.125°N緯度上, 季節(jié)內(nèi)海表面高度變化的周期自東向西逐漸變小。

    圖3 A(a)、B(b)站點(diǎn)SLA’的時(shí)滯空間相關(guān)系數(shù)

    圖4 A(a)、B(b)站點(diǎn)的SLA’的功率譜(藍(lán)線)和95%置信區(qū)間(紅線)

    圖5 20.125°N上SLA’的功率譜

    2.2 EOF分析

    為研究整個(gè)海域季節(jié)內(nèi)變化的時(shí)空分布, 我們選取115°~123°E, 19°~22°N和123°~140°E, 18°~24°N兩矩形海域分別代表呂宋海峽附近海域、大洋內(nèi)部海域, 對(duì)季節(jié)內(nèi)海表面高度進(jìn)行經(jīng)驗(yàn)正交分解, 結(jié)果分別如圖6、圖7所示。在大洋內(nèi)部海域第一模態(tài)(圖7a)和第二模態(tài)的(圖7b)的空間模態(tài)相似, 呈經(jīng)向帶狀分布, 方差解釋率分別是16.7%和15.2%。功率譜分析表明第一模態(tài)和第二模態(tài)時(shí)間序列(圖7c、圖7d)的周期一致, 均是92 d(圖7e、圖7f), 這也符合之前眾多學(xué)者的研究結(jié)論[8-9, 16, 23]。為了更加直觀的看到兩模態(tài)對(duì)應(yīng)的標(biāo)準(zhǔn)化的時(shí)間序列的關(guān)系, 本文只觀察近五年時(shí)間序列(圖8a), 可以明顯看到二者存在相位差, 進(jìn)一步做時(shí)滯相關(guān)分析(圖8b), 當(dāng)?shù)诙B(tài)的時(shí)間序列滯后19 d時(shí), 二者相關(guān)性系數(shù)達(dá)到最大, 為0.89。為此本文猜測(cè)季節(jié)內(nèi)海表面高度經(jīng)過EOF分解后得到的第一模態(tài)和第二模態(tài)可能為同一個(gè)模態(tài)在前后兩個(gè)時(shí)間節(jié)點(diǎn)下的表征結(jié)果, 兩模態(tài)的空間模態(tài)分布(圖7a、圖7b)存在約1.625°的相位差, 結(jié)合19 d的滯后天數(shù), 計(jì)算得對(duì)應(yīng)的傳播速度為10.3 cm/s, 方向?yàn)橄蛭鱾鞑?。相比之下呂宋海峽附近海域第一、二模態(tài)的空間模態(tài)分布未曾顯出與大洋內(nèi)部海域一致的結(jié)果(圖6a、圖6b), 時(shí)間序列的功率譜分析也未呈現(xiàn)兩者相似的周期(圖6e、圖6f), 但是我們可以看到其第一模態(tài)時(shí)間序列的周期在62 d左右(圖6e), 這與上文在A站點(diǎn)處的結(jié)果較為一致(圖4a、圖5)。

    圖6 呂宋海峽附近海域的SLA’的EOF分析

    注: a. 第一模態(tài)的空間模態(tài), 方差解釋率19.8%; b. 第二模態(tài)的空間模態(tài), 方差解釋率17.5%; c. a對(duì)應(yīng)的時(shí)間序列; d. b對(duì)應(yīng)的時(shí)間序列; e: c對(duì)應(yīng)的功率譜分析; f. d對(duì)應(yīng)的功率譜分析

    圖7 大洋內(nèi)部的SLA’的EOF分析

    注: a. 第一模態(tài)的空間模態(tài), 方差解釋率16.7%; b. 第二模態(tài)的空間模態(tài), 方差解釋率15.2%; c. a對(duì)應(yīng)的時(shí)間序列; d. b對(duì)應(yīng)的時(shí)間序列; e: c對(duì)應(yīng)的時(shí)間序列的功率譜分析; f. d對(duì)應(yīng)的時(shí)間序列的功率譜分析

    圖8 圖7c、圖7d中的時(shí)間序列關(guān)系(a)以及時(shí)滯相關(guān)分析(b)

    2.3 路徑分析

    為了進(jìn)一步驗(yàn)證2.2中的猜測(cè), 本文計(jì)算了季節(jié)內(nèi)海表面高度的時(shí)滯分析相關(guān)性來探究其信號(hào)西傳的路徑。選取B點(diǎn)處進(jìn)行研究, 結(jié)合其63、73 d周期信號(hào)(圖4b), 分析不同時(shí)間超前滯后相關(guān)性的空間分布(圖9)。季節(jié)內(nèi)信號(hào)從超前36天(=-36 d)到超前0天(=0 d)沿緯度平行西傳(圖9a~圖9e), 在滯后0天(=0 d)到滯后36天(=36 d)沿西北方向傳播(圖9e~圖9i), 表明季節(jié)內(nèi)海表面高度變化信號(hào)從大洋內(nèi)部沿同緯度西傳, 到達(dá)黑潮流軸(呂宋海峽東側(cè))附近轉(zhuǎn)向西北方向, 推測(cè)與局地背景流場(chǎng)有關(guān)。

    結(jié)合在20.125°N周期處于50~100 d之間(圖5), 繼而對(duì)海表面高度進(jìn)行50~100 d的帶通濾波, 圖10顯示的是2014~2016年間的50~100 d帶通濾波后的海表面高度的時(shí)間-經(jīng)度分布圖, 藍(lán)線代表A站點(diǎn), 紅線代表B站點(diǎn), 可見西傳的季節(jié)內(nèi)海表面高度信號(hào)傳播路徑在A、B站點(diǎn)之間發(fā)生“間斷”。且速度發(fā)生改變(圖3), 因此我們有理由相信在20°N附近海域黑潮對(duì)于季節(jié)內(nèi)海表面高度信號(hào)西傳有一定影響, 位置約在123°E。

    圖9 B點(diǎn)處(黑色星號(hào))SLA’的時(shí)滯相關(guān)分布

    3 討論與結(jié)論

    3.1 物理機(jī)制討論

    根據(jù)自由Rossby波的頻散關(guān)系, 只考慮緯向傳播且在長(zhǎng)波近似下, 本文計(jì)算Rossby波西傳相速度(圖11a)顯示總體呈現(xiàn)緯向分布, 在海陸邊界和海底存在較大地形處發(fā)生彎曲變形。在20.125°N緯度帶上, 呂宋海峽東側(cè)第一斜壓Rossby波波速明顯大于西側(cè)速度, 分別是9.5 cm/s和6.5 cm/s, 與圖3具有一致的趨勢(shì)。根據(jù)公式(1), 當(dāng)=1/1時(shí)求得周期表明在大洋內(nèi)部20°N上具有90~100 d的理論周期(圖11b), 與上文計(jì)算求得大洋內(nèi)部準(zhǔn)90 d信號(hào)周期結(jié)果一致(圖7e、圖7f)。第一斜壓Rossby波的波速和周期能夠很好的解釋季節(jié)內(nèi)海表面高度傳播的速度和周期, 因此我們推斷海表面高度的季節(jié)內(nèi)變化主要是受第一斜壓Rossby波的影響。

    圖10 50~100 d帶通濾波后海表面高度異常的時(shí)間經(jīng)度分布

    注: 藍(lán)色線代表A點(diǎn), 紅色線代表B點(diǎn)

    對(duì)于第一斜壓Rossby波在呂宋海峽西側(cè)的周期(圖11b)與之前的計(jì)算(圖4、圖6)存在較大差異, 理論周期遠(yuǎn)大于本文計(jì)算的周期, 這可能是與呂宋海峽附近復(fù)雜的海域海洋動(dòng)力過程(如中尺度渦旋, 渦流相互作用[21, 24-26])有關(guān), 有待進(jìn)一步探究。

    3.2 結(jié)論

    利用20~120 d帶通濾波后的海表面高度數(shù)據(jù), 分別對(duì)研究海域內(nèi)單獨(dú)的站點(diǎn)(A, B)和整體海域(呂宋海峽附近海域, 大洋內(nèi)部海域)進(jìn)行分析, 本文得到以下結(jié)論:

    (1) 海表面高度的季節(jié)內(nèi)變化在以20°N為中心的緯度帶內(nèi)(16°~24°N)信號(hào)最強(qiáng), 自呂宋海峽向東逐漸減弱, 海表面高度的季節(jié)內(nèi)變化強(qiáng)度在6~8月達(dá)到一年中的最大值;

    (2) 海表面高度的季節(jié)內(nèi)變化在呂宋海峽兩側(cè)存在差異, 東側(cè)周期和傳播速度均大于西側(cè)。東側(cè)周期約70 d, 傳播速度約10.7~12.7 cm/s; 西側(cè)周期約60 d, 傳播速度6.5~7.8 cm/s。

    圖11 第一斜壓Rossby波波速(a)和周期(b)的分布

    (3) 在大洋內(nèi)部, 123°~140°E, 18°~24°N矩形海域內(nèi)海表面高度變化存在準(zhǔn)90 d的季節(jié)內(nèi)變化周期, 西傳速度約10.3 cm/s, 與第一斜壓Rossby波速基本一致。

    (4) 海表面高度的季節(jié)內(nèi)信號(hào)西傳路徑約在123°E處受到黑潮的影響發(fā)生改變, 由沿緯度西傳轉(zhuǎn)向向西北方向傳播。第一斜壓Rossby波對(duì)季節(jié)內(nèi)海表面高度變化的周期和傳播速度具有很好的解釋性。

    [1] Hu D X, Cui M C. The western boundary current of the Pacific and its role in the climate[J]. Chinese Journal of Oceanology and Limnology, 1991, 9(1): 1-14.

    [2] Fine R A, Lukas R, Bingham F M, et al. The Western Equatorial Pacific - a Water Mass Crossroads[J]. Journal of Geophysical Research-Oceans, 1994, 99(C12): 25063- 25080.

    [3] Lukas R, Yamagata T, Mccreary J P. Pacific low- latitude western boundary currents and the Indonesian throughflow[J]. Journal of Geophysical Research-Oceans, 1996, 101(C5): 12209-12216.

    [4] Hu D X, Wu L X, Cai W J, et al. Pacific western boun-dary currents and their roles in climate[J]. Nature, 2015, 522(7556): 299-308.

    [5] Qu T D, Lukas R. The bifurcation of the North Equatorial Current in the Pacific[J]. Journal of Physical Oceanogra-phy, 2003, 33(1): 5-18.

    [6] Wang Q Y, Cao R X, Zhang S W, et al. Bifurcation of Pacific North Equatorial Current at the surface[J]. Sciencein China Series D: Earth Sciences, 2009, 52(2): 227-231.

    [7] 胡石建, 胡敦欣. 西太平洋暖池研究綜述[J]. 海洋科學(xué)集刊, 2016, (51): 37-48. Hu Shijian, Hu Dunxin. Review on Western Pacific Warm Pool Study[J]. Studia Marina Sinica, 2016, (51): 37-48.

    [8] 劉秦玉, 王啟. 熱帶太平洋海平面高度季節(jié)內(nèi)振蕩的空間分布特征[J]. 青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 1999, 29(4): 549-555. Liu Qinyu, Wang Qi. Spatial distribution of the sea surface height Intraseasonal oscillation in the Tropical Pacific[J]. Periodical of Ocean University of China, 1999, 29(4): 549-555.

    [9] 胡瑞金, 劉秦玉. 熱帶太平洋海平面高度年變化與季節(jié)內(nèi)變化特征[J]. 海洋與湖沼, 2002, 33(3): 303-313. Hu Ruijin, Liu Qinyu. Annual and intraseasonal variation in sea surface height over the tropical Pacific[J]. Ocea-nologia et Limnologia Sinica, 2002, 33(3): 303- 313.

    [10] 劉秦玉, 潘愛軍, 劉征宇. 太平洋北赤道流區(qū)上層海洋季節(jié)內(nèi)振蕩及斜壓穩(wěn)定性[J]. 海洋與湖沼, 2003, 34(1): 94-100. Liu Qinyu, Pan Aijun, Liu Zhengyu. Intraseasonal oscillation and baroclinic instability of upper layer ocean in the North Equator Current[J]. Oceanologia et Limnologia Sinica, 2003, 34(1): 94-100.

    [11] Qiao F L, Tal E, Yuan Y L. Zonal distribution features of high frequency planetary waves in the oceans derivedfrom satellite altimeter data[J]. Acta Oceano-lo-gica Sinica, 2004, 23(1): 91-96.

    [12] Prasanna Kumar S, Nuncio M, Ramaiah N, et al. Eddy- mediated biological productivity in the Bay of Bengal during fall and spring intermonsoons[J]. Deep Sea Research Part I: Oceanographic Research Papers, 2007, 54(9): 1619-1640.

    [13] Cheng X H, Xie S P, Mccreary J P, et al. Intraseasonal variability of sea surface height in the Bay of Bengal[J]. Journal of Geophysical Research: Oceans, 2013, 118(2): 816-830.

    [14] Wu C R, Chang C W J. Interannual variability of the South China Sea in a data assimilation model[J]. Geophysical Research Letters, 2005, 32: L17611.

    [15] Zhuang W, Xie S P, Wang D X, et al. Intraseasonal variability in sea surface height over the South China Sea[J]. Journal of Geophysical Research-Oceans, 2010, 115: C04010.

    [16] Polito P S, Liu W T. Global characterization of Rossby waves at several spectral bands[J]. Journal of Geophy-sical Research-Oceans, 2003, 108(C1): 3018.

    [17] Lin X P, Yang J Y, Wu D X, et al. Explaining the global distri-bution of peak-spectrum variability of sea surface height[J]. Geophysical Research Letters, 2008, 35(14): L14602.

    [18] 翟萍. 全球海洋季節(jié)內(nèi)振蕩分布特征及機(jī)制研究[D]. 青島: 中國(guó)海洋大學(xué), 2008. Zhai Ping. The distribution characteristic and mechanism of intraseasonal oscillations in Global Oceans[D]. Qingdao: Ocean University of China, 2008.

    [19] Chelton D B, Deszoeke R A, Schlax M G, et al. Geogra-phical variability of the first baroclinic Rossby radius of deformation[J]. Journal of Physical Oceanography, 1998, 28(3): 433-460.

    [20] 姚菊香, 王盤興, 李麗平. 季節(jié)內(nèi)振蕩研究中兩種數(shù)字濾波器的性能對(duì)比[J]. 南京氣象學(xué)院學(xué)報(bào), 2005, 28(2): 248-253. Yao Juxiang, Wang Panxing, Li Liping. Performance contrast between two filters in Madden-Julian Oscillations analysis[J]. Transactions of Atmospheric Sciences, 2005, 28(2): 248-253.

    [21] Qiu B. Seasonal eddy field modulation of the North Pacific subtropical countercurrent: TOPEX/Poseidon observations and theory[J]. Journal of Physical Oceanography, 1999, 29(10): 2471-2486.

    [22] Zhang Z W, Zhao W, Tian J W, et al. Spatial structure and temporal variability of the zonal flow in the Luzon Strait[J]. Journal of Geophysical Research-Oceans, 2015, 120(2): 759-776.

    [23] 潘愛軍, 劉秦玉, 胡瑞金, 等. 北太平洋副熱帶逆流區(qū)海流的準(zhǔn)90天振蕩[J]. 青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2002, 32(1): 18-24. Pan Aijun, Liu Qinyu, Hu Ruijin, et al. Quasi 90-day oscillation of currents in the area of subtropical Countercurrent in the North Pacific[J]. Periodical of Ocean University of China, 2002, 32(1): 18-24.

    [24] Liu Y, Dong C M, Guan Y P, et al. Eddy analysis in the subtropical zonal band of the North Pacific Ocean[J]. Deep- Sea Research Part I-Oceanographic Research Papers, 2012, 68: 54-67.

    [25] Yang G, Wang F, Li Y L, et al. Mesoscale eddies in the northwestern subtropical Pacific Ocean: Statistical chara-cteristics and three-dimensional structures[J]. Journal of Geophysical Research-Oceans, 2013, 118(4): 1906-1925.

    [26] 律明坤, 臧楠, 王凡. 熱帶西太平洋北赤道逆流區(qū)渦旋統(tǒng)計(jì)分析[J]. 海洋科學(xué), 2017, 41(10): 67-76. Lv Mingkun, Zang Nan, Wang Fan. Statistical analysis of North Equatorial Countercurrent area in tropical we-stern Pacific[J]. Marine Sciences, 2017, 41(10): 67-76.

    Intraseasonal variability of sea surface height near the Luzon Strait

    YUAN Xin1, 2, 3, WANG Qing-ye1, 2

    (1. Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 2. Key Laboratory of Ocean Circulation and Wave, Chinese Academy of Sciences, Qingdao 266071, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China)

    Based on the sea surface height anomaly dataset from 1993 to 2017, the intraseasonal variability (ISV, period of 20~120 days) in the Northwestern Pacific was investigated. Results showed strong ISV occurred in the latitude band centered on 20°N (16°~24°N), and was stronger during June~August. In the eastern part of the Luzon Strait (point B, 123.875°E, 20.125°N), the signal period (70 days) and propagation velocity (10.7~12.7 cm/s) were larger than the corresponding values (60 days, 6.5~7.8 cm/s) on the western side of the Luzon Strait (point A, 119.625°E, 20.125°N). In the open ocean (123°~140°E, 18°~24°N), ISV had a typical period of 90 days and a typical propagation speed of about 10.3 cm/s. The propagation pathway of the ISV is associated with the Kuroshio Current. The ISV signal in the open ocean propagates westward, then turns northwestward when it is approaching the coast. Theory values of the first baroclinic Rossby wave phase period and velocity were consistent with those of the ISV calculated from sea surface height.

    Northwestern Pacific; Luzon Strait; sea surface height; intraseasonal variability

    Mar. 19, 2019

    P731.2

    A

    1000-3096(2020)03-0015-08

    10.11759/hykx20190319003

    2019-03-19;

    2019-05-08

    國(guó)家自然科學(xué)基金項(xiàng)目(41576014)

    [National Natural Science Foundation of China, No. 41576014]

    袁欣(1995-), 男, 河北石家莊人, 碩士在讀, 主要從事物理海洋學(xué)研究, 電話: 0532-82898963, Email: yuanxin17@mails.ucas.ac.cn; 王慶業(yè),通信作者, 研究員, 主要從事物理海洋學(xué)研究, 電話: 0532- 82896260, E-mail: wangqy@qdio.ac.cn

    (本文編輯: 李曉燕)

    猜你喜歡
    斜壓呂宋海峽
    譜松弛動(dòng)力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    海峽
    心聲歌刊(2021年6期)2021-02-16 01:12:32
    郁南縣河口鎮(zhèn)芒果喜獲豐收
    源流(2020年6期)2020-08-03 01:59:35
    以新聞史視角探析晚清《申報(bào)》中呂宋票廣告鏡像
    傳播力研究(2019年1期)2019-03-20 11:28:14
    斜壓渦度的變化與臺(tái)風(fēng)暴雨的關(guān)系研究
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    魔鬼海峽
    海峽可以很遠(yuǎn),但是心可以很近
    海峽姐妹(2016年2期)2016-02-27 15:15:43
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運(yùn)的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    渤海夏季第一斜壓羅斯貝變形半徑的計(jì)算與分析
    欧美成人一区二区免费高清观看| 日韩大尺度精品在线看网址| 亚洲国产欧洲综合997久久,| 在线免费观看的www视频| 久久久久久伊人网av| 日日啪夜夜撸| 91av网一区二区| 午夜免费成人在线视频| 日韩欧美国产一区二区入口| www日本黄色视频网| 丰满的人妻完整版| 国产91精品成人一区二区三区| 国产欧美日韩精品亚洲av| 国产精品女同一区二区软件 | 国产aⅴ精品一区二区三区波| 99精品久久久久人妻精品| 久久久国产成人精品二区| 淫妇啪啪啪对白视频| 一本精品99久久精品77| 麻豆一二三区av精品| 黄色日韩在线| 999久久久精品免费观看国产| www日本黄色视频网| 国产精品日韩av在线免费观看| 丝袜美腿在线中文| 真人一进一出gif抽搐免费| 日本撒尿小便嘘嘘汇集6| 99热这里只有精品一区| 国产探花在线观看一区二区| 国产亚洲精品av在线| 小说图片视频综合网站| 国产精品一区二区性色av| 少妇人妻一区二区三区视频| 国产男人的电影天堂91| 欧美日韩乱码在线| 久久久久九九精品影院| 麻豆国产av国片精品| 成年女人毛片免费观看观看9| 午夜福利成人在线免费观看| 人妻丰满熟妇av一区二区三区| 免费av观看视频| 波野结衣二区三区在线| 变态另类成人亚洲欧美熟女| 不卡一级毛片| 人人妻人人澡欧美一区二区| avwww免费| 久久精品国产清高在天天线| 波多野结衣高清作品| 深爱激情五月婷婷| 性插视频无遮挡在线免费观看| 亚洲七黄色美女视频| 久久午夜福利片| 18禁黄网站禁片午夜丰满| 日韩欧美在线二视频| 联通29元200g的流量卡| 日日摸夜夜添夜夜添av毛片 | 国产高清激情床上av| 精品国产三级普通话版| 国产精品av视频在线免费观看| 深夜精品福利| 欧美性感艳星| 最近在线观看免费完整版| 干丝袜人妻中文字幕| 老师上课跳d突然被开到最大视频| 久久久久久久久久成人| 日韩中文字幕欧美一区二区| 亚洲欧美日韩无卡精品| 亚洲国产欧洲综合997久久,| 亚洲av一区综合| 国产人妻一区二区三区在| 国产单亲对白刺激| 日本三级黄在线观看| 精品久久国产蜜桃| 十八禁网站免费在线| 国产私拍福利视频在线观看| 国产精品一及| 搡女人真爽免费视频火全软件 | 亚洲精品456在线播放app | 国产免费一级a男人的天堂| 一级黄色大片毛片| 欧美激情国产日韩精品一区| 小蜜桃在线观看免费完整版高清| 最近最新中文字幕大全电影3| 成人亚洲精品av一区二区| 精品无人区乱码1区二区| 可以在线观看的亚洲视频| 毛片女人毛片| 毛片一级片免费看久久久久 | 精品人妻熟女av久视频| 午夜激情福利司机影院| 国产人妻一区二区三区在| 99久久精品国产国产毛片| 精品久久久久久久末码| 国产一区二区亚洲精品在线观看| 久久香蕉精品热| 一级黄色大片毛片| 欧美不卡视频在线免费观看| 白带黄色成豆腐渣| 又黄又爽又免费观看的视频| 熟妇人妻久久中文字幕3abv| 联通29元200g的流量卡| 在线播放国产精品三级| 蜜桃久久精品国产亚洲av| 国产在视频线在精品| 中文在线观看免费www的网站| 国产高清三级在线| 久久久久久久久久久丰满 | 欧美国产日韩亚洲一区| 欧美3d第一页| 又爽又黄a免费视频| 极品教师在线视频| 国产精品一区二区免费欧美| 欧美日韩国产亚洲二区| 亚洲色图av天堂| 国产精品1区2区在线观看.| 国产女主播在线喷水免费视频网站 | 日本-黄色视频高清免费观看| 日韩欧美精品v在线| 国产蜜桃级精品一区二区三区| 黄色女人牲交| 在线免费观看的www视频| 久久精品影院6| 色哟哟·www| 九色成人免费人妻av| 色哟哟·www| 亚洲狠狠婷婷综合久久图片| 琪琪午夜伦伦电影理论片6080| 国产伦人伦偷精品视频| 麻豆av噜噜一区二区三区| 99热这里只有精品一区| 别揉我奶头 嗯啊视频| 最近最新中文字幕大全电影3| 又爽又黄a免费视频| av黄色大香蕉| 人妻丰满熟妇av一区二区三区| 亚洲中文日韩欧美视频| 欧美日韩国产亚洲二区| 久久久久久久久久成人| 亚洲va日本ⅴa欧美va伊人久久| 一个人看视频在线观看www免费| 又黄又爽又免费观看的视频| 床上黄色一级片| 国产精品三级大全| 麻豆国产97在线/欧美| 亚洲精品国产成人久久av| 精品久久国产蜜桃| 男女边吃奶边做爰视频| 国产一区二区三区在线臀色熟女| 岛国在线免费视频观看| 天天躁日日操中文字幕| 欧美zozozo另类| 欧美xxxx黑人xx丫x性爽| 亚洲最大成人手机在线| 一本一本综合久久| 国产麻豆成人av免费视频| 最近在线观看免费完整版| 三级毛片av免费| 日本撒尿小便嘘嘘汇集6| 成年女人毛片免费观看观看9| 十八禁国产超污无遮挡网站| 精品不卡国产一区二区三区| 精品人妻偷拍中文字幕| 国产一级毛片七仙女欲春2| 老熟妇仑乱视频hdxx| 日韩中字成人| 丰满人妻一区二区三区视频av| 国产亚洲91精品色在线| 哪里可以看免费的av片| 女生性感内裤真人,穿戴方法视频| 国产一区二区三区在线臀色熟女| av视频在线观看入口| 欧美zozozo另类| videossex国产| 一本久久中文字幕| 国产精华一区二区三区| 免费观看的影片在线观看| 在线免费十八禁| 美女高潮的动态| 九色国产91popny在线| 免费av观看视频| 欧美不卡视频在线免费观看| 午夜福利在线观看免费完整高清在 | 午夜福利高清视频| 少妇丰满av| 女生性感内裤真人,穿戴方法视频| 99久久无色码亚洲精品果冻| 欧美xxxx性猛交bbbb| 久久久久久久久久久丰满 | 成人av一区二区三区在线看| 婷婷色综合大香蕉| 久久人人爽人人爽人人片va| 69人妻影院| 久9热在线精品视频| 欧美丝袜亚洲另类 | 婷婷六月久久综合丁香| 亚洲自偷自拍三级| 国产精品国产三级国产av玫瑰| 丰满人妻一区二区三区视频av| 欧美日韩中文字幕国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 久久久久久久精品吃奶| 欧美3d第一页| 欧美一级a爱片免费观看看| 国产老妇女一区| 国产精品久久久久久久电影| 国产免费一级a男人的天堂| 自拍偷自拍亚洲精品老妇| 免费观看人在逋| 人妻少妇偷人精品九色| 赤兔流量卡办理| 国产精品久久久久久亚洲av鲁大| 久久久成人免费电影| 深夜精品福利| 国产精品爽爽va在线观看网站| 噜噜噜噜噜久久久久久91| 欧美在线一区亚洲| 毛片女人毛片| 午夜福利成人在线免费观看| 十八禁国产超污无遮挡网站| 少妇的逼水好多| 在线a可以看的网站| 男人舔奶头视频| 美女黄网站色视频| 俄罗斯特黄特色一大片| 99久久精品国产国产毛片| 日本黄色片子视频| 窝窝影院91人妻| 小说图片视频综合网站| 亚洲av中文av极速乱 | 日本欧美国产在线视频| 国产亚洲欧美98| 久久久久久久久久久丰满 | 免费观看在线日韩| 国产午夜精品久久久久久一区二区三区 | 久99久视频精品免费| 成年女人毛片免费观看观看9| 国产又黄又爽又无遮挡在线| 成人精品一区二区免费| 国产精品1区2区在线观看.| 亚洲专区中文字幕在线| 一a级毛片在线观看| АⅤ资源中文在线天堂| 亚洲欧美日韩卡通动漫| 五月玫瑰六月丁香| 国产精品三级大全| а√天堂www在线а√下载| 成年人黄色毛片网站| 欧美中文日本在线观看视频| 久久久久久久亚洲中文字幕| 在线观看舔阴道视频| 黄色日韩在线| av天堂在线播放| 精品一区二区免费观看| 日韩强制内射视频| 18禁裸乳无遮挡免费网站照片| 狂野欧美白嫩少妇大欣赏| 在线播放国产精品三级| 中文资源天堂在线| 日本黄色视频三级网站网址| 亚洲自拍偷在线| 亚洲18禁久久av| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区www在线观看 | 久久精品人妻少妇| 亚洲av.av天堂| 免费观看人在逋| 好男人在线观看高清免费视频| 91精品国产九色| a级一级毛片免费在线观看| 国产老妇女一区| 99热网站在线观看| 精品久久久久久久末码| 十八禁国产超污无遮挡网站| 我的老师免费观看完整版| 国产黄a三级三级三级人| 日韩亚洲欧美综合| 中文在线观看免费www的网站| 国产精品1区2区在线观看.| 欧美bdsm另类| 亚洲av日韩精品久久久久久密| 久久亚洲真实| 人妻久久中文字幕网| 亚洲国产日韩欧美精品在线观看| 夜夜看夜夜爽夜夜摸| 国产91精品成人一区二区三区| 最近最新免费中文字幕在线| 国产在视频线在精品| 国产精品久久久久久久电影| 久久久久久久久久黄片| 高清毛片免费观看视频网站| 又爽又黄a免费视频| 日本成人三级电影网站| 丰满的人妻完整版| av视频在线观看入口| 日本 欧美在线| 日韩欧美在线二视频| 婷婷精品国产亚洲av| 婷婷六月久久综合丁香| 欧美成人a在线观看| 亚洲天堂国产精品一区在线| 午夜福利在线观看免费完整高清在 | 国国产精品蜜臀av免费| 午夜精品在线福利| 麻豆精品久久久久久蜜桃| 中文字幕高清在线视频| 国产真实伦视频高清在线观看 | 免费看光身美女| 3wmmmm亚洲av在线观看| 日韩欧美免费精品| 色视频www国产| 欧美高清成人免费视频www| 一区二区三区激情视频| 亚洲美女视频黄频| 国产高清有码在线观看视频| 美女高潮的动态| av福利片在线观看| 五月玫瑰六月丁香| 亚洲av免费高清在线观看| 国产激情偷乱视频一区二区| 日韩欧美 国产精品| 在线免费观看的www视频| 亚洲四区av| netflix在线观看网站| 男人舔女人下体高潮全视频| 禁无遮挡网站| 国产激情偷乱视频一区二区| 亚洲成人精品中文字幕电影| 国产精品久久久久久av不卡| 精品久久久久久久久久免费视频| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 内地一区二区视频在线| 一区二区三区激情视频| 男女边吃奶边做爰视频| 成人鲁丝片一二三区免费| 欧美+日韩+精品| 日韩一本色道免费dvd| 国产一区二区激情短视频| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 我要看日韩黄色一级片| 免费在线观看成人毛片| 日韩,欧美,国产一区二区三区 | 午夜福利18| 噜噜噜噜噜久久久久久91| 午夜激情欧美在线| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 床上黄色一级片| 老司机午夜福利在线观看视频| 动漫黄色视频在线观看| 久久久久久国产a免费观看| 久久久国产成人精品二区| 国产精品1区2区在线观看.| 哪里可以看免费的av片| 性欧美人与动物交配| 免费av不卡在线播放| 国产一区二区在线av高清观看| 中出人妻视频一区二区| 一区二区三区激情视频| 国产精品爽爽va在线观看网站| 91麻豆av在线| 欧美三级亚洲精品| a级毛片免费高清观看在线播放| 欧美中文日本在线观看视频| 精品人妻偷拍中文字幕| 欧美成人性av电影在线观看| 国产色爽女视频免费观看| 亚洲无线观看免费| 久久这里只有精品中国| 久久久国产成人免费| 一个人看视频在线观看www免费| 88av欧美| 99热这里只有是精品在线观看| h日本视频在线播放| 99在线人妻在线中文字幕| 一级av片app| 99视频精品全部免费 在线| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产99精品国产亚洲性色| 色综合婷婷激情| 欧美日本亚洲视频在线播放| 一夜夜www| 欧洲精品卡2卡3卡4卡5卡区| 成年版毛片免费区| 精品久久久久久久久av| 少妇人妻精品综合一区二区 | 少妇被粗大猛烈的视频| 亚洲aⅴ乱码一区二区在线播放| 亚洲av成人精品一区久久| av在线观看视频网站免费| 嫁个100分男人电影在线观看| 国内久久婷婷六月综合欲色啪| 又黄又爽又免费观看的视频| 搡女人真爽免费视频火全软件 | 99久久中文字幕三级久久日本| 日本一二三区视频观看| 欧美一级a爱片免费观看看| 一区二区三区高清视频在线| 最近视频中文字幕2019在线8| 色视频www国产| 少妇高潮的动态图| 波多野结衣高清作品| 亚洲成a人片在线一区二区| 国产精品国产高清国产av| 国产精品一及| 可以在线观看的亚洲视频| 国产精品亚洲美女久久久| 日韩精品有码人妻一区| 美女cb高潮喷水在线观看| 国产精品野战在线观看| 欧美日本亚洲视频在线播放| 少妇猛男粗大的猛烈进出视频 | 99精品在免费线老司机午夜| 国产爱豆传媒在线观看| 免费观看人在逋| 蜜桃久久精品国产亚洲av| 999久久久精品免费观看国产| 亚洲无线在线观看| 极品教师在线视频| 欧美区成人在线视频| 99在线人妻在线中文字幕| 夜夜夜夜夜久久久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av成人av| 国产视频内射| 国产午夜福利久久久久久| 免费高清视频大片| 国产麻豆成人av免费视频| 伊人久久精品亚洲午夜| 九九久久精品国产亚洲av麻豆| 男人的好看免费观看在线视频| 窝窝影院91人妻| 精品人妻视频免费看| 午夜免费激情av| 色5月婷婷丁香| 美女高潮喷水抽搐中文字幕| 桃色一区二区三区在线观看| 亚洲中文字幕日韩| 精品免费久久久久久久清纯| 久久久久精品国产欧美久久久| 悠悠久久av| 日本a在线网址| 欧洲精品卡2卡3卡4卡5卡区| 99riav亚洲国产免费| 国产精品电影一区二区三区| 久久热精品热| 两个人视频免费观看高清| 大又大粗又爽又黄少妇毛片口| 亚洲欧美日韩无卡精品| 亚洲va在线va天堂va国产| 久久久久久久久久成人| 精品久久国产蜜桃| 乱码一卡2卡4卡精品| 日韩中字成人| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 亚洲精品成人久久久久久| 男人的好看免费观看在线视频| 真人做人爱边吃奶动态| 男人舔女人下体高潮全视频| 精品一区二区三区视频在线| 老司机深夜福利视频在线观看| 久久婷婷人人爽人人干人人爱| 99在线视频只有这里精品首页| 男人和女人高潮做爰伦理| 尾随美女入室| 精品午夜福利在线看| 一级av片app| ponron亚洲| 男女视频在线观看网站免费| 亚洲精华国产精华液的使用体验 | 可以在线观看的亚洲视频| 在线观看66精品国产| 超碰av人人做人人爽久久| 精品国内亚洲2022精品成人| 日韩欧美精品v在线| 久久国内精品自在自线图片| 国产免费男女视频| 88av欧美| 三级毛片av免费| 在线观看66精品国产| 亚洲狠狠婷婷综合久久图片| 欧美国产日韩亚洲一区| av天堂在线播放| 男人和女人高潮做爰伦理| 国产精品1区2区在线观看.| 午夜亚洲福利在线播放| 18+在线观看网站| 99在线人妻在线中文字幕| 国产综合懂色| 国产伦精品一区二区三区视频9| 成人鲁丝片一二三区免费| 成人国产一区最新在线观看| 18禁裸乳无遮挡免费网站照片| 白带黄色成豆腐渣| 国产黄色小视频在线观看| 欧美日本视频| 国国产精品蜜臀av免费| 人妻丰满熟妇av一区二区三区| 久久亚洲精品不卡| 久99久视频精品免费| 老熟妇乱子伦视频在线观看| 国产一区二区三区在线臀色熟女| 国产精品一区二区免费欧美| 午夜a级毛片| 欧美3d第一页| 国产乱人视频| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 婷婷精品国产亚洲av| 欧美性感艳星| 成人一区二区视频在线观看| 高清日韩中文字幕在线| 精品久久久久久成人av| 亚洲欧美日韩无卡精品| 乱系列少妇在线播放| 九九在线视频观看精品| 日本爱情动作片www.在线观看 | 天堂网av新在线| 男人舔女人下体高潮全视频| 欧美3d第一页| 国产男靠女视频免费网站| 久久午夜亚洲精品久久| 日本熟妇午夜| 国产精品日韩av在线免费观看| 亚洲不卡免费看| 国产 一区精品| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 日本欧美国产在线视频| 在线观看免费视频日本深夜| 欧美在线一区亚洲| 久久精品国产亚洲av涩爱 | 亚洲精华国产精华液的使用体验 | 欧美最黄视频在线播放免费| 日日撸夜夜添| 内地一区二区视频在线| 久久国产乱子免费精品| 嫩草影视91久久| 久久久国产成人免费| 成年女人看的毛片在线观看| 中国美女看黄片| 看免费成人av毛片| 制服丝袜大香蕉在线| 亚洲精华国产精华液的使用体验 | 男女那种视频在线观看| 日韩欧美一区二区三区在线观看| 如何舔出高潮| 久久人人精品亚洲av| 久久精品国产亚洲网站| 91午夜精品亚洲一区二区三区 | 99视频精品全部免费 在线| 99久久无色码亚洲精品果冻| 国产精品乱码一区二三区的特点| 18禁在线播放成人免费| 亚洲精华国产精华液的使用体验 | 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 乱码一卡2卡4卡精品| 久久99热这里只有精品18| 久久精品国产亚洲网站| 国产白丝娇喘喷水9色精品| 国内毛片毛片毛片毛片毛片| 国产探花在线观看一区二区| 黄片wwwwww| 91麻豆av在线| 精品久久久久久,| 一级黄片播放器| 美女免费视频网站| 精品99又大又爽又粗少妇毛片 | 一卡2卡三卡四卡精品乱码亚洲| 97人妻精品一区二区三区麻豆| 色视频www国产| 老熟妇仑乱视频hdxx| av女优亚洲男人天堂| av在线蜜桃| 少妇人妻精品综合一区二区 | 亚洲av第一区精品v没综合| 少妇高潮的动态图| 最近最新免费中文字幕在线| 看十八女毛片水多多多| 91在线精品国自产拍蜜月| 黄色配什么色好看| 啪啪无遮挡十八禁网站| 老熟妇乱子伦视频在线观看| 色精品久久人妻99蜜桃| 三级男女做爰猛烈吃奶摸视频| 日本 欧美在线| 日本黄色视频三级网站网址| 日本一二三区视频观看| 身体一侧抽搐| 亚洲av免费高清在线观看| 日韩欧美 国产精品| 国产人妻一区二区三区在| 久久热精品热| 国产精品伦人一区二区| 91av网一区二区| 在线播放国产精品三级| 99热网站在线观看| 国产aⅴ精品一区二区三区波| 国产av不卡久久| 在线观看舔阴道视频| 国产高清有码在线观看视频| 国内精品一区二区在线观看| 一区二区三区免费毛片| 免费观看在线日韩| 乱系列少妇在线播放| 村上凉子中文字幕在线| 国产单亲对白刺激| 久久精品夜夜夜夜夜久久蜜豆| 日本三级黄在线观看| 久久中文看片网|