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.

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 (CO2) sink. However, the scarcity of in situ observations in remote areas prevents the deciphering of processes that force the CO2 flux variability. In this study, by examining atmospheric CO2 measurements from satellites in the period 2009-2018, we find recurrent end-of-dry-season CO2 pulses over the Australian continent. These pulses largely control the year-to-year variability of Australia's CO2 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.

Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

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.

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 CO2 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 CO2 (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.

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

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.

Future bioenergy expansion could alter carbon sequestration potential and exacerbate water stress in the United States.

The maximum future projected bioenergy expansion potential, in scenarios limiting warming to 2°C or below, is equivalent to half of present-day croplands. We quantify the impacts of large-scale bioenergy expansion against re/afforestation, which remain elusive, using an integrated human-natural system modeling framework with explicit representation of perennial bioenergy crops. The end-of-century net carbon sequestration due to bioenergy deployment coupled with carbon capture and storage largely depends on fossil fuel displacement types, ranging from 11.4 to 31.2 PgC over the conterminous United States. These net carbon sequestration benefits are inclusive of a 10 PgC carbon release due to land use conversions and a 2.4 PgC loss of additional carbon sink capacity associated with bioenergy-driven deforestation. Moreover, nearly one-fourth of U.S. land areas will suffer severe water stress by 2100 due to either reduced availability or deteriorated quality. These broader impacts of bioenergy expansion should be weighed against the costs and benefits of re/afforestation-based strategies.

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.

Increasing the spatial and temporal impact of ecological research: A roadmap for integrating a novel terrestrial process into an Earth system model.

Terrestrial ecosystems regulate Earth's climate through water, energy, and biogeochemical transformations. Despite a key role in regulating the Earth system, terrestrial ecology has historically been underrepresented in the Earth system models (ESMs) that are used to understand and project global environmental change. Ecology and Earth system modeling must be integrated for scientists to fully comprehend the role of ecological systems in driving and responding to global change. Ecological insights can improve ESM realism and reduce process uncertainty, while ESMs offer ecologists an opportunity to broadly test ecological theory and increase the impact of their work by scaling concepts through time and space. Despite this mutualism, meaningfully integrating the two remains a persistent challenge, in part because of logistical obstacles in translating processes into mathematical formulas and identifying ways to integrate new theories and code into large, complex model structures. To help overcome this interdisciplinary challenge, we present a framework consisting of a series of interconnected stages for integrating a new ecological process or insight into an ESM. First, we highlight the multiple ways that ecological observations and modeling iteratively strengthen one another, dispelling the illusion that the ecologist's role ends with initial provision of data. Second, we show that many valuable insights, products, and theoretical developments are produced through sustained interdisciplinary collaborations between empiricists and modelers, regardless of eventual inclusion of a process in an ESM. Finally, we provide concrete actions and resources to facilitate learning and collaboration at every stage of data-model integration. This framework will create synergies that will transform our understanding of ecology within the Earth system, ultimately improving our understanding of global environmental change, and broadening the impact of ecological research.

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

Surface ozone (O3) 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 O3 measurements, observational-based assessments of O3 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 O3 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 O3 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 O3 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".

Our study suggests that the global CO2 fertilization effect (CFE) on vegetation photosynthesis has declined during the past four decades. The Comments suggest that the temporal inconsistency in AVHRR data and the attribution method undermine the results' robustness. Here, we provide additional evidence that these arguments did not affect our finding and that the global decline in CFE is robust.

Greening drylands despite warming consistent with carbon dioxide fertilization effect.

The rising atmospheric CO2 concentration leads to a CO2 fertilization effect on plants-that is, increased photosynthetic uptake of CO2 by leaves and enhanced water-use efficiency (WUE). Yet, the resulting net impact of CO2 fertilization on plant growth and soil moisture (SM) savings at large scale is poorly understood. Drylands provide a natural experimental setting to detect the CO2 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 CO2 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 CO2 concentration.

Interannual and Seasonal Drivers of Carbon Cycle Variability Represented by the Community Earth System Model (CESM2).

Earth system models are intended to make long-term projections, but they can be evaluated at interannual and seasonal time scales. Although the Community Earth System Model (CESM2) showed improvements in a number of terrestrial carbon cycle benchmarks, relative to its predecessor, our analysis suggests that the interannual variability (IAV) in net terrestrial carbon fluxes did not show similar improvements. The model simulated low IAV of net ecosystem production (NEP), resulting in a weaker than observed sensitivity of the carbon cycle to climate variability. Low IAV in net fluxes likely resulted from low variability in gross primary productivity (GPP)-especially in the tropics-and a high covariation between GPP and ecosystem respiration. Although lower than observed, the IAV of NEP had significant climate sensitivities, with positive NEP anomalies associated with warmer and drier conditions in high latitudes, and with wetter and cooler conditions in mid and low latitudes. We identified two dominant modes of seasonal variability in carbon cycle flux anomalies in our fully coupled CESM2 simulations that are characterized by seasonal amplification and redistribution of ecosystem fluxes. Seasonal amplification of net and gross carbon fluxes showed climate sensitivities mirroring those of annual fluxes. Seasonal redistribution of carbon fluxes is initiated by springtime temperature anomalies, but subsequently negative feedbacks in soil moisture during the summer and fall result in net annual carbon losses from land. These modes of variability are also seen in satellite proxies of GPP, suggesting that CESM2 appropriately represents regional sensitivities of photosynthesis to climate variability on seasonal time scales.

Solar geoengineering can alleviate climate change pressures on crop yields.

Solar geoengineering (SG) and CO2 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 CO2, 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 CO2 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.

Recent global decline of CO2 fertilization effects on vegetation photosynthesis.

The enhanced vegetation productivity driven by increased concentrations of carbon dioxide (CO2) [i.e., the CO2 fertilization effect (CFE)] sustains an important negative feedback on climate warming, but the temporal dynamics of CFE remain unclear. Using multiple long-term satellite- and ground-based datasets, we showed that global CFE has declined across most terrestrial regions of the globe from 1982 to 2015, correlating well with changing nutrient concentrations and availability of soil water. Current carbon cycle models also demonstrate a declining CFE trend, albeit one substantially weaker than that from the global observations. This declining trend in the forcing of terrestrial carbon sinks by increasing amounts of atmospheric CO2 implies a weakening negative feedback on the climatic system and increased societal dependence on future strategies to mitigate climate warming.

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 CO2 (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 CO2 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 CO2 fertilization. This challenges the view that CO2 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 CO2 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.

Towards a multiscale crop modelling framework for climate change adaptation assessment.

Predicting the consequences of manipulating genotype (G) and agronomic management (M) on agricultural ecosystem performances under future environmental (E) conditions remains a challenge. Crop modelling has the potential to enable society to assess the efficacy of G × M technologies to mitigate and adapt crop production systems to climate change. Despite recent achievements, dedicated research to develop and improve modelling capabilities from gene to global scales is needed to provide guidance on designing G × M adaptation strategies with full consideration of their impacts on both crop productivity and ecosystem sustainability under varying climatic conditions. Opportunities to advance the multiscale crop modelling framework include representing crop genetic traits, interfacing crop models with large-scale models, improving the representation of physiological responses to climate change and management practices, closing data gaps and harnessing multisource data to improve model predictability and enable identification of emergent relationships. A fundamental challenge in multiscale prediction is the balance between process details required to assess the intervention and predictability of the system at the scales feasible to measure the impact. An advanced multiscale crop modelling framework will enable a gene-to-farm design of resilient and sustainable crop production systems under a changing climate at regional-to-global scales.