Search NASA⌕ Search

SEARCH · Search NASA

Results for “shallow water”

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

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantification of Off-Channel Inundated Habitat for Pacific Chinook Salmon ( Oncorhynchus tshawytscha ) along the Sacramento River, California, Using Remote Sensing Imagery

Off-channel areas are one of the most impacted aquatic habitats by humans globally, as extensive agricultural and urban development has limited them to roughly 10% of historical extent. This is also true for California’s Sacramento River Valley, where historically frequent widespread inundation has been reduced to a few off-channel water bodies along the mid-Sacramento River. This remaining shallow-water habitat provides crucial ecological benefits to multiple avian and fish species, but especially to floodplain-adapted species such as Chinook salmon (Oncorhynchus tshawytscha). Characterizing spatiotemporal off-channel dynamics, including inundation extent and residence time, is fundamental to better understanding the intrinsic value of such habitats and their potential to support recovery actions. Remote sensing techniques have been increasingly used to map surface water at regional and local scales, with improved resolutions. As such, this study maps off-channel inundation areas and describes their temporal dynamics by analyzing pixel-based time- series of multiple water indices, modified Normalized Difference Water Index (mNDWI) and the Automated Water Extraction Index (AWEI), generated from LandSat-8 and Sentinel-2 data between 2013–2021. Quantified off-channel area was similar with each water index and method used, but improved performance was associated with Sentinel-2 products and AWEI index to identify wetted areas under lower mainstem discharges. Results indicate an uneven distribution of off-channel habitat in the study area, with limited inundated areas in upstream reaches (<16% of total off-channel area for greater flows). In addition, much less habitat exists for flows under 400 m 3 /s, an important migration cue for endangered winter-run Chinook salmon, limiting juvenile access to areas with enhanced rearing conditions. Off-channel habitat residence times averaged between 7 and 16 days, primarily defined by the rate of receding flows, with rapid flow recession providing marginal off-channel habitat. This study shows reasonable performance of moderate resolution LandSat-8 and Sentinel-2 remote sensing imagery to characterize shallow-water inundated habitat in higher-order rivers, and as a method to inform restoration and native fish recovery efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multiple biomarkers highlight the importance of water column processes in treatment wetland organic matter cycling

A suite of biomarkers, including amino acids, pigments, and lignin phenols coupled with high resolution mass spectrometry were used to evaluate differences in the sources and fate of organic matter (OM) in Everglades treatment wetlands as a model for OM cycling in shallow water wetlands. Five components of the system (water column particulate organic matter, vertical traps, flocculent material, periphyton, and surface soil) were assessed for OM transformations down-profile (i.e. water column to soil) and between treatment cells dominated by emergent aquatic vegetation (EAV) and submerged aquatic vegetation (SAV), with comparisons to reference sites that resemble historic Everglades conditions. We found that OM cycling is fundamentally different between EAV and SAV wetlands, and that SAV wetlands have some shared characteristics with similar habitats in the Everglades. Other than locations densely populated by Typha spp., water column particulate organic C was derived from microbial/cryptomonad sources, rather than macroscopic sources (vascular plants and algal mats). Bacterial amino acid biomarkers were positively correlated with amino acid degradation indices and organic P (Po), respectively suggesting that microbial abundance is associated with "fresher" OM, and that further investigation into relationships between microbial biomass and Po is warranted. Overall, this multi-biomarker approach can evaluate the relative "freshness" of OM pools, identify sources of OM, and highlight the importance of water column processes in shallow water wetlands.

treatment wetlands, organic matter, Everglades, am↗
DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗
DOE OSTI2024

Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai

The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Light Regulation of Phytoplankton Growth in San Francisco Bay Studied Using a 3D Sediment Transport Model

In San Francisco Bay (SFB), light availability is largely determined by the concentration of suspended particulate matter (SPM) in the water column. SPM exhibits substantial variation with time, depth, and location. To study how SPM influences light and phytoplankton growth, we coupled a sediment transport model with a hydrodynamic model and a biogeochemical model. The coupled models were used to simulate conditions for the year of 2011 with a focus on northern SFB. For comparison, two simulations were conducted with ecosystem processes driven by SPM concentrations supplied by the sediment transport model and by applying a constant SPM concentration of 20 mg l –1 . The sediment transport model successfully reproduced the general pattern of SPM variation in northern SFB, which improved the chlorophyll-a simulation resulting from the biogeochemical model, with vertically integrated primary productivity varying greatly, from 40 g[C] m –2 year –1 over shoals to 160 g[C] m –2 year –1 in the deep channel. Primary productivity in northern SFB is influenced by euphotic zone depth ( Ze ). Our results show that Ze in shallow water regions (<2 m) is mainly determined by water depth, while Ze in deep water regions is controlled by SPM concentration. As a result, Ze has low (high) values in shallow (deep) water regions. Large (small) differences in primary productivity exist between the two simulations in deep (shallow) water regions. Furthermore, we defined a new parameter F light for “averaged light limitation” in the euphotic zone. The averaged chlorophyll-a concentration in the euphotic zone and F light share a similar distribution such that both have high (low) values in shallow (deep) water regions. Our study demonstrates that light is a critical factor in regulating the phytoplankton growth in northern SFB, and a sediment transport model improves simulation of light availability in the water column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Large isotopic variability at the micron-scale in ‘Shuram’ excursion carbonates from South Australia

Ediacaran-aged (635–541 million years ago) marine sediments contain a large negative carbon isotope (δ 13 C) excursion, in which carbonate δ 13 C values reach -12‰ (VPDB). Known as the ‘Shuram’ excursion, many workers have interpreted this δ 13 C record as an unprecedented perturbation to the global carbon cycle, leading to speculation about a causal connection to the broadly contemporaneous rise of animal life. Others have interpreted the δ 13 C signal as a product of diagenesis, thereby minimizing its relevance for understanding the evolution of metazoans. In this work, we present SEM imaging and in-situ δ 13 C and δ 18 C values measured by secondary ion mass spectrometry (SIMS) to assess these competing hypotheses in the Wonoka Formation of South Australia. Our results from the minimum of the excursion show that rounded sedimentary grains of calcite have δ 13 C values between -12.8 to -10.6‰ and δ 18 C values between 17.8 to 15.5‰ (VPDB). Euhedral dolomite that appears to have grown unimpeded in open sedimentary pore spaces also is present. These early-stage dolomites are interpreted as early authigenic in origin and have δ 13 C values that reach +5‰, requiring a formation fluid with a substantially different δ 13 C composition from basin waters or bulk sediment. Together, these results provide little evidence for the hypothesis that a late diagenetic overprint has generated the ‘Shuram’ excursion in the Wonoka. Instead, they suggest the presence of a large carbon isotopic gradient in the surface environment, with shallow waters capable of precipitating carbonates with very low δ 13 C (down to 12‰) and deeper shelf and/or marine pore waters generating carbonates with positive carbon isotope values (up to +5‰). Because negative isotope excursions of similar magnitude are found in widely dispersed Ediacaran basins, it is likely that this gradient was characteristic of shelf environments of this period and that a still-unknown global process led to the ‘Shuram’ excursion in shallow water carbonates.

58 GEOSCIENCES↗
DOE OSTI2026

Amapari Marker Band Metal‐Enrichments: Potential Mechanisms and Implications for Surface and Subsurface Water and Weathering in Gale Crater

NASA's Curiosity rover is exploring a 5 km tall sedimentary mound that is hypothesized to record the transition from a warm and wet (phyllosilicate-rich) to a cold and drier (sulfate-rich) Mars. Evidence of magnesium sulfate-bearing rock has shown that Curiosity has crossed through this phyllosilicate-sulfate transition. Recently, Curiosity arrived at the Amapari Marker Band, a darker, indurated unit that can be traced laterally for tens of kilometers in orbiter images. Here, Curiosity found evidence for a very broad lake, and bedforms interpreted as wave-ripple laminated sedimentary rock that likely was deposited in shallow water in the explored location, before becoming a deeper lake. These rocks are enriched in Fe, Mn, and Zn which has major implications for groundwater paleohydrology in Gale crater. Three formation hypotheses are considered: concretion formation during early diagenetic alteration of shallow lake sediments, laterization or leaching of the sediments, and addition of Fe, Mn, and Zn by a mildly acidic and reducing groundwater interacting with a redox and/or pH front in a stratified lake. The preferred interpretation of the metal enrichments within the Amapari Marker band sedimentary rocks is that they formed in a shallow water environment at a redox and/or pH front within the ripple unit, which drove precipitation and concentration of metals. If the enrichments are due to groundwater alteration, these processes could link subsurface and surface environments. Water and the presence of high amounts of redox sensitive elements and other metals are favorable indicators for habitability.

58 GEOSCIENCES↗
DOE OSTI2020

Performance Evaluation of a Two-Dimensional Flood Model on Heterogeneous High-Performance Computing Architectures

This paper describes the implementation of a two-dimensional hydrodynamic flood model with two different numerical schemes on heterogeneous high-performance computing architectures. Both schemes were able to solve the nonlinear hyperbolic shallow water equations using an explicit upwind first-order approach on finite differences and finite volumes, respectively, and were conducted using MPI and CUDA. Four different test cases were simulated on the Summit supercomputer at Oak Ridge National Laboratory. Both numerical schemes scaled up to 128 nodes (768 GPUs) with a maximum 98.2x speedup of over 1 GPU. The lowest run time for the 10 day Hurricane Harvey event simulation at 5 meter resolution (272 million grid cells) was 50 minutes. GPUDirect communication proved to be more convenient than the standard communication strategy. Both strong and weak scaling are shown.

Sharif, Md Bulbul↗
DOE OSTI2020

Umbilical Fatigue Analysis for a Wave Energy Converter

Unique umbilical designs for wave energy converters (WECs), including the ability to handle significantly larger motions and loads over long deployments, are often required when conventional marine umbilical designs for offshore oil and gas and offshore wind may not meet the design and cost needs of wave energy technologies. This study details a fatigue analysis of a dynamic power umbilical attached to a two-body floating point absorber WEC system, using the sea states provided for the PacWave testing facilities. The 6 degrees of freedom motion time history for the WEC was simulated, and the motions of the attachment point for the umbilical on the WEC and respective sea states were used to analyze the dynamic motions and fatigue of the connected power umbilical to predict the fatigue life. The results show that the fatigue damage observed is more significant in shallow water, and extensive fatigue damage may occur because of the larger curvature response of the umbilical. The umbilical configurations departing at 90 deg off incoming waves were found to have the highest fatigue life attributed to less extension or compression of the umbilical. However, additional bend stiffener/limiter features may need to be incorporated into the buoyancy section and touchdown regions to minimize curvature-induced fatigue.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2021

On steady alternate bars forced by a localized asymmetric drag distribution in erodible channels

Studying the effect of different in-stream fluvial turbines siting on river morphodynamics allowed us to witness the onset of a time-averaged, large-scale, alternate distortion of bed elevations, which could not be exclusively related to the turbine rotor blockage. The longitudinal profiles of this two-dimensional bathymetric perturbation resemble those of steady fluvial bars. In this contribution we generalize the problem addressing a spatially impulsive, asymmetric distribution of drag force in the channel cross-section. Furthermore, this is experimentally investigated through the deployment of differently sized grids perpendicular to the flow, and analytically explored as a finite perturbation of an open channel flow over an erodible sediment layer, as described by a coupled flow–sediment shallow water equation. The steady solutions of this fluvial morphodynamic problem, physically represented by alternate bars scaling with the channel width, highlight the importance of the resonant conditions in defining the spatial extent of the bed deformation. The equations further suggest that in very shallow flows any asymmetric obstruction may lead to an upstream propagation of the steady bars, consistent with previous studies on the effects of channel curvature. In broad terms, this study provides the preliminary framework to control the onset of river meandering through imposed finite perturbations of the cross-section. In a more applied sense, it provides a tool to predict non-local scour–deposition patterns associated with the deployment of energy converters or other flow obstructions.

42 ENGINEERING↗
DOE OSTI2026

Quantum mechanical closure of partial differential equations with symmetries

We develop a statistical framework for the dynamical closure of spatiotemporal dynamics governed by partial differential equations. Employing the mathematical framework of quantum mechanics to embed the original classical dynamics into a quantum mechanical representation, we use the space of quantum density operators to model the unresolved degrees of freedom of the original dynamics in a statistical sense, and the framework of quantum measurement to predict their contributions to the resolved dynamics. The embedded dynamics is discretized by a positivity preserving process, leading to a compressed representation that is invariant under the dynamical symmetries of the resolved dynamics. We present a data based formulation of the closure scheme and apply it to a closure problem for the shallow water equations. The numerical results demonstrate that our closure model can accurately predict the main features of the true dynamics, including for out of sample initial conditions.

Delay embedding↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Analysis of the performance of different culvert boundary conditions in 2D shallow flow models

Culverts allow roads to safely traverse small streams or drainage ditches, and their proper design is critical to ensure a safe and reliable transportation network. A correct modelization of these hydraulic structures becomes crucial in the assessment of flood footprints or discharge peak estimation in a risk evaluation plan. The question of how to include culverts comes up frequently when assembling a hydraulic model that requires the presence of as many singular elements as possible. In this work, three different culvert integrations with the surface domain are studied and compared in the context of a 2D shallow water (SW) model. All of them are based on the Federal Highway Administration (FHWA) formulation for the culvert discharge estimation but differ in complexity and in the interaction with the numerical model for surface flow, some of them as internal boundary conditions. Several steady and unsteady validation test cases are presented and the numerical results are compared with the predictions from HEC-RAS 1D and HY-8 software. Furthermore, the culvert area, shape and their sensitivity to the 2D computational mesh is also analyzed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Parallel exponential time differencing methods for geophysical flow simulations

Two ocean models are considered for geophysical flow simulations: the multilayer shallow water equations and the multilayer primitive equations. For the former, we investigate the parallel performance of exponential time differencing (ETD) methods, including exponential Rosenbrock–Euler, ETD2wave, and B-ETD2wave. For the latter, we take advantage of the splitting of barotropic and baroclinic modes and propose a new two-level method in which an ETD method is applied to solve the fast barotropic mode. Furthermore, these methods could improve the computational efficiency of numerical simulations because ETD methods allow for much larger time step sizes than traditional explicit time-stepping techniques that are commonly used in existing computational ocean models. Several standard benchmark tests for ocean modeling are performed and comparison of the numerical results demonstrates a great potential of applying the parallel ETD methods for simulating real-world geophysical flows.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A conservative discontinuous-Galerkin-in-time (DGiT) multirate time integration framework for interface-coupled problems with applications to solid–solid interaction and air–sea models

In this paper we extend the DGiT multirate framework, developed in Connors and Sockwell (2022) for scalar transmission problems, to a solid–solid interaction (SSI) problem involving two coupled elastic solids and a coupled air–sea model with the rotating, thermal shallow water equations. In so doing we aim to demonstrate the broad applicability of the mathematical theory and governing principles established in Connors and Sockwell (2022) to coupled problems characterized by subproblems evolving at different temporal scales. Further, multirate time integration algorithms employing different time steps, optimized for the dynamics of each subproblem, can significantly improve simulation efficiency for such coupled problems. However, development of multirate algorithms is a highly non-trivial task due to the coupling, which can impact accuracy, stability or other desired properties such as preservation of system invariants. DGiT provides a general template for multirate time integration that can achieve these properties. To elucidate the manner in which DGiT accomplishes this task, we fully detail each step in the application of the framework to the SSI and air–sea coupled problems. Numerical examples illustrate key properties of the resulting multirate schemes for both problems.

42 ENGINEERING↗
DOE OSTI2021

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗