Kilauea's Summit Research Articles

3‐D Modeling of Irregular Volcanic Sources Using Sparsity‐Promoting Inversions of Geodetic Data and Boundary Element Method

AbstractGeodetic observations of surface deformation associated with volcanic activities can be used to constrain volcanic source parameters and their kinematics. Simple analytical models, such as point and spherical sources, are widely used to model deformation data. The inherent nature of oversimplified model geometries makes them unable to explain fine details of surface deformation. Current nonparametric, geometry‐free inversion approaches resolve the distributed volume change, assuming it varies smoothly in space, which may detect artificial volume change outside magmatic source regions. To obtain a physically meaningful representation of an irregular volcanic source, we devise a new sparsity‐promoting modeling scheme assuming active magma bodies are well‐localized melt accumulations, namely, outliers in the background crust. First, surface deformation data are inverted using a hybridL1‐ andL2‐norm regularization scheme to solve for sparse volume change distributions. Next, a boundary element method is implemented to solve for the displacement discontinuity distribution of the reservoir, which satisfies a uniform pressure boundary condition. The inversion approach is thoroughly validated using benchmark and synthetic tests, of which the results show that source dimension, depth, and shape can be recovered appropriately. We apply this modeling scheme to deformation observed at Kilauea summit for periods of uplift and subsidence leading to and following the 2007 Father's Day event. We find that the magmatic source geometries for these periods are statistically distinct, which may be an indicator that magma is released from isolated compartments due to large differential pressure leading to the rift intrusion.

Very long period conduit oscillations induced by rockfalls at Kilauea Volcano, Hawaii

Eruptive activity at the summit of Kilauea Volcano, Hawaii, beginning in 2010 and continuing to the present time is characterized by transient outgassing bursts accompanied by very long period (VLP) seismic signals triggered by rockfalls from the vent walls impacting a lava lake in a pit within the Halemaumau pit crater. We use raw data recorded with an 11‐station broadband network to model the source mechanism of signals accompanying two large rockfalls on 29 August 2012 and two smaller average rockfalls obtained by stacking over all events with similar waveforms to improve the signal‐to‐noise ratio. To determine the source centroid location and source mechanism, we minimize the residual error between data and synthetics calculated by the finite difference method for a point source embedded in a homogeneous medium that takes topography into account. We apply a new waveform inversion method that accounts for the contributions from both translation and tilt in horizontal seismograms through the use of Green's functions representing the seismometer response to translation and tilt ground motions. This method enables a robust description of the source mechanism over the period range 1–1000 s. The VLP signals associated with the rockfalls originate in a source region ∼1 km below the eastern perimeter of the Halemaumau pit crater. The observed waveforms are well explained by a simple volumetric source with geometry composed of two intersecting cracks including an east striking crack (dike) dipping 80° to the north, intersecting a north striking crack (another dike) dipping 65° to the east. Each rockfall is marked by a similar step‐like inflation trailed by decaying oscillations of the volumetric source, attributed to the efficient coupling at the source centroid location of the pressure and momentum changes induced by the rock mass impacting the top of the lava column. Assuming a simple lumped parameter representation of the shallow magmatic system, the observed pressure and volume variations can be modeled with the following attributes: rockfall volume (200–4500 m3), length of magma column (120–210 m), diameter of pipe connecting the Halemaumau pit crater to the subjacent dike system (6 m), average thickness of the two underlying dikes (3–6 m), and effective magma viscosity (30–210 Pa s). Most rockfalls occur during episodes of sustained deflation of the Kilauea summit. The mass loss rate in the shallow magmatic system is estimated to be 1400–15,000 kg s−1 based on measurements of the temporal variation of VLP period in the two large rockfalls that occurred on 29 August 2012.

Forward Modeling to Assess and Improve Gravity Network Geometry at Kilauea Volcano, Hawai`i

Scientists have been using campaign gravity surveys to monitor volcanic activity at Kilauea's summit for decades, yet we have a poor understanding of the ability of the existing network to resolve sources of magma accumulation with different mass changes and depths. We also do not yet have a fully quantified measure of the relative importance of the stations in the network. This research tests the network using a simple forward modeling approach over a range of likely source volumes and depths. The analysis determines network sensitivity to three different likely source locations, calculates the relative importance of stations in the network, and examines the problem of signal distortion imposed by network geometry. This work finds that the current network is least sensitive to south caldera sources, and investigates the location and number of stations that will resolve this problem most effectively.

Seismicity and deformation induced by magma accumulation at three basaltic volcanoes

We analyzed the evolution of volcano‐tectonic (VT) seismicity and deformation at three basaltic volcanoes (Kilauea, Mauna Loa, Piton de la Fournaise) during phases of magma accumulation. We observed that the VT earthquake activity displays an accelerating evolution at the three studied volcanoes during the time of magma accumulation. At the same times, deformation rates recorded at the summit of Kilauea and Mauna Loa volcanoes were not accelerating but rather tend to decay. To interpret these observations, we propose a physical model describing the evolution of pressure produced by the accumulation of magma into a reservoir. This variation of pressure is then used to force a simple model of damage, where damage episodes are equivalent to earthquakes. This model leads to an exponential increase of the VT activity and to an exponential decay of the deformation rate during accumulation phases. Seismicity and deformation data are well fitted by such an exponential model. The time constant, deduced from the exponential increase of the seismicity, is in agreement with the time constant predicted by the model of magma accumulation. This VT activity can thus be a direct indication of the accumulation of magma at depth, and therefore can be seen as a long‐term precursory phenomenon, at least for the three studied basaltic volcanoes. Unfortunately, it does not allow the prediction of the onset of future eruptions, as no diverging point (i.e., critical time) is present in the model.

Small Explosion From New Vent at Kilauea's Summit

At 0258 Hawaii‐Aleutian Standard Time (HST) on 19 March 2008, a small explosion scattered altered and fresh lithic debris across a 40‐hectare area at the summit of Kilauea volcano. This explosion, the first recorded there since 1924, issued from a vent about 35 meters wide along the east wall of Halema'uma'u Crater. Ballistic fragments—the largest measuring nearly 1 meter across—were propelled upward more than 70 meters onto the Halema'uma'u crater rim. Coarse ash and centimeter‐size lithic debris covered part of Crater Rim Drive, and fine ash was deposited farther than 30 kilometers to the southwest.

Comparative velocity structure of active Hawaiian volcanoes from 3-D onshore–offshore seismic tomography

We present a 3-D P-wave velocity model of the combined subaerial and submarine portions of the southeastern part of the Island of Hawaii, based on first-arrival seismic tomography of marine airgun shots recorded by the onland seismic network. Our model shows that high-velocity materials (6.5–7.0 km/s) lie beneath Kilauea's summit, Koae fault zone, and the upper Southwest Rift Zone (SWRZ) and upper and middle East Rift Zone (ERZ), indicative of magma cumulates within the volcanic edifice. A separate high-velocity body of 6.5–6.9 km/s within Kilauea's lower ERZ and upper Puna Ridge suggests a distinct body of magma cumulates, possibly connected to the summit magma cumulates at depth. The two cumulate bodies within Kilauea's ERZ may have undergone separate ductile flow seaward, influencing the submarine morphology of Kilauea's south flank. Low velocities (5.0–6.3 km/s) seaward of Kilauea's Hilina fault zone, and along Mauna Loa's seaward facing Kao'iki fault zone, are attributed to thick piles of volcaniclastic sediments deposited on the submarine flanks. Loihi seamount shows high-velocity anomalies beneath the summit and along the rift zones, similar to the interpreted magma cumulates below Mauna Loa and Kilauea volcanoes, and a low-velocity anomaly beneath the oceanic crust, probably indicative of melt within the upper mantle. Around Kilauea's submarine flank, a high-velocity anomaly beneath the outer bench suggests the presence of an ancient seamount that may obstruct outward spreading of the flank. Mauna Loa's southeast flank is also marked by a large, anomalously high-velocity feature (7.0–7.4 km/s), interpreted to define an inactive, buried volcanic rift zone, which might provide a new explanation for the westward migration of Mauna Loa's current SWRZ and the growth of Kilauea's SWRZ.

Rapid passage of a small-scale mantle heterogeneity through the melting regions of Kilauea and Mauna Loa Volcanoes

Recent Kilauea and Mauna Loa lavas provide a snapshot of the size, shape, and distribution of compositional heterogeneities within the Hawaiian mantle plume. Here we present a study of the Pb, Sr, and Nd isotope ratios of two suites of young prehistoric lavas from these volcanoes: (1) Kilauea summit lavas erupted from AD 900 to 1400, and (2) 14C-dated Mauna Loa flows erupted from ∼ 2580–140 yr before present (relative to AD 1950). These lavas display systematic isotopic fluctuations, and the Kilauea lavas span the Pb isotopic divide that was previously thought to exist between these two volcanoes. For a brief period from AD 250 to 1400, the 206Pb/ 204Pb and 87Sr/ 86Sr isotope ratios and ε Nd values of Kilauea and Mauna Loa lavas departed from values typical for each volcano (based on historical and other young prehistoric lavas), moved towards an intermediate composition, and subsequently returned to typical values. This is the only known period in the eruptive history of these volcanoes when such a simultaneous convergence of Pb, Sr, and Nd isotope ratios has occurred. The common isotopic composition of lavas erupted from both Kilauea and Mauna Loa during this transient magmatic event was probably caused by the rapid passage of a small-scale compositional heterogeneity through the melting regions of both volcanoes. This heterogeneity is thought to have been either a single body (∼ 35 km long based on the distance between the summits of these volcanoes) or the plume matrix itself (which would be expected to be present beneath both volcanoes). The time scale of this event (centuries) is much shorter than previously noted for variations in the isotopic composition of Hawaiian lavas due to the upwelling of heterogeneities within the plume (thousands to tens of thousands of years). Calculations based on the timing of the isotopic convergence suggest a maximum thickness for the melting region (and thus, the heterogeneity) of ∼ 5–10 km. The small size of the heterogeneity indicates that melt can be extracted from small regions within the Hawaiian plume with minimal subsequent chemical modification (beyond the effects of crystal fractionation). This would be most effective if melt transport in the mantle beneath Hawaiian shield volcanoes occurs mostly in chemically isolated channels.

Deep magma transport at Kilauea volcano, Hawaii

The shallow part of Kilauea's magma system is conceptually well-understood. Long-period and short-period (brittle-failure) earthquake swarms outline a near-vertical magma transport path beneath Kilauea's summit to 20 km depth. A gravity high centered above the magma transport path demonstrates that Kilauea's shallow magma system, established early in the volcano's history, has remained fixed in place. Low seismicity at 4–7 km outlines a storage region from which magma is supplied for eruptions and intrusions. Brittle-failure earthquake swarms shallower than 5 km beneath the rift zones accompany dike emplacement. Sparse earthquakes extend to a decollement at 10–12 km along which the south flank of Kilauea is sliding seaward. This zone below 5 km can sustain aseismic magma transport, consistent with recent tomographic studies. Long-period earthquake clusters deeper than 40 km occur parallel to and offshore of Kilauea's south coast, defining the deepest seismic response to magma transport from the Hawaiian hot spot. A path connecting the shallow and deep long-period earthquakes is defined by mainshock–aftershock locations of brittle-failure earthquakes unique to Kilauea whose hypocenters are deeper than 25 km with magnitudes from 4.4 to 5.2. Separation of deep and shallow long-period clusters occurs as the shallow plumbing moves with the volcanic edifice, while the deep plumbing is centered over the hotspot. Recent GPS data agrees with the volcano-propagation vector from Kauai to Maui, suggesting that Pacific plate motion, azimuth 293.5° and rate of 7.4 cm/yr, has been constant over Kilauea's lifetime. However, volcano propagation on the island of Hawaii, azimuth 325°, rate 13 cm/yr, requires southwesterly migration of the locus of melting within the broad hotspot. Deep, long-period earthquakes lie west of the extrapolated position of Kilauea backward in time along a plate-motion vector, requiring southwesterly migration of Kilauea's magma source. Assumed ages of 0.4 my for Kilauea and 0.8 my for Mauna Loa are consistent with this model. Younger ages would apply if Kilauea began its growth south of the locus of maximum melting, as is true for Loihi seamount. We conclude that Kilauea is fed from below the eastern end of the zone of deep long-period earthquakes. Magma transport is vertical below 30 km, then sub-horizontal, following the oceanic mantle boundary separating plagioclase- and spinel-peridotite, then near-vertical beneath Kilauea's summit. The migration of the melting region within the hotspot and Kilauea's sampling of different sources within the melting region can explain (1) the long-term geochemical separation of Kilauea from neighboring volcanoes Mauna Loa and Loihi, and (2) the short-term changes in trace-element and isotope signatures within Kilauea.

Location of long‐period events below Kilauea Volcano using seismic amplitudes and accurate relative relocation

We present methods for improving the location of long‐period (LP) events, deep and shallow, recorded below Kilauea Volcano by the permanent seismic network. LP events might be of particular interest to understanding eruptive processes as their source mechanism is assumed to directly involve fluid transport. However, it is usually difficult or impossible to locate their source using traditional arrival time methods because of emergent wave arrivals. At Kilauea, similar LP waveform signatures suggest the existence of LP multiplets. The waveform similarity suggests spatially close sources, while catalog solutions using arrival time estimates are widely scattered beneath Kilauea's summit caldera. In order to improve estimates of absolute LP location, we use the distribution of seismic amplitudes corrected for station site effects. The decay of the amplitude as a function of hypocentral distance is used for inferring LP location. In a second stage, we use the similarity of the events to calculate their relative positions. The analysis of the entire LP seismicity recorded between January 1997 and December 1999 suggests that a very large part of the LP event population, both deep and shallow, is generated by a small number of compact sources. Deep events are systematically composed of a weak high‐frequency onset followed by a low‐frequency wave train. Aligning the low‐frequency wave trains does not lead to aligning the onsets indicating the two parts of the signal are dissociated. This observation favors an interpretation in terms of triggering and resonance of a magmatic conduit. Instead of defining fault planes, the precise relocation of similar LP events, based on the alignment of the high‐energy low‐frequency wave trains, defines limited size volumes.

The 12 September 1999 Upper East Rift Zone dike intrusion at Kilauea Volcano, Hawaii

Deformation associated with an earthquake swarm on 12 September 1999 in the Upper East Rift Zone of Kilauea Volcano was recorded by continuous GPS receivers and by borehole tiltmeters. Analyses of campaign GPS, leveling data, and interferometric synthetic aperture radar (InSAR) data from the ERS‐2 satellite also reveal significant deformation from the swarm. We interpret the swarm as resulting from a dike intrusion and model the deformation field using a constant pressure dike source. Nonlinear inversion was used to find the model that best fits the data. The optimal dike is located beneath and slightly to the west of Mauna Ulu, dips steeply toward the south, and strikes nearly east‐west. It is approximately 3 by 2 km across and was driven by a pressure of ∼15 MPa. The total volume of the dike was 3.3 × 106 m3. Tilt data indicate a west to east propagation direction. Lack of premonitory inflation of Kilauea's summit suggests a passive intrusion; that is, the immediate cause of the intrusion was probably tensile failure in the shallow crust of the Upper East Rift Zone brought about by persistent deep rifting and by continued seaward sliding of Kilauea's south flank.

Volcano monitoring using the Global Positioning System: Filtering strategies

Permanent Global Positioning System (GPS) networks are routinely used for producing improved orbits and monitoring secular tectonic deformation. For these applications, data are transferred to an analysis center each day and routinely processed in 24‐hour segments. To use GPS for monitoring volcanic events, which may last only a few hours, real‐time or near real‐time data processing and subdaily position estimates are valuable. Strategies have been researched for obtaining station coordinates every 15 min using a Kalman filter; these strategies have been tested on data collected by a GPS network on Kilauea Volcano. Data from this network are tracked continuously, recorded every 30 s, and telemetered hourly to the Hawaiian Volcano Observatory. A white noise model is heavily impacted by data outages and poor satellite geometry, but a properly constrained random walk model fits the data well. Using a borehole tiltmeter at Kilauea's summit as ground‐truth, solutions using different random walk constraints were compared. This study indicates that signals on the order of 5 mm/h are resolvable using a random walk standard deviation of 0.45 cm/√h. Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1–6 hours.

Constraints on dike propagation from continuous GPS measurements

The January 1997 East Rift Zone eruption on Kilauea volcano, Hawaii, occurred within a network of continuous Global Positioning System (GPS) receivers. The GPS measurements reveal the temporal history of deformation during dike intrusion, beginning ∼8 hours prior to the onset of the eruption. The dike volume as a function of time, estimated from the GPS data using elastic Green's functions for a homogeneous half‐space, shows that only two thirds of the final dike volume accumulated prior to the eruption and the rate of volume change decreased with time. These observations are inconsistent with simple models of dike propagation, which predict accelerating dike volume up to the time of the eruption and little or no change thereafter. Deflationary tilt changes at Kilauea summit mirror the inferred dike volume history, suggesting that the rate of dike propagation is limited by flow of magma into the dike. A simple, lumped parameter model of a coupled dike magma chamber system shows that the tendency for a dike to end in an eruption (rather than intrusion) is favored by high initial dike pressures, compressional stress states, large, compressible magma reservoirs, and highly conductive conduits linking the dike and source reservoirs. Comparison of model predictions to the observed dike volume history, the ratio of erupted to intruded magma, and the deflationary history of the summit magma chamber suggest that most of the magma supplied to the growing dike came from sources near to the eruption through highly conductive conduits. Interpretation is complicated by the presence of multiple source reservoirs, magma vesiculation and cooling, as well as spatial variations in dike‐normal stress. Reinflation of the summit magma chamber following the eruption was measured by GPS and accompanied a rise in the level of the Pu'u O'o lava lake. For a spheroidal chamber these data imply a summit magma chamber volume of ∼20 km3, consistent with recent estimates from seismic tomography. Continuous deformation measurements can be used to image the spatiotemporal evolution of propagating dikes and to reveal quantitative information about the volcanic plumbing systems.

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.

A Rapid Fluctuation in the Mantle Source and Melting History of Kilauea Volcano Inferred from the Geochemistry of its Historical Summit Lavas (1790-1982)

The geochemical variations of Kilauea's historical summit lavas (1790–1982) document the short-term magmatic evolution of one of the Earth's most active volcanoes. Most of these lavas are thought to have erupted directly from the shallow (2–4 km deep) magma reservoir that underlies the volcano's summit region. This paper details a remarkable variation of lava chemistry that spans nearly the entire known compositional range of the volcano in only 200 years. The Pb, Sr, and Nd isotope and incompatible trace element ratios (e.g. La/Yb or Nb/Y) of the lavas vary systematically over time with an abrupt reversal after 1924. This rapid geochemical fluctuation records the temporal changes in the parental magma composition delivered to Kilauea's summit reservoir since 1790. The isotope and incompatible trace element ratio systematics suggest that the source region of historical Kilauea magma is both isotopically and chemically heterogeneous. These source variations can be explained by the melting of small-scale heterogeneities within the Hawaiian mantle plume. Model calculations suggest that the degree of partial melting decreased from the early 19th century until the mid-20th century, which correlates with a lower eruption rate (and presumably a lower magma supply rate) for the volcano between 1840 and 1959. This interval of declining output from the Hawaiian plume culminated with an explosive summit eruption in 1924 and the longest quiescent period in Kilauea's historical record (1934–1952). Lavas erupted just after 1924 are geochemically anomalous and may have been contaminated by the assimilation of country rock into the volcano's magma reservoir during the explosions. Subsequently, the inferred degree of partial melting and the volcano's eruption rate have increased, with the highest values since the early 19th century observed during the current Puu Oo rift zone eruption.

The size and shape of Kilauea Volcano's summit magma storage reservoir: a geochemical probe

One of the most important components of the magmatic plumbing system of Kilauea Volcano is the shallow (2–4 km deep) magma storage reservoir that underlies the volcano's summit region. Nevertheless, the geometry (shape and size) of Kilauea's summit reservoir is controversial. Two fundamentally different models for the reservoir's shape have been proposed based on geophysical observations: a plexus of dikes and sills versus a single, `spherical' magma body. Furthermore, the size of the reservoir is poorly constrained with estimates ranging widely from 0.08 to 40 km 3. In this study, we use the temporal variations of Pb, Sr, and Nd isotope and incompatible trace element (e.g., La/Yb and Nb/Y) ratios of Kilauea's historical summit lavas (1790–1982) to probe the geometry of the volcano's summit reservoir. These lavas preserve a nearly continuous, 200-year record of the changes in the composition of the parental magma supplied to the volcano. The systematic temporal variations in lava chemistry at Kilauea since the early 19th century suggest that the shape of the volcano's summit reservoir is relatively simple. Residence time analysis of these rapid geochemical fluctuations indicates that the volume of magma in Kilauea's summit reservoir is only ∼2–3 km 3, which is smaller than most geophysical estimates (2–40 km 3). This discrepancy can be explained if the volume calculated from lava chemistry represents the hotter, molten core of the reservoir in which magma mixing occurs, whereas the volumes estimated from geophysical data also include portions of the reservoir's outer crystal-mush zone and a hot, ductile region that surrounds the reservoir. Although our volume estimate is small, the amount of magma stored within Kilauea's summit reservoir since the early 19th century is an order of magnitude larger than the magma body supplying Piton de la Fournaise Volcano, another frequently active ocean-island volcano.

Stabilization of volcanic flanks by dike intrusion: an example from Kilauea

Dike propagation and dilation increases the compression of adjacent rocks. On volcanoes, especially oceanic shields, dikes are accordingly thought to be structurally destabilizing. As compression is incremented, volcanic flanks are driven outward or downslope and thus increase their susceptibility to destructive earthquakes and giant landslides. We show, however, that the 2-m-thick dike emplaced along the east rift zone of Kilauea in 1983 actually stabilized that volcano's flank. Specifically, production of flank earthquakes dropped more than twofold after 1983 as maximum downslope motion slowed to 6 cm·year–1 from approximately 40 cm·year–1 during 1980–1982. As much as 65 cm of deflationary subsidence above Kilauea's summit and upper rift zones accompanied the dike intrusion. According to recent estimates, this deflation corresponds to a reduction in magma-reservoir pressure of approximately 4 MPa, probably about as much as the driving pressure of the 1983 dike. The volume of the dike, approximately 0.10–0.15 km3, is orders of magnitude less than the estimated 200- to 250-km3 volume of Kilauea's reservoir of magma and nearby hot, mushy rock. Thus, deflation of that reservoir reduces the compressional load on the flank over a much larger area than intrusion of the dike adds to it, particularly at the dominant depth of seismicity, 8–9 km. A Coulomb block model for flank motion during intervals between major earthquakes requires the low-angle fault beneath Kilauea's flank to exhibit slip weakening, conducive to earthquake instability. Accordingly, the triggering mechanism of destructive earthquakes, several of which have struck Hawaii during the past 150 years, need not require stresses accumulated by dike intrusions.

Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea volcano, Hawaii

We use data from broadband seismometers deployed around the summit of Kilauea Volcano to quantify the mechanism associated with a transient in the flow of magma feeding the east rift eruption of the volcano. The transient is marked by rapid inflation of the Kilauea summit peaking at 22 μrad 4.5 hours after the event onset, followed by slow deflation over a period of 3 days. Superimposed on the summit inflation is a series of sawtooth displacement pulses, each characterized by a sudden drop in amplitude lasting 5–10 s followed by an exponential recovery lasting 1–3 min. The sawtooth waveforms display almost identical shapes, suggesting a process involving the repeated activation of a fixed source. The particle motion associated with each sawtooth is almost linear, and its major swing shows compressional motion at all stations. Analyses of semblance and particle motion are consistent with a point source located 1 km beneath the northeast edge of the Halemaumau pit crater. To estimate the source mechanism, we apply a moment tensor inversion to the waveform data, assuming a point source embedded in a homogeneous half‐space with compressional and shear wave velocities representative of the average medium properties at shallow depth under Kilauea. Synthetic waveforms are constructed by a superposition of impulse responses for six moment tensor components and three single force components. The origin times of individual impulses are distributed along the time axis at appropriately small, equal intervals, and their amplitudes are determined by least squares. In this inversion, the source time functions of the six tensor and three force components are determined simultaneously. We confirm the accuracy of the inversion method through a series of numerical tests. The results from the inversion show that the waveform data are well explained by a pulsating transport mechanism operating on a subhorizontal crack linking the summit reservoir to the east rift of Kilauea. The crack acts like a buffer in which a batch of fluid (magma and/or gas) accumulates over a period of 1–3 min before being rapidly injected into a larger reservoir (possibly the east rift) over a timescale of 5–10 s. The seismic moment and volume change associated with a typical batch of fluid are approximately 1014 N m and 3000 m3, respectively. Our results also point to the existence of a single force component with amplitude of 109 N, which may be explained as the drag force generated by the flow of viscous magma through a narrow constriction in the flow path. The total volume of magma associated with the 4.5‐hour‐long activation of the pulsating source is roughly 500,000 m3 in good agreement with the integrated volume flow rate of magma estimated near the eruptive site.

Effects of eruption and lava drainback on the H 2 O contents of basaltic magmas at Kilauea Volcano

Lava drainback has been observed during many eruptions at Kilauea Volcano: magma erupts, degasses in lava fountains, collects in surface ponds, and then drains back beneath the surface. Time series data for melt inclusions from the 1959 Kilauea Iki picrite provide important evidence concerning the effects of drainback on the H2O contents of basaltic magmas at Kilauea. Melt inclusions in olivine from the first eruptive episode, before any drainback occurred, have an average H2O content of 0.7±0.2 wt.%. In contrast, many inclusions from the later episodes, erupted after substantial amounts of surface degassed lava had drained back down the vent, have H2O contents that are much lower (≥0.24 wt.% H2O). Water contents in melt inclusions from magmas erupted at Pu'u 'O'o on the east rift zone vary from 0.39–0.51 wt.% H2O in tephra from high fountains to 0.10–0.28 wt.% H2O in spatter from low fountains. The low H2O contents of many melt inclusions from Pu'u 'O'o and post-drainback episodes of Kilauea Iki reveal that prior to crystallization of the enclosing olivine host, the melts must have exsolved H2O at pressures substantially less than those in Kilauea's summit magma reservoir. Such low-pressure H2O exsolution probably occurred as surface degassed magma was recycled by drainback and mixing with less degassed magma at depth. Recognition of the effects of low-pressure degassing and drainback leads to an estimate of 0.7 wt.% H2O for differentiated tholeiitic magma in Kilauea's summit magma storage reservoir. Data for MgO-rich submarine glasses (Clague et al. 1995) and melt inclusions from Kilauea Iki demonstrate that primary Kilauean tholeiitic magma has an H2O/K2O mass ratio of ∼1.3. At transition zone and upper mantle depths in the Hawaiian plume source, H2O probably resides partly in a small amount of hydrous silicate melt.

A dynamic balance between magma supply and eruption rate at Kilauea volcano, Hawaii

The dynamic balance between magma supply and vent output at Kilauea volcano is used to estimate both the volume of magma stored within Kilauea volcano and its magma supply rate. Throughout most of 1991 a linear decline in volume flux from the Kupaianaha vent on Kilauea's east rift zone was associated with a parabolic variation in the elevation of Kilauea's summit as vent output initially exceeded then lagged behind the magma supply to the volcano. The correspondence between summit elevation and tilt established with over 30 years of data provided daily estimates of summit elevation in terms of summit tilt. The minimum in the parabolic variation in summit tilt and elevation (or zero elevation change) occurs when the magma supply to the reservoir from below the volcano equals the magma output from the reservoir to the surface, so that the magma supply rate is given by vent flux on that day. The measurements of vent flux and tilt establish that the magma supply rate to Kilauea volcano on June 19, 1991, was 217,000±10,000 m3/d (or 0.079±0.004 km3/yr). This is close to the average eruptive rate of 0.08 km3/yr between 1958 and 1984. In addition, the predictable response of summit elevation and tilt to each east rift zone eruption near Puu Oo since 1983 shows that summit deformation is also a measure of magma reservoir pressure. Given this, the correlation between the elevation of the Puu Oo lava lake (4 km uprift of Kupaianaha and 18 km from the summit) and summit tilt provides an estimate for magma pressure changes corresponding to summit tilt changes. The ratio of the change in volume to the change in reservoir pressure (dV/dP) during vent activity may be determined by dividing the ratio of volume erupted to change in summit tilt (dV/dtilt) by the ratio of pressure change to change in summit tilt (dP/dtilt). This measure of dV/dP, when combined with laboratory measurements of the bulk modulus of tholeitic melt, provides an estimate of 240±50 km3 for the volume of Kilauea's magma reservoir. This estimate is much larger than traditional estimates but consistent with seismic tomographic imaging and geophysical modeling of Kilauea's magma system.

Fluid flow and water‐rock interaction in the East Rift Zone of Kilauea Volcano, Hawaii

The East Rift Zone of Kilauea Volcano in Hawaii represents a major area of geothermal activity. Fluid inclusion and stable isotope analyses of secondary hydrothermal minerals in core samples from three scientific observation holes (SOH) drilled into the rift zone indicate that the geothermal system is dominated by meteoric waters to depths of as much as 1500 m below sea level. Calculated δ18O and δD values for fluids on the north side of the rift zone indicate that the deep meteoric fluids may be derived from precipitation on the upper slopes of Mauna Loa Volcano. In the interior of the rift zone, recharge is dominated by seawater mixed with local meteoric water. Water/rock ratios in the rift area are approximately 2, but strongly 18O‐enriched fluids in the deeper parts of the SOH‐2 and SOH‐4 drill holes (on the north side of the rift) indicate that the fluids underwent extensive interaction with rocks prior to reaching this part of the rift zone. Marine carbonates at the subaerial to submarine transition (between 1700 and 1780 m depth) in SOH‐4 have not fully equilibrated with the fluids, suggesting that the onset of hydrothermal activity in this area was relatively recent (<2000 years). This may represent increased volcanic activity along the rift after the end of the Ai La'au phase of eruptive activity at the Kilauea summit approximately 1000 years ago, or it may reflect progressive evolution of the hydrothermal system in response to southward migration of intrusive activity within the rift.