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

    黑潮延伸體區(qū)域50~100 km渦旋分布特征

    2021-12-09 08:33:08錢思佳于方杰
    海洋科學(xué) 2021年11期
    關(guān)鍵詞:黑潮中尺度旋渦

    錢思佳, 于方杰, 陳 戈

    黑潮延伸體區(qū)域50~100 km渦旋分布特征

    錢思佳, 于方杰, 陳 戈

    (中國海洋大學(xué) 信息科學(xué)與工程學(xué)院, 山東 青島 266100)

    本文使用基于熱成風(fēng)速度的渦旋識別拓展方法, 通過海表面溫度數(shù)據(jù)對黑潮延伸體區(qū)域50~ 100 km渦旋進(jìn)行研究, 發(fā)現(xiàn)50~100 km渦旋主要分布在黑潮延伸體流軸兩側(cè), 氣旋渦和反氣旋渦的壽命、半徑分布具有一致性。氣旋渦多出現(xiàn)在35°N以北, 反氣旋渦在35°N以南比較集中, 與尺度較小的中尺度渦旋分布特征較為相似。冬夏兩季渦旋地理分布存在一定差異, 主要與不同季節(jié)該區(qū)域海表溫度梯度及風(fēng)應(yīng)力旋度的變化有關(guān)。35°N以南50~100 km渦旋數(shù)量的季節(jié)性變化與風(fēng)速大小的季節(jié)性變化存在明顯的正相關(guān)性。35°N以南50~100 km渦旋三倍半徑內(nèi)風(fēng)速異常和風(fēng)應(yīng)力旋度歸一化表明, 氣旋渦對應(yīng)風(fēng)速負(fù)異常而反氣旋渦對應(yīng)風(fēng)速正異常, 反氣旋渦的產(chǎn)生依賴于風(fēng)應(yīng)力負(fù)旋度, 氣旋渦的生成與風(fēng)應(yīng)力正旋度有關(guān)。

    黑潮延伸體; 50~100 km渦旋; 海表面溫度; 渦旋識別拓展方法; 渦旋-風(fēng)速

    黑潮延伸體(140°~180°E, 28°~40°N)是伴隨大振幅灣流和大量脫落渦旋的帶向射流, 存在豐富的多尺度動力過程。前人的研究已經(jīng)揭示了該區(qū)域中尺度渦旋的分布特征、季節(jié)性變化和生成機(jī)制[1-3]。由于渦旋運(yùn)動通常會引起局部海平面高度變化, 因此人們普遍認(rèn)為渦旋比較容易被衛(wèi)星高度計(jì)捕獲。然而, 現(xiàn)有的衛(wèi)星高度計(jì)具有相對較高的沿軌采樣分辨率但留有較大的跨軌間隙[4], 使得盡管多顆同時在軌的衛(wèi)星高度計(jì)的融合產(chǎn)品的分辨率已經(jīng)有了較大提高, 但也不足以準(zhǔn)確描述100 km以下海洋渦旋的動態(tài)變化[5]。因此, 100 km以下渦旋目前處于觀測盲點(diǎn), 其分布特征仍存在較大的認(rèn)知空白, 且目前對100 km以下渦旋的尺度劃分還存在一定爭議??紤]到黑潮延伸體區(qū)域第一斜壓羅斯貝變形半徑約為30~40 km[6], 而亞中尺度要小于所在緯度的第一斜壓羅斯貝變形半徑, 因而本文以50 km來大致區(qū)分該區(qū)域的中尺度和亞中尺度特征, 已有研究中也提出過50 km以下為亞中尺度這一觀點(diǎn)[7-8]。作為典型的中尺度(100 km以上)與亞中尺度的過渡, 50~100 km渦旋存在著重要意義。一方面, 對黑潮延伸體區(qū)域中尺度渦旋的統(tǒng)計(jì)分析表明, 半徑越小的中尺度渦旋數(shù)量越多[2, 9], 因而推斷出該區(qū)域存在豐富的50~100 km渦旋, 但卻無法被衛(wèi)星高度計(jì)準(zhǔn)確觀測。另一方面, 渦旋相互作用可以實(shí)現(xiàn)能量的正向和逆向級聯(lián), 使能量通過各種尺度在時間和空間上重新分布, 因此, 50~100 km渦旋充當(dāng)著海洋中大小尺度之間能量級聯(lián)的關(guān)鍵橋梁。

    盡管許多數(shù)值模式已經(jīng)被用來研究100 km以下渦旋在海洋熱傳輸及物質(zhì)輸運(yùn)等方面的作用[8, 10], 但相關(guān)結(jié)論仍需通過現(xiàn)場觀測加以驗(yàn)證。因此, 如何進(jìn)行有效觀測一直是近年來的研究熱點(diǎn)。水文測量儀器和浮標(biāo)[11], 岸基高頻雷達(dá)[12], 漂流浮標(biāo)[13], 水下滑翔機(jī)[14]等已用于觀測100 km以下渦旋。然而, 這些研究主要集中在局部區(qū)域的個例分析上, 所得到的渦旋特征很大程度上取決于觀測設(shè)備的數(shù)量, 而海洋設(shè)備的部署是昂貴且耗時的。與衛(wèi)星測高數(shù)據(jù)相比, 海面溫度(sea surface temperature, SST)數(shù)據(jù)具有更高的空間分辨率, 且大部分上層海洋動力過程都與SST直接或間接相關(guān)。因此, 從理論上講, SST是在全球尺度上研究50~100 km渦旋的一種可行方法。盡管SST存在易受天氣過程影響等局限性, 但其較高的時空分辨率使得在基于海表面高度(sea surface height, SSH)的渦旋檢測算法無法正常工作的情況下, 它可以作為檢測渦旋的有力補(bǔ)充。目前使用人工神經(jīng)網(wǎng)絡(luò)、準(zhǔn)輪廓跟蹤和聚類、隨機(jī)橢圓擬合與圖像分割、與渦旋結(jié)構(gòu)相關(guān)的iso-SST模式, 小波變換分析海表溫度鋒面等方法均能基于SST成功檢測出渦旋[15-19]。前人提出了一種基于熱成風(fēng)速度的渦旋識別方法[20], 僅基于渦旋幾何形狀來確定渦旋。該方法在許多地區(qū)得到了運(yùn)用, 均體現(xiàn)了不錯的效果[21-23]。

    本文使用基于熱成風(fēng)速度的渦旋識別拓展方法, 基于高分辨率SST數(shù)據(jù)對黑潮延伸體區(qū)域50~100 km渦旋分布特征進(jìn)行研究并探究其與風(fēng)速之間的相關(guān)性。

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

    1.1 數(shù)據(jù)

    本文使用Remote Sensing Systems發(fā)布的SST數(shù)據(jù)(http://www.remss.com/measurements/sea-surface- temperature/)來檢測50~100 km渦旋。該數(shù)據(jù)是結(jié)合微波的穿云功能和紅外SST數(shù)據(jù)高分辨率及近海岸觀測能力的近實(shí)時產(chǎn)品。數(shù)據(jù)從2002年6月1日發(fā)布至今, 其空間分辨率為9 km, 時間分辨率為1 d, 每天可以覆蓋全球所有的海洋[20]。

    本文使用的風(fēng)場數(shù)據(jù)來自美國國家航空航天局(NASA)的QuikSCAT衛(wèi)星(https://coastwatch.pfeg. noaa.gov/data.html#)。該衛(wèi)星上的海風(fēng)觀測儀是一種雙波束微波散射儀, 用于測量全球海洋上風(fēng)力的大小和方向, 風(fēng)速測量的參考高度為地表以上10 m, 可以覆蓋全球90%的海洋, 被認(rèn)為是海表風(fēng)速的準(zhǔn)確呈現(xiàn)[24]。數(shù)據(jù)集時間范圍為1999年7月21日至2009年11月21日, 空間分辨率為0.125°, 時間分辨率為1 d。NASA對海風(fēng)數(shù)據(jù)進(jìn)行了處理, 得到與風(fēng)速同分辨率的風(fēng)應(yīng)力數(shù)據(jù)。

    本文使用的海平面異常(sea level anomaly, SLA)數(shù)據(jù)為AVISO發(fā)布的all-sat融合產(chǎn)品(http://www. aviso.altimetry.fr)[25]。該數(shù)據(jù)融合了所有可用的衛(wèi)星高度計(jì)觀測數(shù)據(jù), 空間分辨率為0.25°, 時間分辨率為1 d, 能有效反映出全球海平面高度信息。

    本文使用的drifter數(shù)據(jù)為全球漂流浮標(biāo)計(jì)劃(global drifter program, GDP)每小時數(shù)據(jù)集的1.01版本(https://www.aoml.noaa.gov/phod/gdp/hourly_data. php), 該數(shù)據(jù)集中共包含985個由全球定位系統(tǒng)跟蹤(GPS-tracked)的drifter軌跡文件, 每個文件中記錄了每小時該drifter的位置和速度。

    1.2 方法

    本文對基于熱成風(fēng)速度的渦旋識別方法進(jìn)行了拓展, 命名為基于熱成風(fēng)速度的渦旋識別拓展方法, 該方法通過以下步驟定義渦旋:

    1) SST降噪。計(jì)算出海表面溫度異常(sea surface temperature anomaly, SSTA)以最大程度地降低由天氣和其他偶然因素引起的每日溫度變化的影響。

    SSTA=SST–DC, (1)

    其中, DC (daily climatology)為每日氣候態(tài)SST, 由18年間(2003—2020年)的每日SST數(shù)據(jù)得到。第日的DC值為第日所在月份的平均SST與每年第日的平均SST的均值(閏年最后一天的DC值用第365天的值代替)。

    2) 依據(jù)基于熱成風(fēng)速度的渦旋識別方法[20]將SSTA轉(zhuǎn)化為對應(yīng)速度場并進(jìn)行渦旋識別得到初始渦旋的渦心和極性, 渦旋邊界定義為SSTA水平梯度最大處溫度異常等值線。

    3) 空間濾波。對于所有初始渦旋, 依據(jù)其空間尺度大小(渦旋內(nèi)包含的最小網(wǎng)格數(shù)量需滿足min=8[26])進(jìn)行篩選, 去除由SST隨機(jī)噪聲產(chǎn)生的小的“偽渦旋”。

    4) 時間濾波。由于SST易受干擾, 短期存在的渦旋為噪聲引起的“偽渦旋”的可能性較大。依據(jù)相鄰時間幀的渦旋極性和位置對所有渦旋進(jìn)行連續(xù)追蹤觀測, 根據(jù)渦旋生命周期對其進(jìn)行時間濾波, 去除短期波動的影響。生命周期≥的渦旋被認(rèn)為是真實(shí)渦旋。

    =L≥2, (4)

    其中,L是第個初始渦旋的壽命(此處不考慮僅存在一天的初始渦旋),為初始渦旋總數(shù),表示初始渦旋的平均壽命,為標(biāo)準(zhǔn)差,取最接近整數(shù)。由于不同地理位置渦旋壽命和數(shù)量可能存在一定差異, 因而在不同區(qū)域計(jì)算得到的值可能會存在不同(綜合考慮保留的渦旋數(shù)量和準(zhǔn)確性, 本文中黑潮延伸體區(qū)域值選取為4)。

    2 50~100 km渦旋分布

    本文采用2002年6月1日至2018年5月31日的SST數(shù)據(jù), 使用基于熱成風(fēng)速度的渦旋識別拓展方法, 成功檢測出了大部分中尺度渦旋和50~100 km渦旋。圖1顯示了2006年1、4、7、10月份某日渦旋分布情況。在同一時間幀下, 將渦旋與對應(yīng)的SSTA及SLA在空間上進(jìn)行疊加, 結(jié)果表明大多數(shù)溫度異常都能被算法檢測為渦旋。由于本文采用的SST數(shù)據(jù)分辨率為9 km而SLA數(shù)據(jù)分辨率為0.25°, 因此一些尺度較小的渦旋不能在SLA上體現(xiàn)出來, 但算法檢測出來的半徑較大(>40 km)的渦旋附近的等高線均較為密集并伴隨著SLA極值點(diǎn)的產(chǎn)生。圖1半徑較大的渦旋中, 約75%的渦旋的渦心與SLA中心并不重合, 但這些渦旋均與SLA極值點(diǎn)最外圍等值線確定的區(qū)域有所重疊。在不同月份, 半徑較大的渦旋均能與附近的SLA匹配, 算法體現(xiàn)了不錯的檢測效果。

    圖1 渦旋分布示意圖

    注: 圖中渦旋半徑均大于20公里, 底圖為SSTA, 金線為SLA, 紅色(藍(lán)色)曲線表示反氣旋渦(氣旋渦)

    為檢驗(yàn)算法在黑潮延伸體區(qū)域的準(zhǔn)確性, 本文結(jié)合SLA和drifter數(shù)據(jù)在該區(qū)域進(jìn)行了個例驗(yàn)證, 圖2顯示了一個氣旋渦在2015年12月15日至26日的位置變化。在觀測期間內(nèi)目標(biāo)渦旋附近出現(xiàn)了一個drifter(ID: dd61504040), 從15日開始逐漸接近目標(biāo)渦旋并一直圍繞渦旋外圍運(yùn)動。由于SST易受干擾, 觀測期間內(nèi)渦旋出現(xiàn)過短暫消失, 但均能在下一時間節(jié)點(diǎn)上被檢測出來。渦旋所在位置附近等高線較為密集, 且在渦旋西北方向有閉合SLA等值線出現(xiàn)(黃色三角形), 表明SLA在該位置也能檢測出相應(yīng)渦旋。由連續(xù)時間內(nèi)SLA等值線的變化可以看到兩個渦旋逐漸分離的過程(圖西北角), 而這一過程也能在SSTA上顯示出來。渦旋與drifter和SLA的一致性表明算法在該區(qū)域的確具有較好的檢測效果。

    2002年6月1日至2018年5月31日間黑潮延伸體區(qū)域共識別出24 289個50~100 km渦旋, 氣旋渦和反氣旋渦的數(shù)量分別為13 385和10 904, 比值約為1.2。50~100 km渦旋生命周期分布如圖3a所示, 渦旋最大生命周期約為14 d, 最短生命周期為4 d(此處僅統(tǒng)計(jì)4 d及以上的渦旋), 與前人研究得到的10~100 km海洋動力過程時間變化通常持續(xù)數(shù)天到數(shù)周結(jié)論較為一致[27]。渦旋平均壽命約為5 d, 95%的渦旋生命周期在7 d以內(nèi), 相同壽命的渦旋中氣旋渦數(shù)量略多于反氣旋渦, 但兩者生命周期的分布趨勢基本一致。50~100 km渦旋半徑分布如圖3b所示, 隨半徑增大渦旋數(shù)量逐漸減少, 氣旋渦和反氣旋渦半徑的分布趨勢具有一致性。

    對每個1°×1°網(wǎng)格內(nèi)渦旋數(shù)量進(jìn)行統(tǒng)計(jì), 得到50~100 km渦旋地理分布如圖4所示。50~100 km渦旋總體上呈緯向分布, 大多沿等溫線分布于黑潮延伸體南北兩端, 32°~36°N之間渦旋數(shù)量相對較少。對該區(qū)域50~100 km渦旋極性分布(圖4b)進(jìn)行分析表明, 氣旋渦和反氣旋渦具有不明顯的斑塊狀交替分布特征。相對而言, 氣渦旋主要分布在黑潮延伸體北側(cè), 而反氣旋渦在黑潮延伸體南側(cè)相對密集。

    圖2 某氣旋渦2015年12月15日至26日的位置變化

    注: 顏色表示SSTA, 靛藍(lán)色細(xì)線表示SLA, 藍(lán)色粗線表示目標(biāo)渦旋, 黑線為當(dāng)天drifter運(yùn)動軌跡, 黃色三角形為SLA閉合等值線中心

    圖3 50~100 km渦旋生命周期和半徑分布

    盡管黑潮延伸體流軸形態(tài)年際變化較大, 但多年間其流軸位置基本位于32°~36°N之間[28]。黑潮延伸體流軸附近50~100 km渦旋數(shù)量較少而流軸兩側(cè)渦旋數(shù)量較多, 與前人利用SSH得到的渦旋分布特征并不完全相符。已有研究表明黑潮延伸體區(qū)域渦旋生成機(jī)制主要有兩種: 一是由黑潮延伸體流軸彎曲脫落生成, 二是由水平剪切力不穩(wěn)定產(chǎn)生。流軸彎曲脫落形成的渦旋主要位于流軸附近, 大多為壽命較長半徑較大攜帶能量較強(qiáng)的渦旋, 其中反氣旋渦主要分布在流軸北側(cè), 氣旋渦主要分布在流軸南側(cè)[29]。而由水平剪切力不穩(wěn)定產(chǎn)生的渦旋壽命和半徑相對較小, 這一部分渦旋中, 氣旋渦主要分布在35°N以北, 反氣旋渦主要分布在35°N以南[2]。圖4結(jié)果表明50~100 km渦旋地理分布特征與由水平剪切力不穩(wěn)定產(chǎn)生的中尺度渦旋較為相似。

    圖4 1°×1°網(wǎng)格內(nèi)50~100 km渦旋數(shù)量和極性分布

    注: 渦旋極性表示一個點(diǎn)位于氣旋渦或反氣旋渦內(nèi)的概率, 計(jì)算公式為(ae–ce)/(ae+ce)[30], 其中ae(ce)表示該點(diǎn)出現(xiàn)反氣旋渦(氣旋渦)的時間占統(tǒng)計(jì)時間的百分比, 極性為正(負(fù))表示在這一點(diǎn)上很可能會出現(xiàn)反氣旋渦(氣旋渦)

    3 50~100 km渦旋數(shù)量與風(fēng)速的相關(guān)性

    高分辨率衛(wèi)星觀測表明, 在熱帶和亞熱帶地區(qū), SST和海表風(fēng)速在中尺度上呈正相關(guān)[31], 因此本文探究50~100 km尺度上SST和海表風(fēng)速是否也存在相關(guān)性??紤]到黑潮延伸體區(qū)域50~100 km渦旋的極性分布大致以35°N為界, 圖5顯示了2002年6月1日至2009年5月31日黑潮延伸體區(qū)域35°N南北兩側(cè)的平均風(fēng)速與50~100 km渦旋生成數(shù)量的季節(jié)性變化。在35°N以北, 風(fēng)速夏季較小而冬季較大, 氣旋渦和反氣旋渦數(shù)量冬季較小而春季較大(圖5a, 5b)。風(fēng)速大小與氣旋渦和反氣旋渦數(shù)量的相關(guān)性分別為–0.2和–0.02, 表明35°N以北渦旋數(shù)量和風(fēng)速大小季節(jié)性變化幾乎不存在相關(guān)性。在35°N以南, 冬春季節(jié), 渦旋數(shù)量較多, 風(fēng)速較大; 夏秋季節(jié), 渦旋數(shù)量較少, 風(fēng)速也相對較小, 渦旋數(shù)量和風(fēng)速季節(jié)性變化具有一致性(圖5c, 5d)。風(fēng)速大小與氣旋渦和反氣旋渦數(shù)量的相關(guān)性分別為0.74和0.8, 表明35°N以南渦旋數(shù)量和風(fēng)速大小季節(jié)性變化存在較強(qiáng)的正相關(guān)。

    圖5 35°N南北兩側(cè)50~100 km渦旋數(shù)量與風(fēng)速的季節(jié)性變化

    由于渦旋數(shù)量與風(fēng)速的季節(jié)性變化一致, 進(jìn)一步分析渦旋地理分布的季節(jié)性變化與風(fēng)速之間是否存在相關(guān)性。選取夏季(6月—8月)和冬季(12月—2月)對渦旋地理分布進(jìn)行分析, 結(jié)果如圖6所示。夏季渦旋主要集中在黑潮延伸體北側(cè)(36°N以北), 南側(cè)(32°N以南)渦旋數(shù)量較少, 氣旋渦和反氣旋渦在(144°E, 39°N)附近較為集中。冬季時黑潮延伸體南側(cè)渦旋數(shù)量相對夏季明顯增多, 且整個區(qū)域渦旋生成數(shù)量更多, 位置覆蓋范圍相對更廣, 氣旋渦和反氣旋渦分布特征變化基本一致。

    圖6 夏季和冬季50~100 km渦旋地理分布

    為探究黑潮延伸體區(qū)域冬夏兩季渦旋地理分布差異的原因, 對冬夏兩季SST水平梯度及風(fēng)應(yīng)力旋度進(jìn)行分析如圖7所示。夏季SST范圍為17~27 ℃, 等溫線分布較為稀疏(圖7a), 對應(yīng)的SST水平梯度在(143°E, 36°N)和(146°E, 39°N)附近較大, 約為3.5 ℃/100 km, 其他區(qū)域較小, 35°N以南更小, 在1 ℃/100 km以下(圖7c)。對于風(fēng)應(yīng)力旋度而言, 夏季時35°N以南風(fēng)應(yīng)力為負(fù)旋度, 約為–1×10–7Pa/m, 35°N以北除SST水平梯度較大處的風(fēng)應(yīng)力旋度為2×10–7Pa/m外, 其他區(qū)域風(fēng)應(yīng)力旋度均接近于0 (圖7e)。在冬季, SST范圍為8~22 ℃, 等溫線分布較為密集(圖7b), 35°N以北SST水平梯度相較于夏季有明顯增強(qiáng), 最大值為5 ℃/100 km, 35°N以南SST水平梯度與夏季基本一致(圖7d)。冬季時風(fēng)應(yīng)力旋度相比于夏季有明顯增強(qiáng), 水平梯度較強(qiáng)(>3.5 ℃/100 km)的區(qū)域?qū)?yīng)的風(fēng)應(yīng)力旋度基本為正, 其他區(qū)域風(fēng)應(yīng)力為負(fù)旋度, 正負(fù)旋度的最大值均接近5×10–7Pa/m (圖7f)。

    渦旋冬夏兩季地理分布的差異與該區(qū)域SST水平梯度及風(fēng)應(yīng)力旋度有關(guān)。夏季時風(fēng)應(yīng)力旋度較小, SST水平梯度強(qiáng)度較弱, 上方大氣的斜壓較弱, 使得向渦動能轉(zhuǎn)換的斜壓能量較弱, 渦旋在此區(qū)域的生成較少。渦旋生成主要集中在水平梯度強(qiáng)度較大處(約為142°~146°E, 36°~40°N), 因此夏季時這一區(qū)域的渦旋數(shù)量相較于其他區(qū)域較多。由于這一區(qū)域風(fēng)應(yīng)力為正旋度, 風(fēng)向海水輸入正的渦度, 有助于氣旋渦的生成, 使得這一區(qū)域生成的渦旋中氣旋渦數(shù)量多于反氣旋渦(圖6a, b)。冬季時黑潮延伸體區(qū)域風(fēng)應(yīng)力旋度大, SST水平梯度較大, 該區(qū)域海表溫度鋒強(qiáng)度較強(qiáng), 增強(qiáng)了其上空的大氣斜壓, 使得向渦動能轉(zhuǎn)化的斜壓能量增強(qiáng), 渦旋在此區(qū)域生成更頻繁[32], 因此冬季黑潮延伸體區(qū)域渦旋數(shù)量較夏季明顯增多。較大的風(fēng)應(yīng)力負(fù)旋度使高緯度海水向南運(yùn)動, 導(dǎo)致黑潮延伸體區(qū)域南側(cè)渦旋數(shù)量增加。風(fēng)應(yīng)力負(fù)旋度給海水輸送負(fù)的渦度, 有助于反氣旋渦生成, 使得冬季反氣旋渦數(shù)量較夏季明顯增多, 相對集中在風(fēng)應(yīng)力旋度為負(fù)的黑潮延伸體南側(cè)(圖6d)。同時, 冬季較強(qiáng)的西風(fēng)使得渦旋生成后向下游運(yùn)動, 因而下游區(qū)域渦旋數(shù)量相對于上游較多。

    圖7 夏季和冬季平均SST、SST水平梯度和風(fēng)應(yīng)力旋度

    之前的分析表明黑潮延伸體區(qū)域35°N以南50~ 100 km渦旋數(shù)量與風(fēng)速大小和風(fēng)應(yīng)力旋度之間存在明顯的相關(guān)性, 為進(jìn)一步探究35°N以南50~100 km渦旋與風(fēng)速和風(fēng)應(yīng)力旋度之間的關(guān)系, 對冬夏兩季50~ 100 km渦旋3倍半徑內(nèi)風(fēng)速異常(風(fēng)速異常為每日風(fēng)速的大小與該日所在月份10 a間平均風(fēng)速值大小的差異)和渦旋生成時3倍半徑內(nèi)風(fēng)應(yīng)力旋度進(jìn)行歸一化如圖8所示。夏季反氣旋渦附近風(fēng)速異常為正值, 隨著與渦心距離的增大, 風(fēng)速異常逐漸減小; 氣旋渦附近風(fēng)速異常為負(fù)值, 渦心處風(fēng)速異常值最小, 遠(yuǎn)離渦心處風(fēng)速異常逐漸增大。冬季反氣旋渦3倍半徑內(nèi)均為強(qiáng)烈的風(fēng)速正異常, 氣旋渦3倍半徑內(nèi)均為強(qiáng)烈的風(fēng)速負(fù)異常(圖8上)。相較于夏季而言, 冬季氣旋渦和反氣旋渦附近風(fēng)速異常均明顯增強(qiáng), 這可能是由于冬季時風(fēng)向較為穩(wěn)定且SST水平梯度較大, 使得SST與風(fēng)之間的響應(yīng)增強(qiáng)。氣旋渦對應(yīng)風(fēng)速負(fù)異常而反氣旋渦對應(yīng)風(fēng)速正異常, 表明氣旋渦的溫度冷異??赡芟魅醯孛骘L(fēng)而反氣旋渦的溫度暖異常會增強(qiáng)地面風(fēng)。

    渦旋3倍半徑內(nèi)風(fēng)應(yīng)力旋度歸一化結(jié)果表明, 夏季反氣旋渦生成時附近的風(fēng)應(yīng)力幾乎均為負(fù)旋度, 氣旋渦生成時渦旋外圍出現(xiàn)了明顯的風(fēng)應(yīng)力正旋度。冬季時反氣旋渦3倍半徑內(nèi)風(fēng)應(yīng)力旋度為負(fù)且相較于夏季明顯增強(qiáng); 氣旋渦內(nèi)出現(xiàn)風(fēng)應(yīng)力正旋度, 渦旋外圍其他區(qū)域風(fēng)應(yīng)力旋度為負(fù), 正負(fù)旋度均有明顯增強(qiáng)(圖8下)。由于夏季和冬季黑潮延伸體區(qū)域35°N以南平均風(fēng)應(yīng)力旋度均為負(fù)值, 因而歸一化后氣旋和反氣旋渦附近均出現(xiàn)了風(fēng)應(yīng)力負(fù)旋度。不同的是反氣旋渦附近風(fēng)應(yīng)力幾乎均為負(fù)旋度而氣旋渦附近出現(xiàn)了明顯的正旋度, 這表明35°N以南50~100 km渦旋生成可能與風(fēng)應(yīng)力旋度有關(guān), 負(fù)旋度有助于反氣旋渦的生成, 而正旋度有助于氣旋渦的產(chǎn)生。

    圖8 50~100 km渦旋3倍半徑內(nèi)風(fēng)速異常和風(fēng)應(yīng)力旋度歸一化

    4 結(jié)論

    本文采用了基于熱成風(fēng)速度的渦旋識別拓展方法對黑潮延伸體區(qū)域50~100 km渦旋分布特征及其與風(fēng)速的相關(guān)性進(jìn)行研究, 得到如下結(jié)論:

    1) 50~100 km渦旋壽命大多在一周以內(nèi), 其生命周期及半徑分布不依賴于渦旋的極性, 其空間分布與該區(qū)域半徑較小的中尺度渦旋分布較為相似, 即氣旋渦(反氣旋渦)多分布于黑潮延伸體北(南)側(cè)。

    2) 35°N以南50~100 km渦旋數(shù)量分布具有明顯的季節(jié)性變化, 表現(xiàn)為春季較多而夏季較少; 風(fēng)速大小和氣旋渦數(shù)量、反氣旋渦數(shù)量之間均存在較強(qiáng)的正相關(guān)。在35°N以北, 渦旋數(shù)量和風(fēng)速大小之間不存在相關(guān)性。

    3) 冬夏兩季50~100 km渦旋地理分布存在一定差異, 主要與不同季節(jié)風(fēng)應(yīng)力旋度和SST水平梯度的不同有關(guān)。夏季渦旋主要分布于SST梯度較大的黑潮延伸體北側(cè)上游區(qū)域, 南側(cè)下游區(qū)域渦旋數(shù)量較少。冬季黑潮延伸體南側(cè)渦旋數(shù)量增多, 南側(cè)下游區(qū)域渦旋數(shù)量相較于夏季有明顯增加。

    4) 黑潮延伸體區(qū)域35°N以南50~100 km渦旋與風(fēng)速異常和風(fēng)應(yīng)力旋度存在相關(guān)性。氣旋渦對應(yīng)風(fēng)速負(fù)異常而反氣旋渦對應(yīng)風(fēng)速正異常, 反氣旋渦生成時其附近風(fēng)應(yīng)力旋度為負(fù)值, 氣旋渦的生成與風(fēng)應(yīng)力正旋度有關(guān)。

    [1] Itoh S, Yasuda I. Characteristics of mesoscale eddies in the Kuroshio–Oyashio Extension region detected from the distribution of the sea surface height anomaly[J]. Journal of Physical Oceanography, 2010, 40(5): 1018-1034.

    [2] Ji J, Dong C, Zhang B, et al. Oceanic eddy characteristics and generation mechanisms in the Kuroshio Extension region[J]. Journal of Geophysical Research: Oceans, 2018, 123(11): 8548-8567.

    [3] Yang Y, San Liang X. On the seasonal eddy variability in the Kuroshio Extension[J]. Journal of Physical Oceanography, 2018, 48(8): 1675-1689.

    [4] Fu L, Ubelmann C. On the transition from profile altimeter to swath altimeter for observing global ocean surface topography[J]. Journal of Atmospheric and Oceanic Technology, 2014, 31(2): 560-568.

    [5] Durand M, Fu L, Lettenmaier D P, et al. The surface water and ocean topography mission: Observing terrestrial surface water and oceanic submesoscale eddies[J]. Proceedings of the IEEE, 2010, 98(5): 766-779.

    [6] Chelton D B, Deszoeke R A, Schlax M G, et al. Geographical variability of the first baroclinic Ros-sby radius of deformation[J]. Journal of Physical Ocea-nography, 1998, 28(3): 433-460.

    [7] Zhang Z, Qiu B. Evolution of submesoscale ageos-trophic motions through the life cycle of oceanic mesoscale eddies[J]. Geophysical Research Letters, 2018, 45(21): 11847-11855.

    [8] Su Z, Wang J, Klein P, et al. Ocean submesoscales as a key component of the global heat budget[J]. Nature Communications, 2018, 9: 775.

    [9] 胡冬, 陳希, 毛科峰, 等. 黑潮延伸體鄰近區(qū)域中尺度渦特征統(tǒng)計(jì)分析[J]. 海洋與湖沼, 2018, 49(3): 497- 511.

    Hu Dong, Chen Xi, Mao Kefeng, et al. Statistical characteristics of mesoscale eddies near the Kuroshio Extension region[J]. Oceanology et Limnologia Sinica, 2018, 49(3): 497-511.

    [10] Resplandy L, Lévy M, McGillicuddy D J. Effects of eddy‐driven subduction on ocean biological carbon pump[J]. Global Biogeochemical Cycles, 2019, 33(8): 1071-1084.

    [11] Kasajima Y, Olsson K A, Johannessen T, et al. A submesoscale coherent eddy in the Greenland Sea in 2003[J]. Journal of Geophysical Research, 2006, 111(C7): C07013.

    [12] Kim S Y. Observations of submesoscale eddies using high-frequency radar-derived kinematic and dynamic quantities[J]. Continental Shelf Research, 2010, 30(15): 1639-1655.

    [13] Zheng S, Du Y, Li J, et al. Eddy characteristics in the South Indian Ocean as inferred from surface drifters[J]. Ocean Science, 2015, 11(3): 361-371.

    [14] Kolodziejczyk N, Testor P, Lazar A, et al. Subsurface fine-scale patterns in an anticyclonic eddy off Cap-Vert Peninsula observed from glider measurements[J]. Journal of Geophysical Research: Oceans, 2018, 123(9): 6312-6329.

    [15] Zhang C, Li H, Liu S, et al. Automatic detection of oceanic eddies in reanalyzed SST images and its application in the East China Sea[J]. Science China Earth Sciences, 2015, 58(12): 2249-2259.

    [16] Jin H, Yang X, Gong J, et al. Automatic eddy extraction from SST imagery using artificial neural network[C]//Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Beijing, China: ISPRSArchives, 2008, XXXVII (B6b): 279-282.

    [17] Fernandes A. Automatic water eddy detection in SST maps using random ellipse fitting and vectorial fields for image segmentation[M]. Berlin, Heidelberg: Springer Berlin Heidelberg, 2006: 77-88.

    [18] D’Alimonte D. Detection of mesoscale eddy-related structures through Iso-SST patterns[J]. IEEE Geoscience and Remote Sensing Letters, 2009, 6(2): 189-193.

    [19] Karoui I, Chauris H, Garreau P, et al. Multi- resolution eddy detection from ocean color and sea surface temperature images[C]. OCEANS'10 IEEE SYDNEY, Sydney, NSW, 2010, 1-6.

    [20] Dong C, Nencioli F, Liu Y, et al. An automated approach to detect oceanic eddies from satellite remo-tely sensed sea surface temperature data[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(6): 1055- 1059.

    [21] Liu Y, Dong C, Guan Y, 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.

    [22] Lin X, Dong C, Chen D, et al. 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, 2015, 99: 46-64.

    [23] Sun W, Dong C, Wang R, et al. Vertical structure anomalies of oceanic eddies in the Kuroshio Extension region[J]. Journal of Geophysical Research-Oceans, 2017, 122(2): 1476-1496.

    [24] Fan X, Krieger J R, Zhang J, et al. Assimilating QuikSCAT ocean surface winds with the weather research and forecasting model for surface wind-field simulation over the Chukchi/Beaufort Seas[J]. Boundary-Layer Meteorology, 2013, 148(1): 207-226.

    [25] 劉穎潔, 田豐林, 陳戈. 南海中尺度渦旋海表溫度特征統(tǒng)計(jì)研究[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2020, 50(5): 146-156.

    Liu Yingjie, Tian Fenglin, Chen Ge. Statistical characterization of sea surface temperature over meso-scale eddies in the South China Sea[J]. Periodical of Ocean University of China, 2020, 50(5): 146-156.

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

    [27] Morrow R, Fu L, Ardhuin F, et al. Global observations of fine-scale ocean surface topography with the surface water and ocean topography (SWOT) mission[J]. Frontiers in Marine Science, 2019, 6: 232.

    [28] Qiu B, Chen S. Eddy-mean flow interaction in the decadally modulating Kuroshio Extension system[J]. Deep-Sea Research Part II-Topical Studies in Oceanography, 2010, 57(13/14): 1098-1110.

    [29] 張笑, 賈英來, 沈輝, 等. 黑潮延伸體區(qū)域海洋渦旋研究進(jìn)展[J]. 氣候變化研究快報(bào), 2013, 2(1): 1-8.

    Zhang Xiao, Jia Yinglai, Shen Hui, et al. Review on mesoscale eddy studies in the Kuroshio Extension region[J]. Climate Change Research Letters, 2013, 2(1): 1-8.

    [30] Chen G, Hou Y, Chu X. Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure[J]. Journal of Geophysical Research, 2011, 116: C06018.

    [31] Chelton D B, Esbensen S K, Schlax G, et al. Observations of coupling between surface wind stress and sea surface temperature in the eastern tropical Pacific[J]. Journal of Climate, 2001, 14(7): 1479-1498.

    [32] 劉明洋, 李崇銀, 譚言科, 等. 黑潮延伸體區(qū)域海表溫度鋒的季節(jié)變化對北太平洋風(fēng)暴軸的影響[J]. 氣象, 2017, 43(4): 443-449.

    Liu Mingyang, Li Chongyin, Tan Yanke, et al. Seasonal variation of Kuroshio Extension SST front and it influence on the Pacific storm track[J]. Meteorological Monthly, 2017, 43(4): 443-449.

    Distribution characteristics of eddies with a scale of 50~100 km in the Kuroshio Extension region

    QIAN Si-jia, YU Fang-jie, CHEN Ge

    (College of Information Science and Engineering, Ocean University of China, Qingdao 266100, China)

    This paper utilizes a high-resolution sea surface temperature (SST) and an extended eddy detection algorithm based on the thermal-wind velocity to study the characteristics of eddies with a scale of 50~100 km in the Kuroshio Extension region. Results show that eddies with a scale of 50–100 km are mainly distributed on both sides of the Kuroshio Extension region. The lifetime and radius distribution of cyclonic eddies (CEs) and anticyclonic eddies (AEs) are consistent. CEs are mostly distributed in the north of 35°N while AEs are more concentrated in the south of 35°N; this distribution is similar to the distribution characteristic of mesoscale eddies at small scales. There is a slight difference in the eddy distribution characteristics in winter and summer, which is mainly due to the difference in the SST horizontal gradient and wind stress curl in different seasons. The seasonal variation of the number of CEs and AEs in the south of 35°N is positively correlated with the seasonal variation of the wind speed. The normalization of the wind speed anomaly and wind stress curl within three times the radius of eddies with a scale of 50~100 km in the south of 35°N show that AEs correspond to the positive anomaly of the wind speed while CEs correspond to the negative wind speed anomaly. AEs and CEs are generated due to the negative and positive curls of the wind stress, respectively.

    Kuroshio Extension; eddy with a scale of 50~100 km; sea surface temperature; extended eddy detection algorithm; eddy-wind speed

    Dec. 14, 2020

    P731.21

    A

    1000-3096(2021)11-0010-10

    10.11759/hykx20201214002

    2020-12-14;

    2021-03-08

    國家重點(diǎn)研發(fā)計(jì)劃(2016YFC1400904, 2016YFC1402608)

    [National key research and development program of China, No. 2016YFC1400904, No. 2016YFC1402608]

    錢思佳(1997—), 女, 湖北荊州人, 碩士生, 主要從事中小尺度渦旋分布研究, E-mail: sijiaqiancug@163.com; 于方杰(1979—),通信作者, 男, 山東青島人, 副教授, 主要從事海洋移動觀測研究, E-mail: yufangjie@ouc.edu.cn

    (本文編輯: 趙衛(wèi)紅)

    猜你喜歡
    黑潮中尺度旋渦
    譜松弛動力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動研究
    小心,旋渦來啦
    基于深度學(xué)習(xí)的中尺度渦檢測技術(shù)及其在聲場中的應(yīng)用
    大班科學(xué)活動:神秘的旋渦
    旋渦笑臉
    山間湖
    2016年7月四川持續(xù)性強(qiáng)降水的中尺度濾波分析
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    久久精品国产亚洲av天美| 亚洲av日韩在线播放| 国产精品一区www在线观看| av天堂久久9| 亚洲电影在线观看av| 99国产精品免费福利视频| 丁香六月天网| 中文字幕人妻熟人妻熟丝袜美| 99热这里只有精品一区| 精品久久久精品久久久| 18禁在线无遮挡免费观看视频| 国产精品久久久久久av不卡| 免费不卡的大黄色大毛片视频在线观看| 少妇人妻久久综合中文| 国产精品偷伦视频观看了| 午夜福利网站1000一区二区三区| 日韩电影二区| 两个人免费观看高清视频 | 最近手机中文字幕大全| 久久亚洲国产成人精品v| 国产日韩一区二区三区精品不卡 | 欧美另类一区| 99久久精品热视频| xxx大片免费视频| 蜜桃在线观看..| 美女cb高潮喷水在线观看| av国产精品久久久久影院| 国产精品一区二区三区四区免费观看| 国产亚洲欧美精品永久| 国产淫语在线视频| 亚洲伊人久久精品综合| 亚洲第一av免费看| 国产一区二区三区av在线| 国产男女超爽视频在线观看| 中文字幕av电影在线播放| 欧美国产精品一级二级三级 | 亚洲精品乱码久久久久久按摩| 亚洲av在线观看美女高潮| 六月丁香七月| 美女国产视频在线观看| 纵有疾风起免费观看全集完整版| 国产伦理片在线播放av一区| 99热这里只有精品一区| 大香蕉97超碰在线| 9色porny在线观看| 中文资源天堂在线| 99久久精品国产国产毛片| www.av在线官网国产| 美女xxoo啪啪120秒动态图| 18禁裸乳无遮挡动漫免费视频| 日韩人妻高清精品专区| 日韩,欧美,国产一区二区三区| 女的被弄到高潮叫床怎么办| 老司机影院成人| 欧美国产精品一级二级三级 | 夫妻午夜视频| 自拍偷自拍亚洲精品老妇| 免费大片黄手机在线观看| 国产一区有黄有色的免费视频| 久久久久精品性色| 丰满迷人的少妇在线观看| 最近最新中文字幕免费大全7| 熟女电影av网| 91在线精品国自产拍蜜月| 精品国产一区二区久久| 香蕉精品网在线| 在线观看av片永久免费下载| 人妻系列 视频| 婷婷色综合www| 人人妻人人澡人人看| 欧美日韩一区二区视频在线观看视频在线| 成人美女网站在线观看视频| 国产成人aa在线观看| 久久6这里有精品| 午夜激情久久久久久久| 精品国产一区二区三区久久久樱花| 日本91视频免费播放| 久久人人爽av亚洲精品天堂| 99国产精品免费福利视频| 婷婷色综合www| 男人舔奶头视频| 亚洲成人一二三区av| 3wmmmm亚洲av在线观看| 乱人伦中国视频| 国产成人91sexporn| 日韩电影二区| 视频中文字幕在线观看| 少妇高潮的动态图| 一区二区三区免费毛片| 九九在线视频观看精品| 建设人人有责人人尽责人人享有的| 91午夜精品亚洲一区二区三区| 人妻系列 视频| 丰满人妻一区二区三区视频av| 久久国产精品大桥未久av | 国产黄色视频一区二区在线观看| 亚洲欧洲国产日韩| 亚洲成人一二三区av| 天堂8中文在线网| 麻豆成人av视频| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 欧美97在线视频| 另类精品久久| 国产日韩欧美视频二区| 黄色毛片三级朝国网站 | 最近手机中文字幕大全| 一级片'在线观看视频| 日韩中文字幕视频在线看片| 99热国产这里只有精品6| 国内少妇人妻偷人精品xxx网站| 欧美97在线视频| 中文字幕亚洲精品专区| 99热6这里只有精品| 2021少妇久久久久久久久久久| 亚洲欧洲日产国产| 搡女人真爽免费视频火全软件| 自拍偷自拍亚洲精品老妇| 大话2 男鬼变身卡| 国产熟女欧美一区二区| 国产精品.久久久| 一级毛片电影观看| 一级片'在线观看视频| 国产真实伦视频高清在线观看| 久久国内精品自在自线图片| 久久韩国三级中文字幕| 亚洲精品国产成人久久av| 国产一区二区三区av在线| 亚洲美女黄色视频免费看| 大香蕉久久网| 日韩一区二区视频免费看| 亚洲国产毛片av蜜桃av| 午夜激情福利司机影院| 亚洲欧美清纯卡通| 国产黄片视频在线免费观看| 国产综合精华液| 一级毛片久久久久久久久女| 国产精品女同一区二区软件| av在线播放精品| 男女国产视频网站| 黄色一级大片看看| 亚洲av成人精品一区久久| 免费黄频网站在线观看国产| 观看免费一级毛片| 涩涩av久久男人的天堂| 高清毛片免费看| 国产精品麻豆人妻色哟哟久久| 国产伦理片在线播放av一区| av视频免费观看在线观看| 久久久久久久久久久丰满| 人妻系列 视频| 人妻少妇偷人精品九色| 五月天丁香电影| 热re99久久国产66热| 日韩强制内射视频| 99热网站在线观看| 亚洲综合色惰| 亚洲美女视频黄频| 欧美bdsm另类| 亚洲真实伦在线观看| 91午夜精品亚洲一区二区三区| a级毛片在线看网站| 人妻夜夜爽99麻豆av| 你懂的网址亚洲精品在线观看| 久久婷婷青草| 国产日韩欧美在线精品| 国产精品国产三级国产专区5o| 国产精品偷伦视频观看了| 黑丝袜美女国产一区| av不卡在线播放| 久久99热这里只频精品6学生| 久久国产精品大桥未久av | 成年人免费黄色播放视频 | 亚洲精品国产成人久久av| 97在线视频观看| 亚洲av电影在线观看一区二区三区| 七月丁香在线播放| 少妇人妻 视频| av在线观看视频网站免费| 国产精品99久久久久久久久| 国产精品一区www在线观看| 日韩av免费高清视频| 精品少妇黑人巨大在线播放| 亚洲四区av| 最近最新中文字幕免费大全7| 亚洲经典国产精华液单| 久久精品国产鲁丝片午夜精品| 日日摸夜夜添夜夜爱| 亚洲av成人精品一二三区| 国产成人aa在线观看| 大码成人一级视频| 寂寞人妻少妇视频99o| 美女中出高潮动态图| 一级av片app| 激情五月婷婷亚洲| 欧美精品一区二区大全| 成人综合一区亚洲| 国产色婷婷99| 极品少妇高潮喷水抽搐| 久久精品国产亚洲av天美| 中文字幕亚洲精品专区| videos熟女内射| 欧美一级a爱片免费观看看| 一个人免费看片子| 黑人巨大精品欧美一区二区蜜桃 | 大香蕉久久网| av天堂中文字幕网| av福利片在线观看| 欧美精品人与动牲交sv欧美| 寂寞人妻少妇视频99o| 欧美精品一区二区大全| 寂寞人妻少妇视频99o| 成人毛片60女人毛片免费| 成人黄色视频免费在线看| 国产亚洲5aaaaa淫片| av专区在线播放| 中文资源天堂在线| 秋霞伦理黄片| 国产老妇伦熟女老妇高清| 精品亚洲成a人片在线观看| 青青草视频在线视频观看| 啦啦啦视频在线资源免费观看| 久久精品久久久久久噜噜老黄| 看十八女毛片水多多多| 久久女婷五月综合色啪小说| 久久久午夜欧美精品| 视频区图区小说| 亚洲av.av天堂| 日韩av免费高清视频| 国产视频首页在线观看| 欧美性感艳星| 观看美女的网站| 成人无遮挡网站| 女性被躁到高潮视频| 国产乱来视频区| 日本黄色片子视频| 成人二区视频| 午夜精品国产一区二区电影| 十八禁高潮呻吟视频 | 性高湖久久久久久久久免费观看| 中国美白少妇内射xxxbb| 免费看光身美女| 国产精品一区二区在线观看99| 成人亚洲欧美一区二区av| 国产成人精品无人区| 汤姆久久久久久久影院中文字幕| 观看av在线不卡| 国产精品国产三级国产av玫瑰| 亚洲内射少妇av| √禁漫天堂资源中文www| 欧美少妇被猛烈插入视频| 久久青草综合色| 精品少妇内射三级| 日本av手机在线免费观看| 22中文网久久字幕| 日韩一本色道免费dvd| 午夜影院在线不卡| 中文字幕免费在线视频6| 麻豆乱淫一区二区| 97在线视频观看| 一本大道久久a久久精品| 久久久久网色| 婷婷色综合大香蕉| 国产色爽女视频免费观看| 中文字幕制服av| 久久久久久人妻| 中国美白少妇内射xxxbb| 中国美白少妇内射xxxbb| 成人亚洲精品一区在线观看| 亚洲国产精品一区二区三区在线| 欧美精品人与动牲交sv欧美| 国产精品成人在线| 99久久精品热视频| 一级毛片电影观看| 91精品伊人久久大香线蕉| 日本欧美国产在线视频| 女性被躁到高潮视频| 亚洲av综合色区一区| 中文乱码字字幕精品一区二区三区| 国产高清有码在线观看视频| tube8黄色片| av女优亚洲男人天堂| 久久精品久久精品一区二区三区| 精品一区二区免费观看| 成人黄色视频免费在线看| 亚洲国产av新网站| 多毛熟女@视频| 少妇猛男粗大的猛烈进出视频| 在线观看免费高清a一片| 亚洲四区av| 中文字幕免费在线视频6| 久久av网站| 看免费成人av毛片| 两个人的视频大全免费| 男女边摸边吃奶| 黑人高潮一二区| 少妇猛男粗大的猛烈进出视频| 免费观看无遮挡的男女| 成人黄色视频免费在线看| 多毛熟女@视频| 老司机影院成人| 欧美日韩精品成人综合77777| 伦理电影大哥的女人| 中国美白少妇内射xxxbb| av有码第一页| 寂寞人妻少妇视频99o| av免费观看日本| 2018国产大陆天天弄谢| 国产免费视频播放在线视频| 丝袜脚勾引网站| 久久久a久久爽久久v久久| 赤兔流量卡办理| 全区人妻精品视频| 深夜a级毛片| 夫妻午夜视频| 精品少妇内射三级| 欧美日韩av久久| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 制服丝袜香蕉在线| 亚洲天堂av无毛| 久久久久视频综合| 91久久精品国产一区二区三区| 91精品伊人久久大香线蕉| 亚洲精品日本国产第一区| 如何舔出高潮| 男女无遮挡免费网站观看| 欧美区成人在线视频| 精品人妻偷拍中文字幕| 国产极品粉嫩免费观看在线 | 欧美精品亚洲一区二区| 国产精品一区www在线观看| 亚洲不卡免费看| 91成人精品电影| 人妻一区二区av| 国产精品人妻久久久久久| 熟女电影av网| 国产精品秋霞免费鲁丝片| 青春草亚洲视频在线观看| 久久久久久久久久久久大奶| 久久精品久久久久久噜噜老黄| 日本色播在线视频| 99久久精品一区二区三区| 免费看光身美女| 免费黄网站久久成人精品| 高清视频免费观看一区二区| 久久久久久久久久人人人人人人| 欧美 日韩 精品 国产| 91成人精品电影| 亚洲精品国产av成人精品| 欧美 日韩 精品 国产| 免费人成在线观看视频色| 久久人妻熟女aⅴ| 少妇被粗大猛烈的视频| 黑人猛操日本美女一级片| 亚洲精品国产av成人精品| 女的被弄到高潮叫床怎么办| 少妇的逼水好多| 制服丝袜香蕉在线| 亚洲综合色惰| 性色avwww在线观看| 国产无遮挡羞羞视频在线观看| 欧美变态另类bdsm刘玥| 日本免费在线观看一区| 日韩精品有码人妻一区| 青青草视频在线视频观看| 亚洲欧美中文字幕日韩二区| 亚洲人成网站在线播| 熟女av电影| 蜜臀久久99精品久久宅男| 自拍欧美九色日韩亚洲蝌蚪91 | 超碰97精品在线观看| 亚洲人成网站在线观看播放| 永久免费av网站大全| 成人国产av品久久久| 久久久欧美国产精品| 精品人妻一区二区三区麻豆| 寂寞人妻少妇视频99o| 99久久中文字幕三级久久日本| 少妇人妻精品综合一区二区| 久久精品熟女亚洲av麻豆精品| 精品熟女少妇av免费看| 精品一区在线观看国产| 欧美成人午夜免费资源| 交换朋友夫妻互换小说| 2018国产大陆天天弄谢| 男人爽女人下面视频在线观看| 国产乱来视频区| 免费黄网站久久成人精品| 国产在线免费精品| 最黄视频免费看| 91久久精品国产一区二区三区| av不卡在线播放| 中文资源天堂在线| 精品国产乱码久久久久久小说| av国产久精品久网站免费入址| 精品久久久精品久久久| 男人和女人高潮做爰伦理| 国产免费又黄又爽又色| 人妻制服诱惑在线中文字幕| 午夜精品国产一区二区电影| 成人免费观看视频高清| 青春草视频在线免费观看| 亚洲国产日韩一区二区| 亚洲美女黄色视频免费看| 日本黄色日本黄色录像| 99视频精品全部免费 在线| 99久久精品一区二区三区| 欧美最新免费一区二区三区| 亚洲av男天堂| 五月天丁香电影| 午夜视频国产福利| 我的老师免费观看完整版| 高清av免费在线| 一级黄片播放器| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 一本久久精品| 国产高清三级在线| 久久99一区二区三区| 2018国产大陆天天弄谢| 亚洲人与动物交配视频| kizo精华| 欧美bdsm另类| 大片免费播放器 马上看| 日本91视频免费播放| 熟女电影av网| 一级片'在线观看视频| 久久人人爽人人片av| 亚洲无线观看免费| 一区在线观看完整版| 午夜91福利影院| 日本欧美视频一区| 国产极品天堂在线| 久久97久久精品| 国产爽快片一区二区三区| 国产精品一区二区在线观看99| 亚州av有码| 最近的中文字幕免费完整| 99re6热这里在线精品视频| 亚洲av男天堂| 在线 av 中文字幕| 国产成人午夜福利电影在线观看| 亚洲天堂av无毛| 亚洲精品国产av成人精品| 99热这里只有是精品50| 国产免费又黄又爽又色| 国产极品天堂在线| 最后的刺客免费高清国语| 国产精品不卡视频一区二区| 亚洲一区二区三区欧美精品| 女人久久www免费人成看片| 成人国产av品久久久| 国产 一区精品| 嫩草影院入口| 黄色欧美视频在线观看| 2022亚洲国产成人精品| 国产黄片视频在线免费观看| 老熟女久久久| 国产女主播在线喷水免费视频网站| 久久ye,这里只有精品| 女的被弄到高潮叫床怎么办| 免费av不卡在线播放| 国产亚洲欧美精品永久| 免费看av在线观看网站| 午夜福利,免费看| av天堂久久9| 九九爱精品视频在线观看| 中文字幕人妻丝袜制服| 亚洲精品国产色婷婷电影| 亚洲经典国产精华液单| 观看av在线不卡| 99热6这里只有精品| 久久久久视频综合| 永久免费av网站大全| 久久精品久久精品一区二区三区| 美女视频免费永久观看网站| 亚洲av日韩在线播放| 日韩欧美精品免费久久| 日韩,欧美,国产一区二区三区| 国产成人精品婷婷| 国产深夜福利视频在线观看| 哪个播放器可以免费观看大片| 777米奇影视久久| 菩萨蛮人人尽说江南好唐韦庄| 久久人人爽人人片av| 大话2 男鬼变身卡| 黄片无遮挡物在线观看| 老司机影院成人| 成人毛片60女人毛片免费| 女性被躁到高潮视频| 美女大奶头黄色视频| 人妻少妇偷人精品九色| 精品人妻偷拍中文字幕| 精品久久久精品久久久| 日本免费在线观看一区| 热99国产精品久久久久久7| 永久网站在线| 大香蕉97超碰在线| 国产免费福利视频在线观看| 一本一本综合久久| 成年美女黄网站色视频大全免费 | 亚洲精品色激情综合| 亚洲国产毛片av蜜桃av| 亚洲自偷自拍三级| 日韩视频在线欧美| 啦啦啦在线观看免费高清www| av在线app专区| 视频区图区小说| 亚洲美女搞黄在线观看| 26uuu在线亚洲综合色| 18禁动态无遮挡网站| 观看av在线不卡| 国产精品久久久久久精品古装| 亚洲精品,欧美精品| 亚洲成人一二三区av| 观看av在线不卡| 蜜桃久久精品国产亚洲av| 日韩精品有码人妻一区| 中文字幕免费在线视频6| 免费播放大片免费观看视频在线观看| 久热久热在线精品观看| 3wmmmm亚洲av在线观看| 中文字幕亚洲精品专区| 日本wwww免费看| 自线自在国产av| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 日韩中文字幕视频在线看片| 国产欧美日韩精品一区二区| 精品一区在线观看国产| 两个人免费观看高清视频 | 人妻 亚洲 视频| av天堂久久9| 欧美激情国产日韩精品一区| 观看美女的网站| 成人美女网站在线观看视频| 中文欧美无线码| 麻豆乱淫一区二区| 日本色播在线视频| 亚洲国产欧美日韩在线播放 | 国精品久久久久久国模美| 亚洲中文av在线| 久久韩国三级中文字幕| 国产高清不卡午夜福利| 国产免费一区二区三区四区乱码| 大香蕉久久网| 国产av国产精品国产| 精品酒店卫生间| 亚洲婷婷狠狠爱综合网| 波野结衣二区三区在线| 精品国产国语对白av| 99视频精品全部免费 在线| 黄色怎么调成土黄色| 久久久精品免费免费高清| 一本—道久久a久久精品蜜桃钙片| 如何舔出高潮| 最新的欧美精品一区二区| 日韩亚洲欧美综合| 制服丝袜香蕉在线| 成年人午夜在线观看视频| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影小说| 亚洲第一区二区三区不卡| 午夜免费男女啪啪视频观看| 亚洲av男天堂| 久久精品国产a三级三级三级| 最新中文字幕久久久久| 免费看日本二区| 国产成人精品久久久久久| 亚洲欧美日韩卡通动漫| 五月开心婷婷网| 99热6这里只有精品| 99re6热这里在线精品视频| 偷拍熟女少妇极品色| 麻豆成人av视频| 七月丁香在线播放| 久久国产精品男人的天堂亚洲 | 王馨瑶露胸无遮挡在线观看| 亚洲成色77777| 丰满少妇做爰视频| 国产欧美亚洲国产| 免费观看无遮挡的男女| 久久精品国产亚洲av涩爱| 男女边吃奶边做爰视频| av视频免费观看在线观看| 免费久久久久久久精品成人欧美视频 | 一区在线观看完整版| 啦啦啦啦在线视频资源| 亚洲精品亚洲一区二区| 久久精品夜色国产| 天天躁夜夜躁狠狠久久av| 亚洲国产色片| 纵有疾风起免费观看全集完整版| 国产亚洲精品久久久com| 你懂的网址亚洲精品在线观看| 一区二区三区免费毛片| 免费av不卡在线播放| 91精品国产国语对白视频| 伦理电影大哥的女人| 91精品国产国语对白视频| av天堂中文字幕网| 久久午夜福利片| 99精国产麻豆久久婷婷| 国产免费一区二区三区四区乱码| 免费av不卡在线播放| 国模一区二区三区四区视频| 精品国产一区二区三区久久久樱花| 精品少妇内射三级| 久久久久久久久久久免费av| 国产成人a∨麻豆精品| 久久女婷五月综合色啪小说| 大片免费播放器 马上看| 国产白丝娇喘喷水9色精品| 国产精品久久久久久av不卡|