A review of flux-profile relationships

Flux-profile relationships in the constant flux layer are reviewed. The preferred relationships are found to be those of Dyer and Hicks (1970), namely, _H = _W =(1–16(z/L))^–1/2, _M =(1–16(z/L))^–1/4 for the unstable region, and _H = _W = _M = 1+5(z/L) for the stable region.The carefully determined results of Businger et al. (1971) remain a difficulty which calls for considerable clarification.     

Analytical solutions for the Ekman layer

The PBL equation that governs the transition from the constant-stress surface layer to the geostrophic wind in a neutrally stratified atmosphere for which the eddy viscosityK(z) is assumed to vary smoothly from the surface-layer value U _*z (0.4,U _*=friction velocity,z=elevation) to the geostrophic asymptoteK _GU _*d forzd is solved through an expansion in fd/U _*1 (f=Coriolis parameter). The resulting solution is separated into Ekman's constant-K solution an inner component that reduces to the classical logarithmic form forzd and isO() relative to the Ekman component forzd. The approximationKU _*d is supported by the solution of Nee and Kovasznay's phenomenological transport equation forK(z), which yieldsK–U _*d exp(–z/d), where is an empirical constant for which observation implies, 1. The parametersA andB in Kazanskii and Monin's similarity relation forG/U _* (G=geostrophic velocity) are determined as functions of . The predicted values ofG/U _* and the turning angle are in agreement with the observed values for the Leipzig wind profile. The predicted value ofB based on the assumption of asymptotically constantK is 4.5, while that based on the Nee-Kovasznay model is 5.1; these compare with the observed value of 4.7 for the Leipzig profile. A thermal wind correction, an asymptotic solution for arbitraryK(z) and 1, and an exact (unrestricted ) solution forK(z)=U _*d[1–exp(–z/d)] are developed in appendices.     

Numerical modeling of the nocturnal urban boundary layer

A numerical case study with a second-order turbulence closure model is proposed to study the role of urban canopy layer (UCL) for the formation of the nocturnal urban boundary layer (UBL). The turbulent diffusion coefficient was determined from an algebraic stress model. The concept of urban building surface area density is proposed to represent the UCL. Calculated results were also compared with field observation data. The height of the elevated inversion above an urban center was simulated and found to be approximately twice the average building height. The turbulent kinetic energy k, energy dissipation rate , and turbulence intensities u ² and w ² increase rapidly at the upwind edge of the urban area. The Reynolds stress uw displayed a nearly uniform profile inside the UBL, and the vertical sensible heat flux w had a negative value at the inversion base height. This indicates that the downward transport of sensible heat from the inversion base may play an important role in the formation of the nocturnal UBL.     

Standard deviation of vertical two-point longitudinal velocity differences in the atmospheric boundary layer

The standard deviation of vertical two-point longitudinal velocity fluctuation differences is analyzed experimentally with eleven sets of turbulence measurements obtained at the NASA 150-m ground-winds tower site at Cape Kennedy, Florida. It is concluded that /u _*0 is proportional to (fz/u _*0)^0.22, where the coefficient of proportionality is a function of fz/u _*0 and u _*0/fL ₀. The quantities f and L₀ denote the Coriolis parameter and the surface Monin-Obukhov stability length, respectively; u _*0 is the surface friction velocity; z is the vertical distance between the two points over which the velocity difference is calculated; and zz is the mean height of the mid-point of the interval z above natural grade. The results of the analysis are valid for 20<-u _*0/fL ₀<2000.     

Radiometrically determined skin temperature and scalar roughness to estimate surface heat flux part II: Performance of parameterized scalar roughness for the determination of sensible heat

Several formulations and proposals to determine the value of the radiometric scalar roughness for sensible heatz _0h,r are tested with respect to their performance in the estimation of the sensible heat flux by means of the profile equations derived from Monin-Obukhov similarity theory. The equations are applied to the data set of spatially averaged surface skin temperature and profiles of wind speed and temperature observed in a pasture field during a growing season. The use of a physical model developed for a dense canopy to estimate scalar roughness for sensible heatz _0h,r produced sensible heat fluxH with a correlation coefficientr=0.884, the ratio of means being H _s /H=1.19 in a comparison with reference values ofH _s . In comparison, a proposal for a fixed value ofz _0h yieldedr=0.887, H _s /H=0.879. In both cases, the validity ofz _0h =z _0h,r was assumed. All expressions derived to estimatez _0h,r from a multiple linear regression with such predictors as leaf area index, solar radiation and the ratio of solar radiation to extraterrestrial radiation, were found to produce a better result, withr better than 0.90 and H _s /H around 1.0. However, when the constantsc andf of a linear regression equationHs=cH+f are used to evaluate the equations, a marked difference in performance of each formulation appeared. In general, equations with smaller numbers of predictors tend to produce a biased result, i.e., an overestimation ofH at largeH _s . These values ofH are used in conjunction with the energy balance equation to derive values of the latent heat fluxLE, which are shown to be in good agreement with the reference valuesLE _s , withr greater than 0.97.     

A study of the seasonal migration of ITCZ and the quasi-periodic oscillations in a simple monsoon system using an energy balance model

Summary A zonally averaged global energy balance model with feedback mechanisms was constructed to simulate (i) the poleward limits of ITCZ over the continent and over the ocean and (ii) a simple monsoon system as a result of differential heating between the continent and the ocean. Three numerical experiments were performed with lower boundary as (1) global continent, (2) global ocean and (3) continent-ocean, with freezing latitudes near the poles. Over the continent, midlatitude deserts were found and the ITCZ migrates 25° north and south with seasons. Over a global swamp ocean results do not show migration of ITCZ with time but once the ocean currents are introduced the ITCZ migrates 5° north and south with seasons. It was found that the seasonal migration of ITCZ strongly depends on the meridional distribution of the surface temperature. It was also found that continent influences the location of the oceanic ITCZ. In the tropics northward progression of quasi-periodic oscillations called events are found during the pre- and post-monsoon periods with a period of 8 to 15 days. This result is consistent with the observed quasi-periodic oscillations in the tropical region. Northward propagation of the surface temperature perturbation appears to cause changes in the sensible heat flux which in turn causes perturbations in vertical velocity and latent heat flux fields.List of Symbols vertical average - ₀ zonal average - vertical mean of the zonal average - _0s zonal average at the surface - _0a zonal average at 500 mb level - latitude We now define the various symbols used in the model rate of atmospheric heating due to convective cloud formation (K/sec) - dp/dt (N/m²/sec) - density - potential temperature (K) - rate of rotation of the earth (rad/sec) - empirical constant - humidity mixing ratio - ^* saturated humidity mixing ratio - opacity of the atmosphere - ₁,₂ factors for downward and upward effective black body long wave radiation from the atmosphere - Stefan-Boltzmann constant - emissivity of the surface - _D subsurface temperature (K) - a specific volume - _0xs ,_0ys eastward and northward components of surface frictional stress - ^* vertical velocity at the top of the boundary layer (N/m²/sec) - P Thickness of the boundary layer (mb) - nondimensional function of pressure - P pressure - P _a pressure of the model atmosphere (N/m²) - P _s pressure at the surface (N/m²) - t time (sec) - U eastward wind speed (m/sec) - V northward wind speed (m/sec) - surface water availability - T absolute temperature (K) - heat addition due to water phase changes - g acceleration due to gravity (m²/sec) - a radius of the earth (m) - R gas constant for dry air (J/Kg/K) - C _p specific heat of air at constant pressure (J/Kg/K) - k R/C _p - L latent heat of condensation (J/Kg) - f coriolis parameter (rad/sec) - H _s H _0s ⁽¹⁾ +H _0s ⁽²⁾ +H _0s ⁽³⁾ +H _0s ⁽⁴⁾ +H _0s ⁽⁵⁾ (J/m²/Sec)=sum of the rates of vertical heat fluxes per unit surface area, directed toward the surface - H _a H _0a ⁽¹⁾ +H _0a ⁽²⁾ +H _0a ⁽³⁾ +H _0a ⁽⁴⁾ (J/m²/Sec)=sum of the rates of heat additions to the atmospheric column per unit horizontal area by all processes - H _0s ⁽¹⁾ ,H _0a ⁽¹⁾ heat flux due to short wave radiation - H _0s ⁽²⁾ ,H _0a ⁽²⁾ heat flux due to long wave radiation - H _0s ⁽³⁾ ,H _0a ⁽³⁾ heat flux due to small scale convection - H _0s ⁽⁴⁾ heat flux due to evaporation - H _0a ⁽⁴⁾ heat flux due to condensation - H _0s ⁽⁵⁾ heat flux due to subsurface conduction and convection - e ^* saturation vapor pressure - R solar constant (W/m²) - r _a albedo of the atmosphere - r _s albedo of the surface - b ₂ empirical constant (J/m²/sec) - c ₂ empirical constant (J/m²/sec) - e ₂ nondimensional empirical constant - f ₂ empirical constant (J/m₂/sec) - factor proportional to the conductive capacity of the surface medium - a _s constant used in Sellers model - b _s positive constant of proportionality used in the Sellers model (kg m²/J/sec²) - K _HT coefficient for eddy diffusivity of heat (m²/sec) - K _HE exchange coefficient for water vapor (m²/sec) - h depth of the water column (m) - z height (m) - V _0ws meridional component of surface current (m/sec) - n cloud amount - G _0,n long wave radiation form the atmosphere for cloud amount n (W/m²) - B ₀ long wave radiation from the surface (W/m²) - S _0,n short wave radiation from the atmosphere for cloud amount n (W/m²) - A _n albedo factor for a cloud amount n - R _f1 large scale rainfall (mm/day) - R _f2 small scale rainfall (mm/day) With 22 Figures     

Near-field dispersion from instantaneous sources in the surface layer

This paper considers the near-field dispersion of an ensemble of tracer particles released instantaneously from an elevated source into an adiabatic surface layer. By modelling the Lagrangian vertical velocity as a Markov process which obeys the Langevin equation, we show analytically that the mean vertical drift velocity w(t) is w()=bu _*(1–e ^–(1+)), where is time since release (nondimensionalized with the Lagrangian time scale at the source), b Batchelor's constant, and u _*, the friction velocity. Hence, the mean height and mean depth of the ensemble are calculated. Although the derivation is formally valid only when 1, the predictions for w, mean height and mean depth are consistent in the downstream limit ( 1) with surface-layer Lagrangian similarity theory and with the diffusion equation. By comparing the analytical predictions with numerical, randomflight solutions of the Langevin equation, the analytical predictions are shown to be good approximations at all times, both near-field and far-field.     

Free convection frequency spectra of atmospheric turbulence over the sea

Frequency spectra of atmospheric turbulenceS (f) in the inertial subrange are considered in the free convection regime over the sea surface in a case of motionless instrument measurements (Eulerian frequency spectra). The frequency spectra formulaef _* S (f)/ ² =c (f _*/f)^5/3 for wind velocity (=1–3), temperature (=t) and humidity (=e) fluctuations are derived on the basis of similarity theory and the –5/3 law. These relations also can be derived from a consideration of convective large-scale advection of small eddies. The frequency scalef _* = (N ^{1 2}/)^1/2 (H/z ²)^1/3 is the lower bound of the inertial subrange and it is of order 10^–2 Hz.The spectra formulae are compared with direct measurements of atmospheric turbulence from the fixed research tower in the coastal zone of the Black Sea in calm weather. It is shown that these formulae are realized at least over two to three decades of the frequency range (approximately from 10^–2 to 10 Hz) and values of the numerical coefficients are found. The derived formulae can be used for calculations of sensible and latent heat fluxes by measuring the high-frequency range of spectra at a fixed point at low wind speeds when the conventional inertial dissipation method is not applicable.     

Equatorial disturbances,monsoons, and 1982–1983 ENSO

Summary Interannual modes are described in terms of three-month running mean anomaly winds (u,v), outgoing longwave radiation (OLR), and sea surface temperature (T _* ). Normal atmospheric monsoon circulations are defined by long-term average winds (u _n,v _n) computed every month from January to December. Daily winds are grouped into three frequency bands, i.e., 30–60 day filtered winds (u _L,v _L); 7–20 day filtered winds (u _M,v _M); and 2–6 day filtered winds (u _S,v _S). Three-month running mean anomaly kinetic energy (signified asK _L , K _M , andK _S , respectively) is then introduced as a measure of interannual variation of equatorial disturbance activity. Interestingly, all of theseK _L , K _M , andK _S perturbations propagate slowly eastward with same phase speed (0.3 ms^–1) as ENSO modes. Associated with this eastward propagation is a positive (negative) correlation between interannual disturbance activity (K _L , K _M , K _S ) and interannualu (OLR) modes. Namely, (K _L , K _M , K _S ) becomes more pronounced than usual nearly simultaneously with the arrival of westerlyu and negativeOLR (above normal convection) perturbutions. In these disturbed areas with (K _L , K _M , K _S >0), upper ocean mixing tends to increase, resulting in decreased sea surface temperature, i.e.T _* 0. Thus, groups (not individual) of equatorial disturbances appear to play an important role in determiningT _* variations on interannual time scales. HighestT _* occurs about 3 months prior to the lowestOLR (convection) due primarily to radiational effects. This favors the eastward propagation of ENSO modes. The interannualT _* variations are also controlled by the prevailing monsoonal zonal windsu _n, as well as the zonal advection of sea surface temperature on interannual time scales. Over the central Pacific, all of the above mentioned physical processes contribute to the intensification of eastward propagating ENSO modes. Over the Indian Ocean, on the other hand, some of the physical processes become insignificant, or even compensated for by other processes. This results in less pronounced ENSO modes over the Indian Ocean.With 10 FiguresContribution No. 89-6, Department of Meteorology, University of Hawaii, Honolulu, Hawaii.     

Effect of finite sampling on atmospheric spectra

The effect of a finite averaging time on variances is well known, but its effect on power spectra is less clearly understood. We present numerical solutions for the spectral distortion arising from sampling over a finite time interval T and show that the commonly used filter function (1 – sinc²f T), valid for variances, is a reasonable approximation for power spectra only when T 10 _m , where f is the cyclic frequency, and _m is the dominant time scale of the process. Our results exhibit an increasingly steeper low-frequency roll-off as T decreases relative to _m , indicating that the measured spectrum is subject to a greater suppression of the lower frequencies (f > 1/T) than predicted by (1 – sinc²f T). This suppression is, in a sense, compensated by an overestimation of spectral estimates in the frequency range f 1/T.     

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.     

Nearly steady convection and the boundary-layer budgets of water vapor and sensible heat

The budgets of water vapor and sensible heat in the convective atmospheric boundary (mixed) layer are analyzed by means of a simple slab approach adapted to steady large-scale advective conditions with radiation and cloud activity. The entrainment flux for sensible heat is assumed to be a linear function of the surface flux. The flux of water vapor at the top of the mixed layer is parameterized by extending the first-order Betts-Deardorff approach, i.e., by adopting linear changes for both the specific humidity and the flux across the mixed layer and across the inversion layer of finite thickness. In this way the dissimilarity of sensible heat and water vapor transport in the mixed layer can be taken into account. The experimental data were obtained from the Air Mass Transformation Experiment (AMTEX). The entrainment constant for sensible heat at the top of the mixed layer was found to have values similar to those observed in other weakly convective situations, i.e., around 0.4 to 0.6. This appears to indicate that the effect of mechanical turbulence was not negligible; however, the inclusion of this effect in the formulation did not improve the correlation. In contrast to the first-order approach, the zero-order approach, i. e., the jump equation commonly used for the flux of a scalar at the inversion, (ovwc ) _h = w_e c (where w _e is the entrainment velocity and c the concentration jump across the inversion), was found to be invalid and incapable of describing the data.     

Boundary-layer parameterizations for applied dispersion modeling over urban areas

Many applied dispersion models require the knowledge of boundary-layer parameters such as sensible heat flux,Q _H , friction velocity,u _*, and turbulent energy components, _w and _v . Formulas are suggested for calculating these parameters over a wide variety of types of ground surfaces, based on simple observations of wind speed near the ground and fractional cloud cover, and specification of constants such as roughness length, albedo, and soil moisture availability. Observations ofu _*,Q _H , _w , and _v during field experiments in St. Louis and Indianapolis are used to test the formulas for urban sites. Relative errors of about ±20% in the predictions are seen to occur whenu _*,Q _H , _w , and _v are large. However, when these quantities are small (e.g.,u _* < 0.2 m/s), the errors in the predictions are as large as the mean value of the quantity itself.In addition, it is concluded from studies of available field data and theories that the magnitude of _w is not well-known at elevations above about 100m during the late afternoon and night. Some simple parameterizations for _w . are suggested that are consistent with the observed steady decrease in ground-level concentration in the afternoon and the sudden increase in concentration that can occur a few hours after sunset due to wind shears associated with a low-level jet, for continuous plumes emitted from moderate to tall stacks.     

Reactions of alkoxy radicals under atmospheric conditions: The relative importance of decomposition versus reaction with O2

The reactions of alkoxy radicals determine to a large extent the products formed during the atmospheric degradations of emitted organic compounds. Experimental data concerning the decompositions, 1,5-H shift isomerizations and reactions with O₂ of several classes of alkoxy radicals are inconsistent with literature estimations of their absolute or relative rate constants. An alternative, although empirical, method for assessing the relative importance under atmospheric conditions of the reactions of alkoxy radicals with O₂ versus decomposition was derived. This estimation method utilizes the differences in the heats of reaction, (H)=(H_{decomposition}–H_{O 2 reaction}), between these two reactions pathways. For (H)[22–0.5(H_{O 2 reaction})], alkoxy radical decomposition dominates over the reaction with O₂ at room temperature and atmospheric pressure of air, while for (H)[25-0.5(H_{O 2 reaction})], the O₂ reaction dominates over decomposition (where the units of H are in kcal mol^–1). The utility and shortcomings of this approach are discussed. It is concluded that further studies concerning the reactions of alkoxy radicals are needed.     

Flux determination over a smooth surface under strongly unstable conditions

Careful micrometeorological measurements on an empty parking lot allowed determination of the surface fluxes of sensible heatH and of momentum by applying profile equations derived from Monin-Obukhov similarity theory with two sets of the stability correction function for momentum _m and sensible heat _h . These fluxes were compared with reference values ofH independently determined by means of an eddy correlation technique. In general, better agreement was found betweenH values derived from profiles with the stability functions of Brutsaert (1992) and referenceH values, than when the Businger-Dyer functions were used to deriveH. The disagreement in the latter comparison was especially serious under strongly unstable conditions, with the value ofy=–z/L (wherez is the height andL is the Obukhov length) larger than 10. A closer look at the procedure for calculatingH from the profiles revealed that the large differences between theH values derived with these two different versions of the stability correction functions were caused by the small differences of the _h values, and not by the larger differences of the _m values. This result stems from the strong sensitivity of the resultingH values on the choice of _h .     

Turbulence Reynolds number and the turbulent kinetic energy balance in the atmospheric surface layer

The relation between the turbulence Reynolds numberR and a Reynolds numberz* based on the friction velocity and height from the ground is established using direct measurements of the r.m.s. longitudinal velocity and turbulent energy dissipation in the atmospheric surface layer. Measurements of the relative magnitude of components of the turbulent kinetic energy budget in the stability range 0 >z/L 0.4 indicate that local balance between production and dissipation is maintained. Approximate expressions, in terms of readily measured micrometeorological quantities, are proposed for the Taylor microscale and the Kolmogorov length scale .     

The Askervein Hill Project: Mean wind variations at fixed heights above ground

This is one of a series of papers on the Askervein Hill Project. It presents results on the variations in mean wind speed at fixed heights (z) above the ground from linear arrays of anemometer posts and towers. Most of the data are for z = 10 m but some are for z = 3 m. Selected and directionally grouped data from the 55 Mean Flow runs are presented together with mean flow data from Askervein '83 Turbulence runs. Comparisons are made between the data and guideline estimates of fractional speed-up ratio at hilltop locations and between the data and MS3DJH/3 model predictions along the tower lines. There is good agreement in most cases.     

Diurnal variation of horizontal wind direction fluctuations in complex terrain at Geysers,Cal.

Horizontal diffusion in the surface layer is dependent on the standard deviation of wind direction fluctuations . Diurnal variation of this parameter in complex terrain was studied for the July 1979 Geysers, Cal., experiment using data from a network of 11 short meteorological towers in the 25 km² Anderson Creek watershed Valley side slopes are roughly 20 ° and maximum terrain difference is about 1 km.Values of for wind directions sampled for one hour at a height of 10 m are about 35 ° during the daytime. They slowly decrease to about 20 ° by 8 to 10 p.m. as stability increases but wind speeds are still relatively high. After 10 p.m. the drainage flow sets in at most stations, with speeds of 1 to 2 m s^-1, and average increases to about 30° during the period 11 p.m. to 6 a.m. In general, highest values of at night are associated with lowest values of wind speed and greatest static stability. This enhancement of by the terrain suggests that horizontal diffusion at night always conforms to that expected during nearly neutral stabilities. That is, Pasquill class D diffusion applies to the horizontal component all night in complex terrain.     

Flux–Profile Relationships Over a Fetch Limited Beech Forest

The influence of an internal boundary layer and a roughness sublayer on flux–profile relationships for momentum and sensible heat have been investigated for a closed beech forest canopy with limited fetch conditions. The influence was quantified by derivation of local scaling functions for sensible heat flux and momentum (_h and _m) and analysed as a function of atmospheric stability and fetch. For heat, the influences of the roughness sublayer and the internal boundary layer were in agreement with previous studies. For momentum, the strong vertical gradient of the flow just above the canopy top for some wind sectors led to an increase in _m, a feature that has not previously been observed. For a fetch of 500 m over the beech forest during neutral atmospheric conditions, there is no height range at the site where profiles can be expected to be logarithmic with respect to the local surface. The different influence of the roughness sublayer on _h and _m is reflected in the aerodynamic resistance for the site. The aerodynamic resistance for sensible heat is considerably smaller than the corresponding value for momentum.     

Similarity of atmospheric Reynolds shear stress and heat flux fluctuations over a rough surface

Horizontal u and vertical w velocity fluctuations have been measured together with temperature fluctuations in the atmospheric surface layer, at a small height above a wheat crop canopy. Marginal probability density functions are presented for both individual fluctuations u, w, and for the instantaneous Reynolds stress uw, and heat fluxes w and u. Probability density functions of the velocity fluctuations deviate less significantly from the Gaussian form than the probability density of temperature. There appears to be closer similarity between statistics of the instantaneous heat fluxes than between the momentum flux and either of the heat fluxes investigated. The mean momentum flux receives equal contributions from the events referred to as ejections and sweeps in laboratory boundary layers. Sweeps provide the largest contribution to the heat fluxes.