rho 1e3[kg/m^3] Fluid density thetas_1 0.339 Porosity, layer 1 thetas_2 0.399 Porosity, layer 2 thetar_1 0.001 Residual saturation, layer 1 thetar_2 0.001 Residual saturation, layer 2 Ss_1 (thetas_1-thetar_1)/(1[m]*g_const*rho) Storage coefficient, layer 1 Ss_2 (thetas_2-thetar_2)/(1[m]*g_const*rho) Storage coefficient, layer 2 Ks_1 0.454[m/d] Saturated hydraulic conductivity, layer 1 Ks_2 0.298[m/d] Saturated hydraulic conductivity, layer 2 alpha_1 1.39[1/m] Parameter alpha, layer 1 alpha_2 1.74[1/m] Parameter alpha, layer 2 n_1 1.6 Parameter n, layer 1 n_2 1.38 Parameter n, layer 2 Hp0 0.01[m] Pressure head in the ring N0 0.01*Ks_1*rho Leak from the base rhob 1400[kg/m^3] Bulk density kp 1e-4[m^3/kg] Partition coefficient Dl 3.74e-3[m^2/d] Diffusion coefficient phil 0.05[1/d] Decay rate in liquid phip 0.01[1/d] Decay rate in solid alphar 0.005[m] Longitudinal dispersivity alphaz 0.001[m] Transverse dispersivity c0 1[kg/m^3] Solute concentration in the ring