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

    南海上層海洋次中尺度過程空間差異和季節(jié)變化特征

    2021-11-02 01:08:14楊瀟霄曹海錦經(jīng)志友
    熱帶海洋學報 2021年5期
    關鍵詞:斜壓鋒面渦度

    楊瀟霄, 曹海錦, 經(jīng)志友

    1. 河海大學海洋學院, 江蘇 南京 210024;

    2. 熱帶海洋環(huán)境國家重點實驗室(中國科學院南海海洋研究所), 廣東 廣州 510301

    海洋環(huán)流中包含復雜的多尺度動力過程。以往對于渦旋的研究主要側重于大尺度和中尺度過程(地轉平衡條件)(Ferrari et al, 2009), 尤其近幾十年來, 衛(wèi)星高度計的成功應用大力推動了海洋中尺度過程的研究(Stammer, 1997; Chelton et al, 2007)。然而由于目前觀測能力在空間分辨率上的限制, 海洋中的次中尺度[空間尺度O(0.1~10km), 時間尺度O(1d)]過程難以被準確的觀測。盡管近年來, 海洋的水色衛(wèi)星、合成孔徑雷達(Synthetic Aperture Radar,SAR)圖像能夠捕捉到次中尺度過程的表面形態(tài)結構(Zeng et al, 2014; Zheng, 2017; Yu et al, 2018), 但僅通過瞬時的圖像難以用于研究其復雜的動態(tài)特征和形成機制等問題。隨著觀測和數(shù)值模擬水平的提高, 大量研究表明次中尺度過程對海洋中能量正向串級有著重要作用(McWilliams et al, 2001; Capet et al, 2008), 同時由于其非地轉過程產(chǎn)生的垂向速度較中尺度過程大一個量級, 對海洋中熱量、鹽度、營養(yǎng)鹽等物質(zhì)的垂向交換也起到了至關重要的作用(e.g., Rudnick et al, 1996; Spall et al, 2000; Lapeyre et al, 2006; Mahadevan et al, 2006)。

    次中尺度流具有較強的局地渦度(接近甚至超過行星渦度)以及較小的Richardson數(shù)[O(1)], 表明其同時具有地轉與非地轉的特征, 且易于發(fā)生剪切不穩(wěn)定和斜壓不穩(wěn)定。次中尺度流的形成主要受到背景流和強迫場的影響, 動力學上的主要生成機制包括: 1) 鋒生作用(Thomas et al, 2008); 2) 非強迫作用, 比如非地轉斜壓不穩(wěn)定(Molemaker et al,2005; Boccaletti et al, 2007); 3) 外強迫作用, 如浮力損失(Haine et al, 1998), 因此往往具有較大的水平浮力梯度, 且隨著斜壓性的增強, 尤其容易引起位渦的變化。然而, 在實際海洋中往往難以辨別和診斷其形成機制。

    南海是西太平洋最大的邊緣海, 由于其獨特的地形特征、受季風影響顯著、黑潮入侵等特點, 有著豐富的中尺度渦旋和鋒面(或渦絲)結構。這些特征是次中尺度過程形成的背景條件, 因此南海是天然的研究次中尺度過程的實驗海域(Liu et al, 2010;羅士浩 等, 2016; Cao et al, 2019; 黃小龍 等, 2020;Lin et al, 2020; Zhang et al, 2020; Zheng et al, 2020)。其中, 部分學者重點關注了南海次中尺度過程的季節(jié)性差異, 尤其在南海北部, 表現(xiàn)出典型的“冬強夏弱”的特征(Dong et al, 2018; Lin et al, 2020; Zhang et al, 2020)。而在空間上, 20°N以北地區(qū)的次中尺度過程普遍強于南部海域(Li et al, 2019)。該區(qū)域冬季豐富的次中尺度過程與混合層不穩(wěn)定以及冬季高應變密切相關(Dong et al, 2018; Zhang et al, 2020), 同時該海域次中尺度渦旋與黑潮和島嶼的相互作用也存在聯(lián)系(Lin et al, 2020)。這些次中尺度過程能夠有效地從地轉流中汲取能量, 并對垂向速度產(chǎn)生影響,增強上層海洋混合, 使得混合層深度加深(Liu et al,2010; 羅士浩 等, 2016; 鄭瑞璽 等, 2018)。通過對該區(qū)域中尺度渦旋實驗觀測到的反氣旋渦的橫斷面進行微觀結構觀測, Zhong等(2017)發(fā)現(xiàn)次中尺度過程可能導致垂向輸送且強度至少比中尺度渦旋大一個量級。Zhang等(2016)通過對該反氣旋渦進行能量收支分析, 提出次中尺度過程的生成可能是該區(qū)域反氣旋渦耗散的主要機制。更進一步研究指出,不平衡的次中尺度過程可能是該渦旋周邊高混合層湍流耗散率的主要成因(Yang et al, 2017)。而針對南海西部夏季上升流海域, 黃小龍等(2020)通過數(shù)值模擬結合高分辨率衛(wèi)星資料發(fā)現(xiàn)該海域存在豐富次中尺度現(xiàn)象, 并指出沿鋒面射流風應力引起的浮力損失可能是導致鋒面發(fā)生次中尺度對稱不穩(wěn)定的重要機制之一。Zheng等(2020)綜述了關于南海次中尺度過程的研究, 并將整個南海的次中尺度過程分為次中尺度波(長內(nèi)波、內(nèi)潮、近慣性波、不穩(wěn)定波等)、次中尺度渦旋、次中尺度陸架過程(如鋒面、離岸射流等)和次中尺度湍流4個大類, 同時發(fā)現(xiàn)針對南海北部次中尺度過程的研究較多, 而針對南部海域的研究較少, 目前對于南海整個海域的次中尺度過程的特征和機制尚不明確。本文基于高分辨率的區(qū)域海洋模型ROMS(Regional Ocean Modeling System)模式結果, 分析了南海4個子區(qū)域(北部海域、中部海域、西部海域、南部海域)次中尺度過程的季節(jié)變化特征, 并對比了不同海域的差異。通過對位渦和能量轉化的分析, 初步揭示了不同海區(qū)上層海洋次中尺度過程的控制因素和形成機制。

    1 模式配置和研究方法

    1.1 模式配置

    為了研究南海海域次中尺度特征的季節(jié)性變化,采用ROMS模式對南海海域進行模擬, 模式使用網(wǎng)格優(yōu)化的單向在線嵌套方案。低水平分辨率模式的空間分辨率約為7.5km, 覆蓋了整個西北太平洋海域(圖1a), 經(jīng)過20年啟動計算達到統(tǒng)計穩(wěn)定狀態(tài)后,為子區(qū)域高分辨率模式提供初始場和邊界條件, 高分辨率模式水平網(wǎng)格約為1.5km (圖1a中的黑色實線區(qū)域)。模擬區(qū)域為整個南海海域(2°12′—25°36′N,105°—127°E), 時間長度為1年。模式從7.5km嵌套到1.5km (5: 1)相對于多數(shù)3: 1的嵌套原則稍大, 但仍在ROMS手冊的建議范圍之內(nèi)。經(jīng)過對模式結果的檢驗發(fā)現(xiàn), 在嵌套邊界處有一定的影響, 因而在子區(qū)域模式中, 采用了較大的外邊界, 以避免邊界帶來的影響(圖1b)(Dauhajre et al, 2017)。模式在垂直方向上使用sigma坐標, 共分60層, 在上層海洋垂向加密, 上層200m為第25層。模式使用的表層大氣強迫(風應力、熱通量、淡水通量)來源于日平均的快速散射計衛(wèi)星(Quick Scatterometer,QuikSCAT)氣候態(tài)數(shù)據(jù)集以及綜合海洋大氣數(shù)據(jù)集(International Comprehensive Ocean Atmosphere Data Set, ICOADS)數(shù)據(jù)。最大區(qū)域模擬的邊界和初始條件來自于月平均的簡單海洋數(shù)據(jù)同化的(Simple Ocean Data Assimilation, SODA)海洋氣候態(tài)數(shù)據(jù)。次網(wǎng)格垂向通量和混合采用垂直混合參數(shù)化(K-profile parameterization, KPP)方案。模擬結果(溫度、表層流以及中尺度渦旋能量水平等)已經(jīng)通過衛(wèi)星測量和有限的現(xiàn)場觀測進行檢驗(黃小龍 等, 2020)。

    圖1 水平分辨率為7.5km的最大模擬區(qū)域的水深分布情況(a)和水平分辨率1.5km的模擬子區(qū)域(b)圖a中的方框為圖b區(qū)域。圖b中的方框表示不同的研究區(qū)域; 紅色虛線表示南海諸島歸屬范圍線; 黑色虛線表示200m等深線。R1表示北部研究海域; R2表示中部研究海域; R3表示西部研究海域; R4表示南部研究海域。該圖基于國家測繪地理信息局標準地圖服務網(wǎng)站下載的審圖號為GS(2016)1665的標準地圖制作Fig. 1 Water depth in the largest model domain with a horizontal resolution of 7.5 km (a), and in the sub-region with a horizontal resolution of 1.5 km (b)

    由于本文主要關注上層海洋的次中尺度過程(Tandon et al, 1995; Boccaletti et al, 2007), 因此本文主要研究上層200m的次中尺度過程, 子區(qū)域的劃分如圖1b中黑框所示。需要指出, 本文在模式中并未考慮潮汐強迫, 因而無法對內(nèi)波和次中尺度過程的相互作用進行討論。這是因為內(nèi)波在一定程度上會加強次中尺度過程的形成, 并且由于部分內(nèi)波和次中尺度的時空尺度接近, 會給模擬數(shù)據(jù)的分析帶來一定的難度(Torres et al, 2018; Cao et al, 2019)。

    1.2 次中尺度特征參數(shù)

    垂向相對渦度量級達到甚至超越行星渦度是次中尺度過程的重要特征, 其計算公式為:

    式中,u、v為水平背景速度場(單位: m·s–1)。垂向相對渦度(ζ)的絕對值越大, 次中尺度渦旋特征也就越明顯, 次中尺度過程也就越活躍(Qiu et al, 2014)。此外, 流場應變率(S)公式為:

    以及水平方向散度(δ)公式為:

    以上變量都能在一定程度上反應次中尺度過程中局地流場的強烈變化。水平應變率能夠表征流場水平變形的強弱, 強烈的水平流場拉伸往往暗示了強烈的鋒生過程(Thomas et al, 2008), 加速傾斜的鋒面容易導致次中尺度不穩(wěn)定。水平散度對于滿足體積守恒的不可壓縮流體而言與垂向速度梯度-相匹配, 式中w表示垂向速度(單位: m·s–1), 表征了垂向速度場在垂直方向的變化, 同時也表現(xiàn)了流場的輻聚/輻散。為方便比較, 這些特征參數(shù)均被除以科氏參數(shù)f進行標準化處理。

    位勢渦度(Potential vorticity, PV)作為重要的動力學示蹤劑, 是判斷海洋環(huán)流以及多尺度動力過程穩(wěn)定性的重要依據(jù)之一(Hoskins, 1974; Thomas et al,2013)。對于科氏參數(shù)f為正的北半球而言, 位渦為負時容易引發(fā)多種不穩(wěn)定, 特別是對稱不穩(wěn)定(Thomas et al, 2013)。這種不穩(wěn)定與上層海洋混合層強烈的水平浮力梯度以及在外強迫作用下產(chǎn)生的邊界浮力損失密切相關。本文采用旋轉分層流體內(nèi)的完全厄特爾(Ertel)位渦公式:

    并將其分為水平斜壓分量和垂向分量:

    此外, 本文還將從正壓、斜壓兩個方面對次中尺度動能的轉換效率隨時間、深度的變化規(guī)律進行探討。為此, 需定義次中尺度的尺度范圍。由于在低緯度地區(qū)(如南部海域)科氏參數(shù)f較小, 經(jīng)典的混合層斜壓變形半徑L=NH/f(Boccaletti et al, 2007)較大(式中H表示混合層深度), 由此定義的次中尺度過程的空間尺度偏大, 且會隨著混合層變化而變化。因此, 本文采用了基于平衡運動和不平衡運動動能相對關系的譜方法(Cao et al, 2021)定義次中尺度的范圍, 簡述如下。根據(jù)公式(7)可以將動能譜密度分解為渦度貢獻部分與散度貢獻部分(Torres et al,2018):

    2 結果

    2.1 相對渦度、水平散度以及應變率

    為了進一步分析各區(qū)域次中尺度過程的季節(jié)性變化, 圖3和表1給出了不同區(qū)域上層200m水平區(qū)域平均的相對渦度ζ、水平散度δ以及應變率S的標準化結果隨時間、深度變化情況以及它們在不同區(qū)域上層50m的季節(jié)平均值(由于ζ/f與δ/f存在正負值, 因而取其均方根值表征強度)。

    表1 不同區(qū)域特征參數(shù)上層50m各季節(jié)平均的結果Tab. 1 Seasonal averages of characteristic parameters of different regions in the upper 50 m

    圖3 R1 (a, e, i)、R2 (b, f, j)、R3 (c, g, k)和R4 (d, h, l)各區(qū)域內(nèi)平均水平相對渦度ζ/f均方根(a~d), 水平散度δ/f均方根(e~h), 水平應變率S/f (i~l)隨時間、深度變化情況圖中黑線表示該區(qū)域平均混合層深度隨時間變化情況Fig. 3The variations of averaged ζ/f, δ/f, and S/f over time and depth in different regions. The black curve shows the variation of the averaged mixed-layer depth over time

    整體來看, 相對渦度ζ/f均方根存在著明顯的季節(jié)性變化, 且隨深度明顯變小, 高值區(qū)域主要在混合層以內(nèi), 表明南海上層的次中尺度過程主要集中在混合層內(nèi)。散度δ/f均方根與應變率S/f隨時間變化情況與ζ/f均方根基本一致。不同的是,δ/f均方

    根與S/f在垂直方向上存在一個雙極值結構, 尤其是散度, 其在表層與混合層底部附近分別存在一個極值, 最大值出現(xiàn)在表層, 這可能與表層與混合層底接近為0的垂向速度有關。表明在南海次中尺度過程活躍的時期, 這些過程引起的海水輻聚(或輻散)的加強或減弱以及流場的水平變形均主要集中在很淺的表層和混合層底部, 而非混合層中間。

    不同區(qū)域結果對比來看, 在R1區(qū)域,ζ/f均方根的最大值出現(xiàn)在冬季(2月中旬), 最小值出現(xiàn)在夏季(6月末—7月初)。隨時間變化情況與混合層類似,特別是在混合層較深的冬季, 說明混合層不穩(wěn)定可能是這個區(qū)域次中尺度過程的重要控制機制, 但是兩者之間存在約半個月的滯后。此外, 由于該區(qū)域位于呂宋海峽附近, 冬季受黑潮入侵影響很強(趙偉等, 2007), 容易激發(fā)較強的流場側向剪切從而使相對渦度變強。自5月末到夏季結束, R1區(qū)域ζ/f均方根始終小于0.3, 且隨深度變化很小, 顯然在該時間段內(nèi)次中尺度過程明顯較弱。整體變化情況與Zhang等(2020)的結果一致。

    R2區(qū)域大部分特征與R1相似, 但ζ/f均方根的最大值出現(xiàn)在1月末—2月初, 最小值則是在6月中旬。此外, 不同于R1, 在5月末開始到夏季結束的時間段內(nèi), R2區(qū)域次中尺度過程雖然同樣弱于冬季, 但不同于始終保持一個很低水平的R1區(qū)域, R2區(qū)域在6—10月強度明顯大于R1區(qū)域, 且在8月中旬有一個較弱的極大值。該區(qū)域的混合層深度變化情況存在相似特征, 但是整體位相相比ζ/f均方根提前約一個月, 表明該海域次中尺度的季節(jié)性變化同樣受到明顯的混合層的滯后影響。

    R3區(qū)域不同于前兩個區(qū)域, 該區(qū)域次中尺度過程在6—10月份最為顯著。ζ/f均方根的第一個極大值出現(xiàn)在7月中旬, 這似乎與夏季西南季風的爆發(fā)有關; 第二個極大值出現(xiàn)在9月末, 主要成因可能是中南半島沿岸部分冷渦受到南沙群島影響, 產(chǎn)生的地形尾渦和西南季風的共同作用。在12月下旬到3月中下旬稍弱, 這一特征與該時間段較深的混合層相吻合, 表明該時段可能受混合層不穩(wěn)定的影響比較明顯, 之后在4—5月份出現(xiàn)極小值。值得注意的是R3區(qū)域次中尺度過程在夏秋季節(jié)活躍, 但混合層相對較淺, 說明該時段混合層不穩(wěn)定的影響并不占主導地位。

    R4區(qū)域包含南沙群島, 與其他海域明顯不同,其混合層冬季最淺而夏季最深, 8月中旬達到最大值。次中尺度過程與混合層深度變化的關系并不顯著, 有兩個相對高值分別出現(xiàn)在10月上旬和1月下旬, 最小值則在12月, 季節(jié)性差異相對較小。該區(qū)域相比于其他區(qū)域整體羅斯貝數(shù)(ζ/f)偏大, 這主要與該地區(qū)豐富的島嶼地形以及較小的科氏參數(shù)有關。島嶼地形引起的流場水平方向剪切, 易于產(chǎn)生地形尾渦(Dong et al, 2007), 激發(fā)次中尺度過程, 也使得該區(qū)域流場應變率、散度明顯強于其他區(qū)域。另外, 與其他區(qū)域不同的是, R4區(qū)域δ/f均方根、S/f的變化情況與ζ/f均方根的變化并不一致, 兩者大值均出現(xiàn)在1月和8月。

    表1列出了標準化后的相對渦度、水平散度的區(qū)域均方根值及應變率在不同季節(jié)的平均結果。對比發(fā)現(xiàn): 相對渦度ζ/f均方根在R1區(qū)域表現(xiàn)出了冬強夏弱的特點, 在冬季達到夏季的兩倍以上, 其他兩個參數(shù)的季節(jié)變化也類似; R2區(qū)域ζ/f均方根的最大平均值出現(xiàn)在冬季, 最小值出現(xiàn)在秋季, 而δ/f均方根、S/f的最小值則出現(xiàn)在夏季, 3個參數(shù)在夏秋季的差距不大; R3區(qū)域ζ/f均方根在秋季最大, 夏季略低于秋季, 春季最小, 該區(qū)域S/f分布情況與ζ/f均方根一致, 但δ/f的均方根最大值出現(xiàn)在夏季; R4區(qū)域ζ/f均方根在各季節(jié)的平均值均保持在0.7左右, 普遍高于其他3個區(qū)域,S/f的特征與相對渦度一致, 但是散度δ/f均方根在冬季的平均值明顯高于其他3個季節(jié)。下文將對各區(qū)域次中尺度特征參數(shù)強度的季節(jié)差異的原因做進一步探討。

    由于除了R4區(qū)域, 其他3個區(qū)域的次中尺度過程均在冬、夏兩季表現(xiàn)出極大或是極小值, 表明次中尺度過程在冬夏季差異顯著。至于R4區(qū)域,ζ/f均方根隨時間變化不明顯, 且δ/f均方根與S/f的兩個極大值分別出現(xiàn)在冬夏兩季, 因此本文研究重點主要集中在冬夏兩季。為了進一步量化對比各參數(shù)冬夏季的差異, 考察了不同區(qū)域各參數(shù)在1月和7月的概率密度分布函數(shù)(Probability density distribution function, PDF), 如圖4。

    圖4 R1 (a, e, i)、R2 (b, f, j)、R3 (c, g, k)和R4 (d, h, l)各區(qū)域1月和7月份水平相對渦度ζ/f (a~d), 水平散度δ/f (e~h), 水平應變率S/f (i~l)的概率密度分布曲線相對渦度圖像中的數(shù)字為不同月份的PDF的偏態(tài)系數(shù)Fig. 4 PDF curves of ζ/f, δ/f, and S/f in January and July.The numbers in the relative vorticity images are the skewness coefficients of PDF in different months

    PDF分布顯示, 各區(qū)域的δ/f均很好地符合了高斯分布, 且在各自次中尺度過程較強的月份相對零線更為分散。應變率S/f的PDF圖像近似滿足χ分布, 在次中尺度過程相對活躍的月份, 曲線峰值位置距離零點更遠, 表明此時S/f的眾數(shù)較大, 表層流體應變普遍更強。相對渦度ζ/f, 平均值在0附近,整體上看近似滿足高斯分布, 但存在不對稱的情況,偏態(tài)系數(shù)結果表明次中尺度過程活躍時期這種不對稱現(xiàn)象更為明顯, 各區(qū)域偏態(tài)系數(shù)均大于1, 此時可以較明顯地看到曲線在正半軸存在一個較長的“尾巴”, 圖像在0線附近的集中程度也較弱, 表明存在稀疏的強正相對渦度結構。

    圖5為ζ/f的偏態(tài)系數(shù)與標準差隨深度變化情況, 在混合層以上各區(qū)域均以正偏態(tài)為主。不同區(qū)域偏態(tài)系數(shù)隨深度變化存在差異, 但整體上先減小后增大, 至200m深度以下開始趨于0。ζ/f的標準差則隨深度減小, 不同月份逐漸趨于一致, 在混合層底部附近存在一個比較明顯的突變, 次中尺度過程活躍時期尤為明顯。不同月份標準差在表層的大小分布與均方根一致, 進一步證明次中尺度過程活躍時期ζ/f的分布更為分散。

    圖5 R1 (a, e)、R2 (b, f)、R3 (c, g)和R4 (d, h)各區(qū)域相對渦度ζ/f的偏態(tài)系數(shù)與標準差隨深度變化情況Fig. 5 Variations of skewness coefficient and standard deviation of ζ/f with depth in each region

    2.2 位渦分析

    基于位渦在冬季(以1月15日為代表)、夏季(以7月15日為代表) 10m水深處的水平分布(圖7)和典型斷面S1、S2、S3 (S1斷面跨越了R1區(qū)域冬季的強烈渦旋結構; S2斷面跨越了R3區(qū)域夏季強鋒面結構; S3斷面位于冬季R4區(qū)域, 受到地形影響顯著)的垂向分布(圖8), 對不同區(qū)域兩季的位渦分布特征進行分析。并基于公式(5)和公式(6)將位勢渦度q分成垂向和水平斜壓兩個分量, 對比其季節(jié)性差異(垂向分量qvert由垂向絕對渦度和浮力頻率控制, 斜壓分量qbc則受浮力的水平梯度和流場的垂向剪切影響)。

    圖6 R1 (a, e)、R2 (b, f)、R3 (c, g)和R4 (d, h)各區(qū)域1月(a~d)和7月(e~h)水平應變率S/f和相對渦度ζ/f的聯(lián)合概率密度分布圖中灰色虛線表示S /f=Fig. 6 Joint probability density distributions of S/f, ζ/f in different regions in January and July

    圖7 1月15日(a~c)和7月15日(d~f)位渦q (a, d)及其分量qvert (b, e)、qbc (c, f)在表層10m深度的水平分布情況圖中黑色實線S1、S2、S3為所選典型斷面, 圖中紅色虛線表示南海諸島歸屬范圍線。審圖號為GS(2016)1665Fig. 7 Horizontal distributions of q and its components in different seasons at 10 m. The time in (a~c) and (d~f) is consistent with that in Fig. 2. S1, S2, and S3 are the typical sections selected

    圖8 典型斷面S1 (a~c)、S2 (d~f)、S3 (h~j)處位渦q (a, d, h)及其分量qvert (b, e, i)、qbc (c, f, j)的分布情況圖中黑色虛線為該斷面混合層深度, 黑色實線為等密度線; 圖h~j中的灰色部分表示地形Fig. 8 Distributions of q and its components at S1, S2, and S3. The dashed line shows the mixed-layer depth at each section, and the solid contours are isopyrdensity lines

    如圖7所示, 整體上各區(qū)域表層位渦q空間分布存在差異, 在鋒面中尺度渦邊緣等次中尺度過程活躍海域更容易出現(xiàn)負值, 容易產(chǎn)生不穩(wěn)定現(xiàn)象。同時, 垂向分量qvert在這些區(qū)域有著較大的絕對值。斜壓分量qbc則存在明顯的正負變化, 這與北半球流體在準地轉平衡態(tài)下qbc始終為負的特點不符, 表明此時局地非地轉作用不可忽視, 流體垂向速度的水平剪切十分重要。與qvert類似,qbc在鋒面、次中尺度渦旋、渦絲及中尺度渦邊緣海域絕對值明顯偏大。同時, 這些區(qū)域有著較大的水平浮力梯度, 所以更容易使qbc<0, 有利于產(chǎn)生對稱不穩(wěn)定。此外, 對于有明顯渦旋特征的區(qū)域, 冬季qbc正值主要分布在渦旋南部, 負值則多位于北部, 夏季分布情況相反。qbc普遍強于qvert, 表明在這些區(qū)域表層, 斜壓分量對于q的貢獻占主導地位。

    對不同區(qū)域進行分析, R1區(qū)域表層q在冬季以負為主, 夏季以正為主。冬季區(qū)域內(nèi)qvert同樣以負值為主, 如圖2所示, 該區(qū)域除東北角渦旋中心部分外, 負相對渦度絕對值普遍小于科氏參數(shù), 因此混合層內(nèi)密度梯度受東北季風攪拌以及黑潮入侵出現(xiàn)正值, 引發(fā)重力不穩(wěn)定可能是導致qvert<0的主要原因。與羅士浩等(2016)的結果不同, R1區(qū)域冬季q在渦旋南側邊緣存在明顯正值, 這可能與本文對于斜壓分量qbc沒有采用準地轉近似有關; 夏季由于季風的攪拌作用減弱, 重力不穩(wěn)定情況減弱,qvert多為正, 但靠近呂宋海峽海域受黑潮影響, 部分區(qū)域仍有為負的情況。而斜壓分量qbc在冬季季風對流場剪切的影響下強度也明顯高于夏季。

    圖2 研究區(qū)域10m層水平相對渦度ζ/f, 水平散度δ/f, 水平應變率S/f在1月15日(a~c)和7月15日(d~f)水平分布的瞬時圖像圖中黑色方框表示不同的研究區(qū)域; 紅色虛線表示南海諸島歸屬范圍線; 黑色虛線表示200m等深線。審圖號為GS(2016)1665Fig. 2 Horizontal distributions ofζ/f, δ/f, and S/f in the 10 m layer of the model. An instantaneous image taken on January 15 (a~c), and on July 15 (d~f)

    R2區(qū)域q及其分量的季節(jié)性特征與R1區(qū)域一致, 但由于受黑潮影響弱于R1區(qū)域,qvert<0的情況相對微弱, 特別是在夏季, 此時除了北部小部分區(qū)域外,qvert基本為正。

    類似地, R3區(qū)域q和qvert為負的情況在冬季出現(xiàn)更為頻繁, 但由于冬季風影響進一步減弱, 強度較前兩者更低。特別的是, 在夏季10°—14°N存在的強鋒面區(qū)域內(nèi),q為負的現(xiàn)象同樣十分顯著。主要是因為該區(qū)域在夏季受西南季風影響, 存在強烈的東向離岸上升流鋒面, 跨鋒面水平浮力梯度很強,同時在沿鋒面風應力的影響下引起的Ekman輸運會將冷水向暖水側跨鋒面輸運, 進一步增強水平浮力梯度并削弱垂向?qū)咏Y, 使得qbc<0, 從而克服正qvert產(chǎn)生負位渦, 引發(fā)對稱不穩(wěn)定。

    R4區(qū)域冬夏兩季q沒有明顯的正負偏向, 但冬季存在更多的q高值區(qū), 這可能是因為該海域上層流場在冬季整體強度更高, 在島嶼地形豐富的條件下容易產(chǎn)生更強的垂向流速剪切, 產(chǎn)生較大的局地位渦水平分量絕對值qbc。qvert以正值為主, 這主要歸因于該地區(qū)受季風影響較小, 風的攪拌作用較弱。但由于南沙群島的存在, 受地形尾渦等作用, 該海域全年次中尺度過程較為活躍qvert會有零星的負值分布。

    根據(jù)不同特點的典型斷面考察q及其分量的垂向分布特點, 如圖8所示。整體上, 混合層以上q較小, 正負差異并不明顯, 僅在約10m以上存在較大的絕對值, 而在混合層底及以下, 則以正值為主,表明對于整個南海, 混合層以上流體穩(wěn)定性相對更差。對于垂向分量qvert, 其絕對值在混合層內(nèi)很小,僅在表層比較顯著, 這主要是因為混合層內(nèi)的層結較弱。而對于層化較強的混合層以下區(qū)域, 由于ζ變?nèi)? 所以這部分深度范圍qvert大小主要由垂向?qū)咏Y決定, 表現(xiàn)為顯著正值, 并隨深度減小。斜壓分量qbc以負值為主, 僅在表層存在小部分正值區(qū)域, 這與準地轉平衡態(tài)下, 北半球qbc<0的特點一致?;旌蠈拥撞? 溫躍層頂, 由于等密度線分布密集, 當其發(fā)生彎曲時會產(chǎn)生較大的水平浮力梯度, 易使qbc出現(xiàn)顯著負值, 削弱流體穩(wěn)定性。與表層不同, 混合層以下,q主要貢獻來自于qvert,qbc的影響幾乎可以忽略, 表明混合層以下流體受斜壓性影響減弱, 流體的不穩(wěn)定性是主要來自于水平切變產(chǎn)生的正壓不穩(wěn)定。

    對不同典型斷面進行分析, 發(fā)現(xiàn)它們各自有著獨特的特征。S1斷面橫跨了R1區(qū)域冬季顯著的反氣旋渦, 其等密度線在渦旋位置明顯下凹。該斷面表層q在渦旋邊緣高于中心區(qū)域, 等密度線露頭處有明顯極大值存在, 表層位渦水平分量絕對值qbc有著同樣的分布特點。此外, 在混合層底部, 渦旋邊緣區(qū)域由于水平密度梯度極大,qbc有明顯負值; 位于夏季R3區(qū)域的S2斷面跨越了12°—13°N范圍內(nèi)的強上升流鋒面, 受該鋒面影響, 等密度線傾斜,水平浮力梯度增大, 部分等密度線甚至可以從表層延伸至50m以下, 因而在鋒面附近存在從表層沿等密度線延伸到混合層底的強qbc負值, 有助于削弱上層流體的穩(wěn)定性, 引發(fā)對稱不穩(wěn)定; S3斷面, 由于其所在R4區(qū)域包括了南沙群島且冬季流速相對較強,因而受地形尾渦的影響在部分區(qū)域混合層底部可能會出現(xiàn)強烈的負渦度, 導致qvert<0, 引發(fā)慣性不穩(wěn)定。

    2.3 能量轉換

    圖10 R1 (a, e)、R2 (b, f)、R3 (c, g)和R4 (d, h)各區(qū)域正壓轉換項BSK的水平分量HRS (a~d)與垂向分量VRS (e~h)隨時間變化情況圖中黑色實線表示混合層深度隨時間變化情況Fig. 10 Variations of HRS and VRS in different regions over time. Panels from top to bottom correspond to regions R1~R4,respectively. The black curve shows the variation of regionally averaged mixed-layer depth over time

    R1區(qū)域(圖9a), 秋末和冬季混合層加深的時間段內(nèi), 斜壓轉換項PKE主要存在于混合層內(nèi), 且振幅與混合層深度變化情況一致, 最大值位于混合層中間位置, 表明這一階段上層海洋PKE主要由混合層不穩(wěn)定控制, 高應變在混合層中容易產(chǎn)生鋒生作用, 激發(fā)混合層發(fā)生再分層, 使APE向SMKE轉換,這一過程主要發(fā)生在混合層中部。自2月中旬起, 混合層深度開始衰減, 混合層不穩(wěn)定產(chǎn)生的APE源和混合層內(nèi)的斜壓能量轉換開始減弱, 到了夏季與混合層以下的PKE項大致處于同一量級。R2區(qū)域的情況與R1基本一致(圖9b), 但強度較弱。R3區(qū)域(圖9c)在冬、春季與R1、R2區(qū)域一樣, 上層PKE強度與混合層深度變化情況一致。不同的是, R3區(qū)域夏季受西南季風影響, 出現(xiàn)顯著東向離岸上升流, 由于鋒面海區(qū)強水平梯度、弱層結的特點, 極易產(chǎn)生鋒面不穩(wěn)定, 引起APE向SMKE的轉換, 因此表層PKE的最大值出現(xiàn)在7—8月。至于R4區(qū)域(圖9d),PKE全年變化并不明顯, 僅在冬夏兩季稍強, 可見該區(qū)域PKE并不完全由混合層不穩(wěn)定控制。另外值得一提的是, 由于R1、R2區(qū)域冬季受到黑潮入侵的限制影響(特別是R1), PKE明顯高于其他海區(qū)。在混合層以下可能出現(xiàn)PKE為負的情況, 但相對較弱,并隨深度逐漸趨于0。

    圖9 R1 (a, e)、R2 (b, f)、R3 (c, g)和R4 (d, h)各區(qū)域平均的斜壓轉換項PKE (a~d)、正壓轉換項BSK (e~h)隨時間變化情況圖中黑線為區(qū)域平均混合層深度隨時間變化情況Fig. 9 Variations of averaged PKE and BSK in each region over time. Panels from top to bottom correspond to regions R1~R4, respectively. The black curve shows the variation of regionally averaged mixed-layer depth over time

    各區(qū)域BSK隨時間變化情況與PKE相近, 且主要存在于混合層內(nèi), 但在混合層以下也可能存在明顯正負值。相對PKE而言較弱, 尤其是對于R1 (圖9e)、R2 (圖9f)區(qū)域, 表明這兩個區(qū)域次中尺度動能的主要來源是混合層內(nèi)斜壓不穩(wěn)定引起的APE釋放, 特別是在冬季, 由于混合層較深, 斜壓性更強, 這一特征更為明顯, 該結論與Zhang等(2020)一致。R3區(qū)域(圖9g)夏季由于有強烈的離岸上升流鋒面存在,BSK與PKE的貢獻相當, 表明夏季中尺度流場剪切引起的動能轉換是該區(qū)域次中尺度動能的重要來源。正壓轉換項BSK主要由剪切拉伸等引起, 其強度與背景流場剪切密切相關, 因而在垂直方向上表現(xiàn)出與應變率相似的雙極值結構, 尤其是在R4區(qū)域, 這一特征十分明顯(圖9h)。

    不同于Zhang 等(2020)僅考慮了背景流水平剪切的作用, 本文將BSK分為水平分量(HRS)和垂向分量(VRS)(圖10)。其中HRS僅在表層20m以上有明顯正值分布, 隨深度減小甚至出現(xiàn)負值, 混合層以下開始逐漸趨于0, 且明顯的正負值往往同時出現(xiàn)。這表明, 南海大部分區(qū)域由背景流水平剪切和雷諾應力引發(fā)的背景流動能向次中尺度動能的正壓能量轉換僅發(fā)生在表面邊界層, 而在表層以下則表現(xiàn)為能量的逆向串級。

    垂向分量VRS的貢獻普遍大于水平分量HRS,表明研究區(qū)域內(nèi), BKE向SMKE的轉化更多是通過背景流垂向剪切以及垂向雷諾應力實現(xiàn)。VRS在混合層內(nèi)多為正, 混合層以下逐漸趨于0。需要強調(diào)的幾個特點是: 冬末R1區(qū)域(圖10e)在100~200m深度范圍出現(xiàn)明顯正值, 這主要是該地區(qū)北部陸架結構向外延伸以及東沙群島部分地形共同作用的結果; R3區(qū)域內(nèi)(圖10g), 與HRS一樣的, 受強烈上升流鋒面影響, VRS最大值出現(xiàn)在夏季; R4區(qū)域(圖10h)由于島嶼分布密集、地形復雜, 在部分海區(qū)會產(chǎn)生較強的背景流垂向剪切, 有可能使得等密度線密集的混合層底部、密躍層頂部區(qū)域出現(xiàn)VRS絕對值的高值區(qū), 甚至有可能大于混合層內(nèi)VRS的絕對值, 其他區(qū)域也存在這種現(xiàn)象,但相對來說比較微弱。

    次中尺度動能SMKE隨時間變化情況與各區(qū)域ζ/f均方根值一致, 但是與中尺度動能EKE的存在明顯差異(圖11)。僅R1區(qū)域由于受黑潮影響顯著,次中尺度過程活躍程度與中尺度過程活躍程度密切相關, 其他區(qū)域EKE與SMKE相關性均不顯著。背景流能量向次中尺度動能的轉換項(Total kinetic energy transform, TKE, TKE=BSK+PKE)在各區(qū)域隨時間變化情況與SMKE一致, 表明次中尺度過程的活躍程度與背景流向次中尺度動能的能量轉換密切相關, 但不同區(qū)域之間有一定的差別, R1區(qū)域一致性較好, 而R2區(qū)域8月中下旬、R3區(qū)域9月中下旬和R4區(qū)域10月上旬的SMKE圖像上均存在一個峰值, 表明此時SMKE可能受大氣強迫作用或是海洋外部能量輸送的影響更為顯著。

    圖11 不同區(qū)域上層50m區(qū)域平均的中尺度動能EKE (a)、次中尺度動能SMKE (b)與次中尺度動能的正壓轉換項BSK、斜壓轉換項PKE之和(c)隨時間變化情況Fig. 11 Variations of averaged EKE, SMKE, and the sum of BSK and PKE in the upper 50 m depth of different regions

    4 結論

    本文基于高分辨率ROMS結果, 將整個南海海域分成4個子區(qū)域, 分析了不同子區(qū)域上層200m海域次中尺度過程的季節(jié)性變化特征。研究表明, 南海次中尺度過程在不同子區(qū)域的季節(jié)性特征存在顯著差異, 與區(qū)域的位置、底部地形和受季風影響情況等條件有關。特征參量ζ/f的均方根值隨時間變化情況表明(圖3): 靠近呂宋海峽的R1區(qū)域和位于中部偏北的R2區(qū)域次中尺度過程冬季最強、夏秋季較弱, 這與Zhang等(2020)的結論一致; 靠近中南半島東部的R3區(qū)域夏秋季最強, 冬季次之, 春季最弱, 這一特征與黃小龍等(2020)的結果相吻合;位于南沙群島附近的R4區(qū)域各季節(jié)差異相對較小,僅在秋末冬初存在明顯偏小的情況。能量轉換的結果表明, TKE主要貢獻來源于斜壓不穩(wěn)定PKE (R3區(qū)域除外), 受背景流水平浮力梯度影響顯著。斜壓不穩(wěn)定帶來的勢能轉換主要出現(xiàn)在上混合層, 以混合層斜壓不穩(wěn)定為主, 強弱變化與混合層深度變化情況基本一致。特別的是, R3區(qū)域受夏季西南風影響, 存在顯著東向離岸上升流, 此時PKE的主要來源為鋒面不穩(wěn)定, 并伴隨著強烈的動能正壓轉換(BSK)。背景流場剪切和雷諾應力通過正壓不穩(wěn)定過程引起的BKE向SMKE的轉換BSK僅發(fā)生在表面邊界層, 受鋒生作用和地形導致的流場剪切(如R4區(qū)域)的影響比較顯著。

    猜你喜歡
    斜壓鋒面渦度
    熱聲耦合燃燒振蕩中火焰鋒面識別分析
    煤氣與熱力(2022年2期)2022-03-09 06:29:20
    2019年夏季長江口及鄰近海域鋒面控制下葉綠素a的分布特征及其環(huán)境影響因素分析
    海洋通報(2021年5期)2021-12-21 06:19:52
    南海相對渦度的時空變化特征分析*
    基于核心素養(yǎng)的高中地理“問題式教學”——以“鋒面氣旋”為例
    斜壓渦度的變化與臺風暴雨的關系研究
    2017年11月16日至17日博州降雪天氣分析
    科技資訊(2018年30期)2018-03-11 09:30:58
    Rossby波的線性穩(wěn)定性
    科技資訊(2016年23期)2016-05-30 20:25:49
    2005—2009年、2011年和2013年南海東北部120°E斷面秋季體積輸運的年際變化*
    海洋與湖沼(2016年1期)2016-01-15 03:50:46
    渤海夏季第一斜壓羅斯貝變形半徑的計算與分析
    一次東北南部特大暴雪的動力機制分析
    桃红色精品国产亚洲av| 国产aⅴ精品一区二区三区波| 一级毛片女人18水好多| 熟女电影av网| 美女午夜性视频免费| 久久久精品欧美日韩精品| 十八禁人妻一区二区| 国产69精品久久久久777片 | 少妇粗大呻吟视频| 亚洲精品美女久久av网站| 一边摸一边抽搐一进一小说| 国产人伦9x9x在线观看| 国产精品久久久久久久电影 | 老熟妇仑乱视频hdxx| 亚洲美女黄片视频| 欧美成人性av电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 9191精品国产免费久久| 91字幕亚洲| 久久香蕉国产精品| 午夜老司机福利片| 国产熟女xx| 亚洲一区二区三区色噜噜| 一进一出好大好爽视频| 两性夫妻黄色片| 欧美国产日韩亚洲一区| 国产精品一区二区三区四区免费观看 | 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 国产99白浆流出| 在线永久观看黄色视频| a在线观看视频网站| 欧美国产日韩亚洲一区| 在线观看www视频免费| 久久精品国产清高在天天线| 国产精华一区二区三区| 免费观看精品视频网站| 亚洲av成人一区二区三| 国产亚洲精品第一综合不卡| 午夜亚洲福利在线播放| 欧美中文综合在线视频| 在线a可以看的网站| 在线观看一区二区三区| 母亲3免费完整高清在线观看| 免费看日本二区| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| 免费看美女性在线毛片视频| 亚洲avbb在线观看| 免费无遮挡裸体视频| 搡老岳熟女国产| 日韩欧美免费精品| 欧美中文综合在线视频| 在线观看66精品国产| 国产主播在线观看一区二区| 国产精华一区二区三区| 成人精品一区二区免费| 成人三级黄色视频| 午夜亚洲福利在线播放| 性欧美人与动物交配| 在线观看舔阴道视频| 日韩欧美国产一区二区入口| 岛国视频午夜一区免费看| 午夜福利视频1000在线观看| 国产激情偷乱视频一区二区| 国产精品影院久久| 国内毛片毛片毛片毛片毛片| 中文字幕精品亚洲无线码一区| 国产激情偷乱视频一区二区| 女人高潮潮喷娇喘18禁视频| cao死你这个sao货| 久久久久精品国产欧美久久久| 日韩免费av在线播放| 少妇被粗大的猛进出69影院| 亚洲成av人片免费观看| 怎么达到女性高潮| 午夜影院日韩av| 91av网站免费观看| 色尼玛亚洲综合影院| 黄色成人免费大全| 88av欧美| 亚洲成a人片在线一区二区| 91在线观看av| 亚洲国产看品久久| 日韩高清综合在线| 久久人妻av系列| 精品国产超薄肉色丝袜足j| 一二三四社区在线视频社区8| 亚洲欧美日韩高清专用| 国产亚洲精品av在线| 村上凉子中文字幕在线| 夜夜爽天天搞| 看免费av毛片| 午夜影院日韩av| 51午夜福利影视在线观看| 欧美极品一区二区三区四区| 欧美精品亚洲一区二区| 亚洲国产欧洲综合997久久,| 亚洲人成电影免费在线| 国产伦在线观看视频一区| www日本黄色视频网| 色播亚洲综合网| 欧洲精品卡2卡3卡4卡5卡区| 国产一区二区三区视频了| 国产av不卡久久| 午夜老司机福利片| 亚洲午夜精品一区,二区,三区| 亚洲精品一卡2卡三卡4卡5卡| av国产免费在线观看| 免费搜索国产男女视频| 成人18禁高潮啪啪吃奶动态图| 一区二区三区激情视频| 久久久久久大精品| 亚洲男人的天堂狠狠| 一二三四在线观看免费中文在| 1024香蕉在线观看| 一本一本综合久久| 午夜日韩欧美国产| 欧美黑人精品巨大| 两性夫妻黄色片| 中文字幕久久专区| 午夜影院日韩av| 亚洲精品中文字幕一二三四区| 久久久久国产一级毛片高清牌| 国产精品一区二区三区四区免费观看 | 欧美日韩国产亚洲二区| 婷婷六月久久综合丁香| 午夜福利视频1000在线观看| 真人做人爱边吃奶动态| 人人妻人人澡欧美一区二区| 欧美成人一区二区免费高清观看 | 在线观看www视频免费| 国产一区二区在线av高清观看| 波多野结衣高清无吗| 一个人免费在线观看电影 | 欧美日韩精品网址| 变态另类丝袜制服| 亚洲avbb在线观看| 亚洲全国av大片| 久久中文字幕人妻熟女| 级片在线观看| 午夜免费成人在线视频| 国产精品久久久久久人妻精品电影| 欧美日韩亚洲国产一区二区在线观看| 两个人的视频大全免费| 日日干狠狠操夜夜爽| 国产片内射在线| 色哟哟哟哟哟哟| 欧美性猛交黑人性爽| 性欧美人与动物交配| 国产精品99久久99久久久不卡| 18禁国产床啪视频网站| 久久久久久亚洲精品国产蜜桃av| 777久久人妻少妇嫩草av网站| 国产精品,欧美在线| 可以在线观看的亚洲视频| 一本综合久久免费| 亚洲男人天堂网一区| 免费在线观看黄色视频的| 久久久久免费精品人妻一区二区| 国产私拍福利视频在线观看| 最新在线观看一区二区三区| av福利片在线观看| 精品久久久久久久毛片微露脸| 欧美一级毛片孕妇| 亚洲精品在线美女| 日本免费a在线| www.精华液| 老司机靠b影院| 欧美日韩亚洲综合一区二区三区_| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 亚洲av日韩精品久久久久久密| 精品欧美一区二区三区在线| 日韩欧美 国产精品| 最近视频中文字幕2019在线8| 人妻久久中文字幕网| 国产亚洲精品久久久久久毛片| 久久精品影院6| 国产区一区二久久| 国产午夜福利久久久久久| 淫秽高清视频在线观看| 黄色丝袜av网址大全| 大型黄色视频在线免费观看| 最近最新免费中文字幕在线| 一本综合久久免费| 久久精品影院6| 成年女人毛片免费观看观看9| 国产精品一区二区精品视频观看| 国产精品爽爽va在线观看网站| 精品午夜福利视频在线观看一区| 国内少妇人妻偷人精品xxx网站 | 天堂影院成人在线观看| 久久欧美精品欧美久久欧美| 男女视频在线观看网站免费 | 一夜夜www| 久久午夜综合久久蜜桃| 亚洲国产精品999在线| 免费看十八禁软件| 黄色毛片三级朝国网站| 久久这里只有精品中国| 久久久久性生活片| 国产一区二区三区视频了| 国产亚洲精品久久久久5区| 亚洲欧洲精品一区二区精品久久久| 国内少妇人妻偷人精品xxx网站 | 欧美日韩亚洲国产一区二区在线观看| 不卡一级毛片| 99久久无色码亚洲精品果冻| 中文资源天堂在线| 中文字幕人妻丝袜一区二区| 法律面前人人平等表现在哪些方面| 黄色 视频免费看| 性色av乱码一区二区三区2| 观看免费一级毛片| 午夜a级毛片| 日韩有码中文字幕| 中文亚洲av片在线观看爽| 国产精品自产拍在线观看55亚洲| 国产探花在线观看一区二区| 久久婷婷人人爽人人干人人爱| 国产精品美女特级片免费视频播放器 | 国语自产精品视频在线第100页| 国产一区二区在线av高清观看| 蜜桃久久精品国产亚洲av| 一进一出抽搐gif免费好疼| 黄频高清免费视频| 午夜日韩欧美国产| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 身体一侧抽搐| 久久精品成人免费网站| 9191精品国产免费久久| 国产午夜精品久久久久久| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 十八禁人妻一区二区| 一本综合久久免费| 欧美日韩中文字幕国产精品一区二区三区| 久久人妻av系列| 国产成人系列免费观看| aaaaa片日本免费| 一边摸一边做爽爽视频免费| 黄色视频,在线免费观看| 国产成人啪精品午夜网站| 亚洲av成人不卡在线观看播放网| 欧美成狂野欧美在线观看| 国产精品亚洲美女久久久| 免费电影在线观看免费观看| 欧美另类亚洲清纯唯美| www.自偷自拍.com| 久久午夜亚洲精品久久| 一二三四在线观看免费中文在| 极品教师在线免费播放| 成人一区二区视频在线观看| 两个人免费观看高清视频| 国产成人aa在线观看| 女同久久另类99精品国产91| 亚洲一码二码三码区别大吗| 国产午夜精品论理片| АⅤ资源中文在线天堂| xxx96com| 日韩欧美一区二区三区在线观看| 亚洲欧美精品综合久久99| 午夜福利高清视频| 久9热在线精品视频| 手机成人av网站| 露出奶头的视频| 欧美大码av| 97人妻精品一区二区三区麻豆| 在线观看免费午夜福利视频| 国产精品亚洲一级av第二区| 久久天躁狠狠躁夜夜2o2o| 欧美乱色亚洲激情| 99国产综合亚洲精品| 啦啦啦观看免费观看视频高清| 草草在线视频免费看| xxxwww97欧美| 欧美+亚洲+日韩+国产| 男插女下体视频免费在线播放| 久久久久国产精品人妻aⅴ院| a级毛片a级免费在线| 精品第一国产精品| 欧美午夜高清在线| www.熟女人妻精品国产| 男人舔女人下体高潮全视频| 少妇的丰满在线观看| 伊人久久大香线蕉亚洲五| 国产亚洲av嫩草精品影院| 亚洲最大成人中文| 免费观看人在逋| 一级毛片精品| 久久精品91无色码中文字幕| 无限看片的www在线观看| 女生性感内裤真人,穿戴方法视频| 在线播放国产精品三级| 一区福利在线观看| 国产亚洲欧美98| 亚洲精品中文字幕在线视频| 久久九九热精品免费| 老熟妇乱子伦视频在线观看| 日本a在线网址| 国产av麻豆久久久久久久| 久久中文看片网| 欧美极品一区二区三区四区| 日韩 欧美 亚洲 中文字幕| 久久久精品国产亚洲av高清涩受| 又粗又爽又猛毛片免费看| 给我免费播放毛片高清在线观看| 777久久人妻少妇嫩草av网站| 日韩欧美免费精品| 成人永久免费在线观看视频| 极品教师在线免费播放| 人人妻人人看人人澡| 十八禁网站免费在线| 热99re8久久精品国产| 两性午夜刺激爽爽歪歪视频在线观看 | 99精品久久久久人妻精品| 熟妇人妻久久中文字幕3abv| 少妇裸体淫交视频免费看高清 | 欧美黑人欧美精品刺激| 两个人看的免费小视频| 国产精品一及| 法律面前人人平等表现在哪些方面| 国产精品久久久av美女十八| 免费看日本二区| 亚洲精品一区av在线观看| 国产成人aa在线观看| 麻豆久久精品国产亚洲av| 丝袜人妻中文字幕| 全区人妻精品视频| 男插女下体视频免费在线播放| 变态另类成人亚洲欧美熟女| 丝袜人妻中文字幕| 国产精品久久久久久久电影 | 性色av乱码一区二区三区2| 亚洲成av人片在线播放无| 国产亚洲欧美在线一区二区| 人人妻人人看人人澡| 国产一区二区在线av高清观看| 50天的宝宝边吃奶边哭怎么回事| 国产欧美日韩精品亚洲av| 国产精品一区二区免费欧美| 亚洲av成人av| 黄色毛片三级朝国网站| 国产一区二区在线av高清观看| 久久久精品大字幕| 国产精品一区二区三区四区久久| 人妻丰满熟妇av一区二区三区| 一进一出好大好爽视频| 久久国产精品人妻蜜桃| 欧美黑人欧美精品刺激| 一级毛片女人18水好多| av有码第一页| 午夜日韩欧美国产| 亚洲性夜色夜夜综合| 91在线观看av| 一个人免费在线观看电影 | 一级a爱片免费观看的视频| 国产69精品久久久久777片 | 国产成人精品久久二区二区91| 一本大道久久a久久精品| 黄色视频,在线免费观看| 亚洲成av人片在线播放无| 国产精品爽爽va在线观看网站| 男男h啪啪无遮挡| 久久久国产精品麻豆| 精品高清国产在线一区| 国产一区二区三区在线臀色熟女| 香蕉久久夜色| 久久人妻福利社区极品人妻图片| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久精品大字幕| 变态另类丝袜制服| 日韩欧美在线乱码| 少妇裸体淫交视频免费看高清 | 色综合站精品国产| 亚洲av片天天在线观看| 免费观看人在逋| 熟妇人妻久久中文字幕3abv| cao死你这个sao货| 在线观看免费视频日本深夜| 久久性视频一级片| 国内少妇人妻偷人精品xxx网站 | 好看av亚洲va欧美ⅴa在| 亚洲狠狠婷婷综合久久图片| 亚洲自偷自拍图片 自拍| 亚洲av第一区精品v没综合| 99国产精品一区二区三区| 亚洲av第一区精品v没综合| 在线十欧美十亚洲十日本专区| 国产黄色小视频在线观看| 在线十欧美十亚洲十日本专区| 日日干狠狠操夜夜爽| 亚洲午夜精品一区,二区,三区| 啪啪无遮挡十八禁网站| 在线观看免费视频日本深夜| 桃红色精品国产亚洲av| 少妇粗大呻吟视频| 在线看三级毛片| 午夜福利18| 精品久久蜜臀av无| 久久这里只有精品中国| 成人18禁高潮啪啪吃奶动态图| 亚洲精品色激情综合| 国产精品一区二区免费欧美| 久久久久久国产a免费观看| 男男h啪啪无遮挡| 90打野战视频偷拍视频| 香蕉久久夜色| 亚洲黑人精品在线| 久久精品91蜜桃| 久久久国产成人免费| 黄色 视频免费看| 99精品久久久久人妻精品| 人妻夜夜爽99麻豆av| 久99久视频精品免费| 久久婷婷成人综合色麻豆| 国产亚洲av嫩草精品影院| 黄色 视频免费看| 黄色片一级片一级黄色片| 国产精品久久久久久亚洲av鲁大| 亚洲欧美精品综合久久99| 久久精品aⅴ一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 日韩欧美三级三区| 波多野结衣巨乳人妻| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 丁香六月欧美| 久久中文字幕一级| 777久久人妻少妇嫩草av网站| 18禁裸乳无遮挡免费网站照片| 久久精品91蜜桃| 欧美日韩福利视频一区二区| 日韩欧美精品v在线| 亚洲avbb在线观看| 男女床上黄色一级片免费看| 一级毛片高清免费大全| 男女之事视频高清在线观看| 黑人巨大精品欧美一区二区mp4| 91国产中文字幕| 国产精品综合久久久久久久免费| 欧美丝袜亚洲另类 | 午夜亚洲福利在线播放| 亚洲欧美精品综合一区二区三区| 俺也久久电影网| 老司机靠b影院| 中文字幕精品亚洲无线码一区| 国产野战对白在线观看| 欧美国产日韩亚洲一区| 哪里可以看免费的av片| 欧美乱码精品一区二区三区| 国产成人精品无人区| 在线a可以看的网站| 久久性视频一级片| 日韩欧美国产一区二区入口| 校园春色视频在线观看| 性欧美人与动物交配| 99re在线观看精品视频| 中文在线观看免费www的网站 | 好看av亚洲va欧美ⅴa在| 大型av网站在线播放| 中文亚洲av片在线观看爽| 在线观看日韩欧美| 无遮挡黄片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 色哟哟哟哟哟哟| 亚洲人与动物交配视频| 美女午夜性视频免费| 欧美大码av| 亚洲中文av在线| 久久精品国产综合久久久| 成人午夜高清在线视频| 不卡一级毛片| 久久久久久亚洲精品国产蜜桃av| 国产精品久久久久久精品电影| 夜夜躁狠狠躁天天躁| 亚洲狠狠婷婷综合久久图片| 欧美精品亚洲一区二区| av福利片在线| 啦啦啦免费观看视频1| 日韩欧美国产在线观看| 90打野战视频偷拍视频| 亚洲色图av天堂| 亚洲色图 男人天堂 中文字幕| 此物有八面人人有两片| 午夜成年电影在线免费观看| 亚洲精品美女久久av网站| 国产精品亚洲美女久久久| 亚洲精华国产精华精| 一进一出抽搐gif免费好疼| 女人被狂操c到高潮| 久久精品国产清高在天天线| 在线国产一区二区在线| 国产精品 国内视频| 国产一区二区三区视频了| 桃色一区二区三区在线观看| 婷婷六月久久综合丁香| 免费一级毛片在线播放高清视频| 国产成人精品无人区| 不卡一级毛片| 国产单亲对白刺激| 88av欧美| 老司机深夜福利视频在线观看| 精品久久久久久久末码| 三级国产精品欧美在线观看 | 很黄的视频免费| 五月伊人婷婷丁香| 国产成人av激情在线播放| 国产精品av久久久久免费| 亚洲专区国产一区二区| 午夜激情福利司机影院| 国产蜜桃级精品一区二区三区| 免费一级毛片在线播放高清视频| 国产成人精品无人区| 亚洲中文字幕一区二区三区有码在线看 | 日本一区二区免费在线视频| 美女午夜性视频免费| 夜夜躁狠狠躁天天躁| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲一级av第二区| 99久久精品国产亚洲精品| 男女做爰动态图高潮gif福利片| 啦啦啦观看免费观看视频高清| 午夜亚洲福利在线播放| 999久久久精品免费观看国产| 亚洲一区高清亚洲精品| 国产精品,欧美在线| 在线观看舔阴道视频| 亚洲五月天丁香| www.自偷自拍.com| 黄色丝袜av网址大全| 欧美性猛交黑人性爽| 国产精品99久久99久久久不卡| 久99久视频精品免费| 男女下面进入的视频免费午夜| 可以在线观看毛片的网站| 国产成人精品久久二区二区91| 每晚都被弄得嗷嗷叫到高潮| 精品日产1卡2卡| 国产精品一区二区三区四区久久| 日本熟妇午夜| 国产熟女午夜一区二区三区| www日本黄色视频网| 亚洲 欧美一区二区三区| 久久久久久久精品吃奶| 男人舔女人的私密视频| 成人18禁在线播放| 舔av片在线| 在线a可以看的网站| 日本撒尿小便嘘嘘汇集6| 午夜免费观看网址| 亚洲专区字幕在线| 亚洲欧美日韩高清专用| 国产精品香港三级国产av潘金莲| 麻豆一二三区av精品| 色精品久久人妻99蜜桃| 波多野结衣高清无吗| 一本精品99久久精品77| 日本免费a在线| 麻豆久久精品国产亚洲av| 成人三级做爰电影| 麻豆国产97在线/欧美 | 精品免费久久久久久久清纯| 12—13女人毛片做爰片一| 亚洲av电影不卡..在线观看| 亚洲中文av在线| 在线a可以看的网站| 国产97色在线日韩免费| 99精品欧美一区二区三区四区| 日韩国内少妇激情av| 99精品久久久久人妻精品| 欧美黑人精品巨大| 99久久99久久久精品蜜桃| 午夜免费激情av| 天天躁狠狠躁夜夜躁狠狠躁| 精品国产乱子伦一区二区三区| 欧美zozozo另类| 日韩精品青青久久久久久| 91大片在线观看| 黄频高清免费视频| 国产精品综合久久久久久久免费| 欧美国产日韩亚洲一区| 国产精品av久久久久免费| 亚洲一区二区三区色噜噜| 亚洲男人的天堂狠狠| 亚洲狠狠婷婷综合久久图片| 精品国产超薄肉色丝袜足j| 亚洲欧美精品综合久久99| 久久久久久久午夜电影| 午夜两性在线视频| 久久久久国产一级毛片高清牌| 国产三级黄色录像| 欧美成人午夜精品| 成人国语在线视频| 亚洲精品美女久久久久99蜜臀| 欧美日韩精品网址| 夜夜夜夜夜久久久久| 99精品欧美一区二区三区四区| 激情在线观看视频在线高清| 精品一区二区三区四区五区乱码| 欧美黄色片欧美黄色片| 亚洲一区中文字幕在线| 免费在线观看亚洲国产| 亚洲国产欧洲综合997久久,| 欧美精品亚洲一区二区| 欧美黄色淫秽网站| 国产午夜精品久久久久久| 精品少妇一区二区三区视频日本电影| 又爽又黄无遮挡网站| 亚洲国产欧美人成| 国产视频一区二区在线看| 亚洲欧美精品综合久久99| 一个人免费在线观看的高清视频|