Experiments On Buoyant Plume Dispersion In A Laboratory Convection Tank

Plume dispersion in the convective boundary layer (CBL) is investigated experimentally in a laboratory convection tank. The focusis on highly-buoyant plumes that loft near or become trapped in the CBL capping inversion and resistdownward mixing. Such plumes are defined by dimensionless buoyancy fluxes F_* 0.1, where F_* = F_b/(U w_* ² z_i), F_b is the stack buoyancy flux,U is the mean wind speed, w_* is the convective velocity scale, and z_i is the CBL depth. The aim is to obtain statistically-reliable mean (C) and root-mean-square (rms, _c) concentration fields as a function of F_* and the dimensionless distance X = w_*x/(U z_i), where x is the distance downstream of the source.The experiments reveal the following mainresults: (1) For 3 X 4and F_* 0.1, the crosswind-integrated concentration (CWIC) fields exhibit distinctly uniform profiles below z_i with a CWIC maximum aloft, in contrast to the nonuniform profiles obtained earlier by Willis and Deardorff. (2) The lateral dispersion (_y) variation with X is consistent with Taylor's theory for _* 0.1 and a buoyancy-enhanced dispersion, _y/z_i F_* ^1/3X^2/3, forF_* = 0.2 and 0.4. (3) The entrapment, the plume fraction above z_i, has a mean (E) that follows a systematic variationwith X and F_*, and a variability (_e/E) that is broad ( 0.3 to 2) near the source but subsides to 0.25 far downstream. (4) Vertical profiles of the concentration fluctuation intensity (c/C) are uniform for z < z_i and X > 1.5, but exhibit significant increases: (a) at the surface and close to the source (X 1.5), and(b) in the entrainment zone. (5) The cumulative distribution functions (CDFs) of the scaled concentration fluctuations (c/_c) separate into mixed-layer and entrainment-layer CDFs for X 2, with the mixed-layer group collapsing to a single distribution independent of z.These are the first experiments to obtain all components of the lateral and vertical dispersion parameters (rms meander, relative dispersion, total dispersion) for continuous buoyant releases in a convection tank. They also are the first tank experiments to demonstrate agreement with field observations of: (1) the scaled ground-level concentration along the plume centreline, and (2) the dimensionless lateral dispersion _y/z_i of buoyant plumes.     

Interactions between coherent eddies in the lower convective boundary layer

The processes of interaction between the atmospheric surface and mixed layers in daytime convective conditions over land are studied using a data set obtained during flights by an instrumented aircraft. Profiles of partitioned run-averaged statistics and examples of time series plots are discussed in the light of results from a recently published study by the authors, in which the average structure and flow within coherent eddies was reconstruced using a compositing technique. This evidence is used to support a conceptual model of the mechanisms of interaction between surface-layer plumes and mixed-layer thermal columns. The divergent flow created near the surface by the downdraft arms of the large-scale mixed-layer circulation patterns, forces the development of lines of convergence in the surface layer (the so-called thermal walls), which channel air into the bases of the mixed-layer thermals. Plumes progressively group and merge together with height in the surface and free convection layers, and move along these convergence lines toward large collector plumes at the intersection points, or hubs. Above the hubs are the thermals, and air parcels originating from plumes and their environment are strongly mixed as they rise, leading to an increased difficulty of the conditional sampling method to distinguish between them. The observed influence of mixed-layer convective processes far down into the surface layer, and the form of the averaged profiles, supports recent refinements of the theory of surface-layer structure suggested in Kader and Yaglom (1990).Notation CBL convective boundary layer - SL surface layer - FCL free convection layer - ML mixed layer     

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.     

The Interference of Higher-Order Statistics of the Concentration Field Produced by Two Point Sources According to a Generalized Fluctuating Plume Model

The higher-order correlation functions for the concentrationfluctuations arising from a two-point-source configuration have beencalculated analytically within the context of the phenomenology of afluctuating plume model (viz., a meandering plume model that explicitlyincorporates internal fluctuations). Explicit expressions for thesecond-, third-, and fourth-order correlationfunctions between the concentrationfluctuations produced by two point sources are given in terms of the sourceseparation d and the five physically based parameters that define thegeneralized fluctuating plume model: namely, the absolute plume dispersion,_a, which determines the outer plume length scale; the relative plume dispersion, _r, which determines the inner plume length scale; the fluctuation intensity, i_r, in relative coordinates, which determines the internal concentration fluctuation level; the correlation coefficient, r,between the positions of the centroids of the two interfering plumes; and,the correlation coefficient, r^*, between the concentration fluctuationsof the two plumes in relative coordinates, which determines the degree ofinternal mixing of the two scalars. Furthermore, the form of the totalconcentration probability density function arising from the interferenceproduced by two point sources is presented. Predictions for the second-ordercorrelation function, , and for the total concentration probabilitydensity function have been compared with some new experimental data fora two-point-source configuration in grid turbulence generated in awater-channel simulation. These results are in good agreement with the dataand suggest that the analytical model for the second-order correlationfunction and the total concentration probability density function canreproduce many qualitative trends in the interaction of plumes from twosources.     

Turbulence characteristics and organized motion in a suburban roughness sublayer

To investigate tubulence characteristics and organized motion within and above an urban canopy, field observations were conducted in July 1991 and November 1992 in Sapporo, Japan. The measurement heights were 5.4, 10.3, 18, 35 and 45 m above ground; the canopy height was 7 m. The profiles of _u peaked slightly above the canopy, while _v and _w had nearly uniform profiles. Vertical profiles of Reynolds stress - peaked slightly at 1.5 times the canopy height and decreased slowly with height thereafter. A four-quadrant analysis showed that sweep and ejection motions caused high-velocity fluid from above moves downward toward the surface and low-velocity fluid from below moves upward. An ensemble-averaging technique was used to isolate typical features of the flow and temperature fields. A time-height cross-section of velocity vectors and temperature contours showed details of the flow structures associated with temperature ramps. It has been noted that the organized motions play important roles in the transport of heat near the urban canopy, where the sweep motion causes negative temperature fluctuations and the ejection motion causes positive temperature fluctuations.     

An Alternative Function For The Wind And Temperature Gradients In Unstable Surface Layers

The function ()=(1+|z/L|^2/3)^1/2,where z is the height, L the Obukhov length, and a constant,is proposed for the nondimensional wind speed and temperaturegradients (flux-profile relationships) in anunstable surface layer. This function agrees quite well withboth wind speed and temperature data,has the theoretically correct behaviour in convective conditions,and leads to simple results when integrated to produce the mean profiles.     

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.     

Transport-Dissipation Analytical Solutions to the E-∈Turbulence Model and their Role in Predictions of the Neutral ABL

E- turbulence model predictions of the neutralatmospheric boundary layer (NABL) are reinvestigated to determine thecause for turbulence overpredictions found in previous applications. Analytical solutions to the coupled E and equations for the case of steady balance between transport and dissipation terms, the dominant balance just below the NABL top, are derived. It is found that analytical turbulence profiles laminarizeat a finite height only for values of closure parameter ratio c ₂ /_e equal toor slightly greater than one, with laminarization as z for greater . The point = 2 is additionally foundthat where analytical turbulent length scale (l) profilesmade a transition from ones ofdecreasing ( < 2) to increasing ( > 2)values with height. Numerically predicted profiles near the NABL topare consistent with analytical findings. The height-increasingvalues of l predicted throughout the NABL with standard values ofclosure parameters thus appear a consequence of 2.5(> 2), implied by these values (c ₂ = 1.92, _{= 1.3}, _{e = 1}). Comparison of numericalpredictions with DNS data shows that turbulence overpredictions obtained with standard-valued parameters are rectifiedby resetting and _e to 1.1 and 1.6, respectively, giving, with c ₂ = 1.92, 1.3, and laminarization of the NABL's cappingtransport-dissipation region at a finite height.     

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.     

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.     

Numerische Integrationen einer Gleichung für die turbulente Diffusion in der Atmosphäre

Zusammenfassung Es werden numerische Integrationen der dreidimensionalen Fickschen Diffusionsgleichung mit höhenabhängigen Diffusionskoeffizienten und Windgeschwindigkeit für den stationären Fall und eine kontinuierlich emittierende Punktquelle durchgeführt. Für die drei typischen Ausbreitungsverhältnisse Starke Durchmischung, Bodeninversion und Inversion über der Quelle, charakterisiert durch die Diffusionskoeffizienten und die Windgeschwindigkeit, werden die Einflüsse unterschiedlicher Quellhöhe, Sinkgeschwindigkeit und Ablagerungsgeschwindigkeit studiert.

---

Summary Numerical integrations of the threedimensional Fickean diffusion equation with height dependent diffusion coefficients and wind velocity are performed in the stationary case with a continuous point source. For the three typical cases strong mixing, surface inversion and inversion above the source characterized by the diffusion coefficients and wind velocity, the effects of source height, sinking velocity and deposition velocity are discussed.

---

Résumé Les auteurs ont procédé à l'intégration numérique de l'équation à trois dimensions de la diffusion, équation établie par Fick. Ils ont, pour cela, utilisé des coefficients de diffusion variant avec l'altitude et la vitesse du vent et distingué l'état stationnaire et le cas d'une source d'émission ponctuelle et continue. Ils étudient en outre les influences de hauteurs différentes de la source, de la rapidité de subsidence et de la vitesse de sédimentation. Ils ont tenu compte pour cela des trois conditions typiques de diffusion mélange très accentué, inversion au sol et inversion au-dessus de la source.

---

---

Mit 7 Textabbildungen     

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.     

Occurrence of roll circulations in a shallow boundary layer

Meteorological measurements taken at the Näsudden wind turbine site during slightly unstable conditions have been analyzed. The height of the convective boundary layer (CBL) was rather low, varying between 60 and 300 m. Turbulence statistics near the ground followed Monin-Obukhov similarity, whereas the remaining part of the boundary layer can be regarded as a near neutral upper layer. In 55% of the runs, horizontal roll vortices were found. Those were the most unstable runs, with -z _i/L > 5. Spectra and co-spectra are used to identify the structures. Three roll indicators were identified: (i) a low frequency peak in the spectrum of the lateral component at low level; (ii) a corresponding increase in the vertical component at mid-CBL; (iii) a positive covariance {ovvw} together with positive wind shear in the lateral direction (V/z) in the CBL. By applying these indicators, it is possible to show that horizontal roll circulations are likely to be a common phenomenon over the Baltic during late summer and early winter.     

Infrared observations and numerical modelling of grassland fires in the Northern Territory, Australia

Summary This paper reports on a small-scale pilot experiment held early in the dry season near Darwin, Australia, in which fine-scale observations of several prescribed fires were made using infrared digital video. Infrared imaging is used routinely to locate fires as infrared radiation suffers little attenuation as it propagates through the smoke that normally obscures visible imagery. However, until now, little use has been made of digital video imagery in analyzing the convective-scale structure of prescribed (or wild) fires. The advantage of digital video imagery is that the individual frames can be objectively analyzed to determine the convective motion in the plane viewed by the camera. The infrared imagery shows mostly rising plumes, much like convective clouds. The flow is highly convective, and the vertical transport of heat is confined to relatively narrow thermals. The updrafts range from a few ms^–1 to around 15ms^–1. A numerical model is used to simulate one of the prescribed fires at very high-resolution. For the most part, the model predictions compare well to the observations. The model produces plumes that are around 7m high, and spaced around 5m apart, which is similar to that observed. The model correctly predicts the mean rate of spread of the fire to be 1.3ms^–1. Perhaps the most serious limitations to using infrared observations of the type presented here are the difficulties in interpreting precisely the relationship between the observed infrared temperature field and the air temperature calculated by the model, and the exact connection between the infrared camera derived flow field and that calculated by the model.     

Application of a “Big-Tree” model to regional climate modeling: a sensitivity study

Summary ¶In order to better understand land-atmosphere interactions and increase the predictability of climate models, it is important to investigate the role of forest representation in climate modeling. Corresponding to the big-leaf model commonly employed in land surface schemes to represent the effects of a forest, a so called big-tree model, which uses multi-layer vegetation to represent the vertical canopy heterogeneity, was introduced and incorporated into the National Center for Atmospheric Research (NCAR) regional climate model RegCM2, to make the vegetation model more physically based. Using this augmented RegCM2 and station data for China during 1991 Meiyu season, we performed 10 experiments to investigate the effects of the application of the big-tree model on the summer monsoon climate.With the big-tree model incorporated into the regional climate model, some climate characteristics, e.g. the 3-month-mean surface temperature, circulation, and precipitation, are significantly and systematically changed over the model domain, and the change of the characteristics differs depending on the area. Due to the better representation of the shading effect in the big-tree model, the temperature of the lower layer atmosphere above the plant canopy is increased, which further influences the 850hPa temperature. In addition, there are significant decreases in the mean latent heat fluxes (within 20–30W/m²) in the three areas of the model domain.The application of the big-tree model influences not only the simulated climate of the forested area, but also that of the whole model domain, and its impact is greater on the lower atmosphere than on the upper atmosphere. The simulated rainfall and surface temperature deviate from the originally simulated result and are (or seem to be) closer to the observations, which implies that an appropriate representation of the big-tree model may improve the simulation of the summer monsoon climate.We also find that the simulated climate is sensitive to some big-tree parameter values and schemes, such as the shape, height, zero-plane displacement height and mixing-length scheme. The simulated local/grid differences may be very large although the simulated areal-average differences may be much lower. The area-average differences in the monthly-mean surface temperature and heat fluxes can amount to 0.5°C and 4W/m², respectively, which correspond to maximum local/grid differences of 3.0°C and 40W/m² respectively. It seems that the simulated climate is most sensitive to the parameter of the zero-plane displacement among the parameters studied.     

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.     

Turbulence in a Katabatic Flow

Turbulence measurements performed in a stable boundary layer over the sloping ice surface of the Vatnajökull in Iceland are described. The boundary layer, in which katabatic forces are stronger than the large-scale forces, has a structure that closely resembles that of a stable boundary layer overlying a flat land surface, although there are some important differences. In order to compare the two situations the set-up of the instruments on an ice cap in Iceland was reproduced on a flat grass surface at Cabauw, the Netherlands. Wind speed and temperature gradients were calculated and combined with flux measurements made with a sonic anemometer in order to obtain the local stability functions _m and _h as a function of the local stability parameter z/L. Unlike the situation at Cabauw, where _m was linear as a function of z/L, in the katabatically forced boundary layer, the dependence of _m on stability was found to be non-linear and related to the height of the wind maximum. Thermal stratification and the depth of the stable boundary layer however seem to be rather similar under these two different forcing conditions.Furthermore, measurements on the ice were used to construct the energy balance. These showed good agreement between observed melt and components contributing to the energy balance: net radiation (supplying 55% of the energy), sensible heat flux (30%) and latent heat flux (15%).Local sources and sinks in the turbulent kinetic energy budget are summed and indicate a reasonable balance in near-neutral conditions but not in more stable situations. The standard deviation of the velocity fluctuations _u, _v, and _w, can be scaled satisfactorily with the local friction velocity u_* and the standard deviation of the temperature fluctuation with the local temperature scale _*.     

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.     

Statistical properties of optical refractive index fluctuations in the marine boundary layer

The optical refractive index fluctuation has been determined from measurements of turbulent temperature and humidity fluctuations in the atmospheric surface layer over the ocean. Probability density, spectral density and even- and odd-order structure functions of the refractive index fluctuation are presented. The spectral density exhibits a significant -5/3 inertial subrange as a result of the existence of a -5/3 subrange in the spectra of temperature and humidity fluctuations and in the temperature-humidity co-spectrum. The behaviour in the inertial subrange of fourth- and sixth-order structure functions of the refractive index is in reasonable agreement with that predicted by the analysis of Antonia and Van Atta (1975, 1978). The third-order structure functions of the refractive index exhibit an approximately linear dependence on separation in the inertial subrange, in agreement with a ramp model for the large coherent structure of the motion.     

On the similarity in atmospheric fluctuations of carbon dioxide,water vapor and temperature over vegetated fields

This paper describes the similarity between atmospheric fluctuations of carbon dioxide, water vapor and temperature using data which cover a wide range of instability (0.02 < < 10). The is the Monin-Obukhov stability parameter including the humidity effect.The spectral analysis shows that the coherency between fluctuations of carbon dioxide and water vapor or temperature is very close to unity, and the phase difference is basically out of phase for whole frequency ranges analyzed. The stability dependence of the normalized standard deviation of carbon dioxide is very similar to those of water vapor and temperature. The normalized standard deviation is about 2.5 under near neutral conditions, and it decreases with increasing instability following the -1/3; power law as (-)^-1/3. The skewness factors of carbon dioxide, water vapor and temperature show a systematic departure with increasing instabilities for 0.02 < s- < 1, and level off at high instabilities for 1 < -\s < 10. The stability dependence of the flatness factors is not so clear as that noted in the standrard deviation and skewness factors. Dissipation rates of carbon dioxide, water vapor and temperature variance are well related to the spectral peak wavelength. This seems to be real since the local production and local dissipation rates are the main terms, almost balancing one another in the variance budget equations for scalar entities.