1. bookAHEAD OF PRINT
Journal Details
License
Format
Journal
First Published
04 Jul 2007
Publication timeframe
1 time per year
Languages
English
access type Open Access

Component resolved equivalent dose estimation using TL glow curves of travertine samples from Anatolia, Turkey

Published Online: 29 Apr 2020
Page range: -
Received: 19 Jan 2019
Accepted: 12 Jul 2019
Journal Details
License
Format
Journal
First Published
04 Jul 2007
Publication timeframe
1 time per year
Languages
English
Abstract

This study provides methodological aspects on the equivalent dose estimation for travertine samples, namely heated calcium carbonate, using the TL multiple-aliquot additive-dose approach. Large equivalent doses (EDs), within the range 750–1300 Gy, were calculated using the plateau method based on the NTL glow curve. Moreover, a component resolved TL glow-peak analysis was carried out, using the integrated intensity of the NTL glow peaks for the ED estimation after deconvolution based on the OTOR model. Three different TL peaks were used, termed P4, P5 and P6. The integrated intensity of TL glow peak P4 resulted in age overestimation (15–26%), compared to the age provided using the plateau method. This overestimation could be attributed to the fact that the temperature range of P4 does not coincide with the plateau region of each sample. Milder overestimation (8–14%) was noticed using the integrated intensity of TL glow peak P6, mostly due to the poor deconvolution resolution. Only the integrated intensity of TL glow peak P5 after deconvolution provides ED values compatible with those yielded using the plateau method, with good accuracy. The present study suggests not using the TL intensity (neither in terms of integrated intensity nor of peak height intensity) for ED estimation; instead it is highly recommended to use either the plateau method, or alternatively integrated intensity of TL peak P5 after deconvolution. Unfortunately, using the peak height of TL P5 is not recommended, due to overlapping with P6.

Keywords

Introduction

Calcium carbonate, with the chemical formula CaCO3, stands among the most widespread minerals of earth’s crust. The importance of dating calcium carbonate based formations for setting age limits has been well demonstrated by both the uranium series disequilibrium method as well as Electron Paramagnetic Resonance (EPR) techniques. As far as luminescence techniques are concerned, Wintle (1977, 1978) has undertaken some of the initial work to assess whether the TL signal from calcite could be used to date geologically relevant samples. The reader may also refer to McDougal (1968) for a review of earlier achievements on this related topic. Due to the problematic properties, or even the lack, of optically stimulated luminescence (OSL) signal, thermoluminescence (TL) stands as the only luminescence technique available for dating applications for these formations.

Despite the numerous reports dealing with the EPR dating applications of calcium carbonates (for example Ikeya, 1993 and references therein), a limited number of papers deal with the corresponding luminescence applications of this mineral within the last 35 years (Debenham, 1983; Debenham and Aitken, 1984; Down et al., 1985; Ninagawa, 1987; Ninagawa et al., 1988, 1992; Liritzis, 1989; 2010; 2011; Duller et al., 2009; Stirling et al., 2014; Polymeris et al., 2016), both successful and unsuccessful. Throughout this aforementioned literature, the TL properties of various types of CaCO3 samples have been reported and the ages were in most cases lower than 500 ka. The present study provides methodological aspects on the ED estimation using the multiple-aliquot additive-dose procedure in the TL of calcium carbonate samples. Two different alternative approaches were adopted, namely a TL peak resolved analysis following deconvolution and the plateau method without deconvolution. Four different samples from the Quaternary were studied.

Materials and methods
Apparatus

All luminescence measurements were conducted using a commercial Risø TL/OSL system (model TL/OSL-DA-20), equipped with blue LEDs (470 nm, FWHM 20 nm) and a 90Sr/90Y beta particle source, delivering a nominal dose rate of 0.1083 Gy/s. An EMI 9235QB photomultiplier tube was used for the detection of all luminescence signals. For the cases of the few OSL and TA – OSL measurements, a Hoya (U-340) filter was used for light detection (340 nm, FWHM 80nm) while for all TL measurements only a BG-39 glass filter was used to restrict the wavelengths detected, similar to the studies of Duller et al., (2009) as well as of Polymeris et al. (2016). This filter has a transmission window that ranges from 350 to 600 nm and it was used for all TL measurements. All heatings and TL measurements were performed in a nitrogen atmosphere using a low constant heating rate of 2°C/s in order to avoid significant temperature lag (Kitis et al., 2015); the samples were heated up to the maximum temperature of 500°C. All luminescence measurements were conducted at the Institute of Nuclear Sciences of Ankara University.

For the structural characterization of the materials, the X-ray powder diffraction (XRPD) technique was applied; the minerals were analyzed to verify their crystalline structure using a Philips diffractometer with monochromatic CuKa (λ = 0.15406 nm) radiation. For a description of the elemental-chemical composition of the samples, Polarized Energy Dispersive X-ray fluorescence (XRF) analysis was performed using an X-ray Fluorescence spectrometer Spectro XLAB-2000. Both PXRD and XRF measurements and analysis were performed at the Earth Sciences Application and Research Center of Ankara University (YEBIM).

Materials and treatment

In this work, four different travertine samples were collected from the Bor district of Niğde city, Central Anatolia, Turkey. Travertine deposits are the product of calcium carbonate precipitation under cool water (Ford and Pedley, 1996). The thickness of the deposits differs between 15 and 25 meters. These deposits are considered belonging to the Quaternary, based on indirect proxies as well as preliminary luminescence dating results (Öztürk et al., 2018). The Southern part of study area was bordered with Bor Segment of Tuzgölü (Saltlake) Fault Zone, Turkey.

Based on Electron dispersion X-ray (EDX) spectroscopy, all samples contain predominantly Ca (38.85%), O (37.87%), while to a lesser extent C (21.03%); all percentage values indicate average values over all samples of the present study. This observation is in excellent agreement with the results of X-ray powder diffraction (XRPD) patterns, indicating the presence of the phase of CaCO3 solely. The sample codes as well as the concentration of each main element are presented in Table 1 for each sample. Typically, travertines are formed where cool carbonate supersaturated freshwater degasses and this is where karst water reaches the surface. Thus, various impurities are expected into the otherwise clean CaCO3 precipitate. Table 1 also presents the XRF results for other available elements as well, including K, Na, Fe, Mg, Mn, Si, Al, Ti and Cl. In some cases, the corresponding concentrations are below the detection limit of the PXRD technique.

Content of main elements. The detection limit for the (%) concentration is 0.01%. BDL=Below Detection Limits.

A/ASampleCodesCacontent(%)Ccontent(%)O2content(%)Kcontent(%)Nacontent(%)Fecontent(%)Mgcontent(%)Mncontent(%)Sicontent(%)Alcontent(%)Ticontent(%)Clcontent(%)
1SB140.1537.4420.53BDL0.050.510.260.030.130.760.09BDL
2SB237.4138.8421.23BDL0.070.410.18BDL0.320.710.710.06
3SB339.0538.1520.28BDLBDL0.610.770.070.210.330.510.03
4SB438.7937.0522.07BDL0.030.530.580.080.140.190.480.05

Before luminescence measurements, treatments and preparation were under-taken in subdued red filtered light conditions. An almost 1 cm thick, outer layer was removed from each sample in the laboratory, to eliminate the light-subjected parts of the samples. Finally, grains with dimensions in the range 4–12 μm were extracted by gently crushing in an agate pestle and mortar, then suspended in acetone and finally precipitated onto 1 cm diameter aluminum discs (Fleming, 1979). As PXRD indicated the presence of solely CaCO3, no further chemical procedure was applied.

Methods of analysis

For the equivalent dose estima tion the multiple aliquot additive dose procedure (MAAD) in TL was applied. A detailed description of the procedure can be found in Aitken (1985; 1998), Wagner (1998) and Liritzis et al. (2013). Mass reproducibility for all discs-aliquots was kept within ±3%. The additive doses applied were 0 (NTL), 150, 300, 450, 600 and 900 Gy. Four fresh aliquots were used for each additive dose. The equivalent dose (ED) was calculated using two different approaches. For the first approach, the ED was calculated using the plateau method (Aitken, 1985; Wagner, 1998) according to the following equation (Liritzis, 2010; Liritzis et al., 2015; Aidona et al., 2018):

ED=NTL(NTL+βi)NTLβi=NTLSi$$ED=\frac{NTL}{\left( NTL+{{\beta }_{i}} \right)-NTL}\cdot {{\beta }_{i}}=\frac{NTL}{{{S}_{i}}}$$

where NTL is the natural TL glow curves, ( NTL + βi) denotes the glow curves after the beta doses were added on top of the natural signal and βi the corresponding administered beta dose in units of Gy. Si denotes the slope of the (linear) fitting of the dose response curve and is defined according to the following equation (Polymeris et al., 2012; Liritzis et al., 2015):

Si=(NTL+βi)NTLβi$${{S}_{i}}=\frac{\left( NTL+{{\beta }_{i}} \right)-NTL}{{{\beta }_{i}}}$$

Eqs. 2.1 and 2.2 were applied to the entire NTL and NTL +βi glow curves. Since all TL curves were measured using 500 data points and heating rate of 1°C/s up to a maximum temperature of 500°C, the preceding calculations take place for each data point of the glow curves. Eventually Eq. 2.2 gives five different curves, each one corresponding to one among the additive doses applied with each data point corresponding to a different temperature in steps of 1°C. In the case of linearity, these five curves should coincide, at least in a selected glow curve temperature range. For this temperature range, the ED is independent on the selection of the integration limits of the TL signal; therefore, if one plots the ED versus the TL glow curve temperature range, a plateau region is yielded. So the ED is calculated within this temperature range, which is specified experimentally.

However, this latter approach does not provide TL peak resolved analysis, but just a plateau within a temperature range within the glow curve. In order to get TL peak-resolved analysis, a computerized glow curve deconvolution analysis was applied to all experimental TL glow curves, using the OTOR model (Halperin and Braner, 1960) and the analytical solution of the corresponding differential equations. ED was calculated using the integral of each TL peak after deconvolution analysis. The equations applied for the deconvolution analysis were the following (Sadek et al., 2014):

I(T)=Imexp(EkTTmTTm)W(zm)+W(zm)2W(z)+W(z)2$$I\left( T \right)={{I}_{m}}exp\left( -\frac{E}{kT}\cdot \frac{{{T}_{m}}-T}{{{T}_{m}}} \right)\cdot \frac{W\left( {{z}_{m}} \right)+W{{\left( {{z}_{m}} \right)}^{2}}}{W\left( z \right)+W{{\left( z \right)}^{2}}}$$

where W(z) is the Lambert function; z, being the independent variable of the equations, is approximated by the expression:

z=exp(R1Rln(1RR)+Eexp(EkTm)kTm2(11.05R1.26)F(T,E))$$z=\exp \left( \frac{R}{1-R}-\ln \left( \frac{1-R}{R} \right)+\frac{E\cdot \exp \left( \frac{E}{k{{T}_{m}}} \right)}{kT_{m}^{2}\left( 1-1.05\cdot {{R}^{1.26}} \right)}\cdot F\left( T,E \right) \right)$$

and zm is the value of z from the above equation, for T = Tm.

These Eqs. 2.2 and 2.3 are the analytical equations used to fit the experimental TL glow curves, with the two adjustable fitting parameters being the activation energy E and the ratio R (with R < 1) (Kitis and Vlachos, 2012). The value of the parameter R indicates the order of kinetics; R value close to 0 indicates first order of kinetics, while a value close to 1 indicates second order of kinetics (Şahiner, 2017). The function F(T,E) in this expression is the exponential integral appearing in TL models and is expressed in terms of the exponential integral function Ei[− E/kT] as (Chen and McKeever, 1997; Kitis et al., 2006):

F(T,E)=T0TeEkTdT=Texp(EkT)+EkEi[EkT]$$F\left( T,E \right)=\int_{{{T}_{0}}}^{T}{{{e}^{-\frac{E}{kT}}}}dT=T\cdot exp\left( -\frac{E}{kT} \right)+\frac{E}{k}\cdot {{E}_{i}}\left[ -\frac{E}{kT} \right]$$

The software package Microsoft Excel was used for the deconvolution analysis, using the Solver utility (Afouxenidis et al., 2012), while the goodness of fit was tested using the Figure Of Merit (FOM) of Balian and Eddy (1977) given by:

FOM(%)=100*p|TLexpTLfit|/pTLfit$$\text{FOM}\left( \text {%} \right)=100*{{\sum }_{\text{p}}}{\left| {\rm T{{L}}_{\exp }}-{\rm T{{L}}_{\text{fit}}} \right|}/{{{\sum }_{\text{p}}}}\;{\rm T{{L}}_{\text{fit}}}$$

The FOM index value provides a measure for the goodness of fit; the lowest its value, the best fit it is. Therefore, every fitting attempt results in minimizing the FOM index value, which was achieved by changing the set of the parameter values of each glow peak. FOM values in all cases were less than 0.9%, indicating the fitting quality.

Results and discussion
OSL and TA-OSL signals

Fig. 1 presents typical OSL as well as thermally assisted optically stimulated luminescence (TA – OSL, Polymeris, 2016) signals from one CaCO3 sample of the present study following natural irradiation, namely without any artificial irradiation in the laboratory. The OSL signal was measured at 110°C, while the TA – OSL at 200°C. It is obvious from this figure that both signals consist of one component with shape rather peculiar and extremely flat. In terms of intensity levels, while the OSL indicates a signal similar to background, the TA – OSL yields intensity as large as ~18 times the OSL background signal. As this figure reveals, a fast decaying component is missing from both signals. Neither the intensity nor the shapes of both OSL signals do change, even after illumination up to 1 ks using blue LEDs. This result indicates that these signals are extremely hard to bleach, at least within relative short time intervals. These features stand as indirect indications regarding the lack of significant quantities of either quartz or feldspar inclusions inside the clear CaCO3 precipitate. The lack of both quartz and feldspar inclusions could be further supported by the contents of Si and Na/K respectively, according to Table 1. In agreement with the previous reports by both Liritzis (1989; 2000; 2011) and Galloway (2002) we conclude that both OSL signals are not useful for luminescence dating. This is the reason why, for dating applications, the luminescence applications to calcium carbonates so far were limited to TL of mostly biogenic calcite (Duller et al., 2009; Liritzis, 2010; 2011; Polymeris et al., 2016). The TL signal of CaCO3 samples is known to be bleachable, however, prolonged bleaching times are required (hours or even days, McDougal, 1968; Liritzis, 1989; 2000). For a recent, indicative example, the authors could refer to Fig. 6 of Polymeris et al. (2016).

Fig. 1

Representative natural OSL and thermally assisted OSL (TA – OSL) signals corresponding to the CaCO3 sample with code SB3. These signals are not useful for dating applications.

ED determination using the plateau method

The natural TL glow curves of all samples exhibit the same main characteristics, namely a glow curve that has the form of a continuum with two distinct but overlapping TL peaks at 325 and at 450°C (Fig. 2a). Figs. 2b and 2c present the MAAD TL analysis results for the sample SB1, which are representative for all the samples that were subject to the present study.

Fig. 2

Typical results of the MAAD plateau method applied for sample SB1. Plot (a) presents a NTL glow curve (as an average of three independently measured glow curves), plot (b) presents the representative Si curves for a selection of additive doses according to equation (2) while plot (c) presents an ED plateau versus glow curve temperature according to equation (1). Error bars correspond to 1σ.

Fig. 2b presents the curves corresponding to the dose–response slopes (Si) which are plotted versus temperature for a selection of different additive doses applied for the sample SB1. Si curves were calculated according to Eq. 2.2 and the preceding subtraction takes place for each data point of the glow curves. According to Fig. 2b, these Si curves coincide, throughout a wide temperature region of the glow curve, namely between 250 and 400°C, confirming thus the dose response linearity within this aforementioned temperature range, despite the large additive doses. Consequently, the ED plateau is expected within the temperature range between 250°C and 400°C. Similar linearity Si features were also indicated for the rest of the samples.

Equivalent doses were calculated using Eq. 2.1. Since (a) all TL curves were measured using 500 data points and heating rate of 1°C/s up to a maximum temperature of 500°C and (b) the calculations take place for each data point of the glow curves, an ED value correspond to each TL glow curve temperature. If one plots these ED values versus the TL temperature, an ED plateau is expected to occur. A typical plot of ED plateau against glow curve temperature is presented in Fig. 2c. Errors, which are derived mainly from the uncertainties in curve fitting, are ±1σ and were calculated by standard error propagation analysis. In all cases, ED plateaus are wide enough, over 100°C wide. The equivalent doses were obtained as the mean values of the plateaus for each sample. Quite large ED values were obtained, ranging between 781(±71) and 1458(±166) Gy. Only linear fittings were performed to the dose response curves. Details of this ED calculation, such as the length, the value and the error of each ED plateau are presented for all samples in Table 2. It is worth emphasizing that for the estimation of the ED plateau region, all reheats (background TL glow curves) were subtracted. Moreover, for all samples the temperature range corresponding to the ED plateau lies within the temperature range of linearity according to the Si curves. Low dose supralinearity was not observed.

Data related to the ED estimation for all four samples.

A/ASampleCodesPlateau ΔΤ(°C)Plateau ED(Gy)P4 ED(Gy)P5 ED(Gy)P6 ED(Gy)
1SB1270–3701258 ± 1551578 ± 1711312 ± 1371422 ± 158
2SB2280–3851458 ± 1661682 ± 1491489 ± 1551358 ± 185
3SB3260–3801077 ± 911236 ± 1331113 ± 1071193 ± 128
4SB4270–380781 ± 70932 ± 98812 ± 86850 ± 88
Deconvolution results and component resolved ED estimation

Fig. 3 presents a NTL + β (450 Gy) glow curve, deconvolved into its individual TL peaks according to Eqs. 2.32.5. For all NTL + βi glow curves, six individual TL glow peaks were used in order to achieve fit of excellent quality (FOM <0.9%). Nevertheless, for the case of the NTL glow curves, only three TL glow peaks are required (namely P4, P5 and P6) in order to achieve fit of excellent quality. A TL background curve has also been included within the deconvolution procedure, as Fig. 3 indicates. Table 3 presents an outline of the deconvolution analysis results, including the values of Tmax, activation energy E and the parameter R which gives information on the order of kinetics. It should be emphasized that the results obtained for all four samples studied were similar. The reproducibility of the Tmax values with ranging additive doses stands in excellent agreement with the R values of <0.1. Moreover, based on the deconvolution results, the lifetimes of the TL traps included in the NTL glow curves are adequate enough for providing ages of the order of 1 Ma.

Fig. 3

Representative deconvolution analysis of a TL glow curve corresponding to the sample SB1, after a dose of 450 Gy (open data points). Six different TL glow peaks and a black body radiation curve (BGK) were used in order to fit the experimental glow curve. Each TL glow peak is presented in solid line and is denoted as Pi. Orange solid line indicates the fitting curve. Table 2 indicates the activation energies of the TL peaks.

Deconvolution results for fitting parameters of all four samples.

A/ASampleCodesP4 Tmax(°C)P4 RP4 E(eV)P5 Tmax(°C)P5 RP5 E(eV)P6 Tmax(°C)P6 RP6 E(eV)
1SB1328 ± 40.061.34 ± 0.12399 ± 110.091.59 ± 0.18453 ± 220.011.81 ± 0.21
2SB2326 ± 50.041.37 ± 0.11402 ± 100.051.63 ± 0.17457 ± 210.011.86 ± 0.19
3SB3328 ± 30.091.34 ± 0.13404 ± 130.071.64 ± 0.16463 ± 220.021.88 ± 0.24
4SB4327 ± 50.051.41 ± 0.12402 ± 130.051.61 ± 0.17459 ± 240.011.79 ± 0.20

Besides the ED plateau method, equivalent doses were also calculated using the integrated TL intensity of each one of the three TL glow peaks which were used to deconvolve the NTL glow curves. An outline of the TL peak resolved ED estimation results could be found in Table 2. TL glow peak P4 is quite prominent, indicating very intense signal. The fitting parameters corresponding to this specific TL peak yield the lowest error values. However, the equivalent dose values calculated using the deconvolved integrated specific TL peak are overestimated when compared to the ED plateau value. Fig. 4 presents a representative additive dose response curve corresponding to sample SB1, plotted for the case of the integrated intensity of TL glow peak P4 after deconvolution.

Fig. 4

Representative additive dose growth curve in terms of the integrated intensity of TL glow peak P4 after deconvolution, for the sample with code SB1. The arrow indicates the ED value.

Despite the fact that the additive doses extend up to 900 Gy, linearity of the TL signal is characteristic for all samples. ED is overestimated at 1578(±177) Gy, while the corresponding plateau ED was calculated as 1258(±155) Gy. Similar overestimation was also yielded while using solely the integral of P6 TL peak. Component resolved overestimation ranges between 8 and 26%, depending on the TL peak used. Nevertheless, the most prominent TL glow peak P4 yields the higher overestimation percentage (15–26%). This experimental overestimation could be attributed to the fact that temperature range of TL glow peak P4 does not coincide with the plateau region for all samples. The overestimation yielded using TL glow peak P6, ranging between 8–14%, could be attributed to the fact that the TL glow curve should be measured up to temperatures higher than 500°C. As this is not the case, the resolution of the deconvolution technique is not the optimum. On the contrary, the additive dose procedure using the TL glow peak P5 yields the most comparable ED results with the ED plateau methodology. This similarity is quite prominent from the results of Table 1 and it is attributed to the fact that TL peak P5 lies within the corresponding plateau temperature range. Once again, it is worth noting that for all four samples and all three TL peaks, the linearity is quite impressive, despite the high additive doses applied.

TL ages

Estimation of the uranium, thorium and potassium geochemical content of the samples in the present study is currently pending. According to previous related luminescence studies using CaCO3, this geochemical content is quite low, indicating annual dose rate values lower than 1Gy/ka (for example Liritzis, 2010; 2011; Polymeris et al., 2016). Based on this rationale, the TL ages of these samples could reach values at least of 1 Ma, or even higher, namely in the Quaternary era. Nevertheless, as the present manuscript aims at presenting methodological aspects, as it was the case of the report by Duller et al., (2009), the detailed presentation of an integrated study dealing with the ages and their geological implications will be presented elsewhere.

Conclusions

The present study provides methodological aspects on the ED estimation for travertine samples, namely calcium carbonate samples, using the multiple-aliquot additive-dose procedure in TL.

Heated CaCO3 can be effectively used for calculating equivalent doses within the range 750–1300 Gy.

OSL and TA – OSL curves were measured and proven inappropriate for age assessment.

The integrated intensity of TL glow peak P4 after deconvolution, provides age overestimation, compared to the age provided using the plateau method. This overestimation could be attributed to the fact that temperature range of P4 does not coincide with the plateau region of each sample.

Using the integrated intensity of TL glow peak P6 after deconvolution is not suggested, as the resolution for resolving this peak is not good. As a result of this poor resolution, ED overestimation is also monitored.

The integrated intensity of TL glow peak P5 after deconvolution provides ED values compatible with those yielded using the plateau methodology within the entire TL glow curve, with better accuracy, namely lower error values.

The present study suggests not using the TL intensity (neither in terms of integrated intensity nor of peak height intensity) for ED estimation; instead it is highly recommended to use either the plateau methodology, or alternatively integrated intensity of TL peak P5 after deconvolution. Unfortunately, using the peak height of TL P5 is not recommended, due to overlapping with P6.

TL age results indicate a formation during the Quaternary.

Fig. 1

Representative natural OSL and thermally assisted OSL (TA – OSL) signals corresponding to the CaCO3 sample with code SB3. These signals are not useful for dating applications.
Representative natural OSL and thermally assisted OSL (TA – OSL) signals corresponding to the CaCO3 sample with code SB3. These signals are not useful for dating applications.

Fig. 2

Typical results of the MAAD plateau method applied for sample SB1. Plot (a) presents a NTL glow curve (as an average of three independently measured glow curves), plot (b) presents the representative Si curves for a selection of additive doses according to equation (2) while plot (c) presents an ED plateau versus glow curve temperature according to equation (1). Error bars correspond to 1σ.
Typical results of the MAAD plateau method applied for sample SB1. Plot (a) presents a NTL glow curve (as an average of three independently measured glow curves), plot (b) presents the representative Si curves for a selection of additive doses according to equation (2) while plot (c) presents an ED plateau versus glow curve temperature according to equation (1). Error bars correspond to 1σ.

Fig. 3

Representative deconvolution analysis of a TL glow curve corresponding to the sample SB1, after a dose of 450 Gy (open data points). Six different TL glow peaks and a black body radiation curve (BGK) were used in order to fit the experimental glow curve. Each TL glow peak is presented in solid line and is denoted as Pi. Orange solid line indicates the fitting curve. Table 2 indicates the activation energies of the TL peaks.
Representative deconvolution analysis of a TL glow curve corresponding to the sample SB1, after a dose of 450 Gy (open data points). Six different TL glow peaks and a black body radiation curve (BGK) were used in order to fit the experimental glow curve. Each TL glow peak is presented in solid line and is denoted as Pi. Orange solid line indicates the fitting curve. Table 2 indicates the activation energies of the TL peaks.

Fig. 4

Representative additive dose growth curve in terms of the integrated intensity of TL glow peak P4 after deconvolution, for the sample with code SB1. The arrow indicates the ED value.
Representative additive dose growth curve in terms of the integrated intensity of TL glow peak P4 after deconvolution, for the sample with code SB1. The arrow indicates the ED value.

Data related to the ED estimation for all four samples.

A/ASampleCodesPlateau ΔΤ(°C)Plateau ED(Gy)P4 ED(Gy)P5 ED(Gy)P6 ED(Gy)
1SB1270–3701258 ± 1551578 ± 1711312 ± 1371422 ± 158
2SB2280–3851458 ± 1661682 ± 1491489 ± 1551358 ± 185
3SB3260–3801077 ± 911236 ± 1331113 ± 1071193 ± 128
4SB4270–380781 ± 70932 ± 98812 ± 86850 ± 88

Deconvolution results for fitting parameters of all four samples.

A/ASampleCodesP4 Tmax(°C)P4 RP4 E(eV)P5 Tmax(°C)P5 RP5 E(eV)P6 Tmax(°C)P6 RP6 E(eV)
1SB1328 ± 40.061.34 ± 0.12399 ± 110.091.59 ± 0.18453 ± 220.011.81 ± 0.21
2SB2326 ± 50.041.37 ± 0.11402 ± 100.051.63 ± 0.17457 ± 210.011.86 ± 0.19
3SB3328 ± 30.091.34 ± 0.13404 ± 130.071.64 ± 0.16463 ± 220.021.88 ± 0.24
4SB4327 ± 50.051.41 ± 0.12402 ± 130.051.61 ± 0.17459 ± 240.011.79 ± 0.20

Content of main elements. The detection limit for the (%) concentration is 0.01%. BDL=Below Detection Limits.

A/ASampleCodesCacontent(%)Ccontent(%)O2content(%)Kcontent(%)Nacontent(%)Fecontent(%)Mgcontent(%)Mncontent(%)Sicontent(%)Alcontent(%)Ticontent(%)Clcontent(%)
1SB140.1537.4420.53BDL0.050.510.260.030.130.760.09BDL
2SB237.4138.8421.23BDL0.070.410.18BDL0.320.710.710.06
3SB339.0538.1520.28BDLBDL0.610.770.070.210.330.510.03
4SB438.7937.0522.07BDL0.030.530.580.080.140.190.480.05

Afouxenidis D, Polymeris GS, Tsirliganis NC and Kitis G, 2012. Computerised curve deconvolution of TL/OSL curves using a popular spreadsheet program. Radiation Protection Dosimetry 149(3): 363–370, DOI 10.1093/rpd/ncr315.AfouxenidisDPolymerisGSTsirliganisNCKitisG2012Computerised curve deconvolution of TL/OSL curves using a popular spreadsheet programRadiation Protection Dosimetry149336337010.1093/rpd/ncr315Open DOISearch in Google Scholar

Aidona E, Polymeris GS, Camps P, Kondopoulou D, Ioannidis N and Raptis K, 2018. Archaeomagnetic versus luminescence methods: the case of an Early Byzantine ceramic workshop in Thessaloniki, Greece. Archaeological and Anthropological Sciences 10: 725– 741, DOI 10.1007/s12520-017-0494-5.AidonaEPolymerisGSCampsPKondopoulouDIoannidisNRaptisK2018Archaeomagnetic versus luminescence methods: the case of an Early Byzantine ceramic workshop in Thessaloniki, GreeceArchaeological and Anthropological Sciences1072574110.1007/s12520-017-0494-5Open DOISearch in Google Scholar

Aitken MJ, 1985. Thermoluminescence dating. Academic Press, London.AitkenMJ1985Thermoluminescence datingAcademic PressLondonSearch in Google Scholar

Aitken MJ, 1998. An introduction to optical dating: the dating of quaternary sediments by the use of photon-stimulated luminescence. Oxford University Press.AitkenMJ1998An introduction to optical dating: the dating of quaternary sediments by the use of photon-stimulated luminescenceOxford University PressSearch in Google Scholar

Balian HG and Eddy NW, 1977. Figure Of Merit (FOM): an improved criterion over the normalized chi-square test for assign the goodness-of-fit of gamma ray spectral peaks. Nuclear Instruments and Methods 145, 389–395, DOI 10.1016/0029-554X(77)90437-2.BalianHGEddyNW1977Figure Of Merit (FOM): an improved criterion over the normalized chi-square test for assign the goodness-of-fit of gamma ray spectral peaksNuclear Instruments and Methods14538939510.1016/0029-554X(77)90437-2Open DOISearch in Google Scholar

Chen R and McKeever SWS, 1997. Theory of Thermoluminescence and Related Phenomena. World Scientific.ChenRMcKeeverSWS1997Theory of Thermoluminescence and Related PhenomenaWorld ScientificSearch in Google Scholar

Debenham NC, 1983. Reliability of thermos-luminescence dating of stalagmitic calcite. Nature 304: 154-156, DOI 10.1038/304154a0.DebenhamNC1983Reliability of thermos-luminescence dating of stalagmitic calciteNature30415415610.1038/304154a0Open DOISearch in Google Scholar

Debenham NC and Aitken MJ, 1984. Thermoluminescence dating of stalagmitic calcite. Archaeometry 26: 155–170, DOI 10.1111/j.1475-4754.1984.tb00330.x.DebenhamNCAitkenMJ1984Thermoluminescence dating of stalagmitic calciteArchaeometry2615517010.1111/j.1475-4754.1984.tb00330.xOpen DOISearch in Google Scholar

Down JS, Flower R, Strain JA and Townsend PD, 1985. Thermoluminescence emission spectra of calcite and Iceland spar. Nuclear Tracks and Radiation Measurements 10: 581–589, DOI 10.1016/0735-245X(85)90061-4.DownJSFlowerRStrainJATownsendPD1985Thermoluminescence emission spectra of calcite and Iceland sparNuclear Tracks and Radiation Measurements1058158910.1016/0735-245X(85)90061-4Open DOISearch in Google Scholar

Duller GAT, Penkman KEH and Wintle AG, 2009. Assessing the potential for using biogenic calcites as dosemeters for luminescence dating. Radiation Measurements 44: 429–433, DOI 10.1016/j.radmeas.2009.02.008.DullerGATPenkmanKEHWintleAG2009Assessing the potential for using biogenic calcites as dosemeters for luminescence datingRadiation Measurements4442943310.1016/j.radmeas.2009.02.008Open DOISearch in Google Scholar

Fleming S, 1979. Thermoluminescence techniques in archaeology. Clarendon Press.FlemingS1979Thermoluminescence techniques in archaeologyClarendon PressSearch in Google Scholar

Ford TD and Pedley HM, 1996. A review of tufa and travertine deposits of the world. Earth-Science Reviews 41(3–4): 117–175, DOI 10.1016/S0012-8252(96)00030-X.FordTDPedleyHM1996A review of tufa and travertine deposits of the worldEarth-Science Reviews413–411717510.1016/S0012-8252(96)00030-XOpen DOISearch in Google Scholar

Galloway RB, 2002. Does limestone show useful optically stimulated luminescence? Ancient TL 20: 1–7.GallowayRB2002Does limestone show useful optically stimulated luminescence?Ancient TL2017Search in Google Scholar

Halperin A and Braner AA, 1960. Evaluation of thermal activation energies from glow curves. Physical Review 117: 408–415, DOI 10.1103/PhysRev.117.408.HalperinABranerAA1960Evaluation of thermal activation energies from glow curvesPhysical Review11740841510.1103/PhysRev.117.408Open DOISearch in Google Scholar

Ikeya M, 1993. New Applications of Electron Spin Resonance: Dating, Dosimetry and Microscopy. World Scientific, Singapore.IkeyaM1993New Applications of Electron Spin Resonance: Dating, Dosimetry and MicroscopyWorld ScientificSingaporeSearch in Google Scholar

Kitis G, Chen R, Pagonis V, Carinou E and Kamenopoulou V, 2006. Thermoluminescence under an exponential heating function: I. Theory. Journal Of Physics D: Applied Physics 39: 1500–1507, DOI 10.1088/0022-3727/39/8/008.KitisGChenRPagonisVCarinouEKamenopoulouV2006Thermoluminescence under an exponential heating function: I. TheoryJournal Of Physics D: Applied Physics391500150710.1088/0022-3727/39/8/008Open DOISearch in Google Scholar

Kitis G and Vlachos ND, 2012. General semi-analytical expressions for TL, OSL and other luminescence stimulation modes derived from the OTOR model using the Lambert W-function. Radiation Measurements 48: 47–54, DOI 10.1016/j.radmeas.2012.09.006.KitisGVlachosND2012General semi-analytical expressions for TL, OSL and other luminescence stimulation modes derived from the OTOR model using the Lambert W-functionRadiation Measurements48475410.1016/j.radmeas.2012.09.006Open DOISearch in Google Scholar

Kitis G, Kiyak NG and Polymeris GS, 2015. Temperature lags of luminescence measurements in a commercial luminescence reader. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms 359: 60–63, DOI 10.1016/j.nimb.2015.07.041.KitisGKiyakNGPolymerisGS2015Temperature lags of luminescence measurements in a commercial luminescence readerNuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms359606310.1016/j.nimb.2015.07.041Open DOISearch in Google Scholar

Liritzis I, 2000. Advances in thermo- and opto- luminescence dating of environmental materials (sedimentary de-posits): part II: applications. Global Nest 2(1): 29–49.LiritzisI2000Advances in thermo- and opto- luminescence dating of environmental materials (sedimentary de-posits): part II: applicationsGlobal Nest212949Search in Google Scholar

Liritzis I, 2010. Strofilas (Andros Island, Greece): new evidence for the cycladic final neolithic period through novel dating methods using luminescence and obsidian hydration. Journal of Archaeological Science 37: 1367–1377, DOI 10.1016/j.jas.2009.12.041.LiritzisI2010Strofilas (Andros Island, Greece): new evidence for the cycladic final neolithic period through novel dating methods using luminescence and obsidian hydrationJournal of Archaeological Science371367137710.1016/j.jas.2009.12.041Open DOISearch in Google Scholar

Liritzis I, 2011. Surface dating by luminescence: an overview. Geochronometria 38: 292–302, DOI 10.2478/s13386-011-0032-7.LiritzisI2011Surface dating by luminescence: an overviewGeochronometria3829230210.2478/s13386-011-0032-7Open DOISearch in Google Scholar

Liritzis Y, 1989. Dating of calcites: some aspects of radiation survey in caves and dose-rates. Bulletin Geologique des Pays Helleniques 34(1): 123–136.LiritzisY1989Dating of calcites: some aspects of radiation survey in caves and dose-ratesBulletin Geologique des Pays Helleniques341123136Search in Google Scholar

Liritzis I, Singhvi AK, Feathers JK, Wagner GA, Kadereit A, Zacharias N and Li SH, 2013. Luminescence Dating in Archaeology,Anthropology and Geoarchaeology: an Overview. Springer Briefs inEarth System Sciences. Springer, Heidelberg. DOI 10.1007/978-3-319-00170-8.LiritzisISinghviAKFeathersJKWagnerGAKadereitAZachariasNLiSH2013Luminescence Dating in Archaeology,Anthropology and Geoarchaeology: an OverviewSpringer Briefs inEarth System Sciences. Springer, Heidelberg10.1007/978-3-319-00170-8Open DOISearch in Google Scholar

Liritzis I, Aravantinos V, Polymeris GS, Zacharias N, Fappas I, Agiamarniotis G, Sfampa IK, Vafiadou A and Kitis G, 2015. Witnessing prehistoric Delphi by luminescence dating. Comptes Rendus Palevol 14: 219–232, DOI 10.1016/j.crpv.2014.12.007.LiritzisIAravantinosVPolymerisGSZachariasNFappasIAgiamarniotisGSfampaIKVafiadouAKitisG2015Witnessing prehistoric Delphi by luminescence datingComptes Rendus Palevol1421923210.1016/j.crpv.2014.12.007Open DOISearch in Google Scholar

McDougal DJ, 1968. Thermoluminescence of Geological Materials. London University Press.McDougalDJ1968Thermoluminescence of Geological MaterialsLondon University PressSearch in Google Scholar

Ninagawa K, 1987. Thermoluminescence dating of fossil calcite shells. Japanese Journal of Applied Physics 26: 2127-2133, DOI 10.1143/JJAP.26.2127.NinagawaK1987Thermoluminescence dating of fossil calcite shellsJapanese Journal of Applied Physics262127213310.1143/JJAP.26.2127Open DOISearch in Google Scholar

Ninagawa K, Adachi K, Uchimura N, Yamamoto I, Wada T, Yamashita Y, Takashima I, Sekimoto K and Hasegawa H, 1992. Thermoluminescence dating of calcite shells in the Pectinidae family. Quaternary Science Reviews 11: 121–126, DOI 10.1016/0277-3791(92)90052-A.NinagawaKAdachiKUchimuraNYamamotoIWadaTYamashitaYTakashimaISekimotoKHasegawaH1992Thermoluminescence dating of calcite shells in the Pectinidae familyQuaternary Science Reviews1112112610.1016/0277-3791(92)90052-AOpen DOISearch in Google Scholar

Ninagawa K, Takahashi N, Wada T, Yamamoto I, Yamashita N and Yamashita Y, 1988. Thermoluminescence measurements of a calcite shell for dating. Quaternary Science Reviews 7: 367–371, DOI 10.1016/0277-3791(88)90031-5.NinagawaKTakahashiNWadaTYamamotoIYamashitaNYamashitaY1988Thermoluminescence measurements of a calcite shell for datingQuaternary Science Reviews736737110.1016/0277-3791(88)90031-5Open DOISearch in Google Scholar

Öztürk MZ, Şener MF and Şahiner E, 2018. Quaternary slip-rates of the Bor segment of Tuzgölü fault zone. Omer Halisdemir University Journal of Engineering Sciences 7(3): 1049–1053 (In Turkish).ÖztürkMZŞenerMFŞahinerE2018Quaternary slip-rates of the Bor segment of Tuzgölü fault zoneOmer Halisdemir University Journal of Engineering Sciences7310491053In TurkishSearch in Google Scholar

Polymeris GS, 2016. Thermally assisted OSL (TA-OSL) from various luminescence phosphors; an overview. Radiation Measurements 90: 145–152, DOI 10.1016/j.radmeas.2016.01.035.PolymerisGS2016Thermally assisted OSL (TA-OSL) from various luminescence phosphors; an overviewRadiation Measurements9014515210.1016/j.radmeas.2016.01.035Open DOISearch in Google Scholar

Polymeris GS, Erginal AE and Kiyak NG, 2012. A comparative morphological, compositional and TL study of Tenedos (Bozcaada) and Şile aeolianites, Turkey. Mediterranean Archaeology & Archaeometry 12(2): 117–131.PolymerisGSErginalAEKiyakNG2012A comparative morphological, compositional and TL study of Tenedos (Bozcaada) and Şile aeolianites, TurkeyMediterranean Archaeology & Archaeometry122117131Search in Google Scholar

Polymeris GS, Kitis G, Kiyak NG, Theodosoglou E, Tsirliganis NC, Ertek A and Erginal AE, 2016. Dating fossil root cast (Black Sea coast, Turkey) using thermoluminescence: Implications for wind-blown drift of shelf carbonates during MIS 2. Quaternary International 401: 184–193, DOI 10.1016/j.quaint.2015.05.060.PolymerisGSKitisGKiyakNGTheodosoglouETsirliganisNCErtekAErginalAE2016Dating fossil root cast (Black Sea coast, Turkey) using thermoluminescence: Implications for wind-blown drift of shelf carbonates during MIS 2Quaternary International40118419310.1016/j.quaint.2015.05.060Open DOISearch in Google Scholar

Sadek AM, Eissa HM, Basha AM and Kitis G, 2014. Development of the peak fitting and peak shape methods to analyze the thermoluminescence glow-curves generated with exponential heating function. Nuclear Instruments & Methods in Physics Research Section B-Beam Interactions With Materials And Atoms 330: 103–107, DOI 10.1016/j.nimb.2014.04.006.SadekAMEissaHMBashaAMKitisG2014Development of the peak fitting and peak shape methods to analyze the thermoluminescence glow-curves generated with exponential heating functionNuclear Instruments & Methods in Physics Research Section B-Beam Interactions With Materials And Atoms33010310710.1016/j.nimb.2014.04.006Open DOISearch in Google Scholar

Şahiner E, 2017. Deconvolution analysis of thermoluminescent glow curves in various commercial dosimeters using two different approaches in the framework of the one-trap, one-recombination model. Turkish Journal of Physics 41: 477–490, DOI 10.3906/fiz-1704-25.ŞahinerE2017Deconvolution analysis of thermoluminescent glow curves in various commercial dosimeters using two different approaches in the framework of the one-trap, one-recombination modelTurkish Journal of Physics4147749010.3906/fiz-1704-25Open DOISearch in Google Scholar

Stirling RJ, Duller GAT and Roberts HM, 2014. Developing a single-aliquot protocol for measuring equivalent dose in biogenic carbonates. Radiation Measurements 47: 725–731, DOI 10.1016/j.radmeas.2012.01.010.StirlingRJDullerGATRobertsHM2014Developing a single-aliquot protocol for measuring equivalent dose in biogenic carbonatesRadiation Measurements4772573110.1016/j.radmeas.2012.01.010Open DOISearch in Google Scholar

Wagner GA, 1998. Age determination of young rocks and artifacts: physical and chemical clocks in quaternary geology and archaeology. Springer-Verlag, Berlin-Heidelberg.WagnerGA1998Age determination of young rocks and artifacts: physical and chemical clocks in quaternary geology and archaeologySpringer-VerlagBerlin-HeidelbergSearch in Google Scholar

Wintle AG, 1977. Thermoluminescence dating of minerals – traps for the unwary. Journal of Electrostatics 3: 281–288, DOI 10.1016/0304-3886(77)90100-0.WintleAG1977Thermoluminescence dating of minerals – traps for the unwaryJournal of Electrostatics328128810.1016/0304-3886(77)90100-0Open DOISearch in Google Scholar

Wintle AG, 1978. A thermoluminescence dating study of some Quaternary calcite: potential and problems. Canadian Journal of Earth Sciences 15: 1977–1986, DOI 10.1139/e78-208.WintleAG1978A thermoluminescence dating study of some Quaternary calcite: potential and problemsCanadian Journal of Earth Sciences151977198610.1139/e78-208Open DOISearch in Google Scholar

Recommended articles from Trend MD

Plan your remote conference with Sciendo