When do plant hydraulics matter in terrestrial biosphere modelling?

The ascent of water from the soil to the leaves of vascular plants, described by the study of plant hydraulics, regulates ecosystem responses to environmental forcing and recovery from stress periods. Several approaches to model plant hydraulics have been proposed. In this study, we introduce four different versions of plant hydraulics representations in the terrestrial biosphere model T&C to understand the significance of plant hydraulics to ecosystem functioning. We tested representations of plant hydraulics, investigating plant water capacitance, and long-term xylem damages following drought. The four models we tested were a combination of representations including or neglecting capacitance and including or neglecting xylem damage legacies. Using the models at six case studies spanning semiarid to tropical ecosystems, we quantify how plant xylem flow, plant water storage and long-term xylem damage can modulate overall water and carbon dynamics across multiple time scales. We show that as drought develops, models with plant hydraulics predict a slower onset of plant water stress, and a diurnal variability of water and carbon fluxes closer to observations. Plant water storage was found to be particularly important for the diurnal dynamics of water and carbon fluxes, with models that include plant water capacitance yielding better results. Models including permanent damage to conducting plant tissues show an additional significant drought legacy effect, limiting plant productivity during the recovery phase following major droughts. However, when considering ecosystem responses to the observed climate variability, plant hydraulic modules alone cannot significantly improve the overall model performance, even though they reproduce more realistic water and carbon dynamics. This opens new avenues for model development, explicitly linking plant hydraulics with additional ecosystem processes, such as plant phenology and improved carbon allocation algorithms.

Machine learning for accelerating process-based computation of land biogeochemical cycles.

Global change ecology nowadays embraces ever-growing large observational datasets (big-data) and complex mathematical models that track hundreds of ecological processes (big-model). The rapid advancement of the big-data-big-model has reached its bottleneck: high computational requirements prevent further development of models that need to be integrated over long time-scales to simulate the distribution of ecosystems carbon and nutrient pools and fluxes. Here, we introduce a machine-learning acceleration (MLA) tool to tackle this grand challenge. We focus on the most resource-consuming step in terrestrial biosphere models (TBMs): the equilibration of biogeochemical cycles (spin-up), a prerequisite that can take up to 98% of the computational time. Through three members of the ORCHIDEE TBM family part of the IPSL Earth System Model, including versions that describe the complex interactions between nitrogen, phosphorus and carbon that do not have any analytical solution for the spin-up, we show that an unoptimized MLA reduced the computation demand by 77%-80% for global studies via interpolating the equilibrated state of biogeochemical variables for a subset of model pixels. Despite small biases in the MLA-derived equilibrium, the resulting impact on the predicted regional carbon balance over recent decades is minor. We expect a one-order of magnitude lower computation demand by optimizing the choices of machine learning algorithms, their settings, and balancing the trade-off between quality of MLA predictions and need for TBM simulations for training data generation and bias reduction. Our tool is agnostic to gridded models (beyond TBMs), compatible with existing spin-up acceleration procedures, and opens the door to a wide variety of future applications, with complex non-linear models benefit most from the computational efficiency.

Environmental controls on the light use efficiency of terrestrial gross primary production.

Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf-level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy-level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C3 species; however, additional assumptions are required to "scale up" from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed-effects model. The same model design was fitted to site-based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple "goodness of fit" comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

Assessing dynamic vegetation model parameter uncertainty across Alaskan arctic tundra plant communities.

As the Arctic region moves into uncharted territory under a warming climate, it is important to refine the terrestrial biosphere models (TBMs) that help us understand and predict change. One fundamental uncertainty in TBMs relates to model parameters, configuration variables internal to the model whose value can be estimated from data. We incorporate a version of the Terrestrial Ecosystem Model (TEM) developed for arctic ecosystems into the Predictive Ecosystem Analyzer (PEcAn) framework. PEcAn treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We examined how variation in 21 parameters in the equation for gross primary production influenced model sensitivity and uncertainty in terms of two carbon fluxes (net primary productivity and heterotrophic respiration) and two carbon (C) pools (vegetation C and soil C). We set up different parameterizations of TEM across a range of tundra types (tussock tundra, heath tundra, wet sedge tundra, and shrub tundra) in northern Alaska, along a latitudinal transect extending from the coastal plain near Utqiagvik to the southern foothills of the Brooks Range, to the Seward Peninsula. TEM was most sensitive to parameters related to the temperature regulation of photosynthesis. Model uncertainty was mostly due to parameters related to leaf area, temperature regulation of photosynthesis, and the stomatal responses to ambient light conditions. Our analysis also showed that sensitivity and uncertainty to a given parameter varied spatially. At some sites, model sensitivity and uncertainty tended to be connected to a wider range of parameters, underlining the importance of assessing tundra community processes across environmental gradients or geographic locations. Generally, across sites, the flux of net primary productivity (NPP) and pool of vegetation C had about equal uncertainty, while heterotrophic respiration had higher uncertainty than the pool of soil C. Our study illustrates the complexity inherent in evaluating parameter uncertainty across highly heterogeneous arctic tundra plant communities. It also provides a framework for iteratively testing how newly collected field data related to key parameters may result in more effective forecasting of Arctic change.

Multi-hypothesis comparison of Farquhar and Collatz photosynthesis models reveals the unexpected influence of empirical assumptions at leaf and global scales.

Mechanistic photosynthesis models are at the heart of terrestrial biosphere models (TBMs) simulating the daily, monthly, annual and decadal rhythms of carbon assimilation (A). These models are founded on robust mathematical hypotheses that describe how A responds to changes in light and atmospheric CO2 concentration. Two predominant photosynthesis models are in common usage: Farquhar (FvCB) and Collatz (CBGB). However, a detailed quantitative comparison of these two models has never been undertaken. In this study, we unify the FvCB and CBGB models to a common parameter set and use novel multi-hypothesis methods (that account for both hypothesis and parameter variability) for process-level sensitivity analysis. These models represent three key biological processes: carboxylation, electron transport, triose phosphate use (TPU) and an additional model process: limiting-rate selection. Each of the four processes comprises 1-3 alternative hypotheses giving 12 possible individual models with a total of 14 parameters. To broaden inference, TBM simulations were run and novel, high-resolution photosynthesis measurements were made. We show that parameters associated with carboxylation are the most influential parameters but also reveal the surprising and marked dominance of the limiting-rate selection process (accounting for 57% of the variation in A vs. 22% for carboxylation). The limiting-rate selection assumption proposed by CBGB smooths the transition between limiting rates and always reduces A below the minimum of all potentially limiting rates, by up to 25%, effectively imposing a fourth limitation on A. Evaluation of the CBGB smoothing function in three TBMs demonstrated a reduction in global A by 4%-10%, equivalent to 50%-160% of current annual fossil fuel emissions. This analysis reveals a surprising and previously unquantified influence of a process that has been integral to many TBMs for decades, highlighting the value of multi-hypothesis methods.

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

An optimality-based model explains seasonal variation in C3 plant photosynthetic capacity.

The maximum rate of carboxylation (Vcmax ) is an essential leaf trait determining the photosynthetic capacity of plants. Existing approaches for estimating Vcmax at large scale mainly rely on empirical relationships with proxies such as leaf nitrogen/chlorophyll content or hyperspectral reflectance, or on complicated inverse models from gross primary production or solar-induced fluorescence. A novel mechanistic approach based on the assumption that plants optimize resource investment coordinating with environment and growth has been shown to accurately predict C3 plant Vcmax based on mean growing season environmental conditions. However, the ability of optimality theory to explain seasonal variation in Vcmax has not been fully investigated. Here, we adapt an optimality-based model to simulate daily Vcmax,25C (Vcmax at a standardized temperature of 25°C) by incorporating the effects of antecedent environment, which affects current plant functioning, and dynamic light absorption, which coordinates with plant functioning. We then use seasonal Vcmax,25C field measurements from 10 sites across diverse ecosystems to evaluate model performance. Overall, the model explains about 83% of the seasonal variation in C3 plant Vcmax,25C across the 10 sites, with a medium root mean square error of 12.3 µmol m-2  s-1 , which suggests that seasonal changes in Vcmax,25C are consistent with optimal plant function. We show that failing to account for acclimation to antecedent environment or coordination with dynamic light absorption dramatically decreases estimation accuracy. Our results show that optimality-based approach can accurately reproduce seasonal variation in canopy photosynthetic potential, and suggest that incorporating such theory into next-generation trait-based terrestrial biosphere models would improve predictions of global photosynthesis.

Contrasting Regional Carbon Cycle Responses to Seasonal Climate Anomalies Across the East-West Divide of Temperate North America.

Across temperate North America, interannual variability (IAV) in gross primary production (GPP) and net ecosystem exchange (NEE) and their relationship with environmental drivers are poorly understood. Here, we examine IAV in GPP and NEE and their relationship to environmental drivers using two state-of-the-science flux products: NEE constrained by surface and space-based atmospheric CO2 measurements over 2010-2015 and satellite up-scaled GPP from FluxSat over 2001-2017. We show that the arid western half of temperate North America provides a larger contribution to IAV in GPP (104% of east) and NEE (127% of east) than the eastern half, in spite of smaller magnitude of annual mean GPP and NEE. This occurs because anomalies in western ecosystems are temporally coherent across the growing season leading to an amplification of GPP and NEE. In contrast, IAV in GPP and NEE in eastern ecosystems is dominated by seasonal compensation effects, associated with opposite responses to temperature anomalies in spring and summer. Terrestrial biosphere models in the MsTMIP ensemble generally capture these differences between eastern and western temperate North America, although there is considerable spread between models.

How are nitrogen availability, fine-root mass, and nitrogen uptake related empirically? Implications for models and theory.

Understanding the effects of global change in terrestrial communities requires an understanding of how limiting resources interact with plant traits to affect productivity. Here, we focus on nitrogen and ask whether plant community nitrogen uptake rate is determined (a) by nitrogen availability alone or (b) by the product of nitrogen availability and fine-root mass. Surprisingly, this is not empirically resolved. We performed controlled microcosm experiments and reanalyzed published pot experiments and field data to determine the relationship between community-level nitrogen uptake rate, nitrogen availability, and fine-root mass for 46 unique combinations of species, nitrogen levels, and growing conditions. We found that plant community nitrogen uptake rate was unaffected by fine-root mass in 63% of cases and saturated with fine-root mass in 29% of cases (92% in total). In contrast, plant community nitrogen uptake rate was clearly affected by nitrogen availability. The results support the idea that although plants may over-proliferate fine roots for individual-level competition, it comes without an increase in community-level nitrogen uptake. The results have implications for the mechanisms included in coupled carbon-nitrogen terrestrial biosphere models (CN-TBMs) and are consistent with CN-TBMs that operate above the individual scale and omit fine-root mass in equations of nitrogen uptake rate but inconsistent with the majority of CN-TBMs, which operate above the individual scale and include fine-root mass in equations of nitrogen uptake rate. For the much smaller number of CN-TBMs that explicitly model individual-based belowground competition for nitrogen, the results suggest that the relative (not absolute) fine-root mass of competing individuals should be included in the equations that determine individual-level nitrogen uptake rates. By providing empirical data to support the assumptions used in CN-TBMs, we put their global climate change predictions on firmer ground.

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity.

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). We evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. These results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

New insights on plant phenological response to temperature revealed from long-term widespread observations in China.

Constraints of temperature on spring plant phenology are closely related to plant growth, vegetation dynamics, and ecosystem carbon cycle. However, the effects of temperature on leaf onset, especially for winter chilling, are still not well understood. Using long-term, widespread in situ phenology observations collected over China for multiple plant species, this study analyzes the quantitative response of leaf onset to temperature, and compares empirical findings with existing theories and modeling approaches, as implemented in 18 phenology algorithms. Results show that the growing degree days (GDD) required for leaf onset vary distinctly among plant species and geographical locations as well as at organizational levels (species and community), pointing to diverse adaptation strategies. Chilling durations (CHD) needed for releasing bud dormancy decline monotonously from cold to warm areas with very limited interspecies variations. Results also reveal that winter chilling is a crucial component of phenology models, and its effect is better captured with an index that accounts for the inhomogeneous effectiveness of low temperature to chilling rate than with the conventional CHD index. The impact of spring warming on leaf onset is nonlinear, better represented by a logistical function of temperature than by the linear function currently implemented in biosphere models. The optimized base temperatures for thermal accumulation and the optimal chilling temperatures are species-dependent and average at 6.9 and 0.2°C, respectively. Overall, plants' chilling requirement is not a constant, and more chilling generally results in less requirement of thermal accumulation for leaf onset. Our results clearly demonstrate multiple deficiencies of the parameters (e.g., base temperature) and algorithms (e.g., method for calculating GDD) in conventional phenology models to represent leaf onset. Therefore, this study not only advances our mechanistic and quantitative understanding of temperature controls on leaf onset but also provides critical information for improving existing phenology models.

Limited effect of ozone reductions on the 20-year photosynthesis trend at Harvard forest.

Ozone (O3 ) damage to leaves can reduce plant photosynthesis, which suggests that declines in ambient O3 concentrations ([O3 ]) in the United States may have helped increase gross primary production (GPP) in recent decades. Here, we assess the effect of long-term changes in ambient [O3 ] using 20 years of observations at Harvard forest. Using artificial neural networks, we found that the effect of the inclusion of [O3 ] as a predictor was slight, and independent of O3 concentrations, which suggests limited high-frequency O3 inhibition of GPP at this site. Simulations with a terrestrial biosphere model, however, suggest an average long-term O3 inhibition of 10.4% for 1992-2011. A decline of [O3 ] over the measurement period resulted in moderate predicted GPP trends of 0.02-0.04 µmol C m-2  s-1  yr-1 , which is negligible relative to the total observed GPP trend of 0.41 µmol C m-2  s-1  yr-1 . A similar conclusion is achieved with the widely used AOT40 metric. Combined, our results suggest that ozone reductions at Harvard forest are unlikely to have had a large impact on the photosynthesis trend over the past 20 years. Such limited effects are mainly related to the slow responses of photosynthesis to changes in [O3 ]. Furthermore, we estimate that 40% of photosynthesis happens in the shade, where stomatal conductance and thus [O3 ] deposition is lower than for sunlit leaves. This portion of GPP remains unaffected by [O3 ], thus helping to buffer the changes of total photosynthesis due to varied [O3 ]. Our analyses suggest that current ozone reductions, although significant, cannot substantially alleviate the damages to forest ecosystems.

The fate of Amazonian ecosystems over the coming century arising from changes in climate, atmospheric CO2, and land use.

There is considerable interest in understanding the fate of the Amazon over the coming century in the face of climate change, rising atmospheric CO2 levels, ongoing land transformation, and changing fire regimes within the region. In this analysis, we explore the fate of Amazonian ecosystems under the combined impact of these four environmental forcings using three terrestrial biosphere models (ED2, IBIS, and JULES) forced by three bias-corrected IPCC AR4 climate projections (PCM1, CCSM3, and HadCM3) under two land-use change scenarios. We assess the relative roles of climate change, CO2 fertilization, land-use change, and fire in driving the projected changes in Amazonian biomass and forest extent. Our results indicate that the impacts of climate change are primarily determined by the direction and severity of projected changes in regional precipitation: under the driest climate projection, climate change alone is predicted to reduce Amazonian forest cover by an average of 14%. However, the models predict that CO2 fertilization will enhance vegetation productivity and alleviate climate-induced increases in plant water stress, and, as a result, sustain high biomass forests, even under the driest climate scenario. Land-use change and climate-driven changes in fire frequency are predicted to cause additional aboveground biomass loss and reductions in forest extent. The relative impact of land use and fire dynamics compared to climate and CO2 impacts varies considerably, depending on both the climate and land-use scenario, and on the terrestrial biosphere model used, highlighting the importance of improved quantitative understanding of all four factors - climate change, CO2 fertilization effects, fire, and land use - to the fate of the Amazon over the coming century.

Confronting model predictions of carbon fluxes with measurements of Amazon forests subjected to experimental drought.

Considerable uncertainty surrounds the fate of Amazon rainforests in response to climate change. Here, carbon (C) flux predictions of five terrestrial biosphere models (Community Land Model version 3.5 (CLM3.5), Ecosystem Demography model version 2.1 (ED2), Integrated BIosphere Simulator version 2.6.4 (IBIS), Joint UK Land Environment Simulator version 2.1 (JULES) and Simple Biosphere model version 3 (SiB3)) and a hydrodynamic terrestrial ecosystem model (the Soil-Plant-Atmosphere (SPA) model) were evaluated against measurements from two large-scale Amazon drought experiments. Model predictions agreed with the observed C fluxes in the control plots of both experiments, but poorly replicated the responses to the drought treatments. Most notably, with the exception of ED2, the models predicted negligible reductions in aboveground biomass in response to the drought treatments, which was in contrast to an observed c. 20% reduction at both sites. For ED2, the timing of the decline in aboveground biomass was accurate, but the magnitude was too high for one site and too low for the other. Three key findings indicate critical areas for future research and model development. First, the models predicted declines in autotrophic respiration under prolonged drought in contrast to measured increases at one of the sites. Secondly, models lacking a phenological response to drought introduced bias in the sensitivity of canopy productivity and respiration to drought. Thirdly, the phenomenological water-stress functions used by the terrestrial biosphere models to represent the effects of soil moisture on stomatal conductance yielded unrealistic diurnal and seasonal responses to drought.

The Grain-for-Green project offsets warming-induced soil organic carbon loss and increases soil carbon stock in Chinese Loess Plateau.

The dynamics of soil organic carbon (SOC) stock is a vital element affecting the climate, and ecological restoration is potentially an effective measure to mitigate climate change by enhancing vegetation and soil carbon stocks and thereby offsetting greenhouse gas emissions. The Grain-for-Green project (GFGP) implemented in Chinese Loess Plateau (LP) since 1999 is one of the largest ecological restoration projects in the world. However, the contributions of ecological restoration and climate change to ecosystem soil carbon sequestration are still unclear. In this study, we improved a soil carbon decomposition framework by optimizing the initial SOC stock based on full spatial simulation of SOC and incorporating the priming effect to investigate the SOC dynamics across the LP GFGP region from 1982 through 2017. Our results indicated that SOC stock in the GFGP region increased by 20.18 Tg C from 1982 through 2017. Most portion (15.83 Tg C) of the SOC increase was accumulated when the GFGP was initiated, with a SOC sink of 16.12 Tg C owing to revegetation restoration and a carbon loss of 0.29 Tg C due to warming during this period. The relationships between SOC and forest canopy height and investigations on the SOC dynamics after afforestation revealed that the accumulation rate of SOC could be as high as 24.68 g C m-2 yr-1 during the 70 years following afforestation, and that SOC could decline thereafter (-8.89 g C m-2 yr-1), which was mainly caused by warming. This study provides a new method for quantifying the contribution of ecological restoration to SOC changes, and also cautions the potential risk of LP SOC loss in the mature forest soil under future warming.

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.