A planetary boundary-layer numerical model containing third-order derivative terms

In this paper, the third-order derivative of velocity with respect to height is included in the traditional motion equations of the neutral PBL. The nonlinear equations are solved numerically to obtain the vertical distribution of wind in the PBL and some PBL characteristic parameters. Reasonable simulations of the Leipzig wind profile using these parameters show the success of this kind of nonlocal closure in a real PBL simulation.     

The Value Of The β Coefficient In The Relaxed Eddy Accumulation Method In Terms Of Fourth-Order Moments

The commonly measured value of in the relaxed eddy accumulationmethod of about 0.56is shown to arise from the non-Gaussiannature of turbulence. Fourth-orderGram–Charlier functions forthe two-dimensional probability distributionsof variation in the horizontal component of wind velocityand concentrations of water vapour, carbondioxide and methane with respect to thevertical component of wind velocity are used to examinethe value of .An analytical solution for ispresented in terms of fourth-order moments.Under mean conditions, this solution givesa value for of0.557. Variation of is shown to be controlledprimarily by the ratio of the mean ofc'w3 (where c'is relevant to the entity of interest andw' is vertical component of windvelocity) to the correlationcoefficient between the entity concentrationand vertical component of wind velocity.     

Wind Stress Structure in the Unstable Marine Surface Layer Detected by Sar

The wind stress in the marine surface layer under unstable conditions and low wind speed has been studied using a Synthetic Aperture Radar (SAR) image of the sea surface and time series of the horizontal and vertical wind velocities and of the wind stress recorded on board the C.N.R. research platform, in the northern Adriatic Sea, during a SAR overflight.A conditional sampling technique has been used on the wind stress time series and on the SAR image to detect downward (sweep) and upward (ejection) bursts of the momentum flux, as well as the two-dimensional structure of the radar backscatter.From the ensemble average of both the wind stress and the backscatter structures, it has been possible to estimate the mean duration of the upward (11 s) and the downward (15 s) wind stress bursts and the mean size of the bright patches of the SAR image (120 m). The front of the mean backscatter structure, associated with the downward wind stress bursts, has been related to the time length of the mean sweep stress structure to get, after accounting for a threshold of the wind stress for the generation of the sea surface wavelets, the translation velocity U_t of the mean wind stress of sweep, very close to the mean wind speed. The vertical coherence of the wind stress structures has permited to refer the translation velocity to a level very close to the sea surface, but above the viscous sublayer. The variability of U_t with height has been studied through comparison with the mean wind speed at different heights z calculated by a boundary-layer model. Accounting for the results reported in the literature, there is an indication that U_t is constant with height in the range 0.5 m z 15 m.The two-dimensional pattern of the wind stress structures has been derived from the SAR image. The structures appear elongated crosswind, as with microfronts, with an average cross- to down-wind ratio of 4. The area covered by the downward wind stress structures represents 13% of the total area.     

A study of a fair weather boundary layer in TOGA-COARE: Parameterization of surface fluxes in large scale and regional models for light wind conditions

The study focuses on a way to parameterize the effect of subgrid scale convective motions on surface fluxes in large scale and regional models for the case of light surface winds. As previously proposed, these subgrid effects are assumed to scale with the convection intensity through the relationship: where is the mean velocity of the wind, U₀ the velocity of the mean wind, w_* the free convection velocity, and an empirical coefficient to be determined. Both observations and numerical simulation are presently used to determine the free convection coefficient .Large eddy simulation of a fair weather convective boundary layer case observed during TOGA-COARE is performed. Comparisons between observations and the simulation of surface properties and vertical profiles in the planetary boundary layer are presented. The simulated vertical turbulent fluxes of heat, moisture and buoyancy range well within estimates from aircraft measurements.The most important result is that the true free convection coefficient , directly estimated from simulation, leads to a value of 0.65, smaller than the ones estimated from temporal and spatial variances. Using observations and simulation, estimates of from temporal and spatial variances are obtained with similar values 0.8. From both theoretical derivations and numerical computations, it is shown that estimates of the true from variances are possible but only after applying a correction factor equal to 0.8. If this correction is not used, is overestimated by about 25%. The time and space sampling problem is also addressed in using numerical simulations.     

On the turbulence structure in the stable boundary layer over the Greenland ice sheet

Turbulence measurements from a 30 m tower in the stably stratifiedboundary layer over the Greenland ice sheet are analyzed. The observationsinclude profile and eddy-correlation measurements at various levels. Atfirst, the analysis of the turbulence data from the lowest level (2 m aboveground) shows that the linear form of the non-dimensional wind profile(m) is in good agreement with the observations for z/L <0.4, whereL represents the Obukhov length. A linear regression yieldsm=1+5.8z/L. The non-dimensional temperature profile (h) at the2m level shows no tendency to increase with increasing stability. The datafrom the upper levels of the tower are analyzed in terms of both localscaling and surface-layer scaling. The m and the h values show atendency to level off at large stability (z/>0.4) where represents the local Obukhov length. Hence, the linear form of the functions is no longer appropriate under such conditions. The bestcorrespondence to the data can be achieved when using the expression ofBeljaars and Holtslag for m and h. The vertical profiles of theturbulent fluxes, the wind velocity variances and temperature variance arealso determined. The momentum flux profile and the profiles of wind speedvariances are in general agreement with other observations if a welldeveloped low-level wind maximum occurs, and the height of this maximum isused as a height scale.     

Parameterizing turbulent diffusion through the joint probability density

The convective mass flux parameterization often used in meteorological modeling expresses the vertical flux of a transported scalar as proportional to the product of the difference in mean values of the scalar in updrafts and downdrafts and their characteristic velocity. The proportionality factor is a constant to be specified. We show that this proportionality factor also appears in the relaxed eddy accumulation technique of Businger and Oncley. That associates the surface-layer flux of a scalar with the product of the standard deviation of vertical velocity and the mean concentration difference between updrafts and downdrafts.We show that this constant (b) is determined uniquely by the joint probability density (jpd) of vertical velocity and the scalar. Using large-eddy simulation, we generate this jpd for a conservative scalar diffusing through a convective boundary layer. It has quite different forms in top-down and bottom-up diffusion geometries. The bottom-up jpd is fairly well represented by a jointly Gaussian form and implies b ~ 0.6, in good agreement with the surface-layer value reported by Businger and Oncley. The top-down jpd is strikingly non-Gaussian and gives b ~ 0.47. Updrafts carry the bulk of the scalar flux - 70% in the bottom-up case, 60% in the top-down case.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Remote sensing of the dynamic stability of the atmospheric boundary layer

Summary The relative strength of the stabilizing effect of buoyancy and the destabilizing effect of velocity shear in a stratified shear flow, such as a stable atmospheric boundary layer, is measured by the gradient Richardson number, Ri_g. The boundary layer static stability, as described by the buoyancy frequency, N, can be calculated from the virtual potential temperature gradient derived from RASS temperature profiles. The mean wind profiles from a sodar can be used to calculate the mean vertical velocity shear. In combination these profilers are potentially a powerful tool for the remotely sensing the dynamic stability of the boundary layer. However, experience shows that the combinations of two experimentally derived quantities, like N and shear, may give highly variable results. On the other hand, a simple sensitivity analysis shows that reasonable estimates of Ri_g are achievable over a range of conditions in the stable nocturnal boundary layer. To test this conclusion, high spatial and temporal resolution temperature and velocity soundings were obtained above 50m in the stable nocturnal boundary layer using a 920MHz continuous wave Radio Acoustic Sounding System (RASS) and 1.875kHz and 5.00kHz Doppler sodars. Examples of the evolution of Ri_g are presented from 24 hours of observations of the boundary layer in Canberra, on the tablelands in south- eastern Australia. Most of the boundary layer had Ri_g between 0.1 and 1. Thus, it was marginally dynamically stable, even with the gradient Richardson number calculated from finite differences over a vertical interval of 68m. A comparison of the results from the two sodars showed that the velocity shear increased significantly when the vertical differencing interval was decreased from 68m to 20m.     

Methodological Development of the Conditional Sampling Method. Part I: Sensitivity to Statistical and Technical Characteristics

The relaxed eddy accumulation (REA), method based on the conditional sampling concept, has received increasing attention over the past few years as it can be used to measure surface fluxes of a wide variety of trace gases for which fast response analysers are not available. In the REA method, a turbulent flux is simply expressed as the product of the standard deviation of vertical wind velocity, the difference between mean scalar concentration in the updrafts and downdrafts and an empirical coefficient, (about 0.63 as based on simulations with a Gaussian distribution, and 0.58 as derived from experimental data). A simulation technique is developed here to evaluate the performance of a ground-based REA system. This analysis uses generated series whose internal structure can be controlled to a large extent. They are stationary and their characteristics are similar to those of physical turbulence. In a first step the influence of some statistical characteristics of vertical velocity and scalar concentration series is investigated. The effect of the third- and fourth-order moments can explain to some degree the difference between calculated and measured values. The impact of a threshold on the vertical velocity is then considered, and its effect on is quantified. The influence of the time lag between and the effective scalar sampling, and the consequences of lowpass filtering of the signal are also investigated. The simulation technique presented in this study can be used to develop elaborate algorithms for near real-time conditional sampling, based on the statistical characteristics of the previous sample.     

Interaction between the atmospheric and oceanic boundary layers

The two-layer system of an atmosphere over water bodies is reduced to a single-layer problem. Values of the interfacial quantities, such as the friction velocity, the surface velocity, the angles, and , between the surface shear stress and the geostrophic wind velocity and the surface wind velocity, respectively, and the surface roughness, all of which depend upon external parameters, such as the geostrophic wind and stratifications, are obtained. The geostrophic drag coefficient C _d, the geostrophic wind coefficient C _f, and the angles , and , of the turbulent flow at the sea-air interface are functions of a dimensionless number, mfG/kg, with S ₁ and S ₂ as two free stratification parameters. The surface roughness is uniquely determined from the geostrophic wind rather than from the wind profile in the boundary layer.Formerly Visiting Research Associate, Applied Physics Branch, Earth Observations Division, NASA-Manned Spacecraft Center, Houston, Texas.     

Temperature distribution and current system in a convectively mixed lake

During spring and autumn, many lakes in temperate latitudes experience intensive convective mixing in the vertical, which leads to almost isothermal conditions with depth. Thus the regime of turbulence appears to be similar with that characteristic of convective boundary layers in the atmosphere. In the present paper a simple analytical approach, based on boundary-layer theory, is applied to convective conditions in lakes. The aims of the paper are firstly to analyze in detail the temperature distribution during these periods, and secondly to investigate the current system, created by the horizontal temperature gradient and wind action. For these purposes, simple analytical solutions for the current velocities are derived under the assumption of depth-constant temperatures. The density-induced current velocities are shown to be small, in the order of a few mm/sec. The analytical model of wind-driven currents is compared with field data. The solution is in good qualitative agreement with observed current velocities under the condition that the wind field is steady for a relatively long time and that residual effects from former wind events are negligible.The effect of the current system on an approximately depth-constant temperature distribution is then checked by using the obtained current velocity fields in the heat transfer equation and deriving an analytical solution for the corrected temperature field. These temperature corrections are shown to be small, which indicates that it is reasonable to describe the temperature distribution with vertical isotherms.Notation T temperature - t time - x, y, z cartesian coordinates - molecular viscosity - _h , _v horizontal and vertical turbulent viscosity - K _h ,K _v horizontal and vertical turbulent conductivity - Q heat flux through the water surface - D depth - u, v, w average current velocity components inx, y andz directions - f Coriolis parameter - p pressure - density - g gravity acceleration - a constant in the freshwater state equation - h _s deviation from the average water surface elevation - L _*,H _* length and depth scale - U _*,W _* horizontal and vertical velocity scale - T temperature difference scale - bottom slope - u _* friction velocity at the water surface - von Karman constant - L Monin-Obukhov length scale - buoyancy parameter - l turbulence length scale - C ₁,C ₂,C ₃ dimensionless constants in the expressions for the vertical turbulent viscosity - , dimensionless vertical coordinate and dimensionless local depth - angle between surface stress direction andx-axis - T _bx ,T _by bottom stress components - c bottom drag coefficient     

Improving the diagnostic relation for the nocturnal boundary-layer height

A one-dimensional model of the nocturnal boundary layer (NBL) has been used to investigate the time variation of the NBL height for stationary and horizontally homogeneous synoptic conditions. The time variation of the well known quantity = hflu _* has been shown to be related to the wind variation at the top of the NBL. For the simple simulated conditions, this variation depends only on the roughness length and the Coriolis parameter. The value of averaged over the whole night is a function of the friction velocity. An expression is proposed for which is compared with observations. Under stationary external conditions, the new relation improves the determination of the NBL height if compared with the classical relation using a constant value of .     

A comparative study of the performance of various vertical discretization schemes

Summary Most finite-difference numerical weather prediction models employ vertical discretizations that are staggered, and are low-order (usually second-order) approximations for the important terms such as the derivation of the geopotential from the hydrostatic equation, and the calculation of the vertically integrated divergence. In a sigma-coordinate model the latter is used for computing both the surface pressure change and the vertical velocity. All of the above-mentioned variables can diminish the accuracy of the forecast if they are not calculated accurately, and can have an impact on related quantities such as precipitation.In this study various discretization schemes in the vertical are compared both in theory and in practice. Four different vertical grids are tested: one unstaggered and three staggered (including the widely-used Lorenz grid). The comparison is carried out by assessing the accuracy of the grids using vertical numerics that range from second-order up to sixth-order.The theoretical part of the study examines how faithfully each vertical grid reproduces the vertical modes of the governing equations linearized with a basic state atmosphere. The performance of the grids is evaluated for 2nd, 4th and 6th-order numerical schemes based on Lagrange polynomials, and for a 6th-ordercompact scheme.Our interpretation of the results of the theoretical study is as follows. The most important result is that the order of accuracy employed in the numerics seems to be more significant than the choice of vertical grid. There are differences between the grids at second-order, but these differences effectively vanish as the order of accuracy increases. The sixth-order schemes all produce very accurate results with the grids performing equally well, and with the compact scheme significantly outperforming the Lagrange scheme. A second major result is that for the number of levels typically used in current operational forecast models, second-order schemes (which are used almost universally) all appear to be relatively poor, for other than the lowest modes.The theoretical claims were confirmed in practice using a large number (100) of forecasts with the Australian Bureau of Meteorology Research Centre's operational model. By comparing test model forecasts using the four grids and the different orders of numerics with very high resolution control model forecasts, the results of the theoretical study seem to be corroborated.With 8 Figures     

SODAR-based estimation of TKE and momentum flux profiles in the atmospheric boundary layer: Test of a parameterization model

Summary A simple parameterization for the estimation of turbulent kinetic energy (TKE) and momentum flux profiles under near-neutral stratification based on sodar measurements of the vertical velocity variance has been tested using data from the LINEX-2000 experiment. Measurements included operation of a phased-array Doppler sodar DSDPA.90 and of a sonic anemometer USA-1 mounted at a meteorological tower at a height of 90m. Good agreement has been found between the TKE and momentum flux values derived from the sonic and sodar data (with correlation coefficients r>0.90 and a slope of the regression lines of about 1.01.1) suggesting the possible use of sodar measurements of _w ² to derive turbulence parameter profiles above the tower range.     

Modification Of The Standard ∈-Equation For The Stable Abl Through Enforced Consistency With Monin–Obukhov Similarity Theory

A condition is derived for consistency of the standard-equation with Monin–Obukhov (MO) similarity theory of thestably-stratified surface layer. The condition is derivedby extending the procedure used to derive the analogous condition forneutral theory to stable stratification. It is shown that consistencywith MO theory requires a function of flux Richardson number, Ri_f, to be absorbed into either of two closure parameters, c ₁ or c ₂.Inconsistency, on the other hand, results if constant values of these are maintained for all Ri_f, as is done in standardapplication of the equation, and the large overpredictions ofturbulence found in such application to the one-dimensionalstable atmospheric boundary layer (1D-SBL) are traced to thisinconsistency. Guided by this, we formulate a MO-consistent-equation by absorbing the aforementioned function intoc ₁, and combine this with a Level-2.5 second-orderclosure model for vertical eddy viscosity and diffusivities.Numerical predictions of the 1D-SBL by the modified model converge to a quasi-steady state, rectifying the predictive failure of the standard -equation for the case.Quasi-steady predictions of non-dimensional variables agree stronglywith Nieuwstadt's theory. Qualitative accuracy of predictionsis inferred from comparisons to field data, large-eddy simulationresults and Rossby-number similarity relationships.     

The structure of the stably stratified internal boundary layer in offshore flow over the sea

Observations obtained mainly from a research aircraft are presented of the mean and turbulent structure of the stably stratified internal boundary layer (IBL) over the sea formed by warm air advection from land to sea. The potential temperature and humidity fields reveal the vertical extent of the IBL, for fetches out to several hundred of kilometres, geostrophic winds of 20–25 m s^–1, and potential temperature differences between undisturbed continental air and the sea surface of 7 to 17 K. The dependence of IBL depth on these external parameters is discussed in the context of the numerical results of Garratt (1987), and some discrepancies are noted.Wind observations show the development of a low-level wind maximum (wind component normal to the coast) and rotation of the wind to smaller cross-isobar flow angles. Potential temperature () profiles within the IBL reveal quite a different structure to that found in the nocturnal boundary layer (NBL) over land. Over the sea, profiles have large positive curvature with vertical gradients increasing monotonically with height; this reflects the dominance of turbulent cooling within the layer. The behaviour is consistent with known behaviour in the NBL over land where curvature becomes negative (vertical gradients of decreasing with height) as radiative cooling becomes dominant.Turbulent properties are discussed in terms of non-dimensional quantities, normalised by the surface friction velocity, as functions of normalised height using the IBL depth. Vertical profiles of these and the normalised wavelength of the spectral maximum agree well with known results for the stable boundary layer over land (Caughey et al., 1979).     

A verification of some methods to determine the fluxes of momentum,sensible heat,and water vapour using standard deviation and structure parameter of scalar meteorological quantities

A set of micro-meteorological data collected over a horizontal, uniform terrain (the plain of La Crau, France) in June 1987 is analysed. Conditions were predominantly sunny and arid, while due to the Mistral the wind speed could exceed 10 m/s. Verification of several methods to evaluate surface fluxes of heat, momentum and water vapour from the standard deviation of temperature, wind and specific humidity is presented. Also, a similar approach using the structure parameter of temperature is considered. These methods are all based on Monin-Obukhov (M-O) similarity theory. It is found that the standard deviation of temperature, vertical and horizontal wind speed as well as the structure parameter for temperature behave according to M-O similarity. It is shown that the sensible heat flux and friction velocity can be determined from a fast response thermometer and a cup anemometer. Also, it appears that the analytic solution of the set of governing equations as derived by the first author yields good results. M-O theory does not appear to work for the standard deviation of specific humidity. This may be due to the relative importance of large eddies.     

Observational aspects on a dominant spectral frequency in an airflow over a weeping-lovegrass stand

Mean and fluctuating wind velocities were measured above a flexible stand (weeping-lovegrass). A waving phenomenon Honami appeared over the stand during the observation period. Some spectral parameters were derived from the vertical wind fluctuations. A dependency of frequency on mean horizontal wind velocity was found. The result, n _m = 0.66u, was obtained under the range of wind speeds from 0.9 m s^-1 to 3.1 m s^-1 just above the canopy.     

Markov chain simulations of vertical dispersion in the neutral surface layer for surface and elevated releases

Vertical dispersion in the neutral surface layer is investigated using a Markov Chain simulation procedure. The conceptual basis of the procedure is discussed and computation procedures outlined. Wind and turbulence parameterizations appropriate to the neutral surface layer are considered with emphasis on the Lagrangian time scale. Computations for a surface release are compared with field data. Good agreement is found for the variation of surface concentration and cloud height to distances 500 m downwind of the source. The functional form of the vertical concentration profile is examined and an exponential with exponent 1.6 is found to give the best fit with simulations.For elevated releases, it is demonstrated that an initial dip of the mass mean height from the simulation can be normalized for various release heights using a non-dimensionalized downwind coordinate incorporating advective wind speed and wind shear. The vertical distribution standard deviation ( _z ), as employed in Gaussian models, shows a fair degree of independence with source height but close examination reveals an optimum source height for maximum _z at a given downwind distance,x. This source height increases with downwind distance. Also the simulations indicate that vertical wind shear is more important than vertical variation of Lagrangian time scale close to the source, with a reverse effect farther downwind.     

Effects of heat and water vapor transport on eddy covariance measurement of CO2 fluxes

Flux densities of carbon dioxide were measured over an arid, vegetation-free surface by eddy covariance techniques and by a heat budget-profile method, in which CO₂ concentration gradients were specified in terms of mixing ratios. This method showed negligible fluxes of CO₂, consistent with the bareness of the experimental site, whereas the eddy covariance measurements indicated large downward fluxes of CO₂. These apparently conflicting observations are in quantitative agreement with the results of a recent theory which predicts that whenever there are vertical fluxes of sensible or latent heat, a mean vertical velocity is developed. This velocity causes a mean vertical convective mass flux (= _cw for CO₂, in standard notation). The eddy covariance technique neglects this mean convective flux and measures only the turbulent flux _c w. Thus, when the net flux of CO₂ is zero, the eddy covariance method indicates an apparent flux which is equal and opposite to the mean convective flux, i.e., _c w = – _c w. Corrections for the mean convective flux are particularly significant for CO₂ because _cw and _c w are often of similar magnitude. The correct measurement of the net CO₂ flux by eddy covariance techniques requires that the fluxes of sensible and latent heat be measured as well.