USGS ScienceSearch

SEARCH · USGS Science

Results for “Computers and Fluids”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 19 recordsLinked to original sources

Development of computational fluid dynamics--habitat suitability (CFD-HSI) models to identify potential passage--Challenge zones for migratory fishes in the Penobscot River

A two-dimensional computational fluid dynamics-habitat suitability (CFD–HSI) model was developed to identify potential zones of shallow depth and high water velocity that may present passage challenges for five anadromous fish species in the Penobscot River, Maine, upstream from two existing dams and as a result of the proposed future removal of the dams. Potential depth-challenge zones were predicted for larger species at the lowest flow modeled in the dam-removal scenario. Increasing flows under both scenarios increased the number and size of potential velocity-challenge zones, especially for smaller species. This application of the two-dimensional CFD–HSI model demonstrated its capabilities to estimate the potential effects of flow and hydraulic alteration on the passage of migratory fish.

Maine

Computational fluid dynamics-habitat suitability index (CFD-HSI) modelling as an exploratory tool for assessing passability of riverine migratory challenge zones for fish

We developed two-dimensional computational fluid hydraulics-habitat suitability index (CFD-HSI) models to identify and qualitatively assess potential zones of shallow water depth and high water velocity that may present passage challenges for five major anadromous fish species in a 2.63-km reach of the main stem Penobscot River, Maine, as a result of a dam removal downstream of the reach. Suitability parameters were based on distribution of fish lengths and body depths and transformed to cruising, maximum sustained and sprint swimming speeds. Zones of potential depth and velocity challenges were calculated based on the hydraulic models; ability of fish to pass a challenge zone was based on the percent of river channel that the contiguous zone spanned and its maximum along-current length. Three river flows (low: 99.1 m 3 sec -1 ; normal: 344.9 m 3 sec -1 ; and high: 792.9 m 3 sec -1 ) were modelled to simulate existing hydraulic conditions and hydraulic conditions simulating removal of a dam at the downstream boundary of the reach. Potential depth challenge zones were nonexistent for all low-flow simulations of existing conditions for deeper-bodied fishes. Increasing flows for existing conditions and removal of the dam under all flow conditions increased the number and size of potential velocity challenge zones, with the effects of zones being more pronounced for smaller species. The two-dimensional CFD-HSI model has utility in demonstrating gross effects of flow and hydraulic alteration, but may not be as precise a predictive tool as a three-dimensional model. Passability of the potential challenge zones cannot be precisely quantified for two-dimensional or three-dimensional models due to untested assumptions and incomplete data on fish swimming performance and behaviours.

River Research and Applications

Benchmarking computational fluid dynamics models of lava flow simulation for hazard assessment, forecasting, and risk management

Numerical simulations of lava flow emplacement are valuable for assessing lava flow hazards, forecasting active flows, designing flow mitigation measures, interpreting past eruptions, and understanding the controls on lava flow behavior. Existing lava flow models vary in simplifying assumptions, physics, dimensionality, and the degree to which they have been validated against analytical solutions, experiments, and natural observations. In order to assess existing models and guide the development of new codes, we conduct a benchmarking study of computational fluid dynamics (CFD) models for lava flow emplacement, including VolcFlow, OpenFOAM, FLOW-3D, COMSOL, and MOLASSES. We model viscous, cooling, and solidifying flows over horizontal planes, sloping surfaces, and into topographic obstacles. We compare model results to physical observations made during well-controlled analogue and molten basalt experiments, and to analytical theory when available. Overall, the models accurately simulate viscous flow with some variability in flow thickness where flows intersect obstacles. OpenFOAM, COMSOL, and FLOW-3D can each reproduce experimental measurements of cooling viscous flows, and OpenFOAM and FLOW-3D simulations with temperature-dependent rheology match results from molten basalt experiments. We assess the goodness-of-fit of the simulation results and the computational cost. Our results guide the selection of numerical simulation codes for different applications, including inferring emplacement conditions of past lava flows, modeling the temporal evolution of ongoing flows during eruption, and probabilistic assessment of lava flow hazard prior to eruption. Finally, we outline potential experiments and desired key observational data from future flows that would extend existing benchmarking data sets.

Journal of Applied Volcanology

A computational fluid dynamics modeling study of guide walls for downstream fish passage

A partial-depth, impermeable guidance structure (or guide wall) for downstream fish passage is typically constructed as a series of panels attached to a floating boom and anchored across a water body (e.g. river channel, reservoir, or power canal). The downstream terminus of the wall is generally located nearby to a fish bypass structure. If guidance is successful, the fish will avoid entrainment in a dangerous intake structure (i.e. turbine intakes) while passing from the headpond to the tailwater of a hydroelectric facility through a safer passage route (i.e. the bypass). The goal of this study is to determine the combination of guide wall design parameters that will most likely increase the chance of surface-oriented fish being successfully guided to the bypass. To evaluate the flow field immediately upstream of a guide wall, a parameterized computational fluid dynamics model of an idealized power canal was constructed in © ANSYS Fluent v 14.5 (ANSYS Inc., 2012). The design parameters investigated were the angle and depth of the guide wall and the average approach velocity in the power canal. Results call attention to the importance of the downward to sweeping flow ratio and demonstrate how a change in guide wall depth and angle can affect this important hydraulic cue to out-migrating fish. The key findings indicate that a guide wall set at a small angle (15° is the minimum in this study) and deep enough such that sweeping flow dominant conditions prevail within the expected vertical distribution of fish approaching the structure will produce hydraulic conditions that are more likely to result in effective passage.

Ecological Engineering

Computational fluid dynamics simulations of the Late Pleistocene Lake Bonneville flood

At approximately 18.0 ka, pluvial Lake Bonneville reached its maximum level. At its northeastern extent it was impounded by alluvium of the Marsh Creek Fan, which breached at some point north of Red Rock Pass (Idaho), leading to one of the largest floods on Earth. About 5320 km 3 of water was discharged into the Snake River drainage and ultimately into the Columbia River. We use a 0D model and a 2D non-linear depth-averaged hydrodynamic model to aid understanding of outflow dynamics, specifically evaluating controls on the amount of water exiting the Lake Bonneville basin exerted by the Red Rock Pass outlet lithology and geometry as well as those imposed by the internal lake geometry of the Bonneville basin. These models are based on field evidence of prominent lake levels, hypsometry and terrain elevations corrected for post-flood isostatic deformation of the lake basin, as well as reconstructions of the topography at the outlet for both the initial and final stages of the flood. Internal flow dynamics in the northern Lake Bonneville basin during the flood were affected by the narrow passages separating the Cache Valley from the main body of Lake Bonneville. This constriction imposed a water-level drop of up to 2.7 m at the time of peak-flow conditions and likely reduced the peak discharge at the lake outlet by about 6%. The modeled peak outlet flow is 0.85·10 6 m 3 s −1 . Energy balance calculations give an estimate for the erodibility coefficient for the alluvial Marsh Creek divide of ∼0.005 m y −1 Pa −1.5 , at least two orders of magnitude greater than for the underlying bedrock at the outlet. Computing quasi steady-state water flows, water elevations, water currents and shear stresses as a function of the water-level drop in the lake and for the sequential stages of erosion in the outlet gives estimates of the incision rates and an estimate of the outflow hydrograph during the Bonneville Flood: About 18 days would have been required for the outflow to grow from 10% to 100% of its peak value. At the time of peak flow, about 10% of the lake volume would have already exited; eroding about 1 km 3 of alluvium from the outlet, and the lake level would have dropped by about 10.6 m.

Lake Bonneville

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona

Predictive techniques for river channel evolution and maintenance

Predicting changes in alluvial channel morphology associated with anthropogenic and natural changes in flow and/or sediment supply is a critical part of the management of riverine systems. Over the past few years, advances in the understanding of the physics of sediment transport in conjunction with rapidly increasing capabilities in computational fluid dynamics have yielded now approaches to problems in river mechanics. Techniques appropriate for length scales ranging from reaches to bars and bedforms are described here. Examples of the use of these computational approaches are discussed for three cases: (1) the design of diversion scenarios that maintain channel morphology in steep cobble-bedded channels in Colorado, (2) determination of channel maintenance flows for the preservation of channel islands in the Snake River in Idaho, and (3) prediction of the temporal evolution of deposits in lateral separation zones for future assessment of the impacts of various dam release scenarios on lateral separation deposits in the Colorado River in Grand Canyon. With continued development of their scientific and technical components, the methodologies described here can provide powerful tools for the management of river environments in the future.

Water, Air, & Soil Pollution

Real-time lava flow forecasting during the 2022 Mauna Loa eruption response

On November 27, 2022, Mauna Loa (Hawai‘i) erupted for the first time in 38 years, initially producing lava flows that covered the floor of its summit caldera, Moku‘āweoweo. Over the first 12 h following the summit eruption, four main fissures opened on Mauna Loa’s Northeast Rift Zone, with “fissure 3” quickly becoming the dominant source of lava flows. For the next 12 days, fissure 3 produced a 19-km-long lava flow to the north, crossing the Mauna Loa Weather Observatory access road and coming within 2.8 km of inundating the Daniel K. Inouye Highway (Saddle Road). Within 40 min of fissure 3 opening, inundation modeling efforts had begun. For the duration of the eruption, the computational fluid dynamics model Lava2d was run in real time, using flow front locations and other field observations to make sequential forecast improvements, eventually producing a set of models which accurately predicted the routing and arrival times of lava from fissure 3. These models were used to inform timing estimates of possible future inundation of the Saddle Road. As the eruption progressed, almost 4000 Lava2d models were made using high-performance computing resources, providing critical information on uncertainty in multi-week forecasts. To our knowledge, this was the first ever real-time physics-based ensemble lava flow modeling and forecasting effort. Here, we present a chronology of these real-time efforts, focusing on the successes and limitations of this approach.

Hawaii

Hydrologic models of modern and fossil geothermal systems in the Great Basin: Genetic implications for epithermal Au-Ag and Carlin-type gold deposits

The Great Basin region in the western United States contains active geothermal systems, large epithermal Au-Ag deposits, and world-class Carlin-type gold deposits. Temperature profiles, fluid inclusion studies, and isotopic evidence suggest that modern and fossil hydrothermal systems associated with gold mineralization share many common features, including the absence of a clear magmatic fluid source, discharge areas restricted to fault zones, and remarkably high temperatures (>200 ??C) at shallow depths (200-1500 m). While the plumbing of these systems varies, geochemical and isotopic data collected at the Dixie Valley and Beowawe geothermal systems suggest that fluid circulation along fault zones was relatively deep (>5 km) and comprised of relatively unexchanged Pleistocene meteoric water with small (<2.5%) shifts from the meteoric water line (MWL). Many fossil ore-forming systems were also dominated by meteoric water, but usually exhibit ??18O fluid-rock interactions with larger shifts of 5???-20??? from the MWL. Here we present a suite of two-dimensional regional (100 km) and local (40-50 km) scale hydrologic models that we have used to study the plumbing of modern and Tertiary hydrothermal systems of the Great Basin. Geologically and geophysically consistent cross sections were used to generate somewhat idealized hydrogeologic models for these systems that include the most important faults, aquifers, and confining units in their approximate configurations. Multiple constraints were used, including enthalpy, ??18O, silica compositions of fluids and/or rocks, groundwater residence times, fluid inclusion homogenization temperatures, and apatite fission track anomalies. Our results suggest that these hydrothermal systems were driven by natural thermal convection along anisotropic, subvertical faults connected in many cases at depth by permeable aquifers within favorable lithostratigraphic horizons. Those with minimal fluid ?? 18O shifts are restricted to high-permeability fault zones and relatively small-scale (???5 km), single-pass flow systems (e.g., Beowawe). Those with intermediate to large isotopic shifts (e.g., epithermal and Carlin-type Au) had larger-scale (???15 km) loop convection cells with a greater component of flow through marine sedimentary rocks at lower water/rock ratios and greater endowments of gold. Enthalpy calculations constrain the duration of Carlin-type gold systems to probably <200 k.y. Shallow heat flow gradients and fluid silica concentrations suggest that the duration of the modern Beowawe system is <5 k.y. However, fluid flow at Beowawe during the Quaternary must have been episodic with a net duration of ???200 k.y. to account for the amount of silica in the sinter deposits. In the Carlin trend, fluid circulation extended down into Paleozoic siliciclastic rocks, which afforded more mixing with isotopically enriched higher enthalpy fluids. Computed fission track ages along the Carlin trend included the convective effects, and ranged between 91.6 and 35.3 Ma. Older fission track ages occurred in zones of groundwater recharge, and the younger ages occurred in discharge areas. This is largely consistent with fission track ages reported in recent studies. We found that either an amagmatic system with more permeable faults (10-11 m2) or a magmatic system with less permeable faults (10-13 m2) could account for the published isotopic and thermal data along the Carlin trend systems. Localized high heat flow beneath the Muleshoe fault was needed to match fl uid inclusion temperatures at Mule Canyon. However, both magmatic and amagmatic scenarios require the existence of deep, permeable faults to bring hot fluids to the near surface. ?? 2008 Geological Society of America.

Geosphere

Comparison of hydraulics and particle removal efficiencies in a mixed cell raceway and Burrows pond rearing system

We compared the hydrodynamics of replicate experimental mixed cell and replicate standard Burrows pond rearing systems at the Dworshak National Fish Hatchery, ID, in an effort to identify methods for improved solids removal. We measured and compared the hydraulic residence time, particle removal efficiency, and measures of velocity using several tools. Computational fluid dynamics was used first to characterize hydraulics in the proposed retrofit that included removal of the traditional Burrows pond dividing wall and establishment of four counter rotating cells with appropriate drains and inlet water jets. Hydraulic residence time was subsequently established in the four full scale test tanks using measures of conductivity of a salt tracer introduced into the systems both with and without fish present. Vertical and horizontal velocities were also measured with acoustic Doppler velocimetry in transects across each of the rearing systems. Finally, we introduced ABS sinking beads that simulated fish solids then followed the kinetics of their removal via the drains to establish relative purge rates. The mixed cell raceway provided higher mean velocities and a more uniform velocity distribution than did the Burrows pond. Vectors revealed well-defined, counter-rotating cells in the mixed cell raceway, and were likely contributing factors in achieving a relatively high particle removal efficiency-88.6% versus 8.0% during the test period. We speculate retrofits of rearing ponds to mixed cell systems will improve both the rearing environments for the fish and solids removal, improving the efficiency and bio-security of fish culture. We recommend further testing in hatchery production trials to evaluate fish physiology and growth.

Idaho

Combining field observations and high-resolution numerical modeling to demonstrate the effect of coral reef roughness on turbulence and its implications for reef restoration design

Coral reefs are effective natural barriers that protect adjacent coastal communities from hazards such as erosion and storm-induced flooding. However, the degradation of coral reefs compromises their ability to protect against these hazards, making degraded reefs a target for restoration. There have been limited field and numerical modeling studies conducted to understand how an increase in coral reef roughness, as would occur due to restoration, can affect wave energy dissipation for a range of real-world wave and water level conditions. To address this knowledge gap, field measurements were collected over adjacent low-roughness and high-roughness reefs off Molokaʻi, Hawaiʻi, USA, subjected to the same oceanographic forcing. Those field data were then used to calibrate and validate OpenFOAM computational fluid dynamics models of the reef. These calibrated models were then used to explore energy dissipation for a range of wave conditions based on measurements from a suite of existing datasets and values from the literature. In general, wave dissipation scales with incident wave conditions, where greater dissipation occurred for shallow depths and shorter-period waves. This tendency for short-period waves to be more readily attenuated is supported by wave energy dissipation factors in the range of 0.1–5, which decline with increasing wave period. Near-bed turbulent kinetic energy dissipation also scales with incident wave conditions, where the greatest difference in dissipation between low and high relief cases occurs for short wave periods. Turbulence becomes less affected by bottom roughness as the wave period increases. Based on this study, wave attenuation and turbulent energy dissipation could be enhanced by 0.5–1 order of magnitude (45% per across-shore meter) if the seabed roughness at the field site were increased by 13%, an achievable goal in coral reef restoration.

Coastal Engineering

A model for evaluation of sediment exposure and burial for freshwater mussels from heavy particle sedimentation

Freshwater mussels (Bivalvia: Unionida) are an ecologically important faunal group. Excessive sediments, both in suspended and deposited formats, are believed to have negative effects on survival of freshwater mussels. However, there is a lack of quantitative tools for assessing the impact of abrupt and excessive sedimentation on freshwater mussel habitats. This gap in knowledge poses challenges for construction planning that necessitates evaluating sedimentation effects on mussels. In this paper, we present a simple Lagrangian particle tracking (LPT) model designed to investigate the downstream distances from the sediment release point where mussels may face risks of sediment exposure and burial during episodic sedimentation events. We validated the model in predicting the deposition of sands ranging from 125 to 625 &#x3BC; m "> 625μm in diameter within an open channel. Subsequently, we applied the LPT model to a representative stretch of mussel habitats in the lower Osage River, Missouri, USA. Implementing the LPT model using computational fluid dynamics analysis of river hydrodynamics, we investigated the downstream distances impacted by a hypothetical sedimentation event. The results indicate that the mussels within tens and hundreds of meters may be at risk of burial under low and high flow conditions, respectively. Mussels within hundreds and potentially thousands of meters may be at risks of sediment exposure under low and high flow conditions, respectively. The impacted distance is highly dependent on the sediment diameters and flow conditions. These results contribute valuable insights into the potential impacts of sedimentation on freshwater mussel habitats, aiding in informed decision-making for construction projects near these critical aquatic ecosystems.

Ecological Modelling

Ecological inferences on invasive carp survival using hydrodynamics and egg drift models

Bighead carp ( Hypophthalmichthys nobilis ), silver carp ( H. molitrix ), black carp ( Mylopharyngodon piceus ), and grass carp ( Ctenopharyngodon idella ), are invasive species in North America. However, they hold significant economic importance as food sources in China. The drifting stage of carp eggs has received great attention because egg survival rate is strongly affected by river hydrodynamics. In this study, we explored egg-drift dynamics using computational fluid dynamics (CFD) models to infer potential egg settling zones based on mechanistic criteria from simulated turbulence in the Lower Missouri River. Using an 8-km reach, we simulated flow characteristics with four different discharges, representing 45–3% daily flow exceedance. The CFD results elucidate the highly heterogeneous spatial distribution of flow velocity, flow depth, turbulence kinetic energy (TKE), and the dissipation rate of TKE. The river hydrodynamics were used to determine potential egg settling zones using criteria based on shear velocity, vertical turbulence intensity, and Rouse number. Importantly, we examined the difference between hydrodynamic-inferred settling zones and settling zones predicted using an egg-drift transport model. The results indicate that hydrodynamic inference is useful in determining the ‘potential’ of egg settling, however, egg drifting paths should be taken into account to improve prediction. Our simulation results also indicate that the river turbulence does not surpass the laboratory-identified threshold to pose a threat to carp eggs.

Scientific Reports

Errors in acoustic doppler profiler velocity measurements caused by flow disturbance

Acoustic Doppler current profilers (ADCPs) are commonly used to measure streamflow and water velocities in rivers and streams. This paper presents laboratory, field, and numerical model evidence of errors in ADCP measurements caused by flow disturbance. A state-of-the-art three-dimensional computational fluid dynamic model is validated with and used to complement field and laboratory observations of flow disturbance and its effect on measured velocities. Results show that near the instrument, flow velocities measured by the ADCP are neither the undisturbed stream velocity nor the velocity of the flow field around the ADCP. The velocities measured by the ADCP are biased low due to the downward flow near the upstream face of the ADCP and upward recovering flow in the path of downstream transducer, which violate the flow homogeneity assumption used to transform beam velocities into Cartesian velocity components. The magnitude of the bias is dependent on the deployment configuration, the diameter of the instrument, and the approach velocity, and was observed to range from more than 25% at 5cm from the transducers to less than 1% at about 50cm from the transducers for the scenarios simulated. ?? 2007 ASCE.

Journal of Hydraulic Engineering

Behavioral responses of sea lamprey (Petromyzon marinus) and white sucker (Catostomus commersonii) to turbulent flow during fishway passage attempts

An understanding of how undesirable and desirable fish species respond behaviorally to turbulent flow in fishways would guide development of selective fish passage techniques. We applied high-resolution computational fluid dynamics modeling and competing risks analysis towards the development of predictive selective passage models. Sea lamprey ( Petromyzon marinus ; an invasive fish in the Great Lakes Basin, North America) upstream passage probability declined from 0.73 to 0.03 as flow conditions became increasingly turbulent, while declines in white sucker ( Catostomus commersonii , a native fish in the region) upstream passage probability were less substantial (0.53 to 0.44). Deploying a sea lamprey trap in the fishway did not effectively reduce sea lamprey upstream passage probability, though capture rate increased during trials with cooler water temperature and low total kinetic energy. Bifurcated fishways that maintain low turbulent flow in the entrapment route and high turbulent flow in the upstream passage route could increase the effectiveness of trapping sea lamprey in fishways as a means to advance selective passage goals.

Wisconsin

Evaluation of the fish passage effectiveness of the Bonneville I prototype surface collector using three-dimensional ultrasonic fish tracking - Final Report

This report describes tests conducted at Bonneville Dam on the Columbia River in the spring of 2000. The studies used three-dimensional (3D) acoustic telemetry and computational fluid dynamics (CFD) hydraulic modeling techniques to evaluate the response of outmigrating juvenile steelhead (Oncorhynchus mykiss) and yearling chinook (O. tshawytscha) to the Prototype Surface Collector (PSC) installed at Powerhouse I of Bonneville Dam in 1998 to test the concept of using a deep-slot surface bypass collector to divert downstream migrating salmon from turbines. The study was conducted by Pacific Northwest National Laboratory (PNNL), the Waterways Experiment Station of the U.S. Army Corp of Engineers (COE), Asci Corporation, and the U.S. Geological Survey (USGS), and was sponsored by COE’s Portland District. The goal of the study was to observe the three-dimensional behavior of tagged fish (fish bearing ultrasonic micro-transmitters) within 100 meters (m) of the surface flow bypass structure to test hypotheses about the response of migrants to flow stimuli generated by the presence of the surface flow bypass prototype and its operation. Research was done in parallel with radio telemetry studies conducted by USGS and hydroacoustic studies conducted by WES & Asci to evaluate the prototype surface collector.

Report

Analysis of alternative weir designs for improved passage of select fish at the U.S. Geological Survey streamgaging weir at Blackwells Mills, New Jersey

As the population of New Jersey continues to remain dense, the need for water supply will likely continue to be high, which can lead to water managers needing to make difficult decisions about managing drinking-water supply. Streamgaging weirs like the ones used by the U.S. Geological Survey (USGS) play a critical role in providing accurate and stable streamflow data, but their presence can affect the passage of diadromous fish species such as river herring ( Alosa pseudoharengus [alewife], Alosa aestivalis [blueback herring], and Alosa sapidissima [American shad]). In some situations, weirs existing in rivers and streams are no longer used because they were part of a farm irrigation system or some type of industrial operation. The weir at the USGS streamgage 01402000 Millstone River at Blackwells Mills, New Jersey, was purposefully built as a hydraulic-control structure that provides a precise and stable control for the measurement of stage and computation of continuous streamflow. To satisfy the dual need of maintaining accurate streamflow data and providing improved fish passage for select species of fish during migration season, the USGS proposed the development and evaluation of two alternative weir designs that would meet the criteria established for successful passage of American shad, alewife, and blueback herring during their yearly migration. The designs were also required to maintain adequate control of the upstream pool elevation necessary for the precise computation of streamflow used by State agencies for municipal water-supply purposes for surrounding communities. Two alternative weir design modifications were incorporated at the center of the Blackwells Mills weir and modeled using two-dimensional hydraulic modeling software and three-dimensional computational fluid-dynamics software to simultaneously evaluate conditions for passage of the target fish species and effects to streamflow computations at the streamgage. The models were calibrated to existing conditions around the weir location using surveyed-elevation data and recorded stage, streamflow, and velocity in the Millstone River. The alternative weir designs lowered the weir crest by 1.02 feet (ft) and the resulting simulations showed an effective increase in depth of 0.98 ft at the median streamflow of 251 cubic feet per second (ft 3 /s) and 0.96 ft at the 95-percent exceedance streamflow of 98 ft 3 /s. The alternative weir designs were also found to increase streamflow depth across the shallowest portions of the weir structure at the downstream anti-scour skirt by lowering the skirt about 4 inches, allowing for two or more body depths of water for American shad, alewife, and blueback herring at the median migration streamflow of 251 ft 3 /s. The alternative weir designs also reduced the highest stream velocities across the downstream weir sill and anti-scour skirt from about 9 to 10 feet per second, and the depth-averaged velocity to about 7 to 8 feet per second. The sensitivity of the weir with respect to the computation of streamflow was increased from about 1.8 cubic feet per second per hundredth foot to 1.6 cubic feet per second per hundredth foot for streamflows of about 10–100 cubic feet per second.

New Jersey