連 展，魏澤勛＊，方國洪，孫寶楠

（1.國家海洋局 第一海洋研究所，山東 青島266061；2.海洋環(huán)境科學(xué)和數(shù)值模擬國家海洋局重點實驗室，山東 青島266061）

近年來，由于溫室氣體排放量有增無減，全球氣溫變暖的態(tài)勢日趨明顯［1］。其帶來的一系列后果越來越受到人們的重視。全球變暖可對地球環(huán)境系統(tǒng)產(chǎn)生多種影響，其中海面高度升高是對人類影響相對較大的一項?？赡軐?dǎo)致海平面變化的因素有很多，如海水鹽度、溫度、質(zhì)量變化、氣象和地殼運動等都有可能對海平面帶來影響。

有研究表明，相對于海水熱比容變化，海水鹽比容變化導(dǎo)致的海平面變化較不顯著（Antonov et al.2005［2］）。一般情況下，因氣象和地殼運動導(dǎo)致的海平面變化幅度遠小于其它因素導(dǎo)致的海平面變化。對于陸地冰川開展的觀測工作較為有限，并且其資料可靠性較低［3］。因此，現(xiàn)今有關(guān)海平面對于氣候變化響應(yīng)的工作較關(guān)注海水熱比容變化導(dǎo)致的海平面變化?；跓崦浝淇s的基本原理，氣溫上升必然導(dǎo)致海水體積的膨脹，從而造成海平面的升高。政府間氣候?qū)I(yè)委員會（IPCC）綜合世界各地驗潮站數(shù)據(jù)與衛(wèi)星測量數(shù)據(jù)，1961—1993年間，因氣候變化造成的全球海平面平均的上升速度為每年1.0 mm［3］。此結(jié)果與蔡怡和李海［4］應(yīng)用不考慮兩極冰蓋融化增水的非Boussinesq近似模式的模擬結(jié)果、CABANES et al，［5］應(yīng)用熱力診斷得到的結(jié)果較為接近。這顯示了全球氣候變化導(dǎo)致的海水體積膨脹是全球海平面升高的主要因素。

在對海平面變化問題的研究中，可依據(jù)的研究方法有多種，主要分為分析觀測資料［4-8］以及數(shù)值模擬［9-10］。但衛(wèi)星觀測資料時間較短，驗潮站等定點觀測又受地殼運動等因素影響，而基于觀測資料驗證的數(shù)值模式可以較好地刻畫海平面變化情況的時空分布特征，同時在對未來情況預(yù)測上具有優(yōu)勢。以往對全球海平面變化的數(shù)值模擬中，應(yīng)用的模式大部分采用了Boussinesq近似，即在數(shù)值積分中考慮體積守恒而不是質(zhì)量守恒。這顯然無法體現(xiàn)海水因受熱而產(chǎn)生的體積膨脹效應(yīng)。少數(shù)非Boussinesq模式多僅進行理想模型試驗等［9］。蔡怡和李海［10］應(yīng)用改進的POP模式模擬了非Boussinesq近似下北太平洋的海平面狀況，但其模式覆蓋區(qū)域較小，結(jié)果未能反映全球海平面變化的整體情況，并且其應(yīng)用的改進型POP模式垂向坐標(biāo)為基于距離的Z坐標(biāo)，不利于精確模擬海水密度的變化［11］。本研究應(yīng)用壓力坐標(biāo)下的MITgcm模式（MIT General Circulation Model），建立了一個非Boussinesq近似下的全球海洋模式，模擬了近156 a來全球海平面的變化情況。并且根據(jù)對未來氣候變化趨勢的預(yù)測，對未來80 a內(nèi)全球海平面變化進行了模擬。

1 模式各項參數(shù)設(shè)置和優(yōu)化

MITgcm模式是由麻省理工學(xué)院開發(fā)的大氣-海洋通用環(huán)流模式。它的全稱為MIT General Circulation Model，是一個被設(shè)計用來研究大氣、海洋和氣候的數(shù)值模式［12］，并且得到了廣泛的應(yīng)用。在垂向坐標(biāo)的選擇上MITgcm提供Z坐標(biāo)與壓力坐標(biāo)兩種格式，其中應(yīng)用壓力坐標(biāo)可以實現(xiàn)對非Boussinesq近似流體的模擬［13］。

模式水平分辨率：1°×1°，水平覆蓋范圍：-80°S～80°N，-180°W～180°E（圖1）。

圖1 模式水平網(wǎng)格分布Fig.1 The horizontal grids of the model

垂向分為15層，每層厚度見表1。

表1 模式垂向分層Table 1 The vertically layering of water column for the model

模式地形數(shù)據(jù)選用的是ETOPO5數(shù)據(jù)集，溫鹽初始條件選用WOA（World Ocean Atlas）2005氣候態(tài)年平均數(shù)據(jù)，風(fēng)場數(shù)據(jù)選用ECWMF（European Centre for Medium-Range Weather Forecasts）氣候態(tài)年平均數(shù)據(jù)，SSS數(shù)據(jù)采用WOA2005氣候態(tài)年平均表層數(shù)據(jù)。SST數(shù)據(jù)的來源有兩種，分別為：1）ERSST（ex-tended reconstructed sea surface temperature data set）1855—2010年平均數(shù)據(jù)；2）IPCC第四次綜合報告給出的2010—2090年年平均數(shù)據(jù)。

通過與各類實測數(shù)據(jù)進行對比，開展了模式對于各類參數(shù)的敏感性研究，對關(guān)鍵參數(shù)進行了重點優(yōu)化，最終確定溫鹽垂向擴散系數(shù)相同且同為9.552 58×103Pa／s2≈9.0×10-5m／s2，模式海面溫鹽松弛系數(shù)和時間步長分別選為2 592 000 s和600 s。

2 模式的模擬流程

本研究建立的全球海平面變化數(shù)值模式共包含了3個主要事件過程，模擬流程見圖2，過程分別為：

過程1：模式由靜止啟動，SST強迫采用ERSST1855年年平均數(shù)據(jù)。為保證結(jié)果穩(wěn)定，模式共運行3 000 a。

過程2：以過程1為初始條件，以ERSST1855—2010年平均數(shù)據(jù)作為海表面強迫，對此段時期內(nèi)全球海平面高度進行模擬。

過程3：以過程2為初始條件，以IPCC第四次報告中2010—2090年年間年平均數(shù)據(jù)作為海表面強迫，對此段時期內(nèi)全球海平面高度進行模擬。

圖2 模擬時間流程圖Fig.2 Flow chart of simulation time

3 模擬結(jié)果驗證

圖3為模式過程一的結(jié)果中表層流場，從圖中我們可以看出黑潮、灣流和赤道流系等主要環(huán)流結(jié)構(gòu)在模擬結(jié)果中均已得到了較好的再現(xiàn)，可見本模式對于全球海洋環(huán)流結(jié)果具有一定的再現(xiàn)能力。

為了驗證模式對于全球海平面變化過程的模擬精度，我們將驗潮站水位觀測數(shù)據(jù)與模式過程二得到的水位數(shù)據(jù)進行了對比，驗潮站數(shù)據(jù)來源為PSMSL（Permanent Service for Mean Sea Level？）數(shù)據(jù)集。我們從中選取了觀測有效數(shù)據(jù)較多的且分別位于開闊海區(qū)（A 站點：22°8′W，21°19′N）和近岸海區(qū)（B站點：121°17′E，28°5′N）的兩個站點為代表站點（圖4）。圖5為二站點觀測數(shù)據(jù)與模擬結(jié)果的對比，可以看出本模式對于過程二期間的海平面變化情況具有一定的模擬能力。

圖3 模式表層流場Fig.3 The simulated surface current field

圖4 A和B站點站位圖Fig.4 The location of stations A and B

圖5 驗潮站觀測海平面變化與模擬結(jié)果對比Fig.5 Comparison of the sea-level changes observed at tidal station with those simulated

Church等在2011年應(yīng)用統(tǒng)計分析的方法，綜合采用驗潮站數(shù)據(jù)和衛(wèi)星高度計資料數(shù)據(jù)，對1880—2010年全球海平面變化情況進行了分析［14］，圖6為該統(tǒng)計結(jié)果與本模式模擬結(jié)果的對比。圖6中虛線為實測數(shù)據(jù)統(tǒng)計分析結(jié)果，實線為數(shù)值模擬結(jié)果，二者絕均差為19.724 6 mm，相關(guān)系數(shù)為0.992 8。

圖6 1880—2010年間全球海平面變化觀測結(jié)果與數(shù)值模擬結(jié)果對比Fig.6 Comparison of the global sea-level changes observed during 1880—2010 with those simulated

表2為模擬結(jié)果與我國國家海洋局發(fā)布的《2010年中國海平面公報》結(jié)果對比，可見本模式在中國近海海區(qū)具有較高的可信度。

表2 中國近海模擬結(jié)果與《2010年中國海平面公報》結(jié)果對比Table 2 Comparison of the simulated sea-level changes in the coastal waters of China with those published in“2010 Chinese Sea Level Bulletin”

4 模擬結(jié)果分析

通過對模擬結(jié)果的分析，得到2010年全球平均海平面相對于常年（即1975—1993年）上升43.8 mm，中國近海平均海平面相對于常年上升54.8 mm?？梢娔M時間段內(nèi)全球和中國近海均存在海平面上升的現(xiàn)象，且中國近海海平面上升速度略高于全球平均速度。

圖7為模擬得到的2010年全球海平面變化情況（相對于1975—1993即常年情況），全球海平面上升幅度較大的區(qū)域主要為：西北太平洋、北太平洋中部、南太平洋中部和北大西洋；東北太平洋和南極周邊海域海平面上升幅度較小，甚至存在海平面高度下降的區(qū)域。中國近海鄰近的西北太平洋地區(qū)存在海平面上升到情況，且上升速度由西向東逐漸增大。全球海平面變化幅度最大的地方位于格陵蘭島附近，可能與此地區(qū)劇烈的溫度變化有關(guān)［8］。圖8為中國近海和鄰近海域2010年海平面相對于常年上升幅度，從圖中可見渤海、黃海、東海和南海四個海區(qū)均存在海平面上升的趨勢，上升幅度略有差別。四個海區(qū)中海平面上升幅度最慢的為南海，黃海和東海則上升較快。中國近海海平面上升速度慢于鄰近的西北太平洋海區(qū)。海平面上升的幅度存在由西太平洋海平面上升高值區(qū)向中國近海逐漸減小的趨勢。在中國近海又存在由北至南逐漸減小的趨勢。

圖7 模擬得到的2010年全球海平面變化情況Fig.7 The simulated global sea-level changes in 2010

圖8 中國近海和鄰近海域2010年海平面相對于常年上升幅度（單位：m）Fig.8 The simulated range of the sea-level rising relatively to the perennial in the coastal waters of China in 2010（isoline unit：m）

5 未來海平面變化預(yù)測

為了對未來不同情境下全球和中國近海平均海平面高度變化情況進行可靠的預(yù)測，我們收集了IPCC綜合提供的全球各大主要氣候模式對未來的預(yù)測結(jié)果。所有模式均按照《IPCC排放情景特別報告》（SRES，2000）中所描述的情景。SRES情景分為探索可替代發(fā)展路徑的四個情景族（A1，A2，B1和B2），涉及一系列人口、經(jīng)濟和技術(shù)驅(qū)動力以及由此產(chǎn)生的溫室氣體排放。A1情景假定這樣一個世界：經(jīng)濟增長非?？?，全球人口數(shù)量峰值出現(xiàn)在本世紀中葉，新的和更高效的技術(shù)被迅速引進。A2情景描述了一個很不均衡的世界：人口快速增長、經(jīng)濟發(fā)展、技術(shù)進步緩慢。

綜合了各類模式對于未來不同情境的預(yù)測模擬之后，從中選取了對于2090年全球平均表面溫度升高最快預(yù)測結(jié)果（A2，最快情景）和最慢預(yù)測結(jié)果（A1，最慢情景），并且進而評估得到一般情況下全球平均表面溫度變化趨勢。圖9和圖10分別為最樂觀和最悲觀情境下2090年全球海平面相對于常年變化幅度，從圖中可以發(fā)現(xiàn)在最樂觀情境中，南大洋海區(qū)為海平面上升最強的區(qū)域，北大西洋海平面變化幅度也較劇烈，北太平洋東側(cè)相對于西側(cè)海平面上升幅度較大。在最悲觀情境中，北太平洋海平面上升較大區(qū)域位于海區(qū)中心偏西側(cè)。在最悲觀和最樂觀情境下，全球海平面平均上升幅度分別為和0.48和0.27 m，中國近海海平面平均上升幅度分別為0.42和0.18 m。綜合全球情況，北大西洋和南大洋仍為海平面變化劇烈海區(qū)，澳大利亞以南海區(qū)存在海平面下降的趨勢。

圖9 全球平均表面溫度升高最慢情境下2090年全球海平面高度相對于常年變化幅度Fig.9 The changing range of the global sea-level height relatively to the perennial in 2090 simulated under the condition of the slowest rising of global surface temperature

圖10 全球平均表面溫度升高最快情境下2090年全球海平面高度相對于常年變化幅度Fig.10 The changing range of the global sea-level height relatively to the perennial in 2090 simulated under the condition of the fastest rising of global surface temperature

表3為模擬得到的不同情境下2090年中國近海相對于常年海平面變化情況，可見至2090年渤海、黃海、東海和南海海平面變化趨勢各有不同，在上升最快情境下南海海區(qū)海平面上升幅度明顯小于其它海區(qū)，可能與此情景設(shè)計的氣溫變化情況有關(guān)。

6 結(jié) 論

通過對數(shù)值模擬結(jié)果的分析，我們可以得到如下結(jié)論：

1）全球海平面高度存在上升的現(xiàn)象。2010年全球平均海平面相對于常年（即1975—1993年）上升43.8 mm，中國近海平均海平面相對于常年上升54.8 mm。

2）2010年全球海平面變化情況，全球海平面上升幅度較大的區(qū)域主要為：西北太平洋、北太平洋中部、南太平洋中部和北大西洋。

3）渤海、黃海、東海和南海四個海區(qū)海平面上升幅度各不相同，黃東海海平面上升幅度略高，南海海平面上升幅度最小。中國近海海平面上升速度慢于鄰近的西北太平洋海區(qū)。

4）若全球溫度變化情況與預(yù)測相同，那么未來80 a后，全球平均海平面上升速度將遠遠超過現(xiàn)今水平。最極端情況下中國近海海平面上升幅度明顯大于全球平均水平。

在模式的建立和調(diào)整過程中，我們進行了多種敏感性實驗。實驗結(jié)果顯示模擬結(jié)果對于溫鹽垂向擴散系數(shù)有著較高的敏感性（圖在本文中未畫出），這顯示在考慮溫度長期變化趨勢的海面數(shù)值模式中，如何選取合理有效的溫鹽垂向擴散系數(shù)對于模式結(jié)果的可靠性至關(guān)重要。受計算時間的限制，本工作中選取的溫鹽垂向擴散系數(shù)為常數(shù)，下一步工作計劃采用更有實際代表性的參數(shù)化方案代替。同時，本研究的結(jié)果顯示，即使忽略海水質(zhì)量的變化，考慮海水熱脹冷縮的數(shù)值模式仍然可以較好地再現(xiàn)過去一段時期內(nèi)海平面的變化過程，并為未來變化趨勢提出有意義的預(yù)測。

另一方面，我們發(fā)現(xiàn)根據(jù)不同氣候變化情景預(yù)測得到的未來全球海平面變化的幅度均遠大于近期。根據(jù)黃瑞新［15］的研究結(jié)果，若采用Boussinesq近似對海洋進行模擬將會導(dǎo)致虛假的重力位能。在對未來海洋進行預(yù)測模擬中，這部分重力位能的差異將被進一步放大，這提示我們?nèi)舨捎肂oussinesq近似進行較長時間尺度的海洋數(shù)值模式對未來情景進行預(yù)測，Boussinesq近似造成的系統(tǒng)誤差將被越來越突出的表現(xiàn)出來，相對于低緯度地區(qū)，這部分誤差在北太平洋、北大西洋和南大洋海區(qū)格外顯著。

：

［1］ RICHARD A K.GLOBAL CHANGE：The Hottest Year，By a Hair［J］.Science，1998，279（5349）：315-316.

［2］ ANTONOV J I，LEVITUS S，BOYER T P.Steric sea level variations during 1957—1994：Importance of salinity［J］.J.Geophys Res，Oceans（1978—2012），2002，107（C12）：8013，doi：101029／2001JC000964.

［3］ IPCC.Climate Change 2007：the Fourth Assessment Report of the Intergovernmental Panel on Climate Change［M］.Cambridge，U K：Cambridge University Press，2007.

［4］ CAI Y，SONG X J.The analysis for sea surface temperature and sea surface height change in based on global warming［J］.Atomspheric science research and application，2004，26（1）：42-48.蔡怡，宋學(xué)家.全球增暖背景下太平洋海溫、海面氣溫和海面高度變化分析［J］.大氣科學(xué)研究與應(yīng)用，2004，26（1）：42-48.

［5］ CABANES C，CAZENAVE A，PROYOST C L.Sea level rise during past 40 years determined from satellite and in situ observations［J］.Science，2001，294（5543）：840-842.

［6］ MUNK W.Twentieth century sea level：An enigma［J］.Proc.Natl.Acad.Sci.U.S.A，2002，99：6550-6555.

［7］ DOUGLAS B C，KEARNEY M S，LEATHERMAN S P.Sea Level Rise：History and Consequences［M］.New York：Academic，2001.

［8］ CHURCH J A，GREGORY J M，HU YBRECHTYS P，et al.Climate Change 2001：The Scientific Basis［M］.Cambridge，U K：Cambridge University.Press，2001：639-694.

［9］ HUANG R X，JIN X Z，ZHANG X H.An oceanic general circulation model in pressure coordinates［J］.Advances in Atmospheric Sciences.2001，18（1）：1-22.

［10］ CAI Y，LI H.Simulation of sea level and heat content of Pacific Ocean based on global warming during last 40 years［J］.Journal of tropical oceanography，2006，25（4）：14-19.蔡怡，李海.全球變暖背景下最近40年太平洋海面高度和熱容量變化數(shù)值模擬［J］.熱帶海洋學(xué)報，2006，25（4）：14-19.

［11］ ZHANG X H，LI W，JIN X Z，et al.Development and Prospect of Ocean general circulation model in pressure coordinate［J］.Progress in Natural Science，2002，10（12）：1015-1020.張學(xué)洪，李薇，金向澤，等.壓力坐標(biāo)海洋環(huán)流模式的發(fā)展和前景［J］.自然科學(xué)進展，2002，10（12）：1015-1020.

［12］ MAROTZKE J，GIERING R，ZHANG K Q，et al.Construction of the adjoint MIT ocean general circulation model and application to Atlantic heat transport variability［J］.J.Geophysical Res.，1999，104（C12）：29529-29547.

［13］ DESZOEKE R A，SAMELSON R M.The duality between the Boussinesq and Non-Boussinesq hydrostatic equations of motion［J］.J.Phys.Oceanogr.，2002，32（8）：2194-2203.

［14］ CHURCH J A，WHITE N J.A 20th century acceleration in global sea-level rise［J］.Geophys.Res.Lett，2006，33（1）：L01602，doi：10.1029／2005GL024826.

［15］ HUANG R X.On the Balance of Energy in the Oceanic General Circulation［J］.Scientia atm ospherica sinica，1998，22（4）：562-574.黃瑞新.論大洋環(huán)流的能量平衡［J］.大氣科學(xué)，1998，22（4）：562-574.