Challenging terrestrial biosphere models with data from the long-term multifactor Prairie Heating and CO2 Enrichment experiment

Annual precipitation explains variability in dryland vegetation greenness globally but not locally

Dryland vegetation productivity is strongly modulated by water availability. As precipitation patterns and variability are altered by climate change, there is a pressing need to better understand vegetation responses to precipitation variability in these ecologically fragile regions. Here we present a global analysis of dryland sensitivity to annual precipitation variations using long-term records of normalized difference vegetation index (NDVI). We show that while precipitation explains 66% of spatial gradients in NDVI across dryland regions, precipitation only accounts for <26% of temporal NDVI variability over most (>75%) dryland regions. We observed this weaker temporal relative to spatial relationship between NDVI and precipitation across all global drylands. We confirmed this result using three alternative water availability metrics that account for water loss to evaporation, and growing season and precipitation timing. This suggests that predicting vegetation responses to future rainfall using space-for-time substitution will strongly overestimate precipitation control on interannual variability in aboveground growth. We explore multiple mechanisms to explain the discrepancy between spatial and temporal responses and find contributions from multiple factors including local-scale vegetation characteristics, climate and soil properties. Earth system models (ESMs) from the latest Coupled Model Intercomparison Project overestimate the observed vegetation sensitivity to precipitation variability up to threefold, particularly during dry years. Given projections of increasing meteorological drought, ESMs are likely to overestimate the impacts of future drought on dryland vegetation with observations suggesting that dryland vegetation is more resistant to annual precipitation variations than ESMs project.

Using climate-driven leaf phenology and growth to improve predictions of gross primary productivity in North American forests

Forest ecosystems are an important sink for terrestrial carbon sequestration. Hence, accurate modeling of the intra- and interannual variability of forest photosynthetic productivity remains a key objective in global biology. Applying climate-driven leaf phenology and growth in models may improve predictions of the forest gross primary productivity (GPP). We used a dynamic non-structural carbohydrates (NSC) model (FORCCHN2) that couples leaf development and phenology to investigate the relationships among photosynthesis and environmental factors. FORCCHN2 simulates spring and autumn phenological events from heat and chilling, respectively. Leaf area index data from satellites along with climate data estimated localized phenological parameters. NSC limitation, immediate temperature, accumulated heat, and growth potential comprised a daily leaf-growth model. Functionally, leaf growth was decoupled from photosynthesis. Leaf biomass determined overall photosynthetic production. We compared this model with outputs of the other six terrestrial biospheric models and with observations from the North American Carbon Program Site Interim Synthesis in 18 forest sites. This model improved the predicted performance of yearly GPP with a 57%-210% increase in correlation (median) and up to a 102% reduction in biases (median), compared to three prognostic models and three prescribed models. At the North America continental scale, the model predicted the average annual GPP of 7.38 Pg C/year from forest ecosystems during 1985-2016. The results showed an increasing trend of GPP in North America (1.0 Pg C/decade). The inclusion of climate-driven phenology and growth has a significant potential for improving dynamic vegetation models, and promotes a further understanding of the complex relationship between environment and photosynthesis.

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily-monthly) timescales and reduces on longer (seasonal-annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. The relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

Plant profit maximization improves predictions of European forest responses to drought

Knowledge of how water stress impacts the carbon and water cycles is a key uncertainty in terrestrial biosphere models. We tested a new profit maximization model, where photosynthetic uptake of CO₂ is optimally traded against plant hydraulic function, as an alternative to the empirical functions commonly used in models to regulate gas exchange during periods of water stress. We conducted a multi-site evaluation of this model at the ecosystem scale, before and during major droughts in Europe. Additionally, we asked whether the maximum hydraulic conductance in the soil-plant continuum k_max (a key model parameter which is not commonly measured) could be predicted from long-term site climate. Compared with a control model with an empirical soil moisture function, the profit maximization model improved the simulation of evapotranspiration during the growing season, reducing the normalized mean square error by c. 63%, across mesic and xeric sites. We also showed that k_max could be estimated from long-term climate, with improvements in the simulation of evapotranspiration at eight out of the 10 forest sites during drought. Although the generalization of this approach is contingent upon determining k_max , it presents a mechanistic trait-based alternative to regulate canopy gas exchange in global models.

Towards better representations of carbon allocation in vegetation: a conceptual framework and mathematical tool

The representation of carbon allocation (CA) in ecosystem differs tremendously among models, resulting in diverse responses of carbon cycling and storage to global change. Several studies have highlighted discrepancies between empirical observations and model predictions, attributing these differences to problems of model structure. We analyzed the mathematical representation of CA in models using concepts from dynamical systems theory; we reviewed a representative sample of models of CA in vegetation and developed a model database within the Python package bgc-md. We asked whether these representations can be generalized as a linear system, or whether a more general framework is needed to accommodate nonlinearities. Some of the vegetation systems simulated with the reviewed models have a fixed partitioning of photosynthetic products, independent of environmental forcing. Vegetation is often represented as a linear system without storage compartments. Yet, other structures with nonlinearities have also been proposed, with important consequences on the temporal trajectories of ecosystem carbon compartments. The proposed mathematical framework unifies the representation of alternative CA schemes, facilitating their classification according to mathematical properties as well as their potential temporal behaviour. It can represent complex processes in a compact form, which can potentially facilitate dialog among empiricists, theoreticians, and modellers.

Desiccation time during drought is highly predictable across species of Eucalyptus from contrasting climates

Catastrophic failure of the water transport pathway in trees is a principal mechanism of mortality during extreme drought. To be able to predict the probability of mortality at an individual and landscape scale we need knowledge of the time for plants to reach critical levels of hydraulic failure. We grew plants of eight species of Eucalyptus originating from contrasting climates before allowing a subset to dehydrate. We tested whether a trait-based model of time to plant desiccation t_crit , from stomatal closure g_s90 to a critical level of hydraulic dysfunction Ψ_crit is consistent with observed dry-down times. Plant desiccation time varied among species, ranging from 96.2 to 332 h at a vapour-pressure deficit of 1 kPa, and was highly predictable using the t_crit model in conjunction with a leaf shedding function. Plant desiccation time was longest in species with high cavitation resistance, strong vulnerability segmentation, wide stomatal-hydraulic safety, and a high ratio of total plant water content to leaf area. Knowledge of t_crit in combination with water-use traits that influence stomatal closure could significantly increase our ability to predict the timing of drought-induced mortality at tree and forest scales.

Joint structural and physiological control on the interannual variation in productivity in a temperate grassland: A data-model comparison

Given the important contributions of semiarid region to global land carbon cycle, accurate modeling of the interannual variability (IAV) of terrestrial gross primary productivity (GPP) is important but remains challenging. By decomposing GPP into leaf area index (LAI) and photosynthesis per leaf area (i.e., GPP_leaf), we investigated the IAV of GPP and the mechanisms responsible in a temperate grassland of northwestern China. We further assessed six ecosystem models for their capabilities in reproducing the observed IAV of GPP in a temperate grassland from 2004 to 2011 in China. We observed that the responses to LAI and GPP_leaf to soil water significantly contributed to IAV of GPP at the grassland ecosystem. Two of six models with prescribed LAI simulated of the observed IAV of GPP quite well, but still underestimated the variance of GPP_leaf, therefore the variance of GPP. In comparison, simulated pattern by the other four models with prognostic LAI differed significantly from the observed IAV of GPP. Only some models with prognostic LAI can capture the observed sharp decline of GPP in drought years. Further analysis indicated that accurately representing the responses of GPP_leaf and leaf stomatal conductance to soil moisture are critical for the models to reproduce the observed IAV of GPP_leaf. Our framework also identified that the contributions of LAI and GPP_leaf to the observed IAV of GPP were relatively independent. We conclude that our framework of decomposing GPP into LAI and GPP_leaf has a significant potential for facilitating future model intercomparison, benchmarking and optimization should be adopted for future data-model comparisons.