基于多源數據的全球比容海平面變化研究

李洪超，文漢江，師軍良

(1.黃河水利職業技術學院，河南 開封 454003；2.中國測繪科學研究院，北京 100830)

全球約有10%的人口分布在海拔較低的沿海地帶，海平面的上升將對這些地區的居民生存產生極大的影響。在過去幾十年間，由氣候變暖等因素導致全球海平面持續上升。IPCC(Intergovernmental Panel on Climate Change)第四次評估報告表明：20世紀全球海平面以1.7 mm/a的速度上升。

全球海平面變化主要有兩方面因素構成：一是海水溫度和鹽度變化導致的海水體積變化，這部分稱為比容海平面變化；一是陸地冰川、極地冰蓋等消融導致的海水質量變化。以前由于監測手段的限制，監測海平面變化只能依靠驗潮站進行，精度低且成本較高。衛星測高、衛星重力和海洋浮標技術的發展為監測海平面變化提供新的方式。當前，獲取比容海平面變化的方式主要有兩種：①聯合衛星測高和衛星重力數據。利用衛星測高技術可以獲取總體海平面變化，利用衛星重力數據可以獲取海水質量變化引起的海平面變化，兩者差值即為比容海平面變化。②利用海洋溫度和鹽度的剖面數據直接計算得出比容海平面變化。

1 數 據

1.1 衛星測高數據

衛星數據選取法國海洋衛星數據處理中心AVISO提供的海面高異常數據產品，該數據融合了Topex/Poseidon、Envisat、ERS-1、ERS-2、Jason-1、Jason-2等多顆測高衛星，較單顆衛星數據相比空間分辨率和時間分辨率具有較大的提高。該海面高異常數據的空間分辨率為0.25°×0.25°，時間分辨率為7 d。時間跨度為2005-01—2009-12共60個月數據。

1.2 衛星重力數據

利用GRACE衛星重力數據反演海水質量變化，GRACE數據采用美國德克薩斯空間研究中心(CSR)發布的GRACE RL05 Level-2版本重力場球諧系數模型，RL05版本數據對RL04版本數據中個別錯誤的地方進行了替換如2004-07、2004-10、2005-03和2006-02的數據，RL05版本數據與RL04版本數據相比噪聲水平明顯降低[6]， C20項的數據質量也有提高。RL05版本數據不僅提供60階還提供96階的重力場球諧系數模型。采用60階的重力場球諧系數模型，時間跨度為2005-01—2009-12共60個月的數據。

1.3 海洋數據

作為全球海洋觀測業務系統計劃(Global Ocean Observing System，GOOS)中的一個針對深海區溫鹽結構觀測的子計劃(Array for Real-time Geostrophic Oceanography，Argo)計劃目前在全球海洋布設超過3 900個Argo浮標，每年約可提供10萬條0～2 000 m海水剖面溫、鹽數據。我國于2002年加入Argo計劃，截止目前已經在我國周邊海域布設374個Argo浮標。采用Argo實時資料中心(http://www.argo.org.cn)提供的Argo浮標溫、鹽格網化產品，其分辨率為1°×1°；時間間隔為2005-01—2009-12共60個月的數據。

2 比容海平面計算方法

2.1 海水溫、鹽數據計算比容海平面的方法

海水任意位置處的比容海平面變化可沿剖面方向對海水密度變化積分[7]得到

(1)

2.2 聯合衛星測高、衛星重力數據計算比容海平面變化

2.2.1 衛星重力數據反演海水質量變化

由時變地球重力場球諧系數殘差(ΔClm(t),ΔSlm(t))計算地球表面質量變化(以等效水高表示)的算式[7]為

WlPlm(sinφ){ΔClm(t)cosmλ+ΔSlm(t)cosmλ}.

(2)

式中:ρE為地球的平均密度(5517 kg·m-2)；ρW為純水的密度(1 000 kg·m-2)；kl為負荷勒夫數(load Lover numbers)；Wl為空間平滑函數。

在利用GRACE重力場模型球諧系數計算海水質量變化時，對數據進行以下處理[7]：首先，采用Chen等[10]提出的地心改正模型對一階項球諧系數進行改正；用衛星激光測距(SLR)得到的C20項代替GRACE球諧系數的C20項；將GAC文件的球諧系數添加到GSM文件的球諧系數中，以恢復AOD1B模型影響；從GSM文件的重力場球諧系數中扣除重力場球諧系數的平均值，獲得重力場球諧系數的殘差；對重力場球諧系數殘差進行去相關誤差濾波處理，消除高階項系統誤差；對重力場球諧系數殘差進行空間平均(平滑半徑為500 km)得到海水質量變化。利用GRACE時變重力場模型反演海水質量變化時，冰川均衡調整(Glacial Isostatic Adjustment，GIA)模型的選擇是一個重要的影響因素，本文采用Paulson2007[10]模型。

2.2.2 衛星測高數據的處理

首先應將空間分辨率為0.25°×0.25°、時間分辨率為7 d衛星測高數據加權平均以獲取空間分辨率為1°×1°、時間分辨率為1月的海面高異常數據。

同時，為與GRACE的處理過程保持一致性，在聯合GRACE和衛星測高數據計算比容海平面時，應將對衛星測高數據進行同GRACE數據相同的平滑處理。本文采用與Chambers(2006)相同的處理方法：

(3)

式中：Ω表示對地球整個球面進行積分，利用式(3)計算時要求數據在全球范圍內是連續的，因此將陸地和沒有數據的海洋點的值設為零。

2)將上步得到的球諧系數殘差代入式(4)，并進行同GRACE數據一樣的平滑處理，則

(4)

式中：Δηalti(φ,λ,t)表示經過平滑處理后的海面高異常網格數據值。

3 比容海平面變化分析

利用Argo溫、鹽數據和聯合衛星測高、GRACE衛星數據兩種相對獨立的方法分別計算了全球66°S～66°N范圍內2005—2009年共5年的月比容海平面變化，如圖1所示(限于篇幅，此處僅列出2006年3月和9月的比容海平面變化)。圖1中可以看出全球比容海平面變化基本在-8～8 cm，局部海域比容海平面變化幅度較大如太平洋的中南部、西北部和中東部區域。比較兩種方法得到的比容海平面變化可以發現二者的整體變化趨勢較為一致，但局部存在著較大差異。GRACE時變重力場數據僅能探測出2～3 cm等效水高的海平面變化[3]以及數據處理過程中引入的各種誤差都可能使聯合衛星測高數據得到的比容海平面變化同利用Argo數據所得的比容海平面存在著空間分布上的差異。通過比較2006-03和2006-09的比容海平面變化的空間分布特點，可以發現2006-03的比容海平面北半球較5年的平均比容海平面為負值，南半球則為正值；2006-09則是北半球為正值，南半球為負值。3月北半球海水溫度較平均值低，使得比容海平面較平均值為負，南半球情況相反。

圖1 全球比容海平面變化

將兩種方法計算得到的空間分布的月比容海平面變化，按緯度進行全球加權平均可以得到各月的全球平均比容海平面變化，如圖2所示。從圖中可以看出比容海平面表現出強烈的年周期變化，且利用兩種方法得到的全球平均比容海平面變化的周期吻合較好，極大值一般在2月，極小值一般在8月；二者具有較強的相關性，相關系數為0.86。對全球平均比容海平面變化進行諧波分析可得，基于Argo數據得到的全球平均比容海平面變化振幅為4.3 mm，聯合衛星測高和GRACE得到的全球平均比容海平面變化振幅為10.5 mm。

圖2 全球平均月比容海平面變化

對全球月比容海平面變化數據進行時間序列分析，可以得出利用兩種不同方法獲得的比容海平面長期變化趨勢，如圖3所示。圖中可以看出全球大部分海域的比容海平面呈上升狀態，局部海域呈下降趨勢如太平洋中部。比較圖3(a)、圖3(b)兩幅圖可以發現二者的整體變化趨勢較為一致，但聯合衛星測高和GRACE圖3(a)得到的比容海平面的年變化速率明顯大于基于Argo數據圖3(b)計算得到的比容海平面年變化速率。按緯度進行全球加權平均可以得到聯合衛星測高和GRACE得到的比容海平面全球平均變化速率是 1.63 mm/a；基于Argo數據的比容海平面全球平均變化速率是0.32 mm/a。造成二者差別的部分原因可能是由GRACE所能探測到的海水質量變化的精度(僅能探測到2～3 cm等效水高的海水質量變化)、數據處理誤差等引起的，具體原因仍需做進一步的研究。

圖3 2005—2009年全球比容海平面長期變化趨勢

4 結束語

1)聯合衛星測高、GRACE和Argo數據，通過兩種獨立的方法計算得到的比容海平面變化具有較好的一致性，均呈現明顯的年周期變化；驗證了聯合衛星測高和GRACE重力衛星計算比容海平面變化的可行性。

2)聯合衛星測高、GRACE數據得到的比容海平面變化振幅(10.5 mm)較Argo數據得到的比容海平面變化振幅(4.3 mm)大，具體原因仍需做進一步研究。

3)兩種方法獲得的比容海平面變化整體趨勢較一致，但局部存在較大差異。

參考文獻：

[1] 文漢江，章傳銀．由ERS-2和TOPEX衛星測高數據推算的海面高異常的主成分分析[J]．武漢大學學報(信息科學版)，2006，31(3)：221-223．

[2] 文漢江，李洪超，蔡艷輝，等．聯合Argo浮標、衛星測高和GRACE數據的海平面變化研究[J]．測繪學報，2012，41(5):696-702．

[3] 鐘玉龍，鐘敏，馮偉．近十年全球平均海平面變化成因的衛星重力監測研究以及與ENSO現象的相關分析[J]．地球物理學報，2016,31(2):643-648．

[4] 陳美香，常曼，張雯皓，等．全球比容海平面低頻變化特征研究[J]．海洋科學進展，2016，34(2) :162-174．

[5] CHAMBERS D P.Observing seasonal steric sea level variations with GRACE and satellite altimetry[J]. JOUTNAL OF GEOPHYSICAL RESEARCH, 2006, 111: C03010.

[6] 鞠曉蕾，沈云中，張子占．基于GRACE衛星RL05數據的南極冰蓋質量變化分析[J]．地球物理學報，2013，56(9) :2918-2927．

[7] LLOVER W, GUINEHUT S, CAZENAVE A. Regional and interannual variability in sea level over 2002-2009 based on satellite altimetry, Argo float data and GRACE ocean mass [J]. Ocean Dynamic, 2010, 60: 1193-1204.

[8] FOFONOFF P, MILLARD R C Jr (1983). Algorithms for computation of fundamental properties of seawater[J]. UNESCO Technical Paper in Marine Science UNESCO,Paris,1983,44(53):17-18.

[9] DHOMPS A L, GUINEHUT S, TRAON P Y Le, et al. A global comparison of Argo and satellite altimetry observations[J].Ocean Science Discussions, 2010, 7: 995-1015.

[10] PAULSON A，ZHONG S J，WAHR J．Inference of mantle viscosity from GRACE and relative sea level data．Geophysical Journal International，2007，171(2):497-508．