speciation: ionic forms as a function of pH

Description Usage Arguments Value Author(s) References See Also Examples

Description

Estimates the concentration of the various ionic forms of a molecule as a function of pH

Usage

1
speciation(K1=K1(), K2=NULL, K3=NULL, pH, conc=1)

Arguments

K1

First dissociation constant

K2

Second dissociation constant, default is NULL

K3

Third dissociation constant, default is NULL

pH

pH value, default is 8

conc

concentration of molecule in mol/kg, default is 1 mol/kg

Value

The function returns a data frame containing the following concentrations (in mol/kg if conc is given in mol/kg):

C1

ionic form 1, univalent, bivalent and trivalent molecules

C2

ionic form 2, univalent, bivalent and trivalent molecules

C3

ionic form 3, bivalent and trivalent molecules

C4

ionic form 4, trivalent molecules

Author(s)

Karline Soetaert K.Soetaert@nioo.knaw.nl

References

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

See Also

bjerrum.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
## Speciation of divalent species; example to estimate the various ionic forms
## of dissolved inorganic carbon  (DIC = 0.0021 mol/kg) at a salinity of 35,
## a temperature of 25oC and an hydrostatic pressure of 0:
spec <- speciation (K1(35, 25, 0), K2(35, 25, 0), pH=8, conc=0.0021)
## where (spec\$C1=[CO2], spec\$C2=[HCO3-], spec\$C3=[CO3--])

## Speciation of trivalent species (e.g.,  H3PO4, H2PO4-, HPO4--, PO4---)
speciation(K1p(), K2p(), K3p(), conc=0.001)

## Effect of temperature on pCO2 - Figure 1.4.18 of Zeebe and Wolf-Gladrow (2001)
Tseq <- seq(0, 30, by=0.5)
pHseq <- carb(flag=15, var1=2300e-6, var2=1900e-6, S=35, T=Tseq, P=0)$pH
CO2  <- speciation(K1(T=Tseq), K2(T=Tseq), conc=1900, pH=pHseq)$C1
pCO2 <- CO2/Kh(T=Tseq)
plot(Tseq, pCO2, xlab="Temperature (oC)", ylab="pCO2 (uatm)", type="l", 
	main="effect of temperature on pCO2")
legend("topleft", c(expression(sum(CO[2])==1900~umol~kg^"-1"), 
	expression(TA==2300~umol~kg^"-1")))

seacarb documentation built on May 31, 2017, 3:42 a.m.

Search within the seacarb package
Search all R packages, documentation and source code