基于GIS和RS的山東日照海岸帶1988—2018年間演化分析

苗苗，謝小平

(曲阜師范大學(xué)地理與旅游學(xué)院，日照 276826)

0 引言

海岸帶作為海、陸、氣交互作用的地帶，不僅是海洋生態(tài)系統(tǒng)向陸地生態(tài)系統(tǒng)的過(guò)渡區(qū)域，也是全球最重要的生態(tài)交錯(cuò)帶，各種過(guò)程耦合多變，機(jī)制復(fù)雜，疊加人類活動(dòng)，使得海岸帶成為響應(yīng)全球變化最敏感的區(qū)域[1-4]。隨著改革開(kāi)放的發(fā)展，我國(guó)的海岸帶開(kāi)發(fā)加速，海岸帶局部侵蝕問(wèn)題日益嚴(yán)重，海岸侵蝕已經(jīng)與風(fēng)暴潮和巨浪災(zāi)害并列為我國(guó)三大海洋災(zāi)害，海岸帶變化的研究也日趨受到研究學(xué)者們的重視[5-7]。中國(guó)海岸分布主要集中于遼東半島、山東半島和華南海岸3個(gè)區(qū)域，國(guó)內(nèi)學(xué)者對(duì)山東半島海岸變化情況研究較多，但是對(duì)日照海岸帶時(shí)空變化的研究卻很少。山東省日照市作為山東半島藍(lán)色經(jīng)濟(jì)區(qū)的重要組成部分和我國(guó)東部沿海地區(qū)的新興港口城市，近年來(lái)經(jīng)濟(jì)發(fā)展迅速，海岸帶地區(qū)面臨日趨嚴(yán)峻的生態(tài)環(huán)境風(fēng)險(xiǎn)。因此，對(duì)日照海岸帶進(jìn)行系統(tǒng)研究，了解其發(fā)展趨勢(shì)是極其重要的。

遙感(remote sensing,RS)和地理信息系統(tǒng)(geographic information system,GIS)技術(shù)作為新的監(jiān)測(cè)手段在海岸帶動(dòng)態(tài)變化監(jiān)測(cè)上得到了廣泛應(yīng)用[8-9]，國(guó)內(nèi)學(xué)者基于RS和GIS技術(shù)對(duì)我國(guó)自北向南多個(gè)區(qū)域海岸線變化進(jìn)行了研究，如姚曉靜等[10]利用衛(wèi)星遙感資料與現(xiàn)場(chǎng)實(shí)測(cè)資料對(duì)海南島30 a海岸線時(shí)空變化進(jìn)行了研究；陳金月[11]結(jié)合GIS和RS技術(shù)分析了珠江三角洲40 a間海岸線變遷及驅(qū)動(dòng)因素；康波等[12]基于遙感和GIS技術(shù)對(duì)長(zhǎng)島南五島30 a來(lái)海岸線時(shí)空變遷進(jìn)行詳細(xì)分析。近年來(lái)，一些學(xué)者對(duì)日照海岸帶地貌和地質(zhì)景觀、海岸線形態(tài)、海岸帶景觀格局變化過(guò)程進(jìn)行了研究[13-15]，但只是對(duì)日照海岸帶的局部特性做了研究。因此，本文綜合利用RS與GIS技術(shù)，對(duì)日照海岸帶1988—2018年間近30 a的時(shí)空變化特征及其演變的驅(qū)動(dòng)力進(jìn)行了探究。

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

1.1 研究區(qū)概況

山東省日照市海岸帶位于黃海中部，北起日照市東港區(qū)兩城河口，與青島膠南市海岸接壤，南至嵐山區(qū)繡針河入?？?，與江蘇贛榆縣相連，西至204國(guó)道，東至20 m水深線附近(圖1)。沿岸發(fā)育平直的基巖砂礫質(zhì)海岸，中部的奎山、南部的嵐山兩段為剝蝕丘陵，形成突出的岬角濱海陡崖和水下巖灘；南北部的剝蝕平原和現(xiàn)代海岸線之間為寬廣的沙壩潟湖沉積體系。區(qū)域斷裂構(gòu)造以NE向的日照—青島斷裂和NW向的梭羅樹(shù)斷裂為主，海岸呈NE-SW走向。境內(nèi)地表水水體發(fā)育，河流為內(nèi)陸河入海段，流向多為NW-SE向，自北而南主要有白馬河、兩城河、潮白河、付疃河、巨峰河、龍王河和繡針河，除傅疃河等為永久性河流外，其他多屬季節(jié)性河流。由于該地區(qū)地理位置優(yōu)越，形成了眾多的天然良港，海岸線上的主要港灣包括石臼灣、佛手灣與日照港、嵐山港組成日照港群。日照海岸帶是日照經(jīng)濟(jì)、文化發(fā)展中心，現(xiàn)已發(fā)展成集臨海工業(yè)、港口貿(mào)易、商貿(mào)、金融、交通運(yùn)輸、旅游度假等多功能、外向型、綜合性、現(xiàn)代化的海岸帶[16-18]。

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

1.2 數(shù)據(jù)源及其預(yù)處理

本文選取的數(shù)據(jù)包括：研究區(qū)1988年、1998年、2008年和2018年4期遙感影像數(shù)據(jù)，來(lái)自地理空間數(shù)據(jù)云網(wǎng)站(http://www.gscloud.cn/)，其中2018年影像數(shù)據(jù)為L(zhǎng)andsat8 OLI，2008年、1998年和1988年影像數(shù)據(jù)為L(zhǎng)andsat5 TM；研究區(qū)各年潮汐數(shù)據(jù)來(lái)源于中國(guó)港口網(wǎng)(http://www.chinaports.com)；參考數(shù)據(jù)包括日照市行政區(qū)劃、1∶50 000基礎(chǔ)地理信息數(shù)據(jù)、數(shù)字高程模型(digital elevation model,DEM)數(shù)據(jù)，以10 a為間隔來(lái)研究日照海岸帶時(shí)空變化特征。采用ENVI5.2軟件數(shù)據(jù)處理平臺(tái)，對(duì)研究區(qū)4期遙感影像依次進(jìn)行輻射校正、幾何糾正、鑲嵌、裁剪等數(shù)據(jù)預(yù)處理操作，以消除空間上和輻射物理量上的畸變或噪聲，得到研究區(qū)各時(shí)相遙感影像圖[19-20]。

2 研究方法

2.1 海岸線類型的劃分及提取

海岸線是指平均大潮高潮時(shí)水陸分界的痕跡線，本文根據(jù)海岸形態(tài)、成因、物質(zhì)組成和發(fā)展階段等特征將日照市海岸線劃分為基巖海岸線、砂質(zhì)海岸線、淤泥質(zhì)海岸線、人工海岸線4類，使用紅光、近紅外和短波紅外1波段合成遙感影像，根據(jù)不同海岸線的紋理、形狀、色調(diào)等特征建立各海岸線的解譯標(biāo)志(表1)[21]。

表1 日照市海岸線類型Tab.1 Coastline type in Rizhao City

在提取海岸線前，先利用濾波器突出各類型海岸線的水陸分界線[22-23]。采用修正歸一化差異水體指數(shù)(modified normalized difference water index，MNDWI)提取Landsat5遙感影像的水陸分界線，表達(dá)式為：

MNDWI=(G-MIR)/(G+MIR)，

(1)

式中：MIR為中紅外波段的反射率；G為綠光波段的反射率。采用自動(dòng)提取水體指數(shù)(automatic extraction water index，AWEI)提取Landsat8遙感影像的水陸分界線，表達(dá)式為：

AWEI=B+2.5G-1.5(NIR+S1)-0.25S2，

(2)

式中：B為藍(lán)光波段的反射率；NIR為近紅外波段的反射率；S1和S2分別為2個(gè)短波紅外波段的反射率。運(yùn)用式(1)—(2)獲得4個(gè)時(shí)相的水體指數(shù)影像，分別設(shè)定 MNDWI和AWEI的取值范圍來(lái)提取水陸分界線[24-26]。由于衛(wèi)星過(guò)境時(shí)刻獲取的海岸線影像不能保證其為平均大潮高潮線，而海岸線的位置又易受潮汐等因素的影響，故需要對(duì)瞬時(shí)水邊線進(jìn)行潮位校正[27-29]。人工海岸和基巖海岸水陸分明、坡度大而陡，潮汐影響小，不需要做潮位校正，僅對(duì)砂質(zhì)海岸和淤泥質(zhì)海岸做潮位校正。本文采用外牽引法[30]進(jìn)行潮位校正，根據(jù)衛(wèi)星過(guò)境時(shí)潮位和日最高潮位計(jì)算出潮位差，再由潮位差和海岸平均坡度計(jì)算外牽引距離(L)，計(jì)算公式為：

L=(H2-H1)/tanα，

(3)

式中：H1為瞬時(shí)潮位高度；H2為最高潮位高度；α為海岸平均坡度(°)。海岸線位置計(jì)算的原理如圖2所示，將需要校正的水邊線向前移動(dòng)距離L就可以得到真正意義的海岸線位置。

圖2 海岸線位置計(jì)算原理Fig.2 Theory of counting coastline position

2.2 海岸帶土地利用時(shí)空變化研究方法

海岸帶是指海陸之間相互作用的地帶，是海岸線向陸海兩側(cè)擴(kuò)展一定寬度的帶狀區(qū)域，包括陸域與近岸水域，對(duì)于其范圍，至今尚無(wú)統(tǒng)一的界定。本文規(guī)定的海岸帶調(diào)查范圍為：由海岸線向陸地方向延伸10 km，向海方向延伸2 km。利用各時(shí)相海岸線向陸區(qū)域做距離為10 km的緩沖區(qū)，向海區(qū)域做距離為2 km的緩沖區(qū)，以緩沖區(qū)為掩模裁剪各時(shí)相遙感影像。根據(jù)《土地利用現(xiàn)狀分類標(biāo)準(zhǔn)》的中的一級(jí)分類體系和研究區(qū)特征，將研究區(qū)土地利用類型分為水體、林地、耕地/草地和建設(shè)用地(主要為鹽田和養(yǎng)殖池)。在ENVI5.2軟件支持下，對(duì)各時(shí)相影像采用最大似然法與傳統(tǒng)目視解譯相結(jié)合的方法進(jìn)行監(jiān)督分類，以Google地圖與2018年9—10月實(shí)地考察數(shù)據(jù)為驗(yàn)證數(shù)據(jù)，對(duì)分類結(jié)果進(jìn)行精度驗(yàn)證[31-33]。

2.3 研究方法

在ArcMap10.2軟件支持下，基于研究區(qū)DEM數(shù)據(jù)，運(yùn)用水文分析工具，依次進(jìn)行洼地填充、流向分析、匯流累積量計(jì)算、提取河流網(wǎng)絡(luò)柵格，進(jìn)行流域分析，生成河網(wǎng)。把集水閾值設(shè)為10 000的河網(wǎng)圖層與研究區(qū)各時(shí)相海岸線圖層疊加，分析河口區(qū)海岸線變化特征。再結(jié)合各時(shí)相入海河流遙感影像解譯標(biāo)志，分析1988—2018年30 a間研究區(qū)入海河流河口區(qū)演化特征。

3 結(jié)果與分析

3.1 海岸線演變分析

3.1.1 海岸線潮位校正及精度驗(yàn)證

海水在天體引潮力作用下產(chǎn)生的潮汐現(xiàn)象具有規(guī)律性，日照海岸帶潮汐屬于半日潮型，即一日有2次高潮和低潮[34]。本文獲取了研究區(qū)2017—2019年6月的大潮日的潮高數(shù)據(jù)，并用2017—2019年各個(gè)時(shí)刻潮位數(shù)據(jù)擬合潮位曲線(圖3)。

圖3 2017—2019年潮位擬合曲線Fig.3 The fit curve of tidal level in 2017—2019

由圖3可知，不同年份潮位曲線具有相似的走勢(shì)，對(duì)高低潮時(shí)刻和潮高做相關(guān)性分析，相關(guān)系數(shù)均達(dá)到0.99以上。計(jì)算圖4中任意2條曲線同一時(shí)刻的潮位高度差，可以得到2017—2019年時(shí)刻的平均潮位高差為30.23 cm，根據(jù)三角函數(shù)關(guān)系計(jì)算此誤差在水平方向上為7.5 m，而衛(wèi)星分辨率為30 m，誤差在合理范圍內(nèi)，故可利用潮汐規(guī)律可以對(duì)缺失數(shù)據(jù)的年份做潮位修正。利用研究區(qū)DEM圖像生成坡度圖(圖4)，計(jì)算可得從水邊線向陸地推80 m的范圍內(nèi)平均坡度為2.3°，通過(guò)計(jì)算可得潮位校正外牽引距離(表2)。

圖4 研究區(qū)坡度圖Fig.4 Slope map of the study area

表2 潮位數(shù)據(jù)以及計(jì)算結(jié)果Tab.2 Tidal level data and calculation results

由表2可知，為得到研究區(qū)各時(shí)相正確的海岸線，需將獲得的1988年、1998年、2008年、2018年的未利用的淤泥質(zhì)海岸和砂質(zhì)海岸的瞬時(shí)水邊線分別向前移動(dòng)88 m，83 m，8 m和70 m，獲得各時(shí)相完整的海岸線。為了驗(yàn)證提取效果，首先在各時(shí)相原圖像不同類型海岸線上分別隨機(jī)取40個(gè)像素點(diǎn)，然后在提取出的海岸線圖像上找出相對(duì)應(yīng)的像素點(diǎn)，比較二者之間的位置變化。如果2個(gè)像素點(diǎn)的位置重疊則表示該點(diǎn)海岸線提取效果好；如果2個(gè)像素點(diǎn)相鄰則表示該點(diǎn)海岸線提取效果在誤差允許范圍之內(nèi)；否則，則表示該點(diǎn)海岸線提取失敗。經(jīng)計(jì)算可獲得不同時(shí)相海岸線像素點(diǎn)位置的檢驗(yàn)結(jié)果(表3)[35-36]。

表3 像素點(diǎn)位置檢驗(yàn)結(jié)果Tab.3 The veracity of pixel position

從檢驗(yàn)結(jié)果可以看出基巖海岸和人工海岸的提取效果最好，砂質(zhì)海岸提取效果較差，這是因?yàn)樯百|(zhì)海岸海岸線平直，易產(chǎn)生點(diǎn)的位移；而基巖海岸、砂質(zhì)海岸和未開(kāi)發(fā)淤泥質(zhì)海岸的海岸線相對(duì)曲折，海岸線位置變化不大。

3.1.2 海岸線動(dòng)態(tài)演變分析

把各時(shí)相海岸線疊加顯示(圖5)，通過(guò)ArcMap的數(shù)據(jù)統(tǒng)計(jì)功能提取出各類型海岸線的長(zhǎng)度(表4)，綜合分析不同類型海岸線時(shí)空變化情況。結(jié)合圖5和表4分析可知，1988—2019年30 a間日照海岸線變化劇烈，海岸線演變整體趨勢(shì)為向海推進(jìn)。海岸線長(zhǎng)度及增速均呈現(xiàn)逐年增加的趨勢(shì)，海岸線長(zhǎng)度總體增加達(dá)49.17 km，其中，1988—1998年增加12.1 km，年變化率為1.21 km/a；1998—2008年增加22.4 km，年變化率為2.24 km/a；2018年比2008年增加14.8 km，年變化率為1.48 km/a，3個(gè)時(shí)期均保持較快增長(zhǎng)，1998—2008年增速最快。1988—1998年間，除日照港附近的海岸線向海推進(jìn)外，其他地區(qū)海岸線均后退；1998—2008年，海岸線變遷顯著，隨著日照港到傅疃河口沿岸建設(shè)幅度增大，海岸線向海方向增加明顯；2008—2018年，韓家營(yíng)子到繡針河口段為嵐山港建設(shè)開(kāi)發(fā)區(qū)，此段海岸線變化劇烈。在不同類型海岸線中，人工海岸線所占比重呈逐年上升趨勢(shì)，30 a內(nèi)凈增加85.72 km，到2018年比重已達(dá)70.04%，海岸工程、養(yǎng)殖池是人工海岸線的主要形式，如中部的奎山和南部的嵐山兩段目前改造為碼頭、養(yǎng)殖池，南部的沙壩外普遍修筑了人工堤壩，沙壩間的潟湖區(qū)被改造為養(yǎng)殖池；砂質(zhì)海岸線分布較多，受海水侵蝕作用小，使沿海地區(qū)向陸面積變化小，但是1988年以后，砂質(zhì)海岸線長(zhǎng)度呈逐年遞減趨勢(shì)，30 a間減少12.84 km；基巖海岸線相間分布所占比例很少，主要分布在任家臺(tái)咀、石臼咀、嵐山咀和奎山咀，近幾年基巖海岸線變化很?。挥倌噘|(zhì)海岸線保持緩慢下降趨勢(shì)，主要分布在河流入海口處。

(a)1988—1998年 (b)1998—2008年

表4 各時(shí)期各類型海岸線長(zhǎng)度統(tǒng)計(jì)Tab.4 Length and proportion statistics of coastline in each period

3.2 海岸帶土地利用變化過(guò)程

3.2.1 海岸帶土地利用分類精度驗(yàn)證

在ENVI5.2軟件支持下，以Google地圖與2018年9—10月實(shí)地考察數(shù)據(jù)為驗(yàn)證數(shù)據(jù)，對(duì)進(jìn)行監(jiān)督分類后的分類結(jié)果進(jìn)行精度驗(yàn)證，以總體分類精度和Kappa系數(shù)來(lái)表示分類精度(表5)。由于本文分類系統(tǒng)簡(jiǎn)單，主要分為了4類，且各地類在光譜上有明顯區(qū)別，故而分類精度較高，從表5中可以看出其Kappa系數(shù)均在0.9以上，表明分類結(jié)果可以滿足研究目的。利用Post Classification工具消除噪聲，獲得研究區(qū)各時(shí)相海岸帶土地利用圖(圖6)[37-38]。

表5 分類精度Tab.5 Classification accuracy

3.2.2 海岸帶土地利用變化結(jié)果分析

得到研究區(qū)各時(shí)相土地利用分類結(jié)果后，利用分類后處理工具統(tǒng)計(jì)并計(jì)算各地類所占面積及其所占百分比(表6)，分析其結(jié)構(gòu)組成的時(shí)空變化。由表6可知，1988—2018年30 a間建設(shè)用地面積呈現(xiàn)逐年增加的趨勢(shì)，由1988年的213.77 km2增加到2018年的413.93 km2，其所占總地類面積的比重由21.73%增加到39.43%；30 a間耕地和草地面積則逐年減少，由1988年的445.50 km2減少到2018年的287.03 km2，其所占總地類面積的比重從45.29%減少到27.34%，嵐山區(qū)、東港區(qū)附近的大部分耕地轉(zhuǎn)換為建設(shè)用地；水體和林地分布時(shí)空變化起伏不大。

(a)1988年 (b)1998年

表6 1988—2018年日照海岸帶各土地利用類型面積及比例Tab.6 Land cover type area and proportion from 1988 to 2018 in Rizhao coastal zone

用2018年海岸帶緩沖區(qū)作為感興趣區(qū)裁剪1988年和2018年遙感影像，進(jìn)行監(jiān)督分類，制作1988—2018年轉(zhuǎn)移矩陣表(表7)，分析日照海岸帶30 a間各地類轉(zhuǎn)化情況。由表7可知，1988—2018年各地類相互轉(zhuǎn)化情況如下：有60.82 km2建設(shè)用地轉(zhuǎn)化為耕地/草地，占其轉(zhuǎn)化總面積的27.69%，相對(duì)應(yīng)的有185.51 km2耕地/草地轉(zhuǎn)化為建設(shè)用地，占其轉(zhuǎn)化總面積的比重最大，為39.60%，兩者間凈轉(zhuǎn)移量為124.69 km2；其次，水體轉(zhuǎn)化為建設(shè)用地的面積占其轉(zhuǎn)化為其他地類的總面積的比重最大，轉(zhuǎn)為建設(shè)用地的面積84.63 km2，占其轉(zhuǎn)化總面積的12.97%；林地轉(zhuǎn)為建設(shè)用地的面積也最大，為33.8 km2，占其轉(zhuǎn)化總面積的31.97%?？傮w來(lái)看，建設(shè)用地在各地類使用中的重要程度逐漸加大。

表7 1988—2018年日照海岸帶土地利用面積轉(zhuǎn)移矩陣Tab.7 Transition matrix of land over types in Rizhao coastal zone from 1988 to 2018

3.3 入海河流對(duì)海岸帶的影響

研究區(qū)沿岸的主要入海河流有兩城河、付疃河和繡針河，根據(jù)入海河流河口區(qū)特征(表8)和河口區(qū)海岸線位置變化(圖7)對(duì)日照海岸帶河口區(qū)演化進(jìn)行分析。結(jié)合表8和圖7分析可知，1988—1998年間，兩城河口段海岸線處于蝕退狀態(tài)，1988年后兩城河河口南側(cè)的沙壩潟湖水域出現(xiàn)了大面積圍海和向陸開(kāi)挖養(yǎng)殖池的活動(dòng)，導(dǎo)致潟湖消失；1998—2008年間，在兩城河下游和河道內(nèi)建設(shè)蝦池、育苗廠，大面積占用河道；2008—2018年間，對(duì)占用河道的設(shè)施進(jìn)行清除，并建立地下水庫(kù)攔蓄和調(diào)節(jié)地下水流，向海輸出的泥沙持續(xù)淤積，在河口區(qū)重新堆積砂體[39]。

表8 河口區(qū)特征Tab.8 The characteristics of the estuary

(a)1988—1998年兩城河河口區(qū) (b)1988—1998年付疃河河口區(qū) (c)1988—1998年繡針河河口區(qū)

(g)2008—2018年兩城河河口區(qū) (h)2008—2018年付疃河河口區(qū) (i)2008—2018年繡針河河口區(qū)

1988—1998年間，付疃河有一半泥沙被日照水庫(kù)攔淤，河流供沙大幅度減少，淤進(jìn)速度減緩，河口附近海岸線還被用于建設(shè)養(yǎng)殖池，轉(zhuǎn)變?yōu)槿斯ず０叮?998—2008年間，在付疃河流域境內(nèi)分布著多處采沙點(diǎn)，嚴(yán)重影響了付疃河的泄洪能力；2008年以來(lái)，對(duì)付疃河進(jìn)行了綜合治理，使其向海輸出的泥沙持續(xù)淤積，向海伸出形成突灘[40-41]。

1988—1998年間，嵐山港相繼建成防波堤和多個(gè)噸級(jí)泊位，攔阻了大量泥沙，使河口區(qū)蝕退；1998—2008年間，又在河口灣潟湖中修建了人工砌石護(hù)岸，岸灘泥沙被阻擋，人工開(kāi)采海灘沙和興修水庫(kù)也使入海泥沙大為減少，導(dǎo)致河口東岸的長(zhǎng)條狀沙嘴消失；2008—2018年間，由于政策保護(hù)河口濕地，河口三角洲處轉(zhuǎn)蝕為淤，河口處形成向陸微凹的小型三角洲沖積平原[42-43]。

4 討論

研究區(qū)海岸線變化與沿海土地利用之間是有直接聯(lián)系的，海岸線劇烈變化的區(qū)域也是海岸帶土地利用變化劇烈的區(qū)域。30 a間，研究區(qū)海岸線變化劇烈，北岸海岸線較為完整，變化不大；中岸和南岸海岸線變化較大，海岸線破碎。海岸線長(zhǎng)度持續(xù)增加，自然海岸線持續(xù)減少，人工海岸線持續(xù)增加，其中變化較大的區(qū)域主要集中在日照港、嵐山港及河流入海口附近，如日照港和嵐山港相繼建成防波堤和多個(gè)噸級(jí)港作泊位和生產(chǎn)泊位，作為建設(shè)開(kāi)發(fā)區(qū)；在河口區(qū)出現(xiàn)了大面積圍海和向陸開(kāi)挖養(yǎng)殖池的活動(dòng)、分布著多處采沙點(diǎn)。30 a間，研究區(qū)主要土地利用類型為建設(shè)用地、耕地/草地，建設(shè)用地面積呈現(xiàn)逐年增加的趨勢(shì)，養(yǎng)殖池和鹽田增加面積所占比重較大，區(qū)域內(nèi)嵐山區(qū)、東港區(qū)附近的大部分耕地轉(zhuǎn)換為建設(shè)用地，且建設(shè)用地?cái)U(kuò)張呈現(xiàn)中心擴(kuò)張趨勢(shì)；隨著入?？诤恿鲾y帶的泥沙逐漸減少及大面積圍海和向陸開(kāi)挖養(yǎng)殖池、采沙活動(dòng)的增加，沿岸水體減少，相應(yīng)建設(shè)用地增多；在海岸線變化較小的任家臺(tái)咀、石臼咀、嵐山咀、奎山咀地區(qū)，沿岸土地利用變化也相應(yīng)較小。由此可知，在海岸線變化劇烈的時(shí)間段內(nèi)沿岸土地利用速率也相應(yīng)變化。

5 結(jié)論

本文系統(tǒng)總結(jié)了日照海岸帶影像解譯標(biāo)志，并結(jié)合前人研究成果，利用RS和GIS技術(shù)，分析了日照海岸帶時(shí)空變化特征及其驅(qū)動(dòng)力，得出如下結(jié)論：

1)1988—2018年的30 a間日照海岸線變化劇烈，海岸線演變整體趨勢(shì)為明顯的向海推進(jìn)。海岸線長(zhǎng)度及增速均呈現(xiàn)逐年增加的趨勢(shì)，3個(gè)時(shí)期均保持較快增長(zhǎng)。在不同類型海岸線中，人工海岸線所占比重呈逐年上升趨勢(shì)；砂質(zhì)海岸線逐年減少；基巖海岸線相間分布于砂質(zhì)海岸之間，所占比例很少，幾乎不發(fā)生變化；淤泥質(zhì)海岸線始終變化不大，保持緩慢下降趨勢(shì)。

2)30 a間日照海岸帶各土地利用類型中，建設(shè)用地面積及其所占比重逐年遞增，且養(yǎng)殖池和鹽田所占比重較大；耕地和草地面積及其所占比重逐年減少；林地和水體面積時(shí)空變化不明顯。轉(zhuǎn)移矩陣表顯示各地類轉(zhuǎn)化以耕地/草地轉(zhuǎn)化為建設(shè)用地為主要趨勢(shì)。

3)30 a間河口區(qū)海岸線動(dòng)態(tài)變化明顯，1988—1998年間河口區(qū)進(jìn)入蝕退階段，沙嘴前端變淺，甚至部分解體，河口附近岸段基本上轉(zhuǎn)為人工海岸線；1998—2008年間各入海河流河口地貌保持相對(duì)穩(wěn)定，說(shuō)明河口區(qū)水動(dòng)力對(duì)沉積物的沖蝕和淤積基本保持動(dòng)態(tài)平衡；2008—2018年間河口區(qū)恢復(fù)以堆積狀態(tài)為主，在河流的兩側(cè)分布有規(guī)模較小的帶狀沖積平原。