基于HYCOM的风生大洋环流模拟及季节变化分析

基于垂向混合坐标系统的海洋模式HYCOM建立了全球大洋气候态环流场.在此基础上与前人研究工作进行对比,分析和讨论了全球风生大洋环流场的季节变化情况.从模拟结果看该模式具有较好的模拟能力,可以合理地模拟南极绕极流、赤道流系、黑潮和湾流等世界各主要大洋流系.从断面温度场、流函数分布和断面流量场等分析显示:南极绕极流堪称世界最强流,湾流整体强于黑潮,3者都具有夏季增强、冬季减弱的特点.HYCOM模式在国外的研究方兴未艾,而在国内的应用尚处在起步阶段.本文通过对该模式的介绍和结果分析,向读者推荐使用该模式.     

一种改进的HYCOM模式垂向坐标分层方法和模拟试验

针对22层的HYCOM海洋模式的不足,本文根据北太平洋位势密度的时空分布特征,建立了一种28层的垂向分层方法,并比较了28层HYCOM模式与22层HYCOM模式在积分过程中能量的变化,以及表层流场和赤道断面流系的模拟效果。结果显示:在模式积分过程中,本文采用的28层垂向坐标模式比传统22层模式对强迫场的响应更快,对动能的摄入更大,并且更加容易达到稳定状态;对于流场的刻画,28层模式对赤道表层流系和赤道潜流相比22层模式有明显改进。     

马尔代夫海域水位、海流数值预报研究

基于无结构有限体积法海洋模式(FVCOM),建立了马尔代夫双重嵌套的水位、海流预报模式,并实现了业务化运行。利用三角网格提高重点区域(马尔代夫大桥及岛屿附近海域)的分辨率,最高网格分辨率达到45 m。垂向分层采用σ-s混合坐标的方式划分,分为31层,分别在表层和底层进行加密。采用GFS预报的风场、气压场和热通量结果制作模式表面强迫场文件。在开边界处与HYCOM预报结果进行嵌套,在斜压条件下,采用热启动的方式,业务化模拟了马尔代夫海域2020年的水位流场过程。结果表明,模式能够较好地再现计算海域内天文潮和综合水位的预报,模式预报的水位值与潮位站实测值非常接近。     

浪致混合对亚热带冬季海洋混合强度的影响

上层海洋在全球气候系统中起着至关重要的作用。对上层海洋层结及混合的模拟偏差一直是海洋和气候数值模式发展中悬而未决的问题。本文首先评估了CMIP5中45个模式对上层海洋层结模拟的偏差,确认了冬季亚热带地区海洋模式垂向混合偏强。随后,基于自然资源部第一海洋研究所地球系统模式(FIO-ESM v1.0),分别开展了1986-2005年期间包含和关闭海浪垂向混合情况下的数值实验,分析浪致混合对亚热带冬季海洋混合强度模拟的影响及机制。发现浪致混合使得气候模式中亚热带海域冬季的海洋层结增强,增强的层结使上层海洋更加稳定。首次揭示了增加浪致混合反而降低了海洋总体的垂向混合率:浪致混合使北半球冬季亚热带海域混合率从无浪实验的227 cm²/s降低到有浪实验的178 cm²/s,降低了21.6%;南半球冬季亚热带海域混合率从无浪实验的189 cm²/s降低到有浪实验的165 cm²/s,降低了12.7%。进一步分析发现,浪致混合主要是通过增加冬季亚热带海域上层海洋的热含量从而强化了海洋的层结,最终改善了气候模式对上层...     

台风“巴蓬”对南海上层海洋水文特征与生态环境的影响

研究台风过境影响下海洋生态环境的变化是系统认识强天气过程下海洋响应问题的重要内容,对于防灾减灾、远洋运输等具有积极的意义。基于混合坐标海洋环流模式（Hybrid Coordinate Ocean Model, HYCOM）数据,结合Aqua卫星搭载的中分辨率成像光谱仪（MODIS）传感器数据,研究了慢速移动的台风“巴蓬”（1929号）过境对南海海表温度、叶绿素浓度以及水体垂向温盐结构和流场特征的影响。结果显示,台风“巴蓬”过境期间,波高由3 m增加至9 m,海表流速由0.3 m/s增加至0.7 m/s,台风周围最大流速甚至超过1.8 m/s;台风过境后,海表温度由原25.8℃下降至24.9℃。在台风影响下,水体垂向混合增加,混合层深度增加近15 m,从23 m到38 m。混合层内上部温度降低、盐度略微增加,下部温度升高、盐度降低;上层水体流速增加,方向偏转近90°,流速变化幅度由表及底逐渐减小。台风过境后3天内,由于温度降低、光照减弱,并且垂向混合增强补充的营养盐仍未被浮游植物充分利用,因此,叶绿素浓度并未快速升高,而是呈现短期内浓度微增的现象。     

夏季黄河入海径流对黄河口及附近海域环流影响的数值研究

基于区域海洋模式ROMS,建立了一个三维非线性斜压浅海模式,考虑了包括径流、风场、海面热交换以及黄渤海环流等因素,研究了夏季8月份黄河入海径流量对黄河口及附近海域环流结构的影响。数值实验较好地佐证了黄河冲淡水的"北偏"现象,并很好的体现了冲淡水对河口附近海域环流结构的影响。数值研究表明:黄河入海径流量对河口附近海域环流结构有显著影响,径流越大冲淡水向北-西北方向偏转越明显,同时流轴中心余流流速也显著增大;莱州湾顺时针环流受黄河入海径流影响显著,径流量越大越不利于该环流的发育和维持,而径流量越小环流越稳定;径流量越大导致河口附近海域表层余流加大,余流垂向梯度得到加强,底部补偿流增强,河口垂向环流越明显。     

黄海、东海夏季环流的数值模拟

基于ECOM -si,建立了一个垂向取σ坐标、水平取球面坐标的三维高分辨率湍流闭合的黄海、东海海洋环流数值模式 .考虑实际岸线和水深、边界通量、密度梯度力和风应力 ,应用该模式较成功地模拟出了夏季黄海、东海环流 (黑潮、台湾暖流、对马暖流、黄海冷水团环流和沿岸流 ) ,摸拟结果与以往实际观测和研究结果较为一致 ,表明所建立的模式能较成功地应用于该海区环流的研究     

海洋模式中潮致混合调和分析参数化方案的建立

潮致混合对海洋环流的调整起着重要作用。陆架环流的数值模拟中如果不考虑潮汐作用,往往不能得到与观测相符的垂向温盐结构。本文基于调和分析方法,建立了一套潮致混合参数化方案。该方案通过对垂向混合系数进行调和分析,从而得到随时间变化的潮致混合系数。将该方案用于黄海冷水团数值模拟的结果显示,其能够得到与在数值模式开边界直接加入潮汐强迫相当的冷水团温盐结构。和直接引入潮汐强迫相比,这一潮致混合参数化方案的优势在于,它能够大大节省数值模拟计算机时,因此有望显著提高大规模高分辨率的海洋环流及气候模式的模拟能力和计算效率。     

S坐标下西北太平洋环流和海温数值模拟

综合考虑侧边界水通量、海表面风应力、热通量、蒸发和降水等,建立了一个采用垂向S坐标的西北太平洋海洋模式,模拟和分析西北太平洋环流和海温。S坐标具有使海洋上层趋于同一水深的物理平面和同时底层拟合海底变化的优点。与日本2008年西北太平洋实测温度断面资料比较,模式计算的海温在分布势态和量值上较为一致,模式能较好地验证海洋温度分布情况。与σ坐标相比,S坐标能更好地验证海温的变化。模式较好地模拟出了西北太平洋气候态环流和海温变化。由于在σ坐标系下水深的水平变化体现在σ层上,海洋表层的分辨率比S坐标系下低,以及初始海温的水平分布含有水深水平变化的信息,S坐标海洋模式能更好地模拟环流和海温的变化。     

海洋模式的水平分辨率对吕宋海峡深层环流数值模拟结果的影响

基于南海观测得到的垂向混合率修改KPP垂向混合方案,利用大洋环流模式HYCOM首次模拟得到了吕宋海峡深层环流的空间分布特征,并通过3个不同水平分辨率的实验(1/6(°)、1/12(°)、1/24(°))讨论水平分辨率对模拟吕宋海峡深层环流空间结构的影响。模拟结果表明:(1)1/6(°)水平分辨率过于粗糙,无法分辨巴士海峡和台东海峡地形特征,无法得到吕宋海峡深层环流的空间分布特征;(2)1/12(°)的水平分辨率可以很好的分辨吕宋海峡和台东海峡的地形特征,模拟得到吕宋海峡深层环流的空间分布和流场特征;(3)1/24(°)水平分辨率的模拟结果表明,更高的水平分辨率不会改变吕宋海峡深层环流的空间分布和主要出入口的垂向结构,只是会显示更细节的环流结构。1/12(°)和1/24(°)水平分辨率的模式结果都表明,西北太平洋深层水通过巴士海峡和台东海峡进入吕宋海沟,年平均流量分别为1.1和0.4Sv,然后沿吕宋海沟向南海方向流动,最后主要通过位于恒春海脊上的2个缺口进入南海,年平均流量分别为0.5和0.9Sv。     

A hybrid coordinate ocean model for shelf sea simulation

The general circulation in the North Sea and Skagerrak is simulated using the hybrid coordinate ocean model (HYCOM). Although HYCOM was originally developed for simulations of the open ocean, it has a design which should make it applicable also for coastal and shallow shelf seas. Thus, the objective of this study has been to examine the skills of the present version of HYCOM in a coastal shelf application, and to identify the areas where HYCOM needs to be further developed. To demonstrate the capability of the vertical coordinate in HYCOM, three experiments with different configurations of the vertical coordinate were carried out. In general, the results from these experiments compares quite well with in situ and satellite data, and the water masses and the general circulation in the North Sea and Skagerrak is reproduced in the simulations. Differences between the three experiments are small compared to other errors, which are related to a combined effect of model setup and properties of the vertical mixing scheme. Hence, it is difficult to quantify which vertical coordinate configuration works best for the coastal region. It is concluded that HYCOM can be used for simulations of coastal and shelf seas, and further suggestions for improving the model results are given. Since HYCOM also works well in open ocean and basin scale simulations, it may allow for a realistic modelling of the transition region between the open ocean and coastal shelf seas.     

Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

Comparison of gravity current mixing parameterizations and calibration using a high-resolution 3D nonhydrostatic spectral element model

In light of the pressing need for development and testing of reliable parameterizations of gravity current entrainment in ocean general circulation models, two existing entrainment parameterization schemes, K-profile parameterization (KPP) and one based on Turner’s work (TP), are compared using idealized experiments of dense water flow over a constant-slope wedge using the HYbrid Coordinate Ocean Model (HYCOM). It is found that the gravity current entrainment resulting from KPP and TP differ significantly from one another. Parameters of KPP and TP are then calibrated using results from the high-order nonhydrostatic spectral element model Nek5000. It is shown that a very good agreement can be reached between the HYCOM simulations with KPP and TP, even though these schemes are quite different from each other.     

基于ROMS模式的南海SST与SSH四维变分同化研究

卫星遥感观测获得了大量高分辨率的海面实时信息,包括海面温度（SST）和海面高度（SSH）等,同化进入数值模式可有效提升模拟精度。本文基于ROMS模式与四维变分同化方法（4DVAR）,使用AVHRR SST和AVISO SSH数据,开展了南海区域同化实验。为检验同化的效果,分别利用HYCOM再分析资料和Argo温盐实测数据分析了同化结果的海面高度、流场及温盐剖面的精度。对比结果表明,SST和SSH的同化能够改善ROMS的模拟结果:同化后海面高度场能够更为准确地捕捉海洋的中尺度特征,与HYCOM海面高度再分析资料相比,平均绝对偏差和均方根误差分别为0.054 m和0.066 m;与HYCOM 10 m层流场相比,东向与北向流速平均绝对偏差分别为0.12 m/s和0.11 m/s,相比未同化均提升约0.01 m/s;温盐同化结果与Argo温盐实测具有较高的一致性,温度和盐度平均绝对偏差为0.45℃、0.077,均方根误差为0.91℃、0.11,单个的温盐廓线对比说明,同化结果与HYCOM再分析资料精度相当。     

Impact of Ocean Mean Dynamic Topography on Satellite Data Assimilation

The response of an eddy-permitting ocean model to changes imposed by the use of different mean dynamic topographies (MDT) is analyzed in a multivariate assimilation context, allowing the evaluation of this impact, not only on the surface circulation, but also on the interior ocean representation. The assimilation scheme is a reduced-order sequential Kalman filter (SEEK). In a first set of experiments, high resolution sea surface temperature, along-track sea surface height and sea surface salinity from climatology are assimilated into a 1/3° resolution North and Tropical Atlantic version of the HYCOM model. In a second experiment, in situ profile data are assimilated in addition to the surface measurements.

The first set of experiments illustrates that important differences in the representation of the horizontal model circulation pattern are related to differences in the MDT used. The objective of assimilation is to improve the representation of the 3D ocean state. However, the imperfect representation of the mean dynamic topography appears to be an important limiting factor with regard to the degree of realism obtained in the simulated flow.

Vertical temperature and salinity profiles are key observations to drive a general circulation ocean model toward a more realistic state. The second set of experiments shows that assimilating them in addition to sea surface measurements is a far from trivial exercise. A specific difficulty is due to inconsistencies between the dynamic topography diagnosed from in situ observations and that diagnosed from sea surface height. These two fields obtained from different data sources do not contain exactly the same information. In order to overcome this difficulty, a strategy is proposed and validated.     

Impact of Ocean Mean Dynamic Topography on Satellite Data Assimilation

The response of an eddy-permitting ocean model to changes imposed by the use of different mean dynamic topographies (MDT) is analyzed in a multivariate assimilation context, allowing the evaluation of this impact, not only on the surface circulation, but also on the interior ocean representation. The assimilation scheme is a reduced-order sequential Kalman filter (SEEK). In a first set of experiments, high resolution sea surface temperature, along-track sea surface height and sea surface salinity from climatology are assimilated into a 1/3° resolution North and Tropical Atlantic version of the HYCOM model. In a second experiment, in situ profile data are assimilated in addition to the surface measurements.

The first set of experiments illustrates that important differences in the representation of the horizontal model circulation pattern are related to differences in the MDT used. The objective of assimilation is to improve the representation of the 3D ocean state. However, the imperfect representation of the mean dynamic topography appears to be an important limiting factor with regard to the degree of realism obtained in the simulated flow.

Vertical temperature and salinity profiles are key observations to drive a general circulation ocean model toward a more realistic state. The second set of experiments shows that assimilating them in addition to sea surface measurements is a far from trivial exercise. A specific difficulty is due to inconsistencies between the dynamic topography diagnosed from in situ observations and that diagnosed from sea surface height. These two fields obtained from different data sources do not contain exactly the same information. In order to overcome this difficulty, a strategy is proposed and validated.     

热带西太平洋潜流模拟:(Ⅰ)模式配置与模拟结果验证

利用高分辨率准全球海洋模式HYCOM(HYbrid Coordinate Ocean Model),采用局部加密网格的方法,模式积分30年,模拟了热带西太平洋环流的年平均状态和季节变化.通过与实测的海表面高度、测流资料、温盐资料等对比分析,验证了模拟结果的有效性,认为模式结果可以用于研究热带西太平洋潜流及其季节变化.     

Can an atmospherically forced ocean model accurately simulate sea surface temperature during ENSO events?

The performance of an atmospherically forced ocean general circulation model (OGCM) in simulating daily and monthly sea surface temperature (SST) is examined during the historical El Niño Southern Oscillation (ENSO) events during the time period 1993–2003. For this purpose, we use the HYbrid Coordinate Ocean Model (HYCOM) configured for the North Pacific north of 20°S at a resolution of ≈9 km. There is no assimilation of (or relaxation to) SST data and no date-specific assimilation of any data type. The ability of the model in simulating temporal variations of SST anomalies is discussed by comparing model results with two satellite-based SST products. The HYCOM simulation gives a basin-averaged monthly mean bias of 0.3 °C and rms difference of 0.6 °C over the North Pacific Ocean during 1993–2003. While the model is able to simulate SST anomalies with mean biases  <0.5 °C  in comparison to observations during most of the ENSO events, limitations in the accuracy of atmospheric forcing (specifically, net short-wave radiation) have some influence on the accuracy of simulations. This is specifically demonstrated during the 1998 transition period from El Niño to La Niña, when a record large SST drop of  ≈7 °C  occurred in the eastern equatorial Pacific Ocean.     

涡分辨率的南海区域海洋环流模式构建技术研究

针对复杂海洋中尺度现象模拟仿真的技术难题,利用区域海洋环流模式ROMS,通过模式的网格构建、地形处理以及模式的初始场和强迫场处理技术,构建出一套涡分辨率南海区域海洋环流模式。通过模式模拟结果与卫星遥感实测资料等对比,发现该模式能够较好地模拟出南海涡旋及其引发的海温异常等海洋中小尺度过程,说明该模式可作为研究复杂海洋中尺度现象影响海军武器装备效能的环境数值仿真手段。