周相君，李曉敏，馬 毅，吳培強(qiáng)

（國(guó)家海洋局 第一海洋研究所，山東 青島 266061）

0 引言

海岸線標(biāo)識(shí)了水陸分界線，不僅如此，其附近地理、生態(tài)區(qū)位特殊，是海洋經(jīng)濟(jì)和社會(huì)發(fā)展的重要依托。中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)《海洋學(xué)術(shù)語(yǔ)·海洋地質(zhì)學(xué)》（GB／T18190－2000）給出的海岸線定義是“海岸線是海陸分界線，在我國(guó)系指多年大潮高潮位時(shí)的海陸界線”［1］?！段覈?guó)近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)海岸線修測(cè)技術(shù)規(guī)程》［2］中定義的是“海岸線為平均大潮高潮時(shí)水陸分界的痕跡線”。海岸線的變化直接引起海岸帶環(huán)境的變化，與人類的生活、生產(chǎn)活動(dòng)息息相關(guān)。因此，進(jìn)行系統(tǒng)的海岸線動(dòng)態(tài)監(jiān)視監(jiān)測(cè)與研究，增強(qiáng)對(duì)海岸線變化機(jī)制深入的了解十分必要。

遙感技術(shù)具有大面積、同步觀測(cè)、時(shí)效性強(qiáng)等優(yōu)勢(shì)，已廣泛應(yīng)用到海岸線變遷研究中，國(guó)內(nèi)外在這方面的研究甚多，陳正華 等［3］曾經(jīng)利用4期衛(wèi)星資料研究了1986—2009年浙江省大陸的海岸線變遷情況，并進(jìn)行了岸線分維數(shù)分析；李猷 等［4］分析了深圳市1978—2005年海岸線時(shí)空動(dòng)態(tài)演變特征，并探討了其驅(qū)動(dòng)因素；劉鑫［5］基于RS和GIS技術(shù)提取海岸線，研究了鐵山港地區(qū)1987—2006年的海岸線變化特征；黃鵠 等［6］基于遙感和GIS平臺(tái)全面探討了廣西海岸近幾十年來(lái)岸線的時(shí)空演變特征；SHEIK et al［7］利用 DSAS（Digital Shoreline Analysis System，數(shù)字岸線分析系統(tǒng)）分析了印度Kanyakumari和Tuticorin之間的岸線變化，依據(jù)計(jì)算結(jié)果分區(qū)域分別分析了岸線淤蝕情況；DEWIDAR［8］用2種技術(shù)計(jì)算海岸線蝕退，得出EPR（End Point Rate，端點(diǎn)變化率），研究了埃及紅海沿岸MarsaAlam和Hamata兩市海岸線在自然因素作用下的位置變化。趙宗澤 等［9］利用基線法和面積法對(duì)湄洲灣的海岸線變遷情況進(jìn)行了分析，并計(jì)算了反映岸線曲折度的岸線曲率。曾慶留［10］利用1986—2011年遙感影像，提取了5期山東省大陸岸線，從岸線的位置、長(zhǎng)度和類型三方面對(duì)山東省大陸岸線變遷進(jìn)行分析。

隨著人口增長(zhǎng)、城市化和工業(yè)化的迅速發(fā)展以及對(duì)外貿(mào)易的日漸加強(qiáng)，防城灣沿岸建設(shè)需求加大，從20世紀(jì)70年代以來(lái)，特別是1987年5月南防鐵路通車后，加強(qiáng)了漁澫島與大陸的聯(lián)系，位于島上的防城港進(jìn)入加速建設(shè)時(shí)期；1993年5月，國(guó)務(wù)院正式批準(zhǔn)設(shè)立地級(jí)防城港市，防城灣沿岸建設(shè)進(jìn)一步飛速發(fā)展。陸續(xù)地一系列填海造地、海洋工程建設(shè)等人類活動(dòng)都導(dǎo)致海岸線發(fā)生了不同程度的變化。本文使用1973，1979，1990和2000年共4景Landsat和1景2010年的HJ－1B衛(wèi)星遙感影像，采用人機(jī)交互解譯方式結(jié)合GIS軟件提取了防城灣的海岸線，并分析了本區(qū)近40a來(lái)的海岸線變遷特征。

1 研究區(qū)域與研究方法

1.1 研究區(qū)概況

防城灣地處廣西壯族自治區(qū)沿海西部、北部灣西北岸，西起江山半島的沙尾咀附近，東至企沙半島的炮臺(tái)附近，范圍為21°32′30″～21°43′00″N，108°17′30″～108°28′35″E。防城灣灣口朝南，北、東、西三面為丘陵地，形勢(shì)隱蔽，水域平靜，防風(fēng)條件較好。我國(guó)西南第一大港——防城港位于灣內(nèi)漁澫島之上。從防城港港口出發(fā)，海路東距北海港62海里，東南距?？诟?63海里，西南距越南海防港151海里。港灣深入內(nèi)陸達(dá)11～14km，陸域?qū)拸V，是中國(guó)大陸海岸線最南端的深水海港，是全國(guó)20個(gè)樞紐港之一［11］。整個(gè)防城灣被漁澫島分為兩部分，以西為防城江口，以東為暗埠口江。因此，將研究區(qū)分為3個(gè)岸段：防城江口沿岸、暗埠口江沿岸以及漁澫島，其區(qū)位示意圖如圖1所示。

1.2 數(shù)據(jù)處理

本文以防城灣1973，1979，1990，2000與2010年共5個(gè)時(shí)期的遙感影像為主要數(shù)據(jù)源，其中有2景MSS影像，2景TM影像，1景HJ－1B衛(wèi)星遙感影像（表1）。以“908”專項(xiàng)海島岸線專題數(shù)據(jù)及1景2007年的ALOS正射校正影像為輔助數(shù)據(jù)源，其中海島岸線專題數(shù)據(jù)主要用于對(duì)岸線提取結(jié)果進(jìn)行精度評(píng)價(jià)；ALOS正射校正影像主要用于對(duì)5期遙感影像進(jìn)行配準(zhǔn)。

圖1 研究區(qū)位置示意圖Fig.1 Location of study area

表1 遙感影像表Tab.1 List of the remote sensing images

為了保證海岸線提取位置的準(zhǔn)確性和各期岸線的可對(duì)比性，對(duì)各期影像進(jìn)行了幾何校正。以2007年的ALOS正射校正影像為基準(zhǔn)，選擇無(wú)明顯變化的、易從影像上識(shí)別的路口、基巖島、碼頭等地物作為控制點(diǎn)，采用二次多項(xiàng)式變換對(duì)所用的5期遙感影像進(jìn)行地理配準(zhǔn)，效果如圖2所示。其中，1973和1979年影像進(jìn)行幾何校正時(shí)直接利用雙線性內(nèi)插法作30m重采樣以統(tǒng)一影像的空間分辨率，所有影像的投影和坐標(biāo)系均采用高斯－克呂格投影和WGS－84坐標(biāo)系。

圖2 影像配準(zhǔn)效果（上：ALOS衛(wèi)星遙感影像，左下：TM衛(wèi)星遙感影像，右下：HJ－1B衛(wèi)星遙感影像）Fig.2 Registration effect of remote sensing images（upper：ALOS image；lower left：TM image；lower right：HJ－1Bimage）

為了能更好地從影像上解譯海岸線，本文對(duì)不同影像開(kāi)展了相應(yīng)的增強(qiáng)處理，使用了不同波段進(jìn)行組合，組合的方式為：MSS影像選取6、5、4波段分別賦予紅、綠、藍(lán)色制作假彩色合成圖，對(duì)于水陸分界的岸線解譯較為敏感；TM影像采用4、3、2波段合成假彩色圖像，水陸邊界明顯；至于HJ－1B衛(wèi)星遙感影像，則選取4、3、2波段，并利用公式（1）［12］的計(jì)算結(jié)果替換第3波段以增強(qiáng)色彩效果：

1.3 海岸線提取

基于現(xiàn)場(chǎng)踏勘和當(dāng)?shù)氐膶?shí)際情況，本文將防城灣海岸線分為4類：基巖岸線、淤泥質(zhì)岸線、生物岸線和人工岸線，參考《我國(guó)近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)海岸線修測(cè)技術(shù)規(guī)程》［2］中的界定原則，結(jié)合各類海岸線的地貌特征，確定各類型海岸線的位置?；鶐r岸線確定在基巖質(zhì)陡崖的基部，如圖3a；淤泥海岸潮間帶上常有一條耐鹽植物生長(zhǎng)狀況明顯變化的界線，以此作為淤泥質(zhì)岸線的位置，如圖3b；生物岸線在本區(qū)域主要指紅樹(shù)林岸線，其位置確定在紅樹(shù)林內(nèi)邊界上，海岸線以下為紅樹(shù)林，以上為陸地［13］，如圖3c；本區(qū)人工構(gòu)筑物向陸一側(cè)皆不存在平均大潮高潮時(shí)海水能達(dá)到水域，所以以永久性人工構(gòu)筑物的向海側(cè)的平均大潮時(shí)水陸分界的痕跡線作為人工岸線的位置，如圖3d。

圖3 海岸線的遙感解譯標(biāo)志Fig.3 Remote sensing interpretation signs of shoreline

基于2010年的HJ－1B衛(wèi)星遙感影像，利用所建立的各類型海岸線遙感解譯標(biāo)志，對(duì)2010年的海岸線進(jìn)行信息提取。利用“908”專項(xiàng)海島岸線專題數(shù)據(jù)對(duì)海岸線提取結(jié)果進(jìn)行精度評(píng)價(jià)，其誤差在0.5～1.0個(gè)像元以內(nèi)，評(píng)價(jià)結(jié)果表明本文所建立的各類型海岸線遙感解譯標(biāo)志可行。基于此，對(duì)另外4期的海岸線以人機(jī)交互解譯方式統(tǒng)一在1∶10 000窗口下進(jìn)行提取，結(jié)果見(jiàn)圖4，海岸線長(zhǎng)度統(tǒng)計(jì)情況見(jiàn)圖5。為保證相鄰時(shí)相兩期岸線沒(méi)有變化的位置及屬性信息保持嚴(yán)格一致，每期海岸線在相鄰時(shí)相的海岸線提取結(jié)果基礎(chǔ)上進(jìn)行修邊得到。

圖4 各期防城灣海岸線圖Fig.4 Shoreline of the Fangcheng Bay in different periods

圖5 各期防城灣海岸線長(zhǎng)度統(tǒng)計(jì)Fig.5 Statistics of shoreline length of the Fangcheng Bay in different periods

2 海岸線變遷分析

通過(guò)5期衛(wèi)星遙感影像解譯結(jié)果疊加可以看出：1973—2010年間防城灣海岸線變遷情況較為復(fù)雜，近40a來(lái)沿岸不斷有新增面積出現(xiàn)，共增加了約4 216.61hm2，其中1973－1990年增加了736.64 hm2，1990－2010年增加了3 479.76hm2（圖6）；海岸線總長(zhǎng)度則持續(xù)遞減，只在2010年稍有回升，與1973年相比減少了36.96km（圖7）。

圖6 防城灣增長(zhǎng)面積變化趨勢(shì)圖Fig.6 Changing tendency of increased area along the Fangcheng Bay

圖7 防城灣海岸線長(zhǎng)度變化趨勢(shì)圖Fig.7 Changing tendency of shoreline length along the Fangcheng Bay

從空間上看，漁澫島的海岸線變遷程度最為劇烈，新增長(zhǎng)的岸段及部分原始岸段皆轉(zhuǎn)為人工岸線類型，而基巖與淤泥質(zhì)岸線類型幾乎消失殆盡；其次是暗埠口江沿岸，以淤泥質(zhì)為主，近40a間，海岸不斷向海推進(jìn)，許多原本曲折多變的鹿角灣被逐漸填充，岸線變得日趨平直；防城江口沿岸則相對(duì)變化不大，各類型岸線變化較為穩(wěn)定（圖5）。

從時(shí)間上看，防城灣1973—1979年海岸線變遷程度很小，1979—1990年開(kāi)始加速，1990—2000年以及2000—2010年2個(gè)10a變化最為劇烈（圖5），究其原因人為因素占主導(dǎo)。

2.1 防城江口沿岸海岸線變遷

防城江口沿岸海岸線在1973—2010年近40a間岸線變遷的程度很小，變遷總面積為74.58hm2。其中，1973—1990年間提取的3期岸線幾乎完全重疊（圖4）。2000年時(shí)，僅在防城江河口附近有人工圍填以加固河堤的跡象，圍填面積約為36.75hm2，此處岸線長(zhǎng)度比1990年減少1.35km，而新增人工岸線0.5km；2000年海岸線長(zhǎng)度比1990年縮短了2.09 km。至2010年，繼續(xù)在2000年基礎(chǔ)上擴(kuò)大了河口附近的圍填面積，又增加了約37.83hm2，海岸線總長(zhǎng)度則變化不大，引起這種變化的原因主要是人工圍填使得深入內(nèi)陸的一些較小的溺谷灣消失，之前曲折的岸線變得更加平直；另外，灣口沙尾咀處向海推進(jìn)了約0.1km。

2.2 暗埠口江沿岸海岸線變遷

暗埠口江沿岸海岸線類型以淤泥質(zhì)為主，近40a間，本區(qū)海岸不斷向海推進(jìn)，許多原本曲折多變的鹿角灣被逐漸填充，岸線變得日趨平直，這種變化的速率較為穩(wěn)定，只是在2010年時(shí)，由于防城港企沙港區(qū)的修建，出現(xiàn)大面積的圍海造陸活動(dòng)。

總體來(lái)看，1973—1979年間，本區(qū)海岸線無(wú)太大變化，較明顯的區(qū)域只在風(fēng)流嶺江灣口南岸有2個(gè)較小的海灣以及內(nèi)側(cè)北岸一個(gè)海灣淤積成陸，面積合約51.45hm2；海岸線總長(zhǎng)度則由1973年的126.59km縮短到1979年的121.55km。

1979—1990年，人工岸線突增20.24km。云約江沿岸有3個(gè)港灣消失，其中南岸2個(gè)港灣是由人工圍海養(yǎng)殖造成，面積合約170.18hm2，岸線長(zhǎng)度共減少3.29km；東岸1個(gè)港灣則是自然淤積成陸，面積為33.39hm2，岸線長(zhǎng)度減少1.22km。風(fēng)流嶺江兩岸也有不同程度的人工圍填的痕跡，面積總計(jì)104.04hm2，岸線長(zhǎng)度則減少了6.03km。公車附近養(yǎng)殖堤壩的修建更是與灣內(nèi)2個(gè)海島6451和6452相接，圍出了32.44hm2的養(yǎng)殖池塘，使得岸線長(zhǎng)度減少了1.03km。

2000年與1990年相比，本區(qū)海岸線變遷更甚。主要表現(xiàn)為，淤泥岸段的淤長(zhǎng)以及人工圍填用以養(yǎng)殖或工程、道路建設(shè)。從潭頭到公車段的岸線幾乎都有向海推進(jìn)，風(fēng)流嶺江南岸向海最大推進(jìn)了近1km 云約江南岸新建了84.95hm2的養(yǎng)殖池塘，東岸繼續(xù)淤積，面積為22.14hm2，岸線長(zhǎng)度共縮減2.01km；另外，潭頭附近也有新增的大面積淤積區(qū)，約達(dá)119.65 hm2，岸線長(zhǎng)度縮短了1.9km。

至2010年，海岸線變遷態(tài)勢(shì)還是以向海淤擴(kuò)為主，但是程度明顯不及前1個(gè)10a。大墩咀與潭頭附近分別有22.09hm2與9.54hm2的人工圍填，岸線長(zhǎng)度合計(jì)增加0.8km；云約江東岸持續(xù)淤積，新增面積51.23hm2，岸線長(zhǎng)度縮減1.89km；松柏港附近由于紅樹(shù)林保護(hù)加強(qiáng)，促使沿岸加速淤積，凈增面積12.6hm2，岸線長(zhǎng)度減少0.29km；企沙大道的修建將風(fēng)流嶺江灣底截?cái)?，使得岸線縮減3.51km，靠陸一側(cè)面積約42.88hm2，已基本被圍填。這一時(shí)期，本區(qū)變化最大的當(dāng)屬炮臺(tái)角外圍的人工圍填區(qū)域，面積達(dá)940.19hm2，用作防城港企沙港區(qū)的興建。

2.3 漁澫島海岸線變化

1973—2010年，漁澫島面積同樣呈不斷增加趨勢(shì)（圖8），近40a來(lái)共增加2 027.41hm2，總體增長(zhǎng)了近2倍（圖9），但增長(zhǎng)率有變化；海岸線總長(zhǎng)度則由46.49km變?yōu)?4.03km，縮短了2.46km，但其趨勢(shì)不是遞減的，期間有增長(zhǎng)的“抬升”期?？梢钥闯?，整個(gè)研究時(shí)段中可分為2個(gè)不同的岸線變遷時(shí)期，分別是1973—1990年和1990—2010年（圖10）。

圖8 漁澫島海岸線變遷圖Fig.8 Shoreline change on Yuwan Island

圖9 漁澫島增長(zhǎng)面積變化趨勢(shì)圖Fig.9 Changing tendency of increased area on Yuwan Island

圖10 漁澫島海岸線長(zhǎng)度變化趨勢(shì)圖Fig.10 Changing tendency of shoreline length on Yuwan Island

從海岸線類型來(lái)看，1973年以基巖及淤泥質(zhì)岸線為主，約占97.3%，人工岸線很少；到2010年，人工岸線與生物岸線占絕對(duì)優(yōu)勢(shì)，約為92.4%，新增長(zhǎng)的岸段及部分原始岸段皆轉(zhuǎn)為人工岸線類型，而基巖與淤泥質(zhì)岸線類型幾乎消失殆盡（圖5）。

漁澫島海岸線在1973—1990年間和1990—2010年間變化趨勢(shì)明顯不同。1973—1990年間，漁澫島面積逐漸增長(zhǎng)，但增速較緩，岸線總長(zhǎng)度則呈不斷遞減趨勢(shì)；1990—2010年間，漁澫島面積擴(kuò)速陡增，岸線總長(zhǎng)度也轉(zhuǎn)而遞增。

2.3.1 1973 —1990年漁澫島海岸線變遷

1973—1990年間，漁澫島面積增長(zhǎng)速率較小，約為18.82hm2／a，這主要是由于該時(shí)期港區(qū)雖然在建，但沒(méi)有太大規(guī)模的圍填工程，且速度較為穩(wěn)定，到1986年為止，防城港第一期工程只有7個(gè)萬(wàn)噸級(jí)深水泊位峻工投入使用。

1973—1979年間，整個(gè)漁澫島總體面積變化不大，僅增加了約94.63hm 岸線長(zhǎng)度縮減2.87km此期間處于防城港建設(shè)開(kāi)工初期，西南端部分區(qū)域的圍海造陸占據(jù)了較大變化份額，形成了防城港的雛形。另外，漁澫島北端有部分小海灣也被圍填。

1979—1990年間，隨著進(jìn)一步的人工圍填，本區(qū)面積增加了225.25hm2，岸線總計(jì)縮短2.32km，這主要是由于東南端圍海建造碼頭。北端的所有海灣也基本被圍填，若干近岸小島（茶壺墩、捕魚(yú)嶺等）都被圍入人工岸線以內(nèi)（圖8），岸線曲率變小，原本曲折的自然岸線日趨平直，使得海岸線總長(zhǎng)度呈縮減狀態(tài)，平均縮減幅度約為0.63km／a。1987年南防鐵路建成，漁澫島與大陸聯(lián)系加強(qiáng)，鐵路通過(guò)島體的中段東岸，使其也有小規(guī)模外擴(kuò)。

2.3.2 1990 —2010年漁澫島海岸線變遷

1990—2010年間，受政策影響，漁澫島上建設(shè)提速，面積不斷擴(kuò)大，年均增長(zhǎng)約85.38hm2；更多原始的自然岸線轉(zhuǎn)為人工岸線，海岸線總長(zhǎng)度轉(zhuǎn)而增加，2010年達(dá)到44.03km，比1990年時(shí)的最小值增長(zhǎng)了約6.9km，增長(zhǎng)速率約為0.35km／a。

1993年，防城港被批準(zhǔn)正式設(shè)立地級(jí)市，同時(shí)借助于國(guó)務(wù)院給予的享受沿海開(kāi)放城市的優(yōu)惠政策，島上的基礎(chǔ)設(shè)施建設(shè)陡然提速，總面積于1990—2000年間增長(zhǎng)了約772.89hm2，岸線長(zhǎng)度增加4.62km；由于防欽高速公路的修建，本區(qū)北端再一次與大陸相通，西部沿岸大面積海域被圍入人工岸線以內(nèi)，使得長(zhǎng)山尾島也與漁澫島貫連；期間，東岸及南岸分別有262hm2與79hm2的新增碼頭泊位。2000—2010年間，防城港擴(kuò)建勢(shì)頭不減，新增面積達(dá)934.63hm2，新增岸線長(zhǎng)達(dá)5.84km，集中在南端新建港區(qū)，5萬(wàn)、7萬(wàn)、10萬(wàn)甚至20萬(wàn)噸級(jí)泊位相繼建成投產(chǎn)（圖11）。至此，防城港已有6個(gè)港區(qū)，漁澫島總面積達(dá)到1973年時(shí)的3倍。

2.3.3 漁澫島分維數(shù)計(jì)算

近年來(lái)越來(lái)越多的學(xué)者嘗試?yán)枚糠椒ɑ騾?shù)來(lái)描述與分析海岸線變遷特征，其中，海岸線分維數(shù)是表述海岸線變遷的一種參數(shù)。

采用網(wǎng)格法計(jì)算海岸線的分維數(shù)。根據(jù)分形理論Nk（εk）∝ε－Dk，則D就是被測(cè)海岸線的分維數(shù)。按照網(wǎng)格長(zhǎng)度ε1，ε2，ε3，…，εk，分別統(tǒng)計(jì)被覆蓋海岸線的網(wǎng)格數(shù)目N（ε1），N（ε2），N（ε3），…，N（εk），對(duì)公式

進(jìn)行回歸（圖12），可計(jì)算出海岸線分維數(shù)D，式（2）中A為待定常數(shù)［14］。

首先在ArcGIS中將已提取的1973—2010年5期漁澫島海岸線矢量層，按照48163264128和256m分別進(jìn)行柵格化，獲得不同尺度下漁澫島海岸線長(zhǎng)度，尺度越大海岸線越短［15］。然后統(tǒng)計(jì)每個(gè)空間分辨率網(wǎng)格下海岸線的網(wǎng)格數(shù)量，依照分維數(shù)計(jì)算公式對(duì)每組數(shù)據(jù)進(jìn)行回歸分析，從而得出時(shí)間序列上漁澫島海岸線的分維數(shù)D（表2）。

圖11 衛(wèi)星圖像顯示的漁澫島南端排尾附近海岸變遷情況Fig.11 Coast change nearby Paiwei on south of Yuwan Island displayed on remote sensing images

表2 漁澫島不同時(shí)期海岸線分維數(shù)計(jì)算結(jié)果Tab.2 Fractal dimensions of shoreline in different periods of Yuwan Island

由表2可知，漁澫島分形維數(shù)范圍介于1.023～1.062之間。1973年分形維數(shù)為最大值1.062，表明當(dāng)時(shí)海岸線結(jié)構(gòu)復(fù)雜，人工改造的痕跡較小，接近于自然狀態(tài)。從1973年到2010年呈遞減趨勢(shì)，表明海岸結(jié)構(gòu)的復(fù)雜度降低。1979年與1990年分維數(shù)差別較大，從漁澫島海岸線變遷圖（圖8）可以看出，之前一些較小的港灣到1990年已經(jīng)圍填，岸線曲率變小，總體趨于平緩。1990年與2000年分維數(shù)又相差較大，原因從圖5可以看出，到2000年為止，海岸線幾乎已全部為人工岸線，變得更為平直。而2010年的分維數(shù)為最小值，表明海岸線結(jié)構(gòu)最簡(jiǎn)單，由人工改造的幅度更大，變?yōu)楦?guī)則、平直的幾何形狀。

圖12 分維數(shù)回歸曲線Fig.12 Regression curves of fractal dimensions

3 結(jié)論

防城灣海岸線在1973—2010年間變化較為明顯，本文利用Landsat與HJ－1B衛(wèi)星遙感影像，通過(guò)人機(jī)交互解譯方式提取海岸線，對(duì)防城灣海岸線動(dòng)態(tài)變化進(jìn)行了定量分析，結(jié)果表明：

（1）近40a來(lái)，防城灣海岸線曲折程度不斷減小，日趨平直化，總長(zhǎng)度呈遞減趨勢(shì)，僅在2010年時(shí)稍有回升。不同時(shí)期陸域面積增加的幅度不同：1973—1990年，政策平穩(wěn)，增幅較緩，僅增加了736.64hm2；1990年后，受改革開(kāi)放政策驅(qū)動(dòng)，擴(kuò)張速度陡增，增加了3 479.96hm2。研究區(qū)變化幅度最大的區(qū)域?yàn)闈O澫島西岸和南岸以及企沙半島的炮臺(tái)沿岸，均為防城港港區(qū)的所在；暗埠口江沿岸變化次之；防城江口海岸線變遷幅度最小。

（2）防城江口沿岸以基巖岸線為主，海岸線在1973—2010年間變遷的程度很小，僅在防城江河口沿岸增加了74.58hm2的面積，引起這種變化的原因主要是人工圍填固岸使得深入內(nèi)陸的一些較小的溺谷灣消失。

（3）暗埠口江沿岸以淤泥岸線為主，近40a間不斷向海推進(jìn)，許多原本曲折多變的鹿角灣被逐漸填充，伴隨著養(yǎng)殖業(yè)的發(fā)展，人為修筑堤壩，岸線變得日趨平直，這種變化的速率較為穩(wěn)定。2010年時(shí)由于防城港企沙港區(qū)的修建，出現(xiàn)大面積的圍海造陸活動(dòng)。

（4）1973—2010年，漁澫島面積同樣呈不斷增加趨勢(shì)，近40a來(lái)共增加2 027.41hm2，總體增長(zhǎng)了近2倍，占整個(gè)研究區(qū)陸域增長(zhǎng)面積的48%；海岸線總長(zhǎng)度則由1973年的46.49km減少到2010年的44.03km。漁澫島海岸線在1973—1990年間及1990—2010年間變化趨勢(shì)明顯不同。1973—1990年間，漁澫島面積逐漸增長(zhǎng)，但增速較緩，岸線總長(zhǎng)度則呈不斷遞減趨勢(shì)；1990—2010年間，漁澫島面積擴(kuò)速陡增，岸線總長(zhǎng)度也轉(zhuǎn)而遞增。

本文所用影像的空間分辨率皆不高于30m，這對(duì)分析結(jié)果可能有一定影響，造成與實(shí)際情況有些偏差，今后將嘗試?yán)酶呖臻g分辨率的遙感影像進(jìn)行分析，以獲取更精確的結(jié)果。

致謝 本文所用的Landsat衛(wèi)星遙感影像從NASA網(wǎng)站免費(fèi)下載；所用的環(huán)境數(shù)據(jù)由中國(guó)資源衛(wèi)星應(yīng)用中心免費(fèi)提供，謹(jǐn)致謝忱！

（References）：

［1］GB／T18190—2000Oceanological terminology：Marine geology［S］.2000.

GB／T18190—2000海洋學(xué)術(shù)語(yǔ)·海洋地質(zhì)學(xué)［S］.2000.

［2］908Project Office of State Oceanic Administration.Chinese off－shore investigation and assessment technical regulation of coastline measurement：Trial implementation［M］.Beijing：Ocean Press，2007.

國(guó)家海洋局908專項(xiàng)辦公室.我國(guó)近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)海岸線修測(cè)技術(shù)規(guī)程：試行本［M］.北京：海洋出版社，2007.

［3］CHEN Zheng－h(huán)ua，MAO Zhi－h(huán)ua，CHEN Jian－yu.Coastline change monitoring using 4periods remote sensing data in Zhejiang Province from 1986to 2009［J］.Remote Sensing Technology and Application，2011，26（1）：68－73.

陳正華，毛志華，陳建裕.利用4期衛(wèi)星資料監(jiān)測(cè)1986～2009年浙江省大陸海岸線變遷［J］.遙感技術(shù)與應(yīng)用，2011，26（1）：68－73.

［4］LI You，WANG Yang－lin，PENG Jian，et al.Research on dynamic changes of coastline in Shenzhen City based on Landsat image［J］.Resources Science，2009，31（5）：875－883.

李猷，王仰麟，彭建，等.深圳市1978年至2005年海岸線的動(dòng)態(tài)演變分析［J］.資源科學(xué)，2009，31（5）：875－883.

［5］LIU Xin.Analyzing shoreline changes of Tieshan Port in Guangxi Province using remote sensing technology ［J］.Geospatial Information，2012，10（1）：102－106.

劉鑫.應(yīng)用遙感方法的廣西鐵山港區(qū)海岸線變遷分析［J］.地理空間信息，2012，10（1）：102－106.

［6］HUANG Hu，HU Zi－ning，CHEN Xin－geng，et al.Analyses on spatial and temporal changes of Guangxi shoreline based on remote sensing and GIS［J］.Journal of Tropical Oceanography，2006，25（1）：66－70.

黃鵠，胡自寧，陳新庚，等.基于遙感和GIS相結(jié)合的廣西海岸線時(shí)空變化特征分析［J］.熱帶海洋學(xué)報(bào)，2006，25（1）：66－70.

［7］SHEIK M，CHANDRASEKAR.A shoreline change analysis along the coast between Kanyakumari and Tuticorin，India，using digital shoreline analysis system ［J］.Arabia Journal of Geoscience，2011，14（4）：282－293.

［8］DEWIDAR K.Changes in the shoreline position caused by natural processes for coastline of Marsa Alam and Hamata，Red Sea，E－gypt［J］.International Journal of Geosciences，2011，2：523－529.

［9］ZHAO Zong－ze，LIU Rong－jie，MA Yi，et al.Remote sensing monitoring and analysis of coastline changes in the Meizhou Bay since the last 30years［J］.Coastal Engineering，2013，32（1）：19－27.

趙宗澤，劉榮杰，馬毅，等.近30年來(lái)湄洲灣海岸線變遷遙感監(jiān)測(cè)與分析［J］.海岸工程，2013，32（1）：19－27.

［10］ZENG Qing－liu.Shoreline change analysis of Shandong Province on remote sensing since 1986［D］.Qingdao：Shandong University of Science and Technology，2012：1－60.

曾慶留.1986年以來(lái)山東省海岸線變遷遙感分析［D］.青島：山東科技大學(xué)，2012：1－60.

［11］ The Navigation Guarantee Department of the Chinese Navy Headquarters.Guide to Chinese ports：South China Sea［M］.Tianjin：China Navigation Publication Press，2005.

中國(guó)人民解放軍海軍司令部航海保障部.中國(guó)港口指南：南海海區(qū)［M］.天津：中國(guó)航海圖書(shū)出版社，2005.

［12］LI Jun－jie.Resources and environment satellite image processing methods ［EB／OL］.http：//www.cresda.com／n16／n1115／n1522／n149706／n149771.files／n149794.pdf.

李俊杰.資源、環(huán)境衛(wèi)星影像加工處理方法［EB／OL］.http：∥www.cresda.com／n16／n1115／n1522／n149706／n149771.files／n149794.pdf.

［13］SUN Wei－fu，MA Yi，ZHANG Jie，et al.Study of remote sensing interpretation keys and extraction technique of different types of shoreline［J］.Bulletin of Surveying and Mapping，2011，3：41－44.

孫偉富，馬毅，張杰，等.不同類型海岸線遙感解譯標(biāo)志建立和提取方法研究［J］.測(cè)繪通報(bào)，2011，3：41－44.

［14］ZHU Xiao－h(huán)ua，PAN Ya－juan.Study on fractal determination of coast based on GIS［J］.Marine Science Bulletin，2002，21（2）：49－54.

朱曉華，潘亞娟.GIS支持的海岸類型分形判定研究［J］.海洋通報(bào)，2002，21（2）：49－54.

［15］ZHANG Hua－guo，HUANG Wei－gen.Study on spatial scale of shoreline remote sensing information based on fractal theory［J］.Journal of Remote Sensing，2006，10（4）：463－468.

張華國(guó)，黃韋艮.基于分形的海岸線遙感信息空間尺度研究［J］.遙感學(xué)報(bào)，2006，10（4）：463－468.