ModifyRecord: Model temporal change of firn profile

In EarthSystemDiagnostics/TrenchR: Analyse and Plot Trench-like Proxy Records

Model temporal change of firn profile

Description

Modify a proxy firn profile (e.g., an isotope record) to simulate changes over time given values of downward advection, compression by densification, and diffusional smoothing.

Usage

ModifyRecord(
  rec.in,
  res,
  depth.hires,
  depth.lores = NULL,
  SIGMA = NULL,
  STRETCH = NULL,
  ADV = NULL
)

Arguments

rec.in	numeric vector with the firn proxy (isotope) record to be modified.
res	the (high) depth resolution of the firn record on which the the modification processes are calculated.
depth.hires	numeric vector with the (high-resolution) depth scale on which rec.in is recorded; must be of the same length as rec.in and in units of res.
depth.lores	numeric vector with a lower resolution depth scale of rec.in; in the units of res. This is optional: A sufficiently high depth resolution is needed for proper simulation of the modification processes; if rec.in and depth.hires provide an interpolated version of the original record to achieve such a higher resolution, depth.lores can be used to provide the original lower resolution depth scale, and the modified version of the record is then returned on both the interpolated (high) and the original (low) depth resolution.
SIGMA	diffusion length value to smooth the record; must be in the same units as res. If NULL, no smoothing by diffusion is modelled.
STRETCH	value of the amount of compression due to densification of the original depth scale of record over time; must be in the same units as res. If NULL, no compression by densification is modelled.
ADV	value of downward advection of the record over time; must be in the same units as res. If NULL, no downward advection is modelled.

Value

A list with one or two elements:

HiRes:

the modified record on the high depth resolution res;

LoRes:

if depth.lores is not NULL, the modified record on this lower depth resolution.

Author(s)

Thomas Münch

See Also

DiffuseRecord; CompressRecord

EarthSystemDiagnostics/TrenchR documentation built on Sept. 2, 2023, 3:44 a.m.