The Influence of Advection on the Short Term CO 2 -Budget in and Above a Forest Canopy

An experimental micrometeorological set-up was established at the CARBOEURO-FLUX site in Tharandt, Germany, to measure all relevant variables for the calculation of the vertical and horizontal advective fluxes of carbon dioxide. The set-up includes two auxiliary towers to measure horizontal and vertical CO₂ and H₂O gradients through the canopy, and to make ultrasonic wind measurements in the trunk space. In combination with the long-term flux tower an approximately even-sided prism with a typical side-length of 50 m was established. It is shown that under stable (nighttime) conditions the mean advective fluxes have magnitudes on the same order as the daily eddy covariance (EC) flux, which implies that they play a significant, but not yet fully understood, role in the carbon budget equation. The two advective fluxes are opposite and seem to cancel each other at night (at least for these measurements). During the day, vertical advection tends to zero, while horizontal advection is still present implying a flow of CO₂ out of the control volume. From our measurements, a mean daily gain of 2.2 gC m^–2 d^–1 for the horizontal advection and a mean daily loss of 2.5 gC m^–2d^–1 for the vertical advection is calculated for a period of 20 days. However the large scatter of the advective fluxes has to be further investigated. It is not clear yet whether the large variability is natural or due to measurement errors and conceptual deficiencies of the experiment. Similar results are found in the few comparable studies.     

Estimating the components of the sensible heat budget of a tall forest canopy in complex terrain

Ultrasonic wind measurements, sonic temperature and air temperature data at two heights in the advection experiment MORE II were used to establish a complete budget of sensible heat including vertical advection, horizontal advection and horizontal turbulent flux divergence. MORE II took place at the long-term Carbo-Europe IP site in Tharandt, Germany. During the growing period of 2003 three additional towers were established to measure all relevant parameters for an estimation of advective fluxes, primarily of CO₂. Additionally, in relation to other advection experiments, a calculation of the horizontal turbulent flux divergence is proposed and the relation of this flux to atmospheric stability and friction velocity is discussed. In order to obtain a complete budget, different scaling heights for horizontal advection and horizontal turbulent flux divergence are tested. It is shown that neglecting advective fluxes may lead to incorrect results. If advective fluxes are taken into account, the sensible heat budget based upon vertical turbulent flux and storage change only, is reduced by approximately 30%. Additional consideration of horizontal turbulent flux divergence would in turn add 5–10% to this sum (i.e., the sum of vertical turbulent flux plus storage change plus horizontal and vertical advection). In comparison with available energy horizontal advection is important at night whilst horizontal turbulent flux divergence is rather insignificant. Obviously, advective fluxes typically improve poor nighttime energy budget closure and might change ecosystem respiration fluxes considerably.     

Comparing CO2 Storage and Advection Conditions at Night at Different Carboeuroflux Sites

Anemometer and CO₂ concentration data from temporary campaigns performed at six CARBOEUROFLUX forest sites were used to estimate the importance of non-turbulent fluxes in nighttime conditions. While storage was observed to be significant only during periods of both low turbulence and low advection, the advective fluxes strongly influence the nocturnal CO₂ balance, with the exception of almost flat and highly homogeneous sites. On the basis of the main factors determining the onset of advective fluxes, the ‘advection velocity’, which takes net radiation and local topography into account, was introduced as a criterion to characterise the conditions of storage enrichment/depletion. Comparative analyses of the six sites showed several common features of the advective fluxes but also some substantial differences. In particular, all sites where advection occurs show the onset of a boundary layer characterised by a downslope flow, negative vertical velocities and negative vertical CO₂ concentration gradients during nighttime. As a consequence, vertical advection was observed to be positive at all sites, which corresponds to a removal of CO₂ from the ecosystem. The main differences between sites are the distance from the ridge, which influences the boundary-layer depth, and the sign of the mean horizontal CO₂ concentration gradients, which is probably determined by the source/sink distribution. As a consequence, both positive and negative horizontal advective fluxes (corresponding respectively to CO₂ removal from the ecosystem and to CO₂ supply to the ecosystem) were observed. Conclusive results on the importance of non-turbulent components in the mass balance require, however, further experimental investigations at sites with different topographies, slopes, different land covers, which would allow a more comprehensive analysis of the processes underlying the occurrence of advective fluxes. The quantification of these processes would help to better quantify nocturnal CO₂ exchange rates.     

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.     

Horizontal and Vertical Turbulent Fluxes Forced by a Gravity Wave Event in the Nocturnal Atmospheric Surface Layer Over the Amazon Forest

A nocturnal gravity wave was detected over a south-western Amazon forest during the Large-Scale Biosphere–Atmosphere experiment in Amazonia (LBA) in the course of the dry-to-wet season campaign on October 2002. The atmospheric surface layer was stably stratified and had low turbulence activity, based on friction velocity values. However, the passage of the wave, an event with a period of about 180–300 s, caused negative turbulent fluxes of carbon dioxide (CO₂) and positive sensible heat fluxes, as measured by the eddy-covariance system at 60 m (≈30 m above the tree tops). The evolution of vertical profiles of air temperature, specific humidity and wind speed during the wave movement revealed that cold and drier air occupied the sub-canopy space while high wind speeds were measured above the vegetation. The analysis of wind speed and scalars high frequency data was performed using the wavelet technique, which enables the decomposition of signals in several frequencies allowed by the data sampling conditions. The results showed that the time series of vertical velocity and air temperature were −90° out of phase during the passage of the wave, implying no direct vertical transport of heat. Similarly, the time series of vertical velocity and CO₂ concentration were 90° out of phase. The wave was not directly associated with vertical fluxes of this variable but the mixing induced by its passage resulted in significant exchanges in smaller scales as measured by the eddy-covariance system. The phase differences between horizontal velocity and both air temperature and CO₂ concentration were, respectively, zero and 180°, implying phase and anti-phase relationships. As a result, the wave contributed to positive horizontal fluxes of heat and negative horizontal fluxes of carbon dioxide. Such results have to be considered in nocturnal boundary-layer surface-atmosphere exchange schemes for modelling purposes.     

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.     

Aspects of CO2 advection measurements

Observations of vegetation–atmosphere exchange of carbon dioxide (CO₂) by the eddy covariance (EC) technique are limited by difficult conditions such as nighttime and heterogeneous terrain. Thus, advective flux components are included into the net ecosystem exchange (NEE) budget. However, advection measurements are experimentally challenging and do not always help to solve the night flux problem of the EC technique. This study investigates alternative methods for the observation of horizontal advection, in particular horizontal concentration gradients, as well as different approaches to coordinate rotation and vertical advection. Continuous high-frequency measurements of the horizontal CO₂ concentration field are employed and compared to the often used discontinuous sequential sampling. Significant differences were found in the case of 30-min mean concentration values between the conventional discontinuous sampling approach and the complete observation of the time series by continuous sampling. Estimates of vertical advection rely on accurate estimates of vertical wind velocity ( $\emph{w}$ ). Therefore, different approaches to the planar fit coordinate rotation have been investigated. Sector-wise rotation was able to eliminate directional dependencies of mean $\emph{w}$ . Furthermore, the effect of the data set length used for rotation (window length) was investigated and was found to have significant impact on estimates of vertical advection, with larger window lengths yielding about 50% larger vertical advection. A sequential planar fit with controlled window length is proposed to give reproducible results. The different approaches to the measurement and calculation of horizontal and vertical advection presented are applied to data obtained during the exchange processes in mountainous region experiment at the FLUXNET site Waldstein–Weidenbrunnen (DE-Bay). Estimates of NEE including advection are compared to NEE from turbulent and storage flux alone without advection. NEE including vertical advection with sector-wise planar fit rotation and controlled window length and including horizontal advection from continuous gradient measurements, which were comprehensively bias corrected by a new approach, did compare well with the expected night flux error, with meteorological drivers of the fluxes and with soil chamber measurements. Unrealistically large and noisy values of horizontal advection from the conventional discontinuous sampling approach, which lead to unrealistic values of NEE, could be eliminated by the alternative approaches presented. We therefore suggest the further testing of those approaches at other sites in order to improve the accuracy of advection measurements and, subsequently, estimates of NEE.     

Seasonal and diurnal variations of coherent structures over a deciduous forest

Coherent structures in turbulent flow above a midlatitude deciduous forest are identified using a wavelet analysis technique. Coupling between motions above the canopy (z/h=1.5, whereh is canopy height) and within the canopy (z/h=0.6) are studied using composite velocity and temperature fields constructed from 85 hours of data. Data are classified into winter and summer cases, for both convective and stable conditions. Vertical velocity fluctuations are in phase at both observation levels. Horizontal motions associated with the structures within the canopy lead those above the canopy, and linear analysis indicates that the horizontal motions deep in the canopy should lead the vertical motions by 90°. On average, coherent structures are responsible for only about 40% of overall turbulent heat and momentum fluxes, much less than previously reported. However, our large data set reveals that this flux fraction comes from a wide distribution that includes much higher fractions in its upper extremes. The separation distanceL _s between adjacent coherent structures, 6–10h, is comparable to that obtained in previous observations over short canopies and in the laboratory. Changes in separation between the summer and winter (leafless) conditions are consistent withL _s being determined by a local horizontal wind shear scale.     

The scaling of flow in vegetative structures

Sensible heat flux data obtained from the U.S. Environmental Protection Agency's Regional Air Pollution Study in St. Louis, Missouri are presented and discussed. Large spatial variations exist in heat flux on both a land-use scale and the urban scale. Arguments based upon empirical data and sampling theory show that estimates of heat flux representative of an upwind fetch l _x require a minimum observation height proportional to z ^3/4. The influence of advection on the magnitude of the heat flux is also explored for both the urban and sub-urban or land-use scales. The data clearly indicate that advection can strongly modulate and even dominate the vertical heat flux above surfaces in areas which maintain large horizontal temperature gradients. The advection contribution is positive for cold air advection and negative for warm air advection, and may result from either the urban heat island or land-use mesoscale features. The depth of advective influence is directly proportional to the horizontal scale of the phenomenon and inversely proportional to horizontal temperature gradient.     

Some methodological questions concerning advection measurements: a case study

A dataset from two campaigns conducted at the Vielsalm experimental site in Belgium was used as a basis for discussing some methodological problems and providing intermediate results on estimating CO₂ advection. The analysis focused on the horizontal [CO₂] gradient and on the vertical velocity w, the variables most affected by uncertainty. The sampling error for half-hourly horizontal [CO₂] gradients was estimated to be 1.3 μmol mol⁻¹. Despite this important random error for half-hour estimations of [CO₂], the mean horizontal [CO₂] gradients in advective conditions were shown to be representative at the ecosystem scale and to extend only to the lowest part of a drainage sub-layer, which developed in the trunk space. By contrast, under daytime conditions, this gradient was shown to be more sensitive to local source heterogeneities. The estimation of the short-term averaged vertical velocity ( was the greater source of error when computing advection terms. The traditional correction methods used to obtain are discussed and a (co)sine correction is tested to highlight the instrumental origin of the offset in w. A comparison of measurements by sonic anemometers placed close together above the canopy showed that the uncertainty on was 0.042 m s⁻¹, which is of the same order of magnitude as the velocity itself. In addition, as the drainage sub-layer is limited to the lowest part of the canopy, the representativeness of is questionable. An alternative computation using the divergence of the horizontal wind speed in the trunk space produced a estimation that was four times lower than the single-point measurement. However, this value gives a more realistic estimate of the vertical advection term and improves the CO₂ budget closure at the site.     

An Alternative Approach for CO2 Flux Correction Caused by Heat and Water Vapour Transfer

Energy and CO₂ fluxes are commonly measured above plant canopies using an eddy covariance system that consists of a three-dimensional sonic anemometer and an H₂O/CO₂ infrared gas analyzer. By assuming that the dry air is conserved and inducing mean vertical velocity, Webb et al. (Quart. J. Roy. Meteorol. Soc. 106, 85-100, 1980) obtained two equations to account for density effects due to heat and water vapour transfer on H₂O/CO₂ fluxes. In this paper, directly starting with physical consideration of air-parcel expansion/compression, we derive two alternative equations to correct for these effects that do not require the assumption that dry air is conserved and the use of the mean vertical velocity. We then applied these equations to eddy flux observations from a black spruce forest in interior Alaska during the summer of 2002. In this ecosystem, the equations developed here led to increased estimates of CO₂ uptake by the vegetation during the day (up to about 20%), and decreased estimates of CO₂ respiration by the ecosystem during the night (approximately 4%) as compared with estimates obtained using the Webb et al. approach.     

Turbulent Flux Measurements Above and Below the Overstory of a Boreal Aspen Forest

Turbulent flux measurements both above and beneath the canopy of a boreal aspen forest are described. Velocity skewness showed that, beneath the aspen canopy, turbulence was dominated by intermittent, downward penetrating gusts. Eulerian horizontal length scales calculated from integration of the autocorrelation function or spectral peaks were 9.0 and 1.4 times the mean aspen height of 21.5 m respectively. Above-canopy power spectral slopes for all velocity components followed the -2/3 power law, whereas beneath-canopy slopes were closer to -1 and showed a spectral short cut in the horizontal and vertical components. Cospectral patterns were similar both above and beneath the canopy. The Monin–Obukhov similarity function for the vertical wind velocity variance was a well-defined function of atmospheric stability, both above and beneath the canopy. Nocturnal flux underestimation and departures of this similarity function from that expected from Monin–Obukhov theory were a function of friction velocity. Energy balance closure greater than 80% was achieved at friction velocities greater than 0.30 and 0.10 m s-1, above and below the aspen canopy, respectively. Recalculating the latent heat flux using various averaging periods revealed a minimum of 15 min were required to capture 90% of the 30-min flux. Linear detrending reduced the flux at shorter averaging periods compared to block averaging. Lack of energy balance closure and erratic flux behaviour led to the recalculation of the latent and sensible heat fluxes using the ratio of net radiation to the sum of the energy balance terms.     

Estimating Co2 Source/Sink Distributions Within A Rice Canopy Using Higher-Order Closure Model

Source/sink strengths and vertical fluxdistributions of carbon dioxide within and above arice canopy were modelled using measured meanconcentration profiles collected during aninternational rice experiment in Okayama, Japan (IREX96). The model utilizes an Eulerian higher-orderclosure approach that permits coupling of scalar andmomentum transport within vegetation to infer sourcesand sinks from mean scalar concentration profiles; theso-called `inverse problem'. To compute the requiredvelocity statistics, a Eulerian second-order closuremodel was considered. The model well reproducedmeasured first and second moment velocity statisticsinside the canopy. Using these modelled velocitystatistics, scalar fluxes within and above the canopywere computed and compared with CO<sub>2</sub>eddy-correlation measurements above the canopy. Goodagreement was obtained between model calculations offluxes at the top of the canopy and measurements. Close to the ground, the model predicted higherrespiratory fluxes when the paddy was drained comparedto when it was flooded. This is consistent with thefloodwater providing a barrier to diffusion ofCO<sub>2</sub> from the soil to the atmosphere. TheEulerian sources and flux calculations were alsocompared to source and flux distributions estimatedindependently using a Lagrangian Localized Near Fieldtheory, the first study to make such a comparison.Some differences in source distributions werepredicted by these analyses. Despite this, thecalculated fluxes by the two approaches compared wellprovided a closure constant, accounting for theinfluence of `near-field' sources in the Eulerian fluxtransport term, was given a value of 1.5 instead ofthe value of 8 found in laboratory studies.     

The Influence of Hilly Terrain on Canopy-Atmosphere Carbon Dioxide Exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO₂ (N_ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO₂ sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO₂ exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO₂. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N_ee = F_c, the vertical turbulent flux of CO₂. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F_c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F_c-at canopy top differs by more than 20% compared to the flat-terrain case.     

Intra-field variability of scalar flux densities across a transition between a desert and an irrigated potato field

This paper reports on measurements of sensible and latent heat and CO₂ fluxes made over an irrigated potato field, growing next to a patch of desert. The study was conducted using two eddy correlation systems. One measurement system was located within the equilibrium boundary layer 800 m downwind from the edge of the potato field. The other measurement system was mobile and was placed at various downwind positions to probe the horizontal transition of vertical scalar fluxes. Latent (LE) and sensible (H) heat fluxes, measured at 4 m above the surface, exhibited marked variations with downwind distance over the field. Only after the fetch to height ratio exceeded 75 to 1 didLE andH become invariant with downwind distance. When latent and sensible heat fluxes were measured upwind of this threshold, significant advection of humidity-deficit occurred, causing a vertical flux divergence ofH andLE.The measured fluxes of momentum, heat, and moisture were compared with predictions from a second-order closure two-dimensional atmospheric boundary layer model. There is good agreement between measurements and model predictions. A soil-plant-atmosphere model was used to examine nonlinear feedbacks between humidity-deficits, stomatal conductance and evaporation. Data interpretation with this model revealed that the advection of hot dry air did not enhance surface evaporation rates near the upwind edge of the potato field, because of negative feedbacks among stomatal conductance, humidity-deficits, andLE. This finding is consistent with results from several recent studies.     

On the flow structure and turbulence during sweep and ejection events in a wind-tunnel model canopy

Particle image velocimetry (PIV) data obtained in a wind-tunnel model of a canopy boundary layer is used to examine the characteristics of mean flow and turbulence. The vector spacing varies between 1.7 and 2.5 times the Kolmogorov scales. Conditional sampling based on quadrants, i.e. based on the signs of velocity fluctuations, reveals fundamental differences in flow structure, especially between sweep and ejection events, which dominate the flow. During sweeps, the downward flow generates a narrow, highly turbulent, shear layer containing multiple small-scale vortices just below canopy height. During ejections, the upward flow expands this shear layer and the associated small-scale flow structures to a broad region located above the canopy. Consequently, during sweeps the turbulent kinetic energy (TKE), Reynolds stresses, as well as production and dissipation rates, have distinct narrow peaks just below canopy height, whereas during ejections these variables have broad maxima well above the canopy. Three methods to estimate the dissipation rate are compared, including spectral fits, measured subgrid-scale (SGS) energy fluxes at different scales, and direct measurements of slightly underresolved instantaneous velocity gradients. The SGS energy flux is 40–60% of the gradient-based (direct) estimates for filter sizes inside the inertial range, while decreasing with scale, as expected, within the dissipation range. The spectral fits are within 5–30% of the direct estimates. The spectral fits exceed the direct estimates near canopy height, but are lower well above and below canopy height. The dissipation rate below canopy height increases with velocity magnitude, i.e. it has the highest values during sweep and quadrant 1 events, and is significantly lower during ejection and quadrant 3 events. Well above the canopy, ejections are the most dissipative. Turbulent transport during sweep events acts as a source below the narrow shear layer within the canopy and as a sink above it. Transport during ejection events is a source only well above the canopy. The residual term in the TKE transport equation, representing mostly the effect of pressure–velocity correlations, is substantial only within the canopy, and is dominated by sweeps.     

Heat and moisture budgets from airborne measurements and high-resolution model simulations

High-resolution simulations with a mesoscale model are performed to estimate heat and moisture budgets of a well-mixed boundary layer. The model budgets are validated against energy budgets obtained from airborne measurements over heterogeneous terrain in Western Germany. Time rate of change, vertical divergence, and horizontal advection for an atmospheric column of air are estimated. Results show that the time trend of specific humidity exhibits some deficiencies, while the potential temperature trend is matched accurately. Furthermore, the simulated turbulent surface fluxes of sensible and latent heat are comparable to the measured fluxes, leading to similar values of the vertical divergence. The analysis of different horizontal model resolutions exhibits improved surface fluxes with increased resolution, a fact attributed to a reduced aggregation effect. Scale-interaction effects could be identified: while time trends and advection are strongly influenced by mesoscale forcing, the turbulent surface fluxes are mainly controlled by microscale processes.     

Limitations of an eddy-correlation technique for the determination of the carbon dioxide and sensible heat fluxes

A modified infrared CO₂ gas analyzer, a small thermocouple assembly, a heated-thermocouple anemometer for horizontal wind, and a propeller-type vertical wind sensor were used to measure the eddy fluxes of heat and CO₂ above a corn crop. Experimental results of these fluxes are discussed. The main sources of errors of the eddy fluxes using these instruments were estimated:

1. Sensors with a time constant of 0.5 s appear to be fast enough to detect most of the vertical CO₂ transfer as long as the sensors are located at least one meter above the crop surface.
2. The deviation from steady-state conditions for 10-min periods was found to have a significant effect on the eddy flux estimates.
3. Temperature fluctuations of the air sample passing through the CO₂ infrared gas analyzer were found to be non-negligible but could be easily corrected.
4. A 1° misalignment of the vertical anemometer affected these eddy fluxes by less than 10% under all circumstances studied.

     

Forest-Air Fluxes Of Carbon, Water And Energy Over Non-Flat Terrain

A field study of surface-air exchange of carbon, water, and energy was conducted at a mid-latitude, mixed forest on non-flat terrain to investigate how to best interpret biological signals from the eddy flux data that may be subject to advective influences. It is shown that during periods of Southwest winds (sector with mild topography), the eddy fluxes are well-behaved in terms of energy balance closure, the existence of a constant flux layer, consistency with chamber observations and the expected abiotic controls on the fluxes. Advective influences are evident during periods with wind from a steep (15%) slope to the Northeast of the tower. These influences appear more severe on CO₂ flux, particularly in stable air, than on the energy fluxes. Large positive flux of CO₂ (> 23 mol m^-2 s^-1) occurs frequently at night. The annual sum of the carbon flux is positive, but the issue about whether the forest is a source of atmospheric carbon remains inconclusive.Attempts are made to assess vertical advectionusing the data collected on a single tower. Over the Southwestsector, vertical advection makes a statistically significant but small contribution to the 30-min energy imbalance and CO₂ flux variations. Contributions by horizontal advection may be larger but cannot be verified directly by the current experimental method.     

A Case Study of CO<Subscript>2</Subscript>, CO and Particles Content Evolution in the Suburban Atmospheric Boundary Layer Using a 2-μm Doppler DIAL,a 1-μm Backscatter Lidar and an Array of In-situ Sensors

A network of remote and in-situ sensors was deployed in a Paris suburb in order to evaluate the mesoscale evolution of the daily cycle of CO₂ and related tracers in the atmospheric boundary layer (ABL) and its relation to ABL dynamics and nearby natural and anthropogenic sources and sinks. A 2-μm heterodyne Doppler differential absorption lidar, which combines measurements of, (1) structure of the atmosphere, (2) radial velocity, and (3) CO₂ differential absorption was a particularly unique element of the observational array. We analyse the differences in the diurnal cycle of CO, CO₂, lidar reflectivity (a proxy for aerosol content) and H₂O using the lidar, airborne measurements in the free troposphere and ground-based measurements made at two sites located few kilometres apart. We demonstrate that vertical mixing dominates the early morning drawdown of CO and aerosol content trapped in the former nocturnal layer but not the H₂O and CO₂ mixing ratio variations. Surface fluxes, vertical mixing and advection all contribute to the ABL CO₂ mixing ratio decrease during the morning transition, with the relative importance depending on the rate and timing of ABL rise. We also show evidence that when the ABL is stable, small-scale (0.1-km vertical and 1-km horizontal) gradients of CO₂ and CO are large. The results illustrate the complexity of inferring surface fluxes of CO₂ from atmospheric budgets in the stable boundary layer.