Technical note: Quantifying uranium-series disequilibrium in natural samples for dosimetric dating – Part 1: gamma spectrometry

. Dosimetric dating techniques rely on accurate and precise determination of environmental radioactivity. Gamma spectrometry is the method of choice for determining the activity of 238 U, 232 Th, and 40 K. With the aim to standardize gamma-spectrometric procedures for the purpose of determining accurate parent nuclide activities in natural samples, we outline the basics of gamma spectrometry and practical laboratory procedures here. This includes gamma radiation and instrumentation, sample preparation, ﬁnding the suitable measurement geometry and sample size for a given detector, and using the most suitable energy peaks in a gamma spectrum. The issue of correct efﬁciency calibration is highlighted. The procedures outlined are required for estimating contemporary parent nuclide activity. For estimating changing activities during burial speciﬁc data analyses are required, and these are also highlighted.


Introduction
Dosimetric dating is a method for determining the age of objects from the absorbed dose to their constituent crystals caused by exposure to radiation, typically from natural radionuclides such as uranium, thorium, and 40 K present in most soils and sediments. It is of critical importance to the dating technique that the quantity of the radiation dose is estimated accurately, and this is usually the case when its rate is constant. Because of radioactive decay the dose rate will in many situations vary over time, and the most important source for this variation is the 238 U-series radionuclides, in particular the short-lived daughters of the intermediate radionuclide 226 Ra. Resolution of any initial disequilibrium between 226 Ra and its parents will normally take place on millennial timescales and therefore result in significant changes to the dose rate over time. The simplest and most straightforward method for determining activities of 226 Ra, its parents, and its short-lived daughters is gamma spectrometry (e.g. Murray and Aitken, 1988). Although the initial set-up can be quite demanding, because of its operational simplicity this is increasingly the method of choice (e.g. Murray et al., 2015), particularly where there is evidence of radiometric disequilibrium.
Despite operational simplicity, the interlaboratory comparison (Murray et al., 2015) reveals discrepancies between laboratories that are of great concern. While the exact reasons are unknown, we suspect that problems arise from data analysis, peak interference and self-absorption, and the use of suitable reference material (e.g. Murray et al., 2018) for efficiency calibration. We address these problems in the first part of our paper wherein we outline the basic principles of gamma spectrometry for estimating contemporary dose rates. We show how the degree of isotopic fractionation can be quantified using the approach published in Abdualhadi et al. (2018). For estimating the historical dose rate the radioactive decay equations (the Bateman equations) for the intervening period must be solved (e.g. Degering and Degering, 2020). Analytical solutions of these equations, as well as software procedures for calculating historical dose rates, will be the subject of the second part of our paper.

Gamma radiation and gamma spectrum
Ionizing radiation from atoms undergoing radioactive decay includes both subatomic particles (alpha rays, beta rays) and electromagnetic waves from the higher end of the energy spectrum (gamma rays, X-rays). Gamma rays are emitted when the decay product undergoes de-excitation from a higher energy level to a lower one or to the ground state and are characterized by photon energies in the range 20 keV to 3 MeV. Most of the radionuclides relevant for dosimetric dating emit gamma rays of energies < 1 MeV. Gamma rays can penetrate material objects typically to a depth of several centimetres depending on the photon energy and characteristics of the material (in particular its attenuation coefficient), though with an intensity that attenuates exponentially with distance. For example, the intensity of a 1 MeV photon penetrating SiO 2 is reduced to 50 % at a penetration depth of 4.7 cm; the intensity is reduced to 63 % at 6.8 cm. Values of the attenuation coefficients at other photon energies or for other materials or mixtures can be found in the literature or from online databases (for details see Appendix A).
Although gamma rays make only a small contribution to the radiation dose compared to alpha and beta particles, they play an important role in determining the source and strength of the dose. In contrast to beta particles, gamma photon energies are discrete and specific to particular nuclides. This allows determination of gamma-emitting radionuclides contained within a source from peaks in a gamma-ray energy spectrum ( Fig. 1) obtained from that source.

Detector and instrumentation
Gamma-ray emissions from environmental samples are readily analysed using semiconductor gamma detectors. Electrons displaced by photons traversing the crystal at the heart of the detector are gathered using an applied high voltage and the total charge collected related to the photon energy. These devices are most commonly constructed using crystals made of high-purity germanium (HPGe). Nearly perfect single crystals of Ge can be produced in a range of different sizes and configurations suitable for analysing environmental samples (for details see Gilmore, 2008). Ge detectors are operated at low temperature to reduce the background current, i.e. charge carriers that are generated due to the relatively low band gap (0.7 eV) of germanium. The typical installation of an HPGe detector includes the detector and preamplifier within the cryostat housing and the nitrogen dewar ( Fig. 2A), the shield surrounding the cryostat (Fig. 2B), the amplifier, and the multichannel analyser; the latter two are currently being increasingly replaced by a digital system. The energy range is typically 10-2600 keV. Table 1 compares detector specifications. The BEGe (broad-energy germanium) detector has a short length and a different electrode geometry compared to the nand p-type coaxial detectors. It ensures good charge collection and excellent low-energy resolution. On all detectors types, samples are placed close to the end cap to ensure good interception of gamma emissions from the sample. At low energy (< 200 keV; Abdualhadi et al., 2018) they have a similar efficiency for the same detector diameter, though the better energy resolution of the BEGe detector allows a more precise determination of the peak areas than the coaxial detector. This is important when the low-energy gamma rays are a vital part of the analysis. At higher energies the coaxial detector has a significantly better efficiency due to its longer length. The BEGe detector, on the other hand, typically has a lower Compton background in the spectrum as there is less material to scatter high-energy gamma rays such as 1461 keV from 40 K and 2614 keV from the thorium decay series. . Typical configuration of a gamma spectrometer with (a) detector housing, nitrogen dewar, and PC. (b) View in the detector housing; the shield is composed of 10 cm low-activity lead, 1 mm tin (not visible), and 1 mm copper. Amplifier and preamplifier not shown (adopted from Abdualhadi, 2016).

Interaction between detector and sample
There are interactions occurring between the detector and the sample's gamma rays that are independent of detector configuration and measurement geometry but disturb the recording of the sample's activity. These interactions occur across the energy range of interest with variable significance depending on the average atomic number (Z) of a given sample.
In the low-energy region a photon may be absorbed, and its energy ejects an electron from its host atom (photoelectric effect). In the intervening energy region the incoming photon is scattered by an electron, resulting in a decrease in the photon's energy (Compton scattering). Photons with energies of 1022 keV can cause the formation of an electron and a positron (pair production).

Full-energy peak efficiency
Gamma photons entering the detector crystal may be fully absorbed, partially absorbed, or escape interaction altogether.
Only those photons that are fully absorbed will contribute to the relevant photopeak. Full-energy peak efficiency is defined as the ratio of the number of counts detected in a peak to the number emitted by the source. Its value for any given photon energy will depend on the properties of the detector and the configuration of the sample. Figure 3 shows plots of efficiency versus photon energy for two of the detectors listed in Table 1. These were determined using standard sources of known activity. While both detectors have similar efficiency at low energy, the coaxial detector has a significantly higher efficiency as the gamma-ray energy increases due to its longer length. In both detectors, the practical efficiency at low gamma-ray energies can be reduced by other effects such as self-absorption of photons within the sample (see Sect. 6.2).

Sample preparation and measurement geometries
Samples are normally dried and then placed in a suitable sample holder. There are three common measurement geometries: Marinelli beaker, cylindrical pot, and in the case of a well-type detector, sample tubes designed to fit snugly in the well (see Appendix B, Fig. B1). To obtain the texture suitable for the relevant measurement geometry samples can be pulverized, homogenized, and placed in a geometry that matches that of the standards used to calibrate the detector. One method for achieving this is to mix the dried and powdered sample with wax before casting it in a shape matching the required geometry (Fig. B2). Since some radionuclides may be preferentially selected by certain components within the sample, any treatment should carefully maintain the original composition. Removing a component, for instance or-  Table 1 for a description of BEGe and coaxial (n) of this study.
ganic matter, may generate inaccurate data because different components may be in different stages of equilibrium depending on their geochemical history during burial. If 226 Ra is being determined via its short-lived daughter 214 Pb, samples will need to be sealed to prevent escape of the intermediate radionuclide 222 Rn (e.g. Tudyka et al., 2021) and stored for ∼ 25 d to ensure 226 Ra-222 Rn-214 Pb radioactive equilibrium.
In practice, the sample configuration will be dictated by the detector type, the amount of available material, and the need to minimize self-absorption (see Sect. 6.2). Although near 4π geometry makes Marinelli beakers (ca 300 ml volume; see Fig. B1c) the preferred configuration, this will only be practicable when there is sufficient material to completely fill the beaker. When the amount of material is limited (< 100 g), for a coaxial or BEGe detector the most suitable geometry will be cylindrical32 ( Fig. B1a), with a diameter comparable to that of the detector. Well-type tubes (Fig. B1b) can normally only accommodate very small samples, typically no more than a few grams, and are only suitable for use on well-type detectors and when larger amounts of material are not available.

Data analysis
Calculations of radionuclide activities need to take into account emission probabilities (see Table 2 for values of 238 Useries radionuclides and Table 3 for the 232 Th series), fullenergy peak efficiency (e.g. Fig. 3), geometry-dependent true coincidence summing (see Gilmore, 2008, for details), background characteristics of the detector (see Appendix C for details on background), interactions between sample and de-tector, and measurement geometry. Table 2 lists photon energies and emission probabilities for key members of the 238 U series (see Guibert et al., 1994, for comparison). Table 3 lists corresponding values for key members of the 232 Th series. For radionuclides with multiple energy lines, a best value of the activity may be obtained by using a weighted sum of activities determined from each of the most suitable lines. The other main natural source, 40 K, emits 1460.8 keV photons with a probability of 10.66 ± 0.16 %. The peak area given by most software will include background counts as well as those due to emissions from the sample itself. The number of counts in a peak is most easily measured when the peak is well resolved, stands alone free from the interference of any neighbouring peak, and has a small background contribution. Methods for calculating the peak area may include a simple summation of counts in channels within the region of interest or a procedure such as the Levenberg-Marquardt non-linear least squares fit algorithm for fitting the peak to a prescribed shape.
In routine measurements peak interferences and non-linear peak efficiencies due to self-absorption can be minimized by (i) limiting sample size (i.e. height of fill in the sample holder, Fig. B1), (ii) using a reference material of known activity (see Sect. 5.1 and Appendix C) having exactly the same geometry as the unknown sample and similar self-absorption characteristics, and (iii) using the most reliable energy lines for determining activities (Tables 2 and 3).

Efficiency calibration
For environmental samples the calibration method of choice is comparing the unknown sample with a sample of known activity concentration. It is of vital importance that known and unknown samples are similar in terms of atomic number, mass, and density and that the measurement geometry is kept constant. This straightforward approach is compromised insofar as available certified reference material (e.g. IAEA, NRCAN; Murray et al., 2018) is typically not ideal in terms of its properties. As a consequence, non-certified internal standard material specifically established for dosimetric dating is often in use (e.g. De Corte et al., 2007;Preusser and Kasper, 2001). Alternatively, certified material is deployed in concert with a number of variables that account for emission probability, full-energy peak efficiency, and, eventually, density or chemical composition. This procedure is specific to the detector and can be established through a series of experiments, or it is carried out through software packages provided by the manufacturer (e.g. LabSOCS, ANGLE). For details see Appendix C.
The interlaboratory comparison study published by Murray et al. (2015) shows considerable differences for gammaspectrometric results obtained for 238 U and 232 Th. Because U and Th constitute around 30 % of the total annual dose rate (Aitken, 1985), these differences are of great concern. While some of the differences may arise from procedures carried Table 2. Key gamma-emitting 238 U-series radionuclides, their emission energies (Eγ ), and emission probabilities (P γ ). The 234 Th peak at 63.3 keV is a doubleton that includes emissions at two different energy levels, 62.9 keV (0.016 %) and 63.3 keV (3.7 %). The peak at 92.6 keV includes emission at three different energy levels: 92.3 keV (0.017 %), 92.4 keV (2.12 %), and 92.8 keV (2.10 %). Data are from the National Nuclear Data Centre (https://www.nndc.bnl.gov/, last access: 10 August 2021). out in individual laboratories, the efficiency calibration is an issue affecting all laboratories. An empirical study (Mauz et al., 2021) suggests that the calibration method contributes < 10 % difference for the activity data of 238 U, 232 Th, and 40 K when unsuitable energy peaks (e.g. 186 keV) are excluded.

Uncertainties
The uncertainty in the overall result depends on a range of factors. These include counting statistics, detection efficiency, nuclear decay data, sample composition, and geometry-defined true coincidence summing. For the typical activities to be measured and the sample sizes used the uncertainty arising from counting statistics is typically 2 %-5 % for a counting time of around 1 d. This depends on the intensities of the various gamma rays used in the analysis. The detection efficiency depends on the counting geometry and the detector used, and it also varies with gamma-ray energy. By using a series of different reference materials to fully characterize the systems used an uncertainty of 2 %-3 % can be achieved. The emission probabilities obtained from the nuclear decay data have uncertainties. Tables 2 and 3 show that the majority of these are small. The exception is the decay of 234 Th. Here the two gamma rays have relatively large uncertainties. If both can be used in the analysis then the uncertainty for 234 Th is around 5 %. The sample composition can lead to self-absorption, particularly for low-energy gamma rays. When the density of the sample is known as is its elemental composition, the amount of self-absorption can be determined accurately with an uncertainty of only a few percent. The effect of true coincidence summing also needs to be taken into account. This depends on the decay scheme of the nuclei involved and the detection efficiency of the system being used. For most of the nuclei in the uranium and thorium decay chains this effect is very small for the efficiencies that are used in the majority of counting geometries, and it is not a major contributor to un-certainty. For more details see Abdualhadi (2016). Overall, uncertainties of 3 %-6 % can be achieved.

Quantifying U-series radionuclides for the detection of secular disequilibrium
In samples that have effectively been sealed for long periods of time, members of the 238 U series will all be in radioactive equilibrium: that is, they will all have the same activity. In this case the activity can in principle be determined using any of the nine photopeaks listed in Table 2. In practice, mobilization of intermediate members is likely to result in the U series being in disequilibrium. When this involves short-lived daughters, disequilibrium with the nearest long-lived parent will also be short-lived. For example, 226 Ra will usually be in equilibrium with its short-lived daughters 214 Pb and 214 Bi. Disequilibria between long-lived members such as 230 Th (half-life 75380 years) and 226 Ra (half-life 1600 years) will take many thousands of years to resolve. In these cases the activity must be measured using photopeaks from the radionuclide itself or its short-lived daughters. Because secular disequilibrium is typically caused by the mobilization of 234 U and 226 Ra, the following energies are key to detect radioactive disturbance: 63 and 92 keV (to infer 234 U through 234 Th), 186 keV (to quantify 226 Ra), and 46 keV (to quantify the endmember 210 Pb). There are nontrivial issues associated with these photopeaks that need to be solved in order to determine the activity of 234 Th, 226 Ra, and 210 Pb (see Sect. 6.1 and 6.2). A major obstacle in determining the complete 238 U series is the absence of significant gamma emissions from 230 Th. In samples that have been isolated on timescales of tens of thousands of years, 230 Th will be in secular equilibrium with its daughter radionuclide 226 Ra and can thus be determined via the 226 Ra emission. On shorter timescales, potential disequilibria between 230 Th and 226 Ra are a significant source of uncertainty.

Correcting peak interference
At 92 keV the gamma peak of 234 Th is interfered by X-rays and by other gamma peaks including the 93.3 keV thorium line (Huy and Luyen, 2005; see Appendix D for details). In addition, for Ge detectors there is a contribution below the photopeak from multiple Compton scattering whereby the interaction probability ratio for the photoelectric effect and Compton scattering is detector-dependent. There is also background from Comptons from higher-energy peaks so that, effectively, the background at 92-93 keV is a continuum with a small step at each peak (Fig. 4). It is good practice to first collect some spectra from Th-rich samples to define the position and characteristics of the peaks near 92 keV so that energies and FWHMs can be fixed when deconvoluting the peak region (Fig. 4). All these non-trivial issues may rule out the use of the 92 keV peak for quantifying 234 Th, leaving us with its 63 keV peak that, however, over-  laps with a 232 Th decay peak at 63.8 keV (P γ = 0.263 %; see also https://www.nndc.bnl.gov/, last access: 10 August 2021). For the purpose of data robustness both peaks are, ideally, analysed and compared (Fig. 5). This process can be more straightforward if a detector with excellent low-energy resolution is used. At 186 keV the gamma peaks of 226 Ra and 235 U overlap. There are three ways to handle this.
(1) The first is subtracting the 235 U portion from the total activity deduced from the 186 keV peak. This is achieved by calculating the natural isotopic ratio of 238 U and 235 U by taking into account atomic mass, Avogadro constant, decay rate, and half-life of each radionuclide. This results in 58 % of the counts in the 186 keV peak emitted by 226 Ra and 42 % emitted by 235 U. Assuming secular equilibrium between 238 U and 226 Ra these percentages are constant and can be used as a nominal factor to correct the interference at 186 keV.
(2) The second is calculating the 226 Ra activity by including other gamma peaks (e.g. de Corte et al., 2005). With this method the 235 U activity is directly inferred from its peak at 144 keV or indirectly from the 234 Th peak at 63 keV. The method requires correction of peak interference at 63 keV and of true coincidence summing as well as adoption of the natural isotopic ratio of 238 U and 235 U. Its advantage is that secular equilibrium is assumed only between 238 U and its immediate daughter 234 Th.
(3) The third is bypassing correction procedures by sealing the sample against 222 Rn escape (as indicated in Sect. 4) and determining 226 Ra via the emissions from its shortlived daughters 214 Pb (352 keV) and/or 214 Bi (609 keV). Note, however, that 214 Bi emissions can be significantly affected by coincidence summing.
We have compared approaches (2) and (3). The results are ambiguous (Fig. C1), suggesting that correction (2) is insufficient to accurately determine the 226 Ra activity (see Appendix D for details).

Correcting self-absorption
A fraction of gamma photons generated by radioactive decay are lost by absorption within the sample before reaching the detector. The extent of the losses will be controlled by the sample mass, measurement geometry, and the gammaray attenuation coefficient of materials within the sample. They will generally be relatively small at high photon energies but much more significant at low photon energies. A semi-empirical formula for estimating the losses is where N is the number of photons actually detected, N 0 the number of photons that would have been detected in the absence of self-absorption, m the mass of the sample, the material gamma-ray mass attenuation coefficient, and k a parameter characterizing the geometry of the sample (Appleby et al., 1992). Values of the attenuation coefficient can be determined empirically using the Beer-Lambert law (Eq. A1) by measuring the reduction in intensity when a collimated beam of gamma rays is passed through the sample (transmission method). Table 4 lists empirically determined intensity reduction values (I /I 0 ) for a number of samples with different densities and chemical properties at photon energies in the range 32 to 121 keV. Differences between samples are relatively small at photon energies greater than 60 keV but may be significant at low photon energies. More advanced approaches are available using the Monte Carlo N-particle (MCNP) transport code, which effectively simulates the interaction of gamma rays emitted in a random direction by a sample of known density and chemical composition with the detector (for details see Abdualhadi et al., 2018). When the sample has a known chemical composition estimates of the attenuation coefficient can be made using information available in published databases (see Appendix A).

Summary
This is Part 1 of our paper on the quantification of uraniumseries disequilibrium in natural samples. It focuses on basic principles for estimating contemporary dose rates using gamma spectrometry. We describe the basics of gamma radiation and the penetration depth of this radiation in matter as a function of photon energy and material properties. The detectors suitable for detecting gamma-ray emissions in the typical range of 10-2600 keV (Fig. 1) are high-purity germanium (HPGe) crystals operated at low temperature (Table 1). The fully absorbed gamma photons entering the detector crystal generate discrete full-energy peaks in the gamma spectrum that are specific to 40 K and the radioisotopes of the 238 U and 232 Th series. The analysis of these full-energy peaks accounts for the gamma emission probability of the corresponding radioisotope, the detector's efficiency and background at that energy, and accompanying interactions between the sample and detector. All these factors together result in the choice of a number of photopeaks suitable for estimating the parent nuclide activities (Tables 2 and 3). For these peaks the simple comparison with a reference material of known activity (e.g. Volkegem loess) delivers accurate activities of individual radioisotopes and hence the parent nuclide, with an uncertainty of 3 %-6 %.
The essential prerequisite of this approach is radioactive equilibrium: that is, all members of the 238 U series have the same activity. Since only samples that have been isolated for millions of years will be in secular equilibrium, most samples that are subject to dosimetric dating are likely to exhibit disequilibrium. The significance of this disequilibrium for the dose-rate estimate depends on the radioisotope involved and the duration of the process that leads to loss or gain of a 238 U-series member. If disequilibrium is suspected, the photopeaks at 46, 63, and 92 keV as well as at 186 keV have to be analysed, and this involves accounting for peak interference and self-absorption, both being well-known issues of gamma spectrometry. Advanced approaches are currently available using the Monte Carlo N-particle (MCNP) transport code, Table 4. Transmitted intensities of gamma rays for sandy samples with different densities (and chemical composition). I /I 0 is the ratio between experimentally measured gamma-ray intensity (I ) and gamma-ray intensity without a sample (I 0 ). Gamma energies (E γ ; keV) of point sources are 32 ( 137 Cs), 39 ( 152 Eu), 59 ( 241 Am), 80 ( 133 Ba), and 121 ( 152 Eu) (Abdualhadi, 2016). which effectively simulates the interaction between the detector and gamma rays emitted in a random direction by a sample of known density and chemical composition. This allows quantifying the activities of key members of the 238 U series required to calculate the historical dose rate by solving the Bateman equation. This will the subject of Part 2 of our paper.

Appendix A: Penetration depth and attenuation of gamma rays
Depending on the photon energy and characteristics of the material, gamma rays can penetrate objects to a depth of several centimetres, though the intensity (I ) attenuates exponentially with distance (x) in accordance with the Beer-Lambert law.
The (linear) attenuation coefficient µ is a function of photon energy and material characteristics. The penetration depth is defined as the distance over which the intensity is reduced to a fraction 1/e (37 %) of its initial value I 0 (that is, a 63 % reduction) and is given by The intensity is reduced by 50 % over a depth ln(2)/µ = 0.69 δ. Values of the mass attenuation coefficients given in the literature are usually expressed in terms of the mass attenuation coefficient, defined as where ρ is the mass density. For example, the mass attenuation coefficient for 1 MeV photons in SiO 2 is 0.0637 cm 2 g −1 ; since the density is 2.32 g cm −3 , the linear attenuation coefficient is 0.148 cm −1 . The penetration depth (63 % reduction) is thus 6.8 cm and is 4.7 cm for a 50 % reduction. Values of at other photon energies or for other materials can be found in the literature (e.g. Hubbell, 1983) or from online databases (Hubbell and Seltzer, 2004).

Appendix C: Background, efficiency calibration, and reference material
The gamma spectrum of a given sample records emissions from both the sample and the background environment (including the detector itself). They include (a) discrete fullenergy peaks at photon energies specific to individual radionuclides (whether from the sample or from the background environment) and (b) a more or less continuous spectrum of partial-energy interactions arising from various scattering processes. Most software will calculate both the total number of interactions within the boundaries of each photopeak (the gross area, including those due to scattering) and the number of full-energy interactions (the net area) by subtracting those due to scattering from the gross area. The contribution of the background environment to the net peak area can be determined from "background counts" carried out without a sample being present. Writing N T for the net counts in the peak with the sample present and N B for the background contribution, the number attributable to the sample itself will be N = N T − N B . The activity (A) of the radionuclide can be calculated using the formula where m is the mass of the sample, T the count time, P γ the gamma emission probability, and η a practical efficiency that takes account of not only the full energy efficiency of the detector at that energy, but also the effects of factors such as self-absorption and, where necessary, true coincidence summing (Gilmore, 2008). The effects of self-absorption can be minimized by limiting the sample size. A procedure that avoids many difficulties is using a reference material of known activity with similar self-absorption characteristics as that of the sample, having exactly the same geometry. Activities in the sample can then be calculated by direct comparison with the reference material using the equation Equation (C2) is valid on a peak-by-peak basis, which requires A reference to be known for each peak of interest. Alternatively, if A reference is only known for 238 U, all members of the 238 U series in both the reference material and unknown sample have to be in secular equilibrium. In this case the activity can in principle be determined using any of the photopeaks listed in Table 2. In practice some peaks are more suitable then others depending on the sensitivity of the peak to the sample's Z and density as well as on detector-specific parameters such as energy resolution and full-energy peak efficiency.

Non-certified internal standard material
Two loess samples, namely Volkegem (De Corte et al., 2007) and Nussi (Preusser and Kasper, 2001), have previously been proposed as suitable standard material. Loess is a sedimentary deposit resulting from the accumulation of wind-blown dust. Typically, a loess deposit is several metres thick and composed of 20-40 µm grains. The grains comprise minerals such as feldspars, quartz, and heavy minerals in percentages representative for the provenance areas. Loess is mostly transported over long distances and over large areas, hence carrying the signature of a significant proportion of the Earth's surface, which is then homogeneously distributed in the accumulation. In fact, while grain size may vary between outcrops, radionuclide activity concentrations are strikingly similar (see e.g. Scheidt et al., 2021, for the Carpathian loess; Roberts et al., 2003, for North American loess;Sauer et al., 2016, for central European loess;Lü et al., 2020, for Chinese loess).
Volkegem -this material is derived from a loess deposit in Belgium (De Corte et al., 2007).
Nussi -this material is derived from a loess deposit in southern Germany (Preusser and Kasper, 2001). The material is registered with the International Association of Geoanalysts (http://www.geoanalyst.org, last access: 22 February 2018) as "Loess-1 GeoPT13" with assigned chemical values for all low-mass elements (Potts et al., 2003).
Both samples have been tested in several luminescence laboratories. Here we provide details for the Volkegem material -see Table C1. Table C2 lists a number of materials commonly used for efficiency calibration. It should be noted that some of these materials entail problems; for instance, some radionuclides included in the QCYK multi-nuclide reference material exhibit true coincidence summing when measured close to the detector (Gilmore, 2008), and this limits the usage of some energy lines.

Certified reference material
We compared two different efficiency calibration methods using the coaxial (n) detector of this study described in Table 1: one uses the multi-nuclide reference solution QCYB41 and follows Eq. (C1). The other one is based on the direct comparison with the Volkegem reference material using Eq. (C2). A total of 81 samples with different activities originating from geochemically variable environments were used for the comparative study. Almost perfect agreement exists for energy lines 1460, 338, 352, 609, 911, and 46 keV. Systematic differences appear for energy lines 63 and 238 keV (Mauz et al., 2021). These results suggest that accurate activity data are obtained from the simple comparison with the Ortec ANGLE user guide: https://www.ortec-online. com/ (last access: 3 April 2022).

Appendix D: Peak interference
Overlaps between neighbouring peaks will occur where photopeaks are separated by less than around 2.5 times the FWHM value: that is, between 4 and 5 keV. An example is the 226 Ra peak at 186.2 keV, which is separated from the nearby 235 U peak at 185.7 keV by just 0.5 keV. Most modern high-resolution detectors are incapable of separating them and will record them as a single peak. When the separation is greater and two overlapping peaks can be distinguished, it may be possible to apportion the total area between the two. One method is to assign fractions based on the areas below the lower peak and above the upper peak. Figure D1. Comparing activity data (Bq kg −1 ) obtained from the photopeak at 186 keV. Each blue dot represents a sample in relation to the 1 : 1 black line. 226 Ra 186 activity is derived from Eq. (C1) corrected for the activity of 235 U directly inferred from its peak at 144 keV and plotted versus the weighted mean of activities derived from the 226 Ra daughters emitting 352 and 609 keV.
Although 238 U does not have significant gamma emissions itself, it will normally be in radioactive equilibrium with its short-lived daughter 234 Th (half-life 24.1 d). 234 Th has two significant photopeaks at around 63.3 and 92.6 keV, though both are in fact composite peaks. The 92.6 keV gamma line includes two separate 234 Th peaks at 92.4 and 92.8 keV, a small 235 U peak at 93.4 keV, two 231 Th peaks at 92.3 and 93.0 keV, and X-rays (Huy and Luyen, 2005). 226 Ra activity is best determined using one of the peaks of its short-lived daughter radionuclides 214 Pb or 214 Bi. 235 U may either be determined via its peaks at 144 and 205 keV or, in natural samples, from the established activity ratio with 238 U. See also De Corte et al. (2005) and https://www.nndc.bnl.gov/ (last access: 10 August 2021). In practice, the methods available for correcting the 186 keV emission deliver ambiguous results (e.g. Fig. C1). The correction procedure should therefore be replaced by the use of software in which a detailed detector characterization is combined with the MCNP Monte Carlo modelling code (see Appendix C for software references).
Author contributions. BM designed the work and wrote the text; PGA and PJN oversaw text and figures dealing with gamma radiation, instruments, and associated equations.
Competing interests. The contact author has declared that neither they nor their co-authors have any competing interests.
Disclaimer. Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.