Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.     

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.     

Crustal segments in the North Patagonian Massif,Patagonia: An integrated perspective based on Sm–Nd isotope systematics

New insights on the Paleozoic evolution of the continental crust in the North Patagonian Massif are presented based on the analysis of Sm–Nd systematics. New evidence is presented to constrain tectonic models for the origin of Patagonia and its relations with the South American crustal blocks. Geologic, isotopic and tectonic characterization of the North Patagonian Massif and comparison of the Nd parameters lead us to conclude that: (1) The North Patagonian Massif is a crustal block with bulk crustal average ages between 2.1 and 1.6 Ga T_DM (Nd) and (2) At least three metamorphic episodes could be identified in the Paleozoic rocks of the North Patagonian Massif. In the northeastern corner, Famatinian metamorphism is widely identified. However field and petrographic evidence indicate a Middle to Late Cambrian metamorphism pre-dating the emplacement of the ca. 475 Ma granitoids. In the southwestern area, are apparent 425–420 Ma (?) and 380–360 Ma metamorphic peaks. The latter episode might have resulted from the collision of the Antonia terrane; and (3) Early Paleozoic magmatism in the northeastern area is coeval with the Famatinian arc. Nd isotopic compositions reveal that Ordovician magmatism was associated with attenuated crust. On the southwestern border, the first magmatic recycling record is Devonian. Nd data shows a step by step melting of different levels of the continental crust in the Late Palaeozoic. Between 330 and 295 Ma magmatism was likely the product of a crustal source with an average 1.5 Ga T_DM (Nd). Widespread magmatism represented by the 295–260 Ma granitoids involved a lower crustal mafic source, and continued with massive shallower-acid plutono volcanic complexes which might have recycled an upper crustal segment of the Proterozoic continental basement, resulting in a more felsic crust until the Triassic. (4) Sm–Nd parameters and detrital zircon age patterns of Early Paleozoic (meta)-sedimentary rocks from the North Patagonian Massif and those from the neighboring blocks, suggest crustal continuity between Eastern Sierras Pampeanas, southern Arequipa-Antofalla and the northeastern sector of the North Patagonian Massif by the Early Paleozoic. This evidence suggests that, at least, this corner of the North Patagonian Massif is not allochthonous to Gondwana. A Late Paleozoic frontal collision with the southwestern margin of Gondwana can be reconcilied in a para-autochthonous model including a rifting event from a similar or neighbouring position to its post-collision location. Possible Proterozoic or Early Paleozoic connections of the NPM with the Kalahari craton or the western Antartic blocks should be investigated.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

An adaptive finite element model for wellbore breakout simulation

A finite element formulation is proposed and implemented for analysing the stability of excavated wells using the DiMaggio-Sandler constitutive elastoplastic model with a typical carbonate reservoir configuration. The quality of the finite element approximation is ensured by applying smooth curved elements adapted to the wellbore geometry, and h − p adaptive finite element meshes in the plastic zone. General purpose procedures are defined to transfer the elastoplastic deformation history to newly created integration points. A breakout damage criterion is proposed based on the second invariant of the deviatoric plastic deformation tensor. This damage criterion is used to apply a mesh movement algorithm to represent material collapse. The automatic successive application of the breakout damage criterion results in elliptical realistically looking geometries obtained in experiments reported in the literature.     

Response of suspended sediment particle size distributions to changes in water chemistry at an Andean mountain stream confluence receiving arsenic rich acid drainage

Acid drainage is an important water quality issue in Andean watersheds, affecting the sustainability of urban, agricultural, and industrial activities. Mixing zones receiving acid drainage are critical sites where changes in pH and chemical environment promote the formation and dissolution of iron and aluminum oxy/hydroxides. These particles can significantly change the speciation of toxic metals and metalloids throughout drainage networks via sorption, desorption, and settling processes. However, little is known about the behavior of particle size distributions (PSDs) in streams affected by acid drainage and their relationship to metal speciation. This work studied: (a) the PSDs for a wide range of mixing ratios found at a fluvial confluence affected by acid drainage, and (b) the response of PSDs and arsenic speciation to environmental changes found when the particles approach complete mixing conditions. The confluence between the Azufre River (pH ~ 2, high concentration of dissolved metals) and Caracarani River (pH = 8.6, low concentration of dissolved metals) was used as a representative model for study. Field measurements show a bimodal PSD with modal diameters of ~50 and 300 μm. At shorter distances from the junction, the smaller modes with smaller particle volumes were dominant across the stream cross‐sections. A systematic shift towards larger particle sizes and larger particle volumes occurred downstream. The analysis of laboratory PSDs for Azufre/Caracarani mixing ratios between 0.01 and 0.5 (pH from 6.2 to 2.3) showed a bimodal trend with ~15 and 50 μm characteristic diameters; larger particles formed at pH>4. When particle suspensions were transferred in laboratory experiments from very low pH to full mixing conditions (pH ~ 2.8 and mixing ratio ~ 0.25) particle sizes varied, and the dissolved arsenic concentration decreased. The observed reaction kinetics were slow compared to the time scale of advective transport, creating opportunities for engineered controls for arsenic. This work contributes to a better understanding of the chemical‐hydrodynamic interactions in watersheds affected by mining, and identifying opportunities to improve water quality at points of use.