Ground water depth formula (Subsidence Overlay)

At the start of a simulation, the ground water depth is initialized with the ground water depth geotiff (if provided) and is optionally overwritten by (managed) water level areas' water level. Furthermore, the ground water level can be managed with drainages (provided as underground buildings), either actively or passively.

Additionally the terrain height can change due to subsidence that occurred in previous years and due to actions taken that raised the terrain. Managed water level areas can react to these changes when indexation is configured. For indexation, see indexation formula.

Ground water level managed by drainages

$d_{gw,y}={\begin{cases}d_{d,active},&{\text{if }}n{\text{ drainage is larger than 0}}\\max(0.0,d_{a}-\Delta {d_{a,y}}+d_{d,passive},&{\text{if }}n{\text{ drainage is smaller than 0}}\\max(0.0,d_{gw}+\Delta {d_{a,s,y}})&{\text{ otherwise }}\end{cases}}$

where

• $d_{gw,y}$ is the calculated ground water depth of a grid cell at year y.
• $d_{d,active}$ is the ground water depth actively maintained by a drainage
• $d_{a}$ is the ground water depth of the managed water level area.
• $\Delta {d_{a,y}}$ is the area water level adjustment for year y.
• $d_{d,passive}$ is the ground water level passively maintained by drainages.
• $\Delta {d_{a,s,y}}$ is the ground water level increase calculated based on the are water level adjustment due to indexation and the subsidence in the previous year.

Ground water level managed by water level areas

The effect of changes of managed surface water levels (by water level areas) on the ground water level. The x-axis indicates the relative depth of the managed surface water level. The y-axis indicates the meters of change to the ground water level, per meter change in surface water level. The smaller the relative depth is, the less the ground water level equalizes with changes in the managed surface water level.

The ground water level increase is calculated based on the indexation on the managed water level of the surface water and the calculated subsidence of the previous year. This adjustment is rarely the adjustment of the managed water level. Therefore, the following formula is used to calculate the increase of the ground water level.

$\Delta {d_{a,s,y}}=\Delta {d_{low}}+\Delta {d_{mid}}+\Delta {d_{high}}$

$\Delta {d_{inc}}=-\Delta {d_{a,y}}-s_{y-1}$

$d_{a,adj}=d_{a}+\Delta {d_{inc}}$

$\Delta {d_{low}}=0.6\cdot (\max(0.0,min(0.6,d_{a,adj})-max(0.0,min(0.6,d_{a}))$

$\Delta {d_{mid}}={\frac {\max(0.6,min(1.0,d_{a,adj}))+max(0.6,min(1.0,d_{a}))}{2.0}}\cdot (\max(0.6,min(1.0,d_{a,adj}))-max(0.6,min(1.0,d_{a})))$

$\Delta {d_{high}}=1.0\cdot (\max(0.0,d_{a,adj}-max(0.0,d_{a})$