where K is the hydraulic conductivity; W represents the source‒sink
term; D represents the extent of the simulation area; H represents the
groundwater level; B is the elevation of the phreatic aquifer bottom;
H1 is the initial groundwater level; q is the flow rate;
μ is the specific yield; t is the time; τ1 is the
Dirichlet condition; and τ2 is the Neumann condition.
The model is split into 50 rows and 100 columns, and according to the
Zone Budget module in Visual MODFLOW, the wetland and piedmont plain
aquifers are divided into two zones, zone 1 and zone 2, as shown in
Figure 3. This is done to quantitatively calculate the changes in
reduced base flow, induced infiltration and total water reduction due to
pumping during the wetland–aquifer interaction.