impulseResponse {astrochron} | R Documentation |
Calculate the analytical response function from an impulse forcing using the 1-D advection diffusion model proposed in Schink and Guinasso (1975) to model the bioturbation impact on climate time series.
impulseResponse(G, ML = NULL, v = NULL, nt = 500, genplot = FALSE,
verbose = FALSE)
G |
Bioturbation parameter. G = D/ML/v |
ML |
Mix layer depth (cm) |
v |
Sedimentation rate (cm/kyr) |
nt |
Number of steps after the signal is deposited. |
genplot |
Generate summary plots? (T or F) |
verbose |
Verbose output? (T or F) |
fc Impulse response function
Guinasso, N.L. and Schinck, D.R., 1975, Quantitative estimates of biological mixing rates in abyssal sediments, J. Geophys. Res., 80, 3032-3043.
Goreau, T.J., 1977, Quantitative effects of sediment mixing on stratigraphy and biogeochemistry: a signal theory approach, Nature, 256, 730-732.
Goreau, T.J., 1980, Frequency sensitivity of the deep-sea climatic record, Nature, 287, 620-622.
Liu, H., Meyers, S.R., and Marcott, S.A., 2021, Unmixing dee-sea paleoclimate records: A study on bioturbation effects through convolution and deconvolution, Earth and Planetary Science Letters, 564, 116883.
G <- 4
ML <- 10
v <- 1
# take a look at the IRF
impulseResponse(G=4, ML = 10, v = 1, genplot = TRUE)