• <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ù)的黑潮變異特征
    久久久久久久大尺度免费视频| 在线天堂中文资源库| 久久久精品94久久精品| videos熟女内射| bbb黄色大片| 免费在线观看影片大全网站 | 欧美人与性动交α欧美精品济南到| 一二三四在线观看免费中文在| 一级片免费观看大全| 我要看黄色一级片免费的| 久久99一区二区三区| 亚洲精品在线美女| 中文字幕人妻丝袜制服| 精品久久久精品久久久| 国产熟女欧美一区二区| 婷婷色综合www| 满18在线观看网站| 脱女人内裤的视频| 大香蕉久久网| 亚洲av片天天在线观看| 最近中文字幕2019免费版| netflix在线观看网站| 国产亚洲欧美在线一区二区| 在线观看免费高清a一片| 成人亚洲精品一区在线观看| av国产久精品久网站免费入址| 成在线人永久免费视频| 夜夜骑夜夜射夜夜干| 永久免费av网站大全| 久久国产精品影院| 青草久久国产| 欧美日韩黄片免| 国产精品九九99| 精品久久久久久电影网| 少妇精品久久久久久久| 国产亚洲欧美精品永久| 久久久精品免费免费高清| 成人18禁高潮啪啪吃奶动态图| 成年人免费黄色播放视频| 在线观看免费高清a一片| 久久人妻福利社区极品人妻图片 | 少妇被粗大的猛进出69影院| 黑人欧美特级aaaaaa片| 国产一区二区激情短视频 | 婷婷丁香在线五月| 国产欧美日韩一区二区三 | 国产精品av久久久久免费| 欧美性长视频在线观看| 亚洲精品日韩在线中文字幕| 丝袜在线中文字幕| 两人在一起打扑克的视频| 国产淫语在线视频| 久久影院123| 五月开心婷婷网| 欧美人与性动交α欧美软件| 欧美 亚洲 国产 日韩一| 黄网站色视频无遮挡免费观看| 99久久综合免费| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人精品一二三区| 亚洲激情五月婷婷啪啪| 18禁国产床啪视频网站| 91九色精品人成在线观看| 最新的欧美精品一区二区| 亚洲图色成人| 在线观看国产h片| 视频区欧美日本亚洲| 国产精品一区二区精品视频观看| 日本色播在线视频| 日本欧美视频一区| 高潮久久久久久久久久久不卡| svipshipincom国产片| 午夜91福利影院| 免费看av在线观看网站| 欧美97在线视频| 99精国产麻豆久久婷婷| 99香蕉大伊视频| 国产成人啪精品午夜网站| 啦啦啦中文免费视频观看日本| 秋霞在线观看毛片| 日韩视频在线欧美| 真人做人爱边吃奶动态| 色播在线永久视频| 黄色 视频免费看| 美女午夜性视频免费| 欧美日韩综合久久久久久| 亚洲人成电影观看| bbb黄色大片| 久久久国产一区二区| 国产免费视频播放在线视频| 狂野欧美激情性bbbbbb| 天天躁狠狠躁夜夜躁狠狠躁| 黄色一级大片看看| 日韩av在线免费看完整版不卡| 桃花免费在线播放| 亚洲,一卡二卡三卡| 亚洲一区中文字幕在线| 韩国精品一区二区三区| 久久精品国产亚洲av涩爱| 国产色视频综合| 中文字幕最新亚洲高清| 97人妻天天添夜夜摸| 日韩制服丝袜自拍偷拍| 91精品三级在线观看| 国产色视频综合| 久久ye,这里只有精品| a级毛片黄视频| 日韩视频在线欧美| 国产男女超爽视频在线观看| 后天国语完整版免费观看| 久久ye,这里只有精品| 人人妻人人澡人人看| av天堂久久9| 亚洲欧美激情在线| 日本五十路高清| 国产亚洲欧美在线一区二区| 精品熟女少妇八av免费久了| 国产99久久九九免费精品| 亚洲熟女毛片儿| 久久久精品免费免费高清| 大话2 男鬼变身卡| 国产精品亚洲av一区麻豆| 国产在视频线精品| 黑人巨大精品欧美一区二区蜜桃| 看免费成人av毛片| 成年人黄色毛片网站| 激情五月婷婷亚洲| a级毛片黄视频| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 欧美黄色淫秽网站| 超碰成人久久| 老汉色∧v一级毛片| 国产91精品成人一区二区三区 | 丁香六月欧美| 女人爽到高潮嗷嗷叫在线视频| 热re99久久精品国产66热6| 午夜福利一区二区在线看| 美女大奶头黄色视频| 大码成人一级视频| 欧美国产精品va在线观看不卡| 国产成人欧美| 蜜桃在线观看..| 精品免费久久久久久久清纯 | 久久人人爽av亚洲精品天堂| 建设人人有责人人尽责人人享有的| 久久毛片免费看一区二区三区| 高清黄色对白视频在线免费看| 日本wwww免费看| 亚洲国产精品一区三区| 五月天丁香电影| 在线观看人妻少妇| 看免费成人av毛片| 亚洲人成电影观看| 国产激情久久老熟女| 精品久久久精品久久久| 精品人妻1区二区| 叶爱在线成人免费视频播放| 久久ye,这里只有精品| 欧美日韩亚洲高清精品| 亚洲国产欧美一区二区综合| 一区二区三区四区激情视频| 精品视频人人做人人爽| 午夜av观看不卡| 精品少妇内射三级| 亚洲成人国产一区在线观看 | 国产免费又黄又爽又色| 久9热在线精品视频| 精品高清国产在线一区| 热99久久久久精品小说推荐| av片东京热男人的天堂| 侵犯人妻中文字幕一二三四区| 岛国毛片在线播放| 一级片'在线观看视频| 国产亚洲一区二区精品| 久久99一区二区三区| 青青草视频在线视频观看| 午夜免费成人在线视频| 嫁个100分男人电影在线观看 | 亚洲黑人精品在线| 天堂8中文在线网| av在线app专区| 久久天躁狠狠躁夜夜2o2o | 国产免费福利视频在线观看| 黄片播放在线免费| 美女中出高潮动态图| 亚洲美女黄色视频免费看| 欧美乱码精品一区二区三区| 高清视频免费观看一区二区| 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 日韩,欧美,国产一区二区三区| 欧美日韩视频精品一区| 亚洲精品av麻豆狂野| 一边摸一边做爽爽视频免费| 久久av网站| 性少妇av在线| 国产一级毛片在线| 1024视频免费在线观看| 午夜福利影视在线免费观看| 中文字幕人妻熟女乱码| 男女高潮啪啪啪动态图| 亚洲欧美日韩另类电影网站| 国产高清不卡午夜福利| 七月丁香在线播放| 秋霞在线观看毛片| 大型av网站在线播放| 久久午夜综合久久蜜桃| 午夜影院在线不卡| 女警被强在线播放| 亚洲三区欧美一区| 纵有疾风起免费观看全集完整版| 亚洲av片天天在线观看| 国产不卡av网站在线观看| 极品人妻少妇av视频| 侵犯人妻中文字幕一二三四区| 夫妻性生交免费视频一级片| 久久热在线av| h视频一区二区三区| 操出白浆在线播放| 一级a爱视频在线免费观看| 女人久久www免费人成看片| √禁漫天堂资源中文www| 自线自在国产av| 国产精品久久久久久人妻精品电影 | 男女午夜视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 免费日韩欧美在线观看| 少妇精品久久久久久久| a级片在线免费高清观看视频| 国产高清视频在线播放一区 | 欧美成人午夜精品| 日韩 欧美 亚洲 中文字幕| 亚洲av欧美aⅴ国产| 日韩av免费高清视频| 美女国产高潮福利片在线看| 日韩视频在线欧美| 99精品久久久久人妻精品| 手机成人av网站| 成人亚洲精品一区在线观看| 免费在线观看日本一区| 日韩 欧美 亚洲 中文字幕| 国产极品粉嫩免费观看在线| 国产精品久久久久久精品古装| 在线看a的网站| av在线播放精品| 成人午夜精彩视频在线观看| 少妇人妻 视频| 亚洲欧美日韩另类电影网站| 国产精品久久久久久人妻精品电影 | 制服诱惑二区| 国产精品国产三级国产专区5o| 成人免费观看视频高清| www.熟女人妻精品国产| 极品少妇高潮喷水抽搐| 久久免费观看电影| 亚洲一区中文字幕在线| 黑丝袜美女国产一区| 人人妻人人添人人爽欧美一区卜| 宅男免费午夜| 亚洲精品国产色婷婷电影| videos熟女内射| 99香蕉大伊视频| 亚洲五月色婷婷综合| 久久女婷五月综合色啪小说| 日日爽夜夜爽网站| 一二三四社区在线视频社区8| 婷婷丁香在线五月| 一个人免费看片子| 这个男人来自地球电影免费观看| 校园人妻丝袜中文字幕| 亚洲国产欧美网| 国产日韩一区二区三区精品不卡| 一本综合久久免费| 亚洲国产精品国产精品| 国产亚洲精品久久久久5区| 国产男人的电影天堂91| 97人妻天天添夜夜摸| 欧美激情 高清一区二区三区| 最黄视频免费看| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 国产片内射在线| 欧美av亚洲av综合av国产av| 最新在线观看一区二区三区 | 又黄又粗又硬又大视频| 国产免费一区二区三区四区乱码| 午夜精品国产一区二区电影| 久久国产精品大桥未久av| e午夜精品久久久久久久| 亚洲国产欧美网| 国产免费又黄又爽又色| 亚洲av国产av综合av卡| av天堂久久9| 久久久久久久大尺度免费视频| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 亚洲熟女毛片儿| 国产成人精品在线电影| 亚洲精品国产一区二区精华液| 亚洲国产欧美日韩在线播放| 国产精品九九99| 亚洲综合色网址| 国产精品一二三区在线看| videosex国产| 亚洲国产欧美在线一区| 在线 av 中文字幕| 人妻 亚洲 视频| 久久天躁狠狠躁夜夜2o2o | 在线亚洲精品国产二区图片欧美| 午夜福利免费观看在线| 操美女的视频在线观看| 少妇 在线观看| 国产精品九九99| 亚洲综合色网址| 免费看av在线观看网站| 美女高潮到喷水免费观看| 国产一区亚洲一区在线观看| 亚洲熟女毛片儿| 精品国产乱码久久久久久男人| 欧美中文综合在线视频| 欧美 日韩 精品 国产| 伊人亚洲综合成人网| 黄色毛片三级朝国网站| 欧美97在线视频| 日本色播在线视频| 欧美精品啪啪一区二区三区 | 成人免费观看视频高清| 视频在线观看一区二区三区| 日日摸夜夜添夜夜爱| 精品亚洲成a人片在线观看| av网站免费在线观看视频| 多毛熟女@视频| 中文字幕亚洲精品专区| 成在线人永久免费视频| 又粗又硬又长又爽又黄的视频| 亚洲国产欧美日韩在线播放| 大话2 男鬼变身卡| 捣出白浆h1v1| 国产精品 欧美亚洲| 日本欧美视频一区| av网站在线播放免费| 在线观看一区二区三区激情| 十分钟在线观看高清视频www| av天堂在线播放| 日本一区二区免费在线视频| 男人舔女人的私密视频| 少妇人妻 视频| 久久人人97超碰香蕉20202| 最黄视频免费看| 悠悠久久av| 如日韩欧美国产精品一区二区三区| 久久久精品区二区三区| 日本色播在线视频| 天天躁夜夜躁狠狠躁躁| 菩萨蛮人人尽说江南好唐韦庄| 成年av动漫网址| 99热国产这里只有精品6| 欧美成人精品欧美一级黄| 日本色播在线视频| 亚洲 国产 在线| 国产欧美日韩综合在线一区二区| 伊人久久大香线蕉亚洲五| 视频区图区小说| 亚洲专区中文字幕在线| 亚洲天堂av无毛| a级毛片黄视频| 国产伦理片在线播放av一区| 各种免费的搞黄视频| 在线观看人妻少妇| 国产男人的电影天堂91| 亚洲av在线观看美女高潮| 这个男人来自地球电影免费观看| 老司机影院成人| 亚洲国产精品一区二区三区在线| a级毛片黄视频| 欧美精品一区二区大全| 狂野欧美激情性bbbbbb| 这个男人来自地球电影免费观看| 妹子高潮喷水视频| 欧美日韩一级在线毛片| 人人妻,人人澡人人爽秒播 | 热re99久久精品国产66热6| 久久精品亚洲av国产电影网| 悠悠久久av| 操出白浆在线播放| av在线app专区| 日韩大片免费观看网站| 久久精品熟女亚洲av麻豆精品| 国产精品三级大全| 亚洲,欧美精品.| 亚洲欧美日韩高清在线视频 | 国产欧美日韩一区二区三 | 亚洲欧美一区二区三区久久| 51午夜福利影视在线观看| 男人爽女人下面视频在线观看| 美女大奶头黄色视频| 97精品久久久久久久久久精品| 这个男人来自地球电影免费观看| 999精品在线视频| 在线观看免费午夜福利视频| 97在线人人人人妻| 自拍欧美九色日韩亚洲蝌蚪91| 巨乳人妻的诱惑在线观看| 美女扒开内裤让男人捅视频| a级片在线免费高清观看视频| 欧美日韩视频高清一区二区三区二| 国产一区二区三区综合在线观看| 丝袜美足系列| 大话2 男鬼变身卡| 777久久人妻少妇嫩草av网站| 久久天躁狠狠躁夜夜2o2o | 亚洲欧洲国产日韩| 欧美av亚洲av综合av国产av| 欧美日韩亚洲国产一区二区在线观看 | 成在线人永久免费视频| 啦啦啦在线免费观看视频4| 另类亚洲欧美激情| 亚洲激情五月婷婷啪啪| 国产精品成人在线| 人妻人人澡人人爽人人| 高潮久久久久久久久久久不卡| 老熟女久久久| 在线观看免费日韩欧美大片| 午夜福利影视在线免费观看| 国产一区二区激情短视频 | 2018国产大陆天天弄谢| 女人高潮潮喷娇喘18禁视频| 后天国语完整版免费观看| 麻豆乱淫一区二区| 99热国产这里只有精品6| 女人高潮潮喷娇喘18禁视频| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 精品国产国语对白av| 国产成人精品久久久久久| 亚洲欧美一区二区三区国产| 国产黄色免费在线视频| 久久久久久久大尺度免费视频| 日韩一卡2卡3卡4卡2021年| 一个人免费看片子| 国产视频首页在线观看| 国产精品九九99| 国精品久久久久久国模美| 久久av网站| 久久亚洲精品不卡| 国产成人一区二区三区免费视频网站 | 丰满迷人的少妇在线观看| 80岁老熟妇乱子伦牲交| 国产成人欧美在线观看 | 激情视频va一区二区三区| 人妻一区二区av| 国产99久久九九免费精品| 色网站视频免费| 五月天丁香电影| 男女免费视频国产| 丰满少妇做爰视频| 九草在线视频观看| 我的亚洲天堂| 亚洲精品一二三| xxxhd国产人妻xxx| 免费观看人在逋| 国产精品免费视频内射| 99香蕉大伊视频| 日本91视频免费播放| 大片电影免费在线观看免费| 亚洲久久久国产精品| 大片免费播放器 马上看| 久久99一区二区三区| 午夜两性在线视频| 国产极品粉嫩免费观看在线| 日韩视频在线欧美| 黄色a级毛片大全视频| 老司机影院毛片| 久久久久网色| 人成视频在线观看免费观看| 一本久久精品| 久久女婷五月综合色啪小说| 国产99久久九九免费精品| 在线精品无人区一区二区三| 久久影院123| 久久av网站| 人人澡人人妻人| 久久久精品免费免费高清| 婷婷丁香在线五月| 精品人妻一区二区三区麻豆| 亚洲国产av影院在线观看| 国产高清视频在线播放一区 | av福利片在线| 国产爽快片一区二区三区| 亚洲成国产人片在线观看| 日韩精品免费视频一区二区三区| 精品一区在线观看国产| 久久av网站| 大片电影免费在线观看免费| 999精品在线视频| 亚洲精品一区蜜桃| 国产国语露脸激情在线看| 天天躁夜夜躁狠狠躁躁| 一区二区av电影网| 午夜精品国产一区二区电影| 亚洲美女黄色视频免费看| 亚洲中文字幕日韩| 亚洲人成电影观看| 一边摸一边做爽爽视频免费| 久久天躁狠狠躁夜夜2o2o | 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 最近最新中文字幕大全免费视频 | 久久久久久久久久久久大奶| 欧美成人午夜精品| av片东京热男人的天堂| 国产不卡av网站在线观看| 国产成人精品久久二区二区免费| 久久久亚洲精品成人影院| 久久久久久人人人人人| 久久久久国产一级毛片高清牌| 亚洲av美国av| 亚洲国产最新在线播放| 这个男人来自地球电影免费观看| 午夜免费成人在线视频| 男女高潮啪啪啪动态图| 夫妻性生交免费视频一级片| 一级黄片播放器| 少妇人妻久久综合中文| 久久毛片免费看一区二区三区| 亚洲中文日韩欧美视频| 国产精品 国内视频| 久久久精品区二区三区| 亚洲精品乱久久久久久| 亚洲久久久国产精品| 免费在线观看视频国产中文字幕亚洲 | 日韩av不卡免费在线播放| 久久久久久久国产电影| 欧美国产精品一级二级三级| 国产精品一区二区在线观看99| 久热爱精品视频在线9| 日本wwww免费看| 免费看不卡的av| av国产久精品久网站免费入址| 丝袜人妻中文字幕| 新久久久久国产一级毛片| 在线观看国产h片| 老司机影院毛片| 97人妻天天添夜夜摸| 亚洲五月婷婷丁香| 欧美 日韩 精品 国产| 精品久久久久久电影网| 国产一区有黄有色的免费视频| 大码成人一级视频| 亚洲精品久久成人aⅴ小说| 国产人伦9x9x在线观看| 人人妻人人爽人人添夜夜欢视频| 人人妻人人添人人爽欧美一区卜| 国产精品九九99| 欧美日韩国产mv在线观看视频| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 亚洲激情五月婷婷啪啪| 纯流量卡能插随身wifi吗| 国产精品欧美亚洲77777| 精品人妻熟女毛片av久久网站| 久久青草综合色| 90打野战视频偷拍视频| 免费日韩欧美在线观看| 最新的欧美精品一区二区| 纯流量卡能插随身wifi吗| 国产精品 欧美亚洲| 天天添夜夜摸| 妹子高潮喷水视频| 国产亚洲欧美精品永久| 色94色欧美一区二区| 侵犯人妻中文字幕一二三四区| 两人在一起打扑克的视频| 深夜精品福利| 美女国产高潮福利片在线看| 亚洲av成人精品一二三区| 欧美老熟妇乱子伦牲交| 又黄又粗又硬又大视频| 亚洲一区二区三区欧美精品| 老汉色∧v一级毛片| 国产成人免费观看mmmm| 国产淫语在线视频| 两个人看的免费小视频| 啦啦啦 在线观看视频| 咕卡用的链子| 伊人亚洲综合成人网| 亚洲av综合色区一区| 日本vs欧美在线观看视频| 日韩精品免费视频一区二区三区| 亚洲国产精品一区三区| videosex国产| 五月天丁香电影| 亚洲五月婷婷丁香| 亚洲av成人不卡在线观看播放网 | 亚洲国产欧美一区二区综合| 高潮久久久久久久久久久不卡| 欧美激情极品国产一区二区三区| 亚洲av片天天在线观看| 香蕉丝袜av| 777久久人妻少妇嫩草av网站| 在线观看免费日韩欧美大片| 成人国产一区最新在线观看 | 欧美日韩亚洲高清精品| 国产97色在线日韩免费| 日韩欧美一区视频在线观看| 天堂8中文在线网| 90打野战视频偷拍视频| 精品少妇一区二区三区视频日本电影| 可以免费在线观看a视频的电影网站| av有码第一页| 性色av一级| 熟女少妇亚洲综合色aaa.| 夫妻性生交免费视频一级片| 国产一级毛片在线|