李陽(yáng)東 ，李仁虎，常亮

1.上海海洋大學(xué)海洋科學(xué)學(xué)院，上海 201306；2.農(nóng)業(yè)農(nóng)村部大洋漁業(yè)開(kāi)發(fā)重點(diǎn)實(shí)驗(yàn)室，上海 201306；3.大洋漁業(yè)資源可持續(xù)開(kāi)發(fā)教育部重點(diǎn)實(shí)驗(yàn)室，上海 201306；4.上海市河口海洋測(cè)繪工程技術(shù)研究中心，上海 201306

近年來(lái)浙江省經(jīng)濟(jì)發(fā)展迅速，海洋開(kāi)發(fā)的速度和廣度都有很大程度的提高。赤潮是浙江近岸海域的主要海洋災(zāi)害，中國(guó)最早有記錄的赤潮就是1933年發(fā)生在浙江鎮(zhèn)海、臺(tái)州、象山石浦一帶的夜光藻暴發(fā)。進(jìn)入21世紀(jì)，浙江海域赤潮頻發(fā)，并呈現(xiàn)范圍不斷擴(kuò)大的趨勢(shì)，浙江近岸的舟山、南麂等海域更是高發(fā)區(qū)域。浙江東部沿海赤潮在1—10月均可能發(fā)生，多發(fā)于5—8月，大部分赤潮事件集中高發(fā)于5—6月。因此，對(duì)浙江近岸海域海洋赤潮災(zāi)害研究的要求程度也越來(lái)越高。

受制于水質(zhì)監(jiān)測(cè)起步時(shí)間和技術(shù)發(fā)展程度，早期的赤潮監(jiān)測(cè)主要借助于肉眼、氣味的表觀識(shí)別或者人工對(duì)監(jiān)測(cè)海域的水樣進(jìn)行采集，然后帶回實(shí)驗(yàn)室分析各種理化指標(biāo)，再根據(jù)相應(yīng)的閾值對(duì)當(dāng)前赤潮階段以及未來(lái)時(shí)間段內(nèi)赤潮暴發(fā)趨勢(shì)進(jìn)行預(yù)測(cè)。但這種方法費(fèi)時(shí)費(fèi)力，且僅能夠了解到當(dāng)時(shí)監(jiān)測(cè)海域的赤潮狀況。而近年來(lái)，遙感技術(shù)的發(fā)展為赤潮監(jiān)測(cè)和預(yù)警預(yù)報(bào)提供了一種更為宏觀、及時(shí)的新途徑?；谶b感方法可以通過(guò)衛(wèi)星獲取大面積范圍內(nèi)的時(shí)間序列影像來(lái)對(duì)海洋環(huán)境狀況進(jìn)行監(jiān)測(cè)和反演。就赤潮而言，由于時(shí)間序列遙感影像是以宏觀的尺度來(lái)展現(xiàn)赤潮特性，因此能夠清楚地反映出對(duì)應(yīng)區(qū)域赤潮發(fā)生的現(xiàn)狀、空間分布特征以及其遷移擴(kuò)散趨勢(shì)等情況，獲得其它監(jiān)測(cè)手段難以發(fā)現(xiàn)的信息。

與MODIS相比，GOCI具有高時(shí)間分辨率的特點(diǎn)，國(guó)內(nèi)外許多學(xué)者利用GOCI影像也進(jìn)行了赤潮監(jiān)測(cè)的相關(guān)研究。Choi et al.（2014）基于GOCI數(shù)據(jù)使用CHL-a濃度法對(duì)韓國(guó)東海有害藻類暴發(fā)進(jìn)行了研究；Lou et al.（2014）基于GOCI遙感反射率數(shù)據(jù)使用三波段（442、490、555nm）赤潮指數(shù)法對(duì)中國(guó)東海海域的有害藻類暴發(fā)事件進(jìn)行了研究，并進(jìn)一步研究了有害藻類晝夜垂直遷移的現(xiàn)象；江彬彬等（2017）基于GOCI數(shù)據(jù)使用555 nm波段歸一化離水輻亮度閾值法提取的中國(guó)東海海域的赤潮信息，并以2011年5月、2011年7月和2014年5月的3次東海赤潮事件為例進(jìn)行了研究；張賀等（2017）基于GOCI瑞利散射校正的反射率數(shù)據(jù)使用三波段（442、490、555 nm）赤潮指數(shù)法，通過(guò)對(duì)2014年5月15日、26日和28日的渤海赤潮事件為例，驗(yàn)證了該方法能有效地監(jiān)測(cè)渤海海域的赤潮情況。

浙江近岸海域水體屬于高渾濁度的二類水體，在高渾濁度水體中，由于碎屑和有色可溶性有機(jī)物的吸收作用，相較于可見(jiàn)光和近紅外波段，紫外波段的離水輻亮度可被忽略，因此本文采用He et al.（2012）提出的紫外大氣校正方法對(duì)GOCI的L1B數(shù)據(jù)進(jìn)行大氣校正得到各波段的遙感反射率數(shù)據(jù)，然后進(jìn)一步基于遙感反射率數(shù)據(jù)反演得到總懸浮物濃度（total suspended particulate matter，TPM），并且計(jì)算赤潮指數(shù)（Red Index，RI）。最后通過(guò)設(shè)置不同的TPM和RI的閾值獲得不同閾值下的赤潮提取結(jié)果，通過(guò)統(tǒng)計(jì)分析確定TPM和RI的最佳閾值，并以2015年6月16日發(fā)生在溫州南麂海域和2016年5月58日發(fā)生在蒼南鰲江口附近海域的赤潮事件為例，驗(yàn)證本文方法對(duì)于提取浙江海域赤潮信息的有效性和準(zhǔn)確性。最后，以2011年5月29日溫州海域和2014年5月23日舟山海域的赤潮事件為例，使用本文算法對(duì)赤潮的日變化情況進(jìn)行研究。

1 數(shù)據(jù)與處理方法

1.1 GOCI數(shù)據(jù)

GOCI是搭載在韓國(guó)COMS地球靜止軌道衛(wèi)星上的海洋水色傳感器（Yang et al.，2012），是世界首個(gè)地球靜止軌道海洋水色傳感器，其發(fā)射于2010年6月27日。GOCI設(shè)置有8個(gè)波段，每個(gè)波段可以提供不同的信息（詳細(xì)信息見(jiàn)表1），GOCI的時(shí)間分辨率和空間分辨率分別為1 h和500 m，影像覆蓋區(qū)域面積為2 500 km×2 500 km（Ryu et al.，2012）。

1.2 總懸浮物濃度

浙江近岸海域的海水屬于典型的二類水體，海水較為渾濁，TPM較高。相關(guān)研究表明，TPM高的水體會(huì)在紅色光波段和近紅外波段出現(xiàn)強(qiáng)反射的現(xiàn)象，會(huì)對(duì)赤潮水體的準(zhǔn)確提取產(chǎn)生影響，同時(shí)，TPM高的水體中不易暴發(fā)赤潮（丘仲鋒等，2011）。所以提高赤潮水體提取的準(zhǔn)確性，本研究在提取浙江海域赤潮水體時(shí)需要先找到合適的總懸浮物濃度閾值來(lái)對(duì)高渾濁水體進(jìn)行掩膜。

表1 GOCI各波段信息1)Table1 GOCI bands information

本文使用He et al.（2013）基于GOCI數(shù)據(jù)針對(duì)杭州灣和長(zhǎng)江口海域TPM反演提出的一種算法來(lái)進(jìn)行浙江沿岸海域TPM的反演。該算法是一種經(jīng)驗(yàn)算法，算法公式如下：

其中，Rrs7和Rrs3分別是GOCI第7波段和第3波段的遙感反射率。

1.3 赤潮指數(shù)

赤潮指數(shù)RI最早是由Ahn et al.（2006）提出的，其使用SeaWiFS的443、510、555 nm波段的離水輻亮度組合比值的赤潮指數(shù)方法來(lái)對(duì)韓國(guó)近岸海域水體赤潮信息進(jìn)行提取。之后，Zhao et al.（2016）對(duì)Ahn和Shanmugam的RI模型進(jìn)行了改進(jìn)，使用MODIS離水輻亮度數(shù)據(jù)對(duì)提取阿拉伯灣海域的赤潮信息進(jìn)行了研究。張賀等（2017）基于適用于渤海海域的RI模型提取了渤海海域的赤潮信息，該模型利用了GOCI瑞利校正的反射率數(shù)據(jù)。王芮（2018）基于GOCI數(shù)據(jù)使用CHL-a濃度閾值法和RI閾值法對(duì)于東海海域赤潮監(jiān)測(cè)進(jìn)行了研究，結(jié)果表明RI閾值法比CHL-a閾值法更優(yōu)。

丘仲峰等（2011）通過(guò)對(duì)比正常水體和赤潮水體的實(shí)測(cè)遙感光譜曲線發(fā)現(xiàn)，正常水體與赤潮水體的遙感反射率存在明顯差異，赤潮水體有兩個(gè)反射峰（570—600、680—700 nm），而正常水體只有單峰（500—600 nm）；正常水體在490—555 nm波段遙感反射率的斜率明顯小于赤潮水體。所以，本文使用GOCI第3和第4波段對(duì)應(yīng)的遙感反射率組合比值來(lái)計(jì)算RI，具體如下：

依據(jù)《河北省土地整治項(xiàng)目管理辦法》，對(duì)責(zé)任主體滅失的礦山廢棄地進(jìn)行修復(fù)綠化時(shí)，對(duì)適宜復(fù)墾的礦山廢棄地，實(shí)施土地復(fù)墾工程。復(fù)墾新增的耕地可用于耕地占補(bǔ)平衡，增加的建設(shè)用地指標(biāo)可調(diào)整到異地使用。此外，指標(biāo)流轉(zhuǎn)的收益可用于礦山環(huán)境恢復(fù)治理。在此過(guò)程中，將工礦廢棄地復(fù)墾與耕地占補(bǔ)平衡、土地轉(zhuǎn)用、項(xiàng)目資金整合等多方面政策結(jié)合，適用了相關(guān)政策和支持政策，整合了資金和力量，有效推進(jìn)了示范區(qū)建設(shè)進(jìn)程。

2 結(jié)果與分析

本文首先通過(guò)對(duì)浙江省赤潮數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，討論了浙江省赤潮發(fā)生的時(shí)間規(guī)律，進(jìn)而可以根據(jù)浙江省赤潮發(fā)生的年際和月際規(guī)律，提取典型時(shí)間段的赤潮信息。在提取浙江赤潮信息時(shí)，通過(guò)統(tǒng)計(jì)不同TPM和RI閾值下得到的赤潮信息并進(jìn)一步分析確定浙江海域赤潮提取的TPM和RI最佳閾值組合。以2015年6月和2016年5月浙江沿海區(qū)域發(fā)生的赤潮事件為例，驗(yàn)證了本文赤潮提取算法的效果。

2.1 浙江赤潮時(shí)間特征

浙江近海是赤潮高發(fā)海域，為了解赤潮發(fā)生的頻次和范圍情況，我們對(duì)來(lái)源于2011—2017年浙江省海洋災(zāi)害公報(bào)上該海域的赤潮數(shù)據(jù)按年和月進(jìn)行了統(tǒng)計(jì)。由2011—2017年浙江海洋災(zāi)害公報(bào)可獲得2006—2017年浙江省赤潮發(fā)生的年累計(jì)次數(shù)、面積信息和2011—2017年浙江省赤潮月累計(jì)次數(shù)、面積信息。圖1為2006—2017年浙江省年赤潮發(fā)生的次數(shù)和累計(jì)海域面積。圖2是2011—2017年多年的月赤潮發(fā)生累積次數(shù)和累計(jì)面積。

圖1 2006—2017年浙江省年赤潮發(fā)生次數(shù)與累計(jì)海域面積Fig.1 Interannual accumulative times and area for the red tides of Zhejiang coast between 2006 and 2017

圖2 2011—2017年浙江省各月赤潮發(fā)生次數(shù)與累計(jì)海域面積Fig.2 Monthly accumulative times and area for the red tides of Zhejiang coast between 2011 and 2017

由圖2中可以看出，除2012—2015年間出現(xiàn)波動(dòng)外，浙江省赤潮發(fā)生次數(shù)與累計(jì)海域面積在2006—2012年呈現(xiàn)下降趨勢(shì)，而在2015—2017年又有所上升；圖3顯示，浙江海域的赤潮最早出現(xiàn)在2月，多發(fā)于5、6、7月這3個(gè)月，峰值出現(xiàn)于5月，5月赤潮發(fā)生的累計(jì)次數(shù)和累計(jì)海域面積分別為49次和5 209.9 km2，5月之后赤潮的累計(jì)次數(shù)和累計(jì)海域面積一直下降，赤潮最晚出現(xiàn)在9月。根據(jù)浙江省赤潮發(fā)生的時(shí)間規(guī)律可以看出，春末和夏季是浙江赤潮的高發(fā)期，應(yīng)重點(diǎn)加強(qiáng)這段時(shí)間的赤潮監(jiān)測(cè)和防治。

圖3 2015年6月16日11:30浙江海域赤潮提取結(jié)果Fig.3 The result of red tide extraction for Zhejiang coast in June 16, 2015

2.2 赤潮提取結(jié)果

本文首先使用紫外大氣校正方法獲得GOCI的遙感反射率數(shù)據(jù)，進(jìn)而根據(jù)公式（1）（2）反演得到TPM數(shù)據(jù)，同時(shí)根據(jù)公式（3）計(jì)算得到RI。通過(guò)設(shè)置不同的TPM和RI閾值提取了2011—2016年浙江沿海的赤潮信息，與赤潮公報(bào)中的赤潮信息對(duì)比，發(fā)現(xiàn)對(duì)于浙江沿岸海域TPM的最佳閾值是17 mg·L-1，RI的最佳閾值是0.23，表2是該閾值下的幾次赤潮事件統(tǒng)計(jì)結(jié)果。在這些赤潮事件發(fā)生的過(guò)程中，對(duì)應(yīng)海域GOCI影像的質(zhì)量較好，云覆蓋范圍較少，有效像元數(shù)較多。從表2中可以看出，本文算法提取赤潮面積的誤差范圍為20.45%—127.34%，平均面積誤差為64.86%。

表2 浙江海域赤潮提取結(jié)果統(tǒng)計(jì)Table 2 Statistic for the extraction results of red tides occurred in Zhejiang coast

前面對(duì)提取的赤潮面積進(jìn)行了定性統(tǒng)計(jì)評(píng)價(jià)，為進(jìn)一步驗(yàn)證本文赤潮提取算法的有效性，分別以2015年6月和2016年5月浙江沿海區(qū)域發(fā)生的赤潮事件為例對(duì)所提取赤潮發(fā)生區(qū)域的準(zhǔn)確性進(jìn)行了驗(yàn)證。根據(jù)赤潮公報(bào)的資料，2015年6月12—21日溫州南麂海域（赤潮事件1）發(fā)生最大面積為390 km2的赤潮，赤潮優(yōu)勢(shì)種為多紋膝溝藻；2016年05月22—30日蒼南鰲江口至霞關(guān)海域（赤潮事件2）發(fā)生最大面積為100 km2的赤潮，優(yōu)勢(shì)種為東海原甲藻。使用本文前述赤潮提取方法，對(duì)這兩個(gè)赤潮事件的赤潮范圍進(jìn)行了提取，其結(jié)果分別如圖3和圖4所示，圖中白色代表無(wú)效值（云覆蓋區(qū)域），藍(lán)色代表正常水體，紅色代表赤潮水體。由圖可見(jiàn)，針對(duì)這兩個(gè)赤潮事件所提取得到的赤潮發(fā)生區(qū)域與公報(bào)中的位置一致。

在現(xiàn)有基于GOCI影像的浙江海域赤潮監(jiān)測(cè)的研究中，與赤潮公報(bào)中的赤潮面積對(duì)比時(shí)，提取面積一般都會(huì)存在偏大的情況，如江彬彬等（2017）基于GOCI數(shù)據(jù)提取的浙江海域的赤潮面積誤差最大達(dá)到了525.00%；王芮（2018）使用GOCI數(shù)據(jù)以赤潮指數(shù)法進(jìn)行東海海域赤潮面積提取時(shí)，其面積誤差范圍為59.00%—105.00%，平均誤差為73%。如表2所示，本文算法提取的赤潮面積誤差范圍為20.45%—127.34%（平均面積誤差為64.86%）。雖然對(duì)于面積較小的赤潮事件（如2011年5月29日—6月1日的赤潮事件），其結(jié)果誤差相對(duì)較大（面積誤差為127.34%），但仍在可接受的范圍內(nèi)?？傮w上，通過(guò)赤潮面積和赤潮發(fā)生海域的提取結(jié)果分析可知，本文使用的赤潮提取方法能有效地提取浙江海域的赤潮信息，可用于浙江海域的赤潮監(jiān)測(cè)。

圖4 2016年5月28日11:30浙江海域赤潮提取結(jié)果Fig.4 The result of red tide extraction for Zhejiang coast in May 28, 2016

3 赤潮日變化特征

為了探究浙江海域赤潮的日變化特征，分別提取了2011年5月29日溫州海域和2014年5月23日舟山海域的每日8幅GOCI影像數(shù)據(jù)的赤潮信息，結(jié)果分別如圖5、6所示。這兩日赤潮發(fā)生海域的遙感影像質(zhì)量較好，有利于赤潮信息的提取和赤潮日變化情況的分析。圖7是以2011年5月29日發(fā)生在溫州海域和2014年5月23日發(fā)生在舟山海域?yàn)榇淼恼憬Ｓ蚋鲿r(shí)刻赤潮面積的統(tǒng)計(jì)信息。

圖5 2011年5月29日溫州海域赤潮日變化Fig.5 Red tide daily change of Wenzhou coastal area in May 29, 2011

圖6 2014年5月23日舟山海域赤潮日變化Fig.6 Daily change of red tide occurred in Zhoushan coastal area in May 23, 2014

圖7 浙江海域赤潮面積日變化Fig.7 Red tide daily area change of Zhejiang coastal area

如圖5所示，2011年5月29日發(fā)生于溫州海域的赤潮在09:30出現(xiàn)于洞頭西南部海域，隨后主要向南擴(kuò)散移動(dòng)，在10:30面積有所減少，但之后一直到13:30面積都在增加，并保持向南的擴(kuò)散，到14:30時(shí)面積下降，并且赤潮影響范圍縮減至霞關(guān)附近海域，在08:30和15:30時(shí)未提取到赤潮信息。圖6是2014年5月23日發(fā)生在舟山海域的赤潮提取結(jié)果，在08:30時(shí)刻提取到小范圍的赤潮區(qū)域，09:30赤潮在枸杞至岱山海域零星分布，之后赤潮向南北方向擴(kuò)散，影響面積逐漸增加，在12:30達(dá)到最大值，之后赤潮面積減少，到15:30時(shí)僅有少量出現(xiàn)于東極島和舟山本島以東海域。圖7中，虛線是2011年5月29日發(fā)生于溫州海域的赤潮面積日變化曲線，在09:30—11:30曲線出現(xiàn)波動(dòng)，但在08:30—13:30整體趨勢(shì)是增加的，在13:30達(dá)到最大值，最大面積為114.25 km2，13:30—15:30曲面一直下降；實(shí)線是2014年5月23日發(fā)生在舟山海域的赤潮面積日變化曲線，可以看出在08:30—12:30赤潮面積持續(xù)增加，在12:30面積達(dá)到最大值，為475.25 km2，特別是在09:30—11:30之間，赤潮面積增長(zhǎng)速度較高，隨后在12:30—15:30赤潮面積又以較快的速度減少。

圖7中的兩條曲線表明兩天的赤潮發(fā)生區(qū)域?qū)?yīng)的面積都呈明顯的先增后減的變化規(guī)律且在12:30—13:30（可認(rèn)為是正午時(shí)刻）左右達(dá)到最大值；進(jìn)一步觀察圖中兩條曲線的變化速度和赤潮面積大小發(fā)現(xiàn)，曲線在08:30—11:30面積增長(zhǎng)速度較快，在14:30—15:30面積又以較大的幅度減小。圖7反映出的赤潮面積的日變化特征與浮游植物晝夜垂直遷移的規(guī)律一致，赤潮面積的日變化可能是浮游植物晝夜垂直遷移導(dǎo)致的（Lou et al.，2014；張賀等，2017）。需要指出的是，在08:30和15:30都有提取的赤潮面積為0的情況，說(shuō)明本文算法在這兩個(gè)時(shí)段提取赤潮面積偏小，具體原因有待進(jìn)一步討論驗(yàn)證。

4 總結(jié)與討論

本文主要利用赤潮統(tǒng)計(jì)資料和具有高時(shí)空分辨率的GOCI影像對(duì)浙江海域赤潮發(fā)生時(shí)空規(guī)律和赤潮信息提取展開(kāi)了研究，并基于赤潮指數(shù)模型提出了一種考慮TPM的赤潮提取方法。首先，通過(guò)統(tǒng)計(jì)2011—2017年浙江海洋災(zāi)害公報(bào)中的赤潮災(zāi)害數(shù)據(jù)，分析了浙江赤潮災(zāi)害發(fā)生的年際和月際變化規(guī)律。在基于GOCI數(shù)據(jù)提取浙江海域赤潮信息時(shí)，在利用紫外大氣校正方法對(duì)GOCI數(shù)據(jù)進(jìn)行大氣校正后，利用得到的遙感反射率反演計(jì)算TPM和RI。然后通過(guò)對(duì)比根據(jù)TPM和RI組合提取的赤潮信息和對(duì)應(yīng)的海洋公報(bào)上的赤潮數(shù)據(jù)，對(duì)浙江海域赤潮提取的TPM和RI的最佳閾值組合進(jìn)行確定并選用案例進(jìn)行驗(yàn)證。最后，使用本文方法以分別發(fā)生在溫州和舟山的兩次赤潮事件為例研究了浙江海域赤潮的日變化情況。通過(guò)研究，形成以下結(jié)論。

（1）與2002—2015年趨勢(shì)相反，2015—2017年來(lái)浙江省赤潮發(fā)生呈上升的趨勢(shì)。同時(shí)具有明顯的月際變化規(guī)律。浙江海域赤潮最早會(huì)出現(xiàn)在2月，多發(fā)于5、6、7月這3個(gè)月，峰值出現(xiàn)于5月，春末和夏季是浙江赤潮的高發(fā)期。

（2）通過(guò)設(shè)置TPM和RI閾值提取浙江海域赤潮信息，實(shí)驗(yàn)結(jié)果表明TPM和RI的最佳閾值分別為17 mg·L-1和0.23。對(duì)于2011—2016年5次典型赤潮事件提取結(jié)果的面積誤差為20.45%—127.34%，平均面積誤差為64.86%。同時(shí)以2015年6月和2016年5月溫州海域的兩個(gè)赤潮事件為例對(duì)赤潮發(fā)生位置提取結(jié)果的準(zhǔn)確性進(jìn)行了驗(yàn)證，結(jié)果表明本文方法能較好地提取浙江海域的赤潮信息。

（3）以2011年5月29日溫州海域和2014年5月23日舟山海域的赤潮事件為例研究了赤潮的日變化情況，發(fā)現(xiàn)赤潮面積變化存在明顯的日變化特征，赤潮面積先增后減，在正午時(shí)刻左右達(dá)到最大值，赤潮面積日變化很可能與浮游植物的晝夜垂直遷移有關(guān)。

本文在提取赤潮信息時(shí)使用了TPM和RI這兩個(gè)參數(shù)，下一步將會(huì)結(jié)合其它與赤潮相關(guān)的海洋環(huán)境參數(shù)提高赤潮監(jiān)測(cè)的準(zhǔn)確性，進(jìn)而在赤潮監(jiān)測(cè)的基礎(chǔ)上建立赤潮預(yù)警預(yù)報(bào)模型。同時(shí)，也可進(jìn)一步研究該方法對(duì)于中國(guó)其他海域赤潮監(jiān)測(cè)的適用性。此外，由于海洋水色傳感器獲取的遙感數(shù)據(jù)易受云的影響，接下來(lái)將會(huì)結(jié)合一些云檢測(cè)算法，減小云對(duì)提取結(jié)果準(zhǔn)確性的影響。本文方法在提取08:30和15:30這兩個(gè)時(shí)段的赤潮信息時(shí)存在面積偏小的情況，具體原因有待進(jìn)一步討論驗(yàn)證并進(jìn)行改進(jìn)。