Step 2: Algebraic decomposition of all CW-OSL curves

This report was automatically created at 2021-02-16 09:09:20 by the function RLum.OSL_decomposition() of the R package OSLdecomposition (21.2.15). Visit luminescence.de to learn more.

---

Results

Step 2 of a component-resolved dose calculation is the separation of the signal components in each individual measurement. Mathematically this process is called decomposition. The number of components \(K\) and their signal decay rates were determined in Step 1:

• \(K=3\)
• \(\lambda_1=1.6514\) \(s^{-1}\) (Fast)
• \(\lambda_2=0.5231\) \(s^{-1}\) (Medium)
• \(\lambda_3=0.0094\) \(s^{-1}\) (Slow2)

These parameters lead to the following CW-OSL model:

\[I(t) = n_1e^{-1.6514t}+n_2e^{-0.5231t}+n_3e^{-0.0094t} \]

The signal component intensities \(n_1\) to \(n_K\) were calculated for each of the 384 CW-OSL measurements. First, 3 signal bin intervals were defined. The value of each bin was set equal to the formula above, integrated over the bin interval. In consequence, each measurement was described by 3 equations with 3 unknown intensity values \(n\). These equations systems were solved by an algebraic approach, see section Methods for details.

First two records

The plots in figure 1 show the component separation results for the first two measurements of each aliquot.

#1

Figure 2.1: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#2

Figure 2.2: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#3

Figure 2.3: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#4

Figure 2.4: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#5

Figure 2.5: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#6

Figure 2.6: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#7

Figure 2.7: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#8

Figure 2.8: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#9

Figure 2.9: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#10

Figure 2.10: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#11

Figure 2.11: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#12

Figure 2.12: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#13

Figure 2.13: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#14

Figure 2.14: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#15

Figure 2.15: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#16

Figure 2.16: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#17

Figure 2.17: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#18

Figure 2.18: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#19

Figure 2.19: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#20

Figure 2.20: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#21

Figure 2.21: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#22

Figure 2.22: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#23

Figure 2.23: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

#24

Figure 2.24: Signal components and residual curve of first two OSL records.
Upper left: First CW-OSL measurement (index = 1), usually natural or recovered dose signal (\(L_n\)).
Upper right: Second CW-OSL measurement (index = 3), usually test-dose signal (\(T_n\)).
Middle left and right: Residual curves between decomposition model and measurement
Lower left and right: Estimated type of quartz OSL component (colored), decay rates \(\lambda\), signal intensity \(n\) and signal contribution to first channel \(I_0\).

Intensity over sequence

The behaviour of the signal component intensities in the course of measurement sequences contains information about irradiation and bleaching response and thermal- or photo-transfer effects of the associated defect states. Fig. 2 display the component intensity variations in the course of the measurement sequences. Just components which were bleached to at least 95 % during the analysed time window are displayed, to prevent scaling issues and data misinterpretation due to background correction issues.

#1

Figure 1.1: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#2

Figure 1.2: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#3

Figure 1.3: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#4

Figure 1.4: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#5

Figure 1.5: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#6

Figure 1.6: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#7

Figure 1.7: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#8

Figure 1.8: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#9

Figure 1.9: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#10

Figure 1.10: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#11

Figure 1.11: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#12

Figure 1.12: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#13

Figure 1.13: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#14

Figure 1.14: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#15

Figure 1.15: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#16

Figure 1.16: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#17

Figure 1.17: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#18

Figure 1.18: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#19

Figure 1.19: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#20

Figure 1.20: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#21

Figure 1.21: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#22

Figure 1.22: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#23

Figure 1.23: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

#24

Figure 1.24: OSL signal component intensities over SAR sequence steps.
Dashed line: \(T_x\) values; Dotted line: \(L_x\) values (dose-response curve only).
Just components which are bleached to at least 95 % at the end of the curve are displayed.

Intensity tables

The following tables summarize all in Step 2 calculated values for each RLum.Analysis object (= aliquot). The column record index refers to the index needed to plot a single measurement by calling: plot_OSLcurve(data_set[[aliquot_index]]@records[[data_index]])

#1

Table 1.1: CW-OSL signal intensities of all records for aliquot #1.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	187389 ± 875	42243 ± 800	166194 ± 1357	92.5 %	7.0 %	0.5 %
3	25	40517 ± 394	9175 ± 394	156914 ± 1216	90.9 %	6.9 %	2.2 %
4	50	81051 ± 607	16904 ± 585	169695 ± 1414	92.4 %	6.4 %	1.2 %
6	25	40932 ± 336	9371 ± 399	170573 ± 1436	90.7 %	7.0 %	2.3 %
7	85	137636 ± 677	27915 ± 754	221334 ± 1617	92.8 %	6.3 %	0.9 %
9	25	41808 ± 351	9592 ± 475	207082 ± 1476	90.3 %	6.9 %	2.8 %
10	155	244967 ± 1175	49848 ± 1173	322545 ± 2196	92.9 %	6.3 %	0.8 %
12	25	41734 ± 497	10351 ± 479	269803 ± 1798	89.0 %	7.4 %	3.6 %
13	190	305138 ± 1514	61258 ± 1611	415334 ± 2488	93.0 %	6.2 %	0.8 %
15	25	43466 ± 487	11504 ± 555	330728 ± 1849	88.0 %	7.8 %	4.1 %
16	50	89342 ± 561	22517 ± 557	316843 ± 2028	90.4 %	7.6 %	2.0 %
18	25	44156 ± 393	12448 ± 424	303196 ± 1700	88.0 %	8.3 %	3.7 %
19	0	792 ± 148	-464 ± 282	245172 ± 1476	36.8 %	-7.2 %	70.4 %
21	25	46725 ± 372	11525 ± 338	268216 ± 1612	89.4 %	7.4 %	3.2 %
23	50	92046 ± 718	22556 ± 608	277262 ± 1667	90.8 %	7.4 %	1.7 %
25	50	93946 ± 438	24686 ± 692	324176 ± 1882	90.1 %	7.9 %	1.9 %

#2

Table 1.2: CW-OSL signal intensities of all records for aliquot #2.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	183987 ± 666	37409 ± 744	168497 ± 1374	93.1 %	6.3 %	0.5 %
3	25	42383 ± 315	10561 ± 334	151080 ± 1230	90.5 %	7.6 %	2.0 %
4	50	84465 ± 459	21384 ± 558	164325 ± 1389	91.2 %	7.7 %	1.1 %
6	25	43062 ± 355	11457 ± 357	170856 ± 1524	89.8 %	8.0 %	2.2 %
7	85	145608 ± 705	36404 ± 543	231252 ± 1499	91.4 %	7.6 %	0.9 %
9	25	43983 ± 368	13415 ± 434	220878 ± 1536	88.2 %	9.0 %	2.7 %
10	155	268132 ± 1601	63429 ± 1226	364104 ± 2210	92.0 %	7.3 %	0.8 %
12	25	46147 ± 405	14868 ± 476	304295 ± 1801	87.1 %	9.4 %	3.5 %
13	190	337003 ± 1509	85686 ± 1065	492000 ± 2847	91.4 %	7.8 %	0.8 %
15	25	48355 ± 415	17376 ± 526	387514 ± 2005	85.5 %	10.3 %	4.2 %
16	50	101995 ± 632	34724 ± 585	389570 ± 2033	87.9 %	10.0 %	2.1 %
18	25	50761 ± 573	18052 ± 500	367673 ± 2017	85.9 %	10.2 %	3.8 %
19	0	562 ± 185	1477 ± 333	299669 ± 1828	19.3 %	17.0 %	63.7 %
21	25	53809 ± 614	20054 ± 529	324737 ± 1992	86.0 %	10.7 %	3.2 %
23	50	106989 ± 860	38235 ± 907	359339 ± 2000	87.7 %	10.5 %	1.8 %
25	50	110343 ± 951	38480 ± 667	423267 ± 2299	87.7 %	10.2 %	2.1 %

#3

Table 1.3: CW-OSL signal intensities of all records for aliquot #3.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	151590 ± 862	36605 ± 751	132240 ± 1258	92.1 %	7.4 %	0.5 %
3	25	32961 ± 290	8051 ± 358	120849 ± 1228	90.5 %	7.4 %	2.0 %
4	50	65990 ± 452	15049 ± 526	134229 ± 1127	91.8 %	7.0 %	1.2 %
6	25	33396 ± 253	8295 ± 305	132612 ± 1069	90.3 %	7.5 %	2.2 %
7	85	111916 ± 473	25711 ± 473	174663 ± 1528	92.0 %	7.1 %	0.9 %
9	25	33970 ± 262	9167 ± 377	160790 ± 1289	89.3 %	8.1 %	2.6 %
10	155	2e+05 ± 841	45734 ± 772	253066 ± 1741	92.2 %	7.1 %	0.7 %
12	25	34588 ± 417	9890 ± 359	213193 ± 1353	88.2 %	8.4 %	3.4 %
13	190	252196 ± 990	56851 ± 1000	323978 ± 2018	92.3 %	7.0 %	0.7 %
15	25	35796 ± 432	11140 ± 424	260706 ± 1612	87.0 %	9.1 %	3.9 %
16	50	74986 ± 732	20673 ± 511	247957 ± 1573	89.9 %	8.3 %	1.8 %
18	25	37303 ± 436	11173 ± 485	237769 ± 1708	87.7 %	8.8 %	3.5 %
19	0	448 ± 136	341 ± 253	190893 ± 1361	25.7 %	6.6 %	67.8 %
21	25	38975 ± 395	11829 ± 390	206696 ± 1420	88.1 %	9.0 %	2.9 %
23	50	76130 ± 451	22220 ± 530	220549 ± 1617	89.6 %	8.8 %	1.6 %
25	50	79655 ± 664	22354 ± 526	254297 ± 1594	89.8 %	8.4 %	1.8 %

#4

Table 1.4: CW-OSL signal intensities of all records for aliquot #4.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	173447 ± 983	43822 ± 918	164460 ± 1451	91.7 %	7.8 %	0.5 %
3	25	38086 ± 316	10941 ± 385	148177 ± 1252	89.3 %	8.6 %	2.1 %
4	50	78339 ± 599	18375 ± 500	163128 ± 1423	91.6 %	7.2 %	1.2 %
6	25	38879 ± 328	10367 ± 378	163691 ± 1391	89.7 %	8.0 %	2.3 %
7	85	129573 ± 593	32761 ± 721	211563 ± 1309	91.3 %	7.7 %	0.9 %
9	25	38968 ± 367	10917 ± 421	201493 ± 1469	88.8 %	8.3 %	2.8 %
10	155	233399 ± 920	56368 ± 838	311034 ± 2011	91.8 %	7.4 %	0.8 %
12	25	40636 ± 321	11594 ± 405	266795 ± 1611	88.0 %	8.4 %	3.6 %
13	190	295257 ± 1362	69866 ± 1096	391927 ± 2226	92.0 %	7.3 %	0.8 %
15	25	42112 ± 309	13084 ± 417	321353 ± 1810	86.9 %	9.0 %	4.1 %
16	50	87827 ± 487	23628 ± 698	305322 ± 1774	90.0 %	8.1 %	1.9 %
18	25	43478 ± 473	13963 ± 488	288650 ± 1883	87.1 %	9.4 %	3.6 %
19	0	631 ± 156	-127 ± 306	236882 ± 1689	30.7 %	-2.1 %	71.3 %
21	25	46152 ± 418	13158 ± 551	256120 ± 1651	88.5 %	8.4 %	3.0 %
23	50	89145 ± 878	27122 ± 726	260371 ± 1819	89.3 %	9.1 %	1.6 %
25	50	93179 ± 511	27023 ± 606	306330 ± 1839	89.5 %	8.7 %	1.8 %

#5

Table 1.5: CW-OSL signal intensities of all records for aliquot #5.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	124411 ± 495	29736 ± 521	118997 ± 1125	92.1 %	7.4 %	0.5 %
3	25	25941 ± 244	6957 ± 303	110056 ± 1093	89.6 %	8.0 %	2.4 %
4	50	52353 ± 262	13116 ± 323	119007 ± 1128	91.1 %	7.6 %	1.3 %
6	25	26543 ± 278	6593 ± 375	119997 ± 1239	90.0 %	7.5 %	2.5 %
7	85	89043 ± 449	21162 ± 584	154965 ± 1398	91.7 %	7.3 %	1.0 %
9	25	26946 ± 258	6863 ± 306	143367 ± 1189	89.4 %	7.6 %	2.9 %
10	155	158721 ± 684	36764 ± 897	223906 ± 1639	92.1 %	7.1 %	0.8 %
12	25	27656 ± 416	7701 ± 341	188355 ± 1481	88.1 %	8.2 %	3.7 %
13	190	202005 ± 1202	41739 ± 1272	290893 ± 2309	92.8 %	6.4 %	0.8 %
15	25	28197 ± 311	8722 ± 399	228757 ± 1613	86.7 %	9.0 %	4.3 %
16	50	57979 ± 342	16300 ± 495	218975 ± 1491	89.5 %	8.4 %	2.1 %
18	25	29346 ± 359	8449 ± 334	210567 ± 1472	87.7 %	8.4 %	3.9 %
19	0	344 ± 146	226 ± 261	169196 ± 1369	23.4 %	5.2 %	71.4 %
21	25	30910 ± 337	8273 ± 384	183923 ± 1350	88.8 %	8.0 %	3.3 %
23	50	58844 ± 400	16777 ± 446	193649 ± 1384	89.6 %	8.6 %	1.8 %
25	50	60637 ± 306	18058 ± 424	219397 ± 1468	89.1 %	8.9 %	2.0 %

#6

Table 1.6: CW-OSL signal intensities of all records for aliquot #6.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	188620 ± 648	39907 ± 913	158668 ± 1451	92.9 %	6.6 %	0.5 %
3	25	41236 ± 423	8387 ± 363	147741 ± 1156	91.7 %	6.2 %	2.0 %
4	50	80523 ± 386	17452 ± 468	157120 ± 1259	92.2 %	6.7 %	1.1 %
6	25	40911 ± 411	9252 ± 455	160372 ± 1178	90.9 %	6.9 %	2.2 %
7	85	136991 ± 782	28829 ± 793	210722 ± 1566	92.6 %	6.5 %	0.9 %
9	25	41781 ± 356	9699 ± 430	201628 ± 1439	90.3 %	7.0 %	2.7 %
10	155	246516 ± 1235	49024 ± 1174	316237 ± 2155	93.1 %	6.2 %	0.7 %
12	25	42829 ± 473	9465 ± 452	264767 ± 1698	89.9 %	6.6 %	3.4 %
13	190	309498 ± 1774	58789 ± 1508	411612 ± 2764	93.3 %	5.9 %	0.8 %
15	25	44439 ± 409	10968 ± 441	322082 ± 1844	88.7 %	7.3 %	4.0 %
16	50	91431 ± 554	20878 ± 565	310767 ± 1786	91.1 %	7.0 %	1.9 %
18	25	46254 ± 372	11056 ± 466	292319 ± 1793	89.4 %	7.1 %	3.5 %
19	0	756 ± 161	-144 ± 241	238800 ± 1602	34.6 %	-2.2 %	67.6 %
21	25	47968 ± 276	11036 ± 442	260555 ± 1589	90.0 %	6.9 %	3.0 %
23	50	92345 ± 583	23499 ± 699	267948 ± 1625	90.6 %	7.7 %	1.6 %
25	50	96542 ± 492	23678 ± 574	319311 ± 2007	90.7 %	7.4 %	1.9 %

#7

Table 1.7: CW-OSL signal intensities of all records for aliquot #7.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	177716 ± 861	56533 ± 839	144018 ± 1708	90.0 %	9.6 %	0.4 %
3	25	38919 ± 425	13622 ± 420	141108 ± 1221	87.7 %	10.3 %	2.0 %
4	50	77106 ± 557	25559 ± 643	151157 ± 1505	89.0 %	9.9 %	1.1 %
6	25	39586 ± 420	13196 ± 386	156394 ± 1342	88.0 %	9.8 %	2.1 %
7	85	134036 ± 693	39720 ± 656	210410 ± 1555	90.2 %	9.0 %	0.9 %
9	25	39918 ± 379	14383 ± 463	194737 ± 1430	86.9 %	10.5 %	2.6 %
10	155	242283 ± 924	72610 ± 898	303514 ± 2062	90.2 %	9.1 %	0.7 %
12	25	42012 ± 464	15481 ± 435	260664 ± 1510	86.1 %	10.6 %	3.3 %
13	190	305152 ± 902	91830 ± 1039	398837 ± 2230	90.2 %	9.1 %	0.7 %
15	25	43086 ± 523	16473 ± 560	327923 ± 1811	85.1 %	10.9 %	4.0 %
16	50	91188 ± 903	34252 ± 654	312355 ± 1730	87.2 %	11.0 %	1.8 %
18	25	45992 ± 588	17066 ± 595	305130 ± 1929	85.8 %	10.7 %	3.5 %
19	0	932 ± 181	-190 ± 338	247019 ± 1719	38.9 %	-2.6 %	63.8 %
21	25	47835 ± 601	18545 ± 512	268480 ± 1654	85.9 %	11.2 %	3.0 %
23	50	93552 ± 1003	37322 ± 838	281089 ± 2171	86.8 %	11.6 %	1.6 %
25	50	95316 ± 691	38201 ± 762	325610 ± 1969	86.6 %	11.6 %	1.8 %

#8

Table 1.8: CW-OSL signal intensities of all records for aliquot #8.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	158184 ± 894	53321 ± 893	161618 ± 1760	89.3 %	10.1 %	0.6 %
3	25	35843 ± 350	11566 ± 322	149098 ± 1279	88.2 %	9.5 %	2.3 %
4	50	71413 ± 433	21556 ± 523	162533 ± 1366	89.7 %	9.1 %	1.3 %
6	25	35792 ± 225	11070 ± 282	166721 ± 1462	88.3 %	9.1 %	2.5 %
7	85	121380 ± 672	33706 ± 655	204252 ± 1741	90.6 %	8.4 %	0.9 %
9	25	37086 ± 298	11263 ± 399	197231 ± 1469	88.1 %	9.0 %	2.9 %
10	155	220631 ± 1085	57600 ± 890	304136 ± 2180	91.2 %	8.0 %	0.8 %
12	25	37021 ± 397	12605 ± 394	253269 ± 1784	86.5 %	9.9 %	3.7 %
13	190	271439 ± 1125	76259 ± 779	370745 ± 2185	90.7 %	8.5 %	0.8 %
15	25	38823 ± 305	12854 ± 385	311881 ± 1877	86.2 %	9.6 %	4.3 %
16	50	79705 ± 676	26759 ± 585	289622 ± 1872	88.1 %	9.9 %	2.0 %
18	25	40622 ± 361	13103 ± 452	284407 ± 1464	86.9 %	9.4 %	3.8 %
19	0	776 ± 130	-398 ± 251	230695 ± 1649	37.5 %	-6.4 %	68.9 %
21	25	43070 ± 297	12812 ± 468	247981 ± 1736	88.1 %	8.8 %	3.1 %
23	50	81770 ± 604	28379 ± 693	254395 ± 1850	88.1 %	10.2 %	1.7 %
25	50	84003 ± 752	29176 ± 662	297456 ± 1790	87.9 %	10.2 %	1.9 %

#9

Table 1.9: CW-OSL signal intensities of all records for aliquot #9.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	189573 ± 1019	65287 ± 781	157371 ± 1892	89.2 %	10.3 %	0.5 %
3	25	39392 ± 457	16307 ± 483	146109 ± 1288	86.1 %	11.9 %	2.0 %
4	50	81808 ± 856	27504 ± 627	156175 ± 1487	88.9 %	10.0 %	1.0 %
6	25	41521 ± 433	14152 ± 470	163207 ± 1286	87.8 %	10.0 %	2.1 %
7	85	137979 ± 887	49660 ± 938	200926 ± 1900	88.5 %	10.7 %	0.8 %
9	25	41648 ± 480	16821 ± 500	196365 ± 1613	85.9 %	11.6 %	2.5 %
10	155	251479 ± 1382	85883 ± 1197	304198 ± 2614	89.1 %	10.2 %	0.7 %
12	25	43229 ± 561	17937 ± 594	265141 ± 1698	85.0 %	11.8 %	3.2 %
13	190	321547 ± 1250	105919 ± 1261	394510 ± 2871	89.5 %	9.9 %	0.7 %
15	25	46768 ± 588	18688 ± 491	326665 ± 1801	85.0 %	11.4 %	3.7 %
16	50	95532 ± 856	38519 ± 890	316947 ± 2255	86.5 %	11.7 %	1.8 %
18	25	47577 ± 510	20010 ± 602	303685 ± 1771	84.7 %	11.9 %	3.4 %
19	0	786 ± 158	176 ± 293	247215 ± 1407	33.1 %	2.5 %	64.4 %
21	25	50998 ± 485	21236 ± 664	261671 ± 1936	85.4 %	11.9 %	2.7 %
23	50	98823 ± 862	39282 ± 836	278300 ± 2067	86.9 %	11.6 %	1.5 %
25	50	101879 ± 873	42287 ± 913	321687 ± 2291	86.3 %	12.0 %	1.7 %

#10

Table 1.10: CW-OSL signal intensities of all records for aliquot #10.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	140535 ± 808	48849 ± 914	123927 ± 1598	89.1 %	10.4 %	0.5 %
3	25	31029 ± 257	11478 ± 336	127127 ± 1195	87.0 %	10.8 %	2.2 %
4	50	61991 ± 341	21013 ± 435	137560 ± 1261	88.7 %	10.1 %	1.2 %
6	25	31200 ± 247	10702 ± 374	145911 ± 1319	87.4 %	10.0 %	2.5 %
7	85	105136 ± 661	33970 ± 762	181437 ± 1612	89.4 %	9.7 %	1.0 %
9	25	31951 ± 400	10978 ± 424	178124 ± 1395	87.0 %	10.0 %	3.0 %
10	155	189443 ± 681	57796 ± 924	267690 ± 1975	90.0 %	9.2 %	0.8 %
12	25	31820 ± 417	12616 ± 467	232717 ± 1472	84.9 %	11.3 %	3.8 %
13	190	230408 ± 1327	76925 ± 1162	327945 ± 2303	89.2 %	10.0 %	0.8 %
15	25	32947 ± 481	13353 ± 513	286816 ± 1778	84.1 %	11.4 %	4.5 %
16	50	67257 ± 537	26140 ± 602	266281 ± 1894	86.6 %	11.3 %	2.1 %
18	25	33427 ± 478	13997 ± 565	255391 ± 1869	84.2 %	11.8 %	4.0 %
19	0	533 ± 126	-292 ± 253	215595 ± 1597	30.1 %	-5.5 %	75.4 %
21	25	34482 ± 411	14212 ± 486	227173 ± 1569	84.8 %	11.7 %	3.5 %
23	50	67486 ± 711	28052 ± 741	231632 ± 1761	86.2 %	12.0 %	1.8 %
25	50	70861 ± 478	26714 ± 616	274373 ± 1708	86.9 %	11.0 %	2.1 %

#11

Table 1.11: CW-OSL signal intensities of all records for aliquot #11.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	145265 ± 631	47709 ± 789	113548 ± 1486	89.7 %	9.9 %	0.4 %
3	25	32071 ± 273	10526 ± 376	117277 ± 1194	88.3 %	9.7 %	2.0 %
4	50	63627 ± 531	21509 ± 549	124115 ± 1139	88.9 %	10.1 %	1.1 %
6	25	32473 ± 324	10740 ± 398	132623 ± 1297	88.0 %	9.8 %	2.2 %
7	85	108445 ± 991	34715 ± 668	166247 ± 1402	89.6 %	9.6 %	0.8 %
9	25	32038 ± 539	12470 ± 529	162836 ± 1474	86.1 %	11.2 %	2.7 %
10	155	196762 ± 910	59477 ± 745	242173 ± 1896	90.2 %	9.1 %	0.7 %
12	25	33258 ± 420	12281 ± 529	215127 ± 1749	85.9 %	10.6 %	3.4 %
13	190	245924 ± 917	73337 ± 876	307633 ± 2110	90.3 %	9.0 %	0.7 %
15	25	35249 ± 387	12632 ± 434	259055 ± 1655	85.8 %	10.3 %	3.9 %
16	50	71767 ± 625	27167 ± 571	242558 ± 1692	87.1 %	11.0 %	1.8 %
18	25	35718 ± 408	13919 ± 428	238994 ± 1819	85.3 %	11.1 %	3.5 %
19	0	447 ± 116	159 ± 241	197036 ± 1450	26.0 %	3.1 %	70.9 %
21	25	37921 ± 365	13565 ± 420	209629 ± 1665	86.7 %	10.4 %	3.0 %
23	50	74025 ± 551	27924 ± 586	214284 ± 1635	87.4 %	11.0 %	1.6 %
25	50	75852 ± 685	28529 ± 656	249395 ± 1661	87.2 %	11.0 %	1.8 %

#12

Table 1.12: CW-OSL signal intensities of all records for aliquot #12.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	175649 ± 743	42110 ± 705	150671 ± 1460	92.1 %	7.4 %	0.5 %
3	25	40827 ± 381	10437 ± 377	137525 ± 1286	90.4 %	7.7 %	1.9 %
4	50	82545 ± 418	21615 ± 488	161638 ± 1381	90.9 %	8.0 %	1.1 %
6	25	42588 ± 270	11174 ± 363	164816 ± 1309	90.0 %	7.9 %	2.1 %
7	85	144448 ± 431	38115 ± 599	242554 ± 1639	91.0 %	8.0 %	1.0 %
9	25	44132 ± 396	12513 ± 484	225772 ± 1509	88.8 %	8.4 %	2.8 %
10	155	268901 ± 1204	70280 ± 1374	406235 ± 2928	91.2 %	8.0 %	0.8 %
12	25	46425 ± 435	17103 ± 504	327423 ± 1932	85.7 %	10.6 %	3.7 %
13	190	342383 ± 1975	101706 ± 1407	576214 ± 3698	90.1 %	9.0 %	0.9 %
15	25	50344 ± 432	18931 ± 573	441750 ± 2133	84.7 %	10.7 %	4.6 %
16	50	103839 ± 1146	42383 ± 781	454716 ± 2463	85.9 %	11.7 %	2.3 %
18	25	51782 ± 626	23579 ± 579	421287 ± 2107	83.1 %	12.7 %	4.2 %
19	0	552 ± 202	2410 ± 365	344145 ± 1950	15.8 %	23.1 %	61.0 %
21	25	55757 ± 578	24527 ± 642	389283 ± 2037	84.0 %	12.4 %	3.6 %
23	50	110076 ± 1171	49865 ± 882	444946 ± 2427	85.0 %	12.9 %	2.1 %
25	50	113014 ± 1002	50224 ± 819	517640 ± 3073	85.0 %	12.6 %	2.4 %

#13

Table 1.13: CW-OSL signal intensities of all records for aliquot #13.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	162637 ± 833	33060 ± 1037	148241 ± 1588	93.1 %	6.3 %	0.5 %
3	25	36644 ± 335	7235 ± 381	130157 ± 1076	91.9 %	6.1 %	2.0 %
4	50	73182 ± 384	14855 ± 492	141768 ± 1291	92.6 %	6.3 %	1.1 %
6	25	37192 ± 301	8099 ± 359	141919 ± 1156	91.2 %	6.6 %	2.1 %
7	85	123630 ± 846	26196 ± 709	185064 ± 1359	92.6 %	6.6 %	0.9 %
9	25	37705 ± 435	8408 ± 444	172976 ± 1270	90.7 %	6.8 %	2.6 %
10	155	224143 ± 1058	45405 ± 970	275857 ± 1936	93.0 %	6.3 %	0.7 %
12	25	39210 ± 292	9499 ± 378	230207 ± 1437	89.5 %	7.3 %	3.2 %
13	190	280278 ± 1212	57171 ± 1477	353518 ± 2395	92.9 %	6.4 %	0.7 %
15	25	39969 ± 559	11692 ± 540	279132 ± 1686	87.6 %	8.6 %	3.8 %
16	50	83163 ± 467	20839 ± 652	272131 ± 1750	90.6 %	7.6 %	1.8 %
18	25	42248 ± 334	10527 ± 395	257295 ± 1604	89.2 %	7.4 %	3.4 %
19	0	539 ± 130	244 ± 244	207748 ± 1394	28.3 %	4.3 %	67.4 %
21	25	43187 ± 421	12227 ± 469	224709 ± 1585	88.7 %	8.4 %	2.9 %
23	50	85045 ± 642	22625 ± 666	240548 ± 1710	90.4 %	8.0 %	1.6 %
25	50	88793 ± 655	22288 ± 527	276521 ± 1750	90.6 %	7.6 %	1.8 %

#14

Table 1.14: CW-OSL signal intensities of all records for aliquot #14.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	174169 ± 1166	21171 ± 1595	191793 ± 2213	95.5 %	3.9 %	0.6 %
3	25	38285 ± 414	4980 ± 421	155908 ± 1261	93.6 %	4.1 %	2.4 %
4	50	76548 ± 571	8873 ± 739	171745 ± 1493	95.0 %	3.7 %	1.3 %
6	25	38603 ± 441	5011 ± 499	167250 ± 1244	93.4 %	4.1 %	2.5 %
7	85	130012 ± 1116	14643 ± 1153	226687 ± 2188	95.4 %	3.6 %	1.0 %
9	25	38713 ± 395	6591 ± 468	202543 ± 1338	91.8 %	5.2 %	3.0 %
10	155	232822 ± 1747	26447 ± 2009	332717 ± 3077	95.5 %	3.6 %	0.8 %
12	25	40189 ± 354	6619 ± 486	265860 ± 1507	91.2 %	5.0 %	3.7 %
13	190	292092 ± 2539	33351 ± 2792	434726 ± 4175	95.5 %	3.6 %	0.9 %
15	25	41423 ± 393	7829 ± 515	329363 ± 2117	89.9 %	5.7 %	4.4 %
16	50	87403 ± 710	13595 ± 715	326729 ± 1967	93.0 %	4.8 %	2.1 %
18	25	43891 ± 465	7998 ± 534	302009 ± 1800	90.6 %	5.5 %	3.9 %
19	0	759 ± 175	253 ± 326	243697 ± 1460	32.3 %	3.6 %	64.1 %
21	25	45765 ± 505	8347 ± 478	272235 ± 1636	91.1 %	5.6 %	3.4 %
23	50	89432 ± 610	16561 ± 657	296239 ± 1741	92.4 %	5.7 %	1.9 %
25	50	91941 ± 725	15998 ± 817	341644 ± 2047	92.5 %	5.4 %	2.1 %

#15

Table 1.15: CW-OSL signal intensities of all records for aliquot #15.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	122229 ± 390	33685 ± 437	110417 ± 1208	91.1 %	8.4 %	0.5 %
3	25	26427 ± 306	7724 ± 327	104554 ± 1163	89.1 %	8.7 %	2.2 %
4	50	52849 ± 360	14925 ± 383	115147 ± 1551	90.2 %	8.5 %	1.2 %
6	25	26087 ± 306	8836 ± 328	116346 ± 1051	87.6 %	9.9 %	2.4 %
7	85	90377 ± 696	25579 ± 533	151384 ± 1296	90.5 %	8.6 %	0.9 %
9	25	27563 ± 292	8761 ± 331	142875 ± 1160	87.8 %	9.4 %	2.8 %
10	155	163869 ± 574	44238 ± 688	226137 ± 1691	91.0 %	8.2 %	0.8 %
12	25	27913 ± 390	10009 ± 348	191442 ± 1437	86.0 %	10.3 %	3.6 %
13	190	206392 ± 968	55469 ± 797	293362 ± 1813	91.0 %	8.2 %	0.8 %
15	25	28776 ± 231	10575 ± 379	240865 ± 1638	85.1 %	10.5 %	4.4 %
16	50	61742 ± 251	19523 ± 376	231971 ± 1406	88.6 %	9.4 %	2.1 %
18	25	31094 ± 392	10147 ± 405	222830 ± 1410	86.7 %	9.5 %	3.8 %
19	0	486 ± 141	90 ± 244	180287 ± 1375	29.8 %	1.8 %	68.3 %
21	25	31957 ± 433	11979 ± 437	192477 ± 1543	86.0 %	10.8 %	3.2 %
23	50	62861 ± 588	22251 ± 536	202792 ± 1470	87.8 %	10.4 %	1.8 %
25	50	64353 ± 560	21859 ± 443	238476 ± 1549	88.0 %	10.0 %	2.0 %

#16

Table 1.16: CW-OSL signal intensities of all records for aliquot #16.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	147278 ± 522	44232 ± 760	134352 ± 1476	90.4 %	9.1 %	0.5 %
3	25	33447 ± 217	10355 ± 308	135376 ± 1141	88.6 %	9.2 %	2.2 %
4	50	65092 ± 520	21328 ± 489	140964 ± 1382	89.0 %	9.8 %	1.2 %
6	25	32858 ± 379	11390 ± 340	144692 ± 1227	87.5 %	10.2 %	2.4 %
7	85	112236 ± 755	33620 ± 694	184969 ± 1476	90.0 %	9.0 %	0.9 %
9	25	33379 ± 319	11411 ± 363	175024 ± 1358	87.2 %	10.0 %	2.8 %
10	155	200041 ± 1068	60044 ± 774	260207 ± 1950	90.2 %	9.1 %	0.7 %
12	25	34475 ± 334	12131 ± 462	228245 ± 1629	86.3 %	10.2 %	3.5 %
13	190	250492 ± 869	73323 ± 900	332005 ± 1961	90.4 %	8.9 %	0.7 %
15	25	36002 ± 318	13007 ± 437	276307 ± 1725	85.6 %	10.3 %	4.1 %
16	50	73165 ± 664	27609 ± 586	260729 ± 1733	87.1 %	11.0 %	1.9 %
18	25	37028 ± 522	13901 ± 532	252961 ± 1696	85.6 %	10.8 %	3.6 %
19	0	714 ± 147	-288 ± 273	213147 ± 1376	36.9 %	-5.0 %	68.1 %
21	25	39043 ± 517	14268 ± 531	220479 ± 1651	86.4 %	10.6 %	3.0 %
23	50	76242 ± 698	26281 ± 592	233747 ± 1588	88.1 %	10.2 %	1.7 %
25	50	78919 ± 482	27247 ± 546	267142 ± 1708	88.0 %	10.2 %	1.8 %

#17

Table 1.17: CW-OSL signal intensities of all records for aliquot #17.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	180313 ± 1113	32057 ± 1063	158191 ± 1678	93.9 %	5.6 %	0.5 %
3	25	38177 ± 349	8377 ± 361	134908 ± 1256	91.3 %	6.7 %	2.0 %
4	50	75477 ± 636	17301 ± 561	149820 ± 1294	91.8 %	7.0 %	1.1 %
6	25	39203 ± 294	7488 ± 454	151469 ± 1194	91.9 %	5.9 %	2.2 %
7	85	130532 ± 727	25331 ± 886	203859 ± 1508	93.0 %	6.0 %	0.9 %
9	25	39282 ± 458	8663 ± 483	188900 ± 1371	90.6 %	6.7 %	2.7 %
10	155	237157 ± 1552	43197 ± 1722	309506 ± 2334	93.5 %	5.7 %	0.8 %
12	25	40077 ± 458	10278 ± 486	249534 ± 1742	88.9 %	7.6 %	3.4 %
13	190	292156 ± 1534	56668 ± 1931	392290 ± 2737	93.2 %	6.0 %	0.8 %
15	25	41432 ± 426	11421 ± 440	304938 ± 1580	87.9 %	8.1 %	4.0 %
16	50	87286 ± 552	21389 ± 648	293393 ± 1781	90.7 %	7.4 %	1.9 %
18	25	43296 ± 500	12123 ± 808	279668 ± 1770	88.2 %	8.3 %	3.5 %
19	0	541 ± 145	10 ± 236	224636 ± 1537	28.0 %	0.2 %	71.9 %
21	25	45193 ± 630	12543 ± 452	244958 ± 1568	88.8 %	8.2 %	3.0 %
23	50	89775 ± 551	22130 ± 586	260386 ± 1645	90.9 %	7.5 %	1.6 %
25	50	90873 ± 608	24845 ± 556	297270 ± 1775	90.0 %	8.2 %	1.8 %

#18

Table 1.18: CW-OSL signal intensities of all records for aliquot #18.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	147370 ± 567	38343 ± 588	125611 ± 1258	91.5 %	8.0 %	0.5 %
3	25	30214 ± 265	9033 ± 267	114312 ± 1213	89.0 %	8.9 %	2.1 %
4	50	61854 ± 394	16279 ± 448	120333 ± 1143	90.9 %	8.0 %	1.1 %
6	25	31519 ± 418	8056 ± 331	125737 ± 1073	90.1 %	7.7 %	2.2 %
7	85	105488 ± 398	26855 ± 494	166341 ± 1469	91.3 %	7.8 %	0.9 %
9	25	31797 ± 375	9106 ± 437	154319 ± 1323	88.8 %	8.5 %	2.7 %
10	155	189648 ± 1075	46397 ± 916	241334 ± 1637	91.8 %	7.5 %	0.7 %
12	25	32510 ± 330	9854 ± 327	202479 ± 1469	87.7 %	8.9 %	3.4 %
13	190	236834 ± 1159	60179 ± 947	307755 ± 1851	91.5 %	7.8 %	0.7 %
15	25	34696 ± 409	10271 ± 444	252424 ± 1590	87.4 %	8.7 %	3.9 %
16	50	70036 ± 372	21672 ± 463	239502 ± 1661	88.9 %	9.2 %	1.9 %
18	25	35121 ± 540	12250 ± 540	228015 ± 1390	86.4 %	10.1 %	3.5 %
19	0	768 ± 144	-397 ± 262	185893 ± 1445	43.0 %	-7.4 %	64.4 %
21	25	36330 ± 454	12985 ± 494	200350 ± 1514	86.7 %	10.4 %	3.0 %
23	50	70583 ± 767	25094 ± 609	208066 ± 1624	87.9 %	10.5 %	1.6 %
25	50	72909 ± 664	25654 ± 634	244890 ± 1597	87.8 %	10.3 %	1.8 %

#19

Table 1.19: CW-OSL signal intensities of all records for aliquot #19.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	150560 ± 968	36162 ± 662	127198 ± 1293	92.1 %	7.4 %	0.5 %
3	25	32101 ± 412	9043 ± 326	117807 ± 1187	89.5 %	8.4 %	2.0 %
4	50	64614 ± 403	17507 ± 411	128910 ± 1264	90.7 %	8.2 %	1.1 %
6	25	33388 ± 264	8057 ± 383	133370 ± 1335	90.5 %	7.3 %	2.2 %
7	85	110371 ± 836	26440 ± 634	171223 ± 1266	91.8 %	7.4 %	0.9 %
9	25	32962 ± 509	9649 ± 377	160030 ± 1192	88.6 %	8.7 %	2.7 %
10	155	200224 ± 997	45399 ± 913	251567 ± 1824	92.3 %	7.0 %	0.7 %
12	25	33828 ± 384	10468 ± 425	208516 ± 1512	87.6 %	9.1 %	3.3 %
13	190	249181 ± 956	57759 ± 1019	318339 ± 2150	92.1 %	7.1 %	0.7 %
15	25	36359 ± 383	11638 ± 468	250241 ± 1665	87.0 %	9.3 %	3.7 %
16	50	74334 ± 582	22183 ± 560	240521 ± 1771	89.3 %	8.9 %	1.8 %
18	25	37824 ± 408	11106 ± 492	233144 ± 1424	88.0 %	8.6 %	3.4 %
19	0	796 ± 139	-340 ± 278	191635 ± 1365	42.6 %	-6.1 %	63.5 %
21	25	39172 ± 316	11294 ± 455	201865 ± 1536	88.6 %	8.6 %	2.8 %
23	50	76317 ± 463	22528 ± 597	213574 ± 1751	89.6 %	8.8 %	1.6 %
25	50	79341 ± 783	23852 ± 598	246486 ± 1480	89.3 %	9.0 %	1.7 %

#20

Table 1.20: CW-OSL signal intensities of all records for aliquot #20.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	167390 ± 624	44944 ± 606	128427 ± 1276	91.4 %	8.2 %	0.4 %
3	25	37336 ± 259	10384 ± 366	123714 ± 1051	89.8 %	8.4 %	1.8 %
4	50	74859 ± 437	20377 ± 541	136458 ± 1313	90.7 %	8.3 %	1.0 %
6	25	37739 ± 389	10868 ± 381	142081 ± 1259	89.3 %	8.6 %	2.1 %
7	85	126198 ± 876	36020 ± 648	182164 ± 1598	90.5 %	8.6 %	0.8 %
9	25	38729 ± 430	11033 ± 459	176769 ± 1386	89.0 %	8.5 %	2.5 %
10	155	227454 ± 737	62495 ± 730	274358 ± 1978	91.0 %	8.4 %	0.7 %
12	25	39975 ± 327	12280 ± 429	238381 ± 1372	87.7 %	9.0 %	3.2 %
13	190	289583 ± 1104	72291 ± 1406	369099 ± 2344	91.6 %	7.7 %	0.7 %
15	25	40134 ± 471	14976 ± 513	292526 ± 1709	85.5 %	10.7 %	3.8 %
16	50	84937 ± 608	27904 ± 650	282222 ± 1828	88.4 %	9.7 %	1.8 %
18	25	42993 ± 472	14400 ± 515	269936 ± 1709	86.9 %	9.7 %	3.4 %
19	0	680 ± 163	-62 ± 311	220448 ± 1607	33.6 %	-1.0 %	67.4 %
21	25	44650 ± 498	15410 ± 524	237833 ± 1634	87.1 %	10.1 %	2.9 %
23	50	87530 ± 1010	30835 ± 618	250825 ± 1717	88.0 %	10.4 %	1.6 %
25	50	90060 ± 870	31441 ± 590	290252 ± 1916	88.0 %	10.3 %	1.8 %

#21

Table 1.21: CW-OSL signal intensities of all records for aliquot #21.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	171232 ± 1085	30551 ± 938	166321 ± 1572	93.8 %	5.6 %	0.6 %
3	25	37850 ± 386	7407 ± 364	134334 ± 1148	92.0 %	6.0 %	2.0 %
4	50	76012 ± 543	12943 ± 739	149923 ± 1346	93.5 %	5.3 %	1.1 %
6	25	38099 ± 287	8242 ± 349	142803 ± 1134	91.3 %	6.6 %	2.1 %
7	85	129458 ± 1015	23362 ± 783	194634 ± 1527	93.5 %	5.6 %	0.9 %
9	25	39821 ± 375	8352 ± 418	175580 ± 1424	91.1 %	6.4 %	2.5 %
10	155	234410 ± 1275	40825 ± 1423	293262 ± 2345	93.8 %	5.5 %	0.7 %
12	25	40445 ± 384	9711 ± 493	233687 ± 1514	89.6 %	7.2 %	3.2 %
13	190	293336 ± 1830	52622 ± 1723	375474 ± 2735	93.6 %	5.6 %	0.7 %
15	25	42634 ± 407	9799 ± 473	291433 ± 1617	89.3 %	6.9 %	3.8 %
16	50	87036 ± 654	21076 ± 658	275955 ± 1680	90.8 %	7.4 %	1.8 %
18	25	44079 ± 487	11714 ± 498	260105 ± 1554	88.8 %	7.9 %	3.2 %
19	0	465 ± 168	219 ± 312	213227 ± 1282	25.0 %	3.9 %	71.0 %
21	25	45789 ± 402	12203 ± 486	226840 ± 1628	89.3 %	8.0 %	2.7 %
23	50	90269 ± 560	22286 ± 671	246756 ± 1699	90.9 %	7.5 %	1.5 %
25	50	93050 ± 607	23932 ± 616	283759 ± 1882	90.5 %	7.8 %	1.7 %

#22

Table 1.22: CW-OSL signal intensities of all records for aliquot #22.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	185921 ± 1158	26881 ± 1165	170210 ± 1767	94.9 %	4.6 %	0.5 %
3	25	40092 ± 388	6753 ± 410	149186 ± 1151	92.6 %	5.2 %	2.1 %
4	50	80705 ± 604	10581 ± 703	165726 ± 1408	94.6 %	4.2 %	1.2 %
6	25	40583 ± 421	7123 ± 464	164312 ± 1281	92.3 %	5.4 %	2.3 %
7	85	138621 ± 1019	17782 ± 1291	230141 ± 2030	95.0 %	4.1 %	1.0 %
9	25	41286 ± 413	8054 ± 445	199877 ± 1445	91.3 %	6.0 %	2.7 %
10	155	250763 ± 1922	31401 ± 2196	347253 ± 3551	95.2 %	4.0 %	0.8 %
12	25	42739 ± 450	8664 ± 485	267756 ± 1637	90.4 %	6.1 %	3.5 %
13	190	320110 ± 2831	38191 ± 3134	454371 ± 4627	95.4 %	3.8 %	0.8 %
15	25	45285 ± 365	9958 ± 417	330650 ± 1695	89.4 %	6.6 %	4.0 %
16	50	96294 ± 675	16776 ± 850	330392 ± 2026	92.6 %	5.4 %	2.0 %
18	25	48330 ± 383	10319 ± 413	309637 ± 1854	90.0 %	6.4 %	3.6 %
19	0	802 ± 162	137 ± 315	246476 ± 1796	33.8 %	1.9 %	64.3 %
21	25	50871 ± 501	10637 ± 576	272711 ± 1665	90.6 %	6.4 %	3.0 %
23	50	99878 ± 625	20872 ± 616	298428 ± 2015	91.9 %	6.4 %	1.7 %
25	50	103292 ± 864	20939 ± 778	346255 ± 1963	91.9 %	6.2 %	1.9 %

#23

Table 1.23: CW-OSL signal intensities of all records for aliquot #23.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	174639 ± 847	38662 ± 1015	151953 ± 1389	92.6 %	6.9 %	0.5 %
3	25	38179 ± 396	8945 ± 434	141381 ± 1140	90.8 %	7.1 %	2.1 %
4	50	75187 ± 679	17793 ± 582	153654 ± 1353	91.6 %	7.3 %	1.2 %
6	25	38240 ± 502	9166 ± 472	157560 ± 1345	90.4 %	7.3 %	2.3 %
7	85	127366 ± 492	29269 ± 759	202558 ± 1501	92.0 %	7.1 %	0.9 %
9	25	39068 ± 404	9701 ± 406	188362 ± 1361	89.8 %	7.5 %	2.7 %
10	155	232335 ± 1014	46459 ± 1057	303900 ± 2192	93.0 %	6.2 %	0.8 %
12	25	39962 ± 371	10977 ± 397	247968 ± 1581	88.5 %	8.1 %	3.4 %
13	190	292910 ± 1855	58065 ± 1301	389165 ± 2684	93.1 %	6.2 %	0.8 %
15	25	41090 ± 334	12203 ± 506	303344 ± 1816	87.3 %	8.7 %	4.0 %
16	50	86522 ± 721	21705 ± 570	296063 ± 1542	90.5 %	7.6 %	1.9 %
18	25	43342 ± 392	12224 ± 401	278231 ± 1801	88.2 %	8.3 %	3.5 %
19	0	487 ± 167	317 ± 323	226720 ± 1643	24.4 %	5.3 %	70.3 %
21	25	45641 ± 490	12222 ± 397	245004 ± 1660	89.1 %	8.0 %	3.0 %
23	50	89267 ± 477	23536 ± 651	267155 ± 1693	90.4 %	8.0 %	1.7 %
25	50	91857 ± 657	24875 ± 507	304594 ± 1729	90.0 %	8.2 %	1.8 %

#24

Table 1.24: CW-OSL signal intensities of all records for aliquot #24.

		Component intensity \(n ± \sigma\)			Contribution to first channel
record index	dose (s)	Fast	Medium	Slow2	Fast	Medium	Slow2
1	?	177289 ± 672	42978 ± 653	144948 ± 1350	92.1 %	7.5 %	0.5 %
3	25	38223 ± 308	9604 ± 357	138285 ± 1214	90.4 %	7.6 %	2.0 %
4	50	77605 ± 649	19122 ± 537	151762 ± 1149	91.4 %	7.5 %	1.1 %
6	25	39519 ± 333	9885 ± 379	155755 ± 1182	90.2 %	7.6 %	2.2 %
7	85	134701 ± 860	30335 ± 861	216274 ± 1634	92.1 %	7.0 %	0.9 %
9	25	40017 ± 404	11040 ± 384	200494 ± 1550	89.0 %	8.2 %	2.8 %
10	155	242481 ± 1124	58216 ± 862	323506 ± 2252	91.9 %	7.4 %	0.8 %
12	25	42113 ± 497	12348 ± 444	269570 ± 1732	87.9 %	8.6 %	3.5 %
13	190	306582 ± 1664	74174 ± 1245	431432 ± 2854	91.8 %	7.4 %	0.8 %
15	25	43946 ± 433	14192 ± 489	340201 ± 1693	86.5 %	9.4 %	4.1 %
16	50	90112 ± 696	28732 ± 619	337604 ± 1841	88.5 %	9.4 %	2.0 %
18	25	45919 ± 378	15435 ± 461	316810 ± 1707	86.6 %	9.7 %	3.7 %
19	0	719 ± 159	344 ± 301	260553 ± 1651	29.4 %	4.7 %	65.9 %
21	25	48555 ± 418	15179 ± 505	287226 ± 1669	87.6 %	9.2 %	3.2 %
23	50	94122 ± 1022	31794 ± 703	311614 ± 1904	88.2 %	10.0 %	1.8 %
25	50	96938 ± 1034	33641 ± 816	361473 ± 2058	87.8 %	10.2 %	2.0 %

Methods

Basics

The data analysis approach assumes that all CW-OSL signal curves in a data set can be sufficently described by a sum of exponential decays with first order kinetics:

\[I(t) = \sum_{i=1}^K n_ie^{-\lambda_it}\]

Here, \(I(t)\) is the CW-OSL signal, \(K\) is the number of signal components. Each component is defined by its intensity \(n\) and its decay rate \(\lambda\).

To achieve component-seperated dose information, the data analysis process is divided into multiple steps:

Step 2 workflow

The decompositon process in Step 2 calculates the signal component intensities \(n_1\) to \(n_K\) and their statistical errors \(\sigma\). An intensity value \(n\) is equal to the areas under the signal component curve, including the ongoing virtual signal beyond the stimulation end. The algorithm can be outlined as following:

Detailed explanations on step 2.1 and step 2.2 can be found in Mittelstraß (2019). Practical advice can be found in the R documentations ?optimise_OSLintervals and ?decompose_OSLcurve.

Signal bin intervals

In step 2.1, optimum signal bin intervals for the decomposition process are defined. optimise_OSLintervals() determines iteratively those integration intervals which asure the maximum mathematical independence between the signal components. This minimizes the impact of eventual systematic errors from the Step 1 data analysis. The optimise_OSLintervals() output depends just on the decay rates and the channel settings, not on the specific CW-OSL measurements.

Table 2: Time intervals for the signal bin integration \(I_{bin}=\int_{t_{start}}^{t_{end}}I(t)dt\). If the components decay rates were determined by a global average curve, the global average bin values are given by \(\overline{I}_{bin}\)

signal bin	\(t_{start}\)	\(t_{end}\)	channels	\(\overline{I}_{bin}\)
1	0.0	1.0	10	85485
2	1.0	5.6	46	41308
3	5.6	20.0	144	24181

Step 2.1 used the differential evolution approach of Price et al. (2006)), implemented by Ardia et al. (2020) in the function DEoptim from the package DEoptim.

Decomposition

In Step 2.2, the decomposition process is performed by decompose_OSLcurve() for each record saperately and can be summarized as following:

1. Set value of each signal bin equal to the integrated basic equation. Obtain a system of \(K\) equations with just the component intensities \(n_1\) to \(n_K\) unknown.
2. Solve equation system by Cramer’s rule and get \(n_1\) to \(n_K\).
3. (Optional) Refine \(n_1\) to \(n_K\) by a Levenberg-Marquardt fitting performed by the function nlsLM from the package minpack.lm (Elzhov et al. 2016). This step enhances the result reliability in presence of an uncorrected signal background level but decreases the error estimation accuracy.
4. Calculate the standard deviation for each signal bin by comparing the modelled CW-OSL curve with the actual measurement data.
5. Apply the propagation of uncertainty for independent variables method onto Cramer’s rule and calculate the standard deviation \(\sigma_1\) to \(\sigma_K\) for each component intensity.

This approach (despite step 3) was first presented by Mittelstraß et al. (2014) and further developed by Mittelstraß (2019). An introduction in the mathematical concepts can be found for example in Bevington and Robinson (1992) or Bronštejn et al. (1989).

Error estimation

The intensity error values \(\sigma\) account just for statistical noise (e. g. shot noise and dark current noise). Other sources of error are not sufficiently accounted for. But they might be detectable through not-statistical scattering in the residual plot (lower plot in figure 1). Is the residual plot dominated by statistical noise, the CW.OSL model can be expected as valid. If not, one of the following systematic errors is likely:

• Incorrect decay constants
• Not-first-order decay kinetics
• Insufficient number of signal components (under-fitting)
• Unstable light stimulation
• Inadequate or missing background correction

General remarks

This report was automatically generated by functions of the R package OSLdecomposition written and maintained by Dirk Mittelstraß (dirk.mittelstrass@luminescence.de). This package is preferably used together with the R package Luminescence (link) by Kreutzer et al. (2012). For the dynamic creation of this HTML report, the R packages knitr and rmarkdown are used, see Xie (2015) and Xie et al. (2018). All diagrams are drawn with ggplot2 (Wickham 2016).

The data analysis method covered with this report was developed for CW-OSL SAR protocol measured quartz samples, see Murray and Wintle (2000). It might be also useful for the analysis of Al2O3 or feldspar samples measured with SAR-like protocols. You can use, share and publish this report and the containing results at will. But we demand to refer to the R package OSLdecomposition including its version number (21.2.15) if you publish the results in a peer-reviewed journal. Please include the following reference to your publication:

---

Mittelstraß, D., 2019. Decomposition of weak optically stimulated luminescence signals and its application in retrospective dosimetry at quartz, Master thesis, TU Dresden, Dresden.
https://iktp.tu-dresden.de/IKTP/pub/19/Dirk_Mittelstrass_Master.pdf
(reference will be replaced as soon as peer-reviewed publication is accepted)

---

We allow and encourage you to add the HTML file of this report to the electronical supplement of your publication.

References

Ardia, D., Mullen, K.M., Peterson, B.G., Ulrich, J., 2020. DEoptim: Differential Evolution in R.
https://CRAN.R-project.org/package=DEoptim

Bevington, P.R., Robinson, D.K., 1992. Data reduction and error analysis for the physical sciences, 2nd ed. ed. McGraw-Hill, New York.
http://hosting.astro.cornell.edu/academics/courses/astro3310/Books/Bevington_opt.pdf

Bronštejn, I.N., Semendjaev, K.A., Grosche, G., 1989. Taschenbuch der Mathematik, 24. Aufl. ed. Verl. Nauka, Moskau.

Elzhov, T.V., Mullen, K.M., Spiess, A.-N., Bolker, B., 2016. minpack.lm: R Interface to the Levenberg-Marquardt Nonlinear Least-Squares Algorithm Found in MINPACK, Plus Support for Bounds.
https://CRAN.R-project.org/package=minpack.lm

Kreutzer, S., Schmidt, C., Fuchs, M.C., Dietze, M., Fuchs, M., 2012. Introducing an R package for luminescence dating analysis. Ancient TL 30.
http://ancienttl.org/ATL_30-1_2012/ATL_30-1_Kreutzer_p1-8.pdf

Mittelstraß, D., Kreutzer, S., Richter, D., 2013. An analytical method to separate signal components of CW-OSL measurements of quartz. talk presented at the Central European Conference on Luminescence and Trapped-Charge dating: DLED 2013, Freiberg, Germany

Mittelstraß D., Schmidt C., Beyer J., and Straessner A., 2019. Automated identification and separation of quartz CW-OSL signal components with R. talk presented at the Central European Conference on Luminescence and Trapped-Charge dating: DLED 2019, Bingen, Germany
http://luminescence.de/OSLdecomp_talk.pdf

Murray, A. S. and Wintle, A. G., 2000. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiation Measurements 32
https://doi.org/10.1016/S1350-4487(99)00253-X

Price, K.V., Storn, R.M., Lampinen, J.A., 2006. Differential Evolution - A Practical Approach to Global Optimization, Natural Computing. Springer-Verlag.

Wickham, H., 2016. ggplot2: elegant graphics for data analysis, Second edition. ed, Use R! Springer, Cham.
https://ggplot2.tidyverse.org/

Xie, Y., 2015. Dynamic documents with R and Knitr, Second edition. ed. CRC Press/Taylor & Francis, Boca Raton.
https://yihui.org/knitr/

Xie, Y., Allaire, J.J., Grolemund, G., 2018. R Markdown: the definitive guide. Taylor & Francis, CRC Press, Boca Raton.
https://bookdown.org/yihui/rmarkdown/

Script parameters

Input data
Analyzed data set	data_set
Analyzed record type	OSL
No. of measurements	384
No. of aliquots	24
Algorithm parameter
Applied algorithm	det+nls
Error estimation approach	empiric
With background component?	FALSE

Last major change in the script code: 2021-02-15 by D. Mittelstraß

Session info:

R version 4.0.3 (2020-10-10) Platform: i386-w64-mingw32/i386 (32-bit) Running under: Windows 10 x64 (build 19041)

Matrix products: default

attached base packages: [1] stats graphics grDevices utils datasets methods base

other attached packages: [1] scales_1.1.1 gridExtra_2.3 ggplot2_3.3.3
[4] knitr_1.31 kableExtra_1.3.1 Luminescence_0.9.10
[7] magrittr_2.0.1 OSLdecomposition_21-02-15

loaded via a namespace (and not attached): [1] Rcpp_1.0.6 lattice_0.20-41 tidyr_1.1.2 digest_0.6.27
[5] cellranger_1.1.0 R6_2.5.0 backports_1.2.1 evaluate_0.14
[9] httr_1.4.2 highr_0.8 pillar_1.4.7 rlang_0.4.10
[13] readxl_1.3.1 curl_4.3 rstudioapi_0.13 data.table_1.13.6 [17] car_3.0-10 raster_3.4-5 rmarkdown_2.6 labeling_0.4.2
[21] webshot_0.5.2 stringr_1.4.0 foreign_0.8-80 munsell_0.5.0
[25] broom_0.7.3 compiler_4.0.3 xfun_0.20 pkgconfig_2.0.3
[29] htmltools_0.5.1.1 tidyselect_1.1.0 tibble_3.0.5 rio_0.5.16
[33] codetools_0.2-16 viridisLite_0.3.0 crayon_1.3.4 dplyr_1.0.3
[37] withr_2.4.1 ggpubr_0.4.0 grid_4.0.3 DEoptim_2.2-5
[41] gtable_0.3.0 lifecycle_0.2.0 zip_2.1.1 carData_3.0-4
[45] stringi_1.5.3 farver_2.0.3 ggsignif_0.6.0 sp_1.4-5
[49] xml2_1.3.2 ellipsis_0.3.1 generics_0.1.0 vctrs_0.3.6
[53] cowplot_1.1.1 openxlsx_4.2.3 tools_4.0.3 forcats_0.5.0
[57] glue_1.4.2 purrr_0.3.4 hms_1.0.0 abind_1.4-5
[61] parallel_4.0.3 yaml_2.2.1 colorspace_2.0-0 rstatix_0.6.0
[65] rvest_0.3.6 minpack.lm_1.2-1 haven_2.3.1