This paper describes a chemical model that calculates (solid + liquid) equilibria in the { m 1FeCl 2 + m 2FeCl 3}(aq), { m 1FeSO 4 + m 2Fe 2(SO 4) 3}(aq), { m 1NaCl + m 2FeCl 3}(aq), { m 1Na 2SO 4 + m 2FeSO 4}(aq), { m 1NaCl + m 2FeCl 2}(aq), { m 1KCl + m 2FeCl 3}(aq), { m 1K 2SO 4 + m 2Fe 2(SO 4) 3}(aq), { m 1KCl + m 2FeCl 2}(aq), { m 1K 2SO 4 + m 2FeSO 4}(aq), and { m 1MgCl 2 + m 2FeCl 2}(aq) systems, where m denotes molality at T=298.15 K. The Pitzer ion-interaction model has been used for thermodynamic analysis of the experimental activity data in binary FeCl 2(aq) and FeCl 3(aq) solutions, and ternary solubility data, presented in the literature. The thermodynamic functions needed (binary and ternary parameters of ionic interaction, thermodynamic solubility products) have been calculated and the theoretical solubility isotherms have been plotted. The mixed solution model parameters { θ(MN) and ψ(MNX)} have been chosen on the basis of the compositions of saturated ternary solutions and data on the pure water solubility of the K 2SO 4 · FeSO 4 · 6H 2O double salt. The standard chemical potentials of four ferrous {FeCl 2 · 4H 2O, Na 2SO 4 · FeSO 4 · 4H 2O, K 2SO 4 · FeSO 4 · 6H 2O, and MgCl 2 · FeCl 2 · 8H 2O} and three ferric {FeCl 3 · 6H 2O, 2KCl · FeCl 3 · H 2O, and 2K 2SO 4 · Fe 2(SO 4) 3 · 14H 2O} solid phases have been determined. Comparison of solubility predictions with experimental data not used in model parameterization is given. The component activities of the saturated { m 1MgSO 4 + m 2FeSO 4}(aq) and in the mixed crystalline phase were determined and the change of the molar Gibbs free energy of mixing Δ mix G ∘ m( s) of crystals was determined as a function of the solid phase composition. It is established that at T=298.15 K the mixed (Mg,Fe)SO 4 · 7H 2O and (Fe,Mg)SO 4 · 7H 2O crystals show small positive deviations from the ideal mixed crystals. Limitations of the {Fe(II) + Fe(III)} model due to data insufficiencies are discussed.