Although SWAT+ does not currently print groundwater height in the output files, the water table height is updated daily by the model. Groundwater height is related to groundwater flow by equation 2:4.2.6.
Qgw​=Lgw2​8000∗Ksat​​∗hwtbl​=108000∗μ​∗μ∗Lgw2​10∗Ksat​​∗hwtbl​=800∗μ∗αgw​∗hwtbl​ 2:4.2.19
where Qgw​ is the groundwater flow into the main channel on day i (mm H2​O), Ksat​ is the hydraulic conductivity of the aquifer (mm/day), Lgw​ is the distance from the ridge or subbasin divide for the groundwater system to the main channel (m), hwtbl​ is the water table height (m), μ is the specific yield of the shallow aquifer (m/m), and αgw​ is the baseflow recession constant. Substituting this definition for Qgw​ into equation 2:4.2.9 gives
hwtbl,i​=hwtbl,i−1​∗exp[−αgw​∗Δt]+800∗μ∗αgw​wrchrg​∗(1−exp⌊−αgw​∗Δt⌋)​ 2:4.2.20
where is the water table height on day (m), is the water table height on day (m), is the baseflow recession constant, is the time step (1 day), is the amount of recharge entering the aquifer on day (mm HO), and is the specific yield of the shallow aquifer (m/m).
Table 2:4-1: SWAT+ input variables used in shallow aquifer calculations.
.gw
: Threshold water level in shallow aquifer for revap (mm HO)
REVAPMN
.gw
: Revap coefficient
GW_REVAP
.gw
: Aquifer percolation coefficient
RCHRG_DP
.gw
: Specific yield of the shallow aquifer (m/m)
GW_SPYLD
.gw
δgw​: Delay time for aquifer recharge (days)
GW_DELAY
.gw
aqshthr,q​: Threshold water level in shallow aquifer for base flow (mm H2​O)
GWQMN
.gw
αgw​: Baseflow recession constant
ALPHA_BF