Graphe de Gran
# Graphe de Gran

Un graphe de Gran (aussi méthode de Gran et titrage de Gran) est un moyen usuel d'étalonner un titrate ou un titrant en estimant le volume d'équivalence pendant un titrage d'un acide fort avec une base forte ou pendant un titrage potentiométrique. De tels graphes peuvent aussi servir à estimer le contenu en carbonate d'une solution aqueuse, ainsi que les constantes d'acidité $K_a$ d'acides ou bases faibles à partir de données de titrage.

Les graphes de Gran mettent en œuvre des approximations linéaires de relations a priori non-linéaires entre la quantité mesurée, le pH ou le potentiel électromoteur (f.e.m.), et le volume de titrant ajouté. D'autres types de mesures de concentration, tels que l'absorbance spectrophotométrique ou les déplacements chimiques en RMN, peuvent en principe aussi servir. Ces approximations ne sont valides que près du point d'équivalence, mais non pas au point d'équivalence, et donc la méthode diffère des autres estimations de point d'équivalence utilisant les premiers ou seconds dérivés, méthodes qui requièrent des données au point d'équivalence. À l'origine, c'est-à-dire avant l'arrivée d'ordinateurs, les graphes de Gran ont été conçus comme détermination graphique, où un tracé sur papier quadrillé serait manuellement extrapolé pour estimer l'intersection sur l'abscisse (l'axe du volume). Ce procédé a été depuis remplacé par l'analyse plus précise par la méthode des moindres carrés, surtout à l'aide de logiciels de feuilles de calcul qui possèdent une fonction intégrée appliquant la méthode des moindres carrées.

## Titrage d'acide fort par base forte

Lors d'un titrage pH-métrique d'une solution d'acide fort, nous avons à chaque ième point,

    $K_{w^}=[H^{+}]_{i}[OH^{-}]_{i}$

où $K_w$ est la constante de autoprotolyse de l'eau.

Si le volume de départ de l'acide est $v_{0^}$ et sa concentration $[H^{+}]_{0^}$, et si le titrage utilise une base de concentration $[OH^{-}]_{0^}$, à tout ième point pendant le titrage avec un volume cumulé de titrant $v_{i^}$,

${\frac {v_{0}[H^{+}]_{0}-v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [H^{+}]_{i}{\text{ ou }}10^{-pH_{i}}&{\text{ si }}v_{0^}[H^{+}]_{0}>v_{i}[OH^{-}]_{0}\mathrm {~(r{\acute {e}}gion~acide)} \\=0&{\text{ si }}v_{0^}[H^{+}]_{0}=v_{i}[OH^{-}]_{0}\mathrm {~(point~d'{\acute {e}}quivalence)} \\\approx -[OH^{-}]_{i}{\text{ ou }}-K_{w}10^{pH_{i}}&{\text{ si }}v_{0^}[H^{+}]_{0}<v_{i}[OH^{-}]_{0}\mathrm {~(r{\acute {e}}gion~alcaline)} \end{cases}}$

Au point d'équivalence, le volume d'équivalence $v_{e^}=v_{i^}$.

Or,

| ${\color {Blue}\bullet }$ un graphe de | $({v_{0}+v_{i}})10^{-pH_{i}}$ vs. $v_{i^}$ | aura une région linéaire | avant | l'équivalence, avec pente | $-[OH^{-}]_{0^}$ |
| - | - | - | - | - | - |
| ${\color {Blue}\bullet }$ et un graphe de | $({v_{0}+v_{i}})10^{pH_{i}}$ vs. $v_{i^}$ | aura une région linéaire | après | l'équivalence, avec pente | $[OH^{-}]_{0}/K_{w^}$ |
| ${\color {Blue}\bullet }$ les deux graphes auront | $v_{e}=v_{0}[H^{+}]_{0}/[OH^{-}]_{0^}$ | comme intersection de l'abscisse |

Le volume d'équivalence $v_{e}$ est utilisé pour calculer l'inconnu, soit $[H^{+}]_{0^}$ ou $[OH^{-}]_{0^}$.

Avant d'entamer le titrage, le pH-mètre est d'habitude calibré à l'aide de solutions tampons aux pH connus. La puissance ionique peut être maintenue par un choix judicieux de l'acide et de la base. Par exemple, un titrage de l'acide HCl avec la base NaOH d'à-peu-près la même concentration remplacera $H^+$ par un ion de même charge $(Na^+)$ à plus ou moins la même concentration, pour maintenir la puissance ionique assez constante. Autrement, on peut utiliser une concentration relativement élevée d'électrolyte de fond, ou le coefficient d'activité peut être calculé.

## Titrage de base forte par acide fort

Des graphes illustrants les précédents sont obtenus avec un titrage de base forte par acide fort, et les pentes sont inversées.

${\frac {v_{0}[OH^{-}]_{0}-v_{i}[H^{+}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [OH^{-}]_{i}{\text{ ou }}K_{w}10^{pH_{i}}&{\text{ si }}v_{0^}[OH^{-}]_{0}>v_{i}[H^{+}]_{0}\mathrm {~(r{\acute {e}}gion~alcaline)} \\=0&{\text{ si }}v_{0^}[OH^{-}]_{0}=v_{i}[H^{+}]_{0}\mathrm {~(point~d'{\acute {e}}quivalence)} \\\approx -[H^{+}]_{i}{\text{ ou }}-10^{-pH_{i}}&{\text{ si }}v_{0^}[OH^{-}]_{0}<v_{i}[H^{+}]_{0}\mathrm {~(r{\acute {e}}gion~acide)} \end{cases}}$

Donc,

${\color {Blue}\bullet }$ $({v_{0}+v_{i}})10^{pH_{i}}$ vs. $v_{i^}$ aura une région linéaire avant l'équivalence avec pente $-[H^{+}]_{0}/K_{w^}$

${\color {Blue}\bullet }$ et $({v_{0}+v_{i}})10^{-pH_{i}}$ vs. $v_{i^}$ aura une région linéaire après l'équivalence avec pente $[H^{+}]_{0^}$

${\color {Blue}\bullet }$ les deux graphes auront $v_{e}=v_{0}[OH^{-}]_{0^}/[H^{+}]_{0}$ comme intersection de l'abscisse

La Figure 1 présente des graphes de Gran à la suite d'un titrage d'une base forte par un acide fort.

## Concentrations et constantes de dissociation d'acides faibles

La méthode peut estimer les constantes de dissociation d'acides faibles, ainsi que leurs concentrations (Gran, 1952). Si l'acide est représenté par HA, où

    $K_{a}={\frac {[H^{+}]_{i}[A^{-}]_{i}}{[HA]_{i}}}$,

on aura à chaque ième point d'un titrage d'un volume $v_{0}$ d'un acide à concentration $[HA]_{0}$ par une base à concentration $[OH^{-}]_{0}$ dans la région linéaire loin d'équivalence

    $[HA]_{i}\approx {\frac {v_{0}[HA]_{0}-v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}$ et
    $[A^{-}]_{i}\approx {\frac {v_{i}[OH^{-}]_{0}}{v_{0}+v_{i}}}$

qui sont des approximations valides, d'où

    $K_{a}\approx {\frac {10^{-pH_{i}}v_{i}[OH^{-}]_{0}}{v_{0}[HA]_{0}-v_{i}[OH^{-}]_{0}}}$ , ou
    $K_{a}(v_{0}{\frac {[HA]_{0}}{[OH^{-}]_{0}}}-v_{i})\approx 10^{-pH_{i}}v_{i}$ ou, parce que $v_{e}=v_{0}{\frac {[HA]_{0}}{[OH^{-}]_{0}}}$ ,
    $K_{a}(v_{e}-v_{i})\approx 10^{-pH_{i}}v_{i}$.

Un graphe de $10^{-pH_{i}}v_{i}$ contre $v_{i^}$ aura comme pente $-K_{a^}$ sur la région acide linéaire, et une intersection avec l'abscisse $v_{e^}$, duquel soit $[HA]_{0^}$ ou $[OH^{-}]_{0^}$ peut être calculé. La région alcaline est traitée de la même manière que pour une titrage d'un acide fort par une base forte. La Figure 2 donne un exemple; ici, les deux intersections diffèrent d'à-peu-près 0,2 mL, mais ceci n'est qu'un petit écart (0,5 %), vu le gros volume d'équivalence.

Des équations semblables peuvent être dérivées pour un titrage d'une base faible par un acide fort. (Gran, 1952; Harris, 1998).

## Contenu en carbonate

Martell et Motekaitis (1992) ont utilisé les régions les plus linéaires pour leurs graphes, puis exploité la différence entre le volume d'équivalence mesuré du côté acide et celui mesuré du côté alcalin pour ainsi estimer la quantité de $CO_2$ absorbée par la solution alcaline. Ce processus est illustré dans la Figure 1. Dans ce cas, l'acide supplémentaire qui sert à neutraliser le carbonate, par double protonation, dans un volume $v_{0^}$ de titrate est $(v_{e}-v_{e}^{\prime })[H^{+}]_{0^}=2v_{0}[CO_{2}]_{0}$. Dans le cas contraire d'un titrage d'acide par de la base, la teneur en carbonate est calculée de façon similaire avec $(v_{e}^{\prime }-v_{e})[OH^{-}]_{0^}=2v_{e}^{\prime }[CO_{2}]_{0}$, où cette fois $v_{e}^{\prime }$ représente le volume d'équivalence mesuré du côté alcalin (d'après Martell et Motekaitis).

Quand le contenu total de $CO_2$ absorbé est important, tel que dans les eaux naturelles et les effluents alcalins, deux ou trois inflections peuvent être décelées dans la courbe pH-volume dues aux effets tampons qu'ont le bicarbonate et le carbonate en concentrations importantes. Comme l'ont discuté Stumm et Morgan (1981), l'analyse de telles eaux peut mettre en œuvre jusqu'à six graphes de Gran à partir d'une seule titrage pour estimer les multiples points d'équivalence et mesurer l'alcalinité totale et les concentrations de carbonate et/ou de bicarbonate.

## Contrôle potentiométrique de $H^+$

Pour utiliser des mesures potentiométriques $E_{i^}$ au lieu de $pH_{i^}$ pour contrôler la concentration de $H^{+_}$, on peut de façon triviale laisser $-log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^}$ et appliquer les équations ci-dessus, où $b_{0^}$ est la correction du décalage de la réponse potentiométrique $nFE_{0^}/RT$ et $b_{1^}$ est une correction de la pente de la réponse potentiométrique $nF^{_}/RT$ (1/59.2 unités de pH/mV à 25 °C), de sorte que $-b_{1}E_{i^}$ remplace $pH_{i^}$.

Ainsi, comme vu plus tôt pour un titrage d'acide fort par une base forte,

${\color {Blue}\bullet }$ $({v_{0}+v_{i}})10^{b_{1}E_{i}}$ vs. $v_{i^}$ aura une région linéaire avant l'équivalence, avec pente $-10^{b_{0}}[OH^{-}]_{0^}$

${\color {Blue}\bullet }$ $({v_{0}+v_{i}})10^{-b_{1}E_{i}}$ vs. $v_{i^}$ aura une région linéaire après l'équivalence, avec pente $10^{-b_{0}}[OH^{-}]_{0}/K_{w^}$

${\color {Blue}\bullet }$ les deux auront $v_{e}=v_{0}[H^{+}]_{0}/[OH^{-}]_{0^}$ comme intersection sur x et, comme auparavant, le volume d'équivalence provenant du côté acide servira à étalonner la concentration inconnue, et la différence entre les deux volumes d'équivalence servira à estimer la teneur en carbonate

Des graphiques analogues peuvent être construits avec les données d'un titrage d'une base par de l'acide.

## Calibrage d'électrode

L'analyse de données potentiométriques nécessite la connaissance préalable de $b_{0^}$ et de $b_{1^}$.

Si l'électrode est mal calibré, une correction de décalage peut être calculée in situ à partir de la pente en milieu acide:

* Pour un titrage d'un acide par de la base, la pente du côté acide est $-[OH^{-}]_{0^}10^{b}$, d'où $b_{^}$ peut être calculé à l'aide de la valeur connue de $[OH^{-}]_{0^}$ ou celle provenant du volume d'équivalence. $K_{w^}$ peut ensuite être calculé à partir de la pente du côté alcalin.
* Pour un titrage d'une base par de l'acide, telle qu'illustré dans les exemples, la pente du côté acide, $[H^{+}]_{0^}10^{b}$, sert pareillement à calculer $b_{0^}$, et la pente du côté alcalin, $-[H^{+}]_{0^}10^{-b_{0}}/K_{w}$ sert à calculer $K_{w^}$. Dans l'exemple de la Figure 1, le décalage n'était pas conséquent, car $b_{0^}$ s'est avéré être -0.054 unités de pH.

La valeur de $b_{1^}$, par contre, peut dévier de la valeur théorique et ne peut être évaluée que par le calibrage de l'électrode. Ceci est souvent effectué avec des solutions tampons aux pH connus, ou en titrant un acide fort avec une base forte, dans lequel cas une force ionique constante peut être maintenue, et $[H^{+_}]_{i}$ sera connu à tous les points du titrage si $[H_^{+}]_{0}$ et $[OH_^{-}]_{0}$ sont tous deux connus (et qui devraient être directement reliés à des étalons primaires). Par exemple, Martell et Motekaitis (1992) ont calculé le pH prévu au début du titrage, ayant au préalable titré les solutions d'acide et de base contre des étalons primaires, puis ajusté le pH ainsi relevé, mais ceci ne donne pas de correction de la pente si l'une est nécessaire.

Suivant le travail antérieur de McBryde (1969), Gans et O'Sullivan (2000) décrivent une approche itérative pour calculer $b_{0^}$ et $b_{1^}$ dans la relation $-log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^}$, à partir d'un titrage d'acide fort avec une base forte:

* $b_{0^}$ est d'abord estimé à l'aide des données acides selon Rossotti et Rossotti (1965), et $b_{1^}$ est en premier lieu donné sa valeur théorique;
* on trace des fonctions Gran modifiées, soit $({v_{0}+v_{i}})10^{b_{1}E_{i}-b_{0}}$ vs. $v_{i^}$ du côté acide de l'équivalence et $({v_{0}+v_{i}})K_{w}10^{-b_{1}E_{i}+b_{0}}$ vs. $v_{i^}$ du côté alcalin, et les volumes d'équivalence $v_{e^}$ et $v_{e}^{\prime }$ en sont calculés, comme auparavant;
* comme auparavant, la différence entre ces volumes d'équivalence servent à estimer le contenu en carbonate mais aussi pour calculer une 'concentration effective' de la base $v_{e}[OH^{-}]_{0^}/v_{e}^{\prime }$ du côté alcalin de l'équivalence;
* les valeurs approximatives de $[H^{+}]_{i^}$ sont calculées avec $(v_{0^}[H^{+}]_{0}-v_{i}[OH^{-}]_{0})/(v_{0}+v_{i})$ du côté acide et avec $(v_{e}[H^{+}]_{0}-v_{i}v_{e}[OH^{-}]_{0}/v_{e}^{\prime })/(v_{0}+v_{i})$ du côté alcalin;
* la définition $-log_{10}[H^{+}]_{i}=b_{0}-b_{1}E_{i^}$ est réécrite $E_{i^}=b_{0}/b_{1}+(1/b_{1})log_{10}[H^{+}]_{i}$, et les données de potentiel $E_{i^}$ sont tracées contre $log_{10}[H^{+}]_{i^}$, en utilisant les valeurs de $[H^{+}]_{i^}$ qui correspondent aux valeurs de pH entre 2.5 et 4.5, et entre 10.7 et 11.5; la méthode linéaire des moindres carrés procure alors la pente $1/b_{1^}$ et l'intersection $b_{0}/b_{1^}$;
* les étapes 2 et 3 sont répétées avec les nouvelles valeurs de $b_{0^}$ et $b_{1^}$ pour atteindre une plus grande précision dans le calcul des volumes d'équivalence et du contenu en $CO_2.$

Cette procédure pourrait en principe être modifiée pour des titrages de base par de l'acide. Le logiciel GLEE (nom dérivé de 'GLass Electrode Evaluation') met en œuvre cette approche pour étalonnage d'un électrode à partir d'un titrage d'acide fort par une base forte. Ce logiciel peut aussi calculer (par un procédé non-linéaire de moindres carrés séparé) une « correction » de la concentration de la base. L'avantage de cette méthode est que le calibrage peut être conduit dans le même milieu à force ionique constante que l'on utilisera pour la détermination des constantes d'équilibre.

Il est à noter que les fonctions de Gran habituelles donnent aussi les volumes d'équivalence requis et que, comme $b_{1^}$ porte au départ sa valeur théorique, la valeur estimée de $b_{0^}$ à la première étape aurait pu provenir de la pente du côté acide, tel que détaillé plus tôt. Aussi, ce procédé calcule la teneur en $CO_2$ et peut donc être combiné avec un étalonnage complet de la base, en utilisant la définition de $v_{e^}$ pour obtenir $[OH^{-}]_{0^}$. Finalement, la plage de pH utilisée pourrait être étendue en trouvant $[H^{+}]_{i^}$ par la quadratique $(v_{0{^}}[H^{+}]_{0}-v_{i}[OH^{-}]_{0})/(v_{0}+v_{i})=[H^{+}]_{i}-K_{w}/[H^{+}]_{i}$.

## Contrôle potentiométrique d'autres espèces

Lorsque l'on contrôle une autre espèce, disons $S^{1_}$, par la potentiométrie, on peut appliquer le même formalisme avec $-log_{10}[S^{1}]_{i}=b_{0}-b_{1}E_{i^}$. Alors, un titrage d'une solution d'une seconde espèce $S^{2_}$ par $S^{1_}$ devient analogue à un titrage d'une base par de l'acide avec contrôle du pH, d'où un graphe soit de $({v_{0}+v_{i}})10^{b_{1}E_{i}}$ ou de $({v_{0}+v_{i}})10^{-b_{1}E_{i}}$ contre $v_{i^}$ aura un intersection sur x de $v_{0}[S^{0}]_{0}/[S^{1}]_{0^}$. Dans un titrage contraire de $S^{1_}$ par $S^{2_}$, le volume d'équivalence sera $v_{0}[S^{1}]_{0}/[S^{0}]_{0^}$. La signification des pentes dépendra alors du type d'interaction entre les deux espèces, qu'il s'agisse d'une association ou d'une précipitation (Gran, 1952), mais, d'habitude, le seul résultat désiré est le point d'équivalence. Toutefois, la pente qui précède le point d'équivalence peut en principe servir à mesurer le produit de solubilité $K_{sp^}$ de la même façon que l'on évalue le $K_{w^}$ avec un titrage acide-base, bien que d'autres interactions entre paires d'ions peuvent aussi opérer.

Pour illustrer ceci, prenons l'exemple d'un titrage potentiométrique de $Cl^−$ par $Ag^+:$ ${\frac {v_{0}[Cl^{-}]_{0}-v_{i}[Ag^{+}]_{0}}{v_{0}+v_{i}}}{\begin{cases}\approx [Cl^{-}]_{i}{\text{ ou }}K_{sp}10^{-b_{1}E_{i}+b_{0}}&{\text{ lorsque }}v_{0^}[Cl^{-}]_{0}>v_{i}[Ag^{+}]_{0}\mathrm {~(avant~l'{\acute {e}}quivalence)} \\=0&{\text{ lorsque }}v_{0^}[Cl^{-}]_{0}=v_{i}[Ag^{+}]_{0}\mathrm {~({\grave {a}}~l'{\acute {e}}quivalence)} \\\approx -[Ag^{+}]_{i}{\text{ ou }}-10^{b_{1}E_{i}-b_{0}}&{\text{ lorsque }}v_{0^}[Cl^{-}]_{0}<v_{i}[Ag^{+}]_{0}\mathrm {~(apr{\grave {e}}s~l'{\acute {e}}quivalence)} \end{cases}}$

Donc,

${\color {Blue}\bullet }$ $({v_{0}+v_{i}})10^{-b_{1}E_{i}}$ vs. $v_{i^}$ aura une région linéaire avant l'équivalence, avec pente $-[Ag^{+}]_{0}10^{-b_{0}}/K_{sp^}$

${\color {Blue}\bullet }$ et $({v_{0}+v_{i}})10^{b_{1}E_{i}}$ vs. $v_{i^}$ aura une région linéaire après l'équivalence, avec pente $[Ag^{+}]_{0^}10^{b_{0}}$

${\color {Blue}\bullet }$ dans les deux cas, l'intersection avec x sera $v_{e}=v_{0}[Cl^{-}]_{0^}/[Ag^{+}]_{0}$

La figure 3 donne un exemple de données potentiométriques d'un tel titrage.

## Comportement non-idéal

Dans tout titrage sans espèces tampons, les graphes avant et après l'équivalence devraient idéalement croiser l'abscisse au même point, mais un comportement non-idéal peut être occasionné par des erreurs de mesures (e.g. avec un électrode mal calibré, un temps d'équilibration insuffisant avant d'enregistrer le potentiel, une dérive de la force ionique), des erreurs d'échantillonnage (e.g. une densité de données faible dans les régions linéaires) ou par un modèle chimique incomplet (e.g. avec la présence d'impuretés titrables telle que le carbonate dans la base, ou une précipitation incomplète lors d'un titrage potentiométrique à cause de solutions trop diluées, pour lequel cas Gran et al. (1981) propose des approches alternatives). Buffle et al. (1972) discutent de plusieurs sources de telles erreurs.

Parce que les fonctions en $10^{pH_{i}}$ ou $10^{-pH_{i}}$ de Gran ne tendent qu'asymptotiquement vers zéro, sans jamais l'atteindre, la courbure qui résulte en approche de l'équivalence est attendue dans tous les cas. Toutefois, les praticiens ne s'accordent pas sur quelles données utiliser pour déterminer l'équivalence, qu'il s'agisse d'utiliser les données avant ou après l'équivalence ou des deux côtés, et qu'il s'agisse d'inclure les données les plus proches de l'équivalence ou seulement celles dans les portions les plus $linéaires^,:$ Les données les plus proches de l'équivalence feront en sorte que les deux estimations de l'intersection sur l'axe de x seront plus coïncidants et coïncideront mieux aussi aux estimés provenant des premiers ou seconds dérivés, alors que l'on peut présumer que les données à pH acide seront moins assujetties aux interférences provenant d'impuretés tamponneuses et titrables, telles que le bicarbonate et le carbonate dans la base (voir Contenu en carbonate), ainsi que de l'effet d'une force ionique dérivante. Dans les graphiques donnés ici en exemples, les régions les plus linéaires (les données représentées par les cercles emplis) étaient celles choisies pour le traitement par la méthode des moindres carrés donnant les pentes et intersections. La sélection des données à retenir est toujours subjective.

## Littérature citée

* Buffle, J., Parthasarathy, N. et Monnier, D. (1972): Errors in the Gran addition method. Part I. Theoretical Calculation of Statistical Errors; Anal. Chim. Acta 59, 427-438; Buffle, J. (1972): Anal. Chim. Acta 59, 439.
* Butler, J. N. (1991): Carbon Dioxide Equilibria and Their Applications; CRC Press: Boca Raton, FL.
* Butler, J. N. (1998): Ionic Equilibrium: Solubility and pH Calculations; Wiley-Interscience. Chap. 3.
* Gans, P. et O'Sullivan, B. (2000): GLEE, a new computer program for glass electrode calibration; Talanta, 51, 33–37.
* Gran, G. (1950): Determination of the equivalence point in potentiometric titrages, Acta Chem. Scand., 4, 559-577.
* Gran, G. (1952): Determination of the equivalence point in potentiometric titrages-- Part II, Analyst, 77, 661-671.
* Gran, G., Johansson, A. et Johansson, S. (1981): Automatic titrage by Stepwise Addition of Equal Volumes of Titrant Part VII. Potentiometric Precipitation titrages, Analyst, 106, 1109-1118.
* Harris, D. C.: Quantitative Chemical Analysis, 5th Ed.; W.H. Freeman & Co., New. York, NY, 1998.
* Martell, A. E. et Motekaitis, R. J.: The determination and use of stability constants, Wiley-VCH, 1992.
* McBryde, W. A. E. (1969): Analyst, 94, 337.
* Rossotti, F. J. C. et Rossotti, H. (1965): J. Chem. Ed., 42, 375
* Skoog, D. A., West, D. M., Holler, F. J. et Crouch, S. R. (2003): Fundamentals of Analytical Chemistry: An Introduction, 8th Ed., Brooks and Cole, Chap. 37.
* Stumm, W. et Morgan, J. J. (1981): Aquatic chemistry, 2nd Ed.; John Wiley & Sons, New York.
