carb {seacarb} | R Documentation |
Returns parameters of the seawater carbonate system.
carb(df)
df |
a data frame containing the input parameters (see below) |
The input data frame should contain the following variables in the right order:
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
Note that the results output is NO LONGER stored in a file.
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) |
Aur'elien Proye and Jean-Pierre Gattuso gattuso@obs-vlfr.fr
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.
## 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)