R/calc_CobbleDoseRate.R
calc_CobbleDoseRate.RdCalculates the dose rate profile through the cobble based on Riedesel and Autzen (2020).
Corrects the beta dose rate in the cobble for the grain size following results of Guérin et al. (2012). Sediment beta and gamma dose rates are corrected for the water content of the sediment using the correction factors of Aitken (1985). Water content in the cobble is assumed to be 0.
calc_CobbleDoseRate(input, conversion = "Guerinetal2011")data.frame (required): A table containing all relevant information for each individual layer. For the table layout see details.
Which dose rate conversion factors to use. For accepted values see BaseDataSet.ConversionFactors
The function returns an RLum.Results object for which the first element
is a matrix (DataIndividual) that gives the dose rate results for each slice
for each decay chain individually, for both, the cobble dose rate and the sediment
dose rate. The second element is also a matrix (DataComponent) that gives
the total beta and gamma-dose rates for the cobble and the adjacent sediment
for each slice of the cobble.
The input table layout
| COLUMN | DATA TYPE | DESCRIPTION |
Distance | numeric | distance from the surface of the cobble to the top of each rock slice in mm. The distance for each slice will be listed in this column |
DistanceError | numeric | Error on the distance in mm |
Thickness | numeric | Thickness of each slice in mm |
TicknessError | numeric | uncertainty of the thickness in mm. |
Mineral | character | 'FS' for feldspar, 'Q' for quartz, depending which mineral in the cobble is used for dating |
Cobble_K | numeric | K nuclide content in % of the bulk cobble |
Cobble_K_SE | numeric | error on K nuclide content in % of the bulk cobble |
Cobble_Th | numeric | Th nuclide content in ppm of the bulk cobble |
Cobble_Th_SE | numeric | error on Th nuclide content in ppm of the bulk cobble |
Cobble_U | numeric | U nuclide content in ppm of the bulk cobble |
CobbleU_SE | numeric | error on U nuclide content in ppm of the bulk cobble |
GrainSize | numeric | average grain size in µm of the grains used for dating |
Density | numeric | Density of the cobble. Default is 2.7 g cm^-3 |
CobbleDiameter | numeric | Diameter of the cobble in cm. |
Sed_K | numeric | K nuclide content in % of the sediment matrix |
Sed_K_SE | numeric | error on K nuclide content in % of the sediment matrix |
Sed_Th | numeric | Th nuclide content in ppm of the sediment matrix |
Sed_Th_SE | numeric | error on Th nuclide content in ppm of the sediment matrix |
Sed_U | numeric | U nuclide content in ppm of the sediment matrix |
Sed_U_SE | numeric | error on U nuclide content in ppm of the sediment matrix |
GrainSize | numeric | average grain size of the sediment matrix |
WaterContent | numeric | mean water content of the sediment matrix in % |
WaterContent_SE | numeric | relative error on water content |
Water content The water content provided by the user should be calculated according to:
$$(Wet_weight - Dry_weight) / Dry_weight * 100$$
0.1.0
Riedesel, S., Autzen, M., 2023. calc_CobbleDoseRate(): Calculate dose rate of slices in a spherical cobble. Function version 0.1.0. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J., Mercier, N., Philippe, A., Riedesel, S., Autzen, M., Mittelstrass, D., Gray, H.J., Galharret, J., 2023. Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.9.23. https://CRAN.R-project.org/package=Luminescence
Riedesel, S., Autzen, M., 2020. Beta and gamma dose rate attenuation in rocks and sediment. Radiation Measurements 133, 106295.
## load example data
data("ExampleData.CobbleData", envir = environment())
## run function
calc_CobbleDoseRate(ExampleData.CobbleData)
#>
#> [RLum.Results-class]
#> originator: calc_CobbleDoseRate()
#> data: 3
#> .. $DataIndividual : matrix
#> .. $DataComponent : matrix
#> .. $input : data.frame
#> additional info elements: 1