Tracing regional flow paths to major springs in Trans-Pecos Texas using geochemical data and geochemical models

San Solomon, Giffin, and Phantom Lake Springs, located in Trans-Pecos Texas, have a high TDS, Na–Cl–SO4 baseflow component derived from a regional flow system and a low TDS, mixed cation–mixed anion stormflow component derived from local precipitation events. The hypothesis that the regional flow system maintaining baseflow spring discharge originates in the Salt Basin and flows through the Apache Mountains towards the springs is tested with historical geochemical data from wells and springs. Data from over 1400 wells in the study area over a 50-year period were analyzed and used to delineate 11 hydrochemical facies based on the predominant ions. mical data from samples along the hypothesized regional flow path indicate a trend of increasing dissolved solids and Cl–HCO3 ratios and decreasing Na–Cl ratios. These are consistent with evolution of groundwater in an unconfined regional system dominated by carbonates and evaporites. In the bicarbonate facies, the waters represent recent recharge modified by mineral dissolution and cation exchange. In the sulfate zones, the hydrochemical facies are controlled by gypsum, anhydrite, and halite dissolution, cation exchange, and mixing with Na–Cl waters. In the chloride zones, the hydrochemical facies are controlled by halite dissolution and irrigation return flow. Spring discharge chemistry is most similar to chloride zone waters; Na–Cl and Ca–SO4 ratios suggest that baseflow is derived from the chloride zone waters upgradient along the hypothesized flow path. C modeled groundwater evolution along the hypothesized flow path and spring discharge under stormflow and baseflow conditions. Results indicate that: (1) hydrochemistry along the regional flow path is controlled by dissolution of halite, gypsum, dolomite, and CO2 and by precipitation of calcite; (2) baseflow spring discharge is derived primarily from this regional flow system; and (3) spring discharge after major storm events can constitute as much as 72% local recharge that is further modified by dissolution of calcite, gypsum, and CO2. Data analysis and model results suggest that cave formation in this system is occurs during major storm events.