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.     

Distribution of Indus born mica along Gulf of Kachchh coast: Implications in understanding current dynamics

Indus is one of the major sources of sediments to the Gulf of Kachchh. Yet only its <63 micron fraction is studied in detail with regards to the offshore current dynamics. Hence here we present our study on characteristic signature of the Indus sediment load (i.e. mica minerals) in >63 micron size fraction along the coast of Gulf of Kachchh. The spatial distribution of mica minerals along the Gulf of Kachchh coast was studied which showed in general decreasing trend as we move along the northern and southern coast of the Gulf of Kachchh but, an increase in amount near the southern mouth at Okha. The study shows that the earlier proposed tidal barrier is ineffective in restricting movement of mica across the mouth of the gulf due to its characteristic transport mechanism. Also the presence of mudflats along the gulf of Kachchh coast plays a vital role as sediment receptors in the active sediment transport processes and mica minerals prove to be a promising simple tracer in studying the Indus born sediments in the region.     

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.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

Classifying zones of suitability for manual drilling using textural and hydraulic parameters of shallow aquifers: a case study in northwestern Senegal

A method is proposed that uses analysis of borehole stratigraphic logs for the characterization of shallow aquifers and for the assessment of areas suitable for manual drilling. The model is based on available borehole-log parameters: depth to hard rock, depth to water, thickness of laterite and hydraulic transmissivity of the shallow aquifer. The model is applied to a study area in northwestern Senegal. A dataset of boreholes logs has been processed using a software package (TANGAFRIC) developed during the research. After a manual procedure to assign a standard category describing the lithological characteristics, the next step is the automated extraction of different textural parameters and the estimation of hydraulic conductivity using reference values available in the literature. The hydraulic conductivity values estimated from stratigraphic data have been partially validated, by comparing them with measured values from a series of pumping tests carried out in large-diameter wells. The results show that this method is able to produce a reliable interpretation of the shallow hydrogeological context using information generally available in the region. The research contributes to improving the identification of areas where conditions are suitable for manual drilling. This is achieved by applying the described method, based on a structured and semi-quantitative approach, to classify the zones of suitability for given manual drilling techniques using data available in most African countries. Ultimately, this work will support proposed international programs aimed at promoting low-cost water supply in Africa and enhancing access to safe drinking water for the population.     

Space Weather: Physics,Effects and Predictability

The time varying conditions in the near-Earth space environment that may affect space-borne or ground-based technological systems and may endanger human health or life are referred to as space weather. Space weather effects arise from the dynamic and highly variable conditions in the geospace environment starting from explosive events on the Sun (solar flares), Coronal Mass Ejections near the Sun in the interplanetary medium, and various energetic effects in the magnetosphere–ionosphere–atmosphere system. As the utilization of space has become part of our everyday lives, and as our lives have become increasingly dependent on technological systems vulnerable to the space weather influences, the understanding and prediction of hazards posed by these active solar events have grown in importance. In this paper, we review the processes of the Sun–Earth interactions, the dynamic conditions within the magnetosphere, and the predictability of space weather effects on radio waves, satellites and ground-based technological systems today.     

Deep inside a neoproterozoic intra-oceanic arc: growth,differentiation and exhumation of the Amalaoulaou complex (Gourma,Mali)

We show here that the Amalaoulaou complex, in the Pan-African belt of West Africa (Gourma, Mali), corresponds to the lower and middle sections of a Neoproterozoic intra-oceanic arc. This complex records a 90–130-Ma-long evolution of magmatic inputs and differentiation above a subducting oceanic slab. Early c. 793 Ma-old metagabbros crystallised at lower crustal or uppermost mantle depths (25–30 km) and have geochemical characteristic of high-alumina basalts extracted from a depleted mantle source slightly enriched by slab-derived sedimentary components ((La/Sm)_N < 1; ε_Nd: +5.4–6.2; ⁸⁷Sr/⁸⁶Sr: 0.7027–0.7029). In response to crustal thickening, these mafic rocks were recrystallised into garnet-granulites (850–1,000°C; 10–12 kbar) and subject to local dehydration–melting reactions, forming trondhjemititic leucosomes with garnet–clinopyroxene–rutile residues. Slightly after the granulitic event, the arc root was subject to strong HT shearing during partial exhumation (detachment faults/rifting or thrusting), coeval with the emplacement of spinel- and garnet-pyroxenite dykes crystallised from a high-Mg andesitic parental magma. Quartz and hornblende-gabbros (700–660 Ma) with composition typical of hydrous volcanic rocks from mature arcs ((La/Sm)_N: 0.9–1.8; ε_Nd: +4.6 to +5.2; ⁸⁷Sr/⁸⁶Sr: 0.7028–0.7031) were subsequently emplaced at mid-arc crust levels (~15 km). Trace element and isotopic data indicate that magmas tapped a depleted mantle source significantly more enriched in oceanic sedimentary components (0.2%). Exhumation occurred either in two stages (700–660 and 623 Ma) or in one stage (623 Ma) with a final exhumation of the arc root along cold P-T path (550°C, 6–9 kbar; epidote–amphibolite and greenschist facies conditions) during the main Pan-African collision event (620–580 Ma). The composition of magmas forming the Cryogenian Amalaoulaou arc and the processes leading to intra-arc differentiation are strikingly comparable to those observed in the deep section of exposed Mezosoic oceanic arcs, namely the Kohistan and Talkeetna complex. This evolution of the Amalaoulaou oceanic arc and its accretion towards the West African craton belong to the life and closure of the Pharusian Ocean that eventually led to the formation of the Greater Gondwana supercontinent, a similar story having occurred on the other side of the Sahara with the Mozambique Ocean.     

The Amazon river breeze and the local boundary layer: I. Observations

Observations of the diurnal evolution of the planetary boundary layer over the Amazon rain forest, made at sites close to the confluence of the Solimões and Negro rivers (approximately at 3°S, 60°W) near Manaus, Amazonas, Brazil, show the existence of a diurnal rotation of the wind near the surface and the frequent presence of low-level nocturnal wind maxima. These circulations are shown to be plausibly explained as elements of a river and land breeze circulation induced by the thermal contrast between the rivers and the adjacent forest.     

Temperature dependence of the cation distribution in nickel aluminate (NiAl2O4) spinel: a powder XRD study

NiAl₂O₄ is a largely inverse spinel, which in detail shows increasing randomisation with temperature of Ni and Al between the octahedral and tetrahedral cation sites of the spinel structure. We have used powder XRD to determine this cation distribution in various samples of NiAl₂O₄ quenched after annealing between 700 and 1400° C. The inversion parameter (x) can be measured with a precision of ± 0.004 (one standard deviation), and a comparison of different methods of synthesis, X-ray diffraction and refinement techniques, suggests a probable accuracy of better than 0.01. The results are supported by some preliminary single crystal refinements on flux-grown samples.Below 800° C the rate of cation ordering becomes very slow, and, despite reaching an apparently steady state, it is doubtful if our samples attained complete internal equilibrium. Above 1250° C the cation redistribution becomes so fast that the quenching method becomes unreliable. Between 800 and 1250° C inclusive, the degree of inversion changes smoothly from 0.87 at 800° C to 0.79 at 1250° C, and is accompanied by linear changes in u, the oxygen parameter, from 0.2555 to 0.2563 (±0.0002), and a₀, the lattice parameter, from 8.0462 to 8.0522 Å (±0.0002 Å).