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

    熱帶西太平洋北赤道逆流區(qū)渦旋統(tǒng)計分析

    2017-02-09 09:19:21律明坤
    海洋科學(xué) 2017年10期

    律明坤, 臧 楠, 王 凡, 4, 5

    ?

    熱帶西太平洋北赤道逆流區(qū)渦旋統(tǒng)計分析

    律明坤1, 2, 3, 臧 楠1, 2, 王 凡1, 2, 4, 5

    (1. 中國科學(xué)院海洋研究所, 山東 青島 266071; 2. 中國科學(xué)院海洋環(huán)流與波動重點(diǎn)實(shí)驗(yàn)室, 山東 青島 266071; 3. 中國科學(xué)院大學(xué), 北京 100049; 4. 青島海洋科學(xué)與技術(shù)國家實(shí)驗(yàn)室, 山東 青島 266071; 5. 海洋動力過程與氣候功能實(shí)驗(yàn)室, 山東 青島 266100)

    隨著海洋技術(shù)的發(fā)展, 中尺度的海洋過程越來越多的被揭示, 中尺度渦旋作為重要的中尺度海洋過程, 已經(jīng)被大量的研究。但對于熱帶西北太平洋海區(qū), 中尺度渦旋特征的空間分布、季節(jié)變化以及移動規(guī)律等方面的研究還有所欠缺。本文使用Chelton提供的渦旋數(shù)據(jù)集, 統(tǒng)計分析了熱帶西北太平洋海區(qū)渦旋特征的空間分布, 發(fā)現(xiàn)以往研究較少的北赤道逆流(North Equatorial Countercurrent, NECC)區(qū)(A海區(qū), 120°~180°E, 4°~6°N)較臨近海域生成渦旋數(shù)量更多, 渦旋半徑、振幅、生命周期及非線性強(qiáng)度更大, 移動距離更遠(yuǎn), 并且A海區(qū)渦旋經(jīng)向移動距離服從伽馬分布。渦旋在靠近西邊界的區(qū)域更易向南移動, 而在西邊界以東的區(qū)域更易向北移動。A海區(qū)渦旋的生成數(shù)量具有明顯的季節(jié)變化, 主要受到流場剪切強(qiáng)度的影響。同時ENSO會對該區(qū)渦旋生成產(chǎn)生影響, 其影響機(jī)制需要進(jìn)一步的研究。

    中尺度渦; 北赤道逆流; 特征分布; 季節(jié)變化; ENSO

    熱帶西北太平洋擁有著復(fù)雜的表層環(huán)流系統(tǒng)[1-3]。北赤道流(North Equatorial Current, NEC)在抵達(dá)菲律賓沿岸后分叉形成向北流動的黑潮(Kuroshio Current, KC)和向南流動的棉蘭老海流(Mindanao Current, MC)[2, 4-5]。大部分KC沿著中國東部陸坡向北流動, 小部分通過呂宋海峽入侵南海[6-7]。MC在向南流動的過程中在棉蘭老島東南部發(fā)生分叉, 大部分MC流入蘇拉威西海匯入印尼貫穿流(Indonesian Throughflow, ITF), 其余部分匯入北赤道逆流(North Equatorial Countercurrent, NECC)[5, 7-10](圖1, 背景底色和黑色矢量箭頭使用AVISO的1993年1月1日—2015年12月31日的海表面高度數(shù)據(jù)和地轉(zhuǎn)流數(shù)據(jù)繪制)。

    熱帶西北太平洋海區(qū)普遍存在著渦旋活動。其中有兩個重要的渦旋, 一個是棉蘭老島以東的氣旋式的棉蘭老渦(Mindanao Eddy, ME), 中心位于7°N, 129°~130°E附近(圖1)。另一個是反氣旋式的哈馬黑拉渦(Halmahera Eddy, HE), 它是一個季節(jié)性渦旋, 位于ME東南, 中心位于(131°E, 4°N)附近, 是連接新幾內(nèi)亞沿岸流(New Guinea Coastal Current, NGCC)和NECC的通道[3]。除了兩個大型的渦旋, 熱帶西北太平洋海區(qū)還存在著中尺度渦旋活動, Heron等[11]利用海洋表面流實(shí)時分析場數(shù)據(jù)(Ocean Surface Currents Analyses—Real-time, OSCAR), 實(shí)測TRITON (TRIangle Trans-Ocean buoy Network)浮標(biāo)資料, 海軍研究實(shí)驗(yàn)室分層海洋模式(Naval Research Laboratory Layered Ocean Model, NLOM)資料, 和1985—2000年31個航次的船載聲學(xué)多普勒海流計(Joint Archive for Shipboard ADCP Data, JASADCP)資料研究了帕勞(Palau)海域環(huán)流的季節(jié)變化, 發(fā)現(xiàn)了之前沒有被人注意到的季節(jié)性渦旋——帕勞渦、加羅琳渦和密克羅尼西亞渦。Chen等[12]利用衛(wèi)星高度計數(shù)據(jù)研究了棉蘭老穹頂(Mindanao Dome,MD)海區(qū)的中尺度渦, 發(fā)現(xiàn)此海區(qū)中尺度渦多具有非線性, 直徑小于300 km,由于地形阻礙多為局地產(chǎn)生。Chen等[13]計算熱帶西北太平洋的中尺度渦動能(Eddy Kinetic Energy, EKE)的季節(jié)變化, 發(fā)現(xiàn)此海區(qū)EKE夏強(qiáng)冬弱, 且局限在上層200 m內(nèi)。此外, 中尺度渦旋可以帶來水體的拉格朗日輸運(yùn)[14-15], 其緯向質(zhì)量輸運(yùn)可與風(fēng)驅(qū)和熱鹽環(huán)流相比擬[15]。

    在熱帶西北太平洋海區(qū), 大尺度環(huán)流除西邊界流外, 主要是緯向的NEC和NECC, 其經(jīng)向輸運(yùn)的量級要遠(yuǎn)小于緯向輸運(yùn)的量級。而中尺度渦引起的經(jīng)向輸運(yùn)就可能在這一海區(qū)的經(jīng)向輸運(yùn)中占有很大比例, 從而影響這一海區(qū)的溫鹽結(jié)構(gòu)和環(huán)流結(jié)構(gòu)。這些重要問題的解決必然要建立在了解熱帶西北太平洋海區(qū)渦旋的空間分布特征及其時間變化規(guī)律的基礎(chǔ)上, 但是目前這方面的統(tǒng)計分析工作還比較欠缺。因此, 本文使用Chelton的數(shù)據(jù), 詳細(xì)地分析了熱帶西北太平洋渦旋的主要統(tǒng)計特征, 并討論了渦旋數(shù)量的季節(jié)變化以及ENSO循環(huán)對渦旋生成數(shù)量的影響。

    圖1 熱帶西北太平洋氣候態(tài)流場及海表面高度

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

    本文使用的渦旋數(shù)據(jù)是Chelton通過AVISO (Archiving, Validation and Interpretation of Satellite Oceanographic data, https: //www.aviso.altimetry.fr)的海表面高度(sea surface height, SSH)數(shù)據(jù)識別的渦旋數(shù)據(jù)集, 數(shù)據(jù)集的時間范圍是1993年1月~2015年4月, 空間范圍是全球海洋(除去2°S~2°N)(http: //wombat. coas.oregonstate.edu/eddies/)。渦旋的產(chǎn)生伴隨著流場的輻聚輻散, 從而產(chǎn)生SSH的變化, 而流線近似和海表面高度異常的等值線重合, 那么最外側(cè)的等高線就可以被定義為渦旋邊界。Chelton就是基于這一事實(shí)識別出氣旋渦和反氣旋渦。本文的研究海區(qū)為120°~ 180°E, 2°~20°N, 將使用此海區(qū)內(nèi)全部的渦旋數(shù)據(jù)。

    本文還使用AVSIO的海表面地轉(zhuǎn)流流速數(shù)據(jù)(https: //www.aviso.altimetry.fr)來計算NECC和NEC的流軸。使用的地轉(zhuǎn)流數(shù)據(jù)的時間范圍為1993年1月1日~2014年12月30日, 空間范圍為120°~180°E, 2°~20°N, 空間分辨率為0.25°×0.25°。

    在文章中討論渦旋的生成位置時需要計算流軸的位置[16], 其計算公式為

    另外, 在研究中我們使用流場剪切強(qiáng)度和流場強(qiáng)度的概念來分析渦旋生成的季節(jié)變化與周圍環(huán)境的關(guān)系。流場強(qiáng)度的計算公式為

    它是通過流場對水體的輸運(yùn)量來定義的, 流場對水體的輸運(yùn)量越大意味著流場的強(qiáng)度越強(qiáng)。

    根據(jù)剪切的物理意義可認(rèn)為兩支海流的剪切強(qiáng)度等于兩支海流的流場強(qiáng)度之差, 所以海流的剪切強(qiáng)度計算公式為

    分別定義NECC, NEC, SEC的緯度范圍為1°~6°N, 6°~ 10°N, 7°S~1°N, 并由此計算三支海流之間的剪切強(qiáng)度。

    2 熱帶西北太平洋渦旋的空間分布與基本特征統(tǒng)計

    本文統(tǒng)計研究了熱帶西北太平洋海區(qū)1°×1°網(wǎng)格內(nèi)產(chǎn)生的渦旋數(shù)量、振幅、半徑、生命周期的空間分布及隨緯度的變化(圖2a~圖2h)。整個研究海域內(nèi)共有渦旋5 115個, 其中氣旋渦2 517個, 反氣旋渦2 598個(圖2a)。進(jìn)一步統(tǒng)計渦旋個數(shù)在不同緯度帶上的分布, 渦旋產(chǎn)生數(shù)量的兩個極大值分別位于北赤道逆流區(qū)(4°~6°N)和副熱帶逆流區(qū)(15°~17°N), 分別定義為A區(qū)和B區(qū)(圖1)。在Chelton渦旋數(shù)據(jù)集中, 渦旋的振幅定義為渦旋邊界處的基礎(chǔ)高度和渦旋中心SSH的極值之間的差值[17]。渦旋近似處于地轉(zhuǎn)平衡(即沿著流的方向, 右側(cè)的壓力高, 水位高)[17],因此對于氣旋渦(反氣旋渦), 渦旋中心SSH低于(高于)渦旋邊界的基礎(chǔ)高度。整個海區(qū)渦旋生成時的平均振幅為3.07 cm, 其中小于5 cm的占84.22%, 大于10 cm的占2.74 %。振幅較大的渦旋主要分布在流速較大海區(qū), 包括NECC西部、MC和KC海區(qū)(圖2c), 這主要是由于地轉(zhuǎn)近似下, 流速越大, 海表面高度起伏越大(即振幅越大)。渦旋半徑在南部較大, 北部較小, 12°N以南平均半徑為110.70 km, 12°N以北平均半徑為79.75 km(圖2(e))。這種渦旋半徑隨緯度的變化是由于科氏參數(shù)隨緯度變化引起的。隨著緯度的升高, 科氏參數(shù)變大, 羅斯貝變形半徑將變小[18], 相應(yīng)的渦旋的半徑也將減小。渦旋生命周期存在兩個極大值區(qū), 分別位于16oN以北和6°N以南, 平均生命周期分別為83.78 d和54.85 d。6°~16°N的渦旋的平均生命周期為53.80 d。

    圖2 1993年1月~2015年4月渦旋各變量的空間分布(a、c、e、g)及緯向平均值(b、d、f、h)

    熱帶西北太平洋生成渦旋的個數(shù)隨半徑、振幅和生命周期的變化也有明顯的規(guī)律(圖3)。氣旋渦和反氣旋渦表現(xiàn)出了相同的變化規(guī)律。氣旋渦和反氣旋渦數(shù)量皆隨振幅的增大而減小, 95.61%的氣旋渦振幅小于10 cm(圖3a), 98.97%的反氣旋渦振幅小于10 cm。氣旋渦最大半徑為329.23 km, 94.84%的氣旋渦半徑在30~170 km(圖3b); 反氣旋渦最大半徑為358.84 km, 95.43%的反氣旋渦半徑在30~170 km。90.84%的氣旋渦生命周期小于120 d, 91.34%的反氣旋渦生命周期小于120 d, 生命周期≥38 d的氣旋渦和反氣旋渦的數(shù)量隨著生命周期的增加而減少(圖3c)。

    圖3 渦旋振幅(a)半徑(b)以及生命周期(c)的頻率分布圖

    3 熱帶西太平洋北赤道逆流區(qū)渦旋的經(jīng)向移動與輸運(yùn)

    熱帶西北太平洋的大尺度環(huán)流除西邊界流之外, 主要是緯向流, 這種大尺度的環(huán)流的質(zhì)量、鹽度和熱量的輸運(yùn)都有過相應(yīng)的研究。而對于中尺度及更小尺度的輸運(yùn)研究的很少。更重要的一點(diǎn)是對于大尺度的環(huán)流, 緯向的輸運(yùn)要遠(yuǎn)遠(yuǎn)大于經(jīng)向的輸運(yùn)。那么, 中尺度渦旋的經(jīng)向輸運(yùn)就有可能對質(zhì)量熱量的經(jīng)向輸運(yùn)起到重要作用。另外, 這一海區(qū)鹽度的經(jīng)向梯度遠(yuǎn)遠(yuǎn)大于緯向梯度[19], 中尺度渦旋的經(jīng)向鹽度輸運(yùn)對這一海區(qū)的鹽度結(jié)構(gòu)很可能也有重要影響。為研究這些問題, 就需要對渦旋的經(jīng)向移動有一個初步的認(rèn)識。因此, 本節(jié)我們對渦旋經(jīng)向移動規(guī)律做了統(tǒng)計分析。

    渦旋因不穩(wěn)定而產(chǎn)生后, 因?yàn)槭芨鞣N條件的影響, 具有了不可預(yù)測性, 渦旋的消亡位置在經(jīng)向上可以看成是隨機(jī)的。因此在大樣本的前提條件下, 通過統(tǒng)計得到的頻率分布可以近似的作為概率分布。將A海區(qū)詳細(xì)劃分為A1海區(qū)(120°~140°E, 4°~6°N), A2海區(qū)(140°~160°E, 4°~6°N), A3海區(qū)(160°~180°E, 4°~6°N), 分別統(tǒng)計A1、A2、A3海區(qū)渦旋的經(jīng)向移動距離的頻率分布并對其進(jìn)行擬合, 擬合曲線皆通過K-S伽馬分布檢驗(yàn)(圖4)。三個海區(qū)的形狀參數(shù)分別為2.915 0、5.589 6、4.881 6, 尺度參數(shù)=0.834 0、0.588 9、0.723 4。經(jīng)過統(tǒng)計A1區(qū)共168個渦旋, 其中65.72%的渦旋向南移動; A2區(qū)共193個渦旋, 其中63.26%的渦旋向北移動; A3區(qū)共165個渦旋, 其中67.39%的渦旋向北移動。由年平均地轉(zhuǎn)流數(shù)據(jù)計算的A1、A2、A3三個海區(qū)的平均經(jīng)向流速分別為: –0.003 4 m/s, 0.016 5 m/s, 0.006 9 m/s。因此在A1海區(qū)中向南運(yùn)動的渦旋較多, 在A2, A3海區(qū)中向北移動的渦旋較多的現(xiàn)象是由背景流的平流作用引起的。在B海區(qū)產(chǎn)生的渦旋的經(jīng)向移動距離不符合任何常見的分布(圖略)。

    圖4 A1(a)、A2(b)和A3(c)區(qū)渦旋經(jīng)向移動距離的概率密度分布擬合圖

    移動距離的正值: 渦旋向北移動; 移動距離負(fù)值: 渦旋向南移動

    Positive value of movement distance: eddy moving northward; negative value of movement distance: eddy moving southward

    為了進(jìn)一步分析渦旋在不同緯度帶上的經(jīng)向移動距離的遠(yuǎn)近, 分別給出渦旋的經(jīng)向移動距離分布(圖5a)和經(jīng)向移動距離(絕對值)的緯向平均分布(圖5b),發(fā)現(xiàn)4o~6oN生成的渦旋普遍能夠在經(jīng)向上移動的更遠(yuǎn), 這與該海區(qū)的背景流有關(guān)。圖5c為該海區(qū)的AVISO地轉(zhuǎn)流經(jīng)向流速的緯向平均分布, 可以看出, 經(jīng)向地轉(zhuǎn)流在4o~6oN達(dá)到最強(qiáng), 強(qiáng)的背景流能把渦旋帶到更遠(yuǎn)的緯度位置。

    圖5 渦旋經(jīng)向移動距離分布(a)、距離(絕對值)的緯向平均分布(b)以及AVISO地轉(zhuǎn)流經(jīng)向流速(絕對值)的緯向平均分布(c)

    以往的研究表明, 并不是所有識別出的渦旋都能夠有效地輸運(yùn)水體, 渦旋的非線性作用使得渦旋可以攜帶水體移動[14-15]。渦旋能否進(jìn)行直接的質(zhì)量輸運(yùn)取決于渦旋是否具有非線性。渦旋的非線性可以用非線性度=/來估計[14], 其中為渦旋的最大平均旋轉(zhuǎn)速度(一般為渦旋邊界處的平均速度),為渦旋移動速度, 當(dāng)>1時, 認(rèn)為渦旋具有非線性。分別統(tǒng)計三個海區(qū)(120°~140°E, 140°~160°E, 160°~ 180°E)的渦旋非線性特性, 由渦旋平均非線性度和非線性渦旋所占比例隨緯度變化曲線可知, 三個海區(qū)渦旋的平均非線性度和非線性渦旋所占比例皆在10°N附近達(dá)到最小值(圖6a、圖6b)。隨著緯度的升高, 渦旋的非線性度和非線性渦旋所占比例在量值增加的同時發(fā)生收束現(xiàn)象, 不同經(jīng)度區(qū)間差異變小。而越往南, 渦旋的非線性度和非線性渦旋所占比例在不同經(jīng)度區(qū)間的差異變大, 10°N以南的三個區(qū)域非線性度分別為4.08、1.78和1.28, 非線性渦旋所占比例分別為0.73、0.48和0.36。因此, 在低緯度地區(qū), 越靠近西邊界渦旋的非線性越強(qiáng), 非線性渦旋比例越大。低緯度的西邊界區(qū)的非線性強(qiáng)度明顯高于其它海區(qū), 可能是由于此海區(qū)內(nèi)流場和風(fēng)場復(fù)雜多變, 渦旋和海流之間, 渦旋和風(fēng)場之間的非線性相互作用更加頻繁而劇烈導(dǎo)致的。

    4 北赤道逆流區(qū)(A海區(qū))中渦旋生成數(shù)量的季節(jié)變化及其與副熱帶逆流區(qū)(B海區(qū))的比較

    由渦旋的空間分布特征得知A、B海區(qū)(圖2a) 生成渦旋較多, 進(jìn)一步分析這兩個海區(qū)內(nèi)生成渦旋的季節(jié)變化。統(tǒng)計得到渦旋生成個數(shù)的氣候態(tài)月平均曲線(圖7), A海區(qū)的氣旋渦和反氣旋渦均具有明顯的季節(jié)變化。氣旋渦夏季和冬季多, 春季和秋季少(圖7b)。而反氣旋渦在夏季數(shù)量最少, 冬季數(shù)量最多(圖7a)。B海區(qū)氣旋渦和反氣旋渦均沒有明顯的季節(jié)變化(圖8a、圖8b)。

    圖6 三個區(qū)域(120°~140°E, 140°~160°E, 160°~180°E)渦旋非線性度(a)、非線性渦旋(U/c>1, b)所占比例的緯向平均分布

    圖7 A區(qū)氣旋渦(a)反氣旋渦(b)數(shù)量的氣候態(tài)月平均分布

    圖8 B區(qū)(a)氣旋渦(b)反氣旋渦數(shù)量的氣候態(tài)月平均分布

    為了解釋A區(qū)渦旋生成數(shù)量季節(jié)變化的成因, 我們統(tǒng)計了NECC流軸附近不同區(qū)域的渦旋生成數(shù)量(表1), 其中流軸位置通過公式(1)計算。NECC是東向流, NEC和SEC是西向流, 如果渦旋是海流的剪切作用產(chǎn)生的, 那么在NECC以北會多出現(xiàn)氣旋, NECC以南會多出現(xiàn)反氣旋。經(jīng)統(tǒng)計實(shí)際情況也是如此(表1), 氣旋渦大部分生成在了NECC以北, 反氣旋渦大部分生成在了NECC以南。

    表1 NECC流軸附近不同區(qū)域的渦旋數(shù)量(個)

    圖9 氣候態(tài)月平均的NEC和NECC(a)、NECC和SEC(b)的流場剪切強(qiáng)度變化

    5 北赤道逆流區(qū)(A海區(qū))渦旋生成與氣候異常事件的關(guān)聯(lián)

    根據(jù)氣候預(yù)測中心(Climate Prediction Center, http: //www.cpc.ncep.noaa.gov/data/indices)發(fā)布的Nino3.4指數(shù), 定義Nino3.4指數(shù)≥0.5并持續(xù)5個月的事件為一次厄爾尼諾事件, Nino3.4指數(shù)≤–0.5并持續(xù)5個月的事件為一次拉尼娜事件。則1993年1月~2014年12月一共發(fā)生了6次厄爾尼諾事件(分別為: 1994~1995、1997~1998、2002~2003、2004~2005、2006~2007、2009~2010), 一共發(fā)生了5次拉尼娜事件(分別為: 1995~1996、1998~1999、2000~2001、2007~2008、2010~2012)。我們將所有厄爾尼諾事件和拉尼娜事件中生成的渦旋個數(shù)和對應(yīng)的Nino3.4指數(shù)分別合成序列, 得到圖10??梢钥闯鯝的三個海區(qū)氣旋渦大部分產(chǎn)生在厄爾尼諾成熟期及成熟期之前, 反氣旋渦大部分產(chǎn)生在厄爾尼諾成熟期和成熟期之后, 而在拉尼娜事件中并沒有明顯的規(guī)律。

    為了分析渦旋生成個數(shù)在厄爾尼諾事件中的變化原因, 同樣合成得到流場剪切強(qiáng)度在厄爾尼諾事件中的變化曲線(圖11)??梢? NECC和NEC剪切強(qiáng)度的變化以及NECC和SEC剪切強(qiáng)度的變化滯后于Nino3.4指數(shù)。而NECC和NEC的剪切產(chǎn)生氣旋渦, NECC和SEC的剪切產(chǎn)生反氣旋渦。因此, 可以通過NECC和SEC的剪切強(qiáng)度的變化滯后Nino3.4指數(shù)的變化來解釋反氣旋渦滯后厄爾尼諾成熟期大量生成的現(xiàn)象。而氣旋渦超前于厄爾尼諾成熟期大量生成的現(xiàn)象則不能單獨(dú)通過NECC和NEC的剪切作用解釋。因此, 推測在年際尺度上其他影響氣旋渦生成的因素(比如風(fēng)場)的作用將會顯現(xiàn), 而這有待進(jìn)一步研究。

    圖10 A海區(qū)氣旋渦在厄爾尼諾事件(a)以及拉尼娜事件(b)中的數(shù)量變化, 反氣旋渦在厄爾尼諾事件(c)以及拉尼娜事件(d)中數(shù)量的變化

    圖11 A海區(qū)NECC和NEC的流場剪切強(qiáng)度在厄爾尼諾事件中的變化(a)以及NECC和SEC的流場剪切強(qiáng)度在厄爾尼諾事件中的變化(b)

    6 總結(jié)與討論

    本文使用Chelton提供的渦旋數(shù)據(jù)集, 統(tǒng)計分析了熱帶西北太平洋海區(qū)渦旋特征的空間分布, 發(fā)現(xiàn)以往研究較少的NECC區(qū)(A海區(qū), 120°~180°E, 4°~6°N)較臨近海域生成渦旋數(shù)量更多, 渦旋半徑、振幅、生命周期及非線性強(qiáng)度更大, 移動距離更遠(yuǎn), 并且A海區(qū)渦旋經(jīng)向移動距離服從伽馬分布。渦旋在靠近西邊界的區(qū)域更易向南移動, 而在西邊界以東的區(qū)域更易向北移動, 這是由此海域背景流場的平流作用引起的。A海區(qū)氣旋渦和反氣旋渦的生成數(shù)量存在明顯的季節(jié)變化: 氣旋渦夏季和冬季多, 春季和秋季少; 而反氣旋渦在夏季數(shù)目少, 冬季數(shù)目多。B海區(qū)生成的渦旋數(shù)量沒有明顯的季節(jié)變化。A區(qū)氣旋渦和反氣旋渦的季節(jié)變化分別由NECC與NEC、以及NECC和SEC的流場剪切的季節(jié)變化引起。分別將研究期間的厄爾尼諾事件和拉尼娜事件中的渦旋數(shù)量和Nino3.4指數(shù)進(jìn)行合成, 對比發(fā)現(xiàn)在厄爾尼諾期間, 三個海區(qū)氣旋渦超前于厄爾尼諾成熟期大量產(chǎn)生, 反氣旋渦滯后于厄爾尼諾成熟期大量產(chǎn)生, 而在拉尼娜期間并未發(fā)現(xiàn)明顯的規(guī)律。通過分析發(fā)現(xiàn), NECC和NEC剪切強(qiáng)度的變化以及NECC和SEC剪切強(qiáng)度的變化均滯后于Nino3.4指數(shù)的變化。反氣旋渦滯后于厄爾尼諾成熟期的大量產(chǎn)生可以通過NECC和SEC的剪切強(qiáng)度峰值滯后于Nino3.4指數(shù)峰值來解釋, 而氣旋渦超前于厄爾尼諾成熟期的大量產(chǎn)生則可能由于受到其他因素(比如風(fēng)場)的影響, 這有待進(jìn)一步研究。

    [1] Fine R A, Lukas R, Bingham F M, et al. The western equatorial Pacific: A water mass crossroads [J]. Journal of Geophysical Research Atmospheres, 1994, 99(C12): 25063-25080.

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

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

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

    [5] Wang F, Song L, Li Y, et al. Semi‐annually alternating exchange of intermediate waters east of the Philippines[J]. Geophysical Research Letters, 2016, 43(13): 7059-7065.

    [6] Lee S K. On the structure of supercritical western boundary currents[J]. Dynamics of Atmospheres & Oceans, 2001, 33(4): 303-319.

    [7] Wyrtki K. The thermohaline circulation in relation to the general circulation in the oceans[J]. Deep Sea Research, 1961, 8(1): 39-64.

    [8] Lukas R, Firing E, Hacker P, et al. Observations of the Mindanao Current during the western equatorial Pacific Ocean circulation study[J]. Journal of Geophysical Research Atmospheres, 1991, 96(C4): 7089-7104.

    [9] Zhao J, Li Y, Wang F. Dynamical responses of the west Pacific North Equatorial Countercurrent (NECC) system to El Ni?o events[J]. Journal of Geophysical Research Oceans, 2013, 118(6): 2828-2844.

    [10] Zhao J, Li Y, Wang F. Seasonal variation of the surface North Equatorial Countercurrent (NECC) in the western Pacific Ocean[J]. Chinese Journal of Oceanology and Limnology, 2016, 34(6): 1332-1346.

    [11] Heron S F, Metzger E J, Skirving W J. Seasonal variations of the ocean surface circulation in the vicinity of Palau[J]. Journal of Oceanography, 2006, 62(4): 413-426.

    [12] Chen L, Jia Y, Liu Q. Mesoscale eddies in the Mindanao Dome region[J]. Journal of Oceanography, 2014, 71(1): 133-140.

    [13] Chen X, Qiu B, Chen S, et al. Seasonal eddy kinetic energy modulations along the North Equatorial Countercurrent in the western Pacific[J]. Journal of Geophysical Research Oceans, 2015, 120(9): 6351-6362.

    [14] Early J J, Samelson R M, Chelton D B. The Evolution and Propagation of Quasigeostrophic Ocean Eddies[J]. Journal of Physical Oceanography, 2011, 41(8): 1535- 1555.

    [15] Zhang Z, Wang W, Qiu B. Oceanic mass transport by mesoscale eddies[J]. Science, 2014, 345(6194): 322-324.

    [16] Hsin Y C, Qiu B. Seasonal fluctuations of the surface North Equatorial Countercurrent (NECC) across the Pacific basin[J]. Journal of Geophysical Research Oceans, 2012, 117(C6): 6001.

    [17] Chelton D B, Schlax M G, Samelson R M. Global observations of nonlinear mesoscale eddies[J]. Progress in Oceanography, 2011, 91(2): 167-216.

    [18] Pedlosky J. Geophysical fluid dynamics[M]. 2nd ed. New York: Springer-Verlag, 1987: 710.

    [19] Melnichenko O, Hacker P, Maximenko N, et al. Optimum interpolation analysis of Aquarius sea surface salinity[J]. Journal of Geophysical Research Oceans, 2016, 121(1): 602-616.

    (本文編輯: 李曉燕)

    Statistical analysis of North Equatorial Countercurrent area in tropical western Pacific

    Lü Ming-kun1, 2, 3, ZANG Nan1, 2, WANG Fan1, 2, 4, 5

    (1. Institute of Oceanology, the 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; 4. Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China; 5. Marine Dynamics Process and Climate Functional Laboratory, Qingdao 266100, China)

    With the development of ocean technology, more and more mesoscale ocean processes have been revealed. Mesoscale eddies are an important mesoscale ocean process that has been the subject of much research in recent decades. However, research on the spatial distribution, seasonal variation, and movement characteristics of mesoscale eddies in the northwestern Pacific is lacking. In this paper, we statistically analyzed the spatial distribution of eddies in the northwestern tropical Pacific using the eddy dataset provided by Chelton. We found that more eddies are generated in the North Equatorial Countercurrent (NECC) area (A zone, 120°–180°E, 4°–6°N) where less research effort has been focused. Moreover, the eddy radii, amplitudes, lifecycles, and nonlinear intensities are larger, the movement distances are farther, and the meridional eddy movement distances obey a gamma distribution in the A zone. Most eddies move southward in the region close to the western boundary, and most move northward far away from the western boundary. The quantity of eddies generated in the A zone exhibits obvious seasonal variation, which is mainly affected by the shear intensity of the currents. At the same time, eddy formation is affected by the ENSO in this area, and its mechanism requires further study.

    mesoscale eddies; North Equatorial Countercurrent; characteristic distribution; seasonal variation; ENSO

    Aug. 5, 2017

    P731.2

    A

    1000-3096(2017)10-0067-10

    10.11759/hykx20170508002

    2017-05-08;

    2017-06-17

    國家自然科學(xué)基金項(xiàng)目(41406014); 中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項(xiàng)(XDA11010204); “全球變化與海氣相互作用”專項(xiàng)(GASI-03- 01-01-05); 山東省科技發(fā)展計劃(2013GRC31503)

    [National Natural Science Foundation of China, No. 41406014; the Strategic Priority Research Program of Chinese Academy of Sciences, No. XDA11010201; the Global Change and Air-Sea Interaction Program, No. GASI-03-01-01-05; the Scientific and Technological Development Plan Project of Shandong Province, No. 2013GRC31503]

    律明坤(1990-), 男, 山東膠州人, 碩士在讀, 主要從事中尺度渦旋研究, Email: qdlvmingkun@126.com; 王凡,通信作者, 研究員, 主要從事物理海洋學(xué)研究, Email: fwang@qdio.ac.cn

    国产精品,欧美在线| av超薄肉色丝袜交足视频| 精品熟女少妇八av免费久了| 久久午夜综合久久蜜桃| 亚洲国产精品合色在线| 日本成人三级电影网站| 久久国产精品男人的天堂亚洲| 琪琪午夜伦伦电影理论片6080| 最近在线观看免费完整版| 久久热在线av| 妹子高潮喷水视频| 亚洲av成人av| 丝袜人妻中文字幕| 99热这里只有精品一区 | 国产一区二区三区在线臀色熟女| 欧美日韩一级在线毛片| 亚洲欧美日韩高清在线视频| 夜夜看夜夜爽夜夜摸| 久久天堂一区二区三区四区| 一个人观看的视频www高清免费观看 | 亚洲精品一区av在线观看| 欧美黑人精品巨大| 99国产综合亚洲精品| 国内精品久久久久久久电影| 色综合婷婷激情| 国内少妇人妻偷人精品xxx网站 | 巨乳人妻的诱惑在线观看| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| 最近最新中文字幕大全电影3 | 色综合站精品国产| 欧美日韩一级在线毛片| 国产aⅴ精品一区二区三区波| 美国免费a级毛片| 老汉色∧v一级毛片| 99热只有精品国产| 久久久久九九精品影院| 校园春色视频在线观看| 欧美 亚洲 国产 日韩一| 999久久久国产精品视频| АⅤ资源中文在线天堂| tocl精华| 国产在线精品亚洲第一网站| 久久久久久免费高清国产稀缺| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全免费视频| 亚洲第一欧美日韩一区二区三区| 一本大道久久a久久精品| 神马国产精品三级电影在线观看 | 亚洲色图 男人天堂 中文字幕| 亚洲国产精品久久男人天堂| 侵犯人妻中文字幕一二三四区| 国产不卡一卡二| 欧美日韩乱码在线| 成人三级做爰电影| 搞女人的毛片| 欧美av亚洲av综合av国产av| 国产成人精品无人区| 人人妻人人澡人人看| 免费一级毛片在线播放高清视频| 精品人妻1区二区| 亚洲狠狠婷婷综合久久图片| 一进一出抽搐动态| 欧美日本视频| 亚洲精品在线观看二区| 国产精品自产拍在线观看55亚洲| 夜夜夜夜夜久久久久| 亚洲av成人一区二区三| 午夜福利18| 日韩欧美免费精品| 日本撒尿小便嘘嘘汇集6| 精品国产国语对白av| 老熟妇仑乱视频hdxx| 久久中文看片网| 一本久久中文字幕| 久久国产亚洲av麻豆专区| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 免费在线观看黄色视频的| 啦啦啦观看免费观看视频高清| 亚洲三区欧美一区| 亚洲av中文字字幕乱码综合 | 久久久精品国产亚洲av高清涩受| 国产熟女xx| 久久久久亚洲av毛片大全| 欧美黑人精品巨大| 日韩大码丰满熟妇| 级片在线观看| 在线观看免费视频日本深夜| 国产精品一区二区三区四区久久 | 又紧又爽又黄一区二区| 国产欧美日韩精品亚洲av| 久久香蕉激情| 久久精品国产综合久久久| 欧美又色又爽又黄视频| 国产一区二区三区在线臀色熟女| 欧美中文日本在线观看视频| 757午夜福利合集在线观看| 一区福利在线观看| 俺也久久电影网| 制服人妻中文乱码| 日本免费一区二区三区高清不卡| 淫秽高清视频在线观看| 精品午夜福利视频在线观看一区| 一级a爱视频在线免费观看| 久热这里只有精品99| 母亲3免费完整高清在线观看| 国产精品久久电影中文字幕| 午夜两性在线视频| 亚洲精品粉嫩美女一区| 在线av久久热| av免费在线观看网站| 亚洲国产看品久久| 免费女性裸体啪啪无遮挡网站| xxxwww97欧美| www日本在线高清视频| 亚洲av成人不卡在线观看播放网| 亚洲av熟女| 亚洲avbb在线观看| 亚洲中文字幕日韩| av有码第一页| 首页视频小说图片口味搜索| 久久国产精品影院| 日韩三级视频一区二区三区| 成人国语在线视频| 国产视频一区二区在线看| 亚洲国产欧洲综合997久久, | 男人舔女人下体高潮全视频| 91国产中文字幕| 怎么达到女性高潮| 欧美av亚洲av综合av国产av| 又大又爽又粗| 国产1区2区3区精品| www.自偷自拍.com| 色精品久久人妻99蜜桃| 精品国内亚洲2022精品成人| 91成人精品电影| 久久久久久久精品吃奶| 高潮久久久久久久久久久不卡| 日韩欧美 国产精品| 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 国产精品亚洲一级av第二区| 女性生殖器流出的白浆| 亚洲精品国产一区二区精华液| 国产乱人伦免费视频| 草草在线视频免费看| 国产成人影院久久av| 精品第一国产精品| 女人高潮潮喷娇喘18禁视频| 性色av乱码一区二区三区2| 首页视频小说图片口味搜索| 国语自产精品视频在线第100页| 亚洲成人国产一区在线观看| 美女 人体艺术 gogo| 国产熟女xx| 亚洲片人在线观看| tocl精华| 成人永久免费在线观看视频| 宅男免费午夜| 国产真实乱freesex| 在线国产一区二区在线| 国产精品久久久人人做人人爽| 大香蕉久久成人网| 久久久久久久久中文| 97人妻精品一区二区三区麻豆 | 香蕉久久夜色| 国产精品久久久av美女十八| 国产亚洲av嫩草精品影院| 波多野结衣巨乳人妻| 亚洲欧美日韩无卡精品| 丝袜人妻中文字幕| 国产伦在线观看视频一区| 欧美黑人精品巨大| 99久久精品国产亚洲精品| 国产精品国产高清国产av| 中文字幕人妻丝袜一区二区| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 欧美激情极品国产一区二区三区| 日本三级黄在线观看| 悠悠久久av| 国产精品精品国产色婷婷| 亚洲avbb在线观看| 97碰自拍视频| 国产精品 国内视频| 可以在线观看毛片的网站| 久久人人精品亚洲av| 一个人观看的视频www高清免费观看 | av在线天堂中文字幕| 午夜福利成人在线免费观看| 国产高清视频在线播放一区| 亚洲一码二码三码区别大吗| 中文字幕久久专区| 草草在线视频免费看| 在线观看日韩欧美| 老司机在亚洲福利影院| 亚洲av日韩精品久久久久久密| 久久精品国产99精品国产亚洲性色| 亚洲成人免费电影在线观看| 久久久久久大精品| 成人18禁高潮啪啪吃奶动态图| 日韩欧美 国产精品| 精品免费久久久久久久清纯| 男男h啪啪无遮挡| 午夜a级毛片| 日韩高清综合在线| 天堂影院成人在线观看| 成人国产综合亚洲| 大型黄色视频在线免费观看| 久久亚洲精品不卡| 亚洲av第一区精品v没综合| 日本一区二区免费在线视频| 欧美久久黑人一区二区| 久久久久久免费高清国产稀缺| 99re在线观看精品视频| 国产精品爽爽va在线观看网站 | 精品无人区乱码1区二区| 色综合欧美亚洲国产小说| 50天的宝宝边吃奶边哭怎么回事| 搞女人的毛片| 免费在线观看影片大全网站| 国产精品香港三级国产av潘金莲| 久久天堂一区二区三区四区| 欧美激情高清一区二区三区| 给我免费播放毛片高清在线观看| www.熟女人妻精品国产| 亚洲自偷自拍图片 自拍| 国产在线观看jvid| 男人操女人黄网站| 亚洲色图av天堂| 女性生殖器流出的白浆| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av在线| 丝袜人妻中文字幕| 曰老女人黄片| 亚洲狠狠婷婷综合久久图片| 免费观看精品视频网站| 国产av一区二区精品久久| 真人一进一出gif抽搐免费| 亚洲专区字幕在线| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 亚洲熟女毛片儿| 久久中文字幕一级| 我的亚洲天堂| 国产精品久久久人人做人人爽| 97超级碰碰碰精品色视频在线观看| 亚洲三区欧美一区| 久久亚洲真实| 国产精品乱码一区二三区的特点| 长腿黑丝高跟| 美女扒开内裤让男人捅视频| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡| 黑人操中国人逼视频| 久久狼人影院| bbb黄色大片| 欧美日韩精品网址| 90打野战视频偷拍视频| 人妻丰满熟妇av一区二区三区| 亚洲国产看品久久| 欧美黄色淫秽网站| 欧美+亚洲+日韩+国产| 亚洲精品久久国产高清桃花| 色婷婷久久久亚洲欧美| 99国产综合亚洲精品| 精品久久久久久久久久免费视频| 久久99热这里只有精品18| 首页视频小说图片口味搜索| 亚洲五月色婷婷综合| 一二三四社区在线视频社区8| 欧美色视频一区免费| 精品久久久久久久久久久久久 | 久热这里只有精品99| 后天国语完整版免费观看| 一本精品99久久精品77| 欧美性长视频在线观看| 97人妻精品一区二区三区麻豆 | 热99re8久久精品国产| 人人澡人人妻人| 国产成人av激情在线播放| 女同久久另类99精品国产91| 午夜免费成人在线视频| 国产欧美日韩一区二区三| 国产精品,欧美在线| 成人免费观看视频高清| 欧美人与性动交α欧美精品济南到| 亚洲天堂国产精品一区在线| 久久国产精品人妻蜜桃| 国产精品乱码一区二三区的特点| 此物有八面人人有两片| 国产亚洲精品久久久久5区| 一级毛片女人18水好多| 久久精品国产清高在天天线| 国产高清videossex| 欧美一区二区精品小视频在线| 一级片免费观看大全| 久久久久久久久中文| 在线免费观看的www视频| 免费看日本二区| 男人舔女人的私密视频| 国产亚洲av嫩草精品影院| 亚洲国产精品合色在线| 嫩草影院精品99| 91老司机精品| 久久99热这里只有精品18| 国产区一区二久久| 大型黄色视频在线免费观看| 1024手机看黄色片| 久久久久久免费高清国产稀缺| 久久亚洲真实| 国产乱人伦免费视频| 欧美激情高清一区二区三区| 国产熟女xx| 久久香蕉国产精品| 我的亚洲天堂| 日本精品一区二区三区蜜桃| 级片在线观看| 国产精品久久久久久精品电影 | 亚洲av美国av| 观看免费一级毛片| 欧美激情 高清一区二区三区| 欧美成人免费av一区二区三区| 自线自在国产av| 国产亚洲av嫩草精品影院| 又黄又爽又免费观看的视频| 嫩草影院精品99| 免费在线观看日本一区| 国产欧美日韩一区二区三| 精品久久蜜臀av无| 亚洲美女黄片视频| 99久久精品国产亚洲精品| 丝袜在线中文字幕| 久久性视频一级片| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 中文字幕精品免费在线观看视频| 一个人免费在线观看的高清视频| 一级毛片高清免费大全| 露出奶头的视频| 99精品在免费线老司机午夜| 一二三四社区在线视频社区8| 男女那种视频在线观看| 一进一出好大好爽视频| 国产片内射在线| 国产一区在线观看成人免费| 特大巨黑吊av在线直播 | 中文在线观看免费www的网站 | 精品久久久久久久末码| 婷婷亚洲欧美| 一个人观看的视频www高清免费观看 | 国产精品久久视频播放| 少妇 在线观看| 国产成人精品久久二区二区91| 国产1区2区3区精品| АⅤ资源中文在线天堂| 一级毛片精品| 19禁男女啪啪无遮挡网站| 国产aⅴ精品一区二区三区波| 午夜a级毛片| 可以免费在线观看a视频的电影网站| 男人舔奶头视频| 久久香蕉激情| 亚洲成人久久性| 91国产中文字幕| 麻豆成人av在线观看| 国产精品爽爽va在线观看网站 | 国产私拍福利视频在线观看| 伦理电影免费视频| 久久中文看片网| 亚洲天堂国产精品一区在线| 黄片小视频在线播放| 久久热在线av| 亚洲精品美女久久久久99蜜臀| 色综合婷婷激情| 欧美黑人欧美精品刺激| 可以在线观看毛片的网站| 国产午夜精品久久久久久| 精品一区二区三区av网在线观看| 亚洲,欧美精品.| 伦理电影免费视频| 国产私拍福利视频在线观看| 亚洲天堂国产精品一区在线| 99热只有精品国产| 成人欧美大片| 久久久水蜜桃国产精品网| 高潮久久久久久久久久久不卡| 久久中文看片网| 久久久久精品国产欧美久久久| 久久精品夜夜夜夜夜久久蜜豆 | 91大片在线观看| 亚洲成人国产一区在线观看| 黄色a级毛片大全视频| 男女床上黄色一级片免费看| 日韩大尺度精品在线看网址| 啦啦啦 在线观看视频| 91老司机精品| 免费在线观看完整版高清| 午夜久久久久精精品| 亚洲精品久久成人aⅴ小说| 色综合欧美亚洲国产小说| 亚洲国产精品999在线| 亚洲成国产人片在线观看| 久久人人精品亚洲av| 桃色一区二区三区在线观看| 日韩成人在线观看一区二区三区| 老司机午夜福利在线观看视频| 无遮挡黄片免费观看| 国产亚洲欧美98| 天天躁夜夜躁狠狠躁躁| 久久久国产成人免费| 老熟妇仑乱视频hdxx| 99热6这里只有精品| 嫩草影视91久久| 18禁黄网站禁片午夜丰满| 久久久久久九九精品二区国产 | 欧美最黄视频在线播放免费| 亚洲精品美女久久av网站| xxx96com| 长腿黑丝高跟| 99久久99久久久精品蜜桃| 午夜亚洲福利在线播放| 亚洲色图av天堂| 一级片免费观看大全| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 国产精品免费视频内射| 18美女黄网站色大片免费观看| 色婷婷久久久亚洲欧美| 一级片免费观看大全| 变态另类成人亚洲欧美熟女| 51午夜福利影视在线观看| 不卡av一区二区三区| 欧美成人性av电影在线观看| 国产精品1区2区在线观看.| 国产成人影院久久av| 色哟哟哟哟哟哟| 国产真实乱freesex| 韩国精品一区二区三区| 夜夜夜夜夜久久久久| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9| 18禁观看日本| 夜夜夜夜夜久久久久| 97人妻精品一区二区三区麻豆 | 色婷婷久久久亚洲欧美| av视频在线观看入口| 色婷婷久久久亚洲欧美| 男女视频在线观看网站免费 | 国产一卡二卡三卡精品| 国产成人av激情在线播放| 99在线人妻在线中文字幕| 在线国产一区二区在线| 性色av乱码一区二区三区2| 亚洲色图av天堂| 黄片大片在线免费观看| 国产一卡二卡三卡精品| 久久久久精品国产欧美久久久| 国产精品免费一区二区三区在线| 91老司机精品| 视频区欧美日本亚洲| 午夜日韩欧美国产| 18禁美女被吸乳视频| 白带黄色成豆腐渣| 在线观看免费日韩欧美大片| 日韩欧美 国产精品| 午夜激情福利司机影院| 天天躁狠狠躁夜夜躁狠狠躁| 成熟少妇高潮喷水视频| 久久精品国产清高在天天线| 久久 成人 亚洲| 90打野战视频偷拍视频| 又黄又爽又免费观看的视频| 男人的好看免费观看在线视频 | 90打野战视频偷拍视频| 欧美 亚洲 国产 日韩一| 啪啪无遮挡十八禁网站| 精品国产亚洲在线| 午夜精品在线福利| 神马国产精品三级电影在线观看 | 手机成人av网站| 久久天躁狠狠躁夜夜2o2o| 色哟哟哟哟哟哟| 三级毛片av免费| 日韩成人在线观看一区二区三区| 日本在线视频免费播放| 免费高清视频大片| 色综合婷婷激情| 男女之事视频高清在线观看| 熟女电影av网| 一进一出抽搐动态| 中文字幕人成人乱码亚洲影| 中文字幕最新亚洲高清| 亚洲av日韩精品久久久久久密| 一本综合久久免费| 亚洲黑人精品在线| 午夜两性在线视频| 日本一区二区免费在线视频| 精品久久久久久,| 国产成人啪精品午夜网站| 黄片小视频在线播放| 黄色毛片三级朝国网站| 一夜夜www| 亚洲熟女毛片儿| 美女 人体艺术 gogo| avwww免费| 怎么达到女性高潮| 99国产综合亚洲精品| 成人国产综合亚洲| 成年人黄色毛片网站| 日韩欧美国产一区二区入口| 欧美成人免费av一区二区三区| 丰满的人妻完整版| 人妻丰满熟妇av一区二区三区| 国产精品av久久久久免费| 一二三四在线观看免费中文在| 久久久国产成人免费| 亚洲国产毛片av蜜桃av| 国产精品一区二区免费欧美| 叶爱在线成人免费视频播放| 黄网站色视频无遮挡免费观看| 99热6这里只有精品| av欧美777| 亚洲 欧美一区二区三区| 丝袜美腿诱惑在线| 中文亚洲av片在线观看爽| 亚洲成人久久爱视频| 国产又色又爽无遮挡免费看| 久久久久久久久免费视频了| 亚洲性夜色夜夜综合| 国产亚洲精品久久久久久毛片| 中文字幕精品亚洲无线码一区 | 欧美色欧美亚洲另类二区| 男女做爰动态图高潮gif福利片| 成人欧美大片| 色播亚洲综合网| 18禁观看日本| 午夜成年电影在线免费观看| 美女高潮到喷水免费观看| 成年版毛片免费区| 后天国语完整版免费观看| 亚洲成人精品中文字幕电影| 夜夜躁狠狠躁天天躁| 欧美日韩乱码在线| 妹子高潮喷水视频| 精品日产1卡2卡| АⅤ资源中文在线天堂| 欧美激情久久久久久爽电影| 国产极品粉嫩免费观看在线| 久久亚洲精品不卡| 成人精品一区二区免费| 丝袜人妻中文字幕| 最好的美女福利视频网| 法律面前人人平等表现在哪些方面| 亚洲av日韩精品久久久久久密| 老熟妇乱子伦视频在线观看| 欧美日韩福利视频一区二区| 久久精品国产亚洲av香蕉五月| 精品久久久久久久人妻蜜臀av| 久久久久久久精品吃奶| 伊人久久大香线蕉亚洲五| 午夜成年电影在线免费观看| 成在线人永久免费视频| 亚洲成人免费电影在线观看| 亚洲一区高清亚洲精品| АⅤ资源中文在线天堂| 亚洲一区中文字幕在线| 美女 人体艺术 gogo| 国产高清激情床上av| 热re99久久国产66热| 亚洲精品中文字幕一二三四区| 亚洲成人免费电影在线观看| 精品久久蜜臀av无| 亚洲熟妇中文字幕五十中出| 欧美国产精品va在线观看不卡| 国产色视频综合| 国产又色又爽无遮挡免费看| 欧美最黄视频在线播放免费| 国产精品久久久久久人妻精品电影| 老汉色∧v一级毛片| 精品午夜福利视频在线观看一区| 亚洲国产欧洲综合997久久, | 亚洲一卡2卡3卡4卡5卡精品中文| 草草在线视频免费看| 搡老熟女国产l中国老女人| 男人的好看免费观看在线视频 | 久久亚洲真实| 久久久久亚洲av毛片大全| 制服人妻中文乱码| 欧美日韩乱码在线| 亚洲专区中文字幕在线| 精品久久久久久久久久久久久 | 亚洲av片天天在线观看| 高清毛片免费观看视频网站| 欧美av亚洲av综合av国产av| 免费一级毛片在线播放高清视频| 午夜福利在线观看吧| 天天躁狠狠躁夜夜躁狠狠躁| 国产99久久九九免费精品| 母亲3免费完整高清在线观看| 国产三级在线视频| 亚洲欧美日韩高清在线视频| 色播亚洲综合网| 欧美日本亚洲视频在线播放| 精品国产乱子伦一区二区三区| 日韩av在线大香蕉| 国产三级在线视频| 亚洲人成网站高清观看| 欧美成人性av电影在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲欧洲精品一区二区精品久久久| 淫秽高清视频在线观看| 亚洲国产高清在线一区二区三 | 2021天堂中文幕一二区在线观 |