Hydrogeochemical evolution and processes based on multivariable statistical and inverse simulation modeling: A coal mine in the Northern Coalfield, China

Considering groundwater from the aquifers overlying the bedrock is an important water source for drinking purposes. As such, the investigation of its property is essential. Based on the spatial structure of aquifers in the study area, the aquifers in the Cenozoic strata are divided into three groups. The multivariate statistical approaches are employed to identify the hydrogeochemical processes and hydro-chemical types, and hydrogeochemical inverse modeling is applied to further validate and elucidate the hydrogeochemical process and water–rock interaction. The results are as follows: (1) The hydro-chemical type of the upper aquifer is dominated by the K + Na-HCO3 type, while others have similar water quality types, which are dominated by the K + Na-Cl type and the K + Na-SO4 type. (2) The saturation index of anhydrite, gypsum, halite, and CO2(g) is below zero in three aquifers, indicating that they are unsaturated. While aragonite, calcite, and dolomite in the middle aquifer remain in the unsaturated–saturated state. (3) The cation exchange process accelerating the reduction of Ca2+ concentration and the increase of SO42− concentration occurs in three aquifers, and the dissolution of calcite and dolomite minerals occurs in most cases. This study supports the fundamental evidence for the hydrogeochemical processes and water resource utilization and has a certain practical significance. HIGHLIGHTS Hydrogeochemical characteristics of the Cenozoic aquifers are investigated.; Relationship between ions reveals the hydrogeochemical processes and water–rock interaction.; Cation exchange process accelerates the reduction of Ca2+ concentration and the increase of SO42+ concentration.; A hydro-chemical boundary is identified.; Reverse simulation and cluster analysis further validate the water–rock interaction processes.;