Groundwater/Base Flow

The shallow aquifer contributes base flow to the main channel or reach within the subbasin. Base flow is allowed to enter the reach only if the amount of water stored in the shallow aquifer exceeds a threshold value specified by the user, aqshthr,qaq_{shthr,q}.

The steady-state response of groundwater flow to recharge is (Hooghoudt, 1940):

Qgw=8000KsatLgw2hwtblQ_{gw}=\frac{8000*K_{sat}}{L_{gw}^2}*h_{wtbl} 2:4.2.6

where QgwQ_{gw} is the groundwater flow, or base flow, into the main channel on day ii (mm H2_2O), KsatK_{sat} is the hydraulic conductivity of the aquifer (mm/day), LgwL_{gw} is the distance from the ridge or subbasin divide for the groundwater system to the main channel (m), and hwtblh_{wtbl} is the water table height (m).

Water table fluctuations due to non-steady-state response of groundwater flow to periodic recharge is calculated (Smedema and Rycroft, 1983):

dhwtbldt=wrchrg,shQgw800μ\frac{dh_{wtbl}}{dt}=\frac{w_{rchrg,sh}-Q_{gw}}{800*\mu} 2:4.2.7

where dhwtbldt\frac{dh_{wtbl}}{dt} is the change in water table height with time (mm/day), wrchrg,shw_{rchrg,sh} is the amount of recharge entering the shallow aquifer on day ii (mm H2_2O), QgwQ_{gw} is the groundwater flow into the main channel on day ii (mm H2_2O), and μ\mu is the specific yield of the shallow aquifer (m/m).

Assuming that variation in groundwater flow is linearly related to the rate of change in water table height, equations 2:4.2.7 and 2:4.2.6 can be combined to obtain:

dQgwdt=10KsatμLgw2(wrchrg,shQgw)=αgw(wrchrg,shQgw)\frac{dQ_{gw}}{dt}=10*\frac{K_{sat}}{\mu *L^2_{gw}}*(w_{rchrg,sh}-Q_{gw})=\alpha_{gw}*(w_{rchrg,sh}-Q_{gw})

2:4.2.8

where QgwQ_{gw} is the groundwater flow into the main channel on day ii (mm H2_2O), KsatK_{sat} is the hydraulic conductivity of the aquifer (mm/day),μ\mu is the specific yield of the shallow aquifer (m/m), LgwL_{gw} is the distance from the ridge or subbasin divide for the groundwater system to the main channel (m), wrchrg,shw_{rchrg,sh} is the amount of recharge entering the shallow aquifer on day ii (mm H2_2O) and αgw\alpha_{gw} is the baseflow recession constant or constant of proportionality. Integration of equation 2:4.2.8 and rearranging to solve for QgwQ_{gw} yields:

Qgw,i=Qgw,i1expαgwΔt+wrchrg,sh(1expαgwΔt)Q_{gw,i}=Q_{gw,i-1}*exp\lfloor-\alpha_{gw}*\Delta t\rfloor+w_{rchrg,sh}*(1-exp\lfloor-\alpha_{gw}*\Delta t\rfloor)

if aqsh>aqshthr,qaq_{sh}>aq_{shthr,q} 2:4.2.9

Qgw,i=0Q_{gw,i}=0 if aqshaqshthr,qaq_{sh} \le aq_{shthr,q} 2:4.2.10

where Qgw,iQ_{gw,i} is the groundwater flow into the main channel on day ii (mm H2_2O), Qgw,i1Q_{gw,i-1} is the groundwater flow into the main channel on day i1i-1 (mm H2_2O), αgw\alpha_{gw} is the baseflow recession constant, Δt\Delta t is the time step (1 day), wrchrg,shw_{rchrg,sh} is the amount of recharge entering the shallow aquifer on day ii (mm H2_2O), aqshaq_{sh} is the amount of water stored in the shallow aquifer at the beginning of day ii (mm H2_2O) and aqshthr,qaq_{shthr,q} is the threshold water level in the shallow aquifer for groundwater contribution to the main channel to occur (mm H2_2O).

The baseflow recession constant, αgw\alpha_{gw}, is a direct index of groundwater flow response to changes in recharge (Smedema and Rycroft, 1983). Values vary from 0.1-0.3 for land with slow response to recharge to 0.9-1.0 for land with a rapid response. Although the baseflow recession constant may be calculated, the best estimates are obtained by analyzing measured streamflow during periods of no recharge in the watershed.

When the shallow aquifer receives no recharge, equation 2:4.2.9 simplifies to:

Qgw=Qgw,0expαgwtQ_{gw}=Q_{gw,0}*exp\lfloor-\alpha_{gw}*t\rfloor if aqsh>aqshthr,qaq_{sh}>aq_{shthr,q} 2:4.2.11

Qgw,i=0Q_{gw,i}=0 if aqshaqshthr,qaq_{sh} \le aq_{shthr,q} 2:4.2.12

where QgwQ_{gw} is the groundwater flow into the main channel at time tt (mm H2_2O), Qgw,0Q_{gw,0} is the groundwater flow into the main channel at the beginning of the recession (time tt=0) (mm H2_2O), αgw\alpha_{gw} is the baseflow recession constant, and t is the time lapsed since the beginning of the recession (days), aqshaq_{sh} is the amount of water stored in the shallow aquifer at the beginning of day ii (mm H2_2O) and aqshthr,qaq_{shthr,q} is the threshold water level in the shallow aquifer for groundwater contribution to the main channel to occur (mm H2_2O). The baseflow recession constant is measured by rearranging equation 2:4.2.11.

αgw=1Nln[Qgw,NQgw,0]\alpha_{gw}=\frac{1}{N}*ln[\frac{Q_{gw,N}}{Q_{gw,0}}] 2:4.2.13

where αgw\alpha_{gw} is the baseflow recession constant, NN is the time lapsed since the start of the recession (days), Qgw,NQ_{gw,N} is the groundwater flow on day NN (mm H2_2O), Qgw,0Q_{gw,0} is the groundwater flow at the start of the recession (mm H2_2O).

It is common to find the baseflow days reported for a stream gage or watershed. This is the number of days for base flow recession to decline through one log cycle. When baseflow days are used, equation 2:4.2.13 can be further simplified:

αgw=1Nln[Qgw,NQgw,0]=1BFDln[10]=2.3BFD\alpha_{gw}=\frac{1}{N}*ln[\frac{Q_{gw,N}}{Q_{gw,0}}]=\frac{1}{BFD}*ln[10]=\frac{2.3}{BFD} 2:4.2.14

where αgw\alpha_{gw} is the baseflow recession constant, and BFDBFD is the number of baseflow days for the watershed.

Last updated