東中國(guó)海HYCOM模式垂向坐標(biāo)對(duì)比研究

陳曉斌 , 周 林 陳 璇, 鄭崇偉 吳炎成

(1. 解放軍理工大學(xué) 氣象海洋學(xué)院, 江蘇 南京 211101; 2. 92538部隊(duì)氣象臺(tái), 遼寧 大連 116041; 3. 75822部隊(duì), 廣東 廣州 510510)

隨著研究的深入, 越來(lái)越多的學(xué)者開(kāi)始嘗試?yán)没旌献鴺?biāo)系統(tǒng)克服傳統(tǒng)單一坐標(biāo)系在模擬海洋環(huán)流和溫、鹽分布上的缺陷。目前較為流行的混合坐標(biāo)大洋模式(hybrid ocean model)為 HYCOM[1]、NCOM[2]模式, 其中 HYCOM 采用三坐標(biāo)混合, 即 z坐標(biāo)、σ坐標(biāo)、密度坐標(biāo)。HYCOM模式在垂向坐標(biāo)的設(shè)置上具有靈活性, 為此該模式不但能夠應(yīng)用于深海大洋, 更可以利用其做近岸模擬, 國(guó)外很多學(xué)者已經(jīng)開(kāi)始了這個(gè)工作。

Winther等[3]利用HYCOM模式在北海及丹麥的斯卡格拉克海峽進(jìn)行了近岸垂向坐標(biāo)敏感性試驗(yàn),試驗(yàn)結(jié)果表明該模式不但可以用于準(zhǔn)確描述深海大洋及海盆尺度的模擬, 而且也可以用于沿海和陸架海區(qū)的模擬; Chassignet等[4]利用HYCOM模式在大西洋區(qū)域分別探討了垂向坐標(biāo)、參考?xì)鈮杭皽貕盒?yīng)對(duì)模擬結(jié)果的影響, 結(jié)果指出 HYCOM 模式坐標(biāo)選取具有相當(dāng)?shù)撵`活性, 研究的深度不同選用的參考?xì)鈮翰煌约澳Ｊ街锌紤]溫壓效應(yīng)會(huì)有效減小近表層及近底層的氣壓梯度力誤差。朱江、謝基平[5]在2012年利用HYCOM模式的2.1及2.2版本對(duì)中國(guó)海區(qū)進(jìn)行了敏感性試驗(yàn), 指出了模擬中國(guó)海近岸及陸架海結(jié)構(gòu)需要準(zhǔn)確的風(fēng)場(chǎng)、全球及局地的高分辨率地形數(shù)據(jù), 能夠準(zhǔn)確刻畫(huà)出近岸射流與渦旋的高水平分辨率(1～3 km)及包含潮汐過(guò)程的大洋模式。

由于東中國(guó)海的地形復(fù)雜多變[6], 混合坐標(biāo)模式能在一定程度上解決這個(gè)問(wèn)題。本文的目的在于探討 HYCOM模式中不同垂向坐標(biāo)設(shè)置對(duì)東中國(guó)海近岸區(qū)域的影響以及黑潮流速及路徑對(duì)不同坐標(biāo)設(shè)置的響應(yīng), 期望對(duì) HYCOM 模式更深入的研究提供參考, 為軍地海洋建設(shè)提供科學(xué)依據(jù)[7-8]。

1 模式介紹

1.1 模式簡(jiǎn)介

NERSC-HYCOM 模式[9]是由挪威的南森環(huán)境遙感中心基于標(biāo)準(zhǔn)版 HYCOM 模式并于 2009年發(fā)布的。在網(wǎng)格選取上, 為了使模擬區(qū)域網(wǎng)格加密, 進(jìn)行了正交投影坐標(biāo)轉(zhuǎn)換[10]; 模式中的兩種邊界條件,分別是牛頓張弛邊界條件和開(kāi)邊界條件。在海綿邊界帶(即有一定寬度的緩沖帶, 以減少側(cè)邊界條件的反射), 應(yīng)用了一個(gè)簡(jiǎn)單的牛頓張弛邊界條件; 溫度、鹽度和垂直坐標(biāo)的壓強(qiáng)都隨時(shí)間步長(zhǎng)進(jìn)行更新。當(dāng) HYCOM 采用等密度面坐標(biāo)時(shí), 溫度和鹽度僅在非等密度面混合層(第一層)進(jìn)行張弛運(yùn)算, 下面的深層僅對(duì)鹽度進(jìn)行張弛運(yùn)算; 當(dāng)采用混合坐標(biāo)時(shí),溫度和鹽度在上面的混合坐標(biāo)層進(jìn)行, 下面的深層也是僅對(duì)鹽度進(jìn)行張弛運(yùn)算; 而壓強(qiáng)則是對(duì)整個(gè)垂直坐標(biāo)層進(jìn)行張弛運(yùn)算。

1.2 模式基本設(shè)置

為減少模式邊界取值誤差對(duì)東中國(guó)海區(qū)模擬結(jié)果的影響, 本文采用了單向嵌套方案, 即將東中國(guó)海計(jì)算區(qū)域(以下稱為小區(qū))嵌套于西北太平洋計(jì)算區(qū) 域(以 下 稱 為 大 區(qū))。 大 區(qū) 范 圍 95°E～90°W,28°S～51°N, 分辨率為 0.5°×0.5°, 經(jīng)向格點(diǎn) 370 個(gè),緯向格點(diǎn) 178 個(gè); 小區(qū)范圍 118°～150°E, 23.5°～44°N,分辨率為0.125°, 經(jīng)向格點(diǎn)256個(gè), 緯向格點(diǎn)200個(gè)。垂直分層 23層, 最小水深 10 m。地形數(shù)據(jù)采用ETOPO5數(shù)據(jù), 模式強(qiáng)迫場(chǎng)資料分別采用了ECWMF提供的ERA40氣候態(tài)資料與2008年的0.5°分辨率的ERA-I高頻強(qiáng)迫資料并提供了長(zhǎng)江、黃河的流量數(shù)據(jù);溫鹽初始場(chǎng)替換一貫采用的levitus資料而改用質(zhì)量較高的PHC資料; 四個(gè)敏感性試驗(yàn)都采用了KPP湍封閉方案; 大小區(qū)域模式均采用并行計(jì)算方案以提高模式的運(yùn)行效率。

由于本文研究的是海洋上層的物理過(guò)程, 因此采用海表參考?jí)簭?qiáng), 位勢(shì)密度值的前5層值小于任何海水的位密, 目的是使這5層固定在海表面, 確保上混合層的精度[11]。HYCOM模式在氣候態(tài)風(fēng)場(chǎng)驅(qū)動(dòng)下在大區(qū)域穩(wěn)定積分10 a后利用2008年高頻(0.5°水平分辨率)的 ERA-I風(fēng)場(chǎng)等資料繼續(xù)積分模式1 a整, 然后利用其模擬結(jié)果提供邊界條件驅(qū)動(dòng)HYCOM 模式在小區(qū)運(yùn)行1 a, 最后對(duì)小區(qū)模擬結(jié)果進(jìn)行分析討論。

模式采用了WOA09資料、ARGO及XBT資料用于進(jìn)行模式氣候態(tài)運(yùn)行以及高頻強(qiáng)迫實(shí)驗(yàn)的對(duì)比。

2 垂向坐標(biāo)配置試驗(yàn)方案

2.1 位勢(shì)密度分析

為便于研究東中國(guó)海的垂向坐標(biāo)設(shè)置對(duì)模式模擬結(jié)果的影響并考慮到 HYCOM模式是基于等密度坐標(biāo)模式改進(jìn)的, 為此有必要對(duì)該區(qū)域的位勢(shì)密度進(jìn)行分析。針對(duì)WOA09氣候態(tài)的季平均資料進(jìn)行了位勢(shì)密度的分析, 在 29.5°N斷面(圖略)及 123.5°E 斷面(圖略)中東中國(guó)海區(qū)密度的鉛直向分布是與海水層結(jié)相應(yīng)的, 一般是隨深度的增加而增大。近岸淺水區(qū)位密存在著十分明顯的季節(jié)性變化, 冬季位密大、夏季位密小, 其主要是受海溫的季節(jié)性變化影響;而靠近長(zhǎng)江口附近夏季的位密較其他區(qū)域更低, 這是由于夏季長(zhǎng)江口附近的河流入海淡水通量較冬季明顯增多引起的。由于東中國(guó)海區(qū)的水深較淺致使該區(qū)混合層的深度淺且位密的變化十分明顯, 可以認(rèn)為在東中國(guó)海區(qū)上層并不適于采用等密度坐標(biāo)方案來(lái)描述混合層的特點(diǎn), 應(yīng)該采用z坐標(biāo)或σ坐標(biāo)用以準(zhǔn)確表征此處混合層的季節(jié)性變化特征, 為此本文據(jù)此進(jìn)行了如下的敏感性試驗(yàn)。

2.2 垂向坐標(biāo)設(shè)置試驗(yàn)方案參數(shù)設(shè)計(jì)

考慮到模式靈活的垂向坐標(biāo)設(shè)置[3-4,12], 并為了比較在高頻強(qiáng)迫下不同混合坐標(biāo)的設(shè)置對(duì)模式模擬結(jié)果的影響, 本文建立的四個(gè)垂向坐標(biāo)試驗(yàn)為z-only(40 層)、σ-only(40 層)、σ-z-iso(23 層)以及 z-iso(23層), 其中對(duì)比資料采用Argo溫鹽廓線資料、XBT溫度廓線資料。表1為模式的方案設(shè)計(jì), 對(duì)于含有密度坐標(biāo)的試驗(yàn), 位密的選取為[13]: 0.10, 0.20, 0.30,0.40, 0.50, 23.20, 23.94, 24.64, 25.22, 25.70, 26.11,26.44, 26.72, 26.95, 27.14, 27.30, 27.44, 27.56, 27.66,27.74, 27.80, 27.84, 27.87。

表1 垂向坐標(biāo)設(shè)置試驗(yàn)方案參數(shù)設(shè)計(jì)Tab.1 Parameter design of testing programs for vertical coordinate setting

其中δpmin, δpmax, fp指在深水區(qū)中z坐標(biāo)分層的每層最小厚度及最大厚度、比例因子; δfmin, δfmax, fs指在淺水區(qū)中 z坐標(biāo)分層的每層最小厚度及最大厚度、比例因子; Nσ代表采用σ坐標(biāo)的分層數(shù)目。

試驗(yàn)1: 將模式設(shè)計(jì)為z-only坐標(biāo)模式。

試驗(yàn)2: 將模式設(shè)計(jì)為σ-only坐標(biāo)模式。

試驗(yàn)3: 將模式設(shè)計(jì)為σ-z-iso坐標(biāo)混合模式, 考慮到在淺水區(qū)域有一個(gè)σ坐標(biāo)的過(guò)渡。

試驗(yàn)4: 將模式設(shè)計(jì)為z-iso坐標(biāo)混合模式。

3 垂向坐標(biāo)配置對(duì)東中國(guó)海溫鹽垂直廓線和表層環(huán)流模擬的影響

3.1 溫鹽廓線與誤差分析

東海大陸架區(qū)域水深地形十分復(fù)雜, 近海部分水深100 m左右、陸架區(qū)深度200～300 m, 而日本九州島以西水深達(dá)到2 000 m左右, 九州島以東水深平均達(dá)到5 500 m; 該區(qū)海洋水文環(huán)境主要受黑潮的影響, 它沿著陸架坡折處的200～300 m等深線流動(dòng), 首先進(jìn)入臺(tái)灣島以東, 然后在日本九州島以南流出。

圖1中, 選取大陸架區(qū)域 2008年 2月的 5個(gè)ARGO觀測(cè)點(diǎn)(A～E)及2008年8月的5個(gè)XBT觀測(cè)點(diǎn)(F～J), 點(diǎn) A(130.69°E, 28.27°N)、點(diǎn) B(128.89°E,24.51°N)、點(diǎn) C(130.03°E, 27.61°N)、點(diǎn) D(130.56°E,25.49°N)、點(diǎn) E(131.04°E, 25.24°N), 點(diǎn) F(128.87°E,29.00°N)、點(diǎn) G(128.17°E, 29.56°N)、點(diǎn) H(128.42°E,29.33°N)、點(diǎn) I(128.83°E, 30.00°N)、點(diǎn) J(128.78°E,29.67°N)。

圖1 浮標(biāo)點(diǎn)分布Fig.1 Distribution of buoys points

圖2分別繪制了A～E 5個(gè)站點(diǎn)的Argo與4個(gè)數(shù)值試驗(yàn)的溫鹽垂直廓線的分布圖, 這5個(gè)點(diǎn)所在的區(qū)域深度可達(dá)幾千米, 可以發(fā)現(xiàn)試驗(yàn)1, 試驗(yàn)3, 試驗(yàn)4與Argo溫鹽廓線的趨勢(shì)最為接近, 且隨著深度的增加, 模擬效果越好。溫度方面, 三個(gè)試驗(yàn)在200 m以上的海洋上層以及1 000 m以下的海洋深層的模擬較好, 而在200～600 m這個(gè)范圍與Argo溫度廓線差異較大, 試驗(yàn)2在整個(gè)深度范圍模擬的溫度廓線與 Argo差異都十分明顯; 鹽度方面, 試驗(yàn)1, 試驗(yàn)3, 試驗(yàn)4三個(gè)試驗(yàn)基本模擬出了 Argo鹽度廓線的趨勢(shì), 但是500 m以上區(qū)域并沒(méi)準(zhǔn)確刻畫(huà)出Argo觀測(cè)的結(jié)果, 存在0.3左右的負(fù)偏差, 即模擬值低于觀測(cè)值, 試驗(yàn)2模擬的鹽度廓線則根本無(wú)法體現(xiàn)Argo鹽度廓線的趨勢(shì)。

為了評(píng)估 NERSC-HYCOM 模式中四個(gè)試驗(yàn)?zāi)M溫鹽場(chǎng)的效果, 表2計(jì)算了站點(diǎn)A～E四個(gè)試驗(yàn)相對(duì)于 Argo溫鹽廓線資料的平均誤差(ME)、均方根差(RMS)、相關(guān)系數(shù)(R), 表中的T, S分別表示溫度和鹽度。結(jié)果指出: 在陸架地區(qū)的深水區(qū), 可以得出試驗(yàn)3(23層混合坐標(biāo)試驗(yàn))與試驗(yàn)1(40層 z-only試驗(yàn))刻畫(huà)的溫鹽場(chǎng)差異不大, 都能夠比較好地反應(yīng)實(shí)測(cè)的溫鹽場(chǎng)信息。

圖3繪制了集中在東海陸架及陸架坡處5個(gè)點(diǎn)的 XBT觀測(cè)的溫度廓線(鹽度觀測(cè)資料缺失)與試驗(yàn)1～試驗(yàn)4的模擬的溫度廓線。由于陸架區(qū)地形崎嶇復(fù)雜, 溫度廓線的分布并不平滑, 0～500 m范圍內(nèi)四個(gè)試驗(yàn)方案模擬的溫度都與XBT資料存在一定的偏差[14-15]。同前述深水區(qū)A～E五個(gè)觀測(cè)點(diǎn)的分析相似,通過(guò)表3的數(shù)據(jù)分析, 可以發(fā)現(xiàn)試驗(yàn)3還是能夠較好地?cái)M合陸架區(qū)的溫度趨勢(shì)。試驗(yàn)2模擬的溫度場(chǎng)比較平滑而且與XBT資料差異最大。這也說(shuō)明在陸架及陸架坡這些地形較為復(fù)雜的地區(qū) HYCOM模式的三重混合坐標(biāo)模擬能力的優(yōu)越性得到了很好的體現(xiàn);試驗(yàn)4的效果總體不如試驗(yàn)3這也恰恰解釋了HYCOM 模式在地形復(fù)雜陸架坡處引用 σ坐標(biāo)混合的合理性; 40層的σ-only坐標(biāo)試驗(yàn)與XBT資料存在較大差異, 則體現(xiàn)了單一的σ坐標(biāo)用于海溫計(jì)算, 在地形比較陡峭的陸架過(guò)渡區(qū)、地形陡峭的陸架坡有很多缺點(diǎn)(最大的缺點(diǎn)是斜壓梯度力的計(jì)算誤差)[16]。

3.2 東海黑潮表面流場(chǎng)分析

黑潮作為世界上最強(qiáng)的西邊界流之一, 在日本南部流速最大可達(dá) 1.5～2.0 m/s, 黑潮也能發(fā)生大彎曲[17-18]。圖4給出了2008年以CORA資料表征的研究海域的冬季(以2月為代表)、夏季(以8月為代表)的表層流場(chǎng)。圖中顯示冬夏兩季黑潮流十分明顯, 在親潮處可以看到多個(gè)中尺度渦。夏季最大流速可達(dá)到1.6470 m/s, 冬季最大流速可達(dá)1.0514 m/s; 冬夏兩季可以看到顯著的對(duì)馬暖流, 夏季流速?gòu)?qiáng)于冬季。

圖2 5個(gè)站點(diǎn)觀測(cè)的溫鹽垂直廓線與對(duì)應(yīng)的數(shù)值模擬結(jié)果對(duì)比Fig.2 Profile of thermohaline between Argo and Model for five points

表2 模擬陸架地區(qū)深水區(qū)溫度、鹽度精度檢驗(yàn)Tab.2 Precision of simulation temperature and salinity between Argo and Model in the deep water area of the continental shelf

圖3 5個(gè)站點(diǎn)上的XBT- MODEL溫度廓線Fig.3 Profile of temperature between XBT and Model for five points

表3 模擬陸架地區(qū)地形復(fù)雜區(qū)溫度精度檢驗(yàn)Tab.3 Precision of simulation temperature between XBT and Model in the complex terrain area of the continental shelf

圖5為試驗(yàn)1～試驗(yàn)4模擬的東太平洋海域表層流場(chǎng)。圖中可以看出試驗(yàn)3與試驗(yàn)4在132°～135°E之間模擬出了非常清晰的黑潮“大彎曲”路徑, 黑潮流軸在30°N附近東轉(zhuǎn), 并沿著日本島南岸及東南岸向東北方向流動(dòng), 在(40°N, 142°E)形成了黑潮的延伸體, 與 CORA資料冬夏季流場(chǎng)相比, 整個(gè)黑潮路徑的趨勢(shì)及日本島以南的反氣旋渦模擬較好, 但是可以看到兩個(gè)試驗(yàn)?zāi)M的黑潮主軸偏南, 兩個(gè)試驗(yàn)流速最大值偏高, 尤其日本東北部的渦旋強(qiáng)度較強(qiáng),且渦旋路徑模擬的偏北; 試驗(yàn)1模擬的黑潮主軸位置與 CORA 資料較為接近, 和試驗(yàn)3、試驗(yàn)4相比,試驗(yàn)1對(duì)日本島以南黑潮的模擬路徑有所改善。試驗(yàn)2模擬的黑潮路徑偏差最大, 效果最差。雖然模擬黑潮路徑受模式的水平分辨率的影響較大[19], 本文認(rèn)為在水平分辨率相同的情況下z坐標(biāo)的垂直分辨率對(duì)表層的黑潮路徑影響也很大, 即通過(guò)增加上混合層內(nèi)z坐標(biāo)的分層可以使黑潮路徑的模擬達(dá)到更好的效果。

4 結(jié)論

本文利用 NERSC-HYCOM 模式, 選取了 z-iso,σ-z-iso, z-only及 σ-only四個(gè)不同垂向坐標(biāo)結(jié)構(gòu), 采用單向嵌套方式, 對(duì)東中國(guó)海區(qū)進(jìn)行數(shù)值模擬研究,結(jié)論如下。

1) 東中國(guó)海區(qū)上層并不適于采用等密度坐標(biāo)方案, 應(yīng)采用z坐標(biāo)或σ坐標(biāo)以表征此處混合層的季節(jié)性變化特征。

2) 數(shù)值模擬結(jié)果與 Argo溫鹽廓線、XBT溫度廓線對(duì)比分析發(fā)現(xiàn), 對(duì)于溫度場(chǎng)而言, 試驗(yàn) 3、試驗(yàn)4和試驗(yàn)1與Argo溫度廓線、XBT溫度廓線的趨勢(shì)最為接近, 且隨著深度的增加, 模擬效果越好; 對(duì)于鹽度場(chǎng), 四個(gè)試驗(yàn)在表層區(qū)域的鹽度都低于Argo資料且存在較大的偏差, 但是試驗(yàn)3、試驗(yàn)4和試驗(yàn)1則能夠很好吻合實(shí)際的鹽度分布態(tài)勢(shì)。在所選取的10個(gè)觀測(cè)點(diǎn)中, 試驗(yàn)2的結(jié)果與實(shí)測(cè)資料偏差最大,效果最差。

圖4 源于CORA資料的表層流場(chǎng)Fig.4 Surface flow field from CORA data

圖5 試驗(yàn)1～試驗(yàn)4模擬的東中國(guó)海區(qū)域表層流場(chǎng)Fig.5 Surface flow field from four tests in the East China Sea

3) 對(duì)黑潮及周邊流場(chǎng)的模擬結(jié)果表明: 試驗(yàn) 1模擬的黑潮主軸位置與CORA資料最為接近; 試驗(yàn)2模擬的黑潮路徑偏差最大, 效果最差; 試驗(yàn)3與試驗(yàn) 4對(duì)整個(gè)黑潮路徑的趨勢(shì)及日本島以南的反氣旋渦模擬相對(duì)較好, 但是模擬的黑潮主軸偏南, 流速最大值普遍偏高, 尤其日本東北部的渦旋強(qiáng)度較強(qiáng)。

4) NERSC-HYCOM 模式可以模擬東中國(guó)海水深地形崎嶇復(fù)雜海域的溫鹽垂直分布和表層流場(chǎng)。本文采用密度坐標(biāo)的混合坐標(biāo)實(shí)驗(yàn)方案模擬的黑潮流系有待改進(jìn), 提高上混合層z坐標(biāo)分辨率的方式可能會(huì)使黑潮的流場(chǎng)模擬更為合理。

[1] 鄭沛楠.常用海洋數(shù)值模式簡(jiǎn)介[J].海洋預(yù)報(bào), 2008,25(4): 108-120.

[2] 鄭沛楠.黑潮對(duì)日本海邊界環(huán)流的影響[D].青島: 中國(guó)海洋大學(xué), 2009.

[3] Nina G, Winther G E. A Hybrid Coordinate Ocean Model for shelf sea simulation[J].Ocean Modelling,2006, 1: 221-237.

[4] Chassignet E P, Smith L T, Halliwell G R, et al. North Atlantic simulation with the Hybrid Coordinate Ocean Model (HYCOM): Impact of the vertical coordinate choice, reference density, and the Rmobaricity[J]. Phys Oceanogr, 2003, 33: 2504-2526.

[5] Jiang Z, Xie J P. Modelling and data assimilation in Marginal Seas around China[OE/BL].[2013-03-17].http: //www.godae.org/...TT.../2.9-Zhu-Jiang-2012-Miami.pdf.

[6] 陳敏, 侯一筠, 趙保仁.冬季東中國(guó)海環(huán)流中的中尺度渦數(shù)值模擬[J].海洋科學(xué), 2003, 27(1): 53-60.

[7] 鄭崇偉, 游小寶, 潘靜, 等.釣魚(yú)島、黃巖島海域風(fēng)能及波浪能開(kāi)發(fā)環(huán)境分析[J].海洋預(yù)報(bào), 2014, 31(1):49-57.

[8] 鄭崇偉, 潘靜, 黃剛.利用WW3模式實(shí)現(xiàn)中國(guó)海擊水概率數(shù)值預(yù)報(bào)[J].北京航空航天大學(xué)學(xué)報(bào), 2014, 40(3):314-320.

[9] Francois K. NERSC-HYCOM 2.2[OE/BL]. [2013-03-17].https: //svn.nersc.no/hycom.

[10] 于慶龍, 王輝, 萬(wàn)莉穎. 南海海表溫度時(shí)空分布特征的數(shù)值模擬[J].海洋預(yù)報(bào), 2010, 27(4): 59-66.

[11] George R H. Evaluation of vertical coordinate and vertical mixing algorithms in the Hybrid—Coordinate Ocean[J]. Ocean Modelling, 2004, 7: 285-322.

[12] 王驥鵬, 聞斌, 耿再興, 等. 混合坐標(biāo)大洋環(huán)流模式Hycom的坐標(biāo)選取與參數(shù)化設(shè)置敏感性試驗(yàn)[J].海洋預(yù)報(bào), 2008, 25(3): 93-102.

[13] 白志鵬, 高松, 王海棠. HYCOM 模式對(duì)東海黑潮的氣候態(tài)模擬[J]. 海洋通報(bào), 2010, 29(2): 121-129.

[14] Thacker W C, Esenkov O E．Assimilating XBT data into HYCOM [J]．Journal of Atmospheric and Oceanic Technology, 2002, 19: 709-724.

[15] Thacker W C, Lee S K, Halliwell G R．Assimilating 20 years of Atlantic XBT data into HYCOM: a first look[J].Ocean Modeling, 2004, 7: 183-210．

[16] 馮興如, 楊德周, 尹寶樹(shù). FVCOM 在龍口海域潮汐潮流模擬中的應(yīng)用研究[J].海洋科學(xué), 2012, 36(4):107-112.

[17] 許雪峰, 羊天柱, 孫志林, 等. 杭州灣南岸淺灘的漫灘潮流數(shù)值模擬[J].海洋科學(xué), 2010, 34(6): 94-99.

[18] 劉增宏. Argos漂流浮標(biāo)的若干觀測(cè)結(jié)果[J].熱帶海洋學(xué)報(bào), 2005, 24(1): 67-76.

[19] Guo X Y. A triply nested Ocean Model for simulating the Kuroshio—Roles of horizontal resolution on JEBAR[J].Journal of Physical Oceanography, 2003, 33: 146-169.