影响海上军事活动的海洋水文要素主要包括：海浪、风暴潮、潮汐潮流、海流、海温、盐度、内波、水声、海冰等。由于海洋环境资料获取困难，极大地制约了海洋水文预报水平，海洋水文预报主要依赖数值预报。影响海浪、风暴潮、潮汐潮流的物理机制相对来说比较简单，它们的预报技术目前已很成熟，基本能够满足业务化应用的要求。影响海流、海温、盐度的物理机制则要复杂很多，其预报难度也要大很多。海流、海温、盐度预报又是开展内波、水声、海冰预报的基础。此外，海洋中尺度现象对海洋动力环境、军事活动也有重要的影响，针对海洋中尺度现象的研究近年来越来越得到重视。海洋中尺度现象主要包括中尺度涡、海温锋面、盐度锋面、跃层、上升流、内
波等。海洋中尺度现象预报水平也主要取
决于海流、海温、盐度预报。关于海流、海温、盐度数值预报的模式通常称为海洋(大洋)环流模式，它一直是海洋科学研究的核心问题。
本文主要介绍被美海军青睐的新一代混合坐标海洋环流模式H Y C O M 。
1 混合坐标的提出及HYCOM 模式简介
海洋环流数值模式从单一坐标到混合坐标、从二维模型到三维模型，数值模式的发展越来越快。尽管单一坐标模式的应用比较广泛，但是国外的C HA SS IG N ET EP [1]chainsaw
进行的大洋环流比较试验表明，单一垂向坐标用来模拟海洋环流并不总是最佳选择，因此混合坐标的概念就被提出了。考虑到海底复杂的地形，国内外很多学者都开展了多种混合坐标模式的研究工作，本文采用的H Y C O M 三坐标混合模式越来越受到重视，国内越来越多的学者也将其作为主流模式。1.1标准版HYCOM
标准版的混合坐标大洋环流模式HYCOM [2,3](Hybrid Coordinate Oceanic Circulation Model)是在美国迈阿密大学原有等密度面海洋模式M I C O M 的基础上发展改进的新一代原始方程海洋环流模式(H all iwel l,2004)。模式的水平网格与M I C O M 一样采用C 网格。
HYCOM模式[4]可以实现三种坐标的自适应，即在开放的、层结的海洋中垂直坐标取等密度坐标；在弱层结的上层海洋混合层中，垂直坐标平滑地过渡到z 坐标；在浅水区域(地形变化强烈的海域)，过渡到随底垂向坐标；而在更浅的海域坐标又过渡到z 坐标。混合坐标扩展了传统的等密度面坐标海洋模式的应用范围，弥补了M IC OM 数值模式在浅海海域垂向分层过薄的缺点(Wallcraft,2003)，使得模式的适应能力更好。
模式提供了两种边界条件[4]：一种是牛顿松弛边界条件；另一种为开边界条件。
H Y C O M 数值模式还提供了多种湍封闭方案[4]，来解决上混合层和层化相对较弱区域的跨等密度面混合问题。如K -T (K r a u s -T u r n e r ),K P P (K -P r o f i l e -Parameterization),PWP(Price-Weller-
Pinkel),MY-2.5(Mellor-Yamada level 2.5),GISS(NASA Goddard Institute for Space Studies level 2)(Large,et a1．,1994 and 1997;Eric P,2003;Halliwell Jr．,2004)方案，通过实验结果表明不同的方案在气候态的模拟中存在较大的差异，其中非局地K P P 参数化方案、GISS湍封闭方案和MY湍封闭方案的模拟效果较好。1.2NERSC-HYCOM
挪威版N E R S C -H Y C O M 模式[5]是由挪威的南森环境遥感中心基于标准版H Y C O M 模式而改进的。在网格选取上，为了使模拟区域网格加密，进行了正交投影坐标转换[6]；虽然该模式基于标准版但在嵌套、潮汐边界、海冰模式、强迫场的处理中采用了不同的方法；对于数据的读取而言，该模式可直接处理N E T C D F 格式数据，这就使得模式的操作较为简化；挪威版模式的开发为用户提供了一个更为友好的界面，大量的脚本使得用户不必深究模式代码却能较容易的建立模式；为了能够及时对新版本的模式代码做出相应的升级更新，由南森中心改进的代码已从H YC OM 代码中单独分离出来。N E R S C -H Y C O M 模式
①作者简介：陈晓斌(1985—)，男，辽宁沈阳人，硕士研究生，主要从事海气相互作用研究。
混合坐标大洋环流模式H Y C O M 及其研究进展①
陈晓斌 刘志宏 林刚 李柄更
(92538部队气象台 辽宁大连 116041)
摘 要：传统单一坐标由于自身的优缺点，使得利用单一坐标的海洋模式的模拟效果受到了限制。混合坐标的出现，正是为了综合z坐标、σ坐标、等密度坐标的优势并尝试解决这一难题而提出的重要手段。本文基于美国海军研发的HYCOM模式以及挪威南森环境遥感中心改进的NERSC-HYCOM模式，针对当前HYCOM模式的应用现状进行了综述。关键词:HYCOM NERSC-HYCOM 混合坐标中图分类号:P731.27文献标识码:
A 文章编号:1672-3791(2013)12(a)-0207-04
图1 HYCOM 模式混合坐标示意(V.Garnier and E.P.Chassignet,2004)
.com. All Rights Rerved.
tequila现在的版本为2.2.37，现有的版本都是并行版本并被广泛用到海洋环境预报当中。
2 垂向坐标的调整
2.1垂向坐标在开阔大洋的调整
垂向坐标在开阔大洋的调整一般是指等密度坐标的调整
[7,8]
。在HYCOM中 为比
容，开阔大洋中应该满足 1k k (下层密度比上层大)，如果取三个连续的等密度层0，层1，层2，而 1
不等于模式所选取的位密
参考值为了恢复等密度条件，水体必须
重新离散以保证水体的总高度。质量守恒关系使得一层或多层界面被重新分配到不同的层。在调整层1密度的时候，如果层1密
度太大(
1 移动使得在层0转移较轻的水到层1；相反如
果层1密度太小( 1
就会向下使得在层2转移较重的水到层1。
H Y C O M 中，通过指定 n hyb 确定混合坐标选用的层数。由 n hyb 所确定的较高层(upperlayer)的层界面通过下面的算法进行调整，对于更深的模式层则以等密度坐标的形式存在，调整大致分为两种情况:
情况1:(
南京教育网1 团会在层0到层1之间发生移动与改变。由
dp 的质量守恒可得到:
10121p p p p
010121p p p p (1)
其中1p 是格点重新调整之后层1上层层界面处压强,
1012110
p p p
(2)
因为 2p 的权重为负值，会得出1p 0p ，为了保证层0的最小厚度，将1p 替换为
1100
max(,)p p p (3)
这里 0 是一个指定的最小层厚度。当然，利用这种替换关系也使得不会再出现等密度的情况了。根据Bleck和Benjamin (1993)的相关工作，可以选择一个绝对最小厚度参数 0 ，那样 0 的计算就利用了模式提供的连续的分段缓冲(cushion)函数来实现，即当 0
3
p ≤时,
2
000133p p
；而当 03p 时, 0 p (其中 p 1p p
牛津小学英语1a110p p )即 1p 1p 1
p ，从而实现了最小厚度 0 会按照等密度坐标与非等密度坐标之间平滑过渡的方式进行计算。实际上，如果温盐场
以通量的形式穿越了重新分
配的层
界面，那么等密度坐标的恢复就会变得不可能还有一些极少发生的特殊情况,
1 的值会有不同的符号，例如：在地中海盐
舌之下的模式层中，就会发生了抬升了层1界面而减少 1 值等情况。
重复的使用垂向坐标调整算法会使 1 的值过度的偏离层1的参考密度值进而产生很大的垂向坐标移动；尤其是当温盐以平流方式输运的时候这个问题会更加严重。因此模式提供了代码以阻止 1 的值会有不同的符号的情况发生。尽管这种方法可以避免 1 的值会迅速偏离参考密度值，但是参考值的恢复是不可能的，唯一避免这个问题的方法是调整温度以及密度，或者盐度以及密度。
在模式中由于温、盐的重新调整依然会存在偏离层参考密度值的现象，为了尽
量解决这个问题模式中增加了一个迭代程
序确保在考虑使用垂向坐标算法的时
候
通过(2)估算出 1p 的时候，温盐以通量的形式穿过重新调整的垂向分层后所产生的温度、盐度改变会被估算，然后通过状态方程
估算出内，那么 1p 与新计算，直至符
合要求。情况1( 1
1
情况:(1) 1p 1p 表示如果上层界面抬升，使得 1 增加为于参考密度值 1对于 1(1)的基础上将倒三角符号替换为波浪号:
12101121
()()
p p
p p p p (4)
如果选择了温度、盐度的混合，那么下一步就在规定的范围内判
定 的值是否会发生符号的改变。当抬升接口时所产生的新温度值为
1110121
()p p T T T p p
(5)
新的盐度值 1S 采用同样的方式估算.如果 1 的符号错误，那么 1p 就不被调整；否 1
p 的值就会通过下面的式子重新计算:
11112
1
0p p p p (14)
其中 1p 不允许超过 1p 的值。(2) 11
p p 这种情况下，上层层界面被允许降低从而保证层0的最小厚度，利用 dp 的守恒性导
出
10
通过增加层0厚度的方式来
确保最小层厚度的保守性使得等密度坐标的恢复成为不可能，所以该种情况并没有用到任何
1p 1
1
p p 时提到的算法。
情况2:( 1 下移动，水团会从层2进入层1，由 dp 的质量守恒可得到:
121232()p p p p 1
2
12
3
2
朗文词典p
p p p
(8)
低层层界面需要被重新调整为
111
21222
)()p p p
(9)
而为了保证层2的最小厚度值, 2p 的值被取代为:
2232
min(,)p p p (10)
其中 2 通过缓冲分段函数来确定：如
果 p ≤
23 ,
22222133p p
；如果 2
3p
, 2
p
(其中 32
p p p )。当温
盐发生调整的时候，迭代算法会被程序调用以保证等密度坐标的恢复。情况2中还有一种特殊情形需要考虑，就是当层1为底层时。该种情况下，层界面2是不能向下移动
的，因此非常有必要使用(1)式并且向下移动层界面1去恢复层1的等密度情况。为了达到这样的目的，层1的水体必须被重新限制为两个子层以便于上部分的子层密度等于层0的密度，而第二个子层的密度值应该接近于参考密度值，并且两个子层的平均密度应该等于原始的层密度值。因此假定:
110
q 层1被重新分配为:
1p 1121+()p p q p p (12)
在第二个子层中的温度表示为:
11101q T T T T q
(13)
该层的盐度和密度则通过类似的方程进行计算得到。另外，热动力学变量中的两个量可以从(13)式诊断得到，而第三个量则通过状态方程得到。这样做可以达到两个目的:(1)为了避免层1的厚度减小超过50%，采用了:
1122
pollyanna
min(,)p p p (14)
其中 221()/2p p ;
(2)为了避免温盐场在河床的变化，需
.com. All Rights Rerved.
满足条件
1001T T T T ≤ (15)
2.2垂向坐标在近岸区域的调整
H Y C O M 模式中垂向坐标在近岸区域的调整可以认为是σ坐标及z 坐标的自适应过程。在模式中存在一个最小厚度项被用于实现从开阔大洋(其非密度坐标存在于大洋内部等密度层上部)向浅水区域σ坐标以及更浅的海域z坐标的转换。在模式的每一层中都指定一个最小厚度 n ，那样就会在等密度层上产生z坐标；为确保在浅水区有一个适当的坐标转换，通过 s N 及
min 去分别指定σ层数与绝对最小厚度。而
σ坐标的最小厚度表示为
s s D N
，
其中D是
总水深。然后再重新进行计算每个模式层的 min max ,min ,n s 。
在给定的模式分层中，当水深变得足够浅的时
候，向sigma坐标转换的满足条件为 s ＜s n ;
而模式坐标转换回z坐标的满足条件是 s ＜min s
。替 n 考虑坐标转换才是适当的。
3.3温、盐、密度及动量的调整
热力学变量以及动量的调整理论上必须遵守垂向平均值守恒、对于热力学变量应该尽可能的调整密度值与等密度参考值
接近,H YCO M中提供的算法会按照上述两点原则将每个变量从原先的格点重新映射到新格点上。模式中可以调整热力学变量中的任何两个并通过状态方程诊断出第三个量。当采用混合坐标设置时，会分别在每个计算格点上进行调整，对于热力学变量会优先在气压格点上调整。如果采用全K-T混合层方案，含有混合层底的模式层会被分为两个次层。混合层基带会被暂时认作另一种垂向坐标，基带上下的两个次层也会被暂时看做两个模式层，在两个次层中的热力学变量会根据K-T方案中的非混合算法进行估算。在进行垂向坐标调整之前，初始的一维温度、盐度、密度廓线以及一维层界面处压力数组会被储存。一旦廓线数组被储存，在气压格点上的垂向坐标的垂向调整就开始进行，然后每个变量会被重新映射到新的垂向坐标上来。下面用温度场来说明这种重新映射的关系:
1
11
1max ,min ,Mempress
n m
n
m n m n n
T T p p p p p
1min ,max ,n m n p p p
(16)
其中在原始格点上的层界面处压力为
he
m p (m取值范围为1,M)；对于新格点上的层
界面处压力为 n p (n取值范围为1,N)，在这里N 为模式的层数，现在讨论的情况中为M =N 。需要指出的是，如果采用全K -T 混合层方案，则用M =N +1来表示在含有混合层基带的模式层中的两个次层。
经过热力学变量调整后，动量部分将会进行调整以保证每个格点上的垂向平均动量的守恒，而通过压力格点获得的新旧垂向坐标会通过上述方程被依次差值到U 、V 格点上。
3 HYCOM 模式的研究进展
需要说明的是挪威版模式刚刚公开代码，因此国内外对N E R S C -H Y C O M 模式的研究较少，下面的内容更多是基于标准版的，可以说随着标准版模式的升级，挪威版也相应做出了调整。3.1模式可行性研究情况
在NRL(Naval Rearch Laboratory)2009年的回顾报告中指出自从2006年基于H Y C O M 模式的第二代预报系统投入到运营后其预报能力通过了观测数据的测试证明该模式有效[9]。3.2海气耦合应用
海气相互作用的深入研究使得大气与海洋模式、大气与海浪模式、甚至是大气、海浪以及大洋模式的耦合得到前所未有的发展。早期的海流模式选用P O M ，但是H Y C O M 模式自身的优势以及P O M 模式本身发展的受限使得耦合模式的研究侧重点转移到了日渐成熟的H Y C O M 模式上面,Hendrik L,Tolman [10]指出HYCOM模式代替传统的P O M 模式加入到与H W R F 、W W 3模式的耦合是完全可行的，而且在2011年他的团队提高了模式的并行效率;Jo h n Michalakes [11]对目前可用的耦合器ESMF,MC EL,MCT ,HR F
Cou ple r在WR F上的使
用效率进行了对比，并详细给出了利用M C E L 耦合器将H Y C O M 与W R F 进行耦合的整套流程。3.3国内外研究进展
国内外学者及相关研究机构利用模式研究了深海大洋、沿海和陆架海区等不同区域的水文特征并做了很多前瞻性的工作。George Halliwell [12～14]在2002年率先利用H Y C O M 模式对美国西福罗里达州海岸的嵌套技术应用、垂向坐标选取及混合算法的应用进行了初步研究，并指出采用嵌套技术可以明显提高模拟能力，垂向坐标的选取及混合算法的选取对近岸区域的流场会产生一定的影响以及采用
底部边界层参
数化的K P P 方案的有效性;G e o r g e
R.
Halliwell [15]在2004年又对HYCOM模式的垂直混合层算法及垂向坐标进行了评估，得到不同的方案在气候态中存在较大的差异，其中非局地KPP参数化方案、GISS湍封闭方案和MY 湍封闭方案的模拟结果较好，这为后来的研究奠定了基础;Hogan [16]在2006年针对已开展的H Y C O M 模式全球及近岸区域的嵌套试验指出了高分辨率嵌套试验能够提供较好的边界条件，准确的反应近岸的流场;Cécile Renaudie [17]在2009年利用H Y C O M 模式就混合层的垂直分辨率的适应性进行了分析，指出通过对混合层的深度进行适当的定义去控制模式垂向坐标的尺度因子及分层虽然能够给出较为合理的垂向分辨率设置提高了模拟效果，但也指出该方法受混合层深度的定义、采用的强迫场的衰减率及边界条件因素影响较大;Nina Gjerde Winther [18]等人利用HYCOM 模式在北海及丹麦的斯卡格拉克海峡进行了近岸垂向坐标敏感性试验，试验结果表明该模式不但可以用于准确描述深海大洋及海盆尺度的模拟，而且也可以用于沿海和陆架海区的模拟;ERIC P.CHASSIGNET,L I N D A T.S M I T H ,A N D G E O R G E R.H A L L IW E L L [19]利用H Y C O M 模式在大西洋区域分别探讨了垂向坐标、参考气压及温压效应对模拟结果的影响，结果指出HYC OM模式坐标选取具有相当的灵活性，研究的深度不同选用的
参考气压不同以及模式中考虑温压效应会有效减小近表层及近底层的气压梯度力误差。Pa tr ic k J.H o g a n [20]在N R L 利用N C O M 模式与H Y C O M 模式就近岸区域进行了高分辨率数值比较试验，指出了H YC OM 模式的趋势是发展它的近岸模拟能力；白志鹏，高松等[21]利用H Y C O M 进行三重嵌套对黑潮的气候态进行了模拟，得出了高水平分辨率模式对黑潮特征量的模拟会有明显的改进；孙双文[22]等在2009年利用H Y C O M 模式对黑潮对中国近海热盐输运的季节和年际变化进行了研究并指出明，黑潮对东海的热量输运与该区域的海表面热通量是同量级的，且呈反向变化，热量输运冬季最大，夏季最小；朱江、谢基平[23]在2012年利用HYCOM模式的2.1及2.2版本对中国海区进行了以嵌套技术为基础的敏感性试验，明确指出了模拟中国海近岸及陆架海结构所需的条件为准确的风场、全球及局地的高分辨率地形数据、能够准确刻画出近岸射流与涡旋的高水平分辨率(1～3km)、包含潮汐过程的
.com. All Rights Rerved.
大洋模式。Metzger[24]等给南海模式的使用提出了一个难题，他们利用N L O M模式对南海进行一系列模拟后认为，黑潮在吕宋海峡的路径模拟对吕宋海峡内的小岛屿和浅滩的表征极其敏感，水平分辨率至少需要达到1/16°才能保证模拟的黑潮平均路径和观测一致，然而整个南海区域均达到1/16°分辨率对计算资源的要求仍然过高，后经卢著敏，尚晓东，陈桂英[25]采用混合坐标模式(HYCOM)模拟南海环流得到即使0. 2°的分辨率,E R A-15气候数据驱动的H Y C O M模式对南海环流已经能够获得相当
合理的结果，为未来建立基于H YC OM模式的南海预报系统奠定了坚实的基础。
邀请函英文
4 结论与展望
本文基于美国海军研发的H Y C O M模式以及挪威南森环境遥感中心改进的N E R S C-H Y C O M模式，将国内外学者及相关研究机构利用该模式就深海大洋、沿海和陆架海区等不同区域所做的很多前瞻性的工作进行了总结。可以看到HY CO M模式由于采用了三个坐标的混合，因此无论是作为大洋模式还是区域模式都有着较多的应用，更为重要的是国外的工作指出了H Y C O M模式近岸模拟的可行性。
目前国内使用该模式的学者也越来越多，但国内的侧重点主要是大洋及深度较深的水域，对于近岸的模拟工作开展相对很少，而国外早已实现H YC OM模式的业务化运行并取得了很好的效果，为了弥补这一工作的空白应该利用H Y C O M模式发展国内的近岸模式从而提高对近岸海域温、盐、流场的业务化保障能力。
参考文献
[1]Chassignet E.P.,H.E. Hurlburt,O.
M. Smedstad ,et al.Ocean Prediction
with the HYbrid Coordinate Ocean
Model (HYCOM)[M].In“Ocean Weather
Forecasting:An Integrated View of
Oceanography”,E.P. Chassignet and
J.Verron(Eds.),Springer,2006:413-
426.
[2]郑沛楠，朱军，张芳苒，等.常用海洋数
值模式简介[J].海洋预报，2008,25(4):
108-120.
[3]郑沛楠.黑潮对日本海边界环流的影响
[D].中国海洋大学,2009.
[4]郑沛楠，吴德馨，陈学恩，等.基于
HY CO M的风生大洋环流模拟及季节变
化分析[J].中国海洋大学学报：自然科
学版,2009,39(1):7-12,114.
[5]NERSC-HYCOM 2.2[EB/OL].Avail-
able sc.no/hycom.
[6]于庆龙，王辉，万莉颖.南海海表温度时
空分布特征的数值模拟[J].海洋预报,
2010,27(4):59-66.
[7]Hycom ur’s manual[EB/OL].Avail-
able at panoramix.rsmas.
miami. edu/hycom/documentation.
html.
[8]Rainer Bleck.An oceanic general cir-
culation model framed in hybrid isopy-
cnic-Cartesian coordinates[J].Ocean
Modelling,2002(37):55-88.
[9] E.J.Metzg er,H.E.Hurl burt,A.J.
Wallcraft.Predicting “Ocean Weather”
Using the HYbrid Coordinate Ocean
Model (HYCOM)[J].Ocean science and
technology,2009.
[10]H en d ri k L.T ol m an．H YC O M a n d
WAVEWATCH III for advancing hur-
ricane Guidance[R]．Marine Modeling and
Analysis Branch,2010,2.
[11]John Michalakes.Coupling WRF to
Other Models[R].NREL,2010,6.
[12]George Halliwell.Initial Progress on H
YCOM Nested West Florida Shelf Sim
ulations[EB/OL].Available at
hycom/attachments/086_HYCOM
-Halliwell.pdf.
[13]George Halliwell. Sensitivity of Wes
t Florida Shelf Simulations to Initial a
nd Boundary Conditions Provided by
HYCOM Data-Assimilative Ocean Hi
ndcasts[EB/OL].Available at
hycom/attachments/086_halliwell.
pdf.
[14]George Halliwell Evaluation of
HYCOM Performance on the West
Florida Shelf[EB/OL].Available at
ht tp://hy com/att ach ment s/
086_5_Halliwell.pdf.
[15]George R.Halliwell.Evaluation of ver-
tical coordinate and vertical mixing
algorithms in the Hybrid-Coordinate
Ocean[J].Ocean Modelling,2004,7:285-
322.
[16]Hogan.An Overview of Nested Re-
gions Using HYCOM[EB/OL].Avail-
able at hycom/attach-
ments/086_Hogan.pdf.
[17]Cécile Renaudie.Adaptation of the
vertical resolution in the mixed layer
for HYCOM[J].Ocean Modelling,2009,
3(2-3):179-189.
[18]Nina Gjerde Winther,Geir Evenn.A
Hybrid Coordinate Ocean Model for
shelf a simulation[J],2006(1):221-237.
[19]Chassignet,E.P.,L.T.Smith,G.R.
英语诗歌朗诵Halliwell,and R.Bleck,North Atlantic
simulation with the HYbrid Coordi-
nate Ocean Model(HYCOM):Impact of
the vertical coordinate choice[J].ref-
erence density and thermobaricity. J.
Phys.Oceanogr.2003(33),2504-2526.
[20]Patrick J.Hogan,Nesting studies with
hycom at NRL[EB/OL].Available at
ht tp://hy com/att ach ment s/
086_HYCOM-Hogan_1.pdf[EB/OL].
[21]白志鹏，高松，王海棠.H Y C O M模式对
东海黑潮的气候态模拟[J].海洋通报,
2010,29(2):121-129.
[22]孙双文，万彪.黑潮对中国近海热盐输
运的季节和年际变化[J].华东师范大学
学报：自然科学版,2009(3):98-107.
[23]Jiang Zhu,Jiping Xie.Modelling and
data assimilation in marginal as around
China[EB/OL].Available at
w w w.g o d a e.o r g/.../2.9-
Zhu-Jiang-2012-Miami.pdf.
[24]M ET Z G ER E J,H UR L B UR T H E.
The nondeterministic nature of Kuroshio
penetration and eddy shedding in the
South China Sea[J].J Phys Oceanogr,
2001(31):1712-1732.
[25]卢著敏，尚晓东，陈桂英.混合坐标模式
H Y C O M模拟C O A D S强迫下的南海平
均环流[J].热带海洋学报,2008,27(4).
23-31.
.com. All Rights Rerved.
