袁超，張靖宇，肖潔，傅明珠，張學(xué)雷*，崔廷偉，王宗靈

( 1. 自然資源部第一海洋研究所 海洋生態(tài)環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室，山東 青島 266061；2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室 海洋生態(tài)與環(huán)境科學(xué)功能實(shí)驗(yàn)室，山東 青島 266237)

1 引言

自2007 年以來，以滸苔（Ulva prolifera）為肇事種的黃海綠潮連年周期性暴發(fā)，對(duì)山東南岸的海洋生態(tài)環(huán)境、水產(chǎn)養(yǎng)殖和濱海旅游造成重大損失[1-2]。黃海綠潮在4-5 月份發(fā)生于蘇北淺灘，隨季風(fēng)和海流向北漂移并快速生長(zhǎng)，覆蓋和分布面積逐漸擴(kuò)大，在6 月份到達(dá)山東南岸堆積成災(zāi)。滸苔具有高生長(zhǎng)率[3]、強(qiáng)漂浮力[4]和多樣化的繁殖方式[5-6]，在適宜條件下會(huì)大規(guī)模增殖[1,7]。漂浮綠藻生物量在蘇北淺灘僅為數(shù)千噸[8-9]，而到達(dá)山東南岸后可達(dá)數(shù)百萬噸。在綠潮形成規(guī)模后打撈，防控成本和難度將大大增加。因此，應(yīng)重視綠潮的早期防控工作。

現(xiàn)場(chǎng)調(diào)查、遙感追溯和數(shù)值模擬等手段均證實(shí)黃海綠潮發(fā)源于蘇北淺灘[1,9-14]。蘇北淺灘紫菜養(yǎng)殖筏架上的附生綠藻脫落入海是南黃海大規(guī)模綠潮的重要來源[12-13]。在早期發(fā)生過程中，漂浮綠藻由零星綠藻、小斑塊逐漸發(fā)展為條帶狀的大斑塊[1,7,13]。然而，由于蘇北淺灘水文環(huán)境復(fù)雜且面積較大，船測(cè)調(diào)查手段難以實(shí)時(shí)和大范圍地對(duì)漂浮綠藻進(jìn)行監(jiān)測(cè)。對(duì)于漂浮綠藻在蘇北淺灘海域的分布、發(fā)生和發(fā)展過程仍缺乏精細(xì)的刻畫，制約了綠潮早期防控工作的研究。

衛(wèi)星遙感方法具有時(shí)效性好和覆蓋范圍廣的優(yōu)點(diǎn)，是綠潮動(dòng)態(tài)觀測(cè)的重要手段。目前，國內(nèi)外學(xué)者基于多源衛(wèi)星，發(fā)展出一系列的綠潮提取算法，開展了綠潮覆蓋面積、影響面積和漂移路徑的研究，并嘗試進(jìn)行綠潮生物量的估算[14-22]。由于漂浮綠藻與陸地植被具有相似的光譜特征，光學(xué)遙感的綠潮提取算法主要基于漂浮綠藻和海水在紅光和近紅外波段的光譜差異。常見的綠潮提取算法包括歸一化植被指數(shù)（Normalized Difference Vegetation Index，NDVI）、增強(qiáng)型植被指數(shù)（Enhanced Vegetation Index，EVI）、漂浮藻類指數(shù)（Floating Algae Index，F(xiàn)AI）和針對(duì)GOCI 衛(wèi)星的漂浮綠藻指數(shù)（Index of Floating Green Algae for GOCI，IGAG）等。其中，NDVI 由于表現(xiàn)穩(wěn)定可靠而被廣泛采用[15-17]。在數(shù)據(jù)源方面，前期研究多基于中低分辨率衛(wèi)星（MODIS、GOCI 等）進(jìn)行綠潮的跟蹤監(jiān)測(cè)，研究區(qū)域集中于鹽城以北海域[14,20]。由于蘇北淺灘漂浮綠藻斑塊面積較小[8]，中低空間分辨率衛(wèi)星難以對(duì)其進(jìn)行識(shí)別且可能會(huì)造成較大誤差[17,22]。近年來，隨著國內(nèi)外高分辨率衛(wèi)星，如哨兵2 號(hào)（Sentinel-2）、高分（GF）系列衛(wèi)星等的運(yùn)行，使研究漂浮綠藻在蘇北淺灘時(shí)空分布規(guī)律成為可能。

本文利用高分辨率衛(wèi)星（以哨兵2 號(hào)為主）數(shù)據(jù)，提取了2018 年黃海綠潮發(fā)生期間蘇北淺灘的漂浮綠藻信息，研究了漂浮綠藻在蘇北淺灘的時(shí)空分布特征，結(jié)合環(huán)境數(shù)據(jù)分析了地形、水溫和風(fēng)向?qū)ζ【G藻分布、發(fā)生和發(fā)展的影響，以期為黃海綠潮的早期防控提供參考。

2 數(shù)據(jù)與方法

2.1 研究海區(qū)

研究區(qū)域位于南黃海西部的蘇北淺灘及其北部（32.4°～34.4°N，120°～122°E）。蘇北淺灘具有似三角洲狀的地形地貌—輻射沙洲（圖1）。輻射沙洲以江蘇東臺(tái)弶港為中心，向東北和東南方向延伸出若干沙脊和潮溝。蘇北淺灘的潮流受輻射沙洲影響顯著。漲、落潮的主流向與潮溝走向一致，潮溝流速大于沙脊[23]。蘇北淺灘是中國條斑紫菜筏式養(yǎng)殖的重要區(qū)域，而養(yǎng)殖筏架上的附生綠藻脫落入海被認(rèn)為是黃海綠潮的重要來源[24]。

2.2 研究數(shù)據(jù)

哨兵2 號(hào)衛(wèi)星是歐洲航天局哥白尼計(jì)劃發(fā)射的多光譜對(duì)地觀測(cè)衛(wèi)星。哨兵2 號(hào)衛(wèi)星的高度為786 km，幅寬為290 km，重訪周期為5 d。在13 個(gè)光譜波段中，可見光（B2，B3，B4）和近紅外（B8）波段的空間分辨率為10 m。哨兵2 號(hào)衛(wèi)星的高時(shí)空分辨率特征為研究漂浮綠藻在蘇北淺灘的分布、發(fā)生和發(fā)展規(guī)律提供了數(shù)據(jù)支撐。

圖 1 蘇北淺灘水深和地形概況Fig. 1 Bathymetric topography of the Subei Shoal

哨兵2 號(hào)影像為免費(fèi)開放獲取，下載于歐洲航天局網(wǎng)站（https://scihub.copernicus.eu/dhus/）。本文篩選了2018 年黃海綠潮暴發(fā)期的少云影像（5 月14 日、5 月24 日、6 月3 日、6 月8 日、6 月13 日、6 月23 日、7 月3 日和7 月13 日），研究漂浮綠藻在蘇北淺灘的時(shí)空分布特征。哨兵2 號(hào)L1C 級(jí)影像是經(jīng)過正射校正和亞像元級(jí)幾何精校的大氣表觀反射率產(chǎn)品。L1C 級(jí)數(shù)據(jù)可基于歐洲航天局發(fā)布的Sen2cor 插件進(jìn)行輻射定標(biāo)和大氣校正，得到各波段的大氣底層反射率（即L2A 級(jí)產(chǎn)品）。5 月14 日，哨兵2 號(hào)影像未檢測(cè)到漂浮綠藻。5 月24 日，漂浮綠藻條帶在蘇北淺灘南部可見，但該區(qū)域被厚云覆蓋。因此，本研究還選取了2018 年5 月23 日的高分1 號(hào)WFV 影像，其空間分辨率為16 m。該數(shù)據(jù)來源于國家衛(wèi)星海洋應(yīng)用中心數(shù)據(jù)分發(fā)系統(tǒng)（dds.nsoas.org.cn/mainIndex.do）。

為了驗(yàn)證哨兵2 號(hào)影像提取漂浮綠藻信息的可靠性，本研究于2018 年6 月13 日在蘇北淺灘南部海域開展了漂浮綠藻分布狀態(tài)的現(xiàn)場(chǎng)觀測(cè)。

利用美國Remote Sensing Systems（RSS）網(wǎng)站提供的微波+紅外融合海表溫度（Microwave+IR Signal Fusion SST， MISST） 數(shù) 據(jù) 和CCMP（ Cross Calibrated Multi-Platform）海面風(fēng)場(chǎng)數(shù)據(jù)，研究影響漂浮綠藻時(shí)空分布的主要環(huán)境因子。MISST 和CCMP 海面風(fēng)場(chǎng)數(shù)據(jù)的空間分辨率分別為9 km 和0.25°。

2.3 漂浮綠藻提取方法

利用NDVI 指數(shù)進(jìn)行漂浮綠藻信息提?。?/p>

式中，NIR、R為近紅外和紅外波段的反射率。本研究中，NIR采 用哨兵2 號(hào)的B8 波 段（842 nm），R采 用B4 波段（665 nm）。

通過設(shè)定閾值T來確定漂浮綠藻的存在與否。本研究通過NDVI 閾值分割和假彩色合成圖像進(jìn)行目視解譯來共同確定T值。閾值T的理論值為0，但受各種因素（淺灘、水深等）的影響而存在波動(dòng)。本文中T值的范圍為[-0.10，0.05]。

通過計(jì)算綠潮的覆蓋面積（S）和漂浮綠藻分布區(qū)域的外包絡(luò)面的面積（D）來表征綠潮的暴發(fā)規(guī)模和影響范圍。S為NDVI＞T的像元數(shù)N與像元面積的乘積：

式中，r為衛(wèi)星傳感器的空間分辨率。

3 結(jié)果與討論

3.1 哨兵2 號(hào)影像中漂浮綠藻信息的現(xiàn)場(chǎng)驗(yàn)證

通過現(xiàn)場(chǎng)同步調(diào)查，對(duì)2018 年6 月13 日哨兵2 號(hào)影像提取的漂浮綠藻信息進(jìn)行驗(yàn)證（圖2，表1）。結(jié)果表明，哨兵2 號(hào)提取的漂浮綠藻信息能夠與現(xiàn)場(chǎng)數(shù)據(jù)相匹配。漂浮綠藻位于海水表面，在渾濁水體中可以與背景區(qū)分開。1 號(hào)站點(diǎn)存在微小的綠藻斑點(diǎn)，面積小于0.5 m2，哨兵2 號(hào)影像未探測(cè)到漂浮綠藻信息。2 號(hào)站點(diǎn)漂浮綠藻分布較為稀疏，覆蓋度約為10%，影像檢測(cè)到微弱漂浮綠藻信息。3 號(hào)和4 號(hào)站點(diǎn)存在顯著的漂浮綠藻條帶，可以清晰地從影像中看到。

圖 2 2018 年6 月13 日哨兵2 號(hào)衛(wèi)星影像中漂浮綠藻信號(hào)的現(xiàn)場(chǎng)驗(yàn)證Fig. 2 Field verification of floating green algae signal from Sentinel-2 images on June 13, 2 018

表 1 2018 年6 月13 日不同站位哨兵2 號(hào)影像中漂浮綠藻信號(hào)的現(xiàn)場(chǎng)驗(yàn)證Table 1 Field verification of floating green algae signal by Sentinel image at different stations on June 13, 2 018

3.2 漂浮綠藻在蘇北淺灘的時(shí)空分布特征

本文采取NDVI 閾值分割與目視解譯法相結(jié)合的方法，對(duì)2018 年蘇北淺灘海域漂浮綠藻信息進(jìn)行了提?。▓D3）；圖4 反映了漂浮綠藻外包絡(luò)面、覆蓋面積和分布面積隨時(shí)間的變化規(guī)律。

2018 年蘇北淺灘漂浮綠藻的發(fā)展過程可大體劃分為發(fā)生—發(fā)展—維持—消失4 個(gè)階段。5 月23 日，零星漂浮綠藻在蘇北淺灘南部的太陽島附近首次被發(fā)現(xiàn)，覆蓋面積為1.78 km2。6 月3 日，漂浮綠藻的分布范圍迅速擴(kuò)大至整個(gè)蘇北淺灘，覆蓋面積擴(kuò)大了32 倍，分布面積達(dá)到17 830 km2。6 月中下旬，漂浮綠藻發(fā)展進(jìn)入維持階段。6 月13 日的漂浮綠藻的分布特征與6 月3 日相似，但在淺灘以北海域的分布面積顯著增加。6 月23 日，漂浮綠藻在蘇北淺灘以北的影響范圍向東擴(kuò)展到122°E，覆蓋面積（199.28 km2）達(dá)到最大，比5 月23 日約高兩個(gè)數(shù)量級(jí)（圖4）。7 月3 日，蘇北淺灘大部分區(qū)域被厚云覆蓋，但在無云區(qū)仍然可觀測(cè)到漂浮綠藻條帶。7 月13 日，研究區(qū)域已無法檢測(cè)到漂浮綠藻。整個(gè)過程持續(xù)約50 d，與歷年黃海綠潮生消過程相一致[16,25]。

在蘇北淺灘，漂浮綠藻聚集成條帶狀，自紫菜養(yǎng)殖筏架區(qū)向外沿各潮溝呈輻射狀分布。如圖3 所示，漂浮綠藻最早可追溯至淺灘中心筏架養(yǎng)殖區(qū)邊緣，而后沿潮溝延伸擴(kuò)展。綠藻條帶的寬度在10～200 m之間，沿潮溝可斷續(xù)綿延數(shù)十千米。出蘇北淺灘后，漂浮綠藻仍呈條帶狀分布，生物量與覆蓋面積顯著增加，但其分布無明顯規(guī)律，呈散布狀（圖3d，圖3e，圖3f）。6 月3 日，綠藻條帶已覆蓋整個(gè)北部淺灘，其中部分漂浮綠藻從淺灘西北沿岸和東部離岸潮溝漂移出淺灘。淺灘以北海域，漂浮綠藻在6 月3 日至6 月13 日明顯向江蘇近岸聚集，并在34°N 以北出現(xiàn)東西跨越半個(gè)經(jīng)度的綠潮鋒面；6 月23 日漂浮綠藻則向外海擴(kuò)散，從近岸至122°E 連成高覆蓋率、高生物量的大面積綠潮。

3.3 漂浮綠藻時(shí)空分布的環(huán)境驅(qū)動(dòng)機(jī)制

適宜的環(huán)境條件是黃海綠潮暴發(fā)的必要條件。蘇北淺灘及周邊海域較高的營養(yǎng)鹽水平是支撐滸苔快速生長(zhǎng)和繁殖的重要物質(zhì)基礎(chǔ)[26-28]。圖5 展示了研究期間蘇北淺灘及周邊海域的海表溫度和海面風(fēng)場(chǎng)的變化過程。溫度是影響滸苔微觀繁殖體萌發(fā)和附生綠藻群落演替的重要環(huán)境因子[1]。室內(nèi)實(shí)驗(yàn)結(jié)果表明，滸苔的最適生長(zhǎng)溫度范圍為14～26℃[3]，而其微觀繁殖體在10～25℃萌發(fā)率較高[5,6,29]。本研究中，研究區(qū)域的海表溫度在各時(shí)間的空間分布較為一致，自南至北、由近岸向外海逐漸降低。蘇北淺灘的海表溫度最高，通常比同緯度的黃海外海高2℃左右。本研究首次探測(cè)到漂浮綠藻時(shí)蘇北淺灘的海表溫度約為18℃，與往年船測(cè)調(diào)查中該區(qū)域大面積漂浮綠藻出現(xiàn)時(shí)的海表溫度相一致[13]，也與歷年衛(wèi)星遙感首次探測(cè)到漂浮綠藻時(shí)的海表溫度相符[18,21,25]。在維持階段，海表溫度雖然持續(xù)上升，但是仍處于滸苔的最適生長(zhǎng)溫度區(qū)間內(nèi)[3]。

圖 3 2018 年不同日期哨兵2 號(hào)影像的漂浮綠藻（紅點(diǎn)）分布（白色區(qū)域?yàn)樵蒲谀ぃ〧ig. 3 Time-series distribution of the floating green algae (red dots) in the Subei Shoal in 2 018 based on the Sentinel-2 images (the thick clouds are marked by white dots)

在蘇北淺灘，漂浮綠藻自紫菜養(yǎng)殖筏架區(qū)沿潮溝向外輻射分布。這表明蘇北淺灘是黃海綠潮的源頭。前期研究表明，蘇北淺灘的紫菜養(yǎng)殖筏架為滸苔的微觀繁殖體提供了附著基，是黃海綠潮的孕育場(chǎng)所[5-8]。附生綠藻在紫菜養(yǎng)殖筏架回收過程中脫落入海，由附生轉(zhuǎn)變?yōu)槠顟B(tài)，是黃海綠潮的重要來源[12-13,24]。蘇北淺灘獨(dú)特的潮流系統(tǒng)和輻射沙洲地貌特征有利于漂浮綠藻的帶狀聚集。該區(qū)域潮溝與沙脊相間，潮溝由外海向核心區(qū)變窄變淺。相同流量下，潮溝內(nèi)海水流速大于沙脊，有利于漂浮綠藻聚集并向外輸運(yùn)[22-23]。漂浮綠藻在潮溝中匯聚還可能與谷物圈效應(yīng)（Cheerios Effect）有關(guān)[30]。當(dāng)兩個(gè)能夠產(chǎn)生液面變形的物體相互靠近時(shí)，由于系統(tǒng)重力勢(shì)能和表面能傾向于最低，物體之間就會(huì)相互靠近聚集。

相關(guān)研究表明，綠潮漂移路徑主要由風(fēng)場(chǎng)驅(qū)動(dòng)下的表層流場(chǎng)決定[11,18,21]。從對(duì)應(yīng)日期的CCMP 海面風(fēng)場(chǎng)圖（圖5）可以看出，漂浮綠藻條帶在蘇北淺灘以北的走向與盛行風(fēng)向較為一致。5 月中旬，海表風(fēng)場(chǎng)由北風(fēng)向南風(fēng)轉(zhuǎn)變，有利于漂浮綠藻的向北漂移。6 月3-13 日研究區(qū)盛行東南風(fēng)，漂浮綠藻主要在江蘇近岸堆積。6 月23 日，當(dāng)風(fēng)向轉(zhuǎn)變?yōu)槲髂巷L(fēng)，綠藻條帶則向外海擴(kuò)展。在較小空間尺度上，漂浮綠藻仍呈帶狀分布。這可能與拉彌爾環(huán)流有關(guān)，小的拉彌爾環(huán)流逐步合并會(huì)引起滸苔或其他漂浮物的條帶狀聚集[11,31]。

圖 4 2018 年蘇北淺灘海域漂浮綠藻的外包絡(luò)面（a）、覆蓋面積和分布面積（b）變化Fig. 4 Temporal changes of envelops (a), coverage and distribution area (b) of the floating green algae in the Subei Shoal in 2 018

3.4 對(duì)黃海綠潮早期防控的啟示

黃海綠潮具有異地起源、大規(guī)模和長(zhǎng)距離漂移的特征，對(duì)山東半島南岸造成重大的經(jīng)濟(jì)和社會(huì)損失[1]。對(duì)漂浮綠藻在蘇北淺灘分布和發(fā)展過程的充分認(rèn)識(shí)能夠?yàn)橹贫ê蛯?shí)施黃海綠潮早期防控措施提供科學(xué)依據(jù)。

減少綠藻附著和控制附生綠藻入海是控制黃海綠潮暴發(fā)規(guī)模的第一步。通過研發(fā)新型材料或?qū)ΜF(xiàn)有材料進(jìn)行處理以減少綠藻附著的相關(guān)研究已有報(bào)道，但由于成本和操作等問題而難以得到大規(guī)模應(yīng)用[28]。目前，將紫菜養(yǎng)殖網(wǎng)簾直接運(yùn)回陸地處理可以減少附生綠藻入海。然而，仍有大量附著于梗繩上的綠藻通過機(jī)械作業(yè)去除入海，成為附生綠藻脫落入海的主要方式[12,24]。附生綠藻入海后，滸苔由于具有高比生長(zhǎng)率[3]、強(qiáng)漂浮能力[4]和多樣化的繁殖方式，迅速成為漂浮綠藻群落的單一優(yōu)勢(shì)種[1]。據(jù)估計(jì)，漂浮綠藻的初始入海生物量?jī)H為幾千噸，但到達(dá)山東沿岸時(shí)可達(dá)數(shù)百萬噸[8-9,32]。因此，黃海綠潮的早期防控可以起到事半功倍的效果。

本研究中，哨兵2 號(hào)衛(wèi)星影像首次探測(cè)到漂浮綠藻的時(shí)間為5 月23 日，比船測(cè)調(diào)查手段首次發(fā)現(xiàn)零星漂浮綠藻的時(shí)間（4 月25 日）晚了近1 個(gè)月[33]。這可能是由于早期零星漂浮綠藻面積較小，只有當(dāng)匯聚成較大斑塊才能夠從衛(wèi)星影像中解譯出來[13]。因此，在漂浮綠藻的早期監(jiān)測(cè)和預(yù)警過程中，高分辨率衛(wèi)星遙感仍需與船舶和無人機(jī)等手段相結(jié)合。

高分辨率衛(wèi)星（哨兵2 號(hào)和高分系列衛(wèi)星等）能夠提供準(zhǔn)實(shí)時(shí)、精細(xì)的漂浮綠藻分布特征，為漂浮綠藻的前期防控提供數(shù)據(jù)支撐。本研究結(jié)果表明蘇北淺灘的漂浮綠藻受輻射沙洲的約束而在潮溝中聚集成條帶狀，且可能對(duì)黃海綠潮有持續(xù)貢獻(xiàn)。綠潮在潮溝內(nèi)的帶狀聚集有利于開展船舶打撈工作。在出淺灘后，綠藻條帶走向會(huì)隨風(fēng)和流的作用而發(fā)生改變，增加了打撈難度。在向北漂移的過程中，漂浮綠藻會(huì)進(jìn)入快速生長(zhǎng)期，增加了打撈的經(jīng)濟(jì)成本。因此，在蘇北淺灘進(jìn)行漂浮綠藻打撈可能是一種切實(shí)可行、經(jīng)濟(jì)且高效的黃海綠潮防控手段。

圖 5 2018 年5 月到7 月南黃海西部海表溫度（等值線，℃）和海面風(fēng)場(chǎng)（矢量）Fig. 5 Sea surface temperature (contour) and wind field (vectors) in the southwestern Yellow Sea during May to July, 2 018

4 結(jié)論

本文主要采用空間分辨率為10 m×10 m 的哨兵2 號(hào)影像，通過NDVI 閾值分割和目視解譯法相結(jié)合的方法提取2018 年蘇北淺灘的漂浮綠藻信息，結(jié)合地形、MISST 海表溫度和CCMP 海面風(fēng)場(chǎng)數(shù)據(jù)分析影響漂浮綠藻時(shí)空變化的主要環(huán)境因子。結(jié)論如下：

（1）2018 年，漂浮綠藻最早于5 月23 日在淺灘南部被探測(cè)到，在6 月份逐漸向北擴(kuò)大并維持，在7 月中旬消失。

（2）蘇北淺灘的漂浮綠潮聚集成條帶狀，自紫菜養(yǎng)殖筏架區(qū)沿潮溝向外輻射分布，并在黃海綠藻發(fā)展過程中持續(xù)向北輸出。出淺灘后，漂浮綠藻的分布和走向與盛行風(fēng)向密切相關(guān)。

（3）漂浮綠藻在蘇北淺灘的帶狀聚集有利于開展船舶打撈工作。