Research
ANDESITE FLOWS AT MT NGAURUHOE, NEW ZEALAND, AND THE IMPLICATIONS FOR POTASSIUM-ARGON "DATING"
ANDREW A. SNELLING, Ph.D. - 1998
Presented at the Fourth International Conference on Creationism
Pittsburgh, PA, August 3-8, 1998
Copyright 1998 by Creation Science Fellowship, Inc.
Pittsburgh, PA USA - All Rights Reserved
KEYWORDS
Andesite, 1949-1975 flows, Mt Ngauruhoe, New Zealand, potassium-argon dating, anomalous model "ages", excess 40Ar*, excess 40Ar* in rocks and minerals, upper mantle, geochemical reservoirs, mantle-crust domains, crustal mixing, magma genesis
ABSTRACT
New Zealand's newest and most active volcano, Mt Ngauruhoe in the Taupo Volcanic Zone, produced andesite flows in 1949 and 1954, and avalanche deposits in 1975. Potassium-argon "dating" of five of these flows and deposits yielded K-Ar model "ages" from <0.27 Ma to 3.5 - 0.2 Ma. "Dates" could not be reproduced, even from splits of the same samples from the same flow, the explanation being variations in excess 40Ar* content. A survey of anomalous K-Ar "dates" indicates they are common, particularly in basalts, xenoliths and xenocrysts such as diamonds that are regarded as coming from the upper mantle. In fact, it is now well established that there are large quantities of excess 40Ar* in the mantle, which in part represent primordial argon not produced by in situ radioactive decay of 40K and not yet outgassed. And there are mantle-crust domains between, and within, which argon circulates during global tectonic processes, magma genesis and mixing of crustal materials. This has significant implications for the validity of K-Ar and 40Ar/39Ar "dating".
INTRODUCTION
Mt Ngauruhoe is an andesite stratovolcano of 2,291 m elevation, rising above the Tongariro volcanic massif within the Tongariro Volcanic Center of the Taupo Volcanic Zone, North Island, New Zealand (Figure 1) [38, 97]. Though not as well publicised as its neighbor, Mt Ruapehu (about 12 km to the south), Ngauruhoe is an imposing, almost perfect cone that rises more than 1,000 m above the surrounding landscape. Eruptions from a central 400 m diameter crater have constructed the steep (33-), outer slopes of the cone [46, 97].
GEOLOGIC SETTING
The Taupo Volcanic Zone, a volcanic arc and marginal basin of the Taupo-Hikurangi arc-trench (subduction) system [13], is a southward extension on the Tonga-Kermadec arc into the continental crustal environment of New Zealand's North Island. It has been interpreted as oblique subduction of the Pacific plate beneath the Australian plate. The zone extends approximately 300 km north-northeast across the North Island from Ohakune to White Island (Figure 1) and is up to 50 km wide in the central part, narrowing northwards and southwards. This volcano-tectonic depression (Taupo-Rotorua depression [42]) comprises four rhyolitic centers (Rotorua, Okataina, Maroa and Taupo), plus the calc-alkaline Tongariro Volcanic Center, part of a young (<0.25 Ma) andesite-dacite volcanic arc with no associated rhyolitic volcanism extending along the eastern side of the zone [38].
The Tongariro Volcanic Center extends for 65 km south-southwest from Lake Taupo at the southern end of the Taupo Volcanic Zone (Figure 1) and consists of four large predominantly andesite volcanoes - Kakaramea, Pihanga, Tongariro and Ruapehu (Figure 2); two smaller eroded centers at Maungakatote and Hauhungatahi; a satellite cone and associated flows at Pukeonake and four craters at Ohakune (Figure 2) [15, 46].
figure 1 - Mt Ngauruhoe in the Taupo Volcanic Zone
Figure 1. The location of Mt Ngauruhoe in the Taupo Volcanic Zone (TVZ), New Zealand, showing the main structural features. The shaded area is the andesite arc, and the inset shows the major components of the boundary between the Australian and Pacific Plates in the New Zealand region (arrows indicate relative motions). Solid triangles are basalt-andesite volcanoes [38].
Most vents lie close to the axis of a large graben in which Quaternary volcanic rocks overlie a basement of Mesozoic greywacke and Tertiary sediments [41, 69]. North-northeast-trending normal faults with throws up to 30 m cut the volcanoes within the graben. Nearly all vents active within the last 10 ka lie on a gentle arc which extends 25 km north-northeast from the Rangataua vent on the southern slopes of Ruapehu through Ruapehu summit and north flank vents, Tama Lakes, Ngauruhoe, Red Crater, Blue Lake and Te Mari craters. None of the young vents lie on the mapped faults, which mostly downthrow towards the axis of the graben. The vent lineation lies above this axis, which is considered to mark a major basement fracture [41, 43, 69] that allows the intrusion of andesite dikes.
The Tongariro volcanics unconformably overlie late Miocene marine siltstones beneath Hauhungatahi, and a minimum age for the onset of volcanism is measured by the influx of andesite pebbles in early Pleistocene conglomerates of the Wanganui Basin to the south [15, 43]. The oldest dated lavas from the Tongariro massif are horn-blende andesites exposed at Tama Lakes between Ngauruhoe and Ruapehu, at 0.26 - 0.003 Ma; from Ruapehu, 0.23 � 0.006 Ma; and from Kakoramea, 0.22 - 0.001 Ma (potassium-argon dates) [91].
Tongariro itself is a large volcanic massif that consists of at least 12 composite cones, the youngest and most active of which is Ngauruhoe. A broad division has been made into older (>20 ka) and younger (<20 ka) lavas [15, 92]. There is a north-northeast alignment of the younger vents of Tongariro, particularly evident between Te Mari and Ngauruhoe.
NGAURUHOE
Ngauruhoe is the newest cone of the Tongariro massif and has been active for at least 2.5 ka [43, 69, 97]. It has been one of the most active volcanoes in New Zealand, with more than 70 eruptive episodes since 1839, when the first steam eruption was recorded by European settlers [41, 69, 97]. Prior to European colonisation the Maoris witnessed many eruptions from the mountain [41]. The first lava eruption seen by European settlers occurred between April and August 1870, with two or three flows witnessed spilling down the north-western flanks of the volcano on July 7 [41, 69]. Following that event there have been pyroclastic (ash) eruptions every few years [69], with major explosive activity in April-May 1948.
figure 2
Figure 2. Location and deposits of the Tongariro Volcanic Center [15, 46].
The next lava extrusion was in February 1949, beginning suddenly with ejection of incandescent blocks, and a series of hot block and ash flows down the north-western slopes on February 9 [41, 69]. The southern sub-crater filled with lava, which by late on February 10 had flowed over the lowest part of the rim and down the north-west slopes of the cone. By February 12 the flow had ceased moving, subsequent mapping placing its volume at about 575,000 m3 (Figure 3) [7, 41]. Further explosive pyroclastic (ash) eruptions followed, reaching a maximum about February 19-21. The eruptions ended on March 3.
The eruption from May 13, 1954 to March 10, 1955 began with explosive ejection of ash and blocks, although red-hot lava had been seen in the crater five months previously [41, 69]. The eruption was remarkable for the estimated large volume of almost 8 million m3 of lava that then flowed from the crater from June through September 1954, and was claimed to be the largest flow of lava observed in New Zealand (that is, by the European settlers) [41, 97]. The lava was actually expelled from the crater in a series of 17 distinct flows on the following dates [40, 41]:
June 4, 30
July 8, 9, 10, 11, 13, 14, 23, 28, 29, 30
August 15 (?), 18
September 16, 18, 26
figure 3 - Mt Ngauruhoe
Figure 3. Map of the north-western slopes of Mt Ngauruhoe showing the lava flows of 1949 and 1954, and the 1975 avalanche deposits [7, 40, 41, 68, 69]. The location of samples collected for this study are marked.
Figure 3 shows the distribution of those 1954 lava flows that are still able to be distinguished on the north-western and western slopes of Ngauruhoe. All flows were of aa lava (as was the February 1949 flow), typified by rough, jagged, clinkery surfaces made up of blocks of congealed lava. The lava flows were relatively viscous, some being observed at close quarters slowly advancing at a rate of about 20 cm per minute [40, 41, 97]. The August 18 flow was more than 18 m thick and still warm almost a year after being erupted. Intermittent explosive eruptions and spectacular lava fountaining during June and July 1954 built a spatter-and-cinder cone around the south sub-crater, modifying the western summit of the mountain. Activity decreased for two months after the last of the lava flows on September 26, but increased again during December 1954 and January 1955 with lava fountaining and many highly explosive pyroclastic (ash) eruptions. The last ash explosion was reported on March 10, 1955, but red-hot lava remained in the crater until June 1955 [40, 41].
After the 1954-1955 eruption, Ngauruhoe steamed semi-continuously, with numerous small eruptions of ash derived from comminuted vent debris. Incandescent ejecta were seen in January 1973, and ash erupted in December 1973 contained juvenile glassy andesite shards [69]. Cannon-like, highly explos-ive eruptions in January and March 1974, the largest since 1954-1955, threw out large quantities of ash and incandescent blocks, one of which was reported as weighing 3000 tonnes and thrown 100 m [69, 97]. Pyroclastic avalanches flowed from the base of large convecting eruption columns, down the west and north slopes of the cone, and the crater became considerably shallower [69, 70].
A series of similar but more violent explosions occurred on February 19, 1975, accompanied by clearly visible atmospheric shock waves and condensation clouds [67, 69, 97]. Ash and blocks up to 30 m across were ejected and scattered within a radius of 3 km from the summit. The series of nine, cannon-like, individual eruptions followed a 1.5 hour period of voluminous gas-streaming emission, which formed a convecting eruption plume between 11 km and 13 km high [68, 69, 97]. The explosions took place at 20-60 minute intervals for more than five hours. Numerous pyroclastic avalanches were also generated by fallback from the continuous eruption column, the avalanches consisting of a turbulent mixture of ash, bombs and larger blocks which rolled swiftly down Ngauruhoe's sides at about 60 km per hour [69, 97]. The deposits from these avalanches and the later explosions accumulated as sheets of debris in the valley at the base of the cone, but did not extend beyond 2 km from the summit. It is estimated that a minimum bulk volume of 3.4 million m3 of pyroclastic material was erupted in the 7-hour eruption sequence on that day [68]. Figure 3 shows the location of these avalanche deposits.
There have been no eruptions since February 1975. A plume of steam or gas is still often seen above the summit of the volcano, as powerful fumaroles in the bottom of the crater discharge hot gases. However, the temperature of these fumaroles in the crater floor has steadily cooled significantly since 1979, suggesting that the main vent is becoming blocked.
SAMPLE COLLECTION
Field work and collection of samples was undertaken in January 1996. The Ngauruhoe area was accessed from State Highway 47 via Mangateopopo Road. From the parking area at the end of the road, the Mangateopopo Valley walking trail was followed to the base of the Ngauruhoe cone, from where the darker-colored recent lava flows were clearly visible and each one easily identified on the north-western slopes against the lighter-colored older portions of the cone (Figure 3).
Eleven 2-3 kg samples were collected - two each from the February 11, 1949, June 4, 1954 and July 14, 1954 lava flows and from the February 19, 1975 avalanche deposits, and three from the June 30, 1954 lava flows. The sample locations are marked on Figure 3. Care was taken to ensure correct identification of each lava flow and that the samples collected were representative of each flow and any variations in textures and phenocrysts in the lavas.
LABORATORY WORK
All samples were sent first for sectioning - one thin section from each sample for petrographic analysis. A set of representative pieces from each sample (approximately 100 g) was then despatched to the AMDEL Laboratory in Adelaide, South Australia, for whole-rock major, minor and trace element analyses. A second representative set (50-100 g from each sample) was sent progressively to Geochron Laboratories in Cambridge (Boston), Massachusetts, for whole-rock potassium-argon (K-Ar) dating - first a split from one sample from each flow, then a split from the second sample from each flow after the first set of results was received, and finally, the split from the third sample from the June 30, 1954 flow.
At the AMDEL Laboratory each sample was crushed and pulverised. Whole-rock analyses were undertaken by total fusion of each powdered sample and then digesting them before ICP-OES for major and minor elements, and ICP-MS for trace and rare earth elements. Fe was analysed for amongst the major elements by ICP-OES as Fe2O3 and reported accordingly, but separate analyses for Fe as FeO were also undertaken via wet chemistry methods. The detection limit for all major element oxides was 0.01%. For minor and trace elements the detection limits varied between 0.5 and 20 ppm, and for rare earth elements between 0.5 and 1 ppm.
The potassium and argon analyses were undertaken at Geochron Laboratories under the direction of Richard Reesman, the K-Ar laboratory manager. No specific location or expected age information was supplied to the laboratory. However, the samples were described as andesites that probably contained "low argon" and therefore could be young, so as to ensure the laboratory took extra care with the analytical work.
Because the sample pieces were submitted as whole rocks, the K-Ar laboratory undertook the crushing and pulverising preparatory work. The concentrations of K2O (weight %) were then measured by the flame photometry method [18], the reported values being the averages of two readings for each sample. The 40K concentrations (ppm) were calculated from the terrestrial isotopic abundance using the measured concentrations of K2O. The concentrations in ppm of 40Ar*, the supposed "radiogenic" 40Ar, were derived using the conventional formula from isotope dilution measurements on a mass spectrometer by correcting for the presence of atmospheric argon whose isotopic composition is known [18]. The reported concentrations of 40Ar* are the averages of two values for each sample. The ratios 40Ar*/Total Ar and 40Ar/36Ar are also derived from measurements on the mass spectrometer and are also the averages of two values for each sample.
PETROGRAPHY AND CHEMISTRY
Clark [11] reported that most of the flows from Ngauruhoe are labradorite-pyroxene andesite with phenocrysts of plagioclase (labradorite), hypersthene and rare augite in a hyalopilitic (needle-like microlites set in a glassy mesostasis) groundmass containing abundant magnetite. However, all lava, lapilli and incandescent blocks that have been analysed from eruptions this century also contain olivine; chemically they may be classed as low-silica (or basaltic) andesites (using the classification scheme of Gill [36]). The published analyses in Table 1 show only trivial changes in composition between 1928 and 1975. In fact, the 1954 and 1974 andesites are so similar that Nairn et al. [70] suggested that a solid plug of 1954-andesite was heated to incandescence and partially remobilised on top of a rising magma column in 1974. This plug was disrupted and blown from the vent as ejecta ranging in texture from solid blocks, through expanded scoria to spatter bombs.
Table 2 lists the whole-rock major element analyses of the eleven samples collected in this study. Comparison of the data for each flow with the corresponding data in Table 1 indicates that in their bulk chemistries all the samples analyzed (and thus all the flows) are virtually identical to one another, the trivial differences being attributable to the statistics of analytical errors, sampling and natural variations. Thus it is not unreasonable to conclude that these basaltic andesites are cogenetic, coming from the same magma and magma chamber, even as they have been observed to flow from the same volcano.
  1 2 3 4 5 6 7 8
SiO2 56.63 57.24 55.90 56.22 56.2 55.83 55.73 56.05
TiO2 0.81 0.81 0.76 0.76 0.8 0.76 0.77 0.84
Al2O3 16.71 16.75 16.90 16.63 16.6 17.03 17.23 16.08
Fe2O3 1.16 1.54 2.10 2.37 1.4 2.14 2.13 3.24
FeO 7.00 6.44 6.30 6.14 7.0 6.36 6.47 5.57
MnO 0.16 0.12 0.15 0.15 0.1 0.16 0.16 0.15
MgO 4.85 4.58 5.20 5.24 5.2 4.79 4.89 5.03
CaO 8.16 7.95 8.40 8.31 8.3 8.37 8.52 7.93
Na2O 2.85 2.74 2.60 3.14 3.1 2.93 2.86 2.82
K2O 1.16 1.49 1.00 1.14 1.2 1.20 1.08 1.23
P2O5 0.16 0.17 0.10 0.17 0.2 0.15 0.15 0.15
H2O 0.31 0.24 0.06 0.19 N/A 0.14 0.10 0.38
TOTAL 99.95 100.07 99.47 100.46 100.1 99.86 100.09 99.47
 
1. Ejecta from March 1928 eruption [11,46]
2. Lava from February 1949 eruption [11, 46, 90]
3. Lava from June 30, 1954 flow [11, 27, 90]
4. 1954 lava (VU 29250) [38]
5. 1954 lava (VU 29250) [15, 69]
6. Average of four lava flows from 1954 eruptions [46, 70]
7. Average of five blocks and bombs from January and March 1974 eruptions [46, 70]
8. Lapilli from February 19, 1975 eruption [46] 		 
 
Table 1. Whole-rock, major-element oxide analyses of recent lava flows at Mt Ngauruhoe, New Zealand, as reported in the literature.
Nevertheless, Nairn et al.[70] suggested that even though the 1949 and 1954 lavas were both olivine-bearing andesite, the chemical analyses (Table 1) showed the 1954 lava to be slightly more basic than the 1949 lava, with slightly higher MgO, CaO and total iron oxides, but lower SiO2 and alkalis. However, these trends are not duplicated with any statistical significance by the analytical results of this study (Table 2). At least they found that their analyses of the 1974 lava blocks and bombs were identical within the limits of error with the 1954 lava (Table 1), which was also substantiated in this study with respect to the 1975 avalanche material and the 1954 lava (Table 2).
Clark [11] and Cole [12] recognized five lava types in the Tongariro Volcanic Center based on the modal proportions of the phenocryst minerals. Graham [37] modified this scheme to six types based on a combination of mineralogy and chemistry, but given their uniform bulk chemistry and petrology, these Ngauruhoe lava flows group together as plagioclase-pyroxene andesite within Graham's -Type 1-. Cole et al. [15] have described -Type 1- lavas as volumetrically dominant within the Tongariro Volcanic Center and as exhibiting coherent chemical trends with increasing silica content. They are relatively Fe-rich and follow a typical calc-alkaline trend on the AFM diagram.
  1A 1B 2A 2B 3A 3B 3C 4A 4B 5A 5B
SiO2 56.7 56.2 55.3 55.8 56.3 55.9 55.6 56.1 55.6 56.0 55.4
TiO2 0.79 0.85 0.74 0.77 0.76 0.75 0.74 0.75 0.84 0.79 0.78
Al2O3 17.2 17.3 16.5 17.3 17.0 16.9 16.7 16.9 17.5 17.0 16.5
Fe2O3* 9.10 9.63 9.26 9.23 9.11 9.17 9.59 9.29 9.61 9.25 9.43
MnO 0.15 0.16 0.15 0.15 0.15 0.15 0.16 0.15 0.16 0.15 0.16
MgO 4.28 3.84 5.21 4.71 4.75 5.00 5.09 4.71 3.84 4.31 5.27
CaO 7.61 7.93 8.22 8.29 7.95 8.16 8.17 8.00 8.17 7.83 8.56
Na2O 3.08 3.19 2.91 3.03 3.06 2.98 2.95 3.02 3.11 3.08 2.86
K2O 1.15 1.01 1.05 1.00 1.10 1.08 1.06 1.12 1.04 1.10 1.09
P2O5 0.13 0.13 0.12 0.12 0.13 0.13 0.13 0.13 0.13 0.14 0.14
L0I 0.42 0.48 0.51 0.37 0.38 0.50 0.53 0.62 0.42 0.70 0.41
TOTAL 100.51 100.72 99.97 100.77 100.69 100.72 100.72 100.79 100.42 100.35 100.60
*Total Fe as Fe2O3 1. February 11, 1949 flow, samples A and B
2. June 4, 1954 flow, samples A and B
3. June 30, 1954 flow, samples A, B and C
4. July 14, 1954 flow, samples A and B
5. February 19, 1975 flow, samples A and B
 
Table 2. Whole-rock, major-element oxide analyses of five recent lava flows at Mt Ngauruhoe, New Zealand (Analyst: AMDEL, Adelaide; April 1996).
Adapting the terminology of Gill [36], the Ngauruhoe lavas are described as basic andesites (53-58 wt% SiO2) [15]. Their designation as plagioclase-pyroxene andesites is based on the predominant phenocrysts present, with plagioclase greater than or equal to pyroxene. Two modal analyses are listed in Table 3 which very closely resemble the samples collected for this study.
Component 1 2
Plagioclase
Augite
Orthopyroxene
Olivine
Iron Oxide
Xenoliths
Groundmass		 22.6
2.6
6.0
0.2
-
2.6
66.0		 21.6
2.6
5.8
0.2
g*
4.5
65.3
TOTAL 100.0 100.0

g* in groundmass
1. Ngauruhoe VU 29250 [15], a 1954 flow.
2. Olivine-bearing low-Si andesite, June 30, 1954 Ngauruhoe flow [12].
Table 3. Modal analyses of two recent lava flows at Mt Ngauruhoe, New Zealand, as reported in the literature.
All samples of the five lava flows examined in this study exhibited a porphyritic texture, with phenocrysts (up to 3 mm across) consistently amounting to 35-40% by volume. The phenocryst assemblage is dominated (2:1) by plagioclase, but orthopyroxene and augite (clinopyroxene) are always major components, while olivine and magnetite are only present in trace amounts. This POAM phenocryst assemblage is a typical anhydrous mineralogy [15]. The groundmass consists of microlites of plagioclase, orthopyroxene and clinopyroxene, and is crowded with minute granules of magnetite and/or Fe-Ti oxides. Small amounts (9-10%) of brown transparent (acid-residuum) glass are also present, and the overall texture is generally pilotaxitic.
Steiner [90] stressed that xenoliths are a common constituent of the 1954 Ngauruhoe lava, but also noted that Battey [7] reported the 1949 Ngauruhoe lava was rich in xenoliths. All samples in this study contained xenoliths, including those from the 1975 avalanche material. However, many of these aggregates are more accurately described as glomerocrysts and mafic (gabbro, websterite) nodules [39]. They are 3-5 mm across, generally have hypidiomorphic-granular textures, and consist of plagioclase, orthopyroxene and clinopyroxene in varying proportions, and very occasionally olivine. The true xenoliths are often rounded and invariably consist of fine quartzose material. Steiner [90] also described much larger xenoliths of quartzo-feldspathic composition and relic gneissic structure.
The plagioclase phenocrysts have been reported as ranging in composition from An89 to An40 (andesine to bytownite), but in Ngauruhoe lavas are usually labradorite (An68-55). They are subhedral and commonly exhibit complex oscillatory zoning with an overall trend from calcic cores to sodic rims [12, 14, 15]. Thin outer rims are usually compositionally similar to groundmass microlites. Twinning and hourglass structures are common.
Orthopyroxene predominates (>2:1) over clinopyroxene. Subhedral-euhedral orthopyroxene is typically pleochroic and sometimes zoned. Compositions range from Ca4 Mg74 Fe22 to Ca3 Mg47 Fe50 [14, 15], but representative bulk and partial analyses of Ngauruhoe orthopyroxenes [12, 26, 38] indicate a hypersthene composition predominates, which is confirmed by optical determinations [11, 12]. Euhdral-subhedral clinopyroxene is typically twinned and zoned, but compositions show a restricted range of Ca43 Mg47 Fe10 to about Ca35 Mg40 Fe25, all of which is augite [15, 39].
The olivine present is strongly magnesian, analyses indicating some compositional zoning from Fo88 to Fo78. The magnetite present in the groundmass is titanomagnetite, judging from the amount of TiO2 present in whole-rock analyses (Tables 1 and 2), but some ilmenite is likely to occur sporadically in association with it [15, 39].
K-Ar RESULTS
All analytical results received from Geochron Laboratories are listed in Table 4, grouped in chronological order according to the historic date of each flow. The 40Ar* quantity refers to the amount of radiogenic 40Ar measured in each sample. All other quantities are self-explanatory, some of them being calculated from the analytical results supplied by the laboratory.
The "age" of each sample is calculated from the analytical results using the general model-age equation [21, 28]:-
formula 1(1)
where:
t=the "age"
l =the decay constant of the parent isotope
Dt=the number of daughter atoms in the rock presently
Do=the number of daughter atoms initially in the rock
Pt=the number of parent atoms presently in the rock
To date a rock, Dt and Pt are measured, and equation (1) can then be used if an assumption about the original quantity of daughter atoms (Do) is made. Applied specifically to K-Ar dating, equation (1) thus
becomes:-
formula 2 (2)
where:
t=the "age" in Ma (millions of years)
5.543 x 10-10=the current estimate for the decay constant of 40K
0.1048=the estimated fraction of 40K decays producing 40Ar
40Ar*/40K=the calculated mole ratio of radiogenic 40Ar to 40K in the sample
It should be noted that to make equation (2) equivalent to equation (1), 40Ar* is assumed to be equal to (Dt - Do), which thus means the 40Ar* measurement has included within it an assumption concerning the initial quantity of 40Ar in the rock, namely, no radiogenic argon is supposed to have existed when the rock formed (that is, Do=0). Thus equation (2) yields a "model age" assuming zero radiogenic argon in the rock when it formed.
The model ages listed in Table 4 range from <0.27 Ma to 3.5 - 0.2 Ma. However, it should be noted that the samples, one from each flow, that yielded model ages of <0.27 Ma and <0.29 Ma (that is, below the detection limits of the equipment for 40Ar*) were all processed at the K-Ar laboratory in the same batch, suggesting the possibility of a systematic problem with the analytical procedure and equipment (in particular, the gas extraction "line"). When this question was raised with the laboratory manager, Richard Reesman, he kindly rechecked his equipment and then re-ran several of the samples, producing similar results and thus ruling out a systematic laboratory "error".
Click here to view Table 4
Table 4. K-Ar analytical results and model ages for five recent lava flows at Mt Ngauruhoe, New Zealand (Analyst: Geochron Laboratories, Boston; July 1996, December 1996, September 1997 and January 1998). Constants used: 40K/K=1.193 x 10-4g/g; Fraction of 40K decays to 40 Ar*=0.1048; Decay constant of 40K=5.543 x 10-10yr-1; Atmospheric 40Ar/36Ar=295.5
However, an independent blind check was then made, by submitting to the K-Ar laboratory duplicate splits from two samples already analysed, to establish if results really were reproducible. The samples chosen were the A and B samples of the June 30, 1954 flow, because their first splits had produced the lowest and highest model ages, <0.27 Ma and 3.5 - 0.2 Ma respectively. The results of these additional analyses are shown in Table 4 as A#2 and B#2, and yielded model ages of 1.3 - 0.3 Ma and 0.8 - 0.2 Ma respectively. Clearly, reproducibility was not obtained, but this is not surprising given the analytical uncertainties at such low to negligible levels of 40Ar*, which are at the detection limits of the laboratory's equipment [84].
DISCUSSION
In spite of the wide variations in model "ages" obtained between and within these recent lava flows, and of the difficulties obtaining analytical reproducibility, it is apparent that the cause of the anomalous K-Ar model "ages" is excess argon in the lavas, that is, non-zero concentrations of radiogenic argon (40Ar*). This of course is contrary to the assumption of zero radiogenic argon in equation (2) for calculating the model "ages". When analysed the oldest of the lavas was less than 50 years old, so there has been insufficient time since cooling for measurable quantities of 40Ar* to have accumulated within the lavas due to the slow radioactive decay of 40K. Thus the measurable 40Ar* can't be from in situ radioactive decay since cooling, and therefore must have been present in the molten lavas when extruded from Mt Ngauruhoe.
No Radiogenic Argon Assumption Violated by Many Anomalous "Ages"
The assumption of no radiogenic argon (40Ar*) when the rocks formed is usually stated as self-evident. For example, Geyh and Schleicher [35, p.56] state:
What is special about the K-Ar method is that the daughter nuclide is a noble gas, which is not normally incorporated into minerals and is not bound in the mineral in which it is found.
Similarly, Dalrymple and Lanphere [18, p.46] state:
"a silicate melt will not usually retain the 40Ar that is produced, and thus the potassium-argon clock is not "set" until the mineral solidifies and cools sufficiently to allow the 40Ar to accumulate in the mineral lattice."
Dalrymple [17, p.91] has recently put the argument more strongly:
"The K-Ar method is the only decay scheme that can be used with little or no concern for the initial presence of the daughter isotope. This is because 40Ar is an inert gas that does not combine chemically with any other element and so escapes easily from rocks when they are heated. Thus, while a rock is molten the 40Ar formed by decay of 40K escapes from the liquid."
However, these dogmatic statements by Dalrymple are inconsistent with even his own work on historic lava flows [16], some of which he found had non-zero concentrations of 40Ar* in violation of this key assumption of the K-Ar dating method. He does go on to admit that "Some cases of initial 40Ar remaining in rocks have been documented but they are uncommon" [17], but then refers to his study of 26 historic, subaerial lava flows [16]. Five (almost 20%) of those flows contained 'excess argon', but Dalrymple still then says "that 'excess' argon is rare in these rocks"! The flows and their "ages" were [16]:  "
Hualalai basalt, Hawaii (AD1800-1801) 1.6 - 0.16 Ma
1.41 -0.08 Ma
Mt Etna basalt, Sicily (122BC) 0.25 -0.08 Ma
Mt Etna basalt, Sicily (AD1792) 0.35 - 0.14 Ma
Mt Lassen plagioclase, California (AD1915) 0.11 - 0.03 Ma
Sunset Crater basalt, Arizona (AD1064�1065) 0.27 - 0.09 Ma
0.25 - 0.15 Ma
Far from being rare, there are numerous examples reported in the literature of excess 40Ar* in recent or young volcanic rocks producing excessively old whole-rock K-Ar "ages":
Akka Water Fall flow, Hawaii (Pleistocene) 32.3 - 7.2 Ma [58]
Kilauea Iki basalt, Hawaii (AD1959) 8.5 - 6.8 Ma [58]
Mt Stromboli, Italy, volcanic bomb (Sept. 23, 1963) 2.4 - 2 Ma [58]
Mt Etna basalt, Sicily (May 1964) 0.7 - 0.01 Ma [58]
Medicine Lake Highlands obsidian, Glass Mountains, California (<500 years old) 12.6 - 4.5 Ma [58]
Hualalai basalt, Hawaii (AD1800�1801) 22.8 - 16.5 Ma [58]
Rangitoto basalt, Auckland, NZ (<800 yrs old) 0.15 - 0.47 Ma [65]
Alkali basalt plug, Benue, Nigeria (<30 Ma) 95 Ma [30]
Olivine basalt, Nathan Hills, Victoria Land, Antarctica (<0.3 Ma) 18.0 - 0.7 Ma [1]
Anorthoclase in volcanic bomb, Mt Erebus, Antarctica (1984) 0.64 - 0.03 Ma [25]
Kilauea basalt, Hawaii (<200 yrs old) 21 - 8 Ma [72]
Kilauea basalt, Hawaii (<1,000 yrs old) 42.9 - 4.2 Ma [19]
30.3 - 3.3 Ma [19]
East Pacific Rise basalt (<1 Ma) 690 - 7 Ma [34]
Seamount basalt, near East Pacific Rise (<2.5 Ma) 580 - 10 Ma [33]
700 - 150 Ma [31]
East Pacific Rise basalt (<0.6 Ma) 24.2 - 1.0 Ma [23]
Other studies have also reported measurements of excess 40Ar* in lavas. Fisher [29] investigated submarine basalt from a Pacific seamount and found "the largest amounts of excess 4He and 40Ar ever recorded" (at that time). McDougall [64] not only found "extraneous radiogenic argon present in three of the groups of basalt flows" on the young volcanic island of R�union in the Indian Ocean, but "extraneous argon" was also "detected in alkali feldspar and amphibole in hyperbyssal drusy syenites that are exposed in the eroded core of Piton des Neiges volcano." Significant quantities of excess 40Ar* have also been recorded in submarine basalts, basaltic glasses and olivine phenocrysts from the currently active Hawaiian volcanoes, Loihi Seamount and Kilauea, as well as on the flanks of Mauna Loa and Hualalai volcanoes, also part of the main island of Hawaii [51, 95], and in samples from the Mid-Atlantic Ridge, East Pacific Rise, Red Sea, Galapagos Islands, McDonald Seamount and Manus Basin [62, 89]. Patterson et al. [76] claimed that some of the initial Loihi analytical results were due to atmospheric contamination of the magma either during intrusion or eruption, but subsequent work [51, 95] has confirmed that the excess 40Ar* is not from atmospheric contamination at all.
 
Excess 40Ar* Occluded in Minerals
Austin [2] has investigated the 1986 dacite lava flow from the post-October 26, 1980 lava dome within the Mt St Helens crater, and has established that the 10-year-old dacite yields a whole-rock K-Ar model "age" of 0.35 - 0.05 Ma due to excess 40Ar* in the rock. He then produced concentrates of the constituent minerals, which yielded anomalous K-Ar model "ages" of 0.34 - 0.06 Ma (plagioclase), 0.9 - 0.2 Ma (hornblende), 1.7 - 0.3 Ma (pyroxene) and 2.8 - 0.6 Ma (pyroxene ultra-concentrate). While these mineral concentrates were not ultra-pure, given the fine-grained glass in the groundmass and some Fe-Ti oxides, it is nonetheless evident that the excess 40Ar* responsible for the anomalous K-Ar "ages" is retained within the different constituent minerals in different amounts. Furthermore, the whole-rock "age" is very similar to the "age" of the plagioclase concentrate because plagioclase is the dominant constituent of the dacite.
That the excess 40Ar* can be occluded in the minerals within lava flows, rather than between the mineral grains, has been established by others also. Laughlin et al. [60] found that the olivine, pyroxene and plagioclase in Quaternary basalts of the Zuni-Bandera volcanic field of New Mexico contained very significant quantities of excess 40Ar*, as did the olivine and clinopyroxene phenocrysts in Quaternary flows from New Zealand volcanoes [78]. Similarly, Potts et al. [83] separated olivine and clinopyroxene phenocrysts from young basalts from New Mexico and Nevada and then measured "ubiquitous excess argon" in them. Damon et al. [20] have reported several instances of phenocrysts with K/Ar "ages" 1�7 million years greater than that of the whole rocks, and one K/Ar "date" on olivine phenocrysts of greater than 110 Ma in a recent (<13,000 year old) basalt. Damon et al. thus suggested that large phenocrysts in volcanic rocks contain the excess 40Ar* because their size prevents them from completely degassing before the flows cool, but Dalrymple [16] concluded that there does not appear to be any correlation of excess 40Ar* with large phenocrysts or with any other petrological or petrographic parameter.
Most investigators have come to the obvious conclusion that the excess 40Ar* had to have been present in the molten lavas when extruded, which then did not completely degas as they cooled, the excess 40Ar* becoming 'trapped' in the constituent minerals, and in some instances, the rock fabrics themselves. Laboratory experiments have tested the solubility of argon in synthetic basalt melts and their constituent minerals near 1300-C at one atmosphere pressure in a gas stream containing argon [8,9]. When quenched, synthetic olivine in the resultant material was found to contain 0.34 ppm 40Ar*. Broadhurst et al. [8] commented that "The solubility of Ar in the minerals is surprisingly high", and concluded that the argon is held primarily in lattice vacancy defects within the minerals.
In a different experiment, Karpinskaya et al. [54] heated muscovite to 740--860-C under high argon pressures (2,800-5,000 atmospheres) for periods of 3 to 10.5 hours. The muscovite absorbed significant quantities of argon, producing K/Ar "ages" of up to 5 billion years, and the absorbed argon appeared like ordinary radiogenic argon (40Ar*). Karpinskaya subsequently [53] synthesized muscovite from a colloidal gel under similar argon pressures and temperatures, the resultant muscovite retaining up to 0.5 wt% argon at 640-C and a vapor pressure of 4,000 atmospheres. This is approximately 2,500 times as much argon as is found in natural muscovite. These experiments show that under certain conditions argon can be incorporated into minerals and rocks that are supposed to exclude argon when they crystallize.
 
Applications to the Mt Ngauruhoe Andesite Flows
Therefore, the analytical results from the very recent (1949-1975) andesite flows at Mt Ngauruhoe, New Zealand, that yield anomalous K-Ar model "ages" because of excess 40Ar*, are neither unique nor an artifact of poor analytical equipment or technique. This realization that the presence of the excess 40Ar* in these rocks is both real and measurable, and has not been derived from radioactive decay of 40K in situ, leads to the obvious questions as to whether there is any pattern in the occurrences of excess 40Ar*, and from whence came this excess 40Ar*?
It is clear that the excess 40Ar* was in the lavas when they flowed from the Mt Ngauruhoe volcano and were trapped in the andesite as it cooled. That there were gases in the lavas is readily evident from the copious "frozen" bubble holes now in the rock, implying that much of the gas content escaped as the lavas flowed and cooled. When choosing samples, care was taken to select pieces from each flow that were different from one another (e.g., copious "frozen" gas bubble holes compared with virtually no such holes). It is hardly surprising, therefore, that the 40Ar* measurements on four of the five flows were consistent with such differences - the samples from each flow which had very few or virtually no "frozen" gas bubble holes yielded excess 40Ar* and thus anomalous K-Ar model "ages", whereas the other samples from each of these flows which contained copious "frozen" gas bubble holes failed to yield detectable 40Ar* (<0.27 Ma and <0.29 Ma in Table 4).
The exception was the June 30, 1954 flow - not only was this expected relationship between excess 40Ar* and lack of "frozen" gas bubble holes not duplicated, but analyses on duplicate splits off the same samples yielded widely divergent results (<0.27 Ma versus 1.3 - 0.3 Ma and 3.5 - 0.2 Ma versus 0.8 - 0.2 Ma, see Table 4). Thus the presence (or absence) of excess 40Ar* must also depend on which portion of a rock sample is being analysed, which in turn implies dependence on the mineral constituents present, including the glass in the groundmass. As already noted, Austin [2] found widely different amounts of excess 40Ar* in the mineral separates concentrated from Mt St Helens 1986 dacite, while numerous other studies [60, 78, 83, 95] have located excess 40Ar* in phenocrysts.
 
Cooling Rates, Pressures and Potassium Alteration
Another factor is the rate of cooling of lavas. Dalrymple and Moore [19] found that the 1 cm thick glassy rim of a pillow in a Kilauea submarine basalt had greater than 40 times more excess 40Ar* than the basalt interior just 10 cm below. The glassy pillow rim is, of course, produced by rapid quenching of the hot basalt lava immediately it contacts the cold ocean water, so the excess 40Ar* in the lava is rapidly trapped and retained. Dymond [23] obtained similar results on four deep-sea basalt pillows from near the axis of the East Pacific Rise. Dalrymple and Moore [19] also found that the excess 40Ar* contents of the glassy rims of basalt pillows increased systematically with water depth, leading them to conclude that the amount of excess 40Ar* is a direct function of both the hydrostatic pressure and the rate of cooling. In a parallel study, Noble and Naughton [72] reported K-Ar "ages" from zero to 22 Ma with increasing sample depth for submarine basalts probably less than 200 years old, also from the active Kilauea volcano.
Seidemann [86] has reported yet another intriguing relationship. He analysed deep-sea basalt samples obtained from DSDP drillholes in the floor of the Pacific Ocean basin and found K-Ar "ages" increased with increasing K contents of the basalts, a relationship he noted also appeared in similar data published by DSDP staff [86, Figure 1]. In basalt pillows the K content increases from the margin to a maximum at an intermediate distance into the pillows, whereas holocrystalline basalts show a decrease of K inward from the margin [49]. Seidemann [86] concluded, as had others before him, that submarine weathering adds K to the basalts, as does alteration at the time of formation, whereas the glassy pillow margins are largely impervious to seawater. The net result, however, is unreliable K-Ar "dates", because the measured 40Ar* was probably not derived by radioactive decay of the measured 40K contents. Seidemann also determined that sediment cover is not a significant barrier to the diffusion of K into basalt.
It is possible that some of these factors are relevant to the pattern of excess 40Ar* measured in the samples from the recent Mt Ngauruhoe andesite flows. For example, the surfaces of the flows would have cooled more rapidly to crusts on top of the still molten flow interiors, which would certainly have been the case with the August 18, 1954 flow that was reported as being 18m thick. Furthermore, the overburden pressure within the deep interiors of such thick flows would likewise inhibit degassing of the lava as it cooled. However, so long as the crusts on the tops of the flows remained unbroken and intact they would have sealed in the molten lava and its contained gases, including excess 40Ar*. But this sealing was probably short-lived, because today the flows mostly outcrop as piles of pieces of andesite that look like rubble (typical aa lavas). The continued flow of the lavas down the sides of the volcano would have broken up the crusts as soon as they congealed, as would contraction with cooling, thus enabling the molten interiors to degas as they cooled. So the speed of cooling was likely the most relevant factor, and this would have varied laterally and vertically within the flows, even at localised scales of a few centimeters.
Any effects of weathering on the K contents of these flows can be discounted. On the one hand these are subaerial flows that do not appear to have been subjected to leaching or addition of K, or to any K-rich alteration for that matter, while on the other they have very uniform K contents (see Tables 1, 2 and 4). This is not unexpected, given the fact they flowed from the same magma source/chamber close together timewise. The K-Ar data in Table 4 do not reveal any discernible relationship between K contents and K-Ar model "ages" of these flows, unlike the negative correlation found in the Middle Proterozoic Cardenas Basalt upper member flows in Grand Canyon, Arizona [3].
Perhaps the key issues, though, are where this excess 40Ar* has come from, and whether it has been derived from radioactive decay of 40K. One possibility is that the excess 40Ar* can be accounted for by radioactive decay during long term residence of magmas in chambers before eruption. Esser et al. [25] discounted this option for the Mt Erebus anorthoclase phenocrysts. Dalrymple [16] found that whereas the Mt Lassen (1915) plagioclase phenocrysts yielded excess 40Ar* and an anomalous K-Ar model "age", a plagioclase from the 1964 eruption of Surtsey only had argon whose isotopic composition matched that of air. Because phenocrysts usually crystallize from lavas after eruption, they may arbitrarily trap excess 40Ar* during lava cooling, 40Ar* that will thus not be from in situ 40K radioactive decay.
 
Negative K-Ar Model "Ages" and Atmospheric Argon
Another relevant consideration bearing on these issues is the observation noted by Dalrymple [16] that some modern lava samples actually yield negative K-Ar model "ages", apparently due to excess 36Ar. Air has an 40Ar/36Ar ratio of 295.5, but some of Dalrymple's samples had ratios less than 295.5 (and hence negative "ages"). Some of the Mt Ngauruhoe samples in this study also yielded 40Ar/36Ar ratios less than 295.5 (see Table 4). According to the straightforward interpretation of the K-Ar dating methodology, this should be impossible. Dalrymple was not willing to attribute these anomalous ratios to experimental error, and neither was Richard Reesman of Geochron Laboratories.

Dalrymple [16] suggested three possible explanations that might account for the excess 36Ar: (1) incorporation of "primitive argon", (2) production of 36Ar by the radioactive decay of 36Cl, or (3) fractionation of atmospheric argon by diffusion. He rejected the possibility of significant 36Ar formation in situ from nuclear reactions [option (2)] because the Cl content of basalts and the production rate of 36Cl by cosmic-ray neutrons both are too low to account for any significant amount of 36Ar. Instead, Dalrymple seemed to favor option (3), that when atmospheric argon diffused back into lavas as they cooled, 36Ar diffused in preferentially. However, he also recognized the weakness of this argument -- it is difficult to explain why some lavas are enriched in 36Ar while others are not (as at Mt Ngauruhoe also). To be consistent, if fractionation of atmospheric argon occurred during diffusion, then this would mean that even supposedly "zero age" lavas actually have an apparent age, and that most lavas do not degas upon eruption. In fact, depending on how strong the fractionation of 36Ar was during diffusion, it could even be that all lavas do not completely degas.

This only leaves Dalrymple's option (1), that the lavas with the anomalously high 36Ar come from areas of the mantle (and perhaps also the crust) which have primordial argon that has not been diluted with radiogenic 40Ar, and have not completely degassed. However, this means that there is no reason to assume that lavas whose argon matches that in the atmosphere have degassed either, because they may have simply started with argon which matches atmospheric argon. Nevertheless, Dalrymple is convinced that "much of the volatile juvenile content may still be present in volcanic rocks quenched on the ocean floor" [19]. Indeed, Dalrymple has specifically defined excess 40Ar* as 40Ar that is not attributed to atmospheric argon or in situ radioactive decay of 40K [18]. Krummenacher [58] is more cautious, attributing anomalous 40Ar/36Ar ratios and excess 40Ar* to the "mass fractionation effect on argon of atmospheric isotopic composition" trapped in the lavas, as well as to the presence of "magmatic" argon different in isotopic composition.

The Role of Xenoliths

Is the excess 40Ar* simply "magmatic" argon, that is, argon that collects in the magma and then is inherited by the lavas from it? Funkhouser and Naughton [32] found that the excess 40Ar* in the 1800?1801 Hualalai flow, Hawaii, resided in fluid and gaseous inclusions in olivine, plagioclase and pyroxene in ultramafic xenoliths in the basalt. The quantities of excess 40Ar* were sufficient to yield K-Ar model "ages" from 2.6 Ma to 2960 Ma. However, Dalrymple [16] subsequently only used the presence of the ultramafic xenoliths and their excess 40Ar* contained in inclusions to explain partly the excess 40Ar* and anomalous K-Ar model "ages" he obtained from the same 1800?1801 Hualalai flow, suggesting instead that the large single inclusions are not directly responsible for the excess argon in the flows and that the 40Ar* is distributed more uniformly throughout the rocks. Nevertheless, those K-Ar and Ar-Ar geochronologists who are concerned about the excess 40Ar* in their samples undermining their "dating" are careful to check for xenoliths, and xenocrysts. Esser et al. [25] did so and discounted xenocrystic contamination.

Xenoliths are present in the Ngauruhoe andesite flows (Table 3), but they are minor and less significant as the location of the excess 40Ar* residing in these flows than the plagioclase and pyroxene phenocrysts, and the much larger glomerocrysts of plagioclase, pyroxene, or plagioclase and pyroxene that predominate. The latter are probably the early-formed phenocrysts that accumulated together in the magma within its chamber prior to eruption of the lava flows. Nevertheless, any excess 40Ar* they might contain had to have been supplied to the magma from its source. The xenoliths that are in the andesite flows have been described by Steiner [90] as gneissic, and are therefore of crustal origin, presumably from the basement rocks through which the magma passed on its way to eruption.

Noble Gases from the Mantle

With the advent of the necessary technology, the isotopic concentrations and ratios of noble gases (including argon) in rock and mineral samples are now obtainable. Honda et al. [51] have reported such analyses on submarine basalt glass samples from Loihi Seamount and Kilauea, Hawaii, and concluded that helium and neon isotopic ratios in particular, being uniquely different from atmospheric isotopic ratios, are indicative of the mantle source area of the plume responsible for the Hawaiian volcanism rather than from atmopsheric contamination of the magma [76]. The 40Ar/36Ar ratios are consistent with excess 40Ar* having also come with the magma from the mantle. A subsequent study [95], in which a larger suite of basalt glass samples, and also samples of olivine phenocrysts, from the same and additional Hawaiian volcanoes were analysed, concluded that the isotopic systematics indicate that the helium and neon have been derived from the mantle and have not been preferentially affected by secondary processes. Consequently, the excess 40Ar* also in these samples would have been also carried from the upper mantle source area of these basalts by the magma plume responsible for the volcanism. Moreira et al. [66] have suggested, based on new experimental data from single vesicles in mid-ocean ridge basalt samples dredged from the North Atlantic, that the excess 40Ar* in the upper mantle may be almost double previous estimates [89] (that is, almost 150 times more than the atmospheric content relative to 36Ar), and represents a primordial mantle component not yet outgassed. Burnard et al. [10] obtained similar results on the same samples, but maintained that because some of the 36Ar is probably surface-adsorbed atmospheric argon, the upper mantle content of excess 40Ar* could be even ten times higher.

Similar results [94] have been obtained from ultramafic mantle xenoliths in basalts from the Kerguelen Archipelago in the southern Indian Ocean, and the considerable excess 40Ar* measured concluded to be a part of the mantle source signature of this hotspot volcanism. However, it has not only been the suboceanic mantle that has thus been sampled for its excess 40Ar* via such magma plumes. Matsumoto et al. [63] have reported high 40Ar/36Ar ratios in spinel-lherzolites from five eruption centers in the youthful (<7 Ma) Newer Volcanics of southeastern Australia. These anhydrous lherzolites have compositions representative of the upper lithospheric mantle, and the significant excess 40Ar* in them indicates the presence of a subcontinental mantle reservoir with a very high 40Ar/36Ar ratio, and thus substantial excess 40Ar*, similar to that found in mid-ocean ridge and plume/hotspot basalts. Another example is the Cardenas Basalt and associated diabase (Middle Proterozoic) of eastern Grand Canyon, regarded as part of the pervasive mafic mid-continental magmatism of the southwestern United States and thus also sourced from the subcontinental mantle. Austin and Snelling [3] have found that the 40Ar/36Ar?40K/36Ar isochrons for 14 and six samples of these rocks respectively yield initial 40Ar/36Ar ratios of 787 ? 118 and 453 ? 42, indicative of some initial excess 40Ar*.

Sampling the Mantle with Diamonds and Their Inclusions

Another means of "sampling" the mantle is the study of diamonds and their micro-inclusions. It is now firmly established that diamonds are thermodynamically stable in the pressure-temperature regime in the mantle at depths greater than 150 km, and their origin is believed to extend back to the Archean and the early crust of the Earth [55, 56]. Diamonds are formed in a number of processes associated with two rock types, ecologite and peridotite, xenoliths of which are also brought up into the upper crust with diamonds from the upper mantle below continental Precambrian shields (cratons) by kimberlite and lamproite "pipe" eruptions [44, 45, 56]. Even though the host kimberlite or lamproite may be relatively young (even in conventional terms), many diamonds date back to the Archean and thus the early history of the Earth [56, 85]. To account for all this evidence, it is postulated that the formation of most diamonds was closely associated with subduction of the Archean oceanic crust into the mantle [55, 56], the required carbon, which was originally thought to be primordial carbon already in the mantle, now believed to derive from sedimentary marine carbonates and biogenic carbon from bacteria/algae in the sediments subducted with the oceanic crust [24, 56, 57].

The noble gas contents of diamonds are consistent with their ancient and mantle origin, high helium isotopic ratios (290 times the atmospheric ratio) being regarded as primordial and rivalling those measured for the Sun today [44, 73]. Of significance here is the postulation that He, Ar, K, Pb, Th and U are added to the convecting upper mantle circulation, and the proportions and isotopic compositions are strongly determined by entrainment from the lower mantle (below 670 km) [45, 50]. This is reflected in those Ar isotopic measurements that have been made on diamonds and their micro-inclusions.

Rather than focus on attempting to date only diamond micro-inclusions as others had done, Zashu et al. [98] carefully selected 10 Zaire diamonds and examined them for purity before undertaking K-Ar dating analyses of the diamonds themselves. However, at the outset they noted that there had been almost no direct radiometric dating of diamonds except for conventional K-Ar dating, and the results had been questionable due to the possible presence of excess 40Ar*. To avoid this problem, they used the K-Ar isochron dating method. Their experimental data showed good linear correlations, but these isochrons yielded an age of 6.0 ? 0.3 Ga, which of course was unacceptable because these diamonds would be older than the Earth itself. Mistakes in the experimental procedure were easily discounted, so they were forced to conclude that excess 40Ar* was responsible, and that it needed to be in a fluid state to ensure the homogenization necessary to give such a constant 40Ar/K ratio. Alternately, they speculated that the diamonds might differ in K isotopic composition from common potassium, but this was discounted in a follow-up study [81] in which it was found that 40K was present in these diamonds in normal abundance. Because 40Ar/39Ar analyses yielded the same unacceptable "age", it was concluded that the excess 40Ar* was not generated in situ, but was an inherited or "trapped" component from the mantle reservoir when and where the diamonds formed.

These Zaire diamonds are not the only ones which have yielded excess 40Ar*. Phillips et al. [79] used a laser-probe to 40Ar/39Ar date eclogitic clinopyroxene inclusions in diamonds from the Premier kimberlite, South Africa, and found moderate 40Ar/36Ar ratios indicative of much less excess 40Ar* than in the Zaire diamonds. The "age" of these eclogitic diamonds was thus determined to be 1.198 ? 0.014 Ga, much younger than the 3.3 Ga peridotitic diamonds at Kimberley and Finsch [56], also in South Africa, so Phillips et al. [79] interpreted the moderate excess 40Ar* as characteristic of mantle conditions prevailing at the time and in the region of Premier eclogitic diamond formation.

Zashu et al. [98] postulated that the excess 40Ar* in the Zaire diamonds needed to be in a fluid state. Though Navon et al. [71] did not analyse for argon when they investigated fluids in micro-inclusions in diamonds from Zaire and Botswana, they found a high content of volatiles and incompatible elements in the uniform average composition of the micro-inclusions, with the amounts of water and CO2 (in carbonates) almost an order of magnitude higher than the volatile contents of kimberlites and lamproites (host rocks to diamonds). At 1?3 wt%, the chlorine levels were also much higher than those of kimberlites (<0.1%), although the bulk composition of the micro-inclusions, including the high K2O content (up to 29.7 wt%), resembled that of such potassic magmas. They concluded that these micro-inclusions represent the volatile-rich (~40% volatiles) fluid or melt from the upper mantle in which the diamonds grew, and that because of the high volatile content in this hydrous mantle fluid, high levels of rare gases may also be expected and explain the high 40Ar/K ratios (the excess 40Ar*) and anomalous "ages".

As a result of continued investigation of the Zaire cubic diamonds, which produced 40Ar/39Ar age spectra yielding a ~5.7 Ga isochron, Ozima et al. [74] discovered that just as there was an excellent correlation between their potassium contents and 40Ar/36Ar ratios, there is also a correlation between their chlorine contents and 40Ar. They concluded from their data "that the 40Ar is an excess component which has no age significance, and that the 40Ar and its associated potassium are contained in sub-micrometer inclusions of mantle-derived fluid." Turner et al. [93] also used the 40Ar/39Ar technique through correlations with K, Cl and 36Ar to unscramble the mixtures of radiogenic and parentless (excess) Ar components in fluid inclusions in "coated" Zaire diamonds and in olivine from an East African mantle xenolith. Their results proved conclusively that 40Ar is present in a widespread chlorine-rich component, which implies the existence of H2O/CO2-rich phases with 40Ar/Cl ratios that are "remarkably uniform over large distances", with enrichments of these two incompatible elements by almost four orders of magnitude relative to bulk upper-mantle values. Clearly, excess 40Ar* is abundant in the mantle and can be easily transported up into the crust.

Crustal Excess 40Ar*

Is there only evidence for excess 40Ar* in the mantle, gleaned from rocks (basalts and ultramafic xenoliths) and minerals (olivine, pyroxene, plagioclase and diamonds) that were formed in, or ascended from, the mantle? Patterson et al. [77] envisage noble gases from the mantle (and the atmosphere) migrating and circulating through the crust, so there should be evidence of excess 40Ar* in crustal rocks and minerals. In fact, noble gases in CO2-rich natural gas wells support such migration and circulation ?that is, the isotopic signatures clearly indicate a mantle origin for the noble gases, including amounts of excess 40Ar* in some CO2-rich natural gas wells exceeding those in the mantle-derived mid-ocean ridge basalts [6, 10, 66, 88, 89]. Staudacher [88] also notes that the quantities of excess 40Ar* in the continental crust can be as much as five times that found in such mantle-derived mid-ocean ridge basalts [89], strongly suggesting that excess 40Ar* in crustal rocks and their constituent minerals could well be the norm rather than the exception, thus making all K-Ar (and Ar-Ar) dating questionable.

It has now been established that some diamonds can form in the crust ? during high-grade metamorphism [22, 87] and via shock metamorphism during meteorite or asteroid impact [52]. The pressures and temperatures of high-grade metamorphism had been regarded as insufficient to produce diamonds, but the key ingredient was found to be volatile N2-CO2-rich fluids. Noble gas data on these diamonds are not yet available, due to their size and rarity, but such data have been definitive in establishing the crustal origin of carbonado diamonds [75]. Nevertheless, they still contain excess 40Ar*.

Dalrymple [17], referring to metamorphism and anatexis of rocks in the crust, commented, "If the rock is heated or melted at some later time, then some or all the 40Ar may escape and the K-Ar clock is partially or totally reset". In other words, 40Ar* escapes to migrate in the crust where it may then be incorporated in other minerals as excess 40Ar*, just as 40Ar* degassing from the mantle does. Thus, for example, excess 40Ar* has been recorded in many minerals (some of which contain no 40K) in crustal rocks, such as quartz, plagioclase, pyroxene, hornblende, biotite, olivine, beryl, cordierite, tourmaline, albite and spodumene [33, 61] ? in pegmatites, metamorphic rocks, and lavas. And it is not just K-Ar dating analyses that detect excess 40Ar*, as Lanphere and Dalrymple [59] used the 40Ar/39Ar method to confirm the presence of excess 40Ar* in feldspars and pyroxenes. Indeed, in a recent study [80], 128 Ar isotopic analyses were obtained from ten profiles across biotite grains in amphibolite-granulite facies metamorphic rocks, and apparent 40Ar/39Ar "ages" within individual grains ranged from 161 to 514 Ma. The investigators concluded that these observations cannot be solely due to radiogenic build-up of 40Ar*, but must be the result of incorporation by diffusion of excess 40Ar* from an external source, namely, 40Ar* from the mantle and other crustal rocks and minerals. Indeed, Harrison and McDougall [47] were able to calculate a well-defined law for 40Ar diffusion from hornblende in a gabbro due to heating. They also found that the excess 40Ar* which had developed locally in the intergranular regions of the host gabbro reached partial pressures in some places of at least 10-2 atm.

This crustal migration of 40Ar* is known to cause grave problems in attempted regional geochronology studies. In the Middle Proterozoic Musgrave Block of northern South Australia, Webb [96] found a wide scatter of K-Ar mineral ages ranging from 343 Ma to 4493 Ma due to inherited (or excess) 40Ar*, so that no meaningful interpretation could be drawn from the rocks (granulite, gneiss, pseudotrachylite, migmatite, granite and diabase). Of the diabase dikes which gave anomalous ages, he concluded that "The basic magmas probably formed in or passed through zones containing a high partial pressure of 40Ar*, permitting inclusion of some of the gas in the crystallizing minerals." Likewise, when Baski and Wilson [5] attempted to argon date Proterozoic granulite-facies rocks in the Fraser Range (Western Australia) and Strangways Range (central Australia), they found that garnet, sapphirine and quartz in those rocks contained excess 40Ar* that rendered their argon dating useless because of "ages" higher than expected. They also concluded that the excess 40Ar* was probably incorporated at the time of formation of the minerals, and their calculations suggested a partial pressure of ~0.1 atm Ar in the Proterozoic lower crust of Australia, which extends over half the continent.

In a detailed 40Ar/39Ar dating study of high-grade metamorphic rocks in the Broken Hill region of New South Wales (Australia), Harrison and McDougall [48] found evidence of widely distributed excess 40Ar*. The minerals most affected were plagioclase and hornblende, with step heating 40Ar/39Ar "age" spectra yielding results of up to 9.588 Ga. Such unacceptable "ages" were produced by excess40Ar* release, usually at temperatures of 350?650?C and/or 930?1380?C, suggesting the excess 40Ar* is held in sites within the respective mineral lattices with different heating requirements for its release. There are three principal trapping sites for Ar in solids ? structural holes, edge dislocations and lattice vacancies. (Argon is also known to be held sometimes in some minerals in fluid inclusions.) Clearly, this study shows that at crustal temperatures, which are less than 930?C, some excess 40Ar* will always be retained in those trapping sites in minerals where it is obviously "held" more tightly, thus rendering K-Ar and 40Ar/39Ar dating questionable. Harrison and McDougall [48] were only able to produce a viable interpretation of the data because they made assumptions about the expected age of the rocks and of a presumed subsequent heating event (based on Pb-Pb and Rb-Sr dating), the latter being the time when they conjecture that accumulated 40Ar* was released from minerals causing a significant regional Ar partial pressure of ~3x10-4 atm to develop.

Mantle-Crust Domains and Excess 40Ar*

Harte and Hawkesworth [50] have identified domains within the mantle and crust and described the interaction between them, all of which is relevant to the migration and circulation of argon (and thus excess 40Ar*) from the lower mantle to the crust and to lavas extruded on the Earth?s surface. The six domains are physically distinct units which show wide differences in average physical and chemical properties, as well as apparent age, structure, and tectonic behavior. They are the lower mantle (below 670 km), upper mantle, continental mantle lithosphere, oceanic mantle lithosphere, continental crust and oceanic crust, and each is a distinct geochemical reservoir. Each domain may provide material for magmatic rocks, and particular geochemical features of magmas may be associated with particular domains. Thus the convecting upper mantle which comes to the surface at mid-ocean ridges may be identified as the source of most geochemical features of mid-ocean ridge basalts, including their excess 40Ar* content. Similarly, the convecting lower mantle is regarded as the primordial or bulk Earth geochemical reservoir, which may also contribute excess 40Ar* to mid-ocean ridge basalts, but is more important for its contribution to ocean island basalts (e.g. Hawaii) and other plume-related basalts (continental alkali basalts and continental flood basalts). However, considerable complexity may be added to the deeper mantle geochemical structure as a result of localized accumulation of subducted oceanic lithosphere.

Porcelli and Wasserburg [82] have proposed a steady-state upper mantle model for mass transfer of rare gases, including argon. The rare gases in the upper mantle are derived from mixing of rare gases from the lower mantle, subducted rare gases, and radiogenic nuclides produced in situ. Porcelli and Wasserburg claim that all of the 40Ar in the closed-system lower mantle has been produced by 40K decay in the lower mantle, but this claim is based on the assumption of a 4.5 Ga Earth. In any case, they contradict themselves, because they also state [82, p. 4924], "The lower mantle is assumed to have evolved isotopically approximately as a closed system with the in situ decay of 129I, 244Pu, 238U, 232Th, and 40K adding to the complement of initial rare gases." In other words, they admit that some of the 40Ar must be primordial and not dervied from radioactive 40K. They then go on to claim that in the upper mantle, 40K decay further increases the radiogenic 40Ar from the lower mantle by a factor of ~3, but again this presupposes a 4.5 Ga Earth and doesn?t allow for primordial 40Ar that could well be also in the upper mantle if it?s admitted to be in the lower mantle.

In the case of the continental and oceanic lithospheric domains, the lack of convective stirring means that different geological processes and events may implant in each domain a variety of geochemically distinct materials which will remain isolated from one another. Therefore, these domains do not have a single set of geochemical characteristics; thus identification of geochemically defined "sources" with particular physically defined crust-mantle domains is complex, and the geochemical definition of particular reservoirs cannot be regarded as simply definition of major physical entities. Nevertheless, excess 40Ar* will be added to these domains by the passage of basaltic magma plumes from the upper mantle to the Earth?s surface.

Furthermore, the processes of oceanic lithosphere formation from the convecting upper mantle in association with mid-ocean ridge activity mean that its isotopic characteristics everywhere will be largely similar to those of the convecting upper mantle and mid-ocean ridge basalts, including the addition of excess 40Ar*. The corollary to this is that the oceanic crust is formed as part of these same processes. However, the oceanic crust generally has a thin veneer of sediments over it, and thick wedges of sediments adjacent to the domains of continental crust, whereas sections of oceanic crust are hydrothermally altered. The compositions of these components of the oceanic crust may, therefore, include a considerable contribution from continental detritus and ocean water, so that this oceanic crustal material may give rise to a distinct geochemical reservoir, the fate of which during subduction back into the upper mantle becomes critically important if it contributes to island arc volcanics, plume-related intra-place magmas and mantle-derived xenoliths.

The complexity of continental crustal material is well known through direct observation, and the mantle lithosphere attached to it may be expected to show a similar complexity. Nevertheless, it is evident that excess 40Ar* also resides in the continental mantle lithosphere, as indicated by xenoliths [63]. Likewise, there is evidence of excess 40Ar* in crustal magmatic rocks (e.g., gabbros [47], pegmatites [61]), migrating through metamorphic terrains [5, 48, 96], and in natural gas in sedimentary reservoirs [88].

Mt Ngauruhoe in its Tectonic Framework

The presence, therefore, of excess 40Ar* in the recent andesite flows at Mt Ngauruhoe is to be expected. The Taupo Volcanic Zone is a volcanic arc and marginal basin of the Taupo-Hikurangi arc-trench system (see Figure 1 again), which is a southward extension of the Tonga-Kermadec arc into the continental crustal environment of New Zealand?s North Island [13]. Geophysical investigations indicate that the Pacific Plate is being obliquely subducted beneath the Australian Plate on which most of New Zealand?s North Island sits, and that the volcanoes of the Taupo Volcanic Zone, including Ngauruhoe in the Tongariro Volcanic Center, are about 80 km directly above the subducting Pacific Plate, a zone of earthquakes revealing where the movement is taking place [97]. Friction along the plane of contact is believed to cause melting to produce pockets of magma, which then feed via conduits to the volcanoes above. Thus the recent andesite flows at Mt Ngauruhoe are calc-alkaline island arc volcanics.

The tectonic and geochemical framework of the Ngauruhoe andesite flows within the mantle-crust domains of Harte and Hawkesworth [50] is that of subducting oceanic crust (derived from the convecting upper mantle), carrying with it the wedge of continental sedimentary detritus which has accumulated at the continental margin and in the adjacent trench to the east of the coastline. Attached beneath the subducting oceanic crust is its associated oceanic mantle lithosphere, and together they are being thrust downwards into the upper mantle. Above the subducting plate are the continental crust and continental mantle lithosphere of the overriding plate, the continental crust being at the contact plane at shallow depths near the trench, and then the attached continental mantle lithosphere beneath at a depth of about 35 km [38]. Thus the geochemical reservoir from which the Ngauruhoe andesite magma has been drawn is potentially a mixture of melted oceanic crust, continental sedimentary detritus and continental crust, and possibly continental mantle lithosphere, or even upper mantle.

Genesis of the Mt Ngauruhoe Andesite Magma and its Excess 40Ar*

One of the easier investigations of the petrogenesis of these volcanic rocks of the Taupo Volcanic Zone was that of Stipp [91], and Ewart and Stipp [27]. They analysed samples that had been systematically collected, including not only the lavas and the pyroclastics, but also the Permian to Jurassic interbedded greywackes, siltstones and shales (the potential crustal source rocks) which are spatially related to, and underlie, the volcanics. Of primary interest were Sr, Rb and K contents, and 87Sr/86Sr and 87Rb/86Sr ratios. Three possibilities for the origin of the calc-alkaline andesite magma were under investigation ? fractional crystallization of a basalt magma under oxidizing conditions; some form of hybridization between basaltic and acidic magmas, possibly followed by fractional crystallization; and derivation of a primary andesite magma from the upper mantle. Ewart and Stipp [27] regarded their Sr isotopic data as more consistent with the production of the andesites by partial assimilation of sedimentary material by basaltic magma (derived from the upper mantle), the adjacent greywackes, siltstones and shales being the most likely sedimentary material, and the unassimilated gneissic xenoliths probably representing the basement rocks to those sediments. However, they admitted that the data did not exclude the possibility of a primary andesitic magma derived directly from the upper mantle, provided that some assimilation of crustal material modified it prior to eruption.

Subsequent investigations by Cole [12] favored the alternate petrogenetic model of a primary andesitic magma. He suggested that the subducting oceanic crust assimilated the greywacke-siltstone-shale and overlying sediments east of the Taupo Volcanic Zone to produce amphibolite, which subsequently broke down to produce phlogopite eclogite below 90 km. This in turn partially melted at 150?200 km, and the resultant magma fractionated in the upper mantle or lower crust to produce andesite. However, based on rare-earth element geochemistry, Cole et al. [14] modified that petrogenetic model, suggesting that while the andesite magma genesis was probably associated in the upper mantle with the downgoing slab and some crustal contamination occurred, the andesite does not appear to have had an eclogite parent. This would then suggest that the melting associated with the subducting slab to generate the andesite magma occurred at a depth of less than 90 km.

Graham and Hackett [38] agreed with this conclusion, demonstrating from geophysical evidence that the top of the subducting slab is at a depth of about 80 km below Ngauruhoe, and that the crust there is probably less than 20 km thick. Thus the upper mantle wedge between would consist only of plagioclase-peridotite and spinel-peridotite. At 80 km depth the hydrated amphibolite assemblage of the upper portion of the subducting slab of oceanic crust and oceanic mantle lithosphere would have started to dehydrate, thus liberating water and possibly other volatile constituents into the overlying upper mantle wedge, significantly lowering its melting point. Graham and Hackett [38] then showed that the geochemical evidence requires the andesite magma for the Ngauruhoe lava flows to have been generated from an original low-alumina basalt magma produced in the upper mantle wedge by anatexis of the asthenosphere (uppermost mantle) and/or subcontinental mantle lithosphere probably catalysed by hydrous, metasomatic fluids from the subducting slab.

Some specific geochemical enrichment then appears to have occurred as a result of this mantle metasomatism and continental crustal contamination during ascent and storage of the magma. Graham and Hackett [38] used least squares geochemical modelling to show how the andesite magma could be generated from such a parent basalt magma by a process of combined assimilation of crustal material (addition of 6 per cent assimilant) and fractional crystallization (30 per cent removal of crystals). Furthermore, the presence of xenoliths in the Ngauruhoe andesite flows, particularly the vitrified meta-greywacke and gneissic xenoliths, indicate conclusively that the assimilant was most likely a partial melt of gneiss, originally the adjacent greywacke-siltstone-shale sediments [38, 39].

These processes responsible for the generation of the andesite magma did not diminish the excess 40Ar* content of the resultant flows. Though the amount of excess 40Ar* is not high when compared with that found in mid-ocean ridge basalts, it is nonetheless significant that the excess 40Ar* was still present in the lavas upon eruption and cooling. The evidence indicates that the parent basaltic magma was generated in the upper mantle where the excess 40Ar* in the geochemical reservoir is now known to be upwards of 150 times more than the atmospheric content, relative to 36Ar. The subsequent crustal contamination and fractional crystallization to form the andesite magma during ascent, and the degassing of the magma during eruption and lava flow and cooling, did not remove all the excess 40Ar*, a small portion of which was left to be trapped in the congealed lava and its constituent minerals.

This model for the generation of the andesite magma in the post-Flood world is, of course, based on the plate tectonics model for global tectonics through Earth history. Even though the postulated plate movements today are extremely slow, and thus extrapolated back over millions of years by uniformitarians, a catastrophic model for plate tectonics in the context of the Flood is entirely compatible with both Scripture and the scientific data [4]. Plate movements are regarded as occurring catastrophically during the Flood and then rapidly slowing down to today?s rates in the post-Flood era.

CONCLUSIONS

The fact that there is even some excess 40Ar* in these recent andesite flows, and that it appears to have ultimately come from the upper mantle geochemical reservoir, where it is regarded as leftover primordial argon not yet fully expelled by the process of outgassing that is supposed to have occurred since the initial formation of the Earth, has very significant implications.

First, this is clearly consistent with a young Earth, where the very short time-scale since the creation of the Earth has been insufficient for all the primordial argon to be released yet from the Earth?s deep interior. Furthermore, it would also seem that even the year-long global catastrophic Flood, when large-scale convection and turdecer occurred in the mantle [4], was insufficient to expel all the deep Earth?s primordial argon.

Second, this primordial argon is, in part, "excess" 40Ar not generated by radioactive decay of 40K, which has then been circulated up into crustal rocks where it may continue migrating and building up to partial pressure status regionally. Because the evidence clearly points to this being the case, then when samples of crustal rocks are analysed for K-Ar "dating" the investigators can never really be sure that whatever 40Ar* is in the samples is from in situ radioactive decay of 40K since the formation of the rocks, or whether some or all of it is from the "excess 40Ar*" geochemical reservoirs in the lower and upper mantles. This could even be the case when the K-Ar analyses yield "dates" compatible with other radioisotopic "dating" systems and/or with fossil "dating" based on evolutionary assumptions. And there would be no way of knowing because the 40Ar* from radioactive decay of 40K cannot be distinguished analytically from primordial 40Ar not from radioactive decay, except of course by external assumptions about the ages of the samples.

Therefore, these considerations call into question all K-Ar "dating", whether "model ages" or "isochron ages", and all 40Ar/39Ar "dating", as well as "fossil dating" that has been calibrated against K-Ar "dates". Although seemingly insignificant in themselves, the anomalous K-Ar "model ages" for these recent andesite flows at Mt Ngauruhoe, New Zealand, lead to deeper questions. Why is there excess40Ar* in these rocks? From where did it come? Answers to these questions in turn point to significant implications that totally undermine such radioactive "dating" and that are instead compatible with a young Earth.

FUTURE RESEARCH

Further research is very definitely warranted. The most pressing need is to attempt to quantify how much primordial 40Ar there is today in the upper mantle. Also, how much has circulated into crustal rocks, how much is in natural gas reservoirs, and how much might have escaped into the atmosphere during 6,000?7,000 years, including accelerated rates during the Flood. It might then be possible to quantify how much primordial 40Ar there was in the mantle at the time of the Earth?s creation. From these calculations and associated modelling exercises there might develop quantifiable evidence for the Earth?s youthfulness.

Additionally, further research is needed to quantify how much "excess 40Ar*" is in all the crustal rocks and minerals that have been, and are, subject to K-Ar and 40Ar/39Ar "dating". This would include what are regarded as mantle xenoliths and xenocrysts (e.g. diamonds). It is helpful to show on the one hand that such "dating" is questionable, but on the other hand there are still many "dates" that are concordant, that is, they agree with other uniformitarian dating systems and schemes. So ultimately we need to explain why that happens when other "dates" are discordant and anomalous. There may, in fact, be some pattern or systematic way in which "excess40Ar" has been trapped in rocks and occluded in minerals at different levels (depths and relative ages) in the geological record. If so, then K-Ar and 40Ar/39Ar "dating" would irrevocably be discredited.

REFERENCES

[1] Armstrong, R.L., K-Ar Dating: Late Cenozoic McMurdo Volcanic Group and Dry Valley Glacial History, Victoria Land, Antarctica, New Zealand Journal of Geology and Geophysics, 21:6 (1978), pp. 685-698.
[2] Austin, S.A., Excess Argon within Mineral Concentrates from the New Dacite Lava Dome at Mount St Helens Volcano, Creation Ex Nihilo Technical Journal, 10:3 (1996), pp. 335-343.
[3] Austin, S.A. and Snelling, A.A., Discordant Potassium-Argon Model and Isochron "Ages" for Cardenas Basalt (Middle Proterozoic) and Associated Diabase of Eastern Grand Canyon, Arizona, this volume
[4] Austin, S.A., Baumgardner, J.R., Humphreys, D.R., Snelling, A.A., Vardiman, L. and Wise, K.P., Catastrophic Plate Tectonics: A Global Flood Model of Earth History, Proceedings of the Third International Conference on Creationism, R.E. Walsh, Editor, 1994, Creation Science Fellowship, Pittsburgh, PA, pp. 609-621.
[5] Baksi, A.K. and Wilson, A.F., An Attempt at Argon Dating of Two Granulite-Facies Terranes, Chemical Geology, 30 (1980), pp. 109-120.
[6] Ballentine, C.J., Resolving the Mantle He/Ne and Crustal 21Ne/22Ne in Well Gases, Earth and Planetary Science Letters, 152 (1997), pp. 233-249.
[7] Battey, M.H., The Recent Eruption of Ngauruhoe, Records of the Auckland Institute and Museum, 3 (1949), pp. 389-395.
[8] Broadhurst, C.L., Drake, M.J., Hagee, B.E. and Benatowicz, T.J., Solubility and Partitioning of Ar in Anorthite, Diopside, Fosterite, Spinel, and Synthetic Basaltic Liquids, Geochimica et Cosmochimica Acta, 54 (1990), pp. 299-309.
[9] Broadhurst, C.L., Drake, M.J., Hagee, B.E. and Benatowicz, T.J., Solubility and Partitioning of Ne, Ar, Kr and Xe in Minerals and Synthetic Basaltic Melts, Geochimica et Cosmochimica Acta, 56 (1992), pp. 709-723.
[10] Burnard, P., Graham, D. and Turner, G., Vesicle-Specific Noble Gas Analyses of "Popping Rock": Implications for Primordial Noble Gases in Earth, Science, 276 (1997), pp. 568-571.
[11] Clark, R.H., Petrology of the Volcanic Rocks of Tongariro Subdivision, The Geology of the Tongariro Subdivision, D.R. Gregg, New Zealand Geological Survey Bulletin n.s.40 (1960), pp. 107-123.
[12] Cole, J.W., Andesites of the Tongariro Volcanic Centre, North Island, New Zealand, Journal of Volcanology and Geothermal Research, 3 (1978), pp. 121-153.
[13] Cole, J.W. and Lewis, K.B., Evolution of the Taupo-Hikurangi Subduction System, Tectonophysics, 72 (1981), pp. 1-21.
[14] Cole, J.W., Cashman, K.V. and Rankin, P.C., Rare-Earth Element Geochemistry and the Origin of Andesites and Basalts of the Taupo Volcanic Zone, New Zealand, Chemical Geology, 38 (1983), pp. 255-274.
[15] Cole, J.W., Graham, I.J., Hackett, W.R. and Houghton, B.F., Volcanology and Petrology of the Quaternary Composite Volcanoes of Tongariro Volcanic Centre, Taupo Volcanic Zone, Late Cenozoic Volcanism in New Zealand, I.E.M. Smith, Editor, Royal Society of New Zealand Bulletin, 23 (1986), pp. 224-250.
[16] Dalrymple, G.B., 40Ar/36Ar Analyses of Historic Lava Flows, Earth and Planetary Science Letters, 6 (1969), pp. 47-55.
[17] Dalrymple, G.B., The Age of the Earth, 1991, Stanford University Press, Stanford.
[18] Dalrymple, G.B. and Lanphere, M.A., Potassium-Argon Dating: Principles, Techniques and Applications to G
ochronology, 1969, W.H. Freeman, San Francisco.
[19] Dalrymple, G.B. and Moore, J.G., Argon 40: Excess in Submarine Pillow Basalts from Kilaueau Volcano, Hawaii, Science, 161 (1968), pp. 1132-1135.
[20] Damon, P.E., Laughlin, A.W. and Precious, J.K., Problem of Excess argon-40 in Volcanic Rocks, Radioactive Dating Methods and Low-Level Counting, 1967, International Atomic Energy Agency, Vienna, pp. 463-481.
[21] Dickin, A.P., Radiogenic Isotope Geology, 1995, Cambridge University Press, Cambridge.
[22] Dobrzhinetskaya, L.F., Eide, E.A., Larsen, R.B., Sturt, B.A., Tr-, R.G., Smith, D.C., Taylor, W.R. and Posukhova, T.V., Microdiamond in High-Grade Metamorphic Rocks in the Western Gneiss Region, Norway, Geology, 23 (1995), pp. 597-600.
[23] Dymond, J., Excess Argon in Submarine Basalt Pillows, Geological Society of America Bulletin, 81 (1970), pp. 1229?1232.
[24] Eldridge, C.S., Compston, W., Williams, I.S., Harris, J.W. and Bristow, J.W., Isotope Evidence for the Involvement of Recycled Sediments in Diamond Formation, Nature, 353 (1991), pp. 649-653.
[25] Esser, R.P., McIntosh, W.C., Heizler, M.T. and Kyle, P.R., Excess Argon in Melt Inclusions in Zero-Age Anorthoclase Feldspar from Mt Erebus, Antarctica, as Revealed by the 40Ar/39Ar Method, Geochimica et Cosmochimica Acta, 61 (1997), pp. 3789-3801.
[26] Ewart, A., Notes on the Chemistry of Ferromagnesian Phenocrysts from Selected Volcanic Rocks, Central Volcanic Region, New Zealand Journal of Geology and Geophysics, 14:2 (1971), pp. 323-340.
[27] Ewart, A. and Stipp, J.J., Petrogenesis of the Volcanic Rocks of the Central North Island, New Zealand, as Indicated by a Study of Sr87/Sr86 Ratios, and Sr, Rb, K, U and Th Abundances, Geochimica et Cosmochimica Acta, 32 (1968), pp. 699-736.
[28] Faure, G., Principles of Isotope Geology, second edition, 1986, John Wiley and Sons, New York.
[29] Fisher, D.E., Heavy Rare Gases in a Pacific Seamount, Earth and Planetary Science Letters, 9 (1970), pp. 331-335.
[30] Fisher, D.E., Excess Rare Gases in a Subaerial Basalt from Nigeria, Nature, 232 (1971), pp. 60?61.
[31] Fisher, D.E., U/He Ages as Indicators of Excess Argon in Deep Sea Basalts, Earth and Planetary Science Letters, 14 (1972), pp. 255-258.
[32] Funkhouser, J.G. and Naughton, J.J., Radiogenic Helium and Argon in Ultramafic Inclusions from Hawaii, Journal of Geophysical Research, 73 (1968), pp. 4601-4607.
[33] Funkhouser, J.G., Barnes, I.L. and Naughton, J.J., Problems in the Dating of Volcanic Rocks by the Potassium-Argon Method, Bulletin of Volcanology, 29 (1966), pp. 709-717.
[34] Funkhouser, J.G., Fisher, D.E. and Bonatti, E., Excess Argon in Deep-Sea Rocks, Earth and Planetary Science Letters, 5 (1968), pp. 95100.
[35] Geyh, M.A. and Schleicher, H., Absolute Dating Methods, 1990, Springer-Verlag, Berlin.
[36] Gill, J.B., Orogenic Andesites and Plate Tectonics, 1981, Springer-Verlag, Berlin.
[37] Graham, I.J., Petrochemical and Sr Isotopic Studies of Lavas and Xenoliths from Tongariro Volcanic Centre - Implications for Crustal Contamination in Calc-Alkaline Magmas, Ph.D. Thesis (unpublished), Victoria University of Wellington (1985).
[38] Graham, I.J. and Hackett, W.R., Petrology of Calc-Alkaline Lavas from Ruapehu Volcano and Related Vents, Taupo Volcanic Zone, New Zealand, Journal of Petrology, 28:3 (1987), pp. 531-567.
[39] Graham, I.J., Cole, J.W., Briggs, R.M., Gamble, J.A. and Smith, I.E.M., Petrology and Petrogenesis of Volcanic Rocks from the Taupo Volcanic Zone: A Review, Journal of Volcanology and Geothermal Research, 68 (1995), pp. 59-87.
[40] Gregg, D.R., Eruption of Ngauruhoe 1954-55, New Zealand Journal of Science and Technology, B37 (1956), pp. 675-688.
[41] Gregg, D.R., The Geology of Tongariro Subdivision, New Zealand Geological Survey Bulletin, n.s.40 (1960).
[42] Grindley, G.W., Geological Map of New Zealand 1:250,000 Sheet 8 - Taupo, New Zealand Geological Survey (1960).
[43] Grindley, G.W., Tongariro National Park, Stratigraphy and Structure, New Zealand Department of Scientific and Industrial Research Information Series, 50 (1965), pp. 79-86.
[44] Gurney, J.J., Diamonds, Kimberlites and Related Rocks, Vol. 2: Their Mantle/Crust Setting, Diamonds and Diamond Exploration, Proceedings of the Fourth International Kimberlite Conference, 1986, Geological Society of Australia Special Publication No. 14, pp. 935?965.
[45] Gurney, J.J., The Diamondiferous Roots of Our Wandering Continent, South African Journal of Geology, 93 (1990), pp. 423-437.
[46] Hackett, W.R. and Houghton, B.F., Active Composite Volcanoes of Taupo Volcanic Zone, Central North Island Volcanism, New Zealand Geological Survey Record 21 (1987), pp. 61-114.
[47] Harrison, T.M. and McDougall, I., Investigations of an Intrusive Contact, Northwest Nelson, New Zealand - II. Diffusion of Radiogenic and Excess 40Ar in Hornblende Revealed by 40Ar/39Ar Age Spectrum Analysis, Geochimica et Cosmochimica Acta, 44 (1980), pp. 2005-2020.
[48] Harrison, T.M. and McDougall, I., Excess 40Ar in Metamorphic Rocks from Broken Hill, New South Wales: Implications for 40Ar/39Ar Age Spectra and the Thermal History of the Region, Earth and Planetary Science Letters, 55 (1981), pp. 123-149.
[49] Hart, S.R., K, Rb, Cs Contents and K/Rb, K/Cs Ratios of Fresh and Altered Submarine Basalts, Earth and Planetary Science Letters, 6 (1969), pp. 295-303.
[50] Harte, B. and Hawkesworth, C.J., Mantle Domains and Mantle Xenoliths, Kimberlites and Related Rocks, Vol. 2: Their Mantle/Crust Setting, Diamonds and Diamond Exploration, Proceedings of the Fourth International Kimberlite Conference, 1986, Geological Society of Australia Special Publication No. 14, pp. 649-686.
[51] Honda, M., McDougall, I., Patterson, D.B., Doulgeris, A. and Clague, D.A., Noble Gases in Submarine Pillow Basalt Glasses from Loihi and Kilauea, Hawaii: a Solar Component in the Earth, Geochimica et Cosmochimica Acta, 57 (1993), pp. 859?874.
[52] Hough, R.M., Gilmour, I., Pillinger, C.T., Langenhorst, F. and Montanari, A., Diamonds from the Iridium-Rich K-T Boundary Layer at Arroyo El Mimbral, Tamaulipas, Mexico, Geology, 25 (1997), pp. 1019-1022.
[53] Karpinskaya, T.B., Synthesis of Argon Muscovite, International Geology Review, 9 (1967), pp. 1493-1495.
[54] Karpinskaya, T.B., Ostrovskiy, I.A. and Shanin, L.L., Synthetic Introduction of Argon into Mica at High Pressures and Temperatures, Isv Akad Nauk S.S.S.R. Geology Series, 8 (1961), pp. 87-89.
[55] Kesson, S.E. and Ringwood, A.E., Slab-Mantle Interactions 2. The Formation of Diamonds, C
emical Geology, 78 (1989), pp. 97-118.
[56] Kirkley, M.B., Gurney, J.J. and Levinson, A.A., Age, Origin and Emplacement of Diamonds: A Review of Scientific Advances in the Last Decade, Canadian Institute of Mining Bulletin, 84 (1992), pp. 48-57.
[57] Kirkley, M.B., Gurney, J.J., Otter, M.L., Hill, S.J. and Daniels, L.R., The Application of C Isotope Measurements to the Identification of the Sources of C in Diamonds: A Review, Applied Geochemistry, 6 (1991), pp. 477-494.
[58] Krummenacher, D., Isotopic Composition of Argon in Modern Surface Volcanic Rocks, Earth and Planetary Science Letters, 8 (1970), pp. 109-117.
[59] Lanphere, M.A. and Dalrymple, G.B., Identification of Excess 40Ar by the 40Ar/39Ar Age Spectrum Technique, Earth and Planetary Science Letters, 12 (1976), pp. 359-372.
[60] Laughlin, A.E., Poths, J., Healey, H.A., Reneau, S. and WoldeGabriel, G., Dating of Quaternary Basalts Using the Cosmogenic 3He and 14C Methods with Implications for Excess 40Ar, Geology, 22 (1994), pp. 135-138.
[61] Laughlin, A.W., Excess Radiogenic Argon in Pegmatite Minerals, Journal of Geophysical Research, 74 (1969), pp. 6684-6690.
[62] Marty, B. and Humbert, F., Nitrogen and Argon Isotopes in Oceanic Basalts, Earth and Planetary Science Letters, 152 (1997), pp. 101-112.
[63] Matsumoto, T., Honda, M., McDougall, I. and O?Reilly, S.Y., Noble Gases in Anhydrous Lherzolites from the Newer Volcanics, Southeastern Australia: a MORB-Like Reservoir in the Subcontinental Mantle, Research School of Earth Sciences Annual Report 1997, Australian National University, Canberra, pp. 156-157.
[64] McDougall, I., The Geochronology and Evolution of the Young Volcanic Island of Reunion, Indian Ocean, Geochimica et Cosmochimica Acta, 35 (1971), pp. 261-288.
[65] McDougall, I., Polach, H.A. and Stipp, J.J., Excess Radiogenic Argon in Young Subaerial Basalts from the Auckland Volcanic Field, New Zealand, Geochimica et Cosmochimica Acta, 33 (1969), pp. 1485-1520.
[66] Moreira, M., Kunz, J. and Allegre, C., Rare Gas Systematics in Popping Rock: Isotopic and Elemental Compositions in the Upper Mantle, Science, 279 (1998), pp. 1178-1181.
[67] Nairn, I.A., Atmospheric Shock Waves and Condensation Clouds from Ngauruhoe Explosive Eruptions, Nature, 259 (1976), pp. 190-192..
[68] Nairn, I.A. and Self, S., Explosive Eruptions and Pyroclastic Avalanches from Ngauruhoe in February 1975, Journal of Volcanology and Geothermal Research, 3 (1978), pp. 39-60.
[69] Nairn, I.A. and Wood, C.P., Active Volcanoes of Taupo Volcanic Zone, Active Volcanoes and Geothermal Systems, Taupo Volcanic Zone, New Zealand Geological Survey Record 22 (1987), pp. 5-84.
[70] Nairn, I.A., Hewson, C.A.Y., Latter, J.H. and Wood, C.P., Pyroclastic Eruptions of Ngauruhoe Volcano, Central North Island, New Zealand, 1974 January and March, Volcanism in Australasia, R.W. Johnson, Editor, 1976, Elsevier, Amsterdam, pp. 385-405
[71] Navon, O., Hulcheon, I.D., Rossman, G.R. and Wasserburg, G.J., Mantle-Derived Fluids in Diamond Micro-Inclusions, Nature, 335 (1988), pp. 784-789.
[72] Noble, C.S. and Naughton, J.J., Deep-Ocean Basalts: Inert Gas Content and Uncertainties in Age Dating, Science, 162 (1968), pp. 265-267.
[73] Ozima, M. and Zashu, S., Primitive Helium in Diamonds, Science, 219 (1983), pp. 1067-1068.
[74] Ozima, M., Zashu, S., Takigami, Y. and Turner, G., Origin of the Anomalous 40Ar-39Ar Age of Zaire Cubic Diamonds: Excess 40Ar in Pristine Mantle Fluids, Nature, 337 (1989), pp. 226-229.
[75] Ozima, M., Zashu, S., Tomura, K. and Matsuhisa, Y., Constraints from Noble-Gas Contents on the Origin of Carbonado Diamonds, Nature, 351 (1991), pp. 472-474.
[76] Patterson, D.B., Honda, M. and McDougall, I., Atmospheric Contamination: a Possible Source for Heavy Noble Gases in Basalts from Loihi Seamount, Hawaii, Geophysical Research Letters, 17 (1990), pp. 705-08.
[77] Patterson, D.B., Honda, M. and McDougall, I., The Noble Gas Cycle Through Subduction Systems, Research School of Earth Sciences, Annual Report 1992, Australian National University, Canberra, pp. 104-106.
[78] Patterson, D.B., Honda, M. and McDougall, I., Noble Gases in Mafic Phenocrysts and Xenoliths from New Zealand, Geochimica et Cosmochimica Acta, 58 (1994), pp. 4411-4427.
[79] Phillips, D., Onstott, T.C. and Harris, J.W., 40Ar/39Ar Laser-Probe Dating of Diamond Inclusions from the Premier Kimberlite, Nature, 340 (1989), pp. 460-462.
[80] Pickles, C.S., Kelley, S.P., Reddy, S.M. and Wheeler, J., Determination of High Spatial Resolution Argon Isotope Variations in Metamorphic Biotites, Geochimica et Cosmochimica Acta, 61 (1997), pp. 3809-3833.
[81] Podosek, F.A., Pier, J., Nitoh, O., Zashu, S. and Ozima, M., Normal Potassium, Inherited Argon in Zaire Cubic Diamonds, Nature, 334 (1988), pp. 607-609.
[82] Porcelli, D. and Wasserburg, G.J., Transfer of Helium, Neon, Argon, and Xenon Through a Steady-State Upper Mantle, Geochimica et Cosmochimica Acta, 59 (1995), pp. 4921-4937.
[83] Poths, J., Healey, H. and Laughlin, A.W., Ubiquitous Excess Argon in Very Young Basalts, Geological Society of America Abstracts with Programs, 25 (1993), p. A-462.
[84] Reesman, R., letters dated February 19, 1997 and January 26, 1998.
[85] Richardson, S.H., Gurney, J.J., Erlank, A.J. and Harris, J.W., Origin of Diamonds in Old Enriched Mantle, Nature, 310 (1984), pp. 198-202.
[86] Seidemann, D.E., Effects of Submarine Alteration on K-Ar Dating of Deep-Sea Igneous Rocks, Geological Society of America Bulletin, 88 (1977), pp. 1660-1666.
[87] Sobolev, N.V. and Shatsky, V.S., Diamond Inclusions in Garnets from Metamorphic Rocks: A New Environment for Diamond Formation, Nature, 343 (1990), pp. 742-746.
[88] Staudacher, T., Upper Mantle Origin of Harding County Well Gases, Nature, 325 (1987), pp. 605-607.
[89] Staudacher, T., Sarda, P., Richardson, S.H., Allegre, C.J., Sagna, I and Dimitrieu, L.V., Noble Gases in Basalt Glasses from a Mid-Atlantic Ridge Topographic High at 14?N: Geodynamic Consequences, Earth and Planetary Science Letters, 96 (1989), pp. 119-133.
[90] Steiner, A., Petrogenetic Implications of the 1954 Ngauruhoe Lava and Its Xenoliths, New Zealand Journal of Geology and Geophysics, 1 (1958), pp. 325-363.
[91] Stipp, J.J., The Geochronology and Petrogenesis of the Cenozoic Volcanics of the North Island, New Zealand, Ph.D. Thesis (unpublished), Australian National University, Canberra (1968).
[92] Topping, W.W., Tephrostratigraphy and Chronology of Late Quaternary Eruptives from the Tongariro Volcanic Centre, New Zealand, New Zealand Journal of Geology and Geophysics, 16 (1973), pp. 397-423.
[93] Turner, G., Burgess, R. and Bannon, M., Volatile-Rich Mantle Fluids Inferred from Inclusions in Diamond and Mantle Xenoliths, Nature, 344 (1990), pp. 653-655.
[94] Valbracht, P.J., Honda, M., Matsumoto, T., Mattielli, N., McDougall, I., Ragettli, R. and Weis, D., Helium, Neon and Argon Isotope Systematics in Kerguelen Ultramafic Xenoliths: Implications for Mantle Source Signatures, Earth and Planetary Science Letters, 138 (1996), pp. 29-38.
[95] Valbracht, P.J., Staudigel, H., Honda, M., McDougall, I. and Davies, G.R., Isotopic Tracing of Volcanic Source Regions from Hawaii: Decoupling of Gaseous from Lithophile Magma Components, Earth and Planetary Science Letters, 144 (1996), pp. 185-198.
[96] Webb, A.W., Geochronology of the Musgrave Block, Mineral Resources Review, South Australia, 155 (1985), pp. 23-27.
[97] Williams, K., Volcanoes of the South Wind: A Field Guide to the Volcanoes and Landscape of the Tongariro National Park, Tongariro Natural History Society, Turangi, New Zealand (1994).
[98] Zashu, S., Ozima, M. and Nitoh, O., K-Ar Isochron Dating of Zaire Cubic Diamonds, Nature, 323 (1986), pp. 710-712.

Additional Resources:

Radioisotopes and the Age of the Earth by Vardiman, Snelling, Austin, Chaffin, DeYoung, Baumgardner (2001)

Impact #309 Potassium-Argon and Argon-Argon Dating of Crustal Rocks
and the Problem of Excess Argon by Andrew A. Snelling (Mar. 1999)

Impact #307 "Excess Argon": The "Achilles' Heel" of Potassium-Argon and
Argon-Argon "Dating" of Volcanic Rocks by Andrew A. Snelling, Ph.D. (Jan. 1999)

The Mythology of Modern Dating Methods by John Woodmorappe (1999)

Excess Argon within Mineral Concentrates from the New Dacite Lava Dome at
Mount St. Helens Volcano by STEVEN A. AUSTIN, Creation Ex Nihilo Technical Journal - Vol. 10 Part 3 (1996)