Episode 49 of the Pu'u 'O'o-Kupaianaha eruption of Kilauea volcano-breakdown of a steady-state eruptive era

The Pu'u 'O'o-Kupaianaha eruption (1983-present) is the longest lived rift eruption of either Kilauea or neighboring Mauna Loa in recorded history. The initial fissure opening in January 1983 was followed by three years of episodic fire fountaining at the Pu'u 'O'o vent on Kilauea's east rift zone 19km from the summit (episodes 4–47). These spectacular events gave way in July 1986 to five and a half years of nearcontinuous, low-level effusion from the Kupaianaha vent, 3km to the cast (episode 48). A 49th episode began in November 1991 with the opening of a new fissure between Pu'u 'O'o and Kupaianaha. this three week long outburst heralded an era of more erratic eruptive behavior characterized by the shut down of Kupaianaha in February 1992 and subsequent intermittent eruption from vents on the west flank of Pu'u 'O'o (episodes 50 and 51). The events occurring over this period are due to progressive shrinkage of the rift-zone reservoir beneath the eruption site, and had limited impact on eruption temperatures and lava composition.     

Volume of magma accumulation or withdrawal estimated from surface uplift or subsidence,with application to the 1960 collapse of Kilauea volcano

An elastic point source model proposed by Mogi for magma chamber inflation and deflation has been applied to geodetic data collected at many volcanoes. The volume of ground surface uplift or subsidence estimated from this model is closely related to the volume of magma injection into or withdrawal from the reservoir below. The analytical expressions for these volumes are reviewed for a spherical chamber and it is shown that they differ by the factor 2(1-v), where v is Poisson's ratio of the host rock. For the common estimate v=0.25, as used by Mogi and subsequent workers, the uplift volume is 3/2 the injection volume. For highly fractured rocks, v can be even less and the uplift volume can approach twice the injection volume. Unfortunately, there is no single relation between the inflation of magma reservoirs and the dilation or contraction of host rocks. The inflation of sill-like bodies, for instance, generates no overall change in host rock volume. Inflation of dike-like bodies generates contraction such that, in contrast with Mogi's result, the uplift volume is generally less than the injection volume; for v=0.25, the former is only 3/4 of the latter. Estimates of volumes of magma injection or withdrawal are there-fore greatly dependent on the magma reservoir configuration. Ground surface tilt data collected during the 1960 collapse of Kilauea crater, one of the first events interpreted with Mogi's model and one of the largest collapses measured at Kilauea, is not favored by any one of a variety of deformation models. These models, however, predict substantially different volumes of both magma withdrawal and ground surface subsidence.     

Geochemical modeling of magma mixing and magma reservoir volumes during early episodes of Kīlauea Volcano’s Pu‘u ‘Ō‘ō eruption

Geochemical modeling of magma mixing allows for evaluation of volumes of magma storage reservoirs and magma plumbing configurations. A new analytical expression is derived for a simple two-component box-mixing model describing the proportions of mixing components in erupted lavas as a function of time. Four versions of this model are applied to a mixing trend spanning episodes 3–31 of Kilauea Volcano’s Puu Oo eruption, each testing different constraints on magma reservoir input and output fluxes. Unknown parameters (e.g., magma reservoir influx rate, initial reservoir volume) are optimized for each model using a non-linear least squares technique to fit model trends to geochemical time-series data. The modeled mixing trend closely reproduces the observed compositional trend. The two models that match measured lava effusion rates have constant magma input and output fluxes and suggest a large pre-mixing magma reservoir (46±2 and 49±1 million m³), with little or no volume change over time. This volume is much larger than a previous estimate for the shallow, dike-shaped magma reservoir under the Puu Oo vent, which grew from ∼3 to ∼10–12 million m³. These volumetric differences are interpreted as indicating that mixing occurred first in a larger, deeper reservoir before the magma was injected into the overlying smaller reservoir. Electronic Supplementary Material Supplementary material is available at and is accessible for authorized users.     

Estimating Crustal Heterogeneity from Double-difference Tomography

Seismic velocity parameters in limited, but heterogeneous volumes can be inferred using a double-difference tomographic algorithm, but to obtain meaningful results accuracy must be maintained at every step of the computation. MONTEILLER et al. (2005) have devised a double-difference tomographic algorithm that takes full advantage of the accuracy of cross-spectral time-delays of large correlated event sets. This algorithm performs an accurate computation of theoretical travel-time delays in heterogeneous media and applies a suitable inversion scheme based on optimization theory. When applied to Kilauea Volcano, in Hawaii, the double-difference tomography approach shows significant and coherent changes to the velocity model in the well-resolved volumes beneath the Kilauea caldera and the upper east rift. In this paper, we first compare the results obtained using MONTEILLER et al.'s algorithm with those obtained using the classic travel-time tomographic approach. Then, we evaluated the effect of using data series of different accuracies, such as handpicked arrival-time differences (``picking differences'), on the results produced by double-difference tomographic algorithms. We show that picking differences have a non-Gaussian probability density function (pdf). Using a hyperbolic secant pdf instead of a Gaussian pdf allows improvement of the double-difference tomographic result when using picking difference data. We completed our study by investigating the use of spatially discontinuous time-delay data.     

Inflation along Kilauea's Southwest Rift Zone in 2006

We report on InSAR and GPS results showing the first crustal inflation along the southwest rift zone at Kilauea volcano in over 20 years. Two independent interferograms (May 2–August 2, 2006 and June 22–Nov 7, 2006) from the ALOS PALSAR instrument reveal domal uplift located southwest of the main caldera. The uplift is bounded on the northeast by the caldera and follows the southwest rift zone for about 12 km. It is approximately 8 km wide. We use data derived from permanent GPS stations to calibrate the InSAR displacement data and estimate uplift of 7.7 cm during the first interferogram and 8.9 cm during the second with line-of-sight volumes of 2.8 × 10⁶ m³ and 3.0 × 10⁶ m³ respectively. The earthquake record for the periods before, during, and after inflation shows that a swarm of shallow earthquakes (z < 5 km) signaled the beginning of the uplift and that elevated levels of shallow seismicity along the rift zones occurred throughout the uplift period. GPS data indicate that the inflation occurred steadily over nine months between mid-January and mid-October, 2006 making injection of a sill unlikely. We attribute the inflation to recharge of a shallow ductile area under the SWRZ.     

Petrology of gabbroic xenoliths in 1960 Kilauea basalt: crystalline remnants of prior (1955) magmatism

The 1960 Kapoho lavas of Kilauea’s east rift zone contain 1–10 cm xenoliths of olivine gabbro, olivine gabbro-norite, and gabbro norite. Textures are poikilitic (ol+sp+cpx in pl) and intergranular (cpx+pl±ol±opx). Poikilitic xenoliths, which we interpret as cumulates, have the most primitive mineral compositions, Fo_82.5, cpx Mg# 86.5, and An_80.5. Many granular xenoliths (ol and noritic gabbro) contain abundant vesicular glass that gives them intersertal, hyaloophitic, and overall ‘open’ textures to suggest that they represent ‘mush’ and ‘crust’ of a magma crystallization environment. Their phase compositions are more evolved (Fo_80–70, cpx Mg# 82–75, and An_73–63) than those of the poikilitic xenoliths. Associated glass is basaltic, but evolved (MgO 5 wt%; TiO₂ 3.7–5.8 wt%). The gabbroic xenolith mineral compositions fit existing fractional crystallization models that relate the origins of various Kilauea lavas to one another. FeO/MgO crystal–liquid partitioning is consistent with the poikilitic ol-gabbro assemblage forming as a crystallization product from Kilauea summit magma with ∼8 wt% MgO that was parental to evolved lavas on the east rift zone. For example, least squares calculations link summit magmas to early 1955 rift-zone lavas (∼5 wt% MgO) through ∼28–34% crystallization of the ol+sp+cpx+pl that comprise the poikilitic ol-gabbros. The other ol-gabbro assemblages and the olivine gabbro-norite assemblages crystallized from evolved liquids, such as represented by the early 1955 and late 1955 lavas (∼6.5 wt% MgO) of the east rift zone. The eruption of 1960 Kapoho magmas, then, scoured the rift-zone reservoir system to entrain portions of cumulate and solidification zones that had coated reservoir margins during crystallization of prior east rift-zone magmas. Received: January 7, 1993/Accepted: November 23, 1993     

The emplacement of pahoehoe lavas on Kilauea and in the Deccan Traps

There is a growing interest in deciphering the emplacement and environmental impact of flood basalt provinces such as the Deccan, India. Observations of active volcanism lead to meaningful interpretations of now-extinct volcanic systems. Here, I illustrate and discuss the morphology and emplacement of the modern and active lava flows of Kilauea volcano in Hawaii, and based on them, interpret the compound pahoehoe lavas of the Deccan Traps. The latter are vastly larger (areally extensive and voluminous) than Kilauea flows, and yet, their internal architecture is the same as that of Kilauea flows, and even the sizes of individual flow units often identical. Many or most compound flows of the Deccan Traps were emplaced in a gentle, effusive, Kilauea-like fashion. Bulk eruption rates for the Deccan province are unknown, and were probably high, but the local eruption rates of the compound flows were no larger than Kilauea’s. Large (≥ 1000 km³) individual compound pahoehoe flows in the Deccan could have been emplaced at Kilauea-like local eruption rates (1 m³/sec per metre length of fissure) in a decade or less, given fissures of sufficient length (tens of kilometres), now exposed as dyke swarms in the province.     

Emplacement of subaerial pahoehoe lava sheet flows into water: 1990 Kūpaianaha flow of Kilauea volcano at Kaimū Bay, Hawai`i

Episode 48 of the ongoing eruption of Kilauea, Hawai`i, began in July 1986 and continuously extruded lava for the next 5.5 years from a low shield, Kūpaianaha. The flows in March 1990 headed for Kalapana and inundated the entire town under 15–25 m of lava by the end of August. As the flows advanced eastward, they entered into Kaimū Bay, replacing it with a plain of lava that extends 300 m beyond the original shoreline. The focus of our study is the period from August 1 to October 31, 1990, when the lava buried almost 406,820 m² of the 5-m deep bay. When lava encountered the sea, it flowed along the shoreline as a narrow primary lobe up to 400 m long and 100 m wide, which in turn inflated to a thickness of 5–6 m. The flow direction of the primary lobes was controlled by the submerged delta below the lavas and by damming up lavas fed at low extrusion rates. Breakout flows through circumferential and axial inflation cracks on the inflating primary lobes formed smaller secondary lobes, burying the lows between the primary lobes and hiding their original outlines. Inflated flow lobes eventually ruptured at proximal and/or distal ends as well as mid-points between the two ends, feeding new primary lobes which were emplaced along and on the shore side of the previously inflated lobes. The flow lobes mapped with the aid of aerial photographs were correlated with daily observations of the growing flow field, and 30 primary flow lobes were dated. Excluding the two repose periods that intervened while the bay was filled, enlargement of the flow field took place at a rate of 2,440–22,640 square meters per day in the bay. Lobe thickness was estimated to be up to 11 m on the basis of cross sections of selected lobes measured using optical measurement tools, measuring tape and hand level. The total flow-lobe volume added in the bay during August 1–October 31 was approximately 3.95 million m³, giving an average supply rate of 0.86 m³/s.     

Pit crater formation on Kilauea volcano, Hawaii.

Most of the known pit craters in Hawaii occur along the East and Southwest Rift Zones of Kilauea volcano. The pit craters typically are either astride a single rift zone fracture or between a pair of rift zone fractures. These fractures are prominent in the pit crater walls. The pit craters are elliptical in plan view, with their major diameters ranging from 8 to 1140 m. They range in depth from 6 m to 186 m. They typically develop with initially steep, locally overhanging walls, but as the walls collapse, the craters fill with talus and become shaped like inverted elliptical cones. None of the craters apparently formed as eruptive vents, although some have been subsequently filled by lava. Devil's Throat is the best-exposed pit crater along the East Rift Zone. It is sited at a `waist' between two east-striking zones of ground cracks; the spacing between the crack zones decreases towards Devil's Throat. East-striking fractures are also prominent in the pit crater walls. Pit craters along the Southwest Rift Zone typically are elongate in plan view along the direction of the rift, have large caves at their bases along the long axes of the craters, and are smaller than those of the East Rift Zone. Some closely spaced pits there have coalesced to form a trough. Based on our observations and mechanical considerations, we infer that pit craters form by stoping over an underlying large-aperture rift zone fracture, and not by piston-like collapse over broad magma bodies or voids. Flow of magma along the underlying fracture may remove stoped blocks and prevent the fracture from being choked with debris. This mechanism is consistent with pit crater location, ground crack patterns, the preferred orientation of fractures in pit crater walls, and pit crater geometry (both in map view and cross-section). The mechanism also fits with observations of stoping into a gaping rift fracture that conducted lava from Kilauea caldera during the 1920s. Additionally, the ratio of pit crater width to depth of 0.5 to 2 is consistent with pit craters forming over a nearly vertical opening mode fracture.     

The relationship between earthquake swarms and magma transport: Kilauea volcano,Hawaii

A detailed investigation of earthquake locations and focal mechanisms for swarms associated with intrusive events at Kilauea volcano, Hawaii, further illuminates the relationships among stress state, faulting, and magma transport. We determine the earthquake locations and mechanisms using a three-dimensional crustal model to improve their accuracy and consistency. Swarms in Kilauea's upper east and southwest rift zones, from the years 1980 through 1982, provide clear evidence for the propagation and/or dilation of dikes. Focal mechanisms are predominantly strike-slip, and the faulting and inferred dike orientations can be interpreted quite consistently in terms of the model ofHill (1977). Stresses induced by the summit magma reservoir system strongly control faulting and magma transport in the rift zones close to the summit.