李娟，左軍成，譚偉，陳美香，趙雪

(1.國家海洋局 南海預(yù)報(bào)中心，廣東 廣州 510300；2.河海大學(xué) 港口海岸與近海工程學(xué)院，江蘇 南京 210098)

21世紀(jì)格陵蘭冰川融化速率對海平面變化的影響

李娟1，2*，左軍成2，譚偉2，陳美香2，趙雪1

(1.國家海洋局 南海預(yù)報(bào)中心，廣東 廣州 510300；2.河海大學(xué) 港口海岸與近海工程學(xué)院，江蘇 南京 210098)

本文利用大洋環(huán)流模式POP研究RCP4.5情景下21世紀(jì)格陵蘭冰川不同的融化速率對全球及區(qū)域海平面變化的影響。結(jié)果顯示：當(dāng)格陵蘭冰川的融化速率以每年1%增加時(shí)，全球大部分海域的動(dòng)力和比容海平面變化基本不變，主要是由于格陵蘭冰川在低速融化時(shí)并不會(huì)導(dǎo)致大西洋經(jīng)向翻轉(zhuǎn)流減弱。當(dāng)格陵蘭冰川的融化速率以每年3%和每年7%增加時(shí)，動(dòng)力海平面在北大西洋副極地、大西洋熱帶、南大西洋副熱帶和北冰洋海域呈現(xiàn)出顯著的上升趨勢，這是因?yàn)楦窳晏m冰川快速融化導(dǎo)致大量的淡水輸入附近海域，造成該上層海洋層化加強(qiáng)和深對流減弱，導(dǎo)致大西洋經(jīng)向翻轉(zhuǎn)流顯著減弱；與此同時(shí)，熱比容海平面在北冰洋、格陵蘭島南部海域和大西洋副熱帶海域顯著下降，而在熱帶大西洋和灣流海域明顯上升；此時(shí)鹽比容海平面的變化與熱比容海平面是反相的，這是由于大量的低溫低鹽水的輸入，造成北大西洋副極地海域變冷變淡、大西洋經(jīng)向翻轉(zhuǎn)流和熱鹽環(huán)流顯著減弱，引起了太平洋向北冰洋的熱通量和淡水通量減少，導(dǎo)致了北冰洋海水變冷變淡，同時(shí)熱帶大西洋滯留了更多的高溫高鹽水，隨著灣流被帶到北大西洋，北大西洋副極地海域低溫低鹽的海水，被風(fēng)生環(huán)流輸運(yùn)到副熱帶海域。

海平面變化；21世紀(jì)；格陵蘭冰川；RCP4.5；經(jīng)向翻轉(zhuǎn)流

1 引言

由于溫室氣體引起的全球變暖，格陵蘭冰川質(zhì)量減少的速率從20世紀(jì)90年代中期開始加速，在不到10年的時(shí)間里就增加了1倍[1]，到21世紀(jì)初格陵蘭冰川質(zhì)量的減少速率可能已達(dá)到了歷史最大值約為300 Gt/a[2—4]。模式結(jié)果顯示，如果溫室氣體的濃度保持在1 000×10-6，格陵蘭冰川將會(huì)在3 000年內(nèi)全部融化[5]，而冰川在第一個(gè)千年就會(huì)大部分融化，可使海平面上升3 m以上，這就意味著巨大的淡水將流入北大西洋，將會(huì)顯著影響北大西洋的經(jīng)向翻轉(zhuǎn)流。氣候變暖對于海洋和陸地冰的貢獻(xiàn)已經(jīng)包含了區(qū)域影響，其形式有海洋動(dòng)力導(dǎo)致的區(qū)域海水密度的變化和因陸地冰川融化導(dǎo)致的海水質(zhì)量的重新分布[6]。陸地冰川融化導(dǎo)致的淡水徑流和冰蓋的崩裂，能引起海平面的變化并會(huì)潛在的影響海洋的熱鹽環(huán)流，熱鹽環(huán)流的變化也會(huì)影響到海平面的變化[2，7]。目前一些模式結(jié)果顯示：如果氣候繼續(xù)變暖將會(huì)導(dǎo)致冰川快速融化，大量淡水輸入大洋，進(jìn)而引起經(jīng)向翻轉(zhuǎn)流急劇減弱[8—10]。

格陵蘭冰川逐漸減小，導(dǎo)致海水質(zhì)量逐年增加，從而影響全球及區(qū)域海平面變化。在21世紀(jì)，格陵蘭和南極冰川融化可能是對全球海平面變化貢獻(xiàn)最大的因子[11]。Gregory和Oerlemans[12]預(yù)估高山冰源和格陵蘭冰蓋的融化對海平面變化的貢獻(xiàn)在21世紀(jì)大約為20 cm。就目前的研究進(jìn)展來看，諸如冰川面積、高度以及形狀特征和質(zhì)量平衡等都有待更好地了解和深入研究。

陸地冰融化將會(huì)導(dǎo)致大量的淡水輸入到北大西洋，進(jìn)而引起大西洋經(jīng)向翻轉(zhuǎn)流的減弱[13—14]。大西洋經(jīng)向翻轉(zhuǎn)流(或者熱鹽環(huán)流)是全球大尺度海洋環(huán)流，把海洋上層暖而咸的水帶到北大西洋副極地海域，在那里失熱下沉，然后向南流到世界大洋的其他海域上升[15]。

對于未來氣候變化，利用模式研究陸地冰融化到目前為止，大多數(shù)的研究焦點(diǎn)放在了格陵蘭冰川的融冰水流入了北大西洋發(fā)生深層對流的海域。IPCC第四次評估報(bào)告指出，目前的耦合模式中沒有包含陸地冰模塊，主要的假設(shè)是未來氣候變暖會(huì)產(chǎn)生一個(gè)比較穩(wěn)定的層結(jié)海洋[16—17]。近年來的模式結(jié)果顯示，如果把格陵蘭冰川融化的融冰水理想的加入到模式模擬中去，如果格陵蘭冰川的融冰水足夠多，大西洋的經(jīng)向翻轉(zhuǎn)流將會(huì)減弱[11，18—19]。區(qū)域海平面偏離全球平均海平面因?yàn)閯?dòng)力海平面和靜力平衡海平面結(jié)構(gòu)的綜合影響，動(dòng)力海平面主要受海洋、大氣環(huán)流和溫度、鹽度的重新分配的影響；靜力平衡海平面結(jié)構(gòu)受重力、彈力和質(zhì)量重新分配旋轉(zhuǎn)效應(yīng)的影響。在目前格陵蘭冰蓋的融化速率下，局域海平面的變化主要受動(dòng)力海平面變化的影響；更高格陵蘭冰川融化速率下，動(dòng)力海平面趨勢在西北大西洋最強(qiáng)，當(dāng)融化冰相當(dāng)于20 cm等價(jià)海平面時(shí)，靜力平衡效應(yīng)占主導(dǎo)[20]。

IPCC為第五次報(bào)告開發(fā)了4種以穩(wěn)定性濃度為特征的新情景，它們分別為RCP8.5、RCP6、RCP4.5和RCP2.6。本文選擇的是中等排放情景RCP4.5。RCP4.5情景下，2100年預(yù)計(jì)升溫2.4～5.5℃，2100年后輻射強(qiáng)度穩(wěn)定在4.5 W/m2，大氣溫室氣體濃度穩(wěn)定在526×10-6CO2當(dāng)量，是一個(gè)比較穩(wěn)定的情景。本文關(guān)注的焦點(diǎn)是格陵蘭冰川不同的融化速率是如何通過影響大西洋經(jīng)向翻轉(zhuǎn)流來影響全球及區(qū)域海平面變化。在CMIP5模式結(jié)果的基礎(chǔ)上，本文利用三維海洋環(huán)流POP模式模擬RCP4.5情景下21世紀(jì)格陵蘭冰川不同的融化速率對全球及區(qū)域海平面變化的影響。

2 模式和數(shù)值試驗(yàn)介紹

2.1 模式介紹

海洋環(huán)流模式POP是美國Los Alamos國家實(shí)驗(yàn)室受能源部CHAMMP計(jì)劃資助開發(fā)的由早期的Bryan，Cox，Semtner和Chervin模型發(fā)展而來的一個(gè)三維海洋環(huán)流模式，2001年P(guān)OP被NCAR的共同氣候系統(tǒng)模式CCSM(Community Climate System Model)采用，作為其中的海洋模式。POP模式采用的坐標(biāo)系為球坐標(biāo)系，垂直方向采用z坐標(biāo)，模型控制方程采用了靜壓近似和Boussinesq近似。海表面采用的是自由表面邊界條件，該條件認(rèn)為海氣交界面可以自由發(fā)展，這就使得海表面高度成為可以預(yù)測的變量。

本文所使用的網(wǎng)格與CCSM4.0中海洋模式的網(wǎng)格是一致的，該網(wǎng)格在沿經(jīng)度的方向上有320個(gè)格點(diǎn)，在沿緯度方向上有384個(gè)點(diǎn)，南邊界在82°S的位置封閉，北邊界大約在80°N的位置封閉，網(wǎng)格的分辨率沿緯度的方向約為(1/2)°、沿經(jīng)度的方向大約為1°。從海表到深度為5 500 m的底層，垂直方向上共分為40層，最上層的間隔大約5 m，然后逐漸加大過渡到了深層的250 m左右。

2.2 數(shù)值試驗(yàn)

本文首先模擬20世紀(jì)的全球水文場作為21世紀(jì)模擬的初始場，模式的初始溫鹽場是CCSM4.0模擬的1850-1870年的氣候態(tài)月均數(shù)據(jù)，模式運(yùn)行200年穩(wěn)定后開始輸入逐月的變強(qiáng)迫場，模式從1850年運(yùn)行到2005年。

本文設(shè)置了4組試驗(yàn)，所有試驗(yàn)均開始于模式在2005年12月的restart場。試驗(yàn)1是控制試驗(yàn)(contralexp)，即沒有考慮格陵蘭冰川融化的影響；試驗(yàn)2～4是敏感性試驗(yàn)，即考慮了格陵蘭冰川不同的融化速率對全球及區(qū)域的海平面變化的影響。在敏感性試驗(yàn)2～4中，假設(shè)格陵蘭冰川融化只在北半球的夏半年融化，且融冰水全部均勻地流進(jìn)格陵蘭島南部的附近海域(58°～72°N，12°～60°W，見圖1)，其中冰川融化的影響是以等效鹽通量的形式加入到模式試驗(yàn)中的。格陵蘭冰川最初的平均融化速率為0.01 Sv (1 Sv=106m3/s，或者是1 mm/a全球等效海平面)，這個(gè)融化速度與目前觀測到的比較接近[2—3]。在敏感性試驗(yàn)2～4中，假設(shè)格陵蘭冰川融化在最初融化速率的基礎(chǔ)上增長率分別為每年1%(1%exp)、3%(3%exp)和7%(7%exp)，其中每年1%、3%的增長速率持續(xù)到2100年，每年7%試驗(yàn)中持續(xù)到2055年，然后保持不變到2100年。到2100年時(shí)，增長速率為1%的試驗(yàn)格陵蘭冰川的融化速率為0.03 Sv，增長速率為每年3%的試驗(yàn)中為0.19 Sv，每年7%試驗(yàn)為0.29 Sv，其中每年7%試驗(yàn)比較接近過去10年觀測得到的格陵蘭冰川質(zhì)量減少的增長率[2，4]。

本文所進(jìn)行的數(shù)值試驗(yàn)采用的表面強(qiáng)迫包括表面淡水通量、表面熱通量和表面風(fēng)應(yīng)力場，強(qiáng)迫場均來自CMIP5中CCSM4在RCP4.5情景下的模式結(jié)果。熱通量、淡水通量和風(fēng)應(yīng)力等強(qiáng)迫場是逐時(shí)的月均場。本文設(shè)置的4組試驗(yàn)，控制試驗(yàn)1作為敏感性試驗(yàn)的參考。

圖1 21世紀(jì)格陵蘭冰川融化的融冰水流入的海域Fig.1 The sea area of the input Greenland ice sheet melting water during the 21st century

3 21世紀(jì)格陵蘭冰川融化對區(qū)域動(dòng)力海平面變化的影響

海平面變化直接反映三維大尺度海洋環(huán)流變化。動(dòng)力海平面可簡單地通過水平流速和海平面梯度(或水平壓力梯度)之間的地轉(zhuǎn)平衡關(guān)系反映出來，流動(dòng)沿海平面等值線方向。區(qū)域海平面不僅受比容效應(yīng)的影響，還會(huì)受到海洋動(dòng)力的影響，即海平面變化還會(huì)受區(qū)域外海洋環(huán)流的控制。經(jīng)向翻轉(zhuǎn)流和風(fēng)生環(huán)流的變化引起局地海水堆積，從而導(dǎo)致局地海平面變化，而對全球而言，平均動(dòng)力海平面是個(gè)常數(shù)。

21世紀(jì)動(dòng)海平面上升最快的海域位于南大洋的43°S緯度帶海域，最大可達(dá)35 cm(見圖2)，而與其相鄰的南極繞極流海域動(dòng)海平面急劇下降，最大可為-50 cm，這與南大洋高緯度海域的海平面下降造成整個(gè)南大洋海面梯度增大和南極繞極流輸運(yùn)增強(qiáng)有關(guān)。北冰洋大部分海域、北大西洋的副極地環(huán)流海域、北太平洋北赤道流和副熱帶環(huán)流的西側(cè)海域動(dòng)海平面上升的比較快約為10～20 cm，這與大部分的模擬結(jié)果相似[21]，與該海域氣溫上升幅度大及淡水輸入有關(guān)。北大西洋副熱帶中心海域、南半球的副熱帶海域和北太平洋副極地海域動(dòng)海平面存在不同程度的下降，下降的幅度為1～20 cm(見圖2)，與Yin等[22]結(jié)果一致，這可能是由于該海域Sverdrup輸運(yùn)增強(qiáng)引起的動(dòng)力海平面下降。

研究表明，密度大的深層水在北大西洋生成，它是經(jīng)向翻轉(zhuǎn)流的一部分，因此，北大西洋的動(dòng)力海平面相對于北太平洋是較低的[23]。一般情況下，當(dāng)北大西洋深層水的生成減弱時(shí)，北大西洋的動(dòng)力海平面比其他海域上升的快[28]。格陵蘭冰川不同的融化速率會(huì)影響大西洋經(jīng)向翻轉(zhuǎn)流的變化(見圖3)，當(dāng)格陵蘭冰川的融化速率比較低時(shí)，經(jīng)向翻轉(zhuǎn)流不會(huì)減弱(見圖3a)，隨著格陵蘭冰川融化速率的增大，大西洋經(jīng)向翻轉(zhuǎn)流顯著減弱，特別是20°N以北1 500 m以淺的海域，最大可減弱0.5 Sv(見圖3b)，這與Hu等[11]的研究結(jié)果一致。

當(dāng)大西洋經(jīng)向翻轉(zhuǎn)流明顯減弱時(shí)，動(dòng)力海平面在北冰洋、北大西洋副極地、熱帶大西洋和南大西洋副熱帶海域顯著上升，而在印度洋、太平洋和南大西洋海域則是下降的(見圖4)。當(dāng)格陵蘭冰川以每年1%的加速度融化時(shí)，只有北大西洋副極地海域動(dòng)海平面呈現(xiàn)出微弱的上升趨勢(見圖4a)，這是由于冰川的融化速度較低時(shí)大西洋經(jīng)向翻轉(zhuǎn)流幾乎不會(huì)減弱(見圖3a)。隨著格陵蘭冰川融化速率的增大，動(dòng)力海平面在北大西洋副極地環(huán)流海域、大西洋熱帶、南大西洋副熱帶和北冰洋呈現(xiàn)出較快的上升趨勢(見圖4b、c)，到21世紀(jì)末，格陵蘭島南部海域和北美沿岸動(dòng)海平面最大可上升5 cm(見圖4c)；在歐洲沿岸、北大西洋副熱帶、太平洋、印度洋和南大洋呈現(xiàn)出不同程度的下降(見圖5b、c)，這是由于經(jīng)向翻轉(zhuǎn)流在3%exp和7%exp試驗(yàn)中顯著地減弱導(dǎo)致了灣流和北大西洋暖流減弱、拉布拉多寒流加強(qiáng)(見圖3)。

圖2 2081-2100年平均動(dòng)力海平面和1981-2000年平均動(dòng)力海平面的差值Fig.2 Dynamic sea level anomalies in the late 21st century relative to that in the late 20th century

圖3 敏感性試驗(yàn)2081-2100年大西洋經(jīng)向流函數(shù)與控制試驗(yàn)的差值Fig.3 Mean Atlantic meridional streamfunction anomaly between the 1%exp (a)，the 7%exp (b) and the contralexp at the end of 2081 to 2100

圖4 2081-2100年平均動(dòng)力海平面與控制試驗(yàn)的差值Fig.4 Dynamic sea level anomalies relative to that in the contralexp test of 2081 to 2100

4 格陵蘭冰川融化對全球平均比容海平面的影響

大氣中溫室氣體濃度的增加和格陵蘭冰川融化顯著影響全球和區(qū)域海平面變化。溫室氣體濃度的增加會(huì)導(dǎo)致海洋吸收較多的熱量，海洋中熱量的增加將會(huì)引起海水體積的增加，繼而影響熱比容海平面的上升。陸地冰融化、降水和徑流等還會(huì)引起海水的鹽度產(chǎn)生變化，這會(huì)導(dǎo)致鹽比容海平面的變化[19]。

21世紀(jì)全球平均比容海平面上升主要是由于熱比容效應(yīng)，鹽比容對比容海平面上升的貢獻(xiàn)相對比較小，這主要是由于海冰的融化、蒸發(fā)降水或者是海洋中鹽度的重新分布。敏感性試驗(yàn)中的比容海平面變化在21世紀(jì)與控制試驗(yàn)相似，在21世紀(jì)的前半個(gè)世紀(jì)，敏感性試驗(yàn)和控制試驗(yàn)的全球平均比容海平面幾乎是重合的(見圖5)，在21世紀(jì)中后期敏感性試驗(yàn)中的全球平均比容海平面變化相對于控制試驗(yàn)中的比容海平面有所上升，在21世紀(jì)末上升最大可達(dá)0.6 cm(見圖5)。

5 格陵蘭冰川融化對全球及區(qū)域比容海平面變化的影響

格陵蘭冰川融化產(chǎn)生的大量低溫低鹽水流入北大西洋，會(huì)導(dǎo)致經(jīng)向翻轉(zhuǎn)流的強(qiáng)度顯著減弱[13—14]。格陵蘭冰川融化會(huì)使北大西洋副極地海域的表層變冷和變淡，加上該海域的海水熱量不斷向大氣散熱，

因而密度變大，接著密度變大的海水下沉到深層，向南流動(dòng)，在大洋的其他海域上升。大量的模式結(jié)果顯示，經(jīng)向翻轉(zhuǎn)流的顯著減弱會(huì)引起北大西洋和附近海域變冷[11，14，24]，也會(huì)影響該海域的熱比容和鹽比容海平面的變化。

相對于20世紀(jì)的熱比容海平面，21世紀(jì)熱比容海平面在全球大部分海域呈現(xiàn)出顯著的上升趨勢，尤其是在20°N以南的大西洋比同緯度的太平洋上升的快(見圖6)，這是由于該海域2 000 m以淺的溫度在大西洋比太平洋增加的快[11]。熱比容海平面在北大西洋的20°～40°N海域呈現(xiàn)出顯著的下降趨勢(見圖6)，這是因?yàn)殡m然該海域的表層溫度是上升的，但是在深度為500～2000 m處急劇下降，最大下降了2.5℃[11]。熱比容海平面在南大洋表現(xiàn)出微弱的下降，這是由于雖然該海域表層是升溫的，而在500 m以深是降溫的[11]。

圖5 2000-2100年全球平均比容海平面變化Fig.5 Global mean steric sea level change of 2000 to 2100

圖6 控制試驗(yàn)2081-2100年平均熱比容海平面與1981-2000年的差值Fig.6 Thermosteric sea level anomalies in 2081 to 2100 relative to that in 1981 to 2000

格陵蘭冰川不同的融化速率對熱比容海平面變化的影響存在顯著的差異(見圖7)。當(dāng)格陵蘭冰川的融化速率較小時(shí)，全球大部分海域的熱比容海平面相對于控制試驗(yàn)來說變化不大，受影響的海域主要位于熱帶大西洋和北大西洋(見圖7a)，隨著格陵蘭冰川融化速率的增大，熱比容海平面在北冰洋、格陵蘭島南部海域和大西洋副熱帶海域表現(xiàn)出顯著的下降，在歐洲北部沿岸和格陵蘭島南部海域下降的最快，最大可達(dá)-2 cm，這是由于當(dāng)格陵蘭島冰川的融化速率比較大時(shí)，大量的低溫低鹽水輸入到格陵蘭南部海域，引起北大西洋副極地海域的深對流減弱，導(dǎo)致大西洋經(jīng)向翻轉(zhuǎn)流和熱鹽環(huán)流減弱(見圖3)，進(jìn)而導(dǎo)致太平洋向北冰洋的熱通量輸送減少，北冰洋海水變冷[11]；大西洋45°N附近海域熱比容海平面顯著上升，主要是受灣流的影響，大量的格陵蘭冰川融冰水輸入格陵蘭島南部海域?qū)е铝讼虮钡臒猁}環(huán)流減弱，從而使熱帶大西洋滯留了更多的熱量，隨著灣流被帶到北大西洋，引起了北大西洋45°N附近海域的次表層變暖(見圖8)。北大西洋副極地海域的熱比容海平面顯著下降的原因主要是由于經(jīng)向翻轉(zhuǎn)流的減弱導(dǎo)致了流向該海域的熱輸送減少(見圖8)；北大西洋副熱帶海域的熱比容海平面顯著下降主要是由于副極地的冷水被風(fēng)生環(huán)流輸運(yùn)到副熱帶海域(見圖8)。

圖7 2081-2100年平均熱比容海平面與1981-2000年的差Fig.7 Thermosteric sea level anomalies in 2081 to 2100 relative to that in 1981 to 2000

圖8 7%exp試驗(yàn)2081-2100年大西洋緯向平均溫度與控制試驗(yàn)的差值Fig.8 Zonal mean oceanic temperature anomaly between the 7%exp and the contralexp at the end of the 21st century in the Atlantic

RCP4.5情景下21世紀(jì)鹽比容海平面相對于20世紀(jì)末變化明顯(圖9)。鹽比容海平面在北大西洋、太平洋和東印度洋海域呈現(xiàn)出不同程度的上升趨勢(圖9)，這主要是由于該海域500～2 000 m的鹽度急劇減小[11]。熱帶大西洋、南大西洋和南大洋鹽比容海平面下降的比較快，北冰洋呈現(xiàn)出微弱的下降趨勢，這與緯向平均鹽度的分布有很好的對應(yīng)關(guān)系，即南極繞極流海域的鹽度在2 000 m以上是增加的，大西洋的20°S～35°N海域700 m以上鹽度明顯上升[11]。

圖9 控制試驗(yàn)2081-2100年平均鹽比容海平面與1981-2000年的差Fig.9 Halosteric sea level anomalies in 2081 to 2100 relative to that in 1981 to 2000

圖10 7%exp試驗(yàn)2081-2100年大西洋緯向平均鹽度與控制試驗(yàn)的差值Fig.10 Zonal mean oceanic salinity anomaly between the 7%exp and the contralexp of 2081 to 2100 in the Atlantic

圖11 2081-2100年平均鹽比容海平面與1981-2000年的差Fig.11 Halosteric sea level anomalies in 2081 to 2100 relative to that in 1981 to 2000

當(dāng)格陵蘭冰融化速率較小時(shí)，只會(huì)影響北大西洋的海平面變化，隨著融化速率的增大，繼而影響北冰洋、南大西洋、太平洋和印度洋海域的海平面變化(見圖11)。隨著格陵蘭冰川融化速率的增大，鹽比容海平面在北冰洋、格陵蘭島南部和大西洋副熱帶東側(cè)海域表現(xiàn)出顯著的上升，歐洲北部沿岸和格陵蘭島南部上升的最快，最大可達(dá)6 cm(見圖11c)，這是由于格陵蘭島冰川快速融化導(dǎo)致大量的淡水流入北冰洋和北大西洋(圖10)；大西洋45°N附近海域和熱帶大西洋鹽比容海平面顯著下降(見圖11c)，這是由于大西洋經(jīng)向翻轉(zhuǎn)流減弱導(dǎo)致了向北和向南的熱鹽環(huán)流減弱，從而使熱帶大西洋滯留了較多的高溫高鹽水，有一部分隨著灣流被帶到北大西洋副熱帶西側(cè)(圖10)。

6 結(jié)論

本文利用CCSM在RCP4.5情景下的模擬結(jié)果作為強(qiáng)迫場，模擬RCP4.5情景下21世紀(jì)格陵蘭冰川不同的融化速率對全球及區(qū)域海平面變化的影響。結(jié)果顯示當(dāng)格陵蘭冰川的融化速率以每年1%增加時(shí)，全球大部分海域的動(dòng)力海平面幾乎沒有變化，只有格陵蘭島南部海域和北美沿岸海域的動(dòng)力海平面有所上升，南大洋的南極繞極流和赤道太平洋海域呈現(xiàn)出微弱的下降；熱比容海平面和鹽比容海平面只在北大西洋發(fā)生微弱的變化，這是因?yàn)楦窳晏m冰川在低速融化時(shí)并不會(huì)導(dǎo)致經(jīng)向翻轉(zhuǎn)流減弱，即拉布拉多寒流、灣流和北大西洋暖流不會(huì)發(fā)生改變。

當(dāng)格陵蘭冰川的融化速率以每年3%和7%增加時(shí)，隨著格陵蘭冰川融化速率的增加，大西洋的經(jīng)向翻轉(zhuǎn)流顯著減弱，動(dòng)力海平面在北大西洋副極地環(huán)流海域、大西洋熱帶、南大西洋副熱帶和北冰洋呈現(xiàn)出顯著的上升趨勢，上升最快的海域是格陵蘭島南部海域和北美沿岸，這是因?yàn)楦窳晏m冰川融化大量的淡水輸入附近海域，造成該海域的上層海洋層化加強(qiáng)，深對流減弱導(dǎo)致大西洋經(jīng)向翻轉(zhuǎn)流顯著減弱，即拉布拉多寒流加強(qiáng)、灣流和北大西洋暖流減弱。

當(dāng)格陵蘭冰川的融化速率比較快時(shí)，熱比容海平面在北冰洋、格陵蘭島南部海域和大西洋副熱帶海域表現(xiàn)出顯著的下降，熱帶大西洋和灣流海域則上升明顯；與此同時(shí)鹽比容海平面的變化與熱比容海平面是反相的，這是由于格陵蘭島冰川的融化速率比較大時(shí)，大量的較冷較淡的水輸入到附近海域，造成北大西洋副極地海域的海水變冷變淡、大西洋經(jīng)向翻轉(zhuǎn)流和熱鹽環(huán)流顯著減弱，導(dǎo)致了太平洋向北冰洋的熱通量和淡水通量減少，北冰洋海水變冷變淡，熱帶大西洋滯留了更多的較暖較咸的海水，隨著灣流被帶到北大西洋，引起了灣流海域熱比容海平面上升和鹽比容海平面下降，北大西洋副極地海域較冷較淡的海水，被風(fēng)生環(huán)流輸運(yùn)到副熱帶海域。

[1] Rignot E，Kanagaratnam P. Changes in the velocity structure of the Greenland Ice Sheet[J]. Science，2006，311(5763): 986-990．

[2] Velicogna I. Increasing rates of ice mass loss from the Greenland and Antarctic ice sheets revealed by GRACE[J]. Geophysics Research Letter，2009，36: L19503．

[3] Rignot E I，Velicogna M R，van den Broeke，et al. Acceleration of the contribution of the Greenland and Antarctic ice sheets to sea level rise [J]. Geophysics Research Letter，2011，38: L05503．

[4] Li Juan，ZuoJuncheng，Cheng Meixiang，et al. Assessing the global averaged sea level budget over 2003-2010[J]. Acta Oceanologica Sinica，2013，32(10):16-23．

[5] Alley R B，Clark P U，Huybrechts P，et al. Ice-sheet and sea-level changes[J]. Science，2005，310(5747): 456-460．

[6] Katsman C A，Hazeleger W，Drijfhout S S，et al.Climate scenarios of sea level rise for the northeast Atlantic Ocean: a study including the effects of ocean dynamics and gravity changes induced by ice melt[J]. Climatic Change，2007，91(3/4): 351-374．

[7] Pritchard H D，Arthern R J，Vaughan D G，et al. Extensive dynamic thinning on the margins of the Greenland and Antarctic ice sheets[J].Nature，2009，461: 971-975．

[8] Huybrechts P，Janssens I，Poncin C，et al. The response of the Green-land Ice Sheet to climate changes in the twenty-first century by interactive coupling of an AOGCM with a thermomechnical ice-sheet model[C]//International symposium on Ice Cores and Climate，Kangerlussuaq (Greenland).2001．

[9] Jungclaus J H，Coauthors P. Will Greenland melting halt the thermohaline circulation？ [J]. Geophysics Research Letter，2006，33: L17708．

[10] Vizcaíno M，Mikolajewicz U，Gr?ger M，et al. Long-term ice sheet-climate interactions under anthropogenic greenhouse forcing simulated with a complex Earth System Model[J]. Climate Dynamics，2008，31(6): 665-690．

[11] Hu A，Meehl G A，Han W，et al. Effect of the potential melting of the Greenland Ice Sheet on the meridional overturning circulation and global climate in the future[J]. Deep Sea Research II，2011，58: 1914-1926．

[12] Gregory J M，Oerlemans J. Simulated future sea-level rise due to glacier melt based on regionally and seasonally resolved temperature changes[J]. Nature，1998，391: 474-476．

[13] Stouffer R J，Yin J，Gregpry J M，et al. Investigating the causes of the response of the thermohaline circulation to past and future climate changes[J]. Journal of Climate，2006，19(8): 1365-1387．

[14] Stouffer R J，Seidov D，Haupt B J. Climate response to external sources of freshwater: North Atlantic versus the Southern Ocean[J].Journal of Climate，2007，20(3): 436-448

[15] Stocker T F，Broecker W S. Observation and modeling of North Atlantic deep water formation and its variability: Introduction[J]. Journal of Geophysical Research，1994，99(12): 317．

[16] Schmittner A，Latif M，Schneider B. Model projections of the North Atlantic thermohaline circulation for the twenty-first century assessed by observations[J]. Geophysical Research Letter，2005，32: L23710．

[17] Meehl G A，Stocker T F. Global climate projections[R]//IPCC Fourth Assessment Report: Climate Change 2007. Cambridge: Cambridge University Press，2007: 747-845．

[18] Hu A，Meehl G A，Han W Q，et al.Transient response of the MOC and climate to potential melting of the Greenland Ice Sheet in the twenty-first century[J].Geophysical Research Letter，2009，36(10): L10707．

[19] Hu A，Meehl G A，Han W Q，et al. Influence of continental ice retreat on future global climate[J]. Journal of Climate，2013，26(3): 3087-3111．

[20] Kopp R E，Mitrovica J X，Griffies S M，et al. The impact of Greenland melt on local sea levels: a partially coupled analysis of dynamic and static equilibrium effects in idealized water-hosing experiments[J]. Climatic Change，2010，103: 619-625．

[21] Gregory J M，Church J A，Boer G J，et al. Comparison of results from several AOGCMs for global and regional sea-level change 1900-2100[J].Climate Dynamics，2001，18(3/4): 225-240．

[22] Yin J J，Schlesinger M E，Stouffer R J. Model projections of rapid sea-level rise on the northeast coast of the United States[J]. Nature Geoscience，2009，2(2): 262-266．

[23] Rio M H，Hernandez F. A mean dynamic topography computed over the world ocean from altimetry，in situ measurements，and a geoid model[J]. Journal of Geophysical Research，2004, 109: C12032.

[24] Sakai K，Peltier W R. Dansgaard Oeschger Oscillation in a coupled atmosphere——Ocean climate model [J]. Journal of Climate，1997，10(2): 949-970．

The influence of the melting rate of the Greenland ice sheet on the sea level variation for 21st century

Li Juan1，2，Zuo Juncheng2，Tan Wei2，Chen Meixiang2，ZhaoXue1

(1.SouthChinaSeaMarinePredictionCenter，StateOceanicAdministration，Guangzhou510300，China；2.CoastalandoffshoreEngineeringHohaiUniversity，Nanjing210098，China)

The influence of the different melting rate of Greenland ice sheet on sea level of the global and regional ocean are studied based on POP model under the RCP4.5 scenario for the 21st century. In the 21st century，the low melting rate of Greenland ice sheet has small influence on the global and regional sea level variations，which mainly due to the meridional overturning flow doesn’t weak under the low melting rate of Greenland ice sheet in the Atlantic. The high melting rate will have notably effect on the dynamic and steric sea level in the Atlantic，the Arctic and Antarctic Circumpolar Current (ACC). The dynamic sea level will significantly rise in the north Atlantic subpolar Ocean，tropical Atlantic Ocean，the south Atlantic subtropical and Arctic Ocean，which will rise fast in ocean of southern Greenland and along the coast of North America，that is due to notable weakened of the meridional overturning flow caused by strengthen stratification of the upper ocean and deep convection weakened because of a lot of fresh water input Greenland nearby ocean in the high rate of melting of Greenland ice sheet，namely the labrad or current will strengthen，the Gulf Stream and the north Atlantic current weakened. The thermosteric sea level will significantly decline in the Arctic ocean，the ocean of the southern Greenland and the Atlantic subtropical ocean，while which obviously rise the tropical Atlantic ocean and the Gulf Stream. Meanwhile the halosteric sea level change will be opposite to the thermosteric sea level. The steric sea level showed change due to the water is cooling and thin in the north Atlantic subpolar ocean，the weakened of the Atlantic meridional overturning flow and thermohaline circulation caused by input of lot of low temperature and low salt water in the 21st century．

sea level variation; 21st century; Greenland ice sheet; RCP4.5; meridional overturning flow

10.3969/j.issn.0253-4193.2015.07.003

2014-07-18；

2015-03-19。

國家自然科學(xué)基金(41276018)；海洋公益性項(xiàng)目(201005019)；海洋公益性項(xiàng)目(201305020-8)。

李娟(1982—)，女，山東省菏澤市人，博士，從事氣候與海平面變化和潮波動(dòng)力學(xué)研究。E-mail：twlj0419@126.com

P731.23

A

0253-4193(2015)07-0022-11

李娟，左軍成，譚偉，等. 21世紀(jì)格陵蘭冰川融化速率對海平面變化的影響[J]. 海洋學(xué)報(bào)，2015，37(7): 22-32，

Li Juan，Zuo Juncheng，Tan Wei，et al. The influence of the melting rate of the Greenland ice sheet on the sea level variation for 21st century[J]. Haiyang Xuebao，2015，37(7): 22-32，doi：10.3969/j.issn.0253-4193.2015.07.003