Experimental and numerical evaluation of single-layer covers placed on acid-generating tailings

This study focuses on the reclamation work being performed on two former acid-generating tailings sites, located in Quebec, Canada. At both sites, the tailings were partially oxidized due to extended exposure, and the pore water is acidic. The reclamation solution applied to control acid mine drainage is a monolayer cover, made of non acid-generating tailings in one case and a till in the other. The goal of this project was to assess the response of the tailings-cover systems under various conditions. Tailings samples were collected in situ and characterized in the laboratory. Large column tests were conducted to evaluate the hydrogeological and geochemical behaviour of the covered tailings following wetting and drying cycles. The instrumented columns were designed to reproduce some of the existing site conditions and provide representative results for longer term analyses. Volumetric water content, suction, and oxygen concentrations were also monitored over time. The experimental data were used to validate different numerical models including those constructed with Vadose/W (GeoSlope Int.). Additional simulations were conducted under field conditions to evaluate the effect of various influence factors such as the depth of the water table, climatic conditions and the thickness of the cover. The combination of experimental and numerical results show how the behaviour and efficiency of a monolayer cover placed over reactive tailings depend on these factors, highlighting the importance of hydrogeological properties and water table depth. In many instances, the cover materials were prone to desaturation, especially when the water table was deeper than about 2 m below the tailings-cover interface. These results tend to indicate that relatively thin monolayer covers would not be able to prevent oxygen ingress under some of the conditions observed in the field. Desaturation of the cover and/or reactive tailings is due to a combination of drainage and evaporation. In such cases, increasing the thickness of the monolayer cover only has a limited effect. A discussion follows on the practical implications of these laboratory experiments and the numerical simulations for field design of reclamation measures at these two tailings impoundments.     

新疆夏季变湿的大气环流异常特征

利用1961-2003年NCEP/NCAR再分析资料和中国气象局整编的新疆夏季(6~8月)月降水量资料,分析了新疆夏季1971-1986年干旱期和1987-2003年湿润期的大气环流变化异常特征.结果表明:在平均环流场上,中亚-巴尔喀什湖槽及上下游地区脊的增强是新疆地区夏季变湿的环流场特征之一,且中亚-巴尔喀什湖槽随高度增加强度明显增强,200 hPa达到最强;源于低纬阿拉伯海向北直至中亚对流层低层偏南的强气流是湿润期环流异常、降水增多的又一特征,也是水汽重要来源之一;中亚上空高空急流轴南压是新疆夏季湿润期与干旱期高空急流最重要的差异特征.     

Fully coupled simulation of a hydraulic fracture interacting with natural fractures with a hybrid discrete‐continuum method

A hybrid discrete‐continuum numerical scheme is developed to study the behavior of a hydraulic fracture crossing natural fractures. The fully coupled hybrid scheme utilizes a discrete element model for an inner domain, within which the hydraulic fracture propagates and interacts with natural fractures. The inner domain is embedded in an outer continuum domain that is implemented to extend the length of the hydraulic fracture and to better approximate the boundary effects. The fracture is identified to propagate initially in the viscosity‐dominated regime, and the numerical scheme is calibrated by using the theoretical plane strain hydraulic fracture solution. The simulation results for orthogonal crossing indicate three fundamental crossing scenarios, which occur for various stress ratios and friction coefficients of the natural fracture: (i) no crossing, that is, the hydraulic fracture is arrested by the natural fracture and makes a T‐shape intersection; (ii) offset crossing, that is, the hydraulic fracture crosses the natural fracture with an offset; and (iii) direct crossing, that is, the hydraulic fracture directly crosses the natural fracture without diversion. Each crossing scenario is associated with a distinct net pressure history. Additionally, the effects of strength contrast and stiffness contrast of rock materials and intersection angle between the hydraulic fracture and the natural fracture are also investigated. The simulations also illustrate that the level of fracturing complexity increases as the number and extent of the natural fractures increase. As a result, we can conclude that complex hydraulic fracture propagation patterns occur because of complicated crossing behavior during the stimulation of naturally fractured reservoirs. Copyright © 2017 John Wiley & Sons, Ltd.     

Assessing the impacts of sea-level rise and precipitation change on the surficial aquifer in the low-lying coastal alluvial plains and barrier islands,east-central Florida (USA)

A three-dimensional variable-density groundwater flow and salinity transport model is implemented using the SEAWAT code to quantify the spatial variation of water-table depth and salinity of the surficial aquifer in Merritt Island and Cape Canaveral Island in east-central Florida (USA) under steady-state 2010 hydrologic and hydrogeologic conditions. The developed model is referred to as the ‘reference’ model and calibrated against field-measured groundwater levels and a map of land use and land cover. Then, five prediction/projection models are developed based on modification of the boundary conditions of the calibrated ‘reference’ model to quantify climate change impacts under various scenarios of sea-level rise and precipitation change projected to 2050. Model results indicate that west Merritt Island will encounter lowland inundation and saltwater intrusion due to its low elevation and flat topography, while climate change impacts on Cape Canaveral Island and east Merritt Island are not significant. The SEAWAT models developed for this study are useful and effective tools for water resources management, land use planning, and climate-change adaptation decision-making in these and other low-lying coastal alluvial plains and barrier island systems.     

Stochastic groundwater modeling of a sedimentary aquifer: evaluation of the impacts of abstraction scenarios under conditions of reduced recharge

A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m³/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.     

Control of sea-water intrusion by salt-water pumping: Coast of Oman

A shallow alluvial coastal aquifer in the Batinah area of Oman, with sea-water intrusion that extends several kilometres inland, has been studied experimentally, analytically and numerically. The water table is proved to have a trough caused by intensive pumping from a fresh groundwater zone and evaporation from the saline phreatic surface. Resistivity traverses perpendicular to the shoreline indicated no fresh groundwater recharge into the sea. Using an analytical Dupuit-Forchheimer model, developed for the plain part of the catchment, explicit expressions for the water table, sharp interface location and stored volume of fresh water are obtained. It is shown that by the pumping of salt water from the intruded part of the aquifer, this intrusion can be mitigated. Different catchment sizes, intensities of fresh groundwater pumping, evaporation rates, water densities, sea level, incident fresh water level in the mountains and hydraulic conductivity are considered. SUTRA code is applied to a hypothetical case of a leaky aquifer with line sinks modeling fresh water withdrawal and evaporation. The numerical code also shows that pumping of saline water can pull the dispersion zone back to the shoreline.     

Analysis of flow behavior in a landfill with cover soil of low hydraulic conductivity

 This paper presents the results of field tests of hydrologic parameters in a landfill and the results of numerical simulation to find the efficiency of the pumping method to reduce leachate levels in the landfill. The field hydraulic conductivity and storativity of waste and buried cover soils in the landfill are measured by pumping and slug tests. The hydrologic condition inside the landfill is first calibrated using the drawdown-time curve obtained from the pumping test, and the flow behavior of leachate during pumping in the landfill, when various layers of waste and buried cover soil exist, is analyzed through three-dimensional numerical simulation of flow. The results of the field investigation show that the buried cover soil of low hydraulic conductivity forms an impermeable layer preventing the downward flow of leachate and upward flow of landfill gas. The hydraulic conductivities of the pumping test and slug tests were quite close on the same order of magnitude. It was also possible to match the drawdown-time data of the field tests with those of the model using input data close to the hydrologic property obtained from the field tests. The numerical flow analysis showed that pumping was possible up to 120 tons/day for a single well without a drain, while the pumping rate could be increased to 300 tons/day for the same well with the drain. From the vertical section of the flow vector with a horizontal drain, the barrier role of buried cover soil is identified, which was proposed by examining the water contents of the disposed cover soil and waste in the field. Received: 15 May 1998 · Accepted: 4 January 1999     

Arsenic mobility in weathered gold mine tailings under a low-organic soil cover

Historical gold mining operations in Nova Scotia, Canada, resulted in numerous deposits of publicly accessible, arsenic (As)-rich mine waste that has weathered in situ for 75–150 years, resulting in a wide range of As-bearing secondary minerals. The geochemical heterogeneity of this mine waste creates a challenge for identifying a single remediation approach that will limit As mobility. A 30-cm-thick, low-organic content soil cover was evaluated in a laboratory leaching experiment where, to simulate natural conditions, the equivalent of 2 years of synthetic rainwater was leached through each column and two dry seasons were incorporated into the leaching protocol. Each column was a stratigraphic representation of the four major tailings types found at the historical Montague and Goldenville gold mine districts: hardpan tailings, oxic tailings, wetland tailings, and high Ca tailings. Hardpan tailings released acidic, As-rich waters (max 12 mg/L) under the soil cover but this acidity was buffered by surrounding oxic tailings. Leachate from the oxic tailings was circumneutral, with average As concentrations between 4.4 and 9.7 mg/L throughout the experiment. The presence of carbonates in the high Ca tailings resulted in near-neutral to weakly alkaline leachate pH values and average As concentrations between 2.1 and 6.1 mg/L. Oxidation of sulfides in the wetland tailings led to acidic leachate over time and a decrease in As concentrations to values that were generally less than 1 mg/L. This study shows that the use of a low-organic content soil cover does not create reducing conditions that would destabilize oxidized, As-bearing secondary phases in these tailings. However, oxygen penetration through the cover during dry seasons would continue to release As to tailings pore waters via sulfide oxidation reactions.     

Modeling of seawater intrusion in a coastal aquifer of Karaburun Peninsula,western Turkey

Seawater intrusion is a major problem to freshwater resources especially in coastal areas where fresh groundwater is surrounded and could be easily influenced by seawater. This study presents the development of a conceptual and numerical model for the coastal aquifer of Karareis region (Karaburun Peninsula) in the western part of Turkey. The study also presents the interpretation and the analysis of the time series data of groundwater levels recorded by data loggers. The SEAWAT model is used in this study to solve the density-dependent flow field and seawater intrusion in the coastal aquifer that is under excessive pumping particularly during summer months. The model was calibrated using the average values of a 1-year dataset and further verified by the average values of another year. Five potential scenarios were analyzed to understand the effects of pumping and climate change on groundwater levels and the extent of seawater intrusion in the next 10 years. The result of the analysis demonstrated high levels of electrical conductivity and chloride along the coastal part of the study area. As a result of the numerical model, seawater intrusion is simulated to move about 420 m toward the land in the next 10 years under “increased pumping” scenario, while a slight change in water level and TDS concentrations was observed in “climate change” scenario. Results also revealed that a reduction in the pumping rate from Karareis wells will be necessary to protect fresh groundwater from contamination by seawater.     

Energy‐work‐based numerical manifold seepage analysis with an efficient scheme to locate the phreatic surface

A major challenge in seepage analysis is to locate the phreatic surface in an unconfined aquifer. The phreatic surface is unknown and assumed as a discontinuity separating the seepage domain into dry and wet parts, thus should be determined iteratively with special schemes. In this study, we systematically developed a new numerical manifold method (NMM) model for unconfined seepage analysis. The NMM is a general numerical method for modeling continuous and discontinuous deformation in a unified mathematical form. The novelty of our NMM model is rooted in the NMM two‐cover‐mesh system: the mathematical covers are fixed and the physical covers are adjusted with iterations to account for the discontinuity feature of the phreatic surface. We developed an energy‐work seepage model, which accommodates flexible approaches for boundary conditions and provides a form consistent with that in mechanical analysis with clarified physical meaning of the potential energy. In the framework of this energy‐work seepage model, we proposed a physical concept model (a pipe model) for constructing the penalty function used in the penalty method to uniformly deal with Dirichlet, Neumann, and material boundaries. The new NMM model was applied to study four example problems of unconfined seepage with varying geometric shape, boundary conditions, and material domains. The comparison of our simulation results to those of existing numerical models for these examples indicates that our NMM model can achieve a high accuracy and faster convergence speed with relatively coarse meshes. This NMM seepage model will be a key component of our future coupled hydro‐mechanical NMM model. Copyright © 2014 John Wiley & Sons, Ltd.     

Effects of climate change on the groundwater system in the Grote-Nete catchment,Belgium

The effects of climate change on the groundwater systems in the Grote-Nete catchment, Belgium, covering an area of 525 km², is modeled using wet (greenhouse), cold or NATCC (North Atlantic Thermohaline Circulation Change) and dry climate scenarios. Low, central and high estimates of temperature changes are adopted for wet scenarios. Seasonal and annual water balance components including groundwater recharge are simulated using the WetSpass model, while mean annual groundwater elevations and discharge are simulated with a steady-state MODFLOW groundwater model. WetSpass results for the wet scenarios show that wet winters and drier summers are expected relative to the present situation. MODFLOW results for wet high scenario show groundwater levels increase by as much as 79 cm, which could affect the distribution and species richness of meadows. Results obtained for cold scenarios depict drier winters and wetter summers relative to the present. The dry scenarios predict dry conditions for the whole year. There is no recharge during the summer, which is mainly attributed to high evapotranspiration rates by forests and low precipitation. Average annual groundwater levels drop by 0.5 m, with maximum of 3.1 m on the eastern part of the Campine Plateau. This could endanger aquatic ecosystem, shrubs, and crop production.     

Reconstructing Phreatic Palaeogroundwater Levels in a Geoarchaeological Context: A Case Study in Flanders,Belgium

The complex debate on prehistoric settlement decisions is no longer tackled from a purely archaeological perspective but from a more landscape‐oriented manner combined with archaeological evidence. Therefore, reconstruction of several components of the former landscape is needed. Here, we focus on the reconstruction of the groundwater table based on modeling. The depth of the phreatic aquifer influences, for example, soil formation processes and vegetation type. Furthermore, it directly influences settlement by the wetness of a site. Palaeogroundwater modeling of the phreatic aquifer was carried out to produce a series of full‐coverage maps of the mean water table depth between 12.7 ka and the middle of the 20th century (1953) in Flanders, Belgium. The research focuses on the reconstruction of the input data and boundary conditions of the model and the model calibration. The model was calibrated for the 1924–1953 time period using drainage class maps. Archaeological site data and podzol occurrence data act as proxies for local drainage conditions over periods in the past. They also served as a control on the simulated phreatic palaeogroundwater levels. Model quality testing on an independent validation data set showed that the model predicts phreatic water table levels at the time of soil mapping well (mean error of 1.8 cm; root mean square error of 65.6 cm). Simulated hydrological conditions were in agreement with the occurrence of archaeological sites of Mesolithic to Roman age at 96% of the validation locations, and also with the occurrence of well‐drained podzols at 97% of the validation locations.     

Hydraulic Conductivity of Fly Ash-Amended Mine Tailings

The objective of this study was to evaluate the effect of fly ash addition on hydraulic conductivity (k) of mine tailings. Mine tailings used in this study were categorized as synthetic tailings and natural tailings; two synthetic tailings were developed via blending commercially-available soils and natural tailings were collected from a garnet mine located in the U.S. Two fly ashes were used that had sufficient calcium oxide (CaO) content (17 and 18.9 %) to generate pozzolanic activity. Hydraulic conductivity was measured on pure tailings and fly ash-amended tailings in flexible-wall permeameters. Fly ash was added to mine tailings to constitute 10 % dry mass of the mixture, and specimens were cured for 7 and 28 days. The influence of fly ash-amendment on k of mine tailings was attributed to (1) molding water content and (2) plasticity of the mine tailings. Tailings that classified as low-plasticity silts with clay contents less than 15 % exhibited a decrease in k when amended with fly ash and prepared wet of optimum water content (w _opt ). Tailings that classified as low-plasticity clay exhibited a one-order magnitude increase in k with addition of fly ash for materials prepared dry or near w _opt . The decrease in k for silty tailings was attributed to formation of cementitious bonds that obstructed flow paths, whereas the increase in k for clayey tailings was attributed to agglomeration of clay particles and an overall increase in average pore size. The results also indicated that the effect of curing time on k is more pronounced during the early stages of curing (≤7 days), as there was negligible difference between k for 7 and 28-days cured specimens.     

Numerical modeling the role of rubber dams on groundwater recharge and phreatic evaporation loss in riparian zones

Rubber dams have been widely built for their advantages in increase of flooding resources utilization in the north arid and semiarid plain regions of China. Rise in river water stage by the dams, particularly during the drought periods, increases lateral seepage of river water into groundwater, and thus groundwater table and phreatic evaporation loss in the riparian zones. In this study, a riparian area of the Baihe River in Nanyang of Henan Province, China was selected for investigation of influences of the river dams on the groundwater recharge and evaporation loss. A hydraulic model, HEC-RAS, was used for simulation of the river stage variations along the Baihe River, and a numerical groundwater model, MODFLOW, was applied for simulation of groundwater dynamics and estimation of river flow seepage into aquifer and evaporation loss. The results show that the dams increase river stages of 2–3 m during January 2000–December 2002. The increase in the captured groundwater recharge was 7.15–34.06 million m³/a and the increased phreatic evaporation loss occupies 10% of the increased recharge when four rubber dams were built.     

The mobility of radium-226 and trace metals in pre-oxidized subaqueous uranium mill tailings

The exchange of ²²⁶Ra and trace metals across the tailings-water interface and the mechanisms governing their mobility were assessed via sub-centimetre resolution profiling of dissolved constituents across the tailings–water interface in Cell 14 of the Quirke Waste Management Area at Rio Algom's Quirke Mine, near Elliot Lake, Ontario, Canada. Shallow zones (<1.5 m water depth) are characterized by sparse filamentous vegetation, well-mixed water columns and fully oxygenated bottom waters. Profiles of dissolved O₂, Fe and Mn indicate that the tailings deposits in these areas are sub-oxic below tailings depths of ∼3 cm. These zones exhibit minor remobilization of Ra in the upper 5 cm of the tailings deposit; ²²⁶Ra fluxes at these sites are relatively small, and contribute negligibly to the water column activity of ²²⁶Ra. The shallow areas also exhibit minor remobilization of Ni, As, Mo and U. The release of these elements to the water cover is, however, limited by scavenging mechanisms in the interfacial oxic horizons. The presence of thick vegetation (Chara sp.) in the deeper areas (>2 m water depth) fosters stagnant bottom waters and permits the development of anoxia above the benthic boundary. These anoxic tailings are characterized by substantial remobilization of ²²⁶Ra, resulting in a relatively large flux of ²²⁶Ra from the tailings to the water column. The strong correlation between the porewater profiles of ²²⁶Ra and Ba (r²=0.99), as well as solubility calculations, indicate that the mobility of Ra is controlled by saturation with respect to a poorly ordered and/or impure barite phase [(Ra,Ba)SO₄]. In the anoxic zones, severe undersaturation with respect to barite is sustained by microbial SO₄ reduction. Flux calculations suggest that the increase in ²²⁶Ra activity in the water cover since 1995 (from <0.5 to 2.5 Bq l⁻¹) can be attributed to an increase in the spatial distribution of anoxic bottom waters caused by increased density of benthic flora. The anoxic, vegetated areas also exhibit minor remobilization with respect to dissolved As, Ni and Zn. The removal of trace metals in the anoxic bottom waters appears to be limited by the availability of free sulphide. Collectively, the data demonstrate that while the water cover over the U mill tailings minimizes sulphide oxidation and metal mobility, anoxic conditions which have developed in deeper areas have led to increased mobility of ²²⁶Ra.     

Assessing the impact of future climate change on groundwater recharge in Galicia-Costa, Spain

Climate change can impact the hydrological processes of a watershed and may result in problems with future water supply for large sections of the population. Results from the FP5 PRUDENCE project suggest significant changes in temperature and precipitation over Europe. In this study, the Soil and Water Assessment Tool (SWAT) model was used to assess the potential impacts of climate change on groundwater recharge in the hydrological district of Galicia-Costa, Spain. Climate projections from two general circulation models and eight different regional climate models were used for the assessment and two climate-change scenarios were evaluated. Calibration and validation of the model were performed using a daily time-step in four representative catchments in the district. The effects on modeled mean annual groundwater recharge are small, partly due to the greater stomatal efficiency of plants in response to increased CO₂ concentration. However, climate change strongly influences the temporal variability of modeled groundwater recharge. Recharge may concentrate in the winter season and dramatically decrease in the summer–autumn season. As a result, the dry-season duration may be increased on average by almost 30 % for the A2 emission scenario, exacerbating the current problems in water supply.     

Assessing the impacts of climate and land use change on streamflow,water quality and suspended sediment in the Kor River Basin,Southwest of Iran

This research addressed the separate and combined impacts of climate and land use change on streamflow, suspended sediment and water quality in the Kor River Basin, Southwest of Iran, using (BASINS–WinHSPF) model. The model was calibrated and validated for hydrology, sediment and water quality for the period 2003–2012. The model was run under two climate changes, two land use changes and four combined change scenarios for near-future period (2020–2049). The results revealed that projected climate change impacts include an increase in streamflow (maximum increases of 52% under RCP 2.6 in December and 170% under RCP 8.5). Projected sediment concentrations under climate change scenarios showed a monthly average decrease of 10%. For land use change scenarios, agricultural development scenario indicated an opposite direction of changes in orthophosphate (increases in all months with an average increase of 6% under agricultural development scenario), leading to the conclusion that land use change is the dominant factor in nutrient concentration changes. Combined impacts results indicated that streamflows in late fall and winter months increased while in summer and early fall decreased. Suspended sediment and orthophosphate concentrations were decreased in all months except for increases in suspended sediment concentrations in September and October and orthophosphate concentrations in late winter and early spring due to the impact of land use change scenarios.     

某磷石膏尾矿库堆积坝渗透稳定性分析

随着尾矿库堆积坝高度的增加,库内的地下水渗流场将发生改变,当尾矿堆积坝处的实际水力梯度大于临界水力梯度时,尾矿库堆积坝将产生渗透变形或失稳破坏。为了预测某已建磷石膏尾矿库堆积坝的渗透稳定性,本文运用地下水位计水位恢复试验去确定磷石膏尾矿渗透系数并计算水力梯度,同时采用Geo-studio软件中的SEEP/W模块分析磷石膏堆积坝在不同工况下的渗透稳定性。计算和分析结果表明当堆积坝处于945m堆积高程,在无干滩和100m干滩情况下,实际水力梯度分别为0.784和0.583;当堆积坝处于960m堆积高程,无干滩情况下,实际水力梯度为0.692,在以上两种情况下,堆积坝渗出点高程高于初始坝高程,堆积坝将出现流土变形,坝体处于不稳定状态。当堆积坝处于960m堆积高程100m干滩时,堆积坝上无渗出点,堆积坝处于稳定状态。分析研究结果对该尾矿库的安全运行和维护管理提供技术依据,同时为同类型尾矿库堆积坝渗透稳定性评价提供参考。     

Mobilization of arsenic from mine tailings through reductive dissolution of goethite influenced by organic cover

Mine tailings at the former Delnite gold mine in northern Ontario were characterized to assess the impact of a biosolids cover on the stability of As species and evaluate options for long-term management of the tailings. Arsenic concentrations in the tailings range from 0.15 to 0.36 wt% distributed among goethite, pyrite and arsenopyrite. Pyrite and arsenopyrite occur as small and liberated particles that are enveloped by goethite in the uncovered tailings and the deeper portions of the biosolids-covered tailings. Sulfide particles in the shallower portions of the biosolids-covered tailings are free of goethite rims. Arsenic occurs predominantly as As⁵⁺ with minor amount of As¹⁻ in the uncovered tailings. Coinciding with the disappearence of goethite rims on sulfide particles, the biosolids-covered tailings have As³⁺ species gradually increasing in proportion towards the cover. Leaching tests indicated that the As concentrations in the leachate gradually increase from less than 0.085 to 13 mg/L and Fe from 28 to 179 mg/L towards the biosolids cover. These are in sheer contrast to the leachate concentrations of less than 0.085 mg/L As and 24–64 mg/L Fe obtained from the uncovered tailings confirming the role of biosolids-influenced reduction and mobilization of As in the form of As³⁺ species. The evidence suggests that reductive dissolution of goethite influenced by the biosolids-cover caused the mobilization of As as As³⁺ species.     

Microcosm assessment of the biogeochemical development of sulfur and oxygen in oil sands fluid fine tailings

Bitumen recovery from Alberta oil sands generates fluid fine tailings, which are retained in tailings ponds where solids settle and release process water. The recovered water is recycled for bitumen extraction, while the resulting tailings are incorporated into various landforms for reclamation, with one option being conversion of tailings basins to viable end pit lakes. Tailings ponds commonly host diverse microbial communities, including SO₄-reducing prokaryotes. The highly reducing nature of the hydrogen sulfide produced by these prokaryotes may impact the biogeochemical cycling of key nutrients. However, the behavioral dynamics of hydrogen sulfide production in ponds containing fluid fine tailings remain to be clearly explained. In this study, microcosms are used as analogues of the sediment–water interface of a tailings pond undergoing reclamation to determine sulfide generation patterns and the behavior of O₂. In the microcosms, hydrogen sulfide fluxes correlated positively with biotic activity, reaching levels of over 2 × 10³ nmol cm⁻² s⁻¹, leading to Fe sulfide formation. Depth-related hydrogen sulfide profiles in the microcosms were comparable to those encountered in situ, in Syncrude’s West In-Pit, an active tailing pond. Oxygen diffusion across the fluid fine tailing sediment–water interface was controlled to different degrees by both biotic and abiotic processes. The results have implications for quantitatively estimating the impact of hydrogen sulfide production, O₂ availability, and biogeochemical cycling of key nutrients important for the success of life in fluid fine tailings-affected ecosystems. This paper shows that this production of hydrogen sulfide may be a self-limiting process, which will begin to decrease after a period of time.