Impact Of Dominant Waves On Sea Drag

The impact of air-flow separation from breaking dominant waves is analyzed.This impact results from the correlation of the pressure drop with theforward slope of breaking waves. The pressure drop is parameterized via thesquare of the reference mean velocity. The slope of breaking waves isrelated to the statistical properties of the wave breaking fronts describedin terms of the average total length of breaking fronts. Assuming that thedominant waves are narrow and that the length of breaking fronts is relatedto the length of the contour of the breaking zone it is shown that theseparation stress supported by dominant waves is proportional to thebreaking probability of dominant waves. The breaking probability of dominantwaves, in turn, is defined by the dominant wave steepness. With thedominant wave steepness increasing, the breaking probability is increasedand so does the separation stress. This mechanism explains wave age (youngerwaves being steeper) and finite depth (the spectrum is steeper in shallowwater) dependence of the sea drag. It is shown that dominant waves support asignificant fraction of total stress (sea drag) for young seas due to theair-flow separation that occurs when they break. A good comparison of themodel results for the sea drag with several data sets is reported.     

荒漠戈壁大气总体曳力系数和输送系数观测研究

利用“我国西北干旱区陆一气相互作用观测试验”在甘肃敦煌进行的陆面过程野外试验的观测资料,依据三种不同方法确定了干旱戈壁区动量输送的曳力系数Cd、感热和潜热交换的总体输送系数Ch和Cq。结果表明：尽管这三种方法计算的曳力系数和总体输送系数有一定的差别,但在量级上相当,尤其是Cd和Ch的平均值比较接近。本文还通过对风向的分析,剔除了附近建筑物干扰和来自绿洲湿平流的影响,得到了荒漠戈壁总体输送系数的特征及其与理查孙数的关系。     

Effects of Drag Coefficients on Surface Heat Flux during Typhoon Kalmaegi (2014)

The lack of in situ observations and the uncertainties of the drag coefficient at high wind speeds result in limited understanding of heat flux through the air-sea interface and thus inaccurate estimation of typhoon intensity in numerical models. In this study, buoy observations and numerical simulations from an air-sea coupled model are used to assess the surface heat flux changes and impacts of the drag coefficient parameterization schemes on its simulations during the passage of Typhoon Kalmaegi (2014). Three drag coefficient schemes, which make the drag coefficient increase, level off, and decrease, respectively, are considered. The air-sea coupled model captured both trajectory and intensity changes better than the atmosphere-only model, though with relatively weaker sea surface cooling (SSC) compared to that captured by buoy observations, which led to relatively higher heat flux and thus a stronger typhoon. Different from previous studies, for a moderate typhoon, the coupled simulation with the increasing drag coefficient scheme outputted an intensity most consistent with the observation because of the strongest SSC, reasonable ratio of latent and sensible heat exchange coefficients, and an obvious reduction in the overestimated surface heat flux among all experiments. Results from sensitivity experiments showed that surface heat flux was significantly determined by the drag coefficient-induced SSC rather than the resulting wind speed changes. Only when SSC differs indistinctively (<0.4°C) between the coupled simulations, heat flux showed a weak positive correlation with the drag coefficient-impacted 10-m wind speed. The drag coefficient also played an important role in decreasing heat flux even a long time after the passage of Kalmaegi because of the continuous upwelling from deeper ocean layers driven by the impacted momentum flux through the air-sea interface.     

Drag Partition for Regularly-Arrayed Rough Surfaces

Vegetation and other roughness elements distributed across a surface can providesignificant protection against wind erosion by extracting momentum from the flowand thereby reducing the shear stress acting at the surface. A theoretical model haspreviously been presented to specify the partition of drag forces for rough surfacesand to predict required vegetation density to suppress wind erosion. However, themodel parameters have not yet been constrained and the predictive capacity of themodel has remained uncertain. A wind-tunnel study was conducted to measure thedrag partition for a range of roughness densities and to parameterise the model inorder to improve its range of potential applicability. The drag forces acting on bothan array of roughness elements and the intervening surface were measured independentlyand simultaneously using new drag balance instrumentation. A detailed measure of thespatial heterogeneity of surface shear stresses was also made using Irwin sensors. Thedata agreed well with previous results and confirmed the general form of the model.Analysis of the drag partition confirmed the parameter definition = C_R/C_S(where C_R and C_S are roughness element and surface drag coefficients,respectively) and a constant proportional difference between the mean and maximumsurface shear stress was found. The results of this experiment suggest that the definitionfor m, the surface shear stress inhomogeneity parameter, should be revised, although thetheoretical and physical reasons for including this parameter in the model appear to bevalid. Best-fit values for m ranged from 0.53 to 0.58.     

Determination Of The Surface Drag Coefficient

This study examines the dependence of the surface drag coefficienton stability, wind speed, mesoscale modulation of the turbulent flux and method of calculation of the drag coefficient. Data sets over grassland, sparse grass, heather and two forest sites are analyzed. For significantly unstable conditions, the drag coefficient does not depend systematically on z/L but decreases with wind speed for fixed intervals of z/L, where L is the Obukhov length. Even though the drag coefficient for weak wind conditions is sensitive to the exact method of calculation and choice of averaging time, the decrease of the drag coefficient with wind speed occurs for all of the calculation methods. A classification of flux calculation methods is constructed, which unifies the most common previous approaches.The roughness length corresponding to the usual Monin–Obukhovstability functions decreases with increasing wind speed. This dependence on wind speed cannot be eliminated by adjusting the stability functions. If physical, the decrease of the roughness length with increasing wind speed might be due to the decreasing role of viscous effectsand streamlining of the vegetation, although these effects cannot be isolated from existing atmospheric data.For weak winds, both the mean flow and the stress vector often meander significantly in response to mesoscale motions. The relationship between meandering of the stress and wind vectors is examined. For weak winds, the drag coefficient can be sensitive to the method of calculation, partly due to meandering of the stress vector.     

The Impact Of Air-Flow Separation On The Drag Of The Sea Surface

An approach that allows assessment ofthe impact of air-flow separation (AFS) fromwave breaking fronts on the sea-surface drag is presented. Wave breaking fronts are modelled by the discontinuities of the sea-surface slope. It is assumedthat the dynamics of the AFS from wave breaking crests is similar to thatfrom the backward facing step. The form drag supported by an individualbreaker is described by the action of the pressure drop distributed alongthe forward face of the breaking front. The total stress due to the AFS isobtained as a sum of contributions from breaking fronts of different scales.Outside the breaking fronts the drag of the sea surface is supported by theviscous surface stress and the wave-induced stress. To calculate the stressdue to the AFS and the wave-induced stress a physical model of the wind-wavespectrum is used. Together with the model of the air flow described in termsof surface stresses it forms a self-consistent dynamical system for the seasurface-atmosphere where the air flow and wind waves are strongly coupled.Model calculations of the drag coefficient agree with measurements. It is shownthat the dimensionless Charnock parameter (roughness length normalized onthe square of the friction velocity and the acceleration of gravity)increases with the increase of the wind speed in agreement with fieldmeasurements. The stress due to the AFS normalized on the square of thefriction velocity is proportional to the cube of wind speed. At low windsthe viscous surface stress dominates the drag. The role of the form drag,which is the sum of the stress due to the AFS and the wave-induced stress, isnegligible. At moderate and high winds the form drag dominates. At windspeeds higher than 10 m s^-1 the stress supported by the AFS becomescomparable to the wave-induced stress and supports up to 50% of the totalstress.     

A Note on the Drag of the Sea Surface at Hurricane Winds

Based on the solution of the turbulent kinetic energy balance equation for the airflow in the regime of limited saturation by suspended sea-spray droplets, some experimental evidence, and simple arguments, a resistance law of the sea surface at hurricane winds is derived. It predicts the reduction of the drag coefficient for the wind speed exceeding hurricane values of 30–40 m s ^-1 in agreement with field data.     

Numerical Simulation of Drag Partition over Rough Surfaces

We present a numerical simulation of drag partition over rough surfaces. A computational fluid dynamics model is applied with high resolution to simulatingturbulent flows over arrays of roughness elements positioned on asmooth surface. The skin drag on the surface and the pressure drag on the roughnesselements are computed. The simulated drag partition compares well with wind-tunnelmeasurements and theoretical estimates for similar rough surfaces. This confirms that the computational approach offers an alternative to wind-tunnel and field experiments in studying drag and drag partition. The model is then applied to studying drag partition over rough surfaces with various roughness configurations. It is shown that drag partition depends not only on the magnitude of the roughness frontal area but also on the sizes and arrangement of roughness elements, because (1) the pressure drag coefficient is sensitive to roughness-element dimensions and (2) the arragement of roughness elements lead to different interferences of turbulent wakes. The impact ofthe latter factor is not insignificant.     

1979—2013年ERA-Interim资料的青藏高原低涡活动特征分析

为了更好地了解青藏高原多尺度地形的动力作用,并为改良数值模式中地形的表示方法奠定基础,通过采用2010年青藏高原西南部6个地面台站的观测资料以及4种不同分辨率的分析(再分析)资料,分别估算了冈底斯山及整个青藏高原主体范围内的地表气压拖曳,得出了青藏高原可能存在的拖曳类型,并且分析了青藏高原气压拖曳的一些特征。得出如下主要结论:由罗斯贝波产生的波动拖曳作为行星尺度的拖曳对青藏高原地区总拖曳的贡献最大;同时,青藏高原范围内存在着大量与天气过程密切相关的天气尺度的拖曳;对于冈底斯山对气流的中尺度动力作用的进一步分析可知,夏季基本全为气流分离,冬季500 hPa以下为气流分离,500—200 hPa为气流分离和波动破碎的混合区,而200 hPa以上的平流层则为重力波的产生及其破碎区域;冈底斯山地区的地表气压拖曳主要集中在3000—5000 m高度,并且,冈底斯山总拖曳的方向近乎与山脊垂直;地表气压和地形高度资料的分辨率越高,所能分辨出的更小波长的气压拖曳也越多,估算出的高原主体范围内的拖曳值也越大;变压梯度和地形梯度是影响气压拖曳的基本因子,但地形梯度对拖曳的影响最终是通过气压梯度来实现的。     

A Note on a Parameterization of the Sea Drag

A new parameterization of the sea drag is based on a wind-over-wavescoupling theory. The parameterization accounts for the wind speed, wave ageand finite depth dependencies of the sea drag. The latter two are introducedthrough the integral parameters of the wind-wave field: the dominant waveheight and the wavenumber at the spectral peak, and the water depth. Theparameterization is checked against the wind-over-waves model results andtwo field datasets obtained in a wide range of the wind speed and wave age.The comparison is encouraging. The parameterization is aimed for use inoperational ocean-state and atmosphere models.     

Laboratory Studies Of Wind Stress Over Surface Waves

Simultaneous laboratory observations of wind speed, wind stress, and surfacewind-wave spectra are made under a variety of wind forcing patterns using cleanwater as well as water containing an artificial surfactant. Under typical experimentalconditions, more than half of the total stress is supported by the wave-induced stressrather than by the surface viscous stress. When the surfactant reduces the shortwind-wave spectra, the wind stress also decreases by as much as 20–30% at agiven wind speed. When the wind forcing is modulated in time, the wind stresstends to be higher under decreasing wind than under increasing wind at a givenwind speed, mainly because the response of short wind-wave spectra to varyingwind forcing is delayed in time. These examples clearly demonstrate that therelationship between the wind speed and the wind stress can be significantlymodified if the surface wave field is not in equilibrium with the wind forcing.Next, we examine whether the wind stress is estimated accurately if the wave-inducedstress by all surface wave components is explicitly evaluated by linear superpositionand is added to the surface viscous stress. It is assumed that the surface viscous stressis uniquely related to the wind speed, and that the wind input rate is determined by thelocal, reduced turbulent stress rather than the total stress. Our wind stress estimatesincluding the wave contributions agree well with observed wind stress values, evenif the surface wave field is away from its equilibrium with the wind in the presenceof surface films and/or under time-transient wind forcing. These observations stronglysuggest that the wind stress is accurately evaluated as a sum of the wave-induced stressand the surface viscous stress. At very high winds, our stress estimates tend to be lowerthan the observations. We suspect that this is because of the enhancement of wind stressover very steep (or breaking) short wind-waves.     

‘Surface Drag in the Arctic Marginal Sea-ice Zone: A Comparison of Different Parameterisation Concepts’

Two parameterisation schemes for the turbulent surface fluxes and drag coefficients over the Arctic marginal sea-ice zone (MIZ) are (further) developed, and their results are compared with each other. Although the schemes are based on different principles (flux averaging and parameter averaging), the resulting drag coefficients differ only slightly in the case of neutral and stable stratification. For unstable stratification and sea-ice conditions being typical for the north-eastern Fram Strait, the drag coefficients resulting from the parameter-averaging concept are 5–10% larger than those of the flux-averaging concept. At a sea-ice concentration of 45%, the parameter-averaging method overestimates the heat fluxes by a factor of 1.2. An inclusion in the schemes of form drag caused by floe edges and ridges has a much larger effect on the drag coefficient, and on the momentum fluxes, than the choice between the parameter-averaging or flux-averaging methods. Based on sensitivity studies with the flux-averaging scheme, a simple formula for the effective drag coefficient above the Arctic MIZ is derived. It reduces the computational costs of the more complex parameterisations and could also be used in larger scale models. With this simple formula, the effective drag coefficient can be calculated as a function of the sea-ice concentration and skin drag coefficients for water and ice floes. The results obtained with this parameterisation differ only slightly from those using the more complex schemes. Finally, it is shown that in the MIZ, drag coefficients for sea-ice models may differ significantly from the effective drag coefficients used in atmospheric models.     

Nonlinear Formulation of the Bulk Surface Stress over Breaking Waves: Feedback Mechanisms from Air-flow Separation

Historically, our understanding of the air-sea surface stress has been derived from engineering studies of turbulent flows over flat solid surfaces, and more recently, over rigid complex geometries. Over the ocean however, the presence of a free, deformable, moving surface gives rise to a more complicated drag formulation. In fact, within the constant-stress turbulent atmospheric boundary layer over the ocean, the total air-sea stress not only includes the traditional turbulent and viscous components but also incorporates surface-wave effects such as wave growth or decay, air-flow separation, and surface separation in the form of sea-spray droplets. Because each individual stress component depends on and alters the sea state, a simple linear addition of all stress components is too simplistic. In this paper we present a model of the air-sea surface stress that incorporates air-flow separation and its effects on the other stress components, such as a reduction of the surface viscous stress in the separated region as suggested by recent measurements. Naturally, the inclusion of these effects leads to a non-linear stress formulation. This model, which uses a variable normalized dissipation rate of breaking waves and normalized length of the separation bubble, reproduces the observed features of the drag coefficient from low to high wind speeds despite extrapolating empirical wave spectra and breaking wave statistics beyond known limits. The model shows the saturation of the drag coefficient at high wind speeds for both field and laboratory fetches, suggesting that air-flow separation over ocean waves and its accompanying effects may play a significant role in the physics of the air-sea stress, at least at high wind speeds.     

Optimized Computation of Transfer Coefficients in Surface Layer with Different Momentum and Heat Roughness Lengths

A new method for computing the surface transfer coefficients is proposed, based on state-of-the-art empirical flux-profile relationships. The influence of the roughness length ratio is first demonstrated with the classical iterative calculation method. Then a non-iterative algorithm is developed, taking into account the difference between momentum and heat roughness lengths.The new method is validated by comparison with the reference iterative computation. The large gain-in computer processing time (CPU) time gain for the calculation of surface fluxes in Eulerian grid models is finally assessed.     

与ENSO相联系的热带太平洋典型风应力异常场结构

利用多年逐月海温距平和风应力距平观测资料,运用线性回归和EOF分析方法,分析了与ENSO相联系的热带太平洋典型风应力异常场结构。结果显示,与ENSO线性相关的风应力异常场在时间尺度上表现为低频变化,在水平结构上主要表现为四个典型分布。其中,分布型1主要表现为日界线以东赤道地区东风异常和赤道风应力场辐散;分布型2主要表现为日界线以东赤道地区西风异常和经向异常风应力向赤道气流;分布型3主要表现为日界线以东赤道偏南地区西风异常和风应力场辐合,日界线以西为东风异常;分布型4主要表现为160°W以东的弱东风异常和160°W以西的西风异常。与ENSO线性无关的风应力场主要表现为高频过程,在水平空间结构上其典型场主要位于赤道外地区。还用与ENSO有关的那部分回归风应力异常场强迫海洋距平模式,成功地再现了ENSO的主要信号。这表现观测揭示的典型风应力异常型对于El Nino的产生是根本性的。      

Heat Flux in the Coastal Zone

Various difficulties with application of Monin–Obukhov similarity theory are surveyed including the influence of growing waves, advection and internal boundary-layer development. These complications are normally important with offshore flow. The transfer coefficient for heat is computed from eddy correlation data taken at a mast two kilometres off the Danish coast in RASEX. For these coastal zone data, the thermal roughness length shows no well-defined relation to the momentum roughness length or roughness Reynolds number, in contrast to previous theories. The variation of the momentum roughness length is dominated by wave state. In contrast, the thermal roughness length shows significant dependence on wave state only for small values of wave age where the mixing is apparently enhanced by wave breaking. The development of thin internal boundary layers with offshore flow substantially reduces the heat transfer and thermal roughness length but has no obvious influence on momentum roughness length. A new formulation of the thermal roughness length based on the internal boundary-layer depth is calibrated to the RASEX data. For the very stable case, the turbulence is mainly detached from the surface and existing formulations do not apply.As an alternative to adjusting the thermal roughness length, the transfer coefficient is related directly to the stability and the internal boundary-layer depth. This avoids specification of roughness lengths resulting from the usual integration of the non-dimensional temperature function. The resulting stability function is simpler than previous ones and satisfies free convection similarity theory without introduction of the gustiness factor. The internal boundary layer also influences the moisture transfer coefficient.     

Vertical Structure Of Turbulence In Offshore Flow During Rasex

The adjustment of the boundary layer immediately downstream froma coastline is examined based on two levels of eddy correlation data collected on a mast at the shore and six levels of eddy correlation data and profiles of mean variables collected from a mast 2 km offshore during the Risø Air-Sea Experiment. The characteristics of offshore flow are studied in terms of case studies and inter-variable relationships for the entire one-month data set. A turbulent kinetic energy budget is constructed for each case study.The buoyancy generation of turbulence is small compared to shear generation and dissipation. However, weakly stable and weakly unstable cases exhibit completely different vertical structure. With flow of warm air from land over cooler water, modest buoyancy destruction of turbulence and reduced shear generation of turbulence over the less rough sea surface cause the turbulence to rapidly weaken downstream from the coast. The reduction of downward mixing of momentum by the stratification leads to smaller roughness lengths compared to the unstable case. Shear generation at higher levels and advection of stronger turbulence from land often lead to an increase of stress and turbulence energy with height and downward transport of turbulence energy toward the surface.With flow of cool air over a warmer sea surface, a convective internal boundary layer develops downstream from the coast. An overlying relatively thick layer of downward buoyancy flux (virtual temperature flux) is sometimes maintained by shear generation in the accelerating offshore flow.     

温度,盐度和风应力对南海海流模拟的影响

用美国普林斯顿大学海洋模式（ＰＯＭ）对南中国海的年平均海流进行了数值模拟,对温盐结构和风应力在海流形成中的作用进行了较详细的讨论。结果表明,仅有温盐水平不均匀分布也可以驱动海水而生成南海海流,但此种海流的结构较乱,最大流速只有３０～４０ｃｍ·ｓ－１。若温盐无水平结构,则在风应力驱动下,南海海流的结构较为有序,且最大流速可增至６０～７０ｃｍ·ｓ－１。在温盐水平分布不均匀并有风应力的作用时,生成的南海海流与仅有风应力作用时的海流场较相似,说明在南海海流的形成中,风应力的作用更为重要。海面自由高度的分析也证明了上述结论。