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

    孟加拉灣海表鹽度變化特征的衛(wèi)星遙感應(yīng)用研究

    2016-06-01 06:57:42林新宇邱云杜巖周喜武宣莉莉許金電
    海洋學(xué)報(bào) 2016年5期
    關(guān)鍵詞:孟加拉灣環(huán)流

    林新宇,邱云 *,杜巖,周喜武,宣莉莉,許金電

    (1.國家海洋局第三海洋研究所,福建廈門361005;2.中國科學(xué)院南海海洋研究所熱帶海洋環(huán)境國家重點(diǎn)實(shí)驗(yàn)室,廣東廣州510301;3.國家海洋信息中心,天津300171)

    ?

    孟加拉灣海表鹽度變化特征的衛(wèi)星遙感應(yīng)用研究

    林新宇1,邱云1 *,杜巖2,周喜武1,宣莉莉3,許金電1

    (1.國家海洋局第三海洋研究所,福建廈門361005;2.中國科學(xué)院南海海洋研究所熱帶海洋環(huán)境國家重點(diǎn)實(shí)驗(yàn)室,廣東廣州510301;3.國家海洋信息中心,天津300171)

    摘要:本文利用2011年8月至2014年3月Aquarius衛(wèi)星鹽度產(chǎn)品結(jié)合Argo等實(shí)測鹽度資料,探討了孟加拉灣海表鹽度的季節(jié)及年際變化特征。結(jié)果顯示,Aquarius與Argo鹽度呈顯著線性正相關(guān),總體較Argo鹽度值低,偏差為-0.13,其中在孟加拉灣北部海域負(fù)偏差值比南部海域更大,分別為-0.28和-0.10。Aquarius衛(wèi)星與Argo浮標(biāo)在表層鹽度觀測深度上的差別是造成此系統(tǒng)偏差的主因。Aquarius鹽度資料清晰顯示了孟加拉灣海表鹽度具有明顯的季節(jié)變化特征,包括阿拉伯海高鹽水的入侵引起灣南部海域鹽度的變化以及灣北部淡水羽分布范圍的季節(jié)性遷移等主要特征。此外,分析還揭示了2011(2012)年春季整個(gè)灣內(nèi)出現(xiàn)異常高鹽(低鹽)現(xiàn)象。研究表明,2010(2011)年灣北部夏季降雨減少(增加)導(dǎo)致該海域海水鹽度偏高(偏低),并通過表層環(huán)流向南輸運(yùn)引起次年春季灣內(nèi)表層鹽度出現(xiàn)異常高鹽(低鹽)現(xiàn)象,春季風(fēng)應(yīng)力旋度正(負(fù))距平通過影響鹽度垂直混合過程對同期表層鹽度異常高鹽(低鹽)變化也有影響。

    關(guān)鍵詞:孟加拉灣;Aquarius衛(wèi)星;表層鹽度;環(huán)流

    1 引言

    孟加拉灣位于東印度洋北部(圖1),是受南亞季風(fēng)影響的典型海域。南亞夏季風(fēng)給孟加拉灣帶來了豐沛的降水,而且大大超過蒸發(fā)作用,加之灣周邊諸多地表徑流的注入,整個(gè)海域年平均淡水凈收支高達(dá)4.05×1012m3[1]。大量淡水的輸入,使得孟加拉灣成為印度洋鹽度最低的海域,其上層鹽度的變化不僅影響了灣內(nèi)環(huán)流結(jié)構(gòu)[2],還通過引發(fā)障礙層現(xiàn)象[3—4]對南亞季風(fēng)的爆發(fā)與變異有重要的調(diào)制作用[5]。

    已有研究主要針對孟加拉灣表層鹽度的季節(jié)變化規(guī)律,且主要來自于少數(shù)現(xiàn)場水文觀測或再分析資料及數(shù)值模擬結(jié)果[6—8]。早期的研究[6]顯示,受降水、陸地徑流及阿拉伯海高鹽水等因素的影響,孟加拉灣表層鹽度分布呈北低南高的態(tài)勢,灣頂鹽度遠(yuǎn)低于灣口(圖1),在夏季風(fēng)期間二者鹽度差達(dá)到最大[9]。W u等[10]利用鹽度再分析資料給出了灣內(nèi)淡水羽的季節(jié)分布特征,并指出表層環(huán)流的平流輸運(yùn)過程決定了淡水羽的分布范圍。近期的研究[11]認(rèn)為夏季風(fēng)期西南季風(fēng)漂流帶來的次表層阿拉伯海高鹽水經(jīng)Ekman抽吸過程進(jìn)入混合層是維持灣內(nèi)混合層鹽量收支平衡的一個(gè)重要機(jī)制。然而,受制于鹽度資料的限制,對孟加拉灣鹽度年際變化的研究很少[12—13]。

    當(dāng)前,衛(wèi)星鹽度遙感技術(shù)的成熟使得監(jiān)測海洋表層鹽度的變化成為可能,已有的研究[14—15]表明其數(shù)據(jù)質(zhì)量是可靠的。因此,本文嘗試?yán)肁-quarius/S A C-D衛(wèi)星的鹽度遙感資料在孟加拉灣開展海表鹽度變化特征的應(yīng)用研究。為了便于說明,下文的春、夏、秋、冬4個(gè)季節(jié)分別對應(yīng)冬夏季風(fēng)過渡期(3-4月)、夏季季風(fēng)期(5-9月),夏冬季風(fēng)過渡期(10月)以及冬季季風(fēng)期(11月-翌年2月)。

    圖1 孟加拉灣海底地形(等值線,單位:m),圖中疊加了氣候態(tài)表層年平均鹽度分布圖(填色圖),鹽度數(shù)據(jù)來自美國國家海洋數(shù)據(jù)中心世界海洋圖集2013版(W O A 2013)Fig.1 Topography of the Bay of Bengal(B O B,contour in m)with climatological surface-layer salinity(shaded)superimposed.The salinity data comes from W orld Ocean Atlas 2013(W O A 2013)provided by A merican National Oceanographic Data Center

    2 數(shù)據(jù)

    海表鹽度遙感資料來自美國國家航空航天局噴氣推進(jìn)實(shí)驗(yàn)室(N ASA's Jet Propulsion Laboratory,JPL)發(fā)布的Aquarius 2.0版產(chǎn)品。Aquarius/SA C-D衛(wèi)星于2011年6月進(jìn)入太空,其所搭載的Aquarius傳感器是由N ASA研制,通過被動(dòng)式L波段微波輻射計(jì)及主動(dòng)式L波段散射計(jì)觀測的海表面亮溫及海表面風(fēng)速等要素反演得到海表鹽度[16]。文中使用的是L3級A DPS(Aquarius Data Processing System)月平均產(chǎn)品,數(shù)據(jù)已經(jīng)過大氣校正、海表鹽度定標(biāo)及海面粗糙度修正等。鹽度月平均產(chǎn)品空間分辨率為1° ×1°,時(shí)間跨度為2011年8月至2014年3月。

    Argo月平均格點(diǎn)數(shù)據(jù)來自美國斯克里普斯海洋研究所(Scripps Institution of Oceanography)[17],時(shí)間跨度為2004年1月至2013年12月,水平空間分辨率為1°×1°,垂向上共有58層,本文采用其表層(2.5 m層)的鹽度數(shù)據(jù)。Argo鹽度剖面數(shù)據(jù)由法國海洋開發(fā)研究院(French Research Institute for Exploitation of the Sea,IFR E M E R)提供,時(shí)間跨度為2011年8月至2014年3月。本文還使用了“非洲—亞洲—澳洲季風(fēng)分析與預(yù)測研究錨定觀測陣列”(Research M oored Array for African-Asian-Australian M onsoon Analysis and Prediction,R A M A)R1浮標(biāo)獲取的20 m層月平均鹽度時(shí)間序列,該浮標(biāo)位于灣中部海域(15°N,90°E,圖1),觀測時(shí)段為2007年12月至2014年8月,只有20 m層有比較連續(xù)的鹽度資料。上述鹽度觀測資料用于檢驗(yàn)研究海域鹽度遙感數(shù)據(jù)的質(zhì)量,并用于分析鹽度的變化特征。

    海面風(fēng)場資料采用法國海洋開發(fā)研究院提供的ASC AT(the Advanced Scatterometer)月平均風(fēng)場數(shù)據(jù),空間分辨率為0.25°×0.25°,時(shí)間覆蓋范圍為2007年5月至2014年4月。用于計(jì)算研究海域海面淡水通量的降雨及蒸發(fā)資料分別來自T R M M衛(wèi)星(Tropical Rainfall M easuring Mission)月平均3B43-V7版降雨資料[18]及美國伍茲霍爾海洋研究所發(fā)布的海氣通量可觀分析資料集(Objectively Analyzed air-sea Fluxes,O Aflux)[19],其水平分辨率分別為0.25°×0.25°及1°×1°,時(shí)間覆蓋范圍均為2004年1月至2013年12月。河流流域的空間信息(圖1)來自日本東京大學(xué)發(fā)布的T RIP(Total Runoff Integrating Pathways)[20]全球河流流域信息數(shù)據(jù)集。

    3 衛(wèi)星遙感鹽度產(chǎn)品評估

    圖2a給出遙感海表鹽度對Argo鹽度的散點(diǎn)分布。圖中顯示,除低鹽端(鹽度小于32)的部分點(diǎn)外,各數(shù)據(jù)點(diǎn)基本沿y = x周圍呈軸對稱分布,且越靠近對稱軸,樣本密度越大,意味著遙感與Argo實(shí)測鹽度總體吻合較好。從統(tǒng)計(jì)結(jié)果看,二者呈顯著線性正相關(guān),相關(guān)系數(shù)(R)為0.93,均方根誤差(R M SE)為0.48,二者間存在的負(fù)偏差值(-0.13)說明研究海域Aquarius遙感鹽度比Argo鹽度低,而且負(fù)偏差值在低鹽端似乎更大(圖2a)。圖2b為資料點(diǎn)的分布情況,圖2c、d分別為遙感海表鹽度與Argo的平均偏差以及均方根誤差。從平均偏差分布可以看出,在印度半島東西兩側(cè)近岸海域及灣口遙感鹽度略高于Argo鹽度,其余大部分海域遙感鹽度較Argo低,且負(fù)偏差值較大的點(diǎn)主要位于灣北部(15°N以北),在灣頂區(qū)域兩者的鹽度平均偏差可高達(dá)1(圖2c),相應(yīng)地,均方根誤差在灣北部也較大(圖2d)。為進(jìn)一步說明負(fù)偏差隨鹽度變化的特征,表1給出各鹽度值區(qū)間內(nèi)衛(wèi)星與Argo鹽度的比較結(jié)果,可見遙感鹽度在各鹽度值區(qū)間內(nèi)均較Argo鹽度值低,并以鹽度值32為界,在低鹽端(高鹽端)負(fù)偏差值均較大(?。?,為-0.27~-0.29 (-0.08~-0.15),相應(yīng)地,二者的均方差在低鹽端也遠(yuǎn)大于高鹽端,說明在表層鹽度較高的海域這兩種數(shù)據(jù)吻合更好(系統(tǒng)偏差值為-0.10),但在低鹽區(qū)遙感鹽度顯著低于Argo鹽度(系統(tǒng)偏差值為-0.28)。

    圖2 Aquarius與Argo表層鹽度匹配結(jié)果(a),顏色代表以每個(gè)配對點(diǎn)為中心1°×1°方格內(nèi)的配對點(diǎn)個(gè)數(shù),其中Argo只選取有觀測樣本的格點(diǎn)鹽度數(shù)據(jù);匹配點(diǎn)的分布情況(b),顏色代表對應(yīng)網(wǎng)格的配對點(diǎn)個(gè)數(shù);Aquarius與Argo的平均偏差(c);Aquarius與Argo的均方根誤差(d)Fig.2 The collocated data pairs of Aquarius and Argo sea surface salinity(SSS)(a),colors represent the number of dots in the 1°×1°box with its center located at every data pair.Spatial distribution of collocated data pairs(b),colors represent the number of collocated data pairs.M ean deviation between Aquarius and Argo(c).Root mean square error between Aquarius and Argo(d)

    表1 Aquarius海表鹽度與Argo表層鹽度比較的統(tǒng)計(jì)結(jié)果Tab.1 The statistics of the comparison between Aquarius and Argo SSS

    與上述分析結(jié)果一致,已有研究[21—22]顯示A-quarius遙感鹽度在熱帶海域與實(shí)測鹽度資料或Argo觀測鹽度相比存在系統(tǒng)性的偏低現(xiàn)象。Tang等[21]研究指出Aquarius觀測的是海表面鹽度,而Argo大部分剖面只觀測到5 m以深的鹽度,于是海表與5 m層之間的鹽度差別構(gòu)成了Aquarius與Argo鹽度之間系統(tǒng)偏差的一個(gè)來源,并且他們強(qiáng)調(diào)在熱帶海域近表層(5 m以淺)強(qiáng)鹽度層化情形下,會(huì)出現(xiàn)較大的負(fù)偏差。就孟加拉灣而言,淡水(降水和陸地徑流)主要從灣北部注入,導(dǎo)致灣北部近表層鹽度層化最為顯著[6]。從圖1看,低鹽水(鹽度小于32)分布區(qū)域?qū)?yīng)于灣北部的強(qiáng)鹽度層化區(qū),因此不難理解在孟加拉灣低鹽區(qū)遙感海表鹽度較Argo鹽度顯著偏低的現(xiàn)象(圖2,表1)。此外,由于孟加拉灣低鹽水比較貼近岸邊,故上述低鹽區(qū)內(nèi)遙感鹽度顯著偏低還可能與近岸的強(qiáng)射頻干擾(radio frequency interference,RFI)的影響有關(guān)[14]。

    需要指出的是,Aquarius衛(wèi)星的反演目標(biāo)為在時(shí)間分辨率為30 d,考慮所有傳感器和物理模型的隨機(jī)偏差的前提下,鹽度觀測預(yù)期精度達(dá)0.20[23]。但實(shí)際上,衛(wèi)星鹽度觀測精度在不同海區(qū)之間差異較大,在中低緯度大洋的大部分海區(qū)精度較高(精度為0.2 ~0.3),而在近岸海域、強(qiáng)降雨區(qū)及邊緣海鹽度精度較低[21,24—25],比如在南海Aquarius鹽度與實(shí)測鹽度的均方根誤差高達(dá)0.62[25]。類似于南海,孟加拉灣氣象條件異常復(fù)雜,顯著影響了衛(wèi)星鹽度的反演精度,從研究結(jié)果看,Aquarius鹽度在該海域的精度與在南海的大體相當(dāng)。

    4 衛(wèi)星觀測的表層鹽度季節(jié)變化特征

    圖3給出了衛(wèi)星觀測期間1、4、7、10月(分別代表冬、春、夏及秋季)的衛(wèi)星和Argo觀測的鹽度空間分布及二者的差異。圖中顯示二者鹽度差異時(shí)空分布的特點(diǎn)有:(1)Aquarius鹽度在大部分海域較Argo偏低;(2)大致15°N為界,15°N以南海域Aquarius與Argo差值較小,基本界于±0.3之間,局部海域因Argo樣本數(shù)不足,差值稍大,15°N以北海域二者差值較大,并以負(fù)偏差為主,尤其是10月負(fù)偏差最為顯著,但1月在灣西北部Aquarius鹽度明顯高于Argo,可能是因?yàn)樵摵S駻rgo樣本數(shù)不足,導(dǎo)致其網(wǎng)格化的鹽度缺乏代表性。

    雖然上述結(jié)果表明灣內(nèi)部分海域Aquarius鹽度與Argo差值較大,但二者所揭示的鹽度季節(jié)變化特征基本一致(圖3)。具體而言,灣內(nèi)各月海表鹽度分布均總體呈北低南高的態(tài)勢。7月和10月灣南部海域鹽度較高,可見一高鹽水舌(鹽度值為34左右)從阿拉伯海南部經(jīng)斯里蘭卡島東部向?yàn)硟?nèi)延伸。該高鹽水一直持續(xù)至冬季(1月),反映的是西南季風(fēng)漂流攜帶的阿拉伯海高鹽水對灣南部的影響[26]。灣北部鹽度分布的變化顯示了淡水羽(鹽度小于等于32)位置的季節(jié)性遷移。7月受降水和陸地徑流的影響,淡水羽首先在灣頂部出現(xiàn)。從Aquarius鹽度看西部淡水羽更為明顯,隨后淡水羽影響范圍沿著灣東、西岸逐漸向南擴(kuò)大,至10月份影響范圍達(dá)到最大。之后1 月Aquarius鹽度顯示灣西北部淡水羽范圍開始縮小,至4月西北部淡水羽已不明顯,但東北部淡水羽仍存在。不同于Aquarius結(jié)果,Argo鹽度顯示1月灣西北部淡水羽依然清晰可見,且影響范圍甚至比10月還大??紤]到Argo在該海域樣本數(shù)不足,因此其結(jié)果可能無法反映實(shí)際情況。Aquarius鹽度所揭示灣北部淡水羽分布范圍的季節(jié)性變化與W u等[10]的結(jié)果一致,作者并認(rèn)為灣北部表層環(huán)流對淡水羽的演變起重要作用。

    5 衛(wèi)星觀測的表層鹽度年際變化特征

    圖4給出R1浮標(biāo)以及浮標(biāo)位置Aquarius、Argo觀測的逐月鹽度時(shí)間序列。與前文結(jié)果一致,Aquarius鹽度在大部分時(shí)段略低于Argo。由于R1浮標(biāo)觀測的是20 m層的鹽度,因而其大部分鹽度值明顯高于Argo和Aquarius。Aquarius所顯示的鹽度季節(jié)變化特征與其余二者較為一致,均大致表現(xiàn)為西南季風(fēng)期(5-9月)鹽度較高,東北季風(fēng)期(10月-翌年4月)鹽度較低。除季節(jié)變化外,Aquarius還顯示浮標(biāo)位置春季(3-4月)鹽度有顯著的年際差異:2012年春季浮標(biāo)位置的鹽度值較低,僅為31.04,明顯低于2013年的32.25和2014年的32.27。Argo和R1浮標(biāo)也觀測到此異常低鹽現(xiàn)象,鹽度值分別為30.98和31.23,均為各自觀測時(shí)段內(nèi)同期最低,比各自同期平均值(各為32.08和32.33)均分別低1.10。Argo時(shí)間序列還顯示2011年春季鹽度(32.73)為歷年最高,比其平均值高0.65。

    圖3 衛(wèi)星觀測期間(2011年8月至2014年3月)孟加拉灣Aquarius(等值線)及Argo(填色)鹽度合成(a-d)以及Aquarius與Argo鹽度差異(e-h)Fig.3 M onthly SSS of Aquarius(contour)and Argo(shaded)(a-d)and their differences(e-h)over August 2011 to M arch 2014

    圖5為各年春季鹽度距平場,用于進(jìn)一步說明鹽度異常變化的空間分布特征。其中Aquarius鹽度距平是參照Zeng等[14]的做法,將Aquarius逐月鹽度值扣除其與Argo鹽度系統(tǒng)偏差值(-0.13,圖2)后,再與Argo對應(yīng)月份多年均值相減得到。Aquarius結(jié)果顯示2012年春季,除灣西北部近岸局部海域外,幾乎整個(gè)灣出現(xiàn)異常低鹽現(xiàn)象,尤其是在灣中北部(12°~18°N,80°~92°E)鹽度負(fù)距平值大都低于-0.8,最低可達(dá)-1.0。Argo證實(shí)了此異常低鹽現(xiàn)象,并為自有Argo觀測以來春季鹽度最低的年份。與2012年恰好相反,2011年為Argo觀測時(shí)段內(nèi)春季鹽度最高的年份,幾乎整個(gè)灣均出現(xiàn)異常高鹽信號,且也以灣中北部最為突出。雖然2010年異常高鹽信號也較顯著,但因Argo觀測樣本數(shù)太少,較難進(jìn)行相互驗(yàn)證。從Aquarius及Argo鹽度的年際方差分布(圖5k)可看出,在衛(wèi)星觀測期間(2011年8月至2014年3月),兩種資料鹽度年際變化空間分布特征大體一致,灣北部鹽度年際變化比灣南部顯著得多。對比2004-2013年及衛(wèi)星觀測期間Argo海表鹽度的方差(圖5k、5l)顯示,衛(wèi)星觀測期間灣北部的鹽度年際變化方差更大,體現(xiàn)了2012年灣北部異常低鹽現(xiàn)象對該時(shí)間段內(nèi)鹽度年際變化方差的貢獻(xiàn)。

    圖4 R1浮標(biāo)位置(15°N,90°E)的月平均鹽度時(shí)間序列Fig.4 Time series of monthly salinity at the location of R1 float(15°N,90°E)

    6 2011、2012年春季表層鹽度異常變化的成因

    上文分析業(yè)已表明,Aquarius衛(wèi)星、Argo和R1浮標(biāo)均觀測到2012年春季灣內(nèi)出現(xiàn)異常低鹽現(xiàn)象。春季低鹽水是引起同期海溫及夏季風(fēng)爆發(fā)的重要因素[5],因此對春季低鹽水的年際變化的研究有重要的意義。為探討其成因,選取春季鹽度異常高的年份(即2011年)進(jìn)行對比分析。由于孟加拉灣表層鹽度變化主要與海面淡水通量(包括蒸發(fā)與降水過程)、入海河流徑流以及環(huán)流的輸運(yùn)作用(包括水平平流及垂直卷挾過程)等因素有關(guān)[6,10—11],因此本文通過比較這些強(qiáng)迫因子在2011年和2012年的異同,找出引起春季表層鹽度年際差別的可能原因。

    由圖6a和圖6b可見,2011(2012)年春季,幾乎整個(gè)灣的降雨較往年增多(減少),相應(yīng)地,海面凈淡水通量較往年增加(減少),會(huì)導(dǎo)致表層鹽度降低(升高)。顯然,同期海面凈淡水通量的變化均難以解釋這兩年春季表層鹽度的實(shí)際變化情況(圖4,圖5)。此外,從灣頂主要流域降雨看(圖6a、b),2011(2012)年春季較往年減少(增加)。參照前人研究[14,27],河流流域降雨的積分近似為該河口入海徑流量,依此方法計(jì)算得到的這兩年春季灣北部入海徑流量分別較往年減少和增加,也就分別有利于灣北部海域鹽度增加和減少。然而,同期表層鹽度變化顯著區(qū)位于灣中部海域(圖5),因而河口徑流的變化也難以解釋這兩年春季鹽度的差別。

    從前文的結(jié)果看,同期海面淡水通量及河口徑流均無法解釋2011年和2012年春季鹽度的異常變化,這也就意味著鹽度平流輸運(yùn)過程可能為關(guān)鍵原因。由圖7可見,Aquarius鹽度距平分布清晰顯示2012年春季灣南部的低鹽距平信號來自于灣北部,此低鹽距平從2011年夏季(8月)出現(xiàn)在灣北部(16°N以北)并持續(xù)至2012年1月前后,然后向南輸運(yùn),至4月延伸至12°N以南海域,導(dǎo)致春季整個(gè)灣出現(xiàn)異常低鹽的現(xiàn)象。Argo鹽度分布也有大體一致的現(xiàn)象,還顯示了2011年春季高鹽距平信號與2012年相似的擴(kuò)展過程,亦即該高鹽信號源于2010年夏季灣北部海域并在隨后1月向南輸運(yùn)而形成的。很顯然,2010年(2011年)夏季灣北部的高(低)鹽距平是由同期該區(qū)域降水異常變化引起的(圖6c和圖6d)。

    Pant等[28]也注意到了2011年夏季至冬季期間(2011年6月至2012年2月)表層鹽度的淡化現(xiàn)象,并利用鹽度收支方程進(jìn)行了鹽度平衡診斷分析(見其文之圖7),他們分析結(jié)果指出正IO D事件引起6-9月的降雨增多導(dǎo)致灣北部表層鹽度降低,同期夏季風(fēng)環(huán)流以及秋季期間受正IO D事件影響而減弱甚至消失的東印度沿岸流使得低鹽水一直滯留在灣北部,之后冬季風(fēng)期灣內(nèi)反氣旋環(huán)流東支將灣北部低鹽水向南部輸運(yùn),使得冬季灣內(nèi)南部表層鹽度降低。從我們的結(jié)果看(圖5和圖7),此次(2011年)正偶極子事件引起灣北部表層鹽度大幅降低通過環(huán)流輸運(yùn)過程的影響可從冬季一直延續(xù)至翌年(2012年)春季,并導(dǎo)致春季灣內(nèi)出現(xiàn)海盆尺度的低鹽異常現(xiàn)象;類似地,2010年負(fù)偶極子事件期間,灣北部夏季降水顯著減少,導(dǎo)致灣北部表層鹽度升高并同樣通過鹽度平流過程引起翌年(2011年)春季幾乎整個(gè)灣表層鹽度升高。

    圖5 2004-2013年孟加拉灣春季(3-4月)Aquarius(等值線)及Argo(填色)鹽度距平(a-j);衛(wèi)星觀測期間(2011年8月至2014年3月)Argo(填色)和Aquarius(等值線)的年際標(biāo)準(zhǔn)差分布(k);2004-2013年Argo(填色)鹽度年際標(biāo)準(zhǔn)差分布(l)Fig.5 Spring(M arch to April)SSS anomaly from Aquarius(contour)and Argo(shaded)in the B O B during 2004 to 2013(a-j),square deviation of Argo SSS(shaded)and Aquarius SSS(contour)over August 2011 to M arch 2014(k),square deviation of Argo SSS(shaded)during 2004 to 2013(l)

    圖6 孟加拉灣降雨(填色,單位:m m/d)及海面凈淡水通量(蒸發(fā)量-降雨量,即E-P,等值線,單位:m m/d)距平空間分布Fig.6 Spacial distribution of precipitation anomalies(shaded,unit:m m/d)and net fresh water flux anomalies (evaporation minus precipitation,E-P,contour,unit:m m/d)in the B O B

    圖7 85.5°~90.5°E之間平均Argo(填色)及Aquarius(等值線)表層鹽度距平的時(shí)間—緯度分布圖Fig.7 Time-latitude plots of SSS anomaly from Aquarius(contour)and from Argo(shaded)averaged over 85.5°~90.5°E for the period 2010-2013

    參照Pant等[28]的做法,使用Argo鹽度格點(diǎn)數(shù)據(jù)、T R M M 3B43-V7逐月降雨資料、O A FL U X蒸發(fā)資料及G E K C O海流資料利用鹽度收支方程進(jìn)行了鹽度平衡診斷,結(jié)果與Pant等類似(圖略),顯示灣北部的異常高鹽海水在2011年1月份通過平流作用向南延伸形成2011年3-4月的灣內(nèi)大范圍的異常高鹽現(xiàn)象。2012年過程同2011年,但為異常低鹽信號。

    除了上述鹽度水平平流輸運(yùn)過程外,有研究指出春季海洋垂直混合過程對表層鹽度的有重要影響[8]。在多年平均意義下,孟加拉灣春季為海盆尺度反氣旋風(fēng)場的控制[29]。2011年春季灣內(nèi)以東北風(fēng)或東風(fēng)距平為主,并大體呈氣旋式旋轉(zhuǎn),幾乎整個(gè)灣風(fēng)應(yīng)力旋度距平為正值(圖8a),其引起的向上Ek man抽吸將次表層高鹽水卷入混合層有利于混合層鹽度升高,因而對該季節(jié)表層鹽度升高可能有貢獻(xiàn)。2012年海面風(fēng)場及風(fēng)應(yīng)力旋度距平的特征與2011年大體相反,顯示了幾乎整個(gè)灣受負(fù)風(fēng)應(yīng)力旋度距平控制(圖8b),會(huì)通過抑制鹽度垂向混合從而對同期鹽度負(fù)距平的形成有促進(jìn)作用。為進(jìn)一步了解這兩年春季鹽度垂直混合過程的差異,這里挑選了這兩年春季相對應(yīng)的Argo浮標(biāo)鹽度剖面進(jìn)行比較,將距離在1°以內(nèi),相對時(shí)間(基準(zhǔn)時(shí)間為3月1日)相差在±5 d內(nèi)為一組配對的剖面,共計(jì)21組鹽度剖面入選(圖8c和圖8d)??梢钥闯鲈?50 m層以上,2011年春季鹽度垂直混合較好,2012年鹽度層化明顯,反映了這兩年風(fēng)應(yīng)力旋度的差異對上層鹽度垂直混合過程的影響。

    圖8 2011年(a)和2012年(b)孟加拉灣春季(3-4月)海面風(fēng)場(矢量,單位:m/s)及風(fēng)應(yīng)力旋度(填色圖,向上為正,單位:N/m3)距平。c和d分別同a和b,但分別為Argo浮標(biāo)鹽度剖面Fig.8 The distribution of spring(M arch to April)surface wind anomalies(vectors,unit:m/s)and wind stress curl anomalies(shaded,upward is positive,unit:N/m3)in the B O B during 2011(a)and 2012(b).c and d are the same period as a and b,respectively,but for salinity profiles from Argo

    7 結(jié)論

    通過Argo浮標(biāo)鹽度資料對2011年8月至2014 年3月孟加拉灣Aquarius 2.0版L3級月平均遙感鹽度數(shù)據(jù)進(jìn)行對比驗(yàn)證,其結(jié)果表明:在本研究海域A-quarius海表鹽度與Argo浮標(biāo)表層鹽度呈顯著性正相關(guān)(相關(guān)系數(shù)為0.93),均方根差異為0.48,并比Argo鹽度偏低(總體偏差為-0.13),其中孟加拉灣北部低鹽區(qū)比南部高鹽區(qū)偏低更多(分別為-0.28 和-0.10)。其偏低主因在于遙感測量的是海表鹽度,而大部分Argo剖面最淺只觀測到5 m左右深度,于是二者觀測層次上的差異構(gòu)成了此系統(tǒng)偏差的一個(gè)來源。在灣北部受降水及陸地徑流注入的影響,近表層的鹽度層化較灣南部更強(qiáng),偏低更顯著。此外,灣北部近岸的強(qiáng)射頻干擾可能也有一定影響。

    分析結(jié)果還表明,Aquarius數(shù)據(jù)清晰呈現(xiàn)了孟加拉灣表層鹽度的季節(jié)變化特征,與已有的認(rèn)識基本一致,包括阿拉伯海高鹽水的入侵引起孟加拉灣南部海域鹽度的變化以及北部淡水羽分布范圍的季節(jié)性遷移等主要特征。Aquarius還揭示了2012年春季整個(gè)灣內(nèi)出現(xiàn)顯著異常低鹽現(xiàn)象,得到了Argo及R A M A等多種實(shí)測資料的證實(shí),同時(shí)Argo鹽度資料表明2011年春季灣內(nèi)存在異常高鹽的現(xiàn)象。分析表明,這兩年春季鹽度變化分別與各自前一個(gè)夏季降雨異常變化引起灣北部表層鹽度出現(xiàn)相應(yīng)的變化有關(guān),北部鹽度變化信號通過反氣旋環(huán)流向南輸運(yùn)最終引起春季出現(xiàn)海盆尺度的鹽度異常變化。此外,2011 (2012)年春季風(fēng)應(yīng)力旋度正(負(fù))異常通過影響鹽度垂直混合過程對該異常高(低)鹽現(xiàn)象也有貢獻(xiàn)。

    致謝:感謝國家海洋局第三海洋研究所郭小鋼研究員給出的寶貴意見。

    參考文獻(xiàn):

    [1]Sengupta D,Bharath Raj G N,Shenoi S S C.Surfacefreshwaterfrom Bay of Bengalrunoff and Indonesian throughflow in the tropicalIndian Ocean [J].Geophysical Research Letters,2006,33(22):L22609.

    [2]Dube S K,Rao A D,Sinha P C,et al.Implications of climatic variationsin the fresh water outflow on the wind-induced circulation of the Bay of Bengal[J].Atmospheric Environ ment,1995,29(16):2133-2138.

    [3]Vinayachandran P N,M urty V S N,Ramesh Babu V.Observations of barrier layer formation in the Bay of Bengal during su m mer monsoon[J].Journal of Geophysical Research:Oceans(1978-2012),2002,107(C12):SR F 19-1-SR F 19-9.

    [4]Thadathil P,M uraleedharan P M,Rao R R,et al.Observed seasonal variability of barrier layer in the Bay of Bengal[J].Journal of Geophysical Research:Oceans(1978-2012),2007,112(C2):C02009.

    [5]M asson S,Luo J J,M adec G,et al.Impact of barrierlayer on winter-spring variability of the southeastern Arabian Sea[J].Geophysical research letters,2005,32(7):L07703.

    [6]Varkey M J,M urty V S N,Suryanarayana A.Physical oceanography ofthe Bay of Bengal and Andaman Sea[J].Oceanography and M arine Biology:An Annual Review,1996,34:1-70.

    [7]張玉紅,杜巖,徐海明,等.赤道印度洋中部斷面東西水交換的季節(jié)變化及其區(qū)域差異[J].海洋學(xué)報(bào),2012,34(2):30-38.Zhang Yuhong,Du Yan,Xu Haiming,et al.Seasonal variability of the west-east water mass exchange on the section of central equatorial Indian Ocean and its regional difference[J].Haiyang Xuebao,2012,34(2):30-38.

    [8]Akhil V P,Durand F,Lengaigne M,et al.A modeling study ofthe processes of surface salinity seasonal cyclein the Bay of Bengal[J].Journal of Geophysical Research:Oceans,2014,119(6):3926-3947.

    [9]M urty V S N,Sarma Y V B,Rao D P,et al.W ater characteristics,mixing and circulation in the Bay of Bengal during southwest monsoon[J].Journal of M arine Research,1992,50(2):207-228.

    [10]W u L J,W ang F,Yuan D L,et al.Evolution offreshwater plu mes and salinity frontsin the northern Bay of Bengal[J].Journal of Geophysical Research:Oceans(1978—2012),2007,112(C08):C08017.

    [11]Vinayachandran P N,Shankar D,Vernekar S,et al.A su m mer monsoon pu mp to keep the Bay of Bengal salty[J].Geophysical Research Letters,2013,40(9):1777-1782.

    [12]Zhang Y H,Du Y,Zheng S J,et al.Impact of Indian Ocean Dipole on the salinity budgetin the equatorialIndian Ocean[J].Journal of Geophysical Research:Oceans,2013,118(10):4911-4923.

    [13]Subrah manyam B,M urty V S N,Heffner D M.Sea surface salinity variability in the tropical Indian Ocean[J].Remote Sensing of Environ ment,2011,115(3):944-956.

    [14]Zeng L L,Timothy Liu W,Xue H J,et al.Freshening in the South China Sea during 2012 revealed by Aquarius and in situ data[J].Journal of Geophysical Research:Oceans,2014,119(12):8296-8314.

    [15]Du Y,Zhang Y H.Satellite and Argo observed surface salinity variationsin the tropical Indian Ocean and their association with the Indian Ocean Dipole mode[J].Journal of Climate,2015,28(2):695-713.

    [16]Le Vine D M,Lagerloef G S E,Colomb F R,et al.Aquarius:An instru ment to monitor sea surface salinity from space[J].Geoscience and Remote Sensing,IEEE Transactions on,2007,45(7):2040-2050.

    [17]Roem mich D,Gilson J.The 2004—2008 mean and annualcycle oftemperature,salinity,and steric heightin the global ocean from the Argo Program[J].Progress in Oceanography,2009,82(2):81-100.

    [18]H uffman G J,Bolvin D T,Nelkin E J,et al.The T R M M multisatellite precipitation analysis(T M P A):Quasi-global,multiyear,combined-sensor precipitation estimates at fine scales[J].Journal of H ydrometeorology,2007,8(1):38-55.

    [19]Yu L S,W eller R A.Objectively analyzed air-sea heat fluxes for the globalice-free oceans(1981-2005)[J].Bulletin of the A merican M eteoro-logical Society,2007,88(4):527-539.

    [20]Oki T,Sud Y C.Design of Total Runoff Integrating Pathways(T RIP)-A global river channel network[J].Earth interactions,1998,2(1):1-37.

    [21]Tang W Q,Yueh S H,F(xiàn)ore A G,et al.Validation of Aquarius sea surface salinity with in situ measurementsfrom Argo floats and moored buoys [J].Journal of Geophysical Research:Oceans,2014,119(9):6171-6189.

    [22]Reagan J,Boyer T,Antonov J,et al.Comparison analysis between Aquarius sea surface salinity and W orld Ocean Database in situ analyzed sea surface salinity[J].Journal of Geophysical Research:Oceans,2014,119(11):8122-8140.

    [23]Le Vine D M,Lagerloef G S E,Pellerano F,et al.The Aquarius/SA C-D mission and status of the Aquariusinstru ment[C]//Microwave Radiometry and Remote Sensing of the Environ ment,2008.MIC R O R A D 2008.IEEE,2008:1-4.

    [24]Ebuchi N,Abe H.Evaluation of sea surface salinity observed by Aquarius[C]//Geoscience and Remote Sensing Symposiu m(IG A RSS),2014 IEEE International.IEEE,2014:4427-4430.

    [25]王新新,楊建洪,趙冬至,等.Aquarius/SA C-D衛(wèi)星鹽度數(shù)據(jù)在中國南海的準(zhǔn)確度評估[J].熱帶海洋學(xué)報(bào),2013,32(5):23-28.W ang Xinxin,Yang Jianghong,Zhao Dongzhi,et al.Assessment of Aquarius/SA C-D salinity data accuracyin the South China Sea[J].Journal of Tropical Oceanography,2013,32(5):23-28.

    [26]Xuan L L,Qiu Y,Xu J D,et al.H ydrography and circulation in the eastern tropicalIndian Ocean during April-M ay 2011[J].Atmospheric and Oceanic Science Letters,2012,5(4):284-289.

    [27]Dai A,Trenberth K E.Estimates offreshwater discharge from continents:Latitudinal and seasonal variations[J].Journal of hydrometeorology,2002,3(6):660-687.

    [28]Pant V,Girishku mar M S,Udaya Bhaskar T V S,et al.Observed interannual variability of near-surface salinity in the Bay of Bengal[J].Journal of Geophysical Research:Oceans,2015,120(5):3315-3329.

    [29]邱云,李立.孟加拉灣上層地轉(zhuǎn)環(huán)流周年變化的遙感研究[J].海洋學(xué)報(bào),2007,29(3):39-46.Qiu Yun,Li Li.Annual variation of geostrophic circulation in the Bay of Bengal observed with T O PE X/Poseidon altimeter data[J].Haiyang Xuebao,2007,29(3):39-46.

    林新宇,邱云,杜巖,等.孟加拉灣海表鹽度變化特征的衛(wèi)星遙感應(yīng)用研究[J].海洋學(xué)報(bào),2016,38(5):46-56,doi:10.3969/ j.issn.0253-4193.2016.05.005

    Lin Xinyu,Qiu Yun,Du Yan,et al.Study on the utilization of satellite remote sensing for variation characteristic of sea surface salinity in the Bay of Bengal[J].Haiyang Xuebao,2016,38(5):46-56,doi:10.3969/j.issn.0253-4193.2016.05.005

    Study on the utilization of satellite remote sensing for variation characteristic of sea surface salinity in the Bay of Bengal

    Lin Xinyu1,Qiu Yun1,Du Yan2,Zhou Xiwu1,Xuan Lili3,Xu Jindian1
    (1.Third Institute of Oceanography,State Oceanic Administration,Xiamen 361005,China;2.State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology,Chinese Academy of Sciences,Guangzhou 510301,China;3.National Marine Data and Information Service,Tianjin 300171,China)

    Abstract:Sea surface salinity(SSS)from Aquarius during August 2011 and M arch 2014,together with insitu salinity data(e.g.Argo)are used to explore seasonal and interannual variations of the SSS in the Bay of Bengal.It shows that Aquarius SSS has a significant positive correlation with Argo SSS.The Aquarius SSS is relatively lower than Argo SSS,and there exists a negative bias of-0.13 between these two data.The negative bias is larger in the northern bay than thatin the southern bay(-0.28 vs-0.10).The bias mainly comes from the different depth between Aquarius and Argo floatsin measuring SSS.Aquarius SSS captures main seasonal variations of the SSS in the bay,such as the intrusion of saltier water from Arabian Sea and seasonal migrations of the freshwater plu mes in the northern bay.Besides,the SSS reveals anomalous salting(freshening)in the bay during the spring 2011(2012).Our results suggestthatthe reduced(enhanced)precipitation during the su m mer 2010(2011)in the northern bay reduced(enhanced)freshening in that region,and the anomalous salting(freshening)water was advected southward by the surface-layer circulation causing anomalous salting(freshening)in the whole bay.M eanwhile,the positive(negative)wind stress curl anomaly also contribute to the anomalous salting(freshening)event during spring 2011(2012)in the bay via impacting the vertical mixing processes.

    Key words:Bay of Bengal;Aquarius;sea surface salinity;ocean circulation

    *通信作者:邱云(1979—),博士,研究員,主要研究熱帶邊緣海及熱帶大洋海洋環(huán)流動(dòng)力學(xué)與海氣相互作用。E-mail:qiuyun@tio.org.cn

    作者簡介:林新宇(1990—),女,福建省福鼎市人,主要從事熱帶邊緣海海洋環(huán)流動(dòng)力學(xué)研究。E-mail:linxinyu@tio.org.cn

    基金項(xiàng)目:海洋三所基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(海三科2009046);國家自然科學(xué)基金項(xiàng)目(41276034,40806014,41306027);熱帶海洋環(huán)境國家重點(diǎn)實(shí)驗(yàn)室(中國科學(xué)院南海海洋研究所)開放課題(L T O1402);海洋公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目(201005033-4)。

    收稿日期:2015-06-29;

    修訂日期:2015-10-26。

    中圖分類號:P731.12;P716+.14

    文獻(xiàn)標(biāo)志碼:A

    文章編號:0253-4193(2016)05-0046-11

    猜你喜歡
    孟加拉灣環(huán)流
    內(nèi)環(huán)流控溫技術(shù)應(yīng)用實(shí)踐與發(fā)展前景
    溫暖的墨西哥灣
    孟加拉灣東海岸波浪特征分析
    2017年8月9日~11日林芝暴雨過程分析
    孟加拉灣春季小型暖池對熱帶氣旋的影響研究
    一種風(fēng)電變流器并聯(lián)環(huán)流抑制方法
    熱鹽環(huán)流方程全局弱解的存在性
    謎底大揭秘
    孟加拉灣北部逆溫現(xiàn)象的觀測特征及機(jī)制分析
    兩相坐標(biāo)系下MMC環(huán)流抑制策略
    尤物成人国产欧美一区二区三区| 久久精品国产亚洲av涩爱 | 观看免费一级毛片| a在线观看视频网站| 丰满乱子伦码专区| 久久精品国产自在天天线| 在线观看av片永久免费下载| 国产视频一区二区在线看| 欧美日韩综合久久久久久 | 日本黄大片高清| 欧美日韩综合久久久久久 | 欧美zozozo另类| or卡值多少钱| 免费av毛片视频| 国产精品久久久久久精品电影| 欧美黑人欧美精品刺激| 国产高清激情床上av| www日本黄色视频网| 国产免费一级a男人的天堂| 日韩中文字幕欧美一区二区| 日本黄色视频三级网站网址| 国产熟女xx| 国产精品电影一区二区三区| 国产探花在线观看一区二区| 亚洲av成人精品一区久久| 国产在线精品亚洲第一网站| 免费一级毛片在线播放高清视频| 免费观看精品视频网站| 天堂av国产一区二区熟女人妻| 可以在线观看毛片的网站| 波多野结衣巨乳人妻| 中文资源天堂在线| 欧美最黄视频在线播放免费| 看黄色毛片网站| 香蕉久久夜色| 免费av不卡在线播放| 97人妻精品一区二区三区麻豆| 久久久国产精品麻豆| 最近在线观看免费完整版| 国产高清videossex| 少妇高潮的动态图| 九九久久精品国产亚洲av麻豆| 在线观看美女被高潮喷水网站 | 高清毛片免费观看视频网站| 亚洲五月婷婷丁香| 麻豆久久精品国产亚洲av| 国产高清videossex| 日本三级黄在线观看| 日韩精品中文字幕看吧| 中亚洲国语对白在线视频| 性色av乱码一区二区三区2| 岛国在线免费视频观看| 天天躁日日操中文字幕| 亚洲精华国产精华精| 天天躁日日操中文字幕| 老鸭窝网址在线观看| 一边摸一边抽搐一进一小说| 日本一二三区视频观看| 美女cb高潮喷水在线观看| 久久精品91无色码中文字幕| 欧美在线一区亚洲| 深爱激情五月婷婷| 88av欧美| 婷婷精品国产亚洲av在线| 国产黄片美女视频| 国产激情欧美一区二区| 婷婷丁香在线五月| 一边摸一边抽搐一进一小说| 一本久久中文字幕| 国产精品嫩草影院av在线观看 | 久久九九热精品免费| 熟女人妻精品中文字幕| 色视频www国产| 精品日产1卡2卡| 天堂影院成人在线观看| 免费观看精品视频网站| 国产精品 欧美亚洲| 国产在视频线在精品| 中文在线观看免费www的网站| 最近最新中文字幕大全电影3| 亚洲在线自拍视频| 可以在线观看毛片的网站| 国产免费男女视频| 亚洲av二区三区四区| 精品国产亚洲在线| 老熟妇仑乱视频hdxx| 欧美最新免费一区二区三区 | 国产精品久久久久久亚洲av鲁大| 91字幕亚洲| 九九久久精品国产亚洲av麻豆| 色哟哟哟哟哟哟| 久99久视频精品免费| 99久久综合精品五月天人人| 国产av一区在线观看免费| 亚洲人成网站高清观看| 国产色爽女视频免费观看| 日日夜夜操网爽| 亚洲不卡免费看| 国产一区二区激情短视频| 国产国拍精品亚洲av在线观看 | 欧美高清成人免费视频www| 精品无人区乱码1区二区| 亚洲精品粉嫩美女一区| 悠悠久久av| 老司机在亚洲福利影院| 真实男女啪啪啪动态图| 最近最新免费中文字幕在线| www日本黄色视频网| 老司机在亚洲福利影院| 啦啦啦观看免费观看视频高清| 久久中文看片网| 草草在线视频免费看| 身体一侧抽搐| 嫩草影视91久久| 国产主播在线观看一区二区| 有码 亚洲区| а√天堂www在线а√下载| 最近在线观看免费完整版| 国产v大片淫在线免费观看| 国产精品一区二区免费欧美| 久久精品国产综合久久久| 久久久久国产精品人妻aⅴ院| 久久香蕉精品热| 91字幕亚洲| 婷婷六月久久综合丁香| 久久久久久国产a免费观看| 午夜亚洲福利在线播放| 男插女下体视频免费在线播放| bbb黄色大片| 亚洲av不卡在线观看| 国产又黄又爽又无遮挡在线| 国产激情偷乱视频一区二区| 成人高潮视频无遮挡免费网站| 色哟哟哟哟哟哟| 中文字幕av成人在线电影| 色噜噜av男人的天堂激情| 国产三级黄色录像| 欧美黑人巨大hd| 亚洲欧美激情综合另类| 日本在线视频免费播放| 国产伦精品一区二区三区四那| 欧美区成人在线视频| 精品人妻1区二区| 国产精品,欧美在线| 日本免费a在线| 午夜精品久久久久久毛片777| 熟女少妇亚洲综合色aaa.| 日韩欧美在线乱码| 午夜两性在线视频| 精华霜和精华液先用哪个| 国产精品一区二区三区四区久久| 男女之事视频高清在线观看| 亚洲最大成人手机在线| 非洲黑人性xxxx精品又粗又长| 中文字幕人妻熟人妻熟丝袜美 | 1000部很黄的大片| 国产高清激情床上av| 免费大片18禁| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 真实男女啪啪啪动态图| 亚洲一区二区三区不卡视频| 日本黄大片高清| 男人舔女人下体高潮全视频| 亚洲熟妇熟女久久| 日日干狠狠操夜夜爽| 俺也久久电影网| 免费看日本二区| 精品久久久久久久人妻蜜臀av| 国产精品亚洲一级av第二区| 亚洲精品日韩av片在线观看 | 一级作爱视频免费观看| 法律面前人人平等表现在哪些方面| 天堂av国产一区二区熟女人妻| 中文亚洲av片在线观看爽| 我的老师免费观看完整版| 国产一区二区在线观看日韩 | 美女cb高潮喷水在线观看| 亚洲av电影不卡..在线观看| 国产熟女xx| 人妻夜夜爽99麻豆av| 3wmmmm亚洲av在线观看| 国产精品三级大全| 久久精品91无色码中文字幕| 欧美3d第一页| 国产aⅴ精品一区二区三区波| 岛国在线观看网站| 一级毛片女人18水好多| 丁香六月欧美| 午夜两性在线视频| 一a级毛片在线观看| 听说在线观看完整版免费高清| 毛片女人毛片| 欧洲精品卡2卡3卡4卡5卡区| 欧美绝顶高潮抽搐喷水| 一区二区三区高清视频在线| 岛国视频午夜一区免费看| 国产精品久久久久久人妻精品电影| 老司机午夜十八禁免费视频| 中文字幕久久专区| 99视频精品全部免费 在线| 亚洲片人在线观看| 亚洲av一区综合| 99热6这里只有精品| 免费一级毛片在线播放高清视频| 一个人看视频在线观看www免费 | 岛国在线免费视频观看| 一级黄片播放器| 中文字幕久久专区| 波野结衣二区三区在线 | 琪琪午夜伦伦电影理论片6080| 国产成人影院久久av| 欧美日韩亚洲国产一区二区在线观看| 香蕉久久夜色| 亚洲精品亚洲一区二区| 亚洲av美国av| 两个人视频免费观看高清| 日本撒尿小便嘘嘘汇集6| 丰满人妻熟妇乱又伦精品不卡| 韩国av一区二区三区四区| 人妻夜夜爽99麻豆av| 他把我摸到了高潮在线观看| 亚洲av成人av| 久久久精品欧美日韩精品| av视频在线观看入口| 免费观看的影片在线观看| 亚洲国产精品久久男人天堂| 国产亚洲精品一区二区www| ponron亚洲| 日本精品一区二区三区蜜桃| 亚洲精华国产精华精| 国产成人欧美在线观看| 午夜日韩欧美国产| 国产爱豆传媒在线观看| 宅男免费午夜| 欧美成人a在线观看| 亚洲久久久久久中文字幕| 又黄又粗又硬又大视频| 男女午夜视频在线观看| 九九热线精品视视频播放| 真人一进一出gif抽搐免费| 观看免费一级毛片| 手机成人av网站| 国产真人三级小视频在线观看| 免费大片18禁| 又爽又黄无遮挡网站| 在线播放无遮挡| 久久欧美精品欧美久久欧美| 露出奶头的视频| 18禁黄网站禁片午夜丰满| 久久久精品大字幕| 日本成人三级电影网站| 国产蜜桃级精品一区二区三区| 日本 欧美在线| 美女高潮的动态| 国产伦在线观看视频一区| 免费看a级黄色片| 国产一区二区激情短视频| 欧美日韩一级在线毛片| 亚洲欧美日韩高清在线视频| 成年女人看的毛片在线观看| 国产私拍福利视频在线观看| 久久久久久久精品吃奶| 女人被狂操c到高潮| 在线播放无遮挡| 日韩精品青青久久久久久| 亚洲av日韩精品久久久久久密| 一级作爱视频免费观看| 夜夜躁狠狠躁天天躁| 身体一侧抽搐| 99热只有精品国产| 国产精华一区二区三区| 久久久久久久亚洲中文字幕 | 国产精品一区二区三区四区免费观看 | 性欧美人与动物交配| 欧美日本视频| 男人舔奶头视频| 国产精品女同一区二区软件 | 国产精品爽爽va在线观看网站| 女人被狂操c到高潮| 又黄又粗又硬又大视频| 国内揄拍国产精品人妻在线| 国产欧美日韩一区二区精品| 国产亚洲精品一区二区www| 日本 av在线| 国产精品女同一区二区软件 | 亚洲熟妇中文字幕五十中出| 给我免费播放毛片高清在线观看| 一本精品99久久精品77| 国产精品久久久久久久电影 | 午夜老司机福利剧场| 丰满人妻一区二区三区视频av | 乱人视频在线观看| 18美女黄网站色大片免费观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲欧美激情综合另类| 97超级碰碰碰精品色视频在线观看| 亚洲av成人精品一区久久| 97超级碰碰碰精品色视频在线观看| 看免费av毛片| 精品一区二区三区人妻视频| 男女床上黄色一级片免费看| 精品欧美国产一区二区三| 成人欧美大片| 99久久久亚洲精品蜜臀av| 一二三四社区在线视频社区8| 午夜激情欧美在线| 香蕉久久夜色| avwww免费| 久久久久九九精品影院| 一级毛片女人18水好多| 日韩免费av在线播放| 国产一区二区亚洲精品在线观看| 亚洲成av人片免费观看| 一级毛片女人18水好多| 欧美一区二区精品小视频在线| 亚洲国产欧洲综合997久久,| 成年人黄色毛片网站| 免费观看精品视频网站| 亚洲中文字幕日韩| 欧美性猛交黑人性爽| 狂野欧美激情性xxxx| 欧美日韩精品网址| 国产免费av片在线观看野外av| 女人十人毛片免费观看3o分钟| 国产av一区在线观看免费| 国产亚洲精品久久久久久毛片| 国产国拍精品亚洲av在线观看 | 蜜桃亚洲精品一区二区三区| 欧美+日韩+精品| 午夜福利在线在线| 日韩欧美免费精品| 一卡2卡三卡四卡精品乱码亚洲| 日本黄色视频三级网站网址| 日本精品一区二区三区蜜桃| 国产麻豆成人av免费视频| 99久久99久久久精品蜜桃| 3wmmmm亚洲av在线观看| 嫁个100分男人电影在线观看| 岛国在线免费视频观看| 日韩av在线大香蕉| 久久久久精品国产欧美久久久| 色老头精品视频在线观看| 国产一区二区在线观看日韩 | 久久久久九九精品影院| 国产精品电影一区二区三区| 亚洲成人免费电影在线观看| 亚洲电影在线观看av| 精品一区二区三区人妻视频| or卡值多少钱| 欧美大码av| 午夜福利18| 成人无遮挡网站| 亚洲国产欧美网| 波多野结衣高清无吗| 国产精品久久久久久人妻精品电影| 夜夜躁狠狠躁天天躁| 神马国产精品三级电影在线观看| 一个人观看的视频www高清免费观看| 精品人妻一区二区三区麻豆 | 国产精品久久久久久亚洲av鲁大| 亚洲内射少妇av| 操出白浆在线播放| 91九色精品人成在线观看| 色综合亚洲欧美另类图片| 可以在线观看的亚洲视频| 在线十欧美十亚洲十日本专区| 久久久久久久久大av| 精品无人区乱码1区二区| 欧美不卡视频在线免费观看| 久久久成人免费电影| 亚洲最大成人手机在线| 99久久99久久久精品蜜桃| 少妇人妻一区二区三区视频| 日日干狠狠操夜夜爽| 国产精品永久免费网站| 此物有八面人人有两片| 久久亚洲真实| 午夜久久久久精精品| 亚洲最大成人中文| 精品国产美女av久久久久小说| 亚洲人成电影免费在线| 国产精品乱码一区二三区的特点| 免费av毛片视频| 三级毛片av免费| 国产精品98久久久久久宅男小说| 亚洲,欧美精品.| 成人鲁丝片一二三区免费| 国产老妇女一区| 深爱激情五月婷婷| 中出人妻视频一区二区| 亚洲成人免费电影在线观看| 19禁男女啪啪无遮挡网站| 在线看三级毛片| 小说图片视频综合网站| 国产黄a三级三级三级人| xxxwww97欧美| 国产黄片美女视频| 床上黄色一级片| 狂野欧美激情性xxxx| 欧美日韩黄片免| 香蕉丝袜av| 免费在线观看日本一区| 亚洲精品日韩av片在线观看 | 亚洲男人的天堂狠狠| 免费在线观看影片大全网站| 18禁美女被吸乳视频| 亚洲av成人精品一区久久| 18禁黄网站禁片免费观看直播| 免费观看的影片在线观看| 欧美日韩国产亚洲二区| 国产欧美日韩一区二区三| 亚洲欧美激情综合另类| 小说图片视频综合网站| 欧美日韩中文字幕国产精品一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲成人久久性| 日本a在线网址| 精品一区二区三区av网在线观看| 国产男靠女视频免费网站| 日韩高清综合在线| 99久久99久久久精品蜜桃| 国产91精品成人一区二区三区| 在线看三级毛片| 成年女人看的毛片在线观看| 超碰av人人做人人爽久久 | 国产真人三级小视频在线观看| 一夜夜www| 啦啦啦韩国在线观看视频| 日韩 欧美 亚洲 中文字幕| 啪啪无遮挡十八禁网站| 国产精品久久久久久亚洲av鲁大| 亚洲色图av天堂| www日本黄色视频网| 国产真实乱freesex| 少妇熟女aⅴ在线视频| 国产精品综合久久久久久久免费| 国内精品美女久久久久久| 国产精品自产拍在线观看55亚洲| 国产美女午夜福利| 在线天堂最新版资源| 最近视频中文字幕2019在线8| 日本一二三区视频观看| 国内精品美女久久久久久| 中文字幕熟女人妻在线| 国产成人啪精品午夜网站| 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 国产成人a区在线观看| 午夜免费观看网址| 日本三级黄在线观看| 少妇丰满av| 91麻豆av在线| 2021天堂中文幕一二区在线观| 99热6这里只有精品| 国产免费一级a男人的天堂| 美女cb高潮喷水在线观看| 亚洲 欧美 日韩 在线 免费| 十八禁人妻一区二区| 国产精品一区二区三区四区免费观看 | av中文乱码字幕在线| 别揉我奶头~嗯~啊~动态视频| 人妻夜夜爽99麻豆av| 久久精品91无色码中文字幕| 麻豆一二三区av精品| 国产精品一及| 成熟少妇高潮喷水视频| 日本五十路高清| 成人av在线播放网站| 精品国产亚洲在线| 亚洲精品影视一区二区三区av| 精品久久久久久久人妻蜜臀av| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 欧美日韩精品网址| 久久久久久大精品| 国产精品永久免费网站| av中文乱码字幕在线| 国产亚洲av嫩草精品影院| 99国产精品一区二区三区| 天堂√8在线中文| 欧美另类亚洲清纯唯美| 国产免费一级a男人的天堂| 欧美zozozo另类| 人人妻人人看人人澡| 宅男免费午夜| 中出人妻视频一区二区| 午夜福利高清视频| 天天躁日日操中文字幕| 好看av亚洲va欧美ⅴa在| 十八禁网站免费在线| 亚洲久久久久久中文字幕| av中文乱码字幕在线| 日韩亚洲欧美综合| 激情在线观看视频在线高清| 国产一区二区三区在线臀色熟女| 国产 一区 欧美 日韩| 日韩成人在线观看一区二区三区| 国产淫片久久久久久久久 | 国产精品永久免费网站| 亚洲片人在线观看| 18禁裸乳无遮挡免费网站照片| 日韩有码中文字幕| 午夜两性在线视频| 啪啪无遮挡十八禁网站| 99国产精品一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲av中文字字幕乱码综合| 亚洲18禁久久av| 两个人看的免费小视频| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 精品国产三级普通话版| av黄色大香蕉| 天堂动漫精品| 亚洲欧美日韩高清在线视频| 99在线人妻在线中文字幕| 日日夜夜操网爽| 亚洲内射少妇av| 亚洲国产欧美网| 欧美一区二区亚洲| 内射极品少妇av片p| 女人十人毛片免费观看3o分钟| 色老头精品视频在线观看| 亚洲精品影视一区二区三区av| 在线天堂最新版资源| 亚洲欧美日韩高清专用| 久久精品人妻少妇| 很黄的视频免费| 亚洲成人精品中文字幕电影| e午夜精品久久久久久久| 亚洲精品日韩av片在线观看 | 国产成人欧美在线观看| 色噜噜av男人的天堂激情| 国内久久婷婷六月综合欲色啪| 国产主播在线观看一区二区| 在线播放国产精品三级| 久久精品国产亚洲av香蕉五月| 国产高清videossex| 黄色女人牲交| www.熟女人妻精品国产| 国产私拍福利视频在线观看| 人妻丰满熟妇av一区二区三区| 精品乱码久久久久久99久播| 人妻丰满熟妇av一区二区三区| 久久久久久国产a免费观看| 中文字幕高清在线视频| 99国产综合亚洲精品| 国产免费男女视频| 制服丝袜大香蕉在线| 国产精品99久久久久久久久| 欧美zozozo另类| 波野结衣二区三区在线 | 久久国产精品人妻蜜桃| 日本精品一区二区三区蜜桃| 久久久精品大字幕| 成年女人永久免费观看视频| 中文字幕av在线有码专区| 叶爱在线成人免费视频播放| 国产探花极品一区二区| 久久6这里有精品| 又紧又爽又黄一区二区| 美女cb高潮喷水在线观看| 我的老师免费观看完整版| 少妇人妻精品综合一区二区 | 夜夜看夜夜爽夜夜摸| 日本 欧美在线| 国产亚洲精品久久久久久毛片| 久久久久九九精品影院| 国产极品精品免费视频能看的| 欧美乱妇无乱码| 不卡一级毛片| 亚洲人成电影免费在线| 波多野结衣高清无吗| 久久久国产成人精品二区| 3wmmmm亚洲av在线观看| 亚洲avbb在线观看| 一区福利在线观看| 国产又黄又爽又无遮挡在线| 色吧在线观看| 两个人看的免费小视频| 免费高清视频大片| 窝窝影院91人妻| 中文字幕熟女人妻在线| 亚洲一区二区三区不卡视频| 岛国在线免费视频观看| 久久久久免费精品人妻一区二区| 亚洲av第一区精品v没综合| 亚洲国产精品合色在线| 嫩草影院入口| 亚洲精品在线美女| 亚洲人与动物交配视频| 男女午夜视频在线观看| 亚洲av电影不卡..在线观看| 日韩欧美一区二区三区在线观看| 男女那种视频在线观看| 亚洲精品粉嫩美女一区| 午夜日韩欧美国产| 亚洲七黄色美女视频| 久久国产精品影院| 午夜精品在线福利| 亚洲专区中文字幕在线| 欧美中文日本在线观看视频| 搡女人真爽免费视频火全软件 | 精品99又大又爽又粗少妇毛片 | 亚洲天堂国产精品一区在线| 男女床上黄色一级片免费看| 国产激情偷乱视频一区二区| 18禁在线播放成人免费| 欧美不卡视频在线免费观看| 一本久久中文字幕| 精品人妻偷拍中文字幕| 久久久国产成人精品二区| 热99在线观看视频| 尤物成人国产欧美一区二区三区|