• <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í)
    国产色婷婷99| 亚洲男人天堂网一区| 久久国产精品大桥未久av| 国产精品偷伦视频观看了| 爱豆传媒免费全集在线观看| 一级毛片我不卡| 国产一级毛片在线| 国产又色又爽无遮挡免| 久久99一区二区三区| 国产成人91sexporn| xxxhd国产人妻xxx| 卡戴珊不雅视频在线播放| 国产男女内射视频| 九色亚洲精品在线播放| 99久国产av精品国产电影| 在线 av 中文字幕| 欧美成人精品欧美一级黄| 国产在线一区二区三区精| 亚洲av欧美aⅴ国产| 国产精品无大码| 日日撸夜夜添| 色吧在线观看| 日韩中文字幕视频在线看片| 欧美精品av麻豆av| 亚洲精品久久午夜乱码| 99九九在线精品视频| 老司机影院成人| 国产日韩欧美视频二区| 热re99久久精品国产66热6| 极品少妇高潮喷水抽搐| 亚洲国产看品久久| 亚洲男人天堂网一区| 五月开心婷婷网| 狠狠婷婷综合久久久久久88av| 国产精品 国内视频| 别揉我奶头~嗯~啊~动态视频 | 侵犯人妻中文字幕一二三四区| 激情视频va一区二区三区| 成年av动漫网址| 国产亚洲av片在线观看秒播厂| 在线观看一区二区三区激情| 精品一品国产午夜福利视频| 人妻 亚洲 视频| 欧美日本中文国产一区发布| 亚洲欧美成人精品一区二区| 人妻一区二区av| 国产在线免费精品| 天堂中文最新版在线下载| 这个男人来自地球电影免费观看 | 久久久久国产精品人妻一区二区| 欧美成人午夜精品| 亚洲专区中文字幕在线 | 爱豆传媒免费全集在线观看| 亚洲精品国产色婷婷电影| 国精品久久久久久国模美| 一本久久精品| 亚洲五月色婷婷综合| 中文字幕制服av| 欧美av亚洲av综合av国产av | 欧美少妇被猛烈插入视频| 国产高清国产精品国产三级| 久久久久久久久免费视频了| 亚洲国产欧美网| 下体分泌物呈黄色| 悠悠久久av| 午夜福利视频在线观看免费| 精品一品国产午夜福利视频| 日本黄色日本黄色录像| 叶爱在线成人免费视频播放| 一级,二级,三级黄色视频| 波多野结衣一区麻豆| 国产精品人妻久久久影院| 丰满迷人的少妇在线观看| 日日爽夜夜爽网站| 亚洲自偷自拍图片 自拍| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠躁躁| 亚洲欧美色中文字幕在线| 国产精品久久久人人做人人爽| 日韩,欧美,国产一区二区三区| 亚洲七黄色美女视频| 91国产中文字幕| 七月丁香在线播放| 国产高清国产精品国产三级| 99久久99久久久精品蜜桃| 国产精品99久久99久久久不卡 | 免费观看a级毛片全部| 最新在线观看一区二区三区 | 黄片无遮挡物在线观看| 成人国产麻豆网| 国产亚洲欧美精品永久| 黄色怎么调成土黄色| 中文字幕精品免费在线观看视频| 欧美成人午夜精品| 男女午夜视频在线观看| 考比视频在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲成色77777| 最新的欧美精品一区二区| 最近最新中文字幕免费大全7| 亚洲av电影在线观看一区二区三区| 操美女的视频在线观看| 国产熟女欧美一区二区| 在线观看免费日韩欧美大片| 精品久久蜜臀av无| 免费久久久久久久精品成人欧美视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品国产av蜜桃| 成年人免费黄色播放视频| av网站免费在线观看视频| 日韩中文字幕视频在线看片| 国产在视频线精品| 日韩一区二区三区影片| 日韩欧美一区视频在线观看| 免费观看人在逋| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕色久视频| 国产在视频线精品| 中文字幕亚洲精品专区| 午夜福利视频在线观看免费| 欧美日本中文国产一区发布| 大话2 男鬼变身卡| 亚洲 欧美一区二区三区| 天天操日日干夜夜撸| 丝袜人妻中文字幕| 性少妇av在线| 精品一区在线观看国产| 午夜福利视频精品| 精品国产超薄肉色丝袜足j| 自线自在国产av| 热99久久久久精品小说推荐| 午夜福利在线免费观看网站| 亚洲av中文av极速乱| 色婷婷av一区二区三区视频| 久久午夜综合久久蜜桃| 无限看片的www在线观看| 免费看av在线观看网站| 亚洲国产毛片av蜜桃av| 久久鲁丝午夜福利片| 麻豆av在线久日| 美女主播在线视频| www.精华液| netflix在线观看网站| 韩国av在线不卡| av线在线观看网站| 亚洲精华国产精华液的使用体验| 国产伦人伦偷精品视频| 麻豆av在线久日| 日韩不卡一区二区三区视频在线| 久久影院123| 午夜91福利影院| 麻豆乱淫一区二区| 极品少妇高潮喷水抽搐| 91精品三级在线观看| 久久久亚洲精品成人影院| 少妇人妻久久综合中文| 人人澡人人妻人| 九色亚洲精品在线播放| av在线app专区| 人人妻人人爽人人添夜夜欢视频| 纵有疾风起免费观看全集完整版| 老鸭窝网址在线观看| 久久亚洲国产成人精品v| 欧美激情高清一区二区三区 | av网站在线播放免费| 日韩免费高清中文字幕av| 999精品在线视频| 尾随美女入室| 电影成人av| 亚洲国产av影院在线观看| 久久综合国产亚洲精品| 中文乱码字字幕精品一区二区三区| 97在线人人人人妻| 成人免费观看视频高清| 国产精品一国产av| 国产一区二区 视频在线| 国产精品久久久久久精品电影小说| 97在线人人人人妻| 色94色欧美一区二区| 观看av在线不卡| av天堂久久9| 免费在线观看视频国产中文字幕亚洲 | 蜜桃国产av成人99| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 亚洲国产av影院在线观看| 精品少妇久久久久久888优播| 亚洲精品成人av观看孕妇| 亚洲国产毛片av蜜桃av| 亚洲成色77777| 日本午夜av视频| 国产人伦9x9x在线观看| 老司机影院成人| 无遮挡黄片免费观看| 免费观看av网站的网址| 精品国产一区二区久久| av国产精品久久久久影院| 韩国高清视频一区二区三区| 久久亚洲国产成人精品v| 麻豆av在线久日| 日韩av免费高清视频| 欧美人与性动交α欧美软件| 亚洲av在线观看美女高潮| 黄色视频不卡| 亚洲人成电影观看| 国产精品三级大全| 老司机亚洲免费影院| 日日啪夜夜爽| 五月开心婷婷网| 少妇被粗大猛烈的视频| 午夜免费鲁丝| 青春草国产在线视频| 天天影视国产精品| 黄频高清免费视频| 波多野结衣一区麻豆| 啦啦啦在线免费观看视频4| 国产一级毛片在线| 黄色视频不卡| 日本色播在线视频| 在线观看国产h片| 国产乱来视频区| 亚洲欧美成人精品一区二区| 视频在线观看一区二区三区| 一级毛片电影观看| 国产日韩欧美视频二区| 亚洲人成网站在线观看播放| 久久久亚洲精品成人影院| 免费av中文字幕在线| 国产精品免费视频内射| 久热这里只有精品99| 欧美精品人与动牲交sv欧美| 啦啦啦在线观看免费高清www| 曰老女人黄片| 久热这里只有精品99| 日韩一本色道免费dvd| 亚洲av男天堂| 欧美日韩国产mv在线观看视频| 亚洲国产最新在线播放| 乱人伦中国视频| 国产av码专区亚洲av| 色吧在线观看| 国产精品久久久久成人av| 极品少妇高潮喷水抽搐| 国产福利在线免费观看视频| www.精华液| 老司机在亚洲福利影院| 考比视频在线观看| 国产精品成人在线| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 久久精品人人爽人人爽视色| 国产极品粉嫩免费观看在线| 91精品国产国语对白视频| 亚洲成人免费av在线播放| 人人妻人人爽人人添夜夜欢视频| 别揉我奶头~嗯~啊~动态视频 | 最近最新中文字幕大全免费视频 | 国产激情久久老熟女| 国产日韩欧美视频二区| 国产成人一区二区在线| 国产熟女午夜一区二区三区| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 五月开心婷婷网| 国产精品秋霞免费鲁丝片| 美国免费a级毛片| 青春草视频在线免费观看| 国产成人精品无人区| 国产一区二区三区av在线| 日本色播在线视频| 欧美久久黑人一区二区| 午夜日韩欧美国产| 蜜桃国产av成人99| 久久久久久久久久久免费av| 日韩精品有码人妻一区| 午夜福利视频精品| 国产亚洲精品第一综合不卡| 免费女性裸体啪啪无遮挡网站| 操美女的视频在线观看| 一区二区三区乱码不卡18| 亚洲少妇的诱惑av| 国产福利在线免费观看视频| 亚洲国产中文字幕在线视频| 97精品久久久久久久久久精品| 日韩欧美一区视频在线观看| 亚洲一区二区三区欧美精品| 搡老乐熟女国产| 欧美日韩亚洲高清精品| 久久久国产精品麻豆| 国产成人欧美在线观看 | 色婷婷av一区二区三区视频| 日韩免费高清中文字幕av| 秋霞伦理黄片| 国产精品偷伦视频观看了| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 香蕉丝袜av| av网站在线播放免费| 黄网站色视频无遮挡免费观看| 久久精品国产亚洲av涩爱| 日韩不卡一区二区三区视频在线| netflix在线观看网站| 日韩电影二区| 男女免费视频国产| 一边亲一边摸免费视频| 韩国av在线不卡| 国产在线视频一区二区| av福利片在线| 欧美日韩国产mv在线观看视频| 美女午夜性视频免费| 最新的欧美精品一区二区| 国产99久久九九免费精品| av在线app专区| 欧美精品一区二区大全| 国产亚洲欧美精品永久| 一级a爱视频在线免费观看| 精品人妻在线不人妻| 香蕉国产在线看| 色播在线永久视频| 日本欧美视频一区| 欧美日韩成人在线一区二区| 美女福利国产在线| 男女午夜视频在线观看| 国产成人精品久久二区二区91 | 国产精品一区二区精品视频观看| 两个人免费观看高清视频| 两性夫妻黄色片| 色播在线永久视频| 久热这里只有精品99| 婷婷色综合大香蕉| 国产免费一区二区三区四区乱码| 一级爰片在线观看| 视频在线观看一区二区三区| 久久99热这里只频精品6学生| 日本91视频免费播放| 人成视频在线观看免费观看| 丰满迷人的少妇在线观看| 日韩制服丝袜自拍偷拍| 男女下面插进去视频免费观看| 青青草视频在线视频观看| 黑人巨大精品欧美一区二区蜜桃| 精品免费久久久久久久清纯 | 日韩大码丰满熟妇| 久久精品久久精品一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产女主播在线喷水免费视频网站| 久久人人爽av亚洲精品天堂| 97在线人人人人妻| 国产成人精品久久久久久| 国产成人啪精品午夜网站| 18禁动态无遮挡网站| 黄频高清免费视频| 久久午夜综合久久蜜桃| 丝瓜视频免费看黄片| 纵有疾风起免费观看全集完整版| 国产亚洲最大av| 新久久久久国产一级毛片| 男女边摸边吃奶| 人人澡人人妻人| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 丰满少妇做爰视频| 国产不卡av网站在线观看| 国产欧美日韩综合在线一区二区| 亚洲精品国产av成人精品| 亚洲欧美成人精品一区二区| 亚洲国产av新网站| 欧美在线一区亚洲| 在线观看国产h片| 精品亚洲成a人片在线观看| a级毛片在线看网站| 亚洲精品av麻豆狂野| 成人午夜精彩视频在线观看| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 欧美黑人欧美精品刺激| 中文字幕色久视频| 久久精品久久久久久久性| 肉色欧美久久久久久久蜜桃| 亚洲国产精品一区三区| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 咕卡用的链子| 亚洲成人国产一区在线观看 | 精品一区二区免费观看| 老司机影院成人| 国产一区二区三区综合在线观看| av又黄又爽大尺度在线免费看| 一级黄片播放器| 欧美日韩亚洲综合一区二区三区_| 亚洲欧洲精品一区二区精品久久久 | 男人爽女人下面视频在线观看| 国产在线视频一区二区| 亚洲综合色网址| 一区二区三区激情视频| 国产有黄有色有爽视频| 亚洲情色 制服丝袜| 免费少妇av软件| 母亲3免费完整高清在线观看| 高清不卡的av网站| 亚洲精品久久成人aⅴ小说| 99香蕉大伊视频| 老司机在亚洲福利影院| 女人久久www免费人成看片| 一边亲一边摸免费视频| 悠悠久久av| 欧美日韩亚洲综合一区二区三区_| 欧美少妇被猛烈插入视频| 精品亚洲乱码少妇综合久久| 男女免费视频国产| 看免费av毛片| 成年动漫av网址| 日日撸夜夜添| 777久久人妻少妇嫩草av网站| 老汉色av国产亚洲站长工具| www.精华液| 国产亚洲精品第一综合不卡| 十八禁人妻一区二区| netflix在线观看网站| 国产探花极品一区二区| 午夜久久久在线观看| 伊人久久国产一区二区| 成人黄色视频免费在线看| av女优亚洲男人天堂| 国产成人精品在线电影| 久久久久久久久免费视频了| 天天躁狠狠躁夜夜躁狠狠躁| 日本av免费视频播放| 色婷婷久久久亚洲欧美| 美女高潮到喷水免费观看| svipshipincom国产片| 日本91视频免费播放| 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 大香蕉久久网| 搡老乐熟女国产| 巨乳人妻的诱惑在线观看| 国产国语露脸激情在线看| 国产乱来视频区| 日韩人妻精品一区2区三区| 亚洲自偷自拍图片 自拍| 午夜福利在线免费观看网站| 久久综合国产亚洲精品| 亚洲第一av免费看| 少妇被粗大猛烈的视频| 久久人人爽人人片av| 丝袜美腿诱惑在线| 久久99热这里只频精品6学生| 国产又色又爽无遮挡免| 爱豆传媒免费全集在线观看| 嫩草影视91久久| 波多野结衣一区麻豆| 最近手机中文字幕大全| 丝袜在线中文字幕| 日韩人妻精品一区2区三区| 丰满少妇做爰视频| 久久精品国产亚洲av涩爱| 中国三级夫妇交换| 丝袜在线中文字幕| 国产精品 国内视频| 国产精品嫩草影院av在线观看| av有码第一页| 亚洲国产精品一区三区| 一区二区日韩欧美中文字幕| 狂野欧美激情性xxxx| 老汉色∧v一级毛片| 1024视频免费在线观看| 成年人午夜在线观看视频| 你懂的网址亚洲精品在线观看| 久久国产亚洲av麻豆专区| 亚洲男人天堂网一区| 亚洲精品国产色婷婷电影| 日韩欧美一区视频在线观看| 一本—道久久a久久精品蜜桃钙片| 狠狠精品人妻久久久久久综合| 制服诱惑二区| 一区二区三区激情视频| 婷婷色麻豆天堂久久| 九色亚洲精品在线播放| 在线看a的网站| 亚洲综合精品二区| 日韩熟女老妇一区二区性免费视频| 精品少妇久久久久久888优播| a级片在线免费高清观看视频| 高清av免费在线| 巨乳人妻的诱惑在线观看| 一级毛片 在线播放| 中文字幕亚洲精品专区| 男人操女人黄网站| 在线观看www视频免费| 亚洲av国产av综合av卡| 国产 精品1| 最近最新中文字幕大全免费视频 | 午夜91福利影院| 成人亚洲精品一区在线观看| 亚洲色图综合在线观看| 中文字幕人妻丝袜一区二区 | 免费看av在线观看网站| 1024视频免费在线观看| 亚洲精品日本国产第一区| 波野结衣二区三区在线| 午夜福利视频精品| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 老汉色av国产亚洲站长工具| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 国产人伦9x9x在线观看| 午夜福利乱码中文字幕| 又粗又硬又长又爽又黄的视频| 波多野结衣av一区二区av| 亚洲av男天堂| 国产在线一区二区三区精| 欧美人与性动交α欧美精品济南到| 少妇 在线观看| 天天操日日干夜夜撸| 五月天丁香电影| 亚洲av日韩在线播放| 日日啪夜夜爽| 免费观看a级毛片全部| 亚洲三区欧美一区| 久久久亚洲精品成人影院| 老司机影院成人| 午夜福利网站1000一区二区三区| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 大片免费播放器 马上看| av有码第一页| 少妇猛男粗大的猛烈进出视频| 精品少妇黑人巨大在线播放| 中文字幕最新亚洲高清| 国产又爽黄色视频| 国产精品亚洲av一区麻豆 | 精品国产乱码久久久久久小说| 啦啦啦啦在线视频资源| 一区二区三区精品91| 久久人人爽人人片av| 一级,二级,三级黄色视频| 97人妻天天添夜夜摸| a级毛片在线看网站| 免费在线观看视频国产中文字幕亚洲 | 精品一区二区免费观看| 99久久人妻综合| 99热网站在线观看| 伊人久久国产一区二区| 久久久精品国产亚洲av高清涩受| av线在线观看网站| 日本vs欧美在线观看视频| 国产日韩欧美亚洲二区| 亚洲av成人精品一二三区| 成年av动漫网址| 午夜老司机福利片| 久久人人97超碰香蕉20202| 日本wwww免费看| 精品亚洲乱码少妇综合久久| 日韩一卡2卡3卡4卡2021年| 亚洲国产精品国产精品| 性少妇av在线| 一个人免费看片子| 欧美日韩亚洲国产一区二区在线观看 | av有码第一页| 91精品伊人久久大香线蕉| 性高湖久久久久久久久免费观看| 嫩草影院入口| 免费观看人在逋| 国产片特级美女逼逼视频| 国产免费福利视频在线观看| 99久国产av精品国产电影| 精品福利永久在线观看| 最近的中文字幕免费完整| 婷婷色麻豆天堂久久| 热re99久久精品国产66热6| 日本午夜av视频| 黑人巨大精品欧美一区二区蜜桃| av片东京热男人的天堂| 在线观看人妻少妇| 亚洲综合色网址| 精品视频人人做人人爽| 久久久精品94久久精品| 亚洲成人国产一区在线观看 | 国产1区2区3区精品| av有码第一页| 精品久久久精品久久久| 国产免费福利视频在线观看| 丝袜脚勾引网站| 日韩大码丰满熟妇| 悠悠久久av| 99久久精品国产亚洲精品| 欧美av亚洲av综合av国产av | 亚洲国产欧美日韩在线播放| 精品一区二区三区四区五区乱码 | 一二三四在线观看免费中文在| 美女大奶头黄色视频| 嫩草影院入口| 性高湖久久久久久久久免费观看| 美国免费a级毛片| 男女免费视频国产| 丝袜脚勾引网站| 又大又爽又粗| 亚洲国产欧美日韩在线播放| 亚洲国产av新网站| 亚洲欧美色中文字幕在线| 亚洲国产中文字幕在线视频| 国产精品国产av在线观看| 天堂8中文在线网| 国产一区二区三区综合在线观看| 午夜91福利影院| 成年美女黄网站色视频大全免费| 黄色视频在线播放观看不卡| 大香蕉久久成人网| 999精品在线视频| 一区二区av电影网| 美女视频免费永久观看网站| 大香蕉久久成人网| 国产片特级美女逼逼视频|