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

    大尺度環(huán)流異常對南極印度洋扇區(qū)海表面高鹽異常的影響

    2021-08-21 14:06:38倪旭彬杜凌史荒原
    海洋學(xué)報(bào) 2021年7期
    關(guān)鍵詞:海冰

    倪旭彬,杜凌*,史荒原

    (1.中國海洋大學(xué) 海洋與大氣學(xué)院,山東 青島 266100)

    1 引言

    氣候變化背景下,南極表層鹽度的變化特征不僅體現(xiàn)了全球水循環(huán)的改變,也與海洋和大氣中的大尺度環(huán)流緊密聯(lián)系,對全球水循環(huán)和未來氣候變化都將產(chǎn)生深遠(yuǎn)影響[1]。與南極海冰變化類似,南極表層鹽度也表現(xiàn)出復(fù)雜的變化特征,如不同鹽度數(shù)據(jù)集呈現(xiàn)出的20 世紀(jì)后半葉的趨勢變化,在60°S 以南海域的差別很顯著[2]。南極印度洋扇區(qū)是南大洋變化的重要海區(qū)[3]。近年來,南極印度洋扇區(qū)的高鹽異常現(xiàn)象十分明顯,咸化趨勢在60°S 以南、500 m 以淺的上層海洋體現(xiàn)得最為顯著[4]。該海域包含我國南大洋調(diào)查的重點(diǎn)研究海域——普里茲灣,眾多學(xué)者對該海域水團(tuán)也開展了廣泛深入的研究,發(fā)現(xiàn)存在3 種高鹽水體,均為該海域的永久性水團(tuán),分別為67°S 以南至陸架區(qū)域的普里茲灣高鹽陸架水、埃默里冰架前緣的高鹽冰架水,以及位于63°~64°S 的繞極深層水[5-7]。徐智昕等[7]指出2011 年普里茲灣灣口與陸架區(qū)分布著一條顯著高鹽帶,鹽度超過34.4。林麗娜等[8]的研究結(jié)果也表明,2009 年普里茲灣陸架的大部分水體出現(xiàn)高鹽異常,比2004 年高約0.75,且高鹽陸架水向北擴(kuò)張。最近,Masuda[9]使用EN4 數(shù)據(jù)探究南極印度洋扇區(qū)上層水體的鹽度變化特征,發(fā)現(xiàn)表層與次表層水體因受海冰等過程影響,鹽度季節(jié)變化較為明顯,其中100 m 以淺的表層水體鹽度在冬季最大,夏季最小。

    南極上層海洋鹽度變化與大尺度環(huán)流異常聯(lián)系密切。已有的研究結(jié)果表明,南極濤動(dòng)(Antarctic Oscillation,AAO)與印度洋偶極子(Indian Ocean Dipole,IOD)對南極海表面鹽度變化的影響不可忽視。AAO正位相期間,南半球中緯度西風(fēng)增強(qiáng),造成向北的??寺斔停约按罅康暮1蒣10],為上層海洋貢獻(xiàn)更多的鹽通量[11],促進(jìn)了表層與次表層水體鹽度的增大[9]。另一方面,AAO 正位相下南極沿岸東風(fēng)減弱,加深了沿岸的密度躍層,促使高鹽的繞極深層水入侵陸架,有利于維持上層海洋的高鹽特征[12]。由于AAO與IOD 強(qiáng)正位相的共同影響,南極印度洋扇區(qū)的海冰生成將顯著增多,從而為海表面提供更多鹽通量[13]。同時(shí),IOD 的正位相能夠增強(qiáng)南極印度洋扇區(qū)沿岸的海表面蒸發(fā)[14],同樣有利于表層海水出現(xiàn)高鹽異常。

    局地過程對海表面鹽度也有重要影響,有貢獻(xiàn)的局地因素主要包括蒸發(fā)增加與降水減少[14]、淡水平流減弱[15],以及南極沿岸海冰生成的增強(qiáng)[16]。局地風(fēng)和風(fēng)場旋度變化影響蒸發(fā)與降水的強(qiáng)弱,是影響海表面鹽度的淡水通量的重要原因。南極印度洋扇區(qū)及其以北的海域盛行下降風(fēng)且風(fēng)速較大,使得蒸發(fā)增強(qiáng)、降水減弱[17],也促使海表面鹽度增加。局地風(fēng)也可以通過??寺槲龠M(jìn)混合層加深,有利于下層的高鹽水與表層水之間的對流混合[18],增大上層水體的鹽度并減弱層化。西風(fēng)帶增強(qiáng)南移,沿岸東風(fēng)南退并減弱,由局地東風(fēng)驅(qū)動(dòng)的沿岸淡水平流減少[15],也將使得沿岸的表層水體鹽度增大。

    目前,有關(guān)南極印度洋扇區(qū)海洋水團(tuán)性質(zhì)和劃分方面已取得重要進(jìn)展,而在鹽度長期變化與大尺度環(huán)流異常的聯(lián)系等方面仍需深入。本文探討了南極印度洋扇區(qū)長期持續(xù)的海表面高鹽異?,F(xiàn)象與AAO、IOD 等環(huán)流異常的聯(lián)系,主要利用實(shí)測CTD 和EN4數(shù)據(jù)詳細(xì)分析了該海區(qū)高鹽異常的時(shí)空分布特征,進(jìn)而從AAO、IOD 影響下的大尺度環(huán)流和局地過程的角度,討論這一海表面高鹽異常的影響機(jī)理。

    2 數(shù)據(jù)

    2.1 鹽度數(shù)據(jù)

    EN4 數(shù)據(jù)[19]為經(jīng)過質(zhì)量控制的實(shí)測海洋溫鹽剖面再分析數(shù)據(jù),包含Argo、WOD09、ASBO(Arctic Synoptic Basin-wide Oceanography)以及美國GTSPP(Global Temperature and Salinity Profile Programme)實(shí)測廓線。由英國氣象局哈德利氣候變化中心(Met Office Hadley Centre for Climate Change)提供,目前版本為EN4.2.1。本文采用的是其全球海洋溫鹽場月均數(shù)據(jù),空間分辨率為1°×1°。根據(jù)Durack 和Wijffels[3]計(jì)算全球溫鹽趨勢的研究,并考慮了數(shù)據(jù)網(wǎng)站對2019-2020 年EN4鹽度數(shù)據(jù)的使用提示,本文選定南極印度洋扇區(qū)(60°~71°S,30°~120°E)為研究海域,選取1950-2018 年的鹽度資料進(jìn)行分析。

    WOD18(World Ocean Database 2018)[20]數(shù)據(jù)由政府間海洋學(xué)委員會(huì)(IOC)下設(shè)的“國際海洋學(xué)數(shù)據(jù)和信息交換”計(jì)劃(IODE)和美國國家海洋和大氣管理局(NOAA)的國家環(huán)境信息中心(NCEI)共同提供(以下稱WOD 數(shù)據(jù)),經(jīng)月平均處理后,得到該海域1981-2018 年間的鹽度時(shí)間序列。

    海豹數(shù)據(jù)[21]由MEOP(Marine Mammals Exploring the Oceans Pole to Pole)收集并發(fā)布(以下稱MEOP 數(shù)據(jù)),根據(jù)安置在海豹身上的CTD 等儀器,獲得各海域的物理、生物、地質(zhì)、化學(xué)等海洋變量的實(shí)測資料。本文采用澳大利亞海洋集成觀測系統(tǒng)(Integrated Marine Observing System,IMOS)、德國PANGAEA 數(shù)據(jù)中心以及法國提供的南極印度洋扇區(qū)海豹數(shù)據(jù)資料,選取南極印度洋扇區(qū)中溫鹽結(jié)構(gòu)較為合理、具有代表性的溫鹽廓線,得到該海區(qū)2004-2018 年溫鹽數(shù)據(jù)。

    2.2 大氣數(shù)據(jù)

    本文采用的大氣再分析數(shù)據(jù)為哥白尼氣候變化服務(wù)中心提供的ERA5 全球再分析產(chǎn)品[22],數(shù)據(jù)包括1978 年至今的大氣變量,比ERA-Interim 覆蓋了更長的時(shí)段、提供更多的變量且具有更高的時(shí)空分辨率,數(shù)據(jù)網(wǎng)址為:https://cds.climate.copernicus.eu/。本文使用月均的10 m 風(fēng)速、海平面氣壓、蒸發(fā)和降水?dāng)?shù)據(jù),數(shù)據(jù)空間分辨率為0.25°×0.25°,時(shí)間范圍為1979-2018 年。此外,大尺度環(huán)流指數(shù)包括AAO 指數(shù)和IOD 指數(shù)[23],它們分別由NOAA 下設(shè)的氣候預(yù)測中心(Climate Prediction Center,CPC)、全球氣候觀測系統(tǒng)(Global Climate Observing System,GCOS)提供。

    2.3 海冰數(shù)據(jù)

    海冰密集度數(shù)據(jù)[24]由美國國家冰雪數(shù)據(jù)中心(NSIDC,https://nsidc.org/)提供,本文采用Passive Microwave Sea Ice Concentration v3 的月均數(shù)據(jù)。該數(shù)據(jù)的空間分辨率為25 km×25 km,時(shí)間范圍為1980 年1 月至2018 年12 月。目前,與其他兩種海冰密集度算法的結(jié)果相比,基于該版本數(shù)據(jù)得到的南極海冰密集度與實(shí)測數(shù)據(jù)相比,誤差最小[25]。本文選取了1980-2018年的海冰密集度月均數(shù)據(jù)。

    3 南極印度洋扇區(qū)鹽度變化

    3.1 表層分布

    印度洋扇區(qū)是南極海表面鹽度長期變化的關(guān)鍵區(qū)之一。60°S 以南的南極海表面鹽度EOF 分析的第一、第二模態(tài)及時(shí)間序列,均體現(xiàn)了印度洋扇區(qū)高鹽現(xiàn)象的出現(xiàn)與持續(xù)(圖1)。第一模態(tài)方差貢獻(xiàn)率約23.1%,1950-2018 年南極海表面鹽度變化基本一致,總體而言東南極比西南極變化更顯著。其中,南極印度洋扇區(qū)(30°~120°E)的鹽度變化最為顯著,南極半島西側(cè)的別林斯高晉海與威德爾海以西海域的鹽度變化次之。結(jié)合時(shí)間序列的變化可以發(fā)現(xiàn),在2008年之前,60°S 以南的海表面鹽度變化較小,季節(jié)變化并不明顯;此后,南極表層鹽度由負(fù)異常逐漸轉(zhuǎn)變?yōu)轱@著的正異常。2008-2018 年,時(shí)間序列則表現(xiàn)為持續(xù)鹽度正異常疊加以顯著的季節(jié)變化(圖1a)。濾掉季節(jié)變化后,這種持續(xù)顯著的高鹽異?,F(xiàn)象仍然存在,可能與大尺度環(huán)流異常有關(guān)(4.1 節(jié))。

    圖1 1950-2018 年60°S 以南的EN4 海表面鹽度EOF 分析結(jié)果Fig.1 The first two EOFs of EN4 sea surface salinity in the south of 60°S from 1950 to 2018

    海表面鹽度的偶極子分布是南極鹽度變化的重要特征之一,主要出現(xiàn)在第二模態(tài),其方差貢獻(xiàn)率為10.5%。這種分布特征可能與南極海冰偶極子(Antarctic Dipole,ADP)[26]密切相關(guān),也可能與兩個(gè)海區(qū)存在深層水的上涌有關(guān)。這種鹽度偶極子分布與EOF1的關(guān)鍵區(qū)明顯不同,表現(xiàn)為南極印度洋扇區(qū)與羅斯海的海表面鹽度呈同相變化,而與南極半島周邊海域的海表面鹽度則為反相變化。2008-2018 年印度洋扇區(qū)與羅斯海均呈現(xiàn)顯著的增鹽趨勢,也有利于南極海表面的高鹽異常持續(xù)存在。EOF 前兩個(gè)模態(tài)均顯示,2008-2018 年60°S 以南海表面鹽度表現(xiàn)出明顯的季節(jié)振蕩增強(qiáng)。

    南極印度洋扇區(qū)的實(shí)測資料也顯示2008 年后持續(xù)顯著的海表面高鹽異常(圖2a)。EN4 數(shù)據(jù)與實(shí)測資料的表現(xiàn)基本一致,低鹽特征在2008 年以后幾乎發(fā)生了逆轉(zhuǎn)。依據(jù)這種鹽度變化特征,本文以1979-2002 年、2008-2018 年兩個(gè)時(shí)段(通過置信度為99%的均值檢驗(yàn)),重點(diǎn)探討該海域鹽度的長期變化。小波分析的結(jié)果同樣也顯示,表層鹽度在前一時(shí)間段主要表現(xiàn)為無明顯周期性的低鹽特征,而在近年來(2008-2018 年)表現(xiàn)出季節(jié)性和年際變化相疊加的高鹽特征,兩個(gè)時(shí)段間表層鹽度增加了0.2。2008 年開始持續(xù)性地出現(xiàn)12~16 個(gè)月、32 個(gè)月的顯著周期(通過置信度為95%的顯著性檢驗(yàn),圖2b)。兩個(gè)時(shí)段的WOD 和MEOP 資料,也證實(shí)了該海域前期低鹽、后期高鹽的特征,其鹽度變化可達(dá)0.8。我們同時(shí)也發(fā)現(xiàn),海豹MEOP 數(shù)據(jù)所表現(xiàn)出的高鹽現(xiàn)象更為顯著,這可能是由于海豹記錄時(shí)間與生活習(xí)性有關(guān)。而2008 年后的WOD 數(shù)據(jù)主要集中在低鹽的夏季,數(shù)據(jù)連續(xù)性明顯不如海豹MEOP 數(shù)據(jù),故該時(shí)段的實(shí)測數(shù)據(jù)主要參考了MEOP 數(shù)據(jù)。

    圖2 南極印度洋扇區(qū)表層鹽度變化(a)及其小波分析(b)Fig.2 The low frequency change of sea surface salinity derived from the observed and reanalysis data (a) and the wavelet analysis (b)in the Indian sector of the Southern Ocean (ISO)

    南極印度洋扇區(qū)的表層高鹽異常主要集中在南極沿岸,MEOP、WOD 數(shù)據(jù)實(shí)測數(shù)據(jù)更為清晰的呈現(xiàn)出近年來的高鹽異常特征(圖3)。2008-2018 年表層鹽度為33.9±0.09,在埃默里冰架、沙克爾頓冰架及鄰近冰間湖附近海域出現(xiàn)顯著的高鹽特征,如沙克爾頓冰架(Shackleton Ice Shelf,SIS)及其以北海域(95°~115°E)出現(xiàn)鹽度高于34.5 的高鹽冰架水;達(dá)恩利冰間湖(Darnley Polynya,DP)附近(62°~81°E)也存在這一高鹽水,但分布范圍較前者而言更為集中。相較而言,前一時(shí)間段內(nèi)的絕大部分海域表層海水鹽度低于33.9,埃默里冰架前緣出現(xiàn)鹽度低至33.4 的低溫低鹽的普里茲灣表層水,僅在沿岸的達(dá)恩利冰間湖附近海域存在高鹽水體,且較2008-2018 年更為集中也鹽度偏低[11,18]。

    圖3 MEOP、WOD 與EN4 數(shù)據(jù)給出的1979-2002 年(a)、2008-2018 年(b)的南極印度洋扇區(qū)表層鹽度分布Fig.3 Surface salinity distribution derived from MEOP,WOD and EN4 in the Indian sector of the Antarctic during 1979-2002 (a) and 2008-2018 (b) respectively

    對比兩個(gè)時(shí)段EN4 數(shù)據(jù)鹽度的空間分布可以看出,南極印度洋扇區(qū)的海表面鹽度在沿岸明顯較大。與前一個(gè)時(shí)間段相比,2008-2018 年33.9 等鹽線由南極沿岸向北退縮,高鹽水體基本占據(jù)了65°S 以南的南極印度洋扇區(qū)沿岸海域(圖3b中的灰實(shí)線)。95°E 以西、65°S 以北的深海平原是低鹽表層水僅存的影響區(qū)域。前一個(gè)時(shí)間段的等值線以60°E 的沿岸為中心向外海延伸,鹽度逐漸變小。在此期間,海表面鹽度基本低于33.9。在后一個(gè)時(shí)間段內(nèi),等鹽線以60°~70°E、110°E 的沿岸海域?yàn)橹行南蛲夂Q由?,中心沿岸海域鹽度高于33.9,體現(xiàn)出南極印度洋扇區(qū)的海表面鹽度增大、高鹽范圍增加。

    3.2 垂直結(jié)構(gòu)

    鹽度垂直結(jié)構(gòu)的變化主要集中在200 m 以淺的上層海洋(圖4a)??紤]到2008-2018 年表層鹽度季節(jié)變化增強(qiáng)(圖1),本文選擇海表面鹽度年變化顯著且與海冰季節(jié)性生消密切相關(guān)的冬季(7 月)、夏季(2 月)的實(shí)測WOD 與MEOP 數(shù)據(jù)的鹽度廓線,兼顧了觀測數(shù)據(jù)點(diǎn)盡可能多、覆蓋的時(shí)間也較長,并且能夠反映鹽度垂直廓線在1979-2018 年的逐年變化。鹽度在不同季節(jié)的年際變化表明,不論是海冰生成、海表面鹽度較高的冬季(圖4b),還是在海冰消融、海表面鹽度較低的夏季(圖4c),南極印度洋扇區(qū)表層及次表層的鹽度在2008 年后均表現(xiàn)出明顯的高鹽異常。

    2008-2018 年的上層海洋鹽度比1979-2002 年平均鹽度明顯增大(圖4a),增鹽現(xiàn)象在50 m 以淺尤為顯著。鹽躍層深度約為100 m,躍層以上水體的鹽度隨深度變化較小。主要的增鹽區(qū)域集中在62°~81°E之間的沿岸海域(圖3),這里分布著達(dá)恩利冰間湖、麥肯齊冰間湖和普里茲灣冰間湖,高鹽異常與冰間湖的結(jié)冰析鹽過程密切相關(guān)[16]。在海冰生成期間,100 m以淺的海水均存在鹽度明顯增大現(xiàn)象。此時(shí)躍層較深,上混合層厚度可達(dá)120 m。表層水體在海冰融化期間則表現(xiàn)為鹽度降低(圖4c),上層海洋對流明顯減弱,躍層甚至淺于50 m,并且150 m 以深的鹽度亦有所增加。年均鹽度廓線清晰的呈現(xiàn)出這一高鹽異常的垂直結(jié)構(gòu)變化,與1979-2002 年相比,2008-2018年的平均鹽度的高鹽異常主要發(fā)生在0~50 m 的表層、150~500 m 的深層,而在次表層50~150 m 深度內(nèi)略有變淡。其中,次表層的變淡可能與史久新[27]所概述的變性繞極深層水?dāng)y帶的熱量入侵冰腔所導(dǎo)致的底部融化有關(guān)??傮w說來,躍層在2008-2018 年明顯減弱,這可能與深層水變化有關(guān)。

    圖4 南極印度洋扇區(qū)兩個(gè)時(shí)段的實(shí)測年均鹽度(a)、各年7 月(b)和2 月(c)鹽度的垂直廓線Fig.4 The mean salinity profile observations in ISO during two time periods (a),and mean salinity profile in July (b),February(c) of each year

    南極印度洋扇區(qū)鹽度經(jīng)向斷面在前后兩個(gè)時(shí)間段的對比,體現(xiàn)了高鹽陸架水的向北擴(kuò)張,以及繞極深層水涌升的增強(qiáng)。圖5給出了兩個(gè)時(shí)段經(jīng)向斷面的鹽度(陰影)、溫度(等值線)結(jié)構(gòu)。黑色實(shí)線(34.2 等鹽線)指示的是南極印度洋扇區(qū)冰間湖結(jié)冰析鹽產(chǎn)生對流,使得表層水與深層高鹽水發(fā)生混合的等鹽面[16]。在前一時(shí)間段內(nèi),低溫低鹽的水體廣泛存在于60°S 以南、50 m 以淺的近表層,50 m 以淺的水體分層明顯。南極幅散帶附近、鹽度大于34.5 的繞極深層水集中于150 m 以深,其北側(cè)的水體鹽度較低而南側(cè)高。2008-2018 年期間近表層水體增鹽明顯。繞極深層水上涌以63°S 為核心,更為集中并向南入侵,出現(xiàn)高鹽的繞極深層水顯著入侵陸架。34.2 等鹽線爬升至50 m 以淺,最淺可達(dá)20 m。Guo等[28]在普里茲灣的研究也發(fā)現(xiàn)了變性繞極深層水涌升入侵陸架,上層低鹽水明顯向北堆積,鹽度低于34.2 的水層明顯變薄,表明表層水與深層水混合可以發(fā)生在更淺深度上。100 m 深度附近的繞極深層水向上層海洋的入侵更明顯。此外,北擴(kuò)的高鹽陸架水與南侵的繞極深層水之間的溫鹽梯度也變得更大。

    圖5 1979-2002 年(a)和2008-2018 年(b)南極印度洋扇區(qū)EN4 經(jīng)向溫(等值線,單位:°C)鹽斷面Fig.5 Meridional mean salinity and temperature(isoline,unit:°C) sections in ISO (30°-120°E) during 1979-2002 (a) and 2008-2018 (b)

    從整個(gè)斷面來看,2008-2018 年間高鹽陸架水體積顯著增大,繞極深層水上涌增強(qiáng),近岸(65°S 以南)上層水體混合更加均勻。北擴(kuò)的高鹽陸架水直接影響了印度洋扇區(qū)高鹽異常的水平分布與垂直結(jié)構(gòu),而繞極深層水的涌升對表層水體增鹽也有很大促進(jìn)作用[28]。以上高鹽異?,F(xiàn)象出現(xiàn)及持續(xù)時(shí)間與近年來AAO、IOD 正位相密切相關(guān),為深入探究與兩者的可能聯(lián)系,下文著重從大尺度環(huán)流異常和局地過程的角度進(jìn)行討論和分析。

    4 大尺度環(huán)流的作用

    4.1 AAO、IOD 的影響

    AAO、IOD 正位相期間南極印度洋扇區(qū)海表面表現(xiàn)出顯著的高鹽異常(圖6)。依據(jù)AAO 和IOD 指數(shù),確定大尺度環(huán)流異常(超出1 倍方差)的顯著年份(表1)。大尺度環(huán)流異常年份下表層鹽度異常(SSSa)的合成分析與EOF 分析結(jié)果(圖1)相似,南極印度洋扇區(qū)、南極半島-威德爾海附近海區(qū)在AAO、IOD 正位相期間均出現(xiàn)顯著的高鹽異常,羅斯海及其鄰近海域則表現(xiàn)為鹽度負(fù)異常。AAO 與IOD 負(fù)位相期間,南極印度洋扇區(qū)也呈現(xiàn)出顯著的鹽度負(fù)異常。環(huán)流指數(shù)的時(shí)間序列顯示,自2008 年以來AAO 與IOD 主要表現(xiàn)為正位相(表1),是導(dǎo)致南極印度洋扇區(qū)的高鹽異常長期存在(圖2)的主要大尺度環(huán)流背景。為探究大尺度環(huán)流系統(tǒng)對高鹽異常的影響,本文分析了AAO、IOD 指數(shù)與南極印度洋扇區(qū)SSSa 時(shí)間序列的超前滯后相關(guān)(圖略),結(jié)果表明對于SSSa 低頻變化,特別是年際尺度上的變化,AAO 超前1~2 年時(shí)正相關(guān)系數(shù)最大(r≥0.7),而IOD 超 前1~2 年和6~10 年有顯著正相關(guān)(r=0.6~0.7),相關(guān)系數(shù)均通過95%的置信檢驗(yàn),這在某種程度上表明AAO 較IOD 的影響更顯著。

    圖6 1979-2018 年間AAO、IOD 指數(shù)序列與環(huán)流異常期間的南極海表面鹽度異常合成分析結(jié)果Fig.6 The AAO and IOD indices and the corresponding composite analysis of sea surface salinity anomaly in the Antarctic from 1979 to 2018

    表1 1979-2018 年間AAO、IOD 指數(shù)的環(huán)流異常年份Table 1 The anomaly years of AAO and IOD indices from 1979 to 2018

    大尺度環(huán)流異常(AAO 與IOD 正位相)主要通過風(fēng)場旋度、海冰生消影響南極印度洋扇區(qū)的海表面鹽度變化。本文利用多年逐月的大氣要素(包括風(fēng)、氣壓、蒸發(fā)和降水)和海冰密集度進(jìn)行合成分析用以探究AAO 和IOD 正、負(fù)位相的不同響應(yīng)過程。圖7等值線表示海表面氣壓異常(實(shí)線代表正異常,虛線代表負(fù)異常),陰影表示風(fēng)場旋度異常。AAO 與IOD正位相時(shí)60°S 以南的南大洋表現(xiàn)為較為一致的氣壓場和風(fēng)場旋度的響應(yīng),體現(xiàn)為強(qiáng)勁的風(fēng)場旋度負(fù)異常與低壓異常,這種影響在冬季(6-8 月)體現(xiàn)得更加明顯(圖7),以AAO 正位相時(shí)的南極印度洋扇區(qū)最為顯著。風(fēng)場旋度異常通過增強(qiáng)埃克曼抽吸,促進(jìn)繞極深層水涌升,有利于對流混合過程產(chǎn)生海表面的高鹽異常。這種作用在靜力穩(wěn)定度較差的南極印度洋扇區(qū)影響更加明顯[29]。海冰異常的響應(yīng)主要發(fā)生在秋季(3-5 月)海冰生成期間,因秋冬季風(fēng)場異常無明顯變化,故圖8將秋季海冰和冬季風(fēng)場異常進(jìn)行合成分析。AAO、IOD 正位相期間,西風(fēng)帶向極移動(dòng),南極沿岸東風(fēng)向南退縮,南極印度洋扇區(qū)沿岸出現(xiàn)強(qiáng)勁的西風(fēng)異常,促進(jìn)海表面氣溫降低,有利于海冰的生成、結(jié)冰析鹽產(chǎn)生的高鹽水的向北輸運(yùn)[30]??v觀60°S以南的南大洋,威德爾海對AAO、IOD 正位相的響應(yīng)過程與印度洋扇區(qū)相近,在羅斯海則明顯不同,尤其是IOD 正位相期間的秋季海冰異常。

    圖7 南極冬季海表面氣壓異常與風(fēng)場旋度異常的合成分析Fig.7 The composite analysis of sea surface pressure anomaly and wind curl anomaly in the Antarctic in winter

    圖8 南極冬季風(fēng)場異常與秋季海冰密集度異常的合成分析Fig.8 The composite analysis of winter wind anomaly and autumn sea ice concentration anomaly in the Antarctic

    比較AAO 和IOD 不同位相期間的大氣和海冰變化,也表明大尺度環(huán)流異常的響應(yīng)過程是影響印度洋扇區(qū)鹽度變化的重要因素。如AAO 與IOD 負(fù)位相期間,南極印度洋扇區(qū)的風(fēng)場與氣壓場呈顯著的正異常(圖7),不利于繞極深層水的上涌與表層結(jié)冰析鹽產(chǎn)生的鹵水進(jìn)行對流混合[11],上層水體增鹽受到阻礙。IOD 負(fù)位相期間僅在南極60°~90°E 范圍內(nèi)表現(xiàn)出弱的海冰正異常(圖8),60°S 以南的西風(fēng)也較弱,主要表現(xiàn)為東風(fēng)異常,水體向岸積聚,不利于暖而咸的深層水上涌以及混合增密過程的持續(xù)[31]。

    4.2 局地過程的響應(yīng)

    在大尺度環(huán)流異常的背景下,局地風(fēng)場剪切與蒸發(fā)的增強(qiáng)也是高鹽異常的重要影響因素。根據(jù)前文分析得到的南極印度洋扇區(qū)的典型高鹽異常區(qū)(圖3的SIS 和DP 附近海域),著重探討了SSSa 與局地風(fēng)場(緯向風(fēng)、風(fēng)場旋度)異常的滑動(dòng)相關(guān)關(guān)系(圖9)。圖中虛線指示95%的顯著性檢驗(yàn)域。

    圖9 高鹽異常區(qū)(圖3)的SSSa 與局地風(fēng)場的滑動(dòng)相關(guān)Fig.9 The running-mean correlation coefficients of zonal wind anomaly and wind curl anomaly with SSSa in the positive salinity anomaly (Figure 3 )

    滑動(dòng)相關(guān)的結(jié)果表明,2008 年后,緯向風(fēng)異常、風(fēng)場旋度異常共同促進(jìn)了高鹽異常區(qū)(SIS、DP 區(qū))的高鹽異常。在SIS 附近海域,緯向風(fēng)與高鹽異常的正相關(guān)在前一時(shí)段短暫存在,而在2008-2018 年間則顯著增強(qiáng),與此同時(shí)局地風(fēng)旋度負(fù)異常與高鹽異常表現(xiàn)出持續(xù)的負(fù)相關(guān),這種影響主要表現(xiàn)為局地風(fēng)旋度負(fù)異常下,增強(qiáng)的埃克曼抽吸對高鹽深層水涌升的促進(jìn)作用。西風(fēng)異常與局地風(fēng)旋度負(fù)異常的影響在DP附近海域也有明顯體現(xiàn),2008-2018 年間,緯向風(fēng)與SSSa 的相關(guān)性也由負(fù)變正,對高鹽異?,F(xiàn)象具有促進(jìn)作用。這是由于近年來AAO、IOD 正位相下的西風(fēng)增強(qiáng)有關(guān)。西風(fēng)風(fēng)速增大且向極移動(dòng),增強(qiáng)了水體離岸的埃克曼輸送。由于補(bǔ)償作用,沿岸高鹽深層水上涌,與表層高鹽水的混合得到促進(jìn)。同時(shí),沿岸的東風(fēng)減弱[12],減少了表層的湍流熱通量使海表面的溫度降低[32],使表層水體在對流混合過程中仍具有持續(xù)生成海冰并產(chǎn)生高鹽水的能力。

    此外,AAO 與IOD 正位相顯著增強(qiáng)了沿岸海域的蒸發(fā),對該海區(qū)海表面鹽度的增大有直接影響(圖10),相較而言降雨的變化影響不大。圖中陰影表示蒸發(fā)異常,實(shí)線(虛線)表示降水正(負(fù))異常。與1979-2002年相比,2008-2018 年南極印度洋扇區(qū)的海表面蒸發(fā)明顯增強(qiáng),尤其集中在沿岸冰間湖、SIS 附近以及65°S 以北、55°E 以東的海域。這種特征在AAO、IOD正位相下體現(xiàn)得更加明顯。

    圖10 蒸發(fā)異常與降水異常在1979-2002 年、2008-2018 年兩個(gè)時(shí)間段的平均態(tài)與AAO、IOD正位相下的合成分析Fig.10 The climatological evaporation and precipitation anomaly during the 1979-2002,2008-2018,and the composite analysis of AAO,IOD positive phases

    當(dāng)IOD 處于較強(qiáng)正位相時(shí),緯向風(fēng)的經(jīng)向梯度增強(qiáng),向極移動(dòng)的西風(fēng)與減弱的東風(fēng)使沿岸的風(fēng)場剪切增大;另一方面,從熱帶到印度洋東部的羅斯貝(Rossby)波列能夠調(diào)節(jié)海冰帶附近的氣旋活動(dòng)[13]。這兩種作用都能在南極印度洋扇區(qū)激發(fā)更多的氣旋式環(huán)流[32],促進(jìn)海表面蒸發(fā)的增強(qiáng)。同時(shí),AAO 正位相也有利于氣旋頻率增大[33],這種蒸發(fā)的增強(qiáng)應(yīng)是AAO 與IOD 的共同影響。

    5 總結(jié)與討論

    本文采用實(shí)測MEOP、WOD 資料和EN4 再分析數(shù)據(jù)對60°S 以南的南極海表面高鹽異?,F(xiàn)象進(jìn)行分析,探討大尺度環(huán)流(AAO 和IOD)異常的影響。EOF分析表明,印度洋扇區(qū)是南極海表面鹽度長期變化的關(guān)鍵區(qū)之一。2008-2018 年間,該海區(qū)海表面鹽度變化表現(xiàn)為持續(xù)的高鹽異常,并疊加以顯著的12~16個(gè)月、32 個(gè)月的變化。這種高鹽異常集中在SIS 和DP 附近海域。高鹽異常的影響深度加深,100 m 以淺的水體均存在鹽度明顯增大的現(xiàn)象,其中冬季是最顯著的?;旌蠈觾?nèi)的低鹽水向北退縮,65°S 以南的水體明顯咸化,致使高鹽陸架水的體積增大,繞極深層水涌升增強(qiáng)且顯著入侵,兩種高鹽水體間的混合變得更加均勻。

    此高鹽異常現(xiàn)象與近年來大尺度環(huán)流AAO、IOD持續(xù)的正位相及其在局地產(chǎn)生的影響有關(guān)。在AAO 與IOD 的正位相下,南極印度洋扇區(qū)表現(xiàn)出顯著的風(fēng)場旋度負(fù)異常與低壓異常,增強(qiáng)了??寺槲c高鹽繞極深層水的上涌,有利于對流混合過程產(chǎn)生海表面的高鹽異常。同時(shí),西風(fēng)顯著增強(qiáng)并向極移動(dòng),有利于海冰大量生成,發(fā)生在冰間湖與冰架底部冷冰腔的結(jié)冰析鹽過程能夠提供大量高鹽水。另一方面,增強(qiáng)的西風(fēng)使離岸的水體輸送增強(qiáng),繞極深層水在近岸涌升并入侵陸架,與表層結(jié)冰析鹽的高鹽水發(fā)生對流混合,且埃默里冰架的冷冰腔有利于高鹽陸架水的生成,并且阻礙了高溫、高鹽的繞極深層水融化冰架底部[27],有利于該高鹽現(xiàn)象的維持。AAO 與IOD 的局地影響在SIS 和DP 兩個(gè)高鹽異常區(qū)表現(xiàn)得最為顯著。2008 年后緯向風(fēng)異常與高鹽異常呈持續(xù)顯著的正相關(guān),體現(xiàn)沿岸西風(fēng)異常促使SIS 與DP 附近水體的離岸輸送,由于補(bǔ)償作用,沿岸高鹽深層水與表層高鹽水的混合得到促進(jìn)。2008-2018 年期間,AAO 與IOD 持續(xù)的正位相增大了沿岸風(fēng)場剪切,增強(qiáng)該海區(qū)的海表面蒸發(fā),對高鹽異常有直接的促進(jìn)作用。這種影響在達(dá)恩利冰間湖、麥肯齊冰間湖、普里茲灣冰間湖、沙克爾頓冰架附近海區(qū)體現(xiàn)得更加顯著。

    在南極印度洋扇區(qū),持續(xù)的強(qiáng)風(fēng)維持了海表面較長時(shí)間的低溫環(huán)境,使下層暖水上涌帶來的熱量不足以抑制海冰的生成,使結(jié)冰析鹽作用與深對流過程能夠持續(xù)為海表面提供高鹽水。然而與南極印度洋扇區(qū)變化明顯不同,羅斯海并未出現(xiàn)如此顯著的高鹽異常,這可能是由于IOD 正位相下羅斯海產(chǎn)生顯著的北風(fēng)異常,暖平流使海冰生成減弱,與AAO 正位相在此的作用相反。同時(shí),本文還探索了厄爾尼諾-南方濤動(dòng)與南極印度洋扇區(qū)的海表面鹽度異常的相關(guān)性,結(jié)果發(fā)現(xiàn)AAO、IOD 仍是主要影響因素。此外,本文發(fā)現(xiàn)海表面鹽度異常還表現(xiàn)出羅斯海與威德爾海的偶極子,這種分布特征可能與AAO 的變化有關(guān)[34],及其在全球變暖背景下其對南極表層鹽度的影響仍待深入研究。

    致謝:感謝國際MEOP 聯(lián)合會(huì)及其做出貢獻(xiàn)的國家收集并免費(fèi)提供海洋哺乳動(dòng)物數(shù)據(jù)(http://www.meop.net);感謝哥白尼氣候變化服務(wù)中心(Copernicus Climate Change Service)提供的ERA5 數(shù)據(jù)資料。

    猜你喜歡
    海冰
    基于Argo浮標(biāo)的南極海冰范圍變化分析
    末次盛冰期以來巴倫支海-喀拉海古海洋環(huán)境及海冰研究進(jìn)展
    近三十年以來熱帶大西洋增溫對南極西部冬季海冰變化的影響
    南極海冰融化致帝企鵝減少
    1979—2015年羅斯海海冰運(yùn)動(dòng)速度變化
    極地研究(2018年2期)2018-06-27 09:09:34
    基于SIFT-SVM的北冰洋海冰識(shí)別研究
    海冰,來年再見啦!
    海洋世界(2016年4期)2016-06-15 01:51:28
    累積海冰密集度及其在認(rèn)識(shí)北極海冰快速變化的作用
    應(yīng)用MODIS數(shù)據(jù)監(jiān)測河北省近海海域海冰
    河北遙感(2014年4期)2014-07-10 13:54:59
    基于TerraSAR-X全極化數(shù)據(jù)的北極地區(qū)海冰信息提取
    亚洲色图av天堂| 人体艺术视频欧美日本| 亚洲欧美成人精品一区二区| www.av在线官网国产| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久久久久久| 亚洲自拍偷在线| 神马国产精品三级电影在线观看| 99久久中文字幕三级久久日本| 舔av片在线| 欧美色视频一区免费| 午夜精品国产一区二区电影 | a级一级毛片免费在线观看| 人人妻人人澡人人爽人人夜夜 | 波多野结衣巨乳人妻| 永久网站在线| 国产精品国产三级国产av玫瑰| 超碰av人人做人人爽久久| 午夜激情福利司机影院| 99精品在免费线老司机午夜| 国产熟女欧美一区二区| 国产精品av视频在线免费观看| 我要看日韩黄色一级片| 免费观看的影片在线观看| 少妇裸体淫交视频免费看高清| avwww免费| 中文字幕av成人在线电影| 22中文网久久字幕| 12—13女人毛片做爰片一| 日本五十路高清| 26uuu在线亚洲综合色| 国产精品免费一区二区三区在线| 色吧在线观看| av又黄又爽大尺度在线免费看 | 日韩欧美一区二区三区在线观看| 五月伊人婷婷丁香| 免费电影在线观看免费观看| 欧美一区二区国产精品久久精品| 少妇裸体淫交视频免费看高清| 在线观看美女被高潮喷水网站| 一边摸一边抽搐一进一小说| 在线免费观看不下载黄p国产| av在线播放精品| av卡一久久| 欧美xxxx性猛交bbbb| 又爽又黄a免费视频| 在线播放国产精品三级| 日韩大尺度精品在线看网址| 特级一级黄色大片| 国产精品麻豆人妻色哟哟久久 | 久久久午夜欧美精品| 久久精品国产亚洲av香蕉五月| 国产蜜桃级精品一区二区三区| 免费看a级黄色片| av天堂在线播放| 日韩制服骚丝袜av| 国产国拍精品亚洲av在线观看| 18禁在线播放成人免费| 99riav亚洲国产免费| 日韩一区二区三区影片| 亚洲精品乱码久久久v下载方式| 真实男女啪啪啪动态图| 尾随美女入室| 99精品在免费线老司机午夜| 亚洲第一区二区三区不卡| 亚州av有码| 男人的好看免费观看在线视频| 91精品国产九色| 日本熟妇午夜| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 麻豆精品久久久久久蜜桃| 亚洲成av人片在线播放无| 男人舔奶头视频| 国国产精品蜜臀av免费| 在线免费观看不下载黄p国产| 丰满乱子伦码专区| 国产精品一及| 高清毛片免费观看视频网站| 中文字幕人妻熟人妻熟丝袜美| 伦理电影大哥的女人| 精品久久久久久久久亚洲| 亚洲精品自拍成人| 看非洲黑人一级黄片| 午夜精品一区二区三区免费看| 国内精品宾馆在线| 久久午夜福利片| 午夜精品国产一区二区电影 | 国产精品一区二区三区四区久久| 国产熟女欧美一区二区| 国产高潮美女av| 国产伦精品一区二区三区视频9| 夜夜看夜夜爽夜夜摸| 成人欧美大片| 性欧美人与动物交配| 国产在线精品亚洲第一网站| 国产精品人妻久久久久久| .国产精品久久| 神马国产精品三级电影在线观看| 少妇人妻一区二区三区视频| 99在线人妻在线中文字幕| 日日摸夜夜添夜夜爱| 51国产日韩欧美| 午夜久久久久精精品| 熟妇人妻久久中文字幕3abv| 91久久精品电影网| 寂寞人妻少妇视频99o| 国产精品三级大全| 久久久久性生活片| 亚洲在久久综合| 少妇被粗大猛烈的视频| 欧美xxxx黑人xx丫x性爽| 国产在线精品亚洲第一网站| 日韩av不卡免费在线播放| 中国国产av一级| 国产精品伦人一区二区| 网址你懂的国产日韩在线| 高清日韩中文字幕在线| 精品欧美国产一区二区三| 国产精品野战在线观看| 自拍偷自拍亚洲精品老妇| 两个人视频免费观看高清| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久免费av| 99热只有精品国产| 一进一出抽搐动态| 欧美日韩国产亚洲二区| 看非洲黑人一级黄片| 成人无遮挡网站| 国内精品久久久久精免费| 天天躁夜夜躁狠狠久久av| 一级毛片久久久久久久久女| 国产久久久一区二区三区| av.在线天堂| 国产成人福利小说| 99热这里只有是精品在线观看| 成人永久免费在线观看视频| 日韩精品有码人妻一区| 99久国产av精品国产电影| 久久久久久久久中文| 国产精品无大码| 亚洲婷婷狠狠爱综合网| 亚洲三级黄色毛片| 国产午夜精品久久久久久一区二区三区| 国产精品精品国产色婷婷| 特大巨黑吊av在线直播| 日本成人三级电影网站| 有码 亚洲区| 亚洲一区高清亚洲精品| av天堂中文字幕网| 哪里可以看免费的av片| 亚洲精品456在线播放app| 国产精品1区2区在线观看.| 精品午夜福利在线看| 少妇人妻精品综合一区二区 | 久久中文看片网| 国内精品宾馆在线| 色5月婷婷丁香| 婷婷六月久久综合丁香| 久久午夜亚洲精品久久| 亚洲无线观看免费| 精品一区二区免费观看| 下体分泌物呈黄色| 伊人久久精品亚洲午夜| 国产精品国产三级国产av玫瑰| 国产精品熟女久久久久浪| 日韩成人伦理影院| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 能在线免费看毛片的网站| av专区在线播放| 午夜av观看不卡| 成人无遮挡网站| 少妇的逼水好多| 国产一区二区在线观看日韩| 免费日韩欧美在线观看| 亚洲久久久国产精品| 成人国语在线视频| a级毛色黄片| 女人久久www免费人成看片| 七月丁香在线播放| 日韩人妻高清精品专区| 亚洲精品久久午夜乱码| 中文字幕人妻熟人妻熟丝袜美| 中文字幕人妻丝袜制服| 黑人猛操日本美女一级片| 999精品在线视频| 亚洲人成77777在线视频| 啦啦啦在线观看免费高清www| 亚洲精品日本国产第一区| 亚洲无线观看免费| 菩萨蛮人人尽说江南好唐韦庄| 嫩草影院入口| 天天影视国产精品| 免费人妻精品一区二区三区视频| 观看美女的网站| av线在线观看网站| 国产免费一级a男人的天堂| 国产精品女同一区二区软件| 国产乱人偷精品视频| 纵有疾风起免费观看全集完整版| 精品一区在线观看国产| 欧美成人精品欧美一级黄| 肉色欧美久久久久久久蜜桃| 国产精品三级大全| 另类精品久久| 免费观看在线日韩| 搡老乐熟女国产| 内地一区二区视频在线| 欧美精品人与动牲交sv欧美| 91在线精品国自产拍蜜月| 丁香六月天网| 汤姆久久久久久久影院中文字幕| 人人妻人人添人人爽欧美一区卜| 丝瓜视频免费看黄片| 两个人的视频大全免费| 色视频在线一区二区三区| av黄色大香蕉| 亚洲激情五月婷婷啪啪| 亚洲婷婷狠狠爱综合网| 欧美变态另类bdsm刘玥| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| kizo精华| 九九久久精品国产亚洲av麻豆| 热99国产精品久久久久久7| 欧美激情 高清一区二区三区| 成人国产麻豆网| 丝袜脚勾引网站| 亚洲在久久综合| 十分钟在线观看高清视频www| 国产成人aa在线观看| 欧美激情极品国产一区二区三区 | 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 亚洲精品,欧美精品| 插阴视频在线观看视频| 久久亚洲国产成人精品v| 国产一区二区在线观看av| 免费黄色在线免费观看| 亚洲成色77777| 日韩av免费高清视频| 中文字幕制服av| 美女大奶头黄色视频| 一二三四中文在线观看免费高清| 日本91视频免费播放| 亚洲精品久久久久久婷婷小说| 一区二区三区四区激情视频| 美女中出高潮动态图| 少妇被粗大猛烈的视频| 水蜜桃什么品种好| 国产高清不卡午夜福利| 99九九线精品视频在线观看视频| 十八禁高潮呻吟视频| 亚洲激情五月婷婷啪啪| 精品久久久精品久久久| 尾随美女入室| 大码成人一级视频| 18禁动态无遮挡网站| 成人毛片60女人毛片免费| 亚洲精品av麻豆狂野| 国产成人午夜福利电影在线观看| 久久精品国产自在天天线| 日韩 亚洲 欧美在线| 一级毛片黄色毛片免费观看视频| 18禁在线播放成人免费| 精品视频人人做人人爽| 成人国产av品久久久| 亚洲av中文av极速乱| 我要看黄色一级片免费的| 国产欧美亚洲国产| 日韩强制内射视频| 美女中出高潮动态图| 99热网站在线观看| 国产亚洲一区二区精品| .国产精品久久| xxx大片免费视频| 在现免费观看毛片| 日韩欧美一区视频在线观看| 国产熟女午夜一区二区三区 | 考比视频在线观看| 国产精品偷伦视频观看了| 精品少妇久久久久久888优播| 国产69精品久久久久777片| 免费大片黄手机在线观看| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| 国产淫语在线视频| 只有这里有精品99| 美女福利国产在线| 国产日韩欧美在线精品| 男女边摸边吃奶| 高清不卡的av网站| 国产色爽女视频免费观看| 日韩人妻高清精品专区| 中文字幕av电影在线播放| 久久国产精品男人的天堂亚洲 | 亚洲美女黄色视频免费看| 水蜜桃什么品种好| 男女免费视频国产| 精品久久久精品久久久| 亚洲av二区三区四区| 日韩成人av中文字幕在线观看| 青春草视频在线免费观看| 日本wwww免费看| 999精品在线视频| 18在线观看网站| 美女国产高潮福利片在线看| av.在线天堂| 美女福利国产在线| 国产欧美亚洲国产| xxxhd国产人妻xxx| 高清视频免费观看一区二区| 在线观看国产h片| 美女中出高潮动态图| 国产精品国产三级国产av玫瑰| 美女视频免费永久观看网站| 大香蕉久久成人网| 一边摸一边做爽爽视频免费| 男的添女的下面高潮视频| 99久久精品国产国产毛片| 久久久国产精品麻豆| av在线播放精品| 人人妻人人爽人人添夜夜欢视频| 国产欧美另类精品又又久久亚洲欧美| 香蕉精品网在线| 久久99精品国语久久久| 亚洲精品第二区| 永久免费av网站大全| 免费观看的影片在线观看| 老司机影院毛片| 91久久精品国产一区二区成人| 一级毛片 在线播放| 国产精品蜜桃在线观看| 国产熟女午夜一区二区三区 | 日本午夜av视频| 久久久精品区二区三区| 免费高清在线观看视频在线观看| 亚洲美女黄色视频免费看| 美女cb高潮喷水在线观看| 夫妻性生交免费视频一级片| 热re99久久国产66热| 日韩,欧美,国产一区二区三区| 丰满乱子伦码专区| 国产在线一区二区三区精| 久久久午夜欧美精品| 国产精品人妻久久久久久| 成人亚洲欧美一区二区av| 制服人妻中文乱码| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| 人妻人人澡人人爽人人| 91精品三级在线观看| 日韩三级伦理在线观看| 精品久久久精品久久久| 人妻制服诱惑在线中文字幕| 免费久久久久久久精品成人欧美视频 | 日韩 亚洲 欧美在线| 久久这里有精品视频免费| 久久久久精品久久久久真实原创| 国产精品女同一区二区软件| 亚洲av在线观看美女高潮| 欧美亚洲 丝袜 人妻 在线| 成人漫画全彩无遮挡| 高清毛片免费看| 亚洲激情五月婷婷啪啪| 久久97久久精品| 日韩电影二区| 日本午夜av视频| av不卡在线播放| 老司机影院毛片| av天堂久久9| 成人综合一区亚洲| 国产精品偷伦视频观看了| 国产免费一区二区三区四区乱码| 国产深夜福利视频在线观看| 日本-黄色视频高清免费观看| 精品久久久久久久久av| 美女视频免费永久观看网站| 国产成人精品久久久久久| 九九久久精品国产亚洲av麻豆| 少妇猛男粗大的猛烈进出视频| 中文字幕制服av| 国产探花极品一区二区| 日韩三级伦理在线观看| 国产欧美日韩一区二区三区在线 | 国产在视频线精品| 欧美人与善性xxx| a级毛片在线看网站| 国产免费一区二区三区四区乱码| 亚洲精品久久成人aⅴ小说 | 欧美日韩av久久| 久久久久久久久久久久大奶| 少妇人妻精品综合一区二区| 欧美精品一区二区免费开放| 国产无遮挡羞羞视频在线观看| 精品酒店卫生间| 日日撸夜夜添| 亚洲中文av在线| 日韩制服骚丝袜av| 中文字幕亚洲精品专区| 国产精品无大码| 青春草国产在线视频| 18在线观看网站| 在线免费观看不下载黄p国产| 国产有黄有色有爽视频| 好男人视频免费观看在线| 性高湖久久久久久久久免费观看| 久久久久人妻精品一区果冻| 两个人免费观看高清视频| 有码 亚洲区| 国产淫语在线视频| 久久人人爽人人爽人人片va| 国产免费又黄又爽又色| 日韩制服骚丝袜av| 国产深夜福利视频在线观看| 一区二区av电影网| 高清毛片免费看| 精品国产一区二区久久| 午夜精品国产一区二区电影| 精品国产乱码久久久久久小说| a级毛片黄视频| 亚洲国产日韩一区二区| 亚洲综合色网址| 久久久久久久大尺度免费视频| 亚洲av日韩在线播放| 99九九在线精品视频| 亚洲性久久影院| 一级a做视频免费观看| 午夜福利视频在线观看免费| 久久久久精品久久久久真实原创| 青春草国产在线视频| 亚洲精品久久久久久婷婷小说| 插阴视频在线观看视频| 黄色配什么色好看| 在线观看国产h片| 成年女人在线观看亚洲视频| 国产精品国产三级国产av玫瑰| 国产无遮挡羞羞视频在线观看| 久久99一区二区三区| 能在线免费看毛片的网站| 黄色配什么色好看| 免费大片18禁| 少妇被粗大的猛进出69影院 | 精品少妇黑人巨大在线播放| 国产av国产精品国产| 国产在线视频一区二区| 美女xxoo啪啪120秒动态图| 夜夜爽夜夜爽视频| 亚洲一级一片aⅴ在线观看| 十八禁高潮呻吟视频| av福利片在线| 丝袜美足系列| 麻豆乱淫一区二区| 精品人妻熟女毛片av久久网站| 最黄视频免费看| 最后的刺客免费高清国语| 麻豆成人av视频| 中文精品一卡2卡3卡4更新| 高清毛片免费看| 成人手机av| 高清欧美精品videossex| 久久精品人人爽人人爽视色| 色视频在线一区二区三区| 9色porny在线观看| 一级爰片在线观看| 日韩精品有码人妻一区| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 国产视频内射| 亚洲av中文av极速乱| 欧美少妇被猛烈插入视频| 精品国产国语对白av| 亚洲国产日韩一区二区| 一级二级三级毛片免费看| 国产一区有黄有色的免费视频| 亚洲av欧美aⅴ国产| 亚洲成人手机| 不卡视频在线观看欧美| 一边亲一边摸免费视频| 成年av动漫网址| 街头女战士在线观看网站| 美女福利国产在线| 亚洲精品美女久久av网站| 欧美日韩视频高清一区二区三区二| 亚洲国产精品国产精品| 五月伊人婷婷丁香| 妹子高潮喷水视频| 免费高清在线观看视频在线观看| 能在线免费看毛片的网站| 国产免费现黄频在线看| 亚洲精品久久午夜乱码| 日韩在线高清观看一区二区三区| 国产一区亚洲一区在线观看| 国产精品不卡视频一区二区| 亚洲高清免费不卡视频| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲网站| 九九久久精品国产亚洲av麻豆| 天天影视国产精品| 视频中文字幕在线观看| 精品熟女少妇av免费看| 丰满饥渴人妻一区二区三| 亚洲美女视频黄频| 久久这里有精品视频免费| 亚洲色图综合在线观看| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的| 成人黄色视频免费在线看| av在线老鸭窝| 欧美性感艳星| 黑丝袜美女国产一区| 免费黄色在线免费观看| 久久久久精品久久久久真实原创| 天天影视国产精品| 久久久久久久精品精品| 久久99一区二区三区| 美女视频免费永久观看网站| 国产男人的电影天堂91| 国产成人freesex在线| 国产精品熟女久久久久浪| 欧美日韩在线观看h| 欧美精品高潮呻吟av久久| 欧美性感艳星| 免费观看在线日韩| 亚洲精品456在线播放app| 黑人高潮一二区| 99国产精品免费福利视频| 日韩免费高清中文字幕av| 久久99蜜桃精品久久| 久久久a久久爽久久v久久| 久久综合国产亚洲精品| 成年美女黄网站色视频大全免费 | 亚洲av日韩在线播放| 99热全是精品| 99久久中文字幕三级久久日本| 国产av一区二区精品久久| 国产精品女同一区二区软件| 国产欧美亚洲国产| 99九九在线精品视频| 婷婷色麻豆天堂久久| 国产乱人偷精品视频| 日韩三级伦理在线观看| 精品人妻在线不人妻| 日韩成人av中文字幕在线观看| 亚洲成人av在线免费| 亚洲精品第二区| 久久久久久伊人网av| 熟女av电影| 最近手机中文字幕大全| 欧美日韩视频精品一区| 男女无遮挡免费网站观看| 国产免费视频播放在线视频| 亚洲精品自拍成人| 精品亚洲成a人片在线观看| 亚洲情色 制服丝袜| 夜夜看夜夜爽夜夜摸| 嘟嘟电影网在线观看| 久久久国产一区二区| 大香蕉久久网| 99热6这里只有精品| 十分钟在线观看高清视频www| 亚洲精华国产精华液的使用体验| 丝袜美足系列| 免费人妻精品一区二区三区视频| 97在线人人人人妻| 在线观看美女被高潮喷水网站| 国产熟女欧美一区二区| 超碰97精品在线观看| 一边亲一边摸免费视频| 久久女婷五月综合色啪小说| 韩国高清视频一区二区三区| 中文精品一卡2卡3卡4更新| 伦理电影大哥的女人| 在线免费观看不下载黄p国产| 国产伦精品一区二区三区视频9| 午夜免费观看性视频| videossex国产| 久久99蜜桃精品久久| 美女中出高潮动态图| 丝袜美足系列| 亚洲美女黄色视频免费看| 精品视频人人做人人爽| 王馨瑶露胸无遮挡在线观看| 中文欧美无线码| 久久婷婷青草| 日韩欧美精品免费久久| 亚洲天堂av无毛| 国产av一区二区精品久久| 国产在视频线精品| 3wmmmm亚洲av在线观看| 国产又色又爽无遮挡免| 国产亚洲精品久久久com| 少妇丰满av| 久久亚洲国产成人精品v| 女性被躁到高潮视频| 亚洲美女视频黄频| 亚洲性久久影院| 欧美日韩av久久| 我的女老师完整版在线观看| 亚洲精品美女久久av网站| 少妇人妻精品综合一区二区| 久久久久视频综合| 性高湖久久久久久久久免费观看| 下体分泌物呈黄色| av免费在线看不卡| 日本黄色日本黄色录像| 亚洲精品成人av观看孕妇| 亚洲精品色激情综合| 久久久久久久久久久免费av| 欧美3d第一页| 欧美另类一区| 亚洲欧美一区二区三区黑人 |