Thermodynamic models are essential to facilitate the advancement of process design, optimization, and operation of electrolyte systems. In this work, a comprehensive thermodynamic framework based on the Electrolyte Nonrandom Two-Liquid model is developed to calculate phase equilibria behavior and salt solubility of aqueous LiCl, LiBr, LiI, and LiNO3 solutions. The model describes the non-ideality of the electrolyte solutions by using two binary interaction parameters for each molecule-electrolyte pair in the system. Each binary interaction parameter is further expressed with up to three temperature coefficients which are regressed from experimental data. To take into account the hydration of lithium ion, two separate chemistries for the dissociation of lithium salts are investigated. In the first case, the lithium ion is considered as a bare ion, Li+, while in the second case hydration of the lithium ion from Li+ to Li(H2O)+ is considered. The calculated thermodynamic properties compare adequately with the experimental data for both sets of chemistries for concentrations up to saturation and temperatures from 273.15 K up to 623.15 K. Moderate to significant improvements are observed with the incorporation of the hydration chemistry for aqueous LiCl, LiBr, and LiI solutions when compared to the non-hydrated lithium ion model results.