Hydrochemical Formation Mechanisms and Source Apportionment in Multi-Aquifer Systems of Coastal Cities: A Case Study of Qingdao City, China

This study systematically unravels the hydrochemical evolution mechanisms and driving forces in multi-aquifer systems of Qingdao, a coastal economic hub. Integrated hydrochemical analysis of porous, fissured, and karst water, combined with PHREEQC modeling and Positive Matrix Factorization (PMF), deciphers water–rock interactions and anthropogenic perturbations. Groundwater exhibits weak alkalinity (pH 7.2–8.4), with porous aquifers showing markedly higher TDS (161.1–8203.5 mg/L) than fissured (147.7–1224.8 mg/L) and karst systems (361.1–4551.5 mg/L). Spatial heterogeneity reveals progressive hydrochemical transitions (HCO3-Ca → SO4-Ca·Mg → Cl-Na) in porous aquifers across the Dagu River Basin. While carbonate (calcite) and silicate weathering govern natural hydrochemistry, evaporite dissolution and seawater intrusion drive severe groundwater salinization in the western Pingdu City and the Dagu River Estuary (localized TDS up to 8203.5 mg/L). PMF source apportionment identifies acid deposition-enhanced dissolution of carbonate/silicate minerals, with nitrate contamination predominantly sourced from agricultural runoff and domestic sewage. Landfill leachate exerts pronounced impacts in Laixi and adjacent regions. This study offering actionable strategies for salinity mitigation and contaminant source regulation, thereby providing a scientific framework for sustainable groundwater management in rapidly urbanizing coastal zones.