Using environmental tracers to assess the extent of river–groundwater interaction in a quarried area of the English Chalk

The Swanscombe area of Kent, SE England represents a typical example of a heavily quarried Chalk area currently undergoing re-development. Because the Chalk is also an important aquifer, a good understanding of groundwater movement is required if environmental impacts are to be minimised and the water resource maximised. In particular, the nature of the relationship between the River Darent and groundwater in the Swanscombe Chalk Block requires better characterisation. Here, ‘environmental tracers’ in the form of ambient concentrations of stable isotopes, chlorofluorocarbons (CFCs), sulphur hexafluoride (SF₆) and tritium (³H) are used to investigate this and other aspects of groundwater movement in the vicinity of the quarries. Stable isotopic contrasts indicate little evidence for widespread river infiltration to the regional Chalk aquifer, although stable isotope and ³H data suggest that 20–35% of the abstraction by river-valley public water supply boreholes may be derived from the river. The CFCs, while present at above-modern concentrations in almost all groundwaters, can be used as tracers, indicating basically S–N flowpaths in the area south of the quarries, though sub-karstic conduits associated with areas of Palaeogene cover add a level of uncertainty at the local scale. Simple piston flow residence times based on SF₆ range from 1 to 17 a, but the data are probably better interpreted in terms of mixing between varying amounts of modern recharge derived from the south and deeper stored groundwater. The information gained from environmental tracers can therefore contribute to effective resource management.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

Towards hydrological and geochemical understanding of an ephemeral palustrine perched water table “wetland” (Lanseria Gneiss,Midrand, South Africa)

Wetland delineation is commonly based on terrain unit, soil form, soil wetness and vegetation indicators which along with the shallow groundwater or proximate surface water are often absent in many ephemeral inland South African areas due to, for instance, prolonged dry periods (seasons to years) which mask these indicators, as well as disruption of surface materials due to construction, agricultural activities and field fires. Furthermore, many “wetlands” in South Africa comply with the four indicators, but the notable absent requirement for wetlands is missing, namely the shallow groundwater table, as many of these systems form in hillslope seeps, catenas or from perched water tables. A 200-m-long excavation through one such a system is being studied in Midrand (Gauteng, South Africa). The site is underlain by Lanseria Gneiss and is waterlogged after long and intense periods of intense rainfall. Frequent downslope soil profiling, horizon-specific laboratory analyses for grading and Atterberg Limits, X-Ray Diffraction and X-Ray Fluorescence Spectroscopy data are interpreted together with soil percolation tests to generate a conceptual model of the system. The results are discussed in terms of the need to consider these temporary systems that do not have a shallow groundwater table and that are not in direct contact with surface drainage features as a possible special type of wetlands, notably in arid regions where they play a very important role in biodiversity and should, therefore, be protected.     

Strain and deformation mechanisms in the Variscan nappes of Vendée,South Brittany,France

The Variscan nappes of Vendée (South Brittany, France) were displaced westward through a mechanism of regional subhorizontal ductile shear involving the entire body of the nappes. Ductile shear also affected the upper part of the Silurian footwall, including the rhyolite of Chapelle-Hermier, which was mylonitized and sheared with upward-increasing intensity, as well as the calcalkaline Vendée porphyroids of volcanic and volcanoclastic origin which form a mylonite sheet outlining the major basal thrust zone of the nappes. Nappe displacement was accompanied by a severe extensional strain, estimated at up to 300–400% in both the porphyroids and the rhyolite, and indicating that the deformational history may have involved a gravitational spreading component. Within the mylonitic rhyolite, accommodation of the bulk strain was mainly by solution-assisted mass-transfer, as shown by (1) growth of quartz fibres in low-pressure zones, (2) development of polycrystalline ribbons through continuous growth of quartz grains in pressure shadows and (3) crystallization of muscovite at the expense of the felsitic groundmass. A similar deformation process was probably active in the mylonitic porphyroids at the sole of the nappe, suggesting that pressure-solution flow was an important mechanism in nappe displacement.     

Integration of geochemical and isotopic tracers for elucidating water sources and salinization of shallow aquifers in the sub-Saharan Drâa Basin,Morocco

In the arid sub-Saharan of southern Morocco, groundwater salinization poses a direct threat for agricultural production in six oases’ basins that are irrigated by water imported from the High Atlas Mountains. Here the geospatial distribution of salinity is evaluated in shallow groundwater, springs and surface waters in the Drâa Basin, integrating major and trace element geochemistry and isotopic tracers (O, H, Sr and B). The data show that water discharge from the High Atlas Mountains to the Upper section of the Drâa Basin is characterized by both low and high salinity, a distinctive low δ¹⁸O and δ²H composition (as low as −9‰ and −66‰, respectively), typical for meteoric water from high elevation, a ⁸⁷Sr/⁸⁶Sr range of 0.7078–0.7094, and δ¹¹B of 12–17‰. The Ca–Mg–HCO₃, Na–Cl–SO₄, and Ca–SO₄ compositions as well as the Br/Cl, ⁸⁷Sr/⁸⁶Sr, and δ¹¹B values of the saline water suggest dissolution of Lower Jurassic carbonates and evaporite rocks in the High Atlas Mountain catchment. Storage and evaporation of the imported water in a man-made open reservoir causes an enrichment of the stable isotope ratios with a δ¹⁸O/δ²H slope of <8 but no change in the Sr and B isotope fingerprints. Downstream from the reservoir, large salinity variations were documented in the shallow groundwater from the six Drâa oases, with systematically higher salinity in the three southern oases, up to 12,000 mg/L. The increase of the salinity is systematically associated with a decrease of the Br/Cl ratio, indicating that the main mechanism of groundwater salinization in the shallow aquifers in the Drâa oases is via salt dissolution (gypsum, halite) in the unsaturated zone. Investigation of shallow groundwater that flows to the northern Drâa oases revealed lower salinity (TDS of 500–4225) water that is characterized by depleted ¹⁸O and ²H (as low as −9‰ and −66‰, respectively) and higher ⁸⁷Sr/⁸⁶Sr ratios (∼0.7107–0.7115) relative to irrigation water and groundwater flow from the Upper Drâa Basin. This newly-discovered low-saline groundwater with a different isotopic imprint flows from the northeastern Anti-Atlas Jabel Saghro Mountains to the northern oases of the Lower Drâa Basin. This adjacent subsurface flow results in a wide range of Sr isotope ratios in the shallow oases groundwater (0.7084–0.7131) and appears to mitigate salinization in the three northern Drâa oases. In contrast, in the southern oases, the higher salinity suggests that this mitigation is not as affective and increasing salinization through cycles of water irrigation and salt dissolution appears inevitable.     

Effects of the October 25, 2010, Mentawai tsunami in La Réunion Island (France): observations and crisis management

On October 25, 2010, a Mw = 7.7 magnitude earthquake was recorded in the Kepulauan Mentawai archipelago, in Indonesia. Following the earthquake, the Pacific Tsunami Warning Center issued a tsunami bulletin to the countries of the Indian Ocean informing them that a local tsunami watch was in effect for Indonesia. La Réunion Island, a French territory located 5,000 km southwest of the seismic source, did not trigger its tsunami early warning protocol, but did suffer the impact of this tsunami 7h20 post-earthquake. This paper presents the post-tsunami survey conducted in La Réunion Island a few days after the event, providing run-up measurements that reach 1.72 m locally as well as information about the material losses and the crisis management. It provides a time sequence of the observations along the coast and the reactions of the population and authorities in dealing with this crisis. This paper also shows the local authorities’ lack of action in managing this event—thus increasing human and material vulnerabilities—the consequences for victims and the impact on the credibility of the local early warning system. These results provide geophysicists with information to better describe the source of this tsunami through far-field inversion and provide local decision makers with a comprehensive assessment of the weaknesses of the local early warning system and protocols.     

High-resolution cyclostratigraphic analysis from magnetic susceptibility in a Lower Kimmeridgian (Upper Jurassic) marl–limestone succession (La Méouge,Vocontian Basin,France)

High-resolution magnetic susceptibility (MS) analysis was carried out on a Lower Kimmeridgian alternating marl–limestone succession of pelagic origin that crops out at La Méouge (Vocontian Basin, southeastern France). The aim of the study was to characterize the strong, dm-scale sedimentary cyclicity of the succession at a very high resolution, and to analyze the cycles for evidence of astronomical forcing. From marl to limestone, MS varies progressively and closely tracks the highest frequency cyclicity corresponding to the basic marl–limestone couplets. Long-term wavelength cycling modulates the high-frequency cyclicity (couplets), and appears to be controlled by clay content. Spectral analysis of the MS record reveals the presence of the complete suite of orbital frequencies in the precession, obliquity, and eccentricity (95–128 ka and 405 ka) bands with very high amplitude of the precession index cycles originating from dm-scale couplets. 405 ka-eccentricity cycles are very pronounced in the MS maxima of the marl members of the couplets, suggesting eccentricity-driven detrital input to the basin. 405 ka-orbital tuning of the MS maxima further sharpens all of the orbital frequencies present in the succession. These results are similar to those of previous studies at La Méouge that used carbonate content observed in field. Our results are also in accordance with cyclostratigraphic studies in Spain and Canada that report dominant precession index forcing. By contrast, in the Kimmeridge Clay (Dorset, UK), obliquity forcing dominates cyclic sedimentation, with weaker influence from the precession index. Ammonite zone duration estimates are made by counting the interpreted precession cycles, and provide an ultra-high resolution assessment of geologic time. In sum, this study demonstrates the power of the MS as a proxy in characterizing the high-resolution cyclostratigraphy of Mesozoic sections, particularly in alternating marl–limestone successions, and for high-resolution correlation and astronomical calibration of the geologic time scale.     

Groundwater flow in an ‘underfit’ carbonate aquifer in a semiarid climate: application of environmental tracers to the Salt Basin,New Mexico (USA)

The Salt Basin is a semiarid hydrologically closed drainage basin in southern New Mexico, USA. The aquifers in the basin consist largely of Permian limestone and dolomite. Groundwater flows from the high elevations (～2,500 m) of the Sacramento Mountains south into the Salt Lakes, which are saline playas. The aquifer is ‘underfit’ in the sense that depths to groundwater are great (～300 m), implying that the aquifer could transmit much more water than it does. In this study, it is speculated that this characteristic is a result of a geologically recent reduction in recharge due to warming and drying at the end of the last glacial period. Water use is currently limited, but the basin has been proposed for large-scale groundwater extraction and export projects. Wells in the basin are of limited utility for hydraulic testing; therefore, the study focused on environmental tracers (major-ion geochemistry, stable isotopes of O, H, and C, and ¹⁴C dating) for basin analysis. The groundwater evolves from a Ca–HCO₃ type water into a Ca–Mg (Na) – HCO₃–Mg (Cl) water as it flows toward the center of the basin due to dedolomitization driven by gypsum dissolution. Carbon-14 ages corrected for dedolomitization ranged from less than 1,000 years in the recharge area to 19,000 years near the basin center. Stable isotopes are consistent with the presence of glacial-period recharge that is much less evaporated than modern. This supports the hypothesis that the underfit nature of the aquifer is a result of a geologically recent reduction in recharge.     

Trace elements and environmental isotopes as tracers of surface water–groundwater interaction: a case study at Xin’an karst water system,Shanxi Province,Northern China

In this case study, silica concentration, oxygen and strontium isotopes of water samples were used to study surface water–groundwater interaction at the Xin’an karst water system. The silica concentration in rain water is commonly less than 1 mg/l. In the areas around the south tributary of the Zhuozhang River, silica concentrations in the groundwater in Quaternary aquifers range between 4.04 and 7.66 mg/l while that of the surface water varies from 1.49 to 6.9 mg/l. Silica concentrations of most surface water samples increase with TDS, indicating the effect of groundwater recharge on river water chemistry. On the contour map of silica concentration of groundwater in Quaternary aquifers, samples located close to surface water often have lower silica concentrations as a result of surface water recharge. Both overland flow and surface water have impact on karst water according to our hydrogeochemical study of stable oxygen isotope, Sr isotope and strontium contents. Calculation results of three end member mixing model show that the contribution of karst water, surface water and overland flow water is 45, 28 and 27%, respectively.     

Geochronological and geochemical characterization of magmatic-hydrothermal events within the Southern Variscan external domain (Cévennes area,France)

Geochronological, mineralogical, and geochemical analyses have been focussed on the Mont-Lozère-Borne plutonic complex and surrounding rocks (Cévennes, French Massif Central) in which B–W-Sn and As-Au-Sb mineralization is encountered. Two main results are highlighted: (1) the existence of a 301–306 Ma magmato-hydrothermal event unrelated to the emplacement of the Pont-de-Montvert-Borne plutonic body at 316 Ma; (2) the magmatic and hydrothermal features are strongly associated, both in time and in space, thus demonstrating an intimate connection between mineralizing processes and magmatism in this part of the French Massif Central. We also show that mineralization and associated hydrothermal occurrences do not correspond to a simple and single geochemical signature and that a contamination model must be invoked in order to account for the complexity of isotopic results. This study demonstrates that the application of the O and H isotopic signatures as tracers of the source and nature of fluids in an orogenic context requires some specific care. Finally, a model of the tectonic-magmatic-hydrothermal evolution of the study area is suggested in which we discuss two alternative scenarios. The first one implies the existence of two different hydrothermal/mineralizing events (Bo-W-Sn and As-Au-Sb ones). The second one suggests the same source for all hydrothermal and mineralized structures.     

Permeable Reactive Barrier for the treatment of contaminated surface water in Katedan industrial development area, Hyderabad, India

Permeable Reactive Barrier （PRB） is an emplacement of inert material （s） in the subsurface, designed to intercept a contaminated plume, provides a preferential flow path through the reactive media, and transforms the contaminant into environmentally acceptable forms to attain concentration remediation goals at the discharge of the barrier. The phenomena, which help in remediation within PRB, are adsorption/sorption, precipitation, oxidation/reduction and biodegradation. Various materials like zero-valent iron, zero-valent bi-metals, natural zeolites, organic carbon, fly ash, zeolites, limestone, activated alumina, apatites, etc. have been tried by many researchers to remove organic and inorganic contaminants. In USA, Canada, and many European countries commercial full-scale and pilot scale PRBs are successfully working. The design and installation of full scale PRBs needs laboratory treatability and dynamic flow column experiments？ The concept of PRB is being applied to treat contaminated surface water in the Katedan industrial area, Hyderabad, India. National Geophysical Research Institute （NGRI）, Hyderabad, India, conducted systematic studies in collaboration with Norwegian Geotechnical Institute （NGI）, Norway, to develop PRB technique to decontaminate the surface water pollution due to industrial effluent. A site assessment study in the Katedan Industrial Area, were carried out and water, soil and sediment from the lakes of the area were found to be polluted with high concentrations of heavy metals like As, Pb, Cr, Cd, Ni, etc. Adsorption studies at NGRI with synthetic samples and in-situ industrial effluent using different reactive media for removing contaminants like arsenic, chromium, cadmium, copper, nickel, lead and zinc have been carried out and yielded satisfactory results. The performance of zero-valent iron and limestone is encouraging in removing As,     

Influence of car body repair activities on particulate matter pollution in a suburban area (Bouzaréah City,Algiers)

This work was carried out in order to study the impact of car body repair activities in Bouzaréah City (suburban area located at the plateau of Algiers) on the air quality. Atmospheric particulate matter was targeted using high-volume sampler. Precisely, the particulate organic matter such as organophosphate esters, especially triphenyl phosphate and tributyl phosphate, was evaluated by the implementation of accelerated solvent extraction. The results of the major extractable organic components revealed that tributyl phosphate and triphenyl phosphate were found at 5.46 and 1.78 ng m^?3, respectively. On the other hand, the classical pollutants such as n-alkanes, n-alkanoic monocarboxylic acids and polycyclic aromatic hydrocarbons were also measured, respectively, at average total concentrations of 35.8, 107.4 and 9.8 ng m^?3. This work confirms that sanding and painting car bodies may contribute to air pollution, especially after detection of diethylnitrosamine which remains the most surprising compound in the particulate matter due to its volatility and can be considered as emerging contaminant.     

Hydrogeochemical characterization and evolution of a regional karst aquifer in the Cuatrociénegas area,Mexico

The Cuatrociénegas area is useful for the investigation of the effect of groundwater extraction in the Chihuahuan freshwater xeric ecoregion. It has been investigated at this time using a selection of geochemical indicators (major, minor and trace elements) and δ³⁴S data, to characterize the origin of groundwater, the main geochemical processes and the mineral/groundwater interactions controlling the baseline geochemistry. The area is composed of limestones of Mesozoic age, with a composite thickness of about 500 m, overlaid by basin fill (poorly consolidated young sediments). Substantial water extraction and modification of natural discharges from the area along the last century have produced a detrimental impact on ecosystem structure and function. Water–rock interactions, mixing and evaporative processes dominate the baseline groundwater quality. Natural recharge is HCO₃–Ca type in equilibrium with calcite, low salinity (TDS?<?500 mg/L), Cl^? lower than 11 mg/L and average Li⁺ concentration of 0.005 mg/L. Along the groundwater flow systems, δ³⁴S evidence and mass transfer calculations indicate that Cretaceous gypsum dissolution and dedolomitization reactions adjust water composition to the SO₄–Ca type. The increase of water–rock interaction is reflected by Cl^? values increase (average 68 mg/L), TDS up to about 1500 mg/L and an average Li⁺ concentration of 0.063 mg/L. Calculations with chemical geothermometers indicate that temperature at depth could be at maximum of 15–20 °C higher than field-measured temperature for pozas. After groundwater is discharged to the surface, chemical evolution continues; water evaporation, CO₂ degassing and precipitation of minerals such as gypsum, calcite and kaolinite represent the final processes and reactions controlling water chemical composition.     

Melt percolation and reaction atop a plume: evidence from the poikiloblastic peridotite xenoliths from Borée (Massif Central, France)

Poikiloblastic harzburgite xenoliths (P-type) from Borée, France are characterised by large (>1 cm), essentially unstrained olivines and high equilibrium temperatures (>1200 °C). Mineralogical data, trace element abundances and Sr-Nd-O isotopes of the constituent minerals are consistent with formation as a result of melt percolation-reactions in a lherzolite precursor during lithospheric erosion by an upwelling plume. This petrogenetic model contrasts with previous models involving isochemical recrystallisation from a granular lherzolite precursor (G-type) or derivation as metacumulates from tholeiitic magmas. Numerical simulation of percolation reactions at the lithosphere-plume boundary using the plate model of Vernières et al. (1997) indicates that the different textured xenoliths may represent mantle from different levels in a percolation-reaction column. If correct then the P-type harzburgites resulted from pyroxene-dissolving and olivine-producing reactions at increasing melt fraction (>3%) at the lower part of column (base of the lithosphere), whereas the G-type lherzolites were located within the low-porosity domain (<0.1%) above a permeability barrier, and are formed through a melt-rock reaction at decreasing melt mass. Given the very low melt fraction, the REE fractionation in this zone is controlled by chromatographic effects coupled with source effects of reaction. The variations in porosity, melt/rock ratio and melt-rock reaction mechanism are believed to be responsible for the diversity of REE patterns and striking correlation between REE abundance and texture in Borée xenoliths. Received: 15 June 1997 / Accepted: 7 January 1998     

Application of the content–area (C–A) fractal model and prediction–area (P–A) plot for mineral prospectivity modeling in the Luchun area of Yunnan Province,China

This method of assigning weights based on expert opinion introduces bias when we are evaluating the relative importance of evidence values. In this paper, we used a prediction–area (P–A) plot method and content–area (C–A) fractal model to estimate the weight of each evidence map. In this paper, we used the content region (C–A) fractal model to divide the evidence maps to the threshold of the corresponding dimensions. The P–A plot approach is an objective data-driven approach for evaluating map weights. Using geochemical layer and remote sensing, hydroxyl layers as weight evidence maps are the highlights of this study. We use the P–A method from which we can evaluate the predictive ability of each evidence map with respect to the known ore occurrences. We used the P–A plot for weighting each evidence map and choosing the appropriate threshold for predictor maps in the Luchun area of Yunnan Province, China. The method adopted in this paper can improve the prediction efficiency of ore prospecting.     

Multi-tracer tests to evaluate the hydraulic setting of a complex aquifer system (Brévilles spring catchment, France)

For good management of groundwater resources, and to comply with European and national regulations, a detailed understanding of an aquifer’s hydraulic setting is required. In order to better characterize a sandy aquifer that is affected by diffuse pollution (Brévilles spring catchment, Val d’Oise, France), and to quantify the transfer time in the saturated zone, a multi-tracer test involving a new technique, the ‘finite volume point dilution method’, has been performed in natural flow conditions. In November 2005, injections of four different tracers took place in four piezometers involving different locations and depths in the aquifer. Recovery of the tracers was observed at two different places near the aquifer outlet. A particularly long and unusual monitoring exercise (27 months) demonstrated the existence of several different velocities within the sandy layer, which seems to be linked to the decrease of hydraulic conductivity with depth. The new insight and parameter quantification brought by interpretation of these tests contribute to a better characterization of the saturated zone. The particularly long-term monitoring exercise also gives new information to understand and forecast the trend and persistence of groundwater contamination by pesticides in the catchment.     

Petrogenesis of Fe–Ti oxides in amphibole-rich veins from the Lherz orogenic peridotite (Northeastern Pyrénées, France)

Accessory, homogeneous ilmenite and rutile are important oxide phases in amphibole-rich high-pressure cumulate veins which crosscut the Lherz orogenic lherzolite massif. Those veins crystallized from alkaline melts at P = 1.2–1.5 GPa within the uppermost lithospheric mantle. Transitional basalts contaminated by peridotitic wall-rocks and then uncontaminated alkali basalts (basanites) reused the same vein conduits. Petrographic observations give evidence that Fe–Ti oxide saturation depends on the silica contents of each parental melt. The water-poor silica-rich transitional melts that generated websterites and plagioclase-rich clinopyroxenites reached early Ti-oxide saturation (1,200°C; 1.5 GPa). Rutile is as abundant as ilmenite. It is enriched with Nb–Zr–Hf by a factor of 10–100 relative to either amphibole or ilmenite. The amphibole pyroxenites and hornblendites crystallized from basanites reached late Fe–Ti oxide saturation after precipitation of amphibole, with ilmenite crystallizing along with phlogopite in the latter. The Lherz ilmenites are devoid of exsolution and contain very little trivalent iron. This compositional feature indicates more reducing crystallization conditions than usually inferred for alkali lavas and their megacrysts (FMQ ± 1). The veins incompletely equilibrated for redox conditions with their wall-rock peridotites which record more oxidizing conditions (FMQ ± 1). The veins also exchanged magnesium and chromium, as suggested by Cr-bearing, Mg-rich ilmenite (up to 44 mol% MgTiO₃₎ in veins less than 3–4 cm thick. Mg-rich ilmenite megacrysts occurring in alkali basalts could be actually xenocrysts from veins similar in thickness to those occurring at the Lherz massif, although crystallized from more oxidized magmas.