carb {seacarb}R Documentation

Parameters of the seawater carbonate system

Description

Returns parameters of the seawater carbonate system.

Usage

carb(df)

Arguments

df a data frame containing the input parameters (see below)

Details

The input data frame should contain the following variables in the right order:

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 pH and pCO2 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
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)
k1k2
"r" for using K1 and K2 from Roy et al. and "m" for using K1 and K2 from Mehrbach,
phflag
The pH scale is either the total scale (0) or the free scale (1), default is 0

Note that the results output is NO LONGER stored in a file.

Value

The function returns a dataframe containing the following columns:

S Salinity
T Temperature in degrees Celsius
P Pressure in bar
pH pH
CO2 CO2 concentration (mol/kg)
pCO2 pCO2, CO2 partial pressure (μatm)
fCO2 fCO2, CO2 fugacity (μ atm)
HCO3 HCO3 concentration (mol/kg)
CO3 CO3 concentration (mol/kg)
DIC DIC concentration (mol/kg)
ALK ALK, total alkalinity (mol/kg)
Oa Omega aragonite, aragonite saturation state
Oc Omega calcite, calcite saturation state
PhiD PhiD, chemical buffer factor (dpH/d[DIC]); input/output of dissolved CO2 (unit pH per mol/kg)
BetaD BetaD, homogeneous or Revelle buffer factor (dln(pCO2)/dln[DIC]); input/output of dissolved CO2
PiD PiD, chemical buffer factor (dpCO2/d[DIC]); input/output of dissolved CO2 (μ atm per mol/kg)
PhiB PhiB, chemical buffer factor (dpH/d[DIC]); from input/output of bicarbonate (unit pH per mol/kg)
BetaB BetaB, homogeneous buffer factor (dln(pCO2)/dln[DIC]); input/output of bicarbonate
PiB PiB, chemical buffer factor (dpCO2/d[DIC]); input/output of dissolved CO2 (μ atm per mol/kg)
PhiC PhiC, chemical buffer factor (dpH/d[DIC]); input/output of carbonate (unit pH per mol/kg)
BetaC BetaC, homogeneous buffer factor (dln(pCO2)/dln[DIC]); input/output of carbonate
PiC PiC, chemical buffer factor (dpCO2/d[DIC]); input/output of carbonate (μ atm per mol/kg)
PhiH PhiH, chemical buffer factor (dpH/d[ALK]); input/output of strong acid (unit pH per mol/kg)
PiH PiH, chemical buffer factor (dpCO2/d[ALK]); input/output of strong acid (μ atm per mol/kg)

Author(s)

Aur'elien Proye and Jean-Pierre Gattuso gattuso@obs-vlfr.fr

References

DOE, 1994 Handbook of methods for the analysis of the various parameters of the carbon dioxide system in sea water. ORNL/CDIAC-74. Oak Ridge,Tenn.: Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory.

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

Zeebe, R. E. and Wolf-Gladrow, D. A., 2001 CO2 in seawater: equilibrium, kinetics, isotopes. Amsterdam: Elsevier, 346 pp.

Examples

## Using a data set provided with the package:
data(seacarb_test)
carb(seacarb_test)

## Using another data set
x1 <- data.frame(flag=8,var1=8.2,var2=2400e-6,S=35,T=25,P=0,k1k2='r',phflag=0)
x2 <- data.frame(flag=8,var1=8.2,var2=2400e-6,S=30,T=30,P=0,k1k2='r',phflag=0)
x <- rbind(x1,x2)
carb(x)

[Package seacarb version 1.2.3 Index]