风浪场波生横向彻体切应力的分析

使用风浪经验公式和风浪传入近岸的小振幅波动理论计算并分析了地转条件下波生横向彻体切应力随风速、风区、水深等因素的变化规律。结果表明，彻体切应力随风速、风区的增大而增大，直到波浪达到完全成长状态，但增长速度并不规则；为更加直观的显示彻体切应力的海流驱动作用，将其与定常Ekman漂流中的风生湍粘性力进行了有意义的比较，结果表明，在中、高纬度的通常海况下两者具有相同量级，因此，在有波浪场存在的条件下，由风、浪共生海流的总体强度将比理论Ekman漂流大，并且不会完全符合Ekman漂流理论的转向规律，在表层尤其明显。由此说明，在一定情况下在风生漂流研究中考虑由地转引起的波生横向彻体切应力是非常必要的。     

Stokes漂流近似公式对海洋表层流场估算的影响

本文采用波浪订正的Ekman模型,研究分析了三种Stokes漂流近似公式(单波公式、e指数公式、Phillips谱近似公式)对海洋表层流场估算的影响。海表总流场由海表面高度(SSH)数据计算的地转流和海浪模式WAVEWATCH Ⅲ输出结果计算的非地转流组成,并采用拉格朗日浮标观测数据对计算结果进行了验证。研究表明,随着Stokes漂流近似公式精度的提高,其计算的拉格朗日流速更接近于谱积分公式的计算结果,更贴近拉格朗日浮标观测数据。与谱积分公式计算的海表拉格朗日流速相比,单波公式的平均相对偏差为0.0834,e指数公式的平均相对偏差为0.0392,Phillips谱近似公式的平均相对偏差为0.0101,说明Phillips谱近似公式在不同风速下均能对谱积分公式有良好的近似效果。在低风速条件下,由Stokes漂流近似公式精度引起的海洋表层流场估算误差可以忽略不计,但随着风速增加,由近似公式精度引起的偏差逐渐变大,此时应该选择Phillips谱近似公式计算Stokes漂流,来减小误差。     

随机波作用下海床动态响应分析

真实的海洋波浪是随机的,而前人对海床的动态响应分析大都是选用线性波或者Stokes波理论,对海床的模拟大都采用Biot拟静力模型,忽略了流体速度及土体位移加速度的影响。联合使用Longuet-Higgins随机波模型(采用Jonswap谱)以及动力u-p形式的海床响应计算模型,使用COMSOL Multiphysics多场耦合软件的PDE模块输入方程进行有限元计算,得到随机波作用下整体海床动态响应结果。将随机波结果与一阶Stokes波和椭圆余弦波结果进行对比,并对渗透系数和饱和度进行参数分析,研究表明渗透系数和饱和度对于随机波作用下海床动态响应影响显著。     

深海畸形波振幅与增益的关系

依据深海表面波的非线性薛定谔方程,推导深海表面波的增益谱,研究在调制不稳定条件下的非线性薛定谔方程的一阶解与增益之间的关系,又结合呼吸孤子产生畸形波的现象,得到畸形波的最大振幅为初始波振幅与增益最大值的叠加。这个结果对研究畸形波的产生具有一定的参考价值。     

一种改进的基于逆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漂流对于海面流场的影响不可忽视。     

波-流相互作用源函数及其对海浪模拟的影响

当前,基于第三代海浪模式的波-流相互作用研究,通常仅考虑背景流场对海浪群速度大小、方向的改变及其所引起的波浪折射效应,忽略了流场对波浪能谱的影响。本文以波-流相互作用源函数为切入点,利用MASNUM海浪模式进行理想实验,阐释了水平方向上的波-流相互作用对海浪能谱的改变以及对海浪模拟可能造成的影响。结果表明：水平流速梯度量级、水深、波浪成长状态以及流速梯度方向均可以影响波-流相互作用源函数的作用效果,使得有效波高、谱峰波向的模拟结果发生变化。在更符合实际的水平流速梯度量级下,波-流相互作用源函数对有效波高和谱峰波向模拟值的影响非常小;但考虑到涌浪的长距离传播,其对谱峰波向的影响仍可能明显改变涌浪的传播形态。     

台风过程中浪致Stokes漂流对海洋上层温度影响的数值模拟

Stokes漂流对海洋上层温度变化具有重要影响。本文以"麦德姆(Matmo)"台风过程为例,基于浪流耦合模型,通过对比分析考虑和未考虑Stokes漂流的模拟结果,研究了台风过程中浪致Stokes漂流及其效应对海洋上层温度变化的影响。研究表明Stokes漂流及其效应与海浪大小的分布密切相关,海浪越大, Stokes漂流、Stokes输运和Ekman-Stokes数相对越大。Stokes漂流在台风过程中起降低海表面温度的作用,台风路径处的Stokes漂流及其效应较大,降温较明显,最大降温约2°C。产生降温的原因是Stokes漂流造成海表流场改变,以及Stokes输运引起海水辐散等作用加强了上下层海水质量和能量的交换。利用Argo资料进行验证,发现考虑了Stokes漂流作用的海洋上层温度模拟结果与Argo测量结果更接近。     

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

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

深水风浪的风区指数律

对已有根据观测提出的幂函数形式风浪成长关系进行了分析。发现这些风浪成长关系在消去无因次风区后一致地与3／2指数律相协调，尽管它们原来存在较大的不协调性。发现Jeffreys，Sverdrup和Munk以及Platit的风能输入源函数在谱积分意义下具有相似性，而Tsikunov，Hasselmann和Phillips的破波耗散源函数在谱积分意义下也具有相似性，尽管这些源函数的原始形式和物理背景显著地不同。利用有效波能量平衡方程，将3／2指数律和发现的风能输入及破波耗散源函数相似性相结合，提出了深水风浪随风区成长的分式指数律，以得到的分式指数律拟合已有基于观测提出的风浪成长关系提出了半经验的风浪成长关系，与已有观测数据符合。     

大洋细结构的一种可能的机制Ⅲ.对混合过程动能耗散的估计

继第部分之后研究了惯性内波和近惯性内波由f～的作用所致的剪切不稳定引起的破碎机制。物理上,该机制很象存在由风应力所致薄表面涡旋漂流层时表面波的破碎与饱和过程。惯性内波和近惯性内波的破碎产物与小尺度湍流一起形成了混合块,它与Gregg等人(1986)的持久混合观测结果一致。依据Thorpe(1973)实验的结果作者提出了一个估计湍流动能耗散率和消衰时间的方法。结果表明,在剪切不稳定中近惯性内波在湍动耗散中起了关键作用,而惯性内波引起非常弱的湍动耗散。使用内波能量谱的标准总能量密度估计出的近惯性内波的耗散率和消衰时间与PATCHEX测量结果非常一致。文中还讨论了几个与此破碎机制有关的问题。     

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42, 134–149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313–351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313–351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971–5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift.     

An approximate solution of wave-modified Ekman current for gradually varying eddy viscosity

An approximate steady solution of the wave-modified Ekman current is presented for gradually varying eddy viscosity by using the WKB method with the variation of parameters technique. The parameters involved in the solution can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water. The solution reduces to the exact solution when the eddy viscosity is taken as a constant. As illustrative examples, for a fully developed wind-generated sea with different wind speeds and a few proposed gradually varying eddy viscosities, the current profiles calculated from the approximate solutions are compared with those of the exact solutions or numerical ones by using the Donelan and Pierson wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. It is shown that the approximate solution presented has an elegant form and yet would be valid for any given gradually varying eddy viscosity. The applicability of the solution method to the real ocean is discussed following the comparisons with published observational data and with the results from a large eddy simulation of the Ekman layer.     

Improving the estimate of wind energy input into the Ekman layer within the Antarctic Circumpolar Current

Based on the data and method offered by Liu et al. (2009), the direct wind and Stokes drift-induced energy inputs into the Ekman layer within the Antarctic Circumpolar Current (ACC) area are reestimated since the results of the former have been proved to be underestimated. And the result shows that the total rate of energy input into the Ekman-Stokes layer within the ACC area is 852.41 GW, including 649.75 GW of direct wind energy input (76%) and 202.66 GW of Stoke drift-induced energy input (24%). Total increased energy input, due to wave-induced Coriolis-Stokes forcing added to the classical Ekman model, is 52.05 GW, accounting for 6.5% of the wind energy input into the classical Ekman layer. The long-term variability of direct wind and Stokes drift-induced energy inputs into the Ekman layer within the ACC is also investigated, and the result shows that the Stokes drift hinders the decadal increasing trend of direct wind energy input. Meanwhile, there is a period of 4-5 a in the energy spectrums, as same as the Antarctic circumpolar wave.     

Global estimates of wind energy input to subinertial motions in the Ekman-Stokes layer

By incorporating the wave-induced Coriolis-Stokes forcing into the classical Ekman model, the wind energy input to the Ekman-Stokes layer is investigated, with an emphasis on the surface wave effects when the direction of Stokes drift deviates from that of wind stress. Theoretical analysis of the kinetic energy balance of the Ekman-Stokes layer shows that the total wind energy input consists of the direct wind energy input and the wave-induced energy input. Details of the direct wind and wave-induced energy input are discussed. Based on the ECMWF ERA-40 Re-Analysis wind stress and surface wave data, the global total wind energy input to subinertial motions in the Ekman-Stokes layer is estimated at 2.19 TW, including 0.26 TW (12%) wave-induced energy input and 1.93 TW (88%) direct wind energy input. The effect of sea-ice coverage on the energy input to the Ekman-Stokes layer is also considered. It is shown that the global total energy input could be overestimated by 0.08 TW (about 4%) without taking the sea-ice coverage into account.     

A note on surface wind-driven flow

A simple operationally oriented model of surface wind-driven currents is presented in which Lagrangian surface drift is assumed to be composed of a linear combination of a wave-induced Stokes drift plus a wind-driven Ekman drift. Using this approach, Stokes drift accounts for as much as half the total surface current magnitude. The Lagrangian current is predicted to be about 3.5% of the 10 m wind magnitude directed in the sense of an Ekman spiral about a 20° deviation angle. For comparison to this model, a second model is proposed that accounts for the interaction of Stokes current and Coriolis force. An inference drawn from this model is that there is only weak coupling between Coriolis force and Stokes drift. Such a conclusion, if correct, leads one to focus attention on the Lagrangian model for operationally oriented current estimates. Results of the Lagrangian model agree with observations of investigators for currents at the air-sea interface and may have application in the movement f oil slicks or surface drifters at sea under fetch or duration limited sea states.     

Vertical momentum transfer by internal waves when eddy viscosity and diffusion are taken into account

Free internal waves are considered in a Boussinesq approximation in the situation when horizontal eddy viscosity and diffusion in a vertically inhomogeneous flow are taken into account. The dispersion relation and wave damping factor are found in a linear approximation. The Stokes drift velocity is determined in the second order of smallness based on the wave amplitude. It has been indicated that the Stokes drift velocity, transverse with respect to the wave propagation direction, differs from zero if the flow-rate transverse component depends on the vertical coordinate. Vertical momentum fluxes differ from zero and can be comparable with or exceed the corresponding turbulent fluxes if eddy viscosity and diffusion are taken into account.     

波浪辐射应力对海流模式影响的数值分析

将波浪辐射应力,特别是地转意义下的波浪辐射应力引入海流数值模式POM(princeton ocean model),在渤海海域进行了初步的数值研究。在目前的数值分析中仅考虑了波浪辐射应力的横向分量(也是最重要的分量)。在POM模式中引入非地转和地转意义下的波浪辐射应力两种方案,并与原模式直接运行(即不考虑波浪辐射应力)的结果进行比较。比较显示,波浪辐射应力,特别是地转意义下的波浪辐射应力对海流模式结果的影响不容忽略。在海浪场存在的条件下,由风应力和地转意义下浪致作用力共同作用产生的海流强度应比理论上Ekman漂流的强度大,尤其是在浪致作用力显著的表层,表层流将明显增强,且不会完全符合Ekman漂流理论的转向规律。     

A note on ocean surface drift with application to surface velocities measured with HF Radar

The ocean drift current consists of a (local) pure drift current generated by the interaction of wind and waves at the sea surface, to which the surface geostrophic current is added vectorially. We present (a) a similarity solution for the wave boundary layer (which has been validated through the prediction of the 10-m drag law), from which the component of pure drift current along the direction of the wind (and hence the speed factor) can be evaluated from the 10-m wind speed and the peak wave period, and (b) a similarity solution for the Ekman layers of the two fluids, which shows that under steady-state neutral conditions the pure drift current lies along the direction of the geostrophic wind, and has a magnitude 0.034 that of the geostrophic wind speed. The co-existence of these two similarity solutions indicates that the frictional properties of the coupled air-sea system are easily evaluated functions of the 10-m wind speed and the peak wave period, and also leads to a simple expression for the angle of deflection of the pure drift current to the 10 m wind. The analysis provides a dynamical model for global ocean drift on monthly and annual time scales for which the steady-state neutral model is a good approximation. In particular, the theoretical results appear to be able to successfully predict the mean surface drift measured by HF Radar, which at present is the best technique for studying the near surface velocity profile.     

Energy budget of surface waves in the global ocean

Mechanical energy input from atmosphere and losses from wave-breaking dissipation of sea surface waves are estimated by a direct scheme. This scheme is based on the integration in the wavenumber space of the wind input and breaking dissipation source functions of the MASNUM wave model. The global amount of wind energy input, averaged in 2005, is about 57 TW, and the wave-breaking dissipation summed in deep-water is about 33 TW, over a half of the wind energy input. The residual may be dissipated by beach processes. Global distributions of the energy input and breaking dissipation concentrate in the westerlies of the Southern Hemisphere.     

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.