Figure 1 should be here
2.2.2 Governing equations
The model coupled the variable saturated flow (Voss, 1984) based on
Richards’ Equation and solute transport and reaction process based on
Advection-Dispersion-Reaction Equation (1). Richards’ Equation
considered the poroelastic response of aquifer media that caused by
changing hydrostatic loading, which was the same as that of previous
studies (e.g., Reeves, 2000; Gardner and Wilson, 2006; Boutt, 2010).
\(\frac{\partial}{\partial t}\left(\theta C_j\right)-\nabla\bullet\left(\theta D\nabla C_j\right)+\nabla\bullet\left(\mathbf{q}C_j\right)=R_{j}\) (1)