The mechanisms governing tree drought mortality and recovery remain a subject of inquiry and active debate given their role in the terrestrial carbon cycle and their concomitant impact on climate ...change. Counter‐intuitively, many trees do not die during the drought itself. Indeed, observations globally have documented that trees often grow for several years after drought before mortality. A combination of meta‐analysis and tree physiological models demonstrate that optimal carbon allocation after drought explains observed patterns of delayed tree mortality and provides a predictive recovery framework. Specifically, post‐drought, trees attempt to repair water transport tissue and achieve positive carbon balance through regrowing drought‐damaged xylem. Furthermore, the number of years of xylem regrowth required to recover function increases with tree size, explaining why drought mortality increases with size. These results indicate that tree resilience to drought‐kill may increase in the future, provided that CO2 fertilisation facilitates more rapid xylem regrowth.
The impacts of climate extremes on terrestrial ecosystems are poorly understood but important for predicting carbon cycle feedbacks to climate change. Coupled climate–carbon cycle models typically ...assume that vegetation recovery from extreme drought is immediate and complete, which conflicts with the understanding of basic plant physiology. We examined the recovery of stem growth in trees after severe drought at 1338 forest sites across the globe, comprising 49,339 site-years, and compared the results with simulated recovery in climate-vegetation models. We found pervasive and substantial "legacy effects" of reduced growth and incomplete recovery for 1 to 4 years after severe drought. Legacy effects were most prevalent in dry ecosystems, among Pinaceae, and among species with low hydraulic safety margins. In contrast, limited or no legacy effects after drought were simulated by current climate-vegetation models. Our results highlight hysteresis in ecosystem-level carbon cycling and delayed recovery from climate extremes.
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 CO
, 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 CO
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 CO
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.
Terrestrial biosphere models (TBMs) have become an integral tool for extrapolating local observations and understanding of land-atmosphere carbon exchange to larger regions. The North American Carbon ...Program (NACP) Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) is a formal model intercomparison and evaluation effort focused on improving the diagnosis and attribution of carbon exchange at regional and global scales. MsTMIP builds upon current and past synthesis activities, and has a unique framework designed to isolate, interpret, and inform understanding of how model structural differences impact estimates of carbon uptake and release. Here we provide an overview of the MsTMIP effort and describe how the MsTMIP experimental design enables the assessment and quantification of TBM structural uncertainty. Model structure refers to the types of processes considered (e.g., nutrient cycling, disturbance, lateral transport of carbon), and how these processes are represented (e.g., photosynthetic formulation, temperature sensitivity, respiration) in the models. By prescribing a common experimental protocol with standard spin-up procedures and driver data sets, we isolate any biases and variability in TBM estimates of regional and global carbon budgets resulting from differences in the models themselves (i.e., model structure) and model-specific parameter values. An initial intercomparison of model structural differences is represented using hierarchical cluster diagrams (a.k.a. dendrograms), which highlight similarities and differences in how models account for carbon cycle, vegetation, energy, and nitrogen cycle dynamics. We show that, despite the standardized protocol used to derive initial conditions, models show a high degree of variation for GPP, total living biomass, and total soil carbon, underscoring the influence of differences in model structure and parameterization on model estimates.
Given the magnitude of soil carbon stocks in northern ecosystems, and the vulnerability of these stocks to climate warming, land surface models must accurately represent soil carbon dynamics in these ...regions. We evaluate soil carbon stocks and turnover rates, and the relationship between soil carbon loss with soil temperature and moisture, from an ensemble of eleven global land surface models. We focus on the region of NASA's Arctic-Boreal vulnerability experiment (ABoVE) in North America to inform data collection and model development efforts. Models exhibit an order of magnitude difference in estimates of current total soil carbon stocks, generally under- or overestimating the size of current soil carbon stocks by greater than 50 PgC. We find that a model's soil carbon stock at steady-state in 1901 is the prime driver of its soil carbon stock a hundred years later-overwhelming the effect of environmental forcing factors like climate. The greatest divergence between modeled and observed soil carbon stocks is in regions dominated by peat and permafrost soils, suggesting that models are failing to capture the frozen soil carbon dynamics of permafrost regions. Using a set of functional benchmarks to test the simulated relationship of soil respiration to both soil temperature and moisture, we find that although models capture the observed shape of the soil moisture response of respiration, almost half of the models examined show temperature sensitivities, or Q10 values, that are half of observed. Significantly, models that perform better against observational constraints of respiration or carbon stock size do not necessarily perform well in terms of their functional response to key climatic factors like changing temperature. This suggests that models may be arriving at the right result, but for the wrong reason. The results of this work can help to bridge the gap between data and models by both pointing to the need to constrain initial carbon pool sizes, as well as highlighting the importance of incorporating functional benchmarks into ongoing, mechanistic modeling activities such as those included in ABoVE.
Inter‐comparison and similarity analysis to gauge consensus among multiple simulation models is a critical visualization problem for understanding climate change patterns. Climate models, ...specifically, Terrestrial Biosphere Models (TBM) represent time and space variable ecosystem processes, like, simulations of photosynthesis and respiration, using algorithms and driving variables such as climate and land use. While it is widely accepted that interactive visualization can enable scientists to better explore model similarity from different perspectives and different granularity of space and time, currently there is a lack of such visualization tools.
In this paper we present three main contributions. First, we propose a domain characterization for the TBM community by systematically defining the domain‐specific intents for analyzing model similarity and characterizing the different facets of the data. Second, we define a classification scheme for combining visualization tasks and multiple facets of climate model data in one integrated framework, which can be leveraged for translating the tasks into the visualization design. Finally, we present SimilarityExplorer, an exploratory visualization tool that facilitates similarity comparison tasks across both space and time through a set of coordinated multiple views. We present two case studies from three climate scientists, who used our tool for a month for gaining scientific insights into model similarity. Their experience and results validate the effectiveness of our tool.
•We compare 14 runoff products with USGS runoff (depleted and corrected).•TBMs “best” match USGS runoff, coupled reanalysis are “worst”.•Land-based reanalysis and water budget estimates have ...above-average skill.•Runoff, despite the range in TBM skill, is better predicted than carbon cycling.•Resolving runoff mismatch causes requires a constrained ensemble at fine(r) scales.
Significant changes in the water cycle are expected under current global environmental change. Robust assessment of present-day water cycle dynamics at continental to global scales is confounded by shortcomings in the observed record. Modeled assessments also yield conflicting results which are linked to differences in model structure and simulation protocol. Here we compare simulated gridded (1° spatial resolution) runoff from six terrestrial biosphere models (TBMs), seven reanalysis products, and one gridded surface station product in the contiguous United States (CONUS) from 2001 to 2005. We evaluate the consistency of these 14 estimates with stream gauge data, both as depleted flow and corrected for net withdrawals (2005 only), at the CONUS and water resource region scale, as well as examining similarity across TBMs and reanalysis products at the grid cell scale. Mean runoff across all simulated products and regions varies widely (range: 71 to 356mmyr−1) relative to observed continental-scale runoff (209 or 280mmyr−1 when corrected for net withdrawals). Across all 14 products 8 exhibit Nash–Sutcliffe efficiency values in excess of 0.8 and three are within 10% of the observed value. Region-level mismatch exhibits a weak pattern of overestimation in western and underestimation in eastern regions—although two products are systematically biased across all regions—and largely scales with water use. Although gridded composite TBM and reanalysis runoff show some regional similarities, individual product values are highly variable. At the coarse scales used here we find that progress in better constraining simulated runoff requires standardized forcing data and the explicit incorporation of human effects (e.g., water withdrawals by source, fire, and land use change).
Earth system processes exhibit complex patterns across time, as do the models that seek to replicate these processes. Model output may or may not be significantly related to observations at different ...times and on different frequencies. Conventional model diagnostics provide an aggregate view of model-data agreement, but usually do not identify the time and frequency patterns of model-data disagreement, leaving unclear the steps required to improve model response to environmental drivers that vary on characteristic frequencies. Wavelet coherence can quantify the times and timescales at which two time series, for example time series of models and measurements, are significantly different. We applied wavelet coherence to interpret the predictions of 20 ecosystem models from the North American Carbon Program (NACP) Site-Level Interim Synthesis when confronted with eddy-covariance-measured net ecosystem exchange (NEE) from 10 ecosystems with multiple years of available data. Models were grouped into classes with similar approaches for incorporating phenology, the calculation of NEE, the inclusion of foliar nitrogen (N), and the use of model-data fusion. Models with prescribed, rather than prognostic, phenology often fit NEE observations better on annual to interannual timescales in grassland, wetland and agricultural ecosystems. Models that calculated NEE as net primary productivity (NPP) minus heterotrophic respiration (HR) rather than gross ecosystem productivity (GPP) minus ecosystem respiration (ER) fit better on annual timescales in grassland and wetland ecosystems, but models that calculated NEE as GPP minus ER were superior on monthly to seasonal timescales in two coniferous forests. Models that incorporated foliar nitrogen (N) data were successful at capturing NEE variability on interannual (multiple year) timescales at Howland Forest, Maine. The model that employed a model-data fusion approach often, but not always, resulted in improved fit to data, suggesting that improving model parameterization is important but not the only step for improving model performance. Combined with previous findings, our results suggest that the mechanisms driving daily and annual NEE variability tend to be correctly simulated, but the magnitude of these fluxes is often erroneous, suggesting that model parameterization must be improved. Few NACP models correctly predicted fluxes on seasonal and interannual timescales where spectral energy in NEE observations tends to be low, but where phenological events, multi-year oscillations in climatological drivers, and ecosystem succession are known to be important for determining ecosystem function. Mechanistic improvements to models must be made to replicate observed NEE variability on seasonal and interannual timescales.
Gross primary production (GPP) is often expressed as the product of absorbed photosynthetically active radiation and the efficiency (ε) with which a plant community uses absorbed radiation in biomass ...production. Light-use efficiency is affected by environmental stresses, and varies diurnally and seasonally. Uncertainty about ε can be a serious limitation when modeling GPP. An important determinant of ε is the amount and type of solar radiation incident on a canopy, because an abundance of light can trigger a photo-protective reaction, diminishing GPP. The radiation regime in a forest canopy is determined by the predominant sky conditions and by mutual shading of tree crowns. Shading effects, producing shifts in the amount of incident direct and diffuse solar radiation, have been largely ignored, however, because they depend on forest structure and are difficult to measure. We describe a new approach for estimating changes in mutual canopy shading throughout the day and year based on a canopy structure model derived from light detection and ranging (LiDAR). Proportions of canopy shading were then combined with eddy covariance data to assess the explanatory power for variance in ε by regression tree analysis over half-hourly, daily and weekly time scales. The approach explained between 75 and 97% of variance in ε, representing an increase of between 5 and 16% compared with models driven solely by meteorological variables.
The El Niño Southern Oscillation is the dominant year-to-year mode of global climate variability. El Niño effects on terrestrial carbon cycling are mediated by associated climate anomalies, primarily ...drought, influencing fire emissions and biotic net ecosystem exchange (NEE). Here we evaluate whether El Niño produces a consistent response from the global carbon cycle. We apply a novel bottom-up approach to estimating global NEE anomalies based on FLUXNET data using land cover maps and weather reanalysis. We analyze 13 years (1997-2009) of globally gridded observational NEE anomalies derived from eddy covariance flux data, remotely-sensed fire emissions at the monthly time step, and NEE estimated from an atmospheric transport inversion. We evaluate the overall consistency of biospheric response to El Niño and, more generally, the link between global CO2 flux anomalies and El Niño-induced drought. Our findings, which are robust relative to uncertainty in both methods and time-lags in response, indicate that each event has a different spatial signature with only limited spatial coherence in Amazônia, Australia and southern Africa. For most regions, the sign of response changed across El Niño events. Biotic NEE anomalies, across 5 El Niño events, ranged from -1.34 to +0.98 Pg C yr-1 , whereas fire emissions anomalies were generally smaller in magnitude (ranging from -0.49 to +0.53 Pg C yr-1 ). Overall drought does not appear to impose consistent terrestrial CO2 flux anomalies during El Niños, finding large variation in globally integrated responses from -1.15 to +0.49 Pg C yr-1 . Despite the significant correlation between the CO2 flux and El Niño indices, we find that El Niño events have, when globally integrated, both enhanced and weakened terrestrial sink strength, with no consistent response across events.