• <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
    乱码一卡2卡4卡精品| 天堂中文最新版在线下载| 亚洲欧美精品专区久久| 精品久久久久久久久av| 久久午夜福利片| 丝袜在线中文字幕| 中文资源天堂在线| 国产精品无大码| 国产男女内射视频| 国产精品99久久99久久久不卡 | 亚洲在久久综合| 三级经典国产精品| 伊人久久精品亚洲午夜| 99re6热这里在线精品视频| 亚洲性久久影院| 亚洲,欧美,日韩| 国产av码专区亚洲av| 中文天堂在线官网| 日韩精品有码人妻一区| 成人特级av手机在线观看| 99久久精品国产国产毛片| 国产无遮挡羞羞视频在线观看| 国产成人精品一,二区| 99re6热这里在线精品视频| av网站免费在线观看视频| 国产精品无大码| 久久久久久伊人网av| 久久6这里有精品| 丰满饥渴人妻一区二区三| 人人妻人人看人人澡| 成人毛片a级毛片在线播放| 亚洲精品色激情综合| 日韩欧美精品免费久久| 国产高清三级在线| 两个人的视频大全免费| 色哟哟·www| 下体分泌物呈黄色| 午夜免费男女啪啪视频观看| 自拍偷自拍亚洲精品老妇| 成人特级av手机在线观看| a 毛片基地| 中文字幕久久专区| 久久免费观看电影| 中国国产av一级| 麻豆精品久久久久久蜜桃| 赤兔流量卡办理| 在线观看av片永久免费下载| 我要看日韩黄色一级片| 超碰97精品在线观看| 欧美+日韩+精品| 国产黄频视频在线观看| 天堂8中文在线网| 建设人人有责人人尽责人人享有的| 成人午夜精彩视频在线观看| 自拍偷自拍亚洲精品老妇| 欧美性感艳星| 超碰97精品在线观看| 99re6热这里在线精品视频| 亚洲熟女精品中文字幕| 亚洲高清免费不卡视频| 久久久久久久亚洲中文字幕| 在线观看av片永久免费下载| 亚洲精品国产成人久久av| 51国产日韩欧美| 久热这里只有精品99| 天天躁夜夜躁狠狠久久av| 2018国产大陆天天弄谢| 久久精品国产亚洲网站| 国产欧美亚洲国产| 亚洲综合色惰| 国产 精品1| 久热这里只有精品99| 91午夜精品亚洲一区二区三区| 一级毛片电影观看| 婷婷色综合www| 高清欧美精品videossex| 黄片无遮挡物在线观看| 免费高清在线观看视频在线观看| 十分钟在线观看高清视频www | 久久久国产一区二区| 中文精品一卡2卡3卡4更新| 在线观看美女被高潮喷水网站| 黄色毛片三级朝国网站 | 99re6热这里在线精品视频| 夫妻性生交免费视频一级片| av.在线天堂| 免费人妻精品一区二区三区视频| 99久久精品热视频| 欧美少妇被猛烈插入视频| 天堂8中文在线网| 老司机亚洲免费影院| 亚洲欧美中文字幕日韩二区| 免费看av在线观看网站| 乱码一卡2卡4卡精品| 色5月婷婷丁香| 精品熟女少妇av免费看| 看非洲黑人一级黄片| 亚洲精品乱久久久久久| 丰满少妇做爰视频| h视频一区二区三区| 一级,二级,三级黄色视频| 我要看日韩黄色一级片| 精品亚洲乱码少妇综合久久| 中文字幕人妻丝袜制服| 黑人猛操日本美女一级片| 99热6这里只有精品| 亚洲国产精品999| av免费观看日本| 国产免费视频播放在线视频| 国产色婷婷99| 街头女战士在线观看网站| 久久综合国产亚洲精品| 精品亚洲乱码少妇综合久久| 亚洲精品亚洲一区二区| 久久精品久久精品一区二区三区| 亚洲av成人精品一区久久| 人妻夜夜爽99麻豆av| 一级毛片电影观看| av线在线观看网站| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 国产日韩欧美在线精品| 欧美日韩在线观看h| 91精品国产国语对白视频| 亚洲国产精品国产精品| 日韩一本色道免费dvd| 丰满迷人的少妇在线观看| 久热久热在线精品观看| 免费观看av网站的网址| 国产精品免费大片| 国产精品伦人一区二区| 欧美xxⅹ黑人| av天堂久久9| 丝袜在线中文字幕| 国产免费一级a男人的天堂| 一级毛片电影观看| 亚洲精品一区蜜桃| 日日啪夜夜撸| 日韩欧美精品免费久久| 日韩在线高清观看一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 男人和女人高潮做爰伦理| 成年女人在线观看亚洲视频| 五月伊人婷婷丁香| 亚洲精品日本国产第一区| 午夜免费男女啪啪视频观看| 美女国产视频在线观看| 免费黄色在线免费观看| 欧美最新免费一区二区三区| 男女国产视频网站| 久久久久久人妻| 欧美xxⅹ黑人| 国精品久久久久久国模美| 丁香六月天网| 自拍偷自拍亚洲精品老妇| 美女福利国产在线| 赤兔流量卡办理| 最黄视频免费看| 亚洲va在线va天堂va国产| 日韩一区二区视频免费看| 免费高清在线观看视频在线观看| 中文字幕精品免费在线观看视频 | 久久久国产欧美日韩av| 国内少妇人妻偷人精品xxx网站| 一级二级三级毛片免费看| 亚洲美女搞黄在线观看| 22中文网久久字幕| 六月丁香七月| 18+在线观看网站| 免费播放大片免费观看视频在线观看| 在现免费观看毛片| 久久久久国产网址| 男女边吃奶边做爰视频| 蜜桃久久精品国产亚洲av| 午夜福利在线观看免费完整高清在| 免费av中文字幕在线| 亚洲精品456在线播放app| 欧美日韩综合久久久久久| 精品久久久久久久久av| 日日啪夜夜爽| 91aial.com中文字幕在线观看| 能在线免费看毛片的网站| 国产在线男女| 丰满人妻一区二区三区视频av| 99热全是精品| 制服丝袜香蕉在线| 久久韩国三级中文字幕| 国产伦理片在线播放av一区| 美女国产视频在线观看| 亚洲国产欧美日韩在线播放 | 免费在线观看成人毛片| 在线天堂最新版资源| 国产亚洲最大av| 国产爽快片一区二区三区| 亚洲精品aⅴ在线观看| 亚洲情色 制服丝袜| 国产深夜福利视频在线观看| 日韩av在线免费看完整版不卡| 日日摸夜夜添夜夜爱| 99热这里只有是精品50| av福利片在线| 秋霞在线观看毛片| 精品熟女少妇av免费看| 精品少妇内射三级| 亚洲无线观看免费| 免费久久久久久久精品成人欧美视频 | 欧美精品一区二区大全| 久久人人爽av亚洲精品天堂| 男女无遮挡免费网站观看| 日本色播在线视频| 日韩 亚洲 欧美在线| 性高湖久久久久久久久免费观看| 日韩强制内射视频| 日韩视频在线欧美| 日本wwww免费看| 国产美女午夜福利| 肉色欧美久久久久久久蜜桃| 国产成人免费无遮挡视频| 亚洲丝袜综合中文字幕| 视频区图区小说| 最近的中文字幕免费完整| 性高湖久久久久久久久免费观看| 欧美另类一区| 国产色婷婷99| 女人精品久久久久毛片| 成年av动漫网址| 日韩av不卡免费在线播放| 王馨瑶露胸无遮挡在线观看| 国产在线视频一区二区| 噜噜噜噜噜久久久久久91| 久久精品国产鲁丝片午夜精品| 亚洲av中文av极速乱| 美女内射精品一级片tv| 国产探花极品一区二区| 男女国产视频网站| 久久狼人影院| 色婷婷av一区二区三区视频| 晚上一个人看的免费电影| 人妻夜夜爽99麻豆av| 久久国产亚洲av麻豆专区| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄| 狂野欧美白嫩少妇大欣赏| 最近中文字幕2019免费版| 国产综合精华液| 日本欧美视频一区| 桃花免费在线播放| 免费黄网站久久成人精品| 2022亚洲国产成人精品| 简卡轻食公司| 亚洲av日韩在线播放| 日日爽夜夜爽网站| 欧美精品高潮呻吟av久久| 日韩中字成人| 成人二区视频| 只有这里有精品99| 免费黄色在线免费观看| 欧美精品人与动牲交sv欧美| 2021少妇久久久久久久久久久| 伊人久久精品亚洲午夜| 亚洲情色 制服丝袜| 一个人免费看片子| 极品人妻少妇av视频| 国产69精品久久久久777片| 天堂中文最新版在线下载| 久久 成人 亚洲| 国产精品国产av在线观看| 日日撸夜夜添| 纯流量卡能插随身wifi吗| a级毛片在线看网站| 国产在视频线精品| 国产高清不卡午夜福利| 一级毛片久久久久久久久女| 久久毛片免费看一区二区三区| 久久久久视频综合| 涩涩av久久男人的天堂| 国产一区二区三区av在线| 人妻少妇偷人精品九色| 男女免费视频国产| 亚洲精品国产色婷婷电影| 丝袜在线中文字幕| 老女人水多毛片| 伊人久久国产一区二区| 日日啪夜夜撸| 久久 成人 亚洲| 人妻 亚洲 视频| 观看av在线不卡| 2021少妇久久久久久久久久久| 老熟女久久久| 国产在线视频一区二区| 久久国产乱子免费精品| 伦理电影免费视频| 中文字幕制服av| 亚洲欧美一区二区三区黑人 | 妹子高潮喷水视频| 色网站视频免费| 精品一区在线观看国产| 久久精品国产鲁丝片午夜精品| 欧美精品人与动牲交sv欧美| 久久 成人 亚洲| 久久久欧美国产精品| 午夜av观看不卡| 一本久久精品| 日韩 亚洲 欧美在线| 激情五月婷婷亚洲| 日韩一区二区三区影片| 成人亚洲欧美一区二区av| 久久久国产一区二区| 性色av一级| 亚洲国产日韩一区二区| 日日撸夜夜添| 好男人视频免费观看在线| 精品亚洲成国产av| 国产精品一区二区性色av| 亚洲精品国产色婷婷电影| 国产精品蜜桃在线观看| 国产精品一区二区在线观看99| 亚洲怡红院男人天堂| 精品久久久精品久久久| 免费人成在线观看视频色| 欧美亚洲 丝袜 人妻 在线| 美女xxoo啪啪120秒动态图| 香蕉精品网在线| 看十八女毛片水多多多| 十分钟在线观看高清视频www | 一个人免费看片子| 熟女电影av网| 国模一区二区三区四区视频| av福利片在线观看| 丰满少妇做爰视频| 最近中文字幕2019免费版| 在线观看免费高清a一片| 亚洲性久久影院| 久久久久人妻精品一区果冻| 久久亚洲国产成人精品v| 午夜久久久在线观看| 亚洲av成人精品一区久久| 亚洲av.av天堂| 欧美成人午夜免费资源| 极品人妻少妇av视频| 自线自在国产av| 婷婷色综合大香蕉| 欧美日韩亚洲高清精品| 视频区图区小说| 久久久久久久国产电影| 草草在线视频免费看| 国产精品久久久久久精品古装| 蜜臀久久99精品久久宅男| 少妇人妻精品综合一区二区| 亚洲精品自拍成人| 久久人人爽人人片av| 久久久亚洲精品成人影院| 人妻少妇偷人精品九色| 国产永久视频网站| 久久精品国产亚洲网站| 欧美精品一区二区免费开放| 久热久热在线精品观看| 日韩精品有码人妻一区| 高清在线视频一区二区三区| 免费不卡的大黄色大毛片视频在线观看| videos熟女内射| 日日撸夜夜添| 久久国产精品男人的天堂亚洲 | 最新中文字幕久久久久| 成年av动漫网址| 日韩不卡一区二区三区视频在线| 国产乱人偷精品视频| 国产伦精品一区二区三区视频9| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 中文字幕制服av| 久久精品国产自在天天线| 极品教师在线视频| 久久精品久久久久久噜噜老黄| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| 成年av动漫网址| 亚洲欧美成人综合另类久久久| 高清毛片免费看| 免费播放大片免费观看视频在线观看| 精品人妻一区二区三区麻豆| 一本色道久久久久久精品综合| 少妇精品久久久久久久| 啦啦啦在线观看免费高清www| 亚洲第一av免费看| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 久久久久人妻精品一区果冻| 国产成人freesex在线| 久久精品久久精品一区二区三区| 波野结衣二区三区在线| 欧美 亚洲 国产 日韩一| 亚洲第一区二区三区不卡| 日韩亚洲欧美综合| 欧美精品人与动牲交sv欧美| 青春草国产在线视频| 蜜臀久久99精品久久宅男| 欧美丝袜亚洲另类| 深夜a级毛片| 久久久久久久大尺度免费视频| 日韩精品有码人妻一区| 观看免费一级毛片| 精品久久久久久久久av| 亚洲精华国产精华液的使用体验| 人人妻人人添人人爽欧美一区卜| 成人亚洲精品一区在线观看| 中文字幕精品免费在线观看视频 | 国产av精品麻豆| 最近手机中文字幕大全| 成人毛片60女人毛片免费| 亚洲欧美一区二区三区国产| 欧美精品一区二区免费开放| 高清毛片免费看| 色吧在线观看| 久久狼人影院| 欧美性感艳星| 国产黄色免费在线视频| 成人二区视频| www.色视频.com| 久久久久久久亚洲中文字幕| 欧美另类一区| 免费看光身美女| 热re99久久国产66热| 亚洲成人av在线免费| 亚洲精品国产av成人精品| 国产免费一级a男人的天堂| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 自线自在国产av| 十分钟在线观看高清视频www | 久久久久久伊人网av| 免费少妇av软件| 色哟哟·www| 尾随美女入室| 卡戴珊不雅视频在线播放| 国产熟女欧美一区二区| 欧美精品亚洲一区二区| 如日韩欧美国产精品一区二区三区 | 青青草视频在线视频观看| 国产欧美亚洲国产| 中文字幕精品免费在线观看视频 | 久久久久国产网址| 一本久久精品| 亚洲精品自拍成人| 纯流量卡能插随身wifi吗| 国产成人免费观看mmmm| 国产欧美日韩精品一区二区| 偷拍熟女少妇极品色| av有码第一页| 中文字幕精品免费在线观看视频 | 国产成人精品久久久久久| 狂野欧美白嫩少妇大欣赏| 日本av免费视频播放| 亚洲人成网站在线播| 日日撸夜夜添| 少妇的逼水好多| 在线观看美女被高潮喷水网站| h日本视频在线播放| 免费观看的影片在线观看| 中文天堂在线官网| 亚洲成色77777| 成人无遮挡网站| 男女免费视频国产| 国产乱来视频区| 久久精品久久久久久噜噜老黄| 国产亚洲av片在线观看秒播厂| 一级a做视频免费观看| 亚洲国产精品999| 毛片一级片免费看久久久久| 久久99一区二区三区| 亚洲欧美日韩另类电影网站| 少妇被粗大猛烈的视频| 国产日韩欧美亚洲二区| 日本猛色少妇xxxxx猛交久久| 午夜福利,免费看| 午夜福利视频精品| 我的老师免费观看完整版| 久久99热6这里只有精品| 少妇丰满av| 免费少妇av软件| 成人漫画全彩无遮挡| 亚洲经典国产精华液单| 国产精品久久久久久精品电影小说| 日产精品乱码卡一卡2卡三| 国产91av在线免费观看| 国产成人精品婷婷| 中文资源天堂在线| 王馨瑶露胸无遮挡在线观看| 美女大奶头黄色视频| 婷婷色av中文字幕| 女人精品久久久久毛片| 午夜免费男女啪啪视频观看| 美女福利国产在线| 成人毛片a级毛片在线播放| 一级av片app| 最近最新中文字幕免费大全7| 久热这里只有精品99| 国产高清有码在线观看视频| 久久人妻熟女aⅴ| 少妇熟女欧美另类| 国产精品秋霞免费鲁丝片| 免费人成在线观看视频色| 午夜免费鲁丝| 全区人妻精品视频| 99热6这里只有精品| 一区二区av电影网| 久热这里只有精品99| 中文字幕av电影在线播放| 日本-黄色视频高清免费观看| 久久人人爽人人片av| 日本wwww免费看| 99久久精品一区二区三区| 久久6这里有精品| 久久99一区二区三区| 国产精品久久久久久精品古装| 中国三级夫妇交换| 日韩在线高清观看一区二区三区| 精品卡一卡二卡四卡免费| 国内少妇人妻偷人精品xxx网站| 女的被弄到高潮叫床怎么办| av又黄又爽大尺度在线免费看| 国产精品不卡视频一区二区| 一级毛片我不卡| 国产极品天堂在线| 又爽又黄a免费视频| 欧美少妇被猛烈插入视频| 国产白丝娇喘喷水9色精品| 曰老女人黄片| 3wmmmm亚洲av在线观看| 成年女人在线观看亚洲视频| 久久国产亚洲av麻豆专区| 男的添女的下面高潮视频| 亚洲精品国产av蜜桃| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲va在线va天堂va国产| 国产亚洲5aaaaa淫片| 97在线人人人人妻| 成人漫画全彩无遮挡| 高清午夜精品一区二区三区| 欧美3d第一页| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av免费高清在线观看| 亚洲精品日韩av片在线观看| 中文字幕久久专区| 国产免费视频播放在线视频| 日韩av不卡免费在线播放| 国产在线一区二区三区精| 激情五月婷婷亚洲| 色网站视频免费| 青春草国产在线视频| 97超视频在线观看视频| 亚洲在久久综合| 国产老妇伦熟女老妇高清| 少妇被粗大猛烈的视频| 久久久精品免费免费高清| 亚洲婷婷狠狠爱综合网| 欧美精品人与动牲交sv欧美| 国产深夜福利视频在线观看| 久久影院123| 久久人妻熟女aⅴ| 晚上一个人看的免费电影| 一级a做视频免费观看| 久久女婷五月综合色啪小说| 亚洲精品第二区| 最近最新中文字幕免费大全7| 国产精品99久久久久久久久| 国语对白做爰xxxⅹ性视频网站| 黄片无遮挡物在线观看| 欧美97在线视频| 日日爽夜夜爽网站| 亚洲精品自拍成人| 丝袜喷水一区| 国精品久久久久久国模美| 男女啪啪激烈高潮av片| 涩涩av久久男人的天堂| 黄色怎么调成土黄色| 高清黄色对白视频在线免费看 | 精品久久国产蜜桃| 久久久久久久久久久免费av| 国产日韩欧美视频二区| av又黄又爽大尺度在线免费看| 九草在线视频观看| 午夜激情福利司机影院| 亚洲国产精品专区欧美| 久久亚洲国产成人精品v| 热99国产精品久久久久久7| 看非洲黑人一级黄片| 免费看不卡的av| 2021少妇久久久久久久久久久| 精品国产乱码久久久久久小说| 久久久欧美国产精品| 国产亚洲精品久久久com| 日韩伦理黄色片| 国产成人午夜福利电影在线观看| 我要看黄色一级片免费的| 最黄视频免费看| 亚洲精华国产精华液的使用体验| 2021少妇久久久久久久久久久| 国产欧美日韩综合在线一区二区 | 一边亲一边摸免费视频| 国产高清不卡午夜福利| 亚洲成人一二三区av| 亚洲内射少妇av| 欧美最新免费一区二区三区| 亚洲成人一二三区av| 爱豆传媒免费全集在线观看| 国产成人午夜福利电影在线观看| 91在线精品国自产拍蜜月| 色网站视频免费| 亚洲国产精品成人久久小说| 久久99蜜桃精品久久| 久久久久久久久久久丰满| 国产精品免费大片| 黑丝袜美女国产一区| 如何舔出高潮|