6.1. Verification Modeling
The ability of the VADOSE/W for predicting the water balance was investigated by Walter and Dubreuilh (2007) [
47]. They compared the model outputs with the field measurements. The results revealed that, in some instances, VADOSE/W predicted the water balance very close to the field measurements. However, in other instances, the model results did not correlate well to the field conditions. The reasons for the numerical code not being able to predict the field conditions could be due to uncertainty in the hydraulic properties, vegetation water use parameters, heterogeneous soil conditions, inaccurate representation of rainfall events, and/or unexpected generation of high surface runoff. Considering the uncertainty in the modeling results, it is generally suggested to verify/calibrate the model prior to its use for future modelling considerations.
Field measurements to verify/calibrate the model for the mine site are not available. Therefore, an attempt was made to assess the suitability of the current modeling effort by reproducing the results of previous modeling studies carried out at Detour Lake Mine. Dobchuk (2002) [
30] predicted the water balance of the tailings covers at the site using the SoilCover numerical code. The modelling results were never quantitatively compared to the site condition, but were reported to provide a qualitative assessment of the tailings covers. The objective of the current study is to assess the effects of climate change on tailings covers; therefore, verifying the model with Dobchuk (2002) [
30] modeling results will provide a consistent approach to the previous research carried out at the site. To accomplish this, water balances for coarse and fine tailings covers are compared in
Figure 6.
Figure 6a shows a comparison of the water balance for current and previous models for the homogenous coarse tailings cover (CC profile). The coarse tailings cover shows the annual net infiltration (
) of about 2% less in the current model to that of Dobchuk (2002) [
30]. The current model predicts the annual surface runoff (
) of 11 mm (2.3%) whereas the previous results showed minimal surface runoff. The difference in annual actual evaporation (
) between both models is insignificantly small (~0.01%). About 5% increase in
is recorded in the current model as compared to the model by Dobchuk (2002) [
30].
The water balance comparison for the fine tailings cover is shown in
Figure 6b. For the fine cover, both models show that the predicted annual
is significantly more than the annual precipitation; therefore, the annual
is negative, indicating net water loss conditions at the ground surface. A difference of 6% in
and 7% in
was observed between the current modeling results and those reported by Dobchuk (2002) [
30]. The comparisons made in
Figure 6a,b reveal that the water balance for the current models (coarse and fine tailings covers) are generally comparable to that of Dobchuk (2002) [
30] models. The minor variations in the results could be due to the use of different numerical codes.
6.3. Water Balance at the Ground Surface
The water balance at the ground surface quantifies the availability of water and evaporative demand together with quantities of water that enters or exits the ground surface. In this research, the components of water balance which enter the ground surface are considered as positive while the components which exit the ground surface are considered negative. The exception is runoff, which is also considered positive. The deep percolation, which is the flux across the cover–tailings interface, is considered positive if the flow is from the tailings into the cover. Similarly, is expressed as a negative quantity if it flows downwards from the cover to the tailings.
Figure 7 shows the water balance for the homogeneous coarse profile (CC) for dry and wet year climates for the base climate (BC) and future climates (FC-1 and FC-2). Dry and wet years were identified based on the lowest and highest
for the respective 30-year climate datasets. The annual cumulative precipitation for dry years increases 19% and 13% for future climates FC-1 and FC-2 respectively in comparison to the base climate (
Figure 7a). However, for the wet years, the annual cumulative precipitation for FC-1 is 21% lower than the baseline condition, while for FC-2 it is 9% greater than the baseline value (
Figure 7b).
Potential evaporation is the total quantity of water that can be evaporated if the unlimited supply of water is available at the ground surface. The potential evaporation for dry years is lowest for BC followed by FC-1 and FC-2. The increase in the potential evaporation relative to the base climate is 12% and 30% for FC-1 and FC-2, respectively (
Figure 7a). However, for wet years, the lowest
can be observed for FC-1, followed by BC and FC-2.
Actual evaporation is related to
but also depends on the availability of water (due to precipitation) and potential moisture difference between earth surface and atmosphere. Although cumulative precipitation and the
of homogeneous coarse profile for dry years varies for base and future climates, the cumulative actual evaporation is almost similar (within 1% difference). This is due to the low retention and high conduction behavior of the coarse material. Conversely, the actual evaporation varies considerably (up to 31%) during wet year for base and future climates. The reason for larger variations in cumulative
during wet years could be due to the variations in the amount of precipitation. For example, looking at the precipitation data in
Figure 7b, the FC-2 has more extreme events especially at the start. This implies that there is more availability of water in a shorter period of time, although the evaporative demand for that period is not too different among the three. The lowest cumulative actual evaporation was observed during FC-1, followed by BC and FC-2 and this trend is in line with the precipitation trends of the wet year.
Net infiltration quantifies the amount of water entering or leaving the ground surface. It is calculated by subtracting the sum of
and
from the precipitation. The cumulative net infiltration during a dry year generally follows the precipitation trends for the homogeneous coarse profile (
Figure 7a) with lowest during base climate, followed by FC-2 and then FC-1. The homogeneous coarse cover profile shows positive net infiltration during dry and wet years of representative future climates. The profile (CC) also shows water gain conditions in wet year of BC. However, water deficit conditions were observed during the dry year of BC. The lowest cumulative net infiltration of wet year was observed during climate FC-2, followed by BC and FC-1, which inversely follows the precipitation trends of the wet years of their respective climate ensembles. The reason for such trends (which inversely follow the precipitation) is due to actual evaporation values, which are highest during FC-2, followed by BC and FC-1.
The infiltrated water flow through gravity and entering the tailings is generally termed as deep percolation (
). In addition to gravity, another competing process is the capillary action, in which water from the saturated portion of the tailings can rise and enter the cover. Since the homogenous coarse profile has low air entry value with little expected capillary rise; therefore, a major portion of infiltrating water during all climates flow through the cover and enters the tailings (
Figure 7a,b).
The water balance for the homogeneous fine material profile of dry and wet years during base and future climates is shown in
Figure 8. The cumulative precipitation and potential evaporation of dry and wet years for the homogeneous fine profile during all the climates are identical to their respective climate ensembles of the coarse homogeneous profile, as these (
and
) are independent of the cover type.
However, actual evaporations are significantly higher for the homogenous fine profile as compared to the coarse profile for dry and wet years for all the climate ensembles. The relative higher cumulative for the fine profile is due to its high retention and low drainage behavior. The variation in actual evaporation of dry year during for all climates is relatively small (within 5% of BC). However, the difference in cumulative for wet years are high, with the highest of 818 mm of FC-2 (33% more than BC), followed by BC (611 mm) and FC-1 (443 mm). Again, the reasons for the variations in the could be the variations in the amount of precipitation in wet years for base and future climates.
The net infiltration of the homogeneous fine profile for dry years showed water deficit conditions for all the representative design climate ensembles. The water deficit conditions are due to increasing high actual evaporation and lower cumulative precipitation. Additionally, the higher air entry value (50 kPa) of the fine profile keeps the saturation at a higher level above the groundwater table due to capillarity which results in generation of higher
(
Figure 8a). During wet years, the net annual precipitation is more than the annual actual evaporation; therefore, the wet years of base and future climates show water gain conditions (
Figure 8b).
For the case of a homogeneous fine profile, a greater capillary rise due to the higher air entry value of the fine tailings is expected. At the beginning of the active year, the flow of water in the cover is predominantly controlled by gravity drainage. As time passes, water flows across the cover–tailings interface in the upward direction due to capillary rise. During dry years, due to higher and lower , the water flows from the tailings to the cover due to capillary rise will result in net gain to the cover storage. Conversely, during wet years, the water flows across the interface from cover to the tailings for BC and FC-1. This is due to the higher amount of precipitation ( > ( = )) during wet years. However, for wet years during FC-2, the cover experiences net gain in its storage. This is due to the high (comparable to ), which creates potential continuous water deficit conditions at the ground surface which in turn were compensated by water flow from the tailings to the cover due to capillary action.
The water balance results for coarse cover profile over fine tailings (CF) showed similar trends to that of profile CC, and the fine cover over coarse tailings profile shows water balance similar to that of profile FF; therefore, they are not discussed any further. In addition to the above mentioned water balance results for FC-1 and FC-2, two additional representative climate ensembles were also analyzed, and their water balance results were found to be similar to those reported here. Details can be found in [
37].
6.4. Annual Variations in Water Storage and Oxygen Flux
The box and whisker plots for change in annual storage, Δ
S (difference in
and
) and oxygen flux are presented in
Figure 9 for the both fine and coarse tailing covers for all the climates. The change in storage of the cover layer is determined by estimating the quantity of water entering or leaving the top (ground surface) and bottom boundary (cover tailings interface) of the cover. As mentioned earlier, the meteoric water that enters or leaves the cover surface at the top boundary is called
.
A portion of this water can make its way into the tailings and is known as . In some instances, for shallow groundwater table depth and tailings with high AEV, water from the tailings can also enter the cover via capillary action. The difference between the and predicts the changes in the cover storage (saturation). The positive value of ΔS represents the gain in storage, while the negative predicts the loss of in cover storage. The value of ΔS is related to the ingress of the oxygen in the sense that it affects the cover saturation which is the controlling parameter for oxygen transport. It should be noted that the cumulative value of ΔS indicates the net water loss or gain during the period in which cover is thawed and moisture transfer can take place across the soil atmosphere boundary.
Figure 9a shows the box and whisker plots of the annual Δ
S and oxygen flux for all the cover profiles. The results indicate that for the CC profile, the mean annual Δ
S (30 years) are −24, −22 and −42 mm for the base, FC-1 and FC-2 future climates, respectively. Therefore, it can be concluded that the difference in mean values for base and future climates is minimal. However, the year to year variation in Δ
S for the future climates is much larger than the base climate. For example, the difference in maximum and minimum annual Δ
S for the CC profile is 219, 256 and 403 mm for base, FC-1 and FC-2 climates, respectively. This increasing trend (in the range) is reflective of more year to year variation and the existence of extreme weather years in the future climate ensembles. More variation can be seen for future climate ensemble FC-2. Similar observations regarding the changes in Δ
S can also be made for the CF profile and Δ
S values for the CC and CF profiles are very similar.
Figure 9a also shows the box and whisker plots for the annual oxygen flux for the coarse tailing covers. For cover profile CC, it can be observed that changes in cover saturation result in variations in yearly flux values. Larger year to year variation can be observed for FC-1, while larger mean and median increases can be observed for future climate FC-2. For cover profile CF, it can be observed that oxygen flux values are much smaller than the values for CC profile. This result seems surprising, as Δ
S values for both profiles are quite similar.
Figure 10 shows the time history of saturation for the CC and CF profiles for the base climate. This figure shows that the saturation profiles for CC and CF are quite different. For example, for the CF cover profile, the reactive tailings remain at saturations higher than 90%. This is in contrast to the CC profile where most of the reactive tailings remain at saturations levels between 40 and 50%. Additionally, it can also be observed that the cover of the CF profile has much higher saturations than for the CC profiles. The higher tailings and cover saturation in the CF profile is due to higher
AEV of the fine reactive tailings and which is responsible for lower oxygen flux in the CF profile. Therefore, it can be concluded that the performance of coarse cover lying over fine reactive tailings is less dependent on the climate and will be less prone to future climate changes.
Figure 9b shows the changes in saturation and oxygen flux values for the FF and FC profiles. The change in saturation for these profiles are similar to those predicted for the coarse cover profiles. It can be observed that with the changing climate, larger variations in year to year in saturation can be observed. It can also be observed that this variation is more pronounced for fine cover lying over fine tailings. The oxygen flux values for the FF and FC profiles shown in
Figure 9b indicate that oxygen flux increases considerably for the future climates FC-1 and FC-2 in comparison to the base climate. This increase is greater for FC-2 as compared to FC-1. It can also be observed that FC profiles show more sensitivity to future climates than the FF profile. This is consistent with our observation from coarse cover profiles, where cover overlying fine tailings showed less dependence on climate. Based on the observations made for all cover profiles, it can be concluded that larger quantities of oxygen will enter the fine covers with changing climates.
6.5. Cumulative Oxygen Fluxes
Figure 11 shows a comparison of the 30-year cumulative oxygen flux values for all profiles predicted at the cover surface and the cover–tailings interface. The coarse tailing covers (CC and CF profiles) reveal a smaller variation in cumulative oxygen flux for all climates. The variations in the surface oxygen flux for future climates are shown as a percentage change relative to the base climate in
Figure 12. The results indicate that these changes are minimal. These further reinforce earlier observations that performance of coarse covers will not be impacted by the changing climate. It should be noted that based on the results presented in the previous section for CC cover profile some variation in yearly oxygen flux values is expected for changing climate. However, in the cumulative sense, these only account for 2% and 5% change in relation to the base climate.
The 30-year cumulative oxygen flux values for the fine cover profiles in
Figure 11 indicate considerable differences compared to the coarse cover profiles. For fine cover profiles, the oxygen flux at the tailings–cover interface is negligible. This is not a surprising result, as the presence of fine material in tailings and/or cover is expected to maintain higher levels of saturation resulting in minimal oxygen ingress.
Figure 13 shows the saturation time history of the FF and CF profiles for 30 years of base climate. The saturation time history of the FF profile shows that saturation of the fine acid generating tailings never falls below 85%. This is due to the location of the groundwater table depth and high
AEV of the tailings. It can also be observed that the cover comprising of the desulfurized fine tailings is also very effective in maintaining higher saturations and the bottom 30 cm of the cover always remains at saturation in excess of 85%. The time history of saturation for the FC profile is also shown in the same figure. It can be observed that for most of the coarse acid generating tailings, the saturation varies from 20% to 60%. However, the cover saturation near the cover tailings interface remains above 85%. There are only a few instances during the 30-year period when saturation at the interface is 70% or slightly lower.
Figure 11 and
Figure 12 also show the oxygen flux values at the cover surface for fine covers. The flux values increase with changing climate. Larger changes in the 30-year cumulative flux values for FF and FC profiles can be observed for the FC-2 climate ensemble. In relation to the base climate, the flux values for the FF cover profile increase by 18 and 53% for the FC-1 and FC-2 climates, respectively. Similarly, for the FC cover profile the oxygen flux at the cover surface increases by 24% and 65% for the FC-1 and FC-2 climates, respectively. The increase in surface oxygen flux for fine cover profiles is related to the higher effective reactive coefficient (
) of the fine tailings as compared to the coarse. The
value of fine tailings is an order of the magnitude larger than for the coarse tailings material (
for fine 44 year
−1 and
for coarse is 3.44 year
−1). The other reason could be due to the increased evaporative demand because of the changing climate. The fine cover will retain more water in the near surface layer, increasing actual evaporation, resulting in a decrease in cover saturations with increased oxygen ingress at the surface.
The results reveal that the FC profile has the largest variation in the cumulative water and the oxygen flux and is observed as more sensitive to the climate variations (
Figure 9b and
Figure 11). The profile CF represents the most accurately as-built configuration of tailings cover over reactive tailings at Detour Lake Mine [
10]. Therefore, the CF and
FC profiles are selected further to assess the temporal effects of climate changes on these profiles.
Time histories of saturation and oxygen concentration versus depth for the CF profile for all design climates are shown in
Figure 14.
Figure 14a–c show the saturation depth plots whereas the
Figure 14d–f show concentration depth plots for the three climates.
Figure 14a–c clearly show that the tailings are at nearly 100% saturation, except the region close to the interface, for all climates. The saturations at the interface are consistently between 60 and 70% for CF profile for all climates. There are also some instances when the saturation becomes more than 85% at the interface. The time history of saturation with depth in the cover is generally in the range between 30 and 50% during BC, FC-1 and FC-2, with some dry spells when the saturation in the covers drops in the range between 10% and 20%. These low saturations in the cover are more obvious near the cover surface. These low saturation conditions are more frequent in FC-1 and FC-2 (every year) and the depth of these low saturation contours can reach a depth of up to 1.75 m below the cover surface, especially during FC-2. It can be concluded that the cover saturation with depth is decreasing with change in climate and the saturation starts decreasing from the surface, and as the climate becomes drier, the low saturation contour progresses further down.
Oxygen concentration inversely follows the saturation plots (
Figure 14d–f). As the tailings remain nearly saturated, the oxygen ingress in the reactive tailings is negligible for all climates. Almost one-fourth of the upper portion of the cover receives high oxygen concentration (~280 g/m
2) during BC. The high oxygen concentration contours almost reach the upper half of the cover for FC-1 climate and go beyond this depth that for the climate FC-2. However, there are instances when the oxygen concentration drops in the range between 150 and 250 g/m
2 near the cover surface. The oxygen concentration in the cover near the interface is generally in the range between 100 and 200 g/m
2 for base climate and between 200 to 250 g/m
2 for future climates, respectively. The reason for the high oxygen concentration in the cover near the interface during the future climate is due to lower saturation near the interface. It can be concluded that no sufficient oxygen will pass through the cover in to the tailings, but with climate change conditions, a greater presence of oxygen near the interface is expected.
The 30-year time histories for saturation and oxygen concentration versus depth for the FC profile of all the design climates are shown in
Figure 15.
Figure 15a–c show the saturation depth plots, whereas
Figure 15d–f show concentration depth plots for BC, FC-1 and FC-2, respectively. In this figure, it can be observed that the saturation in the tailings increases with depth, and that the tailings become saturated at a depth of 4.0 m (which is groundwater table). The time saturation plot also shows that the saturation in the tailings near the interface is between 20 and 30% for all climates (
Figure 15a–c). However, there are instances when the saturation increases in the range between 30% to 50% and these cases are more obvious for BC. The fine cover of the profile FC is expected to have high saturation in the cover layer. The upper portion of the cover up to about 0.6 m experiences a wide range of saturation conditions (between 20% and 95%). Beyond this depth, the cover saturation is normally in the range between 60% and 95 %. The saturation in the cover near the interface is generally higher than 85%. However, there are instances when the saturation over the interface is below 85% and these instances are more common during FC-2.
As the saturation over the interface remains 85% in most of the instances, therefore, the oxygen concentration in the tailings is extremely low (
Figure 15d–f). However, in the instances when the saturation at the interface is below 85%, oxygen concentrations between 0 and 1 g/m
2 can be observed in the tailings for all climates. For FC-1 and FC-2, there are instances when the concentration in the tailings increases to a range between 5 and 10 g/m
2. Moreover, there are few years during FC-2 (year 3, 7, and 23) when the concentration goes past this range up to 20 g/m
2. This is due to the lower saturation near the interface during this time of the year. Similar to the saturation plots, the oxygen concentration in the upper portion (0 to 0.6 m approximately) is highly variable, falling in the range between 100 and 280 g/m
2. Beyond this depth, the oxygen concentration is generally in the range between 5 and 100 g/m
2. Overall, the fine cover maintains high saturation during 30-year time period near the interface, which results in the low ingress of oxygen into the tailings in all climates.