SEARCH · Engineering Papers

Results for “topography”

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Timeseries Photos of a Variably Inundated Stream: Umtanum Creek, Washington, United States

This dataset is associated with a broader study using game camera timeseries photos collected to evaluate stream variable inundation via changes in width (i.e. wet fraction). Four game cameras were deployed along Umtanum Creek (Washington, United States) to track changes in stream inundation over time. Drone imagery was collected at the same location on October 18, 2024 which was used to construct a digital elevation model (DEM) of the streambed topography. The associated paper and data can be found at https://doi.org/10.1016/j.envsoft.2025.106715 (Bao et al., 2025a)) and https://doi.org/10.15485/2589885 (Bao et al., 2025b), respectively. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to this readme, this data package also includes a file-level metadata (FLMD) files that describes each file and a data dictionaries (DD) that describe all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) readme; (4) field metadata; (5) field protocol; and (5) folders containing game camera photos. Game camera photos are organized into folders for each camera (CDL, CUL, CDR, CUR; see readme for information on camera naming) by the month photos were collected. All files are .csv, .jpg, or .pdf.

AI image segmentation↗

DOE OSTI2024

Growth, mortality, wood density, biomass data from BIONTE inventories in Manaus, Brazil

BIONTE (BIOmass and NuTrient Experiment) is a selective logging experiment established at the Experimental Station of Tropical Forestry (EEST, aka “ZF2”) field research station in the mid 1980s in the central Amazon (Higuchi et al. 1997, Amaral et al. 2019). Led by the National Institute for Amazon Research (INPA) in Brazil, the project aimed at assessing the effects of logging intensity on forest dynamics and enabling the creation of a model of forest management for the Central Amazon. The experiment included three levels of increasing selective logging intensity and controls, with 1 hectare sample plots (12 total) located at the center of 4 hectare treatment plots. The site’s Köppen classification is tropical rainforest (Af), characterized by high temperatures and humidity, with mean annual temperatures around 27 ℃ and mean annual precipitation around 2200 mm of rain. The vegetation has a high floristic diversity, the soils of the region are poor in nutrients, and the topography is characterized by plateaus (where BIONTE is located), and also valley bottoms and slopes. The inventory (growth and mortality) and biomass data included here covers the 1990 to 2019 period, with wood density being averaged from existing datasets. This dataset includes a data file in .csv file format and a .txt file, BIONTE_mortality-rates_headers.txt, that provides descriptions for the data file headers.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

AmeriFlux US-xNW NEON Niwot Ridge Mountain Research Station (NIWO)

This is the AmeriFlux version of the carbon flux data for the site US-xNW NEON Niwot Ridge Mountain Research Station (NIWO). Site Description - The Niwot Ridge sits approximately 27 km west of Boulder, Colorado, and 6 km east of the Continental Divide. Topography, climate, and biota of the site are representative of Rocky Mountain alpine ecosystems, including extensive alpine tundra (mostly herbs, some shrubs and scree) and subalpine coniferous forests (Abies lasciocarpa and Picea engelmanii at higher elevations), talus slopes, wetlands and a variety of glacial landforms. Characterized by cold and relatively long winters, Niwot Ridge has an average annual temperature of 1.5°C and average annual precipitation of 800 mm. Most precipitation falls as snow and summer precipitation falls primarily during afternoon thunderstorms. Located on the eastern side of the Continental Divide at 3,000-3,500 m elevation, the site best captures chemical inputs produced along the Front Range and is well situated to observe other east/west flows across the Southern Rockies in conjunction with other NEON sites.

Network), NEON (National Ecological Observatory↗

DOE OSTI2022

AmeriFlux FLUXNET-1F CA-Gro Ontario - Groundhog River, Boreal Mixedwood Forest

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site CA-Gro Ontario - Groundhog River, Boreal Mixedwood Forest. This is the FLUXNET version of the carbon flux data for the site CA-Gro Ontario - Groundhog River, Boreal Mixedwood Forest produced by applying the standard ONEFlux (1F) software. Site Description - Groundhog River (FCRN or CCP site "ON-OMW") is situated in a typical boreal mixedwood forest in northeastern Ontario (48.217 degrees north and 82.156 degrees west) about 80 km southwest of Timmins in Reeves Twp. near the Groundhog River. Rowe (1972) places the site in the Missinaibi-Cabonga Section of the Boreal Forest Region. In terms of ecoregion and ecozone, the site is in the Lake Timiskaming Lowlands of the Boreal Shield. The forest developed after high-grade logging in the 1930's. The average age in 2013 is estimated at beteen 75 and 80 years. The forest is dominated by five species characteristic of Ontario boreal mixedwoods: trembling aspen (Populus tremuloides Michx.), black spruce (Picea mariana (Mill.) B.S.P.), white spruce (Picea glauca (Moench.) Voss.), white birch (Betula papyrifera Marsh.), and balsam fir (Abies balsamea (L.) Mill.). The surficial geology is a lacustrine deposit of varved or massive clays, silts and silty sands. The soil is an orthic gleysol with a soil moisture regime classified as fresh to very fresh. Plonski (1974) rates it as a site class 1. The topography is simple and flat with an overall elevation of 340 m ASL.

McCaughey, Harry↗

DOE OSTI2025

AmeriFlux FLUXNET-1F US-xNW NEON Niwot Ridge Mountain Research Station (NIWO)

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-xNW NEON Niwot Ridge Mountain Research Station (NIWO). This is the FLUXNET version of the carbon flux data for the site US-xNW NEON Niwot Ridge Mountain Research Station (NIWO) produced by applying the standard ONEFlux (1F) software. Site Description - The Niwot Ridge sits approximately 27 km west of Boulder, Colorado, and 6 km east of the Continental Divide. Topography, climate, and biota of the site are representative of Rocky Mountain alpine ecosystems, including extensive alpine tundra (mostly herbs, some shrubs and scree) and subalpine coniferous forests (Abies lasciocarpa and Picea engelmanii at higher elevations), talus slopes, wetlands and a variety of glacial landforms. Characterized by cold and relatively long winters, Niwot Ridge has an average annual temperature of 1.5°C and average annual precipitation of 800 mm. Most precipitation falls as snow and summer precipitation falls primarily during afternoon thunderstorms. Located on the eastern side of the Continental Divide at 3,000-3,500 m elevation, the site best captures chemical inputs produced along the Front Range and is well situated to observe other east/west flows across the Southern Rockies in conjunction with other NEON sites.

Network), NEON (National Ecological Observatory [N↗

DOE OSTI2026

AmeriFlux FLUXNET-1F US-RC4 Moscow Mountain on-farm site

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-RC4 Moscow Mountain on-farm site. This is the FLUXNET version of the carbon flux data for the site US-RC4 Moscow Mountain on-farm site produced by applying the standard ONEFlux (1F) software. Site Description - The Moscow Mountain On-farm site operated 2012-2016 as part of a cluster of 5 towers (RC1 to RC5) operated for the Regional Approaches to Climate Change (REACCH) USDA-supported research project. The site is located in the high precipitation agroecological zone of the Columbia Plateau’s dryland cropping region, near Moscow, Idaho. The site was managed with conventional (reduced) tillage. The crop rotation was spring barley, spring field peas, winter wheat, and spring wheat. Soils are silt loam mollisols including the Latahco, Thatuna, Southwick, and Larkin soil series. The site topography is sloping and is situated within the rolling hills of the Palouse region.

Chi, Jinshu [The Hong Kong University of Science a↗

DOE OSTI2026

AmeriFlux FLUXNET-1F US-RC5 Moses Lake on-farm site

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-RC5 Moses Lake on-farm site. This is the FLUXNET version of the carbon flux data for the site US-RC5 Moses Lake on-farm site produced by applying the standard ONEFlux (1F) software. Site Description - The Moses Lake On-farm site operated 2013-2015 as part of a cluster of 5 towers (RC1 to RC5) operated for the Regional Approaches to Climate Change (REACCH) USDA-supported research project. The field was a half-circle irrigated plot with the irrigation pivot located on the edge of the field and the flux tower located next to the center of the pivot. The field was in wheat from tower establishment in June 2013 to harvest in August 2013. A cover crop of arugula and mustard was grown from August to October 2013. Potatoes were grown April to August 2014 and wheat (a spring cultivar planted in fall) was grown October 2014 to June 2015. Soils are coarse sandy loam mollisols in the Timmerman series. The site topography is flat.

Chi, Jinshu [The Hong Kong University of Science a↗

DOE OSTI2026

AmeriFlux FLUXNET-1F US-RC3 WSU Lind Dryland Research Station

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-RC3 WSU Lind Dryland Research Station. This is the FLUXNET version of the carbon flux data for the site US-RC3 WSU Lind Dryland Research Station produced by applying the standard ONEFlux (1F) software. Site Description - RC3 operated from 2011-2016 at the Washington State University Lind Dryland Experimental Station, as part of a cluster of 5 towers (RC1 to RC5) operated for the Regional Approaches to Climate Change (REACCH) USDA-supported research project. Lind is in the low precipitation region of the Columbia Plateau dryland cropping region. The area is characterized as a stable crop-fallow agroecological zone, with crops typically grown every second year to allow soil water to accumulate during intervening fallow years. Winter wheat was grown August 2012-August 2013, and August 2014-July 2015. The growing seasons of 2012 and 2014 were fallow years. Soils at Lind are predominantly silt loam texture Mollisols in the Shano and Ritzville series. Site topography is flat.

Chi, Jinshu [The Hong Kong University of Science a↗

DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

DOE OSTI2023

SKS Polarization Anomalies Due to the Coriolis Force

The Earth’s Coriolis force has been well–known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear–wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear–wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well–isolated phases. Here, we recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

58 GEOSCIENCES↗

DOE OSTI2024

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗

DOE OSTI2024

Impact of the Earthquake Rupture on Ground-Motion Variability of the 24 August 2016 M w 6.2 Amatrice, Italy, Earthquake

Here, the devastating 24 August 2016 M w 6.2 earthquake that struck Amatrice, Italy, marked the beginning of a prolonged seismic sequence dominated by three subsequent M w ≥6.0 events in the central Apennines region. The earthquake destroyed Amatrice’s historic center, claiming the lives of 299 individuals and causing widespread damage in the neighboring villages. The severity of the ground shaking, with a recorded maximum acceleration of 850 cm/s 2 on the east–west component at the Amatrice station, was far greater than the predicted acceleration based on the Italian ground-motion model (GMM). As pointed out by several investigations, the observed ground-motion amplitude and its spatial variability during the earthquake can be linked to specific rupture characteristics, including slip distribution and rupture directivity effects revealed by the observed data (Tinti et al., 2016; Pischiutta et al., 2021). In this study, we conducted physics-based 3D numerical simulations of ground motion for the Amatrice earthquake for frequencies up to 3 Hz. We employed a series of kinematic rupture models and a well-constrained local 3D velocity model incorporating surface topography. The kinematic rupture realizations were generated using multiscale hybrid and fully stochastic models, following the technique proposed by Graves and Pitarka (2016). We focused on assessing the sensitivity of near-fault ground-motion amplitudes to earthquake rupture characteristics, in particular, the spatial slip pattern. To evaluate the quality of our simulations, we employed goodness-of-fit measurements performed in comparisons of simulated and recorded ground motions. The simulated ground motions compare well with the recorded data and predictions from GMMs for Italy, ITA18 (Lanzano et al., 2019). However, we found that the simulated interevent ground-motion variability (randomness in the source process) of peak ground velocity, σ (PGV) is higher than the constant σ (PGV) predicted by conventional GMMs. Our simulations using several rupture scenarios demonstrate that the near-fault ground-motion amplification pattern is directly related to the slip distribution pattern.

3D Ground Motions Simulations↗

DOE OSTI2025

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗

DOE OSTI2022

Finite-Difference Simulation for Infrasound Generated by Finite-Extent Ground Motions

Underground explosions can produce infrasound in the atmosphere, and the wavefield characteristics are often governed by the ground surface motions. Finite-difference methods are popular for infrasound simulation as their generality and robustness allow for complex atmospheric structures and surface topography. A simple point-source approximation is often used because infrasound wavelengths tend to be large relative to the source dimensions. However, this assumption may not be able to capture the complexity of explosion-induced ground motions if the surface area is not compact, and appropriate source models must be incorporated into the finite-difference simulations for accurate infrasound prediction. In this study, we develop a point source representation of the complex ground motions for infrasound sources. Instead of a single point source, we use a series of point sources distributed over the source area. These distributed point sources can be equivalent to air volume changes produced by the ground motions in the atmosphere. Further, we apply the distributed point-source method to a series of buried chemical explosions conducted during the Source Physics Experiment Phase I. Epicentral ground-motion measurements during the experiments provide a way to calculate accurate distributed point sources. We validate and evaluate the accuracy of distributed point source approach for infrasound simulations by direct comparison with acoustic observations in the field experiment.

58 GEOSCIENCES↗

DOE OSTI2025

Simulating High-Frequency Seismograms in Realistic Earth Models to Better Understand Source Discrimination Based on Differential Magnitudes ( M L− M c)

Discriminating low-yield underground nuclear explosions from small earthquakes is a key task in monitoring nuclear test ban treaties. P/S amplitude ratios have been an effective discriminant for moderate-sized events recorded at regional distances, but it is unclear if they are as effective in discriminating small seismic events recorded at local distances (<150 km). The difference between local magnitude (M L ) and coda duration magnitude (M c ) has been proposed as a new discriminant that may complement P/S amplitude ratios at local distances. Here, in this work, we calculate high-frequency (up to ∼4 Hz) synthetic seismograms at epicentral distances of 0–30 km in realistic models of the Salt Lake basin (Utah, United States) to better understand how variations in source type and depth affect M L −M c values. The Earth models incorporate simplified 1D and deterministic 3D structures, small-wavelength stochastic velocity perturbations, and surface topography. Coda waves are enhanced for the more complicated models compared to the base 1D model, but still underpredict observed durations by about a factor of two, which results in overprediction of amplitude to duration ratios (i.e., M L −M c values) for a near-surface explosion and a 7 km deep earthquake. For both source types, the predicted M L and M c values decrease as source depth increases, and M L −M c shows only minor variation with depth; however, M L −M c is on average ∼0.5 units smaller for explosions than earthquakes. This finding may imply that M L −M c has sensitivity to source type, in addition to being a depth discriminant, but more modeling is needed given the limitations of the current study. Future modeling should incorporate higher-frequency (≳5 Hz) simulations over a larger distance range (0–150 km), where M L and M c are commonly measured, while honoring low shear velocities (<300 m/s) near the surface and sampling a wider range of earthquake and explosion source mechanisms.

Hutchings, Sean J. [Univ. of Utah, Salt Lake City,↗

DOE OSTI2025

Acoustic Observations of the OSIRIS-REx Sample Return Capsule Re-Entry from Wendover Airport

The Origins, Spectral Interpretation, Resource Identification, and Security‐Regolith Explorer sample return capsule (SRC) re‐entered the Earth’s atmosphere at hypersonic speeds from interplanetary space on 24 September 2023. The current work reports on 18 ground‐based acoustic sensors deployed at Wendover Airport, the same location that the Genesis and Stardust SRC re‐entries were recorded. Four different sensors (Chaparral Physics, Gem, Wilson Engineering Research and Development [WERD], and RedVox) were deployed in close proximity to compare their performance. All the sensors captured an N‐wave signal associated with the SRC re‐entry shock wave followed by a broadband coda. The Chaparral Physics array served as the high‐fidelity reference measurement. The N‐wave signal had a peak‐to‐peak amplitude of 4.07 Pa with a fundamental frequency of 4.98 Hz from 167.5° measured clockwise from north, nearly perpendicular to the SRC trajectory. In addition, high coherence in the coda was shown to be associated reflections from the surrounding mountains. In general, the more economical sensors (Gem, WERD, and RedVox) produced results that were consistent with these observations and sensor specifications. Beamforming with these single sensors arranged as an array showed agreement with the high‐fidelity array to within a couple of degrees. Furthermore, the current high‐fidelity results were compared with the measurements during the Genesis and Stardust SRC re‐entries. All three entries produced a broadband fundamental peak at a frequency that was inversely related to the SRC diameter as well as evidence of reflections from the surrounding topography.

KC, Real J. [Oklahoma State University, Stillwater↗

DOE OSTI2023

Topographic Evolution of Nb Surfaces Subjected to Electropolishing

Surface quality plays an important role in the performance of SRF cavities. We report the topographic evolution of H degassed, N and O doped Nb surfaces subjected to electropolishing using atomic force microscopy. N doped Nb shows a substantial increase in roughness due to formation of niobium nitrides into the surface. Implications of the topography on magnetic field enhancement and vortex entry are discussed.

Lechner, E.M.↗

DOE OSTI2020

A GPU-based 2D shallow water quality model

In this paper, a 2D shallow water flow solver integrated with a water quality model is presented. The interaction between the main water quality constituents included is based on the Water Quality Analysis Simulation Program. Efficiency is achieved by computing with a combination of a Central Processing Unit (CPU) and a Graphics Processing Unit (GPU) device. This technique is intended to provide robust and accurate simulations with high computation speedups with respect to a single-core CPU in real events. The proposed numerical model is evaluated in cases that include the transport and reaction of water quality components over irregular bed topography and dry–wet fronts, verifying that the numerical solution in these situations conserves the required properties (C-property and positivity). The model can operate in any steady or unsteady form allowing an efficient assessment of the environmental impact of water flows. The field data from an unsteady river reach test case are used to show that the model is capable of predicting the measured temporal distribution of dissolved oxygen and water temperature, proving the robustness and computational efficiency of the model, even in the presence of noisy signals such as wind speed.

54 ENVIRONMENTAL SCIENCES↗