Geology ReportsSearch

USGS · 70274549

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

Abstract

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.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Scott W. Anderson. 2026-02-18. Channel change and sediment transport in the Puyallup River watershed through 2022. https://doi.org/10.31223/x5hr0n

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Investigation of drivers and sources of potential harmful algal blooms in Lake Anna, Virginia, 2023-24

Algal blooms in freshwater reservoirs are an increasing concern because of the risk to aquatic ecosystem health, recreational use, and water quality. Among the taxa that comprise an algal bloom, cyanobacteria are of particular concern due to the potential to produce toxic compounds which can cause acute and chronic illness in humans if ingested through contaminated water or shellfish. When toxins are above a toxin-specific threshold in the water column, the Commonwealth of Virginia categorizes the bloom as a harmful algal bloom (HAB), distinguishing it from an algal bloom without toxins present. Toxins are of particular concern in bodies of water which serve as major recreational areas, such as Lake Anna in central Virginia. Lake Anna is a 27 km-long and 47 km2 reservoir in central Virginia. Since the Commonwealth of Virginia revised the HAB monitoring framework in 2018, algal bloom advisories for potentially toxigenic blooms have been issued in the western, riverine zone of the reservoir. There were detections of microcystin and anatoxin production but no detections of cylindrospermopsin and saxitoxin production in Lake Anna; however, toxin concentrations were low, at values less than 0.65 µg/L. In 2023, the U.S. Geological Survey, in cooperation with Virginia Department of Environmental Quality, initiated a 19-month intensive and extensive study of western Lake Anna and its tributaries (North Anna River and Pamunkey Creek) to (1) characterize the algal community, biomass, and toxin production; (2) identify potential drivers leading to algal bloom initiation, persistence, and decline; and (3) assess the primary sources of the drivers contributing to algal blooms. Continuous and discrete monitoring of water quality within the lake, watershed inputs, meteorological conditions, and algal communities were monitored. Multiple modeling scenarios identified that the formation and persistence of algal blooms in Lake Anna are driven by interactions among chemical, macro- and micronutrient, and physical factors. Key drivers include total nitrogen, total phosphorus, water temperature, suspended sediment concentration, wind speed, copper, iron, molybdenum, nickel, sodium, and alkalinity. During this study, winter stormflow events delivered the highest loading of macro- and micronutrients from the tributaries to the lake. The lake is a sink for these nutrients in the water column or in the lake-bed sediments. During high-occupancy and recreation periods on the lake in the summer months, boating activity can promote internal loading by the resuspension of the lake-bed sediments, thereby increasing the nutrient availability for algal uptake. Thermal stratification and hypoxic conditions were observed, which can contribute to the passive release of nutrients into the water column, further facilitating algal bloom proliferation. These processes may create a positive feedback loop that supports the formation, persistence, and annual recurrence of algal blooms in the western region of Lake Anna. Continued water-quality monitoring of early warning indicators may improve the detection and forecasting of algal growth in the lake and help managers proactively manage water quality and algal growth.

EarthArXiv

An automatable method for developing preliminary stream crossing designs using lidar and hydraulic modeling in Massachusetts

Upgrading aging and undersized stream crossings (places where culverts and bridges intersect a stream) with replacements designed for aquatic organism passage can improve infrastructure resiliency, geomorphic stability, and maintenance costs. The goal of this work is to develop preliminary culvert designs that convey hydraulic design floods and meet the Massachusetts Stream Crossings Standards for select existing stream crossing sites across Massachusetts. We derive stream and infrastructure data from lidar and other geospatial datasets as input for one-dimensional hydraulic models. Culvert dimensions are iteratively increased until the structures can convey the hydraulic design flood flows with adequate freeboard. The approach was evaluated by comparing lidar-derived designs to designs produced with surveyed data at 30 stream crossings across Massachusetts. The automatable methodology was then applied to 343 stream crossing sites in the Housatonic River watershed in southwestern Massachusetts. For each stream crossing site, preliminary designs were developed to meet the Massachusetts Stream Crossing Standards and convey the 10-, 4-, 2-, and 1-percent annual exceedance probability flood flows. The U.S. Geological Survey StreamStats web application hosts the 3-sided box, 3-sided arch, and pipe culvert preliminary designs. Provided supporting information includes site information, habitat quality assessments, and modeled water surface elevations. The publicly available preliminary designs can be used by municipalities, engineers, and others to evaluate design options, prioritize replacements, and analyze ecological benefits. These designs also may help municipalities pursue grant opportunities and estimate replacement costs prior to investing in field assessments and engineering design.

EarthArXiv

Ecosystem metabolism as an early warning indicator of lake algal blooms

Algal blooms represent ecosystem state shifts that degrade drinking water, restrict recreation, threaten public health, and lower property values. Detecting blooms in advance on management relevant timescales of days to weeks can support proactive intervention. Early warning statistics derived from indicator time series offer a framework for detecting state shifts, but the use of lake metabolic estimates (gross primary production, ecosystem respiration, and net ecosystem production) in this context has been limited. Here, we modeled lake metabolism using water-quality monitoring data from eight U.S. lakes spanning a trophic gradient to test the hypothesis that metabolic estimates serve as effective early warning indicators of algal blooms. Time series of early warning statistics were analyzed using the conventional Kendall’s tau approach and a complementary threshold-based approach. The threshold-based approach yielded higher bloom detection rates, longer early warning lead times, and fewer false alarms than the Kendall’s tau approach. Across lakes, metabolic metrics detected 86% of blooms, with the strongest early warning signals in oligotrophic, eutrophic, and hypereutrophic systems. However, directly measured water-quality parameters, particularly water temperature and chlorophyll concentration, were the primary early warning indicators of blooms. Overall, these results demonstrate that metabolic estimates can provide early warning signals of bloom development but are secondary to directly measured water-quality parameters, indicating a complementary rather than primary role in early bloom detection.

EarthArXiv