表面波浪对海洋上层混合的作用

A new three-dimensional numerical model is derived through a wave average on the primitive N-S equations, in which both the＂Coriolis-Stokes forcing＂ and the＂Stokes-Vortex force＂ are considered. Three ideal experiments are run using the new model applied to the Princeton ocean model （POM）. Numerical results show that surface waves play an important role on the mixing of the upper ocean. The mixed layer is enhanced when wave effect is considered in conjunction with small Langmuir numbers. Both surface wave breaking and Stokes production can strengthen the turbulent mixing near the surface. However, the influence of wave breaking is limited to a thin layer, but Stokes drift can affect the whole mixed layer. Furthermore, the vertical mixing coefficients clearly rise in the mixed layer, and the upper ocean mixed layer is deepened especially in the Antarctic Circumpolar Current when the model is applied to global simulations. It indicates that the surface gravity waves are indispensable in enhancing the mixing in the upper ocean, and should be accounted for in ocean general circulation models.     

热力学效应及斯托克斯漂流对混合层影响的研究

本文通过理想化的外部强迫以及海洋站点实测数据驱动普林斯顿海洋模式来研究海洋热力学效应和斯托克斯漂流对上混合层数值模拟的影响。在Mellor-Yamada湍流闭合方案中，经常出现夏季海表面温度偏暖和混合层深度偏浅的模拟误差。实验表明，斯托克斯漂流在冬季和夏季均能增强湍流动能，加深混合层深度。这种效应可以改善夏季的模拟结果，但与观测数据相比，将增大冬季混合层深度的模拟误差。斯托克斯漂流可以通过增强湍动能来加深混合层深度。结果表明，将斯托克斯漂流与冷皮层和暖层对上部混合层的热效应相结合，可以正确地模拟混合层深度。在夏季，海洋冷皮层和暖层通过“阻挡结构”和双温跃层结构模拟出更真实的上混合层变化。在冬季，海洋热力学效应通过增强上层海洋层结平衡了斯托克斯漂流的影响，并且由斯托克斯漂流引起的过度混合被校正。     

全球Stokes 漂流的时空分布特征研究

利用海浪模式WWIII(Wave Watch III)2008年的模拟结果对海面Stokes漂流、Stokes输运、Stokes深度以及全球Langmuir数的年平均分布特征和季节平均分布特征分别进行了详细的研究与分析。结果表明,海面Stokes漂流和Stokes输运均呈现高纬度偏大的特征,以南极绕极流海域最为突出。全球大部分海域Stokes漂流影响深度在20 m以内,呈现大洋东部偏大,西部偏小的分布特征。全球大部分海域的混合作用是剪切不稳定性和Langmuir湍效应并存的状态,甚至有些海域是以Langmuir湍效应为主。因此,在进行大尺度的海洋数值模拟时,应该考虑波浪导致的混合效应。     

台风条件下朗缪尔环流对上层海洋混合的影响研究进展

回顾了近10年来台风条件下朗缪尔环流影响上层海洋混合的研究进展，朗缪尔致湍流对海洋上混合层的形成和加深的重要作用已形成了基本共识，但对于朗缪尔致湍流对海洋上混合层的混合作用机制和程度仍然存在诸多不确定性。观测表明台风条件下台风眼附近的混合层平均湍流动能受到了较强的抑制，可能与台风不同位置朗缪尔致湍流的特征变异有关；台风条件下，现有的朗缪尔致湍流参数化方案在上层混合过程模拟中还有显著误差。在今后研究中，通过改进斯托克斯漂流剖面的计算方法，优化表征台风条件下海面状况的朗缪尔致湍流参数化计算方案，是进一步揭示台风条件下朗缪尔环流对海洋上层混合的影响机理的必要途径。     

Effect of Langmuir circulation on upper ocean mixing in the South China Sea

Effect of Langmuir circulation (LC) on upper ocean mixing is investigated by a two-way wave-current coupled model. Themodel is coupled of the ocean circulationmodel ROMS (regional ocean modeling system) to the surface wave model SWAN (simulating waves nearshore) via the model-coupling toolkit. The LC already certified its importance by many one-dimensional (1D) research andmechanismanalysis work. This work focuses on inducing LC’s effect in a three-dimensional (3-D) model and applying it to real field modeling. In ROMS, theMellor-Yamada turbulence closuremixing scheme is modified by including LC’s effect. The SWAN imports bathymetry, free surface and current information fromthe ROMS while exports significant wave parameters to the ROMS for Stokes wave computing every 6 s. This coupled model is applied to the South China Sea (SCS) during September 2008 cruise. The results show that LC increasing turbulence and deepening mixed layer depth (MLD) at order of O (10 m) in most of the areas, especially in the north part of SCS where most of our measurements operated. The coupled model further includes wave breaking which will bringsmore energy into water. When LC works together with wave breaking,more energy is transferred into deep layer and accelerates the MLD deepening. In the north part of the SCS, their effects aremore obvious. This is consistent with big wind event in the area of the Zhujiang River Delta. The shallow water depth as another reasonmakes themeasy to influence the oceanmixing as well.     

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.     

在大洋环流模型中对应于惯性内波破碎的垂直混合的一个参数化方案

Based on the theoretical spectral model of inertial internal wave breaking(fine structure) proposed previously, in which the effects of the horizontal Coriolis frequency component f-tilde on a potential isopycnal are taken into account, a parameterization scheme of vertical mixing in the stably stratified interior below the surface mixed layer in the ocean general circulation model(OGCM) is put forward preliminarily in this paper. Besides turbulence, the impact of sub-mesoscale oceanic processes(including inertial internal wave breaking product) on oceanic interior mixing is emphasized. We suggest that adding the inertial internal wave breaking mixing scheme(F-scheme for short) put forward in this paper to the turbulence mixing scheme of Canuto et al.( T-scheme for short) in the OGCM, except the region from 15°S to 15°N. The numerical results of F-scheme by using WOA09 data and an OGCM(LICOM, LASG/IAP climate system ocean model) over the global ocean are given. A notable improvement in the simulation of salinity and temperature over the global ocean is attained by using T-scheme adding F-scheme, especially in the mid- and high-latitude regions in the simulation of the intermediate water and deep water. We conjecture that the inertial internal wave breaking mixing and inertial forcing of wind might be one of important mechanisms maintaining the ventilation process. The modeling strength of the Atlantic meridional overturning circulation(AMOC) by using T-scheme adding F-scheme may be more reasonable than that by using T-scheme alone, though the physical processes need to be further studied, and the overflow parameterization needs to be incorporated. A shortcoming in F-scheme is that in this paper the error of simulated salinity and temperature by using T-scheme adding F-scheme is larger than that by using T-scheme alone in the subsurface layer.     

波浪破碎下湍流混合的实验研究

波浪破碎过程产生的湍流动量和能量垂向输运对于加快海洋上混合层中垂向混合具有显著效果。采用二维实验室水槽中对波浪破碎过程进行模拟。对采集的波浪振幅时间序列采用希尔伯特变换定位破碎波位置,波浪的破碎率随有效波高的增加而增大,波浪谱分析得到的波浪基本周期与有效周期结果相似。实验中采用粒子图像测速技术(particle image velocimetry, PIV)计算波浪破碎过程中湍动能耗散率的空间分布。湍流强度与波浪的相位密切相关,波峰位置处湍流活动最为剧烈,而且波峰位置处湍流混合区内湍动能耗散率量值的垂向分布基本保持不变,即出现"湍流饱和"现象,湍流影响深度可以达到波高的70%—90%。计算湍流扩散系数的垂向分布发现,湍流扩散在混合区上部随深度的增大以指数函数的形式增加,在混合区下部趋于稳定。作为对比,在相同位置处对声学多普勒流速测量仪(acoustic Doppler velocimeter, ADV)测量的单点流速做频谱分析,发现与该位置处PIV湍动能耗散率结果量级处于同一水平,进一步验证了实验结果的准确性。     

Langmuir环流在雷诺平均型模式中的参数化研究进展

Langmuir环流影响着海洋上层的能量输入,对海洋上混合层的形成和加深起着重要作用,对于海洋上混合层具有重要意义。近年来许多学者采用大涡模拟(LES)方法对Langmuir环流进行机制研究,并通过在雷诺平均模型中参数化Langmuir环流效应,将Langmuir环流过程引入到三维海洋环流或海洋耦合模式中,提出了一系列混合参数化方案。本文回顾了Langmuir环流在雷诺平均模式参数化中的研究进展,主要可分为以下几种方案:一种方法是用Langmuir数在KPP垂直混合参数化方案中引入湍流特征速度增强因子,并不断发展Langmuir数的定义;一种是在Mellor-Yamada2.5湍流闭合模型中增加斯托克斯漂流剪切效应项,此外还有通过修改模式中混合长方程来加入Langmuir效应等。通过在雷诺平均模式中应用的结果来看,现有的参数化方案在一定程度上改善了混合层深度和SST的模拟,肯定了Langmuir环流在加深混合等方面的作用,但仍存在一些问题需要在今后的研究中进一步改进。     

A global wave parameter database for geophysical applications. Part 1: Wave-current–turbulence interaction parameters for the open ocean based on traditional parameterizations

Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift U_ss, volume Stokes transport T_S, a wave breaking height scale H_swg, and the flux of energy from waves to ocean turbulence Φ_oc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φ_oc is 68 TW, and the yearly mean value of T_S is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of U_ss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, U_ss ≃ 0.013U₁₀, where U₁₀ is the wind speed at 10 m height.     

Effect of surface wave breaking on the surface boundary layer of temperature in the Yellow Sea in summer

Near-surface enhancement of turbulent mixing and vertical mixing coefficient for temperature owing to the effect of surface wave breaking is investigated using a two-dimensional (2-D) ocean circulation model with a tidal boundary condition in an idealized shelf sea. On the basis of the 2-D simulation, the effect of surface wave breaking on surface boundary layer deepening in the Yellow Sea in summer is studied utilizing a 3-D ocean circulation model. A well-mixed temperature surface layer in the Yellow Sea can be successfully reconstructed when the effect of surface wave breaking is considered. The diagnostic analysis of the turbulent kinetic energy equation shows that turbulent mixing is enhanced greatly in the Yellow Sea in summer by surface wave breaking. In addition, the diagnostic analysis of momentum budget and temperature budget also show that surface wave breaking has an evident contribution to the turbulent mixing in the surface boundary layer. We therefore conclude that surface wave breaking is an important factor in determining the depth of the surface boundary layer of temperature in the Yellow Sea in summer.     

随机波浪下斯托克司漂移速度是否引起水面漂移物离散的研究

漂浮于自由水面的污染物的的迁移、扩散会受到天然随机海浪的影响。之前的研究（以Herterich和Hasselmann（1982）为代表）普遍认为，随机波浪作用下的斯托克司漂移速度会引起水面污染物的离散，这个离散甚至有可能跟风和海流引起的离散同一量级。本研究就随机波浪作用下的斯托克司漂移速度是否会引起水面漂移物的离散进行理论和试验探讨。从理论推导可知，随机波浪下的质量输移速度是个定常分量，因此它不会随时间变化而引起水面漂移物的离散。随后我们在实验室水槽中进行了漂移物在随机波浪（P-M谱）作用下的漂移过程的测量。试验结果也印证了随机波浪作用下的斯托克司漂移速度不会引起水面漂移物离散的结论。     

FIO-ESM v2.0 CORE2-forced experiment for the CMIP6 Ocean Model Intercomparison Project

We introduced the Coupled Model Intercomparison Project Phase 6 (CMIP6) Ocean Model Intercomparison Project CORE2-forced (OMIP-1) experiment by using the First Institute of Oceanography Earth System Model version 2.0 (FIO-ESM v2.0), and comprehensively evaluated the simulation results. Unlike other OMIP models, FIO-ESM v2.0 includes a coupled ocean surface wave component model that takes into account non-breaking surface wave-induced vertical mixing in the ocean and effect of surface wave Stokes drift on air-sea momentum and heat fluxes in the climate system. A sub-layer sea surface temperature (SST) diurnal cycle parameterization was also employed to take into account effect of SST diurnal cycle on air-sea heat ?uxes to improve simulations of air-sea interactions. Evaluations show that mean values and long-term trends of significant wave height were adequately reproduced in the FIO-ESM v2.0 OMIP-1 simulations, and there is a reasonable fit between the SST diurnal cycle obtained from in situ observations and that parameterized by FIO-ESM v2.0. Evaluations of model drift, temperature, salinity, mixed layer depth, and the Atlantic Meridional Overturning Circulation show that the model performs well in the FIO-ESM v2.0 OMIP-1 simulation. However, the summer sea ice extent of the Arctic and Antarctic is underestimated.     

Effects of Stokes production on summer ocean shelf dynamics

A two-dimensional numerical model,which is configured on the basis of Princeton ocean model（POM）,is used to study the effect of Stokes production（SP） of the turbulent kinetic energy on a density profile and Ekman transport in an idealized shelf region in summer.The energy input from SP is parameterized and included into the Mellor-Yamada turbulence closure submodel.Results reveal that the intensity of wind-driven upwelling fronts near the sea surface is weakened by the SP-associated turbulent kinetic energy input.The vertical eddy viscosity coefficient in the surface boundary layer is enhanced greatly owing to the impact of SP,which decreases the alongshore velocity and changes the distribution of upwelling.In addition,the SP-induced mixing easily suppresses the strong stratification and significantly increases the depth of the upper mixed layer（ML） under strong winds.     

A global wave parameter database for geophysical applications. Part 1: Wave-current–turbulence interaction parameters for the open ocean based on traditional parameterizations

Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift U_ss, volume Stokes transport T_S, a wave breaking height scale H_swg, and the flux of energy from waves to ocean turbulence Φ_oc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φ_oc is 68 TW, and the yearly mean value of T_S is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of U_ss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, U_ss  0.013U₁₀, where U₁₀ is the wind speed at 10 m height.     

A linear model for the structure of turbulence beneath surface water waves

The structure of turbulence in the ocean surface layer is investigated using a simplified semi-analytical model based on rapid-distortion theory. In this model, which is linear with respect to the turbulence, the flow comprises a mean Eulerian shear current, the Stokes drift of an irrotational surface wave, which accounts for the irreversible effect of the waves on the turbulence, and the turbulence itself, whose time evolution is calculated. By analysing the equations of motion used in the model, which are linearised versions of the Craik–Leibovich equations containing a ‘vortex force’, it is found that a flow including mean shear and a Stokes drift is formally equivalent to a flow including mean shear and rotation. In particular, Craik and Leibovich’s condition for the linear instability of the first kind of flow is equivalent to Bradshaw’s condition for the linear instability of the second. However, the present study goes beyond linear stability analyses by considering flow disturbances of finite amplitude, which allows calculating turbulence statistics and addressing cases where the linear stability is neutral. Results from the model show that the turbulence displays a structure with a continuous variation of the anisotropy and elongation, ranging from streaky structures, for distortion by shear only, to streamwise vortices resembling Langmuir circulations, for distortion by Stokes drift only. The TKE grows faster for distortion by a shear and a Stokes drift gradient with the same sign (a situation relevant to wind waves), but the turbulence is more isotropic in that case (which is linearly unstable to Langmuir circulations).     

一种改进的基于逆Omega-K算法的海面场景SAR原始数据仿真方法

合成孔径雷达(SAR)海面场景原始数据仿真是研究海洋动力参数(表面波浪、风矢量和洋流)的有效工具。目前海面场景原始数据仿真方法已经基于逆Omega-K算法实现了海洋运动参数的空间变化。但是目前仅仅讨论了正侧视情况下的海面场景仿真,应用范围有限,同时没有考虑Stokes漂流以及Bragg相速度的影响,而这两者都是存在于真实海面的。通常情况下为了反演得到海面流场的二维速度矢量,雷达需要从两个不同的方位方向观察海面的同一个区域,因此这就需要考虑大斜视的雷达波束,同时Stokes漂流和Bragg相速度是SAR海表面流场观测不容忽视的两种运动。本文在不改变原有正侧视逆Omega-K算法的情况下,通过增加重新计算零方位时刻的斜视波束中心位置坐标,并据此确定SAR原始数据在多普勒域的位置来将其扩展到大斜侧视逆Omega-K算法,并通过时域Stokes漂流公式到频域内离散化Stokes漂流公式的推导来加入Stokes漂流,以及根据Bragg散射机制加入了Bragg相速度。仿真结果表明,经过聚焦成像后的SAR图像很好的体现了真实海面波浪场的形状,同时能够很好地反演出设定的雷达径向流场速度,且流速精度误差控制在6%以内。最后也证明了Bragg相速度以及Stokes漂流对于海面流场的影响不可忽视。     

近惯性内波破碎混合方案对MOM4模式的影响

风生近惯性内波破碎引起的跨等密度面混合在海洋内部混合中起重要作用。然而其参数化对海洋模式的模拟影响仍有待进一步认识。本文给出的是在模块化海洋模式（MOM）中海洋表面边界层以下引入一个考虑风驱动近惯性内波破碎引起的跨等密度面混合参数化方案的研究工作。模拟结果显示，该方案有效改善MOM4模拟的上层1 000 m以上的温盐偏差，特别是在北太平洋和北大西洋的通风地区。数值试验表明，风生近惯性内波破碎有可能是维持海洋通风过程的重要机制之一，它使得海洋通风区的位温变冷，盐度变淡，整层等位密面加深。维持的通风过程使得北太平洋副极地大涡的影响延伸到副热带大涡。从而模拟的北太平洋中层水源头及其副热带大涡东侧的温盐更接近观测实际。同时，模拟的北大西洋经圈翻转环流强度也更为合理。     

Modeling of the eddy viscosity by breaking waves

Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing,heat and gases exchange across air-sea interface.The energy loss from wave breaking and the bubble plume penetration depth are estimated.As a consequence,the vertical distribution of the turbulent kinetic energy(TKE),the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided.It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface.High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than υwl(=κuwz),the value predicted for the wall layer scaling close to the surface,where uw is the friction velocity in water,κ with 0.4 is the von Kármán constant,and z is the water depth,and the strength of the eddy viscosity depends both on wind speed and sea state,and decays rapidly through the depth.This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer,well above the classical values expected from the similarity theory.Deeper down,however,the effects of wave breaking on the vertical mixing become less important.