氣候變化對南極冰面湖的影響研究：以埃默里和拉森A冰架為例

王輝 盧善龍 丁俊 邱玉寶 唐海龍,3 閆強

(1 湘潭大學,土木工程與力學學院,湖南 湘潭 411105;2 中國科學院遙感與數字地球研究所,中國科學院數字地球重點實驗室,北京 100101;3 成都理工大學地球科學學院,四川 成都 610059;4 中國地質科學院礦產資源研究所,北京 100037)

提要 冰面湖作為區域氣候變化的靈敏指示器,對研究全球氣候變化背景下南極冰凍圈的穩定性具有重要意義。本文基于2000—2019年區域氣象數據、Landsat ETM+/OLI 影像與LIMA(Landsat Image Mosaic of Antarctica)鑲嵌影像,選取東南極的埃默里冰架和西南極的拉森A冰架作為研究區,采用歸一化水體指數(Normalized Difference Water Index,NDWI)提取冰面湖面積信息,對冰面湖面積與氣象數據開展遙相關分析,探索冰面湖與區域氣候變化的聯系。結果表明:埃默里冰面湖總面積較為穩定并有減小趨勢,拉森A冰面湖總面積則呈增加趨勢;冰面湖與氣溫和海溫之間存在正向反饋特征;兩個冰架區域的反饋特征規律明顯不同,埃默里冰面湖面積受氣溫影響顯著,但受海溫影響并不明顯,而拉森A冰面湖則對海溫變化敏感,其面積與海溫變化趨勢一致。

0 引言

聯合國氣候變化專門委員會(Intergovernmental Panel on Climate Change,IPCC)《第五次氣候變化評估報告》指出,以全球氣候變暖為主要特征的氣候變化已成事實,在1880—2012年間,全球陸地平均氣溫上升了0.65～1.06℃,2016—2035年預計將上升0.3～0.7℃,2081—2100年將上升0.3～0.7℃[1]。在氣候變化模式RCP 8.5(Representative Concentration Pathway)情景下,未來100年全球平均氣溫將升高3.7℃[2]。預計2300年南極冰蓋融化對海平面上升的貢獻為3～10 m[3-4]。Bronselaer 等[5]認為,在RCP 8.5 情景下,若加入冰蓋融化的影響,超過全球大氣變暖平均值(1.5～2.0℃)的時間將往后推遲10 a以上。據世界冰川監測服務處(World Glacier Monitoring Service,WGMS)公布的數據顯示,伴隨全球氣溫上升,從2000年開始,全球所有冰川冰蝕率平均值已達0.5 m·a–1,與20世紀80年代相比提升了2倍。美國宇航局下屬噴氣推進實驗室對1996—2005年格陵蘭島進行了研究,結果指出冰川融化速度逐年加快,在2005年一年中格陵蘭島的冰川融化量是10年前的兩倍[6]。世界自然基金會(World Wide Fund for Nature,WWF)的一份報告顯示,隨著氣候變暖不斷加速,青藏高原的變暖速率也在不斷加快,自20世紀50年代以來,年平均氣溫有顯著上升趨勢,在1961—2007年,年平均氣溫的上升速率為0.37℃·(10a)–1,而且冷季溫度上升加速度要大于暖季[7],1990—2010年,青藏高原冰湖數目增加了1 099個,面積從553.9 km2增加到682.4 km2[8],這與氣候變化趨勢相符。

南極獨特的自然條件使其成為對全球氣候變化最為敏感的指示器,科學研究發現從20世紀50年代到20世紀末這40 多年里,南極地面平均溫度上升了2.5℃,這與全球氣候變暖同步[9-10]。冰架是南極冰蓋邊緣最為顯著的物質交換區,同時也是南極冰蓋邊緣最為活躍的特征區,它們的動態變化對南極乃至全世界都有著深刻影響。冰架表面融水流經冰蓋會形成冰面湖,而冰面湖可能引發冰架的崩解[11],這使得冰面湖成為冰蓋響應氣候變化的一個日益重要的因素[12]。拉森B 冰架(Larsen B Ice Shelf)崩解后其所對應的補給冰川流速明顯加快,該冰架表面的冰面湖已經成為儲存融水破壞性的象征[13]。在氣候變暖情況下,隨之增加的冰面湖加速了南極未來冰量的損失[14]。隨著全球范圍的氣候變暖,越來越多的冰面湖將在南極冰架表面形成[15],南極冰面湖的關注度也在不斷提高。Bell 等[16]研究了冰面湖與南極冰蓋質量平衡間的關系,DeConto和Pollard[17]的研究模型顯示南極冰面湖有很大可能會造成冰架邊緣的坍塌,Stokes 等[18]利用2017年的遙感影像分析得到東南極冰面湖的大致分布,Langley 等[19]對南極東部朗霍夫德冰川冰面湖進行了研究,結果顯示冰面湖是東南極冰川動力學的重要部分,LaBarbera和MacAyeal[20]則通過構建模型探討了西南極喬治六世冰架與接地線運動之間的聯系。總的來說,南極冰面湖的研究多集中在冰川動力學和冰架崩解等方面,對冰面湖關于氣候變化的響應關注較少。

東南極最大的埃默里冰架(Amery Ice Shelf)是區域內重要的冰流入海口,也是與南大洋主要的熱量與物質交換區。南大洋氣候變化會對埃默里冰架產生重大影響,同時,埃默里冰架的融冰狀態也可能會潛在影響海平面上升[21],而冰面湖作為冰蓋反饋氣候變化的重要組成部分,會影響冰蓋上的冰融化和冰流,并可能導致冰架解體。拉森A冰架(Larsen A Ice Shelf)位于西南極的南極半島上,是西南極冰蓋的重要組成部分。目前南極冰蓋融化的主要因素依然存在爭議,本文以埃默里冰架和拉森A冰架為例,研究南極冰面湖的變化特征及其影響因素,這對于理解南極冰川與區域氣候以及南極冰川與海洋環境間的相互作用關系具有重要意義。

1 研究區域與數據

埃默里冰架(63°E～74°E,68°S～73°S)位于東南極北查爾斯王子山與拉斯曼丘陵之間,地勢較高,平均海拔約2 500 m。冰架長約500 km,寬在50～100 km 之間,面積約 71 260 km2,厚度在300～2 500 m 之間[22],與蘭伯特冰川(Lambert Glacier)組成的蘭伯特冰川-埃默里冰架系統(Lambert Glacier-Amery Ice Shelf System)是東南極地區最大的冰流冰架系統。埃默里冰架外延長度只占整個南極洲海岸線的1.7%,但冰流量卻占東南極的14%[23]。由于夏季埃默里冰架冰面湖集中于冰架西部,故本文冰面湖數據主要取自該冰架西部。拉森A冰架(59°W～61°W,64°S～65°S)位于威德爾海西北部(Northwest Weddell Sea),接近于南極半島的頂端,除極少數裸露的巖石之外,幾乎完全被永久性積雪和冰川所覆蓋。拉森A冰架于1995年1月發生崩解,研究區為拉森A冰架崩解后的區域范圍,面積約700 km2(圖1)。

圖1 研究區分布圖.底圖為南極洲陸地衛星鑲嵌圖LIMA,長城站位于拉森A 西北部南設得蘭群島喬治王島上,中山站位于埃默里冰架拉斯曼丘陵上(圖中紅點位置),放大圖中的紅色方框部分為精度驗證樣本分布的位置Fig.1.Location of study area.The base image is LIMA,a Landsat Image Mosaic of Antarctica.Great Wall Station is located on King George Island in the South Shetland Islands northwest of Larsen A,and Zhongshan Station is located on the Larsemann Hills of the Amery Ice Shelf (red dot).The red box in the enlarged figure shows the location of the sample distribution for accuracy verification

研究數據包括衛星影像、氣溫以及海溫。影像數據用于提取冰面湖面積信息,影像包括南極洲陸地衛星鑲嵌圖LIMA(Landsat Image Mosaic of Antarctica)[24]和Landsat ETM+(Enhanced Thematic Mapper)/OLI(Operational Lander Imager) 15 m分辨率影像。LIMA 影像由超過1 000 景成像于1999—2003年間的陸地衛星ETM+影像鑲嵌而成,這些影像在南極洲上空拼接成一幅幾乎無云的影像。LIMA 有真彩色組合與假彩色組合兩種,空間分辨率均為15 m。本文因使用歸一化水體指數(Normalized Difference Water Index,NDWI)方法提取水體,故下載的影像為假彩色組合。冰面徑流特征通常比其周圍冰的顏色更深更藍,利用LIMA 影像可明確識別冰面湖徑流。Landsat為美國國家航天局(NASA)的陸地衛星,Landsat-7 衛星上的Landsat 增強型專題制圖儀(ETM+)傳感器以16天為重復周期自1999年7月以來幾乎連續獲得了全球的影像。Landsat-7 ETM+影像由8個光譜段構成,1～7波段的空間分辨率為30 m,全色波段8的分辨率為15 m。所有的頻段均可以采集兩種增益設置(高或低)中的一種,以達到提高輻射敏感度和動態范圍的目的,其中第6波段可以采集所有場景的高低增益。每景影像的圖幅近似大小為南北長170 km,東西長183 km。Landsat-8在 2013年 2月 11日發射升空,攜帶有OLI(Operational Land Imager,陸地成像儀)和TIRS(Thermal Infrared Sensor,熱紅外傳感器)兩個主要載荷。其中OLI 由科羅拉多州的鮑爾航天技術公司研制,TIRS 由NASA的戈達德太空飛行中心研制。

冰架表面融化多發生在南半球的夏季,為避免云、極夜等不利因素的影響,數據獲取時間為每年的1—2月。另因2000—2003年、2011年和2013年Landsat 受云量影響較大,缺少有效的影像數據,導致這6年的數據缺失,所以本文搜集了埃默里冰架2004—2010年、2012年Landsat-7 15 m 分辨率的數據影像,以及2014—2019年Landsat-8 15 m 分辨率的數據影像。受云量等因素影響,本文只搜集到拉森A冰架2003—2008年、2012年和2016年的影像,其中2003年數據是利用LIMA 影像提取獲得,2005—2008年和2012年為Landsat-7 影像,2016年為Landsat-8 影像,為方便研究,本文將 2003年的冰面湖面積看做1999—2003年間冰面湖面積的平均值。

一年中若出現多景質量較好的影像,選擇冰面湖面積最大的影像為研究數據。Landsat 系列影像數據通過美國地質調查局(United States Geological Survey,USGS)的官方網站(http://earthexplorer.usgs.gov/)獲得,LIMA 影像通過網站(https://lima.usgs.gov/)獲得。具體如表1所示。

表1 遙感影像數據Table1.Remote sensing image data

氣溫與海溫在氣候變化中扮演著重要角色,在當前全球變暖的趨勢下,氣溫的升高正導致南北極冰川融化加劇[25],而海溫的變化則可反映海洋對大氣和陸地的影響[26]。本研究選用氣溫與海溫兩項指標來分析氣候變化背景下冰面湖的變化與大氣系統和海洋系統之間的關系。其中,氣溫數據選用2000—2019年中山站(海拔18 m)與長城站(海拔10 m)的地面氣象站氣溫數據,用于區域氣候變化分析。中山站的氣象數據最初是人工觀測,在2002年之后有線遙測自動觀測系統被用于氣溫等氣象數據的觀測。其中,月平均數據由日平均數據計算得到,而日平均數據采用國際標準時,即通過每天的00:00、06:00、12:00 以及18:00的觀測數據得到。長城站自建立之后,其數據就被一些學者作為研究南極半島氣候變化的重要參考[27-28]。文中所用氣象站數據通過由南極科學委員會資助的READER數據庫獲得[29]。海溫數據選用美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)地球系統研究實驗室(Earth System Research Laboratory)物理科學部(Physical Sciences Division)的COBE-SST2 數據集(2000—2017年兩個冰架區域的海表面溫度數據)。COBE-SST2 數據源為日本海洋資料中心(Japanese Oceanographic Data Center)的監控數據,空間分辨率約100 km。該數據集通過估算全球平均海表溫度長時間序列的偏差,對海表面溫度進行了修正[30-31]。數據使用前根據研究區范圍進行了裁剪和均值處理。

2 研究方法

2.1 冰面湖邊界提取

在每年的1—2月的南極冰蓋消融期,冰架表面分布著大量的消融中的冰、雪,以及被水體滲透而產生的藍冰。在近紅外波段,水體有極低的反射率,而在藍綠波段水體卻有極高的反射率,所以本文使用歸一化水體提取指數NDWI 提取水體,公式如下:

式中,Band1為綠色波段,Band2為近紅外波段,兩者在經過全色波段與多光譜波段的融合后,空間分辨率由原先的30 m 提高到了15 m。閾值的選擇通過參考相關文獻資料和提取效果最終確定[18],每年的具體閾值見表1。

因研究區域沒有更高分辨率衛星遙感影像可供使用,本文采用目視對比方法對遙感提取的水體結果進行了對比。樣本以分布在冰架不同區域為原則,分別從兩個冰架隨機抽取了不同形狀冰面湖12 組。直觀目視比較結果表明,利用水體指數提取的冰面湖冰水邊界與人工目視解譯的結果具有較好的一致性(圖2)。面積統計對比結果表明,遙感提取的冰面湖面積的總體精度為0.948(表2),其中:

目視解譯是指憑借人的眼睛或光學儀器,依靠解譯者的經驗、知識和相關資料,通過分析、推理、判斷,提取遙感圖像中有用的信息的遙感圖像解譯手段,是地學研究和遙感應用的一項基本技能,該方法優點在于能夠獲得較高精度的數據,缺點是操作繁瑣,耗時長。本文中總體精度是指12 組樣本重合區域之和與目視解譯結果之和的比值。

圖2 NDWI 結果與目視解譯結果對比Fig.2.NDWI results were compared with visual interpretation results

表2 NDWI 結果與目視解譯面積對比Table2.Comparison of NDWI results with visual interpretation area

2.2 冰面湖面積統計

對二值化后(這里是指圖像的二值化)的研究區域影像進行預處理,得到冰架冰面湖區域矢量圖。格式轉換后,通過統計冰架研究區域的像元面積和像元數量,就可以得到每一年冰架的冰面湖面積。通過分析統計每一年的冰面湖面積值,可以獲得冰架表面冰面湖的動態變化情況。

2.3 相關性分析

分別對埃默里冰架與拉森A冰架冰面湖數據、氣溫數據、海溫數據進行時間序列分析,總結兩區域冰面湖變化特征,分析區域氣候變化規律,并利用皮爾遜相關系數將冰面湖數據分別與氣溫數據、海溫數據進行遙相關分析,探討氣候變化與冰面湖之間的聯系。皮爾遜相關系數計算如下:

式中,X、Y代表輸入變量,E為數學期望,D為方差,Cov(X,Y)為兩個變量的協方差。

3 氣候變化對冰面湖的影響分析

3.1 氣溫變化

冰面湖的形成多發生于1—2月,本文著重分析1—2月的氣溫變化。利用中山站與長城站2000—2019年1—2月的地面站氣溫數據,計算得到每年1—2月的平均溫度值并繪制曲線圖。圖3顯示,2000—2019年間,中山站月均氣溫最低的年份發生于2000年,為–3.1℃,月均氣溫最高的年份則為2003年,溫度0.05℃。2000—2003年,氣溫有明顯上升趨勢,2003年后氣溫開始下降,直至2008年,氣溫達到–2.45℃,下降趨勢消失,并出現為期4年的波動。2011—2015年,氣溫逐漸上升,直至2015年,氣溫再次出現波動,但在2011年之后,最低氣溫下限一直在上升,這表明該區域氣溫總體呈上升趨勢。圖中趨勢線為該區域氣溫的二次多項式擬合,二次項系數為負,開口向下。趨勢線顯示該區域在2000—2019年,氣溫存在上升趨勢,但增速趨緩。需要說明的是,文中圖3—圖8中的趨勢線都為數據的二次多項式擬合。

長城站氣溫1—2月月均氣溫有輕微降低趨勢(圖4),但下降速度趨緩。2000—2009年氣溫呈階段性上升,2010年可能經歷了極端天氣,氣溫出現短暫的急劇下降,達到這 20年來最低值0.35℃,在2011年回暖之后,氣溫再次出現波動。趨勢線二次項系數為正,開口向上,同樣在2010年出現最低值,趨勢線總體顯示該區域氣溫存在先減小后上升的趨勢。

圖3 中山站2000—2019年1—2月月均氣溫變化Fig.3.Average temperature change at Zhongshan Station from January to February of 2000—2019

圖4 長城站2000—2019年1—2月月均氣溫變化Fig.4.Average temperature change at Great Wall Station from January to February of 2000—2019

3.2 海溫變化

2000—2017年埃默里冰架海域1—2月月均海溫呈上升趨勢(圖5),2007年為月均海溫最低的一年,只有–1.843℃,最高的月均海溫發生在2017年,為–1.712℃。2000—2007年,氣溫一直在上下波動,周期約為3—4年。2008—2017年,除2011和2014年海溫有所下降外,其余年份皆上升。圖5中趨勢線二次項系數為正,且整體同樣呈上升趨勢。拉森A區域海溫上升趨勢也較為明顯(圖6),最低海溫出現在2000年,為–0.217℃,在經歷2000—2006年的海溫快速上升后,2006年達到18年來最高月均海溫0.643℃,期間只有2004年有小幅度下降。2006—2012年,海溫經歷了兩次下降,周期為3—4年。2013—2017年,海溫在經歷2014年的短暫下降后再次回升。趨勢線二次項系數為負,開口向下,趨勢線顯示,在2009年前,該區域海溫增加趨勢較為明顯,增速較快,但在2009年后,海溫增速明顯減緩。

圖5 埃默里冰架海域2000—2017年1—2月月均海溫變化Fig.5.Average SST change on Amery Ice Shelf from January to February of 2000—2017

圖6 拉森A冰架2000—2017年1—2月月均海溫變化Fig.6.Average SST change on Larsen A Ice Shelf from January to February of 2000—2017

3.3 冰面湖面積變化

埃默里冰架2004—2019年冰面湖面積年際變化呈減少趨勢(圖7),其中2005年冰面湖面積最大,達1 219 km2,2005—2007年,面積出現明顯下降,2007—2012年,冰面湖面積一直維持在較低的水平,并有一定程度的下降。2012—2019年,冰面湖面積呈現一定的增加趨勢,2019年冰面湖面積再次突破1 000 km2,達到1 144 km2。圖7趨勢線二次項系數為正,開口向上,R2為0.4534,說明擬合較好,趨勢線顯示埃默里冰面湖面積先減后增,且在2012 后,變化速度有增加趨勢。2003—2016年,拉森A冰面湖面積整體呈增加趨勢(圖8)。2003—2006年,冰面湖面積有明顯的增加趨勢,2006年后,則開始下降,直到2012年下降趨勢趨緩,之后再次快速增加。圖8趨勢線二次項系數為正,開口向上,R2為0.332 6,趨勢線顯示拉森A冰面湖面積增加趨勢明顯,且增速可能會有所提升。

圖7 埃默里冰架2004—2019年冰面湖年際變化Fig.7.Annual variation of supraglacial lakes on Amery Ice Shelf from 2004 to 2019

圖8 拉森A冰架2003—2016年冰面湖年際變化Fig.8.Annual variation of supraglacial lakes on Larsen A Ice Shelf from 2003 to 2016

3.4 冰面湖面積變化與區域氣溫和海溫相關性分析

3.4.1 埃默里冰架

埃默里冰面湖面積與區域氣溫的相關性分析顯示,2004—2019年兩者變化趨勢并不完全相符(圖9),相關系數為0.377 52(圖10a)。2004—2012年,冰面湖與氣溫相關系數為0.214 36(圖10b),2010年后,冰面湖面積與溫度的變化高度一致,呈正相關關系,相關系數0.786 87,達顯著水平(P=0.04<0.05)(圖10c)。2004—2017年,埃默里冰面湖面積與區域海溫一致性較差(圖11),相關系數只有0.004(圖12a),2004—2012年,兩者相關系數為0.051 36(圖12b),2012年后相關系數也只有0.003(圖12c),均未超過0.1。綜合以上分析,說明該冰架冰面湖對氣溫變化敏感,但與海溫間的關系并不明顯。

上述觀點也與Golledge 等[32]的認識一致,區域尺度的風型和局部冰反照率可以很好地解釋這一現象[16,33]。在冰架接地線附近低海拔地區,經常受到來自冰蓋內部的持續的下降風影響,當它向下吹時會使空氣變暖并與之混合,導致夏季近地表的溫度可能比冰上或冰下的地區高出3℃,接近地面的融水產量比其他地區高出一倍[18];另一方面低反照率的藍冰、冰原島峰和冰屑因太陽能吸收的增加而加速融化[34]。藍冰來源于雪的升華、融化以及風蝕作用[35],并且通常與冰原島峰相鄰[36]。崎嶇的地形會產生疾風,而低反照率的巖石會使氣溫升高[37],冰融和風蝕作用擴大了裸巖的面積,導致了升溫、冰融、裸巖和藍冰之間的耦合。但形成在基巖地形和冰流相互作用的地表洼地上的冰面湖[38],面積受非氣候變化因素影響,這可能使得2010年前后冰湖對氣溫變化的響應產生差異。

海溫作用不能有效穿越冰架基底和邊緣可能是海溫與埃默里冰面湖相關性較差的重要原因。首先,埃默里厚厚的冰層會限制基底融化對冰面湖的影響,海水在冰架基底冷卻的過程中,隨著熱量與鹽分的析出會產生高鹽低溫的冰架水[39],這種水團溫度在–1.9℃以下,低于區域海溫[40],作為南極第三大灣的普里茲灣,其海灣地形雖然限制了冰架水與開闊海洋之間的交換,但也會促進上述過程[41];其次,因接地線附近的冰面湖沿排水系統向下流向冰架內部,融冰的大量產生會使得冰架邊緣溫度下降[42],這會阻礙海溫透過冰架邊緣向冰架內部傳導。

圖9 2004—2019年埃默里冰架冰面湖面積變化與區域內1—2月月均氣溫變化對比Fig.9.Comparison of the area changes of supraglacial lakes on Amery Ice Shelf and the average temperature in the region from January to February of 2004 to 2019

圖10 埃默里冰架冰面湖面積與區域內1—2月月均氣溫相關性分析.a)2004—2019;b)2004—2012;c)2012—2019Fig.10.Correlation analysis of the area change of supraglacial lakes on Amery Ice Shelf with average temperature in the region from January to February.a)2004—2019;b)2004—2012;c)2012—2019

圖11 2004—2017年埃默里冰架冰面湖面積變化與區域內1—2月月均海溫變化對比Fig.11.Comparison of the area changes of supraglacial lakes on Amery Ice Shelf and the average SST change in the region from January to February of 2004—2017

3.4.2 拉森A冰架

拉森A冰面湖與長城站氣溫數據分析結果表明(圖13),2005—2008年,兩者雖存在短暫的一致,但兩條曲線變化趨勢總體存在較大差異,相關系數只有0.0140(圖14a)。2016年冰面湖數據與之前數據相差四年,所以在去除2016 數據后,再次進行相關分析,相關系數為0.370 56(圖14b),依舊為弱相關。拉森A冰面湖面積與1—2月月均海溫進行對比(圖15),在2003—2008年,兩者表現高度一致,雖然2009年冰面湖數據缺失，但可以發現2008—2012年海溫與冰面湖面積皆經歷了先下降后上升的過程，同樣具有一致性，因2012—2016年，冰面湖面積數據存在缺失較多，所以該時段二者趨勢無法判斷。2003—2016年兩者相關系數為0.290 37(圖16a),但在2003—2012年兩者相關系數達到0.952 37(極顯著相關,P=0.003<0.01,圖16b),說明拉森A冰面湖面積的變化主要受海溫變化的影響。

圖12 埃默里冰架冰面湖面積與區域內1—2月月均海溫相關性分析.a)2004—2017；b)2004—2012；c)2012—2017Fig.12.Correlation analysis of the area change of supraglacial lakes on Amery Ice Shelf with average SST in the region from January to February.a)2004—2017;b)2004—2012;c)2012—2017

拉森A 周邊海域水團眾多[43],而拉森A冰架作為海洋性冰蓋,會受到水團的劇烈影響。有研究表明拉森冰架融化與1970年以來威德爾海深層水變暖有關[44]。Nicholls 等[45]通過觀察冰架底部的流出物,推導出冰架水是改進的威德爾海深層水與大氣相互作用后冷卻至表面冰點的結果,這證明了溫暖的洋流可以進入冰架腔洞,進而加劇拉森冰架的融化。Nicholls 等觀察的流出物來自整個冰架,所以他們的推論是普遍適用的,從而在拉森冰架A和威德爾海洋條件之間建立起聯系,Holland 等[46]的模型也證明了這一點。除此之外,由地表融水融化冰架產生的熱鹽環流也會促進冰架的融化[47]。

圖13 2003—2019年拉森A冰架冰面湖面積變化與區域內1—2月月均氣溫變化對比Fig.13.Comparison of the area changes of supraglacial lakes on Larsen A Ice Shelf and the average SST change in the region from January to February of 2003—2019

圖14 拉森A冰架冰面湖面積與區域1—2月月均氣溫相關性分析.a)2003—2019;b)2003—2012Fig.14.Correlation analysis of the area change of supraglacial lakes on Larsen A Ice Shelf with average temperature in the region from January to February.a)2003—2019;b)2003—2012

圖15 拉森A冰架冰面湖面積變化與區域內1—2月月均海溫變化對比Fig.15.Comparison of the area changes of supraglacial lakes on Larsen A Ice Shelf and the average SST change in the region from January to February of 2003—2017

南極半島周圍的冰架在最近幾十年不斷地后退,人們普遍認為與快速的大氣變暖有關[48],但拉森A冰面湖對區域氣溫變化并不敏感,這與南極西部松島冰川(the Pine Island Glacier)和斯韋茨冰川(Thwaites Glacier)的研究成果一致[16]。King等[49]就焚風對拉森C冰架表面融化影響進行了研究,發現焚風對冰架表面融化具有重要影響,而焚風的強度由氣溫和相對濕度決定。說明氣溫對拉森A冰架表面直接影響較弱,但可與其他氣候因素綜合作用于冰面湖的形成。

圖16 拉森A冰架冰面湖面積與區域內1—2月月均海溫相關性分析.a)2003—2016;b)2003—2012Fig.16.Correlation analysis of the area change of supraglacial lakes on Larsen A Ice Shelf with average SST in region from January to February.a)2003—2016;b)2003—2012

4 結論與討論

本文利用中山站與長城站的2000—2019年氣溫、NOAA 2000—2017年海表面溫度數據以及遙感影像提取的冰面湖數據,對埃默里冰架與拉森A冰架區域內的氣溫、海溫、冰面湖變化進行了趨勢分析,在此基礎上分析了氣溫、海溫變化對冰面湖變化的影響。研究發現埃默里冰面湖對氣溫變化更為敏感,海溫在埃默里冰面湖形成中的影響卻并不顯著,這可能與南極東西部冰蓋下基巖地形的幾何形狀、冰層厚度的區域差異以及冰架水和融冰的綜合作用有關。拉森A冰面湖則相反,因為威德爾海深層水的存在而與區域海溫聯系較為密切,但其與區域氣溫的關系并不顯著。

埃默里冰架區域隨著氣溫上升,涌入普里茲灣的冰流和冰川融水將會增加,這很大程度會影響普里茲灣的熱鹽環流,使海溫下降,普里茲灣一直被認為是南極底層水的潛在生成地,普里茲灣海溫下降,很可能影響南極環流溫度,緩解氣候變暖。兩冰架區域海溫升高是整個南極周邊海溫升高的縮影,升高的海溫一方面會進一步促進海溫敏感區的冰川融化,加速海平面上升;另一方面通過與南大洋的物質交換,對世界范圍的熱鹽平衡產生深遠影響。

由于南極影像數據的缺乏,本研究主要采用了中等分辨率的Landsat 影像和產品的冰面湖面積信息,未來將會綜合空間分辨率更高的Sentinel 10 m影像對冰面湖數據集進行延長,并嘗試尋找更多的遙相關證據,探索氣候變化對冰面湖更為全面具體的影響機制。因為氣象觀測資料長度有限,區域氣候變化更為長期的變化規律以及對冰面湖變化的影響需要搜集更多資料做進一步研究,后續研究將使用更多的氣象資料,特別是再分析資料,并選擇更多的樣區來開展進一步工作。