Wave-dependence of friction velocity,roughness length,and drag coefficient over coastal and open water surfaces by using three databases

The parameterization of friction velocity, roughness length, and the drag coefficient over coastal zones and open water surfaces enables us to better understand the physical processes of air-water interaction. In context of measurements from the Humidity Exchange over the Sea Main Experiment (HEXMAX), we recently proposed wave-parameter dependent approaches to sea surface friction velocity and the aerodynamic roughness by using the dimensional analysis method. To extend the application of these approaches to a range of natural surface conditions, the present study is to assess this approach by using both coastal shallow (RASEX) and open water surface measurements (Lake Ontario and Grand Banks ERS-1 SAR) where wind speeds were greater than 6.44 m s^-1. Friction velocities, the surface aerodynamic roughness, and the neutral drag coefficient estimated by these approaches under moderate wind conditions were compared with the measurements mentioned above. Results showed that the coefficients in these approaches for coastal shallow water surface differ from those for open water surfaces, and that the aerodynamic roughness length in terms of wave age or significant wave height should be treated differently for coastal shallow and open water surfaces.     

On the deep‐water fetch laws for wind‐generated surface gravity waves

Abstract

With the object of providing an accurate set of open‐sea wave spectra in a variety of conditions, we deployed, in conjunction with CASP, an array of 9 wave buoys (3 directional, 6 non‐directional) along a 30‐km line offshore from Martinique Beach, N.S. A large set of high‐quality wave spectra was collected in conjunction with extensive meteorological information. The data set is unique in the sense that a large onshore swell component was normally present.

Offshore‐wind cases for three windows: ±5°, ±15° and ±30° with respect to the shore normal, have been considered. Wind speed was found to be a strong function of fetch, and attempts were made to allow for this in the analysis. Power‐law regressions have been produced of dimensionless sea energy, peak frequency and high‐frequency spectral level (the Kitaigorodskii “alpha” parameter) vs dimensionless fetch and wind speed (inverse wave age). The regressions are compared with earlier work: the Joint North Sea Waves Project (Jonswap) and the Canada Centre for Inland Waters (CCIW) Lake Ontario study.

The comparisons indicate that dimensionless wave energies, peak frequencies and alpha values in this experiment are comparable with those from earlier experiments; in spite of different wind analysis methods, the CASP and CCIW fetch‐limited growth laws are consistent within the contexts of the two experiments. Differences among the estimated parameters are as large within the analyses of the three windows as they are among the three experiments we compare.     

TESTING DRAG COEFFICIENT APPROACHES BY USING THE BUOY DATA COLLECTED IN MODERATE TO HIGH WIND UNDER FOLLOWING, CROSSING AND OPPOSING SWELL CONDITIONS

Hurricane intensity and track are strongly affected by air-sea interactions. Classified as following swells, crossing swells, and opposing swells, the observed wave height was parameterized by using the 10-m wind speed collected on 5 buoys by the National Buoy Data Center during 13 hurricane events. The path information of these 13 hurricanes was obtained from the National Hurricane Center Best Track (NHC-BT). Results show that the wave height increases exponentially with the 10-m wind speed, and the wave height reaches the maximum value, 11.2 m (8.1 m), when 10-m wind speed is 40 m s-1 under the following and crossing (opposing) swell conditions. We find that the wave steepness (the ratio of wave height to wave length) is proportional to the -2/3 power of the wave age (the ratio of wave phase velocity to 10-m wind speed). The parameterizations of friction velocity and drag coefficient are tested using the buoy data collected in moderate to high wind under following, crossing and opposing swell conditions. A wave age dependent equation for drag coefficient is found more accurate and suggested for future usage in numerical models. Further, these algorithms also suggest that wind-swell orientation needs to be considered to retrieve accurate surface drag under high winds and strong swells.     

ANALYSIS OF SEA-SURFACE DRAG PARAMETERIZATIONS IN OPEN OCEAN CONDITIONS

Data from the Surface Waves and Processes Program (SWAPP) are employed to test current sea-surface drag parameterizations in open ocean conditions. General trends in the data indicate that drag increases with increasing wind speed and wave height, and decreases with wave age. However, scatter in the data limits the use of these parameters and other wave dependent parameterizations for modelling efforts. Upon close inspection, it is found that during the onset of three wind events analyzed separately, each of these parameters correlate well with the drag coefficient. However, the dependence of the drag coefficient on each of these parameters varies markedly from event to event. The disparity appears most closely linked to the turning rate of the wind, indicating that temporal and directional effects may play an important role. A temporal lag of O(4) hours between the rise of the wind and subsequent rise in the drag coefficient is also noticed, further pointing out the complexity of the wind-stress system.     

Airborne surveys of the atmospheric boundary layer above the marginal ice zone on the Newfoundland shelf

Abstract

Airborne measurements in the atmospheric boundary layer (ABL) above the marginal ice zone (MIZ) on the Newfoundland Shelf reveal strong lateral variations in mean wind, temperature and the vertical fluxes of heat and momentum under conditions of cold, off‐ice wind. Flux measurements in (and near) the surface layer indicate that the neutral 10‐m drag coefficient depends on ice concentration, ranging from 2 × 10^‐3 at 10% coverage to 5 × 10^‐3 at 90%. Furthermore, cross‐ice‐edge transects consistently show increasing wind speed, temperature and heat flux in the off‐ice direction, but the momentum flux may either increase or decrease, depending on the relative importance of surface buoyancy flux and roughness. For the conditions encountered in this experiment, it appears surface wave maturity does not have a significant influence on the drag coefficient in fetch‐limited regimes near the ice edge.     

Drag coefficient modeling for the near coastal zone

Using the JONSWAP spectrum for describing the surface wave state in the near coastal zone, models for the roughness length and the drag coefficient are used to simulate the dependence of the wind stress on fetch and depth. The results of each model are then compared with a compiled set of past investigations of the neutral drag coefficient over a variety of conditions. It is found that the models of Donelan, Hsu, and Kitaigorodskii correctly predict the trends in the drag coefficient with fetch and depth. Although it did not account for all the observed variations in the neutral drag coefficient. Kitaigorodskii's model, when incorporating the JONSWAP spectrum, more accurately simulated the slopes of the various C_DN regressions against windspeed.     

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.     

One-dimensional theory of the wave boundary layer

Results obtained in a 2-D modeling of the statistical structure of the wave boundary layer (WBL) are used for elaboration of the general approach to 1-D modeling taking into account the spectral properties of wave drag for an arbitrary wave field. In the case of the wave field described by the JONSWAP spectrum, the momentum and energy spectral density exchange, vertical profiles of the wave-induced momentum flux and dependence of total roughness parameter and drag coefficient on peak frequency are given. The reasons that the total roughness parameter increases with decreasing fetch are explained. The role of wind waves as an active element of the ocean-atmosphere dynamic system is also discussed.     

Wind profiles,momentum fluxes and roughness lengths at Cabauw revisited

We describe the results of an experiment focusing on wind speed and momentum fluxes in the atmospheric boundary layer up to 200 m. The measurements were conducted in 1996 at the Cabauw site in the Netherlands. Momentum fluxes are measured using the K-Gill Propeller Vane. Estimates of the roughness length are derived using various techniques from the wind speed and flux measurements, and the observed differences are explained by considering the source area of the meteorological parameters. A clear rough-to-smooth transition is found in the wind speed profiles at Cabauw. The internal boundary layer reaches the lowest k-vane (20 m) only in the south-west direction where the obstacle-free fetch is about 2 km. The internal boundary layer is also reflected in the roughness lengths derived from the wind speed profiles. The lower part of the profile (< 40 m) is not in equilibrium and no reliable roughness analysis can be given. The upper part of the profile can be linked to a large-scale roughness length. Roughness lengths derived from the horizontal wind speed variance and gustiness have large footprints and therefore represent a large-scale average roughness. The drag coefficient is more locally determined but still represents a large-scale roughness length when it is measured above the local internal boundary layer. The roughness length at inhomogeneous sites can therefore be determined best from drag coefficient measurements just above the local internal boundary layers directly, or indirectly from horizontal wind speed variance or gustiness. In addition, the momentum and heat fluxes along the tower are analysed and these show significant variation with height related to stability and possibly surface heterogeneity. It appears that the dimensionless wind speed gradients scale well with local fluxes for the variety of conditions considered, including the unstable cases.     

étude expérimentale de la couche limite de montagne

The characteristics of the boundary layer over complex terrain (Lannemezan - lat.: 43.7° N and, long.: 0.7 ° E) are analyzed for various scales, using measurements obtained during the COCAGNE Experiment. In this first part, the dynamic characteristics of the flow are studied with respect to atmospheric stability and the relief at small (~20 km) and medium scales (~100 km). These relief scales depend on the topographical profile of the Lannemezan Plateau along the dominant axis of the wind (E-W) and the Pyrénées Mountains located at the south of the experimental site. The terrain heterogeneities have a standard deviation of ~48 m and a wavelength of ~2 km.The averaged vertical profiles of wind speed and direction over the heterogeneous terrain are analyzed. The decrease of wind speed within the boundary layer is greater than over flat terrain (WANGARA Experiment). However, a comparison between ETTEX (complex terrain) and COCAGNE vertical wind speed profiles shows good agreement during unstable conditions. In contrast, during neutral conditions a more rapid increase with normalized height is found with COCAGNE than with ETTEX and WANGARA data. The vertical profiles of wind direction reveal an influence of the Pyrénées Mountains on the wind flow. The wind rotation in the BL is determined by the geostrophic wind direction-Pyrénées axis angle (negative deviation) as the geostrophic wind is connected with the Mountain axis.When the geostrophic wind does not interact with the Pyrénées axis, the mean and turbulent wind flow characteristics (drag coefficient C _D, friction velocity u _*) depend on the topography of the plateau. When the wind speed is strong (>6 m s ^-1), an internal boundary layer is generated from the leading edge of the Plateau.     

Equilibrium Atmospheric Boundary-Layer Flow: Computational Fluid Dynamics Simulation with Balanced Forces

Forcing relationships in steady, neutrally stratified atmospheric boundary-layer (ABL) flow are thoroughly analyzed. The ABL flow can be viewed as balanced between a forcing and a drag term. The drag term results from turbulent stress divergence, and above the ABL, both the drag and the forcing terms vanish. In computational wind engineering applications, the ABL flow is simulated not by directly specifying a forcing term in the ABL but by specifying boundary conditions for the simulation domain. Usually, these include the inflow boundary and the top boundary conditions. This ‘boundary-driven’ ABL flow is dynamically different from its real counterpart, and this is the major reason that the simulated boundary-driven ABL flow does not maintain horizontal homogeneity. Here, first a dynamical approach is proposed to develop a neutrally stratified equilibrium ABL flow. Computational fluid dynamics (CFD) software (Fluent 6.3) with the standard \(k\) – \(\varepsilon \) turbulence model is employed, and by applying a driving force profile, steady equilibrium ABL flows are simulated by the model. Profiles of wind speed and turbulent kinetic energy (TKE) derived using this approach are reasonable in comparison with the conventional logarithmic law and with observational data respectively. Secondly, the equilibrium ABL profiles apply as inflow conditions to simulate the boundary-driven ABL flow. Simulated properties between the inlet and the outlet sections across a fetch of 10 km are compared. Although profiles of wind speed, TKE, and its dissipation rate are consistently satisfactory under higher wind conditions, a deviation of TKE and its dissipation rate between the inlet and outlet are apparent (7–8 %) under lower wind-speed conditions (2 m s \(^{-1}\) at 10 m). Furthermore, the simulated surface stress systematically decreases in the downwind direction. A redistribution of the pressure field is also found in the simulation domain, which provides a different driving pattern from the realistic case in the ABL.     

Wind‐induced microseisms from Lake Ontario

Abstract

The characteristics of microseisms measured in four vaults of the Southern Ontario Seismic Network within 30 km of the shore of Lake Ontario are analysed. It is shown that the rms values in the 1–3 Hz band are coherent between the stations, indicating a common generative mechanism. A distinct onshore intermittent flux of Rayleigh‐like wave energy was detected at a site near the shore. Microseismic energy in this band is distinctly correlated with the wind speed. The incremental microseismic energy above an absolute minimum activity as a function of wind direction, for a given fixed wind speed, correlates with the average fetch of the wind over the lake, indicating that the source of microseisms is the lake itself. The sensitivity to fetch effects is similar for both onshore and offshore stations indicating that shoaling is probably not a source. Niagara Falls, which also can have a wind‐dependent flow from Lake Erie, causes a measurable effect to at least 25 km but does not significantly affect stations at a distance of 150 km.     

New evidence for a relation between wind stress and wave age from measurements during ASGAMAGE

Data from the 1996 ASGAMAGE experiment, performed in the southern North Sea at research platform Meetpost Noordwijk (MPN), are analysed for the parameters affecting the momentum flux. The stress turns out to be quadratically related to the 10-m wind speed and linearly to the wind speed at a wavelength related level. The Charnock parameter (dimensionless roughness length) shows a pronounced correlation with wave age. This implies, due to a coupling between wave age and the steepness of the waves, a connection between the stress and the steepness. We find that our North Sea results are consistent withopen ocean observations. For a given wind speed the mean stress at MPN turns out to be higher because the wave age there is in general lower. We define and give an expression for a drag coefficient at a wavelength related level that can be calculated straightforwardly from the wave age and then reduced to a standard level.     

南海北部风浪特征的计算

探讨南海北部海域风浪成长时有效波高与风速、风时、风区之间的关系,同时分析了5种风浪要素的推算方法,探讨其在南海北部海域的适用性。结果表明：1）在南海北部,风速和风时呈现线性增长的关系,风速越大,风浪从过度状态成长到充分成长状态所需风时就越长;风速大小和风区长度之间满足平方关系,风速越大,风浪充分成长所需风区长度就越长。2）在南海北部,有效波高的大小与风速的大小、风时的长短、风区的长度3者密切相关。3）SMB方法、W ilson IV方法和青岛方法,在计算南海北部的风浪关系中体现出了一定的稳定性和适应性。     

Drag of the sea surface

It is shown how the drag of the sea surface can be computed from the wind speed and the sea state. The approach, applicable both for fully developed and for developing seas, is based on conservation of momentum in the boundary layer above the sea, which allows one to relate the drag to the properties of the momentum exchange between the sea waves and the atmosphere.The total stress is split into two parts: a turbulent part and a wave-induced part. The former is parameterized in terms of mixing-length theory. The latter is calculated by integration of the wave-induced stress over all wave numbers. Usually, the effective roughness is given in terms of the empirical Charnock relation. Here, it is shown how this relation can be derived from the dynamical balance between turbulent and wave-induced stress. To this end, the non-slip boundary conditions is assigned to the wave surface, and the local roughness parameter is determined by the scale of the molecular sublayer.The formation of the sea drag is then described for fully developed and developing seas and for light to high winds.For the Charnock constant, a value of about 0.018–0.030 is obtained, depending on the wind input, which is well within the range of experimental data.It is shown that gravity-capillary waves with a wavelength less than 5 cm play a minor role in the momentum transfer from wind to waves. Most of the momentum is transferred to decimeter and meter waves, so that the drag of developing seas depends crucially on the form of the wave spectrum in the corresponding high wavenumber range.The dependence of the drag on wave age depends sensitively on the dependence of this high wavenumbertail on wave age. If the tail is wave-age independent, the sea drag appears to be virtually independent of wave age. If the tail depends on wave age, the drag also does. There is contradictory evidence as to the actual dependence. Therefore, additional experiments are needed.The investigation was in part supported by the Netherlands Geosciences Foundation (GOA) with financial aid from the Netherlands Organization for Scientific Research (NWO).     

Cross-Spectra Over the Sea from Observations and Mesoscale Modelling

Cospectra and quadrature spectra are calculated for six pairs of tall offshore measurement masts near the Horns Rev I wind farm in the Danish North Sea and the Nysted wind farm in the Baltic sea. The mast-pairs are separated from one another by horizontal distances of 2.13–12.4 km. Cospectra and quadrature spectra for the two sites are classified in terms of the angle between the mean wind direction and the line connecting each pair of masts. The frequency axes of the spectra are normalized to remove the effect of mean wind speed and separation distance. Results indicate a larger contribution to the quadrature spectrum for flow from the sea than for flow from the land, and the patterns in the spectra are clearer and better defined for Horns Rev I (which has a long uninterrupted sea-fetch from the west) than for Nysted (which is surrounded by a more complicated coastline). The analysis is replicated based on 3-month simulations using the weather research and forecasting (WRF) numerical model with a horizontal grid spacing of 2 km. For the sea-fetch directions, good agreement in spectral properties between the model and observations is found. Analytical expressions based on the properties of the cross-correlation function and an exponentially decaying coherence function are fitted to the normalized cospectra and quadrature spectra. The expressions are shown to be a good fit to the spectra calculated from the WRF simulations and to the observed spectra for directions with a long sea-fetch, which suggests that to a good approximation, the average cospectra and quadrature spectra over the sea can be written as functions of frequency, mean wind speed, separation distance and the angle between the wind direction and the orientation of the masts.     

Impact of waves on the sea drag: measurements in the Baltic sea and a model interpretation

Measurements from the Baltic Sea and a wind-over-wave coupled model are used to study the wave impact on the sea drag. The study has been carried out for different wave conditions, namely a pure wind-sea, following-swell/ mixed sea and cross-swell/ mixed sea. Measurements reveal the fact that the sea drag is dependent on the sea-state. In stationary conditions and in the absence of severe cross-swell, swell reduces drag compared to wind-sea at the same wind speed. The cross-swell enhances the drag as compared to the following-swell case and the magnitude of the drag coefficient is increased with increasing the angle of swell propagation to the wind. It is shown that the agreement between the model results and measurements is good for pure wind-sea and stationary mixed-sea cases. Discrepancies occur at light winds, where most of the data represent pure swell conditions. During these pure swell conditions the data are characterized by a large variation of the drag coefficient. The variation is caused by mesoscale variability in the stress co-spectra, wind-cross-swell effects and nonstationarity in the wave and wind fields not represented in the model.     

Momentum and sensible heat fluxes calculated by the dissipation technique during the Ocean Storms Project

High frequency measurements of wind velocity and temperature were made during the Ocean Storms Project in November 1987. The dissipation method was applied to the resulting time series in order to determine friction velocities,u _*, and the characteristic temperature scale,t _*, at 1-min intervals. These values were then compared to the 1-min mean wind speed and air-sea temperature differences to determine relationships for the drag coefficient (C _d ) and Stanton number (C _h ). The drag coefficient was comparable to other values reported in the literature, although the variation with wind speed was greater than reported by other investigators. An examination of the residual time series indicated a systematic low frequency periodicity of about 2-hr duration which was attributed to a fluctuating wind interacting with the surface gravity wave field. The temperature fluctuations did not produce meaningful estimates ofC _h for stable conditions. For unstable conditions, a value of 1.09±0.02×10^–3 was found.     

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.     

On the growth rate of wind‐generated waves

Abstract

A new approach to fetch‐limited wave studies is taken in this paper. Using data from five towers arranged along a line from the eastern shore of Lake St Clair, the differential growth between towers is explored as a function of local wave age. It is argued that this method avoids the usual fetch‐limited pitfall of inhomogeneity over long fetches and, in particular, the changes in wind speed downfetch of an abrupt roughness change. It is found that the growth rate decreases uniformly downfetch as the waves approach full development. This differential method leads to a smooth transition from rapidly growing short fetch waves to the asymptotic invariant state of full development. When the variation in wind speed after an abrupt (land to water) roughness change is taken into account, the idea of a universal fetch‐limited growth curve is called into question.