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

    南海中南部上層的鹽指與湍流混合過程*

    2014-03-19 12:17:25魯遠(yuǎn)征周生啟岑顯榮郭雙喜尚曉東
    海洋與湖沼 2014年6期
    關(guān)鍵詞:湍流海域南海

    魯遠(yuǎn)征 周生啟 岑顯榮 郭雙喜 尚曉東

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

    跨越等密度面的混合對(duì)海洋中物理特征的形成有極其重要的作用, 是控制海洋中物理特征分布的有效動(dòng)力因素, 它影響著營(yíng)養(yǎng)物、溶解氣體、顆粒物的垂向輸移, 驅(qū)動(dòng)著邊緣海與大洋之間的水體交換,甚至影響著大洋熱鹽環(huán)流, 進(jìn)而影響全球氣候(MacKinnonet al, 2003; Rahmstorf, 2003; Tianet al,2009)。海洋內(nèi)部混合造成的浮力(熱量)平衡及渦動(dòng)力作用對(duì)海洋上升流強(qiáng)度及水平環(huán)流至關(guān)重要(Polzinet al, 1997)。Munk等(1998)估算了維持全球徑向翻轉(zhuǎn)環(huán)流所需的平均跨越等密度面的混合率應(yīng)達(dá)到10–4m2/s的量級(jí), 而大量海洋觀測(cè)發(fā)現(xiàn)海洋內(nèi)部混合較弱, 約為10–5m2/s的量級(jí)。因此, 邊緣海較強(qiáng)的混合不僅影響到了局地, 甚至對(duì)全球大洋環(huán)流與熱量的輸運(yùn)都起到了重要的作用。近年來, 邊緣?;旌涎芯康玫搅藦V泛的關(guān)注。中國(guó)南海是熱帶西太平洋最大的邊緣海之一, 也是廣泛被認(rèn)為混合較強(qiáng)的海域(Laurent, 2008; Tianet al, 2009; Liuet al, 2012; Yanget al, 2014 )。Laurent(2008)在東沙群島附近的坡折區(qū)發(fā)現(xiàn)耗散率可高達(dá)10–6W/kg。Yang等(2014)在南海北部陸架區(qū)發(fā)現(xiàn)平均耗散率可達(dá)2.8×10–7W/kg, 而平均混合率為1.6×10–3m2/s。Tian等(2009)根據(jù)南海熱量守恒估算南海的混合率約為1.8×10–3m2/s, 而在呂宋海峽底海實(shí)際觀測(cè)甚至高達(dá)10–2m2/s的量級(jí)。

    南海強(qiáng)烈的混合被認(rèn)為主要是由于內(nèi)波的破碎所造成。地形的劇烈變化與正壓潮的相互作用產(chǎn)生了大量的內(nèi)波, 而內(nèi)波的破碎為混合提供了充足的能量。目前有關(guān)南?;旌系难芯恐饕杏谀虾5谋辈筷懠軈^(qū)及呂宋海域, 而這些海域也被認(rèn)為是具有高內(nèi)潮能量的海域(Laurent, 2008; Liuet al, 2012)。另外,北太平洋熱帶水(NPTW), 北太平洋中層水(NPIW)及北太平洋深層水(NPDW)通過呂宋海峽進(jìn)入南海, 形成南海次表層高溫高鹽水團(tuán)和南海中層低溫低鹽水團(tuán)的溫鹽結(jié)構(gòu)(Tianet al, 2009)。這種穩(wěn)定的垂向溫度梯度及不穩(wěn)定的垂向鹽度梯度使得該深度范圍內(nèi)存在鹽指生成的環(huán)境。鹽指也同樣可以產(chǎn)生有效的跨越等密度面的混合(Laurentet al, 1999; Schmittet al,2005), 特別是會(huì)影響水團(tuán)與水團(tuán)之間的相互作用,改變水體的溫度與鹽度結(jié)構(gòu), 進(jìn)而影響大尺度環(huán)流甚至氣候的變化(Johnsonet al, 2009)。Lambert等(1977)指出由鹽指造成的有效混合率可達(dá)到5×10–4m2/s, 而Schmitt等(2005)利用示蹤劑及微結(jié)構(gòu)剖面儀同時(shí)觀測(cè)的結(jié)果表明鹽指所引起的混合均可達(dá)到10–4m2/s的量級(jí)。雖然有學(xué)者在南海也發(fā)現(xiàn)了存在鹽指的溫鹽臺(tái)階結(jié)構(gòu)(楊麗芬等, 2008), 但鹽指對(duì)南?;旌系呢暙I(xiàn)有多大尚不清楚。

    為了解夏季南海中南部上層混合特征, 評(píng)估鹽指和湍流對(duì)南?;旌系呢暙I(xiàn)及探究產(chǎn)生該混合特征的可能機(jī)制, 本文利用2012年8月12日至9月5日南海海洋湍流微結(jié)構(gòu)剖面儀(TurboMap)和溫鹽深剖面儀(CTD)觀測(cè)資料, 對(duì)南海中南部海域上層500m以淺的混合進(jìn)行了分析。

    1 現(xiàn)場(chǎng)觀測(cè)

    2012年8月12日至9月5日, 中國(guó)科學(xué)院南海海洋研究所“實(shí)驗(yàn)3號(hào)”科考船開展了南海海洋斷面科學(xué)考察。觀測(cè)范圍(6°—18.25°N; 108.5°—119.33°E)內(nèi)共有75個(gè)常規(guī)測(cè)站, 均有CTD觀測(cè), 其中在53個(gè)站位進(jìn)行了TurboMap微結(jié)構(gòu)觀測(cè), 觀測(cè)區(qū)域和站位見圖1。CTD為美國(guó)Sea-Bird公司生產(chǎn), 型號(hào)為SBE911plus, 采樣頻率為24Hz, 最大觀測(cè)水深為1500m。自由沉降式微結(jié)構(gòu)剖面儀TurboMap-L由日本Alec Electronics公司生產(chǎn), 可獲得垂向高分辨率的流速剪切(?u/?z)、溫度梯度(?T/?z)、快速溫度(FP07)、鹽度、壓力及葉綠素等信息。TurboMap-L的最大采樣頻率為512Hz, 最大測(cè)量水深為500m,觀測(cè)時(shí)自由下降速度為0.5—0.7m/s。剪切探頭噪聲主要由儀器的振動(dòng)產(chǎn)生, 其量級(jí)約10–10W/kg(Wolket al, 2002)。

    將所有測(cè)站CTD數(shù)據(jù)綜合可得到該海域的溫鹽圖(圖2), 其顏色表示不同的緯度, 由于南海次表層高溫高鹽水團(tuán)和南海中層低溫低鹽水團(tuán)存在, 南海溫鹽結(jié)構(gòu)呈明顯的倒“S”形狀。Tian等(2009)指出北太平洋熱帶水團(tuán)(NPTW)和北太平洋中層水團(tuán)(NPIW)溫鹽圖呈倒“S”形狀, 其倒“S”形狀從太平洋到呂宋海峽再到南海逐漸變?nèi)?。北太平洋熱帶水團(tuán)鹽度減小了0.3, 而同時(shí)北太平洋中層水團(tuán)鹽度增加了0.3, 表明跨越等密度面的混合可能起到了重要的作用。而由圖2可見, 在最靠近呂宋海峽的測(cè)站, 溫鹽圖中次表層出現(xiàn)了鹽度為34.8的高鹽水, 反映了太平洋高鹽水的入侵。另外密度在25—27之間的水體明顯體現(xiàn)出由北向南倒“S”形狀也在逐漸變?nèi)? 也表明混合對(duì)南海內(nèi)區(qū)水體溫鹽結(jié)構(gòu)也起到了重要的作用。前人研究推測(cè)造成該混合因素可能有鹽指和湍流作用(楊麗芬等,2008; Laurent, 2008; Tianet al, 2009; Liuet al, 2012)。然而南海中南部混合的研究很少, 其混合特征尚不清楚。

    圖1 2012年南海海洋斷面科學(xué)考察站位分布(2012/8/12—9/5), 共75個(gè)站位, 各站均有CTD觀測(cè), 其中實(shí)心點(diǎn)(53個(gè))為同時(shí)有TurboMap微結(jié)構(gòu)觀測(cè)站點(diǎn), 圖中等值線表示水深(m)Fig.1 Sampling site of the cruise (2012/8/12—9/5) in the southern South China Sea. CTD measurements were taken at all stations, and TurboMap measurements are marked as solid diamonds, black contours indicate the depth(m)

    圖2 南海1500m以淺海水溫鹽圖, 顏色表示緯度(°N),小圖為徑度為113°E斷面站位且位密度在25—27之間的溫鹽圖Fig.2 Diagram of potential temperature vs salinity above 1500m in the South China Sea; the color bar indicates the latitude(°N) of stations, the θ-S diagram along the section 113°E is enlarged in the inset

    2 鹽指混合

    雙擴(kuò)散作用普遍存在于海洋中, 全球約44%的海洋發(fā)生雙擴(kuò)散現(xiàn)象, 其中全球30%海洋發(fā)生的雙擴(kuò)散類型為鹽指(You, 2002)。特別是在亞熱帶海域,水體的特定垂向溫鹽結(jié)構(gòu)更易形成鹽指, 這主要是由于上層高溫高鹽的中央水團(tuán)(central water)與其下層來自高緯的冷而淡的水團(tuán)相互作用的結(jié)果, 例如亞熱帶大西洋及印度洋密度躍層區(qū)(Schmitt, 1981,2003; You, 2002; Inoueet al, 2008)。

    2.1 鹽指環(huán)境診斷

    You(2002)將鹽指形成的原因歸納為三種: 1, 海洋的蒸發(fā)量大于降雨量(中央水團(tuán)); 2, 海洋環(huán)流造成的水團(tuán)的影響, 例如: 北大西洋深層水(NADW)與繞極深層水(CDW); 3, 邊緣海水注入的影響(地中海及紅海等)。其中因素1為全球鹽指形成的主要的原因。南海的緯度范圍是4°—12°N, 為低緯熱帶地區(qū)非中央水團(tuán)形成區(qū)域。中國(guó)南海海域30年平均蒸發(fā)減降雨量如圖3所示。降水資料采用美國(guó)國(guó)家海洋大氣局(NOAA)氣候預(yù)報(bào)中心(CPC)根據(jù)雨量器資料、衛(wèi)星資料和NCEP/NCAR 再分析資料等, 通過融合算法生成的全球逐月平均(CMAP)降水, 蒸發(fā)量是采用美國(guó)伍茲霍爾海洋研究所OAflux研究項(xiàng)目分析的海氣通量數(shù)據(jù), 時(shí)間從1979年1月到2008年7月, 空間分辨率2.5°×2.5°。由圖3可見降雨量普遍大于蒸發(fā)量,表明在南海海域自身無(wú)法形成高溫高鹽的次表層水,即不滿足鹽指形成的因素1。另外, 南海沒有大量的陸源徑流及邊緣水體注入的影響, 即不滿足因素3。但由于亞熱帶的太平洋水通過呂宋海峽進(jìn)入南海,衍生出高溫高鹽的次表層水。這樣的垂向溫鹽結(jié)構(gòu)可能為鹽指生長(zhǎng)提供必要的環(huán)境, 即可能有因素2的影響。國(guó)內(nèi)已有學(xué)者在南海北部觀測(cè)到有類似溫鹽臺(tái)階的結(jié)構(gòu)(楊麗芬等, 2008), 因此要首先診斷該研究海域是否有鹽指成長(zhǎng)的環(huán)境。

    圖3 1979年1月到2008年7月南海平均蒸發(fā)減降雨量分布Fig.3 Contours of the annual-mean field of evaporation(E)minus precipitation(P) from January 1979 to July 2008

    Turner(1973)指出雙擴(kuò)散的強(qiáng)弱可以用溫度與鹽度的梯度比值來衡量, 該比值稱為穩(wěn)定度比或密度比Rρ。密度比越接近1, 鹽指越有可能生成。發(fā)生鹽指不穩(wěn)定的條件為:

    Ruddick (1983)首次引用Turner角的診斷方法來判斷鹽指生長(zhǎng)的強(qiáng)弱, 其表達(dá)式如方程(2)所示:

    利用Turner角來診斷鹽指的信號(hào)比用密度比更清晰且更容易界定和解釋, 不易產(chǎn)生混淆。Turner角45°—90°為鹽指發(fā)生區(qū)間, –45°—–90°為擴(kuò)散對(duì)流發(fā)生區(qū)間(雙擴(kuò)散的別一種類型), –45°—45°為水體雙穩(wěn)定區(qū)間, 其他為水體不穩(wěn)定區(qū)間。Turner角越接近90°越易于生成鹽指, 越接近45°越不易生成鹽指。You(2002)對(duì)Turner角的物理意義做了更仔細(xì)的闡述,將整個(gè)鹽指發(fā)生區(qū)間(45°—90°)又分別用75°與60°進(jìn)一步分為強(qiáng)、中、弱三個(gè)區(qū)間。Turner角與密度比的關(guān)系是Rρ=-tan(Tu+45), 即45°與90°分別對(duì)應(yīng)的密度比為±∞和1。

    圖4(a)為緯向18°N, 10°N, 6°N和徑向113°E斷面的Turner角分布, 可見在200—500m的密度躍層區(qū)均存在鹽指形成的環(huán)境, 但其Turner角度都偏小。其垂向溫、鹽及Turner角的垂向廓線如圖4(b)(117.4937°E, 17.9929°N), 可見在由次表層高溫高鹽水團(tuán)到中層低溫低鹽水團(tuán)過渡的區(qū)間, Turner角大于45°, 但都小于60°, 屬于弱鹽指發(fā)生區(qū)間(其他站點(diǎn)垂向廓線結(jié)構(gòu)相似)。

    圖4 Turner角的觀測(cè)斷面分布及垂向廓線分布Fig.4 The distributions of the Turner angle along the sections and the profiles

    采用Schmitt(1981)的方法, 對(duì)發(fā)生鹽指的相關(guān)水團(tuán)的溫鹽圖進(jìn)行線性擬合, 可得到大尺度水團(tuán)間的平均密度比。計(jì)算得到的南海次表層水團(tuán)與南海中層水團(tuán)之間的密度比為11.27, 即Turner角約為50°,遠(yuǎn)大于全球大洋中央水團(tuán)的密度比平均值2.55。Schmitt(1981; 2003)對(duì)不同大洋的中央水團(tuán)進(jìn)行了分析得出, 在密度比大于1.7時(shí)或?qū)咏Y(jié)指數(shù)(Layer Index,即溫度梯度小于0.01°C/m的概率)較小時(shí)(小于0.3),其溫鹽臺(tái)階結(jié)構(gòu)便不明顯。將所有觀測(cè)水深在200至500m的數(shù)據(jù)在4m的間隔計(jì)算層結(jié)指數(shù)、密度比及Turner角。這三個(gè)表征鹽指生成的強(qiáng)弱程度的指標(biāo)統(tǒng)計(jì)如圖5所示。層結(jié)指數(shù)為0.161, 明顯小于0.3; 密度比主要集中在30左右, 遠(yuǎn)大于1.7; 而Turner角基本都分布在60°以下, 主要集中在48°。以上參數(shù)均表明, 鹽指發(fā)育非常弱。而在該研究區(qū)域溫、鹽廓線證實(shí)有類似于臺(tái)階狀的結(jié)構(gòu), 但其臺(tái)階特征很不明顯(圖6b)。

    由以上分析可以看出, 在南海的中南部海域, 由于水團(tuán)之間特定的溫鹽結(jié)構(gòu)使得該海域存在鹽指生成的環(huán)境, 但鹽指診斷參量結(jié)果表明鹽指發(fā)育非常弱, 和實(shí)際觀測(cè)的溫鹽臺(tái)階結(jié)構(gòu)不明顯相一致。

    2.2 鹽指混合作用

    既然南海存在鹽指的形成環(huán)境, 那么鹽指所起到的混合作用有多大呢?St. Laurent等(1999)指出即使在湍流混合發(fā)生的環(huán)境下, 鹽指信號(hào)也可能被識(shí)別出來。如果能將鹽指的信號(hào)清晰地識(shí)別出來, 就可以定量地評(píng)估鹽指對(duì)混合的貢獻(xiàn)。Schmitt(2003)和Inoue等(2007)指出鹽指具有高的溫度耗散率(χ)、高的混合效率(Γ)、低湍動(dòng)能耗散率(ε)及低浮性雷諾數(shù)(Rε)等特征, 即“高χ,Γ; 低ε,Rε”, 而湍流則具有“高χ,ε,Rε; 低?!碧卣?。

    圖5 表征鹽指信號(hào)發(fā)生強(qiáng)弱的參數(shù)統(tǒng)計(jì)Fig.5 The distribution of the parameters characterizing the strength of the salt finger signal

    圖6 鹽指發(fā)生區(qū)的剪切、溫度、溫度梯度、湍動(dòng)能耗散率及溫度耗散率廓線(113.5015°E, 9.9953°N)Fig.6 Typical profiles of shear, temperature, temperature gradient, turbulent kinetic energy dissipation rate and thermal dissipation rate in the finger-favorable region(113.5015°E, 9.9953°N)

    由于微結(jié)構(gòu)剖面儀(TurboMap)可以獲得高分辨率的流速剪切, 因此可直接通過速度剪切譜φ(k)積分的方法計(jì)算湍動(dòng)能耗散率(Oakey, 1982; Wolket al,2002)。

    其中,ν為運(yùn)動(dòng)黏性系數(shù),k為垂向波數(shù), 積分下限k1取1cpm, 積分上限k2為噪聲影響的最大波數(shù)。

    利用TurboMap的快速溫度探頭(FP07)可獲得高分辨率的溫度垂向剪切, 從而可積分溫度垂向剪切譜ψ(k)計(jì)算溫度耗散率(Oakey, 1982)。

    其中,κθ為分子熱擴(kuò)散系數(shù)。

    圖6 為微結(jié)構(gòu)剖面儀觀測(cè)所得到的南海次表層水團(tuán)與中層水團(tuán)之間流速剪切、溫度及溫度梯度廓線,根據(jù)式(3)和(4)計(jì)算的動(dòng)能耗散率及溫度耗散率, 其中積分上限為儀器震蕩噪聲干擾的最大波數(shù), 下限波數(shù)取1cpm, 平均間隔為4m。局部放大(圖6b)可見位溫有類似臺(tái)階狀結(jié)構(gòu), 但并不明顯。且從圖6中可以看出, 湍動(dòng)能耗散率及溫度耗散率垂向變化較為均勻, 沒有呈現(xiàn)出鹽指混合為主導(dǎo)的水體環(huán)境的特征——溫度耗散率在臺(tái)階界面處明顯的極大值而在臺(tái)階混合區(qū)極小值交替出現(xiàn)的特征(Schmittet al, 2005)。

    對(duì)于湍流的情況, 熱、鹽與密度的有效擴(kuò)散系數(shù)認(rèn)為是相同的, 即κT=κS=κρ。采用Osborn(1980)的模型來計(jì)算有效密度擴(kuò)散系數(shù)為:

    其中,N2為浮力頻率的平方, 利用CTD獲得。Γ為混合效率, 在湍流中Γ=Jb/ε(Jb為浮力通量), 即混合的能量一部分轉(zhuǎn)化為勢(shì)能, 另外的基本被耗散掉了,混合效率表征二者的能量之比。湍混合模型中一般取Rf為通量理查森數(shù)(Osborn, 1980;Schmitt, 2003; Inoueet al, 2007)。而實(shí)際水體的混合效率在得到ε后可以由觀測(cè)獲得, 即:

    在有鹽指發(fā)生的區(qū)域, 觀測(cè)得到的耗散率為湍混合與鹽指混合共同作用的結(jié)果。但由于鹽指與湍流所造成的密度通量是相反的, 即鹽指所造成的密度通量垂直向下, 而湍流所造成的密度通量垂直向上。所以觀測(cè)得到的混合效率并不是湍流與鹽指混合效率簡(jiǎn)單相加, 而表示為整個(gè)混合過程中所產(chǎn)生的凈浮力通量和耗散量之比。對(duì)于鹽指, 熱、鹽與密度的有效擴(kuò)散系數(shù)是各不相同的(Schmitt, 2003)。

    其中γ為浮力通量比, 即熱通量與鹽通量的比值:T',S',w'分別為溫度、鹽度與垂向速度的脈動(dòng)。而參數(shù)化鹽指的密度有效擴(kuò)散系數(shù)時(shí), 混合效率應(yīng)為(St. Laurentet al, 1999):

    但需要提到的是, 該混合效率是在密度比在1到10范圍內(nèi)得到的。

    由式(3)、(4)和(7)可以求得湍動(dòng)能耗散率、溫度耗散率及有效熱擴(kuò)散系數(shù)。將其參量在200—500m水深進(jìn)行平均后得到其按緯度的分布(圖7)。18°N斷面的動(dòng)能耗散率明顯小于其他地區(qū), 約小一個(gè)量級(jí),而溫度耗散率及熱有效擴(kuò)散系數(shù)不同緯度間相差不大。即18oN斷面呈現(xiàn)出相對(duì)的“高χ; 低ε”特征, 而南海南部呈現(xiàn)出相對(duì)“高χ;ε”特征。

    圖7 南海不同緯度次表層水團(tuán)與中層水團(tuán)之間的平均動(dòng)能, 溫度耗散率及有效熱擴(kuò)散系數(shù)Fig.7 The average kinetic energy, thermal dissipation rate and effective thermal diffusivity between subsurface water and intermediate water of South China Sea

    為進(jìn)一步區(qū)別雙擴(kuò)散與湍流, 通常采用一個(gè)最為常用的參數(shù), 即浮性雷諾數(shù):

    Ruddick等(1997)發(fā)現(xiàn)在10120, 基本可以認(rèn)為湍流占主導(dǎo)地位。當(dāng)雙擴(kuò)散占主導(dǎo)時(shí), 浮力通量基本可以與耗散達(dá)到平衡(Jb≈ε), 即浮力雷諾數(shù)就等于了Stern數(shù)Jb/(νN2), 在雙擴(kuò)散溫鹽臺(tái)階界面處Rε<25(Inoueet al, 2007)。

    由式(6)和(10)可分別求得南海次表層水團(tuán)到南海中層水團(tuán)間的混合效率和浮性雷諾數(shù)。由圖7可知18°N斷面與南海南部海域混合有明顯的不同, 因此分別將全部數(shù)據(jù)、18°N斷面和南海南部進(jìn)行統(tǒng)計(jì), 計(jì)算的這兩個(gè)參數(shù)見圖8。南海中部及南部整體混合效率小于湍流特征值0.2, 南海南部海域更明顯小于0.2,而在18°N斷面接近0.2, 略有偏小。在湍流環(huán)境中,混合效率0.2為Osborn(1980)模型中混合效率的上限,因此可知本文所研究的海域整體呈現(xiàn)為湍流的環(huán)境。南海中部及南部整體浮性雷諾數(shù)大于臨界值20, 南海南部海域更明顯大于20, 但18°N斷面主要集中在小于20的區(qū)間內(nèi)。因此可知18°N斷面存在鹽指的信號(hào), 呈現(xiàn)出相對(duì)的“高Γ; 低Rε”特征, 而南海南部呈現(xiàn)出相對(duì)“低Γ; 高Rε”特征。其研究海域整體呈現(xiàn)出“低Γ; 高Rε”的湍流占主導(dǎo)的特征。

    綜上, 由計(jì)算得到的湍動(dòng)能耗散率、溫度耗散率、熱有效擴(kuò)散系數(shù)、混合效率及浮性雷諾數(shù)可知,18°N斷面呈現(xiàn)出相對(duì)的“高χ,Γ; 低ε,Rε”特征, 而南海南部呈現(xiàn)出相對(duì)“高χ,ε,Rε; 低?!钡奶卣? 即表明18°N度斷面明顯存在鹽指信號(hào), 其鹽指的作用明顯強(qiáng)于南海南部。但整體海域仍然呈現(xiàn)出“低Γ; 高Rε”的湍流占主導(dǎo)的特征。

    3 跨越等密面混合特征

    南海被普遍認(rèn)為是混合比較強(qiáng)的邊緣海, 很多學(xué)者已經(jīng)發(fā)現(xiàn)了南海東北部呂宋海峽及北部陸架海域存在較強(qiáng)的混合耗散區(qū)(Laurent, 2008; Tianet al.,2009; Liuet al, 2012; Yanget al, 2014), 但南海中部及南部的混合特征鮮有研究。通過上文的分析可知, 南海18°N斷面明顯與南部海域的混合特征有所不同,需要進(jìn)一步探究南海中部及南部上層混合的特征及可能造成混合的原因。

    圖9為整個(gè)航次觀測(cè)的緯向18°N, 10°N, 6°N和徑向113°E四個(gè)主要斷面的湍動(dòng)能耗散率及鹽度廓線, 從圖中可知18°N斷面的耗散率由表層至中層500m左右均較弱, 相對(duì)10°N與6°N斷面小一個(gè)量級(jí)。圖10為整個(gè)航次觀測(cè)的四個(gè)主要斷面的溫度耗散率及位溫廓線, 從圖中可知18°N斷面的溫度耗散率雖然與10°N和6°N斷面的值相差不大, 但仍有偏小。由此可知, 18°N斷面的混合明顯弱于南部海域。

    圖8 南海次表層水團(tuán)與中層水團(tuán)之間的混合效率及浮性雷諾數(shù)統(tǒng)計(jì), 從上至下分別是南海所有測(cè)站、18°N斷面測(cè)站及18°N以南測(cè)站Fig.8 The distribution of the mixing efficiency and the buoyancy Reynolds number Rε between subsurface water and intermediate water of South China Sea. The upper panel represents those of the all stations, middle panel for those of the station at18°N section, and the lower panel for those of the station south of the 18°N section.

    圖9 南海緯向18°N, 10°N, 6°N和徑向113oE斷面觀測(cè)的動(dòng)能耗散率廓線(彩色)及鹽度廓線(黑線)Fig.9 Turbulent kinetic energy dissipation across section of 18°N, 10°N, 6°N and 113°E in the South China Sea, respectively (color lines). The salinity profiles at each station are plotted in black lines

    為進(jìn)一步探究其混合特征, 分別選取18°N及10°N斷面中部的兩個(gè)深水站位進(jìn)行比較。將這兩個(gè)站位有微結(jié)構(gòu)觀測(cè)深度以內(nèi)的溫度、鹽度、Turner角、浮性頻率、動(dòng)能耗散率及溫度耗散率廓線繪于圖11(115.9988°E, 18.0063°N)和圖12(114.5010°E,10.1165°N)。由圖7a、圖9和圖11可知18°N斷面強(qiáng)的混合主要位于上部混合層, 而混合層以下隨著深度的增加而減小, 但平均約為10–9W/kg量級(jí)。由圖7a、圖9和圖12可知南部混合層依然是強(qiáng)的混合區(qū),混合層以下混合雖然略有下降但仍然較強(qiáng), 明顯不同于18°N斷面。南海南部500m以淺全水柱均呈現(xiàn)出較強(qiáng)的混合, 平均約為10–8W/kg量級(jí)。

    圖10 南海緯向18°N, 10°N, 6°N和徑向113°E斷面觀測(cè)的溫度耗散率廓線(彩色)及位溫廓線(黑線)Fig.10 Temperature dissipation rate across the section of 18°N, 10°N, 6°N and 113°E in the South China Sea, respectively (color lines).The potential temperature profiles at each station are plotted in black lines

    圖11 南海18°N上層500m以淺的溫度、鹽度、Turner角、浮性頻率、動(dòng)能耗散率及溫度耗散率廓線(115.9988°E, 18.0063°N)Fig.11 Characteristic profiles of temperature, salinity, Turner angle, buoyancy frequency, kinetic energy dissipation rate and thermal dissipation rate of the lower 500m at 18°N section(115.9988°E, 18.0063°N)

    由浮性頻率N可知從南海次表層水團(tuán)至南海中層水團(tuán), 18°N處的層結(jié)強(qiáng)于10°N, 因此18oN的鹽指生成強(qiáng)度應(yīng)該大于10°N(圖11, 12)。由Turner角可確定有鹽指生成的環(huán)境的深度范圍, 而在該深度范圍內(nèi)10°N的混合卻明顯強(qiáng)于18°N, 鹽指所造成的混合并不能影響背景湍流混合的特征。那么南海南部的強(qiáng)湍流混合是什么原因造成的呢?

    近年來, 很多學(xué)者發(fā)現(xiàn)在粗糙地形處, 局地跨越等密面混合率往往大于10–5m2/s的量級(jí), 甚至是海洋內(nèi)部混合的100—1000倍, 例如海山, 洋脊, 海盆間通道等(Polzinet al, 1997; Tianet al, 2009)。而從圖9和圖10的地形上看, 18°N斷面主要位于海盆區(qū), 而南部的站位地形較淺且較為粗糙復(fù)雜。南海南部較強(qiáng)的混合很可能是由于內(nèi)潮波與局地地形相互作用造成內(nèi)波破碎而產(chǎn)生的。由Jan等(2007)和Shang等(2014)通過溫鹽數(shù)據(jù)、正壓潮流和地形數(shù)據(jù)計(jì)算的整個(gè)南海的內(nèi)潮強(qiáng)迫項(xiàng)分布可知, 在呂宋海峽、西沙群島及南沙群島海域均是內(nèi)潮的多發(fā)海域。由于這些地形復(fù)雜且變化較快, 也是內(nèi)潮能量耗散強(qiáng)的區(qū)域(Laurent, 2008; Tianet al, 2009; Liuet al, 2012)。

    圖12 南海10°N上層500m以淺的溫度、鹽度、Turner角、浮性頻率、動(dòng)能耗散率及溫度耗散率廓線(114.5010°E, 10.1165°N)Fig.12 Characteristic profiles of temperature, salinity, Turner angle, buoyancy frequency, kinetic energy dissipation rate and thermal dissipation rate of the lower 500m at 10°N section(114.5010°E, 10.1165°N)

    綜上, 南海中部及南海上層混合的特征表現(xiàn)為,18°N斷面的中部海域整個(gè)上層海洋平均耗散率約10–10—10–9W/kg的量級(jí), 而南海南部海域混合率可達(dá)到10–8W/kg的量級(jí), 18°N整體弱于南海南部海域。南海南部的上層海洋500m以淺全水柱均呈現(xiàn)出較強(qiáng)的混合耗散, 隨深度的增加基本沒有減弱的特征。結(jié)合地形及南海內(nèi)潮強(qiáng)迫項(xiàng)分布推測(cè), 南海南部較強(qiáng)的混合可能是由于內(nèi)潮波在南部較淺且復(fù)雜的地形作用下的破碎造成的。

    4 結(jié)論

    南海次表層高溫高鹽水團(tuán)和南海中層低溫低鹽水團(tuán)的特殊溫鹽結(jié)構(gòu), 導(dǎo)致南海存在生成鹽指的環(huán)境。通過計(jì)算層結(jié)指數(shù)、密度比及Turner角診斷指數(shù), 發(fā)現(xiàn)鹽指發(fā)育較弱, 從溫鹽廓線上也沒有發(fā)現(xiàn)明顯的溫鹽臺(tái)階結(jié)構(gòu)。相對(duì)于南海南部海域, 18°N斷面水體呈現(xiàn)出“高χ,Γ; 低ε,Rε”特征, 而南海南部呈現(xiàn)出相對(duì)“高χ,ε,Rε; 低?!钡奶卣? 表明18°N斷面處鹽指的作用強(qiáng)于南海南部。從湍流及鹽指對(duì)混合的貢獻(xiàn)來看, 鹽指所造成的混合作用非常有限, 南海中部及南部的上層混合還是以湍流混合為主導(dǎo)。從上層水體的混合特征來看, 呈現(xiàn)出北部弱(18°N斷面)、南部強(qiáng)的特征。其垂向結(jié)構(gòu)上看, 18°N斷面表層混合較強(qiáng)、次表層及中層較弱, 而南部的斷面呈現(xiàn)出500m以淺全水柱混合均較強(qiáng)的特征。結(jié)合南海地形及南海內(nèi)潮強(qiáng)迫項(xiàng)分布, 可推測(cè)南部較強(qiáng)的混合可能是由于內(nèi)潮與地形相互作用造成的。

    楊麗芬, 田紀(jì)偉, 謝玲玲, 2008. 南海北部陸坡區(qū)混合過程觀測(cè). 海洋科學(xué), 32(12): 10—16

    Inoue R, Kunze E, St Laurent Let al, 2008. Evaluating salt-fingering theories. Journal of Marine Research, 66(4):413—440

    Inoue R, Yamazaki H, Wolk Fet al, 2007. An estimation of buoyancy flux for a mixture of turbulence and double diffusion. Journal of Physical Oceanography, 37(3): 611—624

    Jan S, Chern C S, Wang Jet al, 2007. Generation of diurnal K1 internal tide in the Luzon Strait and its influence on surface tide in the South China Sea. Journal of Geophysical Research: Oceans, 112(C6): C06019

    Johnson G C, Kearney K A, 2009. Ocean climate change fingerprints attenuated by salt fingering? Geophysical Research Letters, 36(21), L21603

    Lambert Jr R B, Sturges W, 1977. A thermohaline staircase and vertical mixing in the thermocline. Deep Sea Research,24(3): 211—222

    Laurent L S, Schmitt R W, 1999. The contribution of salt fingers to vertical mixing in the north Atlantic tracer release experiment. Journal of Physical Oceanography, 29(7): 1404—1424

    Laurent L S, 2008. Turbulent dissipation on the margins of the South China Sea. Geophysical Research Letters, 35(23),L23615

    Liu Z, Lozovatsky I, 2012. Upper pycnocline turbulence in the northern South China Sea. Chinese Science Bulletin, 57(18):2302—2306

    MacKinnon J, Gregg M, 2003. Mixing on the late-summer New England Shelf-solibores, shear, and stratification. Journal of Physical Oceanography, 33(7): 1476—1492

    Munk W, Wunsch C, 1998. Abyssal recipes II: energetics of tidal and wind mixing. Deep-Sea Research Part I, 45(12): 1977—2010

    Oakey N, 1982. Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shear microstructure measurements. Journal of Physical Oceanography, 12(3): 256—271

    Osborn T, 1980. Estimates of the local rate of vertical diffusion from dissipation measurements. Journal of Physical Oceanography, 10(1): 83—89

    Polzin K, Toole J, Ledwell Jet al, 1997. Spatial variability of turbulent mixing in the abyssal ocean. Science, 276(5309):93—96

    Rahmstorf S, 2003. Thermohaline circulation: The current climate. Nature, 421(6924): 699—699

    Ruddick B, 1983. A practical indicator of the stability of the water column to double-diffusive activity. Deep Sea Research Part I: Oceanographic Research Papers, 30(10):1105—1107

    Ruddick B, Walsh D, Oakey N, 1997. Variations in apparent mixing efficiency in the North Atlantic Central Water.Journal of Physical Oceanography, 27(12): 2589—2605

    Schmitt R W, 1981. Form of the temperature-salinity relationship in the central water: Evidence for double-diffusive mixing.Journal of Physical Oceanography, 11(7): 1015—1026

    Schmitt R W, 2003. Observational and laboratory insights into salt finger convection. Progress in oceanography, 56(3):419—433

    Schmitt R W, Ledwell J, Montgomery Eet al, 2005. Enhanced diapycnal mixing by salt fingers in the thermocline of the tropical Atlantic. Science, 308(5722): 685—688

    Shang X D, Liu Q, Xie X Het al, 2014. Source and seasonal variability of internal tides in the southern South China Sea.Deep-Sea Research Part I: Oceanographic Reseach Papers,submitted

    Tian J, Yang Q, Zhao W, 2009. Enhanced diapycnal mixing in the South China Sea. Journal of Physical Oceanography, 39(12):3191—3203

    Turner J S, 1973. Buoyancy effects in fluids. Cambridge:Cambridge University Press, 367

    Wolk F, Yamazaki H, Seuront Let al, 2002. A new free-fall profiler for measuring biophysical microstructure. Journal of Atmospheric and Oceanic Technology, 19(5): 780—793

    Yamazaki H, 1990. Stratified turbulence near a critical dissipation rate. Journal of Physical Oceanography, 20(10): 1583—1598

    Yang Q, Tian J, Zhao Wet al, 2014. Observations of turbulence on the shelf and slope of northern South China Sea. Deep Sea Research Part I: Oceanographic Research Papers, 87: 43—52

    You Y, 2002. A global ocean climatological atlas of the Turner angle: implications for double-diffusion and water-mass structure. Deep Sea Research Part I: Oceanographic Research Papers, 49(11): 2075—2093

    猜你喜歡
    湍流海域南海
    南海明珠
    遺落海域
    北海北、南海南
    黃河之聲(2021年10期)2021-09-18 03:07:18
    埕島海域海上獨(dú)立樁拆除方案探討
    重氣瞬時(shí)泄漏擴(kuò)散的湍流模型驗(yàn)證
    廣東省海域使用統(tǒng)計(jì)分析
    南海的虎斑貝
    南海隨筆
    草原(2016年1期)2016-01-31 21:21:51
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    国产亚洲91精品色在线| 国产精品综合久久久久久久免费| 岛国在线免费视频观看| 亚洲国产高清在线一区二区三| 国产精品野战在线观看| 国产91av在线免费观看| 国产麻豆成人av免费视频| av在线天堂中文字幕| 床上黄色一级片| 亚洲av.av天堂| 美女大奶头视频| 国产三级中文精品| 国产精品99久久久久久久久| 亚洲欧美日韩东京热| 国产在视频线精品| 国产亚洲午夜精品一区二区久久 | 亚洲va在线va天堂va国产| 国产视频首页在线观看| 99热网站在线观看| 哪个播放器可以免费观看大片| 午夜福利高清视频| 亚洲国产高清在线一区二区三| 国产探花极品一区二区| 日本一本二区三区精品| 亚洲av电影不卡..在线观看| 免费黄色在线免费观看| 亚洲四区av| 少妇猛男粗大的猛烈进出视频 | 久久精品国产鲁丝片午夜精品| 韩国av在线不卡| 男插女下体视频免费在线播放| 一级毛片久久久久久久久女| 九九热线精品视视频播放| 在线免费观看的www视频| 在线观看一区二区三区| 免费观看a级毛片全部| 国内少妇人妻偷人精品xxx网站| 黄色日韩在线| 91狼人影院| 久久亚洲国产成人精品v| 国产午夜精品一二区理论片| 啦啦啦观看免费观看视频高清| 观看免费一级毛片| 午夜久久久久精精品| 中文字幕制服av| av线在线观看网站| 中文字幕亚洲精品专区| 极品教师在线视频| 啦啦啦啦在线视频资源| 日韩 亚洲 欧美在线| 中文乱码字字幕精品一区二区三区 | 久久久久久久久中文| 老女人水多毛片| av黄色大香蕉| 国产综合懂色| 尤物成人国产欧美一区二区三区| 久久久久久久国产电影| 国产黄色视频一区二区在线观看 | 亚洲欧美日韩东京热| 午夜免费男女啪啪视频观看| 99热这里只有精品一区| 特大巨黑吊av在线直播| 亚洲精品456在线播放app| 午夜激情欧美在线| 麻豆乱淫一区二区| 亚洲国产色片| 99在线视频只有这里精品首页| 亚洲婷婷狠狠爱综合网| 亚洲av电影在线观看一区二区三区 | 久久综合国产亚洲精品| 69人妻影院| av视频在线观看入口| 国产精品久久视频播放| 性色avwww在线观看| 99久国产av精品国产电影| 亚洲欧美日韩高清专用| 成人性生交大片免费视频hd| 日本欧美国产在线视频| 国产亚洲av嫩草精品影院| 日韩高清综合在线| 精品一区二区三区视频在线| 亚洲人与动物交配视频| 亚洲精品456在线播放app| 成年免费大片在线观看| 色视频www国产| 超碰97精品在线观看| 午夜精品国产一区二区电影 | 国语自产精品视频在线第100页| 成年版毛片免费区| 天堂中文最新版在线下载 | 一边亲一边摸免费视频| 99在线人妻在线中文字幕| 国产亚洲91精品色在线| 我要看日韩黄色一级片| 精品免费久久久久久久清纯| 狠狠狠狠99中文字幕| 欧美成人一区二区免费高清观看| 午夜福利视频1000在线观看| 亚洲欧美清纯卡通| 国产麻豆成人av免费视频| 国产精品蜜桃在线观看| 亚洲四区av| 免费黄色在线免费观看| 久久精品国产亚洲网站| 卡戴珊不雅视频在线播放| 国产精华一区二区三区| 99久国产av精品国产电影| 成人av在线播放网站| 久久99热6这里只有精品| 久久草成人影院| 亚洲av电影不卡..在线观看| 亚洲天堂国产精品一区在线| 秋霞伦理黄片| 91狼人影院| 欧美区成人在线视频| 国产精品一区二区在线观看99 | 中文在线观看免费www的网站| 成人美女网站在线观看视频| 久久久久久久久中文| 中文字幕久久专区| 联通29元200g的流量卡| 国产欧美日韩精品一区二区| 免费av观看视频| 午夜免费男女啪啪视频观看| 免费看光身美女| 国产精品国产高清国产av| 在线天堂最新版资源| 国产爱豆传媒在线观看| 日本与韩国留学比较| 少妇人妻一区二区三区视频| 中文欧美无线码| 美女大奶头视频| 日韩欧美三级三区| 天堂√8在线中文| 高清日韩中文字幕在线| 日韩大片免费观看网站 | 亚洲精品国产成人久久av| 在线观看美女被高潮喷水网站| 嫩草影院新地址| 亚洲四区av| 国产精品国产三级国产专区5o | 婷婷色麻豆天堂久久 | 国产美女午夜福利| 国产亚洲一区二区精品| 99久久人妻综合| 亚洲av成人精品一区久久| 国产亚洲精品av在线| 特级一级黄色大片| 国产精品一区二区在线观看99 | 欧美成人午夜免费资源| 乱码一卡2卡4卡精品| 成人欧美大片| 特级一级黄色大片| 小说图片视频综合网站| 夜夜看夜夜爽夜夜摸| 久久99精品国语久久久| 久久婷婷人人爽人人干人人爱| 国产成人福利小说| 免费看美女性在线毛片视频| 国产精品久久久久久av不卡| 一区二区三区乱码不卡18| 久久精品夜色国产| 直男gayav资源| 国产高清三级在线| 91精品一卡2卡3卡4卡| 在线播放无遮挡| 夫妻性生交免费视频一级片| 国产高潮美女av| 精华霜和精华液先用哪个| 丰满人妻一区二区三区视频av| 久久久久久久午夜电影| 好男人视频免费观看在线| 三级毛片av免费| 精品99又大又爽又粗少妇毛片| 欧美97在线视频| 免费看光身美女| 91久久精品电影网| 热99在线观看视频| 淫秽高清视频在线观看| 亚洲国产精品成人综合色| 人妻少妇偷人精品九色| 国产久久久一区二区三区| 成人午夜精彩视频在线观看| 黄色日韩在线| 天天躁夜夜躁狠狠久久av| 精品久久久久久电影网 | 在线播放无遮挡| 国产精品久久久久久久电影| 亚洲经典国产精华液单| 日本-黄色视频高清免费观看| 国产精品一区二区三区四区免费观看| 少妇熟女aⅴ在线视频| 在线天堂最新版资源| 日本午夜av视频| 精品久久国产蜜桃| 亚洲国产最新在线播放| 久久99热6这里只有精品| 欧美97在线视频| 欧美日韩综合久久久久久| 亚洲精品,欧美精品| 五月玫瑰六月丁香| 一二三四中文在线观看免费高清| 在线免费观看不下载黄p国产| 精品久久久噜噜| 中文字幕av在线有码专区| 免费无遮挡裸体视频| 亚洲激情五月婷婷啪啪| 男人舔奶头视频| 日本与韩国留学比较| 亚洲av成人av| 午夜福利在线观看吧| 建设人人有责人人尽责人人享有的 | 亚洲成人精品中文字幕电影| 色综合亚洲欧美另类图片| 淫秽高清视频在线观看| 一本久久精品| 亚洲欧美精品综合久久99| 欧美潮喷喷水| 波野结衣二区三区在线| 欧美日韩精品成人综合77777| 国产欧美另类精品又又久久亚洲欧美| 久久韩国三级中文字幕| 非洲黑人性xxxx精品又粗又长| 中文在线观看免费www的网站| 天天一区二区日本电影三级| 日日啪夜夜撸| 国产欧美另类精品又又久久亚洲欧美| 欧美zozozo另类| 桃色一区二区三区在线观看| 欧美一级a爱片免费观看看| 国产成人freesex在线| 婷婷色麻豆天堂久久 | 久久精品国产鲁丝片午夜精品| 色5月婷婷丁香| 永久免费av网站大全| 日日撸夜夜添| 国产黄片视频在线免费观看| 最近2019中文字幕mv第一页| av在线老鸭窝| 亚洲成色77777| 91av网一区二区| 最近视频中文字幕2019在线8| 久久国内精品自在自线图片| 亚洲图色成人| 国产探花极品一区二区| 汤姆久久久久久久影院中文字幕 | 欧美又色又爽又黄视频| 国产黄a三级三级三级人| 亚洲va在线va天堂va国产| 久久久久久久久久久丰满| 日本熟妇午夜| 国产三级在线视频| a级毛片免费高清观看在线播放| 午夜免费激情av| 亚洲不卡免费看| 色5月婷婷丁香| 波多野结衣高清无吗| 成年女人永久免费观看视频| 色播亚洲综合网| 亚洲欧美精品综合久久99| 久久亚洲国产成人精品v| 久久热精品热| 亚洲av福利一区| 亚洲va在线va天堂va国产| 国产综合懂色| 成人亚洲欧美一区二区av| 亚洲精品日韩av片在线观看| 欧美性猛交黑人性爽| 亚洲国产色片| 中文字幕久久专区| 国产熟女欧美一区二区| 少妇丰满av| 免费看美女性在线毛片视频| a级毛色黄片| 国产又色又爽无遮挡免| 午夜福利高清视频| 一边亲一边摸免费视频| 国产精品嫩草影院av在线观看| 在线播放国产精品三级| 亚洲最大成人中文| av播播在线观看一区| 国产视频首页在线观看| 欧美xxxx性猛交bbbb| av专区在线播放| 成人毛片a级毛片在线播放| 久久久久久久久久黄片| 18禁裸乳无遮挡免费网站照片| 亚洲不卡免费看| 亚洲经典国产精华液单| 少妇人妻一区二区三区视频| 免费观看性生交大片5| 久久草成人影院| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 听说在线观看完整版免费高清| 免费在线观看成人毛片| 欧美激情国产日韩精品一区| 日韩av不卡免费在线播放| 天堂√8在线中文| 麻豆成人午夜福利视频| 免费人成在线观看视频色| 九九在线视频观看精品| 亚洲av男天堂| 嫩草影院入口| 国产乱人偷精品视频| 精品不卡国产一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品日韩av片在线观看| 深爱激情五月婷婷| 国产亚洲一区二区精品| 欧美人与善性xxx| 中文欧美无线码| 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| 伊人久久精品亚洲午夜| 少妇人妻一区二区三区视频| 18禁在线播放成人免费| 欧美成人精品欧美一级黄| 亚洲,欧美,日韩| 日本午夜av视频| 精品久久久久久久末码| 精品人妻一区二区三区麻豆| 欧美三级亚洲精品| or卡值多少钱| 99久国产av精品| 欧美+日韩+精品| 日本与韩国留学比较| 两个人视频免费观看高清| 狂野欧美白嫩少妇大欣赏| 午夜免费男女啪啪视频观看| 少妇人妻一区二区三区视频| 99久久无色码亚洲精品果冻| 午夜久久久久精精品| 国产精品久久久久久久电影| 亚洲国产欧美在线一区| 好男人在线观看高清免费视频| 免费人成在线观看视频色| 国产极品天堂在线| 日韩成人av中文字幕在线观看| 久久精品国产99精品国产亚洲性色| 日本熟妇午夜| 婷婷六月久久综合丁香| 一本一本综合久久| 国产精品人妻久久久久久| 村上凉子中文字幕在线| 搡老妇女老女人老熟妇| 午夜精品国产一区二区电影 | 国产一区亚洲一区在线观看| 国产淫片久久久久久久久| 久久亚洲国产成人精品v| 日本午夜av视频| 国产亚洲精品久久久com| 成人欧美大片| 黄片无遮挡物在线观看| 亚洲内射少妇av| 七月丁香在线播放| 在线观看av片永久免费下载| 精品人妻一区二区三区麻豆| 国产亚洲最大av| www.av在线官网国产| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 91在线精品国自产拍蜜月| 免费黄网站久久成人精品| 高清午夜精品一区二区三区| 国产一区二区亚洲精品在线观看| 韩国高清视频一区二区三区| 看免费成人av毛片| 色吧在线观看| 成人性生交大片免费视频hd| 波多野结衣巨乳人妻| 亚洲精品成人久久久久久| 国语对白做爰xxxⅹ性视频网站| 99热精品在线国产| 久久国产乱子免费精品| av国产久精品久网站免费入址| 成人综合一区亚洲| 免费看a级黄色片| 99国产精品一区二区蜜桃av| 高清在线视频一区二区三区 | 亚洲自偷自拍三级| 麻豆乱淫一区二区| 网址你懂的国产日韩在线| 精品久久久久久久久av| 亚洲一区高清亚洲精品| h日本视频在线播放| 国产精品一区二区三区四区久久| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费| 婷婷六月久久综合丁香| 国产真实乱freesex| 18+在线观看网站| 在线观看一区二区三区| 亚洲欧美日韩高清专用| 国产成人freesex在线| 中文精品一卡2卡3卡4更新| 又爽又黄无遮挡网站| 欧美日韩一区二区视频在线观看视频在线 | 色综合亚洲欧美另类图片| 熟女人妻精品中文字幕| 国产精品人妻久久久久久| 色噜噜av男人的天堂激情| 国产淫片久久久久久久久| 高清视频免费观看一区二区 | 黑人高潮一二区| 国产精品伦人一区二区| 日本黄色视频三级网站网址| 深爱激情五月婷婷| 国产在视频线在精品| 欧美日本亚洲视频在线播放| 一二三四中文在线观看免费高清| 大又大粗又爽又黄少妇毛片口| 亚洲精品日韩av片在线观看| 日韩亚洲欧美综合| 人妻制服诱惑在线中文字幕| 国产av不卡久久| 韩国高清视频一区二区三区| 亚洲成色77777| or卡值多少钱| 蜜桃亚洲精品一区二区三区| 高清午夜精品一区二区三区| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 亚洲av日韩在线播放| 亚洲国产精品成人综合色| 午夜视频国产福利| 国产精品国产三级国产av玫瑰| 人妻系列 视频| 插逼视频在线观看| 国产亚洲最大av| 国产大屁股一区二区在线视频| 国产三级在线视频| 欧美成人一区二区免费高清观看| 亚洲精华国产精华液的使用体验| 黑人高潮一二区| 中文字幕制服av| 日韩欧美 国产精品| 一二三四中文在线观看免费高清| 国产伦精品一区二区三区四那| 哪个播放器可以免费观看大片| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 国产高清视频在线观看网站| 成人二区视频| 亚洲欧美精品自产自拍| 少妇裸体淫交视频免费看高清| 男女那种视频在线观看| 成人毛片60女人毛片免费| 99国产精品一区二区蜜桃av| 国产午夜精品一二区理论片| av免费在线看不卡| 麻豆精品久久久久久蜜桃| 韩国高清视频一区二区三区| 一个人看的www免费观看视频| 国产国拍精品亚洲av在线观看| 欧美成人午夜免费资源| 久久欧美精品欧美久久欧美| 丰满少妇做爰视频| 亚洲在久久综合| 国产黄a三级三级三级人| 亚洲久久久久久中文字幕| 日产精品乱码卡一卡2卡三| 2021少妇久久久久久久久久久| 少妇熟女aⅴ在线视频| 一级毛片电影观看 | 全区人妻精品视频| 校园人妻丝袜中文字幕| 黄色欧美视频在线观看| 九九久久精品国产亚洲av麻豆| 三级毛片av免费| 精品熟女少妇av免费看| 卡戴珊不雅视频在线播放| kizo精华| 国产av不卡久久| av黄色大香蕉| 精华霜和精华液先用哪个| 非洲黑人性xxxx精品又粗又长| 日本黄色视频三级网站网址| 国产高清有码在线观看视频| 日韩在线高清观看一区二区三区| av免费在线看不卡| 特大巨黑吊av在线直播| 男人舔奶头视频| 日本一二三区视频观看| www.av在线官网国产| 乱人视频在线观看| 国产精品电影一区二区三区| 99久久精品热视频| 免费看日本二区| 91午夜精品亚洲一区二区三区| 一本一本综合久久| 精品国产一区二区三区久久久樱花 | 国产免费男女视频| 一边亲一边摸免费视频| 老女人水多毛片| 国产一区二区在线观看日韩| 青春草国产在线视频| 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 国产精品一区二区三区四区免费观看| 国产精品爽爽va在线观看网站| 久久精品熟女亚洲av麻豆精品 | 在线天堂最新版资源| 性插视频无遮挡在线免费观看| 精品少妇黑人巨大在线播放 | 免费人成在线观看视频色| 美女黄网站色视频| 精品人妻一区二区三区麻豆| 国产伦精品一区二区三区视频9| 国产一级毛片在线| 国产极品天堂在线| 亚洲国产色片| 插逼视频在线观看| 亚洲国产精品成人综合色| 亚洲欧美成人精品一区二区| 国产男人的电影天堂91| 国产大屁股一区二区在线视频| 日本wwww免费看| 国产精品久久电影中文字幕| 亚洲自拍偷在线| 麻豆av噜噜一区二区三区| 永久免费av网站大全| 日韩亚洲欧美综合| 九九在线视频观看精品| 好男人视频免费观看在线| 波野结衣二区三区在线| 一个人免费在线观看电影| 国产精品99久久久久久久久| 青春草亚洲视频在线观看| 中文字幕av在线有码专区| 国产一区二区在线av高清观看| 国产大屁股一区二区在线视频| 久久精品久久久久久噜噜老黄 | 国产精品国产三级专区第一集| 亚洲国产色片| 天堂中文最新版在线下载 | 看黄色毛片网站| 色综合站精品国产| 大又大粗又爽又黄少妇毛片口| 日韩欧美三级三区| 久久人妻av系列| 国产精品一及| 久久人妻av系列| 久久久久久伊人网av| 久久精品国产亚洲av天美| eeuss影院久久| 国产v大片淫在线免费观看| 亚洲自拍偷在线| 国产精品伦人一区二区| 中文欧美无线码| 蜜桃亚洲精品一区二区三区| 国产精品福利在线免费观看| 99热这里只有是精品在线观看| 国产精品国产三级专区第一集| 麻豆av噜噜一区二区三区| 日本与韩国留学比较| 婷婷色av中文字幕| av免费观看日本| 中文字幕免费在线视频6| 亚洲av成人精品一二三区| 亚洲美女搞黄在线观看| 国产精品av视频在线免费观看| 色尼玛亚洲综合影院| 99久久精品热视频| 久久精品国产亚洲av天美| 日韩精品有码人妻一区| 欧美成人午夜免费资源| 日日撸夜夜添| videos熟女内射| 国产精品一及| 亚洲国产精品国产精品| 国产精品一及| 又黄又爽又刺激的免费视频.| 日韩大片免费观看网站 | 最近中文字幕2019免费版| 18禁裸乳无遮挡免费网站照片| 高清午夜精品一区二区三区| 最后的刺客免费高清国语| 欧美日韩一区二区视频在线观看视频在线 | 美女cb高潮喷水在线观看| 大话2 男鬼变身卡| 99在线视频只有这里精品首页| 国产极品精品免费视频能看的| 亚洲欧美日韩高清专用| 国产精品综合久久久久久久免费| 国产精品久久久久久精品电影小说 | 免费av观看视频| 一级毛片久久久久久久久女| 欧美高清成人免费视频www| 国产亚洲最大av| 久久久久久久久久成人| 久久99热这里只频精品6学生 | 日本五十路高清| 超碰97精品在线观看| 日本午夜av视频| 久久久久久久久久久丰满| 久久久久久久久久黄片| 青春草视频在线免费观看| 麻豆国产97在线/欧美| 国产大屁股一区二区在线视频| 黑人高潮一二区| 免费一级毛片在线播放高清视频| 国产高清国产精品国产三级 | 国内精品一区二区在线观看| 国产精品久久久久久av不卡| 成人亚洲欧美一区二区av| 一区二区三区乱码不卡18| 你懂的网址亚洲精品在线观看 | 精品国产一区二区三区久久久樱花 | 久久久成人免费电影| 日本三级黄在线观看| 久久这里只有精品中国| 乱码一卡2卡4卡精品| 精品久久国产蜜桃| 国产免费男女视频|