Search NASA⌕ Search

SEARCH · Search NASA

Results for “Land subsidence”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rising concerns of climate extremes and land subsidence impacts

A recent article in Reviews of Geophysics explores land subsidence drivers, rates, and impacts across the globe. It also discusses the need for improved process representations and the inclusion of the interplay among land subsidence and climatic extremes, including their effects in models and risk assessments. Here, we asked the lead author to explain the concept of land subsidence, its impacts, and future directions needed for improved mitigation.

Earth science↗
DOE OSTI2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics↗
DOE OSTI2023

Predicting Short-Term Deformation in the Central Valley Using Machine Learning

Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers’ integrity.

58 GEOSCIENCES↗
DOE OSTI2025

Sediment transport and morphodynamic changes due to upstream variation and sea level rise on the Western Mekong Delta, Vietnam

The hydrological and morphological dynamics across the Vietnam Mekong Delta (VMD) coastline are increasingly shaped by human interventions and climate change, affecting the sustainability of socio-economic environments. This study employs the TELEMAC-MASCARET modeling system (TELEMAC-2D, TOMAWAC, and SISYPHE modules) to simulate hydrodynamics, sediment transport, and morphological changes for a 10-year period under baseline (SC1) and future climate (SC2) scenarios, considering wave climatology of Northeast (NEM) and Southwest monsoons (SWM). The results of sedimentation and erosion were analyzed across five key coastal zones from Ca Mau Cape (KV1) to Ha Tien (KV5). Model results highlight distinct seasonal patterns of sediment redistribution shaped by monsoonal flows, wave action, and seabed morphology. During the NEM, sediment deposition is concentrated near estuarine zones, while the SWM favors offshore sediment accumulation due to perpendicular wave actions, and shallow, gently sloping bathymetry that dissipates wave energy. Coastal erosion predominantly occurs at convex shoreline features subjected to stronger coastal currents, particularly in fine-grained sediment areas. Under future sea level rise and land subsidence scenarios, reduced flow velocities lead to lower bottom shear stress, promoting deposition and mitigating erosion problems despite increased wave energy. Besides that, sediment transport analyses across 14 cross-sections (C-S) reveal a southward trend during the NEM and a northward reversal during the SWM, with upstream sediment reductions primarily influencing near-estuarine zones and exerting minimal impact beyond the mid-coastal regions. Furthermore, these findings advance the understanding of sediment redistribution processes under climate change and seasonal change, and provide a scientific basis for adaptive coastal management and sustainable development.

Climate change↗
DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗
DOE OSTI2023

A Fusion of Geothermal and InSAR Data with Machine Learning for Enhanced Deformation Forecasting at the Geysers

The Geysers geothermal field in California is experiencing land subsidence due to the seismic and geothermal activities taking place. This poses a risk not only to the underlying infrastructure but also to the groundwater level which would reduce the water availability for the local community. Because of this, it is crucial to monitor and assess the surface deformation occurring and adjust geothermal operations accordingly. In this study, we examine the correlation between the geothermal injection and production rates as well as the seismic activity in the area, and we show the high correlation between the injection rate and the number of earthquakes. This motivates the use of this data in a machine learning model that would predict future deformation maps. First, we build a model that uses interferometric synthetic aperture radar (InSAR) images that have been processed and turned into a deformation time series using LiCSBAS, an open-source InSAR time series package, and evaluate the performance against a linear baseline model. The model includes both convolutional neural network (CNN) layers as well as long short-term memory (LSTM) layers and is able to improve upon the baseline model based on a mean squared error metric. Then, after getting preprocessed, we incorporate the geothermal data by adding them as additional inputs to the model. This new model was able to outperform both the baseline and the previous version of the model that uses only InSAR data, motivating the use of machine learning models as well as geothermal data in assessing and predicting future deformation at The Geysers as part of hazard mitigation models which would then be used as fundamental tools for informed decision making when it comes to adjusting geothermal operations.

Yazbeck, Joe (ORCID:0000000302235260)↗
DOE OSTI2023

High–Resolution Maps of Near–Surface Permafrost for Three Watersheds on the Seward Peninsula, Alaska Derived From Machine Learning

Permafrost soils are a critical component of the global carbon cycle and are locally important because they regulate the hydrologic flux from uplands to rivers. Furthermore, degradation of permafrost soils causes land surface subsidence, damaging infrastructure that is crucial for local communities. Regional and hemispherical maps of permafrost are too coarse to resolve distributions at a scale relevant to assessments of infrastructure stability or to illuminate geomorphic impacts of permafrost thaw. Here we train machine learning models to generate meter–scale maps of near–surface permafrost for three watersheds in the discontinuous permafrost region. The models were trained using ground truth determinations of near–surface permafrost presence from measurements of soil temperature and electrical resistivity. We trained three classifiers: extremely randomized trees (ERTr), support vector machines (SVM), and an artificial neural network (ANN). Model uncertainty was determined using k–fold cross validation, and the modeled extents of near–surface permafrost were compared to the observed extents at each site. At–a–site near–surface permafrost distributions predicted by the ERTr produced the highest accuracy (70%–90%). However, the transferability of the ERTr to the sites outside of the training data set was poor, with accuracies ranging from 50% to 77%. The SVM and ANN models had lower accuracies for at–a–site prediction (70%–83%), yet they had greater accuracy when transferred to the non–training site (62%–78%). These models demonstrate the potential for integrating high–resolution spatial data and machine learning models to develop maps of near–surface permafrost extent at resolutions fine enough to assess infrastructure vulnerability and landscape morphology influenced by permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impact of Global Warming on U.S. Summertime Mesoscale Convective Systems: A Simple Lagrangian Parcel Model Perspective

Abstract Mesoscale convective systems (MCSs) are the dominant rainfall producer over the United States during the warm season, causing natural disasters and severe weather every year. Global climate models have large uncertainty in projecting precipitation changes in the future climate. Here, a simple Lagrangian parcel model is used to investigate the impact of global warming on MCS initiation and growth. The single-column parcel model projects a mean precipitation decrease over the central United States and an increase to its east, in agreement with the CMIP5 model projection. It also highlights the crucial role of current climate mean-state model bias in exaggerating the change in future mean precipitation projection by 25%. As for convective population, the model captures the decreased occurrence frequency of weak to moderate convection and increased frequency of strong convection due to the increased CAPE and CIN, in agreement with convection-permitting regional simulations. Novel parameterizations of gust-front propagation speed and subsidence strength are developed as guided by cloud-resolving simulations. The multicolumn parcel model employing those parameterizations captures readily the cold pool–induced upscale growth feature. It simulates smaller mesoscale clusters over the central United States under global warming due to gust-front slowdown and subsidence strength enhancement, which are further attributed to land aridity–induced weakening of initial accumulated precipitation and strengthening of updraft speed, respectively. That said, mesoscale clusters could become bigger under more favorable conditions in future climate, including boundary layer moistening, convection lifetime lengthening, and cold pool mechanical-lifting enhancement, which require further investigations to improve mechanistic understanding of future MCS changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Response of the Current Climate to Land‐Ocean Contrasts in Parameterized Cumulus Entrainment

Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temporal Patterns in Soil Redox Potential Vary Across a Freshwater Coastal Delta

Widespread and persistent flooding is submerging coastal ecosystems, particularly along the Louisiana Gulf coast, where flooding results from the compounding effects of ground subsidence and rising sea level. Restoration projects aim to mitigate land loss by diverting sediment loads from rivers into degraded areas to increase ground elevation. To predict how coastal ecosystems will change over time in response to projected changes in relative sea level and restoration, it is necessary to understand how subsurface biogeochemical processes respond to dynamic hydrologic forcings. Here, this study evaluates how environmental parameters that integrate biogeochemical processes vary with water table fluctuations in the freshwater Wax Lake Delta (WLD) in Louisiana, USA, where water diversions have formed one of the only active deltas along the coast. High‐frequency observations of water level, soil redox potential, specific conductance and pH were made for 1 year along elevation transects located on the older, proximal and younger, distal ends of a deltaic island. Redox responded rapidly to changing water tables, with fluctuations occurring primarily in shallow soils (< 20 cm) and at higher elevations. Deeper soils and those at lower elevation remained inundated and reduced. Semi‐diurnal tidal fluctuations were pronounced in younger, distal soils, presumably due to rapid groundwater exchange with the river channel. Tidal signals were muted in older soils that instead exhibited seasonal variability associated with river discharge and evapotranspiration. Although much of the delta sediments are persistently reducing and anoxic, redox fluctuations in the natural levees that border the deltaic islands likely drive high rates of biogeochemical activity. Evaluating how hydrology drives the frequency and duration of redox fluctuations provides a basis for understanding how biogeochemical processes might vary with complex hydrological interactions in coastal systems.

Herndon, Elizabeth [Oak Ridge National Laboratory ↗
DOE OSTI2025

AmeriFlux US-YNS North Star Yedoma

This is the AmeriFlux version of the carbon flux data for the site US-YNS North Star Yedoma. Site Description - Dry hillslope underlain by ice rich permafrost and now covered with native grasses and managed as a rustic golf course. The land clearing performed decades ago has caused permafrost thawing, surface subsidence, and development of thermokarst mounds or baidjarahks. No foreign fill materials, pesticides or fertilizers have been used, and the area is not irrigated. The area outside the flux footprint is second growth deciduous forest with similar soil conditions.

Anthony, Katey Walter [University of Alaska Fairba↗
DOE OSTI2024

Characterizing Land-Atmosphere Interactions During the Afternoon-to-Evening Transition Using ARM SGP Observations

The land-atmosphere (LA) system consists of the soil, vegetation, and lowest layer of the atmosphere. The energy and water transport processes in this system contribute to the structure and the daily cycle of the atmospheric boundary layer (ABL), which is the layer between the land surface and the free troposphere that controls the exchange of momentum, energy, and water. The state of the ABL is the result of complex two-way interactions between local processes such as surface fluxes and entrainment at the ABL top as well as large-scale and mesoscale processes such as advection and subsidence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Changing Interactions Between Trace Gas Fluxes, Belowground Chemistry, and Plant Traits Across an Arctic Thermokarst Landscape

Arctic permafrost soils are increasingly subject to thermokarst that is, abrupt ground subsidence caused by thaw. Wetlands can form within these depressions, leading to changes in organic matter decomposition and gas fluxes (CO 2 , CH 4 , N 2 O, NH 3 ). Thermokarst wetlands tend to be dominated by graminoids, while surrounding upland tussock tundra tends to be dominated by mixed communities of shrubs and graminoids. Here, to investigate how thermokarst alters the land-atmosphere exchange of C and N gases in Arctic tundra, we analyzed soil, porewater, above- and belowground biomass, and measured gas fluxes across dominant plant functional types (PFTs) within a lowland thermokarst wetland and adjacent upland tussock tundra. Both locations were overall sinks of CO 2 , sources of CH 4 , and sources of both N 2 O and NH 3 . We found that thermokarst wetlands emitted enough CH 4 to generate a positive radiative forcing in CO 2 equivalents (+1.2 μmol m −2 s −1 CO 2 -eq), counteracting the high CO 2 uptake. In contrast, the upland tussock tundra had a net negative radiative forcing (−1.2 μmol m −2 s −1 CO 2 -eq). Differences in gas flux and soil chemistry between upland and lowland are primarily driven by flooded conditions present in thermokarst wetland. Additionally, root biomass from graminoids across both lowlands and uplands significantly correlated with CH 4 fluxes, supporting previous observations of plant-mediated transport of CH 4 . Graminoid cover was correlated with increases in low molecular weight dissolved organic carbon, possibly associated with root exudates that fuel methanogenesis. Forb cover in the upland tussock tundra was significantly correlated with nine soil chemical variables, indicating that forbs may influence local soil chemistry or conversely, that soil chemistry controls where forbs grow. Overall, our findings indicate the variability in gas fluxes in the upland tussock tundra is partially controlled by PFT cover, while thermokarst wetlands emit enough CH 4 to counteract CO 2 uptake, with implications for carbon budget changes in Arctic systems.

Land-atmosphere feedback↗
DOE OSTI2025

Large Eddy Simulation of the Diurnal Cycle of Shallow Convection in the Central Amazon

Climate models often face challenges in accurately simulating the daily precipitation cycle over tropical land areas, particularly in the Amazon. One contributing factor may be the incomplete representation of the diurnal evolution of shallow cumulus (ShCu) clouds. This study aimed to enhance the understanding of the diurnal cycles of ShCu clouds—from formation to maturation and dissipation—over the Central Amazon (CAMZ). Using observational data from the Green Ocean Amazon 2014 (GoAmazon) campaign and large eddy simulation (LES) modeling, we analyzed the diurnal cycles of six selected pure ShCu cases and their composite behavior. Our results revealed a well-defined cycle, with cloud formation occurring between 10 and 11 local time (LT), maturity from 13 to 15 LT, and dissipation by 17–18 LT. The vertical extent of the liquid water mixing ratio and the intensity of the updraft mass flux were closely associated with increases in turbulent kinetic energy (TKE), enhanced buoyancy flux within the cloud layer, and reduced large-scale subsidence. We further analyzed the diurnal cycles of the convective available potential energy (CAPE), the convective inhibition (CIN), the Bowen ratio (BR), and the vertically integrated TKE in the mixed layer (ITKE-ML), exploring their relationships with the cloud base mass flux (Mb) and cloud depth across the six ShCu cases. ITKE-ML and Mb exhibited similar diurnal trends, peaking at approximately 14–15 LT. However, no consistent relationships were found between CAPE (or BR) and Mb. Similarly, comparisons of the cloud depth with CAPE, BR, ITKE-ML, CIN, and Mb revealed no clear relationships. Smaller ShCu clouds were sometimes linked to higher CAPE and lower CIN. It is important to emphasize that these findings are preliminary and based on a limited sample of ShCu cases. Further research involving an expanded dataset and more detailed analyses of the TKE budget and synoptic conditions is necessary. Such efforts would yield a more comprehensive understanding of the factors influencing ShCu clouds’ vertical development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Subsidence Measurements at the Teller 27 Field Site, Seward Peninsula, Alaska, 2021

In September 2021, differential GPS (dGPS) points were collected in areas of visual subsidence in the Teller watershed near mile marker 27 of the Bob Blodgett Nome-Teller Memorial Highway on the Seward Peninsula of Alaska. This study aimed to understand how structural permafrost loss is affecting meter-scale ground elevation. As ice thaws, heaving and slumping of the landscape creates new microtopographical features on the landscape that may alter surface hydrology and soil moisture and thus plant community composition and biogeochemical cycling. The Teller 27 field site is underlain with discontinuous permafrost, so the landscape features a variety of permafrost features. We collected ground elevations from eight subsidence areas in two forms: 1) Transects were sampled laterally through areas of ground subsidence to capture the high edges and low spots where the earth slumped. 2) The circumference of the subsidence area was sampled in order to measure the ground area affected by subsidence. These data were collected with an Emlid Reach RS2 dGPS and a base station. This dataset includes one *.csv of dGPS locations of the eight subsidence locations and one *.kml of measurement locations.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw: Modeling Archive

This Modeling Archive is in support of a NGEE Arctic publication in review "Drying of tundra landscapes will limit subsidence-induced acceleration of permafrost thaw". The study used a cryohydrology model to assess the potential risk of abrupt permafrost thaw triggered by melting ground ice, a key open question associated with permafrost response to a warming Arctic. The spatially resolved simulations are for a small catchment 465 ice-wedge polygons in polygonal tundra near Utqiagvik, Alaska in the high-emissions RCP8.5 climate scenario. The simulations are compared to runoff, evapotranspiration and subsidence in the current climate and agree well. The study used the ATS code configured as an intermediate-scale cryohydrology model (Advanced Terrestrial Simulator, Version v2). The archive includes input files for spinup (1985 to 2005) and projections (2006 to 2100 - the manuscript reports 2006 to 2098). Three of the projections include the effects of subsidence and microtopography change with different depth profiles of ice content corresponding to the median, 20th percentile and 80th percentile. The fourth projection has subsidence turnoff and uses the reference case (median) ice content. Spatially averaged or aggregated output variables are included for the projections. Selected checkpoint files for the projections and postprocessing scripts are also included. Files included are *.xml; *.h5; *.exo; *.py; *.sh, *. nb, mesh files, and one file in the original Excel format plus *.csv and *.pdf to conserve formatting contained in the original file. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

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↗