劉伊格 柯長青 張杰

(1南京大學(xué)地理與海洋科學(xué)學(xué)院, 江蘇 南京 210023;2國家海洋局第一海洋研究所, 山東 青島 266061;3中國南海研究協(xié)同創(chuàng)新中心, 江蘇 南京 210023)

0 引言

海冰主要分布在地球的兩極[1],其存在阻礙了風(fēng)場對海水動量的直接輸入, 減緩了海洋和大氣間的能量交換, 影響海洋和大氣間的熱量、能量、動力的傳遞, 進(jìn)而影響大氣和海洋的循環(huán)[2]。南極海冰的變化不僅影響該地區(qū)局部的能量平衡,同時也對全球大氣環(huán)流和氣候變化產(chǎn)生影響。自20世紀(jì)70年代末有連續(xù)衛(wèi)星觀測記錄以來, 國內(nèi)外學(xué)者對南極海冰進(jìn)行了許多研究[3-8], 然而南極整體海冰變化與南極大陸邊緣海的海冰變化不同[7], 研究南極大陸邊緣海域的海冰變化有助于加深對南極整體海冰變化的認(rèn)識[9]。

羅斯海(Ross Sea)是南大洋五大海區(qū)中緯度最高的海域, 位于羅斯陸緣冰之北, 維多利亞地與瑪麗伯德地之間, 是南極洲西南極太平洋扇區(qū)內(nèi)的一個大海灣, 地理坐標(biāo)處于60°S—85°S、140°W—150°E(圖1)。覆蓋面積約為96萬km2, 全年覆蓋冰層, 多冰山分布。羅斯海是南極周邊的海區(qū)中海冰輸出最強(qiáng)的海區(qū), 所以該海區(qū)海冰變化會通過海洋、大氣等各種形式向低緯度傳播, 進(jìn)而影響同緯度和較低緯度的海洋和大氣變化[10]。

國內(nèi)外學(xué)者針對羅斯海海冰變化特征進(jìn)行過很多相關(guān)研究。鄔曉冬[10]利用ICESat衛(wèi)星數(shù)據(jù)集,首次給出了羅斯海海冰厚度的時空變化特征, 同時結(jié)合海冰密集度數(shù)據(jù), 分析了海冰面積的變化特征, 探討了兩者與大氣要素變化之間的聯(lián)系。劉帥斌等[1]在對羅斯海與普里茲灣海域海冰范圍進(jìn)行時間序列分析時, 發(fā)現(xiàn)羅斯海地區(qū)表現(xiàn)出“快速縮小、迅速擴(kuò)大”的特性, 2003—2014年海冰整體變化趨勢呈現(xiàn)(1.39±1.12)×104km2·a–1的擴(kuò)大趨勢, 而夏季的年際變化為減少趨勢。趙杰臣等[11]利用衛(wèi)星海冰密集度資料和船基海冰走航觀測數(shù)據(jù)分析了2012年12月—2013年3月羅斯海海冰密集度、厚度和浮冰尺寸等參數(shù)的時空變化特征。Cavalieri和Parkinson[7]研究了整個南極和五個分區(qū)海域1979—2006年海冰范圍和海冰面積變化, 發(fā)現(xiàn)羅斯海海冰范圍每年以(11.4±4.6)×103km2增長,海冰面積以(8.9±3.6)×103km2增長。Comiso等[12]研究了1978—2008年羅斯海海冰季節(jié)和年際變化特征, 發(fā)現(xiàn)羅斯海海冰范圍在整個南極海域中是增長最快的; 1992—2008年海冰漂移數(shù)據(jù)顯示, 整個羅斯冰架海冰輸出的凈增長率為30 000 km2·a–1。對羅斯海海冰變化的研究, 有助于更加深入了解這一海區(qū)的海冰變化特征, 探究其與更低緯度的大氣和海洋之間變化的關(guān)聯(lián), 加深對整個南極海冰變化的認(rèn)識。但是目前對羅斯海海冰的研究主要集中在海冰范圍、面積和厚度上[1,7,9-12], 對海冰運動速度的研究還較少。大氣和海洋對海冰的熱力和動力作用是海冰變化的主要原因, 熱力作用會直接影響海冰的凍結(jié)和融化[13-15], 動力作用主要影響海冰運動、破碎等, 進(jìn)而影響海冰分布以及熱力變化, 所以海冰動力作用在海冰總量的動態(tài)平衡過程中起著非常重要的作用[16-18]。

本文利用被動微波遙感衛(wèi)星資料分析了1979—2015年羅斯海海冰運動年際和季節(jié)變化特征以及海冰運動速度與海冰范圍之間的聯(lián)系, 最后對可能影響羅斯海海冰運動速度的因素進(jìn)行了探究。

圖1 羅斯海區(qū)域地理位置(箭頭表示海冰運動方向)Fig.1. Location of the Ross Sea (arrows indicate the direction of sea ice motion)

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

1.1 數(shù)據(jù)

美國冰雪數(shù)據(jù)中心(National Snow and Ice Data Center, NSIDC)提供1979年至今的每日和月平均南北極海冰運動矢量數(shù)據(jù), 該數(shù)據(jù)由AVHRR(Advanced Very High Resolution Radiometer)、SMMR(Scanning Multichannel Microwave Radiometer)和SSM/I(Special Sensor Microwave Imagers)傳感器獲得[19]。海冰密集度數(shù)據(jù)來自NSIDC發(fā)布的海冰索引數(shù)據(jù)集, 該數(shù)據(jù)集提供了1978年11月至今的南北極每日和月平均海冰密集度影像和數(shù)據(jù)值,數(shù)據(jù)集來自不同的衛(wèi)星微波輻射, 包括搭載在Nimbus-7衛(wèi)星上的SMMR (Scanning Multichannel Microwave Radiometer),美國國防氣象衛(wèi)星(Defense Meteorological Satellite Program,DMSP)F8、F11和F13搭載的SSM/I (Special Sensor Microwave Imagers)以及F17攜載的SSMIS(Special Sensor Microwave Imager Sounder)。兩種數(shù)據(jù)的分辨率都是25 km, 極地等面積方位投影。

海表面風(fēng)速月平均數(shù)據(jù)來自美國氣候中心發(fā)布的NCEP/NCAR再分析數(shù)據(jù)集, 該數(shù)據(jù)以NetCDF(Network Common Data Format)格式存儲1948年至今的資料, 網(wǎng)格點分辨率2.5°×2.5°。

1.2 方法

NSIDC通過最大互相關(guān)法獲得海冰運動速度, 即由相鄰時刻兩幅遙感影像中海冰的位移進(jìn)行計算得出。該數(shù)據(jù)計算海冰運動是基于被動微波連續(xù)影像上亮溫的變化, 使用Nimbus-7衛(wèi)星上的SMMR和DMSP衛(wèi)星上的SSM/I、SSMIS被動微波數(shù)據(jù)集得到的海冰密集度數(shù)據(jù)來表征海冰是否存在(海冰密集度大于15%被認(rèn)為存在)。海冰運動矢量數(shù)據(jù)產(chǎn)品最終以EASE-Grids網(wǎng)格形式發(fā)布, 每個網(wǎng)格包含分量u、v和第三個變量(像素值為0代表該處無海冰)。本文根據(jù)計算第三個變量為非0的網(wǎng)格點處海冰運動速度,不考慮海冰運動矢量的方向, 累加之后再平均即為月平均海冰運動速度。海冰范圍是海冰邊緣以內(nèi)所包圍的海域總面積, 本文通過NSIDC的海冰密集度數(shù)據(jù)集計算, 將密集度15%作為區(qū)分海冰和海域的分界線, 海冰密集度大于15%的網(wǎng)格點面積全部相加起來即獲得海冰范圍面積[20]。在研究區(qū)海域, 計算1979年1月—2015年12月海冰運動速度、海冰范圍以及風(fēng)速的月平均、季節(jié)平均和年平均, 在此基礎(chǔ)上利用最小二乘法對序列進(jìn)行回歸估計趨勢分析, 以確定數(shù)據(jù)的最佳線性擬合直線, 且都通過了顯著性為0.05的F檢驗。

2 結(jié)果與分析

2.1 海冰運動速度年際變化特征

1979—2015年羅斯海年平均海冰運動速度隨時間變化整體呈現(xiàn)增加趨勢(圖2), 其中1979—2010年保持速度波動增加的變化, 2010年海冰運動速度達(dá)到最大值234 m·h–1后開始呈下降趨勢, 之后幾年一直在144 m·h–1左右波動。海冰運動速度月變化曲線(圖3)則詳細(xì)反應(yīng)了每月的變化, 海冰運動速度總體上在波動中增加, 在2010年3月達(dá)到353 m·h–1的最大值, 之后運動速度明顯減緩, 與年平均變化(圖2)整體一致。每年海冰運動速度變化周期性并不明顯, 各年份海冰運動速度最大值的變化較大, 在83—353 m·h–1之間波動, 多出現(xiàn)在秋冬季。1979—2010年基本一直在增加, 2010年后下降穩(wěn)定在230 m·h–1左右。海冰運動速度最小值相比最大值變化波動較小, 在11—71 m·h–1波動, 大部分出現(xiàn)在每年夏季。

圖2 1979—2015年羅斯海海冰年平均運動速度變化Fig.2 Variation of yearly mean sea ice motion velocity in the Ross Sea from 1979 to 2015

用某月海冰運動數(shù)據(jù)減去所有數(shù)據(jù)中相同月份的平均值得到月距平變化數(shù)據(jù), 利用月距平變化可以詳細(xì)分析羅斯海海冰運動速度年際變化波動規(guī)律, 同時最大程度上減少季節(jié)變化帶來的影響[20]。總體來看也可以得到海冰運動速度在浮動范圍中保持增長的趨勢(圖4), 海冰運動速度平均每年增加0.3 m·h–1, 月距平變化大致為47 m·h–1。

圖3 1979—2015年羅斯海海冰運動速度月變化Fig.3 Variation of monthly sea ice motion velocity in the Ross Sea from 1979 to 2015

圖4 1979—2015年羅斯海海冰運動速度月距平變化Fig.4. Variation of monthly deviations about sea ice motion velocity in the Ross Sea from 1979 to 2015

2.2 海冰運動速度季節(jié)變化特征

本文按照南半球傳統(tǒng)季節(jié)劃分方法將南極分為春季(10—12月)、夏季(1—3月)、秋季(4—6月)和冬季(7—9月)進(jìn)行研究。海冰運動速度季節(jié)變化(圖5)表明羅斯海海冰運動速度增加趨勢最快的季節(jié)為秋季, 增加速率為4.4 m·h–1·a–1, 2007年秋季海冰運動速度更是達(dá)到了最大值287 m·h–1;其次是冬季, 增加速率為3.9 m·h–1·a–1, 1998年冬季海冰運動速度達(dá)到了次大值265 m·h–1, 但冬季是海冰平均運動速度最快的季節(jié); 春季為3.1 m·h–1·a–1,與秋冬季增加速率相差不大; 增加趨勢最小的是夏季, 同時也是四個季節(jié)中平均運動速度最慢的季節(jié), 夏季每個年份變化趨勢相對其他三個季節(jié)較為平緩, 但整體上還是略有增加, 增加速率為0.8 m·h–1·a–1。

圖5 1979—2015年羅斯海海冰運動速度季節(jié)變化Fig.5. Variation of seasonal mean sea ice motion velocity in the Ross Sea from 1979 to 2015

從月平均變化(圖6)中可以看到, 海冰運動速度在2月達(dá)到最低值后, 3、4月開始迅速加快, 5月運動速度有所下降, 之后月份速度緩慢增加, 8月達(dá)到峰值后又開始慢慢下降, 也剛好對應(yīng)羅斯海海冰平均運動速度冬季最快, 其次是秋季、春季,最后是夏季。

圖6 1979—2015年羅斯海海冰運動速度月平均變化Fig.6. Variation of monthly mean sea ice motion velocity in the Ross Sea from 1979 to 2015

2.3 海冰運動速度與海冰范圍之間的關(guān)系

上文結(jié)果得到1979—2015年羅斯海海冰平均運動速度具有加快的趨勢, 對海冰平均運動速度和海冰范圍的相關(guān)分析來探究海冰運動與海冰范圍之間可能存在的聯(lián)系。

圖7的回歸分析表明海冰運動速度與范圍均呈現(xiàn)上升趨勢, 海冰運動速度增長速率為0.3 m·h–1,而海冰范圍以每年1.11×103km2的速度擴(kuò)張。相關(guān)性研究發(fā)現(xiàn)兩者在0.01置信水平上顯著正相關(guān),相關(guān)系數(shù)為0.687。表明隨著羅斯海海冰運動速度的提高, 海冰范圍不斷地增加。羅斯海海冰運動模式的特點是沿海岸東部海冰流入羅斯海相對較弱, 西北向流出更強(qiáng)[12]。羅斯海海冰平均運動速度37年間具有加快的趨勢, 表明此區(qū)域常年盛行的西北向海冰環(huán)流速度逐年加快, 致使羅斯海區(qū)域海冰向低緯度移出, 海冰范圍逐漸擴(kuò)大, 從而對南極整體海冰分布產(chǎn)生影響。

海冰運動速度與海冰范圍月平均變化趨勢整體上一致(圖8), 都在2月達(dá)到最小值, 然后快速增加, 海冰運動速度在8月達(dá)到最大值而海冰范圍在9月達(dá)到最大值, 隨后相繼下降。不同的是, 海冰運動速度在5月達(dá)到極值后有所減慢,海冰范圍則是從2—9月一直保持?jǐn)U張趨勢。海冰范圍變化相對海冰運動速度有所滯后, 滯后1—2個月。

圖7 1979—2015年羅斯海年平均海冰運動速度與海冰范圍變化Fig.7 Relationship between yearly mean sea ice motion velocity and sea ice extent from 1979 to 2015

圖8 1979—2015年羅斯海海冰運動速度與海冰范圍月平均變化Fig.8. Variation of monthly mean sea ice motion velocity and sea ice extent in the Ross Sea from 1979 to 2015

3 海冰運動速度變化影響因素

3.1 風(fēng)場對海冰運動速度變化的影響

在0.01的置信水平下, 1979—2015年羅斯海區(qū)域年平均海冰運動速度與年平均風(fēng)速的相關(guān)系數(shù)是0.455, 兩者存在顯著的正相關(guān)關(guān)系(圖9)。1979年1月—2015年12月逐月海冰運動速度與逐月風(fēng)速的相關(guān)系數(shù)是0.492, 大致規(guī)律為風(fēng)速每加快1 km·h–1, 海冰運動速度增加8.4 m·h–1(圖10)。風(fēng)速加強(qiáng)是海冰運動速度加快的重要原因。

風(fēng)的作用直接給予海冰拖曳力, 使得海冰發(fā)生運動。Thorndike等[21]和Serreze等[22]發(fā)現(xiàn)風(fēng)速和海冰運動速度的關(guān)系是線性的, 只有在風(fēng)速很弱的情況下才變?yōu)榉蔷€性關(guān)系。國內(nèi)外很多文獻(xiàn)同樣證明風(fēng)場對海冰有持續(xù)性的影響[16,23], 與本文研究結(jié)果相一致。

圖9 1979—2015年羅斯海海冰運動速度與風(fēng)速變化對比Fig.9. The relationship between sea ice motion velocity and wind speed in the Ross Sea from 1979 to 2015

圖10 1979—2015年逐月海冰運動速度與風(fēng)速的相關(guān)關(guān)系Fig.10. The relationship between the monthly sea ice motion velocity and wind speed from 1979 to 2015

3.2 其他影響因素

除了風(fēng)場對海冰運動作用顯著外, 其他海洋環(huán)境因素也會對其產(chǎn)生影響。Comiso等[12]研究了1992—2008年羅斯海區(qū)域海冰運動情況, 認(rèn)為海冰運動的強(qiáng)度和空間位置分布取決于大氣中一直存在的低壓模式, 即氣壓場會影響海冰運動。Spreen等[23]認(rèn)為海冰運動在較小程度上會受到洋流場的影響, 同時也受到大氣、海洋阻力系數(shù)以及海冰內(nèi)部之間相互作用擠壓的影響。

4 結(jié)論和討論

使用1979—2015年NSIDC海冰運動和海冰密集度數(shù)據(jù), 分析了羅斯海海冰運動速度的年際和季節(jié)變化特征及其與海冰范圍可能存在的聯(lián)系,最后用海表面風(fēng)速數(shù)據(jù)探討海冰運動與風(fēng)場的關(guān)系, 得到如下結(jié)論。

1. 1979—2015年羅斯海海冰運動速度總體呈現(xiàn)加快趨勢。海冰運動速度增加趨勢最快的季節(jié)為秋季, 其次是冬季、春季, 增加趨勢最小的為夏季; 冬季海冰平均運動速度最大, 接著是秋季、春季和夏季。

2. 回歸分析表明37年間海冰運動速度與范圍均呈現(xiàn)上升趨勢, 海冰范圍月平均變化滯后海冰運動速度1—2個月。其中海冰運動速度增長速率為0.3 m·h–1, 而海冰范圍以每年1.11×103km2的速度擴(kuò)張。相關(guān)性研究發(fā)現(xiàn)兩者在0.01水平上顯著正相關(guān), 表明隨著羅斯海海冰運動速度的加快, 西北向海冰環(huán)流逐年增強(qiáng), 海冰向低緯度運動, 致使羅斯海海冰范圍不斷地增加, 進(jìn)而影響南極海冰整體分布。

3. 羅斯海年平均海冰運動速度與風(fēng)速、逐月海冰運動速度與風(fēng)速都呈現(xiàn)顯著正相關(guān)關(guān)系, 風(fēng)速對海冰的運動發(fā)揮著重要作用, 除此之外其他海洋環(huán)境因素也會影響海冰運動。

4. 海冰運動是一個氣-冰-海相互作用的復(fù)雜過程, 本文僅對羅斯海海冰運動時間序列演變進(jìn)行了初步研究, 空間地理位置的差異也會導(dǎo)致不同區(qū)域海冰運動的不同, 因此海冰運動空間分布的變化將是下一步研究方向之一。另外, 表征羅斯海海冰變化的因子有很多, 除了海冰運動、海冰范圍之外, 還有海冰厚度、海冰體積等, 這些因子對羅斯海海冰的變化分析也很重要, 在今后也有待進(jìn)一步研究。

1 劉帥斌, 周春霞, 王澤民. 羅斯海和普里茲灣海域海冰范圍變化對比分析[J]. 極地研究, 2016, 28(2): 228—234.

2 柯長青, 彭海濤, 孫波, 等. 2002年—2011年北極海冰時空變化分析[J]. 遙感學(xué)報, 2013, 17(2): 452—466.

3 PARKINSON C L. Global sea ice coverage from satellite data: annual cycle and 35-Yr trends[J]. Journal of Climate, 2014, 27(24): 9377—9382.

4 GLOERSEN P, CAMPBELL W J, CAVALIERI D J, et al. Satellite passive microwave observations and analysis of Arctic and Antarctic sea ice, 1978—1987[J]. Annals of Glaciology, 1993, 17: 149—154.

5 解思梅, 魏立新, 郝春江, 等. 南極海冰和陸架冰的變化特征[J]. 海洋學(xué)報, 2003, 25(3): 32—46.

6 卞林根, 林學(xué)椿. 近30年南極海冰的變化特征[J]. 極地研究, 2005, 17(4): 233—244.

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

8 TURNER J, COMISO J C, MARSHALL G J, et al. Non-annular atmospheric circulation change induced by stratospheric ozone depletion and its role in the recent increase of Antarctic sea ice extent[J]. Geophysical Research Letters, 2009, 36(8): L08502.

9 牟龍江, 趙進(jìn)平. 格陵蘭海海冰外緣線變化特征分析[J]. 地球科學(xué)進(jìn)展, 2013, 28(6): 709—717.

10 鄔曉冬. 南極羅斯海海冰厚度和面積變化研究[D]. 青島: 中國海洋大學(xué), 2012.

11 趙杰臣, 張林, 田忠翔, 等. 南極羅斯海2012年夏季海冰特征分析[J]. 極地研究, 2014, 26(3): 342—351.

12 COMISO J C, KWOK R, MARTIN S, et al. Variability and trends in sea ice extent and ice production in the Ross Sea[J]. Journal of Geophysical Research, 2011, 116(C4): C04021.

13 JACOB T, WAHR J, PFEFFER W T, et al. Recent contributions of glaciers and ice caps to sea level rise[J]. Nature, 2012, 482(7386):514—518, doi: 10.1038/nature10847.

14 CHAPMAN W L, WALSH J E. Recent variations of sea ice and air temperature in high latitudes[J]. Bulletin of the American Meteorological Society, 1993, 74(1): 33—48.

15 EBERT E E, CURRY J A. An intermediate one-dimensional thermodynamic sea ice model for investigating ice-atmosphere interactions[J]. Journal of Geophysical Research, 1993, 98(C6): 10085—10109.

16 田忠翔. 北極海冰運動變化特征分析[D]. 北京: 國家海洋環(huán)境預(yù)報研究中心, 2012.

17 左正道, 高郭平, 程靈巧, 等. 1979—2012年北極海冰運動學(xué)特征初步分析[J]. 海洋學(xué)報, 2016, 38(5): 57—69.

18 田忠翔, 李春花, 張林, 等. 氣候態(tài)下的北極海冰運動特征[J]. 海洋預(yù)報, 2012, 29(6): 66—73.

19 TSCHUDI M, FOWLER C, MASLANIK J, et al. Polar pathfinder daily 25 km EASE-grid sea ice motion vectors, version 3[R]. Boulder,Colorado, USA: NASA National Snow and Ice Data Center Distributed Active Archive Center, 2016, doi: 10.5067/O57VAIT2AYYY.

20 沈校熠, 柯長青, 張杰. 2002—2011年南極海冰變化的遙感分析[J]. 極地研究, 2017, 29(1): 163—171.

21 THORNDIKE A S, COLONY R. Sea ice motion in response to geostrophic winds[J]. Journal of Geophysical Research, 1982, 87(C8):5845—5852.

22 SERREZE M C, WALSH J E, CHAPIN III F S, et al. Observational evidence of recent change in the northern high-latitude environment[J]. Climatic Change, 2000, 46(1—2): 159—207.

23 SPREEN G, KWOK R, MENEMENLIS D. Trends in Arctic sea ice drift and role of wind forcing: 1992—2009[J]. Geophysical Research Letters, 2011, 38(9): L19501, doi: 10.1029/2011GL048970.