Evaluation of the simulation capability of the Wavewatch Ⅲ model for Pacific Ocean wave

Wave climate analysis and other applications for the Pacific Ocean require a reliable wave hindcast. Five source and sink term packages in the Wavewatch III model(v3.14 and v4.18) are compared and assessed in this study through comprehensive observations, including altimeter significant wave height, advanced synthetic aperture radar swell, and buoy wave parameters and spectrum. In addition to the evaluation of typically used integral parameters, the spectra partitioning method contributes to the detailed wave system and wave maturity validation. The modified performance evaluation method(PS) effectively reduces attribute numbers and facilitates the overall assessment. To avoid possible misleading results in the root mean square error-based validations, another indicator called HH(indicating the two authors) is also calculated to guarantee the consistency of the results. The widely used Tolman and Chalikov(TC) package is still generally efficient in determining the integral properties of wave spectra but is physically deficient in explaining the dissipation processes. The ST4 package performs well in overall wave parameters and significantly improves the accuracy of wave systems in the open ocean. Meanwhile, the newly published ST6 package is slightly better in determining swell energy variations. The two packages(ACC350 and BJA) obtained from Wavewatch III v3.14 exhibit large scatters at different sea states. The three most ideal packages are further examined in terms of reproducing waveinduced momentum flux from the perspective of transport. Stokes transport analysis indicates that ST4 is the closest to the NDBC-buoy-spectrum-based transport values, and TC and ST6 tend to overestimate and underestimate the transport magnitude, respectively, in swell mixed areas. This difference must be considered,particularly in air–wave–current coupling research and upper ocean analysis. The assessment results provide guidance for the selection of ST4 for use in a background Pacific Ocean hindcast for high wave climate research and China Sea swell type analysis.     

海浪波面极大值分布的非线性影响

依非线性海浪模型,在二阶近似下,利用特征函数展开技术和直接求矩的方法,导出了定点波面位移及其对时间的一阶和二阶导数的联合分布。由此联合分布,导出了二阶非线性近似下的波面极大值分布。它由线性意义下的海浪频谱及所考虑的二阶级。波相互作用所确定。当忽略非线性影响时,文中给出的波面极大值分布退化为Cartwright和Longuet-Higgins所导出的分布。     

吕宋海峡及周边海域湍流混合的时空分布特征研究

利用1992—2002年的温盐深数据与2012—2016年的Argo数据,基于细尺度参数化方法研究了吕宋海峡及周边海域(12°—30°N,115°—129°E)湍流混合的时空分布特征,并分析了地形粗糙度、内潮以及风输入的近惯性能通量对湍流混合的影响。结果表明,吕宋海峡和东海陆坡处具有强混合的特征,扩散率高达4×10~(-3) m~2/s,主要是由内潮产生导致的,其中吕宋海峡主要是M2、K1和O1内潮的贡献,而东海陆坡处主要是M_2内潮的贡献;南海北部也呈现较强的混合,且陆坡处的混合比海盆高1—2个量级;南海中央海盆和离岸的菲律宾海混合较弱,扩散率为O (10-5 m2/s)。此外,在研究区域内,湍流混合的年际变化和季节变化均不明显,且混合扩散率与风输入的近惯性能通量未表现出明显的季节相关。     

利用TOPEX卫星高度计资料分析东中国海的风、浪场特征

利用TOPEX卫星高度计和日本气象厅浮标观测资料,对东中国海的有效波高和风速进行比较,分析了卫星高度计资料的有效性。利用有效波高和风速的3种概率密度函数分布,结合TOPEX卫星高度计资料,并采用最大似然方法对统计分布参数进行估计,结果表明,有效波高的对数-正态概率密度分布与观测资料的直方图在有效波高的整个范围内符合较好,风速的直方图与Weibul概率密度分布符合较好。同时,分析了有效波高大于4 m的巨浪在东中国海的时空分布特征,表明巨浪多出现在冬、秋两季,平均有效波高最大值出现在夏季,且主要分布在东中国海东南部。     

三层流体系统非线性界面内波传播理论的研究

基于波陡很小的假设, 利用摄动法, 讨论了任意深度的二维不可压缩、无黏性、无旋的三层流体系统. 在刚性上边界、平底不可渗透条件下, 给出了界面内波传播的统一理论以及描述其波剖面的近似非线性演化方程(NEEs). 最后讨论了几种特殊情形下的近似NEEs.结果表明文献导出的理论结果为本文的特殊情形.     

Strongly nonlinear long interfacial waves in uniform basic flows and their dispersion relationship

In this paper, long interfacial waves of finite amplitude in uniform basic flows are considered with the assumption that the aspect ratio between wavelength and water depth is small. A new model is derived using the velocities at arbitrary distances from the still water level as the velocity variables instead of the commonly used depth-averaged velocities. This significantly improves the dispersion properties and makes them applicable to a wider range of water depths. Since its derivation requires no assumption on wave amplitude, the model thus can be used to describe waves with arbitrary amplitude.     

Modeling vertical heat transport in the wave affected surface layer of the ocean

In considering the vertical heat transport problems in the upper ocean, the flat upper boundary approximation for the free surface and the horizontal homogenous hypothesis are usually applied. However, due to the existence of the wave motion, the application of this approximation may result in some errors to the solar irradiation since it decays quickly in respect to the actual thickness of the water layer below the surface; on the other hand, due to the fluctuation of the water layer depth, it is improper to neglect the effects of the horizontal advection and turbulent diffusion since they also contribute to the vertical heat transport. A new model is constructed in this study to reflect these effects. The corresponding numerical simulations show that the wave motion may remarkably accelerate the vertical heat transferring process and the variation of the temperature in the wave affected layer appears in an oscillating manner. Supported by the National High Technology Research and Development Program of China (863 Program, No. 2006AA09A309); China Postdoctoral Science Foundation (No. 20070411111) and the Fund of Shandong Province for the Excellent Post-Doctors (No. 200603056)     

塔基海-气界面动量通量观测：向岸风和离岸风情形比较

本文使用塔基直接观测法研究海洋大气边界层中的海-气界面动量通量。首先,我们收集数据并和前人观测结果比对,其比对结果符合一致。其次,在低风速至中等风速条件下,我们发现海-气界面动量通量的交换系数（又称拖曳系数）对于向岸风和离岸风两种情形有所差异。为此,我们使用一个考虑表面波的参数化方案解释海-气界面动量通量对于表面波的依赖关系。这些结果一方面证实表面波对于海-气界面动量通量的影响,另一方面验证一个考虑表面波参数化方案的有效性。     

低风速涌浪条件下浪致应力的研究

本文基于2012年2月15日至5月8日的广东省茂名市附近海上铁塔实测数据,分析了海气边界层的湍流运动。结果证实,在涌浪条件下,协谱和Ogive曲线存在明显受涌浪影响的区域。本文进一步研究,发现浪致雷诺应力的正负符号变化与波龄相关,并给出了顺风向和侧风向浪致雷诺应力正负符号反转时波龄的统计分布。此外,本文基于协谱法提取了浪致雷诺应力τ_swell,并与传统方法得到的浪致雷诺应力τ_*swell做比较。结果表明,整体上τ_swell比τ_*swell大2～3个量级,这意味着传统方法显著低估了涌浪的浪致雷诺应力,传统的浪致雷诺应力计算方法还有待改进。     

Large eddy simulation of the rotation effect on the ocean turbulence kinetic energy budget in the surface mixed layer

A non-hydrostatic, Boussinesq, and three-dimensional large eddy simulation(LES) model was used to study the impact of the Earth's rotation on turbulence and the redistribution of energy in turbulence kinetic energy(TKE) budget. A set of numerical simulations was conducted,(1) with and without rotation,(2) at different latitudes(10°N, 30°N, 45°N, 60°N, and 80°N),(3) with wave breaking and with Langmuir circulation, and(4) under different wind speeds(5, 10, 20, and 30 m/s). The results show that eddy viscosity decreases when rotation is included, indicating that rotation weakens the turbulence strength. The TKE budget become tight with rotation and the effects of rotation grow with latitude. However, rotation become less important under Langmuir circulation since the transport term is strong in the vertical direction. Finally, simulations were conducted based on field data from the Boundary Layer and Air-Sea Transfer Low Wind(CBLAST-Low) experiment. The results, although more complex, are consistent with the results obtained from earlier simulations using ideal numerical conditions.