Search NASA⌕ Search

SEARCH · Search NASA

Results for “TiDE 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Low-Reynolds-number oscillating boundary layers on adiabatic slopes

We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.

42 ENGINEERING↗
DOE OSTI2023

Barotropic tides in MPAS-Ocean (E3SM V2): impact of ice shelf cavities

Abstract. Oceanic tides are seldom represented in Earth system models (ESMs) owing to the need for high horizontal resolution to accurately represent the associated barotropic waves close to coasts. This paper presents results of tides implemented in the Model for Prediction Across Scales–Ocean or MPAS-Ocean, which is the ocean component within the U.S. Department of Energy developed Energy Exascale Earth System Model (E3SM). MPAS-Ocean circumvents the limitation of low resolution using unstructured global meshing. We are at this stage simulating the largest semidiurnal (M2, S2, N2) and diurnal (K1, O1) tidal constituents in a single-layer version of MPAS-O. First, we show that the tidal constituents calculated using MPAS-Ocean closely agree with the results of the global tidal prediction model TPXO8 when suitably tuned topographic wave drag and bottom drag coefficients are employed. Thereafter, we present the sensitivity of global tidal evolution due to the presence of Antarctic ice shelf cavities. The effect of ice shelves on the amplitude and phase of tidal constituents are presented. Lower values of complex errors (with respect to TPXO8 results) for the M2 tidal constituents are observed when the ice shelf is added in the simulations, with particularly strong improvement in the Southern Ocean. Our work points towards future research with varying Antarctic ice shelf geometries and sea ice coupling that might lead to better comparison and prediction of tides and thus better prediction of sea-level rise and also the future climate variability.

58 GEOSCIENCES↗
DOE OSTI2025

Applying Transfer Learning for Street-Scale Nuisance Flood Forecasting in Coastal-Urban Cities

An important challenge with Machine Learning (ML) is its transferability; that is, whether a ML model trained on one set of data can be applied to a second set of data without requiring a full re-training of the model. Transfer Learning (TL) addresses this challenge by transferring knowledge learned in the source domain (the data it was trained on) to the target domain (a second set of data that is statistically different but related, which the model was not trained on). This study investigates the use of TL for street-scale nuisance flood forecasting by exploring whether a ML model trained on data collected for one set of streets can effectively forecast flooding for another set of streets in the same city using TL. The envisioned use case is a city deploying a new flood depth monitoring sensor on a street and using TL to apply a ML model, trained on sensor data from an existing flood depth sensor network, to this new street. Eventually, the new flood depth sensor will have a sufficient dataset for training its own ML model, but TL can be used to fill the gap in time while this new dataset is being generated. This method is explored using a Long Short-Term Memory (LSTM) model trained on data for the flood-prone streets of Norfolk City, Virginia. The data used for training includes environmental time series (rainfall, tide), topographic features (Digital Elevation Model (DEM), Topographic Wetness Index (TWI), Depth To Water (DTW)), and street-scale flood depth time series obtained from a high-fidelity physics-based model, acting as a synthetic street-scale stream depth sensor dataset since actual stream depth sensor data is generally unavailable for most cities. A set of 180 flood-prone streets was used to train a base model, while another set of 180 flood-prone streets was used to re-train that model using different TL strategies. The results show that full-weight re-training proved most effective and minimal re-training of only the output layer was insufficient. The advantage of TL was most pronounced when target data was limited, meaning data collected at the new water depth sensor location included generally less than 18 flood events. As target data increased beyond 18 flood events, the benefit of TL diminished relative to training a ML model directly on the local flood events. These findings can assist cities as they implement street-scale flood sensing systems to create accurate forecasts for new sensing locations that do not yet have sufficient data records to train a local ML model.

Roy, Binata [Univ. of Virginia, Charlottesville, V↗
DOE OSTI2024

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tidal energy extraction modifies tidal asymmetry and transport in a shallow, well-mixed estuary

Tidal energy extraction is increasingly being studied as a potential renewable energy resource in estuaries worldwide. Although it is understood that energy extraction via tidal stream turbines can modify currents and transport within estuaries, it is not clear how the underlying nonlinear physical mechanisms dictating tidal hydrodynamics are modulated. This research investigates the influence of a hypothetical tidal stream turbine array on barotropic tidal processes in a shallow, well-mixed system: the Piscataqua River – Great Bay (PRGB) estuary, using a numerical model. The modeled turbine farm includes 180 turbines which would extract an estimated 44.7 GWh of energy, annually. The tidal hydrodynamic model for the existing condition is validated with in-situ observations of currents and water level before analyzing tidal asymmetry and transport with and without tidal turbines. Results indicate that the tidal turbine array will decrease tidal elevation and current magnitudes system-wide, but generally reduce ebb currents and transport more than flood over most of the estuary footprint, thereby diminishing tidal asymmetry. The smaller asymmetric distortion compared to the no-turbine case is attributed to reductions in the storage volume of water over the estuary’s extensive tidal flat regions between low and high waters which decreases the associated nonlinear intertidal storage mechanism up to 25%. This leads to weakened ebb dominance over estuary sections from the mouth to mid-reaches, where depths are deep enough to keep the combined nonlinear shallow water and frictional effects from asserting control over the storage mechanism. Even in upstream shallow regions where depth-dependent friction controls asymmetry in both cases, the frictional mechanism is reduced only by 10% with turbines. Some environmental considerations of this work are discussed, with focus on sediment transport, water quality, and ecology.

16 TIDAL AND WAVE POWER↗
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↗
DOE OSTI2021

A database of global storm surge reconstructions

Storm surges are among the deadliest coastal hazards and understanding how they have been affected by climate change and variability in the past is crucial to prepare for the future. However, tide gauge records are often too short to assess trends and perform robust statistical analyses. Here we use a data-driven modeling framework to simulate daily maximum surge values at 882 tide gauge locations across the globe. We use five different atmospheric reanalysis products for the storm surge reconstruction, the longest one going as far back as 1836. The data that we generate can be used, for example, for long-term trend analyses of the storm surge climate and identification of regions where changes in the intensity and/or frequency of storms surges have occurred in the past. It also provides a better basis for robust extreme value analysis, especially for tide gauges where observational records are short. The data are made available for public use through an interactive web-map as well as a public data repository.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The Pristine dwarf galaxy survey – IV. Probing the outskirts of the dwarf galaxy Boötes I

ABSTRACT We present a new spectroscopic study of the dwarf galaxy Boötes I (Boo I) with data from the Anglo-Australian Telescope and its AAOmega spectrograph together with the Two Degree Field multi-object system. We observed 36 high-probability Boo I stars selected using Gaia Early Data Release 3 proper motions and photometric metallicities from the Pristine survey. Out of those, 27 are found to be Boo I stars, resulting in an excellent success rate of 75 per cent at finding new members. Our analysis uses a new pipeline developed to estimate radial velocities and equivalent widths of the calcium triplet lines from Gaussian and Voigt line profile fits. The metallicities of 16 members are derived, including 3 extremely metal-poor stars ([Fe/H] < −3.0), which translates into a success rate of 25 per cent at finding them with the combination of Pristine and Gaia. Using the large spatial extent of our new members that spans up to 4.1 half-light radii and spectroscopy from the literature, we find a systemic velocity gradient of 0.40 ± 0.10 km s−1 arcmin−1 and a small but resolved metallicity gradient of −0.008 ± 0.003 dex arcmin−1. Finally, we show that Boo I is more elongated than previously thought with an ellipticity of ϵ = 0.68 ± 0.15. Its velocity and metallicity gradients as well as its elongation suggest that Boo I may have been affected by tides, a result supported by direct dynamical modelling.

Longeard, Nicolas↗
DOE OSTI2023

Shallow-Water Habitat in the Lower Columbia River Estuary: A Highly Altered System

Abstract Decreases in shallow-water habitat area (SWHA) in the Lower Columbia River and Estuary (LCRE) have adversely affected salmonid populations. We investigate the causes by hindcasting SWHA from 1928 to 2004, system-wide, based on daily higher high water (HHW) and system hypsometry. Physics-based regression models are used to represent HHW along the system as a function of river inflow, tides, and coastal processes, and hypsometry is used to estimate the associated SWHA. Scenario modeling is employed to attribute SWHA losses to levees, flow regulation, diversion, navigational development, and climate-induced hydrologic change, for subsidence scenarios of up to 2 m, and for 0.5 m fill. For zero subsidence, the system-wide annual-average loss of SWHA is 55 ± 5%, or 51 × 10 5 ha/year; levees have caused the largest decrease ( $${54}_{-14}^{+5}$$ 54 - 14 + 5 %, or ~ 50 × 10 5 ha/year). The loss in SWHA due to operation of the hydropower system is small, but spatially and seasonally variable. During the spring freshet critical to juvenile salmonids, the total SWHA loss was $${63}_{-3}^{+2}$$ 63 - 3 + 2 %, with the hydropower system causing losses of 5–16% (depending on subsidence). Climate change and navigation have caused SWHA losses of $${5}_{-5}^{+16}$$ 5 - 5 + 16 % and $${4}_{-6}^{+14}$$ 4 - 6 + 14 %, respectively, but with high spatial variability; irrigation impacts have been small. Uncertain subsidence causes most of the uncertainty in estimates; the sum of the individual factors exceeds the total loss, because factors interact. Any factor that reduces mean or peak flows (reservoirs, diversion, and climate change) or alters tides and along-channel slope (navigation) becomes more impactful as assumed historical elevations are increased to account for subsidence, while levees matter less.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Alternating Conditional Expectations: Introducing a Non‐Parametric Statistical Method to Interpret Long‐Term Greenhouse Gas Flux Measurements Over Semi‐Arid and Wetland Ecosystems

Abstract We explore the potential of using a non‐parametric statistical method called Alternating Conditional Expectations, ACE, to quantify functional relationships in biogeosciences. Here, ACE is used to quantify the non‐linear and multi‐faceted responses of greenhouse gas fluxes to a set of biophysical forcings, when the shapes of those response surfaces are unknown. We evaluated the statistical method over two contrasting ecosystems and two contrasting time steps. One case involved quantifying the biophysical controls of water vapor and carbon dioxide (CO 2 ) fluxes over a semi‐arid oak savanna using daily integrated fluxes. The other case evaluated the responses of CO 2 and methane (CH 4 ) flux measurements to a set of biophysical forcings at a restored tidal wetland using thirty‐minute averages. The statistical model, based on 4 independent variables, explained up over 90% of the variation in daily integrated flux densities of water vapor and net carbon dioxide exchange at the savanna site. This fit was defined by distinct non‐linear responses to such drivers as gross primary production, photosynthetically active radiation, air temperature, vapor pressure deficit and soil moisture. At the tidal wetland site, we evaluated net carbon dioxide and methane fluxes with short‐term measurements to capture the influence of rising and falling tides and seasonality in biological activity. The statistical model defined the shape of the forcing of fluxes due to the roles of carbon exudates, water table depth, oxygen level in the water column, temperature and vegetation status. The statistical fits of the greenhouse gas fluxes were less precise than the savanna case. The fetch varies on a run‐to‐run basis as it is comprised of a heterogeneous mosaic of open water and vegetation. Furthermore, it is difficult to monitor the environmental conditions of the archaea and bacteria in the sediments that produce methane and carbon dioxide.

Environmental Sciences & Ecology↗
DOE OSTI2025

Forecasting Multi-Step-Ahead Street-Scale Nuisance Flooding using a seq2seq LSTM Surrogate Model for Real-Time Application in a Coastal-Urban City

In coastal-urban cities facing an elevated risk of nuisance flooding (by rain and tide) due to increased heavy rainfall, sea level rise, urbanization, and aging drainage systems, real-time flood forecasting at the street-scale can provide useful information to transportation decision-makers. Physics-Based Models (PBMs) that offer high accuracy come with high computational runtimes and costs that limit their application for real-time flood forecasting. To address this challenge, Machine Learning (ML) surrogate models trained from PBMs have been proposed to provide street-scale flood forecasts. Previous related studies have focused on using Long Short-Term Memory (LSTM) architectures to model hourly flood depth on streets. While LSTM models can capture input sequences effectively, they fall short in accurately preserving output sequences, limiting their suitability for multi-step-ahead forecasts. The seq2seq LSTM architecture offers a key advantage here by capturing the full sequence of input–output, making it potentially more suitable for multi-step-ahead flood forecasts compared to traditional LSTM models. However, seq2seq LSTM has not been tested for street-scale flood forecasting, particularly for rapidly fluctuating nuisance flooding events which require special attention to its temporal sequences. Hence, in this study, we applied the seq2seq LSTM model to explore multi-step-ahead street-scale nuisance flooding and compared its results to the traditional LSTM model as a benchmark model. LSTM and seq2seq LSTM surrogate models were applied to 22 flood-prone streets in Norfolk, Virginia, as a case study with a 4-hr (short-term) and 8-hr (long-term) lead time. The models were trained with environmental (rainfall and tide) and topographic (elevation, Topographic Wetness Index, and Depth-To-Water) features along with PBM-derived water depths for different storm events. The results demonstrated satisfactory performance of both LSTM and seq2seq LSTM surrogate models throughout the forecast period compared to the PBM. However, the seq2seq LSTM showed lower Mean Absolute Error (MAE)/ Root Mean Square Error (RMSE) and higher Nash–Sutcliffe Efficiency (NSE)/ correlation than the LSTM across most lead times, particularly for long-term forecasting due to its supremacy in handling both input–output sequences together, which is missing in the traditional LSTM. For example, in the long-term, the average RMSE ranges were 0.0268–0.0373 m for LSTM and 0.0226–0.0319 m for seq2seq LSTM, while in the short-term, they were 0.0263–0.0293 m and 0.0261–0.0283 m, respectively. Additionally, while both models exhibited similar performance in distinguishing flooded and non-flooded streets for flood depth ≥ 0.1 m, the seq2seq LSTM model demonstrated superior performance for higher flood depths (such as ≥ 0.2 m and ≥ 0.3 m). Once trained, inference took only 0.09 to 0.11 s (short-term) and 0.30 to 0.35 s (long-term) per storm event for the 22 streets, making the application highly suitable for real-time decision-making during nuisance flood events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The dynamical tides of spinning Newtonian stars

ABSTRACT We carefully develop the framework required to model the dynamical tidal response of a spinning neutron star in an inspiralling binary system, in the context of Newtonian gravity, making sure to include all relevant details and connections to the existing literature. The tidal perturbation is decomposed in terms of the normal oscillation modes, used to derive an expression for the effective Love number which is valid for any rotation rate. In contrast to previous work on the problem, our analysis highlights subtle issues relating to the orthogonality condition required for the mode-sum representation of the dynamical tide and shows how the prograde and retrograde modes combine to provide the overall tidal response. Utilizing a slow-rotation expansion, we show that the dynamical tide (the effective Love number) is corrected at first order in rotation, whereas in the case of the static tide (the static Love number) the rotational corrections do not enter until second order.

Pnigouras, P. (ORCID:0000000318959431)↗
DOE OSTI2022

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗
DOE OSTI2024

Mixing of the Connecticut River Plume During Ambient Flood Tides: Spatial Heterogeneity and Contributions of Bottom–Generated and Interfacial Mixing

The Connecticut River plume is influenced by energetic ambient tides in the Long Island Sound receiving waters. The objectives of this modeling study are (a) characterizing the spatial heterogeneity of turbulent buoyancy fluxes, (b) partitioning turbulent buoyancy fluxes into bottom-generated and interfacial shear contributions, and (c) quantifying contributions to plume-integrated mixing within the tidal plume. The plume formed during ambient flood tides under low river discharge, spring tides, and no winds is analyzed. Turbulent buoyancy fluxes (B) and depth-integrated B through the plume (Bd) are characterized by pronounced spatial heterogeneity. Strong mixing (Bd ~ 10 –5 -10 –4 m 3 /s 3 ) occurs near the mouth, in the nearfield plume turning region, over shoals, and nearshore shallow areas. Low to moderate mixing (Bd ~ 10 –8 -10 –6 m 3 /s 3 ) occupies half the plume. Buoyancy fluxes are first partitioned based on the depth of the shear stress minimum between plume-generated and bottom-generated shear maxima. Four other tested partitioning methods are based on open channel flow and stratified shear flow parameterizations. Interfacial and bottom-generated shear contribute to different areas of intense and moderate mixing. All methods indicate a significant plume mixing role for bottom-generated mixing, but interfacial mixing is a bigger contributor. Plume-integrated total and interfacial mixing peak at max ambient flood and the timing of peak bottom-generated mixing varies among partitioning methods. Two-thirds of the mixing occurs in concentrated intense mixing areas. A parameter space with the ambient tidal Froude number and plume thickness to depth ratio as axes indicates many tidally modulated plumes are moderately to dominantly influenced by bottom-generated tidal mixing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Comparative Study of Physics‐Informed and Data‐Driven Neural Networks for Compound Flood Simulation at River‐Ocean Interfaces: A Case Study of Hurricane Irene

Simulating compound flooding (CF) at the river-ocean interface within large-scale Earth System Models (ESMs) presents significant challenges due to complex interactions between river discharge, storm surge, and tides. This study assesses the comparative advantages of physics-informed and data-driven machine learning (ML) approaches for enhancing local ESM performance. We systematically compare data-driven neural network models (i.e., CNNs, U-Net, Long Short-Term Memory (LSTM), Gated Recurrent Unit), and physics-informed neural network (PINN) models, including vanilla PINN and a finite-difference-based PINN (FD-PINN). Specifically, FD-PINN is introduced to enhance computational efficiency, accelerating vanilla PINNs by ∼6.5 times while improving accuracy. To enhance data-driven model training, a new data-generation approach is developed to sample historical fluvial and coastal flood events, which ensures a robust data set for extreme event prediction. The models are evaluated using a realistic one-dimensional river domain extracted from an ESM's river mesh and the Hurricane Irene event as an independent test case. Results show that FD-PINN achieves accurate predictions with significantly reduced computational costs relative to vanilla PINNs. Among data-driven models, the best overall performance is achieved by a CNN-LSTM hybrid, which balances accuracy and efficiency. While a fully connected CNN (CNN-FC) provides the best accuracy, it incurs high computational cost. Architectures lacking strong temporal modeling tend to underperform on unseen events. These findings highlight the importance of sequence-aware designs for robust generalization. This study reveals the trade-offs between physics-informed and data-driven models and proposes an adaptive hybrid framework for integrating ML into ESMs to enhance local flood simulations.

Earth Systems Modeling↗
DOE OSTI2024

Simulation of Compound Flooding Using River‐Ocean Two‐Way Coupled E3SM Ensemble on Variable‐Resolution Meshes

Abstract Coastal zone compound flooding (CF) can be caused by the interactive fluvial and oceanic processes, particularly when coastal backwater propagates upstream and interacts with high river discharge. The modeling of CF is limited in existing Earth System Models (ESMs) due to coarse mesh resolutions and one‐way coupled river‐ocean components. In this study, we present a novel multi‐scale coupling framework within the Energy Exascale Earth System Model (E3SM), integrating global atmosphere and land with interactively coupled river and ocean models using different meshes with refined resolutions near the coastline. To evaluate this framework, we conducted ensemble simulations of a CF event (Hurricane Irene in 2011) in a Mid‐Atlantic estuary. The results demonstrate that the novel E3SM configuration can reasonably reproduce river discharge and sea surface height variations. The two‐way river‐ocean coupling improves the representation of coastal backwater effects at the terrestrial‐aquatic interface that are caused by the combined actions of tide and storm surge during the CF event, thus providing a valuable modeling tool for better understanding the river‐estuary‐ocean dynamics in extreme events under climate change. Notably, our results show that the most significant CF impacts occur when the highest storm surge generated by a tropical cyclone meets with a moderate river discharge. This study highlights the state‐of‐the‐art advancements developed within E3SM for simulating multi‐scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗