Dispersion of a Passive Scalar Within and Above an Urban Street Network

The transport of a passive scalar from a continuous point-source release in an urban street network is studied using direct numerical simulation (DNS). Dispersion through the network is characterized by evaluating horizontal fluxes of scalar within and above the urban canopy and vertical exchange fluxes through the canopy top. The relative magnitude and balance of these fluxes are used to distinguish three different regions relative to the source location: a near-field region, a transition region and a far-field region. The partitioning of each of these fluxes into mean and turbulent parts is computed. It is shown that within the canopy the horizontal turbulent flux in the street network is small, whereas above the canopy it comprises a significant fraction of the total flux. Vertical fluxes through the canopy top are predominantly turbulent. The mean and turbulent fluxes are respectively parametrized in terms of an advection velocity and a detrainment velocity and the parametrization incorporated into a simple box-network model. The model treats the coupled dispersion problem within and above the street network in a unified way and predictions of mean concentrations compare well with the DNS data. This demonstrates the usefulness of the box-network approach for process studies and interpretation of results from more detailed numerical simulations.     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

Nocturnal Intermittent Coupling Between the Interior of a Pine Forest and the Air Above It

A 1-year set of measurements of CO₂ and energy turbulent fluxes above and within a 25-m pine forest in southern Brazil is analyzed. The study focuses on the coupling state between two levels and its impact on flux determination by the eddy-covariance method. The turbulent series are split in their typical temporal scales using the multiresolution decomposition, a method that allows proper identification of the time scales of the turbulent events. Initially, four case studies are presented: a continually turbulent, a continually calm, a calm then turbulent, and an intermittent night. During transitions from calm to turbulent, large scalar fluxes of opposing signs occur at both levels, suggesting the transference of air accumulated in the canopy during the stagnant period both upwards and downwards. Average fluxes are shown for the entire period as a function of turbulence intensity and a canopy Richardson number, used as an indicator of the canopy coupling state. Above the canopy, CO₂ and sensible heat fluxes decrease in magnitude both at the neutral and at the very stable limit, while below the canopy they increase monotonically with the canopy Richardson number. Latent heat fluxes decrease at both levels as the canopy air becomes more stable. The average temporal scales of the turbulent fluxes at both levels approach each other in neutral conditions, indicating that the levels are coupled in that case. Average CO₂ fluxes during turbulent periods that succeed very calm ones are appreciably larger than the overall average above the canopy and smaller than the average or negative within the canopy, indicating that the transfer of air accumulated during calm portions at later turbulent intervals affects the flux average. The implications of this process for mean flux determination are discussed.     

Investigating a Hierarchy of Eulerian Closure Models for Scalar Transfer Inside Forested Canopies

Modelling the transfer of heat, water vapour, and CO₂ between the biosphere and the atmosphere is made difficult by the complex two-way interaction between leaves and their immediate microclimate. When simulating scalar sources and sinks inside canopies on seasonal, inter-annual, or forest development time scales, the so-called well-mixed assumption (WMA) of mean concentration (i.e. vertically constant inside the canopy but dynamically evolving in time) is often employed. The WMA eliminates the need to model how vegetation alters its immediate microclimate, which necessitates formulations that utilize turbulent transport theories. Here, two inter-related questions pertinent to the WMA for modelling scalar sources, sinks, and fluxes at seasonal to inter-annual time scales are explored: (1) if the WMA is to be replaced so as to resolve this two-way interaction, how detailed must the turbulent transport model be? And (2) what are the added predictive skills gained by resolving the two-way interaction vis-à-vis other uncertainties such as seasonal variations in physiological parameters. These two questions are addressed by simulating multi-year mean scalar concentration and eddy-covariance scalar flux measurements collected in a Loblolly pine (P. taeda L.) plantation near Durham, North Carolina, U.S.A. using turbulent transport models ranging from K-theory (or first-order closure) to third-order closure schemes. The multi-layer model calculations with these closure schemes were contrasted with model calculations employing the WMA. These comparisons suggested that (i) among the three scalars, sensible heat flux predictions are most biased with respect to eddy-covariance measurements when using the WMA, (ii) first-order closure schemes are sufficient to reproduce the seasonal to inter-annual variations in scalar fluxes provided the canonical length scale of turbulence is properly specified, (iii) second-order closure models best agree with measured mean scalar concentration (and temperature) profiles inside the canopy as well as scalar fluxes above the canopy, (iv) there are no clear gains in predictive skills when using third-order closure schemes over their second-order closure counterparts. At inter-annual time scales, biases in modelled scalar fluxes incurred by using the WMA exceed those incurred when correcting for the seasonal amplitude in the maximum carboxylation capacity (V _{cmax, 25}) provided its mean value is unbiased. The role of local thermal stratification inside the canopy and possible computational simplifications in decoupling scalar transfer from the generation of the flow statistics are also discussed.

“The tree, tilting its leaves to capture bullets of light; inhaling, exhaling; its many thousand stomata breathing, creating the air”. Ruth Stone, 2002, In the Next Galaxy

     

A lagrangian random-walk model for simulating water vapor,CO2 and sensible heat flux densities and scalar profiles over and within a soybean canopy

An integrated canopy micrometeorological model is described for calculating CO₂, water vapor and sensible heat exchange rates and scalar concentration profiles over and within a crop canopy. The integrated model employs a Lagrangian random walk algorithm to calculate turbulent diffusion. The integrated model extends previous Lagrangian modelling efforts by employing biochemical, physiological and micrometeorological principles to evaluate vegetative sources and sinks. Model simulations of water vapor, CO₂ and sensible heat flux densities are tested against measurements made over a soybean canopy, while calculations of scalar profiles are tested against measurements made above and within the canopy. The model simulates energy and mass fluxes and scalar profiles above the canopy successfully. On the other hand, model calculations of scalar profiles inside the canopy do not match measurements.The tested Lagrangian model is also used to evaluate simpler modelling schemes, as needed for regional and global applications. Simple, half-order closure modelling schemes (which assume a constant scalar profile in the canopy) do not yield large errors in the computation of latent heat (LE) and CO₂ (F _c ) flux densities. Small errors occur because the source-sink formulation of LE andF _c are relatively insensitive to changes in scalar concentrations and the scalar gradients are small. On the other hand, complicated modelling frames may be needed to calculate sensible heat flux densities; the source-sink formulation of sensible heat is closely coupled to the within-canopy air temperature profile.     

Scalar Fluxes Near a Tall Building in an Aligned Array of Rectangular Buildings

Scalar dispersion from ground-level sources in arrays of buildings is investigated using wind-tunnel measurements and large-eddy simulation (LES). An array of uniform-height buildings of equal dimensions and an array with an additional single tall building (wind tunnel) or a periodically repeated tall building (LES) are considered. The buildings in the array are aligned and form long streets. The sensitivity of the dispersion pattern to small changes in wind direction is demonstrated. Vertical scalar fluxes are decomposed into the advective and turbulent parts and the influences of wind direction and of the presence of the tall building on the scalar flux components are evaluated. In the uniform-height array turbulent scalar fluxes are dominant, whereas the tall building produces an increase of the magnitude of advective scalar fluxes that yields the largest component. The presence of the tall building causes either an increase or a decrease to the total vertical scalar flux depending on the position of the source with respect to the tall building. The results of the simulations can be used to develop parametrizations for street-canyon dispersion models and enhance their capabilities in areas with tall buildings.     

Intermittency and the Exchange of Scalars in the Nocturnal Surface Layer

The determination of nocturnal surface fluxes in low wind conditions is a major problem for micrometeorological studies. The eddy correlation technique, extensively used in field measurements, becomes inappropriate if not enough turbulent activity exists. At the same time, the phenomenon of turbulence intermittency is responsible for the existence of localized events of short duration within which a large fraction of the total nighttime scalar exchange occurs. The scalar flux within a certain intermittent event varies considerably depending on the window used for the flux calculation. In many cases, events with very different time durations occur in the same night, and therefore, the proper determination of the surface flux would require averaging within data windows of different sizes for each event. In this work, the surface exchanges of temperature, moisture and carbon dioxide are analysed at a micrometeorological tower at southern Brazil. Intermittent turbulence is a common occurrence at the location. The analysis shows that the fluxes vary with turbulence intensity and the estimation technique. A variable-window size method for flux estimation is suggested and shown to cause an increase in the magnitude of the nocturnal surface fluxes     

Testing of the bandpass eddy covariance method for a long-term measurement of water vapour flux over a forest

The bandpass eddy covariance method has been used to measure the turbulent flux of scalar quantities using a slow-responsescalar sensor. The method issimilar in principle to the traditional eddy correlation method but includes the estimation of high-frequency components of the flux on the basis of cospectral similarity in the atmospheric surface layer. In order to investigate the performance of the method, measurements of the water vapour flux over a forest with the bandpass eddy covariance method and the direct eddy correlation method were compared. The flux obtained by the bandpass eddy covariance method agreed with that by the eddy correlation method within ±20% for most cases, in spite of a rather slow sensor-response of the adopted hygrometer. This result supports its relevance to a long-term continuous operation, since a stable, low-maintenance,general-purpose sensor canbe utilized for scalar quantities. Oneweak point of the method isits difficulty in principle to measure the correct flux when the magnitude of the sensible heat flux is very small, because the method uses the sensible heat flux as a standard reference for the prediction of undetectable high-frequency components of the scalar flux. An advanced method is then presented to increase its robustness. In the new method, output signals from a slow-response sensor are corrected using empirical frequency-responsefunctions for the sensor,thereby extending the width of the bandpass frequency region where components of the flux are directly measured (not predicted). The advanced method produced correct fluxes for all cases including the cases of small sensible heat flux. The advanced bandpass eddy covariance method is thus appropriate for along-term measurement of the scalar fluxes.     

Local advection of momentum,heat, and moisture in micrometeorology

The local advection of momentum, heat and moisture in micrometeorology due to a horizontal inhomogeneity in surface conditions is numerically investigated by a higher-order turbulence closure model which includes equations for the mean quantities, turbulent fluxes, and the viscous dissipation rate. The application of the two-dimensional model in this paper deals with the simulation of the flow from an extensive smooth dry area to a grassy wet terrain. The mean wind speed, temperature, and humidity distributions in the resulting internal boundary layer downstream of the surface discontinuity are determined such that the energy and moisture balances at the Earth's surface are satisfied.Numerical calculations of the mean temperature and humidity profiles are compared with available observed ones. The results include the advective effects on turbulent flux distributions, surface energy balance, evaporation rate, and Bowen ratio. The sensitivity of the predicted mean profiles and turbulent flux distributions to the surface relative humidity, thermal stratification, and the roughness change is discussed.NRC-NAS Resident Research Associate at AFCRL.     

Coherent Structures and the Dissimilarity of Turbulent Transport of Momentum and Scalars in the Unstable Atmospheric Surface Layer

Atmospheric stability effects on the dissimilarity between the turbulent transport of momentum and scalars (water vapour and temperature) are investigated in the neutral and unstable atmospheric surface layers over a lake and a vineyard. A decorrelation of the momentum and scalar fluxes is observed with increasing instability. Moreover, different measures of transport efficiency (correlation coefficients, efficiencies based on quadrant analysis and bulk transfer coefficients) indicate that, under close to neutral conditions, momentum and scalars are transported similarly whereas, as the instability of the atmosphere increases, scalars are transported increasingly more efficiently than momentum. This dissimilarity between the turbulent transport of momentum and scalars under unstable conditions concurs with, and is likely caused by, a change in the topology of turbulent coherent structures. Previous laboratory and field studies report that under neutral conditions hairpin vortices and hairpin packets are present and dominate the vertical fluxes, while under free-convection conditions thermal plumes are expected. Our results (cross-stream vorticity variation, quadrant analysis and time series analysis) are in very good agreement with this picture and confirm a change in the structure of the coherent turbulent motions under increasing instability, although the exact structure of these motions and how they are modified by stability requires further investigation based on three-dimensional flow data.     

An Analytical Model for the Two-Scalar Covariance Budget Inside a Uniform Dense Canopy

The two-scalar covariance budget is significant within the canopy sublayer (CSL) given its role in modelling scalar flux budgets using higher-order closure principles and in estimating the segregation ratio for chemically reactive species. Despite its importance, an explicit expression describing how the two-scalar covariance is modified by inhomogeneity in the flow statistics and in the vertical variation in scalar emission or uptake rates within the canopy volume remains elusive even for passive scalars. To progress on a narrower version of this problem, an analytical solution to the two-scalar covariance budget in the CSL is proposed for the most idealized flow conditions: a stationary and planar homogeneous flow inside a uniform and dense canopy with a constant leaf area density distribution. The foliage emission (or uptake) source strengths are assumed to vary exponentially with depth while the forest floor emission is represented as a scalar flux. The analytical solution is a superposition of a homogeneous part that describes how the two-scalar covariance at the canopy top is transported and dissipated within the canopy volume, and an inhomogeneous part governed by local production mechanisms of the two-scalar covariance. The homogeneous part is primarily described by the canopy adjustment length scale, and the attenuation coefficients of the turbulent kinetic energy and the mean velocity. Conditions for which the vertical variation of the two-scalar covariance is controlled by the rapid attenuation in the mean velocity and turbulent kinetic energy profiles, vis-à-vis the vertical variation of the scalar source strength, are explicitly established. This model also demonstrates how dissimilarity in the emissions from the ground, even for the extreme binary case with one scalar turned ‘on’ and the other scalar turned ‘off’, modifies the vertical variation of the two-scalar covariance within the CSL. To assess its applicability to field conditions, the analytical model predictions were compared with observations made at two different forest types—a sparse pine forest at the Hyytiälä SMEAR II-station (in Finland) and a dense alpine hardwood forest at Lavarone (in Italy). While the model assumptions do not represent the precise canopy morphology, attenuation properties of the turbulent kinetic energy and the mean velocity, observed mixing length, and scalar source attenuation properties for these two forest types, good agreement was found between measured and modelled two scalar covariances for multiple scalars and for the triple moments at the Hyytiälä site.     

An Analytical Model for the Distribution of CO2 Sources and Sinks, Fluxes, and Mean Concentration Within the Roughness Sub-Layer

A one-dimensional analytical model that predicts foliage CO₂ uptake rates, turbulent fluxes, and mean concentration throughout the roughness sub-layer (RSL), a layer that extends from the ground surface up to 5h, where h is canopy height, is proposed. The model combines the mean continuity equation for CO₂ with first-order closure principles for turbulent fluxes and simplified physiological and radiative transfer schemes for foliage uptake. This combination results in a second-order ordinary differential equation in which soil respiration (R) and CO₂ concentration well above the RSL are imposed as lower and upper boundary conditions, respectively. An inverse version of the model was tested against datasets from two contrasting ecosystems: a tropical forest (h = 40m) and a managed irrigated rice canopy (h = 0.7m), with good agreement noted between modelled and measured mean CO₂ concentration profiles within the entire RSL. Sensitivity analysis on the model parameters revealed a plausible scaling regime between them and a dimensionless parameter defined by the ratio between external (R) and internal (stomatal conductance) characteristics controlling the CO₂ exchange process. The model can be used to infer the thickness of the RSL for CO₂ exchange, the inequality in zero-plane displacement between CO₂ and momentum, and its consequences on modelled CO₂ fluxes. A simplified version of the solution is well suited for being incorporated into large-scale climate models. Furthermore, the model framework here can be used to a priori estimate relative contributions from the soil surface and the atmosphere to canopy-air CO₂ concentration, thereby making it synergetic to stable isotopes studies.     

Impact of Eddy Characteristics on Turbulent Heat and Momentum Fluxes in the Urban Roughness Sublayer

Eddy-covariance observations above the densely built-up Centre of Nanjing were made from December 2011 to August 2012. Separate eddy-covariance systems installed at two levels on a 36-m tower located on a rooftop were operated simultaneously, and observations grouped into two sectors (A, B) according to the prevalent wind directions. For sector A, where the nearby buildings are all below the lower measurement level, the sensible heat and momentum fluxes are generally greater at the upper level. For sector B, where several high-rise buildings are located upwind, the sensible heat and momentum fluxes at the upper level are close to those at the lower level. The analysis shows that the turbulent eddy characteristics differ between the two wind sectors, leading to a different behaviour of turbulent exchange between the two levels. A hypothesis is proposed that addresses the vertical variation of turbulent fluxes in the urban roughness sublayer (RSL). For sector A, the buildings block the flow, change the trajectory of scalars, and distort the footprint of scalar fluxes; this ‘blocking effect’ is believed to lead to a smaller sensible heat flux above the canopy layer. Such an effect should decrease with height in the RSL, explaining the increase of the observed turbulent heat flux with height. In addition, the presence of non-uniform building heights adversely affects turbulence organization around the canopy top, and likely elevates the inflection point of the mean flow to a higher elevation close to the upper measurement level, where larger shear results in a larger momentum flux. For sector B, wake effects from the nearby high-rise buildings strongly reduce turbulence organization at higher elevations, leading to similar sensible heat and momentum fluxes at both measurement levels.     

Boundary Conditions for Scalar (Co)Variances over Heterogeneous Surfaces

The problem of boundary conditions for the variances and covariances of scalar quantities (e.g., temperature and humidity) at the underlying surface is considered. If the surface is treated as horizontally homogeneous, Monin–Obukhov similarity suggests the Neumann boundary conditions that set the surface fluxes of scalar variances and covariances to zero. Over heterogeneous surfaces, these boundary conditions are not a viable choice since the spatial variability of various surface and soil characteristics, such as the ground fluxes of heat and moisture and the surface radiation balance, is not accounted for. Boundary conditions are developed that are consistent with the tile approach used to compute scalar (and momentum) fluxes over heterogeneous surfaces. To this end, the third-order transport terms (fluxes of variances) are examined analytically using a triple decomposition of fluctuating velocity and scalars into the grid-box mean, the fluctuation of tile-mean quantity about the grid-box mean, and the sub-tile fluctuation. The effect of the proposed boundary conditions on mixing in an archetypical stably-stratified boundary layer is illustrated with a single-column numerical experiment. The proposed boundary conditions should be applied in atmospheric models that utilize turbulence parametrization schemes with transport equations for scalar variances and covariances including the third-order turbulent transport (diffusion) terms.     

A semi-analytical solution for the three-dimensional advection–diffusion equation considering non-local turbulence closure

Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection–diffusion equation. Moreover, large eddies are able to mix scalar quantities in a manner that is counter to the local gradient. In this work we present a semi-analytical solution for the three-dimensional steady-state advection–diffusion equation, considering non-local turbulence closure using the Generalized Integral Advection Diffusion Multilayer Technique (GIADMT). We report some examples of applications of the new solution for two different datasets and for a water tank experiment.     

Latent and sensible heat flux predictions from a uniform pine forest using surface renewal and flux variance methods

A surface renewal model that links organized eddy motion to the latent and sensible heat fluxes is tested with eddy correlation measurements carried out in a 13m tall uniform Loblolly pine plantation in Duke Forest, Durham, North Carolina. The surface renewal model is based on the occurance of ramp-like patterns in the scalar concentration measurements. To extract such ramp-like patterns from Eulerian scalar concentration measurements, a newly proposed time-frequency filtering scheme is developed and tested. The time-domain filtering is carried out using compactly-supported orthonormal wavelets in conjunction with the Universal Wavelet Thresholding approach of Donoho and Johnstone, while the frequency filtering is carried out by a band-pass sine filter centered around the ramp-occurrence frequency as proposed by other studies. The method was separately tested for heat and water vapour with good agreement between eddy correlation flux measurements and model predictions. The usefulness of the flux-variance method to predict sensible and latent heat fluxes is also considered. Our measurements suggest that the simple flux-variance method reproduces the measured heat and momentum fluxes despite the fact that the variances were measured within the roughness sublayer and not in the surface layer. Central to the predictions of water vapour fluxes using the flux-variance approach is the similarity between heat and water vapour transport by the turbulent air flow. This assumption is also investigated for this uniform forest terrain.     

Surface energy balance,parameterizations of boundary-layer heights and the application of resistance laws near an Antarctic Ice Shelf front

A study of the surface energy balance with turbulent fluxes obtained by the Monin-Obukhov similarity theory and a comparison with results for resistance laws are presented for the strong baroclinic conditions in the vicinity of the Filchner/Ronne Ice Shelf front. The data are taken from a field experiment in the Antarctic summer season 1983/84. For the first time in the coastal Antarctic region, this data set comprises synchronous energy balance measurements over the polynya and the ice shelf together with soundings of the boundary layer, yielding vertical profiles of the wind velocity and temperature over the ice shelf, at the ice shelf front and over the polynya.Over the ice shelf, the radiation balance is the largest component of the energy fluxes and is mainly compensated by the subsurface energy flux and the turbulent heat flux in the daily mean. Over the polynya, turbulent fluxes of sensible and latent heat lead to large energy losses of the water surface in the night-time and in situations of very low air temperatures.Different parameterizations for boundary-layer height are compared using tethered sonde and energy balance measurements. With the height of the inversion base over the polynya and the height of the critical bulk Richardson number over the ice shelf, external parameters for the application of resistance laws were determined. The comparison of turbulent surface fluxes obtained by the energy balance measurements and by the resistance laws shows good agreement for the convective conditions over the polynya. For the stably stratified boundary layer over the ice shelf with small amounts of the turbulent heat flux, the deviation is large for the case of a cold air outflow with a superposed inertial oscillation.     

A Lagrangian Stochastic Model for Sea-Spray Evaporation in the Atmospheric Marine Boundary Layer

The dispersion of heavy particles subjected to a turbulent forcing is often simulated with Lagrangian stochastic models. Although these models have been employed successfully over land, the implementation of traditional LS models in the marine boundary layer is significantly more challenging. We present an adaptation of traditional Lagrangian stochastic models to the atmospheric marine boundary layer with a particular focus on the representation of the scalar turbulence for temperature and humidity. In this new model, the atmosphere can be stratified and the bottom boundary is represented by a realistic wavy surface that moves and deforms. Hence, the correlation function for the turbulent flow following a particle is extended to the inhomogenous, anisotropic case. The results reproduce behaviour for scalar Lagrangian turbulence in a stratified airflow that departs only slightly from the expected behaviour in isotropic turbulence. When solving for the surface temperature and the radius of evaporating heavy water droplets in the airflow, the modelled turbulent forcing on the particle also behaves remarkably well. We anticipate that this model will prove especially useful in the context of sea-spray dispersion and its associated sensible heat, latent heat, and gas fluxes between spray droplets and the atmosphere.     

Scaling Properties of Biologically Active Scalar Concentration Fluctuations in the Atmospheric Surface Layer over a Managed Peatland

The higher-order scalar concentration fluctuation properties are examined in the context of Monin–Obukhov similarity theory for a variety of greenhouse gases that have distinct and separate source/sink locations along an otherwise ideal micrometeorological field site. Air temperature and concentrations of water vapour, carbon dioxide and methane were measured at high frequency (10 Hz) above a flat and extensive peat-land soil in the San Joaquin–Sacramento Delta (California, USA) area, subjected to year-round grazing by beef cattle. Because of the heterogeneous distribution of the sources and sinks of CO₂ and especially CH₄ emitted by cattle, the scaling behaviour of the higher-order statistical properties diverged from predictions based on a balance between their production and dissipation rate terms, which can obtained for temperature and H₂O during stationary conditions. We identify and label these departures as ‘exogenous’ because they depend on heterogeneities and non-stationarities induced by boundary conditions on the flow. Spectral analysis revealed that the exogenous effects show their signatures in regions with frequencies lower than those associated with scalar vertical transport by turbulence, though the two regions may partially overlap in some cases. Cospectra of vertical fluxes appear less influenced by these exogenous effects because of the modulating role of the vertical velocity at low frequencies. Finally, under certain conditions, the presence of such exogenous factors in higher-order scalar fluctuation statistics may be ‘fingerprinted’ by a large storage term in the mean scalar budget.     

Application of transilient turbulent theory to study interactions between the atmospheric boundary layer and forest canopies

The new Forest-Land-Atmosphere ModEl called FLAME is presented. The first-order, nonlocal turbulence closure called transilient turbulence theory (Stull, 1993) is applied to study the interactions between a forested land-surface and the atmospheric boundary layer (ABL). The transilient scheme is used for unequal vertical grid spacing and includes the effects of drag, wake turbulence, and interference to vertical mixing by plant elements. Radiation transfer within the vegetation and the equations for the energy balance at the leaf surface have been taken from Norman (1979). Among others, the model predicts profiles of air temperature, humidity and wind velocity within the ABL, sensible and latent heat fluxes from the soil and the vegetation, the stomata and aerodynamic resistances, as well as profiles of temperature and water content in the soil. Preliminary studies carried out for a cloud free day and idealized initial conditions are presented. The canopy height is 30 m within a vertical domain of 3 km. The model is able to capture some of the effects usually observed within and above forested areas, including the relative wind speed maximum in the trunk space and the counter gradient-fluxes in the lower part of the plant stand. Of special interest is the determination of the location and magnitude of the turbulent mixing between model layers, which permits one to identify the effects of large eddies transporting momentum and scalar quantities into the canopy. A comparison between model simulations and field measurements will be presented in a future paper.