Feb 10, 1994 - We assumed that a complex pressure time history could cause ...... network on Mount St. Helena Volcano and observations of long-period ...
JOURNAL OF GEOPHYSICAL
RESEARCH,
VOL. 99, NO. B2, PAGES 2637-2652, FEBRUARY
10, 1994
Array measurements and modelingof sourcesof shallow volcanic tremor at Kilauea Volcano, Hawaii Peter Goldstein
EarthSciences Division,LawrenceLivermoreNationalLaboratory, Livermore,California Bernard
Chouet
U.S. Geological Survey,MenloPark, California
We showthat resonances in near-surfacelayers (path effects),or a combinationof near-surfaceresonances and sourceresonances(sourceeffects),could be responsiblefor the sharp spectral peaks found in the records of gas-piston activity and volcanic tremor recorded near Puu Oo crater on Kilauea Volcano, Hawaii. Two types of sources are found to be compatible with observations and can not be distinguished using our data. In one model, impulsive, explosion point sources are embedded in the Puu Oo structure. In the other model, impulsive, point sources of pressure excite the resonances of a fluid-filled
crackembeddedin that structure. Both modelsrequireshallowsources(z < 100 m) and display a strong dependence of the radiated spectrum on source depth over a depth range comparable to the overall thickness of the surficial layers. Based on these results, it may be possible to track changes in source depth through temporal changes in observed spectra, a potentially useful tool for monitoring volcanic activity at this site. We estimate the depth and spatial extent of volcanic tremor and gas-piston activity using two dense arrays with respective apertures of 800 and 120 rn located near Puu Oo. Measurements of slowness
(ray parameter)and azimuth as a function of time clearly indicate that the sourcesof volcanic tremor and gas-piston activity are located beneath or in close proximity to the Puu Oo crater at depths shallower than approximately 1 km. Based on slowness and particle motion analyses we find that the records of volcanic tremor and gas-piston events at Puu Oo consist of a complex combination of body and surface waves.
INTRODUCTION
Over the last 10 years observations of seismic activity in the vicinity of and prior to volcanic eruptionshave provided an extremely useful tool for the mitigationof volcanichazards[e.g.,Malone et aI., 1981; Koyanagi et al., 1988; Watanabe,1988; Okada et aI., 1990; Chouet et aI., 1993]. Unfortunately, existingpredictivecapabilitiesare primarily empirical and may be misleadingin new environments.Therefore an important goalshouldbe an improvedunderstanding of the physics of volcanic sourcesand the signals they produce. The primary aim of this study is to improve our understanding of volcanic sources in the vicinity of Puu Oo crater, Kilauea Volcano, Hawaii, through detailed measurementsof spatial and temporal variations in near-sourceseismicwavefields. An additional objectiveis to determine how well we can constrain the spatial extent of a volcanic sourceusing
magma movements, they would provide a useful tool for hazard mitigation. Models of Volcanic Tremor
At present, there are two most commonly cited models for volcanic tremor. Those whose spectral features are controlled by resonating sources and those whose spectral features are controlled by path and site effects. Resonating source models have been the
more popularof these [e.g., Aki et al., 1977; Aki and Koyanagi, 1981; Crosson and Bame, 1985; Chouet, 1985, 1988, 1992; Feblet and Chouet, 1982; Feblet,
1983; Ferrazziniand Aki, 1992; Chouetet al., 1993] because they can explain the common spectral peaks of tremor observed at widely separated seismicstations
[e.g., Aki and Koyanagi, 1981; Chouetet aI., 1993] and can produce temporal variations in tremor spectra
[e.g., Shimozuruet aI., 1966; Kubotera,1974; Kamo
near-source dense-array recordings of seismic signals. et aI., 1977; Aki and Koyanagi, 1981; Schick et al., If accurate images could be obtained and related to 1982;Feblet,1983]. However,a numberof studieshave
presentedevidencesuggestingthat path and/or site effects could also be.responsiblefor some of the observed
Copyright 1994 by the American Geophysical Union.
spectralfeaturesof shallowvolcanictremor[e.g.,Omer,
Paper number 93JB02639.
1950; $v[inakami, 1960; Malone et aI., 1981; Gordeev
0148-0227/94/93JB-02639 $05.00
et al., 1990; Gordeev,1992]. Evidenceof path or 2637
2638
GOLDSTEIN AND CHOUET: SOURCESOF VOLCANIC TREMOR AT KILAUEA
site effects include
observations
of similarities
between
We use three-component array data to show that
shallowtremor and shallowexplosionwaveforms[e.g., the tremor waveforms consistof a complex combination of body and surface waves, but note that the surface McNutt, 1986, Gotdeer,1992]. We present new evidencethat suggeststhat many of the spectral peaks found in tremor at Puu Oo could be
waves have the largest amplitudes and durations. By itself, this observation restricts a significant portion due to resonances in near-surface layers(path effects). of the tremor source to shallow depths. Previous In particular we show that shallow impulsive sources studiesat Puu Oo [e.g.,Ferrazziniet al., 1991]did not can generate time seriesand spectra that are consistent find significant evidence of body waves because they with those that are observed. However, we also show averaged data fr.om long time windows. We also note that spectral peaks due to resonancesof shallow fluid- that the timing of gas-pistonactivity (the fountaining filled cracks could be hidden by such path effects. and drainback of lava associatedwith the bursting of Furthermore, we show that changesin sourcedepth can large gas slugs that have risen through the magma explain observed temporal variations in the time series conduit and break through ponded lava at the surface and spectra. [Swanson et al., 1979]is wellcorrelated with the largest We present a variety of additional evidence that amplitudetremor signals[Dietel et al., 1989; Chouet suggests that path effects contribute significantly to and Shaw, 1991]. This also suggeststhat shallow the waveforms and spectra of tremor at Puu Oo. sources are an important contributor to the tremor For example, we demonstrate that sources of tremor wavefield. at Puu Oo are shallow and confined to the vicinity Puu 0o Crater of Puu Oo by locating them using high-resolution The Puu Oo crater was formed shortly after the frequency-wavenumbertechniques. Evidence of shallow sources of volcanic tremor has also been found in a initial phase of the January 1983 Kilauea eruption number of previous studies. For example, McNutt and became the focus of degassing and fountaining
[1986]notedsimilaritiesbetweenB typeevents,volcanic activitiesuntil July 1986 [Wolfe et al., 1988]. The explosions, and tremor at Pavlov volcano, Alaska that activity then evolved into a more sustained, low-level suggest a shallow source. Waveform cross-correlation outpouring of lava from a new vent at Kupaianaha, studiesby Furumotoet al. [1990]and Yamaokaet al. 3 km downrift from Puu Oo. The activity at Puu Oo [1991]alsoindicatea shallowsourceof volcanictremor simultaneously changed to low-level degassingmarked by intermittent bursts associated with stronger lava beneath Izu-Oshima volcano, Japan.
-.--180 rn-I Oldlavaflow
-•:.,• ...... ..... .-.:...:.:.....:.:...:.:.:.:.:.:. ..:.:.: ..:...:.....:.:.... ... .....:......... • .
Dense Arrays
::::-::....•::-'-.:: ..-..:.:.:....:?• :.:..-.-:.........:.....
. .... .... -.. :.•:?::::-....:::-:•:;-::.:..;:::::-::.-::::::.::::.:.:: .::'• .• -' .. ............. . .....::..: ..............•...:•
. . . *•,
....•
.......:•27;• .....
30rn
lOOm
Fig. 1. Schematicillustrationof topography, magmatransportsystem,andearthstructurein the vicinity of Puu Oo crater.Potentialsourcemechanisms includeexplodingbubbles(grayovals)and/or a resonating fluid-filledcrack(rectangle)excitedby explodingbubbles.Tke upperboundaries of the magmatransport
system(thickdashed lines)areconstrained to beconsistent withsurface ventsandreliefinsideandoutside the crater.Interfacedepthsand layervelocitiesin the one-dimensional modelof Ferrazziniet al. [1991] are given at the left and right of the model, respectively.
GOLDSTEIN
AND CHOUET:
SOURCES OF VOLCANIC
TREMOR
AT KILAUEA
2639
movements confined within the vent. Continuous aperture array because of an apparent lack of signal weak tremor has characterized the persistent pressure coherence observed in the records at this spacing. The second array was deployed in a semicircular fluctuations in the magmatic system under Puu Oo since that time. configuration with a station spacing of 10 m and an Successive episodes of fountaining at Puu Oo apertureof 120m (Figure3). In additionto improving gradually built a spatter cone, which at the time of broadband coherence between adjacent elements, the the array experiments stood some 230 m above the size and geometry of this array was chosen so that surrounding lava fields and contained a cylindrical the statisticalanalysismethodof Aki [1957]couldbe crater with a diameter of 180 m and estimated applied and so that the array would sample roughly depth of 150 m. Two active vents were present at half a wavelength of Love wave at 5 Hz. The resulting the northeastern and southwestern perimeters of the small aperture array was very useful for measurements craterfloor (Figure 1). Visualobservations conducted of surface wave dispersion and modeling of the structure from
the
crater
rim
indicated
that
both
vents
beneathPuu Oo [Ferrazziniet al., 1991]. However,
acted
interchangeably as inlet and outlet to lava movement
the reduction in the array aperture may have been
during episodesof gas-pistonactivity [Dietel et al., unwarranted from the perspective of frequency-slowness 1989;ChouetandShaw,1991;FerrazziniandAki, 1992]. analysesfor two reasons. First, as indicated by bandDATA, INSTRUMENTATIONAND ARRAY DESIGN
pass-filtered records and narrow-band cross-correlation function from two stations separated by approximately
950 m (Figure4), the lackof coherence in the broadband During January-February1988 two densearrays of records was misleading; there were clearly coherent seismometers were deployedin the vicinity of Puu Oo arrivals propagating across the first array. We found crateron the eastrift of KilaueaVolcano[Dietelet al., similax correlation coefficients in a number of frequency
1989](Figure2). The primarygoalof thisdeploymentbands, especially in the vicinity of spectral peaks. was to obtain near-source dense-array recordings of Second, reducing the aperture of the array reduced seismicsignalsgeneratedby ongoingeruptive activity at its spatial resolution capabilities which are inversely Puu Oo. These experiments consisted of two separate proportional to the aperture. In other words, the main deploymentslocatedon a fiat expenseof denselycracked peak in the beam pattern of the second array is much pondedpahoehoelava flowscoveringan areaof 0.5km2 wider than that of the first array. Based on these I km west of Puu Oo. In the first deployment an observations, we suggest that naxrow-band or low-passarray was configuredas a mix of linear segmentsand filtered estimates of coherence should be used when grid-shaped patterns with an aperture of 800 m and determining array design. The array instrumentation consisted of 12 Mark approximatestationspacingof 100 m (Figure 3). After 10 daysthe arraywasreconfigured into a muchsmaller- Productsthree-component L-22-JD seismometers with 155ø20 ' 19ø30 '
155o10 '
I KILAUEA
-o
I CALDERA
PuuOø--Array (;yCrater
,x,,,
19o20 '
I
i
I
I
Fig. 2 Map of the southeastflank of Kilauea volcanoshowinglocationsof Puu Oo crater, densearrays, and earthquakesusedin this study. The earthquakesare indicatedby solid diamondsand are labeled with their correspondingJulian day. The island of Hawaii is shown as an inset.
2640
GOLDSTEIN
AND CHOUET'
SOURCES OF VOLCANIC
TREMOR
AT KILAUEA
ß ,K6 ß
ARRAY
,
ARRAY 2 /
I
ß •-'.• ;• K2,L4 ß
Puu 0o Crater
amo
lOO METERS I
Fig. 3. Map of Puu Oo craterwith densearrays1 (solidcircles)and 2 (dashedlines).StationsK2 of array 1 and L4 of array 2 were approximatelycolocatedand were used to comparespectra from both arrays.
-
i
i
i
[
i
i
I
i
I
i
i
i
i
I
i
i
i
i
I
i
i
i
i
o
>-• -2--
' I ' ' ' I ' ' •'jI ', ' ' I '"' '
0•• 2 •
' '"'
•-2
0
'"'""""'•11" "'*'"'•
20
40
"' 60
"1
80
l;me(s)
- ' I ' ' ' I ' ' ' _•1. ' ' ' I ' ' ' IAfK6.cc ' --
.• 0.5 •
.
•,,
•,,- -,,,..... . .......... ,,, ..,..,,,,..• .I; ..... ; .,•,.,-•
• -0.5 -40
-20
0
20
40
•me Log(s)
Yig.4. (Top)Band-pass-filtered and(bo[•om)cross-correlated records ofvolcanic •remoronJanuary22, 1988. The band passis from 2.7 [o 3.3 Hs and con•Mnsoneof •he prominen•peaksin mos• of [he Iremot
andgas-piston spectra.Theserecords arefroms•a•ions in array1 wi•h•hegrea•es• separation (s•a[ions
A6 andK6 wi•h m950m separation; seeYigure3). They show•ha[ •hereis significan[ correla[ionbetween •he signal, within the band passof in[eres•,across•he wholearray.
GOLDSTEIN
0
5
AND CHOUET'
lb
SOURCES OF VOLCANIC
115
20
TREMOR
25
AT KILAUEA
50
2641
55
40
Time(s) Fig. 5. Typical vertical component seismograms of tremor, gas-piston activity and earthquakes recorded on both arrays. Except for the earthquake on January 31, all the records are from stations K2 and L4
which are nearly colocaied(seeFigure 3). The Julian day/event type, local magnitude,and depth of each event is listed at the right of each trace (VT -- volcanictremor, GP -- gas-pistonactivity, EQ earthquake)(seeTable 1).
sensitivity of 0.5 V/cm/s and 36 Mark Products gas-pistonactivity, and a shallow(depth 1.1 km) and vertical component L-4A seismometerswith sensitivity relatively deep (depth 9.3 km) earthquake(Table 1).
of I V/cm/s.
All the instrumentshad a natural The locations of these earthquakes are indicated in
frequency of 2 Hz. Recording was performed using 12 Figure 2. GEOS digital seismographs,each consisting of a six-
channel,16-bit recorder[Botchertitet al., 1985]. All GEOS recordeddata at 200 samples/s/channel usinga commontime basewith an accuracyof I ms among all
SHALLOW
POINT
FLUID-FILLED
SOURCES
OR
CRACKS?
the channels.
In this section we describe our modeling of the Gas-piston activity and volcanic tremor were two ground motions and spectra of volcanic tremor using persistentsourcesof seismicactivity at Puu Oo during point sources and resonating fluid-filled cracks. Since the array experiments. Small earthquakes to the much of our modeling was motivated and constrained by south of Puu Oo were recorded as well. Typical our array measurements and spectra] analysis, we begin vertical componentseismograms for tremor, gas-piston by summarizing the most important results of those activity and earthquakes recordedat two approximately measurements and their implications for our modeling. colocated sensors are shown in Figure 5. Events Details of these measurementswi]] be presentedin later occurring prior to January 30 were recorded on the sections. larger aperture array. These include an example Our spectral analyses can be summarize(] by the
of volcanictremor and three earthquakes(Table 1). followingobservations:(1) individualstationsiteeffects Unfortunately,there was no gas-pistonactivity during can not explain the sharp spectra] peaks in volcanic the period of deployment of the first array. Data tremorspectraat Puu Oo, and (2) there are significant from the smaller array include an exampleof tremor, variations in tremor spectra as a function of time.
2642
GOLDSTEIN
TABLE
AND
CHOUET:
SOURCES
OF VOLCANIC
TREMOR
AT KILAUEA
1. Source Parameters of Earthquake, Tremor, and Gas-Piston Events Discussed in the Text.
Date*
Origin Time, UT
Latitude,
Longitude, M L
deg
Depth, km
deg
Earthquakes
Jan. 23, 1988 (023) Jan. 26, 1988 (026) Jan. 27, 1988 (027) Jan. 31, 1988 (031) Feb. 10, 1988 (041)
0936:37 0728:56 0124:36 0028:54 1326:58 '•
19.3395 19.3453 19.3248 19.3457 19.3517 VolcanicEvents
-155.1190 -155.1283 -155.0572 -155.0978 -155.0415
Jan. 22, 1988 (022/VT) Jan. 30, 1988 (030/VT) Feb. 2, 1988 (033/GP)
0359 0259 0449
19.3915 19.3915 19.3915
-155.0992 -155.0992 -155.0992
1.8 2.2 3.9 1.8 3.4
1.72 6.19 5.22 1.08 9.27
*Each date is followedby its correspondingJulian day.
Results from our array measurementsindicate that resonances (path effects)couldbe responsible for many (1) mostof the coherentenergyin the tremorwavefields of the sharp peaks in tremor spectra. In this section we is comingfrom shallow(z < I km) sourcesin close compare results from this model with those predicted proximity to Puu Co crater, and (2) the observed by the resonating crackmodel[e.g.,Chouet,1992].We wavefieldscontainsignificantamountsof both body and focusour modelingon shallow(z < I km), nearby(r • surface waves. I km), sourcesbut haveconsidered a few examplesof The aboveobservationsand our experiencemodeling deeper more distant sources. Our simulations include
explosionwaveforms[e.g., Goldsteinet al., 1992] both body and surface waves. Based on the above motivated us to test the hypothesis that near-surface observationsand the followingsimulations,we conclude
- Surfoce --
--
- Z=100
iSurfoc'e ' Z=100 m
-
ß
- ou•
-_DATA - o•/•P DATA
I
0
2
,
I
4
,
I
$
,
I
,
8
Time(s)
-
1
4
0
2
4
$
8
Frequency (Hz)
Fig. 6. Comparisonof groundvelocitiesand powerdensityspectraof shallow,impulsive,explosionpoint sourceswith those of observedtremor and gas-pistonevents.Each syntheticis a vertical component calculatedat a distanceof I km usingthe structureof Ferrazziniet al. [1991].The top two traces correspondto a singlesourceat the surfaceand 100 rn depth, respectively.The middle three traceswere obtainedby delayingand superposing multiplecopiesof a singlewaveformfrom a sourceat the surface, 100m, and 300 m depth,respectively. We usedrandomtime delaysbetween0 an.d2 s. The bottomtraces are examplesof tremor on the first array and gas-pistonactivity on the secondarray.
GOLDSTEIN
AND CHOUET:
SOURCES OF VOLCANIC the
TREMOR
data
for
AT KILAUEA
sources
between
2643
that resonancesin the surficial layers are a plausible
with
the
surface
and
explanationfor the observedspectralfeaturesbut that
100 m depth (Figure 6), and suggestthat a complex
time-dependent excitation of shallow point sourcesin a a resonating source couldalsobe present. We tried to determineif the spectralpegksin layered medium might explain the observations.
volcanic tremor and gas-pistonrecords from Puu Oo Shallow,resonating,fluid-filledcrackmodels[Chou•t, were sourceor path related using forward modeling. 1988, 1992]can alsoexcitenear surfaceresonances and We simulatedpath effectswith impulsive,point-source generate spectra that are consistent with the observed
explosionreflectivity syntheticseismograms [•f•'ller, spectral peaks. For example, Figure 7 shows synthetic 1985]in thelayeredvelocitystructureof Ferrazzini•t al. seismogramsand spectra of the vertical component of [1991]IFigure 1). We computedverticalcomponent ground velocity calculated in the near field of a vertisyntheticsfor a variety of sourcedepthsbetween0 and cal fluid-filled crack of rectangular shape embedded in 1 km and for epicentraldistancesbetween800 and 1200 a homogeneousmedium and the layered structure dem. We found that both time series and spectra were
rived by Ferrazziniet al. [1991](Figure 1). This cal-
quite sensitiveto sourcedepthwith muchlesssensitivity to distance. Spectra of sourcescloseto and within the uppermostlow-velociWlayersvariedrapidly with depth and have spectral features that are similar to those of
culation was performed using the discrete wavenumber
velocities in the homogeneous half space are the same
the data (Figure6).
as those found
We modeled the nearly continuoustremor and gaspiston data by assuming that a number of point sourceswere occurring in a random time sequence.We superposedidentical syntheticwaveformswith random time delaysbetween0 and 2 s. Waveformsand spectra of the resulting time seriesare in qualitative agreement
(Figure 1). For the fluid in the crackwe assumedan
0
2
4
6
8
1•)
method [Bouchon,1979]and propagator-based formalism of Chouet[1987]. Compressional and shearwave at the bottom
of Ferrazzini
et al's model
acoustic wave speed a- 1.25 km/s, a value compatible with a basaltic magma containing void fractions of gas
rangingup to a few percent[Chouet,1992]. A sketch of the source-receivergeometry used in our model is given in Figure 7. Crack excitation is provided
4
Fig. 7. Comparisonof ground velocities and power density spectra of a vertical fluid-filled crack with those of observedvolcanic tremor and gas-pistonevents.Each synthetic is a vertical component calculated at a distance of 1 km and an azimuth of 25 ø. The top two traces correspondto a single step in pressure at the top of a crack that intersectsthe free surfacein a homogeneoushalf-spaceand the layered structure
of Ferrazzir•iet al. [1991],respectively. The next three traceswereobtainedby delayingand superposing
multiple copiesof crack syntheticsin a layered medium. The first two correspondto sourcesat the top
(z = 0 m) and middle(z -- 100 m), respectively, of a crackthat intersects the free surface(h -- 0 m). The third tracecorresponds to a crackwith a sourceat its midpoint(z -- 400 m) whosetop is at
a depth of 300 m. Medium type, number of sources,and source and crack depth are identified at the left of the spectra. Time delays were the same as those used to superposepoint sources. The bottom traces are examplesof volcanictremor on the first array and gas-pistonactivity on the secondarray. The source-receivergeometry is shown in Cartesian coordinates as an inset. The dark patches on the crack axis represent the locations and extent of the pressuretransients triggering crack resonance.
2644
GOLDSTEIN AND CHOUET' SOURCES OF VOLCANIC TREMOR AT KILAUEA
by a suddenstep in pressureappliedover a small area surficiallayersagreequalitatively with our observations of the crack wall located on the main crack axis near (Figure7}. Again,we find that the spectrafor sources the top, center,or bottomof the crack.The top of the within the uppermost low-velocity layers depend on crack is buried at a depth/• and the receiveris located source depth and structure. Differencesbetween the at a distancer from the epicenterand azimuth • from shallow sourcesynthetics and observedspectra could be accountedfor by changing the crack geometry or the vertical crack plane. We assumed a crack with horizontal and vertical
dimensionsof 100 and 200 m, respectively,on a plane
acoustic velocity in the fluid. In contrast, spectra correspondingto resonatingcracksbelow the surficial
alignedwith the rift axis as definedby the eruption layers(a=105, 300, and 500 m} are fairly insensitive
fissuresin Figure 1, and choser = 1000m and • = 25ø to source depth and have much less bandwidth and for the receiver. This receiver position coincideswith complexity than those of the data. These results that shallow(a < 105m} resonating cracksmay the location K2 of the large array or location L4 of suggest the small array (seeFigure 2). We considered four be a component of a complex sourceand structure at crack depths, /• = 0 m, /• = 105m, /• = 300m, and Puu Oo. In combination with results from the following array h = 500 m, and three sourcedepths along the main and spectral analysis, the above results suggestthat crackaxis,z = h, z = h+ 100m, and z = h+ 200 m. The volcanic tremor and gas-piston activity at Puu 0o syntheticswereobtainedovera bandwidthof 10 Hz or crater are generatedby a near-surfacesourceand that greaterusinga windowof 6 s that includesthe entire path effects or a combinationof sourceand path effects historyof groundmotionat thereceiver(e.g.,Figure7). could explain the observedspectral complexity. We assumedthat a complex pressuretime history could cause sequentialexcitation of the crack. We SPECTRAL ANALYSIS modeledthis complexsourceexcitationby superposing identical waveforms with the same random time delays In this section we compare array-averagedspecused to superposepoint sourcesynthetics.Waveforms tra from a variety of seismicevents and show that and spectra of resonating cracks that penetrate the individual-station site effects are an unlikely explana-
Hdrmoni'c Tre•or
LargeArray I
Hdrmoni'c Tre•nor SmallArray
'
' I ' I ' I ' I ' i ' i '022/v 1
i
i
i
' /•'
'
'
'
'
'
'
'
'
'
'
L
4x,
'
I
'
I
'
I
'
I
Ods-Pi•on E•,ent ' ' /•' -
mo,, ^rro, -'
' -•
q-'--
I
'
I
Earthquake
^'
I
/ •/'
I
'
I
'
[ Earthquake
'
'
x_.•
I_Large Array•_•_.••,•/ '
i
I
vvuv
'
I
023/œo '
I
'
I
'
I
'
I
'
V
.............. Eo'rthqu(•ke ........
o•Vm •=1,8
SmallArray _
'"'- '
I
'
__
I
Earthquake
Ii
A
o•V•o
Small Array 2
3
4
5
6
7
8
9
Frequency (Hz) :Fi•. 8. A•a,y-avera,•ed power density spectra of ve•tica,] components of events sbow• i• :Fi•u.re 5. "['be spectrum of each station is obtained using It, overlapping, 30-s long windows and 40 s of data.
GOLDSTEIN
AND
CHOUET:
SOURCES
OF VOLCANIC
TREMOR
AT KILAUEA
2645
tion for mostof the largespectralpeaksin the volcanic array (event 022). If a site effect were responsible tremorandgas-piston eventsobserved at Puu 0o. We for the sharp peaks in the tremor spectra, we would alsoshowthat there are significanttemporalvariations expect to see it in the spectra of earthquakesas well. in tremorspectraandsuggest that thesevariationsmay This argument is lesspersuasivefor tremor recordedon the secondarray (events030 and 033) becausethere be dueto temporalvariationsin sourcedepth. With the exception of an earthquake on January 23, are some notable peaks in the spectra of earthquakes spectralestimateswere derivedby computinga power recordedon that array (events031 and041). However, densityspectrum[Jenkinsand Watts,1968]for each the locations of spectral peaks of tremor recorded on vertical component of motion using 11 successive30- the first and second arrays are well correlated, while s windows with 1-s overlap, and then averaging the the peaksin the spectraof earthquakesrecordedon the resultingspectraacrossthe array. The sameprocedure secondarray appear to be more randomly distributed. Evidence of the temporal variability of tremor was applied to the January 23 earthquakeexceptthat spectra at Puu Oo is given'in Figure 9. In this figure we the overlappingwindowswerelimited to a length of 10 s see the tremor spectrum evolve from one with a single because of the smaller signal duration of that event. fairly sharp dominant peak to a complex spectrum We looked at temporal variations in tremor spectraby containing multiple peaks in the 2.5 to 4.0 Hz band. comparingresultsfrom consecutive40-s windowswith Variations similar to these could easily be replicated
20 s of overlap.
Our primary evidence that sharp spectral peaks by superposingspectra from successiveexplosionpoint in tremor spectra are not site effects is indicated sourceswith varying depth (e.g., Figure 6), or from in Figure 8. In this figure we see that the array- successiveimpulseresponsestriggered at variousdepths averagedspectra of earthquakesrecorded on the first in a shallow fluid-filled crack embedded in the Puu Oo
array (events023, 026, and 027) are fairly smoothin structure(e.g.,Figure7). comparisonto those of tremor recorded on the first
'
I
'
I
'
I
'
I
'
I
'
I
'
I
'
I
0,4 O2 0.0 0.4
02
0.40,2o. o.
o. o 0.t
O•0.2L .
0.• 0.1• 0.4-
0,2
0,4'
0 0
Fig. 9. Array-averaged power density spectra of consecutive, 40-s long, overlapping windows of volcanic tremor recorded on January 22, 1988. Each spectrum is obtained using 11, overlapping, 30-s long windows. There is 20 s of overlap between successivespectra.
2646
GOLDSTEIN
AND CHOUET'
SOURCES OF VOLCANIC
TREMOR
AT KILAUEA
Multiple SignalClassification technique[Schmidt,1986; Goldstein and Archuleta, 1987], and seismogram alignment[Goldstein,1988; Goldsteinand Archuleta, In this section we use frequency-slownessanalyses, 1991a]. The advantagesof these techniquesover SPATIAL
EXTENT
OF TREMOR
AND DEPTHS
SOURCES
more conventional array signal processing techniques
particle motions, and ray tracing to estimate the spatial extent and depths of the sources of volcanic tremor and gas-piston activity beneath Puu Oo. We perform similar calculations with data from small earthquakes to identify any trends in ray parameter and azimuth that might be due to path or site effects.
has been demonstratedin detail by Goldstein[1988] and Goldsteinand Archuleta[1991ab]. We beganby
aligning the peak of a prominent arrival in the vertical component waveforms of each event. For example, the first peak in each waveform was used for most of the We measuredtotal slowness(ray parameter)and earthquakes. Volcanic tremor and gas-piston records azimuth as a function of time using MUSIC, the were scanned in record section format to identify
'
I
'
'
I
'
'
--
0.5-
'
I
'
-
--0.5 •
)
--
_
_
,
-1.0 ' f=Z7 •-'-...'d•b= 1_2 • '
1.0 , • ,
, • ,
.
I , , I , f= 3.1 db•b= 1,5
sx=-0320 sy=-OJ'20
sx=0.000sy=-0240
f-- 3g dbd>=
f= 3.5 db•b= 1
sy= sx=-0280 sy=-O.•20 sx=-0',360
.
-
0.0
'
-
0.••••.•
0.5
I
-
--
(
-0.5 -
,
I
,
,
1.0 , I '
I
,
,
,
,
I
,
,
I
,
SX=
' I '--
.
I
,,
--
0.5 ..
ß
0.0 .
_
-0.5 .
-
,
-1.0
4--.,
,
b= 4.7 dbob=2.2
sx=-0.360sy=-•1-60 -
0.5
,
,
0,0
I
0,5
,
1.0-1.0
-0,5
0,0
0.5
...
-
0,0
I
-1,0 -0.5
'
•
•
-0,5 .
-1:-01,0 -0,5 0.0 0,5
-
'
1,0-1,0
-0,5
0,0
0,5
1,0
x componen[ of slowness Fig. 10. Narrow-band frequency-slownessspectra and the correspondingstacked spectrum from a single
windowof volcanictremor on January22, 1988.The centerfrequency(/), peak spectralamplitudein dB abovethe baselinespectrallevel (dbab;seetext), and corresponding components of slowness at the peak (sx, sy) are listedat the top of eachnarrow-bandspectrum.The stackedspectrumis at the bottomof the secondcolumn.It is labeledby bandwidth(b) insteadof centerfrequency.Sourcedepthsand azimuths were estimated from slownessescorresponding to peaks in the stacked spectra.
1,0
GOLDSTEIN
AND
CHOUET:
SOURCES
OF VOLCANIC
TREMOR
AT KILAUEA
2647
prominentarrivals. The slowness/rayparameterof a
slownesswas converted to ray parameter and azimuth
reference arrival was determined by fitting a plane wave to its time delays across the array. After alignment, we computed slowness spectra during successive,2.5s long time windows with 2.0 s of overlap between adjacent windows. Slowness spectra were estimated over a bandwidth from 2 to 6 Hz using slownessstacking
using
q•- 9/2- arctan(Sy/Sx)
whereP is theray parameter,$x and$y arethe z (east) of slowness, respectively, and [Spudichand Oppenheimer, 1986]. A typical set of andy (north)components • is the apparent direction of propagation or signal narrow-band(0.4 Hz bandwidth)slowness spectraand azimuth
their corresponding stacked spectrum, for one time window, are shown in Figure 10. Vector slownesswas estimated from the locations of peaks in the stacked spectra provided their amplitude was at least i dB
The baseline spectral level is defined as two standard deviations above the mean spectral level. Based on both simulations and experience with data, this choice filters our peaks due to noise and signals with poorly constrained slowness. Both components of slowness
were initially allowedto range over a 4 s/km window on
the
reference
arrivals
slowness.
clockwise from north.
SourceAzimuths and Spatial Extent
higher than our baselinespectra level (dbab _• 1).
centered
as measured
These
The azimuthal distribution of earthquake,volcanic tremor, and gas-pistonsignals,measuredduring a 40s time window, are indicated in Figure 11. This plot shows that volcanic tremor and gas-piston signal azimuths are concentrated in a cone approximately 25ø wide pointing in the direction of Puu Oo with a slight trend in the down rift direction towards Kupaianaha. In contrast, earthquake azimuths axe more evenly distributed and typically span a much larger range than
windowswere narrowedto 2 s/km windowswhen we those of gas-piston and tremor activities. The cluster found them sufficient to capture all the significant in azimuths near 270ø for the January 31 earthquake arrivals. The initial range was chosen to be large is due to the superposition of this small event on enough to include both body and surfacewaves. Vector a background of volcanic tremor. The number of 2,0
'
'
'
'
I
'
'
'
'
o•- La, rqe, Ar?y I ....
2,0,
I
'
I
....
'
'
'
I
'
I
....
'
'
'
I
i
.
-
....
1-'Z_
•m-Harmohio Tremor 0,5S,maJJ Ar[ay
•
••••,•••)k • •/HT-_
(;as-Piston Even'[
-
•
all Ar,r ay Sin,
-
-,-'
Earthquake
-
!,_
"
'
'""
I
'
'
'
'
{ - '
'
'
'
{
'
'
'
'
I
""
'
r'L - Earthquake •c•
- ' •k,a, rqe, Ar[ay I '
'
m7/EO-: '
'
I
'
'
'
'
I
'
J '
'
I
'
'
'
lm- E•hquake
2,0 -
,Ar[ay• ....
i ....
• ....
•
-
om-,Smali, Array•, , , , • , , , , • , , , .•
1•
1.•
!
•
i
20o
Azimuth (deg) Fig. 11. Plots of ray parameteras a functionof azimuthfor mostof the eventsshownin Figure5. The first 40 s of eacheventweredividedinto consecutive 2.5-slongwindowswith 2 s of overlap.Azimuth is measuredclockwise from north. The azimuthfrom eacheventis indicatedby the vertical]Line through
each corresponding data set.
2648
GOLDSTEIN AND CHOUET: SOURCES OF VOLCANIC TREMOR AT KILAUEA
Ray Parameter(s/km) 0.5 1.0 1.5
coherent arrivals in the volcanic tremor and gas-piston records and their clustering in azimuth around the direction of Puu 0o indicate that sources localized in
closeproximity to Puu 0o crater are responsiblefor
I
the observed volcanic tremor and gas-piston activity.
0.2'
Assumingthe tremor sourcesare at a distanceof 1 km, the approximatedistancebetweenthe arrays and Puu 0o crater, the horizontal extent of the source
0.4-
region(œ----re) is approximately 400m. This is roughly twicethe sizeof Puu Oo crater (Figure1). Basedon (ee of GoldsteinandArchuleta[1991a]),the uncertaintyin
0.8'
this estimate of spatial extent is approximately 100 m for the larger array and 300 m for the secondsmaller array. Differencesin precisionof the arrays are due to differencesin aperture, averagestation spacing,and uncertainties in delay times used to align the data.
I
I
I'
I ,
0.6F
1.0"
Fig. 12. UpgoingP and S waveray parametersas a function
of sourcedepth, in the model of Ferrazziniet al. [1991]. S wave ray parameters are indicated by the thin curves. These curvesprovide a bound on the maximum sourcedepth because reflections from sources within the top two layers could have a smaller ray parameter than a direct wave.
Ray Parameters and SourceDepths
Most of the volcanictremor ray parameters(Fig- or in the vicinity of the shallowsurficial layers. We base ure 11), especiallythosewith azimuthscorrespondingthis conclusion on our observation that only shallow to the direction of Puu Oo, range betweenslownesses of sources are expected to generate significant surface about0.6 and 2.0 s/km. The low endof thisrange,less waves and that surfacewaves are a large componentof
than approximately1.4 s/km, is incompatible with sur- the observed tremor wavefield. However, we also show face waves in low-velocity layers such as the onesfound that a significant amount of body waves are present by Ferrazziniet al. [1991].Basedon this observation,and note that their correlation with the smallest ray we assumethat the low slownessarrivals are body waves parameterswe observesuggeststhat there may be some and use this range of ray parameters to constrain the contributions to the tremor wavefields from sources as maximum sourcedepth as follows. First, we use Snell's deep as I km. We began by looking for body or surface wave law to relate the observed ray parameter $ to a corpolarizations in individual three-componentrecordings respondingdepth (z) for an upgoingbodywavesource but found significantvariability from station to station; as not suprising in light of the complex signal properties we have identified thus far. In order to identify those features that were common to most or all of the stations, we compared record section plots of the product of band-pass-filtered vertical and radial componentswith where ci and hi are the layer velocitiesand thicknesses the product of band-pass-filtered radial and transverse (Figure1), and r is the epicentraldistancebetweenthe components. In these plots, P and S V waves should source and receiver, which is fixed at I km. We only appear primarily on the vertical[radialrecordsas a considerdirect upgoingwaves becausethey correspond pulse or series of pulses that are consistently positive to the greatest sourcedepth for a given ray parameter. or consistently negative. Rayleigh waves should also Sourcedepth is plotted as a function of P and $ wave appear primarily on the vertical[radial recordsbut ray parameters in Figure 12. Given the lower cutoff should oscillate about zero. SH and Love waves should
in observedray parametersof about 0.6 s[km, this be observedon the radial[transverserecordswith
restricts a P wave source to a depth shallower than about 250 m. The maximum depth for an S wave source is about I km. Uncertainties in maximum sourcedepth, based on uncertainties in slownessat the large array, are approximately 100 m for P waves and 300 m for S waves. These results clearly show that the source of volcanic
tremor
at Puu
Oo is close to or within
the
upper kilometer of the crust.
consistently positive or consistently negative motions. We differentiate $H and Love wavesbased on expected differencesin their apparent velocities. We found that the polarization of a few arrivals were consistent acrossmost of the array. For example, in Figures 13 and 14 we show 2-Hz low-pass-filtered
vertical[radial and radial[transverserecord section plots for the first 20 s of the gas-piston event on February 2. Based on its oscillatory waveform, and
Particle Motions and Wave Types
large slowness (• 1.5 s/km), the arrival identifiedby
In this sectionwe use array measurementsof particle motions to identify wave types and argue that most of
the subvertical
the tremor
arrival marked by the line between 12 and 14 s in
source or sources must be confined
to within
line in the vertical-radial
record
section
(Figure13) is a Rayleighwave.Similarly,the impulsive
GOLDSTEIN
AND CHOUET:
SOURCES OF VOLCANIC
TREMOR
AT KILAUEA
2649
the most vigorous bursting of gasbubbles at the surface
of Puu Oo crater [Dictelet al., 1989;Ohouetand Shaw, 1991]. This suggests that point sourcesare associated with the explosionof large gas bubbles, someof which have risen through the magma conduit to the surface. This model may have important implications for monitoring. Our observation that the spectra of point
sourcesare sensitiveto depth (Figure6) suggests that it may be possibleto determine the depth of the most prominent source of energy from a single record. It
mightevenbe possible to trackchanges in sourcedepth through temporal changesin observedspectra. If such changes in source depth are correlated with and occur prior to changes in volcanic activity they would be a useful tool for monitoring. Even if correlations between activity and source depth do not exist, observationssuch as these provide useful indications of the subsurface plumbing and magma dynamics. For example, our observation that most of the sourcesare located in the upper few hundred meters beneath or in close proximity to Puu Oo with sources outside this region trending downrift toward Kupaianaha suggeststhat most of the magma dynamics is occurring in shallow conduits or chambers beneath Puu Oo with some activity extended in the direction of Kupaianaha. This is consistentwith visual observations 0
/5
10
1/5
20
Time (s)
Fig. 13. Record sectionplot of the product of the vertical and radial components of motion. Each component was
band-pass filtered below 2 Hz before multiplication. The solid line near 6 s indicates the start of a Rayleigh wave.
the radial-transverse recordsection(Figure 14) has a slowness of approximately0.9 (muchtoo low to be a Love wave) and is an S/-/wave. Althoughvolcanic tremor records from both arrays were quite variable, they both had consistentphasesthat could be identified as body or surfacewaves;with surfacewavesbeing more common. Based on this analysis, volcanic tremor at Puu Oo appears to be dominated by sourceswithin or in the vicinity of the surficial layers, although sources as deep as 1 km may also be contributing. DISCUSSION
In this section we discuss plausible mechanisms for and implications of point source models of shallowvolcanic tremor. We compare advantagesand disadvantagesof thesemodelswith those of resonating
fluid-filledcrackmodels[e.g.,Chouet,1988,1992].We provide a few additional observations that should be
useful constraints on further modeling and conclude with some suggestionsfor future work.
I
0
/5
10
I
[
I
I
I
15
I
20
Time (s)
Mechanism and Implications of Point Sources
Fig. 14. Record section plot of the product of the radial and transverse components of motion. Each component was Visual observationsindicate that the large ampli- band-pass filtered below 2 Hz before multiplication. The
tudes in gas-pistonrecordingsare well correlated with
solid line near
12 s indicates
the start
of an $H wave.
2650
GOLDSTEIN
AND CHOUET:
SOURCES OF VOLCANIC
TREMOR
AT KILAUEA
of the Puu Oo crater and Kupaianaha lava pond made at the time of this experiment. A point sourcemodel for volcanicsignMsat Puu Oo also has implicationsfor the rheology and composition of the magma. For example, if gas bubbles are the primary source, their frequency and magnitudes suggest a significant percentage of gas in the magma composition. The ability of the magma system to transport these gas bubbles should provide some constraint on magma viscosity.
spectra. These models are also consistent with visual observationsof magmatic activity in Puu Oo crater. Based on the above observations, either a model consistingof point sourcesembeddedin shallowsurficial layers or a model consistingof point sourcesin a fluid-
Point
Additional
Sources Versus Fluid-Filled
Cracks
The primary advantage of the point sourcemodel of
shallow tremor is its simplicity and ability to explain the observeddata. A potential disadvantageis that it does not provide a formal descriptionof the coupling between processesin the gas, liquid and solid. Such a description would aid in our understanding of the volcanicsourceprocessbut may not be identifiablefrom this data. This model may also be inappropriate when well-correlated spectral peaks are observedat stations with large separations/ncethey would most likely have different path effects;but this is clearly not the casefor
filled crack embeddedin shallowsurficiallayers appear to be appropriate for observations of shallow volcanic tremor in the open vent of Puu Oo. The selectionof one of thosemodelswould require a more extendedcoverage of seismic stations
than
Observations
available
at Puu Oo.
and Future
Work
During the processof modeling shallow tremor and gas-piston activity we found that shallow attenuation
mustbe fairly high (Q, • 20) to explainthe relativesize
of spectralamplitudes above6 Hz (e.g.,Figure6). We also found that long-term averagesprovide a different perspective than our analysis because surface waves,
whichhavelargeramplitudes,are emphasized.This who I0180 s) windows,found that surfacewavesdominated
their observations. We looked at 2.5-s long windows and found a significant amount of body waves. Based on the above observations,further comparthe Puu Oo arraydata (Figure4). Strictly speaking, the fluid-driven crack model is a isons of impulsive point source and fluid-filled crack composite of a triggerelement(the energysource)and models in a variety of volcanic settings are warranted. a resonator(the fluid-filledcrack). In the modelswe Additional observationswith large aperture arrays or considered,the trigger is again a point sourcewith the multiple arrays would be most usefulfor examining the essentialdifferencefrom the earlier model being that spatial as well as temporal characteristicsof volcanic the point source is now embedded in the fluid rather sources.The developmentof a model of time-dependent than the elastic solid. Thus the fluid becomes an active gas bubble generation and transport that is consistent participantin the generation of the elasticwavefieldwith the short-term rise and fall in activity within an via the process of acoustic resonance in the crack. ifidividual gas-pistonevent and the longer-termperiodThe primary advantage of this model is the direct link icities thai were observedbetweenmultiple gas-piston
betweenthe observedwavefieldand the fluid dynamics. Another advantageis the ability of the model to explain commonpeaks in the spectra of tremor seenat widely spacedstations[e.g.,Aki and Koyanagi,1981; Chouet
events is also needed. CONCLUSIONS
We have suggestedthat point sources,within and
et al., 1993]. When the crackpenetratessurficiallow- in the vicinity of shallow surficial layers, can explain velocity layers,path resonances in thoselayerscombine observationsof shallowvolcanictremor and gas-piston with sourceresonancesto produce a broadband multi- activity at Puu Oo crater, Kilauea volcano, Hawaii.
peakedspectrum.Whenthe sourceisremoved fromthe Resonating fluid-filledcracks[e.g., Chouet1933, 1992] surficiallayers,lessenergyis pumpedinto thoselayers that penetrate the surficial layers can also excite nearand the spectrumdisplaysa narrowerband dominated surfaceresonancesand would be difficult to distinguish
by sourceresonances. Thus a shallow(h < 100m, z < from point sources. Waveforms and spectra of both 100 m) resonatingcrackcouldbe presentat Puu Oo of these models are very sensitive to source depth, an but probably could not be distinguishedfrom shallow observation that may be useful for volcano monitoring. impulsive point sourcesbecause our arrays did not These observationsalso suggestthat path effectsshould providea wideenoughcoverageto answerunequivocally be understood and accounted for before empirical the questionof whetheror not sourcepeaksare present models for predicting the size and timing of eruptions are used in new environments. in the observedspectra. Measurements of ray parameter and source azimuth The point sourceand fluid-drivencrack modelsin a layered structure can both explain the waveforms, indicate that both tremor and gas-piston events are frequencycontent,ray parameters,azimuths,and wave localized beneath or in closeproximity to Puu Oo crater types of the observedtremor and gas-pistonsignals and are due to a nearly continuoussuccessionof sources. at Puu Oo. Assuming the position of the source We show that most of the sources must be located
and/or geometryof the crackcan vary, both models within or in the vicinity of the suficial layers to explain can also explain time-dependent variations in these the composition of the observedwavefields. Body wave
GOLDSTEIN AND CHOUET: SOURCES OF VOLCANIC TREMOR AT KILAUEA
2651
Chouet, B., A seismic model for the source of long-period events and harmonic tremor, in IAVCEI Proceedings in Volcanology,vol. 3, Volcanic Seisinology, edited by Based on ray parameter and particle motion P. Gasparini, R. Scarpa, and K. Aki, pp. 133-156, Springer-Verlag, New York, 1992. analyses,tremor and gas-pistonrecordingsat Puu Oo Chouet, B., and H. R. Shaw, Fractal properties of tremor consistof a complex combinationof body and surface and gas-piston events observed at Kilauea Volcano,
ray parametersconstrainthe maximumsourcedepthto about
i km.
waves. Analysis of long-term averagesof these signals Hawaii, J. Geophys. Res., 96, 10,177-10,189, 1991. obscuresthe body wave contributionsbecausethere are Chouet, B. A., R. A. Page, C. D. Stephens, J. C. Lahr, and J. A. Power, Precursory swarms of long-period events at more surfacewaves and they have greater amplitudes. Redoubt Volcano(1989-1990), Alaska:Their origin and Based on these observations,we suggest that more use as a forecasting tool, J. Volcanol. Geotherm. Res., in press, 1993. detailed measurementsof the spatial extent of similar R. S., and D. A. Bame, A spherical source model volcanic sources are possible using larger aperture Crosson, for low-frequency volcanic earthquakes, J. Geophys.Res., arrays or pairs of arrays and that the choiceof array 90, 10,237-10,247, 1985. aperture shouldbe basedon narrow-bandor low-pass- Dietel, C., B. Chouet, K. Aki, V. Ferrazzini, P. Roberts, filtered coherence estimates.
If such measurements
can
be used to estimate the extent of volcanic sources, they
may improveestimatesof volcanichazardsby helpingto constrain the amount of magma and available pathways to the surface.
and R. Koyanagi, Data summary for dense GEOS array observations of seismic activity associated with magma transport at Kilauea Volcano, Hawaii, U.S. Geol. Surv. Open File Rep., 39-113, 1-171, 1989. Fehler, M., Observations of volcanic tremor at Mount St. Helens Volcano, J. Geophys. Res., 88, 3476-3484, 1983. Fehler, M., and B. Chouet, Operation of a digital seismic network
Acknowledgments. Paul Kasameyerdeservescredit for motivating the first author to emphasizehypothesis testing in his work and for carefully reviewing this manuscript. We are also grateful to Phil Harben, Dave Harris, Howard
Patton, Steve Jarpe, Mike Fehler, and Yoshiaki Ida for their thoughtful reviews of this manuscript. Charles Carrigan, Bill Walter, and George Zandt provided many helpful comments. Christopher Dietel was especiallyhelpful in dearchiving data and answering questions about the field experiments. We are grateful to Keiiti Aki, Valerie Ferrazzini, Christopher Dietel, Peter Roberts, and Robert Koyanagi for their participation in the field experiments. We also thank David Okita of Hilo, Hawaii, for providing superb helicopter assistanceand a constant helping hand in the field. This work was supported in part by the National Science Foundation under grant EAR-8618107, through funding from the University of Southern California, and under the auspicesof the U.S. Department of Energy by the Lawrence Livermore National
Laboratory
under contract
W-7405-ENG-48. REFERENCES
Aki, K., Space and time spectra of stationary stochastic waves, with special reference to microtremors, Bull. Earthquake Res. Inst. Univ. Tokyo, 25, 415-457, 1957.
Aki, K., and R. Y. Koyanagi, Deep volcanic tremor and magma ascent mechanism under Kilauea, Hawaii, J. Geophys. Res., 86, 7095-7109, 1981. Aki, K., M. Fehler, and S. Das, Source mechanism of volcanic tremor: Fluid-driven crack models and their application to the 1963 Kilauea eruption, J. Volcanol. Geotherm. Res., œ, 259-287, 1977.
Borcherdt, R. D., J. Fletcher, E. Jensen, G. Maxwell, J. Van Schaack, R. Warrick, E. Cranswick, M. Johnston, and R. McClearn, A General Earthquake Observation System
(GEOS), Bull. Setstool.Soc.Am., 75, 1783-1823,1985. Bouchon, M., Discrete wave number representation of elastic wave fields in three-space dimensions, J. Geophys.Res., 8•, 3609-3614, 1979. Chouet, B., Excitation of a buried magmatic pipe: A seismic source model for volcanic tremor, J. Geophys. Res., 90, 1881-1893, 1985.
Chouet, B., Representation of an extended seismic source in a propagator-based formalism, Bull. Setstool. Soc. Am., 77, 14-27, 1987. Chouet, B., Resonance of a fluid-driven crack: Radiation properties and implications for the source of long-period events and harmonic tremor, J. Geophys. Res., 93, 43734400, 1988.
on Mount
St. Helena
Volcano
and
observations
of long-period seismic events that originate under the volcano, Geophys. Res. Lett., 9, 1017-1020, 1982. Ferrazzini, V., and K. Aki, Preliminary results from a field experiment on volcanic events at Kilauea using an array of digital seismographs, in IAVCEI Proceedings in Volcanology, vol. 3, Volcanic Seisinology, edited by P. Gasparini, R. Scarpa, and K. Aki, pp. 168-189, Springer-Verlag, New York, 1992. Ferrazzini, V., K. Aki, and B. Chouet, Characteristics of seismic waves composing Hawaiian volcanic tremor and gas-piston events observed by a near-source array, J. Geophys. Res., 96, 6199-6209, 1991. Furumoto, M., T. Kunitomo, H. Inoue, I. Yamada, K. Yamaoka, A. Ikami, and Y. Fukao, Twin sources of high-frequency volcanic tremor of Izu-Oshima Volcano, Japan, Geophys. Res. Lett., 17, 25-27, 1990. Goldstein, P., Array measurements of earthquake rupture, Ph.D. thesis, Univ. of Calif. Santa Barbara, 1988. Goldstein, P., and R. J. Archuleta, Array analysis of seismic signals, Geophys. Res. Lett., 1•, 13-16, 1987. Goldstein, P., and R. J. Archuleta, Deterministic frequencywavenumber techniques and direct measurements of rupture propagation during earthquakes using a dense array-theory and methods, J. Geophys. Res., 96, 61736185, 1991a.
Goldstein, P., and R. J. Archuleta, Deterministic frequencywavenumber techniques and direct measurements of rupture propagation during earthquakes using a dense array-data analysis, J. Geophys. Res., 96, 6187-6198, 1991b.
Goldstein, P., W. R. Walter, G. Zandt, Upper mantle structure beneath the Soviet platform from broadband regional recordings of an array of nuclear explosions, J. Geophys. Res., 97, 14,097-14,113, 1992. Gordeev, E. I., Modelling of volcanic tremor wave fields, J. Volcanol. Geotherm. Res., , 51, 145-160, 1992. Gordeev, E. I., V. A. Saltykov, V. I. Sintrain, and V. N. Chebrov, Temporal and spatial characteristics of volcanic tremor wavefields, J. Volcanol. Geotherm. Res., •0, 89-101, 1990.
Jenkins, G. M. and D. G. Watts, Spectral Analysis and Its Applications, pp. 209-257, Holden-Day, Oakland, Calif., 1968.
Kamo, K., T. Furuzawa, and J. Akamatsu, Some natures of volcanic micro-tremors at the Sakurajima Volcano, Bull. Volcanol. Soc. Jpn., œœ,41-58, 1977. Koyanagi, R. Y., W. R. Tanigawa, and J. S. Nakata, Seismicity associated with the eruption, in The Puu Oo Eruption of Kilauea Volcano, Hawaii: Episodes I Through 20, January 3, 1983, Through June 8, 1984, edited by E. W. Wolfe, U.S. Geol. Surv. Prof. Pap., 1•63, 183-235, 1988.
2652
GOLDSTEIN
AND
CHOUET:
SOURCES
Kubotera, A., Volcanic tremors at Aso volcano, in Physical Volcanology,edited by L. Civetta, P. G asparini, G. Luongo, and A. Rapolla, pp. 29-47, Elsevier, New York, 1974.
Malone, S. D., E. T. Endo, C. S. Weaver, and J. W. Ramey, Seismic monitoring for eruption prediction, in The 1980 Eruptions of Mount St. Helens, Washington, edited by P. W. Lipman and D. R. Mullineaux, U.S. Geol. Surv. Prof. Pap., 1œ50,803-813, 1981. McNutt, S. R., Observations and analysis of B-type earthquakes, explosions, and volcanic tremor at Pavlov volcano, Alaska, Bull. Seismol. Soc. Am., 76, 153-175, 1986.
OF VOLCANIC
TREMOR
AT KILAUEA
eruption, Kilauea Volcano, Hawaii, it U.S. Prof. Pap., 1056, 55 pp., 1979.
Geol. Surv.
Shimozuru, D., K. Kamo, and W. Kinoshita, Volcanic tremor of Kilauea Volcano,Hawaii duringJuly-December, 1963, Bull. EarthquakeRes. Inst. TokyoUniv., 44, 10931133, 1966.
Spudich, P., and D. Oppenheimer, Dense seismograph array observationsof earthquake rupture dynamics, in EarthquakeSourceMechanics,Geophys.Monogr.Set., vol. 37, edited by S. Das, J. Boatwright, and C. H. Scholz,pp. 285-296, AGU, Washington, D.C.,
1986.
Watanabe,H., Seismicactivity,in The 1986-1987Eruption of Izu-Oshima Volcano,edited by S. Aramaki, pp. 12-21,
Minakami, T., Fundamental research for predicting volcanic Earthq. Res. Inst. Univ. Tokyo, Tokyo, 1988. eruptions, I, Earthquakes and crustal deformations Wolfe, E. W., C. A. Neal, N. G. Banks, and T. J. Duggan, originating from volcanic activities, Bull. Earthquake Res. Geologicobservationsand chronologyof eruptive events, Inst. Univ. Tokyo, $8, 497-544, 1960. in The Puu Oo Eruption of Kilauea Volcano, Hawaii: Mfiller, G., The reflectivity method: A tutorial, J. Geophys., EpisodesI Through 20, January 3, 1983, Through June 8, 58, 153-174, 1985.
Okada, H., Y. Nishimura, H. Miyamachi, H. Mori, and K. Ishihara, Geophysical significance of the 1988-1989 explosive eruptions of Mt. Tokachi, Hokkaido, Japan, Bull. Volcanol. Soc. Jpn., $5, 175-203, 1990. Omer, G. C., Jr., Volcanic tremor, 2, The theory of volcanic tremor, Bull. Seismol. Soc. Am., 40, 175-194, 1950. Schick, R., G. Lombardo, and G. Patane, Volcanic tremors
1984, edited by E. W. Wolfe, U.S. Geol.Surv. Prof. Pap., 1463, 1-97, 1988.
Yamaoka, K., J. Oikawa, and Y. Ida, An isotropicsource of
volcanic
tremor--observation
with
a
dense
mic network at Izu-Oshima Volcano, Japan, canol. Geotherm. Res., 47, 329-336, 1991.
seis-
J. Vol-
B. Chouet, U.S. GeologicalSurvey,345 Middlefield Road, and shocksassociatedwith eruptionsat Etna (Sicily), Menlo Park, CA 94025. September 1980, Y. Volcanol. Geotherm. Res., 14, 261-
279, 1982. Schmidt, R. O., Multiple emitter location and signal parameter estimation, IEEE Trans. Antennas Propag., AP, $4, 276-280, 1986. Swanson, D. A., W. Duffield, D. Jackson, and D. Peterson, Chronological narrative of the 1969-71 Mauna Ulu
P. Goldstein, Earth SciencesDepartment, L-205, Lawrence Livermore National Laboratory, P.O. Box 808, Livermore, CA, 94550.
(Received May 17, 1993;revisedSeptember 7, 1993; acceptedSeptember13, 1993.)