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

    黑潮延伸體海域典型渦旋的次中尺度特征分析*

    2021-12-04 00:49:40張旭經(jīng)志友鄭瑞璽黃小龍曹海錦
    熱帶海洋學(xué)報(bào) 2021年6期
    關(guān)鍵詞:黑潮鋒面中尺度

    張旭, 經(jīng)志友, 鄭瑞璽, 黃小龍, 曹海錦

    1. 熱帶海洋環(huán)境國家重點(diǎn)實(shí)驗(yàn)室(中國科學(xué)院南海海洋研究所), 廣東 廣州 510301;

    2. 中國科學(xué)院大學(xué), 北京 100049;

    3. 河海大學(xué)海洋學(xué)院, 江蘇 南京 210098

    黑潮是全球較強(qiáng)的西邊界流之一, 其在吐噶喇海峽處(30°N)與親潮匯合轉(zhuǎn)向東流動, 形成黑潮延伸體(140°—180°E, 30°—40°N, 如圖1 所示)(Ma et al, 2015; Sasaki et al, 2017; 胡冬 等, 2018)。在強(qiáng)勁的東向射流和大尺度不穩(wěn)定作用下, 中尺度過程在黑潮延伸體海域較為活躍(Ma et al, 2015; Liu et al 2017; Sasaki et al, 2017; Ji et al, 2018; Jing et al,2019)?;跍u旋識別的方法, 胡冬 等(2018)研究發(fā)現(xiàn), 黑潮延伸體海域在20 年間形成了近6000 個中尺度渦旋。前人研究結(jié)果表明, 中尺度渦旋對上層海洋顆粒有機(jī)碳、浮游生物、位渦等水平輸運(yùn)與垂向沉降都有重要影響(Zhang et al, 2014; Hosegood et al, 2017; Abernathey et al, 2018; Jing et al, 2021)。不同極性的中尺度渦旋對?!?dú)鉄嵬俊⒑!獨(dú)膺吔鐚铀粨Q等也有明顯的調(diào)制作用(Jiang et al, 2019;Wang et al, 2019; Zhang et al, 2020a)。然而, 在傳統(tǒng)的地轉(zhuǎn)理論框架下, 中尺度渦旋海域地轉(zhuǎn)流貢獻(xiàn)的垂向速度較小, 其動能也趨于逆向串級(Scott et al,2005), 難以很好地解釋渦旋場物質(zhì)垂向收支平衡和能量正向串級等問題(Chelton et al, 2011; Lévy et al,2012; Abernathey et al, 2018)。近年來, 高分辨率觀測和理論研究表明, 介于中尺度和小尺度之間的次中尺度過程可能是中尺度能量正向串級和物質(zhì)能量垂向輸運(yùn)的主要途徑之一(Thomas et al, 2013;Brannigan, 2016; Mahadevan, 2016)。

    圖1 黑潮延伸體海域(虛線框)2009 年5 月2 日海表面渦動能與地轉(zhuǎn)流(箭頭)的水平分布Fig. 1 Eddy kinetic energy (shading) and geostrophic velocity (vector) of the sea surface in the Kuroshio Extension (black dashed box), based on the AVISO data

    次中尺度過程具有相對較小的時空尺度[O(1d),O(10km)]、O(1)的羅斯貝數(shù)(Ro)和瑞查德森數(shù)(Ri),是同時具有地轉(zhuǎn)和非地轉(zhuǎn)特征的三維運(yùn)動(Capet et al, 2008; Thomas et al, 2008; Mcwilliams, 2016)。理論研究和現(xiàn)場觀測表明, 次中尺度過程及其不穩(wěn)定能在上層海洋引發(fā)非地轉(zhuǎn)次級環(huán)流并形成較強(qiáng)的垂向速度, 從而顯著增強(qiáng)上層海洋物質(zhì)與能量的垂向交換(Lévy et al, 2001; Mcgillicuddy et al, 2003; Klein et al, 2009; D'Asaro et al, 2011), 對上層海洋熱和碳的垂向通量的貢獻(xiàn)超過50%(Omand et al, 2015); 此外, 次中尺度過程還具有再層化(restratification)作用, 可以調(diào)節(jié)周期性變化的混合層深度(Rocha et al,2016; Yu et al, 2019)。另一方面, 通過多種不穩(wěn)定的次中尺度過程, 能有效釋放中尺度地轉(zhuǎn)動能和儲存在鋒面內(nèi)的有效位能, 并將能量正向串級至小尺度耗散, 對中尺度能量的正向串級和耗散有著重要的意義(Capet et al, 2008; Thomas et al, 2013; Brannigan et al, 2017; Qiu et al, 2019)。

    黑潮延伸體海域次中尺度過程相關(guān)研究結(jié)果表明, 寬約1km 的次中尺度鋒面將黑潮延伸體海域?!?dú)膺吔鐚觾?nèi)的能量耗散率提高了1 到2 個數(shù)量級(D’Asar et al, 2011); 該海域活躍的次中尺度過程對上層海洋再層化有著重要的貢獻(xiàn), 且是導(dǎo)致中尺度渦旋季節(jié)性變化的重要原因之一(Sasaki et al, 2017)。雖然次中尺度過程的重要性逐漸被認(rèn)知, 但由于缺少足夠的高時空分辨率觀測或模擬數(shù)據(jù), 關(guān)于黑潮延伸體海域中尺度渦旋場的次中尺度特征研究仍不充分,其三維結(jié)構(gòu)以及在中尺度渦旋生命周期內(nèi)的演化過程仍有待于深入研究。本文將基于高分辨率衛(wèi)星觀測資料和ROMS(Regional Ocean Modeling System)模式數(shù)據(jù), 對黑潮延伸體典型反氣旋渦旋海域次中尺度過程的時空分布及其動力學(xué)特征進(jìn)行分析。

    1 數(shù)據(jù)來源與計(jì)算方法

    1.1 衛(wèi)星觀測資料

    本文使用的衛(wèi)星遙感資料包括海表葉綠素濃度、海表溫度(sea surface temperature, SST)和海表高度異常(sea level anomaly, SLA)數(shù)據(jù)。海表溫度數(shù)據(jù)采用了美國國家海洋數(shù)據(jù)中心發(fā)布的高分辨率融合產(chǎn)品GHRSST(Group High Resolution Sea Surface Temperature, http://data.nodc.noaa.gov/ghrsst/L4/GL OB/UKMO/OSTIA/), 空間分辨率為0.05°×0.05°, 時間分辨率為1day, 反演精度約為0.57℃。日平均的SLA 數(shù)據(jù)來源于由法國國家空間研究中心提供的AVISO(Archiving Validation and Interpolation of Satellite Oceanographic)網(wǎng)格化產(chǎn)品(ftp://ftp.aviso.oceanobs.com/global/), 其空間分辨率為0.25°×0.25°,反演精度為0.02m。海表葉綠素濃度數(shù)據(jù)來自中等分辨率成像光譜儀(MEdium Resolution Imaging Spectrometer, MERIS)提供的 2 級產(chǎn)品數(shù)據(jù)(ftp://merisfrs-fts-ds.eo.esa.int), 其水平分辨率最高約為300m, 絕對精確度為10-3mg·m-3。

    1.2 數(shù)值模擬資料

    本文利用ROMS 數(shù)值模式對黑潮延伸體海域進(jìn)行了高分辨率嵌套模擬。其中, 最外層模型覆蓋區(qū)域?yàn)檎麄€西太平洋(R0, 水平分辨率7.5km), 第一層嵌套模型涵蓋黑潮延伸體海域(R1, 水平分辨率1.5km), 第二層嵌套模型聚焦在黑潮延伸體主軸附近(R2, 水分辨率500m, 圖2b); 模型垂向分層均為60 層, 并在上層和近底層進(jìn)行了加密處理。初始邊界條件以及風(fēng)場強(qiáng)迫分別采用SODA(simple ocean data assimilation data)海洋數(shù)據(jù)集提供的氣候態(tài)月平均(1990—2010 年)(Carton et al, 2000)和氣候態(tài)日平均(Risien et al, 2008)QuikSCAT 風(fēng)場的數(shù)據(jù)。熱通量和淡水通量等氣候態(tài)月平均?!?dú)馔縼碓从趪H海洋大氣綜合數(shù)據(jù)集(International Comprehensive Ocean-Atmosphere Data Set, ICOADS), 空間分辨率為1°×1°。地形數(shù)據(jù)采用美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)提供的ETOPO2 數(shù)據(jù)。湍流混合則使用了KPP 參數(shù)化方案(K-profile parameterization)進(jìn)行參數(shù)化(Large et al, 1994)。最外層模型R0 計(jì)算20 年后達(dá)到穩(wěn)定, 而后在線嵌套R1 和R2 模型繼續(xù)運(yùn)行并輸出第21 年的日平均模擬結(jié)果。本文所用的模擬數(shù)據(jù)皆來自500m 分辨率的R2 模型, 第二層嵌套的R2 模型由第一層嵌套的R1 模型提供每天一次的側(cè)邊界強(qiáng)迫, 表面強(qiáng)迫包括風(fēng)場、熱通量和淡水通量,底邊界受底摩擦約束。

    為驗(yàn)證ROMS 數(shù)值模擬結(jié)果, 本文對比了模擬結(jié)果與衛(wèi)星觀測的黑潮延伸體海域海表面溫度場、流場以及海表高度異常場(圖2)。對比結(jié)果顯示, 500m 分辨率的R2 模擬結(jié)果對于黑潮主軸的刻畫較為準(zhǔn)確, 主軸位置與主軸南北向彎曲的幅度基本一致, 且SST、SLA 模態(tài)分布也較為一致。這表明R2 模擬結(jié)果能較好地刻畫黑潮延伸體鋒面和中尺度渦旋, 且相比于~25km分辨率的AVISO衛(wèi)星觀測結(jié)果, 500m分辨率的模擬數(shù)據(jù)能進(jìn)一步刻畫次中尺度鋒面等過程(黃小龍等, 2020; Luo et al, 2020; 周霄雯 等, 2020; Jing et al,2021)。因此, 本文選取R2 模擬結(jié)果對黑潮延伸體典型中尺度渦旋海域的次中尺度特征進(jìn)行分析。

    圖2 衛(wèi)星觀測的5 月3 日黑潮延伸體氣候態(tài)SST 和水平流速(箭頭)(a), 以及SLA(c)和模式模擬的第21 年5 月3 日黑潮延伸體SST、水平流速(b)及SLA(d)分布圖c、d 中黑色曲線為氣候態(tài)平均海面高度等值線, 用以表示黑潮延伸體主軸的大致位置Fig. 2 Spatial distributions of SSH (shading) and surface currents (vector), and SLA of the Kuroshio Extension provided by remote sensing satellite (a, c) and ROMS mode (b, d). The black contours in (c) and (d) are SSH, which indicate the location of the Kuroshio Extension axis.

    1.3 相關(guān)參數(shù)計(jì)算

    下列公式中, “′”均表示對數(shù)據(jù)進(jìn)行10km 高通濾波后得到的次中尺度空間異常, 如u′、v′、w′分別為對水平流速東西分量u、水平流速南北分量v以及垂向流速w進(jìn)行10km 高通濾波后得到的次中尺度流速異常。衛(wèi)星高度計(jì)資料與模擬結(jié)果顯示,該渦旋海域特征流速約為1m·s-1, 當(dāng)羅斯貝數(shù)值為O(1)時, 表明該海域次中尺度過程活躍(Thomas et al, 2008; Yang et al, 2017; Jing et al, 2021)。因此根據(jù)尺度分析法估算, 10km 約為該渦旋海域次中尺度過程的特征尺度(L=U fRo , 特征流速U=1m·s-1,Ro=1,f=10-4s-1, 計(jì)算得L=10km)(Boccaletti et al,2007)。

    為分析渦旋生命周期內(nèi)次中尺度過程的演變特征, 本文計(jì)算了該渦旋海域的海表次中尺度動能(submesoscale kinetic energy, SKE)(Qiu et al, 2014)。

    其中,u′、v′分別為水平流速東西分量、南北分量的次中尺度流速異常。

    表征流場拉伸切變幅度的水平剪切速率(horizontal strain rate,St)與表征流場切變和相對渦度之間相對強(qiáng)度的Okubo—Weiss 參數(shù)(OW)分別定義為:

    其中,ζ=v x-uy為垂向相對渦度, 羅斯貝數(shù)oR=ζf,u、v分別為水平全流的東西分量和南北分量, 下標(biāo)表示對該方向的偏導(dǎo)數(shù)。渦旋場背景OW參數(shù)為O0W =0.2wσ, 其中wσ為渦旋場OW 參數(shù)的標(biāo)準(zhǔn)差(Bracco et al, 2000)。當(dāng)OW>OW0時, 表示該區(qū)域以拉伸切變?yōu)橹鲗?dǎo)動力過程, 強(qiáng)流場剪切會進(jìn)一步增強(qiáng)水平浮力梯度, 有利于鋒面、渦絲等結(jié)構(gòu)的形成; 當(dāng)OW<OW0時, 強(qiáng)渦旋流動主導(dǎo)該區(qū)域(Bracco et al, 2000; Pasquero et al, 2001;Isern-Fontanet et al, 2004)。

    為了分析鋒生作用對渦旋場次中尺度過程的影響, 本文還計(jì)算了鋒面強(qiáng)度(frontal sharpness,Fs)與鋒生函數(shù)(front function,Ff), 二者分別表征了流場的水平浮力梯度與鋒面強(qiáng)度的變化速率(Sullivan et al, 2018; Zhang et al, 2019):

    ?hb為水平浮力梯度。

    浮力通量(buoyancy flux, Bflux)表示垂向上的湍流浮力輸運(yùn), 當(dāng)Bflux 為正數(shù)時, 非地轉(zhuǎn)擾動使浮力向上輸運(yùn), 反之則使浮力向下輸運(yùn)。Fb定義為:

    其中,ω′為垂向次中尺度流速異常,b′為次中尺度浮力異常。

    2 結(jié)果分析

    2009 年5 月2 日, Meris 水色遙感衛(wèi)星在黑潮延伸體主軸北側(cè)捕捉到了一個半徑約為100km 的中尺度反氣旋渦(圖3), 在該中尺度渦旋海域分布著側(cè)向?qū)挾燃s為O(10km)的葉綠素濃度高值螺旋帶, 與其他文獻(xiàn)揭示的開闊海域水色遙感觀測結(jié)果類似(Eldevik et al, 2002; McWilliams, 2016; Munk et al,2020), 表明次中尺度過程引起的垂向次級環(huán)流增強(qiáng)可能是貢獻(xiàn)渦旋海域葉綠素濃度增長的重要貢獻(xiàn)者之一(Lévy et al, 2001; D'Asaro et al, 2011; Gula et al,2014)。Zhang 等(2020b)通過分析全球范圍中尺度渦旋海域葉綠素濃度高值帶的時空分布特征, 發(fā)現(xiàn)葉綠素螺旋帶是渦旋場典型的次中尺度結(jié)構(gòu)之一。

    圖3 2009 年5 月2 日MERIS 海表葉綠素濃度與AVISO高度計(jì)觀測的地轉(zhuǎn)流速(黑色箭頭)的水平分布Fig. 3 Spatial distributions of chlorophyll a concentration(shading) and geostrophic velocity (vector), based on Meris and AVISO data

    由于衛(wèi)星遙感資料僅能反映海表信息, 本文選用了高分辨率的R2 模擬結(jié)果, 并選取了與遙感觀測相對應(yīng)海域內(nèi)(151°—154°30′E, 36°30′—38°30′N)的渦旋, 以分析該渦旋海域次中尺度過程的演變特征及其動力學(xué)特征。值得一提的是, 本文僅就該海域(151°—154°30′E, 36°30′—38°30'N) 的渦旋個例進(jìn)行了分析, 但類似的中尺度渦旋在黑潮延伸體海域非常普遍, 因此本文將所分析的渦旋個例稱為“典型渦旋”, 以便區(qū)分該渦旋個例與黑潮延伸體內(nèi)普遍存在的其他渦旋。

    2.1 渦旋及次中尺度過程演變特征

    為分析渦旋演變不同階段次中尺度過程的變化情況, 本文基于Yang 等(2017)及Chu 等(2014)中利用流速及SLA 最大(最小)閉合等值線界定中尺度渦旋邊界的方法, 以渦旋場SKE 最大時刻的最小SLA閉合等值線(SLA=0.14m)是否閉合為依據(jù), 將渦旋演變過程大致分為了前期、中期、后期三個階段(圖4)。并以流速最大閉合等值線為渦旋外邊界, SLA=0.28m 等值線為區(qū)分渦旋邊緣與中心的分界線, 其可以將渦旋近似分為面積相等的內(nèi)外兩部分。其中,Day1—7(Day1 為模擬第21 年4 月24 日)為渦旋演變前期, Day8—29 為渦旋演變中期, Day30—37 為渦旋演變后期。

    圖4 渦旋發(fā)展前期(a—c)、中期(d—f)和后期(g—i)表層SKE 的水平分布黑色線為渦旋外邊界, 紅色線為渦旋中心外邊界Fig. 4 Surface submesoscale kinetic energy distributions in early (a-c), mid (d-f), and late (g-i) development stages of the eddy. Eddy boundary and eddy core are marked by black and red curves in each panel, respectively

    渦旋區(qū)域海表平均SKE 以及鋒生函數(shù)的時間演變曲線 (圖5)可以更直觀地展現(xiàn)渦旋邊緣與渦旋中心SKE 在不同階段的演變情況。結(jié)果顯示, 渦旋演變前期(圖4a—c), 渦旋邊緣SKE 隨背景場鋒生作用在迅速增強(qiáng), 意味著該階段背景流場對渦旋發(fā)展起促進(jìn)作用。在渦旋演變中期, 渦旋邊界閉合, 其結(jié)構(gòu)較為穩(wěn)定, 該階段渦旋邊緣SKE 與背景場鋒生函數(shù)演變趨勢相似且二者峰值對應(yīng)良好, 其相關(guān)系數(shù)為0.86, 表明背景場鋒生過程與渦旋邊緣次中尺度動能高度相關(guān), 可能是該海域次中尺度過程生成的重要機(jī)制之一(Nikurashin et al, 2013; McWilliams,2016); 該階段渦旋中心鋒生作用始終趨向于0, 背景流場對渦旋中心影響極弱, 渦旋中心SKE 仍呈現(xiàn)下降的趨勢, 渦旋中心與渦旋邊緣平均SKE 之比由0.45 逐漸下降至 0.08, 與渦旋羅斯貝波(Vortex Rossby Waves, VRWs)理論所預(yù)測的結(jié)果相似(Koszalka et al, 2009; Rodríguez-Marroyo et al, 2009;Zhang et al, 2020b)。當(dāng)渦旋演變至后期, 渦旋無法繼續(xù)維持穩(wěn)定閉合結(jié)構(gòu)(圖 4g—i), 渦旋中心鋒生函數(shù)短暫增強(qiáng)后再次迅速減弱, 此階段渦旋迅速消亡并匯入背景場平流之中, SKE 隨渦旋消亡而逐漸消失。

    結(jié)合葉綠素濃度的水平分布圖(圖3)與SKE 演變曲線(圖5), 推測圖3 所示渦旋剛發(fā)展至中期, 此時渦旋邊緣鋒面剛剛閉合, 渦旋中心葉綠素螺旋帶結(jié)構(gòu)(次中尺度過程)由渦旋中心水體受前期背景場鋒生作用影響形成, 且會隨著渦旋演變逐漸向渦旋邊緣移動并減弱。為探討該典型渦旋海域的次中尺度特征, 本文選取了與觀測渦旋相同階段SKE 最大時刻(Day9)的模擬結(jié)果做進(jìn)一步分析。

    圖5 渦旋邊緣(藍(lán)線, 單位: 10-3m2·s-2)與渦旋中心(紅線,單位: 10-3m2·s-2)的平均次中尺度動能及渦旋邊緣(綠線,單位: 10-16s-5)與渦旋中心(紫線, 單位: 10-16s-5)鋒生函數(shù)的時間演變曲線Fig. 5 Time series of SKE in eddy periphery (blue curve),eddy center (red curve), and front function in eddy periphery(green curve) and eddy center (purple curve)

    2.2 次中尺度空間結(jié)構(gòu)特征

    模擬的海表溫度(圖6a)結(jié)果顯示, 該反氣旋渦的基本結(jié)構(gòu)與觀測相似, 在中尺度渦旋海域中存在著寬約10km 的次中尺度鋒面結(jié)構(gòu), 且其在渦旋邊緣分布更為密集, 渦內(nèi)溫度最高可達(dá)20 ℃, 較渦旋外高2~4℃。地轉(zhuǎn)流在渦旋邊緣達(dá)到1m·s-1, 而在渦旋中心較弱, 不均勻的流速分布會顯著改變該海域的相對渦度, 為次中尺度過程的發(fā)生提供有利的動力學(xué)條件。渦旋海域Ro 水平空間分布結(jié)果顯示, 由于渦旋邊緣較強(qiáng)的拉伸作用, 其正值在渦旋邊緣接近甚至大于1 且呈現(xiàn)10~20km 寬的帶狀分布, 表明次中尺度過程在渦旋邊緣具有渦絲狀結(jié)構(gòu)特征; 而渦旋中心Ro 整體為負(fù), 其最小值為-0.8, 正值主要分布在渦旋中心附近的狹窄條帶上, 表明渦旋中心也存在次中尺度過程, 但遠(yuǎn)弱于渦旋邊緣, 這一結(jié)果與SST 的水平空間分布相匹配。

    圖6 模擬結(jié)果第21 年5 月3 日渦旋海域的海表面溫度、水平流速(箭頭)(a)與羅斯貝數(shù)Ro(b)的水平分布a 中白色線為等溫線; b 中黑色實(shí)線為下文所分析的37°30'N 斷面位置Fig. 6 Maps of SST (a) and Ro (b) from the R2 simulation. The vectors and white contours in (a) are for surface currents and isotherms, respectively. The black line in (b) is the location of the 37°30'N section analyzed later in the paper

    跨渦旋斷面(37°30′N 斷面, 圖6 中黑線所示位置)進(jìn)一步展示了該中尺度渦旋的三維結(jié)構(gòu)。溫度和位勢密度斷面結(jié)果(圖7)顯示, 渦旋中心溫度最高達(dá)到20℃, 而密度最低約為1024.5kg·m-3, 且高溫、低密水體主要存在于100m 以淺。渦旋邊緣露頭的等溫線(等密線)向下彎曲深度可達(dá)300m 以下, 較為密集的等溫線(等密線)分布表明渦旋邊緣存在較強(qiáng)的鋒面, 其水平溫度梯度(dT/dx)與密度梯度在渦旋邊緣最大分別可達(dá)到 0.35 ℃·km-1和1.2×10-4kg·m-4,與D’Asaro 等(2011)在黑潮延伸體海域所觀測的溫度與密度鋒面量值相近。此外, 由于水體混合不均勻, 渦旋中心同樣存在類似渦旋邊緣的等溫線(等密線)露頭現(xiàn)象, 其主要分布在渦旋中心上100m, 渦旋中心溫度與密度水平梯度最大值分別為6.1×10-2℃·km-1和2.3×10-5kg·m-4。以上分析結(jié)果表明, 整個中尺度渦旋海域內(nèi)均存在強(qiáng)度不同的次中尺度鋒面結(jié)構(gòu), 渦旋邊緣鋒面比中心約強(qiáng)5 倍。

    圖7 37°30'N 斷面的溫度(a)和位勢密度(b)分布a、b 中灰色線分別為等溫線和等密度線Fig. 7 Vertical profiles of temperature (a) and potential density (b). Gray contours in (a) and (b) are isotherms and isopycnals,respectively

    2.3 次中尺度動力學(xué)特征

    為探討2.2 節(jié)渦旋海域次中尺度結(jié)構(gòu)空間特征可能的形成原因, 本文進(jìn)一步診斷計(jì)算了與次中尺度過程相關(guān)的動力學(xué)參量。渦旋海域OW 參數(shù)標(biāo)準(zhǔn)差為σw=1.2×10-9s-2, 即背景場OW 參數(shù)OW0=2.4× 10-10s-2。渦旋海域St、OW 參數(shù)水平空間分布(圖8b) 顯示, 渦旋邊緣存在明顯的流場剪切和正 OW 參數(shù)(均值為5.0×10-10s-2, 大于OW0), 表明渦旋邊緣流動以剪切變形為主導(dǎo), 有利于鋒面的加強(qiáng), 易引起鋒生過程。Ro、St、OW 參數(shù)的診斷結(jié)果表明, 在該渦旋邊緣具有次中尺度特征的流場側(cè)向拉伸剪切有助于水平浮力梯度的增強(qiáng)和鋒生過程的發(fā)生。而渦旋中心St較弱,且OW 參數(shù)以負(fù)值為主(負(fù)值均值為-3.2×10-10s-2, 小于-OW0), 表明渦旋中心以旋轉(zhuǎn)流動為主動。渦旋中心存在少量渦絲狀OW 正值結(jié)構(gòu)表明渦旋中心存在較弱的次中尺度過程, 但背景流場剪切較弱, 不利于次中尺度過程發(fā)展。

    圖8 模擬結(jié)果第21 年5 月3 日渦旋海域的流場水平剪切率(a)與Okubo—Weiss 參數(shù)(b)的水平分布Fig. 8 Spatial distributions of the strain rate (a) and Okubo-Weiss parameter (b)

    Fs與Ff的斷面分布結(jié)果(圖9a、b)顯示, 渦旋邊緣區(qū)域存在較強(qiáng)的鋒面, 鋒面強(qiáng)度最大值為9.6×10-14s-4。渦旋邊緣區(qū)域鋒生趨勢為正(圖9b), 表明背景變形流場會擠壓鋒面, 使鋒面區(qū)域水平浮力梯度增強(qiáng), 為次中尺度不穩(wěn)定的發(fā)生提供了有利條件。而渦旋中心鋒面較弱, 鋒面強(qiáng)度最大值僅為1.1×10-14s-4,比渦旋邊緣約小一個量級(圖9a), 且鋒生作用在該區(qū)域整體較弱, 其均值趨向于0(圖5、圖9b), 意味著渦旋中心背景流場對次中尺度過程的形成與發(fā)展無明顯地促進(jìn)作用。

    渦旋邊緣沿鋒面區(qū)域Ro 值接近于1, 且其分布與鋒面強(qiáng)度和鋒生函數(shù)對應(yīng)良好, 進(jìn)一步表明受鋒生作用的影響, 渦旋邊緣形成了較強(qiáng)的次中尺度過程(圖9c)。同時, 渦旋邊緣存在明顯的垂向次級環(huán)流, 其在渦旋外側(cè)下降, 在渦旋內(nèi)側(cè)上升,顯著增強(qiáng)了鋒面區(qū)域的垂向流速, 最大可達(dá)100m·d-1(圖9a), 比中尺度過程引發(fā)的垂向流速強(qiáng)約一個量級(McWilliams, 2016; Sullivan et al,2018)。垂向次中尺度流速異常與浮力通量的斷面分布(圖9a、d)顯示, 次中尺度過程能夠引起較大的垂向流速, 尤其在渦旋邊緣, 可以影響至混合層以下, 達(dá)到200m 深度處, 顯著增強(qiáng)了該區(qū)域物質(zhì)、能量等的垂向輸送, 從而有效促進(jìn)了海洋表層與內(nèi)部的物質(zhì)交換。

    圖9 鋒面強(qiáng)度(填色)及垂向次中尺度流速異常ω’(箭頭)(a)、鋒生函數(shù)(b)、羅斯貝數(shù)(c)和垂向浮力通量(d)的37°30'N斷面分布a 中灰色曲線為混合層深度(采用密度閾值法計(jì)算, σ=0.03kg·m-3); a—d 中黑色豎直線為渦旋邊緣與中心的分界Fig. 9 Vertical profiles of frontal sharpness (shading), vertical velocity anomaly (black arrows) (a), front function (b), Ro (c),and bouncy flux (d). The grey curve in (a) represents the mixed-layer depth, and the black lines in each panel separate the eddy center and periphery

    3 總結(jié)與討論

    本文利用高分辨率ROMS 模式模擬結(jié)果結(jié)合高分辨率衛(wèi)星觀測資料, 對黑潮延伸體典型中尺度渦旋海域的次中尺度特征進(jìn)行了分析。衛(wèi)星觀測和高分辨率模式模擬結(jié)果均顯示, 次中尺度過程在渦旋邊緣具有顯著的渦絲狀結(jié)構(gòu); 進(jìn)一步的動力學(xué)參數(shù)診斷分析結(jié)果顯示, 背景流場剪切引起的鋒生過程可能是渦旋邊緣次中尺度動能增強(qiáng)的重要機(jī)制。模擬分析結(jié)果顯示, 渦旋邊緣背景流場剪切、強(qiáng)水平密度梯度以及弱垂向?qū)咏Y(jié)有利于次中尺度不穩(wěn)定的發(fā)生。垂向結(jié)構(gòu)分析結(jié)果表明, 次中尺度過程可以引起跨鋒面的次級環(huán)流并顯著增強(qiáng)其所在區(qū)域的垂向流速, 最大可達(dá)100m·d-1, 尤其在渦旋邊緣, 能影響到數(shù)百米深度處, 可以有效地促進(jìn)海洋表層與內(nèi)部的物質(zhì)交換, 從而顯著增強(qiáng)黑潮延伸體上層海洋物質(zhì)、能量等的垂向輸運(yùn)。

    本文僅基于衛(wèi)星觀測與數(shù)值模擬結(jié)果探討了黑潮延伸體典型渦旋海域次中尺度過程的動力學(xué)特征并初步解釋了其可能的形成原因, 關(guān)于次中尺度過程的形成與演變機(jī)制, 以及多尺度過程的能量交換還有待于結(jié)合現(xiàn)場觀測開展更深入的研究。

    猜你喜歡
    黑潮鋒面中尺度
    譜松弛動力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    熱聲耦合燃燒振蕩中火焰鋒面識別分析
    煤氣與熱力(2022年2期)2022-03-09 06:29:20
    2019年夏季長江口及鄰近海域鋒面控制下葉綠素a的分布特征及其環(huán)境影響因素分析
    南海中尺度渦的形轉(zhuǎn)、內(nèi)轉(zhuǎn)及平移運(yùn)動研究
    基于深度學(xué)習(xí)的中尺度渦檢測技術(shù)及其在聲場中的應(yīng)用
    基于核心素養(yǎng)的高中地理“問題式教學(xué)”——以“鋒面氣旋”為例
    2016年7月四川持續(xù)性強(qiáng)降水的中尺度濾波分析
    黃淮地區(qū)一次暖區(qū)大暴雨的中尺度特征分析
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計(jì)數(shù)據(jù)的黑潮變異特征
    国产精品99久久久久久久久| 午夜免费激情av| 少妇人妻精品综合一区二区 | 欧美人与善性xxx| 色哟哟哟哟哟哟| 亚洲欧洲国产日韩| 免费大片18禁| 人妻系列 视频| 亚洲人成网站高清观看| 亚洲电影在线观看av| 少妇被粗大猛烈的视频| 综合色丁香网| 成年免费大片在线观看| 内射极品少妇av片p| 22中文网久久字幕| 精品久久国产蜜桃| 乱系列少妇在线播放| 中文字幕人妻熟人妻熟丝袜美| 两性午夜刺激爽爽歪歪视频在线观看| 久久99精品国语久久久| 18禁在线播放成人免费| 国产中年淑女户外野战色| 精品久久久久久久久久免费视频| 国产极品精品免费视频能看的| 亚洲欧美成人综合另类久久久 | 成人性生交大片免费视频hd| 青春草视频在线免费观看| 日本成人三级电影网站| 免费av毛片视频| 亚洲欧美日韩卡通动漫| 国产91av在线免费观看| 九草在线视频观看| 热99在线观看视频| 18禁在线无遮挡免费观看视频| 一级av片app| av在线亚洲专区| 亚洲精品亚洲一区二区| 国产伦精品一区二区三区视频9| 亚洲国产欧美在线一区| 美女国产视频在线观看| 亚洲国产欧洲综合997久久,| 九九久久精品国产亚洲av麻豆| 午夜a级毛片| 婷婷六月久久综合丁香| 久久久久久久久大av| 99久久中文字幕三级久久日本| 久久久久久久久久黄片| 国产成人福利小说| 久久久精品欧美日韩精品| 最近视频中文字幕2019在线8| 一区二区三区四区激情视频 | 白带黄色成豆腐渣| 人妻夜夜爽99麻豆av| 久久99热这里只有精品18| 一个人免费在线观看电影| 久久99热6这里只有精品| 波多野结衣巨乳人妻| 成人毛片a级毛片在线播放| 日韩三级伦理在线观看| 麻豆国产av国片精品| 又黄又爽又刺激的免费视频.| 免费av不卡在线播放| 欧美bdsm另类| 午夜老司机福利剧场| 欧洲精品卡2卡3卡4卡5卡区| 亚洲色图av天堂| 黄片wwwwww| 久久6这里有精品| 成人二区视频| 日本熟妇午夜| 精品久久久久久久久亚洲| 国产精品女同一区二区软件| 亚洲av电影不卡..在线观看| 欧美不卡视频在线免费观看| 免费电影在线观看免费观看| 久久午夜亚洲精品久久| 亚洲av中文字字幕乱码综合| 又粗又爽又猛毛片免费看| 给我免费播放毛片高清在线观看| 亚洲图色成人| 美女脱内裤让男人舔精品视频 | 亚洲成人久久爱视频| 国内精品久久久久精免费| 欧美日韩一区二区视频在线观看视频在线 | 国产女主播在线喷水免费视频网站 | 不卡视频在线观看欧美| 在线观看一区二区三区| 国产精品久久视频播放| 别揉我奶头 嗯啊视频| 亚洲国产精品合色在线| 99久久中文字幕三级久久日本| 97热精品久久久久久| 直男gayav资源| 色综合色国产| 国产精品,欧美在线| 欧美性猛交黑人性爽| 国产老妇伦熟女老妇高清| 欧美最黄视频在线播放免费| 国产av在哪里看| 成人特级av手机在线观看| 亚洲人与动物交配视频| 亚洲无线观看免费| www.色视频.com| 国产成人影院久久av| 观看免费一级毛片| av在线蜜桃| 哪里可以看免费的av片| 久久久国产成人精品二区| 久久鲁丝午夜福利片| avwww免费| 夜夜爽天天搞| 啦啦啦观看免费观看视频高清| 中文字幕av成人在线电影| 乱人视频在线观看| 免费黄网站久久成人精品| 麻豆成人av视频| 男人舔女人下体高潮全视频| av女优亚洲男人天堂| 麻豆一二三区av精品| 国产三级在线视频| 嫩草影院新地址| 男人狂女人下面高潮的视频| 国产黄片美女视频| 亚洲精品色激情综合| 日本欧美国产在线视频| 亚洲成av人片在线播放无| 亚洲国产精品sss在线观看| 一区二区三区四区激情视频 | 国国产精品蜜臀av免费| 成年av动漫网址| 日本成人三级电影网站| 久久久精品欧美日韩精品| 国产一级毛片在线| 精品一区二区三区人妻视频| 18禁在线播放成人免费| 国产 一区 欧美 日韩| 日本黄色视频三级网站网址| 乱人视频在线观看| 午夜久久久久精精品| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站| 精品无人区乱码1区二区| 日韩国内少妇激情av| 日韩欧美一区二区三区在线观看| 精品欧美国产一区二区三| 亚洲av中文字字幕乱码综合| 一区福利在线观看| 欧美色视频一区免费| 看黄色毛片网站| 亚洲国产精品sss在线观看| 亚洲婷婷狠狠爱综合网| 少妇的逼水好多| 亚洲自拍偷在线| 性插视频无遮挡在线免费观看| 十八禁国产超污无遮挡网站| 最好的美女福利视频网| 一区二区三区免费毛片| 波多野结衣巨乳人妻| 婷婷色综合大香蕉| av女优亚洲男人天堂| 欧美+日韩+精品| 成人国产麻豆网| 久久久久九九精品影院| 亚洲无线观看免费| 国产男人的电影天堂91| 久久草成人影院| 日韩国内少妇激情av| 12—13女人毛片做爰片一| 日本色播在线视频| 一区福利在线观看| 亚洲精品国产成人久久av| 午夜老司机福利剧场| 国产69精品久久久久777片| av专区在线播放| 午夜福利高清视频| 少妇的逼好多水| 好男人视频免费观看在线| 桃色一区二区三区在线观看| 十八禁国产超污无遮挡网站| www.av在线官网国产| 亚洲精品日韩在线中文字幕 | 观看免费一级毛片| 免费av不卡在线播放| www日本黄色视频网| 国产精品精品国产色婷婷| 99热6这里只有精品| 午夜精品一区二区三区免费看| 亚洲国产高清在线一区二区三| 久久中文看片网| 国产亚洲av嫩草精品影院| 欧美另类亚洲清纯唯美| 亚洲三级黄色毛片| 在线观看免费视频日本深夜| 亚洲欧美日韩无卡精品| a级毛片a级免费在线| 日本欧美国产在线视频| 一本精品99久久精品77| 国产毛片a区久久久久| 2022亚洲国产成人精品| 国产精品爽爽va在线观看网站| 亚洲精品国产av成人精品| 免费黄网站久久成人精品| 久久鲁丝午夜福利片| 久久99蜜桃精品久久| 国产一区二区激情短视频| 亚洲欧美中文字幕日韩二区| 精品少妇黑人巨大在线播放 | 欧美一区二区亚洲| 人妻系列 视频| 三级男女做爰猛烈吃奶摸视频| 九色成人免费人妻av| 精品日产1卡2卡| 寂寞人妻少妇视频99o| 免费大片18禁| 一卡2卡三卡四卡精品乱码亚洲| 老熟妇乱子伦视频在线观看| 我的女老师完整版在线观看| 少妇猛男粗大的猛烈进出视频 | 观看免费一级毛片| 国内精品一区二区在线观看| 国产午夜精品久久久久久一区二区三区| 中文字幕熟女人妻在线| 国产精品久久久久久精品电影小说 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 校园人妻丝袜中文字幕| 国产精品伦人一区二区| 亚洲,欧美,日韩| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 在线免费观看不下载黄p国产| 国产探花在线观看一区二区| 久久精品国产自在天天线| 亚洲国产色片| 国产亚洲精品久久久久久毛片| 我要看日韩黄色一级片| 黄色欧美视频在线观看| 高清日韩中文字幕在线| 特大巨黑吊av在线直播| 波多野结衣高清无吗| 亚洲精品影视一区二区三区av| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 久久精品国产99精品国产亚洲性色| 亚洲国产高清在线一区二区三| 天美传媒精品一区二区| 久久久欧美国产精品| 美女国产视频在线观看| 欧美成人免费av一区二区三区| ponron亚洲| 99国产精品一区二区蜜桃av| 欧美激情国产日韩精品一区| 久久久精品欧美日韩精品| 国产高清三级在线| 亚洲精品成人久久久久久| 国产精品一区www在线观看| 网址你懂的国产日韩在线| 日韩强制内射视频| 日韩一区二区视频免费看| www.色视频.com| 欧美丝袜亚洲另类| 12—13女人毛片做爰片一| 欧美成人a在线观看| 日日摸夜夜添夜夜添av毛片| 成年女人永久免费观看视频| 国内精品久久久久精免费| 日韩精品有码人妻一区| 特大巨黑吊av在线直播| 亚洲欧美日韩无卡精品| 青春草国产在线视频 | 两个人视频免费观看高清| 亚洲av中文字字幕乱码综合| 国产av不卡久久| 成熟少妇高潮喷水视频| 最好的美女福利视频网| 中文字幕精品亚洲无线码一区| 国语自产精品视频在线第100页| a级毛片a级免费在线| 精品久久久久久久久久免费视频| 国产精品一及| 高清毛片免费看| 日本与韩国留学比较| 在现免费观看毛片| 美女脱内裤让男人舔精品视频 | av卡一久久| 国产高潮美女av| 国产伦精品一区二区三区视频9| 午夜老司机福利剧场| 非洲黑人性xxxx精品又粗又长| 亚洲精华国产精华液的使用体验 | 最近视频中文字幕2019在线8| 蜜桃久久精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 黑人高潮一二区| 国产黄色视频一区二区在线观看 | av福利片在线观看| 国产私拍福利视频在线观看| 久久久国产成人免费| 小蜜桃在线观看免费完整版高清| 日本与韩国留学比较| 天天躁日日操中文字幕| 蜜桃亚洲精品一区二区三区| 网址你懂的国产日韩在线| 哪里可以看免费的av片| 成人鲁丝片一二三区免费| 不卡一级毛片| 校园人妻丝袜中文字幕| 99久久成人亚洲精品观看| av.在线天堂| 国产 一区 欧美 日韩| 人体艺术视频欧美日本| 国产成年人精品一区二区| 精品欧美国产一区二区三| 大又大粗又爽又黄少妇毛片口| 成人毛片a级毛片在线播放| 偷拍熟女少妇极品色| 欧美高清性xxxxhd video| 插逼视频在线观看| 国产视频首页在线观看| 91精品国产九色| 三级男女做爰猛烈吃奶摸视频| 变态另类丝袜制服| 欧美3d第一页| 变态另类成人亚洲欧美熟女| 黑人高潮一二区| 干丝袜人妻中文字幕| 免费av不卡在线播放| 久久久久久久午夜电影| 一级av片app| 高清午夜精品一区二区三区 | 欧美另类亚洲清纯唯美| 中文字幕av成人在线电影| 日韩亚洲欧美综合| 女同久久另类99精品国产91| 六月丁香七月| 国产毛片a区久久久久| 国产黄色小视频在线观看| 日日摸夜夜添夜夜添av毛片| 欧美一区二区亚洲| 内射极品少妇av片p| 深爱激情五月婷婷| 亚洲无线观看免费| 嫩草影院精品99| 免费观看精品视频网站| 少妇的逼好多水| 欧美3d第一页| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 亚洲精品亚洲一区二区| 最近的中文字幕免费完整| 22中文网久久字幕| 一级黄片播放器| 国产午夜福利久久久久久| 亚洲欧美清纯卡通| 国内精品美女久久久久久| 给我免费播放毛片高清在线观看| 身体一侧抽搐| 人妻少妇偷人精品九色| 男人舔奶头视频| 亚洲人与动物交配视频| av天堂在线播放| 国产美女午夜福利| 中文字幕av成人在线电影| 男人的好看免费观看在线视频| 国产精品一区二区三区四区免费观看| 国产精品久久久久久av不卡| 美女脱内裤让男人舔精品视频 | 国产成人一区二区在线| 嫩草影院新地址| 伊人久久精品亚洲午夜| а√天堂www在线а√下载| 内射极品少妇av片p| 国产黄色小视频在线观看| 亚洲在久久综合| 一级毛片久久久久久久久女| 国产伦在线观看视频一区| 青青草视频在线视频观看| 欧美+日韩+精品| 一个人免费在线观看电影| 欧美激情国产日韩精品一区| 亚洲aⅴ乱码一区二区在线播放| 日本熟妇午夜| 嫩草影院入口| 国产精品久久久久久久电影| 国内精品宾馆在线| 在线观看av片永久免费下载| 国产精品一二三区在线看| 村上凉子中文字幕在线| 男人狂女人下面高潮的视频| 精品国产三级普通话版| 国产乱人偷精品视频| 久久精品国产亚洲av香蕉五月| 国产精品永久免费网站| 91麻豆精品激情在线观看国产| 最近最新中文字幕大全电影3| 99热只有精品国产| 午夜精品一区二区三区免费看| 长腿黑丝高跟| 国产成人午夜福利电影在线观看| 亚洲av一区综合| 国产成人一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久久久久免费av| 一本久久精品| 91久久精品电影网| 日本撒尿小便嘘嘘汇集6| av在线观看视频网站免费| 欧美3d第一页| 3wmmmm亚洲av在线观看| 亚洲国产精品久久男人天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产午夜精品久久久久久一区二区三区| 久久久国产成人精品二区| 欧美日韩精品成人综合77777| 精品99又大又爽又粗少妇毛片| 九九久久精品国产亚洲av麻豆| 亚洲天堂国产精品一区在线| 九色成人免费人妻av| 久久久久久久久中文| 国产精品女同一区二区软件| 国产精品.久久久| 成人无遮挡网站| 免费观看在线日韩| 一进一出抽搐动态| 亚洲欧美日韩无卡精品| 别揉我奶头 嗯啊视频| 久久中文看片网| 日韩视频在线欧美| 国产精品日韩av在线免费观看| 中文资源天堂在线| 日日干狠狠操夜夜爽| 夫妻性生交免费视频一级片| 久久久久网色| 久久99蜜桃精品久久| 中文字幕制服av| 亚洲美女视频黄频| 国产高清三级在线| 亚洲精品国产成人久久av| 欧美+日韩+精品| 久久韩国三级中文字幕| 在线天堂最新版资源| 日韩亚洲欧美综合| 国产av在哪里看| 99热只有精品国产| 成人欧美大片| 亚洲av一区综合| 久久久久久久久久久免费av| 99久国产av精品国产电影| 麻豆久久精品国产亚洲av| 国产精品一区www在线观看| 2022亚洲国产成人精品| 99久国产av精品| 久久久成人免费电影| 黄色欧美视频在线观看| 99九九线精品视频在线观看视频| 久久久久久久久久久免费av| 精品欧美国产一区二区三| 亚洲精品日韩在线中文字幕 | 爱豆传媒免费全集在线观看| 午夜免费男女啪啪视频观看| 日本爱情动作片www.在线观看| 三级男女做爰猛烈吃奶摸视频| 国产精品嫩草影院av在线观看| 日韩av不卡免费在线播放| 国产精品伦人一区二区| 色吧在线观看| 久久久欧美国产精品| 狠狠狠狠99中文字幕| 人人妻人人澡人人爽人人夜夜 | 国产亚洲91精品色在线| 日日摸夜夜添夜夜添av毛片| 国产精品久久视频播放| 成人毛片a级毛片在线播放| 啦啦啦啦在线视频资源| 亚洲av.av天堂| 国产伦精品一区二区三区四那| 日本免费a在线| 一级黄色大片毛片| av女优亚洲男人天堂| 亚洲成人精品中文字幕电影| 亚洲久久久久久中文字幕| 亚洲精品色激情综合| 国产成人午夜福利电影在线观看| 99久久无色码亚洲精品果冻| 不卡一级毛片| 干丝袜人妻中文字幕| 国产成人91sexporn| 国产欧美日韩精品一区二区| а√天堂www在线а√下载| 中文亚洲av片在线观看爽| 亚洲欧美中文字幕日韩二区| 18禁黄网站禁片免费观看直播| 国产又黄又爽又无遮挡在线| 舔av片在线| 岛国毛片在线播放| 久久久久网色| www日本黄色视频网| 男女那种视频在线观看| 一级二级三级毛片免费看| 波野结衣二区三区在线| 亚洲成人久久爱视频| 3wmmmm亚洲av在线观看| 亚洲三级黄色毛片| 亚洲av第一区精品v没综合| 全区人妻精品视频| 麻豆乱淫一区二区| 丰满的人妻完整版| 在线国产一区二区在线| 久久午夜福利片| 日韩av不卡免费在线播放| 中文字幕制服av| 欧美精品国产亚洲| 色5月婷婷丁香| videossex国产| 国产女主播在线喷水免费视频网站 | 噜噜噜噜噜久久久久久91| 黄色配什么色好看| 在线免费十八禁| 免费观看a级毛片全部| 国产精品麻豆人妻色哟哟久久 | 久久精品影院6| 国产在线精品亚洲第一网站| 亚洲国产欧洲综合997久久,| 亚洲成av人片在线播放无| 国产精品爽爽va在线观看网站| 久久久欧美国产精品| 欧美在线一区亚洲| 身体一侧抽搐| 国产精品蜜桃在线观看 | 亚洲第一电影网av| 综合色av麻豆| 在线观看免费视频日本深夜| 人妻系列 视频| 久久久久久久久大av| kizo精华| 亚洲va在线va天堂va国产| 亚洲精品亚洲一区二区| 亚洲国产精品国产精品| 国产精品爽爽va在线观看网站| 欧美日韩一区二区视频在线观看视频在线 | 亚洲丝袜综合中文字幕| 69av精品久久久久久| 麻豆久久精品国产亚洲av| 如何舔出高潮| 麻豆国产97在线/欧美| 国产精品一区二区三区四区免费观看| 人妻少妇偷人精品九色| 热99在线观看视频| 在线国产一区二区在线| 一级毛片我不卡| 国产av在哪里看| 久久人妻av系列| 免费观看人在逋| 国产一级毛片七仙女欲春2| 性欧美人与动物交配| 黑人高潮一二区| 不卡视频在线观看欧美| 在线免费十八禁| 成人国产麻豆网| 看黄色毛片网站| 日韩欧美精品免费久久| 国产乱人视频| 高清在线视频一区二区三区 | 搡老妇女老女人老熟妇| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 亚洲熟妇中文字幕五十中出| 国产精品国产高清国产av| 最近2019中文字幕mv第一页| 久久久久久伊人网av| 欧美日本视频| 久久九九热精品免费| 看非洲黑人一级黄片| 国产精品野战在线观看| 最近视频中文字幕2019在线8| 伦精品一区二区三区| 亚洲国产精品成人综合色| 男的添女的下面高潮视频| 久久精品国产亚洲av天美| 一级毛片电影观看 | 人人妻人人澡人人爽人人夜夜 | 最新中文字幕久久久久| 久久精品国产鲁丝片午夜精品| 亚州av有码| 日韩国内少妇激情av| 国产精华一区二区三区| 一区二区三区高清视频在线| 日韩一区二区三区影片| 男人舔女人下体高潮全视频| 18禁在线播放成人免费| 天堂av国产一区二区熟女人妻| 日本免费a在线| 欧美日韩一区二区视频在线观看视频在线 | 我要搜黄色片| 看十八女毛片水多多多| 午夜福利视频1000在线观看| 精品人妻偷拍中文字幕| 欧美成人免费av一区二区三区| 搡女人真爽免费视频火全软件| 大又大粗又爽又黄少妇毛片口| 亚洲精品国产成人久久av| 亚洲最大成人中文| 男人的好看免费观看在线视频| 神马国产精品三级电影在线观看| 欧美成人精品欧美一级黄| av又黄又爽大尺度在线免费看 | 国产成人a区在线观看| or卡值多少钱| 夫妻性生交免费视频一级片| 久久久国产成人免费| 国产午夜精品久久久久久一区二区三区| 国产在视频线在精品| 舔av片在线| 熟女电影av网| a级毛色黄片| 国产精品久久久久久精品电影|