Search

Model Intercomparison Projects (MIPs) are fundamental to our understanding of how the land surface responds to changes in climate. However, MIPs are challenging to conduct, requiring the organization of multiple, decentralized modeling teams throughout the world running common protocols. We explored centralizing these models on a single supercomputing system. We ran nine offline terrestrial biosphere models through the Terrestrial Biosphere Model Farm: CABLE, CENTURY, HyLand, ISAM, JULES, LPJ-GUESS, ORCHIDEE, SiB-3, and SiB-CASA. All models were wrapped in a software framework driven with common forcing data, spin-up, and run protocols specified by the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) for years 1901-2100. We ran more than a dozen model experiments. We identify three major benefits and three major challenges. The benefits include: (a) processing multiple models through a MIP is relatively straightforward, (b) MIP protocols are run consistently across models, which may reduce some model output variability, and (c) unique multimodel experiments can provide novel output for analysis. The challenges are: (a) technological demand is large, particularly for data and output storage and transfer; (b) model versions lag those from the core model development teams; and (c) there is still a need for intellectual input from the core model development teams for insight into model results. A merger with the open-source, cloud-based Predictive Ecosystem Analyzer (PEcAn) ecoinformatics system may be a path forward to overcoming these challenges.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Terrestrial vegetation removes CO2 from the atmosphere; an important climate regulation service that slows global warming. This 119 Pg C per annum transfer of CO2 into plants-gross primary productivity (GPP)-is the largest land carbon flux globally. While understanding past and anticipated future GPP changes is necessary to support carbon management, the factors driving long-term changes in GPP are largely unknown. Here we show that 1901 to 2010 changes in GPP have been dominated by anthropogenic activity. Our dual constraint attribution approach provides three insights into the spatiotemporal patterns of GPP change. First, anthropogenic controls on GPP change have increased from 57% (1901 decade) to 94% (2001 decade) of the vegetated land surface. Second, CO2 fertilization and nitro  gen deposition are the most important drivers of change, 19.8 and 11.1 Pg C per annum (2001 decade) respectively, especially in the tropics and industrialized areas since the 1970's. Third, changes in climate have functioned as fertilization to enhance GPP (1.4 Pg C per annum in the 2001 decade). These findings suggest that, from a land carbon balance perspective, the Anthropocene began over 100 years ago and that global change drivers have allowed GPP uptake to keep pace with anthropogenic emissions.

Clarifying how increased atmospheric CO2 concentration (eCO2) contributes to accelerated land carbon sequestration remains important since this process is the largest negative feedback in the coupled carbon–climate system. Here, we constrain the sensitivity of the terrestrial carbon sink to eCO2 over the temperate Northern Hemisphere for the past five decades, using 12 terrestrial ecosystem models and data from seven CO2 enrichment experiments. This constraint uses the heuristic finding that the northern temperate carbon sink sensitivity to eCO2 is linearly related to the site-scale sensitivity across the models. The emerging data-constrained eCO2 sensitivity is 0.64 ± 0.28 PgC yr−1 per hundred ppm of eCO2. Extrapolating worldwide, this northern temperate sensitivity projects the global terrestrial carbon sink to increase by 3.5 ± 1.9 PgC yr−1 for an increase in CO2 of 100 ppm. This value suggests that CO2 fertilization alone explains most of the observed increase in global land carbon sink since the 1960s. More CO2 enrichment experiments, particularly in boreal, arctic and tropical ecosystems, are required to explain further the responsible processes.

The ability to accurately predict ecosystem drought response and recovery is necessary to produce reliable forecasts of land carbon uptake and future climate. Using a suite of models from the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), we assessed modeled net primary productivity (NPP) response to, and recovery from, drought events against a benchmark derived from tree ring observations between 1948 and 2008 across forested regions of the US and Europe. We find short lag times (0-6 months) between climate anomalies and modeled NPP response. Although models accurately simulate the direction of drought legacy effects (i.e. NPP decreases), projected effects are approximately four times shorter and four times weaker than observations suggest. This discrepancy between observed and simulated vegetation recovery from drought reveals a potential critical model deficiency. Since productivity is a crucial component of the land carbon balance, models that underestimate drought recovery time could overestimate predictions of future land carbon sink strength and, consequently, underestimate forecasts of atmospheric CO2.

Water use efficiency (WUE), defined as the ratio of gross primary productivity and evapotranspiration at the ecosystem scale, is a critical variable linking the carbon and water cycles. Incorporating a dependency on vapor pressure deficit, apparent underlying WUE (uWUE) provides a better indicator of how terrestrial ecosystems respond to environmental changes than other WUE formulations. Here we used 20th century simulations from four terrestrial biosphere models to develop a novel variance decomposition method. With this method, we attributed variations in apparent uWUE to both the trend and interannual variation of environmental drivers. The secular increase in atmospheric CO2 explained a clear majority of total variation (66 ± 32%: mean ± one standard deviation), followed by positive trends in nitrogen deposition and climate, as well as a negative trend in land use change. In contrast, interannual variation was mostly driven by interannual climate variability. To analyze the mechanism of the CO2 effect, we partitioned the apparent uWUE into the transpiration ratio (transpiration over evapotranspiration) and potential uWUE. The relative increase in potential uWUE parallels that of CO2, but this direct CO2 effect was offset by 20 ± 4% by changes in ecosystem structure, that is, leaf area index for different vegetation types. However, the decrease in transpiration due to stomatal closure with rising CO2 was reduced by 84% by an increase in leaf area index, resulting in small changes in the transpiration ratio. CO2 concentration thus plays a dominant role in driving apparent uWUE variations over time, but its role differs for the two constituent components: potential uWUE and transpiration.

Soil carbon (C) is a critical component of Earth system models (ESMs), and its diverse representations are a major source of the large spread across models in the terrestrial C sink from the third to fifth assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Improving soil C projections is of a high priority for Earth system modeling in the future IPCC and other assessments. To achieve this goal, we suggest that (1) model structures should reflect real-world processes, (2) parameters should be calibrated to match model outputs with observations, and (3) external forcing variables should accurately prescribe the environmental conditions that soils experience. First, most soil C cycle models simulate C input from litter production and C release through decomposition. The latter process has traditionally been represented by first-order decay functions, regulated primarily by temperature, moisture, litter quality, and soil texture. While this formulation well captures macroscopic soil organic C (SOC) dynamics, better understanding is needed of their underlying mechanisms as related to microbial processes, depth-dependent environmental controls, and other processes that strongly affect soil C dynamics. Second, incomplete use of observations in model parameterization is a major cause of bias in soil C projections from ESMs. Optimal parameter calibration with both pool- and flux-based data sets through data assimilation is among the highest priorities for near-term research to reduce biases among ESMs. Third, external variables are represented inconsistently among ESMs, leading to differences in modeled soil C dynamics. We recommend the implementation of traceability analyses to identify how external variables and model parameterizations influence SOC dynamics in different ESMs. Overall, projections of the terrestrial C sink can be substantially improved when reliable data sets are available to select the most representative model structure, constrain parameters, and prescribe forcing fields.

Author(s): Christopher R. Schwalm (corresponding author) [1, 2]; William R. L. Anderegg [3]; Anna M. Michalak [4]; Joshua B. Fisher [5]; Franco Biondi [6]; George Koch [2]; Marcy Litvak [7]; [...]