# Optically stimulated luminescence thermochronometry

Optically stimulated luminescence (OSL) thermochronometry is a dating method used to determine the time since quartz and/or feldspar began to store charge as it cools through the effective closure temperature. The closure temperature for quartz and Na-rich K-feldspar is 30-35 °C and 25 °C respectively. When quartz and feldspar are beneath the earth, they are hot. They cool when any geological process e.g. focused erosion causes their exhumation to the earth surface. As they cool, they trap electron charges originating from within the crystal lattice. These charges are accommodated within crystallographic defects or vacancies in their crystal lattices as the mineral cools below the closure temperature.
During detrapping of these electrons, luminescence is produced. The luminescence or light emission from the mineral is assumed to be proportional to the trapped electron charge population. The age recorded in standard OSL method is determined by counting the number of these trapped charges in an OSL detection system. The OSL age is the cooling age of the quartz and/or feldspar. This cooling history is a record of the mineral's thermal history, which is used to reconstruct the geological event.
The sub-Quaternary period (10$^4$ to 10$^5$ years) is the geological age where OSL is a favourable dating technique because of low closure temperature of quartz and feldspar used in this technique. The Quaternary period is marked by intense crustal erosion particularly within active mountain ranges, leading to high exhumation rate of crustal rocks and formation of sub-Quaternary sediments. Previous techniques (e.g. Apatite Fission Track, Zircon Fission Track, and (Uranium-Thorium)/ Helium dating) could not adequately track the geological age records particularly in the last ~300 thousand years. OSL dating is currently the only dating method that has been successfully applied to understand the cooling ages of the geological events.

## Theoretical concepts of electron trapping and detrapping for OSL measurement

In natural environment, crystal lattices of quartz and/or feldspar are bombarded with radiation released from radiogenic source such as in -situ radioactive decay. As the crystals are irradiated, charges are stored up in their crystallographic defects. The charge trapping process involves atomic-scale ionic substitution of both electron and hole within the crystal lattices of quartz and feldspar. The electron diffusion happens in response to ionizing radiation as the minerals cools below their closure temperature.
If quartz or feldspar gains are exposed to natural light source such as the sun, the trapped charges will be evicted in form of luminescence. This natural process is called bleaching. Any other process that could heat up the sample will also cause the trapped electrons to escape from the crystal lattice known as thermal bleaching. Optical bleaching of the mineral leads to eviction of trapped charges in the minerals, hence, careful sampling and handling must be followed to avoid using bleached sample for OSL thermochronometry. To artificially produce luminescence in the laboratory for luminescence study of the mineral, these two processes are adopted.

## Kinetic or rate equations for trapping and detrapping processes

A wide range of kinetic models have been developed to explain trapping and detrapping processes in quartz and feldspar crystals. Two of these models are particularly useful in determining the cooling histories of quartz or feldspar These models are known as the general order kinetic model and band tail model. The two models consider three major processes to characterize the mineral luminescence, which are: trapping process, thermal detrapping process and athermal detrapping process. Each of the processes are guided by different equations discussed below. These models are useful for the determining of cooling history of the mineral, which involves subtracting the differential sum of thermal detrapping and athermal detrapping from the trapping process (i.e. Trapping – (Thermal detrapping + Athermal detrapping).

### Rate equations

| Trapping process The process describes the rate at which a mineral acquires charge from its natural environment. The process is defined by the equation below. ${\frac {d{\tilde {n}}}{dt}}={\frac {D_{R}}{D_{o}}}(1-{\tilde {n}})^{\alpha }$ where ${\tilde {n}}$ is the ratio of electrons that are trapped (n) and the full storage capacity (N) in the mineral crystal. $D_{R}$ is the charge stored per unit time, known as dose rate (unit is Gy/ky). $D_{o}$ is the dose that correspond to the maximum storage capacity of the mineral. The probability that the mineral will get filled defines the ${\dot {D}}$ and $D_{o}$ relationship. $\alpha$ is an exponential factor greater than or equal to 1. | Thermal detrapping The equation describes how charges are evicted in form of luminescence from the mineral system when the mineral absorbed heat energy. This is where the general order kinetic model and the band tail model equations differs. 1. For the general order kinetic Model 2. For the band tail Model ${\frac {d{\tilde {n}}}{dt}}=-s{\tilde {n}}^{\beta }\exp \left({\frac {-E}{kT}}\right)$ ${\frac {d{\tilde {n}}}{dt}}=-s{\tilde {n}}^{\beta }\exp \left({\frac {-(E_{t}-E_{b})}{kT}}\right)$ Where $s$ frequency factor measured in per second (s$^{−1}$), $\beta$ is the kinetic order, which is usually less than or equal to 1; $E$ is the activation energy measured in (eV), $E_{t}$ and $E_{b}$ are the activation energies at the conduction band and the valence band respectively in the crystal lattices of the mineral. $k$ is the Boltzmann constant (eV/K) and $T$ is the absolute temperature (K). | Athermal detrapping This equation describes the luminescence recorded as the signal fades out known as 'anomalous fading'. This is determined without heating the mineral during fading experiment The relationship is defined by this equation. ${\frac {d{\tilde {n}}}{dt}}=-s{\tilde {n}}^{\beta }\exp \left(-p'^{-{\tfrac {1}{3}}}r'\right)$ where $p'$ represents the density of the crystal's recombination centers and r is the distance between the recombination centers and the electron traps. |

### Determination of cooling history from the kinetic equations

By combining the four equations above, a single differential equation is developed to convert the luminescence into cooling rate. We have:
${\frac {d{\tilde {n}}}{dt}}={\frac {D_{R}}{D_{o}}}(1-{\tilde {n}})^{\alpha }-s{\tilde {n}}^{\beta }\exp \left({\frac {-E}{kT}}\right)-s{\tilde {n}}^{\beta }\exp \left(-p'^{-{\tfrac {1}{3}}}r'\right)$ for the general order kinetic model; and
${\frac {d{\tilde {n}}}{dt}}={\frac {D_{R}}{D_{o}}}(1-{\tilde {n}})^{\alpha }-s{\tilde {n}}^{\beta }\exp \left({\frac {-(E_{t}-E_{b})}{kT}}\right)-s{\tilde {n}}^{\beta }\exp \left(-p'^{-{\tfrac {1}{3}}}r'\right)$ for the band tail model.
Any of the models can be used because the same series of laboratory experiments are followed for the estimation of all the parameters involved in the equations. The inversion of measured ${\tilde {n}}$ for a range of temperature -time history or T-t path can be used to determine the cooling rate. Sufficient number of T-t paths conducted in the laboratory is used to build a probability density function, which will help to determine the most likely cooling histories undergone by the mineral.

## Sample preparation

Bedrock samples from earth surface or boreholes are required earth materials for OSL dating. Minerals (quartz and/or feldspar) are usually separated from the rock or sediment samples under regulated laboratory lighting system similar to procedures used in archaeological OSL dating. The light source is usually a controlled red light condition to avoid luminescence signal resetting.$^$ Crushing of sample are gently carried out to avoid generating heat that is strong enough to reset OSL signal in the minerals. Crushed samples are separated by means of a sieve to get fine-grained. A range of values varying from 90 – 125 microns, 100 – 200 microns and 180 – 212 microns can be used for OSL measurement. The selected grains are chemically treated with HCl to digest carbonates and with H$_2$O$_2$ to remove organic materials that can contaminate the sensitivity of OSL signal during measurement. Feldspar and quartz with densities of < 2.62 g cm$^{−3}$ and < 2.68 g cm$^{−3}$ respectively are separated from other heavier minerals by density separation. Inclusions of zircon, apatite and feldspar in quartz as well as alpha-particles irradiated grain edges that can contaminate OSL signal are removed by etching in hydrofluoric acid (HF).

## OSL signal detection system

OSL ages are commonly measured using an automated Ris$\varnothing$Thermal Luminescence Reader (e.g. TL-DA-20). It contains an internal beta-source (e.g. $^{90}$Sr/$^{90}$Y) with optical stimulation emitted through laser diodes (LEDs). The reader also has a detection filter for transmission of stimulated luminescence signals. During this measurement, the mineral grain (quartz or feldspar) is glued on a heater strip (stainless-steel discs) using adhesive (commonly silicone spray). The mineral grain is stimulated with the light source. This light is the series of light emitting diode. This bombardment stimulates the electrons, which are trapped and begin to recombine in the crystal. During this process, they give the OSL signal, which is collected or recorded in the ray sensitive photomultiplier tube. The photomultiplier tube converts all the incident photons (i.e. light) to electronic charge. This is the basic principle of how the luminescence (light) emission from the minerals under investigation is measured.

## OSL age determination

To determine the OSL age of the sample, the dose rate, ($D_{R}$) and the equivalent dose ($D_{E}$). A dose is the quantity of natural radiation or energy absorbed by a mineral. The dose rate is the effective radiation absorbed from naturally occurring ionizing source per unit time.
The age is calculated by determining the ratio of equivalent dose ($D_{E}$) and the dose rate ($D_{R}$) using the equation below.
$A={\frac {D_{E}}{D_{R}}}$
where $A$ is the age (yr), $D_{E}$ is measured in Gray (Gy). Note that 1 Gy is equivalent to 1 J.kg$^{−1}$ (Joule per kilogram) and $D_{R}$ is Gy year$^{−1}$

### Dose rate determination

For a single grain of mineral, the dose rate ($D_{R}$) can be determined by measuring the concentrations of uranium, potassium and thorium by direct mass spectrometric analysis of quartz or feldspar grains. Ge-Gamma, INAA, X-ray flourescnce and ICP-MS or ICP-OES are spectrometers that can could be used. Other methods for the determination the dose rate include: (1) overburden cosmic dose rate estimation, (2) water content attenuation method, and (3) disequilibrium dose rate correction method. An average dose rate is usually calculated as representative of the dose rate.

### Equivalent dose determination

The equivalent dose ($D_{E}$) is also known as the dose response is determined from the dose response curve (see Plot B). The single-aliquot regenerative (SAR) protocol is a commonly used method for the determination of the equivalent dose. The protocol involves series of laboratory measurement of OSL signal (see Plot A), which is emitted by the aliquot after it has been optically stimulated at a known beta dose within a given time in seconds. The beta-source may be $^{90}$Sr/$^{90}$Y in an automated Ris$\varnothing$Thermal Luminescence Reader . During SAR protocol, the difference in the measurement for quartz and feldspar is mainly on the degree of heat required per time and the source of stimulation.
The first stage involves determination of the natural dose (see Plot B) preheating the aliquot to about 160 -130 °C (for feldspar) for 10 s or 160-300 °C (for quartz) when the natural luminescence signal (i.e. natural dose) is still intact. This is done to remove unstable signals in the mineral. After preheating, the aliquot is optically stimulated by Infrared light emitting diode (for feldspar) or Blue light emitting diode (for quartz) depending on which mineral (see OSL detection system) for 40 s at 125 °C (for feldspar) or 100 s at 125 °C (for quartz) and the natural OSL signal (N$_L$) is measured and recorded in the photomultiplier tube. For the second stage, the aliquot is irradiated with a fixed known test dose (beta dose). The aliquot is preheated at temperature less than 160 °C. The IRSL signal measurement is taken as a test dose IRSL response (N$_T$) after it has been optically stimulated for 40 s at 125 °C (for feldspar) or 100 s at 125 °C (for quartz). At this stage the aliquot is completely bleached. A regenerative test dose is then started after bleaching.
The same procedure as described above is followed but a range of regenerative dose is given at different temperature for sensitivity correction of OSL signal (See Plot B). For the regenerative dose measurement, the aliquot is irradiated with a known dose before preheating at 160-130 °C for 10 s or 160-300 °C for feldspar or quartz respectively while the signal response (R$_i$) is measured. A fixed test dose is by irradiating the aliquot and a preheating of the aliquot is carried out at a temperature less than 160 °C. The aliquot is optically stimulated at the same rate and the IRSL signal (R$_T$) is measured. The steps are repeated for range of different regenerative dose including zero test dose. During each of the tests, all OSL signals are recorded in the photomultiplier tube and the OSL counts are plotted against the OSL exposure time in seconds as shown in the OSL signal curve (first graph).
For sensitivity correction, N$_L$ is plotted against N$_T$ representing the natural OSL signal while the plot of R$_i$ against R$_T$ representing regenerative dose test (see Plot B). The natural dose is along the vertical axis because no laboratory dose is given at the stage. The regenerative dose measurement will vary with respect to the given dose at each stage. The equivalent dose ($D_{E}$) is determined by drawing a line (red discontinuous line in Plot B) from the natural dose to intercept with the regenerative dose curve. The point of interception with the curve represent the equivalent dose by reading its value on the horizontal axis (See Plot B). The corresponding dose value at the horizontal axis is recorded for the equivalent dose ($D_{E}$).

## Applications.
### General applications

OSL finds application in all low-temperature (<50 °C) tectonic and sedimentary processes. These studies are mainly captured within the sub-Quaternary period including, but not limited to focused fluvial and/or glacial erosion, rock exhumation and evolution of topography in active tectonic regions. Other applications include glaciation deposits, lagoon deposits, storm surge and tsunami deposits, lake deposits including shoreline migration history, fluvial erosion deposits, loess deposit records. For example, the rate of slip on a normal faults plane can also be modelled, the rate of glacial or fluvial erosion of the earth surface can also be modelled as well as when sedimentary deposits are found within the sub-Quaternary period.
In active tectonics regions, the application of OSL dating is very useful in tracking the thermal history and rate of rock exhumation towards the Earth's surface. The closer the cooling ages, the higher the rate of erosion and/or exhumation of the rock unit under investigation. When the OSL age of quartz or feldspar is known, the obtained ages are coupled with the existing thermal-mechanical equations e.g. Pecube to reconstruct the thermal-mechanical history.
The OSL ages (see diagram), cooling ages, elevation data are plotted against the horizontal distance where samples and elevation data were collected to interpret the exhumation rate of rock or the evolution of the relief system through time. For example, OSL dating has been applied in determined the cooling histories of some rapidly eroding active regions at sub-Quaternary time-scale (i.e. 10$^4$ to 10$^5$ years). These examples are Whataroa-Perth catchment area in the Southern Alps of New Zealand and Namche Barwa-Gyala Peri dome in eastern Himalaya. In the Namche Barwa-Gyala Peri dome, river erosion was prevalent while glacial erosion was the main active process in the Whataroa-Perth catchment area. In both studies, the rate of exhumation and evolution of the relief systems were estimated by inversion of OSL thermochronological ages.