DensificationRate: Firn densification rate

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

Firn densification rate

Description

Function to return the relative densification rate in % per metre over a specified depth, given a record of firn densities. The densification rate is found from linear regression of the firn density profile against depth.

Usage

DensificationRate(dat = NULL, lim)

Arguments

dat	a list providing a firn density profile with two elements: depth, a vector with the depth scale of the density profile in metre; and density, a vector of firn densities. If NULL (the default), the stacked B41/B42 density profile, provided with the package, is used for the calculation (see Münch et al., 2017).
lim	the depth (in m) over which the densification rate shall be calculated.

Value

the relative linear densification rate (in % per metre) from zero depth to the depth lim.

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.

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