• <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
    国产成人一区二区三区免费视频网站| 久久香蕉激情| 国产精品自产拍在线观看55亚洲| av在线天堂中文字幕| 久久久久久国产a免费观看| 欧美一区二区国产精品久久精品 | 欧美日韩福利视频一区二区| 夜夜躁狠狠躁天天躁| 麻豆成人av在线观看| 国产伦在线观看视频一区| 日韩中文字幕欧美一区二区| 一进一出好大好爽视频| 久久精品成人免费网站| 久久热在线av| 身体一侧抽搐| 亚洲精品美女久久久久99蜜臀| 国模一区二区三区四区视频 | 欧美日韩瑟瑟在线播放| 国产一区在线观看成人免费| 欧美成人性av电影在线观看| 成年女人毛片免费观看观看9| 午夜精品一区二区三区免费看| 亚洲精品一卡2卡三卡4卡5卡| 女人高潮潮喷娇喘18禁视频| 亚洲自拍偷在线| 中文在线观看免费www的网站 | 国产av麻豆久久久久久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲av成人一区二区三| 国产成人aa在线观看| 色综合站精品国产| 亚洲精品中文字幕一二三四区| a级毛片在线看网站| 国产成人一区二区三区免费视频网站| 欧美极品一区二区三区四区| 制服人妻中文乱码| 特级一级黄色大片| 日韩欧美国产一区二区入口| 久久这里只有精品中国| 嫩草影院精品99| 1024手机看黄色片| 国产av一区在线观看免费| 看黄色毛片网站| 99国产精品一区二区三区| 又黄又粗又硬又大视频| av免费在线观看网站| 亚洲精品一区av在线观看| 久久中文看片网| а√天堂www在线а√下载| bbb黄色大片| 欧洲精品卡2卡3卡4卡5卡区| 精品午夜福利视频在线观看一区| 亚洲天堂国产精品一区在线| 亚洲欧美日韩东京热| 国内精品一区二区在线观看| 老司机靠b影院| 国产精品av视频在线免费观看| 久久久水蜜桃国产精品网| 亚洲国产欧美一区二区综合| 精品欧美国产一区二区三| 国产区一区二久久| 欧美国产日韩亚洲一区| 久久午夜综合久久蜜桃| 国产1区2区3区精品| 757午夜福利合集在线观看| 欧美一区二区精品小视频在线| 91麻豆精品激情在线观看国产| 性色av乱码一区二区三区2| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看 | 国产精品九九99| 18美女黄网站色大片免费观看| 两人在一起打扑克的视频| 丁香六月欧美| 婷婷亚洲欧美| 国产99白浆流出| 国产主播在线观看一区二区| 搡老妇女老女人老熟妇| 女警被强在线播放| 99热这里只有是精品50| 国产精品一区二区精品视频观看| 欧美精品啪啪一区二区三区| 国产一区二区三区视频了| 少妇被粗大的猛进出69影院| 琪琪午夜伦伦电影理论片6080| 母亲3免费完整高清在线观看| 免费在线观看黄色视频的| 99久久综合精品五月天人人| 久久人妻av系列| 亚洲一区中文字幕在线| 女人高潮潮喷娇喘18禁视频| 天堂√8在线中文| www.www免费av| 亚洲一区高清亚洲精品| 久久久国产成人精品二区| 性欧美人与动物交配| 精品久久久久久成人av| 观看免费一级毛片| 99国产精品99久久久久| 亚洲色图av天堂| 天天躁夜夜躁狠狠躁躁| 久久精品夜夜夜夜夜久久蜜豆 | 久久精品国产亚洲av高清一级| 好男人在线观看高清免费视频| 亚洲男人的天堂狠狠| 精品久久久久久久人妻蜜臀av| 白带黄色成豆腐渣| 12—13女人毛片做爰片一| av欧美777| 国产亚洲精品久久久久5区| 男插女下体视频免费在线播放| 精品无人区乱码1区二区| 九色国产91popny在线| 免费在线观看影片大全网站| 欧美性猛交╳xxx乱大交人| 首页视频小说图片口味搜索| 三级男女做爰猛烈吃奶摸视频| 国产精品av视频在线免费观看| 18禁美女被吸乳视频| 91麻豆av在线| 在线a可以看的网站| 黄片小视频在线播放| 搞女人的毛片| 久久婷婷人人爽人人干人人爱| 一本久久中文字幕| 久久精品国产综合久久久| 国产真人三级小视频在线观看| 中亚洲国语对白在线视频| 久久热在线av| 村上凉子中文字幕在线| 一级片免费观看大全| 久久伊人香网站| 999久久久精品免费观看国产| 久久久久久国产a免费观看| 久久久久久亚洲精品国产蜜桃av| 亚洲 欧美 日韩 在线 免费| a级毛片a级免费在线| 精品欧美一区二区三区在线| 亚洲中文字幕一区二区三区有码在线看 | 男女午夜视频在线观看| cao死你这个sao货| 免费一级毛片在线播放高清视频| 亚洲国产精品999在线| 免费搜索国产男女视频| 99精品在免费线老司机午夜| av国产免费在线观看| 99riav亚洲国产免费| 日日夜夜操网爽| 国产男靠女视频免费网站| 大型av网站在线播放| 观看免费一级毛片| 久久人人精品亚洲av| cao死你这个sao货| 人人妻人人看人人澡| 巨乳人妻的诱惑在线观看| 日韩大码丰满熟妇| 一本大道久久a久久精品| 日本 av在线| 久久午夜亚洲精品久久| 欧美乱码精品一区二区三区| 国产精品,欧美在线| 在线观看免费午夜福利视频| 国产乱人伦免费视频| 天堂动漫精品| 精品久久久久久成人av| 男人舔奶头视频| 99国产精品一区二区蜜桃av| 精品久久久久久久人妻蜜臀av| 丁香欧美五月| 桃红色精品国产亚洲av| 欧美高清成人免费视频www| 国产精品久久久av美女十八| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 国产精品一及| 性欧美人与动物交配| 亚洲最大成人中文| 久久精品人妻少妇| 深夜精品福利| 男女视频在线观看网站免费 | 国产精品一区二区三区四区免费观看 | 18禁观看日本| 婷婷六月久久综合丁香| 久久人妻福利社区极品人妻图片| 在线观看免费午夜福利视频| 国产av一区在线观看免费| 91九色精品人成在线观看| 麻豆久久精品国产亚洲av| 男人舔女人下体高潮全视频| 国产精品av久久久久免费| 成年人黄色毛片网站| 在线播放国产精品三级| 激情在线观看视频在线高清| 亚洲国产精品999在线| 亚洲九九香蕉| 免费在线观看成人毛片| 精品久久久久久久人妻蜜臀av| 12—13女人毛片做爰片一| 美女高潮喷水抽搐中文字幕| 久久精品成人免费网站| 日韩欧美在线二视频| 免费无遮挡裸体视频| 天天躁狠狠躁夜夜躁狠狠躁| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| 一级毛片精品| 宅男免费午夜| 精品日产1卡2卡| 亚洲中文字幕一区二区三区有码在线看 | 天天添夜夜摸| 日本在线视频免费播放| 午夜福利在线在线| 91国产中文字幕| 韩国av一区二区三区四区| 99re在线观看精品视频| 久久这里只有精品19| 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 国产97色在线日韩免费| 欧美三级亚洲精品| 午夜福利免费观看在线| 午夜久久久久精精品| 女人高潮潮喷娇喘18禁视频| 国产av又大| 国产伦在线观看视频一区| 国产私拍福利视频在线观看| 欧美丝袜亚洲另类 | 国产免费男女视频| 国产黄a三级三级三级人| 欧美一级毛片孕妇| 伊人久久大香线蕉亚洲五| 午夜亚洲福利在线播放| 日韩欧美在线乱码| 男女床上黄色一级片免费看| 精品欧美一区二区三区在线| 国产一级毛片七仙女欲春2| 俄罗斯特黄特色一大片| 日韩大尺度精品在线看网址| 97人妻精品一区二区三区麻豆| 91国产中文字幕| 特级一级黄色大片| 亚洲精品美女久久av网站| 成年女人毛片免费观看观看9| 999久久久精品免费观看国产| 狠狠狠狠99中文字幕| 免费在线观看日本一区| 日韩 欧美 亚洲 中文字幕| 91国产中文字幕| 老司机深夜福利视频在线观看| 中文字幕人妻丝袜一区二区| 精品午夜福利视频在线观看一区| 2021天堂中文幕一二区在线观| 黄片大片在线免费观看| 手机成人av网站| 午夜福利免费观看在线| 精品国产乱子伦一区二区三区| 精华霜和精华液先用哪个| 欧美色视频一区免费| 欧美日韩精品网址| 少妇裸体淫交视频免费看高清 | 黄色毛片三级朝国网站| 亚洲激情在线av| 亚洲一区高清亚洲精品| 久久久久亚洲av毛片大全| 1024视频免费在线观看| 久久精品国产综合久久久| 色综合亚洲欧美另类图片| 亚洲人与动物交配视频| 精品少妇一区二区三区视频日本电影| 舔av片在线| 久久午夜综合久久蜜桃| 午夜福利欧美成人| 成年免费大片在线观看| 亚洲国产高清在线一区二区三| 黄色视频,在线免费观看| 国产成人啪精品午夜网站| 国产精品一区二区三区四区免费观看 | 日韩国内少妇激情av| 波多野结衣高清无吗| 亚洲国产欧美一区二区综合| 国产伦一二天堂av在线观看| 色综合站精品国产| 亚洲欧洲精品一区二区精品久久久| 日韩欧美 国产精品| 在线观看美女被高潮喷水网站 | 一区福利在线观看| 少妇的丰满在线观看| 我的老师免费观看完整版| 99国产精品一区二区蜜桃av| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 成人三级做爰电影| 99久久久亚洲精品蜜臀av| 欧美zozozo另类| 99久久99久久久精品蜜桃| 波多野结衣高清作品| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美精品v在线| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 亚洲一码二码三码区别大吗| 日韩欧美在线乱码| 日韩三级视频一区二区三区| 久久热在线av| 日韩欧美三级三区| 蜜桃久久精品国产亚洲av| 午夜福利在线在线| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 欧美成狂野欧美在线观看| 757午夜福利合集在线观看| videosex国产| 麻豆一二三区av精品| 久久99热这里只有精品18| 中文字幕精品亚洲无线码一区| 日韩国内少妇激情av| 日韩欧美三级三区| 最新美女视频免费是黄的| 日韩欧美精品v在线| 99久久国产精品久久久| 两个人的视频大全免费| 国产精华一区二区三区| 国产精品亚洲一级av第二区| 亚洲国产欧美一区二区综合| 成人高潮视频无遮挡免费网站| 亚洲狠狠婷婷综合久久图片| 99精品在免费线老司机午夜| 深夜精品福利| 国产人伦9x9x在线观看| 久久天堂一区二区三区四区| www.熟女人妻精品国产| 国产精品av视频在线免费观看| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| x7x7x7水蜜桃| 国产精品影院久久| 美女 人体艺术 gogo| 我的老师免费观看完整版| 国产91精品成人一区二区三区| 国产亚洲欧美98| 97碰自拍视频| 99在线视频只有这里精品首页| 一二三四在线观看免费中文在| 最近最新中文字幕大全免费视频| 少妇熟女aⅴ在线视频| 99国产综合亚洲精品| 成人精品一区二区免费| 午夜福利免费观看在线| 韩国av一区二区三区四区| 精品久久久久久久毛片微露脸| 黄色毛片三级朝国网站| 91成年电影在线观看| 国产激情欧美一区二区| 久久久国产精品麻豆| 看片在线看免费视频| 十八禁网站免费在线| 精品第一国产精品| 天堂动漫精品| 国产精品野战在线观看| 天堂√8在线中文| 怎么达到女性高潮| 久久久久久国产a免费观看| 亚洲精品av麻豆狂野| 欧美三级亚洲精品| 午夜免费激情av| 国产亚洲精品久久久久久毛片| a级毛片在线看网站| 婷婷精品国产亚洲av| 国产精品久久久久久人妻精品电影| 欧美日韩瑟瑟在线播放| 免费人成视频x8x8入口观看| 欧美精品亚洲一区二区| 这个男人来自地球电影免费观看| 亚洲国产高清在线一区二区三| 高清毛片免费观看视频网站| 身体一侧抽搐| 无人区码免费观看不卡| 免费看a级黄色片| 成年免费大片在线观看| 五月伊人婷婷丁香| 国产成人啪精品午夜网站| 成人三级黄色视频| 老司机午夜福利在线观看视频| 国产精品久久久av美女十八| 人妻久久中文字幕网| 在线观看一区二区三区| www.熟女人妻精品国产| 久久中文看片网| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久亚洲av鲁大| 日韩大码丰满熟妇| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 国产亚洲欧美98| 免费无遮挡裸体视频| 精品高清国产在线一区| 女同久久另类99精品国产91| 亚洲狠狠婷婷综合久久图片| 国产人伦9x9x在线观看| 欧美又色又爽又黄视频| 麻豆成人午夜福利视频| 欧美激情久久久久久爽电影| 国产午夜精品久久久久久| 亚洲男人的天堂狠狠| 日韩欧美国产一区二区入口| 色精品久久人妻99蜜桃| 看片在线看免费视频| 国产蜜桃级精品一区二区三区| 成年免费大片在线观看| 国产一级毛片七仙女欲春2| 性色av乱码一区二区三区2| 中文字幕久久专区| 99久久精品热视频| 亚洲欧美激情综合另类| 日韩中文字幕欧美一区二区| 日本a在线网址| 亚洲成av人片在线播放无| 亚洲性夜色夜夜综合| 久久精品国产综合久久久| 日本免费a在线| 在线十欧美十亚洲十日本专区| 成人午夜高清在线视频| 一二三四在线观看免费中文在| 亚洲五月婷婷丁香| 麻豆久久精品国产亚洲av| 欧美日韩一级在线毛片| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 欧美黄色淫秽网站| 亚洲国产精品sss在线观看| 91麻豆av在线| 免费av毛片视频| 天堂影院成人在线观看| 九色成人免费人妻av| 成人精品一区二区免费| 天天添夜夜摸| 曰老女人黄片| 18美女黄网站色大片免费观看| 在线观看66精品国产| 两性夫妻黄色片| 麻豆久久精品国产亚洲av| 麻豆成人午夜福利视频| av中文乱码字幕在线| 中出人妻视频一区二区| 草草在线视频免费看| 亚洲一区高清亚洲精品| 亚洲全国av大片| www.www免费av| 成人高潮视频无遮挡免费网站| 久久久久久大精品| 看免费av毛片| 久久天堂一区二区三区四区| av片东京热男人的天堂| 国产精品99久久99久久久不卡| 国产免费男女视频| 他把我摸到了高潮在线观看| 亚洲第一欧美日韩一区二区三区| 一进一出抽搐gif免费好疼| 久久精品成人免费网站| 欧美日韩国产亚洲二区| 免费高清视频大片| 欧美最黄视频在线播放免费| 国产欧美日韩精品亚洲av| 亚洲av成人精品一区久久| 亚洲激情在线av| 久久久久精品国产欧美久久久| 国产蜜桃级精品一区二区三区| 国产精品久久电影中文字幕| 国产蜜桃级精品一区二区三区| 日韩欧美在线二视频| 最好的美女福利视频网| 99精品久久久久人妻精品| 国产精品1区2区在线观看.| 欧美性猛交黑人性爽| 国产在线精品亚洲第一网站| 美女午夜性视频免费| 国产69精品久久久久777片 | 怎么达到女性高潮| 国产欧美日韩一区二区三| 日韩欧美国产在线观看| 少妇被粗大的猛进出69影院| www.精华液| 婷婷六月久久综合丁香| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 亚洲免费av在线视频| 亚洲美女黄片视频| 18禁黄网站禁片午夜丰满| 久久精品91蜜桃| 久久久精品欧美日韩精品| 亚洲欧美日韩高清在线视频| 少妇被粗大的猛进出69影院| 国产久久久一区二区三区| 男女做爰动态图高潮gif福利片| 精品欧美一区二区三区在线| 搡老岳熟女国产| 国产单亲对白刺激| 日本五十路高清| 免费在线观看成人毛片| 久久人人精品亚洲av| 无人区码免费观看不卡| www日本在线高清视频| 欧美一区二区国产精品久久精品 | 亚洲国产欧美网| 他把我摸到了高潮在线观看| 首页视频小说图片口味搜索| 国产精品综合久久久久久久免费| 十八禁人妻一区二区| 手机成人av网站| 欧美日韩亚洲综合一区二区三区_| 国产99久久九九免费精品| 国产精品一区二区三区四区久久| 欧美大码av| 精品国产美女av久久久久小说| 亚洲精品国产一区二区精华液| 亚洲欧美日韩无卡精品| 亚洲欧美激情综合另类| 国内少妇人妻偷人精品xxx网站 | 国产成人aa在线观看| 三级国产精品欧美在线观看 | 久久九九热精品免费| 成人国产一区最新在线观看| 特大巨黑吊av在线直播| 中文字幕人妻丝袜一区二区| 在线观看舔阴道视频| 色播亚洲综合网| 香蕉丝袜av| 曰老女人黄片| 男人的好看免费观看在线视频 | 老司机福利观看| 99精品欧美一区二区三区四区| 国产精品亚洲美女久久久| 欧美黑人巨大hd| 色哟哟哟哟哟哟| 精品欧美一区二区三区在线| 亚洲成人国产一区在线观看| 国产午夜精品论理片| 中国美女看黄片| 国产精品久久久久久亚洲av鲁大| 校园春色视频在线观看| 亚洲av美国av| 老司机深夜福利视频在线观看| 国产成人精品无人区| 久久精品亚洲精品国产色婷小说| 久久中文字幕人妻熟女| 成人特级黄色片久久久久久久| 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| 毛片女人毛片| 亚洲成人精品中文字幕电影| 精品第一国产精品| 舔av片在线| 亚洲中文字幕日韩| 精品一区二区三区四区五区乱码| 男女做爰动态图高潮gif福利片| 男男h啪啪无遮挡| 亚洲aⅴ乱码一区二区在线播放 | 精品久久蜜臀av无| 欧美黄色片欧美黄色片| 免费在线观看视频国产中文字幕亚洲| 国产亚洲av嫩草精品影院| 精品人妻1区二区| 91大片在线观看| av在线天堂中文字幕| 熟女电影av网| 亚洲免费av在线视频| 香蕉av资源在线| 国产亚洲精品综合一区在线观看 | 我的老师免费观看完整版| 老司机午夜福利在线观看视频| 亚洲精品美女久久av网站| 日韩国内少妇激情av| 99riav亚洲国产免费| 亚洲欧美日韩无卡精品| 两个人看的免费小视频| 午夜亚洲福利在线播放| 国产精品综合久久久久久久免费| 91麻豆av在线| 九色成人免费人妻av| 国产精品1区2区在线观看.| 可以在线观看的亚洲视频| 色综合亚洲欧美另类图片| 91在线观看av| 91国产中文字幕| 午夜福利欧美成人| 香蕉国产在线看| 成人国产综合亚洲| 成年人黄色毛片网站| 无遮挡黄片免费观看| 欧美一级a爱片免费观看看 | 欧美成人午夜精品| 日本撒尿小便嘘嘘汇集6| 亚洲av美国av| 少妇熟女aⅴ在线视频| svipshipincom国产片| 成人三级黄色视频| 淫妇啪啪啪对白视频| 国产伦人伦偷精品视频| 悠悠久久av| 国产av一区在线观看免费| 精品久久久久久久毛片微露脸| 午夜两性在线视频| 国产成人欧美在线观看| 深夜精品福利| 老司机在亚洲福利影院| 很黄的视频免费| 国产精品一及| 岛国在线观看网站| 91麻豆av在线| 免费高清视频大片| 757午夜福利合集在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| av欧美777| 一本综合久久免费| 一本一本综合久久|