This paper describes a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the Li-H-Na-K-Cl-OH-H 2 O system from dilute to high solution concentration within the 0 to 250°C temperature range. The model coherently extends to Li the temperature-variable H-Na-K-OH-Cl-H 2 O model of Christov and Møller (2004b, p. 1309). The solubility modeling approach based on Pitzer's (1973, p. 268) specific interaction equations is used. All binary (LiCl-H 2 O and LiOH-H 2 O) and ternary (LiCl-HCl-H 2 O, LiCl-NaCl-H 2 O, LiCl-KCl-H 2 O, LiOH-NaOH-H 2 O, LiOH-KOH-H 2 O, and LiOH-LiCl-H 2 O) lithium subsystems are included in the model parameterization. The model for the LiCl-H 2 O system is parameterized using two different approaches: (1) with 4 ion interaction binary parameters ( (0) ,  (1) ,  (2) , and C ), and (2) with 3 ion interaction binary parameters ( (0) ,  (1) , and C ) and including neutral aqueous LiCl 0 (aq) species. Approach (2) provides a better fit of activity data in unsaturated binary solutions and accurately predicts solid solubilities up to 40 mol.kg ؊1 and up to 250°C. Therefore, this approach was used to parameterize lithium chloride mixed systems. Temperature functions for the thermodynamic solubility product (as log K o sp ) of 5 simple lithium salts (LiCl.2H 2 O(cr), LiCl. H 2 O(cr), LiCl(cr), LiOH. H 2 O(cr), and LiOH(cr)) are determined. The log K o sp values of 3 double lithium basic salts precipitating in the LiOH-LiCl-H 2 O system at 50°C (3LiOH. LiCl(cr), LiOH.LiCl(cr), and LiOH.3LiCl(cr)) are also estimated using solubility data.