Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl2 and CaCl2) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K1, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl2 solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.