Compute constant depth and denudation rate curves array for two-nuclides plots

tnp_soil_mixing(h, E, rhos, rhob, p1, p2, L, S, n = 100)

Arguments

h

soil thickness vector (cm)

E

erosion rate vector (m/Ma)

rhos

soil density (g/cm3)

rhob

bedrock density (g/cm3)

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

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 soil depth 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
rhob = 2.65 # bedrock density (g/cm3)
rhos = rhob/2 # soil density  (g/cm3)
N1 = "Be10" # longer half-life
N2 = "Al26" # shorter half-life
res = tnp_curves(prm[,N1],prm[,N2],Lambda,S,rhob)
plot(NA,xlim=c(0.75,5),ylim=c(3,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
h = seq(100,1000,by = 100) # increments in soil depth (cm)
E = c(0.1,0.2,0.5,1,2,5,10) # increments in denudation rate (m/Ma)
res = tnp_soil_mixing(h,E,rhos,rhob,prm[,N1],prm[,N2],L,S,n=100) # compute array
lines(res[[1]]$C1/1e6,res[[1]]$C2/res[[1]]$C1,col="grey") # constant depth
lines(res[[2]]$C1/1e6,res[[2]]$C2/res[[2]]$C1,col="black") # constant denudation