diff --git a/ebos/eclproblem.hh b/ebos/eclproblem.hh
index eb4c913ea..034f3b725 100644
--- a/ebos/eclproblem.hh
+++ b/ebos/eclproblem.hh
@@ -1748,6 +1748,9 @@ public:
     EclWellModel& wellModel()
     { return wellModel_; }
 
+    EclAquiferModel& mutableAquiferModel()
+    { return aquiferModel_; }
+
     // temporary solution to facilitate output of initial state from flow
     const InitialFluidState& initialFluidState(unsigned globalDofIdx) const
     { return initialFluidStates_[globalDofIdx]; }