The objective of this work is to validate and compare PHREEQC and Reaktoro and some of the available databases distributed freely with the software. We do this by systematically comparing the models implemented in both software programs against gas solubility measurements reported in the literature. The models are later used to assess the potential and extent of geothermal fluids degassing at different pressure and temperature conditions. Being able to predict the bubble point pressure of geothermal brine is pivotal for either selecting a sufficiently high operating and design pressure of the surface facilities or appropriately designing the surface facilities to handle any gases formed during operation. For the latter, quantifying the amount of gas formed is essential as well as knowing the composition of the formed gas phase.