Pyrite formation driven by MSW landfill leachate in the Madrid Basin, Spain

The role of municipal solid waste (MSW) landfill leachate on the genesis of minor amounts of pyrite associated with gypsum in an otherwise predominantly evaporitic sequence was studied in geological and geochemical terms. The potential association between landfill leachate and the conditions required for bacterial reduction of sulfate and fixation of H₂S as pyrite were examined. The lithological column was generally found to contain little or no Fe. The δ³⁴S values for sulfates were consistent with previously reported data; however, the measured δ¹⁸O values were slightly higher. Sulfides disseminated in the marl/lutite exhibited higher δ³⁴S values (≈−8‰) than gypsum-coating pyrite crystals (δ³⁴S < −30‰). Dissolution of gypsum to sulfate and the supply of metabolizable organic matter and Fe required for H₂S fixation as sulfides may have originated from landfill leachate. Intermittent availability of leachate, a result of the precipitation regime, can facilitate sulfur disproportionation and lead to fractionations as high as      

The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.     

Investigation about the dependence of spectral diffuse-to-direct-beam irradiance ratio on atmospheric turbidity and solar zenith angle

Summary The modifications of the solar spectral diffuse and direct-beam irradiances as well as the diffuse-to-direct-beam ratio, E_dλ/E_bλ, as a function of the aerosol optical depth, AOD, and solar zenith angle, SZA, is investigated. The E_dλ/E_bλ ratios decrease rapidly with wavelength and exponential curves in the form E_dλ/E_bλ = aλ^−b can be fitted with a great accuracy. These curves are strongly modified by the solar spectrum distribution, which is affected by the aerosol loading, aerosol optical properties and SZA. The spectral dependence of the above E_dλ/E_bλ ratios in logarithmic coordinates does not yield a straight line, while a significant departure from the linearity is revealed. The reasons for this departure are investigated in detail and it is established that the aerosol physical properties such as single scattering albedo and size distribution along with the effect of SZA are responsible. These parameters strongly affect the scattering processes in the atmosphere and as a consequence the diffuse spectral distribution. The E_dλ/E_bλ ratio, which is an indicator of the atmospheric transmittance (King, 1979), exhibits a strong wavelength and aerosol-loading dependence. The observed differences between turbid and clear atmospheres constitute a manifestation of contrasting air properties and influence solar irradiance spectra. The present work aims at investigating the effect of atmospheric turbidity and SZA on the E_dλ/E_bλ ratio. For this reason, two distinct cases are examined: one having different atmospheric turbidity conditions but same SZA and a second having different SZAs and same atmospheric turbidity levels.     

Climate change and summer mass tourism: the case of Spanish domestic tourism

This paper investigates the impact of climate change on destination choice decisions in a context of domestic coastal tourism in Spain. Destinations are characterized in terms of travel cost and coastal ‘attractors’, such as temperature and beach-related attributes. By means of a discrete choice model based on the random utility theory, these variables are used to explain the observed pattern of interprovincial domestic trips, showing trade-offs between temperature and attractiveness in the probability of a particular destination being chosen. The model is used to investigate the impact of two climate change scenarios on the allocation of domestic tourism within Spain. The findings show that while Spain’s northern colder provinces would benefit from rising temperatures, provinces in the south would experience a decrease in the frequency of trips.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

Pennsylvanian carbonate platforms adjacent to deltaic systems in an active marine foreland basin (Escalada Fm., Cantabrian Zone,NW Spain)

The Pennsylvanian marine foreland basin of the Cantabrian Zone (NW Spain) is characterized by the unique development of kilometre‐size and hundred‐metre‐thick carbonate platforms adjacent to deltaic systems. During Moscovian time, progradational clastic wedges fed by the orogen comprised proximal alluvial conglomerates and coal‐bearing deltaic sequences to distal shelfal marine deposits associated with carbonate platforms (Escalada Fm.) and distal clay‐rich submarine slopes. A first phase of carbonate platform development (Escalada I, upper Kashirian‐lower Podolskian) reached a thickness of 400 m, nearly 50 km in width and developed a distal high‐relief margin facing a starved basin, nearly 1000‐m deep. Carbonate slope clinoforms dipped up to 30° and consisted of in situ microbial boundstone, pinching out downslope into calciturbidites, argillaceous spiculites and breccias. The second carbonate platform (Escalada II, upper Podolskian‐lower Myachkovian) developed beyond the previous platform margin, following the basinward progradation of siliciclastic deposits. Both carbonate platforms include: (1) a lower part composed of siliciclastic‐carbonate cyclothems characterized by coated‐grain and ooid grainstones; and (2) a carbonate‐dominated upper part, composed of tabular and mound‐shaped wackestone and algal‐microbial boundstone strata alternating at the decametre scale with skeletal and coated‐grain grainstone beds. Carbonate platforms initiated in distal sectors of the foreland marine shelf during transgressions, when terrigenous sediments were stored in the proximal part, and developed further during highstands of 3rd‐order sequences in a high‐subsidence context. During the falling stage and lowstand systems tracts, deltaic systems prograded across the shelf burying the carbonate platforms. Key factors involved in the development of these unique carbonate platforms in an active foreland basin are: (1) the large size of the marine shelf (approaching 200 km in width); (2) the subsidence distribution pattern across the marine shelf, decreasing from proximal shoreline to distal sectors; (3) Pennsylvanian glacio‐eustacy affecting carbonate lithofacies architecture; and (4) the environmental conditions optimal for fostering microbial and algal carbonate factories.     

Controls on the morphological evolution of embayed beaches: Morphometry versus external forcing

The morphological evolution of embayed beaches on a microtidal coast is assumed to largely respond to the degree of exposure to wave conditions, decreasing the mobility with increasing beach indentation (and vice versa). However, the number of sediment arrivals at the beach or the impact of extreme storms can modify this relationship. Here, we present an analysis of 10 embayed beaches along the Catalan coast with different morphometric and sedimentary characteristics to identify the most relevant parameters controlling the morphological evolution of these embayed beaches at the inter-annual and decadal scales. The study was mostly based on LiDAR topographic data collected from 2012 to 2017, aerial photographs from 1945 to 2021, sediment sampling and a long-term series analysis of the forcing parameters (waves, sea level, precipitation and land-use changes). The results show a net loss of volume on all the studied beaches at an inter-annual scale and a general shoreline retreat during the last few decades, suggesting the influence of common processes on the evolution of the studied beaches. Smaller pocket beaches with medium-to-high indentations are more sensitive to changes induced by local factors and show higher variability in the volume of the emerged beach and shoreline position than larger beaches. The most relevant factors influencing the evolution of the studied beaches on a decadal scale were identified as changes in sea level and the reduction in sediment inputs provided by streams due to land-use changes in the drainage basin. At the inter-annual scale, the impact of extreme events is the main factor controlling beach behaviour. These general trends can be opposite locally for beaches that receive large amounts of sediment via longshore transport from adjacent beaches.     

Long-term deforestation in NW Spain: linking the Holocene fire history to vegetation change and human activities

The Holocene fire regime is thought to have had a key role in deforestation and shrubland expansion in Galicia (NW Spain) but the contribution of past societies to vegetation burning remains poorly understood. This may be, in part, due to the fact that detailed fire records from areas in close proximity to archaeological sites are scarce. To fill this gap, we performed charcoal analysis in five colluvial soils from an archaeological area (Campo Lameiro) and compared the results to earlier studies from this area and palaeo-ecological literature from NW Spain. This analysis allowed for the reconstruction of the vegetation and fire dynamics in the area during the last ca 11 000 yrs. In the Early Holocene, Fabaceae and Betula sp. were dominant in the charcoal record. Quercus sp. started to replace these species around 10 000 cal BP, forming a deciduous forest that prevailed during the Holocene Thermal Maximum until ～5500 cal BP. Following that, several cycles of potentially fire-induced forest regression with subsequent incomplete recovery eventually led to the formation of an open landscape dominated by shrubs (Erica sp. and Fabaceae). Major episodes of forest regression were (1) ～5500–5000 cal BP, which marks the mid-Holocene cooling after the Holocene Thermal Maximum, but also the period during which agropastoral activities in NW Spain became widespread, and (2) ～2000–1500 cal BP, which corresponds roughly to the end of the Roman Warm Period and the transition from the Roman to the Germanic period. The low degree of chronological precision, which is inherent in fire history reconstructions from colluvial soils, made it impossible to distinguish climatic from human-induced fires. Nonetheless, the abundance of synanthropic pollen indicators (e.g. Plantago lanceolata and Urtica dioica) since at least ～6000 cal BP strongly suggests that humans used fire to generate and maintain pasture.     

The multi-component Hekla Ö Tephra,Iceland: a complex widespread mid-Holocene tephra layer

Large Plinian eruptions from Hekla volcano, Iceland, produce compositionally zoned tephra used as key markers in tephrochronology. However, spatial variations in chemical composition of a tephra layer may complicate its identification. An example is the 5950–6180 cal a bp Hekla Ö tephra layer, which shows compositional spread from rhyolite, dacite and andesite to basalt. In soil sections north of Hekla, the SiO₂ content of the tephra glass reaches 76 wt% in the lowest unit of the Hekla Ö deposit and decreases to 62–63 wt% in the uppermost unit. Intermingled within the whole deposit are basalt tephra grains having 46–47 wt% SiO₂. The composition of the basalt glass includes primitive basalt and a more evolved basalt (MgO >6 and <6 wt%, respectively). Together with literature data, the Hekla Ö tephra and the so-called T-Tephra/Hekla-T are most likely from contemporaneous eruptions of different vents on the Hekla volcanic system, forming a single important marker tephra (Hekla ÖT) deposited over 80% of Iceland. Identification is complicated by its spatial compositional heterogeneity, such as systematic decrease in SiO₂ content from the east to the west of Hekla volcano. Consequently, an individual tephra layer from a large explosive eruption can have different composition at different locations. © 2020 John Wiley & Sons, Ltd.