楊靜，蔡文博，李海

（國家海洋環(huán)境預(yù)報(bào)中心，北京 100081）

1 引言

近年來，我國近海赤潮發(fā)生的頻率和面積均居高不下。2010年我國沿海共發(fā)現(xiàn)赤潮69次，累計(jì)面積10892 km2，赤潮災(zāi)害造成直接經(jīng)濟(jì)損失達(dá)2.06億元【1】。赤潮已成為頻繁侵害我國的主要海洋災(zāi)害之一。目前業(yè)務(wù)化的赤潮預(yù)測工作主要應(yīng)用經(jīng)驗(yàn)、統(tǒng)計(jì)或統(tǒng)計(jì)結(jié)合數(shù)值的方法來對(duì)赤潮發(fā)生的可能性進(jìn)行判斷，但由于赤潮的發(fā)生是物理、化學(xué)、生物等多種因素綜合作用的結(jié)果，直接對(duì)赤潮的發(fā)生時(shí)間和范圍進(jìn)行預(yù)測往往不能滿足赤潮災(zāi)害預(yù)警的要求。赤潮爆發(fā)后常伴有一定范圍轉(zhuǎn)移、擴(kuò)散的發(fā)生，例如2000年5月12—16日，浙江中部臺(tái)州列島附近海域發(fā)生面積為1000 km2的赤潮，18日再次發(fā)現(xiàn)赤潮面積約4560 km2，赤潮生物以具齒原甲藻為主，20日赤潮區(qū)域擴(kuò)展至5800 km2，24日赤潮區(qū)域較5月20日有所北移，面積進(jìn)一步擴(kuò)大[2]。2012年5月17日平潭島附近海域發(fā)現(xiàn)帶狀、塊狀赤潮，5月18—26日漂移至福鼎、霞浦近岸，最大發(fā)生面積約130 km2。通過衛(wèi)星遙感或現(xiàn)場觀測確定其藻種及分布特征，然后在觀測的基礎(chǔ)上由數(shù)值模式對(duì)赤潮的可能影響位置及范圍提前進(jìn)行判斷，能最大程度地減輕赤潮災(zāi)害可能造成的損失。

本文將赤潮漂移預(yù)報(bào)模式與業(yè)務(wù)化海面風(fēng)場和中國海流數(shù)值預(yù)報(bào)相結(jié)合，基于C/S架構(gòu)建立赤潮漂移預(yù)報(bào)系統(tǒng)，通過數(shù)值計(jì)算對(duì)赤潮漂移軌跡和擴(kuò)散狀況進(jìn)行預(yù)測，將會(huì)進(jìn)一步完善現(xiàn)有的赤潮預(yù)報(bào)系統(tǒng)，更有效地服務(wù)于赤潮減災(zāi)工作，為赤潮災(zāi)害的損失評(píng)估提供參考和依據(jù)。

2 赤潮漂移擴(kuò)散預(yù)報(bào)系統(tǒng)介紹

2.1 模塊功能組成

從圖1上看，赤潮漂移擴(kuò)散預(yù)報(bào)系統(tǒng)由赤潮發(fā)生區(qū)域提取模塊，赤潮預(yù)報(bào)輸入輸出模塊，赤潮漂移數(shù)值計(jì)算模塊和赤潮預(yù)報(bào)可視化模塊四部分組成。

圖1 赤潮漂移擴(kuò)散數(shù)值預(yù)報(bào)系統(tǒng)模塊組成

赤潮發(fā)生區(qū)域提取模塊用于加載衛(wèi)星觀測資料，通過遙感解析、人機(jī)交互等方式確定赤潮發(fā)生的多邊形區(qū)域；赤潮預(yù)報(bào)輸入輸出模塊使用戶通過交互配置輸入輸出數(shù)據(jù)文件路徑等參數(shù)，同時(shí)提供讀取輸入示范海域的海流、風(fēng)場和赤潮發(fā)生區(qū)域等數(shù)據(jù)，并將預(yù)報(bào)產(chǎn)品輸出為指定的文件格式；赤潮漂移數(shù)值計(jì)算模塊是赤潮漂移擴(kuò)散預(yù)報(bào)的核心，主要應(yīng)用數(shù)值方法計(jì)算赤潮漂移的軌跡、區(qū)域范圍等數(shù)值預(yù)報(bào)結(jié)果；赤潮預(yù)報(bào)可視化模塊基于GIS開發(fā)圖形顯示系統(tǒng)，實(shí)現(xiàn)赤潮漂移擴(kuò)散預(yù)報(bào)結(jié)果的展示。

2.2 赤潮漂移擴(kuò)散預(yù)報(bào)模式簡介

利用衛(wèi)星遙感、航空遙感、船舶觀測等手段確定監(jiān)測海域赤潮發(fā)生位置和范圍等初始信息，通過赤潮信息數(shù)字化處理輸入赤潮漂移路徑預(yù)報(bào)模塊。模型不考慮赤潮的生態(tài)繁殖擴(kuò)展過程，而只考慮動(dòng)力環(huán)境對(duì)赤潮漂移的影響，由國家海洋環(huán)境預(yù)報(bào)中心現(xiàn)有業(yè)務(wù)化運(yùn)行系統(tǒng)和資料同化分析系統(tǒng)提供風(fēng)場、海流等動(dòng)力環(huán)境，利用獲取的風(fēng)速、風(fēng)向和流速、流向等數(shù)據(jù)計(jì)算赤潮生物團(tuán)的漂移過程，應(yīng)用粒子拉格朗日隨機(jī)游走模式來模擬計(jì)算赤潮生物團(tuán)的擴(kuò)散過程。

赤潮藻團(tuán)在風(fēng)和流共同作用下，只考慮其在水平方向的物理過程，運(yùn)動(dòng)方程如下所示：

式中，va為海流速度，包括潮流、環(huán)流，vd為風(fēng)速，R為經(jīng)驗(yàn)系數(shù)，R×vd(xi,t)為風(fēng)對(duì)赤潮藻團(tuán)的拖拽速度。方程求解方法采用了一階求解，公式如下：

圖2 赤潮漂移擴(kuò)散預(yù)報(bào)業(yè)務(wù)流程圖

3 業(yè)務(wù)運(yùn)行流程

赤潮漂移擴(kuò)散預(yù)報(bào)模塊是基于C/S架構(gòu)建立的，其流程見圖2。C/S即Client/Server（客戶機(jī)/服務(wù)器），是較為流行和通用的軟件架構(gòu)。其原理是通過將任務(wù)合理分配到Client端和Server端，大大降低了系統(tǒng)的通訊開銷。C/S架構(gòu)服務(wù)器端程序主要提供數(shù)據(jù)管理、數(shù)據(jù)共享，而用戶的程序則部署在客戶端，主要完成用戶的各種具體的數(shù)據(jù)處理任務(wù)。

客戶端：預(yù)報(bào)員用客戶端程序加載研究區(qū)遙感影像，通過人機(jī)交互的方式勾出赤潮發(fā)生區(qū)域矢量信息。同時(shí)，通過客戶端程序界面配置風(fēng)場、海流等初始場數(shù)據(jù)文件的網(wǎng)絡(luò)地址和相關(guān)參數(shù)信息。赤潮發(fā)生區(qū)域矢量信息和初始化配置參數(shù)信息存放在xml文件中，并通過ftp服務(wù)上傳到服務(wù)器端。通過服務(wù)器端的高性能計(jì)算機(jī)計(jì)算赤潮漂移擴(kuò)散數(shù)值預(yù)報(bào)結(jié)果后，下載標(biāo)準(zhǔn)的GIS柵格數(shù)據(jù)，實(shí)現(xiàn)赤潮漂移擴(kuò)散可視化并生成圖片、動(dòng)畫等預(yù)報(bào)產(chǎn)品。

服務(wù)器端：響應(yīng)預(yù)報(bào)員客戶端的預(yù)報(bào)服務(wù)請(qǐng)求。首先解析上傳的xml文件獲取風(fēng)場、海流初始場的網(wǎng)絡(luò)地址相關(guān)參數(shù)信息；獲取赤潮發(fā)生的空間信息，并以多組多邊形的形式表示。加載必要的輸入數(shù)據(jù)，在服務(wù)器端完成赤潮漂移擴(kuò)散的數(shù)值模式計(jì)算，并生成netcdf格式的預(yù)報(bào)結(jié)果。通過數(shù)值預(yù)報(bào)后處理功能，將netcdf格式的預(yù)報(bào)結(jié)果轉(zhuǎn)換為標(biāo)準(zhǔn)的GIS柵格數(shù)據(jù)。

圖3 2007年7月24—27日長江口附近海域發(fā)生赤潮的連續(xù)遙感監(jiān)測衛(wèi)星圖

表1 赤潮藻團(tuán)24 h和48 h的漂移擴(kuò)散預(yù)測

4 模型的后報(bào)試驗(yàn)

以2007年7月22日—8月6日長江口附近發(fā)生赤潮為例，進(jìn)行后報(bào)試驗(yàn)。利用建立的赤潮漂移與擴(kuò)散數(shù)值預(yù)報(bào)模型，根據(jù)衛(wèi)星遙感影像選取赤潮藻團(tuán)位置點(diǎn)，組成赤潮范圍包絡(luò)線，在不考慮赤潮的生物、化學(xué)轉(zhuǎn)化過程的前提下，通過模型計(jì)算對(duì)赤潮藻團(tuán)未來24－48 h的漂移軌跡和分布范圍進(jìn)行預(yù)測。

根據(jù)東海海洋環(huán)境保護(hù)通報(bào)，2007年7月24日上午，發(fā)現(xiàn)在 29°54′N 、122°28.92′E，29°54.42′N 、122°36.83′E ，30°02′N 、122°37′E，30°1.92′N 、122°29.17′E四點(diǎn)連線范圍內(nèi)均有赤潮發(fā)生，面積200多平方公里，海水呈紅褐色，呈大面狀分布。同日，經(jīng)現(xiàn)場采樣分析，赤潮生物主要優(yōu)勢(shì)種為扁面角毛藻，無毒性，密度為5.08×106個(gè)/升，次優(yōu)勢(shì)種為膝溝藻屬，無毒性，密度為5.12×105個(gè)/升，另外伴有旋鏈角毛藻、中肋骨條藻、圓篩藻等。7月26日，發(fā)現(xiàn)該海域赤潮依然存在，估計(jì)面積超過200 km2，赤潮位置在朱家尖島東部一側(cè)（122°26′E以東海域），與前幾日相比，赤潮更接近朱家尖島東部沿岸。7月27日，發(fā)現(xiàn)在29°56′N 、122°38′E—29°51′N 、122°46′E，29°34′N 、122°21′E—29°27′N 、122°21′E沿線范圍內(nèi)均有赤潮，累計(jì)面積約400 km2，近海海域顏色為深褐色，離岸海域?yàn)榈稚?。赤潮海域已接近朱家尖島東部岸線，等步島、桃花島及普陀山海濱浴場等沿岸海域已受到赤潮影響。同日，對(duì)朱家尖東部海域的赤潮進(jìn)行應(yīng)急跟蹤監(jiān)測，發(fā)現(xiàn)大面積赤潮依然存在，面積超過200 km2。衛(wèi)星遙感反演得到的赤潮發(fā)生范圍見圖3，圖中紅色部分表示監(jiān)測到的赤潮區(qū)域。

圖4顯示了赤潮藻團(tuán)24 h和48 h漂移擴(kuò)散的位置，黑色實(shí)線表示衛(wèi)星遙感赤潮發(fā)生的位置，紅色虛線表示24 h預(yù)報(bào)位置，藍(lán)色虛線表示48 h預(yù)報(bào)位置。根據(jù)2007年7月24日衛(wèi)星遙感赤潮發(fā)生位置預(yù)測，24 h漂移方向?yàn)槠狈较?，漂移距離為9.61 km，漂移速度0.40 km/h，48 h漂移距離為8.19 km，漂移速度0.34 km/h。與7月25日衛(wèi)星遙感圖對(duì)比，由于計(jì)算沒有考慮赤潮藻種的生物化學(xué)變化，預(yù)報(bào)結(jié)果無法完全刻畫藻類生長繁殖導(dǎo)致赤潮藻團(tuán)的面積增加，預(yù)報(bào)落區(qū)與觀察結(jié)果還存在一定差異，這有可能是赤潮藻在海域內(nèi)增殖所致，但總體上預(yù)報(bào)赤潮藻團(tuán)的主要漂移擴(kuò)散方向與衛(wèi)星圖片反映的實(shí)際情況是一致的。從2007年7月25—27日的衛(wèi)星遙感圖可知，如短時(shí)間內(nèi)忽略赤潮藻類繁殖生長，赤潮的漂移擴(kuò)散趨勢(shì)是逐漸向北方向，與各時(shí)刻起報(bào)的預(yù)報(bào)結(jié)果一致，且24 h與48 h預(yù)測的漂移距離和速度相對(duì)較小，不利于赤潮藻團(tuán)的消散，故赤潮面積有進(jìn)一步發(fā)展趨勢(shì)。據(jù)東海海洋環(huán)境保護(hù)通報(bào)，7月28日，跟蹤監(jiān)測發(fā)現(xiàn)朱家尖東部海域赤潮依然存在，顏色為紅褐色、深褐色，赤潮生物優(yōu)勢(shì)種仍為扁面角毛藻。7月29日的跟蹤監(jiān)視顯示，舟山朱家尖-普陀山東部海域以及桃花島東南部海域的赤潮依然存在，赤潮區(qū)主要位于29°59′N 、122°24′E與29°51′N 、122°25′E兩點(diǎn)連線以東海域，觀測所及的赤潮海域面積超過200 km2，顏色為紅色、褐色。另外，監(jiān)測飛機(jī)在朱家尖東北方向發(fā)現(xiàn)大面積赤潮，赤潮主要位于距離朱家尖20 km以內(nèi)的海域，其中近岸海域附近顏色為褐色，較遠(yuǎn)海域?yàn)榈稚?。由此可見，從不同時(shí)刻遙感圖選定赤潮范圍，進(jìn)行24 h和48 h變化趨勢(shì)預(yù)報(bào)是可行的。

5 總結(jié)

在結(jié)合實(shí)例進(jìn)行準(zhǔn)業(yè)務(wù)化應(yīng)用預(yù)報(bào)試驗(yàn)中，能在較短時(shí)間內(nèi)實(shí)現(xiàn)赤潮藻團(tuán)的漂移擴(kuò)散預(yù)報(bào)，提供赤潮藻團(tuán)24—72 h的主要漂移擴(kuò)散方向、距離和速度等信息，具備數(shù)據(jù)預(yù)處理、數(shù)值模式計(jì)算及預(yù)報(bào)結(jié)果可視化的功能，預(yù)測結(jié)果與實(shí)際情況吻合。

圖4 2007年7月24日起預(yù)報(bào)赤潮藻團(tuán)24 h和48 h的漂移擴(kuò)散位置

圖5 2007年7月25日起預(yù)報(bào)赤潮藻團(tuán)24 h和48 h的漂移擴(kuò)散位置

圖6 2007年7月26日起預(yù)報(bào)赤潮藻團(tuán)24 h和48 h的漂移擴(kuò)散位置

模塊基于C/S架構(gòu)，具有響應(yīng)速度快的優(yōu)點(diǎn)，在客戶端的設(shè)計(jì)上，可以根據(jù)客戶的需求進(jìn)行個(gè)性化的設(shè)計(jì)，滿足客戶復(fù)雜的個(gè)性化需求。同時(shí)，C/S架構(gòu)充分利用服務(wù)器端硬件資源，將大量的數(shù)據(jù)處理工作在服務(wù)器端完成，，數(shù)據(jù)處理能力比較強(qiáng)大，在很大程度上提高了赤潮災(zāi)害的應(yīng)急處置能力。目前該模塊還在不斷完善中，數(shù)值模擬將赤潮的發(fā)生發(fā)展過程與環(huán)境條件進(jìn)行耦合計(jì)算，最終形成赤潮漂移擴(kuò)散預(yù)測預(yù)警系統(tǒng)，該系統(tǒng)的業(yè)務(wù)化應(yīng)用將進(jìn)一步完善現(xiàn)有條件下的赤潮預(yù)測系統(tǒng)，為赤潮防災(zāi)減災(zāi)及災(zāi)害的損失評(píng)估提供參考和依據(jù)。

[1]國家海洋局.2010年中國海洋災(zāi)害公報(bào)[R].2011.

[2]國家海洋局.2000年中國海洋環(huán)境質(zhì)量公報(bào)[R].2001.

[3]周為峰，樊偉.應(yīng)用MODIS進(jìn)行赤潮遙感監(jiān)測的研究進(jìn)展[J].遙感技術(shù)與應(yīng)用，2007，22(6):768-772.

[4]齊雨藻.中國沿海赤潮[M].北京:科學(xué)出版社，2003.

[5]唐峰華,伍玉梅,樊偉,等.長江口浮游植物分布情況及與徑流關(guān)系的初步探討[J].生態(tài)環(huán)境學(xué)報(bào),2010,19（12）：2934-2940.

[6]范麗,程金平,郭光勇.我國東海海域赤潮發(fā)生年際變化趨勢(shì)及其影響因素分析[J].上海環(huán)境科學(xué),2009,28(1):15-17.

[7]申力,徐惠平,吳平.長江口及東海赤潮海洋環(huán)境特征綜合研討[J].海洋環(huán)境科學(xué),2010,29(5):631-635.