Search

The Amazon rainforest plays a crucial role in global carbon storage, but a minor destabilization of these forests could result in considerable carbon loss. Among the external factors affecting vegetation, elevated CO2 (eCO2) levels have long been anticipated to have positive impacts on vegetation, including the direct enhancement of both photosynthesis and productivity and increasing water use efficiency. However, the overall impact of eCO2 on the net carbon balance, especially concerning tree-mortality-induced carbon loss and recovery following extreme drought events, has remained elusive. Here, we use a process-based model that couples physiological CO2 effects with demography and both drought mortality and resistance processes. The model was previously calibrated to reproduce observed drought responses of Amazon forest sites. The model results, based on factorial simulations with and without eCO2, reveal that eCO2 enhances forest growth and promotes competition between trees, leading to more natural self-thinning of forest stands. This occurs following a growth–mortality trade-off response, although the growth outweighs the tree loss. Additionally, eCO2 provides water-saving benefits, reducing the risk of tree mortality during drought episodes. However, extra carbon losses could still occur due to an eCO2-induced increase in background biomass density, leading to “more carbon available to lose” when severe droughts happen. Furthermore, we found that eCO2 accelerates drought recovery and enhances drought resistance and resilience. By delving into the less-explored aspect of tree mortality response to eCO2, the model improvements advance our understanding of how carbon balance responds to eCO2, particularly regarding mechanisms of continuous competition-induced carbon loss vs. pulses of drought-induced carbon loss. These findings provide valuable insights into the intricate ways in which rising CO2 influences forest carbon dynamics and vulnerability, offering a critical understanding of the Amazon rainforest's evolution amidst more frequent and intense extreme climate events.

In early April 2021 several days of harsh frost affected central Europe. This led to very severe damage in grapevine and fruit trees in France, in regions where young leaves had already unfolded due to unusually warm temperatures in the preceding month (March 2021). We analysed with observations and 172 climate model simulations how human-induced climate change affected this event over central France, where many vineyards are located. We found that, without human-caused climate change, such temperatures in April or later in spring would have been even lower by 1.2 ∘C (0.75 to 1.7 ∘C). However, climate change also caused an earlier occurrence of bud burst that we characterized in this study by a growing degree day index value. This shift leaves young leaves exposed to more winter-like conditions with lower minimum temperatures and longer nights, an effect that overcompensates the warming effect. Extreme cold temperatures occurring after the start of the growing season such as those of April 2021 are now 2 ∘C colder (0.5 to 3.3 ∘C) than in preindustrial conditions, according to observations. This observed intensification of growing-period frosts is attributable, at least in part, to human-caused climate change with each of the five climate model ensembles used here simulating a cooling of growing-period annual temperature minima of 0.41 ∘C (0.22 to 0.60 ∘C) since preindustrial conditions. The 2021 growing-period frost event has become 50 % more likely (10 %–110 %). Models accurately simulate the observed warming in extreme lowest spring temperatures but underestimate the observed trends in growing-period frost intensities, a fact that yet remains to be explained. Model ensembles all simulate a further intensification of yearly minimum temperatures occurring in the growing period for future decades and a significant probability increase for such events of about 30 % (20 %–40 %) in a climate with global warming of 2 ∘C.

Extreme drought events in Amazon forests are expected to become more frequent and more intense with climate change, threatening ecosystem function and carbon balance. Yet large uncertainties exist on the resilience of this ecosystem to drought. A better quantification of tree hydraulics and mortality processes is needed to anticipate future drought effects on Amazon forests. Most state-of-the-art dynamic global vegetation models are relatively poor in their mechanistic description of these complex processes. Here, we implement a mechanistic plant hydraulic module within the ORCHIDEE-CAN-NHA r7236 land surface model to simulate the percentage loss of conductance (PLC) and changes in water storage among organs via a representation of the water potentials and vertical water flows along the continuum from soil to roots, stems and leaves. The model was evaluated against observed seasonal variability in stand-scale sap flow, soil moisture and productivity under both control and drought setups at the Caxiuanã throughfall exclusion field experiment in eastern Amazonia between 2001 and 2008. A relationship between PLC and tree mortality is built in the model from two empirical parameters, the cumulated duration of drought exposure that triggers mortality, and the mortality fraction in each day exceeding the exposure. Our model captures the large biomass drop in the year 2005 observed 4 years after throughfall reduction, and produces comparable annual tree mortality rates with observation over the study period. Our hydraulic architecture module provides promising avenues for future research in assimilating experimental data to parameterize mortality due to drought-induced xylem dysfunction. We also highlight that species-based (isohydric or anisohydric) hydraulic traits should be further tested to generalize the model performance in predicting the drought risks.

Over the last decades, the Amazon rainforest has been hit by multiple severe drought events. Here, we assess the severity and spatial extent of the extreme drought years 2005, 2010 and 2015/16 in the Amazon region and their impacts on the regional carbon cycle. As an indicator of drought stress in the Amazon rainforest, we use the widely applied maximum cumulative water deficit (MCWD). Evaluating nine state-of-the-art precipitation datasets for the Amazon region, we find that the spatial extent of the drought in 2005 ranges from 2.2 to 3.0 (mean =2.7) ×106 km2 (37 %–51 % of the Amazon basin, mean =45 %), where MCWD indicates at least moderate drought conditions (relative MCWD anomaly <-0.5). In 2010, the affected area was about 16 % larger, ranging from 3.0 up to 4.4 (mean =3.6) ×106 km2 (51 %–74 %, mean =61 %). In 2016, the mean area affected by drought stress was between 2005 and 2010 (mean =3.2×106 km2; 55 % of the Amazon basin), but the general disagreement between datasets was larger, ranging from 2.4 up to 4.1×106 km2 (40 %–69 %). In addition, we compare differences and similarities among datasets using the self-calibrating Palmer Drought Severity Index (scPDSI) and a dry-season rainfall anomaly index (RAI). We find that scPDSI shows a stronger and RAI a much weaker drought impact in terms of extent and severity for the year 2016 compared to MCWD. We further investigate the impact of varying evapotranspiration on the drought indicators using two state-of-the-art evapotranspiration datasets. Generally, the variability in drought stress is most dependent on the drought indicator (60 %), followed by the choice of the precipitation dataset (20 %) and the evapotranspiration dataset (20 %). Using a fixed, constant evapotranspiration rate instead of variable evapotranspiration can lead to an overestimation of drought stress in the parts of Amazon basin that have a more pronounced dry season (for example in 2010). We highlight that even for well-known drought events the spatial extent and intensity can strongly depend upon the drought indicator and the data sources it is calculated with. Using only one data source and drought indicator has the potential danger of under- or overestimating drought stress in regions with high measurement uncertainty, such as the Amazon basin.

In 2018 and 2019, central Europe was affected by two consecutive extreme dry and hot summers (DH18 and DH19). The DH18 event had severe impacts on ecosystems and likely affected vegetation activity in the subsequent year, for example through depletion of carbon reserves or damage from drought. Such legacies from drought and heat stress can further increase vegetation susceptibility to additional hazards. Temporally compound extremes such as DH18 and DH19 can, therefore, result in an amplification of impacts due to preconditioning effects of past disturbance legacies. Here, we evaluate how these two consecutive extreme summers impacted ecosystems in central Europe and how the vegetation responses to the first compound event (DH18) modulated the impacts of the second (DH19). To quantify changes in vegetation vulnerability to each compound event, we first train a set of statistical models for the period 2001–2017, which are then used to predict the impacts of DH18 and DH19 on enhanced vegetation index (EVI) anomalies from MODIS. These estimates correspond to expected EVI anomalies in DH18 and DH19 based on past sensitivity to climate. Large departures from the predicted values can indicate changes in vulnerability to dry and hot conditions and be used to identify modulating effects by vegetation activity and composition or other environmental factors on observed impacts. We find two regions in which the impacts of the two compound dry and hot (DH) events were significantly stronger than those expected based on previous climate–vegetation relationships. One region, largely dominated by grasslands and crops, showed much stronger impacts than expected in both DH events due to an amplification of their sensitivity to heat and drought, possibly linked to changing background CO2 and temperature conditions. A second region, dominated by forests and grasslands, showed browning from DH18 to DH19, even though dry and hot conditions were partly alleviated in 2019. This browning trajectory was mainly explained by the preconditioning role of DH18 on the impacts of DH19 due to interannual legacy effects and possibly by increased susceptibility to biotic disturbances, which are also promoted by warm conditions. Dry and hot summers are expected to become more frequent in the coming decades, posing a major threat to the stability of European forests. We show that state-of-the-art process-based models could not represent the decline in response to DH19 because they missed the interannual legacy effects from DH18 impacts. These gaps may result in an overestimation of the resilience and stability of temperate ecosystems in future model projections.

The Cretaceous evolution of angiosperm leaves towards higher vein densities enables unprecedented leaf stomatal conductance. Still, simulating and quantifying the impact of such change on plant productivity and transpiration in the specific environmental conditions of the Cretaceous remain challenging. Here, we address this issue by combining a paleo proxy-based model with a full atmosphere–vegetation model that couples stomatal conductance to carbon assimilation. Based on the fossil record, we build and evaluate three consistent proto-angiosperm vegetation parameterizations under two end-member scenarios of pCO2 (280 and 1120 ppm) for the mid-Cretaceous: a reduction of leaf hydraulic or photosynthetic capacity and a combination of both, supported by a likely coevolution of stomatal conductance and photosynthetic biochemistry. Our results suggest that decreasing leaf hydraulic and/or photosynthetic capacities always generates a reduction of transpiration that is predominantly the result of plant productivity variations modulated by light, water availability in the soil, atmospheric evaporative demand and pCO2. The high pCO2 acts as a fertilizer on plant productivity that strengthens plant transpiration and water-use efficiency. However, we show that proto-angiosperm physiology does not allow vegetation to grow under low pCO2 because of a positive feedback between leaf stomatal conductance and leaf area index. Our modeling approach stresses the need to better represent paleovegetation physiological traits. It also confirms the hypothesis of a likely evolution of angiosperms from a state of low leaf hydraulic and photosynthetic capacities at high pCO2 to a state of high leaf hydraulic and photosynthetic capacities linked to leaves with more and more veins together with a more efficient biochemistry at low pCO2.

Estimating the risk of forest collapse due to extreme climate events is one of the challenges of adapting to climate change. We adapt a concept from ruin theory, which is widely used in econometrics and the insurance industry, to design a growth–ruin model for trees which accounts for climate hazards that can jeopardize tree growth. This model is an elaboration of a classical Cramer–Lundberg ruin model that is used in the insurance industry. The model accounts for the interactions between physiological parameters of trees and the occurrence of climate hazards. The physiological parameters describe interannual growth rates and how trees react to hazards. The hazard parameters describe the probability distributions of the occurrence and intensity of climate events. We focus on a drought–heatwave hazard. The goal of the paper is to determine the dependence of the forest ruin and average growth probability distributions on physiological and hazard parameters. Using extensive Monte Carlo experiments, we show the existence of a threshold in the frequency of hazards beyond which forest ruin becomes certain to occur within a centennial horizon. We also detect a small effect of the strategies used to cope with hazards. This paper is a proof of concept for the quantification of forest collapse under climate change.

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4 yr−1 (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4 yr−1 or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4 yr−1 or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4 yr−1 larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4 yr−1 larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4 yr−1, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget, < 30∘ N) compared to mid-latitudes (∼ 30 %, 30–60∘ N) and high northern latitudes (∼ 4 %, 60–90∘ N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters. Some of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 Tg CH4 yr−1 lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 Tg CH4 yr−1 by 8 Tg CH4 yr−1, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5 % compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GCP-CH4-2019 (Saunois et al., 2020) and from the Global Carbon Project.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (EFF) are based on energy statistics and cement production data, while emissions from land use and land-use change (ELUC), mainly deforestation, are based on land use and land-use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (GATM) is computed from the annual changes in concentration. The ocean CO2 sink (SOCEAN) and terrestrial CO2 sink (SLAND) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (BIM), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2008–2017), EFF was 9.4±0.5 GtC yr−1, ELUC 1.5±0.7 GtC yr−1, GATM 4.7±0.02 GtC yr−1, SOCEAN 2.4±0.5 GtC yr−1, and SLAND 3.2±0.8 GtC yr−1, with a budget imbalance BIM of 0.5 GtC yr−1 indicating overestimated emissions and/or underestimated sinks. For the year 2017 alone, the growth in EFF was about 1.6 % and emissions increased to 9.9±0.5 GtC yr−1. Also for 2017, ELUC was 1.4±0.7 GtC yr−1, GATM was 4.6±0.2 GtC yr−1, SOCEAN was 2.5±0.5 GtC yr−1, and SLAND was 3.8±0.8 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2 concentration reached 405.0±0.1 ppm averaged over 2017. For 2018, preliminary data for the first 6–9 months indicate a renewed growth in EFF of +2.7 % (range of 1.8 % to 3.7 %) based on national emission projections for China, the US, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. The analysis presented here shows that the mean and trend in the five components of the global carbon budget are consistently estimated over the period of 1959–2017, but discrepancies of up to 1 GtC yr−1 persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations show (1) no consensus in the mean and trend in land-use change emissions, (2) a persistent low agreement among the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models, originating outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018, 2016, 2015a, b, 2014, 2013). All results presented here can be downloaded from https://doi.org/10.18160/GCP-2018.

Field measurements of aboveground net primary productivity (ANPP) in temperate grasslands suggest that both positive and negative asymmetric responses to changes in precipitation (P) may occur. Under normal range of precipitation variability, wet years typically result in ANPP gains being larger than ANPP declines in dry years (positive asymmetry), whereas increases in ANPP are lower in magnitude in extreme wet years compared to reductions during extreme drought (negative asymmetry). Whether the current generation of ecosystem models with a coupled carbon–water system in grasslands are capable of simulating these asymmetric ANPP responses is an unresolved question. In this study, we evaluated the simulated responses of temperate grassland primary productivity to scenarios of altered precipitation with 14 ecosystem models at three sites: Shortgrass steppe (SGS), Konza Prairie (KNZ) and Stubai Valley meadow (STU), spanning a rainfall gradient from dry to moist. We found that (1) the spatial slopes derived from modeled primary productivity and precipitation across sites were steeper than the temporal slopes obtained from inter-annual variations, which was consistent with empirical data; (2) the asymmetry of the responses of modeled primary productivity under normal inter-annual precipitation variability differed among models, and the mean of the model ensemble suggested a negative asymmetry across the three sites, which was contrary to empirical evidence based on filed observations; (3) the mean sensitivity of modeled productivity to rainfall suggested greater negative response with reduced precipitation than positive response to an increased precipitation under extreme conditions at the three sites; and (4) gross primary productivity (GPP), net primary productivity (NPP), aboveground NPP (ANPP) and belowground NPP (BNPP) all showed concave-down nonlinear responses to altered precipitation in all the models, but with different curvatures and mean values. Our results indicated that most models overestimate the negative drought effects and/or underestimate the positive effects of increased precipitation on primary productivity under normal climate conditions, highlighting the need for improving eco-hydrological processes in those models in the future.

Simulation of vegetation–climate feedbacks in high latitudes in the ORCHIDEE land surface model was improved by the addition of three new circumpolar plant functional types (PFTs), namely non-vascular plants representing bryophytes and lichens, Arctic shrubs and Arctic C3 grasses. Non-vascular plants are assigned no stomatal conductance, very shallow roots, and can desiccate during dry episodes and become active again during wet periods, which gives them a larger phenological plasticity (i.e. adaptability and resilience to severe climatic constraints) compared to grasses and shrubs. Shrubs have a specific carbon allocation scheme, and differ from trees by their larger survival rates in winter, due to protection by snow. Arctic C3 grasses have the same equations as in the original ORCHIDEE version, but different parameter values, optimised from in situ observations of biomass and net primary productivity (NPP) in Siberia. In situ observations of living biomass and productivity from Siberia were used to calibrate the parameters of the new PFTs using a Bayesian optimisation procedure. With the new PFTs, we obtain a lower NPP by 31 % (from 55° N), as well as a lower roughness length (−41 %), transpiration (−33 %) and a higher winter albedo (by +3.6 %) due to increased snow cover. A simulation of the water balance and runoff and drainage in the high northern latitudes using the new PFTs results in an increase of fresh water discharge in the Arctic ocean by 11 % (+140 km3 yr−1), owing to less evapotranspiration. Future developments should focus on the competition between these three PFTs and boreal tree PFTs, in order to simulate their area changes in response to climate change, and the effect of carbon–nitrogen interactions.

The use of dynamic global vegetation models (DGVMs) to estimate CO2 emissions from land-use and land-cover change (LULCC) offers a new window to account for spatial and temporal details of emissions and for ecosystem processes affected by LULCC. One drawback of LULCC emissions from DGVMs, however, is lack of observation constraint. Here, we propose a new method of using satellite- and inventory-based biomass observations to constrain historical cumulative LULCC emissions (ELUCc) from an ensemble of nine DGVMs based on emerging relationships between simulated vegetation biomass and ELUCc. This method is applicable on the global and regional scale. The original DGVM estimates of ELUCc range from 94 to 273 PgC during 1901–2012. After constraining by current biomass observations, we derive a best estimate of 155 ± 50 PgC (1σ Gaussian error). The constrained LULCC emissions are higher than prior DGVM values in tropical regions but significantly lower in North America. Our emergent constraint approach independently verifies the median model estimate by biomass observations, giving support to the use of this estimate in carbon budget assessments. The uncertainty in the constrained ELUCc is still relatively large because of the uncertainty in the biomass observations, and thus reduced uncertainty in addition to increased accuracy in biomass observations in the future will help improve the constraint. This constraint method can also be applied to evaluate the impact of land-based mitigation activities.

The global methane (CH4) budget is becoming an increasingly important component for managing realistic pathways to mitigate climate change. This relevance, due to a shorter atmospheric lifetime and a stronger warming potential than carbon dioxide, is challenged by the still unexplained changes of atmospheric CH4 over the past decade. Emissions and concentrations of CH4 are continuing to increase, making CH4 the second most important human-induced greenhouse gas after carbon dioxide. Two major difficulties in reducing uncertainties come from the large variety of diffusive CH4 sources that overlap geographically, and from the destruction of CH4 by the very short-lived hydroxyl radical (OH). To address these difficulties, we have established a consortium of multi-disciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate research on the methane cycle, and producing regular (∼ biennial) updates of the global methane budget. This consortium includes atmospheric physicists and chemists, biogeochemists of surface and marine emissions, and socio-economists who study anthropogenic emissions. Following Kirschke et al. (2013), we propose here the first version of a living review paper that integrates results of top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models, inventories and data-driven approaches (including process-based models for estimating land surface emissions and atmospheric chemistry, and inventories for anthropogenic emissions, data-driven extrapolations). For the 2003–2012 decade, global methane emissions are estimated by top-down inversions at 558 Tg CH4 yr−1, range 540–568. About 60 % of global emissions are anthropogenic (range 50–65 %). Since 2010, the bottom-up global emission inventories have been closer to methane emissions in the most carbon-intensive Representative Concentrations Pathway (RCP8.5) and higher than all other RCP scenarios. Bottom-up approaches suggest larger global emissions (736 Tg CH4 yr−1, range 596–884) mostly because of larger natural emissions from individual sources such as inland waters, natural wetlands and geological sources. Considering the atmospheric constraints on the top-down budget, it is likely that some of the individual emissions reported by the bottom-up approaches are overestimated, leading to too large global emissions. Latitudinal data from top-down emissions indicate a predominance of tropical emissions (∼ 64 % of the global budget, 4 yr−1, range 51–72, −14 %) and higher emissions in Africa (86 Tg CH4 yr−1, range 73–108, +19 %) than bottom-up values used as prior estimates. Overall, uncertainties for anthropogenic emissions appear smaller than those from natural sources, and the uncertainties on source categories appear larger for top-down inversions than for bottom-up inventories and models. The most important source of uncertainty on the methane budget is attributable to emissions from wetland and other inland waters. We show that the wetland extent could contribute 30–40 % on the estimated range for wetland emissions. Other priorities for improving the methane budget include the following: (i) the development of process-based models for inland-water emissions, (ii) the intensification of methane observations at local scale (flux measurements) to constrain bottom-up land surface models, and at regional scale (surface networks and satellites) to constrain top-down inversions, (iii) improvements in the estimation of atmospheric loss by OH, and (iv) improvements of the transport models integrated in top-down inversions. The data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (http://doi.org/10.3334/CDIAC/GLOBAL_METHANE_BUDGET_2016_V1.1) and the Global Carbon Project.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2006–2015), EFF was 9.3 ± 0.5 GtC yr−1, ELUC 1.0 ± 0.5 GtC yr−1, GATM 4.5 ± 0.1 GtC yr−1, SOCEAN 2.6 ± 0.5 GtC yr−1, and SLAND 3.1 ± 0.9 GtC yr−1. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9 ± 0.5 GtC yr−1, showing a slowdown in growth of these emissions compared to the average growth of 1.8 % yr−1 that took place during 2006–2015. Also, for 2015, ELUC was 1.3 ± 0.5 GtC yr−1, GATM was 6.3 ± 0.2 GtC yr−1, SOCEAN was 3.0 ± 0.5 GtC yr−1, and SLAND was 1.9 ± 0.9 GtC yr−1. GATM was higher in 2015 compared to the past decade (2006–2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4 ± 0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2 % (range of −1.0 to +1.8 %) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Niño conditions of 2015–2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565 ± 55 GtC (2075 ± 205 GtCO2) for 1870–2016, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015b, a, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2016).

The high-resolution CO2 record from Law Dome ice core reveals that atmospheric CO2 concentration stalled during the 1940s (so-called CO2 plateau). Since the fossil-fuel emissions did not decrease during the period, this stalling implies the persistence of a strong sink, perhaps sustained for as long as a decade or more. Double-deconvolution analyses have attributed this sink to the ocean, conceivably as a response to the very strong El Niño event in 1940–1942. However, this explanation is questionable, as recent ocean CO2 data indicate that the range of variability in the ocean sink has been rather modest in recent decades, and El Niño events have generally led to higher growth rates of atmospheric CO2 due to the offsetting terrestrial response. Here, we use the most up-to-date information on the different terms of the carbon budget: fossil-fuel emissions, four estimates of land-use change (LUC) emissions, ocean uptake from two different reconstructions, and the terrestrial sink modelled by the TRENDY project to identify the most likely causes of the 1940s plateau. We find that they greatly overestimate atmospheric CO2 growth rate during the plateau period, as well as in the 1960s, in spite of giving a plausible explanation for most of the 20th century carbon budget, especially from 1970 onwards. The mismatch between reconstructions and observations during the CO2 plateau epoch of 1940–1950 ranges between 0.9 and 2.0 Pg C yr−1, depending on the LUC dataset considered. This mismatch may be explained by (i) decadal variability in the ocean carbon sink not accounted for in the reconstructions we used, (ii) a further terrestrial sink currently missing in the estimates by land-surface models, or (iii) LUC processes not included in the current datasets. Ocean carbon models from CMIP5 indicate that natural variability in the ocean carbon sink could explain an additional 0.5 Pg C yr−1 uptake, but it is unlikely to be higher. The impact of the 1940–1942 El Niño on the observed stabilization of atmospheric CO2 cannot be confirmed nor discarded, as TRENDY models do not reproduce the expected concurrent strong decrease in terrestrial uptake. Nevertheless, this would further increase the mismatch between observed and modelled CO2 growth rate during the CO2 plateau epoch. Tests performed using the OSCAR (v2.2) model indicate that changes in land use not correctly accounted for during the period (coinciding with drastic socioeconomic changes during the Second World War) could contribute to the additional sink required. Thus, the previously proposed ocean hypothesis for the 1940s plateau cannot be confirmed by independent data. Further efforts are required to reduce uncertainty in the different terms of the carbon budget during the first half of the 20th century and to better understand the long-term variability of the ocean and terrestrial CO2 sinks.

We examined the net terrestrial carbon flux to the atmosphere (FTA) simulated by nine models from the TRENDY dynamic global vegetation model project for its seasonal cycle and amplitude trend during 1961–2012. While some models exhibit similar phase and amplitude compared to atmospheric inversions, with spring drawdown and autumn rebound, others tend to rebound early in summer. The model ensemble mean underestimates the magnitude of the seasonal cycle by 40 % compared to atmospheric inversions. Global FTA amplitude increase (19 ± 8 %) and its decadal variability from the model ensemble are generally consistent with constraints from surface atmosphere observations. However, models disagree on attribution of this long-term amplitude increase, with factorial experiments attributing 83 ± 56 %, −3 ± 74 and 20 ± 30 % to rising CO2, climate change and land use/cover change, respectively. Seven out of the nine models suggest that CO2 fertilization is the strongest control – with the notable exception of VEGAS, which attributes approximately equally to the three factors. Generally, all models display an enhanced seasonality over the boreal region in response to high-latitude warming, but a negative climate contribution from part of the Northern Hemisphere temperate region, and the net result is a divergence over climate change effect. Six of the nine models show that land use/cover change amplifies the seasonal cycle of global FTA: some are due to forest regrowth, while others are caused by crop expansion or agricultural intensification, as revealed by their divergent spatial patterns. We also discovered a moderate cross-model correlation between FTA amplitude increase and increase in land carbon sink (R2 = 0.61). Our results suggest that models can show similar results in some benchmarks with different underlying mechanisms; therefore, the spatial traits of CO2 fertilization, climate change and land use/cover changes are crucial in determining the right mechanisms in seasonal carbon cycle change as well as mean sink change.

The Land Surface, Snow and Soil Moisture Model Intercomparison Project (LS3MIP) is designed to provide a comprehensive assessment of land surface, snow and soil moisture feedbacks on climate variability and climate change, and to diagnose systematic biases in the land modules of current Earth system models (ESMs). The solid and liquid water stored at the land surface has a large influence on the regional climate, its variability and predictability, including effects on the energy, water and carbon cycles. Notably, snow and soil moisture affect surface radiation and flux partitioning properties, moisture storage and land surface memory. They both strongly affect atmospheric conditions, in particular surface air temperature and precipitation, but also large-scale circulation patterns. However, models show divergent responses and representations of these feedbacks as well as systematic biases in the underlying processes. LS3MIP will provide the means to quantify the associated uncertainties and better constrain climate change projections, which is of particular interest for highly vulnerable regions (densely populated areas, agricultural regions, the Arctic, semi-arid and other sensitive terrestrial ecosystems). The experiments are subdivided in two components, the first addressing systematic land biases in offline mode (“LMIP”, building upon the 3rd phase of Global Soil Wetness Project; GSWP3) and the second addressing land feedbacks attributed to soil moisture and snow in an integrated framework (“LFMIP”, building upon the GLACE-CMIP blueprint).

Grassland management type (grazed or mown) and intensity (intensive or extensive) play a crucial role in the greenhouse gas balance and surface energy budget of this biome, both at field scale and at large spatial scale. However, global gridded historical information on grassland management intensity is not available. Combining modelled grass-biomass productivity with statistics of the grass-biomass demand by livestock, we reconstruct gridded maps of grassland management intensity from 1901 to 2012. These maps include the minimum area of managed vs. maximum area of unmanaged grasslands and the fraction of mown vs. grazed area at a resolution of 0.5° by 0.5°. The grass-biomass demand is derived from a livestock dataset for 2000, extended to cover the period 1901–2012. The grass-biomass supply (i.e. forage grass from mown grassland and biomass grazed) is simulated by the process-based model ORCHIDEE-GM driven by historical climate change, rising CO2 concentration, and changes in nitrogen fertilization. The global area of managed grassland obtained in this study increases from 6.1 × 106 km2 in 1901 to 12.3 × 106 km2 in 2000, although the expansion pathway varies between different regions. ORCHIDEE-GM also simulated augmentation in global mean productivity and herbage-use efficiency over managed grassland during the 20th century, indicating a general intensification of grassland management at global scale but with regional differences. The gridded grassland management intensity maps are model dependent because they depend on modelled productivity. Thus specific attention was given to the evaluation of modelled productivity against a series of observations from site-level net primary productivity (NPP) measurements to two global satellite products of gross primary productivity (GPP) (MODIS-GPP and SIF data). Generally, ORCHIDEE-GM captures the spatial pattern, seasonal cycle, and interannual variability of grassland productivity at global scale well and thus is appropriate for global applications presented here.

The response of crops to changing climate and atmospheric CO2 concentration ([CO2]) could have large effects on food production, and impact carbon, water, and energy fluxes, causing feedbacks to the climate. To simulate the response of temperate crops to changing climate and [CO2], which accounts for the specific phenology of crops mediated by management practice, we describe here the development of a process-oriented terrestrial biogeochemical model named ORCHIDEE-CROP (v0), which integrates a generic crop phenology and harvest module, and a very simple parameterization of nitrogen fertilization, into the land surface model (LSM) ORCHIDEEv196, in order to simulate biophysical and biochemical interactions in croplands, as well as plant productivity and harvested yield. The model is applicable for a range of temperate crops, but is tested here using maize and winter wheat, with the phenological parameterizations of two European varieties originating from the STICS agronomical model. We evaluate the ORCHIDEE-CROP (v0) model against eddy covariance and biometric measurements at seven winter wheat and maize sites in Europe. The specific ecosystem variables used in the evaluation are CO2 fluxes (net ecosystem exchange, NEE), latent heat, and sensible heat fluxes. Additional measurements of leaf area index (LAI) and aboveground biomass and yield are used as well. Evaluation results revealed that ORCHIDEE-CROP (v0) reproduced the observed timing of crop development stages and the amplitude of the LAI changes. This is in contrast to ORCHIDEEv196 where, by default, crops have the same phenology as grass. A halving of the root mean square error for LAI from 2.38 ± 0.77 to 1.08 ± 0.34 m2 m−2 was obtained when ORCHIDEEv196 and ORCHIDEE-CROP (v0) were compared across the seven study sites. Improved crop phenology and carbon allocation led to a good match between modeled and observed aboveground biomass (with a normalized root mean squared error (NRMSE) of 11.0–54.2 %), crop yield, daily carbon and energy fluxes (with a NRMSE of ∼ 9.0–20.1 and ∼ 9.4–22.3 % for NEE), and sensible and latent heat fluxes. The simulated yields for winter wheat and maize from ORCHIDEE-CROP (v0) showed a good match with the simulated results from STICS for three sites with available crop yield observations, where the average NRMSE was ∼ 8.8 %. The model data misfit for energy fluxes were within the uncertainties of the measurements, which themselves showed an incomplete energy balance closure within the range 80.6–86.3 %. The remaining discrepancies between the modeled and observed LAI and other variables at specific sites were partly attributable to unrealistic representations of management events by the model. ORCHIDEE-CROP (v0) has the ability to capture the spatial gradients of carbon and energy-related variables, such as gross primary productivity, NEE, and sensible and latent heat fluxes across the sites in Europe, which is an important requirement for future spatially explicit simulations. Further improvement of the model, with an explicit parameterization of nutritional dynamics and management, is expected to improve its predictive ability to simulate croplands in an Earth system model.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates as well as consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover-change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2, and land-cover change (some including nitrogen–carbon interactions). We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2005–2014), EFF was 9.0 ± 0.5 GtC yr−1, ELUC was 0.9 ± 0.5 GtC yr−1, GATM was 4.4 ± 0.1 GtC yr−1, SOCEAN was 2.6 ± 0.5 GtC yr−1, and SLAND was 3.0 ± 0.8 GtC yr−1. For the year 2014 alone, EFF grew to 9.8 ± 0.5 GtC yr−1, 0.6 % above 2013, continuing the growth trend in these emissions, albeit at a slower rate compared to the average growth of 2.2 % yr−1 that took place during 2005–2014. Also, for 2014, ELUC was 1.1 ± 0.5 GtC yr−1, GATM was 3.9 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and SLAND was 4.1 ± 0.9 GtC yr−1. GATM was lower in 2014 compared to the past decade (2005–2014), reflecting a larger SLAND for that year. The global atmospheric CO2 concentration reached 397.15 ± 0.10 ppm averaged over 2014. For 2015, preliminary data indicate that the growth in EFF will be near or slightly below zero, with a projection of −0.6 [range of −1.6 to +0.5] %, based on national emissions projections for China and the USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the global economy for the rest of the world. From this projection of EFF and assumed constant ELUC for 2015, cumulative emissions of CO2 will reach about 555 ± 55 GtC (2035 ± 205 GtCO2) for 1870–2015, about 75 % from EFF and 25 % from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set (Le Quéré et al., 2015, 2014, 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2015).

Processes that describe the distribution of vegetation and ecosystem succession after disturbance are an important component of dynamic global vegetation models (DGVMs). The vegetation dynamics module (ORC-VD) within the process-based ecosystem model ORCHIDEE (Organizing Carbon and Hydrology in Dynamic Ecosystems) has not been updated and evaluated since many years and is known to produce unrealistic results. This study presents a new parameterization of ORC-VD for mid- to high-latitude regions in the Northern Hemisphere, including processes that influence the existence, mortality and competition between tree functional types. A new set of metrics is also proposed to quantify the performance of ORC-VD, using up to five different data sets of satellite land cover, forest biomass from remote sensing and inventories, a data-driven estimate of gross primary productivity (GPP) and two gridded data sets of soil organic carbon content. The scoring of ORC-VD derived from these metrics integrates uncertainties in the observational data sets. This multi-data set evaluation framework is a generic method that could be applied to the evaluation of other DGVM models. The results of the original ORC-VD published in 2005 for mid- to high-latitudes and of the new parameterization are evaluated against the above-described data sets. Significant improvements were found in the modeling of the distribution of tree functional types north of 40° N. Three additional sensitivity runs were carried out to separate the impact of different processes or drivers on simulated vegetation distribution, including soil freezing which limits net primary production through soil moisture availability in the root zone, elevated CO2 concentration since 1850, and the effects of frequency and severity of extreme cold events during the spin-up phase of the model.

The land and ocean absorb on average just over half of the anthropogenic emissions of carbon dioxide (CO2) every year. These CO2 "sinks" are modulated by climate change and variability. Here we use a suite of nine dynamic global vegetation models (DGVMs) and four ocean biogeochemical general circulation models (OBGCMs) to estimate trends driven by global and regional climate and atmospheric CO2 in land and oceanic CO2 exchanges with the atmosphere over the period 1990–2009, to attribute these trends to underlying processes in the models, and to quantify the uncertainty and level of inter-model agreement. The models were forced with reconstructed climate fields and observed global atmospheric CO2; land use and land cover changes are not included for the DGVMs. Over the period 1990–2009, the DGVMs simulate a mean global land carbon sink of −2.4 ± 0.7 Pg C yr−1 with a small significant trend of −0.06 ± 0.03 Pg C yr−2 (increasing sink). Over the more limited period 1990–2004, the ocean models simulate a mean ocean sink of −2.2 ± 0.2 Pg C yr−1 with a trend in the net C uptake that is indistinguishable from zero (−0.01 ± 0.02 Pg C yr−2). The two ocean models that extended the simulations until 2009 suggest a slightly stronger, but still small, trend of −0.02 ± 0.01 Pg C yr−2. Trends from land and ocean models compare favourably to the land greenness trends from remote sensing, atmospheric inversion results, and the residual land sink required to close the global carbon budget. Trends in the land sink are driven by increasing net primary production (NPP), whose statistically significant trend of 0.22 ± 0.08 Pg C yr−2 exceeds a significant trend in heterotrophic respiration of 0.16 ± 0.05 Pg C yr−2 – primarily as a consequence of widespread CO2 fertilisation of plant production. Most of the land-based trend in simulated net carbon uptake originates from natural ecosystems in the tropics (−0.04 ± 0.01 Pg C yr−2), with almost no trend over the northern land region, where recent warming and reduced rainfall offsets the positive impact of elevated atmospheric CO2 and changes in growing season length on carbon storage. The small uptake trend in the ocean models emerges because climate variability and change, and in particular increasing sea surface temperatures, tend to counter\-act the trend in ocean uptake driven by the increase in atmospheric CO2. Large uncertainty remains in the magnitude and sign of modelled carbon trends in several regions, as well as regarding the influence of land use and land cover changes on regional trends.

Ecosystems are important and dynamic components of the global carbon cycle, and terrestrial biospheric models (TBMs) are crucial tools in further understanding of how terrestrial carbon is stored and exchanged with the atmosphere across a variety of spatial and temporal scales. Improving TBM skills, and quantifying and reducing their estimation uncertainties, pose significant challenges. The Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal multi-scale and multi-model intercomparison effort set up to tackle these challenges. The MsTMIP protocol prescribes standardized environmental driver data that are shared among model teams to facilitate model–model and model–observation comparisons. This paper describes the global and North American environmental driver data sets prepared for the MsTMIP activity to both support their use in MsTMIP and make these data, along with the processes used in selecting/processing these data, accessible to a broader audience. Based on project needs and lessons learned from past model intercomparison activities, we compiled climate, atmospheric CO2 concentrations, nitrogen deposition, land use and land cover change (LULCC), C3 / C4 grasses fractions, major crops, phenology and soil data into a standard format for global (0.5° × 0.5° resolution) and regional (North American: 0.25° × 0.25° resolution) simulations. In order to meet the needs of MsTMIP, improvements were made to several of the original environmental data sets, by improving the quality, and/or changing their spatial and temporal coverage, and resolution. The resulting standardized model driver data sets are being used by over 20 different models participating in MsTMIP. The data are archived at the Oak Ridge National Laboratory Distributed Active Archive Center (ORNL DAAC, http://daac.ornl.gov) to provide long-term data management and distribution.

Fire is an important global ecological process that influences the distribution of biomes, with consequences for carbon, water, and energy budgets. Therefore it is impossible to appropriately model the history and future of the terrestrial ecosystems and the climate system without including fire. This study incorporates the process-based prognostic fire module SPITFIRE into the global vegetation model ORCHIDEE, which was then used to simulate burned area over the 20th century. Special attention was paid to the evaluation of other fire regime indicators such as seasonality, fire size and fire length, next to burned area. For 2001–2006, the simulated global spatial extent of fire agrees well with that given by satellite-derived burned area data sets (L3JRC, GLOBCARBON, GFED3.1), and 76–92% of the global burned area is simulated as collocated between the model and observation, depending on which data set is used for comparison. The simulated global mean annual burned area is 346 Mha yr−1, which falls within the range of 287–384 Mha yr−1 as given by the three observation data sets; and is close to the 344 Mha yr−1 by the GFED3.1 data when crop fires are excluded. The simulated long-term trend and variation of burned area agree best with the observation data in regions where fire is mainly driven by climate variation, such as boreal Russia (1930–2009), along with Canada and US Alaska (1950–2009). At the global scale, the simulated decadal fire variation over the 20th century is only in moderate agreement with the historical reconstruction, possibly because of the uncertainties of past estimates, and because land-use change fires and fire suppression are not explicitly included in the model. Over the globe, the size of large fires (the 95th quantile fire size) is underestimated by the model for the regions of high fire frequency, compared with fire patch data as reconstructed from MODIS 500 m burned area data. Two case studies of fire size distribution in Canada and US Alaska, and southern Africa indicate that both number and size of large fires are underestimated, which could be related with short fire patch length and low daily fire size. Future efforts should be directed towards building consistent spatial observation data sets for key parameters of the model in order to constrain the model error at each key step of the fire modelling.

Climate change is leading to a disproportionately large warming in the high northern latitudes, but the magnitude and sign of the future carbon balance of the Arctic are highly uncertain. Using 40 terrestrial biosphere models for the Alaskan Arctic from four recent model intercomparison projects – NACP (North American Carbon Program) site and regional syntheses, TRENDY (Trends in net land atmosphere carbon exchanges), and WETCHIMP (Wetland and Wetland CH4 Inter-comparison of Models Project) – we provide a baseline of terrestrial carbon cycle uncertainty, defined as the multi-model standard deviation (σ) for each quantity that follows. Mean annual absolute uncertainty was largest for soil carbon (14.0 ± 9.2 kg C m−2), then gross primary production (GPP) (0.22 ± 0.50 kg C m−2 yr−1), ecosystem respiration (Re) (0.23 ± 0.38 kg C m−2 yr−1), net primary production (NPP) (0.14 ± 0.33 kg C m−2 yr−1), autotrophic respiration (Ra) (0.09 ± 0.20 kg C m−2 yr−1), heterotrophic respiration (Rh) (0.14 ± 0.20 kg C m−2 yr−1), net ecosystem exchange (NEE) (−0.01 ± 0.19 kg C m−2 yr−1), and CH4 flux (2.52 ± 4.02 g CH4 m−2 yr−1). There were no consistent spatial patterns in the larger Alaskan Arctic and boreal regional carbon stocks and fluxes, with some models showing NEE for Alaska as a strong carbon sink, others as a strong carbon source, while still others as carbon neutral. Finally, AmeriFlux data are used at two sites in the Alaskan Arctic to evaluate the regional patterns; observed seasonal NEE was captured within multi-model uncertainty. This assessment of carbon cycle uncertainties may be used as a baseline for the improvement of experimental and modeling activities, as well as a reference for future trajectories in carbon cycling with climate change in the Alaskan Arctic and larger boreal region.

Agro-land surface models (agro-LSM) have been developed from the integration of specific crop processes into large-scale generic land surface models that allow calculating the spatial distribution and variability of energy, water and carbon fluxes within the soil–vegetation–atmosphere continuum. When developing agro-LSM models, particular attention must be given to the effects of crop phenology and management on the turbulent fluxes exchanged with the atmosphere, and the underlying water and carbon pools. A part of the uncertainty of agro-LSM models is related to their usually large number of parameters. In this study, we quantify the parameter-values uncertainty in the simulation of sugarcane biomass production with the agro-LSM ORCHIDEE–STICS, using a multi-regional approach with data from sites in Australia, La Réunion and Brazil. In ORCHIDEE–STICS, two models are chained: STICS, an agronomy model that calculates phenology and management, and ORCHIDEE, a land surface model that calculates biomass and other ecosystem variables forced by STICS phenology. First, the parameters that dominate the uncertainty of simulated biomass at harvest date are determined through a screening of 67 different parameters of both STICS and ORCHIDEE on a multi-site basis. Secondly, the uncertainty of harvested biomass attributable to those most sensitive parameters is quantified and specifically attributed to either STICS (phenology, management) or to ORCHIDEE (other ecosystem variables including biomass) through distinct Monte Carlo runs. The uncertainty on parameter values is constrained using observations by calibrating the model independently at seven sites. In a third step, a sensitivity analysis is carried out by varying the most sensitive parameters to investigate their effects at continental scale. A Monte Carlo sampling method associated with the calculation of partial ranked correlation coefficients is used to quantify the sensitivity of harvested biomass to input parameters on a continental scale across the large regions of intensive sugarcane cultivation in Australia and Brazil. The ten parameters driving most of the uncertainty in the ORCHIDEE–STICS modeled biomass at the 7 sites are identified by the screening procedure. We found that the 10 most sensitive parameters control phenology (maximum rate of increase of LAI) and root uptake of water and nitrogen (root profile and root growth rate, nitrogen stress threshold) in STICS, and photosynthesis (optimal temperature of photosynthesis, optimal carboxylation rate), radiation interception (extinction coefficient), and transpiration and respiration (stomatal conductance, growth and maintenance respiration coefficients) in ORCHIDEE. We find that the optimal carboxylation rate and photosynthesis temperature parameters contribute most to the uncertainty in harvested biomass simulations at site scale. The spatial variation of the ranked correlation between input parameters and modeled biomass at harvest is well explained by rain and temperature drivers, suggesting different climate-mediated sensitivities of modeled sugarcane yield to the model parameters, for Australia and Brazil. This study reveals the spatial and temporal patterns of uncertainty variability for a highly parameterized agro-LSM and calls for more systematic uncertainty analyses of such models.

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and a methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates, consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil-fuel combustion and cement production (EFF) are based on energy statistics, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated for the first time in this budget with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models forced by observed climate, CO2 and land cover change (some including nitrogen–carbon interactions). All uncertainties are reported as ±1σ, reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. For the last decade available (2003–2012), EFF was 8.6 ± 0.4 GtC yr−1, ELUC 0.9 ± 0.5 GtC yr−1, GATM 4.3 ± 0.1 GtC yr−1, SOCEAN 2.5 ± 0.5 GtC yr−1, and SLAND 2.8 ± 0.8 GtC yr−1. For year 2012 alone, EFF grew to 9.7 ± 0.5 GtC yr−1, 2.2% above 2011, reflecting a continued growing trend in these emissions, GATM was 5.1 ± 0.2 GtC yr−1, SOCEAN was 2.9 ± 0.5 GtC yr−1, and assuming an ELUC of 1.0 ± 0.5 GtC yr−1 (based on the 2001–2010 average), SLAND was 2.7 ± 0.9 GtC yr−1. GATM was high in 2012 compared to the 2003–2012 average, almost entirely reflecting the high EFF. The global atmospheric CO2 concentration reached 392.52 ± 0.10 ppm averaged over 2012. We estimate that EFF will increase by 2.1% (1.1–3.1%) to 9.9 ± 0.5 GtC in 2013, 61% above emissions in 1990, based on projections of world gross domestic product and recent changes in the carbon intensity of the economy. With this projection, cumulative emissions of CO2 will reach about 535 ± 55 GtC for 1870–2013, about 70% from EFF (390 ± 20 GtC) and 30% from ELUC (145 ± 50 GtC). This paper also documents any changes in the methods and data sets used in this new carbon budget from previous budgets (Le Quéré et al., 2013). All observations presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_2013_V2.3).

This paper describes a new model for the calculation of daily, high-resolution (up to 1 km) fire emissions, developed in the framework of the APIFLAME (Analysis and Prediction of the Impact of Fires on Air quality ModEling) project. The methodology relies on the classical approach, multiplying the burned area by the fuel load consumed and the emission factors specific to the vegetation burned. Emissions can be calculated on any user-specified domain, horizontal grid, and list of trace gases and aerosols, providing input information on the burned area (location, extent), and emission factors of the targeted species are available. The applicability to high spatial resolutions and the flexibility to different input data (including vegetation classifications) and domains are the main strength of the proposed algorithm. The modification of the default values and databases proposed does not require any change in the core of the model. The code may be used for the calculation of global or regional inventories. However, it has been developed and tested more specifically for Europe and the Mediterranean area. A regional analysis of fire activity and the resulting emissions in this region is provided. The burning season extends from June to October in most regions, with generally small but frequent fires in eastern Europe, western Russia, Ukraine and Turkey, and large events in the Mediterranean area. The resulting emissions represent a significant fraction of the total yearly emissions (on average amounting to ~ 30% of anthropogenic emissions for PM2.5, ~ 20% for CO). The uncertainty regarding the daily carbon emissions is estimated at ~ 100% based on an ensemble analysis. Considering the large uncertainties regarding emission factors, the potential error on the emissions for the various pollutants is even larger. Comparisons with other widely used emission inventories show good correlations but discrepancies of a factor of 2–4 in the amplitude of the emissions, our results being generally on the higher end.

Stand-replacing fires are the dominant fire type in North American boreal forests. They leave a historical legacy of a mosaic landscape of different aged forest cohorts. This forest age dynamics must be included in vegetation models to accurately quantify the role of fire in the historical and current regional forest carbon balance. The present study adapted the global process-based vegetation model ORCHIDEE to simulate the CO2 emissions from boreal forest fire and the subsequent recovery after a stand-replacing fire; the model represents postfire new cohort establishment, forest stand structure and the self-thinning process. Simulation results are evaluated against observations of three clusters of postfire forest chronosequences in Canada and Alaska. The variables evaluated include: fire carbon emissions, CO2 fluxes (gross primary production, total ecosystem respiration and net ecosystem exchange), leaf area index, and biometric measurements (aboveground biomass carbon, forest floor carbon, woody debris carbon, stand individual density, stand basal area, and mean diameter at breast height). When forced by local climate and the atmospheric CO2 history at each chronosequence site, the model simulations generally match the observed CO2 fluxes and carbon stock data well, with model-measurement mean square root of deviation comparable with the measurement accuracy (for CO2 flux ~100 g C m−2 yr−1, for biomass carbon ~1000 g C m−2 and for soil carbon ~2000 g C m−2). We find that the current postfire forest carbon sink at the evaluation sites, as observed by chronosequence methods, is mainly due to a combination of historical CO2 increase and forest succession. Climate change and variability during this period offsets some of these expected carbon gains. The negative impacts of climate were a likely consequence of increasing water stress caused by significant temperature increases that were not matched by concurrent increases in precipitation. Our simulation results demonstrate that a global vegetation model such as ORCHIDEE is able to capture the essential ecosystem processes in fire-disturbed boreal forests and produces satisfactory results in terms of both carbon fluxes and carbon-stock evolution after fire. This makes the model suitable for regional simulations in boreal regions where fire regimes play a key role in the ecosystem carbon balance.

This study describes how management of grasslands is included in the Organizing Carbon and Hydrology in Dynamic Ecosystems (ORCHIDEE) process-based ecosystem model designed for large-scale applications, and how management affects modeled grassland–atmosphere CO2 fluxes. The new model, ORCHIDEE-GM (grassland management) is enabled with a management module inspired from a grassland model (PaSim, version 5.0), with two grassland management practices being considered, cutting and grazing. The evaluation of the results from ORCHIDEE compared with those of ORCHIDEE-GM at 11 European sites, equipped with eddy covariance and biometric measurements, shows that ORCHIDEE-GM can realistically capture the cut-induced seasonal variation in biometric variables (LAI: leaf area index; AGB: aboveground biomass) and in CO2 fluxes (GPP: gross primary productivity; TER: total ecosystem respiration; and NEE: net ecosystem exchange). However, improvements at grazing sites are only marginal in ORCHIDEE-GM due to the difficulty in accounting for continuous grazing disturbance and its induced complex animal–vegetation interactions. Both NEE and GPP on monthly to annual timescales can be better simulated in ORCHIDEE-GM than in ORCHIDEE without management. For annual CO2 fluxes, the NEE bias and RMSE (root mean square error) in ORCHIDEE-GM are reduced by 53% and 20%, respectively, compared to ORCHIDEE. ORCHIDEE-GM is capable of modeling the net carbon balance (NBP) of managed temperate grasslands (37 ± 30 gC m−2 yr−1 (P < 0.01) over the 11 sites) because the management module contains provisions to simulate the carbon fluxes of forage yield, herbage consumption, animal respiration and methane emissions.

Accurate assessments of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO2 emissions from fossil fuel combustion and cement production (EFF) are based on energy statistics, while emissions from Land-Use Change (ELUC), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms. For the last decade available (2002–2011), EFF was 8.3 &pm; 0.4 PgC yr−1, ELUC 1.0 &pm; 0.5 PgC yr−1, GATM 4.3 &pm; 0.1PgC yr−1, SOCEAN 2.5 &pm; 0.5 PgC yr−1, and SLAND 2.6 &pm; 0.8 PgC yr−1. For year 2011 alone, EFF was 9.5 &pm; 0.5 PgC yr−1, 3.0 percent above 2010, reflecting a continued trend in these emissions; ELUC was 0.9 &pm; 0.5 PgC yr−1, approximately constant throughout the decade; GATM was 3.6 &pm; 0.2 PgC yr−1, SOCEAN was 2.7 &pm; 0.5 PgC yr−1, and SLAND was 4.1 &pm; 0.9 PgC yr−1. GATM was low in 2011 compared to the 2002–2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Niña conditions in the Pacific Ocean. The global atmospheric CO2 concentration reached 391.31 &pm; 0.13 ppm at the end of year 2011. We estimate that EFF will have increased by 2.6% (1.9–3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as &pm;1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future. All data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_V2013).

This REgional Carbon Cycle Assessment and Processes regional study provides a synthesis of the carbon balance of terrestrial ecosystems in East Asia, a region comprised of China, Japan, North and South Korea, and Mongolia. We estimate the current terrestrial carbon balance of East Asia and its driving mechanisms during 1990–2009 using three different approaches: inventories combined with satellite greenness measurements, terrestrial ecosystem carbon cycle models and atmospheric inversion models. The magnitudes of East Asia's terrestrial carbon sink from these three approaches are comparable: −0.293±0.033 PgC yr−1 from inventory–remote sensing model–data fusion approach, −0.413±0.141 PgC yr−1 (not considering biofuel emissions) or −0.224±0.141 PgC yr−1 (considering biofuel emissions) for carbon cycle models, and −0.270±0.507 PgC yr−1 for atmospheric inverse models. Here and in the following, the numbers behind ± signs are standard deviations. The ensemble of ecosystem modeling based analyses further suggests that at the regional scale, climate change and rising atmospheric CO2 together resulted in a carbon sink of −0.289±0.135 PgC yr−1, while land-use change and nitrogen deposition had a contribution of −0.013±0.029 PgC yr−1 and −0.107±0.025 PgC yr−1, respectively. Although the magnitude of climate change effects on the carbon balance varies among different models, all models agree that in response to climate change alone, southern China experienced an increase in carbon storage from 1990 to 2009, while northern East Asia including Mongolia and north China showed a decrease in carbon storage. Overall, our results suggest that about 13–27% of East Asia's CO2 emissions from fossil fuel burning have been offset by carbon accumulation in its terrestrial territory over the period from 1990 to 2009. The underlying mechanisms of carbon sink over East Asia still remain largely uncertain, given the diversity and intensity of land management processes, and the regional conjunction of many drivers such as nutrient deposition, climate, atmospheric pollution and CO2 changes, which cannot be considered as independent for their effects on carbon storage.

We analyze how biases of meteorological drivers impact the calculation of ecosystem CO2, water and energy fluxes by models. To do so, we drive the same ecosystem model by meteorology from gridded products and by meteorology from local observation at eddy-covariance flux sites. The study is focused on six flux tower sites in France spanning across a climate gradient of 7–14 °C annual mean surface air temperature and 600–1040 mm mean annual rainfall, with forest, grassland and cropland ecosystems. We evaluate the results of the ORCHIDEE process-based model driven by meteorology from four different analysis data sets against the same model driven by site-observed meteorology. The evaluation is decomposed into characteristic time scales. The main result is that there are significant differences in meteorology between analysis data sets and local observation. The phase of seasonal cycle of air temperature, humidity and shortwave downward radiation is reproduced correctly by all meteorological models (average R2 = 0.90). At sites located in altitude, the misfit of meteorological drivers from analysis data sets and tower meteorology is the largest. We show that day-to-day variations in weather are not completely well reproduced by meteorological models, with R2 between analysis data sets and measured local meteorology going from 0.35 to 0.70. The bias of meteorological driver impacts the flux simulation by ORCHIDEE, and thus would have an effect on regional and global budgets. The forcing error, defined by the simulated flux difference resulting from prescribing modeled instead of observed local meteorology drivers to ORCHIDEE, is quantified for the six studied sites at different time scales. The magnitude of this forcing error is compared to that of the model error defined as the modeled-minus-observed flux, thus containing uncertain parameterizations, parameter values, and initialization. The forcing error is on average smaller than but still comparable to model error, with the ratio of forcing error to model error being the largest on daily time scale (86%) and annual time scales (80%). The forcing error incurred from using a gridded meteorological data set to drive vegetation models is therefore an important component of the uncertainty budget of regional CO2, water and energy fluxes simulations, and should be taken into consideration in up-scaling studies.

Atmospheric CO2 drives most of the greenhouse effect increase. One major uncertainty on the future rate of increase of CO2 in the atmosphere is the impact of the anticipated climate change on the vegetation. Dynamic Global Vegetation Models (DGVM) are used to address this question. ORCHIDEE is such a DGVM that has proven useful for climate change studies. However, there is no objective and methodological way to accurately assess each new available version on the global scale. In this paper, we submit a methodological evaluation of ORCHIDEE by correlating satellite-derived Vegetation Index time series against those of the modeled Fraction of absorbed Photosynthetically Active Radiation (FPAR). A perfect correlation between the two is not expected, however an improvement of the model should lead to an increase of the overall performance. We detail two case studies in which model improvements are demonstrated, using our methodology. In the first one, a new phenology version in ORCHIDEE is shown to bring a significant impact on the simulated annual cycles, in particular for C3 Grasses and C3 Crops. In the second case study, we compare the simulations when using two different weather fields to drive ORCHIDEE. The ERA-Interim forcing leads to a better description of the FPAR interannual anomalies than the simulation forced by a mixed CRU-NCEP dataset. This work shows that long time series of satellite observations, despite their uncertainties, can identify weaknesses in global vegetation models, a necessary first step to improving them.

This paper is a modelling study of crop management impacts on carbon and water fluxes at a range of European sites. The model is a crop growth model (STICS) coupled with a process-based land surface model (ORCHIDEE). The data are online eddy-covariance observations of CO2 and H2O fluxes at five European maize cultivation sites. The results show that the ORCHIDEE-STICS model explains up to 75 % of the observed daily net CO2 ecosystem exchange (NEE) variance, and up to 79 % of the latent heat flux (LE) variance at five sites. The model is better able to reproduce gross primary production (GPP) variations than terrestrial ecosystem respiration (TER) variations. We conclude that structural deficiencies in the model parameterizations of leaf area index (LAI) and TER are the main sources of error in simulating CO2 and H2O fluxes. A number of sensitivity tests, with variable crop variety, nitrogen fertilization, irrigation, and planting date, indicate that any of these management factors is able to change NEE by more than 15 %, but that the response of NEE to management parameters is highly site-dependent. Changes in management parameters are found to impact not only the daily values of NEE and LE, but also the cumulative yearly values. In addition, LE is shown to be less sensitive to management parameters than NEE. Multi-site model evaluations, coupled with sensitivity analysis to management parameters, thus provide important information about model errors, which helps to improve the simulation of CO2 and H2O fluxes across European croplands.

Dynamic Vegetation Models (DVMs) simulate energy, water and carbon fluxes between the ecosystem and the atmosphere, between the vegetation and the soil, and between plant organs. They also estimate the potential biomass of a forest in equilibrium having grown under a given climate and atmospheric CO2 level. In this study, we evaluate the Above Ground Woody Biomass (AGWB) and the above ground woody Net Primary Productivity (NPPAGW) simulated by the DVM ORCHIDEE across Amazonian forests, by comparing the simulation results to a large set of ground measurements (220 sites for biomass, 104 sites for NPPAGW). We found that the NPPAGW is on average overestimated by 63%. We also found that the fraction of biomass that is lost through mortality is 85% too high. These model biases nearly compensate each other to give an average simulated AGWB close to the ground measurement average. Nevertheless, the simulated AGWB spatial distribution differs significantly from the observations. Then, we analyse the discrepancies in biomass with regards to discrepancies in NPPAGW and those in the rate of mortality. When we correct for the error in NPPAGW, the errors on the spatial variations in AGWB are exacerbated, showing clearly that a large part of the misrepresentation of biomass comes from a wrong modelling of mortality processes. Previous studies showed that Amazonian forests with high productivity have a higher mortality rate than forests with lower productivity. We introduce this relationship, which results in strongly improved modelling of biomass and of its spatial variations. We discuss the possibility of modifying the mortality modelling in ORCHIDEE, and the opportunity to improve forest productivity modelling through the integration of biomass measurements, in particular from remote sensing.

European ecosystems are thought to take up large amounts of carbon, but neither the rate nor the contributions of the underlying processes are well known. In the second half of the 20th century, carbon dioxide concentrations have risen by more that 100 ppm, atmospheric nitrogen deposition has more than doubled, and European mean temperatures were increasing by 0.02 °C yr−1. The extents of forest and grasslands have increased with the respective rates of 5800 km2 yr−1 and 1100 km2 yr−1 as agricultural land has been abandoned at a rate of 7000 km2 yr−1. In this study, we analyze the responses of European land ecosystems to the aforementioned environmental changes using results from four process-based ecosystem models: BIOME-BGC, JULES, ORCHIDEE, and O-CN. The models suggest that European ecosystems sequester carbon at a rate of 56 TgC yr−1 (mean of four models for 1951–2000) with strong interannual variability (±88 TgC yr−1, average across models) and substantial inter-model uncertainty (±39 TgC yr−1). Decadal budgets suggest that there has been a continuous increase in the mean net carbon storage of ecosystems from 85 TgC yr−1 in 1980s to 108 TgC yr−1 in 1990s, and to 114 TgC yr−1 in 2000–2007. The physiological effect of rising CO2 in combination with nitrogen deposition and forest re-growth have been identified as the important explanatory factors for this net carbon storage. Changes in the growth of woody vegetation are suggested as an important contributor to the European carbon sink. Simulated ecosystem responses were more consistent for the two models accounting for terrestrial carbon-nitrogen dynamics than for the two models which only accounted for carbon cycling and the effects of land cover change. Studies of the interactions of carbon-nitrogen dynamics with land use changes are needed to further improve the quantitative understanding of the driving forces of the European land carbon balance.

We are comparing spatially explicit process-model based estimates of the terrestrial carbon balance and its components over Africa and confront them with remote sensing based proxies of vegetation productivity and atmospheric inversions of land-atmosphere net carbon exchange. Particular emphasis is on characterizing the patterns of interannual variability of carbon fluxes and analyzing the factors and processes responsible for it. For this purpose simulations with the terrestrial biosphere models ORCHIDEE, LPJ-DGVM, LPJ-Guess and JULES have been performed using a standardized modeling protocol and a uniform set of corrected climate forcing data. While the models differ concerning the absolute magnitude of carbon fluxes, we find several robust patterns of interannual variability among the models. Models exhibit largest interannual variability in southern and eastern Africa, regions which are primarily covered by herbaceous vegetation. Interannual variability of the net carbon balance appears to be more strongly influenced by gross primary production than by ecosystem respiration. A principal component analysis indicates that moisture is the main driving factor of interannual gross primary production variability for those regions. On the contrary in a large part of the inner tropics radiation appears to be limiting in two models. These patterns are partly corroborated by remotely sensed vegetation properties from the SeaWiFS satellite sensor. Inverse atmospheric modeling estimates of surface carbon fluxes are less conclusive at this point, implying the need for a denser network of observation stations over Africa.

Globally, the year 2003 is associated with one of the largest atmospheric CO2 rises on record. In the same year, Europe experienced an anomalously strong flux of CO2 from the land to the atmosphere associated with an exceptionally dry and hot summer in Western and Central Europe. In this study we analyze the magnitude of this carbon flux anomaly and key driving ecosystem processes using simulations of seven terrestrial ecosystem models of different complexity and types (process-oriented and diagnostic). We address the following questions: (1) how large were deviations in the net European carbon flux in 2003 relative to a short-term baseline (1998–2002) and to longer-term variations in annual fluxes (1980 to 2005), (2) which European regions exhibited the largest changes in carbon fluxes during the growing season 2003, and (3) which ecosystem processes controlled the carbon balance anomaly . In most models the prominence of 2003 anomaly in carbon fluxes declined with lengthening of the reference period from one year to 16 years. The 2003 anomaly for annual net carbon fluxes ranged between 0.35 and –0.63 Pg C for a reference period of one year and between 0.17 and –0.37 Pg C for a reference period of 16 years for the whole Europe. In Western and Central Europe, the anomaly in simulated net ecosystem productivity (NEP) over the growing season in 2003 was outside the 1σ variance bound of the carbon flux anomalies for 1980–2005 in all models. The estimated anomaly in net carbon flux ranged between –42 and –158 Tg C for Western Europe and between 24 and –129 Tg C for Central Europe depending on the model used. All models responded to a dipole pattern of the climate anomaly in 2003. In Western and Central Europe NEP was reduced due to heat and drought. In contrast, lower than normal temperatures and higher air humidity decreased NEP over Northeastern Europe. While models agree on the sign of changes in simulated NEP and gross primary productivity in 2003 over Western and Central Europe, models diverge in the estimates of anomalies in ecosystem respiration. Except for two process models which simulate respiration increase, most models simulated a decrease in ecosystem respiration in 2003. The diagnostic models showed a weaker decrease in ecosystem respiration than the process-oriented models. Based on the multi-model simulations we estimated the total carbon flux anomaly over the 2003 growing season in Europe to range between –0.02 and –0.27 Pg C relative to the net carbon flux in 1998–2002.

Three terrestrial biosphere models (LPJ, Orchidee, Biome-BGC) were evaluated with respect to their ability to simulate large-scale climate related trends in gross primary production (GPP) across European forests. Simulated GPP and leaf area index (LAI) were compared with GPP estimates based on flux separated eddy covariance measurements of net ecosystem exchange and LAI measurements along a temperature gradient ranging from the boreal to the Mediterranean region. The three models capture qualitatively the pattern suggested by the site data: an increase in GPP from boreal to temperate and a subsequent decline from temperate to Mediterranean climates. The models consistently predict higher GPP for boreal and lower GPP for Mediterranean forests. Based on a decomposition of GPP into absorbed photosynthetic active radiation (APAR) and radiation use efficiency (RUE), the overestimation of GPP for the boreal coniferous forests appears to be primarily related to too high simulated LAI - and thus light absorption (APAR) – rather than too high radiation use efficiency. We cannot attribute the tendency of the models to underestimate GPP in the water limited region to model structural deficiencies with confidence. A likely dry bias of the input meteorological data in southern Europe may create this pattern. On average, the models compare similarly well to the site GPP data (RMSE of ~30% or 420 gC/m2/yr) but differences are apparent for different ecosystem types. In terms of absolute values, we find the agreement between site based GPP estimates and simulations acceptable when we consider uncertainties about the accuracy in model drivers, a potential representation bias of the eddy covariance sites, and uncertainties related to the method of deriving GPP from eddy covariance measurements data. Continental to global data-model comparison studies should be fostered in the future since they are necessary to identify consistent model bias along environmental gradients.