“…The MODFE approximates the following two‐dimensional partial differential equation of groundwater flow in the UFA, subject to appropriate boundary and initial conditions that represent flow through porous media (Jones & Torak, ): where t is time, ( x , y ) are the cartesian coordinate directions, H ( x , y , t ) is the hydraulic head of the USCU, h ( x , y , t ) is the aquifer hydraulic head, R ( x , y , t ) is the vertical hydraulic conductance (vertical hydraulic conductivity divided by thickness) of USCU, W ( x , y , t ) is the unit areal recharge or discharge rate (infiltration), P ( x , y , t ) is point source or sink, S ( x , y , t ) is the storage coefficient, and symmetric transmissivity is written in matrix form as …”