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.     

A GLUE‐based uncertainty assessment framework for tritium‐inferred transit time estimations under baseflow conditions

The last decade has seen major technical and scientific improvements in the study of water transfer time through catchments. Nevertheless, it has been argued that most of these developments used conservative tracers that may disregard the oldest component of water transfer, which often has transit times greater than 5 years. Indeed, although the analytical reproducibility of tracers limits the detection of the older flow components associated with the most dampened seasonal fluctuations, this is very rarely taken into account in modelling applications. Tritium is the only environmental tracer at hand to investigate transfer times in the 5‐ to 50‐year range in surface waters, as dissolved gases are not suitable due to the degassing process. Water dating with tritium has often been difficult because of the complex history of its atmospheric concentration, but its current stabilization together with recent analytical improvements open promising perspectives. In this context, the innovative contribution of this study lies in the development of a generalized likelihood uncertainty estimation‐based approach for analysing the uncertainties associated with the modelling of transit time due to both parameter identification and tracer analytical precision issues. A coupled resampling procedure allows assessment of the statistical significance of the transfer time differences found in diverse waters. This approach was developed for tritium and the exponential‐piston model but can be implemented for virtually any tracer and model. Stream baseflow, spring and shallow aquifer waters from the Vallcebre research catchments, analysed for tritium in different years with different analytical precisions, were investigated by using this approach and taking into account other sources of uncertainty. The results showed three groups of waters of different mean transit times, with all the stream baseflow and spring waters older than the 5‐year threshold needing tritium. Low sensitivity of the results to the model structure was also demonstrated. Dual solutions were found for the waters sampled in 2013, but these results may be disambiguated when additional analyses will be made in a few years. Copyright © 2016 John Wiley & Sons, Ltd.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Exploring lava‐flow hazards at Pico Island,Azores Archipelago (Portugal)

Pico, the youngest island of the Azores Archipelago (Portugal), is characterized by a central volcano and a 30‐km‐long fissure zone. Its eruption rate is the highest of the Azores islands, with more than 35 eruptions in the last 2000 years. Here, we estimate the lava‐flow hazard for Pico Island by combining the vent opening probability derived from the spatial distribution of eruptive fissures, the classes of expected eruptions inferred from the physical and chemical characteristics of historical eruptions, and the lava‐flow paths simulated by the MAGFLOW model. The most likely area to host new eruptions is along a WNW–ESE trend centred on the central volcano, with the highest hazard affecting the two main residential zones of Lajes do Pico and Madalena. Our analysis is the first attempt to assess the lava‐flow hazard for Pico Island, and may have important implications for decision‐making in territorial management and future land‐use planning.     

Significance of post‐peak metamorphic reaction microstructures in the ultrahigh temperature Eastern Ghats Province,India

Ultrahigh temperature (UHT) granulites in the Eastern Ghats Province (EGP) have a complex P–T–t history. We review the P–T histories of UHT metamorphism in the EGP and use that as a framework for investigating the P–T–t history of Mg–Al‐rich granulites from Anakapalle, with the express purpose of trying to reconcile the down‐pressure‐dominated P–T path with other UHT localities in the EGP. Mafic granulite that is host to Mg–Al‐rich metasedimentary granulites at Anakapalle has a protolith age of c. 1,580 Ma. Mg–Al‐rich metasedimentary granulites within the mafic granulite at Anakapalle were metamorphosed at UHT conditions during tectonism at 960–875 Ma, meaning that the UHT metamorphism was not the result of contact metamorphism from emplacement of the host mafic rock. Reworking occurred during the Pan‐African (c. 600–500 Ma) event, and is interpreted to have produced hydrous assemblages that overprint the post‐peak high‐T retrograde assemblages. In contrast to rocks elsewhere in the EGP that developed post‐peak cordierite, the metasedimentary granulites at Anakapalle developed post‐peak, generation ‘2’ reaction products that are cordierite‐absent and nominally anhydrous. Therefore, rocks at Anakapalle offer the unique opportunity to quantify the pressure drop that occurred during so‐called M2 that affected the EGP. We argue that M2 is either a continuation of M1 and that the overall P–T path shape is a complex counter‐clockwise loop, or that M1 is an up‐temperature counter‐clockwise deviation superimposed on the M2 path. Therefore, rather than the rocks at Anakapalle having a metamorphic history that is apparently anomalous from the rest of the EGP, we interpret that other previously studied localities in the EGP record a different part of the same P–T path history as Anakapalle, but do not preserve a significant record of pressure decrease. This is due either to the inability of refractory rocks to extensively react to produce a rich mineralogical record of pressure decrease, or because the earlier high‐P part of the rocks history was erased by the M1 loop. Irrespective of the specific scenario, models for the tectonic evolution of the EGP must take the substantial pressure decrease during M2 into account, as it is probable the P–T record at Anakapalle is a reflection of tectonics affecting the entire province.     

Detecting the stepwise propagation of the Eastern Sicily thrust belt (Italy): insight from thermal and thermochronological constraints

An extensive dataset of vitrinite reflectance, FTIR parameters on organic matter, illite content in mixed layers illite‐smectite, apatite fission tracks and U‐Th/He dating has been used to reconstruct the stepwise propagation of the Eastern Sicily fold‐and‐thrust belt during Late Palaeogene and Neogene times. The results indicate that the fold‐and‐thrust belt is divisible into two levels of thermal maturity. These levels consist of a less evolved level of thermal maturity that records limited sedimentary burial and minor heating, and a more evolved level of thermal maturity that indicates tectonic burial and exhumation at different times. Deformation and exhumation of shallowly buried units are linked to wedge forward propagation by low‐angle thrusts, whereas the evolution of deeply buried units is associated with tectonic imbrications by duplex formation and steep thrusts. The two tectonic styles alternate during evolution of the fold‐and‐thrust belt under low erosion rates.