章 志，陳 鵬，冒士鳳，常 曼

(江蘇省海涂研究中心，江蘇 南京 210036)

引 言

海岸線對于海岸帶生態(tài)系統(tǒng)興衰、環(huán)境演變及海岸的侵蝕等均起到重要的指示作用[1]?？焖贉蚀_地測定海岸線的動態(tài)變化，對于海域使用管理具有十分重要的意義[2]。江蘇中部地區(qū)，受1855年黃河北歸的影響，以河流輸沙作用而引起海岸線東進，轉(zhuǎn)為以海洋動力作用為主導的侵蝕，海岸線進入調(diào)整時期。1855—1921年，“射陽河口—埒子河口”海岸經(jīng)過了海岸塑造及巨量泥沙重新分配搬運過程，南段海岸一直沿襲著歷史海岸淤進狀態(tài)。20世紀70年代至今，江蘇海岸在海洋水動力作用影響下逐漸變得更加順直流暢，但進退變化特點仍基本沿襲著歷史海岸線的演變蹤跡[3]。2009年江蘇沿海開發(fā)上升為國家戰(zhàn)略以后，受入海河流和海洋的雙重動力作用，以及多來大規(guī)模的灘涂圍墾對海岸地貌及海岸線變化的影響，以往射陽河口為南北沖淤的拐點，現(xiàn)有向南移動的趨勢[4-6]。海岸線監(jiān)測是了解海岸沖淤變化的基礎(chǔ)，對于岸灘沖淤變化具有重要的指示意義。

海岸線測量的傳統(tǒng)方法是設(shè)置測點連接成海岸線，需進行野外實地測量，費時費力。遙感具有快速、大面積、動態(tài)、經(jīng)濟的特點，尤其是高分辨率遙感具有較高的空間分辨率和定位精度，可彌補常規(guī)海岸線測量方法的不足，應用遙感影像探測海岸帶地形信息成為海洋遙感技術(shù)發(fā)展的重點方向[7]。不同學者先后對海岸線遙感提取方法進行了研究，一種是通過目視解譯方式，利用解譯標志，如地物的色調(diào)、紋理、植被生態(tài)位等；另一種是基于潮位校正方法提取嚴格定義岸線[8]。對于海岸變遷研究目前主要有多重緩沖區(qū)法[9]、動態(tài)分割法[10]、剖面計算法[11]、面積法[12]等。本文基于海岸線遙感解譯的結(jié)果，運用海岸線變遷模型，分析了20年來鹽城海岸變遷的過程，探討引起海岸線變化的驅(qū)動力，以期為海岸線保護、整治修復及開發(fā)利用提供一定的決策依據(jù)和參考。

1 研究區(qū)及數(shù)據(jù)源

研究區(qū)位于江蘇省鹽城海岸區(qū)域(圖1)。海岸線北起灌河口團港南側(cè)的“響灌線”陸域分界，南至“安臺線”陸域分界，總長377.88 km，占全省海岸線的42.51%。境內(nèi)灘涂面積683萬畝，占江蘇灘涂面積的70%，并且以每年30～75萬畝左右的成陸速度向大海延伸。依據(jù)岸線變化情況，將研究區(qū)從北到南分為5個部分：“灌河口—廢黃河口”、“廢黃海口—射陽河口”、“射陽河口—新洋港”、“新洋港—斗龍港”和“斗龍港—弶港”。研究利用1997—2017年的landsat7、landsat5、landsat8數(shù)據(jù)，數(shù)據(jù)行列號為120/036、119/037，影像分辨率為30 m。

圖1 研究區(qū)位置圖

2 研究方法

2.1 海岸線分類與提取

海岸線分類目前沒有統(tǒng)一的標準，參考相關(guān)學者的研究成果[13-15]，將鹽城海岸線分為人工和自然岸線。自然岸線以粉砂淤泥質(zhì)為主。遙感解譯特征如下：粉砂淤泥質(zhì)岸線，地質(zhì)類型為粉砂淤泥，灘面坡度平緩，向陸一側(cè)一般植被生長茂盛，呈紅色或暗紅色，向海一側(cè)植被較為稀疏呈淺紅色或沒有植被，裸露潮灘上多有樹枝狀潮溝發(fā)育。影像中，植被茂盛與稀疏程度明顯差異處為淤泥質(zhì)岸線所在位置。河口岸線以最靠近河口的道路橋梁或防潮閘作為河海分界線或以河口突然展寬處的突出點連線作為河海分界線[16]。人工岸線由永久性構(gòu)筑物組成的，向陸一側(cè)不存在平均大潮高潮時海水可達的水域，其向海一側(cè)水陸分界線即是海岸線，一般包括養(yǎng)殖圍堤、鹽田圍堤、碼頭岸線、建設(shè)與交通圍堤。

2.2 岸線變遷分析方法

采用基線法，即以平行于所有歷史岸線基本走向的要素為基準線，構(gòu)建垂直于基線相交的剖面，計算岸線的變化速率，剖面與岸線的交點構(gòu)成岸線的位置序列。基線法主要包括端點變化率(EPR)、線性回歸率(LRR)、最小中位數(shù)平方回歸率(LMS)3 種算法。EPR方法相對于LRR、LMS兩種方法，更能反映岸線的劇烈活動，更能夠反映整個岸線變遷過程中的凈變遷速率[17]。因此選擇EPR算法來定量計算岸線變遷。具體過程如下：

(1)

式中，D1、D2分別為垂線與最早、晚一期岸線的交點到基準線的距離；T1、T2為最早一期岸線時間與最晚一期岸線時間。

具體步驟如下：首先確定基線，繪制大致平行于所提取的海岸線、岸灘走向的分割基線(圖2)。其次，以間隔自北向南沿基線共生成條垂線，采樣間隔為100 m，共產(chǎn)生2486條斷面，刪除其中部分與岸線相交時不合格的垂線，使用ARCGIS軟件運用EPR算法分析岸線變遷過程及變化速率[18]。

圖2 基線及斷面分布

3 結(jié)果與分析

3.1 岸線長度及類型變化

利用上述數(shù)據(jù)與方法，提取了1997、2002、2007、2012、2017年鹽城海岸線長度及岸線類型信息(圖3)。

圖3 1997-2017年岸線分布Fig.3 Coastline distribution from 1997 to 2017

提取結(jié)果表明，鹽城海岸線由人工岸線和自然岸線組成。自然岸線以粉砂淤泥質(zhì)岸線為主，人工岸線以港口、養(yǎng)殖圍堤和建設(shè)圍堤為主。20年間鹽城海岸線長度整體變化不大，但岸線類型和位置變化較大。1997年岸線類型以粉砂淤泥質(zhì)岸線為主，隨著岸線人工開發(fā)干預，自然岸線逐步萎縮，自然岸線比重由1997年75%減少到62%，與之對應人工岸線呈加速增長的態(tài)勢，并于2017達到1997—2017年來的最大值156.13 km。從岸線空間分布上來看，人工岸線主要分布在濱海港、射陽港、大豐港、王港河口和條子泥，以及四卯酉河南側(cè)、王港河南側(cè)區(qū)域。

3.2 岸線進退變化及變遷速率

計算結(jié)果如圖4和表1所示。從結(jié)果來看，海岸線向海推進強度總體由北向南逐步增加，不同岸段向海推進的強度不同，平均為93.01 m/a，向海推進最大距離為10780 m，位于條子泥墾區(qū)，向岸退縮最大距離為1710 m，位于雙洋河口附近。岸段A和岸段B位于射陽河口北側(cè)，平均變化率為-0.46 m/a，為侵蝕區(qū)域。岸段C和岸段D，侵蝕淤積交替發(fā)生，呈現(xiàn)過度變化的特征，岸線變化較小，為侵蝕淤長過渡區(qū)域。岸段E向海推進287.12 m/a，為淤積岸段。

從海岸線變遷變化趨勢來看，岸段A侵蝕速率越來越小，由1997—2002年-5.44 m/a的侵蝕轉(zhuǎn)為2012—2017年12.23 m/a的弱淤積；岸段B侵蝕越來越強，由1997—2002年淤積 28.14m/a轉(zhuǎn)為2012—2017年-26.25 m/a的侵蝕；岸段C侵蝕越來越強，由1997—2002年淤積86.96 m/a到2012—2017年-7.76 m/a的侵蝕。岸段D同樣侵蝕越來越強，由1997—2002年淤積243.71 m/a到2012—2017年-117.19 m/a的侵蝕。岸段E淤積較強，但是幅度在減小，由1997—2002年384.81 m/a下降到2012—2017年227.81 m/a，下降了40.8%。

分析1997—2017年海岸線斷面變化情況， 2002—2007年岸線向海推進65.82 m/a，2012—2017年海岸線變化率為-7.76 m/a，在斷面號56處(圖4)海岸線變化轉(zhuǎn)為負數(shù)，據(jù)此，鹽城海岸南北沖淤的拐點南移至射陽河口以南7 km射陽鹽場附近。

總之，鹽城海岸整體穩(wěn)定性在增加，強侵蝕、強淤積岸段在減少；部分岸段采取了人工防護，海岸穩(wěn)定。廢黃河口北側(cè)侵蝕速率趨緩，廢黃河口南側(cè)到扁擔港之間侵蝕加劇了，雙洋港到射陽河口岸段淤積逐漸減緩。射陽河口往南岸線整體向海前進，但是前進幅度在逐年減小。

注：“-”表示侵蝕，“+”表示淤積

4 岸線變遷影響因素

4.1 自然因素

岸線變化的自然因素主要是泥沙供給和水動力沖刷，海岸的侵蝕淤積是泥沙供給與動力沖刷之間平衡的結(jié)果，泥沙供給能補償動力沖刷的泥沙損失，岸線便處于穩(wěn)定，否則會海岸發(fā)生會侵蝕淤積。鹽城廢黃河口海岸原為古黃河入?？冢?855年黃河北歸，陸域供沙斷絕，沖淤平衡被打破，海岸在潮汐、波浪和風暴潮的沖刷導致缺沙的海岸發(fā)生侵蝕。

目前鹽城海岸線自然因素為主引起的變化主要位于廢黃河口到射陽河口岸段和新洋港至龍港岸段。海岸侵蝕最嚴重的區(qū)域位于廢黃河口南側(cè)的雙洋港和扁擔河口，該岸段雖然進行了養(yǎng)殖圍堤，但受波浪風暴潮影響，大量圍堤被沖垮，海岸線不斷改變后退，侵蝕造成的最大后退距離達到1700 m。扁擔河口北側(cè)(圖5a)，2017年相對1997年侵蝕后退750 m。雙洋港北側(cè)(圖5b)1997—2017年侵蝕導致岸線變化十分劇烈，2017年相比1997年岸線后退187 m；雙洋港南側(cè)(圖5c)2017年相對2012年侵蝕后退220 m。位于鹽城國家級珍禽自然保護區(qū)的新洋港至斗龍港岸段，由于禁止任何形式的建設(shè)開發(fā)，主要是自然因素引起的海岸線來回擺動，擺動幅度在-130～130 m之間。

圖5 典型岸段侵蝕嚴重年份的海岸線變化

4.2 人為因素

岸線變化的人為因素表現(xiàn)為農(nóng)業(yè)圍墾、港口與海岸防護工程的建設(shè)等人工構(gòu)筑物建設(shè)導致的岸線向海前進和向陸后退。2017年鹽城海岸線人工岸線的比例達到38%，人工岸線類型以港口岸線，圍墾和防波堤為主。鹽城海岸的人類活動以射陽河口為界，北側(cè)以海岸防護工程為主，南側(cè)是以農(nóng)業(yè)圍墾為主，還包括濱海港、射陽港、大豐港為中心港口建設(shè)岸線利用為主導致岸線的人工化。

廢黃河口北側(cè)修建了濱海防波堤，岸線保持穩(wěn)定，比如六合莊岸段護坡工程建設(shè)[19]，歷經(jīng)多次臺風、風暴潮，岸線仍然穩(wěn)定。港口建設(shè)主要為濱海港、射陽港、大豐港港口的建設(shè)導致海岸線持續(xù)向海前進，比如大豐港建設(shè)導致1997—2017年岸線向海推進9.5 km，陸域面積增加59.09 km2(圖6)。

此外，由于鹽城南側(cè)為淤積型海岸，為圍墾提供了有利條件，圍墾加速了圍墾區(qū)潮灘淤積速率，改變了潮灘剖面形態(tài)，導致海岸線格局發(fā)生變化[20]。2010年由于受《江蘇沿海灘涂圍墾開發(fā)利用規(guī)劃綱要》規(guī)劃的270萬畝圍墾計劃的影響，該時期為填海造地最快的一個時期，工程建設(shè)導致人工岸線的快速向海推進和岸線的平直化。比如條子泥圍墾工程自2002年開始建設(shè)，圍墾增加面積增加到2017年的125 km2，自然岸線轉(zhuǎn)為人工岸線，人工岸線增加了22 km，岸線向海推進最大10.78 km，圍墾速度達到6.85 km2/a(圖6)。

圖6 1997—2017年大豐港、梁垛河、條子泥圍墾情況

5 結(jié)論

利用1997、2002、2007、2012和2017年遙感影像數(shù)據(jù)提取了海岸線，結(jié)合數(shù)字海岸線分析系統(tǒng)，利用EPR模型定量計算了海岸線變遷速率，分析了岸線變化的驅(qū)動力。結(jié)論如下：

(1)1997到2017年，鹽城海岸線長度變化不大，但位置變化較大。20年來鹽城海岸線人工化程度在逐步增強。人工岸線由1997年264.10 km增加到2017年303.10 km。自然岸線以粉砂淤泥質(zhì)為主，人工岸線以港口岸線、養(yǎng)殖圍堤岸線為主，主要分布在濱海港、射陽港、大豐港、王港河口和條子泥，以及四卯酉河南側(cè)、王港河南側(cè)區(qū)域。

(2)鹽城海岸線向岸后退與向海推進并存。向岸后退岸段位于灌河口到射陽河口之間，以廢黃河口為中心，南到射陽河口，20年間平均變化率分別為-0.56 m/a和-0.35 m/a，廢黃河北側(cè)海岸侵蝕趨緩，南側(cè)侵蝕加劇。向海推進岸段分布于射陽河口南側(cè)，20年間平均變化率為140.6 m/a，每5年變化率為305.16 m/a、160.79 m/a、121.02 m/a和102.86 m/a，海岸淤長速度在減緩，岸線向海推進幅度在逐年減小，南北沖淤拐點向南擴展至射陽鹽場附近。

(3)鹽城海岸線變化是人為和自然因素的共同影響的結(jié)果。以射陽河口為界，泥沙和動力作用導致的海岸侵蝕是射陽河口北側(cè)岸線向岸退縮的主要原因，農(nóng)業(yè)圍墾、港口建設(shè)和海岸防護工程的建設(shè)，是岸線人工化和向海推進的主要原因。