Reclaiming Wetlands after Oil Sands Mining in Alberta, Canada: The Changing Vegetation Regime at anExperimental Wetland

Surface mining for oil sand results in the formation of large pits that must be reclaimed. Some of these pits are backfilled with a myriad of substrates, including tailings rich in cations and anions, to form a solid surface. Experimental reclamation of the East in-pit located on the Syncrude Canada Ltd. mine lease was initiated in 2011 with Sandhill Wetland. Here, we report on monitoring (between 2015 and 2021) of Sandhill Wetland plant communities and significant environmental features, including base cations and water tables. Multivariate analyses demonstrated that the three dominant plant communities established in 2013 have continued to be dominated by the same species nine years after reclamation was initiated, but with reduced species richness. Plant communities have shifted across the wetland in response to water table changes and increases in sodium concentrations. The stoichiometry of base cations is unlike the natural wetlands of the region, and the surficial water chemistry of the wetland is unique. In response to variability in precipitation events coupled with wetland design, water tables have been highly variable, creating shifting water regimes across the wetland. Plant community responses to these shifting conditions, along with increases in base cation concentrations, especially sodium, provide background data for future in-pit reclamation designs. The plant responses underscore the need to develop reclamation designs for landscapes disturbed by mining that alleviate extreme water table fluctuation events and decrease cation concentrations to levels that approach natural wetlands.


Introduction
Peatlands are an integral part of the boreal forest landscape, with approximately 24% of the boreal forest worldwide covered by peatlands [1]. Globally, peatlands store between 20-30% of the world's soil carbon [2,3] and 10% of soil nitrogen [4]. Peat-forming wetlands (bogs and fens) are an important part of the western Canadian landscape, covering approximately 365,000 km 2 or 40% of the western boreal forest of Canada [5]. In western Canada, the Oil Sands Region (OSR) occupies approximately 142,000 km 2 and has been estimated to hold 1.65 billion barrels of crude oil [6]. The Athabasca Oil Sands Area (AOSA) is the largest deposit in the region (66% of the OSR) and includes a mineable surface area for bitumen of 488,970 ha. Wetlands cover 43% of the AOSA landscape and include 28% fens and bogs and 15% mineral wetlands (swamps and marshes) [7].
Peatlands develop in specific landscape settings and respond to changes in hydrological characteristics of the surrounding landscape, as well as internal autogenic interactions [8]. Three ecological gradients have long been recognized as important in influencing vegetation patterning and peatland development: (1) pH and base cation (and associated anion) concentrations provide the basis for different bog and fen site types, each with specific floral elements [9,10]; (2) water tables and seasonal fluctuations set boundary conditions for the structural aspects of peatland site-types; and (3) ecosystem structure and function, especially decomposition and plant growth, and thus peat accumulation, are influenced significantly by nutrient inputs (N, P) that originate from the surrounding landscape and precipitation [11,12]. Changes in these parameters can affect the functioning of peatlands over short periods resulting in changes in vegetation (or structural aspects) and flora (including species diversity) of individual peatlands.
First recognized by Weber in 1902 [13], bogs are hydrologically isolated from surrounding terrestrial inputs, while fens are closely connected to and influenced by hydrological inputs from associated uplands. Changes in water qualities and quantities, especially those from the inflowing surface and groundwater, can dramatically alter carbon storage, species diversity, and nutrient retention, eventually altering the direction of plant succession and peatland site type. Increases in water table fluctuation and salinity have a considerable influence on whether individual sites will develop as marshes, brackish wetlands, saline wetlands, or riparian meadows [14].
Open pit mining causes extensive disturbance to the landscape, including clearing forests and wetlands, draining or diverting surface and groundwater, and removing tens of meters of overburden. After extracting the bitumen, there remain large, deep 'in-pits' filled with composite tailings containing high concentrations of base cations and anions originating either from natural sources or from process waters with high salinity levels. Provincial legislation mandates the reclamation of previously mined areas back to equivalent land capability following mining closure [15]. Hence, appropriate protocols for reclaiming wetlands are a focus for many oil sand mining companies. Among the many challenges that exist in these reclamation efforts are providing suitable water quality and quantities that support appropriate wetland plant species [16] and function from an integrated hydrological perspective [17,18]. To test the geotechnical feasibility of constructing a watershed capable of supporting a wetland placed on a composite tailings deposit, an experimental watershed containing a wetland (Sandhill Wetland-SHW) was constructed, with water and plants introduced in 2012 [19].
The importance of monitoring the vegetation and hydrological changes at Sandhill Wetland is that future wetland reclamation sites will be constructed above previously mined in-pits that contain saline-sodic substrates, wherein the diffusion and attenuation of sodium (~1100 mg Na L −1 ) have not been shown to lessen by overburden depth in the first four years [20,21]. During the first nine years of post-wet-up development, the vegetation and hydrology at SHW have changed. Previous research has centered on hydrological processes and has generated considerable information on early aspects of the watershed. For example, Nichols et al. [22] showed the hydrological fluctuations of SHW during two normal rainfall years, 2013 and 2014. The inflow/outflow pumping of water contributed significantly to the hydrological fluxes in 2013; however, the artificial addition of water was negligible in 2014. The evapotranspiration rates at SHW were similar between years but were lower than older nearby upland systems. After 2014, hydrological management efforts declined, and Spennato et al. [23] found that while the water tables at SHW were on average greater than those observed in natural analogues, the lowland features at SHW exhibited similar functions of lateral flow reversals to natural systems. Among the trends reported over the past seven years has been an increase in concentrations of sodium supporting the continued upward transport of underlying process-affected water [18,24].
The predicted high concentrations of sodium at Sandhill Wetland were a critical component in selecting appropriate wetland plants to form the foundational vegetation of the developing wetland. Among potential candidate plant species, several appeared to be of interest. The broad tolerances of Carex aquatilis (water sedge) to a series of natural environmental gradients [25] suggested that it was a species of interest for in-pit reclamations, and its tolerances to sodium were examined in greenhouse studies [26,27]. Additionally, tolerances to sodicity were examined for a widespread, early successional species, Beckmannia syzigachne (slough grass) [28], and for Typha latifolia (broad-leaved cattail) [29], showing de-creased morphological performance with increased sodium concentrations, notably above 850 mg L −1 and 300 mg L −1 , respectively. Early colonization from local sources by species of bryophytes was confirmed by Vitt and House (2015) [30]. After three years post-wet-up, the wetland was composed of four distinct plant assemblages with distributions mostly related to wetness [31]. Hartsock et al. [32] evaluated the porewater chemistry from decadal monitoring stations present at SHW, with comparisons to provincial water quality guidelines and the chemistry of natural wetlands. Additionally, Hartsock et al. [33] compared the vegetation and surficial water chemistry present in 2018 at SHW to that of regional natural wetlands, concluding that SHW was most similar to brackish marshes. Similarly, Biagi et al. [34] reported the lack of peat profile stratification and the suggested absence of water-conserving feedback mechanisms were likely to facilitate marsh development.
Here we examine the development of vegetation on Sandhill Wetland over a period of three to nine years post-wet-up and how vegetation change has affected species richness and plant occurrences on the wetland. We also provide summaries of salinity and water table variation changes that may have influenced the temporal changes in vegetation responses since wet-up. In particular, from our annual monitoring data obtained 3-9 years post-wetup, we examine plant abundances, base cation concentrations, and water tables at SHW to: (1) examine the development of plant communities on the wetland; (2) determine what components of these early communities changed; (3) follow the changes in spatial shifting distributions of these early plant communities; and (4) understand how plant community evolution on the wetland has been associated with variations in water chemistry and water tables. This work follows that presented in Vitt et al. [31] and Hartsock et al. [33], which assessed earlier vegetation patterns, while this more comprehensive study assesses the temporal and spatial evolution of plant communities that originated on the wetland after wet-up.  [35]. The 10 m tailings sand cap was overlain by 0.5 m of clay till and 0.5 m of fresh peat harvested from a nearby peatland slated for mining with mixed characteristics, including microsites having bog and rich fen vegetation. Peat placement was completed in January 2011. The site consists of a 52 ha watershed consisting of constructed upland hills and a central 17 ha wetland.

Water Management
Sandhill Wetland was designed to maintain wetness and mitigate salinity levels by controlling the inflow and outflow of water through a system of underdrains and pumps and landform design features [19,36]. The system allows freshwater to be supplied from a nearby undisturbed lake, while outflows can be introduced to the surrounding reclaimed areas. Outflows can only occur when the outflow pumps are activated [19,34]. The water inflow only occurred during 2013, with a small amount in 2014. During these two years, outflow volumes were similar to inflows [34]. In 2015, a dry year, outflow pumps released 55 mm of water over 60 h. In 2016, a wet year with two rain events over 50 mm on 9 July and 3 September, outflow pumps remained off, increasing water tables in the late season [34]. An unusually wet spring in 2017 with several rain events between 20-40 mm increased water tables. Throughout this period, outflow pumps remained off, creating extremely high water tables throughout the wetland. In mid-summer 2017, outflow pumps were turned on, releasing some 17,000+ m 3 of water and lowering water tables by about 102 mm [34]. In 2018, another year with several spring rain events above 20 mm and the 21 July event of 92 mm again raised water tables that were lowered some 210 mm with outflow pumping [34]. Although data on outflow pumping are unavailable for 2019-2021, high water tables were maintained over the wetland in 2019-2020 due to high precipitation totals. In 2021, water tables were much reduced to the lowest levels observed on the wetland, leaving many areas with water below the soil surface.

Early Vegetation Structure: Years 1-3
In November 2011, a native seed mix was spread across the wetland. The seed mixture contained the following species: Carex aquatilis, Carex diandra, Carex utriculata, Scirpus atrocinctus, Carex bebbii, Carex paupercula, Scirpus microcarpus, Carex lasiocarpa, Carex rostrata, Carex limosa, Carex interior, and Juncus tenuis. Eighty percent of the seed mixture was comprised of Carex aquatilis [31]. In the spring of 2012, 16 plant introduction plots were established and planted with a variety of greenhouse-grown wetland plants. In the spring of 2013, additional species were introduced to these plots. The site remained dry until August 2012, when water was introduced from a natural lake through a series of pipes and a water storage pond [19], with continued additional water inputs in 2013 [19,34]. Thus, 2013 was the first full year of adequate water available for plant growth.
In 2013, plant establishment was limited, with only scattered individuals of graminoids, occasional annual upland weeds, and cattails present on the site. By 2014, plant abundances had increased, with mean vascular plant cover in the wetland area reaching 65% [30]. Vascular plant cover increased in 2015 to between 90-100% cover, with three distinct plant assemblages zonally distributed on the wetland. The main three zones were characterized by different dominant plant species and different water tables. The wettest zone, with persistent standing water, was dominated by Typha latifolia; the second zone with water tables close to or at the peat surface for most of the growing season was dominated by Carex aquatilis, and the third zone with water tables well below the peat surface was dominated by Calamagrostis canadensis [31].

Vegetation and Plant Communities
Using a grid of 78-87 (depending on annual access) permanent plots centered at 40 m intervals across SHW, we conducted annual vegetation surveys in early August from 2015-2019 and 2021, representing years 3 through 7 and 9 following the introduction of water to the wetland. No sampling occurred in 2020. We present plant community and environmental data for four of these seven years, beginning in 2015 and continuing through 2017, 2019, and 2021. At each grid point we identified all plants (vascular plants and bryophytes) within a circular 8 m 2 plot and estimated abundances to the nearest 5% using the same personnel and following the identical protocols in Vitt et al. [31].
Plant species abundances from each plot were utilized in a Bray-Curtis nonmetric multidimensional scaling (NMDS) ordination with 25 restarts using field-derived abundance data. Group average cluster analyses to delineate plant communities were performed each year. Similarities for the four major clusters ranged from 30-45% (Supplementary Materials Figure S1). PERMANOVA was used to test whether a priori groups (e.g., plant communities) differed among years [37]. A two-way PERMANOVA was used, testing groups (plant communities) determined by the cluster analyses and year. Both were fixed effects. Pair-wise tests were then used to determine which groups and years differed. Abundances for the three most abundant plant species were overlain as segmented bubbles on the yearly ordinations.
Species richness was evaluated through four measures: (1) alpha diversity is the number of species recorded in individual plots; (2) gamma diversity is the total number of species in a group; (3) beta diversity is a measure of the amount of variation in species composition or species turnover between plots calculated as gamma divided by alpha [38]; and (4) additionally, we used Pielou's J as a measure of evenness [39]. Since bryophytes occurred only in low abundances, the frequencies of bryophyte species were estimated by summing the number of occurrences of all species in the sampled plots.
Linear regression was used to analyze the abundance of dominant species (Calamagrostis canadensis, Carex aquatilis, and Typha latifolia) across the entire wetland [40].

Water Chemistry
In 2011, 10 water samplers that allowed wetland water to be collected in 10 cm intervals to 50 cm depth were installed in the wetland area. In 2015, we added 16 additional samplers across the site to better capture spatial water chemistry patterns. Each sampler consisted of 5 sections of thinly slotted 2.5 cm-diameter PVC pipe. We report here on water collected in early August of each year from the surficial 10 cm segment. In 2019, we sampled nearsurface water from each of the gridded vegetation plots, either from standing surface water or from shallow pits dug in the peat. Samples were filtered immediately through Whatman 541 filter paper into acid-washed Nalgene bottles and refrigerated. In the laboratory, water samples were analyzed for pH using an Accumet AB15 pH meter, electrical conductance (EC) using an Orion 4 Star EC meter, and Na + , Ca 2+ , and Mg 2+ using a Varian 220 FS atomic absorption spectrophotometer.

Water Tables
Water tables vary between years [34] and over the growing season, with seasonal highs in May and early June due to snowmelt and influxes from surrounding uplands. Water tables decrease over the growing season and stabilize in July and August, periodically recharged from summer storm events [34]. We measured water tables during vegetation surveys in 2016-2019 in early August, during a period without precipitation during the previous week. Our data reflect a growing season relative water table for mid-summer each year. Water tables at each gridded plot were measured relative to the soil surface with a 30 cm ruler. We placed the measurements into five categories representing: deep water (40.1-100.0 cm) above the soil surface, shallow water (10.1-40.0 cm) above the soil surface, water just above the soil surface (0.1-10.0 cm), moist soil with water near to just beneath the soil surface (0.0-−9.9 cm), and dry soil with water −10-−50 cm beneath the soil surface.

Spatial Analyses
Distributions of plant communities and dominant species, sodium, and water tables were mapped using data from the gridded plots and/or water sampler locations. Georeferenced points were created from coordinates recorded at the surveyed locations using a handheld Garmin Montana 680 GPS unit. Maps were generated in ArcMap version 10.6.

Plant Communities on Sandhill Wetland
Although we sampled vegetation at Sandhill Wetland annually between 2015-2019 and in 2021, we present data for alternate years-2015, 2017, 2019, and 2021-as representative across the seven years. Across all years, four plant communities are different from one another (Group average cluster analysis at 30-45% similarity (Supplementary Materials Figure S1) and PERMANOVA p = 0.001). The four communities also differ between years (PERMANOVA p = 0.001), and there is a significant interaction between plant community and year (PERMANOVA p = 0.001). The four communities ordinate along the x-axis correlated to wetness and base cation gradients (Figure 1a

Plant Communities on Sandhill Wetland
Although we sampled vegetation at Sandhill Wetland annually between 20 and in 2021, we present data for alternate years-2015, 2017, 2019, and 2021-a sentative across the seven years. Across all years, four plant communities are from one another (Group average cluster analysis at 30-45% similarity (Supple Materials Figure S1) and PERMANOVA p = 0.001). The four communities also d tween years (PERMANOVA p = 0.001), and there is a significant interaction betwe community and year (PERMANOVA p = 0.001). The four communities ordinate a x-axis correlated to wetness and base cation gradients (Figure 1a   . DTW = distance of the water table from soil surface; water tables increase above the soil surface to the left, water tables increase below the soil surface to the right; Ca 2+ , Mg 2+ , and Na + concentrations increase to the lower right. Plant community 1 consists of plots having high abundances of Typha latifolia, with Carex aquatilis and Utricularia minor occurring in low abundances. Overall, these plots are species-poor. Plant Community 2 consists of plots with high abundances of Carex aquatilis. Other abundant species in this plant community include Carex utriculata (in ±50% of the plots), Calamagrostis canadensis (in ±25% of the plots), Salix petiolaris, and Typha latifolia (in 40-45% of the plots in low abundance).
Plant Community 3 has high abundances of Calamagrostis canadensis. Other important species in this community include Carex aquatilis (occurring in low abundances in about two-thirds of the plots), Carex atherodes (occurring in 10-15% of the plots), Carex utriculata (occurring in one-third of the plots), Populus balsamifera, Salix monticola, and Salix petiolaris.
In 2015, five plots were dominated by other species. One plot was characterized by the presence of Carex atherodes, one plot was characterized by Carex utriculata, and three plots were characterized by a suite of species found in dry areas of the wetland: Poa pratensis, Epilobium angustifolium, Melilotus alba, and Sonchus arvensis (Figure 1a). In 2017, these outlier plots comprised three plots. Two of these plots were characterized by the presence of Carex atherodes (Figure 1b-on the lower right of the ordination), and the third plot had low species cover overall. In 2019, three plots were characterized by the presence of Carex atherodes (Figure 1c). In 2021, outliers consisted of five plots: Two plots characterized by Carex atherodes, two plots characterized by Carex utriculata, and the fifth characterized by the presence of both Carex atherodes and Carex utriculata ( Figure 1d).
Overall, the most abundant plants on SHW are Typha latifolia, Carex aquatilis, and Calamagrostis canadensis. We plotted the abundances of these dominant species as recorded in each of the four years of vegetation surveys on the yearly ordinations (Figure 1a-d). These three species consistently remained dominant in the wetland throughout the study period, with the three plant communities remaining distinct. However, the spatial distributions of the plant communities and abundances of dominant species in individual plots changed over the years.

Spatial Changes of Plant Communities
Plots in Plant Community 1, with high abundance of Typha latifolia, increased from 11 plots in 2015 to 26 in 2021 (Table 1), with the greatest increases in 2017 (11 to 17) and 2021 (16 to 26) ( Table 1). These plots occupy the central portions of the wetland, and over the seven years have expanded outward into areas dominated by Carex aquatilis (Figure 2a-d). Over the seven years of sampling, overall abundance of T. latfolia on the wetland significantly increased from 6% cover in 2015 to 26% cover in 2021 (p = 0.018, r 2 = 0.791). The number of plots in Plant Community 2, with high abundance of Carex aquatilis, has remained relatively stable over the seven years (Table 1); however, in the central portions of the wetland, individual plots dominated by C. aquatilis in 2015 have been replaced by those dominated by Typha latifolia, while plots located in drier areas dominated by Calamagrostis canadensis in 2015 have been replaced by those dominated by C. aquatilis in 2021 (Figure 2e-h). Overall, the abundance of C. aquatilis remained the same over the seven years (p = 0.693, r 2 = 0.043).
In 2015, Plant Community 3, with a high abundance of Calamagrostis canadensis, occurred in just over one-third of the plots (37 total), including much of the marginal areas on the northern and southern sides of the wetland. By 2021, the number of plots with a high abundance of C. canadensis decreased to 22 plots, or about one-fourth of the plots (Table 1). Many C. canadensis-dominated plots transitioned into Carex aquatilis-dominated plots in later years. Especially noteworthy is the loss of C. canadensis-dominated plots along the northern boundary of the wetland (Figure 2i-l); however, the overall abundance of C. canadensis on the wetland did not significantly decrease (p = 0.063, r 2 = 0.621). The number of plots in Plant Community 2, with high abundance of Carex aquatilis, has remained relatively stable over the seven years (Table 1); however, in the central portions of the wetland, individual plots dominated by C. aquatilis in 2015 have been replaced by those dominated by Typha latifolia, while plots located in drier areas dominated by Calamagrostis canadensis in 2015 have been replaced by those dominated by C. aquatilis in 2021 (Figure 2e-h). Overall, the abundance of C. aquatilis remained the same over the seven years (p = 0.693, r 2 = 0.043).
In 2015, Plant Community 3, with a high abundance of Calamagrostis canadensis, occurred in just over one-third of the plots (37 total), including much of the marginal areas on the northern and southern sides of the wetland. By 2021, the number of plots with a high abundance of C. canadensis decreased to 22 plots, or about one-fourth of the plots (Table 1). Many C. canadensis-dominated plots transitioned into Carex aquatilis-dominated plots in later years. Especially noteworthy is the loss of C. canadensis-dominated plots

Bryophyte Species Richness
Over the entire wetland, bryophyte species richness declined over the seven years, with the number of occurrences decreasing from 316 in 2015 to 74 in 2021 ( Table 2, Supplementary  Materials Table S1). In 2015, 64 (74%) of the 87 sampled plots had at least one species, and 20% of the plots had more than five species. Comparatively, in 2021, 18% of the plots had at least one species, and only 3% of the plots had more than five species. Along with the decrease in the number of occurrences, alpha diversity decreased from 4.5 to 3.3 species, and gamma diversity decreased from 24 to 13 bryophyte species found on the wetland. In 2015, bryophyte gamma diversity varied from zero in T. latifolia-dominated plots to 18 species in C. aquatilis-dominated plots and 19 in C. canadensis-dominated plots. This diversity declined over the seven years by 50% to 9 species in C. aquatilis plots and by 32% to 13 species in C. canadensis plots. The greatest decrease in the number of species was in 2017 ( Table 2). Bryophyte beta diversity (Table 1) could only be calculated across the years for C. canadensis-dominated plots. Here, the bryophyte beta diversity gradually increased from 3.9 in 2015 to 6.5 in 2021, indicating a gradual increase in species turnover between plots. In 2014, only two years after the wet-up, a diverse set of bryophyte species was present in the wetland. Sixty-six of the 87 sampled plots had at least one species, and 18% of the plots had more than five species [30]. This high diversity and frequent occurrence of species continued through 2015-2016. Beginning in 2017, the number of plots with no bryophytes increased from 29 to 58, with species occurrence decreasing by 30%; however, overall species richness at both the plot and site levels remained steady. The number of plots with high diversity also remained similar to past years. In 2018, site species richness declined from 25 to 17 species, occurrences from 182 to 124, and high diversity plots from 19 to 8. This decline continued through 2021, and by 2021, 71 (or 82%) of the plots contained no bryophytes. During the first 3-4 years, the developing bryophyte ground layer contained a highly diverse set of fen species. In 2017 and 2018 and associated with high water tables that created more areas of standing water, along with higher sodicity, the number of bryophyte occurrences decreased substantially; however, site diversity was less affected, suggesting that there are a low number of local habitats on the wetland that have remained suitable for bryophytes and may continue to support a healthy, but limited bryophyte flora.

Vascular Plant Species Richness
In 2015, 96 species of vascular plants occurred on the wetland, declining to 44 in 2021 (Supplementary Materials Table S2). However, this marked decline was not uniform across the three plant communities. In T. latifolia-dominated areas, species richness (gamma diversity) declined from nine species to five in 2021. Comparatively, in C. aquatilis-dominated areas, gamma diversity declined from 68 species in 2015 to only 17 species in 2021 (a 75% decrease in the number of species). In C. canadensis-dominated areas, richness declined from 73 species in 2015 to 38 species in 2021 (a decrease of 48%). Alpha diversity followed similar trends. Overall, species richness declined 45-46% in T. latifoliaand C. canadensis-dominated plots and by 75% in C. aquatilis-dominated plots. This large decrease is especially evident in 2017 in C. aquatilis-dominated plots when richness decreased from 86 species to 37 (60%) ( Table 1).
Beta diversity for vascular plants, a measure of species turnover between plots, varied across years within communities but was consistently different between communities, with the least species turnover evident in T. latifolia-dominated plots and the most turnover in C. aquatilis-dominated plots (Table 1).
Pielou's measure of species evenness declined in two of the three plant communities. In T. latifolia-dominated plots, evenness declined from 0.46 in 2015 to 0.08 in 2021. Evenness also declined (but less so) in C. aqutilis-dominated plots declining from 0. 45 (2015) to 0.37 (2021), while C. canadensis-dominated plots remained relatively stable, varying from 0.50 in 2021 to 0.65 in 2019 (Table 1).

Stoichiometric Comparisons of Base Cations at Sandhill Wetland to Natural Sites
Natural occurring rich fens have surficial (0-10 cm depth) water chemistry dominated by divalent cations (Ca 2+ and Mg 2+ ), whereas brackish wetlands (marshes on shallow peat and fens on deep peat) have waters dominated by Na + (Figure 3). In general, divalent cations at natural rich fen sites are present in concentrations less than 2 mmol L −1 , while brackish sites have divalent concentrations up to 4 mmol L −1 . Comparably, natural rich fens have Na + concentrations mostly below 1 mmol L −1 but up to 2 mmol L −1 compared to brackish wetlands with higher concentrations of Na + (4-10 mmol L −1 ) (Figure 3). Wetland sites in Alberta traditionally classified as saline wetlands (or saline fens or saline marshes) have surficial waters dominated by Na + (>20 mmol L −1 ) and are better recognized as sodic wetlands. These cation-rich wetlands also have somewhat higher concentrations of divalent cations (<20 mmol L −1 ) but with Na + concentrations about 10 times those of the divalent cations ( Figure 3).  Like sodic wetlands and unlike natural rich fens, the surficial waters at Sandhill land have high concentrations of Na + . Still, the stoichiometric ratio to divalent cati quite different (calculated on a mass basis-mmol L −1 ). At SHW, the ratios of Na + to + Mg 2+ range from 1.7 (2017-2019) to 3.7:1 (2021) ( Table 3), with most sites having centered on a ratio of 3:1 in 2021, compared to ratios between 6-10:1 for natural wetlands (Figure 4). When SHW waters are compared to natural rich fens and bra wetlands, they differ in higher divalent cation concentrations (5-20 mmol L −1 com Like sodic wetlands and unlike natural rich fens, the surficial waters at Sandhill Wetland have high concentrations of Na + . Still, the stoichiometric ratio to divalent cations is quite different (calculated on a mass basis-mmol L −1 ). At SHW, the ratios of Na + to Ca 2+ + Mg 2+ range from 1.7 (2017-2019) to 3.7:1 (2021) ( Table 3), with most sites having ratios centered on a ratio of 3:1 in 2021, compared to ratios between 6-10:1 for natural sodic wetlands (Figure 4). When SHW waters are compared to natural rich fens and brackish wetlands, they differ in higher divalent cation concentrations (5-20 mmol L −1 compared to <2 mmol L −1 for natural rich fens) (Figure 4).

Variation of Base Cations and pH of Surficial Waters
The pH of surface water at Sandhill Wetland has varied little over the years, with means varying from 7.9 to 8.1. In contrast, surface water Na + concentrations have increased over the years (Table 3). In 2015, Na + concentrations ranged from 89 to 476 mg L −1 (n = 14, mean: 379 mg L −1 ). Comparatively, in 2021, Na + ranged from 64 to 1950 mg L −1 (n = 83, mean: 705 mg L −1 ). During these seven years, Na + concentrations at sites in the central portion of the wetland have remained at the lowest levels, while sites in other parts of the wetland have increased concentrations, expanding the variation present on the wetland (standard deviation of 267 in 2015 and 319 in 2021). In all four years of data presented here (Figure 5a-d (Table 3), a 1.2 fold increase in mean values. In 2015 on a mass basis, the ratio of mean Na + to divalent cations was 2.6 compared to 3.7 in 2021. On a charge basis (meq L −1 ), Na + has increased in proportion to divalent cations, representing 57% in 2015, compared to 66% in 2021.

Variation in Water Tables
Water tables in early August varied from 92.5 cm above the soil surface to 50 cm below the surface (Figure 6a-c). In 2016 (our first year of water table data), 42% of the sampled plots had water tables 10 cm or more below the soil surface (dry plots), while 34% of the plots had water 10 cm or more above the soil surface (wet plots). Twenty-four  (Table 3), a 1.2 fold increase in mean values. In 2015 on a mass basis, the ratio of mean Na + to divalent cations was 2.6 compared to 3.7 in 2021. On a charge basis (meq L −1 ), Na + has increased in proportion to divalent cations, representing 57% in 2015, compared to 66% in 2021.

Variation in Water Tables
Water tables in early August varied from 92.5 cm above the soil surface to 50 cm below the surface (Figure 6a-c). In 2016 (our first year of water table data), 42% of the sampled plots had water tables 10 cm or more below the soil surface (dry plots), while 34% of the plots had water 10 cm or more above the soil surface (wet plots). Twenty-four percent of the plots had water near the soil surface (10 cm above to 10 cm below the surface-moist plots). Increases in 2017 water tables are quite clear, whereby dry plots decreased to 10%, wet plots increased to 68%, while moist water tables remained similar (22%). Both 2018 and 2019 were wet years with relatively high precipitation totals. In 2018, outflow pumps decreased water tables to some degree. Dry plots increased to 23%, wet plots decreased to 39%, and moist plots increased to 38%. Overall, sites with water tables above the soil surface increased from 46% in 2016 to 82% in 2017 and 71% in 2019. In 2017, and continuing into 2018-2019, the increase in wet habitats favorable to T. latifolia-dominated plant communities allowed this plant community to expand, while the decrease in dry habits favorable to the C. canadensis-dominated community allowed the C. aquatilis-dominated plant community to expand, albeit with attendant losses to the expanding T. latifolia-dominated community. surface-moist plots). Increases in 2017 water tables are quite clear, whereby dry plots decreased to 10%, wet plots increased to 68%, while moist water tables remained similar (22%). Both 2018 and 2019 were wet years with relatively high precipitation totals. In 2018, outflow pumps decreased water tables to some degree. Dry plots increased to 23%, wet plots decreased to 39%, and moist plots increased to 38%. Overall, sites with water tables above the soil surface increased from 46% in 2016 to 82% in 2017 and 71% in 2019. In 2017, and continuing into 2018-2019, the increase in wet habitats favorable to T. latifolia-dominated plant communities allowed this plant community to expand, while the decrease in dry habits favorable to the C. canadensis-dominated community allowed the C. aquatilisdominated plant community to expand, albeit with attendant losses to the expanding T. latifolia-dominated community.

Changing Plant Communities
Plant communities 1, 2, and 3 dominated by Typha latifolia, Carex aquatilis, and Calamagrostis canadensis, respectively, remained distinct in ordination space over the seven years of assessment. A few plots dominated by additional species also persisted in the wetland. The number of plots where Carex atherodes (a potential valuable species for saline habitat reclamation [27]) increased over time. Each of the three main plant communities continued to be dominated by the same species; however, many individual plots exhibited plant composition changes. For example, in 2015, 11 plots located in plant community 1 (dominated by T. latifolia) increased to 26 in 2021. Likewise, in 2015, plant community 3 (dominated by C. canadensis) had 36 plots which decreased to 20 in 2021. Spatial changes between 2015 and 2021 were: (1) C. canadensis-dominated areas decreased, becoming restricted to drier habitats at the margins of the wetland, (2) C. aquatilis-dominated areas overall remained steady, but migrated toward the drier margins-losing plots to the wetter interior and gaining plots towards the site margins, and (3) T. latifolia increased in areas with higher water depths, replacing plots with C. aquatilis. These community changes were coupled with a large reduction in species richness and a general loss of infrequent species on the wetland. In particular, alpha diversity in the Carex aquatilis-dominated plots decreased from 14.6 to 3.6 species at the plot level and from 86 to 26 at a whole site level, a loss of over two-thirds of the wetland flora. Alpha diversity declines were also evident in T. latifolia-dominated sites (9 species to 5). In the C. canadensis-dominated areas, total site richness decreased from 92 to 51 species. Additionally, between 2015 and 2021, there was a loss of between 39-48% of the vascular plant flora and 30-50% of the bryophytes from the wetland. These species extinctions have increased species turnover between plots and decreased evenness. Current levels of species richness at Sandhill Wetland are well below those reported for fens in Alberta [41][42][43].
Most concerning was the dramatic loss of bryophytes, with decreases at the plot level of 3.3 and 4.9 in 2015 to 0.7 and 2.0 in 2021 for C. aquatilisand C. canadensis-dominated communities, respectively. Bryophytes are an important part of peatland ecosystems, and both bogs and fens have ground layers with 90-100% cover of bryophytes [44]. In a review of primary production for wetlands of western Canada, Campbell et al. [45] estimated that for fens and bogs, the bryophyte ground layer provides approximately 41% of the total annual plant production. Comparatively, bryophytes are absent or rare in wetlands with shallow open water, marshes, and saline wetlands, with essentially no contributions to primary production [44,46]. Additionally, the developing peat column of fens and bogs is composed of un-decomposed fragments of bryophytes along with vascular plant litter and roots. Analysis of macrofossils from a 700 cm-long peat core from a rich fen in northeastern Alberta showed bryophytes composed between 40 and 80% of the identifiable macrofossils [46]. In comparison, marshes, swamps, and occasional sedge-dominated rich fens have peat deposits composed of roots and litter, with higher bulk density and high amounts of debris [47].

Changing Water Chemistry
Until 2018, porewaters at Sandhill Wetland had higher concentrations of divalent cations compared to Na + with, in general, ratios of Na + 0.3-0.8 to 1 Ca 2+ + Mg 2+ (data for 2018 not shown). After 2018, Na + concentrations have been somewhat higher (1.1-1.2:1.0) than those of Ca 2+ + Mg 2+ . These changing stoichiometric base cation ratios may play an important role in limiting species occurrences on the wetland.
Comparing cation surficial water concentrations at SHW to those from brackish and saline (sodic) wetlands and natural rich fens reveals that SHW has unique water chemistry. In the early years, porewaters at SHW had higher concentrations of divalent cations compared to Na + [24]. In 2015 and onward, Na + concentrations have increased at a greater rate than those of divalent cations, increasing the sodicity of the wetland. Natural sodic and brackish wetlands are dominated by Na + with comparatively low concentrations of diva-lent cations, whereas natural rich fen waters are characteristically dominated by Ca 2+ and Mg 2+ [48], but at lower concentrations compared to brackish and sodic wetlands [33,49]. In comparison, Na + concentrations at SHW currently are lower than most sodic wetlands, yet divalent cation concentrations are higher than brackish and natural rich fens. The stoichiometric ratios are quite different from sodic, brackish, and natural rich fens. Therefore, on a charge basis (meq L −1 ), divalent ions have a similar effect on salinity (and expressed as electrical conductance) as do those of Na + . Greenhouse studies have demonstrated that Carex aquatilis [26,27,50], C. atherodes [27], and Typha latifolia [29,51] may have reduced performance at Na + concentrations currently present in some areas of SHW; however, no one to our knowledge has examined wetland plant responses to high Na + in the presence of high divalent cations. Tables   Water table maintenance across Sandhill Wetland has been challenging, with water tables controlled by a balance of precipitation and evaporation coupled with variations in annual precipitation, along with sporadic management of removing water during high water events by outflow pumps [23,34]. Water tables at our monitoring sites were extremely variable, especially in 2016, when 54% of the sites had water tables below the peat surface and 46% above the surface. A late-season rain event in 2016, combined with a wet spring in 2017 and a delay in outflow pump activation in June 2017, resulted in unusually high, thought to be temporary, water tables throughout the wetland. These conditions resulted in 82% of our sites with water tables above the soil surface during the 2017 growing season. Pumping excess water reduced the water tables in late 2017. In 2018, a year with relatively large amounts of precipitation [34], outflow pumping reduced water tables somewhat, but in 2019, 71% of the sites remained under water. Many bryophytes and vascular plants sensitive to sodium and prolonged inundation consequently were lost from the site from 2017 to 2019.

Linkages between Environmental Factors and Plant Community Change
We contemplated several correlative environmental changes that could affect plant community responses. Firstly, there were increasing Na + concentrations in many portions of the wetland. Secondly, variable water tables may have affected vegetation, especially the 2017 water table increase. Thirdly, excess water table variation and Na + increases occurring together exposed the plant communities to a highly complex set of evolving conditions. With only a few exceptions, bryophytes are non-halophytes [52] and cannot occupy habitats with high Na + concentrations. Sodic wetlands in Alberta have few or no bryophyte species [33]. Increases in Na + on SHW may have contributed to bryophyte species extinctions in the wetland except for a few species with tolerance to Na + [53]. Increases in Na + on SHW may have contributed to bryophyte species extinctions in the wetland. Another factor potentially contributing to a loss of bryophytes is areas inundated by water, habitats similar to those found in marshes, where few bryophytes can survive. The temporary flooding and subsequent exposure of bryophytes to sodic water in 2017 may have contributed to these extinctions. Under controlled phytotron conditions, Carex aquatilis is tolerant of Na + concentrations up to 1650 mg L −1 ; however, this species cannot compete with Typha latifolia in deep water habitats, no matter the concentration of Na + . As a result, C. aquatilis was eliminated from areas with especially deep water in 2017, and those patterns continued into 2019. Furthermore, these wet conditions flooded many areas occupied by C. canadensis before the 2017 high water event, allowing C. aquatilis to expand towards the wetland margins. The loss of species richness in the C. aquatilis-dominated communities may also have been affected by the increased productivity of C. aquatilis that produced copious amounts of litter, effectively smothering smaller, less abundant vascular plants and bryophytes.

Conclusions
Among the unforeseen developments at Sandhill Wetland has been the extreme variability of water chemistry and water tables. Associated with this hydrological and chemical variability is a suite of originally diverse plant communities that, over the seven years of our monitoring, have also undergone considerable change. The plant communities have changed by the loss of species richness, namely for ground layer bryophytes, plants critical for future fen development. Cationic water chemistry has changed, with Na + increasing at marginal sites more than at centrally located sites. Ca 2+ and Mg 2+ have continued to be co-dominant cations with Na + and contribute to the site salinity. Water tables have varied spatially from more than 70 cm above the peat surface to more than 24 cm below the surface. Due to an early season flood event in 2017, water tables remained above the peat surface through 2019 across much of the site. Plant community changes associated with water table changes and increases in Na + have led to a reduction in species richness by more than 50%.
We conclude that fluctuating water tables were the main contributing factor to the movement of dominant species within the three plant communities. In contrast, excess Na + and high-water tables provided conditions unfavorable for bryophytes and smaller vascular plants. Habitat loss under highly productive C. aquatilis in areas with a high abundance of this species also contributed. Although C. aquatilis was replaced in areas of standing water by T. latifolia, areas, where water tables increased to just above the soil surface favored C. aquatilis, not C. canadensis.
Sandhill Wetland has provided the first test case of whether a peat-forming wetland could be established on an in-pit filled with composite tailings. Much has been learned over the nine years of development. These learnings have included acknowledgment that water table and salinity management are the two critical environmental factors to be reckoned with. Initial dry conditions for the first growing season (2012) provided habitats for species such as Calamagrostis canadensis, a species that has remained dominant in the drier portions of the wetland, with a reduced distribution due to water table change events. Early on, extreme wet conditions and ponding in the central portion of the wetland allowed the invasion of Typha latifolia, and it too has continued to survive and multiply in the inundated parts of the wetland. In 2011, seeds of common plant species from nearby wetlands were introduced to the wetland, and several experimental plots populated with seedlings of a diversity of wetland species were initiated [31]. The seed mix was dominated by Carex aquatilis, a species at the time considered a prime candidate for early reclamation success. This species now dominates at sites with water tables just above the soil surface, and in 2021, it was abundant and highly productive, producing copious biomass with no effect of increased salt concentrations (Vitt, pers. com.). The combined high concentrations of divalent cations, plus the contribution from Na + , SHW, is unlike any natural wetland site presently known in the province. The high concentrations of divalent cations may alleviate some of the negative effects of high Na + on plant species' health, and plant responses to the increasing salinity have been minimal. The calculation of SAR values may be useful in future monitoring.

Implications
After nine years, there have been consequences at Sandhill Wetland that provide guidance for future reclamation efforts. Among these consequences are: (1) The combination of high Na + with high divalent cations provides unique water chemistry that the dominant species in the wetland appear to tolerate; however, the plant responses to extremely high concentrations of divalent cations have not been explored. (2) The three plant communities [31] closely associated with the water table position appear to have stable dominant plant species and provide the basis for vegetation that most closely resembles that of lakeshores, sedge marshes, and wet riparian meadows, depending on the local water tables. (3) Carex aquatilis has been successfully established on the wetland and may provide suitable biomass production in the future; however, its aggressive colonization behavior may impede the development of plant communities with high species richness. (4) Although individual species have continued to dominate plant communities, the local distributions of the plant communities have been responsive to water table changes, with changing distributions and reduction of species richness. Bryophytes were almost certainly affected by sites being flooded with Na + -rich water. In 2021, the species richness had not recovered from the high water tables of 2017. (5) The periodic high water, in combination with high Na + , has all but eliminated the ground layer (bryophytes) from the wetland and greatly limited active peat accumulation processes and the future development of fen vegetation.
Future wetland reclamation must limit excessive water table fluctuations, especially those related to seasonal flooding. In hindsight, the engineered outlet that allows spring freshet flooding has been difficult to manage. Future engineering designs should allow a more natural outflow design. Proposals concerning site planning could consider those proposed by Biagi et al., 2021 [34] to better control water tables and chemistry. The unique water chemistry of the wetland also must be examined more closely to determine how overall high cationic (and associated anionic species such as sulfate and chloride) concentrations can be constrained and how these affect plant community development. Without constraints on water table fluctuations and chemistry, reclamation of fens on in-pit surfaces may be difficult; however, other wetland site types may provide hope for placing wetlands on the landscape.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.