The role of interannual,seasonal, and synoptic climate on the carbon isotope ratio of ecosystem respiration at a semiarid woodland

The terrestrial carbon cycle is influenced by environmental variability at scales ranging from diurnal to interannual. Here, we present 5‐years of growing season (day 131–275) observations of the carbon isotope ratio of ecosystem respiration (δ¹³C_R) from a semiarid woodland. This ecosystem has a large necromass component resulting from 97%Pinus edulis mortality in 2002, is dominated by drought‐tolerant Juniperus monosperma trees, and experiences large variability in the timing and intensity of seasonal and synoptic water availability. Mean growing season δ¹³C_R was remarkably invariant (?23.57±0.4‰), with the exception of particularly enriched δ¹³C_R in 2006 following a winter with anomalously low snowfall. δ¹³C_R was strongly coupled to climate during premonsoon periods (～May to June), including fast (≤2 days) responses to changes in crown‐level stomatal conductance (G_c) and vapor pressure deficit (vpd) following rain pulses. In contrast, δ¹³C_R was relatively decoupled from G_c and environmental drivers during monsoon and postmonsoon periods (July–August and September, respectively), exhibiting only infrequent couplings of δ¹³C_R to vpd and soil water content (SWC) with longer lags (～8 days) and variable response slopes (both positive and negative). Notably, δ¹³C_R exhibited consistent dynamics after rainfall events, with depleted δ¹³C_R occurring within 1 h, progressive hourly δ¹³C_R enrichment over the remainder of the night, and net δ¹³C_R depletions over the multiple nights postevent in monsoon and postmonsoon periods. Overall this ecosystem demonstrated strong dependence of δ¹³C_R on precipitation, with an apparent dominance by the autotrophic δ¹³C signal in premonsoon periods when deep soil moisture is abundant and surface soil moisture is low, and weaker coupling during monsoonal periods consistent with increasing heterotrophic dominance when deep soil moisture has declined and surface moisture is variable.     

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.     

Associations between carbon isotope ratios of ecosystem respiration, water availability and canopy conductance

We tested the hypothesis that the stable carbon isotope signature of ecosystem respiration (δ¹³C_R) was regulated by canopy conductance (G_c) using weekly Keeling plots (n=51) from a semiarid old‐growth ponderosa pine (Pinus ponderosa) forest in Oregon, USA. For a comparison of forests in two contrasting climates we also evaluated trends in δ¹³C_R from a wet 20‐year‐old Douglas‐fir (Pseudotsuga menziesii) plantation located near the Pacific Ocean. Intraannual variability in δ¹³C_R was greater than 8.0‰ at both sites, was highest during autumn, winter, and spring when rainfall was abundant, and lowest during summer drought. The δ¹³C_R of the dry pine forest was consistently more positive than the wetter Douglas‐fir forest (mean annual δ¹³C_R: ?25.41‰ vs. ?26.23‰, respectively, P=0.07). At the Douglas‐fir forest, δ¹³C_R–climate relationships were consistent with predictions based on stomatal regulation of carbon isotope discrimination (Δ). Soil water content (SWC) and vapor pressure deficit (vpd) were the most important factors governing δ¹³C_R in this forest throughout the year. In contrast, δ¹³C_R at the pine forest was relatively insensitive to SWC or vpd, and exhibited a smaller drought‐related enrichment (～2‰) than the enrichment observed during drought at the Douglas‐fir forest (～5‰). Groundwater access at the pine forest may buffer canopy–gas exchange from drought. Despite this potential buffering, δ¹³C_R at the pine forest was significantly but weakly related to canopy conductance (G_c), suggesting that δ¹³C_R remains coupled to canopy–gas exchange despite groundwater access. During drought, δ¹³C_R was strongly correlated with soil temperature at both forests. The hypothesis that canopy‐level physiology is a critical regulator of δ¹³C_R was supported; however, belowground respiration may become more important during rain‐free periods.     

Short-term variation in the isotopic composition of organic matter allocated from the leaves to the stem of Pinus sylvestris: effects of photosynthetic and postphotosynthetic carbon isotope fractionation

We aimed to quantify the separate effects of photosynthetic and postphotosynthetic carbon isotope discrimination on δ¹³C of the fast‐turn‐over carbon pool (water soluble organic carbon and CO₂ emitted from heterotrophic tissues), including their diel variation, along the pathway of carbon transport from the foliage to the base of the stem. For that purpose, we determined δ¹³C in total and water‐soluble organic matter of the foliage plus δ¹³C and δ¹⁸O in phloem organic matter of twigs and at three heights along the stem of Pinus sylvestris over a nine‐day period, including four measurements per day. These data were related to meteorological and photosynthesis parameters and to the δ¹³C of stem‐emitted CO₂. In the canopy (foliage and twigs), the δ¹³C of soluble organic matter varied diurnally with amplitudes of up to 1.9‰. The greatest ¹³C enrichment was recorded during the night/early morning, indicating a strong influence of starch storage and remobilization on the carbon isotope signatures of sugars exported from the leaves. ¹³C enrichment of soluble organic matter from the leaves to the twig phloem and further on to the phloem of the stem was supposed to be a result of carbon isotope fractionation associated with metabolic processes in the source and sink tissues. CO₂ emitted from the stem was enriched by 2.3–5.2‰ compared with phloem organic matter. When day‐to‐day variation was addressed, water‐soluble leaf δ¹³C and twig phloem δ¹⁸O were strongly influenced by c_i/c_a and stomatal conductance (G_s), respectively. These results show that both photosynthetic and postphotosynthetic carbon isotope fractionation influence δ¹³C of organic matter over time, and over the length of the basipetal transport pathway. Clearly, these influences on the δ¹³C of respired CO₂ must be considered when using the latter for partitioning of ecosystem CO₂ fluxes or when the assessment of δ¹³C in organic matter is applied to estimate environmental effects in c_i/c_a.     

Use of decreasing foliar carbon isotope discrimination during water limitation as a carbon tracer to study whole plant carbon allocation

Foliar carbon isotope discrimination (Δ) of C₃ plants decreases in water‐deficit situations as discrimination by the photosynthetic primary carboxylation reaction decreases. This diminished Δ in leaves under water deficit can be used as a tracer to study whole plant carbon allocation patterns. Carbon isotope composition (δ¹³C value) of leaf hot water extracts or leaf tissue sap represents a short‐term integral of leaf carbon isotope discrimination and thus represents the δ¹³C value of source carbon that may be distributed within a plant in water‐deficit situations. By plotting the δ¹³C values of source carbon against the δ¹³C values of sink tissues, such as roots or stems, it is possible to assess carbon allocation to and incorporation into sink organs in relation to already present biomass. This natural abundance labelling method has been tested in three independent experiments, a one‐year field study with the fruit tree species Ziziphus mauritiana and peach (Prunus persica), a medium‐term drought stress experiment with Ziziphus rotundifolia trees in the glasshouse, and a short‐term drought stress experiment with soybean (Glycine max). The data show that the natural abundance labelling method can be applied to qualitatively assess carbon allocation in drought‐stressed plants. Although it is not possible to estimate exact fluxes of assimilated carbon during water deficit the method represents an easy to use tool to study integrated plant adaptations to drought stress. In addition, it is a less laborious method that can be applied in field studies as well as in controlled experiments, with plants from any developmental stage.     

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.     

Variation in the carbon and oxygen isotope composition of plant biomass and its relationship to water‐use efficiency at the leaf‐ and ecosystem‐scales in a northern Great Plains grassland

Measurements of the carbon (δ¹³C_m) and oxygen (δ¹⁸O_m) isotope composition of C₃ plant tissue provide important insights into controls on water‐use efficiency. We investigated the causes of seasonal and inter‐annual variability in water‐use efficiency in a grassland near Lethbridge, Canada using stable isotope (leaf‐scale) and eddy covariance measurements (ecosystem‐scale). The positive relationship between δ¹³C_m and δ¹⁸O_m values for samples collected during 1998–2001 indicated that variation in stomatal conductance and water stress‐induced changes in the degree of stomatal limitation of net photosynthesis were the major controls on variation in δ¹³C_m and biomass production during this time. By comparison, the lack of a significant relationship between δ¹³C_m and δ¹⁸O_m values during 2002, 2003 and 2006 demonstrated that water stress was not a significant limitation on photosynthesis and biomass production in these years. Water‐use efficiency was higher in 2000 than 1999, consistent with expectations because of greater stomatal limitation of photosynthesis and lower leaf c_i/ca during the drier conditions of 2000. Calculated values of leaf‐scale water‐use efficiency were 2–3 times higher than ecosystem‐scale water‐use efficiency, a difference that was likely due to carbon lost in root respiration and water lost during soil evaporation that was not accounted for by the stable isotope measurements.     

基于FORCCHN模型的中国典型森林生态系统碳通量模拟

森林生态系统是陆地碳循环的重要组成部分,其固碳能力显著高于其他陆地生态系统,研究森林生态系统碳通量是认识和理解全球变化对碳循环影响的关键。碳循环模型是研究森林生态系统碳通量有效工具。以长白山温带落叶阔叶林、千烟洲亚热带常绿针叶林、鼎湖山亚热带常绿阔叶林和西双版纳热带雨林等4种中国典型森林生态系统为研究对象,利用涡度相关2003-2012年观测数据,评估FORCCHN模型对生态系统呼吸（ER）,总初级生产力（GPP）,净生态系统生产力（NEP）的模型效果。结果表明：（1） FORCCHN模型能够较好的模拟中国4种典型森林生态系统不同时间尺度的碳通量。落叶阔叶林和常绿针叶林ER和GPP的逐日变化模拟效果较好（ER的相关系数分别为0.94和0.92,GPP的相关系数分别为0.86和0.74）;（2）4种森林生态系统碳通量季节动态模拟值和观测值显著相关（P<0.01）,ER、GPP、NEP的观测值和模拟值的R²分别为0.77-0.93、0.54-0.88和0.15-0.38;模型可以很好地模拟森林生态系统不同季节碳汇（NEP>0）,碳源（NEP<0）的变化规律;（3）4种森林生态系统碳通量模拟值与观测值的年际变化有很好的吻合度,但在数值大小上存在差异,模型高估了常绿阔叶林的ER和GPP,略微低估了其他3种森林生态系统ER和GPP。     

Comparing the performance of different stomatal conductance models using modelled and measured plant carbon isotope ratios (δ13C): implications for assessing physiological forcing

Accurate modelling of long‐term changes in plant stomatal functioning is vital to global climate change studies because changes in evapotranspiration influence temperature via physiological forcing of the climate. Various stomatal models are included in land surface schemes, but their robustness over longer timescales is difficult to validate. We compare the performance of three stomatal models, varying in their degree of complexity, and coupled to a land surface model. This is carried out by simulating the carbon isotope ratio of tree leaves (δ¹³C_leaf) over a period of 53 years, and comparing the results with carbon isotope ratios obtained from tree rings (δ¹³C_stem) measured at six sites in northern Europe. All three stomatal models fail to capture the observed interannual variability in the measured δ¹³C_stem time series. However, the Soil‐Plant‐Atmosphere (SPA) model performs significantly better than the Ball‐Berry (BB) or COX models when tested for goodness‐of‐fit against measured δ¹³C_stem. The δ¹³C_leaf time series simulated using the SPA model are significantly positively correlated (P < 0.05) with measured results over the full time period tested, at all six sites. The SPA model underestimates interannual variability measured in δ¹³C_stem, but is no worse than the BB model and significantly better than the COX model. The inability of current models to adequately replicate changes in stomatal response to rising levels of CO₂ concentrations, and thus to quantify the associated physiological forcing, warrants further investigation.     

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.     

Carbon cycling and storage in world forests: biome patterns related to forest age

Forest age, which is affected by stand‐replacing ecosystem disturbances (such as forest fires, harvesting, or insects), plays a distinguishing role in determining the distribution of carbon (C) pools and fluxes in different forested ecosystems. In this synthesis, net primary productivity (NPP), net ecosystem productivity (NEP), and five pools of C (living biomass, coarse woody debris, organic soil horizons, soil, and total ecosystem) are summarized by age class for tropical, temperate, and boreal forest biomes. Estimates of variability in NPP, NEP, and C pools are provided for each biome‐age class combination and the sources of variability are discussed. Aggregated biome‐level estimates of NPP and NEP were higher in intermediate‐aged forests (e.g., 30–120 years), while older forests (e.g., >120 years) were generally less productive. The mean NEP in the youngest forests (0–10 years) was negative (source to the atmosphere) in both boreal and temperate biomes (?0.1 and –1.9 Mg C ha^?1 yr^?1, respectively). Forest age is a highly significant source of variability in NEP at the biome scale; for example, mean temperate forest NEP was ?1.9, 4.5, 2.4, 1.9 and 1.7 Mg C ha^?1 yr^?1 across five age classes (0–10, 11–30, 31–70, 71–120, 121–200 years, respectively). In general, median NPP and NEP are strongly correlated (R²=0.83) across all biomes and age classes, with the exception of the youngest temperate forests. Using the information gained from calculating the summary statistics for NPP and NEP, we calculated heterotrophic soil respiration (R_h) for each age class in each biome. The mean R_h was high in the youngest temperate age class (9.7 Mg C ha^?1 yr^?1) and declined with age, implying that forest ecosystem respiration peaks when forests are young, not old. With notable exceptions, carbon pool sizes increased with age in all biomes, including soil C. Age trends in C cycling and storage are very apparent in all three biomes and it is clear that a better understanding of how forest age and disturbance history interact will greatly improve our fundamental knowledge of the terrestrial C cycle.     

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.     

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.     

Partitioning net ecosystem carbon exchange and the carbon isotopic disequilibrium in a subalpine forest

We investigate the utility of an improved isotopic method to partition the net ecosystem exchange of CO₂ (F) into net photosynthesis (F_A) and nonfoliar respiration (F_R). Measurements of F and the carbon isotopic content in air at a high‐elevation coniferous forest (the Niwot Ridge AmeriFlux site) were used to partition F into F_A and F_R. Isotopically partitioned fluxes were then compared with an independent flux partitioning method that estimated gross photosynthesis (GEE) and total ecosystem respiration (TER) based on statistical regressions of night‐time F and air temperature. We compared the estimates of F_A and F_R with expected canopy physiological relationships with light (photosynthetically active radiation) and air temperature. Estimates of F_A and GEE were dependent on light as expected, and TER, but not F_R, exhibited the expected dependence on temperature. Estimates of the isotopic disequilibrium D , or the difference between the isotopic signatures of net photosynthesis (δ_A, mean value ?24.6‰) and ecosystem respiration (δ_R, mean value ?25.1‰) were generally positive (δ_A>δ_R). The sign of D observed here is inconsistent with many other studies. The key parameters of the improved isotopic flux partitioning method presented here are ecosystem scale mesophyll conductance (g_m) and maximal vegetative stomatal conductance (g_cmax). The sensitivity analyses of F_A, F_R, and D to g_cmax indicated a critical value of g_cmax (0.15 mol m^?2 s^?1) above which estimates of F_A and F_R became larger in magnitude relative to GEE and TER. The value of D decreased with increasing values of g_m and g_cmax, but was still positive across all values of g_m and g_cmax. We conclude that the characterization of canopy‐scale mesophyll and stomatal conductances are important for further progress with the isotope partitioning method, and to confirm our anomalous isotopic disequilibrium findings.     

Strong seasonal disequilibrium measured between the oxygen isotope signals of leaf and soil CO2 exchange

The oxygen isotope composition (δ¹⁸O) of atmospheric CO₂ is among a very limited number of tools available to constrain estimates of the biospheric gross CO₂ fluxes, photosynthesis and respiration at large scales. However, the accuracy of the partitioning strongly depends on the extent of isotopic disequilibrium between the signals carried by these two gross fluxes. Chamber‐based field measurements of total CO₂ and CO¹⁸O fluxes from foliage and soil can help evaluate and refine our models of isotopic fractionation by plants and soils and validate the extent and pattern of isotopic disequilibrium within terrestrial ecosystems. Owing to sampling limitations in the past, such measurements have been very rare and covered only a few days. In this study, we coupled automated branch and soil chambers with tuneable diode laser absorption spectroscopy techniques to continuously capture the δ¹⁸O signals of foliage and soil CO₂ exchange in a Pinus pinaster Aït forest in France. Over the growing season, we observed a seasonally persistent isotopic disequilibrium between the δ¹⁸O signatures of net CO₂ fluxes from leaves and soils, except during rain events when the isotopic imbalance became temporarily weaker. Variations in the δ¹⁸O of CO₂ exchanged between leaves, soil and the atmosphere were well explained by theory describing changes in the oxygen isotope composition of ecosystem water pools in response to changes in leaf transpiration and soil evaporation.     

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₂.     

Spatial distribution of leaf water-use efficiency and carbon isotope discrimination within an isolated tree crown

The spatial variations in the stable carbon isotope composition (δ¹³C) of air and leaves (total matter and soluble sugars) were quantified within the crown of a well‐watered, 20‐year‐old walnut tree growing in a low‐density orchard. The observed leaf carbon isotope discrimination (Δ) was compared with that computed by a three‐dimensional model simulating the intracanopy distribution of irradiance, transpiration and photosynthesis (previously parameterized and tested for the same tree canopy) coupled to a biophysically based model of carbon isotope discrimination. The importance of discrimination associated with CO₂ gradients encountered from the substomatal sites to the carboxylation sites was evaluated. We also assessed by simulation the effect of current irradiance on leaf gas exchange and the effect of long‐term acclimation of photosynthetic capacity and stomatal and internal conductances to light regime on intracanopy gradients in Δ. The main conclusions of this study are: (i) leaf Δ can exhibit important variations (5 and 8‰ in total leaf material and soluble sugars, respectively) along light gradients within the foliage of an isolated tree; (ii) internal conductance must be taken into account to adequately predict leaf Δ, and (iii) the spatial variations in Δ and water‐use efficiency resulted from the short‐term response of leaf gas exchange to variations in local irradiance and, to a much lesser extent, from the long‐term acclimation of leaf characteristics to the local light regime.     

Seasonal and annual variation of carbon exchange in an evergreen Mediterranean forest in southern France

We present 9 years of eddy covariance measurements made over an evergreen Mediterranean forest in southern France. The goal of this study was to quantify the different components of the carbon (C) cycle, gross primary production (GPP) and ecosystem respiration (R_eco), and to assess the effects of climatic variables on these fluxes and on the net ecosystem exchange of carbon dioxide. The Puéchabon forest acted as a net C sink of ?254 g C m^?2 yr^?1, with a GPP of 1275 g C m^?2 yr^?1 and a R_eco of 1021 g C m^?2 yr^?1. On average, 83% of the net annual C sink occurred between March and June. The effects of exceptional events such the insect‐induced partial canopy defoliation that occurred in spring 2005, and the spring droughts of 2005 and 2006 are discussed. A high interannual variability of ecosystem C fluxes during summer and autumn was observed but the resulting effect on the annual net C budget was moderate. Increased severity and/or duration of summer drought under climate change do not appear to have the potential to negatively impact the average C budget of this ecosystem. On the contrary, factors affecting ecosystem functioning (drought and/or defoliation) during March–June period may reduce dramatically the annual C balance of evergreen Mediterranean forests.     

Assimilation exceeds respiration sensitivity to drought: A FLUXNET synthesis

The intensification of the hydrological cycle, with an observed and modeled increase in drought incidence and severity, underscores the need to quantify drought effects on carbon cycling and the terrestrial sink. FLUXNET, a global network of eddy covariance towers, provides dense data streams of meteorological data, and through flux partitioning and gap filling algorithms, estimates of net ecosystem productivity (F_NEP), gross ecosystem productivity (P), and ecosystem respiration (R). We analyzed the functional relationship of these three carbon fluxes relative to evaporative fraction (EF), an index of drought and site water status, using monthly data records from 238 micrometeorological tower sites distributed globally across 11 biomes. The analysis was based on relative anomalies of both EF and carbon fluxes and focused on drought episodes by biome and climatic season. Globally P was ≈50% more sensitive to a drought event than R. Network‐wide drought‐induced decreases in carbon flux averaged ?16.6 and ?9.3 g C m^?2 month^?1 for P and R, i.e., drought events induced a net decline in the terrestrial sink. However, in evergreen forests and wetlands drought was coincident with an increase in P or R during parts of the growing season. The most robust relationships between carbon flux and EF occurred during climatic spring for F_NEP and in climatic summer for P and R. Upscaling flux sensitivities to a global map showed that spatial patterns for all three carbon fluxes were linked to the distribution of croplands. Agricultural areas exhibited the highest sensitivity whereas the tropical region had minimal sensitivity to drought. Combining gridded flux sensitivities with their uncertainties and the spatial grid of FLUXNET revealed that a more robust quantification of carbon flux response to drought requires additional towers in all biomes of Africa and Asia as well as in the cropland, shrubland, savannah, and wetland biomes globally.     

Evaluating the carbon balance estimate from an automated ground‐level flux chamber system in artificial grass mesocosms

Measuring and modeling carbon (C) stock changes in terrestrial ecosystems are pivotal in addressing global C‐cycling model uncertainties. Difficulties in detecting small short‐term changes in relatively large C stocks require the development of robust sensitive flux measurement techniques. Net ecosystem exchange (NEE) ground‐level chambers are increasingly used to assess C dynamics in low vegetation ecosystems but, to date, have lacked formal rigorous field validation against measured C stock changes. We developed and deployed an automated and multiplexed C‐flux chamber system in grassland mesocosms in order rigorously to compare ecosystem total C budget obtained using hourly C‐flux measurements versus destructive net C balance. The system combines transparent NEE and opaque respiration chambers enabling partitioning of photosynthetic and respiratory fluxes. The C‐balance comparison showed good agreement between the two methods, but only after NEE fluxes were corrected for light reductions due to chamber presence. The dark chamber fluxes allowed assessing temperature sensitivity of ecosystem respiration (R_eco) components (i.e., heterotrophic vs. autotrophic) at different growth stages. We propose that such automated flux chamber systems can provide an accurate C balance, also enabling pivotal partitioning of the different C‐flux components (e.g., photosynthesis and respiration) suitable for model evaluation and developments.