Skip Navigation



Journal of Petrology Advance Access published online on July 19, 2007

Journal of Petrology, doi:10.1093/petrology/egm034
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
49/4/757    most recent
egm034v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Garcia, M. O.
Right arrow Articles by Eiler, J. M.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2007. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org

Oxygen Isotope Evidence for Chemical Interaction of Ki lauea Historical Magmas with Basement Rocks

Michael O. Garcia1,*, Emi Ito2 and John M. Eiler3

1Department of Geology & Geophysics, University of Hawai’i, Honolulu, HI 96822, USA
2Department of Geology & Geophysics, University of Minnesota, Minneapolis, MN 55455, USA
3Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125, USA

Received January 30, 2007; Revised typescript accepted June 13, 2007


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Kilauea historical summit lavas have a wide range in matrix {delta}18OVSMOW values (4·9–5·6{per thousand}) with lower values in rocks erupted following a major summit collapse or eruptive hiatus. In contrast, {delta}18O values for olivines in most of these lavas are nearly constant (5·1 ± 0·1{per thousand}). The disequilibrium between matrix and olivine {delta}18O values in many samples indicates that the lower matrix values were acquired by the magma after olivine growth, probably just before or during eruption. Both Mauna Loa and Kilauea basement rocks are the likely sources of the contamination, based on O, Pb and Sr isotope data. However, the extent of crustal contamination of Kilauea historical magmas is probably minor (< 12%, depending on the assumed contaminant) and it is superimposed on a longer-term, cyclic geochemical variation that reflects source heterogeneity. Kilauea's heterogeneous source, which is well represented by the historical summit lavas, probably has magma {delta}18O values within the normal mid-ocean ridge basalt mantle range (5·4–5·8{per thousand}) based on the new olivine {delta}18O values.

KEY WORDS: Hawaii; Kilauea; basalt; oxygen isotopes; crustal contamination


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Kilauea volcano has undergone dramatic changes over the last ~320 years varying from continuous effusion to explosive eruptions with major summit collapses (>100 m) (Macdonald et al., 1983Go; Holcomb 1987Go; Garcia et al., 2003Go). The 1924 summit explosions and collapse were accompanied by a reversal in the temporal trend in lava composition (including incompatible trace element abundances and Pb, Sr and Nd isotope ratios) that was interpreted to reflect changes in melting processes and compositional heterogeneity in the source (Pietruszka & Garcia, 1999aGo; Garcia et al., 2003Go). The additional effects of crustal contamination were superimposed on the compositions of lavas erupted just after the collapse (Pietruszka & Garcia, 1999aGo). Kilauea's historical period has been punctuated by numerous collapses (1823, 1832, 1840, 1868, 1924; Table 1). Here we evaluate the role of crustal contamination during the entire historical period using O isotopes, which have been shown to be a sensitive indicator of crustal contamination (e.g. Taylor, 1974Go).


View this table:
[in this window]
[in a new window]

 
Table 1: Historical Kilauea summit geological events

 
Crustal contamination has been documented at various oceanic islands (e.g. Pitcairn; Eiler et al. 1995Go; Canary, Thirlwall et al., 1997Go; Hawai’i, Gaffney et al., 2005Go) including the current eruption of Kilauea, where magmas were stored for many years in the rift zone prior to eruption (Garcia et al., 1998Go). It should also be noted that {delta}18O variations could also reflect source heterogeneity (e.g. Eiler 2001Go). To unravel the potentially offsetting influences of source heterogeneity vs crustal contamination, we analyzed O isotopes in coexisting olivine and matrix material from well-characterized historical lavas from Kilauea volcano (20 eruptions during ~400 years; Pietruszka & Garcia, 1999aGo; Garcia et al., 2003Go). These results provide a detailed time series of O isotopes, which are compared with those for another active hotspot volcano, Grimsvötn in Iceland, where 800 years of volcanic activity (including the Laki flood basalt) have been investigated (Bindeman et al., 2006Go). After evaluating the effects of crustal contamination on Kilauea historical lavas, we examine the issue of {delta}18O values for the Hawaiian plume source feeding this volcano.


    KiLAUEA ERUPTIVE HISTORY SUMMARY
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Kilauea's historical activity (~1790 to present) includes prolonged periods of quiescent lava effusion, interrupted by major collapses (≥100 m) of the caldera floor, rare violent explosions, and brief to prolonged rift zone eruptions (Table 1). The most dramatic event occurred in ~1790, when violent explosions killed ~80 Hawaiians (McPhie et al., 1990Go). This event was accompanied by a major lower east rift zone eruption (Macdonald et al., 1983Go). It was preceded by several explosive eruptions dating back to AD ~1670 (Swanson et al., 2006Go) and followed, sometime between 1820 and 1823, by high lava fountaining within the caldera (Sharp et al., 1987Go). In 1823, a ‘great’ crack extending 22 km down the SW rift produced a large lava flow causing the summit lava lake to drain and a summit collapse of >100 m (Brigham, 1909Go; Macdonald et al., 1983Go). Eruptive activity resumed until 1832, when cracks formed on the NE rim of the caldera producing a small lava flow, which was followed by the summit lava lake draining and a >100 m collapse of the caldera floor (Brigham, 1909Go). Lava-lake activity restarted a few months later and persisted until 1840, when the caldera floor sank ~100 m during a major east rift zone eruption (Dana, 1891Go). Formula Crater, in the SW corner of the caldera (Fig. 1), was the focus of subsequent lava lake activity until a M ~7·9 earthquake in 1868 caused the caldera floor to sink ~90 m to ~180 m (Brigham, 1909Go; Wyss & Koyanagi, 1992Go). Vigorous activity resumed immediately in Formula and continued until 1894, when the lava lake drained. For the next 13 years, eruptive activity was episodic with a 5·5 year hiatus (Brigham, 1909Go). Lava-lake activity returned to Formula in 1907, with major overflows and fissure eruptions from 1918 to 1921, covering large portions of the caldera floor.


Figure 1
View larger version (78K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1. Topographic map of the island of Hawai’i showing the five shield volcanoes that form the island [from youngest to oldest; Kilauea, Mauna Loa, Hualalai (Hu), Mauna Kea (MK), Kohala (Ko)] and the young submarine volcano Loihi (Lo). Thick gray lines separate the shield volcanoes. Bold long-dashed line marks the location of the schematic cross-section (Fig. 6). Also shown as bold short-dashed lines are Kilauea's two rift zones, including the east rift zone (ERZ) where the current Figure 1 eruption is occurring (P) and the 2·4 km deep KS-3 well is located (spot). Contour interval is 500 m.

 
In 1924, the Formula lava lake drained, its crater floor collapsed ~300 m, and numerous violent phreatic explosions ensued (Jaggar, 1924Go). This event marked a reversal in Kilauea's lava compositional trend that has persisted for >80 years (Pietruszka & Garcia, 1999aGo; Garcia et al., 2003Go) and a major change in its eruptive style. Summit eruptions became short and infrequent (seven brief eruptions during the next 10 years; Jaggar, 1924Go; Finch, 1940Go). From 1934 to 1952, the volcano remained quiet. The next 30 years (1952–1982) witnessed mostly short-lived eruptions from fissures in and near Formula crater (Fig. 1) and the more frequent east and SW rift zone eruptions, following a 115 year hiatus in rift zone eruptive activity (Macdonald et al., 1983Go). There have been no summit eruptions since the start of the current east rift zone, Formula eruption in 1983 (e.g. Garcia et al., 2000Go).


    SAMPLES
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Well-dated and chemically characterized summit flow and tephra samples were utilized for this study [see Pietruszka & Garcia (1999aGo) and Garcia et al. (2003Go) for locations, petrography and compositional information (all for XRF major and trace elements, and many for Pb, Sr and Nd isotopes)]. This sample suite spans ~310 years of Kilauea's eruptive history (~1670 to September 1982) and includes two young prehistoric samples [dates of 1790 were given by Pietruszka & Garcia (1999aGo) based on the work of McPhie et al. (1990Go); however, subsequent 14C dating shows these samples are ~120 years older (Swanson et al., 2006Go) than previously thought]. Unfortunately, only one sample (1832) is available for the period 1832–1865, during which two major collapse events occurred (1832 and 1840). Kilauea erupted just prior to the 1832 collapse. Subsequent major collapses (1868 and 1924) and eruptive gaps (1897–1902 and 1934–1954) are well bracketed by the sample suite (Table 2). A previous study determined {delta}18O values on whole-rocks for five historical summit samples, including samples from three of the same eruptions examined in our study (1885, 1954 and 1971), obtaining {delta}18O values (all values in this paper are reported relative to Vienna Standard Mean Ocean Water, VSMOW), of 5·5–5·6{per thousand} (Kyser et al., 1982Go). Another study of O isotopes of olivine in 10 SOH-4 drill-core lavas from Kilauea's east rift zone yielded {delta}18O values of 4·9–5·2{per thousand} , averaging 5·1 ± 0·1{per thousand} (Eiler et al., 1996Go), consistent with the Kyser et al. (1982Go) summit results. In contrast, olivines from the current Formula eruption are more variable, ranging in {delta}18O values from 4·4 to 4·9{per thousand} for 16 samples, with low values during the first 2 years of episodic eruptions (4·8 ± 0·1{per thousand}), followed by even lower values for the continuous effusion phase (4·64 ± 0·2{per thousand}; Garcia et al., 1998Go). The matrix for these and other Formula lavas (24 total) also ranges widely in {delta}18O values (4·56–5·25{per thousand}), averaging 4·8 ± 0·2{per thousand} for the early lavas and 5·1 ± 0·1{per thousand} for later lavas (Garcia et al., 1998Go).


View this table:
[in this window]
[in a new window]

 
Table 2: XRF major and trace element analyses of Kilauea rocks

 
Most of the summit samples are simple in mineralogy, consisting of rare to common olivine crystals (0·1–4·7 vol.%) in a glass matrix. Insufficient olivine was present in samples erupted between 1929 and 1934 to allow for O isotope analyses of olivine in these samples. Glass is present in most samples; in those samples where it is absent, matrix material was used for {delta}18O determinations. Olivine elemental compositions for many of the summit samples were reported by Garcia et al. (2003Go). Most are in equilibrium with the whole-rock composition, except a few crystals from prehistoric samples (1790-1) and the early 19th century samples (1820 and 1832). The 1790 and 1820 samples contain a few lower forsterite crystals (Fo 84–85) that are weakly reversely zoned and the 1832 sample has rare higher forsterite olivines (88%; Garcia et al., 2003Go). Glass MgO contents of the summit samples range widely (from 9·3 wt% in the 1820 samples to 6·2 wt% in the September 1982 sample, 1982S-12; Table 3). Whole-rock MgO contents (labeled ‘R’ in Table 3 for samples without glass) range from 12·3 to 6·6 wt%. These values are considered reasonable indicators of matrix MgO content in the olivine-poor samples (< 2·5 vol.%). Sample 1790-1, with 4·7 vol.% olivine, is the exception.


View this table:
[in this window]
[in a new window]

 
Table 3: Oxygen isotope data for Kilauea summit volcanics, accidental blocks and cuttings from east rift zone well

 
To evaluate possible crustal contaminants for Kilauea magmas, new samples were collected from the accidental debris from the 1924 explosions (denoted by the prefix Z followed by a number and either ‘b’ for basalt or ‘g’ for gabbro) and from well cuttings collected every ~3 m from the Puna Geothermal Venture well KS-3 in the lower east rift zone of Kilauea (labeled PGV followed by the depth of sample collection in meters). This well, which reached temperatures of 320°C (D. Thomas, personal communication, 2006), provides deep samples (up to 2237 m) from the interior of Kilauea volcano. The cuttings are 2–4 mm across and most have a gray matrix with sparse to common olivine. Within each sample interval studied, the chips are petrographically similar.

The accidental blocks were derived from the shallow section (< 300 m) overlying the summit magma reservoir (Jagger, 1924Go; Macdonald, 1944Go) and are thought to be possible contaminants for the summit reservoir, especially the 1924 collapse (Pietruszka & Garcia, 1999aGo). The blocks range in diameter from 25 to 150 cm, contain abundant olivine (~30 vol.%) and have a grayish green coating. Sample Z19b is distinctive for its black olivine, which has been interpreted to result from high-temperature metamorphism (>700°C; Haggerty & Baker, 1967Go). Olivines in sample Z19b range from forsterite 94 to 98, the highest ever reported for a Hawaiian rock. These high forsterite contents result from exsolution of magnetite within olivine. Jaggar (1924Go) measured accidental block temperatures to 700°C, a minimum pre-eruption value. Mineral compositions (olivine, magnetite and orthopyroxene) within this block yielded temperatures of 970 ± 170°C using the Quilf95 program (D. Lindsley, personal communication, 1999). Sample Z6g has mildly oxidized olivine (metallic gray color), whereas olivines in sample Z2b are unoxidized. Thus, these blocks show a broad range of metamorphic histories from mild (Z2b) to intense baking (Z19b).


    METHODS AND RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Whole-rock XRF
The new accidental block and well cutting samples were analyzed by X-ray fluorescence (XRF) at the University of Massachusetts for whole-rock major and trace (Nb, Rb, Zr, Sr, Y, V, Zn, Ni, Cr, Ba, La, and Ce) elements using methods described by Rhodes & Vollinger (2004Go). The analytical uncertainties for these analyses have been given by Rhodes & Vollinger (2004Go). Samples were washed in deionized water until the solution was clear. Clean chips were hand picked before crushing in a tungsten carbide mill for less than 2 min.

The KS-3 geothermal well cuttings have typical Hawaiian tholeiite major element compositions with MgO concentrations of 6·3–9·0 wt%. Their compositions resemble Kilauea lavas in having relatively high TiO2, moderate SiO2, and Zr/Nb of 9·6–11 (Table 2). These samples have experienced at least moderate levels of alteration based on their loss on ignition (LOI) of 1· 4–2·6 wt%, their somewhat low K2O/P2O5 compared with magmatic values (< 1·5 vs >1·5; Wright, 1971Go), and their low Rb contents (< 5 ppm, except the deepest sample). There is no correlation of these alteration indicators with the depth at which the cuttings were collected (Table 2). Except for these indicators, the major and trace element compositions of the five chip samples are remarkably similar to those of typical Kilauea lavas (e.g. Nb/Zr of 10–11; Rhodes et al., 1989Go; Garcia et al., 2000Go, 2003Go).

The accidental blocks are high-MgO tholeiites (18·9–21· 6 wt%) with correspondingly high Cr and Ni concentrations (Table 2). The blocks have a wide range in K2O, TiO2 and incompatible trace elements (Table 2), magmatic K2O/P2O5 (1· 6–2·1) and low LOI values (<0·25 wt%). The blocks have Kilauea-like Zr/Nb ratios (12–13; Rhodes et al., 1989Go). They are petrographically and geochemically similar to the suite C picritic basalts from Kilauea's caldera cliffs (see Casadevall & Dzurisin, 1987Go).

Matrix oxygen isotopes
Oxygen isotope analyses of matrix material were made by conventional methods (Clayton & Mayeda, 1963Go) at the University of Minnesota using ClF3 (Borthwick & Harmon, 1982Go). CO2 gas from each extraction carried out during 1993 and 1994 was analyzed on a Finnigan MAT delta E mass spectrometer. CO2 gas extracted during 1997 and 1998 was analyzed with a Finnigan MAT 252 mass spectrometer. Each sample was split into 2–5 aliquots depending on the total sample weight, and at least two aliquots were analyzed to check for analytical reproducibility (see Table 3 for the number of replicates for each). For all but one sample (accident block Z19b, with the black olivines), analytical reproducibility is 0·1{per thousand} or better. For 1993 and 1994 extractions, one laboratory standard (Nain plagioclase) was analyzed for every four matrix samples. The overall reproducibility for Nain plagioclase for the duration of this work was ±0·20{per thousand} (1{sigma}) for 125 analyses. A mid-ocean ridge basalt (MORB) glass from East Pacific Rise (EPRD 10RB), which had previously been analyzed at the Geophysical Laboratory (Ito et al., 1987Go), was analyzed as a laboratory standard during 1997 and 1998. The overall reproducibility for EPRD 10RB for the two years was ±0·17{per thousand} (1 {sigma}) for 56 analyses. The results reported here are relative to {delta}18OVSMOW , normalized against SLAP (Standard Light Antarctic Precipitation, which is defined as –55{per thousand} relative to VSMOW; Coplen, 1988Go).

The matrix {delta}18O values for the summit suite of samples vary from 4·90 to 5·62{per thousand} (Table 2) with no systematic variation with the date of eruption (Fig. 2). PVG K3-3 cuttings show a wide variation {delta}18O from 1·75 to 5·31{per thousand}, generally decreasing with depth (Table 3). The three accidental blocks are remarkably similar in their O isotope values (5·50 ± 0·04{per thousand}) given their wide range in degree of metamorphism.


Figure 2
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2. Temporal variation in matrix oxygen isotope values for Kilauea summit ({blacktriangleup}) and Figure 2 ({circ}) lavas. Gray vertical bars show times of major collapses (dates are listed at top of figure); stippled zones show major gaps in eruptive activity (see Table 1 for details). There are marked decreases in matrix oxygen isotope values following summit collapses in 1868 and 1924, and following eruptive hiatuses in ~1897–1902 and 1934–1952. Figure 2 data are shown only for the continuous phase of activity from 1986 to 1997 (data from Garcia et al., 1998Go). MORB glass values: {square}, Ito et al. (1987Go); {blacksquare}, Eiler (2001Go). MORB glass error bar represents combined range from these two studies. Other bars are the maximum standard variation based on 2–3 analyses of each sample (see Table 3 for specific samples).

 
Olivine oxygen isotopes
Olivines were separated from historical Kilauea lavas and tephra by hand picking from coarsely crushed, but otherwise untreated, samples. Oxygen isotope compositions of c. 2 mg aliquots of these separates were determined by laser fluorination, using a 50W CO2 laser and BrF5 as reagent (Valley et al., 1995Go). Product O2 was converted to CO2 by reaction with hot graphite; CO2 was then analyzed for its isotopic composition by dual-inlet gas source mass spectrometry on a Thermo Finnegan 253 system at California Institute of Technology. Data are reported in units of per mil versus the VSMOW standard. Analyses were standardized by comparison with measurements of Gore Mountain garnet standard (Valley et al., 1995Go) and a San Carlos olivine (Eiler et al., 1995Go, 1996Go; Bindeman et al., 2006Go). Both standards were analyzed between two and five times on each day of analysis, and data for unknowns analyzed on that day were corrected by the average difference between measured and accepted values for these standards. The external precision of repeat measurements of separate splits of unknown samples averaged ±0·05{per thousand} (1{sigma}). This is comparable with the typical external precision for replicate measurements of silicate standards for this laboratory and technique (e.g. Eiler et al., 1995Go, 1996Go; Bindeman et al., 2006Go), suggesting that the olivine separates analyzed in this study are homogeneous in {delta}18O at the scale of c. 2 mg aliquots. Garcia et al. (1998Go) evaluated the interlaboratory differences between Caltech and the University of Minnesota by comparing the results for a matrix glass from the Formula eruption. The sample was run in duplicate in both laboratories. The results were 4·68{per thousand} by conventional methods vs 4·79{per thousand} by laser, which are within the analytical error of the two methods.

Ten of 13 olivine separates yield average {delta}18O values that group tightly around a value of 5·11 ± 0·07{per thousand} (Fig. 3). Members of this population are indistinguishable from one another, within analytical precision, and equal in {delta}18O to olivine in equilibrium with typical MORB (Eiler et al., 2000Go; Cooper et al., 2004Go). This average value is also similar to the {delta}18O values of olivines from Mauna Loa, Lo ’ihi, shield-building Haleakala (i.e. Honomanu series), submarine Mauna Kea lavas, and olivines recovered from the Kilauea SOH4 drill core (all of which have suite-average {delta}18O values from 5·05 to 5·22; Eiler et al., 1996Go; Wang et al., 2003Go), and within the range typical of olivines from most ocean-island basalts and mantle peridotites ({delta}18O values of 5·0–5·4{per thousand}; Eiler, 2001Go).


Figure 3
View larger version (15K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3. Temporal evolution in olivine oxygen isotope values for Kilauea summit ({blacktriangleup}) and Figure 3 ({circ}) lavas. Most of the summit olivines plot within the MORB olivine field (if the errors bars are included), with no temporal variation. The exceptions are the 1961 and 1971 samples and Figure 3 lavas (only data for the continuous phase of activity are shown; Garcia et al., 1998Go). MORB olivine range is from Eiler (2001Go). Error bars are ± maximum standard variation based on replicate analyses (see Table 3).

 
Three historical Kilauea olivine samples are distinguishable from this tightly grouped set of 10 samples. Olivine from sample 1894-2 has a {delta}18O of 4·91{per thousand}, just below the range typical of mantle-derived olivines and similar to previously observed averages for the Kula series of Haleakala. Samples 1971S and 1982-A-20 have the lowest values for summit lavas, 4·70 and 4·60, respectively, overlapping the range of values for olivines from the continuous phase of Kilauea's Formula eruption (Fig. 3) and within the range typical of subaerial Mauna Kea lavas (Eiler et al., 1996Go; Wang et al., 2003Go). Relatively low {delta}18O values in Hawaiian olivines have been previously interpreted to result from lithospheric contamination (Eiler et al., 1996Go), a low-{delta}18O plume component (Lassiter & Hauri, 1998Go) or contamination within the volcanic edifice (Garcia et al., 1998Go; Wang et al., 2003Go).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
The wide range in matrix {delta}18O values (4·9–5·6{per thousand}) during the ~310 year period for Kilauea summit lavas is in striking contrast to the nearly constant matrix {delta}18O values (3·1 ± 0·1{per thousand}) for ~800 years at Iceland's Grimsvötn volcano, including the 15 km3 Laki eruption (Bindeman et al., 2006Go). The variation observed at Kilauea is a magmatic signature rather than a secondary feature because all of the summit samples are fresh, rapidly quenched volcanic rocks (Garcia et al., 2003Go). Could the Kilauea variation be related to mantle source variations, as suggested for oceanic island basalts (e.g. Harmon & Hoefs, 1995Go)? The rapid variations in Kilauea summit {delta}18O values are difficult to reconcile with components in the Hawaiian mantle, especially given the lack of correlation of {delta}18O values with incompatible trace element ratios (Nb/Y; Fig. 4) that are considered a good source indicator (Rhodes & Hart, 1995Go). Furthermore, {delta}18O values of olivines from summit lavas are relatively constant from prehistoric times to 1982 (5·1 ± 0·1{per thousand} ; Fig. 3), whereas Sr, Pb and Nd isotope ratios varied considerably during this period (e.g. 206Pb/204Pb 18·39–18·66; Pietruszka & Garcia, 1999aGo). These results are also distinct from those for subaerial HSDP2 core olivine, which showed lower average values (4·79 ± 0·13{per thousand} ; Wang et al., 2003Go). Variable amounts of peridotite partial melting or crystallization of olivine and minor spinel are likely to have opposite but minimal impacts on the matrix {delta}18O value (+0·1{per thousand} for 20% crystallization vs –0·1{per thousand} for up to 30% partial melting; Eiler 2001Go; Wang et al., 2003Go). Also, there is no correlation in the extent of fractionation (MgO content; Fig. 4) with matrix {delta}18O value (e.g. at 7· 0 MgO, {delta}18O values span the entire range of values; Table 3). Below we examine the temporal variation and discuss the possible causes for the {delta}18O variation in historical Kilauea lavas.


Figure 4
View larger version (14K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4. Matrix oxygen isotope value vs MgO ({blacktriangleup}, glass; {triangleup}, rock) and whole-rock Nb/Y in historical Kilauea summit lavas. The lack of correlations of these geochemical parameters indicates that fractionation (which is dominated by olivine in these rocks; Garcia et al., 2003Go) and source (Nb/Y is strongly influenced by source in Hawaiian basalts; e.g. Pietruszka & Garcia, 1999aGo) have not played an important role in controlling the oxygen isotope values in Kilauea summit lava. Error bars are smaller than the symbol size. Data from Table 3 and Garcia et al. (2003Go).

 
Temporal variation in oxygen isotopes
Two features emerge when the data are viewed in the context of Kilauea summit activities (Table 1): (1) matrix {delta}18O values are lower following major summit collapses (100+ m; 1868 and 1924) or significant gaps in eruptive activity (>5 years; 1898 and 1934; Fig. 2); (2) they increase during periods of nearly continuous eruptive activity between these events (1877–1894 and 1912–1921), except for the period of short, infrequent eruptions from 1952 to 1982 (when very little lava was erupted; Dvorak & Dzurisin, 1993Go) and during the continuing Formula eruption (Fig. 2). The drop in matrix {delta}18O values following summit collapses and during gaps in eruptive activity or reduced levels of effusion may represent times when the summit reservoir magma is contaminated. The second feature, increasing matrix {delta}18O values, may indicate that contaminated magma was progressively diluted with new, uncontaminated magma. New magma was more or less continuously added to the summit reservoir during historical times (e.g. Finch, 1940Go; Macdonald et al., 1983Go; Dvorak & Dzurisin, 1993Go). For the Formula eruption, there was a 0·2–0·3{per thousand} increase in matrix {delta}18O values when the eruptive style switched from episodic to continuous effusion, although there was no change in other geochemical parameters (Garcia et al., 1998Go). These fluctuations in Formula matrix {delta}18O values were interpreted to represent variable amounts of crustal contamination within Kilauea's east rift zone, with less contamination during periods of higher effusion (Garcia et al., 1998Go).

The olivine {delta}18O values for summit historical lavas record a different history. The values are all relatively high, comparable with those for MORB olivine, and nearly constant except for the 1961 and 1971 lavas (Fig. 3). The Formula lavas also have low olivine values. These results are in contrast to the wide variation in olivine isotope values for the Laki eruption (2·3–5·2{per thousand}), which were attributed to the incomplete equilibration of olivine from influxes of mantle-derived magma into a large chamber with {delta}18O-depleted magma (Bindeman et al., 2006Go). The Kilauea summit olivines apparently grew in magmas with {delta}18O values of 5·45–5·75{per thousand}, assuming equilibrium growth and an olivine–melt isotopic fractionation of 0·55{per thousand} [average of 0·7{per thousand} (Ito & Stern, 1986Go) and 0·4{per thousand} (Wang et al., 2003Go)]. The summit olivine data do not provide independent evidence for contamination except for the 1961 and 1971 summit and Formula lavas (Fig. 3).

If the matrix and olivine {delta}18O values are compared (Fig. 4), a record of disequilibrium in Kilauea summit magmas is revealed. The extent of disequilibrium was greatest following the late 17th century major explosions and progressively decreased through the end of the 19th century when equilibrium values were achieved (Fig. 4). The 1924 major summit collapse and explosions, and the 1934–1954 eruptive hiatus marked another major matrix–olivine disequilibrium event, with the magnitude of the disequilibrium decreasing from 1954 to 1971 (Fig. 5). However, the April 1982 lava and many of the Formula lavas (those erupted during periods of lower eruption rate) show disequilibrium values. No event was recorded that might explain contamination in the 1982 lavas. The presence of normal {delta}18O values in olivine (5·0–5·2{per thousand}) from the September 1982 eruption (Fig. 3) suggests that these magmas were contaminated after the olivines grew, probably during eruption. All of the post-1950 lavas have low matrix {delta}18O isotope values indicating that there was a major disruption in the magmatic plumbing system, which may have allowed for more magma contamination. As discussed above, the Formula lavas with disequilibrium matrix–olivine values were erupted during periods of lower eruption rates. In summary, the two large summit events in the late 17th and early 20th centuries had a profound impact on the {delta}18O values of Kilauea magmas. If contamination caused the drop in O isotope values, where did it occur?


Figure 5
View larger version (36K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 5. Temporal evolution in the difference between matrix and olivine oxygen isotope values for Kilauea summit ({blacktriangleup}) and Figure 5 ({circ}) lavas. Only the 1924 collapse and 1934–1952 eruption pause are shown because they mark a major reversal in matrix–ol ratios. Figure 5 data are shown only for the continuous phase of activity from 1986 to 1997 (data from Garcia et al., 1998Go). Estimates for the equilibrium fractionation between matrix and olivine range from 0·7 to 0·4{per thousand} (Ito & Stern, 1986Go; Wang et al., 2003Go). The equilibrium field shown spans these two values and covers the higher values observed for summit lavas. The values below this field are thought to reflect contamination. Error bars are the combined olivine and matrix maximum variation based on replicate analyses (see Table 3).

 

Figure 6
View larger version (22K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 6. O isotope matrix values vs whole-rock 206Pb/204Pb and 87Sr/86Sr ratios for Kilauea summit ({blacktriangleup}), Figure 6 ({circ}), and Mauna Loa submarine rift ({diamond}) lavas. The isotopes show a broad correlation, which may result from contamination of Kilauea magma by high-temperature altered Mauna Loa lavas. O isotope data for Mauna Loa from Garcia et al. (1989Go); oceanic crust upper section from King et al. (1993Go). Pb isotope data for Mauna Loa from Kurz et al. (1995Go) and Rhodes & Hart (1995Go); Kilauea data from Pietruszka & Garcia (1999aGo); Figure 6 data from Garcia et al. (1998Go, 2000Go); oceanic crust data from King et al. (1993Go) and Fekiacova et al. (2007Go). Sr isotope data for Mauna Loa from Gurriet (1988Go). Kilauea field based on samples with mantle-like O isotope values of 5·4–5·6{per thousand}, which are assumed to be uncontaminated.

 
Implications of O, Sr and Pb isotope data for contamination source and location
There is a general correlation of {delta}18O with 87Sr/86Sr ratios, to a lesser extent with 206Pb/204Pb (Fig. 6) but none with {varepsilon}Nd. The Kilauea summit samples with the lowest {delta}18O values tend to have higher Sr and somewhat lower Pb isotope ratios, although there is a broad range at a given {delta}18O value (Fig. 6). These results help to identify potential sources of the contamination that created the lower {delta}18O values in historical Kilauea summit lavas. Shallow-level, low-temperature processes should produce higher {delta}18O values (Gregory & Taylor, 1981Go). For example, values up to 8·7{per thousand} have been reported for rocks from the shallow levels of Kilauea geothermal wells (Smith & Thomas, 1990Go). Deeper in these wells where temperatures are higher (~300°C), the {delta}18O values are lower (3·2{per thousand}; Smith & Thomas, 1990Go). This trend is confirmed with our new data for well KS-3, which generally decrease with depth from 5·32{per thousand} at intermediate depths (1353 m) to 1· 86{per thousand} at 2237 m (Table 3), the lowest value reported for Hawaiian lavas. The Kilauea high-temperature metamorphosed accidental blocks from the 1924 eruption, which were derived from depths <300 m (Jaggar, 1924Go), lack low {delta}18O values (Table 3). Therefore, they are unsuitable contaminants for generating the lower {delta}18O values observed in some summit lavas (Fig. 6).

The Pacific oceanic crust section, which was invoked by Gaffney et al. (2005Go) to explain low {delta}18O values in West Maui lavas, has 206Pb/204Pb ratios too high (18·5–19·1; King et al, 1993Go; Fekiacova et al., 2007Go) to explain the range or correlation in Kilauea matrix O and Pb isotope values (Fig. 6). The only other potential contaminant for Kilauea magmas is Mauna Loa lavas (Fig. 7), which have lower 206Pb/204Pb and higher 87Sr/86Sr ratios than Kilauea lavas (Fig. 6). Glasses from fresh Mauna Loa submarine lavas have a wide range in {delta}18O values (4·7–5·5{per thousand} ; Garcia et al., 1989Go). The lavas with lower values would be suitable contaminants for Kilauea magmas to generate the positive O–Sr and O–Pb isotope correlation (Fig. 6). An alternative contaminant would be moderately hydrothermally altered Mauna Loa lavas.


Figure 7
View larger version (34K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 7. Schematic cross-section of the southern part of the island of Hawai’i (see Fig. 1 for location) showing the basement stratigraphy and magmatic plumbing system for Kilauea volcano based on regional geology and seismic cross-sections (Macdonald et al., 1983Go; Hill & Zucca, 1987Go; Klein et al., 1987Go; Garcia et al., 2000Go). Stippled zone around Kilauea conduit is a region of dunite and gabbroic cumulates. Age of oceanic crust adjusted from ODP Site 843, located 300 km west of the island of Hawai’i (Waggoner, 1993Go).

 
The relatively high {delta}18O values for the summit olivines (Fig. 3), and the disequilibrium between olivine and matrix {delta}18O values in some samples (Fig. 5), indicate that the contamination occurred late, probably just prior to or during eruption. The most likely location for this contamination would be in Kilauea's shallow (3–6 km deep; Klein et al., 1987Go) magma chamber, which is probably located within the flanks of Mauna Loa, or in the conduit during eruption (Fig. 7). The frequent, major collapses of the summit crater floor of 100+ m and the gaps in eruptive activity provide a suitable mechanism for disruption of the roof of the summit magma chamber and the assimilation of Mauna Loa or Kilauea hydrothermally altered rocks into Kilauea magma. The 1924 collapse resulted in significant changes in lava trace element chemistry, which were explained by 3·5% contamination from partially melted, seawater altered Kilauea basalt (Pietruszka & Garcia, 1999aGo). However, the trends of the lava chemistry changes following this event are consistent with a Mauna Loa source (decreasing Pb and Nd and increasing Sr isotopes). The small drop in the {delta}18O value during this period (0·22{per thousand}; Table 3) would require ~5% bulk contamination, if the lowest {delta}18O value from the KS-3 well (1· 86{per thousand}) is used. The results from these two assimilation models are consistent and suggest that the extent of the contamination was small but significant for the 1924 event. The larger drop in matrix {delta}18O value following the 1868 collapse (0·45{per thousand}; Table 3) requires more extensive contamination (~12%) using the same assumptions. Unfortunately, no Pb, Sr or Nd isotope data are available for the 1877 sample to test this hypothesis. Additional testing of contamination scenarios using rock major or trace element data is complicated by the probability that the contaminants were only partially digested (e.g. 1924 eruption; Pietruszka & Garcia, 1999aGo). Mauna Loa rocks have lower incompatible elements than Kilauea rocks, although during partial melting and assimilation, the incompatible element concentrations of the resulting melts will be increased relative to the host Kilauea magma. An alternative process to bulk contamination for the summit rocks with more typical Kilauea geochemistry (higher Pb and lower Sr isotope ratios; e.g. sample 1954-1) is O isotope exchange with groundwater or altered lavas, which would leave no geochemical signature other than lowering the {delta}18O value (e.g. Kilauea groundwater and steam condensates have {delta}18O values of –1 to –14{per thousand}; Hinkley et al., 1995Go). Regardless of the contamination source, the rapid change in matrix {delta}18O values suggests that Kilauea's shallow magma reservoir is small, which is consistent with some geophysical estimates (e.g. Decker, 1987Go) and geochemical modeling (Pietruszka & Garcia, 1999bGo). The Kilauea situation is in marked contrast to the consistency of matrix {delta}18O values for Grimsvötn volcano over the last 800 years. The consistency of Grimsvötn's matrix {delta}18O values was related to its large crustal magma chamber, which is supported by the 15 km3 flood basalt flow it produced in 1783 (Bindeman et al., 2006Go). In contrast, Kilauea's summit reservoir is thought to have been small (only 2–3 km3) over the last 200 years (Pietruszka & Garcia, 1999bGo). The crustal contamination of Kilauea magmas is superimposed on a longer-term geochemical variation that reflects source heterogeneity (Pietruszka & Garcia, 1999aGo). The {delta}18O composition of this source is discussed in the next section.

Kilauea source oxygen isotope values
The origin of the diverse {delta}18O values observed in mantle-derived basalts has been of considerable interest. Although a limited range of {delta}18O values typifies MORB glasses (5·4–5·8{per thousand}), {delta}18O values correlate with source geochemical parameters (Eiler 2001Go) suggesting a mantle origin for at least some of the observed {delta}18O diversity in mantle-derived basalts (Harmon & Hoefs, 1995Go). However, concern has surrounded the origin of lower O isotopes (< 5·4{per thousand}, below normal MORB values), which are observed in many ocean island basalts. Some studies relate lower values to crustal contamination (e.g. Garcia et al., 1998Go; Thirlwall et al., 1997Go), whereas others have suggested they are an inherent feature of the mantle (e.g. Garcia et al., 1993Go; Harmon & Hoefs, 1995Go). For example, Wang et al. (2003Go) argued for a ‘primary’ mantle {delta}18O value of 5·5{per thousand} for HSDP2 Mauna Kea lavas, with lower values resulting from contamination (e.g. Wang et al., 2003Go). The new results presented here for historical Kilauea summit lavas support the contamination explanation for lower {delta}18O values.

The consistency of the olivine {delta}18O values in historical Kilauea summit lavas (5·1 ± 0·1{per thousand} ; Fig. 3), despite the significant range in Pb and Sr isotopes during this period (Table 3), is strong evidence in favor of the interpretation that there is a limited range in {delta}18O values in the mantle source (5·45–5·75{per thousand}; the value depends on the magnitude of the equilibrium matrix–olivine fractionation), despite the heterogeneous source producing Kilauea magmas (e.g. Pietruszka & Garcia, 1999aGo). This range is consistent with the prevailing value for MORB mantle (Ito et al., 1987Go; Eiler 2001Go), suggesting that the {delta}18O value of the Hawaiian mantle plume feeding Kilauea is typical of the upper mantle. However, it should be noted that diversity exists within the plume. In particular, olivines from some Ko’olau volcano lavas have much higher {delta}18O values (5·5–6·0 vs 5·0–5·2{per thousand}), which are inferred to result from melting eclogite blocks within the Hawaiian plume (Eiler et al., 1996Go).


    SUMMARY
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
A wide range in matrix {delta}18O value (4·9–5·6{per thousand}) was determined for Kilauea summit lavas spanning ~310 years of eruptive history. These results are in marked contrast to the essentially constant but much lower matrix {delta}18O values (3·1 ± 0·1{per thousand}) reported for the last 800 years of eruptive activity for the Icelandic hotspot volcano Grimsvötn. Kilauea matrix values do not vary systematically with time, although lower values were found in rocks erupted following major summit collapses or eruptive hiatuses. The difference between the matrix {delta}18O values for these volcanoes probably reflects the size of the magma reservoir system, with the Kilauea system being much smaller. Olivine {delta}18O values in most of the Kilauea lavas are nearly constant (5·1 ± 0·1{per thousand}), whereas those from Grimsvötn volcano show substantial variation (2·3–5·2{per thousand}). The disequilibrium between matrix and olivine {delta}18O values in Kilauea samples indicates that the lower matrix values were acquired by the magma after olivine grew, probably just before or during eruption. The sources of the contamination of Kilauea historical lavas are probably the host rocks for the magma reservoir system. These include Mauna Loa and Kilauea rocks, based on O, Pb and Sr isotope data, whereas oceanic crust contamination is not indicated. The extent of crustal contamination of Kilauea historical magmas was probably minor (~5%) for the 1924 collapse, depending on the assumed contaminant. The contamination signature of Kilauea magma is superimposed on a larger, longer-term cyclic geochemical variation that reflects source heterogeneity. Kilauea's heterogeneous source, which is well represented by the analyzed summit lavas, probably has {delta}18O values within the normal MORB mantle range (5·4–5·8{per thousand}), based on the new olivine {delta}18O values.


    ACKNOWLEDGEMENTS
 
This paper is dedicated to David Green on the occasion of his 70th birthday. David has been an inspiration to all who study basaltic magmas. Thanks go to the many friends who have assisted with this study, including Sorena Sorensen for lending critical samples from the Smithsonian collection; Nancy Baker, Rachel Konishi, Joann Romano, Kate Bridges, Kristina Garcia and Joe Horrell for help in the field; Chad Shishido and Kelly Kolysko for sample preparation; to Rhea Workman for assistance with the laboratory work; and Reed McEwan and John Tacinnelli for matrix sample oxygen isotope analysis. Analytical work in the Caltech stable isotope laboratories was supported, in part, by NSF grant EAR-0337736 and a grant from the EAR Technician support program. Our thanks go to journal reviewers Amy Gaffney, J. M. Rhodes and Ilya Bindeman for their helpful comments. This paper is SOEST Contribution 7149, and was supported by NSF grants (EAR03-36874) to MG and (EAR-0345905) to JE.


*Corresponding author. E-mail: morgarcia{at}hawaii.edu


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 KiLAUEA ERUPTIVE HISTORY SUMMARY
 SAMPLES
 METHODS AND RESULTS
 DISCUSSION
 SUMMARY
 REFERENCES
 
Bindeman IN, Sigmarsson O, Eiler JM. Time constraints on the origin of large volume basalts derived from O-isotope and trace element mineral zoning and U-series disequilibria in the Laki and Grimsvotn volcanic system. Earth and Planetary Science Letters (2006) 245:245–259.

Borthwick J, Harmon RS. A note regarding ClF3 as an alternative to BrF5 for oxygen isotope analysis. Geochimica et Cosmochimica Acta (1982) 46:1665–1668.[CrossRef][ISI]

Brigham WT. The Volcanoes of Kilauea and Mauna Loa on the Island of Hawaii, Their Variously Recorded History to the Present Time. In: Memoirs of the Bernice Pauahi Bishop Museum (1909) 2. Honolulu, HI: Bishop Museum Press. [Reprinted 1974, Kraus Reprint Company, Millwood, NY, 608 pp.].

Casadevall TJ, Dzurisin D. Stratigraphy and petrology of the Uwekahuna Bluff Section, Kilauea Caldera. US Geological Survey, Professional Papers (1987) 1350:351–375.

Clayton RN, Mayeda TK. The use of bromine pentafluoride in the extraction of oxygen from oxides and silicates for isotopic analysis. Geochimica et Cosmochimica Acta (1963) 27:43–52.[CrossRef][ISI]

Cooper K, Eiler JM, Asimow P, Langmuir C. Oxygen isotope evidence for the origin of enriched mantle beneath the mid-Atlantic ridge. Earth and Planetary Science Letters (2004) 220:297–316.[CrossRef][ISI]

Coplen TB. Normalization of oxygen and hydrogen isotope data. Chemical Geology (Isotope Geoscience) (1988) 72:293–297.[CrossRef]

Dana JD. Characteristics of Volcanoes (1891) New York: Dodd, Mead. 399.

Decker RW. Dynamics of Hawaiian volcanoes: an overview. US Geological Survey Professional Paper (1987) 1350:997–1018.

Dvorak JJ, Dzurisin D. Variations in magma supply rate at Kilauea Volcano, Hawaii. Journal of Geophysical Research (1993) 98:22255–22268.

Eiler JM. Oxygen isotope variations of basaltic lavas and upper mantle rocks. In: Stable Isotope Geochemistry. Mineralogical Society of America, Reviews in Mineralogy—Valley JW, Cole DR, eds. (2001) 43:319–364.

Eiler JM, Farley K, Valley JW, Stolper EM, Hauri E, Craig H. Oxygen isotope evidence against bulk recycled sediment in the source of Pitcairn Island lavas. Nature (1995) 377:138–141.[CrossRef]

Eiler JM, Farley KA, Valley JW, Hofmann A, Stolper EM. Oxygen isotope constraints on the sources of Hawaiian volcanism. Earth and Planetary Science Letters (1996) 144:453–468.[CrossRef][ISI]

Eiler JM, Schiano P, Kitchen N, Stolper EM. Oxygen isotope evidence for recycled crust in the sources of mid-ocean ridge basalts. Nature (2000) 403:530–534.[CrossRef]

Fekiacova Z, Abouc