While extant research has predominantly focused on elucidating the mechanisms of fluorine (F⁻) enrichment in groundwater within the North China Plain, the occurrence and formation mechanisms of high-F⁻ groundwater in Xiong’an New Area remain unexplored. Consequently, 365 groundwater samples (172 from shallow groundwater, 193 from deep groundwater) were collected from Xiong’an New Area. Hydrochemical analysis, geochemical modeling, and statistical analysis were used to explore the occurrence and formation mechanisms of high-F⁻ groundwater. The results reveal that the highest F⁻ concentrations in shallow and deep groundwater were up to 3.22 mg/L and 1.79 mg/L, respectively. High-F⁻ groundwater was primarily located at the eastern part of the study area. The distribution area of high-F⁻ shallow groundwater was much greater than that of deep groundwater. F-bearing minerals dissolution and ion exchange were the principal formation mechanisms of high-F⁻ groundwater in both shallow and deep aquifers. Moreover, competitive adsorption, evaporation, and the impacts of Ca²⁺ and Mg²⁺ dissolution equilibrium on F-bearing dissolution were crucial to the formation of high-F⁻ groundwater in shallow aquifers. Desorption in an alkaline environment, evaporites dissolution and salt effects were vital to the formation of high-F⁻ groundwater in deep aquifers. These findings can contribute to the support of local groundwater security and management.

In previous investigations, the demarcation of capture zones within a specific research area predominantly relied on a singular method, leading to pronounced limitations and uncertainties. To address this challenge, an extensive field survey was conducted, focusing on the systematic classification of water sources in the Linfen City region. Building upon this classification, an intricate fusion of a hydrogeological analysis and formulaic methodology was employed. This integrated approach, coupled with independent numerical simulation methods, was applied to delineate recharge areas for both alluvial fan pore water in piedmont regions and exposed karst water in small- to medium-sized water sources. Simultaneously, precise spatial interpolation was carried out on water quality monitoring data from supply wells within the water source area for the year 2020. This meticulous analysis facilitated the determination of the spatial distribution characteristics of hydrochemical elements. To assess the water quality within the capture zone, the class III groundwater quality standards of China were employed as a pivotal tool for validating the results of the delineation of water source recharge areas. In the final analysis, a comparative study between the integrated coupling method and numerical simulation outcomes revealed the successful delineation of the boundaries for the water supply areas of Tumen and Caojiapo in Linfen City, covering areas of 5.5 km² and 22.29 km², respectively. Simultaneously, the combination of the three methods accurately outlined the boundary of the Hexi water supply area, encompassing an area of 2.5224 km². These results vividly illustrate that the amalgamation of various methodologies proves more beneficial for the precise delineation of capture zones, particularly in diverse types and scales of groundwater sources. The synergy exhibited by these three methods underscores their collective efficacy, providing a more comprehensive and intuitive delineation of the recharge areas for small- to medium-sized water sources. Consequently, these findings significantly enhance the practical application value of the study and hold promise in making substantial contributions to local groundwater security and management initiatives.

With the increasing environmental impacts of human activities, the problem of polygenic multipollutants in groundwater has attracted the attention of researchers. Identifying the hydrobiogeochemical characteristics of the surface sewage that replenishes groundwater is crucial to addressing this problem. The input of polygenic multipollutants into groundwater leads to not only the mechanical superposition of pollutants but also the formation of secondary pollutant types. The evolution of polygenic multipollutants is influenced by aquifer characteristics, carbon sources, microbial abundance, etc. Therefore, this study took a sewage leakage point in Northwest China as the research object, carried out a controlled laboratory experiment on the impact of sewage discharge on groundwater, and, combined with long-term field monitoring results, determined the main hydrobiogeochemical processes of polygenic multipollutants and their secondary pollutants. The results showed that the redox environment and the gradient change in pH were identified as the most critical controlling factors. In oxidative groundwater during the early stage of vertical infiltration, sewage carries a substantial amount of NH₄⁺, which is oxidized to form the secondary pollutant NO₃⁻. As O₂ is consumed, the reduction intensifies, and secondary pollutants NO₃⁻, Mn (IV), and Fe(III) minerals are successively reduced. Compared with the natural conditions of rainwater vertical infiltration, the reaction rates and intensities of various reactions significantly increase during sewage vertical infiltration. However, there is a notable difference in the groundwater pH between sewage and rainwater vertical infiltration. In O₂ and secondary pollutant NO₃⁻ reduction, a large amount of CO₂ is rapidly generated. Excessive CO₂ dissolves to produce a substantial amount of H⁺, promoting the acidic dissolution of Mn (II) minerals and generation of Mn²⁺. Sewage provides a higher carbon load, enhancing Mn (II) acidic dissolution and stimulating the activity of dissimilatory nitrate reduction to ammonium, which exhibits a higher contribution to NO₃⁻ reduction. This results in a portion of NO₃⁻ converted from NH₄⁺ being reduced back to NH₄⁺ and retained in the groundwater, reducing the denitrification’s capacity to remove secondary NO₃⁻. This has important implications for pollution management and groundwater remediation, particularly monitored natural attenuation.

Leakage is important in groundwater cycle research and water resource management, the key is to construct the spatial distribution of aquitard. Influenced by natural and anthropogenic factors, the density of boreholes restricts the accuracy of aquitard structure. This paper explores the K-Nearest Neighbor (KNN) method to construct multi-layer aquitard structure in the hinterland of Songnen Plain, Northeast China, with limited boreholes, and compare it with Inverse Distance Weight (IDW) and Ordinary Kriging (OK) methods to evaluate the leakage. The KNN needs to first identify the eigenvalue-k of each layer, which is used to enhance the trainset to obtain the testset with a larger amount of sample data enhanced, and to make a prediction to obtain the 3-dimensional structure of aquitard. KNN is more accurate than the IDW and OK. The aquitrad structure obtained by the KNN is different from the IDW and OK, the thickness difference is 3.53%-54.00%, the area difference is 2.28%-23.91%, and the difference in leakage can be up to 24.51%. KNN also has a high degree of identification of aquitard edges, which can help the prevention of groundwater pollution. The results provide guidance for water balance analysis, borehole engineering design, water resources management and exploitation.