diff --git a/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.hpp b/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.hpp
index 7af31260e..86702f76f 100644
--- a/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.hpp
+++ b/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.hpp
@@ -44,6 +44,8 @@ namespace Opm {
 
         std::unordered_map<size_t, const NumericalAquiferCell*> allAquiferCells() const;
 
+        std::unordered_map<size_t, double> aquiferCellVolumes() const;
+
         std::vector<NNCdata> aquiferNNCs(const EclipseGrid& grid, const FieldPropsManager& fp) const;
 
         std::unordered_map<size_t, AquiferCellProps> aquiferCellProps() const;
diff --git a/src/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.cpp b/src/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.cpp
index 4ecec978d..a1e3a3cce 100644
--- a/src/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.cpp
+++ b/src/opm/parser/eclipse/EclipseState/Aquifer/NumericalAquifer/NumericalAquifers.cpp
@@ -168,4 +168,13 @@ namespace Opm {
         }
         return nncs;
     }
+
+    std::unordered_map<size_t, double> NumericalAquifers::aquiferCellVolumes() const {
+        std::unordered_map<size_t, double> cell_volumes;
+        const auto aquifer_cells = this->allAquiferCells();
+        for (const auto& [global_index, cell] : aquifer_cells) {
+            cell_volumes.insert(std::make_pair(global_index, cell->cellVolume()));
+        }
+        return cell_volumes;
+    }
 }
\ No newline at end of file