RecordCompression: Calculate compression of firn-core record

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

Calculate compression of firn-core record

Description

Calculate the compression due to densification with a constant rate during downward-advection of a firn-core record.

Usage

RecordCompression(z.start = 0, adv, length.in, rate, res = 1)

Arguments

z.start	the depth of the upper end of the firn-core record in [cm].
adv	the amount of downward advection in [cm].
length.in	the initial length of the firn-core record in [cm].
rate	constant rate of firn densification relative to the surface density; in [% per metre].
res	resolution (in [cm]) of the compression calculation, i.e. the size of the individual depth bins of the record ("snow parcels") that are advected downwards.

Value

the length change of the firn-core record in [cm].

Author(s)

Thomas Münch

References

Münch, T., et al., Constraints on post-depositional isotope modifications in East Antarctic firn from analysing temporal changes of isotope profiles, The Cryosphere, 11(5), 2175-2188, doi:10.5194/tc-11-2175-2017, 2017.

See Also

Compression

Examples

# Length change of trench record for mean estimated densification rate
# (see Münch et al., 2017):
RecordCompression(adv = 50, length.in = 100, rate = 4.5)

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