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

    北太平洋副熱帶西部模態(tài)水年代際變化特征及其機(jī)制分析*

    2021-10-11 09:53:30代秋亭劉子洲翟方國(guó)顧艷鎮(zhèn)李元杰
    海洋與湖沼 2021年5期
    關(guān)鍵詞:海表預(yù)先水體

    代秋亭 劉子洲 劉 聰 翟方國(guó) 顧艷鎮(zhèn) 李元杰

    (1. 中國(guó)海洋大學(xué)海洋與大氣學(xué)院 青島 266100; 2. 海南浙江大學(xué)研究院 三亞 572000;3. 中國(guó)人民解放軍75839 部隊(duì) 廣州 510080)

    Masuzawa (1969)發(fā)現(xiàn)并命名了北太平洋副熱帶西部模態(tài)水(subtropical mode water, STMW), 指出STMW 位于北太平洋副熱帶環(huán)流西北部溫躍層內(nèi), 溫度范圍為 16—19 °C, 具有性質(zhì)均一、位勢(shì)渦度(potential vorticity, PV)低等特點(diǎn), 主要形成于黑潮及黑潮延伸體(Kuroshio Extension, KE)南側(cè)深混合層區(qū)。

    從11 月至次年3 月, KE 流域受強(qiáng)東亞季風(fēng)的控制, 巨大的海氣溫差和風(fēng)應(yīng)力在海洋上層引起強(qiáng)烈的垂向混合, 使得該海域南側(cè)混合層快速發(fā)展。至3月初, 混合層深度(mixing layer depth, MLD)達(dá)到最深(Hanawa, 1987; Sugaet al, 1990; Bingham, 1992;Croninet al, 2013)。隨后, 季節(jié)性溫躍層形成并將潛沉的低位勢(shì)渦度水體與大氣環(huán)境相隔離(Sugaet al,1990, 1995)。潛沉后形成的STMW, 被風(fēng)生副熱帶環(huán)流輸送至臺(tái)灣島東側(cè)黑潮流域以及副熱帶逆流(subtropical countercurrent, STCC)北側(cè)(Kubokawaet al, 1999; Liuet al, 2007; Yanet al, 2013; Yuet al,2015), 改變當(dāng)?shù)卮伪韺訙佧}結(jié)構(gòu)。部分STMW 隨黑潮再次回到生成區(qū)通過(guò)浮露進(jìn)入冬季混合層, 改變生成區(qū)冬季混合層水體溫度。因此, STMW 再通風(fēng)也被認(rèn)為是海洋對(duì)太平洋年代際濤動(dòng)(Pacific decadal oscillation, PDO)(Trenberthet al, 1994; Mantuaet al,1997)的負(fù)反饋機(jī)制之一(Liuet al, 2007; Newmanet al, 2016)。此外, 部分STMW 在STCC 北側(cè)堆棧形成低PV 池, 引起溫躍層抬升, 利于東向STCC 的形成和維持(Kubokawaet al, 1999; Yamanakaet al, 2008;Kobashiet al, 2012)。又因, 具有晚冬混合層水體性質(zhì)的模態(tài)水富含溶解氧、可溶性有機(jī)物、無(wú)機(jī)碳、營(yíng)養(yǎng)鹽等, STMW 分布對(duì)海洋生態(tài)系統(tǒng)亦有重要意義(Sukigaraet al, 2009, 2014)。

    綜上, STMW 能夠常年保留生成區(qū)海氣相互作用的信號(hào), 影響北太平洋副熱帶環(huán)流以及氣候變化, 同時(shí)也影響區(qū)域生態(tài)系統(tǒng)。因此, 其年代際變化機(jī)制一直是近幾十年的研究重點(diǎn)。而引起STMW 體積年代際變化的機(jī)制主要有兩種, 即海氣相互作用和海洋動(dòng)力過(guò)程。早期研究普遍認(rèn)為海氣相互作用是調(diào)制STMW 溫度特性以及體積變化的主要機(jī)制(Bingham,1992; Sugaet al, 1995; Yasudaet al, 1999)。而隨后,Qiu 等(2006)利 用 1993—2005 年 期 間 的 CTD(conductivity-temperature-depth) 和 XBT(expendable bathythermograph)數(shù)據(jù)研究了冬季海表冷卻以及KE系統(tǒng)動(dòng)力狀態(tài)對(duì)STMW 形成的影響, 指出KE 動(dòng)力狀態(tài)是STMW 年代際變化的主要影響因素。當(dāng)KE 處于動(dòng)力不穩(wěn)定狀態(tài)時(shí), 氣旋式渦旋活動(dòng)加強(qiáng), 將大量高位勢(shì)渦度性質(zhì)的KE 水體注入副熱帶環(huán)流, 使得海洋層結(jié)加強(qiáng), 不利于深混合層發(fā)展, 由此造成STMW厚度減小。另外, 有不少研究表明KE 動(dòng)力狀態(tài)可能與PDO 誘導(dǎo)的中東太平洋的風(fēng)應(yīng)力旋度異常密切相關(guān): PDO 正相位期間, 中東太平洋海表有正的風(fēng)應(yīng)力旋度異常, 產(chǎn)生負(fù)的海表面高度異常(sea surface height anomaly, SSHA); 這種SSHA 異常信號(hào)將以第一斜壓Rossby 波的速度向西傳輸, 約3—4 a 后抵至KE 流域, 引起KE 系統(tǒng)向不穩(wěn)定狀態(tài)轉(zhuǎn)變(Seageret al, 2001; Schneideret al, 2002; Qiu, 2003; Taguchiet al,2007; Sasakiet al, 2013; Qiuet al, 2014)。此外, Yu 等(2015)分析62 a 的模式模擬數(shù)據(jù)后發(fā)現(xiàn)STMW 生成區(qū)的潛沉率及次表層水體位勢(shì)渦度也都呈現(xiàn)出明顯的PDO 信號(hào)。以上可知STMW 年代際變化實(shí)際上能與PDO 有所聯(lián)系, 且KE 是兩者相聯(lián)系的重要紐帶。同時(shí), 近期 Kim 等(2020)通過(guò)運(yùn)行海洋環(huán)流模式(Ocean General Circulation Model)對(duì)STMW 在1960—2009 年的年代際變化機(jī)制的診斷結(jié)論明確指出,STMW 年代際變化的主要驅(qū)動(dòng)因素在20 世紀(jì)80 年代后期發(fā)生轉(zhuǎn)變: 由海氣相互作用主導(dǎo)轉(zhuǎn)變?yōu)楹Q髣?dòng)力過(guò)程主導(dǎo), 而這主要?dú)w因于西太平洋和中東太平洋的大氣強(qiáng)迫條件的變化。

    綜上可知, 近40 a STMW 年代際變化主要受海洋動(dòng)力過(guò)程的影響, 且主要由海洋預(yù)先層結(jié)調(diào)控。但就目前而言, 對(duì)STMW 年代際變化的具體調(diào)節(jié)方式還尚不明確。為解決這一問(wèn)題, 本文對(duì)混合層形成體積進(jìn)行了量化分析, 指出KE 上游南側(cè)區(qū)域的海洋上層預(yù)先層結(jié)年代際變化(3.1)可以通過(guò)調(diào)節(jié)混合層底部的卷吸作用使冬季混合層體積產(chǎn)生顯著年代際變化(3.2), 以此影響春季潛沉水體體積, 調(diào)控STMW形成體積年代際變化。因此, 本研究通過(guò)混合層收支量化分析將預(yù)先層結(jié)對(duì)STMW 的調(diào)節(jié)過(guò)程用卷吸和潛沉串聯(lián)起來(lái), 更深入地剖析了STMW 年代際變化機(jī)制。另外, 本文指出STMW 年代際變化實(shí)際上受PDO 的遠(yuǎn)場(chǎng)調(diào)控(3.3)。

    1 數(shù)據(jù)和方法

    1.1 數(shù)據(jù)

    為了探究STMW 形成的年代際變化特征及機(jī)制,本文使用來(lái)自大洋環(huán)流模式ECCO2 (the estimating the circulation and climate of the ocean, Phase II project)中Cube 92 的再分析數(shù)據(jù), 數(shù)據(jù)時(shí)間選擇了1992 年 1 月—2019 年 3 月, 水平分辨率為 1/4°(~27 km), 垂向上共有50 層, 最深深度達(dá)6 135 m。ECCO2 通過(guò)使用格林函數(shù)結(jié)合多種觀測(cè)數(shù)據(jù)對(duì)環(huán)流模型進(jìn)行同化, 擁有精細(xì)的分辨率, 能再現(xiàn)全球范圍的海洋。數(shù)據(jù)約束包了高度計(jì)資料中的海表平面高度異常、平均海平面、全球海洋數(shù)據(jù)同化實(shí)驗(yàn)高分辨率海表溫度試點(diǎn)項(xiàng)目(global ocean data assimilation experiment high resolution sea surface temperature pilot project)的海表面溫度數(shù)據(jù)以及來(lái)自海洋環(huán)流實(shí)驗(yàn)(world ocean circulation experiment)、熱帶海洋大氣(tropical atmosphere ocean)、實(shí)時(shí)地轉(zhuǎn)海洋學(xué)陣列(array for real-time geostrophic oceanography, ARGO)、拋棄式深水溫度計(jì)(expendable bathythermograph)等的溫鹽剖面數(shù)據(jù); 控制參數(shù)包括初始溫度和鹽度, 大氣表層邊界條件, Large 等(1994) KPP(K-profile parameterization)方案中的背景場(chǎng)垂向擴(kuò)散及理查森數(shù), 海氣/海冰/氣冰拖曳系數(shù), 冰/海洋/雪反照率, 底部阻力和垂直黏度。詳細(xì)信息可以參閱Menemenlis等(2008)和Fox-Kemper 等(2008)。另外, 本文所使用的ECCO2 數(shù)據(jù)中, 海表熱通量及蒸發(fā)由Large 等(1982)的塊體公式計(jì)算得出; 海表風(fēng)應(yīng)力的計(jì)算使用了E. Vera 的阻力系數(shù)參數(shù)(Largeet al,1995); 降水?dāng)?shù)據(jù)來(lái)自全球降水氣候?qū)W計(jì)劃(global precipitation climatology project)。

    圖1 155°E 剖面處ECCO2 和EN4 氣候態(tài)平均數(shù)據(jù)對(duì)比Fig.1 Comparison between ECCO2 and EN4 climatological data at 155°E section

    1.2 方法

    在副熱帶環(huán)流西北有一個(gè)深MLD 池, 其向東南遞減, 并與南部淺MLD 區(qū)域之間存在狹窄的深度遞減帶, 即MLD 鋒; 低PV 水形成于MLD 鋒和露頭線的交界處, 密度越高的低PV 水形成位置越偏向東北(Kubokawaet al, 1999)。因此, 探討混合層年代際變化機(jī)制對(duì)研究STMW 年代際變化十分重要。在這里,MLD 被定義為密度較10 m 處密度大0.1 kg/m3水體所在深度(Oka, 2009; Guoet al, 2018); 混合層體積守恒公式如下(Nishikawaet al, 2013; Guoet al, 2018):

    方程左邊為混合層形成體積, 是時(shí)間1t至t2時(shí)刻控制體的體積變化量; 右邊各項(xiàng)分別為海氣形成體積、平流輸運(yùn)體積、卷吸體積以及表征越密混合體積的殘余項(xiàng)(R)。MLD 具有明顯的季節(jié)性變化特征, 但只有晚冬的MLD 信號(hào)能夠保留在STMW 之中(Stommel,1979)。因此這里本文選擇MLD 加深時(shí)期進(jìn)行時(shí)間積分。關(guān)于積分時(shí)間的選擇在3.2 節(jié)再詳細(xì)說(shuō)明。

    本次研究將STMW 定義為黑潮南側(cè), 混合層之下,PV≤2×10-10m-1s-1且θσ?[24.6, 25.6] kg/m3的低PV 水體。另外, KE 主軸由0.4 m 的海表高度等值線表征。這里對(duì)STMW 的密度范圍相對(duì)于以往研究更大, 主要是考慮到近年來(lái)KE 南側(cè)表層水體的溫度持續(xù)上升且該區(qū)域鹽度從2008 年開(kāi)始顯著降低, 引起KE 南側(cè)冬季混合層水體以及模態(tài)水密度發(fā)生顯著變化(Sugimotoet al, 2013, 2017)。定義STMW 密度范圍較小時(shí), 會(huì)使得較輕的低位勢(shì)渦度水體被去除, 造成模態(tài)水體積減小的虛假信號(hào)。

    在黑潮南側(cè)MLD 具有顯著的季節(jié)變化, 其在冬季迅速加深, 至3 月達(dá)到最深, 在這一階段MLD 的持續(xù)加深為STMW 形成提供重要的物質(zhì)準(zhǔn)備。隨后,季節(jié)性溫躍層逐漸建立, MLD 迅速變淺, 部分混合層水體潛沉進(jìn)入到溫躍層中形成STMW(Rainvilleet al,2007; Oka, 2009; Daviset al, 2011)。至5 月, 季節(jié)性溫躍層形成, 潛沉結(jié)束。因此本研究將3—5 月劃分為STMW 形成的時(shí)期, 這一時(shí)期STMW 體積在局部迅速增大, 侵蝕和輸運(yùn)對(duì)STMW 體積影響可以忽略,STMW 的形成體積Vf、形成厚度Hf可以5 月及3 月間STMW 體積和厚度差異表征:

    這里的下標(biāo)‘f’表示形成, ‘5’和‘3’分別表示5 月和3 月。STMW 的形成面積則以STMW 形成厚度大于10 m 的區(qū)域面積表征。另外, 本文選擇3—5月(此處以下標(biāo)‘345’表示)對(duì)潛沉進(jìn)入到溫躍層的水體體積進(jìn)行積分, 以得到年潛沉體積:

    式中, 345 表示3、4、5 月。KE 南側(cè), 根據(jù)STMW 定義,θσ?[24.6, 25.6] kg/m3的混合層水體年潛沉體積即為STMW 形成體積。

    2 STMW 體積年代際變化

    為了方便后文對(duì)STMW 形成年代際變化機(jī)制的分析, 在這一部分, 我們?cè)谶@里引入STMW 生成區(qū)(133—170°E、30—36°N; 圖3 中藍(lán)色矩形框框出), 該區(qū)域是具有STMW 密度性質(zhì)水體的3 月露頭區(qū), 也是MLD 以及低PV 水體形成厚度年變化最顯著區(qū)域。其中低 PV 水形成厚度的年變化標(biāo)準(zhǔn)差普遍高于170 m(圖3a), MLD 的年變化相對(duì)較小(圖3b)。而無(wú)論是低PV 水形成厚度還是MLD, 它們的年變化信號(hào)都集中于生成區(qū)北部近KE 處, 特別是KE 向北凸起處南側(cè)。其中在(140°E, 32°N)、(144°E, 34°N)、(150°E,34°N)附近存在MLD 年較差中心。對(duì)應(yīng)的低PV 水體厚度在這3 個(gè)中心附近都較大。東部的兩個(gè)中心對(duì)應(yīng)KE 主軸向北凸起南側(cè), 正好是SSH 抬升的區(qū)域。由圖3 可發(fā)現(xiàn), 顯著的低PV 水形成厚度年變化出現(xiàn)在偏西位置, 本文將其中具體的原因在3.1 機(jī)制分析中進(jìn)行了詳細(xì)的探討。

    圖2 混合層體積收支控制體Fig.2 The control body for mixing layer volume budget

    圖3 低PV 水體形成厚度及3 月MLD 的年變化標(biāo)準(zhǔn)差場(chǎng)Fig.3 Multi-year standard deviation field of low PV water formation thickness and March MLD

    繪制3 月和5 月STMW 的體積, STMW 形成體積、形成厚度、形成面積, 3 月生成區(qū)混合層體積以及其區(qū)域面積及平均如圖4 所示。圖4a 為5 月(紅線)和3 月(藍(lán)線)STMW 體積的年變化曲線, 由此可知,由ECCO2 和EN4 數(shù)據(jù)得到的STMW 體積都存在較為同步且顯著的年代際變化, 且振幅相當(dāng)。往年形成水體經(jīng)過(guò)侵蝕后殘存的STMW 信號(hào)可以在3 月體現(xiàn)出來(lái), 隨后與新形成的STMW(圖4b 黑色實(shí)線)疊加使5 月STMW 體積年代際變化進(jìn)一步放大。在這里需要強(qiáng)調(diào)的是, 3 月STMW 雖存在年代際變化, 但其峰值落后于5 月STMW 以及STMW 形成體積1—2 a的時(shí)間, 顯然3 月STMW 體積只是往年殘存年代際信號(hào)的體現(xiàn), 并不能作為促使STMW 體積產(chǎn)生年代際變化的決定性因素, 唯有STMW 形成體積可決定STMW 體積的年代際變化。

    圖4b 中, STMW 生成體積的年變化特征(黑色加粗曲線), 在年至年代際上顯然都與3 月生成區(qū)混合層水體體積(黑色加粗虛線)相吻合, 也呈現(xiàn)出顯著的年代際變化特征以及整體減小趨勢(shì), 振幅為3×1014—4×1014m3, 其中混合層體積的年代際變化由MLD 的年代際變化決定, 而STMW 形成體積的年代際變化同時(shí)體現(xiàn)在其形成面積和形成厚度上。STMW形成及混合層體積的遞減速率分別為0.057×1014m3/a和0.023×1014m3/a。其中混合層體積的遞減應(yīng)與近年觀測(cè)到的海表升溫有關(guān)(Sugimotoet al, 2017; Yuet al,2020)。

    圖4 3 月和5 月STMW 體積、STMW 形成及3 月混合層的時(shí)間序列Fig.4 Time series of March and May STMW volume, STMW formation, and March mixed layer

    將STMW 生成體積和生成區(qū)混合層體積的整體變化趨勢(shì)剔除, 我們發(fā)現(xiàn)STMW 形成體積及混合層體積呈現(xiàn)明顯正、負(fù)異常(如圖4b 陰影區(qū)和空白區(qū)):正異常分別在 1992—1997 年、2000—2005 年和2011—2017 年; 負(fù)異常對(duì)應(yīng)在 1998—1999 年和2006—2010 年。盡管在正、負(fù)異常期間MLD 差較STMW 形成厚度差明顯更大(圖4c), 但在年代際變化及生成區(qū)MLD 體積與STMW 形成體積在空間分布上十分吻合(圖5)。首先, 對(duì)于3 月生成區(qū)MLD 來(lái)說(shuō), 其在STMW 形成體積正異常期間, 普遍大于240 m, 最大MLD 可超過(guò)300 m(圖5b); 而STMW 形成體積負(fù)異常期間, MLD 厚度相對(duì)減小了60 m 左右, 最大深度只能達(dá)到240 m(圖5d)。此外, 正、負(fù)異常期間的MLD 差呈現(xiàn)出東南-西北向的條狀帶, 在這些條狀帶上MLD 差取顯著正值, 對(duì)應(yīng)其上正SSHa 由外向內(nèi)迅速增加(圖5f)。再對(duì)比STMW 體積正、負(fù)異常期間, KE 主軸的彎曲狀態(tài), 便可以發(fā)現(xiàn)正相位時(shí)期KE主軸相對(duì)有規(guī)律的南北彎曲能夠使其向北凸起處南側(cè)的MLD 加深及SSH 升高, 而向南凹陷處南側(cè)的MLD 及SSH 有弱減小。

    而在STMW 形成體積正異常期間, STMW 形成厚度在黑潮及KE 上游南側(cè)區(qū)域增加顯著, 特別是在138—144°E 之間的區(qū)域(圖5e)。此外, MLD 在西部有更為明顯的加深, 可知生成區(qū)西部是STMW 形成厚度及MLD 年代際變化產(chǎn)生的關(guān)鍵區(qū)域。

    圖5 STMW 形成厚度(a, c)、3 月MLD(b, d)以及它們?cè)赟TMW 形成體積正、負(fù)異常時(shí)期的差(e, f)Fig.5 Multi-year mean field of STMW formation thickness (a, c), March MLD (b, d), and their difference in positive and negative phases of STMW formation volume

    以上的內(nèi)容, 主要探討了正負(fù)相位時(shí)期STMW形成厚度及MLD 的分布差異。但STMW 形成體積的年代際變化同時(shí)展現(xiàn)在形成厚度和面積上, 于此本文進(jìn)一步分析了STMW 形成區(qū)域的緯度分布特征。以155°E 剖面為例, 5 月和3 月PV、STMW 密度邊界、MLD 在STMW 體積正、負(fù)異常期間多年平均結(jié)果如圖6 所示。圖中STMW 密度邊界(圖6 中兩條白色曲線)由5 月STMW 在剖面處的密度性質(zhì)決定, 以準(zhǔn)確地再現(xiàn)不同時(shí)期STMW 所在密度層。雖然這里的密度邊界于STMW 形成體積正、負(fù)異常期間取不同值,但可以準(zhǔn)確定位3 月STMW 露頭區(qū)(圖6 中白色曲線之間的水體在海表的位置)。如圖6a 和c 可見(jiàn), STMW體積正、負(fù)異常期間3 月混合層除了有顯著的深度變化之外, 露頭區(qū)南北邊界也向外擴(kuò)張(圖6b 和6d): 在STMW 體積負(fù)異常期間, 露頭區(qū)南邊界在30°N, 北邊界在36°N; 而在正異常期間, 露頭區(qū)南邊界向外擴(kuò)了0.5°, 北邊界擴(kuò)張不明顯。露頭區(qū)的擴(kuò)張也引起了5 月STMW 形成面積增加。如圖6a 和6c 所示,STMW 厚度在28°N 左右有一個(gè)往北迅速遞增的鋒面,可見(jiàn)當(dāng)年形成的STMW 在5 月主要分布在28°N 以北,相對(duì)于露頭區(qū)偏于南部1°—2°。厚度鋒面受3 月露頭區(qū)影響, 在不同時(shí)期有明顯的位置變化。

    圖6 5 月(a, c)和3 月(b, d)155°E 斷面處的PV 分布圖Fig.6 PV distribution at cross section 155°E in May (a, c) and March (b, d)

    取STMW 上邊界為STMW 厚度頂部, STMW 下邊界為STMW 厚度底部, 繪制STMW 生成區(qū)區(qū)域平均的5 月STMW 上下邊界、3 月MLD、5 月STMW厚度、24.6 和25.6 kg/m3所在深度的年變化曲線, 結(jié)果如圖7 所示。可知, 盡管24.6、25.6 kg/m3的深度存在與STMW 形成體積、厚度和面積類(lèi)似的年代際變化, 但STMW 厚年年代際變化并不由其定義的密度邊界位置決定, 主要由下邊界深度變化引起, 而這進(jìn)一步與3 月MLD 深度有關(guān)。

    圖7 生成區(qū)區(qū)域平均的5 月24.6 和25.6 kg/m3 等密度面所在深度、STMW 上下邊界、STMW 厚度及3 月MLD 的年時(shí)間序列Fig.7 Annual time series of the isopycnic faces of 24.6 kg/m3 and 25.6 kg/m3, the upper and lower boundary, and thickness of STMW in May and MLD in March

    3 機(jī)制

    3.1 影響因素

    結(jié)合(2)式可知, 海氣形成率由單位密度上海表浮力通量隨密度的變化梯度表征。而從圖8 中不難發(fā)現(xiàn)B(t)從黑潮及KE 主軸向南北兩側(cè)遞減, 南側(cè)遞減梯度更大, 可知KE 南側(cè)STMW 形成區(qū)有較大的海氣形成率, 為該區(qū)域晚冬深MLD 的形成提供必要條件。

    如圖9a 所示, 海表浮力通量正異常與正SSHa 區(qū)域重疊, 特別是在圖中紅色框圈定的范圍(144°—163°E, 33°—36°N)。STMW 體積正(負(fù))異常期間, KE 主軸即圖9 紅棕色加粗實(shí)(虛)線向北(南)遷移。而在KE 流域, 越往北海表風(fēng)應(yīng)力越強(qiáng)(圖8a 和8b),因此KE 主軸的向北遷移, 不免造成KE 海表蒸發(fā)失熱加劇, 導(dǎo)致海表浮力通量的顯著增加。

    圖8 海表浮力通量、海表渦動(dòng)能(EKE)和預(yù)先層結(jié)(N)在STMW 體積正異常(左)和負(fù)異常(右)時(shí)期的多年平均空間場(chǎng)Fig.8 Multi-year mean spatial fields of surface buoyancy flux (B), surface eddy kinetic energy (EKE), and pre-stratification (N) at positive (left) and negative (right) anomaly of STMW formation volume

    圖9 海表浮力通量(B)、海表渦動(dòng)能(EKE)和預(yù)先層結(jié)(N)在STMW 體積正負(fù)異常時(shí)期間的差異Fig.9 Differences of surface buoyancy flux (B), surface eddy kinetic energy (EKE), and pre-stratification (N) between positive and negative anomalies of STMW volume

    KE 南側(cè)海域海表EKE 是KE 動(dòng)力狀態(tài)的間接反映: KE 不穩(wěn)定時(shí), KE 路徑快速變化, 區(qū)域EKE 成倍增加(Qiuet al, 2007)。如圖8e—8f 所示, EKE 主要沿著KE 主軸分布, 且于140°—152°E KE 主軸向南彎曲處最大, 且140°—152°E 的KE 路徑變化較顯著。STMW 體積正異常時(shí), EKE 在KE 上游南側(cè)(圖9c 藍(lán)框, 132°—147°E, 30°—35°N)表現(xiàn)出極強(qiáng)的負(fù)異常信號(hào), 對(duì)應(yīng)其上也有有正的SSHa 差(圖9c)。

    預(yù)先層結(jié)的分布如圖8g 和8h 所示。由圖可知預(yù)先層結(jié)在KE 及黑潮附近具有與EKE 相似的分布形態(tài)。而在KE 南側(cè)至28°N, 預(yù)先層結(jié)相對(duì)較弱, 這也是STMW 生成區(qū)具有深MLD 的重要原因(Yuet al,2020)。同時(shí), 相對(duì)于STMW 體積負(fù)異常時(shí)期, STMW正異常時(shí)期, 預(yù)先層結(jié)明顯減弱。STMW 體積正、負(fù)異常時(shí)期間的層結(jié)頻率差與EKE 差相似, 同樣在KE上游南側(cè)藍(lán)色框所圈定的區(qū)域(132°—147°E,30°—35°N)內(nèi)具有顯著的負(fù)異常(圖9d)。

    以上, 我們知道在STMW 體積正異常期間, 在KE 上游南側(cè)區(qū)域?qū)?yīng)有EKE 減小與預(yù)先層結(jié)減弱,在KE 下游近KE 主軸處對(duì)應(yīng)海表浮力通量升高。而在它們顯著變化的區(qū)域上, 都有海表面升高與之相對(duì)應(yīng)。其中海表浮力通量產(chǎn)生顯著變化的區(qū)域位于STMW 生成區(qū)偏東部, 其上的MLD 年變化特征并不顯著。因此, 本文將重點(diǎn)放在了KE 上游南側(cè)區(qū)域, 探究EKE 與預(yù)先層結(jié)對(duì)冬季MLD 加深的影響。

    如圖10 所示, 本文分析了STMW 生成區(qū)、KE上游泳南側(cè)及其下游區(qū)域的MLD 與SSH、EKE 與SSH、預(yù)先層結(jié)與EKE 之間的關(guān)系。分析可知, 同一組變量在不同區(qū)域的相關(guān)性差異很明顯。首先對(duì)于MLD 與SSH 來(lái)說(shuō), 兩者在STMW 生成區(qū)、生成區(qū)中的KE 上游以及KE 下游都為正相關(guān), 相關(guān)性最高的區(qū)域在KE 上游, 但是與其他兩個(gè)區(qū)域相差不大。在KE 上游, SSH 每升高0.1 m, MLD 對(duì)應(yīng)加深35 m。然后, 對(duì)于EKE 與SSH 來(lái)說(shuō), 兩者之間表現(xiàn)為負(fù)相關(guān),不同區(qū)域它們之間的相關(guān)性顯著不同, 在KE 上游EKE 與SSH 之間的相關(guān)系數(shù)可達(dá)–0.53, 而在整個(gè)STMW 生成區(qū)EKE 與SSH 相關(guān)性只有–0.17。上游區(qū)域SSH 每上升0.1 m, EKE 對(duì)應(yīng)減小約0.03 m2/s2。

    圖10 STMW 生成區(qū)、KE 上游南側(cè)及其下游區(qū)域的MLD-SSH、EKE-SSH 和N-EKE 相關(guān)性Fig.10 MLD-SSH, EKE-SSH and N-EKE correlations in the STMW formation region, south side of KE upstream and downstream region of KE

    而受KE 狀態(tài)影其向南分離出來(lái)的渦旋活動(dòng), 海洋層結(jié)也發(fā)生相應(yīng)變化(Qiuet al, 2005, 2006, 2007)。如圖10g—10h, EKE 與預(yù)先層結(jié)表現(xiàn)為正相關(guān), 特別是在KE 上游, 兩者的相關(guān)系數(shù)可達(dá)0.68。在該區(qū)域, EKE每增加0.1 m2/s2, 可促使預(yù)先層結(jié)升高0.1×10–3s–1??梢?jiàn), 相對(duì)生成區(qū)中的其它區(qū)域, KE 上游南側(cè)區(qū)域(132°—147°E, 30°—35°N)的MLD 及SSHA 更易受到海洋動(dòng)力狀態(tài)和層結(jié)強(qiáng)度的調(diào)節(jié): 在STMW 體積正(負(fù))異常時(shí)期, KE 上游南側(cè)出現(xiàn)的渦動(dòng)能和預(yù)先層結(jié)強(qiáng)度負(fù)異常(正異常)將進(jìn)一步地調(diào)控當(dāng)?shù)氐腗LD 和SSHA, 促使深厚(淺薄)的混合層水體在晚冬形成。

    3.2 混合層體積收支

    在3.1 節(jié), 分析了受控于海洋動(dòng)力過(guò)程的預(yù)先層結(jié)(Qiuet al, 2007)和能夠反映海氣相互作用的海表浮力通量?jī)烧叩目臻g分布特性, 初步認(rèn)識(shí)到影響MLD 發(fā)展的不同物理過(guò)程具有區(qū)域差異。本節(jié)將進(jìn)一步量化海洋動(dòng)力環(huán)境和海氣相互作用對(duì)混合層加深的影響, 深入探討引起混合層體積年代際變化的主要原因。同時(shí), 前文分析得出的海表浮力通量變化難以對(duì)STMW 生成區(qū)混合層產(chǎn)生顯著影響的結(jié)論,在本節(jié)會(huì)再一次得到驗(yàn)證。

    本節(jié)研究基于混合層體積收支方程(7), 該方程表明, 混合層體積變化主要受4 個(gè)過(guò)程調(diào)節(jié), 它們分別為海氣形成率、沿著等密度面的側(cè)向輸運(yùn)、混合層底部卷吸以及越密混合。各項(xiàng)作用于圖2 所示的控制體, 是密度和時(shí)間的函數(shù)。其中控制體的密度半徑取0.1 kg/m3(θσΔ =0.1 kg/m3)、中心密度為θσ。若以月為單位進(jìn)行時(shí)間積分再取氣候平均可得到各項(xiàng)對(duì)不同密度范圍控制體作用的月變化如圖11 所示。首先,側(cè)向輸運(yùn)在STMW 生成區(qū)及24.6—25.6 kg/m3密度層上為負(fù)(圖11d), 且整體較其他項(xiàng)小了一個(gè)量級(jí), 可見(jiàn)其對(duì)混合層體積的影響非常微弱。同時(shí), 越密混合項(xiàng)(R)在24.6—25.6 kg/m3密度層上, 對(duì)冬季的混合層體積變化的作用也不顯著(圖11e)。

    在8—11 月, 24.6—25.6 kg/m3之間的混合層水體(黑色虛線之間)的露頭區(qū)主要位于KE 北側(cè), 各月混合層形成體積為正值, 混合層體積持續(xù)增加(圖11a),而這期間海氣形成率為負(fù)值(圖11b), 只有卷吸項(xiàng)能為混合層體積的增大做出貢獻(xiàn)(圖11c)。另外, 圖中紅色的曲線為133°—170°E 經(jīng)向平均的30°N 和36°N 等緯度線。若以此為參照, 可發(fā)現(xiàn)混合層體積年內(nèi)變化強(qiáng)度在不同緯度差異較大。而 8—11 月, 24.6—25.6 kg/m3露頭區(qū)在KE 北部, 海表浮力通量隨密度的向北增加而減小(圖8a 和8b), 因此在這一時(shí)期, 海氣形成體積為負(fù)值。當(dāng)12 月露頭區(qū)回歸STMW 形成區(qū)時(shí), 該區(qū)域具有的弱層結(jié)、高海表浮力通量梯度等特點(diǎn)使得冬季MLD 可以達(dá)到很深的深度, 直至3 月初, 混合層都迅速加深, 混合層形成、海氣形成及卷吸體積都表現(xiàn)出顯著正值。綜上可知, 24.6—25.6 kg/m3露頭區(qū), 在8 月至次年3 月初不斷南移, 混合層形成主要受到混合層底部的卷吸作用調(diào)控; 當(dāng)其露頭區(qū)遷移至STMW 形成區(qū)時(shí)(12 月至3 月), 卷吸作用增強(qiáng)的同時(shí)正的海氣形成也顯著增加, 在兩者共同的作用下混合層迅速加深。

    圖11 133°—170°E 范圍內(nèi)不同密度層上混合層體積收支公式中各項(xiàng)的季節(jié)變化Fig.11 Seasonal variation of each item in the mixed layer volume budget formula at different density layers in 133°—170°E

    為了進(jìn)一步評(píng)估海氣形成作用、卷吸作用等對(duì)生成區(qū)混合層體積變化的影響, 本文計(jì)算了上年10 月至3 月密度處于24.6—25.6 kg/m3之間的各項(xiàng)作用。其中海氣形成體積(~2×1014m3/a)的年代際變化相對(duì)不明顯; 卷吸體積(~3×1014m3/a)與混合層形成體積(~4×1014m3/a)在年變化上幾乎是同步的, 且具有顯著的年代際變化; 側(cè)向輸運(yùn)體積(~0.6×1014m3/a)對(duì)混合層的發(fā)展具有抑制作用, 但相對(duì)其他項(xiàng)小了約一個(gè)量級(jí), 且沒(méi)有明顯的年代際變化。

    3—5 月為潛沉發(fā)生時(shí)期(圖11c), 累積該時(shí)段內(nèi)24.6—25.6 kg/m3密度層上的卷吸率, 得到圖12 藍(lán)色虛線所示的卷吸體積, 其負(fù)值為潛沉體積。由圖12可知潛沉體積與混合層形成體積具有大小相當(dāng)且同步的年代際變化。而潛沉體積年代際變化由混合層形成體積年代際變化決定。綜上, 本次研究進(jìn)一步表明卷吸項(xiàng)主要調(diào)控混合層形成體積的年代際變化, 使生成區(qū)混合層具有顯著的年代際變化, 以此調(diào)控STMW 形成。

    圖12 混合層體積收支中各項(xiàng)及潛沉年變化Fig.12 Annual variation of each item in mixed layer volume budget and spring subduction

    3.3 PDO 調(diào)制

    回顧前文, 影響STMW 年代際變化無(wú)論是MLD還是預(yù)先層結(jié)、EKE 等都與SSHA 緊密聯(lián)系, 因此本文借助ECCO2 的SSHA 數(shù)據(jù)對(duì)其年代際變化進(jìn)行了溯源。

    如圖13c 所示, PDO 指數(shù)呈現(xiàn)出顯著的年代際變化。當(dāng)PDO 處于正(負(fù))相位時(shí)期, 遠(yuǎn)離STMW 形成區(qū)的東太平洋(130°—162°W, 20°—50°N)被激發(fā)出正(負(fù))的風(fēng)應(yīng)力旋度異常(圖13a), 海表為Ekman 輻散(輻聚), 引起海洋上層有向上(下)的Ekman 抽吸速度異常, 導(dǎo)致海表面高度降低(升高), 對(duì)應(yīng)溫躍層抬升(加深); 海表面高度和溫躍層異常信號(hào)再以第一斜壓Rossby 波的速度向西傳輸, 于3—4 a 后抵達(dá)KE 上游南側(cè)STMW 生成區(qū)(圖13b), 可顯著增強(qiáng)(減弱)南側(cè)模態(tài)水生成區(qū)(132°—147°E, 30°—35°N)的預(yù)先層結(jié),從而削弱混合層底部地卷吸過(guò)程, 并最終阻礙(促進(jìn))冬季混合層形成, 使得STMW 形成體積減少(增加)。

    圖13 PDO 指數(shù)及與之有關(guān)的海表風(fēng)應(yīng)力、風(fēng)應(yīng)力旋度及SSHAFig.13 PDO index and the surface wind stress, wind stress curl, and SSHA corresponding to PDO

    4 結(jié)論

    近年來(lái), 許多研究就影響STMW 生成的機(jī)制進(jìn)行了探討, 將20 世紀(jì)80 年代后STMW 的年代際變化歸因于海洋的動(dòng)力過(guò)程, 即由其調(diào)控的海洋預(yù)先層結(jié)引起。然而, 這些動(dòng)力過(guò)程對(duì)STMW 年代際變化的調(diào)節(jié)方式還尚不明確。因此本文采用ECCO2 再分析數(shù)據(jù)對(duì)1992—2019 年STMW 的年代際變化及其機(jī)制進(jìn)行了探討。

    特征分析結(jié)果表明, STMW 形成體積年代際變化完全反映在晚冬生成區(qū)(133°—170°E, 30°—36°N)混合層體積的變化上, STMW 形成厚度和面積均呈現(xiàn)類(lèi)似的年代際變化。其中, STMW 形成厚度的變化由其下邊界深度決定, 而這進(jìn)一步取決于生成區(qū)MLD。另外, 其形成面積的變化與露頭區(qū)的南北收縮、擴(kuò)張有關(guān)。STMW 形成體積在1992—1997 年、2000—2005年和2011—2017 年為正異常, 在1998—1999 年和2006—2010 年為負(fù)異常。正異常期間 MLD 以及STMW 形成厚度顯著增加的區(qū)域?qū)?yīng)正SSHa, 且主要集中于生成區(qū)西部。

    隨后, 本文對(duì)STMW 形成體積正負(fù)異常期間,海洋層結(jié)、海表EKE、海表浮力通量特征進(jìn)行了探討。結(jié)果表明, STMW 形成體積正異常時(shí)期, 在KE 上游南側(cè)(132°—147°E, 30°—35°N)對(duì)應(yīng)有EKE 減小與預(yù)先層結(jié)減弱, 在KE 下游近KE 主軸處對(duì)應(yīng)海表浮力通量升高。而在它們顯著變化的區(qū)域上, 都有正SSHA 與之相對(duì)應(yīng)。KE 上游南側(cè)區(qū)域(132°—147°E,30°—35°N)預(yù)先層結(jié)被削弱(增強(qiáng)), 是促進(jìn)(阻礙)晚冬深MLD 形成主要因素。

    此外, 本文利用了混合層體積收支公式對(duì)混合層形成體積進(jìn)行量化分析。混合層體積變化主要受四個(gè)過(guò)程影響, 它們分別為海氣形成率、沿著等密度面的水體側(cè)向輸運(yùn)、混合層底部卷吸以及越密混合。各密度層上各項(xiàng)引起的體積隨月份的變化顯示,24.6—25.6 kg/m3露頭區(qū), 在8—11 月初不斷南移, 混合層形成主要受到混合層底部的卷吸作用調(diào)控; 12 月開(kāi)始其露頭區(qū)遷移至STMW 形成區(qū)時(shí), 卷吸作用增強(qiáng)的同時(shí)正的海氣形成也顯著增加, 在兩者共同作用下混合層迅速加深。累積上年10 月至3 月24.6—25.6 kg/m3密度層上的各項(xiàng)作用, 得出的各項(xiàng)引起的體積變化量的年時(shí)間序列顯示, 其中海氣形成體積(~2×1014m3/a)以及側(cè)向輸運(yùn)體積(~0.6×1014m3/a)的年代際變化相對(duì)不明顯, 卷吸體積(~3×1014m3/a)與混合層形成體積(~4×1014m3/a)具有大小相當(dāng)且同步的年代際變化。因此, 混合層底部的卷吸作用主要調(diào)控混合層形成體積的年代際變化, 是調(diào)控STMW 形成體積的年代際變化的間接主導(dǎo)因素。

    進(jìn)一步的分析結(jié)果表明, 當(dāng)PDO 處于正(負(fù))相位時(shí)期, 遠(yuǎn)離 STMW 形成區(qū)的中、 東太平洋(130°—162°W, 20°—50°N)被激發(fā)出正(負(fù))的風(fēng)應(yīng)力旋度異常, 海表為Ekman 輻散(輻聚), 引起海洋上層有向上(下)的Ekman 抽吸速度異常, 導(dǎo)致海表面高度降低(升高), 東太平洋產(chǎn)生負(fù)(正)SSHA 信號(hào), 對(duì)應(yīng)溫躍層抬升(加深); 海表面高度和溫躍層異常信號(hào)再以第一斜壓Rossby 波的速度向西傳輸, 于3—4 a 后抵達(dá)KE 上游, 可顯著增強(qiáng)(減弱)KE 南側(cè)模態(tài)水生成區(qū)的預(yù)先層結(jié), 從而削弱混合層底部的卷吸過(guò)程, 并最終阻礙(促進(jìn))冬季混合層形成, 使得STMW 形成體積減少(增加)。

    猜你喜歡
    海表預(yù)先水體
    詞語(yǔ)辨析
    詞語(yǔ)辨析
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    基于無(wú)人機(jī)的海表環(huán)境智能監(jiān)測(cè)系統(tǒng)設(shè)計(jì)與應(yīng)用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對(duì)比分析及初步驗(yàn)證
    太陽(yáng)總輻照度對(duì)熱帶中太平洋海表溫度年代際變化的可能影響
    國(guó)片需自強(qiáng)!2017年12月三大電影排行榜及預(yù)先榜
    廣元:治理黑臭水體 再還水清岸美
    午夜福利视频精品| 涩涩av久久男人的天堂| 免费av中文字幕在线| 日韩大码丰满熟妇| 欧美黑人欧美精品刺激| 日韩视频在线欧美| 精品亚洲成a人片在线观看| 欧美日韩av久久| 最近最新免费中文字幕在线| 久久人人97超碰香蕉20202| 亚洲欧美清纯卡通| 午夜福利视频精品| 一区二区三区四区激情视频| 亚洲七黄色美女视频| 国产免费av片在线观看野外av| 日韩一区二区三区影片| 搡老乐熟女国产| 色综合欧美亚洲国产小说| 久久久国产成人免费| 亚洲av日韩在线播放| 80岁老熟妇乱子伦牲交| 国产精品影院久久| 国产成+人综合+亚洲专区| √禁漫天堂资源中文www| 成人国语在线视频| 久久精品国产亚洲av高清一级| 亚洲欧美日韩另类电影网站| 国产在线视频一区二区| 精品亚洲乱码少妇综合久久| 亚洲自偷自拍图片 自拍| 999久久久精品免费观看国产| 成人黄色视频免费在线看| 久久中文看片网| 中文精品一卡2卡3卡4更新| 国产日韩一区二区三区精品不卡| 99re6热这里在线精品视频| 国产在线一区二区三区精| 十八禁高潮呻吟视频| 亚洲欧美精品综合一区二区三区| 一本综合久久免费| 日韩免费高清中文字幕av| 搡老熟女国产l中国老女人| 美女福利国产在线| netflix在线观看网站| 免费在线观看日本一区| 大片电影免费在线观看免费| 国产区一区二久久| 国产精品成人在线| 纵有疾风起免费观看全集完整版| 十八禁网站网址无遮挡| 国产伦理片在线播放av一区| 久久精品国产a三级三级三级| 国产精品一区二区免费欧美 | 精品欧美一区二区三区在线| 久久久国产一区二区| 夫妻午夜视频| 狠狠狠狠99中文字幕| 久久天堂一区二区三区四区| 中文字幕制服av| 亚洲中文字幕日韩| 在线十欧美十亚洲十日本专区| 一区二区三区四区激情视频| 亚洲欧洲日产国产| 日本av免费视频播放| 亚洲黑人精品在线| 1024香蕉在线观看| 国产成人系列免费观看| 婷婷色av中文字幕| 成年人午夜在线观看视频| 99国产精品99久久久久| 伦理电影免费视频| 多毛熟女@视频| 免费不卡黄色视频| 中文欧美无线码| av网站免费在线观看视频| 国产一区二区三区在线臀色熟女 | 国产日韩欧美亚洲二区| 欧美精品啪啪一区二区三区 | 999久久久精品免费观看国产| 岛国在线观看网站| 色播在线永久视频| 中文字幕色久视频| 中亚洲国语对白在线视频| 黄频高清免费视频| 又大又爽又粗| 超色免费av| 国产精品自产拍在线观看55亚洲 | 自线自在国产av| 9色porny在线观看| 欧美另类亚洲清纯唯美| 人人妻人人澡人人看| 麻豆国产av国片精品| bbb黄色大片| 欧美黄色淫秽网站| 性少妇av在线| 国产一区二区激情短视频 | 亚洲精品日韩在线中文字幕| 亚洲性夜色夜夜综合| 黄频高清免费视频| 国产免费福利视频在线观看| av欧美777| 精品人妻1区二区| 黄色 视频免费看| 又紧又爽又黄一区二区| 国产精品99久久99久久久不卡| 一区二区三区激情视频| 亚洲专区中文字幕在线| 亚洲精品国产区一区二| 国产日韩欧美视频二区| 动漫黄色视频在线观看| 久久久久网色| 免费日韩欧美在线观看| 久热爱精品视频在线9| 超碰97精品在线观看| 老司机在亚洲福利影院| 我的亚洲天堂| 国产免费视频播放在线视频| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三区av网在线观看 | 岛国毛片在线播放| 亚洲中文av在线| videosex国产| 免费高清在线观看日韩| 欧美日韩中文字幕国产精品一区二区三区 | 两个人看的免费小视频| 中文精品一卡2卡3卡4更新| 在线av久久热| 日本撒尿小便嘘嘘汇集6| 日韩有码中文字幕| 久久香蕉激情| 欧美精品一区二区免费开放| 麻豆国产av国片精品| 亚洲美女黄色视频免费看| av网站免费在线观看视频| 汤姆久久久久久久影院中文字幕| 日本一区二区免费在线视频| 日韩大片免费观看网站| 国产极品粉嫩免费观看在线| 一级,二级,三级黄色视频| 日本a在线网址| 超色免费av| 狂野欧美激情性xxxx| 色综合欧美亚洲国产小说| 国产一区二区在线观看av| 香蕉国产在线看| 亚洲国产av影院在线观看| 高清av免费在线| 免费av中文字幕在线| 色精品久久人妻99蜜桃| 一二三四社区在线视频社区8| av视频免费观看在线观看| 国产91精品成人一区二区三区 | 99国产极品粉嫩在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 老司机靠b影院| 亚洲全国av大片| av在线老鸭窝| 久久天堂一区二区三区四区| 亚洲精品美女久久av网站| 捣出白浆h1v1| 亚洲中文av在线| 免费观看av网站的网址| 日韩熟女老妇一区二区性免费视频| 久久国产精品大桥未久av| 精品福利永久在线观看| 麻豆国产av国片精品| 操出白浆在线播放| 黄色视频不卡| 美女高潮喷水抽搐中文字幕| 老司机影院毛片| 1024视频免费在线观看| 最新在线观看一区二区三区| 青草久久国产| 日韩有码中文字幕| 国产成人免费观看mmmm| 久久这里只有精品19| 一区在线观看完整版| 夫妻午夜视频| 日韩欧美国产一区二区入口| 免费在线观看影片大全网站| 久久久久久人人人人人| 男人舔女人的私密视频| 十八禁高潮呻吟视频| 搡老岳熟女国产| 一区二区日韩欧美中文字幕| 宅男免费午夜| 美女高潮喷水抽搐中文字幕| 19禁男女啪啪无遮挡网站| 91成年电影在线观看| 久久人妻福利社区极品人妻图片| 日韩欧美免费精品| 在线看a的网站| 18在线观看网站| 欧美人与性动交α欧美精品济南到| 菩萨蛮人人尽说江南好唐韦庄| av网站免费在线观看视频| 男人操女人黄网站| 免费久久久久久久精品成人欧美视频| 精品高清国产在线一区| 久久久久久久精品精品| 精品少妇久久久久久888优播| 黄片播放在线免费| 曰老女人黄片| 日韩制服丝袜自拍偷拍| 美女中出高潮动态图| 在线观看免费视频网站a站| 欧美在线黄色| 99精品欧美一区二区三区四区| 黄频高清免费视频| 国产一区二区三区av在线| 母亲3免费完整高清在线观看| 捣出白浆h1v1| 国产精品九九99| 亚洲精品自拍成人| 十分钟在线观看高清视频www| 淫妇啪啪啪对白视频 | 免费日韩欧美在线观看| 亚洲欧美一区二区三区久久| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品香港三级国产av潘金莲| 亚洲精品美女久久av网站| 母亲3免费完整高清在线观看| 欧美黄色淫秽网站| 男女国产视频网站| 99国产精品一区二区三区| avwww免费| 亚洲精品中文字幕一二三四区 | 最近最新免费中文字幕在线| 日本vs欧美在线观看视频| 亚洲 欧美一区二区三区| 国产日韩欧美视频二区| 国产成人系列免费观看| 午夜福利视频在线观看免费| 国产成人系列免费观看| 国产在线免费精品| 一区二区三区精品91| 菩萨蛮人人尽说江南好唐韦庄| 少妇人妻久久综合中文| 老司机午夜十八禁免费视频| 亚洲情色 制服丝袜| av免费在线观看网站| 亚洲国产毛片av蜜桃av| 亚洲第一av免费看| 亚洲第一青青草原| 精品久久久久久久毛片微露脸 | 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产av品久久久| 超色免费av| 国产成人一区二区三区免费视频网站| 多毛熟女@视频| 欧美中文综合在线视频| 高潮久久久久久久久久久不卡| 黑人操中国人逼视频| 欧美性长视频在线观看| 婷婷色av中文字幕| 天天影视国产精品| 捣出白浆h1v1| 日韩欧美免费精品| 桃红色精品国产亚洲av| 老司机午夜十八禁免费视频| 伦理电影免费视频| 91九色精品人成在线观看| 97精品久久久久久久久久精品| 男人舔女人的私密视频| 午夜免费观看性视频| 国产三级黄色录像| 亚洲激情五月婷婷啪啪| 日韩制服骚丝袜av| 欧美精品av麻豆av| 中文精品一卡2卡3卡4更新| 亚洲 欧美一区二区三区| 不卡av一区二区三区| 精品久久久精品久久久| 又紧又爽又黄一区二区| 天天影视国产精品| 午夜福利免费观看在线| 少妇裸体淫交视频免费看高清 | 国产深夜福利视频在线观看| 1024香蕉在线观看| 69av精品久久久久久 | 最新的欧美精品一区二区| av又黄又爽大尺度在线免费看| 欧美精品啪啪一区二区三区 | 叶爱在线成人免费视频播放| a在线观看视频网站| 免费人妻精品一区二区三区视频| 精品国产一区二区久久| 欧美另类一区| 无限看片的www在线观看| 精品亚洲乱码少妇综合久久| 国产亚洲欧美在线一区二区| 母亲3免费完整高清在线观看| 一区二区三区精品91| 色综合欧美亚洲国产小说| 丰满迷人的少妇在线观看| 欧美 亚洲 国产 日韩一| 亚洲全国av大片| 久久久水蜜桃国产精品网| 在线观看免费日韩欧美大片| 久久久久国产一级毛片高清牌| 久久ye,这里只有精品| www日本在线高清视频| 纯流量卡能插随身wifi吗| 涩涩av久久男人的天堂| www.999成人在线观看| 日韩大码丰满熟妇| 久久午夜综合久久蜜桃| 久久人人爽人人片av| 精品人妻在线不人妻| 国产欧美亚洲国产| 国产成人影院久久av| 老司机福利观看| 欧美乱码精品一区二区三区| 国产日韩欧美亚洲二区| 黄片小视频在线播放| 国产免费福利视频在线观看| 夜夜夜夜夜久久久久| 热re99久久国产66热| 国产欧美日韩一区二区三区在线| 男女床上黄色一级片免费看| 亚洲激情五月婷婷啪啪| e午夜精品久久久久久久| 不卡av一区二区三区| 在线观看一区二区三区激情| 大香蕉久久网| 如日韩欧美国产精品一区二区三区| 久久影院123| 久久精品国产综合久久久| 国产亚洲av片在线观看秒播厂| 一二三四社区在线视频社区8| 1024视频免费在线观看| 人妻久久中文字幕网| 精品久久久精品久久久| 91字幕亚洲| 国产在线视频一区二区| 69精品国产乱码久久久| 久久青草综合色| 国产高清视频在线播放一区 | 91麻豆av在线| 欧美国产精品va在线观看不卡| 青春草视频在线免费观看| 日韩,欧美,国产一区二区三区| 超色免费av| www.熟女人妻精品国产| 97精品久久久久久久久久精品| 欧美激情 高清一区二区三区| 波多野结衣一区麻豆| 免费一级毛片在线播放高清视频 | 如日韩欧美国产精品一区二区三区| 国产精品一区二区免费欧美 | 免费看十八禁软件| 亚洲av男天堂| 国产高清国产精品国产三级| 国产精品熟女久久久久浪| 国产亚洲av片在线观看秒播厂| 性色av乱码一区二区三区2| av又黄又爽大尺度在线免费看| 手机成人av网站| 亚洲成人国产一区在线观看| 日本a在线网址| 欧美成狂野欧美在线观看| 精品国内亚洲2022精品成人 | 国产欧美日韩一区二区三区在线| 亚洲九九香蕉| 老司机福利观看| 日韩中文字幕欧美一区二区| 岛国毛片在线播放| 日韩制服骚丝袜av| 精品国产国语对白av| 男女午夜视频在线观看| 性色av乱码一区二区三区2| 国产男女内射视频| 成年人午夜在线观看视频| 1024香蕉在线观看| 日韩有码中文字幕| 欧美 亚洲 国产 日韩一| 他把我摸到了高潮在线观看 | 人人妻人人爽人人添夜夜欢视频| 狠狠狠狠99中文字幕| svipshipincom国产片| 亚洲五月色婷婷综合| 国产精品二区激情视频| 少妇人妻久久综合中文| 国产免费一区二区三区四区乱码| 飞空精品影院首页| 男女午夜视频在线观看| 人人妻人人澡人人爽人人夜夜| 高清黄色对白视频在线免费看| 国产成人欧美在线观看 | 亚洲成av片中文字幕在线观看| 成人黄色视频免费在线看| 18在线观看网站| 国产片内射在线| 中文字幕制服av| 狠狠婷婷综合久久久久久88av| 色综合欧美亚洲国产小说| 女性被躁到高潮视频| 久热这里只有精品99| 超色免费av| 免费在线观看日本一区| 成年av动漫网址| 亚洲九九香蕉| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频在线播放观看不卡| 欧美午夜高清在线| 国产一区二区三区av在线| 啦啦啦视频在线资源免费观看| av在线app专区| 黄色视频,在线免费观看| 中文字幕色久视频| 午夜两性在线视频| 成人亚洲精品一区在线观看| 久久久国产一区二区| 成人av一区二区三区在线看 | 美女扒开内裤让男人捅视频| 日韩 亚洲 欧美在线| 欧美午夜高清在线| 久久性视频一级片| 国产成+人综合+亚洲专区| 久久av网站| 国产精品亚洲av一区麻豆| 亚洲中文字幕日韩| 国产免费现黄频在线看| 色综合欧美亚洲国产小说| 丰满饥渴人妻一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品一区三区| 丁香六月天网| 久久性视频一级片| 一本综合久久免费| 1024香蕉在线观看| 波多野结衣一区麻豆| 精品国产乱子伦一区二区三区 | 又大又爽又粗| 久久青草综合色| 天天躁狠狠躁夜夜躁狠狠躁| 国产色视频综合| www日本在线高清视频| 亚洲精品国产精品久久久不卡| av片东京热男人的天堂| 免费在线观看视频国产中文字幕亚洲 | 欧美黄色片欧美黄色片| 日本vs欧美在线观看视频| 国产精品1区2区在线观看. | 欧美成狂野欧美在线观看| 国产精品久久久久成人av| 黄色视频不卡| 精品一区在线观看国产| 999精品在线视频| 欧美精品人与动牲交sv欧美| 精品少妇久久久久久888优播| 亚洲国产毛片av蜜桃av| 一进一出抽搐动态| 久久久久网色| 成人免费观看视频高清| 成人影院久久| 日韩电影二区| 精品国产乱码久久久久久小说| 国产97色在线日韩免费| 国产亚洲午夜精品一区二区久久| 欧美日韩亚洲高清精品| cao死你这个sao货| 亚洲成人手机| 亚洲av美国av| 色综合欧美亚洲国产小说| 伊人久久大香线蕉亚洲五| 老鸭窝网址在线观看| 一级毛片精品| av欧美777| 99热国产这里只有精品6| 久久久久网色| 国产一区二区三区综合在线观看| 十分钟在线观看高清视频www| 久久人人爽人人片av| 大陆偷拍与自拍| 中文字幕最新亚洲高清| 国产在线一区二区三区精| 亚洲激情五月婷婷啪啪| 精品高清国产在线一区| 久热这里只有精品99| av片东京热男人的天堂| 老司机靠b影院| 色精品久久人妻99蜜桃| 无限看片的www在线观看| 国产三级黄色录像| 国产精品一区二区免费欧美 | 久久九九热精品免费| 老司机亚洲免费影院| 老汉色av国产亚洲站长工具| 亚洲国产欧美在线一区| 亚洲国产成人一精品久久久| 这个男人来自地球电影免费观看| 一本大道久久a久久精品| av有码第一页| 国产精品一区二区免费欧美 | 99热网站在线观看| 亚洲精品国产区一区二| 国产精品久久久久久人妻精品电影 | 精品福利永久在线观看| 亚洲欧美一区二区三区久久| 在线天堂中文资源库| 久久久精品94久久精品| 国产精品麻豆人妻色哟哟久久| 国产av又大| 精品福利观看| 天天添夜夜摸| 丰满迷人的少妇在线观看| 日韩欧美国产一区二区入口| 国产av又大| 亚洲精品成人av观看孕妇| 亚洲七黄色美女视频| 欧美日韩国产mv在线观看视频| 国产一区二区三区综合在线观看| 久久久久精品国产欧美久久久 | 黄色视频不卡| 亚洲第一av免费看| 99久久国产精品久久久| 亚洲中文字幕日韩| 桃红色精品国产亚洲av| 狂野欧美激情性bbbbbb| 性色av乱码一区二区三区2| 亚洲国产欧美网| 亚洲精品粉嫩美女一区| 天堂8中文在线网| a级片在线免费高清观看视频| 熟女少妇亚洲综合色aaa.| av在线app专区| 日本一区二区免费在线视频| 另类精品久久| 亚洲精品日韩在线中文字幕| 下体分泌物呈黄色| 欧美xxⅹ黑人| 欧美日韩国产mv在线观看视频| 老汉色av国产亚洲站长工具| 亚洲情色 制服丝袜| 国产不卡av网站在线观看| 国产男女内射视频| 波多野结衣av一区二区av| 性少妇av在线| 丝瓜视频免费看黄片| 久久久国产成人免费| 激情视频va一区二区三区| 日韩人妻精品一区2区三区| 美国免费a级毛片| 久久女婷五月综合色啪小说| 久久国产亚洲av麻豆专区| 黄色视频不卡| 国精品久久久久久国模美| 蜜桃在线观看..| 正在播放国产对白刺激| 丝袜喷水一区| 欧美日韩一级在线毛片| 精品国产超薄肉色丝袜足j| 亚洲精品久久久久久婷婷小说| 美女国产高潮福利片在线看| 成年美女黄网站色视频大全免费| 丁香六月欧美| 亚洲成人国产一区在线观看| www.精华液| 一区二区日韩欧美中文字幕| 大片电影免费在线观看免费| 免费少妇av软件| 欧美一级毛片孕妇| 波多野结衣av一区二区av| 日韩一卡2卡3卡4卡2021年| 亚洲精品国产av蜜桃| 亚洲熟女精品中文字幕| 国产亚洲一区二区精品| 亚洲五月婷婷丁香| 国产日韩欧美视频二区| 日韩人妻精品一区2区三区| 岛国在线观看网站| 国产免费福利视频在线观看| 9191精品国产免费久久| 国产成人精品在线电影| 亚洲精品av麻豆狂野| 国产在线视频一区二区| 国产一区二区 视频在线| 国产免费av片在线观看野外av| 免费久久久久久久精品成人欧美视频| 国产一区有黄有色的免费视频| 十分钟在线观看高清视频www| av天堂在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲伊人色综图| 精品久久久久久电影网| 日韩制服丝袜自拍偷拍| 三级毛片av免费| 人人妻,人人澡人人爽秒播| 久久这里只有精品19| 国产精品影院久久| 夜夜夜夜夜久久久久| 涩涩av久久男人的天堂| 久久久久久亚洲精品国产蜜桃av| 新久久久久国产一级毛片| 国产区一区二久久| 久久精品熟女亚洲av麻豆精品| 国产高清国产精品国产三级| 欧美日本中文国产一区发布| 精品卡一卡二卡四卡免费| 日韩欧美免费精品| 久久精品亚洲熟妇少妇任你| 国产成人精品久久二区二区91| 黄色视频在线播放观看不卡| a在线观看视频网站| 不卡一级毛片| 久久精品久久久久久噜噜老黄| 国产在线一区二区三区精| 两人在一起打扑克的视频| 精品亚洲成国产av| 亚洲精品乱久久久久久| 日本欧美视频一区| 亚洲欧洲日产国产| 亚洲精品久久成人aⅴ小说| 爱豆传媒免费全集在线观看|