UBC Faculty Research and Publications

Hawaiian hot spot dynamics as inferred from the Hf and Pb isotope evolution of Mauna Kea volcano. Weis, Dominique 2003

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata


52383-Weis_AGU_2002GC000340.pdf [ 1.34MB ]
JSON: 52383-1.0075975.json
JSON-LD: 52383-1.0075975-ld.json
RDF/XML (Pretty): 52383-1.0075975-rdf.xml
RDF/JSON: 52383-1.0075975-rdf.json
Turtle: 52383-1.0075975-turtle.txt
N-Triples: 52383-1.0075975-rdf-ntriples.txt
Original Record: 52383-1.0075975-source.json
Full Text

Full Text

Geochemistry Geophysics Geosystems  3  G  Article Volume 4, Number 2 19 February 2003 8704, doi:10.1029/2002GC000340  AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society  ISSN: 1525-2027  Hawaiian hot spot dynamics as inferred from the Hf and Pb isotope evolution of Mauna Kea volcano Janne Blichert-Toft Ecole Normale Supe´rieure de Lyon, 46 Alle´e d’Italie, 69364 Lyon, Cedex 7, France ( jblicher@ens-lyon.fr)  Dominique Weis De´partement des Sciences de la Terre et de l’Environnement, Universite´ Libre de Bruxelles, 1050 Brussels, Belgium Now at Earth and Ocean Sciences, University of British Columbia, Vancouver V6T 1Z4, Canada.  Claude Maerschalk De´partement des Sciences de la Terre et de l’Environnement, Universite´ Libre de Bruxelles, 1050 Brussels, Belgium  Arnaud Agranier and Francis Albare`de Ecole Normale Supe´rieure de Lyon, 46 Alle´e d’Italie, 69364 Lyon, Cedex 7, France  [1] The present work reports multiple collector inductively coupled plasma mass spectrometry (MC-ICPMS) measurements of the isotopic compositions of Hf and Pb in the first 3 km of the deep core retrieved by the Hawaii Scientific Drilling Project. The measurements cover all the samples from the standard geochemical reference set, glasses from the deep hole, and replicates from the pilot hole. Both Hf and Pb are less radiogenic in Mauna Loa compared to Mauna Kea. The transition between Mauna Kea and Mauna Loa lavas in the deep core is progressive for eHf and 208Pb/204Pb, but a sharp discontinuity is observed for 208 Pb*/206Pb*. There is no correlation between the alkalinity of the samples and isotopic composition. In detail, the Hf isotope compositions of samples from the pilot hole are not all identical to those of the HSDP-2 core for samples retrieved from a similar depth, suggesting that steep topography existed at the time of emplacement or that a different eruptive sequence was recorded. The strong correlation between 208 Pb*/206Pb* and 3He/4He (He data from M. D. Kurz et al. (Rapid helium isotopic variability in Mauna Kea shield lavas from the Hawaiian Scientific Drilling Project, submitted to Geochemistry Geophysics Geosystems, 2002)) requires the episodic incorporation of a component that resembles the basalts erupted by either Kilauea or the Loihi eruptive centers (this component is referred to as K/L). The data suggest that some 500 kyr ago, Mauna Kea was tapping a mantle source similar to that tapped by Kilauea today. Isotopic variability of Pb and He cannot be accounted for by radiogenic ingrowth in a closed system, but requires the mixing of mantle source components with distinct outgassing histories. The time series of isotopic and concentration data in Mauna Kea samples spanning about 350,000 years of age indicate the recurrence of geochemical patterns in the melting column. Ignoring the most recent alkalic samples, we find that the dominant fluctuations of eHf and 207Pb/204Pb correspond to a period of 50,000 years. For La/ Yb, Zr/Nb, 87Sr/86Sr, 206Pb/204Pb, 207Pb/206Pb, and 208Pb/206Pb, a dominant period of ca. 18,000 years is obtained. Once provision is made for the existence of harmonics, the consistency between the isotopic spectrum of the pilot hole and the HDSP-2 core is very good. The input of the K/L component does not seem to be periodic. We use these recurrence intervals in conjunction with the upwelling rate deduced from buoyancy flux and seismic evidence of the maximum dimension of scatterers to constrain the radius of the Hawaiian plume conduit to be in the range of 10–50 km and the upwelling velocity to be in the range of 0.13–3 m/yr. Plausible vertical length scales of heterogeneities in the conduit are 6.5–160 km.  Copyright 2003 by the American Geophysical Union  1 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Components: 12,431 words, 16 figures, 3 tables. Keywords: HSDP; Hawaii; Mauna Kea; Mauna Loa; Kilauea; Loihi; Hf isotope ratios; Pb isotope ratios; oceanic island basalts; isotope evolution of the mantle; mantle plumes; mantle heterogeneity; time series. Index Terms: 1025 Geochemistry: Composition of the mantle; 1040 Geochemistry: Isotopic composition/chemistry; 1094 Geochemistry: Instruments and techniques; 8120 Tectonophysics: Dynamics of lithosphere and mantle—general; 8121 Tectonophysics: Dynamics, convection currents and mantle plumes. Received 5 March 2002; Revised 18 November 2002; Accepted 18 November 2002; Published 19 February 2003. Blichert-Toft, J., D. Weis, C. Maerschalk, A. Agranier, and F. Albare`de, Hawaiian hot spot dynamics as inferred from the Hf and Pb isotope evolution of Manua Kea volcano, Geochem. Geophys. Geosyst., 4(2), 8704, doi:10.1029/2002GC000340, 2003. ————————————  Theme: Hawaiian Scientific Drilling Project  Guest Editors: Don DePaolo, Ed Stolper, and Don Thomas  1. Introduction [2] Although numerous important scientific results were obtained from the geochemistry of the HSDP pilot hole lavas, this relatively short core did not make it possible to use the geochemical analyses, in particular those of radiogenic isotopes, as time series because the time window was not sufficiently broad to be of geophysical significance with respect to hot spot dynamics [e.g., DePaolo et al., 2001]. In the present contribution to the HSDP initiative, we use a comprehensive time series of combined highresolution, high-precision Hf and Pb isotope data for the continuous 3 km sequence of Mauna Kea lavas to set constraints on the nature of the Mauna Kea mantle sources and how they varied through time. Before reaching Mauna Kea by drilling, ca 250 m of lavas from the currently active Mauna Loa volcano had to be traversed, and here we also present the Hf and Pb isotope data for this suite of samples. The petrology and geochemistry of modern Hawaiian lava flows offer a two-dimensional view of the hot spot activity which has been used to infer cross sections of melt extraction [Ribe and Christensen, 1999] and mantle heterogeneities in the plume conduit [DePaolo et al., 2001]. One of the most important goals of a high-resolution series of samples covering a significant interval of the lifetime of a volcano is to provide geochemical time series (assuming a reliable timescale for the stratigraphic column) that allow the geochemical cross sections inferred from the modern eruptions to be tested. The present results for the Mauna Kea and Mauna Loa  volcanoes have allowed us to deduce aspects of the chemical and physical structure and the dynamics of the Hawaiian plume beneath Mauna Kea, examine how the plumbing systems of different volcanoes may interact, and make inferences about the relationships of Mauna Kea with its neighboring volcanoes, notably Mauna Loa, Kilauea, and Loihi.  2. Samples and Analytical Techniques 2.1. Hf Isotope Analyses [3] We report 155 Hf isotope analyses of 123 different samples from Mauna Loa (19 samples) and Mauna Kea (104 samples) from the second Hawaii Scientific Drilling Project core (HSDP-2) (Table 1). The Hf isotopic compositions were analyzed by MC-ICP-MS (the VG Plasma 54 instrument in Lyon) during the period of January 2000 to April 2001. Hafnium purification and mass spectrometry followed the procedures described by Blichert-Toft et al. [1997]. To avoid ambiguity between measured isotopic variations and potential chemistry problems or machine drift during long run sessions (as required by the large number of samples analyzed), samples were routinely randomized for both Hf separation chemistry and mass spectrometric analysis. Furthermore, in order to monitor machine performance, the JMC-475 Hf standard was run systematically in between every second or third sample throughout the entire project. Hf total procedural blanks routinely were less then 25 pg. 2 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Table 1. (Representative Sample) Hf and Pb Isotope Compositions of Mauna Loa and Mauna Kea Samples From the HSDP-2 Core [The full Table 1 is available in the HTML version of this article at http://www.g.cubed.org]. 176  Hf/177Hf  eHf  Sample Name  mbsl  SR0008-2.70  9.5  SR0023-2.90  34.0  SR0031-0.50  45.5  SR0036-1.22  53.5  SR0040-1.07  59.5  SR0046-1.15  68.7  0.283051 0.283059 0.283058 0.283064 0.283077 0.283083 0.283086 0.283085 0.283070 0.283087 0.283091  SR0057-3.00 SR0061-0.00 SR0066-0.00  85.9 91.2 98.8  0.283087 ± 6 0.283087 ± 7 0.283096 ± 8  11.1 11.1 11.5  SR0080-0.35 SR0083-7.85 SR0089-1.15 SR0093-6.55 SR0098-2.00  125.4 137.0 149.9 163.3 177.8  SR0104-4.95  197.4  0.283084 ± 6 0.283088 ± 7 0.283087 ± 7 0.283088 ± 5 0.283108 ± 6 0.283115 ± 10 0.283074 ± 6  11.0 11.2 11.1 11.2 11.9 12.1 10.7  SR0109-0.65 SR0113-6.20 SR0117-4.00 SR0120-1.00  209.1 222.5 233.7 242.0  0.283067 0.283092 0.283078 0.283080  10.4 11.3 10.8 10.9  SR0121-4.40  246.2  SR0124-3.90 SR0125-6.25  252.9 256.5  0.283103 ± 5 0.283112 ± 10 0.283100 ± 5 0.283097 ± 5  SR0127-4.75 SR0129-5.20  261.7 267.5  SR0131-6.92  274.4  SR0133-8.20 SR0137-5.98 SR0141-7.90  ± ± ± ± ± ± ± ± ± ± ±  ± ± ± ±  ± ± ± ± ± ± ±  6 7 8 7 9 7 7 8 9 7 6  5 5 6 8  Mauna Loa 9.9 10.1 10.1 10.3 10.8 11.0 11.1 11.1 10.5 11.1 11.3  Mauna Kea 11.7 12.0 11.6 11.5  281.3 293.0  0.283104 0.283126 0.283126 0.283125 0.283123 0.283121 0.283109  4 7 7 5 8 5 5  11.7 12.5 12.5 12.5 12.4 12.3 11.9  305.8  0.283114 ± 5  12.1  206  Pb/204Pb  207  Pb/204Pb  208  Pb/204Pb  18.0923 18.0978 18.0841 18.0728 18.0948  15.4621 15.4673 15.4701 15.4463 15.4592  37.8526 37.8667 37.9204 37.8387 37.8257  18.1248 18.1174 18.1152 18.1262 18.0831 18.0950 18.1662 18.1332 18.2433 18.2316 18.2230 18.2063 18.1591 18.2469 18.2516  15.4827 15.4580 15.4698 15.4683 15.4498 15.4514 15.4647 15.4591 15.4742 15.4635 15.4644 15.4575 15.4605 15.4603 15.4666  37.9123 37.8290 37.8704 37.8483 37.8039 37.7984 37.8855 37.8716 37.9572 37.9419 37.9462 37.9125 37.9147 37.9569 37.9674  18.0803 18.0614 18.1263 18.1798 18.2166 18.2505 18.2357  15.4812 15.4565 15.4561 15.4647 15.4765 15.4564 15.4531  37.9174 37.8375 37.8824 37.9237 37.9954 37.9407 37.9291  18.4111 18.4179 18.4081 18.4076 18.4329 18.4310 18.4362  15.4825 15.4822 15.4774 15.4936 15.4888 15.4849 15.4794  37.9959 37.9853 37.9487 38.0349 38.0079 38.0095 37.9970  18.4413  15.4863  38.0281  18.4301 18.4339 18.4343 18.4465  15.4834 15.4794 15.4779 15.4832  37.9918 37.9807 37.9664 38.0129  mbsl, meters below sea level. Hf and Pb isotopic compositions measured by MC-ICP-MS (VG Plasma 54). 176 Hf/177Hf normalized for mass fractionation to 179Hf/177Hf = 0.7325. 176Hf/177Hf of JMC-475 Hf standard = 0.28216 ± 0.00001 (i.e., 2s external reproducibility = 0.0035% or 35 ppm as also seen from the p sample replicate analyses). Hf standard run every second or third sample. Uncertainties reported on Hf measured isotope ratios are in-run 2s/ n analytical errors in last decimal place, where n is the number of measured isotopic ratios. NBS 981 Pb standard run after every two samples. The 2s standard error on NBS 981 on a typical day of analysis was 142, 120, and 125 ppm, respectively, for 206Pb/204Pb, 207Pb/204Pb, and 208Pb/204Pb. Average individual errors on Pb isotopic analyses (two standard deviations) were 0.0048 for 205Tl/203Tl, 0.0204 for 208Pb/204Pb, 0.0186 for 207Pb/204Pb, and 0.0190 for 206Pb/204Pb. Reproducibility of replicate analyses, including leaching, was 0.6 permil (relative difference) or 0.025 for 208Pb/204Pb, 0.008 for 207Pb/204Pb, and 0.010 for 206Pb/204Pb. a Hand-picked glass as opposed to the rest of the samples, which are powders from the standard geochemical reference sample set. b Original analyses published by Blichert-Toft and Albare`de [1999]. Numbers in parentheses correspond to original values.  3 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  [4] The analyzed sample suite spans the entire depth interval from sea level to current coring depth with the deepest sample coming from 3069 mbsl. One hundred and nineteen of these samples are powders from the 119-sample standard geochemical reference set, while four samples are handpicked glasses from Mauna Kea from the depth interval between 2200 and 2400 mbsl. Based on major element data (J. M. Rhodes and M. J. Vollinger, Composition of basaltic lavas sampled by phase 2 of the Hawaiian Scientific Drilling Project: Geochemical stratigraphy and magma types, submitted to Geochemistry Geophysics Geosystems, 2002; hereinafter referred to as J. M. Rhodes and M. J. Vollinger, submitted manuscript, 2002), this particular depth interval may represent a unique excursion of increased alkalinity in the middle of several kilometers of basalts with a strong tholeiitic character. We considered it important to investigate this feature in more detail by undertaking the additional glass analyses. Nineteen of the samples from the reference suite are from Mauna Loa (sea level down to 246 mbsl) and 100 of the samples are from Mauna Kea (between 246 and 3098 mbsl, the coring depth at the time of writing). [5] In order to verify the external reproducibility for Hf isotopic measurement of 0.35 epsilon units, as indicated by the JMC-475 Hf standard analyses, six samples from Mauna Loa and 20 samples from Mauna Kea were analyzed in duplicate. Some of these samples were chosen for repeat measurement because their in-run errors were greater than the external reproducibility as judged from the analyses of JMC-475. Three additional samples from Mauna Kea were further analyzed in triplicate in order to certify their measured Hf isotopic composition. The excellent reproducibility of all replicate analyses (Table 1) and the consistency of Hf isotopic compositions between consecutive samples, which, as noted above, were not analyzed in successive order, indicate that our error assignment is conservative. In order to further check for potential analytical bias between Hf isotopic measurements for HSDP-1 (the first Hawaii Scientific Drilling Project pilot core) and HSDP-2, four samples from the bottom of HSDP-1 (which overlaps with HSDP-2) were reanalyzed simultaneously with samples from HSDP-2  (Table 1). The results for the HSDP-1 replicates are identical within error to the original analyses carried out by Blichert-Toft and Albare`de [1999], thus allowing direct comparison between the two cores.  2.2. Pb Isotope Analyses [6] All the samples measured for Pb isotopic composition were processed in the clean laboratory at Universite´ Libre de Bruxelles. The samples were leached with 6 N HCl, as described by Weis and Frey [1991], in order to remove a maximum of secondary alteration products. Lead was subsequently separated by a single passage through an anion-exchange column. Total blanks for the entire procedure were 100–200 pg, which is negligible with respect to the Pb of the samples analyzed ($100 – 300 ng). The Pb isotopic compositions were measured on the MC-ICP-MS VG Plasma 54 in Lyon from June 2000 to June 2001. Thallium spike was added to all samples and standards to monitor mass fractionation [White et al., 2000]. The total Pb ion beam intensities varied from 0.6 to 8 V, but only results for samples with an intensity exceeding 2 V of total Pb were used in order to have sufficient precision on the isotopic ratios. Furthermore, because the Tl/Pb ratio appears to be critical to the accuracy of the Pb isotope results, especially for 207Pb/204Pb, only results for samples with Tl/Pb < 0.5 were retained. Almost all the samples from the standard geochemical reference set were analyzed. Among these, 19 are duplicates and two are triplicates. Nineteen samples are from Mauna Loa and 95 samples are from Mauna Kea. As done for Hf, samples were processed randomly for both Pb separation chemistry and mass spectrometric analysis. Because the NBS 981 standard is run systematically after every two samples as a reference to calculate the sample isotopic ratios, no single value of this standard can be reported that would be an independent test of machine accuracy. However, it has been shown by White et al. [2000] that the method used in this study reproduces thermal-ionization Pb isotope results on basaltic samples very accurately. It is worth reemphasizing here that the current technique does not assume that Pb and Tl have similar mass biases, but rather that the biases of these two elements maintain a constant 4 of 27  3  Geochemistry Geophysics Geosystems  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  0  (b)  (a)  Mauna Loa 200 400  Mauna Kea  600 800 1000  subaerial submarine  1200  Depth (mbsl)  1400 1600  triple-spike data  1800 2000 2200 2400 2600 2800 3000 3200 18.0  18.1  18.2  18.3  18.4  18.5  18.6  18.7  206Pb/204Pb  15.44 15.45 15.46 15.47 15.48 15.49 15.50 15.51 15.52 207Pb/204Pb  Figure 1. Downhole Pb isotopic variations for HSDP-2 as measured by the thallium-spike method using MC-ICPMS (Plasma 54). Except for 207Pb/204Pb, the analytical uncertainties are smaller than the size of the symbols. Overall good overlap and comparable down core variations with the triple-spike TIMS data from J. Eisele et al. (submitted manuscript, 2002) (open circles) are demonstrated between the two data sets with the exception of 207Pb/204Pb. 208 Pb*/206Pb* ratios were calculated by removing primordial lead from the 206Pb/204Pb and 208Pb/204Pb ratios.  ratio over the time period of a few runs. The good performance of the Plasma 54 was further substantiated by the measurement of multiple NIST 981NIST 982 mixtures (B. Nelson and F. Albare`de, unpublished data): 207Pb/206Pb and 208Pb/206Pb are precise to better than 50 ppm, while 206Pb/204Pb, 207 Pb/204Pb, and 208Pb/204Pb typically are precise to within 100 ppm per atomic mass unit. [7] The Pb isotope results and a comparison with triple-spike Pb isotope data (J. Eisele et al., The long-term Pb isotope evolution of Mauna Kea lavas from the HDSP-2 drill core, submitted to Geochemistry Geophysics Geosystems, 2002; hereinafter referred to as J. Eisele et al., submitted manuscript, 2002) on a subset of samples from the Mauna Kea section of the deep core (HSDP-2) is presented in  Figure 1. Except for 207Pb/204Pb, for which the Plasma 54 data are clearly significantly more noisy than the triple-spike data, the agreement between the MC-ICP-MS and the triple-spike values is remarkable, especially when taking into account the down core isotopic variations. The cause for the noise on the 207Pb/204Pb measurements by MC-ICP-MS is not yet understood. Although matrix left after elution may alter mass fractionation to some extent, it is our experience that Pb isotope measurements on solutions with a relatively heavy matrix reproduce well within the stated errors. Memory should not be a significant issue since the Plasma 54 of Lyon never sees any 204Pb or 207Pb spike, while the 207Pb/204Pb ratios of the bracketing standard (NIST 981) and the samples are nearly identical. J. Eisele et al. (submitted manuscript, 2002) observe a lack of reprodu5 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  0  (c)  (d)  200 400 600 800 1000  subaerial  1200  submarine  Depth (mbsl)  1400 1600 1800 2000 2200 2400 2600 2800 3000 3200 37.7  37.8  37.9  38.0  38.1  38.2  0.925  0.930  0.935  208Pb/204Pb  0.940  0.945  0.950  0.955  0.960  208Pb*/206Pb*  Figure 1. (continued)  cibility of 207Pb/204Pb measurements on some subaerial basalts from the HSDP pilot core and ascribe it to contamination by particles transported with groundwater. Because our leaching procedure is somewhat less severe than the procedure used in Mainz, we may be seeing a similar effect. A last option is that our HCl-based chemistry fractionates Pb isotopes on anion-exchange resin. Such a process has been previously documented for Zn, Cu [Mare´chal et al., 1999; Mare´chal and Albare`de, 2002] and Fe [Anbar et al, 2000] and has been tested here for Pb by eluting a few micrograms of NIST 981 on an anion-exchange column. The surprising result of this experiment (mass fractionation on 207Pb/204Pb in excess of 2000 ppm across the elution peak) is described in Appendix A. A similar experiment with similar outcome was described by Baker and Waight [2002]. Unfortunately, it has not been possible to assess the yield of our Pb purification chemistry for  the HSDP-2 samples, which were all run prior to this experiment and for which Pb concentrations are not precisely known. However, the observation that the Pb isotope partition coefficients between the anionexchange resin and HCl may change by as little as 40 ppm per mass unit from 204Pb to 208Pb may account for some of the noise observed in our HSDP-2 data. Natural variations in 206Pb/204Pb and 208Pb/204Pb are an order of magnitude larger than the 207Pb/204Pb variations: the impact of potential analytical artifacts is therefore less critical on the first two ratios.  3. Results 3.1. Hf Isotopes [8] Hf isotopic compositions for Mauna Loa lavas from the deep core (HSDP-2) have 176Hf/177Hf that 6 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  0 Mauna Loa Mauna Kea  200 400 600  Depth (mbsl)  800 1000  subaerial submarine  1200 1400  ML MK MK (glass)  1600 1800 2000  Hyaloclastite  2200 2400 2600  error  2800 3000 3200  9  10  11  ε  12  13  14  15  Hf  Figure 2. Variations in eHf with depth for Mauna Loa (ML) and Mauna Kea (MK) HSDP-2 lavas. Replicate analyses as listed in Table 1 are all shown individually.  varies from 0.283051 to 0.283115, corresponding to eHf = +9.9 À +12.1 (Table 1; Figure 2). This is similar to, though slightly more variable than observed for Mauna Loa lavas from the pilot core (HSDP-1) (eHf = +10.5 À +11.6 [Blichert-Toft and Albare`de, 1999]). Samples from the top of the Mauna Loa section are the least radiogenic, with successive samples becoming progressively more radiogenic with depth except for a small wiggle at about 125 mbsl and an abrupt excursion to less radiogenic values at about 200 mbsl. Samples following this excursion gradually become more radiogenic during the remaining 50 m leading to  the contact with Mauna Kea (Figure 2). The range of 176Hf/177Hf for the HSDP-2 Mauna Kea section (0.283087– 0.283152 or eHf = +11.1 À +13.4) likewise is indistinguishable from that of the HSDP-1 Mauna Kea section (eHf = +11.4 À +13.1 [Blichert-Toft and Albare`de, 1999]). However, between 500 and 750 mbsl, Hf in HSDP-2 is offset systematically with about 0.5 epsilon units toward more radiogenic values relative to HSDP-1 (Figure 3b). The discrepancy in Hf isotopic composition for this part of the overlapping depth interval between the two cores cannot be due to an analytical artifact as the four duplicate analyses 7 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  0  (a)  (b)  200 400 HSDP-1 600 800 1000 Mauna Loa  1200  Depth (mbsl)  Mauna Loa break 1400  Mauna Kea  1600  Mauna Kea breaks K/L group  1800 2000 2200 2400 2600 2800 3000 3200 -2.5  -2.0  -1.5  -1.0  -0.5  0.0  0.5  1.0  1.5  Alkalinity  9.0  9.5  10.0 10.5 11.0 11.5 12.0 12.5 13.0 13.5 14.0  ε  Hf  Figure 3. Compared variations in alkalinity (a) and eHf (b) with depth for Mauna Loa and Mauna Kea HSDP-2 lavas. Replicate analyses plotted as single mean values. Alkalinity from J. M. Rhodes and M. J. Vollinger (submitted manuscript, 2002). The K/L group of mid-depth Mauna Kea defined on the basis of its Pb isotopic composition (see text for explanation) is shown with a bull’s-eye symbol. Breaks in Hf isotopic composition that exceed error bars (five in total) are shown with open symbols in panel b. Panel b also traces with solid lines the Hf isotopic variation in the HSDP pilot hole [Blichert-Toft and Albare`de, 1999] for purpose of comparison between the two cores. Between 500 and 750 mbsl, eHf in HSDP-2 is higher by about 0.5 epsilon units relative to the pilot hole. Another difference between the two cores in the Mauna Kea section is observed in the depth interval leading up to the contact with Mauna Loa, where HSDP-2 displays Hf isotopic compositions intermediate between those of Mauna Kea and Mauna Loa, whereas, in the pilot hole, Hf becomes very radiogenic and shows a sharp discontinuity with respect to Mauna Loa. By contrast, the Mauna Loa Hf isotopic signatures are identical for the pilot core and the HSDP-2 core.  from HSDP-1 analyzed together with the new samples from HSDP-2 came out identical within error to the first set of analyses by Blichert-Toft and Albare`de [1999] (Table 1). It would of course have been preferable to duplicate the Hf isotope measurements in the section where the two cores are isotopically distinct. We consider, however, that the reproducibility of HSDP pilot core data acquired on the same samples three years apart and the overall consistency of the multiple replicates of HSDP-2 analyses make it unlikely that an analyt-  ical bias may account for the Hf isotope discrepancy between the two cores. Despite the two holes having been drilled just a few kilometers apart and the similarity of one particular alkaline flow in each core (SR131 and R177 as pointed out by S. Huang and F. A. Frey (Trace element abundances of Mauna Kea basalt from phase 2 of the Hawaii Scientific Drilling Project: Petrogenetic implications of correlations with major element content and isotopic ratios, submitted to Geochemistry Geophysics Geosystems, 2002; hereinafter referred 8 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  to as S. Huang and F. A. Frey, submitted manuscript, 2002)), it seems clear that stratigraphic correlations cannot be retrieved from the depth in the core alone. [9] As observed previously for the pilot core [Blichert-Toft and Albare`de, 1999], Hf isotope compositions of Mauna Kea lavas in the HSDP-2 core (eHf $ +12.3 ± 0.9) are overall more radiogenic than those of Mauna Loa lavas (eHf $ +11.0±0.9) (Figure 2). While the transition from Mauna Kea to Mauna Loa lavas in the pilot core is abrupt with a drop of about one eHf unit [BlichertToft and Albare`de, 1999], it is preceded in the HSDP-2 core by a conspicuous gradual decrease of almost two eHf units over the top 200 meters of Mauna Kea flows (Figure 2). Hafnium in the flows between 500 and 1500 mbsl generally seems to be more radiogenic than in the flows at the bottom of the core. Hf isotope compositions oscillate with only few observed breaks that exceed error bars (Figure 3b). The hyaloclastite sample SR0705-0.15 (at 1823.2 mbsl), analyzed in triplicate to verify its aberrant (with respect to neighboring samples) isotopic composition, is clearly anomalous (Figure 2) and will not be discussed further. The four handpicked glasses from the depth interval of 2200 to 2400 mbsl were leached and also analyzed for their Hf isotopic compositions. Their eHf values are indistinguishable within error from those of the corresponding whole rock powders from the standard geochemical reference set (Figure 2). The glasses define a distinct trend of increasingly radiogenic 176Hf/177Hf (or increasing eHf) with decreasing depth that overlaps the identical trend defined by the whole rock powders.  3.2. Pb Isotopes [10] Lead isotopic compositions for the Mauna Loa part of the deep core (HSDP-2) vary between 18.061 and 18.250 for 206Pb/204Pb and 37.804 and 37.995 for 208 Pb/ 204 Pb (Table 1; Figure 1). Whereas 206Pb/204Pb and 207Pb/204Pb do not show systematic variability down core, 208Pb/204Pb increases more or less systematically between 9.5 and 242 mbsl, the latter depth defining the contact with Mauna Kea. For the Mauna Kea section of HSDP-2, 206Pb/204Pb varies between 18.408 and  18.670 and 208 Pb/ 204 Pb between 37.933 and 38.248. In a 208Pb/204Pb (and 207Pb/204Pb) versus 206 Pb/204Pb plot, there is no overlap between the Mauna Loa and Mauna Kea data (Figure 4). The lavas in the upper part of the Mauna Kea section are characterized by slightly less radiogenic Pb isotopic compositions than typical for deeper Mauna Kea lavas and cluster in the least radiogenic part of the Mauna Kea field in a Pb-Pb diagram. [11] While 208Pb/204Pb, like 176Hf/177Hf, displays some continuity between Mauna Kea and Mauna Loa samples, the transition between the two volcanoes for 206 Pb/ 204 Pb, 207 Pb/ 204 Pb, and 208 Pb*/206Pb* is marked by a sharp discontinuity (Figure 1). Relative to Mauna Loa, Mauna Kea is characterized by more radiogenic Pb and a lower source Th/U (208Pb*/206Pb*) ratio, the latter showing an overall increase with depth. Nine samples from middepth, notably from the interval 1973.8 to 2123.7 mbsl and from 2321.6, 2467.3, 2503.5, and 2654.1 mbsl, which all have a less negative alkalinity than the bulk of Mauna Kea lavas (Figure 3a), have distinctly higher 208Pb*/206Pb* ratios. For reasons to be explained in the Discussion, we will hereafter refer to these samples as the K/L group. As mentioned above, in a 207 Pb/ 204Pb versus 206 Pb/204Pb plot (Figure 4b), the difference between Mauna Kea and Mauna Loa stands out, but in contrast with the results of Abouchami et al. [2000a] on the pilot core, no fine structure can be defined within the data. We ascribe this to the larger errors on 207Pb/204Pb associated with the thalliumspike method used here with respect to the triplespike method employed by Abouchami et al. [2000a]. Another reason may be a variable efficiency of Pb contamination removal by acid leaching (see J. Eisele et al., submitted manuscript, 2002). As also remarked above, the same contrast between Mauna Kea and Mauna Loa shows up in a 208 Pb/204Pb versus 206Pb/204Pb plot with the more alkaline K/L Mauna Kea samples plotting with unusually high 208Pb/204Pb for a given 206Pb/204Pb (Figure 4a). Despite the potential analytical limitations on 207Pb/204Pb, we used a 3-D rotated Pb projection to attempt to demonstrate that most Mauna Kea data fall within a plane, whereas the 9 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  38.30 38.25 38.20  Mauna Loa Mauna Kea K/L group  208Pb/204Pb  38.15 38.10 38.05 38.00 37.95 37.90 37.85 37.80  (a)  37.75  15.50  207Pb/204Pb  15.49  15.48  15.47  15.46  15.45  (b) 15.44 18.0  18.1  18.2  18.3  18.4  18.5  18.6  18.7  206Pb/204Pb  Figure 4. 208Pb/204Pb (a) and 207Pb/204Pb (b) versus 206Pb/204Pb for Mauna Loa and Mauna Kea HSDP-2 lavas. Note the distinct Pb isotopic compositions of Mauna Kea and Mauna Loa volcanoes in both Pb-Pb plots and of the K/L group in 208Pb/204Pb-206Pb/204Pb space.  more alkaline K/L samples plot distinctly above this plane (Figure 5).  4. Discussion 4.1. Interaction Between Mauna Kea Basalts and the Pacific Lithosphere [12] The field of combined Hf and Pb isotope compositions for the HSDP-2 samples is very sim-  ilar to that of the pilot hole and we will not repeat here the discussion presented elsewhere for Hf [Blichert-Toft and Albare`de, 1999; Blichert-Toft et al., 1999; S. Huang and F. A. Frey, submitted manuscript, 2002] and Pb [Abouchami et al., 2000a]. We note, however, that the present Pb and Hf isotope compositions plot on the general hyperbolic array that Blichert-Toft et al. [1999] used to infer that the East Pacific plate is not involved as a significant 10 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  K/L group Intermediate Mauna Kea Mauna Kea Mauna Loa  208Pb/204Pb  206Pb/204Pb  208Pb/204Pb  206Pb/204Pb  207Pb/204Pb  207Pb/204Pb  Figure 5. 3-D rotated projection of 206Pb/204Pb, 207Pb/204Pb, and 208Pb/204Pb shown in two different frozen positions that together illustrate the plane defined by the bulk of Mauna Kea samples and the displacement above this plane of the more alkaline K/L group. We have for better legibility identified intermediate Mauna Kea samples as transitional to the K/L group.  component in Hawaiian basalts (Figure 6). This point is discussed by Z. Wang et al. (An examination of the roles of plume, lithosphere, and edifice in Hawaiian shield-building magmatism based on oxygen isotope date from the second HSDP core, submitted to Geochemistry Geophysics Geosystems, 2002; hereinafter referred to as Z. Wang et al., submitted manuscript, 2002) who emphasize that the eHf of Pacific MORB and Mauna Kea basalts overlap. The histogram shown as an inset in Figure 6 does indeed show some overlap and it may be argued that some Mauna Kea basalts may have incorporated particularly infrequent E-type (enriched) lithospheric material, e.g., similar to MORB from the Siqueiros fracture zone [Sims et al., 2002]. As a whole, however, the Pacific plate remains a very unlikely depleted end-member of Hawaiian basalts.  4.2. Transition Between Mauna Loa and Mauna Kea [13] Although Hf and Pb isotopic data allow a clear distinction to be made between lavas from Mauna Loa and Mauna Kea, there is a tendency for the uppermost 200 meters of Mauna Kea lavas to show isotopic compositions intermediate to those of the two volcanoes (Figures 1–3). The series of 207 Pb/204Pb and 208Pb/204Pb ratios are virtually  continuous across the transition, while the series of 206Pb/204Pb ratios shows a gap but no discontinuity of the trend. In contrast, the 208Pb*/206Pb* ratios show both a gap and an abrupt discontinuity of the trend. Only the most recent (topmost) of these transitional lavas are alkaline (Figure 3a), which confirms the observation from the pilot core that there is no overall correlation between alkalinity and Hf (or Pb) isotopic composition [Blichert-Toft and Albare`de, 1999]. Hafnium in the most alkaline samples from the uppermost Mauna Kea flows in the HSDP-2 core is substantially less radiogenic than in the alkaline samples from the same depth interval in the pilot hole [Blichert-Toft and Albare`de, 1999] (Figure 3b). Were it not for 208 Pb*/206Pb*, the HSDP-2 core Hf and Pb isotope data would support a transitional regime between Mauna Kea and Mauna Loa indicative of a mixed source of interacting batches of magmas, probably reflecting the gradual invasion of the plumbing system of Mauna Kea by that of Mauna Loa. This observation contrasts with the sharp discontinuity between the two volcanoes observed in the Hf and Pb isotope data from the pilot hole [Blichert-Toft and Albare`de, 1999; Abouchami et al., 2000a]. The distinctive properties of the two cores are further illustrated by the significantly more radiogenic character of Hf in 11 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  16  14  12  Mauna Kea (this work) Koolau Haleakala Kilauea Hilina Mauna Loa (this work) Loihi Puu Oo Kahoolawe  10  90  Pacific MORB  8  number of samples  εHf  80  6  4  Mauna Kea  70 60 50 40 30 20 10  2  0 10  0 17.8  18.0  18.2  11  12  18.4 206  Pb/  13  14  15  εHf  18.6  16  17  18  19  18.8  20  19.0  204  Pb  Figure 6. Position of the new Pb and Hf isotopic data on the Hawaiian basalts’ mixing hyperbola of Blichert-Toft et al. [1999]. Although in general Hawaiian basalts require more than two components to account for their geochemical variations, the Hf-Pb correlation is exceptionally good and well curved. The end-members are best represented by the Koolau and Kilauea basalts. The hyperbolic array allows the eHf of the Kilauea end-member and the 206Pb/204Pb of the Koolau end-member to be well constrained. Inset: histogram of eHf of Pacific MORB [Chauvel and Blichert-Toft, 2001; Sims et al., 2002] (Blichert-Toft, unpublished data) and Mauna Kea (this work) which shows that, although EMORB are undoubtedly present occasionally along the East Pacific Rise, the Pacific lithosphere as a whole remains an unlikely end-member of Hawaiian basalts.  the lava flows of the HSDP-2 core between 500 and 750 mbsl relative to those of the pilot core. These discrepancies suggest that, at equivalent depths, the two drill holes did not sample the exact same lava piles, indicating either a lateral shift in stratigraphy between the two cores over the ca 2 km that separate them or the accumulation of different lava flow sequences at the two locations, perhaps caused by steep topography at the time of eruption. In contrast, Mauna Loa lavas have identical eHf versus depth relationships in the two cores. As for Pb isotopes, no samples from the pilot hole were analyzed by MC-ICP-MS, thus preventing direct comparison between the two  cores. However, the agreement between the Pb isotopic compositions obtained for HSDP-2 in the present study by MC-ICP-MS and those obtained by triple spike for the HSDP-1 core [Abouchami et al., 2000a] is good.  4.3. Mauna Kea Source Heterogeneity [14] In order to evaluate which isotopic properties may be correlated, we first ran our own Hf and Pb isotope data together with the Sr and Nd isotope data produced by J. G. Bryce et al. (manuscript in preparation, 2003) and the He isotope data produced by M. D. Kurz et al. (Rapid helium isotopic varia12 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  component 3: 13.0%  1.5 1 87Sr /86Sr  0.5  206Pb /204Pb  εHf  0  208Pb /204Pb 208Pb*/206Pb*  εNd  0.5  3He /4He  1  component 2: 24.2%  1.5 1 0.5 208Pb*/206Pb*  εNd 0 εHf  87Sr /86Sr 3He /4He 208Pb /204Pb  0.5  206Pb /204Pb  1 1.5 1.5  1  0.5  0  0.5  1  1.5  component 1: 46.3% Figure 7. Principal component analysis of isotopic compositions in Mauna Kea basalts. Data from this work for Pb and Hf, from J. G. Bryce et al. (manuscript in preparation, 2003) for Sr and Nd, and from M. D. Kurz et al. (submitted manuscript, 2002) for He. Two components account for more than 70% of the isotopic variability and three components for 84%. The eHf and 206Pb/204Pb axes essentially lie in the plane of the first two components, which justifies the binary mixing shown in Figure 6. Note also that He and 208Pb*/206Pb* are the only coherent systems. This does not reflect a consistency of the U-Th-Pb system with the U-Th-He system (see text) but rather the existence of two separate components with distinct degassing histories.  bility in Mauna Kea shield lavas from the Hawaiian Scientific Drilling Project, submitted to Geochemistry Geophysics Geosystems, 2002; hereinafter referred to as M. D. Kurz et al., submitted manuscript, 2002) on the same samples through a Princi-  pal Component Analysis (PCA) software. PCA is a method used to identify the orthogonal directions of maximum variability and is described in a number of different textbooks [e.g., Johnson and Winchern, 1992; Albare`de, 1995; Reyment and Jo¨reskog, 13 of 27  Geochemistry Geophysics Geosystems  3  G  24  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  (a)  3He/4He  0.8  20  (b)  Nb/Y  0.7  16  0.6  12 0.5 8 208Pb*/206Pb* 0.935  0.940  (c)  SiO2 49.5  0.945  Zr/Y 6.00  0.950  6.75  7.50  (d)  TiO2  K/L group Intermediate Mauna Kea Mauna Kea  8.25  3.75  48.0  3.00  46.5  2.25  45.0  1.50  SiO2  MgO 10  15  20  25  45.0  46.5  48.0  49.5  Figure 8. 3He/4He versus 208Pb*/206Pb* (a), Nb/Y versus Zr/Y (b), SiO2 versus MgO (c), and TiO2 versus SiO2 (d) for Mauna Kea HSDP-2 lavas. Note in particular the positive correlation between 3He/4He and 208Pb*/206Pb*. He isotope data from M. D. Kurz et al. (submitted manuscript, 2002) and major and trace element data from J. M. Rhodes and M. J. Vollinger (submitted manuscript, 2002).  1996]. The three components obtained for 206Pb/ 204 Pb, 207Pb/204Pb, and 208Pb/204Pb, respectively, account for 46.3, 24.2, and 13.0 percent of the variability. There is hardly any correlation between the Pb isotopic ratios and other isotope compositions, since, as demonstrated by Blichert-Toft et al. [1999], such mixing relationships are strongly nonlinear. One exception stands out, however: 208 Pb*/204Pb* and 3He/4He plot next to each other in the space of the first three components (Figure 7) which indicates that they are strongly correlated [Eiler et al., 1998]. We surmise that the odd position of 87Sr/86Sr along the axis of the third component may signal the occasional presence of a seawater component, either as a contaminant in the lava pile or as a real mantle component. The position of the eHf and 206Pb/204Pb axes in the plane of the first two components confirms that, although three compo-  nents are necessary to describe the overall isotopic variability of Hawaiian basalts, two components seem to be enough to describe the relationship between eHf and 206Pb/204Pb [Blichert-Toft et al., 1999]. [15] The systematics of Pb isotopic compositions, particularly in 208Pb/204Pb-206Pb/204Pb space (Figure 4a), require that more than two components be involved in the Mauna Kea lavas. As indicated in section 3, the middepth K/L samples as a group and the rest of the Mauna Kea samples have different Pb isotopic compositions. The K/L samples, which J. M. Rhodes and M. J. Vollinger (submitted manuscript, 2002) refer to as type-3 basalts, are readily identified by their high values of 208Pb/204Pb for a given value of 206Pb/204Pb, but also stand out by their less negative alkalinity, high Ca, Ti, and Nb 14 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  contents, and low silica (Figure 8). Although J. M. Rhodes and M. J. Vollinger (submitted manuscript, 2002) found these lavas very similar to Loihi basalts, Frey and Rhodes [1993] list these features as characteristic of lavas from Kilauea with respect to Mauna Loa. The 3He/4He ratios of the Mauna Kea samples reach values of up to 25 (M. D. Kurz et al., submitted manuscript, 2002) which would be more diagnostic of Loihi basalts if it was not because 3He/4He values in excess of 20 are also known from Kilauea [Kurz et al., 1996]. More extensive sampling of Kilauea may turn up even less radiogenic He. The Sr, Nd [Garcia et al., 1995], and Hf [Blichert-Toft et al., 1999] (Figure 6) isotope compositions of the two volcanoes are not distinguishable from each other. The isotopic ratio for which the overlap between Loihi and Kilauea data is minimal is 208Pb*/204Pb*, with Loihi values being higher. However, evidence from conventional TIMS data [Staudigel et al., 1984; Garcia et al., 1993, 1995] that these two volcanoes do not overlap in Pb isotope space is not very strong and a better case must await full publication of high-precision triple-spike Pb isotope data [Abouchami et al., 2000a]. By referring to a K/L component, it is not our intention to suggest that a common KilaueaLoihi component exists which lumps together the properties of these two distinct volcanoes but rather that, based on the available data, one volcano may be as appropriate a reference as the other. With respect to the chronological relationships discussed below, we simply point out here that assigning the K/L (type-3) basalts to a Loihi stage in the evolution of the Mauna Kea volcano requires strengthening of existing arguments. [16] 3He/4He usually is observed to be decoupled from nonvolatile isotopic systems [e.g., Valbracht et al., 1996] but Eiler et al. [1998] demonstrated that, in Hawaiian basalts, there is a strong correlation between 3He/4He and some combination of Pb isotopic ratios that we wish to elaborate further on here. As expected from the PCA results, there is a strikingly positive correlation between 3He/4He and 208Pb*/206Pb* for Mauna Kea (Figure 8a). The K/L samples represent the high-3He/4He and high-208Pb*/206Pb* end-member of this correlation. It can be shown by writing the equation for  radiogenic ingrowth that this correlation is not due to mere radiogenic ingrowth of 4He in a closed system: 4  He 3 He    4  ¼ 0   238  Á He U À l238U T 0 þ8 e À1 3 He 3 He T0 0 235  Á U À l235U T 0 e À1 þ7 3 He 0 232  Á Th À l232Th T 0 e À1 þ6 3 He 0  in which the subscript 0 refers to modern values, T0 is the age of the Earth and li stands for the decay constant for the nuclide i. Likewise, the radiogenic ingrowth of the 208Pb*/206Pb* ratio is given by 208  Pbà 206 Pbà  232   ¼ 0   Th el232Th T 0 À 1 238 U l238U T 0 À 1 0 e  Upon elimination of the Th/U ratio, the last two equations can be combined as: 4  He 3 He    4  238  Á He U À l238U T 0 ¼ 3 þ 3 e À1 He He 0 T0 0 2  3 208 à  l235U T 0 e À 1 Pb 6 7  l T Á 48 þ  þ 6 206 à 5 Pb 0 e 238U 0 À 1 238 U=235 U 7  0  This relationship indicates that radiogenic ingrowth in a closed system with a broadly uniform 238 3 U/ He ratio but a variable Th/U ratio would lead to a negative correlation between 208Pb*/206Pb* and 3He/4He, which is the opposite of what we observe for the Mauna Kea data (Figure 8a). The trend observed in this figure shows that it is not the variation of the Th/U ratios in the source that creates the observed correlation but mixing between two mantle components with a distinct fractionation/outgassing history: the relatively undegassed K/L component with high 3He/4He (low integrated 238U/3He) happened to have a high source Th/U while the relatively degassed component with a ‘‘normal’’ (7–15) 3He/4He (higher integrated 238U/3He) happened to have a lower Th/U ratio. We are therefore not observing a real correlation between 3He/4He and the U-Th-Pb isotopic system [Eiler et al., 1998] but rather a simple mixing relationship between the 3He/4He and the U/He ratios of two variably outgassed 15 of 27  Geochemistry Geophysics Geosystems  3  observed εHf  G  13  ing average runn 12  11 150  detrended εHf  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  200  250  300  350  400  450  500  550  600  200  250  300  350  400  450  500  550  600  1  0  -1  Age (ka) Figure 9. Principle of secular trend removal (example shown is for Hf isotopes). A running average with Gaussian weights is calculated for each point. In the present case (solid line in the upper panel), we used a standard deviation of 50,000 years. The top panel shows the residues after removal of the secular trend. This procedure tends to remove the low frequencies for which the periodic character cannot safely be established.  mantle components with different Th/U ratios. In other words, the conundrum of He isotope decoupling from nonvolatile systems is yet to be resolved. This decoupling, which may reflect the decoupled dynamics of volatiles and magmas [Valbracht et al., 1996], probably explains why 3 He/4He is so variable in the HSDP-2 samples whereas their eHf is relatively constant. [17] The bulk of Mauna Kea samples do not separate into two distinct groups, but rather form a trend that shows that each Mauna Kea sample contains different proportions of a Kilauea-like mantle source (as an example we have singled out the ‘‘intermediate Mauna Kea’’ group as shown  in Figure 8). The Kilauea ‘‘flavor’’ is dominant in the nine K/L samples that, according to the age model of D. J. DePaolo et al. (manuscript in preparation, 2003), should be 496–532 kyr old. Kilauea is the youngest aerial Hawaiian volcano. The oldest dated Kilauea sample is only 341 kyr [Quane et al., 2000], but, from the growth model established on Mauna Kea, DePaolo and Stolper [1996] inferred that Kilauea is probably as old as 600 kyr. Because essentially all middepth Kea samples contain some proportion of the K/L component, it is unlikely that the HSDP-2 core inadvertently sampled isolated Kilauea lava flows. As first suggested by Abouchami et al. [2000b], this rather indicates that the young Mauna Kea was 16 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Table 2. Most Significant Peaks of the Lomb Periodograms Before Removal of the Secular Trend Ratio  Frequency, kaÀ1  Power  Period T, ka  Confidence Level  Log2T  Data Source  La/Yb La/Yb eHf 87 Sr/86Sr 87 Sr/86Sr 87 Sr/86Sr Zr/Nb 206 Pb/204Pb 207 Pb/206Pb 207 Pb*/206Pb* La/Yb Zr/Nb 208 Pb/206Pb 206 Pb/204Pb 207 Pb/206Pb 207 Pb*/206Pb* 207 Pb/204Pb eHf eHf eNd Zr/Nb 208 Pb/204Pb Zr/Nb eHf 208 Pb*/206Pb* 3 He/4He eNd 208 Pb/204Pb 207 Pb/206Pb 207 Pb*/206Pb* alkalinity La/Yb  0.2241 0.2208 0.2130 0.2027 0.1994 0.1247 0.1202 0.0998 0.0998 0.0997 0.0935 0.0933 0.0522 0.0521 0.0520 0.0519 0.0208 0.0198 0.0087 0.0083 0.0066 0.0052 0.0025 0.0023 0.0014 0.0009 0.0001 0.0001 0.0001 0.0001 0.0001 0.0001  6.8 6.8 5.4 6.5 6.4 5.4 5.6 6.8 6.9 7.0 5.6 6.3 5.4 5.6 6.4 6.9 9.1 6.0 7.8 5.1 9.4 11.6 6.8 12.1 16.0 15.3 9.6 27.8 5.3 5.4 10.5 22.5  4.46 4.53 4.69 4.93 5.02 8.02 8.32 10.0 10.0 10.0 10.7 10.7 19.2 19.2 19.2 19.3 48.1 50.5 114.9 120.5 151.5 192.3 400 435 714 1,111 10,000 10,000 10,000 10,000 10,000 10,000  0.873 0.879 0.574 0.849 0.825 0.602 0.635 0.879 0.896 0.900 0.658 0.803 0.618 0.651 0.823 0.892 0.988 0.741 0.951 0.526 0.990 0.999 0.874 0.999 1.000 1.000 0.993 1.000 0.578 0.605 0.997 1.000  2.16 2.18 2.23 2.30 2.33 3.00 3.06 3.32 3.32 3.33 3.42 3.42 4.26 4.26 4.27 4.27 5.59 5.66 6.84 6.91 7.24 7.59 8.64 8.76 9.48 10.12 13.29 13.29 13.29 13.29 13.29 13.29  a  a S. Huang and F. A. Frey (submitted manuscript, 2002). b This study. c J. G. Bryce et al. (manuscript in preparation, 2003). d J. M. Rhodes and M. J. Vollinger (submitted manuscript, e  a b c c c d b b b a d b b b b b b b b d b d b b e b b b b d a  2002).  M. D. Kurz et al. (submitted manuscript, 2002).  tapping a mantle source similar to the source Kilauea is tapping today.  4.4. Time Series Analysis [18] With reliable age information, a relatively long core, such as the HSDP-2 core, becomes a prime source of information for the length scales of heterogeneities in the mantle and the timescales of melting processes. One limitation, however, is the possible presence of hiatuses in the volcanic activity at the location of the drill core. The survival time of volcano surfaces (resurfacing time) is, however, known to be a Poisson process. It has been demonstrated that the ages of Hawaiian surface elements are exponentially distributed with,  for Kilauea, a mean resurfacing time of about 500 years [Holcomb, 1987], an estimate that also seems to hold for Mauna Loa [Kauahikaua et al., 1995]. The HSDP-2 core counts about 1000 flow units, most of which erupted over a period of about 350 kyr, suggesting an average resurfacing time of 350 years, i.e., not very different from estimates for other volcanoes. Assuming time-invariant resurfacing with an exponential distribution of resurfacing times, the number of lava flows within a given time interval Dt is Poisson-distributed with, for Mauna Kea, a mean value of Dt/350 and a standard deviation of (Dt/350)1/2. An estimate of the maximum error on age assignment at the 1s level for the entire core is therefore on the order of (1000)1/2 * 350 % 10,000 years, an uncertainty smaller than 17 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Table 3. Most Significant Peaks of the Lomb Periodograms After Removal of the Secular Trend by Running a Gaussian Filter With a Standard Deviation of 50,000 Yearsa Ratio  Frequency, kaÀ1  Power  Period T, ka  Confidence Level  Log2T  La/Yb La/Yb eHf 87 Sr/86Sr 87 Sr/86Sr 208 Pb/206Pb 207 Pb/206Pb La/Yb Zr/Nb 206 Pb/204Pb 207 Pb/206Pb 207 Pb*/206Pb* 208 Pb/206Pb alkalinity 207 Pb/204Pb eHf alkalinity eHf 208 Pb*/206Pb* eHf Zr/Nb 208 Pb/204Pb  0.2276 0.2210 0.2131 0.2026 0.1995 0.1598 0.0632 0.0587 0.0585 0.0528 0.0526 0.0524 0.0523 0.0462 0.0208 0.0198 0.0174 0.0129 0.0104 0.0087 0.0068 0.0064  5.41 6.28 5.60 6.87 7.69 5.26 5.11 5.74 5.96 6.73 7.58 8.14 6.60 5.62 9.93 6.38 5.73 5.40 5.37 8.20 6.63 9.88  4.39 4.52 4.69 4.94 5.01 6.26 15.82 17.04 17.09 18.94 19.01 19.08 19.12 21.65 48.08 50.51 57.47 77.52 96.15 114.9 147.1 156.3  0.590 0.803 0.644 0.888 0.949 0.559 0.507 0.686 0.737 0.875 0.944 0.968 0.859 0.650 0.995 0.818 0.680 0.585 0.593 0.968 0.855 0.994  2.13 2.18 2.23 2.30 2.32 2.65 3.98 4.09 4.10 4.24 4.25 4.25 4.26 4.44 5.59 5.66 5.84 6.28 6.59 6.84 7.20 7.29  a  See Table 2 for the source of data.  that associated with K-Ar ages. Consequently, although basalt production by Mauna Kea probably declined over its lifetime, we can argue that random hiatuses are unlikely to introduce a major bias in the time analysis of the geochemical data. [19] The ages in the stratigraphic column have been established by fitting the existing Ar-Ar dates (W. D. Sharp et al., manuscript in preparation, 2003), which are few and have relatively large error bars. Because the samples are irregularly spaced down core, we cannot use standard FFT (Fast Fourier Transform) techniques, as interpolation is known to produce spurious peaks. Instead, we used the method of periodograms specifically designed for unevenly spaced data [Lomb, 1976; Scargle, 1982; Albare`de, 1995; Abouchami et al., 2000a] and which is nothing more but a fit of the data by a periodic function of prescribed frequency. For the purpose of not creating spurious peaks at very low frequencies, an appropriate background (i.e., the secular trend) must first be removed. This background was evaluated here by smoothing the data with a Gaussian filter, which is a running mean with weights normally distributed around the  point to be smoothed (Figure 9). The standard deviation of 50,000 years chosen here for this filter ensures there is no significant bias in the highfrequency range. Other values of the standard deviation were tested with no or only very little consequence to the well-constrained part of the spectrum, but our option was that peaks beyond 50,000 years are insufficiently constrained within the considered time window and were therefore ascribed to the secular trend. Once normalized to the standard deviation of the analyzed variable, the peak heights can be associated with significance levels allowing for the probability that the peak is not a false alarm to be estimated [Scargle, 1982]. [20] Tables 2 and 3 presents the results of the periodogram analysis before and after removal of the secular trend for the isotopic compositions of Hf, Pb, Nd, Sr (Nd and Sr data from J. G. Bryce et al. (manuscript in preparation, 2003)), and He (from M. D. Kurz et al., submitted manuscript, 2002), for the incompatible element ratios La/Yb (from S. Huang and F. A. Frey, submitted manuscript, 2002) and Zr/ Nb, and for alkalinity (two latter both from J. M. Rhodes and M. J. Vollinger (submitted manuscript, 18 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  18 16  εHf  14  207Pb/206Pb  12 10 117 (0.97)  8  19.0 (0.94)  50.5 (0.82) 77.5 (0.58)  6  4.7 (0.64)  Power  4 2  16  206Pb/204Pb  208Pb/204Pb  14 12 10  156 (0.90)  8 18.9 (0.87)  6  6.3 (0.56)  4 2 0 0  0.05  0.10  0.15  0.20  0.05  Frequency  0.10  0.15  0.20  0.25  (ka-1)  Figure 10. Periodograms of the Hf and Pb isotopic signal after removal of a secular trend by running a Gaussian filter with a 50,000 years standard deviation. The peaks are labeled with the period (frequency À1) and, in parentheses, the significance level (levels smaller than 50% are not shown).  2002)) in Mauna Kea flows. The results for Pb and Hf isotopic data after removal of the secular trend are shown in Figure 10. These results report the frequency at the peak, the power (height of the peak), the period (one over the frequency), and the confidence level. Only peaks with confidence levels in excess of 50 percent have been considered. As it is clear that some of the peaks are harmonics of lower frequency (longer period) fluctuations, we plotted the base-two logarithm of the peak heights obtained after applying the Gaussian filter (Figure 11). In this base, two peaks differing by one unit are harmonics of each other (multiplication by two). The strongest peaks characterizing the trace element ratios and alkalinity signal are less significant than those characterizing the isotopic ratios, which strongly suggests that melting blurs the signal of source heterogeneity. Two families of peaks emerge that  differ by a factor of four and therefore are double harmonics of each other. These two groups of peaks (one at 4500–5000 years, the other at 17,000 – 19,000 years) may not represent fundamental periods but the HSDP-2 core is currently too short to allow for the identification of potentially longer periods. These two groups involve the incompatible element ratios La/Yb and Zr/Nb, 87 Sr/ 86 Sr, 206 Pb/204Pb, 207Pb/206Pb, 208Pb/206Pb and at a lower confidence level 176Hf/177Hf. In addition to the stable 19,000 years period, the unfiltered data show strong 10 kyr peaks consistent with the 206Pb/204Pb peak found by Abouchami et al. [2000a] on the pilot core and by J. Eisele et al. (submitted manuscript, 2002) for the HSDP-2 samples, but this feature appears to be an intermediate harmonics that disappears after removal of the secular trend. Slower and very significant fluctuations are observed for 19 of 27  3  Geochemistry Geophysics Geosystems  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  256  8 7  128  >80 confidence >50 confidence ≈50 ka group  5  64 32  ≈19 ka group  16  4 3  8 ≈ 4.8 ka group  f  Pb  04  Pb /2  20 8  εH  Zr /N b  f  Pb *  εH  6  εH f al ka lin ity  Pb  04  06  Pb */ 2 0  20 8  20 7  Pb al ka lin ity 20 7 Pb /2  0  6 Pb Pb */ 2 06 Pb 20 * 8 Pb /2  Pb  04  Pb /2  20 7  Zr /N b  Pb /2  20 6  20 7  Pb /2  04  6  04  20 8  Pb /2  Sr /8  87  f  6  εH  Sr /8  87  La /Y b  1 Pb  0 Pb  2  Sr  1  Sr  4  La /Y b  2  La /Y b  period (ka)  ln2 period (ka)  6  Figure 11. Peaks of the periodograms for isotopic ratios, two critical incompatible element ratios, and alkalinity ordered by increasing period. A secular trend with a period of 50,000 years was removed (Tables 2 and 3). The vertical scale is the base-two logarithm of the period in ka so that harmonics of the same fundamental period differ by exactly one unit. It is highly likely that the groups at 4.8 ka and 19 ka are in a harmonic relationship. The intermediate harmonics at ca. 10,000 a is present in the unfiltered data [see also Abouchami et al., 2000a]. The 50 ka group is present in both the filtered and unfiltered data (see Tables 2 and 3 and Abouchami et al. [2000a]). It is highly significant but longer periods are only marginally sampled.  207  Pb/204Pb (48,000 years), probably in phase and harmonic relationships with those of Hf (50,500 and 115,000 years) and reproduce the strong peak observed at 48,000 years by Abouchami et al. [2000a] on the pilot core. The 208Pb*/206Pb* peak has a low confidence level (0.59) and its position moves from 714,000 to 96,000 upon removal of the secular component (Tables 2 and 3). Neither 3 He/4He nor eNd variations seem to have a periodic component in that time window. Input of the K/L component does not seem to be a periodically recurring phenomenon. [21] DePaolo [1996] discussed the existence of recurring geochemical patterns in the pilot core with the help of a chromatographic model. A particularly interesting application is the possibility of relating the phase separation between the various isotopic patterns to the hydrodynamic dispersivity of the porous flow. The chromatographic model, however, assumes Henry’s law behavior for all the  melt components, and hence a complete decoupling of trace and major elements. It is therefore required that incompatible elements, such as Nd and Hf, move faster than the major elements which make up the melt itself [Albare`de, 1995, p. $517]. Although the assumption that the mobile phase is inert with respect to the stationary phase is acceptable for column chromatography in the laboratory, it is inadequate for basaltic melts that may react with the porous peridotite matrix. Consequently, we consider that the recurring character of the isotopic patterns in the HSDP-2 core reveal source characteristics rather than hydrodynamic features.  4.5. Size of Mantle Heterogeneities [22] With rare exceptions [Nataf and VanDecar, 1993], plume conduits away from ridge axes are too thin to be observed by seismic waves. It is usually assessed that an upper limit of the conduit diameter is the Fresnel zone of body seismic waves, 20 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Figure 12. Relationship between upwelling in the hot spot conduit, plate velocity, and source sampling by melting and eruptions. Models such as those used by Ribe and Christensen [1999] and DePaolo et al. [2001] assume a timeinvariant geochemical cross section of the conduit (left). A disk represents the surface over which melt is collected by the activity of a single volcano. In such a case, time series reveal the horizontal geochemical gradients across the conduit integrated over the collected area. For high upwelling rates, however, the vertical component of heterogeneities cannot be neglected (right). Although the actual situation may be a combination of the two models, the vertical component of plume heterogeneity is likely to be very significant.  i.e., the interval between two objects that scattering at the P and S wavelengths does not allow to separate. It is therefore assumed that plume conduits must be narrower than this wavelength, i.e., on the order of a few tens of kilometers, a diameter confirmed by geodynamic modeling [Loper and Stacey, 1983; Olson et al., 1993].  centered on the summit. They estimated the instantaneous eruption rate Qe as the integral of the melt supply rate q (x, y) of the plume conduit over the area of a disk of radius r as:  [23] The shape and volume of the mantle source in the plume conduits that Hawaiian volcanoes actually tap is unclear. Consequently, the shape of isotopic heterogeneities becomes relevant to the interpretation of their record in volcanic time series. Ribe and Christensen [1999] assumed that melt is extracted vertically at steady state and considered that, when moving across the plume conduit over their lifetime (approximately 1 My), Hawaiian volcanoes tap a disk-shaped volume  where x is the coordinate on the lithosphere along the volcanic chain axis, y the transverse coordinate, V the plate velocity, and rc the basalt density. They determined that a value of r = 33.6 km accounts for the typical volume of a single volcano. We can likewise determine the instantaneous concentration of a particular element in the erupted lava as:  Qe ðtÞ ¼  Ce ðt Þ ¼  1 rc Qe ðt Þ  1 rc  Z qð x þ Vt; yÞdx dy C ðrÞ  Z qð x þ Vt; yÞC ð x þ Vt; yÞdx dy C ðr Þ  21 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  BASE OF THE LITHOSPHERE  ∆V  ∆V  ∆l  vz 2R Figure 13. Cartoon illustrating the relationship between the hot spot conduit and the volcanic islands. R is the radius of the conduit, vz is the upwelling velocity. Dl is the height of the mantle column and DV is the volume of mantle rising in the conduit during the time corresponding to the recurrence of isotopic patterns.  This complex formula expresses the rather trivial statement that a basalt has the composition of its constitutive melt packets weighted by their respective sizes. The model used by DePaolo et al. [2001] uses a geometry similar to that of Ribe and Christensen [1999] and a steady state melting column to constrain the melt supply rate q(x, y) and then inverts the set of geochemical data on Hawaiian volcanoes to produce C(x, y) maps and their isotopic analogs. Both Ribe and Christensen [1999] and DePaolo et al. [2001] therefore assume that the geochemical cross section of the plume conduit is time invariant. A similar assumption was  used by J. Eisele et al. (submitted manuscript, 2002) to model their time series of Pb isotopes: these authors view heterogeneities as very long filaments stretched parallel to the plume axis such that fluctuating isotopic compositions exclusively reflect the density of filaments across the plume conduit. [24] However, whether volcanic successions record a horizontal compositional gradient across the plume conduit at the successive locations of the volcano summit or vertical heterogeneities in the plume conduit is arguable (Figure 12). If we ignore the upward broadening of the plume con22 of 27  3  Geochemistry Geophysics Geosystems  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  1  0.9  1.002  meas 208Pb/204Pb /NIST 981  0.8  fraction eluted  0.7  0.6  0.5  0.4  1.001  1  0.3 0.999  0.2  0.1  0  0  5  10  15  eluted volume (ml)  0.998  7  8  9  10  11  12  eluted volume (ml)  Figure 14. Elution of NIST 981 Pb standard through an anion-exchange resin AG1-X8 (200-400 mesh) in 2.5 N HCl medium. Left: elution data for 208Pb (open symbols) and fit by a Gaussian curve [Mare´chal and Albare`de, 2002]. Right: deviation of the 208Pb/204Pb ratio in each fraction with respect to the mean value. Solid line: fit using a resin/ solution Pb partition coefficient of 30 and a difference of 40 ppm per mass unit of the partition coefficients between the two isotopes.  duit, the equation for advective transport in two dimensions can be formulated as @C @C @C ¼ ÀV À vz @t @x @z  where vz is the upwelling velocity. For a hot spot the strength of Hawaii, modeling predicts centerline vertical velocities of about 1 m/yr [Loper and Stacey, 1983; Olson et al., 1993; Schubert et al., 2001], which indicates that the vertical component of heterogeneity gradients cannot be neglected with respect to the radial components. Two factors will therefore contribute to increase the relative importance of the second term on the right-hand side of the above equation with respect to the first term:  (i) the plate velocity V is likely to be smaller than the upwelling velocity vz, and (ii) cross-sectional heterogeneities (as measured by @C/@x) will be horizontally averaged by magma extraction. [25] The size and shape distribution of heterogeneities is largely unknown since numerical experiments of continuum mechanics, such as those used by Farnetani and Richards [1995], Ribe and Christensen [1999], and Farnetani et al. [2002] must mesh the mantle using elements with a dimension still fairly coarse (10 – 20 km) with respect to the radius of the conduit and the source disk. Stretching in plume conduits is, however, potentially significant, although far less important than in plume heads [Farnetani et al., 2002]. With 23 of 27  3  Geochemistry Geophysics Geosystems  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  3.605  Pb elution on AG1-X8 in HCl 11  y = 1.44 x - 0.35 r = 0.996  3.604  ln 208Pb/204Pb  10 9 7  3.603  6  8  5 4  3  3.602 2 1  3.601  2.739  2.740  2.741  ln 207Pb/204Pb Figure 15. Mass dependence of the mass fractionation upon elution of Pb on AG1-X8 resin in 2.5 N HCl medium. Numbers refer to the order of eluted fractions. The observed slope (1.44) is only slightly higher than the theoretical slope of ln(208/204)/ln(207/204) = 1.33.  stretching simply reflecting velocity gradients, its effect on bulk compositions is minimized by the mere fact that the very slow outer conduit layers will have a very small contribution to the overall material budget of the plume. The reason why heterogeneities stretched along the conduit walls are unlikely to contribute significantly to the magma budget is the same as that accounting for the apparent lack of entrained upper mantle in the source of hot spot basalts [Farnetani and Richards, 1995; Ribe and Christensen, 1999]: temperature dependence of viscosity in a strong temperature gradient confines plume flow to a well-defined high-velocity conduit with a rather sharp, efficiently ‘‘lubricated’’ edge, and with piston flow as the best analogue. We do not claim here that the  plume-conduit cross section is geochemically homogeneous (‘‘striped’’ plume), and the very existence of contrasting Kea and Loa trends is there to illustrate it. However, ignoring the vertical component of heterogeneities along the conduit is bound to result in inadequate description of the plume. [26] If the observed peaks in our periodograms (Figure 11) are fundamental and not harmonic modes of much slower fluctuations, they may be taken as indicating that the diameter of the mantle plume is at least on the order of a few kilometers to a few tens of kilometers. Davies [1988, 1999] and Sleep [1990] argued that the Hawaiian swell, i.e., the area of shallower seafloor depth around 24 of 27  3  Geochemistry Geophysics Geosystems  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  Deviation from NIST 981 values in permil  0.0  -0.2  -0.4  206  Pb/  204  Pb  -0.6  -0.8 207  Pb/  204  Pb  -1.0  208  -1.2  Pb/  204  Pb  -1.4 0.0  0.2  0.4  0.6  0.8  1.0  Yield  Figure 16. Potential mass fractionation resulting from incomplete yield of Pb extraction chemistry.  the hot spot, is due to upwelling of hot and therefore light material. In contrast, Phipps Morgan et al. [1995] argued that a large part of the buoyancy flux may be due to the depletion of the upwelling material. From the geometry of the Hawaiian swell and the assumption that the excess buoyancy is entirely of thermal origin, the flux of plume material Qp was inferred to range from 7.5 km3/yr [Davies, 1988] to 9.5 km3/yr [Sleep, 1990], which Ribe and Christensen [1999] reduced to less than 4 km3/yr. The magma production of the recent Hawaiian volcanoes is of about 0.1–0.2 km3/yr [Swanson, 1972; Clague and Dalrymple, 1987; Moore, 1987; Watts and ten Brink, 1989] but may drop to a mean value of 0.03 km3/yr or less when integrated over the entire volcanic chain [Clague and Dalrymple, 1987; Davies, 1999]. Such a value corresponds to 1–5 percent of the hot spot upwelling rate, a number consistent with the degrees of melting inferred from petrological and geochemical arguments. [27] For isotopic fluctuations modulated by plume flux, it is possible to evaluate the amount of plume material upwelling between two surges of a partic-  ular component and elaborate on characteristic length scales associated with isotopic heterogeneities. Hafnium isotope and 207Pb/204Pb fluctuations come in ‘‘cycles’’ of T = 50,500 years and therefore, using Ribe and Christensen’s [1999] value of 4 km3/yr, are associated with packets of mantle with a volume DV = QpT (Figure 13) of about 202,000 km3. The range of plume conduit radii and upwelling velocities must be consistent with continuum mechanics constraints [Schubert et al., 2001]. If we assume a maximum radius R of the hot spot conduit of about 50 km (as suggested by the width of the Fresnel zone of P and S waves), the layers with a same mean isotope composition must be separated by an interval Dl = DV/4pR2 of about 6.5 km and the upwelling velocity vz = Dl/T must be 13 cm/yr. Such a low velocity is equivalent to the plate velocity and therefore is at odds with the apparent fixity (or near fixity) of the Hawaiian hot spot over most of the Cenozoic. With a conduit radius R = 10 km, the separation Dl would increase to 160 km and the upwelling velocity vz to 3 m/yr, a value more consistent with the hot spot remaining fixed in the upper mantle circulation. All these results are contingent on the 25 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  assumption that we successfully identified the fundamental mode and not the harmonics: if, as suggested by J. Eisele et al. (submitted manuscript, 2002), the strongest Pb isotope period is 100 kyr rather than 50 kyr, the length scale of heterogeneities should be modified accordingly. Only future deepening of the HSDP-2 drill hole may resolve this ambiguity.  Appendix A: Evidence for Pb Isotope Fractionation on Anion-Exchange Columns in HCl Medium [28] This experiment, carried out in the clean laboratory in Lyon, was designed to elucidate whether incomplete yields of Pb elutions can account for significant isotope fractionation during chemical processing of samples. Twenty micrograms of NIST 981 Pb standard were evaporated, loaded in 2.5 N HCl on 1ml AG1-X8 200-400 mesh anion-exchange resin, and then eluted with the same acid. The 208Pb signal and the Pb isotopic composition of each fraction were measured on the Plasma 54. We found that 206Pb/204Pb, 207Pb/204Pb, and 208Pb/204Pb were fractionated by several permil across the elution peak (Figure 14) and that isotopic fractionation was mass dependent (Figure 15). We used the same procedure, essentially a fit of the elution curve of each isotope by an error function, as in the work of Mare´chal and Albare`de [2002], to retrieve the partition coefficient of each isotope. We found that a partition coefficient of 30, which is consistent with the classical results of Kraus and Nelson [1958], and a variation of the partition coefficients of 40 ppm per mass unit account extremely well for the observed results. We then integrated the isotopic ratios of the individual fractions collected up to a certain yield (Figure 16) in order to evaluate the potential isotopic bias introduced by incomplete recovery of sample Pb. Although the conditions used to establish mass-dependent fractionation of Pb isotopes on the column are slightly different from those used for sample elution, notably a larger column volume and a weaker HCl normality, they illustrate that a chemistry yield of less than 100 percent may introduce a mass bias that could be mistaken for a mass spectrometry artifact.  Acknowledgments [29] We thank Philippe Te´louk for assistance with the Plasma 54 and the French Institut National des Sciences de l’Univers through the program Inte´rieur de la Terre and the Belgian Fonds National de la Recherche Scientifique for financial support. Reviews by Wafa Abouchami, Al Hofmann, Don DePaolo, Fred Frey, and Bill White and discussions with Cinzia Farnetani and Yanick Ricard helped improve the manuscript considerably.  References Abouchami, W., S. J. G. Galer, and A. W. Hofmann, High precision lead isotope systematics of lavas from the Hawaiian Scientific Drilling Project, Chem. Geol., 169, 187 – 209, 2000a. Abouchami, W., A. W. Hofmann, J. Eisele, and S. J. G. Galer, Long-term Pb isotope record of sources within the Hawaiian plume, Eos Trans. AGU Fall Meet. Suppl., 81(48), F1355, 2000b. Albare`de, F., Introduction to Geochemical Modeling, 543 pp., Cambridge Univ. Press, New York, 1995. Anbar, A. D., J. E. Roe, J. Barling, and K. H. Nealson, Nonbiological fractionation of iron isotopes, Science, 288, 126 – 128, 2000. Baker, J., and T. Waight, Pb isotope analysis using Tl and a 207 Pb-204Pb spike on a double focusing MC-ICPMS, in Proceedings of the Goldschmidt Conference, Geochim. Cosmochim. Acta, 66, A44, 2002. Blichert-Toft, J., and F. Albare`de, Hf isotopic compositions of the Hawaii Scientific Drilling Project core and the source mineralogy of Hawaiian basalts, Geophys. Res. Lett., 26, 935 – 938, 1999. Blichert-Toft, J., C. Chauvel, and F. Albare`de, Separation of Hf and Lu for high-precision isotope analysis of rock samples by magnetic sector-multiple collector ICP-MS, Contrib. Mineral. Petrol., 127, 248 – 260, 1997. Blichert-Toft, J., F. A. Frey, and F. Albare`de, Hf isotope evidence for pelagic sediments in the source of Hawaiian basalts, Science, 285, 879 – 882, 1999. Chauvel, C., and J. Blichert-Toft, A hafnium isotope and trace element perspective on melting of the depleted mantle, Earth Planet. Sci. Lett., 190, 137 – 151, 2001. Clague, D. A., and G. B. Dalrymple, The Hawaiian-Emperor volcanic chain, U.S. Geol. Surv. Prof. Pap., 1350, 5 – 54, 1987. Davies, G. F., Ocean bathymetry and mantle convection, 1, Large-scale flow and hotspots, J. Geophys. Res., 93, 10,467 – 10,480, 1988. Davies, G. F., Dynamic Earth. Plates, Plumes and Mantle Convection, 458 pp., Cambridge Univ. Press, New York, 1999. DePaolo, D. J., High-frequency isotopic variations in the Mauna Kea tholeiitic basalt sequence: Melt zone dispersivity and chromatography, J. Geophys. Res., 101, 11,855 – 11,864, 1996. DePaolo, D. J., and E. M. Stolper, Models of Hawaiian volcano growth and plume structure: Implications of results from 26 of 27  Geochemistry Geophysics Geosystems  3  G  blichert-toft et al.: hawaiian hot spot dynamics 10.1029/2002GC000340  the Hawaii Scientific Drilling Project, J. Geophys. Res., 101, 11,643 – 11,654, 1996. DePaolo, D. J., J. G. Bryce, A. Dodson, D. L. Shuster, and B. M. Kennedy, Isotopic evolution of Mauna Loa and the chemical structure of the Hawaiian plume, Geochem. Geophys. Geosyst., 2, 2000GC000139, 2001. Eiler, J. M., K. A. Farley, and E. M. Stolper, Correlated helium and lead isotope variations in Hawaiian lavas, Geochim. Cosmochim. Acta, 62, 1977 – 1984, 1998. Farnetani, C. G., and M. A. Richards, Thermal entrainment and melting in mantle plumes, Earth Planet. Sci. Lett., 136, 251 – 267, 1995. Farnetani, C. G., B. Legras, and P. J. Tackley, Mixing and deformations in mantle plumes, Earth Planet. Sci. Lett., 196, 1 – 15, 2002. Frey, F. A., and J. M. Rhodes, Intershield geochemical differences among Hawaiian volcanoes: Implications for source compositions, melting process and magma ascent paths, Philos. Trans. R. Soc. London, Ser. A, 342, 121 – 136, 1993. Garcia, M. O., B. Jorgenson, J. J. Mahoney, E. Ito, and A. J. Irving, An evaluation of temporal geochemical evolution of Loihi Seamount lavas: Results from Alvin submersible dives, J. Geophys. Res., 98, 357 – 380, 1993. Garcia, M. O., D. J. P. Foss, H. B. West, and J. J. Mahoney, Geochemical and isotopic evolution of Loihi volcano, Hawaii, J. Petrol., 36, 1647 – 1674, 1995. Holcomb, R. T., Eruptive history and long-term behavior of Kilauea volcano, U.S. Geol. Surv. Prof. Pap., 1350, 261 – 350, 1987. Johnson, R. A., and D. W. Winchern, Applied Multivariate Statistical Analysis, 594 pp., Prentice-Hall, Old Tappan. N. J., 1992. Kauahikaua, J., S. Margriter, J. Lockwood, and F. A. Trusdell, Applications of GIS to the estimation of lava flow hazards on Mauna Loa volcano, Hawai’i, in Mauna Loa Revealed: Structure, Composition, History, and Hazards, Geophys. Monogr. Ser., vol. 92, edited by J. M. Rhodes and J. P. Lockwood, pp. 315 – 325, AGU, Washington, D. C., 1995. Kraus, K. A., and F. Nelson, Anion exchange studies, 22, Metal separations by anion exchange, in American Society for Testing and Materials, Spec. Tech. Publ., 195, pp. 27 – 59, 1958. Kurz, M. D., T. C. Kenna, J. C. Lassiter, and D. J. DePaolo, Helium isotopic evolution of Mauna Kea volcano: First results from the 1-km drill core, J. Geophys. Res., 101, 11,781 – 11,791, 1996. Lomb, N. R., Least-squares frequency analysis of unequally spaced data, Astrophys. Space Sci., 39, 447 – 462, 1976. Loper, D. E., and F. D. Stacey, The dynamical and thermal structure of deep mantle plumes, Phys. Earth Planet. Inter., 33, 304 – 317, 1983. Mare´chal, C., and F. Albare`de, Ion-exchange fractionation of copper and zinc isotopes, Geochim. Cosmochim. Acta, 66, 1499 – 1509, 2002.  Mare´chal, C. N., P. Te´louk, and F. Albare`de, Precise analysis of copper and zinc isotopic compositions by plasma-source mass spectrometry, Chem. Geol., 156, 251 – 273, 1999. Moore, J. G., Subsidence of the Hawaiian Ridge, U.S. Geol. Surv. Prof. Pap., 1350, 85 – 100, 1987. Nataf, H.-C., and J. VanDecar, Seismological detection of a mantle plume, Nature, 364, 115 – 120, 1993. Olson, P., G. Schubert, and C. Anderson, Structure of axisymmetric mantle plumes, J Geophys. Res., 98, 6829 – 6844, 1993. Phipps Morgan, J., W. J. Morgan, and E. Price, Hotspot melting generates both hotspot volcanism and hotspot swell, J. Geophys. Res., 100, 8045 – 8062, 1995. Quane, S. L., M. O. Garcia, H. Guillou, and T. P. Hulsebosch, Magmatic history of the East Rift Zone of Kilauea Volcano, Hawaii based on drill core from SOH 1, J. Volcanol. Geotherm. Res., 102, 319 – 338, 2000. Reyment, R., and K. G. Jo¨reskog, Applied Factor Analysis in the Natural Sciences, 371 pp., Cambridge Univ. Press, New York, 1996. Ribe, P. H., and U. R. Christensen, The dynamical origin of Hawaiian volcanism, Earth Planet. Sci. Lett., 171, 517 – 531, 1999. Scargle, J. D., Studies in astronomical time series analysis, 2, Statistical aspects of spectral analysis of unevenly spaced data, Astrophys. J., 263, 835 – 853, 1982. Schubert, G., D. L. Turcotte, and P. Olson, Mantle Convection in the Earth and Planets, Cambridge Univ. Press, New York, 2001. Sims, K. W. W., et al., Chemical and isotopic constraints on the generation and transport of magma beneath the East Pacific Rise, Geochim. Cosmochim. Acta, 66, 3481 – 3504, 2002. Sleep, N. H., Hotspots and mantle plumes: Some phenomenology, J. Geophys. Res., 95, 6715 – 6736, 1990. Staudigel, H., A. Zindler, S. R. Hart, T. Leslie, C.-Y. Chen, and D. Clague, The isotope systematics of a juvenile intraplate volcano: Pb, Nd, and Sr isotope ratios of basalts from Loihi Seamount, Hawaii, Earth Planet. Sci. Lett., 69, 13 – 29, 1984. Swanson, D. A., Magma supply rate at Kilauea volcano, 1952 – 1971, Science, 175, 169 – 170, 1972. Valbracht, P. J., H. Staudigel, M. Honda, I. McDougall, and G. R. Davies, Isotopic tracing of volcanic source regions from Hawaii: Decoupling of gaseous and lithophile magma components, Earth Planet. Sci. Lett., 144, 185 – 198, 1996. Watts, A. B., and U. S. ten Brink, Crustal structure, flexure and subsidence history of the Hawaiian islands, J. Geophys. Res., 94, 10,473 – 10,500, 1989. Weis, D., and F. A. Frey, Isotope geochemistry of Ninetyeast Ridge basalts: Sr, Nd, and Pb evidence for the involvement of the Kerguelen hotspot, in Proceedings of the Ocean Drilling Program, Sci. Results, vol. 121, edited by J. Weissel, J. Peirce, E. Taylor, and J. Alt et al., pp. 591 – 610, Ocean Drilling Progam, College Station, Tex., 1991. White, W. M., F. Albare`de, and P. Te´louk, High-precision analysis of Pb isotopic ratios using multi-collector ICPMS, Chem. Geol., 167, 257 – 270, 2000.  27 of 27  


Citation Scheme:


Citations by CSL (citeproc-js)

Usage Statistics



Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            async >
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:


Related Items