Advanced argillic alteration at Cave di Caolino, Lipari, Aeolian Islands (Italy): Implications for the mitigation of volcanic risks and the exploitation of geothermal resources

Four sites in the western sector of Lipari Island with still active hydrothermal activity are here considered. The petrography (mesoscopic observations and XRPD) and geochemistry (major, minor and trace elements chemistry) of ten representative and extremely altered volcanic rocks were characterized. Two types of parageneses of altered rocks are discriminable, one rich in silicate phases (opal/cristobalite, montmorillonite, kaolinite, alunite and hematite) and one in sulphates (gypsum, plus minor amounts of anhydrite or bassanite). The altered silicate-rich rocks are rich in SiO₂, Al₂O₃, Fe₂O₃ and H₂O, and depleted in CaO, MgO, K₂O and Na₂O, while the sulphate-rich ones are extremely enriched in CaO and SO₄ in comparison with local unaltered volcanic rocks. The content of many incompatible elements is similar in altered silicate-rich rocks and lower in sulphate-rich ones with respect to the pristine volcanic rocks; conversely, almost all REEs are markedly enriched in silicate-rich rocks and heavy REEs are enriched in sulphate-rich altered rocks compared to unaltered volcanic rocks. Reaction path modelling of basaltic andesite dissolution in local steam condensate predicts the production of amorphous-silica, anhydrite, goethite, and kaolinite (or smectites and saponites) as stable secondary minerals and alunite, jarosite, and jurbanite as ephemeral minerals. Considering possible post-depositional reactions and admitting that the presence of two distinct parageneses is apparent, since gypsum is prone to form large crystals, it can be concluded that there is an excellent agreement between the alteration minerals occurring in nature and those predicted by geochemical modelling. Consequently, the modelled process is the main responsible for the production of the advanced argillic alteration assemblage of "Cave di Caolino" on Lipari Island. Since rock alteration is sustained by the H₂SO₄ solution produced by hydrothermal steam condensation, there is no need to invoke the involvement of SO₂-HCl-HF-bearing magmatic fluids, in line with the absence of fluoride minerals.

Hydrogeochemical and isotopic characterization of the Gioia Tauro coastal Plain (Calabria - southern Italy): A multidisciplinary approach for a focused management of vulnerable strategic systems

This work pursues the hydro-geochemical and isotopic characterization of the complex groundwater system of the Gioia Tauro Plain, one of the most important industrialized and agricultural coastal areas of southern Italy. The anthropic pressure exposes the water resources at risk of depletion and quality degradation making the plain groundwater a system of high scientific and social interest. The plain is characterized by a shallow aquifer, mostly recharged by local rains and a deep aquifer apparently less influenced by local precipitation. Both aquifers are mainly Ca-HCO₃ waters except for localized sectors where Na-HCO₃, Na-Cl and Ca-SO₄ waters are present. In deep aquifer, both prolonged interaction with sedimentary rocks, mainly deriving from the erosion of crystalline rocks, and direct cation exchange represent the primary factors controlling the formation of Na-HCO₃ waters. Mixing processes between these waters and either connate brine and/or deep thermal waters contribute to the formation of isolated high salinity Na-Cl-rich waters. In shallow aquifer, inputs of N-rich sewage and agriculture-related contaminants, and SOx emissions in proximity of the harbor are responsible of the increasing nitrate and sulphate concentrations, respectively. The Cl/Br and NO₃/Cl ratios highlight contamination mainly linked to agricultural activities and contribution of wastewater. Along the northern boundary, the warmest groundwater (Na-Cl[SO₄]) were found close to a bend of the main strike-slip fault system, locally favouring the rising of B- and Li-rich deep waters, testifying the influence of geological-structural features on deep water circulation. Despite the high-water demand, a direct marine intrusion is localized in a very restricted area, where we observed an incipient groundwater-seawater mixing (seawater contribution ≤7 %). The qualitative and quantitative conditions of the shallow aquifer still have acceptable levels because of the relatively high recharge inflow. A reliable hydrogeochemical conceptual model, able to explain the compositional variability of the studied waters, is proposed.

Use of reaction path modelling to investigate the evolution of water chemistry in shallow to deep crystalline aquifers with a special focus on fluoride

Crystalline aquifers are layered systems in which the hydrogeological path of waters extends from highly weathered, shallow and porous rocks to poorly weathered, deep and fissured rocks. This varying hydrogeological setting influences the water chemistry in different ways. The paper aims to reconstruct the water-rock interaction process in these various environments starting from a solid reactant represented by an average granite rock and several waters from the shallow aquifer. Afterwards, the water-rock interaction processes occurring in the deep environment are reconstructed, varying the geochemical conditions (primary reactants, secondary mineral phases allowed to precipitate, fO₂ and fCO₂), with a special focus on fluoride (F-). The evolution from the F-poor, Ca-HCO₃ facies to the F-rich, Na-HCO₃ water type of high pH was simulated using reaction path modelling. The obtained results show that the theoretical evolution trends well reproduce both shallow and deep water samples providing detailed information on the behavior of fluoride and other relevant constituents (i.e., Na, K, Ca, Mg, SiO₂). The performed model represents a flexible and powerful tool for environmental research, applicable in other areas hosting F-rich groundwater.

A multivariate non-parametric approach for estimating probability of exceeding the local natural background level of arsenic in the aquifers of Calabria region (Southern Italy)

The concept of natural background level (NBL) aims at distinguishing the natural and anthropogenic contributions to concentrations of specific contaminants, as groundwater management and protection tools. This is usually defined as a unique value at a regional scale, even when the hydrogeological and geochemical features of a certain territory are far from homogeneous. The concentration of target contaminants is affected by multiple hydrogeochemical processes. This is the case of arsenic in the Calabria region, where concentrations are definitely variable in groundwater. To overcome the limitation of a traditional approach and to include the intrinsic hydrogeological and geochemical heterogeneity into the definition of the natural contribution to As content in groundwater, an integrated probabilistic approach to the NBL assessment combining aquifer-based preselection criteria and multivariate non-parametric geostatistics was proposed. In detail, different NBL values were selected, based on the aquifer type and/or hydrogeochemical features. Then, these aquifer-based NBL values of arsenic were used in the Probability Kriging method to map the probability of exceedance and to provide contamination risk management tools. This multivariate geostatistical approach that takes advantage of the physico-chemical variables used in the aquifer-based NBL values definition allowed mapping the probability of exceedance of As in a physically-based way. The hydrogeochemical diversity of the study area and all the processes affecting As concentrations in the aquifers have been considered too. As a result, the obtained map was characterized by a short-range and long-range variability due to local hydrogeochemical anomalies and water-rock interaction and/or atmospheric precipitation. By this approach, the NBL exceedance probability maps proved to be less "noisy", because the local hydrogeochemical conditions were filtered, and more capable of pointing out anthropogenic inputs or very anomalous natural contributions, which need to be investigated more in detail and properly managed.

Arsenic polluted waters: Application of geochemical modelling as a tool to understand the release and fate of the pollutant in crystalline aquifers

Arsenic (As) is one of the most investigated elements worldwide due to its negative impact on the natural system. Its geochemical behavior depends on several geogenic processes, which can cause hazardous enrichment into natural waters, even in remote areas, far from anthropogenic sources. In this work the arsenic pollution issue has been addressed by studying water-rock interaction processes and applying reaction path modelling as a tool to understand the rock-to-water release of As and the fate of this natural pollutant in crystalline aquifers. In-depth geochemical characterization of several water samples discharging from crystalline aquifers was performed. The obtained data were used to fix the boundary conditions and validate the modelling outcomes. The performed modelling allowed to reconstruct the water-rock interaction processes which occur (i) in shallow and relatively shallow crystalline aquifers in which no As anomalies were observed and (ii) in As-rich areas, coupling reaction path modelling of granite dissolution with adsorption of dissolved As onto precipitating crystalline and amorphous Fe(III)-oxyhydroxides given the widespread presence of these phases in the studied environment. The results of the geochemical modelling are in agreement with the analytical data and reproduce them satisfactorily. The performed geochemical modelling is of high environmental significance because it is a flexible and powerful tool that correctly defines the water-rock interaction processes occurring in crystalline aquifers, providing valuable data to improve the knowledge on As behavior, not only in the study area, but also in similar geological settings worldwide. Therefore, the present research has broad future perspectives in the environmental field.

Mobility of REE from a hyperacid brine to secondary minerals precipitated in a volcanic hydrothermal system: Kawah Ijen crater lake (Java, Indonesia)

Rare Earth Elements (REE; lanthanides and yttrium) are elements with high economic interest because they are critical elements for modern technologies. This study mainly focuses on the geochemical behavior of REE in hyperacid sulphate brines in volcanic-hydrothermal systems, where the precipitation of sulphate minerals occurs. Kawah Ijen lake, a hyperacid brine hosted in the Ijen caldera (Indonesia), was used as natural laboratory. ∑REE concentration in the lake water is high, ranging from 5.86 to 6.52 mg kg^-1. The REE pattern of lake waters normalized to the average local volcanic rock is flat, suggesting isochemical dissolution. Minerals spontaneously precipitated in laboratory at 25 °C from water samples of Kawah Ijen were identified by XRD as gypsum. Microprobe analyses and the chemical composition of major constituents allow to identify possible other minerals precipitated: jarosite, Al-sulphate and Sr, Ba-sulphate. ∑REE concentration in minerals precipitated (mainly gypsum) range from 59.53 to 78.64 mg kg^-1. The REE patterns of minerals precipitated normalized to the average local magmatic rock show enrichment in LREE. The REE distribution coefficient (K_D), obtained from a ratio of its concentration in the minerals precipitated (mainly gypsum) and the lake water, shows higher values for LREE than HREE. K_D-LREE/K_D-HREE increases in the studied samples when the concentrations of BaO, MgO, Fe₂O₃, Al₂O₃, Na₂O and the sum of total oxides (except SO₃ and CaO) decrease in the solid phase. The presence of secondary minerals different than gypsum can be the cause of the distribution coefficient variations. High concentrations of REE in Kawah Ijen volcanic lake have to enhance the interest on these environments as possible REE reservoir, stimulating future investigations. The comparison of the K_D calculated for REE after mineral precipitation (mainly gypsum) from Kawah Ijen and Poás hyperacid volcanic lakes allow to generalize that the gypsum precipitation removes the LREE from water.

Geochemical modeling of chromium release in natural waters and treatment by RO/NF membrane processes

In this work, a geochemical approach was used as strong-scientific tool for pre-selection of suitable remediation systems to treat Cr-contaminated groundwaters. The geochemical characterization allowed to select Nanofiltration (NF) and Reverse Osmosis (RO) as suitable remediation processes, whereas through a new geochemical modeling, the evolution of water chemistry during the water-rock interaction was also studied. The new reaction path modelling was performed re-evaluating the role of Fe as main oxidant in the system and the analytic concentrations of relevant solutes, including Cr(VI), were reproduced. The spring with the highest Cr(VI) content was treated to lower its concentration below the threshold values. A laboratory-scale set-up was used to carry out both NF and RO experiments. The experiments were conducted on different commercial membranes varying the operating pressures. The results showed high Cr(VI) rejections (around 95%) for all tested membranes, leading to Cr(VI) concentrations below the threshold limits. The high flux, obtained already at lower operating pressures, combined with high selectivity towards Cr(VI) makes NF a favorable remediation option.

Release and fate of Cr(VI) in the ophiolitic aquifers of Italy: the role of Fe(III) as a potential oxidant of Cr(III) supported by reaction path modelling

Reaction path modelling of serpentinite dissolution in meteoric water was performed, varying Fe₂O₃/(FeO + Fe₂O₃) weight ratios of serpentine and reproducing the analytical concentrations of relevant solutes, including Cr(VI), in the Mg-HCO₃ groundwaters hosted in the ophiolitic aquifers of Italy. The occurrence of geogenic Cr(VI) in these groundwaters appears to be potentially controlled by the oxidation of trivalent Cr to the hexavalent redox state, driven by the reduction of trivalent Fe to the divalent redox state. In fact, trivalent Fe is the only oxidant present in suitable amounts in serpentinite rock, and even serpentine contains a high content of trivalent Fe as proven by recent studies. In contrast, the generally accepted hypothesis that geogenic Cr(VI) in waters interacting with serpentinites is driven by the reduction of trivalent and tetravalent Mn is questionable, since serpentinite rock has a low Mn content and it is necessary to invoke adsorption of trivalent Cr ions onto the surface of Mn oxides, which are oxidised as a surface complex to hexavalent Cr. Moreover, Mn oxides are present in the aqueous solution as suspended particles or in rocks as coatings or nodules.

Ultra-trace determination of total mercury in Italian bottled waters

Mercury (Hg) is a widespread, highly toxic persistent pollutant with adverse health effects on humans. So far, concentrations below the method detection limit have always been reported by studies on the concentration of mercury in bottled water when determined using instrumental analytical methods. These are often very expensive and are unaffordable for many laboratories. In this work, a less expensive method based on cold vapour atomic fluorescence spectrometry has been employed to determine total mercury (Hg_T) concentrations in bottled natural mineral waters. In all, 255 waters representing 164 different typologies were analysed. They came from 136 springs located in 18 Italian regions. In all samples, Hg_T concentrations were found in the range of sub-nanogram to a few nanograms per litre, well below the National and European regulatory limit (1 μg L^-1). Differences in Hg_T concentrations were related not only to the environmental characteristics of the springs but also to the extent and impact of human activities. Higher concentrations were found in waters coming from regions with former mining and/or natural thermal and volcanic activity. These data allowed us to estimate the mercury intake by population (adults, children and toddlers) from drinkable mineral waters consumption. The mean mercury daily intake was found to be remarkably lower, not only than the provisional tolerable value (1 μg L^-1 according to European and Italian legislation) but also than the estimated provisional tolerable weekly intake (PTWI) value (4 μg kg^-1 body weight) recommended by the Joint FAO/WHO Expert Committee on Food Additives (JECFA).

Effects of source rocks, soil features and climate on natural gamma radioactivity in the Crati valley (Calabria, Southern Italy)

The study, which represents an innovative scientific strategy to approach the study of natural radioactivity in terms of spatial and temporal variability, was aimed to characterize the background levels of natural radionuclides in soil and rock in the urban and peri-urban soil of a southern Italy area; to quantify their variations due to radionuclide bearing minerals and soil properties, taking into account nature and extent of seasonality influence. Its main novelty is taking into account the effect of climate in controlling natural gamma radioactivity as well as analysing soil radioactivity in terms of soil properties and pedogenetic processes. In different bedrocks and soils, activities of natural radionuclides ((238)U, (232)Th (4) K) and total radioactivity were measured at 181 locations by means of scintillation γ-ray spectrometry. In addition, selected rocks samples were collected and analysed, using a Scanning Electron Microscope (SEM) equipped with an Energy Dispersive Spectrometer (EDS) and an X-Ray Powder Diffraction (XRPD), to assess the main sources of radionuclides. The natural-gamma background is intimately related to differing petrologic features of crystalline source rocks and to peculiar pedogenetic features and processes. The radioactivity survey was conducted during two different seasons with marked changes in the main climatic characteristics, namely dry summer and moist winter, to evaluate possible effects of seasonal climatic variations and soil properties on radioactivity measurements. Seasonal variations of radionuclides activities show their peak values in summer. The activities of (238)U, (232)Th and (4) K exhibit a positive correlation with the air temperature and are negatively correlated with precipitations.