Modeling of volcanic processes is limited by a lack of knowledge of the time scales of storage, mixing, and final ascent of magmas into the shallowest portions of volcanic plumbing systems immediately prior to eruption. It is impossible to measure these time scales directly; however, micro-analytical techniques provide indirect estimates based on the extent of diffusion of species through melts and crystals. We use diffusion in olivine phenocrysts from the A.D. 1959 Kīlauea Iki (Hawai’i, USA) eruption to constrain the timing of mixing events in the crustal plumbing system on time scales of months to years before eruption. The time scales derived from zonation of Fe-Mg in olivines, combined with contemporaneous geophysical data, suggest that mixing occurred on three time scales: (1) as much as 2 yr prior to eruption in the deep storage system; (2) in a shallow reservoir, between incoming hot melts and resident melt for several weeks to months prior to eruption; and (3) in the conduit and summit reservoir, between the resident magma and cooled surface lava, draining back into the vent on time scales of hours to several days during pauses between episodes. Synchronous inflation of the shallow reservoir with deep earthquake swarms and mixing suggests an intermittently open transcrustal magmatic system.
Determining the time scales of pre-eruptive and syn-eruptive magmatic processes is a primary objective of modern volcanology. Time scales of magma movement through the lithosphere prior to eruptions are difficult to resolve unequivocally with geophysical techniques, yet understanding these time scales is critical for eruption forecasting. Swarms of earthquakes at a range of depths often occur prior to eruptions and are inferred to be caused by magma movement. It has not been possible to state with certainty whether these earthquakes are linked with the magma that eventually erupts. However, the chemistry and microtextures of crystals may be used as rate meters of magma mixing (e.g., Costa et al., 2009; Kahl et al., 2011; Martin et al., 2008; Ruprecht and Plank, 2013). Analysis of the chemical stratigraphy of crystals, using time scales from diffusion chronometry (e.g., Costa et al., 2003, 2008, 2009; Costa and Chakraborty, 2004; Costa and Dungan, 2005; Morgan et al., 2004; Morgan and Blake, 2005; Nakamura, 1995), provides constraints on magma ascent rates and crystal residence times that may be compared with contemporaneous geophysical data acquired prior to and during a volcanic eruption (e.g., Kahl et al., 2013).
We apply this method to the 1959 eruption of Kīlauea Iki, Kīlauea Volcano, Hawai’i, USA (Fig. 1), the type example of a powerful Hawaiian fountaining eruption. It included the highest fountaining ever observed at Kīlauea (Richter et al., 1970) and was the first eruption at Kīlauea that was rigorously geophysically monitored (Eaton and Murata, 1960; Eaton et al., 1987). The magma budget is well estimated by tilt measurements (Eaton et al., 1987), direct observation (Richter et al., 1970), and geochemical studies (Helz, 1987; Anderson and Brown, 1993; Sides et al., 2014). The eruption presents an excellent opportunity to retrospectively constrain the timing of magma ascent and mixing with respect to eruption triggering. It was preceded by significant seismicity and ground deformation linked to magma movement (Eaton et al., 1987) and involved 17 episodes of eruption (Richter et al., 1970). Geophysical and geochemical observations suggest that magma stored in a shallow reservoir was mixed, during and between episodes, with primitive magmas from depth as well as with cooler, outgassed lava that had drained back from a lava lake formed in episode 1 (Murata and Richter, 1966; Richter and Murata, 1966; Wright, 1973; Eaton et al., 1987; Helz, 1987; Loewen, 2013; Sides et al., 2014).
Our aim is to constrain the timing of magma mixing events prior to and during the A.D. 1959 eruption of Kīlauea Iki through the characterization of Mg-Fe zonation in olivine crystals. These time scales are then compared with the timings of deformation and seismicity recorded during the eruption. We seek to understand time scales of magma storage and transport to better constrain our understanding of the significance of seismic swarms that traverse the mantle and crust prior to eruptions at Kīlauea Volcano.
When a crystal, in equilibrium with a melt, experiences a change in intrinsic parameters of the host magma (e.g., pressure, temperature, or composition), the system’s chemistry will respond to re-approach equilibrium conditions. If the crystal remains at conditions such that significant diffusion of a particular element is possible, then over time the composition profile of that element at the crystal boundary will relax, creating a smoothed compositional profile whose shape reflects the diffusivity of the element and the relaxation time (Costa and Morgan, 2011). In principle, it is possible to estimate the time elapsed since the change in intrinsic parameters using the geometry of the diffusion profile (Costa and Morgan, 2011) and a diffusion coefficient. For Fe-Mg in olivine, the coefficient of diffusion is dependent upon temperature, oxygen fugacity, pressure, composition, and crystallographic orientation. By comparing the measured composition profile to a set of profiles with specific initial and boundary conditions, assuming a diffusion coefficient, a time scale can be determined between the mixing-transport event and quenching (e.g., Hartley et al., 2016). We use high-resolution backscatter images of 70 crystallographically oriented olivine crystals from 6 episodes of the 1959 eruption. Profiles along 44 olivines were calibrated for Fe-Mg composition by electron probe microanalysis. Of these, 24 compositional profiles reconstructed from backscatter images were used for diffusion chronometry modeling. Several of the 24 modeled profiles had more than one compositional zone; therefore, 29 diffusional zones were modeled. For full details of the samples, methods, and uncertainties, see the GSA Data Repository1.
The profiles were modeled using temperatures calculated by using the Mg thermometer for matrix glass (Helz and Thornber, 1987) and temperatures calculated to be in equilibrium with olivine core compositions. These temperatures provide upper and lower bounds, unique to each crystal, on the temperature (Table 1; see the Data Repository). The oxygen fugacity (fO2) of the system was assumed to be Δlog fO2 = –1.7, relative to the Ni-NiO (NNO) buffer (Roeder et al., 2003); however, using a higher oxygen fugacity of Δlog fO2 = –0.5, relative to the NNO buffer (Gerlach, 1993), reduces the time scales by a factor of ∼1.6. This uncertainty is small compared to the uncertainty due to temperature.
The composition of the olivine cores is typically more magnesium (forsterite, Fo) rich than the rim, i.e., they are normally zoned. Olivine cores typically contain ∼87 mol% Fo but range between 82.7 and 88.2 mol% Fo. Rim and intermediate compositions are more variable, from 80 to 88 mol% Fo, with a broad peak at ∼82.5 mol% Fo. There is a paucity of olivine compositions of 84–86 mol% Fo (Fig. 2).
Figure 3 shows the ages of mixing events as determined from diffusion chronometry for crystals from the six studied episodes of the eruption. The compositional zones may be divided into distinct populations: 15 crystals displayed a single normal zone formed from –600 days; 6 crystals had broad reverse zones (although 2 of these were not modeled due to the effects of nonuniform initial composition), formed from –450 days, and 4 of these had a much narrower normal zone at the crystal rim, formed in <31 days. The remaining three crystals displayed complex zonation patterns with a large range of time scales. Some crystals displayed very thin (<3 µm wide) deviations from the modeled profiles at the crystal rim, likely to represent the eruptive quench event; these were not modeled. Later episodes contain a larger proportion of crystals with more than one diffusional zone and, within these, the youngest zone (at the rim of the crystal) always corresponds to times close to or after the beginning of the eruption.
The time scales calculated from diffusion modeling using the temperatures calculated to be in equilibrium with the olivine cores are approximately one-third of those calculated using the temperatures from the glass compositions. The inferred ages are shown in Figure 4 along with contemporaneous geophysical data and a schematic diagram of the plumbing system beneath Kīlauea Volcano. The crystals record mixing events up to ∼–820 days. Mixing events became increasingly frequent toward the eruption onset; 50% of the mixing events occurred from ∼–50 days onward. Mixing took place quasi-continuously during the 35 days of eruption, generating both reverse and normal zoning.
The distribution of compositions of the olivine cores and rims, combined with the zoning patterns, suggests a complex petrogenetic history for the crystal cargo of this eruption. The wide range of olivine core compositions (Fig. 2) has been linked to crystal mush disaggregation supplying primitive olivines, with fractionation during shallow storage producing the olivines with the lower Fo cores (Clague et al., 1995).
There are several possible causes of the zonation identified here. In all crystals diffusion acts due to disequilibrium across the crystal face, thus the diffusion is caused by changes in the intrinsic properties of temperature, pressure, oxygen fugacity, and/or melt composition. Of these, temperature and melt composition are most significant for changing the equilibrium chemistry of olivine (Roeder and Emslie, 1970; Ford et al., 1983). Cooling, or a decrease in MgO content, causes normal zoning while sub-liquidus heating, or an increase in MgO content, causes reverse zoning.
Approximately 60% of the mixing events recorded by the olivines took place after the onset of shallow (<10 km) seismicity (Fig. 4). The number of mixing events increases rapidly with time up to the eruption onset, paralleling the increasing frequency of upper crustal seismicity from 50 events/day initially to more than 1000 events/day by –14 days (Eaton et al., 1987). The mixing events and seismicity in the upper crust over the few weeks prior to eruption plausibly represent the intrusion of melts from depth into a shallow summit reservoir. We are unable to distinguish between the south caldera and Halema’uma’u reservoirs as the location of mixing, but it is clear that magmas were hosted, albeit briefly, and mixed in one of the reservoirs.
Approximately 40% of the mixing events recorded by the olivines occurred before the onset of shallow seismicity (Fig. 4). During this time period, summit tilt meters recorded quasi-continuous inflation and sporadic deep (45–65 km) swarms of earthquakes (as much as 600 days before eruption) were recorded. This period was marked by an absence of shallow seismicity (Eaton et al., 1987). Although these old zones may reflect mixing during the steady inflation of the summit, there is a period of time, from –350 to –150 days, with no measured inflation. We suggest that the old mixing events, before –150 days, could be associated with the deep seismic swarms recorded at 45–65 km depth. The old zones include both normal and reverse zones, implying the mixture of different composition and/or temperature melts, each bearing their own crystal cargoes, i.e., a hot, mantle-derived, MgO-rich magma bearing high-Fo olivines mixing with cooler, MgO-poor magmas bearing low-Fo olivines. Although the deep earthquake swarms cannot be linked individually to any particular magma mixing event as recorded in the olivine, there is a suggestion that mixing occurred sporadically at depth, associated with seismicity, prior to the Kīlauea Iki eruption. The time scales calculated here are consistent with estimates of time scales between melt segregation from the mantle and eruption based on U-series disequilibria suggesting a maximum of 7 yr (Girard et al., 2012). It is critical, however, that the shallow magma reservoir beneath the summit exhibited periods of inflation from –800 to –350 days, and from –150 days until the eruption (Eaton et al., 1987). This suggests that deep intrusions were supplying magma via an open system (generating little seismicity) throughout this period. This analysis implies that the time scales for magma ascent from 45 to 65 km depth might be significantly less than 800 days and may be as short as weeks to months. These short time scales would be consistent with other studies showing rapid ascent of magmas from mantle depths to the surface (Ruprecht and Plank, 2013). Furthermore, the occurrence of deep swarms of earthquakes after the start of the eruption might indicate that magmas were drawn up from depth on the time scale of the eruption, similar to the downward-propagating seismicity observed during the Eyjafjallajokull (Iceland) eruption in 2010 (Tarasewicz et al., 2012).
The volume of magma intruded into the shallow reservoir declined as the eruption continued (to near zero by episode 8; Eaton et al., 1987), while drainback volumes remained fairly constant throughout the eruption. Several of the diffusion zones indicate mixing events after episode 8; it therefore seems likely that these zones were formed as a result of mixing between the magma in the shallow reservoir and the drainback lava, which may have driven further vesiculation and new episodes of fountaining. The return of drained-back lavas to the shallow storage reservoir is supported by melt inclusion geochemistry (Sides et al., 2014; Wallace and Anderson, 1998).
This study provides the first geochemical constraint on the time scales of magma residence and mixing through the mantle and crust at Kīlauea. The 1959 Kīlauea Iki eruption was triggered by magma mixing in the shallow crust, synchronous with shallow seismicity in the weeks leading up to the eruption. Mixing occurred up to 800 days before eruption, consistent with the timing of sporadic swarms of deep earthquakes, which probably mark the intrusion of mantle-derived melts into the deep Kīlauea plumbing system. Magmas migrated up into the shallow system rapidly and quasi-continuously during the ∼2 yr preceding the eruption. Approximately 60% of the mixing events analyzed took place <100 days prior to and during the eruption, caused by magma movement in the upper crust and by drainback processes during the eruption. The link to the deep mantle was maintained even after the start of the eruption.
We acknowledge Natural Environment Research Council studentship funds (Sides) and a U.S. Geological Survey (USGS) Jack Kleinman grant, which allowed samples for this study to be collected. Don Swanson (USGS) provided invaluable advice in the field. Rae gratefully acknowledges M. Kahl for several valuable discussions. We thank M. Poland, N. Vinet, M. Rhodes, and four anonymous reviewers for their detailed and constructive comments.