Search USGS⌕ Search

Geology topics

Jim E. O'Connor

Publications and source records attributed to Jim E. O'Connor.

54 records · Page 3Linked to original sources

Preliminary assessment of channel stability and bed-material transport in the Tillamook Bay tributaries and Nehalem River basin, northwestern Oregon

This report summarizes a preliminary study of bed-material transport, vertical and lateral channel changes, and existing datasets for the Tillamook (drainage area 156 square kilometers [km 2 ]), Trask (451 km 2 ), Wilson (500 km 2 ), Kilchis (169 km 2 ), Miami (94 km 2 ), and Nehalem (2,207 km 2 ) Rivers along the northwestern Oregon coast. This study, conducted in coopera-tion with the U.S. Army Corps of Engineers and Oregon Department of State Lands to inform permitting decisions regarding instream gravel mining, revealed that: Study areas along the six rivers can be divided into reaches based on tidal influence and topography. The fluvial (nontidal or dominated by riverine processes) reaches vary in length (2.4-9.3 kilometer [km]), gradient (0.0011-0.0075 meter of elevation change per meter of channel length [m/m]), and bed-material composition (a mixture of alluvium and intermittent bedrock outcrops to predominately alluvium). In fluvial reaches, unit bar area (square meter of bar area per meter of channel length [m 2 /m]) as mapped from 2009 photographs ranged from 7.1 m 2 /m on the Tillamook River to 27.9 m 2 /m on the Miami River. In tidal reaches, all six rivers flow over alluvial deposits, but have varying gradients (0.0001-0.0013 m/m) and lengths affected by tide (1.3-24.6 km). The Miami River has the steepest and shortest tidal reach and the Nehalem River has the flattest and longest tidal reach. Bars in the tidal reaches are generally composed of sand and mud. Unit bar area was greatest in the Tidal Nehalem Reach, where extensive mud flats flank the lower channel. Background factors such as valley and channel confinement, basin geology, channel slope, and tidal extent control the spatial variation in the accumulation and texture of bed material. Presently, the Upper Fluvial Wilson and Miami Reaches and Fluvial Nehalem Reach have the greatest abundance of gravel bars, likely owing to local bed-material sources in combination with decreasing channel gradient and valley confinement. Natural and human-caused disturbances such as mass movements, logging, fire, channel modifications for navigation and flood control, and gravel mining also have varying effects on channel condition, bed-material transport, and distribution and area of bars throughout the study areas and over time. Existing datasets include at least 16 and 18 sets of aerial and orthophotographs that were taken of the study areas in the Tillamook Bay tributary basins and Nehalem River basin, respectively, from 1939 to 2011. These photographs are available for future assessments of long-term changes in channel condition, bar area, and vegetation establishment patterns. High resolution Light Detection And Ranging (LiDAR) surveys acquired in 2007-2009 could support future quantitative analyses of channel morphology and bed-material transport in all study areas. A review of deposited and mined gravel volumes reported for instream gravel mining sites shows that bed-material deposition tends to rebuild mined bar surfaces in most years. Mean annual deposition volumes on individual bars exceeded 3,000 cubic meters (m 3 ) on Donaldson Bar on the Wilson River, Dill Bar on the Kilchis River, and Plant and Winslow Bars on the Nehalem River. Cumulative reported volumes of bed-material deposition were greatest at Donaldson and Dill Bars, totaling over 25,000 m 3 per site from 2004 to 2011. Within this period, reported cumulative mined volumes were greatest for the Donaldson, Plant, and Winslow Bars, ranging from 24,470 to 33,940 m 3 . Analysis of historical stage-streamflow data collected by the U.S. Geological Survey on the Wilson River near Tillamook (14301500) and Nehalem River near Foss (14301000) shows that these rivers have episodically aggraded and incised, mostly following high flow events, but they do not exhibit systematic, long-term trends in bed elevation. Multiple cross sections show that channels near bridge crossings in all six study areas are dynamic with many subject to incision and aggradation as well as lateral shifts in thalweg position and bank deposition and erosion. In fluvial reaches, unit bar area declined a net 5.3-83.6 percent from 1939 to 2009. The documented reduction in bar area may be attributable to several factors, including vegetation establishment and stabilization of formerly active bar surfaces, lateral channel changes and resulting alterations in sediment deposition and erosion patterns, and streamflow and/or tide differences between photographs. Other factors that may be associated with the observed reduction in bar area but not assessed in this reconnaissance level study include changes in the sediment and hydrology regimes of these rivers over the analysis period. In tidal reaches, unit bar area increased on the Tillamook and Nehalem Rivers (98.0 and 14.7 percent, respectively), but declined a net 24.2 to 83.1 percent in the other four tidal reaches. Net increases in bar area in the Tidal Tillamook and Nehalem Reaches were possibly attributable to tidal differences between the photographs as well as sediment deposition behind log booms and pile structures on the Tillamook River between 1939 and 1967. The armoring ratio (ratio of the median grain sizes of a bar's surface and subsurface layers) was 1.6 at Lower Waldron Bar on the Miami River, tentatively indicating a relative balance between transport capacity and sediment supply at this location. Armoring ratios, however, ranged from 2.4 to 5.5 at sites on the Trask, Wilson, Kilchis, and Nehalem Rivers; these coarse armor layers probably reflect limited bed-material supply at these sites. On the basis of mapping results, measured armoring ratios, and channel cross section surveys, preliminary conclusions are that the fluvial reaches on the Tillamook, Trask, Kilchis, and Nehalem Rivers are currently sediment supply-limited in terms of bed material - that is, the transport capacity of the channel generally exceeds the supply of bed material. The relation between transport capacity and sediment is more ambiguous for the fluvial reaches on the Wilson and Miami Rivers, but transport-limited conditions are likely for at least parts of these reaches. Some of these reaches have possibly evolved from sediment supply-limited to transport-limited over the last several decades in response to changing basin and climate conditions. Because of exceedingly low gradients, all the tidal reaches are transport-limited. Bed material in these reaches, however, is primarily sand and finer grain-size material and probably transported as suspended load from upstream reaches. These reaches will be most susceptible to watershed conditions affecting the supply and transport of fine sediment. Compared to basins on the southwestern Oregon coast, such as the Chetco and Rogue River basins, these six basins likely transport overall less gravel bed material. Although tentative in the absence of actual transport measurements, this conclusion is supported by the much lower area and frequency of bars and longer tidal reaches along all the northcoast rivers examined in this study. Previous studies suggest that the expansive and largely unvegetated bars visible in the 1939 photographs are primarily associated with voluminous sedimentation starting soon after the first Tillamook Burn fire in 1933. However, USGS studies of temporal bar trends in other Oregon coastal rivers unaffected by the Tillamook Burn show similar declines in bar area over approximately the same analysis period. In the Umpqua and Chetco River basins, historical declines in bar area are associated with long-term decreases in flood magnitude. Other factors may include changes in the type and volume of large wood and riparian vegetation. Further characterization of hydrology patterns in these basins and possible linkages with climate factors related to flood peaks, such as the Pacific Decadal Oscillation, could support inferences of expected future changes in vegetation establishment and channel planform and profile. More detailed investigations of bed-material transport rates and channel morphology would support assessments of lateral and vertical channel condition and longitudinal trends in bed material. Such assessments would be most practical for the fluvial study areas on the Wilson, Kilchis, Miami, and Nehalem Rivers and relevant to several ongoing management and ecological issues pertaining to sand and gravel transport. Tidal reaches may also be logical subjects for indepth analysis where studies would be more relevant to the deposition and transport of fine sediment (and associated channel and riparian conditions and processes) rather than coarse bed material.

Oregon↗

Geomorphic response of the Sandy River, Oregon, to removal of Marmot Dam

The October 2007 breaching of a temporary cofferdam constructed during removal of the 15-meter (m)-tall Marmot Dam on the Sandy River, Oregon, triggered a rapid sequence of fluvial responses as ~730,000 cubic meters (m 3 ) of sand and gravel filling the former reservoir became available to a high-gradient river. Using direct measurements of sediment transport, photogrammetry, airborne light detection and ranging (lidar) surveys, and, between transport events, repeat ground surveys of the reservoir reach and channel downstream, we monitored the erosion, transport, and deposition of this sediment in the hours, days, and months following breaching of the cofferdam. Rapid erosion of reservoir sediment led to exceptional suspended-sediment and bedload-sediment transport rates near the dam site, as well as to elevated transport rates at downstream measurement sites in the weeks and months after breaching. Measurements of sediment transport 0.4 kilometers (km) downstream of the dam site during and following breaching show a spike in the transport of fine suspended sediment within minutes after breaching, followed by high rates of suspended-load and bedload transport of sand. Significant transport of gravel bedload past the measurement site did not begin until 18 to 20 hours after breaching. For at least 7 months after breaching, bedload transport rates just below the dam site during high flows remained as much as 10 times above rates measured upstream of the dam site and farther downstream. The elevated sediment load was derived from eroded reservoir sediment, which began eroding when a meters-tall knickpoint migrated about 200 m upstream in the first hour after breaching. Rapid knickpoint migration triggered vertical incision and bank collapse in unconsolidated sand and gravel, leading to rapid channel widening. Over the following days and months, the knickpoint migrated upstream more slowly, simultaneously decreasing in height and becoming less distinct. Within 7 months, the knickpoint had migrated 2 km upstream from the dam site and became a riffle-like feature approximately 1 m high and a few tens of meters long. Knickpoint migration, vertical incision, and lateral erosion evacuated about 15 percent of the initial reservoir volume (125,000 m 3 ) within 60 hours following breaching, and by the end of the high flows in May 2008, about 50 percent of the volume had been evacuated. Large stormflows in November 2008 and January 2009 eroded another 6 percent of the original volume of impounded sediment. Little additional sediment eroded during the remainder of the second year following breaching. The rapid erosion of sediment by the modest flow that accompanied dam breaching was driven mainly by the steep hydraulic gradient associated with the abrupt change of base level and knickpoint formation and was aided by the unconsolidated and cohesionless character of the reservoir sediment. In the ensuing months, transport competence diminished as channel geometry evolved and the river gradient through the reservoir reach diminished. Changes in profile gradient in conjunction with channel coarsening and widening led to a rapid slowing of the rate of reservoir erosion. Sediment transport and deposition were strongly controlled by channel-gradient discontinuities and valley morphology downstream of the dam site. Those influences led to a strong divergence of sand and gravel transport and to deposition of a sediment wedge, as much as 4 m thick, that tapered to the preremoval channel bed 1.3 km downstream of the dam site. After 2 years, that deposit contained about 25 percent of the total volume of sediment eroded from the reservoir. The balance was distributed among pools within the Sandy River gorge, a narrow bedrock canyon extending 2 to 9 km downstream of the dam site, and along the channel farther downstream. A two-fraction sediment budget for the first year following breaching indicates that most of the gravel eroded from the reservoir reach was deposited within the sediment wedge and within the gorge, whereas eroded sand largely passed through the gorge and was broadly dispersed farther downstream. The sequence of transporting flows affected the specific trajectory of reservoir erosion and downstream sediment transport during the 2 years following breaching. However, because the overall erosion was largely a consequence of knickpoint retreat and channel widening, which in the 2 years after removal had affected most of the reservoir reach, it is unlikely that the specific sequence of flows significantly affected the overall outcome. Because the knickpoint had largely passed through the reservoir within 2 years, and the remaining reservoir sediment is mostly isolated high above armored or bedrock banks, it is unlikely that substantial additional sediment from the reservoir site will enter the system unless very large flows occur. Continued channel evolution downstream of the dam site is probable as deposits formed in the first 2 years are episodically mobilized. Below the Sandy River gorge, detection of effects related to release of reservoir sediment is challenging, especially in areas of sand deposition, because of the high background supply of sand in the river and substantial channel dynamism.

Oregon↗

Estimation of bed-material transport in the lower Chetco River, Oregon, water years 2009-2010

This assessment of bed-material transport uses methods developed in a previous study (Wallick and others, 2010) to estimate bed-material flux at the USGS Chetco River streamflow gaging station located at flood-plain kilometer 15 (14400000). On the basis of regressions between daily mean flow and transport capacity, daily bed-material flux was calculated for the period October 1, 2008 to March 30, 2011. The daily flux estimates were then aggregated by water year (WY) for WY 2009 and WY 2010 and the period April 1-March 31 during 2008-09, 2009-10 and 2010-11. The main findings were: *Estimated bed-material flux for WY 2009 (October 1, 2008 to September 30, 2009) was 87,300 metric tons as calculated by the Parker (1990a, b) equation (hereinafter \'the Parker equation\') and 116,900 metric tons as calculated by the Wilcock and Crowe (2003) equation (hereinafter \'the Wilcock-Crowe equation\'). *Estimated bed-material flux for water year 2010 (October 1, 2008 to September 30, 2009) was 56,800 metric tons as calculated by the Parker equation and 96,700 metric tons as calculated by the Wilcock-Crowe equation. *Estimated bed-material flux for April 1, 2008, to March 31, 2009, was 84,700 metric tons as calculated by the Parker equation and 111,700 metric tons as calculated by the Wilcock-Crowe equation. Flux values from April 1 to September 30, 2008, are from Wallick and others (2010). *Estimated bed-material flux for April 1, 2009, to March 31, 2010, was 45,500 metric tons as calculated by the Parker equation and 79,100 metric tons as calculated by the Wilcock-Crowe equation. *Estimated bed-material flux for April 1, 2010, to March 31, 2011, was 67,100 metric tons as calculated by the Parker equation and 134,300 metric tons as calculated by the Wilcock-Crowe equation. These calculations used provisional flow data for December 29, 2010, to March 31, 2011, and may be subject to revision. *Water years 2009 and 2010 both had less bed-material transport than the average annual transport values of 105,300 and 152,300 metric tons for the period 1970-2010 as calculated by the Parker and Wilcock-Crowe equations, respectively.

Oregon↗

Preliminary assessment of channel stability and bed-material transport along Hunter Creek, southwestern Oregon

This preliminary assessment of (1) bed-material transport in the Hunter Creek basin, (2) historical changes in channel condition, and (3) supplementary data needed to inform permitting decisions regarding instream gravel extraction revealed the following: Along the lower 12.4 km (kilometers) of Hunter Creek from its confluence with the Little South Fork Hunter Creek to its mouth, the river has confined and unconfined segments and is predominately alluvial in its lowermost 11 km. This 12.4-km stretch of river can be divided into two geomorphically distinct study reaches based primarily on valley physiography. In the Upper Study Reach (river kilometer [RKM] 12.4-6), the active channel comprises a mixed bed of bedrock, boulders, and smaller grains. The stream is confined in the upper 1.4 km of the reach by a bedrock canyon and in the lower 2.4 km by its valley. In the Lower Study Reach (RKM 6-0), where the area of gravel bars historically was largest, the stream flows over bed material that is predominately alluvial sediments. The channel alternates between confined and unconfined segments. The primary human activities that likely have affected bed-material transport and the extent and area of gravel bars are (1) historical and ongoing aggregate extraction from gravel bars in the study area and (2) timber harvest and associated road construction throughout the basin. These anthropogenic activities likely have varying effects on sediment transport and deposition throughout the study area and over time. Although assessing the relative effects of these anthropogenic activities on sediment dynamics would be challenging, the Hunter Creek basin may serve as a case study for such an assessment because it is mostly free of other alterations to hydrologic and geomorphic processes such as flow regulation, dredging, and other navigation improvements that are common in many Oregon coastal basins. Several datasets are available that may support a more detailed physical assessment of Hunter Creek. The entire study area has been captured in aerial photographs at least once per decade since the 1940s. This temporally rich photograph dataset would support quantitative analyses of changes in channel planform as well as vegetation cover. Light Detection And Ranging (LiDAR) data collected in 2008 would facilitate hydraulic and sediment-transport modeling and characterization of bar elevations throughout most of the study area. Few studies describing channel morphology and sediment transport exist for the Hunter Creek basin. The most detailed study reported channel incision and bank instability as well as the loss of point bars and pools in the lower 3.9 km of Hunter Creek from slightly downstream of its confluence with Yorke Creek to its mouth (EA Engineering, Sci-ence, and Technology, 1998). Repeat channel cross-sections collected from 1994 to 2010 at four bridges indicate that Hunter Creek is dynamic and subject to channel shifting, aggradation, and incision. Despite this dynamism, the channel at three bridge crossings showed little net change in thalweg elevation during this period. However, the channel thalweg aggraded 0.55 m from 2004 to 2008 near the bridge at RKM 3.5. Systematic delineation of gravel bars from aerial photographs collected in 1940, 1965, 2005, and 2009 indicates a 52-percent reduction in the area of bed-material sediment throughout the study area from 1940 to 2009. Net bar loss was greatest in the Lower Study Reach from RKM 1-4 and mainly is associ-ated with the encroachment of vegetation onto upper-bar surfaces lacking apparent vegetation in 1940. Bar-surface material was approximately equal in size to bar-subsurface material at Conn Creek Bar, whereas it was distinctly coarser than the subsurface material at Menasha Bar. Armoring ratios, which indicate the coarseness of the bar surface relative to the bar subsurface, were calculated as 0.97 for Conn Creek Bar and 1.5 for Menasha Bar. These ratios tentatively show that

Oregon↗

Columbia River Estuary ecosystem classification—Concept and application

This document describes the concept, organization, and application of a hierarchical ecosystem classification that integrates saline and tidal freshwater reaches of estuaries in order to characterize the ecosystems of large flood plain rivers that are strongly influenced by riverine and estuarine hydrology. We illustrate the classification by applying it to the Columbia River estuary (Oregon-Washington, USA), a system that extends about 233 river kilometers (rkm) inland from the Pacific Ocean. More than three-quarters of this length is tidal freshwater. The Columbia River Estuary Ecosystem Classification ("Classification") is based on six hierarchical levels, progressing from the coarsest, regional scale to the finest, localized scale: (1) Ecosystem Province; (2) Ecoregion; (3) Hydrogeomorphic Reach; (4) Ecosystem Complex; (5) Geomorphic Catena; and (6) Primary Cover Class. We define and map Levels 1-3 for the entire Columbia River estuary with existing geospatial datasets, and provide examples of Levels 4-6 for one hydrogeomorphic reach. In particular, three levels of the Classification capture the scales and categories of ecosystem structure and processes that are most tractable to estuarine research, monitoring, and management. These three levels are the (1) eight hydrogeomorphic reaches that embody the formative geologic and tectonic processes that created the existing estuarine landscape and encompass the influence of the resulting physiography on interactions between fluvial and tidal hydrology and geomorphology across 230 kilometers (km) of estuary, (2) more than 15 ecosystem complexes composed of broad landforms created predominantly by geologic processes during the Holocene, and (3) more than 25 geomorphic catenae embedded within ecosystem complexes that represent distinct geomorphic landforms, structures, ecosystems, and habitats, and components of the estuarine landscape most likely to change over short time periods.

Columbia River Estuary↗

Channel change and bed-material transport in the Umpqua River basin, Oregon

The Umpqua River drains 12,103 square kilometers of western Oregon; with headwaters in the Cascade Range, the river flows through portions of the Klamath Mountains and Oregon Coast Range before entering the Pacific Ocean. Above the head of tide, the Umpqua River, along with its major tributaries, the North and South Umpqua Rivers, flows on a mixed bedrock and alluvium bed, alternating between bedrock rapids and intermittent, shallow gravel bars composed of gravel to cobble-sized clasts. These bars have been a source of commercial aggregate since the mid-twentieth century. Below the head of tide, the Umpqua River contains large bars composed of mud and sand. Motivated by ongoing permitting and aquatic habitat concerns related to in-stream gravel mining on the fluvial reaches, this study evaluated spatial and temporal trends in channel change and bed-material transport for 350 kilometers of river channel along the Umpqua, North Umpqua, and South Umpqua Rivers. The assessment produced (1) detailed mapping of the active channel, using aerial photographs and repeat surveys, and (2) a quantitative estimation of bed-material flux that drew upon detailed measurements of particle size and lithology, equations of transport capacity, and a sediment yield analysis. Bed-material transport capacity estimates at 45 sites throughout the South Umpqua and main stem Umpqua Rivers for the period 1951-2008 result in wide-ranging transport capacity estimates, reflecting the difficulty of applying equations of bed-material transport to a supply-limited river. Median transport capacity values calculated from surface-based equations of bedload transport for each of the study reaches provide indications of maximum possible transport rates and range from 8,000 to 27,000 metric tons per year (tons/yr) for the South Umpqua River and 20,000 to 82,000 metric tons/yr for the main stem Umpqua River upstream of the head of tide; the North Umpqua River probably contributes little bed material. A plausible range of average annual transport rates for the South and main stem Umpqua Rivers, based on bedload transport capacity estimates for bars with reasonable values for reference shear stress, is between 500 and 20,000 metric tons/yr. An empirical bed-material yield analysis predicts 20,000-50,000 metric tons/yr on the South Umpqua River and main stem Umpqua River through the Oregon Coast Range, decreasing to approximately 30,000 metric tons/yr at the head of tide. Surveys of individual mining sites in the South Umpqua River indicate minimum local bed-material flux rates that are typically less than 10,000 metric tons/yr but range up to 30,600 metric tons/yr in high-flow years. On the basis of all of these analyses, actual bedload flux in most years is probably less than 25,000 metric tons/yr in the South Umpqua and main stem Umpqua Rivers, with the North Umpqua River probably contributing negligible amounts. For comparison, the estimated annual volume of commercial gravel extraction from the South Umpqua River between 2001 and 2004 ranged from 610 to 36,570 metric tons, indicating that historical in-stream gravel extraction may have been a substantial fraction of the overall bedload flux.

Oregon↗

Flood-frequency analyses from paleoflood investigations for Spring, Rapid, Boxelder, and Elk Creeks, Black Hills, western South Dakota

Flood-frequency analyses for the Black Hills area are important because of severe flooding of June 9-10, 1972, that was caused by a large mesoscale convective system and caused at least 238 deaths. Many 1972 peak flows are high outliers (by factors of 10 or more) in observed records that date to the early 1900s. An efficient means of reducing uncertainties for flood recurrence is to augment gaged records by using paleohydrologic techniques to determine ages and magnitudes of prior large floods (paleofloods). This report summarizes results of paleoflood investigations for Spring Creek, Rapid Creek (two reaches), Boxelder Creek (two subreaches), and Elk Creek. Stratigraphic records and resulting long-term flood chronologies, locally extending more than 2,000 years, were combined with observed and adjusted peak-flow values (gaged records) and historical flood information to derive flood-frequency estimates for the six study reaches. Results indicate that (1) floods as large as and even substantially larger than 1972 have affected most of the study reaches, and (2) incorporation of the paleohydrologic information substantially reduced uncertainties in estimating flood recurrence. Canyons within outcrops of Paleozoic rocks along the eastern flanks of the Black Hills provided excellent environments for (1) deposition and preservation of stratigraphic sequences of late-Holocene flood deposits, primarily in protected slack-water settings flanking the streams; and (2) hydraulic analyses for determination of associated flow magnitudes. The bedrock canyons ensure long-term stability of channel and valley geometry, thereby increasing confidence in hydraulic computations of ancient floods from modern channel geometry. Stratigraphic records of flood sequences, in combination with deposit dating by radiocarbon, optically stimulated luminescence, and cesium-137, provided paleoflood chronologies for 29 individual study sites. Flow magnitudes were estimated from elevations of flood deposits in conjunction with hydraulic calculations based on modern channel and valley geometry. Reach-scale paleoflood chronologies were interpreted for each study reach, which generally entailed correlation of flood evidence among multiple sites, chiefly based on relative position within stratigraphic sequences, unique textural characteristics, or results of age dating and flow estimation. The FLDFRQ3 and PeakfqSA analytical models (assuming log-Pearson Type III frequency distributions) were used for flood-frequency analyses for as many as four scenarios: (1) analysis of gaged records only; (2) gaged records with historical information; (3) all available data including gaged records, historical flows, paleofloods, and perception thresholds; and (4) the same as the third scenario, but ?top fitting? the distribution using only the largest 50 percent of gaged peak flows. The PeakfqSA model is most consistent with procedures adopted by most Federal agencies for flood-frequency analysis and thus was (1) used for comparisons among results for study reaches, and (2) considered by the authors as most appropriate for general applications of estimating low-probability flood recurrence. The detailed paleoflood investigations indicated that in the last 2,000 years all study reaches have had multiple large floods substantially larger than in gaged records. For Spring Creek, stratigraphic records preserved a chronology of at least five paleofloods in approximately (~) 1,000 years approaching or exceeding the 1972 flow of 21,800 cubic feet per second (ft 3 /s). The largest was ~700 years ago with a flow range of 29,300-58,600 ft 3 /s, which reflects the uncertainty regarding flood-magnitude estimates that was incorporated in the flood-frequency analyses. In the lower reach of Rapid Creek (downstream from Pactola Dam), two paleofloods in ~1,000 years exceeded the 1972 flow of 31,200 ft 3 /s. Those occurred ~440 and 1,000 years ago, with flows of 128,000-256,000 and 64,000-128,000 ft 3 /s, respectively. Five smaller paleofloods of 9,500-19,000 ft 3 /s occurred between ~200 and 400 years ago. In the upper reach of Rapid Creek (above Pactola Reservoir), the largest recorded floods are substantially smaller than for lower Rapid Creek and all other study reaches. Paleofloods of ~12,900 and 12,000 ft 3 /s occurred ~1,000 and 1,500 years ago. One additional paleoflood (~800 years ago) was similar in magnitude to the largest gaged flow of 2,460 ft 3 /s Boxelder Creek was treated as having two subreaches because of two tributaries that affect peak flows. During the last ~1,000 years, paleofloods of ~39,000-78,000 ft 3 /s and 40,000-80,000 ft 3 /s in the upstream subreach have exceeded the 1972 peak flow of 30,800 ft 3 /s. One other paleoflood was similar to the second largest gaged flow (16,400 ft 3 /s in 1907). For the downstream subreach, paleofloods of 61,300-123,000 ft 3 /s and 52,500-105,000 ft 3 /s in the last ~1,000 years have substantially exceeded the 1972 flood (50,500 ft 3 /s). Four additional paleofloods had flows between 14,200 and 33,800 ft 3 /s. The 1972 flow on Elk Creek (10,400 ft 3 /s) has been substantially exceeded at least five times in the last 1,900 years. The largest paleoflood (41,500-124,000 ft 3 /s) was ~900 years ago. Three other paleofloods between 37,500 and 120,000 ft 3 /s occurred between 1,100 and 1,800 years ago. A fifth paleoflood of 25,500-76,500 ft 3 /s was ~750 years ago. Considering analyses for all available data (PeakfqSA model) for all six study reaches, the 95-percent confidence intervals about the low-probability quantile estimates (100-, 200-, and 500-year recurrence intervals) were reduced by at least 78 percent relative to those for the gaged records only. In some cases, 95-percent uncertainty intervals were reduced by 99 percent or more. For all study reaches except the two Boxelder Creek subreaches, quantile estimates for these long-term analyses were larger than for the short-term analyses. The 1972 flow for the Spring Creek study reach (21,800 ft 3 /s) corresponds with a recurrence interval of ~400 years. Recurrence intervals are ~500 years for the 1972 flood magnitudes along the lower Rapid Creek reach and the upstream subreach of Boxelder Creek. For the downstream subreach of Boxelder Creek, the large 1972 flood magnitude (50,500 ft 3 /s) exceeds the 500-year quantile estimate by about 35 percent. The recurrence interval of ~100 years for 1972 flooding along the Elk Creek study reach is small relative to other study reaches along the eastern margin of the Black Hills. All of the paleofloods plot within the bounds of a national envelope curve, indicating that the national curve represents exceedingly rare floods for the Black Hills area. Elk Creek, lower Rapid Creek, and the downstream subreach of Boxelder Creek all have paleofloods that plot above a regional envelope curve; in the case of Elk Creek, by a factor of nearly two. The Black Hills paleofloods represent some of the largest known floods, relative to drainage area, for the United States. Many of the other largest known United States floods are in areas with physiographic and climatologic conditions broadly similar to the Black Hills-semiarid and rugged landscapes that intercept and focus heavy precipitation from convective storm systems. The 1972 precipitation and runoff patterns, previous analyses of peak-flow records, and the paleoflood investigations of this study support a hypothesis of distinct differences in flood generation within the central Black Hills study area. The eastern Black Hills are susceptible to intense orographic lifting associated with convective storm systems and also have high relief, thin soils, and narrow and steep canyons-factors favoring generation of exceptionally heavy rain-producing thunderstorms and promoting runoff and rapid concentration of flow into stream channels. In contrast, storm potential is smaller in and near the Limestone Plateau area, and storm runoff is further reduced by substantial infiltration into the limestone, gentle topography, and extensive floodplain storage. Results of the paleoflood investigations are directly applicable only to the specific study reaches and in the case of Rapid Creek, only to pre-regulation conditions. Thus, approaches for broader applications were developed from inferences of overall flood-generation processes, and appropriate domains for application of results were described. Example applications were provided by estimating flood quantiles for selected streamgages, which also allowed direct comparison with results of at-site flood-frequency analyses from a previous study. Several broad issues and uncertainties were examined, including potential biases associated with stratigraphic records that inherently are not always complete, uncertainties regarding statistical approaches, and the unknown applicability of paleoflood records to future watershed conditions. The results of the paleoflood investigations, however, provide much better physically based information on low-probability floods than has been available previously, substantially improving estimates of the magnitude and frequency of large floods in these basins and reducing associated uncertainty.

South Dakota↗

Channel change and bed-material transport in the Umpqua River basin, Oregon

The Umpqua River drains 12,103 km 2 of western Oregon, heading in the Cascade Range and draining portions of the Klamath Mountains and Coast Range before entering the Pacific Ocean. Above the head of tide, the Umpqua River, along with its major tributaries, the North and South Umpqua Rivers, flows on a mixed bedrock and alluvium bed, alternating between bedrock rapids and intermittent, shallow gravel bars composed of gravel to cobble-sized clasts. These bars have been a source of commercial aggregate since the mid-twentieth century. Below the head of tide, the Umpqua River contains large bars composed of mud and sand. Motivated by ongoing permitting and aquatic habitat concerns related to instream gravel mining on the fluvial reaches, this study evaluated spatial and temporal trends in channel change and bed-material transport for 350 km of river channel along the Umpqua, North Umpqua, and South Umpqua Rivers. The assessment produced (1) detailed mapping of the active channel, using aerial photographs and repeat surveys and (2) a quantitative estimation of bed-material flux that drew upon detailed measurements of particle size and lithology, equations of transport capacity, and a sediment yield analysis. Bed-material transport capacity estimates at 45 sites throughout the South Umpqua and mainstem Umpqua Rivers for the period 1951–2008 result in wide-ranging transport capacity estimates, reflecting the difficulty of applying equations of bed-material transport to a supply-limited river. Median transport capacity values calculated from surface-based equations of bedload transport for each of the study reaches provide indications of maximum possible transport rates and range from 8,000 to 27,000 metric tons/yr for the South Umpqua River and 20,000 to 82,000 metric tons/yr for the mainstem Umpqua River upstream of the head of tide; the North Umpqua River probably contributes little bed material. A plausible range of average annual transport rates for the South and mainstem Umpqua Rivers, based on bedload transport capacity estimates for bars with reasonable values for reference shear stress, is between 500 and 20,000 metric tons/yr. An empirical bed-material yield analysis predicts 20,000–50,000 metric tons/yr on the South Umpqua River and mainstem Umpqua River through the Coast Range, decreasing to approximately 30,000 metric tons/yr at the head of tide. Surveys of individual mining sites in the South Umpqua River indicate minimum local bed-material flux rates that are typically less than 10,000 metric tons/yr but that range up to 30,600 metric tons/yr in high-flow years. On the basis of all of these analyses, actual bedload flux in most years is probably less than 25,000 metric tons/yr in the South Umpqua River and Umpqua Rivers, with the North Umpqua River probably contributing negligible amounts. For comparison, the estimated annual volume of commercial gravel extraction from the South Umpqua River between 2001 and 2004 ranged from 610 to 36,570 metric tons, indicating that historical instream gravel extraction may have been a substantial fraction of the overall bedload flux.

Oregon↗

Channel change and bed-material transport in the Lower Chetco River, Oregon

The lower Chetco River is a wandering gravel-bed river flanked by abundant and large gravel bars formed of coarse bed-material sediment. Since the early twentieth century, the large gravel bars have been a source of commercial aggregate for which ongoing permitting and aquatic habitat concerns have motivated this assessment of historical channel change and sediment transport rates. Analysis of historical channel change and bed-material transport rates for the lower 18 kilometers shows that the upper reaches of the study area are primarily transport zones, with bar positions fixed by valley geometry and active bars mainly providing transient storage of bed material. Downstream reaches, especially near the confluence of the North Fork Chetco River, are zones of active sedimentation and channel migration. Multiple analyses, supported by direct measurements of bedload during winter 2008–09, indicate that since 1970 the mean annual flux of bed material into the study reach has been about 40,000–100,000 cubic meters per year. Downstream tributary input of bed-material sediment, probably averaging 5–30 percent of the influx coming into the study reach from upstream, is approximately balanced by bed-material attrition by abrasion. Probably little bed material leaves the lower river under natural conditions, with most net influx historically accumulating in wider and more dynamic reaches, especially near the North Fork Chetco River confluence, 8 kilometers upstream from the Pacific Ocean. The year-to-year flux, however, varies tremendously. Some years may have less than 3,000 cubic meters of bed material entering the study area; by contrast, some high-flow years, such as 1982 and 1997, likely have more than 150,000 cubic meters entering the reach. For comparison, the estimated annual volume of gravel extracted from the lower Chetco River for commercial aggregate during 2000–2008 has ranged from 32,000 to 90,000 cubic meters and averaged about 59,000 cubic meters per year. Mined volumes probably exceeded 140,000 cubic meters per year for several years in the late 1970s. Repeat surveys and map analyses indicate a reduction in bar area and sinuosity between 1939 and 2008, chiefly in the period 1965–95. Repeat topographic and bathymetric surveys show channel incision for substantial portions of the study reach, with local areas of bed lowering by as much as 2 meters. A specific gage analysis at the upstream end of the study reach indicates that incision and narrowing followed aggradation culminating in the late 1970s. These observations are all consistent with a reduction of sediment supply relative to transport capacity since channel surveys in the late 1970s, probably owing to a combination of (1) bed sediment removal and (2) transient river adjustments to large sediment volumes brought by floods such as those in 1964 and, to a lesser extent, 1996.

Oregon↗

Preliminary assessment of vertical stability and gravel transport along the Umpqua River, southwestern Oregon

This report addresses physical channel issues related to instream gravel mining on the Umpqua River and its two primary tributaries, the North and South Umpqua Rivers. This analysis constitutes a “Phase I” investigation, as designated by an interagency team cochaired by the U.S. Army Corps of Engineers, Portland District, and the Oregon Department of State Lands to address instream gravel mining issues across Oregon. Phase I analyses rely primarily on existing datasets and cursory analysis to determine the vertical stability of a channel to ascertain whether a particular river channel is aggrading, degrading, or at equilibrium. Additionally, a Phase I analysis identifies other critical issues or questions pertinent to physical channel conditions that may be related to instream gravel mining activities. This analysis can support agency permitting decisions as well as possibly indicating the need for additional studies. This specific analysis focuses on the mainstem Umpqua River from the Pacific Ocean at River Mile (RM) 0 to the confluence of the North and South Umpqua Rivers (at RM 111.8), as well as the lower 29 mi of the North Umpqua River and the lower 80 mi of the South Umpqua River (fig. 1). It is within these reaches where mining of gravel bars for aggregate has been most prevalent.

Oregon↗

Channel change and bed-material transport in the Lower Chetco River, Oregon

The lower Chetco River is a wandering gravel-bed river flanked by abundant and large gravel bars formed of coarse bed-material sediment. The large gravel bars have been a source of commercial aggregate since the early twentieth century for which ongoing permitting and aquatic habitat concerns have motivated this assessment of historical channel change and sediment transport rates. Analysis of historical channel change and bed-material transport rates for the lower 18 kilometers show that the upper reaches of the study area are primarily transport zones, with bar positions fixed by valley geometry and active bars mainly providing transient storage of bed material. Downstream reaches, especially near the confluence of the North Fork Chetco River, have been zones of active sedimentation and channel migration. Multiple analyses, supported by direct measurements of bedload during winter 2008–09, indicate that since 1970 the mean annual flux of bed material into the study reach has been about 40,000–100,000 cubic meters per year. Downstream tributary input of bed-material sediment, probably averaging 5–30 percent of the influx coming into the study reach from upstream, is approximately balanced by bed-material attrition by abrasion. Probably very little bed material leaves the lower river under natural conditions, with most of the net influx historically accumulating in wider and more dynamic reaches, especially near the North Fork Chetco River confluence, 8 kilometers upstream from the Pacific Ocean. The year-to-year flux, however, varies tremendously. Some years probably have less than 3,000 cubic meters of bed-material entering the study area; by contrast, some high-flow years, such as 1982 and 1997, likely have more than 150,000 cubic meters entering the reach. For comparison, the estimated annual volume of gravel extracted from the lower Chetco River for commercial aggregate during 2000–2008 has ranged from 32,000 to 90,000 cubic meters and averaged about 59,000 cubic meters per year. Mined volumes probably exceeded 140,000 cubic meters per year for several years in the late 1970s. Repeat surveys and map analyses indicate a reduction in bar area and sinuosity between 1939 and 2008, chiefly in the period 1965–95. Repeat topographic and bathymetric surveys show channel incision for substantial portions of the study reach, with local areas of bed lowering by as much as 2 meters. A specific gage analysis at the upstream end of the study reach indicates that incision and narrowing followed aggradation culminating in the late 1970s. These observations are all consistent with a reduction of sediment supply relative to transport capacity since channel surveys in the late 1970s, probably owing to a combination of (1) bed-sediment removal and (2) transient river adjustments to large sediment volumes brought by floods such as those in 1964, and to a lesser extent, 1996.

Open-File Report↗

The Portland Basin: A (big) river runs through it

Metropolitan Portland, Oregon, USA, lies within a small Neogene to Holocene basin in the forearc of the Cascadia subduction system. Although the basin owes its existence and structural development to its convergent-margin tectonic setting, the stratigraphic architecture of basin-fill deposits chiefly reflects its physiographic position along the lower reaches of the continental-scale Columbia River system. As a result of this globally unique setting, the basin preserves a complex record of aggradation and incision in response to distant as well as local tectonic, volcanic, and climatic events. Voluminous flood basalts, continental and locally derived sediment and volcanic debris, and catastrophic flood deposits all accumulated in an area influenced by contemporaneous tectonic deformation and variations in regional and local base level.

Oregon↗

Geologic map of the Camas Quadrangle, Clark County, Washington, and Multnomah County, Oregon

The Camas 7.5' quadrangle is in southwestern Washington and northwestern Oregon approximately 20 km east of Portland. The map area, bisected by the Columbia River, lies on the eastern margin of the Portland Basin, which is part of the Puget-Willamette Lowland that separates the Cascade Range from the Oregon Coast Range. Since late Eocene time, the Cascade Range has been the locus of an episodically active volcanic arc associated with underthrusting of oceanic lithosphere beneath the North American continent along the Cascadia Subduction Zone. Bedrock consists largely of basalt and basaltic andesite flows that erupted during late Oligocene time from one or more vents located outside the map area. These rocks crop out only north of the Columbia River: at the base of Prune Hill in Camas, where they dip southward at about 5°; and east of Lacamas Creek, where they dip to the southeast at 15 to 30°. The volcanic bedrock is unconformably overlain by Neogene sediments that accumulated as the Portland Basin subsided. In the Camas quadrangle, most of these sediments consist of basaltic hyaloclastic debris generated in the volcanic arc to the east and carried into the Portland Basin by the ancestral Columbia River. The dominant structures in the map area are northwest-striking dextral strike-slip faults that offset the Paleogene basin floor as well as the lower part of the basin fill. The Oligocene rocks at Prune Hill and to the east were uplifted in late Pliocene to early Pleistocene time within a restraining bend along one of these dextral faults. In Pleistocene time, basaltic andesite flows issued from a volcano centered on the west side of Prune Hill; another flow entered the map area from the east. These flows are part of the Boring volcanic field, which comprises several dozen late Pliocene and younger monogenetic volcanoes scattered throughout the greater Portland region. In latest Pleistocene time, the Missoula floods of glacial-outburst origin inundated the Portland Basin. The floods deposited huge bars of poorly sorted gravel in the lee of Prune Hill and west of the Sandy River. Volcanic debris from Mount Hood form a prominent delta at the mouth of the Sandy River. This map is a contribution to a program designed to improve geologic knowledge of the Portland Basin region of the Pacific Northwest urban corridor, the densely populated Cascadia forearc region of western Washington and Oregon. More detailed information on the bedrock and surficial geology of the basin and its surrounding area is necessary to refine assessments of seismic risk, ground-failure hazards and resource availability in this rapidly growing region.

Oregon, Washington↗

Spatial distribution of the largest rainfall‐runoff floods from basins between 2.6 and 26,000 km2 in the United States and Puerto Rico

We assess the spatial distribution of the largest rainfall‐generated streamflows from a database of 35,663 flow records composed of the largest 10% of annual peak flows from each of 14,815 U.S. Geological Survey stream gaging stations in the United States and Puerto Rico. High unit discharges (peak discharge per unit contributing area) from basins with areas of 2.6 to 26,000 km 2 (1–10,000 mi 2 ) are widespread, but streams in Hawaii, Puerto Rico, and Texas together account for more than 50% of the highest unit discharges. The Appalachians and western flanks of Pacific coastal mountain systems are also regions of high unit discharges, as are several areas in the southern Midwest. By contrast, few exceptional discharges have been recorded in the interior West, northern Midwest, and Atlantic Coastal Plain. Most areas of high unit discharges result from the combination of (1) regional atmospheric conditions that produce large precipitation volumes and (2) steep topography, which enhances precipitation by convective and orographic processes and allows flow to be quickly concentrated into stream channels. Within the conterminous United States, the greatest concentration of exceptional unit discharges is at the Balcones Escarpment of central Texas, where maximum U.S. rainfall amounts apparently coincide with appropriate basin physiography to produce many of the largest measured U.S. floods. Flood‐related fatalities broadly correspond to the spatial distribution of high unit discharges, with Texas having nearly twice the average annual flood‐related fatalities of any other state.

Water Resources Research↗

The world's largest floods, past and present: Their causes and magnitudes

Floods are among the most powerful forces on earth. Human societies worldwide have lived and died with floods from the very beginning, spawning a prominent role for floods within legends, religions, and history. Inspired by such accounts, geologists, hydrologists, and historians have studied the role of floods on humanity and its supporting ecosystems, resulting in new appreciation for the many-faceted role of floods in shaping our world. Part of this appreciation stems from ongoing analysis of long-term streamflow measurements, such as those recorded by the U.S. Geological Survey's (USGS) streamflow gaging network. But the recognition of the important role of flooding in shaping our cultural and physical landscape also owes to increased understanding of the variety of mechanisms that cause floods and how the types and magnitudes of floods can vary with time and space. The USGS has contributed to this understanding through more than a century of diverse research activities on many aspects of floods, including their causes, effects, and hazards. This Circular summarizes a facet of this research by describing the causes and magnitudes of the world's largest floods, including those measured and described by modern methods in historic times, as well as floods of prehistoric times, for which the only records are those left by the floods themselves.

Circular↗

Large floods in the United States: Where they happen and why

The spatial distribution of large gaged floods throughout the United States shows that the locations of most of the largest flows are related to specific combinations of regional climatology, topography, and basin size. Key factors include the general northward trend of decreasing atmospheric moisture, proximity to oceanic moisture sources such as the Pacific Ocean and the Gulf of Mexico, and orientation of topographic features relative to directions of moisture flow, with the largest flows being at locations where topographically high areas are oriented perpendicular to directions of moisture flow. Additionally, the largest flows in large river basins are primarily caused by persistent climatologic conditions such as seasonal snowfall. In contrast, the largest flows in smaller basins are most commonly the result of intense precipitation due to convective storms.

Circular↗

Origin, extent, and thickness of quaternary geologic units in the Willamette Valley, Oregon

Stratigraphic and chronologic information collected for Quaternary deposits in the Willamette Valley, Oregon, provides a revised stratigraphic framework that serves as a basis for a 1:250,000-scale map, as well as for thickness estimates of widespread Quaternary geologic units. We have mapped 11 separate Quaternary units that are differentiated on the basis of stratigraphic, topographic, pedogenic, and hydrogeologic properties. In summation, these units reflect four distinct episodes in the Quaternary geologic development of the Willamette Valley: 1) Fluvial sands and gravels that underlie terraces flanking lowland margins and tributary valleys were probably deposited between 2.5 and 0.5 million years ago. They are the oldest widespread surficial Quaternary deposits in the valley. Their present positions and preservation are undoubtedly due to postdepositional tectonic deformation - either by direct tectonic uplift of valley margins, or by regional tectonic controls on local base level. 2) Tertiary and Quaternary excavation or tectonic lowering of the Willamette Valley accommodated as much as 500 m (meters) of lacustrine and fluvial fill. Beneath the lowland floor, much of the upper 10 to 50 m of fill is Quaternary sand and gravel deposited by braided channel systems in subhorizontal sheets 2 to 10 m thick. These deposits grade to gravel fans 40 to 100 m thick where major Cascade Range rivers enter the valley and are traced farther upstream as much thinner valley trains of coarse gravel. The sand and gravel deposits have ages that range from greater than 420,000 to about 12,000 years old. A widely distributed layer of sand and gravel deposited at about 12 ka (kiloannum, thousands of years before the present) is looser and probably more permeable than older sand and gravel. Stratigraphic exposures and drillers' logs indicate that this late Pleistocene unit is mostly between 5 and 20 m thick where it has not been subsequently eroded by the Willamette River and its major tributaries. 3) Between 15,000 and 12,700 years ago, dozens of floods from Glacial Lake Missoula flowed up the Willamette Valley from the Columbia River, depositing up to 35 m of gravel, sand, silt, and clay. 4) Subsequent to 12,000 years ago, Willamette River sediment and flow regimes changed significantly: the Pleistocene braided river systems that had formed vast plains of sand and gravel evolved to incised and meandering rivers that are constructing today's fine-grained floodplains and gravelly channel deposits. Sub-surface channel facies of this unit are loose and unconsolidated and are highly permeable zones of substantial groundwater flow that is likely to be well connected to surface flow in the Willamette River and major tributaries. Stratigraphic exposures and drillers' logs indicate that this unit is mostly between 5 and 15 m thick.

Oregon↗

Methods for predicting peak discharge of floods caused by failure of natural and constructed earthen dams

Floods from failures of natural and constructed dams constitute a widespread hazard to people and property. Expeditious means of assessing flood hazards are necessary, particularly in the case of natural dams, which may form suddenly and unexpectedly. We revise statistical relations (derived from data for past constructed and natural dam failures) between peak discharge ( Q p ) and water volume released ( V 0 ) or drop in lake level ( d ) but assert that such relations, even when cast into a dimensionless form, are of limited utility because they fail to portray the effect of breach-formation rate. We then analyze a simple, physically based model of dam-breach formation to show that the hydrograph at the breach depends primarily on a dimensionless parameter η= kV 0 / g l/2 d 7/2 , where k is the mean erosion rate of the breach and g is acceleration due to gravity. The functional relationship between Q p and η takes asymptotically distinct forms depending on whether η ≪ 1 (relatively slow breach formation or small lake volume) or η ≫ 1 (relatively fast breach formation or large lake volume). Theoretical predictions agree well with data from dam failures for which k , and thus η, can be estimated. The theory thus provides a rapid means of predicting the plausible range of values of peak discharge at the breach in an earthen dam as long as the impounded water volume and the water depth at the dam face can be estimated.

Water Resources Research↗