The current paper is addressed mainly to users of luminescence results obtained in the Gliwice Luminescence Laboratory (GLL) and other researchers who would like to set up their own laboratories. We describe procedures necessary for calibration of the measuring devices used in our laboratory. Technical details of such nature are not always easy to find in articles related to the presentation and interpretation of luminescence results themselves.
Currently, numerical dating using luminescence methods is widely applied in geosciences and archaeology in establishing the ages of sediments and archaeological artefacts (Liritzis
We hope that the set of relevant measurement procedures, as well as the reference materials used by the GLL, presented in this manuscript can be an example of a reliable approach to ensuring the quality of the luminescence ages of tested materials.
Currently, the GLL is one of two laboratories working in Poland (the second one is in Toruń). Historically, luminescence dating has been applied in Poland for approximately 40 years, primarily for quaternary research, and was first implemented by M. Prószyński at Warsaw University, J. Butrym at Lublin University, S. Fedorowicz at Gdańsk University and M. Pazdur and A. Bluszcz at the Silesian University of Technology in Gliwice. The first studies addressing key aspects related to using luminescence phenomena for dating were published in Poland in the 1980s (Fedorowicz and Olszak, 1985; Bluszcz, 1986a, 1986b; Pazdur and Bluszcz, 1987a, 1987b). Currently, the laboratory is equipped with four optically stimulated luminescence dating (OSL/TL) readers (two Daybreak and two Risø), three high-resolution gamma spectrometers and three μDose systems (Tudyka
Establishing the luminescence age of a geological deposit requires different types of measurements. At this point, it is necessary to illustrate which laboratory measurements and procedures are vital to obtaining a single result of dating and how it can be interpreted. As a part of scientific services, the GLL focuses on quartz dating for all fractions using the standard SAR protocol (Murray and Wintle, 2000). Therefore, the following description is not relevant to feldspar dating.
Starting from the simplest definition, only two parameters are needed to determine the age of the sample to be tested: Age =
One of the most important aspects of luminescence dating is the appropriate selection of places from which the research material is collected. The elimination of undesirable uncertainties begins with the correct and appropriate sampling (Moska, 2019). The GLL uses the Eijkelkamp system (Eijkelkamp Soil & Water factory; Fig. 1), which meets all the requirements needed for the collection of samples for luminescence dating, including the ease of use, combined with the possibility of collecting material from even very hard sediments, and total protection of the collected material from sunlight. The Eijkelkamp system's design and the materials used to ensure safe and comfortable sampling conditions. This is particularly true about the protection of the collected material from sunlight, which is crucial for the OSL method.
Using metal pipes for collecting material from locations like walls or foundations of ancient buildings is not practicable. In such cases, a cordless drill with a coring saw (Fig. 2) is used to collect cores with a length and diameter of up to several centimetres.
The determination of the dose rate is an integral part of luminescence dating, and the precision and final quality of the luminescence dating strongly depend on the precision and accuracy of the dose rate determination. The ionising radiation to which minerals in the sediment are subjected originates from the 238U, 235U and 232Th decay chains and the isotopes 40K and 87Rb, which are natural components of the Earth's crust; in addition, there are secondary cosmic rays. At present, several methods are used for determining the concentration of the listed radioisotopes to assess dose rate in luminescence laboratories worldwide (Cunningham
In the GLL, three radiometric methods are used for determining 238U, 232Th decay chains and 40K content, namely, high-resolution gamma spectrometry, portable gamma scintillation spectrometry and the μDose system.
In the case of typical geological materials, one may safely assume the availability of several hundred grams of raw material for laboratory research purposes, and high-resolution gamma spectrometry can be an option. However, with regard to archaeologically significant material or sediment cores, a large amount of material is rarely available, and only a small piece of ceramics or a few grams of the cored material may be available for testing. In such cases, the μDose system measurement can be preferred as <3 g of material is needed. An additional problem with testing archaeological samples is possibly a significant difference in the content of natural radioisotopes in the samples themselves in comparison to the surrounding sediment. This can be, for example, assessed using portable gamma scintillation spectrometry. The Laboratory, therefore, always asks for additional samples of surrounding sediment in the case of dating archaeological artefacts as well as for geological, geographical and geometrical descriptions.
Depending on the material type, the most suitable dose rate measurement procedure is selected. In general, the dose rate is modelled considering the best available knowledge and estimates with:
In the GLL, routinely high-resolution gamma spectrometry is used to determine the dose rate, which is calculated from radioisotope specific activities. GLL is equipped with three high-resolution gamma spectrometers (Fig. 3). Two of them are equipped with extended-range germanium detectors. In one case, the germanium detector is characterised mathematically and equipped with the LabSOCS Software (CANBERRA). The gamma spectrometers are additionally used for measurements of the activity of various isotopes for environmental and technical studies (Poręba and Bluszcz, 2007; Tylman
For a typical geological sample, the material is dried and then crushed in a designated room. Next, samples are stored for at least 21 days in the gBEAKERs where >98% of 222Rn emanating from the sample remains in the container (Poręba
To analyse radioactivity in samples, the gamma spectrometer must be calibrated. There are a few methods to achieve this. The two main methods are either through obtaining the total efficiency curve or through the measurement of standards prepared in the same geometry and matrix composition as the investigated samples (Murray
Gamma spectrum lines used to determine the activities of the respective isotopes.
238U decay chain | 234Th | 63.3 keV |
234mPa | 1001.0 keV | |
226Ra | 186.2 keV | |
214Pb | 295.2 keV | |
351.9 keV | ||
214Bi | 609.3 keV | |
768.4 keV | ||
1120.3 keV | ||
1764.5 keV | ||
210Pb | 46.5 keV | |
232Th decay chain | 228Ac | 911.2 keV |
969.0 keV | ||
212Bi | 238.6 keV | |
208Tl | 583.6 keV | |
2614.5 keV | ||
40K | 40K | 1460.8 keV |
An example of a gamma-ray spectrum obtained for soil samples is presented in Fig. 5. In the case of low gamma radiation energy, to obtain the correct value of activity in the sample, self-absorption is included in the sample according to Cutshall
A wide group of samples are beyond the capabilities of conventional gamma spectrometry. This includes but is not limited to: very small archaeological samples, samples where secular equilibrium in decay chains is questionable or internal alpha- and beta-dose rates (Jacobs, 2004; Jacobs
The μDose system is calibrated with the materials as described in Section 2.2.1. Besides, a reference material mix of IAEA-RGU-1, IAEA-RGTh-1 and IAEA-RGK-1 in proportion 1:1:1 allows improving the calibration. Dose rate and its uncertainties are calculated using Monte Carlo and Bayesian methods that are built in the μDose software and take into account correlations (Tudyka
The μDose system requires different sample preparation than the HPGe detectors. First, the material must be dried. Then, 5–20 g of sample is ground down to particles of a diameter of approximately 20 microns using a planetary ball mill. This size is achieved in approximately 45 min at 200 rounds/min. Next, the sample must be placed on the 70 mm diameter sample holder, as shown in Fig. 7. Using a precise balance, exactly 3 g of the dried and ground sample is placed in the holder.
Next, the sample is placed inside the scintillation head, as shown in Fig. 8. Special attention should be paid to flattening the sample uniformly. The system calibration is only valid for the same mass of samples and standards.
The GLL is equipped with a portable gamma scintillation spectrometer manufactured by Canberra (InSpector 1000; Fig. 9) for recording energy spectra of ionising radiation in field conditions at sampling points for OSL dating. The portable gamma spectrometer was calibrated using concrete calibration blocks located at Oxford University (Rhodes and Schwenninger, 2007). Usually, four calibration blocks are used. Three of them are enriched with U, Th and K, and the fourth is pure concrete for background measurement. The portable gamma spectrometer is very useful for dose rate measurements, in particular, in areas where luminescence samples were taken from places close to boundaries between various sedimentation structures (which sometimes show a significantly different concentration of natural isotopes in the material – the best examples are sandy inserts in loess areas).
In addition to the measurement in the spectrometer, it is necessary to determine the dose derived from cosmic rays, which is calculated using the equation proposed by Prescott and Hutton (1994).
Water content is one of the basic parameters necessary for the correct determination of the dose rate. It has to be borne in mind that it is not, in this case, the current water content of the sediment, but a value that is representative for the entire deposition period. Based on the experience of the Gliwice luminescence dating laboratory, we approach this subject as follows:
In the laboratory, the current water content and full saturation level are measured for all sediments. Usually, there is no possibility to precisely estimate sediment changes of water content for the last several thousand years (we can only make some approximations). Based on the collected environmental information for each position and hydrological condition, we determine the most likely humidity range for a given place.
The entire procedure for preparing the material for OSL measurements takes place in a specialised chemical laboratory, where the only red light is used which does not remove the natural luminescence signal from the sample. The GLL's chemical preparation room is presented in Fig. 10. Throughout the entire laboratory, safe (deep golden amber) light is used. For this purpose, standard white fluorescent lamps equipped with an appropriate optical filter (professional certified foil – LEE filters 106 Primary Red) were used.
Chemical treatment is necessary to determine the equivalent dose from the extracted pure quartz. This procedure also depends on the size of the grains. Bigger grains (fraction above 45 μm up to 250 μm) at the beginning are treated with hydrochloric acid (20%) to remove carbonates, and hydrogen peroxide (20%) to remove organic matter. Between each step of the chemical preparation, the sample is rinsed with deionised water. A vital step in the extraction of pure quartz separation is the density separation using sodium polytungstate. During this procedure, a material with a density between 2.62 g/cm3 and 2.75 g/cm3 is separated. The used heavy liquid is collected, filtered and recycled.
The final step is treatment with concentrated (40%) hydrofluoric acid for 1 h to remove the outer layer of the grains of the quartz (about 10 μm) responsible for the absorption of alpha radiation dose (Aitken, 1985, 1998) and other minerals that might be still present after the density separation. After etching, to remove any residues, a final treatment with 20% hydrochloric acid for 20 min is used. For luminescence measurement, quartz extracted in this way is placed on steel disks with a diameter of 1 cm. Silicon oil (silkospray) as an adhesive is applied through a mask, and the grains are deposited as 6 mm circles.
The procedure for preparing fine grains (4–11 μm) is a little more complicated. The first step is to obtain the fraction below 45 μm using sieves. Next, as mentioned above, the material is treated with hydrochloric acid and hydrogen peroxide. Finally, the material is etched in concentrated hydrofluorosilicic acid (34%, H2SiF6) for a few days. Next, grains are washed minimum five times in distilled water, and hydrochloric acid (20%) is added for 2 h. After the following washes and drying, grains are ready for gravitational separation. The sample is placed in a tall test tube with acetone, in which gravitational separation of various fractions occurs over time. This procedure is repeated at least four times to obtain the most accurate material of the desired fraction. After obtaining the appropriate fraction, it is placed in 50 ml of acetone, and the well-mixed solution is then pipetted into 2 ml flat-bottomed tubes. Each tube, before filling, is fitted with an appropriately sized steel disc, on which, after the evaporation of the acetone, the fine grain fraction will be deposited. It is only then that the disks are ready for luminescence measurements.
A similar procedure is necessary to extract quartz from the archaeological artefacts such as ceramics. First, a mini drill with milling cutters is necessary to extract material from the sample (Fig. 11). The outer layer of about 1 mm is sanded off, and the powder is collected from under this layer. The next steps are similar to the procedure, which is associated with the preparation of fine grains.
The most representative fractions for different kinds of sample material can be divided into three main ranges:
4 ÷ 11 μm – ceramics, loess sediments 45 ÷ 63 μm – loess sediments, loess-like sediments 90 ÷ 200 μm – sandy material, dunes, deserts, bricks
The GLL is equipped with four luminescence readers, two Risø TL/OSL DA-20 (Fig. 12) and two Daybreak 2200 TL/OSL automatic luminescence readers (Fig. 13), all readers are equipped with a calibrated 90Sr/90Y beta source. For calibration, in GLL, we use Riso quartz (Hansen
For quartz samples, equivalent doses are determined using the SAR protocol (Murray and Wintle, 2000). The readers are fitted with appropriate optical filters enabling the correct detection of the luminescence signal emitted from the tested samples. Table 2 presents all the necessary steps of the measurement process.
Steps used in the protocol used for determining equivalent doses. For the quartz fraction, the SAR protocol (Murray and Wintle, 2000) is used.
1 | Irradiation with the regenerative beta dose |
2 | Preheat at a temperature of 260°C for 10 s (final temperature is determined after preheat plateau test) |
3 | Blue light stimulation at a temperature of 125°C for 100 s |
4 | Irradiation with the test dose |
5 | Cut-heat at a temperature of 220°C |
6 | Blue light stimulation at a temperature of 125°C for 100 s |
OSL, optically stimulated luminescence dating; SAR, single-aliquot regenerative-dose.
In our laboratory, we focus on single aliquot measurements because single grain measurements are much more time-consuming and the final results do not necessarily deliver more information and precision (Thomsen
The result of the registration of natural luminescence for the sample and luminescence obtained following irradiation with specific regeneration doses is the so-called growth curve, an example of which is shown in Fig. 14 (a typical luminescence decay curve is shown in addition to the growth curve). This curve is created by fitting the appropriate function (usually saturating exponential) to the measurement points. The final result of the analysis of each growth curve is the determination of the equivalent dose.
In the SAR protocol, sensitivity changes that may occur from one measurement cycle to another are measured by the OSL response to a small test dose (Murray and Wintle, 2000). The corrected OSL ratio (regenerated OSL response/OSL response to the fixed test dose) should be independent of prior dose or thermal treatment. This is tested by repeating a particular regenerative dose after various larger values have been used and comparing the ratio of the two regenerated sensitivity-corrected OSL responses (known as a recycling ratio); this ratio should ideally be close to unity from 0.90 to 1.10 (Murray and Wintle, 2000). Preheating the sample can also cause the recuperation of the OSL signal (Aitken, 1985). To test this, a 0 Gy regenerative dose step is incorporated into the SAR protocol (Murray and Wintle, 2000). The luminescence signal should then be zero (this is known as the ‘recuperated’ luminescence signal). Any grains for which this sensitivity-corrected recuperated signal is >5% of the corresponding natural signal are rejected (Murray and Olley, 2002).
The dose recovery test is a basic requirement to determine the suitability of the SAR protocol. This is performed by administering a known dose in the laboratory and estimating it by the same measurement procedure as for dose estimates (Wallinga
For each analysed sample, measurements are usually performed for approximately a dozen individual quartz aliquots leading to obtaining a dozen ages. To determine the final value of the equivalent dose, an appropriate statistical model needs to be implemented – in this case, the central age model (CAM, Galbraith
Within this model, the overdispersion parameter (σOD) is determined. This parameter refers to the distribution of the equivalent dose values for individual samples and increases with their dispersion. Fig. 15 shows two different examples of scattered values obtained for specific samples in the loess profile of Biały Kościół. These examples illustrate how different the results of individual quartz portions may appear in one sample. For the presented samples, the parameters of overdispersion were 5% and 19%, respectively. According to Galbraith's recommendations (Galbraith
The sources of overdispersion are both intrinsic (e.g. counting statistics as described by Adamiec
It is vital not to interpret the multimodal distributions of the probability density, an example of which is the distribution on the right in Fig. 15b, as the possibility of the occurrence of three independent measurement results (there are three maxima of probability density) because this would over-interpret the obtained results and could lead to the falsification of the real age of the sediment. Another crucial piece of information is that we can read from such distributions the efficiency of resetting the luminescent signal in the sample during its redeposition. Fig. 15a shows a graph for a sample of material originating from a typical loess period L1LL1 (lower upper loess, Markovic
In its many years of cooperation with the scientific community in Poland and abroad, the GLL has participated in the implementation of a large number of domestic and five international grants. This level, scale and scope of involvement in the scientific research community is reflected in the quantity and quality of its proprietary publications and success in obtaining grants for leading scientific research.
We would like to mention at least the most interesting research projects implemented in our laboratory:
‘CONTINENT’ – aimed specifically at the high-resolution reconstruction of climatic conditions in the Holocene and Eemian periods. Several gravitational cores have been obtained from different locations on the lake to provide undisturbed sediment records covering the last ca. 150 ky. Different climate parameters were deduced from organic and inorganic components of the deposited matter. OSL was one of the methods selected for dating Lake Baikal sediments and the only capable of providing absolute dating for the Eemian period (Moska ‘Establishing the chronology of loess formation in the eastern and western area of Poland based on high-resolution luminescence dating and geochemical studies of selected loess profiles’ – this project involves the use of a broad spectrum of research aimed at the reconstruction of climate change on the Polish territory over the past 100,000 years by a detailed analysis of the loess material. When we talk about climate change for the whole southern part of the Polish territory (because loess deposits occur only there), we have to remember that from the geographical location and the size of our country, the climate is not uniform. In the western part of Poland, there is a clear influence of the Atlantic climate, while in the east, continental climate, therefore considering the whole climate change in the past. In this project, a stratigraphic model was based on the detailed geochemical studies and a series of luminescence and radiocarbon dates (Moska and Bluszcz, 2013; Moska ‘Settlement strategies, activity and natural environment of late Palaeolithic hunters and gatherers in the upland areas of Poland on the example of a complex of sites from Sowin’ – the project aimed at presenting a new chronological model of recolonisation of the area located north of the Sudetes and the Carpathians after Late Glacial Maximum (LGM). Until recently, it was believed that reoccupation of these areas occurred only due to the Magdalenian people. New chronometric data (radiocarbon and OSL) coming from the Magdalenian and Epigravettian sites together with the application of Bayesian modelling of new records allowed us to present another model. In light of the new research, it seems that the groups classified as the Epigravettian and Magdalenian could have coexisted at the same time in the same areas. The new model states that it is possible to assume the coexistence of these groups during 2300 years between 16,500 and 14,200 years BP (Wiśniewski ‘The disappearance of the last ice sheet in the area of the South Baltic in the light of dating reservoir sediments and glaciofluvial deltas by the OSL method’ – the project was concentrating on detailed OSL dating of glaciofluvial and ice-marginal lake sediments which occur between end moraines of the Słupsk Bank and the Polish coast. As a result of research, it was found that the mode of deglaciation was changed from frontal and aerial into subaqueous after ice sheet retreat from moraines of the Gardno Phase. A large ice-marginal lake existed in front of the ice sheet during the Słupsk Bank Phase. The average age of this episode during the last SIS decay is 15.2 ka with a standard deviation of 0.9 ka (Uścinowicz ‘Life and death of an extinct rhino ( ‘Palaeohydrological Extreme Events – evidence and archives’. Results from geomorphological, sedimentological and geochronological analyses, together with micropaleontological and mineralogical characteristics of lacustrine deposits in five locations within the Kuray intermountain depression, southeast Altai, mountains of south Siberia, support the thesis of repeated formations of ice-dammed lakes during MIS-2 and their draining by high energy floods. OSL was a tool to estimate the age for a sandy layer at the top of the diluvial (i.e. large flood) deposit, revealed in a sedimentary sequence of the 1570 m a.s.l. strandline – one of the lowest preserved strandlines in the western part of the basin (Agatova ‘Multi-proxy investigations and advanced methods of determining numerical time scales in the reconstruction of the evolution of inland dunes in Poland during the end of the last glaciation’ – the European Sand Belt in Central Europe is a remarkably interesting remnant of the Late Glacial aeolian accumulation. Due to its geographical location, Poland is characterised by very extensive dune areas, which occur both in the form of large dune fields and individual dunes scattered throughout the country. Over the last several dozen years, a large volume of scientific and research works have been devoted to the stratigraphy of inland dunes. From the very beginning of the stratigraphic research, some researchers posed a purely hypothetical theory of the earlier formation of aeolian covers in the southern part of Poland (Older Dryas), and only later in the northern part (Younger Dryas). The assumptions of the hypothesis seem to be justified since the width of the European Sand Belt in Poland exceeds 400 km in the north-south direction. Such a significant spread should allow the detection of the temporal and spatial gradient of the formation of the aeolian units (Moska ‘Correlation of primary loess profiles of Ukraine based on luminescence dates’ – project concerns the issues of age as well as lithological and pedological development of the loess-paleosol sequences selected after direct research at the Ukraine territory. For dozens of identical samples, a series of comparative laboratory TL dating (Gdansk, Lublin) and OSL dating (Gliwice) was performed, supplemented by the physicochemical and micromorphological analyses (Fedorowicz ‘Reconstruction of anthropogenic soil erosion on selected loess areas of South Poland by dendrochronology, luminescence and isotopes (Cs-137, Pb-210) methods’. In this project, the SAR-OSL dating method was used to study Holocene water slope sediment simultaneously witch 137Cs and 210Pbex to measure modern soil erosion and sediment accumulation. Research carried out as part of this project clearly showed that loess accumulated loops in the Holocene are suitable for dating by the OSL method and showed why, despite the relatively short transport over the slope, good bleaching of the original luminescent signal is observed (Poręba ‘μDose’ this project was aimed to commercialise innovative solutions for measuring dose rate. The developed system, with the same name, offers new possibilities of increasing dose rate precession and accuracy, measuring small samples including internal dose rate and obtaining reliable dose rate in samples where secular equilibrium is uncertain. A detailed description of the μDose system is given in Section 2.2.2 and by Miłosz
All final luminescence results are placed in the Laboratory's GdTL luminescence database, and a final report from luminescence dating is generated. A part of the report, all necessary information connected with a particular example is presented (see Fig. 16). Each dated sample has its unique laboratory number (e.g. GdTL-4000). The report also contains information about the fraction used for measurement, water content, the activity of 238U and 232Th decay chains and 40K, effective dose rates, the number of investigated aliquots and regenerative doses used in luminescence protocol. Finally, the equivalent dose and the age of the sample are presented.
The age limit of luminescence generally ranges from several years to hundreds of thousands of years. Among all dated samples in the history in the GLL, the youngest dated sample was colluvium loess from Szyszczyce (Poręba