Comparison of eddy-covariance measurements of CO2 fluxes by open- and closed-path CO2 analysers

Eddy fluxes of CO₂ estimated using a sonic anemometer and a closed-path analyser were, on average, 16% lower than those obtained with the same anemometer and an adjacent open-path CO₂ analyser. Covariances between vertical windspeed and CO₂ density from the closed-path analyser were calculated using data points for CO₂ that were delayed relative to anemometer data by the time required for a parcel of air to travel from the tube inlet to the CO₂ sensor. Air flow in the intake tube was laminar. Densities of CO₂ that had been corrected for spurious fluctuations arising from fluctuations in temperature and humidity were used in the flux calculations. Corrections for the cross-sensitivity of CO₂ analysers to water vapour were also incorporated. Spectral analysis of the corrected CO₂ signal from the closed-path analyser showed that damping of fluctuations in the sampling tube at frequencies f > 0.1 Hz caused the apparent loss in flux. The measured losses can be predicted accurately using theory that describes the damping of oscillations in a sampling tube. High-frequency response of the closed-path system can be improved substantially by ensuring turbulent flow in the tube, using a combination of high volumetric flow rate and small tube diameter. The analysis of attenuation of turbulent fluctuations in flow through tubes is applicable to the measurement of fluxes of other minor atmospheric constituents using the eddy covariance method.     

Comparison of flux measurements with open- and closed-path gas analyzers above an agricultural field and a forest floor

Comparison was made of the flux measurements of a closed-path CO₂/H₂O analyzer and an open-path H₂O analyzer above a clover field and the forest floor of a Douglas-fir stand. The attenuation of the gas concentration fluctuations caused by the sampling tube of the closed-path analyzer resulted in underestimation of the H20 flux above both surfaces. The degree of underestimation above the clover field depended on wind speed, but was smaller than that calculated from the transfer function for laminar flow in a circular tube and the scalar cospectrum in the neutral and unstable surface layer. Above the forest floor CO₂ fluctuations led those of H₂O by 0.7s. The implications of this are discussed regarding the determination of the time delay caused by the sampling tube of the closedpath analyzer. The day-time CO₂ efflux from the forest floor, averaged over three days, was 0.043 mg/(m²s).     

Eddy correlation measurement of CO2 flux using a closed-path sensor: Theory and field tests against an open-path sensor

We have examined the potential of using a closed-path sensor to accurately measure eddy fluxes of CO₂. Five inlet tubeflow configurations were employed in the experimental setup. The fluxes of CO₂ were compared against those measured with an open-path sensor. Sampling air through an intake tube causes a loss of flux, due to the attenuation of CO₂ density fluctuations. Adjustments need to be made to correct for this loss and to account for density effects due to the simultaneous transfer of heat and water vapor. Theory quantifying these effects is discussed.The raw CO₂ flux measured with the closed-path sensor was smaller than that measured with the open-path sensor by about 15% (on average) for the turbulent tubeflow configurations with a short (3 m) intake tube, by 31% for turbulent tubeflow with a longer (6 m) intake tube and by 24% for laminar tubeflow. The difference was, in part, caused by tube attenuation of the CO₂ density fluctuations and inadequate sensor time response. The elimination of the flux adjustment for the simultaneous transfer of sensible heat (i.e., the attenuation of ambient temperature fluctuations in the intake tube) generally accounted for the rest of this difference.The raw flux measured with the closed-path sensor was corrected for frequency response and density effects. Except in the case of laminar tubeflow, the corrected closed-path flux agreed consistently with the corrected open-path flux within a few percent (<5%). These results suggest that closed-path sensors, with appropriate corrections, can be used to measure CO₂ flux accurately. Recommendations are included on selecting an optimum flow configuration to minimize the effect of sampling air through a tube.Published as Paper No. 9938, Journal Series, Nebraska Agricultural Research Division.     

An evaluation of aircraft flux measurements of CO2, water vapor and sensible heat

Ground-based flux measurements of carbon dioxide and water vapor integrate physiological processes taking place on a field scale. Aircraft flux measurements have recently been undertaken to attempt to widen the scope of applicability of such measurements. However, because of the intermittency of turbulent transfer, flux measurements must be averaged over long periods of time or long distances to give reproducible results. This requirement makes it difficult to relate aircraft flux measurements to local surface processes. Flux measurements of CO₂, latent and sensible heat obtained from repeated passes in four directions and at three elevations over a homogeneous wheat-growing area are compared with ground-based measurements. Averages based on four runs of 4 km in length gave results consistent with ground-based measurements. The largest percentage differences were in the sensible heat flux. Cospectral analyses showed no significant high frequency losses for the data from flight levels of 25 and 50 m, but an underestimation of approximately 10% resulted at 10 m. Flight direction with respect to wind direction was relatively unimportant at 10 and 25 m but some effects were observed at 50 m. It was also shown that at 25 m, over a relatively smooth and homogeneous surface, the means of either three or four runs 4 km in length were similar to the means of 12–16 km runs. This confirms that at this altitude, most of the flux contribution is contained at wavelengths less than 4 km and that the mean of 3 to 4 passes accounts for most of the intermittency of turbulent transfer.     

Mean flux estimation and cospectra of airborne carbon dioxide and water vapour eddy flux measurements in the planetary surface layer

CO₂ and H₂O eddy-flux measurements derived from airborne vertical wind and concentration data are examined to determine the accuracy of the eddy-flux estimate in moderately unstable conditions within the surface layer. Integral scale estimation is used to determine the minimum length of a sample required to achieve a given accuracy. Cospectral analysis is used to examine the scales responsible for the transport of water vapour and carbon dioxide in the surface layer. Cospectra are found to broaden as the sampling altitude increases, a well known result. At low altitudes, the cospectrum is found to become negative for large scales. A significant fraction of the flux is carried by this range of scales, which suggests that scales larger than the boundary-layer height may play an important role in the transport at low altitudes.     

Accuracy of the relaxed eddy-accumulation technique,evaluated using CO2 flux measurements

A system capable of measuring the fluxes of trace gases was developed. It is based on a simpler version of the eddy-accumulation technique (EA), known as the relaxed eddy-accumulation technique (REA). It accumulates air samples associated with updrafts and downdrafts at a constant flow rate in two containers for later analysis of the trace gas mean concentration. The flux integration is based on the durations of updraft and downdraft events, rather than on the vertical wind velocity (W) as is the case for EA and eddy-correlation (EC) techniques. The flux, calculated by the REA technique, is equal to the difference in the mean concentration of the trace gas of interest between the upward and downward moving eddies, multiplied by the standard deviation of the vertical wind velocity and an empirical coefficient. CO₂ fluxes measured for 162 half-hour periods over a soybean field by both EC and REA techniques showed excellent agreement (coefficient of determination,R ²=0.92). The slope (0.985) and the intercept (–0.042 mg m^–2 s^–1) were not significantly different from 1 and 0, respectively, at the 5% level; and the standard error of estimate was 0.074 mg m^–2 s^–1. It is also shown that the empirical coefficient can be calculated from either latent or sensible heat fluxes. A model describing the effect on this empirical coefficient of not sampling aroundW equal to zero is proposed.Centre for Land and Biological Resources Research Contribution No. 92-212.     

Glacial termination: sensitivity to orbital and CO2 forcing in a coupled climate system model

 To study glacial termination and related feedback mechanisms, a continental ice dynamics model is globally and asynchronously coupled to a physical climate (atmosphere-ocean-sea ice) model. The model performs well under present-day, 11 kaBP (thousand years before present) and 21 kaBP perpetual forcing. To address the ice-sheet response under the effects of both perpetual orbital and CO₂ forcing, sensitivity experiments are conducted with two different orbital configurations (11 kaBP and 21 kaBP) and two different atmospheric CO₂ concentrations (200 ppmv and 280 ppmv). This study reveals that, although both orbital and CO₂ forcing have an impact on ice-sheet maintenance and deglacial processes, and although neither acting alone is sufficient to lead to complete deglaciation, orbital forcing seems to be more important. The CO₂ forcing has a large impact on climate, not uniformly or zonally over the globe, but concentrated over the continents adjacent to the North Atlantic. The effect of increased CO₂ (from 200 ppmv to 280 ppmv) on surface air temperature has its peak there in winter associated with a reduction in sea-ice extent in the northern North Atlantic. These changes are accompanied by an enhancement in the intensity of the meridional overturning and poleward ocean heat transport in the North Atlantic. On the other hand, the effect of orbital forcing (from 21 kaBP to 11 kaBP) has its peak in summer. Since the summer temperature, rather than winter temperature, is found to be dominant for the ice-sheet mass balance, orbital forcing has a larger effect than CO₂ forcing in deglaciation. Warm winter sea surface temperature arising from increased CO₂ during the deglaciation contributes to ice-sheet nourishment (negative feedback for ice-sheet retreat) through slightly enhanced precipitation. However, the precipitation effect is totally overwhelmed by the temperature effect. Our results suggest that the last deglaciation was initiated through increasing summer insolation with CO₂ providing a powerful feedback. Received: 22 February 2000 / Accepted: 17 September 2000     

Carbon and water flux responses to physiology by environment interactions: a sensitivity analysis of variation in climate on photosynthetic and stomatal parameters

Sensitivity of carbon uptake and water use estimates to changes in physiology was determined with a coupled photosynthesis and stomatal conductance (g _s) model, linked to canopy microclimate with a spatially explicit scheme (MAESTRA). The sensitivity analyses were conducted over the range of intraspecific physiology parameter variation observed for Acer rubrum L. and temperate hardwood C₃ (C₃) vegetation across the following climate conditions: carbon dioxide concentration 200–700 ppm, photosynthetically active radiation 50–2,000 μmol m^?2 s^?1, air temperature 5–40 °C, relative humidity 5–95 %, and wind speed at the top of the canopy 1–10 m s^?1. Five key physiological inputs [quantum yield of electron transport (α), minimum stomatal conductance (g ₀), stomatal sensitivity to the marginal water cost of carbon gain (g ₁), maximum rate of electron transport (J _max), and maximum carboxylation rate of Rubisco (V _cmax)] changed carbon and water flux estimates ≥15 % in response to climate gradients; variation in α, J _max, and V _cmax input resulted in up to ~50 and 82 % intraspecific and C₃ photosynthesis estimate output differences respectively. Transpiration estimates were affected up to ~46 and 147 % by differences in intraspecific and C₃ g ₁ and g ₀ values—two parameters previously overlooked in modeling land–atmosphere carbon and water exchange. We show that a variable environment, within a canopy or along a climate gradient, changes the spatial parameter effects of g ₀, g ₁, α, J _max, and V _cmax in photosynthesis-g _s models. Since variation in physiology parameter input effects are dependent on climate, this approach can be used to assess the geographical importance of key physiology model inputs when estimating large scale carbon and water exchange.     

On the scavenging of SO2 by cloud and raindrops: I. A theoretical study of SO2 absorption and desorption for water drops in air

An extention of our previous theory for trace gas absorption into freely-falling cloud and raindrops is presented. This theory describes the convective diffusion of a trace gas through air and into a water drop with internal circulation, the drop falling at its terminal velocity. Using flow fields for the circulating water inside and for the moving air outside the drop, obtained by numerical solutions to the Navier—Stokes equation of motion, we numerically solved the convective diffusion equation to determine the uptake of SO₂ by water drops of various sizes, time exposure to the gas phase, and concentration of SO₂ in the gas phase. It was found that for drops of radius larger than 1 mm and relatively low gas concentrations (10 ppb (v)), resistance to gas diffusion lies mainly in the gas phase; while for drops of radius less than 500 m and gas concentrations larger than those found in the atmosphere (1% (v)), the resistance to diffusion lies primarily in the liquid phase. With drop sizes and gas concentrations between these limits, the rate of SO₂ uptake is controlled by a coupled resistance to diffusion inside and outside the drop. In addition to our general model, a simplified version was formulated which allows considerable savings in computer time for evaluation and improved ease of handling without significant loss of accuracy. A comparison between our simplified model and that of Barrie (1978) shows that the boundary-layer approach of Barrie may be a useful alternate approach to estimating trace gas absorption by water drops, provided appropriate values are chosen for the thickness of the boundary layers involved.     

On the scavenging of SO2 by cloud and raindrops: II. An experimental study of SO2 absorption and desorption for water drops in air

For the purpose of testing our previously described theory of SO₂ scavenging a laboratory investigation was carried out in the UCLA 33 m long rainshaft. Drops with radii between 250 and 2500 m were allowed to come to terminal velocity, after which they passed through a chamber of variable length filled with various SO₂ concentrations in air. After falling through a gas separating chamber consisting of a fluorocarbon gas the drops were collected and analyzed for their total S content in order to determine the rate of SO ₂ absorption.The SO₂ concentration in air studied ranged between 1 and 60% (v). Such relatively large concentrations were necessary due to the short times the drops were exposed to SO₂ in the present setup. The present experimental results were therefore not used to simulate atmospheric conditions but rather to test our previously derived theory which is applicable to any laboratory or atmospheric condition. Comparison of our studies with the results from our theory applied to our laboratory conditions led to predicted values for the S concentration in the drops which agreed well with those observed if the drops had radii smaller than 500 m. In order to obtain agreement between predicted and observed S concentrations in larger drops, an empirically derived eddy diffusivity for SO₂ in water had to be included in the theory to take into account the effect of turbulent mixing inside such large drops.In a subsequent set of experiments, drops initially saturated with S (IV) were allowed to fall through S-free air to determine the rate of SO ₂ desorption. The results of these studies also agreed well with the results of our theoretical model, thus justifying the reversibility assumption made in our theoretical models.In a final set of experiments, the effects of oxidation on SO₂ absorption was studied by means of drops containing various amounts of H₂O₂. For comparable exposure times to SO₂, the S concentration in drops with H₂O₂ was found to be up to 10 times higher than the concentration in drops in which no oxidation occurred.     

Synoptic aspects of wet and dry conditions in central Australia: observations and GCM simulations for 1 × CO2 and 2 × CO2 conditions

Synoptic climatological patterns that produce anomalous wet conditions in central Australia during the period from September to April have been studied. The analysis was done by using observed daily rainfall data at a number of stations, wind and mean sea level pressure from the European Centre for Medium Range Weather Forecasts (ECMWF), Tropical Ocean and Global Atmosphere (TOGA) data from 1985 to 1991, and the CSIRO 9-level (CSIR09) global climate model (GCM) simulated data for 1 × CO₂ and 2 × CO₂ experiments. On the basis of rainfall values above 99.5 percentile in observed and simulated data, wet days have been selected to study the synoptic-scale weather systems that produce anomalous wet events in central Australia. As the vast majority of days in central Australia are dry, the same number of days with no rainfall for both observed and simulated conditions have been selected randomly. The observed synoptic climatological patterns have been compared with the results of the control simulation of CSIRO9. A comparison between CSIRO9 simulated synoptic patterns and observed synoptic patterns reveals that the model fairly well captures the synoptic climatological characteristics which produce anomalous wet and contrasting dry weather conditions during the period from September to April. Under enhanced greenhouse experiments, the main features of the synoptic patterns are intensified both for wet and dry conditions, which result in an increase in extreme weather conditions, an increase in rainfall intensity, a spatial expansion of the heavy rainfall region during wet days, and an expansion of the dry area during dry days. During anomalous wet conditions, the low pressure area is intensified, monsoonal winds and southeasterlies are strengthened and strong wind shear over tropical Australia is simulated. During this condition, the monsoon shear line moves poleward particularly over the Northern Territory. In contrast, during dry conditions, the anticyclonic circulation over the continent is strengthened.     

The importance of population,climate change and CO2 plant physiological forcing in determining future global water stress

Future levels of water stress depend on changes in several key factors including population, climate-change driven water availability, and a carbon dioxide physiological-forcing effect on evaporation and run-off. In this study we use an ensemble of the HadCM3 climate model forced with a range of future emissions scenarios combined with a simple water scarcity index to assess the contribution of each of these factors to the projected population living in water stress over the 21st century.Population change only scenarios increase the number of people living in water stress such that at peak global population 65% of people experience some level of water stress. Globally, the climate model ensemble projects an increase in water availability which partially offsets some of the impacts of population growth. The result is 1 billion fewer people living in water stress by the 2080s under the high end emissions scenarios than if population increased in the absence of climate change.This study highlights the important role plant-physiological forcing has on future water resources. The effect of rising CO₂ is to increase available water and to reduce the number of people living in high water stress by around 200 million compared to climate only projections. This effect is of a similar order of magnitude to climate change.     

Atmospheric and terrestrial water budgets: sensitivity and performance of configurations and global driving data for long term continental scale WRF simulations

Driving data and physical parametrizations can significantly impact the performance of regional dynamical atmospheric models in reproducing hydrometeorologically relevant variables. Our study addresses the water budget sensitivity of the Weather Research and Forecasting Model System WRF (WRF-ARW) with respect to two cumulus parametrizations (Kain–Fritsch, Betts–Miller–Janji?), two global driving reanalyses (ECMWF ERA-INTERIM and NCAR/NCEP NNRP), time variant and invariant sea surface temperature and optional gridded nudging. The skill of global and downscaled models is evaluated against different gridded observations for precipitation, 2 m-temperature, evapotranspiration, and against measured discharge time-series on a monthly basis. Multi-year spatial deviation patterns and basin aggregated time series are examined for four globally distributed regions with different climatic characteristics: Siberia, Northern and Western Africa, the Central Australian Plane, and the Amazonian tropics. The simulations cover the period from 2003 to 2006 with a horizontal mesh of 30 km. The results suggest a high sensitivity of the physical parametrizations and the driving data on the water budgets of the regional atmospheric simulations. While the global reanalyses tend to underestimate 2 m-temperature by 0.2–2 K, the regional simulations are typically 0.5–3 K warmer than observed. Many configurations show difficulties in reproducing the water budget terms, e.g. with long-term mean precipitation biases of 150 mm month^?1 and higher. Nevertheless, with the water budget analysis viable setups can be deduced for all four study regions.     

Aircraft measurements and model simulations of stratospheric ozone and N2O: implications for chemistry and transport processes in the models

Airborne measurements of stratospheric ozone and N₂O from the SCIAMACHY (Scanning Imaging Absorption Spectrometer) Validation and Utilization Experiment (SCIA-VALUE) are presented. The campaign was conducted in September 2002 and February–March 2003. The Airborne Submillimeter Radiometer (ASUR) observed stratospheric constituents like O₃ and N₂O, among others, spanning a latitude from 5°S to 80°N during the survey. The tropical ozone source regions show high ozone volume mixing ratios (VMRs) of around 11 ppmv at 33 km altitude, and the altitude of the maximum VMR increases from the tropics to the Arctic. The N₂O VMRs show the largest value of 325 ppbv in the lower stratosphere, indicating their tropospheric origin, and they decrease with increasing altitude and latitude due to photolysis. The sub-tropical and polar mixing barriers are well represented in the N₂O measurements. The most striking seasonal difference found in the measurements is the large polar descent in February–March. The observed features are interpreted with the help of SLIMCAT and Bremen Chemical Transport Model (CTMB) simulations. The SLIMCAT simulations are in good agreement with the measured O₃ and N₂O values, where the differences are within 1 ppmv for O₃ and 15 ppbv for N₂O. However, the CTMB simulations underestimate the tropical middle stratospheric O₃ (1–1.5 ppmv) and the tropical lower stratospheric N₂O (15–30 ppbv) measurements. A detailed analysis with various measurements and model simulations suggests that the biases in the CTMB simulations are related to its parameterised chemistry schemes.