Geology topics

Scott W. Anderson

Publications and source records attributed to Scott W. Anderson.

At least 19 recordsLinked to original sources

USGS2026

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

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.

EarthArXiv↗

USGS2024

Quantifying fine sediment infiltration in spawning gravel used by Chinook salmon (Oncorhynchus tshawytscha) in the Sauk River Basin, Washington, 2018–21

Fine sediment can infiltrate into river substrate that salmonid fish species ( Oncorhynchus spp.) use to spawn. High levels of sediment infiltration can increase egg-to-fry mortality, which corresponds to the period when salmonids are still residing in the subsurface gravels. This study quantifies fine sediment infiltration of Chinook salmon ( Oncorhynchus tshawytscha ) spawning habitat during the egg-to-fry emergence period over three years in the Sauk River, which has naturally high fine sediment loads and important native salmon populations. Additionally, this study qualitatively assesses how grain size distribution of the riverbed and adjacent gravel bars compare to grain size distribution following fine sediment infiltration to evaluate if riverbed or gravel bar grain size distributions may provide information on the potential for fine sediment infiltration in spawning gravels. Fine sediment infiltration into spawning gravels was quantified using sediment boxes and infiltration bags that were installed in artificial redds constructed at known Chinook salmon spawning locations at three study sites on the Sauk River over the expected egg-to-fry period. Over the three-year study period (August 2018–April 2021), fraction finer of sediment (grain sizes of less than two millimeters), ranged from 0.12 to 0.23 across the three study sites and years. Based on a comparison of field observations from this study and percentage egg-to-fry survival curves found in the literature, the expected survival for Chinook salmon eggs in the Sauk River is roughly 30 percent. Expected survival increases to approximately 90 percent if eggs are eyed and thus farther along in their development. Our field study did not evaluate the progression of infiltration, so it is unknown if observed fine sediment infiltration was at this relatively high rate during the period that corresponded to early egg development, when eggs are more sensitive to fine sediment infiltration. Dissolved oxygen in the gravels is largely above critically low levels (4 milligrams per liter) during sensitive periods corresponding to egg development and is interpreted not to affect egg-to-fry mortality. Active channel morphology in the middle reaches of the Sauk River may pose an additional challenge to pre-emergence survival. Channel change, deposition, and potential scour at the middle Sauk River study site likely contributed to low recovery rates of both sediment boxes and infiltration bags in two of the three study years. In terms of grain size distributions of riverbed sediment and adjacent gravel bars representing the potential for fine sediment infiltration into spawning gravels, both riverbed and gravel bar bulk subsurface sediment samples had higher fraction finer for the representative fine grain size of two millimeters compared to the sediment boxes and infiltration bags. Therefore, riverbed and gravel bar samples may serve as a conservative first order proxy for potential fine sediment infiltration into spawning gravels, with the understanding that these samples may overestimate fine sediment infiltration by up to 15 percent.

Washington↗

USGS2024

River channel response to the removal of the Pilchuck River Diversion Dam, Washington State

In August 2020, the 3-m tall Pilchuck River Diversion Dam was removed from the Pilchuck River, allowing free fish passage to the upper third of the watershed for the first time in over a century. The narrow, 300-m long impoundment behind the dam was estimated to hold 4,000–7,500 m 3 of sand and gravel, representing less than one year's typical bedload flux. Repeat cross section surveys, stage sensors, and time-lapse cameras were used to document the physical channel response over the first year following dam removal. A total of 7,400 m 3 (effectively 100%) of impoundment sediment was eroded in the first year, with 25% accomplished by manual regrading during dam removal. Most river-caused erosion occurred during a sequence of modest flows in October 2020. Downstream deposition totaled 4,300 m 3 , predominately filling in the first 100 m downstream of the dam site. Deposition tapered below detectable levels within 350 m, and most downstream channel adjustments occurred before November 2020. Multiple high flows after December 2020 caused little upstream or downstream change. The physical river response to this dam removal then appears to have been largely accomplished within several months by modest flows, consistent with pre-removal modeling and observations from other regional dam removals. Efficient sediment evacuation was likely aided by the narrow and steep-walled impoundment geometry. Our observations support existing guidance that the physical river response to small dam removals is typically rapid and modest; the benefits of removal may then often be gained with minimal negative downstream geomorphic impacts.

Washington↗

USGS2023

A general approach for evaluating of the coverage, resolution, and representation of streamflow monitoring networks

Streamflow monitoring networks provide information for a wide range of public interests in river and streams. A general approach to evaluate monitoring for different interests is developed to support network planning and design. The approach defines three theoretically distinct information metrics (coverage, resolution, and representation) based on the spatial distribution of a variable of interest. Coverage is the fraction of information that a network can provide about a variable when some areas are not monitored. Resolution is the information available from the network relative to the maximum information possible given the number of sites in the network. Representation is the information that a network provides about a benchmark distribution of a variable. Information is defined using Shannon entropy where the spatial discretization of a variable among spatial elements of a landscape or sites in a network indicates the uncertainty in the spatial distribution of the variable. This approach supports the design of networks for monitoring of variables with heterogeneous spatial distributions (“hot spots” and patches) that might otherwise be unmonitored because they occupy insignificant portions of the landscape. Areas where monitoring will maintain or improve the metrics serve as objective priorities for public interests in network design. The approach is demonstrated for the streamflow monitoring network operated by the United States Geological Survey during water year 2020 indicating gaps in the coverage of coastal rivers and the resolution of low flows.

Colorado↗

USGS2023

Changes in suspended-sediment yields under divergent land-cover disturbance histories: A comparison of two large watersheds, Olympic Mountains, USA

Improvements in timber harvest practices and reductions in harvest volumes over the past half century are commonly presumed to have reduced sediment loads in many western US rivers. However, direct assessments in larger watersheds are relatively sparse. Here, we compare 2019–21 sediment concentrations against those of the late 1970s in the Bogachiel and Calawah River watersheds, adjacent and similarly sized (~300 km 2 ) basins in the western Olympic Mountains of Washington State. The Calawah River watershed has experienced significant land-cover disturbance, including a large 1951 fire, extensive post-fire salvage logging, and relatively high rates of timber harvest through the 1990s. In contrast, the Bogachiel River watershed did not burn, and experienced only modest timber harvest that largely post-dated 1970s sediment monitoring. Channel-width trends suggest the Calawah River was still recovering from 1950s disturbances in the late 1970s. We found that 2019–21 suspended-sediment loads in the Calawah River were 2.3–2.6 times lower than would have been expected based on 1970s sediment rating curves, while recent loads in the Bogachiel River were a factor of 1.4 ± 1.0 lower. We consider the plausibility and possible explanations of declining concentrations in the less-disturbed Bogachiel River. Suspended-sediment yields in the Bogachiel River were two times higher than yields in the Calawah River, which is attributed to a combination of modestly higher precipitation, more efficient runoff generation, and more extensive and erodible Quaternary valley fills in the Bogachiel River. Regional shifts in flood hydrology have also influenced suspended-sediment loads in both watersheds. Our results then document a significant decline in suspended-sediment concentrations in the Calawah River over the past half century. Reduced land-cover disturbance provides the simplest and most likely explanation for this decline, though the wide range of possible concentration changes in the Bogachiel River leaves open possibilities that other processes (human, natural, or methodologic) could be a factor.

Washington↗

USGS2023

Multi-decadal erosion rates from glacierized watersheds on Mount Baker, Washington, USA, reveal topographic, climatic, and lithologic controls on sediment yields

Understanding land surface change in and sediment export out of proglacial landscapes is critical for understanding geohazard and flood risks over engineering timescales and characterizing landscape evolution over geomorphic timescales. We used automated Structure from Motion software to process historical aerial photographs and, with modern lidar data, generated a high-resolution DEM time series with coverage over 10 glacierized watersheds on Mount Baker, Washington, USA for the time period between 1947 and 2015. We measured basin-wide sediment yields and sediment redistribution on hillslopes and in stream channels. Slopes within most measured erosion sites are above theoretical and observed debris-flow thresholds. We observed significant erosion of hillslopes and limited deposition on hillslopes and in stream channels. Sediment delivery ratios during time periods with net erosion averaged 0.73. We determined, consistent with previous field observations, that debris flows originating from moraines are a primary erosion mechanism in proglacial zones on Mount Baker. Time series measurements indicate that temporal variability in erosion rates is associated with climate oscillations, with higher erosion rates during cooler-wetter periods. Basin-wide sediment yield is positively correlated with lithology (r 2 = 0.54), hillslope angle (r 2 = 0.52), drainage area (r 2 = 0.82), and negatively correlated with stream channel slope (r 2 = 0.67). Topographic differences between high and low yielding basins indicate that spatial variability in erosion on Mount Baker is sensitive to Pleistocene and Holocene glacial and volcanic activity. Specific sediment yields in six basins averaged 4600 ton/km 2 /yr, consistent with global measurements in glacierized catchments. Specific sediment yield decreased with increasing basin area, with total loads in the downstream main stem Nooksack River estimated between 480 and 820 ton/km 2 /yr. Proglacial sediment yields account for between 18 and 32 % of total sediment load in the main stem Nooksack River and exceed contributions by bluff and terrace erosion, which account for between 8 and 13 % of total load. Our findings indicate that erosion in glacierized basins is sensitive to decadal climate oscillations and that high proglacial sediment yields provide an important contribution to river systems downstream, particularly in catchments where upland topography and lithology is favorable.

Washington↗

USGS2022

Relative contributions of suspended sediment between the upper Suiattle River Basin and a non-glacial tributary, Washington, May 2016–September 2017

Concentrations of suspended sediment were measured in discrete samples and turbidity was continuously monitored at four U.S. Geological Survey streamgages in western Washington State, including one gage on the Sauk River; two gages on the Suiattle River, a tributary to the Sauk River; and one gage on Downey Creek, a tributary to the Suiattle River. The Suiattle River is a sediment-rich stream with headwaters on Glacier Peak, a glaciated volcano in the northern Cascade Range. Contributions of suspended sediment to the Suiattle River from unglaciated tributaries, represented by Downey Creek, were compared to the contributions from Glacier Peak in the upper Suiattle River watershed. During summer 2017, a period for which complete records of discharge and sediment data were available for all three streamgages in the Suiattle River Basin, the suspended-sediment load from Downey Creek (drainage area [DA] 93 square kilometers [km 2 ]) was 1,400 metric tons, which is equivalent to a sediment yield of about 15 metric tons per km 2 . During the same period, the suspended-sediment load from the upper Suiattle River (DA 176 km 2 ) was 142,000 metric tons, or a sediment yield of about 800 metric tons per km 2 ; and the suspended-sediment load from the lower Suiattle River (DA 733 km 2 ) was 230,000 metric tons, or a sediment yield of about 300 metric tons per km 2 . The Downey Creek Basin accounts for 13 percent of the drainage area of the Suiattle River watershed but contributed only 0.6 percent of the suspended-sediment load over the summer of 2017 and water year 2017 (October 1, 2016–September 30, 2017). In contrast, the upper Suiattle River Basin, which accounts for 24 percent of the entire Suiattle River watershed, contributed 62 percent of the suspended-sediment load during the summer of 2017. Given the short period for which data were collected, it cannot be known with certainty whether the above values are representative of long-term means. The relatively minor contribution of suspended sediment from Downey Creek, however, is consistent with the expectation that the upper Suiattle River, which drains Glacier Peak, is the dominant contemporary source of suspended sediment to the Sauk River. During summer 2016, the suspended-sediment load in the upper Siuattle River (180,000 metric tons) was more than double the estimated load in the lower Sauk River (80,000 metric tons), even though the upper Suiattle River represents only 10 percent of the total contributing area to the lower Sauk River Basin. This ratio of relative contribution is interpreted as an indication of transient storage of sediment along the Suiattle and Sauk Rivers between the two streamgaging stations. In the glaciated upper Suiattle River Basin, sediment is transported by annual glacial-melt processes in spring and summer months, deposited during the summer base-flow period, and then remobilized by fall and winter floods for delivery to the lower Sauk River.

Washington↗

USGS2022

Bedload and suspended-sediment transport in lower Vance Creek, western Washington, water years 2019–20

Vance Creek drains a 24 square mile area of the Olympic Mountains in western Washington. The lower 4 miles of the creek often go dry in discontinuous patches during the summer, limiting salmon rearing success. To better understand sediment transport dynamics in the creek and aid in potential restoration design, bedload and suspended-sediment concentration samples were collected for water years 2019–20 at a site about 2 miles upstream from the creek’s confluence with the South Fork Skokomish River. Fifty bedload samples and 7 suspended-sediment concentration samples were collected over 7 sampling days. These samples were used to develop rating curves relating bedload flux or suspended-sediment concentration to discharge. Mean annual bedload flux was estimated to be 12,200 ± 2,300 tons per year for water years 1930–2020 period of record, based on application of the derived bedload rating curve to an extrapolated daily discharge record. The mean annual suspended-sediment load over the same period was estimated to be 9,000 tons per year with large, but unquantified, uncertainty. Bedload material was predominantly gravel from 0.08 to 2.5 inches (2 to 64 millimeters) in diameter. At the highest sampled discharges, approximately equivalent to a 50 percent annual exceedance probability (2-year peak-flow event), the bedload grain-size distribution was similar to that of the local channel bed. Bedload grain-size distributions generally coarsened as discharge increased. The suspended-sediment load was consistently one-half sand and one-half silt and clay, regardless of discharge. Bedload constituted about 60 percent of the total sediment flux (bedload plus suspended load). This is near the upper limit of values observed in a global compilation of long-term load partitioning data. Sediment transport at the Vance Creek sampling site was compared with sediment-transport data from five other watersheds in the region. To facilitate comparisons, mean annual loads were divided by mean annual runoff volume to obtain an effective average sediment concentration. This normalization accounts for differences in both drainage area and mean runoff depth between the comparison watersheds. At the three comparison watershed sites with relatively complete sediment-transport data, mean bedload concentrations ranged from 44 to 109 milligrams per liter (mg/L) and mean suspended-sediment concentrations ranged from 139 to 374 mg/L; bedload constituted 21 to 29 percent of the total sediment load. The mean bedload concentration at the Vance Creek sampling site (69 mg/L) fell in the middle of the range observed in comparison watersheds, whereas the mean suspended-sediment concentration (50 mg/L) was markedly lower. Bedload samples at the Vance Creek sampling site also were generally less sand rich (sample-average sand fraction was 13 percent at Vance Creek versus 20 to 37 percent for comparison waters). Bedload transport rates at the Vance Creek sampling site appear relatively average for the region, given the drainage basin area and average runoff. In contrast, the supply and transport of finer material, both in the suspended load and the sand fraction of the bedload, are relatively low.

Washington↗

USGS2021

Spatial and temporal controls on proglacial erosion rates: A comparison of four basins on Mount Rainier, 1960 to 2017

The retreat of alpine glaciers since the mid-19th century has triggered rapid landscape adjustments in many headwater basins. However, the degree to which decadal-scale glacier retreat is associated with systematic or substantial changes in overall coarse sediment export, with the potential to impact downstream river dynamics, remains poorly understood. Here, we use repeat topographic surveys to assess geomorphic change in four partly-glaciated basins on a stratovolcano (Mount Rainier) in Washington State at roughly decadal intervals from 1960 to 2017. The proglacial extents of the four basins show distinct geomorphic trajectories, ranging from substantial and sustained net erosion to relatively inactive with net deposition. These different trajectories correspond to differences in initial (1960) valley floor gradients, and can be effectively understood as valley floor grade adjustments. Significant erosion was most often accomplished by debris flows triggered by extreme rainfall or glacial outburst floods, though a single rockfall mobilized more material than all other events combined. Year-to-year runoff events had little measurable geomorphic impact. Exported material tended to accumulate in broad deposits within several kilometers of source areas and largely remained there through the end of the study period. Over 10- to 100-year timescales, we nd that the impact of glacier retreat on coarse sediment yield may then vary substantially according to the geometry and storage trends of the newly-exposed valley floor; the timing of that response may also be dictated, and potentially obscured, by stochastic and/or extreme events.

Washington↗

USGS2021

Channel response to a dam‐removal sediment pulse captured at high‐temporal resolution using routine gage data

In this study, we captured how a river channel responds to a sediment pulse originating from a dam removal using multiple lines of evidence derived from streamflow gages along the Patapsco River, Maryland, USA. Gages captured characteristics of the sediment pulse, including travel times of its leading edge (~7.8 km yr −1 ) and peak (~2.6 km yr −1 ) and suggest both translation and increasing dispersion. The pulse also changed local hydraulics and energy conditions, increasing flow velocities and Froude number, due to bed fining, homogenization and/or slope adjustment. Immediately downstream of the dam, recovery to pre‐pulse conditions occurred within the year, but farther downstream recovery was slower, with the tail of the sediment pulse working through the lower river by the end of the study 7 years later. The patterns and timing of channel change associated with the sediment pulse were not driven by large flow or suspended sediment‐transporting events, with change mostly occurring during lower flows. This suggests pulse mobility was controlled by process‐factors largely independent of high flow. In contrast, persistent changes occurred to out‐of‐channel flooding dynamics. Stage associated with flooding increased during the arrival of the sediment pulse, 1 to 2 years after dam removal, suggesting persistent sediment deposition at the channel margins and nearby floodplain. This resulted in National Weather Service‐indicated flood stages being attained by 3–43% smaller discharges compared to earlier in the study period. This study captured a two‐signal response from the sediment pulse: (1) short‐ to medium‐term (weeks to months) translation and dispersion within the channel, resulting in aggradation and recovery of bed elevations and changing local hydraulics; and (2) dispersion and persistent longer‐term (years) effects of sediment deposition on overbank surfaces. This study further demonstrated the utility of US Geological Survey gage data to quantify geomorphic change, increase temporal resolution, and provide insights into trajectories of change over varying spatial and temporal scales.

Maryland↗

USGS2020

Coarse sediment dynamics in a large glaciated river system: Holocene history and storage dynamics dictate contemporary climate sensitivity

The gravel-bedded White River drains a 1279 km 2 basin in Washington State, with lowlands sculpted by continental glaciation and headwaters on an actively glaciated stratovolcano. Chronic aggradation along an alluvial fan near the river’s mouth has progressively reduced flood conveyance. In order to better understand how forecasted climate change may influence coarse sediment delivery and aggradation rates in this lowland depositional setting, we assessed the contemporary delivery and routing of coarse sediment through the watershed; this assessment was based on a rich set of topographic, sedimentologic, and hydrologic data from the past century, with a focus on repeat high-resolution topographic surveys from the past decade. We found that most of the lower river’s contemporary bed-load flux originates from persistent erosion of alluvial deposits in the lower watershed. This erosion is a response to a drop in local base level caused by a major avulsion across the fan in 1906 and then augmented by subsequent dredging. The 1906 avulsion and modern disequilibrium valley profiles reflect landscape conditioning by continental glaciation and a massive mid-Holocene lahar. In the proglacial headwaters, infrequent large sediment pulses have accomplished most of the observed coarse sediment export, with exported material blanketing downstream valley floors; during typical floods, transported bed material is largely sourced from erosion of these valley floor deposits. Throughout the watershed, we observe decadal-scale coarse sediment dynamics strongly related to the filling or emptying of valley-scale sediment storage over 10 2 −10 4 yr time scales, often in response to major disturbances that either emplace large deposits or influence their redistribution. Paraglacial responses in large watersheds are suggested to be inherently complicated and punctuated as a result of internal landform interactions and stochastic/threshold-dependent events. We argue, in combination, that Holocene disturbance, storage dynamics, and human flow modification make coarse sediment fluxes in the lower White River relatively insensitive to decadal climate variability. Results highlight the degree to which river sensitivity to contemporary disturbance, climatic or otherwise, may be contingent on local and idiosyncratic watershed histories, underscoring the need to unpack those histories while demonstrating the utility of watershed-scale high-resolution topography toward that end.

Washington↗

USGS2019

Sediment storage and transport in the Nooksack River basin, northwestern Washington, 2006–15

The Nooksack River is a dynamic gravel-bedded river in northwestern Washington, draining off Mount Baker and the North Cascades into Puget Sound. Working in cooperation with the Whatcom County Flood Control Zone District, the U.S. Geological Survey studied topographic, hydrologic, and climatic data for the Nooksack River basin to document recent changes in sediment storage, long-term bed elevation trends, rates of sediment transport, and factors influencing surficial drainage in order to support ongoing river management. Differences in elevations between topographic and bathymetric surveys in 2005/06 and 2013/15 indicate the active channel aggraded about 1–2 feet locally near the cities of Ferndale and Everson but was primarily stable between them. The active channel upstream of Nugent’s Corner generally incised. Total incision upstream of Nugent’s Corner to Glacier Creek generated 2.3 ± 1.7 million cubic yards of sediment from 2005/06 to 2013 and likely represented a significant source of coarse sediment to the lower mainstem river over that time. Long-term records of local channel-bed elevation, derived from U.S. Geological Survey streamgage data, show bed-elevation changes of about 1–3 feet. The river bed at most streamgages exhibits long-term trends, with relatively consistent rates of change on the order of 1 foot per decade that persist years to decades. Lagged correlations in bed-elevation trends at all seven streamgages in the North Fork Nooksack and mainstem Nooksack suggest that decadal periods of persistent aggradation and incision originate in the North Fork and translate downstream. The channel-change signal propagates downstream 0.5–2.5 miles per year, with the rate of propagation scaling closely with channel slope. The pattern of incision and aggradation in the North Fork correlates with regional climate, where persistent incision follows extended cold and wet periods, and persistent aggradation follows extended warm and dry periods. Climate-driven variation in coarse-sediment delivery, primarily from the North Fork Nooksack, then appears to be a strong control on long-term vertical channel adjustments at sites downstream. The downstream-translating climate signal generated in the North Fork would account for recently observed aggradation at Everson and Ferndale but not the observed incision in unconfined reaches upstream of Nugent’s Corner from 2005–06 to 2013. This mismatch indicates that understanding how changes in sediment-supply influence those unconfined reaches remains a key uncertainty for predicting future channel change. Continuous turbidity monitoring integrated with suspended sediment and limited bedload sampling were used to calculate annual sediment loads at five sites in the basin. The sediment load in the lower river at Ferndale ranged from 0.78 to 1.17 million tons per year and averaged 0.97 million tons per year for WYs 2012–17. Suspended sediment made up 93 percent of the load, and bedload made up 7 percent. Most of the fine sediment load of the lower river is supplied from headwaters of the North, Middle, and South Fork Nooksack basins, with relatively little net increase in fine sediment loads in the lower mainstem basin. The three forks supply approximately equal proportions of the lower-river fine sediment load. However, the glacially sourced North and Middle Fork Nooksack basins carry a notably sandier suspended-sediment load than the South Fork Nooksack. A comparison of monthly streamflow and precipitation trends since 1981 indicate statistically significant increases in total spring precipitation and the number of spring days with measurable precipitation in much of the basin, as well as increases in mean spring river stage near Ferndale. Since no trends in mean spring discharge are observed, the trends in river stage are attributed primarily to observed changes in bed elevation. Changes in bed elevation and precipitation may then both have plausibly impacted field drainage in the lower river below Ferndale.

Washington↗

USGS2019

Downstream‐propagating channel responses to decadal‐scale climate variability in a glaciated river basin

Regional climate is an important control on the rate of coarse sediment mobilization and transport in alpine river systems. Changes in climate are then expected to cause a cascade of geomorphic responses, including adjustments in downstream channel morphology. However, the mechanics and sensitivity of channel response to short‐term climate variability remain poorly documented. In the Nooksack River, which drains a glaciated stratovolcano in Washington State, bed elevation changes were inferred from shifting stage–discharge relations at seven USGS stream gages. Decadal‐scale elevation trends at most sites can be explained as a downstream‐propagating channel response to regional climate variability, where periods of persistent warm, dry [cool, wet] conditions corresponded to periods of aggradation [incision]. The channel elevation response propagated downstream at a rate of one to four kilometers per year; propagation rate scaled closely with channel slope. Historical trends in glacier extent and flood intensity both show some potential to explain climate–sediment linkages, though assessing causation is complicated by the shared climate signal in both records. Results show the influence of the Pacific Decadal Oscillation, with relatively high coarse sediment yields prior to 1950 and since 1980, and notably lower sediment yields from 1950 to 1980. Measured sediment yields from nearby glaciated basins corroborate this history, suggesting a regional coherence to these climate–sediment linkages. These results document consistent relations between climate, sediment supply and downstream channel response at the basin‐scale, with channel responses propagating downstream over periods of decades with little apparent attenuation.

Washington↗

USGS2019

Uncertainty in quantitative analyses of topographic change: Error propagation and the role of thresholding

Topographic surveys inevitably contain error, introducing uncertainty into estimates of volumetric or mean change based on the differencing of repeated surveys. In the geomorphic community, uncertainty has often been framed as a problem of separating out real change from apparent change due purely to error, and addressed by removing measured change considered indistinguishable from random noise from analyses (thresholding). Thresholding is important when quantifying gross changes (i.e. total erosion or total deposition), which are systematically biased by random errors in stable parts of a landscape. However, net change estimates are not substantially influenced by those same random errors, and the use of thresholds results in inherently biased, and potentially misleading, estimates of net change and uncertainty. More generally, thresholding is unrelated to the important process of propagating uncertainty in order to place uncertainty bounds around final estimates. Error propagation methods for uncorrelated, correlated, and systematic errors are presented. Those equations demonstrate that uncertainties in modern net change analyses, as well as in gross change analyses using reasonable thresholds, are likely to be dominated by low‐magnitude but highly correlated or systematic errors, even after careful attempts to reduce those errors. In contrast, random errors with little to no correlation largely cancel to negligible levels when averaged or summed. Propagated uncertainty is then typically insensitive to the precision of individual measurements, and is instead defined by the relative mean error (accuracy) over the area of interest. Given that real‐world mean elevation changes in many landscape settings are often similar in magnitude to potential mean errors in repeat topographic analyses, reducing highly correlated or systematic errors will be central to obtaining accurate change estimates, while placing uncertainty bounds around those results provides essential context for their interpretation.

Earth Surface Processes and Landforms↗

USGS2017

Characterizing aquatic habitats for long‐term monitoring of a fourth‐order, regulated river in the Pacific Northwest, USA

A pragmatic approach to the long‐term monitoring of rivers leverages available information with targeted field investigations to address key uncertainties relevant to management decisions. An over‐arching management issue for many rivers is how reservoir operation affects the amount and location of in‐channel sediment and the resulting distribution of aquatic habitats. We integrate remotely acquired and field‐survey morphologic data for the Cedar River, Washington, to constitute the current status of aquatic habitats and benchmarks for long‐term monitoring that will inform streamflow management. Four key habitats (river edge, side channels, riffles, and pools) are feasible to monitor with high‐resolution aerial imagery, a longitudinal profile of the river, and a side channel inventory, but full characterization of the functional differences within these habitats requires additional information. Habitat use information such as redd surveys will continue to be important for long‐term monitoring where it cannot be inferred reliably from physical habitat characteristics.

Washington↗

USGS2017

Suspended sediment, turbidity, and stream water temperature in the Sauk River Basin, western Washington, water years 2012-16

The Sauk River is a federally designated Wild and Scenic River that drains a relatively undisturbed landscape along the western slope of the North Cascade Mountain Range, Washington, which includes the glaciated volcano, Glacier Peak. Naturally high sediment loads characteristic of basins draining volcanoes like Glacier Peak make the Sauk River a dominant contributor of sediment to the downstream main stem river, the Skagit River. Additionally, the Sauk River serves as important spawning and rearing habitat for several salmonid species in the greater Skagit River system. Because of the importance of sediment to morphology, flow-conveyance, and ecosystem condition, there is interest in understanding the magnitude and timing of suspended sediment and turbidity from the Sauk River system and its principal tributaries, the White Chuck and Suiattle Rivers, to the Skagit River. Suspended-sediment measurements, turbidity data, and water temperature data were collected at two U.S. Geological Survey streamgages in the upper and middle reaches of the Sauk River over a 4-year period extending from October 2011 to September 2015, and at a downstream location in the lower river for a 5-year period extending from October 2011 to September 2016. Over the collective 5-year study period, mean annual suspended-sediment loads at the three streamgages on the upper, middle, and lower Sauk River streamgages were 94,200 metric tons (t), 203,000 t, and 940,000 t streamgages, respectively. Fine (smaller than 0.0625 millimeter) total suspended-sediment load averaged 49 percent at the upper Sauk River streamgage, 42 percent at the middle Sauk River streamgage, and 34 percent at the lower Sauk River streamgage.

Washington↗

USGS2017

Geomorphic response of the North Fork Stillaguamish River to the State Route 530 landslide near Oso, Washington

On March 22, 2014, the State Route 530 Landslide near Oso, Washington mobilized 8 million cubic meters of unconsolidated Pleistocene material, creating a valley‑spanning deposit that fully impounded the North Fork Stillaguamish River. The river overtopped the 8-meter high debris impoundment within 25 hours and began steadily incising a new channel through the center of the deposit. Repeat topographic surveys, sediment transport measurements, bedload transport models, and observations of downstream channel change were used to document the establishment of that new channel through the landslide and assess the potential for downstream aggradation or channel change that might increase downstream flood hazards. Efficient erosion of the landslide deposit, associated with the steep knickzone formed by the downstream edge of the deposit, resulted in the re-establishment of a 20–40 meters wide, deeply inset channel through the entire deposit by May 2014, 2 months after the landslide. The mean water-surface elevation of the channel through the landslide decreased 7 meters during that 2-month period, and was about 1 meter above the pre-landslide profile in July 2014. The 2014–15 flood season, which included flows near the 0.5 annual exceedance probability discharge (2-year flood), widened the channel tens of meters, and further lowered the water-surface profile 0.5 meter. The planform position evolved slowly as a result of 5–20-meter high banks predominantly composed of clay-rich, cohesive lacustrine material. Erosion of the landslide deposit delivered a total of 820 thousand metric tons of sediment to the North Fork Stillaguamish River over the 18 months following the landslide. The sediment delivery from the deposit was predominantly fine grained: 77 percent (by mass) of the eroded material was silt or clay (less than 0.063 millimeter [mm]), 19 percent sand (0.063–2 mm), and 4 percent pebbles and cobbles (greater than 2 mm). Over the 18 months following the landslide, the bedload at a site 5 kilometers downstream of the landslide was estimated to be 310±65 thousand metric tons, and the suspended load at that same site was estimated to be 990±110 thousand metric tons. These loads represent the combined input from the landslide and ambient upstream sources; over the study interval, landslide sediment made up about 20–40 percent of the bedload, and 65–85 percent of the suspended-sediment load at this site. At a site 70 kilometers downstream of the landslide, near the mouth of the main‑stem Stillaguamish River, suspended sediment loads were estimated to be about 1,440 thousand metric tons, of which about 600 thousand metric tons, or 30 percent, likely was derived from the landslide. The mass of landslide sediment in suspension at the mouth of the river, and the timing of arrival of that sediment, indicates that about 70 percent of the landslide sediment eroded during the study period was quickly transported through the entire basin, exiting into Puget Sound within weeks of initial entrainment. Empirical bedload transport equations, in conjunction with surficial grain-size data and output from a one‑dimensional hydraulic model, were used to estimate spatial trends in bedload transport capacity, highlighting areas where reach-scale conditions would be most likely to promote deposition of coarse landslide sediment. Transport capacities decreased sharply over a reach about 5 kilometers downstream of the landslide and remained relatively low over the next 10 kilometers downstream. However, the magnitude of calculated transport capacities are large relative to the coarse sediment input from the landslide, suggesting that substantial deposition of landslide sediment was not likely to occur. These assessments were corroborated by observations of channel change, which indicated that the downstream channel response to the landslide was modest and short-lived. The most pronounced downstream effects included a wedge of aggradation just downstream of the landslide, about 1 meter high and extending a kilometer downstream, and a 0.3-meter pulse of aggradation observed 5 kilometers downstream of the landslide. In both locations, peak aggradation and channel response occurred within about a month of the landslide, and both sites had largely recovered to pre-landslide conditions by July 2014. No substantial channel change clearly linked to the landslide was observed after July 2014 except for a modest fining of surficial gravel size distributions and continued recovery and incision of the reach just downstream of the landslide. The muted downstream response of the North Fork Stillaguamish River to the State Route 530 Landslide primarily can be attributed to the cohesive, silt- and clay-rich material that bounded most of the new channel. Although the river efficiently incised a new channel through the deposit, subsequent rates of lateral erosion were slowed by the tall, cohesive banks, limiting the total volume of sediment delivery. Once entrained, however, most landslide material was rapidly transported downstream in suspension with little geomorphic effect. Landslide material coarse enough to travel as bedload was predominantly sand and fine gravel, and sediment transport models and observations of downstream change indicated that the rate of coarse sediment delivery from the landslide did not exceed the rivers ability to transport that material. The generally muted downstream response to sediment delivery from the State Route 530 Landslide, as well as the mechanics of that delivery and response, were generally consistent with observations made following the intentional removal of constructed dams. The rate and efficiency of erosion from the landslide decreased over the period of analysis, as the new channel approached a quasi-equilibrium form. In the absence of additional hillslope activity, rates of erosion from the landslide are likely to be small compared to those over the first 18 months after the landslide. The modest channel response to the highest rates of sediment delivery, and rapid recovery thereafter, indicate that the river should be able to convey the continued supply of landslide-derived sediment effectively with little effect on the downstream morphology and flood risks.

Washington↗

USGS2016

Transport and deposition of asbestos-rich sediment in the Sumas River, Whatcom County, Washington

Heavy sediment loads in the Sumas River of Whatcom County, Washington, increase seasonal turbidity and cause locally acute sedimentation. Most sediment in the Sumas River is derived from a deep-seated landslide of serpentinite that is located on Sumas Mountain and drained by Swift Creek, a tributary to the Sumas River. This mafic sediment contains high amounts of naturally occurring asbestiform chrysotile. A known human-health hazard, asbestiform chrysotile comprises 0.25–37 percent, by mass, of the total suspended sediment sampled from the Sumas River as part of this study, which included part of water year 2011 and all of water years 2012 and 2013. The suspended-sediment load in the Sumas River at South Pass Road, 0.6 kilometers (km) downstream of the confluence with Swift Creek, was 22,000 tonnes (t) in water year 2012 and 49,000 t in water year 2013. The suspended‑sediment load at Telegraph Road, 18.8 km downstream of the Swift Creek confluence, was 22,000 t in water year 2012 and 27,000 t in water year 2013. Although hydrologic conditions during the study were wetter than normal overall, the 2-year flood peak was only modestly exceeded in water years 2011 and 2013; runoff‑driven geomorphic disturbance to the watershed, which might have involved mass wasting from the landslide, seemed unexceptional. In water year 2012, flood peaks were modest, and the annual streamflow was normal. The fact that suspended-sediment loads in water year 2012 were equivalent at sites 0.6 and 18.8 km downstream of the sediment source indicates that the conservation of suspended‑sediment load can occur under normal hydrologic conditions. The substantial decrease in suspended-sediment load in the downstream direction in water year 2013 was attributed to either sedimentation in the intervening river reach, transfer to bedload as an alternate mode of sediment transport, or both. The sediment in the Sumas River is distinct from sediment in most other river systems because of the large percentage of asbestiform chrysotile in suspension. The suspended sediment carried by the Sumas River consists of three major components: (1) a relatively dense, largely non-flocculated material that settles rapidly out of suspension; (2) a lighter component containing relatively high proportions of flocculated material, much of it composed of asbestiform chrysotile; and (3) individual chrysotile fibers that are too small to flocculate or settle out, and remain in suspension as wash load (these fibers are on the order of microns in length and tenths of microns in diameter). Whereas the bulk density of the first (heaviest) component of suspended sediment was between 1.5 and 1.6 grams per cubic centimeter (g/cm3), the bulk density of the flocculated material was an order of magnitude lower (0.16 g/cm3), even after 24 hours of settling. Soon after immersion in water, the fresh chrysotile fibers derived from the Swift Creek landslide seem to flocculate readily into large bundles, or floccules, that exhibit settling velocities characteristic of coarse silts and fine sands (30 and 250 micrometers). In quiescent water within this river system, the floccules settle out quickly, but still leave between 2.4 and 19.5 million chrysotile fibers per liter in the clear overlying water. Consistent with the results from previous laboratory research, the amounts of asbestiform chrysotile in the water column in Swift Creek, as well as in the Sumas River close to and downstream of its confluence with Swift Creek, were determined to be directly correlated with pH. This observation offers a possible alternative to either turbidity or suspended‑sediment concentration as a surrogate for the concentration of fresh asbestiform chrysotile in suspension. Continued movement and associated erosion of the landslide through mass wasting and runoff will maintain large sediment loads in Swift Creek and in the Sumas River for the foreseeable future. Given the present channel morphology of the river system, aggradation (that is, sediment accumulation) in Swift Creek and the Sumas River are also likely to continue.

Washington↗