<figure><img src="image_1.png" /></figure>

HAL
open science

# Cosmogenic 3He chronology of postglacial lava flows at Mt Ruapehu, Aotearoa / New Zealand

Pedro Doll, Shaun Robert Eaves, Ben Matthew Kennedy, Pierre-Henri Blard, Alexander Robert Lee Nichols, Graham Sloan Leonard, Dougal Bruce Townsend, Jim William Cole, Chris Edward Conway, Sacha Baldwin, et al.

► To cite this version:

Pedro Doll, Shaun Robert Eaves, Ben Matthew Kennedy, Pierre-Henri Blard, Alexander Robert Lee Nichols, et al.. Cosmogenic 3He chronology of postglacial lava flows at Mt Ruapehu, Aotearoa / New Zealand. Geochronology, 2024, 6 (3), pp.365-395. 10.5194/gchron-6-365-2024. insu-04783358

HAL Id: insu-04783358

https://insu.hal.science/insu-04783358v1

Submitted on 14 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

# Cosmogenic ³He chronology of postglacial lava flows at Mt Ruapehu, Aotearoa / New Zealand

Pedro Doll<fnref n="1" />, Shaun Robert Eaves<fnref n="2" />,<fnref n="3" />, Ben Matthew Kennedy<fnref n="1" />, Pierre-Henri Blard<fnref n="4" />, Alexander Robert Lee Nichols<fnref n="1" />, Graham Sloan Leonard<fnref n="5" />, Dougal Bruce Townsend<fnref n="5" />, Jim William Cole<fnref n="1" />, Chris Edward Conway<fnref n="6" />, Sacha Baldwin<fnref n="1" />,★, Gabriel Fénisse<fnref n="4" />,★, Laurent Zimmermann<fnref n="4" />,★, and Bouchaïb Tibari<fnref n="4" />,★

**Correspondence:** Pedro Doll (pedro.doll@pg.canterbury.ac.nz)

Received: 18 January 2024 – Discussion started: 23 January 2024

Revised: 31 March 2024 – Accepted: 10 May 2024 – Published: 10 July 2024

**Abstract.** Accurate volcanic hazard assessments rely on a detailed understanding of the timing of past eruptions. While radiometric methods like $^{40}$Ar/$^{39}$Ar or K/Ar are by far the most conventional lava flow dating tools, their low resolution for young (< 20 ka) deposits interferes with the development of precise chronologies of recent effusive activity on most volcanoes. Mt Ruapehu (Aotearoa / New Zealand) has produced many lava flows throughout its history, but the precise timing of many recent eruptions remains largely unknown. In this study, we use cosmogenic ³He exposure dating to provide 23 eruption ages of young lava flows at Ruapehu. We then compare our results with existing $^{40}$Ar/$^{39}$Ar and paleomagnetic constraints, highlighting the value of cosmogenic nuclide exposure dating in refining recent eruptive chronologies. Of the 23 sampled flows, 16 provided robust eruption ages (5%–20% internal $2\sigma$; $n \ge 3$) between ca. 20 and 8 ka, except for one lava flow that erupted at around 43 ka, and their age distribution indicates that, during the last 20 kyr, effusive activity at Ruapehu peaked at 17–12 ka and at 9–7.5 ka. Nearly identical eruption ages of lavas located in different flanks of the volcanic edifice suggest concurrent activity from multiple vents during relatively short time intervals

(0–2 kyr) at around 13, 10, and 8 ka. We analysed four individual lava flows previously dated by $^{40}$Ar/$^{39}$Ar, two of which yield eruption ages older than the older limit of the $2\sigma$ interval of the radiometric dates, but the good clustering of individual samples from our sites suggests that our results better represent the real eruption age of these flows. Our ³He-based chronology shows excellent agreement with paleomagnetic constraints, suggesting that production rate uncertainties are unlikely to impact the accuracy of our eruption ages. This study demonstrates how cosmogenic nuclide dating can provide greater detail on the recent effusive chronology of stratovolcanoes, helping to resolve the low resolution of and difficulty in applying radiometric dating methods to young lava flows.

## 1 Introduction

Effusive volcanism is the main mechanism driving edifice growth on stratovolcanoes and poses a great hazard to the infrastructure, the natural environment, and the social fabric

<footnote n="1">¹School of Earth and Environment, University of Canterbury, Private Bag 4800, Ōtautahi / Christchurch 8041, New Zealand</footnote>

<footnote n="2">²Antarctic Research Centre, Victoria University of Wellington, P.O.Box 600, Te Whanganui-a-Tara / Wellington 6140, New Zealand</footnote>

<footnote n="3">³School of Geography, Environment and Earth Sciences, Victoria University of Wellington, P.O. Box 600, Te Whanganui-a-Tara / Wellington 6140, New Zealand</footnote>

<footnote n="4">⁴CRPG, CNRS, Université de Lorraine, 15 Rue Notre Dame des Pauvres, Vandoeuvre-les Nancy 54000, France</footnote>

<footnote n="5">⁵GNS Science, 1 Fairway Drive, Avalon, Te Awa Kairangi ki Tai / Lower Hutt 5011, New Zealand</footnote>

<footnote n="6">⁶Geological Survey of Japan, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Higashi, Tsukuba, Ibaraki 305-8567, Japan</footnote>

<footnote>★These authors contributed equally to this work.</footnote>
and livelihoods of local communities (Trusdell, 1995; Wilson et al., 2014; Harris, 2015; Jenkins et al., 2017; Tsang and Lindsay, 2020). Accurate hazard assessments rely on precise knowledge of recent eruption footprints, magnitudes, and frequencies (Connor et al., 2015) and, hence, accurate dating of eruptive events.

Most chronological studies of lava flows on stratovolcanoes are based on radiometric methods such as $^{40}$Ar/$^{39}$Ar and K/Ar. Recent advances in these methods (Coble et al., 2011; Fleck et al., 2014; Clay et al., 2015) have improved the precision of age determinations for Pleistocene lavas. However, errors in the ages of young ($< 20$ ka) products are still too large to precisely resolve recent eruptive chronologies (e.g. Wijbrans et al., 2011; Conway et al., 2016; Ramos et al., 2016; Calvert et al., 2018; Preece et al., 2018; Pure et al., 2020), hindering our ability to discriminate distinct eruptive episodes or to determine temporal relationships between ef- fusive eruptions and other volcanic processes. If available, radiocarbon dating of burned coal beneath lava flows can pro- vide accurate eruption ages, and it has been used widely in Hawai'i (e.g. Buchanan-Banks et al., 1989; Trusdell, 1995; see also Lockwood and Lipman, 1980) as well as in vari- ous volcanic regions (e.g. Moore and Rubin, 1991; Mishra et al., 2019; Sherrod et al., 2006). However, the use of ra- diocarbon is limited to areas with sufficient vegetation at the time of lava flow emplacement, so it is not applicable at high elevations or in periglacial environments. Alternative meth- ods, such as paleomagnetism or cosmogenic nuclide expo- sure dating, can support radiometric studies in non-vegetated areas and considerably reduce $^{40}$Ar/$^{39}$Ar and K/Ar uncer- tainties for Late Pleistocene and Holocene products (Sher- rod et al., 2006; Parmelee et al., 2015; Wright et al., 2015; Greve et al., 2016) and are therefore important when gener- ating more accurate eruptive histories in a wider spectrum of volcanic environments.

Cosmogenic nuclides are isotopes that originate when primary and secondary cosmic rays interact with atomic nuclei (Leya et al., 2000; Dunai, 2010). Some of them (terrestrial in situ cosmogenic nuclides, or TCNs) are formed in the upper few metres of the Earth's surface and can be used to calculate exposure ages of geological deposits provided they are rare in geological materials, are produced and retained in common minerals, are able to be analysed with reasonable confidence, are stable or have a half-life comparable to the timescales of the studied process, and have a well-understood origin and known relative contributions of their production mechanisms (Dunai, 2010). The number of TCNs that fulfil these requirements and have well-established methodologies developed for Earth science applications is relatively small (see Dunai, 2010), and the production rates and retention efficiency of TCNs vary across different minerals. $^3$He is a stable isotope with the highest production rate of all TCNs and a low detection limit in several geological settings (Blard, 2021), which makes it the ideal nuclide for dating young lava flows (Gosse and Phillips, 2001). This gas suffers diffusion loss in felsic

minerals (e.g. quartz and feldspars, and in volcanic groundmass containing them; Lippolt and Weigel, 1988; Tremblay et al., 2014) at Earth's surface temperatures, and it is therefore not normally used for silicic lithologies, which are better studied using $^{10}$Be or $^{26}$Al (e.g. Klein et al., 1986; Nishizumi et al., 1991; Smith et al., 2005). $^3$He is more efficiently retained in olivines and pyroxenes (Kurz, 1986a; Gosse and Phillips, 2001; Shuster et al., 2004; Blard, 2021), so it is suitable for dating volcanic eruptions (e.g. Kurz et al., 1990; Foeken et al., 2009; Parmelee et al., 2015), reconstructing glacial histories (e.g. Cerling and Craig, 1994; Blard et al., 2007) and fault kinematics (e.g. Fenton et al., 2001), or estimating erosion rates (e.g. Ferrier et al., 2013; Puchol et al., 2017), considering that the studied rocks contain these minerals.

Surface exposure dating using TCNs is applicable to geological deposits that have been brought to the surface and remained exposed to the cosmic ray flux ever since, provided there is no significant erosion or shielding (glacial, snow, debris, soil, tephra, or vegetation cover) that could have affected their cosmogenic nuclide inventory. In temperate climates, suitable sites will lie at elevations between the vegetation limit and where cryogenic processes begin to dominate. In dynamic environments such as stratovolcanoes, original surfaces are more likely to be preserved on younger lava flows, which have had a relatively limited time exposed to erosive and/or depositional processes. In addition, flow interiors with crystalline groundmass necessary for $^{40}$Ar/$^{39}$Ar or K/Ar dating are less likely to be exposed in young lava flows for the same reason. For young lava flows, cosmogenic $^3$He ($^3$He$_{cos}$) has the potential to resolve events down to 100 years under the most favourable conditions (low magmatic He and eruption ages $\le 10$ ka; Niedermann, 2002) and commonly yields ages with uncertainties of 15%–20% (2$\sigma$ including production rate errors), which is significantly more precise than traditional radiometric techniques for lavas $< 20$ ka (e.g. Wijbrans et al., 2011; Calvert et al., 2018; Pure et al., 2020). Thus, $^3$He$_{cos}$ can be used to complement chronological studies by providing greater detail of recent construction histories of volcanic edifices (e.g. Kurz et al., 1990; Foeken et al., 2009; Espanon et al., 2014; Parmelee et al., 2015). However, most of this research is focused on basaltic lava flows in extensional environments (e.g. Kurz et al., 1990; Licciardi et al., 2007; Foeken et al., 2009; Espanon et al., 2014; Marchetti et al., 2014; Medynski et al., 2015), and the application of $^3$He$_{cos}$ on stratovolcanoes (e.g. Parmelee et al., 2015) is still limited. 'A'ā lavas (commonly found in andesitic stratovolcanoes) normally have prominent tumuli standing out from the landscape, which are less likely to accumulate large amounts of snow or tephra compared to flatter primary morphologies (e.g. ropy pāhoehoe surfaces targeted in basic lavas; Kurz et al., 1990; Espanon et al., 2014; Marchetti et al., 2014; Parmelee et al., 2015), making them ideal targets for surface exposure dating (see Licciardi et al., 2007).
<figure><img src="image_6.png" /></figure>

**Figure 1.** Location map of study area. (a) Te Ika-a-Māui/North Island (Aotearoa/New Zealand) with its main active volcanic areas detailed: AVF – monogenetic Auckland Volcanic Field, Wh – Whakaari/White Island, Mo – Motuhara, Pu – Putauaki, TVZ – Taupō Volcanic Zone, Tg – Tongariro, Ru – Ruapehu, and Tk – Taranaki. (b) Detail of the “Central Plateau” at the southern end of the TVZ: Tp – Taupō, WB – Waimarino basalt, Pa – Pihanga, Tg – Tongariro, TM – Te Maari, Nga – Ngauruhoe, RC – Red Crater, Hh – Hauhungatahi, and OB – Ohakune Formation basalt. (c) Study area, with Ruapehu’s postglacial and late synglacial lava units mapped after Townsend et al. (2017), and sampled sites in this study. The maximum glacial extent during the LGM (Last Glacial Maximum; ~20–15 ka) after Barrell (2011) is outlined by the dashed yellow line. Panel (c) shows the NV – northern vent and SV – southern vent. Abbreviations next to sampled sites refer to lava flow names; there is a full list in Table A1. (d) Photo of Ruapehu taken from the south, with the Mangaehuehu Glacier directly beneath Ruapehu’s summit (viewpoint’s location is shown in panel (c) with the white arrow labelled as “d”).
<table><caption>Table 1. Previous chronological studies from lava flows at Ruapehu.</caption><thead><tr><th>Formation</th><th>Member</th><th>Eruption ages (±2σ)</th><th>Methods</th><th><sup>3</sup>He-based eruption ages, this study (ka, ±2σ)</th></tr></thead><tbody><tr><td rowspan="10">Whakapapa<br>(< 15 ka;<br>postglacial)</td><td>Crater Lake (< 5 ka)</td><td>2400–2050 BP<sup>a</sup>, 0.2 ± 2.2 ka<sup>b</sup></td><td>Paleomagnetism (refined from <sup>40</sup>Ar/<sup>39</sup>Ar)</td><td>–</td></tr><tr><td>Iwikau<br>Tawhainui<br>flows (9–7 ka)</td><td>DC: 8200–7900 BP<sup>a</sup>,<br>6.0 ± 2.4 ka<sup>b</sup>;<br>BR: 8800–8500 BP<sup>a</sup></td><td>Paleomagnetism (refined from <sup>40</sup>Ar/<sup>39</sup>Ar and tephra stratigraphy)</td><td>DC: 7.8 ± 1.5; BR: 8.1 ± 2.1;<br>WG: ≥7.8 ± 2.4</td></tr><tr><td>Mangatoetoenui<br>flows<br>(< 17 ka*)</td><td>LC: 0.8 ± 5.6 ka;<br>TSb: 9.2 ± 8.0 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>LC: 11.4 ± 2.3; TSa: 9.4 ± 1.8;<br>TSb: 11.5 ± 2.2; TFt: ≥8.6 ± 2.6</td></tr><tr><td>Taranaki Falls<br>flow (11–9 ka*)</td><td>TFa: 10 800–8900 BP<sup>a</sup>,<br>8.8 ± 2.8 ka<sup>b</sup></td><td>Paleomagnetism (refined from <sup>40</sup>Ar/<sup>39</sup>Ar)</td><td>TFa: 14.6 ± 2.9</td></tr><tr><td>Saddle Cone (10–8 ka*)</td><td>SC: 9850–8650 BP<sup>a</sup></td><td>Paleomagnetism (refined from tephra stratigraphy)</td><td>SC: 9.9 ± 2.0; WP: ≥11.2 ± 2.2</td></tr><tr><td>Pinnacle Ridge (~ 10 ka*)</td><td>PR: ~ 10 ka<sup>c</sup></td><td>Correlation with tephra</td><td>PR: 20.4 ± 4.0</td></tr><tr><td>Tureiti (15–9 ka)</td><td>12.5 ± 2.6; 11.9 ± 2.8 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Rangataua (~ 15–10 ka*)</td><td>~ 15–10 ka<sup>d,e</sup></td><td>Stratigraphy</td><td>RTp: 13.6 ± 2.6; RTm: 15.8 ± 3.0</td></tr><tr><td>Paretetaitonga (~ 15 ka)</td><td>14.8 ± 3.0 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>WT: 13.3 ± 2.6</td></tr><tr><td>Turoa (17–10 ka*)</td><td>15.1 ± 2.4; 11.9 ± 2.2 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>MN: 8.3 ± 1.6; MS: ≥6.6 ± 1.7;<br>CTa: ≥13.6 ± 2.7; CTb: 8.6 ± 1.7;<br>TC: 13.4 ± 2.6</td></tr><tr><td rowspan="10">Mangawhero<br>(50–15 ka;<br>synglacial)</td><td>Makotuku (24–16 ka)</td><td>20.9 ± 2.8; 17.8 ± 2.2 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>MF: 12.6 ± 3.5; NR: 42.9 ± 8.6;<br>MA: ≥54.0 ± 18.0</td></tr><tr><td>Waitonga (25–21 ka)</td><td>23.0 ± 1.6<sup>b</sup>; 23 ± 8 ka<sup>f</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Te Piripiri (~ 21 ka)</td><td>21 ± 6 ka<sup>f</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Horonuku (29–15 ka)</td><td>23 ± 4; 22 ± 7 ka<sup>f</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Whakapapaiti (~ 26 ka)</td><td>25.7 ± 3.8 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Manganuioteao (36–22 ka)</td><td>25.7 ± 2.6; 27.2 ± 4.8;<br>30.7 ± 5.2; 30.9 ± 2.2 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Mananui (42–38 ka)</td><td>40.3 ± 2.2<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Te Kohatu (44–36 ka)</td><td>47.6 ± 1.4; 39.1 ± 1.4;<br>39.2 ± 2.0; 42.6 ± 1.8 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Mangaturuturu (46–36 ka)</td><td>38.4 ± 2.4; 41.3 ± 1.8;<br>43.4 ± 2.4 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Mangaehuehu (47–40 ka)</td><td>42.8 ± 1.0; 43.1 ± 1.4;<br>43.3 ± 1.6; 44.2 ± 1.8;<br>45.4 ± 2.0 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>GR: ≥14.2 ± 2.7</td></tr><tr><td>Ngahuinga (48–35 ka)</td><td>44.8 ± 3.0 ka<sup>b</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Waihianoa (166–80 ka)</td><td>88.1 ± 6.4; 95.9 ± 7.0;<br>120.7 ± 4.0; 121.4 ± 2.8;<br>121.7 ± 4.2; 133.6 ± 6.4<sup>b</sup>;<br>119 ± 12; 129 ± 15; 130 ± 23;<br>131 ± 27; 133 ± 11; 134 ± 12;<br>135 ± 14; 138 ± 14; 147 ± 10;<br>147 ± 12; 154 ± 12 ka<sup>f</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr><tr><td>Te Herenga (200–150 ka)</td><td>158.8 ± 8.2; 169.4 ± 7.8;<br>174.6 ± 3.4; 186.2 ± 6.8;<br>187.9 ± 34.4<sup>b</sup>; 183 ± 13;<br>197 ± 12; 205 ± 27 ka<sup>f</sup></td><td><sup>40</sup>Ar/<sup>39</sup>Ar</td><td>–</td></tr></tbody></table>

<footnote><sup>a</sup> Greve et al. (2016). <sup>b</sup> Conway et al. (2016). <sup>c</sup> Donoghue et al. (1999). <sup>d</sup> Price et al. (2012). <sup>e</sup> Eaves et al. (2016b). <sup>f</sup> Gamble et al. (2003). <sup>*</sup> Age limits redefined in this study based on <sup>3</sup>He<sub>cos</sub> eruption ages. Ages separated by commas represent the same lava flow and semicolons separate dates from different flows.</footnote>
age, varying between 16% and 23% for 20–11 ka deposits and 32% and 1000% for Holocene lavas (see Table 1).

Greve et al. (2016) refined the eruption age for ⁴⁰Ar/³⁹Ar-dated and tephra-constrained Holocene lava flows by comparing characteristic magnetization directions recorded in the lavas with a paleosecular variation record based on lake sediments from Mavora Lakes (Te Waipounamu/South Island, Aotearoa/New Zealand), which were independently calibrated using ¹⁴C (Turner et al., 2015). Dating lava flows using paleomagnetic directions, however, requires a previous eruption age constraint and is limited to the Holocene in Aotearoa/New Zealand due to the extension of the sediment record. Only the ages of five flows were constrained using this method: one from the Crater Lake Member, three from the Iwikau Member (Delta Corner, Bruce Road, and Taranaki Falls flows), and one from the western lobe of the Saddle Cone Member. Eruption ages provided by Greve et al. (2016) for the Crater Lake, Delta Corner, and Bruce Road flows are tightly constrained (age ranges of ca. 300 years), while the preferred ages for the Taranaki Falls flow and Saddle Cone Member span ~2 and ~1.2 kyr, respectively (Table 1).

# 3 Methods

## 3.1 Sampling site selection

The selection of an adequate sampling site is an important step for cosmogenic nuclide exposure dating. Evidence of negligible erosion and confidence that the targeted rock has not been covered by other rocks, soil, ice, volcanic ash, or vegetation for a significant amount of time since formation are essential. For lava flows, effective sampling was achieved by targeting tumuli, spikes, and other features standing above the main flow surface (e.g. Anderson et al., 1994; Licciardi et al., 2006), which preserve characteristic primary cooling morphologies of flow surfaces (Fig. 2). Additional photos of sampled sites and examples of sites that are not suitable can be found in the Supplement file S1.

Using aerial photographs and digital elevation models (DEMs) based on aerial imagery and a newly acquired lidar dataset, we revised the existing maps (Townsend et al., 2017) and identified individual lava flows within each of the different members of the Whakapapa Formation, which we then targeted in our sampling. Lack of adequate lava surface exposures did not permit us to sample lavas from the Tureiti and Crater Lake members. Due to the lack of chronological data of several lavas of the Makotuku Member of the Mangawhero Formation (24–16 ka; Conway et al., 2016, Table 1), we additionally targeted three flows of this unit on outcrops outside of the LGM ice limits (MF, NR, and MA; Fig. 1c). We also sampled a site (GR) that we consider to be postglacial due to the presence of original (non-eroded) lava surfaces and its location inside the LGM ice limit of Barrell (2011, Fig. 1c). Note that this exposure was previously mapped as the Mangawhero Formation (Mangaehuehu Member)

based primarily on its location on the volcano and similarity in appearance to nearby geochemically fingerprinted outcrops.

## 3.2 Sample collection

All samples were collected under a research and collection permit of the Department of Conservation of Aotearoa/New Zealand, which was obtained after a consultation process involving local iwi with rightful claims to guardianship of Ruapehu. We sampled between three and six shallow surfaces (< 6 cm below the flow top) for each targeted flow using a hammer and chisel. For recording the coordinates and altitude of each surface, we used a differential Trimble Geo 7X GPS (vertical precision of 0.1 m) corrected by data of VGMT (Ohakune, Land Information New Zealand) and the Chateau Observatory base (GeoNet) stations. We also measured surface dip and orientation and azimuth–horizon angle pairs to account for topographic shielding. For the CTa, CTb, and TC samples, in situ topographic shielding could not be acquired, so representative azimuth–horizon angle pairs were selected based on observations of DEMs. To test the accuracy of this approach, we compared values derived from DEMs to field-obtained shielding factors from other sites, showing an agreement of 95%–99%.

## 3.3 Mineral separation

For each sample fragment used, the mean thickness was calculated using a caliper in 5–40 points, and then a sample thickness average was obtained, which was weighted by rock fragment mass. Afterwards, samples were crushed and sieved to obtain a 100–1000 µm size fraction, which was then rinsed to eliminate dust and organic matter and dried at 60 °C.

Density separation was done using a 3.0 g cm⁻³ sodium polytungstate solution, after which the heavy concentrates were leached in a 5% HF, 2.5% NaOH bath for 24 h before immersing them in 3 M HCl to remove fluoride precipitates, following Bromley et al. (2014). After checking under a microscope, we leached a second and third time, if necessary, in a 5% HF, 2.5% NaOH, and/or 2.5% HF, 1.25% NaOH solution, until we achieved total removal of groundmass on most crystals. We then carried out magnetic separation of oxides and magnetic groundmass, and then we finally visually removed the remaining impurities, based on colour and texture, to leave pure pyroxenes (olivines and pyroxenes in the GR samples).

## 3.4 Geochemical analyses

For each studied lava flow, major and trace element compositions were analysed at the Service d'Analyse des Roches et Minéraux (SARM) of the Centre de Recherches Pétrographiques et Géochimiques (CRPG, Université de Lorraine, Nancy, France) by inductively coupled plasma optical
<figure><img src="image_7.png" /></figure>

**Figure 2.** Examples of targeted sites. Red arrows point to a 20 cm long GPS. Panels (a), (b), (c), and (d) represent typical sampled surfaces ('a'ā morphologies). (a) A large tumuli standing out 1.5 m above ground level (sample RTp-PD027). (b) A large tumuli standing out 2.5 m above ground level (sample DC-PD330). (c) Detail of lava top with rough, irregular surfaces resembling 'a'ā lava flow morphologies, which is indicative of minimal erosion (sample MN-PD220). (d) 'A'ā block standing out ~40 cm above the ground with a pencil for scale (sample SC-PD001). (e) Surface of the Pinnacle Ridge spatter deposit (sample PR-PD085). (f) Sampled surface of the Waihohonu Plateau blocky flow (sample WP-PD008).

emission spectroscopy (ICP-OES) and inductively coupled plasma mass spectrometry (ICP-MS), respectively, for both bulk rock and pure pyroxenes and olivines. Each analysed sample consisted of 1 g of powdered rock/minerals that was fused at 980 °C for 60 min in Pt crucibles together with ultra-pure LiBO₂ in a 1 : 3 ratio prior to glass dissolution and measurements. The complete procedure is described in detail in Carignan et al. (2001).

## 3.5 Measurement of helium isotope concentrations

We analysed ³He and ⁴He concentrations in pyroxenes and olivines using a GV Instruments Helix Split Flight Tube multi-collector noble gas mass spectrometer attached to a gas line at CRPG (e.g. Schimmelpfennig et al., 2011; Blard et al., 2013, 2015).

Pure minerals were wrapped in tin capsules, loaded in a carousel, and baked for one night at 100 °C under ultra-high vacuum. The samples were heated to > 1300 °C for
15 min in a full metal induction furnace (Zimmermann et al., 2018), and the expelled gases were purified using four activated charcoal traps at 77 K in order to trap large amounts of CO₂, H₂O, N₂, and heavy noble gases (Ar, Kr, and Xe) from the melted samples by physisorption. In parallel, four getters initially activated at 800 °C were used at room temperature to trap all reactive species (e.g. H₂O, CO₂, N₂, and O₂) by chemisorption. After these two steps, He was condensed using a cryogenic trap at 12 K under ultra-low pressure (0.5–1 × 10⁻⁸ mbar) and then released at 75 K towards the mass spectrometer that measured, in static mode, ³He and ⁴He. The source settings were adjusted to get the best compromise between linearity, sensibility, and stability (e.g. ³He sensitivity = 4.30 × 10¹⁸ ± 5 % cps mol⁻¹ and ⁴He sensitivity = 7.45 × 10¹³ ± 2 % mV mol⁻¹). HESJ gas standards (20.63 R/Ra, Matsuda et al., 2002, Ra: atmospheric ³He/⁴He ratio of 1.39, Lupton and Evans, 2013) were measured daily with a reproducibility of 4.7 %, and ⁴He and ³He values were also routinely compared with CRONUS-P standards (Blard et al., 2015; Schaefer et al., 2016, reproducibility of 5.0 %). The main source of background He (measured daily with typical ³He blanks < 5 ± 3.5 × 10³ atoms, typical ⁴He blanks of 1.3 × 10⁹ ± 1.8 × 10⁸ atoms, and ³He/⁴He ratios similar to 1 Ra) was the Ta crucible, which was degassed at 1800 °C for 30 min prior to the analyses.

Crushed-released He isotopic analyses (used for magmatic corrections) were performed in samples with larger crystals (dominant fraction of 500–1000 µm, which were shown to contain larger amounts of magmatic He likely hosted in inclusions; Puchol et al., 2017) using a soft iron slug activated by external solenoids. Samples were crushed for 5 to 7 min at 100 strokes per minute with tube-specific ³He blanks between 3.8 ± 1.1 and 0.6 ± 0.3 × 10⁴ atoms and ⁴He blanks between 3.1 ± 0.1 × 10⁹ and 2.0 ± 1.8 × 10⁸ atoms. For a detailed description of the in vacuo crushing He extraction method, see Puchol et al. (2017).

## 3.6 Surface exposure age determinations

### 3.6.1 Calculation of cosmogenic ³He

To correctly determine the concentration of ³He_cos, it is necessary to consider the non-cosmogenic contributions to total ³He measured when fused in vacuo (³He_tot), which is described as

$$ ^{3}\text{He}_{\text{tot}} = ^{3}\text{He}_{\text{cos}} + ^{3}\text{He}_{\text{atm}} + ^{3}\text{He}_{\text{nuc}} + ^{3}\text{He}_{\text{mag}}, \quad (1) $$

where ³He_atm is the atmospheric ³He hosted at the minerals' surfaces as a contaminant and is time-independent. ³He_nuc is the nucleogenic ³He produced by capture of low-energy neutrons emitted by ⁶Li and dependent on Li concentrations in the mineral, U and Th concentrations in the rock, and the mineral closure age (equivalent to eruption age for pyroxenes and olivines in volcanic rocks, Kurz, 1986a). ³He_mag is

the magmatic ³He contribution (time-independent) present in melt and fluid inclusions and within the matrix of the minerals.

Atmospheric He (both ³He and ⁴He) concentrations are inversely proportional to the mineral grain size and become insignificant for minerals larger than 100 µm (Blard, 2021), so they were considered non-existent in our calculations. ³He_nuc quotas are normally negligible for uneroded lava flows in which the closure and exposure ages are the same (Kurz, 1986a), as shown by our calculations (Table A3) based on the spreadsheet developed by Blard (2021).

The total contribution of ³He_mag was accounted for in Eq. (3) and estimated using a magmatic ratio obtained as an uncertainty-weighted average from isotopic analyses of three samples crushed in vacuo and previous data from pyroxene and olivine phenocrysts in the Waimarino and Ohakune basalts (Sect. 4.2 and Supplement files S2.1 and S2.2).

The total amount of ⁴He measured in each sample (⁴He_tot) is defined by the following equation:

$$ ^{4}\text{He}_{\text{tot}} = ^{4}\text{He}_{\text{mag}} + ^{4}\text{He}_{\text{atm}} + ^{4}\text{He}_{\text{rad}} + ^{4}\text{He}_{\text{cos}}, \quad (2) $$

where ⁴He_mag corresponds to the time-independent magmatic ⁴He quota naturally present in the minerals, while ⁴He_atm accounts for atmospheric ⁴He contaminating the minerals' surfaces (time-independent). ⁴He-rad is generated by the decay of radioactive isotopes present in the minerals (such as U, Th, and Sm) and dependent on the abundance of these elements in the minerals and the closure age. Crystals normally exhibit an enriched ⁴He exterior rim generated by implanted ⁴He-rad from the matrix (Lal, 1989), typically with higher concentrations of U, Th, and Sm. ⁴He_cos refers to the cosmogenic contribution of ⁴He, which is negligible compared to other non-cosmogenic varieties of ⁴He (Blard, 2021) and is therefore also omitted from our calculations.

In this paper, we follow the approach of Blard and Farley (2008), which corrects for the contributions of ⁴He-rad, ⁴He_mag, and ³He_mag for uneroded lava flows using the following equation:

$$ ^{3}\text{He}_{\text{cos}} = \frac{^{3}\text{He}_{\text{tot}} - ^{4}\text{He}_{\text{tot}}\left(\frac{^{3}\text{He}}{^{4}\text{He}}\right)_{\text{mag}}}{R}, \quad (3) $$

where R (or R factor) is defined by

$$ R = 1 - \left(\frac{P_4}{P_3}\right)\left(\frac{^{3}\text{He}}{^{4}\text{He}}\right)_{\text{mag}}, \quad (4) $$

where P₄ and P₃ are the ⁴He-rad and local ³He_cos production rates.

The use of the R factor is essential when using ³He_cos to date uneroded lava flows as it permits the incorporation of a time-dependent ⁴He-rad quota, avoiding the issue of an underestimation or overestimation of the ⁴He_mag (and, hence, ³He_mag) contribution.
Individual values of $P_4$ were calculated for each lava flow using the spreadsheet developed by Blard (2021), neglecting the implanted $ ^{4}He_{rad} $ component to account for the removal of the $ ^{4}He $-enriched crystal rim with HF leaching.

Sample-specific $P_3$ estimates were obtained following the Lal–Stone time-corrected scaling scheme (Lal, 1991; Stone, 2000; Nishiizumi, 1989; Balco et al., 2008) using the online calculator “Cosmic Ray Exposure program” (CREp; https://crep.otelo.univ-lorraine.fr/, last access: 25 May 2024; Martin et al., 2017) and the global $ ^{3}He_{cos} $ production rate database therein.

### 3.6.2 Determination of exposure and eruption ages

To obtain exposure ages, we used the CREp online calculator, which calculated exposure ages based on our $ ^{3}He_{cos} $ concentrations and scaling parameters, the Lal–Stone time-corrected scaling scheme (Lal, 1991; Stone, 2000; Nishiizumi, 1989; Balco et al., 2008), the ERA-40 atmosphere model (Uppala et al., 2005), the geomagnetic framework of Muscheler et al. (2005), and the worldwide mean $ ^{3}He_{cos} $ production rates of $122 \pm 12$ at $g^{-1}yr^{-1}$ at sea level and high latitudes (SLHL).

Exposure ages calculated using the LSD scaling scheme (Lifton et al., 2014) and different atmospheric models and geomagnetic databases are available in the Supplement (file S3), showing variations of 1%–3% compared with the exposure ages calculated using the parameters outlined above. This is, however, not the case of the LSD geomagnetic framework, which provides exposure ages between 8.6% and 3.8% younger. This discrepancy can be explained by a higher spatial variability in the LSD framework than other models and especially by the model’s relative scaling factor high over the Aotearoa/New Zealand region during the Holocene (Lifton, 2016). New paleosecular variation records based on Aotearoa/New Zealand lake sediment cores (Turner et al., 2015; Turner and Corkill, 2023) suggest that this scaling factor high is a spatial artefact caused by the small number of Southern Hemisphere records used to make up the global model of the LSD framework. Thus, we place greater emphasis on results produced using models that do not contain such effects (e.g. Muscheler et al., 2005; Lifton, 2016).

$ ^{3}He_{cos} $ production rates have been shown to be indistinguishable in clinopyroxenes and orthopyroxenes (Delunel et al., 2016), justifying our decision to use a worldwide mean production rate estimate for our exposure age determinations. Additionally, this production rate value is supported by a local calibration test using the radiocarbon-dated debris avalanche deposits of the Murimotu Formation on the outer northwestern slopes of Ruapehu (Eaves et al., 2015). Despite some studies that suggested that olivines concentrate slightly larger amounts of $ ^{3}He_{cos} $ compared to pyroxenes (Ackert et al., 2003; Fenton et al., 2009), the difference was almost statistically insignificant, and, in a more recent study, Fenton and Niedermann (2014; as well as previous data from Blard

et al., 2006) provided results implying that olivine and pyroxenes have similar amounts of $ ^{3}He_{cos} $.

We measured three to five samples per lava flow to counter the possibility that individual samples may be affected by erosion or shielding that would compromise their accuracy for constraining the time of lava flow emplacement. To derive single exposure ages for lava flows from these multiple measurements, we used each sample’s internal age uncertainty ($1\sigma$, not including the external uncertainty from $P_3$) and implemented the summary age statistics and outlier removal routine contained in version 3 of the Balco et al. (2008) online exposure age calculator, described in the documentation (Sect. 4.C, available at https://sites.google.com/a/bgc.org/v3docs/, last access: 25 May 2024). In summary, we used weighted mean summary ages if the samples formed a single population at the 95% confidence interval using the chi-squared statistic (outliers are indicated in italics throughout Sect. 4.4). If this result could not be achieved by incremental outlier removal while maintaining a sample population $\ge$ 3, then we used the mean and standard error as the summary age of the lava flow. We finally propagated the $P_3$ uncertainty into all summary ages (reported with their $2\sigma$ interval), which is necessary when comparing TCN-based eruption ages to those from other geochronological methods (e.g. $ ^{40}Ar/$$^{39}Ar $). In the case of flows for which fewer than three samples passed the single population test (or only two samples were analysed), we considered the summary age to be a minimum eruption age. For those flows with three or more exposure ages passing this test, summary ages were considered robust eruption ages. We used internal $2\sigma$ intervals (INT $2\sigma$, not including $P_3$ errors) to compare intra-site and inter-site age distributions and clustering.

# 4 Results

## 4.1 Bulk rock and mineral geochemistry

Major and trace element concentrations of bulk rocks and minerals from each of the lava flows studied can be found in Table A2.

All bulk rock analyses yielded basaltic andesite to andesitic compositions according to the classification scheme of Le Maitre (2002). Our results indicate that, from the sampled flows, younger flows tend to be less evolved than older flows (Fig. 3a).

Most flows have a bulk geochemistry similar to the reported ranges (Conway et al., 2016) for the respective units they were classified as (Townsend et al., 2017). The only exception is the site here referred to as NR, which shows higher MgO (6.22 wt %) and lower Na₂O (2.95 wt %) than other samples of the Makotuku Member (2 wt %–3 wt % and 3.4 wt %–4 wt %, respectively; Conway et al., 2016). Instead, major element geochemistry of our NR sample matches that of the Mangaehuehu Member (4.7 wt %–7 wt % MgO and
<figure><img src="image_8.png" /><figcaption>**Figure 3.** (a) TAS classification diagram of the sampled lava flows (Le Maitre, 2002). Coloured areas represent geochemical ranges of Whakapapa and Mangawhero lavas. (b) Pyroxene compositions according to the classification scheme of Morimoto et al. (1988). Each triangle represents the average geochemistry of each lava flow's pyroxene population. (c) U, Th, and Sm concentrations in the samples. The x and y axes represent the maximum concentrations in minerals (pyroxenes and olivines) and in bulk rock, respectively.</figcaption></figure>

3 wt %–3.4 wt % Na₂O; Conway et al., 2016), the lavas of which are significantly older (Table 1).

Mineral geochemistry shows that, on average, the pyroxenes are pigeonite (Fig. 3b), although analyses of modal phases of Ruapehu lavas (Hackett, 1985; Conway, 2016) suggest that this represents a combination of augite and enstatite

crystals. MN and TSa yield average compositions of enstatite phases, indicating that the orthopyroxene phase dominates over the clinopyroxene in these flows. The analysed olivines (sample GR-PD023) are magnesium-rich (Fo₆₉; Table A2). Comparing the obtained average compositions with previous ³Hecos studies, our pyroxenes show higher contents of or-
thopyroxene than those analysed by Blard et al. (2006) and higher clinopyroxene contents than samples of Eaves et al. (2015).

In general, trace element concentrations are relatively homogeneous across the sampled sites. Figure 3c shows the concentrations of the main radioactive elements producing $^4$He$_{rad}$ (U, Th, and Sm) in bulk rock and in the mineral phases (pyroxenes and olivines). Bulk rocks contain 0.94–1.74 ppm U, 4.04–6.50 ppm Th, and 2.41–3.25 ppm Sm, while pyroxenes contain 0.01–0.10 ppm U, 0.04–0.36 ppm Th, and 0.44–2.07 ppm Sm (uncertainties < 20% and detection limits of 0.01 ppm). Note that U and Th concentrations in the rock are not involved in the production of the measured $^4$He$_{rad}$, as the external crystal rims were removed before the analyses. GR olivines have lower contents of these elements (with U below the detection limit) and, therefore, larger $P_4$ associated errors. However, element concentrations provided for minerals represent maximum values, as there is a possibility of groundmass and/or melt inclusion contamination that may not be accounted for at the time of measurement. These values indicate (maximum) partition coefficients (Kd) of 0.006–0.085 for U, 0.006–0.080 for Th, and 0.15–0.74 for Sm in pyroxenes and 0.045 for U, 0.045 for Th, and 0.11 for Sm in olivines. The pyroxene maximum Kd values, in general, agree with values from the literature (Dostal et al., 1983; Luhr and Carmichael, 1980; Gallahan and Nielsen, 1992; Nicholls and Harris, 1980). Those values for U and Th in olivines are similar to those reported by Dunn and Sen (1994) and Villemant (1988), while the Kd for Sm in our olivines is an order of magnitude larger than that of Dunn and Sen (1994), which can be explained by the impact of fluid inclusions with higher Sm contents within the olivine crystals.

## 4.2 Local magmatic $^3$He/$^4$He ratio

We measured $^3$He and $^4$He released after in vacuo crushing for samples MA-PD058, WG-PD326, and DC-PD329 (data available in Supplement file S2.1). These values result in $^3$He/$^4$He ratios of 5.5 ± 1.0, 17.9 ± 6.9, and 9.2 ± 6.1 × 10⁻⁶ for each sample, respectively. The large uncertainties associated with the ratios measured in samples WG-PD326 and DC-PD329 are a result of the low total $^4$He values (<5 × 10⁹ at g⁻¹).

We used the three measured $^3$He/$^4$He values and two ratios from Patterson et al. (1994) to constrain the local magmatic $^3$He/$^4$He value. These are 6.5 ± 2.4 × 10⁻⁶ (Ohakune basalt pyroxenes, one sample) and 8.6 ± 3.7 × 10⁻⁶ (Waimarino basalt olivines, mean of three aliquots from one sample), which are comparable to those obtained by in vacuo crushing of our samples. All analyses from Patterson et al. (1994) are from fused samples (and not from in vacuo crushed samples, which is the standard approach to release predominantly magmatic He; Kurz, 1986b), but we assumed that all the measured $^3$He and $^4$He have a magmatic origin, as the samples come from flow interiors

of young flows (i.e. they likely contain low $^4$He$_{rad}$ and minimal to no $^3$He$_{cos}$). With these data, we calculated an uncertainty-weighted mean $^3$He/$^4$He ratio using IsoplotR (Vermeesch, 2018) and obtained a value of 5.9 ± 2.6 × 10⁻⁶ (or 4.2 ± 1.9 Ra; see Supplement file S2.2), which we used for the magmatic corrections in this study. The impact of the obtained magmatic ratio and its uncertainty in our results is described in Sect. 4.3.

## 4.3 Fusion-released helium isotopes and cosmogenic $^3$He concentrations

We analysed a total of 77 samples from 23 individual flows. All fusion $^3$He and $^4$He measurements, calculated $^3$He$_{cos}$ concentrations, and derived exposure and eruption ages are shown in Table 2. Measured $^3$He varies between 2.1 × 10⁶ and 2.4 × 10⁷ at g⁻¹, with 2%–7% of relative associated error (1σ). $^4$He$_{tot}$ values are surprisingly low across most of our samples (possibly due to the repeated HF-leaching steps the samples were exposed to prior to analysis; see Bromley et al., 2014), typically ranging between 0.3 and 9.6 × 10¹⁰ at g⁻¹ with uncertainties between 0.04 and 0.18 × 10¹⁰ at g⁻¹. These values normally result in total $^3$He/$^4$He ratios of 130–800 Ra, although they are lower (50–90 Ra) or higher (1200–1500 Ra) in some cases (see Table 2).

The complete detail of all sources of corrections is available in Table A3 and in Supplement file S4. Calculated $^3$He$_{nuc}$ production rates ($P_{nuc}$) are 4 orders of magnitude below $P_3$ values, making $^3$He$_{cos}$ results insensitive to nucleogenic corrections. $P_4$ ranges between 4 × 10⁴ and 3 × 10⁵ at g⁻¹ yr⁻¹. We assume a 10% error associated with all $P_4$ results except for the site GR, which has lower concentrations of radioactive elements (and, hence, the lowest $P_4$ number within our lavas) with uncertainties of 20%–40%, for which we considered a 25% uncertainty in our $P_4$ estimates. Uncertainties associated with the calculated local magmatic $^3$He/$^4$He ratio represent < 10% of the informed error associated with $^3$He$_{cos}$ results. This ratio, combined with our $P_4$ calculations (3.5 × 10⁴–6.1 × 10⁵ at g⁻¹ yr⁻¹) and local $P_3$ values between 313 and 584 at g⁻¹ yr⁻¹ (el evations between 1288 and 2148 m a.s.l.), yields R factors > 0.99 (Table A3). This indicates that, even if the measured concentrations of radioactive elements in our minerals represent maximum values, corrections for $^4$He$_{rad}$ have a minor (< 1%) impact on our final $^3$He$_{cos}$ values. The $^3$He$_{cos}$/^3$He$_{tot}$ ratios calculated for our samples vary between 0.90 and 0.99, implying that the $^3$He$_{cos}$ quota dominates over magmatic (and nucleogenic) $^3$He.

The used magmatic $^3$He/$^4$He ratio (5.9 ± 2.6 × 10⁻⁶ or 4.2 ± 1.9 Ra) is derived from only three samples from Ruapehu crushed in vacuo (this study) and two samples from the Ohakune and Waimarino basalts (data from fused samples; Patterson et al., 1994), and it is therefore not well constrained. However, this does not significantly affect our final
results due to the low $^4$He$_{tot}$ measured in most of our samples. To demonstrate this, we estimated the resulting $^3$He$_{cos}$ concentration if the magmatic $^3$He/$^4$He ratio of the sample SC-PD001 (which has the smallest measured $^3$He/$^4$He ratio across our samples and is, hence, the most sensitive to this test) was 8.4 and 2.1 Ra (twice and half the mean value of 4.2 Ra used for our calculations and covering most of the range globally observed in subduction zone volcanism; Hilton et al., 2002). This test yields $^3$He$_{cos}$ concentrations of 3.28 ± 0.31 and 3.71 ± 0.31 × 10$^6$ at g$^{-1}$ (resulting in exposure ages of 9.74 ± 0.85 and 10.91 ± 0.83 ka) with magmatic ratios of 8.4 and 2.1 Ra, respectively, both falling within the error of the concentration obtained using a magmatic ratio of 4.2 ± 1.9 Ra for SC-PD001 (3.57 ± 0.31 × 10$^6$ at g$^{-1}$ and an exposure age of 10.53 ± 0.83 ka). This indicates that the potentially variable magmatic $^3$He/$^4$He ratios present in our samples do not significantly impact our results, although they might partially explain small differences between the obtained exposure ages of samples from the same flow.

## 4.4 Lava flows: background and new cosmogenic $^3$He constraints

We obtained 16 eruption ages and 7 minimum eruption ages (Table 2) based on the criteria defined in Sect. 3.6.2.

### 4.4.1 Iwikau Member (Whakapapa Formation)

The Iwikau Member of the Whakapapa Formation covers a large area on the northwestern and eastern flanks of Ruapehu (Fig. 1) and is subdivided into three flow packages: Tawhainui, Mangatoetoenui, and Taranaki Falls flows (Fig. 4a, b), all interpreted to have originated from Ruapehu's northern vent (Townsend et al., 2017).

#### Tawhainui flows

The Tawhainui flows comprise a voluminous sequence of lava flows on the northwestern slopes of the volcano. They have been the most studied unit of Ruapehu due to their accessibility and availability of the fresh exposures of flow interiors, facilitated by the construction of the largest ski field in Te Ika-a-Māui/North Island. We sampled three flows from this unit: the Delta Corner flow (DC samples), the Bruce Road flow (BR samples, both after Greve et al., 2016), and the Whakapapa Glacier flow (WG samples).

The fresh-looking Delta Corner flow was previously dated at 6.0 ± 2.4 ka with $^{40}$Ar/$^{39}$Ar by Conway et al. (2016), and this age was refined to 8200–7900 BP by Greve et al. (2016) based on paleomagnetic data. Analyses from three samples from an area with distinct 'a'ā surface morphologies (see Fig. 2b) yield well-clustered exposure ages, which result in an eruption age of 7.8 ± 1.5 ka. Our results are consistent with the age range of 8200–7900 BP provided by Greve et al.

(2016), suggesting that the flow's true age lies on the upper end of the uncertainty provided by Conway et al. (2016).

The Bruce Road flow is a large 'a'ā flow that underlies the Delta Corner flow and has been constrained to 8800–8500 BP by Greve et al. (2016) using paleomagnetism. Downslope from the BR sample site, the flow has unclear boundaries, as it is covered by vegetation. Based on four individual exposure ages (7.4, 8.1, 7.8, and 9.1 ka), we obtained an eruption age of 8.1 ± 2.1 ka for the Bruce Road flow, which is consistent with its paleomagnetic constraint.

The Whakapapa Glacier flow is one of the youngest lavas of the sequence based on stratigraphic relations, which suggest a comparable age to that of the Delta Corner flow. Due to the highly eroded nature of the Whakapapa Glacier flow's surface, only two WG samples were collected, which yield a minimum eruption age of 7.8 ± 2.4 ka. This result is consistent with the stratigraphy and the age of the Delta Corner flow.

#### Mangatoetoenui flows

This subunit includes a group of lava flows on the eastern slopes of Ruapehu, and its age is poorly constrained (Table 1). We sampled four individual flows classified based on geochemistry and location within the Mangatoetoenui flows: Lava Cascade (LC samples), Tukino Slopes-a (TSa samples), Tukino Slopes-b (TSb samples), and Tukino Flats (TFt samples) flows (Fig. 4b).

The LC sample site is interpreted to be part of an approximately 4 km long lava flow terminating at a 20 m high lava cascade at 1620 m a.s.l. This flow was described in detail by Rhodes (2012) and dated on a cliff at its terminus at 0.8 ± 5.6 ka by Conway et al. (2016). We analysed four individual samples from an outcrop located ~1 km upslope from the lava toe and obtained an eruption age of 11.4 ± 2.3 ka for the Lava Cascade flow (outside the 2σ interval of Conway et al., 2016; see Sect. 5.1), with one young outlier removed (sample LC-PD256). The outlier can be explained by local erosion, shielding from a now collapsed neighbouring lava tumuli (and hence an underestimation of the shielding factor) or a period of tephra cover that could have reduced the $^3$He production on the surface sampled.

The Tukino Slopes-a flow has not been previously dated, but its location and stratigraphic position suggest a similar eruption age to the Lava Cascade and Tukino Slopes-b flows. All measured TSa samples (8.7, 9.5, and 9.9 ka) form a single population and provide an eruption age of 9.4 ± 1.8 ka, in good agreement with the stratigraphy.

The TSb sample site likely corresponds to the same flow dated with $^{40}$Ar/$^{39}$Ar at 9.2 ± 8.0 ka by Conway et al. (2016). We obtained exposure ages of 10.5, 11.9, and 11.9 ka for the TSb samples, which result in a refined eruption age of 11.5 ± 2.2 ka for the Tukino Slopes-b flow. The eruption age we obtained for the Tukino Slopes-b flow is consistent with (and more precise than) the existing radiometric age and, as
<figure><img src="image_9.png" /></figure>

**Figure 4.** Map of dated < 21 ka lava flows on (a) northern, (b) eastern, (c) western, and (d) southern Ruapehu. Polygons redefined from Townsend et al. (2017). Boundaries of the Mangaehuehu Member (Mangawhero Formation), as of Townsend et al. (2017), shown for context of site NR of this study. The shaded-grey areas represent postglacial flows without chronological data.
suggested by the stratigraphy, similar to the age of the Tukino Slopes-a flow.

The TFt sample site is located at a lower elevation (~1515 m a.s.l.), and its stratigraphic position suggests a similar or older age than the rest of the Mangatoetoenui flows. Our results for three TFt samples (7.7, 10.8, and 7.3 ka) do not form a single population and result in a minimum eruption age of 8.6 ± 4.6 ka for the Tukino Flats flow. The older exposure age (10.8 ka) is difficult to explain as an outlier, as the presence of inherited ³He is not justifiable for lava flows, whereas the younger ages may be explained as outliers owing to surface erosion or temporal burial by alluvium or tephra. Lack of additional samples hindered our ability to obtain a robust eruption age. Considering a minimum eruption age of 8.6 ± 4.6 ka, the ages of the other flows from the Mangatoetoenui flows, and their stratigraphic position, our best estimate for the Tukino Flats flow is 12–10 ka.

### Taranaki Falls flow

The Taranaki Falls flow (TFa samples) is a rootless (not continuous towards the vent it would have been erupted from), elongated lava flow that outcrops discontinuously for ~8 km almost directly to the north of the volcano's summit area (Townsend et al., 2017, Fig. 4a) and terminates at the 20 m waterfall after which it is named. The flow was dated with $^{40}$Ar/$^{39}$Ar at 8.8 ± 2.8 ka by Conway et al. (2016). Based on this date, Greve et al. (2016) found two age ranges (10 800–10 200 and 9500–8900 BP) with a better match to the local paleosecular variation record.

We sampled the flow at an outcrop 800 m upstream from the flow terminus and obtained exposure ages of 14.6, 14.2, and 15.0 ka, resulting in an eruption age of 14.6 ± 2.9 ka, which is outside the confidence interval of the radiometric age (see Sect. 5.1).

### 4.4.2 Saddle Cone Member (Whakapapa Formation)

This unit comprises a large, lobate 'a'ā flow originating from a parasitic cone on the north-northeastern side of Ruapehu that is almost disconnected from the main edifice (Fig. 4a). The only available constraint for this flow was provided by Greve et al. (2016), who suggested an age of 9850–8650 BP based on paleomagnetic analyses of samples from the western lobe of the flow. The Saddle Cone Member also includes a smaller blocky lava flow lying between this cone and Ruapehu's summit region (that likely originated from a satellite vent), adjacent to the Waihohonu Ridge and here referred to as the Waihohonu Plateau flow, which was linked to the main Saddle Cone deposits by its geochemical similarity and location. Nairn et al. (1998) suggested that the Waihohonu Plateau ("1990 m lava" therein) might be younger than 5 ka as no deposits from the Papakai tephra were found above the flow.

Individual exposure ages of samples from the main western lobe of the Saddle Cone lavas (SCw samples; 10.5, 10.2, and 9.2 ka, Fig. 2b) show good agreement. Additionally, we analysed a sample from the eastern lobe (SCe, whose surface elevation is more than 100 m below that of the main lobe; see Fig. 4a) to test the hypothesis of a multi-episodic origin. The obtained exposure age of this sample is 9.6 ka (Table 2), indistinguishable from those of the western lobe. We suggest a single eruption age of 9.9 ± 2.0 ka ($n = 4$) for both lobes, which is consistent with the existing paleomagnetic constraint for this flow.

The blocky nature of the Waihohonu Plateau flow made it difficult to find uneroded surfaces, and only two samples were obtained (WP samples, Fig. 2f). Analyses from these samples result in a minimum eruption age of 11.2 ± 2.2 ka.

### 4.4.3 Pinnacle Ridge Member (Whakapapa Formation)

The Pinnacle Ridge Member is a welded spatter deposit (Fig. 2e) linked to a dike on a ridge of the same name on the northern flanks of the volcano (Fig. 4a). Due to its geochemistry and geographic location, Donoghue et al. (1999) linked this isolated spatter-fed lava deposit to the Taurewa pyroclastic unit (ca. 10 ka) described by Topping and Kohn (1973), manifested as a tephra layer with isopachs centred on the northern flanks of Ruapehu.

PR samples yield exposure ages of 20.8, 19.0, and 21.5 ka, resulting in an eruption age of 20.4 ± 4.0 ka for this unit, which is at least 10 kyr prior to the Taurewa eruptive event.

### 4.4.4 Rangataua Member (Whakapapa Formation)

The Rangataua Member includes the longest and most voluminous known lava flow of Ruapehu ($\ge$ 15 km long and ~1.5 km³). It first outcrops ~3.5 km south from the summit, which led to the hypothesis that it is sourced from a satellite vent (Hackett, 1985; Price et al., 2012), although Townsend et al. (2017) suggest initial transport over ice as a possible alternative explanation for its rootless nature. Based on geochemical differences, this unit was first subdivided by Price et al. (2012) into proximal, medial, and distal flows (the latter being the largest flow of the sequence). They suggested eruption ages of 12–10 ka based on underlying and overlying tephra sequences (unpublished data). These lavas overlies left lateral moraines at 1600–1400 m a.s.l., which have been correlated to right lateral moraines of the Mangaehuehu River valley dated at 11–14 ka (Eaves et al., 2019) using ³He/cos dating (Fig. 4d). We sampled the Rangataua Member at two locations: one close to the highest outcrops (RTp, "proximal") and another one approximately 1 km to the south (RTm, "medial"). We did not sample the distal flows, which are interpreted to be older than the medial flows, due to vegetation cover (Fig. 1c).

RTp samples (e.g. Fig. 2a) yield exposure ages of 13.9, 12.4, 13.9, and 14.4 ka and a final eruption age of
13.6 ± 2.6 ka. Results of RTm samples (16.2, 16.0, 15.3, and 8.2 ka) include a young outlier, but the remaining samples are internally consistent and indicate an eruption age of 15.8 ± 3.0 ka, which agrees with the field relationships of the area, as this flow underlies RTp, but not so with previous age estimates (see Sect. 5.1). The ages of the Rangataua proximal and medial flows and their INT 2σ uncertainties (13.6 ± 0.6 and 15.8 ± 0.8, respectively) do not overlap, indicating that they correspond to different eruptive episodes.

### 4.4.5 Paretetaitonga Member (Whakapapa Formation)

The Paretetaitonga Member comprises a series of lava flows that likely originated from the northern summit vent of Ruapehu and emplaced in the headwaters of the Whakapapaiti Stream, northwest of the summit area (Fig. 4a). We sampled one lava flow (Whakapapaiti flow, WT samples) stratigraphically higher than the only flow dated from this unit (14.8 ± 3.0 ka; Conway et al., 2016).

We obtained exposure ages in good agreement with each other (12.8, 13.4, and 13.7 ka), resulting in an eruption age of 13.3 ± 2.6 ka, consistent with the existing chronology.

### 4.4.6 Turoa Member (Whakapapa Formation)

The Turoa Member corresponds to a sequence of numerous flows extending directly west from the edge of Ruapehu's crater rim and reaching the Mangaturuturu valley bottom. Based on the distributions of the flows and two ⁴⁰Ar/³⁹Ar dates (Table 1), this unit is assumed to have been formed by effusive activity from the southern summit vent at ca. 17–10 ka. We sampled five sites, distributed on the northern (MN and MS), central (CTa and CTb), and western (TC) areas (Fig. 4c), covered by this unit.

The Mangaturuturu North flow corresponds to a flow on the headwaters of the Mangaturuturu Stream, and, due to stratigraphic relations and flow morphologies, it was suspected to be the youngest lava on western Ruapehu. We analysed five surfaces of the Mangaturuturu North flow (MN samples, Fig. 2c, with exposure ages of 8.0, 8.9, 6.0, 8.9, and 7.7 ka), and, eliminating the young outlier of 6.0 ka, they yield a robust eruption age of 8.3 ± 1.6 ka.

The Mangaturuturu South flow underlies the Mangaturuturu North flow and extends down ~3 km from the summit area. Poor exposures of original flow surfaces prevented us from collecting more than three samples from the Mangaturuturu South flow (MS samples). Additionally, purification of the minerals in these samples was incomplete due to high (> 50%) mass loss with each HF-leaching cycle, and we suspect an overestimation of measured pyroxene mass for these samples. Sample analyses result in exposure ages that do not pass the single population test (Table 2) but provide a minimum eruption age of 6.1 ± 1.7 ka.

The Central Turoa-a and Central Turoa-b flows are located close to each other and at a similar elevation south of

the MN and MS sample sites. We only collected two samples from the Central Turoa-a flow (CTa) due to a lack of suitable surfaces, which suggest a minimum eruption age of 13.6 ± 2.7 ka. Three out of the four analysed Central Turoa-b flow samples (CTb samples, which yielded exposure ages of 4.9, 8.8, 8.4, and 8.5 ka) show good agreement and yield an eruption age of 8.6 ± 1.7 ka. These results indicate that the Central Turoa-a and Central Turoa-b flows correspond to two different eruptive episodes.

The Turoa Cascades flow (TC samples) is a large flow that reaches the Mangaturuturu River valley floor, and its stratigraphic position indicates that it is likely the oldest flow of the Turoa Member. Individual exposure ages of the TC samples (11.4, 14.1, 13.1, and 13.3 ka) include a young outlier and indicate an eruption age of 13.4 ± 2.6 ka for the Turoa Cascades flow, which is in good agreement with the rest of the ages obtained for the Turoa Member lavas.

### 4.4.7 Makotuku Member (Mangawhero Formation)

We sampled three flows previously mapped as part of the Makotuku Member of the Mangawhero Formation: the Makotuku Flat flow (MF samples) on the southwest and the Ngā Rimutāmaka and Makahikatoa flows (NR and MA samples, named after local site and stream, respectively) on the south of Ruapehu's edifice. The spatial distribution of Makotuku lavas suggests that they originated from the southern summit vent.

The Makotuku Flats flow extends to the west of the edifice, reaching the Makotuku valley bottom (Fig. 4d), and overlies a 11–15 ka moraine (Townsend et al., 2017) at the sampled site. Although results of analyses of MF samples are not particularly well clustered, they behave as a single population and provide an eruption age of 12.6 ± 3.5 ka.

Analyses of NR samples yield well-clustered exposure ages, and we interpret an eruption age of 42.9 ± 8.6 ka, which corresponds to the only date provided for this lava flow so far. It is worth noting that this age and the geochemical composition of this flow match with the ⁴⁰Ar/³⁹Ar ages and the high-MgO and low-Al₂O₃ nature of the Mangaehuehu Member lavas (Table 1).

The small area where the Makahikatoa flow outcrops prevented us from obtaining more than two suitable samples, which result in a minimum eruption age of 54.0 ± 18.0 ka, this being the first age constraint for this flow.

These three eruption ages do not contradict previous chronology or stratigraphy, but they do not match the age ranges indicated by the geochemical affinities for the Makotuku Member lavas as described by Conway et al. (2016) and Townsend et al. (2017; see Sect. 5.2).

### 4.4.8 Mangaehuehu Member (Mangawhero Formation)

We sampled a lava flow (Girdlestone Ridge, GR samples) outcropping on a ridge top ~1.5 km south from Ruapehu's
summit and 800 m southwest from Girdlestone Peak. This site was previously mapped as Mangaehuehu Member lavas (Townsend et al., 2017) based on interpretation of aerial im- agery. However, the uneroded aspect of the flow's surface observed in the field during this study suggests that it could be younger than previously interpreted. The mineral separa- tion process applied to all samples produced the only olivine concentrate (with a minor pyroxene population) of this study.

Analyses of these samples indicate a minimum eruption age of $14.2 \pm 2.7$ ka (mean calculated from the two oldest exposure ages after the elimination of two outliers), which represents the first age constraint for this lava flow.

# 5 Discussion

## 5.1 Comparison of new cosmogenic $^3$He ages with previous age constraints

The new Holocene $^3$He exposure ages yield eruption ages with higher precision than $^{40}$Ar/$^{39}$Ar dates of Conway et al. (2016) for this time range (Fig. 5). Additionally, young ($< 20$ ka) $^{40}$Ar/$^{39}$Ar ages of individual samples have nor- mally weak isochrons, as the R values for their linear fits used to calculate crystallization age (released $^{40}$Ar/$^{36}$Ar vs. $^{39}$Ar/$^{36}$Ar in increasing temperature steps) tend to be rela- tively low (e.g. Harpel et al., 2004; Conway et al., 2016; Preece et al., 2018). Therefore, young $^{40}$Ar/$^{39}$Ar ages are very susceptible to the decisions involved in the selection of steps included (or discarded) in the calculation of weighted mean plateau and isochron ages, and our exposure ages based on multiple samples provide more reliable results.

From the four sampled flows in this study with existing $^{40}$Ar/$^{39}$Ar dates (Conway et al., 2016), two yield eruption ages that agree with the radiometric dates (Delta Corner and Tukino Slopes-b flows) and two not only outside the $2\sigma$ confidence interval of Conway et al. (2016) but older than the $^{40}$Ar/$^{39}$Ar ages, the Lava Cascade ($^3$He$_{cos}$: $11.4 \pm 2.3$ ka and $^{40}$Ar/$^{39}$Ar: $0.8 \pm 5.6$ ka, Mangatoetoenui flows), and the Taranaki Falls ($^3$He$_{cos}$: $14.6 \pm 2.9$ ka and $^{40}$Ar/$^{39}$Ar: $8.8 \pm 2.8$ ka) flows. The imprecise nature of the radiomet- ric age of the Lava Cascade flow and its weak isochron, to- gether with the good agreement between our LC samples and the eruption ages we obtained for the Mangatoetoenui flows, leads us to conclude that our eruption age for the Lava Cas- cade flow is more robust than the date provided by Conway et al. (2016). Based on the good clustering of our results (Ta- ble 2), we suggest that our $^3$He$_{cos}$ eruption age better repre- sents the true eruption age of the Taranaki Falls flow. Addi- tionally, our eruption age would explain the rootless nature of the flow (Townsend et al., 2017), as it is older than the flank collapse event that affected the northern summit area of Ruapehu at ca. 10.5 ka (Eaves et al., 2015) and, hence, the upper section of the Taranaki Falls flow (Fig. 6a).

Our results show, in general, good agreement with the lava flow eruption ages refined by Greve et al. (2016) at Ruapehu

(Fig. 5). The only exception is the Taranaki Falls flow; the refinement by Greve et al. (2016) is based on the $^{40}$Ar/$^{39}$Ar date of Conway et al. (2016), and, thus, it intrinsically agrees with this age and not with our results. Our $^3$He$_{cos}$ erup- tion ages for the Delta Corner ($7.8 \pm 1.5$ ka; INT $2\sigma$ 0.6 ka), Bruce Road ($8.1 \pm 2.1$ ka; INT $2\sigma$ 1.5 ka), and Saddle Cone ($9.9 \pm 2.0$ ka; INT $2\sigma$ 0.7 ka) flows match the respective age ranges of 8200–7900, 8800–8500, and 9850–8650 BP pro- vided by Greve et al. (2016). Moreover, these results suggest that it is unlikely that $P_3$ errors have a significant impact on the accuracy of the eruption ages from this work, which is also supported by the good agreement of the local $^3$He$_{cos}$ production rate calibration test by Eaves et al. (2015) with the worldwide mean production rate used in this study.

Eruption ages obtained for the Rangataua proximal and medial flows ($13.6 \pm 2.6$ and $15.8 \pm 3.0$ ka, respectively) do not agree with a $12–10$ ka constraint suggested by Price et al. (2012) based on tephra stratigraphy (using unpublished data). However, tephra correlation on Ruapehu is complex due to the large number of pyroclastic units emplaced at $20–11$ ka and their broad geochemical ranges (Pardo et al., 2012a). Detailed studies (Donoghue et al., 2007) attempted to sys- tematize tephra correlation in this area without success, in- dicating that the andesitic tephras are highly heterogeneous, displaying wide compositional fluctuations during short time intervals. Hence, our eruption ages are more robust than the estimate of $12–10$ ka by Price et al. (2012). The other existing constraint for the Rangataua flows was given by a right lat- eral moraine of the Mangaehuehu valley dated at $11–14$ ka by Eaves et al. (2019), which was thought to correspond in age to the left lateral moraine overlain by the RTm flow (Fig. 4d). Our eruption age of $15.8 \pm 0.8$ ka (INT $2\sigma$; $P_3$ errors not con- sidered as the moraines were dated using $^3$He$_{cos}$) suggests that the moraine underlying the Rangataua flows is older than the dated right lateral moraine, rather than its equivalent.

Most of the flows dated in this study lack previous age con- straints beyond estimations based on geochemical similarity and geographical proximity to lavas with $^{40}$Ar/$^{39}$Ar dates. The eruption ages obtained for about half of these flows do not agree with these correlations (Fig. 5). Five of them (MN, MS, CTb, MF, and GR flows) yield ages younger than any of the dates informed for the units they were correlated to (i.e. Turoa, Makotuku, and Mangaehuehu members). This can be explained by a sampling bias of Conway et al. (2016) towards older flows that are more likely to have exposed their slowly cooled flow interiors (suitable for $^{40}$Ar/$^{39}$Ar dating) due to their longer periods exposed to erosive processes and the presence of collapsed thick margins in the case of pre- viously ice-impounded flows (Conway et al., 2015). PR and MA deposits are relatively isolated (Fig. 4a, b), so the pre- vious geochemical correlations are weaker. The age previ- ously assigned to the PR deposits (Table 1) was, unlike any other lava in this study, based on a correlation with a pyro- clastic unit, adding another layer of uncertainty. Our results represent the first dates for lavas at the PR and MA sites and
<table><caption>Table 2. Continued.</caption><thead><tr><th>Sample</th><th>Latitude (S)</th><th>Longitude (E)</th><th>Elevation (m.s.l.)</th><th>Shielding factor</th><th><sup>3</sup>He<sub>tot</sub> ±1σ (10<sup>6</sup> at g<sup>-1</sup>)</th><th><sup>4</sup>He<sub>tot</sub> ±1σ (10<sup>10</sup> at g<sup>-1</sup>)</th><th>Total <sup>3</sup>He/<sup>4</sup>He (R/Ra)</th><th><sup>3</sup>He<sub>cos</sub> ±1σ (10<sup>6</sup> at g<sup>-1</sup>)</th><th>Exposure age ±1σ (ka)</th></tr></thead><tbody><tr><td colspan="10"><strong>PR – Pinnacle Ridge spatter deposit</strong></td></tr><tr><td>PR-PD083</td><td>39.2370</td><td>175.5672</td><td>1730.7</td><td>0.979</td><td>9.39 ± 0.44</td><td>6.56 ± 0.20</td><td>103 ± 6</td><td>9.03 ± 0.48</td><td>20.82 ± 1.00</td></tr><tr><td>PR-PD084</td><td>39.2386</td><td>175.5689</td><td>1860.9</td><td>0.988</td><td>9.42 ± 0.44</td><td>7.72 ± 0.24</td><td>88 ± 5</td><td>8.99 ± 0.49</td><td>19.04 ± 0.93</td></tr><tr><td>PR-PD085</td><td>39.2385</td><td>175.5688</td><td>1857.9</td><td>0.997</td><td>10.69 ± 0.49</td><td>5.84 ± 0.18</td><td>132 ± 7</td><td>10.37 ± 0.52</td><td>21.48 ± 1.00</td></tr><tr><td colspan="10"><strong>RTp – Rangataua proximal flow</strong></td></tr><tr><td>RTp-PD027</td><td>39.3140</td><td>175.5509</td><td>1831.4</td><td>0.997</td><td>6.38 ± 0.25</td><td>1.74 ± 0.08</td><td>264 ± 16</td><td>6.28 ± 0.26</td><td>13.88 ± 0.53</td></tr><tr><td>RTp-PD028</td><td>39.3140</td><td>175.5509</td><td>1833.1</td><td>0.996</td><td>5.73 ± 0.26</td><td>2.13 ± 0.12</td><td>194 ± 14</td><td>5.61 ± 0.26</td><td>12.41 ± 0.56</td></tr><tr><td>RTp-PD029</td><td>39.3140</td><td>175.5509</td><td>1832.9</td><td>0.996</td><td>6.42 ± 0.36</td><td>2.24 ± 0.18</td><td>206 ± 20</td><td>6.29 ± 0.37</td><td>13.88 ± 0.75</td></tr><tr><td>RTp-PD030</td><td>39.3143</td><td>175.5512</td><td>1816.4</td><td>0.988</td><td>6.45 ± 0.31</td><td>0.98 ± 0.06</td><td>474 ± 35</td><td>6.40 ± 0.31</td><td>14.35 ± 0.62</td></tr><tr><td colspan="10"><strong>RTm – Rangataua medial flow</strong></td></tr><tr><td>RTm-PD045</td><td>39.3234</td><td>175.5520</td><td>1585.9</td><td>0.991</td><td>6.30 ± 0.38</td><td>1.83 ± 0.07</td><td>248 ± 17</td><td>6.20 ± 0.38</td><td>16.22 ± 0.91</td></tr><tr><td>RTm-PD046 (a)</td><td>39.3249</td><td>175.5508</td><td>1567.6</td><td>0.979</td><td>6.14 ± 0.25</td><td>0.74 ± 0.05</td><td>601 ± 48</td><td>6.10 ± 0.25</td><td></td></tr><tr><td>RTm-PD046 (b)</td><td>39.3249</td><td>175.5508</td><td>1567.6</td><td>0.979</td><td>5.98 ± 0.30</td><td>0.98 ± 0.08</td><td>438 ± 43</td><td>5.93 ± 0.30</td><td></td></tr><tr><td>RTm-PD046 mean</td><td>39.3249</td><td>175.5508</td><td>1567.6</td><td>0.979</td><td></td><td></td><td></td><td>6.04 ± 0.27</td><td>16.03 ± 0.64</td></tr><tr><td>RTm-PD047</td><td>39.3251</td><td>175.5503</td><td>1567.4</td><td>0.997</td><td>5.91 ± 0.29</td><td>1.13 ± 0.06</td><td>376 ± 28</td><td>5.85 ± 0.29</td><td>15.27 ± 0.68</td></tr><tr><td>RTm-PD048</td><td>39.3250</td><td>175.5503</td><td>1567.3</td><td>0.997</td><td>3.08 ± 0.16</td><td>1.64 ± 0.07</td><td>135 ± 9</td><td>2.99 ± 0.16</td><td>8.17 ± 0.40*</td></tr><tr><td colspan="10"><strong>WT – Whakapapaiti flow</strong></td></tr><tr><td>WT-PD073</td><td>39.2569</td><td>175.5428</td><td>1892.4</td><td>0.987</td><td>6.01 ± 0.28</td><td>0.35 ± 0.08</td><td>1230 ± 270</td><td>6.00 ± 0.28</td><td>12.78 ± 0.59</td></tr><tr><td>WT-PD074</td><td>39.2569</td><td>175.5428</td><td>1892.5</td><td>0.991</td><td>6.36 ± 0.26</td><td>0.73 ± 0.04</td><td>624 ± 39</td><td>6.33 ± 0.27</td><td>13.41 ± 0.54</td></tr><tr><td>WT-PD075</td><td>39.2560</td><td>175.5397</td><td>1785.0</td><td>0.990</td><td>6.06 ± 0.26</td><td>1.24 ± 0.05</td><td>352 ± 20</td><td>6.00 ± 0.27</td><td>13.74 ± 0.57</td></tr><tr><td colspan="10"><strong>MN – Mangaturuturu North flow</strong></td></tr><tr><td>MN-PD217</td><td>39.2829</td><td>175.5322</td><td>1815.9</td><td>0.993</td><td>3.49 ± 0.24</td><td>0.30 ± 0.05</td><td>830 ± 150</td><td>3.48 ± 0.24</td><td>7.99 ± 0.50</td></tr><tr><td>MN-PD218</td><td>39.2829</td><td>175.5321</td><td>1813.9</td><td>0.993</td><td>3.82 ± 0.23</td><td>0.37 ± 0.04</td><td>739 ± 97</td><td>3.80 ± 0.23</td><td>8.85 ± 0.51</td></tr><tr><td>MN-PD219</td><td>39.2829</td><td>175.5321</td><td>1812.1</td><td>0.993</td><td>2.47 ± 0.19</td><td>0.42 ± 0.05</td><td>770 ± 200</td><td>2.46 ± 0.19</td><td>5.99 ± 0.39*</td></tr><tr><td>MN-PD220</td><td>39.2829</td><td>175.5322</td><td>1817.5</td><td>0.993</td><td>3.91 ± 0.25</td><td>0.77 ± 0.03</td><td>668 ± 89</td><td>3.89 ± 0.25</td><td>8.89 ± 0.54</td></tr><tr><td>MN-PD221</td><td>39.2829</td><td>175.5325</td><td>1822.8</td><td>0.993</td><td>3.30 ± 0.20</td><td>0.08 ± 0.04</td><td>1340 ± 400</td><td>3.29 ± 0.20</td><td>7.66 ± 0.42</td></tr><tr><td colspan="10"><strong>MS – Mangaturuturu South flow</strong></td></tr><tr><td>MS-PD222</td><td>39.2845</td><td>175.5304</td><td>1750.6</td><td>0.954</td><td>2.58 ± 0.16</td><td>0.50 ± 0.04</td><td>371 ± 35</td><td>2.55 ± 0.16</td><td>6.67 ± 0.35</td></tr><tr><td>MS-PD223</td><td>39.2845</td><td>175.5305</td><td>1751.4</td><td>0.992</td><td>2.51 ± 0.17</td><td>0.12 ± 0.05</td><td>1560 ± 710</td><td>2.51 ± 0.17</td><td>6.32 ± 0.36</td></tr><tr><td>MS-PD224</td><td>39.2845</td><td>175.5305</td><td>1750.9</td><td>0.992</td><td>2.08 ± 0.16</td><td>0.40 ± 0.11</td><td>370 ± 100</td><td>2.06 ± 0.16</td><td>5.33 ± 0.37</td></tr><tr><td colspan="10"><strong>CTa – Central Turoa-a flow</strong></td></tr><tr><td>CTa-PD229</td><td>39.2958</td><td>175.5395</td><td>1924.0</td><td>0.996</td><td>6.57 ± 0.33</td><td>2.35 ± 0.09</td><td>201 ± 13</td><td>6.45 ± 0.34</td><td>13.24 ± 0.67</td></tr><tr><td>CTa-PD230</td><td>39.2959</td><td>175.5396</td><td>1925.1</td><td>0.996</td><td>6.93 ± 0.36</td><td>2.76 ± 0.11</td><td>181 ± 12</td><td>6.78 ± 0.37</td><td>13.89 ± 0.69</td></tr><tr><td colspan="10"><strong>CTb – Central Turoa-b flow</strong></td></tr><tr><td>CTb-PD231</td><td>39.2998</td><td>175.5392</td><td>1877.5</td><td>0.996</td><td>2.11 ± 0.14</td><td>0.66 ± 0.06</td><td>230 ± 25</td><td>2.07 ± 0.14</td><td>4.93 ± 0.30*</td></tr><tr><td>CTb-PD232</td><td>39.3001</td><td>175.5390</td><td>1873.2</td><td>0.991</td><td>4.00 ± 0.24</td><td>0.74 ± 0.05</td><td>390 ± 33</td><td>3.96 ± 0.24</td><td>8.79 ± 0.49</td></tr><tr><td>CTb-PD233</td><td>39.3001</td><td>175.5390</td><td>1872.0</td><td>0.994</td><td>3.86 ± 0.25</td><td>0.93 ± 0.06</td><td>298 ± 28</td><td>3.81 ± 0.25</td><td>8.39 ± 0.52</td></tr><tr><td>CTb-PD234 (a)</td><td>39.3003</td><td>175.5391</td><td>1873.4</td><td>0.996</td><td>3.80 ± 0.24</td><td>0.90 ± 0.07</td><td>283 ± 28</td><td>3.75 ± 0.24</td><td></td></tr><tr><td>CTb-PD234 (b)</td><td>39.3003</td><td>175.5391</td><td>1873.4</td><td>0.996</td><td>3.96 ± 0.27</td><td>0.60 ± 0.05</td><td>472 ± 50</td><td>3.93 ± 0.26</td><td></td></tr><tr><td>CTb-PD234 mean</td><td>39.3003</td><td>175.5391</td><td>1873.4</td><td>0.996</td><td></td><td></td><td></td><td>3.85 ± 0.25</td><td>8.47 ± 0.51</td></tr><tr><td colspan="10"><strong>TC – Turoa Cascades flow</strong></td></tr><tr><td>TC-PD066</td><td>39.3014</td><td>175.5193</td><td>1533.2</td><td>0.997</td><td>4.13 ± 0.20</td><td>1.03 ± 0.06</td><td>288 ± 23</td><td>4.07 ± 0.21</td><td>11.37 ± 0.53*</td></tr><tr><td>TC-PD067</td><td>39.3015</td><td>175.5192</td><td>1533.6</td><td>0.997</td><td>5.24 ± 0.27</td><td>1.14 ± 0.05</td><td>331 ± 21</td><td>5.18 ± 0.26</td><td>14.09 ± 0.65</td></tr><tr><td>TC-PD068</td><td>39.3015</td><td>175.5192</td><td>1533.1</td><td>0.997</td><td>4.74 ± 0.23</td><td>1.09 ± 0.04</td><td>313 ± 18</td><td>4.68 ± 0.23</td><td>13.05 ± 0.62</td></tr><tr><td>TC-PD070</td><td>39.3012</td><td>175.5193</td><td>1528.0</td><td>0.997</td><td>4.89 ± 0.23</td><td>0.82 ± 0.04</td><td>431 ± 28</td><td>4.85 ± 0.23</td><td>13.27 ± 0.61</td></tr><tr><td colspan="10"><strong>TC – Turoa Cascades flow</strong></td></tr><tr><td>TC-PD066</td><td>39.3014</td><td>175.5193</td><td>1533.2</td><td>0.997</td><td>4.13 ± 0.20</td><td>1.03 ± 0.06</td><td>288 ± 23</td><td>4.07 ± 0.21</td><td>11.37 ± 0.53*</td></tr><tr><td>TC-PD067</td><td>39.3015</td><td>175.5192</td><td>1533.6</td><td>0.997</td><td>5.24 ± 0.27</td><td>1.14 ± 0.05</td><td>331 ± 21</td><td>5.18 ± 0.26</td><td>14.09 ± 0.65</td></tr><tr><td>TC-PD068</td><td>39.3015</td><td>175.5192</td><td>1533.1</td><td>0.997</td><td>4.74 ± 0.23</td><td>1.09 ± 0.04</td><td>313 ± 18</td><td>4.68 ± 0.23</td><td>13.05 ± 0.62</td></tr><tr><td>TC-PD070</td><td>39.3012</td><td>175.5193</td><td>1528.0</td><td>0.997</td><td>4.89 ± 0.23</td><td>0.82 ± 0.04</td><td>431 ± 28</td><td>4.85 ± 0.23</td><td>13.27 ± 0.61</td></tr><tr><td colspan="10"><strong>TC – Turoa Cascades flow</strong></td></tr><tr><td>TC-PD066</td><td>39.3014</td><td>175.5193</td><td>1533.2</td><td>0.997</td><td>4.13 ± 0.20</td><td>1.03 ± 0.06</td><td>288 ± 23</td><td>4.07 ± 0.21</td><td>11.37 ± 0.53*</td></tr><tr><td>TC-PD067</td><td>39.3015</td><td>175.5192</td><td>1533.6</td><td>0.997</td><td>5.24 ± 0.27</td><td>1.14 ± 0.05</td><td>331 ± 21</td><td>5.18 ± 0.26</td><td>14.09 ± 0.65</td></tr><tr><td>TC-PD068</td><td>39.3015</td><td>175.5192</td><td>1533.1</td><td>0.997</td><td>4.74 ± 0.23</td><td>1.09 ± 0.04</td><td>313 ± 18</td><td>4.68 ± 0.23</td><td>13.05 ± 0.62</td></tr><tr><td>TC-PD070</td><td>39.3012</td><td>175.5193</td><td>1528.0</td><td>0.997</td><td>4.89 ± 0.23</td><td>0.82 ± 0.04</td><td>431 ± 28</td><td>4.85 ± 0.23</td><td>13.27 ± 0.61</td></tr><tr><td colspan="10"><strong>TC – Turoa Cascades flow</strong></td></tr><tr><td>TC-PD066</td><td>39.3014</td><td>175.5193</td><td>1533.2</td><td>0.997</td><td>4.13 ± 0.20</td><td>1.03 ± 0.06</td><td>288 ± 23</td><td>4.07 ± 0.21</td><td>11.37 ± 0.53*</td></tr><tr><td>TC-PD067</td><td>39.3015</td><td>175.5192</td><td>1533.6</td><td>0.997</td><td>5.24 ± 0.27</td><td>1.14 ± 0.05</td><td>331 ± 21</td><td>5.18 ± 0.26</td><td>14.09 ± 0.65</td></tr><tr><td>TC-PD068</td><td>39.3015</td><td>175.5192</td><td>1533.1</td><td>0.997</td><td>4.74 ± 0.23</td><td>1.09 ± 0.04</td><td>313 ± 18</td><td>4.68 ± 0.23</td><td>13.05 ± 0.62</td></tr><tr><td>TC-PD070</td><td>39.3012</td><td>175.5193</td><td>1528.0</td><td>0.997</td><td>4.89 ± 0.23</td><td>0.82 ± 0.04</td><td>431 ± 28</td><td>4.85 ± 0.23</td><td>13.27 ± 0.61</td></tr></tbody></table>
<table><caption>Table 2. Continued.</caption><thead><tr><th>Sample</th><th>Latitude (S)</th><th>Longitude (E)</th><th>Elevation (m.s.l.)</th><th>Shielding factor</th><th><sup>3</sup>He<sub>tot</sub> ± 1σ (10<sup>6</sup> at g<sup>-1</sup>)</th><th><sup>4</sup>He<sub>tot</sub> ± 1σ (10<sup>10</sup> at g<sup>-1</sup>)</th><th>Total <sup>3</sup>He/<sup>4</sup>He (R/Ra)</th><th><sup>3</sup>He<sub>cos</sub> ± 1σ (10<sup>6</sup> at g<sup>-1</sup>)</th><th>Exposure age ± 1σ (ka)</th></tr></thead><tbody><tr><td colspan="10"><strong>MF – Makotuku Flats flow</strong></td></tr><tr><td>MF-PD061</td><td>39.3169</td><td>175.5143</td><td>1437.1</td><td>0.971</td><td>4.92 ± 0.26</td><td>1.93 ± 0.07</td><td>183 ± 12</td><td>4.82 ± 0.27</td><td>14.21 ± 0.73</td></tr><tr><td>MD-PD063</td><td>39.3168</td><td>175.5146</td><td>1434.8</td><td>0.991</td><td>4.00 ± 0.22</td><td>2.19 ± 0.09</td><td>131 ± 9</td><td>3.88 ± 0.23</td><td>11.45 ± 0.62</td></tr><tr><td>MF-PD064</td><td>39.3167</td><td>175.5146</td><td>1433.8</td><td>0.987</td><td>4.08 ± 0.22</td><td>1.65 ± 0.07</td><td>178 ± 12</td><td>3.99 ± 0.22</td><td>11.80 ± 0.61</td></tr><tr><td>MF-PD065</td><td>39.3167</td><td>175.5147</td><td>1433.3</td><td>0.988</td><td>4.47 ± 0.24</td><td>1.79 ± 0.08</td><td>180 ± 13</td><td>4.37 ± 0.25</td><td>12.80 ± 0.70</td></tr><tr><td colspan="10"><strong>NR – Ngā Rimutāmaka flow</strong></td></tr><tr><td>NR-PD053</td><td>39.3381</td><td>175.5873</td><td>1369.8</td><td>0.996</td><td>16.08 ± 0.67</td><td>2.46 ± 0.08</td><td>474 ± 25</td><td>16.05 ± 0.67</td><td>46.60 ± 2.12</td></tr><tr><td>NR-PD054</td><td>39.3384</td><td>175.5880</td><td>1372.9</td><td>1.000</td><td>15.34 ± 0.63</td><td>1.89 ± 0.07</td><td>582 ± 31</td><td>15.21 ± 0.64</td><td>43.38 ± 1.75</td></tr><tr><td>NR-PD055</td><td>39.3384</td><td>175.5879</td><td>1372.7</td><td>0.999</td><td>14.63 ± 0.62</td><td>1.79 ± 0.08</td><td>587 ± 35</td><td>14.52 ± 0.62</td><td>41.41 ± 1.53</td></tr><tr><td>NR-PD057</td><td>39.3384</td><td>175.5880</td><td>1372.6</td><td>0.995</td><td>14.80 ± 0.62</td><td>2.61 ± 0.10</td><td>408 ± 23</td><td>14.67 ± 0.63</td><td>42.11 ± 1.63</td></tr><tr><td colspan="10"><strong>MA – Makahikatoa flow</strong></td></tr><tr><td>MA-PD058</td><td>39.3125</td><td>175.6116</td><td>1594.8</td><td>0.996</td><td>20.03 ± 0.83</td><td>4.82 ± 0.15</td><td>299 ± 15</td><td>19.75 ± 0.84</td><td>48.96 ± 2.60</td></tr><tr><td>MA-PD059</td><td>39.3125</td><td>175.6116</td><td>1593.4</td><td>0.998</td><td>24.08 ± 1.18</td><td>9.56 ± 0.27</td><td>181 ± 10</td><td>23.58 ± 1.21</td><td>58.98 ± 2.84</td></tr><tr><td colspan="10"><strong>GR – Girdlestone Ridge flow</strong></td></tr><tr><td>GR-PD022 (Ol)</td><td>39.3072</td><td>175.5613</td><td>2148.0</td><td>0.996</td><td>7.88 ± 0.21</td><td>1.24 ± 0.10</td><td>457 ± 37</td><td>7.79 ± 0.21</td><td>14.02 ± 0.35</td></tr><tr><td>GR-PD023 (Ol)</td><td>39.3074</td><td>175.5615</td><td>2147.3</td><td>0.921</td><td>7.61 ± 0.49</td><td>1.19 ± 0.13</td><td>460 ± 59</td><td>7.54 ± 0.49</td><td>14.46 ± 0.85</td></tr><tr><td>GR-PD024 (Ol)</td><td>39.3074</td><td>175.5616</td><td>2145.4</td><td>0.990</td><td>6.19 ± 0.39</td><td>1.28 ± 0.07</td><td>348 ± 30</td><td>6.12 ± 0.39</td><td>11.11 ± 0.64*</td></tr><tr><td>GR-PD025 (Ol)</td><td>39.3078</td><td>175.5615</td><td>2128.1</td><td>0.993</td><td>3.61 ± 0.25</td><td>1.87 ± 0.10</td><td>139 ± 12</td><td>3.50 ± 0.25</td><td>6.81 ± 0.41*</td></tr><tr><td colspan="10"><strong>GR – Girdlestone Ridge flow</strong></td></tr><tr><td colspan="10"><strong>MA – Makahikatoa flow</strong></td></tr><tr><td>MA-PD058</td><td>39.3125</td><td>175.6116</td><td>1594.8</td><td>0.996</td><td>20.03 ± 0.83</td><td>4.82 ± 0.15</td><td>299 ± 15</td><td>19.75 ± 0.84</td><td>48.96 ± 2.60</td></tr><tr><td>MA-PD059</td><td>39.3125</td><td>175.6116</td><td>1593.4</td><td>0.998</td><td>24.08 ± 1.18</td><td>9.56 ± 0.27</td><td>181 ± 10</td><td>23.58 ± 1.21</td><td>58.98 ± 2.84</td></tr><tr><td colspan="10"><strong>GR – Girdlestone Ridge flow</strong></td></tr><tr><td>GR-PD022 (Ol)</td><td>39.3072</td><td>175.5613</td><td>2148.0</td><td>0.996</td><td>7.88 ± 0.21</td><td>1.24 ± 0.10</td><td>457 ± 37</td><td>7.79 ± 0.21</td><td>14.02 ± 0.35</td></tr><tr><td>GR-PD023 (Ol)</td><td>39.3074</td><td>175.5615</td><td>2147.3</td><td>0.921</td><td>7.61 ± 0.49</td><td>1.19 ± 0.13</td><td>460 ± 59</td><td>7.54 ± 0.49</td><td>14.46 ± 0.85</td></tr><tr><td>GR-PD024 (Ol)</td><td>39.3074</td><td>175.5616</td><td>2145.4</td><td>0.990</td><td>6.19 ± 0.39</td><td>1.28 ± 0.07</td><td>348 ± 30</td><td>6.12 ± 0.39</td><td>11.11 ± 0.64*</td></tr><tr><td>GR-PD025 (Ol)</td><td>39.3078</td><td>175.5615</td><td>2128.1</td><td>0.993</td><td>3.61 ± 0.25</td><td>1.87 ± 0.10</td><td>139 ± 12</td><td>3.50 ± 0.25</td><td>6.81 ± 0.41*</td></tr><tr><td colspan="10"><strong>GR – Girdlestone Ridge flow</strong></td></tr><tr><td>GR-PD022 (Ol)</td><td>39.3072</td><td>175.5613</td><td>2148.0</td><td>0.996</td><td>7.88 ± 0.21</td><td>1.24 ± 0.10</td><td>457 ± 37</td><td>7.79 ± 0.21</td><td>14.02 ± 0.35</td></tr><tr><td>GR-PD023 (Ol)</td><td>39.3074</td><td>175.5615</td><td>2147.3</td><td>0.921</td><td>7.61 ± 0.49</td><td>1.19 ± 0.13</td><td>460 ± 59</td><td>7.54 ± 0.49</td><td>14.46 ± 0.85</td></tr><tr><td>GR-PD024 (Ol)</td><td>39.3074</td><td>175.5616</td><td>2145.4</td><td>0.990</td><td>6.19 ± 0.39</td><td>1.28 ± 0.07</td><td>348 ± 30</td><td>6.12 ± 0.39</td><td>11.11 ± 0.64*</td></tr><tr><td>GR-PD025 (Ol)</td><td>39.3078</td><td>175.5615</td><td>2128.1</td><td>0.993</td><td>3.61 ± 0.25</td><td>1.87 ± 0.10</td><td>139 ± 12</td><td>3.50 ± 0.25</td><td>6.81 ± 0.41*</td></tr></tbody></table>

indicate older eruption ages than suggested by geochemical correlations.

## 5.2 Inconsistency with previous unit classification

Most of the eruption ages measured in this study are consistent with the age and geochemical ranges of the units to which they were assigned by Townsend et al. (2017). Here, we discuss the results we obtained that do not agree with the existing classification.

- Donoghue et al. (1999) linked the Pinnacle Ridge spatter-fed lava with the Taurewa pyroclastic unit (ca. 10 ka) based on geochemistry and the concentric nature of the isopachs of the Taurewa deposits around the location of PR. Our results indicate that the Pinnacle Ridge deposit was emplaced at 20.4 ± 4.0 ka, during the LGM and ca. 10 kyr prior to the Taurewa eruptive event, which is consistent with the lack of preservation of a proximal vent, likely associated with a significant erosive period and the retreat of large ice masses. Hence, our eruption age for Pinnacle Ridge suggests that this unit should be included as part of the Mangawhero Formation (50–15 ka) instead of the Whakapapa Formation (< 15 ka).

- MF samples were taken from a large flow considered to be part of the Makotuku Member of the Mangawhero Formation (ca. 24–16 ka; Table 1) based on its geochemistry. Our results show that this lava flow erupted at 12.6 ± 3.5 ka, which suggests that, based on age criteria, it could be classified as part of the Whakapapa Formation (< 15 ka).

- Our NR site was mapped as part of the Makotuku Member in an area dominated by outcrops of Mangaehuehu lavas (Fig. 4d). Our eruption age of 42.9 ± 8.6 ka for this site, together with the geochemical similarity of the NR samples to Mangaehuehu lavas (47–40 ka; Conway et al., 2016; see Table 1), suggests that the sampled outcrop is part of the Mangaehuehu Member.

- The outcrop we collected the MA samples from has, due to its geochemical similarity, been considered part of the Makotuku Member. Two exposure ages indicate that the Makahikatoa flow was emplaced at, or prior to, 50 ka, suggesting that it was formed during the first eruptive stages of the Mangawhero or in the late stages of the Waihianoa Formation (see Table 1) with a geochemical signature common in lavas emplaced at 24–16 ka.

- Exposure ages of GR samples (previously mapped as part of the Mangaehuehu Member) suggest that this lava

<footnote n="3">³He<sub>cos</sub> values were calculated using Eq. (3), with a magmatic ³He/⁴He of 5.9 ± 2.6 × 10⁻⁶ (~4.2 ± 1.9Ra). Individual samples are informed with 1σ for reproducibility using the CREp online calculator. Summary eruption age uncertainties represent 2σ values including production rate errors. Internal (INT) 2σ errors do not include production rate errors. All analysed samples consisted of pure pyroxenes with the exception of the site GR where analysed crystals were olivines with subordinate pyroxenes. For complete data and corrections, see Table A3. Outliers are marked with a * after the calculated exposure age. Two aliquots were measured for samples RTm-PD046 and CTb-PD234 for which we calculated a weighted mean of the ³He<sub>cos</sub> as a sample summary.</footnote>
<figure><img src="image_10.png" /><figcaption>Figure 5. Comparison between eruption ages obtained in this study and previous chronological constraints of the sampled flows. Unit colours correspond to the colours in Fig. 1. (a) Lavas < 20 ka. (b) Lava flows that are, or were thought to be, older than 20 ka.</figcaption></figure>

was emplaced during the last 15 kyr, which is inconsistent with it being part of the Mangawhero Formation. However, its geochemistry differentiates this outcrop from the rest of the Whakapapa lavas (Conway et al., 2016); thus it is likely to be part of a new member within the Whakapapa Formation.

- The results we obtained for flows from the Turoa Member indicate that lava was emplaced on Ruapehu's western flanks at ca. 15–12 ka (Turoa Cascades and Central Turoa-a flows as well as data from Conway et al., 2016) and, after a hiatus of ~4 kyr, again at around 8 ka (Mangaturuturu North and Central Turoa-b flows). Thus, we suggest the extension of the younger limit of the Turoa Member to 8 ka. Similarly, the obtained eruption ages redefine the age limits of the Rangataua Member (17–12 ka), Saddle Cone Member (12–8.5 ka), Taranaki Falls flow (16–13 ka), and Mangatoetoenui flows (12–9 ka).

## 5.3 Postglacial effusive activity of Ruapehu

Our ³Hecos-based eruption ages allow two periods of enhanced effusive activity since the LGM to be identified on Ruapehu (17–12 ka, Fig. 6a and b; and 9–7.5 ka, Fig. 6e), during which lava emplacement on different areas of the volcano occurred nearly simultaneously.

Our results show that, during the last glacial termination (ca. 17–14 ka; Fig. 6a), effusive activity affected the southern (Rangataua medial and, likely, the immense distal Rangataua flows of > 1.5 km³) and northern (Taranaki Falls flow) slopes of Ruapehu, suggesting that its southern and northern vents were active during this period. Radiometric dates published by Conway et al. (2016; see Table 1) suggest that, during this period, lava flows were also emplaced on Ruapehu's western (15.1 ± 2.4 ka, Turoa Member) and northwestern (14.8 ± 3.0 ka, Paretetaitonga Member) flanks. This period of generalized activity across Ruapehu continued until ca. 12 ka (Fig. 6b), with increasing intensity on the western flanks and decreasing intensity on the southern flanks. Eruption ages of the Whakapapaiti (13.3 ± 0.7 ka), Turoa Cascades (13.4 ± 0.7 ka), and Rangataua proximal (13.6 ± 0.6 ka) flows are nearly identical, indicating that lava emplacement occurred nearly simultaneously on different flanks of the volcanic edifice. In the Early Holocene (i.e. 12–10.5 ka, Fig. 6c), activity was focused on the east and northeast of the volcano, generating the first lavas of the Mangatoetoenui flows as well as lavas from satellite vents (Waiho-honu Plateau flow). After a flank collapse that affected part of the northern edifice at ca. 10.5 ka (Eaves et al., 2015), lava flows continued to be emplaced on the eastern flanks from the northern vent and erupted from satellite vents on the northeast in short time lapses (< 2 kyr), generating the large Saddle Cone flow (Fig. 6d). The rate of lava produc-
<figure><img src="image_11.png" /><figcaption>**Figure 6.** Lava flows emplaced at Ruapehu through time after the LGM. Collapse scars corresponds to flank collapse episodes at (d) 10.4–10.6 ka cal BP (Palmer and Neall, 1989; Eaves et al., 2015, Murimotu debris avalanche) and at (f) ca. 5.2 ka cal BP (Mangaio Formation, Donoghue, 1991; Donoghue and Neall, 2001). Lava flows with dotted boundaries in panels (a) and (e) have not been dated; their ages have been assigned based on geochemical and geomorphological similarities with dated flows.</figcaption></figure>

tion (i.e. number of individual lava flows produced) between 9 and 7.5 ka (Fig. 6e) was likely to have been the highest in the last 20 kyr at Ruapehu. Our results suggest that, during this time, most of the flows forming the Tawhainui sequence on north Ruapehu were emplaced from the northern vent, filling a topographic low left by the flank collapse. At a similar time, the last lavas of the Turoa Member (Mangaturuturu North and Central Turoa-b flows) were erupted from the southern vent and flowed to the west of the edifice. Effusive activity then declined, and, after another episode of

flank collapse that modified the topography surrounding the summit southern vent, lava flow emplacement was confined to the current outlet of Ruapehu's crater lake and flowed to the east (Whangaehu valley, Fig. 6f) at 2400–2050 BP (Greve et al., 2016).

Between ~23 and ~10 ka, Ruapehu produced at least five Plinian eruptions (as well as dozens of smaller explosive events) sourced from its northern vent (Pardo et al., 2012b). In contrast, effusive activity occurred from both the southern and northern vents until ~8 ka. Lack of high-resolution
ages of the pyroclastic deposits, however, hinders our ability to precisely compare the timing of these events. After this period of enhanced volcanism (finishing at ~10 ka for explosive events, Pardo et al., 2012b, and at ~8 ka for effusive events), activity at Ruapehu decreased significantly in magnitude and was restricted to the southern vent. However, our data expose time intervals during the last 17 kyr when lavas have been emplaced from both of Ruapehu's summit vents, challenging the assumption that volcanic hazards should be expected from the southern vent but not from the northern vent (e.g. Keys and Green, 2010; Leonard et al., 2021).

## 5.4 Applicability of cosmogenic ³He dating on stratovolcanoes

This study represents the first large-scale application of ³He cosmogenic dating tool for lava flows at stratovolcanoes. We provide ³He cosmogenic eruption age constraints for 20 young lava flows at Ruapehu, contributing to a detailed lava flow eruptive history for Ruapehu during the last 20 kyr (Sect. 5.3). Our data have good intra-flow clustering, inter-flow consistency, and good agreement with previous chronological constraints, demonstrating that robust eruption ages can be obtained for lava flows using ³He cosmogenic not only for basaltic lavas (e.g. Kurz et al., 1990; Licciardi et al., 2007; Foeken et al., 2009; Marchetti et al., 2014; Medynski et al., 2015) but also for andesitic lavas at stratovolcanoes.

Analyses of our samples yielded low ⁴Hetot values (likely influenced by repeated HF-leaching steps of the pyroxenes during sample preparation; Bromley et al., 2014) and low concentrations of radioactive elements (as expected from samples of intermediate compositions), which in turn resulted in small non-cosmogenic corrections and, added to analytical errors, small internal uncertainties in the obtained exposure ages. Like most other ³He cosmogenic ages, however, the ³He cosmogenic production rate uncertainty makes the largest contribution to our errors, imparting an uncertainty of ~10% to all calculated ages, which points out that more high-quality calibration sites are required to reduce these uncertainties and improve the quality of ³He-based exposure ages (Blard, 2021).

Considering these sources of uncertainties, our data show that the resolution of ³He cosmogenic eruption ages can be higher than ⁴⁰Ar/³⁹Ar or K/Ar for young intermediate lavas (see Fig. 5). In older lava flows (> 20 ka), radiometric methods can resolve emplacement ages more precisely (e.g. Lanphere, 2000; Harpel et al., 2004; Conway et al., 2016), whereas cosmogenic exposure ages become less certain due to production rate errors. Consequently, cosmogenic nuclide exposure dating has the potential to yield better results compared to ⁴⁰Ar/³⁹Ar or K/Ar when dating post-LGM lava flows (e.g. Harpel et al., 2004; Parmelee et al., 2015; Conway et al., 2016; Alcalá-Reygosa et al., 2018) and offers a valid alternative to date older lavas when no radiometric dating method can be applied (e.g. the site NR from this

study, whose age matches with higher-precision ⁴⁰Ar/³⁹Ar dates of geochemically similar lavas). Additionally, young lava flows are more likely to have original lava surfaces preserved as they were exposed to erosive and/or depositional processes for a relatively limited time. For the same reason, they are less likely to have exposed flow interiors needed for ⁴⁰Ar/³⁹Ar or K/Ar dating (Calvert and Lanphere, 2006; Fierstein et al., 2011), which makes ³He cosmogenic dating an ideal supplementary technique to radiometric methods when dating young pyroxene-bearing and olivine-bearing lavas at both basaltic volcanic areas and andesitic stratovolcanoes.

# 6 Conclusions

We analysed pyroxene-hosted and olivine-hosted ³He cosmogenic in 77 samples from 23 lava flows on Ruapehu volcano, Aotearoa/New Zealand, and obtained 16 eruption ages (between 7.8 ± 0.6 and 42.9 ± 1.7 ka; analytical $2\sigma$) and 7 minimum eruption ages, refining the chronology of lava flow emplacement at Ruapehu in the last 20 kyr.

Our data expose that weak ⁴⁰Ar/³⁹Ar isochrons led to unreliable eruption ages for two postglacial lavas at Ruapehu and stress the necessity of robust age constraints when using paleomagnetism as an age-refining tool.

Our results show that effusive activity at Ruapehu occurred from different vents during the last 17 kyr, affecting various sectors of the volcanic edifice over short time intervals. Based on our observations, we propose that the number of effusive eruptions during the last 20 kyr peaked at 17–12 and 9–7.5 ka. This represents a significant contribution to the hazard database of Aotearoa/New Zealand and valuable data for investigating temporal links of volcanic activity in the Tāupo Volcanic Zone.

Cosmogenic nuclide exposure dating can provide greater detail on the recent effusive chronology of stratovolcanoes, filling the gap left by the low resolution and challenges in acquiring adequate samples for radiometric dating of young lava flows.
Appendix A

<table><caption>Table A1. Abbreviations list used for sampling sites and samples.</caption><thead><tr><td>Abbreviation</td><td>Lava flow name</td><td>Area</td></tr></thead><tbody><tr><td>BR</td><td>Bruce Road</td><td>North</td></tr><tr><td>CTa</td><td>Central Turoa-a</td><td>West</td></tr><tr><td>CTb</td><td>Central Turoa-b</td><td>West</td></tr><tr><td>DC</td><td>Delta Corner</td><td>North</td></tr><tr><td>GR</td><td>Girdlestone Ridge</td><td>South</td></tr><tr><td>LC</td><td>Lava Cascade</td><td>East</td></tr><tr><td>MA</td><td>Makahikatoa</td><td>Southeast</td></tr><tr><td>MF</td><td>Makotuku Flat</td><td>West</td></tr><tr><td>MN</td><td>Mangaturuturu North</td><td>West</td></tr><tr><td>MS</td><td>Mangaturuturu South</td><td>West</td></tr><tr><td>NR</td><td>Ngā Rimutāmaka</td><td>South</td></tr><tr><td>PR</td><td>Pinnacle Ridge</td><td>North</td></tr><tr><td>RTm</td><td>Rangataua medial</td><td>South</td></tr><tr><td>RTp</td><td>Rangataua proximal</td><td>South</td></tr><tr><td>SC</td><td>Saddle Cone</td><td>Northeast</td></tr><tr><td>SCw</td><td>Saddle Cone – western lobe</td><td>Northeast</td></tr><tr><td>SCe</td><td>Saddle Cone – eastern lobe</td><td>Northeast</td></tr><tr><td>TC</td><td>Turoa Cascades</td><td>West</td></tr><tr><td>TFa</td><td>Taranaki Falls</td><td>North</td></tr><tr><td>TFt</td><td>Tukino Flats</td><td>East</td></tr><tr><td>TSa</td><td>Tukino Slopes-a</td><td>East</td></tr><tr><td>TSb</td><td>Tukino Slopes-b</td><td>East</td></tr><tr><td>WG</td><td>Whakapapa Glacier</td><td>North</td></tr><tr><td>WP</td><td>Waihohonu Plateau</td><td>Northeast</td></tr><tr><td>WT</td><td>Whakapapaiti</td><td>Northwest</td></tr></tbody></table>
<table><caption>Table A2. Normalized major and trace elements of bulk rock and analysed minerals for each sampled lava flow.</caption><thead><tr><th colspan="2">Bulk rock</th><th colspan="10">Normalized wt %</th><th colspan="10">Parts per million</th></tr><tr><th>Site</th><th>SiO<sub>2</sub></th><th>Al<sub>2</sub>O<sub>3</sub></th><th>Fe<sub>2</sub>O<sub>3</sub></th><th>MnO</th><th>MgO</th><th>CaO</th><th>Na<sub>2</sub>O</th><th>K<sub>2</sub>O</th><th>TiO<sub>2</sub></th><th>P<sub>2</sub>O<sub>5</sub></th><th>LOI</th><th>Majot</th><th>Li</th><th>B</th><th>Cr</th><th>Co</th><th>Ni</th><th>Gd</th><th>Sm</th><th>U</th><th>Th</th></tr></thead><tbody><tr><td>BR</td><td>57.01</td><td>15.88</td><td>8.23</td><td>0.13</td><td>5.37</td><td>7.42</td><td>3.14</td><td>1.40</td><td>0.68</td><td>0.18</td><td>0.56</td><td>99.26</td><td>17.9</td><td>20.0</td><td>85.19</td><td>26.68</td><td>45.59</td><td>2.69</td><td>2.75</td><td>1.08</td><td>4.04</td></tr><tr><td>WG</td><td>55.80</td><td>16.49</td><td>7.69</td><td>0.12</td><td>4.74</td><td>6.38</td><td>3.11</td><td>1.55</td><td>0.68</td><td>0.14</td><td>3.31</td><td>100.23</td><td>20.4</td><td>22.8</td><td>101.82</td><td>23.64</td><td>40.78</td><td>2.79</td><td>2.83</td><td>1.30</td><td>4.91</td></tr><tr><td>LC</td><td>54.87</td><td>16.02</td><td>9.25</td><td>0.13</td><td>5.34</td><td>6.59</td><td>2.92</td><td>1.41</td><td>0.72</td><td>0.14</td><td>2.62</td><td>99.46</td><td>17.6</td><td>20.1</td><td>91.32</td><td>26.95</td><td>45.80</td><td>2.75</td><td>2.75</td><td>1.13</td><td>4.36</td></tr><tr><td>TSa</td><td>52.61</td><td>18.94</td><td>8.79</td><td>0.13</td><td>4.21</td><td>5.18</td><td>2.65</td><td>1.52</td><td>0.78</td><td>0.15</td><td>5.04</td><td>99.79</td><td>19.6</td><td>20.5</td><td>45.56</td><td>21.82</td><td>20.55</td><td>2.87</td><td>2.84</td><td>1.46</td><td>5.22</td></tr><tr><td>TSb</td><td>53.97</td><td>17.20</td><td>8.91</td><td>0.13</td><td>4.36</td><td>6.51</td><td>2.81</td><td>1.41</td><td>0.79</td><td>0.16</td><td>3.76</td><td>99.52</td><td>16.3</td><td>19.6</td><td>45.17</td><td>22.39</td><td>14.11</td><td>2.79</td><td>2.81</td><td>1.19</td><td>4.79</td></tr><tr><td>TFt</td><td>52.75</td><td>18.33</td><td>8.68</td><td>0.13</td><td>4.62</td><td>5.87</td><td>2.64</td><td>1.42</td><td>0.77</td><td>0.14</td><td>4.65</td><td>100.47</td><td>19.0</td><td>20.0</td><td>85.21</td><td>24.29</td><td>24.85</td><td>2.79</td><td>2.76</td><td>1.42</td><td>5.03</td></tr><tr><td>TFa</td><td>57.17</td><td>16.80</td><td>7.81</td><td>0.12</td><td>4.09</td><td>6.88</td><td>3.15</td><td>1.52</td><td>0.72</td><td>0.13</td><td>1.61</td><td>100.25</td><td>16.8</td><td>22.3</td><td>60.99</td><td>20.34</td><td>19.91</td><td>2.91</td><td>2.96</td><td>1.34</td><td>5.09</td></tr><tr><td>SC</td><td>56.76</td><td>16.48</td><td>7.88</td><td>0.12</td><td>4.64</td><td>7.02</td><td>3.07</td><td>1.52</td><td>0.72</td><td>0.15</td><td>1.63</td><td>99.87</td><td>21.3</td><td>21.7</td><td>86.23</td><td>23.75</td><td>31.36</td><td>2.77</td><td>2.81</td><td>1.32</td><td>4.92</td></tr><tr><td>WP</td><td>55.42</td><td>16.60</td><td>8.13</td><td>0.12</td><td>4.87</td><td>6.58</td><td>2.77</td><td>1.31</td><td>0.70</td><td>0.13</td><td>3.37</td><td>100.04</td><td>18.8</td><td>20.3</td><td>95.47</td><td>24.80</td><td>36.95</td><td>2.72</td><td>2.60</td><td>1.18</td><td>4.49</td></tr><tr><td>PR</td><td>58.59</td><td>16.12</td><td>6.78</td><td>0.10</td><td>4.42</td><td>5.69</td><td>3.04</td><td>1.79</td><td>0.65</td><td>0.14</td><td>2.68</td><td>100.29</td><td>22.5</td><td>21.9</td><td>183.67</td><td>20.39</td><td>53.61</td><td>2.57</td><td>2.50</td><td>1.45</td><td>5.50</td></tr><tr><td>RFt</td><td>56.56</td><td>17.39</td><td>7.60</td><td>0.12</td><td>3.37</td><td>5.35</td><td>3.15</td><td>1.74</td><td>0.75</td><td>0.15</td><td>3.84</td><td>99.83</td><td>19.2</td><td>23.9</td><td>36.93</td><td>16.88</td><td>12.67</td><td>3.00</td><td>3.11</td><td>1.49</td><td>5.91</td></tr><tr><td>RFm</td><td>56.64</td><td>16.82</td><td>7.67</td><td>0.12</td><td>3.57</td><td>5.72</td><td>3.15</td><td>1.69</td><td>0.75</td><td>0.21</td><td>3.67</td><td>99.53</td><td>18.6</td><td>23.4</td><td>40.90</td><td>17.75</td><td>13.20</td><td>3.15</td><td>3.17</td><td>1.53</td><td>6.03</td></tr><tr><td>WT</td><td>55.32</td><td>17.34</td><td>8.29</td><td>0.12</td><td>4.19</td><td>5.79</td><td>2.85</td><td>1.52</td><td>0.76</td><td>0.14</td><td>3.69</td><td>100.24</td><td>20.3</td><td>20.9</td><td>57.45</td><td>20.84</td><td>18.59</td><td>2.76</td><td>2.80</td><td>1.32</td><td>5.40</td></tr><tr><td>MN</td><td>53.66</td><td>17.72</td><td>8.07</td><td>0.12</td><td>5.01</td><td>6.09</td><td>2.93</td><td>1.47</td><td>0.66</td><td>0.13</td><td>4.14</td><td>100.44</td><td>19.2</td><td>22.5</td><td>93.87</td><td>25.17</td><td>47.08</td><td>2.42</td><td>2.41</td><td>1.34</td><td>4.86</td></tr><tr><td>MS</td><td>55.58</td><td>15.71</td><td>9.54</td><td>0.13</td><td>4.51</td><td>6.18</td><td>2.87</td><td>1.49</td><td>0.84</td><td>0.13</td><td>3.01</td><td>100.20</td><td>13.8</td><td>19.2</td><td>54.80</td><td>23.51</td><td>17.53</td><td>2.45</td><td>2.45</td><td>0.94</td><td>4.99</td></tr><tr><td>Cta</td><td>56.67</td><td>16.49</td><td>8.12</td><td>0.12</td><td>4.36</td><td>6.49</td><td>3.01</td><td>1.50</td><td>0.72</td><td>0.15</td><td>2.37</td><td>99.64</td><td>20.2</td><td>21.8</td><td>88.79</td><td>22.15</td><td>29.15</td><td>2.75</td><td>2.80</td><td>1.35</td><td>5.02</td></tr><tr><td>CTh</td><td>56.01</td><td>16.78</td><td>8.32</td><td>0.12</td><td>4.60</td><td>7.12</td><td>3.01</td><td>1.39</td><td>0.74</td><td>0.13</td><td>1.77</td><td>99.29</td><td>19.3</td><td>20.8</td><td>89.86</td><td>23.46</td><td>27.88</td><td>2.72</td><td>2.75</td><td>1.20</td><td>4.51</td></tr><tr><td>TC</td><td>56.18</td><td>16.49</td><td>8.44</td><td>0.13</td><td>4.60</td><td>7.03</td><td>3.06</td><td>1.48</td><td>0.78</td><td>0.14</td><td>1.66</td><td>99.90</td><td>17.3</td><td>20.0</td><td>71.52</td><td>22.48</td><td>18.63</td><td>2.96</td><td>2.95</td><td>1.22</td><td>4.71</td></tr><tr><td>MF</td><td>58.67</td><td>16.98</td><td>6.52</td><td>0.093</td><td>2.96</td><td>5.91</td><td>3.50</td><td>1.87</td><td>0.84</td><td>0.16</td><td>2.50</td><td>100.13</td><td>26.8</td><td>27.3</td><td>49.22</td><td>16.60</td><td>18.75</td><td>2.96</td><td>3.13</td><td>1.63</td><td>6.12</td></tr><tr><td>NR</td><td>57.92</td><td>14.71</td><td>6.81</td><td>0.11</td><td>6.24</td><td>5.96</td><td>2.96</td><td>1.85</td><td>0.73</td><td>0.19</td><td>2.53</td><td>99.77</td><td>29.7</td><td>23.4</td><td>308.14</td><td>25.08</td><td>109.34</td><td>3.14</td><td>3.25</td><td>1.74</td><td>6.50</td></tr><tr><td>MA</td><td>58.73</td><td>17.28</td><td>6.77</td><td>0.10</td><td>3.09</td><td>5.26</td><td>3.29</td><td>1.55</td><td>0.68</td><td>0.13</td><td>3.12</td><td>99.82</td><td>21.0</td><td>24.4</td><td>34.02</td><td>17.42</td><td>14.05</td><td>2.77</td><td>2.85</td><td>1.31</td><td>5.00</td></tr><tr><td>GR</td><td>55.76</td><td>15.21</td><td>7.40</td><td>0.12</td><td>6.09</td><td>6.77</td><td>2.80</td><td>1.35</td><td>0.66</td><td>0.14</td><td>3.71</td><td>100.13</td><td>20.4</td><td>20.7</td><td>215.77</td><td>26.79</td><td>73.65</td><td>2.71</td><td>2.78</td><td>1.14</td><td>4.50</td></tr><tr><td colspan="12">Minerals</td></tr><tr><td>DC</td><td>52.11</td><td>1.58</td><td>20.53</td><td>0.41</td><td>22.11</td><td>4.08</td><td>0.09</td><td>&lt;DL</td><td>0.31</td><td>&lt;DL</td><td>-1.20</td><td>99.63</td><td>&lt;2</td><td>10.7</td><td>285.89</td><td>98.63</td><td>187.40</td><td>1.07</td><td>0.78</td><td>0.01</td><td>0.040</td></tr><tr><td>BR</td><td>51.76</td><td>1.97</td><td>18.03</td><td>0.36</td><td>19.70</td><td>8.31</td><td>0.20</td><td>0.04</td><td>0.39</td><td>&lt;DL</td><td>-0.76</td><td>99.55</td><td>&lt;2</td><td>8.59</td><td>398.47</td><td>85.29</td><td>172.81</td><td>2.37</td><td>1.97</td><td>0.04</td><td>0.144</td></tr><tr><td>WG</td><td>51.71</td><td>1.61</td><td>20.15</td><td>0.39</td><td>21.71</td><td>4.76</td><td>0.10</td><td>&lt;DL</td><td>0.34</td><td>&lt;DL</td><td>-0.77</td><td>99.85</td><td>&lt;2</td><td>12.8</td><td>371.23</td><td>92.22</td><td>182.44</td><td>1.38</td><td>1.08</td><td>0.01</td><td>0.048</td></tr><tr><td>LC</td><td>51.54</td><td>1.59</td><td>20.41</td><td>0.43</td><td>20.50</td><td>5.74</td><td>0.12</td><td>&lt;DL</td><td>0.39</td><td>&lt;DL</td><td>-0.71</td><td>99.32</td><td>&lt;2</td><td>10.7</td><td>337.72</td><td>85.76</td><td>129.55</td><td>1.97</td><td>1.60</td><td>0.02</td><td>0.056</td></tr><tr><td>TSa</td><td>49.38</td><td>1.51</td><td>24.99</td><td>0.44</td><td>21.19</td><td>2.14</td><td>0.05</td><td>&lt;DL</td><td>0.86</td><td>&lt;DL</td><td>-0.57</td><td>99.82</td><td>&lt;2</td><td>7.35</td><td>404.65</td><td>88.37</td><td>143.62</td><td>0.65</td><td>0.44</td><td>0.02</td><td>0.052</td></tr><tr><td>TSb</td><td>51.60</td><td>1.52</td><td>21.23</td><td>0.45</td><td>21.34</td><td>4.16</td><td>0.08</td><td>&lt;DL</td><td>0.37</td><td>&lt;DL</td><td>-0.76</td><td>100.33</td><td>&lt;2</td><td>10.9</td><td>277.49</td><td>88.05</td><td>117.96</td><td>1.39</td><td>1.08</td><td>0.01</td><td>0.050</td></tr><tr><td>TFt</td><td>51.18</td><td>1.53</td><td>21.60</td><td>0.45</td><td>20.94</td><td>4.64</td><td>0.09</td><td>&lt;DL</td><td>0.49</td><td>&lt;DL</td><td>-0.90</td><td>99.54</td><td>&lt;2</td><td>16.7</td><td>279.53</td><td>88.18</td><td>120.06</td><td>1.59</td><td>1.27</td><td>0.01</td><td>0.047</td></tr><tr><td>TFa</td><td>51.27</td><td>1.51</td><td>20.55</td><td>0.42</td><td>20.94</td><td>5.09</td><td>0.10</td><td>&lt;DL</td><td>0.40</td><td>&lt;DL</td><td>-0.28</td><td>99.10</td><td>&lt;2</td><td>11.0</td><td>367.38</td><td>85.09</td><td>130.22</td><td>1.65</td><td>1.35</td><td>0.02</td><td>0.062</td></tr><tr><td>SC</td><td>51.41</td><td>1.91</td><td>19.70</td><td>0.39</td><td>19.54</td><td>7.18</td><td>0.18</td><td>0.04</td><td>0.42</td><td>&lt;DL</td><td>-0.78</td><td>100.00</td><td>&lt;2</td><td>9.17</td><td>394.61</td><td>80.89</td><td>127.03</td><td>2.56</td><td>2.07</td><td>0.05</td><td>0.185</td></tr><tr><td>WP</td><td>52.19</td><td>2.37</td><td>16.91</td><td>0.33</td><td>19.95</td><td>8.16</td><td>0.28</td><td>0.10</td><td>0.38</td><td>&lt;DL</td><td>-0.67</td><td>99.86</td><td>&lt;2</td><td>10.2</td><td>645.87</td><td>77.50</td><td>164.79</td><td>2.21</td><td>1.80</td><td>0.10</td><td>0.358</td></tr><tr><td>PR</td><td>51.05</td><td>2.01</td><td>20.06</td><td>0.39</td><td>20.02</td><td>6.31</td><td>0.13</td><td>&lt;DL</td><td>0.41</td><td>&lt;DL</td><td>-0.39</td><td>99.57</td><td>&lt;2</td><td>31.0</td><td>632.12</td><td>83.14</td><td>163.52</td><td>1.94</td><td>1.54</td><td>0.03</td><td>0.141</td></tr><tr><td>RFt</td><td>51.08</td><td>1.50</td><td>22.65</td><td>0.43</td><td>21.53</td><td>3.04</td><td>0.07</td><td>&lt;DL</td><td>0.49</td><td>&lt;DL</td><td>-0.80</td><td>100.29</td><td>&lt;2</td><td>5.71</td><td>405.44</td><td>91.30</td><td>153.24</td><td>1.06</td><td>0.80</td><td>0.01</td><td>0.049</td></tr><tr><td>RFm</td><td>51.38</td><td>1.46</td><td>23.10</td><td>0.44</td><td>21.71</td><td>2.91</td><td>0.06</td><td>&lt;DL</td><td>0.47</td><td>&lt;DL</td><td>-1.53</td><td>99.84</td><td>&lt;2</td><td>7.48</td><td>321.50</td><td>91.05</td><td>135.00</td><td>1.07</td><td>0.82</td><td>0.01</td><td>0.045</td></tr><tr><td>WT</td><td>47.28</td><td>1.59</td><td>26.88</td><td>0.43</td><td>20.27</td><td>2.90</td><td>0.06</td><td>&lt;DL</td><td>1.41</td><td>&lt;DL</td><td>-0.82</td><td>99.76</td><td>&lt;2</td><td>10.7</td><td>465.85</td><td>96.87</td><td>142.07</td><td>0.89</td><td>0.68</td><td>0.03</td><td>0.069</td></tr><tr><td>MN</td><td>52.02</td><td>1.51</td><td>21.69</td><td>0.43</td><td>22.68</td><td>2.47</td><td>0.06</td><td>&lt;DL</td><td>0.31</td><td>&lt;DL</td><td>-1.18</td><td>99.37</td><td>&lt;2</td><td>13.4</td><td>269.95</td><td>102.76</td><td>196.34</td><td>0.66</td><td>0.46</td><td>0.01</td><td>0.040</td></tr><tr><td>MS</td><td>51.76</td><td>1.62</td><td>20.40</td><td>0.42</td><td>21.30</td><td>5.00</td><td>0.10</td><td>&lt;DL</td><td>0.38</td><td>&lt;DL</td><td>-0.97</td><td>99.36</td><td>&lt;2</td><td>23.6</td><td>374.34</td><td>88.31</td><td>122.86</td><td>1.52</td><td>1.20</td><td>0.01</td><td>0.046</td></tr><tr><td>Cta</td><td>51.45</td><td>1.64</td><td>20.05</td><td>0.41</td><td>20.42</td><td>5.92</td><td>0.13</td><td>&lt;DL</td><td>0.38</td><td>&lt;DL</td><td>-0.40</td><td>99.89</td><td>&lt;2</td><td>9.23</td><td>409.65</td><td>84.51</td><td>143.50</td><td>2.11</td><td>1.72</td><td>0.03</td><td>0.098</td></tr><tr><td>CTh</td><td>51.08</td><td>1.63</td><td>20.86</td><td>0.41</td><td>20.43</td><td>5.87</td><td>0.12</td><td>&lt;DL</td><td>0.62</td><td>&lt;DL</td><td>-1.03</td><td>99.35</td><td>&lt;2</td><td>10.8</td><td>410.51</td><td>84.63</td><td>144.39</td><td>2.03</td><td>1.67</td><td>0.02</td><td>0.065</td></tr><tr><td>TC</td><td>51.56</td><td>1.62</td><td>20.37</td><td>0.41</td><td>21.10</td><td>5.02</td><td>0.10</td><td>&lt;DL</td><td>0.35</td><td>&lt;DL</td><td>-0.53</td><td>100.21</td><td>&lt;2</td><td>8.01</td><td>391.62</td><td>88.11</td><td>124.68</td><td>1.51</td><td>1.19</td><td>0.01</td><td>0.046</td></tr><tr><td>MF</td><td>50.94</td><td>1.66</td><td>22.07</td><td>0.39</td><td>21.93</td><td>3.35</td><td>0.07</td><td>&lt;DL</td><td>0.82</td><td>&lt;DL</td><td>-1.22</td><td>99.65</td><td>&lt;2</td><td>13.6</td><td>475.60</td><td>95.40</td><td>181.14</td><td>0.94</td><td>0.73</td><td>0.02</td><td>0.058</td></tr><tr><td>NR</td><td>52.64</td><td>1.95</td><td>16.10</td><td>0.28</td><td>22.62</td><td>6.82</td><td>0.15</td><td>&lt;DL</td><td>0.35</td><td>&lt;DL</td><td>-0.90</td><td>99.71</td><td>&lt;2</td><td>11.7</td><td>1184.21</td><td>83.92</td><td>352.04</td><td>1.56</td><td>1.27</td><td>0.02</td><td>0.042</td></tr><tr><td>MA</td><td>49.70</td><td>2.03</td><td>22.77</td><td>0.39</td><td>21.09</td><td>3.88</td><td>0.08</td><td>&lt;DL</td><td>0.89</td><td>&lt;DL</td><td>-0.83</td><td>99.44</td><td>&lt;2</td><td>17.7</td><td>332.48</td><td>92.02</td><td>106.83</td><td>0.88</td><td>0.64</td><td>0.02</td><td>0.081</td></tr><tr><td>GR</td><td>54.07</td><td>1.52</td><td>13.92</td><td>0.27</td><td>27.48</td><td>2.89</td><td>0.05</td><td>&lt;DL</td><td>0.19</td><td>&lt;DL</td><td>-0.39</td><td>99.30</td><td>&lt;5</td><td>6.22</td><td>1649.09</td><td>87.13</td><td>360.96</td><td>0.36</td><td>0.26</td><td>&lt;DL</td><td>0.022</td></tr></tbody></table>

Detection limits (DL) are at 0.03 wt % for K₂O, 0.10 wt % for P₂O₅, 2 ppm for B, and 0.01 ppm for U.
<table><caption>Table A3. Sample data used to compute exposure ages.</caption><thead><tr><th>Sample</th><th>Latitude (S)</th><th>Longitude (E)</th><th>Elevation (m.s.l.)</th><th>Surface dip (°)</th><th>Dip direction (°)</th><th>Shielding factor</th><th>Density (g cm<sup>-3</sup>)</th><th>Thickness (cm)</th><th>P<sub>nuc</sub> (10<sup>2</sup> at g<sup>-1</sup> yr<sup>-1</sup>)</th><th>Closure age (Ma)</th><th>He<sub>nuc</sub> (at)</th><th>P<sub>3</sub> (at g<sup>-1</sup> yr<sup>-1</sup>)</th><th>P<sub>4</sub> (10<sup>5</sup> at g<sup>-1</sup> yr<sup>-1</sup>)</th><th>R factor</th></tr></thead><tbody><tr><td>SC-PD001</td><td>39.2143</td><td>175.6011</td><td>1439.0</td><td>—</td><td>—</td><td>0.998</td><td>1.89</td><td>3.6</td><td>406</td><td>0.010</td><td>406</td><td>231.22</td><td>3.13</td><td>0.9932</td></tr><tr><td>SC-PD002</td><td>39.2143</td><td>175.6010</td><td>1439.3</td><td>—</td><td>—</td><td>0.998</td><td>2.01</td><td>3.2</td><td></td><td></td><td></td><td>245.83</td><td></td><td>0.9932</td></tr><tr><td>SC-PD003</td><td>39.2146</td><td>175.5997</td><td>1443.3</td><td>—</td><td>—</td><td>0.998</td><td>2.01</td><td>2.4</td><td></td><td></td><td></td><td>245.86</td><td></td><td>0.9932</td></tr><tr><td>WP-PD007</td><td>39.2479</td><td>175.5882</td><td>1911.7</td><td>10</td><td>190</td><td>0.996</td><td>2.13</td><td>3.3</td><td>414</td><td>0.010</td><td>415</td><td>487.90</td><td>6.08</td><td>0.9906</td></tr><tr><td>WP 008</td><td>39.2479</td><td>175.5882</td><td>1912.1</td><td>18</td><td>30</td><td>0.995</td><td>2.06</td><td>3.0</td><td></td><td></td><td></td><td>487.90</td><td></td><td>0.9906</td></tr><tr><td>BR-PD014</td><td>39.2201</td><td>175.5406</td><td>1360.0</td><td>—</td><td>—</td><td>0.999</td><td>2.15</td><td>2.7</td><td>565</td><td>0.010</td><td>565</td><td>317.93</td><td>0.75</td><td>0.9982</td></tr><tr><td>BR-PD016</td><td>39.2198</td><td>175.5379</td><td>1359.2</td><td>28</td><td>55</td><td>0.981</td><td>2.32</td><td>2.6</td><td></td><td></td><td></td><td>320.37</td><td></td><td>0.9982</td></tr><tr><td>BR-PD017</td><td>39.2190</td><td>175.5409</td><td>1332.6</td><td>—</td><td>—</td><td>0.998</td><td>2.25</td><td>2.7</td><td></td><td></td><td></td><td>313.04</td><td></td><td>0.9982</td></tr><tr><td>BR-PD018</td><td>39.2190</td><td>175.5411</td><td>1332.4</td><td>—</td><td>—</td><td>0.998</td><td>2.15</td><td>4.3</td><td></td><td></td><td></td><td>317.93</td><td></td><td>0.9982</td></tr><tr><td>GR-PD022</td><td>39.3072</td><td>175.5613</td><td>2148.0</td><td>—</td><td>—</td><td>0.996</td><td>2.15</td><td>2.7</td><td>254</td><td>0.020</td><td>508</td><td>589.39</td><td>0.44</td><td>0.9994</td></tr><tr><td>GR-PD023</td><td>39.3072</td><td>175.5615</td><td>2147.2</td><td>45</td><td>180</td><td>0.921</td><td>2.24</td><td>4.5</td><td></td><td></td><td></td><td>584.50</td><td></td><td>0.9994</td></tr><tr><td>GR-PD024</td><td>39.3074</td><td>175.5616</td><td>2145.4</td><td>—</td><td>—</td><td>0.990</td><td>2.18</td><td>4.6</td><td></td><td></td><td></td><td>574.72</td><td></td><td>0.9994</td></tr><tr><td>GR-PD025</td><td>39.3078</td><td>175.5615</td><td>2128.1</td><td>—</td><td>—</td><td>0.993</td><td>2.80</td><td>4.7</td><td></td><td></td><td></td><td>539.25</td><td></td><td>0.9994</td></tr><tr><td>RT-PD027</td><td>39.3140</td><td>175.5509</td><td>1831.4</td><td>—</td><td>—</td><td>0.997</td><td>1.84*</td><td>5.5</td><td>296</td><td>0.015</td><td>444</td><td>468.33</td><td>0.74</td><td>0.9988</td></tr><tr><td>RT-PD028</td><td>39.3140</td><td>175.5509</td><td>1833.1</td><td>—</td><td>—</td><td>0.996</td><td>1.79</td><td>5.3</td><td></td><td></td><td></td><td>467.11</td><td></td><td>0.9988</td></tr><tr><td>RT 029</td><td>39.3140</td><td>175.5509</td><td>1832.9</td><td>—</td><td>—</td><td>0.996</td><td>1.84*</td><td>5.3</td><td></td><td></td><td></td><td>469.56</td><td></td><td>0.9988</td></tr><tr><td>RT-PD030</td><td>39.3143</td><td>175.5512</td><td>1816.4</td><td>20</td><td>230</td><td>0.988</td><td>1.89</td><td>5.1</td><td></td><td></td><td></td><td>465.89</td><td></td><td>0.9988</td></tr><tr><td>RT-PD045</td><td>39.3234</td><td>175.5520</td><td>1585.9</td><td>20</td><td>80</td><td>0.991</td><td>2.80</td><td>4.2</td><td>401</td><td>0.015</td><td>601</td><td>399.86</td><td>0.75</td><td>0.9986</td></tr><tr><td>RT-PD046</td><td>39.3249</td><td>175.5508</td><td>1567.6</td><td>29</td><td>145</td><td>0.979</td><td>2.07</td><td>3.6</td><td></td><td></td><td></td><td>393.74</td><td></td><td>0.9986</td></tr><tr><td>RT-PD047</td><td>39.3251</td><td>175.5503</td><td>1567.4</td><td>—</td><td>—</td><td>0.997</td><td>2.27</td><td>3.0</td><td></td><td></td><td></td><td>392.52</td><td></td><td>0.9986</td></tr><tr><td>RT-PD048</td><td>39.3250</td><td>175.5503</td><td>1567.3</td><td>—</td><td>—</td><td>0.997</td><td>2.19</td><td>3.4</td><td></td><td></td><td></td><td>374.18</td><td></td><td>0.9985</td></tr><tr><td>NR-PD053</td><td>39.3381</td><td>175.5873</td><td>1369.8</td><td>15</td><td>3</td><td>0.996</td><td>2.39</td><td>4.8</td><td>654</td><td>0.045</td><td>2944</td><td>359.50</td><td>1.01</td><td>0.9979</td></tr><tr><td>NR-PD054</td><td>39.3384</td><td>175.5880</td><td>1372.9</td><td>—</td><td>—</td><td>1.000</td><td>2.06</td><td>4.2</td><td></td><td></td><td></td><td>360.73</td><td></td><td>0.9979</td></tr><tr><td>NR-PD055</td><td>39.3384</td><td>175.5879</td><td>1372.5</td><td>—</td><td>—</td><td>0.999</td><td>2.15</td><td>3.2</td><td></td><td></td><td></td><td>358.28</td><td></td><td>0.9979</td></tr><tr><td>NR-PD057</td><td>39.3384</td><td>175.5880</td><td>1372.6</td><td>—</td><td>—</td><td>0.995</td><td>2.47</td><td>3.5</td><td></td><td></td><td></td><td>359.50</td><td></td><td>0.9979</td></tr><tr><td>MA-PD058</td><td>39.3125</td><td>175.6116</td><td>1570.7</td><td>14</td><td>190</td><td>0.996</td><td>2.45</td><td>2.9</td><td>818</td><td>0.050</td><td>4089</td><td>414.53</td><td>1.29</td><td>0.9977</td></tr><tr><td>MA-PD059</td><td>39.3125</td><td>175.6116</td><td>1569.3</td><td>—</td><td>—</td><td>0.998</td><td>2.21</td><td>3.20</td><td></td><td></td><td></td><td>408.42</td><td></td><td>0.9976</td></tr><tr><td>MF-PD061</td><td>39.3169</td><td>175.5143</td><td>1437.0</td><td>—</td><td>—</td><td>0.971</td><td>2.15</td><td>3.0</td><td>729</td><td>0.013</td><td>948</td><td>353.39</td><td>1.14</td><td>0.9975</td></tr><tr><td>MF-PD063</td><td>39.3168</td><td>175.5146</td><td>1434.8</td><td>—</td><td>—</td><td>0.991</td><td>1.81</td><td>3.8</td><td></td><td></td><td></td><td>348.50</td><td></td><td>0.9975</td></tr><tr><td>MF-PD064</td><td>39.3167</td><td>175.5146</td><td>1433.8</td><td>—</td><td>—</td><td>0.987</td><td>1.95</td><td>3.6</td><td></td><td></td><td></td><td>349.72</td><td></td><td>0.9975</td></tr><tr><td>MF-PD065</td><td>39.3167</td><td>175.5147</td><td>1433.3</td><td>11</td><td>180</td><td>0.988</td><td>1.80</td><td>2.6</td><td></td><td></td><td></td><td>350.94</td><td></td><td>0.9975</td></tr><tr><td>TC-PD066</td><td>39.3014</td><td>175.5193</td><td>1533.2</td><td>—</td><td>—</td><td>0.997</td><td>2.17</td><td>6.2</td><td>292</td><td>0.013</td><td>380</td><td>375.40</td><td>7.59</td><td>0.9985</td></tr><tr><td>TC-PD067</td><td>39.3015</td><td>175.5192</td><td>1533.6</td><td>—</td><td>—</td><td>0.997</td><td>2.11</td><td>4.5</td><td></td><td></td><td></td><td>379.07</td><td></td><td>0.9985</td></tr><tr><td>TC-PD068</td><td>39.3015</td><td>175.5192</td><td>1533.1</td><td>—</td><td>—</td><td>0.997</td><td>2.15</td><td>7.4</td><td></td><td></td><td></td><td>377.85</td><td></td><td>0.9985</td></tr><tr><td>TC-PD070</td><td>39.3012</td><td>175.5193</td><td>1528.0</td><td>—</td><td>—</td><td>0.997</td><td>1.99</td><td>4.3</td><td></td><td></td><td></td><td>376.62</td><td></td><td>0.9985</td></tr><tr><td>WT-PD073</td><td>39.2569</td><td>175.5428</td><td>1892.4</td><td>20</td><td>260</td><td>0.987</td><td>2.09</td><td>3.5</td><td>493</td><td>0.013</td><td>641</td><td>468.33</td><td>0.77</td><td>0.9984</td></tr><tr><td>WT-PD074</td><td>39.2569</td><td>175.5428</td><td>1892.1</td><td>—</td><td>—</td><td>0.988</td><td>2.22</td><td>3.8</td><td></td><td></td><td></td><td>471.99</td><td></td><td>0.9984</td></tr><tr><td>WT-PD075</td><td>39.2560</td><td>175.5397</td><td>1891.2</td><td>—</td><td>—</td><td>0.997</td><td>2.18</td><td>3.3</td><td></td><td></td><td></td><td>471.99</td><td></td><td>0.9984</td></tr><tr><td>PR-PD083</td><td>39.2370</td><td>175.5672</td><td>1730.7</td><td>16</td><td>310</td><td>0.979</td><td>2.28</td><td>3.5</td><td>1557</td><td>0.020</td><td>3115</td><td>453.66</td><td>2.11</td><td>0.9965</td></tr><tr><td>PR-PD084</td><td>39.2386</td><td>175.5689</td><td>1860.9</td><td>24</td><td>180</td><td>0.988</td><td>2.18</td><td>4.7</td><td></td><td></td><td></td><td>492.79</td><td></td><td>0.9968</td></tr><tr><td>PR-PD085</td><td>39.2385</td><td>175.5688</td><td>1857.9</td><td>16</td><td>330</td><td>0.997</td><td>2.12</td><td>4.3</td><td></td><td></td><td></td><td>497.68</td><td></td><td>0.9968</td></tr><tr><td>TFa-PD088</td><td>39.2067</td><td>175.5668</td><td>1308.2</td><td>—</td><td>—</td><td>0.999</td><td>2.39</td><td>3.8</td><td>513</td><td>0.015</td><td>769</td><td>321.60</td><td>1.17</td><td>0.9973</td></tr><tr><td>TFa-PD090</td><td>39.2060</td><td>175.5665</td><td>1290.4</td><td>16</td><td>90</td><td>0.996</td><td>2.31</td><td>5.0</td><td></td><td></td><td></td><td>316.71</td><td></td><td>0.9972</td></tr><tr><td>TFa-PD091</td><td>39.2059</td><td>175.5664</td><td>1288.2</td><td>—</td><td>—</td><td>0.999</td><td>2.23</td><td>4.2</td><td></td><td></td><td></td><td>317.93</td><td></td><td>0.9972</td></tr><tr><td>SC-PD093</td><td>39.2115</td><td>175.6139</td><td>1308.2</td><td>17</td><td>40</td><td>0.993</td><td>2.30</td><td>2.6</td><td>537</td><td>0.010</td><td>537</td><td>313.04</td><td>3.71</td><td>0.9911</td></tr><tr><td>TSa-PD205</td><td>39.2761</td><td>175.6021</td><td>1905.0</td><td>—</td><td>—</td><td>0.983</td><td>2.12</td><td>5.6</td><td>305</td><td>0.010</td><td>305</td><td>476.89</td><td>1.05</td><td>0.9983</td></tr><tr><td>TSa-PD206</td><td>39.2761</td><td>175.6021</td><td>1905.9</td><td>—</td><td>—</td><td>0.997</td><td>2.21</td><td>4.8</td><td></td><td></td><td></td><td>480.56</td><td></td><td>0.9984</td></tr><tr><td>TSa-PD207</td><td>39.2761</td><td>175.6021</td><td>1905.5</td><td>—</td><td>—</td><td>0.997</td><td>2.37</td><td>4.2</td><td></td><td></td><td></td><td>481.78</td><td></td><td>0.9984</td></tr><tr><td>Tsb-PD209</td><td>39.2815</td><td>175.5993</td><td>1932.5</td><td>—</td><td>—</td><td>0.997</td><td>2.20*</td><td>2.4</td><td>506</td><td>0.010</td><td>506</td><td>494.01</td><td>0.75</td><td>0.9989</td></tr><tr><td>Tsb-PD210</td><td>39.2815</td><td>175.5992</td><td>1935.0</td><td>17</td><td>90</td><td>0.989</td><td>2.14</td><td>2.9</td><td></td><td></td><td></td><td>500.13</td><td></td><td>0.9988</td></tr><tr><td>Tsb-PD211</td><td>39.2816</td><td>175.5993</td><td>1929.2</td><td>20</td><td>110</td><td>0.993</td><td>2.26</td><td>3.7</td><td></td><td></td><td></td><td>497.68</td><td></td><td>0.9988</td></tr><tr><td>Tft-PD212</td><td>39.2726</td><td>175.6261</td><td>1521.2</td><td>—</td><td>—</td><td>0.994</td><td>2.07</td><td>5.9</td><td>703</td><td>0.010</td><td>703</td><td>359.50</td><td>0.76</td><td>0.9984</td></tr><tr><td>Tft-PD213</td><td>39.2726</td><td>175.6263</td><td>1522.0</td><td>10</td><td>40</td><td>0.998</td><td>2.14</td><td>6.5</td><td></td><td></td><td></td><td>369.29</td><td></td><td>0.9985</td></tr><tr><td>Tft-PD214</td><td>39.2723</td><td>175.6271</td><td>1506.4</td><td>—</td><td>—</td><td>0.988</td><td>2.23</td><td>6.5</td><td></td><td></td><td></td><td>353.39</td><td></td><td>0.9984</td></tr><tr><td>MN-PD217</td><td>39.2829</td><td>175.5322</td><td>1815.9</td><td>—</td><td>—</td><td>0.993</td><td>2.11</td><td>2.6</td><td>623</td><td>0.008</td><td>498</td><td>446.32</td><td>0.65</td><td>0.9989</td></tr><tr><td>MN-PD218</td><td>39.2829</td><td>175.5321</td><td>1813.9</td><td>13</td><td>220</td><td>0.993</td><td>2.24</td><td>5.2</td><td></td><td></td><td></td><td>448.77</td><td></td><td>0.9989</td></tr><tr><td>MN-PD219</td><td>39.2829</td><td>175.5321</td><td>1812.1</td><td>—</td><td>—</td><td>0.993</td><td>2.22</td><td>4.0</td><td></td><td></td><td></td><td>425.53</td><td></td><td>0.9988</td></tr><tr><td>MN-PD220</td><td>39.2829</td><td>175.5322</td><td>1817.5</td><td>—</td><td>—</td><td>0.993</td><td>2.06</td><td>3.2</td><td></td><td></td><td></td><td>449.99</td><td></td><td>0.9989</td></tr><tr><td>MN-PD221</td><td>39.2829</td><td>175.5325</td><td>1822.8</td><td>—</td><td>—</td><td>0.993</td><td>2.20</td><td>4.0</td><td></td><td></td><td></td><td>446.32</td><td></td><td>0.9989</td></tr><tr><td>MS-PD222</td><td>39.2845</td><td>175.5304</td><td>1750.6</td><td>36</td><td>170</td><td>0.955</td><td>2.23</td><td>5.0</td><td>901</td><td>0.010</td><td>901</td><td>414.53</td><td>0.76</td><td>0.9986</td></tr><tr><td>MS-PD223</td><td>39.2845</td><td>175.5305</td><td>1751.4</td><td>—</td><td>—</td><td>0.992</td><td>2.41</td><td>4.0</td><td></td><td></td><td></td><td>412.08</td><td></td><td>0.9986</td></tr><tr><td>MS-PD224</td><td>39.2845</td><td>175.53005</td><td>1750.9</td><td>—</td><td>—</td><td>0.992</td><td>2.33</td><td>4.1</td><td></td><td></td><td></td><td>401.08</td><td></td><td>0.9986</td></tr><tr><td>CTa-PD229</td><td>39.2958</td><td>175.5395</td><td>1924.0</td><td>7</td><td>300</td><td>0.996</td><td>1.96</td><td>3.2</td><td>415</td><td>0.015</td><td>623</td><td>498.90</td><td>1.81</td><td>0.9972</td></tr><tr><td>CTa-PD230</td><td>39.2959</td><td>175.5396</td><td>1925.1</td><td>—</td><td>—</td><td>0.996</td><td>2.21</td><td>2.9</td><td></td><td></td><td></td><td>500.13</td><td></td><td>0.9973</td></tr><tr><td>CTb-PD231</td><td>39.2998</td><td>175.5392</td><td>1877.5</td><td>—</td><td>—</td><td>0.996</td><td>2.22</td><td>3.4</td><td>439</td><td>0.008</td><td>351</td><td>432.87</td><td>1.26</td><td>0.9978</td></tr><tr><td>CTb-PD232</td><td>39.3001</td><td>175.5390</td><td>1873.2</td><td>20</td><td>190</td><td>0.991</td><td>2.14</td><td>4.1</td><td></td><td></td><td></td><td>467.11</td><td></td><td>0.9979</td></tr><tr><td>CTb-PD233</td><td>39.3001</td><td>175.5390</td><td>1872.0</td><td>15</td><td>240</td><td>0.994</td><td>2.17*</td><td>2.8</td><td></td><td></td><td></td><td>465.89</td><td></td><td>0.9979</td></tr><tr><td>CTb-PD234</td><td>39.3003</td><td>175.5391</td><td>1873.4</td><td>—</td><td>—</td><td>0.996</td><td>2.15</td><td>3.2</td><td></td><td></td><td></td><td>467.11</td><td></td><td>0.9979</td></tr><tr><td>LC-PD254</td><td>39.2718</td><td>175.6052</td><td>1827.1</td><td>—</td><td>—</td><td>0.997</td><td>2.01</td><td>5.4</td><td>506</td><td>0.010</td><td>506</td><td>462.22</td><td>1.18</td><td>0.9981</td></tr><tr><td>LC-PD255</td><td>39.2718</td><td>175.6053</td><td>1826.6</td><td>—</td><td>—</td><td>0.997</td><td>2.07</td><td>6.4</td><td></td><td></td><td></td><td>461.00</td><td></td><td>0.9981</td></tr><tr><td>LC-PD256</td><td>39.2718</td><td>175.6053</td><td>1825.6</td><td>—</td><td>—</td><td>0.996</td><td>2.08</td><td>6.0</td><td></td><td></td><td></td><td>452.44</td><td></td><td>0.9980</td></tr><tr><td>LC-PD257</td><td>39.2718</td><td>175.6053</td><td>1824.7</td><td>16</td><td>330</td><td>0.996</td><td>2.05</td><td>3.7</td><td></td><td></td><td></td><td>462.22</td><td></td><td>0.9981</td></tr><tr><td>WG-PD325</td><td>39.2557</td><td>175.5551</td><td>2079.1</td><td>21</td><td>357</td><td>0.991</td><td>2.25</td><td>4.0</td><td>329</td><td>0.008</td><td>264</td><td>536.81</td><td>2.45</td><td>0.9966</td></tr><tr><td>WG-PD326</td><td>39.2556</td><td>175.5549</td><td>2066.7</td><td>—</td><td>—</td><td>0.995</td><td>2.30</td><td>3.2</td><td></td><td></td><td></td><td>520.91</td><td></td><td>0.9965</td></tr><tr><td>DC-PD327</td><td>39.2346</td><td>175.5515</td><td>1600.4</td><td>—</td><td>—</td><td>0.999</td><td>2.22</td><td>3.5</td><td>401</td><td>0.008</td><td>321</td><td>379.07</td><td>0.67</td><td>0.9987</td></tr><tr><td>DC-PD329</td><td>39.2342</td><td>175.5509</td><td>1591.8</td><td>—</td><td>—</td><td>0.999</td><td>2.37</td><td>4.3</td><td></td><td></td><td></td><td>380.29</td><td></td><td>0.9987</td></tr><tr><td>DC-PD330</td><td>39.2341</td><td>175.5507</td><td>1590.3</td><td>—</td><td>—</td><td>0.999</td><td>2.21</td><td>3.4</td><td></td><td></td><td></td><td>377.85</td><td></td><td>0.9987</td></tr></tbody></table>

Density measures were obtained with the hydrostatic method. Density values marked with * were calculated by averaging the densities of other samples from the same site. $P_{nuc}$, Closure age, He$_{nuc}$, and $P_4$ values are considered equal for all samples of the same flow.
**Data availability.** All used data are available in the Supplement file S4 and Appendix Table A2.

**Supplement.** The supplement related to this article is available online at: https://doi.org/10.5194/gchron-6-365-2024-supplement.

**Author contributions.** PD carried out field sampling, mineral separation, He isotope measurements, data processing and interpretation, and manuscript writing. SRE assisted with sampling, data processing, and manuscript revision. BMK handled the project supervision, obtained resources, and reviewed the manuscript. PHB helped with methodology, data analysis, and manuscript revision. ARLN reviewed and edited the manuscript. GSL helped with resources and data interpretation. DBT assisted with data interpretation. JWC helped with manuscript revision. CEC helped with data interpretation. SB assisted with mineral separation. GF, LZ, and BT helped with He isotope measurements.

**Competing interests.** The contact author has declared that none of the authors has any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.

**Acknowledgements.** We would like to acknowledge Ngāti Rangi, Uenuku, and Ngāti Tūwharetoa iwi; tangata whenua and guardians of Ruapehu. We are grateful to Mark Kurz, David Marchetti, and Eric Portenga, whose comments as referees significantly improved the quality of this paper. We would like to thank New Zealand's Department of Conservation for sampling permission. We would also like to thank Amy Dreyer, Gilles Seropian, and Alexander James Marshall for their assistance in the field; Hollei Gabrielsen for her advice on Māori subjects and aid with the permit process; Chris Grimshaw for help with laboratory procedures; and Mark Henson and Nigel Seebeck for their help in the Tukino Access Road ford.

**Financial support.** This research has been supported by the Resilience to Nature's Challenges (RNC) volcano program (grant no. GNS-RNC047), the SAAFE scholarship 2022 from the Australian Institute of Nuclear Science and Engineering (AINSE), the Mason Trust (SEE, University of Canterbury), the Wellman Research Award 2022 (Royal Society of New Zealand), and the Tongariro Memorial Award 2022 (the Tongariro Natural History Society).

**Review statement.** This paper was edited by Marissa Tremblay and reviewed by Mark Kurz, David Marchetti, and Eric Portenga.

## References

Ackert, R. P., Singer, B. S., Guillou, H., Kaplan, M. R., and Kurz, M. D.: Long-term cosmogenic ³He production rates from ⁴⁰Ar/³⁹Ar and K-Ar dated Patagonian lava flows at 47° S, Earth Planet. Sc. Lett., 210, 119–136, https://doi.org/10.1016/S0012-821X(03)00134-1, 2003.

Alcalá-Reygosa, J., Palacios, D., Schimmelpfennig, I., Vázquez-Selem, L., García-Sancho, L., Franco-Ramos, O., Villanueva, J., Zamorano, J. J., Aumaître, G., Bourlès, D., and Keddadouche, K.: Dating late Holocene lava flows in Pico de Orizaba (Mexico) by means of in situ-produced cosmogenic ³⁶Cl, lichenometry and dendrochronology, Quat. Geochronol., 47, 93–106, https://doi.org/10.1016/j.quageo.2018.05.011, 2018.

Anderson, S. W., Krinsley, D. H., and Fink, J. H.: Criteria for recognition of constructional silicic lava flow surfaces, Earth Surf. Proc. Land., 19, 531–541, https://doi.org/10.1002/esp.3290190606, 1994.

Balco, G., Stone, J. O., Lifton, N. A., and Dunai, T. J.: A complete and easily accessible means of calculating surface exposure ages or erosion rates from ¹⁰Be and ²⁶Al measurements, Quat. Geochronol., 3, 174–195, https://doi.org/10.1016/j.quageo.2007.12.001, 2008.

Barrell, D. J.: Quaternary Glaciers of New Zealand, in: Developments in Quaternary Sciences, 15 edn., chap. 75, Elsevier, 1047–1064, https://doi.org/10.1016/B978-0-444-53447-7.00075-1, 2011.

Barrell, D. J. A., Almond, P. C., Vandergoes, M. J., Lowe, D. J., and Newnham, R. M.: A composite pollen-based stratotype for inter-regional evaluation of climatic events in New Zealand over the past 30,000 years (NZ-INTIMATE project), Quaternary Sci. Rev., 74, 4–20, https://doi.org/10.1016/j.quascirev.2013.04.002, 2013.

Blard, P.-H.: Cosmogenic ³He in terrestrial rocks: A review, Chem. Geol., 586, 120543, https://doi.org/10.1016/j.chemgeo.2021.120543, 2021.

Blard, P.-H. and Farley, K. A.: The influence of radiogenic ⁴He on cosmogenic ³He determinations in volcanic olivine and pyroxene, Earth Planet. Sc. Lett., 276, 20–29, https://doi.org/10.1016/j.epsl.2008.09.003, 2008.

Blard, P. H., Pik, R., Lavé, J., Bourlès, D., Burnard, P. G., Yokochi, R., Marty, B., and Trusdell, F.: Cosmogenic ³He production rates revisited from evidences of grain size dependent release of matrix-sited helium, Earth Planet. Sc. Lett., 247, 222–234, https://doi.org/10.1016/j.epsl.2006.05.012, 2006.

Blard, P. H., Lavé, J., Pik, R., Wagnon, P., and Bourlès, D.: Persistence of full glacial conditions in the central Pacific until 15,000 years ago, Nature, 449, 591–594, https://doi.org/10.1038/nature06142, 2007.

Blard, P.-H., Braucher, R., Lavé, J., and Bourlès, D.: Cosmogenic ¹⁰Be production rate calibrated against ³He in the high Tropical Andes (3800–4900 m, 20–22° S), Earth Planet. Sc. Lett., 382, 140–149, https://doi.org/10.1016/j.epsl.2013.09.010, 2013.

Blard, P.-H., Balco, G., Burnard, P. G., Farley, K. A., Fenton, C. R., Friedrich, R., Jull, A. J., Niedermann, S., Pik, R., Schaefer, J. M., Scott, E. M., Shuster, D. L., Stuart, F. M., Stute, M., Tibari, B., Winckler, G., and Zimmermann, L.: An interlaboratory comparison of cosmogenic ³He and radiogenic ⁴He in
the CRONUS-P pyroxene standard, Quat. Geochronol., 26, 11–19, https://doi.org/10.1016/j.quageo.2014.08.004, 2015.

Bromley, G. R., Winckler, G., Schaefer, J. M., Kaplan, M. R., Licht, K. J., and Hall, B. L.: Pyroxene separation by HF leaching and its impact on helium surface-exposure dating, Quat. Geochronol., 23, 1–8, https://doi.org/10.1016/j.quageo.2014.04.003, 2014.

Buchanan-Banks, J. M., Lockwood, J. P., and Rubin, M.: Radiocarbon dates for lava flows from northeast rift zone of Mauna Loa volcano, Hilo 7 1/2' Quadrangle, island of Hawaii, Radiocarbon, 31, 179–186, 1989.

Calvert, A. T. and Lanphere, M. A.: Argon geochronology of Kilauea's early submarine history, J. Volcanol. Geoth. Res., 151, 1–18, https://doi.org/10.1016/j.jvolgeores.2005.07.023, 2006.

Calvert, A. T., Fierstein, J., and Hildreth, W.: Eruptive history of Middle Sister, Oregon Cascades, USA-Product of a late Pleistocene eruptive episode, Geosphere, 14, 2118–2139, https://doi.org/10.1130/GES01638.1, 2018.

Carignan, J., Hild, P., Mevelle, G., Morel, J., and Yeghicheyan, D.: Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: A study of geochemical reference materials BR, DR-N, UB-N, AN-G and GH, Geostandard. Newslett., 25, 187–198, https://doi.org/10.1111/j.1751-908x.2001.tb00595.x, 2001.

Cerling, T. E. and Craig, H.: Geomorphology and in-situ cosmogenic isotopes, Annu. Rev. Earth Pl. Sc., 273–317, 1994.

Clay, P. L., Busemann, H., Sherlock, S. C., Barry, T. L., Kelley, S. P., and McGarvie, D. W.: $^{40}$Ar/$^{39}$Ar ages and residual volatile contents in degassed subaerial and subglacial glassy volcanic rocks from Iceland, Chem. Geol., 403, 99–110, https://doi.org/10.1016/j.chemgeo.2015.02.041, 2015.

Coble, M. A., Grove, M., and Calvert, A. T.: Calibration of Nu-Instruments Noblesse multicollector mass spectrometers for argon isotopic measurements using a newly developed reference gas, Chem. Geol., 290, 75–87, https://doi.org/10.1016/j.chemgeo.2011.09.003, 2011.

Cole, J. W. and Lewis, K. B.: Evolution of the Taupo-Hikurangi subduction system, Tectonophysics, 72, 1–21, https://doi.org/10.1016/0040-1951(81)90084-6, 1981.

Connor, C., Bebbington, M., and Marzocchi, W.: Probabilistic Volcanic Hazard Assessment, Elsevier Inc., second edn., ISBN 9780123859389, https://doi.org/10.1016/b978-0-12-385938-9.00051-1, 2015.

Conway, C. E.: Studies on the Glaciovolcanic and Magmatic Evolution of Ruapehu Volcano, New Zealand, PhD thesis, Victoria University of Wellington, https://researcharchive.vuw.ac.nz/handle/10063/5152 (last access: 25 May 2024), 2016.

Conway, C. E., Townsend, D. B., Leonard, G. S., Wilson, C. J., Calvert, A. T., and Gamble, J. A.: Lava-ice interaction on a large composite volcano: a case study from Ruapehu, New Zealand, B. Volcanol., 77, 21, https://doi.org/10.1007/s00445-015-0906-2, 2015.

Conway, C. E., Leonard, G. S., Townsend, D. B., Calvert, A. T., Wilson, C. J., Gamble, J. A., Eaves, S. R., Pure, L. R., Leonard, G. S., Townsend, D. B., Wilson, C. J., Calvert, A. T., Cole, R. P., Conway, C. E., Gamble, J. A., and Smith, T. B.: A high-resolution $^{40}$Ar/$^{39}$Ar lava chronology and edifice construction history for Ruapehu volcano, New Zealand, J. Volcanol. Geoth. Res.,

327, 152–179, https://doi.org/10.1016/j.jvolgeores.2016.07.006, 2016.

Delunel, R., Blard, P.-H., Martin, L. C., Nomade, S., and Schlünegger, F.: Long term low latitude and high elevation cosmogenic $^{3}$He production rate inferred from a 107 ka-old lava flow in northern Chile; 22° S–3400 m a.s.l., Geochim. Cosmochim. Ac., 184, 71–87, https://doi.org/10.1016/j.gca.2016.04.023, 2016.

Donoghue, S. L.: Late quaternary volcanic stratigraphy of the southeastern sector of the Mount Ruapehu ring plain New Zealand, PhD thesis, Massey University, https://mro.massey.ac.nz/items/516a0d80-eda3-4a7e-a495-2e13fcb7821c (last access: 25 May 2024), 1991.

Donoghue, S. L. and Neall, V. E.: Late Quaternary constructional history of the southeastern Ruapehu ring plain, New Zealand, New Zealand. J. Geol. Geop., 44, 439–466, https://doi.org/10.1080/00288306.2001.9514949, 2001.

Donoghue, S. L., Neall, V. E., and Palmer, A. S.: Stratigraphy and chronology of late quaternary andesitic tephra deposits, Tongariro Volcanic Centre, New Zealand, J. Roy. Soc. New Zeal., 25, 115–206, https://doi.org/10.1080/03014223.1995.9517487, 1995.

Donoghue, S. L., Stewart, R. B., Neall, V. E., Lecointre, J., Price, R., Palmer, A. S., McClelland, E., and Hobson, K.: The Taurewa Eruptive Episode: Evidence for climactic eruptions at Ruapehu volcano, New Zealand, B. Volcanol., 61, 223–240, https://doi.org/10.1007/s004450050273, 1999.

Donoghue, S. L., Vallance, J. W., Smith, I. E., and Stewart, R. B.: Using geochemistry as a tool for correlating proximal andesitic tephra: case studies from Mt Rainier (USA) and MT Ruapehu (New Zealand), J. Quaternary Sci., 22, 395–410, https://doi.org/10.1002/jqs.1065, 2007.

Dostal, J., Dupuy, C., Carron, J. P., Le Guen de Kerneizon, M., and Maury, R. C.: Partition coefficients of trace elements: Application to volcanic rocks of St. Vincent, West Indies, Geochim. Cosmochim. Ac., 47, 525–533, https://doi.org/10.1016/0016-7037(83)90275-2, 1983.

Dunai, T. J.: Cosmogenic Nuclides. Principles, Concepts and Application in the Earth Surface Sciences, Cambridge University Press, ISBN 9780521873802, 2010.

Dunn, T. and Sen, C.: Mineral/matrix partition coefficients for orthopyroxene, plagioclase, and olivine in basaltic to andesitic systems: A combined analytical and experimental study, Geochim. Cosmochim. Ac., 58, 717–733, https://doi.org/10.1016/0016-7037(94)90501-0, 1994.

Eaves, S. R. and Brook, M. S.: Glaciers and glaciation of North Island, New Zealand, New Zealand. J. Geol. Geop., 64, 1–20, https://doi.org/10.1080/00288306.2020.1811354, 2021.

Eaves, S. R., Winckler, G., Schaefer, J. J. M., Vandergoes, M. J., Alloway, B. V., Mackintosh, A. N., Townsend, D. B., Ryan, M. T., and Li, X.: A test of the cosmogenic $^{3}$He production rate in the south-west Pacific (39° S), J. Quaternary Sci., 30, 79–87, https://doi.org/10.1002/jqs.2760, 2015.

Eaves, S. R., Mackintosh, A. N., Anderson, B. M., Doughty, A. M., Townsend, D. B., Conway, C. E., Winckler, G., Schaefer, J. M., Leonard, G. S., and Calvert, A. T.: The Last Glacial Maximum in the central North Island, New Zealand: palaeoclimate inferences from glacier modelling, Clim. Past, 12, 943–960, https://doi.org/10.5194/cp-12-943-2016, 2016a.
Eaves, S. R., Mackintosh, A. N., Winckler, G., Schaefer, J. M., Alloway, B. V., and Townsend, D. B.: A cosmogenic ³He chronology of late Quaternary glacier fluctuations in North Island, New Zealand (39°S), Quaternary Sci. Rev., 132, 40–56, https://doi.org/10.1016/j.quascirev.2015.11.004, 2016b.

Eaves, S. R., Winckler, G., Mackintosh, A. N., Schaefer, J. M., Townsend, D. B., Doughty, A. M., Jones, R. S., and Leonard, G. S.: Late-glacial and Holocene glacier fluctuations in North Island, New Zealand, Quaternary Sci. Rev., 223, 105914, https://doi.org/10.1016/j.quascirev.2019.105914, 2019.

Espanon, V. R., Honda, M., and Chivas, A. R.: Cosmogenic ³He and ²¹Ne surface exposure dating of young basalts from Southern Mendoza, Argentina, Quat. Geochronol., 19, 76–86, https://doi.org/10.1016/j.quageo.2013.09.002, 2014.

Fenton, C. R. and Niedermann, S.: Surface exposure dating of young basalts (1-200ka) in the San Francisco volcanic field (Arizona, USA) using cosmogenic ³He and ²¹Ne, Quat. Geochronol., 19, 87–105, https://doi.org/10.1016/j.quageo.2012.10.003, 2014.

Fenton, C. R., Webb, R. H., Pearthree, P. A., Cerling, T. E., and Poreda, R. J.: Displacement rates on the Toroweap and Hurricane faults: Implications for Quaternary downcutting in the Grand Canyon, Arizona, Geology, 29, 1035–1038, https://doi.org/10.1130/0091-7613(2001)029<1035:DROTTA>2.0.CO;2, 2001.

Fenton, C. R., Niedermann, S., Goethals, M. M., Schneider, B., and Wijbrans, J.: Evaluation of cosmogenic ³He and ²¹Ne production rates in olivine and pyroxene from two Pleistocene basalt flows, western Grand Canyon, AZ, USA, Quat. Geochronol., 4, 475–492, https://doi.org/10.1016/j.quageo.2009.08.002, 2009.

Ferrier, K. L., Taylor Perron, J., Mukhopadhyay, S., Rosener, M., Stock, J. D., Huppert, K. L., and Slosberg, M.: Covariation of climate and long-term erosion rates across a steep rainfall gradient on the Hawaiian island of Kaua'i, Bull. Geol. Soc. Am., 125, 1146–1163, https://doi.org/10.1130/B30726.1, 2013.

Fierstein, J., Hildreth, W., and Calvert, A. T.: Eruptive history of South Sister, Oregon Cascades, J. Volcanol. Geoth. Res., 207, 145–179, https://doi.org/10.1016/j.jvolgeores.2011.06.003, 2011.

Fleck, R. J., Hagstrum, J. T., Calvert, A. T., Evarts, R. C., and Conrey, R. M.: ⁴⁰Ar/³⁹Ar geochronology, paleomagnetism, and evolution of the Boring volcanic field, Oregon and Washington, USA, Geosphere, 1283–1314, https://doi.org/10.1130/GES00985.1, 2014.

Foeken, J. P., Day, S., and Stuart, F. M.: Cosmogenic ³He exposure dating of the Quaternary basalts from Fogo, Cape Verde: Implications for rift zone and magnetic reorganisation, Quat. Geochronol., 4, 37–49, https://doi.org/10.1016/j.quageo.2008.07.002, 2009.

Gabrielsen, H., Procter, J., Rainforth, H., Black, T., Harmsworth, G., and Pardo, N.: Reflections from an Indigenous Community on Volcanic Event Management, Communications and Resilience, in: Observing the Volcano World. Advances in Volcanology, edited by: Fearnley, C. J., Bird, D. K., Haynes, K., McGuire, W. J., and Jolly, G., Springer, Cham. 463–479, https://doi.org/10.1007/11157_2016_44, 2018.

Gallahan, W. E. and Nielsen, R. L.: The partitioning of Sc, Y, and the rare earth elements between high-Ca pyroxene and natural mafic to intermediate lavas at 1 atmosphere, Geochim.

Cosmochim. Ac., 56, 2387–2404, https://doi.org/10.1016/0016-7037(92)90196-P, 1992.

Gamble, J. A., Price, R. C., Smith, I. E., McIntosh, W. C., and Dunbar, N. W.: ⁴⁰Ar/³⁹Ar geochronology of magmatic activity, magma flux and hazards at Ruapehu volcano, Taupo Volcanic Zone, New Zealand, J. Volcanol. Geoth. Res., 120, 271–287, https://doi.org/10.1016/S0377-0273(02)00407-9, 2003.

Gosse, J. C. and Phillips, F. M.: Terrestrial in situ cosmogenic nuclides: Theory and application, Quaternary Sci. Rev., 20, 1475–1560, https://doi.org/10.1016/S0277-3791(00)00171-2, 2001.

Greve, A., Turner, G. M., Conway, C. E., Townsend, D. B., Gamble, J. A., and Leonard, G. S.: Palaeomagnetic refinement of the eruption ages of Holocene lava flows, and implications for the eruptive history of the Tongariro Volcanic Centre, New Zealand, Geophys. J. Int., 207, 702–718, https://doi.org/10.1093/gji/ggw296, 2016.

Hackett, W. R.: Geology and petrology of Ruapehu volcano and related vents, PhD thesis, Victoria University of Wellington, http://researcharchive.vuw.ac.nz/handle/10063/743 (last access: 25 May 2024), 1985.

Harpel, C. J., Kyle, P. R., Esser, R. P., McIntosh, W. C., and Caldwell, D. A.: ⁴⁰Ar/³⁹Ar dating of the eruptive history of Mount Erebus, Antarctica: Summit flows, tephra, and caldera collapse, B. Volcanol., 66, 687–702, https://doi.org/10.1007/s00445-004-0349-7, 2004.

Harris, A. J. L.: Basaltic Lava Flow Hazard, in: Volcanic Hazards, Risks and Disasters, edited by: Schroder, J. F. and Papale, P., chap. 2, 17–46, Elsevier, https://doi.org/10.1016/C2011-0-07012-6, 2015.

Hilton, D. R., Fischer, T. P., and Marry, B.: Noble gases and volatile recycling at subduction zones, Rev. Mineral. Geochem., 47, 319–370, https://doi.org/10.2138/rmg.2002.47.9, 2002.

Houghton, B. F., Wilson, C. J., McWilliams, M. O., Laphere, M. A., Weaver, S. D., Briggs, R. M., and Pringle, M. S.: Chronology and dynamics of a large silicic magmatic system: central Taupo Volcanic Zone, New Zealand, Geology, 23, 13–16, https://doi.org/10.1130/0091-7613(1995)023<0013:CADOAL>2.3.CO;2, 1995.

Jenkins, S. F., Day, S. J., Faria, B. V., and Fonseca, J. F.: Damage from lava flows: insights from the 2014–2015 eruption of Fogo, Cape Verde, Journal of Applied Volcanology, 6, 6, https://doi.org/10.1186/s13617-017-0057-6, 2017.

Keys, H. and Green, P. M.: Mitigation of volcanic risks at Mt Ruapehu, New Zealand, Proceedings of the Moutain Risks international conference, Firenze, Italy, 24–26 November, 485–490, https://www.researchgate.net/profile/Harry-Keys/publication/281606085_Mitigation_of_volcanic_risks_at_Mt_Ruapehu_New_Zealand/links/581f8fa608aeccc08af3abd6/Mitigation-of-volcanic-risks-at-Mt-Ruapehu-New-Zealand.pdf (last access: 25 May 2024), 2010.

Klein, J., Giegengack, R., Middleton, R., Sharma, P., Underwood, J. R., and Weeks, R. A.: Revealing Histories of Exposure Using in Situ Produced ²⁶Al and ¹⁰Be in Libyan Desert Glass, Radiocarbon, 28, 547–555, 1986.

Kurz, M. D.: Cosmogenic helium in a terrestrial igneous rock, Nature, 320, 435–439, https://doi.org/10.1038/320435a0, 1986a.

Kurz, M. D.: In situ production of terrestrial cosmogenic helium and some applications to geochronology, Geochim. Cosmochim. Ac.,
50, 2855-2862, https://doi.org/10.1016/0016-7037(86)90232-2, 1986b.

Kurz, M. D., Colodner, D., Trull, T. W., Moore, R. B., and O'Brien, K.: Cosmic ray exposure dating with in situ produced cosmogenic $^{3}$He: results from young Hawaiian lava flows, Earth Planet. Sc. Lett., 97, 177-189, https://doi.org/10.1016/0012-821X(90)90107-9, 1990.

Lal, D.: An important source of $^{4}$He (and $^{3}$He) in diamonds, Earth Planet. Sc. Lett., 96, 1-7, https://doi.org/10.1016/0012-821X(89)90118-0, 1989.

Lal, D.: Cosmic ray labeling of erosion surfaces: in situ nuclide production rates and erosion models, Earth Planet. Sc. Lett., 104, 424-439, https://doi.org/10.1016/0012-821X(91)90220-C, 1991.

Lanphere, M. A.: Comparison of conventional K-Ar and $^{40}$Ar/$^{39}$Ar dating of young mafic volcanic rocks, Quaternary Res., 53, 294-301, https://doi.org/10.1006/qres.1999.2122, 2000.

Le Maitre, R. W.: Igneous Rocks: A Classification and Glossary of Terms, Recommendations of the International Union of Geological Sciences, Subcomission of the Systematics of Igneous Rocks, Cambridge University Press, https://doi.org/10.1017/CBO9780511535581, 2002.

Leonard, G. S., Cole, R. P. R., Christenson, B. W., Conway, C. E., Cronin, S. J., Gamble, J. A., Hurst, T., Kennedy, B. M., Miller, C. A., Procter, J. N., Pure, L. R., Townsend, D. B., White, J. D., and Wilson, C. J.: Ruapehu and Tongariro stratovolcanoes: a review of current understanding, New Zeal. J. Geol. Geop., 64, 389-420, https://doi.org/10.1080/00288306.2021.1909080, 2021.

Leya, I., Lange, H. J., Neumann, S., Wieler, R., and Michel, R.: The production of cosmogenic nuclides in stony meteoroids by galactic cosmic-ray particles, Meteorit. Planet. Sci., 35, 259-286, https://doi.org/10.1111/j.1945-5100.2000.tb01775.x, 2000.

Licciardi, J. M., Kurz, M. D., and Curtice, J. M.: Cosmogenic $^{3}$He production rates from Holocene lava flows in Iceland, Earth Planet. Sc. Lett., 246, 251-264, https://doi.org/10.1016/j.epsl.2006.03.016, 2006.

Licciardi, J. M., Kurz, M. D., and Curtice, J. M.: Glacial and volcanic history of Icelandic table mountains from cosmogenic $^{3}$He exposure ages, Quaternary Sci. Rev., 26, 1529-1546, https://doi.org/10.1016/j.quascirev.2007.02.016, 2007.

Lifton, N.: Implications of two Holocene time-dependent geomagnetic models for cosmogenic nuclide production rate scaling, Earth Planet. Sc. Lett., 433, 257-268, https://doi.org/10.1016/j.epsl.2015.11.006, 2016.

Lifton, N., Sato, T., and Dunai, T. J.: Scaling in situ cosmogenic nuclide production rates using analytical approximations to atmospheric cosmic-ray fluxes, Earth Planet. Sc. Lett., 386, 149-160, https://doi.org/10.1016/j.epsl.2013.10.052, 2014.

Lippolt, H. J. and Weigel, E.: $^{4}$He diffusion in $^{40}$Ar-retentive minerals, Geochim. Cosmochim. Ac., 52, 1449-1458, https://doi.org/10.1016/0016-7037(88)90215-3, 1988.

Lockwood, J. P. and Lipman, P. W.: Recovery of datable charcoal beneath young lavas: Lessons from Hawaii, Bulletin Volcanologique, 43, 609-615, https://doi.org/10.1007/BF02597697, 1980.

Luhr, J. F. and Carmichael, I. S. E.: The Colima Volcanic complex, Mexico, Contrib. Mineral. Petr., 71, 343-372, https://doi.org/10.1007/bf00374707, 1980.

Lupton, J. and Evans, L.: Changes in the atmospheric helium isotope ratio over the past 40 years, Geophys. Res. Lett., 40, 6271-6275, https://doi.org/10.1002/2013GL057681, 2013.

Marchetti, D. W., Hynek, S. A., and Cerling, T. E.: Cosmogenic $^{3}$He exposure ages of basalt flows in the northwestern Payún Matru volcanic field, Mendoza Province, Argentina, Quat. Geochronol., 19, 67-75, https://doi.org/10.1016/j.quageo.2012.10.004, 2014.

Martin, L. C., Blard, P.-H., Balco, G., Lavé, J., Delunel, R., Lifton, N., and Laurent, V.: The CREp program and the ICE-D production rate calibration database: A fully parameterizable and updated online tool to compute cosmic-ray exposure ages, Quat. Geochronol., 38, 25-49, https://doi.org/10.1016/j.quageo.2016.11.006, 2017.

Matsuda, J., Matsumoto, T., Sumino, H., Nagao, K., Yamamoto, J., Miura, Y., Kaneoka, I., Takahata, N., and Sano, Y.: The $^{3}$He/$^{4}$He ratio of new internal He Standard of Japan (HESJ), Geochem. J., 36, 191-195, https://doi.org/10.2343/geochemj.36.191, 2002.

Mc Arthur, J. L. and Shepherd, M. J.: Late Quaternary glaciation of Mt. Ruapehu, North Island, New Zealand, J. Roy. Soc. New Zeal., 20, 287-296, https://doi.org/10.1080/03036758.1990.10416823, 1990.

Medynski, S., Pik, R., Burnard, P., Vye-Brown, C., France, L., Schimmelpfennig, I., Whaler, K., Johnson, N., Benedetti, L., Ayelew, D., and Yirgu, G.: Stability of rift axis magma reservoirs: Spatial and temporal evolution of magma supply in the Dabbahu rift segment (Afar, Ethiopia) over the past 30 kyr, Earth Planet. Sc. Lett., 409, 278-289, https://doi.org/10.1016/j.epsl.2014.11.002, 2015.

Mishra, A. K., Placzek, C., Wurster, C., and Whitehead, P. W.: New radiocarbon age constraints for the 120 km-long Toomba flow, north Queensland, Australia, Aust. J. Earth Sci., 66, 71-79, https://doi.org/10.1080/08120099.2019.1523227, 2019.

Moore, R. B. and Rubin, M.: Radiocarbon dates for lava flows and pyroclastic deposits on Sao Miguel, Azores, Radiocarbon, 33, 151-164, https://doi.org/10.1017/S0033822200013278, 1991.

Morimoto, N., Fabries, J., Ferguson, A., Ginzburg, I., Ross, M., Seifert, F., Zussman, J., Aoki, K., and Gottardi, G.: Nomenclature of pyroxenes Subcommittee on Pyroxenes Commission on New Minerals and Mineral Names International Mineralogical Association, Am. Mineral., 73, 1123-1133, 1988.

Muscheler, R., Beer, J., Kubik, P. W., and Synal, H. A.: Geomagnetic field intensity during the last 60,000 years based on $^{10}$Be and $^{36}$Cl from the Summit ice cores and $^{14}$C, Quaternary Sci. Rev., 24, 1849-1860, https://doi.org/10.1016/j.quascirev.2005.01.012, 2005.

Nairn, I. A., Kobayashi, T., and Nakagawa, M.: The ~10 ka multiple vent pyroclastic eruption sequence at Tongariro Volcanic Centre, Taupo Volcanic Zone, New Zealand: Part 1. Eruptive processes during regional extension, J. Volcanol. Geoth. Res., 86, 19-44, https://doi.org/10.1016/S0377-0273(98)00085-7, 1998.

Nicholls, I. A. and Harris, K. L.: Experimental rare earth element partition coefficients for garnet, clinopyroxene and amphibole coexisting with andesitic and basaltic liquids, Geochim. Cosmochim. Ac., 44, 287-308, https://doi.org/10.1016/0016-7037(80)90138-6, 1980.

Niedermann, S.: Cosmic-ray-produced noble gases in terrestrial rocks: Dating tools for surface processes, Rev. Mineral. Geochem., 47, 731-784, https://doi.org/10.2138/rmg.2002.47.16, 2002.
Nishiizumi, K.: Cosmic ray production rates of $^{10}$Be and $^{26}$Al in quartz from glacially polished rocks, J. Geophys. Res., 94, 17907–17915, https://doi.org/10.1029/jb094ib12p17907, 1989.

Nishiizumi, K., Kohl, C. P., Arnold, J. R., Klein, J., Fink, D., and Middleton, R.: Cosmic ray produced $^{10}$Be and $^{26}$Al in Antarctic rocks: exposure and erosion history, Earth Planet. Sc. Lett., 104, 440–454, https://doi.org/10.1016/0012-821X(91)90221-3, 1991.

Palmer, B. A. and Neall, V. E.: The Murimotu Formation – 9500 year old deposits of a debris avalanche and associated lahars, Mount Ruapehu, North Island, New Zealand, New Zealand. J. Geol. Geop., 32, 477–486, https://doi.org/10.1080/00288306.1989.10427555, 1989.

Pardo, N.: Andesitic Plinian Eruptions at Mt. Ruapehu (New Zealand): From Lithofacies to Eruption Dynamics, PhD thesis, Massey University, https://mro.massey.ac.nz/items/c2f61b33-2579-410f-b493-eeab08691375 (last access: 25 May 2024), 2012.

Pardo, N., Cronin, S., Palmer, A., Procter, J., and Smith, I.: Andesitic Plinian eruptions at Mt. Ruapehu: Quantifying the uppermost limits of eruptive parameters, B. Volcanol., 74, 1161–1185, https://doi.org/10.1007/s00445-012-0588-y, 2012a.

Pardo, N., Cronin, S. J., Palmer, A. S., and Németh, K.: Reconstructing the largest explosive eruptions of Mt. Ruapehu, New Zealand: Lithostratigraphic tools to understand subplinian-plinian eruptions at andesitic volcanoes, B. Volcanol., 74, 617–640, https://doi.org/10.1007/s00445-011-0555-z, 2012b.

Parmelee, D. E., Kyle, P. R., Kurz, M. D., Marrero, S. M., and Phillips, F. M.: A new Holocene eruptive history of Erebus volcano, Antarctica using cosmogenic $^3$He and $^{36}$Cl exposure ages, Quat. Geochronol., 30, 114–131, https://doi.org/10.1016/j.quageo.2015.09.001, 2015.

Patterson, D. B., Honda, M., and McDougall, I.: Noble gases in mafic phenocrysts and xenoliths from New Zealand, Geochim. Cosmochim. Ac., 58, 4411–4427, https://doi.org/10.1016/0016-7037(94)90344-1, 1994.

Preece, K., Mark, D. F., Barclay, J., Cohen, B. E., Chamberlain, K. J., Jowitt, C., Vye-Brown, C., Brown, R. J., and Hamilton, S.: Bridging the gap: $^{40}$Ar/$^{39}$Ar dating of volcanic eruptions from the “Age of Discovery”, Geology, 46, 1035–1038, https://doi.org/10.1130/G45415.1, 2018.

Price, R. C., Gamble, J. A., Smith, I. E., Maas, R., Waight, T., Stewart, R. B., and Woodhead, J.: The anatomy of an andesite volcano: A time-stratigraphic study of andesite petrogenesis and crustal evolution at Ruapehu Volcano, New Zealand, J. Petrol., 53, 2139–2189, https://doi.org/10.1093/petrology/egs050, 2012.

Puchol, N., Blard, P.-H., Pik, R., Tibari, B., and Lavé, J.: Variability of magmatic and cosmogenic $^3$He in Ethiopian river sands of detrital pyroxenes: Impact on denudation rate determinations, Chem. Geol., 448, 13–25, https://doi.org/10.1016/j.chemgeo.2016.10.033, 2017.

Pure, L. R., Leonard, G. S., Townsend, D. B., Wilson, C. J., Calvert, A. T., Cole, R. P., Conway, C. E., Gamble, J. A., and Smith, T. B.: A high resolution $^{40}$Ar/$^{39}$Ar lava chronology and edifice construction history for Tongariro volcano, New Zealand, J. Volcanol. Geoth. Res., 403, 106993, https://doi.org/10.1016/j.jvolgeores.2020.106993, 2020.

Ramos, F. C., Heizler, M. T., Buettner, J. E., Gill, J. B., Wei, H. Q., Dimond, C. A., and Scott, S. R.: U-series and $^{40}$Ar/$^{39}$Ar ages of

Holocene volcanic rocks at Changbaishan volcano, China, Geology, 44, 511–514, https://doi.org/10.1130/G37837.1, 2016.

Rhodes, E.: The Draining of an Andesitic Valley-Confined Lava Flow, Mt Ruapehu, Honours thesis, University of Canterbury, 2012.

Schaefer, J. M., Winckler, G., Blard, P.-H., Balco, G., Shuster, D. L., Friedrich, R., Jull, A. J. T., Wieler, R., and Schluechter, C.: Performance of CRONUS-P – A pyroxene reference material for helium isotope analysis, Quat. Geochronol., 31, 237–239, https://doi.org/10.1016/j.quageo.2014.07.006, 2016.

Schimmelpfennig, I., Williams, A., Pik, R., Burnard, P., Niedermann, S., Finkel, R., Schneider, B., and Benedetti, L.: Intercomparison of cosmogenic in-situ $^3$He, $^{21}$Ne and $^{36}$Cl at low latitude along an altitude transect on the SE slope of Kilimanjaro volcano (3° S, Tanzania), Quat. Geochronol., 6, 425–436, https://doi.org/10.1016/j.quageo.2011.05.002, 2011.

Sherrod, D. R., Hagstrum, J. T., Mcgeehin, J. P., Champion, D. E., and Trusdell, F. A.: Distribution, $^{14}$C chronology, and paleomagnetism of latest Pleistocene and Holocene lava flows at Haleakala Island of Maui, Hawaii: A revision of lava flow hazard zones, J. Geophys. Res., 111, B05205, https://doi.org/10.1029/2005JB003876, 2006.

Shuster, D. L., Farley, K. A., Sisterson, J. M., and Burnett, D. S.: Quantifying the diffusion kinetics and spatial distributions of radiogenic $^4$He in minerals containing proton-induced $^3$He, Earth Planet. Sc. Lett., 217, 19–32, https://doi.org/10.1016/S0012-821X(03)00594-6, 2004.

Smith, J. A., Finkel, R. C., Farber, D. L., Rodbell, D. T., and Seltzer, G. O.: Moraine preservation and boulder erosion in the tropical Andes: Interpreting old surface exposure ages in glaciated valleys, J. Quaternary Sci., 20, 735–758, https://doi.org/10.1002/jqs.981, 2005.

Stipp, J.: The Geochronology and Petrogenesis of the Cenozoic Volcanics of the North Island, New Zealand, PhD thesis, Australian National University, https://www.proquest.com/docview/2617247792?pq-origsite=gscholar&fromopenview=true&sourcetype=Dissertations&Thesis (last access: 9 July 2023), 1968.

Stone, J. O.: Air pressure and cosmogenic isotope production, J. Geophys. Res., 105, 753–759, https://doi.org/10.1029/2000JB900181, 2000.

Tanaka, H., Kawamura, K., Nagao, K., and Houghton, B. F.: K-Ar ages and paleosecular variation of direction and intensity from quaternary lava sequences in the Ruapehu Volcano, New Zealand, Earth Planets Space, 49, 587–599, https://doi.org/10.5636/jgg.49.587, 1997.

Topping, W. W.: Some aspects of quaternary history of Tongariro Volcanic Centre, PhD thesis, Victoria University of Wellington, https://openaccess.wgtn.ac.nz/articles/thesis/Some_Aspects_of_Quaternary_History_of_Tongariro_Volcanic_Centre/19252133/1 (last access: 9 July 2023) 1974.

Topping, W. W. and Kohn, B. P.: Rhyolitic tephra marker beds in the Tongariro area, North Island, New Zealand, New Zealand. J. Geol. Geop., 16, 375–395, https://doi.org/10.1080/00288306.1973.10431367, 1973.

Townsend, D., Leonard, G., Conway, C., Eaves, S., and Wilson, C.: Geology of the Tongariro National Park Area, GNS Science, pp. 1 sheet + 109 pp., ISBN 978-1-98-853030-7, 2017.
Tremblay, M. M., Shuster, D. L., and Balco, G.: Diffusion kinetics of ³He and ²¹Ne in quartz and implications for cosmogenic noble gas paleothermometry, Geochim. Cosmochim. Ac., 142, 186–204, https://doi.org/10.1016/j.gca.2014.08.010, 2014.

Trusdell, F. A.: Lava flow hazards and risk assessment on Mauna Loa Volcano, Hawaii, Geoph. Monog. Series, 92, 327–336, https://doi.org/10.1029/GM092p0327, 1995.

Tsang, S. and Lindsay, J.: Lava flow crises in inhabited areas part I: Lessons learned and research gaps related to effective, basaltic eruptions, Journal of Applied Volcanology, 9, 9, https://doi.org/10.1186/s13617-020-00096-y, 2020.

Turner, G. M. and Corkill, R. M.: NZPSV11k.2023 and NZPSV1k.2023: Holocene palaeomagnetic secular variation master records for New Zealand, Phys. Earth Planet. In., 344, 107093, https://doi.org/10.1016/j.pepi.2023.107093, 2023.

Turner, G. M., Howarth, J. D., de Gelder, G. I., and Fitzsimons, S. J.: A new high-resolution record of Holocene geomagnetic secular variation from New Zealand, Earth Planet. Sc. Lett., 430, 296–307, https://doi.org/10.1016/j.epsl.2015.08.021, 2015.

Uppala, S. M., Källberg, P. W., Simmons, A. J., Andrae, U., da Costa Bechtold, V., Fiorino, M., Gibson, J. K., Haseler, J., Hernandez, A., Kelly, G. A., Li, X., Onogi, K., Saarinen, S., Sokka, N., Allan, R. P., Andersson, E., Arpe, K., Balmaseda, M. A., Beljaars, A. C., van de Berg, L., Bidlot, J., Bormann, N., Caires, S., Chevallier, F., Dethof, A., Dragosavac, M., Fisher, M., Fuentes, M., Hagemann, S., Hólm, E., Hoskins, B. J., Isaksen, L., Janssen, P. A., Jenne, R., McNally, A. P., Mahfouf, J. F., Morcrette, J. J., Rayner, N. A., Saunders, R. W., Simon, P., Sterl, A., Trenberth, K. E., Untch, A., Vasiljevic, D., Viterbo, P., and Woollen, J.: The ERA-40 re-analysis, Q. J. Roy. Meteor. Soc., 131, 2961–3012, https://doi.org/10.1256/qj.04.176, 2005.

Vermeesch, P.: IsoplotR: A free and open toolbox for geochronology, Geosci. Front., 9, 1479–1493, https://doi.org/10.1016/j.gsf.2018.04.001, 2018.

Villemant, B.: Trace element evolution in the Phlegrean Fields (Central Italy): fractional crystallization and selective enrichment, Contrib. Mineral. Petr., 98, 169–183, https://doi.org/10.1007/BF00402110, 1988.

Wijbrans, J., Schneider, B., Kuiper, K., Calvari, S., Branca, S., De Beni, E., Norini, G., Corsaro, R. A., and Miraglia, L.: ⁴⁰Ar/³⁹Ar geochronology of Holocene basalts; examples from Stromboli, Italy, Quat. Geochronol., 6, 223–232, https://doi.org/10.1016/j.quageo.2010.10.003, 2011.

Wilson, C. J., Gravley, D. M., Leonard, G. S., and Rowland, J. V.: Volcanism in the central Taupo Volcanic Zone, New Zealand: tempo, styles and controls, in: Studies in Volcanology: The Legacy of George Walker, edited by: Thordarson, T., Special Publications of IAVCEI 2, 225–247, https://doi.org/10.1144/IAVCEI002.12, 2009.

Wilson, G., Wilson, T. M., Deligne, N. I., and Cole, J. W.: Volcanic hazard impacts to critical infrastructure: A review, J. Volcanol. Geoth. Res., 286, 148–182, https://doi.org/10.1016/j.jvolgeores.2014.08.030, 2014.

Wright, H. M., Vazquez, J. A., Champion, D. E., Calvert, A. T., Mangan, M. T., Stelten, M., Cooper, K. M., Herzig, C., and Jr, A. S.: Episodic Holocene eruption of the Salton Buttes rhyolites, California, from paleomagnetic, U-Th, and Ar/Ar dating Heather, Geochem., Geophy. Geosy., 16, 1198–1210, https://doi.org/10.1002/2015GC005714, 2015.

Zimmermann, L., Avice, G., Blard, P.-H., Marty, B., Füri, E., and Burnard, P. G.: A new all-metal induction furnace for noble gas extraction, Chem. Geol., 480, 86–92, https://doi.org/10.1016/j.chemgeo.2017.09.018, 2018.