Search NASA⌕ Search

SEARCH · Search NASA

Results for “Earth Systems Model”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

The Single Column Atmosphere Model Version 6 (SCAM6): Not a Scam but a Tool for Model Evaluation and Development

The Single Column Atmosphere Model (SCAM) is a single column model version of the Community Atmosphere Model (CAM). Here we describe the functionality and features of SCAM6, available as part of CAM6 in the Community Earth System Model, version 2 (CESM2). SCAM6 features a wide selection of standard cases, as well as the ability to easily configure a case specified by the user based on a particular point in a CAM 3-D simulation. This work illustrates how SCAM6 reproduces CAM6 results for physical parameterizations, mostly of moisture and clouds. We demonstrate how SCAM6 can be used for model development through different physics selections, as well as with parameter sweep experiments to highlight the sensitivity of cloud properties to the specification of the vapor deposition process in the cloud microphysics. Furthermore, we use SCAM6 to illustrate the sensitivity of CAM6 cloud radiative properties and precipitation to variable drop number (cloud microphysics properties). Finally, we illustrate how SCAM6 can be used to explore critical emergent processes such as cloud feedbacks and show that CAM6 cloud responses to surface warming in stratus and stratocumulus regimes are similar to those in CAM5. CAM6 has a larger response in the shallow cumulus regime than CAM5. CAM6 cloud feedbacks in the shallow cumulus regime are sensitive to turbulence parameters. SCAM6 is thus a valuable tool for model development, evaluation, and scientific analy sis and an important part of the model hierarchy in Community Earth System Model, version 2.

58 GEOSCIENCES↗
DOE OSTI2025

ILAMBv2.7 benchmarking results comparing E3SMv2.1 land-atmosphere coupled (BGCv2LNDATM) and stand alone land (ELM) simulations with CMIP6 emission driven historical simulations

This dataset contains land model benchmarking results for the Energy Exascale Earth System Model version 2.1 (E3SMv2.1), including outputs from both coupled biogeochemistry simulations and stand-alone land model simulations. These results are compared against several emission-driven historical simulations from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Benchmarking was conducted using the International Land Model Benchmarking (ILAMB) package, version 2.7 (ILAMBv2.7). CMIP6 model outputs were sourced from the Earth System Grid Federation (ESGF), while the E3SMv2.1 results were derived from raw model outputs. These outputs underwent processing steps such as time serialization, conservative regridding, and data standardization to ensure comparability. For spatial interpolation, the Earth System Modeling Framework (ESMF) tool, ESMF_RegridWeightGen, was employed to generate regridding weights, enabling the transformation of E3SM’s native cubed-sphere grid to a regular latitude-longitude grid.

Feng, Sha [PNNL]↗
DOE OSTI2019

Flood Inundation Generation Mechanisms and Their Changes in 1953-2004 in Global Major River Basins

Despite the serious threats posed by floods, the driving mechanisms of floods are still not well understood. Here, we apply a physically-based inundation model coupled with a river routing model (Model for Scale Adaptive River Transport, MOSART) within the Energy Exascale Earth System Model (E3SM) framework to investigate flood inundation dynamics. After calibration using observed streamflow and satellite-derived flood extent, the model is used to simulate global flood inundation from 1953 to 2004. The mean date and seasonality of annual maximum flood, defined based on flood extent, exhibit significant regional differences across 16 major basins. Generally, soil moisture and monthly maximum daily rainfall (MMR) are the dominant drivers of flood in tropical basins while monthly maximum daily snowmelt (MMS) is the dominant driver in high latitude basins. From 1953-1982 to 1975-2004, significant changes in flood generation mechanisms are found in some basins such as Amazon, Lena, Yenisey, and Kolyma. Analysis of the rainfall seasonality and water balance at grid scale reveals a stronger rainfall seasonality in the Amazon during the later period that increases the synchrony between extreme rainfall and wet soil. With high antecedent soil moisture coinciding with rainfall, MMR contributes more to flood in the later period. Fewer extreme rainfall events and increasing soil moisture reduced the contribution of MMR and increased the role of MMS in floods in the Lena and Yenisey basins, respectively. Lastly, increased soil moisture and frequency of large MMS reduced the contribution of the latter to floods in the Kolyma basin.

Mao, Yuna↗
DOE OSTI2022

The DOE E3SM v1.2 Cryosphere Configuration: Description and Simulated Antarctic Ice‐Shelf Basal Melting

Abstract The processes responsible for freshwater flux from the Antarctic Ice Sheet (AIS), ice‐shelf basal melting and iceberg calving, are generally poorly represented in current Earth System Models (ESMs). Here we document the cryosphere configuration of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) v1.2. This includes simulating Antarctic ice‐shelf basal melting, which has been implemented through simulating the ocean circulation within static Antarctic ice‐shelf cavities, allowing for the ability to calculate ice‐shelf basal melt rates from the associated heat and freshwater fluxes. In addition, we added the capability to prescribe forcing from iceberg melt, allowing for realistic representation of the other dominant mass loss process from the AIS. In standard resolution simulations (using a noneddying ocean) under preindustrial climate forcing, we find high sensitivity of modeled ocean/ice shelf interactions to the ocean state, which can result in a transition to a high basal melt regime under the Filchner‐Ronne Ice Shelf (FRIS), presenting a significant challenge to representing the ocean/ice shelf system in a coupled ESM. We show that inclusion of a spatially dependent parameterization of eddy‐induced transport reduces biases in water mass properties on the Antarctic continental shelf. With these improvements, E3SM produces realistic ice‐shelf basal melt rates across the continent that are generally within the range inferred from observations. The accurate representation of ice‐shelf basal melting within a coupled ESM is an important step toward reducing uncertainties in projections of the Antarctic response to climate change and Antarctica's contribution to global sea‐level rise.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Arctic tundra shrubification: a review of mechanisms and impacts on ecosystem carbon balance

Abstract Vegetation composition shifts, and in particular, shrub expansion across the Arctic tundra are some of the most important and widely observed responses of high-latitude ecosystems to rapid climate warming. These changes in vegetation potentially alter ecosystem carbon balances by affecting a complex set of soil–plant–atmosphere interactions. In this review, we synthesize the literature on (a) observed shrub expansion, (b) key climatic and environmental controls and mechanisms that affect shrub expansion, (c) impacts of shrub expansion on ecosystem carbon balance, and (d) research gaps and future directions to improve process representations in land models. A broad range of evidence, including in-situ observations, warming experiments, and remotely sensed vegetation indices have shown increases in growth and abundance of woody plants, particularly tall deciduous shrubs, and advancing shrublines across the circumpolar Arctic. This recent shrub expansion is affected by several interacting factors including climate warming, accelerated nutrient cycling, changing disturbance regimes, and local variation in topography and hydrology. Under warmer conditions, tall deciduous shrubs can be more competitive than other plant functional types in tundra ecosystems because of their taller maximum canopy heights and often dense canopy structure. Competitive abilities of tall deciduous shrubs vs herbaceous plants are also controlled by variation in traits that affect carbon and nutrient investments and retention strategies in leaves, stems, and roots. Overall, shrub expansion may affect tundra carbon balances by enhancing ecosystem carbon uptake and altering ecosystem respiration, and through complex feedback mechanisms that affect snowpack dynamics, permafrost degradation, surface energy balance, and litter inputs. Observed and projected tall deciduous shrub expansion and the subsequent effects on surface energy and carbon balances may alter feedbacks to the climate system. Land models, including those integrated in Earth System Models, need to account for differences in plant traits that control competitive interactions to accurately predict decadal- to centennial-scale tundra vegetation and carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Design and Impacts of Land-Biogenic-Atmosphere Coupling in the NASA-Unified WRF (NU-WRF) Modeling System

Land-Atmosphere coupling is typically designed and implemented independently for physical (e.g. water and energy) and chemical (e.g. biogenic emissions and surface depositions)-based models and applications. Differences in scale, data requirements, and physics thus limit the ability of Earth System models to be fully coupled in a consistent manner. In order for the physical-chemical-biological coupling to be complete, treatment of the land in terms of surface classification, condition, fluxes, and emissions must be considered simultaneously and coherently across all components. In this study, we investigate a coupling strategy for the NASA-Unified Weather Research and Forecasting (NU-WRF) model that incorporates the traditionally disparate fluxes of water and energy through NASA's LIS (Land Information System) and biogenic emissions through BEIS (Biogenic Emissions Inventory System) and MEGAN (Model of Emissions of Gases and Aerosols from Nature) into the atmosphere. In doing so, inconsistencies across model inputs and parameter data are resolved such that the emissions from a particular plant species are consistent with the heat and moisture fluxes calculated for that land cover type. In turn, the response of the atmospheric turbulence and mixing in the planetary boundary layer (PBL) acts on the identical surface type, fluxes, and emissions for each. In addition, the coupling of dust emission within the NU-WRF system is performed in order to ensure consistency and to maximize the benefit of high-resolution land representation in LIS. The impacts of those self-consistent components on' the simulation of atmospheric aerosols are then evaluated through the WRF-Chem-GOCART (Goddard Chemistry Aerosol Radiation and Transport) model. Overall, this ambitious project highlights the current difficulties and future potential of fully coupled. components. in Earth System models, and underscores the importance of the iLEAPS community in supporting improved knowledge of processes and innovative approaches for models and observations.

Tan, Qian↗
DOE OSTI2018

Impact of Physics Parameterization Ordering in a Global Atmosphere Model

Abstract Because weather and climate models must capture a wide variety of spatial and temporal scales, they rely heavily on parameterizations of subgrid‐scale processes. The goal of this study is to demonstrate that the assumptions used to couple these parameterizations have an important effect on the climate of version 0 of the Energy Exascale Earth System Model (E3SM) General Circulation Model (GCM), a close relative of version 1 of the Community Earth System Model (CESM1). Like most GCMs, parameterizations in E3SM are sequentially split in the sense that parameterizations are called one after another with each subsequent process feeling the effect of the preceding processes. This coupling strategy is noncommutative in the sense that the order in which processes are called impacts the solution. By examining a suite of 24 simulations with deep convection, shallow convection, macrophysics/microphysics, and radiation parameterizations reordered, process order is shown to have a big impact on predicted climate. In particular, reordering of processes induces differences in net climate feedback that are as big as the intermodel spread in phase 5 of the Coupled Model Intercomparison Project. One reason why process ordering has such a large impact is that the effect of each process is influenced by the processes preceding it. Where output is written is therefore an important control on apparent model behavior. Application of k ‐means clustering demonstrates that the positioning of macro/microphysics and shallow convection plays a critical role on the model solution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2020

The Earth Surface Mineral Dust Source Investigation: An Earth Science Imaging Spectroscopy Mission

The Earth Surface Mineral Dust Source Investigation, EMIT, is planned to operate from the International Space Station starting no earlier than the fall of 2021. EMIT will use visible to short wavelength infrared imaging spectroscopy to determine the mineral composition of the arid land dust source regions of the Earth to advance our knowledge of the radiative forcing effect of these aerosols. Mineral dust emitted into the atmosphere under high wind conditions is an element of the Earth system with many impacts to the Earth’s energy balance, atmosphere, surface, and oceans. The Earth’s mineral dust cycle with source, transport, and deposition phases are studied with advanced Earth System Models. Because the chemical composition, optical and surface properties of soil particles vary strongly with the mineral composition of the source, these models require knowledge of surface soil mineral dust source composition knowledge to accurately understand dust impacts on the Earth system now and in the future. At present, compositional knowledge of the Earth’s mineral dust source regions from existing data sets is uncertain as a result of limited measurements. EMIT will use spectroscopically derived surface mineral composition to update the prescribe boundary conditions for state-of-the-art Earth System Models. The EMIT-initialized models will be used to investigate the impact of direct radiative forcing in the Earth system that depends strongly on the composition of the mineral dust aerosols emitted into the atmosphere. These new measurements and related products will be used to address the EMIT science objectives and made available to the science community for additional investigations. An overview of the EMIT science, development, and mission is presented in this paper.

Green, Robert O.↗
DOE OSTI2020

Coastal Wetland Carbon Sequestration in a Warmer Climate (Final Report)

Coastal wetlands are global hotspots of carbon storage and locations where carbon and nitrogen cycles have a disproportionately large impact on land, water, and air in comparison to the area they occupy. The extremely high rates of carbon sequestration in these systems are the result of complex feedbacks between vegetation and the physical environment. The plant and microbial interactions that drive these feedbacks are absent in ecosystem- and global-scale Earth System models. The goal of this proposal was to address the uncertain response of the coastal carbon sink to climate change, a critical knowledge gap in coastal carbon research and a significant barrier to incorporating coastal wetlands into Earth System Models. The Salt Marsh Accretion Response to Temperature eXperiment (SMARTX) was designed to understand the ecosystem-scale consequences of warming and elevated CO₂ (eCO₂) in tidal wetlands at the coastal terrestrial-aquatic interface. We successfully designed and built a novel whole-ecosystem warming and eCO₂ experiment in a coastal wetland. The gradient design of the warming treatments (+0, +1.7, +3.4, +5.1 ⁰C) allowed us to discover unexpected non-linear and non-additive responses arising from plant-microbe interactions. We concluded that these non-linear responses to warming were the result of plants and microbes responding to temperature at different thresholds, with important consequences for forecasting terrestrial ecosystem responses to climate change and understanding global trends.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of E3SM Simulated Aerosols and Aerosol‐Cloud Interactions Across GCM and Convection‐Permitting Scales

This paper introduces an Earth system modeling testbed for predicting aerosols and aerosol‐cloud interactions (ACIs) at convection‐permitting scales. Using the Energy Exascale Earth System Model (E3SM) version 2 with a four‐mode Modal Aerosol Module, we conduct simulations at 3.25 km resolution on a regionally refined mesh (RRM) across four regions with distinct aerosol and cloud regimes. Results are compared with the standard 100 km E3SM configuration and evaluated against satellite, aircraft, and ground‐based observations. We find that increasing model resolution improves heavy precipitation simulation but amplifies positive bias in light drizzle at coarse resolution. These resolution‐induced changes affect cloud and aerosol properties to varying degrees across regions. Generally, cloud cover and liquid water path (LWP) show better agreement with satellite retrievals at 3.25 km, though surface‐based comparisons suggest otherwise. Aerosol composition remains poorly represented at both resolutions. The RRM increases Aitken mode aerosol number concentrations via enhanced new particle formation. However, accumulation mode aerosols are decreased at higher resolution as aerosol removals become more efficient. This partially contributes to fewer cloud condensation nuclei (CCN) and lower cloud droplet number concentrations (N d ), which produces larger model biases in some scenarios. These findings suggest that solely increasing horizontal resolution to kilometer scales is insufficient to broadly improve aerosol and cloud predictions without concurrent advancements in physical and chemical process representations. Nonetheless, the RRM moderately improves key ACI relationships such as CCN‐N d correlation, reflecting enhanced aerosol activation representation. The LWP‐N d relationship is also better captured by RRM, suggesting a better characterization of LWP adjustment.

Huang, Meng [Pacific Northwest National Laboratory↗
DOE OSTI2018

Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1)

Abstract. The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1∘ horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

Abstract A significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

58 GEOSCIENCES↗
DOE OSTI2020

Understanding processes that control dust spatial distributions with global climate models and satellite observations

Dust aerosol is important in modulating the climate system at local and global scales, yet its spatiotemporal distributions simulated by global climate models (GCMs) are highly uncertain. In this study, we evaluate the spatiotemporal variations of dust extinction profiles and dust optical depth (DOD) simulated by the Community Earth System Model version 1 (CESM1) and version 2 (CESM2), the Energy Exascale Earth System Model version 1 (E3SMv1), and the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) against satellite retrievals from Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), Moderate Resolution Imaging Spectroradiometer (MODIS), and Multi-angle Imaging SpectroRadiometer (MISR). We find that CESM1, CESM2, and E3SMv1 underestimate dust transport to remote regions. E3SMv1 performs better than CESM1 and CESM2 in simulating dust transport and the northern hemispheric DOD due to its higher mass fraction of fine dust. CESM2 performs the worst in the Northern Hemisphere due to its lower dust emission than in the other two models but has a better dust simulation over the Southern Ocean due to the overestimation of dust emission in the Southern Hemisphere. DOD from MERRA-2 agrees well with CALIOP DOD in remote regions due to its higher mass fraction of fine dust and the assimilation of aerosol optical depth. The large disagreements in the dust extinction profiles and DOD among CALIOP, MODIS, and MISR retrievals make the model evaluation of dust spatial distributions challenging. Our study indicates the importance of representing dust emission, dry/wet deposition, and size distribution in GCMs in correctly simulating dust spatiotemporal distributions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Predicting the earth's future

The development of earth system models that will simulate the past and present and provide predictions of future conditions is essential now that human activities have the potential to induce changes in the planetary environment. Critical aspects of global change include its pervasiveness and ubiquity, its distribution in several distinct time-scale bands, and the interactions between the atmosphere, ocean, land surface, and the terrestrial and marine biospheres. A model of the earth system on the scale of decades to centuries, developed by the Earth System Science Committee (NASA) with the strategy of dividing by time scale rather than discipline, is presented and the requirements for observations to support the implementation of the model are reviewed.

Dutton, J. A.↗
DOE OSTI2017

Impact of numerical choices on water conservation in the E3SM Atmosphere Model Version 1 (EAM V1)

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations for sea level rise projection. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1° horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model is negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in the new model results in a very thin model layer at the Earth’s surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for this model.

Zhang, Kai↗
DOE OSTI2025

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Geometric remapping of particle distributions in the Discrete Element Model for Sea Ice (DEMSI v0.0)

Abstract. A new sea ice dynamical core, the Discrete Element Model for Sea Ice (DEMSI), is under development for use in coupled Earth system models. DEMSI is based on the discrete element method, which models collections of ice floes as interacting Lagrangian particles. In basin-scale sea ice simulations the Lagrangian motion results in significant convergence and ridging, which requires periodic remapping of sea ice variables from a deformed particle configuration back to an undeformed initial distribution. At the resolution required for Earth system models we cannot resolve individual sea ice floes, so we adopt the sub-grid-scale thickness distribution used in continuum sea ice models. This choice leads to a series of hierarchical tracers depending on ice fractional area or concentration that must be remapped consistently. The circular discrete elements employed in DEMSI help improve the computational efficiency at the cost of increased complexity in the effective element area definitions for sea ice cover that are required for the accurate enforcement of conservation. An additional challenge is the accurate remapping of element values along the ice edge, the location of which varies due to the Lagrangian motion of the particles. In this paper we describe a particle-to-particle remapping approach based on well-established geometric remapping ideas that enforces conservation, bounds preservation, and compatibility between associated tracer quantities, while also robustly managing remapping at the ice edge. One element of the remapping algorithm is a novel optimization-based flux correction that enforces concentration bounds in the case of nonuniform motion. We demonstrate the accuracy and utility of the algorithm in a series of numerical test cases.

58 GEOSCIENCES↗
DOE OSTI2023

SERGHEI (SERGHEI-SWE) v1.0: a performance-portable high-performance parallel-computing shallow-water solver for hydrology and environmental hydraulics

The Simulation EnviRonment for Geomorphology, Hydrodynamics, and Ecohydrology in Integrated form (SERGHEI) is a multi-dimensional, multi-domain, and multi-physics model framework for environmental and landscape simulation, designed with an outlook towards Earth system modelling. At the core of SERGHEI's innovation is its performance-portable high-performance parallel-computing (HPC) implementation, built from scratch on the Kokkos portability layer, allowing SERGHEI to be deployed, in a performance-portable fashion, in graphics processing unit (GPU)-based heterogeneous systems. In this work, we explore combinations of MPI and Kokkos using OpenMP and CUDA backends. In this contribution, we introduce the SERGHEI model framework and present with detail its first operational module for solving shallow-water equations (SERGHEI-SWE) and its HPC implementation. This module is designed to be applicable to hydrological and environmental problems including flooding and runoff generation, with an outlook towards Earth system modelling. Its applicability is demonstrated by testing several well-known benchmarks and large-scale problems, for which SERGHEI-SWE achieves excellent results for the different types of shallow-water problems. Finally, SERGHEI-SWE scalability and performance portability is demonstrated and evaluated on several TOP500 HPC systems, with very good scaling in the range of over 20 000 CPUs and up to 256 state-of-the art GPUs.

58 GEOSCIENCES↗