Climate-ecosystem modelling made easy: The Land Sites Platform

Dynamic Global Vegetation Models (DGVMs) provide a state-of-the-art process-based approach to study the complex interplay between vegetation and its physical environment. For example, they help to predict how terrestrial plants interact with climate, soils, disturbance and competition for resources. We argue that there is untapped potential for the use of DGVMs in ecological and ecophysiological research. One fundamental barrier to realize this potential is that many researchers with relevant expertize (ecology, plant physiology, soil science, etc.) lack access to the technical resources or awareness of the research potential of DGVMs. Here we present the Land Sites Platform (LSP): new software that facilitates single-site simulations with the Functionally Assembled Terrestrial Ecosystem Simulator, an advanced DGVM coupled with the Community Land Model. The LSP includes a Graphical User Interface and an Application Programming Interface, which improve the user experience and lower the technical thresholds for installing these model architectures and setting up model experiments. The software is distributed via version-controlled containers; researchers and students can run simulations directly on their personal computers or servers, with relatively low hardware requirements, and on different operating systems. Version 1.0 of the LSP supports site-level simulations. We provide input data for 20 established geo-ecological observation sites in Norway and workflows to add generic sites from public global datasets. The LSP makes standard model experiments with default data easily achievable (e.g., for educational or introductory purposes) while retaining flexibility for more advanced scientific uses. We further provide tools to visualize the model input and output, including simple examples to relate predictions to local observations. The LSP improves access to land surface and DGVM modelling as a building block of community cyberinfrastructure that may inspire new avenues for mechanistic ecosystem research across disciplines.

Plant water use theory should incorporate hypotheses about extreme environments, population ecology, and community ecology

Plant water use theory has largely been developed within a plant-performance paradigm that conceptualizes water use in terms of value for carbon gain and that sits within a neoclassical economic framework. This theory works very well in many contexts but does not consider other values of water to plants that could impact their fitness. Here, we survey a range of alternative hypotheses for drivers of water use and stomatal regulation. These hypotheses are organized around relevance to extreme environments, population ecology, and community ecology. Most of these hypotheses are not yet empirically tested and some are controversial (e.g. requiring more agency and behavior than is commonly believed possible for plants). Some hypotheses, especially those focused around using water to avoid thermal stress, using water to promote reproduction instead of growth, and using water to hoard it, may be useful to incorporate into theory or to implement in Earth System Models.

Soil respiration-driven CO2 pulses dominate Australia's flux variability

The Australian continent contributes substantially to the year-to-year variability of the global terrestrial carbon dioxide (CO₂) sink. However, the scarcity of in situ observations in remote areas prevents the deciphering of processes that force the CO₂ flux variability. In this study, by examining atmospheric CO₂ measurements from satellites in the period 2009-2018, we find recurrent end-of-dry-season CO₂ pulses over the Australian continent. These pulses largely control the year-to-year variability of Australia's CO₂ balance. They cause two to three times larger seasonal variations compared with previous top-down inversions and bottom-up estimates. The pulses occur shortly after the onset of rainfall and are driven by enhanced soil respiration preceding photosynthetic uptake in Australia's semiarid regions. The suggested continental-scale relevance of soil-rewetting processes has substantial implications for our understanding and modeling of global climate-carbon cycle feedbacks.

Physiological trait networks enhance understanding of crop growth and water use in contrasting environments

Plant function arises from a complex network of structural and physiological traits. Explicit representation of these traits, as well as their connections with other biophysical processes, is required to advance our understanding of plant-soil-climate interactions. We used the Terrestrial Regional Ecosystem Exchange Simulator (TREES) to evaluate physiological trait networks in maize. Net primary productivity (NPP) and grain yield were simulated across five contrasting climate scenarios. Simulations achieving high NPP and grain yield in high precipitation environments featured trait networks conferring high water use strategies: deep roots, high stomatal conductance at low water potential ("risky" stomatal regulation), high xylem hydraulic conductivity and high maximal leaf area index. In contrast, high NPP and grain yield was achieved in dry environments with low late-season precipitation via water conserving trait networks: deep roots, high embolism resistance and low stomatal conductance at low leaf water potential ("conservative" stomatal regulation). We suggest that our approach, which allows for the simultaneous evaluation of physiological traits, soil characteristics and their interactions (i.e., networks), has potential to improve our understanding of crop performance in different environments. In contrast, evaluating single traits in isolation of other coordinated traits does not appear to be an effective strategy for predicting plant performance.

Process-oriented analysis of dominant sources of uncertainty in the land carbon sink

The observed global net land carbon sink is captured by current land models. All models agree that atmospheric CO₂ and nitrogen deposition driven gains in carbon stocks are partially offset by climate and land-use and land-cover change (LULCC) losses. However, there is a lack of consensus in the partitioning of the sink between vegetation and soil, where models do not even agree on the direction of change in carbon stocks over the past 60 years. This uncertainty is driven by plant productivity, allocation, and turnover response to atmospheric CO₂ (and to a smaller extent to LULCC), and the response of soil to LULCC (and to a lesser extent climate). Overall, differences in turnover explain ~70% of model spread in both vegetation and soil carbon changes. Further analysis of internal plant and soil (individual pools) cycling is needed to reduce uncertainty in the controlling processes behind the global land carbon sink.

The global net land sink is relatively well constrained. However, the responsible drivers and above/below-ground partitioning are highly uncertain. Model issues regarding turnover of individual plant and soil components are responsible.

Informing Nature-based Climate Solutions for the United States with the best-available science

Nature-based Climate Solutions (NbCS) are managed alterations to ecosystems designed to increase carbon sequestration or reduce greenhouse gas emissions. While they have growing public and private support, the realizable benefits and unintended consequences of NbCS are not well understood. At regional scales where policy decisions are often made, NbCS benefits are estimated from soil and tree survey data that can miss important carbon sources and sinks within an ecosystem, and do not reveal the biophysical impacts of NbCS for local water and energy cycles. The only direct observations of ecosystem-scale carbon fluxes, for example, by eddy covariance flux towers, have not yet been systematically assessed for what they can tell us about NbCS potentials, and state-of-the-art remote sensing products and land-surface models are not yet being widely used to inform NbCS policymaking or implementation. As a result, there is a critical mismatch between the point- and tree-scale data most often used to assess NbCS benefits and impacts, the ecosystem and landscape scales where NbCS projects are implemented, and the regional to continental scales most relevant to policymaking. Here, we propose a research agenda to confront these gaps using data and tools that have long been used to understand the mechanisms driving ecosystem carbon and energy cycling, but have not yet been widely applied to NbCS. We outline steps for creating robust NbCS assessments at both local to regional scales that are informed by ecosystem-scale observations, and which consider concurrent biophysical impacts, future climate feedbacks, and the need for equitable and inclusive NbCS implementation strategies. We contend that these research goals can largely be accomplished by shifting the scales at which pre-existing tools are applied and blended together, although we also highlight some opportunities for more radical shifts in approach.

Representation of Leaf-to-Canopy Radiative Transfer Processes Improves Simulation of Far-Red Solar-Induced Chlorophyll Fluorescence in the Community Land Model Version 5

Recent advances in satellite observations of solar-induced chlorophyll fluorescence (SIF) provide a new opportunity to constrain the simulation of terrestrial gross primary productivity (GPP). Accurate representation of the processes driving SIF emission and its radiative transfer to remote sensing sensors is an essential prerequisite for data assimilation. Recently, SIF simulations have been incorporated into several land surface models, but the scaling of SIF from leaf-level to canopy-level is usually not well-represented. Here, we incorporate the simulation of far-red SIF observed at nadir into the Community Land Model version 5 (CLM5). Leaf-level fluorescence yield was simulated by a parametric simplification of the Soil Canopy-Observation of Photosynthesis and Energy fluxes model (SCOPE). And an efficient and accurate method based on escape probability is developed to scale SIF from leaf-level to top-of-canopy while taking clumping and the radiative transfer processes into account. SIF simulated by CLM5 and SCOPE agreed well at sites except one in needleleaf forest (R ² > 0.91, root-mean-square error <0.19 W⋅m^-2⋅sr^-1⋅μm^-1), and captured the day-to-day variation of tower-measured SIF at temperate forest sites (R ² > 0.68). At the global scale, simulated SIF generally captured the spatial and seasonal patterns of satellite-observed SIF. Factors including the fluorescence emission model, clumping, bidirectional effect, and leaf optical properties had considerable impacts on SIF simulation, and the discrepancies between simulate d and observed SIF varied with plant functional type. By improving the representation of radiative transfer for SIF simulation, our model allows better comparisons between simulated and observed SIF toward constraining GPP simulations.

Legislative and functional aspects of different metrics used for ozone risk assessment to forests

Surface ozone (O₃) is a threat to forests by decreasing photosynthesis and, consequently, influencing the strength of land carbon sink. However, due to the lack of continuous surface O₃ measurements, observational-based assessments of O₃ impacts on forests are largely missing at hemispheric to global scales. Currently, some metrics are used for regulatory purposes by governments or national agencies to protect forests against the negative impacts of ozone: in particular, both Europe and United States (US) makes use of two different exposure-based metrics, i.e. AOT40 and W126, respectively. However, because of some limitations in these metrics, a new standard is under consideration by the European Union (EU) to replace the current exposure metric. We analyse here the different air quality standards set or proposed for use in Europe and in the US to protect forests from O₃ and to evaluate their spatial and temporal consistency while assessing their effectiveness in protecting northern-hemisphere forests. Then, we compare their results with the information obtained from a complex land surface model (ORCHIDEE). We find that present O₃ uptake decreases gross primary production (GPP) in 37.7% of the NH forested area of northern hemisphere with a mean loss of 2.4% year^-1. We show how the proposed US (W126) and the currently used European (AOT40) air quality standards substantially overestimate the extension of potential vulnerable regions, predicting that 46% and 61% of the Northern Hemisphere (NH) forested area are at risk of O₃ pollution. Conversely, the new proposed European standard (POD1) identifies lower extension of vulnerability regions (39.6%).

Response to Comments on "Recent global decline of CO2 fertilization effects on vegetation photosynthesis"

[Figure: see text].

Greening drylands despite warming consistent with carbon dioxide fertilization effect

The rising atmospheric CO₂ concentration leads to a CO₂ fertilization effect on plants-that is, increased photosynthetic uptake of CO₂ by leaves and enhanced water-use efficiency (WUE). Yet, the resulting net impact of CO₂ fertilization on plant growth and soil moisture (SM) savings at large scale is poorly understood. Drylands provide a natural experimental setting to detect the CO₂ fertilization effect on plant growth since foliage amount, plant water-use and photosynthesis are all tightly coupled in water-limited ecosystems. A long-term change in the response of leaf area index (LAI, a measure of foliage amount) to changes in SM is likely to stem from changing water demand of primary productivity in water-limited ecosystems and is a proxy for changes in WUE. Using 34-year satellite observations of LAI and SM over tropical and subtropical drylands, we identify that a 1% increment in SM leads to 0.15% (±0.008, 95% confidence interval) and 0.51% (±0.01, 95% confidence interval) increments in LAI during 1982-1998 and 1999-2015, respectively. The increasing response of LAI to SM has contributed 7.2% (±3.0%, 95% confidence interval) to total dryland greening during 1999-2015 compared to 1982-1998. The increasing response of LAI to SM is consistent with the CO₂ fertilization effect on WUE in water-limited ecosystems, indicating that a given amount of SM has sustained greater amounts of photosynthetic foliage over time. The LAI responses to changes in SM from seven dynamic global vegetation models are not always consistent with observations, highlighting the need for improved process knowledge of terrestrial ecosystem responses to rising atmospheric CO₂ concentration.

Solar geoengineering can alleviate climate change pressures on crop yields

Solar geoengineering (SG) and CO₂ emissions reduction can each alleviate anthropogenic climate change, but their impacts on food security are not yet fully understood. Using an advanced crop model within an Earth system model, we analysed the yield responses of six major crops to three SG technologies (SGs) and emissions reduction when they provide roughly the same reduction in radiative forcing and assume the same land use. We found sharply distinct yield responses to changes in radiation, moisture and CO₂, but comparable significant cooling benefits for crop yields by all four methods. Overall, global yields increase ~10% under the three SGs and decrease 5% under emissions reduction, the latter primarily due to reduced CO₂ fertilization, relative to business as usual by the late twenty-first century. Relative humidity dominates the hydrological effect on yields of rainfed crops, with little contribution from precipitation. The net insolation effect is negligible across all SGs, contrary to previous findings.

Climate-Driven Variability and Trends in Plant Productivity Over Recent Decades Based on Three Global Products

Variability in climate exerts a strong influence on vegetation productivity (gross primary productivity; GPP), and therefore has a large impact on the land carbon sink. However, no direct observations of global GPP exist, and estimates rely on models that are constrained by observations at various spatial and temporal scales. Here, we assess the consistency in GPP from global products which extend for more than three decades; two observation-based approaches, the upscaling of FLUXNET site observations (FLUXCOM) and a remote sensing derived light use efficiency model (RS-LUE), and from a suite of terrestrial biosphere models (TRENDYv6). At local scales, we find high correlations in annual GPP among the products, with exceptions in tropical and high northern latitudes. On longer time scales, the products agree on the direction of trends over 58% of the land, with large increases across northern latitudes driven by warming trends. Further, tropical regions exhibit the largest interannual variability in GPP, with both rainforests and savannas contributing substantially. Variability in savanna GPP is likely predominantly driven by water availability, although temperature could play a role via soil moisture-atmosphere feedbacks. There is, however, no consensus on the magnitude and driver of variability of tropical forests, which suggest uncertainties in process representations and underlying observations remain. These results emphasize the need for more direct long-term observations of GPP along with an extension of in situ networks in underrepresented regions (e.g., tropical forests). Such capabilities would support efforts to better validate relevant processes in models, to more accurately estimate GPP.

Causes of slowing-down seasonal CO2 amplitude at Mauna Loa

Changing amplitude of the seasonal cycle of atmospheric CO₂ (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO₂ record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing-down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing-down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing-down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO₂ fertilization. This challenges the view that CO₂ fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land-atmosphere CO₂ flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.

Enhanced regional terrestrial carbon uptake over Korea revealed by atmospheric CO2 measurements from 1999 to 2017

Understanding changes in terrestrial carbon balance is important to improve our knowledge of the regional carbon cycle and climate change. However, evaluating regional changes in the terrestrial carbon balance is challenging due to the lack of surface flux measurements. This study reveals that the terrestrial carbon uptake over the Republic of Korea has been enhanced from 1999 to 2017 by analyzing long-term atmospheric CO₂ concentration measurements at the Anmyeondo Station (36.53°N, 126.32°E) located in the western coast. The influence of terrestrial carbon flux on atmospheric CO₂ concentrations (ΔCO₂ ) is estimated from the difference of CO₂ concentrations that were influenced by the land sector (through easterly winds) and the Yellow Sea sector (through westerly winds). We find a significant trend in ΔCO₂ of -4.75 ppm per decade (p < .05) during the vegetation growing season (May through October), suggesting that the regional terrestrial carbon uptake has increased relative to the surrounding ocean areas. Combined analysis with satellite measured normalized difference vegetation index and gross primary production shows that the enhanced carbon uptake is associated with significant nationwide increases in vegetation and its production. Process-based terrestrial model and inverse model simulations estimate that regional terrestrial carbon uptake increases by up to 18.9 and 8.0 Tg C for the study period, accounting for 13.4% and 5.7% of the average annual domestic carbon emissions, respectively. Atmospheric chemical transport model simulations indicate that the enhanced terrestrial carbon sink is the primary reason for the observed ΔCO₂ trend rather than anthropogenic emissions and atmospheric circulation changes. Our results highlight the fact that atmospheric CO₂ measurements could open up the possibility of detecting regional changes in the terrestrial carbon cycle even where anthropogenic emissions are not negligible.

State of the science in reconciling top-down and bottom-up approaches for terrestrial CO2 budget

Robust estimates of CO₂ budget, CO₂ exchanged between the atmosphere and terrestrial biosphere, are necessary to better understand the role of the terrestrial biosphere in mitigating anthropogenic CO₂ emissions. Over the past decade, this field of research has advanced through understanding of the differences and similarities of two fundamentally different approaches: "top-down" atmospheric inversions and "bottom-up" biosphere models. Since the first studies were undertaken, these approaches have shown an increasing level of agreement, but disagreements in some regions still persist, in part because they do not estimate the same quantity of atmosphere-biosphere CO₂ exchange. Here, we conducted a thorough comparison of CO₂ budgets at multiple scales and from multiple methods to assess the current state of the science in estimating CO₂ budgets. Our set of atmospheric inversions and biosphere models, which were adjusted for a consistent flux definition, showed a high level of agreement for global and hemispheric CO₂ budgets in the 2000s. Regionally, improved agreement in CO₂ budgets was notable for North America and Southeast Asia. However, large gaps between the two methods remained in East Asia and South America. In other regions, Europe, boreal Asia, Africa, South Asia, and Oceania, it was difficult to determine whether those regions act as a net sink or source because of the large spread in estimates from atmospheric inversions. These results highlight two research directions to improve the robustness of CO₂ budgets: (a) to increase representation of processes in biosphere models that could contribute to fill the budget gaps, such as forest regrowth and forest degradation; and (b) to reduce sink-source compensation between regions (dipoles) in atmospheric inversion so that their estimates become more comparable. Advancements on both research areas will increase the level of agreement between the top-down and bottom-up approaches and yield more robust knowledge of regional CO₂ budgets.

Increased atmospheric vapor pressure deficit reduces global vegetation growth

Atmospheric vapor pressure deficit (VPD) is a critical variable in determining plant photosynthesis. Synthesis of four global climate datasets reveals a sharp increase of VPD after the late 1990s. In response, the vegetation greening trend indicated by a satellite-derived vegetation index (GIMMS3g), which was evident before the late 1990s, was subsequently stalled or reversed. Terrestrial gross primary production derived from two satellite-based models (revised EC-LUE and MODIS) exhibits persistent and widespread decreases after the late 1990s due to increased VPD, which offset the positive CO₂ fertilization effect. Six Earth system models have consistently projected continuous increases of VPD throughout the current century. Our results highlight that the impacts of VPD on vegetation growth should be adequately considered to assess ecosystem responses to future climate conditions.

Impact of the 2015/2016 El Niño on the terrestrial carbon cycle constrained by bottom-up and top-down approaches

Evaluating the response of the land carbon sink to the anomalies in temperature and drought imposed by El Niño events provides insights into the present-day carbon cycle and its climate-driven variability. It is also a necessary step to build confidence in terrestrial ecosystems models' response to the warming and drying stresses expected in the future over many continents, and particularly in the tropics. Here we present an in-depth analysis of the response of the terrestrial carbon cycle to the 2015/2016 El Niño that imposed extreme warming and dry conditions in the tropics and other sensitive regions. First, we provide a synthesis of the spatio-temporal evolution of anomalies in net land–atmosphere CO₂ fluxes estimated by two in situ measurements based on atmospheric inversions and 16 land-surface models (LSMs) from TRENDYv6. Simulated changes in ecosystem productivity, decomposition rates and fire emissions are also investigated. Inversions and LSMs generally agree on the decrease and subsequent recovery of the land sink in response to the onset, peak and demise of El Niño conditions and point to the decreased strength of the land carbon sink: by 0.4–0.7 PgC yr⁻¹ (inversions) and by 1.0 PgC yr⁻¹ (LSMs) during 2015/2016. LSM simulations indicate that a decrease in productivity, rather than increase in respiration, dominated the net biome productivity anomalies in response to ENSO throughout the tropics, mainly associated with prolonged drought conditions.

This article is part of a discussion meeting issue ‘The impact of the 2015/2016 El Niño on the terrestrial tropical carbon cycle: patterns, mechanisms and implications’.

Global vegetation productivity response to climatic oscillations during the satellite era

Climate control on global vegetation productivity patterns has intensified in response to recent global warming. Yet, the contributions of the leading internal climatic variations to global vegetation productivity are poorly understood. Here, we use 30 years of global satellite observations to study climatic variations controls on continental and global vegetation productivity patterns. El Niño-Southern Oscillation (ENSO) phases (La Niña, neutral, and El Niño years) appear to be a weaker control on global-scale vegetation productivity than previously thought, although continental-scale responses are substantial. There is also clear evidence that other non-ENSO climatic variations have a strong control on spatial patterns of vegetation productivity mainly through their influence on temperature. Among the eight leading internal climatic variations, the East Atlantic/West Russia Pattern extensively controls the ensuing year vegetation productivity of the most productive tropical and temperate forest ecosystems of the Earth's vegetated surface through directionally consistent influence on vegetation greenness. The Community Climate System Model (CCSM4) simulations do not capture the observed patterns of vegetation productivity responses to internal climatic variations. Our analyses show the ubiquitous control of climatic variations on vegetation productivity and can further guide CCSM and other Earth system models developments to represent vegetation response patterns to unforced variability. Several winter time internal climatic variation indices show strong potentials on predicting growing season vegetation productivity two to six seasons ahead which enables national governments and farmers forecast crop yield to ensure supplies of affordable food, famine early warning, and plan management options to minimize yield losses ahead of time.