A methodology to describe spatial surface flux boundary conditions for solving tailings impoundment closure water balance problems

The calculation of surface fluxes for an unsaturated tailings impoundment surface is difficult due to the fact that there is a spatially varying phreatic surface that determines the thickness of the vadose zone. This study presents a spatial flux hypothesis, which states that spatial flux boundary conditions for a generalized tailings impoundment cross-section follows a characteristic shape that is governed by the depth to the phreatic surface.; The hypothesis states that evaporation will be a minimum close to the tailings impoundment embankment wall where the depth to the phreatic surface is the greatest, and will increase to a maximum close to the pool. Inversely, infiltration will be a maximum at the embankment and will decrease to a minimum close to the pool. This study presents the methodology to calculate the spatial flux boundary functions proposed in the hypothesis, and shows how these flux boundary functions can be used as direct input for surface flux boundary conditions in multidimensional saturated/unsaturated flow seepage analysis models. The method effectively bridges the gap that currently exists between rigorous coupled soil/atmosphere one-dimensional surface flux boundary numerical models and multidimensional saturated/unsaturated flow seepage analysis. The use of the calculated spatial flux boundary functions is demonstrated though detailed evaluation modeling.; The calculation of the flux boundary function stems from the development of a technique whereby the one-dimensional SoilCover surface flux boundary model can be used to solve a two-dimensional cross section. The technique consists of a generalized non-dimensionalized tailings impoundment cross-section that comprises a beach profile and a phreatic level function. Material properties and the shape functions have been tested and calibrated through an extensive laboratory and field characterization program for the tailings. The generalized cross-section is divided into a number of equal zones and a SoilCover simulation is performed for each zone before being integrated to give a cumulative result.; The cumulative result is tested and calibrated using detailed field measurements for the tailings impoundment water balance at the Kidston Gold Mine, Australia. This cumulative result represents the spatial flux boundary function that is consistent with the spatial flux hypothesis. In summary the thesis presents a quasi-three-dimensional model for calculating surface flux boundary conditions.