Articles | Volume 7, issue 2
https://doi.org/10.5194/gchron-7-199-2025
https://doi.org/10.5194/gchron-7-199-2025
Research article
 | 
05 Jun 2025
Research article |  | 05 Jun 2025

A comparison between in situ monazite Lu–Hf and U–Pb geochronology

Alexander T. De Vries Van Leeuwen, Stijn Glorie, Martin Hand, Jacob Mulder, and Sarah E. Gilbert
Abstract

In complex metamorphic rocks, monazite U–Pb dates can span a wide concordant range, leading to ambiguous geological interpretations (e.g. slow protracted cooling versus multiphase growth). We present in situ monazite Lu–Hf analysis as an independent chronometer to verify U–Pb age interpretations. Monazite Lu–Hf dates were attained via laser ablation inductively coupled plasma mass spectrometry equipped with collision/reaction cell technology (LA-ICP-MS/MS). In situ Lu–Hf dates for potential reference monazites with uncertainties <1.6 % agree with published U–Pb dates, validating the approach. We demonstrate the method on complex metamorphic samples from the Arkaroola region of the northern Flinders Ranges, South Australia, which exhibit protracted thermal and monazite growth histories due to high geothermal gradient metamorphism. In situ Lu–Hf dates reproduce the main U–Pb monazite age populations, demonstrating the ability to reliably resolve multiple age populations from polymetamorphic monazite samples.

Share
1 Introduction

Monazite is a common accessory mineral in a broad range of metamorphic and felsic igneous rocks and forms across wide-ranging pressure–temperature conditions. In metamorphic rocks, monazite can record multiple stages of crystal growth (e.g. Kohn and Malloy, 2004; Rubatto et al., 2013), undergo fluid-mediated dissolution–precipitation reactions (e.g. Harlov et al., 2011; Seydoux-Guillaume et al., 2002), and at high temperatures and/or strain rates undergo recrystallisation (e.g. Erickson et al., 2015; Kelly et al., 2012). This responsiveness to changing physicochemical conditions makes monazite amenable to recording multiple overprinting events and complex episodes of fluid–rock interaction. Consequently, U–Pb dating of monazite has become routine for deciphering the timing and tempo of thermal events in crustal rocks (e.g. Kohn and Malloy, 2004; Larson et al., 2022; Parrish, 1990; Rubatto et al., 2001).

Widely dispersed concordant dates are a common observation in monazite U–Pb data from metamorphic rocks in complex and/or long-lived orogens (e.g. Clark et al., 2024; De Vries Van Leeuwen et al., 2021; Kirkland et al., 2016; Korhonen et al., 2013). There is often ambiguity surrounding what this dispersion represents. Common interpretations consider prolonged, slow cooling and associated volume diffusion or partial dissolution–reprecipitation by overprinting or prolonged thermal events. Detailed microstructural observations and trace element geochemistry play a key role in contextualising these data; however, in their absence or ambiguity, the significance of dispersion in U–Pb dates can be difficult to interpret. As such, it is important to understand the significance of such concordia dispersion, as it can lead to substantially different tectonic interpretations.

With the recent advent of in situ Lu–Hf dating facilitated by laser ablation inductively coupled plasma mass spectrometry equipped with collision/reaction cell technology (LA-ICP-MS/MS), a new frontier of in situ dating opportunities has emerged (e.g. Glorie et al., 2023; Simpson et al., 2021, 2022; Yu et al., 2024). In this contribution, we first appraise in situ Lu–Hf isotopic data from monazite reference materials and in-house secondary reference materials by comparing calculated Lu–Hf dates with published U–Pb dates. We then compare the results of in situ Lu–Hf and U–Pb geochronology from monazites that record a protracted history of fluid-driven dissolution and re-precipitation. Monazite Lu–Hf dating by LA-ICP-MS/MS was recently demonstrated using an iCap TQ instrument (Wu et al., 2024). However, this instrumental approach lacks axial ion acceleration and the ability to set a wait time between isotope jumps. These limitations affect sensitivity and induce undesirable isobaric interferences, hindering the exploration of in situ monazite Lu–Hf dating from reaching its full potential. Here we present monazite Lu–Hf data acquired using an Agilent 8900x mass spectrometer, with demonstrated better performance for heavy ions, and show that even in complex systems with protracted thermal histories, monazite Lu–Hf dating yields robust geochronometric data that can be used to interrogate U–Pb dates. In situ Lu–Hf dating of monazite can resolve multiple age populations from single grains and thus may be useful in cases where the U–Pb system has been compromised by Pb loss, non-radiogenic Pb contamination, excess 206Pb due to 230Th uptake, low U concentration, or a combination of these factors.

2 Methods

2.1 Lu–Hf geochronology and trace element geochemistry

Monazite Lu–Hf geochronological and trace element analysis was undertaken at Adelaide Microscopy, at the University of Adelaide, following Simpson et al. (2021), which we briefly outline here. Analyses of Lu–Hf were acquired across two sessions using a RESOlution-LR 193 nm excimer laser ablation system coupled to an Agilent 8900 ICP-MS/MS. The reaction gas used was NH3, supplied as a mixture of 10 % NH3 in 90 % He. Laser beam diameters were set to either 43 or 67 µm, depending on Lu concentrations and microstructural constraints (e.g. size and shape of monazite compositional domains). The laser repetition rate was 10 Hz with an average on-sample fluence of ∼3.5 J cm−2. The ablated sample material was transported from the laser cell to the ICP-MS by a He carrier gas (380 mL min−1). Data acquisition consisted of (1) 30 s of baseline acquisition; (2) 40 s of continuous ablation, during which data were collected; and (3) ∼25 s of washout. The following isotopes (mass shifts denoted in parentheses) were measured: 27Al, 43Ca, (47+66)Ti, 88Sr, (89+83)Y, (90+83)Zr, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, 157Gd, 159Tb, 163Dy, 165Ho, 166Er, 169Tm, 172Yb, 175Lu, (175+82)Lu, (176+82)Hf, (178+82)Hf, and (232+15)Th. Axial acceleration was set to 2 V, and a wait time offset of 2 ms was set to avoid memory effects when cycling between isotopes. 175Lu was measured as a proxy for 176Lu and 178Hf as a proxy for 177Hf. The calculation of 176Lu and 177Hf was performed assuming present-day 176Lu/175Lu (0.02659) and 177Hf/178Hf (0.682) ratios (De Biévre and Taylor, 1993), following the procedures outlined in Simpson et al. (2021). Isobaric interference from (176+82)Lu on (176+82)Hf was corrected by monitoring (175+82)Lu and subtracting a proportion of this signal from (176+82)Hf based on the present-day 176Lu/175Lu ratio. No corrections were performed for isobaric interferences from (176+82)Yb on (176+82)Hf, as this has been demonstrated to be negligible (Simpson et al., 2021).

Data reduction was performed in LADR (Norris and Danyushevsky, 2018). Background-subtracted isotopic ratios were normalised to NIST 610 glass using the Nebel et al. (2009) isotope dilution multi-collector inductively coupled plasma mass spectrometry (ID-MC-ICP-MS) isotopic compositions of 176Lu/177Hf=0.1379±0.005 and 176Hf/177Hf=0.282111±0.000009. Analyses of NIST 610 were conducted before and after every 40 unknown analyses and were also used to normalise isotopic ratios and correct for instrument drift. No downhole fractionation corrections were applied, as there was no observable downhole fractionation. This is consistent with the results of Simpson et al. (2021), where no downhole fractionation was observed in garnet, apatite, or xenotime, using laser beam diameters between 43 and 120 µm. A subsequent matrix fractionation correction was applied to the calculated 177Lu/176Hf ratios (see Simpson et al., 2021, 2023). Although matrix-matched reference materials are desirable, it has been demonstrated that when using the same method and instrumentation as that outlined in Simpson et al. (2021), correction factors for materials with similar ablation characteristics analysed with the same laser beam conditions are indistinguishable (e.g. Glorie et al., 2023, 2024a). Here, the Bamble-1 and OD-306 apatite reference materials were employed to perform matrix fractionation corrections for sessions 1 and 2, respectively. Bamble-1 (1102±5 Ma; Simpson et al., 2024) yielded an uncorrected inverse Lu–Hf isochron age of 1150±8 Ma (n=20, mean square weighted deviation (MSWD) = 2.0, p=0.00), and OD-306 (1597±7 Ma; Thompson et al., 2016) yielded an uncorrected inverse Lu–Hf isochron age of 1671±15 Ma (n=25, MSWD = 0.94, p=0.55). This resulted in correction factors of 4.40±0.04 % and 4.71±0.05 % for sessions 1 and 2, respectively. Monazite reference materials TS-Mnz (Budzyń et al., 2021) and RW-1 (Ling et al., 2017) were analysed in both sessions to appraise the accuracy of these corrections (discussed below). Additionally, in-house apatite secondary reference material HR-1 (long-term Lu–Hf age of 344±2 Ma; Glorie et al., 2024a) was analysed across both sessions and yielded corrected isochron ages of 348±4 Ma (n=15, MSWD = 1.30, p=0.17) and 342±3 Ma (n=26, MSWD = 0.96, p=0.52) for sessions 1 and 2, respectively. Trace element data were calibrated using NIST 610, with Ce used as the internal standardisation element for calibration. Ce concentrations were set to 21.39 wt % for RW-1, 21.42 wt % for TS-Mnz, and 20 wt % for all other samples. Trace element concentrations were quantified by normalising the weight percent of oxides to 100 % totals.

Inverse Lu–Hf isochron and weighted mean ages were calculated using IsoplotR (Vermeesch, 2018), using a 176Lu decay constant of 0.00001867±0.00000008 Myr−1 (Söderlund et al., 2004). Given the narrow range of initial terrestrial 177Hf/176Hf ratios, anchored regressions were used to calculate inverse isochrons (following the approach of Glorie et al., 2024a, b). All inverse isochron plots presented in this study were anchored to an initial 177Hf/176Hf value of 3.55±0.05, covering the range of plausible terrestrial possibilities (Spencer et al., 2020). This avoids issues which may be encountered when calculating regressions on samples with low 177Hf/176Hf variability, as samples with little spread along the isochron can lead to spurious upper intercepts yielding geologically implausible initial Hf values (e.g. Vermeesch, 2024). The algorithm employed for performing anchored regressions is detailed in Vermeesch (2024). Both analytical and propagated uncertainties are presented following the format t±x [y] Ma, where t is the calculated Lu–Hf date, x the analytical 2 SE uncertainty, and y the propagated uncertainty. Error propagation involved the quadratic addition of uncertainties on the measured sample date, measured mineral reference material date, the known reference material age, the 176Lu decay constant, and the uncertainty associated with matrix fractionation correction. Uncertainties are reported at the 2 SE level unless the quoted p value is <0.05; then the quoted uncertainty accounts for overdispersion following the method outlined in Vermeesch (2018).

2.2 U–Pb geochronology and trace elements

Monazites were analysed in situ by spot-targeting guided by back-scattered electron (BSE) images collected on an FEI Quanta 450 scanning electron microscope (SEM) housed at Adelaide Microscopy, the University of Adelaide. U–Pb and trace element data were collected using a RESOlution-LR 193 nm excimer laser ablation system coupled to an Agilent 8900 ICP-MS/MS at Adelaide Microscopy, the University of Adelaide. Ablation was performed with a laser frequency of 5 Hz employing a 13 µm laser beam diameter with an average on-sample fluence of ∼2.2 J cm−2. The ablated sample material was transported from the laser cell to the ICP-MS by a He carrier gas (380 mL min−1). Data acquisition consisted of (1) 30 s of baseline acquisition; (2) 30 s of continuous ablation, during which data were collected; and (3) ∼25 s of washout.

The isotopes collected were 29Si, 31P, 43Ca, 89Y, 90Zr, 139La, 140Ce, 141Pr, 146Nd, 147Sm, 153Eu, 157Gd, 159Tb, 163Dy, 165Ho, 166Er, 169Tm, 172Yb, 175Lu, 202Hg, 204Pb, 206Pb, 207Pb, 208Pb, 232Th, and 238U. MAdel was used as the primary reference material to correct for elemental fractionation and mass bias (Payne et al., 2008), with 94-222 as a secondary reference material to monitor precision and accuracy (Maidment, 2005). Reference materials were analysed after every 10–15 unknown analyses. For trace element concentrations, NIST 610 (Pearce et al., 1997) was analysed after every 10–15 unknown analyses. U–Pb isotope and trace element data were reduced using the “U–Pb Geochronology” and “Trace Elements” data reduction schemes in Iolite version 4.9.3 (Paton et al., 2011), respectively. Trace element data were calibrated using NIST 610. The internal standard element used was Ce, and trace elements were quantified by normalising the weight percent of oxides to 100 % totals. Error propagation and uncertainty reporting follow the same approach discussed for Lu–Hf data. Secondary reference materials yielded results comparable to published values, with 94-222 yielding a weighted mean 206Pb/238U age of 447±1 Ma (n=41, MSWD = 1.1, p=0.36), within 2 SE uncertainty of the reference age of 450.2±3.4 Ma (Maidment, 2005). Weighted means and concordia plots were generated using IsoplotR (Vermeesch, 2018).

3 In situ Lu–Hf geochronology of candidate monazite reference materials

3.1 RW-1

Ten ∼1 mm fragments of RW-1 mounted in a 25 mm epoxy resin disk were analysed in this study. The crystals are reddish-brown in colour and free of inclusions and cracks. RW-1 is a high-Th monazite that originates from a 20–30 m wide pegmatite dyke located in the Landsverk 1 quarry in the Evje–Iveland district, south Norway (Ling et al., 2017). Ling et al. (2017) present U–Pb isotope dilution thermal ionisation mass spectrometry (ID-TIMS)/ID-MC-ICP-MS isotopic data. These authors recommend the mean 207Pb/235U age of 904.15±0.26 Ma (95 % confidence) as the best estimate for the crystallisation age of the pegmatite hosting the RW-1 monazite (Ling et al., 2017). Electron probe microanalysis (EPMA) compositional data show that RW-1 has a Ce2O3 content of 25.22 wt %, Nd2O3 of 14.47 wt %, ThO2 of 13.5 wt %, Y2O3 of 2.44 wt %, and UO2 of 0.30 wt % (Ling et al., 2017). Additional LA-ICP-MS data show a Lu content of 27±5 (2σ) ppm (Ling et al., 2017), making the sample amenable to in situ Lu–Hf geochronology.

In this study, RW-1 yields inverse Lu–Hf isochron dates of 907±10 [16] Ma (n=26, MSWD = 1.8, p<0.00) and 907±11 [17] Ma (n=30, MSWD = 2.0, p<0.00) and single-spot weighted mean dates of 905±10 [16] Ma (n=26, MSWD = 1.9, p<0.00) and 909±8 [16] Ma (n=29, MSWD = 1.2, p=0.19) for sessions 1 and 2, respectively (Figs. 1a, 2a). Common Hf is very low, with most analyses yielding 177Hf/176Hf<0.05. These dates are within the uncertainty of published U–Pb ID-TIMS/ID-MC-ICP-MS ages (Ling et al., 2017).

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f01

Figure 1Inverse isochron plots for (a) RW-1, (b) TS-Mnz, (c) Storø, and (d) Pilbara. Purple-coloured ellipses correspond to analyses from session 1, and orange-coloured ellipses correspond to analyses from session 2. Individual data-point uncertainties are 2 SE.

Download

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f02

Figure 2Single-spot weighted mean plots for (a) RW-1, (b) TS-Mnz, (c) Storø, and (d) Pilbara. Purple-coloured bars correspond to analyses from session 1, and orange-coloured bars correspond to analyses from session 2. For TS-Mnz and Storø, darker-coloured bars correspond to analyses collected using a 67 µm spot size, while lighter-coloured bars were collected using a 43 µm spot. Individual data-point uncertainties are 2 SE.

Download

3.2 TS-Mnz

A single ∼7 mm fragment of TS-Mnz mounted in a 25 mm epoxy resin disk was analysed in this study. The crystal is reddish-brown in colour with abundant cracks that host thorite inclusions. These cracks were avoided during analysis, with only fresh monazite being analysed. The crystal, originally attained from a mineral dealer, likely originates from the Arendal region of Norway (Budzyń et al., 2021). U–Pb ID-TIMS data yield a 207Pb/235U age of 910.42±0.34 Ma (2σ) (Budzyń et al., 2021). EPMA shows that TS-Mnz has a Ce2O3 content of 25.09 wt %, Nd2O3 of 15.92 wt %, ThO2 of 4.80 wt %–9.44 wt %, Y2O3 of 2.83 wt %, and UO2 of 0.16 wt %–0.29 wt % (Budzyń et al., 2021). LA-ICP-MS data also presented in Budzyń et al. (2021) indicate that TS-Mnz has a Lu content of 28.2±3.7 (2σ) ppm, making the sample amenable to in situ Lu–Hf geochronology.

In this study, TS-Mnz yields inverse Lu–Hf isochron dates of 915±10 [16] Ma (n=26, MSWD = 2.0, p<0.00) and 912±9 [16] Ma (n=26, MSWD = 1.6, p=0.03) and single-spot weighted mean dates of 913±11 [16] Ma (n=26, MSWD = 2.2, p<0.00) and 911±9 [16] Ma (n=26, MSWD = 1.6, p=0.02) for sessions 1 and 2, respectively (Figs. 1b, 2b). In session 1, the laser beam diameter was varied between 43 and 67 µm; however, aside from comparatively smaller uncertainties on analyses collected with the larger 67 µm spot size, no discernible difference in the accuracy of calculated Lu–Hf dates was observed (Fig. 2b). Common Hf is very low, with most analyses yielding 177Hf/176Hf<0.05. These dates are within the uncertainty of published U–Pb ID-TIMS ages (Budzyń et al., 2021).

3.3 Storø

A total of 35 monazite grains ranging from 30 to 170 µm mounted on two epoxy resin disks were analysed in this study. The grains are yellow in colour with few cracks and inclusions. This sample originates from the Storø Quartzite in West Greenland (Gardiner et al., 2023). Existing laser ablation split-stream ICP-MS data yield concordant U–Pb monazite dates between 2600 and 2630 Ma with an overdispersed concordia age of 2619±8 Ma; the authors estimate the crystallisation age of monazite in this sample to be ca. 2620 Ma (Gardiner et al., 2023).

In this study, Storø yields an inverse Lu–Hf isochron date of 2619±9 [35] Ma (n=44, MSWD = 1.4, p=0.06; Fig. 1c) and a single-spot weighted mean date of 2618±11 [36] Ma (n=44, MSWD = 1.4, p=0.04; Fig. 2c). The laser beam diameter was varied between 43 and 67 µm on a subset of analyses from Storø. Similar to the data from TS-Mnz in session 1, there was no discernible difference in the accuracy of calculated Lu–Hf dates observed between the two spot sizes (Fig. 2c). Common Hf is very low, with most analyses yielding 177Hf/176Hf<0.05. These dates are within the uncertainty of the published U–Pb (LASS)-ICP-MS age of 2619±8 Ma (Gardiner et al., 2023).

3.4 Pilbara

Ten ∼1 mm monazite fragments mounted in a 25 mm epoxy resin disk were analysed in this study. The grains are reddish-brown in colour with few cracks and inclusions. Originating from a granitoid suite in the Pilbara Craton, Western Australia, this sample belongs to the Mawson Collection housed at the University of Adelaide. U–Pb dating of this sample yields an approximate age of ca. 2870 Ma (unpublished).

Pilbara yields an inverse Lu–Hf isochron date of 2836±19 [42] Ma (n=19, MSWD = 1.1, p=0.39; Fig. 1d) and a single-spot weighted mean date of 2840±19 [42] Ma (n=17, MSWD = 0.8, p=0.73; Fig. 2d). Common Hf is very low, with most analyses yielding 177Hf/176Hf<0.05.

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f03

Figure 3(a) Geological setting of the Arkaroola region in South Australia, Australia; (b) simplified geological map of the Arkaroola region.

4 Comparing U–Pb and Lu–Hf data from the Arkaroola region

4.1 Geological background

The Arkaroola region of the northern Flinders Ranges, South Australia, hosts some of the highest-heat-producing basement rocks exposed at Earth's surface (De Vries Van Leeuwen et al., 2021; McLaren et al., 2006). These basement rocks comprise Mesoproterozoic granitoids and metasedimentary rocks which are exposed in two inliers, the Mount Painter and Mount Babbage inliers (Fig. 3). Overlying these high-heat-producing basement rocks is a 12–15 km succession of sedimentary rocks which form the Adelaidean stratigraphy of the Adelaide Superbasin (Lloyd et al., 2020; Paul et al., 1999; Preiss, 2000). Deposition of these sedimentary rocks began in the early-to-middle Neoproterozoic and terminated in the early Cambrian (Lloyd et al., 2020; Preiss, 2000).

The accumulation of this thick sedimentary package on high-heat-producing basement rocks leads to the development of steep thermal gradients, resulting in high-temperature, low-pressure metamorphism of the basal portion of the sedimentary succession (De Vries Van Leeuwen et al., 2021; McLaren et al., 2006). This is recorded by the development of cordierite–biotite-bearing assemblages in metapelitic rocks (Fig. 3; De Vries Van Leeuwen et al., 2021; Mildren and Sandiford, 1995). The consequence of this style of high-heat-production-driven “burial” metamorphism is that high-thermal-gradient conditions will persist, providing the rocks are sufficiently deep.

The cordierite-bearing metapelitic rocks at the base of the Adelaidean stratigraphy at Arkaroola record two distinct periods of monazite growth at ca. 580–540 Ma and ca. 450–400 Ma (De Vries Van Leeuwen et al., 2021), which are interpreted to reflect the timing of thermally and hydrothermally catalysed monazite growth. The ca. 580–540 Ma monazite population corresponds to a significant interval of subsidence and sedimentation in the Adelaide Superbasin, where ∼5 km of sediment was deposited between ca. 580–520 Ma (Paul et al., 1999; Preiss, 2000), significantly increasing the burial depth of the basement. Sediment accumulation was associated with the formation of progressively younger monazite ages for incipient metamorphism up-stratigraphy (De Vries Van Leeuwen et al., 2021). The second monazite population at ca. 450–400 Ma is more enigmatic, as sedimentation in the Adelaide Superbasin had terminated by the onset of the ca. 520–490 Ma Delamerian Orogeny (e.g. Foden et al., 2006; Preiss, 2000). However, evidence exists for a significant, regionally widespread hydrothermal–magmatic event between ca. 460–400 Ma (e.g. Elburg et al., 2013; McLaren et al., 2006). Monazite also exhibits increasing heavy rare earth element (HREE) + Y contents between the ca. 580–540 Ma and ca. 450–400 Ma populations (De Vries Van Leeuwen et al., 2021), suggesting the thermal maximum was attained at ca. 400 Ma, supporting the notion that increasing temperatures were a function of increasing burial depth over at least a ∼150 Myr period (De Vries Van Leeuwen et al., 2021).

4.2 Sample descriptions

Two metapelitic samples, ARK 2017-11 and ARK 2017-15, were collected from the Paralana Quartzite, which forms the basal portion of the Adelaidean stratigraphy and occupies the unconformable interface with the high-heat-producing Mesoproterozoic basement rocks of the Mount Painter Inlier (Fig. 3). These samples, previously described in De Vries Van Leeuwen et al. (2021), were derived from discrete metapelitic layers within broadly psammitic-to-quartzitic packages of the Paralana Quartzite. Although mineral modes vary between these two samples, both are mineralogically similar, exhibiting large porphyroblasts (up to 1 cm) wrapped by a strong fabric defined by biotite, plagioclase, and minor quartz. These porphyroblasts comprise fine-grained intergrowths of plagioclase and K-feldspar, biotite, and hematite and are interpreted to represent altered cordierite.

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f04

Figure 4BSE images of monazite grains from samples ARK 2017 15 and ARK 2017 11. Lu–Hf spots are represented by solid-lined circles and are coloured according to their microstructural domains (mnz1: blue; mnz2: green), corresponding to the colour scheme in Fig. 5. U–Pb spots are represented by dash-lined shaded circles and are coloured according to the corresponding age populations shown in Fig. 6. White-coloured dashed circles correspond to U–Pb analyses that display isotopic mixing. U–Pb spot numbers are in a smaller font size and underlined. Grains were re-polished between Lu–Hf and U–Pb analyses; as such, some U–Pb spots were placed beyond the extent of these BSE images.

Download

Monazite in these samples predominantly exists as large (up to ∼500µm) foliation-parallel elongate grains, anhedral grains throughout the matrix, or inclusions within altered cordierite porphyroblasts. BSE images reveal two distinct generations of monazite (Fig. 4). The first generation, mnz1, forms dark BSE response poikiloblastic cores containing rounded inclusions of quartz and rare hematite (Fig. 4). Mnz1 often exhibits chaotic zoning patterns with high-Th monazite intergrowths, with some grains in sample ARK 2017-11 also exhibiting patchy zoning with no clear core–rim relationship (Fig. 4l). The second generation, mnz2, forms as brighter BSE response rims mantling mnz1 or grains with no core–rim relationships (Fig. 4). These domains are inclusion-poor and can be homogeneous or exhibit patchy or wispy zoning patterns (Fig. 4). All analysed grains exhibit embayed margins.

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f05

Figure 5(a, c) Tera–Wasserburg concordia plots for U–Pb analyses from samples (a) ARK 2017 15 and (c) ARK 2017 11; (b, d) chondrite-normalised REE plots for analyses from samples (b) ARK 2017 15 and (d) ARK 2017 11. Blue ellipses and lines correspond to analyses from mnz1, and green ellipses and lines correspond to analyses from mnz2. Grey and unfilled ellipses in panels (a) and (c) and grey lines in panels (b) and (d) represent isotopically mixed analyses or erroneous analyses. Individual data-point uncertainties are 2 SE.

Download

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f06

Figure 6(a, c) Inverse isochron plots for samples (a) ARK 2017 15 and (c) ARK 2017 11; (b, d) chondrite-normalised REE plots for analyses from samples (b) ARK 2017 15 and (d) ARK 2017 11. Blue ellipses and lines correspond to analyses from mnz1, and green ellipses and lines correspond to analyses from mnz2. Unfilled ellipses in panels (a) and (c) and grey lines in panels (b) and (d) represent mixed analyses which were not considered for age calculations. Individual data-point uncertainties are 2 SE.

Download

4.3 In situ U–Pb and Lu–Hf geochronology

A total of 42 U–Pb spot analyses were collected from sample ARK 2017-15, 7 of which belong to mnz1 and 34 to mnz2 (Fig. 5a). An additional analysis, which yields a concordant 206Pb/238U date of 487±13 Ma, is interpreted to reflect isotopic mixing between mnz1 and mnz2 domains (Fig. 5a). Chondrite-normalised rare earth element (REE) data help to delineate data from mnz1 and mnz2 domains, with monazite belonging to the mnz2 population consistently showing elevated HREE contents (Fig. 5b). Given the large spread of dates in both the mnz1 and mnz2 populations, the range of dates within each population is the preferred method of assigning an “age” to each population. However, given that overdispersed dates often reflect underlying processes, we present them in Fig. 5 for completeness. Analyses from mnz1 yield 206Pb/238U dates of 604–571 Ma, whereas those from mnz2 are spread between 444 and 390 Ma (Fig. 5a). These dates replicate the previously published monazite U–Pb data presented in De Vries Van Leeuwen et al. (2021). From these same monazite grains, 42 Lu–Hf spot analyses were collected, of which 11 were from mnz1 domains and 31 were from mnz2 domains (Fig. 6a). Two analyses from mnz1 and one analysis from mnz2 showed signs of isotopic mixing between the two domains and were not further considered for age calculations (Fig. 6a). Chondrite-normalised REE data from these analyses agree with those attained from U–Pb analyses, with mnz2 analyses exhibiting elevated HREE contents (Fig. 6b). Analyses from mnz1 yield an inverse Lu–Hf isochron age of 601±47 [48] Ma (Fig. 6a; n=9, MSWD = 1.3, p=0.21), whereas analyses from mnz2 yield and inverse Lu–Hf isochron age of 441±11 [13] Ma (Fig. 6a; n=30, MSWD = 2.1, p<0.00).

Thirty-nine U–Pb spot analyses were collected from sample ARK 2017-11, 13 of which belong to mnz1 and 21 to mnz2 (Fig. 5c). Five analyses yield intermediate 206Pb/238U dates between 518 and 478 Ma, with two yielding concordant dates of 515±9 and 518±10 Ma (Fig. 5c). Chondrite-normalised REE data show similar patterns to those in sample ARK 2017-15, with analyses from mnz2 domains exhibiting elevated HREE contents (Fig. 5d). These analyses, as in sample ARK 2017-15, are considered to represent mixing between mnz1 and mnz2. Analyses from mnz1 yield 206Pb/238U dates of 602–569 Ma, whereas those from mnz2 are spread from 460–418 Ma (Fig. 5c). A single analysis from mnz2 yields an anomalously young 206Pb/238U date of 379±13 Ma (Fig. 5c). Similar to sample ARK 2017-15, these data accurately replicate the monazite U–Pb data presented in De Vries Van Leeuwen et al. (2021). From the same grains, 22 Lu–Hf spot analyses were collected from sample ARK 2017-11, of which 17 were from mnz1 domains and 5 were from mnz2 domains (Fig. 6c). Chondrite-normalised REE data from these analyses agree with those attained from U–Pb analyses, with mnz2 analyses exhibiting elevated HREE contents compared to those from mnz1 (Fig. 6d). A single mixed analysis from mnz1 was excluded from age calculations (Fig. 6c). Analyses from mnz1 yielded an inverse Lu–Hf isochron age of 589±37 [38] Ma (Fig. 6c; n=16, MSWD = 1.3, p=0.18), while analyses from mnz2 yielded and inverse Lu–Hf isochron age of 467±62 [62] Ma (Fig. 6c; n=5, MSWD = 0.47, p=0.80).

5 Discussion

5.1 Monazite reference materials

The two monazite reference materials with published ID-TIMS data investigated in this study, RW-1 and TS-Mnz, both yield inverse Lu–Hf isochron and weighted mean dates that lie within 2 SE uncertainty and are accurate to <1 % of their published U–Pb ages (Figs. 1, 2; Budzyń et al., 2021; Ling et al., 2017). This demonstrates that the in situ Lu–Hf approach via LA-ICP-MS/MS, corrected for matrix-dependent fractionation to apatite reference materials, faithfully reproduces the published ID-TIMS/ID-MC-ICP-MS U–Pb ages for monazite reference materials RW-1 and TS-Mnz (Budzyń et al., 2021; Ling et al., 2017).

Across two analytical sessions, RW-1 returned uncorrected inverse Lu–Hf isochron dates of 947±11 Ma (n=26, MSWD = 1.6, p=0.02) and 949±11 Ma (n=30, MSWD = 1.9, p<0.00) and TS-Mnz returned uncorrected inverse Lu–Hf isochron ages of 955±11 Ma (n=26, MSWD = 1.9, p<0.00) and 955±7 Ma (n=26, MSWD = 1.4, p=0.07), corresponding to apparent age offsets from their published U–Pb ID-TIMS ages of ∼4.5 %–5.0 %. If these apparent age offsets are converted to matrix fractionation correction factors, RW-1 yields values of 4.78±0.06 % and 5.00±0.06 %, while TS-Mnz yields values of 4.94±0.06 % and 4.94±0.04 % for sessions 1 and 2, respectively. These values deviate slightly (<1 %) from the matrix fractionation correction factors attained from apatite reference materials. However, the resulting <0.5 % age difference for RW-1 and TS-Mnz when correcting to monazite versus apatite is within uncertainty. The similarity of matrix fractionation correction factors, along with the negligible common Hf contents and relatively high Lu contents (Budzyń et al., 2021; Ling et al., 2017), indicates that these reference monazites would be appropriate for calibrating unknown samples. Hence, although RW-1 and TS-Mnz were used here as secondary reference materials, they can reliably be used to calibrate Lu–Hf ratios for matrix-dependent fractionation in future studies. In their recent study, Wu et al. (2024) also measured Lu–Hf ages for RW-1 but did not present the data, precluding a direct comparison between instruments and laboratories.

Although the Storø and Pilbara monazites are not as well-characterised as RW-1 and TS-Mnz, both yield inverse Lu–Hf isochron and weighted mean dates that fall within ∼1 % of their published U–Pb ages (Figs. 1, 2). This suggests they are suitable as secondary reference materials for evaluating the accuracy of post-acquisition calibrations and corrections (see above). Since Storø originates from a metasedimentary rock and Pilbara from a granitoid, it is evident that monazites from diverse rock types can serve as secondary reference materials, provided they meet the following criteria: (1) sufficient Lu content, (2) negligible common Hf, and (3) consistent Lu–Hf and U–Pb dates. In this regard, laboratories routinely performing in situ U–Pb monazite dating via LA-ICP-MS likely possess various in-house monazite reference materials that could also be used for Lu–Hf dating.

5.2 Comparing in situ U–Pb and Lu–Hf data from complex samples

In situ Lu–Hf geochronological data from samples ARK 2017-15 and ARK 2017-11 produce dates that lie within the spread of 206Pb/238U dates for both the mnz1 and mnz2 domains (ca. 600–570 Ma and ca. 460–390 Ma). This highlights that in situ Lu–Hf isotopic data attained via LA-ICP-MS/MS have the capacity to replicate ages attained via U–Pb LA-ICP-MS geochronology in geologically complex samples. Furthermore, it can resolve multiple age populations from samples which exhibit significant intra-grain complexity, provided careful microstructural targeting is performed and companion trace element data are acquired.

U–Pb data from both mnz1 and mnz2 domains in both ARK samples exhibit large dispersion in concordant U–Pb dates. De Vries Van Leeuwen et al. (2021) argue that this dispersion corresponds to prolonged fluid-mediated dissolution–reprecipitation given the thermally energetic environment in which these rocks were metamorphosed. Although excess analytical scatter and/or common Hf incorporation cannot be ruled out, monazite dissolution–reprecipitation may also explain the overdispersion of the Lu–Hf dates for mnz2 in sample ARK 2017-15 (MSWD = 2.1). This would suggest that (partial) dissolution of monazite effectively expels radiogenic Hf and its re-uptake during reprecipitation is limited. This in turn preserves the timing (and/or time span) of fluid–rock interaction, behaving much the same as Pb during the same process (e.g. Harlov et al., 2011; Seydoux-Guillaume et al., 2002).

https://gchron.copernicus.org/articles/7/199/2025/gchron-7-199-2025-f07

Figure 7Plot depicting the concentration of radiogenic 176Hf accumulated for different Lu concentrations (1, 5, 10, 20, and 50 ppm) as a function of age. The dashed line indicates the approximate lower limit of detection on 176Hf achieved in this study.

Download

5.3 Applications and limitations

In situ Lu–Hf dating of monazite via LA-ICP-MS/MS presents an opportunity to gain geochronological information from a separate decay system to the commonly utilised U–Pb series. As demonstrated in this study, in situ Lu–Hf dating can aid in interpreting the ambiguous dispersion of U–Pb dates, which is a well-documented feature of purportedly long-lived metamorphic systems (Clark et al., 2024; De Vries Van Leeuwen et al., 2021; Kirkland et al., 2016; Korhonen et al., 2013). It is debated whether this dispersion is due to the low degrees of analytically unresolvable Pb loss which generate seemingly protracted spreads of concordant dates or if it represents a truly prolonged episode of metamorphism and/or fluid–rock interaction (e.g. Kirkland et al., 2016). However, when two separate isotopic systems both indicate dispersion (as seen in this study), it bolsters the argument that this represents a true geological signal. In addition to this, the Lu–Hf system in monazite may be particularly useful in situations where the U–Pb system has been compromised or obfuscated by processes such as Pb loss, non-radiogenic Pb incorporation, excess 206Pb due to 230Th uptake, or low U concentrations. Furthermore, the diffusivity of Lu and Hf in monazite is currently poorly understood. However, it is likely that their behaviour differs from that of U, Th, and Pb. As such, the Lu–Hf system may have the ability to resolve geological processes that are not preserved (or are highly obscured) by the U–Pb system.

At present, the greatest limiting factor for this technique is the analytical sensitivity of quadrupole instruments. Given the slow rate of radiogenic 176Hf ingrowth, the accumulation of resolvable quantities is highly dependent on initial Lu concentrations (Fig. 7). Although monazite Lu concentrations are highly variable (dictated by numerous factors such as host rock/protolith composition, PT conditions of metamorphism, and magma/fluid chemistry), they often fall in the range of 1 to 50 ppm. Across the two analytical sessions conducted during this study, 176Hf detection limits were typically observed to be ∼1–3 ppb. From the data presented in Fig. 7, it is evident that Lu-rich monazites (50 ppm) can resolve ages of ∼50 Ma, while Lu-poor monazites (1 ppm) would take more than 2000 Myr to reach 176Hf concentrations at the lower limit of detection. This can be compensated for by increasing laser beam diameters (thus increasing count rates), albeit at the expense of spatial resolution.

6 Conclusions

In situ Lu–Hf dating of monazite via LA-ICP-MS/MS faithfully reproduces published U–Pb ages of two monazite reference materials, RW-1 and TS-Mnz. We further demonstrate the approach for monazite from Arkaroola in South Australia, which formed during a complex and protracted geological history. These data replicate U–Pb geochronological data collected from the same grains and demonstrate that the Lu–Hf system within monazite is sensitive to resetting during fluid-mediated dissolution–reprecipitation. In situ Lu–Hf geochronology may find use in scenarios where the U–Pb system in monazite has been compromised (e.g. Pb loss, common Pb contamination) and is unable to provide reasonable geological information.

Data availability

The data used in this article are available at https://doi.org/10.25909/27441327.v4 (De Vries Van Leeuwen, 2024).

Author contributions

ATDVVL: conceptualisation, investigation, writing (original draft), visualisation; SG: conceptualisation, investigation, methodology, writing (review and editing); MH: conceptualisation, writing (review and editing); JM: resources, writing (review and editing); SEG: methodology, investigation, writing (review and editing).

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

The authors acknowledge the instruments and expertise of Microscopy Australia (ROR: 042mm0k03) at Adelaide Microscopy, University of Adelaide, enabled by NCRIS, university, and state government support. Ben Wade and Ken Neubauer from Adelaide Microscopy are thanked for their help with the collection of SEM data used in this study. Stijn Glorie is supported by an Australian Research Council Future Fellowship (FT210100906). Jacob Mulder is supported by an Australian Research Council Discovery Early Career Researcher Award (DE240101283). Kristoffer Szilas and Nick Gardiner are thanked for providing samples of the Storø Quartzite. Ryan Ickert is thanked for the editorial handling of this paper. Nick Roberts and Stephanie Walker are thanked for providing thoughtful and constructive reviews. This study was supported by the Mineral Exploration Cooperative Research Centre whose activities are funded by the Australian Government's Cooperative Research Centre Programme. This is MinEx CRC Document 2025/09.

Financial support

This research has been supported by the Australian Research Council (grant nos. FT210100906 and DE240101283) and the Mineral Exploration Cooperative Research Centre (Project 7).

Review statement

This paper was edited by Ryan Ickert and reviewed by Nick Roberts and Stephanie Walker.

References

Budzyń, B., Sláma, J., Corfu, F., Crowley, J., Schmitz, M., Williams, M. L., Jercinovic, M. J., Kozub-Budzyń, G. A., Konečny, P., Rzepa, G., and Włodek, A.: TS-Mnz – A new monazite age reference material for U-Th-Pb microanalysis, Chem. Geol., 572, 120195, https://doi.org/10.1016/j.chemgeo.2021.120195, 2021. 

Clark, C., Brown, M., Knight, B., Johnson, T. E., Mitchell, R. J., and Gupta, S.: Ultraslow cooling of an ultrahot orogen, Geology, 52, 880–884, https://doi.org/10.1130/G52442.1, 2024. 

De Biévre, P. and Taylor, P. D. P.: Table of the isotopic compositions of the elements, Int. J. Mass Spectrom., 123, 149–166, https://doi.org/10.1016/0168-1176(93)87009-H, 1993. 

De Vries Van Leeuwen, A.: Monazite Lu–Hf and U–Pb data – A comparison between in situ monazite Lu–Hf and U–Pb geochronology, The University of Adelaide [data set], https://doi.org/10.25909/27441327.v4, 2024. 

De Vries Van Leeuwen, A. T., Hand, M., Morrissey, L. J., and Raimondo, T.: Th–U powered metamorphism: Thermal consequences of a chemical hotspot, J. Metamorph. Geol., 39, 541–565, https://doi.org/10.1111/jmg.12590, 2021. 

Elburg, M. A., Andersen, T., Bons, P. D., Simonsen, S. L., and Weisheit, A.: New constraints on Phanerozoic magmatic and hydrothermal events in the Mt Painter Province, South Australia, Gondwana Res., 24, 700–712, https://doi.org/10.1016/j.gr.2012.12.017, 2013. 

Erickson, T. M., Pearce, M. A., Taylor, R. J. M., Timms, N. E., Clark, C., Reddy, S. M., and Buick, I. S.: Deformed monazite yields high-temperature tectonic ages, Geology, 43, 383–386, https://doi.org/10.1130/G36533.1, 2015. 

Foden, J., Elburg, M. A., Dougherty-Page, J., and Burtt, A.: The Timing and Duration of the Delamerian Orogeny: Correlation with the Ross Orogen and Implications for Gondwana Assembly, J. Geol., 114, 189–210, https://doi.org/10.1086/499570, 2006. 

Gardiner, N. J., Mulder, J. A., Szilas, K., Nebel, O., Whitehouse, M., Jeon, H., and Cawood, P. A.: A record of Neoarchaean cratonisation from the Storø Supracrustal Belt, West Greenland, Earth Planet. Sc. Lett., 602, 117922, https://doi.org/10.1016/j.epsl.2022.117922, 2023. 

Glorie, S., Mulder, J., Hand, M., Fabris, A., Simpson, A., and Gilbert, S.: Laser ablation (in situ) Lu-Hf dating of magmatic fluorite and hydrothermal fluorite-bearing veins, Geosci. Front., 14, 101629, https://doi.org/10.1016/j.gsf.2023.101629, 2023. 

Glorie, S., Hand, M., Mulder, J., Simpson, A., Emo, R. B., Kamber, B., Fernie, N., Nixon, A., and Gilbert, S.: Robust laser ablation Lu–Hf dating of apatite: an empirical evaluation, Geological Society, London, Special Publications, 537, 165–184, https://doi.org/10.1144/SP537-2022-205, 2024a. 

Glorie, S., Simpson, A., Gilbert, S. E., Hand, M., and Müller, A. B.: Testing the reproducibility of in situ LuHf dating using Lu-rich garnet from the Tørdal pegmatites, southern Norway, Chem. Geol., 653, 122038, https://doi.org/10.1016/j.chemgeo.2024.122038, 2024b. 

Harlov, D. E., Wirth, R., and Hetherington, C. J.: Fluid-mediated partial alteration in monazite: the role of coupled dissolution–reprecipitation in element redistribution and mass transfer, Contrib. Mineral. Petr., 162, 329–348, https://doi.org/10.1007/s00410-010-0599-7, 2011. 

Kelly, N. M., Harley, S. L., and Möller, A.: Complexity in the behavior and recrystallization of monazite during high-T metamorphism and fluid infiltration, Chem. Geol., 322–323, 192–208, https://doi.org/10.1016/j.chemgeo.2012.07.001, 2012. 

Kirkland, C. L., Erickson, T. M., Johnson, T. E., Danišík, M., Evans, N. J., Bourdet, J., and McDonald, B. J.: Discriminating prolonged, episodic or disturbed monazite age spectra: An example from the Kalak Nappe Complex, Arctic Norway, Chem. Geol., 424, 96–110, https://doi.org/10.1016/j.chemgeo.2016.01.009, 2016. 

Kohn, M. J. and Malloy, M. A.: Formation of monazite via prograde metamorphic reactions among common silicates: implications for age determinations, Geochim. Cosmochim. Ac., 68, 101–113, https://doi.org/10.1016/S0016-7037(03)00258-8, 2004. 

Korhonen, F. J., Clark, C., Brown, M., Bhattacharya, S., and Taylor, R. J. M.: How long-lived is ultrahigh temperature (UHT) metamorphism? Constraints from zircon and monazite geochronology in the Eastern Ghats orogenic belt, India, Precambrian Res., 234, 322–350, https://doi.org/10.1016/j.precamres.2012.12.001, 2013. 

Larson, K. P., Shrestha, S., Cottle, J. M., Guilmette, C., Johnson, T. A., Gibson, H. D., and Gervais, F.: Re-evaluating monazite as a record of metamorphic reactions, Geosci. Front., 13, 101340, https://doi.org/10.1016/j.gsf.2021.101340, 2022. 

Ling, X.-X., Huyskens, M. H., Li, Q.-L., Yin, Q.-Z., Werner, R., Liu, Y., Tang, G.-Q., Yang, Y.-N., and Li, X.-H.: Monazite RW-1: a homogenous natural reference material for SIMS U–Pb and Th–Pb isotopic analysis, Miner. Petrol., 111, 163–172, https://doi.org/10.1007/s00710-016-0478-7, 2017. 

Lloyd, J. C., Blades, M. L., Counts, J. W., Collins, A. S., Amos, K. J., Wade, B. P., Hall, J. W., Hore, S., Ball, A. L., Shahin, S., and Drabsch, M.: Neoproterozoic geochronology and provenance of the Adelaide Superbasin, Precambrian Res., 350, 105849, https://doi.org/10.1016/j.precamres.2020.105849, 2020. 

Maidment, D. W.: Palaeozoic high-grade metamorphism within the Centralian Superbasin, Harts Range region, central Australia, Thesis, https://doi.org/10.25911/5d763319d0e88, 2005. 

McLaren, S., Sandiford, M., Powell, R., Neumann, N., and Woodhead, J. O. N.: Palaeozoic Intraplate Crustal Anatexis in the Mount Painter Province, South Australia: Timing, Thermal Budgets and the Role of Crustal Heat Production, J. Petrol., 47, 2281–2302, https://doi.org/10.1093/petrology/egl044, 2006. 

Mildren, S. D. and Sandiford, M.: Heat refraction and low-pressure metamorphism in the northern Flinders Ranges, South Australia, Aust. J. Earth Sci., 42, 241–247, https://doi.org/10.1080/08120099508728198, 1995. 

Nebel, O., Morel, M. L. A., and Vroon, P. Z.: Isotope Dilution Determinations of Lu, Hf, Zr, Ta and W, and Hf Isotope Compositions of NIST SRM 610 and 612 Glass Wafers, Geostand. Geoanal. Res., 33, 487–499, https://doi.org/10.1111/j.1751-908X.2009.00032.x, 2009. 

Norris, A. and Danyushevsky, L.: Towards estimating the complete uncertainty budget of quantified results measured by LA-ICP-MS, in: Goldschmidt Abstracts, Goldschmidt, Boston, MA, USA, 1894, https://goldschmidtabstracts.info/abstracts/abstractView?id=2018003403 (last access: 31 October 2024), 2018. 

Parrish, R. R.: U–Pb dating of monazite and its application to geological problems, Can. J. Earth Sci., 27, 1431–1450, https://doi.org/10.1139/e90-152, 1990. 

Paton, C., Hellstrom, J., Paul, B., Woodhead, J., and Hergt, J.: Iolite: Freeware for the visualisation and processing of mass spectrometric data, J. Anal. Atom. Spectrom., 26, 2508–2518, https://doi.org/10.1039/c1ja10172b, 2011. 

Paul, E., Flöttmann, T., and Sandiford, M.: Structural geometry and controls on basement-involved deformation in the northern Flinders Ranges, Adelaide Fold Belt, South Australia, Aust. J. Earth Sci., 46, 343–354, https://doi.org/10.1046/j.1440-0952.1999.00711.x, 1999. 

Payne, J. L., Hand, M., Barovich, K. M., and Wade, B. P.: Temporal constraints on the timing of high-grade metamorphism in the northern Gawler Craton: Implications for assembly of the Australian Proterozoic, Aust. J. Earth Sci., 55, 623–640, https://doi.org/10.1080/08120090801982595, 2008. 

Pearce, N. J. G., Perkins, W. T., Westgate, J. A., Gorton, M. P., Jackson, S. E., Neal, C. R., and Chenery, S. P.: A Compilation of New and Published Major and Trace Element Data for NIST SRM 610 and NIST SRM 612 Glass Reference Materials, Geostandards Newsletter, 21, 115–144, https://doi.org/10.1111/j.1751-908X.1997.tb00538.x, 1997. 

Preiss, W. V.: The Adelaide Geosyncline of South Australia and its significance in Neoproterozoic continental reconstruction, Precambrian Res., 100, 21–63, https://doi.org/10.1016/S0301-9268(99)00068-6, 2000. 

Rubatto, D., Williams, I. S., and Buick, I. S.: Zircon and monazite response to prograde metamorphism in the Reynolds Range, central Australia, Contrib. Mineral. Petr., 140, 458–468, https://doi.org/10.1007/PL00007673, 2001. 

Rubatto, D., Chakraborty, S., and Dasgupta, S.: Timescales of crustal melting in the Higher Himalayan Crystallines (Sikkim, Eastern Himalaya) inferred from trace element-constrained monazite and zircon chronology, Contrib. Mineral. Petr., 165, 349–372, https://doi.org/10.1007/s00410-012-0812-y, 2013. 

Seydoux-Guillaume, A.-M., Paquette, J.-L., Wiedenbeck, M., Montel, J.-M., and Heinrich, W.: Experimental resetting of the U–Th–Pb systems in monazite, Chem. Geol., 191, 165–181, https://doi.org/10.1016/S0009-2541(02)00155-9, 2002. 

Simpson, A., Gilbert, S., Tamblyn, R., Hand, M., Spandler, C., Gillespie, J., Nixon, A., and Glorie, S.: In-situ LuHf geochronology of garnet, apatite and xenotime by LA ICP MS/MS, Chem. Geol., 577, 120299, https://doi.org/10.1016/j.chemgeo.2021.120299, 2021. 

Simpson, A., Glorie, S., Hand, M., Spandler, C., Gilbert, S., and Cave, B.: In situ Lu–Hf geochronology of calcite, Geochronology, 4, 353–372, https://doi.org/10.5194/gchron-4-353-2022, 2022. 

Simpson, A., Glorie, S., Hand, M., Spandler, C., and Gilbert, S.: Garnet Lu-Hf speed dating: A novel method to rapidly resolve polymetamorphic histories, Gondwana Res., 121, 215–234, https://doi.org/10.1016/j.gr.2023.04.011, 2023. 

Simpson, A., Glorie, S., Hand, M., Gilbert, S. E., Spandler, C., Dmitrijeva, M., Swain, G., Nixon, A., Mulder, J., and Münker, C.: In situ apatite and carbonate Lu-Hf and molybdenite Re-Os geochronology for ore deposit research: Method validation and example application to Cu-Au mineralisation, Geosci. Front., 15, 101867, https://doi.org/10.1016/j.gsf.2024.101867, 2024. 

Söderlund, U., Patchett, P. J., Vervoort, J. D., and Isachsen, C. E.: The 176Lu decay constant determined by Lu–Hf and U–Pb isotope systematics of Precambrian mafic intrusions, Earth Planet. Sc. Lett., 219, 311–324, https://doi.org/10.1016/S0012-821X(04)00012-3, 2004. 

Spencer, C. J., Kirkland, C. L., Roberts, N. M. W., Evans, N. J., and Liebmann, J.: Strategies towards robust interpretations of in situ zircon Lu–Hf isotope analyses, Geosci. Front., 11, 843–853, https://doi.org/10.1016/j.gsf.2019.09.004, 2020.  

Thompson, J., Meffre, S., Maas, R., Kamenetsky, V., Kamenetsky, M., Goemann, K., Ehrig, K., and Danyushevsky, L.: Matrix effects in Pb/U measurements during LA-ICP-MS analysis of the mineral apatite, J. Anal. Atom. Spectrom., 31, 1206–1215, https://doi.org/10.1039/C6JA00048G, 2016. 

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. 

Vermeesch, P.: Errorchrons and anchored isochrons in IsoplotR, Geochronology, 6, 397–407, https://doi.org/10.5194/gchron-6-397-2024, 2024. 

Wu, S., Niu, J., Yang, Y., Wang, H., Yang, J., and Wu, F.: First demonstration of in situ Lu–Hf dating using LA-ICP-MS/MS applied to monazite, J. Anal. At. Spectrom., 39, 2703–2715, https://doi.org/10.1039/D4JA00258J, 2024. 

Yu, J., Glorie, S., Hand, M., Simpson, A., Gilbert, S., Szilas, K., Roberts, N., Pawley, M., and Cheng, Y.: Laser ablation (in situ) Lu-Hf geochronology of epidote group minerals, Contrib. Mineral. Petr., 179, 62, https://doi.org/10.1007/s00410-024-02143-y, 2024. 

Download
Short summary
In this contribution, we demonstrate in situ monazite lutetium–hafnium dating and compare results with uranium–lead dating. We present data from monazite reference materials and complex samples to demonstrate the viability of this method. We show that in situ lutetium–hafnium dating of monazite can resolve multiple age populations and may find use in scenarios where the uranium–lead system has been compromised.
Share