韓小燕，潘曉東，馬林芳，林順利

（溫州海洋環(huán)境監(jiān)測(cè)中心站，浙江溫州 325027）

溫州沿岸平均海平面變化特征分析

韓小燕，潘曉東，馬林芳，林順利

（溫州海洋環(huán)境監(jiān)測(cè)中心站，浙江溫州 325027）

根據(jù)溫州沿岸龍灣、瑞安、鰲江3個(gè)水文站的潮位長(zhǎng)序列（1971～2008年）觀測(cè)資料，進(jìn)行統(tǒng)計(jì)分析，得出溫州沿岸海平面年際及年內(nèi)變化特征，并且得出了溫州淺灘工程靈霓大堤合攏前后，各驗(yàn)潮站海平面的變化特征。分析結(jié)果表明：海平面年際變化3站呈現(xiàn)上升趨勢(shì)，尤以鰲江站上升速率最大，且3站不同時(shí)期均有各自變化規(guī)律；三站9月份的海平面變化對(duì)年變化的貢獻(xiàn)率最大；通過(guò)功率譜分析，3站平均海平面年變化都存在2～3年及4年左右的振蕩周期。

溫州；平均海平面；分布特征；貢獻(xiàn)率；功率譜

1 引言

目前，人類活動(dòng)對(duì)氣候、海洋的影響正不斷擴(kuò)大，溫室效應(yīng)、海平面上升已經(jīng)成為全球性重大環(huán)境問題。海平面上升給人類的生存環(huán)境產(chǎn)生很大的影響，特別是對(duì)沿海地區(qū)經(jīng)濟(jì)發(fā)展、城市安全及人民生活環(huán)境帶來(lái)巨大的威脅：導(dǎo)致風(fēng)暴潮災(zāi)害加劇，增加了城市排澇難度；導(dǎo)致鹽水入侵加劇，破壞水環(huán)境，農(nóng)田鹽漬化加速；致使城市防汛工程功能降低等[1]。許多專家對(duì)全球海平面的變化進(jìn)行了研究，發(fā)現(xiàn)在過(guò)去的100年，全球海平面每年大約以1.5 mm的速度上升[2]；在1993～2003年間全球平均海平面的上升速率為2.89 mm/a[3]；近50年來(lái)中國(guó)平均海平面上升速率為2.5 mm/a[4]；各海區(qū)中東海海平面上升速率平均為3.1 mm/a[5]。以目前狀況發(fā)展下去，到2050年全球海平面將上升5～32 cm[4]，中國(guó)沿海海平面將上升13～50 cm[6]。在全球海平面上升的大背景下，研究溫州沿岸（特別是河口區(qū)）附近海域海平面變化對(duì)溫州沿海地區(qū)經(jīng)濟(jì)發(fā)展、城市安全及人民生活十分重要。龍灣、瑞安、鰲江三站為溫州地區(qū)城市防潮警戒水位核定的代表站，因此，我們搜集了這3個(gè)驗(yàn)潮站（見圖1）的多年潮位觀測(cè)資料，研究溫州沿岸海平面變化的規(guī)律特征。

本文中的海平面為相對(duì)海平面。

2 資料來(lái)源及預(yù)處理

所用的潮位基礎(chǔ)資料來(lái)自溫州水文站，資料序列長(zhǎng)度龍灣、鰲江兩站為38年（1971～2008年），瑞安站位37年（1972～2008年），平均海平面資料選用月、年平均值和極值。據(jù)了解，水文站潮位資料用的基面有測(cè)站基面及凍結(jié)基面，經(jīng)過(guò)換算統(tǒng)一到85國(guó)家高程，見表1。

在收集潮位資料的過(guò)程中，發(fā)現(xiàn)有部分?jǐn)?shù)據(jù)由于儀器設(shè)備故障而缺失，造成潮位數(shù)據(jù)時(shí)間上的不連續(xù)性，為此需要對(duì)缺測(cè)數(shù)據(jù)進(jìn)行補(bǔ)足訂正。對(duì)缺測(cè)記錄不超過(guò)一天的數(shù)據(jù)，采用二次拋物線擬合方法補(bǔ)齊資料序列，具體計(jì)算公式為：

對(duì)缺測(cè)記錄超過(guò)一天的數(shù)據(jù)，可以用自報(bào)值補(bǔ)齊資料序列。具體方法為先用前后15 d或29 d的觀測(cè)值代替缺測(cè)數(shù)據(jù)，然后進(jìn)行潮位調(diào)和分析；然后用第一組調(diào)和常數(shù)自報(bào)出缺測(cè)時(shí)刻的潮位值，用來(lái)代替原先賦予的數(shù)值；進(jìn)行第二次調(diào)和分析，以此類推，直至調(diào)和常數(shù)不變或變化很小為止。此外也可以根據(jù)潮汐性質(zhì)相似原理，參照相近驗(yàn)潮站同時(shí)間實(shí)測(cè)潮位曲線規(guī)律補(bǔ)齊資料序列。

表1 各站潮位序列統(tǒng)計(jì)表

表2 海平面年際變化趨勢(shì)表

圖1 驗(yàn)潮站分布圖

圖2 3站年平均海平面變化曲線

3 平均海平面的分布特點(diǎn)

3.1 年際變化線性模擬

文章對(duì)溫州沿岸3個(gè)驗(yàn)潮站年平均海平面變化進(jìn)行逐年分析（見圖2），可以看出龍灣、瑞安、鰲江三站呈線性遞增趨勢(shì)，龍灣、鰲江兩站上升速率相近，分別為3.6 mm/a、3.7 mm/a，瑞安站上升速率為2.0 mm/a，具體線性變化趨勢(shì)項(xiàng)見表2。

分析結(jié)果表明，溫州沿岸平均海平面上升速率為2.0～3.7 mm/a，平均為3.1 mm/a，與東海區(qū)海平面上升速率持平，略高于我國(guó)沿海海平面及全球海平面上升速率。

3.2 年際變化多項(xiàng)式模擬

溫州地區(qū)海平面年際變化雖然呈上升趨勢(shì)，但在不同時(shí)期各站的海平面變化規(guī)律各有特點(diǎn)，為此對(duì)年平均海平面進(jìn)行多項(xiàng)式趨勢(shì)模擬，見圖3a～3c。

圖3 三站年均海平面多項(xiàng)式模擬

結(jié)果顯示，龍灣站從上世紀(jì)70年代中前期海平面呈現(xiàn)上升趨勢(shì)；70年代中期至80年代中期呈現(xiàn)下降趨勢(shì)；至2000年左右一直維持上升的趨勢(shì)，之后海平面開始下降，下降幅度較大。瑞安站上世紀(jì)70年代初期海平面呈現(xiàn)上升趨勢(shì)；70年代中期至80年代初期呈現(xiàn)下降趨勢(shì)；至90年代初期維持上升的趨勢(shì)，2000年前后海平面開始下降。鰲江站海平面變化與上述兩站位相存在差異，上世紀(jì)70年代初期呈現(xiàn)下降趨勢(shì)；至80年代前期呈現(xiàn)上升趨勢(shì)；80年后前至2000年維持上升的趨勢(shì)；之后海平面開始下降，2006年前后又開始呈現(xiàn)上升趨勢(shì)。

除此，龍灣、瑞安兩站海平面變化的多項(xiàng)式變化模擬還存在11～12年的周期變化。

3.3 月均海平面年際變化

為了解三站各月份海平面年際變化對(duì)年平均海平面的影響及貢獻(xiàn)率，分別對(duì)各站月海平面年際變化進(jìn)行分析，發(fā)現(xiàn)9月份溫州沿岸的海平面普遍最高，貢獻(xiàn)率居全年之首；其次溫州北部5月份海平面次高，而中部及南部海平面次高的月份分別為1月及12月，為冬季月份。具體各月上升速率見表3。

溫州地區(qū)9月份經(jīng)常受臺(tái)風(fēng)影響，南部的鰲江地形為一喇叭形，極易受外海影響，因此鰲江站9月海平面高于其它兩站；5月份進(jìn)入梅雨期，降水量突增，北部的甌江為該地區(qū)主要河流，沿岸支流較多，使得甌江徑流量高于飛云江及鰲江，因此5月龍灣站海平面高于其它兩站；12月、1月沒有徑流及臺(tái)風(fēng)的影響，海平面較高的原因極可能與全球的暖冬現(xiàn)象有關(guān)。

3.4 功率譜

為了更加客觀定量地檢測(cè)時(shí)間序列中的顯著周期，在此采用功率譜的方法進(jìn)行分析，功率譜的計(jì)算方法采用落后自相關(guān)方法。功率譜密度的計(jì)算公式為：

表3各月海平面年際變化上升速率統(tǒng)計(jì)表

通過(guò)上述計(jì)算公式得出的譜密度值有隨機(jī)噪聲的影響，為消除噪聲影響，得到比較光滑和平穩(wěn)的譜密度值，采用Hanning方法進(jìn)行平滑處理，最終得到平滑的功率譜估計(jì)值。

計(jì)算過(guò)程中，三站平均海平面序列的落后自相關(guān)系數(shù)r(1)均為較大正值，表明序列具有連續(xù)性，所以采用紅噪音標(biāo)準(zhǔn)譜對(duì)其進(jìn)行顯著性檢驗(yàn)。圖4給出了顯著性水平α=0.05的平均海平面功率譜圖，圖中橫坐標(biāo)為周期（年），縱坐標(biāo)為譜密度值。

結(jié)果表明，3個(gè)站在周期長(zhǎng)度為2.7 a處均有振蕩，龍灣、瑞安兩站還存在12 a、8 a的振蕩周期，鰲江站振蕩周期較少，且振蕩周期較短，集中在2～4 a（見表4）。

龍灣站在周期長(zhǎng)度為2.7年處的功率譜估計(jì)值為一峰值且超過(guò)標(biāo)準(zhǔn)譜，為顯著周期；在12 a、8 a、3 a及2.4 a處的功率譜估計(jì)值也超過(guò)標(biāo)準(zhǔn)譜，說(shuō)明龍灣站平均海平面在上述周期處存在相應(yīng)的振蕩。

表4各站振蕩周期統(tǒng)計(jì)表

瑞安站功率譜圖與龍灣站相似，在周期長(zhǎng)度為2.7 a處的功率譜估計(jì)值為一峰值且超過(guò)標(biāo)準(zhǔn)譜，為顯著周期；在12 a、8 a及2.4 a處的功率譜估計(jì)值也超過(guò)標(biāo)準(zhǔn)譜。

鰲江站在周期長(zhǎng)度為4.8 a、2.7 a處的功率譜估計(jì)值為一峰值，同樣超過(guò)標(biāo)準(zhǔn)譜，為顯著周期；在2.4 a處功率譜估計(jì)值也超過(guò)標(biāo)準(zhǔn)譜，表明鰲江站在此周期內(nèi)存在振蕩。

三站平均海平面的周期變化與一些常用的理論天文周期相對(duì)應(yīng)，如2 a大氣振動(dòng)周期[3]、24～36個(gè)月準(zhǔn)兩年周期振蕩、2～7 a厄爾尼諾周期、11 a太陽(yáng)黑子周期[3]，5 a平均雨生周期。

3.5 年內(nèi)變化

圖5a是平均海平面的逐月變化曲線，3個(gè)站月變化規(guī)律較一致，表現(xiàn)為9月份海平面最高，6月份次高，6～9月之間存在一個(gè)波谷，發(fā)生在7月。龍灣、瑞安、鰲江2月海平面最低，1～4月海平面變化較小。

考慮到三個(gè)站均布設(shè)在河口區(qū)，徑流的多少直接影響著海平面的高低，因此必須了解每個(gè)區(qū)域的徑流大小。由于可獲取的徑流資料有限，此處僅收集到甌江代表站圩仁站2005年、2006年兩年的徑流資料，繪制兩年月平均徑流變化曲線圖（見圖6），分析徑流對(duì)平均海平面的影響，飛云江、鰲江徑流影響可以以此作為參考。

圖4 三站平均海平面功率譜圖

圖5 各站平均海平面累年月變化曲線

圖6圩仁站月均徑流變化曲線

導(dǎo)致海平面年內(nèi)變化的原因?yàn)?月、9月的徑流明顯高于其余月份，與溫州地區(qū)的梅雨、臺(tái)風(fēng)期降水密切相關(guān)。溫州地區(qū)的梅雨期基本發(fā)生在5～6月（春季），期間產(chǎn)生的大量連續(xù)降水導(dǎo)致河口區(qū)徑流量猛增；9月（夏季）由于臺(tái)風(fēng)影響帶來(lái)豐沛的降水導(dǎo)致徑流相對(duì)增大，但由于每次臺(tái)風(fēng)影響時(shí)間較短，相對(duì)于梅雨期徑流量較少；冬季溫州地區(qū)進(jìn)入枯水期，徑流最小。

平均海平面不僅受到徑流的影響，同時(shí)也受到天文潮位的影響，天文最大潮通常發(fā)生在農(nóng)歷中秋前后，即陽(yáng)歷9月左右，圖5平均海平面累年月變化曲線9月高于6月，可以說(shuō)明天文潮位的影響是該海域夏季海平面變化的主要原因，而梅雨期降水是春季海平面變化的主因。7月該區(qū)域出梅，降水帶北跳，溫州進(jìn)入枯水期，此時(shí)臺(tái)風(fēng)尚未對(duì)溫州產(chǎn)生影響，天文潮位也未達(dá)到天文大潮，導(dǎo)致平均海平面在這段時(shí)期產(chǎn)生下降。

4 沿岸大型海洋工程的影響

近年，由于溫州地區(qū)土地資源緊缺，人多地少的矛盾日益突出，而溫州沿岸灘涂資源眾多，圍涂造地工程此起彼伏，溫州淺灘工程便是其中一項(xiàng)大型利用灘涂資源造地的用?；顒?dòng)。經(jīng)過(guò)工程前期一系列模型的論證及專家的評(píng)審，淺灘工程于2003年下半年開始動(dòng)工。該工程的建設(shè)將新增加13萬(wàn)多畝土地，諸多島嶼連成一體，其中連通靈昆島及霓嶼島的靈霓北堤于2006年4月29日建成通車。溫州沿岸圍涂工程的建設(shè)或者局部或者大范圍地改變了外來(lái)海水進(jìn)入河口區(qū)的通道。

平均海平面年際變化曲線（見圖3a～3c）可以看出，龍灣、瑞安兩站平均海平面在2004年左右開始下降，下降速率龍灣最大。對(duì)三個(gè)站2001～2002年與2004～2005年的潮汐數(shù)據(jù)進(jìn)行調(diào)和分析，結(jié)果顯示，淺水系數(shù)M4/M2龍灣站減小12%，瑞安、鰲江兩站略有增高，短期來(lái)看，與溫州沿岸尤其是甌江流域周邊地區(qū)圍涂造地等各種海洋工程建設(shè)不無(wú)關(guān)系。

5 小結(jié)

（1）溫州沿岸平均海平面年際變化呈上升趨勢(shì)，上升速率鰲江站最大；按月的年際變化來(lái)看，龍灣、瑞安、鰲江三站9月份海平面變化對(duì)年均變化貢獻(xiàn)率最大；

（2）平均海平面年內(nèi)變化，3個(gè)站9月份最高，6月份次高，受天文大潮及徑流影響明顯。此外，還受到各站所處地理環(huán)境、地質(zhì)沉降影響[3]。

（3）從年變化的功率譜分析結(jié)果看，三站均有一個(gè)相同的振蕩周期，為2.7 a。此外，龍灣站存在12 a、8 a、3 a及2.4 a的振蕩周期；瑞安站存在12 a、8 a、2.4 a的振蕩周期；鰲江站存在4.8 a、2.4 a的振蕩周期。龍灣、瑞安兩站的振蕩周期不僅與大氣振動(dòng)周期、準(zhǔn)兩年周期振蕩、厄爾尼諾周期密切相關(guān)，而且還與11 a太陽(yáng)黑子周期活動(dòng)對(duì)應(yīng)。

[1]黃祖珂,黃磊.潮汐原理與計(jì)算[M].青島:中國(guó)海洋大學(xué)出版社，2005.

[2]刑燦飛,廖成旺,李樹德.我國(guó)東部沿海海平面變化研究現(xiàn)狀及問題[J].大地測(cè)量與地球動(dòng)力學(xué),2005,25(3):129-132.

[3]吳濤,康建成,李衛(wèi)江.中國(guó)近海海平面變化研究進(jìn)展[J].海洋地質(zhì)與第四紀(jì)地質(zhì),2007,27(4):123-130.

[4]顏梅,左軍成,傅深波.全球及中國(guó)海海平面變化研究進(jìn)展[J].海洋環(huán)境科學(xué),2008,27（2）:197-200.

[5]辛立國(guó),李廣雪,李西雙.中國(guó)東海2萬(wàn)年來(lái)海平面變化分析[J].中國(guó)海洋大學(xué)學(xué)報(bào),2006,36（5):699-705.

[6]吳中鼎,李占橋,趙明才.中國(guó)近海近50年海平面變化速度及預(yù)測(cè)[J].海洋測(cè)繪，2003,23（2):17-19.

[7]羊天柱,應(yīng)仁方,林顯鈺.浙江沿岸海平面變化的初步研究[J].東海海洋,1993,11(4):1-7.

P731

A

1003-0239（2011）01-0066-06

2010-01-26

韓小燕（1977-），女，工程師，主要從事海洋預(yù)報(bào)工作。E-mail：xiaoyong0728@163.com