SEARCH · USGS Science

Results for “Journal of Geophysical Research: Earth Surface”

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 145 records · Page 8Linked to original sources

USGS2006

Reconstructing depositional processes and history from reservoir stratigraphy: Englebright Lake, Yuba River, northern California

Reservoirs provide the opportunity to link watershed history with its stratigraphic record. We analyze sediment cores from a northern California reservoir in the context of hydrologic history, watershed management, and depositional processes. Observations of recent depositional patterns, sediment‐transport calculations, and 137 Cs geochronology support a conceptual model in which the reservoir delta progrades during floods of short duration (days) and is modified during prolonged (weeks to months) drawdowns that rework topset beds and transport sand from topsets to foresets. Sediment coarser than 0.25–0.5 mm deposits in foresets and topsets, and finer material falls out of suspension as bottomset beds. Simple hydraulic calculations indicate that fine sand (0.063–0.5 mm) is transported into the distal bottomset area only during floods. The overall stratigraphy suggests that two phases of delta building occurred in the reservoir. The first, from dam construction in 1940 to 1970, was heavily influenced by annual, prolonged >20 m drawdowns of the water level. The second, built on top of the first, reflects sedimentation from 1970 to 2002 when the influence of drawdowns was less. Sedimentation rates in the central part of the reservoir have declined ∼25% since 1970, likely reflecting a combination of fewer large floods, changes in watershed management, and winnowing of stored hydraulic mining sediment.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Suspended sediment transport in an ephemeral stream following wildfire

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2006

Form drag in rivers due to small-scale natural topographic features: 1. Regular sequences

Small-scale topographic features are commonly found on the boundaries of natural rivers, streams, and floodplains. A simple method for determining the form drag on these features is presented, and the results of this model are compared to laboratory measurements. The roughness elements are modeled as Gaussian-shaped features defined in terms of three parameters: a protrusion height, H; a streamwise length scale, ??; and a spacing between crests, ??. This shape is shown to be a good approximation to a wide variety of natural topographic bank features. The form drag on an individual roughness element embedded in a series of identical elements is determined using the drag coefficient of the individual element and a reference velocity that includes the effects of roughness elements further upstream. In addition to calculating the drag on each element, the model determines the spatially averaged total stress, skin friction stress, and roughness height of the boundary. The effects of bank roughness on patterns of velocity and boundary shear stress are determined by combining the form drag model with a channel flow model. The combined model shows that drag on small-scale topographic features substantially alters the near-bank flow field. These methods can be used to improve predictions of flow resistance in rivers and to form the basis for fully predictive (no empirically adjusted parameters) channel flow models. They also provide a foundation for calculating the near-bank boundary shear stress fields necessary for determining rates of sediment transport and lateral erosion.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Intertidal sand body migration along a megatidal coast, Kachemak Bay, Alaska

[1] Using a digital video-based Argus Beach Monitoring System (ABMS) on the north shore of Kachemak Bay in south central Alaska, we document the timing and magnitude of alongshore migration of intertidal sand bed forms over a cobble substrate during a 22-month observation period. Two separate sediment packages (sand bodies) of 1–2 m amplitude and ∼200 m wavelength, consisting of well-sorted sand, were observed to travel along shore at annually averaged rates of 278 m/yr (0.76 m/d) and 250 m/yr (0.68 m/d), respectively. Strong seasonality in migration rates was shown by the contrast of rapid winter and slow summer transport. Though set in a megatidal environment, data indicate that sand body migration is driven by eastward propagating wind waves as opposed to net westward directed tidal currents. Greatest weekly averaged rates of movement, exceeding 6 m/d, coincided with wave heights exceeding 2 m suggesting a correlation of wave height and sand body migration. Because Kachemak Bay is partially enclosed, waves responsible for sediment entrainment and transport are locally generated by winds that blow across lower Cook Inlet from the southwest, the direction of greatest fetch. Our estimates of sand body migration translate to a littoral transport rate between 4,400–6,300 m 3 /yr. Assuming an enclosed littoral cell, minimal riverine sediment contributions, and a sea cliff sedimentary fraction of 0.05, we estimate long-term local sea cliff retreat rates of 9–14 cm/yr. Applying a numerical model of wave energy dissipation to the temporally variable beach morphology suggests that sand bodies are responsible for enhancing wave energy dissipation by ∼13% offering protection from sea cliff retreat.

Alaska↗

USGS2007

Forces on stationary particles in near-bed turbulent flows

In natural flows, bed sediment particles are entrained and moved by the fluctuating forces, such as lift and drag, exerted by the overlying flow on the particles. To develop a better understanding of these forces and the relation of the forces to the local flow, the downstream and vertical components of force on near-bed fixed particles and of fluid velocity above or in front of them were measured synchronously at turbulence-resolving frequencies (200 or 500 Hz) in a laboratory flume. Measurements were made for a spherical test particle fixed at various heights above a smooth bed, above a smooth bed downstream of a downstream-facing step, and in a gravel bed of similarly sized particles as well as for a cubical test particle and 7 natural particles above a smooth bed. Horizontal force was well correlated with downstream velocity and not correlated with vertical velocity or vertical momentum flux. The standard drag formula worked well to predict the horizontal force, but the required value of the drag coefficient was significantly higher than generally used to model bed load motion. For the spheres, cubes, and natural particles, average drag coefficients were found to be 0.76, 1.36, and 0.91, respectively. For comparison, the drag coefficient for a sphere settling in still water at similar particle Reynolds numbers is only about 0.4. The variability of the horizontal force relative to its mean was strongly increased by the presence of the step and the gravel bed. Peak deviations were about 30% of the mean force for the sphere over the smooth bed, about twice the mean with the step, and 4 times it for the sphere protruding roughly half its diameter above the gravel bed. Vertical force correlated poorly with downstream velocity, vertical velocity, and vertical momentum flux whether measured over or ahead of the test particle. Typical formulas for shear-induced lift based on Bernoulli's principle poorly predict the vertical forces on near-bed particles. The measurements suggest that particle-scale pressure variations associated with turbulence are significant in the particle momentum balance. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Methods to assess natural and anthropogenic thaw lake drainage on the western Arctic coastal plain of northern Alaska

Thousands of lakes are found on the Arctic Coastal Plain of northern Alaska and northwestern Canada. Developed atop continuous permafrost, these thaw lakes and associated drained thaw lake basins are the dominant landscape elements and together cover 46% of the 34,570 km 2 western Arctic Coastal Plain (WACP). Lakes drain by a variety of episodic processes, including coastal erosion, stream meandering, and headward erosion, bank overtopping, and lake coalescence. Comparison of Landsat multispectral scanner (MSS) imagery from the mid-1970s to Landsat 7 enhanced thematic mapper (ETM+) imagery from around 2000 shows that 50 lakes completely or partially drained over the approximately 25 year period, indicating landscape stability. The lake-specific drainage mechanism can be inferred in some cases and is partially dependant on geographic settings conducive to active erosion such as riparian and coastal zones. In many cases, however, the cause of drainage is unknown. The availability of high-resolution aerial photographs for the Barrow Peninsula extends the record back to circa 1950; mapping spatial time series illustrates the dynamic nature of lake expansion, coalescence, and drainage. Analysis of these historical images suggests that humans have intentionally or inadvertently triggered lake drainage near the village of Barrow. Efforts to understand landscape processes and identify events have been enhanced by interviewing Iñupiaq elders and others practicing traditional subsistence lifestyles. They can often identify the year and process by which individual lakes drained, thereby providing greater dating precision and accuracy in assessing the causal mechanism. Indigenous knowledge has provided insights into events, landforms, and processes not previously identified or considered.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Suspended-sediment rating curve response to urbanization and wildfire, Santa Ana River, California

[1] River suspended-sediment concentrations provide insights to the erosion and transport of materials from a landscape, and changes in concentrations with time may result from landscape processes or human disturbance. Here we show that suspended-sediment concentrations in the Santa Ana River, California, decreased 20-fold with respect to discharge during a 34-year period (1968−2001). These decreases cannot be attributed to changes in sampling technique or timing, nor to event or seasonal hysteresis. Annual peak and total discharge, however, reveal sixfold increases over the 34-year record, which largely explain the decreases in sediment concentration by a nonlinear dilution process. The hydrological changes were related to the widespread urbanization of the watershed, which resulted in increases in storm water discharge without detectable alteration of sediment discharge, thus reducing suspended-sediment concentrations. Periodic upland wildfire significantly increased water discharge, sediment discharge, and suspended-sediment concentrations and thus further altered the rating curve with time. Our results suggest that previous inventories of southern California sediment flux, which assume time-constant rating curves and extend these curves beyond the sampling history, may have substantially overestimated loads during the most recent decades.

California↗

USGS2007

Sorted bed forms as self-organized patterns: 2. complex forcing scenarios

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2007

Modeling englacial radar attenuation at Siple Dome, West Antarctica, using ice chemistry and temperature data

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2008

Magnetic fabric of sheared till: A strain indicator for evaluating the bed deformation model of glacier flow

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2008

Modeling the evolution of channel shape: Balancing computational efficiency with hydraulic fidelity

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2009

Simulating and understanding sand wave variation: A case study of the Golden Gate sand waves

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2009

Hindcasting of decadal‐timescale estuarine bathymetric change with a tidal‐timescale model

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.

California↗

USGS2010

Complexities in barrier island response to sea level rise: Insights from numerical model experiments, North Carolina Outer Banks

Using a morphological-behavior model to conduct sensitivity experiments, we investigate the sea level rise response of a complex coastal environment to changes in a variety of factors. Experiments reveal that substrate composition, followed in rank order by substrate slope, sea level rise rate, and sediment supply rate, are the most important factors in determining barrier island response to sea level rise. We find that geomorphic threshold crossing, defined as a change in state (e.g., from landward migrating to drowning) that is irreversible over decadal to millennial time scales, is most likely to occur in muddy coastal systems where the combination of substrate composition, depth-dependent limitations on shoreface response rates, and substrate erodibility may prevent sand from being liberated rapidly enough, or in sufficient quantity, to maintain a subaerial barrier. Analyses indicate that factors affecting sediment availability such as low substrate sand proportions and high sediment loss rates cause a barrier to migrate landward along a trajectory having a lower slope than average barrier island slope, thereby defining an “effective” barrier island slope. Other factors being equal, such barriers will tend to be smaller and associated with a more deeply incised shoreface, thereby requiring less migration per sea level rise increment to liberate sufficient sand to maintain subaerial exposure than larger, less incised barriers. As a result, the evolution of larger/less incised barriers is more likely to be limited by shoreface erosion rates or substrate erodibility making them more prone to disintegration related to increasing sea level rise rates than smaller/more incised barriers. Thus, the small/deeply incised North Carolina barriers are likely to persist in the near term (although their long-term fate is less certain because of the low substrate slopes that will soon be encountered). In aggregate, results point to the importance of system history (e.g., previous slopes, sediment budgets, etc.) in determining migration trajectories and therefore how a barrier island will respond to sea level rise. Although simple analytical calculations may predict barrier response in simplified coastal environments (e.g., constant slope, constant sea level rise rate, etc.), our model experiments demonstrate that morphological-behavior modeling is necessary to provide critical insights regarding changes that may occur in environments having complex geometries, especially when multiple parameters change simultaneously.

North Carolina↗

USGS2010

Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration

Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.

Journal of Geophysical Research F: Earth Surface↗

USGS2011

Flow structures and sandbar dynamics in a canyon river during a controlled flood, Colorado River, Arizona

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.

Journal of Geophysical Research F: Earth Surface↗

USGS2010

A universal approximation to grain size from images of non-cohesive sediment

The two-dimensional spectral decomposition of an image of sediment provides a direct statistical estimate, grid-by-number style, of the mean of all intermediate axes of all single particles within the image. We develop and test this new method which, unlike existing techniques, requires neither image processing algorithms for detection and measurement of individual grains, nor calibration. The only information required of the operator is the spatial resolution of the image. The method is tested with images of bed sediment from nine different sedimentary environments (five beaches, three rivers, and one continental shelf), across the range 0.1 mm to 150 mm, taken in air and underwater. Each population was photographed using a different camera and lighting conditions. We term it a “universal approximation” because it has produced accurate estimates for all populations we have tested it with, without calibration. We use three approaches (theory, computational experiments, and physical experiments) to both understand and explore the sensitivities and limits of this new method. Based on 443 samples, the root-mean-squared (RMS) error between size estimates from the new method and known mean grain size (obtained from point counts on the image) was found to be ±≈16%, with a 95% probability of estimates within ±31% of the true mean grain size (measured in a linear scale). The RMS error reduces to ≈11%, with a 95% probability of estimates within ±20% of the true mean grain size if point counts from a few images are used to correct bias for a specific population of sediment images. It thus appears it is transferable between sedimentary populations with different grain size, but factors such as particle shape and packing may introduce bias which may need to be calibrated for. For the first time, an attempt has been made to mathematically relate the spatial distribution of pixel intensity within the image of sediment to the grain size.

Journal of Geophysical Research F: Earth Surface↗

USGS2010

Iceberg calving as a primary source of regional‐scale glacier‐generated seismicity in the St. Elias Mountains, Alaska

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.

Alaska↗