2.4. Groundwater Chemistry

A key point from the groundwater lectures is the connection between groundwater and surface water – that groundwater levels reflect changes in surface water. In the lecture we showed data for a hydrograph – that is the level of water in a river – and the underlying water table. There is a small delay, but the water table follows changes in the flow of the local river.

We have discussed the importance of groundwater systems, how groundwater moves (and why it can move uphill), and how we can test for and figure out the key aspects that influence groundwater flow, which includes pumping tests. Together these give us a good sense of how to calculate groundwater velocity and the energy required to move groundwater. There are several environmental challenges that face Earth’s groundwater systems, not least of which is the fact that we are removing groundwater faster than it can recharge.

Note

To explore the contamination of groundwater we are going to tell the story of arsenic poisoning in Bangladesh, why there are such high levels of arsenic in the groundwater and the various mechanisms that can help remedy this situation. This is just a specific example of an acute groundwater contamination problem and the way the community and scientists are trying to help remedy it. The details of this case study will be in the lecture slides.

How does Groundwater Chemistry Evolve?

To explore this, we consider the case of arsenic poisoning in groundwater wells in Nepal. There is tremendous elevation change in Nepal, and this leads to temperature gradients from the high Himalaya to the low plains, which in turn leads to nearly an order-of-magnitude difference in annual rainfall. In Nepal, groundwater supplies around 50% of the drinking water for the capital, Kathmandu, and 90% of the drinking water for the bulk of the population in the flood plains (Terai). Only 34% of the population has access to clean drinking water.

Fig. 2.28 Arsenic distribution in groundwater in Nepal Terai districts, showing the widespread contamination of shallow groundwater wells.

We can consider what happens to the chemistry of groundwater along the flow path to understand the chemistry of the drinking water in the flood plains (Fig. 2.29).

Fig. 2.29 Schematic showing how groundwater chemistry evolves along the flow path from recharge to discharge zones.

Chemical Evolution Along the Flow Path

How does the chemistry of groundwater change over distance?

1. Precipitation and Atmospheric Equilibration

   Rain is the major source of water to groundwater systems. Rain is in equilibrium with CO₂ in the atmosphere, and is slightly acidic (pH ~5.5-6.5).

2. Infiltration and Vadose Zone Processes

   During infiltration, chemical reactions occur with soil minerals due to leaching of organic carbon, and there is evapotranspiration which leads to concentration of ions. The water that infiltrates through the vadose zone is therefore slightly acidic and with a higher concentration of ions and dissolved organic matter (DOM).

   

   Fig. 2.30 Schematic of the vadose zone showing plant roots, soil layers, and the processes affecting infiltrating water.

   

   Fig. 2.31 Geochemical processes in the vadose zone, including organic matter leaching, ion exchange, and mineral weathering.

3. Mineral Dissolution

   Along the groundwater flow path there is further dissolution of minerals which contribute ions to the groundwater. Ionic strength therefore increases along the flow path. Typical ions that are picked up are shown in Fig. 2.32.

   

   Fig. 2.32 Common ions found in groundwater, including major cations (Ca²⁺, Mg²⁺, Na⁺, K⁺) and anions (HCO₃^-, SO₄^2-, Cl^-).

   Key to this is an understanding of saturation state and where minerals are undersaturated or oversaturated with respect to the fluids in the groundwater system. We express this using a solubility product for a mineral or salt.

4. Secondary Mineral Precipitation

   The precipitation of secondary minerals can occur as groundwater becomes supersaturated as it ages and changes its chemical composition.

   The longer the flow path, the more altered the chemistry of the fluid will be. In the deepest flow paths you can get exceptionally saline waters that form from dissolution of salt (Fig. 2.33).

   

   Fig. 2.33 Schematic summary of the distribution of Total Dissolved Solids (TDS) with depth in groundwater, showing increasing salinity with depth and residence time.

5. Development of Anoxic Conditions

   The longer groundwater is out of contact with the atmosphere, the more likely it is to become anoxic – and lose its oxygen. Remember that redox reactions are the ones that involve the exchange of electrons. In the subsurface there are many things that use up oxygen (respiration, reaction with minerals) but there is no way to make new oxygen.

This evolving chemistry becomes a real problem for arsenic. Arsenic-bearing minerals are easily weathered (chemically dissolved) in anoxic conditions, and in particular when there is a lot of dissolved organic matter which binds to the arsenate ions. Increased human waste at the surface is increasing the flux of organic matter into the subsurface and creating increasingly anoxic conditions that are accelerating the weathering of arsenic.

Annoxia, Redox and Arsenic

We can remember why groundwater becomes anoxic if we consider the electrode potential associated with redox half reactions. The electrode potential measures the voltage associated with the movement of electrons relative to a reference cell. Once you know the potential, you can calculate the total electrode potential (DE) by coupling two redox reactions – one that gives up electrons and one that accepts electrons, as shown in Fig. 2.34.

Fig. 2.34 The redox ladder showing the relationship between the standard reduction potential (E°) and the Gibbs free energy change (ΔG°) for key redox reactions.

Electrode potential is related to the overall energy of the redox reaction through

\[\Delta G = -n F \Delta E\]

Where \(n\) is the number of moles of electrons transferred in the redox reaction. Remember that the more negative the \(\Delta G\), the more energy for that overall reaction. If you have the conditions you can use the Nernst Equation to calculate the true, insitu \(\Delta G\) for a particular reaction.

What we can see is that the most energy is gained by reactions that are coupled to oxygen. Therefore, provided reactions are not kinetically limited, the reactions involving oxygen will go first. This means in groundwater systems, where there are many electron dense materials (dissolved organic matter, minerals), oxygen will react the fastest and therefore be depleted in the groundwater first.

We can image the change in the chemistry of groundwater using a Pourbaix diagram (Fig. 2.35), which shows the acid-base chemistry (pH) against the electron chemistry (pe) to determine which species are stable. As groundwater chemistry evolves you move gradually down on this diagram towards more reducing conditions, and whether you move left or right (more acidic or basic) depends on the minerals that might be dissolving or forming along the flow path.

Fig. 2.35 A Pourbaix diagram showing the stability of arsenic species as a function of pH and pe (electron activity). The diagram illustrates how arsenic speciation changes under different redox and pH conditions.