Trace elements contamination of agricultural soils affected by sulphide exploitation (Iberian Pyrite Belt,Sw Spain)

Agricultural soils of the Riotinto mining area (Iberian Pyrite Belt) have been studied to assess the degree of pollution by trace elements as a consequence of the extraction and treatment of sulphides. Fifteen soil samples were collected and analysed by ICP-OES and INAA for 51 elements. Chemical analyses showed an As–Cu–Pb–Zn association related with the mineralisation of the Iberian Pyrite Belt. Concentrations were 19–994 mg kg⁻¹ for As, 41–4,890 mg kg⁻¹ for Pb, 95–897 mg kg⁻¹ for Zn and of 27–1,160 mg kg⁻¹ for Cu. Most of the samples displayed concentrations of these elements higher than the 90th percentile of the corresponding geological dominium, which suggests an anthropogenic input besides the bedrock influence. Samples collected from sediments were more contaminated than leptosols because they were polluted by leachates or by mining spills coming from the waste rock piles. The weathering of the bedrock is responsible for high concentrations in Co, Cr and Ni, but an anthropogenic input, such as wind-blown dust, seems to be indicative of the high content of As, Cu, Pb and Zn in leptosols. The metal partitioning patterns show that most trace elements are associated with Fe amorphous oxy-hydroxides, or take part of the residual fraction. According to the results obtained, the following mobility sequence is proposed for major and minor elements: Mn, Pb, Cd, > Zn, Cu > Ni > As > Fe > Cr. The high mobility of Pb, Cu and Zn involve an environmental risk in this area, even in soils where the concentrations are not so high.     

Space–time distribution of ancient and active alluvial karst subsidence: examples from the central Ebro Basin,Spain

Space and time variations of alluvial karst subsidence in the central Ebro Basin are analysed in trenches (paleodolines), aerial photographs (historical dolines) and field surveys (present-day active dolines). The measured subsidence rates, as well as a comparison between densities of paleodolines computed in randomly selected scan-lines and density of active dolines computed in maps, suggest that present-day subsidence in favourable areas is more intense than that of Pleistocene times. According to diachronic maps drawn from (a) different aerial photographs taken after 1946 and (b) field surveys on deformations in urbanized areas, subsidence behaved as a nearly steady process while the whole area was used for agriculture. In contrast, rapid changes (through periods of several tens of years) have occurred in urbanized areas, in which dolines expand their boundaries through small marginal collapses, shifting sharply their sinking centres to neighbouring sites. These rapid changes do not have an equivalent in observed paleodolines.     

Environmental risk evaluation of the use of mine spoils and treated sewage sludge in the ecological restoration of limestone quarries

The ecologic restoration criteria in areas degraded from extraction activities require making use of their mine spoils. These materials do not meet fertility conditions to guarantee restoration success and therefore, need the incorporation of organic amendments to obtain efficient substratum. Reducing the deficiencies in the organic material and restoration material nutrients with the contribution of treated sewage sludge is proposed in this work. This experiment was based on a controlled study using columns. The work was conducted with two mine spoils, both very rich in calcium carbonate. The first mineral, of poor quality, came from the formation of aggregates of crushed limestone (Z). The other residual material examined originated in limestone extraction, formed by the levels of interspersed non-limestone materials and the remains of stripped soils (D). Two treatments were undertaken (30,000 and 90,000 kg/ha of sewage sludge), in addition to a control treatment. The water contribution was carried out with a device that simulated either short-duration rain or a flooding irrigation system in order to cover the surface and then percolate through the soil. The collection of leached water took place 24 h after the applications. Different parameters of the leached water were determined, including pH, electrical conductivity, nitrate anions, ammonium, phosphates, sulphates and chlorides. The values obtained for each irrigation application are discussed, and the nitrate values obtained were very elevated.     

A GPR study of subsidence-creep-fault processes in Morelia, Michoacán, Mexico

In 1983, inhabitants of the City of Morelia, Michoacán, Mexico, began to observe a series of differential settlements causing damages to constructions along linear trends parallel to a system of regional faults. The same phenomenon occurs in others cities of the Mexican Volcanic Belt (MVB), such as Celaya, Aguascalientes, and Querétaro, and is linked to a structurally controlled subsidence, caused by groundwater withdrawal, and the presence of geological faults. We define this subsidence type as Subsidence-Creep-Fault Processes (SCFP), based on the necessary elements for their generation, and we studied them through geophysical and geotechnical techniques. In Morelia, the geophysical investigations have been carried out using ground-penetrating radar (GPR). GPR profiles, perpendicular to the axis of the surface fault generated by the SCFP were carried out. The common-offset single-fold profiling was used, with a central frequency of 50 MHz. In all cases it has been possible to visualize a fault plane dividing two blocks, the presence of synthetic and antithetic faults, influence zones from 20 m to 40 m, and a maximum “net throw” of 4 m. Exploration trenches followed the same direction of the profiles obtained with GPR (perpendicular to the axis of the surface fault). These trenches exposed a fault plane dividing two blocks with different lithology, generating a maximum “net throw” of 4.40 m; as well they help in the determination of influence zones that varied from 14 m to 40 m.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Sedimentary petrology and geochemistry of siliciclastic rocks from the upper Jurassic Tordillo Formation (Neuquén Basin, western Argentina): Implications for provenance and tectonic setting

The Upper Jurassic Tordillo Formation is exposed along the western edge of the Neuquén Basin (west central Argentina) and consists of fluvial strata deposited under arid/semiarid conditions. The pebble composition of conglomerates, mineralogical composition of sandstones and pelitic rocks, and major- and trace-element geochemistry of sandstones, mudstones, and primary pyroclastic deposits are evaluated to determine the provenance and tectonic setting of the sedimentary basin. Conglomerates and sandstones derived almost exclusively from volcanic sources. The stratigraphic sections to the south show a clast population of conglomerates dominated by silicic volcanic fragments and a predominance of feldspathic litharenites. This framework composition records erosion of Triassic–Jurassic synrift volcaniclastic rocks and basement rocks from the Huincul arch, which was exhumed as a result of Late Jurassic inversion. In the northwestern part of the study area, conglomerates show a large proportion of mafic and acidic volcanic rock fragments, and sandstones are characterised by a high content of mafic volcanic rock fragments and plagioclase. These data suggest that the source of the sandstones and conglomerates was primarily the Andean magmatic arc, located west of the Neuquén Basin. The clay mineral assemblage is interpreted as the result of a complex set of factors, including source rock, climate, transport, and diagenesis. Postdepositional processes produced significant variations in the original compositions, especially the fine-grained deposits. The Tordillo sediments are characterised by moderate SiO₂ contents, variable abundances of K₂O and Na₂O, and a relatively high proportion of ferromagnesian elements. The degree of chemical weathering in the source area, expressed as the chemical index of alteration, is low to moderate. The major element geochemistry and Th/Sc, K/Rb, Co/Th, La/Sc, and Cr/Th values point to a significant input of detrital volcanic material of calcalkaline felsic and intermediate composition. However, major element geochemistry is not useful for interpreting the tectonic setting. Discrimination plots based on immobile trace elements, such as Ti, Zr, La, Sc, and Th, show that most data lie in the active continental margin field. Geochemical information is not sufficiently sensitive to differentiate the two different source areas recognized by petrographic and modal analyses of conglomerates and sandstones.     

The passivation of calcite by acid mine water. Column experiments with ferric sulfate and ferric chloride solutions at pH 2

Column experiments, simulating the behavior of passive treatment systems for acid mine drainage, have been performed. Acid solutions (HCl or H₂SO₄, pH 2), with initial concentrations of Fe(III) ranging from 250 to 1500 mg L⁻¹, were injected into column reactors packed with calcite grains at a constant flow rate. The composition of the solutions was monitored during the experiments. At the end of the experiments (passivation of the columns), the composition and structure of the solids were measured. The dissolution of calcite in the columns caused an increase in pH and the release of Ca into the solution, leading to the precipitation of gypsum and Fe–oxyhydroxysulfates (Fe(III)–SO₄–H⁺ solutions) or Fe–oxyhydroxychlorides (Fe(III)–Cl–H⁺ solutions). The columns worked as an efficient barrier for some time, increasing the pH of the circulating solutions from 2 to 6–7 and removing its metal content. However, after some time (several weeks, depending on the conditions), the columns became chemically inert. The results showed that passivation time increased with decreasing anion and metal content of the solutions. Gypsum was the phase responsible for the passivation of calcite in the experiments with Fe(III)–SO₄–H⁺ solutions. Schwertmannite and goethite appeared as the Fe(III) secondary phases in those experiments. Akaganeite was the phase responsible for the passivation of the system in the experiments with Fe(III)–Cl–H⁺ solutions.     

Coupled hydrogeological and reactive transport modelling of the Simpevarp area (Sweden)

The Simpevarp area is one of the alternative sites being considered for the deep geological disposal of high level radioactive waste in Sweden. In this paper, a coupled regional groundwater flow and reactive solute transport model of the Simpevarp area is presented that integrates current hydrogeological and hydrochemical data of the area. The model simulates the current hydrochemical pattern of the groundwater system in the area. To that aim, a conceptual hydrochemical model was developed in order to represent the dominant chemical processes. Groundwater flow conditions were reproduced by taking into account fluid-density-dependent groundwater flow and regional hydrogeologic boundary conditions. Reactive solute transport calculations were performed on the basis of the velocity field so obtained. The model was calibrated and sensitivity analyses were carried out in order to investigate the effects of heterogeneities of hydraulic conductivity in the subsurface medium. Results provided by the reactive transport model are in good agreement with much of the measured hydrochemical data. This paper emphasizes the appropriateness of the use of reactive solute transport models when water-rock interaction reactions are involved, and demonstrates what powerful tools they are for the interpretation of hydrogeological and hydrochemical data from site geological repository characterization programs, by providing a qualitative framework for data analysis and testing of conceptual assumptions in a process-oriented approach.     

Flood hazard delineation combining geomorphological and hydrological methods: an example in the Northern Iberian Peninsula

Flood mapping requires the combination and integration of geomorphological and hydrological-hydraulic methods; however, despite this, there is very little scientific literature that compares and validates both methods. Two types of analysis are addressed in the present article. On the one hand, maps of flood plains have been elaborated using geomorphological evidence and historical flood data in the mountainous area of northwestern Spain, covering an area of more then 232 km² of floodplains. On the other hand, a hydrometeorological model has been developed (Clark semidistributed unit hydrograph) in the Sarria River basin (155 km², NW Spain). This basin is not gauged, hence the model was subjected to a goodness-of-fit test of its parameter (curve number) by means of Monte Carlo simulation. The peak flows obtained by means of the hydrological model were used for hydraulic modeling (one-phase, one-dimensional and steady flow) in a 4 km² urban stretch of the river bed. The delineation of surface areas affected by floods since 1918, as well as those analyzed subsequent to the geomorphological study, reveals a high degree of reliability in the delineation of the flooded areas with frequent recurrence intervals (<50 years). If we compare these flooded surface areas with the estimate obtained by the hydrological-hydraulic method we can see that the latter method overestimates the extent of the surface water by 144% for very frequent recurrence intervals (>10 years) and underestimates it as the recurrence interval increases, by up to 80% less floodplain for exceptional events (>500 years). Finally, a management map is put forth combining the most reliable results available by integrating both methods. Originally presented at the Sixth International Conference on Geomorphology.