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

    呂宋海峽西側(cè)海域混合層鹽度低頻變化特征分析

    2014-05-30 08:52:52
    海洋學(xué)研究 2014年1期
    關(guān)鍵詞:呂宋黑潮海峽

    沈 華

    (1.衛(wèi)星海洋環(huán)境動(dòng)力學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,浙江 杭州 310012;2.國(guó)家海洋局 第二海洋研究所,浙江 杭州 310012)

    0 引言

    鹽度是海水含鹽量的一個(gè)標(biāo)度,是海水的重要特性之一,其分布與變化和海洋中發(fā)生的許多現(xiàn)象和過程密切相關(guān),同時(shí),鹽度對(duì)海洋的熱力和動(dòng)力過程有著非常重要的影響。因此,對(duì)海水鹽度的分布及變化規(guī)律的了解有利于更深入地理解海洋系統(tǒng)及海氣相互作用過程。

    目前,關(guān)于南海海域鹽度變化的研究,主要從季節(jié)變化、年際變化及年代際變化3個(gè)方面展開。南海屬于東亞季風(fēng)區(qū),受季風(fēng)影響強(qiáng),以致季節(jié)變化最為顯著。隨著夏季風(fēng)的建立,南海北部和東部的鹽度降低,南部的鹽度增加[1]。在南海北部,由于西太平洋高鹽水的輸入,季節(jié)變化相對(duì)較?。欢谀喜?,由于受季風(fēng)轉(zhuǎn)換的影響較大,季節(jié)變化相對(duì)比較明顯。此外,南海存在多個(gè)高鹽中心,其中主要的1個(gè)出現(xiàn)在南越南灣,存在的時(shí)間能跨越整個(gè)上半年,并在2、3、4月份達(dá)到最強(qiáng),中心鹽度在34.7以上[2]。WYRTKI[3]利用1950—1955年的觀測(cè)資料首次繪制了南海夏季和冬季表層鹽度的分布圖,發(fā)現(xiàn)冬季有一支高鹽水舌從呂宋海峽進(jìn)入,并沿著大陸架一直延伸到越南南部;夏季該高鹽水舌減弱,同時(shí)越南南部有一支低鹽水舌向南海中部延伸。謝華偉[4]利用 WOA01的資料分析了加里曼丹島周邊的低鹽水團(tuán)的季節(jié)變化,其中西北側(cè)的低鹽水其中心位置基本保持在拉讓河河口位置,僅在8—9月及12—翌年1月會(huì)向東北偏移;他還進(jìn)一步指出,低鹽水冬季影響面積大,夏季至秋季影響面積小,且在2—3月和9—10月會(huì)發(fā)生2次鹽度降低的過程。

    在南海鹽度的年際變化方面,謝華偉[4]利用19502005年的OFES模式資料發(fā)現(xiàn)加里曼丹島西北部水域鹽度的年際變化存在2~5a的周期信號(hào);進(jìn)一步研究發(fā)現(xiàn),凈淡水通量與垂向運(yùn)動(dòng)對(duì)該地區(qū)鹽度的年際變化較為重要,而徑流的變化關(guān)系不大,同時(shí)ENSO事件對(duì)低鹽水也有著極其重要的影響。

    在年代際變化方面,南峰[5]利用1990—2010年的WOD以及Argo觀測(cè)資料,繪制了呂宋口附近海域的T-S曲線和散點(diǎn)隨時(shí)間變化的分布圖,發(fā)現(xiàn)2000年后黑潮入侵有減弱的趨勢(shì),隨之而來的,呂宋口西側(cè)海域的鹽度呈下降趨勢(shì)。

    呂宋海峽位于中國(guó)臺(tái)灣島和菲律賓呂宋島之間,是南海與西北太平洋水體交換的主要通道,因此,呂宋海峽鹽度的變化一直是學(xué)術(shù)界十分關(guān)心的問題之一。邱春華 等[6]利用SODA資料分析了年際尺度上黑潮和呂宋海峽鹽度的關(guān)系,指出南海北部深水海域的鹽度場(chǎng)受黑潮的影響較大,在黑潮流量大的年份呂宋海峽處鹽度值較低,反之,在黑潮流量小的年份呂宋海峽處鹽度值較高。但黑潮與南海鹽度低頻變化特征的聯(lián)系仍不清楚。因此,本文采用了時(shí)間比較長(zhǎng)的WOD09剖面觀測(cè)資料以及Argo資料分析南海鹽度的低頻信號(hào),進(jìn)而利用混合層鹽度平衡方程研究了呂宋海峽西側(cè)海域上層鹽度低頻變化特征的形成機(jī)理。

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

    1.1 數(shù)據(jù)

    本文所用的WOD09剖面溫、鹽資料來自NODC(National Oceanographic Data Center,http://www.nodc.noaa.gov/OC5/SELECT/-dbsearch/dbsearch.html),時(shí)間范圍為1984—2010年,空間范圍為整個(gè)南海海域(0°~25°N,100°~125°E)。由于該資料為散點(diǎn)剖面觀測(cè),需要先剔除大于3倍標(biāo)準(zhǔn)方差的數(shù)據(jù)和明顯的奇異值以進(jìn)行質(zhì)量控制。該資料的鹽度誤差范圍為-0.016±0.006。

    Argo剖面溫、鹽資料來自中國(guó)Argo實(shí)時(shí)資料中心(http://www.argo.org.cn/index.html),時(shí)間范圍為2005—2010年,空間范圍覆蓋整個(gè)南海海區(qū)。

    為匹配各種資料,將以上2種資料整合后,利用反距離加權(quán)插值方法,內(nèi)插到經(jīng)緯度1°×1°的網(wǎng)格上并通過空間高斯濾波進(jìn)行平滑,深度則內(nèi)插到標(biāo)準(zhǔn)層上(0,10,20,30,50,75,100,125,150,200,250,300,400,500,600,700,800,900,1 000,1 100,1 200,1 300,1 400,1 500,1 750,2 000,2 500,3 000,3 500,4 0004 5005 000和5 500m 以上2種均為實(shí)測(cè)資料,具有較高的可靠性。

    計(jì)算混合層深度和鹽通量所使用的SODA(Simple Ocean Data Assimilation)資料為來自美國(guó)德州農(nóng)工大學(xué)SODA 數(shù)據(jù)網(wǎng)站(http://soda.tamu.edu/assim/)的月平均資料,該資料由同化分析系統(tǒng)產(chǎn)生。該同化分析系統(tǒng)采用隨機(jī)連續(xù)估計(jì)理論和質(zhì)量控制方法,如卡爾曼濾波、臨近點(diǎn)檢驗(yàn)法、四維變分方法等多種方法來減小誤差,以保證資料的準(zhǔn)確度、可用性和可靠性[7]??臻g分辨率為經(jīng)緯度1°×1°,時(shí)間范圍為1984—2010年,空間范圍為0°~25°N,100°~125°E。

    降水資料來自GPCP(Global Precipitation Climatology Project)[8]的網(wǎng)格化逐月平均數(shù)據(jù),空間分辨率為經(jīng)緯度2.5°×2.5°,為匹配其他資料將其插值到經(jīng)緯度1°×1°的網(wǎng)格上。時(shí)間范圍為1984—2010年,空間范圍為0°~25°N,100°~125°E。

    蒸發(fā)資料來自 WHOI(Woods Hole Oceanographic Institution)的逐月平均數(shù)據(jù),空間分辨率為經(jīng)緯度1°×1°。時(shí)間范圍為1984—2010年,空間范圍為0°~25°N,100°~125°E。

    流場(chǎng)資料采用日本地球模擬器中心的OFES(O-cean General Circulation Model for Earth Simulator)模式結(jié)果產(chǎn)品。OFES模式是基于MOM3模式,為提高準(zhǔn)確度和可靠性,動(dòng)量的水平混合采取一個(gè)阻尼可變的雙諧算子,以控制水平網(wǎng)格間的計(jì)算誤差[9]。該資料覆蓋了全球75°S~75°N的范圍,水平網(wǎng)格分辨率為經(jīng)緯度0.1°×0.1°,為匹配其他資料將其插值到經(jīng)緯度1°×1°的網(wǎng)格上。本文采取的資料為以氣候性平均的NCEP風(fēng)場(chǎng)作為驅(qū)使風(fēng)場(chǎng),模式運(yùn)行50a后取1984—2010年整個(gè)南海海域的月平均數(shù)據(jù)作為氣候態(tài)下的資料,用于分析平流項(xiàng)的年際變化。

    1.2 計(jì)算方法

    本文使用混合層鹽度平衡方程對(duì)鹽度變化進(jìn)行診斷,該方程定義如下[10]:

    式中:E為蒸發(fā)量;P為降水量;S為混合層的鹽度;hm為混合層深度;u代表混合層的流矢量;ΔS表示混合層的鹽度與混合層底的鹽度之差(S-S-h(huán));H是海維塞單位函數(shù),當(dāng)we<0時(shí),H(we)=0,否則H(we)=1;we表示卷夾速度,可用下式表示:式(1 表示混合層內(nèi)鹽度的變化,其主要受到3方面的影響:

    (1)混合層表面的凈淡水通量,由式(1)右邊第1項(xiàng)表示。這部分主要受到海洋表層的降水以及蒸發(fā)作用的影響,與水循環(huán)有關(guān);在近岸河口附近海域,徑流量也要納入凈淡水通量的計(jì)算之中。

    (2)平流引起的鹽度變化,由式(1)右邊第2項(xiàng)表示。由于海洋中不同區(qū)域的鹽度存在差異,平流會(huì)將其他海域的淡水或咸水帶入,從而影響當(dāng)?shù)氐暮K}度。

    (3)混合層深度的變化帶來的卷夾效應(yīng),由式(1)右邊第3項(xiàng)表示。混合層底與下層躍層水體之間的界面變化引起混合層鹽度的變化,即卷夾作用。當(dāng)混合層加深時(shí),鹽度普遍相對(duì)較高的原界面下層的水體會(huì)與原混合層水體混合,使混合層鹽度升高;當(dāng)混合層變淺時(shí),該水體不與混合層的水體進(jìn)行混合,故混合層的鹽度不會(huì)發(fā)生變化。

    由鹽度控制方程可看出,影響混合層鹽度變化的主要因素可分為內(nèi)部調(diào)整(水平以及垂直平流)和外部強(qiáng)迫(降水、蒸發(fā)和徑流構(gòu)成的海面凈淡水通量)兩個(gè)部分。公式中的不同項(xiàng)對(duì)于鹽度平衡都有著各自的影響周期,并且在不同海域各項(xiàng)的影響程度也不同。蒸發(fā)降水強(qiáng)的區(qū)域或河口附近海域,蒸發(fā)降水或徑流處于主導(dǎo)地位,劇烈影響鹽度的變化;而流速比較大的海域,會(huì)形成較強(qiáng)的鹽度梯度,在小規(guī)模范圍內(nèi)會(huì)形成明顯的結(jié)構(gòu)分化;在接近赤道或近岸上升流強(qiáng)烈的區(qū)域,卷夾作用也會(huì)對(duì)局地的鹽度造成較大的影響;當(dāng)鹽度層化非常強(qiáng)烈的時(shí)候,混合層底部的垂向混合就會(huì)起到十分重要的作用,抵消由于河口附近的徑流和陸架上平流造成的鹽度降低。

    式(1)計(jì)算中需要對(duì)整個(gè)鹽度混合層進(jìn)行積分,所以需要先計(jì)算出混合層深度(MLD),確定鹽度混合的深度從而進(jìn)行定量分析。本文使用由鹽度變化場(chǎng)定義的混合層深度,利用狀態(tài)方程計(jì)算得到密度場(chǎng),進(jìn)而得到混合層深度,反映表層鹽度的混合程度。在實(shí)際定義中對(duì)于不同的海域,混合層閥值的選取有不同的定義。本文采用的閥值定義及求解過程參照文獻(xiàn) [11]中使用的方法進(jìn)行計(jì)算:

    (1)由SODA溫、鹽資料利用UNESCO1981海洋狀態(tài)方程算得密度場(chǎng);

    (2)將10m層的密度σ10作為起始參考層;

    (3)計(jì)算閥值Δσ值,閥值定義:相對(duì)參考層(10m)溫度低0.5℃處所算得的位密(P=0),該值與10m層的位密之差,即:

    (4)得到σT-0.5=σ10+Δσ,σT-0.5所在的深度即為混合層的深度D。

    此外,趨勢(shì)系數(shù)Rxt可以用下式計(jì)算[12]:

    式中:{xi}為需要計(jì)算趨勢(shì)系數(shù)的時(shí)間序列,為{xi}的均值,n為{xi}的個(gè)數(shù),= (n+1)/2。Rxt為正值(負(fù)值),則表示{xi}的線性趨勢(shì)是上升(下降)的。

    2 鹽度變化特征與分析

    2.1 南海鹽度數(shù)據(jù)量分布

    圖1為南海東北部(其它海域的數(shù)據(jù)個(gè)數(shù)不足150個(gè),故不顯示)1929—2011年間經(jīng)緯度1°×1°網(wǎng)格混合層深度內(nèi)WOD09和Argo鹽度觀測(cè)數(shù)據(jù)量分布。由圖可知,數(shù)據(jù)多分布在SW—NE走向的斜線上,這主要是由于該線為南海最主要的航線;在南海內(nèi)部主要分布在呂宋海峽及兩側(cè)海域。基于以上觀測(cè)數(shù)據(jù)分布情況以及研究呂宋海峽鹽度在深入了解黑潮與南海的水交換中的重要性,本文將著重研究呂宋海峽西側(cè)海域(18.5°~21.5°N ,117.5°~120.5°E)混合層鹽度的低頻變化特征及機(jī)理,該區(qū)域內(nèi)共包含了871個(gè)Argo剖面觀測(cè)數(shù)據(jù),網(wǎng)格化后的鹽度數(shù)據(jù)均大于250個(gè),最大為416個(gè),平均為327個(gè)。

    圖1 南海混合層內(nèi)WOD09和Argo鹽度數(shù)據(jù)量分布(單位:個(gè))Fig.1 Spatial distribution of WOD09and Argo salinity data in mixed layer of South China Sea(Unit:number)

    2.2 呂宋海峽西側(cè)海域混合層鹽度低頻變化特征

    圖2為1984—2010年呂宋海峽西側(cè)海域混合層鹽度的線性變化趨勢(shì)。區(qū)域內(nèi)的鹽度整體上都呈下降趨勢(shì),其中北部下降最多,能達(dá)到0.030/a,西南部下降最少,僅為0.010/a。

    圖2 1984—2010年呂宋海峽西側(cè)海域混合層鹽度線性變化趨勢(shì)的空間分布Fig.2 Spatial distribution of the trend of mixed layer salinity at the western Luzon Strait during 1984-2010

    為了直觀地了解呂宋海峽西側(cè)海域混合層鹽度的長(zhǎng)期變化趨勢(shì)特征,分析了區(qū)域平均混合層鹽度的異常變化,本文關(guān)注的是低頻變化特征,故2a以下周期的高頻變化已通過高斯濾波方法去除。觀測(cè)數(shù)據(jù)在時(shí)間上較為連續(xù),中間有少量的缺測(cè)(占總數(shù)的11.5%),采用內(nèi)插法將這些數(shù)據(jù)補(bǔ)全(圖3中的灰色柱)。經(jīng)統(tǒng)計(jì),1999年后四季的鹽度剖面?zhèn)€數(shù)均為1999年之前的2倍左右,基本可以排除季節(jié)變化引起的誤差。從圖3中不難發(fā)現(xiàn)鹽度變化的2個(gè)顯著特征:

    (1)圖3a中紅實(shí)線為鹽度趨勢(shì)線,從1984年到2010年,鹽度整體上呈現(xiàn)下降的趨勢(shì),線性趨勢(shì)為-0.020/a;

    (2)以1999年為界,1999年以前的鹽度距平值多數(shù)為正,1999年以后的鹽度距平值則恰好相反,負(fù)值占絕大多數(shù)。其中最大正距平值可達(dá)0.14,平均值為0.05;最大負(fù)距平值為-0.22,平均值為-0.06(圖3a)。圖3b中的鹽度累計(jì)異常曲線也很好地驗(yàn)證了這個(gè)現(xiàn)象,1999年之前鹽度總體呈上升趨勢(shì),特別是1990—1996年這段時(shí)間上升得最明顯;1999年之后鹽度總體呈下降趨勢(shì),2003年以前下降很明顯,但2003—2008年這幾年間有稍許回升,隨后又開始顯著下降。圖3b中紅虛線對(duì)應(yīng)鹽度累計(jì)異常峰值的位置,即1999年。

    圖3 呂宋海峽西側(cè)海域混合層鹽度變化Fig.3 Mixed layer salinity anomaly at the western Luzon Strait

    3 呂宋海峽西側(cè)海域混合層鹽度低頻變化的機(jī)理分析

    3.1 長(zhǎng)期線性趨勢(shì)

    圖4為呂宋海峽以西海域混合層年平均的鹽度和凈淡水通量(WHOI蒸發(fā)量減去GPCP降水量)變化曲線。凈淡水通量整體呈上升趨勢(shì),為0.019cm/a。這表示進(jìn)入海洋的淡水在減少,鹽度變化應(yīng)該是上升趨勢(shì)。但是實(shí)際的鹽度變化呈下降趨勢(shì),下降幅度為-0.020/a,所以我們認(rèn)為呂宋海峽西側(cè)海域混合層的鹽度下降并不是由凈淡水通量引起的。

    圖4 混合層鹽度與凈淡水通量變化和趨勢(shì)Fig.4 Variation and trend of mixed layer salinity and freshwater flux

    圖5顯示了呂宋海峽西側(cè)海域鹽通量的變化(利用SODA流場(chǎng)資料計(jì)算),其中研究區(qū)域東面和南面的鹽通量為正值,表明外界通過這兩側(cè)輸入鹽分;而研究區(qū)域北面和西面的鹽通量為負(fù)值,表明研究海域內(nèi)部通過這兩側(cè)輸出鹽分;研究海域總的鹽通量為負(fù)值,表明總體上研究區(qū)域內(nèi)部的鹽分是流失的。從趨勢(shì)來看,北面鹽通量下降最多,為-0.015a東面的次之,為-0.012/a;而南面的下降得最少,僅為-0.009/a;此外只有西面的鹽通量是呈上升趨勢(shì)的,為0.016/a。呂宋海峽西側(cè)海域總的鹽通量呈下降趨勢(shì),為-0.016/a。由此我們可以推斷,鹽度的長(zhǎng)期下降趨勢(shì)主要是由于研究區(qū)北面和東面的鹽通量下降引起的。

    圖5 呂宋海峽西側(cè)海域鹽通量的變化Fig.5 Variation of salinity flux at the western Luzon Strait

    3.2 年代際變化

    圖6為呂宋海峽西側(cè)海域1984—1998年和1999—2010年2個(gè)時(shí)間段平均的鹽度和凈淡水通量異常分布。從圖中不難看出,1984—1998年的鹽度是正距平,1999—2010年的是負(fù)距平;而凈淡水通量則正好相反,1984—1998年的為負(fù)距平,1999—2010年的為正距平。另外,從峰值的走向來看,鹽度異常呈SE走向,而凈淡水通量異常則呈NE走向。所以我們認(rèn)為鹽度的年代際變化也不是由凈淡水通量引起的。

    圖7為呂宋海峽西側(cè)和東側(cè)海域(122.5°~124.5°,18.5°~21.5°N)的T-S曲線,其中等值線為等位密線。從圖中不難發(fā)現(xiàn),呂宋海峽西側(cè)海域的鹽度,1999年之后比1999年之前有明顯下降;而呂宋海峽東側(cè)海域的鹽度不僅沒有降低,在上層還略微有所增加。因此,呂宋海峽西側(cè)海域鹽度在1999年后相對(duì)較低的現(xiàn)象并不是由黑潮本身海水的變化引起的。

    圖6 凈淡水通量(a、b)和鹽度(c、d)變化率多年平均分布Fig.6 Multi-year mean variation rate distribution of freshwater flux(a、b)and salinity(c、d)

    圖7 呂宋海峽西側(cè)(a)和東側(cè)(b)海域T-S曲線Fig.7 T-Scurve at western Luzon Strait(a)and eastern Luzon Strait(b)

    圖8為呂宋海峽西側(cè)海域的鹽度沿著緯線方向的時(shí)間序列分布,其中圖8a為每個(gè)觀測(cè)剖面的最大鹽度值,圖8b為混合層的鹽度值,等值線是經(jīng)度網(wǎng)格上取2°內(nèi)作年平均后得到的結(jié)果。眾所周知,次表層的高鹽水是能夠較好地表征黑潮入侵南海的示蹤物,且對(duì)混合層的鹽度會(huì)產(chǎn)生不小的影響。從圖8a中可以看出,以1999年為界(圖中黑色橫線),1999年之前黑潮入侵明顯比1999年之后劇烈,從34.68的等鹽線能明顯看出,1999年之前等鹽線相對(duì)入侵得更遠(yuǎn),即跨越121°E經(jīng)線(圖中黑色豎線)的距離更遠(yuǎn)。從圖8b中的34.14等鹽線也可以看出同樣的現(xiàn)象,即混合層的鹽度在1999年之前入侵得更加劇烈。

    圖8 鹽度緯向時(shí)間序列分布Fig.8 Time series distribution of zonal scatter of salinity

    如圖9所示,OFES模式模擬的鹽度與實(shí)際鹽度(WOD鹽度)的符合度相當(dāng)高,可見OFES模式資料適用于研究呂宋海峽海域的鹽度。從模式的流速中得出,黑潮經(jīng)向流速與混合層鹽度的相關(guān)系數(shù)可達(dá)-0.55,緯向流速與鹽度的相關(guān)系數(shù)可達(dá)0.49,且1999年之前經(jīng)向流速以負(fù)距平為主,緯向流速以正距平為主,1999年之后則相反,經(jīng)向流速以正距平為主,緯向流速以負(fù)距平為主,這很好地說明了1999年之前黑潮入侵更為劇烈的事實(shí)。不考慮其他因素的影響,若黑潮沿著經(jīng)線方向(18.5°N斷面)的流速增大(即經(jīng)向流速增大),則黑潮的緯向流速(121°E斷面)減小,即黑潮更不容易入侵南海,反之,若經(jīng)向流速減小,黑潮更容易入侵南海,這也支持了SHEREMET 的觀點(diǎn)[13]。

    圖9 黑潮流速異常與混合層鹽度異常的關(guān)系Fig.9 Relationship between the mixed layer salinity anomaly and Kuroshio velocity anomaly

    考慮到數(shù)據(jù)采集的時(shí)候沒有較大偏差的前提下[5],1999年之前呂宋海峽西側(cè)海域的鹽度為正距平,而1999年之后的為負(fù)距平,這主要可能是由于PDO(Pacific Decadal Oscillation)的年代際變化引起黑潮的徑向流量變化,從而導(dǎo)致黑潮入侵在1999年之后減弱引起的。

    3.3 鹽度平衡方程

    利用WOD和Argo鹽度資料、GPCP降水和WHOI蒸發(fā)資料、OFES流場(chǎng)資料,通過鹽度平衡方程式分別計(jì)算鹽度變化項(xiàng)、凈淡水通量項(xiàng)、平流項(xiàng)和卷夾項(xiàng),并對(duì)呂宋海峽西側(cè)海域取平均,得到的時(shí)間序列如圖10所示。圖10a中黑色點(diǎn)實(shí)線為鹽度變化項(xiàng),黑實(shí)線為等式(1)右邊3項(xiàng)之和,兩者相關(guān)度達(dá)到0.52,兩者數(shù)值有稍許差異,其原因可能是:資料來源不同,數(shù)據(jù)預(yù)處理方法不一致,且模式對(duì)鹽度的模擬還存在一定的誤差。但總體趨勢(shì)較為一致,因此利用混合層鹽度方程來解釋鹽度變化還是可行的。圖10b中黑色點(diǎn)實(shí)線為凈淡水通量項(xiàng),黑色空心點(diǎn)實(shí)線為平流項(xiàng),黑實(shí)線為卷夾項(xiàng),其中平流項(xiàng)和鹽度變化項(xiàng)的關(guān)系達(dá)到0.48,由此可見,鹽度變化主要受平流作用的控制。

    圖10 呂宋海峽西側(cè)海域鹽度診斷結(jié)果Fig.10 Diagnosis of salinity budget at the western Luzon Strait

    4 結(jié)論

    本文利用WOD09和Argo資料計(jì)算得到溫、鹽網(wǎng)格化數(shù)據(jù),同時(shí)利用SODA、OFES、GPCP和 WHOI等資料對(duì)鹽度低頻變化特征進(jìn)行分析,結(jié)果表明:

    1 呂宋海峽西側(cè)海域的鹽度從19842010年整體上呈現(xiàn)下降趨勢(shì),線性趨勢(shì)為-0.020/a,而凈淡水通量整體上呈現(xiàn)上升趨勢(shì),線性趨勢(shì)為0.019cm/a。相比于凈淡水通量,呂宋海峽西側(cè)海域四周的總鹽通量呈現(xiàn)下降趨勢(shì),線性趨勢(shì)為-0.016/a。其中北側(cè)和東側(cè)的鹽通量下降較多,分別為-0.015/a和-0.012/a。因此鹽度長(zhǎng)期下降的趨勢(shì)和凈淡水通量的關(guān)系比較小,主要是由平流作用,特別是黑潮入侵的變化引起的。

    (2)以1999年為界,呂宋海峽西側(cè)海域的鹽度在1999年之前為正距平,在1999年之后為負(fù)距平,其中最大正距平值可達(dá)0.14,平均值為0.05,而最大負(fù)距平值為-0.22,平均值為-0.06。

    (3)通過鹽度平衡方程計(jì)算可知,呂宋海峽西側(cè)海域的鹽度變化主要受到平流項(xiàng)變化的控制,而受凈淡水通量項(xiàng)和卷夾項(xiàng)的作用則相對(duì)較小。

    (References):

    [1]LI Xiu-zhen,LIANG Wei,WEN Zhi-ping,et al.Preliminary study on the salinity characteristics of South China Sea and its response to the summer monsoon[J].Journal of Tropical Oceanography,2011,30(1):29-34.

    李秀珍,梁衛(wèi),溫之平,等.南海鹽度對(duì)南海夏季風(fēng)響應(yīng)的初步分析[J].熱帶海洋學(xué)報(bào),2011,30(1):29-34.

    [2]CHU P C,LI Rong-feng.South China Sea isopycnal-surface circulation[J].Journal of Physical Oceanography,2000,30(9):2 419-2 438.

    [3]WYRTKI K.Physical oceanography of the southeast Asian waters[D].La Jolla,California:The University of California,Scripps Institution of Oceanography,1961.

    [4]XIE Hua-wei.The salinity variation of the waters surrounding the Kalimantan Island and its mechanism analysis[D].Hangzhou:The Second Institute of Oceanography,SOA,2009.

    謝華偉.加里曼丹島周邊海域鹽度變化及機(jī)理分析[D].杭州:國(guó)家海洋局第二海洋研究所,2009.

    [5]NAN Feng.Spatiotemporal evolution of the current-eddy structure southwest of Taiwan[D].Qingdao:Ocean University of China,2012.

    南峰.臺(tái)灣西南部海域流-渦結(jié)構(gòu)及其演變規(guī)律研究[D].青島:中國(guó)海洋大學(xué),2012.

    [6]QIU Chun-h(huán)ua,JIA Ying-lai.Seasonal and inter-annual variations of temperature and salinity in the Northern South China Sea[J].Periodical of Ocean University of China,2009,39(3):375-380.

    邱春華,賈英來.南海北部深水海域溫度以及鹽度的季節(jié)及年際變化特征[J].中國(guó)海洋大學(xué)學(xué)報(bào),2009,39(3):375-380.

    [7]GUO Jing.Analysis of characters and mechanisms of the mixed layer salinity variation in South China Sea[D].Qingdao:The First Institute of Oceanography,SOA,2012.

    郭敬.南?;旌蠈欲}度變異的特征與機(jī)制分析[D].青島:國(guó)家海洋局第一海洋研究所,2012.

    [8]CHEN Ju,SHI Ping,DU Yan.Intercomparison of rainfall products in South China Sea and its adjacent areas[J].Journal of Tropical Oceanography,2004,23(6):40-51.

    陳舉,施平,杜巖.南海及其鄰近地區(qū)幾種常用降雨產(chǎn)品的相互比較[J].熱帶海洋學(xué)報(bào),2004,23(6):40-51.

    [9]MASUMOTO Y,KAGIMOTO T,KOMORI N,et al.A fiftyyear eddy-resolving simulation of the world ocean preliminary outcomes of OFES[J].Journal of the Earth Simulater,2004,1:35-56.

    [10]LI Ren,RISER S C.Seasonal salt budget in the northeast Pacific Ocean[J].J Geophys Res:Oceans,2009,114(C12):doi:10.1029/2009JC005307.

    [11]Clément de Boyer Montégut,MADEC G,F(xiàn)ISCHER A S,et al.Mixed layer depth over the global ocean:An examination of profile data and a profile-based climatology[J].J Geophys Res:O-ceans,2004,109(C12):doi:10.1029/2004JC002378.

    [12]LIN Chuan-lan,SU Ji-lan,XU Bing-rong,et al.Long--term variations of temperature and salinity of the Bohai Sea and their influence on its ecosystem[J].Progress in Oceanography,2001,49(1-4):7-19.

    [13]SHEREMET V A.Hysteresis of a Western Boundary Current Leaping across a Gap[J].J Phys Oceanogr,2001,31:1 247-1 259.

    猜你喜歡
    呂宋黑潮海峽
    譜松弛動(dòng)力降尺度方法及其在呂宋海峽黑潮模擬中的應(yīng)用
    海峽
    心聲歌刊(2021年6期)2021-02-16 01:12:32
    郁南縣河口鎮(zhèn)芒果喜獲豐收
    源流(2020年6期)2020-08-03 01:59:35
    以新聞史視角探析晚清《申報(bào)》中呂宋票廣告鏡像
    傳播力研究(2019年1期)2019-03-20 11:28:14
    黑潮
    海洋世界(2017年1期)2017-02-13 08:31:44
    基于20 a衛(wèi)星高度計(jì)數(shù)據(jù)的黑潮變異特征
    魔鬼海峽
    海峽可以很遠(yuǎn),但是心可以很近
    海峽姐妹(2016年2期)2016-02-27 15:15:43
    東海冷渦對(duì)黑潮熱輸運(yùn)年代際躍變的響應(yīng)
    呂宋海峽上層海洋對(duì)于臺(tái)風(fēng)南瑪都響應(yīng)的觀測(cè)分析與數(shù)值模擬試驗(yàn)
    国产视频内射| 成人欧美大片| av福利片在线观看| 亚洲欧美日韩无卡精品| 五月玫瑰六月丁香| 97热精品久久久久久| 免费看日本二区| 欧美潮喷喷水| 国产一区二区三区视频了| 国产av在哪里看| 男女边吃奶边做爰视频| 国产亚洲av嫩草精品影院| 桃色一区二区三区在线观看| 午夜福利成人在线免费观看| 女的被弄到高潮叫床怎么办 | 久久香蕉精品热| eeuss影院久久| 色综合站精品国产| 一本一本综合久久| 国产 一区 欧美 日韩| 看十八女毛片水多多多| 精华霜和精华液先用哪个| a级毛片a级免费在线| 国产人妻一区二区三区在| 国产色爽女视频免费观看| 欧美最新免费一区二区三区| eeuss影院久久| 校园人妻丝袜中文字幕| 日本黄色片子视频| netflix在线观看网站| 中文字幕av成人在线电影| 99在线人妻在线中文字幕| 国产精华一区二区三区| av黄色大香蕉| h日本视频在线播放| 3wmmmm亚洲av在线观看| 欧美日韩乱码在线| 日韩 亚洲 欧美在线| 女的被弄到高潮叫床怎么办 | 99久久无色码亚洲精品果冻| 午夜爱爱视频在线播放| 国产成人aa在线观看| av黄色大香蕉| 色综合色国产| 日韩欧美国产一区二区入口| 1000部很黄的大片| 国产老妇女一区| 不卡一级毛片| 男人舔女人下体高潮全视频| 亚洲18禁久久av| 99热只有精品国产| 欧美zozozo另类| 久久精品国产99精品国产亚洲性色| 老女人水多毛片| 国内精品一区二区在线观看| 日韩欧美 国产精品| 最好的美女福利视频网| 欧美黑人欧美精品刺激| 免费搜索国产男女视频| 51国产日韩欧美| 免费在线观看日本一区| 国产老妇女一区| 国产蜜桃级精品一区二区三区| 欧美性猛交╳xxx乱大交人| 国产欧美日韩一区二区精品| 99热6这里只有精品| 麻豆av噜噜一区二区三区| 国产色婷婷99| 99热精品在线国产| 麻豆久久精品国产亚洲av| 91av网一区二区| 久久久久久久午夜电影| 嫩草影院新地址| a在线观看视频网站| 国产久久久一区二区三区| 国产精品女同一区二区软件 | 亚洲av中文字字幕乱码综合| 看免费成人av毛片| 在线观看舔阴道视频| 国产视频一区二区在线看| 我的老师免费观看完整版| 亚洲18禁久久av| 日本 欧美在线| 精品99又大又爽又粗少妇毛片 | 日本免费a在线| 午夜精品在线福利| 精品99又大又爽又粗少妇毛片 | 美女xxoo啪啪120秒动态图| 一进一出抽搐动态| 亚洲天堂国产精品一区在线| 老司机深夜福利视频在线观看| 国产极品精品免费视频能看的| 亚洲真实伦在线观看| 国产高清视频在线观看网站| 一个人观看的视频www高清免费观看| 午夜老司机福利剧场| 精品欧美国产一区二区三| 级片在线观看| 搞女人的毛片| 国产极品精品免费视频能看的| 99热网站在线观看| 日日干狠狠操夜夜爽| 亚洲精华国产精华精| 欧美色欧美亚洲另类二区| 少妇人妻精品综合一区二区 | av天堂在线播放| 香蕉av资源在线| 国产av不卡久久| 国产精品精品国产色婷婷| 日韩av在线大香蕉| 亚洲va在线va天堂va国产| 亚洲欧美激情综合另类| 国产单亲对白刺激| 精品一区二区三区av网在线观看| 日本熟妇午夜| 91av网一区二区| 日韩av在线大香蕉| 黄色丝袜av网址大全| 在线观看免费视频日本深夜| 99热只有精品国产| 国产成人aa在线观看| 伦精品一区二区三区| 五月伊人婷婷丁香| 日本免费a在线| 精品99又大又爽又粗少妇毛片 | 一区二区三区免费毛片| 久久国产精品人妻蜜桃| 悠悠久久av| 日韩强制内射视频| 特大巨黑吊av在线直播| 亚洲人成网站在线播放欧美日韩| 五月玫瑰六月丁香| 国产中年淑女户外野战色| 国产伦在线观看视频一区| 国内揄拍国产精品人妻在线| 久久久久九九精品影院| 国产精品野战在线观看| 最好的美女福利视频网| 免费搜索国产男女视频| 日本三级黄在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产国拍精品亚洲av在线观看| 美女 人体艺术 gogo| 小蜜桃在线观看免费完整版高清| 老熟妇仑乱视频hdxx| 乱人视频在线观看| 欧美中文日本在线观看视频| 日韩精品中文字幕看吧| 男女做爰动态图高潮gif福利片| 日本 欧美在线| 又黄又爽又刺激的免费视频.| 成人毛片a级毛片在线播放| 久久久精品欧美日韩精品| 女的被弄到高潮叫床怎么办 | 国产成人av教育| 国产成人一区二区在线| 99久久精品热视频| 久久久色成人| 91在线观看av| 久久精品夜夜夜夜夜久久蜜豆| 桃红色精品国产亚洲av| 看黄色毛片网站| 国产午夜精品久久久久久一区二区三区 | 国产高清不卡午夜福利| 国产欧美日韩精品亚洲av| 久久精品国产亚洲网站| 黄色一级大片看看| 久久精品国产清高在天天线| 久久久色成人| 看片在线看免费视频| 2021天堂中文幕一二区在线观| 免费搜索国产男女视频| 免费观看精品视频网站| 国产精品无大码| 最近最新免费中文字幕在线| 窝窝影院91人妻| 亚洲精品日韩av片在线观看| a在线观看视频网站| 日韩欧美在线乱码| 婷婷六月久久综合丁香| 在现免费观看毛片| 韩国av一区二区三区四区| 国产免费av片在线观看野外av| 久久久久九九精品影院| 久久久久久久亚洲中文字幕| 我的老师免费观看完整版| 国产黄a三级三级三级人| 成人av在线播放网站| 又爽又黄a免费视频| 亚洲av二区三区四区| 神马国产精品三级电影在线观看| 亚洲一级一片aⅴ在线观看| 91在线精品国自产拍蜜月| 中文字幕精品亚洲无线码一区| 久久久久久久久久久丰满 | 精品人妻视频免费看| 女生性感内裤真人,穿戴方法视频| av视频在线观看入口| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 精品久久国产蜜桃| 久久久国产成人免费| 中文在线观看免费www的网站| av福利片在线观看| 久久精品国产清高在天天线| 嫩草影院入口| 亚洲图色成人| 国语自产精品视频在线第100页| 精品久久久久久久久亚洲 | 亚洲色图av天堂| 99热网站在线观看| 男人的好看免费观看在线视频| 真实男女啪啪啪动态图| 动漫黄色视频在线观看| 亚洲av成人av| 桃红色精品国产亚洲av| 天堂网av新在线| 麻豆久久精品国产亚洲av| 高清日韩中文字幕在线| 国产黄色小视频在线观看| 亚洲人成伊人成综合网2020| 亚洲av中文字字幕乱码综合| 一进一出抽搐gif免费好疼| 久久99热这里只有精品18| 超碰av人人做人人爽久久| 日韩 亚洲 欧美在线| 日日摸夜夜添夜夜添小说| 国产女主播在线喷水免费视频网站 | 国产av麻豆久久久久久久| 久久久午夜欧美精品| 中国美女看黄片| 免费av观看视频| 一进一出好大好爽视频| 桃色一区二区三区在线观看| 免费观看在线日韩| 成人午夜高清在线视频| 看片在线看免费视频| www.色视频.com| 韩国av在线不卡| 免费观看在线日韩| 久久久久久久精品吃奶| 人人妻人人看人人澡| 久久精品国产清高在天天线| 又爽又黄无遮挡网站| 免费在线观看影片大全网站| 欧美成人a在线观看| 欧美日本视频| 亚洲av美国av| 亚洲av电影不卡..在线观看| 身体一侧抽搐| 制服丝袜大香蕉在线| 色尼玛亚洲综合影院| 亚洲美女黄片视频| 99久久精品国产国产毛片| 久久久久久久久中文| 日日撸夜夜添| 美女高潮的动态| 亚洲成人免费电影在线观看| 国产伦人伦偷精品视频| 日日夜夜操网爽| 久久久久久伊人网av| 国产精品伦人一区二区| 午夜福利高清视频| 成人无遮挡网站| 欧美人与善性xxx| 亚洲熟妇熟女久久| 欧美高清性xxxxhd video| 日本黄大片高清| 性插视频无遮挡在线免费观看| 色5月婷婷丁香| 久久精品影院6| 国内久久婷婷六月综合欲色啪| 免费看a级黄色片| 午夜精品久久久久久毛片777| 亚洲经典国产精华液单| 免费看美女性在线毛片视频| 啦啦啦韩国在线观看视频| 国产精品av视频在线免费观看| 免费在线观看成人毛片| 国产女主播在线喷水免费视频网站 | 亚洲国产精品sss在线观看| www.www免费av| 亚洲精华国产精华液的使用体验 | 亚洲中文日韩欧美视频| 在线观看av片永久免费下载| а√天堂www在线а√下载| 久久精品国产亚洲av天美| 亚洲 国产 在线| а√天堂www在线а√下载| 美女被艹到高潮喷水动态| 97人妻精品一区二区三区麻豆| 成年人黄色毛片网站| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 97超视频在线观看视频| 国产成人影院久久av| 午夜日韩欧美国产| 国产 一区精品| av视频在线观看入口| 国产在视频线在精品| 熟女电影av网| www.色视频.com| 最近在线观看免费完整版| 波多野结衣巨乳人妻| 国产精品无大码| 99久久无色码亚洲精品果冻| 尤物成人国产欧美一区二区三区| 亚洲中文字幕一区二区三区有码在线看| 亚州av有码| 国产麻豆成人av免费视频| 在线观看舔阴道视频| 2021天堂中文幕一二区在线观| а√天堂www在线а√下载| 男女那种视频在线观看| 嫩草影视91久久| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 久久6这里有精品| 国产伦一二天堂av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 精品久久久噜噜| 亚洲综合色惰| 日韩欧美在线二视频| 黄片wwwwww| 亚洲精品456在线播放app | 亚洲欧美精品综合久久99| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 欧美日韩黄片免| 国产精品久久视频播放| 少妇熟女aⅴ在线视频| 免费观看人在逋| 人人妻人人看人人澡| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 少妇人妻一区二区三区视频| 夜夜爽天天搞| 免费无遮挡裸体视频| а√天堂www在线а√下载| 欧美+日韩+精品| 极品教师在线免费播放| 18+在线观看网站| 九九热线精品视视频播放| 免费一级毛片在线播放高清视频| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 99久久中文字幕三级久久日本| 高清在线国产一区| 国产免费av片在线观看野外av| 内地一区二区视频在线| 欧美日韩精品成人综合77777| 一级av片app| 亚洲自拍偷在线| 韩国av在线不卡| 日韩欧美国产一区二区入口| 亚洲成人久久爱视频| 一区二区三区高清视频在线| 久久午夜亚洲精品久久| 成人特级黄色片久久久久久久| 欧美xxxx黑人xx丫x性爽| 色综合色国产| 久久久久久久精品吃奶| 白带黄色成豆腐渣| 久久久久久久精品吃奶| 色综合色国产| 欧美黑人巨大hd| 黄色日韩在线| 最近最新免费中文字幕在线| 免费av不卡在线播放| 久久亚洲真实| 亚洲色图av天堂| 欧美成人a在线观看| 国产精品1区2区在线观看.| 国产高清有码在线观看视频| 日韩国内少妇激情av| 在线观看美女被高潮喷水网站| 国产伦精品一区二区三区四那| 亚洲av成人av| 天堂网av新在线| 国产成年人精品一区二区| 乱码一卡2卡4卡精品| 国产免费av片在线观看野外av| 亚洲精品乱码久久久v下载方式| 国产一区二区在线观看日韩| 内地一区二区视频在线| 免费av毛片视频| 午夜久久久久精精品| 97超视频在线观看视频| 亚洲三级黄色毛片| 天堂√8在线中文| 中文字幕久久专区| 九色成人免费人妻av| 国产探花在线观看一区二区| 国产黄a三级三级三级人| 真人做人爱边吃奶动态| a级毛片免费高清观看在线播放| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女| 中国美白少妇内射xxxbb| 观看美女的网站| 又黄又爽又免费观看的视频| 欧美在线一区亚洲| 免费av不卡在线播放| 久久久久九九精品影院| av在线蜜桃| 日韩 亚洲 欧美在线| 人人妻,人人澡人人爽秒播| 桃色一区二区三区在线观看| 18禁黄网站禁片免费观看直播| 中文字幕免费在线视频6| 99久久精品热视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲一级一片aⅴ在线观看| 人妻丰满熟妇av一区二区三区| 长腿黑丝高跟| 国产精品野战在线观看| 久久久精品欧美日韩精品| 内射极品少妇av片p| 色播亚洲综合网| 欧美日韩综合久久久久久 | 成人高潮视频无遮挡免费网站| 岛国在线免费视频观看| 在线a可以看的网站| 久久久久九九精品影院| 国内精品美女久久久久久| 别揉我奶头~嗯~啊~动态视频| 久久欧美精品欧美久久欧美| 亚洲成人久久爱视频| 波多野结衣高清无吗| 亚洲成a人片在线一区二区| 亚洲成人中文字幕在线播放| 亚洲男人的天堂狠狠| 亚洲av.av天堂| 亚洲人成网站高清观看| 亚洲精品粉嫩美女一区| 在线观看舔阴道视频| 91在线精品国自产拍蜜月| 黄色日韩在线| 好男人在线观看高清免费视频| 国产乱人视频| 女人被狂操c到高潮| 男人舔女人下体高潮全视频| 老女人水多毛片| 欧美精品国产亚洲| 中文字幕av在线有码专区| 偷拍熟女少妇极品色| 伦理电影大哥的女人| a在线观看视频网站| 成人毛片a级毛片在线播放| 少妇人妻精品综合一区二区 | 久久人妻av系列| 久久精品影院6| 小蜜桃在线观看免费完整版高清| 精品免费久久久久久久清纯| 欧美日韩综合久久久久久 | 夜夜看夜夜爽夜夜摸| 久99久视频精品免费| 亚洲成人精品中文字幕电影| 亚洲男人的天堂狠狠| 亚洲av免费在线观看| 赤兔流量卡办理| 狠狠狠狠99中文字幕| 九九在线视频观看精品| 简卡轻食公司| 搡老岳熟女国产| 国产在视频线在精品| 有码 亚洲区| 特级一级黄色大片| 亚洲五月天丁香| 身体一侧抽搐| 欧美3d第一页| 国产蜜桃级精品一区二区三区| 日韩大尺度精品在线看网址| 噜噜噜噜噜久久久久久91| 国产精品国产高清国产av| 免费观看的影片在线观看| 男女啪啪激烈高潮av片| 久久草成人影院| 国产私拍福利视频在线观看| 国产激情偷乱视频一区二区| 亚洲av免费在线观看| 亚洲图色成人| 久久精品综合一区二区三区| 欧美不卡视频在线免费观看| 1024手机看黄色片| 国产一区二区三区av在线 | 婷婷丁香在线五月| 五月玫瑰六月丁香| 91狼人影院| 成人无遮挡网站| 中文字幕av成人在线电影| 婷婷丁香在线五月| 国产亚洲91精品色在线| 欧美高清成人免费视频www| 免费搜索国产男女视频| 老熟妇仑乱视频hdxx| 欧美3d第一页| 丰满人妻一区二区三区视频av| 国产激情偷乱视频一区二区| 国产高清有码在线观看视频| 身体一侧抽搐| 极品教师在线视频| 十八禁网站免费在线| 人人妻人人看人人澡| 久久亚洲真实| 91午夜精品亚洲一区二区三区 | 少妇的逼好多水| 国产伦在线观看视频一区| 99riav亚洲国产免费| 午夜福利在线在线| 日韩欧美国产在线观看| 黄色视频,在线免费观看| 变态另类成人亚洲欧美熟女| 亚洲精品影视一区二区三区av| 国产老妇女一区| 日本欧美国产在线视频| 日日摸夜夜添夜夜添小说| 婷婷精品国产亚洲av在线| 欧美zozozo另类| 真人做人爱边吃奶动态| 成年女人永久免费观看视频| 成人精品一区二区免费| 久久人人爽人人爽人人片va| 精品人妻一区二区三区麻豆 | 国产精品亚洲一级av第二区| 亚洲欧美日韩卡通动漫| 一个人看视频在线观看www免费| 少妇猛男粗大的猛烈进出视频 | 97超视频在线观看视频| 一卡2卡三卡四卡精品乱码亚洲| www日本黄色视频网| 女的被弄到高潮叫床怎么办 | 免费搜索国产男女视频| 黄色日韩在线| 成年人黄色毛片网站| 日本黄色视频三级网站网址| 成人午夜高清在线视频| 麻豆成人av在线观看| 亚洲精华国产精华精| 国产探花极品一区二区| 女人十人毛片免费观看3o分钟| 亚洲经典国产精华液单| 欧美xxxx黑人xx丫x性爽| 成年女人看的毛片在线观看| 国产精品嫩草影院av在线观看 | 欧美性猛交黑人性爽| 中文字幕av在线有码专区| 亚洲综合色惰| 成人欧美大片| 久久精品综合一区二区三区| 真人做人爱边吃奶动态| 久久久久久大精品| 亚洲精华国产精华液的使用体验 | 国产亚洲精品久久久久久毛片| 久久久久久国产a免费观看| 99热精品在线国产| 国产成年人精品一区二区| 国产 一区 欧美 日韩| 丰满的人妻完整版| 乱码一卡2卡4卡精品| 国产精品一及| 欧洲精品卡2卡3卡4卡5卡区| 一区福利在线观看| 日本三级黄在线观看| 99久久精品热视频| 中文在线观看免费www的网站| 久久久久久伊人网av| 国产aⅴ精品一区二区三区波| 久久人人爽人人爽人人片va| 天堂av国产一区二区熟女人妻| 俄罗斯特黄特色一大片| 国产午夜精品论理片| 丝袜美腿在线中文| 精品不卡国产一区二区三区| 亚洲四区av| 欧美一区二区精品小视频在线| 在线免费观看不下载黄p国产 | 国内精品久久久久久久电影| 色在线成人网| 别揉我奶头 嗯啊视频| 国产成年人精品一区二区| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 99久久中文字幕三级久久日本| 亚洲性久久影院| 1024手机看黄色片| 九九久久精品国产亚洲av麻豆| 国产在线男女| 久久久久国内视频| 亚洲精品在线观看二区| 亚洲在线观看片| av视频在线观看入口| 中文字幕av成人在线电影| 亚洲天堂国产精品一区在线| а√天堂www在线а√下载| 久久久久久九九精品二区国产| 久久久久精品国产欧美久久久| 国产精品自产拍在线观看55亚洲| 亚洲一区二区三区色噜噜| 亚洲avbb在线观看| a级毛片a级免费在线| 999久久久精品免费观看国产| 91麻豆av在线| 日日撸夜夜添| 午夜免费男女啪啪视频观看 | 免费看日本二区| 他把我摸到了高潮在线观看| 免费在线观看日本一区| 免费观看人在逋| 国产高潮美女av| 精品久久久久久成人av| 偷拍熟女少妇极品色| 老熟妇仑乱视频hdxx| 精品人妻一区二区三区麻豆 | 日韩亚洲欧美综合|