金波文，王慧， 李經(jīng)緯， 李文善， 李歡， 徐浩

（1.國家海洋信息中心，天津 300171；2.中國地震局第一監(jiān)測中心，天津 300180）

在全球氣候變暖背景下，海平面上升已成為全球性重大環(huán)境問題。由于溫度升高和南北極冰川、冰蓋融化等環(huán)境因素的影響，在過去100a 內(nèi)全球海平面上升了10～25 cm[1]。20 世紀以來，全球海平面呈現(xiàn)加速上升趨勢，我國沿海海平面上升速率高于全球同期平均水平，面臨的海平面上升影響風險更大[2-3]。目前，中國沿海海平面監(jiān)測的主要方法是通過驗潮站潮位觀測，其觀測得到的是基于各站水尺零點的相對海平面變化，除了受絕對海平面變化影響外，還與局部地區(qū)陸地垂直運動有關(guān)，后者在很多地區(qū)的重要性甚至超過了前者，并且陸地沉降會加劇相對海平面上升[4-5]。此外，由于驗潮站相對海平面變化包含了其本身的陸地垂直運動因素在內(nèi)，在對近岸融合衛(wèi)星高度計絕對海平面資料進行綜合研究分析時，陸地垂直運動是一個不可忽視的因素。因此，準確地監(jiān)測中國沿海驗潮站及其鄰近地區(qū)的陸地垂直運動規(guī)律不僅對于驗潮水尺零點的維護具有重要意義，而且對聯(lián)合衛(wèi)星高度計數(shù)據(jù)的沿海絕對海平面變化研究分析至關(guān)重要。

沿海地區(qū)陸地垂直運動主要包含大尺度范圍的地殼垂直運動和局部地面垂直運動。大尺度范圍內(nèi)的地殼垂直運動主要包括地殼板塊運動和冰川均衡調(diào)整（Glacial Isostatic Adjustment，GIA）等[6-7]，如北歐斯堪的納維亞半島附近由于其陸上冰川融化陸地抬升的GIA 效應引起了海平面下降[8]。局地地面垂直運動主要由大河三角洲沉積壓實效應、人類活動開采地下水和超高超密建筑物修建等因素引起[9-10]。監(jiān)測陸地垂直運動的傳統(tǒng)方法是定期精密水準測量，黃立人等利用精密水準復測資料對中國東部沿海地區(qū)的近代地殼垂直運動分析得出：東部沿海廣大平原地區(qū)以下降為主，遼西山地、膠東半島以及閩、粵、貴山地和丘陵則以緩慢的塊體上升為主[11-13]。隨著現(xiàn)代大地測量GNSS 技術(shù)的發(fā)展，基于連續(xù)GNSS 觀測提供精確地面沉降結(jié)果的方法得到了快速的發(fā)展和應用。為剝離驗潮站地面垂直運動因素影響，Bevies 等利用GNSS 連續(xù)觀測數(shù)據(jù)對西北太平洋、黑海等沿岸的驗潮站陸地垂直運動信息進行了研究[14-17]，國內(nèi)焦文海等利用驗潮站附近的GNSS 研究了驗潮站地面垂直運動后的絕對海平面變化，但在中國沿海地區(qū)直接利用GNSS 觀測數(shù)據(jù)，整體分析全國驗潮站及其附近區(qū)域地面垂直運動的研究仍較少[18-20]。1993 年以來，隨著衛(wèi)星高度計的發(fā)射，許多學者利用驗潮站和附近的衛(wèi)星高度計觀測數(shù)據(jù)就如何推算驗潮站附近的陸地升降開展了研究，該方法最早由Cazenave 等提出，Kuo 等分別在該方法的基礎上進行了改進，并對地中海和大西洋東部等地區(qū)驗潮站的地面垂直運動進行了研究[21-22]。董鴻聞等利用類似方法分別推算了青島驗潮站和黃河三角洲的地面垂直運動，劉首華等利用該方法推算了渤黃海周邊驗潮站地面垂直運動速率，謝書誼等利用全球191 個和中國沿海29個驗潮站資料聯(lián)合衛(wèi)星高度計資料推算了沿岸陸地垂直運動速率[23-25]。

我國于2009 年開始陸續(xù)在全國沿海50 余個驗潮站增加了GNSS 連續(xù)觀測業(yè)務，GNSS 連續(xù)觀測為提取驗潮站的地面垂直運動信息提供了重要的數(shù)據(jù)來源。本文利用2009—2018 年間的驗潮站GNSS 觀測資料，聯(lián)合驗潮站周邊的中國地殼運動觀測網(wǎng)絡（以下簡稱陸態(tài)網(wǎng)絡） GNSS 基準站資料，對全國沿海驗潮站的地面垂直運動進行了分析。同時，利用沿海衛(wèi)星高度計資料和驗潮站海平面觀測資料推算驗潮站的地面垂直運動，與GNSS獲取的垂直運動進行了對比分析，并將分析結(jié)果與冰川均衡調(diào)整模型的預測結(jié)果進行了對比研究。最后給出中國沿海驗潮站及其鄰近區(qū)域最新的陸地垂直運動情況，可為我國沿海相對海平面變化研究、影響評估和預測提供支撐。

1 數(shù)據(jù)來源和處理方法

1.1 數(shù)據(jù)來源

本文研究數(shù)據(jù)主要包括驗潮站觀測數(shù)據(jù)、GNSS觀測數(shù)據(jù)和衛(wèi)星高度計海平面異常月均值數(shù)據(jù)。

驗潮站實測數(shù)據(jù)來源于中國沿海海洋觀測網(wǎng)，通過實測數(shù)據(jù)得到的月均海平面數(shù)據(jù)集已通過全國海洋站基準潮位核定，對因站址變遷、環(huán)境改變、儀器更換、水尺變動和觀測手段改變等引起的水尺零點變動和資料均一性問題進行了訂正。

衛(wèi)星高度計資料來源于哥白尼環(huán)境監(jiān)測中心（www.psmsl.org）提供的多衛(wèi)星融合訂正后的海平面異常月均值數(shù)據(jù)（MSLA，Mean Sea Level Anomalies），參考框架為1993—2012 年平均海平面，數(shù)據(jù)已進行儀器誤差和相應的地球物理影響因素的訂正。

GNSS 連續(xù)觀測數(shù)據(jù)來源于海洋站GNSS 連續(xù)觀測網(wǎng)以及陸態(tài)網(wǎng)絡GNSS 基準站網(wǎng)。海洋站GNSS 資料來源于海洋站GNSS 實時觀測數(shù)據(jù)。陸態(tài)網(wǎng)絡GNSS 基準站資料來源于中國地震局東部形變數(shù)據(jù)分中心（www.eqdsc.com）提供的時間序列數(shù)據(jù)，采用ITRF2014 坐標框架。陸態(tài)網(wǎng)絡在華北地區(qū)站點數(shù)較多，距離驗潮站較近的站點也比較多，東部和南部沿海地區(qū)站點分布較少。

1.2 驗潮站和衛(wèi)星高度計數(shù)據(jù)處理

為了匹配各GNSS 觀測資料的時間跨度，同時考慮海平面季節(jié)變化帶來的影響，本文對各驗潮站月均值數(shù)據(jù)和衛(wèi)星高度計數(shù)據(jù)的時段進行了截取，各站數(shù)據(jù)時段均為起始年的1 月到結(jié)束年的12 月，總體時段為2009—2018 年。在選取沿海衛(wèi)星高度計資料時，一種是把臨近海域的高度計觀測數(shù)據(jù)進行整體月平均，從而避免高度計在某一格點觀測時間的不連續(xù)性帶來的較大誤差[22]。另一種是選取與驗潮站時間序列相關(guān)系數(shù)最高的格點數(shù)據(jù)作為驗潮站絕對海平面時間序列[26]?？紤]到衛(wèi)星高度計在沿岸地區(qū)的精度不如寬闊海域，因此本文搜索各驗潮站1.5°半徑內(nèi)的高度計格點數(shù)據(jù)，并按反距離加權(quán)獲得平均值作為該站衛(wèi)星高度計絕對海平面資料序列，并對各驗潮站海平面序列與提取的對應站點衛(wèi)星高度計海平面序列進行了相關(guān)性分析，各站相關(guān)系數(shù)如表1 所示[24]。除日照、連云港、北海和防城港4 站驗潮站序列與衛(wèi)星高度計序列相關(guān)系數(shù)小于0.6 外，其余站點相關(guān)性均較高（相關(guān)系數(shù)大于0.8的站點占68%），能較好地代表驗潮站區(qū)域的海平面變化，圖1 為三亞站驗潮站與衛(wèi)星高度計海平面觀測序列。

表1 各站驗潮站資料與衛(wèi)星高度計資料相關(guān)系數(shù)

圖1 三亞站驗潮站與衛(wèi)星高度計海平面觀測序列

計算驗潮站和衛(wèi)星高度計海平面變化速率時，首先去除月均值序列的季節(jié)性信號，進而獲取各站點的年均值海平面序列，最后利用最小二乘法對時間序列進行線性回歸求得海平面變化速率[25]，本文在海平面變化趨勢提取中使用的模型為：

式中：MSLt為驗潮站時間為t 的年平均海平面觀測值；a0為該時段內(nèi)的平均海平面值；b0為平均海平面線性變化速率；nt為擬合誤差。

1.3 GNSS 數(shù)據(jù)處理

為了比較中國沿海和大陸的相對運動，并將海洋站GNSS 觀測結(jié)果統(tǒng)一至全球參考框架下，本文將所有海洋站GNSS 觀測數(shù)據(jù)與中國周邊IGS（International GNSS Service）跟蹤站觀測數(shù)據(jù)聯(lián)合解算。采用GAMIT/GLOBK（10.7）軟件進行基線處理和平差解算?；€處理時軌道選用IGS 發(fā)布的精密軌道星歷，衛(wèi)星高度角設為10°，大氣延遲改正采用SAASTAMOINEN 模型，映射函數(shù)選用維也納映射模型，并進行了海洋潮汐、固體潮汐和極移等改正。平差解算時利用GLOBK 對所有時段單天解進行整體平差，得到各站基于ITRF2014 框架的坐標時間序列，時間序列質(zhì)量控制采用相鄰兩個觀測值之差和可視化人工判讀方法對各站時間序列粗差逐個進行剔除，得到最后的時間序列[27]。

1.4 陸地垂直運動速率計算

1.4.1 直接法（GNSS 垂直運動速率） 國內(nèi)外眾多學者對GNSS 高程分量進行過研究，普遍認為全球許多連續(xù)站單天解算坐標分量的時間序列周期主要有年周期和半年周期兩項，因此本文使用與Li 等相同的高程時間序列模型對單站進行建模，并利用最小二乘法求得各站垂直運動速率[28-31]。

1.4.2 間接法（絕對海平面變化速率減去相對海平面變化速率：A-TG） 通過式（1）模型分別計算各驗潮站相對海平面變化速率和衛(wèi)星高度計絕對海平面變化速率，再利用各站對應的絕對海平面變化速率減去相對海平面變化速率，得到的差值即為陸地垂直運動速率。陸地垂直運動速率與驗潮站和衛(wèi)星高度計觀測的海平面變化速率之間的關(guān)系可用下式表示：

式中，Va為衛(wèi)星高度計觀測的絕對海平面變化速率，VTG為驗潮站觀測的相對海平面變化速率，Vh為陸地垂直運動速率。

2 中國沿海陸地垂直運動分析

中國大陸構(gòu)造環(huán)境監(jiān)測網(wǎng)絡（以下簡稱“陸態(tài)網(wǎng)”）在全國范圍內(nèi)建設了260 個GNSS 連續(xù)觀測基準站，本文選取中國沿海地區(qū)43 個陸態(tài)網(wǎng)基準站資料并利用直接法分析了地面運動速率，結(jié)果如圖2 所示，并與選取衛(wèi)星高度計和驗潮站資料采用間接法得到的陸地垂直運動趨勢進行了對比分析。

2.1 沿海地區(qū)陸態(tài)網(wǎng)基準站陸地垂直運動

中國沿海地區(qū)陸態(tài)網(wǎng)基準站觀測結(jié)果顯示，中國沿海省區(qū)市陸地總體以上升為主，局部區(qū)域陸地呈下降趨勢。其中，遼寧省5 個基準站均呈上升趨勢，平均上升速率為1.46± 0.11 mm/a，沈陽站上升速率最大（3.04±0.07 mm/a），丹東站上升速率最?。?.25±0.04 mm/a）。河北省滄縣站位于華北平原沉降區(qū)，沉降速率為-29.71±0.07 mm/a，承德和唐山站均呈上升趨勢。天津市北部山區(qū)的薊縣站和寶坻站呈上升趨勢，上升速率較小，分別為0.94±0.01 mm/a 和0.41±0.04 mm/a，位于華北平原地區(qū)的武清和濱海站沉降較快，沉降速率分 別為-44.72 ± 0.07 mm/a 和-16.91 ± 0.04 mm/a。山東省8 個基準站均呈上升趨勢，平均上升速率為0.91 ± 0.11 mm/a，臨沂站上升速率最大（1.75 ±0.04mm/a），煙臺站上升速率最?。?.43±0.04mm/a）。江蘇省連云港站較為穩(wěn)定，沉降速率為-0.03 ±0.04 mm/a，鹽城和南通站沉降速率分別為-1.59 ±0.06 mm/a 和-0.95±0.04 mm/a，溧水站呈緩慢上升趨勢，上升速率為0.41±0.04 mm/a。上?；鶞收境两邓俾蕿?2.35±0.02 mm/a。浙江省3 個基準站均呈沉降趨勢，平均沉降速率為-0.69±0.07 mm/a。福建省除廈門站以1.59±0.02 mm/a 速率上升外，其余3個基準站均呈沉降趨勢，平潭和霞浦站沉降速率分別為-0.83 ± 0.04 mm/a 和-0.55 ± 0.04 mm/a。廣東省除湛江站以-0.70 ± 0.04 mm/a 速率沉降外，其余3 個基準站均呈上升趨勢。廣西除梧州站以-1.15±0.05 mm/a 速率沉降外，其余2 個基準站均呈上升趨勢。海南?？诤腿齺喺境两邓俾史謩e為-3.82 ± 0.07 mm/a 和-0.01 ± 0.06 mm/a，瓊中和永興島站均呈上升趨勢，上升速率分別為0.67 ±0.02 mm/a 和3.34±0.07 mm/a。

2.2 沿海驗潮站陸地垂直運動

大部分驗潮站GNSS 與附近的陸態(tài)網(wǎng)絡GNSS基準站的垂直運動方向和速率基本一致（圖2），與顧國華等利用陸態(tài)網(wǎng)GNSS 基準站獲取的垂直運動規(guī)律也基本一致，個別驗潮站GNSS 由于位于碼頭，與附近位于基巖的陸態(tài)網(wǎng)基準站得出的陸地垂直運動規(guī)律具有一定的差異[32-34]。中國沿海驗潮站及其鄰近地區(qū)陸地垂直運動總體分布特征如下：

圖2 中國沿海驗潮站及其鄰近區(qū)域ITRF 框架下陸地垂直速度場

（1）遼寧沿海驗潮站及其鄰近區(qū)域陸地垂直運動基本特點是西部上升明顯、東部緩慢上升而局部呈下降趨勢。沿海驗潮站，除東港驗潮站沉降速率為-1.98±0.08 mm/a、鲅魚圈驗潮站沉降速率為-0.19±0.09mm/a 外，其余地區(qū)均呈現(xiàn)不同程度的上升趨勢。遼寧沿海驗潮站平均上升速率約為0.15± 0.16mm/a。陸態(tài)網(wǎng)葫蘆島基準站（1.53 ±0.04 mm/a） 與葫蘆島驗潮站（1.60 ± 0.06 mm/a）距離較近（約11 km），陸地均呈上升趨勢，速率相差0.07 mm/a。

（2）河北沿海的芷錨灣和秦皇島驗潮站均呈上升趨勢，上升速率分別為1.87±0.07 mm/a 和0.73±0.22 mm/a，這與華北平原北側(cè)天津?qū)氎?、天津薊州、河北唐山、河北承德等幾個山地丘陵地區(qū)的陸態(tài)網(wǎng)一致，均呈現(xiàn)上升趨勢。天津沿海的塘沽驗潮站以-8.26±0.10 mm/a 的速率沉降，其位于華北平原濱海地區(qū)，與陸態(tài)網(wǎng)濱海站相距約14 km，均具有較大的沉降速率。華北平原是海河、黃河、淮河等水系在第三紀以來長期下降背景下共同堆積而成的，位于該平原沉降漏斗區(qū)的天津武清、濱海和河北滄縣附近站點（圖2 中黑色方框內(nèi)的站點）區(qū)域沉降最為顯著，主要是由過量開采地下水并疊加構(gòu)造活動引起。

（3）山東和江蘇沿海驗潮站除威海石島站以-3.06±0.05 mm/a 的速率沉降外，其余站點陸地均表現(xiàn)為上升趨勢。從驗潮站和陸態(tài)網(wǎng)基準站可以看出，山東半島整體處于上升趨勢。

（4）上海和浙江沿海驗潮站除朱家尖和小衢山驗潮站緩慢上升外，其余驗潮站均呈下降趨勢。其中，上海蘆潮港驗潮站沉降速率為-3.10±0.05 mm/a，與陸態(tài)網(wǎng)上海站沉降速率相差0.75 mm/a；浙江沿海驗潮站總體以下降為主，與陸態(tài)網(wǎng)基準站趨勢一致。

（5）福建寧德至泉州沿海驗潮站以沉降為主，福建廈門至漳州、廣東和廣西沿海驗潮站均呈現(xiàn)上升與沉降交替出現(xiàn)現(xiàn)象。其中，閘坡驗潮站沉降速率最大，為-3.19±0.09 mm/a，云澳驗潮站上升速率最大，為1.47±0.20 mm/a。

（6）海南島沿海北部秀英驗潮站和南部三亞驗潮站呈下降趨勢，東部清瀾驗潮站、西部東方驗潮站則呈現(xiàn)緩慢上升趨勢。西沙驗潮站上升速率為2.18±0.81 mm/a，與距離約1 km 的陸態(tài)網(wǎng)絡永興島基準站（3.34±0.07 mm/a）上升趨勢一致。

此外，利用間接法得到驗潮站附近的陸地垂直運動情況并與GNSS 速率進行對比，結(jié)果如表2 所示。兩種方法得出的陸地垂直運動速率大小和趨勢方向存在一定差異，其中有31 站的兩種結(jié)果運動方向一致，其余測站運動方向相反。有51%的站點采用直接法和間接法獲得的陸地垂直運動速率差值在-2～2 mm/a 之間。兩種方法結(jié)果產(chǎn)生差異性較大的因素較為復雜，驗潮站周邊高度計數(shù)據(jù)選取方式、高度計在沿岸精度本身受限以及海平面和陸地垂直運動季節(jié)性差異等因素均會對最終的結(jié)果產(chǎn)生影響。

表2 沿海地區(qū)GNSS 直接法和A-TG 間接法計算陸地垂直運動對比情況 單位：mm/a

本文使用的驗潮站和GNSS 站點有33 站是并址運行的，因此GNSS 垂直運動即可認為是驗潮站陸地垂直運動，分析得出有42%的并址站點采用直接法和間接法獲得的陸地垂直運動速率差值在-2～2 mm/a 之間。

4 結(jié)論與討論

沿海陸地垂直運動對于沿海海平面變化分析、海平面變化影響評估以及未來海平面上升預測等具有重要意義，利用GNSS 連續(xù)觀測獲取陸地垂直形變量的技術(shù)已較為成熟，隨著多衛(wèi)星導航系統(tǒng)的聯(lián)合應用，GNSS 連續(xù)觀測精度也將進一步提高。通過與驗潮站并址的GNSS 連續(xù)觀測獲取陸地垂直形變可直接為沿海海平面變化分析提供準確的形變速率。驗潮站GNSS 與陸態(tài)網(wǎng)絡GNSS 基準站獲取的沿海地區(qū)陸地垂直形變速率及方向具有較好的一致性。中國沿海省區(qū)市及沿海驗潮站陸地垂直運動總體表現(xiàn)為遼寧至江蘇沿海上升、上海至福建泉州沿海沉降、福建廈門至廣西沿海升降交替格局，局部濱海平原地區(qū)如華北平原天津南部、河北平原的滄縣則表現(xiàn)為顯著的沉降特征；上海至福建泉州沿海除個別站點上升外，總體呈下降趨勢；福建廈門至廣西沿海、海南島則呈現(xiàn)出上升與沉降交替出現(xiàn)的特征；潿洲島站附近表現(xiàn)為沉降，西沙站區(qū)域呈現(xiàn)出上升趨勢。當然，本文獲得的結(jié)果僅為驗潮站及其附近的陸地垂直運動情況，由于有些驗潮站本身建設在港口碼頭等局部地面不穩(wěn)定的區(qū)域，并不能完全代表東部沿海地殼的垂直運動情況。全面掌握中國東部沿海地殼垂直運動特征仍需進一步利用穩(wěn)定的GNSS 基準站和沿海地區(qū)的水準監(jiān)測成果進行綜合分析。

謝書誼等認為對一個驗潮站附近有多個GNSS站位和多個高度計網(wǎng)格點的情形，分析結(jié)果對GNSS 和高度計數(shù)據(jù)的選擇和處理會非常敏感，甚至會導致陸地沉降的趨勢性改變[25]。本文選用的與驗潮站并址的33 個站點中直接法和間接法得出的結(jié)果差異性也較大，這與GNSS 資料本身的質(zhì)量以及衛(wèi)星高度計在沿岸地區(qū)的精度受限等有關(guān)。隨著多種衛(wèi)星導航系統(tǒng)的聯(lián)合應用，GNSS 連續(xù)觀測獲取陸地垂直形變量的精度和相關(guān)技術(shù)已相對成熟，而衛(wèi)星高度計本身在沿岸地區(qū)的精度尚未得到很好的改善，在GNSS 結(jié)果與A-TG 結(jié)果存在較大差異時，哪種方法更優(yōu)仍有待在遠離大陸的衛(wèi)星高度計精度較高的驗潮站點區(qū)域開展進一步研究。