紅樹林宜林地灘面高程及潮水退干時間時長的測量方法

何碧蘭，馬士龍

(廣州蘢騰園林景觀設計有限公司，廣東 廣州 510520)

1 引言

紅樹林生態系統是生長在潮間帶的特有植物群落，具有防風消浪、促淤護岸、固碳儲碳和維持生物多樣性等重要功能。20世紀，隨著人口增長和社會經濟發展，沿海地區開展的圍墾造田、圍建鹽田、圍海養殖、填海造地等活動，造成我國紅樹林面積大幅萎縮，生態服務功能嚴重衰退。因此，加強紅樹林保護和修復是我國海洋生態文明建設和國土空間生態保護修復的重要內容，自然資源部 國家林業和草原局[1]要求，到2025年，營造紅樹林9050 hm2。其中，廣東5500 hm2、海南2000 hm2、廣西1000 hm2、福建350 hm2、浙江200 hm2。然而，準確掌握造林地海水鹽度、漲退潮的潮水流速、潮間帶所處位置、灘面高程及潮水退干時間、時長等，是紅樹林造林設計的必須內容，尤其灘面高程及潮水退干時間、時長更是預測該立地能否造林的重要制約因素。

眾多研究表明[2～6]，凡高于平均海水面、月平均每天潮水退干時間在12 h以上的潮間帶灘涂均可視為紅樹林的宜林地。然而，目前對紅樹林擬造林地的立地灘面高程及潮水退干時間的測定，大多采用人工直接實測，或參照附近驗潮站的潮汐表進行估測，或按經驗法以現林緣立地灘面為基準點，凡低于20 cm以上的地方可視為紅樹林的宜林地。

然而，人工直接實測至少要持續1個月時間，費工費時，而且約有15 d潮水退干時間在夜間，對實測的技術人員存在著嚴重的安全隱患。參照附近驗潮站的潮汐表進行估測，因常用的RTK測量的高程為CGCS2000大地高程，無法準確測量造林立地灘面的85高程，常常導致估測失誤等。按照經驗法以現有紅樹林林緣立地灘面為基準點，凡低于20 cm以上的地方可視為紅樹林宜林地的方法，估測得到的紅樹林宜林灘面高程是一種經驗值，無法與其他科學數據進行比較，而且在沒有紅樹林的立地灘涂潮間帶無法使用。

針對已有技術的上述缺點，該技術方案通過對湛江紅樹林國際重要濕地保護修復項目(項目編號：GDYD220050)研究探索，并獲得了對擬營造紅樹林的立地灘面潮間帶85高程測定及潮水退干時間、時長估測的簡易方法。實例驗證表明，該技術方案能夠比較準確地測量造林立地灘面的85高程，進而能夠較準確地判斷該灘涂是否屬于紅樹林宜林地的好辦法，對指導生產具有較高的實用價值。

2 研究內容與方法

2.1 擬造林地的平均海水面高程的獲取

(1)因驗潮站的高程系統均為85高程，所以，在網上查找擬造林地附近驗潮站的潮汐表信息，則可以直接獲得該驗潮站近年的平均海水面85高程。

(2)若在網上查找到的擬造林地附近驗潮站的潮汐表信息，沒有平均海水面高程信息，則可以應用該潮汐表22 d時段的每小時潮位高程信息，計算獲得該時段的平均潮位高程，此高程即為該時段的平均海水面85高程[7,8]。

2.2 CGCS 2000大地高程與85高程的異常值的獲取

2.2.1 85高程與CGCS2000大地高程的概念

85高程：1985高程是我國現行的法定基準高程(簡稱85高程)[9]，是采用青島驗潮站1952～1979年的驗潮資料建立的85國家高程基準，并在觀象山上建立了永久水準原點。該水準原點高出黃海平均海水面72.2604 m，并依據“85國家高程基準”推算國家水準網中各水準點的高程，作為全國統一的高程控制系統。85高程基準與1956黃海高程基準的換算：85高程=56黃海高程-0.029(m)。目前所有工程建設控制點的高程基準，幾乎都采用85高程，因此通常稱85高程為正常高程。

CGCS2000大地高程：高程起算自CGCS 2000參考橢球表面[10]。使用RTK接CORS賬號測量，端口配置為8003，RTK測量的高程即為CGCS 2000大地高程。

85高程與CGCS 2000大地高程的區別：85高程是正常高程，CGCS2000大地高是橢圓球面高程[11]。如圖1所示，正常重力線HB為85高程，垂線hB為CGCS 2000大地高程。

85高程與CGCS 2000大地高程的異常值：高程異常值的計算公式為：高程異常值=85高程-CGCS 2000大地高程；或85高程=CGCS 2000大地高程+高程異常值。可設公式：

ζB=HB-hB或HB=hB+ζB

(1)

式(1)中：HB為85高程；hB為CGCS2000大地高程；ζB為高程異常值。

因此，如果有高程異常值，即可把CGCS 2000大地高程轉換為85高程，或將85高程轉換為CGCS 2000大地高程[9]。

必須注意：因橢圓球面位置不同，導致的高程異常值并不是一個固定值[12]，因此，在實際應用時，必須采用RTK CORS測量CGCS2000大地坐標及其高程信息與測區附近的85高程水準點的高程信息進行轉換獲得。

2.2.2 尋找85高程水準點及其坐標和高程信息

通過當地測繪部門找到擬營造紅樹林立地灘面(待測區)附近的任意1～3個(點越多，可靠性越強)85高程水準點及其坐標和高程信息。

2.2.3 測量85高程水準點的CGCS2000大地坐標及其高程

在待測的任意1～3個85高程水準點上，采用網絡RTK的CORS測量的CGCS2000大地坐標和高程信息。

2.2.4 計算2000大地高程與85高程的異常值

把以上測量的CGCS2000大地高程信息與已知的85高程水準點的高程信息進行比較，計算其差值的平均數則為該地的CGCS2000大地與85高程的異常值。

2.3 擬造林地的CGCS2000大地高程測量和宜林地界線確定

在獲取擬造林地的平均海水面85高程數據后，采用網絡RTK的CORS，或大疆經緯M300RTK無人機，對擬造林地的CGCS2000大地高程和相對高程進行測量并找出擬造林地的平均海水面高程位置，同時以該位置為基準點，按等高線測量繪制出該擬造林地的宜林地界線。換言之，該宜林地界線之上潮間帶灘涂均為適宜紅樹林的造林地。

3 結果應用與分析

根據上述方法，通過當地測繪部門可以找到擬造林地附近水準點的85高程信息，并對該點進行測量，可以獲得該點的85高程(正常高)與2000大地高程的異常值。據筍克勤[7]和劉克修等[8]介紹，平均潮位高程就是平均海水面高程，因此，應用附近驗潮站潮汐表的平均潮位高程對照擬造林地潮間帶灘面85高程，可以測量計算擬造林地灘面的平均海水面高程臨界線和潮水退干時間、時長及其日變化情況等，并以此來測量該地的宜林面積。

3.1 實例一

實例一的驗證地點位于硇洲島(北港)驗潮站附近的“湛江紅樹林國際重要濕地保護修復項目”地。

(1)獲取85高程與2000大地高程的異常值。 該驗證地附近水準點的85高程(B85)為1.805 m。應用上述方法，采用RTK接收CORS賬號測量獲得該點的CGCS2000大地高程(B2000)為-13.10 m。根據公式：ζB=HB-hB，計算獲得該點(B點)的85高程(正常高)與2000大地高程的異常值為14.905 m。

(2)平均海水面的潮水退干時間。查看距離最近的硇洲島(北港)驗潮站的潮汐表[13](圖1)計算可知，自2022年5月31日至6月21日的22 d，平均海水面的85高程為180.5 cm(若有對應驗潮站的平均海水面高程則可直接引用)，該地在此高程以上平均每天潮水退干時間為13.0 h。

(3)平均海水面的潮水退干時長及其日變化情況。 查看硇洲島(北港)驗潮站實時潮汐表[13]，獲得該地每天在平均海水面高程以上的潮水退干時間、時長及其日變化情況(圖2)。

(4)平均海水面以上適宜營造紅樹林的面積。 應用高程異常值14.905 m，采用大疆經緯M300RTK無人機測量相對高程并換算為85高程后，繪圖計算獲得“湛江紅樹林國際重要濕地保護修復項目”中的一塊擬造林地在平均海水面85高程180.5 cm(適宜營造紅樹林界線)以上的潮間帶，其灘涂面積為8.5333 hm2。

3.2 實例二

實例二的驗證地點位于湛江港驗潮站附近的潮間帶濕地灘涂擬營造紅樹林的林地。

(1)獲取85高程與2000大地高程的異常值。 該驗證地附近水準點的85高程(B85)為2.066 m。應用上述方法，采用RTK接收CORS賬號測量獲得該點的CGCS 2000大地高程(B2000)為-13.06 m。根據公式：ζB=HB-hB，計算獲得該點(B點)的85高程(正常高)與2000大地高程的異常值為15.126 m。

(2)平均海水面的潮水退干時間。 查看距離最近的湛江港驗潮站的潮汐表[13](圖3)計算可知，自2022年5月15日至2022年6月5日的22 d，平均海水面高程為85高程207.6 cm，該地在此高程以上平均每天潮水退干時間為12.5 h。

(3)平均海水面的潮水退干時長及其日變化情況。 查看湛江港驗潮站實時潮汐表[13]，獲得該地塊每天在平均海水面高程以上的潮水退干時間、時長及其日變化情況(圖3)。

圖3 2022年5月20日湛江港潮汐表示意

(4)平均海水面以上適宜營造紅樹林的面積。 應用高程異常值207.6 cm，采用大疆經緯M300RTK無人機測量相對高程并換算為85高程后，通過繪圖計算獲得上述地塊擬造林地(適宜營造紅樹林界線以上)的潮間帶灘涂面積為5.5113 hm2。

3.2 兩個實例的應用結果分析

(1)兩個實例的平均海水面高程差異較大及其原因。 雖然，硇洲島(北港)驗潮站與湛江港驗潮站直線相距僅60 km，但是，硇洲島(北港)驗潮站的平均海水面高程比湛江港驗潮站低了0.271 m。導致這一現象的原因，主要為硇洲島(北港)驗潮站位于相對比較開闊的洋面邊緣上，而湛江港驗潮站則位于相對比較狹窄的內海灣邊緣上所造成的。這一結果明確提示，今后使用這一技術時，待測的擬造林地距離驗潮站越近，其結果越準確，尤其在洋面狹長、寬度又窄的內海灣與洋面較開闊短淺海灣的差距會更大。這一結果與韓志遠等[14]和章衛勝等[15]的研究結果相似。

(2)經測量轉換的海水面85高程與驗潮站的海水面高程差異不大。 海水面高程在正常氣象狀況下的海水面高程差異較小，在正常平穩的天氣情況下，天文潮潮位與實測潮位相差并不明顯[16]，所以，該實測驗證工作均在正常氣象狀況下進行。

對實例一的實測驗證于2022年5月15日14:00～17:00進行，經4次對林地海水面實測，其結果為林地的海水面85高程與驗潮站的海水面高程差異不大，達到了許永倫等測點高程準確到0.1 m的要求[17]。對實例二的實測驗證于2022年6月10日12:00～15:00進行，經4次對林地海水面的實測，也得到了類似的結果(表1)。

4 結論

4.1 與現有技術相比，該技術方案具有下列優點

(1)與人工實測相比，該技術方案操作更簡便，工作效率更高，實用性更強，安全性更好。

(2)與直接參照附近驗潮站的潮汐表直接估測的結果相比，因該技術方案準確測量了CGCS2000大地高程并換算成85高程，能夠更準確應用附近驗潮站的潮汐表查出平均海水面高程和潮水退干時間、時長及其日變化情況，因此對宜林地的界定更為準確。

(3)與按經驗法以現林緣立地灘面為基準點，凡低于20 cm以上的地方可視為紅樹林宜林地的結果相比，科學性更強，而且可與其他科學數據進行比較，尤其能在沒有紅樹林的潮間帶立地灘涂上使用。

表1 林地實測海水面與驗潮站海水面的高程差異比較

4.2 實踐證明，該技術方案可在生產上參照應用

經2個實例的驗證，在風速3.4～7.9 m/s、浪高0.3～0.5 m的條件下，應用該技術方案測定的85高程與驗潮站的海水面高程誤差在0.1 m以下。但是，兩個實例驗潮站的平均海水面高程差異達0.271 m，究其因主要是由洋面狹長、寬度又窄的內海灣(湛江港)與洋面較開闊短淺的海灣(硇洲島(北港))的區位所造成的。這一結果表明:今后使用這一技術時，待測的擬造林地距離驗潮站越近，其結果越準確。盡管如此，綜合判斷認為該技術方案的科學性較強，可操作性較好，可以在生產上參照應用。