Geology Reports⌕ Search

Geology topics

Richard R. McDonald

Publications and source records attributed to Richard R. McDonald.

32 records · Page 2Linked to original sources
USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗
USGS2015

Coevolution of bed surface patchiness and channel morphology: 2. Numerical experiments

In gravel bed rivers, bed topography and the bed surface grain size distribution evolve simultaneously, but it is not clear how feedbacks between topography and grain sorting affect channel morphology. In this, the second of a pair of papers examining interactions between bed topography and bed surface sorting in gravel bed rivers, we use a two-dimensional morphodynamic model to perform numerical experiments designed to explore the coevolution of both free and forced bars and bed surface patches. Model runs were carried out on a computational grid simulating a 200 m long, 2.75 m wide, straight, rectangular channel, with an initially flat bed at a slope of 0.0137. Over five numerical experiments, we varied (a) whether an obstruction was present, (b) whether the sediment was a gravel mixture or a single size, and (c) whether the bed surface grain size feeds back on the hydraulic roughness field. Experiments with channel obstructions developed a train of alternate bars that became stationary and were connected to the obstruction. Freely migrating alternate bars formed in the experiments without channel obstructions. Simulations incorporating roughness feedbacks between the bed surface and flow field produced flatter, broader, and longer bars than simulations using constant roughness or uniform sediment. Our findings suggest that patches are not simply a by-product of bed topography, but they interact with the evolving bed and influence morphologic evolution.

Journal of Geophysical Research F: Earth Surface↗
USGS2013

Prediction of suspended-sediment concentrations at selected sites in the Fountain Creek watershed, Colorado, 2008-09

In 2008, the U.S. Geological Survey (USGS), in cooperation with Pikes Peak Area Council of Governments, Colorado Water Conservation Board, Colorado Springs City Engineering, and the Lower Arkansas Valley Water Conservancy District, began a small-scale pilot study to evaluate the effectiveness of the use of a computational model of streamflow and suspended-sediment transport for predicting suspended-sediment concentrations and loads in the Fountain Creek watershed in Colorado. Increased erosion and sedimentation damage have been identified by the Fountain Creek Watershed Plan as key problems within the watershed. A recommendation in the Fountain Creek Watershed plan for management of the basin is to establish measurable criteria to determine if progress in reducing erosion and sedimentation damage is being made. The major objective of this study was to test a computational method to predict local suspended-sediment loads at two sites with different geomorphic characteristics in order to evaluate the feasibility of using such an approach to predict local suspended-sediment loads throughout the entire watershed. Detailed topographic surveys, particle-size data, and suspended-sediment samples were collected at two gaged sites: Monument Creek above Woodmen Road at Colorado Springs, Colorado (USGS gage 07103970), and Sand Creek above mouth at Colorado Springs, Colorado (USGS gage 07105600). These data were used to construct three-dimensional computational models of relatively short channel reaches at each site. The streamflow component of these models predicted a spatially distributed field of water-surface elevation, water velocity, and bed shear stress for a range of stream discharges. Using the model predictions, along with measured particle sizes, the sediment-transport component of the model predicted the suspended-sediment concentration throughout the reach of interest. These computed concentrations were used with predicted flow patterns and channel morphology to determine fluxes of suspended sediment for the median particle size and for the measured range of particle sizes in the channel. Three different techniques were investigated for making the suspended-sediment predictions; these techniques have varying degrees of reliance on measured data and also have greatly differing degrees of complexity. Based on these data, the calibrated Rouse method provided the best balance between accuracy and both computational and data collection costs; the presence of substantial washload was the primary factor in eliminating the simpler and the more complex techniques. Based on this work, using the selected technique at additional sites in the watershed to determine relative loads and source areas appears plausible. However, to ensure that the methodology presented in this report yields reasonable results at other selected sites in the basin, it is necessary to collect additional verification data sets at those locations.

Colorado↗
USGS2012

Using computational modeling of river flow with remotely sensed data to infer channel bathymetry

As part of an ongoing investigation into the use of computational river flow and morphodynamic models for the purpose of correcting and extending remotely sensed river datasets, a simple method for inferring channel bathymetry is developed and discussed. The method is based on an inversion of the equations expressing conservation of mass and momentum to develop equations that can be solved for depth given known values of vertically-averaged velocity and water-surface elevation. The ultimate goal of this work is to combine imperfect remotely sensed data on river planform, water-surface elevation and water-surface velocity in order to estimate depth and other physical parameters of river channels. In this paper, the technique is examined using synthetic data sets that are developed directly from the application of forward two-and three-dimensional flow models. These data sets are constrained to satisfy conservation of mass and momentum, unlike typical remotely sensed field data sets. This provides a better understanding of the process and also allows assessment of how simple inaccuracies in remotely sensed estimates might propagate into depth estimates. The technique is applied to three simple cases: First, depth is extracted from a synthetic dataset of vertically averaged velocity and water-surface elevation; second, depth is extracted from the same data set but with a normally-distributed random error added to the water-surface elevation; third, depth is extracted from a synthetic data set for the same river reach using computed water-surface velocities (in place of depth-integrated values) and water-surface elevations. In each case, the extracted depths are compared to the actual measured depths used to construct the synthetic data sets (with two- and three-dimensional flow models). Errors in water-surface elevation and velocity that are very small degrade depth estimates and cannot be recovered. Errors in depth estimates associated with assuming water-surface velocities equal to depth-integrated velocities are substantial, but can be reduced with simple corrections.

Conference Paper↗
USGS2011

Mechanics of flow and sediment transport in delta distributary channels

Predicting the planform and dimensions of a channel downstream from a confluence of two smaller channels with known sediment and water supplies is a fundamental, well-studied problem in geomorphology and engineering. An analogous but less well understood problem is found well downstream of such confluences, where large river channels split into two or more distributary channels on a river delta. In this case, both the flow and sediment supplies in the downstream distributaries are set by the dynamics near the bifurcation of the upstream channel and by downstream boundary conditions. Over time, the pattern of erosion and deposition in the distributary channels gives rise to variations in the amount of water and sediment routed into them. In the simplest case, this results in channel switching on deltas, but in a more general sense these dynamics produce a rich suite of interesting morphologic change contributing both to the evolution of the channel distributary network and the overall evolution of the delta. As part of a study to develop a better understanding of these processes, we conducted a field study measuring the detailed morphology of the Hong-Luoc junction on the Red River downstream of Hanoi, Vietnam. This junction was selected for such a study because it has a 1,000-year history, modern observations suggest that it is currently switching (changing the proportion of sediment and streamflow provided to each of the distributary channels), and hydrologic configuration of the junction allows for the study of two bifurcations and one confluence in a single junction complex. In this paper, our morphologic observations are used in computational flow models to show how flow and sediment transport changes as a function of total discharge upstream of the junction. This is a key component of understanding how the junction attains stability over a range of flows or how imbalances in the distribution of flow and sediment transport lead to destabilization of the channel bifurcation.

Conference Paper↗
USGS2010

Sediment-transport during three controlled-flood experiments on the Colorado River downstream from Glen Canyon Dam, with implications for eddy-sandbar deposition in Grand Canyon National Park

Three large-scale field experiments were conducted on the Colorado River downstream from Glen Canyon Dam in 1996, 2004, and 2008 to evaluate whether artificial (that is, controlled) floods released from the dam could be used in conjunction with the sand supplied by downstream tributaries to rebuild and sustainably maintain eddy sandbars in the river in Grand Canyon National Park. Higher suspended-sand concentrations during a controlled flood will lead to greater eddy-sandbar deposition rates. During each controlled flood experiment, sediment-transport and bed-sediment data were collected to evaluate sediment-supply effects on sandbar deposition. Data collection substantially increased in spatial and temporal density with each subsequent experiment. The suspended- and bed-sediment data collected during all three controlled-flood experiments are presented and analyzed in this report. Analysis of these data indicate that in designing the hydrograph of a controlled flood that is optimized for sandbar deposition in a given reach of the Colorado River, both the magnitude and the grain size of the sand supply must be considered. Because of the opposing physical effects of bed-sand area and bed-sand grain size in regulating suspended-sand concentration, larger amounts of coarser sand on the bed can lead to lower suspended-sand concentrations, and thus lower rates of sandbar deposition, during a controlled flood than can lesser amounts of finer sand on the bed. Although suspended-sand concentrations were higher at all study sites during the 2008 controlled-flood experiment (CFE) than during either the 1996 or 2004 CFEs, these higher concentrations were likely associated with more sand on the bed of the Colorado River in only lower Glen Canyon. More sand was likely present on the bed of the river in Grand Canyon during the 1996 CFE than during either the 2004 or 2008 CFEs. The question still remains as to whether sandbars can be sustained in the Colorado River in Grand Canyon National Park through use of controlled floods in conjunction with typical amounts and grain sizes of sand supplied by the tributaries that enter the Colorado River downstream from Glen Canyon Dam.

Arizona↗
USGS2010

Topographic evolution of sandbars: Flume experiment and computational modeling

Measurements of sandbar formation and evolution were carried out in a laboratory flume and the topographic characteristics of these barforms were compared to predictions from a computational flow and sediment transport model with bed evolution. The flume experiment produced sandbars with approximate mode 2, whereas numerical simulations produced a bed morphology better approximated as alternate bars, mode 1. In addition, bar formation occurred more rapidly in the laboratory channel than for the model channel. This paper focuses on a steady-flow laboratory experiment without upstream sediment supply. Future experiments will examine the effects of unsteady flow and sediment supply and the use of numerical models to simulate the response of barform topography to these influences.

Conference Paper↗
USGS2010

Mechanics and modeling of flow, sediment transport and morphologic change in riverine lateral separation zones

Lateral separation zones or eddies in rivers are critically important features for sediment storage and for a variety of roles they play in riparian and aquatic ecology. As part of a larger effort to predict the morphology of lateral separation zones in the Colorado River in Grand Canyon for a selection of sediment supply and discharge scenarios, we evaluated the performance of two modeling techniques for predicting flow, sediment transport, and morphodynamics in eddies using field data. In order to understand the relative roles of various exchange mechanisms between the main channel and eddies, we applied two-dimensional unsteady and three-dimensional unsteady models in a reach containing a lateral separation zone. Both models were developed, calibrated, and evaluated using detailed field data comprising acoustic-Doppler velocity measurements, water-surface elevations, sediment concentration by size class, and bathymetry measured during a flood event in the Colorado River. Model results and measurements are used to develop a better understanding of the mechanics of water and sediment exchange between the eddy and the mainstem and other factors that control the morphology of the reach.

Arizona↗
USGS2010

Computational modeling of bedform evolution in rivers with implications for predictions of flood stage and bed evolution

Uncertainties in flood stage prediction and bed evolution in rivers are frequently associated with the evolution of bedforms over a hydrograph. For the case of flood prediction, the evolution of the bedforms may alter the effective bed roughness, so predictions of stage and velocity based on assuming bedforms retain the same size and shape over a hydrograph will be incorrect. These same effects will produce errors in the prediction of the sediment transport and bed evolution, but in this latter case the errors are typically larger, as even small errors in the prediction of bedform form drag can make very large errors in predicting the rates of sediment motion and the associated erosion and deposition. In situations where flows change slowly, it may be possible to use empirical results that relate bedform morphology to roughness and effective form drag to avoid these errors; but in many cases where the bedforms evolve rapidly and are in disequilibrium with the instantaneous flow, these empirical methods cannot be accurately applied. Over the past few years, computational models for bedform development, migration, and adjustment to varying flows have been developed and tested with a variety of laboratory and field data. These models, which are based on detailed multidimensional flow modeling incorporating large eddy simulation, appear to be capable of predicting bedform dimensions during steady flows as well as their time dependence during discharge variations. In the work presented here, models of this type are used to investigate the impacts of bedform on stage and bed evolution in rivers during flood hydrographs. The method is shown to reproduce hysteresis in rating curves as well as other more subtle effects in the shape of flood waves. Techniques for combining the bedform evolution models with larger-scale models for river reach flow, sediment transport, and bed evolution are described and used to show the importance of including dynamic bedform effects in river modeling. For example calculations for a flood on the Kootenai River, errors of almost 1m in predicted stage and errors of about a factor of two in the predicted maximum depths of erosion can be attributed to bedform evolution. Thus, treating bedforms explicitly in flood and bed evolution models can decrease uncertainty and increase the accuracy of predictions.

Conference Paper↗
USGS2009

Simulation of streamflow using a multidimensional flow model for white sturgeon habitat, Kootenai River near Bonners Ferry, Idaho – Supplement to Scientific Investigations Report 2005-5230

During 2005, the U.S. Geological Survey (USGS) developed, calibrated, and validated a multidimensional flow model for simulating streamflow in the white sturgeon spawning habitat of the Kootenai River in Idaho. The model was developed as a tool to aid understanding of the physical factors affecting quality and quantity of spawning and rearing habitat used by the endangered white sturgeon (Acipenser transmontanus) and for assessing the feasibility of various habitat-enhancement scenarios to re-establish recruitment of white sturgeon. At the request of the Kootenai Tribe of Idaho, the USGS extended the two-dimensional flow model developed in 2005 into a braided reach upstream of the current white sturgeon spawning reach. Many scientists consider the braided reach a suitable substrate with adequate streamflow velocities for re-establishing recruitment of white sturgeon. The 2005 model was extended upstream to help assess the feasibility of various strategies to encourage white sturgeon to spawn in the reach. At the request of the Idaho Department of Fish and Game, the USGS also extended the two-dimensional flow model several kilometers downstream of the white sturgeon spawning reach. This modified model can quantify the physical characteristics of a reach that white sturgeon pass through as they swim upstream from Kootenay Lake to the spawning reach. The USGS Multi-Dimensional Surface-Water Modeling System was used for the 2005 modeling effort and for this subsequent modeling effort. This report describes the model applications and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to the location of white sturgeon spawning events during 1994-2001. Model simulations also were used to report on the length and percentage of longitudinal profiles that met the minimum criteria during May and June 2006 and 2007 as stipulated in the U.S. Fish and Wildlife Biological Opinion.

Idaho↗
USGS2007

Surface-exposure ages of Front Range moraines that may have formed during the Younger Dryas, 8.2 cal ka, and Little Ice Age events

Surface-exposure ( 10 Be) ages have been obtained on boulders from three post-Pinedale end-moraine complexes in the Front Range, Colorado. Boulder rounding appears related to the cirque-to-moraine transport distance at each site with subrounded boulders being typical of the 2-km-long Chicago Lakes Glacier, subangular boulders being typical of the 1-km-long Butler Gulch Glacier, and angular boulders being typical of the few-hundred-m-long Isabelle Glacier. Surface-exposure ages of angular boulders from the Isabelle Glacier moraine, which formed during the Little Ice Age (LIA) according to previous lichenometric dating, indicate cosmogenic inheritance values ranging from 0 to ∼3.0 10 Be ka. 1 Subangular boulders from the Butler Gulch end moraine yielded surface-exposure ages ranging from 5 to 10.2 10 Be ka. We suggest that this moraine was deposited during the 8.2 cal ka event, which has been associated with outburst floods from Lake Agassiz and Lake Ojibway, and that the large age range associated with the Butler Gulch end moraine is caused by cosmogenic shielding of and(or) spalling from boulders that have ages in the younger part of the range and by cosmogenic inheritance in boulders that have ages in the older part of the range. The surface-exposure ages of eight of nine subrounded boulders from the Chicago Lakes area fall within the 13.0–11.7 10 Be ka age range, and appear to have been deposited during the Younger Dryas interval. The general lack of inheritance in the eight samples probably stems from the fact that only a few thousand years intervened between the retreat of the Pinedale glacier and the advance of the Chicago Lakes glacier; in addition, bedrock in the Chicago Lakes cirque area may have remained covered with snow and ice during that interval, thus partially shielding the bedrock from cosmogenic radiation.

Colorado↗
USGS2006

A broadscale fish-habitat model development process: Genesee Basin, New York

We describe a methodology for developing species-habitat models using available fish and stream habitat data from New York State, focusing on the Genesee basin. Electrofishing data from the New York Department of Environmental Conservation were standardized and used for model development and testing. Four types of predictive models (multiple linear regression, stepwise multiple linear regression, linear discriminant analysis, and neural network) were developed and compared for 11 fish species. Predictive models used as many as 25 habitat variables and explained 35-91% of observed species abundance variability. Omission rates were generally low, but commission rates varied widely. Neural network models performed best for all species, except for rainbow trout Oncorhynchus mykiss , gizzard shad Dorosoma cepedianum , and brown trout Salmo trutta . Linear discriminant functions generally performed poorly. The species-environment models we constructed performed well and have potential applications to management issues.

Book chapter↗
USGS2005

Simulation of flow and sediment mobility using a multidimensional flow model for the White Sturgeon critical-habitat reach, Kootenai River near Bonners Ferry, Idaho

In 1994, the Kootenai River white sturgeon (Acipenser transmontanus) was listed as an Endangered Species as a direct result of two related observations. First, biologists observed that the white sturgeon population in the Kootenai River was declining. Second, they observed a decline in recruitment of juvenile sturgeon beginning in the 1950s with an almost total absence of recruitment since 1974, following the closure of Libby Dam in 1972. This second observation was attributed to changes in spawning and (or) rearing habitat resulting from alterations in the physical habitat, including flow regime, sediment-transport regime, and bed morphology of the river. The Kootenai River White Sturgeon Recovery Team was established to find and implement ways to improve spawning and rearing habitat used by white sturgeon. They identified the need to develop and apply a multidimensional flow model to certain reaches of the river to quantify physical habitat in a spatially distributed manner. The U.S. Geological Survey has addressed these needs by developing, calibrating, and validating a multidimensional flow model used to simulate streamflow and sediment mobility in the white sturgeon critical-habitat reach of the Kootenai River. This report describes the model and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to biological or other habitat data. This study was conducted in cooperation with the Kootenai Tribe of Idaho along a 23-kilometer reach of the Kootenai River, including the white sturgeon spawning reach near Bonners Ferry, Idaho that is about 108 to 131 kilometers below Libby Dam. U.S. Geological Survey's MultiDimensional Surface-Water Modeling System was used to construct a flow model for the critical-habitat reach of the Kootenai River white sturgeon, between river kilometers 228.4 and 245.9. Given streamflow, bed roughness, and downstream water-surface elevation, the model computes the velocity field, water-surface elevations, and boundary shear stress throughout the modeled reach. The 17.5 kilometer model reach was subdivided into two segments on the basis of predominant grain size: a straight reach with a sand, gravel, and cobble substrate located between the upstream model boundary at river kilometer 245.9 and the upstream end of Ambush Rock at river kilometer 244.6, and a meandering reach with a predominately sand substrate located between upstream end of Ambush Rock and the downstream model boundary at river kilometer 228.4. Model cell size in the x and y (horizontal) dimensions is 5 meters by 5 meters along the computational grid centerline with 15 nodes in the z (vertical) dimension. The model was calibrated to historical streamflows evenly distributed between 141.6 and 2,548.9 cubic meters per second. The model was validated by comparing simulated velocities with velocities measured at 15 cross sections during steady streamflow. These 15 cross sections were each measured multiple (7-13) times to obtain velocities suitable for comparison to the model results. Comparison of modeled and measured velocities suggests that the model does a good job of reproducing flow patterns in the river, although some discrepancies were noted. The model was used to simulate water-surface elevation, depth, velocity, bed shear stress, and sediment mobility for Kootenai River streamflows of 170, 566, 1,130, 1,700, and 2,270 cubic meters per second (6,000, 20,000, 40,000, 60,000, and 80,000 cubic feet per second). The three lowest streamflow simulations represent a range of typical river conditions before and since the construction of Libby Dam, and the highest streamflow simulation (2,270 cubic meters per second) is approximately equal to the annual median peak streamflow prior to emplacement of Libby Dam in 1972. Streamflow greater than 566 cubic meters per second were incrementally increased by 570 cubic meters per second. For each

Scientific Investigations Report↗