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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mechanisms of heat flux across the Southern Greenland continental shelf in 1/10° and 1/12° ocean/sea ice simulations

The increased presence of warm Atlantic water on the Greenland continental shelf has been connected to the accelerated melting of the Greenland Ice Sheet, particularly in the southwest and southeast shelf regions. Results from two high-resolution coupled ocean-sea ice simulations that utilized either the 1/10-degree Parallel Ocean Program (POP) or the 1/12-degree HYbrid Coordinate Ocean Model (HYCOM) are used to understand the flux of heat on and off the southern Greenland shelf. The analysis reveals that the region of greatest heat flux onto the shelf is southeast Greenland. On the southwestern shelf, heat is mainly exported from the shelf to the interior basins. We identify differences in the shelf break current structure and on-shelf heat content between the two simulations. Just south of the Denmark strait, there is a seasonally persistent pattern of multi-day variability in the cross-shelf heat flux in both simulations. In the POP simulation, this high-frequency signal results in net on-shore heat flux. In the HYCOM simulation, the signal is weaker and results in net off-shelf heat flux. This variability is consistent with Denmark Strait Overflow eddies traveling along the shelf break.

58 GEOSCIENCES↗
DOE OSTI2020

Validation of the High-Resolution Salish Sea Tidal Hydrodynamic Model

In this study, a tidal hydrodynamic model was developed and validated to simulate tidal currents in Puget Sound, Washington, to support tidal energy resource characterization using the unstructured-grid, Finite Volume Community Ocean Model (FVCOM). The Salish Sea tidal hydrodynamic model was driven by tides along two open boundaries at the entrance of the Strait of Juan de Fuca and north end of Georgia Strait, and river flows from 19 major rivers in the Salish Sea. To simulate the tidal current in Puget Sound, a high-resolution model grid is required to accurately represent the complex coastlines and bathymetry. The spatial resolution of the model grid varies from ~10 m near river boundaries and ~30 m in small tidal channels and estuaries to near 1000 m inside Georgia Strait and at the open boundaries. Model validation was carried out by comparing simulated and observed water levels at 12 tidal stations and currents at 135 Acoustic Doppler Current Profiler stations in the model domain. A set of model performance metrics, including root mean square error, scatter index, bias, and linear correlation coefficient, were used to quantify the model skills in simulating the tidal hydrodynamics in Puget Sound. Error statistics showed an overall good agreement between simulated and observed tidal elevations and currents, which demonstrated that the Puget Sound tidal model can be used to accurately characterize the tidal stream energy resource in Puget Sound.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Impact of a New Sea Ice Thermodynamic Formulation in the CESM2 Sea Ice Component

Abstract The sea ice component of the Community Earth System Model version 2 (CESM2) contains new “mushy‐layer” physics that simulates prognostic salinity in the sea ice, with consequent modifications to sea ice thermodynamics and the treatment of melt ponds. The changes to the sea ice model and their influence on coupled model simulations are described here. Two simulations were performed to assess the changes in the vertical thermodynamics formulation with prognostic salinity compared to a constant salinity profile. Inclusion of the mushy layer thermodynamics of Turner et al. (2013, https://doi.org/10.1002/jgrc.20171 ) in a fully coupled Earth system model produces thicker and more extensive sea ice in the Arctic, with relatively unchanged sea ice in the Antarctic compared to simulations using a constant salinity profile. While this is consistent with the findings of uncoupled ice‐ocean model studies, the role of the frazil and congelation growth is more important in fully coupled simulations. Melt pond drainage is also an important contribution to simulated ice thickness differences as also found in the uncoupled simulations of Turner and Hunke (2015; https://doi.org/10.1002/2014JC010358 ). However, it is an interaction of the ponds and the snow fraction that impacts the surface albedo and hence the top melt. The changes in the thermodynamics and resulting ice state modify the ice‐ocean‐atmosphere fluxes with impacts on the atmosphere and ocean states, particularly temperature.

Bailey, David A.↗
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Export of Ice Sheet Meltwater from Upernavik Fjord, West Greenland

Meltwater from Greenland is an important freshwater source for the North Atlantic Ocean, released into the ocean at the head of fjords in the form of runoff, submarine melt, and icebergs. The meltwater release gives rise to complex in-fjord transformations that result in its dilution through mixing with other water masses. The transformed waters, which contain the meltwater, are exported from the fjords as a new water mass Glacially Modified Water (GMW). Here we use summer hydrographic data collected from 2013 to 2019 in Upernavik, a major glacial fjord in northwest Greenland, to describe the water masses that flow into the fjord from the shelf and the exported GMWs. Using an optimum multi-parameter technique across multiple years we then show that GMW is composed of 57.8% ± 8.1% Atlantic Water (AW), 41.0% ± 8.3% Polar Water (PW), 1.0% ± 0.1% subglacial discharge, and 0.2% ± 0.2% submarine meltwater. We show that the GMW fractional composition cannot be described by buoyant plume theory alone since it includes lateral mixing within the upper layers of the fjord not accounted for by buoyant plume dynamics. Consistent with its composition, we find that changes in GMW properties reflect changes in the AW and PW source waters. Using the obtained dilution ratios, this study suggests that the exchange across the fjord mouth during summer is on the order of 50 mSv (1 Sv ≡ 10 6 m 3 s –1 ) (compared to a freshwater input of 0.5 mSv). This study provides a first-order parameterization for the exchange at the mouth of glacial fjords for large-scale ocean models.

58 GEOSCIENCES↗
DOE OSTI2021

FVCOM-ICM4.0

FVCOM-ICM is a water quality model that simulates the influence of nutrients and carbon pollution on biogeochemical oceanographic processes such as phytoplankton primary productivity, nutrient consumption, mineralization, decay and settling of particulate organic matter and computes resulting dissolved oxygen (DO) concentrations and pH levels. An “offline” approach of coupling Finite Volume Coastal Ocean Model (FVCOM) hydrodynamic solutions to CE-QUAL-ICM biogeochemical model was used. FVCOM and CE-QUAL-ICM are public domain Fortran Codes that were developed independently by US Army Corps of Engineers and University of Massachusetts. PNNL developed FVCOM-ICM that has the ability to conduct water quality calculations using a previously computed hydrodynamic solution. The biogeochemical calculations using CE-QUAL-ICM kinetics are conducted over the same finite-volume mesh and in the FVCOM framework.

Khangaonkar, Tarang↗
DOE OSTI2021

Coupled Ocean Wave-Atmosphere Models for Offshore Wind Energy Applications

Coupling ocean wave models to mesoscale atmospheric models is necessary to represent the effect of waves on wind turbine hub-height winds. In this report we provide a review of the most widely used ocean waves models and the phased-averaged spectral wave modeling paradigm that they are based on. Methodologies used to couple these wave models to mesoscale atmospheric models are described along with details of existing coupled modeling systems. We summarize impacts on offshore wind resource assessments that have been investigated with such coupled modeling systems to-date. Specifically, coupling in the North and Baltic Seas was shown to have a small negative effect on the offshore wind energy resource at weak-to-moderate wind speeds, with little impact at higher wind speeds. Finally, limitations of these existing coupled modeling systems and impact assessments are discussed, including examples of and the potential for improved parameterizations of the air-sea fluxes, and the use of fine-scaled simulations using phased-resolved wave models.

17 WIND ENERGY↗
DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗
DOE OSTI2025

The Imprint of Southern Ocean Storms on Modeled Surface Chlorophyll, Their Drivers and Satellite Biases

Southern Ocean (SO) phytoplankton chlorophyll is highly variable on sub-seasonal time scales. Although the SO is the windiest ocean basin globally, it is not conclusively understood how storms impact SO phytoplankton dynamics. Much of our existing knowledge stems from satellites, but biases due to data gaps from cloud cover and low solar angles remain unquantified. Here, we use ocean–sea-ice simulations with the Community Earth System Model to quantify the climatological 1997–2018 imprint of storms on chlorophyll and phytoplankton dynamics in the ice-free SO. Additionally, by comparing the full-field model output to synthetic satellite observations, we quantify sampling biases in satellite-derived estimates. We find that both the sign and the magnitude of the average surface chlorophyll imprint vary substantially across storms but last for at least 4 days after the storm passing. Based on our analysis, more than one third of the storms explain the majority of local non-seasonal chlorophyll variability, but satellite-derived storm imprints are often too large in magnitude. On the day of the storm passing, changes in vertical mixing predominantly cause surface chlorophyll anomalies, and reduced light availability due to enhanced cloud cover outweighs the enhanced nutrient availability due to entrainment. Interestingly, storms imprint differently on total net primary production than on surface chlorophyll, demonstrating the difficulty to derive carbon-cycle impacts from a surface-chlorophyll assessment. With SO future storm activity projected to increase, complementing satellite observations with other observing technologies, for example, profiling floats, is necessary to better constrain how storms impact biological carbon cycling in the SO.

Nissen, Cara [University of Colorado, Boulder, CO ↗
DOE OSTI2021

Mid-Holocene Antarctic sea-ice increase driven by marine ice sheet retreat

Abstract. Over recent decades Antarctic sea-ice extent has increased, alongside widespread ice shelf thinning and freshening of waters along the Antarctic margin. In contrast, Earth system models generally simulate a decrease in sea ice. Circulation of water masses beneath large-cavity ice shelves is not included in current Earth System models and may be a driver of this phenomena. We examine a Holocene sediment core off East Antarctica that records the Neoglacial transition, the last major baseline shift of Antarctic sea ice, and part of a late-Holocene global cooling trend. We provide a multi-proxy record of Holocene glacial meltwater input, sediment transport, and sea-ice variability. Our record, supported by high-resolution ocean modelling, shows that a rapid Antarctic sea-ice increase during the mid-Holocene (∼ 4.5 ka) occurred against a backdrop of increasing glacial meltwater input and gradual climate warming. We suggest that mid-Holocene ice shelf cavity expansion led to cooling of surface waters and sea-ice growth that slowed basal ice shelf melting. Incorporating this feedback mechanism into global climate models will be important for future projections of Antarctic changes.

58 GEOSCIENCES↗
DOE OSTI2021

Attribution of Snowpack Errors to Simulated Temperature and Precipitation in E3SMv1 Over the Contiguous United States

Abstract Snow water equivalent (SWE), temperature, and precipitation biases and trends are evaluated in the atmosphere‐land simulations of the Energy Exascale Earth System Model version 1 (E3SMv1) in comparison to the Community Earth System Model version 2 (CESM2) and two other models using the ground measurement‐based University of Arizona (UA) snow product. SWE, temperature, and precipitation biases are highest in magnitude in the Western contiguous United States (CONUS). SWE errors are attributed to temperature and precipitation through multiple linear regressions of normalized errors, the coefficients of which represent the sensitivities to temperature and precipitation errors. SWE errors are more sensitive to temperature errors throughout the CONUS. Model SWE and temperature trends are generally opposite from UA product trends in the Western CONUS. SWE trend errors are also attributed to temperature and precipitation trend errors using multiple linear regressions of normalized trend errors. SWE trend errors are more sensitive to those of precipitation at higher elevations (>1,500 m) in the Western CONUS in these simulations. Thus, the sensitivity to temperature and precipitation differ for SWE errors and its trend errors. Furthermore, the SWE trend errors are more sensitive to temperature and precipitation in the atmosphere‐ocean coupled simulations in which the atmosphere‐land is coupled to an active ocean model. These results suggest that both errors in simulated temperature and precipitation contribute to SWE errors.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

The Influence of Ocean Topography on the Upwelling of Carbon in the Southern Ocean

The physical circulation of the Southern Ocean sets the surface concentration and thus air-sea exchange of CO 2 . However, we have a limited understanding of the three-dimensional circulation that brings deep carbon-rich waters to the surface. Here, we introduce and analyze a novel high-resolution ocean model simulation with active biogeochemistry and online Lagrangian particle tracking. We focus our attention on a subset of particles with high dissolved inorganic carbon (DIC) that originate below 1,000 m and eventually upwell into the near-surface layer (upper 200 m). We find that 71% of the DIC-enriched water upwelling across 1,000 m is concentrated near topographic features, which occupy just 33% of the Antarctic Circumpolar Current. Once particles upwell to the near-surface layer, they exhibit relatively uniform pCO 2 levels and DIC decorrelation timescales, regardless of their origin. Furthermore, our results show that Southern Ocean bathymetry plays a key role in delivering carbon-rich waters to the surface.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Data for Post-Access Report: Tidal Energy Resource Characterization and Model Validation via the Assessment of Cost-effective Spatially Distributed Current Measurements from microFloat Swarms

This project evaluated how high-resolution, spatially distributed field data can be used to refine and validate site-scale hydrodynamic simulations of tidal channels. Use of such spatially-distributed field observations or site-scale hydrodynamic simulations will be needed for producing accurate predictions of tidal energy production over larger arrays of tidal turbines. The data contained in this archive includes field observations of water velocity from a series of short surveys in Agate Pass taken on August 20th, 2020 using station-keeping ADCP (acoustic doppler current profiler) measurements, drifting downward-looking ADCPs, and microFloats, gathered by researchers at APL-UW under a previous project. It also includes high-resolution hydrodynamic model runs for the site during that time generated under TEAMER support. An explanatory file (.pdf format) describes the survey and field data structure. Processed data is stored in .mat format. Survey summary: start time, number of floats, water level. microFloat data: - Cleaned data for each float in survey, - Time series of position (lat., lon.), depth (relative to surface), and horizontal velocity (u, v) Drifting ADCP data: - Cleaned data for each drifting ADCP in survey - Time series of position (lat., lon.), sample depth (relative to surface), and observed water velocity (u,v,w) Stationary ADCP data: - Cleaned data for each station-keeping data set in survey - Time series of position (lat., lon.), sample depth (relative to surface), and observed water velocity (u,v,w) FVCOM (Finite Volume Coastal Ocean Model) model outputs: Water level and velocity, as well as derived variables (e.g., power density) from the final refined model simulations for the Agate Pass subdomain. Contained in NetCDF and ASCII files for the duration of August 18th and 20th 2020 field surveys. Additional scripts for processing the data and generating the figures in the report are included.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

INSIGHTS INTO MARINE UNICELLULAR CYANOBACTERIAL AND NON-CYANOBACTERIAL DIAZOTROPHS THROUGH SINGLE-CELL ANALYSES

Nutrient availability affects primary productivity, and nitrogen (N) is the limiting nutrient over large regions of the ocean. Biological N fixation (BNF) is the conversion of atmospheric dinitrogen (N2) to bioavailable forms of N (ammonium and amino acids). A select group of Bacteria and Archaea, known as diazotrophs, are the only organisms capable of this conversion. BNF is a critical source of N in nutrient-depleted (oligotrophic) marine waters. As such, BNF studies are a fundamental component for better understanding primary productivity and carbon cycling, as well as parameterizing predictive future global ocean models. Unicellular diazotrophs are widespread in marine systems and, in some cases, important contributors of BNF. This dissertation investigates unicellular diazotrophs at the single-cell level, including: UCYN-A, a cyanobacterial diazotroph that lives in symbiosis with an algal host, non-cyanobacterial diazotrophs (NCDs), which are putative heterotrophs, as well as a specific NCD known as Gamma A, a commonly occurring organism throughout the tropics and subtropics that belongs to the class Gammaproteobacteria. The single-cell investigations used in these studies can provide high-resolution insights that are not possible from community-level analyses of BNF, such as the N2 fixation rate of individual cells and cell morphology. Singlecell nanoscale secondary ion mass spectrometry (nanoSIMS) analyses were used in chapter 2 and 3 to investigate the cell-specific N2 fixation rates of (2) UCYN-A in the viii cold waters of the Bering and Arctic Seas where N2 fixation was not previously believed to be possible, and (3) NCDs in the North Pacific Subtropical Gyre where we demonstrated the first direct evidence of NCD N2 fixation on open ocean particles. Chapter 4 focused on the uncultivated NCD, Gamma A, using geneFISH to visualize single cells in natural seawater samples. We successfully identified Gamma A by developing a fluorescent in situ hybridization (FISH) probe targeting the diazotrophic functional marker gene nifH, revealing Gamma A cells are most often attached to particles. Together these studies provide new insights into the distributions, activities, and physiologies of geographically widespread, unicellular diazotrophs, thus improving our understanding of their roles as sources of fixed N in the surface ocean.

Harding, K.↗
DOE OSTI2022

High-order multirate explicit time-stepping schemes for the baroclinic-barotropic split dynamics in primitive equations

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta approach, we propose two high-order multirate explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The proposed schemes allow for a large time step to be used for the three-dimensional baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast) mode, in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the baroclinic velocity is first computed by advancing the baroclinic mode and fluid thickness of the system with the large time-step and the assistance of some intermediate approximations of the barotropicmode obtained by substepping with the small time step; then the barotropic velocity is corrected by using the small time step to re-advance the barotropic mode un-der an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting are relieved via a reconciliation process with carefully calculated flux deficits. Here, two benchmark tests from the “MPAS-Ocean” platform are carried out to numerically demonstrate the performance and parallel scalability of the proposed SSPRK-SE schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang↗
DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗
DOE OSTI2025

Decomposing a Compound Flood Event in an Urban Pacific Northwest Estuary: Primary Drivers and Projections for the Future

The Duwamish River Estuary (DRE) of Washington is prone to compound flooding during atmospheric river (AR) events. The processes contributing to such flooding (coastal and fluvial) have remained opaque to municipalities that are increasingly impacted. Here, we conduct a suite of coupled atmosphere-hydrology-ocean model simulations with varying forcing combinations (tide, surge, and/or river discharge) to identify the primary drivers of compound flooding during a recent AR event. We also test year 2,100 climate forcing to project how flooding and drivers may change in the future for the same event. We identify a clear distinction between dynamics in the downstream, engineered portion of the DRE compared to the upstream, “natural” river. Downstream, tides dominate water levels but contributions from storm surge and nonlinear tide-surge interaction elevate tide-only high waters from no flooding to major flooding. Upstream, total water levels during the event are ∼6 cm higher than downstream due to an increasing influence of river discharge over surge and tides. Notably, nonlinear surge-river and tide-river interactions act to reduce upstream water levels up to 50% compared to estimates which linearly sum tides, surge, and river, likely reducing flood vulnerability. Under two future climate scenarios: one with only sea level rise (SLR) and another with SLR plus atmospheric warming, we find little change in mechanism contributions to water levels. Expanded flooding in both cases is largely due to SLR, as a ∼50% increase to river discharge under the warming scenario has no impact downstream and marginally increases (∼3 cm) water level upstream.

atmospheric river↗