diff --git a/ChemistryLib/PhreeqcIOData/AqueousSolution.h b/ChemistryLib/PhreeqcIOData/AqueousSolution.h
index c73d407179895d6a4fb2047f3985743654467b1d..b802479b1c00c2bc9f4ca0b0bdcbefdc2e59b125 100644
--- a/ChemistryLib/PhreeqcIOData/AqueousSolution.h
+++ b/ChemistryLib/PhreeqcIOData/AqueousSolution.h
@@ -16,7 +16,6 @@
 #include <vector>
 
 #include "MathLib/LinAlg/GlobalMatrixVectorTypes.h"
-#include "MathLib/LinAlg/MatrixVectorTraits.h"
 #include "Output.h"
 
 namespace MeshLib
@@ -45,12 +44,9 @@ struct AqueousSolution
     AqueousSolution(double temperature_, double pressure_,
                     MeshLib::PropertyVector<double>* pe_, double const pe0_,
                     std::vector<Component>&& components_,
-                    ChargeBalance charge_balance_,
-                    std::size_t const num_chemical_systems_)
+                    ChargeBalance charge_balance_)
         : temperature(temperature_),
           pressure(pressure_),
-          pH(MathLib::MatrixVectorTraits<GlobalVector>::newInstance(
-              num_chemical_systems_)),
           pe(pe_),
           pe0(pe0_),
           components(std::move(components_)),
diff --git a/ChemistryLib/PhreeqcIOData/CreateAqueousSolution.cpp b/ChemistryLib/PhreeqcIOData/CreateAqueousSolution.cpp
index d6cf042a27a273af8da334c9998fe73bcd842b41..ee5069eca08716268c8c12cfbc9bd268af7ca72a 100644
--- a/ChemistryLib/PhreeqcIOData/CreateAqueousSolution.cpp
+++ b/ChemistryLib/PhreeqcIOData/CreateAqueousSolution.cpp
@@ -41,13 +41,8 @@ std::unique_ptr<AqueousSolution> createAqueousSolution(
 
     auto charge_balance = createChargeBalance(config);
 
-    return std::make_unique<AqueousSolution>(temperature,
-                                             pressure,
-                                             pe,
-                                             pe0,
-                                             std::move(components),
-                                             charge_balance,
-                                             mesh.getNumberOfBaseNodes());
+    return std::make_unique<AqueousSolution>(
+        temperature, pressure, pe, pe0, std::move(components), charge_balance);
 }
 }  // namespace PhreeqcIOData
 }  // namespace ChemistryLib