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.

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.

Terrestrial carbon cycle models have incorporated increasingly more processes as a means to achieve more-realistic representations of ecosystem carbon cycling. Despite this, there are large across-model variations in the simulation and projection of carbon cycling. Several model intercomparison projects (MIPs), for example, the fifth phase of the Coupled Model Intercomparison Project (CMIP5) (historical simulations), Trends in Net Land-Atmosphere Carbon Exchange (TRENDY), and Multiscale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), have sought to understand intermodel differences. In this study, the authors developed a suite of new techniques to conduct post-MIP analysis to gain insights into uncertainty sources across 25 models in the three MIPs. First, terrestrial carbon storage dynamics were characterized by a three-dimensional (3D) model output space with coordinates of carbon residence time, net primary productivity (NPP), and carbon storage potential. The latter represents the potential of an ecosystem to lose or gain carbon. This space can be used to measure how and why model output differs. Models with a nitrogen cycle generally exhibit lower annual NPP in comparison with other models, and mostly negative carbon storage potential. Second, a transient traceability framework was used to decompose any given carbon cycle model into traceable components and identify the sources of model differences. The carbon residence time (or NPP) was traced to baseline carbon residence time (or baseline NPP related to the maximum carbon input), environmental scalars, and climate forcing. Third, by applying a variance decomposition method, the authors show that the intermodel differences in carbon storage can be mainly attributed to the baseline carbon residence time and baseline NPP (>90% in the three MIPs). The three techniques developed in this study offer a novel approach to gain more insight from existing MIPs and can point out directions for future MIPs. Since this study is conducted at the global scale for an overview on intermodel differences, future studies should focus more on regional analysis to identify the sources of uncertainties and improve models at the specified mechanism level.

NASA has launched the decade-long Arctic-Boreal Vulnerability Experiment (ABoVE). While the initial phases focus on field and airborne data collection, early integration with modeling activities is important to benefit future modeling syntheses. We compiled feedback from ecosystem modeling teams on key data needs, which encompass carbon biogeochemistry, vegetation, permafrost, hydrology, and disturbance dynamics. A suite of variables was identified as part of this activity with a critical requirement that they are collected concurrently and representatively over space and time. Individual projects in ABoVE may not capture all these needs, and thus there is both demand and opportunity for the augmentation of field observations, and synthesis of the observations that are collected, to ensure that science questions and integrated modeling activities are successfully implemented.

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.

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]; [...]

The terrestrial biosphere can release or absorb the greenhouse gases, carbon dioxide (C[O.sub.2]), methane (C[H.sub.4]) and nitrous oxide ([N.sub.2]O), and therefore has an important role in regulating atmospheric composition and [...]

The Arctic-boreal zone (ABZ) is experiencing amplified warming, actively changing biogeochemical cycling of vegetation and soils. The land-to-atmosphere fluxes of CO2 in the ABZ have the potential to increase in magnitude and feedback to the climate causing additional large scale warming. The ability to model and predict this vulnerability is critical to preparation for a warming world, but Earth system models have biases that may hinder understanding the rapidly changing ABZ carbon fluxes. Here we investigate circumpolar carbon cycling represented by the Community Land Model 5 (CLM5.0) with a focus on seasonal gross primary productivity (GPP) in plant functional types (PFTs). We benchmark model results using data from satellite remote sensing products and eddy covariance towers. We find consistent biases in CLM5.0 relative to observational constraints: (1) the onset of deciduous plant productivity to be late, (2) the offset of productivity to lag and remain abnormally high for all PFTs in fall, (3) a high bias of grass, shrub, and needleleaf evergreen tree productivity, and (4) an underestimation of productivity of deciduous trees. Based on these biases, we focus model development of alternate phenology, photosynthesis schemes, and carbon allocation parameters at eddy covariance tower sites. Although our improvements are focused on productivity, our final Model Recommendation results in other component CO2 fluxes, e.g. Net Ecosystem Exchange (NEE) and Terrestrial Ecosystem Respiration (TER), that are more consistent with observations. Results suggest that algorithms developed for lower latitudes and more temperate environments can be inaccurate when extrapolated to the ABZ, and that many land surface models may not accurately represent carbon cycling and its recent rapid changes in high latitude ecosystems, especially when analyzed by individual PFTs. [ABSTRACT FROM AUTHOR]