Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Land Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

An Observation‐Driven Approach to Improve Vegetation Phenology in a Global Land Surface Model

An empirical model calibration approach is presented that aims to approximate missing biosphere processes in a global land surface model without the need for substantial model structural changes. The strategy is implemented here using the NASA Catchment‐CN land surface model and Moderate Resolution Imaging Spectroradiometer (MODIS) observations of the fraction of absorbed photosynthetically active radiation (FPAR). Existing plant functional types (PFTs) of the Catchment‐CN model are divided into three subtypes, based on the bias between the model‐simulated and MODIS‐observed FPAR. Separate sets of vegetation parameters for each subtype are then calibrated at a small number of grid cells with homogeneous, single‐PFT land cover, using MODIS FPAR reference observations from 2003 to 2009. The effectiveness of the empirical approach at improving the realism of modeled vegetation dynamics is investigated with two global model simulations for the period 2010–2016, one using the newly calibrated parameter values and the other using the original values. Globally, the calibrated parameters reduce the root mean square error (RMSE) of the modeled FPAR with respect to MODIS by 0.029 (∼10%) on average. In some regions, substantially larger RMSE reductions are achieved. RMSE reductions are primarily driven by model bias reductions, with neutral effects on the temporal correlation skill. While the empirical approach is suitable for achieving consistent model improvements, it is shown to be sensitive to the characteristics of the model error, specifically a dominance of the bias component in the case of Catchment‐CN. Ultimately, more fundamental model structural changes may be required to achieve better improvements.

J. Kolassa↗

ESS-DIVE guidelines for archiving terrestrial model data

This dataset contains supporting documents and images for ESS-DIVE terrestrial model data archiving guidelines.Terrestrial models are broadly defined as numerical models that couple both land dynamics and energy, water, carbon, or nutrient fluxes. We created these guidelines based on input from the U.S. Department of Energy’s Biological and Environmental Research land modeling community. The guidelines are intended to help modelers determine which components of their terrestrial model data associated with publication should be archived. Based on input from the land modeling community, the guidelines recommend archiving both model input and testing data, as well as code, script, and metadata. The guidelines also recommend archiving model data output, depending on the limitations set by data repositories. Lastly, we provide recommendations for bundling data files for publication as well as a discussion about tools that can facilitate model data archiving and reuse.This dataset is an archive of the associated GitHub repository for our model archiving guidelines (https://github.com/ess-dive-community/essdive-model-data-archiving-guidelines). The ‘README.pdf’ file gives a general introduction to the guidelines, and the ‘instructions.pdf’ file provides more detailed steps for following the guidelines. We also provide 2 figures in this data package: 1) a decision tree (model_data_guidelines_decision_tree.png) that can help users determine which components of their model data to archive. and 2) the ‘model_data_guidelines_flmd.png’ file depicts the different files that can be archived in addition to the model data itself. Lastly, we include 3 digitized tables from our associated manuscript and 3 CSV files with anonymized input from DOE scientists about the importance of different aspects of model data archiving from which we developed the guidelines.Dataset updates for v1.1.0: We updated this data package on 2021-11-22 in response to review comments on our related manuscript. In this update we removed one figure so that the model archiving guidelines are conveyed in text rather than an image. We updated the file-level metadata (FLMD) figure to be in accord with the most recent FLMD recommendations. We made minor edits to the README file to update the recommended citation and added two co-authors. We also added 6 new data files (3 are anonymized input from DOE scientists that helped to inform guidelines, and 3 are digitized tables from our manuscript.

54 ENVIRONMENTAL SCIENCES↗

Towards Global Simulation of Irrigation in a Land Surface Model: Multiple Cropping and Rice Paddy in Southeast Asia

Agricultural land use significantly influences the surface water and energy balances. Effects of irrigation on land surface states and fluxes include repartitioning of latent and sensible heat fluxes, an increase in net radiation, and an increase in soil moisture and runoff. We are working on representing irrigation practices in continental- to global-scale land surface simulation in NASA's Global Land Data Assimilation System (GLDAS). Because agricultural practices across the nations are diverse, and complex, we are attempting to capture the first-order reality of the regional practices before achieving a global implementation. This study focuses on two issues in Southeast Asia: multiple cropping and rice paddy irrigation systems. We first characterize agricultural practices in the region (i.e., crop types, growing seasons, and irrigation) using the Global data set of monthly irrigated and rainfed crop areas around the year 2000 (MIRCA2000) dataset. Rice paddy extent is identified using remote sensing products. Whether irrigated or rainfed, flooded fields need to be represented and treated explicitly. By incorporating these properties and processes into a physically based land surface model, we are able to quantify the impacts on the simulated states and fluxes.

Beaudoing, Hiroko Kato↗

A special issue on the RCPs

This paper provides an editors' introduction to the special issue of Climatic Change on the RCPs. Scenarios form a crucial element in climate change research. They allow researchers to explore the long-term consequences of decisions today, while taking account of the inertia in both the socio-economic and physical system. Scenarios also form an integrating element among the different research disciplines of those studying climate change, such as economists, technology experts, climate researchers, atmospheric chemists and geologists. In 2007, the IPCC requested the scientific community to develop a new set of scenarios, as the existing scenarios (published in the Special Report on Emissions Scenarios, (Nakicenovic and Swart 2000), and called the 'SRES scenarios') needed to be updated and expanded in scope (see Moss et al. (2010) for a detailed discussion). Researchers from different disciplines worked together to develop a process to craft these new scenarios, as summarized by Moss, et al. (2010). The Integrated Assessment Modeling Consortium (IAMC), founded in response to the IPCC call, played a key role in this process.1 The scenario development process aims to develop a set of new scenarios that facilitate integrated analysis of climate change across the main scientific communities. The process comprises 3 main phases: (1) an initial phase, developing a set of pathways for emissions, concentrations and radiative forcing, (2) a parallel phase, comprising both the development of new socio-economic storylines and climate model projections, and (3) an integration phase, combining the information from the first phases into holistic mitigation, impacts and vulnerability assessments. The pathways developed in the first phase were called 'Representative Concentration Pathways (RCPs)'. They play an important role in providing input for prospective climate model experiments, including both the decadal and long-term projections of climate change. The RCPs also provide an important reference point for new research within the integrated assessment modeling (IAM) community by standardizing on a common set of year-2100 conditions, and exploring alternative pathways and policies that could produce these outcomes. By design, the RCPs, as a set, cover the range of radiative forcing levels examined in the open literature and contain relevant information for climate model runs. This Special Issue documents the main assumptions and characteristics of the RCPs, and, in particular, the various steps that were involved in their development. A number of collaborative activities were initiated and finalized during the last 2-3 years to develop the RCPs. This required the cooperation of researchers from various disciplines involved in climate research, including emission experts, climate modelers, atmospheric chemistry modelers, land use modelers and experts involved in integrated assessment. The four RCPs together reflect the range of year-2100 radiative forcing values found in the literature, i.e. from 2.6 to 8.5 W/m{sup 2}. The papers in this Special Issue describe the individual RCPs, but also the various integrative steps that were necessary within the RCP development process to provide a harmonized set of pathways, that show a smooth transition from the past and extend far into the future for very long-term experiments. Important outcomes of this process included, for instance, the development of new emission inventories, new methods for the harmonization of spatial land use patterns, as well as extensions of the RCP trends beyond 2100. They briefly discuss the content of the individual papers.

54 ENVIRONMENTAL SCIENCES↗

Simulating groundwater uptake and hydraulic redistribution by phreatophytes in a high-resolution, coupled subsurface-land surface model

Several new functions representing groundwater dependent vegetation were incorporated into a coupled subsurface-land surface model, ParFlow.CLM, in order to adequately describe groundwater water uptake, hydraulic redistribution, and plant water stress. The modified model was used to conduct three-dimensional, stand-scale simulations of a Mediterranean oak savanna in California. It performed well and captured daily, hourly and spatial water and energy dynamics, as well as groundwater evapotranspiration rates. The new model was then compared to various approaches, the original ParFlow.CLM and a version using a root water uptake compensation equation. During the dry season, the modified model closely predicted the measured transpiration rate while the original model predicted that it would become zero and the compensation approach overestimated it by nearly double. The modified model also allowed for analysis of several key ecohydrological processes, namely the hydraulic redistribution when plants were both active and dormant, the leaf water potential, and xylem cavitation.

54 ENVIRONMENTAL SCIENCES↗

Tropical peatland hydrology simulated with a global land surface model

Tropical peatlands are among the most carbon-dense ecosystems on Earth, and their water storage dynamics strongly control these carbon stocks. The hydrological functioning of tropical peatlands differs from that of northern peatlands, which has not yet been accounted for in global land surface models (LSMs). Here, we integrated tropical peat-specific hydrology modules into a global LSM for the first time, by utilizing the peatland-specific model structure adaptation (PEATCLSM) of the NASA Catchment Land Surface Model (CLSM). We developed literature-based parameter sets for natural (PEATCLSM(Trop,Nat)) and drained (PEATCLSM(Trop,Drain)) tropical peatlands. Simulations with PEATCLSM(Trop,Nat) were compared against those with the default CLSM version and the northern version of PEATCLSM (PEATCLSM(North,Nat)) with tropical vegetation input. All simulations were forced with global meteorological reanalysis input data for the major tropical peatland regions in Central and South America, the Congo Basin, and Southeast Asia. The evaluation against a unique and extensive data set of in situ water level and eddy covariance-derived evapotranspiration showed an overall improvement in bias and correlation compared to the default CLSM version. Over Southeast Asia, an additional simulation with PEATCLSM(Trop,Drain) was run to address the large fraction of drained tropical peatlands in this region. PEATCLSM(Trop,Drain) outperformed CLSM, PEATCLSM(North,Nat) and PEATCLSM(Trop,Nat) over drained sites. Despite the overall improvements of PEATCLSM(Trop,Nat) over CLSM, there are strong differences in performance between the three study regions. We attribute these performance differences to regional differences in accuracy of meteorological forcing data, and differences in peatland hydrologic response that are not yet captured by our model.

Peatland↗

Integration of nitrogen dynamics into the Noah-MP land surface model v1.1 for climate and environmental predictions

Climate and terrestrial biosphere models consider nitrogen an important factor in limiting plant carbon uptake, while operational environmental models view nitrogen as the leading pollutant causing eutrophication in water bodies. The community Noah land surface model with multi-parameterization options (Noah-MP) is unique in that it is the next-generation land surface model for the Weather Research and Forecasting meteorological model and for the operational weather/climate models in the National Centers for Environmental Prediction. Here in this study, we add a capability to Noah-MP to simulate nitrogen dynamics by coupling the Fixation and Uptake of Nitrogen (FUN) plant model and the Soil and Water Assessment Tool (SWAT) soil nitrogen dynamics. This model development incorporates FUN's state-of-the-art concept of carbon cost theory and SWAT's strength in representing the impacts of agricultural management on the nitrogen cycle. Parameterizations for direct root and mycorrhizal-associated nitrogen uptake, leaf retranslocation, and symbiotic biological nitrogen fixation are employed from FUN, while parameterizations for nitrogen mineralization, nitrification, immobilization, volatilization, atmospheric deposition, and leaching are based on SWAT. The coupled model is then evaluated at the Kellogg Biological Station – a Long Term Ecological Research site within the US Corn Belt. Results show that the model performs well in capturing the major nitrogen state/flux variables (e.g., soil nitrate and nitrate leaching). Furthermore, the addition of nitrogen dynamics improves the modeling of net primary productivity and evapotranspiration. The model improvement is expected to advance the capability of Noah-MP to simultaneously predict weather and water quality in fully coupled Earth system models.

58 GEOSCIENCES↗

Utilization of Hydrologic Remote Sensing Data in Land Surface Modeling and Data Assimilation: Current Status and Challenges

Recent advances in remote sensing technologies have enabled the monitoring and measurement of the Earth's land surface at an unprecedented scale and frequency. The myriad of these land surface observations must be integrated with the state-of-the-art land surface model forecasts using data assimilation to generate spatially and temporally coherent estimates of environmental conditions. These analyses are of critical importance to real-world applications such as agricultural production, water resources management and flood, drought, weather and climate prediction. This need motivated the development of NASA Land Information System (LIS), which is an expert system encapsulating a suite of modeling, computational and data assimilation tools required to address challenging hydrological problems. LIS integrates the use of several community land surface models, use of ground and satellite based observations, data assimilation and uncertainty estimation techniques and high performance computing and data management tools to enable the assessment and prediction of hydrologic conditions at various spatial and temporal scales of interest. This presentation will focus on describing the results, challenges and lessons learned from the use of remote sensing data for improving land surface modeling, within LIS. More specifically, studies related to the improved estimation of soil moisture, snow and land surface temperature conditions through data assimilation will be discussed. The presentation will also address the characterization of uncertainty in the modeling process through Bayesian remote sensing and computational methods.

Kumar, Sujay V.↗

Land surface model parameter optimisation using in situ flux data: comparison of gradient-based versus random search algorithms (a case study using ORCHIDEE v1.9.5.2)

Land surface models (LSMs), which form the land component of earth system models, rely on numerous processes for describing carbon, water and energy budgets, often associated with highly uncertain parameters. Data assimilation (DA) is a useful approach for optimising the most critical parameters in order to improve model accuracy and refine future climate predictions. In this study, we compare two different DA methods for optimising the parameters of seven plant functional types (PFTs) of the ORCHIDEE LSM using daily averaged eddy-covariance observations of net ecosystem exchange and latent heat flux at 78 sites across the globe. We perform a technical investigation of two classes of minimisation methods – local gradient-based (the L-BFGS-B algorithm, limited memory Broyden–Fletcher–Goldfarb–Shanno algorithm with bound constraints) and global random search (the genetic algorithm) – by evaluating their relative performance in terms of the model–data fit and the difference in retrieved parameter values. We examine the performance of each method for two cases: when optimising parameters at each site independently (“single-site” approach) and when simultaneously optimising the model at all sites for a given PFT using a common set of parameters (“multi-site” approach). We find that for the single site case the random search algorithm results in lower values of the cost function (i.e. lower model–data root mean square differences) than the gradient-based method; the difference between the two methods is smaller for the multi-site optimisation due to a smoothing of the cost function shape with a greater number of observations. The spread of the cost function, when performing the same tests with 16 random first-guess parameters, is much larger with the gradient-based method, due to the higher likelihood of being trapped in local minima. When using pseudo-observation tests, the genetic algorithm results in a closer approximation of the true posterior parameter value in the L-BFGS-B algorithm. We demonstrate the advantages and challenges of different DA techniques and provide some advice on using it for the LSM parameter optimisation.

58 GEOSCIENCES↗

Land surface models systematically overestimate the intensity, duration and magnitude of seasonal-scale evaporative droughts

Land surface models (LSMs) must accurately simulate observed energy and water fluxes during droughts in order to provide reliable estimates of future water resources. We evaluated 8 different LSMs (14 model versions) for simulating evapotranspiration (ET) during periods of evaporative drought (Edrought) across six flux tower sites. Using an empirically defined Edrought threshold (a decline in ET below the observed 15th percentile), we show that LSMs simulated 58 Edrought days per year, on average, across the six sites, ~3 times as many as the observed 20 d. The simulated Edrought magnitude was ~8 times greater than observed and twice as intense. Our findings point to systematic biases across LSMs when simulating water and energy fluxes under water-stressed conditions. The overestimation of key Edrought characteristics undermines our confidence in the models' capability in simulating realistic drought responses to climate change and has wider implications for phenomena sensitive to soil moisture, including heat waves.

54 ENVIRONMENTAL SCIENCES↗

The NASA-Goddard Multi-Scale Modeling Framework - Land Information System: Global Land/atmosphere Interaction with Resolved Convection

The present generation of general circulation models (GCM) use parameterized cumulus schemes and run at hydrostatic grid resolutions. To improve the representation of cloud-scale moist processes and landeatmosphere interactions, a global, Multi-scale Modeling Framework (MMF) coupled to the Land Information System (LIS) has been developed at NASA-Goddard Space Flight Center. The MMFeLIS has three components, a finite-volume (fv) GCM (Goddard Earth Observing System Ver. 4, GEOS-4), a 2D cloud-resolving model (Goddard Cumulus Ensemble, GCE), and the LIS, representing the large-scale atmospheric circulation, cloud processes, and land surface processes, respectively. The non-hydrostatic GCE model replaces the single-column cumulus parameterization of fvGCM. The model grid is composed of an array of fvGCM gridcells each with a series of embedded GCE models. A horizontal coupling strategy, GCE4fvGCM4Coupler4LIS, offered significant computational efficiency, with the scalability and I/O capabilities of LIS permitting landeatmosphere interactions at cloud-scale. Global simulations of 2007e2008 and comparisons to observations and reanalysis products were conducted. Using two different versions of the same land surface model but the same initial conditions, divergence in regional, synoptic-scale surface pressure patterns emerged within two weeks. The sensitivity of largescale circulations to land surface model physics revealed significant functional value to using a scalable, multi-model land surface modeling system in global weather and climate prediction.

cumulus schemes↗

CO2 and CH4 leaf-level fluxes and soil porewater concentrations from common vegetation patches in Louisiana’s coastal wetlands

This dataset contains leaf-level flux and soil porewater concentration measurements of carbon dioxide (CO2) and methane (CH4 ) in plots in the footprint of Ameriflux sites US-LA2 and US-LA3. Leaf fluxes in US-LA2 were measured on patches dominated by Sagittaria lancifolia and co-dominated by Sagittaria lancifolia and Typha latifolia. In US-LA3, fluxes were measured from distinct Juncus roemerianus and Spartina alterniflora patches. The porewater concentrations were collected across a vertical profile (~50 cm depth) at centric locations within 25 m2 plots where we measured the leaf fluxes. US-LA3 included an additional set of measurements in open water spots. We aimed to evaluate differences in leaf fluxes and porewater pools of CO2 and CH4 of representative ecohydrological patches across a salinity gradient. We also used this dataset to help develop ELM-Wet, a more realistic representation of wetland carbon biogeochemical processes within the U.S. Department of Energy’s Energy Exascale Earth System Model (E3SM) Land Model version 1 (ELM v.1). The files can be opened with regular text editors or spreadsheet programs. Version 2.0 (8/26/2025): This is the latest version of this dataset. The update includes additional samples of soil porewater CH4/CO2 concentrations from June-2021 to November-2022, as well as minor adjustments made to V1 samples via changing Henry's solubility to account for porewater salinity. Additionally leaf-level measurments of spectral indices, PSRI, NDVI, and PRI have been added to complement Leaf-level flux measurements. All V1 data sets have been integrated into V2 sheets, ensuring data from the previous version is contained with the additional samples and consistent with V2 metadata.

54 ENVIRONMENTAL SCIENCES↗

Model Structure and Climate Data Uncertainty in Historical Simulations of the Terrestrial Carbon Cycle (1850–2014)

Abstract The divergence among Earth system models in the terrestrial carbon cycle has prompted interest in how to reduce uncertainty. Previous studies have identified model structural uncertainty arising from process parameterizations and parameter values. The current study highlights the importance of climate forcing in generating carbon cycle uncertainty. We use simulations in which three models (Community Land Model version 4 (CLM4), CLM4.5, CLM5) with substantially different carbon cycles are forced with two climate reconstructions (CRUNCEPv7, Global Soil Wetness Project 3 version 1 (GSWP3v1)) to examine the contributions of model structure and climate to uncertainty in the carbon cycle over the period 1850–2014. Climate uncertainty for global annual net biome production exceeds one third of total uncertainty (defined as the sum of climate and model structure uncertainty) in the first half of the twentieth century, but declines after the 1950s. Global annual gross primary productivity, net primary productivity, heterotrophic respiration, and vegetation and soil carbon stocks have substantial climate uncertainty (relative to total uncertainty) throughout the simulation period. Climate forcing contributes more than one half of total uncertainty for these carbon cycle fluxes and stocks throughout boreal North America and Eurasia, some midlatitude regions, and in eastern Amazonia and western equatorial Africa during the decade 2000–2009. Comparison with observationally based data sets of the carbon cycle using model benchmarking methods provides insight into strengths and deficiencies among models and climate forcings, but we caution against overreliance on benchmarking to discriminate among models. The conceptualization of uncertainty arising from this study implies embracing multiple feasible model simulations rather than focusing on which model or simulation is best.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Joint Effects of Vegetation Characteristics and Soil Properties on Ecosystem Dynamics in a Panama Tropical Forest

Abstract In tropical forests, both vegetation characteristics and soil properties are important not only for controlling energy, water, and gas exchanges directly but also determining the competition among species, successional dynamics, forest structure and composition. However, the joint effects of the two factors have received limited attention in Earth system model development. Here we use a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM‐FATES, to explore how plant traits and soil properties affect tropical forest growth and composition concurrently. A large ensemble of simulations with perturbed vegetation and soil hydrological parameters is conducted at the Barro Colorado Island, Panama. The simulations are compared against observed carbon, energy, and water fluxes. We find that soil hydrological parameters, particularly the scaling exponent of the soil retention curve ( B sw ), play crucial roles in controlling forest diversity, with higher B sw values (>7) favoring late successional species in competition, and lower B sw values (1 ∼ 7) promoting the coexistence of early and late successional plants. Considering the additional impact of soil properties resolves a systematic bias of FATES in simulating sensible/latent heat partitioning with repercussion on water budget and plant coexistence. A greater fraction of deeper tree roots can help maintain the dry‐season soil moisture and plant gas exchange. As soil properties are as important as vegetation parameters in predicting tropical forest dynamics, more efforts are needed to improve parameterizations of soil functions and belowground processes and their interactions with aboveground vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗

NGEE Arctic Field-to-Model

This repository contains workshop documentation and shell scripts/computational infrastructure for building, running, and analyzing output from the DOE Energy Exascale Earth System Model (E3SM) Land Model (ELM).

Fiorella, Rich [@lanl]↗

ELM model simulations of Plum Island Ecosystems LTER low marsh site 2018-2020

Model simulations using the Department of Energy's Energy Exascale Earth System Model (E3SM) land model (ELM) with improved capabilities to represent vegetation response to salinity and inundation. The simulations were conducted for a tidal salt marsh at Plum Island Ecosystems Long Term Ecological Research (LTER) site near Rowley, Massachusetts, USA; the site is a low marsh dominated by Spartina alterniflora. The model was forced with site-specific meteorology, salinity and tidal cycles from 2018-2020. Four sets of model simulations are included and described below:1. Parameterization of the salinity response function. These simulations tested different combinations of values for optimal salinity and salinity tolerance.2. Model evaluation. This comparison conducted simulations using the default model, the salinity function only, the submergence function only, and both the salinity and submergence functions. 3. Salinity scenarios. These simulations used the 2018 salinity input data varied by -5 to +10 ppt salinity.4. Water level scenarios. These simulations used the tide height varied by -10 to +50 cm. These simulations were used to demonstrate how the salinity and submergence functions better represent carbon uptake by tidal salt marshes.The data package includes netCDF files used as forcing files for tide height and salinity, one for each year 2018-2020 at observed salinity concentrations, and an additional three forcing files in which salinity concentrations were varied 5 ppt lower, 5 ppt higher, and 10 ppt higher than the measured 2018 time series. Also included are python scripts for creating forcing files, plain text parameter and command files for running simulations, model outputs in netCDF format, and python scripts for visualizing outputs. Code for the modified E3SM model is archived in Sulman et al 2023 at doi:10.15485/1991625. More detail about files is provided in the README.md file.

54 ENVIRONMENTAL SCIENCES↗

Data-Model Needs for Belowground Ecology. A Summary Report from the Terrestrial Ecosystem Science (TES) Mini-Workshop, May 8, 2014

Terrestrial life on Earth depends on belowground cycling and storage of a significant portion of the planet’s biologically-active carbon. The interplay among organic inputs by primary producers, soil organic matter stabilization, and assimilation and mineralization by soil organisms drives belowground ecology. Understanding belowground carbon dynamics, including the stability and vulnerability of soil organic carbon stocks, is necessary to project the responses and feedbacks of terrestrial ecosystems to a changing climate. Although conceptually understood, belowground hydrobiogeochemical systems are not well characterized with quantitative and mechanistic representations of critical functions. As a result, these systems are inadequately represented in modeling efforts, limiting our ability to understand and project terrestrial carbon source/sink relationships at various temporal and spatial scales. The belowground science community lacks robust databases and data libraries against which model sensitivities and performance can be tested and enhanced. Improved understanding of belowground ecosystems, including biotic and abiotic interactions, is needed to ensure that model structures include and provide adequate representations of critical belowground processes, functions, and dynamics. A one-day workshop brought together experimental and modeling representatives from the four Department of Energy (DOE) National Laboratories where the Office of Biological and Environmental Research (BER) Terrestrial Ecosystem Science (TES) program has investments in belowground ecology. The group focused on identifying and characterizing (1) critical experimental and observational datasets for testing and developing existing and future belowground-process models and land models; (2) the limitations and needs for modeling belowground processes within ecosystem models; and (3) key areas for leveraging current research activities across the National Laboratories and BER-funded projects. Discussion results are organized in this report as: (I) over-arching data requirements for advancing understanding of belowground ecosystems; (II) efforts that could rapidly improve understanding and model representations of belowground components of terrestrial ecosystems; and (III) longer-term research needs that would benefit from greater and more sustained support.

54 ENVIRONMENTAL SCIENCES↗

Effects of Rapid Permafrost Thaw on CO2 and CH4 Fluxes in a Warmer and Wetter Future (Final Technical Report)

When ice-rich permafrost thaws, the ground subsides, creating thermokarst landscapes with dramatically different soil conditions and carbon fluxes than the original ecosystem. Roughly 20% of the northern permafrost region is susceptible to thermokarst formation (Olefeldt et al., 2016). Thermokarst formation is often rapid; tens of meters of permafrost can thaw within a few years (Schuur et al., 2015). In topographically low areas, thermokarst thaw converts boreal forest or tundra dry shrub ecosystems into sedge or Sphagnum moss wetlands (Olefeldt et al., 2016). While this type of landscape transformation releases carbon stored in permafrost into the atmosphere, on longer time scales, it facilitates sequestration of atmospheric carbon in plant biomass because permafrost thaw releases plant-available nutrients and wetlands are highly productive (M. C. Jones et al., 2017). However, wetlands also generate methane, which is a potent greenhouse gas. Methane emissions from thermokarst wetlands can cause these carbon-sequestering systems to have a positive global warming potential (Johansson et al., 2006; Turetsky et al., 2007). Our project objective was to improve Earth System and environmental predictability by advancing understanding of how CO 2 and CH 4 flux in permafrost thaw-induced wetlands (thermokarst) will change in the future as temperatures and climate conditions shift. Northern latitudes are expected to get warmer and wetter (IPCC 2013), and initiation and expansion of thermokarst thaw is expected to increase (Jorgenson et al. 2006; Zhang et al. 2017). Given these expected changes, our work sought to address three broad questions: Q1) How will northern latitude CO2 and CH4 emissions respond to warming temperatures? Q2) What is the impact of precipitation on permafrost thaw and carbon emissions? Q3) How do CO2 and CH4 emissions change as wetlands age after permafrost thaw? To answer these questions, we took both a modeling and measurement approach. Modeling work was conducted with DOE’s Energy Exascale Earth System Model (E3SM) land model (ELM). Empirical work took place in two primary locations. The first was a thermokarst site near Fairbanks, AK. The site is part of the Bonanza Creek Long Term Ecological Research program and is well instrumented (Neumann et al., 2019). The second was an isolated thawing permafrost wetland located on Kenai Peninsula — Brown’s Lake bog (B. M. Jones et al., 2016) — where the current climate is representative of what is expected at higher latitudes in the future. This site provides an ideal opportunity to test our hypotheses about how thermokarst wetland dynamics will respond to future environmental conditions. In addition, the project collaborated with researchers asking similar questions who were collecting measurements in high latitude post-glacial lakes located in collaborators tackling similar questions in Sweden (Stordalen Mire) (Emerson et al., 2021). Project efforts directly aligned with the stated goal of the funding opportunity announcement (FOA), which was “to improve the understanding and representation of terrestrial ecosystems in ways that advance Earth system model parameterizations and capabilities... thereby improving the quality of Earth and environmental model projections and providing the scientific foundation needed to support DOE’s science and energy missions.” Specifically, the project improved sophistication and accuracy of the Energy Exascale Earth System Model (E3SM), which is being developed primarily at DOE National Laboratories to support scientific research and decision-making.

54 ENVIRONMENTAL SCIENCES↗