diff --git a/dumux/porousmediumflow/richards/localresidual.hh b/dumux/porousmediumflow/richards/localresidual.hh
index dac5b3d933bef0fd1091c3f281feaa8135064d3d..6e58d391e27214bd8ace91eeb3419cb73d4c73ab 100644
--- a/dumux/porousmediumflow/richards/localresidual.hh
+++ b/dumux/porousmediumflow/richards/localresidual.hh
@@ -196,7 +196,7 @@ public:
         static_assert(!FluidSystem::isCompressible(0),
                       "richards/localresidual.hh: Analytic Jacobian only supports incompressible fluids!");
 
-        const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity();
+        const auto poreVolume = Extrusion::volume(scv)*curVolVars.porosity()*curVolVars.extrusionFactor();
         static const auto rho = curVolVars.density(0);
 
         // partial derivative of storage term w.r.t. p_w