Dual clumped isotope thermometry

temp_D48() calculates carbonate growth temperature from D47 and D48 values.

Usage

temp_D48(
  D47_CDES90,
  D48_CDES90,
  D47_error,
  D48_error,
  ks,
  add = FALSE,
  col = "black",
  pch = 19
)

Arguments

D47_CDES90

Carbonate D47 values expressed on the CDES90 scale (‰).

D48_CDES90

Carbonate D48 values expressed on the CDES90 scale (‰).

D47_error

Error on the D47 value. Optional.

D48_error

Error on the D48 value. Optional.

ks

Kinetic slope. Has to be negative!

add

Add graphics to an already existing plot? Default: FALSE.

col

Graphical parameter. Optional.

pch

Graphical parameter. Optional.

Value

Returns the carbonate growth temperature (°C). If both D47_error and D48_error are specified temp_D48() returns a data frame.

Details

The function calculates a D47 value as an intersect of two curves: the equilibrium D47 vs D48 curve from Fiebig et al. (2021) and the kinetic slope. The resulting D47 value is then converted to temperature using the temp_D47() function and the equilibrium D47_CDES90 vs temperature equation of Fiebig et al. (2021).

Contributors

The source code of this function contains elements from the reconPlots package, available at https://github.com/andrewheiss/reconPlots

References

References are listed at D48() and D47().

See also

D47() calculates the equilibrium carbonate D47 value. D48() calculates the equilibrium carbonate D48 value.

Other thermometry: temp_D47(), temp_d18O()

Examples

temp_D48(0.617, 0.139, ks = -0.6)
#> [1] 38
temp_D48(0.546, 0.277, ks = -1)
#> [1] 30