Simulation of Leachate Generation from a Waste Rock Dump in Kiruna Using HYDRUS-1D

• Main Author: Atmosudirdjo, Aryani
• Format: Bachelor Thesis
• Language: English
• Published: Uppsala universitet, Institutionen för geovetenskaper 2019
• Subjects: Hydrus-1D; matrix flow; leachate; unsaturated flow; waste rock dump; Oceanography; Hydrology and Water Resources; Oceanografi; hydrologi och vattenresurser; Kiruna; Northern Sweden
• Online Access: http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-379225

The percolation of water through waste rock dumps at mine sites can lead to the production of a leachate with high concentrations of dissolved metals, sulfate and nitrogen compounds. It is important to understand how water flows in waste rock dumps in order to predict the environmental impact of this leachate on recipients. The dynamics of percolation and leachate discharge are controlled by climatological conditions at the site, where relatively large flows in northern Sweden correspond to snowmelt during late Spring. Rock dumps are often tens of meters in height, resulting in an unsaturated water flow system through heterogeneous material. Hence, the simulation of leachate generation requires an accurate representation of the subsurface materials as well as the flow processes, where water flow in waste rock dumps is dominated by matrix flow with macropore flow being of secondary importance. Matrix flow is rather slow and may thus potentially yield relatively high concentrations of contaminants in the leachate, in response to precipitation and snow melt. This study uses Hydrus-1D to predict leachate generation from a small-scale waste rock dump in Kiruna in terms of discharge magnitude and timing. The 3-dimensional geometry of the waste rock dump is approximated by summing simulations from 1225 one-dimensional columns of different length, with a surface area of 1 m2 each. There are four output parameters that are compared between the model results and measured data: snow accumulation, water content, temperature, and discharge. There are some discrepancies between the model results and field measurements, most likely due to uncertainties in the input parameters (especially waste rock properties), limitations in the Hydrus-1D model (i.e. freeze-thaw dynamics), and assumptions that are used in constructing the conceptual model. For better agreement between model results and measured data, a new modelling approach is recommended, potentially using a different program than Hydrus-1D.