Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean Models”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Results of the second Ice Shelf–Ocean Model Intercomparison Project (ISOMIP+)

Ocean-driven basal melting of Antarctic ice shelves plays an important role in the mass loss of the Antarctic Ice Sheet. Ice shelf cavity-resolving ocean models are a valuable tool for understanding ice shelf-ocean interactions and for simulating projections of ice shelf and ocean states under future climate. Designed to assess the current state of ice shelf–ocean modelling, the second Ice Shelf–Ocean Model Intercomparison Project, ISOMIP+, consists of 12 ocean model configurations submitted with a common, idealised experimental setup. Here, we focus on the experiments Ocean0–2 (Asay-Davis et al., 2016), which are ocean models with idealised, static ice shelf geometries, but where the ocean reaches a balance with prescribed far-field ocean conditions. Different thermal transfer coefficient values (ranging from 0.011 to 0.2) are used for each model in the melting parameterisation to achieve a common, tuned melt rate since the models cover a range of types of vertical coordinates, ice–ocean boundary layer treatments, and numerical schemes. These model differences lead to spread in the resultant ocean properties, circulation, boundary-layer structure and spatial distribution of melting. We also highlight similarities between models, such as a shared linear relationship across most models between melt rate and overturning and barotropic streamfunctions during the spin-up and spin-down, demonstrating a robust relationship between melt and circulation across models and forcing conditions. The ISOMIP+ results provide a systematic comparison of ice shelf cavity-capable ocean models. However, we also demonstrate the need for realistic ice shelf–ocean model intercomparison projects (some already underway) to assess model biases and inter-model variation against sparse observations. Further research is needed to understand the differences between models and further improve our modelled representations of the ice–ocean boundary layer and ice shelf cavity circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Data-driven global ocean modeling for seasonal to decadal prediction

Accurate modeling of ocean dynamics is crucial for enhancing our understanding of complex ocean circulation processes, predicting climate variability, and tackling challenges posed by climate change. Although great efforts have been made to improve traditional numerical models, predicting global ocean variability over multiyear scales remains challenging. Here, we propose ORCA-DL (Oceanic Reliable foreCAst via Deep Learning), a data-driven three-dimensional ocean model for seasonal to decadal prediction of global ocean dynamics. ORCA-DL accurately simulates the three-dimensional structure of global ocean dynamics with high physical consistency and outperforms state-of-the-art numerical models in capturing extreme events, including El Niño–Southern Oscillation and upper ocean heat waves. Moreover, ORCA-DL stably emulates ocean dynamics at decadal timescales, demonstrating its potential even for skillful decadal predictions and climate projections. Our results demonstrate the high potential of data-driven models for providing efficient and accurate global ocean modeling and prediction.

Science & Technology - Other Topics↗
DOE OSTI2025

Ocean Model for E3SM Global Applications (OMEGA)

This ocean model is the next generation version of the previous LANL developed MPAS-Ocean (Model for Prediction Across Scales) Ocean model. It is specifically designed to be used effectively on high performance computing, in particular GPU enabled architectures.

Van Roekel, Luke↗
DOE OSTI2024

A Verification Suite of Test Cases for the Barotropic Solver of Ocean Models

Abstract The development of any atmosphere or ocean model warrants a suite of test cases (TCs) to verify its spatial and temporal discretizations, order of accuracy, stability, reproducibility, portability, scalability, etc. In this paper, we present a suite of shallow water TCs designed to verify the barotropic solver of atmosphere and ocean models. These include the non‐dispersive coastal Kelvin wave; the dispersive inertia‐gravity wave; the dispersive planetary and topographic Rossby waves; the barotropic tide; and a non‐linear manufactured solution. These TCs check the implementation of the linear pressure gradient term; the linear constant or variable‐coefficient Coriolis and bathymetry terms; and the non‐linear advection terms. Simulation results are presented for a variety of time‐stepping methods as well as two spatial discretizations: a mimetic finite volume method based on the TRiSK scheme, and a high‐order discontinuous Galerkin spectral element method. The experimental procedure for conducting these numerical experiments is detailed. It underscores several key considerations that vary depending on the chosen spatial discretization method. Finally, convergence studies of every TC are conducted with refinement in both space and time, only in space, and only in time. The convergence slopes match the expected theoretical predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Experimental design for the Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2)

The Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2) is a natural progression of previous and ongoing model intercomparison exercises that have focused on the simulation of ice-sheet and ocean processes in Antarctica. The previous exercises motivate the move towards realistic configurations, as well as more diverse model parameters and resolutions. The main objective of MISOMIP2 is to investigate the performance of existing ocean and coupled ice-sheet–ocean models in a range of Antarctic environments through comparisons to observational data. We will assess the status of ice-sheet–ocean modelling as a community and identify common characteristics of models that are best able to capture observed features. As models are highly tuned based on present-day data, we will also compare their sensitivity to prescribed abrupt atmospheric perturbations leading to either very warm or slightly warmer ocean conditions compared to the present day. The approach of MISOMIP2 is to welcome contributions of models as they are, including global and regional configurations, but we request standardized variables and common grids for the outputs. We target the analysis at two specific regions, the Amundsen Sea and the Weddell Sea, since they describe two different ocean environments and have been relatively well observed compared to other areas of Antarctica. An observational “MIPkit” synthesizing existing ocean and ice-sheet observations for a common period is provided to evaluate ocean and ice-sheet models in these two regions.

58 GEOSCIENCES↗
DOE OSTI2026

The ocean model for E3SM global applications: Omega version 0.1.0 – a new high-performance computing code for exascale architectures

This paper introduces Omega, the Ocean Model for E3SM Global Applications. Omega is a new ocean model designed to run efficiently on high performance computing (HPC) platforms, including exascale heterogeneous architectures with accelerators, such as Graphics Processing Units (GPUs). Omega is written in C and uses the Kokkos performance portability library. These were chosen because they are well-supported and will help future-proof Omega for upcoming HPC architectures. Omega will eventually replace the Model for Prediction Across Scales-Ocean (MPAS-Ocean) in the US Department of Energy's (DOE's) Energy Exascale Earth System Model (E3SM). Omega runs on unstructured horizontal meshes with variable-resolution capability and implements the same horizontal discretization as MPAS-Ocean. This work documents the design and performance of Omega Version 0.1.0 (Omega-V0), which solves the shallow water equations with passive tracers and is the first step towards the full primitive equation ocean model. On Central Processing Units (CPUs), Omega-V0 is 1.4 times faster than MPAS-Ocean with the same configuration. Omega-V0 is more efficient on GPUs than CPUs on a per-watt basis – by a factor of 5.3 on Frontier and 3.6 on Aurora, two of the world's fastest exascale computers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗
DOE OSTI2023

Quantification of Physical and Numerical Mixing in a Coastal Ocean Model Using Salinity Variance Budgets

Numerical mixing, the spurious mixing primarily generated by the discretization of advection, is often significant in estuarine and coastal models due to sharp, energetic fronts. In this study, we compare on- and offline estimates of numerical mixing in a submesoscale-resolving realistic simulation of the ocean state over the Texas-Louisiana continental shelf. While offline estimates of numerical mixing differ from online estimates, offline methods may be the only analysis available. This study offers insight into the differences between the on- and offline methods. We use two methods to estimate numerical mixing offline, based on salinity squared $s^2$ and volume-mean salinity variance $s^{\prime^2}$. Numerical mixing estimated from the $s^{\prime^2}$ budget is generally within 60\% of the magnitude for the online method but captures the temporal variability well. However, the $s^2$ budget compares poorly due to larger truncation errors associated with the tendency and advection terms, which can be reduced by increasing the model output frequency. We also investigate the effects of horizontal resolution on numerical mixing using a two-way nested grid. The volume-integrated numerical mixing constitutes 57\% of the bulk physical mixing -- the mixing prescribed by the turbulence closure scheme -- in the coarse model and may exceed the physical mixing by half an order of magnitude. We find that numerical mixing is reduced by 35\% on average in the nested model, less than expected based on scaling of the numerical mixing for an upwind advection scheme, likely due to new dynamical processes that emerge in the nested simulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Estimating Ocean Observation Impacts on Coupled Atmosphere‐Ocean Models Using Ensemble Forecast Sensitivity to Observation (EFSO)

Ensemble Forecast Sensitivity to Observation (EFSO) is a technique that can efficiently identify the beneficial/detrimental impacts of every observation in ensemble-based data assimilation (DA). While EFSO has been successfully employed on atmospheric DA, it has never been applied to ocean or coupled DA due to the lack of a suitable error norm for oceanic variables. This study introduces a new density-based error norm incorporating sea temperature and salinity forecast errors, making EFSO applicable to ocean DA for the first time. We implemented the oceanic EFSO on the CFSv2-LETKF and investigated the impact of ocean observations under a weakly coupled DA framework. By removing the detrimental ocean observations detected by EFSO, the CFSv2 forecasts were significantly improved, showing the validation of impact estimation and the great potential of EFSO to be extended as a data selection criterion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Streamlining Ocean Dynamics Modeling with Fourier Neural Operators: A Multiobjective Hyperparameter and Architecture Optimization Approach

Training an effective deep learning model to learn ocean processes involves careful choices of various hyperparameters. We leverage DeepHyper’s advanced search algorithms for multiobjective optimization, streamlining the development of neural networks tailored for ocean modeling. The focus is on optimizing Fourier neural operators (FNOs), a data-driven model capable of simulating complex ocean behaviors. Selecting the correct model and tuning the hyperparameters are challenging tasks, requiring much effort to ensure model accuracy. DeepHyper allows efficient exploration of hyperparameters associated with data preprocessing, FNO architecture-related hyperparameters, and various model training strategies. We aim to obtain an optimal set of hyperparameters leading to the most performant model. Moreover, on top of the commonly used mean squared error for model training, we propose adopting the negative anomaly correlation coefficient as the additional loss term to improve model performance and investigate the potential trade-off between the two terms. The numerical experiments show that the optimal set of hyperparameters enhanced model performance in single timestepping forecasting and greatly exceeded the baseline configuration in the autoregressive rollout for long-horizon forecasting up to 30 days. Utilizing DeepHyper, we demonstrate an approach to enhance the use of FNO in ocean dynamics forecasting, offering a scalable solution with improved precision.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Effects of Atmosphere and Ocean Horizontal Model Resolution on Tropical Cyclone and Upper-Ocean Response Forecasts in Four Major Hurricanes

A coupled atmosphere–ocean model is necessary for tropical cyclone (TC) prediction to accurately characterize ocean feedback on atmospheric processes within the TC environment. Here, the ECMWF coupled global model is run at horizontal resolutions from 9 to 1.4 km in the atmosphere, as well as 25 and 8 km in the ocean, to identify how resolution impacts forecast accuracy of four observed major TCs in the Atlantic: Irma, Florence, Teddy, and Ida. Most of the resolutions used here are unprecedented for global models. GOES-16 and synthetic aperture radar (SAR) satellite images and best track data are used for atmospheric validation. Salinity and temperature observations from Air-Launched Autonomous Micro-Observer (ALAMO) floats are used to validate modeled upper-ocean response, including mixed layer deepening, sea surface cooling, and near-inertial waves in the wakes of TCs. Increasing atmospheric resolution leads to more realistic TC structure and stronger winds, significantly improving TC intensity forecasts and modestly improving track errors. Ocean resolution impacts the upper-ocean response but does not influence atmospheric forecasts for the fast-moving TCs considered here. Stronger mixing, sea surface cooling, and near-inertial oscillations are found for both higher atmosphere and ocean resolutions, provided the initial upper-ocean state is the same for the two ocean resolutions. Whether this agrees better with the ALAMO observations also depends on the realism of the initial upper-ocean state in the model, emphasizing the importance of ocean initialization for the accurate upper-ocean response. Overall, the model at all resolutions correctly predicts stronger mixing, surface cooling, and near-inertial oscillation amplitudes to the right of a TC center, as observed by ALAMO floats.

Atmosphere-ocean interaction↗
DOE OSTI2023

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model con figurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

58 GEOSCIENCES↗
DOE OSTI2023

Scalable self attraction and loading calculations for unstructured ocean tide models

Self attraction and earth-loading effects are important for accurately modeling global tides. A common approach of handling this forcing is to expand mass anomalies into spherical harmonics, which are scaled by load Love numbers to account for elastic earth deformation. We investigate two different approaches to perform these calculations for ocean models that employ unstructured meshes and distributed memory parallelization. The first approach leverages a highly efficient spherical harmonics library, but requires all-to-one and one-to-all communications and interpolation operations between the unstructured and a structured mesh. This approach is compared to a parallel algorithm that computes the spherical harmonic transformations directly on the unstructured mesh with an all-reduce communication. Here, our results show that although the unstructured mesh calculations are more expensive, the scalability of the unstructured mesh approach allows for more efficient spherical harmonics transforms for high-resolution meshes and large processor counts. This methodology enables the efficient inclusion of tidal dynamics large-scale Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ice-shelf freshwater triggers for the Filchner–Ronne Ice Shelf melt tipping point in a global ocean–sea-ice model

Abstract. Some ocean modeling studies have identified a potential tipping point from a low to a high basal melt regime beneath the Filchner–Ronne Ice Shelf (FRIS), Antarctica, with significant implications for subsequent Antarctic ice sheet mass loss. To date, investigation of the climate drivers and impacts of this possible event have been limited because ice-shelf cavities and ice-shelf melting are only now starting to be included in global climate models. Using a global ocean–sea-ice configuration of the Energy Exascale Earth System Model (E3SM) that represents both ocean circulations and melting within ice-shelf cavities, we explore freshwater triggers (iceberg melt and ice-shelf basal melt) of a transition to a high-melt regime at FRIS in a low-resolution (30 km in the Southern Ocean) global ocean–sea-ice model. We find that a realistic spatial distribution of iceberg melt fluxes is necessary to prevent the FRIS melt regime change from unrealistically occurring under historical-reanalysis-based atmospheric forcing. Further, improvement of the default parameterization for mesoscale eddy mixing significantly reduces a large regional fresh bias and weak Antarctic Slope Front structure, both of which precondition the model to melt regime change. Using two different stable model configurations, we explore the sensitivity of FRIS melt regime change to regional ice-sheet freshwater fluxes. Through a series of sensitivity experiments prescribing incrementally increasing melt rates from the smaller, neighboring ice shelves in the eastern Weddell Sea, we demonstrate the potential for an ice-shelf melt “domino effect” should the upstream ice shelves experience increased melt rates. The experiments also reveal that modest ice-shelf melt biases in a model, especially at coarse ocean resolution where narrow continental shelf dynamics are not well resolved, can lead to an unrealistic melt regime change at downstream ice shelves. Thus, we find that remote connections between melt fluxes at different ice shelves are sensitive to baseline model conditions. Our results highlight both the potential and the peril of simulating prognostic Antarctic ice-shelf melt rates in a low-resolution global model.

54 ENVIRONMENTAL SCIENCES↗