SWAN模型中不同风拖曳力系数对风浪模拟的影响

本文以荷兰哈灵水道海域为实验区域,通过敏感性实验,研究了在14 m/s、31.5 m/s和50 m/s（分别代表一般大风、强热带风暴和强台风的极端条件）定常风速下SWAN模型中不同风拖曳力系数对风浪模拟的影响程度。结果表明,对于近岸浅水区域（水深小于20 m）,风拖曳力系数计算方案的选择对有效波高影响较小,而且当风速增加到一定程度后,波浪破碎成为影响波高值的主要因素;对于深水区域（水深大于30 m）,一般大风条件下风拖曳力系数计算方案的选择对有效波高影响仍然较小,随着风速的继续增大,风拖曳力系数计算方案的选择对有效波高的影响逐渐显著。对于平均周期,风拖曳力系数计算方案的选择和风速的改变对其影响均较小,而由水深变浅导致的波浪破碎对其影响较为显著。根据敏感性实验结果,本文对SWAN模型中风拖曳力系数计算方案的选择做出如下建议:计算近岸浅水区域风浪场或深水区域一般大风条件风浪场时,其风拖曳力系数可以直接采用模型默认选项;而对于深水区域更大风速条件,可首先采用模型默认选项试算,然后结合当地海域实测波浪资料进行修正。     

Dissipation observations of drag coefficients over the open ocean

Forty-three open-ocean observations of drag coefficients observed at Argus Island Tower near Bermuda by the dissipation technique resulted in constant drag coefficients for mean horizontal wind velocities between 7.8 and10.4 m .s^{1}in good agreement with the larger near-neutral data set of DeLeonibus and Simpson [1] and the neutral data of Large and Pond [2], both of whom observed10^{3}C_{10} = 1.2whereC_{10}is the drag coefficient at l0 m. Ratios of vertical-to-horizontal wind velocity spectral densities averaged over an inertial subrange of 0.8 to 1.6 Hz ranged from 0.7 to 1.07 in agreement with the Busch and Panofsky [3] result that isotropy is approached only when the observation height is much greater than the Nyquist wavelength.     

Wind speed scaling and the drag coefficient

[1]Banner M L, Chen W, Walsh E J, et al.     

Form Drag Caused by Topographically Forced Waves in a Barotropic β Channel: Effect of Higher Mode Resonance

The relationship between form drag and the zonal mean velocity of steady states is investigated in a very simple system; a barotropic quasi-geostrophic β channel with sinusoidal topography. When a steady solution is calculated by the modified Marquardt method, keeping the zonal mean velocity constant as a parameter, the characteristic of the solution changes at a phase speed of a wave with a wavenumber higher than that of the bottom topography. For velocities smaller than this critical value, there exists a stable quasi-linear solution similar to the linear solution. For larger velocities, there exist three solutions whose form drag is very large. In addition, the resonant velocity of the mode, whose wavenumber is the same as the bottom topography, has no effect on these solutions. When the quiescent fluid is accelerated by a constant wind stress, acceleration stops around the critical velocity for a wide range of the wind stress. If the wind stress is too large for the acceleration to stop, the zonal current speed continues to increase infinitely. It is implied that the zonal velocity of equilibrium is mainly determined, not by the wind stress, but by the amplitude of the bottom topography and the viscosity coefficient.     

两参量的海面阻力系数模式的探讨

从风浪的能量平衡方程出发，引进若干风要素与波要素以及波要素之间的定性关系，经演算可导出海面阻力系数（Ｃｐ）或是风速（Ｕ）和波龄（β）或是Ｕ和波高（Ｈ）的函数，然后沿用最小二乘法，终将得出４组１２个回归方程。当β（或β）或Ｈ为某一给定值，惟有Ｕ为唯一参量时，所提各式均可简化为非线性方程：ＣＤ＝ａ＋ｂ，Ｕ＋ｃ．Ｕ＾２；式中ａ，ｂ和ｃ为三个经验系数，就所检验的例子而言，本文的结果与实际的符合前人的为好。     

Effect of the drag coefficient on a typhoon wave model

The effect of the drag coefficient on a typhoon wave model is investigated.Drag coefficients for Pingtan Island are derived from the progress of nine typhoons using COARE 3.0 software.The wind parameters are obtained using the Weather Research and Forecasting model.The simulation of wind agrees well with observations.Typhoon wave fields are then simulated using the third-generation wave model SWAN.The wave model includes exponential and linear growths of the wind input,which determine the wave-growth mode.A triple triangular mesh is adopted with spatial resolution as fine as 100 m nearshore.The SWAN model performs better when using the new drag coefficient rather than the original coefficient.     

海-气界面动量通量的估计方法分析与应用

首次将经验模态分解方法引入湍流稳定性分析,与传统的线性和滑动平均去势方法进行了比较,发现经验模态分解方法的去势效果最好。基于"南海平台通量观测计划"(FOPSCS)近两年的连续通量观测数据,得到了22 476个摩擦速度的估算值,结果表明,当风速小于5m/s时,拖曳系数随风速增大而减小,而风速大于5m/s时,拖曳系数随风速增大而增大,两种情形分别反映了黏性表皮摩擦和波浪引起的形状阻力对海面风应力的贡献。同时发现短风区情形的拖曳系数大于长风区情形,说明波浪成长状态会对海-气界面动量交换产生影响。     

Quasi-linear model of interaction of surface waves with strong and hurricane winds

A quasi-linear model for determining the aerodynamic drag coefficient of the sea surface and the growth rate of surface waves under a hurricane wind is proposed. The model explains the reduction (stabilization) in the drag coefficient during hurricane winds. This model is based on the solution of the Reynolds equations in curvilinear coordinates with the use of the approximation of the eddy viscosity, which takes into account the presence of the viscous sublayer. The profile of the mean wind velocity is found with consideration for nonlinear wave stresses (wave momentum flux), whereas wave disturbances induced in air by waves on the water surface are determined in the context of linear equations. The model is verified by comparing the calculation results with experimental data for a wide range of wind velocities. The growth rate and drag coefficient for hurricane winds are calculated both with and without consideration for the shortwave portion of the windwave spectrum. On the basis of calculations with the quasi-linear model, a simple parametrization is proposed for the drag coefficient and the growth rate of surface waves during hurricane winds. This model is convenient for use in models of forecasting winds and waves.     

引入拖曳系数参数化的海冰自由漂流模拟研究

Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient(C_a) and oceanic drag coefficient(C_w).Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions.In the present study,the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model,and the wind factor α and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of C_a to C_w are studied to investigate their dependence on the impact factors such as local drag coefficients,floe and ridge geometry.The results reveal that in an idealized steady ocean,C_a/C_w increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone,while it remains at a steady level(0.2-0.25) for large floes in the central ice zone.The wind factor α increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%.And the deflection angle θ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than20%without a steady level like α.The values of these parameters agree well with the previously reported observations in Arctic.The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient(C_s' /C_s) and the ratio of ridge form drag coefficient(C_r'/C_r) to the value of C_a/C_w,α,and θ,because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces.Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.     

风应力拖曳系数选取对风暴潮数值模拟的影响

在风暴潮的形成中风应力起决定性作用 ,风应力拖曳系数决定了大气与海洋间的动量传输率。观测结果表明 ,风应力拖曳系数随风速而变化 ,与海面粗糙度有关。文中采用几种与风速有关的风应力拖曳系数表达式进行数值模拟 ,与将其视为常数情况相比较 ,计算结果的精度均有较明显提高。对比各表达式模拟结果 ,采用 Smith(1980 )风应力拖曳系数公式的模拟效果为最好     

波流共同作用下珊瑚礁冠层附近平均流分布及阻力特性试验

本文通过物理模型试验研究了波流共同作用下珊瑚礁冠层附近平均流的分布特征以及阻力特性,分析了典型波浪工况下无潮流、正向潮流和反向潮流分别作用下平均流速、摩阻流速、阻力系数的沿礁变化规律。结果表明：无潮流时礁前斜坡及外礁坪上存在海底回流且在礁缘附近回流最强,在礁坪上冠层附近平均流表现为向岸流,且该流沿礁向海岸方向持续增大。相较于无潮流时,正向潮流作用下冠层内外均为向岸流,在礁坪上冠层内外的向岸流显著增大;反向潮流作用下冠层内外均为离岸流且在礁缘处达到最大,该离岸流在礁坪上逐渐减小然后趋于稳定。无潮流时礁坪上摩阻流速呈小幅波动;相较于无潮流时,正向潮流、反向潮流影响下礁坪上摩阻流速显著增大,其中正向潮流影响下增幅更大;无潮流时礁坪上水力粗糙度沿礁减小,正向潮流和反向潮流影响下水力粗糙度普遍有幅度不等的增加。三种工况下礁坪上的阻力系数均沿礁整体呈下降趋势,相较于无潮流时,正向潮流和反向潮流影响下礁坪上的阻力系数显著增大,且正向潮流作用时增幅更大。     

Laboratory study on wave dissipation by vegetation in combined current–wave flow

Coastal wetlands such as salt marshes and mangroves provide valuable ecosystem services including coastal protection. Many studies have assessed the influence of plant traits and wave conditions on vegetation-induced wave dissipation, whereas the effect of tidal currents is often ignored. To our knowledge, only two studies investigated wave dissipation by vegetation with the presence of following currents (current velocity is in the same direction as wave propagation) (Li and Yan, 2007; Paul et al., 2012). However, based on independent experiments, they have drawn contradictive conclusions whether steady currents increase or decrease wave attenuation. We show in this paper that this inconsistency may be caused by a difference in ratio of imposed current velocity to amplitude of the horizontal wave orbital velocity. We found that following currents can either increase or decrease wave dissipation depending on the velocity ratio, which explains the seeming inconsistency in the two previous studies. Wave dissipation in plant canopies is closely related to vegetation drag coefficients. We apply a new approach to obtain the drag coefficients. This new method eliminates the potential errors that are often introduced by the commonly used method. More importantly, it is capable of obtaining the vegetation drag coefficient in combined current–wave flows, which is not possible for the commonly used calibration method. Based on laboratory data, we propose an empirical relation between drag coefficient and Reynolds number, which can be useful for numerical modeling. The characteristics of drag coefficient variation and in-canopy velocity dynamics are incorporated into an analytical model to help understand the effect of following currents on vegetation-induced wave dissipation.     

风浪和海洋飞沫对海表面拖曳系数和风廓线的影响

基于埃克曼理论,本文将波致应力和飞沫应力引入到海-气边界层的界面应力中,来研究海表面风浪和海洋飞沫对海-气边界层动量交换的影响,并得到修改后的埃克曼模型的理论解。波致应力是由风浪谱和波增长函数估计,并得到在中低风速下,波致应力、飞沫应力与湍流应力相比,对海表面拖曳系数和风廓线的影响非常小。当风速高于25米/秒时,海洋飞沫通过飞沫应力对海-气界面应力的作用远高于波致应力,以至于波致应力可以忽略。海表面拖曳系数在高风速下,随着风速的增大而减小。通过采用风浪谱的不同波龄,得到海洋飞沫的产生会导致海-气边界层风速的增加。最后,理论解与现场的观察数据进行了对比。对比后的数据表明,在中高风速下,飞沫对海-气边界层的影响远大于表面风浪。     

Statistical study on the local equilibrium between wind and wind waves by using data from ocean data buoy stations

The local equilibrium between the wind and wind waves, which is defined by a range of the coefficient of the 3/2-power law between the non-dimensional significant wave height and period, is statistically investigated by using wind and wave data obtained at four ocean data buoy stations in the seas near Japan. The friction velocity is calculated from the wind speed measured at one height together with the significant wave period by using formulas of the wave dependent drag coefficient proposed by Tobaet al. (1990). The data for small waves or for weak winds indicate that the waves do not satisfy the criterion for the local equilibrium, because they may be affected by changing winds or remotely generated swells. In the seas near Japan, the data which satisfy the local equilibrium are about 6% through a year. Otherwise swells are dominant in most situations. Changing winds also cause deviations from the local equilibrium. The degree of satisfaction of the local equilibrium can be classified by ranges of the significant wave height. As the significant wave height exceeds 4 m, the local equilibrium is more frequently satisfied.     

OSMAR-S100便携式高频地波雷达海浪和海面风探测性能分析

OSMAR-S系列便携式高频地波雷达系统采用单极子/交叉环紧凑型天线阵,通过单站雷达即可实现有效探测距离约10km内海浪和海面风的单点观测。为了更好地了解OSMAR-S100雷达系统海浪和海面风的综合探测性能,于2013年1月29日至3月7日在台湾海峡西南部海域进行了雷达与浮标观测的对比试验,得到了有效波高、有效波周期、平均风速和平均风向数据。对比结果表明,OSMAR-S100便携式高频地波雷达可有效观测距雷达10km以内有效波高0.5m以上的海浪平均状况和平均风速5m/s以上的海面风,雷达反演有效波高和有效波周期的均方根误差分别为0.60m和1.60s,反演平均风速和平均风向的均方根误差为1.83m/s和16.7°。在未经区域化标定的情况下,此结果说明了该型雷达产品已初步具备了海浪和海面风的业务化观测水平。     

A Practical bi-parameter formula of gas transfer velocity depending on wave states

The parameter that describes the kinetics of the air-sea exchange of a poorly soluble gas is the gas transfer velocity which is often parameterized as a function of wind speed. Both theoretical and experimental studies suggest that wind waves and their breaking can significantly enhance the gas exchange at the air-sea interface. A relationship between gas transfer velocity and a turbulent Reynolds number related to wind waves and their breaking is proposed based on field observations and drag coefficient formulation. The proposed relationship can be further simplified as a function of the product of wind speed and significant wave height. It is shown that this bi-parameter formula agrees quantitatively with the wind speed based parameterizations under certain wave age conditions. The new gas transfer velocity attains its maximum under fully developed wave fields, in which it is roughly dependent on the square of wind speed. This study provides a practical approach to quantitatively determine the effect of waves on the estimation of air-sea gas fluxes with routine observational data.     

Internal flow structure of short wind waves

For wind waves generated in a wind-wave tunnel, the surface pressure and also the pressure distribution along the internal streamlines were calculated from the measured internal velocity field. In distinct waves, with wave height comparable with or larger than the mean, the surface pressure is found to vary drastically in a narrow region around the crest, showing a dominant minimum near the crest. On the other hand, the pressure distribution along the streamline shows systematic variations that are nearly in phase with the streamline profile. It is shown that the occurrence of the pressure in phase with the streamline profile is linked with the internal vorticity distribution, especially with the presence of a high vorticity region below the crest described in Part I of this study. As a result of the occurrence of such pressure variations, the dispersion relation is modified by about 10% from that for linear irrotational waves. It is argued from the present measurements that the dispersion relation and also the energy transfer from wind into wind waves are strongly affected by the internal vortical structure so that the assumption of irrotational gravity waves cannot be applied to the wind waves being studied.     

强风条件下海面拖曳系数和粗糙长度研究

进一步研究强风条件下海-气湍流动量交换以及海浪特征,有助于提高数值天气模式对台风强度演变、移动路径以及恶劣海况的预报能力。依照前人的方法将台风分为风向与浪向(1)相同,(2)相反,和(3)交叉3个扇形区,并结合台风路径数据,得到了浮标数据相对于台风的方位。分别对3种类型的浮标数据进行分析,进而发现了波浪高度和相速度随风速增加而变化的规律。并利用GWW参数化方案计算出摩擦速度(u_*)、拖曳系数(C_DN)和粗糙长度(z₀)。将这些结果与前人代表性的研究论文中所用观测数据和所得研究结论进行比较,结果表明二者有较强的一致性。该研究证明GWW参数化方案在强风条件下依然有很好的适用性。     

Internal flow structure of short wind waves

The minimum value of wind stress under which the flow velocity in short wind waves exceeds the phase speed is estimated by calculating the laminar boundary layer flow induced by the surface tangential stress with a dominant peak at the wave crest as observed in previous experiments. The minimum value of the wind stress is found to depend strongly on, the ratio of the flow velocity just below the boundary layer and the phase speed, but weakly onL, the wavelength. For wind waves previously studied (=0.5,L=10 cm), the excess flow appears when the air friction velocityu _* is larger than about 30 cm sec^–1. The present results confirm that the excess flow found in my previous experiments is associated with the local growth of a laminar boundary layer flow near the wave crest.     

风浪对海-气界面动量通量估计的影响

利用实验室和有代表性的外海观测数据综合分析表明,海面粗糙度对波龄的依赖性与是否将实验室和外海数据一起考虑有关,而风应力拖曳系数与此无关,且随波龄增大而减小.利用Toba-3/2指数律和风浪成长关系的分析表明,风应力拖曳系数为常数或随波龄的增大而增大,与上述结果定性上相矛盾,说明风浪对风应力拖曳系数影响问题需要进一步研究.