Deborah N. Huntzinger

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

Abstract 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.

Disentangling climatic and anthropogenic controls on global terrestrial evapotranspiration trends

We examined natural and anthropogenic controls on terrestrial evapotranspiration (ET) changes from 1982 to 2010 using multiple estimates from remote sensing-based datasets and process-oriented land surface models. A significant increasing trend of ET in each hemisphere was consistently revealed by observationally-constrained data and multi-model ensembles that considered historic natural and anthropogenic drivers. The climate impacts were simulated to determine the spatiotemporal variations in ET. Globally, rising CO2 ranked second in these models after the predominant climatic influences, and yielded decreasing trends in canopy transpiration and ET, especially for tropical forests and high-latitude shrub land. Increasing nitrogen deposition slightly amplified global ET via enhanced plant growth. Land-use-induced ET responses, albeit with substantial uncertainties across the factorial analysis, were minor globally, but pronounced locally, particularly over regions with intensive land-cover changes. Our study highlights the importance of employing multi-stream ET and ET-component estimates to quantify the strengthening anthropogenic fingerprint in the global hydrologic cycle.

Impact of large‐scale climate extremes on biospheric carbon fluxes: An intercomparison based on MsTMIP data

Understanding the role of climate extremes and their impact on the carbon (C) cycle is increasingly a focus of Earth system science. Climate extremes such as droughts, heat waves, or heavy precipitation events can cause substantial changes in terrestrial C fluxes. On the other hand, extreme changes in C fluxes are often, but not always, driven by extreme climate conditions. Here we present an analysis of how extremes in temperature and precipitation, and extreme changes in terrestrial C fluxes are related to each other in 10 state-of-the-art terrestrial carbon models, all driven by the same climate forcing. We use model outputs from the North American Carbon Program Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP). A global-scale analysis shows that both droughts and heat waves translate into anomalous net releases of CO2 from the land surface via different mechanisms: Droughts largely decrease gross primary production (GPP) and to a lower extent total respiration (TR), while heat waves slightly decrease GPP but increase TR. Cold and wet periods have a smaller opposite effect. Analyzing extremes in C fluxes reveals that extreme changes in GPP and TR are often caused by strong shifts in water availability, but for extremes in TR shifts in temperature are also important. Extremes in net CO2 exchange are equally strongly driven by deviations in temperature and precipitation. Models mostly agree on the sign of the C flux response to climate extremes, but model spread is large. In tropical forests, C cycle extremes are driven by water availability, whereas in boreal forests temperature plays a more important role. Models are particularly uncertain about the C flux response to extreme heat in boreal forests.

Evaluation of continental carbon cycle simulations with North American flux tower observations

Terrestrial biosphere models can help identify physical processes that control carbon dynamics, including land-atmosphere CO2 fluxes, and have great potential to predict the terrestrial ecosystem response to changing climate. The skill of models that provide continental-scale carbon flux estimates, however, remains largely untested. This paper evaluates the performance of continental-scale flux estimates from 17 models against observations from 36 North American flux towers. Fluxes extracted from regional model simulations were compared with co-located flux tower observations at monthly and annual time increments. Site-level model simulations were used to help interpret sources of the mismatch between the regional simulations and site-based observations. On average, the regional model runs overestimated the annual gross primary productivity (5%) and total respiration (15%), and they significantly underestimated the annual net carbon uptake (64%) during the time period 2000- 2005. Comparison with site-level simulations implicated choices specific to regional model simulations as contributors to the gross flux biases, but not the net carbon uptake bias. The models performed the best at simulating carbon exchange at deciduous broadleaf sites, likely because a number of models used prescribed phenology to simulate seasonal fluxes. The models did not perform as well for crop, grass, and evergreen sites. The regional models matched the observations most closely in terms of seasonal correlation and seasonal magnitude of variation, but they have very little skill at interannual correlation and minimal skill at interannual magnitude of variability. The comparison of site vs. regional-level model runs demonstrated that (1) the interannual correlation is higher for site-level model runs, but the skill remains low; and (2) the underestimation of year-to-year variability for all fluxes is an inherent weakness of the models. The best-performing regional models that did not use flux tower calibration were CLM-CN, CASA-GFEDv2, and SIB3.1. Two flux tower calibrated, empirical models, EC-MOD and MOD17þ, performed as well as the best process-based models. This suggests that (1) empirical, calibrated models can perform as well as complex, process-based models and (2) combining process- based model structure with relevant constraining data could significantly improve model performance.

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.

North American carbon dioxide sources and sinks: magnitude, attribution, and uncertainty

North America is both a source and sink of atmospheric carbon dioxide (CO2). Continental sources – such as fossil-fuel combustion in the US and deforestation in Mexico – and sinks – including most ecosystems, and particularly secondary forests – add and remove CO2 from the atmosphere, respectively. Photosynthesis converts CO2 into carbon as biomass, which is stored in vegetation, soils, and wood products. However, ecosystem sinks compensate for only ~35% of the continents fossil-fuel-based CO2 emissions; North America therefore represents a net CO2 source. Estimating the magnitude of ecosystem sinks, even though the calculation is confounded by uncertainty as a result of individual inventory- and model-based alternatives, has improved through the use of a combined approach.

Toward “optimal” integration of terrestrial biosphere models

Multimodel ensembles (MME) are commonplace in Earth system modeling. Here we perform MME integration using a 10-member ensemble of terrestrial biosphere models (TBMs) from the Multiscale synthesis and Terrestrial Model Intercomparison Project (MsTMIP). We contrast optimal (skill based for present-day carbon cycling) versus naive (“one model-one vote”) integration. MsTMIP optimal and naive mean land sink strength estimates (−1.16 versus −1.15 Pg C per annum respectively) are statistically indistinguishable. This holds also for grid cell values and extends to gross uptake, biomass, and net ecosystem productivity. TBM skill is similarly indistinguishable. The added complexity of skill-based integration does not materially change MME values. This suggests that carbon metabolism has predictability limits and/or that all models and references are misspecified. Resolving this issue requires addressing specific uncertainty types (initial conditions, structure, and references) and a change in model development paradigms currently dominant in the TBM community.

Uncertainty in the response of terrestrial carbon sink to environmental drivers undermines carbon-climate feedback predictions

Terrestrial ecosystems play a vital role in regulating the accumulation of carbon (C) in the atmosphere. Understanding the factors controlling land C uptake is critical for reducing uncertainties in projections of future climate. The relative importance of changing climate, rising atmospheric CO2, and other factors, however, remains unclear despite decades of research. Here, we use an ensemble of land models to show that models disagree on the primary driver of cumulative C uptake for 85% of vegetated land area. Disagreement is largest in model sensitivity to rising atmospheric CO2 which shows almost twice the variability in cumulative land uptake since 1901 (1 s.d. of 212.8 PgC vs. 138.5 PgC, respectively). We find that variability in CO2 and temperature sensitivity is attributable, in part, to their compensatory effects on C uptake, whereby comparable estimates of C uptake can arise by invoking different sensitivities to key environmental conditions. Conversely, divergent estimates of C uptake can occur despite being based on the same environmental sensitivities. Together, these findings imply an important limitation to the predictability of C cycling and climate under unprecedented environmental conditions. We suggest that the carbon modeling community prioritize a probabilistic multi-model approach to generate more robust C cycle projections.

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.

Response of Water Use Efficiency to Global Environmental Change Based on Output From Terrestrial Biosphere Models

Author(s): Zhou, S; Yu, B; Schwalm, CR; Ciais, P; Zhang, Y; Fisher, JB; Michalak, AM; Wang, W; Poulter, B; Huntzinger, DN; Niu, S; Mao, J; Jain, A; Ricciuto, DM; Shi, X; Ito, A; Wei, Y; Huang, Y; Wang, G | Abstract: ©2017. American Geophysical Union. All Rights Reserved. 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.