A comprehensive thermodynamic model based on the electrolyte NRTL (eNRTL) activity coefficient equation is developed for the NaCl + H2O binary, the Na2SO4 + H2O binary and the NaCl + Na2SO4 + H2O ternary. The NRTL binary parameters for pairs H2O-(Na+, Cl-) and H2O-(Na+, SO42-), and the aqueous phase infinite dilution heat capacity parameters for ions Cl- and SO42- are regressed from fitting experimental data on mean ionic activity coefficient, heat capacity, liquid enthalpy and dissolution enthalpy for the NaCl + H2O binary and the Na2SO4 + H2O binary with electrolyte concentrations up to saturation and temperature up to 473.15K. The Gibbs energy of formation, enthalpy of formation and heat capacity parameters for solids NaCl(s), NaCl center dot 2H(2)O((s)), Na2SO4(s) and Na2SO4 center dot 10H(2)O((s)) are obtained by fitting experimental data on solubilities of NaCl and Na2SO4 in water. The NRTL binary parameters for the (Na+, Cl-)-(Na+, SO42-) pair are regressed from fitting experimental data on dissolution enthalpies and solubilities for the NaCl + Na2SO4 + H2O ternary. (C) 2011 Elsevier B.V. All rights reserved.