We developed an analytical solution for solute transport under steady-state, two-dimensional, unsaturated flow and transport conditions for the investigation of high-level radioactive waste disposal. The two-dimensional, unsaturated flow problem is treated using the quasilinear flow method for a system with homogeneous material properties. Dispersion is modeled as isotropic and is proportional to the effective hydraulic conductivity. This leads to a quasilinear form for the transport problem in terms of a scalar potential that is analogous to the Kirchhoff potential for quasilinear flow. The solutions for both flow and transport scalar potentials take the form of Fourier series. The particular solution given here is for two sources of flow, with one source containing a dissolved solute. The solution method may easily be extended, however, for any combination of flow and solute sources under steady-state conditions. The analytical results for multidimensional solute transport problems, which previously could only be solved numerically, also offer an additional way to benchmark numerical solutions. An analytical solution for two-dimensional, steady-state solute transport under unsaturated flow conditions is presented. A specific case with two sources is solved but may be generalized to any combination of sources. The analytical results complement numerical solutions, which were previously required to …

To better understand the role of groundwater-level changes on rock-slope deformation and damage, a carbonate rock slope (30 m x 30 m x 15 m) was extensively instrumented for mesoscale hydraulic and mechanical measurements during water-level changes. The slope is naturally drained by a spring that can be artificially closed or opened by a water gate. In this study, a 2-hour slope-dewatering experiment was analyzed. Changes in fluid pressure and deformation were simultaneously monitored, both at discontinuities and in the intact rock, using short-base extensometers and pressure gauges as well as tiltmeters fixed at the slope surface. Field data were analyzed with different coupled hydromechanical (HM) codes (ROCMAS, FLAC{sup 3D}, and UDEC). Field data indicate that in the faults, a 40 kPa pressure fall occurs in 2 minutes and induces a 0.5 to 31 x 10{sup -6} m normal closure. Pressure fall is slower in the bedding-planes, lasting 120 minutes with no normal deformation. No pressure change or deformation is observed in the intact rock. The slope surface displays a complex tilt towards the interior of the slope, with magnitudes ranging from 0.6 to 15 x 10{sup -6} rad. Close agreement with model for both slope surface and internal measurements …