Stability dependence of the drag and bulk transfer coefficients over a coastal sea surface

Bulk transfer coefficients were evaluated from eddy correlation flux measurements on a fixed pier during onshore winds. The mean values are C _D = 1.69 × 10^-3, C _H = 2.58 × 10^-3 and C _E = 1.51 × 10^-3. The drag coefficient, C _D, gradually increases with wind speed but C _H and C _E are independent of wind speed. According to theory and empirical formulas based on experimental results over flat grassland, the transfer coefficients should gradually increase with increasing instability. This is confirmed experimentally in the stable region in our case. However, the drag coefficient appears to decrease with increasing instability, which is against the theoretical result. A stability dependence is not clearly observed for C _H or C _E.     

The scintillation method tested over a dry vineyard area

Measurements of a scintillometer device mounted at 4 m above a dry vineyard area in La Mancha, Spain, are used to obtain an average sensible heat flux densityH. Averaging is over a rectangular area determined by the distance between the scintillometer light source and receptor (875 m) and some upwind distance governed by the horizontal wind speed perpendicular to that line. Using similarity relations obtained from La Crau, a good correspondence betweenH measured with the scintillometer and an eddy-correlation device in the centre of a vineyard is obtained. The friction velocityu _* was either measured directly using a sonic anemometer or obtained indirectly from two wind speeds and known values of the roughness length z_o and displacementd. The free convection formulation underestimates the sensible heat flux by about 30%. This is due to a significant contribution of mechanically generated turbulence to the total turbulent transport, which was caused by relatively strong winds and rough terrain.     

Surface inhomogeneity effects on convective diffusion

It is suggested that convective scaling, with appropriate extensions, provides the most useful framework for estimating the effects of urban-scale surface inhomogeneities on diffusion in convective conditions. Strong contrasts in surface heat flux exist between cropland, forests, urban areas, and water or marshland surfaces. It is argued that a typical fetch for convective turbulence to readjust to changed heat (or buoyancy) input from the surface below is 2(U/w ^*)h, where U is the mean wind speed in the mixing layer, w ^* is the convective scaling velocity, and h is the mixing depth. In contrast, the fetch required for wind speed to readjust to new underlying surface roughness is of the order (U/u ^*)²h/2, where u ^* is the friction velocity.The ratio w ^*/U is the best index of diffusion rates in moderately to very unstable conditions. General urban effects on heat flux, h, and U are discussed separately, then their combined effects on w ^*/U are estimated. While this ratio can double over a large city during light winds, its increase is much less for small cities, or during moderate winds. Finally, some examples of heat flux in- homogeneities causing stationary convective features are presented. Steady downdrafts associated with these features are of the order of 0.4w ^*, and could significantly increase surface concentrations from elevated sources.On assignment from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce.This paper is based on a presentation made at the AMS Specialty Conference on Air Quality Modeling of the Urban Boundary Layer, in Baltimore, late 1983.     

A theory for the scalar roughness and the scalar transfer coefficients over snow and sea ice

Although the bulk aerodynamic transfer coefficients for sensible (C _H ) and latent (C _E ) heat over snow and sea ice surfaces are necessary for accurately modeling the surface energy budget, they have been measured rarely. This paper, therefore, presents a theoretical model that predicts neutral-stability values of C _H and C _E as functions of the wind speed and a surface roughness parameter. The crux of the model is establishing the interfacial sublayer profiles of the scalars, temperature and water vapor, over aerodynamically smooth and rough surfaces on the basis of a surface-renewal model in which turbulent eddies continually scour the surface, transferring scalar contaminants across the interface by molecular diffusion. Matching these interfacial sublayer profiles with the semi-logarithmic inertial sublayer profiles yields the roughness lengths for temperature and water vapor. When coupled with a model for the drag coefficient over snow and sea ice based on actual measurements, these roughness lengths lead to the transfer coefficients. C _E is always a few percent larger than CH. Both decrease monotonically with increasing wind speed for speeds above 1 m s^–1, and both increase at all wind speeds as the surface gets rougher. Both, nevertheless, are almost always between 1.0 × 10^–3 and 1.5 × 10^–3.     

Tilt errors and precipitation effects in trivane measurements of turbulent fluxes over open water

For 390 ten-minute samples of turbulent flux, made with a trivane above a lake, the vertical alignment is determined within 0.1 ° through azimuth-dependent averaging. One degree of instrumental misalignment is found to produce an average tilt error of 9 ± 4% for momentum flux, and 4 ± 2% for heat flux. The tilt error in the vertical momentum flux depends mainly ons _u/u_*, and cannot be much diminished with impunity by high-pass pre-filtering of the turbulence signals. The effects of rain on trivane measurements of vertical velocity are shown to be negligible at high wind speeds, and adaptable to correction in any case.The normalized vertical velocity variance,s _w/u_*, appears to be proportional to the square root ofz/L for unstable stratification. For a wind speed range of 2 to 15 m s^–1, the eddy correlation stresses measured at 4- and 8-m heights can be reasonably well estimated by using a constant drag coefficientC _d=1.3 X 10^-3, while cup anemometer profile measurements give an overestimate of eddy stress at high wind speeds. A good stress estimate is also obtained from the elevation variance; it is suggested that trivane measurement of this variance might be made from a mobile platform, e.g., a moderately stabilized spar buoy.     

Large‐eddy simulation of small‐scale surface effects on the convective boundary‐layer structure

Abstract

The effects of small‐scale surface inhomogeneities on the turbulence structure in the convective boundary layer are investigated using a high‐resolution large‐eddy simulation model. Surface heat flux variations are sinusoidal and two‐dimensional, dividing the total domain into a checkerboard‐like pattern of surface hot spots with a 500‐m wavelength in the x and y directions, or 1/4 of the domain size. The selected wind speeds were 1 and 4 m s^‐l, respectively. As a comparison, a simulation of the turbulence structure was performed over a homogeneous surface.

When the wind speed is light, surface heat flux variations influence the horizontally averaged turbulence statistics, including the higher moments despite the small characteristic length of the surface perturbation. Stronger mean wind speeds weaken the effects of inhomogeneous surface conditions on the turbulence structure in the convective boundary layer.

Results from conditional sampling show that when the mean wind speed is small, weak mean circulations occur, with updraft branches above the high heat flux regions and down‐draft branches above the low heat flux regions. The inhomogeneous surface induces significant differences in the turbulence statistics between the high and low heat flux regions. However, the effect of the surface perturbations weaken rapidly when the mean wind speed increases. This research has implications in the explanation of the large‐scale variability commonly encountered in aircraft observations of atmospheric turbulence.     

Numerical simulation of turbulent convective flow over wavy terrain

By means of a large-eddy simulation, the convective boundary layer is investigated for flows over wavy terrain. The lower surface varies sinusoidally in the downstream direction while remaining constant in the other. Several cases are considered with amplitude up to 0.15H and wavelength ofH to 8H, whereH is the mean fluid-layer height. At the lower surface, the vertical heat flux is prescribed to be constant and the momentum flux is determined locally from the Monin-Obukhov relationship with a roughness lengthz _o=10^–4 H. The mean wind is varied between zero and 5w _*, wherew _* is the convective velocity scale. After rather long times, the flow structure shows horizontal scales up to 4H, with a pattern similar to that over flat surfaces at corresponding shear friction. Weak mean wind destroys regular spatial structures induced by the surface undulation at zero mean wind. The surface heating suppresses mean-flow recirculation-regions even for steep surface waves. Short surface waves cause strong drag due to hydrostatic and dynamic pressure forces in addition to frictional drag. The pressure drag increases slowly with the mean velocity, and strongly with /H. The turbulence variances increase mainly in the lower half of the mixed layer forU/w _*>2.     

BOUNDARY LAYER STRUCTURE IN TYPHOON SAOMAI (2006): UNDERSTANDING THE EFFECTS OF EXCHANGE COEFFICIENT

Recent studies have shown that surface fluxes and exchange coefficients are particularly important to models attempting to simulate the evolution and maintenance of hurricanes or typhoons.By using an advanced research version of the Weather Research and Forecasting(ARW)modeling system,this work aims to study the impact of modified exchange coefficient on the intensity and structures of typhoon Saomai(2006)over the western North Pacific.Numerical experiments with the modified and unmodified exchange coefficients are used to investigate the intensity and structure of the storm,especially the structures of the boundary layer within the storm.Results show that,compared to the unmodified experiment,the simulated typhoon in the modified experiment has a bigger deepening rate after 30-h and is the same as the observation in the last 12-h.The roughness is leveled off when wind speed is greater than 30 m/s.The momentum exchange coefficient(CD)and enthalpy exchange coefficient(CK)are leveled off too,and CD is decreased more than CK when wind speed is greater than 30 m/s.More sensible heat flux and less latent heat flux are produced.In the lower level,the modified experiment has slightly stronger outflow,stronger vertical gradient of equivalent potential temperature and substantially higher maximum tangential winds than the unmodified experiment has.The modified experiment generates larger wind speed and water vapor tendencies and transports more air of high equivalent potential temperature to the eyewall in the boundary layer.It induces more and strong convection in the eyewall,thereby leading to a stronger storm.     

On the Determination of the Neutral Drag Coefficient in the Convective Boundary Layer

Based on the idea that free convection can be considered as a particular case of forced convection, where the gusts driven by the large-scale eddies are scaled with the Deardorff convective velocity scale, a new formulation for the neutral drag coefficient, C_Dn, in the convective boundary layer (CBL) is derived. It is shown that (i) a concept of C_Dn can still be used under strongly unstable conditions including a pure free-convection regime even when no logarithmic portion in the velocity profile exists; (ii) gustiness corrections must be applied for rational calculations of C_Dn; and (iii) the stratification function used in the derivation of C_Dn should satisfy the theoretical free-convection limit. The new formulation is compared with the traditional relationship for C_Dn, and data collected over the sea (during the Tropical Ocean-Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) and the San Clemente Ocean Probing Experiment (SCOPE)) and over land (during the BOREX-95 experiment) are used to illustrate the difference between the new and traditional formulations. Compared to the new approach, the traditional formulation strongly overestimates C_Dn and z_o in the CBL for mean wind speed less than about 2 m s^-1. The new approach also clarifies several contradictory results from earlier works. Some aspects related to an alternate definition of the neutral drag coefficient and the wind speed and the stress averaging procedure are considered.     

The Signature of Sea Spray in the Hexos Turbulent Heat Flux Data

The role of sea spray intransferring heat and moisture across the air-sea interface has remained elusive. Some studies have reported that sea spray does not affect the turbulent air-sea heat fluxes for 10-m wind speeds up to at least 25 m s^-1, while others have reported important spray contributions for wind speeds as low as 12 m s^-1. One goal of the HEXOS (Humidity Exchange over the Sea) program was to quantify spray's contribution to the turbulent air-sea heat fluxes, but original analyses of the HEXOS flux data found the spray signal to be too small to be reliably identified amid the scatter in the data. We look at the HEXOS data again in the context of the TOGA-COARE bulk flux algorithm and a sophisticated microphysical spray model. This combination of quality data andstate-of-the-art modelling reveals a distinct spray signature in virtually all HEXOS turbulent heat flux data collected in winds of 15 m s^-1 and higher. Spray effects are most evident in the latent heat flux data, where spray contributes roughly 10% of the total turbulent flux in winds of 10 m s^-1 and between 10 and 40% in winds of 15–18 m s^-1. The spray contribution to the total sensible heat flux is also at least 10% in winds above 15 m s^-1. These results lead to a new, unified parameterization for the turbulent air-sea heat fluxes that should be especially useful in high winds because it acknowledges both the interfacial and spray routes by which the sea exchanges heat and moisture with the atmosphere.     

Generation and atmospheric heat exchange of coastal polynyas in the Weddell Sea

The forcing mechanisms for Antarctic coastal polynyas and the thermodynamic effects of existing polynyas are studied by means of an air-sea-ice interaction experiment in the Weddell Sea in October and November 1986.Coastal polynyas develop in close relationship to the ice motion and form most rapidly with offshore ice motion. Narrow polynyas occur frequently on the lee side of headlands and with strong curvature of the coastline. From the momentum balance of drifting sea ice, a forcing diagram is constructed, which relates ice motion to the surface-layer wind vector v _z and to the geostrophic ocean current vector c _g . In agreement with the data, wind forcing dominates when the wind speed at a height of 3 m exceeds the geostrophic current velocity by a factor of at least 33. This condition within the ocean regime of the Antarctic coastal current usually is fulfilled for wind speeds above 5 m/s at a height of 3 m.Based on a nonlinear parameter estimation technique, optimum parameters for free ice drift are calculated. Including a drift dependent geostrophic current in the ice/water drag yields a maximum of explained variance (91%) of ice velocity.The turbulent heat exchange between sea ice and polynya surfaces is derived from surface-layer wind and temperature data, from temperature changes of the air mass along its trajectory and from an application of the resistance laws for the atmospheric PBL. The turbulent heat flux averaged over all randomly distributed observations in coastal polynyas is 143 W/m². This value is significantly different over pack ice and shelf ice surfaces, where downward fluxes prevail. The large variances of turbulent fluxes can be explained by variable wind speeds and air temperatures. The heat fluxes are also affected by cloud feedback processes and vary in time due to the formation of new ice at the polynya surface.Maximum turbulent fluxes of more than 400 W/m² result from strong winds and low air temperatures. The heat exchange is similarly intense in a narrow zone close to the ice front, when under weak wind conditions, a local circulation develops and cold air associated with strong surface inversions over the shelf ice is heated above the open water.     

Models for estimating offshore winds from onshore meteorological measurements

To understand and estimate wind speed differences across the coastal zone, two models, one theoretical and another semi-empirical, have been developed and verified by available data sets. Assuming that: (1) mean horizontal motion exists across the coastal zone; and (2) the geostrophic wind does not change appreciably at the top of the planetary boundary layer (PBL), the equation of motion in the direction of the wind can be reduced so that 341-01, where U, H, and C _D are wind speed, height of PBL, and drag coefficient over the sea and land, respectively. For practice, C _D SEA has been modified from a formula with U _LAND as the only input. H _SEA may be estimated routinely from known H _{D LAND LAND} and the temperature difference between land and sea, which can be provided by such means as remote sensing from meteorological satellites. For a given coast, Cmay be estimated also. This formula is recommended for weather forecasters. The semiempirical formula is based mainly on the power law wind distribution with height in the PBL. The formula states that 341-02. Simultaneous offshore and onshore wind measurements made at stations ranging from Somalia, near the equator, to the Gulf of Alaska indicated that values of a and b are 2.98 and 0.34 with a correlation coefficient of -0.95. For oceanographic applications, a simplified equation, i.e., 341-03, is also proposed.     

The sea surface temperature deviation and the heat flow at the sea-air interface

The deviation of the sea surface temperature from the water temperature below is calculated as a function of the heat flow through the air-sea interface, using wind tunnel determinations of the effective thermal diffusivity in a boundary layer. The influence ofQ, shortwave radiation, andH, latent and sensible heat transfer plus effective back radiation, and U, wind speed, can be described by:T ₀ –T _w =C ₁ ·H/U +C ₂ ·Q/U. The calculated coefficients vary slightly with reference depth, Tables II and III. They are in good agreement with independent observations.On leave at Department of Oceanography, Oregon State University, Corvallis, Oregon in 1969–70.     

Flux Footprints in the Convective Boundary Layer: Large-Eddy Simulation and Lagrangian Stochastic Modelling

We investigated the flux footprints of receptors at different heights in the convective boundary layer (CBL). The footprints were derived using a forward Lagrangian stochastic (LS) method coupled with the turbulent fields from a large-eddy simulation model. Crosswind-integrated flux footprints shown as a function of upstream distances and sensor heights in the CBL were derived and compared using two LS particle simulation methods: an instantaneous area release and a crosswind linear continuous release. We found that for almost all sensor heights in the CBL, a major positive flux footprint zone was located close to the sensor upstream, while a weak negative footprint zone was located further upstream, with the transition band in non-dimensional upwind distances −X between approximately 1.5 and 2.0. Two-dimensional (2D) flux footprints for a point sensor were also simulated. For a sensor height of 0.158 z _i, where z _i is the CBL depth, we found that a major positive flux footprint zone followed a weak negative zone in the upstream direction. Two even weaker positive zones were also present on either side of the footprint axis, where the latter was rotated slightly from the geostrophic wind direction. Using CBL scaling, the 2D footprint result was normalized to show the source areas and was applied to real parameters obtained using aircraft-based measurements. With a mean wind speed in the CBL of U = 5.1 m s⁻¹, convective velocity of w _* = 1.37 m s⁻¹, CBL depth of z _i = 1,000 m, and flight track height of 159 m above the surface, the total flux footprint contribution zone was estimated to range from about 0.1 to 4.5 km upstream, in the case where the wind was perpendicular to the flight track. When the wind was parallel to the flight track, the total footprint contribution zone covered approximately 0.5 km on one side and 0.8 km on the other side of the flight track.     

Spatio-temporal Variability of Surface-layer Turbulent Fluxes Over the Bay of Bengal and Arabian sea During the ICARB Field Experiment

We report the spatio-temporal variability of surface-layer turbulent fluxes of heat, moisture and momentum over the Bay of Bengal (BoB) and the Arabian Sea (AS) during the Integrated Campaign for Aerosols, gases Radiation Budget (ICARB) field experiment. The meteorological component of ICARB conducted during March – May 2006 onboard the oceanic research vessel Sagar Kanya forms the database for the present study. The bulk transfer coefficients and the surface-layer fluxes are estimated using a modified bulk aerodynamic method, and then the spatio-temporal variability of these air-sea interface fluxes is discussed in detail. It is observed that the sensible and latent heat fluxes over the AS are marginally higher than those over the BoB, which we attribute to differences in the prevailing meteorological conditions over the two oceanic regions. The values of the wind stress, sensible and latent heat fluxes are compared with those obtained for the Indian Ocean Experiment (INDOEX) period. The variation of drag coefficient (C _D ), exchange coefficients of sensible heat and moisture (C _H = C _E ) and neutral drag coefficient (C _DN ) with wind speed is also discussed.

---

     

Characteristic heat transfer coefficients near the ground at night during strong winds

Based on theoretical radiative cooling values and observed temperature changes with time near the surface during the night, the bulk heat transfer coefficient C _H is estimated from standard meteorological observations obtained from stations representative of open rural, small town and large urban areas, for nights with clear skies and relatively strong winds. It is shown that C _H is smaller than the drag coefficient C _M, and that C _H/C _M over urban areas is smaller than that over open countryside.     

Air–Sea Interaction Features in the Baltic Sea and at a Pacific Trade-Wind Site: An Inter-comparison Study

A systematic comparison of wind profiles and momentum exchange at a trade wind site outside Oahu, Hawaii and corresponding data from the Baltic Sea is presented. The trade wind data are to a very high degree swell dominated, whereas the Baltic Sea data include a more varied assortment of wave conditions, ranging from a pure growing sea to swell. In the trade wind region swell waves travel predominantly in the wind direction, while in the Baltic, significant cross-wind swells are also present. Showing the drag coefficient as a function of the 10-m wind speed demonstrates striking differences for unstable conditions with swell for the wind-speed range 2 m s^?1 < U ₁₀ < 7 m s^?1, where the trade-wind site drag values are significantly larger than the corresponding Baltic Sea values. In striking contrast to this disagreement, other features studied are surprisingly similar between the two sites. Thus, exactly as found previously in Baltic Sea studies during unstable conditions and swell, the wind profile in light winds (3 m s^?1) shows a wind maximum at around 7–8 m above the water, with close to constant wind speed above. Also, for slightly higher wind speeds (4 m s^?1 < U ₁₀ < 7 m s^?1), the similarity between wind profiles is striking, with a strong wind-speed increase below a height of about 7–8 m followed by a layer of virtually constant wind speed above. A consequence of these wind-profile features is that Monin–Obukhov similarity is no longer valid. At the trade-wind site this was observed to be the case even for wind speeds as high as 10 m s^?1. The turbulence kinetic energy budget was evaluated for four cases of 8–16 30- min periods at the trade-wind site, giving results that agree very well with corresponding figures from the Baltic Sea.     

Generation mechanisms of convectively induced internal gravity waves in a three-dimensional framework

The generation mechanisms of convective gravity waves in the stratosphere are investigated in a three-dimensional framework by conducting numerical simulations of four ideal storms under different environmental conditions: one un-sheared and three constant low-level sheared basic-state winds with the depth of the shear layer of 6 km and the surface wind speeds (U_s) of 8, 18, and 28 m s^?1, using the Advanced Regional Prediction System (ARPS) model. The storms simulated under the un-sheared (U_s = 0 m s^?1), weakly sheared (U_s = 8 and 18ms^?1), and strongly sheared (U_s = 28ms^?1) basicstate winds are classified into single-cell, multicell, and supercell storms, respectively. For each storm, the wave perturbations in a control simulation, including nonlinearity and microphysical processes, are compared with those in quasi-linear dry simulations forced by diabatic forcing and nonlinear forcing that are obtained from the control simulation. The gravity waves generated by the two forcing terms in the quasi-linear dry simulations are out of phase with each other for all of the storms. The gravity waves in the control simulation are represented by a linear sum of the wave perturbations generated by the nonlinear forcing and diabatic forcing. This result is consistent with the results of previous studies in a two-dimensional framework. This implies that both forcing mechanisms are important for generating the convective gravity waves in the three-dimensional framework as well. The characteristics of the three-dimensional gravity waves in the stratosphere were determined by the spectral combination of the forcing terms and the wave-filtering and resonance factor that is determined from the basic-state wind and stability as well as the vertical structure of the forcing.     

Aerodynamic roughness of the sea surface at high winds

The role of the surface roughness in the formation of the aerodynamic friction of the water surface at high wind speeds is investigated. The study is based on a wind-over-waves coupling theory. In this theory waves provide the surface friction velocity through the form drag, while the energy input from the wind to waves depends on the friction velocity and the wind speed. The wind-over-waves coupling model is extended to high wind speeds taking into account the effect of sheltering of the short wind waves by the air-flow separation from breaking crests of longer waves. It is suggested that the momentum and energy flux from the wind to short waves locally vanishes if they are trapped into the separation bubble of breaking longer waves. At short fetches, typical for laboratory conditions, and strong winds the steep dominant wind waves break frequently and provide the major part of the total form drag through the air-flow separation from breaking crests, and the effect of short waves on the sea drag is suppressed. In this case the dependence of the drag coefficient on the wind speed is much weaker than would be expected from the standard parameterization of the roughness parameter through the Charnock relation. At long fetches, typical for the field, waves in the spectral peak break rarely and their contribution to the air-flow separation is weak. In this case the surface form drag is determined predominantly by the air-flow separation from breaking of the equilibrium range waves. As found at high wind speeds up to 60 m s⁻¹ the modelled aerodynamic roughness is consistent with the Charnock relation, i.e. there is no saturation of the sea drag. Unlike the aerodynamic roughness, the geometrical surface roughness (height of short waves) could be saturated or even suppressed when the wind speed exceeds 30 m s⁻¹.     

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.