Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ocean 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Wave-powered water pump for upwelling in aquaculture: Numerical model and ocean test

Wave-powered upwelling can increase the productivity and survivability of several aquaculture species. This enhancement is due to transporting cold, nutrient-rich ocean water, typically found lower in the water column, to the surface. Macroalgaes, like kelp, exhibit increased growth from these altered conditions. The University of New Hampshire’s (UNH) wave-powered water pump (wave pump) is a point absorber wave energy converter (WEC) that uses ocean waves to create relative motion between a spar buoy and a concentric float which drives an internal pump. A numerical model of the wave pump was developed using WEC-Sim to predict device performance in the ocean. Wave pump performance was evaluated during a five day ocean test near Appledore Island in Maine in March 2023, where volumetric flow rate, relative distance between spar and float, and wave conditions were measured. These data were then used for numerical model validation. The ocean deployment recorded the device’s performance in a variety of sea states, with average significant wave heights up to 0.7 m. The ocean test data were compared to the WEC-Sim numerical model of the device with favorable results. Average values of device stroke period, stroke height, and flow rate agreed between the ocean test and model data to within approximately 16 to 22%. Furthermore, the validated numerical model provides a valuable tool for improving the design and developing a commercial-scale, wave-powered water pump for use in aquaculture.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

REMORA v1.0

REMORA (previously known as ROMS-X) is a GPU-enabled regional ocean modeling code with the same core capabilities as the Regional Ocean Modeling System (ROMS) code.

Almgren, AnnS↗
DOE OSTI2021

Comparing numerical accuracy and stability for different horizontal discretizations in MPAS-Ocean

This manuscript investigates the effectiveness of two possible horizontal discretizations for the global ocean model MPAS-Ocean, both applied to Spherical Centroidal Voronoi Tessellations (SCVTs). The first discretization is TRiSK, a C-grid, finite-volume method, that possesses many desirable mimetic properties, but has a low order accuracy. The second discretization was introduced for the first time by Peixoto (2016), and consists of modifications to the TRiSK scheme designed to achieve at least first-order accuracy in the L ∞ norm, with the loss of some of the mimetic properties. Tests on shallow-water and primitive-equation models show that the scheme due to Peixoto is indeed more accurate, but presents stability issues with respect to TRiSK. Here, TRiSK is indeed found to be often more stable in time and more robust with respect to errors in the geometric properties of the grid.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

An unstructured C-grid type variational formulation for the sea ice dynamics

Historically, B-grid formulations of sea ice dynamics have been dominant because they have matched the grid type used by ocean models. The reason for the grid match is simple – it facilitates penetration of the curl of ice-ocean stress into the deep ocean with minimal numerical diffusivity because sea ice and ocean velocity are co-located. In recent years, as ocean models have increasingly progressed to C-grids, sea ice models have followed suit on quadrilateral meshes, but the implementation of an unstructured C-grid sea ice models is new. We present an unstructured C-grid discretization of the Elastic Viscous Plastic (EVP) rheology, where the velocity unknowns are discretized at the edges of the mesh cells with n-sides, where typically n is greater than or equal to four, rather than at the vertices, as in the B-grid. Our framework of choice is the Model for Prediction Across Scales (MPAS) within E3SM, the climate model of the U.S. Department of Energy, although our approach is general and could be applied to other models as well. While MPAS-Seaice is currently defined on a B-grid, MPAS-Ocean runs on a C-grid, hence interpolation operators are heavily used when coupled simulations are performed. In this work, we describe a mathematical formulation to transition the dynamics of MPAS-Seaice to a C-grid, in order to ultimately facilitate the coupling with MPAS-Ocean and reduce numerical errors associated with this communication. Numerical results are reported to highlight the features of the method.

58 GEOSCIENCES↗
DOE OSTI2023

The Effect of Tropical Pacific Air‐Sea Coupling on the Rainfall Response to Quadrupled CO 2 Forcing

We perform quadrupled CO 2 climate simulations with the Community Earth System Model version 1 (CESM1) to study how air-sea coupling affects the response of tropical rainfall under global warming. We use a hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows directed perpendicular to the wind, and the near-equator frictional flows directed in the same direction as the wind. We show that the Pacific Ocean's enhanced equatorial warming pattern (EEW) and equatorward ITCZ contraction observed in previous climate simulations emerge when the ocean model includes wind-driven Ekman and frictional flows. Furthermore, the near-equator frictional flow contributes more than half of the heat convergence in the equatorial Pacific Ocean. Finally, we show that although Ekman flow and near-equator frictional flow can both result in EEW, their coupled interactions with the Hadley circulation lead to opposite feedbacks on EEW's strength.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗
DOE OSTI2025

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗
DOE OSTI2025

A GPU‐Based Ocean Dynamical Core for Routine Mesoscale‐Resolving Climate Simulations

Abstract We describe an ocean hydrostatic dynamical core implemented in Oceananigans optimized for Graphical Processing Unit (GPU) architectures. On 64 A100 GPUs, equivalent to 16 computational nodes in current state‐of‐the‐art supercomputers, our dynamical core can simulate a decade of near‐global ocean dynamics per wall‐clock day at an 8‐km horizontal resolution; a resolution adequate to resolve the ocean's mesoscale eddy field. Such efficiency, achieved with relatively modest hardware resources, suggests that climate simulations on GPUs can incorporate fully eddy‐resolving ocean models. This removes a major source of systematic bias in current IPCC coupled model projections, the parameterization of ocean eddies, and represents a major advance in climate modeling. We discuss the computational strategies, focusing on GPU‐specific optimization and numerical implementation details that enable such high performance.

Silvestri, Simone [Massachusetts Institute of Tech↗
DOE OSTI2025

High-Resolution Modeling of the Gulf of Mexico using E3SM

Coastal ocean modeling is a high priority in the DOE‘s Energy Exascale Earth System Model (E3SM). The goal is to accurately predict the risk of damage to coastal resources and infrastructure due to a changing climate in the coming decades. North American coastal communities are areas of particular interest, as this fits under the topic of US national security and planning assessments in a changing climate. LANL Institutional Computing time for the Tier 1 allocation ”Coastal Ocean and Sea Ice Modeling” have been used for development and testing of numerical methods needed for E3SM coastal applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Rain-Induced Stratification of the Equatorial Indian Ocean and Its Potential Feedback to the Atmosphere

Surface freshening through precipitation can act to stably stratify the upper ocean, forming a rain layer (RL). RLs inhibit subsurface vertical mixing, isolating deeper ocean layers from the atmosphere. This process has been studied using observations and idealized simulations. The present ocean modeling study builds upon this body of work by incorporating spatially resolved and realistic atmospheric forcing. Fine-scale observations of the upper ocean collected during the Dynamics of the Madden-Julian Oscillation field campaign are used to verify the General Ocean Turbulence Model (GOTM). Spatiotemporal characteristics of equatorial Indian Ocean RLs are then investigated by forcing a 2D array of GOTM columns with realistic and well-resolved output from an existing regional atmospheric simulation. RL influence on the ocean-atmosphere system is evaluated through analysis of RL-induced modification to surface fluxes and sea surface temperature (SST). This analysis demonstrates that RLs cool the ocean surface on time scales longer than the associated precipitation event. A second simulation with identical atmospheric forcing to that in the first, but with rainfall set to zero, is performed to investigate the role of rain temperature and salinity stratification in maintaining cold SST anomalies within RLs. Approximately one third, or 0.1°C, of the SST reduction within RLs can be attributed to rain effects, while the remainder is attributed to changes in atmospheric temperature and humidity. The prolonged RL-induced SST anomalies enhance SST gradients that have been shown to favor the initiation of atmospheric convection. These findings encourage continued research of RL feedbacks to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Machine Learning Eliminates Reanalysis Warm Bias and Reveals Weaker Winter Surface Cooling Over Arctic Sea Ice

The surface energy budget governs Arctic sea-ice growth/melt, yet observations are sparse, and reanalysis data sets suffer from systematic biases. Here, we train a neural network with observational data to bias-correct hourly ERA5 fluxes over Arctic ice-covered regions (≥70°N; sea-ice concentration >80%) for 1994–2024. Training data cover two full seasonal cycles and different sea-ice regimes. The neural network reduces RMSE for net shortwave radiation by ∼40%, downward longwave radiation by ∼16% and the total surface energy budget by ∼55%, eliminating the wintertime warm bias of ∼4 K in ERA5. Wintertime surface cooling is reduced by ∼50%, yielding thermodynamic ice-growth estimates of ∼80–120 cm, consistent with SMOS–CryoSat satellite thickness increases and in contrast to the 150–200 cm growth implied by ERA5. Our bias-corrected data capture the observed clear/cloudy states of the winter boundary layer and can be used to study Arctic climatology, evaluate climate models and drive sea-ice-ocean models.

Hossain, Akil [Alfred Wegener Institute for Polar ↗
DOE OSTI2022

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analytical population dynamics underlying harmful algal blooms triggered by prey avoidance

A modified version of the NPZD ecosystem model is used to analytically examine the effects of predation avoidance, a possible mechanism for triggering harmful algal blooms (HAB). To resolve HAB development caused by predation avoidance, an additional phytoplankton functional group is considered, one that has slower nutrient uptake and better predation avoidance characteristics than the non-harmful phytoplankton group used in traditional NPZD models. Because the two phytoplankton groups (one non-harmful and one HAB) compete for only one resource within the same system, steady state (equilibrium) conditions cannot occur without the presence of zooplankton; only the non-harmful phytoplankton group, which defeats the HAB group in the resource competition, can survive in the equilibrium. The presence of sufficient zooplankton effectively acts to replenish the nutrient pool by consuming the non-harmful phytoplankton. When this occurs, two equilibrium states are found: one with both phytoplankton groups coexisting, and one that only includes the HAB group. Further, the condition required for equilibrium is that the total nitrogen within the system should be larger than a threshold determined by model coefficients. The threshold and feasibility of the equilibrium are sensitive to the relative HAB predation avoidance coefficient. If the coefficient is larger than the ratio of net growth rates between the HAB and non-harmful phytoplankton group, the threshold becomes infinite, and an equilibrium is not feasible. The time scale for the system to reach an equilibrium state that includes a HAB group is determined asymptotically. The dependence of a threshold condition as a controlling factor may explain the regime shift of dominant species causing HABs. The ecosystem model is fully implemented into the Regional Ocean Modeling System and applied to an idealized coastal embayment (with depths and geometry taken from San Francisco Bay) to show numerically the dominance of prey avoidance dynamics in a natural shallow water environment that includes advection and diffusion. The analytical results improve strategies for HAB modeling and provide guidance for setting model coefficients necessary to resolve a HAB event.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Understanding AMOC stability: the North Atlantic Hosing Model Intercomparison Project

Abstract. The Atlantic meridional overturning circulation (AMOC) is an important part of our climate system. The AMOC is predicted to weaken under climate change; however, theories suggest that it may have a tipping point beyond which recovery is difficult, hence showing quasi-irreversibility (hysteresis). Although hysteresis has been seen in simple models, it has been difficult to demonstrate in comprehensive global climate models. Here, we outline a set of experiments designed to explore AMOC hysteresis and sensitivity to additional freshwater input as part of the North Atlantic Hosing Model Intercomparison Project (NAHosMIP). These experiments include adding additional freshwater (hosing) for a fixed length of time to examine the rate and mechanisms of AMOC weakening and whether the AMOC subsequently recovers once hosing stops. Initial results are shown from eight climate models participating in the Sixth Coupled Model Intercomparison Project (CMIP6). The AMOC weakens in all models as a result of the freshening, but once the freshening ceases, the AMOC recovers in half of the models, and in the other half it stays in a weakened state. The difference in model behaviour cannot be explained by the ocean model resolution or type nor by details of subgrid-scale parameterisations. Likewise, it cannot be explained by previously proposed properties of the mean climate state such as the strength of the salinity advection feedback. Instead, the AMOC recovery is determined by the climate state reached when hosing stops, with those experiments where the AMOC is weakest not experiencing a recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estuarine hurricane wind can intensify surge-dominated extreme water level in shallow and converging coastal systems

Based on the projected increase in hurricane landfall frequency on the middle to lower US east coast, we examined the crucial role of the estuarine wind field in exacerbating coastal flooding. A regionally refined atmospheric and two high-resolution hydrology and ocean models are integrated to provide plausible and physically consistent ensembles of hurricane events and the associated flooding inside the Delaware Bay and River, a US mid-Atlantic estuary. Model results show that the hurricane propagation direction, estuarine geometry, remote surge from the open ocean, and direct nearshore upwind stress could magnify the flood magnitude. More specifically, inland-bound tracks that make landfall before reaching the mid-Atlantic coast produce a more significant surge within Delaware Bay than the shore-parallel tracks, where the estuarine wind direction plays the primary role in surge amplification. Ultimately, this study emphasized the need for integrated models to capture the nonlinear dynamics and interactions in flood hazard modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing Effects of Climate Change on Legacy Waste at the Enewetak Atoll

The Republic of the Marshall Islands (RMI) is in the central Pacific Ocean ~4,500 km west of Hawaii. The Enewetak Atoll, located in the northwest part of the RMI, was the site for 43 nuclear weapon tests between 1948 and 1958. Fallout and deposition from the tests contaminated the island surfaces, lagoon waters and sediment, and nearby ocean waters at the atoll. In the 1970s, a cleanup effort collected radioactive waste and placed it in the Cactus Crater on Runit Island (also called the Runit Dome). In December 2021, Congress directed the U.S. Department of Energy to study the impacts of climate change on the Runit Dome nuclear waste disposal site. Pacific Northwest National Laboratory (PNNL) assembled a multidisciplinary team of climate scientists, ocean modelers, environmental scientists, and health physicists to assess the likely effects of remaining radionuclides at the Enewetak Atoll. PNNL’s approach focused on effects of tropical cyclones that were postulated to mobilize and transport contaminated lagoon sediments and result in human and biota exposure. PNNL’s study estimated (1) the radionuclide source term, (2) the effects of climate change on severe storms, (3) mobilization and transport of radionuclides, and (4) radiation dose to humans and biota. Radionuclides in the lagoon and/or ocean waters of the Enewetak Atoll were characterized by the U.S. Atomic Energy Commission (AEC) in 1972, Woods Hole Oceanographic Institution in 2015, and Lawrence Livermore National Laboratory in 2018. The RMI Nationwide Radiological Study was conducted in the early 1990s for radionuclides remaining in island soils. The 1972 AEC survey remains the most comprehensive source of radionuclide data on lagoon sediments. Climate change modeling at a regional scale in the central Pacific Ocean is limited. PNNL climate scientists simulated severe historical storms postulated to occur both in a recent climate (2015) and in the future (2090) using the Advanced Research Weather Research and Forecasting (WRF-ARW) model, employing a pseudo-global-warming technique. A postulated complete, future failure of the Runit Dome was also considered. PNNL developed a high-resolution regional ocean hydrodynamics model covering the entire RMI extended economic zone using the Finite Volume Coastal Ocean Model (FVCOM). The FVCOM model was run using global reanalysis data for current climate and WRF-ARW simulation for the future climate. PNNL also developed a radionuclide fate and transport model using the FVCOM Integrated Compartment Model (FVCOM-ICM) to simulate the current and future mobilization and transport of radionuclides sorbed to lagoon sediments and the exchange of radionuclides between the water and sediment. FVCOM-ICM-predicted radionuclide concentrations were then used to estimate radiation dose to humans and biota at all islands of the Enewetak Atoll. Under current climate conditions, annual radiation exposures for the southern islands including Enewetak (Fred) and Medren (Elmer) were below the current U.S. standards. Radiation doses were somewhat elevated starting at Runit Island northward and westward to Enjebi Island (Janet). The islands in the northwest quadrant, particularly Bokoluo (Alice) and Bokombako (Belle), remain relatively contaminated. The islands in the southwestern quadrant have low contamination. The highest contribution to radiation doses comes from consumption of locally grown foods. Two radionuclides, 90Sr and 137Cs, contributed the greatest fraction for most terrestrial foods. In current climate conditions, the storms temporarily increased radionuclide concentrations in the lagoon waters, increasing the radiation dose slightly. In future conditions, doses are expected to be smaller, primarily because of the radioactive decay of the shorter-lived radioisotopes of 90Sr and 137Cs. This could make all islands in the far northwest of the atoll – except Bokombako (Belle) and perhaps Bokoluo (Alice) – suitable for residency. For the f

Prasad, Rajiv↗