USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Research”

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

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

At least 397 records · Page 22Linked to original sources
USGS2005

Old groundwater influence on stream hydrochemistry and catchment response times in a small Sierra Nevada catchment: Sagehen Creek, California

The relationship between the chemical and isotopic composition of groundwater and residence times was used to understand the temporal variability in stream hydrochemistry in Sagehen basin, California. On the basis of the relationship between groundwater age and [Ca 2+ ], the mean residence time of groundwater feeding Sagehen Creek during base flow is approximately 28 years. [Cl − ]:[Ca 2+ ] ratios in Sagehen Creek can be used to distinguish between two important processes: changes in the apparent age of groundwater discharging into the creek and dilution with snowmelt. The mean residence time of groundwater discharging into the creek is approximately 15 years during snowmelt periods. The results from this study have implications for hydrograph separation studies as groundwater is not a single, well‐mixed chemical component but rather is a variable parameter that predictably depends on groundwater residence time. Most current models of catchment hydrochemistry do not account for chemical and isotopic variability found within the groundwater reservoir. In addition, this study provides valuable insight into the long‐term hydrochemical response of a catchment to perturbations as catchment‐flushing times are related to the mean residence time of water in a basin.

California↗
USGS2004

Assessing denitrification in groundwater using natural gradient tracer tests with 15N: In situ measurement of a sequential multistep reaction

Denitrification was measured within a nitrate‐contaminated aquifer on Cape Cod, Massachusetts, using natural gradient tracer tests with 15 N nitrate. The aquifer contained zones of relatively high concentrations of nitrite (up to 77 μM) and nitrous oxide (up to 143 μM) and has been the site of previous studies examining ground water denitrification using the acetylene block technique. Small‐scale (15–24 m travel distance) tracer tests were conducted by injecting 15 N nitrate and bromide as tracers into a depth interval that contained nitrate, nitrite, nitrous oxide, and excess nitrogen gas. The timing of the bromide breakthrough curves at down‐gradient wells matched peaks in 15 N abundance above background for nitrate, nitrite, nitrous oxide, and nitrogen gas after more than 40 days of travel. Results were simulated with a one‐dimensional transport model using linked reaction kinetics for the individual steps of the denitrification reaction pathway. It was necessary to include within the model spatial variations in background concentrations of all nitrogen oxide species. The model indicated that nitrite production (0.036–0.047 μmol N (L aquifer) −1 d −1 ) was faster than the subsequent denitrification steps (0.013–0.016 μmol N (L aquifer) −1 d −1 for nitrous oxide and 0.013–0.020 μmol N (L aquifer) −1 d −1 for nitrogen gas) and that the total rate of reaction was slower than indicated by both acetylene block tracer tests and laboratory incubations. The rate of nitrate removal by denitrification was much slower than the rate of transport, indicating that nitrate would migrate several kilometers down‐gradient before being completely consumed.

Water Resources Research↗
USGS2005

Nitrogen transport and transformations in a coastal plain watershed: Influence of geomorphology on flow paths and residence times

Nitrogen transport and groundwater‐surface water interactions were examined in a coastal plain watershed in the southeastern United States. Groundwater age dates, calculated using chlorofluorocarbon and tritium concentrations, along with concentrations of nitrogen species and other redox‐active constituents, were used to evaluate the fate and transport of nitrate. Nitrate is stable only in recently recharged (<10 years) water found in the upper few meters of saturated thickness in the upland portion of a surficial aquifer. Groundwater with a residence time between 10 and 30 years typically has low nitrate and elevated excess N 2 concentrations, indications that denitrification has reduced nitrate concentrations. Groundwater older than 30 years also has low nitrate concentrations but contains little or no excess N 2 , suggesting that this water did not contain elevated concentrations of nitrate along its flow path. Nitrate transport to streams varies between first‐ and third‐order streams. Hydrologic, lithologic, and chemical data suggest that the surficial aquifer is the dominant source of flow and nitrate to a first‐order stream. Iron‐reducing conditions occur in groundwater samples from the bed and banks of the first‐order stream, suggesting that direct groundwater discharge is denitrified prior to entering the stream. However, nitrogen from the surficial aquifer is transported directly to the stream via a tile drain that bypasses these reduced zones. In the alluvial valley of a third‐order stream the erosion of a confining layer creates a much thicker unconfined alluvial aquifer with larger zones of nitrate stability. Age dating and chemical information (SiO 2 , Na/K ratios) suggest that water in the alluvial aquifer is derived from short flow paths through the riparian zone and/or from adjacent streams during high‐discharge periods.

Water Resources Research↗
USGS2004

Hydraulic tests in highly permeable aquifers

A semianalytical solution is presented for a mathematical model describing the flow of groundwater in response to a slug or pumping test in a highly permeable, confined aquifer. This solution, which is appropriate for wells of any degree of penetration and incorporates inertial mechanisms at both the test and observation wells, can be used to gain new insights into hydraulic tests in highly permeable settings. The oscillatory character of slug‐ and pumping‐induced responses will vary considerably across a site, even in an essentially homogeneous formation, when wells of different radii, depths, and screen lengths are used. Thus variations in the oscillatory character of responses do not necessarily indicate variations in hydraulic conductivity (K). Existing models for slug tests in partially penetrating wells in high‐K aquifers neglect the storage properties of the media. That assumption, however, appears reasonable for a wide range of common conditions. Unlike in less permeable formations, drawdown at an observation well in a high‐K aquifer will be affected by head losses in the pumping well. Those losses, which affect the form of the pumping‐induced oscillations, can be difficult to characterize. Thus analyses of observation‐well drawdown should utilize data from the period after the oscillations have dissipated whenever possible. Although inertial mechanisms can have a large impact on early‐time drawdown, that impact decreases rapidly with duration of pumping and distance to the observation well. Conventional methods that do not consider inertial mechanisms should therefore be viable options for the analysis of drawdown data at moderate to large times.

Water Resources Research↗
USGS2005

Speciation and transport of newly deposited mercury in a boreal forest wetland: A stable mercury isotope approach

As part of the Mercury Experiment to Assess Atmospheric Loadings in Canada and the United States (METAALICUS) the fate and transport of contemporary mercury (Hg) deposition in a boreal wetland was investigated using an experimentally applied stable mercury isotope. We applied high purity (99.2% ± 0.1) 202 Hg(II) to a wetland plot to determine if (1) the 202 Hg was detectable above the pool of native Hg, (2) the 202 Hg migrated vertically and/or horizontally in peat and pore waters, and (3) the 202 Hg was converted to methylmercury (MeHg) in situ. The 202 Hg was easily detected by ICP/MS in both solid peat and pore waters. Over 3 months, the 202 Hg migrated vertically downward in excess of 15 cm below the water table and traveled several meters horizontally beyond the experimental plot to the lake margin along the dominant vector of groundwater flow. Importantly, at one location, 6% of aqueous 202 Hg was detected as Me 202 Hg after only 1 day. These results indicate that new inorganic Hg in atmospheric deposition can be readily methylated and transported lakeward by shallow groundwater flow, confirming the important role of wetlands as contributors of Hg to aquatic ecosystems.

Water Resources Research↗
USGS2005

Variations in climate and ephemeral channel recharge in southeastern Arizona, United States

Significant variations in interannual and decadal recharge rates are likely in alluvial basins of the semiarid southwestern United States on the basis of decadal variations in climate and precipitation and correlation of El Niño with high rates of winter precipitation and streamflow. A better understanding of the magnitude of recharge variations in semiarid and arid regions would reduce water budget uncertainty. Variability of ephemeral channel recharge with climate in southeastern Arizona was investigated through analysis of hydrologic monitoring near three ephemeral streams in southeastern Arizona during the middle to late 1990s and by relating the results to long‐term hydrologic and climatic trends. The analysis used precipitation, streamflow, water levels in wells, estimates of groundwater storage change from repeat gravity surveys, and two climatic indicators of El Niño–Southern Oscillation (ENSO), Southern Oscillation index, and Pacific Decadal Oscillation (PDO). Results indicate that variations in winter recharge are related to ENSO. El Niño conditions correspond with a greater probability of high rates of winter precipitation, streamflow, and recharge. La Niña conditions are almost exclusively associated with below‐average recharge. Rates of recharge along Rillito Creek near Tucson during 1977–1998, a period of frequent El Niño conditions and positive PDO values, were 3 times recharge rates during 1941–1957, a period dominated by La Niña conditions and low PDO values. Quantification of recharge variability with decadal climate cycles should improve estimates of rates of aquifer drainage and replenishment in the region. Similar methods are applicable to other regions where thick unsaturated zones can accept significant periodic recharge.

Arizona↗
USGS2005

Influence of sediment storage on downstream delivery of contaminated sediment

Sediment storage in alluvial valleys can strongly modulate the downstream migration of sediment and associated contaminants through landscapes. Traditional methods for routing contaminated sediment through valleys focus on in‐channel sediment transport but ignore the influence of sediment exchanges with temporary sediment storage reservoirs outside the channel, such as floodplains. In theory, probabilistic analysis of particle trajectories through valleys offers a useful strategy for quantifying the influence of sediment storage on the downstream movement of contaminated sediment. This paper describes a field application and test of this theory, using 137 Cs as a sediment tracer over 45 years (1952–1997), downstream of a historical effluent outfall at the Los Alamos National Laboratory (LANL), New Mexico. The theory is parameterized using a sediment budget based on field data and an estimate of the 137 Cs release history at the upstream boundary. The uncalibrated model reasonably replicates the approximate magnitude and spatial distribution of channel‐ and floodplain‐stored 137 Cs measured in an independent field study. Model runs quantify the role of sediment storage in the long‐term migration of a pulse of contaminated sediment, quantify the downstream impact of upstream mitigation, and mathematically decompose the future 137 Cs flux near the LANL property boundary to evaluate the relative contributions of various upstream contaminant sources. The fate of many sediment‐bound contaminants is determined by the relative timescales of contaminant degradation and particle residence time in different types of sedimentary environments. The theory provides a viable approach for quantifying the long‐term movement of contaminated sediment through valleys.

Water Resources Research↗
USGS2004

A method for evaluating the importance of system state observations to model predictions, with application to the Death Valley regional groundwater flow system

We develop a new observation‐prediction (OPR) statistic for evaluating the importance of system state observations to model predictions. The OPR statistic measures the change in prediction uncertainty produced when an observation is added to or removed from an existing monitoring network, and it can be used to guide refinement and enhancement of the network. Prediction uncertainty is approximated using a first‐order second‐moment method. We apply the OPR statistic to a model of the Death Valley regional groundwater flow system (DVRFS) to evaluate the importance of existing and potential hydraulic head observations to predicted advective transport paths in the saturated zone underlying Yucca Mountain and underground testing areas on the Nevada Test Site. Important existing observations tend to be far from the predicted paths, and many unimportant observations are in areas of high observation density. These results can be used to select locations at which increased observation accuracy would be beneficial and locations that could be removed from the network. Important potential observations are mostly in areas of high hydraulic gradient far from the paths. Results for both existing and potential observations are related to the flow system dynamics and coarse parameter zonation in the DVRFS model. If system properties in different locations are as similar as the zonation assumes, then the OPR results illustrate a data collection opportunity whereby observations in distant, high‐gradient areas can provide information about properties in flatter‐gradient areas near the paths. If this similarity is suspect, then the analysis produces a different type of data collection opportunity involving testing of model assumptions critical to the OPR results.

Water Resources Research↗
USGS2005

Estimating soil hydraulic parameters from transient flow experiments in a centrifuge using parameter optimization technique

A modified version of the Hydrus software package that can directly or inversely simulate water flow in a transient centrifugal field is presented. The inverse solver for parameter estimation of the soil hydraulic parameters is then applied to multirotation transient flow experiments in a centrifuge. Using time‐variable water contents measured at a sequence of several rotation speeds, soil hydraulic properties were successfully estimated by numerical inversion of transient experiments. The inverse method was then evaluated by comparing estimated soil hydraulic properties with those determined independently using an equilibrium analysis. The optimized soil hydraulic properties compared well with those determined using equilibrium analysis and steady state experiment. Multirotation experiments in a centrifuge not only offer significant time savings by accelerating time but also provide significantly more information for the parameter estimation procedure compared to multistep outflow experiments in a gravitational field.

Water Resources Research↗
USGS2005

Surface water acidification responses and critical loads of sulfur and nitrogen deposition in Loch Vale watershed, Colorado

We evaluated the sensitivity of The Loch, a subalpine lake in Rocky Mountain National Park in Colorado, to acidification in response to increased atmospheric loading of sulfur (S) and nitrogen (N) using the Model of Acidification of Groundwater in Catchments (MAGIC). Lake water acid‐base chemistry was moderately sensitive to changes in both S and N deposition. However, the loads of S deposition that would drive chronic lake water acid neutralizing capacity (ANC) to below 0 or 20 μeq L −1 were estimated to be 11 and 8 kg S ha −1 yr −1 , respectively, assuming constant future N deposition at current levels. Comparable loads for N deposition, assuming constant future S deposition, were estimated to be 21 and 12 kg N ha −1 yr −1 , respectively. Modeling results for Andrews Creek, an alpine tributary to The Loch, suggested critical loads for surface water acidification that averaged about one third lower. Surface water ANC = 50 μeq L −1 was projected to occur in 50 years in The Loch if S or N deposition increased by a moderate amount (<40%) but could not be achieved in Andrews Creek by reducing either S or N deposition to zero. On the basis of the results of synoptic surveys of lake water chemistry, about one fifth of the wilderness lakes in the Colorado Front Range are more acid‐sensitive than The Loch. This modeling exercise suggests the need for a regional analysis of critical loads for the larger population of acid‐sensitive aquatic resources in order to provide part of the scientific foundation for federally mandated land management decisions.

Water Resources Research↗
USGS2005

Use of dissolved and vapor‐phase gases to investigate methanogenic degradation of petroleum hydrocarbon contamination in the subsurface

At many sites contaminated with petroleum hydrocarbons, methanogenesis is a significant degradation pathway. Techniques to estimate CH 4 production, consumption, and transport processes are needed to understand the geochemical system, provide a complete carbon mass balance, and quantify the hydrocarbon degradation rate. Dissolved and vapor‐phase gas data collected at a petroleum hydrocarbon contaminated site near Bemidji, Minnesota, demonstrate that naturally occurring nonreactive or relatively inert gases such as Ar and N 2 can be effectively used to better understand and quantify physical and chemical processes related to methanogenic activity in the subsurface. In the vadose zone, regions of Ar and N 2 depletion and enrichment are indicative of methanogenic and methanotrophic zones, and concentration gradients between the regions suggest that reaction‐induced advection can be an important gas transport process. In the saturated zone, dissolved Ar and N 2 concentrations are used to quantify degassing driven by methanogenesis and also suggest that attenuation of methane along the flow path, into the downgradient aquifer, is largely controlled by physical processes. Slight but discernable preferential depletion of N 2 over Ar, in both the saturated and unsaturated zones near the free‐phase oil, suggests reactivity of N 2 and is consistent with other evidence indicating that nitrogen fixation by microbial activity is taking place at this site.

Water Resources Research↗
USGS2005

Saline tracer visualized with three-dimensional electrical resistivity tomography: Field-scale spatial moment analysis

Cross-well electrical resistivity tomography (ERT) was used to monitor the migration of a saline tracer in a two-well pumping-injection experiment conducted at the Massachusetts Military Reservation in Cape Cod, Massachusetts. After injecting 2200 mg/L of sodium chloride for 9 hours, ERT data sets were collected from four wells every 6 hours for 20 days. More than 180,000 resistance measurements were collected during the tracer test. Each ERT data set was inverted to produce a sequence of 3-D snapshot maps that track the plume. In addition to the ERT experiment a pumping test and an infiltration test were conducted to estimate horizontal and vertical hydraulic conductivity values. Using modified moment analysis of the electrical conductivity tomograms, the mass, center of mass, and spatial variance of the imaged tracer plume were estimated. Although the tomograms provide valuable insights into field-scale tracer migration behavior and aquifer heterogeneity, standard tomographic inversion and application of Archie's law to convert electrical conductivities to solute concentration results in underestimation of tracer mass. Such underestimation is attributed to (1) reduced measurement sensitivity to electrical conductivity values with distance from the electrodes and (2) spatial smoothing (regularization) from tomographic inversion. The center of mass estimated from the ERT inversions coincided with that given by migration of the tracer plume using 3-D advective-dispersion simulation. The 3-D plumes seen using ERT exhibit greater apparent dispersion than the simulated plumes and greater temporal spreading than observed in field data of concentration breakthrough at the pumping well.

Water Resources Research↗
USGS2005

Solute transport and storage mechanisms in wetlands of the Everglades, south Florida

Solute transport and storage processes in wetlands play an important role in biogeochemical cycling and in wetland water quality functions. In the wetlands of the Everglades, there are few data or guidelines to characterize transport through the heterogeneous flow environment. Our goal was to conduct a tracer study to help quantify solute exchange between the relatively fast flowing water in the open part of the water column and much more slowly moving water in thick floating vegetation and in the pore water of the underlying peat. We performed a tracer experiment that consisted of a constant‐rate injection of a sodium bromide (NaBr) solution for 22 hours into a 3 m wide, open‐ended flume channel in Everglades National Park. Arrival of the bromide tracer was monitored at an array of surface water and subsurface samplers for 48 hours at a distance of 6.8 m downstream of the injection. A one‐dimensional transport model was used in combination with an optimization code to identify the values of transport parameters that best explained the tracer observations. Parameters included dimensions and mass transfer coefficients describing exchange with both short (hours) and longer (tens of hours) storage zones as well as the average rates of advection and longitudinal dispersion in the open part of the water column (referred to as the “main flow zone”). Comparison with a more detailed set of tracer measurements tested how well the model's storage zones approximated the average characteristics of tracer movement into and out of the layer of thick floating vegetation and the pore water in the underlying peat. The rate at which the relatively fast moving water in the open water column was exchanged with slowly moving water in the layer of floating vegetation and in sediment pore water amounted to 50 and 3% h −1 , respectively. Storage processes decreased the depth‐averaged velocity of surface water by 50% relative to the water velocity in the open part of the water column. As a result, flow measurements made with other methods that only work in the open part of the water column (e.g., acoustic Doppler) would have overestimated the true depth‐averaged velocity by a factor of 2. We hypothesize that solute exchange and storage in zones of floating vegetation and peat pore water increase contact time of solutes with biogeochemically active surfaces in this heterogeneous wetland environment.

Water Resources Research↗
USGS2005

Tropical cyclones and the flood hydrology of Puerto Rico

Some of the largest unit discharge flood peaks in the stream gaging records of the U.S. Geological Survey (USGS) have occurred in Puerto Rico. Many of these flood peaks are associated with tropical cyclones. Hurricane Georges, which passed directly over the island on 21–22 September 1998, produced record flood peaks at numerous USGS stations in Puerto Rico. The hydrology and hydrometeorology of extreme flood response in Puerto Rico are examined through analyses of rainfall, based on Weather Surveillance Radar–1988 Doppler (WSR‐88D) radar reflectivity observations and USGS rain gage observations and discharge from USGS stream gaging stations. Peak rainfall accumulations of more than 700 mm occurred in the central mountain region of the island. The largest unit discharge flood peaks, however, were located in the eastern portion of the island in areas with smaller storm total rainfall accumulations but markedly larger rainfall rates at 5–60 min timescale. Orographic precipitation mechanisms played an important role in rainfall distribution over the island of Puerto Rico. Amplification of rainfall accumulations was associated with areas of upslope motion. Elevated low‐level cloud water content in regions of upslope motion played an important role in the maximum rainfall accumulations in the central mountain region of Puerto Rico. The largest unit discharge flood peaks, however, were produced by a decaying eye wall mesovortex, which resulted in a 30–45 min period of extreme rainfall rates over the eastern portion of the island. This storm element was responsible for the record flood peak of the Río Grande de Lóiza. The role of terrain in development and evolution of the eye wall mesovortex is unclear but is of fundamental importance for assessing extreme flood response from the storm. Hydrologic response is examined through analyses of rainfall and discharge from five pairs of drainage basins, extending from east to west over the island. These analyses point to the importance of short‐term rainfall rates for extreme flood response. The hydrologic response of Puerto Rico is compared with two other extreme flood environments, the central Appalachians and Edwards Plateau of Texas. These analyses suggest that the high rainfall environment of Puerto Rico is linked to the development of a hydraulically efficient drainage system.

Water Resources Research↗
USGS2005

How snowpack heterogeneity affects diurnal streamflow timing

Diurnal cycles of streamflow in snow‐fed rivers can be used to infer the average time a water parcel spends in transit from the top of the snowpack to a stream gauge in the river channel. This travel time, which is measured as the difference between the hour of peak snowmelt in the afternoon and the hour of maximum discharge each day, ranges from a few hours to almost a full day later. Travel times increase with longer percolation times through deeper snowpacks, and prior studies of small basins have related the timing of a stream's diurnal peak to the amount of snow stored in a basin. However, in many larger basins the time of peak flow is nearly constant during the first half of the melt season, with little or no variation between years. This apparent self‐organization at larger scales can be reproduced by employing heterogeneous observations of snow depths and melt rates in a model that couples porous medium flow through an evolving snowpack with free surface flow in a channel.

Water Resources Research↗
USGS2005

The role of topography on catchment‐scale water residence time

The age, or residence time, of water is a fundamental descriptor of catchment hydrology, revealing information about the storage, flow pathways, and source of water in a single integrated measure. While there has been tremendous recent interest in residence time estimation to characterize watersheds, there are relatively few studies that have quantified residence time at the watershed scale, and fewer still that have extended those results beyond single catchments to larger landscape scales. We examined topographic controls on residence time for seven catchments (0.085–62.4 km 2 ) that represent diverse geologic and geomorphic conditions in the western Cascade Mountains of Oregon. Our primary objective was to determine the dominant physical controls on catchment‐scale water residence time and specifically test the hypothesis that residence time is related to the size of the basin. Residence times were estimated by simple convolution models that described the transfer of precipitation isotopic composition to the stream network. We found that base flow mean residence times for exponential distributions ranged from 0.8 to 3.3 years. Mean residence time showed no correlation to basin area (r 2 < 0.01) but instead was correlated (r 2 = 0.91) to catchment terrain indices representing the flow path distance and flow path gradient to the stream network. These results illustrate that landscape organization (i.e., topography) rather than basin area controls catchment‐scale transport. Results from this study may provide a framework for describing scale‐invariant transport across climatic and geologic conditions, whereby the internal form and structure of the basin defines the first‐order control on base flow residence time.

Water Resources Research↗
USGS2005

Incorporating seepage losses into the unsteady streamflow equations for simulating intermittent flow along mountain front streams

Seepage losses along numerous mountain front streams that discharge intermittently onto alluvial fans and piedmont alluvial plains are an important source of groundwater in the Basin and Range Province of the Western United States. Determining the distribution of seepage loss along mountain front streams is important when assessing groundwater resources of the region. Seepage loss along a mountain front stream in northern Nevada was evaluated using a one‐dimensional unsteady streamflow model. Seepage loss was incorporated into the spatial derivatives of the streamflow equations. Because seepage loss from streams is dependent on stream depth, wetted perimeter, and streambed properties, a two‐dimensional variably saturated flow model was used to develop a series of relations between seepage loss and stream depth for each reach. This method works when streams are separated from groundwater by variably saturated sediment. Two periods of intermittent flow were simulated to evaluate the modeling approach. The model reproduced measured flow and seepage losses along the channel. Seepage loss in the spring of 2000 was limited to the upper reaches on the alluvial plain and totaled 196,000 m 3 , whereas 64% of the seepage loss in the spring of 2004 occurred at the base of the alluvial plain and totaled 273,000 m 3 . A greater seepage loss at the base of the piedmont alluvial plain is attributed to increased streambed hydraulic conductivity caused by less armoring of the channel. The modeling approach provides a method for quantifying and distributing seepage loss along mountain front streams that cross alluvial fans or piedmont alluvial plains.

Water Resources Research↗
USGS2005

Variation in the reference Shields stress for bed load transport in gravel‐bed streams and rivers

The present study examines variations in the reference shear stress for bed load transport (τ r ) using coupled measurements of flow and bed load transport in 45 gravel‐bed streams and rivers. The study streams encompass a wide range in bank‐full discharge (1–2600 m 3 /s), average channel gradient (0.0003–0.05), and median surface grain size (0.027–0.21 m). A bed load transport relation was formed for each site by plotting individual values of the dimensionless transport rate W * versus the reach‐average dimensionless shear stress τ*. The reference dimensionless shear stress τ* r was then estimated by selecting the value of τ* corresponding to a reference transport rate of W * = 0.002. The results indicate that the discharge corresponding to τ* r averages 67% of the bank‐full discharge, with the variation independent of reach‐scale morphologic and sediment properties. However, values of τ* r increase systematically with average channel gradient, ranging from 0.025–0.035 at sites with slopes of 0.001–0.006 to values greater than 0.10 at sites with slopes greater than 0.02. A corresponding relation for the bank‐full dimensionless shear stress τ* bf , formulated with data from 159 sites in North America and England, mirrors the relation between τ* r and channel gradient, suggesting that the bank‐full channel geometry of gravel‐ and cobble‐bedded streams is adjusted to a relatively constant excess shear stress, τ* bf − τ* r , across a wide range of slopes.

Water Resources Research↗