趙宗澤，劉榮杰，馬 毅，孫偉富

(1. 山東科技大學，山東 青島 266590；2. 國家海洋局 第一海洋研究所，山東 青島 266061)

海岸線是海洋與陸地的分界線。中華人民共和國國家標準《海洋學術(shù)語：海洋地質(zhì)學》(GB/T 18190-2000)給出的海岸線定義是：海陸分界線，在我國系指多年大潮平均高潮位時海陸分界線[1]。海岸線是國際地質(zhì)科學聯(lián)合會(IUGS)提出的27個全球“地質(zhì)指標”之一[2]，它不僅標識了沿海地區(qū)的水陸分界線，而且蘊含著豐富的環(huán)境信息，其變化直接改變潮間帶灘涂資源量及海岸帶環(huán)境，可引起海岸帶多種資源與生態(tài)過程的改變，影響沿海人民的生存發(fā)展[3]。因此，研究海岸線變化對于了解海岸帶生態(tài)環(huán)境變化乃至全球變化具有重要意義。

湄洲灣位于福建省海岸中段，三面環(huán)陸，是深入內(nèi)陸的半封閉狹長形海灣、自然條件優(yōu)越[4]。1989年，湄洲灣被交通部列為我國四大遠景發(fā)展規(guī)劃的國際深水中轉(zhuǎn)港之一，是福建省主要建設的三大港口之一。1999-12，國務院批準秀嶼港、東吳港區(qū)為對外輪開放的一類口岸，湄洲灣港口群已與世界上27個國家和地區(qū)的近50個港口通航[5]。港口的發(fā)展推動了湄洲灣的劇烈變化。為了合理規(guī)劃與管理湄洲灣岸線，必須掌握湄洲灣岸線現(xiàn)狀及變化。

Tuncay等利用Landsat TM和ETM+影像的水體指數(shù)進行閾值處理，提取出土耳其5個海岸帶濕地的水邊線，并利用DSAS (Digital Shoreline Analysis System)對岸線進行變遷分析[6]。Sheik等利用6期(1999，2001，2003，2005，2007，2009年)IRS和Landsat數(shù)據(jù)對印度南部科摩林角和杜蒂戈林之間的岸線進行提取，將研究區(qū)分為4塊，每塊再分3個子區(qū)，利用EPR(End Point Rate)、LRR(Linear Regression Rate)、LMS(Least Median Square)以及JKR計算出岸線的變化率，然后根據(jù)岸線的變化率分析岸線的侵蝕與增長[7]。遙感具有范圍廣，多時相等特點，目前國內(nèi)外普遍采用遙感手段對海岸線進行變遷監(jiān)測與分析[8-14]，但未見關(guān)于湄洲灣海岸線遙感變遷的研究。

因此，本文首先利用1983年、1993年、2001年和2010年的衛(wèi)星影像對湄洲灣海岸線進行提取，然后利用基線法和面積法分別對1983-1993年、1993-2001年和2001-2010年海岸線進行變遷分析，最后對湄洲灣四期海岸線進行曲率分析，得出四期海岸線整體曲折度的變化。

1 研究區(qū)域

湄洲灣位于福建省中部沿海，北鄰興化灣，南鄰泉州灣，東面為莆田市，西面為惠安縣，西北面為仙游縣，共屬于泉州市與莆田市管轄，是福建沿海天然優(yōu)良港灣之一，亦是“中國少有，世界不多”的多泊位天然深水良港，以灣口的湄洲島命名。灣內(nèi)主要有秀嶼港區(qū)、東吳港區(qū)、肖厝港區(qū)和斗尾港區(qū)四大港區(qū)。灣內(nèi)島嶼層層阻擋，口內(nèi)有盤嶼、大竹島、大生島，灣內(nèi)又有橫嶼和羅嶼形成兩道天然屏障。海岸線曲折，自然岸線主要由基巖岸線組成，局部出現(xiàn)淤泥質(zhì)、砂質(zhì)和生物岸線，研究區(qū)如圖1所示。

2 數(shù)據(jù)和方法

2.1 數(shù)據(jù)與處理

為了監(jiān)測湄洲灣海岸線近30 a來的變遷情況，選取1983年的Landsat MSS，1993年的Landsat TM，2001年的Landsat ETM+和2010年的Landsat TM四期影像作為數(shù)據(jù)源，數(shù)據(jù)信息見表1。分別對MSS影像綠、紅、近紅外，TM和ETM+影像的藍、綠、紅、近紅外、短波紅外、中紅外進行波段組合，然后以正射校正精度小于兩個像元(20 m以內(nèi))的SPOT-5影像作為參考，對1983年、1993年、2001年和2010年四期影像進行幾何糾正。

圖1 湄洲灣地理位置Fig.1 Geographical position of the Meizhou Bay

序號軌道號成像時間傳感器空間分辨率/m1119/0431983-09-27MSS802119/0431993-06-26TM303119/0432001-03-04ETM+304119/0432010-05-24TM30

注：序號表示4幅影像

2.2 方 法

采用目視解譯的方法對湄洲灣海岸線進行提取。借鑒現(xiàn)場踏勘資料，根據(jù)海岸線類型在影像中的色彩、紋理、地物鄰接關(guān)系特征，建立海岸線類型的解譯標志，從而提取湄洲灣1983年、1993年、2001年和2010年四期海岸線。

在信息提取結(jié)果的基礎(chǔ)上，利用基線法和面積法對1983-1993年、1993-2001年和2001-2010年湄洲灣海岸線進行變遷監(jiān)測與分析。并利用湄洲灣四期海岸線的曲率分布來反映海岸線的曲折程度。

3 結(jié)果與分析

3.1 1983-2010年湄洲灣海岸線變化

由表2可知，近30 a來湄洲灣海岸線長度整體處于增加態(tài)勢，海岸線長度增加了2.99 km，主要原因是養(yǎng)殖池塘的圍填，其次是港口建設。1993年比1983年岸線長度新增1.28 km，2001年比1993年岸線長度新增2.38 km，而2010年海岸線長度比2001年岸線長度減少了0.67 km，其中2000—2010年海岸線長度減少的主要原因是人工建堤。

近30 a來，由于湄洲灣四期海岸線整體處于一直向海增長狀態(tài)，導致湄洲灣海灣面積減少了70 km2，其中1993年較1983年減少了13 km2，2001年較1993年減少了5 km2，2010較2001年減少了52 km2。

表2 1983-2010年湄洲灣四期海岸線長度和海灣面積Tabel 2 The coastline length and the bay areas in four periods (from 1983 to 2010)in the Meizhou Bay

3.2 湄洲灣海岸線變遷分析

DSAS是基于基線法原理的岸線分析系統(tǒng)，本文利用DSAS系統(tǒng)實現(xiàn)基線法對湄洲灣海岸線的變遷分析。在基線上總共生成了1 303個點，從這些點向四期岸線做橫斷線，從而可以計算出各點處兩期岸線之間的增長(或侵蝕)距離。

利用基線法和面積法分別對1983-1993年、1993-2001年和2001-2010年湄洲灣海岸線變遷情況進行分析。

3.2.1 1983-1993年海岸線變遷分析

1983年至1993年湄洲灣海岸線向海最大增長距離為3 456 m，位于泉州市泉港區(qū)，主要是由人工建堤造成；最大侵蝕距離為746 m，位于惠安縣東北部；平均增長距離為183 m；平均侵蝕距離為44 m。如圖2所示，1983年至1993年湄洲灣岸線變遷導致海岸帶陸地面積增加15.97 km2，海岸帶陸地面積減少2.85 km2，陸地面積凈增加13.13 km2，平均每年凈增加1.13 km2。

1983-1993年岸線變遷主要由養(yǎng)殖池塘圍填、人工建堤和肖厝港區(qū)建設引起。莆田市仙游縣楓亭鎮(zhèn)和莆田市秀嶼區(qū)忠門鎮(zhèn)出現(xiàn)了一些養(yǎng)殖池塘圍填，造成了海岸線的變遷；1985-1988年中交水運規(guī)劃設計院在泉州市泉港區(qū)設計和建成了肖厝萬噸級雜貨碼頭，使海岸線向海增長；人工建堤位于泉州市泉港區(qū)南埔鎮(zhèn)，使岸線產(chǎn)生了較大的變遷。

3.2.2 1993-2001年海岸線變遷分析

1993—2001年湄洲灣海岸線向海最大增長距離為2 123 m，位于莆田市秀嶼區(qū)南部，主要由圍填海養(yǎng)殖造成；最大侵蝕距離為118 m，位于莆田市與泉州市交界處；平均增長距離為76 m，平均侵蝕距離為26 m。如圖3可知，1993年至2001年湄洲灣岸線變遷導致海岸帶陸地面積新增6.04 km2；海岸帶陸地面積減少0.52 km2；陸地面積凈增加5.52 km2，平均每年凈增加0.69 km2。

1993-2001年岸線變遷主要由養(yǎng)殖池塘圍填引起，零散分布于莆田市的城廂區(qū)靈川鎮(zhèn)和秀嶼區(qū)東埔鎮(zhèn)。

圖2 1983-1993年由海岸線變遷引起的湄洲灣面積變化Fig.2 Changes of the Meizhou Bay area caused by the coastline changes from 1983 to 1993

圖3 1993-2001年由海岸線變遷引起的湄洲灣面積變化Fig.3 Changes of the Meizhou Bay area caused by the coastline changes from 1993 to 2001

3.2.3 2001-2010年海岸線變遷分析

2001—2010年湄洲灣海岸線向海最大增長距離為3 406 m，位于惠安縣東部，主要由人工建堤造成。最大侵蝕距離為543 m，位于莆田市城廂區(qū)。平均增長距離為877 m，平均侵蝕距離為46 m。1993—2001年湄洲灣海岸線變遷導致海岸帶陸地面積新增53.98 km2，海岸帶陸地面積減少2.45 km2，陸地面積凈增加51.74 km2，平均每年凈增加5.75 km2。

圖4 2001-2010年由海岸線變遷引起的湄洲灣面積變化Fig.4 Changes of the Meizhou Bay area caused by the coastline changes from 2001 to 2010

2001-2010年湄洲灣海岸線變遷主要由外走馬埭圍墾、斗尾港區(qū)建設和養(yǎng)殖池塘的圍填引起。2004-08-15，惠安縣外走馬埭圍墾工程正式開工建設，2009-05，14.011 km的海堤工程整體建設任務全面完成，造成了海岸線的大面積變遷；2006-04-29，位于泉州市惠安縣的斗尾港區(qū)青蘭山原油碼頭開始投入建設，2009-02-16，青蘭山30萬噸級原油碼頭正式完工，引起了海岸線的向海增長；養(yǎng)殖池塘圍填主要分布于泉州市泉港區(qū)南埔鎮(zhèn)和莆田市的仙游縣楓亭鎮(zhèn)和秀嶼區(qū)山亭鄉(xiāng)。

3.3 湄洲灣海岸線曲折度分析

曲率是幾何體不平坦程度的一種衡量，曲率越大，表示曲線的彎曲程度越大。所以,可以利用曲率分布來反映湄洲灣四期海岸線的曲折程度。

計算出湄洲灣海岸線每年在不同離散點數(shù)情況下的曲率分布，并利用離散點的曲率分布來反映湄洲灣海岸線整體的曲折程度，為了去除細節(jié)噪聲的影響，選取閾值M，并設ai和ai+1為兩相鄰離散點，兩點坐標分別為(axi,ayi)和(axi+1,ayi+1)，其中設d為max(|axi+1-axi|, |ayi+1-ayi|)，s為兩點的曲線距離。當s/d>M時，進行曲率計算。當s/d≤M時，連接下個鄰接點，重復以上過程。通過以上迭代方法，生成新的離散點，可以有效地反映湄洲灣海岸線的曲率分布。

根據(jù)分析，為了較好并準確地反映湄洲灣海岸線的曲率分布，根據(jù)經(jīng)驗并開展大量的對比試驗，分析得出閾值為1.6，離散點數(shù)為4 000時的海岸線曲率分布，可以較準確的反映湄洲灣四期海岸線的曲折程度，如圖5所示，利用matlab編程繪制出了1983年、1993年、2001年和2010年湄洲灣的曲率分布。

圖5 湄洲灣曲率分布圖Fig.5 Distribution of the curvature of the Meizhou Bay

從表3中可以看出，中值可以較好的反映1983-2010年湄洲灣海岸線的整體曲折程度。1983-2010年曲率分布的中值在增加，說明整體曲折度一直在增大，與1983-2010年湄洲灣海岸線變遷情況相對應，其中2001-2010年海岸線整體曲折度變化最大，1993-2001年變化最小。這說明養(yǎng)殖池塘圍填和港口建設不僅造成了湄洲灣海岸線向海增長，同時也造成了湄洲灣海岸線整體曲折程度的增加。

由曲率分布圖可知，1983年，湄洲灣海岸線曲折度較大區(qū)域主要位于莆田市秀嶼區(qū)。1993年，由于肖厝港區(qū)建設和人工建堤，泉州市泉港區(qū)海岸線曲折度增大；莆田市仙游縣出現(xiàn)圍填海養(yǎng)殖，仙游縣海岸線曲折度也有所增大。2001年，莆田市秀嶼區(qū)東埔鎮(zhèn)由于圍填海養(yǎng)殖，造成其岸線曲折度增大；2010年，由于秀嶼港區(qū)和斗尾港區(qū)建設，其附近海岸線曲折度增大；泉州市泉港區(qū)出現(xiàn)了大量圍填海養(yǎng)殖，使得泉州市泉港區(qū)整體曲折度增大。

表3 1983-2010年湄洲灣四期海岸線的曲率分布統(tǒng)計數(shù)據(jù)Tabel 3 Statistis of coastline curvature distributions in the Meizhou Bay from 1983 to 2010

4 結(jié) 論

1983-2010年湄洲灣海岸線長度整體處于增長態(tài)勢，海灣面積減少，海岸線一直向海增長，海岸線曲折度增大，主要由養(yǎng)殖池塘圍填、人工建堤和港口建設等人類活動引起。

泉州市惠安縣：1983-2001年，海岸線變遷緩和；2001-2010年，由于外走馬埭圍墾和斗尾港區(qū)青蘭山原油碼頭建設，變遷面積較大，同時也造成了岸線曲折度的增大。

泉州市泉港區(qū)：1983-1993年，由于人工建堤和肖厝港區(qū)建設，使得海岸線有大的變遷，同時造成了岸線曲折度的增大；1993-2001年，海岸線變遷較?。?001-2010年，出現(xiàn)了養(yǎng)殖池塘圍填，造成了海岸線向海增長和曲折度的增大。

莆田市仙游縣：1983-1993年，楓亭鎮(zhèn)出現(xiàn)養(yǎng)殖塘圍填，使海岸線向海增長，海岸線曲折度增大；1993-2001年，海岸線變遷緩和；2001-2010年，楓亭鎮(zhèn)又出現(xiàn)較大面積養(yǎng)殖池塘圍填，使海岸線向海增長。

莆田市城廂區(qū)：1983-2010年，海岸線無較大變遷。

莆田市秀嶼區(qū)：1983-1993年，忠門鎮(zhèn)出現(xiàn)少許圍填海養(yǎng)殖，造成了海岸線向海增長和曲折度增大；1993-2001年，東埔鎮(zhèn)出現(xiàn)了圍填海，造成了海岸線向海增長；2001-2010年，山亭鄉(xiāng)出現(xiàn)大面積的圍填海，使海岸線出現(xiàn)了較大的變遷。

由于所用的MSS影像分辨率較低，對目視解譯造成了一定的困難，使得提取結(jié)果存在偏差，與TM和ETM+影像的分辨率有一定的差別，使得變遷分析結(jié)果存在一定誤差。但總體上，監(jiān)測結(jié)果能真實地反映出岸線變化的規(guī)模和趨勢。

[1] GB/T18190—2000 海洋學術(shù)語:海洋地質(zhì)學[S].北京:中國標準出版社,2000.

[2] LI R, DI K, MA R. 3-D shoreline extraction from IKONOS satellite imagery [J]. Marine Geodesy, 2003,(26):107-115.

[3] 孫才志,李明昱. 遼寧省海岸線時空變化及驅(qū)動因素分析[J]. 地理與地理信息科學,2010,26(3):63-67.

[4] 高東旭, 劉允召. 福建省湄洲灣港口建設[J]. 水運工程,1996,(10):59-63.

[5] 陳明蔚, 陳貴華. 湄洲灣港港口物流在新格局下的發(fā)展對策[J]. 交通企業(yè)管理,2010,8(264):56-57.

[6] TUNCAY K, ABDULAZIZ G, FEVZI K, et al. Automatic detection of shoreline change on coastal Ramsar wetlands of Turkey [J]. Ocean Engineering,2011,(38):1141-1149.

[7] SHEIK M, CHANDRASEKAR. A shoreline change analysis along the coast between kanyakumari and tuticorin, india, using digital shoreline analysis system[J]. Geo-spatial Information Science, 2011, 14(4):282-293.

[8] YU K, HU C, MULLER-KARGER F E, et al. Shoreline changes in west-central Florida between 1987 and 2008 from Landsat observations[J]. International Journal of Remote Sensing, 2011,32(23):8299-8313.

[9] CHEN W W, CHANG H K. Estimation of shoreline position and change from satellite images considering tidal variation [J]. Estuarine, Coastal and Shelf Science,2009,84(1):54-60.

[10] 孫偉富, 馬毅, 張杰, 等. 不同類型海岸線遙感解譯標志建立和提取方法研究[J].測繪通報, 2011,(3):41-44.

[11] VINAYARAJ P, JOHNSON G, DORA G U, et al. Quantitative estimation of coastal changes along selected locations of Karnataka, India: A GIS and remote sensing approach[J]. International Journal of Geosciences,2011,(2):385-393.

[12] ZHANG Y. Coastal environmental monitoring using remotely sensed data and GIS techniques in the Modern Yellow River delta, China[J]. Environ Monitoring Assess,2011,(179):15-29.

[13] 孫云華, 張安定, 王慶. 基于RS和GIS的近30年來人類活動影響下萊州灣東南岸海岸濕地演變[J].海洋通報, 2011,1(30):66-72

[14] 馮永玖, 劉丹, 韓震. 遙感和GIS支持下的九段沙岸線提取及變遷研究[J].國土資源遙感,2012,1(92):65-69.