diff --git a/ProcessLib/HT/MonolithicHTFEM.h b/ProcessLib/HT/MonolithicHTFEM.h
index 613797b4dfe2c009728fcb8756324c2c6c685c61..3b6c268c23d7fd63cefdaaf9a682af43e93b9b0e 100644
--- a/ProcessLib/HT/MonolithicHTFEM.h
+++ b/ProcessLib/HT/MonolithicHTFEM.h
@@ -120,8 +120,8 @@ public:
             // \todo the argument to getValue() has to be changed for non
             // constant storage model
             auto const specific_storage =
-                process_data.porous_media_properties.getSpecificStorage(t, pos)
-                    .getValue(0.0);
+                solid_phase.property(MaterialPropertyLib::PropertyType::storage)
+                    .template value<double>(vars);
 
             auto const& ip_data = this->_ip_data[ip];
             auto const& N = ip_data.N;
diff --git a/ProcessLib/HT/StaggeredHTFEM-impl.h b/ProcessLib/HT/StaggeredHTFEM-impl.h
index 2d81e71ddc98ecf0e3ce80bd9909854c92f59e35..63c3024ae54f6f3d300ebdd02ba674a451ce94ba 100644
--- a/ProcessLib/HT/StaggeredHTFEM-impl.h
+++ b/ProcessLib/HT/StaggeredHTFEM-impl.h
@@ -123,9 +123,8 @@ void StaggeredHTFEM<ShapeFunction, IntegrationMethod, GlobalDim>::
         // \todo the argument to getValue() has to be changed for non
         // constant storage model
         auto const specific_storage =
-            material_properties.porous_media_properties
-                .getSpecificStorage(t, pos)
-                .getValue(0.0);
+            solid_phase.property(MaterialPropertyLib::PropertyType::storage)
+                .template value<double>(vars);
 
         auto const intrinsic_permeability = intrinsicPermeability<GlobalDim>(
             solid_phase