Search USGSSearch

SEARCH · Search USGS

Results for “Journal of Interpretation Research”

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 235 records · Page 13Linked to original sources

Identifying relationships between baseflow geochemistry and land use with synoptic sampling and R-mode factor analysis

The relationship between land use and stream chemistry is often explored through synoptic sampling of rivers at baseflow conditions. However, baseflow chemistry is likely to vary temporally and spatially with land use. The purpose of our study is to examine the usefulness of the synoptic sampling approach for identifying the relationship between complex land use configurations and stream water quality. This study compares biogeochemical data from three synoptic sampling events representing the temporal variability of baseflow chemistry and land use using R - mode factor analysis . Separate R - mode factor analyses of the data from individual sampling events yielded only two consistent factors. Agricultural activity was associated with elevated levels of Ca 2+ , Mg 2+ , alkalinity, and frequently K + , SO 2- 4 , and NO - 3 . Urban areas were associated with higher concentrations of Na + , K + , and Cl - . Other retained factors were not consistent among sampling events, and some factors were difficult to interpret in the context of biogeochemical sources and processes. When all data were combined, further associations were revealed such as an inverse relationship between the proportion of wetlands and stream nitrate concentrations. We also found that barren lands were associated with elevated sulfate levels. This research suggests that an individual sampling event is unlikely to characterize adequately the complex processes controlling interactions between land use and stream chemistry. Combining data collected over two years during three synoptic sampling events appears to enhance our ability to understand processes linking stream chemistry and land use.

Michigan

Petrology of the Paloma Valley ring complex, southern California batholith

The Paloma Valley ring complex is one of the numerous plutons that make up the Cretaceous southern California batholith. The complex is composite, consisting of (1) an older, single ring dike and two subsidiary short-arced inner dikes, and (2) a younger set of thin short-arced dikes largely inside the older ring dike. The older ring dike, composed of granodiorite-quartz monzonite, has nearly vertical walls and is elliptical in plan; its long axis (14 km) is oriented west-northwest. It was emplaced in and contains numerous fragments of gabbro. The more than 200 short-arced granitic pegmatite dikes, mainly ranging from 0.2 to 1 m in thickness, define a domal ring dike set with moderately to steeply dipping outer dikes and nearly horizontal inner dikes. The younger dikes cut both the older ring dike and the gabbro. Spatially associated with the younger dikes are bodies of fine-grained granophyre that contain stringers of granitic pegmatite. The granophyre has an Mg content similar to that of the younger ring-dike rock, but contains less K and more Fe. The older ring dike is interpreted to have been magmatically emplaced in an elliptical zone of ring fracturing in gabbro; the magma made room by stoping the denser gabbro with little assimilation. Upon a release of pressure, a set of domal fractures formed, along which volatile-rich magma was emplaced, forming the younger ring dikes. Granophyre resulted from pressure-quenching through a loss of volatiles. Residual volatiles, or volatiles that were introduced later, recrystallized parts of the granophyre and caused the formation of pegmatite stringers.

California

A 450-year record of environmental change from Castle Lake, California (USA), inferred from diatoms and organic geochemistry

A 39-cm sediment core from Castle Lake, California (USA) spans the last ~ 450 years and was analyzed for diatoms and organic geochemistry (δ 15 N, δ 13 C, and C:N), with the goal of determining sensitivity to natural climate variation and twentieth century anthropogenic effects. Castle Lake is a subalpine, nitrogen-limited lake with ~ 5 months of annual ice cover. Human impacts include light recreational use, past fish stocking, and experimental use by the Castle Lake Research Station. The base of the core (below 32 cm; pre mid-1700s) represents the period of maximum ice cover. In contrast, the end of the Little Ice Age (mid 1700s–early 1800s) is dominated by cyclotelloids (mostly Discostella stelligera ), indicating significant open-water periods, a condition that persisted into the early 1900s. Cyclotelloids began to decline in the 1960s and were replaced by the Fragilaria tenera grp. (peak in 1970s), succeeded by Asterionella formosa (peak ~ 2010), and accompanied by a reduction in δ 15 N values and a decrease in C:N that may represent increased atmospheric nitrogen deposition. Another anthropogenic signal was discerned in the core and was interpreted to be the result of an ammonium nitrate fertilization experiment of the epilimnion that was conducted in 1980 and 1981. This signal was manifested in the core largely by a negative excursion in δ 15 N, possibly caused by fractionation during denitrification in surface sediment. A phytoplankton monitoring dataset collected by the Castle Lake Research Station from 1967 to 1984 corroborates the timing of increased araphid euplanktonic species in the 1970s, and increases in two benthic diatoms ( Staurosirella pinnata and Tabellaria fenestrata) , entrained in the phytoplankton tows during the experimentation years. Both ice cover and nitrogen addition appear to be strong drivers that affected the lake diatoms, although additional drivers, such as fish stocking and associated cascade effects need further exploration. These data will be helpful for interpreting longer core records from Castle Lake, should the opportunity arise, as well as cores from similar systems in the region.

California

The 1989-1990 eruptions of Redoubt Volcano: an introduction

Redoubt Volcano, located on the west side of Cook Inlet in south-central Alaska, erupted explosively on over 20 separate occasions between December 14, 1989 and April 21, 1990. Fourteen lava domes were emplaced in the summit area, thirteen of which were subsequently destroyed. The eruption caused economic losses estimated at over $160,000,000 making this the second most costly eruption in U.S. history. This economic impact provided the impetus for a integrated comprehensive account of an erupting volcano using both modern and classical research and modern techniques which in turn led to advances in eruption monitoring and interpretation. Research on such topics as dome formation and collapse and the resulting pyroclastic flows, elutriated ash, lightning, tephra, and flooding was blended with the rapid communication of associated hazards to a large user group. The seismology successes in predicting and monitoring eruption dynamics were due in part to (1) the recognition of long-period seismic events as indicators of the readiness of the volcano to erupt, and (2) to the development of new tools that allowed the seismicity to be assessed instantaneously. Integrated studies of the petrology of erupted products and volatile content over time gave clues as to the progress of the eruption towards completion.

Alaska

Large-scale variation in lakebed properties interpreted from single-beam sonar in two Laurentian Great Lakes

Acoustic seabed classification (ASC) is an important method for understanding landscape-level physical and biological patterns in the aquatic environment. Bottom habitats in the Laurentian Great Lakes are poorly mapped to date, and will require a variety of contributors and data sources to complete. We repurposed a long-term split-beam echosounder dataset gathered for purposes of fisheries assessment to estimate lakebed properties utilizing unsupervised classification of echo return data. We interpreted first echo properties representing lakebed hardness and roughness to define and map three statistically supported lakebed classes revealed through cluster analysis. Our results indicate coherent and logical class boundaries and suggest that the dataset has promise for expanded use in ASC. To improve inferences using repeated measures, future work should focus on collecting ground truth information for areas previously surveyed and on collecting concurrent ground truth information when sampling acoustic data moving forward.

Michigan, Ontario, Wisconsin

Seismic reflection/refraction mapping of faulting and regional dips in the eastern Alaska Range

We present the results of a Trans‐Alaska Crustal Transect (TACT) investigation of the upper 2–5 km of the eastern Alaska Range in the vicinity of the Denali fault based on seismic reflection/refraction data, laboratory measurements of rock velocities, and structural mapping. The Denali fault is a major dextral slip structure mappable for more than 2000 km separating the Wrangellia and adjacent terranes to the south from the Yukon‐Tanana and adjacent terranes to the north. Geologic mapping suggests over 400 km of dextral slip has occurred on the fault, yet within the upper 1.5 km of the crust along the TACT corridor, basement rocks juxtaposed along the Denali fault reveal no significant seismic velocity differences, although the fault zone itself is associated with a minor lowering of velocity. The lack of seismic velocity contrast adjacent to the fault is in agreement with laboratory measurements of elastic wave velocities of samples from terranes bordering the fault. Laboratory measurements of elastic wave velocities of the metasedimentary mica‐quartz schists comprising the Yukon‐Tanana basement are highly anisotropic because of preferred orientation of mica and predict significant variations in velocity accompanying variations in foliation dip. Although other interpretations are possible, the northward shallowing of foliation dip of basement rocks in the Yukon‐Tanana terrane combined with the strong anisotropy associated with these highly foliated rocks can explain an observed northward increase in seismic velocity within this terrane. Seismic reflections from basement rocks within the Yukon‐Tanana terrane may originate from variations in anisotropy with depth and/or changes in composition reflecting different proportions of sandstone and shale in the protolith.

Alaska

Acoustic measurements of the 1999 basaltic eruption of Shishaldin volcano, Alaska 2. Precursor to the Subplinian phase

The 1999 eruption of Shishaldin volcano (Alaska, USA) displayed both Strombolian and Subplinian basaltic activity. The Subplinian phase was preceded by a signal of low amplitude and constant frequency (??? 2 Hz) lasting 13 h. This "humming signal" is interpreted as the coalescence of the very shallow part of a foam building up in the conduit, which produces large gas bubbles before bursting. The acoustic waveform of the hum event is modelled by a Helmholtz resonator: gas is trapped into a rigid cavity and can only escape through a tiny upper hole producing sound waves. At Shishaldin, the radius of the hole (??? 5 m) is close to that of the conduit (??? 6 m), the cavity has a length of ??? 60 m, and gas presents only a small overpressure between (??? 1.2 ?? 10-3 and 4.5 ?? 10-3 MPa). Such an overpressure is obtained by the partial coalescence of a foam formed by bubbles with a diameter from ??? 2.3 mm at the beginning of the episode towards ??? 0.64 mm very close to the end of the phase. The intermittency between hum events is explained by the ripening of the foam induced by the H2O diffusion through the liquid films. The two extreme values, from 600 to 10 s, correspond to a bubble diameter from 2.2 to 0.3 mm at the beginning and end of the pre-Subplinian phase, respectively. The extremely good agreement between two independent estimates of bubble diameters in the shallow foam reinforces the validity of such an interpretation. The total gas volume lost at the surface during the humming events is at most 5.9 ?? 106 m3. At the very end of the pre-Subplinian phase, there is a single large bubble with an overpressure of ???0.42 MPa. The large overpressure suggests that it comes from significant depth, unlike other bubbles in the pre-Subplinian phase. This deep bubble may be responsible for the entire foam collapse, resulting in the Subplinian phase. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research

Evidence for Milankovitch periodicities in Cenomanian-Turonian lithologic and geochemical cycles, western interior U.S.A.

The limestone/marlstone bedding couplets of the Bridge Creek Limestone Member, Cenomanian-Turonian Greenhorn Formation, were analyzed by applying spectral techniques to high-resolution lithologic and geochemical data from a core. The results suggest that the Bridge Creek contains a complex record of orbital cyclicity. The dominant signal appears to be obliquity, but signals corresponding to precession and eccentricity were also observed. The development of the bedding couplets is interpreted to have resulted from a combination of factors, including insolation-controlled changes in higher-latitude precipitation leading to dilution/redox cycles, and in lower-latitude evaporation, leading to changes in surface water conditions and productivity cycles in the calcareous plankton. The data interpreted to reflect redox cycles appear to be more strongly influenced by obliquity, and show a weak precessional signal. In contrast, trends in the carbonate record show the opposite response. The complex bedding pattern observed in the Bridge Creek Limestone is interpreted to result from the competing influences of different orbital cycles expressed through different pathways of the depositional system, and was also affected by changes in sedimentation rates related to relative sea level fluctuations, aperiodic dilution by volcanic ash, and changes in organic-matter production and redox conditions related to a global "oceanic anoxic event". These factors complicate cycle analysis in the lower part of the member but leave a relatively undisturbed record in the upper Bridge Creek Limestone.

Journal of Sedimentary Research

Using geologic structures to constrain constitutive laws not accessible in the laboratory

In this essay, we explore a central problem of structural geology today, and in the foreseeable future, which is the determination of constitutive laws governing rock deformation to produce geologic structures. Although laboratory experiments provide much needed data and insights about constitutive laws, these experiments cannot cover the range of conditions and compositions relevant to the formation of geologic structures. We advocate that structural geologists address this limitation by interpreting natural experiments, documented with field and microstructural data, using continuum mechanical models that enable the deduction of constitutive laws. To put this procedure into a historical context, we review the founding of structural geology by James Hutton in the late 18th century, and the seminal contributions to continuum mechanics from Newton to Cauchy that provide the tools to model geologic structures. The procedure is illustrated with two examples drawn from recent and on-going field investigations of crustal and mantle lithologies . We conclude by pointing to future research opportunities that will engage structural geologists in the pursuit of constitutive laws during the 21st century.

Journal of Structural Geology

Characteristics of and corrections for core shortening in unconsolidated sediments

Thinning, bypassing, and compaction of shallow unconsolidated sediments during manual coring or vibracoring operations probably cause more sediment deformation and greater stratigraphic displacement than is commonly reported in the wetland literature. We measured core shortening in open-barrel cores from fluvial wetlands, lagoonal flats, and marshes to document the magnitude and characteristics of shortening where sediments may be stiff and require extra mechanical effort to recover a sufficient length of sample for analysis. Results of those measurements indicate that thinning or non- recovery of discrete sediment intervals can range from 0 to 67 percent and cumulative core shortening can be as much as 30 percent even for cores less than one meter long. Detailed open-barrel measurements also show that core shortening is not uniformly distributed throughout the depth of penetration as is often assumed. Analytical data derived from shortened cores can only be properly interpreted if patterns of shortening are established and incorporated into the analysis. Minor artificial displacement of sediment depths can alter plots of physico-chemical parameters and can significantly influence calculated rates of sedimentation and other depth-dependent statistical relationships. This study (1) demonstrates how plots of interval shortening and cumulative shortening can be used to characterize the distribution of shortening at depth and (2) presents a simple equation for stratigraphic restoration so that core observations and analyses are corrected to their original depths.

Journal of Coastal Research

Landsat yesterday and today: An American vision and an old challenge

Since the late 1960s, the United States government has invested more than $1 billion in designing, launching, and operating the Landsat (land satellite) series of Earth-observing satellites. Global change researchers, geologists, and environmental scientists have used images gathered by the satellites for purposes ranging from human health research, energy exploration, and pollution detection to agricultural assessments, urban growth monitoring, and earthquake lineament studies. The earliest data were captured on a digital medium called wide-band video tape (WBVT). However, two decades of unsound media storage conditions and a poorly maintained processing system have left the physically deteriorating WBVTs with no mechanism for interpretation. A national treasure was in jeopardy. With seed money from the National Aeronautics and Space Administration (NASA), the U.S. Geological Survey (USGS) began a project to rescue the data. More than 21,000 tapes from the 1970s have been transcribed to stable, archival media, preserving the data for future studies in Earth System Science.

Journal of Map & Geography Libraries

Variations in the reflectivity of the moho transition zone beneath the Midcontinent Rift System of North America: results from true amplitude analysis of GLIMPCE data

True amplitude processing of The Great Lakes International Multidisciplinary Program on Crustal Evolution seismic reflection data from the Midcontinent Rift System of North America shows large differences in the reflectivity of the Moho transition zone beneath the axial rift, beneath the rift flanks, and outside of the rift. The Moho reflection from the axial rift has a discontinuous, diffractive character marginally stronger (several decibels) than an otherwise transparent lower crust and upper mantle. Beneath the axial rift, Moho is interpreted to be a synrift igneous feature. Beneath the rift flanks, the reflectivity of the Moho transition is generally well developed with two identifiable boundaries, although in places it is weakly reflective to nonreflective, similar to Moho outside the rift. The two boundaries are interpreted as the base of essentially intact, although stretched, prerift Archean crust (upper boundary) and new synrift Moho 1-2 s (6-7 km) deeper (lower boundary). Beneath the rift flanks, the layered reflection Moho transition results from the preexisting crustal composition and fabric modified by synrift igneous processes and extensional tectonic/metamorphic processes. The geologic evidence for extensive basaltic magmatism in the rift is the basis for interpreting the Moho signature as a Keweenawan structure that has been preserved for 1.1 b.y. Extension and magmatism appear to enhance reflectivity in the lower crust and Moho transition zone only where stretching factors are moderate (rift flanks) and not where they are extreme (axial rift). This leads to the prediction that the reflectivity across analogous volcanic passive continental margins should be greatest beneath the moderately stretched continental shelves and should decrease towards the ocean-continent boundary.

Journal of Geophysical Research B: Solid Earth

Interpretation of gravity data in a complex volcano-tectonic setting, southwestern Nevada

This regional gravity study, based on an irregular 2-km data grid, was conducted during the past few years at Yucca Mountain, southern Nye County, Nevada, as part of a program to locate a suitable repository for high-level nuclear waste. About 100 surface rock samples, three borehole gamma-gamma logs, and one borehole gravity study provide excellent density control. A nearly linear increase in density of 0.26 g/cm 3 per kilometer of depth is indicated in the thick tuff sequences that underlie the mountain. Isostatic and 2.0-g/cm 3 Bouguer corrections were applied to the observed gravity values to remove regional gradients and topographic effects, respectively. The Bare Mountain gravity high, with an isostatic anomaly maximum of 48 mGal, is connected with a greater gravity high over the Funeral Mountains, to the southwest; together, these highs result from a continuous block of dense, metamorphosed Precambrian and Paleozoic rocks that stretches across much of the Walker Lane from the east edge of Death Valley to Bare Mountain. The Calico Hills gravity high appears more likely to originate from a northeast trending buried ridge of Paleozoic rocks that extends southwestward beneath Busted Butte, 5 km southeast of the proposed repository, where two- and three-dimensional modeling indicates that the pre-Cenozoic rocks lie less than 1000 m beneath the surface. Tuff, at least 4000 m thick, fills a large steep-sided depression in the pretuff rocks beneath Yucca Mountain and Crater Flat. The gravity low and the thick tuff section lie within a large collapse area that includes the Crater Flat-Timber Mountain-Silent Canyon caldera complexes. Gravity lows in Crater Flat itself are interpreted to coincide with the source areas of the Prow Pass Member, the Bullfrog Member, and the Tram Member of the Crater Flat Tuff; these source areas add nearly 350 km 2 to the previously recognized extent of the local caldera complexes. Southward extension of the broad gravity low associated with Crater Flat into the Amargosa Desert is evidence for sector graben-type collapse segments related to the formation of the Timber Mountain caldera and superimposed on the other volcanic and extensional structures within Crater Flat.

Journal of Geophysical Research Solid Earth

Mineralogy and chemistry of cobbles at Meridiani Planum, Mars, investigated by the Mars Exploration Rover Opportunity

Numerous loose rocks with dimensions of a few centimeters to tens of centimeters and with no obvious physical relationship to outcrop rocks have been observed along the traverse of the Mars Exploration Rover Opportunity. To date, about a dozen of these rocks have been analyzed with Opportunity's contact instruments, providing information about elemental chemistry (Alpha Particle X‐ray Spectrometer), iron mineralogy and oxidation states (Mössbauer Spectrometer) and texture (Microscopic Imager). These “cobbles” appear to be impact related, and three distinct groups can be identified on the basis of chemistry and mineralogy. The first group comprises bright fragments of the sulfate‐rich bedrock that are compositionally and texturally indistinguishable from outcrop rocks. All other cobbles are dark and are divided into two groups, referred to as the “Barberton group” and the “Arkansas group,” after the first specimen of each that was encountered by Opportunity. Barberton group cobbles are interpreted as meteorites with an overall chemistry and mineralogy consistent with a mesosiderite silicate clast composition. Arkansas group cobbles appear to be related to Meridiani outcrop and contain an additional basaltic component. They have brecciated textures, pointing to an impact‐related origin during which local bedrock and basaltic material were mixed.

Journal of Geophysical Research E: Planets

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research

Impact crater outflows on Venus: Morphology and emplacement mechanisms

Many of the 932 impact craters discovered by the Magellan spacecraft at Venus are associated with lobate flows that originate at or near the crater rim. They extend for several to several hundred kilometers from the crater, and they commonly have a strong radar backscatter. A morphologic study of all identifiable crater outflows on Venus has revealed that many individual flows each consist of two areas, defined by distinct morphologic features. These two areas appear to represent two stages of deposition for each flow. The part of the flow that is generally deposited closest to the crater tends to be on the downrange side of the crater, flows in the downrange direction, and it is interpreted to be a late-stage ejecta. In many cases, this proximal part of the flow is too thin to completely bury the large blocks in subjacent ejecta deposits. Dendritic channels, present in many proximal flows, appear to have drained liquid from the proximal part in the downhill direction, and they debouch to feed the outer part of the flows. This distal part flows downhill, fills small grabens, and is ponded by ridges, behavior that mimics that of volcanic lava flows. The meandering and dendritic channels and the relation of the distal flows to topography strongly suggest that the distal portion is the result of coalescence and slow drainage of impact melt from the proximal portion. Impact melt forms a lining to the transient crater and mixes turbulently with solid clasts, and part of this mixture may be ejected to form the proximal part of the flow during the excavation stage of crater development. A statistical study of the Venusian craters has revealed that, in general, large craters produced by impacts with relatively low incidence angles to the surface are more likely to produce flows than small craters produced by higher-angle impacts. The greater flow production and downrange focusing of the proximal flows with decreasing incidence angle indicate a strong control of the flows by the impactor flight direction, and a high downrange velocity imparted to the proximal flow material in lower angle impacts. On the Moon, small flows interpreted to be composed of impact melt are observed atop the ejecta of large, fresh craters; on Earth, melt-rich suevite deposits form the uppermost layer of ejecta of some fresh craters. These features, albeit much smaller, may be analogous to the flows on Venus. Numerical models have predicted that larger volumes of impact melt would be produced on Venus than on the cooler terrestrial bodies due to high atmospheric and target temperatures, perhaps 3 times the volume produced on the Moon for a given crater diameter.

Journal of Geophysical Research E: Planets

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts

Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska

We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.

Alaska