The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Trace elements in groundwater as indicators of anthropogenic impact

The distribution of several minor and trace elements mainly in fresh (dominating TDS 160–400 mg/l) groundwater of Latvia have been investigated by the Inductively Coupled Plasma-Mass Spectrometry (ICP-MS) technique. An evaluation of results of about 700 analyses leads to the conclusion that concentrations of these elements is influenced by: pH–Eh conditions, groundwater residence time and diffuse contamination, whereas the role of water-bearing sediments is of secondary importance. Most trace elements are characterised by low mobility under alkaline and reducing conditions; concentrations in confined aquifers are much smaller than the Maximum Permissible Values for drinking water. The strongest anomalies of REE, Al and P were found in shallow groundwater around the former agrochemical storehouses.     

Numerical modeling assessment of three-gate structures for capturing contaminated groundwater

This modeling study evaluated the capability of alternative funnel-and-gate structures with three gates for capturing contaminated groundwater in a hypothetical unconfined aquifer. Simulated interceptor structures were linear and 45 m wide, consisting of three gates and two funnels (walls). One gate occupied the center and two gates occupied the ends of the interceptor structures. The structures, positioned perpendicular to regional groundwater flow, traversed the entire thickness of the aquifer. A total of four structures were evaluated (numbers designate widths of end, center, and end gates, respectively, in meters): 3-3-3, 2-5-2, 1-7-1, and 4-1-4. Particle tracking and zonal water budgets identified shapes of capture zones and discharge patterns for each interceptor structure. A mass transport model, accounting for advection and hydrodynamic dispersion, tested the capability of each structure for capturing a contaminant plume. Results suggest that: time-dependent capture zones underestimate the amount of time to capture a contaminant plume, wide center gates facilitate plume capture, and wide end gates facilitate lateral containment of contaminants. Of the structures simulated, the 2-5-2 configuration was relatively efficient at processing and containing the simulated contaminant plume.     

Groundwater recharge of carbonate aquifers of the Silesian-Cracow Triassic (southern Poland) under human impact

A Triassic carbonate unit has been intensively drained by zinc and lead ore mines and numerous borehole fields since the nineteenth century. Its groundwater recharge has increased due to: pumping of water from boreholes, mining activity, and urbanization. An approach to determine the amounts of the recharge at a variety of spatial scales is presented in the paper. Different methods were used to identify and quantify recharge components on a regional and local scale: mathematical modelling was performed for four aquifers included in an aquifer system, an analytical estimation based on the assumption that an average recharge is equal to the average discharge of the hydrogeological system—for six man-made drainage centres, and the method of water level fluctuation (WLF) was applied in one observation borehole. Results of modelling have been supplemented by observation of environmental tracers (δ¹⁸O, δ²H, ³H), noble gases temperatures, and ⁴He_exc in groundwater. The regional aquifer’s current recharge according to estimations performed by means of modelling varies from 39 to 101 mm/year on average. Depending on the aquifer site the average precipitation ranges from 779 to 864 mm/year. In the confined part of the aquifer average recharge ranges from 26 to 61 mm/year. Within outcrops average recharge varies from 96 to 370 mm/year. Current recharge estimated by the analytical method for man-made drainage centres varies from 158 up to 440 mm/year. High values are caused by different recharge sources like precipitation, induced leakage from shallow aquifers, and water losses from streams, water mains and sewer systems. Pumping of water, mining and municipal activities constitute additional factors accounting for the intensified recharge.     

Mapping susceptibility of rainfall-triggered shallow landslides using a probabilistic approach

To prepare a landslide susceptibility map is essential to identify hazardous regions, construct appropriate mitigation facilities, and plan emergency measures for a region prone to landslides triggered by rainfall. The conventional mapping methods require much information about past landslides records and contributing terrace and rainfall. They also rely heavily on the quantity and quality of accessible information and subjectively of the map builder. This paper contributes to a systematic and quantitative assessment of mapping landslide hazards over a region. Geographical Information System is implemented to retrieve relevant parameters from data layers, including the spatial distribution of transient fluid pressures, which is estimated using the TRIGRS program. The factor of safety of each pixel in the study region is calculated analytically. Monte Carlo simulation of random variables is conducted to process the estimation of fluid pressure and factor of safety for multiple times. The failure probability of each pixel is thus estimated. These procedures of mapping landslide potential are demonstrated in a case history. The analysis results reveal a positive correlation between landslide probability and accumulated rainfall. This approach gives simulation results compared to field records. The location and size of actual landslide are well predicted. An explanation for some of the inconsistencies is also provided to emphasize the importance of site information on the accuracy of mapping results.     

Experimental study on the deterioration and natural remediation of the Ariake Sea tidal mud caused by the sea laver treatment acid practice and the upward seepage of pore water liquid

This study presents a laboratory study of the following two aspects: (1) the influence of sea laver treatment acid on the geoenvironmental properties of Ariake Sea tidal mud, and (2) the natural remediation effect on the sea laver treatment acid contaminated Ariake Sea tidal mud caused by the upward seepage of pore water liquid in the mud. Firstly, the mechanisms of the transport of sea laver treatment acid in the Ariake Sea tidal mud and the generation mechanisms of the upward seepage flow in the Ariake Sea tidal mud are discussed. Secondly, a series of one-dimensional laboratory infiltration tests were carried out to investigate the deterioration of the Ariake Sea tidal mud caused by the sea laver acid treatment practice. Test results reveal that the acid treatment practice caused considerable change in the geochemical properties of the mud in terms of increase in sulfide content and decrease in pH value. After the treatment by the sea laver treatment acid, the sulfide content of the mud even exceeded the safe limit value for the benthos, which represents undesirable living condition for benthos. Thirdly, series of laboratory fresh seawater infiltration tests for the deteriorated Iida site mud were conducted to illustrate this natural remediation efficiency. It is found that with the infiltration of the fresh seawater, the sulfide content of the Iida site mud was considerably reduced and pH value increased to an acceptable range for benthos living in the tidal flat mud. With the increase in the infiltration time and the hydraulic gradient, the remediation efficiency could be increased.     

Effects of natural attenuation processes on groundwater contamination caused by abandoned waste sites in Berlin

The aim of this research project is to identify, characterize and quantify natural attenuation (NA) processes in groundwater affected by emissions of abandoned waste disposal sites in Berlin-Kladow/Gatow, Germany. It is part of the funding priority called KORA established by the Federal Ministry for Education and Research (BMBF) to explore the extent to which NA can be used for remedial purposes for varied forms of soil and groundwater contamination. Information on the emission behaviour of individual parameters is generated on the basis of hydrogeochemical comparison of 20 years old and new data. Using groundwater-modelling and CFC-analysis, information on the transport and retention of pollutants in groundwater is compiled. The microbial colonization of contaminated aquifers is characterized by molecular biological methods [polymerase chain reaction (PCR) and denaturing gradient gel electrophoresis (DGGE)] to differentiate between contaminated and not contaminated zones.     

The Ca-Eskola component in eclogitic clinopyroxene as a function of pressure, temperature and bulk composition: an experimental study to 15 GPa with possible implications for the formation of oriented SiO2-inclusions in omphacite

Experiments have been conducted in the P-T range 2.5–15 GPa and 850–1,500°C using bulk compositions in the systems SiO₂–TiO₂–Al₂O₃–Fe₂O₃–FeO–MnO–MgO–CaO–Na₂O–K₂O–P₂O₅ and SiO₂–TiO₂–Al₂O₃–MgO–CaO–Na₂O to investigate the Ca-Eskola (CaEs Ca_0.5□_0.5AlSi₂O₆) content of clinopyroxene in eclogitic assemblages containing garnet + clinopyroxene + SiO₂ ± TiO₂ ± kyanite as a function of P, T, and bulk composition. The results show that CaEs_ss in clinopyroxene increases with increasing T and is strongly bulk composition dependent whereby high CaEs-contents are favoured by bulk compositions with high normative anorthite and low diopside contents. In this study, a maximum of 18 mol% CaEs_ss was found at 6 GPa and 1,350°C in a kyanite-eclogite assemblage garnet + clinopyroxene + kyanite + rutile + coesite. By comparison, no significant increase in CaEs_ss with increasing P could be observed. If the formation of oriented SiO₂-rods frequently observed in eclogititc clinopyroxenes is due to the retrogressive breakdown of a CaEs-component then these textures are a cooling rather than a decompression phenomenon and are most likely to be found in kyanite-bearing eclogites cooled from temperatures ≥750°C. The presence of clinopyroxene with approx. 4 mol% CaEs_ss in an experiment conducted at 2.5 GPa/850°C confirms earlier suggestions based on field data that vacancy-rich clinopyroxenes are not necessarily restricted to ultrahigh pressure metamorphic conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Stable isotope study of the Archaean rocks of the Vredefort impact structure,central Kaapvaal Craton,South Africa

The Vredefort dome in the Kaapvaal Craton was formed as a result of the impact of a large meteorite at 2.02 Ga. The central core of Archaean granitic basement rocks is surrounded by a collar of uplifted and overturned strata of the Witwatersrand Supergroup, exposing a substantial depth section of the Archaean crust. Orthogneisses of the core show little variation in whole-rock δ ¹⁸O value, with the majority being between 8 and 10‰, with a mean of 9.2‰ (n = 35). Quartz and feldspar have per mil differences that are consistent with O-isotope equilibrium at high temperatures, suggesting minimal interaction with fluids during subsequent cooling. These data refute previous suggestions that the Outer Granite Gneiss (OGG) and Inlandsee Leucogranofels (ILG) of the core represent middle and lower crust, respectively. Granulite-facies greenstone remnants from the ILG have δ ¹⁸O values that are on average 1.5‰ higher than the ILG host rocks and are unlikely, therefore, to represent the residuum from the partial melting event that formed the host rock. Witwatersrand Supergroup sedimentary rocks of the collar, which were metamorphosed at greenschist-to amphibolite-facies conditions, generally have lower δ ¹⁸O values than the core rocks with a mean value for metapelites of 7.7‰ (n = 45). Overall, through an ∼20 km thick section of crust, there is a general increase in whole-rock δ ¹⁸O value with increasing depth. This is the reverse of what is normal in the crust, largely because the collar rocks have δ ¹⁸O values that are unusually low in comparison with metamorphosed sedimentary rocks worldwide. The collar rocks have δD values ranging from −35 to −115‰ (average −62‰, n = 29), which are consistent with interaction with water of meteoric origin, having a δD of about −25 to −45‰. We suggest that fluid movement through the collar rocks was enhanced by impact-induced secondary permeability in the dome structure. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Li-Sr-Nd isotope signatures of the plume and cratonic lithospheric mantle beneath the margin of the rifted Tanzanian craton (Labait)

Lithium concentrations and isotopic compositions of olivine and ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd of coexisting clinopyroxene from peridotite xenoliths from the Quaternary Labait volcano, Tanzania, document the influence of rift-related metasomatism on the ancient cratonic mantle. Olivines show negative correlations between Fo content and both δ⁷Li and Li concentrations. Olivines in iron-rich peridotites (Fo_85–87) have high Li concentrations (3.2–4.8 ppm) and heavy δ⁷Li (+5.2 to +6.6). In contrast, olivines in ancient, refractory peridotites have lower Li concentrations (∼2 ppm) and relatively light δ⁷Li (+2.6 to +3.5). This reflects mixing between ancient, refractory cratonic lithosphere and asthenosphere-derived rift magmas. A uniquely fertile, deformed, high-temperature garnet lherzolite, interpreted to be from the base of the lithosphere, has a ⁸⁷Sr/⁸⁶Sr of 0.7029 and ¹⁴³Nd/¹⁴⁴Nd of 0.51286, similar to HIMU oceanic basalts. It provides the best estimate of the Sr–Nd isotope composition of the upwelling mantle (i.e., plume, sensu lato) underlying this portion of the East African Rift, and is slightly less radiogenic compared to previous estimates of the plume that were based on rift basalts. Although elevated δ⁷Li are not exclusive to HIMU source regions, the data collectively indicate that the plume beneath Labait has HIMU characteristics in Sr, Nd and Li isotope composition. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.