Eddy energy dissipation rate and puff diffusion during calms

The paper considers a puff diffusion in its inertial stage when particle separation obeys the laws of the inertial subrange and depends only on eddy energy dissipation rate . The can be determined in the surface layer by the turbulent kinetic energy equation. Similarity equations connect with diffusion measure .A simple analytical model has been deduced to estimate pollutants diffusion during calms.     

AN E – ε – ℓ TURBULENCE CLOSURE SCHEME FOR PLANETARY BOUNDARY-LAYER MODELS: THE NEUTRALLY STRATIFIED CASE

The standard E – model generates aplanetary boundary layerthat appears to be much too deep. The cause of theproblem is traced to the equation for the dissipationrate () of turbulent kinetic energy (E), specifically theparameterization of dissipation production anddestruction. In the context of atmosphericboundary-layer modelling, we argue that a part of thedissipation production should be modelled as the inputto the spectral cascade from the energy-containingpart of the spectrum, with a characteristic length , while the equilibrium imbalancebetween local production and destruction ofdissipation is modelled as proportional toE2/E, as in the standard model. Wepropose an E – – turbulence closurescheme, in which both the mixing length, m, and are prescribed. The importance ofthe equation is diminished, though itstill determines the dissipation rate in the Eequation.     

Time series analyses of concentration and wind fluctuations

Analyses of concentration fluctuation (C) spectra from boundary-layer smoke plume experiments at six separate locations show that the spectra from these experiments generally exhibit an inertial subrange at high frequencies with a slope of -5/3 and indicate peak energy at a time period of about 50 to 100 s. These periods of peak energy are a factor of two to five less than those for the peak of the wind speed fluctuation (u or v) spectra. A general spectral formula fits normalized spectra from the U.S. and Australia, where the frequency, n, is made dimensionless by multiplying by the plume dispersion parameter, _y , and dividing by the wind speed, u. Peak energy occurs at a dimensionless frequency of n _y/u equal to about 0.15. The Kolmogorov constant in the inertial subrange is estimated from a set of averaged spectra. Cross-spectra indicate little relation between concentration and wind fluctuations. However, most of the correlation that exists is due to periods larger than about 10 or 20 s.     

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.     

The response of a propeller anemometer to turbulent flow with the mean wind vector perpendicular to the axis of rotation

The system transfer function ¦H(v)¦² at frequencyv (units of Hz) for a vertical velocity propeller anemometer in a statistically stationary and horizontally homogeneous turbulent flow is determined from: (1) experimental estimates of propeller velocity spectra; and (2) estimates of Eulerian vertical velocity spectra based on the hypothesis that degradation of the input vertical velocity Fourier components occurs in the inertial subrange. The experimental estimates of ¦H(v)¦² were adequately summarized with the mathematical expression for the system transfer function of a first-order system with parameterT which has units of time and is analogous to the time constant of a horizontal velocity propeller anemometer. Dimensional analysis techniques and the Monin-Obukhov similarity hypothesis were used to construct a model for the system parameterT which yielded the result that _w /D ₁ ( _w /)^1/3, where _w , andD ₁ denote the standard deviation of the input vertical velocity fluctuations, the horizontal mean wind speed, and the diameter of the propeller, respectively. The system parameterT is interpreted in terms of the time required for the propeller velocity statistics to become asymptotically independent of time upon being released from rest in a statistically stationary turbulent flow.Currently on leave of absence from the Indian Institute of Technology, New Delhi, India.     

Scavenging Efficiency of ‘Aerosol Carbon’ and Sulfate in Supercooled Clouds at Mt. Sonnblick (3106 m a.s.l., Austria)

Cloud water and interstitial aerosol samples collected at Mt. Sonnblick (SBO) were analyzed for sulfate and aerosol carbon to calculate in-cloud scavenging efficiencies. Scavenging efficiencies for sulfate (_SO) ranged from 0.52 to 0.99 with an average of 0.80. Aerosol carbon was scavenged less efficiently with an average value (_AC) of 0.45 and minimum and maximum values of 0.14 and 0.81, respectively. Both _SO and _AC showed a marked, but slightly different, dependence on the liquid water content (LWC) of the cloud. At low LWC, _SO increased with rising LWC until it reached a relatively constant value of 0.83 above an LWC of 0.3 g/m³. In the case of aerosol carbon, we obtained a more gradual increase of _AC up to an LWC of 0.5 g/m³. At higher LWCs, _{_} remained relatively constant at 0.60. As the differences between _SO and _A varied across the LWC range observed at SBO, we assume that part of the aerosol carbon was incorporated into the cloud droplets independently from sulfate. This hypothesis is supported by size classified aerosol measurements. The differences in the size distributions of sulfate and total carbon point to a partially external mixture. Thus, the different chemical nature and the differences in the size and mixing state of the aerosol particles are the most likely candidates for the differences in the scavenging behavior.     

Local Similarity Relationships In The Urban Boundary Layer

To investigate turbulent structures in an urban boundary layer (UBL) with many tallbuildings, a number of non-dimensional variable groups based on turbulent observationsfrom a 325-m meteorological tower in the urban area of Beijing, China, are analyzedin the framework of local similarity. The extension of surface-layer similarity to localsimilarity in the stable and unstable boundary layer is also discussed. According to localsimilarity, dimensionless quantities of variables: e.g., velocity and temperature standarddeviations _i/u_*l (i=u,v,w) and_T/T_*l,correlation coefficients of uw and wT covariance, gradients of wind and temperature_m and _h, and dissipation rates of turbulent kinetic energy (TKE) andtemperature variance and _N can be represented as a functiononly of a local stability parameter z/, where is the local Obukhovlength and z is the height above ground. The average dissipation rates of TKE andtemperature variance are computed by using the u spectrum, and the uw and wTcospectra in the inertial subrange. The functions above were found to be in a goodagreement with observational behaviour of turbulence under unstable conditions, butthere were obvious differences in the stable air.     

The influence of anisotropy on stress estimation by the indirect dissipation method

A comparison of observations by different authors reveals that systematic differences exist between momentum fluxes measured directly, and momentum fluxes determined indirectly by the dissipation method. This discrepancy is attributed to systematic errors due to the indirect determination of energy dissipation from the presumed inertial subrange spectrum of the horizontal wind component. The discrepancy increases with increasing degree of anisotropy, indicated by the ratio (vertical wind spectrum): (horizontal wind spectrum) deviating from % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4baFfea0dXde9vqpa0lb9% cq0dXdb9IqFHe9FjuP0-iq0dXdbba9pe0lb9hs0dXda91qaq-xfr-x% fj-hmeGabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWcaaqaaiaais% daaeaacaaIZaaaaaaa!33E6!\[\frac{4}{3}\]The results support a value of 0.48 for Kolmogoroff's constant.     

Inertial subrange correlation between temperature and humidity fluctuations in the unstable surface layer above vegetated terrains

Observations of the temperature-humidity cospectrum and correlation spectrum were made with a cold platinum wire and a Ly- hygrometer at 3.7 and 10 m above vegetated surfaces during unstable atmospheric conditions. It was found theoretically that a separation between the temperature and humidity sensors causes a drop-off of the correlation spectrum at wavenumbers > 0.3 ^–1. The observed drop-off follows the theoretical one reasonably well. Measurements made with the temperature sensor placed in the center of the Ly- gap reveal a % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaccaGae8NKby% kaaa!37B5!\[f\]^–5/3 dependence of the temperature-humidity cospectrum in the inertial subrange up to frequencies of 20 Hz. The drop-off at higher frequenties is thought to be caused by limitations inherent to the Ly- humidiometer.     

Turbulence mechanisms at an agricultural site

An extensive set of turbulence data from the 3- and 12-m heights taken over an agricultural site (Marsta, Sweden) are analyzed and compared with data from ideal sites.In unstable air, Monin-Obukhov similarity is found to be valid for the non-dimensional gradients of wind, _m , temperature, _h , and humidity, _e , for (only a few data), for _T /|T _*|,/ _E /|E _*| and for the non-dimensionalized inertial subrange spectra of temperature and humidity. Where comparison is possible, the unstable data also agree with those found in the Kansas study, with one remarkable exception, the inertial subrange constant of the temperature spectrum, ₁, being only 0.39, compared to the value 0.80 found at the Kansas site.On the stable side, most similarity predictions break down, with most of the data differing systematically from the corresponding Kansas results, the only exception being . The inertial subrange constants for temperature, ₁, and for humidity, ₁ are found to have the same values, 0.39 and 0.30, respectively, as they do on the unstable side. Remarkable similarity is found for the shape of the stable u- and - and e-spectra. In addition, this shape is found to be identical with that found in Kansas. The peak wavelength of the stable u-, and -spectra is found to be about four times larger than it is for the corresponding Kansas spectra. This is interpreted to be a result of the increased macro-roughness at the Marsta site as compared with that at the Kansas site. A possible explanation for the low ₁-value is discussed, suggesting that ₁ is not a universal constant, but instead dependent on the turbulent structure.     

Observed variations of σθ and σΦ in the nocturnal drainage flow in a deep valley

The variations of and in the drainage flow in the Brush Creek valley of western Colorado are investigated using data from Doppler acoustic sodars and instrumented towers. The data were obtained on two experimental nights during the 1984 ASCOT field study. There is good agreement between the variations derived from low-level observations of the sodars and those derived from the towers located throughout the valley. The observed hourly average and in the nocturnal drainage flow are about 20 ° to 25 ° and 5 °, respectively; these values are much larger than those generally observed over flat terrain during nighttime stable conditions. After sunrise (about 0600 MST), as the valley warms and the flow direction changes to up-valley, these parameters increase sharply to their peak values at about 0800 MST and then decrease to their normal daytime values after about two hours.In the drainage flow, the hourly average varies inversely with wind speed according to the relation u 0.7ms^-1. The vertical standard deviation is much less enhanced by complex terrain than the horizontal standard deviation. The observed values are predicted fairly well by the local similarity theory.Oak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 andOak Ridge Associated Universities (ORAU) Summer Research Participant at ATDD in 1987 and     

The validation of various schemes for parameterizing evaporation from bare soil for use in meteorological models: A numerical study using in situ data

Parameterization of evaporation from a non-plant-covered surface is very important in the hierarchy strategy of modelling land surface processes. One of the representations frequently used in its computation is the resistance formulation. The performance of the evaporation schemes using the , , and their combination resistance approaches to parameterize evaporation from bare soil surfaces is discussed. For that purpose, the nine schemes, based on a different dependence of and on volumetric soil moisture content and its saturated value, are used.The tests of performances of the considered schemes are based on time integrations by the land surface module (BARESOIL) using observed data. The 23 data sets at a bare surface experimental site in Rimski anevi, Yugoslavia on chernozem soil, were used for the resistance algorithm evaluation. The quality of the schemes was compared with the observed values of the latent heat flux using several statistical parameters.     

On Trajectory Curvature as a Selection Criterion for Valid Lagrangian Stochastic Dispersion Models

Recently Wilson and Flesch (Boundary-Layer Meteorology, 84, 411-426, 1997) suggested that the average increment d z to the orientation = arctan(w/u) of the Lagrangian velocity-fluctuation vector can be used to distinguish the better Lagrangian stochastic models within the well-mixed class. Here it is demonstrated that the specification of d z constitutes neither a sufficient or universally applicable criterion to distinguish the better Lagrangian stochastic models within the well-mixed class. The hypothesis made by Wilson and Flesch that Lagrangian stochastic models with /P_E irrotational are zero-spin models, having d z=0, is proven     

A Green Planet Versus a Desert World: Estimating the Maximum Effect of Vegetation on the Land Surface Climate

We quantify the maximum possible influence of vegetation on the global climate by conducting two extreme climate model simulations: in a first simulation (desert world), values representative of a desert are used for the land surface parameters for all non glaciated land regions. At the other extreme, a second simulation is performed (green planet) in which values are used which are most beneficial for the biosphere's productivity. Land surface evapotranspiration more than triples in the presence of the green planet, land precipitation doubles (as a second order effect) and near surface temperatures are lower by as much as 8 K in the seasonal mean resulting from the increase in latent heat flux. The differences can be understood in terms of more absorbed radiation at the surface and increased recycling of water. Most of the increase in net surface radiation originates from less thermal radiative loss and not from increases in solar radiation which would be expected from the albedo change. To illustrate the differences in climatic character and what it would imply for the vegetation type, we use the Köppen climate classification. Both cases lead to similar classifications in the extra tropics and South America indicating that the character of the climate is not substantially altered in these regions. Fundamental changes occur over Africa, South Asia and Australia, where large regions are classified as arid (grassland/desert) climate in the desert world simulation while classified as a forest climate in the green planet simulation as a result of the strong influence of maximum vegetation on the climate. This implies that these regions are especially sensitive to biosphere-atmosphere interaction.     

A Comparative Analysis of Transpiration and Bare Soil Evaporation

Transpiration E^v and bare soil evaporation E^b processes are comparatively analysed assuming homogeneous and inhomogeneous areal distributions of volumetric soil moisture content . For a homogeneous areal distribution of we use a deterministic model, while for inhomogeneous distributions a statistical-deterministic diagnostic surface energy balance model is applied. The areal variations of are simulated by Monte-Carlo runs assuming normal distributions of .The numerical experiments are performed for loam. In the experiments we used different parameterizations for vegetation and bare soil surface resistances and strong atmospheric forcing. According to the results theE^v()-^Eb() differences are great, especially in dry conditions. In spite of this, the available energy flux curves of vegetation A^v() and bare soil A^b() surfaces differ much less than the E^v() and ^Eb() curves. The results suggest that E^v is much more non-linearly related to environmental conditions than ^Eb. Both E^v and ^Eb depend on the distribution of , the wetness regime and the parameterization used. With the parameterizations, ^Eb showed greater variations than E^v. These results are valid when there are no advective effects or mesoscale circulation patterns and the stratification is unstable.     

Effects of water vapour on the structure parameter of the refractive index for near-infrared radiation

The refractive index of air (n) mainly depends on temperature andwater vapour content. For near-infrared radiation, temperature isthe main determining factor. To determinethe structure parameter of temperature ( ) fromthe structure parameter of the refractive index( ), the influence ofwater vapour content on n needs to be taken into account as acorrection. Three levels of approximation are presented. The firstinvolves the standard deviations of T and q (_T and_q) as well as the correlation coefficient betweenT and q (R_Tq). The second approximation involvesR_Tq and the Bowen ratio (), and the last uses only the Bowen ratio.The latter is the classical Bowen ratio correction.Evaluation of the validity of the assumptions used in the derivationreveals that a large error may be introduced (for situations withR_Tq < 0.8) by replacing _T/_qwith the Bowen ratio.Consequently, estimating from using the Bowen ratiocan contain errors of 5–40% for || < 1. Adding R_Tqas extra information reduces the errors only marginally. For || > 1,the correction is small, and all three approximations give errors of lessthan 1% in . When is used to compute the sensible heat flux,the influence of the quality of the correction forwater vapour fluctuations on the measured energy balance is small: forsmall ||, the correction is large, but the absolute value ofthe heat flux is small, whereas for large ||, the correctionis insignificant.     

PHYSICAL IMPACTS OF CLIMATIC CHANGE ON DELTAIC COASTAL SYSTEMS (II): DRIVING TERMS

The evaluation of the impact of climatic change on deltaic coastal systems requires establishing a reference situation. As deltaic coastal zones are highly dynamic areas, it is necessary to define the reference situation as a function of their present evolutive stage. This implies characterizing coastal processes as a complex system defined by the links between driving or forcing terms and the coastal response at different temporal and spatial scales (long/large, medium and episodic scales). Starting from classical previous works on deltaic systems and including all available field and conceptual knowledge to properly consider present conditions (largely regulated rivers, subsidence, etc.) a detailed methodology to analyse forcing or driving agents at the three mentioned scales is presented. This methodology which developed the general approach presented in a companion paper, will be followed by a future paper focusing on the analysis of the coastal response.     

A dynamical analysis of the drag conditions at the sea surface

When applied to a sea surface, shortcomings are noted for the ordinary classification of drag conditions at rigid underlying surfaces according to the Reynolds roughness number Re _s . It is shown that in the case of mobile underlying surfaces, it would be more natural to use the dynamical classification of drag conditions according to the order of magnitude of the ratio ( = /) of the momentum flux toward the waves ( _w) to the viscous momentum flux through the surface ( _w). The relevant estimates of for the main stages of development of the wind waves indicate that the observed values of the drag coefficient of the sea surface correspond to the case of underdeveloped roughness.     

On the inner-layer scale height of boundary-layer flow over low hills

Jackson and Hunt's (1975) equation for the depth of the inner layer of flows over low hills does not depend on any closure assumption as contrarily supposed in literature. This equation contains a constant which can arbitrarily be specified. It is suggested that this inner-layer constant should be determined from experimental data. A preliminary check with some data from the Askervein experiment suggests that Jackson and Hunt's equation fits these data almost as well as Jensen's equation provided that fitted inner-layer constants are used.     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.