Dynamics of microbial populations and diversity in NAPL contaminated peat soil under varying water table conditions.

Despite the risks that hydrocarbon contamination from pipeline leaks or train derailments impose on the health of peatlands in hydrocarbon production areas and transportation corridors, assessing the effect of such contaminations on the health and sustainability of peatlands has received little attention. This study investigates the impacts of hydrocarbons on peat microbial communities. Column experiments were conducted on non-aqueous phase liquid (NAPL) contaminated undisturbed peat core (0-35 cm) under static and fluctuating water table conditions. Water table fluctuations reduced residual NAPL saturation from 8.1-11.3% to 7.7-9.5%. Biodegradation of n-C8 and n-C12 along with oxidation of CH4 together produced high CO2 concentrations in the headspace. Clear patterns in dynamics in the microbial community structure were observed, with a more pronounced population growth. However, a significant loss of microbial richness was observed in contaminated columns. The result indicates that the phylum Proteobacteria benefited most from NAPL; however, their families differed between static and fluctuating water table conditions. This study established strong evidence that peat microbes and water table fluctuation can be an excellent tool for hydrocarbon removal and its control in peatlands.

Aquifer depressurization and water table lowering induces landscape scale subsidence and hydrophysical change in peatlands of the Hudson Bay Lowlands.

The depositional history of the Hudson Bay Lowlands (HBL) in Ontario, Canada has created a low relief, poorly drained landscape, favouring the formation of one of the largest peatland complexes in the world. High volume dewatering associated with resource extraction in this area, such as the De Beers Victor Diamond Mine, tests the ability of the underlying confining layer to limit water losses in the peatlands above. This research quantifies the deepening of water tables and increase in effective stress related to mine dewatering and the resulting changes to bog and fen peatland hydrophysical structure and function. Long-term implications of these impacts are discussed. One impacted and two unimpacted transects were instrumented for meteorological (precipitation and evapotranspiration) and hydrophysical (hydraulic head, hydraulic conductivity (Ksat), and surface elevation) monitoring over a 12-year period in the vicinity of the Victor Mine. Over this study period, the unimpacted peatlands operated within relative hydrological equilibrium, demonstrated through shallow water tables, negligible subsidence, and stable Ksat. Contrastingly, all impacted peatlands experienced deeper watertables, larger downwards gradients, and measurable long-term subsidence (4-15 cm). Hydrological impacts were highest in bogs with a thin underlying confining layer even if they were farther from the point of dewatering, highlighting the need for environmental monitoring programs which incorporate an assessment of aquitard thickness. Where subsidence occurred, associated decreases in Ksat deflected bog-fen-tributary flow-paths deeper, reducing the upwards transport of solute rich water to downgradient fens. The long-term effects of these landscape scale changes should be studied further, particularly since climate change in this region will potentially increase water deficits and further alter peatland connectivity. Peatland studies should be conducted in different landscapes experiencing water table lowering due to drought or depressurization in order to better understand the associated subsidence patterns and hydrophysical changes in varying geological and morphological regimes.

Multiphase flow behavior of diesel in bog, fen, and swamp peats.

Hydrocarbon fate and transport in various categories of peatlands is complicated by the botanical origin, and thus variations in the hydraulic structures and surface chemistry of its peat soils. There has been no systematic evaluation of the role of different peat types on hydrocarbon migration. Thus, two-phase, and three-phase flow experiments were performed for living and partially decomposed peat cores from bog, fen, and swamp peatlands. Numerical simulations of water drainage were performed using HYDRUS-1D, diesel-water and diesel-water-air flow using MATLAB Reservoir Simulation Toolbox (MRST). Five water table (WT) fluctuations were imposed to explore its potential to reduce residual diesel saturation in peat columns. Our results demonstrate a good match of the relative water permeability (krw) - saturation (S) relations estimated using the unsaturated hydraulic conductivity-S relation derived from HYDRUS-1D modeling of two-phase flow, and krw - S from MRST for three-phase flow, for all tested peat columns. Thus, we recommended using two-phase system based krw - S predictions if multiphase data are unavailable for peatland sites' spill management plans. We found the discharge of water and diesel both increase with increasing hydraulic conductivity, while residual water and diesel were within the range of 0.42-0.52 and of 0.04-0.11, respectively. High diesel discharge rates suggest that quick spill-response is required to manage its spread in peatlands. Up to 29% of residual diesel saturation was yielded by the five WT fluctuations, and thus we strongly recommend WT manipulation as a first step towards diesel decontamination progression in peatlands.

Projecting the hydrochemical trajectory of a constructed fen watershed: Implications for long-term wetland function.

Surface mining operations for bitumen have fundamentally altered large areas of boreal forest and fen peatland in the Athabasca Oil Sands Region (AOSR) of Alberta, Canada. Pilot projects intended to assess the feasibility of fen construction as a reclamation option have been designed, built, and are currently undergoing monitoring. Initial assessments of ecohydrologic function have been conducted for these systems but offer limited insight into their evolution and likely successional pathway. Thus, this study projects the hydrologic and geochemical behaviour of a constructed fen watershed to understand whether the system will be capable of supporting peatland processes into the future. A numerical groundwater flow and sodium transport model was calibrated and validated with 7 years of hydraulic head, water flux, and water chemistry data. Based on Monte Carlo simulations, the projected fen water table would be stable and remain close to the surface (<15 cm), indicating that the design of the system can generate sufficient water quantity to meet evaporative demand and maintain surface water discharge. However, water quality was more sensitive to climatic variability, which induced a large range in potential sodium concentrations at the fen surface (450-850 mg L-1). Evapoconcentration of salts across the surface of the fen will likely limit moss establishment for decades following construction. Yet stress-thresholds of salt-tolerant vegetation like sedges will not be exceeded. Ultimately, these projections support the original design principles and philosophy that guided the creation of the watershed. Nonetheless, this work indicates that increasing the area of the fen relative to the upland would not have a detrimental impact on the ability of the system to maintain a high water table. This could allow for the proportion of peatlands on the reclamation landscape to reflect the pre-disturbance environment more faithfully.

Characterizing the hydraulic and transport properties of a constructed coarse tailings sand aquifer.

Millions of tonnes of coarse tailings sand are produced every year as a byproduct of the bitumen extraction process in the Athabasca Oil Sands Region. These tailings materials contain residual quantities of mobile solutes, which can be transported through groundwater to downgradient terrestrial and aquatic ecosystems. The anticipated ubiquity of coarse tailings sand on the post-mined landscape necessitates the characterization of its hydraulic and transport properties. Hydraulic conductivity and dispersivity was evaluated at multiple scales, and included the first field-scale tracer test conducted in a tailings sand aquifer. Average hydraulic conductivity derived using laboratory cores, single-well response tests, and the tracer test were 3.2 m d-1, 2.9 m d-1, and 3.4 m d-1, respectively. These measurements demonstrated close agreement and were consistent with expectations of a material that experiences some grain-size segregation and homogenization due to the oil sands process and the nature of deposition. The field-scale tracer test appeared to obtain the asymptotic dispersivity of the coarse tailings sand aquifer, reaching a maximum value of 0.5 m after 18 m of displacement. Coarse tailings in the oil sands that experience similar processes of segregation, settling, and deposition on the reclamation landscape could be expected to have similar hydraulic properties.

Assessing benzene and toluene adsorption with peat depth: Implications on their fate and transport.

After a hydrocarbon spill in a peatland, dissolution of water-soluble compounds including benzene and toluene introduces a dissolved-phase plume to the peatland groundwater system, while the adsorption of these solutes onto the peat matrix restrains their distribution velocity. The adsorption of benzene and toluene and its dependency on peat depth, thus degree of decomposition, are investigated. The batch adsorption experiments revealed that benzene and toluene adsorption isotherms in peat are linear, with adsorption coefficients ranging from 16.2 to 48.7 L/kg and 31.6-48.7 L/kg, respectively. In a vertical peat profile benzene adsorption decreased with depth, while toluene adsorption increased. Considering toluene adsorption onto cellulose is significantly less than toluene adsorption onto humic substance, the increase in toluene adsorption was attributed to decreasing cellulose and increasing humic substances with depth. Negligible competition for adsorption was observed between benzene and toluene at the measured concentrations. The retardation factors of benzene and toluene ranged respectively from 3.5 to 10.7 and from 5.4 to 17.7, both increasing with depth. Higher retardation in deeper peat coupled with lower hydraulic conductivity will lead to a weaker solute velocity in deeper peat, thus preferential migration of these dissolved-phase contaminants in shallow layers. The results can help predict the behavior of dissolved hydrocarbons in peatlands after a hydrocarbon spill.

Temperature, moisture and freeze-thaw controls on CO2 production in soil incubations from northern peatlands.

Peat accumulation in high latitude wetlands represents a natural long-term carbon sink, resulting from the cumulative excess of growing season net ecosystem production over non-growing season (NGS) net mineralization in soils. With high latitudes experiencing warming at a faster pace than the global average, especially during the NGS, a major concern is that enhanced mineralization of soil organic carbon will steadily increase CO2 emissions from northern peatlands. In this study, we conducted laboratory incubations with soils from boreal and temperate peatlands across Canada. Peat soils were pretreated for different soil moisture levels, and CO2 production rates were measured at 12 sequential temperatures, covering a range from - 10 to + 35 °C including one freeze-thaw event. On average, the CO2 production rates in the boreal peat samples increased more sharply with temperature than in the temperate peat samples. For same temperature, optimum soil moisture levels for CO2 production were higher in the peat samples from more flooded sites. However, standard reaction kinetics (e.g., Q10 temperature coefficient and Arrhenius equation) failed to account for the apparent lack of temperature dependence of CO2 production rates measured below 0 °C, and a sudden increase after a freezing event. Thus, we caution against using the simple kinetic expressions to represent the CO2 emissions from northern peatlands, especially regarding the long NGS period with multiple soil freeze and thaw events.

Modelling the hydrologic effects of vegetation growth on the long-term trajectory of a reclamation watershed.

Reclamation watersheds that integrate fen peatlands into the design require the inclusion of uplands that are capable of supporting forest development while concurrently supplying sufficient groundwater recharge to downgradient wetland ecosystems. This necessitates selecting materials with suitable soil hydraulic properties and identifying the appropriate thickness and layering to fulfill the dual function of uplands as water storage, and water conveyance features. Currently, these systems incorporate tailings sand - a mine waste material - overlain by a cover soil of fine forest-floor material. The developmental pathway of these uplands is currently unknown, and it is unclear whether these landforms will provide enough groundwater recharge once a climax vegetation community establishes. Therefore, this research attempts to estimate the maximum density of vegetation, and associated water balance fluxes of a constructed upland integrated into a peatland watershed. The numerical modelling software HYDRUS-1D simulated soil moisture dynamics using a 65-year meteorological record, and a plant water stress algorithm was used to estimate the maximum sustainable leaf area index that the upland could support. Based on the thickness of the cover soil, the upland could support an average leaf area index of 1.2. Under this vegetation density, average annual groundwater recharge was 83 mm, and predominantly supplied by snowmelt (64%). Given this quantity of recharge, the model indicates that the upland will continue to provide enough groundwater to offset the anticipated water deficit in the downgradient fen ecosystem. However, by altering the design of the upland, specifically the spatial arrangement and thickness of cover soil, the same recharge could be supplied while also allowing for a higher average vegetation density. Such a design could allow for the creation of watersheds with a higher proportion of peatland.

Soil moisture dynamics modelling of a reclaimed upland in the early post-construction period.

Mine reclamation landscapes typically comprise layers of mine waste materials such as tailings sands, capped with a cover soil. In addition to the arrangement and placement of these materials, their hydraulic properties govern the performance of the built system. Soil evolution due to freeze-thaw cycling can result in dramatically altered soil hydraulic properties compared to the as-built material. Therefore, prediction of present and future hydrologic behaviour relies on understanding the nature and magnitude of this change and the elapsed time associated with stabilization. This research quantifies the transient hydraulic properties of mine reclamation materials at a constructed upland within a reclaimed watershed, and models the effect of this evolution on the partitioning of soil moisture between evaporation and groundwater recharge. Soil moisture dynamics were simulated using HYDRUS-1D for the ice-free period two, three, and five years after construction. A capillary barrier between the fine-grained cover soil and coarse-grained tailings sand regulated percolation past the interface. Soil evolution of the cover soil was responsible for an increase in saturated hydraulic conductivity by an order of magnitude, decrease in air-entry pressure by a factor of 4, and decrease in the van Genuchten n parameter by a factor of 2. The altered soil hydraulic properties associated with the weathered cover soil ultimately resulted in a 64% increase in groundwater recharge as a consequence of the capillary barrier weakening. The cover soil exhibited minor spatial heterogeneity in soil hydraulic properties, and did not contribute substantial uncertainty to the estimates of groundwater recharge and evaporation. Cover soil thickness exerted a strong influence on the partitioning of soil moisture. Reclaimed uplands will provide the most recharge to downgradient ecosystems in the period following the completion of soil evolution (~4 years) but preceding substantial vegetation development.

Characterizing the immiscible transport properties of diesel and water in peat soil.

Extensive pipeline and railway corridors crossing Canadian peatlands make them vulnerable to hydrocarbon spills, potentially impairing ecosystem health, so it is important to be able to forecast hydrocarbon fate and transport within and beyond the peatland. The redistribution of hydrocarbon liquids in groundwater systems are controlled by the multiphase flow characteristics of the aquifer material including capillary pressure-saturation-relative permeability (Pc-S-kr) relations. However, these relations have never been characterized for the hydrocarbon-water phases in peat. To address this, the flow and transport of diesel and water in peat soils were examined through a number of one dimensional vertical immiscible displacement tests, in which diesel was percolated into peat pore space displacing peat water, leading to a two-phase flow regime. Inverse modelling simulations using both Brooks and Corey's and power law relative permeability models, matched the data of the immiscible displacement tests well. Irreducible water saturation (Swirr) and the curvature of water relative permeability relation increased with peat bulk density. The residual diesel saturation (SNr) ranged between 0.3% and 17% and increased with bulk density of peat. In a given peat, SNr was a function of saturation history and increased with increasing maximum diesel saturation. The receding contact angles of water in water-air systems and diesel in diesel-air systems, respectively, were 51.7° and 61.2°, showing that the wetting tendency of peat in the air imbibition condition is toward the draining liquid. These experiments advance our knowledge on the behavior of hydrocarbons in peat, and can improve numerical modelling of hydrocarbon transport after a spill.

The distribution and migration of sodium from a reclaimed upland to a constructed fen peatland in a post-mined oil sands landscape.

Post-mine landscape reclamation of the Athabasca Oil Sands Region requires the use of tailings sand, an abundant mine-waste material that often contains large amounts of sodium (Na+). Due to the mobility of Na+ in groundwater and its effects on vegetation, water quality is a concern when incorporating mine waste materials, especially when attempting to construct groundwater-fed peatlands. This research is the first published account of Na+ redistribution in groundwater from a constructed tailings sand upland to an adjacent constructed fen peat deposit (Nikanotee Fen). A permeable petroleum coke layer underlying the fen, extending partway into the upland, was important in directing flow and Na+ beneath the peat, as designed. Initially, Na+ concentration was highest in the tailings sand (average of 232mgL-1) and lowest in fen peat (96mgL-1). Precipitation-driven recharge to the upland controlled the mass flux of Na from upland to fen, which ranged from 2 to 13tons Na+ per year. The mass flux was highest in the driest summer, in part from dry-period flowpaths that direct groundwater with higher concentrations of Na+ into the coke layer, and in part because of the high evapotranspiration loss from the fen in dry periods, which induces upward water flow. With the estimated flux rates of 336mmyr-1, the Na+ arrival time to the fen surface was estimated to be between 4 and 11years. Over the four-year study, average Na+ concentrations within the fen rooting zone increased from 87 to 200mgL-1, and in the tailings sand decreased to 196mgL-1. The planting of more salt-tolerant vegetation in the fen is recommended, given the potential for Na+ accumulation. This study shows reclamation designs can use layered flow system to control the rate, pattern, and timing of solute interactions with surface soil systems.

Competitive transport processes of chloride, sodium, potassium, and ammonium in fen peat.

There is sparse information on reactive solute transport in peat; yet, with increasing development of peatland dominated landscapes, purposeful and accidental contaminant releases will occur, so it is important to assess their mobility. Previous experiments with peat have only evaluated single-component solutions, such that no information exists on solute transport of potentially competitively adsorbing ions to the peat matrix. Additionally, recent studies suggest chloride (Cl-) might not be conservative in peat, as assumed by many past peat solute transport studies. Based on measured and modelled adsorption isotherms, this study illustrates concentration dependent adsorption of Cl- to peat occurred in equilibrium adsorption batch (EAB) experiments, which could be described with a Sips isotherm. However, Cl- adsorption was insignificant for low concentrations (<500 mg L-1) as used in breakthrough curve experiments (BTC). We found that competitive adsorption of Na+, K+, and NH4+ transport could be observed in EAB and BTC, depending on the dissolved ion species present. Na+ followed a Langmuir isotherm, K+ a linear isotherm within the tested concentration range (~10 - 1500 mg L-1), while the results for NH4+ are inconclusive due to potential microbial degradation. Only Na+ showed clear evidence of competitive behaviour, with an order of magnitude decrease in maximum adsorption capacity in the presence of NH4+ (0.22 to 0.02 mol kg-1), which was confirmed by the BTC data where the Na+ retardation coefficient differed between the experiments with different cations. Thus, solute mobility in peatlands is affected by competitive adsorption.

Solute pools in Nikanotee Fen watershed in the Athabasca oil sands region.

Overburden and tailings materials from oil sands production were used as construction materials as part of a novel attempt to create a self-sustaining, peat accumulating fen-upland ecosystem. To evaluate the potential for elemental release from the construction materials, total elemental concentrations in the tailings sand, petroleum coke and peat used to construct a fen ecosystem were determined using microwave-assisted acid digestions and compared to a leaching experiment conducted under environmentally-relevant conditions. A comparison of solid phase to aqueous Na, Ca, S and Mg concentrations showed they were highly leachable in the materials. Given that the concentrations of these elements can affect plant community structure, it is important to understand their leachability and mobility as they migrate between materials used to construct the system. To that end, a mass balance of aqueous Na, Ca, S and Mg was conducted based on leaching experiments and materials analysis coupled with existing data from the constructed system. The data indicate that there is a large pool of leachable Na, Ca, S and Mg in the system, estimated at 27 t of Na, 14 t of Ca, 37.3 t of S and 8.8 t of Mg. Since recharge mainly drives the fen-upland system water regime, and discharge in the fen, evapo-accumulation of these solutes on the surface may occur.

The hydrological functioning of a constructed fen wetland watershed.

Mine reclamation requires the reconstruction of entire landforms and drainage systems. The hydrological regime of reclaimed landscapes will be a manifestation of the processes operating within the individual landforms that comprise it. Hydrology is the most important process regulating wetland function and development, via strong controls on chemical and biotic processes. Accordingly, this research addresses the growing and immediate need to understand the hydrological processes that operate within reconstructed landscapes following resource extraction. In this study, the function of a constructed fen watershed (the Nikanotee Fen watershed) is evaluated for the first two years following construction (2013-2014) and is assessed and discussed within the context of the construction-level design. The system design was capable of sustaining wet conditions within the Nikanotee Fen during the snow-free period in 2013 and 2014, with persistent ponded water in some areas. Evapotranspiration dominated the water fluxes from the system. These losses were partially offset by groundwater discharge from the upland aquifer, which demonstrated strong hydrologic connectivity with the fen in spite of most construction materials having lower than targeted saturated hydraulic conductivities. However, the variable surface infiltration rates and thick placement of a soil-capping layer constrained recharge to the upland aquifer, which remained below designed water contents in much of the upland. These findings indicate that it is possible to engineer the landscape to accommodate the hydrological functions of a fen peatland following surface oil sands extraction. Future research priorities should include understanding the storage and release of water within coarse-grained reclaimed landforms as well as evaluating the relative importance of external water sources and internal water conservation mechanisms for the viability of fen ecosystems over the longer-term.

Landscape analysis of nutrient-enriched margins (lagg) in ombrotrophic peatlands.

Scientific knowledge of the wet zone - the lagg - that tends to form at the edge of ombrotrophic peatlands is surprisingly limited. In this study, we aim to improve the understanding of the ecohydrological functions of this transition by describing the form and abiotic controls of the laggs and margins of bog peatlands. Data collected in wells and piezometers along 10 transects (within 6 bogs), of the New Brunswick Eastern Lowlands are used to analyse the hydraulic and hydrochemical gradients, while airborne LiDAR data provides new insight on the geomorphology and the vegetation patterns of the bog-lagg-mineral transition zone. Based on their geomorphic character, the study transects are placed into 2 categories: confined and unconfined. Laggs of confined transition are found in a topographic depression, between the bog and a mineral slope >1%, while laggs of unconfined transitions are adjacent to a flat (≤1%) or receding mineral slope (sloping away from the lagg). Water level (4 ± 9 cm vs. -3 ± 9 cm), pH (4.8 ± 0.9 vs. 4.2 ± 0.4), electrical conductivity (ECcorr) (105 ± 52 µS cm(-1) vs. 52 ± 28 µS cm(-1)) and peat depth (55 ± 9 cm vs. 30 ± 9 cm) are found to be higher, respectively, for the confined laggs than for the unconfined. Saturated hydraulic conductivity (Ksat) of the lagg's upper peat layer resembles that of bog environments, but quickly reduces with depth, impeding vertical water flow. The greatest abiotic control of the lagg appears to be topography, which affects water flow rates and direction, thus water chemistry, nutrient transport and availability, hence vegetation characteristics. Our results suggest that the features of the transition zone that include the lagg, influence the quantity and variability of water within the adjacent peatland, and should be considered as integral part of the peatland complex.