Modeling dynamics of stable carbon isotopic exchange between a boreal forest ecosystem and the atmosphere

Stable isotopes of CO₂ contain unique information on the biological and physical processes that exchange CO₂ between terrestrial ecosystems and the atmosphere. In this study, we developed an integrated modeling system to simulate dynamics of stable carbon isotope of CO₂, as well as moisture, energy, and momentum, between a boreal forest ecosystem and the atmosphere, as well as their transport/mixing processes through the convective boundary layer (CBL), using remotely sensed surface parameters to characterize the surface heterogeneity. It has the following characteristics: (i) it accounts for the influences of the CBL turbulent mixing and entrainment of the air aloft; (ii) it scales individual leaf‐level photosynthetic discrimination up to the whole canopy (Δ_canopy) through the separation of sunlit and shaded leaf groups; (iii) it has the capacity to examine the detailed interrelationships among plant water‐use efficiency, isotope discrimination, and vapor pressure deficit; and (iv) it has the potential to investigate how an ecosystem discriminates against ¹³C at various time and spatial scales. The monthly mean isotopic signatures of ecosystem respiration (i.e. δ¹³C_R) used for isotope flux calculation are retrieved from the nighttime flask data from the intensive campaigns (1998–2000) at 20 m level on Fraserdale tower, and the data from the growing season in 1999 are used for model validation. Both the simulated CO₂ mixing ratio and δ¹³C of CO₂ at the 20 m level agreed with the measurements well in different phases of the growing season. On a diurnal basis, the greatest photosynthetic discrimination at canopy level (i.e. Δ_canopy) occurred early morning and late afternoon with a varying range of 10–26‰. The diurnal variability of Δ_canopy was also associated with the phases of growing season and meteorological variables. The annual mean Δ_canopy in 1999 was computed to be 19.58‰. The monthly averages of Δ_canopy varied between 18.55‰ and 20.84‰ with a seasonal peak during the middle growing season. Because of the strong opposing influences of respired and photosynthetic fluxes on forest air (both CO₂ and ¹³CO₂) on both the diurnal and seasonal time scales, CO₂ was consistently enriched with the heavier ¹³C isotope (less negative δ¹³C) from July to October and depleted during the remaining months, whereas on a diurnal basis, CO₂ was enriched with the heavier ¹³C in the late afternoon and depleted in early morning. For the year 1999, the model results reveal that the boreal ecosystem in the vicinity of Fraserdale tower was a small sink with net uptake of 29.07 g ¹²C m^?2 yr^?1 and 0.34 g ¹³C m^?2 yr^?1.     

Partitioning net ecosystem carbon exchange with isotopic fluxes of CO2

Because biological and physical processes alter the stable isotopic composition of atmospheric CO₂, variations in isotopic content can be used to investigate those processes. Isotopic flux measurements of ¹³CO₂ above terrestrial ecosystems can potentially be used to separate net ecosystem CO₂ exchange (NEE) into its component fluxes, net photosynthetic assimilation (F_A) and ecosystem respiration (F_R). In this paper theory is developed to partition measured NEE into F_A and F_R, using measurements of fluxes of CO₂ and ¹³CO₂, and isotopic composition of respired CO₂ and forest air. The theory is then applied to fluxes measured (or estimated, for ¹³CO₂) in a temperate deciduous forest in eastern Tennessee (Walker Branch Watershed). It appears that there is indeed enough additional information in ¹³CO₂ fluxes to partition NEE into its photosynthetic and respiratory components. Diurnal patterns in F_A and F_R were obtained, which are consistent in magnitude and shape with patterns obtained from NEE measurements and an exponential regression between night‐time NEE and temperature (a standard technique which provides alternate estimates of F_R and F_A). The light response curve for photosynthesis (F_A vs. PAR) was weakly nonlinear, indicating potential for saturation at high light intensities. Assimilation‐weighted discrimination against ¹³CO₂ for this forest during July 1999 was 16.8–17.1‰, depending on canopy conductance. The greatest uncertainties in this approach lie in the evaluation of canopy conductance and its effect on whole‐canopy photosynthetic discrimination, and thus the indirect methods used to estimate isotopic fluxes. Direct eddy covariance measurements of ¹³CO₂ flux are needed to assess the validity of the assumptions used and provide defensible isotope‐based estimates of the component fluxes of net ecosystem exchange.     

Canopy-scale δ13C of photosynthetic and respiratory CO2 fluxes: observations in forest biomes across the United States

The δ¹³C values of atmospheric carbon dioxide (CO₂) can be used to partition global patterns of CO₂ source/sink relationships among terrestrial and oceanic ecosystems using the inversion technique. This approach is very sensitive to estimates of photosynthetic ¹³C discrimination by terrestrial vegetation (Δ_A), and depends on δ¹³C values of respired CO₂ fluxes (δ¹³C_R). Here we show that by combining two independent data streams – the stable isotope ratios of atmospheric CO₂ and eddy‐covariance CO₂ flux measurements – canopy scale estimates of Δ_A can be successfully derived in terrestrial ecosystems. We also present the first weekly dataset of seasonal variations in δ¹³C_R from dominant forest ecosystems in the United States between 2001 and 2003. Our observations indicate considerable summer‐time variation in the weekly value of δ¹³C_R within coniferous forests (4.0‰ and 5.4‰ at Wind River Canopy Crane Research Facility and Howland Forest, respectively, between May and September). The monthly mean values of δ¹³C_R showed a smaller range (2–3‰), which appeared to significantly correlate with soil water availability. Values of δ¹³C_R were less variable during the growing season at the deciduous forest (Harvard Forest). We suggest that the negative correlation between δ¹³C_R and soil moisture content observed in the two coniferous forests should represent a general ecosystem response to the changes in the distribution of water resources because of climate change. Shifts in δ¹³C_R and Δ_A could be of sufficient magnitude globally to impact partitioning calculations of CO₂ sinks between oceanic and terrestrial compartments.     

Influence of stand structure on carbon-13 of vegetation, soils, and canopy air within deciduous and evergreen forests in Utah, United States

Carbon isotope ratios (δ¹³C) were studied in evergreen and deciduous forest ecosystems in semi-arid Utah (Pinus contorta, Populus tremuloides, Acer negundo and Acer grandidentatum). Measurements were taken in four to five stands of each forest ecosystem differing in overstory leaf area index (LAI) during two consecutive growing seasons. The δ¹³C_leaf (and carbon isotope discrimination) of understory vegetation in the evergreen stands (LAI 1.5–2.2) did not differ among canopies with increasing LAI, whereas understory in the deciduous stands (LAI 1.5–4.5) exhibited strongly decreasing δ¹³C_leaf values (increasing carbon isotope discrimination) with increasing LAI. The δ¹³C values of needles and leaves at the top of the canopy were relatively constant over the entire LAI range, indicating no change in intrinsic water-use efficiency with overstory LAI. In all canopies, δ¹³C_leaf decreased with decreasing height above the forest floor, primarily due to physiological changes affecting c _i/c _a (> 60%) and to a minor extent due to δ¹³C of canopy air (< 40%). This intra-canopy depletion of δ¹³C_leaf was lowest in the open stand (1‰) and greatest in the denser stands (4.5‰). Although overstory δ¹³C_leaf did not change with canopy LAI, δ¹³C of soil organic carbon increased with increasing LAI in Pinus contorta and Populus tremuloides ecosystems. In addition, δ¹³C of decomposing organic carbon became increasingly enriched over time (by 1.7–2.9‰) for all deciduous and evergreen dry temperate forests. The δ¹³C_canopy of CO₂ in canopy air varied temporally and spatially in all forest stands. Vertical canopy gradients of δ¹³C_canopy, and [CO₂]_canopy were larger in the deciduous Populus tremuloides than in the evergreen Pinu contorta stands of similar LAI. In a very wet and cool year, ecosystem discrimination (Δ_e) was similar for both deciduous Populus tremulodies (18.0 ± 0.7‰) and evergreen Pinus contorta (18.3 ± 0.9‰) stands. Gradients of δ¹³C_canopy and [CO₂]_canopy were larger in denser Acer spp. stands than those in the open stand. However, ¹³C enrichment above and photosynthetic draw-down of [CO₂]_canopy below tropospheric baseline values were larger in the open than in the dense stands, due to the presence of a vigorous understory vegetation. Seasonal patterns of the relationship δ¹³C_canopy versus 1/[CO₂]_canopy were strongly influenced by precipitation and air temperature during the growing season. Estimates of Δ_e for Acer spp. did not show a significant effect of stand structure, and averaged 16.8 ± 0.5‰ in 1933 and 17.4 ± 0.7‰ in 1994. However, Δ_e varied seasonally with small fluctuations for the open stand (2‰), but more pronounced changes for the dense stand (5‰). Received: 15 April 1996 / Accepted: 19 October 1996     

18O composition of CO2 and H2O ecosystem pools and fluxes in a tallgrass prairie: Simulations and comparisons to measurements

In this paper we describe measurements and modeling of ¹⁸O in CO₂ and H₂O pools and fluxes at a tallgrass prairie site in Oklahoma. We present measurements of the δ¹⁸O value of leaf water, depth‐resolved soil water, atmospheric water vapor, and Keeling plot δ¹⁸O intercepts for net soil‐surface CO₂ and ecosystem CO₂ and H₂O fluxes during three periods of the 2000 growing season. Daytime discrimination against C¹⁸OO, as calculated from measured above‐canopy CO₂ and δ¹⁸O gradients, is also presented. To interpret the isotope measurements, we applied an integrated land‐surface and isotope model (ISOLSM) that simulates ecosystem H₂¹⁸O and C¹⁸OO stocks and fluxes. ISOLSM accurately predicted the measured isotopic composition of ecosystem water pools and the δ¹⁸O value of net ecosystem CO₂ and H₂O fluxes. Simulations indicate that incomplete equilibration between CO₂ and H₂O within C₄ plant leaves can have a substantial impact on ecosystem discrimination. Diurnal variations in the δ¹⁸O value of above‐canopy vapor had a small impact on the predicted δ¹⁸O value of ecosystem water pools, although sustained differences had a large impact. Diurnal variations in the δ¹⁸O value of above‐canopy CO₂ substantially affected the predicted ecosystem discrimination. Leaves dominate the ecosystem ¹⁸O‐isoflux in CO₂ during the growing season, while the soil contribution is relatively small and less variable. However, interpreting daytime measurements of ecosystem C¹⁸OO fluxes requires accurate predictions of both soil and leaf ¹⁸O‐isofluxes.     

Seasonal variations of CO2 and water vapour exchange rates over a temperate deciduous forest

Long-term and direct measurements of CO₂ and water vapour exchange are needed over forested ecosystems to determine their net annual fluxes of carbon dioxide and water. Such measurements are also needed to parameterize and test biogeochemical, ecological and hydrological assessment models. Responding to this need, eddy covariance measurements of CO₂ and water vapour were made ever a deciduous forest growing near Oak Ridge, TN, between April 1993 and April 1994. Periodic measurements were made of leaf area index, stomatal resistance, soil moisture and pre-dawn leaf water potential to characterize the gas exchange capacity of the canopy. Four factors had a disproportionate influence on the seasonal variation of CO₂ flux densities. These factors were photon flux densities (during the growing season), temperature (during the dormant season), leaf area index and the occurrence of drought The drought period occurred during the peak of the growing season and caused a significant decline in daily and hourly CO₂ flux densities, relative to observations over the stand when soil moisture was plentiful. The annual net uptake of carbon was calculated by integrating flux measurements and filling missing and spurious data with the relations obtained between measured CO₂ fluxes and environmental forcing variables. The net flux of carbon for the period between April 1993 and April 1994 was -525 g C m^?2 y^?1. This value represents a net flux of carbon from the atmosphere and into the forest. The net annual carbon exchange of this southern temperate broadleaved forest exceeded values measured over a northern temperate forest (which experiences a shorter growing season and has less leaf area) by 200 g C m^?2 y^?1 (cf. Wofsy et al 1993). The seasonal variation of canopy evaporation (latent heat flux) was controlled mostly by changes in leaf area and net radiation. A strong depression in evaporation rates was not observed during the drought Over a broadleaved forest large vapour pressure deficits promote evaporation and trees in a mixed stand are able to tap a variety of deep and shallow water sources.     

Future atmospheric CO2 leads to delayed autumnal senescence

Growing seasons are getting longer, a phenomenon partially explained by increasing global temperatures. Recent reports suggest that a strong correlation exists between warming and advances in spring phenology but that a weaker correlation is evident between warming and autumnal events implying that other factors may be influencing the timing of autumnal phenology. Using freely rooted, field‐grown Populus in two Free Air CO₂ Enrichment Experiments (AspenFACE and PopFACE), we present evidence from two continents and over 2 years that increasing atmospheric CO₂ acts directly to delay autumnal leaf coloration and leaf fall. In an atmosphere enriched in CO₂ (by ～45% of the current atmospheric concentration to 550 ppm) the end of season decline in canopy normalized difference vegetation index (NDVI) – a commonly used global index for vegetation greenness – was significantly delayed, indicating a greener autumnal canopy, relative to that in ambient CO₂. This was supported by a significant delay in the decline of autumnal canopy leaf area index in elevated as compared with ambient CO₂, and a significantly smaller decline in end of season leaf chlorophyll content. Leaf level photosynthetic activity and carbon uptake in elevated CO₂ during the senescence period was also enhanced compared with ambient CO₂. The findings reveal a direct effect of rising atmospheric CO₂, independent of temperature in delaying autumnal senescence for Populus, an important deciduous forest tree with implications for forest productivity and adaptation to a future high CO₂ world.     

δ13C of CO2 respired in the dark in relation to δ13C of leaf metabolites: comparison between Nicotiana sylvestris and Helianthus annuus under drought

The variations of δ¹³C in leaf metabolites (lipids, organic acids, starch and soluble sugars), leaf organic matter and CO₂ respired in the dark from leaves of Nicotiana sylvestris and Helianthus annuus were investigated during a progressive drought. Under well‐watered conditions, CO₂ respired in the dark was ¹³C‐enriched compared to sucrose by about 4‰ in N. sylvestris and by about 3‰ and 6‰ in two different sets of experiments in H. annuus plants. In a previous work on cotyledonary leaves of Phaseolus vulgaris, we observed a constant ¹³C‐enrichment by about 6‰ in respired CO₂ compared to sucrose, suggesting a constant fractionation during dark respiration, whatever the leaf age and relative water content. In contrast, the ¹³C‐enrichment in respired CO₂ increased in dehydrated N. sylvestris and decreased in dehydrated H. annuus in comparison with control plants. We conclude that (i) carbon isotope fractionation during dark respiration is a widespread phenomenon occurring in C₃ plants, but that (ii) this fractionation is not constant and varies among species and (iii) it also varies with environmental conditions (water deficit in the present work) but differently among species. We also conclude that (iv) a discrimination during dark respiration processes occurred, releasing CO₂ enriched in ¹³C compared to several major leaf reserves (carbohydrates, lipids and organic acids) and whole leaf organic matter.     

Pan-European δ13C values of air and organic matter from forest ecosystems

We present carbon stable isotope, δ¹³C, results from air and organic matter samples collected during 98 individual field campaigns across a network of Carboeuroflux forest sites in 2001 (14 sites) and 2002 (16 sites). Using these data, we tested the hypothesis that δ¹³C values derived from large‐scale atmospheric measurements and models, which are routinely used to partition carbon fluxes between land and ocean, and potentially between respiration and photosynthesis on land, are consistent with directly measured ecosystem‐scale δ¹³C values. In this framework, we also tested the potential of δ¹³C in canopy air and plant organic matter to record regional‐scale ecophysiological patterns. Our network estimates for the mean δ¹³C of ecosystem respired CO₂ and the related ‘discrimination’ of ecosystem respiration, δ_er and Δ_er, respectively, were ?25.6±1.9‰ and 17.8 ±2.0‰ in 2001 and ?26.6±1.5‰ and 19.0±1.6‰ in 2002. The results were in close agreement with δ¹³C values derived from regional‐scale atmospheric measurement programs for 2001, but less so in 2002, which had an unusual precipitation pattern. This suggests that regional‐scale atmospheric sampling programs generally capture ecosystem δ¹³C signals over Europe, but may be limited in capturing some of the interannual variations. In 2001, but less so in 2002, there were discernable longitudinal and seasonal trends in δ_er. From west to east, across the network, there was a general enrichment in ¹³C (～3‰ and ～1‰ for the 2 years, respectively) consistent with increasing Gorczynski continentality index for warmer and drier conditions. In 2001 only, seasonal ¹³C enrichment between July and September, followed by depletion in November (from about ?26.0‰ to ?24.5‰ to ?30.0‰), was also observed. In 2001, July and August δ_er values across the network were significantly related to average daytime vapor pressure deficit (VPD), relative humidity (RH), and, to a lesser degree, air temperature (T_a), but not significantly with monthly average precipitation (P_m). In contrast, in 2002 (a much wetter peak season), δ_er was significantly related with T_a, but not significantly with VPD and RH. The important role of plant physiological processes on δ_er in 2001 was emphasized by a relatively rapid turnover (between 1 and 6 days) of assimilated carbon inferred from time‐lag analyses of δ_er vs. meteorological parameters. However, this was not evident in 2002. These analyses also noted corresponding diurnal cycles of δ_er and meteorological parameters in 2001, indicating a rapid transmission of daytime meteorology, via physiological responses, to the δ_er signal during this season. Organic matter δ¹³C results showed progressive ¹³C enrichment from leaves, through stems and roots to soil organic matter, which may be explained by ¹³C fractionation during respiration. This enrichment was species dependent and was prominent in angiosperms but not in gymnosperms. δ¹³C values of organic matter of any of the plant components did not well represent short‐term δ_er values during the seasonal cycle, and could not be used to partition ecosystem respiration into autotrophic and heterotrophic components.     

Estimating the contribution of leaf litter decomposition to soil CO2 efflux in a beech forest using 13C-depleted litter

The contribution of leaf litter decomposition to total soil CO₂ efflux (F_L/F) was evaluated in a beech (Fagus sylvatica L.) forest in eastern France. The Keeling‐plot approach was applied to estimate the isotopic composition of respired soil CO₂ from soil covered with either control (?30.32‰) or ¹³C‐depleted leaf litter (?49.96‰). The δ¹³C of respired soil CO₂ ranged from ?25.50‰ to ?22.60‰ and from ?24.95‰ to ?20.77‰, respectively, with depleted or control litter above the soil. The F_L/F ratio was calculated by a single isotope linear mixing model based on mass conservation equations. It showed seasonal variations, increasing from 2.8% in early spring to about 11.4% in mid summer, and decreasing to 4.2% just after leaf fall. Between December 2001 and December 2002, cumulated F and F_L reached 0.98 and 0.08 kg_C m^?2, respectively. On an annual basis, decomposition of fresh leaf litter accounted for 8% of soil respiration and 80% of total C loss from fresh leaf litter. The other fraction of carbon loss during leaf litter decomposition that is assumed to have entered the soil organic matter pool (i.e. 20%) represents only 0.02 kg_C m^?2.     

Diurnal variation of Δ13CO2, ΔC18O16O and evaporative site enrichment of δH218O in Piper aduncum under field conditions in Trinidad

Concurrent measurements of gas exchange, instantaneous isotope discrimination (Δ) against ¹³CO₂ and C¹⁸O¹⁶O, and extent of ¹⁸O enrichment in H₂O at the evaporative sites, were followed in a tropical forest pioneer, Piper aduncum, on two different days in Trinidad during February 1995. Δ¹³CO₂ differed from that predicted from measurements of internal:external CO₂ concentration (C_i/C_a) and showed a wide range of values which decreased throughout the course of the day. Derivation of C_c (the CO₂ concentration at the carboxylation site) was not possible using carbon isotope discrimination under field conditions in situ and was derived assuming a constant value of internal transfer conductance (g_w). Under low rates of assimilation the derived C_c/C_a, like C_i/C_a, remained relatively stable over the course of both days and ΔC¹⁸O¹⁶O followed evaporative demand. Lower values of ΔC¹⁸O¹⁶O on day 2 occurred in response to the indirect effect of increased leaf-to-air vapour pressure deficits (VPD) and reduced stomatal conductance. For the first time, direct determination of the δH₂¹⁸O of transpired water vapour (δ_t) allowed derivation of evaporative site enrichment without the prerequisite of isotopic steady state (ISS) defined in the Craig and Gordon model. Generally, δ_t was less enriched than the source water (δ_s) in the morning and more enriched in the afternoon, which would be predicted from an increase and decrease in ambient VPD, respectively. On both days, leaves of P. aduncum approached ISS (indicated where δ_t≈δ_s) between 1300 and 1500 h. Evaporative site enrichment was maintained into the late afternoon, despite a decrease in ambient VPD. The data presented provide a greater insight into the natural variation in isotopic discrimination under field conditions, which may help to refine models of terrestrial biome discrimination.     

Interseasonal comparison of CO2 concentrations, isotopic composition, and carbon dynamics in an Amazonian rainforest (French Guiana)

Canopy CO₂ concentrations in a tropical rainforest in French Guiana were measured continuously for 5 days during the 1994 dry season and the 1995 wet season. Carbon dioxide concentrations ([CO₂]) throughout the canopy (0.02–38 m) showed a distinct daily pattern, were well-stratified and decreased with increasing height into the canopy. During both seasons, daytime [CO₂] in the upper and middle canopy decreased on average 7–10 μmol mol⁻¹ below tropospheric baseline values measured at Barbados. Within the main part of the canopy (≥ 0.7 m), [CO₂] did not differ between the wet and dry seasons. In contrast, [CO₂] below 0.7 m were generally higher during the dry season, resulting in larger [CO₂] gradients. Supporting this observation, soil CO₂ efflux was on average higher during the dry season than during the wet season, either due to diffusive limitations and/or to oxygen deficiency of root and microbial respiration. Soil respiration rates decreased by 40% after strong rain events, resulting in a rapid decrease in canopy [CO₂] immediately above the forest floor of about 50␣μmol mol⁻¹. Temporal and spatial variations in [CO₂]_canopy were reflected in changes of δ¹³C_canopy and δ¹⁸O_canopy values. Tight relationships were observed between δ¹³C and δ¹⁸O of canopy CO₂ during both seasons (r ² > 0.86). The most depleted δ¹³C_canopy and δ¹⁸O_canopy values were measured immediately above the forest floor (δ¹³C = −16.4‰; δ¹⁸O = 39.1‰ SMOW). Gradients in the isotope ratios of CO₂ between the top of the canopy and the forest floor ranged between 2.0‰ and 6.3‰ for δ¹³C, and between 1.0‰ and 3.5‰ for δ¹⁸O. The δ¹³C_leaf and calculated c _i/c _a of foliage at three different positions were similar for the dry and wet seasons indicating that the canopy maintained a constant ratio of photosynthesis to stomatal conductance. About 20% of the differences in δ¹³C_leaf within the canopy was accounted for by source air effects, the remaining 80% must be due to changes in c _i/c _a. Plotting 1/[CO₂] vs. the corresponding δ¹³C ratios resulted in very tight, linear relationships (r ² = 0.99), with no significant differences between the two seasons, suggesting negligible seasonal variability in turbulent mixing relative to ecosystem gas exchange. The intercepts of these relationships that should be indicative of the δ¹³C of respired sources were close to the measured δ¹³C of soil respired CO₂ and to the δ¹³C of litter and soil organic matter. Estimates of carbon isotope discrimination of the entire ecosystem, Δ_e, were calculated as 20.3‰ during the dry season and as 20.5‰ during the wet season. Received: 3 March 1996 / Accepted: 19 October 1996     

Exposure to an enriched CO2 atmosphere alters carbon assimilation and allocation in a pine forest ecosystem

We linked a leaf-level CO₂ assimilation model with a model that accounts for light attenuation in the canopy and measurements of sap-flux-based canopy conductance into a new canopy conductance-constrained carbon assimilation (4C-A) model. We estimated canopy CO₂ uptake (A_nC) at the Duke Forest free-air CO₂ enrichment (FACE) study. Rates of A_nC estimated from the 4C-A model agreed well with leaf gas exchange measurements (A_net) in both CO₂ treatments. Under ambient conditions, monthly sums of net CO₂ uptake by the canopy (A_nC) were 13% higher than estimates based on eddy-covariance and chamber measurements. Annual estimates of A_nC were only 3% higher than carbon (C) accumulations and losses estimated from ground-based measurements for the entire stand. The C budget for the Pinus taeda component was well constrained (within 1% of ground-based measurements). Although the closure of the C budget for the broadleaf species was poorer (within 20%), these species are a minor component of the forest. Under elevated CO₂, the C used annually for growth, turnover, and respiration balanced only 80% of the A_nC. Of the extra 700 g C m⁻² a⁻¹ (1999 and 2000 average), 86% is attributable to surface soil CO₂ efflux. This suggests that the production and turnover of fine roots was underestimated or that mycorrhizae and rhizodeposition became an increasingly important component of the C balance. Under elevated CO₂, net ecosystem production increased by 272 g C m⁻² a⁻¹: 44% greater than under ambient CO₂. The majority (87%) of this C was sequestered in a moderately long-term C pool in wood, with the remainder in the forest floor–soil subsystem.     

Water savings in mature deciduous forest trees under elevated CO2

Stomatal conductance of plants exposed to elevated CO₂ is often reduced. Whether this leads to water savings in tall forest‐trees under future CO₂ concentrations is largely unknown but could have significant implications for climate and hydrology. We used three different sets of measurements (sap flow, soil moisture and canopy temperature) to quantify potential water savings under elevated CO₂ in a ca. 35 m tall, ca. 100 years old mixed deciduous forest. Part of the forest canopy was exposed to 540 ppm CO₂ during daylight hours using free air CO₂ enrichment (FACE) and the Swiss Canopy Crane (SCC). Across species and a wide range of weather conditions, sap flow was reduced by 14% in trees subjected to elevated CO₂, yielding ca. 10% reduction in evapotranspiration. This signal is likely to diminish as atmospheric feedback through reduced moistening of the air comes into play at landscape scale. Vapour pressure deficit (VPD)‐sap flow response curves show that the CO₂ effect is greatest at low VPD, and that sap flow saturation tends to occur at lower VPD in CO₂‐treated trees. Matching stomatal response data, the CO₂ effect was largely produced by Carpinus and Fagus, with Quercus contributing little. In line with these findings, soil moisture at 10 cm depth decreased at a slower rate under high‐CO₂ trees than under control trees during rainless periods, with a reversal of this trend during prolonged drought when CO₂‐treated trees take advantage from initial water savings. High‐resolution thermal images taken at different heights above the forest canopy did detect reduced water loss through altered energy balance only at <5 m distance (0.44 K leaf warming of CO₂‐treated Fagus trees). Short discontinuations of CO₂ supply during morning hours had no measurable canopy temperature effects, most likely because the stomatal effects were small compared with the aerodynamic constraints in these dense, broad‐leaved canopies. Hence, on a seasonal basis, these data suggest a <10% reduction in water consumption in this type of forest when the atmosphere reaches 540% ppm CO₂.     

Temporal dynamics and spatial variability in the enhancement of canopy leaf area under elevated atmospheric CO2

Increased canopy leaf area (L) may lead to higher forest productivity and alter processes such as species dynamics and ecosystem mass and energy fluxes. Few CO₂ enrichment studies have been conducted in closed canopy forests and none have shown a sustained enhancement of L. We reconstructed 8 years (1996–2003) of L at Duke's Free Air CO₂ Enrichment experiment to determine the effects of elevated atmospheric CO₂ concentration ([CO₂]) on L before and after canopy closure in a pine forest with a hardwood component, focusing on interactions with temporal variation in water availability and spatial variation in nitrogen (N) supply. The dynamics of L were reconstructed using data on leaf litterfall mass and specific leaf area for hardwoods, and needle litterfall mass and specific leaf area combined with needle elongation rates, and fascicle and shoot counts for pines. The dynamics of pine L production and senescence were unaffected by elevated [CO₂], although L senescence for hardwoods was slowed. Elevated [CO₂] enhanced pine L and the total canopy L (combined pine and hardwood species; P<0.050); on average, enhancement following canopy closure was ～16% and 14% respectively. However, variation in pine L and its response to elevated [CO₂] was not random. Each year pine L under ambient and elevated [CO₂] was spatially correlated to the variability in site nitrogen availability (e.g. r²=0.94 and 0.87 in 2001, when L was highest before declining due to droughts and storms), with the [CO₂]‐induced enhancement increasing with N (P=0.061). Incorporating data on N beyond the range of native fertility, achieved through N fertilization, indicated that pine L had reached the site maximum under elevated [CO₂] where native N was highest. Thus closed canopy pine forests may be able to increase leaf area under elevated [CO₂] in moderate fertility sites, but are unable to respond to [CO₂] in both infertile sites (insufficient resources) and sites having high levels of fertility (maximum utilization of resources). The total canopy L, representing the combined L of pine and hardwood species, was constant across the N gradient under both ambient and elevated [CO₂], generating a constant enhancement of canopy L. Thus, in mixed species stands, L of canopy hardwoods which developed on lower fertility sites (～3 g N inputs m^?2 yr^?1) may be sufficiently enhanced under elevated [CO₂] to compensate for the lack of response in pine L, and generate an appreciable response of total canopy L (～14%).     

Measuring and modelling carbon dioxide and water vapour exchange over a temperate broad-leaved forest during the 1995 summer drought

Forests in the south-eastern United States experienced a prolonged dry spell and above-normal temperatures during the 1995 growing season. During this episode, nearly continuous, eddy covariance measurements of carbon dioxide and water vapour fluxes were acquired over a temperate, hardwood forest. These data are used to examine how environmental factors and accumulating soil moisture deficits affected the diurnal pattern and magnitude of canopy-scale carbon dioxide and water vapour fluxes. The field data are also used to test an integrative leaf-to-canopy scaling model (CANOAK), which uses micrometeorological and physiological theory, to calculate mass and energy fluxes. When soil moisture was ample in the spring, peak rates of net ecosystem CO₂ exchange (N_F) occurred around midday and exceeded 20 μmol m^?2 s^?1. Rates of N_K were near optimal when air temperature ranged between 22 and 25°C. The accumulation of soil moisture deficits and a co-occurrence of high temperatures caused peak rates of daytime carbon dioxide uptake to occur earlier in the morning. High air temperatures and soil moisture deficits were also correlated with a dramatic reduction in the magnitude of N_E. On average, the magnitude of N_E decreased from 20 to 7 μmol m^?2 s^?1 as air temperature increased from 24 to 30°C and the soil dried. The CANAOK model yielded accurate estimates of canopy-scale carbon dioxide and water vapour fluxes when the forest had an ample supply of soil moisture. During the drought and heat spell, a cumulative drought index was needed to adjust the proportionality constant of the stomatal conductance model to yield accurate estimates of canopy CO₂ exchange. The adoption of the drought index also enabled the CANOAK model to give improved estimates of evaporation until midday. On the other hand, the scheme failed to yield accurate estimates of evaporation during the afternoon.     

温带混交林碳水通量模拟及其对冠层分层方式的响应——耦合的气孔导度-光合作用-能量平衡模型

利用Leuning建立的耦合的光合作用、气孔导度和能量平衡方程,以将冠层分成多层的方式,包括Gaussian五点积分法、将冠层平均分为多层的方法,逐层计算温带混交林的碳水通量,最后累加至冠层尺度,以模拟CO2和H2O通量。该模型以常规气象观测数据作为驱动变量,计算出冠层与大气之间的碳水交换,与涡动相关系统的通量观测数据进行比较,分析了不同的冠层分层方式对多层模型模拟结果的影响。从3个温带混交林通量站涡动相关系统的能量平衡闭合度来看,中国长白山站CBS、韩国GDK和日本MMF站点的能量平衡比率(EBR)分别为0.76、0.66和1.07,居于国际同类观测范围(0.34—1.2)的中上水平,因此,涡动相关系统的观测数据较为可靠。从碳水通量的日变化来看,用Gaussian五点积分法将冠层分为五层的模型能较好的模拟碳水通量的"单峰形"日变化趋势。夜间Fc为负值且变化趋势较为平缓,表明生态系统进行呼吸作用释放CO2,从日出开始Fc逐渐变为正值,表明生态系统进行光合作用吸收CO2,Fc在中午时分达到最大值,下午Fc逐渐减小,日落之后又回复到夜间的负值。H2O通量的日变化曲线与CO2通量相似,且模拟值与涡动相关实测值具有较好的一致性。在多层模型中,对冠层采用不同的分层方法,对碳水通量模拟结果有一定的影响。以Gaussian五点积分法将冠层分为五层的方法作为对照,分别将冠层平均分为2、5、10、20层的方法得到的碳水通量与其进行比较。从平均值来看,分层越多,H2O通量模拟值越低,而CO2通量模拟值越高。不同的分层方法产生的差异,主要来自于不同层的辐射吸收、温湿度、风速等环境要素的垂直廓线差异,且叶片光合作用对光的响应是非线性的。     

CO2 flux from soil in pastures and forests in southwestern Amazonia

Stocks of carbon in Amazonian forest biomass and soils have received considerable research attention because of their potential as sources and sinks of atmospheric CO₂. Fluxes of CO₂ from soil to the atmosphere, on the other hand, have not been addressed comprehensively in regard to temporal and spatial variations and to land cover change, and have been measured directly only in a few locations in Amazonia. Considerable variation exists across the Amazon Basin in soil properties, climate, and management practices in forests and cattle pastures that might affect soil CO₂ fluxes. Here we report soil CO₂ fluxes from an area of rapid deforestation in the southwestern Amazonian state of Acre. Specifically we addressed (1) the seasonal variation of soil CO₂ fluxes, soil moisture, and soil temperature; (2) the effects of land cover (pastures, mature, and secondary forests) on these fluxes; (3) annual estimates of soil respiration; and (4) the relative contributions of grass‐derived and forest‐derived C as indicated by δ¹³CO₂. Fluxes were greatest during the wet season and declined during the dry season in all land covers. Soil respiration was significantly correlated with soil water‐filled pore space but not correlated with temperature. Annual fluxes were higher in pastures compared with mature and secondary forests, and some of the pastures also had higher soil C stocks. The δ¹³C of CO₂ respired in pasture soils showed that high respiration rates in pastures were derived almost entirely from grass root respiration and decomposition of grass residues. These results indicate that the pastures are very productive and that the larger flux of C cycling through pasture soils compared with forest soils is probably due to greater allocation of C belowground. Secondary forests had soil respiration rates similar to mature forests, and there was no correlation between soil respiration and either forest age or forest biomass. Hence, belowground allocation of C does not appear to be directly related to the stature of vegetation in this region. Variation in seasonal and annual rates of soil respiration of these forests and pastures is more indicative of flux of C through the soil rather than major net changes in ecosystem C stocks.     

Influence of vegetation and soil CO2 exchange on the concentration and stable oxygen isotope ratio of atmospheric CO2 within a Pinus resinosa canopy

Measurements were made of the concentration and stable oxygen isotopic ratio of carbon dioxide in air samples collected on a diurnal basis at two heights within a Pinus resinosa canopy. Large changes in CO₂ concentration and isotopic composition were observed during diurnal time courses on all three symple dates. In addition, there was strong vertical stratification in the forest canopy, with higher CO₂ concentrations and more negative ¹⁸O values observed closer to the soil surface. The observed daily increases in ¹⁸O values of forest CO₂ were dependent on relative humidity consistent with the modelled predictions of isotopic fractionation during photosynthetic gas exchange. During photosynthetic gas exchange, a portion of the CO₂ that enters the leaf and equilibrates with leaf water is not fixed and diffuses back out of the leaf with an altered oxygen isotopic ratio. The oxygen isotope ratio of CO₂ diffusing out of a leaf depends primarily on the ¹⁸O content of leaf water which changes in response to relative humidity. In contrast, soil respiration caused a decline in the ¹⁸O values of forest CO₂ at night, because CO₂ released from the soil has equilibrated with soil water which has a lower ¹⁸O content than leaf water. The observed relationship between diurnal changes in CO₂ concentration and oxygen isotopic composition in the forest environment were consistent with a gas mixing model that considered the relative magnitudes of CO₂ fluxes associated with photosynthesis, respiration and turbulent exchange between the forest and the bulk atmosphere.     

Evaporation and carbonic anhydrase activity recorded in oxygen isotope signatures of net CO2 fluxes from a Mediterranean soil

The oxygen stable isotope composition (δ¹⁸O) of CO₂ is a valuable tool for studying the gas exchange between terrestrial ecosystems and the atmosphere. In the soil, it records the isotopic signal of water pools subjected to precipitation and evaporation events. The δ¹⁸O of the surface soil net CO₂ flux is dominated by the physical processes of diffusion of CO₂ into and out of the soil and the chemical reactions during CO₂–H₂O equilibration. Catalytic reactions by the enzyme carbonic anhydrase, reducing CO₂ hydration times, have been proposed recently to explain field observations of the δ¹⁸O signatures of net soil CO₂ fluxes. How important these catalytic reactions are for accurately predicting large‐scale biosphere fluxes and partitioning net ecosystem fluxes is currently uncertain because of the lack of field data. In this study, we determined the δ¹⁸O signatures of net soil CO₂ fluxes from soil chamber measurements in a Mediterranean forest. Over the 3 days of measurements, the observed δ¹⁸O signatures of net soil CO₂ fluxes became progressively enriched with a well‐characterized diurnal cycle. Model simulations indicated that the δ¹⁸O signatures recorded the interplay of two effects: (1) progressive enrichment of water in the upper soil by evaporation, and (2) catalytic acceleration of the isotopic exchange between CO₂ and soil water, amplifying the contributions of ‘atmospheric invasion’ to net signatures. We conclude that there is a need for better understanding of the role of enzymatic reactions, and hence soil biology, in determining the contributions of soil fluxes to oxygen isotope signals in atmospheric CO₂.