Search

The importance of phosphorus (P) in regulating ecosystem responses to climate change has fostered P-cycle implementation in land surface models, but their CO 2 effects predictions have not been evaluated against measurements. Here, we perform a data-driven model evaluation where simulations of eight widely used P-enabled models were confronted with observations from a long-term free-air CO 2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO 2 effect on ecosystem carbon (C) sequestration, but they generally overestimated the effects on plant C uptake and growth. We identify leaf-to-canopy scaling of photosynthesis, plant tissue stoichiometry, plant belowground C allocation, and the subsequent consequences for plant-microbial interaction as key areas in which models of ecosystem C-P interaction can be improved. Together, this data-model intercomparison reveals data-driven insights into the performance and functionality of P-enabled models and adds to the existing evidence that the global CO 2 -driven carbon sink is overestimated by models.

The leaf economics spectrum 1,2 and the global spectrum of plant forms and functions 3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species 2 . Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities 4 . However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability 4,5 . Here we derive a set of ecosystem functions 6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems 7,8 ., (© 2021. The Author(s).)

In 2018, central and northern parts of Europe experienced heat and drought conditions over many months from spring to autumn, strongly affecting both natural ecosystems and crops. Besides their impact on nature and society, events like this can be used to study the impact of climate variations on the terrestrial carbon cycle, which is an important determinant of the future climate trajectory. Here, variations in the regional net ecosystem exchange (NEE) of CO 2 between terrestrial ecosystems and the atmosphere were quantified from measurements of atmospheric CO 2 mole fractions. Over Europe, several observational records have been maintained since at least 1999, giving us the opportunity to assess the 2018 anomaly in the context of at least two decades of variations, including the strong climate anomaly in 2003. In addition to an atmospheric inversion with temporally explicitly estimated anomalies, we use an inversion based on empirical statistical relations between anomalies in the local NEE and anomalies in local climate conditions. For our analysis period 1999-2018, we find that higher-than-usual NEE in hot and dry summers may tend to arise in Central Europe from enhanced ecosystem respiration due to the elevated temperatures, and in Southern Europe from reduced photosynthesis due to the reduced water availability. Despite concerns in the literature, the level of agreement between regression-based NEE anomalies and temporally explicitly estimated anomalies indicates that the atmospheric CO 2 measurements from the relatively dense European station network do provide information about the year-to-year variations of Europe's carbon sources and sinks, at least in summer. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'.

Interannual variations in the large-scale net ecosystem exchange (NEE) of CO 2 between the terrestrial biosphere and the atmosphere were estimated for 1957-2017 from sustained measurements of atmospheric CO 2 mixing ratios. As the observations are sparse in the early decades, available records were combined into a 'quasi-homogeneous' dataset based on similarity in their signals, to minimize spurious variations from beginning or ending data records. During El Niño events, CO 2 is anomalously released from the tropical band, and a few months later also in the northern extratropical band. This behaviour can approximately be represented by a linear relationship of the NEE anomalies and local air temperature anomalies, with sensitivity coefficients depending on geographical location and season. The apparent climate sensitivity of global total NEE against variations in pan-tropically averaged annual air temperature slowly changed over time during the 1957-2017 period, first increasing (though less strongly than in previous studies) but then decreasing again. However, only part of this change can be attributed to actual changes in local physiological or ecosystem processes, the rest probably arising from shifts in the geographical area of dominating temperature variations.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'., (© 2018 The Authors.)

Evaluating the response of the land carbon sink to the anomalies in temperature and drought imposed by El Niño events provides insights into the present-day carbon cycle and its climate-driven variability. It is also a necessary step to build confidence in terrestrial ecosystems models' response to the warming and drying stresses expected in the future over many continents, and particularly in the tropics. Here we present an in-depth analysis of the response of the terrestrial carbon cycle to the 2015/2016 El Niño that imposed extreme warming and dry conditions in the tropics and other sensitive regions. First, we provide a synthesis of the spatio-temporal evolution of anomalies in net land-atmosphere CO 2 fluxes estimated by two in situ measurements based on atmospheric inversions and 16 land-surface models (LSMs) from TRENDYv6. Simulated changes in ecosystem productivity, decomposition rates and fire emissions are also investigated. Inversions and LSMs generally agree on the decrease and subsequent recovery of the land sink in response to the onset, peak and demise of El Niño conditions and point to the decreased strength of the land carbon sink: by 0.4-0.7 PgC yr -1 (inversions) and by 1.0 PgC yr -1 (LSMs) during 2015/2016. LSM simulations indicate that a decrease in productivity, rather than increase in respiration, dominated the net biome productivity anomalies in response to ENSO throughout the tropics, mainly associated with prolonged drought conditions.This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'., (© 2018 The Authors.)

Determining whether the terrestrial biosphere will be a source or sink of carbon (C) under a future climate of elevated CO 2 (eCO 2 ) and warming requires accurate quantification of gross primary production (GPP), the largest flux of C in the global C cycle. We evaluated 6 years (2007-2012) of flux-derived GPP data from the Prairie Heating and CO 2 Enrichment (PHACE) experiment, situated in a grassland in Wyoming, USA. The GPP data were used to calibrate a light response model whose basic formulation has been successfully used in a variety of ecosystems. The model was extended by modeling maximum photosynthetic rate (A max ) and light-use efficiency (Q) as functions of soil water, air temperature, vapor pressure deficit, vegetation greenness, and nitrogen at current and antecedent (past) timescales. The model fits the observed GPP well (R 2  = 0.79), which was confirmed by other model performance checks that compared different variants of the model (e.g. with and without antecedent effects). Stimulation of cumulative 6-year GPP by warming (29%, P = 0.02) and eCO 2 (26%, P = 0.07) was primarily driven by enhanced C uptake during spring (129%, P = 0.001) and fall (124%, P = 0.001), respectively, which was consistent across years. Antecedent air temperature (Tair ant ) and vapor pressure deficit (VPD ant ) effects on A max (over the past 3-4 days and 1-3 days, respectively) were the most significant predictors of temporal variability in GPP among most treatments. The importance of VPD ant suggests that atmospheric drought is important for predicting GPP under current and future climate; we highlight the need for experimental studies to identify the mechanisms underlying such antecedent effects. Finally, posterior estimates of cumulative GPP under control and eCO 2 treatments were tested as a benchmark against 12 terrestrial biosphere models (TBMs). The narrow uncertainties of these data-driven GPP estimates suggest that they could be useful semi-independent data streams for validating TBMs., (© 2017 John Wiley & Sons Ltd.)

Large interannual variations in the measured growth rate of atmospheric carbon dioxide (CO 2 ) originate primarily from fluctuations in carbon uptake by land ecosystems. It remains uncertain, however, to what extent temperature and water availability control the carbon balance of land ecosystems across spatial and temporal scales. Here we use empirical models based on eddy covariance data and process-based models to investigate the effect of changes in temperature and water availability on gross primary productivity (GPP), terrestrial ecosystem respiration (TER) and net ecosystem exchange (NEE) at local and global scales. We find that water availability is the dominant driver of the local interannual variability in GPP and TER. To a lesser extent this is true also for NEE at the local scale, but when integrated globally, temporal NEE variability is mostly driven by temperature fluctuations. We suggest that this apparent paradox can be explained by two compensatory water effects. Temporal water-driven GPP and TER variations compensate locally, dampening water-driven NEE variability. Spatial water availability anomalies also compensate, leaving a dominant temperature signal in the year-to-year fluctuations of the land carbon sink. These findings help to reconcile seemingly contradictory reports regarding the importance of temperature and water in controlling the interannual variability of the terrestrial carbon balance. Our study indicates that spatial climate covariation drives the global carbon cycle response.

Atmospheric monitoring of high northern latitudes (above 40°N) has shown an enhanced seasonal cycle of carbon dioxide (CO2) since the 1960s, but the underlying mechanisms are not yet fully understood. The much stronger increase in high latitudes relative to low ones suggests that northern ecosystems are experiencing large changes in vegetation and carbon cycle dynamics. We found that the latitudinal gradient of the increasing CO2 amplitude is mainly driven by positive trends in photosynthetic carbon uptake caused by recent climate change and mediated by changing vegetation cover in northern ecosystems. Our results underscore the importance of climate-vegetation-carbon cycle feedbacks at high latitudes; moreover, they indicate that in recent decades, photosynthetic carbon uptake has reacted much more strongly to warming than have carbon release processes., (Copyright © 2016, American Association for the Advancement of Science.)

The response of the forest carbon (C) balance to changes in nitrogen (N) deposition is uncertain, partly owing to diverging representations of N cycle processes in dynamic global vegetation models (DGVMs). Here, we examined how different assumptions about the degree of flexibility of the ecosystem's C : N ratios contribute to this uncertainty, and which of these assumptions best correspond to the available data. We applied these assumptions within the framework of a DGVM and compared the results to responses in net primary productivity (NPP), leaf N concentration, and ecosystem N partitioning, observed at 22 forest N fertilization experiments. Employing flexible ecosystem pool C : N ratios generally resulted in the most convincing model-data agreement with respect to production and foliar N responses. An intermediate degree of stoichiometric flexibility in vegetation, where wood C : N ratio changes were decoupled from leaf and root C : N ratio changes, led to consistent simulation of production and N cycle responses to N addition. Assuming fixed C : N ratios or scaling leaf N concentration changes to other tissues, commonly assumed by DGVMs, was not supported by reported data. Between the tested assumptions, the simulated changes in ecosystem C storage relative to changes in C assimilation varied by up to 20%., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

The growth rate of atmospheric carbon dioxide (CO2) concentrations since industrialization is characterized by large interannual variability, mostly resulting from variability in CO2 uptake by terrestrial ecosystems (typically termed carbon sink). However, the contributions of regional ecosystems to that variability are not well known. Using an ensemble of ecosystem and land-surface models and an empirical observation-based product of global gross primary production, we show that the mean sink, trend, and interannual variability in CO2 uptake by terrestrial ecosystems are dominated by distinct biogeographic regions. Whereas the mean sink is dominated by highly productive lands (mainly tropical forests), the trend and interannual variability of the sink are dominated by semi-arid ecosystems whose carbon balance is strongly associated with circulation-driven variations in both precipitation and temperature., (Copyright © 2015, American Association for the Advancement of Science.)

We analysed the responses of 11 ecosystem models to elevated atmospheric [CO2 ] (eCO2 ) at two temperate forest ecosystems (Duke and Oak Ridge National Laboratory (ORNL) Free-Air CO2 Enrichment (FACE) experiments) to test alternative representations of carbon (C)-nitrogen (N) cycle processes. We decomposed the model responses into component processes affecting the response to eCO2 and confronted these with observations from the FACE experiments. Most of the models reproduced the observed initial enhancement of net primary production (NPP) at both sites, but none was able to simulate both the sustained 10-yr enhancement at Duke and the declining response at ORNL: models generally showed signs of progressive N limitation as a result of lower than observed plant N uptake. Nonetheless, many models showed qualitative agreement with observed component processes. The results suggest that improved representation of above-ground-below-ground interactions and better constraints on plant stoichiometry are important for a predictive understanding of eCO2 effects. Improved accuracy of soil organic matter inventories is pivotal to reduce uncertainty in the observed C-N budgets. The two FACE experiments are insufficient to fully constrain terrestrial responses to eCO2 , given the complexity of factors leading to the observed diverging trends, and the consequential inability of the models to explain these trends. Nevertheless, the ecosystem models were able to capture important features of the experiments, lending some support to their projections., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

Projections of future changes in land carbon (C) storage using biogeochemical models depend on accurately modeling the interactions between the C and nitrogen (N) cycles. Here, we present a framework for analyzing N limitation in global biogeochemical models to explore how C-N interactions of current models compare to field observations, identify the processes causing model divergence, and identify future observation and experiment needs. We used a set of N-fertilization simulations from two global biogeochemical models (CLM-CN and O-CN) that use different approaches to modeling C-N interactions. On the global scale, net primary productivity (NPP) in the CLM-CN model was substantially more responsive to N fertilization than in the O-CN model. The most striking difference between the two models occurred for humid tropical forests, where the CLM-CN simulated a 62% increase in NPP at high N addition levels (30 g N m(-2) yr(-1)), while the O-CN predicted a 2% decrease in NPP due to N fertilization increasing plant respiration more than photosynthesis. Across 35 temperate and boreal forest sites with field N-fertilization experiments, we show that the CLM-CN simulated a 46% increase in aboveground NPP in response to N, which exceeded the observed increase of 25%. In contrast, the O-CN only simulated a 6% increase in aboveground NPP at the N-fertilization sites. Despite the small response of NPP to N fertilization, the O-CN model accurately simulated ecosystem retention of N and the fate of added N to vegetation when compared to empirical (15) N tracer application studies. In contrast, the CLM-CN predicted lower total ecosystem N retention and partitioned more losses to volatilization than estimated from observed N budgets of small catchments. These results point to the need for model improvements in both models in order to enhance the accuracy with which global C-N cycle feedbacks are simulated., (© 2013 John Wiley & Sons Ltd.)

The purpose of this study was to evaluate 10 process-based terrestrial biosphere models that were used for the IPCC fifth Assessment Report. The simulated gross primary productivity (GPP) is compared with flux-tower-based estimates by Jung et al. [Journal of Geophysical Research 116 (2011) G00J07] (JU11). The net primary productivity (NPP) apparent sensitivity to climate variability and atmospheric CO2 trends is diagnosed from each model output, using statistical functions. The temperature sensitivity is compared against ecosystem field warming experiments results. The CO2 sensitivity of NPP is compared to the results from four Free-Air CO2 Enrichment (FACE) experiments. The simulated global net biome productivity (NBP) is compared with the residual land sink (RLS) of the global carbon budget from Friedlingstein et al. [Nature Geoscience 3 (2010) 811] (FR10). We found that models produce a higher GPP (133 ± 15 Pg C yr(-1) ) than JU11 (118 ± 6 Pg C yr(-1) ). In response to rising atmospheric CO2 concentration, modeled NPP increases on average by 16% (5-20%) per 100 ppm, a slightly larger apparent sensitivity of NPP to CO2 than that measured at the FACE experiment locations (13% per 100 ppm). Global NBP differs markedly among individual models, although the mean value of 2.0 ± 0.8 Pg C yr(-1) is remarkably close to the mean value of RLS (2.1 ± 1.2 Pg C yr(-1) ). The interannual variability in modeled NBP is significantly correlated with that of RLS for the period 1980-2009. Both model-to-model and interannual variation in model GPP is larger than that in model NBP due to the strong coupling causing a positive correlation between ecosystem respiration and GPP in the model. The average linear regression slope of global NBP vs. temperature across the 10 models is -3.0 ± 1.5 Pg C yr(-1) °C(-1) , within the uncertainty of what derived from RLS (-3.9 ± 1.1 Pg C yr(-1) °C(-1) ). However, 9 of 10 models overestimate the regression slope of NBP vs. precipitation, compared with the slope of the observed RLS vs. precipitation. With most models lacking processes that control GPP and NBP in addition to CO2 and climate, the agreement between modeled and observation-based GPP and NBP can be fortuitous. Carbon-nitrogen interactions (only separable in one model) significantly influence the simulated response of carbon cycle to temperature and atmospheric CO2 concentration, suggesting that nutrients limitations should be included in the next generation of terrestrial biosphere models., (© 2013 Blackwell Publishing Ltd.)

Interactions between the terrestrial nitrogen (N) and carbon (C) cycles shape the response of ecosystems to global change. However, the global distribution of nitrogen availability and its importance in global biogeochemistry and biogeochemical interactions with the climate system remain uncertain. Based on projections of a terrestrial biosphere model scaling ecological understanding of nitrogen-carbon cycle interactions to global scales, anthropogenic nitrogen additions since 1860 are estimated to have enriched the terrestrial biosphere by 1.3 Pg N, supporting the sequestration of 11.2 Pg C. Over the same time period, CO2 fertilization has increased terrestrial carbon storage by 134.0 Pg C, increasing the terrestrial nitrogen stock by 1.2 Pg N. In 2001-2010, terrestrial ecosystems sequestered an estimated total of 27 Tg N yr(-1) (1.9 Pg C yr(-1)), of which 10 Tg N yr(-1) (0.2 Pg C yr(-1)) are due to anthropogenic nitrogen deposition. Nitrogen availability already limits terrestrial carbon sequestration in the boreal and temperate zone, and will constrain future carbon sequestration in response to CO2 fertilization (regionally by up to 70% compared with an estimate without considering nitrogen-carbon interactions). This reduced terrestrial carbon uptake will probably dominate the role of the terrestrial nitrogen cycle in the climate system, as it accelerates the accumulation of anthropogenic CO2 in the atmosphere. However, increases of N2O emissions owing to anthropogenic nitrogen and climate change (at a rate of approx. 0.5 Tg N yr(-1) per 1°C degree climate warming) will add an important long-term climate forcing.

Changes in temperature and carbon dioxide during glacial cycles recorded in Antarctic ice cores are tightly coupled. However, this relationship does not hold for interglacials. While climate cooled towards the end of both the last (Eemian) and present (Holocene) interglacials, CO2 remained stable during the Eemian while rising in the Holocene. We identify and review twelve biogeochemical mechanisms of terrestrial (vegetation dynamics and CO2 fertilization, land use, wildfire, accumulation of peat, changes in permafrost carbon, subaerial volcanic outgassing) and marine origin (changes in sea surface temperature, carbonate compensation to deglaciation and terrestrial biosphere regrowth, shallow-water carbonate sedimentation, changes in the soft tissue pump, and methane hydrates), which potentially may have contributed to the CO2 dynamics during interglacials but which remain not well quantified. We use three Earth System Models (ESMs) of intermediate complexity to compare effects of selected mechanisms on the interglacial CO2 and δ13CO2 changes, focusing on those with substantial potential impacts: namely carbonate sedimentation in shallow waters, peat growth, and (in the case of the Holocene) human land use. A set of specified carbon cycle forcings could qualitatively explain atmospheric CO2 dynamics from 8 ka BP to the pre-industrial. However, when applied to Eemian boundary conditions from 126 to 115 ka BP, the same set of forcings led to disagreement with the observed direction of CO2 changes after 122 ka BP. This failure to simulate late-Eemian CO2 dynamics could be a result of the imposed forcings such as prescribed CaCO3 accumulation and/or an incorrect response of simulated terrestrial carbon to the surface cooling at the end of the interglacial. These experiments also reveal that key natural processes of interglacial CO2 dynamics - shallow water CaCO3 accumulation, peat and permafrost carbon dynamics - are not well represented in the current ESMs. Global-scale modeling of these long-term carbon cycle components started only in the last decade, and uncertainty in parameterization of these mechanisms is a main limitation in the successful modeling of interglacial CO2 dynamics.

In 2018, central and northern parts of Europe experienced heat and drought conditions over many months from spring to autumn, strongly affecting both natural ecosystems and crops. Besides their impact on nature and society, events like this can be used to study the impact of climate variations on the terrestrial carbon cycle, which is an important determinant of the future climate trajectory. Here, variations in the regional net ecosystem exchange (NEE) of CO2 between terrestrial ecosystems and the atmosphere were quantified from measurements of atmospheric CO2 mole fractions. Over Europe, several observational records have been maintained since at least 1999, giving us the opportunity to assess the 2018 anomaly in the context of at least two decades of variations, including the strong climate anomaly in 2003. In addition to an atmospheric inversion with temporally explicitly estimated anomalies, we use an inversion based on empirical statistical relations between anomalies in the local NEE and anomalies in local climate conditions. For our analysis period 1999–2018, we find that higher-than-usual NEE in hot and dry summers may tend to arise in Central Europe from enhanced ecosystem respiration due to the elevated temperatures, and in Southern Europe from reduced photosynthesis due to the reduced water availability. Despite concerns in the literature, the level of agreement between regression-based NEE anomalies and temporally explicitly estimated anomalies indicates that the atmospheric CO2 measurements from the relatively dense European station network do provide information about the year-to-year variations of Europe's carbon sources and sinks, at least in summer. This article is part of the theme issue 'Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale'. [ABSTRACT FROM AUTHOR]

Interannual variations in the large-scale net ecosystem exchange (NEE) of CO2 between the terrestrial biosphere and the atmosphere were estimated for 1957-2017 from sustained measurements of atmospheric CO2 mixing ratios. As the observations are sparse in the early decades, available records were combined into a 'quasi-homogeneous' dataset based on similarity in their signals, to minimize spurious variations from beginning or ending data records. During El Niño events, CO2 is anomalously released from the tropical band, and a few months later also in the northern extratropical band. This behaviour can approximately be represented by a linear relationship of the NEE anomalies and local air temperature anomalies, with sensitivity coefficients depending on geographical location and season. The apparent climate sensitivity of global total NEE against variations in pan-tropically averaged annual air temperature slowly changed over time during the 1957-2017 period, first increasing (though less strongly than in previous studies) but then decreasing again. However, only part of this change can be attributed to actual changes in local physiological or ecosystem processes, the rest probably arising from shifts in the geographical area of dominating temperature variations. This article is part of a discussion meeting issue 'The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications'. [ABSTRACT FROM AUTHOR]

We modeled the carbon (C) cycle in Mexico with a process-based approach. We used different available products (satellite data, field measurements, models and flux towers) to estimate C stocks and fluxes in the country at three different time frames: present (defined as the period 2000-2005), the past century (1901-2000) and the remainder of this century (2010-2100). Our estimate of the gross primary productivity (GPP) for the country was 2137±1023 TgC yr-1 and a total C stock of 34 506±7483 TgC, with 20 347±4622 TgC in vegetation and 14 159±3861 in the soil. Contrary to other current estimates for recent decades, our results showed that Mexico was a C sink over the period 1990-2009 (+31 TgC yr-1) and that C accumulation over the last century amounted to 1210±1040 TgC. We attributed this sink to the CO2 fertilization effect on GPP, which led to an increase of 3408±1060 TgC, while both climate and land use reduced the country C stocks by -458±1001 and -1740±878 TgC, respectively. Under different future scenarios, the C sink will likely continue over the 21st century, with decreasing C uptake as the climate forcing becomes more extreme. Our work provides valuable insights on relevant driving processes of the C cycle such as the role of drought in drylands (e.g., grasslands and shrublands) and the impact of climate change on the mean residence time of soil C in tropical ecosystems. [ABSTRACT FROM AUTHOR]

We describe the Max Planck Institute Carbon Cycle Data Assimilation System (MPI-CCDAS) built around the tangent-linear version of the land surface scheme of the MPI-Earth System Model v1 (JSBACH). The simulated terrestrial biosphere processes (phenology and carbon balance) were constrained by observations of the fraction of photosynthetically active radiation (TIP-FAPAR product) and by observations of atmospheric CO2 at a global set of monitoring stations for the years 2005-2009. The system successfully, and computationally efficiently, improved average foliar area and northern extra-tropical seasonality of foliar area when constrained by TIP-FAPAR. Global net and gross carbon fluxes were improved when constrained by atmospheric CO2, although the system tended to underestimate tropical productivity. Assimilating both data streams jointly allowed the MPI-CCDAS to match both observations (TIP-FAPAR and atmospheric CO2) equally well as the single data stream assimilation cases, therefore overall increasing the appropriateness of the resultant parameter values and biosphere dynamics. Our study thus highlights the role of the TIP-FAPAR product in stabilising the underdetermined atmospheric inversion problem and demonstrates the value of multiple-data stream assimilation for the simulation of terrestrial biosphere dynamics. The constraint on regional gross and net CO2 flux patterns is limited through the parametrisation of the biosphere model. We expect improvement on that aspect through a refined initialisation strategy and inclusion of further biosphere observations as constraints. [ABSTRACT FROM AUTHOR]

We modelled the carbon (C) cycle in Mexico with a process-based approach. We used different available products (satellite data, field measurements, models and flux towers) to estimate C stocks and fluxes in the country at three different time frames: present (defined as the period 2000-2005), the past century (1901-2000) and the remainder of this century (2010-2100). Our estimate of the gross primary productivity (GPP) for the country was 2137±1023 TgCyr-1 and a total C stock of 34 506±7483 TgC, with 20 347±4622 PgC in vegetation and 14 159±3861 in the soil. Contrary to other current estimates for recent decades, our results showed that Mexico was a C sink over the period 1990-2009 (+31 TgCyr-1) and that C accumulation over the last century amounted to 1210±1040 TgC. We attributed this sink to the CO2 fertilization effect on GPP, which led to an increase of 3408±1060 TgC, while both climate and land use reduced the country C stocks by -458±1001 and -1740±878 TgC, respectively. Under different future scenarios the C sink will likely continue over 21st century, with decreasing C uptake as the climate forcing becomes more extreme. Our work provides valuable insights on relevant driving processes of the C-cycle such as the role of drought in marginal lands (e.g. grasslands and shrublands) and the impact of climate change on the mean residence time of C in tropical ecosystems. [ABSTRACT FROM AUTHOR]

Terrestrial ecosystem models used for Earth system modelling show a significant divergence in future patterns of ecosystem processes, in particular carbon exchanges, despite a seemingly common behaviour for the contemporary period. An in-depth evaluation of these models is hence of high importance to achieve a better understanding of the reasons for this disagreement. Here, we develop an extension for existing benchmarking systems by making use of the complementary information contained in the observational records of atmospheric CO2 and remotely-sensed vegetation activity to provide a firm set of diagnostics of ecosystem responses to climate variability in the last 30 yr at different temporal and spatial scales. The selection of observational characteristics (traits) specifically considers the robustness of information given the uncertainties in both data and evaluation analysis. In addition, we provide a baseline benchmark, a minimum test that the model under consideration has to pass, to provide a more objective, quantitative evaluation framework. The benchmarking strategy can be used for any land surface model, either driven by observed meteorology or coupled to a climate model. We apply this framework to evaluate the offline version of the MPI-Earth system model's land surface scheme JSBACH. We demonstrate that the complementary use of atmospheric CO2 and satellite based vegetation activity data allows to pinpoint specific model failures that would not be possible by the sole use of atmospheric CO2 observations. [ABSTRACT FROM AUTHOR]

Nitrogen (N) availability plays a key role in terrestrial biosphere dynamics. To understand and quantify the role of terrestrial N in the Earth system, we developed an advanced terrestrial biogeochemical model O-CN that mechanistically couples terrestrial energy, water, carbon, and nitrogen fluxes in terrestrial ecosystems. We evaluate this new model against observations from intensive forest monitoring plots at temperate and boreal locations in Europe. O-CN simulates realistic foliage N concentrations and N cycling rates and reproduces observed diurnal and seasonal cycles of C fluxes as well as observed gradients in vegetation productivity with N availability for the forest sites studied. A sensitivity test reveals that these results are reasonably robust against uncertainties in model parameter estimates. Using this model we quantify the likely contribution of anthropogenic N deposition to present ecosystem C sequestration as 36 (range: 2-79) g C g-1 N in agreement with ecosystem manipulation studies. [ABSTRACT FROM AUTHOR]

Global-scale results of the new O-CN terrestrial biosphere model coupling the carbon (C) and nitrogen (N) cycles show that the model produces realistic estimates of present- day C and N stocks and fluxes, despite some regional biases. N availability strongly affects high-latitude foliage area and foliage N, limiting vegetation productivity and present-day high-latitude net C uptake. Anthropogenic N deposition is predicted to have increased net primary productivity due to increases in foliage area and foliage N, contributing 0.2-0.5 Pg C yr-1 to the I 990s global net C uptake. While O-CN's modeled global 1990s terrestrial net C uptake (2.4 Pg C yr-1) is similar to the estimate not accounting for anthropogenic N inputs and N dynamics (2.6 Pg C yr-1), its latitudinal distribution and the sensitivity of the terrestrial C balance to its driving factors are substantially altered by N dynamics, with important implications for future trajectories of the global carbon cycle. [ABSTRACT FROM AUTHOR]

We evaluate three terrestrial biosphere models (LPJ, Orchidee, Biome-BGC) with respect to their capacity to simulate climate related trends in gross primary production (GPP) of forests in Europe. We compare simulated GPP and leaf area index (LAI) with GPP estimates based on flux separated eddy covariance measurements of net ecosystem exchange (NEE) and LAI measurements along a gradient in mean annual temperature from the boreal to the Mediterranean. 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 zone appears to be primarily related to too high simulated LAI - and thus light absorption (APAR) - rather than too high radiation use efficiency. On average, the models compare similarly well to the site GPP data (RMSE of ∼30% or 420 gC/m²/yr) but differences are apparent for different ecosystem types. Given 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, we find the agreement between site data and simulations acceptable, providing confidence in simulations of GPP for European forests. [ABSTRACT FROM AUTHOR]