Calculation of radiocarbon dates

Last updated

The calculation of radiocarbon dates determines the age of an object containing organic material by using the properties of radiocarbon (also known as carbon-14), a radioactive isotope of carbon.

Contents

Radiocarbon dating methods produce data based on the ratios of different carbon isotopes in a sample that must then be further manipulated in order to calculate a resulting "radiocarbon age". Radiocarbon dating is also referred to as carbon dating or carbon-14 dating. Calculations of radiocarbon dates are typically made based on measurements from beta counting devices or from accelerator mass spectrometers (AMS). There are several possible sources of error in both the beta counting and AMS methods.

Calculations

The calculations to be performed depend on the measurements taken based on the technology used, since beta counters measure the sample's radioactivity, whereas accelerator mass spectrometers (AMS) determine the ratio of the three different carbon isotopes in the sample. [1]

Standards

The calculations to convert measured data to an estimate of the age of the sample require the use of several standards. One of these, the standard for normalizing δ13C values, is Pee Dee Belemnite (PDB), a fossil which has a 13
C
/12
C
ratio of 1.12372%. [2] A related standard is the use of wood, which has a δ13C of -25‰, as the material for which radiocarbon ages are calibrated. Since different materials have different δ13C values, it is possible for two samples of different materials, of the same age, to have different levels of radioactivity and different 14
C
/12
C
ratios. To compensate for this, the measurements are converted to the activity, or isotope ratio, that would have been measured if the sample had been made of wood. This is possible because the δ13C of wood is known, and the δ13C of the sample material can be measured, or taken from a table of typical values. The details of the calculations for beta counting and AMS are given below. [3]

Another standard is the use of 1950 as "present", in the sense that a calculation that shows that a sample's likely age is 500 years "before present" means that it is likely to have come from about the year 1450. This convention is necessary in order to keep published radiocarbon results comparable to each other; without this convention, a given radiocarbon result would be of no use unless the year it was measured was also known—an age of 500 years published in 2010 would indicate a likely sample date of 1510, for example. In order to allow measurements to be converted to the 1950 baseline, a standard activity level is defined for the radioactivity of wood in 1950. Because of the fossil fuel effect, this is not actually the activity level of wood from 1950; the activity would have been somewhat lower. [4] The fossil fuel effect was eliminated from the standard value by measuring wood from 1890, and using the radioactive decay equations to determine what the activity would have been at the year of growth. The resulting standard value, Aabs, is 226 becquerels per kilogram of carbon. [5]

Both beta counting and AMS measure standard samples as part of their methodology. These samples contain carbon of a known activity. [6] The first standard, Oxalic Acid SRM 4990B, also referred to as HOxI, was a 1,000 lb batch of oxalic acid created in 1955 by the National Institute of Standards and Technology (NIST). Since it was created after the start of atomic testing, it incorporates bomb carbon, so measured activity is higher than the desired standard. This is addressed by defining the standard to be 0.95 times the activity of HOxI. [5]

All of this first standard has long since been consumed, and later standards have been created, each of which has a given ratio to the desired standard activity. A secondary standard, Oxalic Acid SRM 4990C, also referred to as HOxII, 1,000 lb of which was prepared by NIST in 1977 from French beet harvests, is now in wide use. [7]

Calculations for beta counting devices

To determine the age of a sample whose activity has been measured by beta counting, the ratio of its activity to the activity of the standard must be found. The equation: [8]

gives the required ratio, where As is the true activity of the sample, Astd is the true activity of the standard, Ms is the measured activity of the sample, Mstd is the measured activity of the standard, and Mb is the measured activity of the blank. [8]

A correction must also be made for fractionation. The fractionation correction converts the 14
C
/12
C
ratio for the sample to the ratio it would have had if the material was wood, which has a δ13C value of -25‰. This is necessary because determining the age of the sample requires a comparison of the amount of 14
C
in the sample with what it would have had if it newly formed from the biosphere. The standard used for modern carbon is wood, with a baseline date of 1950. [3]

Correcting for fractionation changes the activity measured in the sample to the activity it would have if it were wood of the same age as the sample. The calculation requires the definition of a 13
C
fractionation factor, which is defined for any sample material as [4]

The 14
C
fractionation factor, Frac14/12, is approximately the square of this, to an accuracy of 1‰: [4]

Multiplying the measured activity for the sample by the 14
C
fractionation factor converts it to the activity that it would have had had the sample been wood: [4]

where Asn is the normalized activity for the sample, and Frac14/12 (s) is the 14
C
fractionation factor for the sample. [4]

The equation for δ13C given earlier can be rearranged to [4]

Substituting this in the 14
C
fractionation factor, and also substituting the value for δ13C for wood of -25‰, gives the following expression: [4]

where the δ13C value remaining in the equation is the value for the sample itself. This can be measured directly, or simply looked up in a table of characteristic values for the type of sample materialthis latter approach leads to increased uncertainty in the result, as there is a range of possible δ13C values for each possible sample material. Cancelling the PDB 13
C
/12
C
ratio reduces this to: [4]

AMS calculations

The results from AMS testing are in the form of ratios of 12
C
, 13
C
, and 14
C
. These ratios are used to calculate Fm, the "fraction modern", defined as

where Rnorm is the 14
C
/12
C
ratio for the sample, after correcting for fractionation, and Rmodern is the standard 14
C
/12
C
ratio for modern carbon. [9]

The calculation begins by subtracting the ratio measured for the machine blank from the other sample measurements. That is:

where Rs is the measured sample 14
C
/12
C
ratio; Rstd is the measured ratio for the standard; Rpb is the measured ratio for the process blank, and Rmb is the measured ratio for the machine blank. The next step, to correct for fractionation, can be done using either the 14
C
/12
C
ratio or the 14
C
/13
C
ratio, and also depends on which of the two possible standards was measured: HOxI or HoxII. R'std is then R'HOxI or R'HOxII, depending on which standard was used. The four possible equations are as follows. First, if the 14
C
/12
C
ratio is used to perform the fractionation correction, the following two equations apply, one for each standard. [9]

If the 14
C
/13
C
ratio is used instead, then the equations for each standard are: [9]

The δ13C values in the equations measure the fractionation in the standards as CO
2
, prior to their conversion to graphite to use as a target in the spectrometer. This assumes that the conversion to graphite does not introduce significant additional fractionation. [9]

Once the appropriate value above has been calculated, Rmodern can be determined; it is [9]

The values 0.95 and 0.7459 are part of the definition of the two standards; they convert the 14
C
/12
C
ratio in the standards to the ratio that modern carbon would have had in 1950 if there had been no fossil fuel effect. [9]

Since it is common practice to measure the standards repeatedly during an AMS run, alternating the standard target with the sample being measured, there are multiple measurements available for the standard, and these measurements provide a couple of options in the calculation of Rmodern. Different labs use this data in different ways; some simply average the values, while others consider the measurements made on the standard target as a series, and interpolate the readings that would have been measured during the sample run, if the standard had been measured at that time instead. [9]

Next, the uncorrected fraction modern is calculated; "uncorrected" means that this intermediate value does not include the fractionation correction. [9]

Now the measured fraction modern can be determined, by correcting for fractionation. As above there are two equations, depending on whether the 14
C
/12
C
or 14
C
/13
C
ratio is being used. If the 14
C
/12
C
ratio is being used: [9]

If the 14
C
/13
C
ratio is being used: [9]

The δ13Cs value is from the sample itself, measured on CO
2
prepared while converting the sample to graphite. [9]

The final step is to adjust Fmms for the measured fraction modern of the process blank, Fmpb, which is calculated as above for the sample. One approach [note 1] is to determine the mass of the measured carbon, Cms, along with Cpb, the mass of the process blank, and Cs, the mass of the sample. The final fraction modern, Fms is then [9]

The fraction modern is then converted to an age in "radiocarbon years", meaning that the calculation uses Libby's half-life of 5,568 years, not the more accurate modern value of 5,730 years, and that no calibration has been done: [10]

Errors and reliability

There are several possible sources of error in both the beta counting and AMS methods, although laboratories vary in how they report errors. All laboratories report counting statistics—that is, statistics showing possible errors in counting the decay events or number of atoms—with an error term of 1σ (i.e. 68% confidence that the true value is within the given range). [11] These errors can be reduced by extending the counting duration: for example, testing a modern benzene sample will find about eight decay events per minute per gram of benzene, and 250 minutes of counting will suffice to give an error of ± 80 years, with 68% confidence. If the benzene sample contains carbon that is about 5,730 years old (the half-life of 14
C
), then there will only be half as many decay events per minute, but the same error term of 80 years could be obtained by doubling the counting time to 500 minutes. [12] [13] Note that the error term is not symmetric, though the effect is negligible for recent samples; for a sample with an estimated age of 30,600 years, the error term might be +1600 to -1300. [11]

To be completely accurate, the error term quoted for the reported radiocarbon age should incorporate counting errors not only from the sample, but also from counting decay events for the reference sample, and for blanks. It should also incorporate errors on every measurement taken as part of the dating method, including, for example, the δ13C term for the sample, or any laboratory conditions being corrected for such as temperature or voltage. These errors should then be mathematically combined to give an overall term for the error in the reported age, but in practice laboratories differ, not only in the terms they choose to include in their error calculations, but also in the way they combine errors. The resulting 1σ estimates have been shown to typically underestimate the true error, and it has even been suggested that doubling the given 1σ error term results in a more accurate value. [11] [12]

The usual presentation of a radiocarbon date, as a specific date plus or minus an error term, obscures the fact that the true age of the object being measured may lie outside the range of dates quoted. In 1970, the British Museum radiocarbon laboratory ran weekly measurements on the same sample for six months. The results varied widely (though consistently with a normal distribution of errors in the measurements), and included multiple date ranges (of 1σ confidence) that did not overlap with each other. The extreme measurements included one with a maximum age of under 4,400 years, and another with a minimum age of over 4,500 years. [14]

It is also possible for laboratories to have systematic errors, caused by weaknesses in their methodologies. For example, if 1% of the benzene in a modern reference sample is allowed to evaporate, scintillation counting will give a radiocarbon age that is too young by about 80 years. Laboratories work to detect these errors both by testing their own procedures, and by periodic inter-laboratory comparisons of a variety of different samples; any laboratories whose results differ from the consensus radiocarbon age by too great an amount may be suffering from systematic errors. Even if the systematic errors are not corrected, the laboratory can estimate the magnitude of the effect and include this in the published error estimates for their results. [15]

The limit of measurability is approximately eight half-lives, or about 45,000 years. Samples older than this will typically be reported as having an infinite age. Some techniques have been developed to extend the range of dating further into the past, including isotopic enrichment, or large samples and very high precision counters. These methods have in some cases increased the maximum age that can be reported for a sample to 60,000 and even 75,000 years. [16] [17]

Notes

  1. McNichol & Burr give two other calculations, one of which can be shown to be equivalent to the one given here. The other depends on the process blank being the same mass as the sample. [9]

Footnotes

  1. McNichol, Jull & Burr, "Converting AMS Data to Radiocarbon Values: Considerations and Conventions", pp. 313.
  2. Dass (2007), p.276.
  3. 1 2 Aitken, Science-based Dating in Archaeology, p. 92–95.
  4. 1 2 3 4 5 6 7 8 Eriksson Stenstrom et al. A Guide to Radiocarbon Units and Calculations, p. 6.
  5. 1 2 L'Annunziata, Radioactivity, p. 528.
  6. Aitken, Science-based Dating in Archaeology, pp. 82-85.
  7. J. Terasmae, "Radiocarbon Dating: Some Problems and Potential Developments", in Mahaney, Quaternary Dating Methods, p. 5.
  8. 1 2 Eriksson Stenström et al., "A guide to radiocarbon units and calculations", p. 3.
  9. 1 2 3 4 5 6 7 8 9 10 11 12 13 McNichol, Jull & Burr, "Converting AMS Data to Radiocarbon Values: Considerations and Conventions", pp. 315-318.
  10. "Radiocarbon Data Calculations: NOSAMS". Woods Hole Oceanographic Institution. 2007. Retrieved August 27, 2013.
  11. 1 2 3 Taylor, Radiocarbon Dating, p. 102−104.
  12. 1 2 Bowman, Radiocarbon Dating, pp. 3839.
  13. Taylor, Radiocarbon Dating, p. 124.
  14. Taylor, Radiocarbon Dating, pp. 125−126.
  15. Bowman, Radiocarbon Dating, pp. 40−41.
  16. Mike Walker, Department of Archaeology and Anthropology, University of Wales, Lampeter, UK Quaternary Dating Methods, John Wiley & Sons, Ltd. 2005, page 23. Online, pdf 9 MByte Archived 2014-07-14 at the Wayback Machine
  17. Shennan, Ian; Long, Antony; Horton, Benjamin (2014). Handbook of Sea-Level Research. John Wiley & Sons. p. 350. ISBN   978-1118452561 . Retrieved 31 October 2018.

Related Research Articles

In a chemical reaction, chemical equilibrium is the state in which both the reactants and products are present in concentrations which have no further tendency to change with time, so that there is no observable change in the properties of the system. This state results when the forward reaction proceeds at the same rate as the reverse reaction. The reaction rates of the forward and backward reactions are generally not zero, but they are equal. Thus, there are no net changes in the concentrations of the reactants and products. Such a state is known as dynamic equilibrium.

pH Measure of the level of acidity or basicity of an aqueous solution

In chemistry, pH, historically denoting "potential of hydrogen", is a scale used to specify the acidity or basicity of an aqueous solution. Acidic solutions are measured to have lower pH values than basic or alkaline solutions.

Radiocarbon dating is a method for determining the age of an object containing organic material by using the properties of radiocarbon, a radioactive isotope of carbon.

In chemistry, an acid dissociation constant is a quantitative measure of the strength of an acid in solution. It is the equilibrium constant for a chemical reaction

In electrochemistry, the Nernst equation is a chemical thermodynamical relationship that permits the calculation of the reduction potential of a reaction from the standard electrode potential, absolute temperature, the number of electrons involved in the redox reaction, and activities of the chemical species undergoing reduction and oxidation respectively. It was named after Walther Nernst, a German physical chemist who formulated the equation.

In chemistry, the molar mass of a chemical compound is defined as the ratio between the mass and the amount of substance of any sample of said compound. The molar mass is a bulk, not molecular, property of a substance. The molar mass is an average of many instances of the compound, which often vary in mass due to the presence of isotopes. Most commonly, the molar mass is computed from the standard atomic weights and is thus a terrestrial average and a function of the relative abundance of the isotopes of the constituent atoms on Earth. The molar mass is appropriate for converting between the mass of a substance and the amount of a substance for bulk quantities.

In physical chemistry, Henry's law is a gas law that states that the amount of dissolved gas in a liquid is directly proportional to its partial pressure above the liquid. The proportionality factor is called Henry's law constant. It was formulated by the English chemist William Henry, who studied the topic in the early 19th century.

In physical organic chemistry, a kinetic isotope effect (KIE) is the change in the reaction rate of a chemical reaction when one of the atoms in the reactants is replaced by one of its isotopes. Formally, it is the ratio of rate constants for the reactions involving the light (kL) and the heavy (kH) isotopically substituted reactants (isotopologues):

The equilibrium constant of a chemical reaction is the value of its reaction quotient at chemical equilibrium, a state approached by a dynamic chemical system after sufficient time has elapsed at which its composition has no measurable tendency towards further change. For a given set of reaction conditions, the equilibrium constant is independent of the initial analytical concentrations of the reactant and product species in the mixture. Thus, given the initial composition of a system, known equilibrium constant values can be used to determine the composition of the system at equilibrium. However, reaction parameters like temperature, solvent, and ionic strength may all influence the value of the equilibrium constant.

Isotope geochemistry is an aspect of geology based upon the study of natural variations in the relative abundances of isotopes of various elements. Variations in isotopic abundance are measured by isotope ratio mass spectrometry, and can reveal information about the ages and origins of rock, air or water bodies, or processes of mixing between them.

Carbon (6C) has 15 known isotopes, from 8
C
to 22
C
, of which 12
C
and 13
C
are stable. The longest-lived radioisotope is 14
C
, with a half-life of 5.70(3)×103 years. This is also the only carbon radioisotope found in nature—trace quantities are formed cosmogenically by the reaction 14
N
+
n
14
C
+ 1
H
. The most stable artificial radioisotope is 11
C
, which has a half-life of 20.3402(53) min. All other radioisotopes have half-lives under 20 seconds, most less than 200 milliseconds. The least stable isotope is 8
C
, with a half-life of 3.5(1.4)×10−21 s.

An isotopic signature is a ratio of non-radiogenic 'stable isotopes', stable radiogenic isotopes, or unstable radioactive isotopes of particular elements in an investigated material. The ratios of isotopes in a sample material are measured by isotope-ratio mass spectrometry against an isotopic reference material. This process is called isotope analysis.

The environmental isotopes are a subset of isotopes, both stable and radioactive, which are the object of isotope geochemistry. They are primarily used as tracers to see how things move around within the ocean-atmosphere system, within terrestrial biomes, within the Earth's surface, and between these broad domains.

In geochemistry, paleoclimatology and paleoceanography δ18O or delta-O-18 is a measure of the ratio of stable isotopes oxygen-18 (18O) and oxygen-16 (16O). It is commonly used as a measure of the temperature of precipitation, as a measure of groundwater/mineral interactions, and as an indicator of processes that show isotopic fractionation, like methanogenesis. In paleosciences, 18O:16O data from corals, foraminifera and ice cores are used as a proxy for temperature.

Equilibrium constants are determined in order to quantify chemical equilibria. When an equilibrium constant K is expressed as a concentration quotient,

<i>δ</i><sup>13</sup>C Measure of relative carbon-13 concentration in a sample

In geochemistry, paleoclimatology, and paleoceanography δ13C is an isotopic signature, a measure of the ratio of the two stable isotopes of carbon—13C and 12C—reported in parts per thousand. The measure is also widely used in archaeology for the reconstruction of past diets, particularly to see if marine foods or certain types of plants were consumed.

The variation in the 14
C
/12
C
ratio in different parts of the carbon exchange reservoir means that a straightforward calculation of the age of a sample based on the amount of 14
C
it contains will often give an incorrect result. There are several other possible sources of error that need to be considered. The errors are of four general types:

Hydrogen isotope biogeochemistry is the scientific study of biological, geological, and chemical processes in the environment using the distribution and relative abundance of hydrogen isotopes. There are two stable isotopes of hydrogen, protium 1H and deuterium 2H, which vary in relative abundance on the order of hundreds of permil. The ratio between these two species can be considered the hydrogen isotopic fingerprint of a substance. Understanding isotopic fingerprints and the sources of fractionation that lead to variation between them can be applied to address a diverse array of questions ranging from ecology and hydrology to geochemistry and paleoclimate reconstructions. Since specialized techniques are required to measure natural hydrogen isotope abundance ratios, the field of hydrogen isotope biogeochemistry provides uniquely specialized tools to more traditional fields like ecology and geochemistry.

δ<sup>34</sup>S Standardized measurement method

The δ34S value is a standardized method for reporting measurements of the ratio of two stable isotopes of sulfur, 34S:32S, in a sample against the equivalent ratio in a known reference standard. Presently, the most commonly used standard is Vienna-Canyon Diablo Troilite (VCDT). Results are reported as variations from the standard ratio in parts per thousand, per mil or per mille, using the ‰ symbol. Heavy and light sulfur isotopes fractionate at different rates and the resulting δ34S values, recorded in marine sulfate or sedimentary sulfides, have been studied and interpreted as records of the changing sulfur cycle throughout the earth's history.

Methane clumped isotopes are methane molecules that contain two or more rare isotopes. Methane (CH4) contains two elements, carbon and hydrogen, each of which has two stable isotopes. For carbon, 98.9% are in the form of carbon-12 (12C) and 1.1% are carbon-13 (13C); while for hydrogen, 99.99% are in the form of protium (1H) and 0.01% are deuterium (2H or D). Carbon-13 (13C) and deuterium (2H or D) are rare isotopes in methane molecules. The abundance of the clumped isotopes provides information independent from the traditional carbon or hydrogen isotope composition of methane molecules.

References