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

    繞極深層水入侵阿蒙森海陸架的路徑及其性質(zhì)

    2021-03-05 02:13:38白雨趙亮李磊
    極地研究 2021年4期
    關(guān)鍵詞:阿蒙森海槽冰架

    白雨 趙亮 李磊

    研究論文

    繞極深層水入侵阿蒙森海陸架的路徑及其性質(zhì)

    白雨1趙亮1李磊2

    (1天津科技大學(xué)海洋與環(huán)境學(xué)院, 天津 300457;2中國海洋大學(xué)海洋與大氣學(xué)院, 山東 青島 266100)

    近幾十年來較暖的繞極深層水(Circumpolar Deep Water, CDW)不斷入侵阿蒙森海陸架, 使冰架底部融化, 導(dǎo)致阿蒙森海冰架質(zhì)量不斷損失。分析CDW入侵阿蒙森海陸架的路徑及性質(zhì)變化, 對(duì)研究冰架變薄和接地線后退具有重要意義。基于GLORYS12V1[Global Ocean (1/12)° Physical Reanalysis]再分析數(shù)據(jù), 對(duì)阿蒙森海西側(cè)、中央、東側(cè)通道的體積輸運(yùn)和熱輸運(yùn)進(jìn)行了計(jì)算, 給出了陸架上CDW溫度和鹽度的變化特征, 并分析了它們與流場間的關(guān)系。結(jié)果表明, Dotson-Getz海槽內(nèi)的CDW主要源自西側(cè)通道, Pine Island海槽內(nèi)的CDW主要源自中央通道和東側(cè)通道。由東側(cè)通道入侵的CDW溫度最高, 西側(cè)入侵的CDW溫度最低。CDW通過西側(cè)通道入侵陸架的體積輸運(yùn)和熱輸運(yùn)在時(shí)間序列上呈現(xiàn)微弱的上升趨勢(shì)。CDW通過中央通道向Pine Island海槽的體積輸運(yùn)和熱輸運(yùn)約是東側(cè)通道的2倍。CDW進(jìn)入海槽時(shí)的溫度主要受兩方面影響, 一方面與入侵陸架時(shí)的溫度有關(guān), 一方面又被后續(xù)的混合過程所控制。

    繞極深層水 阿蒙森海陸架 流場結(jié)構(gòu) 體積輸運(yùn) 熱輸運(yùn)

    0 引言

    近幾十年來南極冰蓋質(zhì)量不斷損失[1], 特別是西南極冰蓋在加速消融, 融化速率為每年97±13 Gt[2], 最大損失出現(xiàn)在阿蒙森海和別林斯高晉海[3]。南極冰蓋的融化對(duì)海平面變化具有重要影響[4-5], 僅阿蒙森海區(qū)域就有使海平面升高約1.5 m的潛力, 如果阿蒙森海的冰川加倍入海將導(dǎo)致海平面每世紀(jì)上升5 cm[6], 對(duì)世界各地的沿海地區(qū)造成威脅[7-8]。此外, 冰架融化向海洋注入的融冰水還會(huì)對(duì)全球氣候造成影響, 延緩全球氣候變暖的趨勢(shì)[9]。

    阿蒙森海范圍內(nèi)的冰架融化主要受到繞極深層水(Circumpolar Deep Water, CDW) 的影響[10-11], CDW隨著南極繞極流(Antarctic Circumpolar Current, ACC)不斷輸運(yùn), 其南部邊界(Southern Boundary of the Antarctic Circumpolar Current, SBACC)在該海域可以接近大陸邊緣[12], 溫暖的CDW沿著大陸架邊緣的潛流向東流動(dòng), 在近岸大陸架坡折處, 通過海底通道不斷上涌入侵陸架[13-14]。目前已觀測(cè)到, 在2007年CDW通過阿蒙森海東部陸架的西側(cè)通道進(jìn)入陸架, 體積通量為2.34× 105m3·s–1, 熱通量為2.8 TW[13]。通過模型模擬也再現(xiàn)了這個(gè)過程, 并且發(fā)現(xiàn)在其他位置也存在類似入侵過程[15]。在2013年通過水文數(shù)據(jù)驗(yàn)證了CDW從東側(cè)和中央通道進(jìn)入陸架, 對(duì)Pine Island冰架(Pine Island Ice Shelf, PIIS)的融化起重要貢獻(xiàn)[16]。CDW上涌入侵陸架的時(shí)間與局部風(fēng)場強(qiáng)迫有關(guān), 在東風(fēng)增強(qiáng)的時(shí)期通常伴隨較大的流入量[17]。這種入侵陸架的過程除了與風(fēng)場強(qiáng)迫有關(guān)還與海流和地形的交互作用有關(guān), 在陸架坡折處, 陸坡流受埃克曼運(yùn)輸作用可以向阿蒙森海東部陸架深處輸送深水[18]。在陸架上的CDW通常溫度較高, 在暖冰腔過程的作用之下不斷融化上覆冰架[19]。入侵性CDW的溫度被其厚度所控制[16], 而厚度呈現(xiàn)出季節(jié)性變化, 冬季比夏季約厚49 dbar[20], 這種季節(jié)性變化可以通過風(fēng)的水平方向變化和海冰漂移所引起的??寺槲F(xiàn)象所解釋[21]。

    隨著南極冰架與海洋的相互作用日益受到關(guān)注[22], 確定CDW流向南極冰架的路徑及其性質(zhì), 對(duì)研究冰架變薄和接地線后退的現(xiàn)象具有重大意義。目前, 相關(guān)研究主要通過建立模型模擬CDW的入侵過程或通過觀測(cè)數(shù)據(jù)分析CDW入侵后的溫鹽變化, 缺乏兩者之間的聯(lián)系。此外由于阿蒙森海是一個(gè)典型的數(shù)據(jù)匱乏區(qū)域[23], 在海冰的限制下幾乎不存在冬季的觀測(cè)數(shù)據(jù), 在夏季也缺乏連續(xù)觀測(cè)數(shù)據(jù), 較難對(duì)CDW入侵陸架的過程進(jìn)行連續(xù)的過程分析。因此本文采用了全球海洋資料同化再分析數(shù)據(jù)[Global Ocean (1/12)° Physical Reanalysis, GLORYS12V1]對(duì)CDW入侵阿蒙森海陸架的路徑及其性質(zhì)進(jìn)行分析與探討。

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

    本文的研究區(qū)域是阿蒙森海陸架海域, 水深數(shù)據(jù)來自GLORYS12V1, 水深分布及研究斷面如圖1所示, 阿蒙森海大陸架從東向西逐漸變窄, 平均深度約為500 m。在研究區(qū)域內(nèi)存在兩個(gè)主要的海槽系統(tǒng), 分別為Dotson-Getz海槽和Pine Island海槽, 最大深度都超過1600 m[24]。Dotson- Getz海槽內(nèi)的CDW主要來自西側(cè)通道, 該海槽內(nèi)的CDW主要影響Getz冰架(GIS)和Dotson冰架(DIS); Pine Island海槽內(nèi)的CDW主要來自中央和東側(cè)通道, 該海槽內(nèi)的CDW主要影響Thwaites冰架(TIS)、PIIS、Cosgrove冰架(CIS)和Abbot冰架(AIS)。TIT和BI分別表示Thwaites冰舌和Burke島, 等值線分別表示了500 m和700 m等深線。相關(guān)研究已證實(shí)CDW可以通過海底通道上涌入侵大陸架, 并且在b、c、f斷面處都發(fā)現(xiàn)了這種過程[14,16,25]。在本文研究中, 根據(jù)地理位置及水深結(jié)構(gòu)補(bǔ)充并分析了a、d、e和g 4個(gè)新的斷面。

    圖1 阿蒙森海CDW入侵路徑及研究斷面示意圖.綠色線段代表研究斷面, 橙色箭頭代表CDW入侵路徑示意圖, 黃色圓點(diǎn)代表驗(yàn)證斷面, 紫色方框代表海槽范圍, GIS、DIS、TIS、PIIS、CIS、AIS、TIT和BI分別表示Getz冰架、Dotson冰架、Thwaites冰架、Pine Island冰架、Cosgrove冰架、Abbot冰架、Thwaites冰舌和Burke島

    Fig.1.Schematic diagram of CDW intrusion pathways and study sections in the Amundsen Sea.The green line segment represents the study section, the orange arrow represents the invasion path of CDW, the yellow dot represents the validation section, and the purple box represents the extent of the trough.GIS, DIS, TIS, PIIS, CIS, AIS, TIT and BI represent the Getz Ice Shelf, Dotson Ice Shelf, Thwaites Ice Shelf, Pine Island Ice Shelf, Cosgrove Ice Shelf, Abbot Ice Shelf, Thwaites Ice Tongue and Burke Island

    本文采用了基于CMEMS(Copernicus Marine Environment Monitoring Service)系統(tǒng)提供的全球海洋資料同化再分析數(shù)據(jù)GLORYS12V1(數(shù)據(jù)來源: http://marine.copernicus.eu/services-portfolio/ access-to-products/), 數(shù)據(jù)水平分辨率為(1/12)°, 垂向分層為50層。該再分析系統(tǒng)使用的數(shù)值模式是NEMO(Nucleus for European Modelling of the Ocean)3.1版本, 同時(shí)耦合了海冰模式LIM2 (Louvain-la-Neuve Sea Ice Model), 其大氣強(qiáng)迫采用ERA-interm(European Centre for Medium- Range Weather Forecasts interim Re-Analysis)[26]。該模型通過降階卡爾曼濾波器同化了高程衛(wèi)星提供的沿軌道高度計(jì)數(shù)據(jù)、AVHRR(Advanced Very High Resolution Radiometer)提供的衛(wèi)星海表面溫度數(shù)據(jù)、CERSAT(Centre de Recherche et d’Exploi-tation Satellitaire)提供的海冰密集度, CMEMS CORAv4.1(Coriolis Ocean Dataset for Reanalysis)數(shù)據(jù)庫提供的現(xiàn)場溫度和鹽度的垂直剖面數(shù)據(jù)等數(shù)據(jù), 并利用3D-VAR方案對(duì)溫度和鹽度的偏差進(jìn)行修正[27]。

    在本研究中通過中性密度(γ)[28]定義了CDW,γ取28 kg·m–3定義為CDW的上邊界,γ取28.27 kg·m–3定義為CDW的下邊界[29], 上、下邊界之間的距離為CDW的厚度, 并以該層厚范圍內(nèi)的溫度、鹽度、流場分析CDW的性質(zhì)變化。在西側(cè)、中央和東側(cè)通道范圍內(nèi)的7個(gè)斷面處, CDW的體積輸運(yùn)()和熱輸運(yùn)()的計(jì)算公式如下:

    式中,表示網(wǎng)格面積,v表示斷面處流速法向分量,表示海水密度,p表示定壓比熱容,和ref分別表示入流或出流海水溫度和參考溫度, 參考溫度為冰點(diǎn)溫度取–1.9℃。

    2 數(shù)據(jù)驗(yàn)證

    本文選取的驗(yàn)證斷面位置如圖1中的黃色圓點(diǎn)所示, 采用的驗(yàn)證數(shù)據(jù)來自號(hào)考察船(ANT-XXVI/3)2010年在南極阿蒙森海獲取的實(shí)測(cè)CTD(Conductivity-Temperature-Depth)數(shù)據(jù)[30], 使用的測(cè)量設(shè)備是Seabird 911+ CTD (SN 561), 附帶有24個(gè)采水瓶[16]。如圖2所示, 阿蒙森海陸架范圍主要存在3種水團(tuán)[25]: 在300 m以淺的深度范圍主要為南極表層水(Antarctic Surface Water, AASW)和冬季水(Winter Water, WW), 在400 m以深的深度范圍主要為CDW。根據(jù)CDW性質(zhì)的差異, 又分為上繞極深層水(Upper Circumpolar Deep Water, UCDW)、下繞極深層水(Lower Circumpolar Deep Water, LCDW)和變性繞極深層水(Modified Circumpolar Deep Water, MCDW)。在實(shí)測(cè)斷面中, CDW的核心溫度為1.18℃, 核心鹽度為34.7[16]。根據(jù)最高溫度和最大鹽度分別定義了UCDW和LCDW[14], 其溫度最大值約為1.57℃,鹽度最大值約為34.73。CDW入侵陸架的過程中, 不斷與上層水團(tuán)發(fā)生混合, 逐漸轉(zhuǎn)變?yōu)橄鄬?duì)低溫的MCDW[31]。相比之下, 在GLORYS12V1中CDW的溫度最大值偏低, 約為1.29℃, 鹽度最大值與CTD數(shù)據(jù)相差不大, 約為34.74。GLORYS12V1與CTD數(shù)據(jù)的相關(guān)性較高, 溫度相關(guān)系數(shù)為0.85, 鹽度相關(guān)系數(shù)為0.92。

    圖2 驗(yàn)證斷面處的溫鹽點(diǎn)聚圖.a) CTD溫鹽點(diǎn)聚圖, b) GLORYS12V1溫鹽點(diǎn)聚圖

    Fig.2.-diagram of the validation section.a)-diagram of CTD; b)-diagram of GLORYS12V1

    如圖3所示, 在GLORYS12V1和CTD數(shù)據(jù)中, CDW上邊界位置都與1℃等值線較為接近, 在西側(cè)都略淺于1℃等值線, 在東側(cè)都略深于1℃等值線, 且CDW的平均溫度分別為(1.19±0.09)℃和(1.24±0.15)℃, 平均差值分別為0.07℃和0.13℃。雖然CDW的溫度在CTD斷面中更高, 但總體上2個(gè)斷面中CDW的溫度都呈現(xiàn)東高西低的分布特征。此外, CDW上邊界位置與34.6等值線較為接近, 在400 m以深的東側(cè)陸架底部呈現(xiàn)高鹽特征。在GLORYS12V1和CTD數(shù)據(jù)中, CDW的平均鹽度分別為34.70±0.04和34.70±0.03, 平均差值分別為0.03和0.02。綜上所述, 雖然CTD數(shù)據(jù)與GLORYS12V1數(shù)據(jù)存在一些差異, 但總體上GLORYS12V1數(shù)據(jù)可以表現(xiàn)出CDW高溫高鹽的性質(zhì)及其分布特征。

    圖3 驗(yàn)證斷面溫鹽分布圖.a) CTD溫鹽分布, b) GLORYS12V1溫鹽分布.黑色虛線代表1℃等值線, 黑色實(shí)線代表鹽度等值線, 白色實(shí)線表示CDW上邊界, γn=28 kg·m–3

    Fig.3.Temperature and salinity distribution of the validation section.a) temperature and salinity distribution of CTD; b) temperature and salinity distribution of GLORYS12V1.The dashed black line indicates the 1℃ isotherm.The solid black lines indicate the salinity contours.The solid white line indicates the upper boundary of CDW,γ=28 kg·m–3

    3 結(jié)果與討論

    根據(jù)1993—2017年多年的平均流場結(jié)構(gòu), 給出西側(cè)、中央和東側(cè)3個(gè)通道斷面處的體積輸運(yùn)及熱輸運(yùn)變化, 討論了流場對(duì)溫鹽性質(zhì)變化的影響。

    3.1 繞極深層水入侵路徑及性質(zhì)

    CDW上涌入侵陸架的物理機(jī)制主要包括在陸架上由冰架融化驅(qū)動(dòng)的深層流, 沿岸流與海槽的相互作用, 羅斯貝波沿著陸架坡折處與海槽的相互作用這三個(gè)方面[32]。此外, 渦流在入侵過程中也起到重要作用, 其誘導(dǎo)的輸運(yùn)主要集中在遠(yuǎn)離摩擦邊界的中間層[33]。CDW的入侵位置與地形特征有關(guān), 其入侵過程主要發(fā)生在陸架坡折處[34]。在阿蒙森海, CDW沿著ACC的路徑向東輸送, 在ACC與羅斯渦流的分叉附近(120°W~110°W之間)向南彎曲靠近阿蒙森海大陸坡, 觀測(cè)結(jié)果表明在110°W存在沿大陸架向東的潛流, 由東向正壓流誘導(dǎo)的??寺斶\(yùn)將115°W以東的暖深層水抽吸送入海槽之中[18]。在南極陸坡流附近的一些細(xì)微的地形變化也會(huì)影響其上涌過程, 如在海床階梯狀下沉的區(qū)域也會(huì)促進(jìn)上涌過程[35], 而這種入侵過程的季節(jié)性變化還與南極陸坡流的季節(jié)性周期有關(guān)[36]。

    CDW入侵陸架的過程如圖1所示。從西側(cè)通道入侵陸架的CDW主要向Dotson-Getz海槽附近輸運(yùn)。從中央通道入侵陸架的CDW在科氏力的作用下[16]會(huì)先向東輸運(yùn), 與東側(cè)通道d斷面入侵陸架的CDW發(fā)生混合, 混合后的CDW經(jīng)過g斷面繼續(xù)向南延伸進(jìn)入Pine Island海槽。從東部通道的e、f斷面入侵陸架的CDW主要向東輸運(yùn), 進(jìn)入別林斯高晉海陸架。CDW進(jìn)入海槽后, 由于其溫度相對(duì)較高, 會(huì)不斷造成冰架底部融化[37]。從圖4a可以看出, CDW上邊界位置在a和g溫度斷面上與1℃等值線基本重疊, 在b~f溫度斷面上深于1℃等值線。從圖4b可以看出, CDW上邊界位置均略深于34.6等值線, 且在b~g斷面存在鹽度大于34.7的高鹽CDW。從圖4c可以看出, CDW在b斷面處的流速最快, 其核心區(qū)域流速超過0.02 m·s–1。在c、d斷面處, CDW主要從斷面中部流入陸架, 從兩側(cè)流出陸架, 在a、f和g斷面處, CDW主要從東側(cè)進(jìn)入陸架, 從西側(cè)流出陸架, 其核心區(qū)域流速均達(dá)到0.01 m·s–1。

    如圖5所示, 自1993—2017年CDW的最大溫度在a~c斷面處呈現(xiàn)先升高后降低的趨勢(shì)變化, 在d~g斷面呈現(xiàn)降低的趨勢(shì)變化。CDW的鹽度變化較小, 但在f斷面處的鹽度明顯增加, 且在2015年存在鹽度最大值, 約為34.82。從圖5h中可以看出, 25年平均的溫鹽性質(zhì)表明CDW在f斷面處的溫度和鹽度最高, 在7個(gè)斷面中CDW的溫度范圍約為1.07~1.43℃, 鹽度范圍約為34.64~34.74。

    3.2 繞極深層水體積輸運(yùn)及熱輸運(yùn)變化

    在西側(cè)通道(a、b斷面)處, 1993至2017年期間CDW經(jīng)過a、b斷面的平均體積輸運(yùn)分別為2.16×104m3·s–1、2.55×104m3·s–1, 在時(shí)間序列上呈現(xiàn)微弱的上升趨勢(shì),值均小于0.01, 表明CDW向Dotson-Getz海槽中的輸運(yùn)量不斷增加。在CDW不斷入侵陸架的過程中攜帶大量的熱量, 經(jīng)過a、b斷面的熱輸運(yùn)分別為0.27 TW、0.33 TW。根據(jù)表1所示, b斷面處的體積輸運(yùn)方差及熱輸運(yùn)方差均為7個(gè)斷面中的最低值, 表明該斷面處的年際差異最小。

    Fig.4.Vertical sections of average temperature, salinity and normal velocity distribution from 1993 to 2017.a) temperature; b) salinity; c) normal velocity distribution.The white contours indicate the upper boundary of CDW,γ=28 kg·m–3.Positive value in c) indicates the inflow to the Amundsen Sea shelf

    圖5 1993—2017年斷面處的溫鹽點(diǎn)聚圖.a) 斷面a; b) 斷面b; c) 斷面c; d) 斷面d; e) 斷面e; f) 斷面f; g) 斷面g; h) 25年平均

    Fig.5.-diagram of the vertical section from 1993 to 2017.a) section a; b) section b; c) section c; d) section d; e) section e; f) section f; g) section g; h) average data of 25 years

    表1 1993—2017年a~g斷面處CDW體積輸運(yùn)、熱輸運(yùn)、溫度及鹽度的平均值及方差

    在中央通道(c斷面)和東側(cè)通道(d~f斷面)處, 1993至2017年期間CDW經(jīng)過c~f斷面的平均體積輸運(yùn)和熱輸運(yùn)分別為4.28×104m3·s–1、2.15× 104m3·s–1、1.83×104m3·s–1、6.49×104m3·s–1和0.57 TW、0.28 TW、0.25 TW、0.87 TW。Pine Island海槽中的CDW主要源自于中央和東側(cè)通道的c、d斷面, 其中c斷面的體積輸運(yùn)和熱輸運(yùn)約是d斷面的2倍, 并且CDW經(jīng)過g斷面向Pine Island海槽輸送的總體積輸運(yùn)和熱輸運(yùn)約是c、d斷面之和, 在輸運(yùn)過程中由c、d斷面向g斷面流動(dòng)過程中會(huì)產(chǎn)生小部分支流損耗, 此外由東側(cè)通道入侵的少量CDW也會(huì)向南輸運(yùn)經(jīng)過g斷面。由于d斷面和g斷面位置較為接近, 在2個(gè)斷面間的輸運(yùn)過程中產(chǎn)生的支流損耗較少, 因此在年際變化上, d斷面與g斷面的相關(guān)性較高, 相關(guān)系數(shù)為0.84。在1993年至2017年期間CDW經(jīng)過g斷面進(jìn)入Pine Island海槽的平均體積輸運(yùn)和熱輸運(yùn)分別為5.38×104m3·s–1和0.70 TW。在2000至2007年期間, 經(jīng)過e斷面處的體積輸運(yùn)及熱輸運(yùn)呈現(xiàn)顯著的下降趨勢(shì),值小于0.01。在1993至2017年期間, 體積輸運(yùn)和熱輸運(yùn)的峰值均出現(xiàn)在2000年的f斷面處, 并且在大多數(shù)年份上, f斷面的體積輸運(yùn)和熱輸運(yùn)均大于其他斷面, 經(jīng)f斷面入侵到陸架的CDW主要流向別林斯高晉海陸架, 為別林斯高晉海的冰川融化提供了較暖的CDW。雖然經(jīng)過西側(cè)、中央和東側(cè)通道的CDW溫度存在一定的差異(圖6c), 但熱輸運(yùn)的變化與體積輸運(yùn)的變化趨勢(shì)一致(圖6a、6b), 相關(guān)系數(shù)均大于0.99。這與Assmann等[38]的結(jié)論一致, 熱輸運(yùn)主要受流速控制。

    圖6 a–g斷面處CDW的體積輸運(yùn)、熱輸運(yùn)及溫度變化情況。a)體積輸運(yùn); b)熱輸運(yùn); c)溫度變化

    Fig.6.Variations of volume transport, heat transport and temperature of CDW in section a to g.a) volume transport; b) heat transport; c) temperature

    3.3 繞極深層水厚度及溫鹽分布特征

    1993年至2017年期間, 阿蒙森海陸架范圍內(nèi)的CDW平均厚度被水深所限制, 其厚度空間分布情況與水深分布較為接近。如圖7所示, 在71°S~72°S附近, CDW經(jīng)過f斷面沿著700 m等深線向西南方向延伸的海域厚度明顯較大, 在該區(qū)域中厚度普遍大于500 m, 最大厚度可達(dá)1200 m, 此外在100°W以東的海域中同樣存在較厚的CDW層, 最大厚度約1000 m。在冰架附近的海槽中, CDW厚度也相對(duì)較大, 在Pine Island海槽中較厚的CDW分布范圍廣, 最厚約1100 m, 在Dotson-Getz海槽中的CDW厚度雖然比Pine Island海槽薄, 但最厚處也超過600 m。

    在大尺度環(huán)流與中緯度CDW性質(zhì)變化的共同控制下[15], CDW入侵陸架時(shí)的溫度存在一定的差異(圖6c和圖8a)。經(jīng)過西側(cè)通道斷面處的平均溫度最低, 約1.18℃; 中央通道斷面處的溫度約1.25℃; 東側(cè)通道斷面處的溫度最高, 約1.32℃。觀測(cè)數(shù)據(jù)也表明從東側(cè)通道入侵的CDW溫度高于中央通道, 這種溫度差異還與陸架邊緣CDW入侵的厚度有關(guān)[16], 而陸坡流的??寺斶\(yùn)過程又會(huì)影響入侵厚度[18]。阿蒙森海陸架范圍內(nèi)CDW的溫度和鹽度分布特征較為相似(圖8), 在大陸架邊緣溫度和鹽度由北向南、由東向西逐漸降低, Dotson-Getz海槽中的溫度和鹽度明顯低于Pine Island海槽, 且最低值均出現(xiàn)在TIS和TIT西側(cè)。東側(cè)較暖的CDW向Pine Island海槽流動(dòng)的過程中, 會(huì)不斷與來自中央通道較冷的CDW發(fā)生混合, 與此同時(shí), CDW在向南輸運(yùn)的過程中也會(huì)與上層較冷的表層水或融冰水發(fā)生混合[16]。在這兩個(gè)過程的影響下, 溫度逐漸降低, 最終經(jīng)過g斷面向Pine Island海槽輸運(yùn)時(shí)的溫度約為1.18℃。

    圖7 1993—2017年阿蒙森海CDW平均厚度分布

    Fig.7.Average thickness distribution of CDW in the Amundsen Sea from 1993 to 2017

    圖8 1993—2017年阿蒙森海CDW平均溫度和鹽度分布(對(duì)CDW層內(nèi)的溫鹽進(jìn)行深度平均).a)溫度分布; b)鹽度分布

    Fig.8.Average temperature and salinity distribution of CDW in the Amundsen Sea from 1993 to 2017 (depth-averaged of the temperature and salinity within the CDW layer).a) temperature distribution; b) salinity distribution

    3.4 海槽內(nèi)繞極深層水的性質(zhì)變化

    本文以2013—2016年的變化為例, 闡述CDW溫鹽性質(zhì)與流場之間的關(guān)系, 由于溫度與鹽度變化過程基本一致(圖10和圖11), 因此僅以溫度為例闡述其與流場間的關(guān)系。如圖9和圖10所示, 在2013年TIT西北側(cè), 低溫的CDW其流向先指向西北, 在112°W附近流向轉(zhuǎn)變?yōu)槲髂? 進(jìn)入Dotson-Getz海槽范圍, 與此同時(shí)在TIT南側(cè)的冷MCDW也向西北流動(dòng), 低溫的CDW與冷MCDW沿著該路徑不斷與Dotson-Getz海槽內(nèi)相對(duì)較暖的CDW發(fā)生混合, 使其溫度不斷下降, 在2014年可以發(fā)現(xiàn)Dotson-Getz海槽內(nèi)的CDW溫度明顯變低。2015年與2016年流場分布結(jié)構(gòu)比較相似, 但2016年的流速明顯較大, 在TIT東側(cè)的冷MCDW沿著TIT邊界不斷向Pine Island海槽流動(dòng), 此外在2014年TIT的東北側(cè)也存在向Pine Island海槽輸運(yùn)的東南流。在2013—2016年期間, TIT西側(cè)的冷MCDW隨著流場不斷向外擴(kuò)散, 逐漸影響兩個(gè)海槽內(nèi)的溫鹽性質(zhì)。

    如圖11所示, 在2001年, Dotson-Getz海槽中的CDW溫度明顯升高, 并且在后續(xù)幾年內(nèi), 溫度不斷上升。在2014年, 兩個(gè)海槽中的CDW溫度明顯降低, 且在Dotson-Getz海槽內(nèi)溫度差異較大。從圖6c中可以看出, 在2001年CDW通過a、b斷面時(shí)的溫度明顯升高, 并且在后續(xù)幾年仍處于較高溫度, 在2014年經(jīng)過斷面入侵陸架的CDW溫度明顯低于其他年份, 并且g斷面的溫度變化較為顯著, 隨著流入Pine Island海槽的CDW溫度變低, Pine Island海槽內(nèi)的整體溫度也明顯變低。CDW除了入侵陸架時(shí)的溫度存在一定的差異, CDW入侵陸架后的性質(zhì)還受到流場結(jié)構(gòu)變化的影響。在溫鹽方面, 隨著流場結(jié)構(gòu)的變化, 較高溫鹽的CDW與較低溫鹽的CDW在輸運(yùn)過程中

    圖9 2013—2016年CDW的流場分布(對(duì)CDW層內(nèi)的流速進(jìn)行深度平均).a)2013年; b)2014年; c)2015年; d)2016年

    Fig.9.Flow field distribution of CDW from 2013 to 2016 (depth-averaged of the velocity within the CDW layer).a) 2013; b) 2014; c) 2015; d) 2016

    圖10 2013—2016年CDW的溫度、鹽度分布(對(duì)CDW層內(nèi)的溫鹽進(jìn)行深度平均).a) 2013年溫度分布; b) 2013年鹽度分布; c) 2014年溫度分布; d) 2014年鹽度分布; e) 2015年溫度分布; f) 2015年鹽度分布; g) 2016年溫度分布; h) 2016年鹽度分布

    Fig.10.Temperature and salinity distribution of CDW from 2013 to 2016 (depth-averaged of the temperature and salinity within the CDW layer).a) temperature distribution in 2013; b) salinity distribution in 2013; c) temperature distribution in 2014; d) salinity distribution in 2014; e) temperature distribution in 2015; f) salinity distribution in 2015; g) temperature distribution in 2016; h) salinity distribution in 2016

    圖11 Dotson-Getz海槽和Pine Island海槽內(nèi)CDW的溫度和鹽度變化.a) 溫度分布; b) 鹽度分布

    Fig.11.Variations in temperature and salinity of CDW in the Dotson-Getz Trough and Pine Island Trough.a) temperature distribution; b) salinity distribution

    不斷混合。綜上所述, CDW沿著輸運(yùn)路徑最終進(jìn)入海槽時(shí)的溫度主要受兩方面影響, 一方面與入侵陸架時(shí)的溫度有關(guān), 一方面又被后續(xù)的混合過程所控制[38]。

    4 結(jié)論

    本文基于GLORYS12V1數(shù)據(jù), 根據(jù)阿蒙森海陸架上流場結(jié)構(gòu)給出CDW入侵路徑, 計(jì)算了西側(cè)、中央和東側(cè)3個(gè)通道斷面處的體積輸運(yùn)及熱輸運(yùn), 并討論了陸架上CDW溫度和鹽度的性質(zhì)變化與流場之間的關(guān)系。

    在1993至2017年期間, CDW主要通過西側(cè)、中央和東側(cè)3個(gè)通道入侵阿蒙森海陸架, 并且入侵陸架時(shí)的溫度存在一定的差異, 由東側(cè)通道入侵的CDW溫度最高, 西側(cè)入侵的CDW溫度最低。CDW主要通過西側(cè)通道向Dotson-Getz海槽輸運(yùn), 通過中央通道和東側(cè)通道向Pine Island海槽輸運(yùn)。CDW通過西側(cè)通道入侵陸架的體積輸運(yùn)和熱輸運(yùn)在年際變化上呈現(xiàn)微弱的上升趨勢(shì), 使CDW向Dotson-Getz海槽的體積輸運(yùn)和熱輸運(yùn)不斷增加。CDW通過中央通道向Pine Island海槽的輸運(yùn)量約是東側(cè)通道的2倍。

    在阿蒙森海陸架范圍內(nèi)CDW的平均厚度空間分布與水深分布較為接近, 在Dotson-Getz海槽與Pine Island海槽厚度較大。CDW溫度的空間分布與鹽度基本一致, Pine Island海槽內(nèi)的溫度和鹽度明顯高于Dotson-Getz海槽, 且整個(gè)陸架上溫度和鹽度的最低值均出現(xiàn)在TIT和TIS西側(cè)。CDW進(jìn)入海槽時(shí)的溫度主要受兩方面影響, 一方面與入侵陸架時(shí)的溫度有關(guān), 一方面又被后續(xù)的混合過程所控制。

    1 王亞鳳, 溫家洪, 劉吉英.南極冰蓋與冰川的快速變化[J].極地研究, 2006, 18(1): 63-74.

    2 WILLIAMS S D P, MOORE P, KING M A, et al.Revisiting GRACE Antarctic ice mass trends and accelerations considering autocorrelation[J].Earth and Planetary Science Letters, 2014, 385: 12-21.

    3 RIGNOT E, JACOBS S, MOUGINOT J, et al.Ice-shelf melting around Antarctica[J].Science, 2013, 341(6143): 266-270.

    4 RIGNOT E, VELICOGNA I, VAN DEN BROEKE M R, et al.Acceleration of the contribution of the Greenland and Antarctic ice sheets to sea level rise[J].Geophysical Research Letters, 2011, 38(5): L05503.

    5 RITZ C, EDWARDS T L, DURAND G, et al.Potential sea-level rise from Antarctic ice-sheet instability constrained by observations[J].Nature, 2015, 528(7580): 115-118.

    6 VAUGHAN D G.West Antarctic Ice Sheet collapse—The fall and rise of a paradigm[J].Climatic Change, 2008, 91(1-2): 65-79.

    7 WOODROFFE C D, MURRAY-WALLACE C V.Sea-level rise and coastal change: The past as a guide to the future[J].Quaternary Science Reviews, 2012, 54: 4-11.

    8 FEAGIN R A, SHERMAN D J, GRANT W E.Coastal erosion, global sea-level rise, and the loss of sand dune plant habitats[J].Frontiers in Ecology and the Environment, 2005, 3(7): 359-364.

    9 BRONSELAER B, WINTON M, GRIFFIES S M, et al.Change in future climate due to Antarctic meltwater[J].Nature, 2018, 564(7734): 53-58.

    10 JACOBS S S, HELLMER H H, JENKINS A.Antarctic Ice Sheet melting in the southeast Pacific[J].Geophysical Research Letters, 1996, 23(9): 957-960.

    11 錢懿德, 周春霞, 陳一鳴, 等.1947—2020年西南極派恩島冰川前緣變化特征分析[J].極地研究, 2020, 32(4): 435-451.

    12 ORSI A H, WHITWORTH T III, NOWLIN W D Jr.On the meridional extent and fronts of the Antarctic Circumpolar Current[J].Deep Sea Research Part I: Oceanographic Research Papers, 1995, 42(5): 641-673.

    13 WALKER D P, BRANDON M A, JENKINS A, et al.Oceanic heat transport onto the Amundsen Sea shelf through a submarine glacial trough[J].Geophysical Research Letters, 2007, 34(2): L02602.

    14 WALKER D P, JENKINS A, ASSMANN K M, et al.Oceanographic observations at the shelf break of the Amundsen Sea, Antarctica[J].Journal of Geophysical Research: Oceans, 2013, 118(6): 2906-2918.

    15 NAKAYAMA Y, MENEMENLIS D, ZHANG H, et al.Origin of Circumpolar Deep Water intruding onto the Amundsen and Bellingshausen Sea continental shelves[J].Nature Communications, 2018, 9(1): 3403.

    16 NAKAYAMA Y, SCHR?DER M, HELLMER H H.From circumpolar deep water to the glacial meltwater plume on the eastern Amundsen Shelf[J].Deep Sea Research Part I: Oceanographic Research Papers, 2013, 77: 50-62.

    17 THOMA M, JENKINS A, HOLLAND D, et al.Modelling Circumpolar Deep Water intrusions on the Amundsen Sea continental shelf, Antarctica[J].Geophysical Research Letters, 2008, 35(18): L18602.

    18 W?HLIN A K, MUENCH R D, ARNEBORG L, et al.Some implications of Ekman layer dynamics for cross-shelf exchange in the Amundsen Sea[J].Journal of Physical Oceanography, 2012, 42(9): 1461-1474.

    19 史久新.南極冰架-海洋相互作用研究綜述[J].極地研究, 2018, 30(3): 287-302.

    20 MALLETT H K W, BOEHME L, FEDAK M, et al.Variation in the distribution and properties of Circumpolar Deep Water in the eastern Amundsen Sea, on seasonal timescales, using seal-borne tags[J].Geophysical Research Letters, 2018, 45(10): 4982-4990.

    21 KIM T W, HA H K, W?HLIN A K, et al.Is Ekman pumping responsible for the seasonal variation of warm circumpolar deep water in the Amundsen Sea?[J].Continental Shelf Research, 2017, 132: 38-48.

    22 JOUGHIN I, ALLEY R B, HOLLAND D M.Ice-sheet response to oceanic forcing[J].Science, 2012, 338(6111): 1172-1176.

    23 MEREDITH M P, DUCKLOW H W, SCHOFIELD O, et al.The interdisciplinary marine system of the Amundsen Sea, Southern Ocean: Recent advances and the need for sustained observations[J].Deep Sea Research Part II: Topical Studies in Oceanography, 2016, 123: 1-6.

    24 NITSCHE F O, JACOBS S S, LARTER R D, et al.Bathymetry of the Amundsen Sea continental shelf: Implications for geology, oceanography, and glaciology[J].Geochemistry, Geophysics, Geosystems, 2007, 8(10): Q10009.

    25 W?HLIN A K, YUAN X, BJ?RK G, et al.Inflow of warm Circumpolar Deep Water in the central Amundsen Shelf[J].Journal of Physical Oceanography, 2010, 40(6): 1427-1434.

    26 JACKSON L C, DUBOIS C, FORGET G, et al.The mean state and variability of the north Atlantic circulation: A perspective from ocean reanalyses[J].Journal of Geophysical Research: Oceans, 2019, 124(12): 9141-9170.

    27 BUONGIORNO NARDELLI B.A multi-year time series of observation-based 3D horizontal and vertical quasi-geostrophic global ocean currents[J].Earth System Science Data, 2020, 12(3): 1711-1723.

    28 JACKETT D R, MCDOUGALL T J.A neutral density variable for the world's oceans[J].Journal of Physical Oceanography, 1997, 27(2): 237-263.

    29 WHITWORTH T III, ORSI A H, KIM S J, et al.Water masses and mixing near the Antarctic slope front[M]//Ocean, Ice, and Atmosphere: Interactions at the Antarctic Continental Margin.Washington, D.C.: American Geophysical Union, 1985: 1-27.

    30 GOHL K.The expedition of the research vessel “Polarstern” to the Amundsen Sea, Antarctica, in 2010 (ANT-XXVI/3) [J].Berichte zur Polar- und Meeresforschung (Reports on Polar and Marine Research), 2010, 617.

    31 FANG L, LEE S, LEE S A, et al.Removal of refractory dissolved organic carbon in the Amundsen Sea, Antarctica[J].Scientific Reports, 2020, 10: 1213.

    32 ST-LAURENT P, KLINCK J M, DINNIMAN M S.On the role of coastal troughs in the circulation of warm Circumpolar Deep Water on Antarctic shelves[J].Journal of Physical Oceanography, 2013, 43(1): 51-64.

    33 THOMPSON A F, HEYWOOD K J, SCHMIDTKO S, et al.Eddy transport as a key component of the Antarctic overturning circulation[J].Nature Geoscience, 2014, 7(12): 879-884.

    34 KLINCK J M, DINNIMAN M S.Exchange across the shelf break at high southern latitudes[J].Ocean Science, 2010, 6(2): 513-524.

    35 LIU C Y, WANG Z M, CHENG C, et al.On the modified circumpolar deep water upwelling over the Four Ladies Bank in Prydz Bay, East Antarctica[J].Journal of Geophysical Research: Oceans, 2018, 123(11): 7819-7838.

    36 LIU C Y, WANG Z M, CHENG C, et al.Modeling modified Circumpolar Deep Water intrusions onto the Prydz Bay continental shelf, East Antarctica[J].Journal of Geophysical Research: Oceans, 2017, 122(7): 5198-5217.

    37 KIMURA S, JENKINS A, REGAN H, et al.Oceanographic controls on the variability of ice-shelf basal melting and circulation of glacial meltwater in the Amundsen Sea Embayment, Antarctica[J].Journal of Geophysical Research: Oceans, 2017, 122(12): 10131-10155.

    38 ASSMANN K M, JENKINS A, SHOOSMITH D R, et al.Variability of Circumpolar Deep Water transport onto the Amundsen Sea continental shelf through a shelf break trough[J].Journal of Geophysical Research: Oceans, 2013, 118(12): 6603-6620.

    Pathways and properties of Circumpolar Deep Water intrusion on the Amundsen Sea shelf

    Bai Yu1, Zhao Liang1, Li Lei2

    (1College of Marine and Environmental Sciences, Tianjin University of Science and Technology, Tianjin 300457, China;2College of Oceanic and Atmospheric Sciences, Ocean University of China, Qingdao 266100, China)

    The relatively warm Circumpolar Deep Water (CDW) has been intruding on the Amundsen Sea continental shelf in recent decades and melting the ice shelves from below, resulting in a continuous ice mass loss in the Amundsen Sea sector of the West Antarctic Ice Shelf.Analysis of the pathways and variations of the CDW intrusion onto the shelf is of great importance to understanding the phenomena of rapid ice shelf thinning and grounding line retreat in this sector.Based on GLORYS12V1 [Global Ocean (1/12)° Physical Reanalysis] data, we calculated the volume and heat transport in the western, central, and eastern channels.We analyzed the relationship between the variations in temperature and salinity of the CDW on the shelf and its flow field.Results showed that the CDW invaded the Dotson–Getz Trough from the western channel, while it invaded the Pine Island Trough from the central and eastern channels.The volume and heat transport of the CDW intruding on the shelf through the western channel showed a weak upward trend over time.The volume and heat transport of the CDW to the Pine Island Trough through the central channel was about twice as much as that through the eastern channel.The temperature of the CDW intrusion into these troughs was controlled by the temperature of the CDW, when it invaded the shelf, and by subsequent processes on the continental shelf.

    Circumpolar Deep Water, Amundsen Sea shelf, flow field structure, volume transport, heat transport

    2021年1月收到來稿, 2021年3月收到修改稿

    國家海洋局極地考察辦公室項(xiàng)目(RFSOCC2020-2022-No.18)資助

    白雨, 男, 1996年生。碩士研究生, 主要從事物理海洋學(xué)研究。E-mail: baiyu@mail.tust.edu.cn

    趙亮, E-mail: zhaoliang@tust.edu.cn

    10.13679/j.jdyj.20210004

    猜你喜歡
    阿蒙森海槽冰架
    南極冰架融化新原因確定
    鄂西-城口地區(qū)“海槽”西側(cè)長興組-飛仙關(guān)組臺(tái)緣礁灘特征與展布
    沖繩海槽構(gòu)造演化及其與巖漿、熱液和沉積作用的關(guān)系:研究進(jìn)展與展望
    最慢的人先到終點(diǎn)
    最慢的人先到終點(diǎn)
    最慢的人最先到終點(diǎn)
    羅斯冰架探險(xiǎn)
    遲到20年的道歉
    日本南海海槽IODP C0004C巖芯樣品熱釋光信號(hào)特征及其對(duì)古氣候變化的響應(yīng)
    沖繩海槽地殼結(jié)構(gòu)與性質(zhì)研究進(jìn)展和新認(rèn)識(shí)
    日本成人三级电影网站| 国产精品爽爽va在线观看网站| av中文乱码字幕在线| 亚洲专区国产一区二区| 人人妻,人人澡人人爽秒播| 色综合亚洲欧美另类图片| 亚洲av日韩精品久久久久久密| 亚洲精品久久国产高清桃花| 两人在一起打扑克的视频| 简卡轻食公司| 人妻少妇偷人精品九色| 亚洲欧美日韩卡通动漫| 日韩精品有码人妻一区| 成年版毛片免费区| 看十八女毛片水多多多| 久久精品91蜜桃| 欧美三级亚洲精品| 国产乱人视频| 黄色女人牲交| 欧美激情在线99| 99热6这里只有精品| 十八禁国产超污无遮挡网站| 亚洲人成网站在线播放欧美日韩| 色综合亚洲欧美另类图片| 国产精品女同一区二区软件 | 欧美绝顶高潮抽搐喷水| 99国产极品粉嫩在线观看| 欧美黑人欧美精品刺激| 禁无遮挡网站| 亚洲人成网站高清观看| 亚洲电影在线观看av| 日本 av在线| 免费看美女性在线毛片视频| 99热这里只有是精品在线观看| 精品久久久久久久久久免费视频| 免费观看在线日韩| 久久久久久久久久黄片| 嫩草影院入口| 国产精品永久免费网站| h日本视频在线播放| 国产午夜精品久久久久久一区二区三区 | 色5月婷婷丁香| 老熟妇仑乱视频hdxx| 免费av毛片视频| 日日干狠狠操夜夜爽| 久久精品国产亚洲网站| 亚洲精华国产精华液的使用体验 | 国产精品一区二区三区四区免费观看 | 无遮挡黄片免费观看| 有码 亚洲区| 高清毛片免费观看视频网站| 91麻豆av在线| 国产探花极品一区二区| 亚洲美女黄片视频| 日本免费a在线| 亚洲人成伊人成综合网2020| 欧美激情国产日韩精品一区| 91av网一区二区| 国产精品女同一区二区软件 | 99久久九九国产精品国产免费| 人人妻人人看人人澡| 真实男女啪啪啪动态图| 国产综合懂色| 成人一区二区视频在线观看| 人妻久久中文字幕网| 久久精品综合一区二区三区| av国产免费在线观看| 精品人妻1区二区| 国产乱人视频| 悠悠久久av| 精品人妻1区二区| 久久精品国产自在天天线| 香蕉av资源在线| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 国产精品爽爽va在线观看网站| 久9热在线精品视频| 真实男女啪啪啪动态图| 精品一区二区三区视频在线观看免费| 久久99热这里只有精品18| 99久久无色码亚洲精品果冻| 国产大屁股一区二区在线视频| 免费无遮挡裸体视频| 精品久久久久久,| 舔av片在线| 三级毛片av免费| 看片在线看免费视频| 97碰自拍视频| 成人午夜高清在线视频| av黄色大香蕉| 夜夜看夜夜爽夜夜摸| 国产精品,欧美在线| 国产成人影院久久av| 国产综合懂色| 日本免费一区二区三区高清不卡| 18禁黄网站禁片免费观看直播| 桃红色精品国产亚洲av| 精华霜和精华液先用哪个| 黄色女人牲交| 九色成人免费人妻av| 99久久无色码亚洲精品果冻| xxxwww97欧美| 国产精品久久电影中文字幕| 99热这里只有是精品50| 久久国产精品人妻蜜桃| 人妻丰满熟妇av一区二区三区| 看黄色毛片网站| 天天一区二区日本电影三级| 国产极品精品免费视频能看的| 91午夜精品亚洲一区二区三区 | 一本一本综合久久| 午夜精品在线福利| 国内精品久久久久久久电影| 校园春色视频在线观看| 九色国产91popny在线| 免费av观看视频| 国产乱人伦免费视频| 亚洲国产高清在线一区二区三| 婷婷精品国产亚洲av在线| eeuss影院久久| 亚洲自偷自拍三级| 亚洲最大成人av| 久久6这里有精品| 久久久久九九精品影院| x7x7x7水蜜桃| 精品午夜福利视频在线观看一区| 日本三级黄在线观看| 久久久精品大字幕| 欧美日韩国产亚洲二区| 欧美日韩综合久久久久久 | 亚洲av熟女| 亚洲av免费高清在线观看| 99热这里只有是精品在线观看| 99热这里只有精品一区| 精品久久久久久久久久免费视频| 国产亚洲精品综合一区在线观看| 在线观看免费视频日本深夜| 欧美一区二区精品小视频在线| 18禁黄网站禁片免费观看直播| 久久久成人免费电影| 狂野欧美激情性xxxx在线观看| 日韩欧美精品免费久久| a级毛片a级免费在线| 免费看美女性在线毛片视频| 日韩高清综合在线| 麻豆av噜噜一区二区三区| 国产又黄又爽又无遮挡在线| 久久精品国产鲁丝片午夜精品 | 欧美xxxx黑人xx丫x性爽| 给我免费播放毛片高清在线观看| 蜜桃亚洲精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 少妇猛男粗大的猛烈进出视频 | 人人妻人人看人人澡| 日韩在线高清观看一区二区三区 | 精品久久久久久久久亚洲 | 人妻久久中文字幕网| 国产av一区在线观看免费| 国产伦一二天堂av在线观看| 亚洲国产高清在线一区二区三| 国产女主播在线喷水免费视频网站 | 一边摸一边抽搐一进一小说| 久久久久久国产a免费观看| 婷婷六月久久综合丁香| 在线观看av片永久免费下载| 亚洲va日本ⅴa欧美va伊人久久| 免费av观看视频| 级片在线观看| 99久国产av精品| 真人一进一出gif抽搐免费| 亚洲欧美日韩卡通动漫| 精品久久久久久久久av| 久久亚洲精品不卡| 神马国产精品三级电影在线观看| 在线观看免费视频日本深夜| 中文资源天堂在线| 国产综合懂色| 日韩欧美在线二视频| 亚洲无线观看免费| 人妻夜夜爽99麻豆av| 久久精品国产清高在天天线| 狂野欧美激情性xxxx在线观看| 男女视频在线观看网站免费| 少妇人妻精品综合一区二区 | 99久久九九国产精品国产免费| 久久热精品热| 黄片wwwwww| 乱人视频在线观看| 一进一出抽搐动态| 在线播放国产精品三级| 国产精品人妻久久久久久| 欧美区成人在线视频| 精品久久久久久,| 欧美日韩中文字幕国产精品一区二区三区| 欧美成人a在线观看| 91麻豆av在线| 欧美日韩国产亚洲二区| 国产午夜福利久久久久久| 干丝袜人妻中文字幕| 校园春色视频在线观看| 亚洲午夜理论影院| 国产精品久久久久久久电影| 欧美又色又爽又黄视频| 国产成年人精品一区二区| 女人十人毛片免费观看3o分钟| 国产成人av教育| 一本精品99久久精品77| 成人高潮视频无遮挡免费网站| 自拍偷自拍亚洲精品老妇| 国产一区二区在线av高清观看| 男女视频在线观看网站免费| 欧美3d第一页| a级一级毛片免费在线观看| 亚洲电影在线观看av| 黄色女人牲交| 久久亚洲真实| 国产v大片淫在线免费观看| 搡老熟女国产l中国老女人| 欧美三级亚洲精品| 欧美bdsm另类| 69人妻影院| 日本黄色片子视频| 国产亚洲精品综合一区在线观看| 身体一侧抽搐| 亚州av有码| 日本一二三区视频观看| 久久久国产成人精品二区| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 他把我摸到了高潮在线观看| 一区二区三区四区激情视频 | 久久精品国产自在天天线| 我要搜黄色片| 别揉我奶头~嗯~啊~动态视频| 老女人水多毛片| 婷婷亚洲欧美| 人人妻人人澡欧美一区二区| 哪里可以看免费的av片| 婷婷六月久久综合丁香| 免费人成视频x8x8入口观看| 亚洲av免费高清在线观看| 免费看av在线观看网站| 大又大粗又爽又黄少妇毛片口| 国产精品国产高清国产av| 国产高潮美女av| 久久久久久九九精品二区国产| 国产蜜桃级精品一区二区三区| 亚洲av日韩精品久久久久久密| 久久久久久久久久黄片| 亚洲av日韩精品久久久久久密| 日日干狠狠操夜夜爽| 桃红色精品国产亚洲av| 国内精品久久久久久久电影| 亚洲精品456在线播放app | 国产在线男女| 伦精品一区二区三区| 成人特级黄色片久久久久久久| 日韩欧美精品v在线| 免费大片18禁| 日韩高清综合在线| 国产欧美日韩一区二区精品| 免费黄网站久久成人精品| 国产精品野战在线观看| 熟女电影av网| 91麻豆精品激情在线观看国产| 国产精品精品国产色婷婷| 亚洲精品久久国产高清桃花| 女的被弄到高潮叫床怎么办 | 日韩精品有码人妻一区| 久久久久久国产a免费观看| 桃红色精品国产亚洲av| 国产在视频线在精品| 赤兔流量卡办理| 国产精品国产高清国产av| 美女大奶头视频| 亚洲欧美日韩卡通动漫| 成年免费大片在线观看| 国产av麻豆久久久久久久| 黄色配什么色好看| 黄片wwwwww| 舔av片在线| а√天堂www在线а√下载| 国产单亲对白刺激| 国产 一区 欧美 日韩| 欧美区成人在线视频| 日韩欧美国产一区二区入口| 国产男人的电影天堂91| 1024手机看黄色片| 干丝袜人妻中文字幕| 欧美日本视频| 国产免费一级a男人的天堂| 91av网一区二区| 麻豆国产97在线/欧美| 日韩高清综合在线| 免费观看的影片在线观看| 亚洲国产精品久久男人天堂| 国产熟女欧美一区二区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av免费在线观看| 午夜精品在线福利| 精品久久久噜噜| 在线观看免费视频日本深夜| 午夜免费成人在线视频| 欧美3d第一页| 国产探花在线观看一区二区| 欧美三级亚洲精品| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人一区二区免费高清观看| 99久久成人亚洲精品观看| 欧美最新免费一区二区三区| 日本黄色片子视频| 国产精品一区二区免费欧美| 国产午夜精品论理片| 成人二区视频| 国产成人影院久久av| 亚洲精品456在线播放app | 国产乱人视频| 久久久久久伊人网av| 亚洲国产精品成人综合色| 亚洲一区高清亚洲精品| 午夜视频国产福利| 亚洲国产高清在线一区二区三| 男女做爰动态图高潮gif福利片| 综合色av麻豆| 欧美色欧美亚洲另类二区| 天堂av国产一区二区熟女人妻| 我的老师免费观看完整版| 有码 亚洲区| 男人的好看免费观看在线视频| 色av中文字幕| 99久久九九国产精品国产免费| 国产精品无大码| 一区二区三区四区激情视频 | 国产在线精品亚洲第一网站| .国产精品久久| 国产黄a三级三级三级人| 久久久久九九精品影院| 亚洲不卡免费看| 白带黄色成豆腐渣| 欧美+亚洲+日韩+国产| 琪琪午夜伦伦电影理论片6080| 熟女人妻精品中文字幕| 日本在线视频免费播放| 看十八女毛片水多多多| 国产精品,欧美在线| 毛片女人毛片| 精品久久久久久久人妻蜜臀av| 在线a可以看的网站| 国内揄拍国产精品人妻在线| 中国美女看黄片| 欧美高清性xxxxhd video| 欧美不卡视频在线免费观看| 91麻豆精品激情在线观看国产| 久久亚洲真实| 内地一区二区视频在线| 在线观看免费视频日本深夜| 国产精品av视频在线免费观看| 亚洲国产精品成人综合色| 国产又黄又爽又无遮挡在线| 亚洲中文字幕一区二区三区有码在线看| 美女高潮喷水抽搐中文字幕| 一本久久中文字幕| 国产精品一区二区免费欧美| 我的老师免费观看完整版| 丰满的人妻完整版| 一区二区三区四区激情视频 | 日韩欧美国产在线观看| 国产精品亚洲美女久久久| 免费看a级黄色片| 岛国在线免费视频观看| 极品教师在线视频| 99久久精品一区二区三区| 99久久久亚洲精品蜜臀av| 五月玫瑰六月丁香| 非洲黑人性xxxx精品又粗又长| 国产精品不卡视频一区二区| 成人毛片a级毛片在线播放| 久久精品国产自在天天线| 少妇高潮的动态图| 一级av片app| 偷拍熟女少妇极品色| 免费在线观看影片大全网站| 国产精品人妻久久久久久| 欧美成人免费av一区二区三区| 日韩av在线大香蕉| 成人鲁丝片一二三区免费| a级一级毛片免费在线观看| 2021天堂中文幕一二区在线观| 欧美激情在线99| 别揉我奶头 嗯啊视频| 在线天堂最新版资源| 精品久久国产蜜桃| 99热这里只有是精品在线观看| 国产综合懂色| 人人妻,人人澡人人爽秒播| 一夜夜www| 久久久久久久亚洲中文字幕| 日本黄大片高清| 亚洲自偷自拍三级| 国模一区二区三区四区视频| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 午夜老司机福利剧场| 亚洲不卡免费看| 22中文网久久字幕| 91久久精品国产一区二区三区| 91麻豆av在线| 日本-黄色视频高清免费观看| 久久久久久国产a免费观看| 亚洲人成网站在线播| 全区人妻精品视频| 亚洲精品国产成人久久av| 久99久视频精品免费| 国产精品不卡视频一区二区| 99热网站在线观看| 国产精品一区www在线观看 | 少妇人妻一区二区三区视频| 国产亚洲精品久久久com| 国产在线精品亚洲第一网站| 欧美三级亚洲精品| 一区福利在线观看| 国产亚洲精品久久久com| 午夜日韩欧美国产| 日日撸夜夜添| 亚洲精品影视一区二区三区av| 精品久久久久久久久久免费视频| 亚洲自拍偷在线| 亚洲国产精品成人综合色| 天堂网av新在线| 婷婷色综合大香蕉| 精品人妻视频免费看| 国产又黄又爽又无遮挡在线| 国产熟女欧美一区二区| 国产精品人妻久久久影院| 日韩一区二区视频免费看| av在线亚洲专区| 久久久久性生活片| 亚洲在线自拍视频| 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 国产精品久久久久久亚洲av鲁大| 一夜夜www| 日本三级黄在线观看| 亚洲一区二区三区色噜噜| 我要看日韩黄色一级片| 成熟少妇高潮喷水视频| 久久香蕉精品热| 狂野欧美激情性xxxx在线观看| 97碰自拍视频| 少妇人妻一区二区三区视频| 亚洲自拍偷在线| www.色视频.com| 久久午夜亚洲精品久久| 少妇熟女aⅴ在线视频| 噜噜噜噜噜久久久久久91| 免费看av在线观看网站| 国产成人一区二区在线| 一a级毛片在线观看| 亚洲欧美清纯卡通| 国产色婷婷99| 啦啦啦观看免费观看视频高清| 久久精品人妻少妇| videossex国产| 午夜视频国产福利| 亚洲精品一区av在线观看| 久久午夜福利片| 麻豆一二三区av精品| eeuss影院久久| 精品人妻熟女av久视频| 亚洲经典国产精华液单| 亚洲18禁久久av| 久久精品国产99精品国产亚洲性色| 色哟哟哟哟哟哟| 最近中文字幕高清免费大全6 | 人妻久久中文字幕网| 国语自产精品视频在线第100页| 露出奶头的视频| 中文字幕免费在线视频6| 三级国产精品欧美在线观看| 88av欧美| 午夜福利成人在线免费观看| 久久香蕉精品热| 国产男人的电影天堂91| 亚洲最大成人手机在线| 麻豆一二三区av精品| 成人毛片a级毛片在线播放| 长腿黑丝高跟| 一进一出抽搐动态| 淫妇啪啪啪对白视频| 日日摸夜夜添夜夜添小说| 国产一级毛片七仙女欲春2| 欧美黑人欧美精品刺激| 18禁黄网站禁片午夜丰满| 日日摸夜夜添夜夜添小说| 人人妻人人澡欧美一区二区| 99在线人妻在线中文字幕| 在线看三级毛片| 在线播放无遮挡| 亚洲成人中文字幕在线播放| 欧美在线一区亚洲| 级片在线观看| 久久精品91蜜桃| 色综合站精品国产| 日韩欧美精品免费久久| 国产av一区在线观看免费| av.在线天堂| 亚洲熟妇中文字幕五十中出| 三级毛片av免费| 午夜日韩欧美国产| 91av网一区二区| 久久6这里有精品| 国产精品一区二区免费欧美| 亚洲av二区三区四区| 亚洲内射少妇av| 99久久成人亚洲精品观看| avwww免费| 亚洲一区高清亚洲精品| 亚洲av中文av极速乱 | 日韩在线高清观看一区二区三区 | 亚洲精品国产成人久久av| 真实男女啪啪啪动态图| 美女大奶头视频| 一级黄片播放器| 国产黄片美女视频| 国产精品国产三级国产av玫瑰| 97超级碰碰碰精品色视频在线观看| av在线蜜桃| 国产亚洲精品久久久com| 欧美日韩国产亚洲二区| 91久久精品国产一区二区三区| 十八禁国产超污无遮挡网站| 午夜福利欧美成人| 99热这里只有是精品在线观看| 日本免费一区二区三区高清不卡| 午夜免费激情av| 网址你懂的国产日韩在线| 偷拍熟女少妇极品色| 尾随美女入室| 色吧在线观看| 中文在线观看免费www的网站| 一a级毛片在线观看| 成人鲁丝片一二三区免费| 亚洲av中文av极速乱 | 成人国产一区最新在线观看| 亚洲黑人精品在线| 又黄又爽又刺激的免费视频.| 国产蜜桃级精品一区二区三区| 如何舔出高潮| 免费在线观看影片大全网站| 久久久久久久精品吃奶| 久久久久久久久中文| 变态另类丝袜制服| 国内精品久久久久精免费| 男人舔女人下体高潮全视频| 毛片女人毛片| 欧美三级亚洲精品| 一区二区三区免费毛片| 老熟妇乱子伦视频在线观看| 国内毛片毛片毛片毛片毛片| 大又大粗又爽又黄少妇毛片口| 亚洲av美国av| 国产视频内射| 此物有八面人人有两片| 一本一本综合久久| 午夜视频国产福利| 99久久中文字幕三级久久日本| 亚洲av电影不卡..在线观看| 精品一区二区三区人妻视频| 成年人黄色毛片网站| 少妇熟女aⅴ在线视频| 国内毛片毛片毛片毛片毛片| 国产伦一二天堂av在线观看| 91久久精品国产一区二区三区| 亚洲中文日韩欧美视频| 久久精品影院6| 欧美最黄视频在线播放免费| 高清日韩中文字幕在线| 成年女人毛片免费观看观看9| 日本色播在线视频| 久99久视频精品免费| 亚洲中文日韩欧美视频| 久久中文看片网| 99久国产av精品| 日本五十路高清| 干丝袜人妻中文字幕| 嫩草影视91久久| 人妻少妇偷人精品九色| 亚洲 国产 在线| 男人的好看免费观看在线视频| av.在线天堂| 性色avwww在线观看| 一级黄色大片毛片| 国产精品久久视频播放| 黄色日韩在线| 亚洲欧美清纯卡通| 欧美不卡视频在线免费观看| 欧美成人a在线观看| 国产免费男女视频| 亚洲欧美日韩卡通动漫| 白带黄色成豆腐渣| 久久午夜福利片| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品综合一区二区三区| 国产亚洲精品综合一区在线观看| 少妇的逼好多水| 国产在视频线在精品| 久久亚洲精品不卡| 干丝袜人妻中文字幕| 久久亚洲精品不卡| 亚洲男人的天堂狠狠| 国产精品综合久久久久久久免费| 国内精品久久久久久久电影| 欧美日本亚洲视频在线播放|