Search

In the summer of 2018, Kīlauea Volcano underwent one of its most significant eruptions in the past few hundred years. The volcano’s summit and East Rift Zone magma system partially drained, resulting in a series of occasionally explosive partial caldera collapses, and widespread lava flows in the lower East Rift Zone. The Hawaiian Volcano Observatory (HVO) operates a robust permanent monitoring network of about 250 stations, recording a variety of real-time data streams: seismic (short-period, broadband, strong-motion), infrasound, Global Navigation Satellite Systems (GNSS), tilt, camera, laser rangefinder, and gas geochemistry. During the eruption, HVO staff quickly established 35 new temporary monitoring stations, to better constrain evolving volcanic hazards. The partial collapses of the caldera threatened to disrupt important telemetry links in the HVO monitoring network, and a major effort was undertaken in the midst of the eruption crisis to reroute radio telemetry and maintain continuity of data flow. In the process, a new data center was established in Hilo, to mitigate a long-standing potential single point of failure at the HVO facility. Over the course of the eruption from May through August, lava, ashfall, wildfire, and cliff collapse destroyed or disabled 36 stations. Thousands of earthquakes damaged the main HVO facility at Uēkahuna Bluff, causing staff to evacuate the building and relocate observatory operations in the midst of the eruption response, adding more complexity to the response effort. Throughout these events, the HVO team maintained the monitoring network, provided timely information to the public and emergency managers, and collected valuable scientific data to better understand Kīlauea Volcano.

Multiple-aperture SAR interferometry (MAI) has demonstrated outstanding measurement accuracy of along-track displacement when compared to pixel-offset-tracking methods; however, measuring slow-moving (cm/year) surface displacement remains a challenge. Stacking of multi-temporal observations is a potential approach to reducing noise and increasing measurement accuracy, but it is difficult to achieve a significant improvement by applying traditional stacking methods to multi-temporal MAI interferograms. This paper proposes an efficient MAI stacking method, where multi-temporal forward- and backward-looking residual interferograms are individually stacked before the MAI interferogram is generated. We tested the performance of this method using ENVISAT data from Kīlauea Volcano, Hawai‘i, where displacement on the order of several centimeters per year is common. By comparing results from the proposed stacking methods with displacements from GPS data, we documented measurement accuracies of about 1.03 and 1.07 cm/year for the descending and ascending tracks, respectively—an improvement of about a factor of two when compared with that from the conventional stacking approach. Three-dimensional surface-displacement maps can be constructed by combining stacked InSAR and MAI observations, which will contribute to a better understanding of a variety of geological phenomena.

We present here a small baseline subset (SBAS) algorithm to extract both transient and secular ground deformations on the order of millimeters in the presence of tropospheric noise on the order of centimeters, when the transient is of short duration and known time and the background deformation is smooth in time. We applied this algorithm to study the 2010 slow slip event as well as the secular motion of Kilauea's south flank using 49 TerraSAR-X images. We also estimate the tropospheric delay variation relative to a given reference pixel using an interferometric synthetic aperture radar (InSAR) SBAS approach. We compare the InSAR SBAS solution for both ground deformation and tropospheric delays with existing GPS measurements and confirm that the ground deformation signal and tropospheric noise in InSAR data are successfully separated. We observe that the coastal region on the south side of the Hilina Pali moves at a higher background rate than the region north side of the Pali. We also conclude that the 2010 slow slip event (SSE) displacement is mainly horizontal and the maximum magnitude of the 2010 SSE vertical component is less than 5 mm.

[1] The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep decollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. Surface dike opening ceased following the eruption, but subsurface opening in the dike continued into May 2011. Dike volumes increased from 15, to 16, to 21 million cubic meters (MCM) after the first day, eruption end, and 2 months following, respectively. Dike shape is distinctive, with a main limb plunging from the surface to 2–3 km depth in the up-rift direction toward Kīlauea's summit, and a lesser projection extending in the down-rift direction toward Pu`u `Ō`ō at 2 km depth. Volume losses beneath Kīlauea summit (1.7 MCM) and Pu`u `Ō`ō (5.6 MCM) crater, relative to dike plus erupted volume (18.3 MCM), yield a dike to source volume ratio of 2.5 that is in the range expected for compressible magma without requiring additional sources. Inflation of Kīlauea's summit in the months before the March 2011 eruption suggests that the Kamoamoa eruption resulted from overpressure of the volcano's magmatic system.

Over the past decade, Kīlauea and Mauna Loa—adjacent volcanoes in Hawai‘i—have inflated and deflated in tandem, yet their shallow volcanic plumbing systems are separate. Numerical modelling of the volcanoes shows that dynamic stress transfer by asthenospheric pore pressure is a viable mechanism for volcano coupling in Hawai‘i.

On January 29±30, 1997, prolonged steady-state effusion of lava from Pu'u'O'o was briefly disrupted by shallow extension beneath Napau Crater, 1±4 km uprift of the active Kilauea vent. A 23-h-long eruption (episode 54) ensued from fissures that were overlapping or en echelon with eruptive fissures formed during episode 1 in 1983 and those of earlier rift zone eruptions in 1963 and 1968. Combined geophysical and petrologic data for the 1994±1999 eruptive interval, including episode 54, reveal a variety of shallow magmatic conditions that persist in association with prolonged rift zone eruption. Near-vent lava samples document a significant range in composition, temperature and crystallinity of pre-eruptive magma. As supported by phenocryst±liquid relations and Kilauea mineral thermometers established herein, the rift zone extension that led to episode 54 resulted in mixture of near-cotectic magma with discrete magma bodies cooled to 1100 C. Mixing models indicate that magmas isolated beneath Napau Crater since 1963 and 1968 constituted 32±65% of the hybrid mixtures erupted during episode 54. Geophysical measurements support passive displacement of open-system magma along the active east rift conduit into closed-system rift-reservoirs along a shallow zone of extension. Geophysical and petrologic data for early episode 55 document the gradual flushing of episode 54 related magma during magmatic recharge of the edifice.

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.

A continuous Global Positioning System (GPS) network on Kilauea Volcano captured the most recent fissure eruption in Kilauea's East Rift Zone (ERZ) in unprecedented spatial and temporal detail. The short eruption drained the lava pond at Pueu Oeo, leading to a two month long pause in its on-going eruption. Models of the GPS data indicate that the intrusion's bottom edge extended to only 2.4 km. Continuous GPS data reveal rift opening 8 hours prior to the eruption. Absence of precursory summit inflation rules out magma storage overpressurization as the eruption's cause. We infer that stresses in the shallow rift created by the continued deep rift dilation and slip on the south flank decollement caused the rift intrusion.

Campaign Global Positioning System (GPS) measurements from 1990 to 1996 are used to calculate surface displacement rates on Kilauea Volcano, Hawaii. The GPS data show that the south flank of the volcano, which has generated several large earthquakes in the past 3 decades, is displacing at up to ;8 cm/yr to the south-southeast. The summit and rift zones are subsiding, with maximum subsidence rates of ;8 cm/yr observed a few kilometers south of the summit caldera. Elastic dislocation modeling of the GPS data suggests that the active sources of deformation include deep rift opening along the upper east and east rift zone, fault slip along a subhorizontal fault near the base of the volcano, and deflation near the summit caldera. A nonlinear optimization algorithm was used to explore the parameter space and to find the best fitting source geometry. There is a broad range of model geometries that fit the data reasonably well. However, certain models can be ruled out, including those that have shallow rift opening or shallow fault slip. Some offshore, aseismic slip on a fault plane that dips between 258 north- northwest and 88 south-southeast is required. Best fitting slip and rift opening rates are 23-28 cm/yr, although rates as low as 10 cm/yr are permitted by the data.

The rift system traversing about 80 km of the subaerial surface of Kilauea volcano has extended continuously since the M 7.2 flank earthquake of November 1975. Widening across the summit has amounted to more than 250 cm, decelerating after 1975 from about 25 to 4 cm yr−1 since 1983. Concurrently, the summit has subsided more than 200 cm, even as the adjacent south flank has risen more than 50 cm. The axes of the upper zones, about 10 km from the summit, subsided before 1983 at average rates of 9 and 4 cm yr−1, respectively, and at rates of 4 and 3 cm yr−1 since. The middle southwest rift zone is also subsiding and, at the other end of Kilauea's subaerial rift system, subsidence along the lower east rift zone has averaged 1–2 cm yr−1. Deformation of Kilauea's south flank has been continuous, although subject as well to displacements caused by major rift zone seismic swarms. Whereas horizontal strains across the subaerial south flank seem to have been generally compressive after 1975, they have been extensional since about 1980 or 1981, interrupted only by the east rift zone dike intrusion of 1983. Because the magnitudes of these contractions and extensions are much less than the extension across the rift system, the subaerial south flank is apparently sliding seaward on its basal decollement more than it is accumulating horizontal strains within the overlying volcanic pile. Kilauea suffers from gravitational spreading made even more unstable by accumulation of magma along the rift system at depths in excess of about 4–5 km in the presence hot rock incapable of withstanding deviatoric stresses. This seismicly quiescent zone decouples the south flank from the rest of Hawaii's volcanic edifice; the rift zones at lesser depths exhibit a more brittle and, therefore, sporadic extensional behavior. Judging from the modern extension record of the summit, which both predates the M 7.2 earthquake of 1975 and has outlived its 10-year period of aftershocks, Kilauea will continue to spread along its rift system as its south flank slips seaward to accommodate the accretion of magma and its relatively dense olivine-rich differentiate.

Following the amalgamation of a collage of pre-Neogene terranes largely by strike-slip and convergence mechanisms to form the Philippine islands, volcanic chains, related to oppositely dipping subduction zones, developed along the eastern and western margins of the archipelago. There is ample field evidence that this volcanic activity, predominantly calc-alkaline in chemical character, had commenced by the Oligocene. Volcanoes resulting from subduction along the Manila-Negros trench in the west (e.g. Taal, Laguna de Bay and Arayat) form a high-angle linear array, trending away from the MORE field on Pb-isotopic covariation diagrams; have the highest Sr- and lowest Nd-isotopic compositions, of the two chains (but nevertheless plotting above bulk earth on the 87 Sr/ 86 Sr versus 143 Nd/ 144 Nd covariation diagram); and exhibit Sm/Nd and Rb/Sr values that are lower and higher, respectively, than the estimated values for bulk earth. While the Sm/Nd and Rb/Sr characteristics are common to both chains, volcanoes associated with the Philippine-East Luzon trench have Pb-isotopic compositions that fall in the Indian Ocean MORB field and that require time-integrated evolution in a high Th/U environment. They also have higher Nd- and lower Sr-isotopic ratios. The source materials of Philippine volcanoes, therefore, have undergone varied recent enrichments in LILE, as indicated by the decoupling of isotopic and elemental ratios. These enrichments, particularly for the western volcanoes, cannot be entirely due to small degrees of partial melting in the mantle wedge, considering that they were accompanied by elevations in radiogenic Pb. Elevated Pb ratios are best explained by the introduction of subducted, continentally derived sediments. The sedimentary component in the western volcanoes is probably the South China Sea sediments derived largely from Eurasia. That this component is not available in the Philippine-East Luzon trench is reflected by the fact that the eastern volcanoes have higher Nd- and lower Sr-isotopic ratios as well as less radiogenic common Pb.

[1] Rapid changes in ground tilt and GPS positions on Kīlauea Volcano, Hawai'i, are interpreted as resulting from a shallow, two-segment dike intrusion into the east rift zone that began at 1217 UTC (0217 HST) on 17 June 2007 and lasted almost 3 days. As a result of the intrusion, a very small volume of basalt (about 1500 m3) erupted on 19 June. Northward tilt at a coastal tiltmeter, subsidence of south flank GPS sites, southeastward displacements at southwestern flank GPS sites, and a swarm of flank earthquakes suggest that a slow slip event occurred on the decollement beneath Kīlauea's south flank concurrent with the rift intrusion. We use 4 min GPS positions that include estimates of time-dependent tropospheric gradients and ground tilt data to study the spatial and temporal relationships between the two inferred shallow, steeply dipping dike segments extending from the surface to about 2 km depth and decollement slip at 8 km depth. We invert for the temporal evolution of distributed dike opening and decollement slip in independent inversions at each time step using a nonnegative least squares algorithm. On the basis of these inversions, the intrusion occurred in two stages that correspond spatially and temporally with concentrated rift zone seismicity. The dike opening began on the western of the two segments before jumping to the eastern segment, where the majority of opening accumulated. Dike opening preceded the start of decollement slip at an 84% confidence level; the latter is indicated by the onset of northward tilt of a coastal tiltmeter. Displacements at southwest flank GPS sites began about 18 h later and are interpreted as resulting from slow slip on the southwestern flank. Additional constraints on the evolution of the intrusion and decollement slip come from inversion of an Envisat interferogram that spans the intrusion until 0822 UTC on 18 June 2007, combined with GPS and tilt data. This inversion shows that up to 0822 UTC on 18 June, decollement slip is only required in a limited region offshore of Ka'ena Point. A similar inversion of the complete event, which includes GPS and tilt data up to 21 June and a second Envisat interferogram spanning the complete intrusion until 21 June, shows decollement slip spread westward across the south flank. This may suggest westward migration of the decollement slip as the event progressed.

Kilauea erupted almost continuously from January 1983 through 1991. Although the summit began subsiding during the rift zone dike intrusion that initiated this eruption, remarkably steady ground surface motions began in late 1983 after a magnitude 6.6 earthquake beneath the slopes of nearby Mauna Loa volcano and continued until the onset of brief upper east rift zone earthquake swarms in late 1990. During these 7 years the summit and upper rift zones subsided up to 10–11 and 4–8 cm yr−1, respectively, and summit baselines contracted up to 6 cm yr−1. Baselines directed northward from the summit to stations on Mauna Loa extended at rates up to 7 cm yr−1, and a baseline from south of the summit to Mauna Loa extended 4 cm yr−1. Much of this extension is inconsistent with deformation caused solely by summit magma reservoir collapse and more likely reflects rifting as the south flank of the volcano moved seaward from the summit and rift zones. Farther from the summit, baselines crossing the south flank extended up to 2 cm yr−1, and a south flank tide gauge rose 2 cm yr−1; the lower east rift zone, 40–50 km from the summit, subsided about 2 cm yr−1. Motion on Kilauea, then, is broadly consistent with slip along low-angle south flank faults, generating subsidence that is focused at the summit and along the rift system behind the faulting and uplift along the coastal south flank ahead of it. Dislocation models that combine these elements show that much of Kilauea's edifice migrated seaward, producing ground surface motions along the south flank of up to about 6 cm yr−1. The magnitude 6.1 earthquake of 1989 punctuated these motions along the eastern south flank, producing more than 25 cm of seaward displacement and, 15 km east of the epicenter, up to 24 cm of subsidence south of the lower east rift zone. Unlike the magnitude 7.2 south flank earthquake of 1975, the 1989 event was preceded neither by summit magma reservoir inflation nor by rift zone dike intrusions and accompanying compression of the south flank. Deformation was probably caused by the weight of the volcanic overburden and by ongoing dilation and slip within the rift system.

[1] A series of complex events at Kīlauea Volcano, Hawaii, 17 June to 19 June 2007, began with an intrusion in the upper east rift zone (ERZ) and culminated with a small eruption (1500 m3). Surface deformation due to the intrusion was recorded in unprecedented detail by Global Positioning System (GPS) and tilt networks as well as interferometric synthetic aperture radar (InSAR) data acquired by the ENVISAT and ALOS satellites. A joint nonlinear inversion of GPS, tilt, and InSAR data yields a deflationary source beneath the summit caldera and an ENE-striking uniform-opening dislocation with ∼2 m opening, a dip of ∼80° to the south, and extending from the surface to ∼2 km depth. This simple model reasonably fits the overall pattern of deformation but significantly misfits data near the western end of an inferred dike-like source. Three more complex dike models are tested that allow for distributed opening including (1) a dike that follows the surface trace of the active rift zone, (2) a dike that follows the symmetry axis of InSAR deformation, and (3) two en echelon dike segments beneath mapped surface cracks and newly formed steaming areas. The en echelon dike model best fits near-field GPS and tilt data. Maximum opening of 2.4 m occurred on the eastern segment beneath the eruptive vent. Although this model represents the best fit to the ERZ data, it still fails to explain data from a coastal tiltmeter and GPS sites on Kīlauea's southwestern flank. The southwest flank GPS sites and the coastal tiltmeter exhibit deformation consistent with observations of previous slow slip events beneath Kīlauea's south flank, but inconsistent with observations of previous intrusions. Slow slip events at Kīlauea and elsewhere are thought to occur in a transition zone between locked and stably sliding zones of a fault. An inversion including slip on a basal decollement improves fit to these data and suggests a maximum of ∼15 cm of seaward fault motion, comparable to previous slow-slip events.

[1] The global positioning system (GPS) is one of the most common techniques, and the current state of the art, used to monitor volcano deformation. In addition to slow (several centimeters per year) displacement rates, GPS can be used to study eruptions and intrusions that result in much larger (tens of centimeters over hours-days) displacements. It is challenging to resolve precise positions using GPS at subdaily time intervals because of error sources such as multipath and atmospheric refraction. In this paper, the impact of errors due to multipath and atmospheric refraction at subdaily periods is examined using data from the GPS network on Kīlauea Volcano, Hawai'i. Methods for filtering position estimates to enhance precision are both simulated and tested on data collected during the June 2007 intrusion and eruption. Comparisons with tiltmeter records show that GPS instruments can precisely recover the timing of the activity.

[1] Several slow slip events beneath the south flank ofKilauea Volcano, Hawaii, have been inferred from transient displacements in daily GPS positions. To search for smaller events that may be close to the noise level in the GPS time series, we compare displacement fields on Kilauea's south flank with displacement patterns in previously identified slow slip events. Matching displacement patterns are found for several new candidate events, although displacements are much smaller than previously identified events. One of the candidates, 29 May 2000, is coincident with a microearthquake swarm, as are all of the previously identified slow slip events. The microearthquakes follow the onset of slow slip, implying that they are triggered by stress changes during slip. The new slow slip event brings the total number of events on Kilauea, between 1997 and 2007, to eight, the smallest having M w = 5.3, and the largest having M w = 6.0. While the recurrence time between the four largest events is 2.11 ± 0.01 years, the repeat time for all eight events is 0.9 ± 0.6 years. We invert for the fault geometry and distribution of slip during the slow slip events. The optimal source depths of 5 km, assuming uniform slip dislocations in an elastic half-space, are considerably shallower than the accompanying swarm earthquakes (6.5-8.5 km), which would place the earthquakes in a zone of decreased Coulomb stress. Inversions including the effects of topography and layered elastic structure in the forward models favor depths comparable to microearthquake depths, such that the earthquakes are located in a region of increased Coulomb stress. We also invert for time-dependent fault slip directly from the 30 s GPS phase observations, constraining the source to the optimal uniform slip geometry. On the basis of these inversions, the larger events last between 1.5-2.2 days. The data are unable to resolve migration of slip along the fault. The temporal pattern of accompanying microearthquakes is consistent with the fault slip history assuming a seismicity rate theory based on rate and state-friction, making the swarm earthquakes coshocks and aftershocks of the slow slip events.

[1] 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.

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/ AAA p . Values lower than this suppress small signals, and values greater than this have significantly higher noise at periods of 1-6 hours.