Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical 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 1,009 records · Page 56Linked to original sources
USGS2021

The seismo-acoustics of submarine volcanic eruptions

Many of the world’s volcanoes are hidden beneath the ocean’s surface where eruptions are difficult to observe. However, seismo‐acoustic signals produced by these eruptions provide a useful means of identifying active submarine volcanism. A literature survey revealed reports of 119 seismo‐acoustically recorded submarine eruptions since 1939. Submarine eruptions have been recorded in all major tectonic settings, with a range of geochemistries, and at a variety of water depths, but the reports are dominated by eruptions in the Pacific and at only a few locations. Many of the reports offer little detail, with over half of the observations made from distances >500 km, and only about half were confirmed as eruptions by non‐seismo‐acoustic evidence. The reported seismo‐acoustic signals cover a wide variety of processes, including earthquakes, explosions, various types of tremor, signals related to lava extrusion, and landslides. Recorded signals can sometimes be difficult to classify or confidently associate with an eruption, although there has been progress in this regard. Real‐time monitoring of submarine eruptions has been on‐going for several decades on regional and global scales with growing interest and effort in local networks. Real‐time networks are complemented by short‐term instrument deployments that often give more detailed insights into the dynamics and processes of submarine eruptions. Thorough seismo‐acoustic monitoring and study has increased our understanding of submarine eruptions, especially of deep‐sea volcanoes and spreading centers. Despite this, there are still many outstanding questions that need to be addressed for submarine volcanoes to be as well understood and monitored as their terrestrial counterparts.

Journal of Geophysical Research↗
USGS2021

Limited mantle hydration by bending faults at the Middle America Trench

Seismic anisotropy measurements show that upper mantle hydration at the Middle America Trench (MAT) is limited to serpentinization and/or water in fault zones, rather than distributed uniformly. Subduction of hydrated oceanic lithosphere recycles water back into the deep mantle, drives arc volcanism, and affects seismicity at subduction zones. Constraining the extent of upper mantle hydration is an important part of understanding many fundamental processes on Earth. Substantially reduced seismic velocities in tomography suggest that outer rise plate‐bending faults provide a pathway for seawater to rehydrate the slab mantle just prior to subduction. Estimates of outer‐rise hydration based on tomograms vary significantly, with some large enough to imply that, globally, subduction has consumed more than two oceans worth of water during the Phanerozoic. We found that, while the mean upper mantle wavespeed is reduced at the MAT outer rise, the amplitude and orientation of inherited anisotropy are preserved at depths >1 km below the Moho. At shallower depths, relict anisotropy is replaced by slowing in the fault‐normal direction. These observations are incompatible with pervasive hydration but consistent with models of wave propagation through serpentinized fault zones that thin to <100‐m in width at depths >1 km below Moho. Confining hydration to fault zones reduces water storage estimates for the MAT upper mantle from ∼3.5 wt% to <0.9 wt% H 2 0. Since the intermediate thermal structure in the ∼24 Myr‐old MAT slab favors serpentinization, limited hydration suggests that fault mechanics are the limiting factor, not temperatures. Subducting mantle may be similarly dry globally.

Journal of Geophysical Research↗
USGS2021

B-positive: A robust estimator of aftershock magnitude distribution in transiently incomplete catalogs

The earthquake magnitude-frequency distribution is characterized by the b -value, which describes the relative frequency of large versus small earthquakes. It has been suggested that changes in b -value after an earthquake can be used to discriminate whether that earthquake is part of a foreshock sequence or a more typical mainshock-aftershock sequence, with a decrease in b -value heralding a larger earthquake to come. However, the measurement of b -value during an active aftershock sequence is strongly biased by short-term incompleteness of the earthquake catalog and by data-windowing, and these biases have the same direction as the proposed signal. Here I develop a new estimator of the b- value that is insensitive to transient changes in catalog completeness and that does not require data windowing. The new estimator “ b -positive” is based on the positive-only subset of the differences in magnitude between successive earthquakes, which are described by a double-exponential (Laplace) distribution with the same b -value as the magnitude distribution itself. The b -positive estimator greatly improves the robustness of continuous b -value measurements during active earthquake sequences, as well as in historical catalogs with unknown or variable completeness. The new estimator confirms some of the observations of Gulia and Wiemer (2019), although at a reduced level, showing a decrease and recovery of the b -value during several recent foreshock sequences that cannot be attributed simply to measurement bias. However, the unbiased b -value changes may be too subtle to use in a real-time earthquake alarm system.

Journal of Geophysical Research: Solid Earth↗
USGS2021

Heat flux from a vapor-dominated hydrothermal field beneath Yellowstone Lake

We report results from 149 heat flux measurements made over n ∼2-year interval at sites in and around a vapor-dominated geothermal field located at water depths of ∼100–120 m in Yellowstone Lake, Wyoming. Measurements of both in situ temperature and thermal conductivity as a function of depth were made with a 1 m probe via a remotely operated vehicle, and are combined to compute the vertical conductive heat flux. Inside the ∼55.5 × 10 3 m 2 bathymetric depression demarcating the vapor-dominated field, the median conductive flux is 13 W m −2 , with a conductive output of 0.72 MW. Outside the thermal field, the median conductive flux is 3.5 W m −2 . We observed 49 active vents inside the thermal field, with an estimated mass discharge rate of 56 kg s −1 , a median exit-fluid temperature of 132°C, and a total heat output of 29 MW. We find evidence for relatively weak secondary convection with a total output of 0.09 MW in thermal area lake floor sediments. Our data indicate that vapor beneath the thermal field is trapped by a low-permeability cap at a temperature of ∼189°C and a depth of ∼15 m below the lake floor. The thermal output of the Deep Hole is among the highest of any vapor-dominated field in Yellowstone, due in part to the high boiling temperatures associated with the elevated lake floor pressures.

Wyoming↗
USGS2021

Indicators of volcanic eruptions revealed by global M4+ earthquakes

Determining whether seismicity near volcanoes is due primarily to tectonic or magmatic processes is a challenging but critical endeavor for volcanic eruption forecasting and detection, especially at poorly monitored volcanoes. Global statistics on the occurrence and timing of earthquakes near volcanoes both within and outside of eruptive periods reveal patterns in eruptive seismicity that may improve our ability to discern magmatically driven seismicity from purely tectonic seismicity. In this paper, we catalog magnitude four and greater (M4+) earthquakes near volcanoes globally and compute statistics on their occurrence with respect to various eruptive and volcanic attributes, evaluating their utility as diagnostic indicators of eruptions. Using a 2‐week time window and a 30 km radius around the volcanoes, we find that 11% of eruptions are preceded by at least one M4+ earthquake, but only 1% of such earthquakes is followed by eruption. However, earthquakes located 5–15 km from the volcano, those with normal faulting mechanisms and/or large nondouble‐couple components, and those occurring as groups are more commonly associated with eruptions, providing significant forecasting utility in some cases. Similarly, certain volcanoes are more likely to exhibit such precursors, such as those with long repose periods. We illustrate the use of these data in eruption forecasting scenarios, including rapid identification of analogous earthquake sequences at other volcanoes. When integrated within the context of multiparametric, multidisciplinary probabilistic assessments of volcanic activity, global earthquake statistics can improve eruption forecasts, and our work provides a model for use on other rapidly expanding global volcanological databases.

Journal of Geophysical Research↗
USGS2021

Geometry of the décollement below eastern Bangladesh and implications for seismic hazard

Eastern Bangladesh sits on the seismically active Chittagong-Myanmar fold and thrust belt (CMFB), a north-trending accretionary wedge on the eastern side of the India-Eurasia collision. Earthquakes on the basal décollement and associated thrusts within the CMFB present a hazard to this densely populated region. In this study, we interpret 28 seismic reflection profiles from both published and unpublished sources to constrain the depth of the basal décollement. To convert profiles from the time domain to the depth domain, we integrate sonic log and seismic stacking velocity data to generate time-velocity relationships for different parts of the CMFB. Our analysis reveals that the décollement is ∼9 km deep in northeast and southeast Bangladesh, but shallows to ∼5 km in east-central Bangladesh. The décollement has an area of 7.25 × 10 4 km 2 (∼150 × 450 km), making it capable of an Mw 8.5 earthquake. However, the warped geometry of this fault might act as a rupture barrier were a large earthquake to occur on the décollement. Our combined velocity and fault model lay the groundwork for future studies to address seismic segmentation, ground shaking, and rupture modeling in the CMFB. Finally, we use our compiled data set to analyze the evolution of fold kinematics in the CMFB. We observe that folding style and failure mode varies, from mainly ductile deformation in the foreland to mainly brittle in the hinterland. The dual-failure modes within the CMFB support the hypothesis that a region with ductile deformation may still be capable of seismic behavior.

Journal of Geophysical Research: Solid Earth↗
USGS2021

Sources of volcanic tremor associated with the summit caldera collapse during the 2018 east rift eruption of Kīlauea Volcano, Hawai'i

Volcanic tremor occurring at the beginning of the 2018 Kīlauea eruption is characterized using both seismic and tilt data recorded at the Kīlauea summit. An automatic seismic network-based approach detects several types of tremor including (a) 0.5–1 Hz long-period tremor preceding the eruption, located at the south-southwest edge of Halema'uma'u Crater and attributed to the quasi-steady radiation from a shallow hydrothermal system and (b) two sequences of gliding tremor at the beginning of the eruption, both with locations on the edges of the crater and within it. The first sequence is attributed to two swarms of low-amplitude regularly repeating earthquakes induced by the jerky motions of a cylindrical rock piston with radius of 325 m, height of 250 m, and mass of 2.07 × 10 11 kg, progressively intruding 12.3 m into the shallow hydrothermal system with volume of 10 8 m 3 and depth extent of 300 m. The second sequence is attributed to a gradual evolution in the properties of a bubbly magma within an east-striking dike below Halema'uma'u Crater, impacted by repeated roof collapses. A fluid-filled crack model points to a decrease in gas volume fraction from 4.22% to 1.6 × 10 −2 % in the magma filling the dike, and a model of gas retro-diffusion within the melt suggests a two orders of magnitude decrease in bubble number density from 7 × 10 8 m −3 down to 4 × 10 6 m −3 . Both models feature a quasi to totally degassed magma by May 26.

Hawaii↗
USGS2021

Earthquake source mechanisms and stress field variations associated with wastewater-induced seismicity in southern Kansas, USA

The strong increase of seismicity rates in the contiguous USA over the last 10 years is linked to the injection of huge amounts of wastewater from oil and gas production in unconventional hydrocarbon reservoirs. We calculated 549 moment tensors of induced earthquakes ( M W ≤ 4.9) in southern Kansas to study their source mechanisms and their relation to injection activity. Seventeen percent of the events analyzed contained significant volumetric (ISO%) components, and these events mostly occurred near the two largest local earthquakes during the 4 months of largest active wastewater disposal. Mapping the local stress field, we determined that most of the region lies within a transtensional stress regime, with a maximum horizontal stress σ Hmax trending N75°E. In the epicentral area of the M W 4.9 Milan earthquake, the σ Hmax trend is rotated to about S80°E. Locally, two areas display a change in the stress field orientation with depth, from transtensional above 5.5 km depth to strike slip deeper in the basement. Relating the resolved fault geometries to the obtained local stress field orientation, we find that most of the activated fault planes were optimally oriented to the current stress field and thus small stress perturbations caused by the water injection could lead to failure.

Kansas↗
USGS2020

Using tracer variance decay to quantify variability of salinity mixing in the Hudson River Estuary

The salinity structure in an estuary is controlled by time‐dependent mixing processes. However, the locations and temporal variability of where significant mixing occurs is not well‐understood. Here we utilize a tracer variance approach to demonstrate the spatial and temporal structure of salinity mixing in the Hudson River Estuary. We run a 4‐month hydrodynamic simulation of the tides, currents, and salinity that captures the spring‐neap tidal variability as well as wind‐driven and freshwater flow events. On a spring‐neap time scale, salinity variance dissipation (mixing) occurs predominantly during the transition from neap to spring tides. On a tidal time scale, 60% of the salinity variance dissipation occurs during ebb tides and 40% during flood tides. Spatially, mixing during ebbs occurs primarily where lateral bottom salinity fronts intersect the bed at the transition from the main channel to adjacent shoals. During ebbs, these lateral fronts form seaward of constrictions located at multiple locations along the estuary. During floods, mixing is generated by a shear layer elevated in the water column at the top of the mixed bottom boundary layer, where variations in the along channel density gradients locally enhance the baroclinic pressure gradient leading to stronger vertical shear and more mixing. For both ebb and flood, the mixing occurs at the location of overlap of strong vertical stratification and eddy diffusivity, not at the maximum of either of those quantities. This understanding lends a new insight to the spatial and time dependence of the estuarine salinity structure.

Hudson River↗
USGS2021

Seasonal controls on sediment delivery and hydrodynamics in a vegetated tidally influenced interdistributary island

River deltas are maintained by a continuous supply of terrestrial sediments that provide critical land building material to help sustain and protect vulnerable ecological communities and serve as natural storm protection barriers. Local hydrodynamics are important in determining the degree to which fluvial sediments are removed from the water column and retained on the delta complex. During 2014, we measured hydrodynamics and sediment transport characteristics at one of the world's most rapidly prograding deltas, the Wax Lake delta in Louisiana, USA. We observed waves to be the dominant source of bottom stress for 70% of our observations. Sediment concentration tended to increase with shear stress, but only after stresses exceeded 0.01–0.02 Pa. Significant wave height and bottom stress were substantially reduced after June, when the emergence of American lotus ( Nelumbo lutea ) formed a dense canopy over the intertidal regions of the island splay. Hydrodynamics during these summer vegetated conditions were much more favorable to floc formation, and by extension particle settling, as shown by trends in the Kolmogorov microscale parameter over the course of the measurement campaign. Together, these findings suggest that the timing between peak river discharge and the emergence of vegetation may have a strong influence on rates of progradation in seasonally vegetated delta splays, whereby sediments delivered by flood events that extend late into summer may be governed by hydrodynamics that favor particle deposition, whereas those delivered prior to the summer may be more prone to remain in suspension and bypass the delta complex.

Louisiana↗
USGS2021

Persistent multidecadal variability since the 15th century in the southern Barents Sea derived from annually resolved shell-based records

In the North Atlantic Ocean, multidecadal variability in sea surface temperatures (SSTs) over the past several centuries has largely been inferred through terrestrial proxies and decadally resolved marine proxies. Annually resolved proxy records from marine archives provide valuable insight into this variability, but are especially rare from high latitude environments, particularly for centennial timescales. We constructed continuous, absolutely dated records of shell growth (1449–2014 CE; 564 years) and oxygen isotope ratios (δ 18 O shell ; 1539–2014 CE; 476 years) from shells of the bivalve Arctica islandica from coastal northern Norway, a location sensitive to large-scale North Atlantic Ocean dynamics. An annual (January–December) SST reconstruction derived from δ 18 O shell for the past five centuries suggests an increase of at least 2°C from the mid-18th century to 2014. The SST reconstruction correlates significantly with instrumental records and with other proxy reconstructions in the southern Barents Sea region. Spectral analysis of the shell growth and isotope records supports evidence for Atlantic multidecadal variability (65–80 year periodicity) extending into polar and subpolar latitudes for the past five centuries. These results provide additional evidence that multidecadal variability in SSTs are a persistent feature of the North Atlantic marine system.

Greenland Sea, Iceland Sea, Norwegian Sea, Barents↗
USGS2020

Ancient Martian aeolian sand dune deposits recorded in the stratigraphy of Valles Marineris and implications for past climates

Aeolian sediment transport, deposition, and erosion have been ongoing throughout Mars's history. This record of widespread aeolian processes is preserved in landforms and geologic units that retain important clues about past environmental conditions including wind patterns. In this study we describe landforms within Melas Chasma, Valles Marineris, that occur in distinct groups with linear to crescentic shapes, arranged with a characteristic wavelength; some possess slope profiles analogous to modern sand dunes yet show evidence for lithification. Based on the features' dimensions, asymmetry, and spatial patterns relative to modern equivalents, we interpret these landforms to be two classes of aeolian bedforms: decameter‐scale megaripples and sand dunes. The presence of superposed erosional features and depositional units indicates that these landforms were cemented and likely ancient. Melas paleodunes are found atop Hesperian‐aged layered deposits, but we estimate them to be younger, likely lithified in the Amazonian period. Although a range of degradation was observed, some paleodunes are >10 m tall and maintain steep lee sides (>25°), an uncommon scenario for terrestrial examples as other geologic processes lead to dune obliteration. The preserved paleobedform geometries are largely consistent with those of modern aeolian indicators, suggesting no major shifts in wind regime or contributing boundary conditions. Finally, we propose that their appearance and context require sequential periods of dune migration, stabilization following catastrophic burial, cementation, differential erosion, exposure, and burial. The presence of wholly preserved duneforms appears to be more common on Mars compared to the Earth and may signal something important about Martian landscape evolution.

Journal of Geophysical Research: Planets↗
USGS2020

Evidence for a diagenetic origin of Vera Rubin Ridge, Gale Crater, Mars: Summary and synthesis of Curiosity's exploration campaign

This paper provides an overview of the Curiosity rover's exploration at Vera Rubin ridge (VRR) and summarizes the science results. VRR is a distinct geomorphic feature on lower Aeolis Mons (informally known as Mount Sharp) that was identified in orbital data based on its distinct texture, topographic expression, and association with a hematite spectral signature. Curiosity conducted extensive remote sensing observations, acquired data on dozens of contact science targets, and drilled three outcrop samples from the ridge, as well as one outcrop sample immediately below the ridge. Our observations indicate that strata composing VRR were deposited in a predominantly lacustrine setting and are part of the Murray formation. The rocks within the ridge are chemically in family with underlying Murray formation strata. Red hematite is dispersed throughout much of the VRR bedrock, and this is the source of the orbital spectral detection. Gray hematite is also present in isolated, gray‐colored patches concentrated toward the upper elevations of VRR, and these gray patches also contain small, dark Fe‐rich nodules. We propose that VRR formed when diagenetic event(s) preferentially hardened rocks, which were subsequently eroded into a ridge by wind. Diagenesis also led to enhanced crystallization and/or cementation that deepened the ferric‐related spectral absorptions on the ridge, which helped make them readily distinguishable from orbit. Results add to existing evidence of protracted aqueous environments at Gale crater and give new insight into how diagenesis shaped Mars' rock record.

Journal of Geophysical Research: Planets↗
USGS2021

A rock record of complex aeolian bedforms in a Hesperian desert landscape: The Stimson formation as exposed in the Murray Buttes, Gale Crater, Mars

Lithified aeolian strata encode information about ancient planetary surface processes and the climate during deposition. Decoding these strata provides insight regarding past sediment transport processes, bedform kinematics, depositional landscape, and the prevailing climate. Deciphering these signatures requires a detailed analysis of sedimentary architecture to reconstruct dune morphology, motion, and the conditions that enabled their formation. Here, we show that a distinct sandstone unit exposed in the foothills of Mount Sharp, Gale crater, Mars, records the preserved expression of compound aeolian bedforms that accumulated in a large dune field. Analysis of Mastcam images of the Stimson formation shows that it consists of cross‐stratified sandstone beds separated by a hierarchy of erosive bounding surfaces formed during dune migration. The presence of two orders of surfaces with distinct geometrical relations reveals that the Stimson‐era landscape consisted of large dunes (draas) with smaller, superimposed dunes migrating across their lee slopes. Analysis of cross‐lamination and subset bounding surface geometries indicate a complex wind regime that transported sediment toward the north, constructing oblique dunes. This dune field was a direct product of the regional climate and the surface processes active in Gale crater during the fraction of the Hesperian Period recorded by the Stimson formation. The environment was arid, supporting a large aeolian dune field; this setting contrasts with earlier humid depositional episodes, recorded by the lacustrine sediments of the Murray formation (also Hesperian). Such fine‐scale reconstruction of landscapes on the ancient surface of Mars is important to understanding the planet’s past climate and habitability.

Journal of Geophysical Research - Planets↗
USGS2021

Mars: Abundant recurring slope lineae (RSL) following the planet-encircling dust event (PEDE) of 2018

Recurring slope lineae (RSL) are dark linear markings on Mars that regrow annually and likely originate from the flow of either liquid water or granular material. Following the great dust storm (or planet-encircling dust event, PEDE) of Mars year (MY) 34, Mars Reconnaissance Orbiter/High Resolution Imaging Science Experiment has seen many more candidate RSL than in typical Mars years. They have been imaged at more than 285 unique locations from August 2018 (when the atmosphere was clearing as the PEDE decayed) to August 2019, about half (157) of which are locations where RSL have not been documented previously. In MY34, 150 active RSL sites were identified in the southern middle latitudes (SML, -60° to -30°), whereas an average of 36 active sites were observed in each previous year (MY28–33). Post-PEDE RSL are also present during southern summer over a wider range of latitude, slope aspect, and L s (areocentric longitude of the sun) than in prior years. These RSL sites usually show evidence for recent dust deposition: obscuration of relatively dark areas, an overall brighter and redder surface than in prior years, and dust devil tracks, which indicate dust lifting by several mechanisms. We speculate that dust-lifting processes may initiate and sustain RSL activity. The RSL may form from flows of dust (perhaps clumped) and/or sand that is destabilized by dust movement or directly mobilized by dust devils. If this is the case, then the otherwise puzzling recurrence and year-to-year variability of RSL activity can be at least partly explained. The dust replenishment varies from year to year, which could explain interannual variations in RSL activity.

Journal of Geophysical Research: Planets↗
USGS2020

Compositional layering in Io driven by magmatic segregation and volcanism

The compositional evolution of volcanic bodies like Io is not well understood. Magmatic segregation and volcanic eruptions transport tidal heat from Io's interior to its surface. Several observed eruptions appear to be extremely high temperature (≥ 1600 K), suggesting either very high degrees of melting, refractory source regions, or intensive viscous heating on ascent. To address this ambiguity, we develop a model that couples crust and mantle dynamics to a simple compositional system. We analyse the model to investigate chemical structure and evolution. We demonstrate that magmatic segregation and volcanic eruptions lead to stratification of the mantle, the extent of which depends on how easily high temperature melts from the more refractory lower mantle can migrate upwards. We propose that Io's highest temperature eruptions originate from this lower mantle region, and that such eruptions act to limit the degree of compositional stratification.

Journal of Geophysical Research↗
USGS2021

Widespread exposures of extensive clean shallow ice in the mid-latitudes of Mars

Although ice in the Martian midlatitudes is typically covered by a layer of dust or regolith, it is exposed in some locations by fresh impact craters or in erosional scarps. In both cases, the exposed ice is massive or excess ice with a low lithic content. We find that erosional scarps occur between 50° and 61° north and south latitude and that they are concentrated in and near Milankovič crater in the northern hemisphere and southeast of the Hellas basin in the southern hemisphere. These may represent locations of particularly thick or clean bodies of ice. Pits created by retreat of the scarps represent sublimation-thermokarst landforms that evolve in a manner distinct from other ice-loss landforms on Mars. New impact craters reveal that clean subsurface ice is widespread at middle- and high-latitudes in both hemispheres at depths less than 1 m. Both the depth to ice and the ice content appear to exhibit significant variability over tens to hundreds of meters. The lowest-latitude exposed ice is near 39°N and is at the edge of a region where impact exposures between 40° and 50°N are common, consistent with other indications of a high ice content. This lowest-latitude ice may be currently unstable and subliming. Impact craters on lineated valley fill excavate ice blocks that may represent the top of debris-covered glacial ice. Together, these landforms indicate widespread, clean subsurface ice at middle-latitudes on Mars. The distribution and properties of this ice could provide information about past climate conditions.

Journal of Geophysical Research: Planets↗
USGS2021

Self-limitation of sand storage in a bedrock-canyon river arising from the interaction of flow and grain size

Bedrock-canyon rivers tend to be supply limited because they are efficient transporters of sediment and not because the upstream supply of sediment is small. A byproduct of this supply limitation is that the finer alluvium stored in these rivers has shorter residence times and smaller volumes than in alluvial rivers. To improve our understanding of disequilibrium sediment transport and its effect on sand storage in bedrock-canyon rivers, we undertook a 20-year study, synthesized herein, of the Colorado River in Grand Canyon. Despite the large loads for which it was renowned, this river exhibited evidence of natural sand-supply limitation and became the perfect natural laboratory for studying sand transport in a bedrock canyon after upstream dam construction exacerbated this supply limitation. During our study, we made and analyzed an unprecedented ∼2.5 million measurements of the suspended and bed sediment. Results indicate that sand storage in this bedrock-canyon river is self-limiting owing to the physical controls of flow and grain size causing negative feedbacks that likely also operate in other bedrock-canyon rivers. Following episodic tributary floods that supply finer sand, sand migrates quickly downstream in the form of a wave in which large systematic changes in bed-sand grain size occur. These grain-size changes cause discharge-independent systematic changes in suspended-sand concentration in excess of a factor of 20. Although the tributary supply of sand increases the amount of sand storage, it also greatly increases the downstream sand transport by causing bed-sand fining, thus limiting the residence time and volume of sand storage.

Arizona↗