Theoretical progress of groundwater chemical evolution based on Tóthian theory: A review

Dong, Hongzhi; Gao, Zongjun

doi:10.3389/fmars.2022.972426

REVIEW article

Front. Mar. Sci., 22 August 2022

Sec. Coastal Ocean Processes

Volume 9 - 2022 | https://doi.org/10.3389/fmars.2022.972426

This article is part of the Research Topic Coastal Ecotones under Anthropogenic Influences: Conservation, Restoration and Realignment View all 14 articles

Theoretical progress of groundwater chemical evolution based on Tóthian theory: A review

Hongzhi Dong

Zongjun Gao^*

College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao, China

Tóthian theory refers to the gravity driven groundwater flow system (GFS) theory represented by Tóth, which mainly expounds the driving and distribution law of groundwater. The establishment and development of this theory not only deepened people’s understanding of the driving and distribution law of groundwater, but also greatly promoted the study of groundwater chemical evolution (GCE). Modern GCE research is mostly based on Tóthian theory, characterized by combining with advanced scientific and technological means. Based on the clue of time, this paper is divided into two parts. The first part mainly summarizes the establishment and development of Tóthian theory, including the exploration of groundwater driving force and distribution form by hydrogeologists before Tóthian theory, and the enrichment, development and application of Tóthian theory by geologists after its establishment. The second part mainly combs the main theories and application progress of GCE mechanism research, including the main theories and findings of GCE research before the emergence of Tóthian theory, as well as the research progresses of GCE after the emergence of Tóthian theory. With the flow of groundwater in GFS, groundwater undergoes continuous chemical evolution, which eventually leads to the transformation of hydrochemical types and the gradual increase of total dissolved solids (TDS). The distribution of GFS and GCE complement each other. The distribution of GFS directly determines the model of GCE, and the results of GCE also play a certain role in the distribution of GFS. GCE mainly includes dissolution, precipitation, cation exchange and adsorption, which is affected by the physical and chemical conditions of permeable media, organic matter content and microorganisms. GCE has the characteristics of universality, sustainability and diversity. With the increasing global population and the progresses of science and technology, the impact of human life, industrial and agricultural production on groundwater is deepening. The aggravation of pollution directly changes the chemical compositions of groundwater, resulting in changes of the law of GCE.

Introduction

Groundwater, a complex solution containing more than 80 elements, is an important geological force, which is widely involved in geological processes such as karst, sedimentation, diagenesis, metamorphism and mineralization (Zhang et al., 2018). Driven by gravity, density difference and heat, groundwater continuously reacts with the surrounding environment to complete the transportation and distribution of materials and energy, and its chemical properties continue to evolve (Tóth, 2009; Jiang et al., 2011; Zhang et al., 2015; Luan et al., 2017; Zhang et al., 2018; Meng, 2021; Zhao et al., 2021). According to Tóthian theory, the distribution pattern of gravity driven GFS is mainly controlled by the terrain. Regional GFS has multi-level nested systems, i.e., local-GFS, intermediate-GFS and regional-GFS. It is recognized that regional-GFS is the reason for the diversity of geological processes (Tóth, 1999; Jiang et al., 2010; Judit and Tóth, 2015; Jia et al., 2020). The discovery of Tóthian theory makes people master the general law of GFS. The combination of Tóthian theory and groundwater evolution research greatly promotes the study of GCE, which makes it possible to reveal the general law of GCE.

The long-term circulation makes the main components of groundwater controlled by the stratum. Under natural conditions, the main anions are bicarbonate ( $H C O_{3}^{-}$ ), sulfate ( $S O_{4}^{2 -}$ ), chloride ion (Cl⁻) and carbonate ( $C O_{3}^{2 -}$ ), and the main cations are calcium (Ca²⁺ n), magnesium (Mg²⁺), sodium (Na⁺) and potassium (K⁺) (Liang et al., 2012; Ren et al., 2014; Zhou et al., 2014; Zhang et al., 2018; Gao et al., 2020; Jia et al., 2020). Combined with the theory of GFS, the TDS of groundwater generally increases from the recharge area (RA) to the discharge area (DA), and the hydrochemical types often show horizontal zoning, gradually transforming from $H C O_{3}^{-} - C a^{2 +} \cdot M g^{2 +}$ to $S O_{4}^{2 -} - M g^{2 +} \cdot C a^{2 +}$ , and eventually evolve to Cl⁻−Na⁺ (Hao et al., 2020; Wang et al., 2021a; Zhao et al., 2021). In some areas with special geological environment, the chemical evolution of groundwater deviates from the general law. For example, due to the different geological backgrounds of each region, the phenomenon of element process in groundwater in some areas leads to various endemic diseases, such as Sri Lanka, India, the United States and China (Gao et al., 2014; Kaur et al., 2019; Imbulana et al., 2020).

With the gradual expansion of the impact of human activities on nature, human activities continue to change the surrounding environment of groundwater, leading to obvious man-made interference characteristics in the chemical evolution of groundwater, resulting in obvious changes in the chemical components of groundwater. Such as nitrate ( $N O_{3}^{-}$ ) is added to the main ions in groundwater, the content of heavy metals in groundwater in some areas exceeds the background values, arsenic plumbum pollution is spread worldwide, and the change of chemical compositions of groundwater poses a threat to human health to a certain extent (Gao et al., 2014a; Nan et al., 2008; Sinha and Prasad, 2020; Paolalsiam and Mohan, 2021; Wang et al., 2021b; Xiao et al., 2021; Gao et al., 2022). In addition, with the development of industry and the growth of population, the natural recharge of groundwater is insufficient to meet people’s increasing demand for groundwater, the natural balance of groundwater is broken, and overexploitation leads to a series of problems such as land subsidence, ground fissures, land collapse, seawater intrusion and other problems (Zhang and Zhang, 2019).

Although the research on the chemical evolution of groundwater combined with Tóthian theory has achieved rich results which were widely recognized, the distribution law and chemical evolution law of groundwater have changed significantly under the influence of human activities. In the context of today’s industrial revolution, combined with Tóthian theory and human activities, mining the new laws of groundwater distribution and GCE under the influence of human activities can deepen our understanding of the diversity of geological processes, explain the chemical change process of groundwater more comprehensively and objectively, and better reveal the interaction relationship between groundwater and environment.

This paper fully integrates the relevant literatures, aiming at tracing the development and reform of the theory, focusing on combing the typical literature that can reflect the process of theoretical development, and analyzing the development of GCE theory and the main ways of GCE.

Establishment and development of Tóthian theory

Tóthian theory is the theory of gravity driven GFS firstly proposed by József Tóth, a Canadian hydrogeologist. Once published, it has been widely concerned. After the improvement and development of many hydrogeologists, it plays a leading role in hydrogeology. The theory was initially called “Groundwater multi-level nesting theory” (Tóth, 1963), and later evolved into “Gravity driven groundwater flow system theory” (Tóth, 2009). Finally, it was summarized as Tóthian theory in his article “The evolutionary concepts and practical utilization of the Tóthian Theory of regional groundwater flow”, published in 2016 (Tóth, 2016). This section mainly reviews the theoretical establishment and development of Tóthian theory.

Establishment of Tóthian theory

The GFS theory proposed by Tóthian comes from a large number of previous scientific researches. The main classical groundwater researches before Tóthian theory are shown in Table 1. In 1738, Daniel Bernoulli discovered the “Bernoulli theorem”, which had been widely used in hydraulic research. It was defined as that the total head of groundwater is the sum of its position head, pressure head and velocity head, and groundwater always flows from the position with high total head to the place with low total head (Maranzoni, 2020). In 1856, French hydraulic scientist Darcy found that the seepage velocity is directly proportional to the power of hydraulic gradient based on the concept of head, which is called Darcy’s Law (Darcy, 1857). In 1885, Chamberlin pointed out that the topographic elevation difference provided a driving force for groundwater flow (Chamberlin, 1885). King put forward the concept of gravity driven groundwater in 1899 (Figure 1A) (Cederstrom, 1946; Gao et al., 2014b), and Fourmarier proposed the concept of multistage GFS in 1939 (Figure 1b) (Zhang et al., 2015). Meinzer, Meinzer and Wenzel respectively discovered the three-dimensional properties of groundwater overflow and hydraulic gradient in 1936 and 1940 (Meinzer, 1936; Meinzer and Wenzel, 1940). Hubbert (1940) found that the driving force of groundwater was the negative value of the first derivative of fluid potential through Darcy’s law, and drew an approximate flow model of homogeneous permeable medium between gas-water interface potential source and valley potential sink - inter River block network flow diagram (Figure 1C) (Hubbert, 1940). In 1962, Tóth found the defects of Hubbert flow model in the actual hydrogeological work. Combined with the general relationship between piezometric tube surface and terrain in the same well depth in the same rock unit, he applied Laplace equation to solve the drainage problems of simple basin (Tóth, 1962). After continuous correction and improvement, he finally put forward the theory of “Nested multistage flow system” in 1963 (Figure 1D) (Tóth, 1963). Tóth believed that under the control of specific physical, geographical and geological conditions, the groundwater in the basin presented a multi-level flow mode (composed of local-GFS, intermediate-GFS and regional -GFS). In each GFS, from the supply source to the discharge sink, the flow changes from downward movement to horizontal movement, and then upward movement (Tóth, 1962; Tóth, 1963; Zhang et al., 2018).

TABLE 1

Table 1 Main classical theories on groundwater before Tóthian theory.

FIGURE 1

Figure 1 The evolution diagram of Tóthian Theory (A) The cross-section drawn of the gravity-driven groundwater flow by King in 1899; (B) A pattern of complex flow in unconfined aquifer by Fourmarier in 1939; (C) The groundwater flow net diagram of rivers block by Hubbert in 1940; (D) The use of upper head boundary for linear function and sine function superimposed conditions to conduct simulated groundwater flow field by Toth in 1963.

Development of Tóthian theory

Since 1961, Tóth has been committed to the research of GFS and continuously promoted Tóthian theory. The improvement of GFS and its environmental effects are shown in Figure 2. In 2009, Tóth systematically combed his scientific researches of nearly half a century and published the book “Graphic Systems of Ground Water Flow” (Tóth, 2009). During this period, Tóth’s achievements in gravity driven GFS could be briefly summarized as follows: (1) the main driving force of groundwater is gravity. The gravity penetration theory is proposed to expand the application scope of GFS theory from homogeneous small basins to heterogeneous large basins (Tóth, 1978; Tóth, 1979; Tóth, 1980; Tóth, 1995); (2) Reveal the role of GFS in material transport, hydrochemical evolution and thermal evolution (Tóth, 1966a; Tóth, 1984; Tóth, 1996); (3) It is considered that groundwater is an ubiquitous geological agent and the main reason for geological diversity (Tóth, 1971; Tóth, 1999; Tóth, 2016); (4) It promotes the application of GFS theory in practical work (Tóth, 1966b; Tóth, 1977; Tóth and Corbet, 1986).

FIGURE 2

Figure 2 Natural conditions and phenomena due to environmental agency of flowing groundwater in drainage basin (modified after Tóth, 1999 and Tóth, 2016).

In addition to Tóth himself, other hydrogeological scholars have made continuous exploration based on Tóthian theory, which played an important role in the enrichment, development and promotion of the theory. Among them, the representative figures are Freeze and Witherspoon, Maclay, Jay, Schwartz, Engelen, Allen, Zijl, etc. Freeze and Witherspoon (1966;1967) extended the application scope of Tóthian theory from special cases to general cases through numerical methods, and explored the influence of water level and permeability coefficient on groundwater flow mode (Freeze and Witherspoon, 1966; Freeze and Witherspoon, 1967). Maclay and Winter (1967) combined Tóthian theory with the study of GCE. They found that the relative concentration of major ions in groundwater changed regularly from RA to DA. They believed that the combination of Tóthian theory and hydrochemical research could make people better understand the relationship between water quality and groundwater movement (Maclay and Winter, 1967). Lehr was one of early hydrogeological experts to carry out the indoor observation experiment of GFS. He used plexiglass to build a transparent box filled with sand and resin materials, simulated sandstone to build a hydraulic model, and clearly observed two main phenomena in the movement of groundwater: First, groundwater always moves from high head to low head and the second is the characteristics of layered movement of groundwater (Lehr, 1968). Schwartz and Domenico (1973) put forward the conceptual model of solute transport in GFS in combination with Tóthian theory, and summarized the law of groundwater evolution involved in the processes of groundwater flow (Schwartz and Domenico, 1973). Engelen and Jones (1986) participated in the compilation of the book “Developments in the Analysis of Ground Water Flow Systems” published by the International Association of Hydrological Sciences (Engelen and Jones, 1986). Later, Engelen and Kloosterman (1996) presided over the compilation of the book “Introduction to flow systems”. The publication of these two books not only affirms Tóthian theory, but also extends Tóthian theory to a wider world to a great extent (Engelen and Kloosterman, 1996). In 1993, Allen et al. explored the movement law of groundwater in convergent aquifer, and the main conclusions obtained are basically the same as Tóthian theory, which further confirms the applicability of Tóthian theory (Allen T. Hjelmfelt and Pi, 1993). Zijl (1999) explored the temporal and spatial problems of GFS theory and found that Tóthian theory was also very useful in analyzing the time and space of GFS and their relationship (Zijl, 1999). Jiang et al. (2010) found that due to the attenuation of hydraulic conductivity and porosity depth, the renewal and aging of groundwater in the basin exist at the same time (Jiang et al., 2010). Tóth (2016) combed the development of groundwater flow evolution theory and believed that the current GFS had broken the traditional concept of aquifer, and most of the current GFS were non-confined cross stratum flow systems (Tóth, 2016). Zhang et al. (2020a) found that the change of salinity in groundwater DA would change the position of stagnation point in GFS, thus changing the migration law of groundwater flow (Zhang et al., 2020a). With the maturity of GFS theory, people actively apply Tóthian theory to practical hydrogeological work. For example, people in countries such as Canada, the United States, China and Japan have widely used GFS theory in hydrogeological survey and mapping, water source and sewage discharge site selection, ecosystem protection and so on; and fruitful theoretical and practical results were achieved (Cloutier et al., 2006; Palmer et al., 2006; Han et al., 2009; Tóth, 2009; Lei et al., 2012; Zhang et al., 2015; Ono et al., 2019). By means of isotopes, it was found that there was not only a recharge relationship between surface water and groundwater, but also a GFS at the bottom of the ocean (Krall et al., 2017; Xu et al., 2018b; Danish et al., 2020). This discovery further promoted the application of GFS to a broader world.

In China, Tóthian theory has also been widely concerned and applied. Chen introduced Tóthian theory into China (Chen, 1987). Wang briefly reviewed the development of Tóthian theory in the Book “Fundamentals of Hydrogeology” and fully affirmed the significance of Tóthian theory in the history of hydrogeology (Wang, 1995). Later, the book “Fundamentals of Hydrogeology” compiled by Zhang et al. fully drew lessons from Tóthian theory (Zhang et al., 2018). The book of Zhang et al. also includes the multi-level GFS observed by Liang through the GFS demonstrator (Liang et al., 2010). Based on Tóthian theory, Jiang et al. (2009; 2011; 2012b) explored the relationship between permeability coefficient and depth, the distribution law of stagnation points in GFS, the spatial division characteristics of groundwater age, etc. It was found that the permeability coefficient in regional GFS decreased exponentially with the increasing depth, which would affect the distribution of stagnation points (Jiang et al., 2009; Jiang et al., 2011; Jiang et al., 2012b). Gao et al. (2013; 2014b; 2014c) successfully observed the multistage distribution of GFS, the refraction of groundwater in aquifers with different permeabilities and the driving of heat to groundwater through simple sand trough experiment (Gao, 2013; Gao et al., 2014b; Gao et al., 2014c). In 2015, Zhang et al. translated and published Tóth’s book “Gravitational Systems of Groundwater Flow”, which greatly improved the systematic dissemination of Tóthian theory in China (Zhang et al., 2015). Zhao et al. (2008) and Zhang et al. (2017) respectively investigated the GFS in Ordos Basin and they found that the regional - GFS was consistent with the GFS of Tóthian theory, which provides an excellent example for the research of Tóthian theory and is of great value to the research and practical application of Tóthian theory (Zhao et al., 2008; Zhang et al., 2017). Wang et al. (2017) proposed three-dimensional groundwater circulation units and combined these units into a three-dimensional GFS to obtain the fine structure characteristics of groundwater circulation (Wang et al., 2017).

Tóthian theory breaks through the concept of aquifer and makes people realize that there is no absolute water resisting body (layer) in nature. In the huge structural unit (groundwater aquifer system or hydrogeological unit), driven by the potential energy given by the terrain of the RA and the energy difference between it and the discharge point (place), the continuously moving groundwater flow runs in all directions of releasing energy, with low resistance and smooth flow speed. The gradient of potential energy reduction is large and so is the energy release, and vice versa (Gao et al., 2014b; Tóth, 1970; Tóth, 1978; Winter et al., 1998; Li and Hao, 1999; Tóth, 1999; Tóth, 2009; Liang et al., 2010; Jiang et al., 2011; Wang et al., 2011; Jiang et al., 2012a; Liang et al., 2013). If we combine the concept of GFS with hydrochemical methods, we can better understand the relationship between groundwater quality and movement, and better study the mechanism of GCE.

Chemical evolution of groundwater

Theoretical progress of groundwater chemical evolution

The theoretical study of GCE can be divided into two stages according to the establishment of Tóthian theory. The first stage is the study of GCE independent from Tóthian theory before 1963; The second stage is the study of GCE combined with Tóthian theory after 1963.

Before the emergence of Tóthian theory, the research on the chemical evolution of groundwater was generally limited to the aquifer, but a series of remarkable achievements were still made. For example, as early as the 1860s, geochemists established the basic principle of thermodynamic equilibrium that could be used to distinguish the types of natural water bodies (Glynn and Plummer, 2005). Early scientists believed that groundwater was partially balanced, and pH and redox conditions could reflect the ionic composition, mineral reaction and water stability of water (Korzhinskii, 1936; Krumbein and Garrels, 1952). There are many chemical classification methods for groundwater. In 1934, Schukarev took six main components in water as the basis of classification and proposed the classification method of natural water chemical analysis data, which is the most widely used classification method at present (Zhang and Zhang, 2019). In order to more conveniently represent the chemical types and characteristics of groundwater, Piper three-line diagram (1944) and Durov diagram (1948) appeared one after another. These two hydrochemical diagrams are still widely used in today’s hydrochemical research. Chadha (1999) proposed a rectangular chart that could be directly classified by Excel based on the above two hydrochemical charts (Chadha, 1999). The replacement order of cations is basically determined. Magnesium is (Mg²⁺ ) preferentially replaced in water containing calcium (Ca²⁺ ) and Mg²⁺ The order of relative replaceability of cations in clay is as follows: Li⁺>Na⁺>H⁺>K⁺>Mg²⁺>Ca²⁺ (Kelly, 1934; Ross, 1943). Cederstrom (1946) established a groundwater chemical model with the concept of aquifer (Cederstrom, 1946). It was found that the mineralization of groundwater in different aquifers decreased with the increase of depth. At the same time, combined with the sampling and analysis data, the groundwater near the coastal plain was divided into four zones from west to east: descending zone, hard water zone, soft water zone and high chlorine zone (Cederstrom, 1946). Back et al. (1960) proposed the concept of hydrochemical facies based on aquifer system to describe the main ions in groundwater. It was found that in shallow sediments of the North Atlantic coastal plain, the distribution of $H C O_{3}^{-}$ was wider than $S O_{4}^{2 -}$ and Cl⁻, whereas $N a C l$ was dominant in deeper sediments (Back, 1960; Back and Hanshaw, 1971b). In the 1950s, the study of groundwater isotopes began to show clues. Munnich and Vogel (1959) used the ¹⁴C method to determine the age of groundwater (Munnich and Vogel, 1959). Craig (1961) established the global atmospheric precipitation linear equation δD=8δ¹⁸O+10 to trace the source of groundwater (Craig, 1961). During this period, great progress were made in the theoretical research of GCE and many of them are still used today. However, the research on GCE is still limited, and the exploration of GCE law in different stages is insufficient.

The emergence of Tóthian theory has greatly promoted the development of GCE. For example, Maclay and Winter (1967) found that the evolution sequence of main hydrochemical types of groundwater from RA to runoff area to DA: $H C O_{3}^{-} \to H C O_{3}^{-} + S O_{4}^{2 -} \to S O_{4}^{2 -} \to S O_{4}^{2 -} + C l^{-} \to C l^{-} + S O_{4}^{2 -} \to C l^{-}$ (Maclay and Winter, 1967). Based on Tóthian Theory, Schwartz and Domenico (1973) proposed a conceptual model of solute transport in GFS, which could be defined as: the substance enters the GFS in the RA, then completes the transportation and distribution, and finally discharges in the DA; The chemical and biological effects continuously modify the dissolved substance throughout the process (as shown in Figure 3) (Schwartz and Domenico, 1973). Thorstenson put forward the concept of redox and its application in geochemistry; Then the sequence of redox reactions in GFS began to be studied (Thorstenson, 1984). Wang (1995) clearly distinguished the two concepts of aquifer system and GFS in his research, within which an ideal model diagram was proposed based on the hydrochemical evolution of GFS (as shown in Figure 4) (Wang, 1995). Due to the rapid development of computer technology, hydrogeochemical simulation based on GFS originated in the 1960s developed continuously, which played an important role in revealing the chemical evolution of groundwater. For example, the granite model proposed by Garrels and Mackenzi (1967) laid the foundation for hydrogeochemical simulation (Garrels and Mackenzie, 1967). In 1971, Back and Hanshaw (1971a) pointed out that the application of the basic principle of irreversible thermodynamics in GFS provided a theoretical basis for the construction of GCE and head distribution prediction model (Back and Hanshaw, 1971a). With the improvement of calculation level, PHREEQC began to be widely used in mass balance calculation based on hydrogen and oxygen in 1990 (Wolery et al., 1990). The technology of computer greatly promoted the research of GCE theory, and provided a new research tool for this theoretical research. In particular, the emergence of hydrochemical simulation promoted the hydrochemical research from summarizing the current situations and laws to predicting gradually.

FIGURE 3

Figure 3 Conceptualization of mass transfer in regional groundwater flow (from Schwartz 1973).

FIGURE 4

Figure 4 An ideal model for hydrochemical evolution of groundwater flow system (modified after Wang, 1995).

Main processes of groundwater chemical evolution based on Tóthian theory

In the local GFS which was widely distributed in the shallow part of the surface, the groundwater actively participating in the modern water cycle is discharged to the surface through short-distance, shallow burial and short-time runoff, and its hydrochemical type remains unchanged in bicarbonate type, mostly fresh water with low TDS. However, the groundwater in the regional GFS needs to go through runoff with long distance, deep burial and long time, and the hydrochemical type of groundwater gradually evolves. For example, the transformation from $H C O_{3}^{-}$ type to $S O_{4}^{2 -}$ type and then to Cl⁻ type mainly goes through certain evolution paths (as shown in Figure 5) (Li, 1991; Liang et al., 1991; Tóth, 1999; Chen et al., 2010; Liu et al., 2010; Jiang, 2013).

FIGURE 5

Figure 5 Three paths of groundwater chemical evolution.

The chemical types of groundwater change with time, and vary from place to place. It is generally considered that groundwater experiences such actions as ion alternating adsorption, precipitation and dissolution of salt substances, microorganisms such as desulfurization bacteria and nitrifying bacteria participate in the action, and groundwater undergoes redox reaction, etc. (Gao et al., 2019; Chen et al., 2020a; Gao et al., 2020; Xia et al., 2020). The main evolution form of cations in groundwater are shown in Figure 6. In the RA, atmospheric precipitation enters the GFS through the aeration zone. In this process, leaching first occurs. Ca²⁺, Mg²⁺, Na⁺ in soil or rock are dissolved into groundwater. Due to the low content of CO₂ in the RA and high content of TDS in atmospheric precipitation in the RA, the dissolution of Ca²⁺, Mg²⁺ in this period is greater than Na⁺. In the runoff area, with the migration of groundwater, TDS continues to increase and the content of CO₂ decreases, along with the decrease of the solubility of Ca²⁺ and Mg²⁺, and the solubility of Na⁺ changes little. At the same time, Ca²⁺. and Mg²⁺ in groundwater exchange and adsorb with Na⁺ in surrounding geological bodies, so that the content of Ca²⁺ and Mg²⁺ in groundwater in the runoff area can either be greater or smaller than Na⁺. The TDS of groundwater in the DA continues to increase, the water pressure decreases during the discharge process, Ca²⁺ and Mg²⁺ precipitate, and Na⁺ continues to dissolve. Finally, the Na⁺ content in groundwar is greater than Ca²⁺ and Mg²⁺ (Ren et al., 2014; Wang, 1995; Zhang et al., 2015).

FIGURE 6

Figure 6 Schematic diagram of cation evolution in groundwater flow system.

The evolution form of anions in GFS can be shown in Figure 7. In the local GFS, the groundwater actively participating in the modern water cycle is discharged to the surface through short-distance, shallow burial and short-time runoff, and its hydrochemical type remains unchanged in $H C O_{3}^{-}$ type, characterized by low TDS of most fresh water. However, the groundwater in the regional GFS needs to go through long-distance, deep burial and long-time runoff, and the hydrochemical type of groundwater gradually evolves from $H C O_{3}^{-}$ type to $S O_{4}^{2 -}$ type and then to Cl⁻ type (Wang, 1995; Ren et al., 2014; Zhang et al., 2018). When the RA is dominated by limestone and dolomite, the main chemical reactions are (Zhang et al., 2018):

FIGURE 7

Figure 7 Schematic diagram of anion evolution in groundwater flow system.

\begin{array}{l} {① CaCO}_{3} + H_{2} O + {CO}_{2} \to 2 {HCO}_{3}^{-} + {Ca}^{2 +} & ; \end{array}

\begin{array}{l} {② MgCO}_{3} + H_{2} O + {CO}_{2} \to 2 {HCO}_{3}^{-} + {Mg}^{2 +} & ; \end{array}

When the RA is dominated by albite and calciclase, the main chemical reactions are:

\begin{array}{l} {③ Na}_{2} {Al}_{2} {Si}_{6} O_{6} + 2 {CO}_{2} + 3 H_{2} O \to 2 {HCO}_{3}^{-} + 2 {Na}^{+} + H_{4} {Al}_{2} {Si}_{2} O_{4} + 4 {SiO}_{4} & ; \end{array}

\begin{array}{l} ④ CaO \cdot {Al}_{2} O_{3} \cdot 2 {SiO}_{2} + 2 {CO}_{2} + 3 H_{2} O \to 2 {HCO}_{3}^{-} + {Ca}^{2 +} + H_{4} {Al}_{2} {Si}_{2} O_{9} & ; \end{array}

Main chemical reactions in RA:

\begin{array}{l} {⑤ CaSO}_{4} \cdot 2 H_{2} O \to {Ca}^{2 +} + {SO}_{4}^{2 -} + 2 H_{2} O & ; \end{array}

\begin{array}{l} {⑥ FeS}_{2} + 7 O_{2} + 2 H_{2} O \to 2 {FeSO}_{4} + 4 H^{+} + 2 {SO}_{4}^{2 -} & (oxidation environment) \end{array}

Main chemical reactions in the DA:

\begin{array}{l} {⑦ SO}_{4}^{2 -} + 2 C + 2 H_{2} O \to H_{2} S + 2 {HCO}_{3}^{-} & (reducing environment) \end{array}

New theories and methods of groundwater chemical evolution

In the current situation of increasingly serious water pollution caused by the development of industry and agriculture, the traditional groundwater theories are facing enormous challenges. Therefore, new theories and wider research range related to groundwater arise. Due to the rapid development of computer level, the innovation of new testing methods like isotopic methods, the hydrogeochemical analysis has become more and more convenient and effective.

New theories of groundwater chemical evolution

(1) The reverse cation exchange mechanism has been gradually revealed. Monsoon climate (Kumar et al., 2009), environmental pollution (Reddy et al., 2012; Guma et al., 2021), water acidification (Zaidi et al., 2015), seawater intrusion (Chen et al., 2020b; Roy and Zahid, 2021) and other problems may become the trigger mechanism of the reverse cation exchange process. Through leaching experiment, David and Dimitrios (2002) found that when the soil leaching time reached 15-100h, reverse cation exchange occurred (David and Dimitrios, 2002). Reverse cation exchange can reduce the content of Ca²⁺ and increase the content of Na⁺ in groundwater (Askri, 2015), which can be expressed as: Ca²⁺(Mg²⁺)−water+2Na⁺−clay=2Na⁺−water+Ca²⁺(Mg²⁺)−clay (Chen et al., 2020b).

(2) The content of $N O_{3}^{-}$ in groundwater increased significantly, and the traditional GCE theory encountered new challenges. Since the mid-20th century, nitrogen fertilizer has been widely used in agricultural production, resulting in the continuous increase of $N O_{3}^{-}$ /Cl⁻ content in groundwater in agricultural areas (Abascal et al., 2021) and affecting hydrochemical types (Xu et al., 2018a). In some areas, nitric acid water has appeared (Zhou and Zhu, 2014). The entry of $N O_{3}^{-}$ into groundwater triggered a series of chemical reactions (Mahaqi et al., 2021), which changed the original chemical evolution of groundwater to a great extent (Zaryab et al., 2022; Zhang et al., 2022). For example, the nitrate could be eliminated from nitrate-polluted groundwater by this chemical reaction: 5 CH₂O + 4 ${NO}_{3}^{-}$ = 2N₂ + $4 {HCO}_{3}^{-}$ + CO₂+3H₂O (Mahaqi et al., 2021). There is a strong correlation between $N O_{3}^{-}$ and Cl⁻, and denitrification occurs when the ratio of $N O_{3}^{-}$ /Cl⁻ decreases (Rezaei et al., 2017; Su et al., 2020). $N O_{3}^{-}$ content was negatively correlated with pH and Mn²⁺, and positively correlated with ORP (Oxidation-Reduction Potential) (Cong et al., 2021; Gómez et al., 2002). Nitrogen cycle will affect the migration and transformation of arsenic (As) and iron (Fe) (Gao et al., 2021b; Smith et al., 2017)

(3) The atmosphere is providing more Cl⁻ for groundwater, resulting in significant changes in the hydrochemical type. Cl⁻ in groundwater mainly comes from precipitation and salt dissolution (Zhang et al., 2018). With the growth of runoff path, Cl⁻ in groundwater gradually occupies a dominant position (He et al., 2011). As mentioned above, most of the groundwater in the vadose zone is $H C O_{3}^{-}$ type r and Cl⁻ type in the DA. However, with the intensification of industrial pollution, the concentration of Cl⁻ in the atmosphere increases due to the massive emission or discharge of chlorinated compounds, and the Cl⁻ supplied to the groundwater through rainfall continues to increase, which makes the groundwater in the vadose zone in some areas only present a Cl⁻ type. For example, it’s found that Cl⁻ of the groundwater in vadose zones in some valleys may come from the atmosphere rather than geological bodies (Gardner et al., 2020). In some coastal areas, the initial type of groundwater may also show Cl⁻ type due to the influence of seawater aerosol (Gao et al., 2017).

(4) The research on the migration and transformation mechanism of heavy metals in groundwater has attracted more and more attention. With the rapid development of urbanization and industrialization, there are many forms of heavy metal pollution worldwide. The migration and transformation mechanism of heavy metals in groundwater have attracted extensive attention (Adimalla, 2019). For example, in recent years, scientists have found that the increase of Cl⁻ concentration is conducive to the migration of Mn²⁺ , and seawater intrusion can promote the release of Mn²⁺ from rocks to groundwater and enrich Mn²⁺ in groundwater (Russak et al., 2016; Gao et al., 2021a). Fe and Mn²⁺ ions can form complexes with other anions in groundwater to increase the dissolution of As (Zhang et al., 2020b). Natural colloids are conducive to the diffusion of heavy metal ions (Cai et al., 2016), and low pH conditions are more conducive to the dissolution of heavy metals (Jacques et al., 2008).

New methods of groundwater chemical evolution

(1) More isotopes and isotope models are involved in the study of hydrochemical evolution.

In addition to the traditional ¹⁴C, ³H and ¹⁸O isotopes widely used in groundwater traceability analysis (Craig, 1961; Zhang et al., 2006), isotopes such as ³²Si, ³⁷Ar, ⁸⁵Kr, ⁸⁷Sr, ²²⁶Ra and ²²²Rn are gaining more and more popularity in terms of groundwater source and age analysis (Knies et al., 2000; Bauer et al., 2001; Krall et al., 2017; Danish et al., 2020; Seltzer et al., 2021), which makes the research on GCE more diversified. Isotope mixing model and Bayesian mixing model are also widely used in the identification of pollutants and the calculation of pollutant contribution rate (Zhang et al., 2022).

(2) Groundwater simulation is gradually moving towards artificial intelligence based on big data.

The traditional hydrogeochemical software such as MINTEQA2, WATEQ, Netpath, EQ3/6 and PHREEQC/QE are constantly upgraded (Zhang and Zhang, 2019), and the functions of groundwater flow numerical simulation software such as GMS, MODFLOW and FEFLOW are more and more perfect (Li et al., 2021; Liu et al., 2022a). Neural network learning, deep learning and machine learning based on big data have been widely applied to solve the problem of high-dimensional uncertain GCE (Mo, 2019; El Bilali et al., 2021). For example, Liu et al. (2022) explored the low-temperature hydrogeological process of frozen soil environment through three-dimensional numerical simulation (Liu et al., 2022b). Meanwhile, it is a significant trend that big data models were employed to study the groundwater pollution. Podgorski and Berg created a global prediction map of groundwater arsenic exceeding 10 mg/L using a random forest machine-learning model based on geospatial environmental parameters and more than 50,000 aggregated data points of measured groundwater arsenic concentration (Podgorski and Berg, 2020). Keesari et al. studied big data and environmental sustainability based integrated framework for isotope hydrology applications in India (Keesari et al., 2021). Correa et al. proposed and Embedded Edge-processing IoT-based Water Quality Monitoring system for the sake of proper decision-making (Correa et al., 2022). Wang et al. applied more than 4600 water samples for the artificial neural network to predict genesis, occurrence and distribution of high iodine groundwater in China (Wang et al., 2022). Using the similar method to Wang et al. (2022), Cao et al. predicted geogenic groundwater fluoride contamination throughout China (Cao et al., 2022). Xiong et al. (2022) applied machine learning to optimize the groundwater pollution monitoring network (Xiong et al., 2022). It is also noted that geological big data sharing service has received great attention from the governments. For example, the British Geological Survey, supported by the British government, is the first country to build the geological database and they are at the forefront of sharing data worldwide through a comprehensive system in this decade (Normile, 2019; Wen et al., 2021). In China, the Chinese Geological Survey also shows the sharing attitude by developing and opening a geological database called “GeoCloud”, within which a wealth of groundwater data was accessible (Lei et al., 2021). By sharing the groundwater database between different regions, hydrologists can use different types of databases from similar geological background for reference and achieve comparative results of GFS and GCE (Normile, 2019).

In short summary, implementation of big-data methodology in terms of isotope hydrology, water pollution and environmental sustainability has drawn more and more attention. The connected and sharing geological database throughout China or other regions in the world makes it more convenient to carry out big-data based groundwater research.

Conclusion

From Darcy’s Law (1856) to GFS theory (1963), then to Tóthian theory (2016), scientists have been exploring the flow law of water for more than 160 years and will continue to do so. In the more than 50 years since the establishment of GFS force theory to finally become Tóthian theory, its connotation has been gradually enriched through the continuous development, enrichment and practice of many hydrogeologists. Tóthian theory mainly covers the distribution and driving mechanism of GFS, the transmission and distribution of matter and energy by GFS, the interaction between groundwater and environment, and the environment associated phenomena brought by groundwater. Tóthian theory is a milestone in the history of hydrogeology, which breaks the concept of aquifer and promotes the development of hydrogeology from small-scale homogeneous basin to large-scale heterogeneous basin. As a universal geological agent, the environmental effect of groundwater is gradually accepted by people. With the progresses of science and technology, the research scopes of GFS have been extended from the original saturated zone to the non-saturated zone, realizing the leap from two-dimensional space to three-dimensional space, and describing the structural characteristics of GFS more finely.

The research on the mechanism of GCE based on GFS theory and Tóthian theory have promoted each other and grown together in the past 50 years. People have mastered some general laws of GCE, such as: (1) the replacement order of cations in the process of cation exchange adsorption, (2) the main water rock reactions and the main hydrochemical types in the RA, runoff area and DA under different geological backgrounds, (3) horizontal zoning of GCE. (4) in GFS, groundwater TDS generally increases from RA to DA, (5) sources of groundwater in GFSs in different regions, (6) location and sequence of redox reaction. With the expansion of the research scopes of GFS to the unsaturated zone, soil microorganisms and CO₂ have become a research hotspot in the study of GCE. The maturity of traditional isotopic methods and the addition of new isotopes have greatly improved the reliability of the research conclusion of hydrochemical evolution. The maturity of hydrochemical evolution software and the application of neural network learning and deep learning have solved the high-dimensional uncertainty problem.

With the deepening of theoretical research, more and more scientists are aware of the interaction between the distribution of GFS and GCE. The theoretical problems of GFS and the research of GCE are interdependent and inseparable. In order to better deal with the problems of GFS distribution change, GCE environment change and groundwater pollution caused by the further intensification of human transformation of nature in the 21st century, it is also a difficult problem for hydrogeologists to comprehensively use the scientific research results of Tóthian theory and GCE mechanism to promote the harmonious coexistence between man and nature.

Author contributions

HD, original draft preparation, software, and visualization. ZG, review. All authors contributed to the article and approved the submitted version.

Funding

This review was supported by the Natural Science Foundation of Shandong Province (ZR2020MD109); National Natural Science Foundation of China (40772145, 41641022).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Abascal E., Gómez-Coma L., Ortiz I., Ortiz A. (2021). Global diagnosis of nitrate pollution in groundwater and review of removal technologies. Sci. Total Environ. 810, 152233. doi: 10.1016/j.scitotenv.2021.152233

PubMed Abstract | CrossRef Full Text | Google Scholar

Adimalla N. (2019). Heavy metals contamination in urban surface soils of medak province, India, and its risk assessment and spatial distribution. Environ. Geochem. Health42 (1), 59–75. doi: 10.1007/s10653-019-00270-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Allen T. Hjelmfelt J., Pi Z. (1993). Movement of ground water in converging aquifer. J. Irrig. Drain. Eng.119 (2), 312–322. doi: 10.1061/(ASCE)0733-9437(1993)119:2(312)