Compute burial age and denudation rate curves array for two-nuclides plots

tnp_burial(A, E, p1, p2, L, S, rho, n = 100)

Arguments

A

burial age vector (a)

E

erosion rate vector (m/Ma)

p1

production and decay parameters for the first nuclide (4 elements vector)

p1[1] unscaled spallation production rate (at/g/a)
p1[2] unscaled stopped muons production rate (at/g/a)
p1[3] unscaled fast muons production rate (at/g/a)
p1[4] decay constant (1/a)

p2

production and decay parameters for the second nuclide (4 elements vector, same as p1)

L

Attenuation length (3 elements vector in g/cm2)

L[1] neutrons
L[2] stopped muons
L[3] fast muons

S

scaling factors (2 elements vector)

S[1] scaling factor for spallation
S[2] scaling factor for muons

rho

material density (g/cm3)

n

number of along-curve evaluation points (optional default 100)

Value

a list with two dataframes containing the concentrations for the two nuclides as a function of burial age and denudation

Examples

data("prm") # production and decay data
p = prm
data("Lambda") # attenuation length data
L = Lambda
altitude = 1000 # site elevation in m
latitude = 20 # site latitude in degrees
P = atm_pressure(alt=altitude,model="stone2000") # atmospheric pressure at site
S = scaling_st(P,latitude) # Stone 2000 scaling parameters
rho = 2.65 # bedrock density (g/cm3)
N1 = "Be10" # longer half-life
N2 = "Al26" # shorter half-life
res = tnp_curves(prm[,N1],prm[,N2],Lambda,S,rho)
plot(NA,xlim=c(0.25,10),ylim=c(0.5,7),log="x",
    xlab=paste(N1,"(x10^6 at/g)"),ylab=paste(N2,"/",N1))
lines(res[[1]]$C1/1e6,res[[1]]$C2/res[[1]]$C1,lty=2,lwd=2,col="khaki4") # constant exposure
lines(res[[2]]$C1/1e6,res[[2]]$C2/res[[2]]$C1,lty=1,lwd=2,col="khaki4") # steady-state erosion
A = c(1,2,3,4)*1e6 # increments in age (a)
E = c(0.1,0.2,0.5,1,2,5,10) # increments in denudation rate (m/Ma)
res = tnp_burial(A,E,prm[,N1],prm[,N2],L,S,rho,n=100) # compute array
lines(res[[1]]$C1/1e6,res[[1]]$C2/res[[1]]$C1,col="grey") # constant age
lines(res[[2]]$C1/1e6,res[[2]]$C2/res[[2]]$C1,col="black") # constant denudation