Fractured Bedrock Research Articles

Hydrogeochemical and isotopic evaluation of VOC commingled plumes in a weathered fractured bedrock aquifer treated with thermal and bioremediation

Chlorinated ethanes and ethenes isotopic analyses in groundwater and hydrogeochemical results from a former industrial area in Sao Paulo (Brazil) were used to confirm the existence and allow further characterization of source areas and their commingled plumes, both before and after thermal and bioremediation treatments. Prior to full scale remediation, a recently identified off-site source area with unknown history and limited access for further intrusive works presented lower δ13C values (−6.5‰ to −1.8‰ for 1,2-DCA) than the downgradient on-site source area (+8.6‰ to +20.0‰). Intermediate δ13C values for 1,2-DCA were identified further downgradient from the sources, within commingled plumes patterns. The isotope and concentration results show the typical degradation patterns associated with biotic reductive dechlorination for chlorinated ethenes and dihaloelimination for 1,2-DCA. Results following remediation treatments show further levels of isotopic enrichment, for chlorinated ethenes and chlorinated ethanes in the tropically weathered and deeper fractured bedrock (gneisses) groundwater. Hydrogeochemical results, isotopic mass balance and Carbon-Chlorine isotope slopes data are coherent with remediation treatment and a complex commingled plume setting. The results of this study confirmed the Temporal Conceptual Model proposed by Hart et al. (2021) and identified the need for further studies to evaluate isotopic dynamics under thermal remediation, including thermal-induced hydrolysis processes.

Dirt cracking as rock fracture-wedging process in the Mediterranean climate of Victoria, British Columbia, Canada

Dirt cracking physically wedges open fractures in desert bedrock via the synergistic processes of carbonate precipitation and the wetting and drying of clays in fracture fill. We use back-scattered and high-resolution electron microscopy along with 87Sr/86Sr analyses to find that dirt cracking also occurs in the hypermaritime Mediterranean climate of Victoria, British Columbia, Canada. 87Sr/86Sr analyses of calcium carbonate from four different sampling locations do not reflect an ocean source, but 87Sr/86Sr ratios correspond with the known composition of the sampled rocks. Carbonates derived from decay of host-rock minerals precipitate in narrow fractures, widening them. In addition to calcium carbonate, iron carbonate, barium carbonate, and calcium phosphate also precipitate in fractures. Dissolution of silicates, in association with carbonate precipitation, aids in the further penetration of carbonates. Fines falling into fissures includes smectites that undergo wetting (expansion) and drying (contraction) that further promotes fracture widening. Lead contamination, probably from mid-20th century automobile emissions, occurs from heavy metal scavenging by iron oxides. As studies of dirt cracking are few in number, and since the carbonate precipitation process has been observed in Antarctica, monsoonal Asia, and now in a Mediterranean climate, it is possible that dirt cracking is a common rock-decay process in biogeochemical settings dry enough for carbonate precipitation in rock fissures.

Rhizosphere water content drives hydraulic redistribution: Implications of pore-scale heterogeneity to modeling diurnal transpiration in water-limited ecosystems

Trees typically survive prolonged droughts by absorbing water from deeper layers. Where soils are shallow, roots may be extract water from the underlying fractured bedrocks. In dry seasons, surface-soil moisture dynamics reflect hydraulic redistribution (HR). HR is usually estimated based on the gradient of mean, or bulk, soil water potential among layers in the rooting zone (HRB). This approach neglects the potential effect of spatial heterogeneity of water content at the millimeter scale between the rhizosphere and bulk soil. We proposed to account for the rhizosphere water balance, estimating HR to the rhizosphere (HRR) of the dry surface soil from the underlying fractured rock. The model was evaluated using a 15-year dataset collected in Sardinia. When the typical approach, based on moisture gradients among bulk soil layers, was used for estimating HRB, tree transpiration was underpredicted in all seasons, especially in spring and summer. Forcing the model with measured tree transpiration, HRB decreased during spring and summer, while the contribution of the underlying rock layer to tree transpiration was threefold that estimated using HRR-based model. The average water content of the bulk surface soil layer was very low, reaching 0.06 in the driest summers while showing little diurnal dynamics; however, concentrating water in roughly estimated rhizosphere volume, produced rhizosphere water content appreciably higher (≈0.16), and much more dynamic. Predicted HRR dominated evapotranspiration (60% - 65%) in dry springs and summers reaching 80% of tree transpiration. Most importantly, the proposed rhizosphere-HR model correctly predicts the diurnal dynamics of tree transpiration year-round, and the grass transpiration in its active spring period. Eco-hydrological models operating at sub-daily scale should consider partitioning the soil to rhizosphere volume, thus allowing both diagnostic and prognostic estimates of diurnal biosphere-atmosphere mass and energy exchanges.

Water Flows in Rock Wall Permafrost: A Numerical Approach Coupling Hydrological and Thermal Processes

AbstractRock wall permafrost is highly sensitive to atmospheric warming. Its degradation is regarded as a primary factor accounting for the recent increase in periglacial rock slope failures. Investigations of rock wall permafrost have so far ignored hydrological processes acting in bedrock fractures. We propose a 2‐D coupled thermal and hydrological modeling approach applied to alpine rock wall permafrost. Simulations are carried out on the scale of the mountain flank for four case studies with variable saturation levels and hydrological forcings. Results show that cold water infiltration from the summit leads to a deeper permafrost body than with no water flows. Ice‐filled fractures first delay permafrost thawing due to latent heat consumption but then accelerate it when the ice starts to melt. Ice may thus subsist in fractures surrounded by thawed bedrock while thawing corridors may form in fractures embedded in frozen bedrock. As a result, temperature gradients are made steeper. When connected fractures thaw, bottom‐up permafrost degradation can occur through upward propagation of thawing wedges delineated by these fractures. High hydraulic head values are associated with the formation of perched water tables above or within the impermeable permafrost body, and correspond to hydrostatic pressures that could reach critical values in terms of rock wall stability. The proposed model setting is more appropriate for investigating long‐term permafrost changes, rock wall permafrost dynamics and hydrogeology across entire mountain flanks, but could be applied on a finer spatial scale to investigate short‐term permafrost response to climate signals and rock wall destabilization.

Interpretation of a network-scale tracer experiment in fractured rock conducted using open wells

The presence of fractures in bedrock allows for rapid aqueous contaminant transport through complex pathways and for diffusion of solutes between the fractures and the matrix. To better understand transport in these settings, tracer experiments are a commonly used tool. The need for expensive multi-level wells to obtain depth-specific concentrations, however, significantly limits the cost efficiency. The primary objective of this study is to develop a method whereby a discrete fracture network approach can be used to simulate the results of a divergent tracer experiment conducted using open observation boreholes in a well-characterized dolostone over distances of 55 m to 242 m. The experiment was conducted using a fluorescent tracer which allowed for continuous concentration measurement with depth in each open observation well. Two numerical models were employed in the interpretation of the experiment. The first was a 1-D finite difference model focused on flow and transport in the observation wells and the second was a 3-D control-volume finite element model capable of simulating the entire fracture network. Through fitting the experimental data to simulations, the most important fractures for transport in the system were identified. The number of fractures that participated in transport was few relative to the number of fractures observed in core and in constant head test results. Heterogeneous distribution of the fracture apertures was determined to be the likely cause of the highly tortuous transport observed at the site. This study demonstrates that tracer experiments conducted using open observation boreholes and a downhole fluorometer can improve our understanding of large-scale transport in fractured rock, especially when analysed with multiple models, and compared to other measured properties such as matrix porosity, hydraulic aperture, and fracture orientation.

Reconstructing the Gorte and Spiaz de Navesele Landslides, NE of Lake Garda, Trentino Dolomites (Italy)

We applied a multi-method approach to reconstruct the Gorte rock avalanche (85–95 Mm3) located at the northeastern end of Lake Garda. The combination of field mapping, characterization of bedrock discontinuities, Dan3D-Flex runout modeling and dating of boulders with cosmogenic 36Cl supports the conclusion that the deposits stem from a single rock avalanche at 6.1 ± 0.8 ka. The Gorte event may have triggered the Spiaz de Navesele–Salto della Capra landslide (3.2 Mm3), whose deposits cover the southern end of the Gorte deposits. First-order controls on detachment were the NNE–SSW- and WNW–ESE-oriented fractures in the limestone bedrock, related to the Giudicarie and Schio-Vicenza fault systems, respectively. Dan3D-Flex runout modeling sufficiently reproduced the Gorte rock avalanche, which involved detachment and sliding of a quasi-intact block, likely along marly interbeds, followed by rapid disintegration. The frictional rheology in the source area and the turbulent frictional rheology (Voellmy) in the remaining part best replicate the observed deposit extent and thickness. Heavy precipitation that occurred at that time may have contributed to failure at Gorte. Nonetheless, its timing overlaps with the nearby (<15 km) Dosso Gardene (6630–6290 cal BP) and Marocca Principale (5.3 ± 0.9 ka) landslides, making a seismic trigger plausible.

Matrix Diffusion as a Mechanism Contributing to Fractal Stream Chemistry and Long‐Tailed Transit Time Distributions

AbstractSolute transit or travel time distributions (TTDs) in catchments are relevant to both hydrochemical response and inference of hydrologic mechanisms. Long‐tailed TTDs and fractal scaling behavior of stream concentration power spectra (∼1/frequency, or 1/frequency to a power <2) are widely observed in catchment studies. In several catchments, a significant fraction of streamflow is derived from groundwater in shallow fractured bedrock, where matrix diffusion significantly influences solute transport. I present frequency and time domain theoretical analyses of solute transport to quantify the influence of matrix diffusion on fractal scaling and long‐tailed TTDs. The theoretical concentration power spectra exhibit fractal scaling, and the corresponding TTDs resemble a gamma distribution. The tails of the TTDs are influenced by accessible matrix width, exhibiting a sustained power‐law (rather than exponential) decline for large matrix widths. Application to an experimental catchment shows that theoretical spectra match previously reported power spectral estimates derived from concentration measurements.

Occurrence, Geochemistry and Speciation of Elevated Arsenic Concentrations in a Fractured Bedrock Aquifer System

The presence of elevated arsenic concentrations (≥ 10 µg L−1) in groundwaters has been widely reported in areas of South-East Asia with recent studies showing its detection in fractured bedrock aquifers is occurring mainly in regions of north-eastern USA. However, data within Europe remain limited; therefore, the objective of this work was to understand the geochemical mobilisation mechanism of arsenic in this geologic setting using a study site in Ireland as a case study. Physicochemical (pH, Eh, d-O2), trace metals, major ion and arsenic speciation samples were collected and analysed using a variety of field and laboratory-based techniques and evaluated using statistical analysis. Groundwaters containing elevated dissolved arsenic concentrations (up to 73.95 µg L−1) were characterised as oxic-alkali groundwaters with the co-occurrence of other oxyanions (including Mo, Se, Sb and U), low dissolved concentrations of Fe and Mn, and low Na/Ca ratios indicated that arsenic was mobilised through alkali desorption of Fe oxyhydroxides. Arsenic speciation using a solid-phase extraction methodology (n = 20) showed that the dominant species of arsenic was arsenate, with pH being a major controlling factor. The expected source of arsenic is sulphide minerals within fractures of the bedrock aquifer with transportation of arsenic and other oxyanion forming elements facilitated by secondary Fe mineral phases. However, the presence of methylarsenical compounds detected in groundwaters illustrates that microbially mediated mobilisation processes may also be (co)-occurring. This study gives insight into the geochemistry of arsenic mobilisation that can be used to further guide research needs in this area for the protection of groundwater resources.

Probabilistic Evaluation of Geoscientific Hypotheses With Geophysical Data: Application to Electrical Resistivity Imaging of a Fractured Bedrock Zone

AbstractAs climate changes and populations grow, groundwater sustainability is becoming increasingly important. Hydrogeologic models, which are based on a conceptual understanding of the subsurface, are crucial tools for informing decisions. Conceptual models of the subsurface incorporate knowledge of geological processes, and, frequently, observations from geophysical data into a common subsurface parameterization where the parameters may still be uncertain. Many methods exist to test how different conceptual subsurface parameterizations compare to geophysical data, but a frequent problem in hydrogeologic model development occurs when multiple geological phenomena could explain a single subsurface parameterization. In this work, we present a framework for testing geological hypotheses in conditions where a geological feature is observed in geophysical data, but its physical characteristics are uncertain. The framework uses Popper‐Bayes methods developed in previous studies, and is applied to study a fractured bedrock zone in a mountainous watershed in southwest Colorado. First, we propose six hypotheses based on the geological history of the watershed. Then, using the proposed Popper‐Bayes approach, we demonstrate that two of the hypotheses are inconsistent with the electrical resistivity tomography data. Finally, we discuss the importance of the prior model, and in what other scenarios the framework can be applied to.

Shear joints and its relations with subsurface structures in Batu Melintang, Jeli, Kelantan, Malaysia

Shear joint is the common rock deformation structures formed in Batu Melintang, Jeli due to its location within Bentong-Raub Suture zone. The structural analysis of shear joint can give information about the direction of maximum and minimum stress exerted on a rock while undergoing deformation as the effect of stress fields in the study area. The subsurface structural analysis is done by using the geophysical resistivity method. It displays the subsurface structure in the area for confirmation of the structure found on the surface whether it is highly fractured, moderately fractured or low fractured. The research area was divided into six grids for systematic field measurement. The shear joints orientation were taken while conducting geological mapping and recorded using rose diagram analysis; while the geophysical resistivity method was carried out with a varied length of survey lines set at 200/100m and 1.25/2.5/5m electrode spacing. The subsurface depth of penetration for each survey line is varied, ranging from 0m to 50m. The data is processed in RES2DINV software to obtain the pseudosection profile of the subsurface. The study area principal stress was identified; the maximum stress force ?1 was directed in NW-SE in direction of S107°E and N287°W, while minimum stress ?3 was directed in NE-SW in direction of N17°E and S197°W. The pseudosection subsurface image also displayed a correlation between surface shear joint structures and subsurface structures. The subsurface investigation; according to the pseudosection found that the study area consists of highly fractured structure displayed as several weak zones and fractures of bedrock.

Flow directions of stream‐groundwater exchange in a headwater catchment during the hydrologic year

AbstractUnderstanding near‐stream groundwater dynamics and flow directions is important for predicting hillslope‐stream connectivity, streamflow generation, and hydrologic controls of streamwater quality. To determine the drivers of groundwater flow in the stream corridor (i.e., the stream channel and the adjacent groundwater in footslopes and riparian areas), we observed the water levels of 36 wells and 7 piezometers along a headwater stream section over a period of 18 months. Groundwater dynamics during events were controlled by the initial position of the groundwater table relative to the subsurface structure. The near‐stream groundwater table displayed a fast and pronounced response to precipitation events when lying in fractured bedrock with low storage capacity, and responded less frequently and in a less pronounced way when lying in upper layers with high storage capacity. Precipitation depth, intensity, regolith thickness above the fractured bedrock, and proximity to and elevation above the stream channel also had an effect on the groundwater dynamics, which varied with hydrologic conditions. Our high‐frequency and spatially dense measurements highlight the competing influence of groundwater inflow from upslope locations, streamwater level and bedrock properties on the spatiotemporal dynamics of flowpaths in the stream corridor. Near‐stream groundwater pointed uniformly towards the stream channel when the stream corridor was hydrologically connected to upslope groundwater. However, local interruptions of the water inflow from upslope locations caused flow reversals towards the footslopes. The direction of near‐stream groundwater followed the local fractured bedrock topography during dry hydrologic conditions on a few occasions after events. The outcomes of this research contribute to a better understanding of the drivers controlling spatiotemporal changes in near‐stream groundwater dynamics and flow directions in multiple wetness states of the stream corridor.

The distribution of glacial meltwater routes and associated murtoo fields in Finland

High-resolution LiDAR (Light detection and Ranging) -based digital elevation models (DEM) have greatly improved the mapping of glacial landforms and revealed new ones such as murtoos. Murtoos have extensive diversity in form, relief and size, and they often appear along meltwater routes. However, not all meltwater routes in the recently glaciated terrains include murtoos. We mapped different types of subglacial meltwater routes and the related distribution of murtoos in the Finnish part of Fennoscandian Ice Sheet (FIS). Subglacial meltwater routes represent a previously unknown extension of the subglacial hydrological system that supplements esker networks. Murtoo deposition along the routes is dictated by the marked concentration and routing of subglacial meltwater in a high-pressure environment outside the subglacial tunnel flow and channelized drainage zone. The main environments of murtoo route genesis include the margins of glacial lineation fields, lateral shear margins of ice streams or in between ice-flow sectors or corridors, confluence zones of ice stream onset areas, lee-sides of bedrock protrusions or thresholds, bedrock fracture valleys and potential subglacial lake inputs and outputs. This study adds to the knowledge of the development of subglacial drainage and emphasizes the role of murtoos as the potential missing link between the channelized and distributed subglacial drainage.

How Do Fractures Influence Hyporheic Exchange in Sedimentary Rock Riverbeds?

AbstractBedrock rivers represent a hydrogeological environment in which surface water flows along an exposed bedrock surface. Studies of hyporheic exchange have exclusively involved rivers composed of unconsolidated fluvial sediments, leaving a critical knowledge gap. This study evaluates the conditions that could support bedform‐scale hyporheic exchange within a fractured sedimentary bedrock river based on field data collected near Guelph, Canada. Hyporheic exchange at the bedform‐scale was evaluated by numerically modeling the migration of a conservative solute tracer through a bedrock riverbed within a two‐dimensional vertical cross‐section along the flow direction. A stochastic discrete fracture‐matrix framework was developed to represent measured subsurface fractured bedrock properties, producing a probabilistic distribution of potential hyporheic exchange pathways. Flow and transport model results indicate that: (a) bedform‐scale hyporheic exchange within a bedrock river exists with high fluid flow velocities in fractures, yet long solute residence times due to diffusion across the fracture‐matrix boundary; (b) the coincidence of fractures and hydrodynamic head gradients across the riverbed controls the spatial extent of bedform‐scale hyporheic exchange; and (c) the potential variability in hyporheic exchange residence time is large (i.e., tens of years) due to the inherent variability in facture network properties. Our field‐based numerical study indicates that the average (median) residence time may not be a good proxy for the potential natural attenuation capacity of a fractured sedimentary bedrock riverbed and that hyporheic exchange has the potential to emplace surface water contaminants within the fractured porous rock matrix that could become a long‐term source of trace contaminants.

Characterization of groundwater circulations in a headwater catchment from an analysis of chemical concentrations, Sr-Nd-U isotope ratios, and CFC, SF6 gas tracers (Strengbach CZO, France)

In order to characterize the chemical composition and the age of the water circulating in the critical zone of the Strengbach catchment (Vosges mountains, France), water samples from springs, 10–15 m deep piezometers and 50–120 m deep boreholes were collected and analyzed in elementary concentrations, in Sr, Nd, and U isotopic ratios, and in chlorofluorocarbon (CFC) and sulfur hexafluoride (SF6) concentrations. The results evidence a clear distinction between surface water (<≈10–15 m) and deeper water. The latter has much higher conductivity and cationic loads and is marked by lower Sr isotopic ratios and higher U activity ratios. Such a water typology suggests that the spring and piezometer water flows within the same shallow subsurface aquifer while deep borehole water belongs to a different circulation system, which flow path is controlled by the bedrock fracture network. The CFC data show that these two circulation systems are marked by contrasted residence times with a short residence time for the surface water and a longer residence time (water ages > 50 years) for the deep water. These results confirm different circulation histories for surface and deep water in the Strengbach catchment. They also suggest that the higher degree of chemical saturation of the deep water compared to the surface water is caused more by longer water residence times in the deep circulation systems than by differences in the primary minerals involved in each of the water-rock interaction systems. Our results also point that in the Strengbach granitic catchment, the SF6 concentrations cannot be used for water dating due to their lithogenic production in granitic bedrocks.

Geoelectrical Investigation of a Proposed Landfill Site in Ogbomoso, Southwestern Nigeria

An integrated geophysical investigation was carried out at a proposed landfill site along the new Ibadan-Ilorin Expressway in Ogbomoso, Southwestern Nigeria in order to characterize the subsurface and determine the suitability or otherwise of the site for landfill. The study location is underlain by rocks of the basement complex. 29 Schlumberger Vertical Electrical Soundings were conducted along seven traverses each about 100m long, about 20m apart, with station spacing of about 30m and half current electrode spacing varying between 1 and 100m. The VES data were interpreted quantitatively by partial curve matching and 1-D computer iteration and the layer parameters were used to generate the geoelectric sections. Wenner electrical resistivity profiling was also conducted along the traverses with station spacing of 25m and electrode spacing varied between 5 and 15m, and the data were interpreted using 2D inversion procedures. The VES results revealed three layers defined by the topsoil, clay layer and fractured/fresh bedrock. Depth to the bedrock is generally shallow with a mean of 9.0 3.7m. The mean thickness of the saprolite (mainly comprising clay) is 8.6 4.0m. 2D inversion of the resistivity profiling data revealed intense fracturing of the bedrock. The water table is shallow, ranging from 2 to 6m and may be rather close to the landfill base or fall within the landfill wastes. The study area is underlain by clayey overburden capable of protecting groundwater in the underlying aquifer unit from contamination by the leachate dischargeable from the landfill wastes. However, the presence of shallow water table and bedrock fractures may promote contamination of groundwater in the underlying aquifer unit and pose substantial health risk to the host community. It can therefore be concluded that the study area is not a suitable location for siting landfill. Keywords: Contamination, Leachate, Clay, Water table, Bedrock fractures DOI: 10.7176/JEES/11-16-04 Publication date: June 30 th 2021

Groundwater bacterial communities evolve over time in response to recharge

Time series analyses are a crucial tool for uncovering the patterns and processes shaping microbial communities and their functions, especially in aquatic ecosystems. Subsurface aquatic environments are perceived to be more stable than surface oceans and lakes, due to the lack of sunlight, the absence of photosysnthetically-driven primary production, low temperature variations, and oligotrophic conditions. However, periodic groundwater recharge should affect the structure and succession of groundwater microbiomes. To disentangle the long-term temporal changes in bacterial communities of shallow fractured bedrock groundwater, and identify the drivers of the observed patterns, we analysed bacterial 16S rRNA gene sequencing data for samples collected monthly from three groundwater wells over a six-year period (n = 230) along a hillslope recharge area. We showed that the bacterial communities in the groundwater of limestone-mudstone alternations were not stable over time and exhibited non-linear dissimilarity patterns which corresponded to periods of groundwater recharge. Further, we observed an increase in dissimilarity over time (generalized additive model P < 0.001) indicating that the successive recharge events result in communities that are increasingly more dissimilar to the initial reference time point. The sampling period was able to explain up to 29.5% of the variability in bacterial community composition and the impact of recharge events on the groundwater microbiome was linked to the strength of the recharge and local environmental selection. Many groundwater bacteria originated from the recharge-related sources (mean = 66.5%, SD = 15.1%) and specific bacterial taxa were identified as being either enriched or repressed during recharge events. Overall, similar to surface aquatic environments, the microbiomes in shallow fractured-rock groundwater vary through time, though we revealed groundwater recharges as unique driving factors for these patterns. The high temporal resolution employed here highlights the dynamics of bacterial communities in groundwater, which is an essential resource for the provision of clean drinking water; understanding the biological complexities of these systems is therefore crucial.

Research on the Settlement Regulation and Stability of Large Building Foundation over Gobs: A Case Study in the Xiangcheng Coal Mine, China

Due to the acceleration of urbanization in China, the land use and new building construction in mining subsidence areas are increasing. However, a large high-rise building (100 m × 80 m × 100 m) was planned to be built above gob areas in the Xiangcheng Mine subsidence area. According to the analysis of the borehole exploration data, borehole TV observations, and theoretical calculations, there were broken bedrock and developmental fractures in the fracture and caved zones above the gob areas. In addition, the depth affected by the building load is less than 9 m from the top of the fracture zone; thus, the stabilities of the gobs are low. Then, through the theoretical analysis of the space stress, compression settlement of foundation soil, and residual settlement of the gob and coal pillar under loads, combining the similar simulation study on settlement characteristics of foundation soil and gobs, the total settlement of the foundation can be regarded as the vector superposition of the compression settlement and residual settlement; then, the maximum settlement is calculated as 555 mm and the maximum tilt is 14.9 mm/m under the building load, which are greater than the permissible values of 200 mm and 2 mm/m, indicating that the foundation’s stability is inadequate. Therefore, the treatment measures of the deep pile group foundations in the soil layer and grouting reinforcement in the gobs were simultaneously put forward to effectively improve the stability of the building foundations.

Mapping Contaminant Plume at a Landfill in a Crystalline Basement Terrain in Ouagadougou, Burkina Faso, Using Self-Potential Geophysical Technique

The delineation of pollution plumes generated by household waste landfills is not easy, particularly in the case of discontinuous or intricately extending water tables, such as those developed in a fractured crystalline bedrock context. In Ouagadougou (Burkina Faso), there are many uncontrolled landfills throughout the urban area. The water table, generally located between 3 and 10 m deep, is likely to be contaminated by the leachate from these landfills. More than 1000 measurements of spontaneous potential (self-potential), referenced by GPS, have been carried out on a landfill and its immediate surroundings to the south of the urban area. The geostatistical processing by analysis of variograms and correlograms highlights an adapted prospecting technique and reliable cartography. The response seems to be mainly due to the electrochemical component with hot spots within the landfill and a plume heading towards the North-East. The distribution of the spontaneous potential seems to be controlled, not by the topography of the site, but by the fracturing of the mother rock of dominant direction 15° N, and by the mother rock/saprolite contact. Thus, the plume does not flow to the market gardening just below the landfill but rather to a residential area where monitoring of the quality of the borehole water is required.

Characterization of subsurface fluxes at the plot scale during flash floods in the Valescure catchment, France

AbstractThis study focuses on a 10‐m2 plot within a granitic hillslope in Cevennes mountainous area in France, in order to study infiltration and subsurface hydrological processes during heavy rainfalls and flash floods. The monitoring device included water content at several depths (0–70 cm for the shallow soil water; 0–10 m for the deep water) during both intense artificial and natural rainfall events, chemical and physical tracers, time‐lapse electrical resistivity tomography. During the most intense events, the infiltrated water was estimated to be some hundreds of millimetres, which largely exceeds the topsoil capacity (≤40 cm deep in most of the cases). The weathered/fractured rock area below the soil clearly has an active role in the water storage and sub‐surface flow dynamics. Vertical flow was dominant in the first 0–10 m, and lateral flow was effective at 8–10 m depth, at the top of the saturated area. The speed of the vertical flow was estimated between 1 and 10 m/hr, whereas it was estimated between 0.1 and 1 m/hr for the lateral flow. The interpretation of the experiments might lead to a local pattern of the 2D‐hydrological processes and profile properties, which could be generic for most of the mountainous catchments under Mediterranean climate. It suggests that fast triggering of floods at the catchment scale cannot be explained by a mass transfer within the hillslope, but should be due to a pressure wave propagation through the bedrock fractures, which allows exfiltration of the water downstream the hillslope.

Extreme sulfur isotope fractionation of hydrothermal auriferous pyrites from the SW fringe of the Taupo Volcanic Zone, New Zealand: Implications for epithermal gold exploration

The Taupo Volcanic Zone and its environs, New Zealand, represent an ideal place for epithermal gold exploration. We have documented structurally-bound gold in low temperature hydrothermal pyrite veins from the Erua Basin area on the outside fringe of the Taupo Volcanic Zone. Extreme 34S/32S fractionation of 122.2‰ from auriferous pyrite is documented in hydrothermal veins hosted in alluvial deposits from the Erua Basin. We document three generations of pyrite based on petrographic and Nano-SIMS microstructural analysis. First generation pyrite (Py assemblage -1) has extremely high δ34S values of +70.3 to +100.0‰ probably resulting from sequential-cyclic bacterial reduction and Rayleigh fractionation of limited sulfate. Second generation pyrite (Py assemblage -2) has relatively lower δ34S values of -39.8 to -3.3‰, interpreted to reflect the influx of magmatic-related oxidized fluids mixed with remobilized, reduced sedimentary sulfides that then formed the host calcite-pyrite veins. Third generation pyrite (Py assemblage -3) incorporated much heavier 34S resulting in positive δ34S values ranging from +6.8 to +72.2‰, indicating greater reduction of sulfate to H2S in an anoxic (reducing) environment. One single composite ~ 10 um pyrite grain encompassing all three types of subgrains has a large δ34S range from -22.2‰ along its rim and +100.0‰ from its core (122.2‰ for a single grain, and 139.8‰ for the system). This is the largest known variation in δ34S ever reported in a single pyrite grain from any system on Earth, but is comparable to δ34S ranges in pyrite systems of up to 186‰ reported from microbial sulfate reduction in deep bedrock fracture systems in Sweden. The growth stages from Py-2 to Py-3 indicate a transformation of the environment from an open to a closed (highly reducing) system, and it is these late-stage vein pyrites that are slightly auriferous. This process is interpreted to be induced by the precipitation of sulfides perhaps near the margin of the hydrothermal system or closed growth of pyrite in hydrothermal veins, which both restrict circulation and convection of sulfate within the veins. Our results suggest that the Erua Basin may sit on the fringe of a low-temperature epithermal gold system; other “low-temperature” areas around the TVZ may also preserve similar, potentially enriched zones on the edge of an epithermal Au belt.