Acid mine drainage at Cerro Rico de Potosí I: unabated high-strength discharges reflect a five century legacy of mining

Intensive mining and processing of Ag, Sn, Pb and Zn ores have occurred in various locations within and around the city of Potosí, Bolivia since 1545. Surface and subsurface waters, stream sediments and soils are contaminated with various ecotoxic metals in the headwaters of the economically vital, yet highly impacted, upper Rio Pilcomayo watershed. Previous studies have documented downstream trace metal contamination, however, not addressed specific sources. The AMD discharges identified in this study help link downstream pollution to primary origins. The majority of AMD would be considered high-strength due to metal concentrations and acidity orders of magnitude greater than typical AMD. Discharges from both operating and abandoned portals as well as tailings-related deposits displayed a high degree of heterogeneity with total metal concentrations ranging from 0.11 to 7480, <0.022 to 889, <0.0006 to 65.3, <0.001 to 310, 0.12 to 72,100, 0.3 to 402, <0.012 to 34.8, and 0.24 to 19,600 mg/L of Al, As, Cd, Cu, Fe, Mn, Pb and Zn, respectively. Net acidity and pH ranged from −10 to 246,000 mg/L as CaCO₃ equivalent and 0.90–6.94 standard units, respectively. Data were gathered during two sampling events centered around the most extreme periods of the dry and wet seasons of one water-year. Loadings to local streams were marginally greater for most metals in the wet season. If observed loadings are historically representative, Cerro Rico AMD has contributed thousands of tonnes of ecotoxic metals to the upper Rio Pilcomayo over the last five centuries. Metal and hydrogen ion concentrations in the majority of AMD sampled were several orders of magnitude above discharge limits set by the Bolivian government, yet no action has historically or contemporarily been taken.     

Anthropogenic relief changes in a long-lasting lignite mining area (‘Ville’, Germany) derived from historic maps and digital elevation models

Humans constitute one of the main geomorphological agents in modern times. As an example, post-mining regions represent a typical landscape of the Anthropocene. Strong relief modifications are particularly obvious with open pit mining. However, many existing mining areas are lacking detailed pre-mining information for the quantification of anthropogenic relief changes, which is a considerable challenge in regions with historic mining activities. Here, the Ville (Rhenish lignite district, Germany) is used to quantify surface mining induced relief changes in one of the oldest and currently largest lignite districts in Europe. Historical maps from first geodetic mapping in 1893 enabled construction of a historic digital elevation model to quantify the relief changes in comparison to elevation data from 2000 and 2015. The vertical accuracy of the historic data is remarkably high, with relief differences < 2 m in areas not affected by mining. In total, 49.2% of the investigated area (184 km²) shows a relief deficit and 14.5% has positive relief differences. Absolute changes account for more than 80 m heightening (dumpsites of overburden) and lowering of the natural relief (pits). Besides these altitudinal changes, overall steeper slopes are significant for the new topography, but levelling exists likewise. The spatial variabilities are discussed in the context of the regional geology and the mining techniques. Undoubtedly, such large-scale anthropogenic relief changes persist for a very long time and will last as a human legacy far into the future. Only the detailed reconstruction of the pre-mining relief offers the ability to clarify the dimension of humans as geomorphological agents and to understand landscape perception. Due to the fact that the impact of open pit mining has such a large vertical and horizontal extension, their consideration as part of anthropogeomorphology can significantly contribute to support future Critical Zone research in the Anthropocene.     

Thermodynamic controls on element partitioning between titanomagnetite and andesitic–dacitic silicate melts

Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO₂) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg²⁺, Mn²⁺, Co²⁺, Ni²⁺, Zn²⁺) and Cu²⁺/Cu⁺ with increasing fO₂ between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO₂ implies a reciprocal relationship between a(Fe²⁺O) and a(Fe³⁺O_1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO₂ in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO₂ in silicic, but less so in mafic bulk systems.     

Estimating coseismic deformations from near source strong motion records: methods and case studies

Digital strong-motion accelerographs have opened up the possibility of extracting ground motion characteristics at much lower frequencies than was offered by analogue instruments. High-quality digital data obtained close to the faults have tempted several efforts to retrieve permanent ground displacements after an earthquake. Such attempts have been partly successful, and somewhat subjective, the main reason being the presence of baseline offsets in the accelerometric data. We review existing methods for such applications, discuss their limitations and propose a more objective and improved scheme to make baseline adjustments and obtain permanent displacements. The proposed technique is applied to 26 digital recordings from the 29 May 2008 Ölfus Earthquake in South Iceland and 9 recordings from the 1999 Chi-Chi, Taiwan, Earthquake, and the permanent displacements obtained are compared with published results and GPS measurements from nearby stations. Our case studies show that the proposed technique, in addition to being simple and objective, is effective in making adjustments for baseline errors in accelerometric data.     

Structure Function Analysis of Two-Scale Scalar Ramps. Part II: Ramp Characteristics and Surface Renewal Flux Estimation

Ramp features in the turbulent scalar field are associated with turbulent coherent structures, which dominate energy and mass fluxes in the atmospheric surface layer. Although finer scale ramp-like shapes embedded within larger scale ramp-like shapes can readily be perceived in turbulent scalar traces, their presence has largely been overlooked in the literature. We demonstrate the signature of more than one ramp scale in structure functions of the turbulent scalar field measured from above bare ground and two types of short plant canopies, using structure-function time lags ranging in scale from isotropic to larger than the characteristic coherent structures. Spectral analysis of structure functions was used to characterize different scales of turbulent structures. By expanding structure function analysis to include two ramp scales, we characterized the intermittency, duration, and surface renewal flux contribution of the smallest (i.e., Scale One) and the dominant (i.e., Scale Two) coherent structure scales. The frequencies of the coherent structure scales increase with mean wind shear, implying that both Scale One and Scale Two are shear-driven. The embedded Scale One turbulent structure scale is ineffectual in the surface-layer energy and mass transport process. The new method reported here for obtaining surface renewal-based scalar exchange works well over bare ground and short canopies under unstable conditions, effectively eliminating the α calibration for these conditions and forming the foundation for analysis over taller and more complex surfaces.     

Application of sodars in the study and monitoring of the environment

Summary Application of acoustic sounders (sodars) to study the atmospheric boundary layer (ABL) began in the early 1970s. During the last two decades the scope of sodars applications enlarged considerably. The proceedings of eleven symposia of the International Society for Remote Acoustic Sensing of the Atmosphere and Oceans (ISARS) are a unique collection of papers, where all directions of the sodar use since 1981 are presented. In this paper, a review of sodar applications to atmospheric research is presented based on materials published in these proceedings in the following fields: conditions of microwave and light propagation; regional climatology of the ABL; air pollution meteorology and weather forecast; mesoscale phenomena under stable and unstable stratification; micrometeorology; peculiarities of the ABL in remote and complex terrain.     

Quantification of the environmental structural risk with spoiling ties: is randomization worthwhile?

Many recent works show that copulas turn out to be useful in a variety of different applications, especially in environmental sciences. Here the variables of interest are usually continuous, being times, lengths, weights, and so on. Unfortunately, the corresponding observations may suffer from (instrumental) adjustments and truncations, and eventually may show several repeated values (i.e., ties). In turn, on the one hand, a tricky issue of identifiability of the model arises, and, on the other hand, the assessment of the risk may be adversely affected. A possible remedy is to adopt suitable randomization procedures: here three different strategies are outlined. The goal of the work is to carry out a simulation study in order to evaluate the effects of the randomization of multivariate observations when ties are present. In particular, it is investigated whether, how, and to what extent, the randomization may change the estimation of the structural risk: for this purpose, a coastal engineering example will be used, as archetypical of a broad class of models and problems in engineering applications. Practical advices and warnings about the use of randomization techniques are hence given.     

Marymia: an Archean, amphibolite facies-hosted, orogenic lode-gold deposit overprinted by Palaeoproterozoic orogenesis and base metal mineralisation, Western Australia

The Marymia gold deposit, comprising two orebodies, Keillor 1 and Keillor 2, is at the northern end of the Plutonic Well greenstone belt in the Marymia Inlier, in the southern Capricorn Orogen, just north of the Yilgarn craton. The Marymia Inlier is a discrete fault-bounded Archean gneiss-granitoid-greenstone domain surrounded by sedimentary basins that were formed and variably metamorphosed and deformed during several Palaeoproterozoic orogenic cycles. The greenstone sequence at Marymia is stratigraphically and geochemically similar to greenstone sequences in the Yilgarn craton, but was subjected to further deformation and metamorphism in the Palaeoproterozoic. Late Archean deformation (D1-D2) was ductile to brittle-ductile in style, whereas Palaeoproterozoic deformation was predominantly brittle. Equilibrium mineral assemblages indicate that peak amphibolite-facies metamorphism (540-575 °C, <3 kb) was overprinted by greenschist-facies metamorphism (300-360 °C). Petrographic textures indicate that prograde metamorphism was coeval with D1-D2, with peak metamorphism early to syn D2. Gold mineralisation at Marymia is hosted in metamorphosed tholeiitic basalts and banded iron formation. On a gross scale, the distribution of gold is controlled by D2 folds and shear zones. Lithological contacts with strong rheological or chemical contrasts provide local controls. Gold-related alteration comprises subtle millimetre- to centimetre-wide zones of silicification with variable amounts of quartz, hornblende, biotite, K-feldspar, plagioclase, calcite/siderite, scheelite, titanite, epidote, sulfide and telluride minerals. Quartz veins are generally narrow and discontinuous with low total volume of quartz. Gold is sited in the wall rock, at vein salvedges or within stringers of wall rock within veins. There are two distinct opaque-mineral assemblages: pyrite-pyrrhotite-chalcopyrite-galena and hessite-petzite-altaite-Bi-telluride-galena. Ore samples are variably enriched in Ag, Te, Pb, W, Cu, S and Fe reflecting heterogeneity of the ore mineralogy. Structural timing and temperature of formation of alteration and ore minerals support deposition of gold during late peak amphibolite-facies metamorphism from neutral to alkaline (pH=5-6), moderately oxidising (log P_O2,-21-22) and CO₂-bearing (X_CO2 Ƹ.2) fluids. The total sulfur content of the fluid is estimated at 1mDS. Lead isotope compositions support derivation of lead from within the local greenstone sequence. Gold lodes were deformed by faults and shear zones in the Palaeoproterozoic, with only limited remobilisation. Subeconomic, carbonate vein- and breccia-hosted base metal mineralisation is locally hosted within Palaeoproterozoic fault zones, which clearly cut gold lodes. Base-metal-related alteration is characterised by intense carbonatisation, chloritisation, and albitisation of the mafic host rocks. Mineral assemblages are consistent with formation at greenschist facies conditions. Lead isotope compositions support crystallisation at ca. 1.7 Ga from lead that is similar in composition to earlier gold-related galena.