TephraNZ: a major- and trace-element reference dataset for glass-shard analyses from prominent Quaternary rhyolitic tephras in New Zealand and implications for correlation

Abstract. Although analyses of tephra-derived glass shards have been undertaken in New
Zealand for nearly four decades (pioneered by Paul Froggatt), our study is
the first to systematically develop a formal, comprehensive, open-access
reference dataset of glass-shard compositions for New Zealand tephras. These
data will provide an important reference tool for future studies to identify
and correlate tephra deposits and for associated petrological and
magma-related studies within New Zealand and beyond. Here we present the
foundation dataset for TephraNZ, an open-access reference dataset for
selected tephra deposits in New Zealand. Prominent, rhyolitic, tephra deposits from the Quaternary were identified,
with sample collection targeting original type sites or reference locations
where the tephra's identification is unequivocally known based on
independent dating and/or mineralogical techniques. Glass shards were
extracted from the tephra deposits, and major- and trace-element geochemical
compositions were determined. We discuss in detail the data reduction
process used to obtain the results and propose that future studies follow a
similar protocol in order to gain comparable data. The dataset contains
analyses of glass shards from 23 proximal and 27 distal
tephra samples characterising 45 eruptive episodes ranging from Kaharoa (636 ± 12 cal yr BP) to the Hikuroa Pumice member (2.0 ± 0.6 Ma)
from six or more caldera sources, most from the central Taupō Volcanic
Zone. We report 1385 major-element analyses obtained by electron microprobe
(EMPA), and 590 trace-element analyses obtained by laser ablation
(LA)-ICP-MS, on individual glass shards. Using principal component analysis (PCA), Euclidean similarity coefficients, and geochemical investigation,
we show that chemical compositions of glass shards from individual eruptions
are commonly distinguished by major elements, especially CaO, TiO2,
K2O, and FeOtt (Na2O+K2O and SiO2/K2O), but not
always. For those tephras with similar glass major-element signatures, some
can be distinguished using trace elements (e.g. HFSEs: Zr, Hf, Nb; LILE: Ba,
Rb; REE: Eu, Tm, Dy, Y, Tb, Gd, Er, Ho, Yb, Sm) and trace-element ratios
(e.g. LILE/HFSE: Ba/Th, Ba/Zr, Rb/Zr; HFSE/HREE: Zr/Y, Zr/Yb, Hf/Y;
LREE/HREE: La/Yb, Ce/Yb). Geochemistry alone cannot be used to distinguish between glass shards from
the following tephra groups: Taupō (Unit Y in the post-Ōruanui
eruption sequence of Taupō volcano) and Waimihia (Unit S); Poronui (Unit
C) and Karapiti (Unit B); Rotorua and Rerewhakaaitu; and
Kawakawa/Ōruanui, and Okaia. Other characteristics, including
stratigraphic relationships and age, can be used to separate and distinguish
all of these otherwise-similar tephra deposits except Poronui and Karapiti.
Bimodality caused by K2O variability is newly identified in Poihipi and
Tahuna tephras. Using glass-shard compositions, tephra sourced from
Taupō Volcanic Centre (TVC) and Mangakino Volcanic Centre (MgVC) can be
separated using bivariate plots of SiO2/K2O vs.
Na2O+K2O. Glass shards from tephras derived from Kapenga
Volcanic Centre, Rotorua Volcanic Centre, and Whakamaru Volcanic Centre have
similar major- and trace-element chemical compositions to those from the
MgVC, but they can overlap with glass analyses from tephras from Taupō and
Okataina volcanic centres. Specific trace elements and trace-element ratios
have lower variability than the heterogeneous major-element and bimodal
signatures, making them easier to fingerprint geochemically.


Trace elements are more strongly partitioned by fractional crystallisation processes that occur during the formation of melt and therefore have the potential to be unique for discrete eruption episodes (e.g. Pearce et al., 2004). Specifically, a number of key trace elements have been identified as important for the correlation of rhyolitic tephras, including the high field strength elements (HFSEs) Zr and Nb; the large ion lithophile elements (LILEs) Rb, Sr, and Ba; the heavy rare-earth elements (HREEs) Gd, Yb, Sc, and Y; and the light rare-earth elements (LREEs) La and Nd. Trace-element ratios are also identified as important, including (1) HFSE/HREE -for example Zr/Y, Nb/Y, Hf/Y; (2) LILE/HFSE -for example Ba/Th; (3) LREE/HFSE -for example Ce/Th, La/Nb; (4) LREE/HREE -for example La/Yb, Ce/Yb; and (5) HFSE/HFSE -for example Zr/Nb, Zr/Th. Some studies have shown that trace elements and trace-element ratios can distinguish between tephra beds that have indistinguishable glass-shard major-element signatures and thus are a robust way of providing accurate correlations (e.g. Westgate et al., 1994;Pearce et al., 1996Pearce et al., , 2002Pearce et al., , 2004Allan et al., 2008;Hopkins et al., 2017).
Tephra correlation is also increasingly being quantified through statistical approaches on geochemical data (Lowe et al., 2017), but many of these approaches (e.g. supervised learning) often require a robust, comprehensive set of "known" reference data against which to test the analyses of "unknown" samples. Statistics can also scale data to make them comparable, but they cannot account or correct for inter-laboratory or historical variance in analyses. Therefore, incomplete datasets, or datasets constructed from a range of data sources, will limit the ability to provide holistic statistical correlations with accurate outputs. Consequently, the formation of reference datasets that are run in one analytical session, in one lab, with a consistent methodology is highly desirable for minimising sources of error. The production of tephra databases is thus being recognised as an exceptionally useful tool internationally (e.g. Lowe et al., 2017), made more obtainable with open-access journals and online, effectively limitless storage, leading to easier publication and maintenance of large data repositories. Ideally, a global tephra database would exist, but at present this is beyond the scope and remit of any individual researcher, research group, or institute(s). Therefore separate, regional databases for volcanically active (and other) regions are becoming increasingly popular, such as TephraKam -Kamchatka (Portnyagin et al., 2020); Tephrabase -Europe (Newton, 1996); AntT tephra database -Antarctic ice cores (Kurbatov et al., 2014); Alaska Tephra Database (Wallace, 2018); Klondike goldfields, Yukon (Preece et al., 2011); VOLCORE -DSDP, ODP, and IODP marine tephra deposits (Mahony et al., 2020). In addition, in an effort to produce comparable global datasets Abbott et al. (2021) have recently presented guidance for best practices, providing recommendations and templates for tephra collection, sample preparation, and physical and geochemical analysis.

Geologic setting
The volcanically active nature of New Zealand (Mortimer and Scott, 2020) and the longevity and consistency of largescale rhyolitic eruptions (Howorth, 1975;Froggatt and Lowe, 1990;Houghton et al., 1995;Wilson et al., 1995aWilson et al., , 2009Jurado-Chichay and Walker, 2000;Carter et al., 2003;Briggs et al., 2005;Wilson and Rowland, 2016;Barker et al., 2021) mean the landscape currently has a very long, detailed, and complex rhyolitic tephrostratigraphic framework that is used for a wide range of applications . However, at present New Zealand tephra studies are lacking a comprehensive reference dataset resource that has been developed in a systematic way.
Because of their pervasive nature, high repose period, and high preservation potential, tephra deposits are a common stratigraphic and chronological aid in many studies in New Zealand (Shane, 2000;Lowe, 2011;Hopkins et al., 2021). For example, the eruption of Kaharoa (636 ± 12 cal yr BP, Hogg et al., 2003) from Mt Tarawera in the Okataina Volcanic Centre (OVC) has been used to date the arrival of Polynesians in northern New Zealand and map their expansion and impact across the country (Newnham et al., 1998;Lowe and Newnham, 2004). The Rerewhakaaitu eruption (17 496 ± 462 cal yr BP; Lowe et al., 2013), sourced from OVC, is used as a marker horizon for the transition between the last glacial and present interglacial , and several other widespread late Quaternary tephra deposits form boundaries or key stratigraphic markers in the New Zealand Climate Event Stratigraphy developed by the NZ-INTIMATE community (e.g. Kawakawa/Oruanui tephra; Barrell et al., 2013;Lowe et al., 2013). Compositions of glass and mineral components from rhyolitic tephra deposits have also been used to reconstruct changes in magmatic systems and give insight into the complexity of caldera-related eruption episodes (e.g. Smith et al., 2002Smith et al., , 2005Cooper et al., 2012;Barker et al., 2016Barker et al., , 2021Wilson and Rowland, 2016). Many of the commonly found rhyolitic tephra horizons in New Zealand are well studied, dated, and geochemically and mineralogically characterised. However, often these studies have been eruption-, source-, or depocentre-specific, and thus only provide a small, effectively piecemeal catalogue of tephra geochemistry that is not necessarily comparable to those of other studies. In addition, compositional data are not usually published in their entirety, or not at all, meaning future studies can neither access nor use the data for correlation techniques. Furthermore, Lowe et al. (1999) identified that differing procedural methods employed at different institutes around New Zealand before and after 1995 produced variable elemental concentrations for the same tephra (post-1995 SiO 2 values were lower by 0.5 wt %-1.0 wt %, and all other elements had slightly higher values). Therefore, it is likely that some of the older tephra compositions that have been relied upon in the past for correlative purposes are no longer appropriate.
It is therefore timely for a comprehensive, systematic, and accessible New Zealand tephra database to be established and curated. In this study we present TephraNZ as a foundation reference dataset of internally consistent, openaccess data for major-and trace-element compositions of glass shards from a selection of the most pervasive Quaternary tephra deposits in New Zealand (Table 1). This is by exact sampling location. We note here that we have not attempted to sample multiple tephra beds from a single eruptive episode in proximal sequences, nor deposits of the same tephra at different azimuths, as has been undertaken in some more localised or petrologically focussed studies (e.g. Shane et al., 2005Shane et al., , 2008. We recognise this limitation but instead have concentrated on analysing a wide range of pervasive rhyolitic tephras, both proximal and distal, in a systematic and well-documented way so that future tephrostratigraphic studies will have a foundation of new, high-quality glassshard compositional data for facilitating robust correlations and applications. Where we have both, we compare proximal and distal analyses of the same tephra and comment on similarities or differences allowing for an increased understanding of the variability in the geochemistry seen in the pyroclastic products of some eruptions. In addition, we have used statistical methods to demonstrate the integrity of our new datasets (and show how such methods can enable unknown tephras to be classified).

Sample preparation
Bulk tephra samples were disaggregated in water for 1-5 min in an ultra-sonic water bath. Clays and ultra-fine sediments (< 5 µm) were rinsed off, and samples were then wet-sieved using disposable sieve cloths to 125-250 µm or, where necessary, 60-125 µm. Samples were then dried for 12-24 h at 50 • C before mounting in epoxy resin. Seven samples were mounted into individual drill holes (4 mm diameter) in 25 mm epoxy round blocks (a 4 : 1 ratio of EpoTek 301 resin [A]: hardener [B]). Individual drill holes were then backfilled using the same epoxy mix (see Lowe, 2011, p. 124, for a schematic illustration). Sample blocks were polished using the following sequence: ∼ 3 min in a figure-eight pattern on 800-grit sandpaper with water lubricant to remove the epoxy and break through to the glass shards, ∼ 1 min on 1200-grit sandpaper with water lubricant to remove any large scratches, and ∼ 1 min on 2500-grit sandpaper with water lubricant to begin to reveal the outline of the shards. Blocks were then moved on to the diamond laps with their appropriate lubricant, all at 280 revolutions per minute rotating the block 90 • every 30 s followed by 2 min of ultrasonic bathing at < 24 • C between each lap stage to remove any loose material on the surface of the blocks: ∼ 3 min on 6 µm, ∼ 1 min at 3 µm, and ∼ 1 min at 1 µm. Blocks were then carbon coated before loading in the electron microprobe system for analysis. Table 1.

EPMA method and data reduction
Major-element analysis of glass shards was undertaken at Victoria University of Wellington (VUW) by wavelength dispersive X-ray spectroscopy (WDS) on a JEOL JXA8230 Superprobe electron probe microanalyser (EPMA). Broadly the method follows that espoused by Kuehn et al. (2011). Backscatter electron images of each sample were taken and used as block maps to allow the location of EPMA analyses to be replicated for trace-element analysis. A defocused circle beam 10 µm in diameter was used at 8 nA and 15 kV to analyse all major elements as oxides (SiO 2 , TiO 2 , Al 2 O 3 , FeO tt , MnO, MgO, CaO, Na 2 O, K 2 O) and Cl. Run duration for each analysis was ∼ 3 min, including online correction. During standardisation, Na 2 O was run twice, the second time skipping the peak search to reduce the volatilisation of the element, with the second standardisation value then used. Supplement Table S1.1 (a and b) shows the EPMA set-up and run times (after Abbott et al., 2021). During the analysis, VG-568 was run as a calibration standard, and VG-A99 and ATHO-G were run as secondary standards (all standard data can be found in Supplement Table S3), with two of each standard (calibration and secondary) analysed between 10 sample analyses to monitor machine drift (no machine drift was identified).
Initial concentrations were determined using the ZAF correction method, with secondary offline data reduction undertaken to all samples and standards to correct for variability in VG-568. Internal correction values were calculated using the GeoREM reference values of VG-568 from Streck and Wacaster (2006; Eq. 1) and applied to all the data (Eq. 2). Following this, samples were corrected for deviations from 100 wt % total; this assumes any variation is due mostly to magmatic water, with a very small amount of minor and trace elements (Froggatt, 1983;Lowe, 2011) that are not analysed by the EPMA (Eq. 3). The difference is reported as "H 2 O D " in all data tables to allow back calculation to original data values including totals. Results with H 2 O D ≥ 8 wt % were removed and are listed at the bottom of the table as "outliers" (Supplement Table 2). Accuracy and analytical precision of the standards were calculated, where accuracy is the offset from the reference value for the secondary standards (Eq. 4), and precision is the standard deviation of all measured secondary standards throughout a run, reported at 2 standard deviations (SD) to represent a 95 % variability.
where X p m = measured concentration of element X of the calibration standard, and X p r = reference concentration for element X of the calibration standard (reference values taken from GeoRem preferred values http://georem.mpch-mainz. gwdg.de/, last access: June 2021).
where X i m = measured concentration for element X of any sample or standard. secondary hydration corrected data = ((corrected data (Eq. 2)/total for that sample) × 100) (3) offset from standard (accuracy) = (absolute value(X s r − averageX s m )), where X s r = reference concentration for element X of the secondary standard (GeoRem preferred value; MPI-DING; Jochum et al., 2006), and averageX s m = average concentration measured for element X of all analyses of the secondary standard).

LA-ICP-MS method and data reduction
In situ trace-element analysis was undertaken at VUW using laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) where a RESOlution S155-SE 193 nm ArF excimer laser system was coupled with an Agilent 7900 quadrupole ICP-MS. Data for 43 trace elements were acquired using a static spot method, with a 25 µm spot size, ablation time of 30 s, and repetition rate of 5 Hz power (method: 10 s background/washout count, cleaning spot of 25 µm for three laser pulses to clean the glass shard surface, 20 s background count, 30 s acquisition, 10 s washout; see Supplement Table S1 for full LA-ICP-MS set-up details; after Abbott et al., 2021). Synthetic glass standards NIST-612 and NIST-610 were used to tune the ICP-MS and obtain the P/A factors at a range of spot sizes and laser powers. During the analysis, a full range of standards was analysed to determine which produced the most accurate and precise results as a calibration standard, including NIST-612, NIST-610, BHVO2-G, and ATHO-G. StHS6/80-G was analysed as a secondary standard throughout (results of which are discussed below in Sect. 2.5 and shown in Fig. 2 and Supplement Table S5), and all standards (calibration and secondary) were analysed twice every 10 samples. All data were reduced offline using Iolite v.3™ software (Paton et al., 2011), using 43 Ca as the internal standard value (index channel) and the Trace_ Elements_IS data reduction scheme (DRS). The data were reduced against ATHO-G as the calibration standard. No post-processing data reduction was necessary for the trace-element data, but outliers were removed; precision and accuracy were calculated on STHS6/80-G as described above (Eq. 4).

Standardisation method
Multiple calibration standards with different trace-element concentrations were analysed to determine which would be most suitable for trace-element data reduction. Potential calibration standards included NIST-612, NIST-610, BHVO2-G, and ATHO-G. These were each run twice every 10 samples, along with secondary standard STHS6/80-G. Figure 2 shows the STHS6/80-G results of a range of selected, commonly used trace elements, including Zn (transition metal), Rb (LILE), Zr (HFSE), La (LREE), and Yb (HREE), normalised using each of the calibration standards. Overall, the results show that for the lighter masses (e.g. Zn) there is a large variability in the measured STHS6/80-G values across the different standards, but all except BHVO2-G sit within error (2 SD) of the reference value (Fig. 2). For the heavier masses (e.g. La, Yb, Fig. 2), the variation from the reference value observed within the analysed values decreases, except for NIST-610, which remains highly variable in the middle masses (Rb, Zr, Fig. 2), with variability reducing in the heavier masses. The data show that the use of ATHO-G as the calibration standard (for data reduction of rhyolites) produces the most accurate and precise data for the secondary standard, for all except the elements with the heaviest masses and smallest concentrations (e.g. Yb).

Principal component analysis
To visualise elements that distinguish the different tephra compositions we have used principal component analysis (PCA). PCA was run in the coding platform R (R core team, 2019) v3.6.2 and RStudio v.1.2.5033 using packages "ggbiplot" (Vu, 2011), "Hotelling" (Curran, 2018), "ggplot2" (Wickham, 2016), "factoextra" (Kassambara and Mundt, 2020) and "vegan" (Oksanen et al., 2019). Data for Tuhua tephra were removed as these would unnecessarily skew the results due to their distinct geochemistry. Non-normalised, average, elemental values were used from each tephra sample, for example Si (in ppm), no oxide values, no ratios (e.g. SiO 2 /K 2 O), or sums (e.g. Na 2 O+K 2 O). All element values were centred using a centred log-ratio transform to deal with closure effect (clr: column mean subtracted from each value) and scaled (value divided by the standard deviation of the column) to compare elements with concentrations that differ by orders of magnitude. PCA was run using the "prcomp" (Venables and Ripley, 2002) function, and PCA contributions were calculated using "fviz_contrib" (Kassambara and Mundt, 2020) function. A template of the coding script used can be found in Supplement Material 1.

Euclidean similarity coefficients
To identify the tephra samples that were most similar, and could therefore pose problems in attempting to obtain unique fingerprinting, we ran Euclidean similarity coefficient (ESC) analysis. ESC was run in R and RStudio using the package "stats" (R core team, 2019). Following the guidelines of Hunt et al. (1995) for ESC analysis, we used non-normalised, mean concentrations of the elements highlighted by the PCA to be the most indicative of variance in the dataset. These values were input as comparison values, and the function "as.matrix.dist" was used to run the "Euclidean" distance measure. This method calculates the similarity of samples based on an infinite number of comparison input values. A template of the coding script used can be found in Supplement Material 2. The output table was manipulated postproduction to provide the colour formatting.

Results
The averages and their standard deviations for all samples are reported in Table 2; the full reference dataset can be found in Supplement Table S2. All reported values in the text and figures (unless stated otherwise) are recalculated (normalised) to 100 % on a volatile-free basis (following Lowe et al., 2017) with the difference between the raw total and 100 % being reported as "H 2 O D " (Table 2). For best correlation results, we recommend that the full dataset is used in order to see the trends in the geochemical data rather than just the means and standard deviations.

Data quality
Standard values for VG-568 and VG-A99 are taken from the GeoREM. The reference values used as a standard (by other publications) are from Jarosewich et al. (1980). However, for the purpose of this research we have chosen alternate values published by Streck and Wacaster (2006). A comparison of the reference values from both publications is shown in Supplement  Jarosewich et al. (1980). We do note, however, that Cl is challenging to analyse accurately on EPMA for glass due to its low concentration and especially as there are few standards that have similar compositions (e.g. Jochum et al., 2006). Our samples have between 0.3 wt % and 0.06 wt % Cl; therefore, VG-568 (with Cl = 0.1 wt %), ATHOG (with 0.04 wt %), and VG-A99 (with 0.02 wt %) attempt to provide a good range for standard comparison.

Figure 2.
Compilation of trace-element standard data produced during the first run of glass-shard analyses. These data show selected element concentrations of secondary standard STHs6/80 normalised using difference calibration standards (NIST-612 -orange; NIST-610 -green, ATHO-G -blue, and BHVO2-G -red). The grey shaded area shows the preferred GeoREM reference value (http://georem.mpch-mainz.gwdg.de, last access: June 2021) error margin reported for each element for STHS6/80. Note that for standard NIST610, two data points were removed as outliers; full data can be found in Table S5. Table 2. Average and standard deviation values (shown in italics) for glass-shard analysis of major and trace elements. Full data can be found in the Supplement Table S2 Table 2.

Continued.
Tephra name   .4) show the variability in the concentrations analysed by the EPMA of the standard data throughout the running of these samples. For the secondary standards ATHO-G ( Fig. S3.3) and VG-A99 (Fig. S3.4), there is some clear variability within the batches of samples run. For example for SiO 2 for ATHO at point 60, there is a clear jump in the values reported, for Na 2 O for ATHO at points 92-101 there are some very low concentrations, and for MgO for VG-A99 there are clear variations in different run sets. The variation observed in all these data is likely due to a number of factors including (1) a change to a different standard shard during the run; (2) re-calibration of the EPMA after a period of down time (note the dates of analyses) and/or (3) day-today variations in machine performance; (4) for the case of Na 2 O, possible volatilisation of Na 2 O due to repeated analysis of the same standard shard; (5) use of an inappropriate primary calibration standard (for example a rhyolite standard (VG568) to calibrate a basaltic glass (VG-A99) which is used in this case as a secondary standard).
Recently, a number of studies have reported difficulty in accurate analysis of Na 2 O concentrations in reference standard ATHO-G: (1) a large range of values are reported from different analytical techniques (3.53 wt %-4.31 wt %, Jochum et al., 2006), and (2) the reported reference value for ATHO-G from Jochum et al. (2006) is too low (reported as 3.75 wt %; e.g. Lowe et al., 2017;Portnyagin et al., 2020). Our use of Steck and Wacaster (2006) reference data for VG-568 as an internal calibration (3.52 wt %) rather than the Jarosewich et al. (1980) value (3.75 wt %) brings our secondary standard data in alignment with the original Jochum et al. (2006) values for ATHO-G (see Supplement  Table S3.3). However, it is possible that because of this, our sample values reported for Na 2 O are too low. Further studies into this community-wide issue will hopefully allow this discrepancy to be resolved.
During LA-ICP-MS tuning oxide production was monitored using the ThO / Th ratio, and this value was tuned to between 1.3 % and 1.8 %, which is considered high by current standards (e.g. Portnyagin et al., 2020, reported values of 0.5 %-0.7 %), but is comparable with values in older studies (e.g. Jochum et al., 2006, report values "< 1 %-2 %"; Pearce et al., 2011, reported values "typically ∼ 1.5 %";Allan et al., 2008, reported values "typically < 1 %, always < 2 %"). This high oxide production value could have had impacts on some elements. For example, it is likely that there was a high addition of SiO into our analyses; SiO can interfere with 45 Sc, or alternatively, BaO can interfere with 153 Eu. However, be-cause the concentration of SiO 2 in our samples is similar to that of our secondary standard (ATHO-G), the data should still be viable. In addition, to monitor the impact of oxides on our elements we analyse and report multiple isotopes of Sr (86 and 88), Zr (90 and 91), Mo (95 and 98), Ba (137 and 138), and Eu (151 and 153). The concentrations of these elements do not show significant variability (e.g. Figs. S6.2.3, and 6.2.4; R 2 Sr = 0.96, R 2 Zr = 0.99). In addition, when plotted together, Ba vs. 153 Eu shows no relationship (Fig. S6.2.5), proving little-to-no oxide interference has impacted the values obtained for these elements. Ti
Of the 45 tephra samples, 22 have a "homogeneous signature", homogeneity being defined here (as an approximation) when the standard deviation of the sample is equal to or less than analytical error (2 SD of secondary standard: for example, for FeO t = ± 0.23 wt %, CaO = ± 0.10 wt %). The majority (∼ 64 %) of the samples that have a homogeneous signature are from OVC (e.g. Whakatāne, Mamaku, Rotoma) or from calderas older than OVC (∼ 32 %), such as (1) Upper Griffins Road tephra, a correlative of the Whakamaru eruptives, Whakamaru Volcanic Centre (WVC), and (2) Mangapipi tephra, a correlative to deposits of Mangakino Volcanic Centre (MgVC; Fig. 5a). Ten samples show a heterogeneous signature (where standard deviations for both FeO t and CaO are greater than analytical errors), with most from a proximal source (∼ 30 %), or from tephras deposited in the Whanganui Basin area (40 %), and with the remainder being from the Mangaone Subgroup eruptives from the OVC: Hauparu, Maketū, and Ngāmotu (Fig. 5b).
Glass shards from four tephra samples show a bimodal signature in some major and trace elements, where specific elements split the populations into two distinct groups. Tephras showing this phenomenon include Rotorua (OVC), Rerewhakaaitu (OVC), Poihipi (TVC), and Tahuna (TVC). The bimodal signatures of Rerewhakaaitu and Rotorua are well documented , whereas those of Poihipi and Tahuna are newly identified here (Fig. 6). All four of these tephra horizons have their glass-shard bimodal signatures produced predominantly by K 2 O concentrations, into high (≥ 3.8 wt %) and low (≤ 3.6 wt %) populations (Fig. 6) linked to the crystallisation of biotite minerals. This relationship has been discussed in previous research (for Rerewhakaaitu and Rotorua) and modelled as two different biotite populations formed through fractional crystallisation in a zoned magma chamber (e.g. Shane et al., 2003;, resulting in the formation of heterogeneity in the magma and hence the formation of different glass compositions. For five of the tephras, we undertook analyses on glass from both proximal and distal samples. These tephras included Whakatāne, Rotoma, Waiohau, Rotorua, and Rerewhakaaitu, which are all derived from OVC (Table 1). For Rotoma, Rerewhakaaitu, and Waiohau, the signatures of the proximal and distal deposits are indistinguishable, whereas for Whakatāne and Rotorua the proximal signature is highly variable, and the distal signature is homogeneous but overlapping with part of the extent of the proximal signature (e.g. Fig. 7). Similar findings are reported and discussed in more detail for Whakatāne tephra in Kobayshi et al. (2005) and Holt et al. (2011) and for Rotorua tephra in Shane et al. (2003a) and Kilgour and Smith (2008). Figure 8 shows a primitive mantle-normalised spider plot of all the trace-element data for the glass shards analysed (after McDonough and Sun, 1995). The majority of the data plot along a common pattern of variable concentrations of HFSE, LILEs, and LREEs, but they show more consistent concentrations of HREEs (Gd to Lu). Of note are peaks in Nd, a negative Sr anomaly relative to LREE, and a positive Zr-Hf anomaly relative to Sm. Sr and Ba show the largest variability in concentrations that is likely caused by a variability in feldspar crystallisation (Pearce et al., 2004). Several different patterns are observable within this full data suite pertaining to individual samples. The obviously different signature is that for glass from Tuhua tephra, which shows a low concentration of Ba (< 10 ppm) and Sr (< 1 ppm) in comparison with values for the rest of the samples, and with high concentrations of all other elements, especially the REEs (Fig. 8). Analyses of glass shards from the Maketū tephra can also be identified by their high concentrations of all elements in comparison to the TVZ trends but mid-range Nb values (between those of Tuhua and the general trend) (Fig. 8). We also note Er and Lu peaks, which pertain to glasses from the Te Rere tephra, that sit at the higher concentration levels of the general trend (these could potentially be analytical artefacts; Fig. 8; Table 2) and samples from Ngāmotu, Rotoehu/Rotoiti, and Earthquake Flat deposits that sit at the lower overall trace-element concentration levels of the general trend (Fig. 8). For the tephras where both proximal and distal samples of glass have been analysed for trace elements, the HFSEs (including Zr, Hf, Th, and Ti) and LILEs (including Rb, Sr, and Cs) may exhibit heterogeneity between the proximal and distal samples, whereas the HREE and the LREE tend to have a lower variability (Fig. 7).

Major and trace elements in general
In many cases, the major-element concentrations in glass are sufficient to allow different tephras to be distinguished (including through the common use of biplots), a result consistent with the findings from much previous work both in New Zealand and elsewhere (e.g. Lowe et al., 2017). However, previous studies have also shown that for some New Zealand tephras more elements from the glass analyses are often required to distinguish between tephras from different eruptions. For this reason we used principal component analysis (PCA) on the dataset to compare multidimensional data rather than an array of traditional biplots. Looking at data in a multidimensional space can allow variations to be more readily distinguished and visualised because all constituent elements are used, not just two.
PCA results for the glass-shard major elements (Fig. 9) show that PC1 and PC2 explain 82.7 % of the variance within the data. Al, K, Si, Na, and Ti make the highest contributions to PC1 (Fig. 9), while Fe, Mn, and Cl have the great-est loadings on PC2 (Fig. 9); therefore, these elements are most appropriate for distinguishing between tephra deposits for the reference dataset as a whole (Fig. 9). These major elements, especially Fe and K (±Ca), have long been recognised as being useful to distinguish many New Zealand late Quaternary tephras from one another (e.g. Lowe, 1988;Shane, 2000;Alloway et al., 2013); however, the inclusion of Ti, Al, Mn, and Cl is somewhat unusual. In a number of cases (discussed below), however, major-element concentrations are shown to overlap for certain tephra horizons, and thus trace elements and trace-element ratios are investigated to provide additional variables to use as discriminants. PCA was also applied to scaled trace elements and major elements together, with the results indicating that PC1 and PC2 could explain 62.8 % of the variability in the full data suite with V, Co, Mg, Cu, Ti, Sr, Sc, Ca, Cs, and Zr being the 10 highest contributors to PC1 and Cu, Mn, Mg, Cs, Sc, Co, Ti, Sr, Th, and Rb highlighted as the 10 highest to PC2 (Fig. 10). Therefore, these elements, and the ratios of these elements, have the highest potential to distinguish individual tephra horizons when using their glass-shard compositions alone.

Source-specific major and trace elements
The central TVZ contains nine recognised calderas, each with different eruption histories but all having produced large-magnitude/volume tephra-producing rhyolitic eruptives. Some of the calderas are attributed to single caldera collapse events (Rotorua, Reporoa, and Ohakuri), whereas others represent composite collapse events that overlap spatially but not temporally (Mangakino and Kapenga). However, the majority reflect multiple collapse events over an extended period of time (Maroa, Okataina, Taupō, and Whakamaru) ( Fig. 1; Wilson et al., 1995aWilson et al., , 2009Barker et al., 2021). Although the calderas are mostly discrete in space, evidence from multiple eruptions has shown their plumbing systems may be linked tectonically (e.g. Wilson et al., 2009;Allan et al., 2012). Hence, the ability to trace a tephra deposit to a caldera source through glass-shard geochemistry alone could be challenging.
The results of the PCA analysis suggest that tephra sourced from the TVC can be distinguished from those of a proposed Mangakino source (MgVC) (Fig. 9). Using SiO 2 /K 2 O vs. Na 2 O+K 2 O, the glass shards of the TVC tephras generally have higher SiO 2 /K 2 O and lower Na 2 O+K 2 O in comparison to those of the equivalent oxides for MgVC-sourced tephra (Fig. 11a). This information is important, but because of the large age differences for the calderas (TVC ∼ 0.32 Ma to present, and MgVC ∼ 1.6 to 1.53 Ma and ∼ 1.2 to 0.95 Ma), the use of this distinction is likely more important for discussions on mantle source dynamics rather than for geochemical correlation of tephra deposits.
Previous studies have suggested that the geochemical characteristics of glass shards from TVC and OVC tephra deposits post-dating the eruption of the Kawakawa/Oruanui (KOT) can be distinguished using fO 2 of Fe-Ti oxides and minerals (Shane, 1998), pumice and lava compositions (Sutton et al., 2000), and glass chemistry (Stokes et al., 1992). Our results also show there is a bimodality in the TVC glassshard data as a whole and that the post-KOT tephra deposits from the TVC and OVC are quite different, whereas the pre-KOT tephras from OVC and TVC are similar ( Fig. 11c  and d). Most glass shards erupted after the KOT event from the TVC have low SiO 2 (≤ 77 wt %), less variable K 2 O (∼ 3 wt %), and higher values for all other major elements in comparison with those of the glass shards erupted from the OVC (Fig. 11c and d). In contrast, tephras erupted from the TVC and OVC prior to, and including the KOT, do show a large amount of overlap in their glass geochemical signatures. For OVC (Fig. 11e and f), there is a high density of samples that have their SiO 2 concentrations at ∼ 78 wt %; however, there is a high variability in SiO 2 overall, with Maketū, Hauparu, and Ngāmotu of the Mangaone Subgroup plotting with SiO 2 concentrations ≤∼ 76 wt % and the remaining Mangaone Subgroup samples (Unit L, Awakeri, and Mangaone) clustering at SiO 2 = 76 wt %-77.5 wt % (FeO t = ∼ 1.2 wt %, K 2 O = ∼ 2.8 wt %, Al 2 O 3 =∼ 13 wt %, CaO = ∼ 1.2 wt %), a finding consistent with those of , who divided the Mangaone Subgroup into "old" and "young" eruptives on the basis of low and high SiO 2 , respectively, unlike the other OVC-sourced samples that plot around SiO 2 =∼ 77.5 wt %-79 wt %, (FeO t = ∼ 0.8 wt %-0.9 wt %, K 2 O = 2.75 wt %-4.5 wt %, Al 2 O 3 =∼ 12 wt %-13 wt %, CaO = ∼ 0.5 wt %-1.0 wt %; Fig. 11c and d). Analyses from the Rotoehu/Rotoiti tephra deposits plot independently from those of other OVC eruptives for this time period. However, they overlap with those of some TVC-tephraderived glass compositions (Poihipi, Tahuna, Okaia, and KOT). The Rotoehu/Rotoiti tephra deposits have a markedly homogeneous geochemical signature and are also much older than TVC eruptions (Table 1). Hence, coupled with the thickness of the deposits, it is likely that a tephra linked to the Rotoehu/Rotoiti eruption would be obvious to distinguish through stratigraphic relationships and age combined with the geochemistry.
The TephraNZ dataset presented here also includes analyses of glass of samples from tephras erupted from the Kapenga Volcanic Centre (KVC; Earthquake Flat eruption), Rotorua Volcanic Centre (RoVC; Ararātā Gully and Kākāriki), and Whakamaru Volcanic Centre (WVC; Fordell, Upper and Lower Griffins Road, Potaka, Rewa, Pakihikura, and Mangapipi). In addition, some older tephra deposits have been recorded in the Whanganui Basin and elsewhere. These are well-known beds, but their caldera sources are not yet defined (Alloway et al., 1993;Pillans et al., 1994Pillans et al., , 2005Shane et al., 1996;Rees et al., 2018Rees et al., , 2019Rees et al., , 2020. Figure 12 shows a plot for the data from KVC, RoVC, and WVC with those regions populated by glass data from samples from the OVC, TVC, and MgVC sources. Overall, the samples plot with a lower SiO 2 /K 2 O ratio (≤∼ 25) similar to that of the MgVC-sourced tephra, which seems to be indicative of samples from older sources in comparison with those from the OVC and TVC. The samples potentially linked to RoVC (Bussell, 1986;Bussell and Pillans, 1997) show different geochemical compositions. For example, Kākāriki-tephraderived glass has slightly higher SiO 2 ≥ 78 wt %) in comparison to that of the Ararātā Gully tephra (SiO 2 ≤ 77 wt %), Figure 5. Examples of major-element bivariate plots for glassshard analyses of tephras (presented on a normalised basis) which show (a) homogeneous signatures, where the standard deviation of the analysis is less than the analytical error (shown as 2σ ), and (b) heterogeneous signatures, where the standard deviation of the analysis is greater than the analytical error. Different colours indicate the differing caldera sources (shown in Fig. 1), and different symbols show the different tephras. P -proximal sample (see Table 1). Total iron is expressed as FeO.
suggesting that they are likely derived from different eruptions but potentially the same source (Mamaku Ignimbrite reportedly has variable geochemical phases; Milner et al., 2003). Glass from the KVC sample (Earthquake Flat tephra) has a very homogeneous signature in the major elements but a more variable signature in the trace elements, both of which overlap with OVC-and TVC-source signatures. There is a very large spread for the data from the unknown samples, precluding the ability to specify their source based simply on major and trace elements alone. Nevertheless, their glass compositional signatures are more similar to those of the older MgVC-sourced tephra, in comparison to those of the younger TVC and OVC deposits, as would be expected based on their known age range (Table 1). Figure 6. Selected major-element biplots of glass analyses (presented on a normalised basis) of samples from Poihipi and Tahuna tephras (both TVC sourced) that exhibit a bimodal signature. This bimodality is identified as being caused by K 2 O concentration (e.g. see Lowe et al., 2008;Shane et al., 2008), and therefore plots with other elements (major or trace) do not show this bimodality. Total iron is expressed as FeO.

Homogeneous, heterogeneous, and bimodal samples
Fingerprinting of glass shards for correlation relies on the ability to distinguish between different deposits, and therefore a homogeneous signature for a single eruptive that is dis- tinct from all other samples is the ideal "fingerprint". However, sometimes there is more complexity in the geochemical data, and heterogeneity can develop in a tephra deposit through a number of mechanisms including the following (Lowe, 2011): tephras. The Rotoehu Ash signature is also similar to that for the Rotoiti Ignimbrite (which are coeval deposits; Nairn and Kohn, 1973), the Earthquake Flat tephra (Kapenga VC; Nairn and Kohn, 1973), and the Ngāmotu tephra (OVC; Jurado-Chichay and Walker, 2000).
For example, the heterogeneous signature identified for glass from the Kaharoa tephra agrees with previous findings for this eruptive. Nairn et al. (2004) and Sahetapy-Engel et al. (2014) reported that glass compositional variability within the Kaharoa deposits shows sequential tapping of a stratified magma body coupled with syn-eruptive changes in dispersal patterns. In general, this is likely one of the reasons why some of the proximal tephra deposits analysed in this study have a more variable geochemical signature in comparison to those of their distal counterparts (Fig. 7). Although the proximal deposits record the detail in the eruption progression, the distal deposits tend to record the very largest phase of the eruption (e.g. Walker, 1980), but differences can be expected to occur also according to the azimuths of wind direction during an eruption and the number and degree of interconnectedness of magma bodies involved in the eruption (e.g. Walker, 1981;Kilgour and Smith, 2008;Sahetapy-Engel et al., 2014;Storm et al., 2014;Rubin et al., 2016). The tephrochronological principle is much more likely to utilise distal unknown deposits, and therefore we suggest that using the distal signature (or signatures) may be more appropriate for correlation in many studies. In general, distal tephras are more chemically homogeneous -but with some notable and well-documented exceptions -and this attribute therefore allows them to be traced over large areas (Manning, 1996). Alternatively, the identification of heterogeneity or bimodality in distal tephras, once recognised, can be an additional useful characteristic for fingerprinting (e.g. Shane et al., 2003aLowe et al., 2017). These statements, however, rely on the tephra being identified as a primary deposit, and not reworked. Reworking is commonly seen in paleofluvial deposits, for example those in the Whanganui Basin, and in other environments thin tephras are prone to mixing such as in surficial soils. This reworking can mix tephra from Figure 9. Results of PCA analysis on all TephraNZ major-element reference data for glass (normalised). Data are scaled to allow comparison within the PCA analysis. Tephra samples are coloured as per their source centre, and ellipses highlight the mean compositional region for each source caldera (KVC -Kapenga Volcanic Centre, MgVC -Mangakino Volcanic Centre, OVC -Okataina Volcanic Centre, RoVC -Rotorua Volcanic Centre, TVC -Taupō Volcanic Centre, and WVC -Whakamaru Volcanic Centre). PCA analysis was performed in R (see Supplement 1 for R script). Bar plots highlight the top elemental contributions for PC1 and PC2. The red dashed lines on the elemental contribution plots indicate the expected average contribution; if the contribution by each element were uniform, the expected value would be 1/no. of variables (e.g. 1/9 =∼ 11 %). Therefore, a variable with a contribution larger than this cut-off line (∼ 11 %) is considered important in contributing to the component. multiple eruptions and can cause highly variable glass chemistry within a single deposit (e.g. Shane et al., 2005Shane et al., , 2006. Fluvial reworking can be commonly identified by sedimentary structures within the deposit, for example, ripples or cross bedding indicative of fluvial transport and deposition (e.g. Shane, 1994;Schneider et al., 2001), over thickening of deposits (e.g. Vucetich and Pullar, 1969;Lowe, 2011), or through shard morphology, for example anomalously large shards or rounding of shards (e.g. Leahy, 1997).
Heterogeneous signatures (defined in approximation of where the standard deviation of the analyses is greater than the analytical error) in major-element compositions were identified for 10 of the tephra deposits: Kaharoa, Taupō Y5 proximal (P), Whakatāne-P, Hauparu, Maketū, Ngāmotu, Fordell, Onepuhi, Birdgrove, and Ototoka. Our data show that for some samples, specific trace elements and traceelement ratios have lower geochemical variability (Fig. 13a). The elements that work best to separate out the individual units within a deposit with a heterogeneous signature reflect the minerals that have formed during fractional crystallisation of the melt. Because of this, different elements or element ratios work for different tephras. For example, for glass from Kaharoa, Sr exhibits little variability (27-79 ppm), whereas for glass from Taupō, Sr compositional range exemplifies the heterogeneity in the sample (62-158 ppm; Fig. 13a).
Bimodality was identified for glass shards derived from four of the tephra horizons analysed: Rotorua (OVC),

Rerewhakaaitu (OVC), Poihipi (TVC), and Tahuna (TVC).
For all four of these, K 2 O concentration in glass exhibits bimodality, and therefore trace elements with similar chemical properties reinforce the bimodality (for example, LILEs Rb, Sr, and Cs; HFSEs Zr, or REE Eu), whereas most other trace elements do not show this bimodal signature (Fig. 13b).

Indistinguishable tephras
Euclidean similarity coefficient (ESC) analysis was used on all glass-shard reference data for tephras from Rotoiti/Rotoehu to Kaharoa in addition to the PCA and geochemical investigation to determine those samples that have indistinguishable element concentrations at similar ages (Table 3). Table 3a shows that similarities (similarity coefficient values (SC)) are observed in the major-element signatures of glass analyses for the following tephras: Waimihia and Unit K (SC = 0.11); Rotoma-P and Whakatāne-D KOT and Unit L (SC = 8.44) come up with significantly low (< 10) similarity coefficients for trace elements also, hence suggesting these samples will be indistinguishable in both major and trace elements. In addition, when simple geochemical assessment is applied, similarities are observed between glass analyses for Taupō and Waimihia; Mamaku and and Waiohau,Rotorua,and Rerewhakaaitu (Table 4). Table 5 outlines the key eruptions that show similar geochemical signatures in their glass chemistry and the ways in which they can be distinguished. Figure 14a shows that for Poihipi and Tahuna the best separation (although some overlap remains) is seen in the ratios La/Yb vs. Ba/Y; in addition, Tahuna also shows a bimodality in Ba/Th ratio which is not seen for Poihipi. For Rotoma and Mamaku, the tephras can be separated (although some overlap remains) using Ba/Th vs. Rb/Sr and Rb/Zr vs. Rb/Sr (Fig. 14b). Rotoma and Rotoehu/Rotoiti are very similar in their glass-shard major elements though they can be distinguished using specific, but a wide range of, trace elements (Fig. 14c). They are also very different in age and hence should not be too difficult to distinguish on the basis of stratigraphic positioning or dating.
Waimihia and Unit K (Taupō Subgroup) tephras are very difficult to distinguish, and their similar Late-to Mid-Holocene ages (3382 ± 50 and 5088 ± 73 cal yr BP, respectively; Lowe et al., 2013) and mineralogy could see them misidentified if dates were unavailable or imprecise. Geochemical investigation beyond the PCA and SC analyses of glass shows that Lu, Sc, Mn, and Co can be used to geochemically distinguish these two tephras (Fig. 14d), indicative of fractional crystallisation of differing amounts of clinopyroxene, plagioclase, and amphibole during the eruptive events. Although not identified by the SC analysis directly, Poronui (11 195 ± 51 cal yr BP) and Karapiti (11 501 ± 104 cal yr BP) tephras also have comparable age, geochemistry, and mineralogy; thus using major element, trace element, and traceelement ratios these two tephras remain indistinguishable. Glass shards from the three Holocene tephras, Waimihia, Poronui, and Karapiti, also have very similar trace element and trace-element ratios, but, as for Waimihia and Unit K, they can be distinguished with Lu, Sc, Mn, and Co, where Waimihia has higher Sc, Lu, and Mn but lower Co in compar- Figure 11. Major-element biplots to distinguish between caldera sources of tephras based on their glass major-element compositions (presented on a normalised basis; total iron is expressed as FeO): (a) and (b) show a comparison between all glass-shard analyses for the TVCand MgVC-sourced tephras; (c) and (d) indicate the distinction in glass compositional signature for the eruptives from the OVC and TVC that post-date the Kawakawa/Oruanui (KOT) eruption; (e) and (f) plots distinguish the glass analyses for tephra from the OVC into component eruptive time periods, with tephra from the Mangaone Subgroup (Jurado-Chichay and Walker, 2000;Smith et al., 2002) distinguishable from all other tephra from the OVC; (g) and (h) show the similarity in the geochemical compositions of glass of the tephras from the OVC and TVC for the eruptions prior to, and including, the KOT eruption. Colours are consistent for each caldera source; symbols are representative of different groups of tephra defined in the keys for each set of plots. ison to those of the Poronui and Karapiti tephras. They can also be distinguished simply with a biplot of FeO t vs. CaO, or Na 2 O+K 2 O or SiO 2 /K 2 O, or SiO 2 , where the Waimihia samples in general have lower FeO t , Na 2 O+K 2 O, and SiO 2 /K 2 O and higher CaO and SiO 2 in comparison to the equivalent values for Poronui and Karapiti samples (Fig. 14e, Table 4).
Geochemical investigation and PCA also highlight the similarity of the glasses of Waiohau, Rotorua, and Rerewhakaaitu tephras. There is added complexity with these samples as we have both proximal and distal deposits to compare, where, as discussed previously, the proximal samples will likely be more heterogeneous. Glass analyses of the Waiohau tephra show it can be distinguished from those for the Rotorua and Rerewhakaaitu tephras using a range of trace elements and trace-element ratios. In addition, the Rotorua and Rerewhakaaitu tephras are observed to be bimodal for some elements. The Waiohau tephra also has different mineralogy from that of Rotorua and Rerewhakaaitu tephras (Froggatt and Lowe, 1990;Lowe et al., 2008). Conversely, the Rotorua and Rerewhakaaitu tephras are indistinguishable in geochemistry and mineralogy, and therefore accurate dating and stratigraphic super-positioning would have to be re-lied upon to distinguish them with certainty (Fig. 14f, Table  4).
KOT, Okaia, and Unit L (Mangaone Subgroup) show indistinguishable major elements in their constituent glass shards and very similar trace elements. The TephraNZ samples have been compared with existing published data and are complementary with the existing data with respect to major elements (e.g. Sandiford et al., 2002;Smith et al., 2002Smith et al., , 2005Lowe et al., 2008;Allan et al., 2008;Molloy, 2008). This is the first time trace-element glass data have been published for Unit L and Okaia tephras. Our results show that Unit L glass shows bimodality in Rb/Zr, Ba/Th, Ce/Th, and Y/Th, and in this way, it can therefore be distinguished from the KOT and Okaia tephras (Table 5).

Proposed future research
This foundation dataset, derived in a formalised way, is unique in New Zealand and provides researchers with new avenues of research. It is our aim that the foundation dataset can be improved and expanded with analyses of other known deposits and that a subsidiary catalogue of accurately correlated geochemical analyses for such deposits can be added to  Froggatt and Lowe (1990). d Rotoehu tephra-fall volume in km 3 from Froggatt and Lowe (1990). e Unit L tephra-fall volume in km 3 from Jurado-Chichay and Walker (2000).
bolster the dataset. As noted earlier, it is beyond the scope of this paper to dive too deeply into the detail of the data, but we feel that it will provide the basis for countless projects in the future. Below we highlight some of the current gaps which we think would benefit from further research.

Further statistical analysis
We have applied simple ordination and statistical analyses to this dataset; however, we believe that further rigorous statistical analysis could be applied. Firstly, the analyses we present in this publication have been applied to mean values for each of the tephra samples (e.g. data from Table 2); there is no reason why these simple tests could not be applied to the full dataset, using all the individual (shard by shard) val- ues analysed for each sample. Secondly, we chose very basic tests (PCA and ESC) to fit with our requirements, but there is likely some more appropriate statistical test that could be applied to get the most out of this exceptional dataset. For example, (extended) canonical variate analysis (CVA): applying CVA to PCA results could determine optimal discrimination between multivariate data for single tephra deposits. This discrimination will increase the ability to identify an unknown tephra based on its similarity to known signatures plotted in multivariate space (e.g. discriminant function analysis; Tyron et al., 2009Tyron et al., , 2010Lowe et al., 2017;Bolton et al., 2020).

Whanganui Basin correlatives
A number of the tephras reported in this research were sampled from the Whanganui Basin, an uplifted Plio-Pleistocene basin margin sequence that preserves as many as 45 superposed cyclothems deposited since ∼ 3 Ma (Naish et al., 1996Naish and Kamp, 1997;Carter and Naish, 1998;Carter et al., 1999;Pillans, 2017;Grant et al., 2018Grant et al., , 2019. The tephra deposits within the basin contribute to the robust chronological framework that has been constructed for this region (Seward, 1976;Beu and Edwards, 1984;Alloway et al., 1993;Naish and Kamp, 1995;Saul et al., 1999;Pillans et al., 1994Pillans et al., , 2005Naish et al., 1996Rees et al., 2018Rees et al., , 2019Rees et al., , 2020Hopkins et al., 2021). These and other tephras (such as those derived from Tauranga Volcanic Centre) also record a critical time in New Zealand's volcanological history -the transfer between activity from the Coromandel Volcanic Zone to the Taupō Volcanic Zone (Briggs et al., 2005;Pittari et al., 2021). Deposits from this period are generally poorly exposed at source, and thus distal tephras could provide an insight into the eruptive history, geochemical evolution, and potentially even caldera evolution during this period (Houghton et al., 1995;Pittari et al., 2021). Most of the tephras reported in this research are well known and well dated, which is why they were included in the study. However, most do not have a known source caldera or source eruptives, or they have only been variably correlated to other deposits in New Zealand (e.g. Lowe et al., 2001;Pearce et al., 2008). There are also a number of tephra deposits in the Whanganui Basin that have yet to be studied, and thus a research project that is tephra focused, rather than using it as an accessory to a different line of enquiry, is timely.

IODP and ODP correlatives
At present there is a wealth of information that has yet to be fully investigated in the tephra record of the ODP Leg 181 Sites 1122, 1123, 1124, and 1125(Carter et al., 2003Alloway et al., 2005;Allan et al., 2008) and IODP Expedition 372 and 375 sites U1517 and U1520 (Pecher et al., 2018;Saffer et al., 2018). Pioneering work includes that undertaken by Watkins and Huang (1977) and Nelson et al. (1985), and findings from more "local" marine coring expeditions include those reported by Shane et al. (2006). The new reference material built by this project will allow more definitive identification and correlation of tephras within these cores, specifically post-2 Ma. However, the reports currently published on these deposits suggest that there are many more tephra deposits to be found in these marine and offshore sites than we have in the TephraNZ dataset (Carter et al., 2003;Alloway et al., 2005;Holt et al., 2010Holt et al., , 2011Hopkins et al., 2021). The TephraNZ dataset can provide a formalised correlation framework from which other unknown deposits can be determined, characterised, and in- Panel (a) shows analyses of glass from Kaharoa and Taupō tephras, both of which show a heterogeneous signature with most major elements (presented on a normalised basis). Sr has a low variability for Taupō but does not for Kaharoa tephra; conversely, Ba has a low variability for Kaharoa but does not for Taupō. Panel (b) shows the bimodal signature created for Tahuna tephra using K 2 O composition; this is also seen for Cs but not for Ba.
tegrated into a holistic tephrostratigraphic reconstruction. Allan (2008) and Allan et al. (2008) reported the major-and trace-element geochemistry of glass shards for tephra deposits dating from ∼ 1.65 Ma in the ODP 1123 core. They also give orbitally tuned ages for these tephras. However, of the 38 identified tephras only 7 were correlated to onshore equivalents. In addition, Alloway et al. (2005) reported over 100 tephra layers in the four ODP Leg 181 cores, dating back through orbital tuning (astrochronology) to 1.81 Ma. Using major-element chemistry of constituent glass shards, 13 tephras were correlated to equivalent onshore tephras including KOT, Omataroa, Rangitawa/Onepuhi, Kaukatea, Kidnappers-B and Kidnappers-A/Potaka, Unit D/Ahuroa, Ongatiti, Rewa, Sub-Rewa, Pakihikura, Ototoka, and Table  Flat. Analyses of glass from some of these are currently not in the TephraNZ database but could be easily added if Figure 14. Biplots for glass analyses for specific tephras which have very similar compositions and similar ages (see text for discussion and Table 5 for alternative elements). Plots show examples of the elements in glass that enable these tephras to be separated: (a) Poihipi and Tahuna (from TVC); (b) Mamaku and Rotoma-D (from OVC; note no trace-element data were obtained for Rotoma-P); (c) Rotoma and Rotoiti/Rotoehu (from OVC); (d) Waimihia and Unit K (from TVC); (e) Waimihia, Poronui and Karapiti -note that Poronui and Karapiti are indistinguishable using glass chemistry; and (f) Waiohau, Rerewhakaaitu and Rotorua -note that Rerewhakaaitu and Rotorua are indistinguishable using glass chemistry. All major-element data are presented on a normalised basis, and total iron is expressed as FeO.
the appropriate reference samples were available along with the capacity to analyse them. Alloway et al. (2005) reported an additional six tephra deposits that are correlated between the cores, but not to onshore equivalents, leaving potentially ∼ 81 tephra horizons within the ODP cores that are uncorrelated. The information that could be derived from their analysis would provide many details about the timing and evo-lution of the TVZ eruptions that are currently unobtainable from onshore deposits.

Mineral compositions
The TephraNZ reference dataset is only populated by glass major-and trace-element analyses at present. This is because glass geochemistry is one of the most frequently used and Table 5. Dominant ferromagnesian. Mineral assemblages for late Quaternary silicic tephra deposits updated from Froggatt and Lowe (1990). Plag -plagioclase feldspar; opx -orthopyroxene; mnt -magnetite; ilm -ilmenite; hyp -hypersthene; hbl -hornblende; bio -biotite; cgtcummingtonite; aug -augite. The mineral assemblages are listed with mineral species in order of abundance; the diagnostic mineral in each assemblage is in bold. Tephras are listed multiple times if their mineral assemblage changes through the eruption sequence, and deposits in brackets are not included in the TephraNZ database.
Assemblage 1 a Assemblage 2 Assemblage 3 Assemblage 4 Assemblage 5 Assemblage 6 Plag ± opx ± mnt ± ilm Hyp + hbl ± aug Hyp + hbl + bio Hyp + cgt ± hbl Hyp + aug ± hbl Aegirine b  Barker et al. (2015). b Assembly 6 aegirine ± riebeckite ± aenigmatite ± olivine ± tuhualite. c Motutere was listed in Froggatt and Lowe (1990) as a single unit with Assemblage 1 mineralogy, but this has subsequently been redefined by Wilson (1993) into two subunits G and H, which do not have their independent assemblages defined. accessible tools for tephra correlation. Aerodynamic sorting of tephra componentry through transportation adds to the favourability of glass shards as the dominant tool because glass shards tend to be the only phase that is found at both proximal and distal sites. However, previous New Zealand-based studies have specified how mineral assemblages and their geochemical compositions can be used to distinguish certain tephras and their source (e.g. Nairn and Kohn, 1973;Lowe, 1988;Froggatt and Lowe, 1990;Froggatt and Rogers, 1990;Shane, 1998;Shane et al., 2003b;Allan et al., 2008;Lowe et al., 2008;Lowe, 2011). For example, the mineral cummingtonite, where predominant, is a known identifier for tephras from the Haroharo complex of the OVC (Whakatāne, Rotoma, Rotoehu/Rotoiti (Table 5); Ewart, 1968;Lowe, 1988;Froggatt and Lowe, 1990). At present, ferromagnesian mineralogical assemblages (following Froggatt and Lowe, 1990;Smith et al., 2005; for all the TephraNZ samples younger than and including Rotoehu/Rotoiti have been published (see Table 5). Extending this tabulation to include the older samples would add another useful criterion to the correlation toolbox for tephras containing ferromagnesian minerals.
Additionally, the fractional crystallisation of plagioclase, biotite, amphibole, zircon, hydrous mineral phases, or Fe-Ti oxides has been shown to be the key impactor on the traceelement chemistry (Shane, 1998;Allan, 2008;Turner et al., 2009Turner et al., , 2011. Thus the prevalence of these minerals is also an important potential fingerprinting tool. The information on the mineralogy of the tephras is not only useful for fingerprinting but also can be used in determining the characteristics of the magma source components and potentially provide estimates for the temperature, pressure, and oxidation states of the magmatic system before eruption (e.g. Lowe, 1988Lowe, , 2011Shane, 1998). Thus, this information can allow hypotheses to be developed on the reactivation and triggering of these large-scale eruptions, an important step for hazard and risk monitoring.

The New Zealand tephra "Bermuda Triangle"
At present the TephraNZ database is very well populated for samples from the Rotoiti/Rotoehu through to Kaharoa eruption. It also has a high number of samples, but not an exhaustive list, from Mamaku ignimbrite (∼ 0.22-0.23 ka) to the Hikuroa Pumice (2 Ma). There is a stark deficit in tephras between the Rotoiti/Rotoehu eruption and Mamaku ignimbrite (Table 1). This ∼ 150 kyr gap in the volcanic record (∼ 220 to 45 ka) is intriguing as there is proximal evidence for activity during this period. For example, Rosenberg et al. (2020) reported the occurrence of volcanic formations in cores forms the Taupō region in the age range of ∼ 168 to 92 ka, including the Huka Falls formations, Racetrack rhyolites, and the Te Mihi rhyolites. Tephra deposits, in some cases strongly weathered successions of multiple units broadly lumped together as a "formation", such as the so-called Hamilton Ash Formation, have been reported during this time period both terrestrially and in marine and lacustrine sediment cores (Ward, 1967;Pain, 1975;Vucetich et al., 1978;Iso et al., 1982;Froggatt, 1983;Manning, 1996;Lowe et al., 2001;Newnham et al., 2004;Allan et al., 2008;Briggs et al., 2006;Lowe, 2019;Benjamin Laeuchli, personal communication, 2020). However, at present the authors are not aware of a detailed, up-to-date study into the primary compositions of these tephra deposits. The key deposits identified during this time period include (but are not limited to) Kaingaroa Ignimbrite (∼ 0.18 Ma;Froggatt, 1983), Tablelands Tephra Formation (∼ 0.21-0.18 Ma; Iso et al., 1982, 0.39-0.34 Ma;Manning, 1996), Hamilton Ash Formation (0.34-0.125 Ma; Lowe, 2019), Kutarere tephra (= Mamaku ignimbrite 0.22-0.23 Ma; Shane et al., 1994;Houghton et al., 1995;Tanaka et al., 1996;Milner et al., 2003), Kukumoa Subgroup (∼ 0.22-0.05 Ma; Manning, 1996), and Tikotiko Ash (∼ 0.125 ka; Lowe, 2019). A number of these studies are outdated, and with improved methodologies (major-and trace-element analysis, potentially of melt inclusions where preserved, dating techniques, and other measures to help construct time frames such as via phytolith studies to determine glacial vs. interglacial periods, and potentially also paleomagnetic measures as shown for the much older and very strongly weathered Kauroa Ash Formation: Hopkins et al., 2021) it could be timely to further investigate this period of (apparent) deficit.

Conclusions
Major-and trace-element geochemical compositions of glass shards for a large suite of prominent, widespread New Zealand rhyolitic tephras have been analysed systematically and published for the first time as TephraNZ. TephraNZ is a foundation dataset for collating geochemical data about New Zealand tephras based on analyses of their glass components. The foundation reference dataset is made up of known deposits that have their ages quantified through independent methods, and/or are from the type sites where tephras were first defined, or well-documented reference sections. Detailed methodology is reported to allow subsequent research to acquire comparable data to those in this database. Principal component analysis of the glass geochemistry indicates that for the TephraNZ foundation dataset, as a whole, major elements Al, K, Si, Ti, Fe, Mn, and Cl are responsible for the spread along PC1 and PC2 space. When the trace elements are run together with the major elements, V, Co, Mg, Cu, Ti, Sr, Sc, Ca, Cs, Zr, Th, and Rb are most responsible for the separation in PC1 and PC2 space. Euclidean similarity coefficients can also be used to distinguish between some geochemically similar glass analyses. However, further detailed geochemical investigation is required to distinguish others. Geochemically indistinguishable tephras (on the basis of both major-and trace-element glass-shard compositions) are identified as Taupō and Waimihia; Poronui and Karapiti; Rotorua and Rerewhakaaitu; and KOT and Okaia. Only Poronui and Karapiti are noted as entirely indistinguishable, with other methods of characterisation listed as alternative options, including mineralogy, age, and stratigraphic relationships.
Data availability. All the data provided in this article are available as Excel files in the Supplement. The data are also available from GNS Science, New Zealand, at Pet Lab (https://pet.gns.cri. nz, GNS Science, 2004), and as a file submission on EarthChem (https://doi.org/10.26022/IEDA/111724, Hopkins et al., 2020b).
Author contributions. JLH and RJW designed the project. DJL and BJP contributed samples from previous field campaigns, and DJL provided guidance on new and existing field locations for sample collection. JLH and JEB undertook the fieldwork, lab work, analysis, and data reduction. ABHR advised on statistical analysis and R coding. LA supervised and helped JEB develop LA-ICP-MS analysis and data reduction. FT supervised and helped JEB develop sample mounting and polishing procedures. JLH wrote the manuscript with contributions from all co-authors.
Britta Jensen, Maxim Portnyagin, Stephen Kuehn, and an anonymous reviewer for their detailed and helpful comments and reviews in the development of this article and the dataset.
Financial support. This research has been supported by the Royal Society of New Zealand (Royal Society Te Apārangi) (grant nos. MFP-VUW1809 and UOW1006) and the Victoria University of Wellington Summer Scholarship programme (grant no. 136 (2019-2020)).
Review statement. This paper was edited by Britta Jensen and reviewed by Maxim Portnyagin, Stephen Kuehn, and one anonymous referee.