Saline aqueous solutions at elevated pressures and temperatures play an important role in processes such as supercritical water oxidation (SCWO) and supercritical water gasification (SCWG), as well as in natural geochemical processes in Earth and planetary interiors. Some solutions exhibit a negative temperature coefficient of solubility at high temperatures, thereby leading to salt precipitation with increasing temperature. Using modulated FT-IR Raman spectroscopy and classical molecular dynamics simulations (MD), we studied the solute speciation in solutions of 10 wt % Na2SO4, at conditions close to the saturation limit. Our experiments reveal that ion pairing and cluster formation are favored as solid saturation is approached, and ionic clusters form prior to the precipitation of solid sulfate. The proportion of such clusters increases as the phase boundary is approached either by decreasing pressure or by increasing temperature in the vicinity of the three-phase (vapor-liquid-solid) curve.