Over the last 40 years, declining spring water flow rates have become a typical feature of karst springs in Northern China. Wavelet analysis, the Mann-Kendall trend test and the mutation test were used to analyze dynamic monitoring data of groundwater levels and atmospheric precipitation in the Jinan karst spring area, from 1956 to 2013, to study hydrological responses to atmospheric precipitation over one-year periods. Results from this analysis show that: (1) Atmospheric precipitation and the spring water level displayed multi-scale change characteristics, having two very similar cycles of change of 16 and 12 years. This finding shows that atmospheric precipitation generates a direct impact on the level of spring water. (2) From 1956 to 2013, the groundwater level in the Jinan spring area had a significant downward trend (0.65 m/10a). Precipitation recorded an increasing trend (12.65 mm/10a), however this was not significant. The weight of the influencing factors of the spring dynamic therefore changed due to the influence of human factors. (3) A mutation of atmospheric precipitation occurred in 1999, after which annual precipitation increased. Results for the mutation of the groundwater level showed an initial change in 1967. After this change the water level continued to decrease before rapidly increasing after 2004. The future trend of the spring water level should be maintained with consistent precipitation (having an upward trend), indicating that atmospheric precipitation is not the only factor affecting the dynamics of the spring. (4) Different periods were identified on the multiple regression model. The main influencing factors on groundwater level over the past 58 years were identified as a transition from precipitation to artificial mining. These results also validate the suitability and reliability of using wavelet analysis and the Mann-Kendall test method to study groundwater dynamics; these results provide a reference for the future protection of the Jinan City spring.

Groundwater chemistry and its potential health risks are as important as water availability in arid and semiarid regions. This study was conducted to determine the contamination and associated health threats to various populations in a semiarid basin of north China. A total of 78 groundwater samples were collected from the shallow unconfined aquifers. The results showed that the phreatic water was slightly alkaline, hard fresh water with ions in the order of Ca²⁺ > Na⁺+K⁺ > Mg²⁺ and HCO₃⁻ > SO₄²⁻ > Cl⁻. Four hydrochemical elements, NO₃⁻, F⁻, Mn and Zn, exceeded the permissible limits. NO₃⁻ and F⁻ contaminants may pose health risks to local residents, while the risks of Mn and Zn are negligible. Dermal exposure is safe for all populations, while the oral pathway is not. Minors (i.e., infants and children) are susceptible to both NO₃⁻ and F⁻ contaminants, and adults only to NO₃⁻. The susceptibility of various populations is in the order of infants > children > adult males > adult females. Anthropogenic activities are responsible for the elevated levels of NO₃⁻, Zn, Total dissolved solids (TDS), while F⁻ and Mn are from geogenic sources. Thus, differential water supplies, strict control of waste, and rational irrigation practices are encouraged in the basin.

Prevention and remediation strategies for groundwater pollution can be successfully carried out if the location, concentration, and release history of contaminants can be accurately identified. This, however, presents a challenge due to complex groundwater systems. To address this issue, a simulation-optimization (S/O) model by integrating MODFLOW and MT3DMS into a shuffled complex evolution (SCE-UA) optimization algorithm was proposed; this coupled model can identify the unknown groundwater pollution source characteristics. Moreover, the Grids Traversal algorithm was used for automatically searching all possible combinations of pollution source location. The performance of the proposed S/O model was tested by three hypothetical scenarios with various combinations of mixed situations (i.e., single and multiple pollution source locations, known and unknown pollution source locations, steady-state flow and transient flow). The field measurement errors was additionally considered and analyzed. Our results showed that this proposed S/O model performed reasonably well. The identified locations and concentrations of contaminants fairly matched with the imposed inputs with average normalized deviations less than 1% after sufficient generations. We further assessed the impact of generation number on the performance of the S/O model. The performance could be significantly improved by increasing generation number, which yet resulted in a heavy computational burden. Furthermore, the proposed S/O model performed more efficiently and robustly than the traditionally used artificial neural network (ANN)-based model. This is due to the internal linkage of numerical simulation in the S/O model that promotes the data exchange from external files to programming variables. This new model allows for solving the source-identification problems considering complex conditions, and thus for providing a platform for groundwater pollution prevention and management.

Pore water in inland plain areas, generally having poor water quality, contain complex hydrochemical properties. In order to examine groundwater chemical composition formation characteristics, groundwater in the Jiyang area of Lubei Plain was studied using stratified monitoring of drilling, analysis of water level and water quality, isotope analysis, ion ratio coefficient and isothermal adsorption experiments, hydrochemical characteristics, and analysis of variations in different shallow depths. Results show that: (1) Numerous hydrochemistry types are present in the diving. Along with the direction of groundwater flow, total dissolved solids (TDS) of diving in the study area generally increases and the hydrochemical type changes from the HCO₃ type to the HCO₃·SO₄ type, Cl·HCO₃ type and the Cl·SO₄ type. (2) Shallow brackish water and freshwater in the horizontal direction are alternately distributed, and shallow brackish water is distributed in the area between old channels, showing sporadic spots or bands, whose hydrochemistry type is predominantly Cl·SO₄-Na·Mg·Ca. (3) Affected by the sedimentary environment, hydrodynamic conditions and other factors; diving, middle brackish water and deep freshwater are vertically deposited in the study area. The dynamics of middle brackish water quality are stable due to the sedimentary environment and clay deposits. The hydrochemistry types of middle brackish water are mainly Cl·SO₄-Mg·Na and SO₄·Cl-Na·Mg, while the deep confined water is dominated by HCO₃. (4) The optimal adsorption isotherms of Na⁺, Ca²⁺ and Mg²⁺ in groundwater from clay, with a thickness raging from 6–112 m, conformed to the Henry equation and the Langmuir equation. The retardation of Na⁺, Ca²⁺ and Mg²⁺ in groundwater differed with differing depths of the clay deposit. The trend of change in retardation strength correlates strongly with the TDS of groundwater. Groundwater in the inland plain area is affected by complicated hydrogeochemical processes; in addition, leaching, cation exchange, evaporative concentration and clay retardation have an important influence on the horizontal and vertical zonation of groundwater chemical components.

In order to increase understanding of the hydrogeochemical effects that influence changes in the quality of salt water, we investigated the distribution of saline and fresh water in an inland plain area and, in particular, the scarcity of fresh water resources. Taking the inland plain in Jiyang County as a specific case study, samples of undisturbed clay and underground saline water from different depths were collected to examine hydrogeological changes. A wide variety of methods was used to analyze the blocking effect of clay on the chemical characteristics of the groundwater. These include real-time monitoring for field water quality, tests for isothermal adsorption, a factor analysis model, physiochemical analysis, and correlation analysis. Our results show that the optimal adsorption isotherm of clay for Na⁺, Ca²⁺ and Mg²⁺ in groundwater conform to the established Henry and Langmuir equations for adsorption isotherms. The influence of clay mineral types and content on the blockage of Na⁺, Ca²⁺ and Mg²⁺ in groundwater samples were evident at different depths, with the clay adsorption capacity increasing in line with increases in the clay mineral content. Clay at different depths was found to have the strongest blocking effect on Na⁺ in groundwater, being systematically greater than its effect on Ca²⁺ and Mg²⁺. It is believed that the blocking effect of clay has an important influence on the hydrochemical zoning of groundwater in inland plains and the formation of saline water in groundwater systems. This study therefore provides concrete evidence in support of this supposed effect.

Solute transport modeling resolves advection, dispersion, and chemical reactions in groundwater systems with its accuracy depending on the resolution of domain at all scales, thus the computational efficiency of a simulator becomes a bottleneck for the wide application of numerical simulations. However, the traditional serial numerical simulators have reached their limits for the prohibitive computational time and memory requirement in solving large-scale problems. These limitations have greatly hindered the wide application of groundwater solute transport modeling. Thus, the development of an efficient method for handling large-scale groundwater solute transport simulation is urgently required. In this study, we developed and assessed a parallelized MT3DMS (Modular Three-Dimensional Multispecies Transport Model) by using OpenMP (Open specifications for Multi-Processing) to accelerate the solute transport simulation process. The parallelization was achieved by adding OpenMP compile directives (i.e., defining various types of parallel regions) into the most time-consuming packages, including the Advection package (ADV), Dispersion package (DSP), and Generalized Conjugate Gradient Solver package (GCG). This allows parallel processing on shared-memory multiprocessors, i.e., both the memory requirement and computing efforts are automatically distributed among all processors. Moreover, we discussed two different parallelization strategies for handling numerical models with either many layers or few layers. The performance of parallelized MT3DMS was assessed by two benchmark numerical models with different model domain sizes via a workstation with two quad-core processors. Results showed that the running time of parallelized MT3DMS can be 4.15 times faster than that using sequential MT3DMS. The effects of using different preconditioners (procedures that transform a given problem into a form that is more suitable for numerical solving methods) in the GCG package were additionally evaluated. The modified strategy for handling numerical models with few layers also achieved satisfactory results with running time two times faster than that via sequential simulation. Thus, the proposed parallelization allows high-resolution groundwater transport simulation with higher efficiency for large-scale or multimillion-cell simulation problems.