diff --git a/watertap/property_models/multicomp_aq_sol_prop_pack.py b/watertap/property_models/multicomp_aq_sol_prop_pack.py
index 2059758a8f..c742c6956d 100644
--- a/watertap/property_models/multicomp_aq_sol_prop_pack.py
+++ b/watertap/property_models/multicomp_aq_sol_prop_pack.py
@@ -905,8 +905,7 @@ def initialize(
                         self[k].enth_mass_phase["Liq"],
                         self[k].eq_enth_mass_phase["Liq"],
                     )
-                else:
-                    self[k].total_dissolved_solids = 0
+
         # Check when the state vars are fixed already result in dof 0
         for k in self.keys():
             dof = degrees_of_freedom(self[k])