Have a personal or library account? Click to login
Analyzing Statistical Age Models to Determine the Equivalent Dose and Burial Age Using a Markov Chain Monte Carlo Method Cover

Analyzing Statistical Age Models to Determine the Equivalent Dose and Burial Age Using a Markov Chain Monte Carlo Method

Geochronometria
Volume 48 (2021): Issue 1 (January 2021)
By: Jun Peng  
Open Access
|Dec 2021

Figures & Tables

Fig. 1

Statements used to conduct MCMC sampling according to the CAM and the MAM3 (log-scale) using De sets from individual samples. The MXAM3 can be applied in a similar manner to the MAM3 by changing the likelihood function.

Fig. 2

Statements used for conducting MCMC simulation simultaneously according to different age models (log-scale) using De sets from three samples whose ages were constrained such that a1<a2<a3. The three samples were analyzed using the MAM3, CAM, and MXAM3, respectively. The protocol can be easily adapted for MCMC sampling without order constrains for burial ages.

Fig. 3

Measured multi-grain De distributions for four aeolian samples from Gulang County at the southern margin of the Tengger Desert, China (Peng et al., 2016b). N denotes the number of measured aliquots. Note that in Peng et al. (2016b) an additional uncertainty of σb=5% was added in the quadrature to the RSEs of the measured De values to account for sources of errors that were not considered in the measurements. However, no additional uncertainty was added to the De distributions shown here.

Fig. 4

Posterior distributions of burial ages obtained from the CAM for the four aeolian samples taken from the same sedimentary section obtained using the MCMC sampling protocol, as shown in Fig. 2. The blue-coloured density plots denote results obtained by constraining the burial ages in ascending order. The grey-coloured density plots denote results obtained by imposing no constraints on the order of burial ages. The value inside parentheses denotes the calculated RSD of burial age.

Fig. 5

Distribution of OSL sensitivity for the 127 grains of a sample taken from Lake Mungo, Australia. The red line denotes the probability density curve obtained by fitting the OSL sensitivity data using a Gamma distribution. The fitted Gamma distribution has a shape parameter of α=0.565 and a rate parameter of β=0.191. Q[x%] denotes the x% sample quantile of OSL sensitivity.

Fig. 6

Simulated De distributions of nine samples from a sedimentary sequence. Each De set contains 100 De values (and associated standard errors). The shaded bars in each plot indicate the burial dose within 2σ error. The first three samples (#1–#3) are partially-bleached, the next three samples (#4–#6) are fully-bleached, and the last three samples (#7–#9) contain fully-bleached grains that have been subsequently mixed with younger, intrusive grains. The nine samples are assumed to be collected from the same section but at different depths, and their estimated burial dose and age are summarized in Fig. 7 and Table 3.

Fig. 7

Posterior distributions of burial ages of the nine simulated samples from a sedimentary section obtained using the MCMC sampling protocol, as shown in Fig. 2. Samples #1–#3, #4–#6, and #7–#9 were analyzed using the MAM3, CAM, and MXAM3, respectively. The blue-coloured density plots are results obtained by constraining burial ages in ascending order. The grey-coloured density plots are results obtained by imposing no constraints on the order of burial ages. The first and second values inside parentheses denote the calculated RSD and RE of the burial age, respectively.

Fig. 8

Convergence diagnostics for posterior samples of the CAM burial age from sample GL2-1. The plots were automatically generated using function mcmcSAM(). (A) denotes the trace plot; (B) denotes the density plot; (C) denotes the autocorrelations plot, and (D) denotes the Gelman-Rubin convergence diagnostic plot generated using three parallel Markov chains.

A summary of burial dose and age estimated from various statistical age models for nine simulated De sets taken from the same sedimentary sequence_ Results obtained using MLE and results obtained using the MCMC sampling protocol as shown in Fig_ 2 without and with order constraints are presented_ The systematic error on the dose rate measurement shared by all simulated samples was assumed to be σdc=0_

SampleModelAge (ka)Burial dose (Gy)MLEMCMC (without constraints)MCMC (with constraints)
μ (Gy)a (ka)μ (Gy)a (ka)μ (Gy)a (ka)
#1MAM371818.35 ± 1.177.95 ± 0.7517.85 ± 1.177.82 ± 0.7616.58 ± 1.026.80 ± 0.44
#27.52121.16 ± 0.636.95 ± 0.5320.96 ± 0.796.96 ± 0.5721.17 ± 0.657.23 ± 0.40
#382424.47 ± 1.228.41 ± 0.7223.64 ± 1.578.22 ± 0.8123.25 ± 1.437.84 ± 0.42
#4CAM8.52726.91 ± 0.338.24 ± 0.5926.92 ± 0.348.32 ± 0.6126.93 ± 0.338.31 ± 0.39
#593030.40 ± 0.388.50 ± 0.6030.41 ± 0.398.59 ± 0.6330.44 ± 0.388.83 ± 0.42
#69.53332.79 ± 0.4610.59 ± 0.7632.81 ± 0.4710.71 ± 0.7932.66 ± 0.459.62 ± 0.42
#7MXAM3103635.77 ± 1.619.13 ± 0.7636.43 ± 1.869.39 ± 0.8437.18 ± 1.819.95 ± 0.45
#810.53939.26 ± 1.379.32 ± 0.7339.70 ± 1.699.53 ± 0.7940.62 ± 1.8710.39 ± 0.54
#9114243.65 ± 2.5510.63 ± 0.9744.61 ± 2.5210.96 ± 1.0045.43 ± 2.5311.48 ± 0.85

A summary of burial dose and age estimated from the CAM for De sets from four measured aeolian samples using the MCMC sampling protocol shown in Fig_ 2 without and with order constraints_ The systematic error for the dose rate measurement shared by all samples was set to σdc=0_1_

SampleDose rate (Gy/ka)MCMC (without constraints)MCMC (with constraints)
μ (Gy)a (ka)μ (Gy)a (ka)
GL2-13.26 ± 0.2534.55 ± 0.9110.75 ± 0.9634.23 ± 0.869.72 ± 0.52
GL2-22.84 ± 0.2132.71 ± 1.0011.68 ± 1.0632.19 ± 0.9210.28 ± 0.49
GL2-33.23 ± 0.2330.58 ± 1.049.59 ± 0.8531.03 ± 1.0010.59 ± 0.52
GL2-43.40 ± 0.2432.24 ± 1.299.60 ± 0.8533.13 ± 1.2611.13 ± 0.67

Comparisons of burial dose and age estimated from MLE and MCMC for the various statistical age models using De sets from four measured aeolian samples_ Quantities estimated using the MCMC sampling protocol shown in Fig_ 1 are marked in bold_ Quantities estimated using the MLE are inside parentheses_

SampleDose rate (Gy/ka)CAMMAM3MXAM3
μ (Gy)a (ka)μ (Gy)a (ka)μ (Gy)a (ka)
GL2-13.26 ± 0.2534.55 ± 0.91 (34.51 ± 0.84)10.73 ± 0.89 (10.58 ± 0.85)33.96 ± 1.17 (32.79 ± 3.64)10.55 ± 0.89 (10.05 ± 1.36)35.36 ± 1.41 (35.98 ± 3.51)10.96 ± 0.97 (11.04 ± 1.37)
GL2-22.84 ± 0.2132.70 ± 0.98 (32.65 ± 0.91)11.64 ± 0.95 (11.50 ± 0.91)30.40 ± 1.41 (30.33 ± 1.67)10.83 ± 0.98 (10.68 ± 0.98)35.53 ± 2.12 (37.12 ± 2.31)12.66 ± 1.26 (13.7 ± 1.26)
GL2-33.23 ± 0.2330.57 ± 1.05 (30.55 ± 0.98)9.56 ± 0.77 (9.46 ± 0.74)27.92 ± 1.55 (26.03 ± 3.56)8.74 ± 0.80 (8.06 ± 1.24)35.19 ± 2.03 (36.61 ± 2.41)11.00 ± 1.01 (11.33 ± 1.10)
GL2-43.40 ± 0.2432.23 ± 1.27 (32.17 ± 1.18)9.58 ± 0.79 (9.46 ± 0.75)30.09 ± 1.55 (30.75 ± 1.04)8.94 ± 0.80 (9.04 ± 0.71)37.77 ± 2.61 (39.16 ± 2.46)11.22 ± 1.12 (11.52 ± 1.09)

A summary of Gelman-Rubin convergence diagnostics for measured and simulated samples obtained using the MCMC sampling protocol, as shown in Fig_ 2 with order constraints_ n_eff is the effective sample size, while Rhat (i_e_, the shrink factor) is a statistic measure of the ratio of the average variance of samples within each chain to the variance of the pooled samples across chains_ If all chains are at equilibrium, the Rhat will be 1_ If these chains have not converged to a common distribution, the Rhat statistic will be greater than 1_

SampleμaSampleμa
n_effRhatn_effRhatn_effRhatn_effRhat
GL2-1160000.99985541.741.0002#113865.371.00019987.931.0000
GL2-2160000.99998220.181.0003#2160001.000111473.551.0000
GL2-3160000.99977595.601.0002#39405.951.000412466.841.0003
GL2-4160000.99988864.011.0003#4160001.000313161.070.9999
#5160000.999912830.440.9999
#6160000.999914048.980.9998
#7160000.999913073.690.9998
#87780.821.000312665.561.0001
#9160000.9999160000.9999
DOI: https://doi.org/10.1515/geochr-2015-0114 | Journal eISSN: 1897-1695
Journal RSS Feed
Language: English
Page range: 147 - 160
Submitted on: Jan 7, 2019
|
Accepted on: Jul 12, 2019
|
Published on: Dec 31, 2021
Published by: Sciendo
In partnership with: Paradigm Publishing Services
Publication frequency: 1 issue per year
Keywords:
OSL dating,
statistical age models,
equivalent dose,
MCMC sampling
Related subjects:
Geosciences,
Geosciences, other

© 2021 Jun Peng, published by Sciendo
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Previous articleVolume 48 (2021): Issue 1 (January 2021)Next article