Convert data from the depth to the time domain

Converts a data set from the depth to the time domain using a tracked curve/cycle to depth domain an assigning a duration (in kyr) set tracked curve/cycle.

Usage

curve2tune(
  data = NULL,
  tracked_cycle_curve = NULL,
  tracked_cycle_period = NULL,
  genplot = FALSE,
  keep_editable = FALSE
)

Arguments

data

Record to convert to the time domain (matrix or data frame, two columns): depth in metres in column 1 and the proxy in column 2, normally the data set that was analysed with analyse_wavelet to track the cycle. Depths outside the tracked curve get the age of its nearest end. Required (Default=NULL gives an error).

tracked_cycle_curve

Tracked period curve of one astronomical cycle (matrix or data frame, two columns): depth in metres in column 1 and the period of the tracked cycle in metres in column 2, e.g. the output of track_period_wavelet (preferably completed with completed_series). Rows with NA are dropped. Required (Default=NULL gives an error).

tracked_cycle_period

Duration of the tracked cycle in kyr (single numeric value, > 0, e.g. 405 for the long eccentricity cycle), used to convert the tracked period to a sedimentation rate that is integrated to time. Required (Default=NULL gives an error).

genplot

Plot three panels stacked on top of each other (logical, TRUE or FALSE): the data in the depth domain, the depth-time curve and the data in the time domain. Default=FALSE.

keep_editable

Keep the graphical parameters of the plot so that extra elements can be added afterwards (logical, TRUE or FALSE); FALSE restores the previous par settings. Default=FALSE.

Value

The output is a matrix with 2 columns. The first column is time. The second column sedimentation proxy value.

If genplot=TRUE then 3 plots stacked on top of each other will be plotted. Plot 1: the original data set. Plot 2: the depth time plot. Plot 3: the data set in the time domain.

References

Routines for astrochronologic testing, astronomical time scale construction, and time series analysis <doi:10.1016/j.earscirev.2018.11.015>

Author

Part of the code is based on the sedrate2time function of the 'astrochron' R package

Examples

# \donttest{
#The example uses the magnetic susceptibility data set of Pas et al., (2018).
# perform the CWT
mag_wt <- analyse_wavelet(data = mag,
dj = 1/100,
lowerPeriod = 0.1,
upperPeriod = 254,
verbose = FALSE,
omega_nr = 10)

#Track the 405 kyr eccentricity cycle in a wavelet spectra

#mag_track <- track_period_wavelet(astro_cycle = 405,
#                                   wavelet=mag_wt,
#                                   n.levels = 100,
#                                   periodlab = "Period (meters)",
#                                   x_lab = "depth (meters)")

#Instead of tracking, the tracked solution data set mag_track_solution is used
mag_track <- mag_track_solution

mag_track_complete <- completed_series(
  wavelet = mag_wt,
  tracked_curve = mag_track,
  period_up = 1.2,
  period_down = 0.8,
  extrapolate = TRUE,
  genplot = FALSE
)

# smooth the tracking of the 405 kyr eccentricity cycle
mag_track_complete <- loess_auto(time_series = mag_track_complete,
genplot = FALSE, print_span = FALSE)

mag_track_time<- curve2tune(data=mag,
                           tracked_cycle_curve=mag_track_complete,
                           tracked_cycle_period=405,
                           genplot = FALSE,
                           keep_editable=FALSE)
# }