Bjerrum plot

Example Bjerrum plot: Change in carbonate system of seawater from ocean acidification.

A Bjerrum plot (named after Niels Bjerrum), sometimes also known as a Sillén diagram (after Lars Gunnar Sillén), or a Hägg diagram (after Gunnar Hägg) is a graph of the concentrations of the different species of a polyprotic acid in a solution, as a function of pH, when the solution is at equilibrium. Due to the many orders of magnitude spanned by the concentrations, they are commonly plotted on a logarithmic scale. Sometimes the ratios of the concentrations are plotted rather than the actual concentrations. Occasionally H⁺ and OH⁻ are also plotted.

Most often, the carbonate system is plotted, where the polyprotic acid is carbonic acid (a diprotic acid), and the different species are dissolved carbon dioxide, carbonic acid, bicarbonate, and carbonate. In acidic conditions, the dominant form is Template:CO2; in basic (alkaline) conditions, the dominant form is Template:Chem; and in between, the dominant form is Template:Chem. At every pH, the concentration of carbonic acid is assumed to be negligible compared to the concentration of dissolved Template:Chem, and so is often omitted from Bjerrum plots. These plots are very helpful in solution chemistry and natural water chemistry. In the example given here, it illustrates the response of seawater pH and carbonate speciation due to the input of man-made Template:Chem emission by the fossil fuel combustion.

The Bjerrum plots for other polyprotic acids, including silicic, boric, sulfuric and phosphoric acids, are other commonly used examples.

Bjerrum plot equations for carbonate system

If carbon dioxide, carbonic acid, hydrogen ions, bicarbonate and carbonate are all dissolved in water, and at chemical equilibrium, their equilibrium concentrations are often assumed to be given by:

$\begin{array}{cc}[]{\left[{\text{CO}}_2\right]}_{\text{eq}}& =\frac{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2+K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}+K_1{K}_2}\times \text{DIC},\\ {\left[{\text{HCO}}_3^{-}\right]}_{\text{eq}}& =\frac{K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2+K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}+K_1{K}_2}\times \text{DIC},\\ {\left[{\text{CO}}_3^{2-}\right]}_{\text{eq}}& =\frac{K_1{K}_2}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2+K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}+K_1{K}_2}\times \text{DIC},\end{array}$

where the subscript 'eq' denotes that these are equilibrium concentrations, K₁ is the equilibrium constant for the reaction Template:Chem + Template:Chem Template:Eqm H⁺ + Template:Chem (i.e. the first acid dissociation constant for carbonic acid), K₂ is the equilibrium constant for the reaction Template:Chem Template:Eqm H⁺ + Template:Chem (i.e. the second acid dissociation constant for carbonic acid), and DIC is the (unchanging) total concentration of dissolved inorganic carbon in the system, i.e. [[[:Template:CO2]]] + [[[:Template:Chem]]] + [[[:Template:Chem]]]. K₁, K₂ and DIC each have units of a concentration, e.g. mol/L.

A Bjerrum plot is obtained by using these three equations to plot these three species against Template:Nowrap, for given K₁, K₂ and DIC. The fractions in these equations give the three species' relative proportions, and so if DIC is unknown, or the actual concentrations are unimportant, these proportions may be plotted instead.

These three equations show that the curves for Template:CO2 and Template:Chem intersect at Template:Nowrap, and the curves for Template:Chem and Template:Chem intersect at Template:Nowrap. Therefore, the values of K₁ and K₂ that were used to create a given Bjerrum plot can easily be found from that plot, by reading off the concentrations at these points of intersection. An example with linear Y axis is shown in the accompanying graph. The values of K₁ and K₂, and therefore the curves in the Bjerrum plot, vary substantially with temperature and salinity.

Chemical and mathematical derivation of Bjerrum plot equations for carbonate system

Suppose that the reactions between carbon dioxide, hydrogen ions, bicarbonate and carbonate ions, all dissolved in water, are as follows: Template:NumBlk Template:NumBlk

Note that reaction Template:EquationNote is actually the combination of two elementary reactions:

Template:Chem + Template:Chem Template:Eqm Template:Chem Template:Eqm H⁺ + Template:Chem

Assuming the mass action law applies to these two reactions, that water is abundant, and that the different chemical species are always well-mixed, their rate equations are

$\begin{array}{cc}\frac{\text{d}\left[{\text{CO}}_2\right]}{\text{d}t}& =-k_1\left[{\text{CO}}_2\right]+k_{-1}\left[{\text{H}}^{+}\right]\left[{\text{HCO}}_3^{-}\right],\\ \frac{\text{d}\left[{\text{H}}^{+}\right]}{\text{d}t}& =k_1\left[{\text{CO}}_2\right]-k_{-1}\left[{\text{H}}^{+}\right]\left[{\text{HCO}}_3^{-}\right]+k_2\left[{\text{HCO}}_3^{-}\right]-k_{-2}\left[{\text{H}}^{+}\right]\left[{\text{CO}}_3^{2-}\right],\\ \frac{\text{d}\left[{\text{HCO}}_3^{-}\right]}{\text{d}t}& =k_1\left[{\text{CO}}_2\right]-k_{-1}\left[{\text{H}}^{+}\right]\left[{\text{HCO}}_3^{-}\right]-k_2\left[{\text{HCO}}_3^{-}\right]+k_{-2}\left[{\text{H}}^{+}\right]\left[{\text{CO}}_3^{2-}\right],\\ \frac{\text{d}\left[{\text{CO}}_3^{2-}\right]}{\text{d}t}& =k_2\left[{\text{HCO}}_3^{-}\right]-k_{-2}\left[{\text{H}}^{+}\right]\left[{\text{CO}}_3^{2-}\right]\end{array}$

where Template:Nowrap denotes concentration, t is time, and K₁ and k₋₁ are appropriate proportionality constants for reaction Template:EquationNote, called respectively the forwards and reverse rate constants for this reaction. (Similarly K₂ and k₋₂ for reaction Template:EquationNote.)

Template:Em, the concentrations are unchanging, hence the left hand sides of these equations are zero. Then, from the first of these four equations, the ratio of reaction Template:EquationNote's rate constants equals the ratio of its equilibrium concentrations, and this ratio, called K₁, is called the equilibrium constant for reaction Template:EquationNote, i.e. Template:NumBlk

where the subscript 'eq' denotes that these are equilibrium concentrations.

Similarly, from the fourth equation for the equilibrium constant K₂ for reaction Template:EquationNote, Template:NumBlk

Rearranging Template:EquationNote gives Template:NumBlk

and rearranging Template:EquationNote, then substituting in Template:EquationNote, gives Template:NumBlk

The total concentration of dissolved inorganic carbon in the system is given by substituting in Template:EquationNote and Template:EquationNote:

$\begin{array}{cc}\text{DIC}& =\left[{\text{CO}}_2\right]+\left[{\text{HCO}}_3^{-}\right]+\left[{\text{CO}}_3^{2-}\right]\\ & ={\left[{\text{CO}}_2\right]}_{\text{eq}}(1+\frac{K_1}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}}+\frac{K_1{K}_2}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2})\\ & ={\left[{\text{CO}}_2\right]}_{\text{eq}}\left(\frac{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2+K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}+K_1{K}_2}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2}\right)\end{array}$

Re-arranging this gives the equation for Template:Chem:

${\left[{\text{CO}}_2\right]}_{\text{eq}}=\frac{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2}{{\left[{\text{H}}^{+}\right]}_{\text{eq}}^2+K_1{\left[{\text{H}}^{+}\right]}_{\text{eq}}+K_1{K}_2}\times \text{DIC}$

The equations for Template:Chem and Template:Chem are obtained by substituting this into Template:EquationNote and Template:EquationNote.

See also

• Charlot equation
• Gran plot (also known as Gran titration or the Gran method)
• Henderson–Hasselbalch equation
• Hill equation (biochemistry)
• Ion speciation
• Fresh water
• Seawater
• Thermohaline circulation