圖們江口及鄰近區(qū)域海岸線遙感監(jiān)測(cè)及變化分析

苗俊偉，馬毅，張杰，孫偉富，胡亞斌，3

(1.山東科技大學(xué) 青島 266590;2.國家海洋局第一海洋研究所 青島 266061; 3.大連海事大學(xué) 大連 116026)

圖們江口及鄰近區(qū)域海岸線遙感監(jiān)測(cè)及變化分析

苗俊偉1，2，馬毅1，2，張杰1，2，孫偉富2，胡亞斌2，3

(1.山東科技大學(xué) 青島 266590;2.國家海洋局第一海洋研究所 青島 266061; 3.大連海事大學(xué) 大連 116026)

文章利用1985、1995、2005、2009年的Landsat衛(wèi)星遙感影像和2015年的GF-1衛(wèi)星遙感影像，分析近30年來圖們江口海岸線的變化特征。結(jié)果表明:近30年來圖們江口海岸線長度整體呈增長態(tài)勢(shì)，累積增加0.66 km，主要由人類開發(fā)活動(dòng)及泥沙堆積等引起;人工岸線比例增大，但1985—1995年因河流改道致使岸線減少3.8 km。1985—2005年海岸帶淤積面積大于侵蝕面積，其中河流改道使1995—2005年淤積面積達(dá)到4.522 7 km2;2005—2015年侵蝕面積大于淤積面積; 2005—2009年徑流量減少使侵蝕面積達(dá)到4.203 7 km2。河口處的岸線遷移變化最大，其他地區(qū)較穩(wěn)定。

遙感;圖們江口;海岸線;岸線變遷

1 引言

中華人民共和國國家標(biāo)準(zhǔn)《海洋學(xué)術(shù)語:海洋地質(zhì)學(xué)》(GB/T18190—2000)對(duì)海岸線的定義是:海陸分界線，在我國系指多年大潮平均高潮位時(shí)海陸分界線［1］。海岸線不僅標(biāo)志沿海地區(qū)的海陸分界線，而且包含豐富的環(huán)境信息，其變化可以直接影響潮間帶的灘涂資源量以及海岸帶的環(huán)境變化，引起海岸帶資源與生態(tài)過程的改變，影響人類生產(chǎn)活動(dòng)。因此，對(duì)海岸線進(jìn)行動(dòng)態(tài)變化監(jiān)測(cè)，有利于了解海岸帶的生態(tài)環(huán)境變化。

傳統(tǒng)的海岸線獲取方法主要依靠常規(guī)測(cè)量手段，近年來Baptista等［2］開始利用GPS結(jié)合沙灘摩托和步行方法進(jìn)行岸線測(cè)量，雖然岸線提取精度很高，但現(xiàn)場(chǎng)踏勘費(fèi)時(shí)費(fèi)力且很多地區(qū)人員無法到達(dá)。與之相比，遙感技術(shù)具有大范圍、時(shí)效性強(qiáng)、節(jié)約資源等優(yōu)勢(shì)，對(duì)于獲取海岸線長時(shí)間的變化具有獨(dú)特優(yōu)勢(shì)。

國內(nèi)外在海岸線變遷遙感分析上成果頗多。在海岸線自動(dòng)提取方面，Chen等［3］利用直方圖均衡化及Canny邊緣檢測(cè)算子提取臺(tái)灣西部外傘頂洲的水邊線;王李娟［4］、晁暉等［5］用Sobel算法和修復(fù)歸一化水體指數(shù)提取海岸線;以馮永玖、劉云朋為代表的學(xué)者利用CV模型［6］、元胞自動(dòng)機(jī)［7－8］和岸線追蹤技術(shù)［9］等方法提取海岸線，并在人工岸線、砂質(zhì)岸線及淤泥質(zhì)岸線進(jìn)行應(yīng)用和分析;Kuleli等［10］利用Landsat TM和ETM+影像的水體指數(shù)進(jìn)行閾值處理，對(duì)二值圖像矢量化后提取土耳其5個(gè)海岸帶濕地的水邊線。

在海岸線人－機(jī)交互解譯提取方面，部分學(xué)者界定遙感圖像中海岸線的確切位置，給出相應(yīng)提取原則，并對(duì)岸線提取結(jié)果進(jìn)行精度驗(yàn)證［11－15］。Misra等［16］基于Landsat衛(wèi)星數(shù)據(jù)，利用不同波段之間的比值區(qū)分海陸分界線，進(jìn)而人工提取海岸線;Kevin White等［17］基于Landsat系列衛(wèi)星遙感影像提取尼羅河三角洲河口區(qū)域的岸線，并研究岸線動(dòng)態(tài)特征;Mujabar［18］通過人－機(jī)交互方法提取印度兩地之間的岸線，并利用 DSAS分析海岸淤蝕狀況;張怡［19］、包萌［20］、趙宗澤［21］、周相君［22］基于衛(wèi)星影像，采用人－機(jī)交互方法提取海岸線，并對(duì)其變化進(jìn)行時(shí)空分析。

遙感具有范圍廣、多時(shí)相等特點(diǎn)，國內(nèi)外采用遙感手段進(jìn)行海岸線變遷監(jiān)測(cè)與分析較廣泛。本文在建立的各岸線類型解譯標(biāo)志的基礎(chǔ)上，利用1985、1995、2005、2009年的 Landsat衛(wèi)星影像和2015年GF-1衛(wèi)星遙感影像對(duì)圖們江口及鄰近區(qū)域海岸線進(jìn)行提取，從岸線長度變遷、岸線類型變化和海岸帶侵蝕淤積狀況等方面對(duì)圖們江口地區(qū)海岸線變遷情況進(jìn)行分析，監(jiān)測(cè)該區(qū)海岸線近30年的變化，為該區(qū)岸線保護(hù)與開發(fā)提供科學(xué)參考。

2 研究區(qū)域

圖們江大部分在中國境內(nèi)，以琿春市為截點(diǎn)出境，注入日本海。圖們江口地處130°22'E—131°02' E、42°14'N—42°44'N，涉及中、朝、俄三國，毗鄰韓、日兩國，是東北亞地區(qū)重要樞紐，具有優(yōu)越的地理位置與重要的戰(zhàn)略地位，未來該地區(qū)將成為我國在東北地區(qū)與其他國家進(jìn)行合作交流的橋頭堡，推動(dòng)?xùn)|北地區(qū)大發(fā)展。

該區(qū)屬溫帶大陸性氣候，臨近日本海，在東南季風(fēng)、地形因素和日本海域的影響下，降水和徑流都比較豐富。圖們江口海岸線類型以自然岸線為主，近年來由于開發(fā)活動(dòng)加劇，河口處岸線類型及長度發(fā)生較大變化。由于地理位置特殊，對(duì)該區(qū)域海岸線提取及變遷的研究尚為空白，本文基于衛(wèi)星遙感影像進(jìn)行研究。

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

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

為監(jiān)測(cè)圖們江口地區(qū)海岸線近30年來的變化，本文應(yīng)用 1985年、1995年、2005年、2009年的Landsat影像和2015年的GF-1影像(表1)。

表1 研究區(qū)遙感數(shù)據(jù)源信息

Landsat5 TM 30 1995年11月 2 Landsat7 ETM+ 30 2005年10、11月 2 Landsat5 TM 30 2009年10月 2 GF-1 PMS 8 2015年3月2

利用拼接、裁剪后的2005年ETM+衛(wèi)星影像數(shù)據(jù)對(duì)TM及經(jīng)過正射校正的GF-1衛(wèi)星數(shù)據(jù)進(jìn)行幾何配準(zhǔn)，配準(zhǔn)誤差控制在1個(gè)像元之內(nèi)，并對(duì)配準(zhǔn)后的圖像進(jìn)行拼接裁剪，獲得研究區(qū)衛(wèi)星影像，利用近紅外、紅、綠波段組合開展海岸線提取。

3.2 方法

采用目視解譯方法，根據(jù)不同海岸線類型在衛(wèi)星影像中的特征建立海岸線類型的解譯標(biāo)志，進(jìn)而提取圖們江口海岸線。根據(jù)當(dāng)?shù)貙?shí)際狀況與衛(wèi)星遙感影像信息，結(jié)合孫偉富等［13］建立的海岸線類型遙感解譯標(biāo)志，將該地區(qū)海岸線分為基巖岸線、砂質(zhì)岸線、粉砂淤泥質(zhì)岸線與人工岸線4種，參考《我國近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)海岸線修測(cè)技術(shù)規(guī)程》［23］中的界定原則，結(jié)合各類海岸線的地貌特征，確定各類型岸線的位置。

基巖岸線位于由巖石組成的基巖海岸上，常有突出的海岬和深入陸地的海灣，岸線曲折。在影像中基巖岸線位于明顯的水陸分界線上(圖1(a))。

砂質(zhì)岸線位于砂質(zhì)海岸上，砂質(zhì)海岸通常較為平直，在大潮潮水搬運(yùn)下常堆積成一條與海岸平行的脊?fàn)畹纳百|(zhì)沉積——灘脊，該處即為海岸線位置。海岸的干燥灘面光譜反射率較高，在影像上表現(xiàn)為白亮的區(qū)域，灘脊痕跡線處堆積有植物碎屑、雜物等，亮度較低，海水的光譜反射率較低，含水量較高的沙灘光譜反射率也較低，在影像上表現(xiàn)略暗。在影像中砂質(zhì)岸線的位置應(yīng)取在灘脊痕跡線上限處(圖1(b))。

粉砂淤泥質(zhì)岸線位于淤泥質(zhì)海岸上，該類型海岸多由潮汐作用形成，受上沖流影響灘面坡度平緩，寬度可達(dá)數(shù)千米甚至更大。粉砂淤泥質(zhì)海岸向陸一側(cè)植被生長茂盛，成紅色或者暗紅色;向海一側(cè)植被較稀疏，呈現(xiàn)為淺紅色或沒有植被，裸露潮灘上多有樹枝狀潮溝發(fā)育，潮灘淹沒時(shí)高潮線處植被稀疏。在影像中植被茂盛與稀疏程度明顯差異處即為粉砂淤泥質(zhì)岸線所在位置(圖1(c))。

人工岸線是人工建筑物形成的岸線，在影像上建筑物一般亮度較高，呈白色，狹長延伸分布，紋理平滑，堤外淤泥質(zhì)灘涂顏色灰暗，人工岸線一般確定在建筑物外緣(圖1(d))。

圖1 海岸線遙感解譯標(biāo)志

在海岸線提取結(jié)果上，從岸線長度、岸線類型、海岸帶侵蝕淤積、遷移速率4個(gè)方面對(duì)5期海岸線進(jìn)行變化監(jiān)測(cè)，并根據(jù)衛(wèi)星圖像及文獻(xiàn)資料分析岸線變化的主要原因。

4 結(jié)果與分析

4.1 岸線時(shí)空變遷總體分析

通過5期衛(wèi)星遙感圖像的海岸線提取結(jié)果可以發(fā)現(xiàn):1985—2015年圖們江口區(qū)域的海岸線變遷主要發(fā)生在河口地區(qū)，其他地區(qū)變化較弱;目前岸線長度為320.78 km，較1985年增加0.66 km，陸域面積增加1.000 4 km2，泥沙堆積使岸線長度增加，港口碼頭建設(shè)也增加岸線長度。

海岸線變遷程度受自身與外部因素影響。一方面，在海岸線各類型中，基巖岸線最為穩(wěn)定，受外界自然因素影響較小，而砂質(zhì)岸線與粉砂淤泥質(zhì)岸線較為脆弱，易受環(huán)境影響;另一方面，自然因素對(duì)海岸線的變遷作用較為緩慢，而人為因素可在極短時(shí)間內(nèi)產(chǎn)生自然因素影響幾十年才能產(chǎn)生的變化，人類活動(dòng)一直在岸線變遷中起主導(dǎo)作用。

4.2 岸線類型轉(zhuǎn)化分析

根據(jù)提取的5期圖們江口區(qū)域海岸線，獲取5期岸線類型信息(表2)。

表2 不同年份各類型海岸線的長度統(tǒng)計(jì)

結(jié)果表明:岸線總長度及各類型岸線長度均存在波動(dòng)變化。1985—2015年由于人類活動(dòng)使基巖岸線減少;砂質(zhì)岸線整體增長，主要是河道內(nèi)泥沙堆積引起，其中1985—1995年砂質(zhì)岸線大量減少的主要原因是河流改道，1995—2005年河流再次改道，砂質(zhì)岸線有所增加;由于修建港口碼頭、防護(hù)堤等人類活動(dòng)，人工岸線持續(xù)增加;粉砂淤泥質(zhì)岸線的變化與河流徑流量及海浪沖刷關(guān)系較大，但1985—2005年粉沙淤泥質(zhì)岸線減少的原因與砂質(zhì)岸線類似，均為河流改道，而2005—2009年河道再次改變，粉砂淤泥質(zhì)岸線長度增加，但與改道前相比長度減少，2009年之后海浪沖刷使河口處粉砂淤泥質(zhì)岸線繼續(xù)減少。整體上自然岸線比例減小、人工岸線比例增大，說明該區(qū)開發(fā)利用進(jìn)一步加強(qiáng)。

4.3 海岸帶侵蝕淤積分析

圖們江口及鄰近區(qū)域1985—2015年由于海岸線變遷導(dǎo)致陸地面積凈增1.000 4 km2，平均每年增加0.033 3 km2(表3)。

表3 海岸帶侵蝕淤積統(tǒng)計(jì) km2

圖們江口區(qū)域1985—2015年4個(gè)時(shí)段的變遷情況不同，主要是由海水侵蝕、河流徑流量變化、人類開發(fā)與泥沙堆積共同作用引起的。1985—1995年雖然河流改道侵蝕面積較大，但河流徑流量充沛，泥沙堆積使岸線加長，加之圍海造陸使淤積面積大于侵蝕面積，海岸帶總面積增加;1995—2005年河流再次改道，使淤積面積遠(yuǎn)大于上一時(shí)期，海岸帶增長面積也比上一時(shí)期大;2005—2009年降水減少、徑流量減小，泥沙補(bǔ)充不足，原本在河口處的部分陸地變成沙洲，導(dǎo)致海岸帶侵蝕加大，該時(shí)期內(nèi)海岸帶面積減小;2009—2015年徑流量增大，但洪水沖刷河道導(dǎo)致河道侵蝕，使海岸帶面積減小。

圖們江口區(qū)域1985—2015年海岸帶地區(qū)侵蝕淤積變化共存。其中岸線減少的主要原因是河流改道、泥沙補(bǔ)充不足導(dǎo)致侵蝕，以及在入海口地區(qū)海浪對(duì)砂質(zhì)岸線的沖刷侵蝕。如，由于圖們江位于松遼流域，根據(jù)松遼流域河流泥沙公報(bào)，2005—2009年該流域多屬于枯水少沙狀態(tài)，水量、沙量明顯減少［24］，泥沙補(bǔ)充不足;同時(shí)在日本海海浪的侵蝕下，河口堆積的泥沙后退，侵蝕面積增大。2009—2015年該流域河流徑流量變大、總輸沙量增大，尤其是2010年發(fā)生洪水災(zāi)害［24］，使河流徑流增大、攜帶泥沙量增大，導(dǎo)致2009—2015年砂質(zhì)岸線略微變長，海岸帶淤積面積較上一時(shí)期增大;但由于發(fā)生洪水時(shí)河流水速較快，對(duì)河道猛烈沖刷產(chǎn)生侵蝕，加之海浪的侵蝕作用，使侵蝕面積大于淤積面積，因此海岸帶面積依舊呈減少態(tài)勢(shì)。

4.4 海岸線遷移分析

目前量化分析海岸線變遷的主要方法較多［25］，本研究選用基線法分析圖們江口地區(qū)30年來的海岸線遷移狀況。

基線法最早由 Doln等［26］提出，經(jīng)過 Thieler等［27］的進(jìn)一步完善，形成基于ArcGIS平臺(tái)的“Digital Shoreline Analysis System，DSAS”功能模塊，用于分析海岸的時(shí)空變化狀況。采用DSAS分析模塊，由海岸線向陸縱深得到基線并使海岸線在基線同側(cè)，然后從基線向海岸線一側(cè)做長度為D的垂線，與基線交點(diǎn)間的距離為基線采樣間隔，每根垂線與各時(shí)期的海岸線相交，將每個(gè)交點(diǎn)至基線的距離相減，可得海岸線的變化距離。

本文選用的基線由手工繪制，基線采樣間隔為50 m，由DSAS分析得到最終時(shí)空變化(圖2)。

圖2 各時(shí)期不同區(qū)域終點(diǎn)變化速率分布*注:加重線段為河口區(qū)域

從海岸線終點(diǎn)變化速率及其分布空間來看，1985—2015年河口地區(qū)(斷面號(hào)大約位于669～1 088)海岸線變化較為明顯，其他地區(qū)變化較少。1985—1995年由于河道部分地區(qū)泥沙淤積使岸線向水域推進(jìn)，導(dǎo)致終點(diǎn)變化速率為正值，而河口地區(qū)由于河流改道使河口邊界向陸地推移，終點(diǎn)變化速率為負(fù)值且變化速率較大。1995—2005年河口位置再次發(fā)生變化，向水域猛烈推進(jìn)，導(dǎo)致終點(diǎn)變化速率值巨大且為正值。2005—2009年河口推進(jìn)情況與上一時(shí)段相反，終點(diǎn)變化速率值為負(fù)值且數(shù)值巨大。2009—2015年由于部分地區(qū)河道位置改變以及河口岸線位置改變，加之岸線推移較短，該時(shí)段終點(diǎn)變化速率較小、范圍較廣。除河口區(qū)域外的其他變化區(qū)域終點(diǎn)變化速率值較小且變化范圍較小，主要是因?yàn)槿祟惢顒?dòng)少，河流寬度較窄、多為細(xì)小河流，河道位置變化較小，河流徑流量對(duì)河口區(qū)域影響較弱。

通過對(duì)比不同時(shí)期的變化可以發(fā)現(xiàn)，發(fā)生變化的斷面越來越多，橫軸右側(cè)處(約在斷面號(hào)2 200后)波動(dòng)區(qū)域增加;根據(jù)斷面號(hào)定位到遙感圖像上，發(fā)現(xiàn)變化地區(qū)向北推進(jìn)且變化率較大的斷面號(hào)多分布在砂質(zhì)岸線與基巖岸線處，由此說明該區(qū)砂質(zhì)岸線受海浪影響加重，查閱資料并將斷面號(hào)定位到遙感影像上可知，該地區(qū)沿海開發(fā)活動(dòng)進(jìn)一步增強(qiáng)。

5 岸線變遷原因

海岸線變遷驅(qū)動(dòng)力復(fù)雜，通常受自然和人為因素雙重影響，前者包括河流輸沙量、風(fēng)暴潮、氣候變化、海浪和潮汐等，后者包括海岸工程采砂、岸堤修建、養(yǎng)殖池和鹽池建設(shè)、圍海造陸等。通過閱讀研究區(qū)的歷史資料及相關(guān)文獻(xiàn)，并參考不同時(shí)期間的遙感衛(wèi)星影像，結(jié)合本文得出的岸線分析結(jié)果，進(jìn)一步分析岸線變遷原因。

5.1 自然因素

從長時(shí)間序列看，海岸線變遷的主要自然因素是地殼運(yùn)動(dòng)和海平面升降;從短時(shí)間序列看，海岸線變遷的主要自然因素為海岸侵蝕、海水入侵、海浪、潮汐和風(fēng)暴潮等。具體到圖們江口，影響其大陸海岸線變遷的自然原因主要有以下幾個(gè)方面。

5.1.1 波浪及潮流的作用

改變海岸地貌的主要外動(dòng)力為波浪和潮流。波浪在海岸塑造中最積極、最活躍，靠近東海的朝鮮海峽(大韓海峽)處的海潮潮差可達(dá)2 m，劇烈的海浪對(duì)該區(qū)岸線的變化起到推動(dòng)作用;潮流是泥沙運(yùn)營的主要?jiǎng)恿Γ貣|岸北流的對(duì)馬暖流和沿西岸南流的利曼寒流組成日本海海流，形成反時(shí)針型的環(huán)流系統(tǒng)［28］。在波浪、潮流等動(dòng)力作用下，近岸泥沙發(fā)生轉(zhuǎn)運(yùn)，以沙嘴的形式形成堆積岸段;隨著沙嘴的遷移，岸段另一端發(fā)生侵蝕，形成相對(duì)侵蝕的海岸。以2005—2009年為例，在圖們江入海口附近，由于形成反時(shí)針型環(huán)流系統(tǒng)，在近海一端波浪從北向南流，在海浪和潮汐作用下海岸發(fā)生侵蝕;而在遠(yuǎn)海一端海浪侵蝕不到，從而發(fā)生淤積。

5.1.2 河流徑流的作用

該區(qū)位于歐亞大陸東岸，屬溫帶大陸氣候，受東南季風(fēng)、地形因素和日本海域的影響，降水和徑流都比較豐富，因此降水對(duì)徑流影響較大，二者分布規(guī)律大體相一致。近年來氣候變化等因素使河流輸沙逐漸減少，通過查閱松遼流域河流公報(bào)［24］可以發(fā)現(xiàn)，2005—2009年由于流域內(nèi)徑流量減少、河流輸沙量減少，導(dǎo)致河流入海口附近的泥沙堆積減少;而2010年發(fā)生洪水災(zāi)害，河流帶來大量泥沙，并在河道及河口堆積，使該流域入海口的岸線又變長。

5.2 人為因素

該區(qū)海岸線變化的主要原因是以填海造地為代表的各種人為活動(dòng)。在人類活動(dòng)的影響下，自然岸線總長度縮短、人工岸線增加、岸線形態(tài)更為曲折。總體看來，以下因素對(duì)于圖們江口地區(qū)海岸線時(shí)空動(dòng)態(tài)變化特征具有不同程度的影響。

5.2.1 港口碼頭的修建

該區(qū)地理位置復(fù)雜重要，各國為提高在該區(qū)的實(shí)力在海灣地區(qū)大力修建港口，促進(jìn)海上經(jīng)濟(jì)發(fā)展。大規(guī)模的建設(shè)改變海灣的位置形狀，如在朝鮮造山灣附近修建港口使岸線位置及類型發(fā)生變化。

5.2.2 防護(hù)堤壩的修建

近年來海平面上升和極端氣象災(zāi)害頻發(fā)，使岸線遭受嚴(yán)重威脅。為加固海岸，人工堤壩逐漸出現(xiàn)，保護(hù)岸線不受侵蝕，如河口南部地區(qū)在原始自然岸線的基礎(chǔ)上修建人工岸線及堤壩。防護(hù)堤壩在人工岸線中所占比例較大，也是該區(qū)人工岸線變化的主要類型。

6 結(jié)論

本文利用圖們江口地區(qū)1985年、1995年、2005年、2009年的Landsat影像和2015年的GF-1影像，基于GIS、RS技術(shù)，對(duì)5期海岸線進(jìn)行信息提取，分析圖們江口地區(qū)海岸線長度、類型、侵蝕淤積、岸線遷移的變化情況。

將研究區(qū)海岸線分為基巖岸線、砂質(zhì)岸線、粉砂淤泥質(zhì)岸線和人工岸線4類，以近紅外、紅、綠波段組合假彩色影像，根據(jù)不同岸線的特征建立圖們江口地區(qū)主要岸線類型的遙感解譯標(biāo)志，并給出4類岸線的提取原則，經(jīng)過目視解譯提取得到最終的海岸線。

從岸線的長度、類型、侵蝕淤積、岸線遷移等方面對(duì)5期岸線的分布狀況分別進(jìn)行統(tǒng)計(jì)分析。結(jié)果表明，30年間圖們江口地區(qū)岸線長度增減共存，但總體上呈增加趨勢(shì);1985—2015年共增加0.66 km，原因是泥沙堆積使砂質(zhì)岸線增長，而且人類開發(fā)活動(dòng)也使岸線有所增長。從整個(gè)研究時(shí)期看，30年間岸線類型的轉(zhuǎn)化較為平緩，其中減少最多的是基巖岸線，共減少0.87 km，減少的岸線多轉(zhuǎn)變?yōu)槿斯ぐ毒€;增長最多的則是人工岸線與砂質(zhì)岸線，均為1.14 km，岸線變遷態(tài)勢(shì)較穩(wěn)定。1985—2005年淤積面積大于侵蝕面積，其中1985—1995年因徑流量充沛和人類開發(fā)活動(dòng)加大，淤積面積為2.238 0 km2、侵蝕面積僅為0.892 9 km2;1995—2005年河流改道使淤積面積達(dá)到4.522 7 km2，侵蝕面積較上一時(shí)間段約減少50%;2005—2015年侵蝕面積大于淤積面積，2005—2009年因徑流量減少，侵蝕面積達(dá)到4.203 7 km2;2009—2015年侵蝕面積又有大幅度減少。30年間因海岸線變遷導(dǎo)致陸域面積累積新增1.000 4 km2。通過對(duì)比4個(gè)時(shí)期的終點(diǎn)變化速率分布可以發(fā)現(xiàn)，整體上岸線波動(dòng)范圍較小，河口地區(qū)岸線變化波動(dòng)較大且變化區(qū)域向北移動(dòng)，變化點(diǎn)多分布于砂質(zhì)岸線與基巖岸線處。

人類活動(dòng)在岸線變遷中作用時(shí)間短、效果明顯，是岸線變遷的主要推動(dòng)力，如港口及防護(hù)堤壩的修建都能在短時(shí)間內(nèi)改變岸線長度及岸線類型。而自然因素在岸線變遷中推動(dòng)緩慢、效果不明顯，如在不發(fā)生特殊氣象災(zāi)害的情況下，海浪及河流徑流量對(duì)岸線的影響效果要在幾年內(nèi)才可顯現(xiàn)，且作用效果較小。

由于本文使用的Landsat影像為中低分辨率影像，對(duì)海岸線的目視解譯造成一定的困難，使海岸線的提取結(jié)果存在偏差，今后將選用分辨率更高的影像進(jìn)一步研究。此外，在岸線精度評(píng)價(jià)中，由于該區(qū)地理位置特殊，研究尚不全面，如缺乏實(shí)測(cè)岸線，不能更好地判斷本文提取岸線的精度。今后將對(duì)該地區(qū)收集實(shí)測(cè)數(shù)據(jù)，獲得實(shí)際海岸線位置。

［1］GB/T18190—2000海洋學(xué)術(shù)語:海洋地質(zhì)學(xué)［S］．北京:中國標(biāo)準(zhǔn)出版社，2000．

［2］BAPTISTA P，CUNHA T，BERNARDES C，et al．A precise and efficient methodology to analyse the shoreline displacement rate［J］．Journal of Coastal Research，2005，27(2):223－232．

［3］CHEN W W，CHANG H K．Estimation of shoreline position and change from satellite images considering tidal variation［J］．Estuarine Coastal＆Shelf Science，2009，84(1):54－60．

［4］王李娟，牛錚，趙德剛，等．基于ETM遙感影像的海岸線提取與驗(yàn)證研究［J］．遙感技術(shù)與應(yīng)用，2010，25(2):235－239．

［5］晁暉，梁娜．遙感技術(shù)支持下海岸線提取算法［J］．河北聯(lián)合大學(xué)學(xué)報(bào):自然科學(xué)版，2014(3):82－85．

［6］劉云朋，王建梅，樓立志．基于ALOS影像與CV模型的瞬時(shí)海岸線提取［J］．測(cè)繪與空間地理信息，2014(2):91－95．

［7］馮永玖，韓震．海岸線遙感信息提取的元胞自動(dòng)機(jī)方法及其應(yīng)用［J］．中國圖像圖形學(xué)報(bào)，2012，17(3):441－446．

［8］FENG Y J，LIU Y，LIU D．Shoreline mapping with cellular automata and shoreline progradation analysis in Shanghai，China from 1979 to 2008［J］．Arabian Journal of Geosciences，2014，DOI:10．1007/s12517－0140－1515－7．

［9］馮永玖，袁佳宇，宋麗君，等．杭州灣海岸線信息的遙感提取及其變遷分析［J］．遙感技術(shù)與應(yīng)用，2015，30(2):345－352．

［10］KULELI T，GUNEROGLU A，KARSLI F，et al．Automatic detection of shoreline change on coastal Ramsar wetlands of Turkey［J］．Ocean Engineering，2011，38(10):1141－1149．

［11］張霞，莊智，張旭凱，等．秦皇島市海岸線遙感提取及變化監(jiān)測(cè)［J］．遙感技術(shù)與應(yīng)用，2014，29(4):625－630．

［12］巢子豪，高一博，李彥輝，等．基于遙感和GIS的海州灣海岸線提取方法研究［J］．淮海工學(xué)院學(xué)報(bào):自然科學(xué)版，2014(4):87－91．

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

［14］WANG Changying，ZHANG Jie，SONG Pingjian．An intelligent coastline interpretation of several types of seacoasts from TM/ ETM+images based on rules［J］．Acta Oceanological Sinica-English Edition，2014，33(7):89－96．

［15］MARFAI M A，ALMOHAMMAD H，DEY S，et al．Coastal dynamic and shoreline mapping:multi-sources spatial data analysis in Semarang Indonesia．［J］．Environmental Monitoring＆ Assessment，2008，142(1－3):297－308．

［16］MISRA A，BALAJI R．A Study on the Shoreline Changes and LAND-use/Land-cover along the South Gujarat Coastline［J］．Procedia Engineering，2015，116(1):381－389．

［17］KEVIN W，HESHAM M．Monitoring changing position of coastlines using thematic mapper imagery，an example from the nile delta［J］．Geomorphology，1999，29:93－105．

［18］MUJABAR S，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．

［19］張怡，李曉敏，馬毅，等．基于遙感的珠江口海岸線變遷分析［J］．海洋測(cè)繪，2014，34(3):52－55．

［20］包萌，孫偉富，馬毅，等．近40年來清瀾灣海岸線及其鄰接地物遙感監(jiān)測(cè)與變遷分析［J］．海岸工程，2014，33(02):66－76．

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

［22］周相君，李曉敏，馬毅，等．基于遙感的廣西防城灣海岸線變遷分析［J］．海洋學(xué)研究，2014，32(1):47－55．

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

［24］水利部松遼水利委員會(huì)．松遼流域河流泥沙公報(bào)2006－2012［Z］．2006－2012．

［25］姚曉靜，高義，杜云艷，等．基于遙感技術(shù)的近30年海南島海岸線時(shí)空變化［J］．自然資源學(xué)報(bào)，2013(1):114－125．

［26］DOLAN R，HAYDEN B，HEYWOOD J．A new photogrammetric method for determining shoreline erosion［J］．Coastal Engineering，1978，78(2):21－39．

［27］THIELER E R，HIMMELSTOSS E A，ZICHICHI J L，et al．The digital shoreline analysis system(DSAS)version 4．0-an ArcGIS extension for calculating shoreline change［J］．USGS Open-File Report，2008－1278 2009．

［28］佚名．日本海區(qū)的自然環(huán)境條件［J］．海洋通報(bào)，1977(6): 28－43．

Coastline Monitoring by Remote Sensing and Change Analysis in Tumen River Estuary and the Adjacent Area

MIAO Junwei1，2，MA Yi1，2，ZHANG Jie1，2，SUN Weifu2，HU Yabin2，3

(1.Shandong University of Science And Technology，Qingdao 266590，China;2.The First Institute of Oceanography，SOA，Qingdao 266061，China;3.Dalian Maritime University，Dalian 116026，China)

Four years'Landsat satellite images of 1985，1995，2005，2009 and one year's GF-1 remote sensing satellite images of 2015 were used to analyze the changing characteristics of the past 30 years in the Tumen River estuary coastline．The results showed that:over the past 30 years，the whole coastline length shows an increasing trend，mainly caused by human development activities and sediment buildup，and the cumulative increase is 0．66 km．The ratio of artificial shoreline increased，but between 1985 and 1995，the coastline length decreased 3．8 km because the river course was diverted．Between 1985 and 2005，the deposition area of the coastal zone was larger than that of the erosion area，and the river course was diverted to make the deposition area increase 4.522 7 km2between 1995 and 2005．From 2005 to 2015，the erosion area was greater than the deposition area．From 2005 to 2009，the erosion area reachedto 4.203 7 km2because the runoff had been reduced．The shoreline at the mouth had the greatest changes in migration and other regions were more stable．

Remote sensing，Tumen River estuary，Coastline，Coastline change

P74

A

1005－9857(2017)02－0039－07

2016-10-15;

2016-12-17

高分海岸帶遙感監(jiān)測(cè)與應(yīng)用示范．

苗俊偉，碩士研究生，研究方向?yàn)楹Ｑ筮b感與應(yīng)用，電子信箱:rsmiaojunwei@163.com

張杰，研究員，博士，研究方向?yàn)楹Ｑ筮b感，電子信箱:zhangjie@fio.org.cn