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

    呂宋海峽黑潮西側(cè)3個(gè)氣旋渦的三維結(jié)構(gòu)與演變特征*

    2023-08-03 07:32:26趙偉喬趙秋蘭
    海洋與湖沼 2023年4期
    關(guān)鍵詞:呂宋黑潮中尺度

    趙偉喬 南 峰 趙秋蘭 于 非

    呂宋海峽黑潮西側(cè)3個(gè)氣旋渦的三維結(jié)構(gòu)與演變特征*

    趙偉喬1, 2南 峰2, 3, 4, 5①趙秋蘭1于 非2, 3, 4, 5

    (1. 山東科技大學(xué)數(shù)學(xué)學(xué)院 山東青島 266590; 2. 中國科學(xué)院海洋環(huán)流與波動(dòng)重點(diǎn)實(shí)驗(yàn)室 山東青島 266071; 3. 中國科學(xué)院海洋大科學(xué)研究中心 山東青島 266071; 4. 中國科學(xué)院大學(xué)地球科學(xué)學(xué)院 北京 100049; 5. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室 山東青島 266237)

    黑潮入侵南海會(huì)誘生中尺度渦, 對(duì)南海與西太平洋物質(zhì)能量交換起著重要作用。前人對(duì)黑潮誘生反氣旋渦研究較多, 而對(duì)氣旋渦研究較少, 對(duì)其三維結(jié)構(gòu)和生消過程也不清楚。利用衛(wèi)星高度計(jì)數(shù)據(jù)和再分析數(shù)據(jù), 選取2018~2019年間南海東北部的3個(gè)氣旋渦(cyclonic eddy, CE) CE1、CE2和CE3, 研究了其三維結(jié)構(gòu)與演變特征, 并初步討論了其生成機(jī)制。研究結(jié)果表明: 3個(gè)氣旋渦生成于呂宋海峽西南部黑潮主軸左側(cè), 半徑約為47~87 km, 生成后都向西移動(dòng), 最長距離可達(dá)255 km, 遠(yuǎn)小于該區(qū)域中尺度渦平均移動(dòng)距離。氣旋渦的最大旋轉(zhuǎn)速度約為0.4~0.6 m/s, 垂向深度可達(dá)1 200~1 600 m。3個(gè)氣旋渦中心水體上涌, 溫度異常均為負(fù)值, 在垂向上呈單核結(jié)構(gòu), 冷核的位置在50~600 m處, 冷異常最大可達(dá)-2 °C; 中心鹽度異常垂向上呈現(xiàn)“正-負(fù)-正”的三核結(jié)構(gòu), 分別位于0~100 m、200~400 m和500 m以下深度, 低鹽異常最大可達(dá)-0.26。黑潮鋒面的正壓不穩(wěn)定性是氣旋渦生成的主要因素, 能量從黑潮動(dòng)能向渦動(dòng)能轉(zhuǎn)移, 是氣旋渦生成所需能量的主要來源。

    黑潮; 呂宋海峽; 氣旋渦; 三維結(jié)構(gòu); 能量收支

    黑潮作為最強(qiáng)勁的西邊界流之一, 是連接西太平洋和中國近海物質(zhì)能量交換的重要紐帶, 對(duì)周邊海域氣候變化有著重要調(diào)制作用。呂宋海峽黑潮兩側(cè)都是中尺度渦高發(fā)區(qū)。西太平洋的中尺度渦對(duì)呂宋海峽黑潮流量和流軸都有重要的調(diào)制作用, 可以影響黑潮入侵的路徑(Yang, 2020; Shi, 2021; Yan, 2022), 而黑潮入侵南海又會(huì)誘生新的中尺度渦(Zhang, 2017)。中尺度渦在南海和西太平洋的水交換中起到了重要作用, 其可以輸運(yùn)熱量和鹽量, 并影響浮游生物的分布(Chelton, 2011; Hu, 2014; Yang, 2020; Sun, 2022)。研究黑潮-中尺度渦相互作用對(duì)研究西邊界流變化的動(dòng)力機(jī)制、中尺度渦生消過程、海洋中不同尺度能量串級(jí)和海洋漁業(yè)資源分布等都具有重要意義。

    南海的中尺度渦非?;钴S, 統(tǒng)計(jì)分析表明南海東北部海域是中尺度渦的高發(fā)區(qū)(程旭華等, 2005; 郭景松等, 2007; 林宏陽等, 2012)。該海域同時(shí)存在反氣旋渦和氣旋渦, 反氣旋渦的振幅、半徑和存在時(shí)間通常大于氣旋渦, 但其數(shù)量卻比氣旋渦要少(Nan, 2011c; 王萌等, 2019)。黑潮入侵南海是南海東部海域中尺度渦活躍的一個(gè)重要因素(李燕初等, 2003; 管秉賢等, 2006; 王鼎琦等, 2017)。黑潮在流經(jīng)呂宋海峽時(shí), 常常會(huì)在南海東北部誘發(fā)中尺度渦(陳春濤等, 2010; Jia, 2013)。冬季, 黑潮入侵南海后, 部分黑潮水在臺(tái)灣島西南部流出回到西太平洋, 形成黑潮流套(Looping), 反氣旋渦從流套中甩出, 之后向西移動(dòng), 反氣旋渦的東側(cè)又容易誘生氣旋渦, 形成“渦旋偶極子”現(xiàn)象。夏季, 黑潮入侵南海的強(qiáng)度較弱, 在菲律賓東北部海域一支流會(huì)從南海流出匯入黑潮, 在其左側(cè)容易誘生氣旋渦(Nan, 2011b; Zhang, 2017)。

    很多學(xué)者都對(duì)南海東北部的渦旋進(jìn)行了研究。Chen等(2010)研究了呂宋島東北部一個(gè)反氣旋渦的垂直結(jié)構(gòu), 發(fā)現(xiàn)該渦的溫鹽高于本地水域, 原因是該反氣旋渦攜帶了高溫高鹽的黑潮水。黑潮入侵南海形成的黑潮流套中渦旋容易脫落, 且脫落的渦旋位置大多都在119.5°~ 120°E之間(Jia, 2004)。Zu等(2019)利用渦旋數(shù)據(jù)集和再分析數(shù)據(jù)集研究了南海中尺度渦熱鹽結(jié)構(gòu)的季節(jié)特征和形成機(jī)制, 發(fā)現(xiàn)冬季渦旋溫度異常在50 m以上為偶極子型, 50 m以下為中心對(duì)稱型, 在夏季則整個(gè)都為中心對(duì)稱型, 鹽度異常的特征與之相似, 但不對(duì)稱性較弱。Zhao等(2023)利用潛標(biāo)觀測陣列在呂宋海峽北部捕獲了兩個(gè)攜帶黑潮水向西北方向移動(dòng)的反氣旋渦, 并分析了其生成機(jī)制, 認(rèn)為局地負(fù)風(fēng)應(yīng)力旋度導(dǎo)致了反氣旋渦的產(chǎn)生。為了更好地了解渦旋的結(jié)構(gòu)和演變特征, 很多學(xué)者開始研究渦旋的三維結(jié)構(gòu), 一些學(xué)者利用合成方法構(gòu)造復(fù)合渦的模型來研究其溫鹽特征, 如He等(2018)和謝旭丹等(2018)利用插值的方法合成得到南海不同區(qū)域中尺度渦的三維結(jié)構(gòu), 研究了其溫鹽的三維結(jié)構(gòu)特征及其區(qū)域性差異。另一些學(xué)者通過現(xiàn)場觀測的溫鹽數(shù)據(jù)分析渦旋的三維結(jié)構(gòu)和生消過程(Lin, 2015; Zhang, 2016; Liu, 2019)。Sun等(2021)利用9 a (2000~2008年)的模式數(shù)據(jù)資料, 對(duì)呂宋海峽西北部的氣旋型暖核渦和反氣旋型冷核渦的三維結(jié)構(gòu)進(jìn)行了分析, 認(rèn)為黑潮入侵是南海東北部一些異常渦旋形成的關(guān)鍵機(jī)制。

    綜上, 前人對(duì)南海東北部黑潮流套甩出的反氣旋渦以及臺(tái)灣島西南部氣旋渦研究較多, 但對(duì)于黑潮在呂宋海峽西南部直接誘生的氣旋渦研究較少, 對(duì)這些氣旋渦的垂直結(jié)構(gòu)和演變特征尚不清楚。本文利用衛(wèi)星高度計(jì)數(shù)據(jù), 結(jié)合哥白尼中心再分析數(shù)據(jù), 在呂宋海峽西側(cè)選取了2018~2019年間的3個(gè)氣旋渦(CE1、CE2、CE3), 研究了這3個(gè)氣旋渦溫度、鹽度和流速的三維結(jié)構(gòu)和演變特征, 初步探討了其生成機(jī)制。本研究豐富了對(duì)于黑潮誘生中尺度渦的認(rèn)識(shí), 有助于進(jìn)一步研究黑潮-中尺度渦相互作用過程和機(jī)理。

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

    為研究氣旋渦的三維結(jié)構(gòu)和演變特征, 本文使用的數(shù)據(jù)主要包括衛(wèi)星高度計(jì)數(shù)據(jù)和哥白尼中心再分析數(shù)據(jù)。

    1.1 衛(wèi)星高度計(jì)數(shù)據(jù)

    本文使用了衛(wèi)星高度計(jì)數(shù)據(jù)產(chǎn)品(網(wǎng)址為: https:// resources.marine.copernicus.eu/product-detail/SEALEVEL_GLO_PHY_L4_MY_008_047/), 包含海面絕對(duì)動(dòng)力高度(Absolute Dynamic Topography, ADT)、海面高度異常(sea level anomaly, SLA)、地轉(zhuǎn)流速(,)和地轉(zhuǎn)流速異常(,)等要素, 空間分辨率為0.25°× 0.25°, 時(shí)間分辨率為1 d。利用衛(wèi)星高度計(jì)數(shù)據(jù), 可以統(tǒng)計(jì)分析中尺度渦的時(shí)空分布和運(yùn)動(dòng)特征, 也可以研究中尺度渦對(duì)黑潮的影響(秦麗娟等, 2015; 程建婷等, 2017; 王萌等, 2019)。本文利用衛(wèi)星高度計(jì)資料主要研究氣旋渦運(yùn)動(dòng)學(xué)特征。

    1.2 哥白尼中心再分析數(shù)據(jù)

    為研究氣旋渦的三維結(jié)構(gòu)和演變特征, 本文使用歐洲哥白尼中心再分析數(shù)據(jù)(Global Ocean Multi- Observation Products)(網(wǎng)址為: https://data.marine. copernicus.eu/product/SEALEVEL_GLO_PHY_L4_MY_008_047), 數(shù)據(jù)深度范圍0~5 500 m, 共50層, 包括溫度、鹽度和流速(,)以及海面動(dòng)力高度數(shù)據(jù), 空間分辨率為0.25°×0.25°, 時(shí)間分辨率為7 d。再分析數(shù)據(jù)同化了衛(wèi)星高度計(jì)資料和歷史溫鹽剖面觀測數(shù)據(jù), 其絕對(duì)動(dòng)力高度和流速分布圖與衛(wèi)星高度計(jì)數(shù)據(jù)基本一致(圖1)。再分析數(shù)據(jù)選取的時(shí)間跨度和研究區(qū)域與衛(wèi)星高度計(jì)數(shù)據(jù)相同, 用于研究渦旋的三維結(jié)構(gòu)特征和能量收支分析。

    2 結(jié)果與討論

    2.1 3個(gè)氣旋渦的運(yùn)動(dòng)學(xué)及演變特征

    呂宋海峽兩側(cè)是中尺度渦的高發(fā)區(qū), 利用1993~ 2020年的衛(wèi)星高度計(jì)數(shù)據(jù)和再分析數(shù)據(jù)繪制了呂宋海峽兩側(cè)平均的渦動(dòng)能(eddy kinetic energy, EKE)分布(圖2), 渦動(dòng)能的計(jì)算公式如下:

    其中, EEK表示渦動(dòng)能EKE, u'和v'分別表示流速異常東西分量和南北分量。流速異常是每個(gè)時(shí)刻的流速減去1993~2020年中當(dāng)月流速數(shù)據(jù)的平均值得到的。由圖2可以看出, 渦動(dòng)能高值區(qū)基本上都集中在呂宋海峽兩側(cè), 這一現(xiàn)象和黑潮流經(jīng)呂宋海峽以及入侵南海誘生中尺度渦有關(guān)(陳春濤等, 2010; Jia et al, 2013)。由于再分析數(shù)據(jù)同化了衛(wèi)星高度計(jì)數(shù)據(jù), 因此呂宋海峽兩側(cè)的平均流場特征與渦動(dòng)能的空間分布基本一致(圖1和圖2)。

    注: a: 衛(wèi)星高度計(jì)數(shù)據(jù), b: 再分析數(shù)據(jù)結(jié)果; 紅色圓點(diǎn)為渦旋起始位置, 三條黑線分別代表3個(gè)渦旋(CE1、CE2、CE3)的運(yùn)動(dòng)軌跡; 陸地周邊的空白是數(shù)據(jù)插值后的正?,F(xiàn)象

    圖2 1993~2020年平均渦動(dòng)能(EKE)分布圖

    注: a: 衛(wèi)星高度計(jì)數(shù)據(jù); b: 再分析數(shù)據(jù)

    利用2018~2019年數(shù)據(jù), 在呂宋海峽西側(cè)選取了較為明顯的3個(gè)氣旋渦, 為了分析這3個(gè)渦旋的演化過程, 根據(jù)-0.05 m的海面高度異常等值線(Wang, 2003), 氣旋渦開始有明顯的閉合等高線時(shí)刻定為生成時(shí)刻, 閉合等高線消失時(shí)刻定為消亡時(shí)刻, 中間時(shí)刻即為生成時(shí)刻和消亡時(shí)刻的中間值(圖3), 表1統(tǒng)計(jì)了它們的運(yùn)動(dòng)學(xué)參數(shù)。結(jié)果顯示: 黑潮流經(jīng)呂宋海峽時(shí)以西北向入侵南海, 3個(gè)氣旋渦都是在呂宋海峽黑潮主軸[位置判定參考了丁睿彬等(2013)]西南側(cè)生成, 生成后自東向西移動(dòng), 西傳一段距離后逐漸消亡。氣旋渦的運(yùn)動(dòng)軌跡如圖1所示。CE1在2018年8月23日生成, 2018年9月15日消亡, 生命周期有24 d; CE2在2018年11月21日生成, 2018年12月8日消亡, 生命周期有18 d; CE3在2019年2月13日生成, 在2019年3月25日消亡, 其生命周期達(dá)42 d, 存在時(shí)間最長。CE2存在時(shí)間最短, 其強(qiáng)度也最弱, 其平均半徑只有47 km, 平均海表流速和最大海表流速均小于另外兩個(gè)渦旋, 西傳距離最短, 僅在呂宋海峽西側(cè)向西北方向移動(dòng)了57 km就消亡了。CE1的平均半徑有51 km, 其強(qiáng)度比CE2略強(qiáng), 相對(duì)于CE3較弱, 西傳距離和軌跡與CE3相似, CE1向西移動(dòng)距離有177 km。CE3的平均半徑有87 km, 生命周期有42 d, 其強(qiáng)度最強(qiáng), 海表流速和海面高度異常都是3個(gè)渦旋中最大的, 西傳的距離也最遠(yuǎn), 可達(dá)255 km。這可能和CE3的生命周期在冬末春初有一定關(guān)系, 這個(gè)時(shí)期呂宋海峽風(fēng)向由東向西, 有利于渦旋向西移動(dòng)(杜天時(shí), 2019)??傮w來說3個(gè)氣旋渦存在時(shí)間、半徑和移動(dòng)距離都遠(yuǎn)小于這個(gè)區(qū)域氣旋渦的平均存在時(shí)間、半徑和移動(dòng)距離, 跟這個(gè)區(qū)域反氣旋渦相比差異更大(Nan, 2011c), 說明這些由黑潮誘生的氣旋渦可能具有相對(duì)較短的生命周期。CE1、CE2、CE3的平均海面高度異常和最大海面高度異常依次遞增(表1), CE1和CE2的海面高度異常也要小于這個(gè)區(qū)域氣旋渦的平均海面高度異常, CE3的海面高度異常與這個(gè)區(qū)域氣旋渦平均海面高度異常相差不大(Nan, 2011c)。

    圖3 3個(gè)氣旋渦存在期間的海面動(dòng)力高度和流場

    注: a~c為CE1, d~f為CE2, g~i為CE3; a, d和g為生成時(shí)刻, b, e和h為中間時(shí)刻, c, f和i為消亡時(shí)刻; a, d和g中的黑色加粗曲線為黑潮流軸, 紅色橫線為圖4和圖5所選取的緯向斷面; 數(shù)據(jù)來自再分析數(shù)據(jù)

    表1 3個(gè)渦旋的運(yùn)動(dòng)學(xué)特征統(tǒng)計(jì)表

    Tab.1 Statistics of kinematic properties of the three eddies

    為研究3個(gè)氣旋渦的垂向結(jié)構(gòu)和演變特征, 利用再分析數(shù)據(jù)繪制了橫穿渦旋斷面的流場以及溫度異常(圖4)和鹽度異常(圖5)。溫度異常的計(jì)算方法是將選定區(qū)域某一層深度上的所有溫度值取平均值, 然后所有溫度值都減去這個(gè)平均值得到的即為選定區(qū)域某一層深度上的溫度異常, 鹽度異常同理。由于渦旋從生成到消亡的過程中不斷移動(dòng), 因此根據(jù)圖3中渦旋3個(gè)時(shí)刻所在的位置選取了不同的緯度, 用紅色橫線標(biāo)記。氣旋渦在生成時(shí)刻和消亡時(shí)刻較弱, 主要用于對(duì)氣旋渦演變特征的研究, 關(guān)于影響深度最大值及溫鹽負(fù)異常的研究以中間時(shí)刻為主。結(jié)果顯示: 3個(gè)渦旋最大旋轉(zhuǎn)速度基本都在表層, 旋轉(zhuǎn)速度隨深度遞減?;?.1 m/s流速等值線可以看到, CE1的垂直深度可以延伸到1 400 m, 旋轉(zhuǎn)流速最大可達(dá)0.5 m/s; CE3的垂直深度最大, 可達(dá)1 600 m, 旋轉(zhuǎn)流速最大可達(dá)0.6 m/s; CE2的垂直深度最小, 只到1 200 m, 旋轉(zhuǎn)流速最大為0.4 m/s。3個(gè)氣旋渦中CE3旋轉(zhuǎn)速度最大, 是對(duì)溫鹽特性影響最大的渦旋, CE1次之, CE2最小。

    圖4 氣旋渦溫度異常(°C)和南北向流速(m/s)的斷面分布圖

    注: a~c為CE1, d~f為CE2, g~i為CE3; a, d和g為生成時(shí)刻, b, e和h為中間時(shí)刻, c, f和i為消亡時(shí)刻; 0 m處的黑色加粗橫線和圖3中的紅色加粗橫線表達(dá)含義一致, 黑色等值線為流速; 數(shù)據(jù)來自再分析數(shù)據(jù)

    圖5 氣旋渦鹽度異常和南北向流速(m/s)的斷面分布圖

    3個(gè)氣旋渦溫度異常的中心均為負(fù)異常, 這是由于氣旋渦造成水體輻散涌升所致。氣旋渦的溫度異常在垂向上是單核結(jié)構(gòu), 冷核大約在50~600 m。CE1、CE2和CE3的冷異常在200~300 m左右達(dá)到最大值, 冷異常峰值分別為-1.7、-1.5和-2 °C。渦旋的冷核都是從生成時(shí)刻到中間時(shí)刻逐漸增強(qiáng), 之后又逐漸減弱。雖然3個(gè)渦旋的影響深度差不多, 但CE3的影響范圍明顯更大一些, 冷異常的峰值也大于另外兩個(gè)渦旋。

    氣旋渦中心鹽度異常在垂向上似有三核結(jié)構(gòu)(圖5), 0~100 m為正異常, 200~400 m為負(fù)異常, 500 m以下鹽度又為正異常。這是由于南海背景鹽度場垂向分布為反”S”型(南峰, 2012), 氣旋渦中心海水上升, 造成鹽度異常為“正-負(fù)-正”三核結(jié)構(gòu)(謝旭丹等, 2018)。500 m以下的鹽度異常不大, 因此500 m以下的正鹽度異常區(qū)域相較于另外兩個(gè)核不是很明顯。渦旋引起的鹽度負(fù)異常峰值深度在200~400 m, CE1、CE2和CE3的低鹽度異常峰值分別為-0.23、-0.21和-0.26。選取斷面東側(cè)的溫度和鹽度正異??赡芘c黑潮水入侵南海有關(guān)。

    2.2 3個(gè)氣旋渦的三維結(jié)構(gòu)

    為更直觀展示3個(gè)渦的結(jié)構(gòu), 進(jìn)一步繪制了這3個(gè)渦旋中間時(shí)刻溫度異常和鹽度異常的三維結(jié)構(gòu)(圖6和圖7)。根據(jù)渦旋的斷面垂直結(jié)構(gòu)選取了6個(gè)深度, 分別是0、200、400、600、800和1 000 m。渦旋的溫度異常(圖6)結(jié)果顯示: 渦旋的冷核在200和400 m處比較明顯, CE3低溫異常的影響范圍最大, 其次是CE1, CE2的影響范圍最小。渦旋溫度負(fù)異常的最大值在200 m左右, 之后負(fù)異常隨著深度的增加逐漸減小。渦旋的鹽度三維異常(圖7)結(jié)果顯示: 低鹽核也是在200~400 m比較明顯, 這一點(diǎn)和冷核一樣, 從600 m開始, 鹽度異常開始出現(xiàn)一些微弱的正異常。200~400 m處, 渦旋鹽度異常結(jié)構(gòu)的中心基本為負(fù)異常。400 m以下, 隨著深度的增加, 渦旋鹽度異常結(jié)構(gòu)的中心反而呈現(xiàn)正異常, 這是氣旋渦引起該區(qū)域海水上升的結(jié)果。從垂直結(jié)構(gòu)上看, 表層流速最大可達(dá)0.6 m/s, 流速隨著深度的增加而急劇減小, 3個(gè)渦旋在0~1 000 m都能清楚地觀察到渦旋結(jié)構(gòu)。根據(jù)不同深度溫鹽異常剖面的渦旋結(jié)構(gòu), 可以發(fā)現(xiàn)3個(gè)氣旋渦的另一個(gè)特征是隨著深度的增加, 渦旋的中心位置沒有發(fā)生明顯的偏移, 這與Zhang等(2016)觀測到的南海北部中尺度渦中心西向傾斜結(jié)構(gòu)不同。

    圖6 氣旋渦中間時(shí)刻溫度異常和流場的三維結(jié)構(gòu)

    圖7 氣旋渦中間時(shí)刻鹽度異常和流場的三維結(jié)構(gòu)

    結(jié)合3個(gè)渦旋溫鹽異常的三維結(jié)構(gòu), 可以發(fā)現(xiàn)CE2和CE3的東北側(cè)有溫鹽異常明顯升高的現(xiàn)象, 這一現(xiàn)象可能與CE2和CE3的中間時(shí)刻所處的時(shí)期是黑潮入侵南海強(qiáng)度最強(qiáng)的季節(jié)有關(guān)(方國洪等, 2002; 王兆毅等, 2016), 高溫高鹽的黑潮水入侵南海, 可以改變流經(jīng)區(qū)域的背景溫鹽特性。CE1的中間時(shí)刻在夏季末, 這個(gè)時(shí)期是黑潮入侵南海強(qiáng)度比較弱的季節(jié), 渦旋東北側(cè)溫鹽異常升高的現(xiàn)象相對(duì)于另外兩個(gè)渦旋較弱可能也與此有關(guān)。以上分析表明, 不同氣旋渦的溫鹽異常主體結(jié)構(gòu)及影響深度雖然相似, 但還是有一些差異, 可能跟渦旋存在期間風(fēng)場和黑潮入侵等因素有關(guān)。

    2.3 渦旋生成時(shí)刻的能量收支分析

    海洋中的能量形式主要有渦動(dòng)能(EK)、平均流動(dòng)能(MK)、渦動(dòng)勢能(EP)和平均流勢能(MP)四種(Xie, 2007), 計(jì)算公式分別為

    其中,是平均位勢密度的垂向梯度,'是密度異常,和是流速異常,為重力加速度, 上劃線表示時(shí)間平均, 點(diǎn)撇表示和平均值的偏差。

    平均動(dòng)能與渦動(dòng)能的相互作用與正壓不穩(wěn)定有關(guān), 正壓不穩(wěn)定過程由水平剪切主導(dǎo)。正壓轉(zhuǎn)換率(BT)為正值時(shí)代表能量從平均流動(dòng)能向渦動(dòng)能轉(zhuǎn)換, 為負(fù)值時(shí)代表渦旋能量衰減。平均勢能與渦勢能的相互作用與斜壓不穩(wěn)定有關(guān), 斜壓不穩(wěn)定過程由水平溫度/密度梯度或者速度垂向剪切主導(dǎo)。斜壓轉(zhuǎn)換率(BC)為正值時(shí)代表能量從平均勢能向渦勢能轉(zhuǎn)換, 最后再轉(zhuǎn)移到渦動(dòng)能,BC為負(fù)值則代表渦旋能量衰減, 能量從渦向平均流動(dòng)能和勢能轉(zhuǎn)化(徐安琪, 2021; Yan, 2022)。正壓轉(zhuǎn)換率BT和斜壓轉(zhuǎn)換率BC的計(jì)算公式如下:

    其中,和分別是緯向和經(jīng)向,是垂直流速異常。再分析數(shù)據(jù)中只有水平流速, 垂向流速可由水平流速和深度等數(shù)據(jù)計(jì)算得到(Greene, 2019), 計(jì)算公式如下:

    方程兩邊對(duì)深度積分后可得:

    基于圖4和圖5, 可以發(fā)現(xiàn)在氣旋渦的生成時(shí)刻, 渦旋的影響深度大致在600 m, 因此在研究渦旋生成機(jī)制時(shí)只選取0~600 m的深度范圍。圖8是3個(gè)氣旋渦生成時(shí)刻0~600 m的深度平均的正壓轉(zhuǎn)換率和斜壓轉(zhuǎn)化率水平分布圖, 從圖8可以看出: 南海東北部正壓轉(zhuǎn)換率和斜壓轉(zhuǎn)換率正負(fù)值交替分布, 說明這個(gè)區(qū)域平均流和渦旋之間動(dòng)能的能量轉(zhuǎn)換非常頻繁。在渦旋生成時(shí)刻, 渦旋位置附近正壓轉(zhuǎn)換率基本都為正值, 表示能量從平均流動(dòng)能向渦動(dòng)能轉(zhuǎn)移; 渦旋位置附近的斜壓轉(zhuǎn)換率基本為負(fù)值, 說明一部分能量從渦勢能轉(zhuǎn)化成平均流勢能。因此, 黑潮鋒面的水平剪切導(dǎo)致的正壓不穩(wěn)定性是這3個(gè)氣旋渦生成所需能量的主要來源。Zhang等(2017)研究發(fā)現(xiàn)黑潮流套北側(cè)的水平剪切導(dǎo)致的正壓不穩(wěn)定性是臺(tái)灣島西南海域氣旋渦的生成所需能量的主要來源。本文研究的3個(gè)氣旋渦在黑潮鋒面左側(cè)生成, 黑潮鋒面不穩(wěn)定性是呂宋海峽西側(cè)氣旋渦生成主要原因, 氣旋渦生成后能量從渦勢能向平均流勢能轉(zhuǎn)化, 可能是這個(gè)區(qū)域氣旋渦存在時(shí)間較短的原因。

    3 結(jié)論

    本文利用衛(wèi)星高度計(jì)數(shù)據(jù)和再分析數(shù)據(jù), 在呂宋海峽西側(cè)選取了2018~2019年間的3個(gè)氣旋渦, 分析統(tǒng)計(jì)了氣旋渦的三維結(jié)構(gòu)和演變特征, 并初步討論了其生成機(jī)制。研究結(jié)果如下:

    圖8 3個(gè)氣旋渦生成時(shí)刻0~600 m深度平均的正壓轉(zhuǎn)換率RBT和斜壓轉(zhuǎn)化率RBC水平分布圖

    注: 黑色圓點(diǎn)為渦生成時(shí)刻渦核所在位置, a和b為CE1, c和d為CE2, e和f為CE3; a, c和e為正壓轉(zhuǎn)換率BT (單位: W/m3), b, d和f為斜壓轉(zhuǎn)化率BC (單位: W/m3)

    (1) 氣旋渦在呂宋海峽西南部黑潮主軸左側(cè)生成, 生成后都向西移動(dòng), 3個(gè)渦存在時(shí)間在半個(gè)月到一個(gè)半月之間, 半徑在47~87 km之間, 西向移動(dòng)距離在57~255 km之間。CE3是3個(gè)渦旋中強(qiáng)度最強(qiáng)的渦, 其生命周期、半徑和西向移動(dòng)距離都要大于另外兩個(gè)渦旋, CE1次之, CE2最弱。3個(gè)氣旋渦具有相對(duì)較短的生命周期, 其存在時(shí)間、半徑和西傳距離都要小于該海域氣旋渦的平均值。

    (2) 氣旋渦垂直深度可達(dá)1 200~1 600 m, CE2、CE1和CE3的垂直深度依次遞增。中心溫度為負(fù)異常, 垂向上呈現(xiàn)單核結(jié)構(gòu), 冷核的位置大約在50~ 600 m處, 200~300 m左右達(dá)到最大值, 峰值在-1.5~ 2 °C之間。氣旋渦鹽度異常的中心也為負(fù)異常, 中心鹽度異常垂向上呈現(xiàn)“正-負(fù)-正”的三核結(jié)構(gòu), 低鹽異常的峰值可達(dá)-0.26。

    (3) 氣旋渦最大旋轉(zhuǎn)速度在表層可以達(dá)到0.4~ 0.6 m/s, 隨著深度的增加, 流速逐漸減小, 渦旋中心隨深度沒有出現(xiàn)明顯傾斜。渦旋東北側(cè)海水溫鹽異常與不同渦旋存在期間黑潮入侵南海的強(qiáng)弱有關(guān)。

    (4) 呂宋海峽黑潮的水平剪切導(dǎo)致的正壓不穩(wěn)定性是氣旋渦生成的主要因素, 能量從平均流動(dòng)能向渦動(dòng)能轉(zhuǎn)移, 為氣旋渦的生成提供了能量。

    本文以2018~2019年黑潮誘生的3個(gè)氣旋渦為例研究了其三維結(jié)構(gòu)特征和生成機(jī)制, 但氣旋渦每年都會(huì)產(chǎn)生, 我們沒有統(tǒng)計(jì)所有年份的氣旋渦, 其季節(jié)和年際變化規(guī)律尚不清楚?因此針對(duì)該區(qū)域更多氣旋渦旋開展研究是我們下一步的工作。

    丁睿彬, 陳大可, 姜良紅, 2013. 黑潮流軸在呂宋海峽的變化分析[J]. 海洋學(xué)研究, 31(2): 16-25.

    王兆毅, 劉桂梅, 王輝, 等, 2016. 呂宋海峽水交換季節(jié)和年際變化特征的數(shù)值模擬研究[J]. 海洋學(xué)報(bào), 38(5): 1-13.

    王萌, 張艷偉, 劉志飛, 等, 2019. 南海北部中尺度渦的時(shí)空分布特征: 基于衛(wèi)星高度計(jì)資料的統(tǒng)計(jì)分析[J]. 地球科學(xué)進(jìn)展, 34(10): 1069-1080.

    王鼎琦, 方國洪, 邱婷, 2017. 呂宋海峽黑潮脫落渦旋的特征分析[J]. 海洋與湖沼, 48(4): 672-681.

    方國洪, 魏澤勛, 王凱, 等, 2002. 中國近海域際水、熱、鹽輸運(yùn): 全球變網(wǎng)格模式結(jié)果[J]. 中國科學(xué)(D輯), 32(12): 969-977.

    杜天時(shí), 2019. 冬季南海表層黑潮的入侵形式以及相關(guān)風(fēng)場的影響[J]. 海洋預(yù)報(bào), 36(4): 60-67.

    李燕初, 蔡文理, 李立, 等, 2003. 南海東北部海域中尺度渦的季節(jié)和年際變化[J]. 熱帶海洋學(xué)報(bào), 22(3): 61-70.

    陳春濤, 劉佳佳, 2010. 東海北部黑潮鋒面渦的快速生成和消散[J]. 遙感信息, 25(3): 103-107, 128.

    林宏陽, 胡建宇, 鄭全安, 2012. 呂宋海峽附近中尺度渦特征的統(tǒng)計(jì)分析[J]. 海洋學(xué)報(bào), 34(1): 1-7.

    南峰, 2012. 臺(tái)灣西南部海域流—渦結(jié)構(gòu)及其演變規(guī)律研究[D]. 青島: 中國海洋大學(xué): 59-61.

    秦麗娟, 董慶, 樊星, 等, 2015. 衛(wèi)星高度計(jì)的北太平洋中尺度渦時(shí)空分析[J]. 遙感學(xué)報(bào), 19(5): 806-817.

    徐安琪, 2021. 西北太平洋次表層中尺度渦特征及其動(dòng)力機(jī)制研究[D]. 青島: 中國科學(xué)院大學(xué)(中國科學(xué)院海洋研究所): 70-73.

    郭景松, 袁業(yè)立, 熊學(xué)軍, 等, 2007. 呂宋海峽兩側(cè)中尺度渦統(tǒng)計(jì)[J]. 海洋科學(xué)進(jìn)展, 25(2): 139-148.

    程旭華, 齊義泉, 王衛(wèi)強(qiáng), 2005. 南海中尺度渦的季節(jié)和年際變化特征分析[J]. 熱帶海洋學(xué)報(bào), 24(4): 51-59.

    程建婷, 楊德周, 尹寶樹, 等, 2017. 臺(tái)灣以東中尺度渦對(duì)黑潮入侵東海路徑的影響[J]. 海洋科學(xué), 41(2): 81-88.

    謝旭丹, 王靜, 儲(chǔ)小青, 等, 2018. 南海中尺度渦溫鹽異常三維結(jié)構(gòu)[J]. 海洋學(xué)報(bào), 40(4): 1-14.

    管秉賢, 袁耀初, 2006. 中國近海及其附近海域若干渦旋研究綜述Ⅰ. 南海和臺(tái)灣以東海域[J]. 海洋學(xué)報(bào), 28(3): 1-16.

    CHELTON D B, GAUBE P, SCHLAX M G,, 2011. The influence of nonlinear mesoscale eddies on near-surface oceanic chlorophyll [J]. Science, 334(6054): 328-332, doi: 10.1126/science.1208897.

    CHEN G X, HOU Y J, CHU X Q,, 2010. Vertical structure and evolution of the Luzon Warm Eddy [J]. Chinese Journal of Oceanology and Limnology, 28(5): 955-961.

    GREENE C A, THIRUMALAI K, KEARNEY K A,, 2019. The climate data toolbox for MATLAB [J]. Geochemistry, Geophysics, Geosystems, 20(7): 3774-3781, doi: 10.1029/ 2019GC008392.

    HE Q Y, ZHAN H G, CAI S Q,, 2018. A new assessment of mesoscale eddies in the South China Sea: surface features, three-dimensional structures, and thermohaline transports [J]. Journal of Geophysical Research: Oceans, 123(7): 4906-4929.

    HU Z F, TAN Y H, SONG X Y,, 2014. Influence of mesoscale eddies on primary production in the South China Sea during spring inter-monsoon period [J]. Acta Oceanologica Sinica, 33(3): 118-128, doi: 10.1007/s13131-014-0431-8.

    JIA Y L, LIU Q Y, 2004. Eddy shedding from the Kuroshio bend at Luzon strait [J]. Journal of Oceanography, 60(6): 1063-1069.

    JIA Y J, ZHANG Y G, LIN M S, 2013. The numerical simulation of the Kuroshio frontal eddies in the East China Sea using a hybrid coordinate ocean mode [J]. Acta Oceanologica Sinica, 32(5): 31-41.

    LIN X Y, DONG C M, CHEN D K,, 2015. Three-dimensional properties of mesoscale eddies in the South China Sea based on eddy-resolving model output [J]. Deep Sea Research Part I: Oceanographic Research Papers, 99: 46-64.

    LIU Z H, CHEN X R, YU J C,, 2019. Kuroshio intrusion into the South China Sea with an anticyclonic eddy: evidence from underwater glider observation [J]. Journal of Oceanology and Limnology, 37(5): 1469-1480.

    NAN F, XUE H J, CHAI F,, 2011b. Identification of different types of Kuroshio intrusion into the South China Sea [J]. Ocean Dynamics, 61(9): 1291-1304.

    NAN F, XUE H J, XIU P,, 2011c. Oceanic eddy formation and propagation southwest of Taiwan [J]. Journal of Geophysical Research: Oceans, 116(C12): C12045, doi: 10. 1029/2011JC007386.

    SHI Q, WANG G H, 2021. Meander response of the Kuroshio in the East China Sea to impinging eddies [J]. Journal of Geophysical Research: Oceans, 126(9): e2021JC017512.

    SUN W J, LIU Y, CHEN G X,, 2021. Three-dimensional properties of mesoscale cyclonic warm-core and anticyclonic cold-core eddies in the South China Sea [J]. Acta Oceanologica Sinica, 40(10): 17-29.

    SUN Z B, ZHANG Z W, QIU B,, 2022. Subsurface mesoscale eddies observed in the northeastern South China Sea: dynamic features and water mass transport [J]. Journal of Physical Oceanography, 52(5): 841-855, doi: 10.1175/ JPO-D-21-0177.1.

    WANG G H, SU J L, CHU P C, 2003. Mesoscale eddies in the South China Sea observed with altimeter data [J]. Geophysical Research Letters, 30(21): 2121.

    XIE L, LIU X M, PIETRAFESA L J, 2007. Effect of bathymetric curvature on gulf stream instability in the vicinity of the Charleston bump [J]. Journal of Physical Oceanography, 37(3): 452-475.

    YAN X M, KANG D J, PANG C G,, 2022. Energetics analysis of the eddy-Kuroshio interaction East of Taiwan [J]. Journal of Physical Oceanography, 52(4): 647-664.

    YANG Q, LIU H L, LIN P F, 2020. The effect of oceanic mesoscale eddies on the looping path of the Kuroshio intrusion in the Luzon Strait [J]. Scientific Reports, 10(1): 636.

    ZHANG Z W, TIAN J W, QIU B,, 2016. Observed 3D structure, generation, and dissipation of oceanic mesoscale eddies in the South China Sea [J]. Scientific Reports, 6: 24349.

    ZHANG Z W, ZHAO W, QIU B,, 2017. Anticyclonic eddy sheddings from Kuroshio loop and the accompanying cyclonic eddy in the northeastern South China Sea [J]. Journal of Physical Oceanography, 47(6): 1243-1259.

    ZHAO R X, ZHU X H, ZHANG C Z,, 2023. Summer anticyclonic eddies carrying Kuroshio waters observed by a large CPIES array west of the Luzon Strait [J]. Journal of Physical Oceanography, 53(1): 341-359, doi: 10.1175/JPO- D-22-0019.1.

    ZU Y C, SUN S W, ZHAO W,, 2019. Seasonal characteristics and formation mechanism of the thermohaline structure of mesoscale eddy in the South China Sea [J]. Acta Oceanologica Sinica, 38(4): 29-38, doi: 10.1007/s13131-018-1222-4.

    THREE-DIMENSIONAL STRUCTURE AND EVOLUTION OF THREE CYCLONIC EDDIES TO THE WEST OF THE KUROSHIO IN THE LUZON STRAIT

    ZHAO Wei-Qiao1, 2, NAN Feng2, 3, 4, 5, ZHAO Qiu-Lan1, YU Fei2, 3, 4, 5

    (1. School of Mathematics, Shandong University of Science and Technology, Qingdao 266590, China; 2. Key Laboratory of Ocean Circulation and Waves, Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China; 3. Center for Ocean Mega-science, Chinese Academy of Sciences, Qingdao 266071, China; 4. College of Earth and Planetary Sciences, University of Chinese Academy of Sciences, Beijing 100049, China; 5. Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China)

    The Kuroshio intrusion into the South China Sea (SCS) can induce mesoscale eddies, which play an important role in the exchange of matter and energy between the SCS and Western Pacific Ocean. Previous studies focus on anticyclonic eddies induced by Kuroshio and eddy shedding processes. However, cyclonic eddies are less investigated, their three-dimensional structure and evolution remain unclear. Three cyclonic eddies (CE1, CE2, and CE3) in the northeastern SCS formed during 2018~2019 were selected. By using satellite altimeter data and reanalysis data, three-dimensional structure and evolution characteristics of the eddies were studied, and their forming mechanism was discussed. Results show that: cyclonic eddies were generated in the southwest of the Kuroshio in the Luzon Strait, their mean radius were between 47~87 km. They all moved westward in the maximum distance of 255 km, which was far less than the average moving distance of regional mesoscale eddies. The cyclonic eddies rotated counterclockwise in the maximum velocity of 0.4~0.6 m/s, and the vertical depth extended to 1 200~1 600 m. The water in the centers of the cyclonic eddies were uplifted, and the temperature in eddy center showed a negative anomaly with a single core in vertical direction. The cold core was located at 50~600 m, and the maximum cold anomaly reached-2 °C. The central salinity anomaly presented a vertical three-core structure of “positive-negative-positive” located at 0~100 m, 200~400 m, and below 500 m in depth, respectively, and the lowest salinity anomaly reached-0.26. The barotropic instability of the Kuroshio front was the controlling factor on the generation of the cyclonic eddies.

    Kuroshio; Luzon Strait; cyclonic eddies; three-dimensional structure; energy budget

    * 國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目, 2022YFB3205300號(hào); 國家自然科學(xué)基金面上項(xiàng)目, 41676005號(hào)。趙偉喬, 碩士研究生, E-mail: qiao20221109@163.com

    南 峰, 博士生導(dǎo)師, E-mail: nanfeng0515@qdio.ac.cn

    2022-11-11,

    2023-01-17

    P731

    10.11693/hyhz20221100295

    猜你喜歡
    呂宋黑潮中尺度
    譜松弛動(dòng)力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動(dòng)研究
    郁南縣河口鎮(zhèn)芒果喜獲豐收
    源流(2020年6期)2020-08-03 01:59:35
    基于深度學(xué)習(xí)的中尺度渦檢測技術(shù)及其在聲場中的應(yīng)用
    以新聞史視角探析晚清《申報(bào)》中呂宋票廣告鏡像
    傳播力研究(2019年1期)2019-03-20 11:28:14
    2016年7月四川持續(xù)性強(qiáng)降水的中尺度濾波分析
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計(jì)數(shù)據(jù)的黑潮變異特征
    東海冷渦對(duì)黑潮熱輸運(yùn)年代際躍變的響應(yīng)
    午夜免费激情av| 国产午夜福利久久久久久| 婷婷精品国产亚洲av| 日本-黄色视频高清免费观看| 老熟妇仑乱视频hdxx| 欧美丝袜亚洲另类 | 亚洲欧美清纯卡通| 亚洲成人久久爱视频| 一区二区三区免费毛片| 久久中文看片网| 精品国产三级普通话版| 91麻豆精品激情在线观看国产| 美女黄网站色视频| 男女边吃奶边做爰视频| 亚洲综合色惰| 麻豆久久精品国产亚洲av| 黄色日韩在线| 国产午夜精品论理片| 亚洲三级黄色毛片| 少妇人妻精品综合一区二区 | 亚洲狠狠婷婷综合久久图片| 国产一区二区激情短视频| 99久久中文字幕三级久久日本| 亚洲av美国av| 日韩欧美三级三区| 一区福利在线观看| 他把我摸到了高潮在线观看| or卡值多少钱| 最近视频中文字幕2019在线8| 欧美黑人欧美精品刺激| 国产精品三级大全| 欧美激情在线99| 成年人黄色毛片网站| 成人鲁丝片一二三区免费| 免费人成在线观看视频色| 欧美最新免费一区二区三区| 大型黄色视频在线免费观看| 国产午夜精品久久久久久一区二区三区 | 午夜福利成人在线免费观看| 十八禁国产超污无遮挡网站| 亚洲av不卡在线观看| 黄色丝袜av网址大全| 欧美zozozo另类| 久久久色成人| 久久久久国内视频| 国产v大片淫在线免费观看| 免费在线观看日本一区| 色综合色国产| 99久久九九国产精品国产免费| 一个人看视频在线观看www免费| 老熟妇乱子伦视频在线观看| 欧美日韩国产亚洲二区| 小说图片视频综合网站| 国产午夜精品久久久久久一区二区三区 | 一卡2卡三卡四卡精品乱码亚洲| 校园春色视频在线观看| 日日夜夜操网爽| 一个人看的www免费观看视频| 国产熟女欧美一区二区| 一级黄色大片毛片| 国内少妇人妻偷人精品xxx网站| 欧美日韩黄片免| 精品久久久久久,| 久久久久久伊人网av| 九九久久精品国产亚洲av麻豆| 亚洲熟妇熟女久久| 黄色视频,在线免费观看| 深爱激情五月婷婷| 亚洲欧美日韩卡通动漫| 欧美性感艳星| 日韩欧美免费精品| 麻豆国产97在线/欧美| 黄色视频,在线免费观看| 最后的刺客免费高清国语| 性插视频无遮挡在线免费观看| 午夜福利在线观看吧| 一本久久中文字幕| 久9热在线精品视频| 久久精品国产亚洲av天美| 又爽又黄a免费视频| 国产激情偷乱视频一区二区| 成人美女网站在线观看视频| 尤物成人国产欧美一区二区三区| 91在线观看av| .国产精品久久| 22中文网久久字幕| 男人舔奶头视频| 女的被弄到高潮叫床怎么办 | 日本成人三级电影网站| 成人国产综合亚洲| 最近在线观看免费完整版| 狂野欧美激情性xxxx在线观看| 九色成人免费人妻av| 欧美黑人欧美精品刺激| 少妇裸体淫交视频免费看高清| 天天躁日日操中文字幕| 给我免费播放毛片高清在线观看| 成年免费大片在线观看| 国产av一区在线观看免费| 夜夜夜夜夜久久久久| 国产亚洲av嫩草精品影院| 欧美成人一区二区免费高清观看| 给我免费播放毛片高清在线观看| 麻豆av噜噜一区二区三区| 亚洲精品久久国产高清桃花| 欧美性感艳星| 免费无遮挡裸体视频| 免费在线观看日本一区| 九九爱精品视频在线观看| 国产探花极品一区二区| 蜜桃久久精品国产亚洲av| 久久精品国产99精品国产亚洲性色| 综合色av麻豆| 国产精品精品国产色婷婷| 美女高潮喷水抽搐中文字幕| 噜噜噜噜噜久久久久久91| 成人三级黄色视频| 亚洲精品粉嫩美女一区| 国产精品一及| 亚洲国产精品成人综合色| 两个人视频免费观看高清| 搡老妇女老女人老熟妇| www.色视频.com| 婷婷精品国产亚洲av在线| 亚洲 国产 在线| 中文字幕av成人在线电影| 亚洲美女黄片视频| 又粗又爽又猛毛片免费看| 亚洲 国产 在线| 精品一区二区免费观看| 久久99热这里只有精品18| 岛国在线免费视频观看| 岛国在线免费视频观看| 亚洲精品影视一区二区三区av| 国产精品三级大全| 久久国内精品自在自线图片| 一进一出抽搐动态| 乱系列少妇在线播放| 欧美在线一区亚洲| 精品久久久久久久末码| 一夜夜www| 两个人视频免费观看高清| 超碰av人人做人人爽久久| 丰满乱子伦码专区| 欧美性感艳星| 两人在一起打扑克的视频| 内地一区二区视频在线| 日韩一区二区视频免费看| 成人欧美大片| 久久久久久久久久久丰满 | 国产一区二区三区av在线 | 日本一本二区三区精品| 在现免费观看毛片| 国产伦人伦偷精品视频| 综合色av麻豆| 久久精品影院6| 亚洲最大成人av| 深夜精品福利| 国产精品一区二区三区四区久久| 国产一区二区在线观看日韩| 人人妻,人人澡人人爽秒播| 在线观看免费视频日本深夜| 亚洲va日本ⅴa欧美va伊人久久| 窝窝影院91人妻| 国产伦一二天堂av在线观看| 久久久午夜欧美精品| 亚洲av第一区精品v没综合| 一级黄片播放器| 亚洲av中文字字幕乱码综合| 国产精品综合久久久久久久免费| 国产黄色小视频在线观看| 欧美高清性xxxxhd video| 久久精品91蜜桃| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 看片在线看免费视频| 国产精品一区二区性色av| 美女高潮喷水抽搐中文字幕| 久久久国产成人免费| 中文资源天堂在线| 美女免费视频网站| 国产精品美女特级片免费视频播放器| 九色国产91popny在线| 大型黄色视频在线免费观看| 99九九线精品视频在线观看视频| 国产av一区在线观看免费| 国语自产精品视频在线第100页| 国产欧美日韩精品亚洲av| 最后的刺客免费高清国语| 小蜜桃在线观看免费完整版高清| 大型黄色视频在线免费观看| 日本 欧美在线| 久久久久久伊人网av| 亚洲av熟女| 女的被弄到高潮叫床怎么办 | 亚洲人成网站高清观看| 2021天堂中文幕一二区在线观| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 免费观看在线日韩| 国产成年人精品一区二区| 日本爱情动作片www.在线观看 | 老师上课跳d突然被开到最大视频| 国产免费一级a男人的天堂| 国产蜜桃级精品一区二区三区| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 国产精品,欧美在线| 亚洲一区二区三区色噜噜| av在线天堂中文字幕| 乱码一卡2卡4卡精品| av福利片在线观看| 少妇裸体淫交视频免费看高清| 日本黄色片子视频| av.在线天堂| 我的女老师完整版在线观看| 99热6这里只有精品| 中国美白少妇内射xxxbb| 丰满的人妻完整版| 久久久色成人| 亚洲经典国产精华液单| 亚洲av美国av| 国产精品永久免费网站| 婷婷精品国产亚洲av| www.色视频.com| 午夜老司机福利剧场| 别揉我奶头 嗯啊视频| 免费搜索国产男女视频| 久久精品综合一区二区三区| 中亚洲国语对白在线视频| 亚洲成人久久爱视频| 美女高潮的动态| 欧美区成人在线视频| 国产精品av视频在线免费观看| 最近在线观看免费完整版| 亚洲人成网站高清观看| 搡女人真爽免费视频火全软件 | av天堂中文字幕网| 国产成人福利小说| 国产视频内射| 国产精品亚洲一级av第二区| 99在线人妻在线中文字幕| 成熟少妇高潮喷水视频| 亚洲电影在线观看av| 超碰av人人做人人爽久久| 91麻豆av在线| 亚洲 国产 在线| 黄色视频,在线免费观看| 国产综合懂色| 此物有八面人人有两片| 久久午夜福利片| 国产又黄又爽又无遮挡在线| 免费一级毛片在线播放高清视频| 日韩中文字幕欧美一区二区| av在线观看视频网站免费| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品伦人一区二区| 国产在视频线在精品| 久久久久精品国产欧美久久久| av.在线天堂| 黄片wwwwww| 1000部很黄的大片| 精品久久国产蜜桃| 精品免费久久久久久久清纯| 日韩亚洲欧美综合| 国模一区二区三区四区视频| 精品久久久久久久末码| 欧美另类亚洲清纯唯美| 国产麻豆成人av免费视频| 三级毛片av免费| 免费观看精品视频网站| 亚洲中文日韩欧美视频| 中文资源天堂在线| 97碰自拍视频| 久久久久久久久久成人| 内射极品少妇av片p| 变态另类成人亚洲欧美熟女| 精品久久久久久久久亚洲 | 亚洲性夜色夜夜综合| 一本精品99久久精品77| 成人国产一区最新在线观看| 久久久久九九精品影院| 午夜精品一区二区三区免费看| 欧美+亚洲+日韩+国产| 最近视频中文字幕2019在线8| 亚洲天堂国产精品一区在线| 成人国产一区最新在线观看| 好男人在线观看高清免费视频| 国产伦精品一区二区三区四那| 少妇被粗大猛烈的视频| 成人无遮挡网站| 国产亚洲精品综合一区在线观看| 丰满的人妻完整版| 老司机午夜福利在线观看视频| 久久香蕉精品热| 国产伦人伦偷精品视频| 亚洲天堂国产精品一区在线| 亚洲美女视频黄频| 日本欧美国产在线视频| 亚洲专区国产一区二区| 久久久久久国产a免费观看| 午夜久久久久精精品| 18+在线观看网站| 男人舔女人下体高潮全视频| 亚洲成人久久爱视频| 国产熟女欧美一区二区| 欧美一区二区亚洲| 人人妻人人看人人澡| 美女免费视频网站| 91精品国产九色| 中文亚洲av片在线观看爽| 欧美日韩乱码在线| 高清在线国产一区| 日韩中字成人| 久久午夜亚洲精品久久| 精品不卡国产一区二区三区| 国产一区二区三区在线臀色熟女| 可以在线观看毛片的网站| 又粗又爽又猛毛片免费看| 日本一二三区视频观看| 免费av观看视频| 久久中文看片网| 亚洲成人久久爱视频| 免费av不卡在线播放| 久久精品影院6| 老师上课跳d突然被开到最大视频| 日韩一本色道免费dvd| 久久久久九九精品影院| 欧美激情国产日韩精品一区| 综合色av麻豆| 最近最新中文字幕大全电影3| 午夜a级毛片| 午夜精品久久久久久毛片777| 在线a可以看的网站| 国产一区二区亚洲精品在线观看| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 欧美xxxx性猛交bbbb| 欧美激情国产日韩精品一区| 偷拍熟女少妇极品色| 国产激情偷乱视频一区二区| 午夜福利欧美成人| 日本一二三区视频观看| 欧美高清成人免费视频www| 精品午夜福利在线看| 国产精华一区二区三区| 日本一本二区三区精品| 黄片wwwwww| 国产单亲对白刺激| 在线a可以看的网站| 亚洲精华国产精华精| 波多野结衣巨乳人妻| 99在线人妻在线中文字幕| 亚洲av不卡在线观看| 亚洲成人精品中文字幕电影| 午夜精品在线福利| 色哟哟哟哟哟哟| 亚洲精品国产成人久久av| 窝窝影院91人妻| 午夜免费成人在线视频| 中文字幕精品亚洲无线码一区| 尾随美女入室| 极品教师在线免费播放| 尤物成人国产欧美一区二区三区| 久久婷婷人人爽人人干人人爱| 最新在线观看一区二区三区| 久久欧美精品欧美久久欧美| 亚洲av熟女| 免费在线观看日本一区| 国产黄a三级三级三级人| 成年女人毛片免费观看观看9| 欧美又色又爽又黄视频| 亚洲精品久久国产高清桃花| 成年版毛片免费区| or卡值多少钱| 久久精品国产亚洲av涩爱 | 国产v大片淫在线免费观看| 国产精品福利在线免费观看| 久久久久九九精品影院| 中文亚洲av片在线观看爽| 精品久久国产蜜桃| 午夜免费男女啪啪视频观看 | 亚洲av.av天堂| 欧美极品一区二区三区四区| 欧美+亚洲+日韩+国产| 午夜激情福利司机影院| 国产极品精品免费视频能看的| videossex国产| 赤兔流量卡办理| 床上黄色一级片| 中出人妻视频一区二区| 色噜噜av男人的天堂激情| 国产69精品久久久久777片| 亚洲国产高清在线一区二区三| 国产av不卡久久| 精品午夜福利视频在线观看一区| 少妇被粗大猛烈的视频| 日本在线视频免费播放| 国产精品国产三级国产av玫瑰| 麻豆av噜噜一区二区三区| 午夜福利在线在线| 男人和女人高潮做爰伦理| 999久久久精品免费观看国产| 成人午夜高清在线视频| 免费av不卡在线播放| 欧美一区二区国产精品久久精品| 国产午夜精品论理片| 丰满的人妻完整版| 精品久久久久久久久久免费视频| 国产三级在线视频| 免费看日本二区| 国产av在哪里看| 男人和女人高潮做爰伦理| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放| 在线播放无遮挡| 精品午夜福利在线看| 大又大粗又爽又黄少妇毛片口| 精品一区二区免费观看| 欧美中文日本在线观看视频| 国产亚洲精品综合一区在线观看| 国产精品福利在线免费观看| 精品人妻熟女av久视频| 亚洲av二区三区四区| 搡老妇女老女人老熟妇| 亚洲中文日韩欧美视频| 一个人免费在线观看电影| 91精品国产九色| 在线观看一区二区三区| 午夜老司机福利剧场| 一区二区三区激情视频| 99在线视频只有这里精品首页| 亚洲va日本ⅴa欧美va伊人久久| 制服丝袜大香蕉在线| 日韩大尺度精品在线看网址| 国产黄a三级三级三级人| 亚洲av中文字字幕乱码综合| 国产精品免费一区二区三区在线| 嫩草影院精品99| 高清在线国产一区| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 午夜福利欧美成人| 日韩欧美在线二视频| 亚洲成人中文字幕在线播放| 欧美高清性xxxxhd video| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 最近最新免费中文字幕在线| 国产伦人伦偷精品视频| 1024手机看黄色片| 欧美日本视频| 亚洲国产欧美人成| av在线老鸭窝| 黄色日韩在线| 久久这里只有精品中国| 最新在线观看一区二区三区| 精品久久久久久,| 欧美区成人在线视频| 九色成人免费人妻av| 亚洲黑人精品在线| 亚洲无线观看免费| 国产午夜精品久久久久久一区二区三区 | 国产成人影院久久av| 日本三级黄在线观看| 亚洲欧美日韩高清专用| 级片在线观看| 国产精品伦人一区二区| 一进一出抽搐gif免费好疼| 国产成人福利小说| 国产69精品久久久久777片| 看十八女毛片水多多多| 亚洲乱码一区二区免费版| 成人av在线播放网站| 久久久久久久亚洲中文字幕| 国产大屁股一区二区在线视频| 国产主播在线观看一区二区| 亚洲黑人精品在线| 国产精品嫩草影院av在线观看 | 黄色视频,在线免费观看| 色尼玛亚洲综合影院| 国语自产精品视频在线第100页| 欧美日本亚洲视频在线播放| 日韩强制内射视频| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久亚洲 | 成人性生交大片免费视频hd| 看免费成人av毛片| 免费大片18禁| 好男人在线观看高清免费视频| 婷婷精品国产亚洲av| 国产黄a三级三级三级人| 午夜精品一区二区三区免费看| 久久久午夜欧美精品| 亚洲七黄色美女视频| 欧美高清成人免费视频www| 草草在线视频免费看| 国内精品宾馆在线| 午夜日韩欧美国产| av天堂在线播放| 99热6这里只有精品| 小说图片视频综合网站| 亚洲精华国产精华精| 麻豆成人av在线观看| 午夜福利视频1000在线观看| 午夜福利在线观看免费完整高清在 | 亚洲四区av| 欧美日韩亚洲国产一区二区在线观看| av福利片在线观看| 国产精品爽爽va在线观看网站| 久久久色成人| 老熟妇乱子伦视频在线观看| 国产一区二区激情短视频| 老司机午夜福利在线观看视频| 搡老岳熟女国产| 桃红色精品国产亚洲av| 国产探花极品一区二区| 国产精品久久久久久精品电影| 高清毛片免费观看视频网站| 婷婷精品国产亚洲av| 亚洲内射少妇av| 亚洲成人久久爱视频| 永久网站在线| 成人三级黄色视频| a在线观看视频网站| 精品久久久久久久久久久久久| 国产免费男女视频| 久久久精品欧美日韩精品| 久久久久精品国产欧美久久久| 色综合亚洲欧美另类图片| 好男人在线观看高清免费视频| 亚洲,欧美,日韩| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 九九久久精品国产亚洲av麻豆| 精品免费久久久久久久清纯| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 干丝袜人妻中文字幕| 搡老熟女国产l中国老女人| 国产成人福利小说| 在线播放国产精品三级| 丰满的人妻完整版| 国内精品久久久久久久电影| 伊人久久精品亚洲午夜| 亚洲最大成人av| 99久国产av精品| 亚洲成人免费电影在线观看| 国产欧美日韩精品亚洲av| 国产精品嫩草影院av在线观看 | 乱码一卡2卡4卡精品| 亚洲av免费在线观看| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av| 床上黄色一级片| 欧美另类亚洲清纯唯美| 精华霜和精华液先用哪个| 午夜激情欧美在线| 国产精品乱码一区二三区的特点| 午夜亚洲福利在线播放| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 午夜福利视频1000在线观看| 不卡视频在线观看欧美| 人人妻人人澡欧美一区二区| 一区二区三区免费毛片| 久久久久久久久中文| 18禁在线播放成人免费| 乱码一卡2卡4卡精品| 成人高潮视频无遮挡免费网站| а√天堂www在线а√下载| 我的老师免费观看完整版| 日韩av在线大香蕉| 久久香蕉精品热| 日本 欧美在线| 国产极品精品免费视频能看的| 嫩草影院新地址| 天天一区二区日本电影三级| 88av欧美| 可以在线观看的亚洲视频| 国产男靠女视频免费网站| 一本久久中文字幕| 亚洲av中文字字幕乱码综合| 少妇被粗大猛烈的视频| 乱系列少妇在线播放| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 两个人的视频大全免费| 精品免费久久久久久久清纯| 又黄又爽又刺激的免费视频.| 国产av不卡久久| 中出人妻视频一区二区| 国语自产精品视频在线第100页| 亚洲内射少妇av| 亚洲第一电影网av| 日日啪夜夜撸| 简卡轻食公司| av福利片在线观看| 国产精品嫩草影院av在线观看 | 天堂影院成人在线观看| 亚洲成人久久性| 在线免费十八禁| 精品福利观看| 国产精品爽爽va在线观看网站| 伊人久久精品亚洲午夜| 舔av片在线| 亚洲性久久影院| 高清日韩中文字幕在线| 国内精品久久久久精免费| 香蕉av资源在线| 啦啦啦观看免费观看视频高清| 久久婷婷人人爽人人干人人爱| 成人无遮挡网站| 亚洲中文字幕日韩| 久久久久久国产a免费观看|