曾銀東

（1.福建省海洋預(yù)報(bào)臺(tái)，福建福州350003；2.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東青島266071）

平潭海域精細(xì)化三維溫鹽流業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)

曾銀東1,2

（1.福建省海洋預(yù)報(bào)臺(tái)，福建福州350003；2.中國(guó)海洋大學(xué)海洋與大氣學(xué)院，山東青島266071）

根據(jù)平潭附近海域?qū)嶋H需求和地理區(qū)位特點(diǎn)，研發(fā)了三維溫、鹽、流高分辨率業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)，對(duì)該系統(tǒng)采用的模式、配置、預(yù)報(bào)結(jié)果驗(yàn)證和業(yè)務(wù)化運(yùn)行等情況進(jìn)行了詳細(xì)介紹。該系統(tǒng)基于ROMS海流模式，采用正交曲線網(wǎng)格和網(wǎng)格嵌套技術(shù)；建立了包含數(shù)據(jù)預(yù)處理、溫鹽流預(yù)報(bào)和預(yù)報(bào)產(chǎn)品后處理等子程序的業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)。結(jié)果表明：模式預(yù)報(bào)產(chǎn)品具有較高精度，整個(gè)業(yè)務(wù)流程運(yùn)行在30 min之內(nèi)，系統(tǒng)運(yùn)行穩(wěn)定，達(dá)到業(yè)務(wù)化預(yù)報(bào)的要求。

平潭海域；精細(xì)化；海流模式；業(yè)務(wù)化；數(shù)值預(yù)報(bào)系統(tǒng)

1 引言

近年來(lái)，隨著沿海各國(guó)大力發(fā)展海洋經(jīng)濟(jì)，如何保障海上活動(dòng)的公共安全成為關(guān)注焦點(diǎn)，其中海上落水人員搜救、海上溢油的保障預(yù)報(bào)技術(shù)成為維護(hù)國(guó)家海洋權(quán)益、保障海上公共安全的研究重點(diǎn)。為此，各國(guó)都在積極發(fā)展全球海洋預(yù)報(bào)系統(tǒng)、區(qū)域嵌套高分辨率海洋預(yù)報(bào)系統(tǒng)，以及針對(duì)海灣和河口等多重嵌套精細(xì)化海洋預(yù)報(bào)系統(tǒng)。例如，法國(guó)基于NEMO（Nucleus for European Models of the Ocean）模式建立了Mercator海洋預(yù)報(bào)系統(tǒng)[1]；意大利同樣基于NEMO模式建立了“地中海海上安全決策支持系統(tǒng)”（Mediterranean Decision Support System for Marine Safety，MEDESS4MS）[2]；美國(guó)建立了HYCOM（Hybrid Coordinate Ocean Mode）/NCODA（Navy Coupled Ocean Data Assimilation）系統(tǒng)和GRTOFS（Global Real-Time Ocean Forecast System）系統(tǒng)[3-5]；國(guó)家海洋環(huán)境預(yù)報(bào)中心基于ROMS（Regional Ocean Modeling System）模式建立了全球業(yè)務(wù)化海洋學(xué)預(yù)報(bào)系統(tǒng)和中國(guó)海三維溫鹽流數(shù)值預(yù)報(bào)系統(tǒng)[4-5]；國(guó)家海洋局北海預(yù)報(bào)中心基于ROMS模式建立了黃渤海近岸精細(xì)化三維海溫、鹽度、海流（以下簡(jiǎn)稱溫、鹽、流）業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)[6]。業(yè)務(wù)化海洋預(yù)報(bào)系統(tǒng)可以提供相對(duì)可靠的海洋環(huán)境預(yù)報(bào)產(chǎn)品，在海洋工程建設(shè)、海上交通運(yùn)輸、海洋權(quán)益維護(hù)、海上搜救和應(yīng)急處置等方面發(fā)揮了重要作用。

平潭位于福建東部，其近岸海域島嶼眾多，海底地形復(fù)雜，海況和水文氣象條件復(fù)雜多變，是臺(tái)灣海峽重要的水動(dòng)力研究區(qū)域[7-9]。夏季受西南季風(fēng)影響，該海域存在明顯上升流，而冬季由長(zhǎng)江、閩江和九龍江等入海徑流與沿岸海水混合形成的富含營(yíng)養(yǎng)鹽、低溫、低鹽的浙閩沿岸流到達(dá)平潭[10]，對(duì)臺(tái)灣海峽的環(huán)流結(jié)構(gòu)、水團(tuán)組成、海洋生產(chǎn)力、沿岸物質(zhì)輸運(yùn)等產(chǎn)生重要的影響。平潭近岸海域也因此成為船舶交通事故和溢油事故多發(fā)區(qū)，海上船只航行、停泊、漁業(yè)生產(chǎn)和海上作業(yè)安全事故時(shí)有發(fā)生。例如，2010年9月30日凌晨，天津籍貨船“惠盈168”輪在平潭島附近海域沉沒(méi)，船上15名船員落水，當(dāng)晚已有3人死亡、3人獲救，另外9人失蹤。目前，針對(duì)臺(tái)灣海峽及其周邊海域已建立了業(yè)務(wù)化三維溫鹽流數(shù)值預(yù)報(bào)系統(tǒng)和海上突發(fā)事故輔助決策系統(tǒng)[11]，但該系統(tǒng)在平潭近岸海域空間分辨率較低，無(wú)法為該海域提供精細(xì)化的三維溫、鹽、流數(shù)值預(yù)報(bào)產(chǎn)品和服務(wù)。因此，本研究將利用當(dāng)前該海域高精度的水深地形資料和大量的海洋觀測(cè)數(shù)據(jù)，基于國(guó)際先進(jìn)的ROMS海洋環(huán)流模式，采用正交曲線網(wǎng)格和網(wǎng)格嵌套技術(shù)，發(fā)展平潭海域精細(xì)化三維溫、鹽、流數(shù)值預(yù)報(bào)系統(tǒng)，并對(duì)模式進(jìn)行驗(yàn)證和業(yè)務(wù)化試運(yùn)行，最終實(shí)現(xiàn)該系統(tǒng)的業(yè)務(wù)化運(yùn)行，為該海域海上活動(dòng)公共安全提供精細(xì)化溫、鹽、流數(shù)值預(yù)報(bào)產(chǎn)品。

2 系統(tǒng)模式介紹與配置

2.1 ROMS模式介紹

ROMS是一個(gè)自由表面、靜水力學(xué)和三維非線性斜壓原始方程模式，由羅格斯大學(xué)和加利福尼亞大學(xué)洛杉磯分校共同開(kāi)發(fā)完成[12-13]，該模式可以模擬多種尺度的海洋環(huán)流運(yùn)動(dòng)，如全球尺度的環(huán)流模擬、中尺度由氣象因素或是天文潮所引起的水位與流場(chǎng)變化，也可以計(jì)算小尺度河流渠道等流體運(yùn)動(dòng)。模式垂直方向采用S坐標(biāo)系，垂向分層隨著地形變化；水平方向采用正交曲線坐標(biāo)系及Arakawa C差分網(wǎng)格；為提高計(jì)算效率，采用了模態(tài)分離法，將帶自由表面的三維流動(dòng)問(wèn)題分成表面波的傳播問(wèn)題（外模態(tài)）和內(nèi)波的傳播問(wèn)題（內(nèi)模態(tài)），對(duì)外模態(tài)采用顯式差分；計(jì)算垂向渦動(dòng)黏性系數(shù)和紊動(dòng)擴(kuò)散系數(shù)由2.5階紊流模式[14]確定。

2.2 模式配置

2.2.1 模式區(qū)域、網(wǎng)格和水深地形

本模式采用曲線正交網(wǎng)格和雙重嵌套方法，大網(wǎng)格區(qū)域范圍為95°～147.9°E，9°S～44.01°N，小網(wǎng)格區(qū)域?yàn)槠教督逗Ｓ颍ㄒ?jiàn)圖1）。大網(wǎng)格開(kāi)邊界網(wǎng)格距約為45 km，逐漸過(guò)渡到臺(tái)灣海峽為1.5 km。整個(gè)模式區(qū)域的水深地形數(shù)據(jù)采用多種數(shù)據(jù)融合后的數(shù)據(jù)：外海區(qū)域水深地形資料采用ETOPO2；近海區(qū)域水深地形數(shù)據(jù)采用數(shù)字化的沿岸海圖水深（共26張）和臺(tái)灣海峽地形數(shù)據(jù)（0.5′網(wǎng)格化水深）；福清灣海域采用高精度的1∶1000水深。經(jīng)融合和插值到數(shù)值模式所需要的網(wǎng)格化水深和岸線如圖1所示。大、小網(wǎng)格垂向均取30層，采用S坐標(biāo)進(jìn)行垂向分層，在海表分層較密，以反映海表混合層及溫、鹽躍層的復(fù)雜垂直結(jié)構(gòu)。小網(wǎng)格在大網(wǎng)格中以1∶7加密，平潭海域達(dá)到200 m的水平分辨率（見(jiàn)圖1）。對(duì)于平潭小區(qū)域，本模式通過(guò)定義臨界水深（取0.2 m）進(jìn)行干濕邊界的轉(zhuǎn)換，以準(zhǔn)確體現(xiàn)該區(qū)域的納潮量。

圖1 模式計(jì)算網(wǎng)格區(qū)域范圍、水深地形和觀測(cè)站點(diǎn)圖

2.2.2 模式上表面和開(kāi)邊界條件

模式在大小網(wǎng)格的海表都考慮了海表熱通量，并在模式中考慮了河流淡水流入。模式上邊界采用福建省海洋預(yù)報(bào)臺(tái)業(yè)務(wù)化運(yùn)行的西北太平洋的WRF大氣模式風(fēng)場(chǎng)和熱通量場(chǎng)等作為大氣強(qiáng)迫場(chǎng)。模式所需的河流流量數(shù)據(jù)采用《中國(guó)海灣志》的主要河流各月及多年平均流量（見(jiàn)表1）。

模式開(kāi)邊界控制參數(shù)含溫、鹽、流及水位，其中水位含與大洋環(huán)流相關(guān)的大洋表面高度及潮汐所引起的海表起伏。模式的開(kāi)邊界采用MyOcean的業(yè)務(wù)化產(chǎn)品，包括海表高度、海溫、鹽度和海流。MyOcean業(yè)務(wù)化產(chǎn)品由ECMWF（European Centre for Medium-range Weather Forecasts）研發(fā)海洋信息服務(wù)支持系統(tǒng)提供，MyOcean是基于觀測(cè)數(shù)據(jù)、數(shù)值模型與數(shù)據(jù)同化系統(tǒng)建立的全球海洋實(shí)時(shí)觀測(cè)和預(yù)報(bào)系統(tǒng)，可以提供全球范圍海洋觀測(cè)和預(yù)報(bào)產(chǎn)品服務(wù)。開(kāi)邊界的潮汐調(diào)和常數(shù)采用日本NAOTIDE（National Astronomical Observatory（Japan）Tide Prediction Model）模式提供的天文潮調(diào)和常數(shù)，包含M2、S2、K1、O1、N2、P1、K2、Q1、M1、J1、OO1、2N2、Mu2、Nu2、L2和T2 16個(gè)主要天文分潮。

模式上表層熱通量控制邊界條件：

式中：為海表熱及降水凈通量，數(shù)值取自采用福建省海洋預(yù)報(bào)臺(tái)業(yè)務(wù)化運(yùn)行的西北太平洋的WRF大氣模式產(chǎn)品獲取。

表1 中國(guó)主要河流多年平均徑流量的數(shù)據(jù)（單位：104m3/s）

3 模式驗(yàn)證

3.1 實(shí)測(cè)數(shù)據(jù)來(lái)源

海洋觀測(cè)資料來(lái)自福建省海洋預(yù)報(bào)臺(tái)布放在平潭牛山島海域的海床基和平潭近海的浮標(biāo)（位置見(jiàn)圖1），分別代表平潭沿岸海域和平潭近海。海床基布放于海底水深約50 m，從海底向海表觀測(cè)共25層剖面流；浮標(biāo)觀測(cè)表層海溫和海表至水下約50 m深處剖面流，每5 m一層，共10層。

3.2 驗(yàn)證與討論

圖2和圖3為模式24 h、48 h和72 h預(yù)報(bào)海流分別與海床基、浮標(biāo)觀測(cè)海流的對(duì)比情況。由圖可見(jiàn)，模式預(yù)報(bào)流速、流向與觀測(cè)值吻合較好。模式預(yù)報(bào)值和觀測(cè)值均反映了夏季7月份盛行西南風(fēng)的情況下，在潮汐及環(huán)流的共同作用下，流場(chǎng)以往復(fù)流為主，一天之內(nèi)有4個(gè)峰值出現(xiàn)，南向流速大于北向流速。表、中、底層的預(yù)報(bào)流速與觀測(cè)值比較表明，流速較為一致，表層流速大，中、底層流速較小，表、中、底層流向基本一致。在2014年7月23—25日之間出現(xiàn)強(qiáng)的連續(xù)南向流是由于“麥德姆”臺(tái)風(fēng)過(guò)程引起，模式預(yù)報(bào)和觀測(cè)值誤差高于其它時(shí)間段，但兩者總體表現(xiàn)比較一致，能夠反應(yīng)臺(tái)風(fēng)的影響過(guò)程。

圖2 牛山島模式不同時(shí)間預(yù)報(bào)海流與海床基觀測(cè)海流對(duì)比

圖2 （續(xù)）

圖3為模式24 h、48 h和72 h預(yù)報(bào)海溫與浮標(biāo)觀測(cè)海溫的對(duì)比情況。從圖上看，模式預(yù)報(bào)海溫與浮標(biāo)實(shí)測(cè)海溫變化趨勢(shì)比較一致。23—25日浮標(biāo)觀測(cè)海溫出現(xiàn)急劇下降，形成“V”字形狀，主要是受1410號(hào)臺(tái)風(fēng)“麥德姆”正面影響，臺(tái)風(fēng)影響過(guò)后海溫迅速回升，臺(tái)風(fēng)的抽吸作用是導(dǎo)致海溫下降的主要原因[15]。

為了進(jìn)一步檢驗(yàn)?zāi)Ｊ筋A(yù)報(bào)的準(zhǔn)確性，本文選用了比較常見(jiàn)的檢驗(yàn)統(tǒng)計(jì)參數(shù)絕對(duì)平均誤差（MAE）和均方根誤差（RMSE）等統(tǒng)計(jì)方法。

（1）絕對(duì)平均誤差：

（2）均方根誤差

對(duì)模式預(yù)報(bào)結(jié)果分別與海床基和浮標(biāo)實(shí)測(cè)數(shù)據(jù)進(jìn)行誤差統(tǒng)計(jì)，得到表2、表3和表4?？傮w而言，模式預(yù)報(bào)結(jié)果與實(shí)測(cè)結(jié)果比較一致，模式在平潭近海預(yù)報(bào)結(jié)果優(yōu)于沿岸海域。

圖3 3號(hào)浮標(biāo)模式不同時(shí)間預(yù)報(bào)海流與浮標(biāo)實(shí)測(cè)海流對(duì)比

圖3 （續(xù)）

圖4 模式24 h、48 h和72 h預(yù)報(bào)海溫與浮標(biāo)觀測(cè)海溫的對(duì)比

表2 模式預(yù)報(bào)流速（單位：m/s）和流向（單位：°）結(jié)果與海床基觀測(cè)數(shù)據(jù)的誤差統(tǒng)計(jì)

表3 模式預(yù)報(bào)流速（單位：m/s）和流向（單位：°）結(jié)果與3號(hào)浮標(biāo)觀測(cè)數(shù)據(jù)的誤差統(tǒng)計(jì)

表2結(jié)果顯示，模式預(yù)報(bào)流速絕對(duì)平均誤差小于0.15 m/s，表層最大，底層最小，預(yù)報(bào)時(shí)效延長(zhǎng)誤差差異不大；預(yù)報(bào)流向絕對(duì)平均誤差小于35.4°，隨著預(yù)報(bào)時(shí)效延長(zhǎng)誤差略有升高，表層最大，底層最??；預(yù)報(bào)流速均方根誤差小于0.20 m/s，表層最大，底層最小，隨著預(yù)報(bào)時(shí)效延長(zhǎng)誤差略有下降；預(yù)報(bào)流向均方根誤差小于43.4°，表層最大，底層最小，隨著預(yù)報(bào)時(shí)效延長(zhǎng)誤差呈下降趨勢(shì)。

表3結(jié)果顯示，模式預(yù)報(bào)流速絕對(duì)平均誤差都小于0.11 m/s，表層最大，底層最小，預(yù)報(bào)時(shí)效延長(zhǎng)誤差沒(méi)有顯著差異；預(yù)報(bào)流向絕對(duì)平均誤差小于27.6°，表層最大，底層最小，預(yù)報(bào)時(shí)效延長(zhǎng)誤差略有增大；預(yù)報(bào)流速均方根誤差小于0.15 m/s，表層最大，底層最小，隨著預(yù)報(bào)時(shí)效的延長(zhǎng)誤差下降；預(yù)報(bào)流向均方根誤差小于28.2°，表層最大，底層最小，隨著預(yù)報(bào)時(shí)效的延長(zhǎng)誤差略有增大。

表4結(jié)果顯示，模式預(yù)報(bào)24 h、48 h和72 h表層海溫與相應(yīng)的實(shí)測(cè)海溫絕對(duì)平均誤差約1℃左右，隨著預(yù)報(bào)時(shí)效延長(zhǎng)誤差略有增大；模式預(yù)報(bào)表層海溫與實(shí)測(cè)海溫均方根誤差約1.3℃左右，并且同樣隨著預(yù)報(bào)時(shí)效延長(zhǎng)誤差略有增大。

3 系統(tǒng)業(yè)務(wù)化流程及其應(yīng)用

平潭海域精細(xì)化三維海流業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)流程圖如圖5所示。該系統(tǒng)依托福建省海洋預(yù)報(bào)臺(tái)曙光40萬(wàn)億次高性能計(jì)算機(jī)，其存儲(chǔ)容量為248TB。該預(yù)報(bào)系統(tǒng)在預(yù)報(bào)結(jié)果與實(shí)測(cè)數(shù)據(jù)對(duì)比驗(yàn)證和業(yè)務(wù)化試運(yùn)行的基礎(chǔ)上，目前已經(jīng)實(shí)現(xiàn)了從資料接收、處理、模式運(yùn)行到產(chǎn)品可視化及發(fā)布的安全、高效、穩(wěn)定業(yè)務(wù)化運(yùn)行。該系統(tǒng)已穩(wěn)定業(yè)務(wù)化運(yùn)行2—3 a。系統(tǒng)起報(bào)時(shí)間為每日20時(shí)（北京時(shí)），提供未來(lái)72 h溫鹽流預(yù)報(bào)產(chǎn)品。系統(tǒng)業(yè)務(wù)具體流程為：

表4 模式預(yù)報(bào)表層海溫（單位：℃）結(jié)果與3號(hào)浮標(biāo)觀測(cè)數(shù)據(jù)的誤差統(tǒng)計(jì)

圖5 平潭海域精細(xì)化海流業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)流程圖

圖6 平潭海域精細(xì)化溫、鹽、流業(yè)務(wù)化數(shù)值預(yù)報(bào)產(chǎn)品展示圖

圖7 預(yù)報(bào)系統(tǒng)應(yīng)用界面

（1）每日9時(shí)高性能計(jì)算機(jī)啟動(dòng)預(yù)報(bào)主程序；

（2）數(shù)據(jù)預(yù)處理：系統(tǒng)自動(dòng)獲取福建省海洋預(yù)報(bào)臺(tái)業(yè)務(wù)化的WRF氣象模式預(yù)報(bào)產(chǎn)品，并代入模式計(jì)算得到系統(tǒng)所需的海表面凈熱通量、水氣通量和動(dòng)量通量作為海洋模式的海表驅(qū)動(dòng)；系統(tǒng)自動(dòng)下載Myocean當(dāng)天的業(yè)務(wù)化海表高度、溫、鹽和流預(yù)報(bào)產(chǎn)品作為該系統(tǒng)海洋模式所需的開(kāi)邊界數(shù)據(jù)；該系統(tǒng)海洋模式采用熱啟動(dòng)方式，并采用前兩天存貯的歷史文件作為該海洋模式的初始場(chǎng)；

（3）進(jìn)行72 h三維溫、鹽、流數(shù)值預(yù)報(bào)，運(yùn)行時(shí)間為30 min。模式預(yù)報(bào)結(jié)果產(chǎn)生后，系統(tǒng)自動(dòng)啟動(dòng)產(chǎn)品后處理子系統(tǒng)，并生成標(biāo)準(zhǔn)數(shù)據(jù)格式，自動(dòng)化制作圖形圖像產(chǎn)品（見(jiàn)圖6），并將逐時(shí)海流場(chǎng)數(shù)據(jù)自動(dòng)推送到福建省海上突發(fā)事故應(yīng)急輔助決策系統(tǒng)和赤潮災(zāi)害預(yù)警系統(tǒng)業(yè)務(wù)平臺(tái)（見(jiàn)圖7），開(kāi)展海上落水人漂移軌跡預(yù)測(cè)和赤潮漂移軌跡預(yù)測(cè)工作，為海上搜救應(yīng)急提供保障服務(wù)，為赤潮漂移對(duì)養(yǎng)殖區(qū)、濱海浴場(chǎng)等潛在風(fēng)險(xiǎn)的應(yīng)急處置決策提供技術(shù)支撐。整個(gè)預(yù)報(bào)系統(tǒng)無(wú)需人工干預(yù)，具有較高的自動(dòng)化程度。

4 結(jié)論

本文基于ROMS海洋模式，采用正交曲線網(wǎng)格和網(wǎng)格雙重嵌套技術(shù)，建立了以平潭附近海域?yàn)橹攸c(diǎn)區(qū)域，涵蓋西北太平洋區(qū)域的精細(xì)化三維溫、鹽、流數(shù)值預(yù)報(bào)模式。模式驗(yàn)證結(jié)果表明，模式預(yù)報(bào)結(jié)果與實(shí)測(cè)結(jié)果比較一致，模式在平潭近海預(yù)報(bào)結(jié)果優(yōu)于沿岸海域。平潭沿岸海域模式預(yù)報(bào)流速絕對(duì)平均誤差和均方根誤差分別小于0.15 m/s和0.20 m/s，流向絕對(duì)平均誤差和均方根誤差分別小于35.4°和43.4°；平潭近海海域模式預(yù)報(bào)流速絕對(duì)平均誤差和均方根誤差分別小于0.11 m/s和0.15 m/s，流向絕對(duì)平均誤差和均方根誤差分別小于27.6°和28.2°。在模式的基礎(chǔ)上，建立了包含數(shù)據(jù)預(yù)處理、溫鹽流預(yù)報(bào)和預(yù)報(bào)產(chǎn)品后處理等子程序的業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)。業(yè)務(wù)化運(yùn)行結(jié)果表明：數(shù)值預(yù)報(bào)產(chǎn)品具有較高精度，整個(gè)業(yè)務(wù)流程運(yùn)行在30 min之內(nèi)，系統(tǒng)運(yùn)行穩(wěn)定，達(dá)到業(yè)務(wù)化預(yù)報(bào)的要求。該預(yù)報(bào)業(yè)務(wù)系統(tǒng)已業(yè)務(wù)化運(yùn)行2—3 a，多次為平潭海域海上落水人員、赤潮漂移軌跡預(yù)測(cè)提供海洋水動(dòng)力數(shù)據(jù)支持。同時(shí)，該系統(tǒng)也可以為海上軍事活動(dòng)、海上運(yùn)動(dòng)、海上交通運(yùn)輸?shù)忍峁┯行У暮Ｑ蟓h(huán)境安全保障，具有較高的應(yīng)用價(jià)值和社會(huì)經(jīng)濟(jì)價(jià)值。

[1]France GODAE OceanView Science Team.GODAE ocean-view national reports 2011:France Mercator[R].Paris:The 3rd meeting of the GODAE OceanView Science Team,2011.

[2]De Dominicis M,Bruciaferri D,Gerin R,et al.Amulti-model assessment of the impact of currents,waves and wind in modelling surface drifters and oil spill[J].Deep Sea Research Part II:Topical Studies in Oceanography,2016,133:21-38.

[3]Rhodes R C,Hurlburt H E,Wallcraft A J,et al.Navy real-time global modeling systems[J].Oceanography,2002,15(1):29-43.

[4]王輝,萬(wàn)莉穎,秦英豪,等.中國(guó)全球業(yè)務(wù)化海洋學(xué)預(yù)報(bào)系統(tǒng)的發(fā)展和應(yīng)用[J].地球科學(xué)進(jìn)展,2016,31(10):1090-1104.

[5]李云,劉欽政,王旭.海上失事目標(biāo)搜救應(yīng)急預(yù)報(bào)系統(tǒng)[J].海洋預(yù)報(bào),2011,28(5):77-81.

[6]吳玲娟,高松,徐江玲.黃渤海近岸精細(xì)化三維溫鹽流業(yè)務(wù)化數(shù)值預(yù)報(bào)系統(tǒng)研發(fā)[J].防災(zāi)科技學(xué)院學(xué)報(bào),2015,17(4):83-91.

[7]陳冠惠.風(fēng)對(duì)福建中部沿岸春、夏季水文結(jié)構(gòu)和上升流的影響[J].海洋科學(xué),1991,(4):48-53.

[8]黃榮祥,翁學(xué)傳.閩中沿岸上升流特征和成因分析[J].熱帶海洋,1996,15(1):52-59.

[9]林培根,陳照章,胡建宇,等.2011年夏季臺(tái)灣海峽上升流的變化特征及其與風(fēng)場(chǎng)的關(guān)系[J].臺(tái)灣海峽,2012,31(3):307-316.

[10]Hu J Y,Kawamura H,Li C Y,et al.Review on current and seawater volume transport through the Taiwan Strait[J].Journal of Oceanography,2010,66(5):591-610.

[11]江毓武,萬(wàn)振文,柴扉,等.臺(tái)灣海峽三維海流現(xiàn)報(bào)及預(yù)報(bào)系統(tǒng)[C]//海峽兩岸第七屆海洋科學(xué)研討會(huì)論文集.杭州:國(guó)家海洋局,2008.

[12]Haidvogel D B,Arango H G,Hedstrom K,et al.Model evaluation experiments in the North Atlantic Basin:simulations in nonlinear terrain-following coordinates[J].Dynamics of Atmospheres and Oceans,2000,32(3-4):239-281.

[13]Warner J C,Sherwood C R,Arango H G,et al.Performance of four turbulence closure models implemented using a generic length scale method[J].Ocean Modelling,2005,8(1-2):81-113.

[14]Mellor G L,Yamada T.Development of a turbulence closure model for geophysical fluid problems[J].Reviews of Geophysics and,1982,20(4):851-875.

[15]Leipper D F.Observed ocean conditions and Hurricane Hilda,1964[J].Journal ofAtmospheric Sciences,1967,24(2):182-196.

An operational three-dimensional numerical forecasting system of the sea temperature,salinity and water currents with fine resolutions for Pingtan coastal region

ZENG Yin-dong1,2
（1.Marine Forecasting Center of Fujian Province,Fuzhou 350003 China;2.College of Oceanic and Atmospheric Sciences,Ocean University of China,Qingdao 266071 China）

Based on geographic and hydrodynamic characteristics of Pingtan coastal region of Fujian Province of China,an operational three-dimensional numerical system for sea temperature,salinity and current forecasting with high horizontal and vertical resolutions have been developed.This paper details the numerical model,configuration,results validation,and operational running.This forecasting system is based on ROMS model,by using orthotropic curvilinear grids and nested grid techniques.This system is composed of pre-processing data module,forecasting module,and forecasting production visualization module.The forecasting system has shown reasonable results by comparing with observational data,and the operational running time of the system is constrained in 30 min.The system is stable,and has reached the requirement and standard to become an operational forecasting system.

Pingtan coastal region；fine resolution；sea current model；operationalization；numerical forecasting system

P731.3

A

1003-0239（2017）06-0039-09

10.11737/j.issn.1003-0239.2017.06.005

2017-06-07；

2017-07-21。

“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目（2013BAB04B00）；福建省科技計(jì)劃項(xiàng)目（2017Y0005）。

曾銀東（1978-），男，高工，博士在讀，主要從事海洋觀測(cè)預(yù)報(bào)工作。E-mail:zydzyd100@163.com