何恩業(yè)，李雪丁，楊靜，高姍，季軒梁

（1.國家海洋環(huán)境預(yù)報中心自然資源部海洋災(zāi)害預(yù)報技術(shù)重點實驗室，北京100081；2.福建省海洋預(yù)報臺，福建福州350003）

1 引言

近年來，以赤潮和綠潮為主要災(zāi)種的海洋生態(tài)災(zāi)害頻發(fā)，成為與風(fēng)暴潮災(zāi)害、海浪災(zāi)害和海冰災(zāi)害并列的幾大海洋災(zāi)害之一。赤潮作為海洋生態(tài)災(zāi)害的主要災(zāi)種，其高頻次的發(fā)生給我國沿海地區(qū)造成了巨大的經(jīng)濟損失，對海洋生態(tài)系統(tǒng)平衡構(gòu)成嚴(yán)重威脅，有毒赤潮的發(fā)生還會給水產(chǎn)品食用安全保障帶來一定的壓力[1-2]。若能對赤潮進(jìn)行有效預(yù)測，為當(dāng)?shù)卣蜐O民提前預(yù)警，并有針對性地提前部署防治措施和處置相關(guān)海域養(yǎng)殖的水產(chǎn)品，可以降低甚至消除赤潮帶來的損失。

當(dāng)前赤潮的預(yù)測預(yù)報還是一個國際難題，總體水平偏低，也未形成統(tǒng)一的預(yù)報檢驗評價標(biāo)準(zhǔn)[2-3]。誘發(fā)赤潮的成因較多，既有生物因素，也涉及水動力、氣象和環(huán)境理化要素等因子，很難通過傳統(tǒng)的統(tǒng)計方法或數(shù)值方法對赤潮過程進(jìn)行準(zhǔn)確模擬和預(yù)測。國內(nèi)學(xué)者從多角度開展了赤潮預(yù)報研究工作[4-9]，但大都是對赤潮個例開展分析或者基于多要素進(jìn)行建模，這需要以大量和連續(xù)的觀測數(shù)據(jù)為基礎(chǔ)，在實際業(yè)務(wù)預(yù)報應(yīng)用中難免遇到監(jiān)測數(shù)據(jù)缺失問題。部分學(xué)者從氣象條件引發(fā)赤潮的角度入手進(jìn)行研究，但多以經(jīng)驗預(yù)報為主，預(yù)報員個人的經(jīng)驗、技巧和判斷對預(yù)報結(jié)果影響較大，主觀性較強，在實際赤潮預(yù)報業(yè)務(wù)中存在較大的局限性。以上問題的存在使得多數(shù)赤潮預(yù)報研究方法難以應(yīng)用于實際預(yù)報業(yè)務(wù)，不利于對赤潮災(zāi)害開展及時準(zhǔn)確的預(yù)測。本文從沿海赤潮業(yè)務(wù)預(yù)報的現(xiàn)實需求出發(fā)，利用穩(wěn)定易獲取的實時氣象數(shù)據(jù)，基于天氣分型的預(yù)報方法對赤潮進(jìn)行建模預(yù)測。以福建沿海為示范區(qū)，利用2001—2012年福建沿海赤潮記錄和美國國家環(huán)境預(yù)報中心（National Centers for Environmental Prediction，NCEP）的再分析數(shù)據(jù)構(gòu)建預(yù)報模型。敏感性試驗和預(yù)報檢驗結(jié)果表明，該模型具有較高的預(yù)報準(zhǔn)確率，可以在實踐中予以應(yīng)用。

2 數(shù)據(jù)來源和研究方法

2.1 數(shù)據(jù)來源

赤潮數(shù)據(jù)來源于《沿海海洋赤潮專報》和《中國海洋災(zāi)害公報》中2001—2012年發(fā)生在福建沿海的赤潮災(zāi)害過程記錄。數(shù)據(jù)包含歷次赤潮的起止時間、發(fā)生海域、成災(zāi)面積（km2）、優(yōu)勢生物種及密度和災(zāi)害損失情況等相關(guān)信息。本文以2001—2010年福建沿海赤潮過程作為建模所需的敏感性試驗樣本，以2011—2012年福建沿海赤潮過程作為預(yù)報檢驗樣本。

氣象數(shù)據(jù)來源于1981—2012年NCEP再分析數(shù)據(jù)資料，時間分辨率為4次/d，水平分辨率為2.5°×2.5°的全球格點數(shù)據(jù)。

2.2 研究方法

基于天氣分型的預(yù)報研究方法多集中于氣象方面[10-14]，赤潮方面較少。任湘湘等[15]曾嘗試將天氣分型思想引入赤潮預(yù)報，通過經(jīng)驗分析方法總結(jié)了5種易引發(fā)珠江口赤潮的天氣形勢，并對其中的1種用相似法進(jìn)行了檢驗。本文在前人工作的基礎(chǔ)上，對基于天氣分型的赤潮預(yù)報方法進(jìn)行優(yōu)化，采用全定量化的客觀方案構(gòu)建模型對容易誘發(fā)赤潮的天氣形勢進(jìn)行分型，利用相似算法對赤潮進(jìn)行預(yù)報。

基于天氣分型的赤潮預(yù)報模型分為3層（見圖1）。第一層為天氣形勢場數(shù)據(jù)信息輸入層，第二層為天氣分型運算層（包括天氣場聚類分析模塊、合成場試驗?zāi)K和典型天氣形勢分型模塊），第三層為結(jié)果輸出層（模型對輸入的天氣場進(jìn)行判別分類）。

圖1 基于天氣分型的赤潮預(yù)報模型結(jié)構(gòu)

聚類方法采用Pearson（皮爾森）相似算法，對樣本集合歷次赤潮首發(fā)日的天氣形勢場進(jìn)行歸類。合成天氣場采用互為極強相關(guān)且為歷次赤潮暴發(fā)首日的天氣場作為合成樣本集。

本研究利用Fortran、Grads、SPSS和Excel等軟件處理分析數(shù)據(jù)及繪圖，Pearson相似算法算式為：

式中：x和y分別為天氣形勢場Ai和Bi的距平場（氣候場基于1981—2010年的NCEP氣象再分析數(shù)據(jù)）；n表示天氣形勢場區(qū)域內(nèi)的格點數(shù)。

3 結(jié)果與分析

3.1 福建沿海赤潮時空分布特征

2001—2012年福建沿海年均赤潮約15.6次，年均累計面積約1 000 km2（見圖2a）。其中，2001年、2009年和2011年赤潮暴發(fā)頻次和規(guī)模較輕，赤潮年發(fā)頻次均未超過6次，年累計面積均小于300 km2；2002年、2003年、2006年和2010年赤潮暴發(fā)頻次和規(guī)模較重，赤潮年發(fā)頻次均高于20次，年累計面積均超過1 500 km2。

福建毗鄰北回歸線，屬亞熱帶海洋性季風(fēng)氣候，各屬地自然特征差異明顯，受海洋、大氣環(huán)流和地形等綜合影響，福建沿海赤潮也呈現(xiàn)出顯著的季節(jié)變化和區(qū)域性分布特征。每年4月福建沿海進(jìn)入赤潮多發(fā)季節(jié)，5—6月達(dá)到高峰時段，7月開始減弱，8月后發(fā)生頻次大幅降低，赤潮進(jìn)入低發(fā)期（見圖2b）。2001—2012年，福建沿海共發(fā)生赤潮187次，累計面積約11 700 km2。其中，4月第6候—7月第1候是赤潮暴發(fā)最盛期，在此期間福建沿海共發(fā)生赤潮153起，累計面積10 386 km2，分別占12a總次數(shù)和總面積的81.8%和88.5%。

圖3a為福建沿海10個預(yù)警海區(qū)劃分。通過對赤潮發(fā)生地點的分析發(fā)現(xiàn)，12a間廈門近海（Ⅷ區(qū)）發(fā)生赤潮次數(shù)最多，為55次，占總次數(shù)的29.4%，年均赤潮4.6次；其次為三沙灣海域（Ⅰ區(qū)和Ⅱ區(qū)），為68次，占比36.4%；平潭海域（Ⅳ區(qū)）27次，占比14.4%；閩江口海域（Ⅲ區(qū)）22次，占比11.8%；Ⅴ區(qū)、Ⅵ、Ⅶ區(qū)、Ⅸ區(qū)和Ⅹ區(qū)赤潮發(fā)生概率較小，合計發(fā)生15次，占比只有8%（見圖3b）。超過100 km2的大面積赤潮共有30次，主要集中在三沙灣海域（14次）、廈門近海（5次）和平潭海域（4次）。

圖2 2001—2012年福建沿海赤潮發(fā)生狀況統(tǒng)計

圖3 福建沿海預(yù)警海區(qū)劃分示意圖和各區(qū)2001—2012年赤潮發(fā)生狀況

3.2 典型天氣形勢場的合成

3.2.1 天氣分型模型參數(shù)初始化

天氣分型模型參數(shù)初始化是對天氣形勢場的時空范圍進(jìn)行界定，包含空間尺度和時間尺度參數(shù)，前者是對天氣場的垂直層次和水平尺度的范圍界定，后者是對天氣分型定義的時間跨度區(qū)間的界定。

眾多學(xué)者在大氣環(huán)流對赤潮的影響方面進(jìn)行了深入研究，不同程度地揭示了高空大氣環(huán)流對誘發(fā)中國沿海赤潮起到的關(guān)鍵性作用[16-19]。本文以高空500 hPa位勢高度場作為天氣分型的垂直尺度參數(shù)。水平尺度參數(shù)的確定考慮到誘發(fā)赤潮需要一個穩(wěn)定的天氣形勢，范圍選的太小不能把握天氣尺度系統(tǒng)的持續(xù)穩(wěn)定性，范圍太大又不具備反映區(qū)域性天氣形勢的典型特征。本文以水平尺度1 000～3 000 km、時間尺度1～3 d的天氣系統(tǒng)作為參考系，水平尺度參數(shù)設(shè)定在15°～35°N，90°～135°E。

依據(jù)福建沿海赤潮高發(fā)時段，本文將建模的時間尺度參數(shù)設(shè)定在4月26日—7月5日。此外，考慮到西太平洋副熱帶高壓（副高）的跳變規(guī)律造成我國天氣系統(tǒng)在6月前后存在時空尺度和水汽輸送上的季節(jié)性差異[19]，本文將模型的時間尺度參數(shù)細(xì)分為period T1（4月26日—5月31日）和period T2（6月1日—7月5日）兩個時段分別作為季節(jié)尺度參數(shù)。兩個時段的季節(jié)場具有不同的特征（見圖4）：period T1時段副高位置相對偏東和偏南，孟加拉灣上空的副熱帶南支西風(fēng)槽結(jié)構(gòu)特征比較穩(wěn)定，福建沿海氣候主要受來自孟加拉灣和南海北部的西-西南向氣流影響；period T2時段副高明顯西伸并北抬，福建沿海主體受副高外圍控制，氣候主要受來自南海的偏南向氣流影響。

圖4 兩個時段500 hPa季節(jié)場

3.2.2 數(shù)據(jù)處理

小面積赤潮受水質(zhì)、營養(yǎng)鹽和水動力等局地因素影響較大[3,20-22]，若將其加入到模型中進(jìn)行統(tǒng)計運算，會對天氣分型結(jié)果造成不利影響；其次，同一日期在不同海域也會發(fā)生多起赤潮，若重復(fù)運算會增加其所占權(quán)重；再者，同一海區(qū)接連發(fā)生赤潮，誘發(fā)這些赤潮的天氣形勢具有一貫性，因此在統(tǒng)計運算時也需要考慮。以上客觀因素的存在若不作處理，會影響模型的執(zhí)行效率及效果。為此，本文采用以下方案篩選敏感性試驗樣本：

（1）成災(zāi)面積超過20 km2的赤潮事件計入樣本集合[20]；

（2）成災(zāi)面積記錄空白的赤潮事件作剔除處理；

（3）同期發(fā)生在不同海域的赤潮事件，因其對應(yīng)的天氣形勢場相同，故取一處理；

（4）同區(qū)海域前后緊連且依次發(fā)生的赤潮，取最早暴發(fā)的赤潮事件作為樣本。

處理后，period T1時段共擇取敏感性試驗樣本36個，period T2時段共擇取敏感性試驗樣本26個。對各赤潮樣本首發(fā)日的500 hPa位勢高度距平場進(jìn)行交互式相似性運算，以其相關(guān)程度作為劃分依據(jù)，將天氣場進(jìn)行聚類：

（1）將互為強相關(guān)（ρ≥0.6）且通過P=0.05顯著性水平檢驗的天氣場歸為一類；

（2）對（1）形成的各類別中任一天氣場，篩選與其相關(guān)度ρ≥0.5的天氣場劃歸此類。

結(jié)果顯示，period T1時段天氣場形成兩類，共計樣本22個，占該時間段樣本集合的60.2%；period T2時段形成兩類，共計樣本19個，占該時間段樣本集合的73.1%（見表1）。其他剩余樣本之間并無顯著性的相關(guān)關(guān)系，不做歸類處理。

表1各赤潮首發(fā)日天氣場的分類

3.2.3 典型天氣形勢場的合成

將各類樣本集合中互為強相關(guān)的天氣場樣本進(jìn)行合成處理，以期獲取誘發(fā)赤潮的天氣形勢共性。初始合成場s1—s4見表2。將生成的4個合成場分別與2001—2010年相應(yīng)時段（period T1和period T2）的逐日天氣形勢場再次進(jìn)行相似運算，若結(jié)果呈現(xiàn)極強相關(guān)（ρ≥0.8）且當(dāng)日確為赤潮事件首發(fā)日時間，則該天氣形勢計入該類別并作為新一次合成場的樣本。如此反復(fù)試驗，當(dāng)合成場的樣本集合不再發(fā)生變化時，則最終的合成場認(rèn)定為最易誘發(fā)赤潮的典型天氣形勢場。

模型結(jié)果顯示（見表2和圖5）：典型場Type1由13個天氣形勢場樣本合成，為西南季風(fēng)擾動型，其主要特征是福建沿海受副高西北部外圍影響，海區(qū)等壓線呈西南-東北向分布，近海以暖濕的海洋性西南氣流為主；典型場Type2由5個樣本合成，為高壓脊控制型，主要特征是陸上高壓脊與副高打通合并控制了我國大部分地區(qū)，福建沿海處于脊線前部，大部海域受來自南中國海的氣流影響；典型場Type3由6個樣本合成，為東亞低槽擾動型，主要特征是副高位置相對偏南和偏東，福建沿海處于東亞大槽后部，影響該海區(qū)的氣流主要來自孟加拉灣海域；典型場Type4由有5個樣本合成，為副高控制型，主要特征是位于孟加拉灣上空的副熱帶南支西風(fēng)槽已經(jīng)消失，福建沿海完全處于強大的副高控制之下，影響近海的氣流主要來自西北太平洋海域。

表2 各類初始合成場和典型天氣形勢場樣本集合

圖5 誘發(fā)福建沿海赤潮的4種典型天氣形勢場

3.3 模型敏感性試驗和預(yù)報應(yīng)用與檢驗

3.3.1 模型敏感性試驗

2001—2010年赤潮天氣分型敏感試驗結(jié)果表明：若天氣形勢場與典型場相關(guān)系數(shù)連續(xù)2 d超過0.5，則發(fā)生赤潮的概率較大。為了對試驗結(jié)果進(jìn)行定量化評價，本研究參考匡曉迪等[23]在中尺度渦預(yù)報檢驗方面提出的評分方案，結(jié)合赤潮預(yù)報業(yè)務(wù)的客觀實際，引進(jìn)擬合準(zhǔn)確率和預(yù)報準(zhǔn)確率概念，設(shè)臨界參數(shù)值為0.5，定義各型典型場擬合準(zhǔn)確率為：

各型典型場擬合準(zhǔn)確率=(準(zhǔn)確擬合次數(shù)+

式中：準(zhǔn)確擬合是指連續(xù)2 d相關(guān)系數(shù)超過臨界值，且該時段包含赤潮個例的首發(fā)日；半對擬合是指相關(guān)系數(shù)連續(xù)超過臨界值的時間區(qū)間不含赤潮個例首發(fā)日，但與區(qū)間緊挨的下一日是赤潮事件首發(fā)日（考慮到赤潮監(jiān)測受客觀條件限制，赤潮實際暴發(fā)日可能早于發(fā)現(xiàn)日，即模型預(yù)報結(jié)果與實際吻合概率較大，故以半對計入統(tǒng)計）；空報是指相關(guān)系數(shù)連續(xù)超過臨界值的時間段不包含赤潮個例首發(fā)日，且排除半對擬合情況。此外，在敏感性試驗中，由于不同分型的典型場只能預(yù)報本類別天氣形勢誘發(fā)的赤潮，對于其他不同類型天氣形勢誘發(fā)的赤潮不具有預(yù)報能力，因此本研究對各單一典型場敏感試驗結(jié)果驗證不考慮漏報情況。評價結(jié)果顯示（見表3），典型場的擬合準(zhǔn)確率為Type1 61.1%、Type2 65.0%、Type3 60.7%和Type4 64.7%。

3.3.2 模型預(yù)報應(yīng)用與檢驗

2011年福建沿海發(fā)生赤潮5次，是12a中次數(shù)最少的一年，其中4次發(fā)生在period T1時段，periodT2時段未有赤潮發(fā)生；2012年福建沿海發(fā)生赤潮17次，其中12次發(fā)生在period T1時段，4次發(fā)生在period T2時段（見表4）?；谀Ｐ蜕傻?種典型場對2011—2012年福建沿海赤潮進(jìn)行預(yù)報檢驗，定義模型綜合預(yù)報準(zhǔn)確率為：

表3敏感性試驗結(jié)果

檢驗結(jié)果顯示（見表4），2011—2012年20次赤潮過程中，該模型成功預(yù)報了其中的12次，半對預(yù)報了5次，空報了6次，漏報了3次，綜合預(yù)報準(zhǔn)確率為55.8%。分時段來看，period T1時段（Type1和Type2）準(zhǔn)確預(yù)報了9次，半對預(yù)報了4次，漏報了3次，空報了3次，綜合預(yù)報準(zhǔn)確率為57.9%；period T2時段（Type3和Type4）準(zhǔn)確預(yù)報了3次，半對預(yù)報了1次，無漏報次數(shù)，空報了3次，綜合預(yù)報準(zhǔn)確率為50.0%。結(jié)果還顯示，典型場Type2準(zhǔn)確預(yù)報了8次赤潮，關(guān)鍵原因是2012年5月24—26日福建沿海連續(xù)發(fā)生了4次赤潮，2012年5月30日當(dāng)天發(fā)生了4次赤潮，這8次赤潮首發(fā)日的天氣形勢都與典型場Type2呈現(xiàn)出較強的相關(guān)性。此外，在2011—2012年period T2時段與Type3和Type4呈現(xiàn)出較高相似度的天氣場次數(shù)明顯偏少，合計只有7次，而實際上period T2時段赤潮發(fā)生次數(shù)也明顯偏少，合計只有4次。模型成功預(yù)報了其中的3次，剩余1次也為半對預(yù)報，顯示出赤潮暴發(fā)頻次的多寡受天氣形勢影響因素較大。該結(jié)果在一定程度上說明天氣形勢在誘發(fā)赤潮過程中起關(guān)鍵性作用，基于天氣分型的赤潮預(yù)測模型能較準(zhǔn)確地抓住誘發(fā)赤潮的天氣形勢主要特征，并且可以通過定量化的方法對赤潮的發(fā)生進(jìn)行預(yù)測。

表4基于天氣分型的赤潮預(yù)報模型后報檢驗結(jié)果

3.4 誘發(fā)福建沿海赤潮的典型天氣形勢個例分析

3.4.1典型場Type1：西南季風(fēng)擾動型

由該型天氣形勢誘發(fā)的福建沿海赤潮主要集中在4月底—5月底。2001—2012年，period T1時段共發(fā)生赤潮92次，其中44次與典型場Type1具有強相關(guān)性，占比達(dá)47.8%，平均每年約有3次赤潮過程與該類天氣形勢誘發(fā)有關(guān)。例如2002年5月4—14日發(fā)生在福鼎市和霞浦縣東部海域面積為500 km2的赤潮。

2002年5月2—4日高空500 hPa位勢高度場與典型場Type1的相關(guān)系數(shù)依次為0.90、0.92和0.94，連續(xù)多日達(dá)到極強相關(guān)程度，說明本次天氣過程與典型場Type1高度相似，且極其穩(wěn)定（見圖6）。500 hPa位勢高度場上，西太副高588線壓在了福建沿海，且588線包絡(luò)了臺灣省、巴士海峽以及菲律賓大部分地區(qū)，與此同時，位于孟加拉灣上空的副熱帶南支西風(fēng)槽持續(xù)發(fā)展加深，受此影響，來自西南季風(fēng)的暖濕氣流源源不斷地輸送到我國東南部區(qū)域；近地面層，5月2—3日福建沿海處于低壓倒槽之前，海區(qū)等壓線連續(xù)多日呈西北-東南向，同時氣壓梯度力小，這種天氣形勢導(dǎo)致海區(qū)以高溫和高濕的東南風(fēng)為主，且風(fēng)力較小，水體較穩(wěn)定，為赤潮的醞釀和暴發(fā)提供了有利的氣象條件。

3.4.2 典型場Type2：高壓脊控制型

由該型天氣形勢誘發(fā)的赤潮主要發(fā)生在5月中旬和下旬。period T1時段有18次與典型場Type2具有強相關(guān)性，占比19.6%，平均每年約有1～2次赤潮過程與該類天氣形勢誘發(fā)有關(guān)，例如2007年5月29日—6月1日發(fā)生在同安灣海域面積為20 km2的旋鏈角毛藻赤潮。

圖6 2002年5月2—4日500 hPa位勢高度場和近地面層氣壓場

2007年5月27—29日天氣形勢場與典型場Type2的相關(guān)系數(shù)依次為0.63、0.64和0.93，相關(guān)程度由強相關(guān)遞增到極強相關(guān)（見圖7）。5月27—28日，我國陸上西部高壓脊開始發(fā)展東移，并逐漸影響我國華南和東南沿海，伴隨西太副高的西伸北進(jìn)，5月29日陸上高壓脊與西太副高打通合并形成強大的阻塞高壓，造成福建沿海地區(qū)出現(xiàn)極其穩(wěn)定的晴好天氣；在近地面層，持續(xù)多日的縱向等壓線分布態(tài)勢，迫使來自南中國海的暖濕氣流持續(xù)不斷地輸送至福建沿海，良好的光照加之能量的持續(xù)積累，最終誘發(fā)了福建沿海水華。

3.4.3典型場Type3：東亞低槽擾動型

由該型天氣形勢誘發(fā)的福建沿海赤潮主要發(fā)生在6月上旬和中旬。2001—2012年，period T2時段共發(fā)生赤潮61次，其中28次與典型場Type3具有強相關(guān)性，占比達(dá)到45.9%，平均每年約有2次赤潮過程與該類天氣形勢誘發(fā)有關(guān)。典型的赤潮事件如2006年6月11—14日發(fā)生在平潭赤潮監(jiān)控區(qū)及周邊海域的中肋骨條藻赤潮，成災(zāi)面積600 km2；2006年6月12—13日，分別發(fā)生在深滬灣海域（成災(zāi)面積100 km2）、崇武半島海域（成災(zāi)面積200 km2）和南日群島海域（成災(zāi)面積150 km2）的赤潮。

2006年6月8—10日天氣形勢場與典型場Type3連續(xù)多日達(dá)到強相關(guān)程度，其相關(guān)系數(shù)依次為0.79、0.68和0.67（見圖8）。在500 hPa位勢高度場上東亞大槽不斷加深，福建沿海處于低壓槽末端影響之下，西太副高588線相對偏南和偏東，大體處于巴士海峽的東南部地區(qū)，受此影響福建沿海上升流比較活躍。對應(yīng)近地面層，6月8—9日溫帶氣旋東移過境，致使福建沿海地區(qū)出現(xiàn)氣壓降低和氣流輻合的天氣現(xiàn)象。該天氣形勢的配置易在近海海域產(chǎn)生上升流，利于水體底層營養(yǎng)鹽向表層輸送，為赤潮生物的繁殖生長創(chuàng)造了有利的營養(yǎng)條件。6月10日后北部弱高壓南下跟進(jìn)控制了福建沿海，大氣運動處于相對穩(wěn)定狀態(tài)，為觸發(fā)水華提供了有利的氣象條件，之后6月11—12日福建沿海發(fā)生多起赤潮。

3.4.4典型場Type4：西太平洋副熱帶高壓控制型

由該型天氣形勢誘發(fā)的赤潮主要發(fā)生在6月下旬和7月初。period T2時段有17次與典型場Type4具有強相關(guān)性，占比27.9%，平均每年約有1～2次赤潮過程與該類天氣形勢誘發(fā)有關(guān)。以2006年6月21—24日發(fā)生在東山灣海域面積為200 km2的日本星桿藻和旋鏈角毛藻赤潮為例，對誘發(fā)這次赤潮的天氣過程進(jìn)行分析。

圖7 2007年5月27—29日500 hPa位勢高度場和近地面層氣壓場

圖8 2006年6月8—10日500 hPa位勢高度場和近地面層氣壓場

2006年6月19—21日的天氣形勢場與典型場Type4連續(xù)多日達(dá)到極強相關(guān)程度，其相關(guān)系數(shù)依次為0.88、0.89和0.88（見圖9）。500 hPa位勢高度場上，西太副高588線西伸并北抬且持續(xù)控制浙江、福建和廣東東部近海，造成沿海海域下沉氣流運動較強，天氣形勢穩(wěn)定，形成晴熱高溫天氣；對應(yīng)近地面層，福建沿海受海上高壓長時間控制，暖濕的東南風(fēng)持續(xù)，天氣悶熱，海溫逐漸升高，這種天氣形勢容易誘發(fā)水體交換不暢的半封閉性海域赤潮。

圖9 2006年6月19—21日500 hPa位勢高度場和近地面層氣壓場

4 結(jié)論

本文以2001—2012年我國沿海赤潮的相關(guān)數(shù)據(jù)，對福建沿海赤潮的時空分布特征進(jìn)行分析。在此基礎(chǔ)上，利用NCEP再分析資料，從大氣環(huán)流誘發(fā)赤潮的角度出發(fā)，基于典型天氣形勢場相似法構(gòu)建了福建沿海赤潮定量化預(yù)報模型，并對該模型進(jìn)行了敏感性試驗和預(yù)報檢驗。結(jié)果表明：

（1）2001—2012年，福建沿海共發(fā)生赤潮187次，累計面積約11 700 km2，年均赤潮約15.6次，年均面積約1 000 km2。其中，4月第6候—7月第1候是福建沿海赤潮暴發(fā)的鼎盛期，期間共發(fā)生赤潮153起，累計面積10 386 km2，占12a總次數(shù)的81.8%，占總面積的88.5%。從區(qū)域分布來看，閩東北海域和廈門近海是赤潮的高發(fā)海域，大面積赤潮主要集中在三沙灣海域、廈門近海和平潭海域。

（2）敏感性試驗和預(yù)報檢驗結(jié)果表明：基于天氣分型的福建沿海赤潮預(yù)報模型能夠在赤潮大面預(yù)報中獲得較好的效果，其后報檢驗準(zhǔn)確率為55.8%。其中，period T1時段（4月25日—5月31日）的綜合預(yù)報準(zhǔn)確率為57.9%；period T2時段（6月1日—7月5日）的綜合預(yù)報準(zhǔn)確率為50.0%。

（3）誘發(fā)福建沿海赤潮的主要天氣形勢分別為西南季風(fēng)擾動型、高壓脊控制型、東亞低槽擾動型和西太平洋副熱帶高壓控制型。4月26日—7月5日期間，在福建沿海153次赤潮過程中，由該4種天氣形勢誘發(fā)的赤潮為107次，占比達(dá)69.9%，說明該模型抓住了誘發(fā)福建沿海赤潮環(huán)流形勢的主要特征，具有較好的預(yù)報效果。

赤潮成因復(fù)雜，既有赤潮生物種自身的生物因素，也涉及到水動力、氣象和環(huán)境理化要素等因子。本文只從氣象條件對赤潮誘發(fā)的角度進(jìn)行建模分析，在大面預(yù)報中效果較好，若在此基礎(chǔ)上加入實時生態(tài)觀測數(shù)據(jù)進(jìn)行綜合分析則可以做到更精細(xì)化的預(yù)報。總之，在當(dāng)前赤潮生態(tài)動力學(xué)模型尚不成熟的階段，基于天氣分型的定量化預(yù)報方法具有較強的可操作性，為赤潮短期業(yè)務(wù)預(yù)報提供了一種新的途徑。