Towards an improved and more flexible representation of water stress in coupled photosynthesis-stomatal conductance models

Coupled photosynthesis-stomatal conductance (A-g_s) models are commonly used in ecosystem models to represent the exchange rate of CO₂ and H₂O between vegetation and the atmosphere. The ways these models account for water stress differ greatly among modelling schemes. This study provides insight into the impact of contrasting model configurations of water stress on the simulated leaf-level values of net photosynthesis (A), stomatal conductance (g_s), the functional relationship among them and their ratio, the intrinsic water use efficiency (A/g_s), as soil dries. A simple, yet versatile, normalized soil moisture dependent function was used to account for the effects of water stress on g_s, on mesophyll conductance (g_m) and on the biochemical capacity. Model output was compared to leaf-level values obtained from the literature. The sensitivity analyses emphasized the necessity to combine both stomatal and non-stomatal limitations of A in coupled A-g_s models to accurately capture the observed functional relationships A vs. g_s and A/g_svs. g_s in response to drought. Accounting for water stress in coupled A-g_s models by imposing either stomatal or biochemical limitations of A, as commonly practiced in most ecosystem models, failed to reproduce the observed functional relationship between key leaf gas exchange attributes. A quantitative limitation analysis revealed that the general pattern of C₃ photosynthetic response to water stress may be well represented in coupled A-g_s models by imposing the highest limitation strength to g_m, then to g_s and finally to the biochemical capacity.     

Water use efficiency in a soybean field: influence of plant water stress

Measurements were made in 1980 over a fully-developed soybean (Glycine max (L.) Merrill) canopy at Mead, Nebraska to determine how crop water status influences photosynthesis, evapotranspiration and water use efficiency. Water use efficiency was calculated in terms of the CO₂—water flux ratio (CWFR). Micrometeorological techniques were used to measure the exchange rates of CO₂ and water vapor above the crop canopy. Crop water status was evaluated by reference to volumetric soil moisture (θ_v), stomatal resistance (r_s), and leaf water potential (ψ) measurements.Stomatal resistance (r_s) was independent of ψ when the latter was greater than ?1.1 MPa. r_s increased sharply as ψ dropped below this threshold. Canopy CO₂ exchange (F_c) decreased logarithmically with increasing r_s under strong irradiance. Although F_c was found to be strongly correlated with r_s, the influence of low values of ψ and of high air temperature cannot be discounted since these factors affect the enzymatic reactions associated with photosynthesis. Stomatal closure also reduced evapotranspiration and influenced the partitioning of net radiation.Under strong irradiance the CO₂ water flux ratio (CWFR) decreased with increasing stomatal resistance. This observation is at variance with predictions of certain early ‘resistance’ models, but substantiates predictions of some recent models in which leaf energy balance considerations are incorporated.     

Long-term (13 years) measurements of SO2 fluxes over a forest and their control by surface chemistry

Long-term fluxes of sulphur dioxide (SO₂) have been measured over a mixed suburban forest subjected to elevated SO₂ concentrations. The net exchange was shown to be highly dynamic with substantial periods of both upward and downward fluxes observed in excellent conditions for flux measurement. Upward fluxes constituted 30% of selected fluxes and appeared more frequently when the canopy was acidic. Upward fluxes were shown to be due to desorption from a drying surface or when ambient levels declined after periods of increased SO₂ exposure.The long term average SO₂ flux (F) was −59 ng SO₂ m⁻² s⁻¹ for the period 1997-2009 corresponding to an average SO₂ concentration of 12.3 μg SO₂ m⁻³ and a deposition velocity υ_d of 5 mm s⁻¹. The smallest deposition fluxes and υ_d were measured in dry conditions (−42 ng m⁻² s⁻¹ and 3.5 mm s⁻¹, resp.), which represented 57% of all cases. Wet canopies were more efficient sinks for SO₂ and a dew-wetted canopy had a smaller υ_d (6 mm s⁻¹) than a rain-wetted canopy (ca 10 mm s⁻¹). Seasonal variability reflected differences in chemical climate or canopy buffering properties. During the summer half-year when surface acidity was low due to higher NH₃/SO₂ ratios, a higher deposition efficiency (υ_d/υ_dmax) and lower non-stomatal resistance (R_w) were observed compared to winter conditions. Comparisons of R_c for different combinations of canopy wetness and surface acidity categories emphasized the importance of both factors in regulating the non-stomatal sinks of SO₂. Increased surface water acidity gradually led to a lower υ_d/υ_dmax and an increased R_c for all considered canopy wetness categories. The smallest υ_d/υ_dmax ratio and highest R_c were obtained for a dry canopy with high surface acidity. Conversely, a rain-wetted canopy was the most efficient sink for SO₂. The canopy sink strength was further enhanced by high friction velocities (u_*), optimizing the mechanical mixing into the canopy. Long-term trends were strongly coupled to changes in the NH₃/SO₂ ratio, which has clearly enhanced the deposition efficiency of SO₂ in recent years.     

Measuring forest floor CO2 fluxes in a Douglas-fir forest

CO₂ exchange was measured on the forest floor of a coastal temperate Douglas-fir forest located near Campbell River, British Columbia, Canada. Continuous measurements were obtained at six locations using an automated chamber system between April and December, 2000. Fluxes were measured every half hour by circulating chamber headspace air through a sampling manifold assembly and a closed-path infrared gas analyzer. Maximum CO₂ fluxes measured varied by a factor of almost 3 between the chamber locations, while the highest daily average fluxes observed at two chamber locations occasionally reached values near 15 μmol C m⁻² s⁻¹. Generally, fluxes ranged between 2 and 10 μmol C m⁻² s⁻¹ during the measurement period. CO₂ flux from the forest floor was strongly related to soil temperature with the highest correlation found with 5 cm depth temperature. A simple temperature dependent exponential model fit to the nighttime fluxes revealed Q₁₀ values in the normal range of 2–3 during the warmer parts of the year, but values of 4–5 during cooler periods. Moss photosynthesis was negligible in four of the six chambers, while at the other locations, it reduced daytime half-hourly net CO₂ flux by about 25%. Soil moisture had very little effect on forest floor CO₂ flux. Hysteresis in the annual relationship between chamber fluxes and soil temperatures was observed. Net exchange from the six chambers was estimated to be 1920±530 g C m⁻² per year, the higher estimates exceeding measurement of ecosystem respiration using year-round eddy correlation above the canopy at this site. This discrepancy is attributed to the inadequate number of chambers to obtain a reliable estimate of the spatial average soil CO₂ flux at the site and uncertainty in the eddy covariance respiration measurements.     

Estimating transpiration and the sensitivity of carbon uptake to water availability in a subalpine forest using a simple ecosystem process model informed by measured net CO2 and H2O fluxes

Modeling how the role of forests in the carbon cycle will respond to predicted changes in water availability hinges on an understanding of the processes controlling water use in ecosystems. Recent studies in forest ecosystem modeling have employed data-assimilation techniques to generate parameter sets that conform to observations, and predict net ecosystem CO₂ exchange (NEE) and its component processes. Since the carbon and water cycles are linked, there should be additional process information available from ecosystem H₂O exchange. We coupled SIPNET (Simple Photosynthesis EvapoTranspiration), a simplified model of ecosystem function, with a data-assimilation system to estimate parameters leading to model predictions most closely matching the net CO₂ and H₂O fluxes measured by eddy covariance in a high-elevation, subalpine forest ecosystem. When optimized using measurements of CO₂ exchange, the model matched observed NEE (RMSE = 0.49 g C m⁻²) but underestimated transpiration calculated independently from sap flow measurements by a factor of 4. Consequently, the carbon-only optimization was insensitive to imposed changes in water availability. Including eddy flux data from both CO₂ and H₂O exchange to the optimization reduced the model fit to the observed NEE fluxes only slightly (RME = 0.53 g C m⁻²), however this parameterization also reproduced transpiration calculated from independent sap flow measurements (r² = 0.67, slope = 0.6). A significant amount of information can be extracted from simultaneous analysis of CO₂ and H₂O exchange, which improved the accuracy of transpiration estimates from measured evapotranspiration. Conversely, failure to include both CO₂ and H₂O data streams can generate results that mask the responses of ecosystem carbon cycling to variation in the precipitation. In applying the model conditioned on both CO₂ and H₂O fluxes to the subalpine forest at the Niwot Ridge AmeriFlux site, we observed that the onset of transpiration is coincident with warm soil temperatures. However, after snow has covered the ground in the fall, we observed significant inter-annual variability in the fraction of evapotranspiration composed of transpiration; evapotranspiration was dominated by transpiration in years when late fall air temperatures were high enough to maintain photosynthesis, but by sublimation from the surface of the snowpack in years when late fall air temperatures were colder and forest photosynthetic activity had ceased. Data-assimilation techniques and simultaneous measurements of carbon and water exchange can be used to quantify the response of net carbon uptake to changes in water availability by using an ecosystem model where the carbon and water cycles are linked.     

Comparison of the eddy covariance and automated closed chamber methods for evaluating nocturnal CO2 exchange in a Japanese cypress forest

Underestimation of nocturnal CO₂ respiration using the eddy covariance method under calm conditions remains an unsolved problem at many flux observation sites in forests. To evaluate nocturnal CO₂ exchange in a Japanese cypress forest, we observed CO₂ flux above the canopy (F_c), changes in CO₂ storage in the canopy (S_t) and soil, and trunk and foliar respiration for 2 years (2003–2004). We scaled these chamber data to the soil, trunk, and foliar respiration per unit of ground area (F_s, F_t, F_f, respectively) and used the relationships of F_s, F_t, and F_f with air or soil temperature for comparison with canopy-scale CO₂ exchange measurements (=F_c + S_t). The annual average F_s, F_t, and F_f were 714 g C m⁻² year⁻¹, 170 g C m⁻² year⁻¹, and 575 g C m⁻² year⁻¹, respectively. At small friction velocity (u_*), nocturnal F_c + S_t was smaller than F_s + F_t + F_f estimated using the chamber method, whereas the two values were almost the same at large u_*. We replaced F_c + S_t measured during calm nocturnal periods with a value simulated using a temperature response function derived during well-mixed nocturnal periods. With this correction, the estimated net ecosystem exchange (NEE) from F_c + S_t data ranged from −713 g C m⁻² year⁻¹ to −412 g C m⁻² year⁻¹ in 2003 and from −883 g C m⁻² year⁻¹ to −603 g C m⁻² year⁻¹ in 2004, depending on the u_* threshold. When we replaced all nocturnal F_c + S_t data with F_s + F_t + F_f estimated using the chamber method, NEE was −506 g C m⁻² year⁻¹ and −682 g C m⁻² year⁻¹ for 2003 and 2004, respectively.     

Relationship between wheat yield and foliage temperature: theory and its application to infrared measurements

A theoretical basis is presented for the relationship between crop yield and one time-of-day measurements of the foliage-ambient temperature differential (T_f — T_a). The theory was used to analyse two contrasting relationships between wheat yield and T_f — T_a. The relationships resulted from a range of irrigation treatments and two times of sowing imposed on 26 plots. This caused yields to vary from 8.3 to 1.7 t ha^?1 when the crop was sown in June and from 5.5 to 2.2 t ha^?1 when sown in August. To explain these variations T_f — T_a and associated micrometeorological data were collected around solar noon during the period from jointing to maturity. From these data transpiration and the associated aerodynamic and canopy stomatal resistances to water vapour transport were predicted. The associated canopy conductances for diffusion of CO₂ were derived and used to predict the corresponding CO₂ assimilation rates.The predicted transpiration and CO₂ assimilation rates were closely related to yield within each year but not between years. However if the rates were normalised for the shorter growing season of the late sown crop the yields from the 26 plots formed a common relationship. The transpiration vs. yield relationship was further improved by normalising for differences in foliage vapour pressure deficit. The good agreement between field data and theory was probably due to the dominating effect of stomatal control on both T_f — T_a and CO₂ assimilation rate. If CO₂ assimilation rate is strongly influenced by factors other than soil water stress then the theory may not hold and a different relationship may exist. It was concluded that infrared thermometry is a useful technique for studying yield variations in agronomic experiments where these variations are due to stomatal control.     

Vertical structure of evapotranspiration at a forest site (a case study)

The components of ecosystem evapotranspiration of a Norway spruce forest (Picea abies L.) as well as the vertical structure of canopy evapotranspiration were analyzed with a combination of measurements and models for a case study of 5 days in September 2007. Eddy-covariance and sap flux measurements were performed at several heights within the canopy at the FLUXNET site Waldstein-Weidenbrunnen (DE-Bay) in the Fichtelgebirge mountains in Germany. Within and above canopy fluxes were simulated with two stand-scale models, the 1D multilayer model ACASA that includes a third-order turbulence closure and the 3D model STANDFLUX. The soil and understory evapotranspiration captured with the eddy-covariance system in the trunk space constituted 10% of ecosystem evapotranspiration measured with the eddy-covariance system above the canopy. A comparison of transpiration measured with the sap flux technique and inferred from below and above canopy eddy-covariance systems revealed higher estimates from eddy-covariance measurements than for sap flux measurements. The relative influences of possible sources of this mismatch, such as the assumption of negligible contribution of evaporation from intercepted water, and differences between the eddy-covariance flux footprint and the area used for scaling sap flux measurements, were discussed. Ecosystem evapotranspiration as well as canopy transpiration simulated with the two models captured the dynamics of the measurements well, but slightly underestimated eddy-covariance values. Profile measurements and models also gave us the chance to assess in-canopy profiles of canopy evapotranspiration and the contributions of in-canopy layers. For daytime and a coupled or partly coupled canopy, mean simulated profiles of both models agreed well with eddy-covariance measurements, with a similar performance of the ACASA and the STANDFLUX model. Both models underestimated profiles for nighttime and decoupled conditions. During daytime, the upper half of the canopy contributed approximately 80% to canopy evapotranspiration, whereas during nighttime the contribution shifted to lower parts of the canopy.     

The sensitivity of modeled SO2 fluxes and profiles to stomatal and boundary layer resistances

A higher-order closure model for canopy/surface exchange is presented and applied to an oak-hickory forest to model SO₂ deposition for typical summertime conditions. The model is then used as a tool to investigate the sensitivity of modeled fluxes (the deposition velocity) and concentration profiles to the parameters used to compute the leaf boundary layer and stomatal resistances. Both the deposition velocity and concentration profiles show little sensitivity to variations in the leaf boundary-layer resistance (r _b) since it generally comprises only a small fraction of the total resistance to diffusion from the air to the sub-stomatal cavity. The deposition velocity (V _d ) is more sensitive to variations in the minimum stomatal resistance (r _s ,min) and light response coefficient (β) than r _b .It was found that 50% variations in β give a maximum difference of 0.2 cm s^?1 in V _d while 30% variations in r _s , min produce a maximum difference of near 0.3 cm s^?1.     

Applying a Lagrangian dispersion analysis to infer carbon dioxide and latent heat fluxes in a corn canopy

Warland and Thurtell (2000) proposed an analytical dispersion Lagrangian analysis (hereafter WT analysis) to relate the mean scalar concentration field to source profiles inside the canopy. The first objective of this study was to evaluate the performance of the WT analysis with existing turbulence statistics parameterizations in a corn canopy, by comparing its inferred net ecosystem CO₂ exchange (NEE) and latent heat flux (λE) with eddy covariance measurements. The second objective was to assess the performance of the WT analysis to infer the soil CO₂ flux. Four parameterizations of turbulence statistics were used to estimate Lagrangian time scale (T_L) and standard deviation of vertical wind velocity (σ_w) profiles. The estimated T_L and σ_w profiles were then corrected for atmospheric stability conditions. The field experiment was carried out in a corn field from August to October 2007 and 2008. Profiles of water vapour and CO₂ mixing ratios were measured using a multiport sampling system connected to an infrared gas analyzer. Wind velocity within and above the canopy and eddy covariance measurements over the canopy were taken. The soil respiration, estimated using the WT analysis, was compared to estimates obtained by an empirical model. WT analysis fluxes showed good correlation (R² = 0.77-0.88) with NEE and λE obtained by the eddy covariance technique, but overestimated net fluxes, especially when corrections for atmospheric stability were applied. The optimization of T_L and σ_w profiles using in-canopy turbulence measurements improved the agreement between measured and modeled NEE and λE. Inferred soil CO₂ fluxes were underestimated and were poorly correlated (R² = 0.02-0.01) with estimates obtained using an empirical model based on soil temperature. This poor performance in estimating the soil respiration is likely caused by the decoupling between inside and above canopy flows.     

Forest–atmosphere carbon dioxide exchange in eastern Siberia

We investigated the daily exchange of CO₂ between undisturbed Larix gmelinii (Rupr.) Rupr. forest and the atmosphere at a remote Siberian site during July and August of 1993. Our goal was to measure and partition total CO₂ exchanges into aboveground and belowground components by measuring forest and understory eddy and storage fluxes and then to determine the relationships between the environmental factors and these observations of ecosystem metabolism. Maximum net CO₂ uptake of the forest ecosystem was extremely low compared to the forests elsewhere, reaching a peak of only ∼5 μmol m⁻² s⁻¹ late in the morning. Net ecosystem CO₂ uptake increased with increasing photosynthetically active photon flux density (PPFD) and decreased as the atmospheric water vapor saturation deficit (D) increased. Daytime ecosystem CO₂ uptake increased immediately after rain and declined sharply after about six days of drought. Ecosystem respiration at night averaged ∼2.4 μmol m⁻² s⁻¹ with about 40% of this coming from the forest floor (roots and heterotrophs). The relationship between the understory eddy flux and soil temperature at 5 cm followed an Arrhenius model, increasing exponentially with temperature (Q₁₀∼2.3) so that on hot summer afternoons the ecosystem became a source of CO₂. Tree canopy CO₂ exchange was calculated as the difference between above and below canopy eddy flux. Canopy uptake saturated at ∼6 μmol CO₂ m⁻² s⁻¹ for a PPFD above 500 μmol m⁻² s⁻¹ and decreased with increasing D. The optimal stomatal control model of Mäkelä et al. (1996) was used as a `big leaf' canopy model with parameter values determined by the non-linear least squares. The model accurately simulated the response of the forest to light, saturation deficit and drought. The precision of the model was such that the daily pattern of residuals between modeled and measured forest exchange reproduced the component storage flux. The model and independent leaf-level measurements suggest that the marginal water cost of plant C gain in Larix gmelinii is more similar to values from deciduous or desert species than other boreal forests. During the middle of the summer, the L. gmelinii forest ecosystem is generally a net sink for CO₂, storing ∼0.75 g C m⁻² d⁻¹.     

Seasonal variation and fire effects on CH4, N2O and CO2 exchange in savanna soils of northern Australia

Tropical savanna ecosystems are a major contributor to global CO₂, CH₄ and N₂O greenhouse gas exchange. Savanna fire events represent large, discrete C emissions but the importance of ongoing soil-atmosphere gas exchange is less well understood. Seasonal rainfall and fire events are likely to impact upon savanna soil microbial processes involved in N₂O and CH₄ exchange. We measured soil CO₂, CH₄ and N₂O fluxes in savanna woodland (Eucalyptus tetrodonta/Eucalyptus miniata trees above sorghum grass) at Howard Springs, Australia over a 16 month period from October 2007 to January 2009 using manual chambers and a field-based gas chromatograph connected to automated chambers. The effect of fire on soil gas exchange was investigated through two controlled burns and protected unburnt areas. Fire is a frequent natural and management action in these savanna (every 1-2 years). There was no seasonal change and no fire effect upon soil N₂O exchange. Soil N₂O fluxes were very low, generally between −1.0 and 1.0 μg N m⁻² h⁻¹, and often below the minimum detection limit. There was an increase in soil NH₄⁺ in the months after the 2008 fire event, but no change in soil NO₃⁻. There was considerable nitrification in the early wet season but minimal nitrification at all other times.Savanna soil was generally a net CH₄ sink that equated to between −2.0 and −1.6 kg CH₄ ha⁻¹ y⁻¹ with no clear seasonal pattern in response to changing soil moisture conditions. Irrigation in the dry season significantly reduced soil gas diffusion and as a consequence soil CH₄ uptake. There were short periods of soil CH₄ emission, up to 20 μg C m⁻² h⁻¹, likely to have been caused by termite activity in, or beneath, automated chambers. Soil CO₂ fluxes showed a strong bimodal seasonal pattern, increasing fivefold from the dry into the wet season. Soil moisture showed a weak relationship with soil CH₄ fluxes, but a much stronger relationship with soil CO₂ fluxes, explaining up to 70% of the variation in unburnt treatments. Australian savanna soils are a small N₂O source, and possibly even a sink. Annual soil CH₄ flux measurements suggest that the 1.9 million km² of Australian savanna soils may provide a C sink of between −7.7 and −9.4 Tg CO₂-e per year. This sink estimate would offset potentially 10% of Australian transport related CO₂-e emissions. This CH₄ sink estimate does not include concurrent CH₄ emissions from termite mounds or ephemeral wetlands in Australian savannas.     

Regional surface flux of CO2 inferred from changes in the advected CO2 column density

A Lagrangian experiment was conducted over Iowa during the daytime (9:00–17:30 LT) on June 19, 2007 as part of the North American Carbon Program's Mid-Continent Intensive using a light-weight and operationally flexible aircraft to measure a net drawdown of CO₂ concentration within the boundary layer. The drawdown can be related to net ecosystem exchange when anthropogenic emissions are estimated using a combination of the Vulcan fossil fuel emissions inventory coupled with a source contribution analysis using HYSPLIT. Results show a temporally and spatially averaged net CO₂ flux of −9.0 ± 2.4 μmol m⁻² s⁻¹ measured from the aircraft data. The average flux from anthropogenic emissions over the measurement area was 0.3 ± 0.1 μmol CO₂ m⁻² s⁻¹. Large-scale subsidence occurred during the experiment, entraining 1.0 ± 0.2 μmol CO₂ m⁻² s⁻¹ into the boundary layer. Thus, the CO₂ flux attributable to the vegetation and soils is −10.3 ± 2.4 μmol m⁻² s⁻¹. The magnitude of the calculated daytime biospheric flux is consistent with tower-based eddy covariance fluxes over corn and soybeans given existing land-use estimates for this agricultural region. Flux values are relatively insensitive to the choice of integration height above the boundary layer and emission footprint area. Flux uncertainties are relatively small compared to the biospheric fluxes, though the measurements were conducted at the height of the growing season.     

Use of a greenhouse as an open chamber for canopy gas exchange measurements: Methodology and validation

Measurements of whole-canopy gas exchange - of CO₂ and H₂O - are important for agricultural and ecological reasons. The objective of this study was to investigate the use of a full-size greenhouse as an open-chamber system for measuring canopy-scale gas exchange. Measurements were validated by comparison with gas exchange scaled up from leaf- and plant-level measurements. Leaf-level measurements used a conventional hand-held cuvette gas exchange system at many points in the greenhouse. The experiments were done in a greenhouse with an area of (15 × 24) m² in which a pepper crop was grown. Within the canopy photosynthetic activity and transpiration changed with height, as expected. In addition, it was shown both theoretically and experimentally that in the absence of air mixing within the chamber, gradients of CO₂ and H₂O developed along the airflow direction. Theoretical estimates of the gradients were in good agreement with measured values. In spite of the gradients, canopy photosynthesis and transpiration could be estimated relatively accurately. For instance, the values of canopy photosynthesis and transpiration, during the course of the day, as measured with the open-chamber approach, were in good agreement with mean values obtained from measurements on individual leaves. However, transpiration values obtained both from open-chamber measurements and from individual leaves were generally a little lower than those obtained with lysimeters.     

Measured deuterium in water vapour concentration does not improve the constraint on the partitioning of evapotranspiration in a tall forest canopy, as estimated using a soil vegetation atmosphere transfer model

Partitioning the evapotanspiration (ET) flux in a forest into its component fluxes is important for understanding the water and carbon budgets of the ecosystem. We use non-linear parameter estimation to determine the vertical profile of the Lagrangian timescale (T_L) and partitioning of ET that simultaneously optimise agreement between modelled and measured vertical profiles of temperature, water vapour, carbon dioxide concentrations, and deuterated water vapour for a two-week period in November 2006. High precision real-time trace gas measurements were obtained by FTIR spectroscopy. Modelled temperature and concentration profiles are generated using a Lagrangian dispersion theory combined with source/sink distributions of HDO, H₂O, sensible heat, and CO₂. These distributions are derived from an isotopically enabled multilayer Soil Vegetation Atmospheric Transfer (SVAT) model subject to multiple constraints. The soil component of the model was tested in isolation using measured deuterium content of soil chamber evaporate, while the leaf component was tested using isotopic analyses of leaf and xylem water, combined with leaf-level gas exchange measurements. Optimisation of T_L and the partition of ET was performed twice: once using only temperature, H₂O and CO₂ profiles and a second time including HDO as well. The modelled vertical concentration profiles resulting from inclusion of HDO in the cost function demonstrate our ability to make consistent estimates of both the scalar source distributions and the deuterium content of the water vapour sources. However, introducing measurements of deuterium in water vapour does not significantly alter resulting estimates of normalised T_L (0.4 ± 0.1 at canopy top) and the partition of ET (85 ± 2% transpiration), suggesting that the additional data and modelling required to use deuterium are not warranted for the purpose of partitioning ET using the framework presented here.     

A model to simulate response of plant stomata to environmental conditions

In order to simulate plant transpiration under different field climatic conditions we have developed and verified a semi-empirical model for predicting the stomatal response to solar global radiation, leaf temperature, vapour pressure difference between the leaf and ambient air, ambient air CO₂ concentration and soil water potential. The transpiration and the stomatal relative conductance of a Nicotania Tabaccum var “samsun” plant leaves were measured in a laboratory apparatus and compared to those predicted by the model: good agreement was obtained between them for the different investigated cases. The model was incorporated in a numerical greenhouse microclimate model and its effects on the canopy microclimate are discussed here.     

Carbon dioxide fluxes in corn-soybean rotation in the midwestern U.S.: Inter- and intra-annual variations, and biophysical controls

Quantifying carbon dioxide (CO₂) fluxes in terrestrial ecosystems is critical for better understanding of global carbon cycling and observed changes in climate. This study examined year-round temporal variations of CO₂ fluxes in two biennial crop rotations during 4 year of corn (Zea mays L.) and soybean [Glycine max (L.) Merr.] production. We monitored CO₂ fluxes using eddy-covariance (EC) and soil chambers in adjacent production fields near Ames, Iowa. Under the non-limiting soil water availability conditions predominant in these fields, diel and seasonal variations of CO₂ fluxes were mostly controlled by ambient temperature and available light. Air temperature explained up to 81% of the variability of soil respiratory losses during fallow periods. In contrast, with full-developed canopies, available light was the main driver of daytime CO₂ uptake for both crops. Furthermore, a combined additive effect of both available light and temperature on enhanced CO₂ uptake was identified only for corn. Moreover, diurnal hysteresis of net CO₂ uptake with available light was also found for both crops with consistently greater CO₂ uptake in the mornings than afternoons perhaps primarily owing to delay in peak of soil respiration relative to the time of maximum plant photosynthesis. Annual cumulative CO₂ exchange was mainly determined by crop species with consistently greater net uptake for corn and near neutral exchange for soybean (−466 ± 38 and −13 ± 39 g C m⁻² year⁻¹). Concomitantly, within growing seasons, CO₂ sink periods were approximately 106 days for corn and 90 days for soybean, and peak rates of CO₂ uptake were roughly 1.7-fold higher for corn than soybean. Apparent changes in soil organic carbon estimated after accounting for grain carbon removal suggested soil carbon depletion following soybean years and neutral carbon balance for corn. Overall, results suggest changes in land use and cropping systems have a substantial impact on dynamics of CO₂ exchange.     

Seasonal changes in the spatial structures of N2O, CO2, and CH4 fluxes from Acacia mangium plantation soils in Indonesia

We evaluated the spatial structures of nitrous oxide (N₂O), carbon dioxide (CO₂), and methane (CH₄) fluxes in an Acacia mangium plantation stand in Sumatra, Indonesia, in drier (August) and wetter (March) seasons. A 60 × 100-m plot was established in an A. mangium plantation that included different topographical elements of the upper plateau, lower plateau, upper slope and foot slope. The plot was divided into 10 × 10-m grids and gas fluxes and soil properties were measured at 77 grid points at 10-m intervals within the plot. Spatial structures of the gas fluxes and soil properties were identified using geostatistical analyses. Averaged N₂O and CO₂ fluxes in the wetter season (1.85 mg N m⁻² d⁻¹ and 4.29 g C m⁻² d⁻¹, respectively) were significantly higher than those in the drier season (0.55 mg N m⁻² d⁻¹ and 2.73 g C m⁻² d⁻¹, respectively) and averaged CH₄ uptake rates in the drier season (−0.62 mg C m⁻² d⁻¹) were higher than those in the wetter season (−0.24 mg C m⁻² d⁻¹). These values of N₂O fluxes in A. mangium soils were higher than those reported for natural forest soils in Sumatra, while CO₂ and CH₄ fluxes were in the range of fluxes reported for natural forest soils. Seasonal differences in these gas fluxes appears to be controlled by soil water content and substrate availability due to differing precipitation and mineralization of litter between seasons. N₂O fluxes had strong spatial dependence with a range of about 18 m in both the drier and wetter seasons. Topography was associated with the N₂O fluxes in the wetter season with higher and lower fluxes on the foot slope and on the upper plateau, respectively, via controlling the anaerobic-aerobic conditions in the soils. In the drier season, however, we could not find obvious topographic influences on the spatial patterns of N₂O fluxes and they may have depended on litter amount distribution. CO₂ fluxes had no spatial dependence in both seasons, but the topographic influence was significant in the drier season with lowest fluxes on the foot slope, while there was no significant difference between topographic positions in the wetter season. The distributions of litter amount and soil organic matter were possibly associated with CO₂ fluxes through their effects on microbial activities and fine root distribution in this A. mangium plantation.     

Emissions of CO2, CH4 and N2O from Southern European peatlands

Peatlands play an important role in emissions of the greenhouse gases CO₂, CH₄ and N₂O, which are produced during mineralization of the peat organic matter. To examine the influence of soil type (fen, bog soil) and environmental factors (temperature, groundwater level), emission of CO₂, CH₄ and N₂O and soil temperature and groundwater level were measured weekly or biweekly in loco over a one-year period at four sites located in Ljubljana Marsh, Slovenia using the static chamber technique. The study involved two fen and two bog soils differing in organic carbon and nitrogen content, pH, bulk density, water holding capacity and groundwater level. The lowest CO₂ fluxes occurred during the winter, fluxes of N₂O were highest during summer and early spring (February, March) and fluxes of CH₄ were highest during autumn. The temporal variation in CO₂ fluxes could be explained by seasonal temperature variations, whereas CH₄ and N₂O fluxes could be correlated to groundwater level and soil carbon content. The experimental sites were net sources of measured greenhouse gases except for the drained bog site, which was a net sink of CH₄. The mean fluxes of CO₂ ranged between 139 mg m⁻² h⁻¹ in the undrained bog and 206 mg m⁻² h⁻¹ in the drained fen; mean fluxes of CH₄ were between −0.04 mg m⁻² h⁻¹ in the drained bog and 0.05 mg m⁻² h⁻¹ in the drained fen; and mean fluxes of N₂O were between 0.43 mg m⁻² h⁻¹ in the drained fen and 1.03 mg m⁻² h⁻¹ in the drained bog. These results indicate that the examined peatlands emit similar amounts of CO₂ and CH₄ to peatlands in Central and Northern Europe and significantly higher amounts of N₂O.