(a) Steady-state temperature distribution and fluid flow field. (b) Calculated temperatures fit well with the observed temperatures in ZR1. (c) Gas saturation (Sg) distribution illustrates the existence of unsaturated zone at the top of the model. (d) The calculated chemical components in discharge area fit well with observations (ZR1, G1, G2, G3, and G4). The value of “HCO3-” is the sum of HCO3- and CO32-.