diff --git a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
index da90cbdc39ed2c1c2f197319959cce67ae1814dd..421bc2331ee3bca9d168a894a6c6990158de57f4 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
@@ -176,7 +176,7 @@ void LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::
                     t, dt);
 
         auto const porosity =
-            medium->property(MaterialPropertyLib::PropertyType::porosity)
+            solid_phase.property(MaterialPropertyLib::PropertyType::porosity)
                 .template value<double>(vars, pos, t, dt);
         auto const storage =
             solid_phase.property(MaterialPropertyLib::PropertyType::storage)
diff --git a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
index 15aa35995eda92fa617ada5e63bc33b9473d57d9..d8afe847f51ad874dcecdd5793fc87216864dc52 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
+++ b/ProcessLib/LiquidFlow/LiquidFlowProcess.cpp
@@ -32,7 +32,6 @@ void checkMPLProperties(
 
     std::array const requiredPropertyMedium = {
         MaterialPropertyLib::reference_temperature,
-        MaterialPropertyLib::PropertyType::porosity,
         MaterialPropertyLib::PropertyType::permeability};
 
     std::array const requiredPropertyLiquidPhase = {
@@ -40,6 +39,7 @@ void checkMPLProperties(
         MaterialPropertyLib::PropertyType::density};
 
     std::array const requiredPropertySolidPhase = {
+        MaterialPropertyLib::PropertyType::porosity,
         MaterialPropertyLib::PropertyType::storage};
 
     for (auto const& element : mesh.getElements())