Beyond ecosystem modeling: A roadmap to community cyberinfrastructure for ecological data-model integration

Toward a coordinated understanding of hydro-biogeochemical root functions in tropical forests for application in vegetation models

Tropical forest root characteristics and resource acquisition strategies are underrepresented in vegetation and global models, hampering the prediction of forest-climate feedbacks for these carbon-rich ecosystems. Lowland tropical forests often have globally unique combinations of high taxonomic and functional biodiversity, rainfall seasonality, and strongly weathered infertile soils, giving rise to distinct patterns in root traits and functions compared with higher latitude ecosystems. We provide a roadmap for integrating recent advances in our understanding of tropical forest belowground function into vegetation models, focusing on water and nutrient acquisition. We offer comparisons of recent advances in empirical and model understanding of root characteristics that represent important functional processes in tropical forests. We focus on: (1) fine-root strategies for soil resource exploration, (2) coupling and trade-offs in fine-root water vs nutrient acquisition, and (3) aboveground-belowground linkages in plant resource acquisition and use. We suggest avenues for representing these extremely diverse plant communities in computationally manageable and ecologically meaningful groups in models for linked aboveground-belowground hydro-nutrient functions. Tropical forests are undergoing warming, shifting rainfall regimes, and exacerbation of soil nutrient scarcity caused by elevated atmospheric CO2. The accurate model representation of tropical forest functions is crucial for understanding the interactions of this biome with the climate.

Net greenhouse gas balance in U.S. croplands: How can soils be part of the climate solution?

Agricultural soils play a dual role in regulating the Earth's climate by releasing or sequestering carbon dioxide (CO2 ) in soil organic carbon (SOC) and emitting non-CO2 greenhouse gases (GHGs) such as nitrous oxide (N2 O) and methane (CH4 ). To understand how agricultural soils can play a role in climate solutions requires a comprehensive assessment of net soil GHG balance (i.e., sum of SOC-sequestered CO2 and non-CO2 GHG emissions) and the underlying controls. Herein, we used a model-data integration approach to understand and quantify how natural and anthropogenic factors have affected the magnitude and spatiotemporal variations of the net soil GHG balance in U.S. croplands during 1960-2018. Specifically, we used the dynamic land ecosystem model for regional simulations and used field observations of SOC sequestration rates and N2 O and CH4 emissions to calibrate, validate, and corroborate model simulations. Results show that U.S. agricultural soils sequestered13.2 ± 1.16 $$ 13.2\pm 1.16 $$ Tg CO2 -C year-1 in SOC (at a depth of 3.5 m) during 1960-2018 and emitted0.39 ± 0.02 $$ 0.39\pm 0.02 $$ Tg N2 O-N year-1 and0.21 ± 0.01 $$ 0.21\pm 0.01 $$ Tg CH4 -C year-1 , respectively. Based on the GWP100 metric (global warming potential on a 100-year time horizon), the estimated national net GHG emission rate from agricultural soils was122.3 ± 11.46 $$ 122.3\pm 11.46 $$ Tg CO2 -eq year-1 , with the largest contribution from N2 O emissions. The sequestered SOC offset ~28% of the climate-warming effects resulting from non-CO2 GHG emissions, and this offsetting effect increased over time. Increased nitrogen fertilizer use was the dominant factor contributing to the increase in net GHG emissions during 1960-2018, explaining ~47% of total changes. In contrast, reduced cropland area, the adoption of agricultural conservation practices (e.g., reduced tillage), and rising atmospheric CO2 levels attenuated net GHG emissions from U.S. croplands. Improving management practices to mitigate N2 O emissions represents the biggest opportunity for achieving net-zero emissions in U.S. croplands. Our study highlights the importance of concurrently quantifying SOC-sequestered CO2 and non-CO2 GHG emissions for developing effective agricultural climate change mitigation measures.

Digital twins: dynamic model-data fusion for ecology

Digital twins (DTs) are an emerging phenomenon in the public and private sectors as a new tool to monitor and understand systems and processes. DTs have the potential to change the status quo in ecology as part of its digital transformation. However, it is important to avoid misguided developments by managing expectations about DTs. We stress that DTs are not just big models of everything, containing big data and machine learning. Rather, the strength of DTs is in combining data, models, and domain knowledge, and their continuous alignment with the real world. We suggest that researchers and stakeholders exercise caution in DT development, keeping in mind that many of the strengths and challenges of computational modelling in ecology also apply to DTs.

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.

Progress and opportunities in advancing near-term forecasting of freshwater quality

Near-term freshwater forecasts, defined as sub-daily to decadal future predictions of a freshwater variable with quantified uncertainty, are urgently needed to improve water quality management as freshwater ecosystems exhibit greater variability due to global change. Shifting baselines in freshwater ecosystems due to land use and climate change prevent managers from relying on historical averages for predicting future conditions, necessitating near-term forecasts to mitigate freshwater risks to human health and safety (e.g., flash floods, harmful algal blooms) and ecosystem services (e.g., water-related recreation and tourism). To assess the current state of freshwater forecasting and identify opportunities for future progress, we synthesized freshwater forecasting papers published in the past 5 years. We found that freshwater forecasting is currently dominated by near-term forecasts of water quantity and that near-term water quality forecasts are fewer in number and in the early stages of development (i.e., non-operational) despite their potential as important preemptive decision support tools. We contend that more freshwater quality forecasts are critically needed and that near-term water quality forecasting is poised to make substantial advances based on examples of recent progress in forecasting methodology, workflows, and end-user engagement. For example, current water quality forecasting systems can predict water temperature, dissolved oxygen, and algal bloom/toxin events 5 days ahead with reasonable accuracy. Continued progress in freshwater quality forecasting will be greatly accelerated by adapting tools and approaches from freshwater quantity forecasting (e.g., machine learning modeling methods). In addition, future development of effective operational freshwater quality forecasts will require substantive engagement of end users throughout the forecast process, funding, and training opportunities. Looking ahead, near-term forecasting provides a hopeful future for freshwater management in the face of increased variability and risk due to global change, and we encourage the freshwater scientific community to incorporate forecasting approaches in water quality research and management.

Development of a data-assimilation system to forecast agricultural systems: A case study of constraining soil water and soil nitrogen dynamics in the APSIM model

As we face today's large-scale agricultural issues, the need for robust methods of agricultural forecasting has never been clearer. Yet, the accuracy and precision of our forecasts remains limited by current tools and methods. To overcome the limitations of process-based models and observed data, we iteratively designed and tested a generalizable and robust data-assimilation system that systematically constrains state variables in the APSIM model to improve forecast accuracy and precision. Our final novel system utilizes the Ensemble Kalman Filter to constrain model states and update model parameters at observed time steps and incorporates an algorithm that improves system performance through the joint estimation of system error matrices. We tested this system at the Energy Farm, a well-monitored research site in central Illinois, where we assimilated observed in situ soil moisture at daily time steps for two years and evaluated how assimilation impacted model forecasts of soil moisture, yield, leaf area index, tile flow, and nitrate leaching by comparing estimates with in situ observations. The system improved the accuracy and precision of soil moisture estimates for the assimilation layers by an average of 42% and 48%, respectively, when compared to the free model. Such improvements led to changes in the model's soil water and nitrogen processes and, on average, increased accuracy in forecasts of annual tile flow by 43% and annual nitrate loads by 10%. Forecasts of aboveground measures did not dramatically change with assimilation, a fact which highlights the limited potential of soil moisture as a constraint for a site with no water stress. Extending the scope of previous work, our results demonstrate the power of data assimilation to constrain important model estimates beyond the assimilated state variable, such as nitrate leaching. Replication of this study is necessary to further define the limitations and opportunities of the developed system.

PERFICT: A Re-imagined foundation for predictive ecology

Making predictions from ecological models—and comparing them to data—offers a coherent approach to evaluate model quality, regardless of model complexity or modelling paradigm. To date, our ability to use predictions for developing, validating, updating, integrating and applying models across scientific disciplines while influencing management decisions, policies, and the public has been hampered by disparate perspectives on prediction and inadequately integrated approaches. We present an updated foundation for Predictive Ecology based on seven principles applied to ecological modelling: make frequent Predictions, Evaluate models, make models Reusable, Freely accessible and Interoperable, built within Continuous workflows that are routinely Tested (PERFICT). We outline some benefits of working with these principles: accelerating science; linking with data science; and improving science‐policy integration.

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.

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.

Bridging implementation gaps to connect large ecological datasets and complex models

Merging robust statistical methods with complex simulation models is a frontier for improving ecological inference and forecasting. However, bringing these tools together is not always straightforward. Matching data with model output, determining starting conditions, and addressing high dimensionality are some of the complexities that arise when attempting to incorporate ecological field data with mechanistic models directly using sophisticated statistical methods. To illustrate these complexities and pragmatic paths forward, we present an analysis using tree-ring basal area reconstructions in Denali National Park (DNPP) to constrain successional trajectories of two spruce species (Picea mariana and Picea glauca) simulated by a forest gap model, University of Virginia Forest Model Enhanced-UVAFME. Through this process, we provide preliminary ecological inference about the long-term competitive dynamics between slow-growing P. mariana and relatively faster-growing P. glauca. Incorporating tree-ring data into UVAFME allowed us to estimate a bias correction for stand age with improved parameter estimates. We found that higher parameter values for P. mariana minimum growth under stress and P. glauca maximum growth rate were key to improving simulations of coexistence, agreeing with recent research that faster-growing P. glauca may outcompete P. mariana under climate change scenarios. The implementation challenges we highlight are a crucial part of the conversation for how to bring models together with data to improve ecological inference and forecasting.

Coding for Life: Designing a Platform for Projecting and Protecting Global Biodiversity

Time is running out to limit further devastating losses of biodiversity and nature's contributions to humans. Addressing this crisis requires accurate predictions about which species and ecosystems are most at risk to ensure efficient use of limited conservation and management resources. We review existing biodiversity projection models and discover problematic gaps. Current models usually cannot easily be reconfigured for other species or systems, omit key biological processes, and cannot accommodate feedbacks with Earth system dynamics. To fill these gaps, we envision an adaptable, accessible, and universal biodiversity modeling platform that can project essential biodiversity variables, explore the implications of divergent socioeconomic scenarios, and compare conservation and management strategies. We design a roadmap for implementing this vision and demonstrate that building this biodiversity forecasting platform is possible and practical.

Global patterns of forest autotrophic carbon fluxes

Carbon (C) fixation, allocation, and metabolism by trees set the basis for energy and material flows in forest ecosystems and define their interactions with Earth's changing climate. However, while many studies have considered variation in productivity with latitude and climate, we lack a cohesive synthesis on how forest carbon fluxes vary globally with respect to climate and one another. Here, we draw upon 1,319 records from the Global Forest Carbon Database, representing all major forest types and the nine most significant autotrophic carbon fluxes, to comprehensively review how annual C cycling in mature, undisturbed forests varies with latitude and climate on a global scale. Across all flux variables analyzed, rates of C cycling decreased continuously with absolute latitude-a finding that confirms multiple previous studies and contradicts the idea that net primary productivity of temperate forests rivals that of tropical forests. C flux variables generally displayed similar trends across latitude and multiple climate variables, with no differences in allocation detected at this global scale. Temperature variables in general, and mean annual temperature or temperature seasonality in particular, were the best single predictors of C flux, explaining 19%-71% of variation in the C fluxes analyzed. The effects of temperature were modified by moisture availability, with C flux reduced under hot and dry conditions and sometimes under very high precipitation. Annual C fluxes increased with growing season length and were also influenced by growing season climate. These findings clarify how forest C flux varies with latitude and climate on a global scale. In an era when forests will play a critical yet uncertain role in shaping Earth's rapidly changing climate, our synthesis provides a foundation for understanding global patterns in forest C cycling.