Creating a Sea-Level-Enabled E3SM: A critical capability for predicting coastal impacts
Explore the source record for details and available documents.
Earth Sciences↗
SEARCH · Search NASA
Results for “sea level”
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.
Polar Processes, Sea Level Rise, and Coastal Impacts
Explore the source record for details and available documents.
Earth Sciences↗
Polar Processes, Sea Level Rise, & Coastal Impacts
Explore the source record for details and available documents.
Earth Sciences↗
Sea Level-Driven Marsh Migration Results in Rapid Net Loss of Carbon
Sea level rise alters coastal carbon cycling by driving the rapid migration of coastal ecosystems, salinization of freshwater systems, and replacement of terrestrial forests with tidal wetlands. Wetland soils accumulate carbon (C) at faster rates than terrestrial soils, implying that sea level rise may lead to enhanced carbon accumulation. Here, we show that carbon stored in tree biomass greatly exceeds carbon stored in adjacent marsh soils so that marsh migration reduces total carbon stocks by ~50% in less than 100 years. Continued marsh soil carbon accumulation may eventually offset forest carbon loss, but we estimate that the time for replacement is similar to estimates of marsh survival (i.e., centuries), which suggests that forest C may never be replaced. These findings reveal a critical C source not included in coastal C budgets driven by migrating ecosystems and rapidly shifting allocations between carbon stored in soils and biomass.
58 GEOSCIENCES↗
Geodynamically corrected Pliocene shoreline elevations in Australia consistent with midrange projections of Antarctic ice loss
The Mid-Pliocene represents the most recent interval in Earth history with climatic conditions similar to those expected in the coming decades. Mid-Pliocene sea level estimates therefore provide important constraints on projections of future ice sheet behavior and sea level change but differ by tens of meters due to local distortion of paleoshorelines caused by mantle dynamics. We combine an Australian sea level marker compilation with geodynamic simulations and probabilistic inversions to quantify and remove these post-Pliocene vertical motions at continental scale. Dynamic topography accounts for most of the observed sea level marker deflection, and correcting for this effect and glacial isostatic adjustment yields a Mid-Pliocene global mean sea level of +16.0 (+10.4 to +21.5) m (50th/16th to 84th percentiles). Recalibration of recent high-end sea level projections using this revised estimate implies a more stable Antarctic Ice Sheet under future warming scenarios, consistent with midrange forecasts of sea level rise that do not incorporate a marine ice cliff instability.
54 ENVIRONMENTAL SCIENCES↗
ISMIP6-based Antarctic projections to 2100: simulations with the BISICLES ice sheet model
The contribution of the Antarctic Ice Sheet is one of the most uncertain components of sea level rise to 2100. Ice sheet models are the primary tool for projecting future sea level contribution from continental ice sheets. The Ice Sheet Model Intercomparison for the Coupled Model Intercomparison Phase 6 (ISMIP6) provided projections of the ice sheet contribution to sea level over the 21st century, quantifying uncertainty due to ice sheet model, climate model, emission scenario, and uncertain parameters. We present simulations following the ISMIP6 framework with the BISICLES ice sheet model and new experiments extending the ISMIP6 protocol to more comprehensively sample uncertainties in future climate, ice shelf sensitivity to ocean melting, and their interactions. These results contributed to the land ice projections of Edwards et al. (2021), which formed the basis of sea level projections for the Sixth Assessment Report of the Intergovernmental Panel on Climate Change (AR6). Our experiments show the important interplay between surface mass balance processes and ocean-driven melt in determining Antarctic sea level contribution. Under higher-warming scenarios, high accumulation offsets more ocean-driven mass loss when sensitivity to ocean-driven melt is low. Conversely, we show that when sensitivity to ocean warming is high, ocean melting drives increased mass loss despite high accumulation. Overall, we simulate a sea level contribution range across our experiments from 2 to 178 mm. Finally, we show that collapse of ice shelves due to surface warming increases sea level contribution by 25 mm relative to the no-collapse experiments, for both moderate and high sensitivity of ice shelf melting to ocean forcing.
54 ENVIRONMENTAL SCIENCES↗
Insights into the vulnerability of Antarctic glaciers from the ISMIP6 ice sheet model ensemble and associated uncertainty
Abstract. The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from −5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty.
54 ENVIRONMENTAL SCIENCES↗
Manuscript model outputs, model source code, and figure scripts: The role of geomorphology in mediating biomass allocation impacts on salt-marsh resilience and carbon accumulation
This data package provides model source code (C++), model outputs, and figure generation scripts (R) needed to reproduce the following manuscript: Bruns, Nicholas E., Genevieve L. Noyce, and Matthew L. Kirwan. "The Role of Geomorphology in Mediating Biomass Allocation Impacts on Salt-Marsh Resilience and Carbon Accumulation." Estuarine, Coastal and Shelf Science 327 (December 2025): 109549. https://doi.org/10.1016/j.ecss.2025.109549. This manuscript investigates how geomorphology mediates the impact of biomass allocation shifts on salt marsh persistence and carbon (C) sequestration under sea level rise. We use a 1-D soil-column model (Kirwan and Mudd 2012) to perform experiments across a range of static root:shoot ratios (RSR = 1-4) spanning observed values. The model explicitly simulates interactions between tidal inundation, productivity, inorganic sediment deposition, and organic matter accumulation. Experiments determine whether geomorphic feedbacks amplify, dampen, or leave unchanged the ecosystem response to biomass allocation shifts. A first experiment uses constant sea level rise (2.5 mm/yr) to examine equilibrium dynamics and their influence on carbon accumulation rates across different suspended sediment concentrations (SSC). A second set of experiments calculates threshold sea level rise rates for marsh drowning across RSR and SSC combinations. Final experiments apply accelerating sea level rise scenarios (2000-2200) derived from NOAA projections (Sweet et al. 2022) to generate an envelope of expected responses, quantifying the importance of biomass allocation shifts on marsh carbon accumulation and persistence. All experiments are repeated across SSC ranging from 5-50 mg/L to investigate how these interactions vary in micro-tidal marshes with different sediment supplies. Package contents: * README.txt with detailed description of package contents and instructions for reproducing manuscript figures and model outputs * R scripts for generating all manuscript figures * C++ baseline model code from Kirwan and Mudd (2012) * C++ source code for 4 experimental model variants used in the manuscript, extending above baseline code * Model input files (.csv, .txt) including sea level rise scenarios * Model output files (.txt) used in manuscript analyses Temporal coverage: Model simulations span years 1900-2200, with accelerating sea level rise scenarios for 2000-2200. Key variables: Root:shoot ratio, suspended sediment concentration, carbon accumulation rate, vertical accretion rate, marsh elevation, inundation depth, threshold sea level rise rate.
54 ENVIRONMENTAL SCIENCES↗
Climate-driven tradeoffs between landscape connectivity and the maintenance of the coastal carbon sink
Abstract Ecosystem connectivity tends to increase the resilience and function of ecosystems responding to stressors. Coastal ecosystems sequester disproportionately large amounts of carbon, but rapid exchange of water, nutrients, and sediment makes them vulnerable to sea level rise and coastal erosion. Individual components of the coastal landscape (i.e., marsh, forest, bay) have contrasting responses to sea level rise, making it difficult to forecast the response of the integrated coastal carbon sink. Here we couple a spatially-explicit geomorphic model with a point-based carbon accumulation model, and show that landscape connectivity, in-situ carbon accumulation rates, and the size of the landscape-scale coastal carbon stock all peak at intermediate sea level rise rates despite divergent responses of individual components. Progressive loss of forest biomass under increasing sea level rise leads to a shift from a system dominated by forest biomass carbon towards one dominated by marsh soil carbon that is maintained by substantial recycling of organic carbon between marshes and bays. These results suggest that climate change strengthens connectivity between adjacent coastal ecosystems, but with tradeoffs that include a shift towards more labile carbon, smaller marsh and forest extents, and the accumulation of carbon in portions of the landscape more vulnerable to sea level rise and erosion.
54 ENVIRONMENTAL SCIENCES↗
Response of salt intrusion in a tidal estuary to regional climatic forcing
Abstract Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels.
54 ENVIRONMENTAL SCIENCES↗
Investigating coastal backwater effects and flooding in the coastal zone using a global river transport model on an unstructured mesh
Abstract. Coastal backwater effects are caused by the downstream water level increase as a result of elevated sea level, high river discharge and their compounding influence. Such effects have crucial impacts on floods in densely populated regions but have not been well represented in large-scale river models used in Earth system models (ESMs), partly due to model mesh deficiency and oversimplifications of river hydrodynamics. Using two mid-Atlantic river basins as a testbed, we perform the first attempt to simulate the backwater effects comprehensively over a coastal region using the MOSART river transport model under an ESM framework, i.e., Energy Exascale Earth System Model (E3SM) configured on a regionally refined unstructured mesh, with a focus on understanding the backwater drivers and their long-term variations. By including sea level variations at the river downstream boundary, the model performance in capturing backwaters is greatly improved. We also propose a new flood event selection scheme to facilitate the decomposition of backwater drivers into different components. Our results show that while storm surge is a key driver, the influence of extreme discharge cannot be neglected, particularly when the river drains to a narrow river-like estuary. Compound flooding, while not necessarily increasing the flood peaks, exacerbates the flood risk by extending the duration of multiple coastal and fluvial processes. Furthermore, our simulations and analysis highlight the increasing strength of backwater effects due to sea level rise and more frequent storm surge during 1990–2019. Thus, backwaters need to be properly represented in ESMs to improve the predictive understanding of coastal flooding.
54 ENVIRONMENTAL SCIENCES↗
Response Of Salt Intrusion In A Tidal Estuary To Regional Climatic Forcing (FVCOM Output for the baseline scenario)
Salinity distribution in a large tidal estuary is subject to estuarine adjustment under the influences of multiple physical drivers such as freshwater pulses and sea level rise, and is crucial to upstream water quality, aquaculture, and ecosystem functions of the estuary. To better understand the estuarine salinity response to climate change, the unstructured-grid Finite Volume Community Ocean Model was implemented to simulate the salt intrusion in the Delaware Bay Estuary. The model was first validated by multiple observational data sets and subsequently applied in an idealized setting to examine the response of salt front to freshwater pulses in high flow conditions, followed by a long-term drought condition supported by a multi-decadal streamflow drought analysis in the estuary. The model results showed that after the freshwater pulses the salt front location moved further upstream with sea level rise. Under the simulated long-term drought condition, the adjustment timescale of salt intrusion varies nonlinearly with sea level rise. With a significant increase in sea level rise, the adjustment timescale starts to decrease. This shift suggests a transition into a different regime where the estuary becomes more stratified, as indicated by an increasing bulk Simpson number with rising sea levels. These findings have been summarized into a journal manuscript. The data included in this information release request are FVCOM model output associated with the numerical simulations described in the journal manuscript.
Liu, Jinliang↗
Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice
We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.
58 GEOSCIENCES↗
Cosmogenic activation of silicon
In the production of 3 H, 7 Be, and 22 Na by interactions of cosmic-ray particles with silicon can produce radioactive backgrounds in detectors used to search for rare events. Through controlled irradiation of silicon CCDs and wafers with a neutron beam that mimics the cosmic-ray neutron spectrum, followed by direct counting, we determined that the production rate from cosmic-ray neutrons at sea level is (112 ± 24) atoms/( kg day ) for 3 H, (8.1 ± 1.9) atoms/( kg day ) for 7 Be, and (43.0 ± 7.2) atoms/( kg day ) for 22 Na . Complementing these results with the current best estimates of activation cross sections for cosmic-ray particles other than neutrons, we obtain a total sea-level cosmic-ray production rate of (124 ± 25) atoms/(kg day) for 3 H, (9.4 ± 2.0) atoms/( kg day ) for 7 Be, and (49.6 ± 7.4) atoms/( kg day ) for 22 Na. These measurements will help constrain background estimates and determine the maximum time that silicon-based detectors can remain unshielded during detector fabrication before cosmogenic backgrounds impact the sensitivity of next-generation rare-event searches.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A fast particle-based approach for calibrating a 3-D model of the Antarctic ice sheet
We consider the scientifically challenging and policy-relevant task of understanding the past and projecting the future dynamics of the Antarctic ice sheet. The Antarctic ice sheet has shown a highly nonlinear threshold response to past climate forcings. Triggering such a threshold response through anthropogenic greenhouse gas emissions would drive drastic and potentially fast sea level rise with important implications for coastal flood risks. Previous studies have combined information from ice sheet models and observations to calibrate model parameters. These studies have broken important new ground but have either adopted simple ice sheet models or have limited the number of parameters to allow for the use of more complex models. These limitations are largely due to the computational challenges posed by calibration as models become more computationally intensive or when the number of parameters increases. Here, we propose a method to alleviate this problem: a fast sequential Monte Carlo method that takes advantage of the massive parallelization afforded by modern high-performance computing systems. We use simulated examples to demonstrate how our sample-based approach provides accurate approximations to the posterior distributions of the calibrated parameters. The drastic reduction in computational times enables us to provide new insights into important scientific questions, for example, the impact of Pliocene era data and prior parameter information on sea level projections. These studies would be computationally prohibitive with other computational approaches for calibration such as Markov chain Monte Carlo or emulation-based methods. We also find considerable differences in the distributions of sea level projections when we account for a larger number of uncertain parameters. For example, based on the same ice sheet model and data set, the 99th percentile of the Antarctic ice sheet contribution to sea level rise in 2300 increases from 6.5 m to 13.1 m when we increase the number of calibrated parameters from three to 11. With previous calibration methods, it would be challenging to go beyond five parameters. Here, this work provides an important next step toward improving the uncertainty quantification of complex, computationally intensive and decision-relevant models.
54 ENVIRONMENTAL SCIENCES↗
Data-model files associated with the manuscript "Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California"
This package contains the data, simulation setups, notebooks and figures used in “Modeling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California” (Xu et al., 2025). In this study, we selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise on coastal hydrology using the process-based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site-specific data. We designed a novel modeling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modeled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected sea level rise over the next few decades. The low-lying Yampah Marsh, the reference site, is likely to be inundated due to future sea level rise when highest tides arrive; while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.Several files can be found from this data package.1. README.md: This file describes the title, journal, co-authors, abstract, repository structure and model version.2. Simulation_Setups.zip: The file contains the model configuration files (XML format) for ATS. 3. Notebooks.zip: The file contains the Jupyter notebooks for generating the pre- and post-restoration meshes and the meshes of future scenarios. 4. Figures.zip: The file contains the figures used in the manuscript.5. Data.zip: The file contains the data used to drive the model simulations, including watershed and wetlands boundaries, mesh files and references to additional datasets (e.g., meteorological forcing, tidal dataset, DEMs, land cover, soil properties). Also, it contains water level observations at the restored wetland.
54 ENVIRONMENTAL SCIENCES↗
Save the Mekong Delta from drowning
Owing to only a few decades of human influence and unsustainable management of the Mekong River basin’s natural resources, the Mekong Delta is receding rapidly. Most of the delta landform, home to 17 million people and an economic powerhouse, could slip below sea level by 2100. Avoiding such a catastrophic impact will require concerted actions that acknowledge root causes for land loss and the global importance of the delta landform. Deltas persist and grow if sediment supply from an upstream river basin builds delta land at the same or greater rates than land is submerged by relative sea level rise and erosion. Further, with more rapid sea level rise, more sediment resources are needed to maintain the current extent of the delta. Only improved coordination of governance and investments, informed by science, will provide the delta with those critical resources.
54 ENVIRONMENTAL SCIENCES↗
A Conterminous USA-Scale Map of Relative Tidal Marsh Elevation
Abstract Tidal wetlands provide myriad ecosystem services across local to global scales. With their uncertain vulnerability or resilience to rising sea levels, there is a need for mapping flooding drivers and vulnerability proxies for these ecosystems at a national scale. However, tidal wetlands in the conterminous USA are diverse with differing elevation gradients, and tidal amplitudes, making broad geographic comparisons difficult. To address this, a national-scale map of relative tidal elevation ( Z * MHW ), a physical metric that normalizes elevation to tidal amplitude at mean high water (MHW), was constructed for the first time at 30 × 30-m resolution spanning the conterminous USA. Contrary to two study hypotheses, watershed-level median Z * MHW and its variability generally increased from north to south as a function of tidal amplitude and relative sea-level rise. These trends were also observed in a reanalysis of ground elevation data from the Pacific Coast by Janousek et al. (Estuaries and Coasts 42 (1): 85–98, 2019). Supporting a third hypothesis, propagated uncertainty in Z * MHW increased from north to south as light detection and ranging (LiDAR) errors had an outsized effect under narrowing tidal amplitudes. The drivers of Z * MHW and its variability are difficult to determine because several potential causal variables are correlated with latitude, but future studies could investigate highest astronomical tide and diurnal high tide inequality as drivers of median Z * MHW and Z * MHW variability, respectively. Watersheds of the Gulf Coast often had propagated Z * MHW uncertainty greater than the tidal amplitude itself emphasizing the diminished practicality of applying Z * MHW as a flooding proxy to microtidal wetlands. Future studies could focus on validating and improving these physical map products and using them for synoptic modeling of tidal wetland carbon dynamics and sea-level rise vulnerability analyses.
Holmquist, James R. (ORCID:0000000325466766)↗