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

    南海西部夏季上升流鋒面的次中尺度特征分析*

    2020-05-25 04:01:40黃小龍經(jīng)志友鄭瑞璽張旭
    熱帶海洋學(xué)報(bào) 2020年3期

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

    海洋水文學(xué)

    南海西部夏季上升流鋒面的次中尺度特征分析*

    黃小龍1, 2, 經(jīng)志友1, 鄭瑞璽1,2, 張旭1,2

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

    利用衛(wèi)星遙感資料和區(qū)域海洋數(shù)值模式ROMS(regional ocean modeling system)高分辨率數(shù)值模擬結(jié)果, 對(duì)南海西部夏季上升流鋒面的次中尺度特征及其非地轉(zhuǎn)過(guò)程進(jìn)行了探討。高分辨率衛(wèi)星遙感觀測(cè)和數(shù)值模擬結(jié)果顯示, 南海西部夏季鋒面海域存在活躍的次中尺度現(xiàn)象, 其水平尺度約為1~10km, 且具有(1)羅斯貝數(shù)(Rossby number, Ro)的典型次中尺度動(dòng)力學(xué)特征。進(jìn)一步的診斷分析表明, 在夏季西南風(fēng)的驅(qū)動(dòng)下, 沿鋒面射流方向的風(fēng)應(yīng)力(down-front wind stress)引起的跨鋒面??寺斶\(yùn)有利于將海水由鋒面冷水側(cè)向暖水側(cè)輸運(yùn), 減小了鋒面海域的垂向?qū)咏Y(jié)和Ertel位渦, 加劇了鋒面的不穩(wěn)定, 并形成跨鋒面的垂向次級(jí)環(huán)流。高分辨率模擬結(jié)果顯示, 鋒面海域最大垂向流速可達(dá)100m?d–1, 顯著增強(qiáng)了上層海洋的垂向物質(zhì)交換。因此, 活躍在鋒面海域的次中尺度過(guò)程可能是增強(qiáng)南海西部上升流海域垂向物質(zhì)交換的重要貢獻(xiàn)者。

    上升流鋒面; 次中尺度過(guò)程; 區(qū)域海洋數(shù)值模式; 衛(wèi)星遙感; 南海西部

    海洋環(huán)流包含了從幾千公里的大尺度環(huán)流, 到幾百公里的中尺度過(guò)程和幾公里的次中尺度過(guò)程, 直到百米以下的小尺度湍流混合過(guò)程。其中關(guān)于大尺度環(huán)流、中尺度過(guò)程和小尺度湍流混合過(guò)程的研究已取得了豐碩的成果(Munk et al, 1998; Mauritzen et al, 2002; Tian et al, 2009; Chelton et al, 2011; Alford et al, 2015; Wunsch et al, 2018)。但受限于資料分辨率等因素, 目前對(duì)介于中尺度和小尺度之間, 水平尺度為(1~10km), 時(shí)間尺度為幾小時(shí)到數(shù)天的次中尺度過(guò)程尚缺乏足夠的理解。近年來(lái)隨著觀測(cè)技術(shù)的革新、高分辨率模式的發(fā)展以及相關(guān)理論認(rèn)知的提高, 對(duì)普遍存在于上層海洋的次中尺度過(guò)程研究也逐漸深入。

    在中尺度背景環(huán)流、非對(duì)稱(chēng)的大氣強(qiáng)迫(海表風(fēng)場(chǎng)和降溫)以及地形誘導(dǎo)等因素的影響下, 次中尺度過(guò)程在強(qiáng)流、中尺度渦旋邊緣以及鋒面海域尤其活躍(Thomas et al, 2008)。次中尺度過(guò)程具有典型的(1)羅斯貝數(shù)(Rossby number, Ro)和理查森數(shù)(Richardson number, Ri), 表現(xiàn)為完全的三維運(yùn)動(dòng)特征(McWilliams, 2016)。一方面, 現(xiàn)場(chǎng)觀測(cè)和高分辨率數(shù)值研究表明, 次中尺度過(guò)程及其不穩(wěn)定能夠誘發(fā)強(qiáng)烈的垂向非地轉(zhuǎn)次級(jí)環(huán)流, 有利于增強(qiáng)上層海洋熱通量、浮力、營(yíng)養(yǎng)鹽和動(dòng)量的垂向輸運(yùn), 是影響區(qū)域海洋生態(tài)環(huán)境變化的重要?jiǎng)恿^(guò)程之一(Mahadevan et al, 2012; Omand et al, 2015; McKinley et al, 2016; Zhang et al, 2019)。另一方面, 理論分析和觀測(cè)研究顯示, 次中尺度過(guò)程也是平衡態(tài)地轉(zhuǎn)能量向小尺度過(guò)程正向串級(jí)的重要媒介過(guò)程之一(Capet et al, 2008a)。次中尺度不穩(wěn)定(例如對(duì)稱(chēng)不穩(wěn)定和混合層斜壓不穩(wěn)定)能夠有效地汲取地轉(zhuǎn)動(dòng)能并釋放鋒面有效位能, 并在次級(jí)不穩(wěn)定的作用下將能量進(jìn)一步串級(jí)至小尺度湍流混合過(guò)程(Boccaletti et al, 2007; D'Asaro et al, 2011; Barkan et al, 2015)。因此, 次中尺度過(guò)程是中、小尺度相互作用的關(guān)鍵動(dòng)力過(guò)程。

    南海地處東亞季風(fēng)區(qū)(Fang et al, 2002), 在大氣強(qiáng)迫、復(fù)雜海底地形等因素的影響下, 具有多尺度的動(dòng)力過(guò)程(Wang et al, 2003; Chen et al, 2011; Zhang et al, 2016), 是研究次中尺度過(guò)程和中、小尺度相互作用的理想海域(Jing et al, 2016; Zhong et al, 2017; 冀承振等, 2017)。其中南海西部是南海中尺度過(guò)程活躍的主要區(qū)域之一(Zhuang et al, 2010), 且已被現(xiàn)場(chǎng)觀測(cè)所證實(shí)。衛(wèi)星遙感觀測(cè)顯示(圖1), 南海西部海面高度變化(圖1中黑色框)顯著高于中部海盆區(qū), 和南海東北部海域同為南海兩個(gè)典型的中尺度渦旋活躍海區(qū)(Hu et al, 2011; Zhang et al, 2016)。在夏季西南風(fēng)的驅(qū)動(dòng)下, 南海西部出現(xiàn)顯著的東向離岸上升流(Chen et al, 2014a; Li et al, 2014)。同時(shí), 沿岸徑流輸入(Chen et al, 2012)、風(fēng)應(yīng)力旋度(Xie et al, 2003, 2007)以及渦度水平輸運(yùn)(Wang et al, 2006; Gan et al, 2008)等均對(duì)該上升流的生成與不同時(shí)間尺度的演變有著重要影響。上升流將下層低溫、高營(yíng)養(yǎng)鹽的海水帶至表層, 并在平流作用下向東輸運(yùn), 在南海西部形成顯著的溫度鋒面和葉綠素濃度高值區(qū)(Kuo et al, 2000; Chen et al, 2014b)。高分辨率水色遙感觀測(cè)結(jié)果顯示, 南海西部夏季鋒面海域廣泛存在著細(xì)絲狀和渦旋狀的次中尺度結(jié)構(gòu)(Liu et al, 2015; Yu et al, 2018), 但關(guān)于該海域次中尺度過(guò)程的動(dòng)力學(xué)特征及其對(duì)上層海洋的影響尚缺乏足夠的理解。本文嘗試?yán)酶叻直媛蕯?shù)值模擬并結(jié)合衛(wèi)星觀測(cè)資料分析, 探討南海西部夏季上升流鋒面的次中尺度特征及其非地轉(zhuǎn)過(guò)程, 以期為南海西部鋒面海域的高葉綠素濃度和中尺度過(guò)程變異提供動(dòng)力學(xué)解釋, 進(jìn)一步加深對(duì)南海西部多尺度相互作用和生態(tài)效應(yīng)的認(rèn)識(shí)。

    圖1 南海海表面高度均方根的空間分布

    黑色矩形框表示本文研究區(qū)域; 地圖來(lái)自Matlab軟件自帶底圖, 下同

    Fig. 1 Spatial distribution of SSH root mean square (rms) in the South China Sea (SCS).The research domain is delineated by the black box

    1 數(shù)據(jù)來(lái)源與分析方法

    1.1 衛(wèi)星遙感資料

    本文使用的衛(wèi)星遙感資料包括高分辨率海表溫度數(shù)據(jù)(sea surface temperature, SST)、海表面高度(sea surface height, SSH)、海表高度異常(sea level anomaly, SLA)、風(fēng)場(chǎng)數(shù)據(jù)和葉綠素濃度數(shù)據(jù)。海表溫度數(shù)據(jù)(2006—2017年)采用英國(guó)氣象局制作并由美國(guó)國(guó)家海洋數(shù)據(jù)中心發(fā)布的高分辨率融合產(chǎn)品(http://data.nodc.noaa.gov/ghrsst/L4/GLOB/UKMO/OSTIA/), 空間分辨率為0.05°×0.05°, 時(shí)間分辨率為1d, 反演精度約為0.57℃(Dong et al, 2014)。日平均的高度計(jì)數(shù)據(jù)(SSH和SLA,1993—2016年)來(lái)源于法國(guó)國(guó)家空間研究中心提供的AVISO(Archiving Validation and Interpolation of Satellite Oceanographic)網(wǎng)格化產(chǎn)品(ftp://ftp.aviso.oceanobs. com/global/), 空間分辨率為0.25°×0.25°, 反演精度約為0.02m。日平均海表風(fēng)場(chǎng)數(shù)據(jù)(1999—2009年)來(lái)自QuikSCAT(Quick Scatterometer)衛(wèi)星散射計(jì), 空間分辨率為0.25°×0.25°, 資料精度為1m?s–1。海表葉綠素濃度數(shù)據(jù)(2002—2012年)來(lái)自中等分辨率成像光譜儀(MEdium Resolution Imaging Spectrometer, MERIS)提供的2級(jí)產(chǎn)品數(shù)據(jù)(ftp:// merisfrs-fts-ds.eo.esa.int), 水平分辨率約為300m, 絕對(duì)精確度達(dá)到10–3mg·m–3。

    1.2 模式配置

    本文利用ROMS(regional oceanic modeling system)數(shù)值模式對(duì)南海及其西部海域進(jìn)行高分辨率的嵌套模擬。其中最外層模擬區(qū)域覆蓋西太平洋(ROMS0, 7.5km分辨率), 第一層嵌套覆蓋南海區(qū)域(ROMS1, 1.5km分辨率, 如圖2所示), 第二層嵌套聚焦南海西部海域(ROMS2, 500m分辨率); 模式垂向分層均為60層, 并在上邊界和底邊界適當(dāng)加密。初始邊界條件以及風(fēng)場(chǎng)強(qiáng)迫分別采用SODA(simple ocean data assimilation data)海洋數(shù)據(jù)集的氣候態(tài)月平均數(shù)據(jù)(1990—2010年, Carton et al, 2008)和氣候態(tài)日平均QuikSCAT(Quick Scatterometer)風(fēng)場(chǎng)數(shù)據(jù)(Risien et al, 2008)。其中, 采用的SODA數(shù)據(jù)集包含溫度、鹽度、水平和垂直速度以及SSH等變量, 該再分析同化資料的空間分辨率為0.5°×0.5°, 對(duì)應(yīng)的數(shù)據(jù)精度約為: 溫度0.41℃、鹽度0.5‰、速度不確定性小于30%和SSH 0.03m(Carton et al, 2000)。熱通量和淡水通量等氣候態(tài)月平均?!?dú)馔縼?lái)源于國(guó)際海洋大氣綜合數(shù)據(jù)集(International Comprehensive Ocean-Atmosphere Data Set, ICOADS) (Woodruff et al, 2011), 空間分辨率為1°×1°。地形數(shù)據(jù)采用美國(guó)國(guó)家海洋和大氣管理局(National Oceanic and Atmospheric Administration, NOAA)提供的ETOPO2數(shù)據(jù)。小尺度湍流混合采用KPP參數(shù)化方案(K-profile parameterization)進(jìn)行參數(shù)化(Large et al, 1994)。ROMS0診斷計(jì)算20年后, 在線嵌套的ROMS1和ROMS2計(jì)算第21年并輸出日平均模擬結(jié)果。本文選取500m分辨率的ROMS2模擬結(jié)果對(duì)南海西部上升流鋒面海域的次中尺度動(dòng)力特征進(jìn)行分析。

    圖2 南海嵌套模擬區(qū)域示意圖

    第一層嵌套(ROMS1)包括南海海域, 第二層嵌套(ROMS2)為南海西部海域(黑色矩形框)

    Fig. 2 The diagram for the nested ROMS model domain. The western SCS is the nested domain for the ROMS2 simulation (black box) while the nested domain of ROMS1 covers the SCS

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

    位渦(potential vorticity, PV)作為一種重要的動(dòng)力學(xué)示蹤物, 反映海洋環(huán)流和多尺度動(dòng)力過(guò)程的穩(wěn)定性(Hoskins, 1974; Thomas et al, 2013)。由于北半球恒為正值, 因此當(dāng)Ertel位渦小于零時(shí), 易誘發(fā)對(duì)稱(chēng)不穩(wěn)定(Thomas et al, 2013)。在地轉(zhuǎn)平衡的假設(shè)條件下, Ertel位渦定義如下(Hoskins, 1974):

    為評(píng)估風(fēng)應(yīng)力引起的艾克曼輸運(yùn)對(duì)改變海表浮力的作用, 本文計(jì)算了埃克曼浮力通量(Ekman buoyancy flux, EBF), 表達(dá)式如下:

    2 結(jié)果分析

    2.1 鋒面的次中尺度動(dòng)力學(xué)特征

    高分辨率MERIS水色遙感資料顯示, 南海西部上升流海域有著顯著的葉綠素濃度高值區(qū), 側(cè)向?qū)挾燃s20~40km的高葉綠素帶狀區(qū)域由近岸109°12'E向東延伸至112°30'E, 并且在空間分布形態(tài)上與衛(wèi)星觀測(cè)的溫度鋒面位置一致(圖3)。同時(shí), 從高分辨率葉綠素濃度圖像可以清晰地看到鋒面兩側(cè)渦旋狀和渦絲狀的次中尺度結(jié)構(gòu), 這些次中尺度現(xiàn)象的側(cè)向尺度約為10km, 該特征尺度與呂宋海峽、南海西部等海域的現(xiàn)場(chǎng)觀測(cè)結(jié)果一致(Zheng et al, 2008; Song et al, 2019)。

    圖3 南海西部2009年7月26日的海表葉綠素濃度、海表風(fēng)速(箭頭)和海表溫度(等值線, 單位: ℃)平面分布

    本文選取11°48'—12°06'N, 110°30'—110°54'E的典型鋒面區(qū)域(圖4a虛線矩形框)對(duì)次中尺度過(guò)程進(jìn)行進(jìn)一步分析, 其水平分布和跨鋒面結(jié)構(gòu)如圖5所示。沿流軸方向露頭等密線顯示, 流軸附近存在寬約5km的狹長(zhǎng)密度鋒面。垂向上, 密度鋒面主要存在于混合層內(nèi)20m以淺, 且南側(cè)等密線相對(duì)北側(cè)更加密集。

    圖4 模式第21年8月4日南海西部數(shù)值模擬海表溫度、表層流(箭頭)、SLA(等值線, 單位: m)(a)和羅斯貝數(shù)Ro(b)分布

    黑色實(shí)線表示流軸位置(相對(duì)渦度為0); a中虛線矩形框?yàn)楸疚倪x取的典型鋒面區(qū)域

    Fig. 4 Maps of SST (a) and Ro (b) from the ROMS2 simulation. The vectors are surface currents. The gray contours represent the SLA, and the black line is the jet axis

    圖5 典型鋒面區(qū)域鋒面強(qiáng)度F、地轉(zhuǎn)流ug(箭頭)和海水密度(灰色等值線, 單位: kg·m–3)的水平分布(a)和斷面分布(b)

    a中紫紅色實(shí)線表示斷面位置, 黑色實(shí)線表示流軸; b中黑色實(shí)線表示混合層深度

    Fig. 5 Spatial distribution of the intensity of front: horizontal distribution (a) and vertical profile (b). The geostrophic flows (vector) and density (grey contour) are also shown. Black line is the jet axis in (a) and the mixed layer depth in (b). The location of the vertical profile is shown by the magenta line in (a)

    2.2 垂向次級(jí)環(huán)流與位渦分布特征

    為了進(jìn)一步討論上升流鋒面區(qū)域垂向環(huán)流結(jié)構(gòu), 本文通過(guò)空間濾波的方法提取次中尺度的流速信號(hào)。圖6給出了次中尺度垂向流速異常的水平和垂向分布??梢源_定, 鋒面附近存在明顯的跨鋒面的非地轉(zhuǎn)次級(jí)環(huán)流, 其表層海水在鋒面區(qū)域輻聚下沉, 而在其兩側(cè)存在上升流, 最大垂向流速出現(xiàn)在15m深度附近(圖6b), 其中最大垂向流速可達(dá)100m?d–1, 顯著增強(qiáng)了上層海洋物質(zhì)的垂向交換。該特征與利用南海西部、南海北部等海域的現(xiàn)場(chǎng)觀測(cè)及高分辨率模擬診斷的結(jié)果一致(Xie et al, 2017; Zhong et al, 2017)。

    理論研究結(jié)果表明:局地非地轉(zhuǎn)次級(jí)環(huán)流與次中尺度不穩(wěn)定密切相關(guān)(Mahadevan et al, 2006; Gula et al, 2014; McWilliams, 2017)。當(dāng)Ertel位渦與反號(hào)時(shí), 有利于促發(fā)多種類(lèi)型的次中尺度不穩(wěn)定, 其中對(duì)稱(chēng)不穩(wěn)定是能量增長(zhǎng)最快的一種模態(tài)(Hoskins et al, 1972)。與傳統(tǒng)的中尺度過(guò)程不同, 鋒面弱垂向?qū)踊?qiáng)水平浮力梯度以及垂向剪切等均使得位渦的水平斜壓分量(EPVh)與其垂向分量(EPVv)量值相當(dāng), 有利于負(fù)Ertel位渦的發(fā)生并易于引發(fā)次中尺度不穩(wěn)定。

    圖7顯示了基于高分辨率數(shù)值模擬診斷計(jì)算得到的Ertel位渦、水平位渦分量和垂向位渦分量的空間分布。負(fù)Ertel位渦主要出現(xiàn)在上層海洋次中尺度鋒面附近, 表明鋒面區(qū)域發(fā)生次中尺度斜壓不穩(wěn)定。通過(guò)比較各位渦分量可以發(fā)現(xiàn), Ertel位渦水平分量絕對(duì)值大于水平分量, 使得Ertel位渦小于零(圖7a—c)。由Ertel位渦的跨鋒面垂向結(jié)構(gòu)分布可知, 受次級(jí)環(huán)流作用, 海表的負(fù)位渦水體沿等密度面自海表進(jìn)入混合層內(nèi); 下層高位渦海水受鋒面兩側(cè)上升流影響抬升至海表, 表明表層和下層海水在次級(jí)環(huán)流的作用下發(fā)生垂向交換且浮力發(fā)生再分配(圖7d)。在此過(guò)程中, 增強(qiáng)的垂向次級(jí)環(huán)流能夠有效地將下層較豐富的營(yíng)養(yǎng)鹽帶進(jìn)入海表, 為鋒面海域的浮游生物生長(zhǎng)提供營(yíng)養(yǎng)鹽供給。

    圖6 10m層垂向流速異常的水平和垂向分布

    a. 垂向流速異常′(填色)、水平流速異?!?箭頭)和海水密度(等值線, 單位: kg·m–3)的水平分布; b. 斷面溫度(填色)、跨鋒面和垂向流速異常(箭頭)、海水密度(等值線, 單位: kg·m–3)分布; a中紫紅色實(shí)線表示斷面位置; 黑色實(shí)線表示流軸; b中黑色實(shí)線表示混合層深度

    Fig. 6 Map (a) and vertical (b) profile of vertical velocity anomaly at 10 m. (a) Black arrows refer to surface currents anomalies. Magenta line indicates the section for the vertical profile, and black line indicates the jet axis. (b) Shading denotes temperature, and black vector represents the cross-front and vertical velocity anomaly. Black line represents mixed layer depth. Grey contour is for density

    圖7 典型鋒面海域Ertel位渦(填色)、表層流u(箭頭)和海水密度(灰色等值線)分布

    a—c分別表示Ertel位渦(EPV)、Ertel位渦水平分量(EPVh)和Ertel位渦垂向分量(EPVv)的水平分布, 黑色實(shí)線表示流軸; d為Ertel位渦斷面分布, 斷面位置如a中紫紅色實(shí)線所示, 黑色實(shí)線表示混合層深度, 灰色等值線為密度(單位: kg·m–3)

    Fig. 7 Snapshot of the Ertel PV. (a-c) are the EPV, EPVh and EPVv, respectively. The black line represents the jet axis. Black arrows denote surface currents. (d) is vertical distribution of EPV. The location of the vertical section is shown by the magenta line in (a). Black line represents mixed layer depth

    海表外強(qiáng)迫是上層海洋位渦發(fā)生改變的主要源之一(Thomas, 2005; Thomas et al, 2005)。在夏季西南風(fēng)強(qiáng)迫引起的浮力損失等非保守過(guò)程影響下, 南海西部鋒面海域易出現(xiàn)負(fù)Ertel位渦。診斷計(jì)算結(jié)果顯示, 在沿鋒面射流風(fēng)應(yīng)力(down-front wind stress)的作用下, 最大??寺×ν砍霈F(xiàn)在流軸附近, 大小約為1×10–5m2?s–3, 表明沿鋒面射流方向的風(fēng)應(yīng)力引起的??寺斶\(yùn)將冷水向暖水側(cè)跨鋒面輸運(yùn)(圖8), 進(jìn)一步改變鋒面浮力分布并減弱垂向?qū)咏Y(jié), 加強(qiáng)鋒面的水平浮力梯度以及減小Ertel位渦水平斜壓分量, 最終導(dǎo)致負(fù)Ertel位渦出現(xiàn), 從而加劇鋒面的不穩(wěn)定。

    圖8 典型鋒面海域的埃克曼浮力通量(填色)、地轉(zhuǎn)流ug(黑色箭頭)、風(fēng)應(yīng)力tstr(紫紅色箭頭)和海水密度(等值線,單位: kg·m–3)的水平分布

    從空間分布形態(tài)來(lái)看, ??寺×ν靠臻g分布與Ertel位渦結(jié)構(gòu)也有著良好的對(duì)應(yīng)關(guān)系。在一定程度上, 這些結(jié)果表明沿鋒面射流風(fēng)應(yīng)力驅(qū)動(dòng)的跨鋒面??寺斶\(yùn)過(guò)程可能是導(dǎo)致夏季南海西部鋒面海域發(fā)生次中尺度對(duì)稱(chēng)不穩(wěn)定的重要原因。

    3 結(jié)論

    本文利用高分辨率衛(wèi)星遙感資料并結(jié)合數(shù)值模擬結(jié)果, 分析了夏季南海西部上升流鋒面的次中尺度特征及其非地轉(zhuǎn)過(guò)程。高分辨率水色遙感觀測(cè)顯示, 南海西部上升流鋒面海域存在豐富的次中尺度現(xiàn)象。高分辨率模擬結(jié)果發(fā)現(xiàn), 上升流鋒面有顯著的非地轉(zhuǎn)特征并伴隨次中尺度過(guò)程發(fā)生。診斷計(jì)算結(jié)果顯示, 跨鋒面的非線性的埃克曼水平凈輸運(yùn)將海水自鋒面冷水側(cè)向暖水側(cè)輸運(yùn), 增強(qiáng)了鋒面水平浮力梯度并減弱其垂向?qū)咏Y(jié), 導(dǎo)致鋒面海域出現(xiàn)負(fù)Ertel位渦, 為次中尺度對(duì)稱(chēng)不穩(wěn)定提供了有利條件。因此, 沿鋒面射流風(fēng)應(yīng)力引起的浮力損失可能是導(dǎo)致南海西部上升流海域鋒面發(fā)生次中尺度對(duì)稱(chēng)不穩(wěn)定的重要機(jī)制之一。同時(shí), 診斷結(jié)果顯示, 次中尺度過(guò)程能夠引起跨鋒面的垂向次級(jí)環(huán)流, 最大垂向流速可達(dá)100m?d–1, 顯著增強(qiáng)了鋒面海域的垂向物質(zhì)交換。本文基于數(shù)值模擬初步探討了南海西部鋒面海域次中尺度過(guò)程, 關(guān)于其形成與演變的動(dòng)力過(guò)程, 還需要結(jié)合現(xiàn)場(chǎng)觀測(cè)資料進(jìn)行深入研究。

    范植松, 2002. 海洋內(nèi)部混合研究基礎(chǔ)[M]. 北京: 海洋出版社: 7–10. FAN ZHISONG, 2002. Research fundamentals of ocean interior mixing[M]. Beijing: China Ocean Press: 7–10 (in Chinese).

    冀承振, 葉瑞杰, 董濟(jì)海, 等, 2017. 南海中尺度渦邊緣亞中尺度過(guò)程模式研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào), 47(1): 1–6. JI CHENGZHEN, YE RUIJIE, DONG JIHAI, et al, 2017. The simulation of submesoscale process at the periphery of a mesoscale eddy in the South China Sea[J]. Periodical of Ocean University of China, 47(1): 1–6 (in Chinese with English abstract).

    ALFORD M H, PEACOCK T, MACKINNON J A, et al, 2015. The formation and fate of internal waves in the South China Sea[J]. Nature, 521(7550): 65–69.

    BARKAN R, WINTERS K B, SMITH S G L, 2015. Energy cascades and loss of balance in a reentrant channel forced by wind stress and buoyancy fluxes[J]. Journal of Physical Oceanography, 45(1): 272–293.

    BOCCALETTI G, FERRARI R, FOX-KEMPER B, 2007. Mixed layer instabilities and restratification[J]. Journal of Physical Oceanography, 37(9): 2228–2250.

    CAPET X, MCWILLIAMS J C, MOLEMAKER M J, et al, 2008a. Mesoscale to submesoscale transition in the California current system. Part I: flow structure, eddy flux, and observational tests[J]. Journal of Physical Oceanography, 38(1): 29–43.

    CAPET X, MCWILLIAMS J C, MOLEMAKER M J, et al, 2008b. Mesoscale to submesoscale transition in the California current system. Part II: frontal processes[J]. Journal of Physical Oceanography, 38(1): 44–64.

    CARTON J A, CHEPURIN G, CAO X H, 2000. A Simple Ocean Data Assimilation analysis of the global upper ocean 1950–95. Part II: results[J]. Journal of Physical Oceanography, 30(2): 311–326.

    CARTON J A, GIESE B S, 2008. A reanalysis of ocean climate using Simple Ocean Data Assimilation (SODA)[J]. Monthly Weather Review, 136(8): 2999–3017.

    CHELTON D B, SCHLAX M G, SAMELSON R M, 2011. Global observations of nonlinear mesoscale eddies[J]. Progress in Oceanography, 91(2): 167–216.

    CHEN CHANGLIN, WANG GUIHUA, 2014a. Interannual variability of the eastward current in the western South China Sea associated with the summer Asian monsoon[J]. Journal of Geophysical Research: Oceans, 119(9): 5745–5754.

    CHEN CHANGSHENG, LAI ZHIGANG, BEARDSLEY R C, et al, 2012. Current separation and upwelling over the southeast shelf of Vietnam in the South China Sea[J]. Journal of Geophysical Research: Oceans, 117(C3): C03033.

    CHEN GENGXIN, HOU YIJUN, CHU XIAOQING, 2011. Mesoscale eddies in the South China Sea: mean properties, spatiotemporal variability, and impact on thermohaline structure[J]. Journal of Geophysical Research: Oceans, 116(C6): C06018.

    CHEN GENGXIN, XIU PENG, CHAI FEI, 2014b. Physical and biological controls on the summer chlorophyll bloom to the east of Vietnam[J]. Journal of Oceanography, 70(3): 323–328.

    D'ASARO E, LEE C, RAINVILLE L, et al, 2011. Enhanced turbulence and energy dissipation at ocean fronts[J]. Science, 332(6027): 318–322.

    DONG CHANGMING, MCWILLIAMS J C, LIU YU, et al, 2014. Global heat and salt transports by eddy movement[J]. Nature Communications, 5: 3294.

    FANG WENDONG, FANG GUOHONG, SHI PING, et al, 2002. Seasonal structures of upper layer circulation in the southern South China Sea from in situ observations[J]. Journal of Geophysical Research: Oceans, 107(C11): 23-1–23-12.

    GAN JIANPING, QU TANGDONG, 2008. Coastal jet separation and associated flow variability in the southwest South China Sea[J]. Deep Sea Research Part I: Oceanographic Research Papers, 55(1): 1–19.

    GULA J, MOLEMAKER M J, MCWILLIAMS J C, 2014. Submesoscale cold filaments in the Gulf Stream[J]. Journal of Physical Oceanography, 44(10): 2617–2643.

    HOSKINS B J, 1974. The role of potential vorticity in symmetric stability and instability[J]. Quarterly Journal of the Royal Meteorological Society, 100(425): 480–482.

    HOSKINS B J, BRETHERTON F P, 1972. Atmospheric frontogenesis models: mathematical formulation and solution[J]. Journal of the Atmospheric Sciences, 29(1): 11–37.

    HU JIANYU, GAN JIANPING, SUN ZHENYU, et al, 2011. Observed three-dimensional structure of a cold eddy in the southwestern South China Sea[J]. Journal of Geophysical Research: Oceans, 116(C5): C05016.

    JING ZHIYOU, QI YIQUAN, FOX-KEMPER B, et al, 2016. Seasonal thermal fronts on the northern South China Sea shelf: satellite measurements and three repeated field surveys[J]. Journal of Geophysical Research: Oceans, 121(3): 1914–1930.

    KUO NAN-JUNG, ZHENG QUANAN, HO C-R, 2000. Satellite observation of upwelling along the western coast of the South China Sea[J]. Remote Sensing of Environment, 74(3): 463–470.

    LARGE W G, MCWILLIAMS J C, DONEY S C, 1994. Oceanic vertical mixing: a review and a model with a nonlocal boundary layer parameterization[J]. Reviews of Geophysics, 32(4): 363–403.

    LI YUANLONG, HAN WEIQING, WILKIN J L, et al, 2014. Interannual variability of the surface summertime eastward jet in the South China Sea[J]. Journal of Geophysical Research: Oceans, 119(10): 7205–7228.

    LIU FENFEN, TANG SHILIN, CHEN CHUQUN, 2015. Satellite observations of the small-scale cyclonic eddies in the western South China Sea[J]. Biogeosciences, 12(2): 299–305.

    MAHADEVAN A, TANDON A, 2006. An analysis of mechanisms for submesoscale vertical motion at ocean fronts[J]. Ocean Modelling, 14(3–4): 241–256.

    MAHADEVAN A, D'ASARO E, LEE C, et al, 2012. Eddy-driven stratification initiates north Atlantic spring phytoplankton blooms[J]. Science, 337(6090): 54–58.

    MAURITZEN C, POLZIN K L, MCCARTNEY M S, et al, 2002. Evidence in hydrography and density fine structure for enhanced vertical mixing over the Mid-Atlantic Ridge in the western Atlantic[J]. Journal of Geophysical Research: Oceans, 107(C10): 11-1–11-19.

    MCKINLEY G A, PILCHER D J, FAY A R, et al, 2016. Timescales for detection of trends in the ocean carbon sink[J]. Nature, 530(7591): 469–472.

    MCWILLIAMS J C, 2016. Submesoscale currents in the ocean[J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 472(2189): 20160117.

    MCWILLIAMS J C, 2017. Submesoscale surface fronts and filaments: secondary circulation, buoyancy flux, and frontogenesis[J]. Journal of Fluid Mechanics, 823: 391–432.

    MUNK W, WUNSCH C, 1998. Abyssal recipes II: energetics of tidal and wind mixing[J]. Deep Sea Research Part I: Oceanographic Research Papers, 45(12): 1977–2010.

    OMAND M M, D'ASARO E A, LEE C M, et al, 2015. Eddy-driven subduction exports particulate organic carbon from the spring bloom[J]. Science, 348(6231): 222–225.

    RISIEN C M, CHELTON D B, 2008. A global climatology of surface wind and wind stress fields from eight years of QuikSCAT Scatterometer data[J]. Journal of Physical Oceanography, 38(11): 2379–2413.

    SONG BO, WANG HUIZAN, CHEN CHANGLIN, et al, 2019. Observed subsurface eddies near the Vietnam coast of the South China Sea[J]. Acta Oceanologica Sinica, 38(4): 39–46.

    THOMAS L N, 2005. Destruction of potential vorticity by winds[J]. Journal of Physical Oceanography, 35(12): 2457–2466.

    THOMAS L N, LEE C M, 2005. Intensification of ocean fronts by down-front winds[J]. Journal of Physical Oceanography, 35(6): 1086–1102.

    THOMAS L N, TANDON A, MAHADEVAN A, 2008. Submesoscale processes and dynamics[M]//HECHT M W, HASUMI H. Ocean modeling in an eddying regime. Washington: American Geophysical Union: 17–38.

    THOMAS L N, TAYLOR J R, FERRARI R, et al, 2013. Symmetric instability in the Gulf Stream[J]. Deep Sea Research Part II: Topical Studies in Oceanography, 91: 96–110.

    TIAN JIWEI, YANG QINGXUAN, ZHAO WEI, 2009. Enhanced diapycnal mixing in the South China Sea[J]. Journal of Physical Oceanography, 39(12): 3191–3203.

    WANG GUIHUA, SU JILAN, CHU P C, 2003. Mesoscale eddies in the South China Sea observed with altimeter data[J]. Geophysical Research Letters, 30(21): 2121.

    WANG GUIHUA, CHEN DAKE, SU JILAN, 2006. Generation and life cycle of the dipole in the South China Sea summer circulation[J]. Journal of Geophysical Research: Oceans, 111(C6): C06002.

    WOODRUFF S D, WORLEY S J, LUBKER S J, et al, 2011. ICOADS release 2.5: extensions and enhancements to the surface marine meteorological archive[J]. International Journal of Climatology, 31(7): 951–967.

    WUNSCH C, FERRARI R, 2018. 100 years of the ocean general circulation[J]. Meteorological Monographs, 59: 7.1–7.32.

    XIE LINGLING, PALLàS-SANZ E, ZHENG QUANAN, et al, 2017. Diagnosis of 3D vertical circulation in the upwelling and frontal zones east of Hainan Island, China[J]. Journal of Physical Oceanography, 47(7): 755–774

    XIE SHANG-PING, XIE QIANG, WANG DONGXIAO, et al, 2003. Summer upwelling in the South China Sea and its role in regional climate variations[J]. Journal of Geophysical Research: Oceans, 108(C8): 3261.

    XIE SHANG-PING, CHANG C-H, XIE QIANG, et al, 2007. Intraseasonal variability in the summer South China Sea: Wind jet, cold filament, and recirculations[J]. Journal of Geophysical Research: Oceans, 112(C10): C10008.

    YU JIE, ZHENG QUANAN, JING ZHIYOU, et al, 2018. Satellite observations of sub-mesoscale vortex trains in the western boundary of the South China Sea[J]. Journal of Marine Systems, 183: 56–62.

    ZHANG ZHENGGUANG, QIU BO, KLEIN P, et al, 2019. The influence of geostrophic strain on oceanic ageostrophic motion and surface chlorophyll[J]. Nature Communications, 10: 2838.

    ZHANG ZHIWEI, TIAN JIWEI, QIU BO, et al, 2016. Observed 3D structure, generation, and dissipation of oceanic mesoscale eddies in the South China Sea[J]. Scientific Reports, 6: 24349.

    ZHENG QUANAN, LIN HUI, MENG JUNMIN, et al, 2008. Sub-mesoscale ocean vortex trains in the Luzon Strait[J]. Journal of Geophysical Research: Oceans, 113(C4): C04032.

    ZHONG YISEN, BRACCO A, TIAN JIWEI, et al, 2017. Observed and simulated submesoscale vertical pump of an anticyclonic eddy in the South China Sea[J]. Scientific Reports, 7: 44011.

    ZHUANG WEI, XIE SHANG‐PING, WANG DONGXIAO, et al, 2010. Intraseasonal variability in sea surface height over the South China Sea[J]. Journal of Geophysical Research: Oceans, 115(C4): C04010.

    Analysis of submesoscale characteristics of summer upwelling fronts in the western South China Sea

    HUANG Xiaolong1,2, JING Zhiyou1, ZHENG Ruixi1,2, ZHANG Xu1,2

    1. State Key Laboratory of Tropical Oceanography (South China Sea Institute of Oceanology, Chinese Academy of Sciences), Guangzhou 510301, China; 2. University of Chinese Academy of Sciences, Beijing 100049, China

    We investigate submesoscale characteristics of summer upwelling fronts in the western South China Sea (WSCS) and associated ageostrophic processes by using satellite measurements and high-resolution ROMS simulations. Active submesoscale filaments with a typical horizontal scale of(1-10) km are detected to be characterized by(1) Rossby number (Ro) from the fine-resolution satellite images and simulation results. The diagnostic analysis shows that down-front wind forcing drives a net cross-front Ekman transport and advects heavy water over light water. This process at submesoscale tends to reduce the stratification and potential vorticity (PV), exacerbates the frontal instabilities, and forms the cross-front secondary circulation. The high-resolution simulation results show that the maximum vertical velocity in the frontal zone can reach 100m?d–1, which significantly enhances vertical material exchange. In this context, active submesoscale processes may contribute to enhanced vertical exchanges of the upper ocean in the summer upwelling front of the western SCS.

    upwelling front; submesoscale process; ROMS; satellite remote sensing; western South China Sea

    P731.21

    A

    1009-5470(2020)03-0001-09

    10.11978/2019086

    http://www.jto.ac.cn

    2019-09-12;

    2019-12-20。

    林強(qiáng)編輯

    中國(guó)科學(xué)院基礎(chǔ)前沿科學(xué)研究計(jì)劃原始創(chuàng)新項(xiàng)目(ZDBS-LY-DQC011); 國(guó)家自然科學(xué)基金項(xiàng)目(41776040); 中國(guó)科學(xué)院南海生態(tài)環(huán)境工程創(chuàng)新研究院自主部署項(xiàng)目(ISEE2018PY05); 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室OCFL功能實(shí)驗(yàn)室開(kāi)放課題(OCFL-201804); 廣州市科學(xué)研究計(jì)劃(201904010420)。

    黃小龍(1992—), 男, 廣東省信宜市人, 碩士研究生, 主要從事上層海洋次中尺度過(guò)程研究。E-mail: huangxiaolong@scsio.ac.cn

    經(jīng)志友。E-mail: jingzhiyou@scsio.ac.cn

    *感謝三位審稿專(zhuān)家提出的寶貴修改意見(jiàn)。

    2019-09-12;

    2019-12-20.

    Editor: LIN Qiang

    Original Innovation Project of Basic Frontier Scientific Research Program of CAS (ZDBS-LY-DQC011); National Natural Science Foundation of China (41776040); Innovation Academy of South China Sea Ecology and Environmental Engineering, CAS (ISEE2018PY05); Laboratory for Ocean Dynamics and Climate, Pilot Qingdao National Laboratory for Marine Science and Technology (OCFL-201804); Guangzhou Science and Technology Project (201904010420).

    JING Zhiyou. E-mail: jingzhiyou@scsio.ac.cn

    97人妻天天添夜夜摸| 80岁老熟妇乱子伦牲交| 亚洲精品一区蜜桃| 精品亚洲成国产av| av线在线观看网站| 自线自在国产av| 午夜精品国产一区二区电影| 国产成人欧美| 亚洲国产欧美日韩在线播放| 悠悠久久av| 精品亚洲成国产av| 国产成人精品福利久久| 国产黄频视频在线观看| xxxhd国产人妻xxx| 久久久久久人妻| 蜜桃在线观看..| 男的添女的下面高潮视频| 国产成人欧美在线观看 | 中文字幕制服av| 国产av一区二区精品久久| 老汉色av国产亚洲站长工具| 国产欧美亚洲国产| 伦理电影大哥的女人| 黄片小视频在线播放| 国产一卡二卡三卡精品 | 国产熟女午夜一区二区三区| 日本wwww免费看| 少妇人妻久久综合中文| 性少妇av在线| 大片电影免费在线观看免费| 亚洲欧洲精品一区二区精品久久久 | 中文天堂在线官网| 亚洲人成电影观看| 午夜激情久久久久久久| 国产 精品1| 国产xxxxx性猛交| 久久国产精品大桥未久av| 91国产中文字幕| 亚洲欧洲精品一区二区精品久久久 | 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 又大又爽又粗| 91老司机精品| 久久人人爽人人片av| 999久久久国产精品视频| 亚洲国产精品999| 视频在线观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利视频在线观看免费| 免费黄频网站在线观看国产| 成年美女黄网站色视频大全免费| av片东京热男人的天堂| 欧美人与性动交α欧美软件| 日韩 欧美 亚洲 中文字幕| 成人亚洲欧美一区二区av| 哪个播放器可以免费观看大片| 久久久久久久久免费视频了| 中文字幕制服av| 黄色一级大片看看| 99热网站在线观看| 在现免费观看毛片| 国产又爽黄色视频| 久久天躁狠狠躁夜夜2o2o | 欧美精品亚洲一区二区| 久久av网站| 久久99精品国语久久久| 久久精品久久精品一区二区三区| 妹子高潮喷水视频| 亚洲中文av在线| 少妇被粗大猛烈的视频| av卡一久久| 日韩制服丝袜自拍偷拍| 国产男女超爽视频在线观看| 一本—道久久a久久精品蜜桃钙片| 人妻人人澡人人爽人人| 国产精品香港三级国产av潘金莲 | 亚洲综合色网址| 欧美国产精品va在线观看不卡| 日本欧美视频一区| 亚洲国产欧美一区二区综合| 欧美久久黑人一区二区| 蜜桃国产av成人99| 老司机影院毛片| 免费不卡黄色视频| 亚洲欧美清纯卡通| 日韩 欧美 亚洲 中文字幕| 欧美精品一区二区大全| 国产精品女同一区二区软件| 三上悠亚av全集在线观看| 欧美激情极品国产一区二区三区| 免费不卡黄色视频| 中文乱码字字幕精品一区二区三区| 国产淫语在线视频| 国产精品偷伦视频观看了| 黄片播放在线免费| 熟妇人妻不卡中文字幕| av.在线天堂| 啦啦啦在线观看免费高清www| 热99国产精品久久久久久7| 亚洲第一av免费看| 十八禁人妻一区二区| 欧美日韩视频高清一区二区三区二| 超碰97精品在线观看| 国产探花极品一区二区| 99久久精品国产亚洲精品| 我要看黄色一级片免费的| 一级爰片在线观看| 亚洲av日韩在线播放| 成人国产av品久久久| 日韩一卡2卡3卡4卡2021年| 国产精品无大码| 日韩制服骚丝袜av| svipshipincom国产片| 久久国产亚洲av麻豆专区| 下体分泌物呈黄色| 久久久久久人妻| 成年动漫av网址| 国产精品.久久久| 中文天堂在线官网| 日本wwww免费看| 久久久精品94久久精品| 国产黄色视频一区二区在线观看| 亚洲国产看品久久| 国产一区二区激情短视频 | 青草久久国产| av卡一久久| 免费黄频网站在线观看国产| 亚洲精品自拍成人| av在线播放精品| 日日啪夜夜爽| 欧美精品一区二区大全| 一区二区三区四区激情视频| videos熟女内射| 国产女主播在线喷水免费视频网站| 国产精品欧美亚洲77777| 一级黄片播放器| 精品少妇黑人巨大在线播放| 国产精品亚洲av一区麻豆 | 日本色播在线视频| 欧美日韩国产mv在线观看视频| 搡老岳熟女国产| a级毛片在线看网站| 91老司机精品| 伊人久久国产一区二区| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 男女午夜视频在线观看| 汤姆久久久久久久影院中文字幕| 中文字幕高清在线视频| 国产精品久久久av美女十八| 这个男人来自地球电影免费观看 | 你懂的网址亚洲精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 免费黄色在线免费观看| 女性生殖器流出的白浆| 丰满乱子伦码专区| 黑人巨大精品欧美一区二区蜜桃| 日韩av不卡免费在线播放| 欧美久久黑人一区二区| 国产野战对白在线观看| 久久性视频一级片| 亚洲欧美一区二区三区久久| 精品一品国产午夜福利视频| 99精品久久久久人妻精品| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精| 日本一区二区免费在线视频| 成人18禁高潮啪啪吃奶动态图| 一二三四在线观看免费中文在| 男人添女人高潮全过程视频| 国产伦人伦偷精品视频| 可以免费在线观看a视频的电影网站 | 成年美女黄网站色视频大全免费| 91老司机精品| 国产日韩欧美在线精品| 别揉我奶头~嗯~啊~动态视频 | 啦啦啦中文免费视频观看日本| av在线app专区| 18禁裸乳无遮挡动漫免费视频| 欧美精品亚洲一区二区| 午夜福利在线免费观看网站| 美女扒开内裤让男人捅视频| 超碰97精品在线观看| 老司机影院毛片| 国产精品免费大片| 免费看av在线观看网站| 亚洲人成77777在线视频| 成年女人毛片免费观看观看9 | 亚洲精品国产区一区二| 亚洲av日韩在线播放| 青青草视频在线视频观看| 日日啪夜夜爽| 国产精品一国产av| av女优亚洲男人天堂| 亚洲男人天堂网一区| 久久这里只有精品19| 国产 精品1| 九草在线视频观看| 男女边吃奶边做爰视频| 高清在线视频一区二区三区| 久久性视频一级片| 老司机影院成人| 赤兔流量卡办理| 精品国产乱码久久久久久小说| 不卡av一区二区三区| 亚洲成av片中文字幕在线观看| e午夜精品久久久久久久| 欧美精品亚洲一区二区| 久久精品人人爽人人爽视色| 国产高清国产精品国产三级| 久久久国产欧美日韩av| 日韩av在线免费看完整版不卡| 精品少妇久久久久久888优播| 精品亚洲成a人片在线观看| 国产精品 国内视频| 校园人妻丝袜中文字幕| 亚洲少妇的诱惑av| 亚洲精品自拍成人| 无限看片的www在线观看| 国产又色又爽无遮挡免| 日日撸夜夜添| 国产一级毛片在线| 国产免费福利视频在线观看| 我的亚洲天堂| 日韩一区二区视频免费看| 国产 一区精品| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 又粗又硬又长又爽又黄的视频| 亚洲欧美激情在线| av.在线天堂| 国产人伦9x9x在线观看| 精品亚洲成a人片在线观看| 美女大奶头黄色视频| 亚洲国产精品一区三区| 久久国产精品男人的天堂亚洲| 国产精品 国内视频| 丁香六月欧美| 欧美变态另类bdsm刘玥| 一级毛片 在线播放| 曰老女人黄片| 亚洲成人av在线免费| 亚洲精品国产色婷婷电影| 男女午夜视频在线观看| 午夜av观看不卡| 女性被躁到高潮视频| 热99久久久久精品小说推荐| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲高清精品| 午夜日本视频在线| 夜夜骑夜夜射夜夜干| 国产日韩欧美视频二区| 人妻一区二区av| 久久久国产一区二区| 激情视频va一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 嫩草影院入口| 18禁国产床啪视频网站| 中文字幕色久视频| 悠悠久久av| 1024视频免费在线观看| 亚洲精品国产av成人精品| 人人妻人人爽人人添夜夜欢视频| 日韩av免费高清视频| 九九爱精品视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区国产| 亚洲欧美精品自产自拍| 丁香六月欧美| 啦啦啦 在线观看视频| www.熟女人妻精品国产| 精品视频人人做人人爽| 成人国产av品久久久| 成人黄色视频免费在线看| 午夜激情av网站| 国产精品av久久久久免费| av国产精品久久久久影院| 热99久久久久精品小说推荐| 国产精品欧美亚洲77777| 极品少妇高潮喷水抽搐| 精品国产一区二区久久| 悠悠久久av| 汤姆久久久久久久影院中文字幕| 久久久久国产精品人妻一区二区| 女人久久www免费人成看片| av天堂久久9| 丝袜喷水一区| 色吧在线观看| 人人妻,人人澡人人爽秒播 | 亚洲免费av在线视频| 一级a爱视频在线免费观看| 你懂的网址亚洲精品在线观看| 精品一区在线观看国产| 日本91视频免费播放| 亚洲七黄色美女视频| 久久青草综合色| 国产精品国产av在线观看| 少妇被粗大猛烈的视频| 激情视频va一区二区三区| 亚洲av电影在线进入| 大码成人一级视频| tube8黄色片| 国产黄色免费在线视频| 9热在线视频观看99| 捣出白浆h1v1| 精品午夜福利在线看| av在线app专区| 亚洲欧美一区二区三区黑人| 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 亚洲精品日韩在线中文字幕| 丁香六月欧美| 亚洲在久久综合| 国产麻豆69| 老汉色∧v一级毛片| 欧美激情高清一区二区三区 | 国产精品香港三级国产av潘金莲 | 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡| 亚洲激情五月婷婷啪啪| 97精品久久久久久久久久精品| 在线看a的网站| 高清不卡的av网站| 美女高潮到喷水免费观看| 久久久国产一区二区| 亚洲国产毛片av蜜桃av| 欧美人与性动交α欧美精品济南到| 一区二区三区激情视频| 国产免费福利视频在线观看| 在线 av 中文字幕| 你懂的网址亚洲精品在线观看| 久久热在线av| 亚洲情色 制服丝袜| 久久精品aⅴ一区二区三区四区| 日韩精品有码人妻一区| 男人舔女人的私密视频| 天天操日日干夜夜撸| 一本大道久久a久久精品| 最近中文字幕2019免费版| 中文字幕最新亚洲高清| 久久av网站| 久久久久视频综合| 国产亚洲av片在线观看秒播厂| 一级片免费观看大全| 欧美国产精品va在线观看不卡| 久久鲁丝午夜福利片| 精品第一国产精品| 精品视频人人做人人爽| 丰满乱子伦码专区| 人妻 亚洲 视频| 精品国产一区二区久久| 精品人妻一区二区三区麻豆| 如日韩欧美国产精品一区二区三区| 日韩中文字幕欧美一区二区 | 亚洲美女黄色视频免费看| 国产精品久久久久成人av| 成人亚洲欧美一区二区av| 国产精品女同一区二区软件| 一本一本久久a久久精品综合妖精| 最近最新中文字幕免费大全7| 国产在线一区二区三区精| 欧美日韩福利视频一区二区| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 水蜜桃什么品种好| 亚洲av中文av极速乱| 91精品国产国语对白视频| 考比视频在线观看| 伦理电影免费视频| 久久国产精品男人的天堂亚洲| 精品国产乱码久久久久久小说| 亚洲精品第二区| 男女免费视频国产| 国产成人精品久久二区二区91 | 国产成人91sexporn| 国精品久久久久久国模美| 亚洲国产欧美网| 九九爱精品视频在线观看| 国产一区二区三区综合在线观看| 亚洲熟女精品中文字幕| 美国免费a级毛片| 亚洲第一青青草原| bbb黄色大片| 国产精品一二三区在线看| 亚洲图色成人| 精品国产乱码久久久久久小说| 下体分泌物呈黄色| 国产欧美亚洲国产| 综合色丁香网| 精品亚洲成国产av| 久久婷婷青草| 视频在线观看一区二区三区| 久久热在线av| 国产淫语在线视频| 国产精品成人在线| 欧美日韩视频高清一区二区三区二| 一边摸一边做爽爽视频免费| 久久av网站| 日韩精品有码人妻一区| 最近的中文字幕免费完整| 一本一本久久a久久精品综合妖精| 久久99一区二区三区| 国产1区2区3区精品| 少妇精品久久久久久久| 男人操女人黄网站| 国产免费一区二区三区四区乱码| 久热这里只有精品99| 成年女人毛片免费观看观看9 | 又粗又硬又长又爽又黄的视频| 日韩电影二区| 亚洲精品乱久久久久久| 秋霞伦理黄片| 国产男人的电影天堂91| 18禁国产床啪视频网站| √禁漫天堂资源中文www| 另类亚洲欧美激情| 我的亚洲天堂| 老司机影院毛片| xxxhd国产人妻xxx| 99精国产麻豆久久婷婷| 国产99久久九九免费精品| 欧美黑人欧美精品刺激| 久久韩国三级中文字幕| 老汉色av国产亚洲站长工具| 丰满少妇做爰视频| 天天添夜夜摸| 亚洲精品日本国产第一区| 国产男女超爽视频在线观看| 在现免费观看毛片| 成年美女黄网站色视频大全免费| 亚洲欧美日韩另类电影网站| 欧美日韩视频精品一区| 亚洲av成人不卡在线观看播放网 | 一本一本久久a久久精品综合妖精| 亚洲av福利一区| 女的被弄到高潮叫床怎么办| 欧美精品一区二区免费开放| 黄网站色视频无遮挡免费观看| 69精品国产乱码久久久| 欧美日韩精品网址| 国产又爽黄色视频| 在线观看免费午夜福利视频| 又黄又粗又硬又大视频| 午夜免费观看性视频| 哪个播放器可以免费观看大片| 九色亚洲精品在线播放| av有码第一页| 成人黄色视频免费在线看| 免费看av在线观看网站| 午夜福利在线免费观看网站| 国产亚洲午夜精品一区二区久久| 精品国产超薄肉色丝袜足j| 深夜精品福利| 免费不卡黄色视频| 天天躁夜夜躁狠狠躁躁| 国产极品粉嫩免费观看在线| 91成人精品电影| 午夜老司机福利片| 黄色 视频免费看| 日韩制服骚丝袜av| 中文字幕制服av| 久久99一区二区三区| 在现免费观看毛片| 国产成人精品无人区| 亚洲免费av在线视频| 青春草国产在线视频| 免费在线观看黄色视频的| 久久久久久人人人人人| 欧美精品一区二区免费开放| 久久久久久久国产电影| 亚洲天堂av无毛| 国产高清国产精品国产三级| 伦理电影大哥的女人| 国精品久久久久久国模美| 欧美激情高清一区二区三区 | av卡一久久| 亚洲国产日韩一区二区| 亚洲国产精品一区三区| 亚洲国产精品999| √禁漫天堂资源中文www| 亚洲精品美女久久久久99蜜臀 | av卡一久久| 日本黄色日本黄色录像| 亚洲四区av| 欧美激情高清一区二区三区 | 免费观看人在逋| 咕卡用的链子| 亚洲精品视频女| 国产精品一区二区精品视频观看| 交换朋友夫妻互换小说| 18禁动态无遮挡网站| 在线观看免费午夜福利视频| 精品国产一区二区三区四区第35| 少妇的丰满在线观看| 18禁国产床啪视频网站| 亚洲婷婷狠狠爱综合网| 精品国产露脸久久av麻豆| 日韩 亚洲 欧美在线| 亚洲国产日韩一区二区| 中文字幕最新亚洲高清| 国产精品二区激情视频| 天堂8中文在线网| 国产成人系列免费观看| 丝袜美腿诱惑在线| 亚洲人成77777在线视频| 美女扒开内裤让男人捅视频| 亚洲少妇的诱惑av| 国产黄频视频在线观看| 9191精品国产免费久久| 成年动漫av网址| 操美女的视频在线观看| 自线自在国产av| 国产极品粉嫩免费观看在线| 男人舔女人的私密视频| 亚洲欧洲国产日韩| 亚洲图色成人| 日本vs欧美在线观看视频| 老司机深夜福利视频在线观看 | 日韩精品免费视频一区二区三区| 在线观看三级黄色| 悠悠久久av| 纵有疾风起免费观看全集完整版| av电影中文网址| 亚洲第一青青草原| 亚洲一级一片aⅴ在线观看| 亚洲成人一二三区av| 国产成人精品无人区| 日日摸夜夜添夜夜爱| 女人久久www免费人成看片| 久久久国产一区二区| 国产深夜福利视频在线观看| 极品人妻少妇av视频| 国产无遮挡羞羞视频在线观看| 街头女战士在线观看网站| 国产精品久久久久成人av| 亚洲人成网站在线观看播放| 精品国产乱码久久久久久小说| 女性被躁到高潮视频| 美女大奶头黄色视频| 国产日韩欧美视频二区| 美女主播在线视频| 精品一品国产午夜福利视频| 国产成人a∨麻豆精品| 久久久久久久久久久久大奶| 黑人巨大精品欧美一区二区蜜桃| 男女高潮啪啪啪动态图| 久久亚洲国产成人精品v| 黄色毛片三级朝国网站| 成人亚洲精品一区在线观看| 97人妻天天添夜夜摸| 久久久亚洲精品成人影院| 国产一级毛片在线| 亚洲精品日韩在线中文字幕| 老司机亚洲免费影院| 天天操日日干夜夜撸| 妹子高潮喷水视频| 丝袜人妻中文字幕| 人人妻,人人澡人人爽秒播 | 久久久久久久久久久久大奶| 国产老妇伦熟女老妇高清| 久久免费观看电影| 日韩中文字幕欧美一区二区 | 最近中文字幕2019免费版| 免费黄色在线免费观看| 欧美日韩视频精品一区| 精品视频人人做人人爽| 久久久精品94久久精品| 叶爱在线成人免费视频播放| 老司机影院成人| 日韩精品有码人妻一区| 亚洲欧美色中文字幕在线| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产欧美日韩在线播放| 免费不卡黄色视频| 91精品三级在线观看| 一二三四中文在线观看免费高清| 国产精品 国内视频| 国产精品国产三级国产专区5o| 七月丁香在线播放| 日韩电影二区| 一级毛片电影观看| 老汉色av国产亚洲站长工具| 国产精品一区二区精品视频观看| 国产一区二区三区av在线| 欧美97在线视频| 国产精品欧美亚洲77777| 看十八女毛片水多多多| 在现免费观看毛片| 久久久久精品久久久久真实原创| 巨乳人妻的诱惑在线观看| 国产精品人妻久久久影院| 五月开心婷婷网| 久久精品国产综合久久久| 大香蕉久久成人网| 日韩精品免费视频一区二区三区| 亚洲情色 制服丝袜| 成人毛片60女人毛片免费| 亚洲伊人色综图| 啦啦啦在线观看免费高清www| 色94色欧美一区二区| 国产精品熟女久久久久浪| 哪个播放器可以免费观看大片| 午夜福利一区二区在线看| 我的亚洲天堂| 在线天堂最新版资源| 久久久久精品性色| 哪个播放器可以免费观看大片| 嫩草影视91久久| 人人妻人人澡人人看| 男人添女人高潮全过程视频| 我的亚洲天堂| 亚洲av电影在线进入| 美女视频免费永久观看网站| 免费av中文字幕在线|