Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical 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 1,657 records · Page 92Linked to original sources
USGS2012

Geology of the Mount St. Helens area: Record of discontinuous volcanic and plutonic activity in the Cascade Arc of southern Washington

The Quaternary edifice of Mount St. Helens volcano was built upon a deeply eroded terrane of gently folded and altered volcanic and plutonic rocks that represent the core of the Tertiary Cascade magmatic arc. These rocks constitute an east dipping homoclinal sequence, several kilometers thick, of subaerially erupted mafic to silicic flows and volcaniclastic strata; K‐Ar ages from this section range from about 28 to 23 Ma (late Oligocene and earliest Miocene), which corresponds to an apparent lull in Cascade volcanism to the north of Mount Rainier. Volcanism was essentially continuous during this period of time, and neither a well‐defined base nor top of the section is exposed within the mapped area. Basalt and basaltic andesite dominate the lower part of the mid‐Tertiary section, whereas andesitic and dacitic rocks comprise most of the upper part. This section was intruded by numerous mafic to silicic dikes, sills, and irregular plutonic bodies, most no more than a few million years younger than their host rocks, and subjected to pervasive burial metamorphism and widespread hydrothermal alteration. Large areas of hornfelsed rock surrounding even relatively small intrusions indicate that the proportion of plutonic rock becomes significantly greater at shallow depth beneath the existing erosion surface. A large granitic pluton intruded the mid‐Tertiary section north of Spirit Lake at about 21 Ma. The Earl porphyry copper deposit occurs within the pluton but appears too young (17 Ma) to be genetically related to it. In contrast to the rather continuous and voluminous Oligocene to early Miocene activity, volcanism since then in the Mount St. Helens area has been localized and Volumetrically minor. Products of three younger eruptive periods have been recognized: a sequence of 15 m.y. old pyroxene andesite flows resting unconformably on mid‐Tertiary strata south of Mount St. Helens, widespread shallow dikes and sills of pyroxene andesite between 10 and 8 m.y. old, and compositionally diverse rocks erupted during the past 3 m.y. The Quaternary lavas are more potassic than the Tertiary lavas and typically contain phenocrysts of hornblende and biotite, which are absent from the older rocks. A number of Tertiary structures define a broad NNE trending zone that may reflect a deep‐seated lithospheric flaw that has controlled the locus of Cascade magmatism in southern Washington for the past 25 m.y. Mount St. Helens lies within this zone at the intersection of the NNW striking St. Helens seismic zone (SHZ) and an ENE trending alignment of Pleistocene silicic plug‐domes. No surface breakage has been detected along the SHZ, which is apparently very young. The linear zone of silicic vents is probably controlled by a fault that has been interpreted from seismic records to occur directly beneath the volcano. This zone parallels the directions of regional maximum horizontal compressive stress and North America/Juan de Fuca plate convergence. Mount St. Helens is an example of a low‐volume tectonically controlled magmatic system in an early stage of development.

Washington↗
USGS2012

Tectonics and conductivity structures in the Southern Washington Cascades

The tectonic setting of the southern Washington Cascades has been studied with the aid of magnetotelluric (MT) and other geophysical data. The main feature of interest in the geophysical data is a broad high-conductivity anomaly mapped with MT and geomagnetic variation (GMV) data. This anomaly is located roughly within the triangle formed by the volcanoes Mount Rainier, Mount St. Helens, and Mount Adams but exceeds beyond Mount Rainier to the northwest. We interpret the cause of the anomaly to be conductive rocks with resistivities of 1–4 ohm m and thicknesses possibly greater than 15 km. These conductive rocks are found 2–8 km beneath the overlying less conductive volcanic and sedimentary rocks at the surface. Two aeromagnetic lows follow the trend of the conductivity anomaly, and linear belts of strike-slip seismicity are coincident with both these magnetic lows. One of the aeromagnetic lows is coincident with the western margin of the conductivity anomaly. The geophysical data appear to outline a suture zone of probable Eocene age caused by accretion of a large seamount complex (Siletzia) and that may contain large thicknesses of compressed forearc basin and accretionary prism sedimentary rocks of Cretaceous to Eocene age. Part of the shallower conductive rocks may be associated with carbonaceous continental and transitional marine sedimentary rocks of the Puget Group. The contact between the hypothesized compressed basin and the accreted terranes to the west may localize the release of shear stresses in this region of oblique subduction. Several possible explanations for the conductivity anomaly are considered in addition to the compressed basin hypothesis.

Washington↗
USGS2012

An aeromagnetic study of Mount St. Helens

Aeromagnetic data from surveys flown by the U.S. Geological Survey over Mount St. Helens, Washington, before and after the climactic May 18, 1980, eruption were used to determine the bulk magnetic properties of the volcano and to delineate a buried source. We assumed that most of the edifice of preeruption and posteruption Mount St. Helens has a magnetization direction near the present earth's field of 69° inclination and 20° declination and calculated its intensity as 4.1 A/m. After subtraction of magnetic anomalies due to topogaphy magnetized with this direction and intensity, the preeruption and posteruption surveys revealed nearly identical residual magnetic highs and lows, indicative that their sources were not altered or removed by the May 1980 eruption. The residual highs were explained by a 200-m-deep source lying mostly within the edifice of Mount St. Helens. The source could be terrain that predates Mount St. Helens, such as a buried ridge or a cone or a valley filled with lava. We calculated the magnetization of the material removed by the May 1980 eruption and found its intensity to be 4.2 A/m in a direction near the present earth's field, similar to that assumed for the volcano as a whole. This similar result confirmed the validity of the magnetization vector assumed for the entire edifice.

Washington↗
USGS2012

Evidence from gravity data for an intrusive complex beneath Mount St. Helens

On the basis of gravity data, aided by aeromagnetic, magnetotelluric, side-looking radar, and geologic information, we tentatively identify a large, shallow intrusive complex beneath Mount St. Helens. The complex is roughly 5–6 km thick and has apparently intruded a buried and compressed Mesozoic forearc sedimentary sequence. The lateral extent of the intrusive complex is uncertain, and we give three alternative models varying from about 18 by 10 km to as much as 18 by 22 km. The western boundary of the inferred Mount St. Helens intrusion abuts several of a number of Tertiary and younger plutons that crop out in the area. The Mount St. Helens intrusion apparently is adjacent to sedimentary rock or shallow volcanic rock along most of the remainder of its boundary. A ringlike drainage pattern around Mount St. Helens suggests subsidence caused by removing magma from deep beneath Mount St. Helens and adding this mass to the Mount St. Helens edifice and its underlying intrusion. Our analysis indicates that the average bulk density of the volcanic rock comprising Mount St. Helens is about 2.15 g/cm 3 .

Washington↗
USGS2012

Geochemistry, mineralogy, and petrogenesis of basalt from the Gorda Ridge

Basalt pillow lava with glassy rims was dredged from 17 sites along the Gorda Ridge between latitudes 43°N and 40.8°N. All samples are low-K 2 O mid-ocean ridge tholeiitic basalt having a narrow compositional range with Mg numbers clustered around 60–62; more primitive and evolved compositions are present but rare. None of the more fractionated samples is as evolved as ferrobasalt from the Juan de Fuca Ridge. Incompatible trace element ratios are typical of normal mid-ocean ridge basalt but show considerable scatter. The 87 Sr/ 86 Sr ratios, ranging from 0.70233 to 0.70267, are near the low end of the range for normal mid-ocean ridge basalt. The phenocrysts present are abundant plagioclase, less abundant olivine ± spinel, and traces of sulfide. Clinopyroxene is absent in the glasses. The mineral chemistry is typical of ocean floor basalt. However, unusually anorthitic plagioclase in some highly plagioclasephyric samples, complex zoning in plagioclase and olivine, and other disequilibrium features indicate complex petrogenetic processes involving magma mixing. Glass inclusions in olivine and plagioclase indicate primitive trapped melts. The chemical diversity of Gorda Ridge lava appears to have been primarily caused by variable degrees of partial melting, but crystal fractionation and magma mixing were also important processes. Compositional variations along the strike of the ridge indicate that lava from the segment south of the offset at latitude 41.6°N was generated by the smallest percentage of partial melting of a source less depleted than that of the two northern segments. Lava from the central and northern segments was generated by variable percentages of partial melting of a uniformly depleted source. Lava from the northern segment shows the greatest diversity with a large range in percentage of melting and both relatively primitive and evolved compositions in spatially and temporally closely associated lava. The diversity of lava from the northern segment may result from thermal perturbations caused by proximity to the Blanco Fracture Zone. Gorda Ridge lava is more similar to that from the Mid-Atlantic Ridge at 22°–25°N than to that from other spreading centers in the Pacific Ocean. Magma chambers under the Gorda Ridge appear to be small, discontinuous in time and space and possibly at greater depth than those beneath other Pacific spreading centers.

Journal of Geophysical Research Solid Earth↗
USGS2012

Effect of crustal layering upon dislocation modeling

Slip distribution at depth on a fault may be inferred from the deformation observed on the surface. In inverting the surface deformation data to obtain the slip distribution, the Earth is generally approximated by an elastic half‐space. Slip distributions inferred from a half‐space model may contain artifacts, including zones of reversed slip, due solely to effects of layering in the real Earth. This effect is demonstrated for a vertical strike‐slip fault in an Earth consisting of an elastic layer overlying an elastic half‐space. Slip on the fault is taken to be independent of the along‐strike coordinate (i.e., antiplane strain is assumed). For a given slip distribution in this model the slip distribution on a similar fault in an elastic half‐space is found that produces the identical surface deformation. Comparison of the two slip distributions reveals structure introduced into the half‐space equivalent slip profile by crustal layering. The comparisons suggest that low‐resolution inversion schemes (e.g., single screw dislocation models) are not drastically affected by Earth structure, but attempts at detailed inversion are likely to produce profiles contaminated by artifacts of Earth structure.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Morphometric variability within the axial zone of the southern Juan de Fuca Ridge: Interpretation from Sea MARC II, Sea MARC I, and deep-sea photography

The morphometric characteristics of the axial regions of oceanic spreading centers are determined by (1) the type of volcanic flows, (2) the relation between primary volcanic relief (on a scale of a few meters to tens of meters) and degree of sediment cover, and (3) the extent of surficial expression and timing of tectonic disruption of the young oceanic crust. Even within a single, continuous, linear spreading-ridge segment with relatively uniform axial valley dimensions over a distance of 50 or more kilometers, such as along the southern Juan de Fuca Ridge, the changes in morphometric characteristics along axis within the youngest crust indicate distinct variation in tectonic and volcanic activity over short distances within short time periods. An integrated analysis of Sea MARC I, Sea MARC II, and photographic data for the southernmost continuous segment of the Juan de Fuca Ridge shows that generalizations about tectonic and volcanic processes at spreading ridges must consider both the temporal scale of processes as well as the physical scales of observations if predictive models are to be successful. Comparison of the morphometric expression within the major hydrothermal vent area and the rest of the southernmost ridge segment suggests that the mapped distribution of hydrothermal vents may reflect the extent of survey effort rather than uniqueness of geologic setting.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Structure of the southern Juan de Fuca Ridge from seismic reflection records

Twenty-four-channel seismic reflection records were obtained from the axial region of the southern Juan de Fuca Ridge. Two profiles are normal to the strike of the spreading center and intersect the ridge at latitude 44°40′N and 45°05′N; a third profile extends south along the ridge axis from latitude 45°20′N and crosses the Blanco Fracture Zone. Processing of the axial portions of the cross-strike lines resolved a weak reflection centered beneath the axis. The reflector is at a depth similar to seismically detected magma chambers on the East Pacific Rise and a Lau Basin spreading center; we suggest that the reflector represents the top of an axial magma chamber. In the migrated sections the top of the probable magma chamber is relatively flat and 1–2 km wide, and the subbottom depth of the chamber is greater where the depth to the ridge axis is greater.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Uniformity and diversity in the composition of mineralizing fluids from hydrothermal vents on the southern Juan de Fuca Ridge

Abundances of Li, Na, K, Rb, Ca, Sr, Ba, Mn, Fe, Zn, and Si have been determined in fluid samples from seven vents located in three areas on the southern Juan de Fuca Ridge. The hydrothermal component estimated from the Mg contents of the samples ranges from 7% to 76%. Concentrations of Fe and Si, among other elements, in acid-stabilized solutions appear to be generally representative of the parental hydrothermal fluids, but some Zn determinations and most Ba values appear to be too low. Thermodynamic calculations indicate that the acidified samples remain supersaturated with respect to silica, barite, and pyrite; unacidified samples are supersaturated, in addition, with respect to ZnS, FeS, and many silicate phases. Within the constraints of limited sampling there appear to be differences in fluid compositions both within and between the three vent areas. Some uniform differences in the elemental abundances predicted for hydrothermal end-member fluids might be due to inmixing of fresh seawater at depth in the hydrothermal system. The Juan de Fuca hydrothermal fluids contain more Fe but otherwise have relative elemental abundances fairly similar to those in 13°N (East Pacific Rise) fluids, albeit at higher levels. In contrast, fluids from 21°N (East Pacific Rise) and Galapagos have lower K/Rb and much lower Sr and Na abundances; these compositional features probably result from interaction of these fluids with a different mineral assemblage, possibly more mature greenstone.

Journal of Geophysical Research Solid Earth↗
USGS2012

Chemistry of hydrothermal solutions from the southern Juan de Fuca Ridge

Fluids from three vent fields on the southern Juan de Fuca Ridge were sampled in September 1984 using the DSRV Alvin. The fluids are uniquely high in both chloride, which ranges up to twice the seawater value, and in metal content. Simple vapor-liquid phase separation could not have produced both the high chlorinity and gas concentrations observed in these fluids. The cause of the elevated chlorinity can not be uniquely identified but may be the result of either or a combination of two processes: (1) subsurface formation of a degassed brine and subsequent mixing of a small amount of this brine with a hydrothermal seawater which has not undergone a phase separation and (2) dissolution of a chloride-rich phase combined with a possible small loss of gas during sampling procedures. Although measured temperatures were all less than 300°C, quartz geothermometry suggests that the fluids have equilibrated at greater than 340°C. Quartz geobarometry is also in agreement with geophysical estimates of depth to the local magma chamber.

Journal of Geophysical Research B: Solid Earth↗
USGS2012

Mineralogical studies of sulfide samples and volatile concentrations of basalt glasses from the southern Juan de Fuca Ridge

Sulfide samples obtained from the U.S. Geological Survey's DSRV Alvin dives on the southern Juan de Fuca Ridge closely resemble those from the same area described by Koski et al. (1984). Major minerals include sphalerite, wurtzite, pyrite, marcasite, isocubanite, anhydrite, and chalcopyrite. Equilibrium, if attained at all, during deposition of most sulfides was a transient event over a few tens of micrometers at most and was perturbed by rapid temperature and compositional changes of the circulating fluid. Two new minerals were found: one, a hydrated Zn, Fe hydroxy-chlorosulfate, and the other, a (Mn, Mg, Fe) hydroxide or hydroxy-hydrate. Both were formed at relatively low temperatures. Lizardite, starkeyite, and anatase were found for the first time in such an environment. Sulfide geothermometry involving the system Cu-Fe-S indicates a vent temperature of <328°C for one sample. Fluid inclusion studies on crystals from the same vicinity of the same sample give pressure-corrected homogenization temperatures of 268° and 285°C. Ice-melting temperatures on inclusions from the same sample are about −2.8°C, indicating that the equivalent salinity of the trapped fluid is about 50% greater than that of seawater. Volatile concentrations from vesicle-free basalt glass from the vent field are about 0.013 wt % CO 2 and 0.16 wt % H 2 O. CO 2 contents in these samples yield an entrapment depth of 2200 m of seawater, which is the depth from which the samples were collected.

Journal of Geophysical Research Solid Earth↗
USGS2012

Fitting degradation of shoreline scarps by a nonlinear diffusion model

The diffusion model of degradation of topographic features is a promising means by which vertical offsets on Holocene faults might be dated. In order to calibrate the method, we have examined present-day profiles of wave-cut shoreline scarps of late Pleistocene lakes Bonneville and Lahontan. It may be assumed that these scarps were initially at least as steep as the angle of repose. Offsets range from 1 to 12 m, and present slope angles range from 9° to 29°. A parameter called apparent diffusion age, defined as half the mean square horizontal extent of the slope function of each profile, is plotted as a function of scarp offset. The points show a clear trend of apparent age increasing nearly linearly with offset. If linear diffusion held and scarps were initially vertical, apparent diffusion age would be the same for all the shoreline profiles. The increasing trend can only partly be explained by nonvertical initial scarp slope, and therefore the rate of transport of material downslope must increase significantly faster than a linear law in the range of slopes spanned by the data. The transport law must become linear at small slope to reduce scatter between profiles with varying ambient slopes. The transport law adopted for the purpose of dating is k 0 (1 + 5 s 2 ), where s is local slope. The transport coefficient k 0 is correlated inversely with fan slope, suggesting that there is a dependence on the particle size distribution. A table is included that allows easy application of the model to scarps with simple initial shape.

Journal of Geophysical Research Solid Earth↗
USGS2012

Anomalous shear wave attenuation in the shallow crust beneath the Coso volcanic region, California

We use seismograms of local earthquakes to image relative shear wave attenuation structure in the shallow crust beneath the region containing the Coso volcanic-geothermal area of eastern California. SV and P wave amplitudes were measured from vertical component seismograms of earthquakes that occurred in the Cososouthem Sierra Nevada region from July 1983 to August 1985. Seismograms of 16 small earthquakes show SV amplitudes which are greatly diminished at some azimuths and takeoff angles, indicating strong lateral variations in S wave attenuation in the area. Three-dimensional images of the relative S wave attenuation structure are obtained from forward modeling and a back projection inversion of the amplitude data. The results indicate regions within a 20 by 30 by 10 km volume of the shallow crust (one shallower than 5 km) that severely attenuate SV waves passing through them. These anomalies lie beneath the Indian Wells Valley, 30 km south of the Coso volcanic field, and are coincident with the epicentral locations of recent earthquake swarms. No anomalous attenuation is seen beneath the Coso volcanic field above about 5 km depth. Geologic relations and the coincidence of anomalously slow P wave velocities suggest that the attenuation anomalies may be related to magmatism along the eastern Sierra front.

Journal of Geophysical Research Solid Earth↗
USGS2012

Volcanic hazards and public response

Although scientific understanding of volcanoes is advancing, eruptions continue to take a substantial toll of life and property. Some of these losses could be reduced by better advance preparation, more effective flow of information between scientists and public officials, and better understanding of volcanic behavior by all segments of the public. The greatest losses generally occur at volcanoes that erupt infrequently where people are not accustomed to dealing with them. Scientists sometimes tend to feel that the blame for poor decisions in emergency management lies chiefly with officials or journalists because of their failure to understand the threat. However, the underlying problem embraces a set of more complex issues comprising three pervasive factors. The first factor is the volcano: signals given by restless volcanoes are often ambiguous and difficult to interpret, especially at long-quiescent volcanoes. The second factor is people: people confront hazardous volcanoes in widely divergent ways, and many have difficulty in dealing with the uncertainties inherent in volcanic unrest. The third factor is the scientists: volcanologists correctly place their highest priority on monitoring and hazard assessment, but they sometimes fail to explain clearly their conclusions to responsible officials and the public, which may lead to inadequate public response. Of all groups in society, volcanologists have the clearest understanding of the hazards and vagaries of volcanic activity; they thereby assume an ethical obligation to convey effectively their knowledge to benefit all of society. If society resists, their obligation nevertheless remains. They must use the same ingenuity and creativity in dealing with information for the public that they use in solving scientific problems. When this falls short, even excellent scientific results may be nullified.

Journal of Geophysical Research Solid Earth↗
USGS2012

The mechanics and three-dimensional internal structure of active magmatic systems: Kilauea volcano, Hawaii

Interpretation of abundant seismic data suggests that Kilauea's primary conduit within the upper mantle is concentrically zoned to about 34-km depth. This zoned structure is inferred to contain a central core region of relatively higher permeability, surrounded by numerous dikes that are in intermittent hydraulic communication with each other and with the central core. During periods of relatively high magma transport, the entire cross section of the conduit is utilized. During periods of relatively low to moderate transport, however, only the central core is active. As the conduit penetrates the oceanic crust and enters the volcanic shield, it simultaneously supplies the deeper sections of the rift zones (6-to 10-km depth) and the roots of the summit reservoir with picritic magma. The rift zones at depth are inferred to be almost wholly molten and to possess a high degree of fluid continuity from Heiheiahulu in the East Rift Zone, 45 km westward through the roots of the summit magma reservoir, and well into the Southwest Rift Zone. Higher in the shield, the subcaldera magma reservoir and the shallow rift zones occupy the 2-to 4-km depth interval. Summit-differentiated olivine tholeiite (ρ ≈ 2.62 g cm −3 ) is periodically injected laterally along a horizon of neutral buoyancy within the rift zones, where the density of the magma is just balanced by the in situ density of the shield (Ryan, 1987 a , b ). Deep rift zone intrusions push seaward the deep tectonic blocks of the volcano's south flank. Shallow rift intrusions build a sheeted dike complex, inferred to be in isostatic equilibrium with the higher-density deep rift cores below. General finite element analyses are presented for the deformation and stress fields surrounding such dikes in the horizontal and vertical planes. The dike tip in two and three dimensions is surrounded by a tubular core of tensile (σ 1 , σ 2 ) and shear stress (τ max ). The displacement field is characterized by counterrotating cells on either side of the dike tip which, in vertical orientation, produce the characteristic subsidence above the dike complex, with uplift on either side, forming a ridge-trough-ridge structure. A finite element model of Kilauea's shield computes the displacement fields and principal stress (σ 1 ) distributions resulting from intrusive activity on each or both of the rift zones. Within the summit region, tensile stress lobes produced by the three-dimensional upward extension of the intrusions superpose constructively to produce calderawide regimes of tensile stress, conducive to caldera development. Parametric studies of (1) intrusion in the East Rift Zone only, (2) intrusion in the Southwest Rift Zone only, and (3) intrusion in both rift zones demonstrate their unique kinematic contributions. For case 1, the caldera undergoes a counterclockwise rotation (torque up state) conducive to the development of rightstepping en echelon eruptive fissures, as exemplified by the August 14, 1971, eruption. For case 2, the caldera undergoes a clockwise rotation (torque down state) conducive to the development of left-stepping eruptive fissures, as occurred during the December 31, 1974, eruption. For case 3, the caldera substructure is driven due southward, producing the southward migration of the upper portions of the summit magma reservoir.

Journal of Geophysical Research Solid Earth↗
USGS2012

The giant submarine alika debris slide, Mauna Loa, Hawaii

A 4000‐km 2 area of submarine slump and slide deposits along the west flank of Mauna Loa volcano has been mapped with GLORIA side‐scan sonar images, seismic reflection profiles, and new bathymetry. The youngest deposits are two debris avalanche lobes that travelled from their breakaway area near the present shoreline as much as 100 km into the Hawaiian Deep at water depths of 4800 m. The two lobes partly overlap and together are designated the Alika slide. They were derived from the same source area and probably formed in rapid succession. Distinction hummocky topography, marginal levees, and other features on lower slopes (0.3°–0.6°) of these deposits resemble subaerial volcanic debris avalanche deposits such as 1980 Mount St. Helens and suggest high emplacement velocities. The breakaway area for the Alika slide (10°–15° slopes) is characterized by large block slumps, bounded by normal faults, that probably represent multiple subsidence events before, during, and after the debris avalanches. Lower slopes of the slide contain distinctive lobate‐terraced deposits that are interpreted as having been emplaced more slowly, prior to the debris avalanches. Estimated thicknesses of 50–200 m suggest volumes of 200–600 km 3 for the two lobes. The combined volume of the entire slide and slump terrane is probably 1500–2000 km 3 . The slide deposits predate a 13‐ka coral reef and probably postdate the block‐faulted Ninole Basalt, roughly dated as a few hundred thousand years old. The Alika slide, or a similar deposit recognized on GLORIA images further north along the Hawaiian Ridge, probably triggered a giant wave that washed 325 m high on Lanai at about 100 ka. Slumping on Mauna Loa has been most intense adjacent to the large arcuate bend in its southwest rift zone, as the rift zone migrated westward away from the growing Kilauea volcano. Slumping events were probably triggered by seismic activity accompanying dike injection along the rift zone. Such massive slumps, landslides, and distal submarine turbidity flows appear to be widespread on the flanks of Hawaiian volcanoes.

Journal of Geophysical Research↗
USGS2012

The mechanics of ground deformation precursory to dome-building extrusions at Mount St. Helens 1981-1982

Detailed monitoring at Mount St. Helens since 1980 has enabled prediction of the intermittent eruptive activity (mostly dome growth) with unprecedented success. During 1981 and 1982, accelerating deformation of the crater floor around the vent (including radial cracks, thrust faults, and ground tilt) was the earliest indicator of impending activity. Numerical experiments using the finite element method to model the mechanics of the crater floor show that all deformation features can be explained by a uniform shear-stress boundary condition along the conduit wall. The magnitude of the shear stress required to match observed displacements (1–7 MPa) is inversely proportional to the conduit diameter (estimated to be 25–100 m). The most probable source of this shear stress is the flow of viscous magma up the conduit and into the lava dome. We propose a model in which the accelerating deformation, beginning as much as 4 weeks before extrusions, is caused by the increasing velocity of ascending magma in the conduit. This model is examined by using deformation data of the dome before four extrusions in 1981 and 1982 to estimate the volumetric flow rate through the conduit. This flow rate and an estimate of the effective viscosity of the magma enable calculation of an ascent velocity and an applied shear stress that, again, depend on the conduit diameter. The results of these calculations are consistent with the finite element experiments and show that the proposed model is feasible. In light of this model, events observed just before or near the time extrusions began, such as reversals of ground tilt direction from outward to inward and the sudden decrease in the number of shallow earthquakes, may indicate an abrupt decrease of shear stress in the conduit. This could be explained by a decrease in either the ascent velocity, or the effective viscosity of the magma ascending through the shallow conduit, or both, near the time of extrusion. Precursory deformation like that measured at Mount St. Helens should be observable at similar volcanoes elsewhere because it is caused by the fundamental process of magma ascent.

Journal of Geophysical Research Solid Earth↗
USGS2012

Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor

A dynamic source model is presented, in which a three-dimensional crack containing a viscous compressible fluid is excited into resonance by an impulsive pressure transient applied over a small area ΔS of the crack surface. The crack excitation depends critically on two dimensionless parameters called the crack stiffness, C = ( b /μ)( L / d ), and viscous damping loss, F = (12η L )/(ρ ƒ d 2 α), where b is the bulk modulus, η is the viscosity, ρ ƒ is the density of the fluid, μ is the rigidity, α is the compressional velocity of the solid, L is the crack length, and d is the crack thickness. The first parameter characterizes the ability of the crack to vibrate and shapes the spectral signature of the source, and the second quantifies the effect of fluid viscosity on the duration of resonance. Resonance is sustained by a very slow wave trapped in the fluid-filled crack. This guided wave, called the crack wave, is similar to the tube wave propagating in a fluid-filled borehole; it is inversely dispersive, showing a phase velocity that decreases with increasing wavelength, and its wave speed is always lower than the acoustic velocity of the fluid, decreasing rapidly as the crack stiffness increases. The source spectrum shows many sharp peaks characterizing the individual modes of vibration of the crack; the variation of spectral shape, both in the number and width of peaks, is surprisingly complex, reflecting the interference between the lateral and longitudinal modes of resonance, as well as nodes for these modes. The far-field spectrum is marked by narrow-band dominant and subdominant peaks that reflect the interaction of the various source modes. The frequency of the dominant spectral peak radiated by the source is independent of the radiation direction. The frequency, bandwidth, and spacing of the resonant peaks are strongly dependent on the crack stiffness, larger values of the stiffness factor shifting these peaks to lower frequencies and decreasing their bandwidth. The excitation of a particular mode depends on the position of the trigger and on the extent of the crack surface affected by the pressure transient. Fluid viscosity decreases the amplitudes of the main spectral peaks, smears out the finer structure of the spectrum, and greatly reduces the duration of the radiated signal. The energy loss by radiation is stronger for high frequencies, producing a seismic signature that is marked by a high-frequency content near the onset of the signal and dominated by a longer-period component of much longer duration in the signal coda. Such signature is in harmony with those displayed by long-period events observed on active volcanoes and in hydrofracture experiments. The very low velocity which is possible in a crack with high stiffness ( C ≥ 100) also provides an attractive explanation for very long period tremor, such as type 2 tremor at Aso volcano, Japan, without the requirement of an unrealistically large magma container. The standing wave pattern set up on the crack surface by the sustained resonance in the fluid is observable in the near field of the crack, suggesting that the location and extent of the source may be estimated from the mapping of the pattern of nodes and antinodes seen in its vicinity. According to the model, the long-period event and harmonic tremor share the same source but differ in the boundary conditions for fluid flow and in the triggering mechanism setting up the resonance of the source, the former being viewed as the impulse response of the tremor generating system and the latter representing the excitation due to more complex forcing functions.

Journal of Geophysical Research Solid Earth↗