Numerical model to simulate groundwater response to ice sheet decline

The resource includes numerical modeling code developed using FlexPDE to assess ice sheet loading on groundwater head and salinity distributions in sub-ice sheet coastal aquifers. This code is based off the foundational code developed and described by Bense & Person (2008) and simultaneously solves coupled partial differential equations describing subsurface, two-dimensional, variable-density fluid flow and solute transport with mechanical loading. Included in this resource are the spin-up files and full simulation code, as well as metadata describing the simulation process.