S_E: Get error components for stochastic climate signal

In EarthSystemDiagnostics/psem: Proxy Spectral Error Model

Get error components for stochastic climate signal

Description

Get error components for stochastic climate signal

Usage

S_E(
  nu,
  tau_s,
  tau_b,
  tau_r,
  T,
  delta_t,
  N,
  n.k,
  clim.spec.fun,
  clim.spec.fun.args,
  n.nu.prime = 1000,
  plot.integral = FALSE
)

Arguments

nu	frequency
tau_s	sediment slice thickness in years (layer.width / sedimentation rate)
tau_b	timescale of bioturbation (bioturbation depth / sedimentation rate) (L/sr)
tau_r	width of moving average filter that represents the "interpreted" resolution of the timeseries
T	length of a finite time series, e.g. 1e04
delta_t	sampling frequency of the sediment core / climate timeseries
N	number of signal carriers (e.g. individual foraminifera)
n.k	the number of aliasers used when estimating the error spectrum of the stochastic climate signal. Defaults to 15, does not normally need to be adjusted.
clim.spec.fun	a function to return climate power spectral density as a function of frequency, nu
clim.spec.fun.args	arguments to the named climate power spectrum function
n.nu.prime	number of discrete frequencies at which to evaluate the PSD of the error
plot.integral	Plot the PSD of the error of the bioturbated climate timeseries. This PSD is numerically integrated to give the variance for the aliasing of climate variation.

Value

list of frequencies and error spectrum components

See Also

ProxyErrorSpectrum

Examples

spec.pars <- GetSpecPars("Mg_Ca", phi_c = pi, delta_phi_c = pi,
 sig.sq_a = 0.1, tau_p = 1/12)
spec.pars$nu <- GetNu(spec.pars$T, spec.pars$delta_t)
do.call(S_E, spec.pars[names(spec.pars) %in% names(formals(S_E))])

EarthSystemDiagnostics/psem documentation built on March 3, 2024, 10:24 p.m.