沈心儀 張瑜,2,3 陳長勝4,,3 胡松,3

研究論文

加拿大北極群島區(qū)域西北航道海冰冰情長期時空變化特征研究

沈心儀1張瑜1,2,3陳長勝4,1,3胡松1,3

(1上海海洋大學(xué)海洋科學(xué)學(xué)院, 上海 201306;2南方海洋科學(xué)與工程廣東省實驗室(珠海), 廣東 珠海 519082;3上海海洋大學(xué)國際海洋研究中心, 上海 201306;4馬薩諸塞大學(xué)達特茅斯分校海洋科學(xué)與技術(shù)學(xué)院, 馬薩諸塞州 新貝德福德市 02744, 美國)

加拿大北極群島西北航道區(qū)域海冰冰情對西北航道的開通及航行安全有著顯著的影響。本研究將西北航道分為南部、北部以及交匯區(qū)3個區(qū)域, 利用美國國家冰雪中心Bootstrap海冰密集度數(shù)據(jù)和CryoSat-2、CS2SMOS觀測以及高分辨率北極海洋-海冰耦合模型(AO-FVCOM)的海冰厚度數(shù)據(jù)對1978年11月—2017年12月加拿大北極群島區(qū)域西北航道海冰冰情長期時空特征進行研究, 并對南線和北線冰情進行評估。研究結(jié)果顯示西北航道冬春季被高密集度的海冰覆蓋, 夏秋季密集度較小, 呈現(xiàn)北高南低的分布特征; 北部、南部、以及交匯區(qū)域覆蓋范圍呈現(xiàn)減少趨勢, 分別為–0.01×105km2·(10a)–1[–0.77%·(10a)–1]、–0.09×105km2·(10a)–1[–3.50%·(10a)–1]以及–0.12×104km2·(10a)–1[–2.81%·(10a)–1]; 海冰厚度在北部和南部區(qū)域呈現(xiàn)減小趨勢, 分別為–0.25 m·(10a)–1和–0.13 m·(10a)–1。西北航道南線冰情優(yōu)于北線, 海冰密集度、覆蓋范圍、厚度均小于北線。對影響海冰密集度和厚度的大氣和海洋熱力學(xué)因素進行分析發(fā)現(xiàn), 研究區(qū)域的表面大氣溫度、海洋表層溫度呈現(xiàn)上升趨勢, 并和海冰密集度呈現(xiàn)顯著負相關(guān), 和北部、南部區(qū)域海冰厚度呈現(xiàn)顯著負相關(guān), 與交匯區(qū)域海冰厚度相關(guān)不顯著。

加大拿北極群島 西北航道 海冰 海表氣溫 海表溫度

0 引言

西北航道是由加拿大東北部戴維斯海峽(Davis Strait)經(jīng)加拿大北極群島(Canadian Arctic Archipelago, CAA)到美國阿拉斯加北岸的航道, 是歐亞大陸東西兩側(cè)連接太平洋和大西洋之間的最短航線, 相較于傳統(tǒng)歐洲通往太平洋的26 000余公里的巴拿馬運河航線, 西北航道航程縮短近一半, 約為14 000公里, 但是歷史上由于冰情嚴重和船舶通航能力限制, 西北航道一直難于通行。

隨著全球變暖, 北極海冰消融。海冰覆蓋范圍[1-2]和厚度[3-4]都呈現(xiàn)減小趨勢, 其中2012年9月13日北極海冰覆蓋范圍和面積降到歷史新低值, 分別是3.4×106km2和3.0×106km2[2], 北極海冰年平均厚度從1975年的3.59 m減少至2012年的1.25 m, 共減少65%[3]。作為北極的重要區(qū)域之一, 加拿大北極群島的海冰面積在年際變化上也呈現(xiàn)減少趨勢, 其減少速率約為(0.7%±0.6%)·(10a)–1,其中, 夏季約(1.4%±2.2%)·(10a)–1的下降速率, 冬季約(1.2%±0.2%)·(10a)–1的增長速率[5]。在此背景下, 加拿大北極群島作為西北航道的關(guān)鍵區(qū)域, 其夏季海冰消融給西北航道通航帶來更多的可能性[6], 隨著航道的開通將會帶來巨大的經(jīng)濟效益。

西北航道在穿過加拿大北極群島時, 有多條可行路線(圖1)?？勘钡囊粭l航線由麥克盧爾海峽(M’Clure Strait)經(jīng)過梅爾維爾子爵海峽(Viscount Melville Sound)、巴羅海峽(Barrow Strait)和蘭開斯特海峽(Lancaster Sound)到達巴芬灣(Baffin Bay)?？磕系暮骄€則有多條路線, 均起始于阿蒙森灣(Amundsen Gulf), 一條經(jīng)過威爾士王子海峽(Prince of Wales Strait)在Viscount Melville Sound與北線匯合; 另外3條航線則穿過科羅內(nèi)申灣(Coronation Gulf), 在通過維多利亞海峽(Victoria Strait)后分開, 其中一條穿過麥克林托克海峽(M’Clintock Channel)在Viscount Melville Sound與北線匯合, 另一條穿過皮爾海峽(Peel Sound)在Barrow Strait與北線匯合, 還有一條穿過利金特王子灣(Prince Regent Inlet)在Lancaster Sound與北線匯合。

圖1 加拿大北極群島區(qū)域地形圖及西北航道路線圖

Fig.1. The domain of Canadian Arctic Archipelago and the routes of the Northwest Passage

近些年來, 隨著世界經(jīng)濟貿(mào)易需求增加和我國冰上絲綢之路倡議的提出, 西北航道受到了各方的不斷關(guān)注, 國內(nèi)外一些針對西北航道冰情以及對通航條件影響的研究已開展。蘇潔等[7]與李春花等[8]通過利用Advanced Microwave Scanning Radiometer for the Earth Observing System, AMSR-E)海冰密集度數(shù)據(jù)研究發(fā)現(xiàn)西北航道南線海冰冰情優(yōu)于北線, 北線冰情西重東輕, 南線冰情東重西輕。影響西北航道北線開通的主要因素包括北部伊麗莎白女王群島(QEI)輸入的海冰、局地氣溫變化和大西洋暖流。西北航道9月海冰冰情都呈現(xiàn)減輕趨勢[8]。2006年Howell等[9]通過分析研究快速散射計(Quick Scatterometer, QuikSCAT)數(shù)據(jù)反演加拿大北極群島2002—2005年海冰融化情況, 指出在過去的5年中, 雖然融冰期延長, 但西北航道通航的時間并沒有顯著延長。Howell等[9-10]2006和2008年的研究指出西北航道內(nèi)存在由加拿大北極群島北部輸入的多年冰, 且多年冰一旦流入可能會連續(xù)存在幾年, 影響西北航道的通行[10]。2009年, Howell等[11]分析1979—2008年掃描式多通道微波輻射計(Scanning Multichannel Microwave Radiometer, SMMR)和專用傳感器微波/成像儀(Special Sensor Microwave/Imager, SSM/I)數(shù)據(jù)發(fā)現(xiàn)加拿大北極群島9月海冰覆蓋面積呈現(xiàn)強烈減少趨勢, 減少速率為–8.7%·(10a)–1, 融冰期持續(xù)延長。2013年Howell等[12]分析1997—2012年5—11月加拿大北極群島與北冰洋之間的海冰面積通量, 發(fā)現(xiàn)2007、2008、2010和2011年西北航道通航可部分歸因于M’Clure Strait海冰輸入?yún)T乏。除了根據(jù)衛(wèi)星數(shù)據(jù)進行研究外, Melling[13]利用1971—1980年鉆孔海冰厚度數(shù)據(jù)研究發(fā)現(xiàn)天氣變暖, 加拿大北極群島低密集度海冰融化, 開闊水域增加, 加拿大海盆的多年冰可以順利通過原本被海冰堵塞的通道進入QEI以及M’Clure Strait進而到達加拿大北極群島南部地區(qū), 導(dǎo)致冰情可能會更加嚴重。此外, 2015年Hass和Howell[14]使用航空電磁方式測量得到的2014和2015年4月和5月一年冰和多年冰海冰厚度, 結(jié)果表明在西北航道內(nèi)觀測到了來自北冰洋平均厚度在3 m以上的海冰, 西北航道內(nèi)的海冰冰情仍然較為嚴峻。

加拿大北極群島岸線和地形復(fù)雜, 海冰厚度觀測數(shù)據(jù)在空間和時間上都有較大的缺失, 前人對于西北航道冰情的研究主要集中在海冰密集度研究上, 這僅針對無破冰能力的船只而言, 對于有破冰能力的船只, 海冰厚度則是另一重要的考慮因素, 在這方面現(xiàn)有的研究還較少。因此本文綜合使用海冰密集度、海冰厚度遙感觀測數(shù)據(jù), 并結(jié)合高分辨率北極海冰耦合模式的海冰厚度數(shù)據(jù)有效填補海冰厚度數(shù)據(jù)在時間上和空間上的空白, 以此對西北航道海冰進行研究。本研究選擇西北航道最為重要的兩條航線, 根據(jù)海洋環(huán)流和海冰環(huán)流的方向, 定義Amundsen Gulf和M’Clure Strait為入口, Lancaster Sound為出口, 其中靠南的一條航線從Amundsen Gulf穿過Coronation Gulf, Victoria Strait, Peel Sound, Barrow Strait和Lancaster Sound最終到達巴芬灣(下文統(tǒng)稱南線), 這條線是在各種定義中使用最廣泛的, 本研究以此作為南線的代表。另一條靠北的航線由M’Clure Strait經(jīng)過Viscount Melville Sound, Barrow Strait和Lancaster Sound 到達巴芬灣(下文統(tǒng)稱北線)。而從Amundsen Gulf穿過Prince of Wales Strait在Viscount Melville Sound和北線匯合的航線, 由于航行路線與北線大多數(shù)重合, 所以本研究將Prince of Wales Strait區(qū)域并入北線討論。

根據(jù)海冰密集度和厚度的季節(jié)性和年際性變化情況, 我們對南線和北線冰情的長期變化進行評估, 在此基礎(chǔ)上, 我們對影響海冰冰情的熱力學(xué)因素包括表面大氣溫度(Surface Air Temperature, SAT)和海洋表層溫度(Sea Surface Temperature, SST)進行分析討論。

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

1.1 海冰觀測數(shù)據(jù)

鑒于西北航道區(qū)域地形較為復(fù)雜, 本研究將南北兩條航線所經(jīng)過的區(qū)域分為3個子區(qū)域來進行討論, 研究其海冰分布以及季節(jié)、年際和月際變化趨勢: 北部區(qū)域(圖1粉色, M’Clure Strait、Viscount Melville Sound、Prince of Wales Strait)、南部區(qū)域(圖1黃色, Amundsen Gulf、Coronation Gulf、Victoria Strait、Peel Sound)、交匯區(qū)域(圖1紫色, Barrow Strait、Lancaster Sound)。本文研究所用到的Bootstrap逐月海冰密集度數(shù)據(jù)來自美國國家冰雪中心(National Sea and Ice Data Center, NSIDC, http://nsidc.org/data/nsidc-0079), 時間范圍為1978年11月—2017年12月, 分辨率為25 km× 25 km。根據(jù)每個網(wǎng)格點的海冰密集度數(shù)據(jù)和每個網(wǎng)格點的控制面積, 計算研究區(qū)域海冰覆蓋范圍(sea ice extent), 具體公式如下:

其中,S為海冰覆蓋范圍,w為權(quán)重系數(shù),q為海冰密集度,A為網(wǎng)格面積。

本研究使用了2種海冰厚度觀測數(shù)據(jù), 其一來自CryoSat-2雷達測高數(shù)據(jù)(https://www. meereisportal.de/en.html), 數(shù)據(jù)時間范圍是2010年11月—2017年12月(缺失5—9月數(shù)據(jù)), 空間分辨率為25 km×25 km, 通過美國國家航空航天局(National Aeronautics and Space Administration, NASA)冰橋行動2011—2013年每年3月冰橋和CryoSat-2海冰厚度數(shù)據(jù)對比, 平均誤差為0.182 m[15]; 其二是CS2SMOS海冰厚度融合數(shù)據(jù)(https:// data.meereisportal.de/data/cs2smos_awi/n/), CS2S-MOS周平均海冰厚度數(shù)據(jù)融合了CryoSat-2和SMOS海冰厚度數(shù)據(jù), 時間范圍為2010年10月—2017年12月, 空間分辨率為25 km×25 km。由于CryoSat-2和CS2SMOS海冰厚度數(shù)據(jù)并不能完全覆蓋研究區(qū)域并且覆蓋的區(qū)域不同, 再加上5—9月海冰厚度數(shù)據(jù)缺失, 導(dǎo)致空間分布以及時間序列不完整, 因此模式數(shù)據(jù)將為西北航道海冰厚度研究提供必不可少的幫助。

1.2 北極海洋-海冰耦合模式海冰厚度數(shù)據(jù)

本研究選用高分辨率北極有限體積海洋-海冰耦合模型(Arctic Ocean-Finite Volume Community Ocean Model, AO-FVCOM)[16-17]的海冰厚度數(shù)據(jù), 來研究更長時間段1978年11月—2017年12月西北航道海冰厚度空間分布、季節(jié)及長期變化特征。模型空間分辨率為2~40 km, 在加拿大北極群島區(qū)域水平分辨率達2 km, 能夠很好地模擬加拿大北極群島復(fù)雜的岸線和地形, 垂向分辨率為45層, 模型驅(qū)動包括潮汐、風應(yīng)力、凈熱通量、海表面大氣壓、蒸發(fā)降水、河流徑流等。2016年Zhang等[1]已將AO-FVCOM海冰厚度數(shù)據(jù)與北極地區(qū)多源海冰厚度數(shù)據(jù)包括衛(wèi)星觀測、現(xiàn)場鉆孔觀測、機載電磁觀測、海冰冰站觀測等進行了詳細驗證對比, 結(jié)果表明AO-FVCOM模式的海冰厚度數(shù)據(jù)能夠很好地捕捉到可觀測的北極海冰厚度數(shù)據(jù)空間分布、季節(jié)以及年際變化特征, 此外, 與其他6個海冰模式數(shù)據(jù)的綜合對比中, AO-FVCOM與多源海冰厚度觀測數(shù)據(jù)在厚度誤差以及季節(jié)和年際變化的相關(guān)性上也均展現(xiàn)了較好的結(jié)果。

1.3 大氣和海洋熱力學(xué)數(shù)據(jù)

本研究使用歐洲中期天氣預(yù)報中心(Euro-pean Centre for Medium-Range Weather Forecasts, ECMWF, https://cds.climate.copernicus.eu/)的月平均表面大氣溫度(SAT)、海洋表層溫度(SST)再分析數(shù)據(jù)資料來進一步分析研究西北航道海冰冰情及其相互影響。表面大氣溫度和海洋表層溫度數(shù)據(jù)的時間范圍為1979—2017年, 空間分辨率為0.25°×0.25°。

2 海冰密集度時空變化特征

2.1 海冰密集度空間分布特征

如圖2所示, 1979—2017年整個西北航道內(nèi)海冰密集度在不同月份呈現(xiàn)明顯的空間分布差異。海冰密集度代表每個網(wǎng)格點內(nèi)海冰覆蓋面積占網(wǎng)格總面積的比率, 當海冰密集度為1時, 海冰覆蓋整個網(wǎng)格點, 當海冰密集度為0時, 網(wǎng)格點內(nèi)無海冰覆蓋。研究區(qū)域海冰密集度存在顯著的季節(jié)變化, 冬春高, 夏秋低。研究區(qū)域在12月至次年4月被高密集度大于0.95的海冰覆蓋, 此時兩條航線都因海冰封閉, 缺少破冰能力的船只不宜通航。5—11月西北航道區(qū)域海冰密集度呈現(xiàn)北高南低的分布特征。5月海冰從3個通道: Amundsen Gulf、Lancaster Sound和M’Clure Strait開始融化, Amundsen Gulf海冰從通道口向東南部融化, 通道口處海冰密集度小于0.8, 靠近Coronation Gulf區(qū)域的海冰密集度仍高達1.0; Lancaster Sound海冰密集度約為0.8; M’Clure Strait通道口處海冰開始融化, 密集度約為0.9。6月研究區(qū)域海冰密集度整體下降, 上述3個通道口海冰融化最為顯著, 其中Lancaster Sound和Amundsen Gulf區(qū)域存在0.6以下的海冰密集度, M’Clure Strait通道口海冰密集度降至0.8以下。隨著海冰持續(xù)融化, 7月僅有Peel Sound和Viscount Melville Sound仍然存在大面積海冰密集度0.8以上的區(qū)域。夏末秋初的8月與9月海冰密集度在整個研究區(qū)域分布呈現(xiàn)一個低值。西北航道北部區(qū)域海冰密集度高, 大范圍存在密集度大于0.6的海冰。西北航道南部區(qū)域以及交匯區(qū)域海冰密集度低于0.6, 其中Amundsen Gulf靠近Prince of Wales Strait區(qū)域存在無冰區(qū)。研究區(qū)域10月開始結(jié)冰, Peel Sound, Prince of Wales Strait, M’Clure Strait, Viscount Melville Sound, Barrow Strait結(jié)冰速度較快, 其10月平均海冰密集度高于0.8。11月除Amundsen Gulf結(jié)冰速度較慢、海冰密集度明顯小于其他區(qū)域以外, 其他區(qū)域已接近冬季數(shù)值。

2.2 海冰覆蓋范圍時間變化特征

2.2.1 季節(jié)變化特征

為了進一步了解西北航道區(qū)域冰情, 本文對區(qū)域內(nèi)的海冰覆蓋范圍開展了研究。圖3分別給出了1979—2017年3個區(qū)域的海冰覆蓋范圍的季節(jié)變化, 北部區(qū)域總面積為1.54×105km2, 南部區(qū)域為2.51×105km2, 交匯區(qū)域為4.46×104km2。由于3個區(qū)域總面積不同, 本文計算了海冰覆蓋率用于研究分析3個區(qū)域海冰覆蓋情況。海冰覆蓋率的大小為海冰覆蓋范圍占區(qū)域總面積比值。冬春季研究區(qū)域基本被海冰覆蓋, 夏秋季海冰覆蓋范圍較小, 9月海冰覆蓋范圍達到最小值。3個區(qū)域中, 北部區(qū)域海冰覆蓋范圍季節(jié)變化較小, 9月最小海冰覆蓋范圍也高達1.43×105km2, 北部區(qū)域92.59%仍被海冰覆蓋; 南部區(qū)域和交匯區(qū)域海冰覆蓋范圍季節(jié)變化較大, 9月海冰最小覆蓋范圍分別為1.33×105km2和3.35×104km2, 其海冰覆蓋率達到52.92%和75.11%。

圖2 加拿大北極群島區(qū)域西北航道1979—2017年海冰密集度季節(jié)分布圖(1—12月)

Fig.2. Seasonal distribution of sea ice concentration in the Northwest Passage of the Canadian Arctic Archipelago over the period 1979—2017 (January—December)

南部區(qū)域最先出現(xiàn)無冰區(qū), 5月海冰覆蓋范圍減少, 而北部區(qū)域和交匯區(qū)域6月海冰覆蓋范圍開始減少。北部區(qū)域結(jié)冰速度較快, 10月海冰再次覆蓋整個北部區(qū)域, 而南部區(qū)域以及交匯區(qū)域則于11月被海冰覆蓋, 南部區(qū)域結(jié)冰速度慢于交匯區(qū)域, 其10月海冰覆蓋范圍分別為2.51×105km2(86.48%)和4.46×104km2(99.14%)。

圖3 1979—2017年北部區(qū)域、南部區(qū)域、交匯區(qū)域海冰覆蓋范圍季節(jié)變化. 圖中藍線表示多年月平均海冰覆蓋范圍、黑色豎棒表示標準差

Fig.3. Seasonal variation of sea ice extent of the northern area, southern area and intersectional area over the period of 1979—2017. The blue lines are the multi- year monthly sea ice extent and the vertical bars are the standard deviations

2.2.2 年際和月際變化特征

從1979—2017年研究區(qū)域海冰覆蓋范圍年變化(圖4)來看, 北部、南部、交匯區(qū)域的多年平均海冰覆蓋范圍分別為1.53×105km2、2.28×105km2以及4.24×104km2。北部區(qū)域海冰覆蓋率最高, 為99.06%; 其次為交匯區(qū)域, 為95.06%; 南部區(qū)域最低, 為90.51%。南部區(qū)域以及交匯區(qū)域的海冰覆蓋范圍下降趨勢大于北部區(qū)域, 3個區(qū)域的下降趨勢分別是–0.09×105km2·(10a)–1[–3.50%·(10a)–1], –0.12×104km2·(10a)–1[–2.81%·(10a)–1]以及–0.01× 105km2·(10a)–1[–0.77%·(10a)–1](以上趨勢通過99%的顯著性檢驗)。在年際尺度上, 北部、南部及交匯區(qū)域的海冰覆蓋范圍最小值分別發(fā)生在2012、1998及2011年, 其海冰覆蓋范圍分別是1.41×105、1.99×105和3.84×104km2, 覆蓋率分別達到91.60%、79.32%和86.10%。

圖4 1979—2017年北部區(qū)域、南部區(qū)域、交匯區(qū)域年平均海冰覆蓋范圍(藍色柱狀). 圖中紅線表示線性回歸趨勢

Fig.4. Yearly sea ice extent (blue bars) of the northern area, southern area and intersectional area over the period of 1979—2017. The red lines are the linear regression trends

根據(jù)1978年11月—2017年12月研究區(qū)域月平均海冰覆蓋范圍的長期變化(圖5)來看, 1979—1997年北部區(qū)域全年幾乎被海冰覆蓋, 多年平均值為1.54×105km2(99.98%), 南部區(qū)域和交匯區(qū)域存在下降趨勢, 但整體海冰覆蓋范圍較大, 多年平均值分別為2.36×105km2(93.73%)和4.38× 104km2(98.30%); 自1998年起3個區(qū)域海冰覆蓋范圍低值頻發(fā), 北部區(qū)域、南部區(qū)域和交匯區(qū)域多年平均值分別為1.51×105km2(98.20%)、2.20× 105km2(87.46%)和4.10×104km2(91.98%), 北部區(qū)域以及南部區(qū)域海冰覆蓋范圍最小值發(fā)生在2012年9月與海冰密集度最小值有很好的對應(yīng), 覆蓋范圍分別為0.59×105km2, 0.00×105km2, 交匯區(qū)域海冰覆蓋范圍最小值發(fā)生于1998年9月, 覆蓋范圍為0.17×104km2。

圖5也顯示了3個區(qū)域海冰覆蓋范圍的季節(jié)變化, 冬春季3個區(qū)域整體被海冰覆蓋, 海冰覆蓋范圍并沒有下降趨勢, 夏秋季海冰覆蓋范圍較小, 并且3個區(qū)域海冰覆蓋范圍呈現(xiàn)顯著的下降趨勢。由于3個區(qū)域總面積不一致, 本文計算了海冰覆蓋范圍變化率來展示每個區(qū)域海冰覆蓋范圍變化相對大小。海冰覆蓋變化率為海冰覆蓋范圍變化趨勢與1978年11月海冰覆蓋范圍的比值。根據(jù)表1, 北部區(qū)域和南部區(qū)域, 秋季海冰覆蓋范圍下降趨勢略大于夏季分別為–0.03× 105km2· (10a)–1[–2.02%·(10a)–1, 通過99%顯著性檢驗], –0.20×105km2·(10a)–1[–8.59%·(10a)–1, 通過99%顯著性檢驗], 夏季下降趨勢分別為–0.02× 105km2·(10a)–1[–1.04%·(10a)–1, 通過95%顯著性檢驗], –0.14×105km2·(10a)–1[–5.81%·(10a)–1, 通過99%顯著性檢驗]; 交匯區(qū)域夏季海冰覆蓋范圍下降趨勢(–0.27×104km2·(10a)–1, –6.16%·(10a)–1, 通過99%顯著性檢驗)大于秋季(–0.23× 104km2·(10a)–1, –5.08%·(10a)–1, 通過99%顯著性檢驗)。北部區(qū)域10月至次年7月海冰覆蓋范圍沒有顯著趨勢, 8—9月呈現(xiàn)顯著減小趨勢, 分別是–0.04×105km2·(10a)–1[–1.04%·(10a)–1, 通過95%的顯著性檢驗]和–0.09×105km2·(10a)–1[–6.07%·(10a)–1, 通過99%的顯著性檢驗]。南線區(qū)域11月至次年6月海冰覆蓋范圍沒有顯著趨勢, 7—10月呈現(xiàn)顯著減小趨勢, 其中9月海冰覆蓋范圍減小趨勢最大為–0.43×105km2·(10a)–1[–22.31%·(10a)–1, 通過99%的顯著性檢驗]。交匯區(qū)域11月至次年5月海冰覆蓋范圍沒有顯著趨勢, 6—10月呈現(xiàn)減小趨勢, 其中9月海冰覆蓋范圍減小趨勢最大為–0.65×104km2·(10a)–1[–14.78%·(10a)–1, 通過99%的顯著性檢驗]。結(jié)合3個區(qū)域面積大小, 南部區(qū)域和交匯區(qū)域夏秋季海冰覆蓋范圍減小趨勢大于北部區(qū)域。根據(jù)西北航道無冰狀態(tài)來看, 南線優(yōu)于北線, 較早存在無冰區(qū), 并且夏季海冰覆蓋率低, 并伴隨著較強的下降趨勢。

圖5 1978年11月—2017年12月北部區(qū)域、南部區(qū)域、交匯區(qū)域月平均海冰覆蓋范圍(藍線).圖中紅線表示線性回歸趨勢

Fig.5. Monthly sea ice extent (blue lines) of the northern area, southern area and intersectional area over the period of November 1978—December 2017. The red lines are the linear regression trends

3 海冰厚度時空變化特征

3.1 海冰厚度觀測數(shù)據(jù)時空變化特征

由于CryoSat-2衛(wèi)星測高海冰厚度數(shù)據(jù)和CS2SMOS融合海冰厚度數(shù)據(jù)在研究區(qū)域內(nèi)每個時間段的有效觀測點個數(shù)并不一致, 并且缺失5—9月海冰厚度數(shù)據(jù), 本文挑選了覆蓋范圍最廣的2017年3月、2016年11月、2016年12月作為春季、秋季和冬季代表的海冰厚度空間分布來對研究區(qū)域內(nèi)海冰厚度在不同季節(jié)的分布情況(圖6)進行研究。由圖可知, CryoSat-2在研究區(qū)域內(nèi)有效觀測數(shù)據(jù)點較少, CS2SMOS海冰厚度數(shù)據(jù)有效觀測點相較更多, 但都沒能完全覆蓋研究區(qū)域。基于現(xiàn)有的數(shù)據(jù)來看, 3月海冰厚度最大、12月其次、11月最小; 11月、12月研究區(qū)域內(nèi)海冰厚度分布呈現(xiàn)北高南低的分布特征, 3月CS2SMOS海冰厚度分布特征為北高南低, 而CryoSat-2海冰厚度數(shù)據(jù)顯示在研究區(qū)域北部(M’Clure Strait)和南部(Amundsen Gulf)都存在較厚的海冰。

表1 1979—2017年海冰覆蓋范圍變化趨勢

注: S在99%, 95%, 90%的置信水平上有統(tǒng)計學(xué)意義。

圖6 加拿大北極群島區(qū)域西北航道CryoSat-2與CS2SMOS數(shù)據(jù)海冰厚度分布圖

Fig.6. The distribution of sea ice thickness from CryoSat-2 and CS2SMOS in the Northwest Passage of the Canadian Arctic Archipelago

根據(jù)2017年3月CryoSat-2和CS2SMOS海冰厚度分布圖, 研究區(qū)域3月海冰厚度較大, 但兩套觀測數(shù)據(jù)顯示的分布特征并不相同: 從僅有的有效海冰數(shù)據(jù)來看, CryoSat-2海冰厚度分布顯示在Amundsen Gulf和M’Clure Strait通道口都存在3 m厚的海冰; CS2SMOS海冰厚度數(shù)據(jù)顯示3月西北航道區(qū)域海冰厚度分布呈現(xiàn)北線高南線低的分布特征, M’Clure Strait海冰厚度較大, 約為2.5 m, Amundsen Gulf海冰厚度較小, 約為1.7 m, 其中北線西高東低, 在Lancaster Sound存在小于0.7 m厚度的海冰, 南線在Amundsen Gulf和Victoria Strait海冰較厚, 海冰厚度大于1.4 m。通過分析2016年11月研究區(qū)域的海冰厚度分布可以得出, 相對于3月CryoSat-2海冰厚度有效數(shù)據(jù)范圍增加了Victoria Strait區(qū)域。兩組觀測數(shù)據(jù)在2016年11月的海冰空間分布整體呈現(xiàn)北高南低, 但東西分布特征并不一致: CS2SMOS海冰厚度分布顯示, M’Clure Strait、Viscount Melville Sound海冰厚度達到1.4 m,東部Barrow Strait、Lancaster Sound和南線海冰較薄, 厚度小于0.7 m。而CryoSat-2海冰厚度數(shù)據(jù)則顯示海冰最厚位于Viscount Melville Sound東部, 高達2.8 m。相較11月, 2016年12月海冰厚度增大, CryoSat-2和CS2SMOS海冰厚度分布同樣顯示北線高南線低的分布特征, 但在Barrow Strait和Lancaster Sound兩組觀測海冰厚度數(shù)據(jù)顯示海冰厚度差異較大, 其中CryoSat-2海冰厚度較大, 存在高達3 m的厚冰; CS2SMOS海冰厚度較小, 約為0.7 m。

我們根據(jù)2010年11月—2017年12月研究區(qū)域內(nèi)CryoSat-2和CS2SMOS海冰厚度的有效數(shù)據(jù)(缺失5—9月數(shù)據(jù)), 計算了研究區(qū)域平均海冰厚度(圖7)。由圖可知, 兩套觀測海冰厚度數(shù)據(jù)都顯示了10月為海冰厚度低值, 其后海冰厚度開始增大, 于3月或4月達到最大值。CryoSat-2海冰厚度數(shù)據(jù)的最大值和最小值均大于CS2SMOS海冰厚度。2010年11月—2017年12月(缺失5—9月)CryoSat-2海冰厚度最大值發(fā)生在2017年4月, 為1.903 m, 最小值發(fā)生在2012年10月, 為0.77 m; CS2SMOS海冰厚度最大值發(fā)生在2015年3月, 為1.74 m, 最小值發(fā)生在2011年10月, 為0.20 m。

圖7 2010年11月—2017年12月加拿大北極群島區(qū)域西北航道CryoSat-2(藍線)和CS2SMOS(紅線)海冰厚度月變化

Fig.7. Monthly sea ice thickness from CryoSat-2 (blue lines) and CS2SMOS (red lines) in the Northwest Passage of the Canadian Arctic Archipelago over the period of November 2010—December 2017

3.2 AO-FVCOM模式海冰厚度數(shù)據(jù)空間分布特征

基于CryoSat-2與CS2SMOS的有效觀測數(shù)據(jù)點在加拿大北極群島較少、時間序列上不連續(xù)的原因, 為了研究更長時間序列的研究區(qū)域海冰厚度變化, 我們使用AO-FVCOM模式海冰厚度數(shù)據(jù)結(jié)果進一步分析。圖8為1979—2017年西北航道海冰厚度季節(jié)分布圖。結(jié)果表明西北航道區(qū)域西北、東南部海冰厚度較大, 西南海冰厚度較小。

1—5月研究區(qū)域海冰厚度增加, 5月除了Amundsen Gulf、Prince of Wales Strait以及Coronation Gulf東部存在厚度小于2 m的海冰外, 研究區(qū)域海冰厚度大于2 m, 并且在Victoria Strait西南部和Barrow Strait西部存在厚度大于3 m的厚冰。6月海冰從Amundsen Gulf、M’Clure Strait 兩岸、Lancaster Sound西部及北岸開始變薄, 9月西北航道區(qū)域海冰厚度達到最小, 僅有M’Clure Strait、Viscount Melville Sound和Victoria Strait區(qū)域存在大于0.7 m的冰; 而Amundsen Gulf、Lancaster Sound、Coronation Gulf西部海冰厚度較薄, 小于0.3 m。10月除Victoria Strait區(qū)域外, 西北航道其他區(qū)域開始結(jié)冰, 海冰厚度增加。

圖8 加拿大北極群島區(qū)域西北航道1979—2017年海冰厚度季節(jié)分布(1—12月)

Fig.8. Seasonal distribution of sea ice thickness in the Northwest Passage of the Canadian Arctic Archipelago over the period of 1979—2017 (January—December)

3.3 AO-FVCOM模式海冰厚度數(shù)據(jù)的時間變化特征

3.3.1 季節(jié)變化特征

圖9顯示的是西北航道區(qū)域1979—2017年多年月變化時間序列圖, 3個區(qū)域海冰厚度都呈現(xiàn)出春季高、夏秋低的變化特征, 與觀測數(shù)據(jù)所得到的10—4月的季節(jié)變化特征相同。在3個區(qū)域中, 交匯區(qū)域海冰厚度具有強烈的季節(jié)變化, 海冰厚度在4月達到最大值2.39 m, 9月厚度達到最小值0.07 m。南部區(qū)域海冰厚度較小, 4月最大值為2.21 m, 9月最小值為0.14 m。北部區(qū)域海冰厚度最大值發(fā)生在5月, 為2.35 m; 9月最小值相對較高, 為0.52 m。

圖9 1979—2017年北部區(qū)域、南部區(qū)域、交匯區(qū)域海冰厚度季節(jié)變化. 圖中藍線表示多年月平均海冰厚度、黑色豎棒表示標準差

Fig.9. Seasonal variation of sea ice thickness of the northern area, southern area and intersectional area over the period of 1979—2017. The blue lines are the multi- year monthly sea ice thickness and the vertical bars are the standard deviations

3.3.2 年際和月際變化特征

根據(jù)1979—2017年研究區(qū)域海冰厚度年際變化(圖10), 北部區(qū)域多年平均海冰厚度高于南部區(qū)域(1.26 m)和交匯區(qū)域(1.25 m), 為1.56 m, 并且北部和南部區(qū)域海冰厚度呈現(xiàn)顯著下降趨勢, 分別為–0.25 m·(10a)–1和–0.13 m·(10a)–1(以上趨勢通過99%的顯著性檢驗), 交匯區(qū)域趨勢并不顯著。39年間北部和南部區(qū)域海冰厚度最大值發(fā)生在1979年, 分別為2.27 m和1.83 m, 交匯區(qū)域海冰厚度最大值發(fā)生在1980年, 達到1.65 m; 3個區(qū)域最小值分別發(fā)生在2016年(0.89 m)、2016年(0.94 m)和1988年(0.97 m)。

圖10 1979—2017年北部區(qū)域、南部區(qū)域、交匯區(qū)域年平均海冰厚度變化(藍色柱狀), 圖中紅線表示線性回歸趨勢

Fig.10. Yearly sea ice thickness (blue bars) of the northern area, southern area and intersectional area over the period of 1979—2017. The red lines are the linear regression trends

根據(jù)1978年11月—2017年12月研究區(qū)域的月平均海冰厚度時間序列圖(圖11), 北部、南部、交匯區(qū)域海冰厚度最小值發(fā)生在2012年9月, 厚度為0.00 m; 最大值分別發(fā)生在1990年5月、1980年5月以及1998年4月, 海冰厚度分別達到3.15 m、2.62 m與3.54 m(圖11中紅點)。3個區(qū)域的海冰厚度的季節(jié)變化特征為春季高、夏秋低。根據(jù)表2, 北部區(qū)域和南部區(qū)域海冰厚度在四季都呈現(xiàn)顯著下降趨勢, 其中夏季下降趨勢最大[–0.34 m·(10a)–1, –0.17 m·(10a)–1], 其次為秋季[–0.27 m·(10a)–1, –0.14 m·(10a)–1]和冬季[–0.19 m·(10a)–1, –0.13 m·(10a)–1], 春季下降趨勢最小[–0.18 m·(10a)–1, –0.10 m·(10a)–1], 以上趨勢都通過99%的顯著性檢驗; 交匯區(qū)域在冬季、夏季和秋季呈現(xiàn)下降趨勢, 其中冬季和夏季下降趨勢并不顯著, 秋季下降趨勢為–0.11 m·(10a)–1(通過99%的顯著性檢驗), 春季呈現(xiàn)上升趨勢為0.08 m·(10a)–1, 但趨勢并不顯著, 不具有統(tǒng)計學(xué)意義。北部區(qū)域海冰厚度在四季的下降趨勢都大于南部區(qū)域, 下降趨勢最大月份為8月[–0.37 m·(10a)–1], 最小月份為2月[–0.18m·(10a)–1]; 南部區(qū)域下降趨勢最大月份為7月[–0.18 m·(10a)–1], 最小月份為4月[–0.09 m·(10a)–1], 以上趨勢都通過99%的顯著性檢驗。交匯區(qū)域海冰厚度在2—6月都呈現(xiàn)上升趨勢, 但趨勢并不顯著, 不具有統(tǒng)計學(xué)意義; 7月至次年1月呈現(xiàn)下降趨勢, 其中1月趨勢并不顯著, 10月下降趨勢最大為–0.13 m·(10a)–1(通過99%的顯著性檢驗)。

4 冰情與熱力學(xué)影響因素分析

4.1 西北航道南北航線冰情比較

從空間分布情況可以發(fā)現(xiàn), 冬春季南線北線都被高密集度海冰覆蓋。夏秋季北線的海冰分布呈現(xiàn)西高東低的分布特征, M’Clure Strait和Viscount Melville Sound區(qū)域海冰冰情影響著北線的開通。南線冰情較重的區(qū)域則是Coronation Gulf、Victoria Strait、Peel Sound和Barrow Strait, 整體呈現(xiàn)東西輕, 中間重的分布特征。

圖11 1978年11月—2017年12月北部區(qū)域、南部區(qū)域、交匯區(qū)域月平均海冰厚度(藍線). 圖中紅線表示線性回歸趨勢, 紅點表示海冰厚度最大值

Fig.11. Monthly sea ice thickness (blue lines) of the northern area, southern area and intersectional area over the period of November 1978—December 2017. The red lines are the linear regression trends and the red dots are the maximum of sea ice thickness

表2 1979—2017年海冰厚度變化趨勢

注: S在99%, 95%, 90%的置信水平上有統(tǒng)計學(xué)意義。

西北航道北線海冰較厚。5月整體海冰厚度大于2 m, 9月海冰厚度特征為西部大、東部小。在跨水道方向, M’Clure Strait中部和Viscount Melville Sound南岸較厚的海冰是影響西北航道開通的重要因素。南線海冰空間差異性較大, 5月海冰厚度呈現(xiàn)西低東高的分布特征, Victoria Strait西南部海冰存在大于3 m的厚冰, 9月海冰厚度分布東西低、中間較大, Victoria Strait南岸的較厚海冰是南線冰情險要處。

從長期變化趨勢來看, 相較于北線, 南線海冰覆蓋范圍以及海冰厚度較小, 且覆蓋范圍下降趨勢較為顯著, 而北線海冰厚度呈現(xiàn)了更為劇烈的下降趨勢。

4.2 溫度變化對西北航道海冰冰情的影響

4.2.1 西北航道表面大氣溫度(SAT)和海洋表層溫度(SST)的時間變化特征

隨著全球氣候變化, 表面大氣溫度和海洋表層溫度都持續(xù)上升, 圖12顯示的是1979—2017年的年平均表面大氣溫度、海洋表層溫度時間序列圖?？傮w來看, 西北航道區(qū)域表面大氣溫度、海洋表層溫度呈現(xiàn)上升趨勢。就表面大氣溫度而言, 南部區(qū)域氣溫較高, 多年平均氣溫為–11.57℃,并以0.68℃·(10a)–1(通過99%的顯著性檢驗)的趨勢升溫, 北部以及交匯區(qū)域氣溫較低, 多年平均氣溫分別達到–13.78℃和–13.33℃, 上升趨勢為0.58℃·(10a)–1和0.59℃·(10a)–1(通過99%的顯著性檢驗)。對于海洋表層溫度, 同樣也是南部區(qū)域具有較高的海表溫, 多年平均值為–0.87℃, 交匯區(qū)域為–1.28℃, 北部區(qū)域海表溫較低, 為–1.57℃,并且南部區(qū)域海表溫上升趨勢更為顯著, 為0.28℃·(10a)–1, 交匯區(qū)域為0.15℃·(10a)–1, 而北部區(qū)域趨勢較弱, 為0.07℃·(10a)–1(以上趨勢通過99%的顯著性檢驗)。表面大氣溫度和海洋表層溫度對西北航道區(qū)域海冰密集度均呈現(xiàn)顯著負相關(guān), 前者的相關(guān)系數(shù)在北部區(qū)域為–0.74, 南部區(qū)域為–0.84, 交匯區(qū)域為–0.63; 后者的相關(guān)系數(shù)在北部區(qū)域為–0.96, 南部區(qū)域為–0.87, 交匯區(qū)域為–0.85(以上相關(guān)系數(shù)通過99%的顯著性檢驗)。表面大氣溫度和海洋表層溫度對北部區(qū)域和南部區(qū)域的海冰厚度呈現(xiàn)顯著負相關(guān), 對交匯區(qū)域的海冰厚度并沒有顯著相關(guān)性。北部區(qū)域的海冰厚度與表面大氣溫度相關(guān)性略高, 為–0.72, 與海洋表層溫度相關(guān)性略低, 為–0.59(以上相關(guān)系數(shù)通過99%的顯著性檢驗)。南部區(qū)域海冰厚度與海洋表層溫度相關(guān)性為–0.75, 與表面大氣溫度相關(guān)性達到–0.70(以上相關(guān)系數(shù)通過99%的顯著性檢驗)。交匯區(qū)域的海冰厚度在研究中與表面大氣溫度和海洋表層溫度沒有顯著相關(guān), 這是由于交匯區(qū)域作為加拿大北極群島海冰最主要的輸出通道之一, 其海冰厚度會受到上游伊麗莎白女王群島和M’Clure Strait的厚冰輸入影響。

圖12 1979—2017年北部區(qū)域、南部區(qū)域、交匯區(qū)域年平均表面大氣溫度(紅線)、海洋表層溫度(藍線). 圖中紅色虛線表示表面大氣溫度線性回歸趨勢; 藍色虛線表示海洋表層溫度線性回歸趨勢

Fig.12. Yearly SAT (red solid lines), SST (blue solid lines) of the northern area, southern area and intersectional area over the period of 1979—2017. The red dash lines are the linear regress trends estimated from the yearly SAT, and the blue dash lines are the linear regress trends estimated from the yearly SST

4.2.2 冰情對表面大氣溫度和海洋表層溫度的響應(yīng)

由于1998年和2012年為西北航道區(qū)域的最輕冰情年份, 我們以該兩年為例, 計算了表面大氣溫度異常(SAT Anomaly, SATA)和海洋表層溫度異常(SST Anoamly, SSTA), 從而進一步分析對冰情的具體影響(圖13)。表面大氣溫度異常和海洋表層溫度異常的計算方法為該年數(shù)據(jù)減去1979—2017年多年平均值?？傮w上來看, 1998年表面大氣溫度明顯高于2012年, Amundsen Gulf、M’Clure Strait、Viscount Melville Sound在1998年氣溫增溫顯著。對于南部區(qū)域, 1998年Amundsen Gulf通道口處氣溫高于氣候態(tài)平均4℃, 2012年增溫幅度雖不如1998年, 但增溫范圍更廣, 延伸至Victoria Strait。對于北部區(qū)域, 1998年M’Clure Strait通道口處氣溫異常達到3.4℃, 并且增溫延伸至Viscount Melville Sound, 2012年僅M’Clure Strait和Viscount Melville Sound交匯處氣溫異常較高處達到2.8℃。交匯區(qū)域氣溫在1998年和2012年增溫并不顯著, 2012年Lancaster Sound氣溫增溫強于1998。從海洋表層溫度空間分布來看, 除了1998年M’Clure Strait通道口處海表溫異常達到2℃, 西北航道其他區(qū)域2012年海表溫異常顯著高于1998年。南部區(qū)域海表溫增溫最為顯著, 其中Coronation Gulf區(qū)域海表溫異常達到3.5℃。2012年北部區(qū)域和交匯區(qū)域增溫范圍更廣, 海表溫增溫高達1℃。隨著未來北極地區(qū)的氣候變化, 對異常溫度的預(yù)報將對航道的通航起到關(guān)鍵作用。

圖13 1998年、2012年加拿大北極群島區(qū)域西北航道表面大氣溫度異常、海洋表層溫度異常分布圖

Fig.13. Distribution of SATA, SSTA in the Northwest Passage of the Canadian Arctic Archipelago in 1998 and 2012

5 總結(jié)

在海冰密集度方面, 加拿大北極群島的西北航道區(qū)域春冬季被高密集度海冰覆蓋, 夏秋季海冰密集度呈現(xiàn)北高南低的分布特征。1979—2017年南北航線覆蓋范圍均呈現(xiàn)下降趨勢, 南部和交匯區(qū)域下降趨勢大于北部區(qū)域, 下降趨勢主要受到夏秋季海冰覆蓋范圍下降影響, 其中9月下降趨勢最大, 1998年后下降趨勢顯著, 北部和南部區(qū)域海冰覆蓋范圍最小值發(fā)生在2012年9月, 交匯區(qū)域最小值發(fā)生在1998年9月。

在海冰厚度方面, 西北航道海冰厚度整體呈現(xiàn)北部、東南高, 西南區(qū)域低的分布特征, 5月海冰厚度達到最大, 9月為最小值。北部和南部區(qū)域的海冰厚度在四季都呈現(xiàn)顯著的下降趨勢, 夏秋季下降趨勢大于冬春季, 其中北部區(qū)域趨勢大于南部區(qū)域, 交匯區(qū)域海冰厚度并沒有顯著的趨勢, 3個區(qū)域的海冰厚度最小值發(fā)生在2012年9月。結(jié)合海冰密集度和海冰厚度的空間分布和變化趨勢來看, 西北航道南線冰情優(yōu)于北線。

對于海冰密集度和厚度, 熱力學(xué)是影響其變化的關(guān)鍵因素。在研究區(qū)域, 表面大氣溫度和海洋表層溫度均呈現(xiàn)上升趨勢, 并與海冰密集度存在緊密聯(lián)系，呈現(xiàn)顯著負相關(guān); 對于海冰厚度, 北部區(qū)域和南部區(qū)域呈現(xiàn)顯著負相關(guān), 而交匯區(qū)域相關(guān)性不顯著, 這可能是由于來自上游大量海冰的輸入影響了厚度變化, 還需未來進一步結(jié)合動力模型進行研究。在海冰密集度低值年份1998年, 表面大氣溫度在Amundsen Gulf和M’Clure Strait區(qū)域的變化最為顯著, 同時, 海洋表層溫度在Amundsen Gulf和M’Clure Strait區(qū)域也有較強的增溫, Amundsen Gulf和M’Clure Strait是影響該年海冰冰情的關(guān)鍵區(qū)域。相較于1998年, 2012年西北航道區(qū)域的表面大氣溫度增溫弱, 而海洋表層溫度整體增溫顯著, 相較于北部區(qū)域和交匯區(qū)域, 南部區(qū)域的表面大氣溫度和海洋表層溫度增溫顯著, 影響2012年冰情的關(guān)鍵區(qū)域為Amundsen Gulf和Coronation Gulf。

1 ZHANG Y, CHEN C S, BEARDSLEY R C, et al. Seasonal and interannual variability of the Arctic sea ice: A comparison between AO-FVCOM and observations[J]. Journal of Geophysical Research: Oceans, 2016, 121(11): 8320-8350.

2 PARKINSON C L, COMISO J C. On the 2012 record low Arctic sea ice cover: Combined impact of preconditioning and an August storm[J]. Geophysical Research Letters, 2013, 40(7): 1356-1361.

3 LINDSAY R, SCHWEIGER A. Arctic sea ice thickness loss determined using subsurface, aircraft, and satellite observations[J]. The Cryosphere, 2015, 9(1): 269-283.

4 KWOK R, CUNNINGHAM G F, WENSNAHAN M, et al. Thinning and volume loss of the Arctic Ocean sea ice cover: 2003—2008[J]. Journal of Geophysical Research: Oceans, 2009, 114(C7): C07005.

5 PARKINSON C L, CAVALIERI D J. Arctic sea ice variability and trends, 1979—2006[J]. Journal of Geophysical Research: Oceans, 2008, 113(C7): C07003.

6 FALKINGHAM J C, CHAGNON R, MCCOURT S. Sea ice in the Canadian Arctic in the 21st century[C]//Proceedings of the International Conference on Port and Ocean Engineering Under Arctic Conditions. Ottawa, Canada: Port and Ocean Engineering under Arctic Conditions (POAC), 2001: 10.

7 蘇潔, 徐棟, 趙進平, 等. 北極加速變暖條件下西北航道的海冰分布變化特征[J]. 極地研究, 2010, 22(2): 104-124.

8 李春花, 李明, 趙杰臣, 等. 近年北極東北和西北航道開通狀況分析[J]. 海洋學(xué)報, 2014, 36(10): 33-47.

9 HOWELL S E L, TIVY A, YACKEL J J, et al. Application of a SeaWinds/QuikSCAT sea ice melt algorithm for assessing melt dynamics in the Canadian Arctic Archipelago[J]. Journal of Geophysical Research: Oceans, 2006, 111(C7): C07025.

10 HOWELL S E L, TIVY A, YACKEL J J, et al. Changing sea ice melt parameters in the Canadian Arctic Archipelago: Implications for the future presence of multiyear ice[J]. Journal of Geophysical Research: Oceans, 2008, 113(C9): C09030.

11 HOWELL S E L, DUGUAY C R, MARKUS T. Sea ice conditions and melt season duration variability within the Canadian Arctic Archipelago: 1979—2008[J]. Geophysical Research Letters, 2009, 36(10): L10502.

12 HOWELL S E L, WOHLLEBEN T, DABBOOR M, et al. Recent changes in the exchange of sea ice between the Arctic Ocean and the Canadian Arctic Archipelago[J]. Journal of Geophysical Research: Oceans, 2013, 118(7): 3595-3607.

13 MELLING H. Sea ice of the northern Canadian Arctic Archipelago[J]. Journal of Geophysical Research: Oceans, 2002, 107(C11): 2-1-2-21.

14 HAAS C, HOWELL S E L. Ice thickness in the Northwest Passage[J]. Geophysical Research Letters, 2015, 42(18): 7673-7680.

15 KURTZ N T, GALIN N, STUDINGER M. An improved CryoSat-2 sea ice freeboard retrieval algorithm through the use of waveform fitting[J]. The Cryosphere, 2014, 8(4): 1217-1237.

16 CHEN C S, GAO G P, ZHANG Y, et al. Circulation in the Arctic Ocean: Results from a high-resolution coupled ice-sea nested Global-FVCOM and Arctic-FVCOM system[J]. Progress in Oceanography, 2016, 141: 60-80.

17 GAO G P, CHEN C S, QI J H, et al. An unstructured-grid, finite-volume sea ice model: Development, validation, and application[J]. Journal of Geophysical Research: Oceans, 2011, 116(C8): C00D04.

LONG-TERM SPATIAL AND TEMPORAL VARIATIONS OF SEA ICE IN THE NORTHWEST PASSAGE OF THE CANADIAN ARCTIC ARCHIPELAGO

Shen Xinyi1, Zhang Yu1,2,3, Chen Changsheng4,1,3, Hu Song1,3

(1College of Marine Sciences, Shanghai Ocean University, Shanghai 201306, China;2Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China;3International Center for Marine Studies, Shanghai Ocean University, Shanghai 201306, China;4School for Marine Science and Technology, University of Massachusetts Dartmouth, New Bedford 02744, Massachusetts, USA)

Navigation through the Northwest Passage of the Canadian Arctic Archipelago is heavily influenced by sea ice conditions. The model domain used in this study was divided into the northern, southern and intermediate areas. We studied the spatial and temporal characteristics of sea ice in the Northwest Passage between November 1978 and December 2017 and evaluated the ice conditions of the northern and southern routes by using Bootstrap sea ice concentrations from the National Sea and Ice Data Center (NSIDC) and sea ice thicknesses from CryoSat-2, CS2SMOS and Arctic Ocean-Finite Volume Community Ocean Model (AO-FVCOM). Over the study period, sea ice concentration in the Northwest Passage was high in winter and spring, low in summer and fall, high in the north, and low in the south. Sea ice extent has decreased by –0.01×105km2/decade (–0.77%/decade) in the northern area, –0.09×105km2/decade (–3.50%/decade) in the southern area, and –0.12×104km2/decade (–2.81%/decade) in the intermediate area. Sea ice thickness has decreased by –0.25 m/decade in the northern area and –0.13 m/decade in the southern area. Between November 1978 and December 2017, sea ice concentration, extent and thickness along the southern route were lower than those along the northern route. Ice conditions along the southern route were safer for navigation. We analyzed atmospheric and oceanic thermodynamic parameters that influence sea ice concentration and thickness. Surface air temperature and sea surface temperature in the study area have increased and are significantly and negatively correlated with the sea ice concentration and thickness of the northern and southern areas. However, there are no significant correlations between temperatures and sea ice thickness in the intermediate area.

Canadian Arctic Archipelago, Northwest Passage, sea ice, SAT, SST

2020年3月收到來稿, 2020年4月收到修改稿

國家自然科學(xué)基金(41706210)、國家重點研發(fā)計劃(2019YFA0607000, 2016YFC1400903)、海洋環(huán)境安全保障國家重點研發(fā)計劃(2018YFC1406801)資助

沈心儀, 女, 1996年生。碩士研究生, 主要從事極地海冰數(shù)模與觀測相關(guān)研究。E-mail:xyshen9676@163.com

張瑜, E-mail: yuzhang@shou.edu.cn

10. 13679/j.jdyj.20200009