Electrical Conductivity (EC)
Three Calculation Methods
The electrical conductivity (EC) or specific conductance is a useful water-quality parameter. There are several methods available to calucalate EC.1 Three of them are used by the program:
- Linear approach (proportional to ionic strength)
- Pseudo-linear approach (Inverse Marion & Babcock)
- Diffusion coefficient-based approach (Appelo 2010)2 ⇐ default method
The first two approaches are simple empirical methods based on the ionic strength. The third approach is more advanced and relies on diffusion coefficients; it is the default algorithm of the program.3
For any given aqueous solution two values are displayed: (i) the calculated EC at the actual water temperature and (ii) EC25 after conversion to the reference temperature of 25°C. In this way, measurements or samplings taken at different T can be compared.
Method 1: Linear Approach based on Ionic Strength
The simplest empirical method relies on a linear relationship between electrical conductivity and the ionic strength I:
|(1)||EC (µS/cm) = 6.2 104 × I (mol/L)|
This equation is equivalent to the common approximation (i.e. the inverse of Eq.(1)):
|(2)||I (mol/L) = 1.6 10-5 × EC (µS/cm)|
People use it to get a rough number for the ionic strength from easy-to-measure EC values. But we go the other way round: The ionic strength I that enters Eq.(1) is strictly determined by the actual water composition/speciation as:
where the sum runs over all ions i with molar concentration ci and charge number zi. In hydrochemistry, the ionic strength I is calculated anyway because it enters the activity model to account for ion-ion interactions in non-ideal solutions.4
TDS. There is also a simple linear relationship between EC value and TDS – see here.
Method 2: Pseudo-linear Approach (Inverse Marion & Babcock)
|(4)||log10 I = 1.159 + 1.009 log10 EC||for I ≤ 0.3 mol/L|
In this equation the units of I are mmol/L (= mM) and the units of EC are dS/m – which differ significantly from the units in Eq.(1). The rearrangement of Eq.(4) into a form similar to Eq.(2) can be done step by step:
|log10 (EC / dS∙m-1)||=||0.991 log10 (I / mM) - 1.149|
|log10 (10-3 EC / µS∙cm-1)||=||0.991 log10 (103 I/M) - 1.149|
|log10 (EC / µS∙cm-1) - 3||=||0.991 [ 3 + log10 (I/M) ] - 1.149|
|log10 (EC / µS∙cm-1)||=||4.824 + log10 (I/M)|
|(5)||EC (µS/cm) = 6.67 104 × [ I (mol/L) ] 0.991|
Because I0.991 ≈ I, this equation is very similar to Eq.(1). Thus, this approach is called ‘pseudo-linear’. Due to a slightly larger prefactor in Eq.(5), the EC of the pseudo-linear approach exceeds the EC of the linear approach a little.
[The alternative name ‘inverse Marion & Babcock’ results from the fact that we rearranged the original equation from I = I(EC) to EC = EC(I).]
Method 3: Approach based on Diffusion Coefficients
The Nernst-Einstein equation establishes a physical relationship between the molar limiting conductivity and the diffusion coefficient Di for a given ion i. Appelo2 adopted this idea for practical use:
The meaning of the symbols and the mathematical derivation of this formula is described here. This approach is used in aqion as the default algorithm.
Temperature Compensation: EC ⇒ EC25
The EC of most natural waters, including seawater, increases with temperature 1-3% per degree Celsius. Measured EC values are usually referred to 25°C – often indicated by EC25. For this purpose the program converts the calculated EC (valid for the given water temperature T) to EC25 at 25°C. Principally, there are two main approaches: (i) a nonlinear model and (ii) its linear approximation.
Nonlinear T Compensation. The nonlinear model is an outcome of the physical relationship between electrical conductivity, diffusion coefficients, and the viscosity of water. The equation is given by:
|(7)||EC25 = 1.125 · 10-A/B · EC|
with the two parameters taken from Atkins:7
|(7a)||A = 1.37023 (T-20) + 8.36·10-4 (T-20)2|
|(7b)||B = 109 + T and T in °C|
The nonlinear compensation model is the standard method used in aqion.
Linear Approximation. Instead of the general approach in Eq.(7), linear formulas are in widespread use. The most common type of a linear expression is obtained from Eq.(7) by Taylor series expansion (as shown here):
|(8)||EC25 = EC / [ 1 + a (T - 25) ]|
with a = 0.020 °C-1 and T in °C.
Program Output. The program displays both values, the calculated EC (based on diffusion coefficients) and the compensated value EC25. This is done in the output tables:
|EC_25||when checkbox Mol is on|
|EC (T)||when checkbox Mol is off|
Conversion of Units
The physical units of EC take a little getting used to. The conversions between µS/cm (micro Siemens per centimeter) and other EC units are:
|1 mS/m||=||10 µS/cm|
|1 dS/m||=||1000 µS/cm|
|1 dS/m||=||1 mS/cm|
|1 µmho/cm||=||1 µS/cm|
where 1 S = 1 Siemens = 1 ohm-1 = 1 mho. The program uses µS/cm as default unit for EC.
Typical Conductivities of Aqueous Solutions
|absolute pure water||0.055||µS/cm|
|rain water||5 – 30||µS/cm|
|potable water||500 – 1000||µS/cm|
|groundwater||30 – 2000||µS/cm|
|industrial wastewater||≥ 5000||µS/cm|
|concentrated acids and bases||up to 1 000 000||µS/cm|
Pure Water. Due to the self ionization of water into H+ und OH- ions, the electrical conductivity of pure water is non-zero: EC = 0.055 µS/cm at 25°C.
R.B. McCleskey, D.K. Nordstrom, and J.N. Ryan: Comparison of electrical conductivity calculation methods for natural waters, Limnol. Oceanogr.: Methods 10, 952–967 (2012) ↩
To change or select the calculation method for EC click on Settings in the upper menu of aqion. ↩
Garrison Sposito: The Chemistry of Soils, 2nd Edition, Oxford University Press, 2008, (see Eq.(4.23) p.111) ↩
G.M. Marion, K.L. Babcock: Predicting specific conductance and salt concentration in dilute aqueous solutions. Soil Sci. 122, 181–187 (1976) ↩
P. Atkins and J. de Paula: Physical Chemistry, 8th Edition, W. H. Freeman and Company New York, 2006, Table 21.4, p. 1019 ↩