psi: Molar ratio of CO2 released vs CaCO3 precipitated

In seacarb: Calculates parameters of the seawater carbonate system

Description

Returns the molar ratio of CO2 released vs CaCO3 precipitated described by Frankignoulle et al. (1994).

Usage

psi(flag, var1, var2, S=35, T=20, P=0, Pt=0, Sit=0, pHscale="T", kf="pf", k1k2="l")

Arguments

flag	select the couple of variables available. The flags which can be used are: flag = 1 pH and CO2 given flag = 2 CO2 and HCO3 given flag = 3 CO2 and CO3 given flag = 4 CO2 and ALK given flag = 5 CO2 and DIC given flag = 6 pH and HCO3 given flag = 7 pH and CO3 given flag = 8 pH and ALK given flag = 9 pH and DIC given flag = 10 HCO3 and CO3 given flag = 11 HCO3 and ALK given flag = 12 HCO3 and DIC given flag = 13 CO3 and ALK given flag = 14 CO3 and DIC given flag = 15 ALK and DIC given flag = 21 pCO2 and pH given flag = 22 pCO2 and HCO3 given flag = 23 pCO2 and CO3 given flag = 24 pCO2 and ALK given flag = 25 pCO2 and DIC given
var1	enter value of the first variable in mol/kg, except for pH and for pCO2 in μatm
var2	enter value of the second variable in mol/kg, except for pH
S	Salinity
T	Temperature in degrees Celsius
P	Hydrostatic pressure in bar (surface = 0)
Pt	Concentration of total phosphate in mol/kg
Sit	Concentration of total silicate in mol/kg
pHscale	choice of pH scale: "T" for the total scale, "F" for the free scale and "SWS" for using the seawater scale, default is "T" (total scale)
kf	"pf" for using Kf from Perez and Fraga (1987) and "dg" for using Kf from Dickson and Goyet (1979), default is "pf"
k1k2	"l" for using K1 and K2 from Lueker et al. and "r" for using K1 and K2 from Roy et al. , default is "l"

Details

Note that the arguments can be given as a unique number or as vectors. If the lengths of the vectors are different, the longer vector is retained and only the first value of the other vectors is used. It can therefore be critical to use vectors of the same length. For instance, to compute parameters from one couple of variable for a range of temperatures, a vector with temperatures required can be given in enter and other arguments can be completed be one variable this variable will be used for each temperatures.

Value

The function returns a data frame containing the following columns:

psi	ratio of CO2 released vs CaCO3 precipitated (mol/mol)

Author(s)

Jean-Pierre Gattuso [e-mail]

References

Frankignoulle, M. 1994 A complete set of buffer factors for acid/base CO2 system in seawater. Journal of Marine Systems 5, 111-118.

Frankignoulle, M., Canon, C. & Gattuso, J.-P., 1994. Marine calcification as a source of carbon dioxide- Positive feedback of increasing atmospheric CO2. Limnology and Oceanography 2, 458-462.

See Also

speciation.

Examples

## Calculation using the numerical example given in Frankignoulle et al. (1994)
psi(flag=24, var1=350, var2=2400e-6, S=35, T=25, P=0, Pt=0, Sit=0, pHscale="T", kf="pf", k1k2="l")

Example output

Loading required package: oce
Loading required package: gsw
Loading required package: testthat
[1] 0.5962306
attr(,"unit")
[1] "mol CO2/mol CaCO3"
attr(,"pH scale")
[1] "total scale"
attr(,"method")
[1] "Luecker et al. (2000)"

seacarb documentation built on May 2, 2019, 5:34 p.m.