Search USGS⌕ Search

SEARCH · Search USGS

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.

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

At least 91 records · Page 5Linked to original sources
USGS2009

Hydrometeor-enhanced tephra sedimentation: Constraints from the 18 May 1980 eruption of Mount St. Helens

Uncertainty remains on the origin of distal mass deposition maxima observed in many recent tephra fall deposits. In this study the link between ash aggregation and the formation of distal mass deposition maxima is investigated through reanalysis of tephra fallout from the Mount St. Helens 18 May 1980 (MSH80) eruption. In addition, we collate all the data needed to model distal ash sedimentation from the MSH80 eruption cloud. Four particle size subpopulations were present in distal fallout with modes at 2.2 Φ, 4.2 Φ, 5.9 Φ, and 8.3 Φ. Settling rates of the coarsest subpopulation closely matched predicted single‐particle terminal fall velocities. Sedimentation of particles <100 μ m was greatly enhanced, predominantly through aggregation of a particle subpopulation with modal diameter 5.9 ± 0.2 Φ (19 ± 3 μ m). Mammatus on the MSH80 cloud provided a mechanism to transport very fine ash particles, with predicted atmospheric lifetimes of days to weeks, from the upper troposphere to the surface in a matter of hours. In this mechanism, ash particles initiate ice hydrometeor formation high in the troposphere. Subsequently, the volcanic cloud rapidly subsides as mammatus develop from increased particle loading and cloud base sublimation. Rapid fallout occurs as the cloud passes through the melting level in a process analogous to snowflake aggregation. Aggregates sediment en masse and form the distal mass deposition maxima observed in many recent volcanic ash fall deposits. This work provides a data resource that will facilitate tephra sedimentation modeling and allow model intercomparisons.

Washington↗
USGS2007

Glacier changes in southeast Alaska and northwest British Columbia and contribution to sea level rise

The digital elevation model (DEM) from the 2000 Shuttle Radar Topography Mission (SRTM) was differenced from a composite DEM based on air photos dating from 1948 to 1987 to detennine glacier volume changes in southeast Alaska and adjoining Canada. SRTM accuracy was assessed at ??5 in through comparison with airborne laser altimetry and control locations measured with GPS. Glacier surface elevations lowered over 95% of the 14,580 km2 glacier-covered area analyzed, with some glaciers thinning as much as 640 in. A combination of factors have contributed to this wastage, including calving retreats of tidewater and lacustrine glaciers and climate change. Many glaciers in this region are particularly sensitive to climate change, as they have large areas at low elevations. However, several tidewater glaciers that had historically undergone calving retreats are now expanding and appear to be in the advancing stage of the tidewater glacier cycle. The net average rate of ice loss is estimated at 16.7 ?? 4.4 km3/yr, equivalent to a global sea level rise contribution of 0.04 ?? 0.01 mm/yr. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Distributed energy balance modeling of South Cascade Glacier, Washington and assessment of model uncertainty

We have developed a physically based, distributed surface energy balance model to simulate glacier mass balance under meteorological and climatological forcing. Here we apply the model to estimate summer ablation on South Cascade Glacier, Washington, for the 2004 and 2005 mass balance seasons. To arrive at optimal mass balance simulations, we investigate and quantify model uncertainty associated with selecting from a range of physical parameter values that are not commonly measured in glaciological mass balance field studies. We optimize the performance of the model by varying values for atmospheric transmissivity, the albedo of surrounding topography, precipitation-elevation lapse rate, surface roughness for turbulent exchange of momentum, and snow albedo aging coefficient. Of these the snow aging parameter and precipitation lapse rates have the greatest influence on the modeled ablation. We examined model sensitivity to varying parameters by performing an additional 103 realizations with parameters randomly chosen over a ??5% range centered about the optimum values. The best fit suite of model parameters yielded a net balance of -1.69??0.38 m water equivalent (WE) for the 2004 water year and -2.10??0.30 m WE up to 11 September 2005. The 2004 result is within 3% of the measured value. These simulations account for 91% and 93% of the variance in measured ablation for the respective years. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2008

Compositional trends in aeolian dust along a transect across the southwestern United States

Aeolian dust strongly influences ecology and landscape geochemistry over large areas that span several desert ecosystems of the southwestern United States. This study evaluates spatial and temporal variations and trends of the physical and chemical properties of dust in the southwestern United States by examining dust deposited in natural depressions on high isolated surfaces along a transect from the Mojave Desert to the central Colorado Plateau. Aeolian dust is recognized in these depressions on the basis of textural, chemical, isotopic, and mineralogical characteristics and comparisons of those characteristics to the underlying bedrock units. Spatial and temporal trends suggest that although local dust sources are important to the accumulated material in these depressions, Mojave Desert dust sources may also contribute. Depth trends in the depressions suggest that Mojave sources may have contributed more dust to the Colorado Plateau recently than in the past. These interpretations point to the important roles of far-traveled aeolian dust for landscape geochemistry and imply future changes to soil geochemistry under changing conditions in far-distant dust source areas. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2016

A synthesis of the basal thermal state of the Greenland Ice Sheet

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Nitrogen and carbon flow from rock to water: Regulation through soil biogeochemical processes, Mokelumne River watershed, California, and Grand Valley, Colorado

Soil denitrification is an ecologically important nitrogen removal mechanism that releases to the atmosphere the greenhouse gas N2O, an intermediate product from the reduction of NO3- to N 2. In this study we evaluate the relationship between soil carbon and denitrification potential in watersheds with bedrock acting as a nonpoint source of nitrogen, testing the hypothesis that nitrate leaching to stream water is in part regulated by denitrification. Two sites, one in a Mediterranean climate and the other in an arid climate, were investigated to understand the interplay between carbon and denitrification potential. Both sites included carbonaceous bedrock with relatively high nitrogen concentrations (> 1,000 mg N kg-1) and had low background nitrogen concentrations in surface and groundwater. There was a net accumulation of carbon and nitrogen in soil relative to the corresponding bedrock, with the exception of carbonaceous shale from the arid site. There the concentration of carbon in the soil (15,620 mg C kg-1) was less than the shale parent (22,460 mg C kg-1), consistent with the bedrock being a source of soil carbon. Rates of denitrification potential (0.5-83 ??g N kg-1 hr-1) derived from laboratory incubations appeared to be related to the ratio of dissolved organic carbon and nitrate extracted from soils. These data indicate that microbial processes such as denitrification can help maintain background nitrogen concentrations to tens of ??M N in relatively undisturbed ecosystems when nitrogen inputs from weathering bedrock are accompanied by sufficient organic carbon concentrations to promote microbial nitrogen transformations.

Journal of Geophysical Research F: Earth Surface↗
USGS2005

Evolving force balance at Columbia Glacier, Alaska, during its rapid retreat

Changes in driving and resistive stresses play an essential role in governing the buoyancy forces that are important controls on the speed and irreversibility of tidewater glacier retreats. We describe changes in geometry, velocity, and strain rate and present a top-down force balance analysis performed over the lower reach of Columbia Glacier. Our analysis uses new measurements and estimates of basal topography and photogrammetric surface velocity measurements made between 1977 and 2001, while assuming depth-independent strain. Sensitivity tests show that the method is robust and insensitive to small changes in the calculation parameters. Spatial distributions of ice speed show little correspondence with driving stress. Instead, spatial patterns of ice speed exhibit a nonlinear correspondence with basal drag. Primary resistance to flow comes from basal drag, but lateral drag becomes increasingly more important throughout the retreat, which may account for observed increases in speed. Maximum basal drag is always located in a prominent constriction located ~12 km upstream from the preretreat terminus. Once the terminus retreated into deep water off the terminal moraine marking the modern maximum extent, the upstream location of this maximum basal drag helped to promote thinning and decrease effective pressure in the lower region by limiting replenishing ice flow from upstream. An increase in both ice velocity and calving resulted, initiating what appears to be an irreversible retreat. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗
USGS2012

Displacement fields from point cloud data: Application of particle imaging velocimetry to landslide geodesy

Acquiring spatially continuous ground-surface displacement fields from Terrestrial Laser Scanners (TLS) will allow better understanding of the physical processes governing landslide motion at detailed spatial and temporal scales. Problems arise, however, when estimating continuous displacement fields from TLS point-clouds because reflecting points from sequential scans of moving ground are not defined uniquely, thus repeat TLS surveys typically do not track individual reflectors. Here, we implemented the cross-correlation-based Particle Image Velocimetry (PIV) method to derive a surface deformation field using TLS point-cloud data. We estimated associated errors using the shape of the cross-correlation function and tested the method's performance with synthetic displacements applied to a TLS point cloud. We applied the method to the toe of the episodically active Cleveland Corral Landslide in northern California using TLS data acquired in June 2005–January 2007 and January–May 2010. Estimated displacements ranged from decimeters to several meters and they agreed well with independent measurements at better than 9% root mean squared (RMS) error. For each of the time periods, the method provided a smooth, nearly continuous displacement field that coincides with independently mapped boundaries of the slide and permits further kinematic and mechanical inference. For the 2010 data set, for instance, the PIV-derived displacement field identified a diffuse zone of displacement that preceded by over a month the development of a new lateral shear zone. Additionally, the upslope and downslope displacement gradients delineated by the dense PIV field elucidated the non-rigid behavior of the slide.

California↗
USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

Instrumental record of debris flow initiation during natural rainfall: Implications for modeling slope stability

The middle of a hillslope hollow in the Oregon Coast Range failed and mobilized as a debris flow during heavy rainfall in November 1996. Automated pressure transducers recorded high spatial variability of pore water pressure within the area that mobilized as a debris flow, which initiated where local upward flow from bedrock developed into overlying colluvium. Postfailure observations of the bedrock surface exposed in the debris flow scar reveal a strong spatial correspondence between elevated piezometric response and water discharging from bedrock fractures. Measurements of apparent root cohesion on the basal ( C b ) and lateral ( C l ) scarp demonstrate substantial local variability, with areally weighted values of C b = 0.1 and C l = 4.6 kPa. Using measured soil properties and basal root strength, the widely used infinite slope model, employed assuming slope parallel groundwater flow, provides a poor prediction of hydrologic conditions at failure. In contrast, a model including lateral root strength (but neglecting lateral frictional strength) gave a predicted critical value of relative soil saturation that fell within the range defined by the arithmetic and geometric mean values at the time of failure. The 3‐D slope stability model CLARA‐W, used with locally observed pore water pressure, predicted small areas with lower factors of safety within the overall slide mass at sites consistent with field observations of where the failure initiated. This highly variable and localized nature of small areas of high pore pressure that can trigger slope failure means, however, that substantial uncertainty appears inevitable for estimating hydrologic conditions within incipient debris flows under natural conditions.

Oregon↗
USGS2010

The perfect debris flow? Aggregated results from 28 large-scale experiments

Aggregation of data collected in 28 controlled experiments reveals reproducible debris-flow behavior that provides a clear target for model tests. In each experiment ∼10 m 3 of unsorted, water-saturated sediment composed mostly of sand and gravel discharged from behind a gate, descended a steep, 95-m flume, and formed a deposit on a nearly horizontal runout surface. Experiment subsets were distinguished by differing basal boundary conditions (1 versus 16 mm roughness heights) and sediment mud contents (1 versus 7 percent dry weight). Sensor measurements of evolving flow thicknesses, basal normal stresses, and basal pore fluid pressures demonstrate that debris flows in all subsets developed dilated, coarse-grained, high-friction snouts, followed by bodies of nearly liquefied, finer-grained debris. Mud enhanced flow mobility by maintaining high pore pressures in flow bodies, and bed roughness reduced flow speeds but not distances of flow runout. Roughness had these effects because it promoted debris agitation and grain-size segregation, and thereby aided growth of lateral levees that channelized flow. Grain-size segregation also contributed to development of ubiquitous roll waves, which had diverse amplitudes exhibiting fractal number-size distributions. Despite the influence of these waves and other sources of dispersion, the aggregated data have well-defined patterns that help constrain individual terms in a depth-averaged debris-flow model. The patterns imply that local flow resistance evolved together with global flow dynamics, contradicting the hypothesis that any consistent rheology applied. We infer that new evolution equations, not new rheologies, are needed to explain how characteristic debris-flow behavior emerges from the interactions of debris constituents.

Journal of Geophysical Research F: Earth Surface↗
USGS2018

Enhancement of a parsimonious water balance model to simulate surface hydrology in a glacierized watershed

The U.S. Geological Survey monthly water balance model (MWBM) was enhanced with the capability to simulate glaciers in order to make it more suitable for simulating cold region hydrology. The new model, MWBMglacier, is demonstrated in the heavily glacierized and ecologically important Copper River watershed in Southcentral Alaska. Simulated water budget components compared well to satellite‐based observations and ground measurements of streamflow, evapotranspiration, snow extent, and total water storage, with differences ranging from 0.2% to 7% of the precipitation flux. Nash Sutcliffe efficiency for simulated and observed streamflow was greater than 0.8 for six of eight stream gages. Snow extent matched satellite‐based observations with Nash Sutcliffe efficiency values of greater than 0.89 in the four Copper River ecoregions represented. During the simulation period 1949 to 2009, glacier ice melt contributed 25% of total runoff, ranging from 12% to 45% in different tributaries, and glacierized area was reduced by 6%. Statistically significant ( p < 0.05) decreasing and increasing trends in annual glacier mass balance occurred during the multidecade cool and warm phases of the Pacific Decadal Oscillation, respectively, reinforcing the link between climate perturbations and glacier mass balance change. The simulations of glaciers and total runoff for a large, remote region of Alaska provide useful data to evaluate hydrologic, cryospheric, ecologic, and climatic trends. MWBM glacier is a valuable tool to understand when, and to what extent, streamflow may increase or decrease as glaciers respond to a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2009

Compositions of modern dust and surface sediments in the Desert Southwest, United States

Modern dusts across southwestern United States deserts are compositionally similar to dust-rich Av soil horizons (depths of 0-0.5 cm and 1-4 cm at 35 sites) for common crustal elements but distinctly different for some trace elements. Chemical compositions and magnetic properties of the soil samples are similar among sites relative to dust sources, geographic areas, and lithologic substrates. Exceptions are Li, U, and W, enriched in Owens Valley, California, and Mg and Sr, enriched in soils formed on calcareous fan gravel in southeast Nevada. The Av horizons are dominated by dust and reflect limited mixing with substrate sediments. Modern dust samples are also similar across the region, except that Owens Valley dusts are higher in Mg, Ba, and Li and dusts both there and at sites to the north on volcanic substrates are higher in Sb and W. Thus, dust and Av horizons consist of contributions from many different sources that are well mixed before deposition. Modern dusts contain significantly greater amounts of As, Cd, Cr, Cu, Ni, Pb, and Sb than do Av horizons, which record dust additions over hundreds to thousands of years. These results suggest that modern dust compositions are influenced by anthropogenic sources and emissions from Owens (dry) Lake after its artificial desiccation in 1926. Both modern dusts and Av horizons are enriched in As, Ba, Cu, Li, Sb, Th, U, and W relative to average crustal composition, which we interpret to indicate that the geologic sources of dust in the southwestern United States are geochemically distinctive.

Journal of Geophysical Research F: Earth Surface↗
USGS2013

Runoff-generated debris flows: observations and modeling of surge initiation, magnitude, and frequency

Runoff during intense rainstorms plays a major role in generating debris flows in many alpine areas and burned steeplands. Yet compared to debris flow initiation from shallow landslides, the mechanics by which runoff generates a debris flow are less understood. To better understand debris flow initiation by surface water runoff, we monitored flow stage and rainfall associated with debris flows in the headwaters of two small catchments: a bedrock-dominated alpine basin in central Colorado (0.06 km 2 ) and a recently burned area in southern California (0.01 km 2 ). We also obtained video footage of debris flow initiation and flow dynamics from three cameras at the Colorado site. Stage observations at both sites display distinct patterns in debris flow surge characteristics relative to rainfall intensity (I). We observe small, quasiperiodic surges at low I; large, quasiperiodic surges at intermediate I; and a single large surge followed by small-amplitude fluctuations about a more steady high flow at high I. Video observations of surge formation lead us to the hypothesis that these flow patterns are controlled by upstream variations in channel slope, in which low-gradient sections act as “sediment capacitors,” temporarily storing incoming bed load transported by water flow and periodically releasing the accumulated sediment as a debris flow surge. To explore this hypothesis, we develop a simple one-dimensional morphodynamic model of a sediment capacitor that consists of a system of coupled equations for water flow, bed load transport, slope stability, and mass flow. This model reproduces the essential patterns in surge magnitude and frequency with rainfall intensity observed at the two field sites and provides a new framework for predicting the runoff threshold for debris flow initiation in a burned or alpine setting.

California;Colorado↗
USGS2012

Volcanism, isostatic residual gravity and regional tectonic setting of the Cascade volcanic province

A technique to locate automatically boundaries between crustal blocks of disparate densities was applied to upward continued isostatic residual gravity data. The boundary analysis delineates a narrow gravitational trough that extends the length of the Pliocene and Quaternary volcanic arc from Mount Baker in northern Washington to Lassen Peak in California. Gravitational highs interrupt the trough at two localities: a northwest trending high in southern Washington and a northeast trending high between Mount Shasta and Lassen Peak. The latter anomaly is one of a set of northeast trending anomalies that, within the Quaternary arc, appear related to volcanic segmentation proposed previously on the basis of spatial and compositional distributions of volcanoes. These northeast trending anomalies extend hundreds of kilometers northeast of the arc, are caused by sources in the upper crust, and in some cases are related to exposed pre-Tertiary rocks. Segmentation models invoke geometric characteristics of the subducting plate as the primary factor controlling location and chemistry of volcanism, and these northeast trending gravity sources also may be a product of disturbance of the upper crust by the subduction process. More likely, the gravity sources may reflect upper crustal structures older than the High Cascades, possibly relicts from earlier accretionary events or more recent crustal deformation, that have actively influenced the spatial location of more recent volcanism. Much of the Pliocene and Quaternary volcanism of the Cascade arc has concentrated on or near contacts between crustal blocks of disparate density. These contacts may promote the ascension of magma to the Earth's surface.

Journal of Geophysical Research Solid Earth↗
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 bluffs. 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 bluff retreat to assess the relative roles of the length of sea ice-free season, sea level, water temperature, nearshore wavefield, 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 bluff 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 wavefield. 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 bluffs will increase rates of coastal erosion. The certain increases in sea level and potential changes in storminess will both contribute to this effect. As water temperature also influences 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↗