The Terrestrial Biosphere Model Farm.

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.

Author Correction: Modeling suggests fossil fuel emissions have been driving increased land carbon uptake since the turn of the 20th Century.

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

Modeling suggests fossil fuel emissions have been driving increased land carbon uptake since the turn of the 20th Century.

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.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Carbon and Water Use Efficiencies: A Comparative Analysis of Ten Terrestrial Ecosystem Models under Changing Climate.

Terrestrial ecosystems carbon and water cycles are tightly coupled through photosynthesis and evapotranspiration processes. The ratios of carbon stored to carbon uptake and water loss to carbon gain are key ecophysiological indicators essential to assess the magnitude and response of the terrestrial plant to the changing climate. Here, we use estimates from 10 terrestrial ecosystem models to quantify the impacts of climate, atmospheric CO2 concentration, and nitrogen (N) deposition on water use efficiency (WUE), and carbon use efficiency (CUE). We find that across models, WUE increases over the 20th Century particularly due to CO2 fertilization and N deposition and compares favorably to experimental studies. Also, the results show a decrease in WUE with climate for the last 3 decades, in contrasts with up-scaled flux observations that demonstrate a constant WUE. Modeled WUE responds minimally to climate with modeled CUE exhibiting no clear trend across space and time. The divergence between simulated and observationally-constrained WUE and CUE is driven by modeled NPP and autotrophic respiration, nitrogen cycle, carbon allocation, and soil moisture dynamics in current ecosystem models. We suggest that carbon-modeling community needs to reexamine stomatal conductance schemes and the soil-vegetation interactions for more robust modeling of carbon and water cycles.

Land carbon models underestimate the severity and duration of drought's impact on plant productivity.

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.

Accelerating rates of Arctic carbon cycling revealed by long-term atmospheric CO2 measurements.

The contemporary Arctic carbon balance is uncertain, and the potential for a permafrost carbon feedback of anywhere from 50 to 200 petagrams of carbon (Schuur et al., 2015) compromises accurate 21st-century global climate system projections. The 42-year record of atmospheric CO2 measurements at Barrow, Alaska (71.29 N, 156.79 W), reveals significant trends in regional land-surface CO2 anomalies (ΔCO2), indicating long-term changes in seasonal carbon uptake and respiration. Using a carbon balance model constrained by ΔCO2, we find a 13.4% decrease in mean carbon residence time (50% confidence range = 9.2 to 17.6%) in North Slope tundra ecosystems during the past four decades, suggesting a transition toward a boreal carbon cycling regime. Temperature dependencies of respiration and carbon uptake suggest that increases in cold season Arctic labile carbon release will likely continue to exceed increases in net growing season carbon uptake under continued warming trends.

Global patterns of drought recovery.

Drought, a recurring phenomenon with major impacts on both human and natural systems, is the most widespread climatic extreme that negatively affects the land carbon sink. Although twentieth-century trends in drought regimes are ambiguous, across many regions more frequent and severe droughts are expected in the twenty-first century. Recovery time-how long an ecosystem requires to revert to its pre-drought functional state-is a critical metric of drought impact. Yet the factors influencing drought recovery and its spatiotemporal patterns at the global scale are largely unknown. Here we analyse three independent datasets of gross primary productivity and show that, across diverse ecosystems, drought recovery times are strongly associated with climate and carbon cycle dynamics, with biodiversity and CO2 fertilization as secondary factors. Our analysis also provides two key insights into the spatiotemporal patterns of drought recovery time: first, that recovery is longest in the tropics and high northern latitudes (both vulnerable areas of Earth's climate system) and second, that drought impacts (assessed using the area of ecosystems actively recovering and time to recovery) have increased over the twentieth century. If droughts become more frequent, as expected, the time between droughts may become shorter than drought recovery time, leading to permanently damaged ecosystems and widespread degradation of the land carbon sink.

The terrestrial biosphere as a net source of greenhouse gases to the atmosphere.

The terrestrial biosphere can release or absorb the greenhouse gases, carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O), and therefore has an important role in regulating atmospheric composition and climate. Anthropogenic activities such as land-use change, agriculture and waste management have altered terrestrial biogenic greenhouse gas fluxes, and the resulting increases in methane and nitrous oxide emissions in particular can contribute to climate change. The terrestrial biogenic fluxes of individual greenhouse gases have been studied extensively, but the net biogenic greenhouse gas balance resulting from anthropogenic activities and its effect on the climate system remains uncertain. Here we use bottom-up (inventory, statistical extrapolation of local flux measurements, and process-based modelling) and top-down (atmospheric inversions) approaches to quantify the global net biogenic greenhouse gas balance between 1981 and 2010 resulting from anthropogenic activities and its effect on the climate system. We find that the cumulative warming capacity of concurrent biogenic methane and nitrous oxide emissions is a factor of about two larger than the cooling effect resulting from the global land carbon dioxide uptake from 2001 to 2010. This results in a net positive cumulative impact of the three greenhouse gases on the planetary energy budget, with a best estimate (in petagrams of CO2 equivalent per year) of 3.9 ± 3.8 (top down) and 5.4 ± 4.8 (bottom up) based on the GWP100 metric (global warming potential on a 100-year time horizon). Our findings suggest that a reduction in agricultural methane and nitrous oxide emissions, particularly in Southern Asia, may help mitigate climate change.

North American terrestrial CO2 uptake largely offset by CH4 and N2O emissions: toward a full accounting of the greenhouse gas budget.

The terrestrial ecosystems of North America have been identified as a sink of atmospheric CO2 though there is no consensus on the magnitude. However, the emissions of non-CO2 greenhouse gases (CH4 and N2O) may offset or even overturn the climate cooling effect induced by the CO2 sink. Using a coupled biogeochemical model, in this study, we have estimated the combined global warming potentials (GWP) of CO2, CH4 and N2O fluxes in North American terrestrial ecosystems and quantified the relative contributions of environmental factors to the GWP changes during 1979-2010. The uncertainty range for contemporary global warming potential has been quantified by synthesizing the existing estimates from inventory, forward modeling, and inverse modeling approaches. Our "best estimate" of net GWP for CO2, CH4 and N2O fluxes was -0.50 ± 0.27 Pg CO2 eq/year (1 Pg = 1015 g) in North American terrestrial ecosystems during 2001-2010. The emissions of CH4 and N2O from terrestrial ecosystems had offset about two thirds (73 %±14 %) of the land CO2 sink in the North American continent, showing large differences across the three countries, with offset ratios of 57 % ± 8 % in US, 83 % ± 17 % in Canada and 329 % ± 119 % in Mexico. Climate change and elevated tropospheric ozone concentration have contributed the most to GWP increase, while elevated atmospheric CO2 concentration have contributed the most to GWP reduction. Extreme drought events over certain periods could result in a positive GWP. By integrating the existing estimates, we have found a wide range of uncertainty for the combined GWP. From both climate change science and policy perspectives, it is necessary to integrate ground and satellite observations with models for a more accurate accounting of these three greenhouse gases in North America.

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions.

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land-atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century-long (1901-2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project and two observation-based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr-1 with a median value of 51 Pg C yr-1 during 2001-2010. The largest uncertainty in SOC stocks exists in the 40-65°N latitude whereas the largest cross-model divergence in Rh are in the tropics. The modeled SOC change during 1901-2010 ranges from -70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble-estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks-even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Visual Reconciliation of Alternative Similarity Spaces in Climate Modeling.

Visual data analysis often requires grouping of data objects based on their similarity. In many application domains researchers use algorithms and techniques like clustering and multidimensional scaling to extract groupings from data. While extracting these groups using a single similarity criteria is relatively straightforward, comparing alternative criteria poses additional challenges. In this paper we define visual reconciliation as the problem of reconciling multiple alternative similarity spaces through visualization and interaction. We derive this problem from our work on model comparison in climate science where climate modelers are faced with the challenge of making sense of alternative ways to describe their models: one through the output they generate, another through the large set of properties that describe them. Ideally, they want to understand whether groups of models with similar spatio-temporal behaviors share similar sets of criteria or, conversely, whether similar criteria lead to similar behaviors. We propose a visual analytics solution based on linked views, that addresses this problem by allowing the user to dynamically create, modify and observe the interaction among groupings, thereby making the potential explanations apparent. We present case studies that demonstrate the usefulness of our technique in the area of climate science.

New Conceptual Model for Soil Treatment Units: Formation of Multiple Hydraulic Zones during Unsaturated Wastewater Infiltration.

Onsite wastewater treatment systems are commonly used in the United States to reclaim domestic wastewater. A distinct biomat forms at the infiltrative surface, causing resistance to flow and decreasing soil moisture below the biomat. To simulate these conditions, previous modeling studies have used a two-layer approach: a thin biomat layer (1-5 cm thick) and the native soil layer below the biomat. However, the effect of wastewater application extends below the biomat layer. We used numerical modeling supported by experimental data to justify a new conceptual model that includes an intermediate zone (IZ) below the biomat. The conceptual model was set up using Hydrus 2D and calibrated against soil moisture and water flux measurements. The estimated hydraulic conductivity value for the IZ was between biomat and the native soil. The IZ has important implications for wastewater treatment. When the IZ was not considered, a loading rate of 5 cm d resulted in an 8.5-cm ponding. With the IZ, the same loading rate resulted in a 9.5-cm ponding. Without the IZ, up to 3.1 cm d of wastewater could be applied without ponding; with the IZ, only up to 2.8 cm d could be applied without ponding. The IZ also plays a significant role in soil moisture distribution. Without the IZ, near-saturation conditions were observed only within the biomat, whereas near-saturation conditions extended below the biomat with the IZ. Accurate prediction of ponding is important to prevent surfacing of wastewater. The degree of water and air saturation influences pollutant treatment efficiency through residence time, volatility, and biochemical reactions.

Carbon dioxide sequestration in cement kiln dust through mineral carbonation.

Carbon sequestration through the formation of carbonates is a potential means to reduce CO2 emissions. Alkaline industrial solid wastes typically have high mass fractions of reactive oxides that may not require preprocessing, making them an attractive source material for mineral carbonation The degree of mineral carbonation achievable in cement kiln dust (CKD) underambienttemperatures and pressures was examined through a series of batch and column experiments. The overall extent and potential mechanisms and rate behavior of the carbonation process were assessed through a complementary set of analytical and empirical methods, including mass change, thermal analysis, and X-ray diffraction. The carbonation reactions were carried out primarily through the reaction of CO2 with Ca(OH)2, and CaCO3 was observed as the predominant carbonation product. A sequestration extent of over 60% was observed within 8 h of reaction without any modifications to the waste. Sequestration appears to follow unreacted core model theory where reaction kinetics are controlled by a first-order rate constant at early times; however, as carbonation progresses, the kinetics of the reaction are attenuated by the extent of the reaction due to diffusion control, with the extent of conversion never reaching completion.

Mineral carbonation for carbon sequestration in cement kiln dust from waste piles.

Alkaline earth metals, such as calcium and magnesium oxides, readily react with carbon dioxide (CO(2)) to produce stable carbonate minerals. Carbon sequestration through the formation of carbonate minerals is a potential means to reduce CO(2) emissions. Calcium-rich, industrial solid wastes and residues provide a potential source of highly reactive oxides, without the need for pre-processing. This paper presents the first study examining the feasibility of carbon sequestration in cement kiln dust (CKD), a byproduct generated during the manufacturing of cement. A series of column experiments were conducted on segments of intact core taken from landfilled CKD. Based on stoichiometry and measured consumption of CO(2) during the experiments, degrees of carbonation greater than 70% of the material's potential theoretical extent were achieved under ambient temperature and pressure conditions. The overall extent of carbonation/sequestration was greater in columns with lower water contents. The major sequestration product appears to be calcite; however, more detailed material characterization is need on pre- and post-carbonated samples to better elucidate carbonation pathways and products.