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.     

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.     

An explicit three-dimensional nonhydrostatic numerical simulation of a tropical cyclone

Summary A nonhydrostatic numerical simulation of a tropical cyclone is performed with explicit representation of cumulus on a meso- scale grid and for a brief period on a meso- scale grid. Individual cumulus plumes are represented by a combination of explicit resolution and a 1.5 level closure predicting turbulent kinetic energy (TKE).The results demonstrate a number of expected and unexpected important scale interaction processes. Within the central core of the developing cyclone, meso- convective regions grow and breakdown into propagating inertiagravity waves throughout the lifecycle of the cyclone. In the early stages, the amplitude of pressure fluctuations associated with the meso- scale convection exceed the central pressure of the cyclone and strongly modulate its intensity. With each meso- scale pulsation, the cyclone core increases in strength, measured by the central pressure deficit. The increasingly strong inertial frequency of the storm core acts to increasingly trap the convection induced heating within the core by balancing the tangential wind against the low central pressure, before the meso- scale convection breaks down and sends the warmth away as a propagating wave. Eventually, the slow manifold's amplitude exceeds the amplitude of the meso- scale oscillations and a stable eye region is formed. As inertial instability increases, increasingly high thermal warmth can be protected in the core, allowing persistent subsidence to form and to clear out the cyclone eye.On the outside of the eye wall, strong inertial stability gradients in the troposphere cause convective warming to split the inflow to the eye wal! and spawn outwardly propagating inertia gravity waves. These waves carry away all of the heating forced by convection that is not inertially trapped by the eye wall and act as a moderating influence on storm intensity.Inertia gravity waves are also spawned in the stratosphere at the top of the eye wall by the revolution of asymmetric cumulus structures. In all instances, the tropospheric waves are coupled to the propagating stratospheric waves which both move at 35 ms^–1, although there are many instances where the stratospheric waves seem to have no tropospheric counterpart. Hence the anvil top forcing and low level breakdown are linked.The outwardly propagating inertia gravity waves act to initiate outer bands of convection. This initiation is with the assistance of low level boundary layer variations of density related to previous convection and to virga falling from the anvil which moistens and destabilizes the mid levels of _e minimum. The convection initiated by these waves does not move substantially outward with the wave, although may appear to develop outward discontinuously.With 12 Figures     

Relation of Richardson number to the curvature of the wind profile

Values of Richardson number (Ri) and the wind profile curvature () in the atmospheric surface layer are evaluated from experimental data. The relationship between Ri and shows some scatter, but gives agreement with the theoretical relation derived by Businger et al. (1971).     

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.     

Examination of the empirical flux-profile models in the atmospheric surface boundary layer

Values of Richardson number (Ri) and the wind profile Deacon number () in the atmospheric surface layer are evaluated from various experimental data sources. The relationship between Ri and shows some scatter, but gives agreement with the relations derived by the several empirical models available for both unstable and inversion conditions.     

Equilibrium evaporation beneath a growing convective boundary layer

Expressions for the equilibrium surface Bowen ratio ( _s ) and equilibrium evaporation are derived for a growing convective boundary layer (CBL) in terms of the Bowen ratio at the top of the mixed layer _i and the entrainment parameter A _R . If A_R is put equal to zero, the solution for _s becomes-that previously obtained for the zero entrainment or closed box model. The Priestley-Taylor parameter is also calculated and plotted in terms ofA _R and _i . Realistic combinations of the atmospheric parameters give values of in the range 1.1 to 1.4.     

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.     

A study of the sensitivity of bare soil evaporation schemes to soil surface wetness,using the coupled soil moisture and surface temperature prediction model,BARESOIL

Summary The performance of evaporation schemes with and approach and their combination within resistance representation of evaporation from bare soil surface is discussed. For this purpose nine schemes, based on different functions of or , on the ratio of the volumetric soil moisture content and its saturated value are used.The quality of the chosen schemes has been evaluated using the results of time integration by the coupled soil moisture and surface temperature prediction model, BARESOIL, using in situ data. A sensitivity analysis was made using two sets of data derived from the volumetric soil moisture content of the top soil layer. One with values below the wilting point (0.17 m³m^–3) and the second with values above 0.20m³m^–3. Data sets were obtained at the experimental site Rimski anevi, Yugoslavia, from the bare surface of a chernozem soil.With 4 Figures     

Gas and Particle Products Distribution from the Reaction of β-Caryophyllene with Ozone

Gas and particulate reaction products from the ozonolysis of -caryophyllene (I) in the presence of atmospheric air were investigated using a combination of gas chromatography-mass spectrometry (GC-MS) and high performance liquid chromatography (HPLC). A Scanning Mobility Particle Sizer system (3936, TSI) and a Condensation Particle Counter (3025A, TSI) were used to study secondary organic aerosol formation. The nighttime oxidation was carried out in a large outdoor smog chamber (190 m³). A wide range of ring retaining and ring opening products in the gas and particle phase are reported over the course of the reaction. On average, measured gas and particle phase products accounted for 64% of the reacted -caryophyllene (I) carbon. Measurements show that a number of reaction products with low vapor pressure (e.g., -caryophyllone aldehyde (IV), -norcaryophyllone aldehyde (V), -caryophyllonic acid (VIII), -14-hydroxycaryophyllonic acid (XIV)) were found in the sample taken during the first 20 min of the reaction and may play an important role in the early formation of secondary organic aerosol. A detailed mechanism is proposed to account for most products observed in this investigation.     

Influence of heterogeneous land surfaces on surface energy and mass fluxes

Summary Land-surface heterogeneity affects surface energy fluxes. The magnitudes of selected land-surface influences are quantified by comparing observations with model simulations of the FIFE (First ISLSCP Field Experiment) domain. Several plausible heterogeneous and homogeneous initial and boundary conditions are examined, although soilmoisture variability is emphasized. It turns out that simple spatial averages of surface variation produced biased flux values. Simulated maximum latent-heat fluxes were approximately 30 to 40 W m^–2 higher, and air temperatures 0.4 °C lower (at noon), when computations were initialized with spatially averaged soil-moisture and leaf-area-index fields. The planetary boundary layer (PBL) height and turbulent exchanges were lower as well. It additionally was observed that (largely due to the nonlinear relationship between initial soil-moisture availability and the evapotranspiration rate), real latent-heat flux can be substantially less than simulated latent-heat flux using models initialized with spatially averaged soil-moisture fields. Differences between real and simulated fluxes also vary with the resolution at which real soil-moisture heterogeneity is discretized.With 8 Figures     

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.     

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.     

Evaporation from a blanket bog in a foggy coastal environment

Frequent fog severely restricts evaporation from blanket bogs in Newfoundland because it more than halves the radiant energy input, and it eliminates the vapor pressure deficit, resulting in evaporation at the equilibrium rate (average = 0.99 during fog). During these periods, there is no surface resistance to evaporation because the bog has been wetted by fog drip, and although the latent heat flux dominates over sensible heat (average = 0.8), both are small. In contrast, the surface dries during clear periods, increasing the surface resistance to evaporation so that sensible heat becomes more important ( = 1.05). When the mosses are dry, evaporation is below the equilibrium rate ( = 0.87), although the higher available energy ensures that actual evaporation is higher. During clear periods, daily evaporation averaged 2.5 mm, compared to 1.1 and 0.7 mm for fog and rain, respectively. The suppressed evaporation at this site is important in maintaining appropriate hydrological conditions for blanket bog development.     

A resistance-law hypothesis for the non-stationary advective planetary boundary layer

The theoretical developments yielding the resistance law for the planetary boundary layer are summarized, including the refinements of the ageostrophic method. This leads to the hypothesis that the resistance law for the simple non-stationary, advective boundary layer is obtained from the one for a stationary, horizontally homogeneous boundary layer, if the modulous of the geostrophic wind in the latter is replaced with a generalized frictionless wind in the former. The generalized frictionless wind includes the inertial terms.Contribution from the Sonderforschungsbereich Meeresforschung Hamburg of the Deutsche Forschungsgemeinschaft.     

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.     

The interaction of binary vortices in a barotropic model

Summary The interaction of binary cyclonic vortices is investigated using the nondivergent barotropic model of Chan and Williams (1987) under two situations: a quiescent environment and a linearly-sheared background flow. It is found that the mutual interaction between the vortices results from a combination of two processes: the advection of symmetric vorticity by the asymmetric flow and the advection of asymmetric vorticity by the symmetric flow. The latter contribution is rather significant. Whether the vortices in a binary system attract or repel each other depends on the asymmetric vorticity distribution associated with the two vortices. Such a distribution is governed by the structure (size) of and the separation between the vortices. In the presence of a sheared flow, the contribution from the advection of asymmetric vorticity by the symmetric flow may also become appreciable depending on the structure and magnitude of the shear. Furthermore, the geographical locations of the vortices in relation to the sheared flow are also important in determining the relative movement of the vortices.In the presence of , the movements of the vortices are modified by the northwestward -drift However, the relative motion between the vortices is almost identical to that on an f-plane. In other words, the mutual interaction between the vortices is largely independent of . Alternatively, the two vortices can be considered to be one system which drifts towards the northwest under the influence of while they interact with each other within the system. Physically, this independence arises because the two relative vorticity advection terms have much larger magnitudes than the planetary vorticity advection term. However, the -effect is still important in that it modifies the asymmetric flow associated with each vortex and hence the asymmetric vorticity. Such modifications change the advection patterns compared with the =0 case and hence lead to different vortex movements. The presence of a linear shear causes the binary system to move as if it was a large (for a cyclonic shear) or smaller (for an anticyclonic shear) vortex under the influence of .With 22 Figures     

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.     

On some new statistical characteristics of the wind “Koshava”

Summary Koshava is a gusty wind of moderate to strong intensity, blowing from the south-eastern direction, over the area of the Republic of Serbia. It is caused by the interaction between the synoptic circulation and the orography of the Carpathian and the Balkan mountains. The Koshava wind can damage buildings, factories and industrial plants or city infrastructure. Therefore it is important to estimate its gust and the gustiness factor on the basis of the measured data.This paper discusses a statistical analysis of wind data in the maximum influence area of the Koshava wind in the periods of maximum duration of Koshava. The focus of the paper is the examination of urban and suburban effects on Koshava wind and the correlation between the instantaneous maximum wind speed and the hourly mean wind speed. The best fitting with various empirical distributions is proposed.With 10 Figures