SEARCH · USGS Science

Results for “Geophysical Journal”

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

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27Linked to original sources

USGS2014

Triggering and modulation of geyser eruptions in Yellowstone National Park by earthquakes, earth tides, and weather

We analyze intervals between eruptions (IBEs) data acquired between 2001 and 2011 at Daisy and Old Faithful geysers in Yellowstone National Park. We focus our statistical analysis on the response of these geysers to stress perturbations from within the solid earth (earthquakes and earth tides) and from weather (air pressure and temperature, precipitation, and wind). We conclude that (1) the IBEs of these geysers are insensitive to periodic stresses induced by solid earth tides and barometric pressure variations; (2) Daisy (pool geyser) IBEs lengthen by evaporation and heat loss in response to large wind storms and cold air; and (3) Old Faithful (cone geyser) IBEs are not modulated by air temperature and pressure variations, wind, and precipitation, suggesting that the subsurface water column is decoupled from the atmosphere. Dynamic stress changes of 0.1−0.2 MPa resulting from the 2002 M-7.9 Denali, Alaska, earthquake surface waves caused a statistically significant shortening of Daisy geyser's IBEs. Stresses induced by other large global earthquakes during the study period were at least an order of magnitude smaller. In contrast, dynamic stresses of >0.5 MPa from three large regional earthquakes in 1959, 1975, and 1983 caused lengthening of Old Faithful's IBEs. We infer that most subannual geyser IBE variability is dominated by internal processes and interaction with other geysers. The results of this study provide quantitative bounds on the sensitivity of hydrothermal systems to external stress perturbations and have implications for studying the triggering and modulation of volcanic eruptions by external forces.

Wyoming↗

USGS2014

Three-dimensional seismic velocity structure of Mauna Loa and Kilauea volcanoes in Hawaii from local seismic tomography

We present a new three-dimensional seismic velocity model of the crustal and upper mantle structure for Mauna Loa and Kilauea volcanoes in Hawaii. Our model is derived from the first-arrival times of the compressional and shear waves from about 53,000 events on and near the Island of Hawaii between 1992 and 2009 recorded by the Hawaiian Volcano Observatory stations. The V p model generally agrees with previous studies, showing high-velocity anomalies near the calderas and rift zones and low-velocity anomalies in the fault systems. The most significant difference from previous models is in V p / V s structure. The high- V p and high- V p / V s anomalies below Mauna Loa caldera are interpreted as mafic magmatic cumulates. The observed low- V p and high- V p / V s bodies in the Kaoiki seismic zone between 5 and 15 km depth are attributed to the underlying volcaniclastic sediments. The high- V p and moderate- to low- V p / V s anomalies beneath Kilauea caldera can be explained by a combination of different mafic compositions, likely to be olivine-rich gabbro and dunite. The systematically low- V p and low- V p / V s bodies in the southeast flank of Kilauea may be caused by the presence of volatiles. Another difference between this study and previous ones is the improved V p model resolution in deeper layers, owing to the inclusion of events with large epicentral distances. The new velocity model is used to relocate the seismicity of Mauna Loa and Kilauea for improved absolute locations and ultimately to develop a high-precision earthquake catalog using waveform cross-correlation data.

Hawaii↗

USGS2014

Suspended particulate layers and internal waves over the southern Monterey Bay continental shelf: an important control on shelf mud belts?

Physical and optical measurements taken over the mud belt on the southern continental shelf of Monterey Bay, California documented the frequent occurrence of suspended particulate matter features, the majority of which were detached from the seafloor, centered 9–33 m above the bed. In fall 2011, an automated profiling mooring and fixed instrumentation, including a thermistor chain and upward-looking acoustic Doppler current profiler, were deployed at 70 m depth for 5 weeks, and from 12 to 16 October a long-range autonomous underwater vehicle performed across-shelf transects. Individual SPM events were uncorrelated with local bed shear stress caused by surface waves and bottom currents. Nearly half of all observed SPM layers occurred during 1 week of the study, 9–16 October 2011, and were advected past the fixed profiling mooring by the onshore phase of semidiurnal internal tide bottom currents. At the start of the 9–16 October period, we observed intense near-bed vertical velocities capable of lifting particulates into the middle of the water column. This “updraft” event appears to have been associated with nonlinear adjustment of high-amplitude internal tides over the mid and outer shelf. These findings suggest that nonlinear internal tidal motions can erode material over the outer shelf and that, once suspended, this SPM can then be transported shoreward to the middle and shallow sections of the mud belt. This represents a fundamental broadening of our understanding of how shelf mud belts may be built up and sustained.

California↗

USGS2014

Interactions between waves, sediment, and turbulence on a shallow estuarine mudflat

Measurements were collected on a shallow estuarine mudflat in northern San Francisco Bay to examine the physical processes controlling waves, turbulence, sediment resuspension, and their interactions. Tides alone forced weak to moderate currents of 10–30 cm s -1 in depths of 0–3 m, and maintained a background suspension of 30–50 mg L21 of fine sediment. In the presence of wind waves, bottom orbital velocities spanned 20–30 cm s -1 , suspended-sediment concentrations (SSC) at 15 and 30 cm above the bed (cmab) increased by 1–2 orders of magnitude, and vertical gradients in SSC were strong enough to produce turbulence-limiting stratification, with gradient Richardson numbers exceeding 0.25. Simultaneously, turbulent stresses (decomposed from wave motions) increased by an order of magnitude. The apparent contradiction of energetic turbulence in the presence of strong stratification was reconciled by considering the turbulent kinetic energy (TKE) budget: in general, dissipation and buoyancy flux were balanced by local shear production, and each of these terms increased during wave events. The classic wave-current boundary layer model represented the observations qualitatively, but not quantitatively since the velocity profile could not be approximated as logarithmic. Rather, the mean shear was elevated by the Stokes drift return flow and wind-generated surface stress, which diffused sediment upward and limited stratification. Our findings highlight a pathway for waves to supply energy to both the production and destruction of turbulence, and demonstrate that in such shallow depths, TKE and SSC can be elevated over more of the water column than predicted by traditional models.

California↗

USGS2014

Testing the accuracy of a 1-D volcanic plume model in estimating mass eruption rate

During volcanic eruptions, empirical relationships are used to estimate mass eruption rate from plume height. Although simple, such relationships can be inaccurate and can underestimate rates in windy conditions. One-dimensional plume models can incorporate atmospheric conditions and give potentially more accurate estimates. Here I present a 1-D model for plumes in crosswind and simulate 25 historical eruptions where plume height H obs was well observed and mass eruption rate M obs could be calculated from mapped deposit mass and observed duration. The simulations considered wind, temperature, and phase changes of water. Atmospheric conditions were obtained from the National Center for Atmospheric Research Reanalysis 2.5° model. Simulations calculate the minimum, maximum, and average values ( M min , M max , and M avg ) that fit the plume height. Eruption rates were also estimated from the empirical formula M empir = 140 H obs 4.14 ( M empir is in kilogram per second, H obs is in kilometer). For these eruptions, the standard error of the residual in log space is about 0.53 for M avg and 0.50 for M empir . Thus, for this data set, the model is slightly less accurate at predicting M obs than the empirical curve. The inability of this model to improve eruption rate estimates may lie in the limited accuracy of even well-observed plume heights, inaccurate model formulation, or the fact that most eruptions examined were not highly influenced by wind. For the low, wind-blown plume of 14–18 April 2010 at Eyjafjallajökull, where an accurate plume height time series is available, modeled rates do agree better with M obs than M empir .

Journal of Geophysical Research D: Atmospheres↗

USGS2014

Igneous mineralogy at Bradbury Rise: the first ChemCam campaign at Gale crater

Textural and compositional analyses using ChemCam Remote Micro Imager (RMI) and Laser Induced Breakdown Spectroscopy (LIBS) have been performed on 5 float rocks and coarse gravels along the first 100 meters of the Curiosity traverse at Bradbury rise. ChemCam, the first LIBS instrument sent to another planet, offers the opportunity to assess mineralogic diversity at grain-size scales (∼100μm) and, from this, lithologic diversity. Depth profiling indicates that targets are relatively free of surface coatings. One type of igneous rock is volcanic and includes both aphanitic (Coronation) and porphyritic (Mara) samples. The porphyritic sample shows dark grains that are likely pyroxene megacrysts in a finegrained mesostasis containing andesine needles. Both types have magnesium-poor basaltic compositions and in this respect are similar to the evolved Jake-Matijevic rock [Stolper et al., 2013] analyzed further along the Curiosity traverse both with APXS and ChemCam instruments. The second rock-type encountered is a coarse-grained intrusive rock (Thor Lake) showing equigranular texture with mm size crystals of feldspars and Fe-Ti oxides. Such a rock is not unique at Gale as the surrounding coarse gravels (such as Beaulieu) and the conglomerate Link are dominated by feldspathic (andesine-bytownite) clasts. Finally, alkali feldspar compositions associated with a silica polymorph have been analyzed in fractured filling material of Preble rock and in Stark, a putative pumice or an impact melt. These observations document magmatic diversity at Gale and describe the first fragments of feldspar-rich lithologies (possibly an anorthosite) that may be ancient crust transported from the crater rim and now forming float rocks, coarse gravel or conglomerate clasts.

Journal of Geophysical Research E: Planets↗

USGS2014

HiRISE observations of new impact craters exposing Martian ground ice

Twenty small new impact craters or clusters have been observed to excavate bright material inferred to be ice at mid and high latitudes on Mars. In the northern hemisphere, the craters are widely distributed geographically and occur at latitudes as low as 39°N. Stability modeling suggests that this ice distribution requires a long-term average atmospheric water vapor content around 25 precipitable microns, more than double the present value, which is consistent with the expected effect of recent orbital variations. Alternatively, near-surface humidity could be higher than expected for current column abundances if water vapor is not well-mixed with atmospheric CO 2 , or the vapor pressure at the ice table could be lower due to salts. Ice in and around the craters remains visibly bright for months to years, indicating that it is clean ice rather than ice-cemented regolith. Although some clean ice may be produced by the impact process, it is likely that the original ground ice was excess ice (exceeding dry soil pore space) in many cases. Observations of the craters suggest small-scale heterogeneities in this excess ice. The origin of such ice is uncertain. Ice lens formation by migration of thin films of liquid is most consistent with local heterogeneity in ice content and common surface boulders, but in some cases nearby thermokarst landforms suggest large amounts of excess ice that may be best explained by a degraded ice sheet.

Journal of Geophysical Research E: Planets↗

USGS2014

Trace element geochemistry (Li, Ba, Sr, and Rb) using Curiosity's ChemCam: early results for Gale crater from Bradbury Landing Site to Rocknest

The ChemCam instrument package on the Mars rover, Curiosity, provides new capabilities to probe the abundances of certain trace elements in the rocks and soils on Mars using the laser-induced breakdown spectroscopy technique. We focus on detecting and quantifying Li, Ba, Rb, and Sr in targets analyzed during the first 100 sols, from Bradbury Landing Site to Rocknest. Univariate peak area models and multivariate partial least squares models are presented. Li, detected for the first time directly on Mars, is generally low (<15 ppm). The lack of soil enrichment in Li, which is highly fluid mobile, is consistent with limited influx of subsurface waters contributing to the upper soils. Localized enrichments of up to ~60 ppm Li have been observed in several rocks but the host mineral for Li is unclear. Bathurst_Inlet is a fine-grained bedrock unit in which several analysis locations show a decrease in Li and other alkalis with depth, which may imply that the unit has undergone low-level aqueous alteration that has preferentially drawn the alkalis to the surface. Ba (~1000 ppm) was detected in a buried pebble in the Akaitcho sand ripple and it appears to correlate with Si, Al, Na, and K, indicating a possible feldspathic composition. Rb and Sr are in the conglomerate Link at abundances >100 ppm and >1000 ppm, respectively. These analysis locations tend to have high Si and alkali abundances, consistent with a feldspar composition. Together, these trace element observations provide possible evidence of magma differentiation and aqueous alteration.

Journal of Geophysical Research E: Planets↗

USGS2014

Sands at Gusev Crater, Mars

Processes, environments, and the energy associated with the transport and deposition of sand at Gusev Crater are characterized at the microscopic scale through the comparison of statistical moments for particle size and shape distributions. Bivariate and factor analyses define distinct textural groups at 51 sites along the traverse completed by the Spirit rover as it crossed the plains and went into the Columbia Hills. Fine-to-medium sand is ubiquitous in ripples and wind drifts. Most distributions show excess fine material, consistent with a predominance of wind erosion over the last 3.8 billion years. Negative skewness at West Valley is explained by the removal of fine sand during active erosion, or alternatively, by excess accumulation of coarse sand from a local source. The coarse to very coarse sand particles of ripple armors in the basaltic plains have a unique combination of size and shape. Their distribution display significant changes in their statistical moments within the ~400 m that separate the Columbia Memorial Station from Bonneville Crater. Results are consistent with aeolian and/or impact deposition, while the elongated and rounded shape of the grains forming the ripples, as well as their direction of origin, could point to Ma'adim Vallis as a possible source. For smaller particles on the traverse, our findings confirm that aeolian processes have dominated over impact and other processes to produce sands with the observed size and shape patterns across a spectrum of geologic (e.g., ripples and plains soils) and aerographic settings (e.g., wind shadows).

Journal of Geophysical Research E: Planets↗

USGS2014

Volcanic glass signatures in spectroscopic survey of newly proposed lunar pyroclastic deposits

Moon Mineralogy Mapper spectroscopic observations are used to assess the mineralogy of five sites that have recently been proposed to include lunar dark mantle deposits (DMDs). Volcanic glasses have, for the first time, clearly been identified at the location of three of the proposed pyroclastic deposits. This is the first time that volcanic glasses have been identified at such a small scale on the lunar surface from remote sensing observations. Deposits at Birt E, Schluter, and Walther A appear to be glassy DMDs. Deposits at Birt E and Schluter show (1) morphological evidence suggesting a likely vent and (2) mineralogical evidence indicative of the presence of volcanic glasses. The Walther A deposits, although they show no morphological evidence of vents, have the spectroscopic characteristics diagnostic of volcanic glasses. The deposits of the Freundlich-Sharonov basin are separated in two areas: (1) the Buys-Ballot deposits lack mineralogical and morphological evidence and thus are found to be associated with mare volcanism not with DMDs and (2) the Anderson crater deposits, which do not exhibit glassy DMD signatures, but they appear to be associated with possible vent structures and so may be classifiable as DMDs. Finally, dark deposits near the crater Kopff are found to be associated with likely mare volcanism and not associated with DMDs. The spectral identification of volcanic glass seen in many of the potential DMDs is a strong indicator of their pyroclastic origin.

Journal of Geophysical Research E: Planets↗

USGS2014

The rock abrasion record at Gale Crater: Mars Science Laboratory results from Bradbury Landing to Rocknest

Ventifacts, rocks abraded by wind-borne particles, are found in Gale Crater, Mars. In the eastward drive from “Bradbury Landing” to “Rocknest,” they account for about half of the float and outcrop seen by Curiosity's cameras. Many are faceted and exhibit abrasion textures found at a range of scales, from submillimeter lineations to centimeter-scale facets, scallops, flutes, and grooves. The drive path geometry in the first 100 sols of the mission emphasized the identification of abrasion facets and textures formed by westerly flow. This upwind direction is inconsistent with predictions based on models and the orientation of regional dunes, suggesting that these ventifact features formed from very rare high-speed winds. The absence of active sand and evidence for deflation in the area indicates that most of the ventifacts are fossil features experiencing little abrasion today.

Journal of Geophysical Research E: Planets↗

USGS2014

Chemistry and texture of the rocks at Rocknest, Gale Crater: Evidence for sedimentary origin and diagenetic alteration

A suite of eight rocks analyzed by the Curiosity Rover while it was stopped at the Rocknest sand ripple shows the greatest chemical divergence of any potentially sedimentary rocks analyzed in the early part of the mission. Relative to average Martian soil and to the stratigraphically lower units encountered as part of the Yellowknife Bay formation, these rocks are significantly depleted in MgO, with a mean of 1.3 wt %, and high in Fe, averaging over 20 wt % FeOT, with values between 15 and 26 wt % FeOT. The variable iron and low magnesium and rock texture make it unlikely that these are igneous rocks. Rock surface textures range from rough to smooth, can be pitted or grooved, and show various degrees of wind erosion. Some rocks display poorly defined layering while others seem to show possible fractures. Narrow vertical voids are present in Rocknest 3, one of the rocks showing the strongest layering. Rocks in the vicinity of Rocknest may have undergone some diagenesis similar to other rocks in the Yellowknife Bay Formation as indicated by the presence of soluble calcium phases. The most reasonable scenario is that fine-grained sediments, potentially a mixture of feldspar-rich rocks from Bradbury Rise and normal Martian soil, were lithified together by an iron-rich cement.

Journal of Geophysical Research E: Planets↗

USGS2014

Linking channel hydrology with riparian wetland accretion in tidal rivers

The hydrologic processes by which tide affects river channel and riparian morphology within the tidal freshwater zone are poorly understood, yet are fundamental to predicting the fate of coastal rivers and wetlands as sea level rises. We investigated patterns of sediment accretion in riparian wetlands along the non-tidal through oligohaline portion of two coastal plain rivers in Maryland, U.S.A., and how flow velocity, water level, and suspended sediment concentration (SSC) in the channel may have contributed to those patterns. Sediment accretion was measured over a one year period using artificial marker horizons, channel hydrology was measured over a one month period using acoustic Doppler current profilers, and SSC was predicted from acoustic backscatter. Riparian sediment accretion was lowest at the non-tidal sites (mean and standard deviation = 8 ± 8 mm yr -1 ), highest at the upstream tidal freshwater forested wetlands (TFFW) (33 ± 28 mm yr -1 ), low at the midstream TFFW (12 ± 9 mm yr -1 ), and high at the oligohaline (fresh-to-brackish) marshes (19 ± 8 mm yr -1 ). Channel maximum flood and ebb velocity was 2-fold faster at the oligohaline than tidal freshwater zone on both tidal rivers, corresponding with the differences in in-channel SSC: the oligohaline zone's SSC was more than double the tidal freshwater zone's, and was greater than historical SSC at the non-tidal gages. The tidal wave characteristics differed between rivers, leading to significantly greater in-channel SSC during floodplain inundation in the weakly convergent than the strongly convergent tidal river. Overall sediment accretion was higher in the embayed river likely due to a single storm discharge and associated sedimentation.

Maryl↗

USGS2014

Modeling erosion of ice-rich permafrost bluffs along the Alaskan Beaufort Sea coast

The Arctic climate is changing, inducing accelerating retreat of ice-rich permafrost coastal bluﬀs. Along Alaska’s Beaufort Sea coast, erosion rates have increased roughly threefold from 6.8 to 19 m yr −1 since 1955 while the sea ice-free season has increased roughly twofold from 45 to 100 days since 1979. We develop a numerical model of bluﬀ retreat to assess the relative roles of the length of sea ice-free season, sea level, water temperature, nearshore waveﬁeld, and permafrost temperature in controlling erosion rates in this setting. The model captures the processes of erosion observed in short-term monitoring experiments along the Beaufort Sea coast, including evolution of melt notches, topple of ice wedge-bounded blocks, and degradation of these blocks. Model results agree with time-lapse imagery of bluﬀ evolution and time series of ocean-based instrumentation. Erosion is highly episodic with 40% of erosion is accomplished during less than 5% of the sea ice-free season. Among the formulations of the submarine erosion rate we assessed, we advocate those that employ both water temperature and nearshore waveﬁeld. As high water levels are a prerequisite for erosion, any future changes that increase the frequency with which water levels exceed the base of the bluﬀs will increase rates of coastal erosion. The certain increases in sea level and potential changes in storminess will both contribute to this eﬀect. As water temperature also inﬂuences erosion rates, any further expansion of the sea ice-free season into the midsummer period of greatest insolation is likely to result in an additional increase in coastal retreat rates.

Alaska↗

USGS2014

Seismicity within a propagating ice shelf rift: the relationship between icequake locations and ice shelf structure

Iceberg calving is a dominant mass loss mechanism for Antarctic ice shelves, second only to basal melting. An important known process involved in calving is the initiation and propagation of through-penetrating fractures called rifts; however, the mechanisms controlling rift propagation remain poorly understood. To investigate the mechanics of ice-shelf rifting, we analyzed seismicity associated with a propagating rift tip on the Amery Ice Shelf, using data collected during the Austral summers of 2004-2007. We investigated seismicity associated with fracture propagation using a suite of passive seismological techniques including icequake locations, back projection, and moment tensor inversion. We confirm previous results that show that seismicity is characterized by periods of relative quiescence punctuated by swarms of intense seismicity of one to three hours. However, even during periods of quiescence, we find significant seismic deformation around the rift tip. Moment tensors, calculated for a subset of the largest icequakes (M W > -2.0) located near the rift tip, show steeply dipping fault planes, horizontal or shallowly plunging stress orientations, and often have a significant volumetric component. They also reveal that much of the observed seismicity is limited to the upper 50 m of the ice shelf. This suggests a complex system of deformation that involves the propagating rift, the region behind the rift tip, and a system of rift-transverse crevasses. Small-scale variations in the mechanical structure of the ice shelf, especially rift-transverse crevasses and accreted marine ice, play an important role in modulating the rate and location of seismicity associated with propagating ice shelf rifts.

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Predictions of barrier island berm evolution in a time-varying storm climatology

Low-lying barrier islands are ubiquitous features of the world's coastlines, and the processes responsible for their formation, maintenance, and destruction are related to the evolution of smaller, superimposed features including sand dunes, beach berms, and sandbars. The barrier island and its superimposed features interact with oceanographic forces (e.g., overwash) and exchange sediment with each other and other parts of the barrier island system. These interactions are modulated by changes in storminess. An opportunity to study these interactions resulted from the placement and subsequent evolution of a 2 m high sand berm constructed along the northern Chandeleur Islands, LA. We show that observed berm length evolution is well predicted by a model that was fit to the observations by estimating two parameters describing the rate of berm length change. The model evaluates the probability and duration of berm overwash to predict episodic berm erosion. A constant berm length change rate is also predicted that persists even when there is no overwash. The analysis is extended to a 16 year time series that includes both intraannual and interannual variability of overwash events. This analysis predicts that as many as 10 or as few as 1 day of overwash conditions would be expected each year. And an increase in berm elevation from 2 m to 3.5 m above mean sea level would reduce the expected frequency of overwash events from 4 to just 0.5 event-days per year. This approach can be applied to understanding barrier island and berm evolution at other locations using past and future storm climatologies.

Louisiana↗

USGS2014

Transport of fine sediment over a coarse, immobile riverbed

Sediment transport in cobble-boulder rivers consists mostly of fine sediment moving over a coarse, immobile bed. Transport rate depends on several interrelated factors: boundary shear stress, the grain size and volume of fine sediment, and the configuration of fine sediment into interstitial deposits and bed forms. Existing models do not incorporate all of these factors. Approaches that partition stress face a daunting challenge because most of the boundary shear is exerted on immobile grains. We present an alternative approach that divides the bed into sand patches and interstitial deposits and is well constrained by two clear end-member cases: full sand cover and absence of sand. Entrainment from sand patches is a function of their aerial coverage. Entrainment from interstices among immobile grains is a function of sand elevation relative to the size of the immobile grains. The bed-sand coverage function is used to predict the ratio of the rate of entrainment from a partially covered bed to the rate of entrainment from a completely sand-covered bed, which is determined using a standard sand transport model. We implement the bed-sand coverage function in a morphodynamic routing model and test it against observations of sand bed elevation and suspended sand concentration for conditions of nonuniform fine sediment transport in a large flume with steady uniform flow over immobile hemispheres. The results suggest that this approach may provide a simple and robust method for predicting the transport and migration of fine sediment through rivers with coarse, immobile beds.

Journal of Geophysical Research F: Earth Surface↗

USGS2014

Rapid reservoir erosion, hyperconcentrated flow, and downstream deposition triggered by breaching of 38 m tall Condit Dam, White Salmon River, Washington

Condit Dam on the White Salmon River, Washington, a 38 m high dam impounding a large volume (1.8 million m 3 ) of fine-grained sediment (60% sand, 35% silt and clay, and 5% gravel), was rapidly breached in October 2011. This unique dam decommissioning produced dramatic upstream and downstream geomorphic responses in the hours and weeks following breaching. Blasting a 5 m wide hole into the base of the dam resulted in rapid reservoir drawdown, abruptly releasing ~1.6 million m 3 of reservoir water, exposing reservoir sediment to erosion, and triggering mass failures of the thickly accumulated reservoir sediment. Within 90 min of breaching, the reservoir's water and ~10% of its sediment had evacuated. At a gauging station 2.3 km downstream, flow increased briefly by 400 m 3  s −1 during passage of the initial pulse of released reservoir water, followed by a highly concentrated flow phase—up to 32% sediment by volume—as landslide-generated slurries from the reservoir moved downstream. This hyperconcentrated flow, analogous to those following volcanic eruptions or large landslides, draped the downstream river with predominantly fine sand. During the ensuing weeks, suspended-sediment concentration declined and sand and gravel bed load derived from continued reservoir erosion aggraded the channel by >1 m at the gauging station, after which the river incised back to near its initial elevation at this site. Within 15 weeks after breaching, over 1 million m 3 of suspended load is estimated to have passed the gauging station, consistent with estimates that >60% of the reservoir's sediment had eroded. This dam removal highlights the influence of interactions among reservoir erosion processes, sediment composition, and style of decommissioning on rate of reservoir erosion and consequent downstream behavior of released sediment.

Washington↗