A method for deriving net primary productivity and component respiratory fluxes from tower-based eddy covariance data: a case study using a 17-year data record from a Douglas-fir chronosequence

Conventional gap‐filling procedures for eddy covariance (EC) data are limited to calculating ecosystem respiration (R_E) and gross ecosystem productivity (P_G) as well as missing values of net ecosystem productivity (F_NEP). We develop additional postprocessing steps that estimate net primary productivity (P_N), autotrophic (R_a), and heterotrophic respiration (R_h). This is based on conservation of mass of carbon (C), Monte Carlo (MC) simulation, and three ratios: C use efficiency (CUE, P_N to P_G), R_a to R_E, and F_NEP to R_E. This procedure, along with the estimation of F_NEP, R_E, and P_G, was applied to a Douglas‐fir dominated chronosequence on Vancouver Island, British Columbia, Canada. The EC data set consists of 17 site years from three sites: initiation (HDF00), pole/sapling (HDF88), and near mature (DF49), with stand ages from 1 to 56 years. Analysis focuses on annual C flux totals and C balance ratios as a function of stand age, assuming a rotation age of 56 years. All six C balance terms generally increased with stand age. Average annual P_N by stand was 213, 750, and 1261 g C m⁻² yr⁻¹ for HDF00, HDF88, and DF49, respectively. The canopy compensation point, the year when the chronosequence switched from a source to a sink of C, occurred at stand age ca. 20 years. HDF00 and HDF88 were strong and moderate sources (F_NEP=−581 and −138 g C m⁻² yr⁻¹), respectively, while DF49 was a moderate sink (F_NEP=294 g C m⁻² yr⁻¹) for C. Differences between sites were greater than interannual variation (IAV) within sites and highlighted the importance of age‐related effects in C cycling. The validity of the approach is discussed using a sensitivity analysis, a comparison with growth and yield estimates from the same chronosequence, and an intercomparison with other chronosequences.     

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.     

Net ecosystem production of a Douglas-fir stand for 3 years following clearcut harvesting

To investigate the variations in annual and seasonal net ecosystem production (F_NEP) during the development of a young forest, 3 years of continuous eddy covariance measurements of carbon dioxide (CO₂) fluxes were collected following clearcut harvesting and replanting of a coastal Douglas‐fir stand on the east coast of Vancouver Island, BC, Canada. The impact of changing weather and stand structure on F_NEP was examined by developing relationships between F_NEP and variables such as light, temperature, soil moisture, and leaf area index (LAI). In all 3 years, the stand was a large source of CO₂ (620, 520, and 600 g C m^?2 yr^?1 in the first, second, and third years, respectively). During this period, the growth of pioneer and understory species resulted in an increase in maximum growing season LAI from 0.2 in the year the seedlings were planted to 2.5 in the third year. The associated increase in annual gross ecosystem production (P=F_NEP?R_e, where R_e is ecosystem respiration) from 220 g C m^?2 yr^?1 in the first year to 640 g C m^?2 yr^?1 in the third year was exceeded by an increase in annual R_e from 840 to 1240 g C m^?2 yr^?1. Seasonal and interannual variations in daytime F_NEP and P were well described by variations in photosynthetically active radiation, temperature, and changes in LAI. Night‐time measurements of R_e exponentially increased with 2 cm soil temperature with an average Q₁₀ of 2 (relative increase in R_e for a 10°C increase in temperature) and R₁₀ (R_e at 10°C) that increased from 2.1 in the first year to 2.5 in the second year to 3.2 μmol m^?2 s^?1 in the third year. Although the re‐establishment of vegetation in this stand had a major impact on both P and R_e, interannual variations in weather also affected annual F_NEP. Drought, in the summer of the third year, resulted in early senescence and reduced both P and R_e. This resulted in more C being lost from the stand in the third year after harvesting than in the second year.     

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.     

Impact of severe dry season on net ecosystem exchange in the Neotropical rainforest of French Guiana

The lack of information on the ways seasonal drought modifies the CO₂ exchange between Neotropical rainforest ecosystems and the atmosphere and the resulting carbon balance hinders our ability to precisely predict how these ecosystems will respond as global environmental changes force them to face increasingly contrasting conditions in the future. To address this issue, seasonal variations in daily net ecosystem productivity (NEP_d) and two main components of this productivity, daily total ecosystem respiration (R_Ed) and daily gross ecosystem productivity (GEP_d), were estimated over 2 years at a flux tower site in French Guiana, South America (5 °16′54″N, 52 °54′44″W). We compared seasonal variations between wet and dry periods and between dry periods of contrasting levels of intensity (i.e. mild vs. severe) during equivalent 93‐day periods. During the wet periods, the ecosystem was almost in balance with the atmosphere (storage of 9.0 g C m^?2). Seasonal dry periods, regardless of their severity, are associated with higher incident radiation and lower R_Ed combined with reduced soil respiration associated with low soil water availability. During the mild dry period, as is normally the case in this region, the amount of carbon stored in the ecosystem was 32.7 g C m^?2. Severe drought conditions resulted in even lower R_Ed, whereas the photosynthetic activity was only moderately reduced and no change in canopy structure was observed. Thus, the severe dry period was characterized by greater carbon storage (64.6 g C m^?2), emphasizing that environmental conditions, such as during a severe drought, modify the CO₂ exchange between Neotropical rainforest ecosystems and the atmosphere and potentially the resulting carbon balance.     

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.     

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.     

Influence of physiological phenology on the seasonal pattern of ecosystem respiration in deciduous forests

Understanding the environmental and biotic drivers of respiration at the ecosystem level is a prerequisite to further improve scenarios of the global carbon cycle. In this study we investigated the relevance of physiological phenology, defined as seasonal changes in plant physiological properties, for explaining the temporal dynamics of ecosystem respiration (R_ECO) in deciduous forests. Previous studies showed that empirical R_ECO models can be substantially improved by considering the biotic dependency of R_ECO on the short‐term productivity (e.g., daily gross primary production, GPP) in addition to the well‐known environmental controls of temperature and water availability. Here, we use a model‐data integration approach to investigate the added value of physiological phenology, represented by the first temporal derivative of GPP, or alternatively of the fraction of absorbed photosynthetically active radiation, for modeling R_ECO at 19 deciduous broadleaved forests in the FLUXNET La Thuile database. The new data‐oriented semiempirical model leads to an 8% decrease in root mean square error (RMSE) and a 6% increase in the modeling efficiency (EF) of modeled R_ECO when compared to a version of the model that does not consider the physiological phenology. The reduction of the model‐observation bias occurred mainly at the monthly time scale, and in spring and summer, while a smaller reduction was observed at the annual time scale. The proposed approach did not improve the model performance at several sites, and we identified as potential causes the plant canopy heterogeneity and the use of air temperature as a driver of ecosystem respiration instead of soil temperature. However, in the majority of sites the model‐error remained unchanged regardless of the driving temperature. Overall, our results point toward the potential for improving current approaches for modeling R_ECO in deciduous forests by including the phenological cycle of the canopy.     

Semiempirical modeling of abiotic and biotic factors controlling ecosystem respiration across eddy covariance sites

In this study we examined ecosystem respiration (R_ECO) data from 104 sites belonging to FLUXNET, the global network of eddy covariance flux measurements. The goal was to identify the main factors involved in the variability of R_ECO: temporally and between sites as affected by climate, vegetation structure and plant functional type (PFT) (evergreen needleleaf, grasslands, etc.). We demonstrated that a model using only climate drivers as predictors of R_ECO failed to describe part of the temporal variability in the data and that the dependency on gross primary production (GPP) needed to be included as an additional driver of R_ECO. The maximum seasonal leaf area index (LAI_MAX) had an additional effect that explained the spatial variability of reference respiration (the respiration at reference temperature T_ref=15 °C, without stimulation introduced by photosynthetic activity and without water limitations), with a statistically significant linear relationship (r²=0.52, P<0.001, n=104) even within each PFT. Besides LAI_MAX, we found that reference respiration may be explained partially by total soil carbon content (SoilC). For undisturbed temperate and boreal forests a negative control of total nitrogen deposition (N_depo) on reference respiration was also identified. We developed a new semiempirical model incorporating abiotic factors (climate), recent productivity (daily GPP), general site productivity and canopy structure (LAI_MAX) which performed well in predicting the spatio‐temporal variability of R_ECO, explaining >70% of the variance for most vegetation types. Exceptions include tropical and Mediterranean broadleaf forests and deciduous broadleaf forests. Part of the variability in respiration that could not be described by our model may be attributed to a series of factors, including phenology in deciduous broadleaf forests and management practices in grasslands and croplands.     

Severe drought effects on ecosystem CO2 and H2O fluxes at three Mediterranean evergreen sites: revision of current hypotheses?

Eddy covariance and sapflow data from three Mediterranean ecosystems were analysed via top‐down approaches in conjunction with a mechanistic ecosystem gas‐exchange model to test current assumptions about drought effects on ecosystem respiration and canopy CO₂/H₂O exchange. The three sites include two nearly monospecific Quercus ilex L. forests – one on karstic limestone (Puéchabon), the other on fluvial sand with access to ground water (Castelporziano) – and a typical mixed macchia on limestone (Arca di Noè). Estimates of ecosystem respiration were derived from light response curves of net ecosystem CO₂ exchange. Subsequently, values of ecosystem gross carbon uptake were computed from eddy covariance CO₂ fluxes and estimates of ecosystem respiration as a function of soil temperature and moisture. Bulk canopy conductance was calculated by inversion of the Penman‐Monteith equation. In a top‐down analysis, it was shown that all three sites exhibit similar behaviour in terms of their overall response to drought. In contrast to common assumptions, at all sites ecosystem respiration revealed a decreasing temperature sensitivity ( Q ₁₀) in response to drought. Soil temperature and soil water content explained 70–80% of the seasonal variability of ecosystem respiration. During the drought, light‐saturated ecosystem gross carbon uptake and day‐time averaged canopy conductance declined by up to 90%. These changes were closely related to soil water content. Ecosystem water‐use efficiency of gross carbon uptake decreased during the drought, regardless whether evapotranspiration from eddy covariance or transpiration from sapflow had been used for the calculation. We evidence that this clearly contrasts current models of canopy function which predict increasing ecosystem water‐use efficiency (WUE) during the drought. Four potential explanations to those results were identified (patchy stomatal closure, changes in physiological capacities of photosynthesis, decreases in mesophyll conductance for CO₂, and photoinhibition), which will be tested in a forthcoming paper. It is suggested to incorporate the new findings into current biogeochemical models after further testing as this will improve estimates of climate change effects on (semi)arid ecosystems' carbon balances.     

Reduction of ecosystem productivity and respiration during the European summer 2003 climate anomaly: a joint flux tower, remote sensing and modelling analysis

The European CARBOEUROPE/FLUXNET monitoring sites, spatial remote sensing observations via the EOS‐MODIS sensor and ecosystem modelling provide independent and complementary views on the effect of the 2003 heatwave on the European biosphere's productivity and carbon balance. In our analysis, these data streams consistently demonstrate a strong negative anomaly of the primary productivity during the summer of 2003. FLUXNET eddy‐covariance data indicate that the drop in productivity was not primarily caused by high temperatures (‘heat stress’) but rather by limitation of water (drought stress) and that, contrary to the classical expectation about a heat wave, not only gross primary productivity but also ecosystem respiration declined by up to more than to 80 gC m⁻² month⁻¹. Anomalies of carbon and water fluxes were strongly correlated. While there are large between‐site differences in water‐use efficiency (WUE, 1–6 kg C kg⁻¹ H₂O) here defined as gross carbon uptake divided by evapotranspiration (WUE=GPP/ET), the year‐to‐year changes in WUE were small (<1 g kg⁻¹) and quite similar for most sites (i.e. WUE decreased during the year of the heatwave). Remote sensing data from MODIS and AVHRR both indicate a strong negative anomaly of the fraction of absorbed photosynthetically active radiation in summer 2003, at more than five standard deviations of the previous years. The spatial differentiation of this anomaly follows climatic and land‐use patterns: Largest anomalies occur in the centre of the meteorological anomaly (central Western Europe) and in areas dominated by crops or grassland. A preliminary model intercomparison along a gradient from data‐oriented models to process‐oriented models indicates that all approaches are similarly describing the spatial pattern of ecosystem sensitivity to the climatic 2003 event with major exceptions in the Alps and parts of Eastern Europe, but differed with respect to their interannual variability.     

Seasonal and annual respiration of a ponderosa pine ecosystem

The net ecosystem exchange of CO₂ between forests and the atmosphere, measured by eddy covariance, is the small difference between two large fluxes of photosynthesis and respiration. Chamber measurements of soil surface CO₂ efflux (F_s), wood respiration (F_w) and foliage respiration (F_f) help identify the contributions of these individual components to net ecosystem exchange. Models developed from the chamber data also provide independent estimates of respiration costs. We measured CO₂ efflux with chambers periodically in 1996–97 in a ponderosa pine forest in Oregon, scaled these measurements to the ecosystem, and computed annual totals for respiration by component. We also compared estimated half-hourly ecosystem respiration at night (F_nc) with eddy covariance measurements. Mean foliage respiration normalized to 10 °C was 0.20 μmol m^–2 (hemi-leaf surface area) s^–1, and reached a maximum of 0.24 μmol m^–2 HSA s^–1 between days 162 and 208. Mean wood respiration normalized to 10 °C was 5.9 μmol m^–3 sapwood s^–1, with slightly higher rates in mid-summer, when growth occurs. There was no significant difference (P > 0.10) between wood respiration of young (45 years) and old trees (250 years). Soil surface respiration normalized to 10 °C ranged from 0.7 to 3.0 μmol m^–2 (ground) s^–1 from days 23 to 329, with the lowest rates in winter and highest rates in late spring. Annual CO₂ flux from soil surface, foliage and wood was 683, 157, and 54 g C m^–2 y^–1, with soil fluxes responsible for 76% of ecosystem respiration. The ratio of net primary production to gross primary production was 0.45, consistent with values for conifer sites in Oregon and Australia, but higher than values reported for boreal coniferous forests. Below-ground carbon allocation (root turnover and respiration, estimated as F_s– litterfall carbon) consumed 61% of GPP; high ratios such as this are typical of sites with more water and nutrient constraints. The chamber estimates were moderately correlated with change in CO₂ storage in the canopy (F_stor) on calm nights (friction velocity u* < 0.25 m s^–1; R² = 0.60); F_stor was not significantly different from summed chamber estimates. On windy nights (u* > 0.25 m s^–1), the sum of turbulent flux measured above the canopy by eddy covariance and F_stor was only weakly correlated with summed chamber estimates (R² = 0.14); the eddy covariance estimates were lower than chamber estimates by 50%.     

Carbon storage and fluxes in ponderosa pine forests at different developmental stages

We compared carbon storage and fluxes in young and old ponderosa pine stands in Oregon, including plant and soil storage, net primary productivity, respiration fluxes, eddy flux estimates of net ecosystem exchange (NEE), and Biome‐BGC simulations of fluxes. The young forest (Y site) was previously an old‐growth ponderosa pine forest that had been clearcut in 1978, and the old forest (O site), which has never been logged, consists of two primary age classes (50 and 250 years old). Total ecosystem carbon content (vegetation, detritus and soil) of the O forest was about twice that of the Y site (21 vs. 10 kg C m^?2 ground), and significantly more of the total is stored in living vegetation at the O site (61% vs. 15%). Ecosystem respiration (R_e) was higher at the O site (1014 vs. 835 g C m^?2 year^?1), and it was largely from soils at both sites (77% of R_e). The biological data show that above‐ground net primary productivity (ANPP), NPP and net ecosystem production (NEP) were greater at the O site than the Y site. Monte Carlo estimates of NEP show that the young site is a source of CO₂ to the atmosphere, and is significantly lower than NEP(O) by c. 100 g C m^?2 year^?1. Eddy covariance measurements also show that the O site was a stronger sink for CO₂ than the Y site. Across a 15‐km swath in the region, ANPP ranged from 76 g C m^?2 year^?1 at the Y site to 236 g C m^?2 year^?1 (overall mean 158 ± 14 g C m^?2 year^?1). The lowest ANPP values were for the youngest and oldest stands, but there was a large range of ANPP for mature stands. Carbon, water and nitrogen cycle simulations with the Biome‐BGC model suggest that disturbance type and frequency, time since disturbance, age‐dependent changes in below‐ground allocation, and increasing atmospheric concentration of CO₂ all exert significant control on the net ecosystem exchange of carbon at the two sites. Model estimates of major carbon flux components agree with budget‐based observations to within ± 20%, with larger differences for NEP and for several storage terms. Simulations showed the period of regrowth required to replace carbon lost during and after a stand‐replacing fire (O) or a clearcut (Y) to be between 50 and 100 years. In both cases, simulations showed a shift from net carbon source to net sink (on an annual basis) 10–20 years after disturbance. These results suggest that the net ecosystem production of young stands may be low because heterotrophic respiration, particularly from soils, is higher than the NPP of the regrowth. The amount of carbon stored in long‐term pools (biomass and soils) in addition to short‐term fluxes has important implications for management of forests in the Pacific North‐west for carbon sequestration.     

Environmental variation is directly responsible for short- but not long-term variation in forest-atmosphere carbon exchange

Tower‐based eddy covariance measurements of forest‐atmosphere carbon dioxide (CO₂) exchange from many sites around the world indicate that there is considerable year‐to‐year variation in net ecosystem exchange (NEE). Here, we use a statistical modeling approach to partition the interannual variability in NEE (and its component fluxes, ecosystem respiration, R_eco, and gross photosynthesis, P_gross) into two main effects: variation in environmental drivers (air and soil temperature, solar radiation, vapor pressure deficit, and soil water content) and variation in the biotic response to this environmental forcing (as characterized by the model parameters). The model is applied to a 9‐year data set from the Howland AmeriFlux site, a spruce‐dominated forest in Maine, USA. Gap‐filled flux measurements at this site indicate that the forest has been sequestering, on average, 190 g C m⁻² yr⁻¹, with a range from 130 to 270 g C m⁻² yr⁻¹. Our fitted model predicts somewhat more uptake (mean 270 g C m⁻² yr⁻¹), but interannual variation is similar, and wavelet variance analyses indicate good agreement between tower measurements and model predictions across a wide range of timescales (hours to years). Associated with the interannual variation in NEE are clear differences among years in model parameters for both R_eco and P_gross. Analysis of model predictions suggests that, at the annual time step, about 40% of the variance in modeled NEE can be attributed to variation in environmental drivers, and 55% to variation in the biotic response to this forcing. As model predictions are aggregated at longer timescales (from individual days to months to calendar year), variation in environmental drivers becomes progressively less important, and variation in the biotic response becomes progressively more important, in determining the modeled flux. There is a strong negative correlation between modeled annual P_gross and R_eco (r=−0.93, P≤0.001); two possible explanations for this correlation are discussed. The correlation promotes homeostasis of NEE: the interannual variation in modeled NEE is substantially less than that for either P_gross or R_eco     

Carbon sequestration in boreal jack pine stands following harvesting

A large area of boreal jack pine (Pinus banksiana Lamb.) forest in Canada is recovering from clear‐cut harvesting, and the carbon (C) balance of these regenerating forests remains uncertain. Net ecosystem CO₂ exchange was measured using the eddy‐covariance technique at four jack pine sites representing different stages of stand development: three postharvest sites (HJP02, HJP94, and HJP75) and one preharvest site (OJP). The four sites, located in the southern Canadian boreal forest, Saskatchewan, Canada, are typical of low productivity jack pine stands and were 2, 10, 29, and 90 years old in 2004, respectively. Mean annual net ecosystem production (NEP) for 2004 and 2005 was ?137±11, 19±16, 73±28, and 22±30 g C m^?2 yr^?1 at HJP02, HJP94, HJP75 and OJP, respectively, showing the postharvest jack pine stands to be moderate C sources immediately after harvesting, weak sinks at 10 years, moderate C sinks at 30 years, then weak C sinks at 90 years. Mean annual gross ecosystem photosynthesis (GEP) for the 2 years was 96±10, 347±20, 576±34, and 583±35 g C m^?2 yr^?1 at HJP02, HJP94, HJP75, and OJP, respectively. The ratio of annual ecosystem respiration (R) to annual GEP was 2.51±0.15, 0.95±0.04, 0.87±0.03, and 0.96±0.03. Seasonally, NEP peaked in May or June at all four sites but GEP and R were highest in July. R at a reference soil temperature of 10 °C, ecosystem quantum yield and photosynthetic capacity were lowest for the 2‐year‐old stand. R was most sensitive to soil temperature for the 90‐year‐old stand. The primary source of variability in NEP over the course of succession of the jack pine ecosystem following harvesting was stand age due to the changes in leaf area index. Intersite variability in GEP and R was an order of magnitude greater than interannual variability at OJP. For both young and old stands, GEP had greater interannual variability than R and played a more important role than R in interannual variation in NEP. Based on year‐round flux measurements from 2000 to 2005, the 10‐year stand had larger interannual variability in GEP and R than the 90‐year stand. Interannual variability in NEP was driven primarily by early‐growing‐season temperature and growing‐season length. Photosynthesis played a dominant role in the rapid rise in NEP early in stand development. Late in stand development, however, the subtle decrease in NEP resulted primarily from increasing respiration.     

Net ecosystem carbon dioxide exchange over grazed steppe in central Mongolia

This paper presents results of 1 year (from March 25, 2003 to March 24, 2004, 366 days) of continuous measurements of net ecosystem CO₂ exchange (NEE) above a steppe in Mongolia using the eddy covariance technique. The steppe, typical of central Mongolia, is dominated by C₃ plants adapted to the continental climate. The following two questions are addressed: (1) how do NEE and its components: gross ecosystem production (GEP) and total ecosystem respiration (R_eco) vary seasonally? (2) how do NEE, GEP, and R_eco respond to biotic and abiotic factors? The hourly minimal NEE and the hourly maximal R_eco were −3.6 and 1.2 μmol m⁻² s⁻¹, respectively (negative values denoting net carbon uptake by the canopy from the atmosphere). Peak daily sums of NEE, GEP, and R_eco were −2.3, 3.5, and 1.5 g C m⁻² day⁻¹, respectively. The annual sums of GEP, R_eco, and NEE were 179, 138, and −41 g C m⁻², respectively. The carbon removal by sheep was estimated to range between 10 and 82 g C m⁻² yr⁻¹ using four different approaches. Including these estimates in the overall carbon budget yielded net ecosystem productivity of −23 to +20 g C m⁻² yr⁻¹. Thus, within the remaining experimental uncertainty the carbon budget at this steppe site can be considered to be balanced. For the growing period (from April 23 to October 21, 2003), 26% and 53% of the variation in daily NEE and GEP, respectively, could be explained by the changes in leaf area index. Seasonality of GEP, R_eco, and NEE was closely associated with precipitation, especially in the peak growing season when GEP and R_eco were largest. Water stress was observed in late July to early August, which switched the steppe from a carbon sink to a carbon source. For the entire growing period, the light response curves of daytime NEE showed a rather low apparent quantum yield (α=−0.0047 μmol CO₂ μmol⁻¹ photons of photosynthetically active radiation). However, the α values varied with air temperature (T_a), vapor pressure deficit, and soil water content.     

Sensitivity of canopy transpiration to altered precipitation in an upland oak forest: evidence from a long-term field manipulation study

Climate‐induced changes in regional precipitation could have important implications for the carbon, water, and nutrient cycles of forest ecosystems. However, few studies have examined the response of deciduous forests to increases or decreases in precipitation. Therefore, the throughfall displacement experiment (TDE) was established in 1993 near Oak Ridge, Tennessee to examine the sensitivity of an upland oak (Quercus spp.) forest to ambient, wet (+33%), and dry (?33%) precipitation regimes. Sap flux measurements on co‐occurring tree species were scaled using species‐specific estimates of stand sapwood area to derive daily and seasonal rates of canopy transpiration (E_C) from 2000 to 2003. With the exception of 2003, which was an extremely wet year, daily E_C in the dry plot, and occasionally during extended droughts in the ambient and wet plots, declined as water potential in the upper 0.35 m soil profile approached ‐3.0 MPa. Seasonal patterns of soil water potential and treatment‐specific differences in E_C were dependent on precipitation frequency and intensity. Supplemental precipitation added to the wet plot increased seasonal E_C on average by 9% (range ?1% to 19%), whereas extended periods of drought on the dry plot in 2000, 2001, and 2002 were sufficient to reduce seasonal E_C by 26–30% compared with the ambient plot. There was a strong correlation between seasonal E_C and the water stress integral, a cumulative index of drought severity and duration. A polynomial fitted to these data indicated that reductions in seasonal E_C on the order of 40% were possible given TDE‐imposed reductions in soil water potential. Application of this equation to all years of the TDE (1994–2003) revealed considerable interannual and treatment‐specific variation in canopy transpiration. In general, a 33% removal of throughfall on the dry plot during 1995, 1998, and 2002 resulted in a 23–32% reduction in seasonal E_C compared with the ambient plot. While droughts in deciduous forests are often limited in duration and tend to occur late in the growing season, soil water deficits of the magnitude observed in this study have the potential to impact local and regional forest water budgets.     

Chlorophyll fluorescence tracks seasonal variations of photosynthesis from leaf to canopy in a temperate forest

Accurate estimation of terrestrial gross primary productivity (GPP) remains a challenge despite its importance in the global carbon cycle. Chlorophyll fluorescence (ChlF) has been recently adopted to understand photosynthesis and its response to the environment, particularly with remote sensing data. However, it remains unclear how ChlF and photosynthesis are linked at different spatial scales across the growing season. We examined seasonal relationships between ChlF and photosynthesis at the leaf, canopy, and ecosystem scales and explored how leaf‐level ChlF was linked with canopy‐scale solar‐induced chlorophyll fluorescence (SIF) in a temperate deciduous forest at Harvard Forest, Massachusetts, USA. Our results show that ChlF captured the seasonal variations of photosynthesis with significant linear relationships between ChlF and photosynthesis across the growing season over different spatial scales (R² = 0.73, 0.77, and 0.86 at leaf, canopy, and satellite scales, respectively; P < 0.0001). We developed a model to estimate GPP from the tower‐based measurement of SIF and leaf‐level ChlF parameters. The estimation of GPP from this model agreed well with flux tower observations of GPP (R² = 0.68; P < 0.0001), demonstrating the potential of SIF for modeling GPP. At the leaf scale, we found that leaf F_q’/F_m’, the fraction of absorbed photons that are used for photochemistry for a light‐adapted measurement from a pulse amplitude modulation fluorometer, was the best leaf fluorescence parameter to correlate with canopy SIF yield (SIF/APAR, R² = 0.79; P < 0.0001). We also found that canopy SIF and SIF‐derived GPP (GPP_SIF) were strongly correlated to leaf‐level biochemistry and canopy structure, including chlorophyll content (R² = 0.65 for canopy GPP_SIF and chlorophyll content; P < 0.0001), leaf area index (LAI) (R² = 0.35 for canopy GPP_SIF and LAI; P < 0.0001), and normalized difference vegetation index (NDVI) (R² = 0.36 for canopy GPP_SIF and NDVI; P < 0.0001). Our results suggest that ChlF can be a powerful tool to track photosynthetic rates at leaf, canopy, and ecosystem scales.     

Air temperature triggers the recovery of evergreen boreal forest photosynthesis in spring

The timing of the commencement of photosynthesis (P^*) in spring is an important determinant of growing‐season length and thus of the productivity of boreal forests. Although controlled experiments have shed light on environmental mechanisms triggering release from photoinhibition after winter, quantitative research for trees growing naturally in the field is scarce. In this study, we investigated the environmental cues initiating the spring recovery of boreal coniferous forest ecosystems under field conditions. We used meteorological data and above‐canopy eddy covariance measurements of the net ecosystem CO₂ exchange (NEE) from five field stations located in northern and southern Finland, northern and southern Sweden, and central Siberia. The within‐ and intersite variability for P^* was large, 30–60 days. Of the different climate variables examined, air temperature emerged as the best predictor for P^* in spring. We also found that ‘soil thaw’, defined as the time when near‐surface soil temperature rapidly increases above 0°C, is not a useful criterion for P^*. In one case, photosynthesis commenced 1.5 months before soil temperatures increased significantly above 0°C. At most sites, we were able to determine a threshold for air‐temperature‐related variables, the exceeding of which was required for P^*. A 5‐day running‐average temperature (T₅) produced the best predictions, but a developmental‐stage model (S) utilizing a modified temperature sum concept also worked well. But for both T₅ and S, the threshold values varied from site to site, perhaps reflecting genetic differences among the stands or climate‐induced differences in the physiological state of trees in late winter/early spring. Only at the warmest site, in southern Sweden, could we obtain no threshold values for T₅ or S that could predict P^* reliably. This suggests that although air temperature appears to be a good predictor for P^* at high latitudes, there may be no unifying ecophysiological relationship applicable across the entire boreal zone.     

Modeling analysis of primary controls on net ecosystem productivity of seven boreal and temperate coniferous forests across a continental transect

Process‐based models are effective tools to synthesize and/or extrapolate measured carbon (C) exchanges from individual sites to large scales. In this study, we used a C‐ and nitrogen (N)‐cycle coupled ecosystem model named CN‐CLASS (Carbon Nitrogen‐Canadian Land Surface Scheme) to study the role of primary climatic controls and site‐specific C stocks on the net ecosystem productivity (NEP) of seven intermediate‐aged to mature coniferous forest sites across an east–west continental transect in Canada. The model was parameterized using a common set of parameters, except for two used in empirical canopy conductance–assimilation, and leaf area–sapwood relationships, and then validated using observed eddy covariance flux data. Leaf Rubisco‐N dynamics that are associated with soil–plant N cycling, and depend on canopy temperature, enabled the model to simulate site‐specific gross ecosystem productivity (GEP) reasonably well for all seven sites. Overall GEP simulations had relatively smaller differences compared with observations vs. ecosystem respiration (RE), which was the sum of many plant and soil components with larger variability and/or uncertainty associated with them. Both observed and simulated data showed that, on an annual basis, boreal forest sites were either carbon‐neutral or a weak C sink, ranging from 30 to 180 g C m^?2 yr^?1; while temperate forests were either a medium or strong C sink, ranging from 150 to 500 g C m^?2 yr^?1, depending on forest age and climatic regime. Model sensitivity tests illustrated that air temperature, among climate variables, and aboveground biomass, among major C stocks, were dominant factors impacting annual NEP. Vegetation biomass effects on annual GEP, RE and NEP showed similar patterns of variability at four boreal and three temperate forests. Air temperature showed different impacts on GEP and RE, and the response varied considerably from site to site. Higher solar radiation enhanced GEP, while precipitation differences had a minor effect. Magnitude of forest litter content and soil organic matter (SOM) affected RE. SOM also affected GEP, but only at low levels of SOM, because of low N mineralization that limited soil nutrient (N) availability. The results of this study will help to evaluate the impact of future climatic changes and/or forest C stock variations on C uptake and loss in forest ecosystems growing in diverse environments.