ESR-thermochronometry of the Hida range of the Japanese Alps: Validation and future potential

The electron spin resonance (ESR) of quartz has previously been shown to have potential for determining rock cooling histories, however this technique remains underdeveloped. In this study, we explore the ESR of a suite of samples from the Hida range of the Japanese Alps. We develop measurement protocols and models to constrain the natural trapped 15 charge concentration as well as the parameters that govern signal growth and signal thermal decay. The thermal stability of the Al and Ti-centres is similar to that of the luminescence of feldspar. Inverting the ESR data for cooling yields similar thermal histories to paired luminescence data from the same samples. However, a series of synthetic inversions show that whereas the luminescence of feldspar can only resolve minimum cooling histories of ~160 °C/Myr over timescales of 103-5 a, quartz ESR may resolve cooling histories as low as 25-50 °C/Myr over timescales of 103-7 a. This difference arises because quartz ESR 20 has a higher dating limit than the luminescence of feldspar. These results imply that quartz ESR will be widely applicable in the constraint of late-stage rock cooling histories, providing new insights into landscape evolution over late-Quaternary timescales.


Introduction
Thermochronometry based on trapped-charge dating allows the constraint of late stage exhumation and/or rock thermal 25 histories at the scale of glacial-interglacial cycles (e.g. Biswas et al., 2018). Following the study of Herman et al. (2010) which applied optically stimulated luminescence (OSL) dating to constrain the exhumation histories of the Southern Alps of New Zealand, there have been a number of both methodological and applied studies that have almost exclusively focused on luminescence dating (see King et al., 2016a andKing, 2018 for reviews). In this study we explore the potential of a second trapped-charge dating method, electron spin resonance (ESR) of quartz, for ultra-low temperature (i.e. < 100°C) 30 thermochronometry.
Electron spin resonance can be used to measure the time-dependent accumulation of unpaired electrons (paramagnetic centres) in minerals such as quartz (cf. Grün, 1989;Ikeya, 1993). As for luminescence dating, when a mineral is exposed to ionizing radiation, electrons are excited from their ground state in the valence band, to the conduction band. Almost immediately most 35 electrons fall back to the valence band, recombining with the "holes" of positive charge created by the electron's excitation.
However, some electrons become trapped within defects in the crystal lattice, caused by element vacancies or substitutions. In this study we specifically target the Al (hole trapping) centre and the Ti (electron trapping) centre, although other defects such as the E' (oxygen vacancy) centre could also be investigated (e.g. Grün et al., 1999). The Al-centre comprises a hole located at AlO - (Nuttall and Weil, 1981) whereas the Ti-centre comprises the substitution of Si 4+ with Ti 3+ + estabilized with H + , Li + or Na + (Rinneberg and Weil, 1972;Isoya et al., 1983). ESR offers a key advantage over luminescence dating, specifically that 5 ESR signals saturate at higher doses (Rink, 1997;Tsukamoto et al., 2018). Within the context of thermochronometry, this means that whilst the application of luminescence thermochronometry remains geographically limited to regions experiencing extremely rapid cooling/exhumation higher than several mm/yr e.g. New Zealand (Herman et al., 2010), eastern Himalayan syntaxis (King et al., 2016b), ESR thermochronometry could be much more widely applied. 10 The potential of ESR for thermochronometry has been recognized previously. Following from an earlier study (Ikeya, 1983), Toyoda and Ikeya (1991) first suggested that the intensity of quartz ESR centres could be used to determine the lowtemperature thermal histories of the host rock. Scherer et al. (1993; investigated changing ESR centre intensities with depth through the known-thermal history KTB borehole in Germany (Coyle et al., 1997), which has also been used to validate the luminescence thermochronometry technique (Guralnik et al., 2015, Biswas et al., 2018. Scherer et al. (1993;1994) 15 recorded a qualitative reduction in signal intensity of the Al-centre with increasing temperature and depth. In contrast, data for the Ti-centre were much more scattered with zero signal intensity recorded for many samples. However, it was Grün et al. (1999) who reported the first quantitative ESR-thermochronometry results from their study of the Eldzhurtinskiy Granite from the Russian Caucasus. Using the Al and Ti-centres of quartz, they obtained cooling rates of between 160 °C/Myr and 600 °C/Myr which correspond to denudation rates of ~2. 5 and 5.5 mm/yr. 20 Despite the potential illustrated by ESR-thermochronometry in these early studies, the technique has not been applied since, in part associated with the difficulties of making ESR measurements (i.e. gamma or X-ray source availability, absence of automated instrumentation). In this study, we investigate the potential of ESR-thermochronometry through applying new measurement protocols (Tsukamoto et al., 2015), which have been facilitated by developments in instrumentation (Oppermann 25 and Tsukamoto, 2015), and that have recently been validated against samples with independent age control (Richter et al., In Press). We propose a kinetic model inspired by recent progress in luminescence thermochronometry (Lambert, 2018) to facilitate the inference of rock thermal histories from ESR laboratory data and perform a series of synthetic inversions to evaluate the range of cooling histories that ESR thermochronometry may be applicable over. We then investigate six rock samples from the Japanese Alps and contrast their ESR thermal histories with those obtained from optically stimulated 30 luminescence thermochronometry of feldspar (e.g. Guralnik et al., 2015;King et al., 2016b,c) of the same samples.

Theoretical basis
The theoretical basis of ESR-thermochronometry is very similar to that of luminescence thermochronometry (cf. King et al., 2016a;Herman and King, 2018 for reviews), with the advantage that unlike feldspar minerals, quartz minerals are not thought to suffer from athermal signal losses. Here we present the kinetic model for ESR-thermochronometry, before discussing how the parameters that describe signal growth and signal thermal decay can be constrained in the laboratory. 5

Kinetic model
We propose the following kinetic models to describe the evolution of ESR signals with temperature. A saturating system may be described by: : ; < [ $( 4 , )] (1) 10 and a non-saturating system can be described by: : ; < [ $( 4 , )] (2) 15 where, and, 20 where $ is the trapped charge population with activation energy, Ea (eV). In the instance of a saturating system $ is expressed as a saturation ratio, but for a non-saturating system it is expressed as absorbed radiation dose (Gy). The first term on the right-25 hand side of Eqs. (1) and (2) describes charge trapping as a first-order process. For a non-saturating system, / is defined by the environmental dose rate ̇ (Gy), whereas for a saturating system, / is defined as ̇/ B where B is the characteristic dose of saturation (Gy). The second term on the right-hand side of Eqs. (1) and (2) describes thermal charge detrapping, and here we benefit from recent advances made in luminescence thermochronometry, and follow Lambert (2018) by describing thermal detrapping using a model that assumes a Gaussian distribution of activation energies, Ea around the mean trap depth, 30 ( ) ) (eV). Thermal detrapping is also described by the frequency factor, s (s -1 ), the Boltzman constant, kB (eV), temperature, T (K) and P(Ea) the probability of thermally evicting electrons (or holes) from the trap (Eq. (4)). An alternative approach could be to use a first or second order kinetic model as has been done previously (Toyoda and Ikeya, 1991;Ikeya, 1993;Grün et al., 1999) and we discuss our model selection more completely in the supplementary material. 5

Constraining charge trapping
The natural trapped charge concentration, which reflects the equilibrium between charge trapping and thermally stimulated charge detrapping, can be measured in the laboratory through the development of a sample specific radiation dose response curve. This comprises measurement of a sample following increasingly large laboratory radiation doses, and interpolation of the natural ESR signal onto the resultant dose response curve. Measurements can either be made on single (e.g. Tsukamoto et 10 al., 2015) or using multiple aliquots of the same sample (e.g. Grün et al., 1999). The former approach has only recently been made practical, following the introduction of X-ray irradiation for regenerative dosing (Oppermann and Tsukamoto, 2015), as opposed to gamma irradiation which is often done at a laboratory separate to the measurement laboratory.

Constraining charge detrapping
Thermal detrapping can be measured following laboratory isothermal decay experiments, whereby aliquots of a sample are 15 given a radiation dose before being heated at different temperatures for different durations. The resultant signal loss is measured and fitted with the kinetic model described in Eqs. (1-4). Previous investigations have suggested that the thermal decay of quartz ESR can be described by first order or second order kinetics. Here, instead we use a density of states model, originally developed for the luminescence of feldspar (Li and Li, 2013;Lambert, 2018; further details of model selection are given in the Supplementary Material). The selected model is based on a Gaussian distribution of activation energies ( ) ), around the 20 mean trap-depth, ( ) ) (Lambert, 2018), and may be applicable for quartz ESR data where electrons can be trapped in a variety of different defects, e.g. Ti 3+ + echarge compensated by H + , Li + or Na + (Tsukamoto et al., 2018).

Assessing the potential of ESR-thermochronometry
Electron spin resonance dating analyses are not automated, meaning that the laboratory measurements required for ESRthermochronometry analyses are considerably more time-consuming than those required for luminescence 25 thermochronometry. It is thus necessary to verify that ESR-thermochronometry offers advantages over luminescence methods.
To achieve this, a series of synthetic tests for known cooling histories were done using the kinetic parameters of sample KRG16-06 (Table 1). These tests first comprised running a forward model, which uses sample-specific kinetic parameters and a rate equation to describe signal growth. Through forward modelling, it is possible to predict the trapped-charge concentration for a particular cooling history. The second stage of the test comprised inverting the trapped-charge concentrations predicted 30 by the forward model, using the same rate equation, to determine if it is possible to recover the cooling history used in the forward model prediction. Further details of the forward and inverse modelling are given below.

Forward modelling
Five different monotonic cooling scenarios were used to test the potential of ESR-thermochronometry in comparison to optically stimulated luminescence (OSL) thermochronometry, comprising cooling with rates of 100 °C/Myr, 75 °C/Myr, 50 5 °C/Myr, 25 °C/Myr and no cooling (i.e. isothermal holding at 0 °C for 2 Myr). All cooling rates were maintained for at least 2 Myr with a starting temperature of 200 °C which is greater than the anticipated closure temperature of the ESR system (cf. Grün et al., 1999;Scherer, 1993;. Using the kinetic model in Eqs. (1) and (2), a trapped charge population, $ TUV was predicted for both the Ti and Al-centres respectively using the kinetic parameters of sample KRG16-06 (Table 1) for the five different scenarios. In addition, the same exercise was carried out for four feldspar multi-OSL-thermochronometry signals of 10 the same sample using the following kinetic model, after King et al., 2016a (see supplementary information for further details on model selection): 15 where the total accumulation of charge with time, i.e. $( ) is obtained by integrating $( _ ,`, ) over the range of band-tail states, `, and an infinite range of dimensionless distances, _ : where (`) is the probability of evicting electrons into band-tail states of energy `+`, defined as: where B is a pre-exponential multiplier, and where ( _ ) is the probability density distribution of the nearest recombination 25 centre defined by Huntley (2006) as: where dimensionless distance _ ≡ p qHe r s f g , the dimensionless density of recombination centres _ ≡ qHe ru g and is a constant related to the Bohr radius of the electron trap (Huntley, 2006;Kars et al., 2008;Tachiya and Mozumder, 1974). Table 1 3.2 Inverse modelling 5 We inverted the five sets of $ TUV values for the ESR and OSL data described above using a similar approach to King et al. (2016a), which we briefly outline here. The trapped-charge populations were modelled for 10,000 randomly generated timetemperature histories (t-T paths), which were constrained to cool monotonically between 200 °C and 0±5 °C, over 2 Ma. We computed the dose response curves by solving the differential equations described above using a semi implicit Euler method (Press, 2007). For each t-T path we calculated a misfit between the final inverted trapped-charge population, $ wxV , and our 10 forward modelled values, $ TUV (Wheelock et al., 2015), from which the misfit, M and likelihood, L are calculated: for m traps, where σ is the uncertainty. An arbitrary uncertainty on $ TUV of 10% was assumed. Cooling histories are then accepted or rejected by contrasting L with a random number between 0 and 1; if L is greater, the cooling history is retained.
The accepted cooling histories are finally combined to construct a time-temperature history probability density function through dividing the time-temperature axis into 50 intervals and summing the number of paths that cross through each of the 20 different cells. The Al, Ti and OSL data were first inverted separately and then the Al and Ti-centres were inverted together.
The results of the forward modelling and the synthetic inversions for the ESR and OSL data are shown in Fig. 1. The OSL signals for all cooling histories reach saturation (Fig. 1c), and this is reflected in the failure of the OSL to recover any of the cooling histories when inverted. This is apparent because the 1s confidence intervals show a broad range, with the highest 25 density of cooling histories concentrated at temperatures < 20 °C over the past 500 ka indicating that the luminescence signals are saturated (as shown in Fig. 1c). The minimum cooling rate that can be resolved using OSL for sample KRG16-06 is ~160 °C/Myr, calculated from 86% of the luminescence signal saturation level. Signal saturation is the key limitation that restricts the application of luminescence thermochronometry to regions undergoing rapid exhumation. In contrast, it is clear that the ESR data are able to resolve the 100 °C/Myr, 75 °C/Myr and 50 °C/Myr synthetic cooling histories, and cooling rates of 25 30 °C/Myr are distinct from isothermal holding at 0 °C over timescales of ~2 Ma. This is apparent because of the coincidence between the prescribed cooling histories (white lines) and the highest density of accepted cooling histories shown by the brightest colours in the probability density functions. These results are significant as they show that ESR-thermochronometry is applicable in a range of geological settings beyond the rapidly exhuming locations that luminescence-thermochronometry is currently restricted to.

Proof of concept -Hida range, Japanese Alps
To further explore the potential of the ESR method we applied it to a suite of samples from the Hida range of the Japanese Alps. The Japanese Alps which reach elevations of up to 3,000 m are thought to have uplifted since the Pliocene or Quaternary (Yonekura et al., 2001;Tokahashi, 2006) in response to E-W compressional tectonic forces (Takahashi, 2006;Townend and 10 Zoback, 2006;Sueoka et al., 2016). Lithology of the Hida range is dominated by granitic intrusions, including the Kurobegawa granite, which is the youngest known intrusion on Earth and which was emplaced between 10-0.8 Ma ago (Ito et al., 2013;2017). Previous efforts to apply apatite fission-track dating on the Kurobegawa granite have been unsuccessful because of the very low fission-track density (Yamada, 1999). Extremely young apatite (0.50 ± 0.04 Ma) and zircon helium ages (0.37 ± 0.10 Ma) have recently been reported (Spencer et al., 2019), indicating that exhumation in this region has remained rapid throughout 15 the Quaternary Period.
Six bedrock samples were taken from the Kurobegawa granite, northern Hida range of the Japanese Alps. Four surface samples were taken and form an elevation transect, whilst a further two samples were taken from a high-temperature tunnel, which has a present-day temperature of ~40-50 °C but which had temperatures of up to 165 °C at the time of excavation in the late 1930s 20 (Yuhara and Yamamoto, 1983). Samples had a minimum size of 15 x 15 x 15 cm, to ensure that a light safe portion could be extracted from their interiors. Sample details are given in the Supplementary Material.

Sample preparation
Bedrock samples were prepared using standard laboratory methods under subdued red light conditions at the University of Lausanne and University of Bern, Switzerland (cf. King et al., 2016c). At least 10 mm was cut from the exterior of the samples 25 using a water-cooled diamond saw, to extract the light safe interior. A thin section was made using a representative sample of the bedrock exterior and a further representative sample was sent to ActLabs, Canada for ICP-MS analysis. Sample interiors were then hand crushed to extract the 180-212 µm grain size fraction, which was treated with HCl and H2O2 to remove any carbonates and organic material respectively. The K-feldspar and quartz fractions were separated from heavy minerals using heavy liquids. The K-feldspars were retained for luminescence dating, whilst the quartz extracts (2.58>ρ<2.70 g cm -3 ) were 30 etched for 40 minutes using 40% HF, before being treated with HCl to remove fluorides that had precipitated during etching.
The etched samples were sieved to >150 µm, to remove any partially dissolved feldspar grains. Aliquots for ESR measurement comprised 60 mg of quartz loaded into glass tubes with interior and exterior diameters of 2 and 3 mm respectively.

Environmental dose rate determination
The grain size distribution of quartz and feldspar minerals within the parent bedrock was estimated from thin section analysis using the software of Buscombe (2013). The environmental dose rate, , was calculated from the sample specific radioisotope 5 concentrations using DRAC v.1.2 (Durcan et al., 2015), the conversion factors of Guérin et al. (2011), the alpha grain size attenuation factors of Bell (1980) and the beta grain size attenuation factors of Guérin et al. (2012). Because the bedrock samples have only been at the surface for a short period of time, no cosmic dose rate was included in the calculation. The water content was estimated at 2 ± 2%. For the quartz extract, an etch depth of 10 µm was assumed and the alpha dose rate adjusted following Bell (1980); an a-value of 0.040 ± 0.005 was used after Rees-Jones (1995) for any residual alpha dose. No internal 10 dose rate was included. In contrast, the feldspar fraction was not etched, and an a-value of 0.15 ± 0.05 was used after Balescu and Lamothe (1994). An internal K-content of 12.5 ± 5.0 % was assumed following Huntley and Baril (1997). The calculated environmental dose rates are summarized in Table 1 and full calculation details are given in the Supplementary Material.

Electron Spin Resonance
Electron Spin Resonance measurements were done at the Leibniz Institute for Applied Geophysics in Hannover, Germany. 15 Measurements were made on a JEOL JES-FA100 spectrometer using 2.0 mW microwave power, 0.1 mT modulation width, a 333.5 ± 15 mT magnetic field, 0.1 s time constant and 60 s scan which was averaged over 3 scans. All spectra were measured 3 times following sample turning by 60° to avoid any anisotropic effects. Measurements were made at -150 °C. The instrumentation detailed in Oppermann and Tsukamoto (2015) was used to facilitate X-ray irradiation and sample preheating, which is described below. The Ti and Al-centre peaks were fitted using V3.3.35 of the JEOL ESR data processing software, 20 and were normalized relative to the intensity of the 6 th hyperfine line of Mn 2+ from the internal MgO standard, doped with MnO. As our measurements were carried out at -150 °C, it was not possible to differentiate between the Ti-H and Ti-Li centres, and consequently they have been treated as a single centre (Fig. 2). All subsequent data fitting was done using MATLAB.   Tsukamoto et al. (2015;2018) recently showed that it is necessary to preheat ESR samples that are measured in a single aliquot protocol to avoid any signal contribution from trapped charge that is unstable over laboratory timescales, similar to luminescence dating (cf. Murray and Wintle, 2000). Within this study, a series of tests were done to select the most appropriate preheat temperature and duration. The signal intensities of five aliquots of samples KRG16-06 and KRG16-104 were measured 30 following different preheat treatments (i.e. one aliquot per temperature; Fig. 3a). Aliquots of KRG16-06 were preheated for two minutes at temperatures of between 160 °C and 240 °C, whereas aliquots of sample KRG16-104 were preheated for four minutes at temperatures of between 120 °C and 200 °C. The signal intensity of a further aliquot of each sample was measured without laboratory preheating. In addition to measuring the ESR signal intensity, the equivalent doses of the Ti and Al-centres of KRG16-06 for each preheat temperature, were measured in a single aliquot method (Tsukamoto et al., 2015;Fig. 3b). The single aliquot protocol comprised measurement of the natural signal, measurement of a single additive dose, annealing at 420 5 °C for two minutes and measurement following zero dose. All irradiations were given using an X-ray source with a dose rate of ~0.3 Gy s -1 (Tsukamoto et al., 2018); aliquots were manually turned once during irradiation to ensure that even dosing was achieved.

Measurement of the trapped charge concentration
The trapped charge population of the different samples was measured using a single aliquot approach. This comprised measurement of the natural signal, a zero-point measurement following annealing of the aliquot at 380 °C for four minutes, and measurement of two or three regenerative doses points. The natural signal was then interpolated onto the dose response curve to determine the equivalent dose; all equivalent dose values were calculated using a linear fit. To confirm that the 15 measurement protocol was appropriate, a dose-recovery experiment was done. Three aliquots of zero-age sample KRG16-112 were given an X-ray dose of 360 Gy, before measurement using the same protocol outlined above. Trapped-charge dating systems usually experience signal saturation, therefore it is also necessary to constrain the form of ESR centre dose response.
Using a new aliquot of each sample, dose response was measured using the same measurement protocol, but omitting the zeropoint measurement step, i.e. in an additive dose response protocol. 20

Measurement of trapped charge thermal decay
Thermal signal losses were measured using an isothermal decay experiment, whereby three aliquots of each sample were irradiated with an additive dose of 4.30 kGy. The aliquots were then preheated at 160 °C for four minutes prior to initial measurement, and were then measured following isothermal holding at between 130 °C and 180 °C for 4, 8,16,32,64,128 and 256 minutes. This experiment was also repeated on three fresh aliquots of sample KRG16-104 using a smaller dose of 25 2.15 kGy.

OSL measurements
OSL measurements of all samples followed the approach of King et al. (2016b,c). Luminescence measurements were made at the University of Bern using a single aliquot regenerative dose multiple-elevated-temperature (MET) infra-red stimulated luminescence (IRSL) measurement protocol (Li and Li, 2011) comprising a preheat at 250 °C for 60 seconds, followed by four 30 IRSL measurements at 50, 100, 150 and 225 °C each of 100 s duration. A test dose of 160 Gy was used, which is c. 30% of the IRSL50 signal equivalent dose value of samples KRG16-06, KRG16-101 and KRG16-104. Each measurement cycle was followed by a high temperature optical wash at 290 °C for 60 s. Regenerative doses up to ~4.50 kGy were given to three small (2 mm diameter) aliquots of each sample using two different Risø TL-DA-20 luminescence readers with dose rates ranging from 0.06 to 0.10 Gy s -1 dependent on instrument (dose rates are provided for each measurement in the Supplementary Materials). Luminescence signals were detected in the blue part of the visible spectrum using a BG39 and BG3 or Corning 7-5 59 filter combination. The suitability of the selected measurement protocol was confirmed using a dose recovery test.
Rates of athermal and thermal charge detrapping were also measured using a single aliquot regenerative dose method on the same aliquots used to measure the luminescence dose response curve. Athermal detrapping rates were quantified at room temperature by measuring the luminescence response to a fixed dose following different delay periods. Aliquots were 10 preheated prior to storage following Auclair et al. (2003) and maximum fading delays were 122 days. Rates of thermal charge detrapping were measured using isothermal holding experiments. The aliquots were given a dose of 50 Gy, and held at temperatures ranging from 170 to 350 °C for delay times of 0 to 10,240 s prior to measurement.

Electron Spin Resonance 15
The signal intensity experiment indicates a plateau for the Ti-centre of sample KRG16-104 up until 160 °C (Fig. 3a). In contrast, the Al-centre for this sample and the Ti-centre of sample KRG16-06 reduce in intensity by ~5% between room temperature and 160 °C whilst the Al-centre of sample KRG16-06 is depleted further, by ~10%. The signal intensity data for KRG16-06 are relatively noisy (non-monotonic signal decay with increasing preheat temperature) in comparison to KRG16-104. However, in spite of this, within the preheat-plateau experiment (Fig. 3b), a plateau in De values between 160 and 220 °C 20 is recorded for this sample following preheating for two minutes. On the basis of these experiments a preheat temperature of 160 °C was selected as this temperature maximises signal intensity (Fig. 3a) whilst remaining within the De value plateau (Fig.   3b).
Preheating for short durations resulted in the heater unit overshooting the target temperature and poor thermal reproducibility. 25 For this reason, a longer duration preheat at 160 °C for four minutes was selected for all measurements. This selected protocol is further validated by the successful recovery of a 360 Gy dose from naturally zero-age sample KRG16-112 for both the Al and Ti-centres, which yield recovered to given dose ratios of 0.83 ± 0.20 and 1.01 ± 0.06 respectively (n=3).
Measurements of the trapped charge population of the Al and Ti-centres were similar between aliquots resulting in 1s 30 uncertainties of ~20%. Equivalent dose values for the Al-and Ti-centres were within uncertainty for all samples, and ages ranged from 291 ± 13 ka for sample KRG16-05 to 37 ± 2 ka for sample KRG16-104 (Table 1). Samples KRG16-111 and KRG16-112 from the high-temperature tunnel yielded zero age; consequently, full dose response and isothermal decay was not measured for sample KRG16-112 and it is not included in Table 1. Whereas it was possible to saturate the Ti-centre of all samples with the maximum given dose of 19 kGy, the Al-centre continued to grow linearly throughout measurement for all samples (Fig. 4a,c). Continued growth of the Al-centre has been reported previously and has been accommodated through fitting dose response with an exponential plus linear function (e.g. Duval, 2012). In contrast, for the KRG samples, the Al-5 centre is best described using a linear regression (Fig. 4a). In contrast, the Ti-centre of all samples showed a reduction in signal intensity at high doses (i.e. >10 kGy), which has also been reported previously (e.g. Duval and Guilarte, 2015) and has been attributed to changing electron capture probabilities (Woda and Wagner, 2007). To characterise the maximum possible trapped-charge population we excluded data points where the ESR signal intensity started to reduce (white data points in Fig.   4c) and fitted the remaining data with a single saturating exponential function (e.g. Grün and Rhodes, 1991) of the form: 10 where I is the natural ESR signal intensity, Isat is the saturation intensity of the ESR signal and D is the given dose (Gy).
Because the Ti-centre experiences saturation, the equivalent dose value, ' , can also be expressed as a saturation ratio, i.e. 15 $ = ( / ). As only a single aliquot of each sample was dosed until saturation, $ values for the Ti-centre were calculated from interpolation of the average De value (n=3) onto the single dose response curve, thus $ values for the Ti-centre are derived from multiple aliquots.
Toyoda and Ikeya (1991) suggested that the thermal decay of the E', Al and Ti-centres follows second order kinetics, however 20 it was not possible to fit our data using either a first or second order kinetic model (Supplementary Material). Instead the isothermal decay data were fitted using a multiple first order kinetic model (Lambert, 2018; Table 1). Whilst the actual physical meaning of a Gaussian distribution of energies requires further investigation within the context of ESR defects, and it is unlikely that both the Al and Ti-centres follow exactly the same process of thermal decay, preliminary fits to the data using this model are promising (Fig. 4b,d). Values of µ(Et) ranged from 1.3 -1.9 eV between samples and centres (Table 1). 25

Optically Stimulated Luminescence
For all of the samples, the measured luminescence signals fulfilled the acceptance criteria (see Supplementary Material for 30 further details). The IRSL50 signals of all samples exhibited very high rates of fading with g2days values ranging from 6-11 %/decade, whereas for post-IR IRSL measurements at 225 °C, fading rates were 2-4 %/decade. The model introduced by Huntley (2006) was used to fit the athermal detrapping data to determine r'. Using this model to fading correct the trapped charge concentrations following Kars et al. (2008) indicates that the IRSL50 signals of samples KRG16-06 and KRG16-101, and all signals for sample KRG16-05 are saturated (see Supplementary Material). All other signals can be used to determine rock-cooling histories. Saturation of the IRSL50 signals relative to the higher temperature signals is a consequence of their relatively high rate of anomalous fading. The luminescence dose response data of all of the samples and signals were fitted 5 with a single saturating exponential fit to determine the characteristic dose of saturation, D0 and the concentration of trapped charge, $. Although for some samples a general order kinetic model (GOK) fit would result in lower deviation from the measured values, GOK fits have been shown to overestimate sample athermal field saturation values , which must be done accurately to evaluate if a sample contains thermal information (cf. Valla et al., 2016). Finally, the isothermal decay data were fitted using the band-tail states model (Poolton et al., 2009;Li and Li, 2013;Eq. 6 and 7) to determine Et, Eu 10 and s. Values of Et ranged from 1.2 -1.5 eV between samples and signals ( Table 2).

Inversion of ESR and OSL data for cooling histories
In order to invert the data into cooling histories, we used the same approach outlined in section 3. We computed OSL and ESR dose response curves from 10,000 randomly generated t-T paths, which were constrained to cool monotonically between 200 15 °C and 15 ± 5 °C, over 2 Ma. Initially the Al-centre, Ti-centre and OSL centres were inverted separately, before being inverted together. The results for all samples with the exception of naturally zero-age samples

Discussion
Trapped-charge thermochronometers offer benefits over other thermochronometry systems because of their low closure temperature and ability to yield precise cooling histories over Quaternary timescales . However, signal saturation has proven a significant barrier to the application of luminescence thermochronometry (cf. Valla et al., 2016).
For ESR thermochronometry to offer a viable alternative it should exhibit later signal saturation but also similar thermal 25 stability. The measurements presented here are promising because whilst the ages measured for the OSL and ESR systems are similar, the maximum possible ages that can be obtained from the ESR Ti-centre are more than four times greater than the maximum possible age that can be obtained from the OSL signals (Tables 1 and 2). Furthermore, the Al-centre of the KRG samples does not exhibit signal saturation up to 19 kGy, which was the maximum dose explored in this study (Fig. 4a).
Although such linear dose response behaviour has, to our knowledge, not been reported previously and thus may be a property 30 of these exceptionally young quartz minerals, it is an exciting observation that warrants further study through the investigation of further quartz samples. Samples KRG16-111 and KRG16-112 from the high temperature tunnel yielded zero, or near-zero ages for both ESR centres and the IRSL signals investigated (Tables 1 and 2). These samples provide an important local control on the thermal stability of these trapped-charge systems, demonstrating that all charge is evicted from the centres at sufficiently high temperatures.
For the remaining samples, the ages obtained from the two ESR centres are within uncertainties, indicating that they may have 5 similar thermal stability. For the OSL data, some variance in age is recorded between the different signals (Table 2); all signals of sample KRG16-05 are in field saturation, and thus only a minimum sample age of ~180 ka can be calculated ( Table 2). The IRSL50 signals of samples KRG16-06 and KRG16-101 are also in field saturation, yielding the highest apparent ages for these samples, and are not considered further. For samples KRG16-06 and KRG16-101, the remaining IRSL signals show a general reduction in age with increasing stimulation temperature, possibly indicating that the ages have been overcorrected for 10 anomalous fading using the Huntley (2006) model (cf. King et al., 2018. The OSL and ESR ages of samples KRG16-06 and KRG16-104 are similar, indicating that for these samples the ESR and OSL signals have similar thermal stability, and thus that ESR-thermochronometry would also be suitable for resolving late stage cooling histories. In contrast, sample KRG16-101 yields OSL ages twice as large as the ESR ages, which could be indicative of a difference in centre thermal stability.

15
To further evaluate the relative thermal stability of the ESR and OSL signals, the isothermal decay of the ESR and feldspar systems was simulated using the experimentally constrained kinetic parameters of the different samples, for isothermal conditions of 20 °C assuming an initial trapped charge concentration of 1 and assuming no charge trapping (Fig. 6). Note that anomalous fading related signal loss has also been included for the OSL signals, as excluding this variable would result in erroneously high apparent signal stabilities. The ESR centres have similar thermal stability to the IRSL centres for all samples, 20 with the exception of sample  where the ESR centres are more thermally stable. The Ti-centre is more thermally stable than the Al-centre for all samples, with the exception of sample KRG16-104 for the measurement in response to 2.15 kGy. This is consistent with the earlier work of Grün et al. (1999) who also extracted quartz from granitic bedrock and observed that the Ti-centre is more thermally stable than the Al-centre, but contrasts with observations from Chinese loess (Tsukamoto et al., 2018). The contrasting behaviour between the two measurements of KRG16-104 in response to doses of 25 4.30 kGy and 2.15 kGy (Fig. 6d) reflects uncertainty in the derivation of ESR kinetic parameters, potentially related to interaliquot variability. Improved measurement protocols and the development of automated instrumentation may alleviate these discrepancies through improving measurement reproducibility, however despite this, the thermal stability determined in both experiments is broadly similar (Fig. 6d). The general trend of ESR signals exhibiting similar thermal stability to IRSL signals indicates that ESR-thermochronometry will record changes in exhumation histories from a similar thermal range as OSL-30 thermochronometry, whilst benefitting from considerably later signal saturation and being unaffected by anomalous fading.

Figure 6
Inverting the Al and Ti-centres of all samples results in broadly similar time-temperature histories between centres, whist the cooling histories of different samples vary (Fig. 5). This is in agreement with the Al and Ti-centres' similar thermal stabilities ( Fig. 6) and measured ages (Table 1). The two different centres can also be effectively combined to produce a single cooling history (Fig. 5), which is similar to that inverted from the OSL data alone for samples KRG16-101 and KRG16-104 (Fig. 5).
For sample KRG16-05, the saturated OSL signals result in a broad cooling history, whereas for sample KRG16-06 the OSL 5 data yield more rapid cooling than the ESR data. The OSL and ESR data can also be inverted together. These data show that for samples beyond the range of OSL dating, ESR-thermochronometry will be able to provide cooling histories over a similar thermal range, allowing late stage exhumation histories to be determined. The data inversions reveal that rates of rock cooling in the Hida range of the Japanese Alps are consistent with previous investigations that indicate rapid rock cooling (Ito et al., 2013;2017;Spencer et al., 2019). Whereas sample KRG16-05 experienced almost no cooling over the past 2 Myr, cooling 10 rates accelerated from ~100 °C/Myr (calculated from the U/Pb ages of Ito et al., 2013) to rates of >400°C/Myr over the past 100 ka for samples KRG16-06, KRG16-101 and KRG16-104.

Conclusions and Outlook
In this study, the potential of ESR thermochronometry for constraining rates of rock cooling has been explored for a suite of samples from the Hida range of the Japanese Alps. Through using the latest ESR measurement protocols (Tsukamoto et al., 15 2015) and instrumentation (Oppermann and Tsukamoto, 2015) the dose response and thermal stability of both the Al and Ticentres has been constrained. Whilst the Ti-centre can be described with a single saturating dose response curve, the Al-centre continues to grow linearly with laboratory irradiation. A multiple-first order model based on a distribution of trap-depths was successfully used to fit isothermal decay data (Lambert, 2018), which do not follow either simple first-order, or second-order decay. Contrasting the thermal stability of the Al and Ti-centres with that of the luminescence centres of feldspar shows that 20 the ESR of quartz has similar thermal stability. The Al and Ti-centres can be successfully inverted together for rock cooling for all of the samples investigated. It was also possible to invert the OSL and ESR data together for all samples analysed, providing further constraints on their thermal histories. Whereas OSL-thermochronometry of sample KRG16-06 can only recover a minimum cooling rate of ~160 °C/Myr, both ESR centres have the potential to recover cooling rates of as low as 50-25 °C/Myr, illustrating the potential of ESR for resolving late-stage cooling histories . 25 Tissoux, H., Falguères, C., Voinchet, P., Toyoda, S., Bahain, J. J., and Despriée, J., 2007. Potential use of Ti-center in ESR dating of fluvial sediment. Quaternary Geochronology, 2(1-4), 367-372. Townend, J., Zoback, M.D., 2006. Stress, strain, and mountain building in central Japan. Journal of Geophysical Research: Solid Earth (1978-2012) 111. Toyoda, S. and Ikeya, M., 1991. Thermal stabilities of paramagnetic defect and impurity centers in quartz: Basis for ESR 5 dating of thermal history. Geochemical Journal, 25(6), pp.437-445. Tsukamoto, S., Toyoda, S., Tani, A., & Oppermann, F., 2015. Single aliquot regenerative dose method for ESR dating using X-ray irradiation and preheat. Radiation Measurements, 81, 9-15. Tsukamoto, S., Long, H., Richter, M., Li, Y., King, G.E., He, Z., Yang, L., Zhang, J. and Lambert, R., 2018. Quartz natural and laboratory ESR dose response curves: A first attempt from Chinese loess. Radiation Measurements. 10 Valla, P.G., Lowick, S.E., Herman, F., Champagnac, J.D., Steer, P. and Guralnik, B., 2016. Exploring