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)