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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Assessing the limitations of commercial sensors and models for supporting marine carbon dioxide removal monitoring: a case study

Several unknowns remain surrounding marine Carbon Dioxide Removal (mCDR) monitoring, reporting, and verification (MRV) practices and capabilities. Current in-situ sensor technology is limited (primarily pH and pCO 2 ), requiring calculations and assumptions to estimate changes in carbonate chemistry parameters, including total alkalinity (TA). Considering that cost, energy consumption, and accuracy of commercial sensors can vary by orders of magnitude, understanding how well existing sensors perform in an mCDR context is important for this emerging community. Likewise, documenting sensor limitations and how relatively simple models can optimize sensor deployments will improve MRV efforts and support protocol development. Here we (1) compare performance a variety of commercially available sensors in a blind mesocosm experiment simulating ocean alkalinity enhancement (OAE), and how sensor performance impacted carbonate chemistry estimates; (2) evaluate if sensors can distinguish the OAE signal from natural variability during a small scale OAE field test in Sequim Bay, WA, USA, and (3) use an idealized ocean biogeochemistry model to explore optimal sensor network design based on (1) and (2). Our mesocosm results indicate that correctly constraining pH uncertainty will be critical for accurate TA estimates with current sensor technology compared to the less impactful variation caused by uncertainty in pCO 2 (pH data that are presented throughout are reported on the total scale (pH T ) unless otherwise noted). Our pilot field test demonstrated that sensors were capable of distinguishing mCDR signatures from natural variability under optimal real-world conditions. Idealized modeling simulations of the field test showed that a range of sparse and dense (3 to 100) sensors sampling areas of detectable increases will underestimate the net change in surface pH by at least 35–55%, at both realistic and highly elevated alkalinity input levels. We also highlight the limitations of current sensing technology for MRV, and the importance of ocean biogeochemistry models as critical tools for predicting when and where mCDR signals will be detectable using available sensors. Overall, our findings suggest that commercially available pCO 2 sensors and some pH sensors will form an important backbone for mCDR MRV tasks, though complete MRV characterization will require these data to be used in combination with other tools.

OAE↗
DOE OSTI2026

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

LANL Institutional Computing Close-out Report for Project t22_ocean_time_step

This year, my team made a very productive use of LANL IC time. Four papers were published that used IC resources, and two more are under review. These publications fall into three categories: 1.) Improving tide modeling in the global ocean by adding self attraction and loading (Barton et. al. 2022), ice shelf cavities (Pal et al. 2023), and local time stepping (Lilly et al. 2023). 2.) A new sea ice numerical formulation (Capodaglio 2023). 3.) Performance comparisons, methods, and test cases for ocean model development, This includes a verification suite for ocean models (Bishnu, submitted) and a comparison between Julia and Fortran (Strauss, submitted). A new time-stepping method was introduced in Calandrini et al. (2022). These publications, and the use of LANL IC resources, go hand-in-hand with our mentoring efforts to train young scientists. Two of the lead authors are graduate students who are conducting their PhD research: Kristin Barton at the University of Michigan and Jeremy Lilly at Oregon State. These are both their first publications, and they both have DOE funding and DOE mentors. In addition, Bishnu is a post-doctoral researcher at LANL; Strauss conducted his research as a senior in high school; and Capodaglio, Calandrini and Pal are all early-career scientists who were converted to staff in 2021 or 2022. Here we highlight two publications on improvements in tidal modeling: Barton et al. 2022 and Lilly et al. 2023.

58 GEOSCIENCES↗
DOE OSTI2021

Simulations With the Marine Biogeochemistry Library (MARBL)

Abstract The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO 2 flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO 2 is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr −1 and particulate export flux at 100 m of 7.1 Pg C yr −1 . The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗
DOE OSTI2023

Reducing impacts of artificial ponding in modeling salt marshes using a conductivity-formulated subgrid model

The landforms of salt marsh systems are rather complex, with meandering channel networks cutting through low-lying, extensive marsh platforms. Some small-scale topographic details are critical to flooding and draining processes but cannot be reflected in the Digital Elevation Model (DEM) due to limited measurement resolution and vegetation bias. Numerical models need a grid with high enough resolution to capture the small-scale connectivity to correctly model the small-scale flow processes, making the model computationally unaffordable. Here, in this study, we derived a subgrid model, which couples two flow components, one is the regular subgrid model for simulating the subgrid-resolved flows, and the other is a conductivity-formulated porous flow model for modeling small-scale gut flows unresolved at the subgrid level. The model validation was performed using an idealized case with a narrow slot that was unresolved in the coupled flow model. Model results were compared with a slot-resolved regular subgrid model where only surface flows are solved and suggested that, with a proper model parameterization, the coupled flow model is capable of reproducing flooding and draining processes similar to the slot-resolved subgrid model of surface flow. The model was applied to simulating flooding and draining processes in a meso-tidal salt marsh. The model results demonstrate that the coupled subgrid model can be used to reduce the impacts of artificial ponding that usually occur in numerical simulations of salt marshes with numerous narrow gullies and creeks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A practical method for modeling temporally-averaged ocean wave frequency-directional spectra for characterizing wave energy climates

Wave energy resource characterization investigations have been hindered because frequency-directional wave spectra are not broadly available for many coastal regions due to the insufficient spatial coverage of buoy observation networks and spectral wave model hindcast outputs. We address this problem by developing and validating a practical method for approximating temporally-averaged frequency-directional wave spectra averaged over periods of months or years from bulk wave parameters, which are available at a much greater coverage and resolution than frequency-directional wave spectra. While the temporally averaged frequency-directional wave spectrum over these periods cannot be used for analyzing a single sea state, it aggregates multiple sea states, identifies dominant wave energy systems representing wave climate, and resolves their spectral characteristics. Therefore, modelling temporally averaged frequency-directional wave spectra is of great value for planning and designing wave energy projects, e.g., resource characterization, site assessment, and conceptual design of wave energy converters. Temporally-averaged frequency-directional wave spectra and related important wave energy parameters approximated using this method are found to be more accurate than commonly used parametric wave spectrum models. This method can be applied to a wide range of wave climates given its universality and high accuracy. Also, the availability of bulk wave parameters from multi-decade high-resolution wave model hindcasts is increasing.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Storm Surge Modeling as an Application of Local Time-Stepping in MPAS-Ocean

This paper presents the first practical application of local time-stepping (LTS) schemes in the Model for Prediction Across Scales-Ocean (MPAS-O). We use LTS schemes in a single-layer, global ocean model that predicts the storm surge around the eastern coast of the United States during Hurricane Sandy. The variable-resolution meshes used are of unprecedentedly high resolution in MPAS-O, containing cells as small as 125 m wide in Delaware Bay. It is shown that a particular, third-order LTS scheme (LTS3) produces sea-surface height solutions that are of comparable quality to solutions produced by the classical four-stage, fourth-order Runge-Kutta method (RK4) with a uniform time step on the same meshes. Furthermore, LTS3 is up to 35% faster in the best cases considered, where the number of cells using the coarse time-step relative to those using the fine time-step is as low as 1:1. This shows that LTS schemes are viable for use in MPAS-O with the added benefit of substantially less computational cost. The results of these performance experiments inform us of the requirements for efficient mesh design and configuration of LTS regions for LTS schemes. In particular, we see that for LTS to be efficient on a given mesh, it is important to have enough cells using the coarse time-step relative to those using the fine time-step, typically at least 1:5 to see an increase in performance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Lab-Scale, Closed-Loop Experimental Characterization, Model Refinement, and Validation of a Hydrokinetic Energy-Harvesting Ocean Kite

This paper presents a study wherein we experimentally characterize the dynamics and control system of a lab-scale ocean kite, and then refine, validate, and extrapolate this model for use in a full-scale system. Ocean kite systems, which harvest tidal and ocean current resources through high-efficiency cross-current motion, enable energy extraction with an order of magnitude less material (and cost) than stationary systems with the same rated power output. However, an ocean kite represents a nascent technology that is characterized by relatively complex dynamics and requires sophisticated control algorithms. In order to characterize the dynamics and control of ocean kite systems rapidly, at a relatively low cost, the authors have developed a lab-scale, closed-loop prototyping environment for characterizing tethered systems, whereby 3D printed systems are tethered and flown in a water channel environment. While this system has been shown to be capable of yielding similar dynamic characteristics to some full-scale systems, there are also fundamental limitations to the geometric scales and flow speeds within the water channel environment, making many other real-world scenarios impossible to replicate from the standpoint of dynamic similarity. Here, to address these scenarios, we show how the lab-scale framework is used to refine and validate a scalable dynamic model of a tethered system, which can then be extrapolated to full-scale operation. In this work, we present an extensive case study of this model refinement, validation, and extrapolation on an ocean kite system intended for operation in the Gulf Stream or similar current environments.

42 ENGINEERING↗
DOE OSTI2021

Generation of low–latitude seamount–trapped waves: a case study of the Seychelles Plateau

Baroclinic seamount-trapped waves are thought to influence their surrounding ecosystem; however, trapped waves are not well-studied in near-equatorial settings, where stratification is strong and Burger numbers are large. Motivated by observations, we use daily output (2005-2009) from the global 0.1-degree Parallel Ocean Program Model (POP) to examine topographically trapped baroclinic waves around the Seychelles Plateau (S > 400) in the tropical Indian Ocean. These trapped waves are associated with velocity and temperature oscillations at periods of 15–16 days, similar to the dominant period of some equatorial Yanai waves. Energy flux maps using POP output suggest that quasi-biweekly equatorial Yanai waves excite trapped waves on the western and south-western flanks of the Seychelles Plateau, near the surface. Here, the anticyclonic energy flux associated with the trapped wave extends vertically throughout the water column and around most of the plateau circumference, diminishing on the eastern flank of the plateau. This work highlights the role that equatorial planetary waves and trapped waves play in facilitating energy redistribution, dissipation, and mixing in the tropical ocean.

58 GEOSCIENCES↗
DOE OSTI2023

Simulating Southern Ocean Aerosol and Ice Nucleating Particles in the Community Earth System Model Version 2

Southern Ocean (SO) low-level mixed phase clouds have been a long-standing challenge for Earth system models to accurately represent. While improvements to the Community Earth System Model version 2 (CESM2) resulted in increased supercooled liquid in SO clouds and improved model radiative biases, simulated SO clouds in CESM2 now contain too little ice. Previous observational studies have indicated that marine particles are major contributor to SO low-level cloud heterogeneous ice nucleation, a process that initiates a number of cloud processes that govern cloud radiative properties. In this study, we utilize detailed aerosol and ice nucleating particle (INP) measurements from two recent measurement campaigns to assess simulated aerosol abundance, number size distributions, and composition and INP parameterizations for use in CESM2. Our results indicate that CESM2 has a positive bias in simulated surface-level total aerosol surface area at latitudes north of 58°S. Measured INP populations were dominated by marine INPs and we present evidence of refractory INPs present over the SO assumed here to be mineral dust INPs. Results highlight a critical need to assess simulated mineral dust number and size distributions in CESM2 in order to adequately represent SO INP populations and their response to long-term changes in atmospheric transport patterns and land use change. Furthermore, we also discuss important cautions and limitations in applying a commonly used mineral dust INP parameterization to remote regions like the pristine SO.

58 GEOSCIENCES↗
DOE OSTI2023

Drivers of Dry Day Sensitivity to Increased CO 2

Persistent precipitation deficits are among the most impactful consequences of global warming. Here we focus on changes in the annual number of dry days (NDD) and in the annual maximum length of dry spells due to a quadrupling of atmospheric CO 2 . We use atmosphere-only simulations to decompose the projected changes into additive contributions. A fast adjustment leads to a global increase in NDD despite notable regional exceptions (e.g., South Asia and Sahel). The effect of the uniform component of the surface ocean warming is model-dependent but shapes the regional distribution of the NDD response in each model. Finally, the ocean warming pattern also contributes to large uncertainties, likely through contrasting changes in large-scale circulation. Our results thus highlight the complexity of the NDD response, with policy-relevant practical implications for mitigation and adaptation strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evaluating the effect of meso/submesoscale current–wave interactions on wave energy resource characterization at northeast U.S. coast

Wave energy is a promising renewable resource, but accurate assessment is difficult in regions with strong currents due to wave–current interactions (WCI). Here, this study develops a two-way coupled WCI model within the Coupled Ocean Atmosphere Wave Sediment Transport (COAWST) framework at 2 km resolution to improve wave energy characterization along the northeastern U.S. coast, including the Mid-Atlantic Bight and Gulf of Maine. The model integrates WaveWatchIII (WWIII) and the Regional Ocean Modeling System (ROMS) to enhance wave hindcasting by accounting for Doppler-shift, refraction, and nonlinear energy exchanges. Validation against buoy and satellite observations confirms model accuracy. Analysis shows that Doppler-shifting can alter wave power density by over 20%, while strong current gradients and shear distort wave crests via focusing/defocusing and stretching/squeezing, modifying wave direction and frequency. These processes together can induce wave power fluctuations of up to 40% on synoptic scales. Applying a 2.5 MW Ocean Energy Converter power matrix shows that WCI may change harvested energy by up to 100% in shallow-waters and 60% in deep-waters. These results underscore the importance of incorporating nonlinear WCI for reliable wave climate predictions and resource assessments in energetic coastal regions.

doppler-shift↗
DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part II: Modeling Study

Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.

Meteorology & Atmospheric Sciences↗