Correction to “Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration”
No abstract is available for this article.
SEARCH · Geology Reports
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.
No abstract is available for this article.
Anthropogenic and climatic forces have modified the geomorphology of tidal wetlands over a range of timescales. Changes in land use, sediment supply, river flow, storminess, and sea level alter the layout of tidal channels, intertidal flats, and marsh plains; these elements define wetland complexes. Diagnostically, measurements of net sediment fluxes through tidal channels are high-temporal resolution, spatially integrated quantities that indicate (1) whether a complex is stable over seasonal timescales and (2) what mechanisms are leading to that state. We estimated sediment fluxes through tidal channels draining wetland complexes on the Blackwater and Transquaking Rivers, Maryland, USA. While the Blackwater complex has experienced decades of degradation and been largely converted to open water, the Transquaking complex has persisted as an expansive, vegetated marsh. The measured net export at the Blackwater complex (1.0 kg/s or 0.56 kg/m 2 /yr over the landward marsh area) was caused by northwesterly winds, which exported water and sediment on the subtidal timescale; tidally forced net fluxes were weak and precluded landward transport of suspended sediment from potential seaward sources. Though wind forcing also exported sediment at the Transquaking complex, strong tidal forcing and proximity to a turbidity maximum led to an import of sediment (0.031 kg/s or 0.70 kg/m 2 /yr). This resulted in a spatially averaged accretion of 3.9 mm/yr, equaling the regional relative sea level rise. Our results suggest that in areas where seaward sediment supply is dominant, seaward wetlands may be more capable of withstanding sea level rise over the short term than landward wetlands. We propose a conceptual model to determine a complex's tendency toward stability or instability based on sediment source, wetland channel location, and transport mechanisms. Wetlands with a reliable portfolio of sources and transport mechanisms appear better suited to offset natural and anthropogenic loss.
Glen Canyon Dam has caused a fundamental change in the distribution of fine sediment storage in the 99-km reach of the Colorado River in Marble Canyon, Grand Canyon National Park, Arizona. The two major storage sites for fine sediment (i.e., sand and finer material) in this canyon river are lateral recirculation eddies and the main-channel bed. We use a combination of methods, including direct measurement of sediment storage change, measurements of sediment flux, and comparison of the grain size of sediment found in different storage sites relative to the supply and that in transport, in order to evaluate the change in both the volume and location of sediment storage. The analysis shows that the bed of the main channel was an important storage environment for fine sediment in the predam era. In years of large seasonal accumulation, approximately 50% of the fine sediment supplied to the reach from upstream sources was stored on the main-channel bed. In contrast, sediment budgets constructed for two short-duration, high experimental releases from Glen Canyon Dam indicate that approximately 90% of the sediment discharge from the reach during each release was derived from eddy storage, rather than from sandy deposits on the main-channel bed. These results indicate that the majority of the fine sediment in Marble Canyon is now stored in eddies, even though they occupy a small percentage (???17%) of the total river area. Because of a 95% reduction in the supply of fine sediment to Marble Canyon, future high releases without significant input of tributary sediment will potentially erode sediment from long-term eddy storage, resulting in continued degradation in Marble Canyon. Copyright 2006 by the American Geophysical Union.
[1] Evolution of a tectonically active coast is driven by geomorphically destructive energy supplied by ocean waves. Wave energy is episodic and concentrated; sea cliffs are battered by the geomorphic wrecking ball every 4-25 s. We measure the response of sea cliffs to wave assault by sensing the ground motion using near-coastal seismometers. Sea cliffs respond to waves in two distinct styles. High-frequency motion (20 Hz) reflects the natural frequency of the sea cliff as it rings in response to direct wave impact. Low-frequency motion in the 0.1-0.05 Hz (10-20 s) band consistently agrees with the dominant nearshore wave period. Integrating microseismic velocities suggests 50 ??m and 10 ??m displacements in horizontal and vertical directions, respectively. Displacement ellipsoids exhibit simultaneous downward and seaward sea cliff motion with each wave. Video footage corroborates the downward sea cliff flex in response to the imposed water load on the wave cut platform. Gradients in displacement amplitudes documented using multiple seismometers suggest longitudinal and shear strain of the flexing sea cliff on the order of 0.5-4 ?? strains during each wave loading cycle. As this sea cliff flexure occurs approximately 3 million times annually, it has the potential to fatigue the rock through cyclical loading. Local sea cliff retreat rates of 10 cm/yr imply that a given parcel of rock is flexed through roughly 109 cycles of increasing amplitude before exposure to direct wave attack at the cliff face. Copyright 2005 by the American Geophysical Union.
[1] Field measurements of debris-bed friction on a smooth rock tablet at the bed of Engabreen, a hard-bedded, temperate glacier in northern Norway, indicated that basal ice containing 10% debris by volume exerted local shear traction of up to 500 kPa. The corresponding bulk friction coefficient between the dirty basal ice and the tablet was between 0.05 and 0.08. A model of friction in which nonrotating spherical rock particles are held in frictional contact with the bed by bed-normal ice flow can account for these measurements if the power law exponent for ice flowing past large clasts is 1. A small exponent (n < 2) is likely because stresses in ice are small and flow is transient. Numerical calculations of the bed-normal drag force on a sphere in contact with a flat bed using n = 1 show that this force can reach values several hundred times that on a sphere isolated from the bed, thus drastically increasing frictional resistance. Various estimates of basal friction are obtained from this model. For example, the shear traction at the bed of a glacier sliding at 20 m a-1 with a geothermally induced melt rate of 0.006 m a-1 and an effective pressure of 300 kPa can exceed 100 kPa. Debris-bed friction can therefore be a major component of sliding resistance, contradicting the common assumption that debris-bed friction is negligible. Copyright 2005 by the American Geophysical Union.
The Colorado River in Grand Canyon has long been known as a "rapids-and-pools" river, with the rapids owing their existence primarily to tributary debris flows. The debris flows deposit subaerial debris fans that constrict the channel laterally and, when they enter the river, raise the bed elevation. The rapids are short-wavelength (???0.1 to ???1 km), small-amplitude (??????5 m) convexities in the river's longitudinal profile, arising from the shallow gradient in the upstream pool and the steep gradient through the rapid itself. Analysis of the entire longitudinal profile through Grand Canyon reveals two long-wavelength (???100 km), large-amplitude (15-30 m) river profile convexities: the eastern canyon convexity between river mile (RM) 30 and RM 80 and the western canyon convexity between RM 150 and RM 250. Convexities of intermediate scale are also identified in the longitudinal profile. These longer-wavelength, larger-amplitude convexities have strong spatial correlations with high rates of debris flow occurrence, high densities of Holocene debris fans, the largest debris fans along the river, and alluvial thicknesses of 10 m or more. River profile convexities are unstable and require an active and powerful geologic process to maintain them, in this case the abundant, frequent, and voluminous Holocene debris flow activity in Grand Canyon. At all wavelengths the most likely cause for these river profile convexities is Holocene aggradation of the riverbed beneath them, driven by the coarse particles of tributary debris flows. Large enough debris flows will slow river flow for kilometers upstream, causing it to drop much of its suspended load. Integrated over time and all of the tributary point source contributions, this process will build short-wavelength convexities into long-wavelength convexities. For most if not all of the Holocene the Colorado River has been dissipating most of its energy in the rapids and expending the remainder in transporting fine sediment through Grand Canyon, with little or no regional incision of bedrock.
Development of a practical technology for rapid quantification of bed load transport in large rivers would represent a revolutionary advance for sediment monitoring and the investigation of fluvial dynamics. Measurement of bed load motion with acoustic Doppler current profiles (ADCPs) has emerged as a promising approach for evaluating bed load transport. However, a better understanding of how ADCP data relate to conditions near the stream bed is necessary to make the method practical for quantitative applications. In this paper, we discuss the response of ADCP bed velocity measurements, defined as the near-bed sediment velocity detected by the instrument's bottom-tracking feature, to changing sediment-transporting conditions in the lower Missouri River. Bed velocity represents a weighted average of backscatter from moving bed load particles and spectral reflections from the immobile bed. The ratio of bed velocity to mean bed load particle velocity depends on the concentration of the particles moving in the bed load layer, the bed load layer thickness, and the backscatter strength from a unit area of moving particles relative to the echo strength from a unit area of unobstructed bed. A model based on existing bed load transport theory predicted measured bed velocities from hydraulic and grain size measurements with reasonable success. Bed velocities become more variable and increase more rapidly with shear stress when the transport stage, defined as the ratio of skin friction to the critical shear stress for particle entrainment, exceeds a threshold of about 17. This transition in bed velocity response appears to be associated with the appearance of longer, flatter bed forms at high transport stages.
We studied soils on high-purity limestones of Quaternary age on the western Atlantic Ocean islands of Barbados, the Florida Keys, and the Bahamas. Potential soil parent materials in this region, external to the carbonate substrate, include volcanic ash from the island of St. Vincent (near Barbados), volcanic ash from the islands of Dominica and St. Lucia (somewhat farther from Barbados), the fine-grained component of distal loess from the lower Mississippi River Valley, and wind-transported dust from Africa. These four parent materials can be differentiated using trace elements (Sc, Cr, Th, and Zr) and rare earth elements that have minimal mobility in the soil-forming environment. Barbados soils have compositions that indicate a complex derivation. Volcanic ash from the island of St. Vincent appears to have been the most important influence, but African dust is a significant contributor, and even Mississippi River valley loess may be a very minor contributor to Barbados soils. Soils on the Florida Keys and islands in the Bahamas appear to have developed mostly from African dust, but Mississippi River valley loess may be a significant contributor. Our results indicate that inputs of African dust are more important to the genesis of soils on islands in the western Atlantic Ocean than previously supposed. We hypothesize that African dust may also be a major contributor to soils on other islands of the Caribbean and to soils in northern South America, central America, Mexico, and the southeastern United States. Dust inputs to subtropical and tropical soils in this region increase both nutrient-holding capacity and nutrient status and thus may be critical in sustaining vegetation. Copyright 2007 by the American Geophysical Union.
We examine the impacts of a stand-clearing wildfire on the characteristics and magnitude of suspended sediment transport in ephemeral streams draining the burn area. We report the results of a monitoring program that includes 2 years of data prior to the Cerro Grande fire in New Mexico, and 3 years of postfire data. Suspended sediment concentration (SSC) increased by about 2 orders of magnitude following the fire, and the proportion of silt and clay increased from 50% to 80%. For a given flow event, SSC is highest at the flood bore and decreases monotonically with time, a pattern evident in every flood sampled both before and after the fire. We propose that the accumulation of flow and wash load at the flow front is an inherent characteristic of ephemeral stream flows, due to amplified momentum losses at the flood bore. We present a new model for computing suspended sediment transport in ephemeral streams (in the presence or absence of wildfire) by relating SSC to the time following the arrival of the flood bore, rather than to instantaneous discharge. Using this model and a rainfall history, we estimate that in the 3 years following the fire, floods transported in suspension a mass equivalent to about 3 mm of landscape lowering across the burn area, 20% of this following a single rainstorm. We test the model by computing fine sediment delivery to a small reservoir in an adjacent watershed, where we have a detailed record of postfire sedimentation based on repeat surveys. Systematic discrepancies between modeled and measured sedimentation rates in the reservoir suggest rapid reductions in suspended sediment delivery in the first several years after the fire.
We employ a numerical model to study the development of sorted bed forms under a variety of hydrodynamic and sedimentary conditions. Results indicate that increased variability in wave height decreases the growth rate of the features and can potentially give rise to complicated, a priori unpredictable, behavior. This happens because the system responds to a change in wave characteristics by attempting to self-organize into a patterned seabed of different geometry and spacing. The new wavelength might not have enough time to emerge before a new change in wave characteristics occurs, leading to less regular seabed configurations. The new seabed configuration is also highly dependent on the preexisting morphology, which further limits the possibility of predicting future behavior. For the same reasons, variability in the mean current magnitude and direction slows down the growth of features and causes patterns to develop that differ from classical sorted bed forms. Spatial variability in grain size distribution and different types of net sediment aggradation/degradation can also result in the development of sorted bed forms characterized by a less regular shape. Numerical simulations qualitatively agree with observed geometry (spacing and height) of sorted bed forms. Also in agreement with observations is that at shallower depths, sorted bed forms are more likely to be affected by changes in the forcing conditions, which might also explain why, in shallow waters, sorted bed forms are described as ephemeral features. Finally, simulations indicate that the different sorted bed form shapes and patterns observed in the field might not necessarily be related to diverse physical mechanisms. Instead, variations in sorted bed form characteristics may result from variations in local hydrodynamic and/or sedimentary conditions.
The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 ?? 5.7 dB km-1) is somewhat lower than the value derived from radar profiles (25.3 ?? 1.1 dB km-1). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 ?? 2.2 dB km-1. This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet. Copyright 2007 by the American Geophysical Union.
Wet-based portions of ice sheets may move primarily by shearing their till beds, resting in high sediment fluxes and the development of subglacial landforms. This model of glacier movement, which requires high bed shear strains, can be tested using till microstructural characteristics that evolve during till deformation. Here we examine the development of magnetic fabric using a ring shear device to defom two Wisconsin-age basal tills to shear strains as high as 70. Hysteresis experiments and the dependence of magnetic susceptibility of these tills on temperature demonstrate that anisotropy of magnetic susceptibility (AMS) develops during shear due to the rotation of primarily magnetite particles that are silt sized or smaller. At moderate shear strains (???6-25), principal axes of maximum magnetic susceptibility develop a strong fabric (S1 eignevalues of 0.83-0.96), without further strengthening at higher strains, During deformation, directions of maximum susceptibility cluster strongly in the direction of shear and plunge 'up-glacier,' consistent with the behavior of pebbles and sand particles studied in earlier experiments. In contrast, the magnitude of AMS does not vary systematically with strain and is small relative to its variability among samples; this is because most magnetite grains are contained as inclusions in larger particles and hence do not align during shear. Although processes other than pervasive bed deformation may result in strong flow parallel fabrics, AMS fabrics provide a rapid and objective means of identifying basal tills that have not been sheared sufficiently to be compatible with the bed deformation model. Copyright 2008 by the American Geophysical Union.
The cross-sectional shape of a natural river channel controls the capacity of the system to carry water off a landscape, to convey sediment derived from hillslopes, and to erode its bed and banks. Numerical models that describe the response of a landscape to changes in climate or tectonics therefore require formulations that can accommodate evolution of channel cross-sectional geometry. However, fully two-dimensional (2-D) flow models are too computationally expensive to implement in large-scale landscape evolution models, while available simple empirical relationships between width and discharge do not adequately capture the dynamics of channel adjustment. We have developed a simplified 2-D numerical model of channel evolution in a cohesive, detachment-limited substrate subject to steady, unidirectional flow. Erosion is assumed to be proportional to boundary shear stress, which is calculated using an approximation of the flow field in which log-velocity profiles are assumed to apply along vectors that are perpendicular to the local channel bed. Model predictions of the velocity structure, peak boundary shear stress, and equilibrium channel shape compare well with predictions of a more sophisticated but more computationally demanding ray-isovel model. For example, the mean velocities computed by the two models are consistent to within ???3%, and the predicted peak shear stress is consistent to within ???7%. Furthermore, the shear stress distributions predicted by our model compare favorably with available laboratory measurements for prescribed channel shapes. A modification to our simplified code in which the flow includes a high-velocity core allows the model to be extended to estimate shear stress distributions in channels with large width-to-depth ratios. Our model is efficient enough to incorporate into large-scale landscape evolution codes and can be used to examine how channels adjust both cross-sectional shape and slope in response to tectonic and climatic forcing. Copyright 2008 by the American Geophysical Union.
In this paper we present a detailed comparison between measured features of the Golden Gate sand wave field and the results of a nonlinear sand wave model. Because the Golden Gate sand waves exhibit large variation in their characteristics and in their environmental physics, this area gives us the opportunity to study sand wave variation between locations, within one well-measured, large area. The nonlinear model used in this paper is presently the only tool that provides information on the nonlinear evolution of large-amplitude sand waves. The model is used to increase our understanding of the coupling between the variability in environmental conditions and the sand wave characteristics. Results show that the model is able to describe the variation in the Golden Gate sand waves well when both the local oscillating tidal current and the residual current are taken into account. Current and water depth seem to be the most important factors influencing sand wave characteristics. The simulation results give further confidence in the underlying model hypothesis and assumptions. Copyright 2009 by the American Geophysical Union.
Hindcasting decadal-timescale bathymetric change in estuaries is prone to error due to limited data for initial conditions, boundary forcing, and calibration; computational limitations further hinder efforts. We developed and calibrated a tidal-timescale model to bathymetric change in Suisun Bay, California, over the 1867–1887 period. A general, multiple-timescale calibration ensured robustness over all timescales; two input reduction methods, the morphological hydrograph and the morphological acceleration factor, were applied at the decadal timescale. The model was calibrated to net bathymetric change in the entire basin; average error for bathymetric change over individual depth ranges was 37%. On a model cell-by-cell basis, performance for spatial amplitude correlation was poor over the majority of the domain, though spatial phase correlation was better, with 61% of the domain correctly indicated as erosional or depositional. Poor agreement was likely caused by the specification of initial bed composition, which was unknown during the 1867–1887 period. Cross-sectional bathymetric change between channels and flats, driven primarily by wind wave resuspension, was modeled with higher skill than longitudinal change, which is driven in part by gravitational circulation. The accelerated response of depth may have prevented gravitational circulation from being represented properly. As performance criteria became more stringent in a spatial sense, the error of the model increased. While these methods are useful for estimating basin-scale sedimentation changes, they may not be suitable for predicting specific locations of erosion or deposition. They do, however, provide a foundation for realistic estuarine geomorphic modeling applications.
In canyon rivers, debris fan constrictions create rapids and downstream pools characterized by secondary flow structures that are closely linked to channel morphology. In this paper we describe detailed measurements of the three‐dimensional flow structure and sandbar dynamics of two pools along the Colorado River in the Grand Canyon during a controlled flood release from Glen Canyon Dam. Results indicate that the pools are characterized by large lateral recirculation zones (eddies) resulting from flow separation downstream from the channel constrictions, as well as helical flow structures in the main channel and eddy. The lateral recirculation zones are low‐velocity areas conducive to fine sediment deposition, particularly in the vicinity of the separation and reattachment points and are thus the dominant flow structures controlling sandbar dynamics. The helical flow structures also affect morphology but appear secondary in importance to the lateral eddies. During the controlled flood, sandbars in the separation and reattachment zones at both sites tended to build gradually during the rising limb and peak flow. Deposition in shallow water on the sandbars was accompanied by erosion in deeper water along the sandbar slope at the interface with the main channel. Erosion occurred via rapid mass failures as well as by gradual boundary shear stress driven processes. The flow structures and morphologic links at our study sites are similar to those identified in other river environments, in particular sharply curved meanders and channel confluences where the coexistence of lateral recirculation and helical flows has been documented.
Since the installation of the Alaska Regional Seismic Network in the 1970s, data analysts have noted nontectonic seismic events thought to be related to glacier dynamics. While loose associations with the glaciers of the St. Elias Mountains have been made, no detailed study of the source locations has been undertaken. We performed a two-step investigation surrounding these events, beginning with manual locations that guided an automated detection and event sifting routine. Results from the manual investigation highlight characteristics of the seismic waveforms including single-peaked (narrowband) spectra, emergent onsets, lack of distinct phase arrivals, and a predominant cluster of locations near the calving termini of several neighboring tidewater glaciers. Through these locations, comparison with previous work, analyses of waveform characteristics, frequency-magnitude statistics and temporal patterns in seismicity, we suggest calving as a source for the seismicity. Statistical properties and time series analysis of the event catalog suggest a scale-invariant process that has no single or simple forcing. These results support the idea that calving is often a response to short-lived or localized stress perturbations. Our results demonstrate the utility of passive seismic instrumentation to monitor relative changes in the rate and magnitude of iceberg calving at tidewater glaciers that may be volatile or susceptible to ensuing rapid retreat, especially when existing seismic infrastructure can be used.
Recent work has demonstrated that the topological properties of real river networks deviate significantly from predictions of Shreve's random model. At the same time the property of mean self-similarity postulated by Tokunaga's model is well supported by data. Recently, a new class of network model called random self-similar networks (RSN) that combines self-similarity and randomness has been introduced to replicate important topological features observed in real river networks. We investigate if the hypothesis of statistical self-similarity in the RSN model is supported by data on a set of 30 basins located across the continental United States that encompass a wide range of hydroclimatic variability. We demonstrate that the generators of the RSN model obey a geometric distribution, and self-similarity holds in a statistical sense in 26 of these 30 basins. The parameters describing the distribution of interior and exterior generators are tested to be statistically different and the difference is shown to produce the well-known Hack's law. The inter-basin variability of RSN parameters is found to be statistically significant. We also test generator dependence on two climatic indices, mean annual precipitation and radiative index of dryness. Some indication of climatic influence on the generators is detected, but this influence is not statistically significant with the sample size available. Finally, two key applications of the RSN model to hydrology and geomorphology are briefly discussed.