timetostratratebin {DAIME}R Documentation

Transform Rates from Time into Stratigraphic Height (Binned)

Description

This function transforms a given binned temporal rate/signal into the corresponding stratigraphic rate/signal. The transformation is based on the age model constructed from the deposition rate.

This can for example be used to predict the patterns of (first/last) fossil occurrences in the outcrop under a given depositional environment.

Usage

timetostratratebin(binborder, depoval, signalval, pos = NULL)

Arguments

binborder

Vector of strictly increasing numerical values. Defines the borders of the bins in time

depoval

Strictly positive vector of length length(binborder)-1. Defines the deposition rate in the bins, i.e. depoval[i] is the deposition rate in the time interval between binborder[i] and binborder[i+1]

signalval

Positive vector of length length(binborder)-1. Defines the temporal rate/signal in the bins, i.e. signalval[i] is the temporal rate between the times binborder[i] and binborder[i+1]

pos

OPTIONAL. Vector of points in time that will be transformed into stratigraphic rate, and where the stratigraphic rate will evaluated. Default setting is a vector spanning the interval of interest with evenly spaced points (no. is chosen by the code)

Value

A list containing the following entries

x

x values, corresponding to the locations of pos, transformed into stratigraphic height

y

function values of the transformed rate at x

The results return NA if the values coincide with a hiatus or both deposition rate and temporal rate are undefined

Author(s)

Niklas Hohmann

References

Hohmann, Niklas. 2018. Quantifying the Effects of Changing Deposition Rates and Hiatii on the Stratigraphic Distribution of Fossils. <doi:10.13140/RG.2.2.23372.51841>

See Also

For an overview of the functions in the DAIME package, see its vignette (available via vignette("DAIME") )

Examples

## Define deposition rate
binborder=1:6 #temporal bins for the deposition rate
depoval=c(5,4,3,1,2) #deposition rate in the bins
#plot deposition rate
depositionrate=approxfun(binborder,c(depoval,tail(depoval,1)),method="constant",yleft=NA,yright=NA)
plot(depositionrate(seq(from=min(binborder),to=max(binborder),length.out=100)),xlab='time',
  ylab='deposition rate', main='Deposition rate')
## Define temporal rate
signalval=c(1,0,5,2,1)
#plot temporal rate
temporalrate=approxfun(binborder,c(signalval,tail(signalval,1)),method="constant",
  yleft=NA,yright=NA)
plot(temporalrate(seq(from=min(binborder),to=max(binborder),length.out=100)),xlab='time',
  ylab='temporal rate', main='Temporal rate')

## Transform temporal rate into stratigraphic rate
reslist=timetostratratebin(binborder,depoval,signalval)
#plot resulting stratigraphic rate
plot(reslist$height,reslist$val,xlab='Stratigraphic Height',ylab='Intensity',
  main='Stratigraphic Rate')
## With removal of sediment (hiatus)
depoval=c(5,4,-3,1,2) #erosion rate is 3 in the middle time bin
reslist=timetostratratebin(binborder,depoval,signalval)
#plot resulting stratigraphic rate
plot(reslist$height,reslist$val,xlab='Stratigraphic Height',ylab='Intensity',
  main='Stratigraphic Rate')

[Package DAIME version 1.1 Index]