USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal 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 361 records · Page 20Linked to original sources

The oligocene Lund Tuff, Great Basin, USA: A very large volume monotonous intermediate

Unusual monotonous intermediate ignimbrites consist of phenocryst-rich dacite that occurs as very large volume (> 1000 km3) deposits that lack systematic compositional zonation, comagmatic rhyolite precursors, and underlying plinian beds. They are distinct from countless, usually smaller volume, zoned rhyolite-dacite-andesite deposits that are conventionally believed to have erupted from magma chambers in which thermal and compositional gradients were established because of sidewall crystallization and associated convective fractionation. Despite their great volume, or because of it, monotonous intermediates have received little attention. Documentation of the stratigraphy, composition, and geologic setting of the Lund Tuff - one of four monotonous intermediate tuffs in the middle-Tertiary Great Basin ignimbrite province - provides insight into its unusual origin and, by implication, the origin of other similar monotonous intermediates. The Lund Tuff is a single cooling unit with normal magnetic polarity whose volume likely exceeded 3000 km3. It was emplaced 29.02 ?? 0.04 Ma in and around the coeval White Rock caldera which has an unextended north-south diameter of about 50 km. The tuff is monotonous in that its phenocryst assemblage is virtually uniform throughout the deposit: plagioclase > quartz ??? hornblende > biotite > Fe-Ti oxides ??? sanidine > titanite, zircon, and apatite. However, ratios of phenocrysts vary by as much as an order of magnitude in a manner consistent with progressive crystallization in the pre-eruption chamber. A significant range in whole-rock chemical composition (e.g., 63-71 wt% SiO2) is poorly correlated with phenocryst abundance. These compositional attributes cannot have been caused wholly by winnowing of glass from phenocrysts during eruption, as has been suggested for the monotonous intermediate Fish Canyon Tuff. Pumice fragments are also crystal-rich, and chemically and mineralogically indistinguishable from bulk tuff. We postulate that convective mixing in a sill-like magma chamber precluded development of a zoned chamber with a rhyolitic top or of a zoned pyroclastic deposit. Chemical variations in the Lund Tuff are consistent with equilibrium crystallization of a parental dacitic magma followed by eruptive mixing of compositionally diverse crystals and high-silica rhyolite vitroclasts during evacuation and emplacement. This model contrasts with the more systematic withdrawal from a bottle-shaped chamber in which sidewall crystallization creates a marked vertical compositional gradient and a substantial volume of capping-evolved rhyolite magma. Eruption at exceptionally high discharge rates precluded development of an underlying plinian deposit. The generation of the monotonous intermediate Lund magma and others like it in the middle Tertiary of the western USA reflects an unusually high flux of mantle-derived mafic magma into unusually thick and warm crust above a subducting slab of oceanic lithosphere. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

The proximal part of the giant submarine Wailau landslide, Molokai, Hawaii

The main break-in-slope on the northern submarine flank of Molokai at -1500 to -1250 m is a shoreline feature that has been only modestly modified by the Wailau landslide. Submarine canyons above the break-in-slope, including one meandering stream, were subaerially carved. Where such canyons cross the break-in-slope, plunge pools may form by erosion from bedload sediment carried down the canyons. West Molokai Volcano continued infrequent volcanic activity that formed a series of small coastal sea cliffs, now submerged, as the island subsided. Lavas exposed at the break-in-slope are subaerially erupted and emplaced tholeiitic shield lavas. Submarine rejuvenated-stage volcanic cones formed after the landslide took place and following at least 400-500 m of subsidence after the main break-in-slope had formed. The sea cliff on east Molokai is not the headwall of the landslide, nor did it form entirely by erosion. It may mark the location of a listric fault similar to the Hilina faults on present-day Kilauea Volcano. The Wailau landslide occurred about 1.5 Ma and the Kalaupapa Peninsula most likely formed 330??5 ka. Molokai is presently stable relative to sea level and has subsided no more than 30 m in the last 330 ka. At their peak, West and East Molokai stood 1.6 and 3 km above sea level. High rainfall causes high surface runoff and formation of canyons, and increases groundwater pressure that during dike intrusions may lead to flank failure. Active shield or postshield volcanism (with dikes injected along rift zones) and high rainfall appear to be two components needed to trigger the deep-seated giant Hawaiian landslides. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Analysis of long-period events recorded at Mount Etna (Italy) in 1992, and their relationship to eruptive activity

Seismic activity recorded at Mount Etna during 1992 was characterized by long-period (LP) events and tremor with fluctuating amplitudes. These signals were associated with the evolution of the eruptive activity that began on December 14, 1991. Following the occurrence of numerous volcano-tectonic earthquakes at the onset of the eruption, LP events dominated the overall seismicity starting in January, 1992. The LP activity occurred primarily in swarms, which were temporally correlated with episodic collapses of the crater floor in the Northeast Crater. Source depths determined for selected LP events suggest a source region located slightly east of Northeast Crater and extending from the surface to a depth of 2000 m. Based on the characteristic signatures of the time series, four families of LP events are identified. Each family shares common spectral peaks independent of azimuth and distance to the source. These spectral features are used to develop a fluid-filled crack model of the source. We hypothesize that the locus of the LP events represents a segment of the magma feeding system connecting a depressurizing magma body with a dike extending in the SSE direction along the western wall of Valle del Bove, toward the site of the Mount Etna eruption. We surmise that magma withdrawal from the source volume beneath Northeast Crater may have caused repeated collapses of the crater floor. Some collapse events may have produced pressure transients in the subjacent dike which acted as seismic wave sources for LP events.

Mount Etna, Sicily↗

High-resolution aeromagnetic mapping of volcanic terrain, Yellowstone National Park

High-resolution aeromagnetic data acquired over Yellowstone National Park (YNP) show contrasting patterns reflecting differences in rock composition, types and degree of alteration, and crustal structures that mirror the variable geology of the Yellowstone Plateau. The older, Eocene, Absaroka Volcanic Supergroup, a series of mostly altered, andesitic volcanic and volcaniclastic rocks partially exposed in mountains on the eastern margin of YNP, produces high-amplitude, positive magnetic anomalies, strongly contrasting with the less magnetic, younger, latest Cenozoic, Yellowstone Plateau Group, primarily a series of fresh and variably altered rhyolitic rocks covering most of YNP. The Yellowstone caldera is the centerpiece of the Yellowstone Plateau; part of its boundary can be identified on the aeromagnetic map as a series of discontinuous, negative magnetic anomalies that reflect faults or zones along which extensive hydrothermal alteration is localized. The large-volume rhyolitic ignimbrite deposits of the 0.63-Ma Lava Creek Tuff and the 2.1-Ma Huckleberry Ridge Tuff, which are prominent lithologies peripheral to the Yellowstone caldera, produce insignificant magnetic signatures. A zone of moderate amplitude positive anomalies coincides with the mapped extent of several post-caldera rhyolitic lavas. Linear magnetic anomalies reflect the rectilinear fault systems characteristic of resurgent domes in the center of the caldera. Peripheral to the caldera, the high-resolution aeromagnetic map clearly delineates flow unit boundaries of pre- and post-caldera basalt flows, which occur stratigraphically below the post-caldera rhyolitic lavas and are not exposed extensively at the surface. All of the hot spring and geyser basins, such as Norris, Upper and Lower Geyser Basins, West Thumb, and Gibbon, are associated with negative magnetic anomalies, reflecting hydrothermal alteration that has destroyed the magnetic susceptibility of minerals in the volcanic rocks. Within Yellowstone Lake, which is mostly within the Yellowstone caldera, aeromagnetic lows also are associated with known hydrothermal activity in the lake. Many of the magnetic lows extend beyond the areas of alteration and hot springs, suggesting a more extensive currently active or fossil hydrothermal system than is currently mapped. Steep magnetic gradients, suggesting faults or fractures, bound the magnetic lows. This implies that fractures localize the hot springs. Magnetic gradient trends reflect the mapped Basin and Range structural trends of north and northwest, as well as northeasterly trends that parallel the regional trend of the Snake River Plain and the track of the Yellowstone hot spot which follow the Precambrian structural grain. These trends are found both at small scales such as in hydrothermal basins and at more regional fault scales, which suggests that the regional stress field and reactivated older structures may exert some control on localization of hydrothermal activity. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Gas geochemistry of the Valles caldera region, New Mexico and comparisons with gases at Yellowstone, Long Valley and other geothermal systems

Noncondensible gases from hot springs, fumaroles, and deep wells within the Valles caldera geothermal system (210-300??C) consist of roughly 98.5 mo1% CO2, 0.5 mol% H2S, and 1 mol% other components. 3He/4He ratios indicate a deep magmatic source (R/Ra up to 6) whereas ??13C-CO2 values (-3 to -5???) do not discriminate between a mantle/magmatic source and a source from subjacent, hydrothermally altered Paleozoic carbonate rocks. Regional gases from sites within a 50-km radius beyond Valles caldera are relatively enriched in CO2 and He, but depleted in H2S compared to Valles gases. Regional gases have R/Ra values ???1.2 due to more interaction with the crust and/or less contribution from the mantle. Carbon sources for regional CO2 are varied. During 1982-1998, repeat analyses of gases from intracaldera sites at Sulphur Springs showed relatively constant CH4, H2, and H2S contents. The only exception was gas from Footbath Spring (1987-1993), which experienced increases in these three components during drilling and testing of scientific wells VC-2a and VC-2b. Present-day Valles gases contain substantially less N2 than fluid inclusion gases trapped in deep, early-stage, post-caldera vein minerals. This suggests that the long-lived Valles hydrothermal system (ca. 1 Myr) has depleted subsurface Paleozoic sedimentary rocks of nitrogen. When compared with gases from many other geothermal systems, Valles caldera gases are relatively enriched in He but depleted in CH4, N2 and Ar. In this respect, Valles gases resemble end-member hydrothermal and magmatic gases discharged at hot spots (Galapagos, Kilauea, and Yellowstone). Published by Elsevier Science B.V.

Journal of Volcanology and Geothermal Research↗

Experimental and textural constraints on mafic enclave formation in volcanic rocks

We have used experiments and textural analysis to investigate the process of enclave formation during magma mixing at Southwest Trident volcano, Alaska. Andesite enclaves are present throughout the four dacite lava flows produced by the eruption, and resemble mafic enclaves commonly found in other volcanic rocks. Our experiments replicate the pressure–temperature–time path taken by enclave-forming andesite magma as it is engulfed in dacite during magma mixing. Pressure and temperature information for the andesite and dacite are from [Coombs et al., Contrib. Mineral. Petrol. 140 (2000) 99–118]. The andesite was annealed at 1000°C, and then cooled to 890°C at rates of 110°C h −1 , 10°C h −1 , and 2°C h −1 . Once cooled to 890°C, andesite was held at this lower temperature from times ranging from 1 to 40 h. The andesite that was cooled at the slower rates of 2°C h −1 and 10°C h −1 most resembles enclave groundmass texturally and compositionally. Based on simple thermal calculations, these rates are more consistent with cooling of the andesite groundmass below an andesite–dacite interface than with cooling of enclave-sized spheres. If enclaves do crystallize as spheres, post-crystallization disaggregation must occur. Calculations using the MELTS algorithm [Ghiorso and Sack, Contrib. Mineral. Petrol. 119 (1995) 197–212] show that for incoming andesite to become less dense than the dacite ∼34 volume % of its groundmass must crystallize to undergo ∼18 volume % vesiculation; these values are similar to those determined for Southwest Trident enclaves. Thus such crystallization may lead to ‘flotation’ of enclaves and be a viable mechanism for enclave formation and dispersal. The residual melt in the cooling experiments did not evolve to rhyolitic compositions such as seen in natural enclaves due to a lack of a decompression step in the experiments. Decompression experiments on Southwest Trident dacite suggest an average ascent rate for the eruption of ∼2–3 MPa h −1 . An andesite experiment that was cooled and then decompressed at this rate contains melt that matches that of the natural enclaves. It is apparent that decompression (ascent)-induced crystallization occurs in enclaves, but not in the form of microlites as happens in the dacite host, due either to insufficient residence time at chamber temperatures or to the pre-existing microphenocrysts which act as sites for new growth.

Alaska↗

Origin and emplacement of the andesite of Burroughs Mountain, a zoned, large-volume lava flow at Mount Rainier, Washington, USA

Burroughs Mountain, situated at the northeast foot of Mount Rainier, WA, exposes a large-volume (3.4 km 3 ) andesitic lava flow, up to 350 m thick and extending 11 km in length. Two sampling traverses from flow base to eroded top, over vertical sections of 245 and 300 m, show that the flow consists of a felsic lower unit (100 m thick) overlain sharply by a more mafic upper unit. The mafic upper unit is chemically zoned, becoming slightly more evolved upward; the lower unit is heterogeneous and unzoned. The lower unit is also more phenocryst-rich and locally contains inclusions of quenched basaltic andesite magma that are absent from the upper unit. Widespread, vuggy, gabbronorite-to-diorite inclusions may be fragments of shallow cumulates, exhumed from the Mount Rainier magmatic system. Chemically heterogeneous block-and-ash-flow deposits that conformably underlie the lava flow were the earliest products of the eruptive episode. The felsic–mafic–felsic progression in lava composition resulted from partial evacuation of a vertically-zoned magma reservoir, in which either (1) average depth of withdrawal increased, then decreased, during eruption, perhaps due to variations in effusion rate, or (2) magmatic recharge stimulated ascent of a plume that brought less evolved magma to shallow levels at an intermediate stage of the eruption. Pre-eruptive zonation resulted from combined crystallization–differentiation and intrusion(s) of less evolved magma into the partly crystallized resident magma body. The zoned lava flow at Burroughs Mountain shows that, at times, Mount Rainier’s magmatic system has developed relatively large, shallow reservoirs that, despite complex recharge events, were capable of developing a felsic-upward compositional zonation similar to that inferred from large ash-flow sheets and other zoned lava flows.

Washington↗

Rejuvenated-stage volcanism after 0.6-m.y. quiescence at West Maui Volcano, Hawaii: New evidence from K-Ar ages and chemistry of Lahaina Volcanics

West Maui's rejuvenated-stage Lahaina Volcanics were erupted from four discrete sites. New KAr ages indicate two pulses of volcanism, the older about 0.6 Ma and the younger about 0.4 Ma. Compositionally the lava flows are entirely basanitic, but each pulse is diverse. The underlying postshield-stage Honolua Volcanics were emplaced by about 1.2 Ma on the basis of previously published ages. Therefore the duration of volcanic quiescence prior to rejuvenation is about 0.6 m.y. at West Maui, much longer than estimated previously. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Excess nitrogen in selected thermal and mineral springs of the Cascade Range in northern California, Oregon, and Washington: Sedimentary or volcanic in origin?

Anomalous N 2 /Ar values occur in many thermal springs and mineral springs, some volcanic fumaroles, and at least one acid-sulfate spring of the Cascade Range. Our data show that N 2 /Ar values are as high as 300 in gas from some of the hot springs, as high as 1650 in gas from some of the mineral springs, and as high as 2400 in gas from the acid-sulfate spring on Mt. Shasta. In contrast, gas discharging from hot springs that contain nitrogen and argon solely of atmospheric origin typically exhibits N 2 /Ar values of 40–80, depending on the spring temperature. If the excess nitrogen in the thermal and mineral springs is of sedimentary origin then the geothermal potential of the area must be small, but if the nitrogen is of volcanic origin then the geothermal potential must be very large. End-member excess nitrogen (δ 15 N) is +5.3‰ for the thermal waters of the Oregon Cascades but is only about +1‰ for fumaroles on Mt. Hood and the acid-sulfate spring on Mt. Shasta. Dissolved nitrogen concentrations are highest for thermal springs associated with aquifers between 120 and 140°C. Chloride is the major anion in most of the nitrogen-rich springs of the Cascade Range, and N 2 /Ar values generally increase as chloride concentrations increase. Chloride and excess nitrogen in the thermal waters of the Oregon Cascades probably originate in an early Tertiary marine formation that has been buried by the late Tertiary and Quaternary lava flows of the High Cascades. The widespread distribution of excess nitrogen that has been generated in low to moderate-temperature sedimentary environments is further proof of the restricted geothermal potential of the Cascade Range.

Cascade Range↗

Slightly thermal springs and non-thermal springs at Mount Shasta, California: Chemistry and recharge elevations

Temperature measurements, isotopic contents, and dissolved constituents are presented for springs at Mount Shasta to understand slightly thermal springs in the Shasta Valley based on the characteristics of non-thermal springs. Non-thermal springs on Mount Shasta are generally cooler than mean annual air temperatures for their elevation. The specific conductance of non-thermal springs increases linearly with discharge temperature. Springs at higher and intermediate elevations on Mount Shasta have fairly limited circulation paths, whereas low-elevation springs have longer paths because of their higher-elevation recharge. Springs in the Shasta Valley are warmer than air temperatures for their elevation and contain significant amounts of chloride and sulfate, constituents often associated with volcanic hydrothermal systems. Data for the Shasta Valley springs generally define mixing trends for dissolved constituents and temperature. The isotopic composition of the Shasta Valley springs indicates that water fell as precipitation at a higher elevation than any of the non-thermal springs. It is possible that the Shasta Valley springs include a component of the outflow from a proposed 210??C hydrothermal system that boils to supply steam for the summit acid-sulfate spring. In order to categorize springs such as those in the Shasta Valley, we introduce the term slightly thermal springs for springs that do not meet the numerical criterion of 10??C above air temperature for thermal springs but have temperatures greater than non-thermal springs in the area and usually also have dissolved constituents normally found in thermal waters. ?? 2002 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Source mechanism of long-period events at Kusatsu-Shirane Volcano, Japan, inferred from waveform inversion of the effective excitation functions

We investigate the source mechanism of long-period (LP) events observed at Kusatsu-Shirane Volcano, Japan, based on waveform inversions of their effective excitation functions. The effective excitation function, which represents the apparent excitation observed at individual receivers, is estimated by applying an autoregressive filter to the LP waveform. Assuming a point source, we apply this method to seven LP events the waveforms of which are characterized by simple decaying and nearly monochromatic oscillations with frequency in the range 1-3 Hz. The results of the waveform inversions show dominant volumetric change components accompanied by single force components, common to all the events analyzed, and suggesting a repeated activation of a sub-horizontal crack located 300 m beneath the summit crater lakes. Based on these results, we propose a model of the source process of LP seismicity, in which a gradual buildup of steam pressure in a hydrothermal crack in response to magmatic heat causes repeated discharges of steam from the crack. The rapid discharge of fluid causes the collapse of the fluid-filled crack and excites acoustic oscillations of the crack, which produce the characteristic waveforms observed in the LP events. The presence of a single force synchronous with the collapse of the crack is interpreted as the release of gravitational energy that occurs as the slug of steam ejected from the crack ascends toward the surface and is replaced by cooler water flowing downward in a fluid-filled conduit linking the crack and the base of the crater lake. ?? 2003 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Exploration and discovery in Yellowstone Lake: Results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (∼1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.

Wyoming↗

Morphology and composition of spinel in Pu'u 'O'o lava (1996-1998), Kilauea volcano, Hawaii

The morphology and composition of spinel in rapidly quenched Pu'u 'O'o vent and lava tube samples are described. These samples contain glass, olivine phenocrysts (3-5 vol.%) and microphenocrysts of spinel (~0.05 vol.%). The spinel surrounded by glass occurs as idiomorphic octahedra 5-50 μm in diameter and as chains of octahedra that are oriented with respect to each other. Spinel enclosed by olivine phenocrysts is sometimes rounded and does not generally form chains. The temperature before quenching was calculated from the MgO content of the glass and ranges from 1150 o C to 1180 o C. The oxygen fugacity before quenching was calculated by two independent methods and the log f O2 ranged from -9.2 to -9.9 (delta QFM=-1). The spinel in the Pu'u'O'o samples has a narrow range in composition with Cr/(Cr+Al)=0.61 to 0.73 and Fe 2+ /(Fe 2+ +Mg) =0.46 to 0.56. The lower the calculated temperature for the samples, the higher the average Fe 2+ /(Fe 2+ +Mg), Fe 3+ and Ti in the spinel. Most zoned spinel crystals decrease in Cr/(Cr+Al) from core to rim and, in the chains, the Cr/(Cr+Al) is greater in the core of larger crystals than in the core of smaller crystals. The occurrence of chains and hopper crystals and the presence of Cr/(Cr+Al) zoning from core to rim of the spinel suggest diffusion-controlled growth of the crystals. Some of the spinel crystals may have grown rapidly under the turbulent conditions of the summit reservoir and in the flowing lava, and the crystals may have remained in suspension for a considerable period. The rapid growth may have caused very local (μm) gradients of Cr in the melt ahead of the spinel crystal faces. The crystals seem to have retained the Cr/(Cr+Al) ratio that developed during the original growth of the crystal, but the Fe 2+ /(Fe 2+ +Mg) ratio may have equilibrated fairly rapidly with the changing melt composition due to olivine crystallization. Six of the samples were collected on the same day at various locations along a 10-km lava tube and the calculated pre-collection temperatures of the samples show a 5 o C drop with distance from the vent. The average Fe2+/(Fe2++Mg) of the spinel in these samples shows a weak positive correlation with decreasing MgO in the glass of these samples. The range in Cr 2 O 3 (0.041-0.045 wt.%) of the glass for these six samples is too small to distinguish a consistent change along the lava tube. The spinel in the Pu'u 'O'o samples shows a zoning trend in a Cr-Al-Fe 3+ diagram almost directly away from the Cr apex. This compares with a zoning trend in rapidly quenched MORB samples away from Cr coupled with decreasing Fe 3+ . The trend away from Cr displayed by spinel in rapidly quenched samples is in marked contrast to the trend of increasing Fe 3+ shown by spinel in slowly cooled lava.

Journal of Volcanology and Geothermal Research↗

Relations between seismicity and deformation during unrest in Long Valley Caldera, California, from 1995 through 1999

Unrest in Long Valley Caldera and the adjacent Sierra Nevada from 1995 through 2000 was dominated by three major episodes: (1) the March-April 1996 earthquake swarm in the east lobe of the south moat; (2) the July 1997-January 1998 caldera-wide unrest; and (3) a sequence of three M>5 earthquakes (9 June 1998, 13 July 1998, and 15 May 1999 UT) located in the Sierra Nevada block immediately south of the caldera. These three unrest episodes each had distinct characteristics with distinct implications for associated hazards. Seismicity developed as earthquake swarms for the 1996 and 1997-98 episodes, both of which were within the caldera. In contrast, the series of three M>5 earthquakes south of the caldera in 1998-99 each developed as a mainshock-aftershock sequence. Marginal deformation within the caldera associated with the 1996 swarm and the 1998-99 M>5 earthquakes is consistent with the cumulative seismic moments for the respective sequences. Deformation associated with the 1997-98 episode, however, was roughly five times larger than can be accounted for by the cumulative seismic moment of the associated earthquake swarm. We conclude that the 1997-98 episode was associated with mass transport (local intrusion of magma or magmatic brine) and that the associated earthquake swarm activity, which had a relatively high b -value of 1.2, was largely driven by the intrusive process. In contrast, the 1996 earthquake swarm and the 1998-99 M>5 mainshock-aftershock sequences, both with 'normal' b -values of ???0.9, represent brittle relaxation to previously accumulated stresses associated with little or no mass transport. These relations emphasize the importance of simultaneous, real-time monitoring of both seismicity and deformation as a basis for judging whether an evolving unrest episode has the potential for culminating in a volcanic eruption. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

The mechanics of unrest at Long Valley caldera, California: 1. Modeling the geometry of the source using GPS, leveling and two-color EDM data

We surveyed 44 existing leveling monuments in Long Valley caldera in July 1999, using dual frequency global positioning system (GPS) receivers. We have been able to tie GPS and leveling to a common reference frame in the Long Valley area and computed the vertical deformation by differencing GPS-based and leveled orthometric heights. The resurgent dome uplifted 74??7 cm from 1975 to 1999. To define the inflation source, we invert two-color EDM and uplift data from the 1985-1999 unrest period using spherical or ellipsoidal sources. We find that the ellipsoidal source satisfies both the vertical and horizontal deformation data, whereas the spherical point source cannot. According to our analysis of the 1985-1999 data, the main source of deformation is a prolate ellipsoid located beneath the resurgent dome at a depth of 5.9 km (95% bounds of 4.9-7.5 km). This body is vertically elongated, has an aspect ratio of 0.475 (95% bounds are 0.25-0.65) and a volume change of 0.086 km3 (95% bounds are 0.06-0.13 km3). Failure to account for the ellipsoidal nature of the source biases the estimated source depth by 2.1 km (35%), and the source volume by 0.038 km3 (44%). ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

The mechanics of unrest at Long Valley caldera, California. 2. Constraining the nature of the source using geodetic and micro-gravity data

We model the source of inflation of Long Valley caldera by combining geodetic and micro-gravity data. Uplift from GPS and leveling, two-color EDM measurements, and residual gravity change determinations are used to estimate the intrusion geometry, assuming a vertical prolate ellipsoidal source. The U.S. Geological Survey occupied the Long Valley gravity network six times from 1980 to 1985. We reoccupied this network twice, in the summer of 1998 (33 stations), and the summer of 1999 (37 stations). Before gravity data can be used to estimate the density of the intrusion, they must be corrected for the effect of vertical deformation (the free-air effect) and changes in the water table. We use geostatistical techniques to interpolate uplift and water table changes at the gravity stations. The inflation source (a vertical prolate ellipsoid) is located 5.9 km beneath the resurgent dome with an aspect ratio equal to 0.475, a volume change from 1982 to 1999 of 0.136 km3 and a density of around 1700 kg/m3. A bootstrap method was employed to estimate 95% confidence bounds for the parameters of the inflation model. We obtained a range of 0.105-0.187 km3 for the volume change, and 1180-2330 kg/m3 for the density. Our results do not support hydrothermal fluid intrusion as the primary cause of unrest, and confirm the intrusion of silicic magma beneath Long Valley caldera. Failure to account for the ellipsoidal nature of the source biases the estimated source depth by 2.9 km (a 33% increase), the volume change by 0.019 km3 (a 14% increase) and the density by about 1200 kg/m3 (a 40% increase). ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Deformation of the Long Valley Caldera, California: Inferences from measurements from 1988 to 2001

Two periods of volcanic unrest occurred between 1989 and 1998 in the Long Valley Caldera, eastern California. Numerous earthquakes were recorded, and these periods of unrest were documented with high-precision geodetic measurements. The first round of unrest started rapidly in late 1989 and slowly decreased in rate through the early 1990s. For this interval there are both leveling and two-color electronic distance meter (EDM) measurements. The second round of unrest started slowly in mid-1997, climaxed in late 1997, and rapidly returned to quiescence by mid-1998. Deformation was recorded by both the two-color EDM and continuous GPS. Both episodes require inflation at 6-7 km beneath the resurgent dome, and both episodes had roughly 0.1 m extension across the resurgent dome. In addition, the data presented here suggest that there is a deeper, 10-20 km, inflation source beneath the south moat of the caldera. For both episodes, the better-resolved inflation beneath the resurgent dome is a near-vertical, prolate spheroid rather than an isotropic source, which suggests that magma came up through vertical cracks. However, the modeling suggests that the location changed with the depth from 6.0 to 6.7 km for the later episode. In contrast to the earlier episode, the 1997-1998 episode has additional deformation in the south moat, where the simplest model is that of a right-lateral slip on a steeply dipping plane that is defined by the location of earthquakes in the south moat. Models of the time-dependent behavior suggest that slip on this fault occurred from late November through December 1997, corresponding to the time of greatest moment release by the earthquake swarm in the south moat. Confounding the interpretation of these data is an active geothermal field near the center of the EDM network and adjacent to the south moat and resurgent dome. Additional modeling of leveling and EDM data within the geothermal field during a period of low rate of inflation of the dome suggests some methods of adjustments to the EDM data during the inflation episodes. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Water-level changes induced by local and distant earthquakes at Long Valley caldera, California

Distant as well as local earthquakes have induced groundwater-level changes persisting for days to weeks at Long Valley caldera, California. Four wells open to formations as deep as 300 m have responded to 16 earthquakes, and responses to two earthquakes in the 3-km-deep Long Valley Exploratory Well (LVEW) show that these changes are not limited to weathered or unconsolidated near-surface rocks. All five wells exhibit water-level variations in response to earth tides, indicating they can be used as low-resolution strainmeters. Earthquakes induce gradual water-level changes that increase in amplitude for as long as 30 days, then return more slowly to pre-earthquake levels. The gradual water-level changes are always drops at wells LKT, LVEW, and CH-10B, and always rises at well CW-3. At a dilatometer just outside the caldera, earthquake-induced strain responses consist of either a step followed by a contractional strain-rate increase, or a transient contractional signal that reaches a maximum in about seven days and then returns toward the pre-earthquake value. The sizes of the gradual water-level changes generally increase with earthquake magnitude and decrease with hypocentral distance. Local earthquakes in Long Valley produce coseismic water-level steps; otherwise the responses to local earthquakes and distant earthquakes are indistinguishable. In particular, water-level and strain changes in Long Valley following the 1992 M7.3 Landers earthquake, 450 km distant, closely resemble those initiated by a M4.9 local earthquake on November 22, 1997, during a seismic swarm with features indicative of fluid involvement. At the LKT well, many of the response time histories are identical for 20 days after each earthquake, and can be matched by a theoretical solution giving the pore pressure as a function of time due to diffusion of a nearby, instantaneous, pressure drop. Such pressure drops could be produced by accelerated inflation of the resurgent dome by amounts too small to be detected by the two-color electronic distance-measuring network. Opening-mode displacement in the south moat, inferred to have followed a M4.9 earthquake on November 22, 1997, could also create extensional strain on the dome and lead to water-level changes similar to those following dome inflation. Contractional strain that could account for earthquake-induced water-level rises at the CW-3 well is inconsistent with geodetic observations. We instead attribute these water-level rises to diffusion of elevated fluid pressure localized in the south moat thermal aquifer. For hydraulic diffusivities appropriate to the upper few hundred meters at Long Valley, an influx of material at temperatures of 300°C can thermally generate pressure of 6 m of water or more, an order of magnitude larger than needed to account for the CW-3 water-level rises. If magma or hot aqueous fluid rises to within 1 km of the surface in the eastern part of the south moat, then hydraulic diffusivities are high enough to allow fluid pressure to propagate to CW-3 on the time scale observed. The data indicate that seismic waves from large distant earthquakes can stimulate upward movement of fluid in the hydrothermal system at Long Valley.

California↗