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 343 records · Page 19Linked to original sources

Deformation and seismic precursors to dome-collapse and fountain-collapse nuées ardentes at Merapi Volcano, Java, Indonesia, 1994-1998

Following the eruption of January 1992, episodes of lava dome growth accompanied by generation of dome-collapse nuées ardentes occurred in 1994–1998. In addition, nuées ardentes were generated by fountain-collapse in January 1997, and the 1998 events also suggest an explosive component. Significant tilt and seismic precursors on varying time scales preceded these events. Deformation about the summit has been detected by electronic tiltmeters since November 1992, with inflation corresponding generally to lava dome growth, and deflation (or decreased inflation) corresponding to loss of dome mass. Strong short-term (days to weeks) accelerations in tilt rate and seismicity occurred prior to the major nuées ardentes episodes, apart from those of 22 November 1994 which were preceded by steadily increasing tilt for over 200 days but lacked short-term precursors. Because of the combination of populated hazardous areas and the lack of an issued warning, about 100 casualties occurred in 1994. In contrast, the strong precursors in 1997 and 1998 provided advance warning to observatory scientists, enabled the stepped raising of alert levels, and aided hazard management. As a result of these factors, but also the fortunate fact that the large nuées ardentes did not quite descend into populated areas, no casualties occurred. The nuée ardente episode of 1994 is interpreted as purely due to gravitational collapse, whereas those of 1997 and 1998 were influenced by gas-pressurization of the lava dome.

Journal of Volcanology and Geothermal Research↗

Nuées ardentes of 22 November 1994 at Merapi volcano, Java, Indonesia

Nuées ardentes associated with dome collapse on 22 November 1994, at Merapi volcano traveled to the south–southwest as far as 6.5 km, and collectively accumulated roughly 2.5–3 million cubic meters of deposits. The damaged area comprises 9.5 km 2 and is covered by two nuée ardente facies, a conventional “Merapi-type”, valley-fill block-and-ash flow facies and a pyroclastic surge facies. The proximal deposits reflect the accumulation of dozens of nuées ardentes, with many subsidiary flow units. The distal deposits are more simply organized, as only a few individual events reached to distances >3.5 km. The stratigraphic relationships north of Turgo hill indicate that the surge deposits are a facies of particularly mobile nuées ardentes that also deposited channeled block-and-ash flow facies. They further suggest that the surge facies beyond the channel margins correlate laterally with a finer-grained sublayer locally developed at the base of the block-and-ash flow facies. Eyewitness reports suggest that the emplacement of the block-and-ash flow facies in the distal part of the Boyong river may have followed, by a short time interval, the destruction and deposition of the surge facies at Turgo village. The stratigraphy is in accord with the eyewitness reports. The surge facies was emplaced by a dilute surge current, detached from the same dome-collapse nuée ardente that, as a separate flow unit, subsequently emplaced the distal block-and-ash deposit in the Boyong valley. The detachment occurred at higher elevations, likely at or above the slope break at about 2000 m elevation. This flow separation enabled the surge current to shortcut over the landscape and to emplace its deposit even as the block-and-ash flow continued its tortuous southward movement in the Boyong channel. Dome-collapse nuée ardente activity formed the bulk of the eruption, which was accompanied by virtually no significant vertical summit explosive activity.

Journal of Volcanology and Geothermal Research↗

Sequential dome-collapse nuées ardentes analyzed from broadband seismic data, Merapi Volcano, Indonesia

During the sequential dome collapse of Merapi Volcano on 22 November 1994, a broadband seismic station on the western slope was the only operational seismic equipment that provided continuous on-scale recording of the event. According to visual and seismic observations, the collapse activity lasted about 10 h. We divide the activity into two phases: (I) a period with fluctuating but generally increasing seismic activity associated with 24 moderate to large dome-collapse nuées ardentes, lasting about 40 min and culminating with the largest event at 10:54 (all times reported as local time); and (II) activity from 11:42 until 20:00 described by scattered clusters of individual rockfalls and 20 nuées ardentes. The broadband data were evaluated using the assumption that avalanches with the same source areas and descent paths exhibit a linear relation between source volume and recorded seismic-amplitude envelope area. A result of the analysis is the determination of the volume of selected individual events. From the field surveys, the total volume of the collapsed dome lava is 2.6 Mm 3 . Discounting the volumetric influence of rockfalls, the average size of the 44 nuées ardentes is therefore about 60,000 m 3 . The largest collapse event at 10:54 is estimated to involve 260,000 m 3 , based on an analysis of the seismicity. The remaining 23 phase I events averaged 60,000 m 3 , with the total volume of all phase I events accounting for 63% of the unstable dome. The 20 phase II events comprised 37% of the total volume and averaged 47,000 m 3 . The methods described here can be put to practical use in real-time monitoring situations. Broadband data were essential in this study primarily because of the wide dynamic range.

Journal of Volcanology and Geothermal Research↗

Emplacement temperatures of the November 22, 1994 nuee ardente deposits, Merapi Volcano, Java

A study of emplacement temperatures was carried out for the largest of the 22 November 1994 nu&eacute;e ardente deposits at Merapi Volcano, based mainly on the response of plastic and woody materials subjected to the hot pyroclastic current and the deposits, and to some extent on eyewitness observations. The study emphasizes the Turgo&ndash;Kaliurang area in the distal part of the area affected by the nu&eacute;e ardente, where nearly 100 casualties occurred. The term nu&eacute;e ardente as used here includes channeled block-and-ash flows, and associated ash-clouds of surge and fallout origins. The emplacement temperature of the 8 m thick channeled block-and-ash deposit was relatively high, &sim;550&deg;C, based mainly on eyewitness reports of visual thermal radiance. Emplacement temperatures for ash-cloud deposits a few cm thick were deduced from polymer objects collected at Turgo and Kaliurang. Most polymers do not display a sharp melting range, but polyethylene terephthalate used in water bottles melts between 245 and 265&deg;C, and parts of the bottles that had been deformed during fabrication molding turn a milky color at 200&deg;C. The experimental evidence suggests that deposits in the Turgo area briefly achieved a maximum temperature near 300&deg;C, whereas those near Kaliurang were <200&deg;C. Maximum ash deposit temperatures occurred in fallout with a local source in the channeled block-and-ash flow of the Boyong river valley; the surge deposit was cooler (&sim;180&deg;C) due to entrainment of cool air and soils, and tree singe-zone temperatures were around 100&deg;C.

Journal of Volcanology and Geothermal Research↗

Instrumental lahar monitoring at Merapi Volcano, Central Java, Indonesia

More than 50 volcanic debris flows or lahars were generated around Mt Merapi during the first rainy season following the nuees ardentes of 22 November 1994. The rainfalls that triggered the lahars were analyzed, using such instruments as weather radar and telemetered rain gauges. Lahar dynamics were also monitored, using new non-contact detection instrumentation installed on the slopes of the volcano. These devices include real-time seismic amplitude measurement (RSAM), seismic spectral amplitude measurement (SSAM) and acoustic flow monitoring (AFM) systems. Calibration of the various systems was accomplished by field measurements of flow velocities and discharge, contemporaneously with instrumental monitoring. The 1994–1995 lahars were relatively short events, their duration in the Boyong river commonly ranging between 30 min and 1 h 30 min. The great majority (90%) of the lahars was recognized at Kaliurang village between 13:00 and 17:30 h, due to the predominance of afternoon rainfalls. The observed mean velocity of lahar fronts ranged between 1.1 and 3.4 m/s, whereas the peak velocity of the flows varied from 11 to 15 m/s, under the Gardu Pandang viewpoint location at Kaliurang, to 8–10 m/s at a section 500 m downstream from this site. River slopes vary from 28 to 22 m/km at the two sites. Peak discharges recorded in various events ranged from 33 to 360 m 3 /s, with the maximum value of peak discharge 360 m 3 /s, on 20 May 1995. To improve the lahar warning system along Boyong river, some instrumental thresholds were proposed: large and potentially hazardous lahars may be detected by RSAM units exceeding 400, SSAM units exceeding 80 on the highest frequency band, or AFM values greater than 1500 mV on the low-gain, broad-band setting.

Journal of Volcanology and Geothermal Research↗

Period doubling and other nonlinear phenomena in volcanic earthquakes and tremor

Evidence of subharmonic period-doubling cascades has recently been recognized in seismograms of volcanic tremor from several volcanoes. This phenomenon occurs only in nonlinear systems, and is the commonest route by which such systems change from periodic to chaotic behavior. It is predicted to occur in a model of volcanic tremor excitation by flow-induced vibration, and it might well also occur in other volcano-seismic source process. If the possibility of period doubling is not taken into account in interpreting spectra of tremor and long-period earthquakes, then low-frequency "sub-harmonic" oscillations may be mis-identified as normal modes of a linear acoustic resonator, leading to errors of an order of magnitude or more in inferred magma-body dimensions. This example illustrates the importance of nonlinear phenomena in attempts to understand volcano-seismic phenomena physically. Linear systems are fundamentally incapable of causing earthquakes or exciting tremor, so nonlinearity is essential to any theory of volcano-seismic phenomena. Nonlinear processes are in many respects qualitatively different from linear ones. A few of their characteristics that might be relevant in volcanoes include the possibility: (1) that damping might increase, rather than decrease, oscillation frequencies; and (2) that these frequencies might be functions of the amplitude of oscillation, so that temporal variations in spectral peak frequencies might not be manifestations of changes of conditions within the magmatic system.

Journal of Volcanology and Geothermal Research↗

Spherical and ellipsoidal volcanic sources at Long Valley caldera, California, using a genetic algorithm inversion technique

We model the second inflation period at Long Valley caldera, California using a genetic algorithm technique and high quality geodetic measurements of elevation changes and baseline extensions. We compare two source inversions for both spherical Mogi point sources and the finite prolate ellipsoid of Yang and Davis. A sensitivity analysis for the genetic algorithm is performed based upon synthetic data set inversions on similar sources in order to better constrain the areal location, orientation, and volume of the potential sources. The spherical sources are well constrained, the larger located at 9.9 km beneath the resurgent dome, with a volume of 0.036 km 3 , while the second, at only 0.008 km 3 , is located at a depth of 7.3 km beneath the south moat. The depths to the ellipsoidal sources are switched, with the larger source at a depth of 9.6 km and the smaller at 11.8 km, with volumes of 0.037 and 0.002 km 3 , respectively.

California↗

Temporal and spectral characteristics of seismicity observed at Popocatepetl volcano, central Mexico

Popocatepetl volcano entered an eruptive phase from December 21, 1994 to March 30, 1995, which was characterized by ash and fumarolic emissions. During this eruptive episode, the observed seismicity consisted of volcano-tectonic (VT) events, long-period (LP) events and sustained tremor. Before the initial eruption on December 21, VT seismicity exhibited no increase in number until a swarm of VT earthquakes was observed at 01:31 hours local time. Visual observations of the eruption occurred at dawn the next morning. LP activity increased from an average of 7 events a day in October 1994 to 22 events per day in December 1994. At the onset of the eruption, LP activity peaked at 49 events per day. LP activity declined until mid-January 1995 when no events were observed. Tremor was first observed about one day after the initial eruption and averaged 10 h per episode. By late February 1995, tremor episodes became more intermittent, lasting less than 5 min, and the number of LP events returned to pre-eruption levels (7 events per day). Using a spectral ratio technique, low-frequency oceanic microseismic noise with a predominant peak around 7 s was removed from the broadband seismic signal of tremor and LP events. Stacks of corrected tremor episodes and LP events show that both tremor and LP events contain similar frequency features with major peaks around 1.4 Hz. Frequency analyses of LP events and tremor suggest a shallow extended source with similar radiation pattern characteristics. The distribution of VT events (between 2.5 and 10 km) also points to a shallow source of the tremor and LP events located in the first 2500 m beneath the crater. Under the assumption that the frequency characteristics of the signals are representative of an oscillator we used a fluid-filled-crack model to infer the length of the resonator.

Journal of Volcanology and Geothermal Research↗

Record of late holocene debris avalanches and lahars at Iliamna Volcano, Alaska

Iliamna Volcano is a 3053-meter high, glaciated stratovolcano in the southern Cook Inlet region of Alaska and is one of seven volcanoes in this region that have erupted multiple times during the past 10,000 yr. Prior to our studies of Iliamna Volcano, little was known about the frequency, magnitude, and character of Holocene volcanic activity. Here we present geologic evidence of the most recent eruptive activity of the volcano and provide the first outline of Late Holocene debris-avalanche and lahar formation. Iliamna has had no documented historical eruptions but our recent field investigations indicate that the volcano has erupted at least twice in the last 300 yr. Clay-rich lahar deposits dated by radiocarbon to ???1300 and ???90 yr BP are present in two major valleys that head on the volcano. These deposits indicate that at least two large, possibly deep-seated, flank failures of the volcanic edifice have occurred in the last 1300 yr. Noncohesive lahar deposits likely associated with explosive pyroclastic eruptions date to 2400-1300,>1500,???300, and <305 yr BP. Debris-avalanche deposits from recent and historical small-volume slope failures of the hydrothermally altered volcanic edifice cover most of the major glaciers on the volcano. Although these deposits consist almost entirely of hydrothermally altered rock debris and snow and ice, none of the recently generated debris avalanches evolved to lahars. A clay-rich lahar deposit that formed <90??60 radiocarbon yr BP and entered the Johnson River Valley southeast of the volcano cannot be confidently related to an eruption of Iliamna Volcano, which has had no known historical eruptions. This deposit may record an unheralded debris avalanche and lahar. ?? 2000 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Anomalously high b-values in the South Flank of Kilauea volcano, Hawaii: Evidence for the distribution of magma below Kilauea's East rift zone

The pattern of b-value of the frequency-magnitude relation, or mean magnitude, varies little in the Kaoiki-Hilea area of Hawaii, and the b-values are normal, with b = 0.8 in the top 10 km and somewhat lower values below that depth. We interpret the Kaoiki-Hilea area as relatively stable, normal Hawaiian crust. In contrast, the b-values beneath Kilauea's South Flank are anomalously high (b = 1.3-1.7) at depths between 4 and 8 km, with the highest values near the East Rift zone, but extending 5-8 km away from the rift. Also, the anomalously high b-values vary along strike, parallel to the rift zone. The highest b-values are observed near Hiiaka and Pauahi craters at the bend in the rift, the next highest are near Makaopuhi and also near Puu Kaliu. The mildest anomalies occur adjacent to the central section of the rift. The locations of the three major and two minor b-value anomalies correspond to places where shallow magma reservoirs have been proposed based on analyses of seismicity, geodetic data and differentiated lava chemistry. The existence of the magma reservoirs is also supported by magnetic anomalies, which may be areas of dike concentration, and self-potential anomalies, which are areas of thermal upwelling above a hot source. The simplest explanation of these anomalously high b-values is that they are due to the presence of active magma bodies beneath the East Rift zone at depths down to 8 km. In other volcanoes, anomalously high b-values correlate with volumes adjacent to active magma chambers. This supports a model of a magma body beneath the East Rift zone, which may widen and thin along strike, and which may reach 8 km depth and extend from Kilauea's summit to a distance of at least 40 km down rift. The anomalously high b-values at the center of the South Flank, several kilometers away from the rift, may be explained by unusually high pore pressure throughout the South Flank, or by anomalously strong heterogeneity due to extensive cracking, or by both phenomena. The major b-value anomalies are located SSE of their parent reservoirs, in the direction of motion of the flank, suggesting that magma reservoirs leave an imprint in the mobile flank. We hypothesize that the extensive cracking may have been acquired when the anomalous parts of the South Flank, now several kilometers distant from the rift zone, were generated at the rift zone near persistent reservoirs. Since their generation, these volumes may have moved seaward, away from the rift, but earthquakes occurring in them still use the preexisting complex crack distribution. Along the decollement plane at 10 km depth, the b-values are exceptionally low (b = 0.5), suggesting faulting in a more homogeneous medium. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Trends in long-period seismicity related to magmatic fluid compositions

Sound speeds and densities are calculated for three different types of fluids: gas-gas mixture; ash-gas mixture; and bubbly liquid. These fluid properties are used to calculate the impedance contrast (Z) and crack stiffness (C) in the fluid-driven crack model (Chouet: J. Geophys. Res., 91 (1986) 13,967; 101 (1988) 4375; A seismic model for the source of long-period events and harmonic tremor. In: Gasparini, P., Scarpa, R., Aki, K. (Eds.), Volcanic Seismology, IAVCEI Proceedings in Volcanology, Springer, Berlin, 3133). The fluid-driven crack model describes the far-field spectra of long-period (LP) events as modes of resonance of the crack. Results from our calculations demonstrate that ash-laden gas mixtures have fluid to solid density ratios comparable to, and fluid to solid velocity ratios lower than bubbly liquids (gas-volume fractions <10%). This difference results in synthetic far-field spectra with higher impedance contrasts and narrower spectral bandwidths for ash-laden gas mixture than spectra for bubbly liquids. Spectral characteristics are described in terms of the quality factor Q-1. Q-1 is measured by the ratio of the frequency of the dominant spectral peak to the bandwidth of the peak measured at one half of its amplitude. This factor expresses the losses of energy due to elastic radiation Q-1r and other dissipative mechanisms Q-1i at the source, Q-1 = Q-1r + Q-1i. Spectra for LP events recorded at active volcanoes such as Galeras in Colombia and Kilauea in Hawaii, have Q-1 factors in the range of 0.1-0.002. The Q-1r factors due to radiation loss calculated for a sphere filled with a H2O-CO2 or H2O-SO2 gas mixture, vary between 0.0015 and 0.0040 with a change in wt% H2O at 800-1600 K and 10-50 MPa. For gas-rich mixtures, Q-1r has a strong dependence on resonator geometry (spherical versus rectangular). The spectra from a resonating sphere filled with gas-rich mixture yields values of Q-1r an order of magnitude smaller than those from a rectangular crack. For a resonating crack filled with an ash-gas mixture (or pseudogas), Q-1r varies parabolically from ???0.006 for an ash-rich mixture, to 0.0015 or 0.0023 for a H2O-rich or CO2-rich mixture at 800 K and 25 MPa. For low (<20%) gas-volume fraction fluids (foams, bubbly fluids and ash-rich pseudogases), the magnitudes for Q-1r are independent of crack geometry. Spectra associated with a foam (gas-volume fractions 10-90%) or bubbly basalt (gas-volume fractions <10%) may have a dominant spectral peak with values of Q-1r on the order of 0.01 and 0.1, respectively. The spectra from a resonating sphere filled with a foam containing >20% gas-volume fraction yields values of Q-1r similar to those for a rectangular crack. As with gas-gas and ash-gas mixtures, an increase in mass fraction narrows the bandwidth of the dominant mode and shifts the spectra to lower frequencies. Including energy losses due to dissipative processes in a bubbly liquid increases attenuation. Attenuation may also be higher in ash-gas mixtures and foams if the effects of momentum and mass transfer between the phases were considered in the calculations. ?? 2001 Elsevier Science B. V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Implications for eruptive processes as indicated by sulfur dioxide emissions from Kilauea Volcano, Hawai'i, 1979-1997

Kı̄lauea Volcano, Hawai‘i, currently hosts the longest running SO 2 emission-rate data set on the planet, starting with initial surveys done in 1975 by Stoiber and his colleagues. The 17.5-year record of summit emissions, starting in 1979, shows the effects of summit and east rift eruptive processes, which define seven distinctly different periods of SO 2 release. Summit emissions jumped nearly 40% with the onset (3 January 1983) of the Pu`u `Ō`ō-Kūpaianaha eruption on the east rift zone (ERZ). Summit SO 2 emissions from Kı̄lauea showed a strong positive correlation with short-period, shallow, caldera events, rather than with long-period seismicity as in more silicious systems. This correlation suggests a maturation process in the summit magma-transport system from 1986 through 1993. During a steady-state throughput-equilibrium interval of the summit magma reservoir, integration of summit-caldera and ERZ SO 2 emissions reveals an undegassed volume rate of effusion of 2.1×10 5 m 3 /d. This value corroborates the volume-rate determined by geophysical methods, demonstrating that, for Kı̄lauea, SO 2 emission rates can be used to monitor effusion rate, supporting and supplementing other, more established geophysical methods. For the 17.5 years of continuous emission rate records at Kı̄lauea, the volcano has released 9.7×10 6 t (metric tonnes) of SO 2 , 1.7×10 6 t from the summit and 8.0×10 6 t from the east rift zone. On an annual basis, the average SO 2 release from Kı̄lauea is 4.6×10 5 t/y, compared to the global annual volcanic emission rate of 1.2×10 7 t/y.

Journal of Volcanology and Geothermal Research↗

Magmatic gas scrubbing: Implications for volcano monitoring

Despite the abundance of SO2(g) in magmatic gases, precursory increases in magmatic SO2(g) are not always observed prior to volcanic eruption, probably because many terrestrial volcanoes contain abundant groundwater or surface water that scrubs magmatic gases until a dry pathway to the atmosphere is established. To better understand scrubbing and its implications for volcano monitoring, we model thermochemically the reaction of magmatic gases with water. First, we inject a 915??C magmatic gas from Merapi volcano into 25??C air-saturated water (ASW) over a wide range of gas/water mass ratios from 0.0002 to 100 and at a total pressure of 0.1 MPa. Then we model closed-system cooling of the magmatic gas, magmatic gas-ASW mixing at 5.0 MPa, runs with varied temperature and composition of the ASW, a case with a wide range of magmatic-gas compositions, and a reaction of a magmatic gas-ASW mixture with rock. The modeling predicts gas and water compositions, and, in one case, alteration assemblages for a wide range of scrubbing conditions; these results can be compared directly with samples from degassing volcanoes. The modeling suggests that CO2(g) is the main species to monitor when scrubbing exists; another candidate is H2S(g), but it can be affected by reactions with aqueous ferrous iron. In contrast, scrubbing by water will prevent significant SO2(g) and most HCl(g) emissions until dry pathways are established, except for moderate HCl(g) degassing from pH < 0.5 hydrothermal waters. Furthermore, it appears that scrubbing will prevent much, if any, SO2(g) degassing from long-resident boiling hydrothermal systems. Several processes can also decrease or increase H2(g) emissions during scrubbing making H2(g) a poor choice to detect changes in magma degassing. We applied the model results to interpret field observations and emission rate data from four eruptions: (1) Crater Peak on Mount Spurr (1992) where, except for a short post-eruptive period, scrubbing appears to have drastically diminished pre-, inter-, and post-eruptive SO2(g) emissions, but had much less impact on CO2(g) emissions. (2) Mount St. Helens where scrubbing of SO2(g) was important prior to and three weeks after the 18 May 1980 eruption. Scrubbing was also active during a period of unrest in the summer of 1998. (3) Mount Pinatubo where early drying out prevented SO2(g) scrubbing before the climactic 15 June 1991 eruption. (4) The ongoing eruption at Popocate??petl in an arid region of Mexico where there is little evidence of scrubbing. In most eruptive cycles, the impact of scrubbing will be greater during pre- and post-eruptive periods than during the main eruptive and intense passive degassing stages. Therefore, we recommend monitoring the following gases: CO2(g) and H2S(g) in precursory stages; CO2(g), H2S(g), SO2(g), HCl(g), and HF(g) in eruptive and intense passive degassing stages; and CO2(g) and H2S(g) again in the declining stages. CO2(g) is clearly the main candidate for early emission rate monitoring, although significant early increases in the intensity and geographic distribution of H2S(g) emissions should be taken as an important sign of volcanic unrest and a potential precursor. Owing to the difficulty of extracting SO2(g) from hydrothermal waters, the emergence of > 100 t/d (tons per day) of SO2(g) in addition to CO2(g) and H2S(g) should be taken as a criterion of magma intrusion. Finally, the modeling suggests that the interpretation of gas-ratio data requires a case-by-case evaluation since ratio changes can often be produced by several mechanisms; nevertheless, several gas ratios may provide useful indices for monitoring the drying out of gas pathways. Published by Elsevier Science B.V.

Journal of Volcanology and Geothermal Research↗

Evolution of the December 14, 1989 precursory long-period event swarm at Redoubt volcano, Alaska

The intermittency pattern and evolution in waveforms of long-period (LP) seismic events during the intense, 23-h swarm that preceded the December 14, 1989 eruption of Redoubt volcano are investigated. Utilizing cross correlation to exploit the high degree of similarity among waveforms, a substantially more complete event catalog is generated than was available from near realtime detection based on short-term/long-term amplitude ratios, which was saturated by the high rate of activity. The temporal magnitude distribution of the predominant LP events is found to have an unusual banded structure in which the average magnitude of each band slowly increases and then decreases through time. A bifurcation that appears in the uppermost band shortly after the peak in magnitudes is characterized by a quasi-periodicity in intermittency and magnitude that is reminiscent of one of the classic routes to chaotic behavior in some non-linear systems. The waveforms of the predominant events evolve slowly but unsteadily through time. These gradual changes appear to result from variations in the relative amplitudes of spectral peaks that remain stable in frequency, which suggests that they are due to differential excitation of a single, resonant source. Two other previously unrecognized, repetitive waveforms are also identified, but the signals from these secondary events are not clearly recorded at distances beyond the closest station. Similarities among the spectra of the predominant and secondary events suggest that the signals from these events also could represent different modes of exciting the same source. Significant changes in the rates and the sizes of the largest of these secondary events appear to coincide with the peak in the size distribution of the predominant LPs. At least some of the non-repetitive LP waveforms in the swarm appear to be the result of the superposition of signals from the rapid repetition of predominant LP source, thus placing a constraint on the repeat time of the triggering mechanism for this source. A lone hybrid event, which has a waveform character intermediate between the predominant LP events and high-frequency volcano-tectonic events, was also identified in the swarm; the occurrence of this event provides important evidence that the low-frequency character of the LP events is a source rather than a path or site effect. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Wavefield properties of a shallow long-period event and tremor at Kilauea Volcano, Hawaii

The wavefields of tremor and a long-period (LP) event associated with the ongoing eruptive activity at Kilauea Volcano, Hawaii, are investigated using a combination of dense small-aperture (300 m) and sparse large-aperture (5 km) arrays deployed in the vicinity of the summit caldera. Measurements of azimuth and slowness for tremor recorded on the small-aperture array indicate a bimodal nature of the observed wavefield. At frequencies below 2 Hz, the wavefield is dominated by body waves impinging the array with steep incidence. These arrivals are attributed to the oceanic microseismic noise. In the 2-6 Hz band, the wavefield is dominated by waves propagating from sources located at shallow depths (<1 km) beneath the eastern edge of the Halemaumau pit crater. The hypocenter of the LP event, determined from frequency-slowness analyses combined with phase picks, appears to be located close to the source of tremor but at a shallower depth (<0.1 km). The wavefields of tremor and LP event are characterized by a complex composition of body and surface waves, whose propagation and polarization properties are strongly affected by topographic and structural features in the summit caldera region. Analyses of the directional properties of the wavefield in the 2-6 Hz band point to the directions of main scattering sources, which are consistent with pronounced velocity contrasts imaged in a high-resolution three-dimensional velocity model of the caldera region. The frequency and Q of the dominant peak observed in the spectra of the LP event may be explained as the dominant oscillation mode of a crack with scale length 20-100 m and aperture of a few centimeters filled with bubbly water. The mechanism driving the shallow tremor appears to be consistent with a sustained excitation originating in the oscillations of a bubbly cloud resulting from vesiculation and degassing in the magma. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Time-variation of hydrothermal discharge at selected sites in the Western United States: Implications for monitoring

We compiled time series of hydrothermal discharge consisting of 3593 chloride- or heat-flux measurements from 24 sites in the Yellowstone region, the northern Oregon Cascades, Lassen Volcanic National Park and vicinity, and Long Valley, California. At all of these sites the hydrothermal phenomena are believed to be as yet unaffected by human activity, though much of the data collection was driven by mandates to collect environmental-baseline data in acticipation of geothermal development. The time series average 19 years in length and some of the Yellowstone sites have been monitored intermittently for over 30 years. Many sites show strong seasonality but few show clear long-term trends, and at most sites statistically significant decadal-scale trends are absent. Thus, the data provide robust estimates of advective heat flow ranging from ~130 MW in the north-central Oregon Cascades to ~6100 MW in the Yellowstone region, and also document Yellowstone hydrothermal chloride and arsenic fluxes of 1740 and 15-20 g/s, respectively. The discharge time series show little sensitivity to regional tectonic events such as earthquakes or inflation/deflation cycles. Most long-term monitoring to date has focused on high-chloride springs and low-temperature fumaroles. The relative stability of these features suggests that discharge measurements done as part of volcano-monitoring programs should focus instead on high-temperature fumaroles, which may be more immediately linked to the magmatic heat source.

Journal of Volcanology and Geothermal Research↗

Predictions of hydrothermal alteration within near-ridge oceanic crust from coordinated geochemical and fluid flow models

Coordinated geochemical and hydrological calculations guide our understanding of the composition, fluid flow patterns, and thermal structure of near-ridge oceanic crust. The case study presented here illustrates geochemical and thermal changes taking place as oceanic crust ages from 0.2 to 1.0 Myr. Using a finite element code, we model fluid flow and heat transport through the upper few hundred meters of an abyssal hill created at an intermediate spreading rate. We use a reaction path model with a customized database to calculate equilibrium fluid compositions and mineral assemblages of basalt and seawater at 500 bars and temperatures ranging from 150 to 400??C. In one scenario, reaction path calculations suggest that volume increases on the order of 10% may occur within portions of the basaltic basement. If this change in volume occurred, it would be sufficient to fill all primary porosity in some locations, effectively sealing off portions of the oceanic crust. Thermal profiles resulting from fluid flow simulations indicate that volume changes along this possible reaction path occur primarily within the first 0.4 Myr of crustal aging. ?? 2001 Elsevier Science B.V. All rights reserved.

Journal of Volcanology and Geothermal Research↗

Density of basalt core from Hilo drill hole, Hawaii

Density measurements of 1600 samples of core from 889 to 3097 m depth below sea level in the Hawaii Scientific Drilling Program hole near Hilo, Hawaii show marked differences between the basaltic rock types and help define stratigraphy in the hole. Water-saturated densities of subaerial lava flows (occurring above 1079 m depth) have the broadest range because of the large density variation within a single lava flow. Water-saturated densities commonly range from 2.0 to 3.0 with an average of 2.55 ?? 0.24 g/cc. Dikes and sills range from 2.8 to 3.1 g/cc). Densities of hyaloclastite commonly range from 2.3 to 2.7, with an overall average of about 2.5 g/cc. The low-density of most hyaloclastite is due primarily to palagonitization of abundant glass and presence of secondary minerals in the interstices between fragments. Four principal zones of pillow lava, separated by hyaloclastite, occur in the drill core. The shallowest (1983-2136 m) is paradoxically the densest, averaging 3.01 ?? 0.10 g/cc. The second (2234-2470 m) is decidedly the lightest, averaging 2.67 ?? 0.13 g/cc. The third (2640-2790 m) and fourth (2918-bottom at 3097 m) are high, averaging 2.89 ?? 0.17 and 2.97 ?? 0.08 g/cc, respectively. The first pillow zone includes degassed pillows i.e. lava erupted on land that flowed into the sea. These pillows are poor in vesicles, because the subaerial, one-atmosphere vesicles were compressed when the flow descended to deeper water and higher pressure. The second (low-density, non-degassed) pillow zone is the most vesicle-rich, apparently because it was erupted subaqueously at a shallow depth. The higher densities of the third and fourth zones result from a low vesicularity of only a few percent and an olivine content averaging more than 5% for the third zone and about 10% for the fourth zone. The uppermost hyaloclastite extending about 400 m below the bottom of the subaerial basalt is poorly cemented and absorbs up to 6 wt% of water when immersed. Progressing downward the hyaloclastite absorbs less water and becomes better cemented. This change is apparently due to palagonitization of glass and addition of secondary minerals in the deeper older hyaloclastite, a process favored by the increase of temperature with depth. The cementation is largely complete at 1800 m depth where the temperature attains about 20??C. The zone of freshest, uncemented hyaloclastite represents the weakest rock in the drill hole and is a likely level for tectonic or landslide disruption. ?? 2001 Published by Elsevier Science B.V.

Journal of Volcanology and Geothermal Research↗