USGS Publications Warehouse

Sampling and analysis plan for the water-quality monitoring program in Lake Koocanusa and upper Kootenai River, Montana, water years 2022–23

Fri, 3 Apr 2026 18:10:53

The U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples and environmental data in Lake Koocanusa (also known as “Koocanusa Reservoir”), the Kootenai River, and the Tobacco River during water years 2022–23. The transboundary Lake Koocanusa is in southeastern British Columbia, Canada, and northwestern Montana, United States. It was formed by constructing Libby Dam on the Kootenai River 26 kilometers upstream from Libby, Montana. One of the lake sites and the Kootenai River site, in the Libby Dam tailwater (the outflow of the lake flow into the Kootenai River), were equipped with automated, high-frequency ServoSipper water samplers. At the lake site, these samplers were mounted to pontoon platforms during the summer, and a submersible ServoSipper sipper was deployed with ice buoys during the winter. Samples were automatically collected from multiple depths. At the Kootenai River site, these samplers were housed in the gage house. In water year 2022, discrete water-quality samples were collected every 4–6 weeks, year round, at all four lake sites in the Kootenai River between April and November. In water year 2023, discrete water-quality samples were collected at three lake sites and the Kootenai and Tobacco River sites every 4–6 weeks. The goal of this project was to collect multidepth, high-frequency vertical and temporal water-quality samples and data to understand the limnological and biological processes that control variations and trends in selenium concentrations and loads throughout Lake Koocanusa and in the Libby Dam tailwater at the southern end of the lake. This sampling and analysis plan documents the organization, sampling and data-collection scheme and design, pre- and post-collection processes, and quality-assurance and quality-control procedures of the Koocanusa/Kootenai water-quality monitoring program during water years 2022–23.

Methods for estimating daily upstream location of the freshwater- saltwater interface along the Maurice and Cohansey Rivers, New Jersey

Closson, Jennifer; Suro, Thomas; Niemoczynski, Lukasz — Fri, 3 Apr 2026 17:27:03

The Delaware River basin (DRB) provides drinking water to 15 million people in the surrounding area. Water is frequently withdrawn from the freshwater reaches of streams, above head of tide, in the DRB for use as public drinking water. During extended periods of low flow, saltwater can move upstream, which can threaten drinking-water supplies in the basin. Due to spatial patterns in bathymetry, tidal influences within the DRB, and varying weather conditions, it can be hard to predict the movement and upstream extent of the freshwater-saltwater interface, often defined as the salt-front. Although there is a relationship that predicts this location in the main stem of the Delaware River, there lacks a relationship for its tributaries, such as the Maurice and Cohansey Rivers in southwestern New Jersey. In this study, a relationship was developed between daily specific conductance (SC) at gage locations along the tidal river reaches of the Maurice and Cohansey Rivers to the daily upstream location of the salt-front. The study augmented existing real-time tide gage data with the collection of water temperature and specific conductance data to develop the relationship. Additionally, longitudinal profiles upstream of the selected tide gages were conducted during a range of high tide conditions to define the location of the salt-front. Equations were then developed that related the daily SC measured at the tide gage to the upstream location of the salt-front. The equations were used to estimate the daily upstream location of the salt-front for the period of July 15, 2021, to July 15, 2024. This work can aid in understanding the propagation of saltwater upstream, which can affect local communities and crop farmers along these tidal reaches of the DRB.

Thickness and other characteristics of overbank sediment deposited during an extreme flood in May 1978 along the Powder River, Montana

Moody, John; Meade, Robert — Fri, 3 Apr 2026 17:29:50

An extreme flood on the Powder River in southeastern Montana in May 1978 inundated its valley and deposited sediment on the floodplains and terraces at multiple heights. The recurrence interval for this flood was less than 1 percent in the reach between Moorhead and Broadus, Montana. Peak discharges at the U.S. Geological Survey streamgages at Moorhead and Broadus were 779 and 711 cubic meters per second (m³/s), respectively, the difference reflecting the water and sediment stored on the valley surfaces. Bankfull discharge depended on the height of the bank at the start of the valley transect and varied from 243 to 713 m³/s. Sediment-thickness and particle-size data were collected and analyzed in the autumn of 1978 by U.S. Geological Survey scientists at about 900 sites along 20 valley transects between Moorhead and Broadus, Mont. These transects were approximately orthogonal to the floodflow across the floodplain from near the edge of the channel to the high-water mark. Estimated maximum flood depths along these transects ranged from 0.9 to 4.2 meters.

Contrary to theory and controlled laboratory experiments, the distribution of sediment thickness and particle sizes along valley transects did not decrease systematically with distance from the main channel but were affected by the distribution of vegetation. Additionally, some water and sediment—primarily muds and silts—were conveyed by subsidiary channels (often connected to the main channel downriver from the valley transect) during the early stages of the flood before water overtopped the banks at the start of the valley transect. The vegetation created natural sediment traps in the recirculation and wake zones in the lee of trees and shrubs. Sediment that accumulated in these traps formed dunes and thus an undulating surface with many local maximums and minimums in sediment thicknesses. Sediment in the traps are referred to as lee dunes, which recorded flow conditions and a predominance of coarsening-upward sequence of particle sizes (mud to silt to sands) starting at the preflood surface. These sequences were associated with the rising limb of the hydrograph, and later as the flood began to recede, the lee dunes recorded a fining-upward sequence associated with the falling limb of the hydrograph.

Rapid seismic and infrasound assessment of large landslides: A case study from Denali National Park and Preserve (Alaska)

Thu, 26 Mar 2026 20:12:07

Large, rapid landslides are a global hazard that can occur in remote, mountainous areas. Eyewitness reports of landslides and satellite imagery can often be limited or delayed, particularly during inclement weather. However, landslide-generated seismic and infrasound (low-frequency atmospheric sound) waves can be remotely detected in near real-time. This information can significantly expedite characterization and possible landslide response activities. Here, we highlight these capabilities using a > 4 million m³ ice–rock avalanche in Denali National Park and Preserve (Alaska). This event was detected via a landslide-specific seismic location and volume estimation algorithm deployed in Alaska, and — notably — by standard earthquake monitoring systems. Following rapid detection of this event, we combined its seismic and infrasound dataset with optical, synthetic aperture radar, and oblique aerial imagery, multitemporal digital elevation models, and a numerical flow model to reconstruct its failure timeline and dynamics. We apply array processing to infrasound signals traveling > 250 km and find that two precursory events occurred minutes prior to the main failure. We use long-period seismic signals to infer the force exerted by the landslide on the Earth and constrain the rheological parameters of our numerical flow simulation with this result and deposit morphology. The main failure produced a steeply-dipping impulsive initial downward force and reached speeds exceeding 60 m/s. This impulsive force generated relatively strong seismic body waves, which contributed to the earthquake system detection. This large, remote Alaska landslide underscores the key value of seismic and infrasound analysis for rapid landslide assessment and motivates efforts to further operationalize these approaches.

Groundwater budget for the Mountain Home area, southern Idaho, 2022–23

Thomas, Paul — Fri, 27 Feb 2026 21:32:48

The U.S. Geological Survey, with funding from the Idaho Department of Water Resources, developed a groundwater budget for the Mountain Home area in southern Idaho for irrigation year 2023 (November 1, 2022–October 31, 2023). This study focused on the water balance across the Cinder Cone Butte Critical Groundwater Area (CGWA), Mountain Home Groundwater Management Area (GWMA), and the rest of the study area (RoSA), compiling data from various sources, including precipitation records, groundwater level measurements, metered groundwater pumpage data, surface water diversions and evapotranspiration (ET) estimates derived from remote sensing satellite imagery, and ground-based reference data. Key inflow components included recharge from applied surface water irrigation (which incorporates incidental recharge from irrigation practices and conveyance losses), estimated tributary streamflow, and estimated mountain block recharge. The key outflow components were groundwater pumpage for irrigation, municipal, industrial, and domestic uses, and ET. Recharge from applied irrigation and mountain block recharge were the largest inflows, and groundwater pumpage for irrigation was the largest outflow.

The CGWA had a positive groundwater budget residual of 2,170 acre-feet (acre-ft), which contrasts with observed long-term groundwater level declines and historical trends of storage depletion. This positive residual is likely associated with unquantified outflows, including lateral groundwater flow out of the subregion, or other complexities, such as overestimated tributary contributions relative to the actual recharge for the 2023 water budget. The GWMA exhibited a positive residual of 56,563 acre-ft, primarily owing to recharge from applied surface water irrigation and areal recharge during a wetter-than-average year, which allowed irrigation entities to deliver more water from in-basin and out-of-basin reservoirs. The RoSA showed a large positive residual of 124,933 acre-ft. The interpretation of these positive residuals must account for significant uncertainties, including estimations of areal recharge, tributary streamflow (particularly losses and diversions), ET, the volume of surface water loss to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and the unquantified groundwater discharge to the Snake River. These uncertainties, in combination with the complex hydrogeologic controls on water movement and limitations of remotely sensed data, directly affect the accuracy of water availability assessments.

Future data collection efforts would help reduce these uncertainties and support water resource management decisions in the Mountain Home area. Key efforts could include installing additional streamflow gaging stations (particularly to quantify tributary losses and gains and surface water losses to the Snake River), improving groundwater pumpage metering, and validating remotely sensed ET data with ground-based measurements. Furthermore, to better quantify unrepresented or highly uncertain fluxes, focused investigations on groundwater discharge to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and a more robust determination of the actual influence and volume of mountain block recharge would help refine future water availability assessments for the Mountain Home area.

Erosion potential and flood vulnerability of streams and stream crossings at Acadia National Park, Maine

Armstrong, Ian; McCallister, Meghan; Hyslop, Kristina; Benthem, Adam — Fri, 27 Feb 2026 21:38:18

Acadia National Park has had increases in the frequency and magnitude of precipitation in recent years, leading to increased flood flows, stream erosion, and costly infrastructure damage. To improve infrastructure management in a changing climate, the U.S. Geological Survey, in cooperation with the National Park Service, has developed multiple datasets that can help natural resource managers identify stream reaches and stream crossings that have the highest potential for erosion and flood damage within Acadia National Park. To develop these datasets, we first created a lidar- derived hydrography based on a 1- meter digital elevation model and then estimated peak flows at stream crossings and along the stream network using regional regression equations for Maine. We assessed the erosion potential of stream reaches by computing channel morphologic and hydrologic metrics associated with erosive power, such as stream steepness, topographic openness, and percent storage in the contributing watershed. Stream crossing flood vulnerability was assessed by comparing estimated peak flows to stream crossing conveyance capacities. Our results indicate that stream reaches in the headwaters of the Acadia National Park highlands such as Sargent, Penobscot, and Cadillac Mountain, have the highest erosion potential and generally coincide with reaches that have had erosion and infrastructure damage in the past. Stream crossings with the highest flood vulnerability are distributed throughout Mount Desert Island and Acadia National Park, especially south of Jordan Pond, north of Sargent Mountain, and surrounding Eagle Lake. Over a quarter of the total stream crossings have insufficient information to compute flood vulnerability and are often on the parts of the stream with the highest potential for erosion. The datasets allow users to identify stream reaches with the highest erosion potential, stream crossings that are most vulnerable to flood damage, and to highlight areas where supplemental field assessments could most effectively be completed.

Methods for estimating selected streamflow statistics at ungaged sites in Wyoming based on data through water year 2021

Taylor, Nicholas; Sando, Roy — Fri, 27 Feb 2026 21:40:16

The U.S. Geological Survey, in cooperation with the Wyoming Water Development Office, developed regional regression equations based on basin characteristics and streamflow statistics for streamgages through water year 2021 (October 1, 2020, to September 30, 2021). The regression equations allow estimates of mean annual maximum, mean annual, mean seasonal, and mean monthly streamflows; frequency statistics for the 7- day mean low flows with 2- year and 10- year recurrence intervals, 14- and 30- day mean low flows with 5- year recurrence intervals, and 60- and 1- day mean high flow with 2- year and 5- year recurrence intervals, respectively; and the 0.1- , 0.2- , 0.5- , 1- , 2- , 4- , 5- , 10- , 20- , 25- , 30- , 50- , 60- , 70- , 75- , 80- , 90- , 95- , 98- , and 99- percent durations for annual streamflows and 0.1- , 0.5- , 10- , 15- , 20- , 25- , 30- , 40- , 50- , 60- , 70- , 75- , 80- , 85- , 90- , 95- , and 99- percent durations for monthly streamflows for most months for ungaged locations in Wyoming that are largely unaltered by diversions or upstream reservoirs.

Regression equations were developed for 243 streamflow statistics. Best- subset selection was used to assess explanatory variables for respective streamflow statistics. Exploratory data analyses determined that, of the 81 basin characteristics evaluated as potential explanatory variables, characteristics such as drainage area and precipitation often produced models with the highest adjusted coefficient of determination and lowest mean squared error, as determined in the best- subset selection. To address heteroskedasticity of model residuals, model variables were regionalized using fixed- effects models; the percentages of the streamgage basins in selected ecoregions were defined as interaction terms, which represent the model slope for specific ecoregions. Most models were determined to be statistically significant for probability values less than or equal to 0.1 for one or more regional explanatory variables. The final regional regression equations defined in this report are available for use in the U.S. Geological Survey’s StreamStats web application at https://streamstats.usgs.gov/ss/.

Treatability study to evaluate bioremediation of trichloroethene at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, 2020–22

Fri, 20 Feb 2026 18:18:35

Executive Summary

Chlorinated solvents, including trichloroethene (TCE) and other chlorinated volatile organic compounds (cVOCs), are widespread contaminants that can be treated by bioremediation approaches that enhance anaerobic reductive dechlorination. Reductive dechlorination can be enhanced either through the addition of an electron donor (biostimulation) or the addition of a known dechlorinating culture (bioaugmentation) along with an electron donor. Although bioremediation has been applied at many TCE- contaminated groundwater sites, application in source zones at sites where residual dense nonaqueous phase liquid (DNAPL) is present is more limited. In this study, laboratory and field treatability tests were completed to evaluate the potential application of anaerobic bioremediation for a shallow groundwater plume containing TCE in a perched alluvial aquifer at Site K, former Twin Cities Army Ammunition Plant, Arden Hills, Minnesota, which was on the National Priorities List as the New Brighton/Arden Hills Superfund site until 2019. In addition to the presence of residual DNAPL at the site, temporal variability in groundwater flow directions and input of oxygenated recharge were possible complicating factors for the application of enhanced anaerobic biodegradation in the shallow plume. The Site K plume extends beneath the footprint of Building 103, which was demolished in 2006, and soil excavations to a maximum depth of 6 feet (ft) below ground surface in 2014 were known to leave some deeper contaminated soil in place in the TCE source area. Groundwater treatment at the site, formalized as part of the 1997 Record of Decision, has been in operation since 1986 and consists of an extraction trench at the downgradient edge of the plume to collect groundwater, which is then pumped to an on- site air stripper. Groundwater concentrations in the plume have been relatively stable since treatment began, indicating a continued source of TCE in the aquifer. The desire for a destructive remedy that would enhance the removal of cVOCs in the aquifer at Site K and shorten the remediation timeframe led the U.S. Army to request that the U.S. Geological Survey conduct a groundwater treatability study to assess bioremediation. This report describes the U.S. Geological Survey bioremediation treatability study conducted during 2020–22, including pre- design site characterization to assist in formulating the bioremediation approach, laboratory experiments to support the design of the field pilot test, and implementation and 1-year performance monitoring results for the pilot test.

Pre- design site characterization included the collection of soil cores for cVOC analysis and lithologic descriptions and the re- installment of three wells to obtain hydrologic measurements and initial groundwater chemistry. Relatively flat head gradients were measured at the site, and substantial decreases in water- level elevations occurred from spring to summer (May–July 2021). Continuous water- level monitoring indicated a rapid response to precipitation. Groundwater flow velocities were consistently less than 0.5 foot per day, and the pilot bioremediation test was therefore designed with short lateral distances (about 5 ft) between injection and individual monitoring points. Soil analyses confirmed that high volatile organic compound contamination was left in place in the source area. The highest concentrations were near or in clay at the base of the perched aquifer. Concentrations of cVOCs measured in the replaced wells were consistent with historical data and had a maximum TCE concentration of 57,700 micrograms per liter (μg/L), indicative of nearby residual DNAPL based on the general rule of observed concentrations exceeding 1 percent of solubility. The primary TCE daughter product detected was 1,2- cis- dichloroethene (cisDCE), which indicated limited reductive dechlorination in the plume. Groundwater in both the source and downgradient areas was relatively reducing during the pre- design characterization, particularly in the source area where methane concentrations greater than 400 μg/L were measured.

Initial laboratory tests conducted using native aquifer microorganisms from the three replacement wells showed that anaerobic TCE biodegradation rates were low when biostimulated with the addition of sodium lactate as an electron donor, also known as a carbon donor, and resulted in the production of only cisDCE. Addition of a known dechlorinating culture, WBC- 2, however, resulted in rapid biodegradation and production of ethene, verifying complete reductive dechlorination of TCE. Microcosms constructed with aquifer soil collected from the site were used to evaluate other electron donors besides lactate to support reductive dechlorination by WBC- 2, including corn syrup as an alternative fast- release compound and whey, soy- based vegetable oil, and 3- D Microemulsion (Regenesis, San Clemente, California) as slow-release compounds. First- order rate constants for total organic chlorine removal in these WBC- 2 amended microcosms were greatest with either lactate or vegetable oil as the donor, ranging between 0.061 and 0.047 per day or corresponding half- lives of 11–15 days. Testing of commercial products in other WBC- 2- bioaugmented microcosms led to selection for the field pilot test of an emulsified vegetable oil product that also contained some sodium lactate as a fast- release donor. Delaying the addition of WBC- 2 relative to the donor in the microcosms resulted in the most rapid overall biodegradation rates.

The selected design for the pilot test utilized three separate test plots, each about 30-ft wide and 60-ft long: plots GS1 and GS2 in the source area of the plume and plot GS3 in the downgradient area of the plume near the excavation trench. Each test plot had one injection well, one monitoring well upgradient from the injection point, and 12 surrounding monitoring wells in a grid to capture variable groundwater flow directions. Donor injections, which included a bromide tracer, were completed in October 2021, immediately following baseline sampling, and the WBC- 2 culture was injected about 40 days later, between November 30 and December 2, 2021. Performance monitoring conducted until December 2022 included hydrologic measurements and analyses of cVOCs, redox- sensitive constituents, dissolved organic carbon, bromide, volatile fatty acids, compound- specific carbon isotopes, and microbial communities.

The biogeochemical data collected during the pilot tests in the three treatment plots showed that enhanced, complete reductive dechlorination of cVOCs in the groundwater was achieved in the GS1 and GS3 plots. In contrast, evidence of distribution of the injected amendments and subsequent biodegradation was limited in GS2, which was in an area of more heterogeneous soil lithology and low water table elevations. The molar composition of volatile organic compounds in the GS1 and GS3 plots was dominated by ethene in wells that were reached by the injected amendments by the end of the monitoring period. In the GS1 and GS3 plots, similar patterns were observed of cVOC concentrations decreasing to near detection levels, or below, at some wells sampled in July and October 2022, whereas ethene became dominant and indicated sustained complete reductive dechlorination. Baseline cVOC concentrations were more than a factor of 10 higher in the groundwater in the GS1 plot than in GS3, but no apparent inhibition of complete dechlorination occurred. As expected from the initial pre- design site data and the laboratory experiments, enhanced dissolution of residual DNAPL coupled to biodegradation was evident in the GS1 plot, where a marked increase in dichloroethene (DCE) above the initial baseline and upgradient TCE and DCE concentrations occurred. DCE concentrations subsequently declined where DNAPL dissolution was evident, concurrent with production of vinyl chloride and then predominantly ethene. Thus, overall biodegradation rates outpaced the DNAPL dissolution and desorption and DCE production in the source area. This success in complete degradation to predominantly ethene was achieved even in areas where the DCE concentrations reached a maximum of about 30,000 μg/L. Compound specific isotope analysis of carbon in TCE, cisDCE, trans- 1,2- dichloroethene, and vinyl chloride was conducted to provide another line of evidence of the occurrence and extent of anaerobic biodegradation. Along a flow path in each plot that was affected by the injected amendments, carbon isotopes in the TCE and daughter cVOCs in the groundwater became isotopically heavier, indicating biodegradation.

Enhanced biodegradation rates calculated from the field tests in GS1 and GS3 showed half- lives of 36.9–75.3 days for DCE degradation and 9.48–38.5 days for ethene production. Notably, these ethene production rates calculated from the field tests are consistent with the results of WBC- 2- bioaugmented microcosms amended with either lactate or vegetable oil, which had half- lives for total organic chlorine removal that ranged from 11 to 15 days. These rates indicated rapid enhanced biodegradation, which is promising for application of a full- scale bioremediation remedy. Ultimately, however, the mass of residual or sorbed TCE in the aquifer that remains accessible for dissolution and biodegradation would likely control the time required for a full- scale bioremediation effort to achieve performance goals for TCE and cisDCE specified in the Record of Decision for Site K.

The field pilot tests showed that the relatively low hydraulic head gradients and temporal changes in groundwater flow directions in the shallow aquifer would add complexity to a full- scale bioremediation effort. The radius of influence (ROI) at GS1 and GS3 (16.3 ft and 12.7 ft, respectively) were close to the design ROI of 15 ft. The estimated ROI at GS2 was about four times the design ROI, but may be less reliable at this location owing to groundwater flow direction. In addition, the low temperatures following WBC- 2 injection in late November to early December 2021, in combination with the low hydraulic head gradients, were probably major factors in the delay observed before the onset of enhanced biodegradation following injection of the culture. Additional test injections could be beneficial to optimize the timing of donor and culture injections with the variable temperatures and hydraulic head in the shallow aquifer.

Channel change and sediment transport in the Puyallup River watershed through 2022

Anderson, Scott — Tue, 31 Mar 2026 13:38:43

The Puyallup River drains a 990 square mile watershed in western Washington, with headwaters on the glacier-covered flanks of Mount Rainier. Major tributaries include the White, Carbon, and Mowich Rivers. In the levee-confined reaches of the lower watershed, loss of flood conveyance due to sand and gravel deposition has been a chronic issue. Over much of the 20th century, flood conveyance was maintained through sediment removal, but this practice ended in the late 1990s. Flood hazard management activities since the 1990s have primarily involved levee removal or setback projects. Assessments of 1984-2009 repeat cross sections suggested that sediment deposition rates were particularly high in reaches with recent levee setbacks. However, there have been no assessments of recent deposition rates since the 2009 surveys. There are also concerns that intensifying flood hydrology or increased sediment delivery from Mount Rainier may exacerbate deposition. However, assessment of those risks has been hindered by limited understanding of watershed-scale sediment delivery and routing, particularly for coarse sand and gravel.

The U.S. Geological Survey, in cooperation with Pierce County, initiated this study to improve understanding of sediment deposition in the lower Puyallup River watershed. This work is primarily based on differencing of multiple aerial lidar datasets collected during 2002–2022, supplemented by early 1990 photogrammetric elevation datasets, geomorphic assessments of streamgage data, historical topographic surveys from 1907, and previously collected sediment transport measurements. Analyses cover the Puyallup, Carbon, and Mowich Rivers, but do not include the White River.

During 2004–2020, repeat aerial lidar indicates that 1.3 ± 0.3 million yd3 of sediment accumulated in the lower 20 valley miles (VMs) of the Puyallup River, averaging 80,000 ± 20,000 cubic yards per year (yd3/yr). Deposition was observed during both 2004–11 and 2011–20 lidar differencing intervals. This continued a long-term depositional trend that extends back to at least 1977. From 2004 to 2011, deposition rates along the Soldiers Home levee setback reach, the only setback project downstream of VM 20 completed prior to 2011, were approximately four times higher than in adjacent unmodified reaches. From 2011 to 2020, two additional setback projects were completed; volumetric deposition rates over all three setback reaches were similar to adjacent unmodified reaches, suggesting elevated setback deposition in the 2004–11 interval may have been influenced by an extreme flood in November 2006. These levee setback projects increased the local cross-sectional area of the floodway, used as a rough proxy for relative flood conveyance, by 50 to 200 percent above 2004 conditions. If deposition continued at recent rates, cross-sectional area over the levee setback reaches would be reduced back to 2004 values by 2050-90.

Deposition also occurred over the lower six VMs of the Carbon River during 2004–20, though volumes (0.15 ± 0.09 million yd3) were an order of magnitude lower than along the Puyallup River. Relatively lower deposition rates in the Carbon River are most likely the combined result of modestly lower incoming sediment loads, modestly steeper channel slope, and the additional sediment transport capacity provided by two large non-glacial tributaries that enter the Carbon River near VM 5.

Upstream of the depositional reaches described above, 2002–22 sediment storage trends along the Puyallup, Carbon, and Mowich Rivers were predominately negative (net erosion) up to the Mount Rainier National Park boundary. Net erosion was the result of bank and bluff erosion exceeding deposition across wetted channel and bare gravel areas, as opposed to uniform vertical downcutting. Net erosion along these river valleys delivered 3.4 ± 0.6 million yd3 to the river system, equivalent to 190,000 ± 35,000 yd3/yr. Most of that volume was supplied by erosion of relatively low (4–10 ft) surfaces along the Puyallup and Mowich Rivers and tall (300 ft) glacial bluffs along the lower Carbon River. Substantial aggradation from 1984 to 2009 reported by Czuba and others (2010) along reaches of the Puyallup River (VM 19–22) where levee confinement has recently been removed was most likely an artifact of methodologic bias.

The Puyallup, Mowich, and Carbon Rivers drain five distinct glaciated watersheds on the flanks of Mount Rainier, four of which were assessed in this study. All four watersheds were impacted by an extreme November 2006 rainstorm. Between 2002 and 2008, debris flows occurred in all four headwater areas, collectively eroding at least 2.1 million yd3 of sediment. These debris flows formed distinct deposits one to two miles downstream of source areas, depositing 30-50 percent of the material eroded upstream. From 2008 to 2022, no headwater debris flows were observed and overall rates of geomorphic change in the headwaters were low. Rivers eroded into debris flow deposits emplaced over the 2002–08 interval, but re-deposited equivalent volumes of material within a half mile downstream.

Stage-discharge relations at five streamgages on upland rivers draining Mount Rainier show either net channel incision or dynamic variability with no long-term trend over the past 60–100 years. Observations of pervasive river valley erosion and stable or incising trends at long-term streamgages in the upper watershed do not support prior claims of widespread and accelerating aggradation of upland rivers draining Mount Rainier.

Erosion and deposition volumes estimated in this report were combined with sediment transport estimates from limited suspended sediment and bedload measurements, estimates of sub-glacial erosion rates, and sediment delivery from non-glacial tributaries to construct watershed-scale sediment budgets for the Puyallup River watershed. During 2004–20, the estimated sediment load entering the depositional lowlands was well balanced by estimated inputs from, in order of relative magnitude, subglacial erosion (33–60 percent of total sediment load), erosion along the major river valleys (25–45 percent), erosion in recently deglaciated headwater areas (7–17 percent) and non-glacial tributaries (3–9 percent). These results are specific to the study period and represent total sediment loads, most of which is fine material carried in suspension. The relative sourcing of sand and gravel may be different than implied by this sediment budget.

Downstream of VM 12, comparison of 1907 and 2009 channel surveys show net lowering of the channel thalweg of 4–12 ft. A long-term gage near VM 22 shows lowering of 4–5 ft through the 1960s. Lowering at both locations was inferred to be a channel response to the substantial straightening, and so steepening, of the river during major phases of levee construction through the early and mid-20th century.

Application of a simple empirical bedload-discharge power-law relation to an ensemble of model-estimated daily mean discharge records in the lower Puyallup River between 1977 and 2100 projects that annual bedload transport capacity in the lower Puyallup River will increase by 20–60 percent by the middle of the 21st century. Actual changes in bedload transport and deposition rates will depend on concurrent changes in sediment supply and local hydraulics governing deposition.

This report presents several key conclusions. First, the persistence and spatial patterns of sand and gravel deposition along the lower Puyallup River support prior claims that deposition is fundamentally caused by decreases in channel slope moving downstream. Given this underlying cause and the abundance of sand and gravel available to be transported downstream, deposition is likely to continue for the foreseeable future. Second, despite continued sediment deposition, recent levee setback projects in the lower Puyallup River will likely provide several decades of flood conveyance benefits relative to a no-action alternative. Third, while the rivers linking Mount Rainier to the Puget Sound lowlands have often been discussed as conduits that either pass or accumulate sediment from Mount Rainier, observations from 2002–22 show these river valleys acting as substantial sediment sources, delivering three times more sediment than recently deglaciated headwater areas on Mount Rainier. While the persistence and underlying cause of recent river valley erosion remain unknown, sediment storage dynamics along these river valleys are likely to be a major control on sand and gravel delivery to the lower watershed.

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River near Kansas City, Missouri, August 8–9, 2023

Huizinga, Richard; Rivers, Benjamin — Mon, 23 Feb 2026 14:47:38

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 8 bridge crossings of the Missouri River near Kansas City, Missouri, on August 8–9, 2023. A multibeam echosounder mapping system was used to obtain channel- bed elevations for river reaches that extended about 1,550 to 1,640 feet longitudinally and generally extended laterally across the active channel from bank to bank during low floodflow to nonflood conditions. These surveys provided the channel geometry and hydraulic conditions of the river at the time of the surveys and provided characteristics of scour holes, which may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. The data collected from the surveys may also be useful to the Missouri Department of Transportation as a record of low floodflow conditions in regards to the stability and integrity of the bridges with respect to bridge scour. Bathymetric data were collected around every in- channel pier. Scour holes were at most piers where bathymetry could be obtained, except for those piers on banks or surrounded by riprap. All the bridge sites in this study were surveyed and documented in previous studies.

The average difference between the bathymetric surfaces ranged from 0.07 to 4.16 feet higher in 2023 than 2019, which indicates overall deposition between the survey dates, as might be expected based purely on streamflow at the time of the survey. However, the average difference between the bathymetric surfaces ranged from 1.44 feet higher to 1.88 feet lower in 2023 than 2015, which indicates a dynamic equilibrium of scour and deposition overall between those surveys, despite the lower flow conditions in 2023. Similarly, the average difference between the bathymetric surfaces ranged from 3.18 feet higher to 5.19 feet lower in 2023 than 2011, which indicates a relative equilibrium between scour and deposition overall, albeit the trend was toward scour as might be expected because of the substantial flood event in 2011.

Riprap blankets and alignment to flow had a substantial effect on the size of the scour hole for a given pier. Piers that were partially or fully surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers with no rock or riprap and effectively mitigated the scour holes historically observed at these piers. Several of the structures had piers that were skewed to primary approach flow. At most of the structures, the scour hole was deeper and longer on the side of the pier with impinging flow than the leeward side, with some amount of deposition on the leeward side, as typically observed at piers skewed to approach flow.