E ﬀ ects of Capping Strategy and Water Balance on Salt Movement in Oil Sands Reclamation Soils

: The success of oil sands reclamation can be impacted by soil salinity depending on the materials used for soil reconstruction and the capping strategies applied. Using both a greenhouse-based column experiment and numerical modeling, we examined the potential pathways of salt migration from saline groundwater into the rooting zone under di ﬀ erent capping strategies (the type and the thickness of the barrier layer) and water balance scenarios. The experimental results showed that there would be salinity issues in the cover soil within several growing seasons if there was a shallow saline groundwater table and if the soil was not properly reconstructed. The thickness of the barrier layer was the most signiﬁcant factor a ﬀ ecting the upward movement of saline groundwater and salt accumulation in the cover soil. The suitable thickness of the barrier layer for preventing the upward movement of saline groundwater and salt accumulation in the cover soil for each material varied. A numerical simulation for a 15-year period further indicates that, when the cover soil was 50 cm of peat-mineral soil mix and when wet, dry, or normal climatic conditions were considered, the minimum barrier thickness to restrain salt intrusion into the cover soil in the long term was about 75 or 200 cm for coarse tailings sand or overburden barrier material, respectively. In view of the above, to minimize salt migration into the rooting zone and ensure normal plant growth, oil sands reclamation should consider salt migration when designing soil capping strategies.


Introduction
In the Athabasca Oil Sands Region (AOSR) of Canada, surface mining for oil sands have completely altered the landscape from boreal forests and wetlands into mine pits, dumps, and tailings ponds. By 2017, the disturbed lands already covered an area of 895 km 2 , which account for about 0.2% of Alberta's boreal forest [1]. As regulated by the government, the disturbed lands need to be reclaimed. Taking restoration of ecological function and bio-diversity as a priority, all the reclamation practices are required to re-establish the functional boreal ecosystems with equivalent capability prior to disturbance [2]. Meanwhile, the mining and processing operations produce large volumes of mining by-products that must be disposed of properly [3]. Soil covers are used to cap disturbed lands to meet equivalent land capability. These reconstructed soil profiles usually have a top layer (also known as cover soil) that provides nutrients for vegetation and an underlying barrier layer (subsoil) to insulate the mining wastes beneath as well as provide more relatively clean media suitable for root growth [3,4].

Experimental Design
A greenhouse-based column experiment was set up at Alberta Innovates-Technology Futures (now InnoTech Alberta) in Vegreville, AB, Canada in summer 2013. It became operational in October the same year and was destructively sampled in May 2015 after 3 simulated growing seasons (October 2013-January 2014, April 2014-September 2014, and January 2015-April 2015). A completely randomized block design was used in the experiment with barrier material type, barrier thickness, and water balance as treatments. It examined 2 types of barrier materials (coarse-textured tailings sand (TS) vs. fine-textured overburden (OB)), barrier thickness at 4 depths (0, 20, 50 and 100 cm), and water balance at 2 levels (simulated dry vs. simulated rainy condition for the oil sands mining region in Canada).
The two barrier materials with contrasting soil textures, TS and OB, and the cover soil material, peat-mineral soil mix (PMM), were all sourced from the Canadian Natural Resources Limited Horizon mine, north of Fort McMurray, Alberta, Canada. The materials all had low salinity and alkalinity (Table 1). TS and OB had low organic carbon contents (<0.5 g kg −1 ), and PMM was a mixture of peat and mineral soil with high porosity and organic carbon content (>5 g kg −1 ). The materials were packed according to pre-determined bulk densities (Table 1) (Figure 1). This layer was submerged into saline water (3 g L −1 of NaCl) to mimic the salinity threat from underneath in the field such that the water table depth ranged from 50 cm to 150 cm ( Table 2). Water levels in the columns were controlled at the interface between the barrier and the bottom layers using Mariotte bottles, and the amounts of saline water added to the bottles were recorded. Two more controls, which had cover soil and TS or OB barrier layer but without saline groundwater, were also included as part of the experiment (Figure 1 and Table 2). With 3 replications for each treatment, a total of 54 columns were set up.
In each column, four jack pine (Pinus banksiana) seedlings were planted at the beginning of the experiment. Water balance treatment was started in October 2013, 3 months after the seedlings were planted. Precipitation was mimicked using tap water sprinkled on the top of the columns. Average water supplies during the entire experiment were controlled at 260 and 385 mm per growing season for dry and wet conditions, respectively. During the growing season, air temperature in the greenhouse was maintained at 22 • C for 16 h of day and at 16 • C for 8 h of night to simulate the summer condition in the oil sands region. . A completely randomized block design was used in the experiment with barrier material type, barrier thickness, and water balance as treatments. It examined 2 types of barrier materials (coarse-textured tailings sand (TS) vs. fine-textured overburden (OB)), barrier thickness at 4 depths (0, 20, 50 and 100 cm), and water balance at 2 levels (simulated dry vs. simulated rainy condition for the oil sands mining region in Canada).
The two barrier materials with contrasting soil textures, TS and OB, and the cover soil material, peat-mineral soil mix (PMM), were all sourced from the Canadian Natural Resources Limited Horizon mine, north of Fort McMurray, Alberta, Canada. The materials all had low salinity and alkalinity (Table 1). TS and OB had low organic carbon contents (<0.5 g kg −1 ), and PMM was a mixture of peat and mineral soil with high porosity and organic carbon content (>5 g kg −1 ). The materials were packed according to pre-determined bulk densities (Table 1) in columns (45 cm in diameter and 180 cm in height) with the following layers ( Figure 1): (1) The top layer (cover soil): all columns were packed with 50 cm of PMM as cover soil. (2) The barrier layer: TS and OB were packed under the cover soils as barrier materials. Their thicknesses varied from 0, 20, 50, to 100 cm in different treatments. And (3) The bottom layer: The same TS was placed under the barrier materials to keep all columns to the same height (180 cm) (Figure 1). This layer was submerged into saline water (3 g L −1 of NaCl) to mimic the salinity threat from underneath in the field such that the water table depth ranged from 50 cm to 150 cm (Table 2). Water levels in the columns were controlled at the interface between the barrier and the bottom layers using Mariotte bottles, and the amounts of saline water added to the bottles were recorded. Two more controls, which had cover soil and TS or OB barrier layer but without saline groundwater, were also included as part of the experiment (Figure 1 and Table 2). With 3 replications for each treatment, a total of 54 columns were set up.
In each column, four jack pine (Pinus banksiana) seedlings were planted at the beginning of the experiment. Water balance treatment was started in October 2013, 3 months after the seedlings were planted. Precipitation was mimicked using tap water sprinkled on the top of the columns. Average water supplies during the entire experiment were controlled at 260 and 385 mm per growing season for dry and wet conditions, respectively. During the growing season, air temperature in the greenhouse was maintained at 22 °C for 16 h of day and at 16 °C for 8 h of night to simulate the summer condition in the oil sands region.

Instrumentation and Regular Monitoring
Soil dielectric permittivity (ε), bulk EC (EC b ), and temperature data in the materials were recorded hourly during the experiment by multi-functional Decagon GS3 probes. In each column, two probes were installed in the cover soil (PMM) at 20 and 40 cm depths, the other probes were installed in OB or TS, and all of them were above the groundwater level ( Figure 1). Soil volumetric water contents (VWC) were calculated from ε values using one of the following calibration equations specific for the media: In the study, EC b was first converted to the EC of pore water (EC p ), and then the EC of saturated paste (EC sp ), and total dissolved solids (TDS) were calculated using VWC, EC p , bulk density and porosity [20].

Final Sampling and Data Analysis
At the end of the third growing season, soil columns were destructively sampled. Soils were sampled by layers, and the final salt concentrations in the cover soil and barrier materials were determined and compared with the initial values in Table 1. Plant growth and root distribution under different treatments were assessed. The treatment effects on water dynamics and salt accumulation in the soil profiles and plant response after 3 growing seasons were analyzed using ANOVA and conditional inference tree (CIT) tests [21].

Model Description
Numerical modeling was conducted to expand our research duration and understand how soil salinity evolves over the longer term and tried to find out the minimum safe barrier thickness to prevent salt intrusion into the rooting zone for the PMM and TS and PMM and OB designs, respectively.
A soil water flow and solute transport model, HYDRUS-1D, which was designed to analyze water and solute movement in unsaturated, partially saturated, or fully saturated porous media [22] was used to complete this part of the work. It included the following aspects: (1) Developing soil hydraulic and solute transport parameters by interpreting/simulating water and salt dynamics in the columns; (2) evaluating the long-term performance of different capping configurations under various climatic conditions; and (3) proposing recommendations for the design of reconstructed soil profiles when using PMM and TS or PMM and OB combinations that will minimize salt migration into the rooting zone.
With the historic data reported in references as detailed above, a Marquardt-Levenberg type parameter estimation technique was applied to determine the parameters based on our experimental data [22]. Soil bulk density (BD), residual water content (Q r ), saturated water content (Q s ), α and n in the van Genuchten soil water retention function, saturated hydraulic conductivity (K s ), and solute transport parameters such as longitudinal dispersivity (Disp), and molecular diffusion coefficient in free water (Diff w ), were all estimated using the above method.
The upper and lower boundaries for water movement modeling were set up as an atmospheric boundary with surface layer and constant water content. The water content at the lower boundary equaled its saturated water content (water potential equals 0). For solute transport, constant flux concentration and variable concentration boundary conditions were used as the upper and lower solute transport conditions. Among them, irrigation, evaporation, transpiration, salt concentration in the irrigation water, and salt concentration in the recharging groundwater were temporally variable. Irrigation water in the model used the real data recorded. Salt concentration in the irrigation water was taken as 0, and salt concentration at the lower boundary was time-variable, and its values were set according to the monitored EC data from GS3 probes. Potential evapotranspiration (PET, mm) was estimated based on temperature at soil surface (T soil ) using the following equation [25]: Potential transpiration (PT, mm) was calculated as: where t is the time (day) after the seedling was planted, and t r-max is the time (day) that the root grows to its max length. Evaporation was the difference between PET and PT. Soil water dynamics were validated by VWCs at depths. Soil salt transport was validated by comparing the observed and simulated TDS in pore water at different depths in the column. Mean error (ME), mean absolute error (MAE), and root mean squared error (RMSE) were used to evaluate the fitness of the model.
where n is the number of validation points, x m is the experimental data, and x s is the simulated data [26].

Scenarios for Salinization Risk Assessment
The validated HYDRUS-1D model was run for a 15-year period under different barrier material, barrier thickness and water balance scenarios. In the simulation, barrier materials used PMM and TS or PMM and OB combinations; barrier thickness of TS or OB was set as 50, 75, 100, 150, or 200 cm; and wet, dry and normal climatic conditions were considered. The normal climatic condition in the simulation used real weather recordings from 2001 to 2015 at the weather station called Fort McMurray A (weather.gc.ca); and the wet and dry conditions were set to have 130% or 70% of the normal precipitation with other meteorological parameters the same as the normal condition.
For each scenario, modeling was performed only for the growing season (1 May to 31 October). It was assumed that there was no evaporation loss and no loss in the form of snow in winter. Snow accumulation in winter would melt in the first 15 days at the beginning of the growing season. The PET was calculated using Equation (4). The PT was set as half of PET, and the root distribution in different materials was based on the final root distribution data in the column experiment. The lower boundary conditions were the same as those in the simulation for the column experiment. Groundwater levels were maintained at the bottom of the barrier layer. Saline groundwater was the only external salt source for the system and TDS in the groundwater was 3 g L −1 .

Treatment Effects on Soil Water Dynamics
During the experiment, barrier materials or thicknesses showed apparent effects on soil water contents above the groundwater (p < 0.01, ANOVA), while the water balance treatment showed slight effects ( Figure 2 and Table 3). Cover soils were generally moister when there was no barrier layer, or the barrier thickness was thin (≤50 cm). As the barrier thickness increased, VWCs in the cover soils generally decreased regardless of the barrier material. Highest VWCs in cover soils (at 20 and 40 cm depths) were found in treatments without a capping layer; these values were significantly higher than most of the other treatments.
Water 2020, 12, 512 8 of 23 Groundwater levels were maintained at the bottom of the barrier layer. Saline groundwater was the only external salt source for the system and TDS in the groundwater was 3 g L −1 .

Treatment Effects on Soil Water Dynamics
During the experiment, barrier materials or thicknesses showed apparent effects on soil water contents above the groundwater (p < 0.01, ANOVA), while the water balance treatment showed slight effects ( Figure 2 and Table 3). Cover soils were generally moister when there was no barrier layer, or the barrier thickness was thin (≤50 cm). As the barrier thickness increased, VWCs in the cover soils generally decreased regardless of the barrier material. Highest VWCs in cover soils (at 20 and 40 cm depths) were found in treatments without a capping layer; these values were significantly higher than most of the other treatments.  Intense up-flux of saline groundwater means a high risk of rooting zone salinization, and the salts coming up with capillary water have been considered an important cause of the salinity in the rooting zone in oil sands reclamation [24,27]. Since soil moisture regimes in the columns remained almost the same at the beginning and the end of saline water application, the amount of saline water added to the Mariotte bottles generally equaled the net amount of upward capillary water. It was the gross capillary water minus deep percolation and represented the intensity of capillarity. In the experiment, the largest net up-flux of groundwater was found in treatments with thin layers of OB barrier material (20OB−, 46.25 cm yr −1 on average, and 20OB+, 42.81 cm yr −1 on average); and strong capillary effects were also found in treatments without a barrier layer (NC+ and NC−), with only 20 cm of TS (20TS− and 20TS+), and with 50 cm of OB and simulated wet condition (50OB+) (Figure 3). The CIT analysis further indicated that barrier thickness was the most significant factor affecting the intensity of capillarity (as the root node in the CIT model in Figure 4). There was on average 39.21 cm yr −1 of saline groundwater up-flux into the unsaturated zone when the barrier layer was less than 20 cm, regardless of the barrier material type and the amount of precipitation (irrigation). It was almost as much as the amount of water irrigated during the same period. When capped with 50 cm, saline groundwater up-flux was reduced to 20.14 cm yr −1 , and when the barrier layer increased to 100 cm it was further reduced to 5.57 cm yr −1 . Nevertheless, it should be noted that the capillary connections between cover layer and saline groundwater were not completely cut off even with as much as 100 cm of barrier material; especially when fine-textured OB was used, there was still significant upward capillary water movement (18.41 cm yr −1 for 100OB− and 12.61 cm yr −1 for 100OB+). Water balance treatments did not lead to statistically significant differences in capillary water up-flux, although the treatments with lower irrigation amounts generally consumed more saline water (Figure 3).

Treatment Effects on Soil Salinity
Cover soil is the main layer for root growth of the seedlings especially the following years after reclamation. Preventing the intrusion of salts from below is an important function of the barrier layer, and salts accumulated in this layer could also possibly get into the cover soil and eventually influence the success of reclamation. In the column experiment, salts accumulated quickly in some soil profiles ( Figure 5). Salinization occurred in both cover soil (PMM) and the barrier layer (TS or OB) after three growing seasons ( Figure 6). Barrier thicknesses showed strong influence on soil EC at all monitored depths, water balance didn't show significant effects on soil salinity in most cases, while the effects of barrier material were more complicated ( Figure 6 and Table 3).

Treatment Effects on Soil Salinity
Cover soil is the main layer for root growth of the seedlings especially the following years after reclamation. Preventing the intrusion of salts from below is an important function of the barrier layer, and salts accumulated in this layer could also possibly get into the cover soil and eventually influence the success of reclamation. In the column experiment, salts accumulated quickly in some soil profiles ( Figure 5). Salinization occurred in both cover soil (PMM) and the barrier layer (TS or OB) after three growing seasons ( Figure 6). Barrier thicknesses showed strong influence on soil EC at all monitored depths, water balance didn't show significant effects on soil salinity in most cases, while the effects of barrier material were more complicated ( Figure 6 and Table 3).

Treatment Effects on Soil Salinity
Cover soil is the main layer for root growth of the seedlings especially the following years after reclamation. Preventing the intrusion of salts from below is an important function of the barrier layer, and salts accumulated in this layer could also possibly get into the cover soil and eventually influence the success of reclamation. In the column experiment, salts accumulated quickly in some soil profiles ( Figure 5). Salinization occurred in both cover soil (PMM) and the barrier layer (TS or OB) after three growing seasons ( Figure 6). Barrier thicknesses showed strong influence on soil EC at all monitored depths, water balance didn't show significant effects on soil salinity in most cases, while the effects of barrier material were more complicated ( Figure 6 and Table 3). In the cover soils, salinization developed quickly after the initiation of the experiment when the cover soils were directly placed above the saline groundwater or the barrier layer was thin (20 cm). 50 cm of barrier layer had successfully prevented salt concentration increase in cover soil; however, when using OB, this thickness should be more than 50 cm to prevent salinization. Data from EC probes installed at 20 and 40 cm depths indicated that salt concentrations in the cover soils in these treatments all increased from less than 2 dS m −1 to 4 dS m −1 or more, from the beginning to the end of the experiment. These saline conditions could have stressed the plants [7,10,28]. In the barrier layers, OB tended to accumulate more salts than TS. Salinization in this layer could be largely mitigated when using 50 cm or more of TS as barrier material; however, the EC level still increased with 100 cm of OB as barrier material. When the barrier layer was not thick (≤50 cm), salt concentrations increased faster in OB than in TS at 60 and 90 cm depths, especially in its upper part (at the 60 cm depth, below the cover soil-barrier material interface, Figure 5). When the barrier thickness increased to 100 cm, although the upper part might not be strongly affected, the lower part still had a salinity problem when OB was used. Based on the experimental data, the minimum barrier thickness recommended to prevent the increase of salt concentration in the cover soil would be between 20 cm and 50 cm for TS and more than 50 cm for OB. In the cover soils, salinization developed quickly after the initiation of the experiment when the cover soils were directly placed above the saline groundwater or the barrier layer was thin (20 cm). 50 cm of barrier layer had successfully prevented salt concentration increase in cover soil; however, when using OB, this thickness should be more than 50 cm to prevent salinization. Data from EC probes installed at 20 and 40 cm depths indicated that salt concentrations in the cover soils in these treatments all increased from less than 2 dS m −1 to 4 dS m −1 or more, from the beginning to the end of the experiment. These saline conditions could have stressed the plants [7,10,28]. In the barrier layers, OB tended to accumulate more salts than TS. Salinization in this layer could be largely mitigated when using 50 cm or more of TS as barrier material; however, the EC level still increased with 100 cm of OB as barrier material. When the barrier layer was not thick (≤50 cm), salt concentrations increased faster in OB than in TS at 60 and 90 cm depths, especially in its upper part (at the 60 cm depth, below the cover soil-barrier material interface, Figure 5). When the barrier thickness increased to 100 cm, although the upper part might not be strongly affected, the lower part still had a salinity problem when OB was used. Based on the experimental data, the minimum barrier thickness recommended to prevent the increase of salt concentration in the cover soil would be between 20 cm and 50 cm for TS and more than 50 cm for OB. The CIT analysis on the final TDS in the cover soil also confirmed that barrier thickness was the most significant factor affecting salt accumulation in the cover soils (as the root nodes) (Figure 7), which was consistent with the result of CIT analysis on cumulative saline groundwater up-flux. Strong correlations were found between soil salinity and the amount of saline water up-flux especially in the upper part of the soil profiles. At the end of the experiment, the coefficients of determination between groundwater up-flux and soil ECsp at 40, 60, 90, and 140 cm depths were 0.74, Figure 6. Salinization in cover soils after 3 growing seasons. Different letters in a row indicate significant differences (p < 0.05, ANOVA) among treatments at the same depth, which were derived by comparing all the values across the NC, OB and TS trials.
The CIT analysis on the final TDS in the cover soil also confirmed that barrier thickness was the most significant factor affecting salt accumulation in the cover soils (as the root nodes) (Figure 7), which was consistent with the result of CIT analysis on cumulative saline groundwater up-flux. Strong correlations were found between soil salinity and the amount of saline water up-flux especially in the upper part of the soil profiles. At the end of the experiment, the coefficients of determination between groundwater up-flux and soil EC sp at 40, 60, 90, and 140 cm depths were 0.74, 0.58, 0.51, and 0.38, respectively. It indicates that not only the soils close to saline groundwater but also the upper cover soils could be affected by upward water and salt movement when a shallow water table is present, and the salinization process in reclaimed soils could develop very quickly and have marked effects on cover soil salinity in just several growing seasons.
Water 2020, 12, 512 13 of 23 0.58, 0.51, and 0.38, respectively. It indicates that not only the soils close to saline groundwater but also the upper cover soils could be affected by upward water and salt movement when a shallow water table is present, and the salinization process in reclaimed soils could develop very quickly and have marked effects on cover soil salinity in just several growing seasons.

Plant Response
Plant growth could be influenced by soil moisture, nutrient availability as well as salinity. In the experiment, plants showed apparent differences in both aboveground biomass and total root biomass after 3 growing seasons. Among the three treatment factors, barrier material was the most significant factor affecting plant growth ( Figure 8 and Table 3). Treatments using OB as barrier material generally had higher aboveground biomass than those using TS, and treatments without barrier layer had the lowest biomass. Water balances also showed positive effects on plant growth, except in those treatments without a barrier layer. In a period of three growing seasons, the barrier thickness exhibited slight effects on plants if the barrier layer was not missing. High aboveground biomass was found in treatments with thin barrier layer (≤50 cm), especially in the OB treatments. The experiment results also clearly showed that roots preferred to grow in OB rather than in TS, and the growth was also facilitated by water condition (Figure 8). Therefore, high root biomass amounts were found in treatments with OB barrier material and high irrigation.
Within three growing seasons, plant response was not consistent with the soil water and salinity results. Soil columns using OB material or thin barrier layer still accumulated high biomass, although they tended to accumulate more salts. The possible reason is that salt concentration in the rooting zone didn't increase high enough due to the limit of experimental time and it couldn't offset the benefit brought by soil moisture (Table 3). Therefore, the treatments should be validated over a longer period. A numerical simulation was done to further evaluate whether the soils become salinized and "unsuitable" for plants over a longer period (15 years).

Plant Response
Plant growth could be influenced by soil moisture, nutrient availability as well as salinity. In the experiment, plants showed apparent differences in both aboveground biomass and total root biomass after 3 growing seasons. Among the three treatment factors, barrier material was the most significant factor affecting plant growth ( Figure 8 and Table 3). Treatments using OB as barrier material generally had higher aboveground biomass than those using TS, and treatments without barrier layer had the lowest biomass. Water balances also showed positive effects on plant growth, except in those treatments without a barrier layer. In a period of three growing seasons, the barrier thickness exhibited slight effects on plants if the barrier layer was not missing. High aboveground biomass was found in treatments with thin barrier layer (≤50 cm), especially in the OB treatments. The experiment results also clearly showed that roots preferred to grow in OB rather than in TS, and the growth was also facilitated by water condition (Figure 8). Therefore, high root biomass amounts were found in treatments with OB barrier material and high irrigation.
Within three growing seasons, plant response was not consistent with the soil water and salinity results. Soil columns using OB material or thin barrier layer still accumulated high biomass, although they tended to accumulate more salts. The possible reason is that salt concentration in the rooting zone didn't increase high enough due to the limit of experimental time and it couldn't offset the benefit brought by soil moisture (Table 3). Therefore, the treatments should be validated over a longer period. A numerical simulation was done to further evaluate whether the soils become salinized and "unsuitable" for plants over a longer period (15 years).

Validation of the Soil Water and Salt Movement Model
Since soil hydraulic properties could change over time dramatically, especially in the first few years after reclamation [29], we used hourly experimental data more than one year after the column setup for the hydraulic and salt transport parameter inversion. The period covers days (12 September 2014-7 October 2014), and has two irrigation events during the period. Soil hydraulic and salt transport parameters used in the model are listed in Table 4. Table 4. Hydraulic and solute transport parameters 1 for the reclamation materials used in the HYDRUS-1D model. The experimental data of the NC+ treatment between October 2013 and May 2015 recorded by GS3 probes were used for validating the parameters of PMM (single layer); data of the 50TS+ and 100OB+ treatments during the same period were used for validating the parameters of PMM and TS, or PMM and OB constructed in layers, respectively. The temporal changes of measured and simulated VWC and TDS in pore water in different treatments during the whole experiment are

Validation of the Soil Water and Salt Movement Model
Since soil hydraulic properties could change over time dramatically, especially in the first few years after reclamation [29], we used hourly experimental data more than one year after the column setup for the hydraulic and salt transport parameter inversion. The period covers days (12 September 2014-7 October 2014), and has two irrigation events during the period. Soil hydraulic and salt transport parameters used in the model are listed in Table 4. The experimental data of the NC+ treatment between October 2013 and May 2015 recorded by GS3 probes were used for validating the parameters of PMM (single layer); data of the 50TS+ and 100OB+ treatments during the same period were used for validating the parameters of PMM and TS, or PMM and OB constructed in layers, respectively. The temporal changes of measured and simulated VWC and TDS in pore water in different treatments during the whole experiment are shown in Figure 9, and the validation results from 100th day after the start to the end are shown in Table 5. In the first 100 days, there were some poor agreements between the measured and simulated values because of the instability of the soil structure [29]. After that, most of the simulated values were in good agreement with the measured data. All simulated VWC values in different treatments had low ME (−0.02, 0.03), MAE (0.01, 0.03), and RMSE (0.01, 0.04), indicating that the model can be a promising tool for predicting water movement. The prediction of TDS in pore water had low ME, MAE and RMSE and was generally acceptable as well. However, the errors in treatments with TS as barrier material were relatively large, and this was probably due to the hydrophobic nature of this material [30].
Water 2020, 12, 512 15 of 23 shown in Figure 9, and the validation results from 100th day after the start to the end are shown in Table 5. In the first 100 days, there were some poor agreements between the measured and simulated values because of the instability of the soil structure [29]. After that, most of the simulated values were in good agreement with the measured data. All simulated VWC values in different treatments had low ME (−0.02, 0.03), MAE (0.01, 0.03), and RMSE (0.01, 0.04), indicating that the model can be a promising tool for predicting water movement. The prediction of TDS in pore water had low ME, MAE and RMSE and was generally acceptable as well. However, the errors in treatments with TS as barrier material were relatively large, and this was probably due to the hydrophobic nature of this material [30]. 50PMM+ 100TS+ Figure 9. Cont.
Water 2020, 12, 512 16 of 23 100OB+ Figure 9. Comparison of the experimental soil moisture and salt concentration with the simulated data using parameters in Table 2 (NC+ treatment) (• average data of the treatment recorded by GS3 probes, -simulated data from HYDRUS-1D).   Table 2 (NC+ treatment) (• average data of the treatment recorded by GS3 probes, -simulated data from HYDRUS-1D).

Modeling of Long-Term Soil Salinity Evolution
Soil water and salt dynamics in a series of capping designs which used TS or OB as barrier material from 0 to more than 2 m thick were simulated under normal, wet and dry conditions. The final salt distribution in profiles and the cumulative water fluxes at lower boundary of the soil profiles were used to identify the reliability of the design.
When PMM and TS were used as reclamation materials, soil profiles using 50 cm of TS as a barrier layer would be seriously salinized after 15 years under the dry condition (green lines in Figures 10 and 11), especially in the lower part of the PMM and the lower part of TS. Under the normal precipitation condition, it is still under saline risk at the bottom of the PMM layer. Using this capping strategy, there would be continuous net upward water flux from the bottom under the dry condition, and periodic up-flux under the normal condition. This confirmed again that design with a thin barrier layer would risk having salts migrate into the rooting zone ( Figure 11). When the barrier thickness was over 75 cm, the salinity level in the cover soil was always lower than the critical EC sp for plant growth (4 dS m −1 ). From this simulation, it can be summarized that the minimum barrier thickness to restrain salt intrusion into the rooting zone for the PMM and TS combination would be 50 cm and 75 cm.
When OB was used as a barrier material, soil capillarity was more intensive, and more salt accumulated over the 15-year simulation period. Layers with high salt concentration or EC values formed in the soil profile, and the thinner the barrier layer, the closer was the cover soil to the saline layer. When the thickness of the barrier layer (OB) was 100 cm, the salts would intrude into the cover soil, especially under dry condition. The continuous upward flux also indicated that this design would have potential risk under dry and normal conditions ( Figure 11). When the barrier thickness was increased to 150 cm, the EC sp in the cover soil might not reach a problematic level. Therefore, the recommended minimum safe barrier thickness to restrain salt intrusion into the rooting zone for the PMM and OB combination would be 50 cm and 200 cm, which was more than double of the thickness of PMM and TS combination. Soil water and salt dynamics in a series of capping designs which used TS or OB as barrier material from 0 to more than 2 m thick were simulated under normal, wet and dry conditions. The final salt distribution in profiles and the cumulative water fluxes at lower boundary of the soil profiles were used to identify the reliability of the design.
When PMM and TS were used as reclamation materials, soil profiles using 50 cm of TS as a barrier layer would be seriously salinized after 15 years under the dry condition (green lines in Figures 10 and 11), especially in the lower part of the PMM and the lower part of TS. Under the normal precipitation condition, it is still under saline risk at the bottom of the PMM layer. Using this capping strategy, there would be continuous net upward water flux from the bottom under the dry condition, and periodic up-flux under the normal condition. This confirmed again that design with a thin barrier layer would risk having salts migrate into the rooting zone ( Figure 11). When the barrier thickness was over 75 cm, the salinity level in the cover soil was always lower than the critical ECsp for plant growth (4 dS m −1 ). From this simulation, it can be summarized that the minimum barrier thickness to restrain salt intrusion into the rooting zone for the PMM and TS combination would be 50 cm and 75 cm.
When OB was used as a barrier material, soil capillarity was more intensive, and more salt accumulated over the 15-year simulation period. Layers with high salt concentration or EC values formed in the soil profile, and the thinner the barrier layer, the closer was the cover soil to the saline layer. When the thickness of the barrier layer (OB) was 100 cm, the salts would intrude into the cover soil, especially under dry condition. The continuous upward flux also indicated that this design would have potential risk under dry and normal conditions ( Figure 11). When the barrier thickness was increased to 150 cm, the ECsp in the cover soil might not reach a problematic level. Therefore, the recommended minimum safe barrier thickness to restrain salt intrusion into the rooting zone for the PMM and OB combination would be 50 cm and 200 cm, which was more than double of the thickness of PMM and TS combination.

Conclusions
After a comparison of the effects of different capping strategies and different hydrologic conditions on salt movement in oil sands reclamation soils, we conclude that barrier thickness was the most important factor that should be considered to minimize salt migration from below into the rooting zone when a shallow water table is present. Further, the risk of salinization in the rooting zone was higher with the OB barrier material than with the TS due to the greater capillarity of the OB material, although the OB material facilitated plant growth within the first a few growing seasons after reclamation. If less than 20 cm of TS or less than 50 cm of OB was used as the barrier layer, salt concentration in the cover soil could increase dramatically under the simulated saline groundwater conditions and some weather conditions. A further conclusion drawn from the 15-year modeling work was that when 50 cm of PMM was used as a cover soil and the water table is at the bottom of the barrier layer, the minimum barrier thickness to prevent salt intrusion into the rooting zone in long term was 75 cm for TS and 200 cm for OB as a barrier material.