裘誠,王岳峰,方略任,張何朕加,張凱嬪,胡進(jìn),陳沈良,袁瑞

(1.上海市海洋監(jiān)測(cè)預(yù)報(bào)中心 上海 200062;2.上海海事大學(xué)海洋科學(xué)與工程學(xué)院 上海 201306;3.華東師范大學(xué)河口海岸學(xué)國家重點(diǎn)實(shí)驗(yàn)室 上海 200241)

0 引言

潮灘是廣泛發(fā)育的典型的海岸地貌單元,其演變特征和規(guī)律是海洋科學(xué)和海岸工程領(lǐng)域重要的研究方向和實(shí)際問題。在人類活動(dòng)密集的海岸區(qū)域,沉積物輸移和鹽沼植被生長的耦合作用對(duì)于潮灘地貌演變有著顯著的影響[1-4]。在河口海岸區(qū)域,潮灘演變是多種時(shí)間尺度因素共同作用下的結(jié)果,如海平面上升(長時(shí)間尺度)、洪枯季轉(zhuǎn)換(季節(jié)尺度)和極端天氣過程(事件尺度)等。因此,獲取更高時(shí)空分辨率的潮灘資料是研究和認(rèn)識(shí)其演變規(guī)律的重要手段。

傳統(tǒng)的海岸帶監(jiān)測(cè)手段一般采用人工現(xiàn)場測(cè)量(GPS-RTK等)和遙感反演[5],人工碎步測(cè)量方法具有可控的精度,但是時(shí)間長、效率低,且涉及人員作業(yè)安全問題,很難在有效的時(shí)間窗口內(nèi)開展大范圍的潮灘地形測(cè)量;衛(wèi)星遙感圖像分辨率低、時(shí)間周期長、波段冗余信息多,對(duì)于特定的興趣點(diǎn)與時(shí)間窗口無法有效地獲得信息,而潮灘地貌變化迅速,潮周期內(nèi)會(huì)發(fā)生較大變化,遙感技術(shù)多用于單次監(jiān)測(cè),難以用于揭示潮灘系統(tǒng)季節(jié)性和短期極端事件的動(dòng)態(tài)變化。

相比傳統(tǒng)監(jiān)測(cè)手段,無人機(jī)作為一種有動(dòng)力、可控制、能攜帶多種任務(wù)設(shè)備、執(zhí)行多種任務(wù)的低空飛行平臺(tái),其與遙感技術(shù)相結(jié)合具有成本低、操作簡單、影像獲取速度快和高空間分辨率等傳統(tǒng)衛(wèi)星遙感技術(shù)所無可比擬的優(yōu)勢(shì),已被廣泛用于國土監(jiān)測(cè)與城市管理、地質(zhì)災(zāi)害、環(huán)境監(jiān)測(cè)以及應(yīng)急保障等重要領(lǐng)域[6],如海洋環(huán)境監(jiān)測(cè)[7]、海岸線監(jiān)測(cè)[8]、海島測(cè)繪[9]等。上述基于無人機(jī)在海洋、海岸帶監(jiān)測(cè)的應(yīng)用,很好地利用了無人機(jī)平臺(tái)并搭載如RTK和LiDAR等設(shè)備,對(duì)關(guān)注區(qū)域進(jìn)行研究。然而,行業(yè)應(yīng)用無人機(jī)或搭載多種探頭無人機(jī)價(jià)格較為昂貴,且在海岸帶等復(fù)雜環(huán)境下存在較大的風(fēng)險(xiǎn),故監(jiān)測(cè)研究的成本較大。僅搭載高清相機(jī)的輕便型無人機(jī)在海岸帶監(jiān)測(cè)中應(yīng)用便可發(fā)揮低成本的優(yōu)勢(shì),做到高頻次的監(jiān)測(cè)。因此,本研究基于低成本的無人機(jī)航拍測(cè)量與影像處理技術(shù),開展臨港新片區(qū)潮灘植被和沙灘演變狀況監(jiān)測(cè)與動(dòng)力成因研究。

1 研究區(qū)域

臨港潮灘,也稱南匯邊灘,位于長江口與杭州灣之間,屬于開敞海岸,地處長江三角洲南翼,近岸水域?qū)俜钦?guī)半日淺海潮,兼具杭州灣潮波性質(zhì),屬于強(qiáng)潮潮灘[10],地貌演變受漲落潮影響。根據(jù)衛(wèi)星圖片顯示以及前人研究,近30年來臨港潮灘整體北沖南淤、向東南方向延伸發(fā)育[11]。潮灘面積變化表現(xiàn)出明顯的夏秋季弱侵蝕、冬春季強(qiáng)淤積的季節(jié)性轉(zhuǎn)化特征,而且潮灘地貌沖刷和沉積具有顯著的空間分布特征[12]。

研究區(qū)是典型的亞熱帶季風(fēng)性氣候,常年溫和濕潤,潮灘濕地生長了許多鹽沼植被,如蘆葦、海三棱藨草、互花米草等。近年來,由于外來物種互花米草的入侵和人類活動(dòng)增加,使得臨港潮灘的植被空間分布發(fā)生變化,在一定程度上改變了潮灘的生態(tài)系統(tǒng),破壞了其原有的生態(tài)平衡。本研究區(qū)域以蘆葦、海三棱藨草和互花米草群落為研究對(duì)象,進(jìn)行臨港潮灘濕地典型植被時(shí)空分布特性研究。

長江三角洲沿岸大多是淤泥質(zhì)海灘。本研究團(tuán)隊(duì)自2018年以來每周對(duì)臨港潮灘進(jìn)行考察,發(fā)現(xiàn)以往臨港潮灘以粉砂、沙質(zhì)粉砂為主,如今潮上帶以砂為主,潮灘有逐漸沙質(zhì)化的趨勢(shì),其海岸的利用類型也隨之發(fā)生變化,成為休閑娛樂的場所。因此,本研究針對(duì)該區(qū)域潮灘短期內(nèi)沙灘的變化開展深入監(jiān)測(cè)研究。

圖1 研究區(qū)概況

2 無人機(jī)作業(yè)流程

無人機(jī)測(cè)量系統(tǒng)以無人駕駛飛行器為飛行平臺(tái),并以高分辨數(shù)字設(shè)備為機(jī)載傳感器,將獲取高分辨率數(shù)據(jù)作為主要應(yīng)用目標(biāo),并實(shí)現(xiàn)航拍影像監(jiān)測(cè)、信息數(shù)據(jù)傳輸以及圖像預(yù)處理等多項(xiàng)功能[13]。本研究選用的無人機(jī)大疆Mavic 2 Pro,其配備了CMOS 2 000萬有效像素的哈蘇鏡頭,可保證在百米級(jí)空中航拍的圖片分辨率,全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)可接收全球定位系統(tǒng)(Global Positioning System,GPS)和格洛納斯(Global Navigation Satellite System,GLONASS)兩種信號(hào)。

無人機(jī)作業(yè)主要包括外業(yè)工作和內(nèi)業(yè)工作。外業(yè)工作主要是進(jìn)行現(xiàn)場調(diào)查,根據(jù)工作需求收集相關(guān)信息,選擇起降場地和參照區(qū)域,規(guī)劃航拍路線,確定空間點(diǎn)位、飛行高度和作業(yè)方式等,并進(jìn)行航拍作業(yè);內(nèi)業(yè)工作主要是對(duì)現(xiàn)場獲取的航拍圖像進(jìn)行檢查、預(yù)處理、拼接和成圖等,再根據(jù)成像區(qū)域色塊進(jìn)行數(shù)據(jù)化分析。

2.1 外業(yè)數(shù)據(jù)獲取

外業(yè)數(shù)據(jù)主要利用旋翼無人機(jī)空中拍攝的方式獲取,操作步驟如下:在飛行軟件地圖上選擇需要航拍的研究區(qū)域,選擇掃描模式自動(dòng)規(guī)劃好航線,設(shè)定飛行高度和飛行速度;拍照模式設(shè)置為等距間隔拍照,相鄰照片的重疊度不低于60%,以保證圖像拼接時(shí)有足夠的同名特征點(diǎn);相機(jī)朝向平行于主航線,保證無人機(jī)航拍過程選擇單個(gè)方向固定不旋轉(zhuǎn);鏡頭設(shè)置以垂直方向?yàn)橹饕较?鏡頭沿垂直方向向下,正對(duì)地面拍攝,根據(jù)航線完成拍攝。

設(shè)定航線需注意:①選擇航拍區(qū)域面積需大于研究區(qū)域,并保存航線,后續(xù)航拍可直接調(diào)用,確保每次航拍圖像范圍一致性;②在航拍設(shè)置區(qū)域內(nèi)選定參照區(qū)域,實(shí)測(cè)計(jì)算其實(shí)際面積,參照區(qū)域選定需滿足在成圖中可清晰選區(qū),確保面積換算的精確度。

2.2 內(nèi)業(yè)數(shù)據(jù)處理

2.2.1 圖像拼接處理

目前在實(shí)際應(yīng)用中,由于專業(yè)全景圖拼接軟件的自動(dòng)化和智能化程度很高,通常40~80張的原始圖像約20 min即可拼接完成。應(yīng)用較多的全景圖拼接軟件包括PTGui、Image Composite Editor、PanoramaStudio和Photoshop等。航拍全景圖拼接主要包括6個(gè)步驟:①導(dǎo)入航拍圖像;②提取圖像特征點(diǎn);③匹配圖像之間的特征點(diǎn),移除錯(cuò)誤匹配;④通過算法和特征點(diǎn)對(duì)圖像進(jìn)行仿射變換;⑤拼接相鄰圖像,重疊部分融合;⑥導(dǎo)出合成的全景圖。對(duì)于拼接錯(cuò)位變形的區(qū)域,進(jìn)行手動(dòng)調(diào)整局部特征點(diǎn),增加人工特征點(diǎn),消除不匹配的特征點(diǎn),以保證全景圖的質(zhì)量[14]。

2.2.2 圖像數(shù)據(jù)化分析

(1)圖像初步處理:所有航拍圖像統(tǒng)一分辨率,通過快速選擇工具(通過查找和追蹤圖像中的邊緣來創(chuàng)建選區(qū))選擇航拍全景圖中的海水和道路等區(qū)域并刪除,從而降低在色彩范圍內(nèi)劃定區(qū)域的誤差。

(2)外業(yè)中在現(xiàn)場選取參照區(qū)域?qū)崪y(cè)其面積為S1(圖2),通過圖片處理軟件獲取圖像中該實(shí)測(cè)區(qū)域的直方圖數(shù)據(jù)記錄其像素值為N1。

圖2 現(xiàn)場選定參照區(qū)

(3)不同區(qū)域的像素:獲取全景圖的直方圖數(shù)據(jù),查看并記錄全景圖的像素值。通過選取單個(gè)位置的顏色或某一色彩范圍,顯示圖像中在此色彩范圍的區(qū)域,從而提取不同色域圖像并記錄其像素值,根據(jù)植被區(qū)域的屬性提取植被區(qū)域,根據(jù)沙質(zhì)區(qū)域的屬性提取沙灘區(qū)域。

(4)通過參照區(qū)實(shí)測(cè)面積和像素值比例的換算,重復(fù)選擇統(tǒng)計(jì)可分別獲取不同色域的面積。例如,植被區(qū)域像素為N2,植被區(qū)域換算后的面積為S2,則S2=N2×S1/N1。

2.3 精度分析與評(píng)估

因現(xiàn)場觀測(cè)受風(fēng)速、空氣質(zhì)量、光照等因素的影響,航拍精度略有不同。為了驗(yàn)證無人機(jī)航拍的換算面積精度,本研究利用岸邊水工建筑物(如,堤壩、石坡等),選取頂部平整,且有垂直角點(diǎn)的區(qū)域作為檢核區(qū)域[15](圖3)。通過圖像換算得出的面積,與實(shí)地測(cè)量的面積進(jìn)行對(duì)比,其面積誤差如表1所示。

圖3 選取航拍圖中檢核區(qū)域的像素值

表1 像素?fù)Q算面積與測(cè)量面積對(duì)比

從表1中可見,對(duì)8個(gè)不同時(shí)期的航拍圖,以實(shí)測(cè)面積為標(biāo)準(zhǔn)計(jì)算了像素?fù)Q算面積方法的誤差率,計(jì)算的總平均誤差率為6.67%,總體來看無人機(jī)精度能夠滿足海岸監(jiān)測(cè)的需要。

3 結(jié)果與討論

3.1 無人機(jī)監(jiān)測(cè)結(jié)果的時(shí)空分布特征

從研究區(qū)連續(xù)監(jiān)測(cè)的不同時(shí)期植被和沙灘圖像,可以直觀地看出兩者的面積變化趨勢(shì)。2019年9月26日,東西兩側(cè)沙灘呈狹長條帶狀沿水邊線分布,中部呈條狀帶垂直于岸線,植被東西兩側(cè)均有分布,東側(cè)植被較淺,西側(cè)植被較深;2019年10月3日,潮灘沿水邊線的沙灘受侵蝕消失,僅部分東側(cè)沙灘分布范圍變窄,西側(cè)植被大面積減少;2019年11月12日,東側(cè)沙灘向沿水邊線收縮,分布形狀更加狹窄,中部和西側(cè)沙灘持續(xù)向岸發(fā)展,西側(cè)植被基本枯萎消失;2019年12月11日,東側(cè)沙灘呈狹長條帶狀沿水邊線分布,中部沙灘受侵蝕而減少,西側(cè)沙灘相比于2019年9月,向岸發(fā)展呈長條狀,水邊線沙灘被侵蝕逐漸模糊,植被分布無顯著變化;2020年1月6日,沙灘出現(xiàn)淤積,又逐漸恢復(fù)至東西兩側(cè)沿水邊線分布的較為分散的條帶狀,中部呈垂直岸線的條狀,植被分布仍無顯著變化;2020年4月15日,沙灘持續(xù)淤積,中部沙灘向岸發(fā)展顯著,西側(cè)近岸植被分散出現(xiàn),東側(cè)植被逐漸向海發(fā)育;2020年6月4日,東側(cè)沙灘呈狹長分布,中部垂直岸線條狀分布,西側(cè)沙灘侵蝕消失,東側(cè)植被持續(xù)向海延伸,西側(cè)植被面積減少向岸收縮;2020年8月3日,沙灘呈現(xiàn)明顯淤積,東側(cè)沙灘沿水邊線條狀面積擴(kuò)寬,西側(cè)沙灘淤積出現(xiàn)沿水邊線分布,潮灘植被大面積增加,基本覆蓋水邊線以上潮灘區(qū)域;2020年9月3日,沙灘西側(cè)受到?jīng)_刷呈分散分布,植被分布無顯著變化;2020年10月13日,沙灘西側(cè)持續(xù)沖刷,加劇分散分布,植被中部區(qū)域開始枯黃。

一年期間植被面積先減少后增加,植被在潮灘西側(cè)向海稍有擴(kuò)張,植被更加濃密,其密度增大;沙灘面積不斷增多?？傮w而言,潮灘的植被秋冬季枯萎,春季生長,夏季茂盛,分布在西側(cè)的植被變化最為顯著,東側(cè)近岸植被常綠;沙灘在潮灘東側(cè)主要分布在潮灘水邊線及上部范圍,中西側(cè)的分布范圍較廣,植被枯萎期間,中西部沙灘裸露且易受沖刷向岸推進(jìn),植被茂盛期間,中西部沙灘沿水邊線分布并向海發(fā)展。

此外,無人機(jī)監(jiān)測(cè)結(jié)合潮灘取樣發(fā)現(xiàn),研究區(qū)正由淤泥質(zhì)潮灘向沙灘轉(zhuǎn)變。長期以來臨港潮灘為淤泥質(zhì)潮灘,如今該潮灘逐漸形成沙質(zhì)灘肩,呈現(xiàn)由淤泥質(zhì)潮灘向沙質(zhì)海灘的轉(zhuǎn)變。從表2可知,一年期間臨港潮灘的沉積物粒度整體變粗,中值粒徑(D50)由2.84~2.9Φ,變?yōu)?.71~2.77Φ。分選性(S)變好,且砂含量增多,結(jié)合一年前后對(duì)灘肩航拍觀察,該區(qū)域的沉積物表現(xiàn)向海擴(kuò)散和向西側(cè)輸移的現(xiàn)象,波流作用沖刷東側(cè)岸線,并使泥沙懸浮隨落潮流向海擴(kuò)散。

表2 一年后潮灘沉積物變化特征

3.2 植被和沙灘對(duì)監(jiān)測(cè)期間內(nèi)事件的響應(yīng)

在風(fēng)暴等極端事件的影響下,潮灘可在短期內(nèi)產(chǎn)生大尺度變形,引起沉積物和植被分布發(fā)生顯著變化。在研究的一年期間,臨港潮灘共經(jīng)歷了6次臺(tái)風(fēng)影響以及多年一遇的洪水等極端事件(表3)。通過事件前后植被和沙灘變化的監(jiān)測(cè)可以揭示極端事件的影響。

表3 不同事件的時(shí)間及水文特征

3.2.1 沙灘變化

從圖4沙灘面積變化可以看出,在2019年9—10月初連續(xù)兩個(gè)臺(tái)風(fēng)影響期間,潮水漫過潮上帶,淹沒潮灘,在波浪卷動(dòng)下泥沙掀動(dòng)并隨落潮離岸輸移,潮灘沿水邊線分布的沙灘受侵蝕消失,沙灘面積減少4 160 m2。臺(tái)風(fēng)過后,沙灘面積逐漸增加,11月沙灘面積逐漸減少,直至翌年1月潮灘沙灘面積開始不斷增加。

圖4 不同時(shí)期植被和沙灘面積時(shí)序變化

由于2020年6月以來長江流域大范圍持續(xù)性強(qiáng)降雨,出現(xiàn)了較為嚴(yán)重的汛情[16-17],洪水期挾帶大量泥沙向海輸移,且波高不大,沒有明顯的增浪(表3),在利于淤積的條件下,從6月4日—8月3日沙灘面積驟增,增加量約為22 060 m2(圖4),潮灘東側(cè)沿水邊線向海延伸淤積,西側(cè)主要沿水邊線淤積。

2020年8—9月,受連續(xù)4個(gè)臺(tái)風(fēng)的影響,8月3日與9月3日對(duì)比沙灘面積減少10 776 m2,在此期間單個(gè)臺(tái)風(fēng)事件導(dǎo)致沙灘面積的最大減少量為20 199 m2,對(duì)比2020年8月和9月的沙灘圖像,一個(gè)月臺(tái)風(fēng)期沙灘面積變化顯著,相比于東側(cè),西側(cè)沖刷較為劇烈。即使受臺(tái)風(fēng)影響,沙灘面積仍呈增長趨勢(shì),在無臺(tái)風(fēng)影響時(shí)沙灘面積逐漸增加,增加量為3 509 m2,臺(tái)風(fēng)過后9月3—9日沙灘面積增加了13 436 m2。

由此可見,一年之中研究區(qū)經(jīng)歷高頻的極端事件,在沙灘面積變化量中,一次臺(tái)風(fēng)的沖刷量和多年一遇的洪水期引起的淤積量最為顯著,且在臺(tái)風(fēng)事件之后該區(qū)域的恢復(fù)周期很長,甚至在恢復(fù)周期內(nèi),連續(xù)的臺(tái)風(fēng)事件導(dǎo)致此區(qū)域持續(xù)沖刷。潮灘不同區(qū)域?qū)κ录捻憫?yīng)也有差異,潮灘西側(cè)的沙灘受臺(tái)風(fēng)影響更為顯著,表現(xiàn)為沙灘向岸推移或者沙灘沖刷消失;洪水期間東側(cè)沙灘的淤積現(xiàn)象更為顯著,表現(xiàn)為沙灘沿水邊線向海延伸。

3.2.2 植被變化

根據(jù)前人對(duì)鹽沼植被生長情況的研究,蘆葦在上海地區(qū)表現(xiàn)為3月出芽,4—6月大量生長,11月結(jié)果,12月后才逐漸枯萎,蘆葦抗寒能力較強(qiáng),冬季也有地面植株生長的可能性[18]。互花米草每年2月開始返青萌芽長出小苗,9月種子成熟期,10—11月后開始逐漸大面積變黃枯萎,至翌年1—2月莖稈大面積倒伏[19-20]。海三棱藨草生長周期主要表現(xiàn)為4月地面植株發(fā)芽開花,8—11月植株結(jié)果后枯萎[21-23]。

結(jié)合各類植被生長周期和空間分布特征,2019年9月—2020年1月區(qū)域的植被主要為海三棱藨草和互花米草,以及蘆葦,11月至翌年1月(互花米草和海三棱藨草枯萎期)植被面積大致為潮灘上的蘆葦面積,平均為16 107 m2。通過不同植被生長階段的植被面積相減,2019年9月臺(tái)風(fēng)前植被面積減去蘆葦植被面積可得出互花米草和海三棱藨草枯萎的面積為34 862 m2,與2020年9月植被面積相減即可得到潮灘植被一年間增加的面積為40 729 m2,同比增長約79.9%,對(duì)比兩年9月的植被分布,植被在潮灘西側(cè)向海稍有擴(kuò)張,潮灘整體的植被密度明顯增大,表明臨港潮的灘植被不斷增加。

海三棱藨草和互花米草2019年9—10月枯萎期間,由于連續(xù)兩個(gè)強(qiáng)臺(tái)風(fēng)潮汐的影響,地面植株被沖走,潮灘植被面積大幅度下降,減少量約為29 515 m2,從分布上來看潮灘西側(cè)植被減少最為顯著;2020年8—9月,在連續(xù)4個(gè)臺(tái)風(fēng)事件下植被面積減少量約為12 217 m2,可知臺(tái)風(fēng)事件對(duì)植被生長具有較大的影響。

2020年7—8月洪水期間,長江口水流攜帶大量泥沙,提供充足的水分和有機(jī)物質(zhì),適逢植被拔節(jié)期間,因此7—8月臨港潮灘植被生長率增快,植被面積驟增,基本覆蓋水邊線以上的潮灘區(qū)域,面積增加97 306 m2。

由此看出,極端事件對(duì)植被的影響十分顯著,適逢植被特定的生長周期,兩者的疊加效應(yīng)為正相關(guān)狀態(tài)時(shí),則對(duì)植被生長和消亡具有促進(jìn)作用。

4 結(jié)論

傳統(tǒng)的人工測(cè)量和遙感衛(wèi)星圖片無法在特定極端天氣獲取數(shù)據(jù)和圖像,臺(tái)風(fēng)事件云層也會(huì)影響衛(wèi)星圖片質(zhì)量和GPS-RTK信號(hào),而無人機(jī)在實(shí)際操作中可采集所需圖像數(shù)據(jù),從而分析不同事件對(duì)于潮灘植被和沙灘的影響。本研究結(jié)合實(shí)測(cè)面積和航拍圖像像素值,計(jì)算不同區(qū)域的面積,并進(jìn)行了精度檢驗(yàn)。結(jié)果能夠滿足潮灘演變監(jiān)測(cè)的需要,達(dá)到了預(yù)期效果,形成了一套有效的快速獲取海灘植被和沙灘面積數(shù)據(jù)的技術(shù)流程,這是無人機(jī)技術(shù)在海洋行業(yè)中應(yīng)用的一次嘗試,表明低成本、高效率且應(yīng)急性強(qiáng)的無人機(jī)航拍對(duì)于高頻事件影響潮灘形態(tài)變化的研究具有很好的時(shí)效性和可推廣性。

本研究基于無人機(jī)技術(shù)開展臨港潮灘監(jiān)測(cè),獲取了一年間潮灘演變的圖像,提取植被和沙灘面積圖像,分別計(jì)算其面積,通過結(jié)合各類植被生長周期發(fā)現(xiàn),臨港潮灘蘆葦主要生長于潮灘東側(cè)近岸,面積約為16 107 m2,互花米草和海三棱藨草主要生長于潮灘中西側(cè),面積約為34 862 m2,潮灘植被的面積在一年間增加了40 729 m2,表明在臨港潮灘的植被分布范圍和密度不斷增加。

一年的監(jiān)測(cè)表明,臨港潮灘植被和沙灘面積在不同事件影響下的變化也不同,在臺(tái)風(fēng)期間對(duì)植被和沙灘的沖刷,導(dǎo)致其面積皆呈現(xiàn)減少現(xiàn)象,而在洪水期間的來沙和降雨量增多導(dǎo)致植被和沙灘面積驟增;潮灘沙灘面積呈現(xiàn)增加的趨勢(shì),并且該潮灘顯示從淤泥質(zhì)潮灘逐漸轉(zhuǎn)變?yōu)樯迟|(zhì)海灘,而沙質(zhì)海灘正是上海稀缺的海岸帶資源,這為在臨港新片區(qū)海洋旅游資源開發(fā)提供了一個(gè)新的選擇。