Hydrochemical Characteristics and Evolution Mechanisms of Shallow Groundwater in the Alluvial–Coastal Transition Zone of the Tangshan Plain, China

To elucidate the hydrochemical characteristics and evolution mechanisms of shallow groundwater in the alluvial–coastal transitional zone of the Tangshan Plain, 76 groundwater samples were collected in July 2022. An integrated approach combining Piper and Gibbs diagrams, ionic ratio analysis, multivariate statistical methods (including Pearson correlation, hierarchical cluster analysis, and principal component analysis), and PHREEQC inverse modeling was employed to identify hydrochemical facies, dominant controlling factors, and geochemical reaction pathways. Results show that groundwater in the upstream alluvial plain is predominantly of the HCO3–Ca type with low mineralization, primarily controlled by carbonate weathering, water–rock interaction, and natural recharge. In contrast, groundwater in the downstream coastal plain is characterized by high-mineralized Cl–Na type water, mainly influenced by seawater intrusion, evaporation concentration, and dissolution of evaporite minerals. The spatial distribution of groundwater follows a pattern of “freshwater in the north and inland, saline water in the south and coastal,” reflecting the transitional nature from freshwater to saline water. Ionic ratio analysis reveals a concurrent increase in Na+, Cl−, and SO42− in the coastal zone, indicating coupled processes of saline water mixing and cation exchange. Statistical analysis identifies mineralization processes, carbonate weathering, redox conditions, and anthropogenic inputs as the main controlling factors. PHREEQC simulations demonstrate that groundwater in the alluvial zone evolves along the flow path through CO2 degassing, dolomite precipitation, and sulfate mineral dissolution, whereas in the coastal zone, continuous dissolution of halite and gypsum leads to the formation of high-mineralized Na–Cl water. This study establishes a geochemical evolution framework from recharge to discharge zones in a typical alluvial–coastal transitional setting, providing theoretical guidance for salinization boundary identification and groundwater management.