First posted September 16, 2022 For additional information, contact: Southwest Biological Science CenterU.S. Geological Survey2255 N. Gemini DriveFlagstaff, AZ 86001 Sandbars are an important resource in the Colorado River corridor in Marble and Grand Canyons, Arizona, downstream from Glen Canyon Dam. Sandbars provide aquatic and riparian habitat and are used as campsites by river runners and hikers. The study area is the Colorado River between Glen Canyon Dam and Diamond Creek, which is about 388 kilometers (241 miles) downstream from the dam. Closure of Glen Canyon Dam in 1963 and subsequent flow regulation reduced the sediment supply, limited the magnitude and frequency of floods, and increased the magnitude of baseflows. The result has been widespread erosion of sandbars and expansion of native and non-native vegetation on previously bare sand deposits in this debris-fan dominated canyon river. This study reports on the on-going long-term measurement program of Northern Arizona University, initiated in 1990 with the Bureau of Reclamation, and now also with the U.S. Geological Survey’s Grand Canyon Monitoring and Research Center. We report on all sandbar measurements made between 1990 and 2020 to demonstrate the multi-decadal response of the sandbar monitoring sites resulting from flow regulation by Glen Canyon Dam. Because only one study site is located in Glen Canyon, the 25 kilometer (15.5 miles) reach just below Glen Canyon Dam, analyses of sandbar response are only made for the next two canyon segments in the down-river direction, Marble Canyon (388 kilometers [99 miles]) and Grand Canyon (265 kilometers [165 miles]), respectively, where the majority of study sites are located.We show that a majority of monitoring sites increased in volume during a period of frequent controlled floods intended to rebuild sandbars. In the period from 2004 to 2020, which included seven controlled floods, a median discharge of 350 cubic meters per second (m3/s), and greater than average tributary sand inputs in more than half of the years, net deposition occurred at 86 percent of long-term monitoring sites. This period was preceded by a period of net erosion (1990–2003) when there was one controlled flood greater than the nominal powerplant capacity of 940 m3/s. During this period the median discharge from Glen Canyon Dam was 376 m3/s and greater than average sand inputs occurred in only 36 percent of those years. At the end of the monitoring period in 2020, 61 percent of the study sites measured since 1990 underwent a net increase in sand volume. For the entire 31-year period, these trends were statistically significant for all six sandbar types studied, indicating that increased frequency of controlled flooding maintained sandbar volume at the majority of sites monitored. These floods, also referred to as high-flow experiments (HFEs), are part of a decision-making protocol approved in 2012 for coordinating dam releases timed to occur following large sand inputs to the Colorado River by a major tributary.These findings are based on digital elevation models (DEMs) derived from approximately (~)1,800 repeat surveys of sandbar and channel bed topography made annually, or more frequently, at the 45 long-term monitoring sites, of which 31 have been monitored since 1990 and 14 were added between 1990 and 2008. This large collection of monitoring sites comprises just 7 to 9 percent of all sandbars in Marble and Grand Canyons, respectively. Nevertheless, when compared with measurements of a larger sample, these sites provide consistent characterization of average sandbar response, despite the local variability in channel and debris fan geometry. We use sand volume and normalized sand volume for tracking geomorphic changes of sandbars, because these metrics are sensitive to both changes in sandbar area and sandbar elevation. Based on checkpoint comparisons and repeat measurements, DEM elevation uncertainty was determined to be ±0.05 meter (m) and this uncertainty was used in a spatially uniform estimate of volume uncertainty. We find that the magnitudes of the topographic changes were substantially greater than the measurement uncertainty.Sandbars of similar type throughout both Marble and Grand Canyons have responded similarly during the period of the HFE protocol, despite variations in sand supply and longitudinal extent of those inputs. It should be noted that tributary-supplied sand to Glen Canyon is negligible, much of the riverbed is now armored with cobbles, and the channel bed degradation is irreversible in the current flow and sediment supply regime. Because all these HFEs have been conducted during periods of sediment enrichment, other factors such as vegetation and geomorphic setting are likely the primary causes of variation among the monitoring sites. A larger percentage of the sandbar population, predominantly located in narrow reaches where stage changes are greater, is composed of sandbar types that remain dynamic and consistently aggrade during HFEs. In contrast, wide reaches of the river corridor where stage change is not as great are characterized by sandbars that have been stabilized by vegetation and progressive aggradation during floods. In the former case, a majority of sandbars are likely to remain dynamic, requiring continued use of HFEs to achieve desired management goals. In the latter case, HFEs can do no better than replace the sediment eroded during normal dam operation between high-flow events, as they become less effective because of a diminishing amount of accommodation space available for deposition. Long-term sandbar trajectory and the continued effectiveness of HFEs are related to the differential vegetation establishment at each bar type. Future sandbar monitoring may need to consider the effects of riparian vegetation removal.

Strandlines of peak-stage indicators (such as driftwood logs, woody debris, and trash) provide valuable data for understanding the maximum stage and extent of inundation during floods. A series of seven strandlines have been preserved along the Colorado River in Grand Canyon National Park, Arizona, USA. A survey and analysis of these strandlines was completed from the Colorado River at Lees Ferry, Ariz., gaging station to the Colorado River near Grand Canyon, Ariz., gaging station. Owing to the longitudinally discontinuous nature of the strandlines, several lines of evidence were used to determine the year of the flood associated with each strandline segment. This evidence included strandline relative vertical position, degree of peak-stage indicator weathering, datable trash drift, and map-view location. The seven distinct strandlines identified were deposited during floods with the following peak discharges (in cubic feet per second [ft 3 /s]) at the Colorado River at Lees Ferry, Ariz., gaging station (year of flood in parentheses): 210,000 ft 3 /s (1884), 170,000 ft 3 /s (1921), 125,000 ft 3 /s (1957), 108,000 ft 3 /s (1958), 97,000 ft 3 /s (1983), 52,500 ft 3 /s (1986), and 45,000 ft 3 /s (multiple events between 1996 and 2012). Stage-discharge relations were developed in areas where all, or most of the strandlines were present, and were compared to predicted stage-discharge relations from a one-dimensional flow model. River width exerted a strong control on these relations, with much greater stage change occurring for a given discharge change in narrower bedrock-dominated reaches than in wider reaches with more extensive channel-margin alluvium. This comprehensive dataset allows for the verification of model-predicted flood stage along the Colorado River in Grand Canyon National Park.

Effects of riparian vegetation on fluvial sediment dynamics depend on morphological traits of the constituent species. Determining the effects of different morphological guilds on sedimentation rates, as influenced by multiple aspects of dam operations, can help identify viable strategies for streamflow and vegetation management to achieve riparian resource goals. Plants of increasing size and branching density or complexity have been found to have greater effects on sedimentation in free-flowing systems; however, this relationship could differ in regulated rivers. We tested the hypothesis that plant guilds of increasing height and branching complexity would be positively associated with sedimentation rates on 23 sandbars deposited in zones of recirculating flow (eddies) along the Colorado River in Grand Canyon. We used an image-based vegetation classification and digital elevation models from annual topographic surveys to track associations between six plant morphological guilds and topographic change over 5 years. Vegetation had significant associations with deposition after accounting for geomorphic setting, but the ordinal guild scale was not positively correlated with deposition magnitude. Instead, low-statured rhizomatous and herbaceous guilds were particularly effective at capturing sediment in the separation zone of sandbars, whereas tall herbs and large shrubs were most effective at capturing sediment in reattachment zones. These nuanced interactions between geomorphic position and morphological guild may be a direct consequence of flow regulation through modifications to physical deposition and erosion processes. Flow regulation may also select for a narrow subset of morphological guilds, reducing the diversity of vegetation feedbacks on sedimentation and emphasizing geomorphic drivers.

Morphological change in river channels is frequently evaluated in the context of mass balance sediment budgets. In a closed sediment budget, measurements of sediment influx and efflux are coupled with measured changes in channel topography to provide both spatial and temporal resolution, and independent estimates of the mass balance. For sediment budgets constructed over long river segments (~102 channel widths or greater) and long periods (~2 years or longer), spatial and temporal accumulation of measurement uncertainty, compounded by inadequate sampling frequency or spatial coverage, may produce indeterminate results. The degree of indeterminacy may be evaluated in the context of a signal‐to‐noise ratio (SNR), which is a function of the magnitude of the mass balance and the magnitudes of potential systematic uncertainties associated with measurements and incomplete sampling. We report on a closed sand budget consisting of measurements of flux and two morphological surveys for a 50‐km segment of a large river over a 3‐year period. Accurate reporting of the magnitude and sign of the change in sand storage was only possible by using state‐of‐the‐art techniques with high temporal frequency and large spatial extent. Together, a sand flux and morphological mass balance revealed that sand evacuation was temporally concentrated (~100% of mass change occurred during 19% of the study period) and highly localized (70% of mass change occurred in 12% of the study segment). A SNR analysis revealed that uncertainty resulting from undersampling may approach or exceed that caused by measurement uncertainty and that daily sampling of suspended‐sand concentration or repeat mapping of at least 50% of the river segment was required to determine the sand budget with SNR > 1. The approach used here to analyze sand budget uncertainty is especially applicable to other river systems with large temporal variability in sediment transport and large spatial variability in erosion and deposition. © 2018 John Wiley & Sons, Ltd.

First posted July 30, 2018 For additional information, contact: SBSC Staff,Southwest Biological Science CenterU.S. Geological Survey2255 N. Gemini DriveFlagstaff, AZ 86001 Sandbars along the Colorado River are used as campsites by river runners and hikers and are an important recreational resource within Grand Canyon National Park, Arizona. Regulation of the flow of river water through Glen Canyon Dam has reduced the amount of sediment available to be deposited as sandbars, has reduced the magnitude and frequency of flooding events, and has increased the magnitude of baseflows. This has caused widespread erosion of sandbars and has allowed native and non-native vegetation to expand on open sand. Previous studies show an overall decline in campsite area despite the use of controlled floods to rebuild sandbars. Monitoring of campsites since 1998 has shown changes in campsite area, but the factors that cause gains and losses in campsite area have not been quantified. These factors include, among others, changes in sandbar volume and slope under different dam flow regimes that include controlled floods, gullying caused by monsoonal rains, vegetation expansion, and reworking of sediment by aeolian processes.Using 4-band aerial imagery and digital elevation models (DEMs) derived from total-station survey data, we analyzed topographic and vegetation change at 35 of 37 long-term monitoring sites (2 sites were excluded because topographic measurements do not overlap with measurements of campsite area) using data collected between 2002 and 2009 to quantify the factors affecting the size of campsite area. Over the course of the study period, there was a net loss in campsite area of 2,431 square meters (m2). We find that (1) 53 percent of the net loss was caused by topographic change associated with controlled floods and erosion of those flood deposits, (2) 47 percent of the net loss was caused by increases in vegetation cover, the majority of which occurred in high-elevation campsite area, and (3) gullying was significant at certain sites but overall was a minor factor.Sites in critical reaches—sections of river where campsites are infrequent or where there is high demand by river runners—were subjected to more erosion and changes in sandbar slope than sites in noncritical reaches, suggesting that campsite area is less stable in those reaches. There was also a greater increase in vegetation cover at sites in noncritical reaches than at sites in critical reaches. Our results show a continuation of sandbar erosion and vegetation encroachment that has been occurring at campsites since construction of the dam.A new campsite survey methodology using a tablet-based geographic information system (GIS) approach was also developed in an effort to map campsite area on digital orthophotographs. Using a series of repeat measurements, we evaluated the inherent uncertainty in mapping campsite area, the accuracy of the new tablet-based method, and if there is any bias between the tablet method and the total-station method that is currently used. We find that uncertainty associated with surveyor judgment while using the total-station method is about 15 percent, which is higher than a previously reported uncertainty of 10 percent. Use of the tablet method adds additional uncertainty; however, the benefits of being able to quantify factors that lead to campsite-area change in the field may outweigh the additional error. Future campsite monitoring may need to consist of a combination of total-station and orthophotograph techniques.

Sandbars are iconic features of the Colorado River in the Grand Canyon, Arizona, U.S.A. Following completion of Glen Canyon Dam in 1963, sediment deficit conditions caused erosion of eddy sandbars throughout much of the 360km study reach downstream from the dam. Controlled floods in 1996, 2004, and 2008 demonstrated that sand on the channel bed could be redistributed to higher elevations, and that floods timed to follow tributary sediment inputs would increase suspended sand concentrations during floods. Since 2012, a new management protocol has resulted in four controlled floods timed to follow large inputs of sand from a major tributary. Monitoring of 44 downstream eddy sandbars, initiated in 1990, shows that each controlled flood deposited significant amounts of sand and increased the size of subaerial sandbars. However, the magnitude of sandbar deposition varied from eddy to eddy, even over relatively short distances where main-stem suspended sediment concentrations were similar. Here, we characterize spatial and temporal trends in sandbar volume and site-scale (i.e., individual eddy) sediment storage as a function of flow, channel, and vegetation characteristics that reflect the reach-scale (i.e., kilometer-scale) hydraulic environment. We grouped the long-term monitoring sites based on geomorphic setting and used a principal component analysis (PCA) to correlate differences in sandbar behavior to changes in reach-scale geomorphic metrics. Sites in narrow reaches are less-vegetated, stage changes markedly with discharge, sandbars tend to remain dynamic, and sand storage change dominantly occurs in the eddy compared to the main channel. In wider reaches, where stage-change during floods may be half that of narrow sites, sandbars are more likely to be stabilized by vegetation, and floods tend to aggrade the vegetated sandbar surfaces. In these locations, deposition during controlled floods is more akin to floodplain sedimentation, and the elevation of sandbar surfaces increases with successive floods. Because many sandbars are intermediate to the end members described above, high-elevation bar surfaces stabilized by vegetation often have a more dynamic unvegetated sandbar on the channel-ward margin that aggrades and erodes in response to controlled flood cycles. Ultimately, controlled floods have been effective at increasing averaged sandbar volumes, and, while bar deposition during floods decreases through time where vegetation has stabilized sandbars, future controlled floods are likely to continue to result in deposition in a majority of the river corridor.

First posted March 23, 2017 For additional information, contact: GCMRC Staff, Southwest Biological Science Center U.S. Geological Survey Grand Canyon Monitoring and Research Center 2255 N. Gemini Drive Flagstaff, Arizona 86001 https://www.gcmrc.gov/ Bathymetric, topographic, and grain-size data were collected in May 2009 along a 33-mi reach of the Colorado River in Grand Canyon National Park, Arizona. The study reach is located from river miles 29 to 62 at the confluence of the Colorado and Little Colorado Rivers. Channel bathymetry was mapped using multibeam and singlebeam echosounders, subaerial topography was mapped using ground-based total-stations, and bed-sediment grain-size data were collected using an underwater digital microscope system. These data were combined to produce digital elevation models, spatially variable estimates of digital elevation model uncertainty, georeferenced grain-size data, and bed-sediment distribution maps. This project is a component of a larger effort to monitor the status and trends of sand storage along the Colorado River in Grand Canyon National Park. This report documents the survey methods and post-processing procedures, digital elevation model production and uncertainty assessment, and procedures for bed-sediment classification, and presents the datasets resulting from this study.

Many archeological sites in the Grand Canyon are being impacted by gully incision. In March 2008, a high-flow experiment (2008 HFE) was conducted with the intention of redistributing fine sediment (sand, silt, and clay) from the bed of the Colorado River to higher elevations along the channel margin. Deposition of fine sediment in gully mouths has been hypothesized to slow gully erosion rates and lessen impacts to archeological sites. The effects of the 2008 HFE on gullies were evaluated by comparing the topographic changes of three gullies at two study sites before and after the 2008 HFE. Comparison results indicated that sediment was deposited in gully mouths during the 2008 HFE, and that the inundated areas nearest to the river can be extensively altered by mainstream flow during high-flow events. Additionally, the history of gully evolution at the two study sites was examined between 1996 and 2010 and indicated that gullies have been subjected to thalweg incision and gully widening processes over a decadal timescale. Although the small sample size precludes extrapolating the results to other gullies, the findings contribute to the understanding of gully erosion in archeologically significant areas and have implications for future monitoring of gully erosion and evaluating the effectiveness of check dams intended to mitigate that erosion at archaeological sites in the Grand Canyon National Park.

Digital elevation models (DEMs) of eleven 2–5 kilometer reaches of the Colorado River ecosystem (CRE) in Grand Canyon were constructed from repeat bathymetric and topographic surveys collected between August 2000 and December 2004. The DEMs will be used by researchers to study the effects of Glen Canyon Dam (GCD) operations on the sediment resources of the CRE in Grand Canyon by quantifying morphological changes and sediment transfer within and among the study reaches. Airborne surveys collected light detection and ranging (lidar) and photogrammetric data, whereas ground topographic and bathymetric data were collected simultaneously on river trips. Surveys were conducted in August 2000, September 2000, May 2002, May 2004, November 2004, and December 2004. The aerial lidar and photogrammetric data were merged with the ground topographic and bathymetric data to create DEMs of the study areas with a grid resolution of 1 meter. For each survey period, the vertical component of uncertainty (specifically, reproducibility or precision) was estimated for each data type (lidar/photogrammetry, ground surveys, bathymetry) and for two different types of bed-surface texture (smooth and rough). The resulting DEMs from this study are a valuable contribution to ongoing efforts in assessing the effects of GCD operations on the CRE. The DEMs can be used to map the spatial characteristics of geomorphic change within the study reaches and to estimate sediment budgets for different time periods by calculating the difference in sediment volume between surveys. In addition, the DEMs provide essential boundary conditions for numerical models of sediment transport and deposition, as well as help define the spatial distribution of habitat for fisheries investigations.

In 2011, a large magnitude flow release from Flaming Gorge Reservoir, Wyoming and Utah, occurred in response to high snowpack in the middle Rocky Mountains. This was the third highest recorded discharge along the Green River downstream of Flaming Gorge Dam, Utah, since its initial closure in November 1962 and motivated a research effort to document effects of these flows on channel morphology and sedimentology at four long-term monitoring sites within the Canyon of Lodore in Dinosaur National Monument, Colorado and Utah. Data collected in September 2011 included raft-based bathymetric surveys, ground-based surveys of banks, channel cross sections and vegetation-plot locations, sand-bar stratigraphy, and painted rock recovery on gravel bars. As part of this surveying effort, Global Navigation Satellite System (GNSS) data were collected at benchmarks on the canyon rim and along the river corridor to establish a high-resolution survey control network. This survey control network allows for the collection of repeatable spatial and elevation data necessary for high accuracy geomorphic change detection. Nearly 10,000 ground survey points and more than 20,000 bathymetric points (at 1-meter resolution) were collected over a 5-day field campaign, allowing for the construction of reach-scale digital elevation models (DEMs). Additionally, we evaluated long-term geomorphic change at these sites using repeat topographic surveys of eight monumented cross sections at each of the four sites. Analysis of DEMs and channel cross sections show a spatially variable pattern of erosion and deposition, both within and between reaches. As much as 5 meters of scour occurred in pools downstream from flow constrictions, especially in channel segments where gravel bars were absent. By contrast, some channel cross sections were stable during the 2011 floods, and have shown almost no change in over a decade of monitoring. Partial mobility of gravel bars occurred, and although in some locations vegetation such as tamarisk (Tamarix ramosissima) was damaged, wholesale bed motion necessary to fully clear these surfaces was not evident. In flow recirculation zones, eddy sandbars aggraded one meter or more, increasing the area of bars exposed during typical dam operations. Yet overall, the 2011 flood resulted in a decrease in reach-scale sand storage because bed degradation exceeded bar deposition. The 2011 response is consistent with that of a similar event in 1999, which was followed by sand-bar erosion and sediment accumulation on the bed during subsequent years of normal dam operational flows. Although the 1999 and 2011 floods were exceptional in the post-dam system, they did not exceed the pre-dam 2-year flood, isolating their effects to the modern active channel with minor erosion or reworking of pre-dam deposits stabilized through vegetation encroachment.

Prior to the construction of large dams on the Green and Colorado Rivers, annual floods aggraded sandbars in lateral flow-recirculation eddies with fine sediment scoured from the bed and delivered from upstream. Flows greater than normal dam operations may be used to mimic this process in an attempt to increase time-averaged sandbar size. These controlled floods may rebuild sandbars, but sediment deficit conditions downstream from the dams restrict the frequency that controlled floods produce beneficial results. Here, we integrate complimentary, long-term monitoring data sets from the Colorado River in Marble and Grand Canyons downstream from Glen Canyon dam and the Green River in the Canyon of Lodore downstream from Flaming Gorge dam. Since the mid-1990s, several controlled floods have occurred in these canyon rivers. These controlled floods scour fine sediment from the bed and build sandbars in eddies, thus increasing channel relief. These changes are short-lived, however, as interflood dam operations erode sandbars within several months to years. Controlled flood response and interflood changes in bed elevation are more variable in Marble Canyon and Grand Canyon, likely reflecting more variable fine sediment supply and stronger transience in channel bed sediment storage. Despite these differences, neither system shows a trend in fine-sediment storage during the period in which controlled floods were monitored. These results demonstrate that controlled floods build eddy sandbars and increase channel relief for short interflood periods, and this response may be typical in other dam-influenced canyon rivers. The degree to which these features persist depends on the frequency of controlled floods, but careful consideration of sediment supply is necessary to avoid increasing the long-term sediment deficit.

Measurements of morphologic change are often used to infer sediment mass balance. Such measurements may, however, result in gross errors when morphologic changes over short reaches are extrapolated to predict changes in sediment mass balance for long river segments. This issue is investigated by examination of morphologic change and sediment influx and efflux for a 100 km segment of the Colorado River in Grand Canyon, Arizona. For each of four monitoring intervals within a 7 year study period, the direction of sand‐storage response within short morphologic monitoring reaches was consistent with the flux‐based sand mass balance. Both budgeting methods indicate that sand storage was stable or increased during the 7 year period. Extrapolation of the morphologic measurements outside the monitoring reaches does not, however, provide a reasonable estimate of the magnitude of sand‐storage change for the 100 km study area. Extrapolation results in large errors, because there is large local variation in site behavior driven by interactions between the flow and local bed topography. During the same flow regime and reach‐average sediment supply, some locations accumulate sand while others evacuate sand. The interaction of local hydraulics with local channel geometry exerts more control on local morphodynamic response than sand supply over an encompassing river segment. Changes in the upstream supply of sand modify bed responses but typically do not completely offset the effect of local hydraulics. Thus, accurate sediment budgets for long river segments inferred from reach‐scale morphologic measurements must incorporate the effect of local hydraulics in a sampling design or avoid extrapolation altogether.