Controls on the distribution of hydrous defects in forsterite from a thermodynamic model