Leonardo Montagnani

Respiration as the main determinant of carbon balance in European forests

Carbon exchange between the terrestrial biosphere and the atmosphere is one of the key processes that need to be assessed in the context of the Kyoto Protocol. Several studies suggest that the terrestrial biosphere is gaining carbon, but these estimates are obtained primarily by indirect methods, and the factors that control terrestrial carbon exchange, its magnitude and primary locations, are under debate. Here we present data of net ecosystem carbon exchange, collected between 1996 and 1998 from 15 European forests, which confirm that many European forest ecosystems act as carbon sinks. The annual carbon balances range from an uptake of 6.6 tonnes of carbon per hectare per year to a release of nearly 1 t C ha -1 yr-1, with a large variability between forests. The data show a significant increase of carbon uptake with decreasing latitude, whereas the gross primary production seems to be largely independent of latitude. Our observations indicate that, in general, ecosystem respiration determines net ecosystem carbon exchange. Also, for an accurate assessment of the carbon balance in a particular forest ecosystem, remote sensing of the normalized difference vegetation index or estimates based on forest inventories may not be sufficient.

Recent decline in the global land evapotranspiration trend due to limited moisture supply

More than half of the solar energy absorbed by land surfaces is currently used to evaporate water. Climate change is expected to intensify the hydrological cycle and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land—a key diagnostic criterion of the effects of climate change and variability—remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network, meteorological and remote-sensing observations, and a machine-learning algorithm. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

Global patterns of land‐atmosphere fluxes of carbon dioxide, latent heat, and sensible heat derived from eddy covariance, satellite, and meteorological observations

We upscaled FLUXNET observations of carbon dioxide, water, and energy fluxes to the global scale using the machine learning technique, model tree ensembles (MTE). We trained MTE to predict site-level gross primary productivity (GPP), terrestrial ecosystem respiration (TER), net ecosystem exchange (NEE), latent energy (LE), and sensible heat (H) based on remote sensing indices, climate and meteorological data, and information on land use. We applied the trained MTEs to generate global flux fields at a 0.5 degrees x 0.5 degrees spatial resolution and a monthly temporal resolution from 1982 to 2008. Cross-validation analyses revealed good performance of MTE in predicting among-site flux variability with modeling efficiencies (MEf) between 0.64 and 0.84, except for NEE (MEf = 0.32). Performance was also good for predicting seasonal patterns (MEf between 0.84 and 0.89, except for NEE (0.64)). By comparison, predictions of monthly anomalies were not as strong (MEf between 0.29 and 0.52). Improved accounting of disturbance and lagged environmental effects, along with improved characterization of errors in the training data set, would contribute most to further reducing uncertainties. Our global estimates of LE (158 +/- 7 J x 10(18) yr(-1)), H (164 +/- 15 J x 10(18) yr(-1)), and GPP (119 +/- 6 Pg C yr(-1)) were similar to independent estimates. Our global TER estimate (96 +/- 6 Pg C yr(-1)) was likely underestimated by 5-10%. Hot spot regions of interannual variability in carbon fluxes occurred in semiarid to semihumid regions and were controlled by moisture supply. Overall, GPP was more important to interannual variability in NEE than TER. Our empirically derived fluxes may be used for calibration and evaluation of land surface process models and for exploratory and diagnostic assessments of the biosphere.

Influence of spring and autumn phenological transitions on forest ecosystem productivity

We use eddy covariance measurements of net ecosystem productivity (NEP) from 21 FLUXNET sites (153 site-years of data) to investigate relationships between phenology and productivity (in terms of both NEP and gross ecosystem photosynthesis, GEP) in temperate and boreal forests. Results are used to evaluate the plausibility of four different conceptual models. Phenological indicators were derived from the eddy covariance time series, and from remote sensing and models. We examine spatial patterns (across sites) and temporal patterns (across years); an important conclusion is that it is likely that neither of these accurately represents how productivity will respond to future phenological shifts resulting from ongoing climate change. In spring and autumn, increased GEP resulting from an ‘extra’ day tends to be offset by concurrent, but smaller, increases in ecosystem respiration, and thus the effect on NEP is still positive. Spring productivity anomalies appear to have carry-over effects that translate to productivity anomalies in the following autumn, but it is not clear that these result directly from phenological anomalies. Finally, the productivity of evergreen needleleaf forests is less sensitive to phenology than is productivity of deciduous broadleaf forests. This has implications for how climate change may drive shifts in competition within mixed-species stands.

Global Convergence in the Temperature Sensitivity of Respiration at Ecosystem Level

Carbon Cycle and Climate Change As climate change accelerates, it is important to know the likely impact of climate change on the carbon cycle (see the Perspective by Reich). Gross primary production (GPP) is a measure of the amount of CO2 removed from the atmosphere every year to fuel photosynthesis. Beer et al. (p. 834, published online 5 July) used a combination of observation and calculation to estimate that the total GPP by terrestrial plants is around 122 billion tons per year; in comparison, burning fossil fuels emits about 7 billion tons annually. Thirty-two percent of this uptake occurs in tropical forests, and precipitation controls carbon uptake in more than 40% of vegetated land. The temperature sensitivity (Q10) of ecosystem respiratory processes is a key determinant of the interaction between climate and the carbon cycle. Mahecha et al. (p. 838, published online 5 July) now show that the Q10 of ecosystem respiration is invariant with respect to mean annual temperature, independent of the analyzed ecosystem type, with a global mean value for Q10 of 1.6. This level of temperature sensitivity suggests a less-pronounced climate sensitivity of the carbon cycle than assumed by recent climate models. The long-standing discrepancy between modeled and empirical measures of ecosystem temperature sensitivity is resolved. The respiratory release of carbon dioxide (CO2) from the land surface is a major flux in the global carbon cycle, antipodal to photosynthetic CO2 uptake. Understanding the sensitivity of respiratory processes to temperature is central for quantifying the climate–carbon cycle feedback. We approximated the sensitivity of terrestrial ecosystem respiration to air temperature (Q10) across 60 FLUXNET sites with the use of a methodology that circumvents confounding effects. Contrary to previous findings, our results suggest that Q10 is independent of mean annual temperature, does not differ among biomes, and is confined to values around 1.4 ± 0.1. The strong relation between photosynthesis and respiration, by contrast, is highly variable among sites. The results may partly explain a less pronounced climate–carbon cycle feedback than suggested by current carbon cycle climate models.

Redefinition and global estimation of basal ecosystem respiration rate

Basal ecosystem respiration rate (BR), the ecosystem respiration rate at a given temperature, is a common and important parameter in empirical models for quantifying ecosystem respiration (ER) globally. Numerous studies have indicated that BR varies in space. However, many empirical ER models still use a global constant BR largely due to the lack of a functional description for BR. In this study, we redefined BR to be ecosystem respiration rate at the mean annual temperature. To test the validity of this concept, we conducted a synthesis analysis using 276 site-years of eddy covariance data, from 79 research sites located at latitudes ranging from similar to 3 degrees S to similar to 70 degrees N. Results showed that mean annual ER rate closely matches ER rate at mean annual temperature. Incorporation of site-specific BR into global ER model substantially improved simulated ER compared to an invariant BR at all sites. These results confirm that ER at the mean annual temperature can be considered as BR in empirical models. A strong correlation was found between the mean annual ER and mean annual gross primary production (GPP). Consequently, GPP, which is typically more accurately modeled, can be used to estimate BR. A light use efficiency GPP model (i.e., EC-LUE) was applied to estimate global GPP, BR and ER with input data from MERRA (Modern Era Retrospective-Analysis for Research and Applications) and MODIS (Moderate resolution Imaging Spectroradiometer). The global ER was 103 Pg C yr (-1), with the highest respiration rate over tropical forests and the lowest value in dry and high-latitude areas.

Implications of the carbon cycle steady state assumption for biogeochemical modeling performance and inverse parameter retrieval

We analyze the impacts of the steady state assumption on inverse model parameter retrieval from biogeochemical models. An inverse model parameterization study using eddy covariance CO 2 flux data was performed with the Carnegie Ames Stanford Approach (CASA) model under conditions of strict and relaxed carbon cycle steady state assumption (CCSSA) in order to evaluate both the robustness of the models structure for the simulation of net ecosystem carbon fluxes and the assessment of the CCSSA effects on simulations and parameter estimation. Net ecosystem production (NEP) measurements from several eddy covariance sites were compared with NEP estimates from the CASA model driven by local weather station climate inputs as well as by remotely sensed fraction of photosynthetically active radiation absorbed by vegetation and leaf area index. The parameters considered for optimization are directly related to aboveground and belowground modeled responses to temperature and water availability, as well as a parameter (η) that relaxed the CCSSA in the model, allowing for site level simulations to be initialized either as net sinks or sources. A robust relationship was observed between NEP observations and predictions for most of the sites through the range of temporal scales considered (daily, weekly, biweekly, and monthly), supporting the conclusion that the model structure is able to capture the main processes explaining NEP variability. Overall, relaxing CCSSA increased model efficiency (21%) and decreased normalized average error (-92%). Intersite variability was a major source of variance in model performance differences between fixed (CCSSA f ) and relaxed (CCSSA r ) CCSSA conditions. These differences were correlated with mean annual NEP observations, where an average increase in modeling efficiency of 0.06 per 100 g Cm -2 a -1 (where a is years) of NEP is observed (α < 0.003). The parameter η was found to be a key parameter in the optimization exercise, generating significant model efficiency losses when removed from the initial parameter set and parameter uncertainties were significantly lower under CCSSA r . Moreover, modeled soil carbon stocks were generally closer to observations once the steady state assumption was relaxed. Finally, we also show that estimates of individual parameters are affected by the steady state assumption. For example, estimates of radiation-use efficiency were strongly affected by the CCSSA f indicating compensation effects for the inadequate steady state assumption, leading to effective and thus biased parameters. Overall, the importance of model structural evaluation in data assimilation approaches is thus emphasized.

Joint control of terrestrial gross primary productivity by plant phenology and physiology

Significance Terrestrial gross primary productivity (GPP), the total photosynthetic CO2 fixation at ecosystem level, fuels all life on land. However, its spatiotemporal variability is poorly understood, because GPP is determined by many processes related to plant phenology and physiological activities. In this study, we find that plant phenological and physiological properties can be integrated in a robust index—the product of the length of CO2 uptake period and the seasonal maximal photosynthesis—to explain the GPP variability over space and time in response to climate extremes and during recovery after disturbance. Terrestrial gross primary productivity (GPP) varies greatly over time and space. A better understanding of this variability is necessary for more accurate predictions of the future climate–carbon cycle feedback. Recent studies have suggested that variability in GPP is driven by a broad range of biotic and abiotic factors operating mainly through changes in vegetation phenology and physiological processes. However, it is still unclear how plant phenology and physiology can be integrated to explain the spatiotemporal variability of terrestrial GPP. Based on analyses of eddy–covariance and satellite-derived data, we decomposed annual terrestrial GPP into the length of the CO2 uptake period (CUP) and the seasonal maximal capacity of CO2 uptake (GPPmax). The product of CUP and GPPmax explained >90% of the temporal GPP variability in most areas of North America during 2000–2010 and the spatial GPP variation among globally distributed eddy flux tower sites. It also explained GPP response to the European heatwave in 2003 (r2 = 0.90) and GPP recovery after a fire disturbance in South Dakota (r2 = 0.88). Additional analysis of the eddy–covariance flux data shows that the interbiome variation in annual GPP is better explained by that in GPPmax than CUP. These findings indicate that terrestrial GPP is jointly controlled by ecosystem-level plant phenology and photosynthetic capacity, and greater understanding of GPPmax and CUP responses to environmental and biological variations will, thus, improve predictions of GPP over time and space.

Climatic Influences on Seasonal and Spatial Differences in Soil CO2 Efflux

The efflux of C02 from the soil is characterized by large seasonal fluctuations due to seasonal changes in root and microbial respiration. Although several biotic and abiotic factors influence root and microbial activity (see Chap. 3), the control exerted by temperature, and in some cases moisture, is usually dominant. In the absence of water stress, variation in soil temperature accounts for most of the seasonal and diurnal variation in soil C02 efflux. Where water stress frequently occurs, soil C02 efflux may not be correlated with soil temperature, but with its moisture content (Rout and Gupta 1989). Thus, C02 release from the soil appears to respond to temperature or moisture, whichever is most limiting at the time of measurement (Schlentner and van Cleve 1985).

Terrestrial gross primary production inferred from satellite fluorescence and vegetation models

Determining the spatial and temporal distribution of terrestrial gross primary production (GPP) is a critical step in closing the Earths carbon budget. Dynamical global vegetation models (DGVMs) provide mechanistic insight into GPP variability but diverge in predicting the response to climate in poorly investigated regions. Recent advances in the remote sensing of solar-induced chlorophyll fluorescence (SIF) opens up a new possibility to provide direct global observational constraints for GPP. Here, we apply an optimal estimation approach to infer the global distribution of GPP from an ensemble of eight DGVMs constrained by global measurements of SIF from the Greenhouse Gases Observing SATellite (GOSAT). These estimates are compared to flux tower data in N. America, Europe, and tropical S. America, with careful consideration of scale differences between models, GOSAT, and flux towers. Assimilation of GOSAT SIF with DGVMs causes a redistribution of global productivity from northern latitudes to the tropics of 7-8 Pg C yr(-1) from 2010 to 2012, with reduced GPP in northern forests (~3.6 Pg C yr(-1) ) and enhanced GPP in tropical forests (~3.7 Pg C yr(-1) ). This leads to improvements in the structure of the seasonal cycle, including earlier dry season GPP loss and enhanced peak-to-trough GPP in tropical forests within the Amazon Basin and reduced growing season length in northern croplands and deciduous forests. Uncertainty in predicted GPP (estimated from the spread of DGVMs) is reduced by 40-70% during peak productivity suggesting the assimilation of GOSAT SIF with models is well-suited for benchmarking. We conclude that satellite fluorescence augurs a new opportunity to quantify the GPP response to climate drivers and the potential to constrain predictions of carbon cycle evolution.