Documentation — 10Be / 26Al exposure age calculator

Version 3 is the current stable version or the exposure age and erosion age calculator.

This page summarizes some of the documentation from hess.ess.washington.edu,
e.g. from here, here and here.

Calculators home
Exposure ages
Erosion rates

Paper describing the calculation methods: PDF

This paper was published in Quaternary Geochronology (v.3, p. 174, 2008) and is intended to be the referenceable documentation for exposure ages and erosion rates calculated here. For basic information about the inner workings, goals, assumptions, and limitations of the calculators, read this first.

Update to version 3:

Excerpt from the official documentation

Significant changes have been made since the original release, and some parts of the original paper are no longer up-to-date. Updates reflecting these changes are described in this PDF (v2.2), in this PDF (v2.2.1) and here (v2.3).

Version 3 is similar in concept to the 2008 online exposure age calculators (that is, up to version 2.3) that are also available at hess.ess.washington.edu and described in a paper by Balco et al. (2008). The design concept for version 3 is to make the following improvements:

1. Do calculations not just for Be-10 and Al-26, but all commonly measured nuclides produced by spallation and muon interactions.

2. Remove obsolete production rate scaling methods and add more recent methods based on particle transport modeling, including (i) the "LSDn" scaling method of Nat Lifton, which is based on the work of Tatsuhiko Sato, and (ii) potentially, a similar method developed by David Argento.

3. Improve and update various ancillary computation schemes and input parameters, including, among others, muon interaction cross-sections, atmosphere models, and figure generation.

4. While doing all this, improve execution speed so that the code runs fast enough to serve as a back end for online databases that require dynamically calculated exposure ages (e.g., ICE-D:ANTARCTICA).

The difference between this and the system described by Marrero and others (2016) is that the Marrero code was designed to carry out the model-fitting experiments described in Borchers and others (2016). The aim of these experiments was to determine whether or not various scaling models could or could not be fit to production rate calibration data. Clearly for this purpose one must make sure that differences between scaling model predictions and observations are not just due to simplifying assumptions or numerical approximations; thus, their code includes a very comprehensive representation of the physics of nuclide production, does not make simplifying assumptions, and is designed to have very good numerical precision. While all these features are necessary for the purposes of the Borchers study, they make the code quite slow.

In contrast, the aim of the version 3 code described here is to make as many simplifying assumptions and numerical shortcuts as possible while retaining acceptable accuracy for the single application of exposure-dating of surface samples (and its reverse, i.e., estimating production rates from independently dated surface samples). "Acceptable accuracy" is approximately defined as 2-4 per mil, which is small compared to (i) measurement precision for cosmogenic nuclides (in nearly all cases) and (ii) uncertainty in production rate scaling methods. To summarize, the overall design goal for this code is to compute exposure ages and/or production rates as fast as possible while maintaining acceptable precision (and also without doing anything that is really badly unphysical).

The official documentation of "The online exposure age calculator formerly known as the CRONUS-Earth online exposure age calculator", version 3 can be found here.


Input format:

Excerpt from the official documentation

Calculators home
Exposure ages
Erosion rates

What the overall format looks like for non-Cl-36 data:

Data to be entered are arranged in blocks, or lines, that are separated by semicolons; so each string of text that ends in a semicolon describes either a sample or a nuclide concentration measurement made on that sample. Here is an example:

                PH-1	41.3567	-70.7348	91	std	4.5	2.65	1	0.00008	1999;
                PH-1	Be-10	quartz	123453	3717	07KNSTD;
                PH-1	Al-26	quartz	712408	31238	KNSTD;
                

This has three lines: one that describes properties of the sample (what the properties are is described in more detail below); one that describes a Be-10 measurement; and one that describes an Al-26 measurement.

The order of the lines doesn't matter, and the nuclide concentration lines don't need to appear right after the sample they pertain to — they are linked to the sample by the fact that the sample name is included in each of the nuclide concentration lines.

So this text block describing two samples with both Be-10 and Al-26 measurements:

                PH-1	41.3567	-70.7348	91	std	4.5	2.65	1	0.00008	1999;
                PH-1	Be-10	quartz	123453	3717	07KNSTD;
                PH-1	Al-26	quartz	712408	31238	KNSTD;
                WR-2	41.3937	-70.6991	54	std	2	2.65	1	0.00008	1999;
                WR-2	Be-10	quartz	122903	3584	07KNSTD;
                WR-2	Al-26	quartz	749102	42198	KNSTD;
                

Works the same as this block:

                PH-1	41.3567	-70.7348	91	std	4.5	2.65	1	0.00008	1999;
                WR-2	41.3937	-70.6991	54	std	2	2.65	1	0.00008	1999;
                PH-1	Be-10	quartz	123453	3717	07KNSTD;
                PH-1	Al-26	quartz	712408	31238	KNSTD;
                WR-2	Be-10	quartz	122903	3584	07KNSTD;
                WR-2	Al-26	quartz	749102	42198	KNSTD;
                

And the same as this block:

                PH-1	41.3567	-70.7348	91	std	4.5	2.65	1	0.00008	1999;
                WR-2	41.3937	-70.6991	54	std	2	2.65	1	0.00008	1999;
                PH-1	Be-10	quartz	123453	3717	07KNSTD;
                WR-2	Be-10	quartz	122903	3584	07KNSTD;
                PH-1	Al-26	quartz	712408	31238	KNSTD;
                WR-2	Al-26	quartz	749102	42198	KNSTD;
                

Although I've shown the data lines as separate lines of text, the separator is actually the semicolon, so this:

                PH-1 41.3567 -70.7348 91 std 4.5 2.65 1 0.00008 1999;
                PH-1 Be-10 quartz 123453 3717 07KNSTD;
                

Is the same as this, which may often be easier to lay out in spreadsheets:

                PH-1 41.3567 -70.73483333 91 std 4.5 2.65 1 0.00008 1999; PH-1 Be-10 quartz 123453 3717 07KNSTD;
                

Tedious and pedantic formatting details for each data line:


A. Sample data line. These have 10 fields, as follows. Here is an example:

                PH-1	41.3567	-70.7348	91	std	4.5	2.65	1	0.00008	1999;
                

A.1. Sample name. Same as version 2 - a text string not exceeding 24 characters. Can include letters, numbers, and dashes (hyphens). May not contain white space or anything that even vaguely resembles an escape character, e.g., slashes of either direction, commas, quotes, colons, or semicolons.

A.2. Latitude. Same as version 2. Decimal degrees: north positive, south negative.

A.3. Longitude. Same as version 2. Decimal degrees: east positive, west negative.

A.4. Elevation/pressure. Same as version 2: meters or hPa, depending on selection below.

A.5. Elevation/pressure handling flag. Same as version 2; this is a three-letter text string. If you have supplied elevations in meters and the default elevation/pressure relationship (currently a spatially variable scheme based on the ERA-40 reanalysis, courtesy of Nat Lifton) is applicable at your site, enter "std" here. If you have supplied elevations in meters, your site is in Antarctica, and you want to use an Antarctic-specific elevation/pressure relationship, enter "ant" here. If you have entered pressure in hPa, enter "pre" here. Any text other than these three options will cause an error.

A.6. Sample thickness. Same as version 2. Centimeters.

A.7. Sample density. Same as version 2. Grams per cubic centimeter.

A.8. Shielding correction. Same as version 2; a shielding factor between 0 and 1.

A.9. Erosion rate inferred from independent evidence. Same as version 2. Centimeters per year.

A.10. Date of sample collection. This is not in version 2. The idea here is so that paleomagnetic reconstructions for the last couple of hundred years are correctly aligned with the date the sample was collected. Of course, this issue is totally irrelevant for samples that are more than a few hundred years old, and it's probably irrelevant for young samples too. In any case, you only need to worry about this at all if you have really young samples. Thus, you can also enter zero and a default date of 2010 will be assumed.

A.11. Then the line ends with a semicolon.



B. Nuclide concentration lines.

In general, all the lines describing nuclide concentration measurements have information in the following order: the name of the corresponding sample; the nuclide measured; the mineral in which it is measured; the nuclide concentration and uncertainty; and then standardization information, which varies by nuclide, at the end. The currently available set of nuclide-mineral pairs is as follows:

                He-3	quartz
                He-3	olivine
                He-3	pyroxene
                Be-10	quartz
                C-14	quartz
                Al-26	quartz
                Ne-21	quartz
                

Exotica such as Ne-21 in sanidine or Be-10 in olivine are not, at present, supported. The following sections detail the format of the data entry lines for Be-10 and Al-26 in quartz; refer to the official documentation for other nuclide-mineral pairs listed above.

B.1. Be-10-in-quartz and Al-26-in-quartz measurement lines. These have 6 fields, as follows. Here are examples of both:

                PH-1	Be-10	quartz	123453	3717	07KNSTD;
                PH-1	Al-26	quartz	712408	31238	KNSTD;
                

B.1.1 Sample name. This must exactly match the sample name of the corresponding sample as entered in the sample data line, or else, obviously, the code will not be able to match samples to nuclide concentration measurements.

B.1.2. and B.1.3. Nuclide-mineral pair. For Be-10 in quartz, enter "Be-10" and "quartz"; for Al-26 in quartz, this is "Al-26" and "quartz."

B.1.4. and B.1.5. Nuclide concentration and uncertainty. Atoms per gram. Standard or scientific notation.

B.1.6. Name of Be-10 or Al-26 standardization. This is familiar from version 2. Acceptable values for this parameter are given below.

B.1.7. Then the line ends with a semicolon.


Standardization:

See the official documentation.

Calculators home
Exposure ages
Erosion rates

This page provides a list of standardizations to which Be-10 and Al-26 measurements submitted to the online exposure age and erosion rate calculators can be referenced. A 'standardization' means the combination of a particular isotope ratio standard material and an assumed isotope ratio for that material. Each standardization described below is associated with a name that appears in the first column. This name must be entered exactly in the appropriate place in the online calculator input forms. Note that two pieces of information are critical here: the identity of the actual standard material, and the nominal isotope ratio that the material is assumed to have. In a number of cases, the same physical standard material has been used with different assumed isotope ratios.

The goal here is that users not have to renormalize their AMS measurements before submitting them to the online calculators. For all the standardizations listed below, there is a known conversion factor that can be used to convert measurements made using that standardization to be compatible with the online calculators. The calculators will do that conversion internally before calculating exposure ages or erosion rates. Thus, users are NOT required to renormalize data to a particular standardization before submitting it to the calculator — the idea is that the user lists the standardization that was used for a particular sample, and the correction is done internally in a consistent fashion.

The calculator multiplies a nuclide concentration measured using a given standardization with a conversion factor to make it consistent with the 07KNSTD (for Be-10) or KNSTD (for Al-26) standardization.
A table with conversion factors last updated February 25, 2016 is available in this PDF.

If you are using a standard material/isotope ratio pair that is not listed here, and you would like it to be added, please contact Greg Balco.


Be standardizations

Name Description
07KNSTD
Any of a dilution series derived from the so-called "ICN solution" by K. Nishiizumi and described in Nishiizumi et al, 2007 (NIM-B v. 258, p. 403), with the revised nominal isotope ratios listed in that publication and in the printed description of the standards available from K. Nishiizumi and dated May 29, 2007. Measurements made at LLNL-CAMS with these standards and nominal isotope ratios will list '07KNSTDX,' where 'X' is a number related to the isotope ratio of the standard. Be-10 measurements made at PRIME Lab after November 14, 2007 were referenced to this set of standards and nominal isotope ratios. This is the standard on which the internal constants and production rates in the online exposure age calculator are based, so measurements made against any other standard will be internally converted to be consistent with this standard.
KNSTD
This refers to the same standard material as above -- the dilution series derived from the ICN solution by K. Nishiizumi -- but with a different nominal isotope ratio that was assumed for these standards before the 2007 revision. Measurements made at LLNL-CAMS with this set of standards and isotope ratios will list 'KNSTDX' as the name of the standard. Measurements made at PRIME Lab between January 12, 2005 and November 14, 2007 were referenced to this set of standards and isotope ratios.
NIST_Certified
This refers to a standard produced by the National Institute of Standards and Technology (NIST), referred to as SRM4325, with the nominal isotope ratio stated on the certificate for this material (2.68 x 10^-11 for the solution as supplied by NIST). Measurements made at PRIME Lab prior to January 12, 2005 were referenced to this standard and isotope ratio.
NIST_30000
NIST SRM4325 standard material with an assumed isotope ratio of 3.0 x 10^-11 rather than the NIST certified value.
NIST_30200
NIST SRM4325 standard material with an assumed isotope ratio of 3.02 x 10^11 rather than the NIST certified value.
NIST_30300
NIST SRM4325 standard material with an assumed isotope ratio of 3.03 x 10^-11 rather than the NIST certified value.
NIST_30600
NIST SRM4325 standard material with an assumed isotope ratio of 3.06 x 10^-11 rather than the NIST certified value.
NIST_27900
NIST SRM4325 standard material with an assumed isotope ratio of 2.79 x 10^-11 rather than the NIST certified value. This standardization is equivalent to 07KNSTD within rounding error, so users can enter either one.
BEST433
ETH-Zurich standard material "BEST433" with an assumed isotope ratio of 93.1 x 10^-12. This standardization was in use at ETH before April 1, 2010. See Kubik and Christl (2010).
S555
ETH-Zurich standard material "S555" with an assumed isotope ratio of 95.5 x 10^-12. This standardization was in use at ETH before April 1, 2010. See Kubik and Christl (2010).
S2007
ETH-Zurich standard material "S2007" with an assumed isotope ratio of 30.8 x 10^-12. This standardization was in use at ETH prior to April 1, 2010. See Kubik and Christl (2010).
BEST433N
ETH-Zurich standard material originally called "BEST433" with a revised isotope ratio of 83.3 x 10^-12. This standardization is equivalent to 07KNSTD, so users can enter either one. This standardization was adopted at ETH on April 1, 2010. See Kubik and Christl (2010).
S555N
ETH-Zurich standard material originally called "S555" with a revised isotope ratio of 87.1 x 10^-12. This standardization is equivalent to 07KNSTD, so users can enter either one. This standardization was adopted at ETH on April 1, 2010. See Kubik and Christl (2010).
S2007N
ETH-Zurich standard material originally called "S2007" with a revised isotope ratio of 28.1 x 10^-12. This standardization is equivalent to 07KNSTD, so users can enter either one. This standardization was adopted at ETH on April 1, 2010. See Kubik and Christl (2010).
STD11
ASTER standard called "STD-11" with an assumed isotope ratio of 1.191 x 10^11. This standard was calibrated by reference to NIST_27900, which is equivalent to 07KNSTD at rounding error, so all three of these should be equivalent.
LLNL31000
An internal LLNL-CAMS standard. Identified as such on CAMS data reports.
LLNL10000
An internal LLNL-CAMS standard. Identified as such on CAMS data reports.
LLNL3000
An internal LLNL-CAMS standard. Identified as such on CAMS data reports.
LLNL1000
An internal LLNL-CAMS standard. Identified as such on CAMS data reports.
LLNL300
An internal LLNL-CAMS standard. Identified as such on CAMS data reports.

Al standardizations

Name Description
KNSTD
Any of a dilution series described in Nishiizumi, 2004 (NIM-B, v. 223-224, p. 388), with the assumed isotope ratios described in this publication. Measurements made at LLNL-CAMS with this standard will list 'KNSTDX,' where 'X' is a number related to the isotope ratio of the standard. This is the standard on which the internal constants and production rates in the online exposure age calculator are based.
ZAL94
ETH-Zurich standard material "ZAL94" with an assumed isotope ratio of 526 x 10^-12. This standardization was in use at ETH prior to April 1, 2010. The University of Cologne "AL09" standard material with an assumed isotope ratio of 1190 x 10^-12 is also consistent with this standardization. See Kubik and Christl (2010).
ZAL94N
ETH-Zurich standard material originally called "ZAL94" with a revised isotope ratio of 490 x 10^-12. This standardization is equivalent to KNSTD, so users can enter either one. This standardization was adopted at ETH on April 1, 2010. See Kubik and Christl (2010).
SMAL11
Internal standard used at ASTER, named "SM-Al-11" and with defined Al-26/Al-27 = 7.401 x 10^-12. This is part of a dilution series that also includes other standards called "SM-Al-10" (9.352 x 10^-11), "SM-Al-12" (7.21 x 10^-13), and "SM-Al-13" (7.30 x 10^-14).
Z92-0222
Al standard originally prepared at Purdue and used at several other labs with a defined isotope ratio of 4.11 x 10^-11. This standardization is equivalent to KNSTD, so users can enter either one.