Opportunistic wetland formation, characterization, and quantification on landforms reclaimed to upland ecosites in the Athabasca Oil Sands Region

Surface mine operators in the Athabasca Oil Sands Region (AOSR) of northeastern Alberta are required by regulation to mitigate habitat impacts resulting from their operations, including impacts to wetlands. To date, most land reclamation efforts have focused on recreating upland forestlands that resemble the surrounding natural (dry) boreal forest. However, the surficial conditions on these reclaimed upland sites can also promote spontaneous wetland development. At Suncor’s Base Plant mine, opportunistic wetlands occurring on reclamation sites have not been formally included in the current inventory of reclaimed wetland areas and remain largely unquantified. We characterized and delineated an estimated 210 ha of opportunistic wetlands (consisting of shallow open water, marshes, and swamps) using aerial photo interpretation and remote sensing analysis in combination with follow-up field verifications. The remote-based (desktop) delineations consistently underestimated actual wetland extents, due mainly to underestimations in the extent of non-inundated vegetation zones (e.g., wet meadow) as well as shrubby swamp. After field corrections, opportunistic wetland habitat was estimated to constitute ~ 17% of the total study area (1209 ha), representing more than a fourfold increase in aerial wetland extent associated with reclaimed landforms over that delineated prior to this study. The interspersion of opportunistic wetlands with upland reclaimed landforms, although unintended, more closely reflects the pre-disturbance landscape, which was characterized by a matrix of forestlands, peatlands, and mineral wetlands (in contrast to the more peatland-dominated lowlands). At Suncor, wetland vegetation composition varied significantly across the study area and was influenced by topographic variation (e.g., in elevation and % slope) in combination with the reclamation substrates (soils) that were placed prior to seeding/planting. Thus, the inclusion of opportunistic wetland delineation in reclamation tracking and closure planning merits consideration as does the opportunity to manipulate current reclamation practices to promote the establishment and persistence of wetlands on reclaimed landforms.

Abstract Surface mine operators in the Athabasca Oil Sands Region (AOSR) of northeastern Alberta are required by regulation to mitigate habitat impacts resulting from their operations, including impacts to wetlands. To date, most land reclamation efforts have focused on recreating upland forestlands that resemble the surrounding natural (dry) boreal forest. However, the surficial conditions on these reclaimed upland sites can also promote spontaneous wetland development. At Suncor's Base Plant mine, opportunistic wetlands occurring on reclamation sites have not been formally included in the current inventory of reclaimed wetland areas and remain largely unquantified. We characterized and delineated an estimated 210 ha of opportunistic wetlands (consisting of shallow open water, marshes, and swamps) using aerial photo interpretation and remote sensing analysis in combination with follow-up field verifications. The remote-based (desktop) delineations consistently underestimated actual wetland extents, due mainly to underestimations in the extent of non-inundated vegetation zones (e.g., wet meadow) as well as shrubby swamp. After field corrections, opportunistic wetland habitat was estimated to constitute * 17% of the total study area (1209 ha), representing more than a fourfold increase in aerial wetland extent associated with reclaimed landforms over that delineated prior to this study. The interspersion of opportunistic wetlands with upland reclaimed landforms, although unintended, more closely reflects the pre-disturbance landscape, which was characterized by a matrix of forestlands, peatlands, and mineral wetlands (in contrast to the more peatland-dominated lowlands). At Suncor, wetland vegetation composition varied significantly across the study area and was influenced by topographic variation (e.g., in elevation and % slope) in combination with the reclamation substrates (soils) that were placed prior to seeding/planting. Thus, the inclusion of opportunistic wetland delineation in reclamation tracking and closure planning merits consideration as does the opportunity to manipulate current reclama-

Introduction
The development of the Athabasca Oil Sands Region (AOSR) in northeastern Alberta brings large-scale and long-term anthropogenic disturbances that alter the mosaic of peatland-forest-mineral wetland complexes common in the Boreal Plains (Vitt and Chee 1989;Nwaishi et al. 2015a, b;Volik et al. 2020). Between 1980 and 2017, the total land area impacted by surface mining and in-situ activities in the Lower Athabasca region of Alberta increased from 290 to 895 km 2 (Government of Alberta 2017a). Within this forest-wetland mosaic, peatland wetlands (e.g., wooded and shrubby fens) are dominant upon the landscape, comprising up to 62% of the surface minable oil sands area, with mineral wetlands (e.g., shallow open water, marshes) constituting another 3% of the surface minable area (Raine et al. 2002). Consequently, oil sands development has disturbed the ecology and hydrology of many wetlands Daly et al. 2012;Nwaishi et al. 2015a, b). Efforts to mitigate impacts to affected forest and wetland habitats are underway as mandated by government Environmental and Protection Act approvals. Specifically, reclamation regulations require oil sands operators to return disturbed lands to ''equivalent land capability,'' a portion of which will be reclaimed to wetland ecosystems (Government of Alberta 2017b).
Land reclamation occurs on active oil sands leases (both open pit mines and in-situ operations) concurrent with bitumen extraction activities. Most reclamation in the AOSR has targeted the re-establishment of upland forest habitat that is similar to the surrounding, naturally occurring boreal forest Hawkes and Gerwing 2019), although work aimed at directly reclaiming/restoring wetlands has also occurred. To date, wetland reclamation efforts have resulted in the establishment of numerous shallow open water and marsh wetlands on several surface mining leases including the first two mines established in the region: Syncrude's Mildred Lake lease and Suncor's Base Plant (Daly et al. 2009;Daly et al. 2010;Daly 2011;Pollard et al. 2012;Daly et al. 2012;Mollard et al. 2015;Ketcheson et al. 2016;Vitt et al. 2016;Scarlett et al. 2017). As of 2018, for example, Suncor had permanently reclaimed over 2000 ha of disturbed land at its Base Plant, of which * 48 ha was constructed mineral wetland. Some attention has also been given to the reclamation of peatland ecosystems, particularly of fens, on these leases (Ketcheson et al. 2016;Wytrykush et al. 2012;Cooper 2015, 2019;Vitt et al. 2016). Compensation efforts for peatland (fen) loss have included constructed fen creation using peat material harvested from natural peatlands (Ketcheson et al. 2016;Borkenhagen and Cooper 2019).
Although unintended, soil reconstruction methods used to recreate upland forestlands can also produce wetland-forming conditions (Little-Devito et al. 2019). Upland reclamation generally involves application of peat-mineral mix or forest floor materials as a capping substrate on reclaimed land (Hemstock et al. 2010;MacKenzie and Quideau 2010;Pinno and Hawkes 2015). Slight irregularities in the placement of different capping materials (topsoil) over large ([ 100 ha) areas of uneven terrain results in a heterogeneous surface that produces patches of very wet to very dry microsites (Gingras-Hill et al. 2018). This heterogeneity creates opportunities for wetland plant communities to establish spontaneously, such as in depressions within peat-mineral mix applications (Trites and Bayley 2009;Gringas-Hill et al. 2018;Little-Devito et al. 2019) but potentially even on flat terrain where fine-textured soils with low water storage potential promote frequent surface saturation (Little-Devito et al. 2019). Under certain conditions, lateral groundwater flow from constructed watersheds can also maintain soil wetness sufficient for wetland development Ketcheson et al. 2017). Such wetlands can enhance the suitability of reclaimed upland landforms for wildlife, increase biodiversity, improve water quality, and facilitate recharge and water storage capacity (Alberta Environment 2000; Ketcheson et al. 2016;Calhoun et al. 2017).
The study of opportunistic wetland development in the AOSR is relatively recent. On Syncrude's Mildred Lake Lease, opportunistic wetlands on three reclaimed landforms were investigated to assess controls on their structure and formation and the role they play in the ecohydrology of constructed landforms in the AOSR (Little-Devito et al. 2019). At Suncor's Base Plant, opportunistic wetlands occurring on reclaimed sites have not been formally included in the current inventory of reclaimed wetland areas and remain largely unquantified. In particular, the total reclaimed area constituted by these features is not well-understood, nor is the contribution they could potentially be making to the reclamation or partial replacement of wetland habitat impacted by bitumen extraction activities in the AOSR.
Here we report the results of a two-phase study of opportunistic wetlands that formed on landforms reclaimed to an upland forest type on Suncor's Base Plant. The objectives were to: (1) use aerial photo interpretation and remote sensing analysis to identify, classify, and delineate opportunistic wetlands that have developed on upland reclamation areas; (2) conduct a similar desktop delineation of intentionally reclaimed (constructed) mineral wetlands; and (3) use follow-up ground surveys to verify the accuracy of the desktop delineations, determine seasonal permanence of wetlands, characterize wetland vegetation communities, and compare these communities across the different reclamation areas. Results from this work are considered in the context of reclamation requirements in the AOSR and the contribution that opportunistic wetlands forming on upland reclaimed landforms make in mitigating impacts to the landscape resulting from bitumen extraction. The role that opportunistic wetlands have in enhancing the suitability of upland reclaimed landforms for wildlife is also briefly discussed along with the potential problems they present for the structural integrity of permanently reclaimed upland landforms.

Study area
Suncor Base Plant is located approximately 40 km north of Fort McMurray, Alberta, Canada in the Central Mixedwood subregion of the Boreal Forest Natural Ecoregion of Alberta (Natural Regions Committee (NRC) 2006) and the Boreal Plains ecoregion of Canada. Of the 2000 ha of permanently reclaimed land at the mine, * 1209 ha were remotely surveyed as part of this study. For purposes of wetland delineation, the study area was divided into 11 focal areas (Fig. 1). Each of these Areas of Interest (AOI) contained at least one of the following: (1) intentionally reclaimed (constructed) wetlands; (2) opportunistic wetlands (e.g., Fig. 2); or (3) topographic wetness information (derived with GIS; see ''Methods''). Except for intentional wetlands, which were included in this study for comparative purposes, wetland reclamation was not an explicit consideration for the upland landforms and all wetlands sampled formed opportunistically on the landforms reclaimed to an upland forest type.

Phase 1: desktop delineation
Using available remote sensing data (satellite imagery, orthophoto mosaics, and digital elevation models), we identified potential opportunistic wetlands on landforms reclaimed to an upland forest type. All suspected wetland boundaries within the AOIs were manually mapped in a GIS (Fig. 1). Desktop delineation followed the Alberta Wetland Identification and Delineation Directive and the Guide for Assessing Permanence of Wetland Basins (Government of Alberta 2015; AEP 2014).
To predict wetland occurrence, we assessed eight orthophoto mosaics from 2007 to 2014 and 67 GeoEye-1 and WorldView-2 satellite images from 2008 to 2015 for the presence of wetlands or places where wetlands might occur. Available satellite images spanned all calendar months but did not include all possible month-year combinations (Novoa and Hawkes 2016). Most images were composed of three colour bands corresponding to the visible spectrum. A few had a fourth, infrared band, which is useful for detecting water content and for discriminating vegetation types. Three digital elevation models were used, including two from 2012 (terrain and surface models) and one from 2014 (terrain model). All elevation models were generated from LiDAR datasets and had \ 1 m of spatial resolution.
Once areas with obvious water accumulation had been manually delineated in GIS, a topographic wetness index (TWI; Grabs et al. 2009) was calculated to improve the visual recognition of areas prone to moisture accumulation. We used TWI because the spectral signature of a wetland in the visible spectrum is not optimal for image classification algorithms (Amani et al. 2018) and produces low accuracy when using classical remote sensing methods. To calculate the TWI, we used a LiDAR digital terrain model (DTM) from October 2014 with an assumed (because metadata were not available) vertical accuracy of B 30 cm, which we compared to colour infrared aerial photography from September 2014. These were the most recent instances where spectral and elevation data were available for the same year and were thus considered the best available data sources for detecting potential wetlands.
The TWI was used to assign three wetness categories (low, medium, high) to areas with DTM coverage (Fig. 1). Category thresholds were determined from the mean and standard deviation of all TWI values. The mean was used to define the lowest wetness threshold (TWI values 4.71-5.12), the mean plus one half of the standard deviation defined the medium threshold (TWI values 5.12-5.54), and mean plus one standard deviation defined the upper (high) threshold (TWI values [ 5.54). TWI areas with values less than the mean were assigned a wetness category of nil. To delineate wetlands, the aerial photograph was displayed on screen using a colour infrared band (CIR) combination to highlight moisture and water bodies. Because the wetness index was used only as a visual aid and not to automate the delineation process, no comprehensive error analysis was undertaken. In areas without DTM coverage, desktop wetland delineation was accomplished using visual interpretation of the on-screen aerial imagery from September 2014; this method took more time but produced results similar to ones obtained for DTM coverage areas.
Opportunistic wetlands were preliminarily assigned to a wetland class of either shallow open water or marsh as per Alberta Environment and Sustainable Resource Development (AESRD 2015). Areas with darker, homogeneous pixels (i.e., visible standing water) were classified as shallow open water and areas with heterogeneous, dark pixels were classified as marshes. Ambiguous cases, where the wetland could not reliably be typed from the aerial imagery, were temporarily mapped as unspecified until they could be field verified. A fourth designation, intentionally reclaimed, was applied to constructed mineral wetlands. The same monthly satellite imagery was then used to determine the water permanence type (Alberta (1) within-year permanence, or the number of weeks per year that a wetland was wet (Table 1); and (2) year-over-year wetted status, which was the number of years, between 2008 and 2015, that the wetland was wet. If a wetland was wet for C 1/3rd of the years assessed, it was considered continuous; if wet for \ 1/3rd of the years assessed, it was considered noncontinuous. If ponding occurred between May and November in year t and continued to be observed in May of the following year, it was assumed that water was also present from December to April, although likely frozen and under snow. The results of the desktop exercise were used to plan and implement the field verification (phase 2).

Phase 2: field verification
From the dataset created in phase 1 (Novoa and Hawkes 2016), we selected a subsample of four intentionally reclaimed, and 100 opportunistic, wetland polygons to sample for field verification. After first stratifying by wetland type and permanence, the  The Alberta Wetland Classification System (AWCS) wetland types were used to modify the classification of opportunistic and intentionally reclaimed wetlands assessed in this study proportion of total polygons per stratum was used to determine the target sample size by stratum. We used a Generalized Random Tessellation Stratified (GRTS; Stevens and Olsen 2004) approach to produce a spatially-balanced sample dataset in which the distribution of sample sites mimics the population spatial distribution while avoiding the creation of clumped sampling sites (Kincaid and Olsen 2016). Using GRTS also supports design-based inferences across the entire study area. Sampling effort was focused on AOIs with the highest concentrations of identified wetlands ( Fig. 1) to minimize travel time while maximizing the number of unique wetland sites sampled. Field verification visits occurred in June and September 2017. A total of 238 phase 1 polygons were assessed, a number that exceeded the 100 initially selected for sampling (Table 2) because many mapped polygons turned out to be part of a single larger wetland complex when delineated in the field. Consequently, multiple polygons were often captured within a single field delineation. At each sample location, the desktop-identified feature was confirmed to be wetland or not wetland based on hydroregime and vegetation indicators. Each confirmed wetland was classified as to class, form, and permanency type using keys in AESRD (2015) and AEP (2014). At all subsample locations, vegetation and (if needed) soil assays were completed to identify the point at which the wetland ended, and the upland habitat began. For wetlands that had high topographic relief or obvious transitions in vegetation from wetland to upland, vegetation alone was used to verify the desktop interpretation of the wetland boundary. In contrast, wetlands with gentle relief and shallow shorelines tended to have transitional plant communities containing a mix of both wetland and upland species, especially towards the outer extent of the wetland. In these instances, a combination of soil and vegetation indicators was used to determine the true wetland extent (Government of Alberta 2015). Vegetation structure and composition were assessed within 5.64 m-radius (100 m 2 ) circular plots at one or more representative points around the suspected wetland perimeter. The number of points sampled was determined by the complexity of the site; heterogeneous wetlands with multiple community types required a higher sampling rate than more homogenous sites. At each point, a series of circular plots extending across both sides of the suspected wetland boundary was rapidly assessed to pinpoint the precise wetland extent using vegetation indicators identified in Government of Alberta (2015).
Assessments were repeated as the surveyor moved upslope until wetland indicators were replaced by predominantly non-wetland indicators, indicating that the landward extent of the wetland had been reached. This transition point was marked with GPS and temporarily flagged to help guide the subsequent GPS delineation. To provide a permanent data record of the wetland versus non-wetland condition near the wetland boundary, detailed structural and compositional information on two plots, one immediately below the demarcation line, and one immediately above the demarcation line, was recorded using a modified version of the field form in Government of Alberta (2015; App. 7). A total of 72 100 m 2 detailed vegetation plots were sampled near the outskirts of delineated wetlands (including intentional reclaimed wetlands) across the entire study area. If necessary, soil properties were examined on both sides of the suspected wetland:upland boundary to aid in pinpointing the precise wetland extent. A separate soil core was taken at the approximate centre of each vegetation plot to look for evidence of primary soil indicators (Government of Alberta 2015). These primary soil indicators included characteristics of water-altered or hydric soils such as gleying or mottling; presence of organic soils or organic horizons; presence of peat; oxidized rhizospheres, mineral accretions, or sulphuric odors; remains of aquatic invertebrates; and other soil properties associated with a prolonged seasonally high water table (Government of Alberta 2015). If any primary soil indicators of a wetland were present, the area was recorded as wetland. If neither vegetation nor soil indicators were present, the area was recorded as non-wetland.
To assess if plant community composition (presence/absence and percent cover) differed between different reclamation areas and between wetland sites within locations both above and below the wetland perimeter, we used permutational multivariate analysis of variance (PERMANOVA) in combination with non-metric multidimensional scaling (nMDS) plots. Wetland community development relative to prevailing soil properties was likewise assessed for that subset of wetlands for which soil information was available. Data analyses were performed in PRIMER with the PERMANOVA add-on (Anderson et al. 2008;Clarke and Gorley 2015). nMDS was used with 100 restarts (Clarke 1993), and the response variables for nMDS plots were resemblance matrices calculated using Bray-Curtis coefficients (Clarke 1993) constructed from vegetation species composition data. All plots had a stress (goodness of fit) B 0.2 and were considered good 2-dimensonal representations of higher dimensional trends (Clarke 1993).
Where ponded water was present, water physicochemical properties (pH, dissolved oxygen [DO; mg/ L], conductivity [lS/cm], and temperature [°C]) were measured 30 cm below the surface. DO, conductivity, and temperature were collected using a YSI Professional Series 2030 Multi-function Metre (Model Pro2030). An Oakton waterproof pHTester 30 was used to measure pH.
Wetland perimeters were georeferenced in the field using a TrimbleÒ GPS unit with sub-metre GPS accuracy. If the field delineation was within ± 5 m of the desktop delineation, the desktop delineation was considered verified (Government of Alberta 2015). To characterize internal vegetation zonation patterns, one to several supplementary transects were run from the centres of a selection of sampled wetlands uphill as far as the upland interface (or, in some cases, oriented so as to traverse a mosaic of wetland/upland habitats). Detailed, running plant species lists were recorded along each transect. Any evident transition points between vegetation zones (e.g., wet meadow and emergent); between different wetland classes (e.g., marsh and swamp); or between wetland and upland habitat, were georeferenced. A total of 188 supplementary transects were sampled. Once mapped, these transitions, in combination with the delineated wetland boundaries, served as a useful visual aid during subsequent orthophoto interpretations (Fig. 3).
Field-verified wetland areas were compared to aerial estimates from the desktop exercise to assess the accuracy of the latter. Observed differentials between the two delineation exercises were then used to extrapolate across areas that had not been fieldverified to derive estimates of total wetland extent for each AOI. Extrapolations were made by weighting the desktop estimate for each AOI (equating to the total area of all mapped polygons) by a correction factor (CF). The CF was calculated by dividing the total wetland delineated area (desktop wetland area verified ? additional wetland area field-delineated) by the desktop identified wetland area. The corrected total wetland area (ha) was calculated by multiplying the desktop-identified wetland (ha) by the CF. The total area of opportunistic wetlands was then estimated by summing the individual corrected values over all AOI (Table 6, ''Results''). The CF provides an estimate of the extent of wetlands in areas where field verification did not occur. As such, field verification is needed to produce a more accurate estimate of the extent of wetlands in those areas.

Phase 1
Prior to this study, Suncor had independently delineated 48.0 ha of intentionally reclaimed wetlands at its Base Plant mine. During the desktop exercise, a total of 687 polygons representing 41.9 ha of opportunistic wetlands (* 3.5% of the area assessed) were mapped (Table 2). Wetlands were assigned to one of two wetland classes, marsh (n = 90) or shallow open water (n = 300), or else left as ''unspecified'' (n = 283). The remaining 14 wetlands were classified as intentionally reclaimed (Table 2). Further typing for water permanence yielded nine different mapping categories for opportunistic wetlands, plus one category for intentionally reclaimed wetlands (Table 2). Based on the monthly satellite imagery, most marshes (n = 70), all shallow open water (n = 283), and the majority of unspecified wetlands (n = 197) were assessed as having a water permanence type of semi-permanent (Table 2). On a year-over year basis, all assessed marsh and shallow open water areas, and most unassigned wetlands (78%), were continuously inundated through the study period (2008-2015; Table 3). Noncontinuous (i.e., those wetlands wetted for \ 1/ 3rd of the years assessed) occupied 2.37 ha or * 6% Fig. 3 Example transect from a field-delineated, opportunistic wetland complex in area of interest (AOI) 4 ( Fig. 1) showing transitions between shallow marsh, shrubby swamp, and upland habitat. Red polygons represent the field-verified perimeters; turquoise aligns with the desktop delineations of the total wetland area assessed for year-over-year continuity (Table 3).

Field delineations
Confirmed wetlands fell into one of three wetland classes (AESRD 2015): shallow open water, marsh (shallow form), or swamp (shrubby form). In cases where shallow open water was identified, this class generally formed a subcomponent within a larger area of marsh. Shallow open water was generally limited in extent and could either support minor amounts of submersed and/or floating aquatic vegetation or be unvegetated. Based on conductivity readings (AESRD 2015), standing water ranged from slightly brackish to moderately brackish, with pH ranging from 7.2 to 8.7. In contrast to the desktop prediction (Table 2), fieldverified shallow open water wetlands tended to be of seasonal duration, while only a few were classed as semi-permanent.
Of the 238 polygons field-verified, seven were reclassified as non-wetland following field inspection, representing 2.9% of the polygons sampled. The misidentified features included dry swales, dried up mud depressions, dry ravines, and wood debris pilesfeatures that, in the desktop analysis, presented as regions of darker pixelation compared to the surrounding area. The remaining 231 locations were verified as classifiable wetlands (or parts of wetlands). Of these, 80.8% of the originally unspecified wetlands were classed as marsh, 7.0% were classed as shallow open water, and 12.3% were classed as swamp. Of those initially mapped as marsh, 4.7% were reclassified as shallow open water and 1.4% were reclassified as swamp. Of those mapped as shallow open water, 37.5% were reclassified as marsh and 1.8% were reclassified as swamp (Table 4).
With respect to water permanency typing, most (88.5%) unspecified wetland area was verified as seasonal during ground surveys, with 3.7% of this area identified as temporary and 3.7% identified as semipermanent (Table 5). Most (86.1%) of the area identified in the desktop exercise as temporary was also retyped as seasonal, along with * 15% of the area initially typed as semi-permanent (Table 5). Overall, 86% of the area covered by water decreased significantly between June and September and in most cases surface water was completely absent by September. Assuming these observations are reflective of normal seasonal variation, they support the prediction that desktop estimations of standing water based on aerial photo interpretation will be influenced by acquisition date (both month and year) of the aerial photos used (Novoa and Hawkes 2016).
The desktop analysis was successful 96% of the time in identifying the presence of marshes and shallow open water wetlands (Miller et al. 2017). However, wetland aerial extents were consistently underestimated by the desktop delineations. The underestimation percentage ranged from * 35% for AOI 9 up to * 98% for AOI 2, with an overall underestimation rate of * 89% for all areas sampled. Polygons delineating non-inundated wetland habitat had a median underestimation of * 78%, while those delineating standing water (including emergent and shallow open water zones) had a median underestimation of * 16%. Differentials were due in part to underestimations in the extent of non-inundated vegetation zones (e.g., wet meadows) associated with marsh and shallow open water wetlands, and in part to the under-detection of shrubby swamps. The latter, of which * 59 ha were identified and mapped during fieldwork, were non-ponded areas of elevated soil moisture content supporting a plant cover dominated by obligate or facultative wetland shrubs (e.g. Salix bebbiana, Alnus viridis ssp. crispa, Populus balsamifera) and grasses (e.g. Calamagrostis canadensis), with a sparse component of obligate wetland sedges (e.g., Carex aquatilis and C. utriculata) and various facultative upland species. This widespread wetland class formed a component of 71 field-verified polygons, representing * 30% of the 238 polygons surveyed. The total area of opportunistic wetlands delineated during the phase 1 desktop exercise per AOI was between 0.65 and 22.93 ha, for a total of 41.9 ha across AOIs (Table 6). Field verifications (phase 2) led to the delineation of an additional 0.09 to 29.82 ha (total of * 65.8 ha) across AOIs. Combining the two results, between 0.19 and 30.43 ha of opportunistic wetlands were delineated per AOI for a total of * 74.0 ha. Extrapolating across unsampled areas, corrected estimates for individual AOIs were between 1.3 ha (AOI 9 in Fig. 1) and 95.5 ha (AOI 11 in Fig. 1), resulting in an overall estimate for Suncor's Base Plant mine of * 210 ha ( Table 6). The relatively high CF obtained for AOI 2 (49.7) was related to the identification in the desktop exercise of what appeared to be several small, distinct wetlands. Field verifications revealed these isolated areas to be mostly patches of open water nested within a much more extensive area of swamp and marsh (as determined by vegetation and soil sampling) not readily detected in the remote imagery. The total area mapped as opportunistic wetland habitat represents * 17.4% of  the overall area considered (1209 ha)-an increase of [ 400% over the total wetland area (including intentionally constructed wetlands) recognized for the mine site prior to this study (i.e., 48 ha). Of the opportunistic wetlands delineated, most (n = 617) were \ 1 ha in area (Fig. 4). Of those [ 1 ha (10,000 m 2 ), only 12 exceeded 5 ha (Fig. 4). The specific plant compositions recorded on either side of the wetland boundary varied among AOIs (Fig. 1). This variation can be seen in the segregation of plant cover data by AOI, when samples from the two habitat types (upland and wetland) are plotted separately with nMDS ordinations (Fig. 5). For example, Wapisiw (AOI 7), MD5 (AOI 4), and NSE (AOI 11) samples form a distinct cluster on the left, right, The proportion of each AOI sampled during fieldwork is shown along with the correction factor (CF), which was calculated by dividing the total Phase 2 delineated area (Phase 1 area sampled ? additional wetland area field-delineated) by the original Phase 1 area. E.g., for AOI 2: (0.61 ? 29.82)/0.61 = 49.71. The corrected total wetland area (ha) was calculated by multiplying the desktopidentified wetland (ha) by CF (e.g., for AOI 2: 0.65 9 49.71 = 32.1) Fig. 4 Area distribution of mapped opportunistic wetlands (desktop-identified ? field-delineated) on habitats reclaimed to an upland landform on Suncor's Base Plant mine in northeastern Alberta and bottom regions of the upland diagram, respectively (Fig. 5a). In the wetland diagram, MD5 and NSE samples cluster at the top and bottom, respectively, with MD2 and MD8 (AOI 3) samples clustering closely (but still separately) near the middle (Fig. 5c). The variation in vegetation species composition among sites was statistically significant ( Table 7), suggesting that wetland development is sensitive to the specific substrates used during the reclamation process.
Differences in wetland community development among the AOIs are associated with topographic variation (e.g., elevation and % slope) in combination with the specific reclamation substrates that were placed over tailings prior to seeding/planting (Fig. 6). Of the various abiotic parameters modeled (e.g., soil class, surface material, horizon depths [upper/lower], drainage, aspect, slope, and elevation), the most influential variables on community composition were soil class and surface material. For example, presence of Rego Gleysol and Rego Humic Gleysol (R.G and  Total 61 -----Variance components represents the percent of the observed community variation accounted for by each factor. Significant variables are denoted in bold placed (Anth-PMM) tend to be associated with more transitional (facultative) wetland plant associations (e.g., Calamagrostis canadensis and Salix bebbiana; Fig. 6).

Discussion
Field surveys in 2017 confirmed the finding of Novoa and Hawkes (2016) that reclamation areas on Suncor Base Plant support an abundance of opportunistic wetland habitat. We estimate that * 210 ha of wetlands (a combination of shallow open water, marshes, and swamps) have spontaneously developed on reclaimed upland landforms at the mine, representing * 17% of the 1,209 ha considered in this study and a greater than four-fold increase over the total wetland area (including intentionally constructed wetlands) recognized for the mine site prior to this study. The proportion of area comprising opportunistic wetlands on Suncor's Base Plant was greater than that reported for Syncrude's Mildred Lake Lease (* 8%) although the methods used to assess wetland occurrence and distribution at Syncrude differed (Little-Devito et al. 2019). While the desktop (remote sensing) exercise used for phase 1 of this study was successful at predicting the existence of wetland conditions on the ground over 96% of the time, it was less successful at predicting the spatial extent of those wetlands. Wetland extents were consistently underestimated by the desktop delineations, due mainly to underestimations in the extent of non-inundated vegetation zones (e.g., wet meadow) as well as shrubby swamp. Polygons indicating noninundated wetland habitat had a median underestimation of * 78%, while those delineating standing water (including emergent and shallow open water zones) had a much smaller, median underestimation of * 16%. The accuracy of desktop delineations may also have been affected by differences in acquisition dates among the aerial photos, digital elevation model (DEM), and field samples, and by the resolution of the existing DEM. The lack of a suitable available method for calibrating the topographic wetness index (TWI) was likely also a factor affecting mapping accuracy (Novoa and Hawkes 2016). Further explorations with GIS approaches could yield improvements to the desktop methodology. For example, the inclusion of other topographic indexes such as the downslope distance gradient (Hjerdt et al. 2004) could increase the accuracy of wetland detection and delineation. A cursory examination of juxtaposed examples of the TWI and the downslope index, applied to the same opportunistic wetlands at MD2, suggest that the two indices could play a complementary role in identifying wetland locations and extents (Miller et al. 2017). Additional investigations are needed to determine if these different indices can be calibrated with each other in a way that improves wetland detectability (Miller et al. 2017).
Terrestrial and emergent plant communities inhabiting opportunistic and intentionally reclaimed wetlands at Suncor's Base Plant mine were comparable, in terms of vegetation structure, species composition, and species richness, to naturally occurring mineral wetlands in the region (AESRD 2015;Miller et al. 2017). However, submersed and floating-leaved plant communities (e.g., Ceratophyllum demersum and Potamogeton/Nuphar communities) were generally poorly developed to absent in the opportunistic shallow open water wetlands we surveyed (albeit from a distance, as we were not boat-equipped). Low percent cover of aquatic macrophytes may be associated with reduced productivity, possibly reflecting low nutrient levels or salinization . Furthermore, at the scale of the study area, wetland communities varied significantly among each area of interest. Successional trajectories appear to be mediated by topographic variation (e.g., in elevation and % slope) in combination with the reclamation substrates that were placed prior to seeding/planting (Fig. 6). For example, presence of Rego Gleysols, Rego Humic Gleysols, and placed overburden was correlated with higher cover of obligate wetland species such as common cattail and great bulrush, while peat-mineral caps tended to be associated with more transitional (facultative) wetland plant associations. More extensive data are required on the reclamation substrates used at various wetland sites to understand how closely wetland plant composition is correlated with substrate type. Nevertheless, our preliminary results suggest that wetland development is sensitive to the specific substrates used at different locations, and that prior knowledge of substrate types will be important when making predictions about wetland development trajectories.
Most of the reclaimed landforms assessed for this study are between 6-and 19-years post-reclamation and can be characterized as early seral. The possible role that opportunistically forming wetlands can play in mitigating impacts to wetland habitats stemming from bitumen extraction in the Athabasca Oil Sands Region has only recently been considered (this study and Little-Devito et al. 2019). As noted by Little-Devito et al. (2019), the presence of opportunistic wetlands could help fulfil some of the reclamation requirements associated with wetland habitat loss mitigation. Their relative ubiquity at Suncor's Base Plant, combined with their seasonal persistence and year-over-year continuity, suggests that opportunistic wetlands may have a greater influence on the closure landscape than previously recognized. The conditions that seem to promote spontaneous wetland development, such as topographical contouring, heterogeneity of soil cover placements (in terms of depth and composition), and elevated water-holding capacity, may be an inevitable if unintended outcome of current land reclamation practices at Suncor's Base Plant mine and Syncrude's Mildred Lake Lease and, by extension, at other AOSR mines as well.
As more areas are reclaimed to upland forest types planted on a peat-mineral or a forest floor mineral mix substrate, more opportunistic wetlands are prone to develop in situ. The spontaneous appearance of wetlands on reclaimed areas presents an opportunity to consider the intentional facilitation of such features (e.g., via surface contouring or targeted plant seeding) as a potential complementary approach to compensating for habitat impacts to the pre-mine boreal landscape. Certainly, the presence of numerous small wetlands (at least 619 at Suncor Base Plant) has the potential to increase the ecological heterogeneity of upland reclaimed landforms, increase wildlife habitat suitability, and provide hydrologic functions (Little-Devito et al. 2019; Hawkes and Gerwing 2019). For example, Hawkes and Gerwing (2019) found that wildlife (i.e., small mammal, large mammal, songbird, amphibian) communities on reclaimed uplands across the AOSR exhibited 31-62% (average 52%) similarity within 33 years of reclamation. It is reasonable to surmise that the time required for wildlife utilization of opportunistic wetlands to begin to converge with that of regional natural wetlands could be even less, given that shallow open water, marshes, and shrubby swamps can form quickly (within a few years) and probably reach functional maturity within a much shorter timeframe than mixedwood forest. One environmentally sensitive, wetland-dependent taxonomic guild that has already succeeded in occupying opportunistic wetland habitat at Suncor is amphibians [including Wood Frog (Lithobates sylvaticus), Boreal Chorus Frog (Pseudacris maculata), and Canadian Toad (Anaxyrus hemiophrys); unpubl. Data]. Colonization by amphibians has occurred despite the numerous and obvious barriers to dispersal (e.g., open mines, roads, tailings) that exist between many reclamation sites and neighbouring occupied natural sites, reinforcing the idea that opportunistic wetlands could have a role in assisting oil sands operators to meet their regulatory requirements and obligations around wildlife habitat mitigation.
Most wetlands sampled (93%) were less than 1 ha in size (Fig. 4). The maximum water table of these wetlands is just below or above the ground surface, and while their presence could pose a threat to the geotechnical stability of constructed upland landforms, this appears unlikely due to their shallow profiles and high evaporative potential (CEMA 2014; Little-Devito et al. 2019). A recent internal study by Suncor corroborates this, at least for two reclaimed landforms on their Base Plant mine (areas 3 and 4 in Fig. 1). Both upland areas are mining waste dumps without tailings. Wetlands were evaluated relative to their location (e.g., position relative to the crest of a slope) and water depth to determine if their presence represented a threat to the structural stability of the landforms. The internal report did not identify any significant concerns stemming from the occurrence or distribution or opportunistic wetlands for the areas assessed but future monitoring was recommended (Suncor, unpubl. data). However, given that many reclaimed uplands are built with mine tailings necessitating that they be capped, the potential short-and long-term risks to landform integrity need further investigation. Likewise, the longevity of opportunistic wetlands relative to forest succession and changing climatic regimes is presently unknown (Little-Devito et al. 2019). As adjacent planted or naturally establishing upland plant communities begin to mature and increase their draw upon the water table through evapotranspiration, many smaller marshes and shallow open water wetlands will likely fill in and succeed first to shrubby swamp before transitioning to upland forest or drier, graminoid-dominated shrubland over time. During field verifications we observed various signs of successional transition already underway, for example, in the form of willows (Salix spp.) establishing over former cattail stands and encroachment of mixed woody shrubs into wet meadow zones of marshes (Fig. 2c). Likewise, large areas of delineated swamp habitat appeared to be situated well above the water table and supported both planted spruce (Picea glauca) as well as stands of regenerating aspen and poplar (Populus tremuloides, P. balsamifera), all common components of adjacent mixedwood stands.
Given that most of the reclamation planned for the AOSR is scheduled to occur post-2035 (Pickard et al. 2013), greater understanding of the interplay between opportunistic wetland development, upland reclaimed landforms, wildlife use, and wetland persistence on those landscapes is required for purposes of reclamation tracking as it pertains to mine closure plans and habitat reclamation.

Conclusions
Remote imaging can be an effective tool for identifying instances of opportunistic wetland formation on reclaimed upland landforms in the Athabasca Oil Sands Region, but additional field truthing will probably be needed to accurately delimit the extents of individual wetland complexes. This applies to both intentionally reclaimed (constructed) and non-intentional (opportunistic) wetlands. Even when ground surveys are used, deriving precise wetland coverage estimates for a single large reclamation area (more so for the entire mine area) can be challenging due to both the complexity of the terrain at the micro-scale and the large aerial extents occupied by certain wetland classes such as shrubby swamp. It is thus hoped that further explorations with remote sensing approaches will yield improvements to desktop-based mapping. For example, the inclusion of other topographic indexes aside from TWI (topographic wetness index) such as the downslope distance gradient could increase the accuracy of wetland detection and delineation.
Wetland and forest ecosystems are affected by bitumen extraction in the AOSR. Oil Sands operators have an obligation to reclaim and this requires the creation of functional habitats with equivalent land capability. Most reclamation that has occurred to date has focused on reclaiming constructed landforms to an upland forest type consistent with the surrounding and naturally occurring boreal forest (Pinno and Hawkes 2015;Hawkes and Gerwing 2019). Although unintended, the interspersion of opportunistic wetlands with upland reclaimed landforms more closely reflects the pre-disturbance upland landscape, which was characterized by a matrix of peatlands, forested ecosystems, and mineral wetlands. The inclusion of opportunistic wetland development in reclamation tracking and closure planning merits consideration as does the opportunity to manipulate current reclamation practices to promote the establishment and persistence of wetlands on reclaimed landforms (Little-Devito et al. 2019). When combined with purposeful (intentional) wetland creation, the formation of opportunistic wetlands on reclaimed upland landforms can contribute to the development of forest-peatlandmineral wetland mosaics that were once abundant in the mineable portion of the AOSR, thereby assisting with the process of mitigating the impacts associated with this particular natural resource extraction.