Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Geophysical Research - Atmospheres”

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 559 records · Page 31Linked to original sources
USGS2007

Three-dimensional geophysical mapping of rock alteration and water content at Mount Adams, Washington: Implications for lahar hazards

[1] Hydrothermally altered rocks, particularly if water saturated, can weaken stratovolcanoes, thereby increasing the potential for catastrophic sector collapses that can lead to far-traveled, destructive debris flows. Evaluating the hazards associated with such alteration is difficult because alteration has been mapped on few active volcanoes and the distribution and intensity of subsurface alteration are largely unknown on any active volcano. At Mount Adams, some Holocene debris flows contain abundant hydrothermal minerals derived from collapse of the altered edifice. Intense hydrothermal alteration significantly reduces the resistivity and magnetization of volcanic rock, and therefore hydrothermally altered rocks can be identified with helicopter electromagnetic and magnetic measurements. Electromagnetic and magnetic data, combined with geological mapping and rock property measurements, indicate the presence of appreciable thicknesses of hydrothermally altered rock in the central core of Mount Adams north of the summit. We identify steep cliffs at the western edge of this zone as the likely source for future large debris flows. In addition, the electromagnetic data identified water in the brecciated core of the upper 100–200 m of the volcano. Water helps alter the rocks, reduces the effective stress, thereby increasing the potential for slope failure, and acts, with entrained melting ice, as a lubricant to transform debris avalanches into lahars. Therefore knowing the distribution of water is also important for hazard assessments. Our results demonstrate that high-resolution geophysical and geological observations can yield unprecedented views of the three-dimensional distribution of altered rock and shallow pore water aiding evaluation of the debris avalanche hazard.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗
USGS2008

Constraining fault constitutive behavior with slip and stress heterogeneity

We study how enforcing self-consistency in the statistical properties of the preshear and postshear stress on a fault can be used to constrain fault constitutive behavior beyond that required to produce a desired spatial and temporal evolution of slip in a single event. We explore features of rupture dynamics that (1) lead to slip heterogeneity in earthquake ruptures and (2) maintain these conditions following rupture, so that the stress field is compatible with the generation of aftershocks and facilitates heterogeneous slip in subsequent events. Our three-dimensional fmite element simulations of magnitude 7 events on a vertical, planar strike-slip fault show that the conditions that lead to slip heterogeneity remain in place after large events when the dynamic stress drop (initial shear stress) and breakdown work (fracture energy) are spatially heterogeneous. In these models the breakdown work is on the order of MJ/m2, which is comparable to the radiated energy. These conditions producing slip heterogeneity also tend to produce narrower slip pulses independent of a slip rate dependence in the fault constitutive model. An alternative mechanism for generating these confined slip pulses appears to be fault constitutive models that have a stronger rate dependence, which also makes them difficult to implement in numerical models. We hypothesize that self-consistent ruptures could also be produced by very narrow slip pulses propagating in a self-sustaining heterogeneous stress field with breakdown work comparable to fracture energy estimates of kJ/M2. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗
USGS2007

A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

[1] We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green's functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake ( Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth.

Tottori↗
USGS2007

Measurements of storm and nonstorm circulation in the northern Adriatic: October 2002 Through April 2003

Fifteen bottom-mounted Acoustic Doppler Current Profilers were deployed from October 2002 through April 2003 in the northern Adriatic Sea. Average transport from the portion of the Western Adriatic Current (WAC) along the Italian slope was 0.1470 ?? 0.0043 Sv, punctuated by bursts of more than twice that amount during storm events. Monthly means were calculated with times of strong wind-driven circulation excluded. These suggest a 2002/2003 seasonal separation consisting of October, December through February, and March through April. An extreme Po River flood influenced November conditions making seasonal categorization difficult. October generally had more kinetic energy and more vertical structure than other months, and near-inertial waves were more frequent in April and October. The Eastern Adriatic Current (EAC)/WAC (i.e. inflow/outflow) system was clearly present in the means for all months. The cyclonic gyre north of the Po River was present October through February. Generally, in the WAC, over 50% of kinetic energy came from vertically uniform monthly mean flows. Elsewhere, eddy kinetic energy was stronger than mean kinetic energy with 10-40% contributions for vertically uniform monthly mean flows, 40-60% for vertically uniform monthly varying flows, and 10-30% for vertically varying monthly varying flows. Mean currents for bora storms indicate enhancement of the EAC/WAC and the cyclonic northern gyre, a shift toward Kvarner Bay in EAC direction, a circulation null point south of the Po, and double-gyre bifurcation of flow at Istria. Strengthening of both the EAC and WAC also occurs during sirocco storms. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2008

Sediment dispersal in the northwestern Adriatic Sea

Sediment dispersal in the Adriatic Sea was evaluated using coupled three-dimensional circulation and sediment transport models, representing conditions from autumn 2002 through spring 2003. The calculations accounted for fluvial sources, resuspension by waves and currents, and suspended transport. Sediment fluxes peaked during southwestward Bora wind conditions that produced energetic waves and strengthened the Western Adriatic Coastal Current. Transport along the western Adriatic continental shelf was nearly always to the south, except during brief periods when northward Sirocco winds reduced the coastal current. Much of the modeled fluvial sediment deposition was near river mouths, such as the Po subaqueous delta. Nearly all Po sediment remained in the northern Adriatic. Material from rivers that drain the Apennine Mountains traveled farther before deposition than Po sediment, because it was modeled with a lower settling velocity. Fluvial sediment delivered to areas with high average bed shear stress was more highly dispersed than material delivered to more quiescent areas. Modeled depositional patterns were similar to observed patterns that have developed over longer timescales. Specifically, modeled Po sediment accumulation was thickest near the river mouth with a very thin deposit extending to the northeast, consistent with patterns of modern sediment texture in the northern Adriatic. Sediment resuspended from the bed and delivered by Apennine Rivers was preferentially deposited on the northern side of the Gargano Peninsula, in the location of thick Holocene accumulation. Deposition here was highest during Bora winds when convergences in current velocities and off-shelf flux enhanced delivery of material to the midshelf. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗
USGS2007

Bed forms created by simulated waves and currents in a large flume

The morphology and evolution of bed forms created by combinations of waves and currents were investigated using an oscillating plate in a 4-m-wide flume. Current speed ranged from 0 to 30 cm/s, maximum oscillatory velocity ranged from 20 to 48 cm/s, oscillation period was 8 s (except for one run with 12 s period), and the median grain size was 0.27 mm. The angle between oscillations and current was 90°, 60°, or 45°. At the end of each run the sand bed was photographed and ripple dimensions were measured. Ripple wavelength was also determined from sonar images collected throughout the runs. Increasing the ratio of current to wave (i.e., oscillatory) velocity decreased ripple height and wavelength, in part because of the increased fluid excursion during the wave period. Increasing the ratio of current to waves, or decreasing the angle between current and waves, increased the three-dimensionality of bed forms. During the runs, ripple wavelength increased by a factor of about 2. The average number of wave periods for evolution of ripple wavelength to 90% of its final value was 184 for two-dimensional ripples starting from a flat bed. Bed form orientations at the end of each run were compared to four potential controlling factors: the directions of waves, current, maximum instantaneous bed shear stress, and maximum gross bed form normal transport (MGBNT). The directions of waves and of MGBNT were equally good predictors of bed form orientations, and were significantly better than the other two factors.

Journal of Geophysical Research C: Oceans↗
USGS2006

Nature and origin of the hematite-bearing plains of Terra Meridiani based on analyses of orbital and Mars Exploration rover data sets

The ∼5 km of traverses and observations completed by the Opportunity rover from Endurance crater to the Fruitbasket outcrop show that the Meridiani plains consist of sulfate‐rich sedimentary rocks that are largely covered by poorly‐sorted basaltic aeolian sands and a lag of granule‐sized hematitic concretions. Orbital reflectance spectra obtained by Mars Express OMEGA over this region are dominated by pyroxene, plagioclase feldspar, crystalline hematite (i.e., concretions), and nano‐phase iron oxide dust signatures, consistent with Pancam and Mini‐TES observations. Mössbauer Spectrometer observations indicate more olivine than observed with the other instruments, consistent with preferential optical obscuration of olivine features in mixtures with pyroxene and dust. Orbital data covering bright plains located several kilometers to the south of the landing site expose a smaller areal abundance of hematite, more dust, and a larger areal extent of outcrop compared to plains proximal to the landing site. Low‐albedo, low‐thermal‐inertia, windswept plains located several hundred kilometers to the south of the landing site are predicted from OMEGA data to have more hematite and fine‐grained olivine grains exposed as compared to the landing site. Low calcium pyroxene dominates spectral signatures from the cratered highlands to the south of Opportunity. A regional‐scale model is presented for the formation of the plains explored by Opportunity, based on a rising ground water table late in the Noachian Era that trapped and altered local materials and aeolian basaltic sands. Cessation of this aqueous process led to dominance of aeolian processes and formation of the current configuration of the plains.

Journal of Geophysical Research E: Planets↗
USGS2006

The rocks of Gusev Crater as viewed by the Mini-TES instrument

The Miniature Thermal Emission Spectrometer (Mini-TES) on board the Mars Exploration Rover Spirit is part of a payload designed to investigate whether a lake once existed in Gusev Crater. Mini-TES has observed hundreds of rocks along the rover's traverse into the Columbia Hills, yielding information on their distribution, bulk mineralogy, and the potential role of water at the site. Although dust in various forms produces contributions to the spectra, we have established techniques for dealing with it. All of the rocks encountered on the plains traverse from the lander to the base of the Columbia Hills share common spectral features consistent with an olivine-rich basaltic rock known as Adirondack Class. Beginning at the base of the West Spur of the Columbia Hills and across its length, the rocks are spectrally distinct from the plains but can be grouped into a common type called Clovis Class. These rocks, some of which appear as in-place outcrop, are dominated by a component whose spectral character is consistent with unaltered basaltic glass despite evidence from other rover instruments for significant alteration. The northwest flank of Husband Hill is covered in float rocks known as Wishstone Class with spectral features that can be attributed uniquely to plagioclase feldspar, a phase that represents more than half of the bulk mineralogy. Rare exceptions are three classes of basaltic "exotics" found scattered across Husband Hill that may represent impact ejecta and/or float derived from local intrusions within the hills. The rare outcrops observed on Husband Hill display distinctive spectral characteristics. The outcrop called Peace shows a feature attributable to molecular bound water, and the outcrop that hosts the rock called Watchtower displays a dominant basaltic glass component. Despite evidence from the rover's payload for significant alteration of some of the rocks, no unambiguous detection of crystalline phyllosilicates or other secondary silicates has been observed by Mini-TES. The mineralogical results supplied by Mini-TES provide no clear evidence that a lake once existed in Gusev Crater. Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2006

Erosion rates at the Mars Exploration Rover landing sites and long-term climate change on Mars

Erosion rates derived from the Gusev cratered plains and the erosion of weak sulfates by saltating sand at Meridiani Planum are so slow that they argue that the present dry and desiccating environment has persisted since the Early Hesperian. In contrast, sedimentary rocks at Meridiani formed in the presence of groundwater and occasional surface water, and many Columbia Hills rocks at Gusev underwent aqueous alteration during the Late Noachian, approximately coeval with a wide variety of geomorphic indicators that indicate a wetter and likely warmer environment. Two-toned rocks, elevated ventifacts, and perched and undercut rocks indicate localized deflation of the Gusev plains and deposition of an equivalent amount of sediment into craters to form hollows, suggesting average erosion rates of ???0.03 nm/yr. Erosion of Hesperian craters, modification of Late Amazonian craters, and the concentration of hematite concretions in the soils of Meridiani yield slightly higher average erosion rates of 1-10 nm/yr in the Amazonian. These erosion rates are 2-5 orders of magnitude lower than the slowest continental denudation rates on Earth, indicating that liquid water was not an active erosional agent. Erosion rates for Meridiani just before deposition of the sulfate-rich sediments and other eroded Noachian areas are comparable with slow denudation rates on Earth that are dominated by liquid water. Available data suggest the climate change at the landing sites from wet and likely warm to dry and desiccating occurred sometime between the Late Noachian and the beginning of the Late Hesperian (3.7-3.5 Ga). Copyright 2006 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2007

Evidence for montmorillonite or its compositional equivalent in Columbia Hills, Mars

During its exploration of the Columbia Hills, the Mars Exploration Rover "Spirit" encountered several similar samples that are distinctly different from Martian meteorites and known Gusev crater soils, rocks, and sediments. Occurring in a variety of contexts and locations, these "Independence class" samples are rough-textured, iron-poor (equivalent FeO ??? 4 wt%), have high Al/Si ratios, and often contain unexpectedly high concentrations of one or more minor or trace elements (including Cr, Ni, Cu, Sr, and Y). Apart from accessory minerals, the major component common to these samples has a compositional profile of major and minor elements which is similar to the smectite montmorillonite, implicating this mineral, or its compositional equivalent. Infrared thermal emission spectra do not indicate the presence of crystalline smectite. One of these samples was found spatially associated with a ferric sulfate-enriched soil horizon, possibly indicating a genetic relationship between these disparate types of materials. Compared to the nearby Wishstone and Watchtower class rocks, major aqueous alteration involving mineral dissolution and mobilization with consequent depletions of certain elements is implied for this setting and may be undetectable by remote sensing from orbit because of the small scale of the occurrences and obscuration by mantling with soil and dust.

Journal of Geophysical Research E: Planets↗
USGS2006

Overview of the Opportunity Mars Exploration Rover Mission to Meridiani Planum: Eagle Crater to Purgatory Ripple

The Mars Exploration Rover Opportunity touched down at Meridiani Planum in January 2004 and since then has been conducting observations with the Athena science payload. The rover has traversed more than 5 km, carrying out the first outcrop‐scale investigation of sedimentary rocks on Mars. The rocks of Meridiani Planum are sandstones formed by eolian and aqueous reworking of sand grains that are composed of mixed fine‐grained siliciclastics and sulfates. The siliciclastic fraction was produced by chemical alteration of a precursor basalt. The sulfates are dominantly Mg‐sulfates and also include Ca‐sulfates and jarosite. The stratigraphic section observed to date is dominated by eolian bedforms, with subaqueous current ripples exposed near the top of the section. After deposition, interaction with groundwater produced a range of diagenetic features, notably the hematite‐rich concretions known as “blueberries.” The bedrock at Meridiani is highly friable and has undergone substantial erosion by wind‐transported basaltic sand. This sand, along with concretions and concretion fragments eroded from the rock, makes up a soil cover that thinly and discontinuously buries the bedrock. The soil surface exhibits both ancient and active wind ripples that record past and present wind directions. Loose rocks on the soil surface are rare and include both impact ejecta and meteorites. While Opportunity's results show that liquid water was once present at Meridiani Planum below and occasionally at the surface, the environmental conditions recorded were dominantly arid, acidic, and oxidizing and would have posed some significant challenges to the origin of life.

Journal of Geophysical Research E: Planets↗
USGS2007

Visible and near-infrared multispectral analysis of rocks at Meridiani Planum, Mars, by the Mars Exploration Rover Opportunity

Multispectral measurements in the visible and near infrared of rocks at Meridiani Planum by the Mars Exploration Rover Opportunity's Pancam are described. The Pancam multispectral data show that the outcrops of the Burns formation consist of two main spectral units which in stretched 673, 535, 432 nm color composites appear buff- and purple-colored. These units are referred to as the HFS and LFS spectral units based on higher and lower values of 482 to 535 nm slope. Spectral characteristics are consistent with the LFS outcrop consisting of less oxidized, and the HFS outcrop consisting of more oxidized, iron-bearing minerals. The LFS surfaces are not as common and appear, primarily, at the distal ends of outcrop layers and on steep, more massive surfaces, locations that are subject to greater eolian erosion. Consequently, the HFS surfaces are interpreted as a weathering rind. Further inherent spectral differences between layer's and between different outcrop map units, both untouched and patches abraded by the rover's Rock Abrasion Tool, are also described. Comparisons of the spectral parameters of the Meridiani outcrop with a set of laboratory reflectance measurements of Fe3+-bearing minerals show that the field of outcrop measurements plots near the fields of hematite, ferrihydrite, poorly crystalline goethite, and schwertmannite. Rind and fracture fill materials, observed intermittently at outcrop exposures, are intermediate in their spectral character between both the HFS and LFS spectral classes and other, less oxidized, surface materials (basaltic sands, spherules, and cobbles). Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2007

Coordinated analyses of orbital and spirit rover data to characterize surface materials on the cratered plains of Gusev Crater, Mars

Comparison of the Mars Exploration Rover Spirit's Pancam (0.4 to 1.0 ??m) and Mars Express Observatoire pour la Mineralogie l'Eau, les Glaces et l'Activite?? (OMEGA) (0.4 to 2.5 ??m) spectral reflectance data over Spirit's traverses shows that Gusev cratered plains are dominated by nanophase ferric-oxide-rich dust covering weakly altered basaltic sands. This interpretation is also consistent with both observations from OMEGA data covering plains beyond the traverse region and interpretations of data from the other payload instruments on the Spirit Rover. OMEGA observations of relatively low albedo regions where dust has presumably been stripped by dust devils show negative spectral reflectance slopes from 1.5 to 2.5 ??m and moderately masked spectral features which are indicative of olivine or pyroxene. High-albedo regions north and south of the Spirit landing site have flat spectral reflectance slopes and few spectral features, although all spectra have a nanophase ferric-oxide absorption edge between 0.4 and 0.75 ??m. Comparison of THEMIS-derived thermal inertia values with OMEGA-derived spectral parameters shows that although the dust cover can be optically thick (0.4 to 2.5 ??m wavelength region) in some areas, it is not thick enough (???1 cm) to mask the thermal inertia of the underlying substrate for areas included in this study. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
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↗
USGS2006

Form drag in rivers due to small-scale natural topographic features: 2. Irregular sequences

The size, shape, and spacing of small-scale topographic features found on the boundaries of natural streams, rivers, and floodplains can be quite variable. Consequently, a procedure for determining the form drag on irregular sequences of different-sized topographic features is essential for calculating near-boundary flows and sediment transport. A method for carrying out such calculations is developed in this paper. This method builds on the work of Kean and Smith (2006), which describes the flow field for the simpler case of a regular sequence of identical topographic features. Both approaches model topographic features as two-dimensional elements with Gaussian-shaped cross sections defined in terms of three parameters. Field measurements of bank topography are used to show that (1) the magnitude of these shape parameters can vary greatly between adjacent topographic features and (2) the variability of these shape parameters follows a lognormal distribution. Simulations using an irregular set of topographic roughness elements show that the drag on an individual element is primarily controlled by the size and shape of the feature immediately upstream and that the spatial average of the boundary shear stress over a large set of randomly ordered elements is relatively insensitive to the sequence of the elements. In addition, a method to transform the topography of irregular surfaces into an equivalently rough surface of regularly spaced, identical topographic elements also is given. The methods described in this paper can be used to improve predictions of flow resistance in rivers as well as quantify bank roughness.

Journal of Geophysical Research F: Earth Surface↗
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↗
USGS2007

Seismic detection and analysis of icequakes at Columbia Glacier, Alaska

Contributions to sea level rise from rapidly retreating marine-terminating glaciers are large and increasing. Strong increases in iceberg calving occur during retreat, which allows mass transfer to the ocean at a much higher rate than possible through surface melt alone. To study this process, we deployed an 11-sensor passive seismic network at Columbia Glacier, Alaska, during 2004–2005. We show that calving events generate narrow-band seismic signals, allowing frequency domain detections. Detection parameters were determined using direct observations of calving and validated using three statistical methods and hypocenter locations. The 1–3 Hz detections provide a good measure of the temporal distribution and size of calving events. Possible source mechanisms for the unique waveforms are discussed, and we analyze potential forcings for the observed seismicity.

Journal of Geophysical Research↗