Skip to contents

Estimate the average organic carbon sequestration rate

Source: R/estimate_seq_rate.R
estimate_seq_rate.Rd

estimate the average organic carbon sequestration rate to the soil in a indicated time frame (by default last 100 years) from the organic carbon concentration and ages obtained from a age-depth or age-accumulated mass model

Usage

estimate_seq_rate(
  df = NULL,
  timeframe = 100,
  core = "core",
  mind = "mind_corrected",
  maxd = "maxd_corrected",
  dbd = "dbd",
  oc = "eoc",
  age = "age"
)

Arguments

df

A data.frame with, at least, columns: core, mind (minimum depth of the sample), maxd (maximum depth of the sample), dbd (dry bulk density), oc (organic carbon %), age (age of the sample obtained from a age-depth or age-accumulated mass model)

timeframe

standardization time frame, by default 100 years

core

Character Name of the column reporting core ID.

mind

Character Name of the column reporting the minimum depth of each sample.

maxd

Character Name of the column reporting the maximum depth of each sample.

dbd

Character Name of the column reporting dry bulk density.

oc

Character Name of the column reporting organic carbon concentrations.

age

Character Name of the column reporting the age of each sample.

Value

data.frame with columns Core.id, F.WC (organic carbon sequestration rates at the whole core), A.Max (maximum age of the core), and seq_rate (average organic carbon sequestration rate at the indicated time frame)

Examples

bluecarbon_decompact <- decompact(bluecarbon_data)
#> Warning: Setting compaction = 0 for these cores: Sm_03_04, Sm_03_04, Sm_03_04, Sm_03_04, Sm_03_04, Sm_03_04, Sm_03_04, Sg_10_02, Sg_10_02, Sg_10_02, Sg_10_02, Sg_10_02, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sg_11_03, Sm_05_01, Sm_05_01, Sm_05_01, Sm_05_01, Sm_05_01, Sm_05_01, Sm_06_01, Sm_06_01, Sm_06_01, Sm_06_01, Sm_06_01, Sm_06_01
oc <- estimate_oc(bluecarbon_decompact)

#> Warning: The following cores had samples with organic carbon values below the organic carbon range used to built the model: Sg_04_01, Sm_04_03, Sm_04_04, Sm_05_01
#> Warning: The following cores had samples with organic carbon values above the organic carbon range used to built the model: Sg_04_01, Sm_03_01, Sm_04_02, Sm_04_03, Sm_04_04, Sm_05_01
out <- estimate_seq_rate(oc[[1]])