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.     

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.     

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.     

Can Plant Competition and Diversity Reduce the Growth and Survival of Exotic <Emphasis Type="Italic">Phragmites australis</Emphasis> Invading a Tidal Marsh?

The rapid proliferation of Phragmites australis in North America has challenged resource managers to curb its expansion and reduce the loss of functional tidal marsh. We investigated whether native plant competition could reduce the ability of Phragmites to invade a tidal marsh, and if plant diversity (species richness, evenness, and composition) altered the competitive outcome. Immature Phragmites shoots and four native halophytes were transplanted to small but dense field plots (~1,200 shoots m⁻²) comprising three community structure types (Phragmites alone, Phragmites + 1 native species, and Phragmites + 4 native species). Interspecific competition significantly reduced Phragmites aboveground biomass, shoot length production, density, and survival by approximately 60%. Additionally, plots planted with greater native diversity contained Phragmites with the lowest growth and survival, potentially indicating diversity-enhanced resource competition. Competition consistently reduced the growth of Phragmites even under favorable conditions: lack of strong tidal flooding stresses as well as elevated nutrient pools.     

Evidence for the occurrence of hydrogeochemical processes in the groundwater of Khoy plain,northwestern Iran,using ionic ratios and geochemical modeling

Rapid population growth, industrialization, and agricultural expansion in the Khoy area (northwestern Iran) have led to its dependence on groundwater and degradation of groundwater quality. This study attempts to decipher the major processes and factors that degrade the groundwater quality of the Khoy plain. For this purpose, 54 groundwater samples from unconfined and confined aquifers of the plain were collected in July 2017 and analyzed for major cations and anions (Na, K, Ca, Mg, HCO₃, SO_4, and Cl), minor ions (NO₃ and F), and Al. Magnesium and bicarbonate were identified as the dominant cation and anion, respectively. Several ionic ratios and geochemical modeling using PHREEQC indicated that the most important hydrogeochemical processes to affect groundwater quality in the plain were weathering and dissolution of evaporitic and silicate minerals, mixing, and ion exchange. There were smaller effects from evaporation and anthropogenic factors (e.g., industries). Results showed that the high salinity of the groundwater in the northeast area of the plain was due to the high solubility of the evaporitic minerals, e.g., halite and gypsum. Reverse ion exchange and the contribution of mineral dissolution were more significant than ion exchange in the northeastern part of the plain. Elevated salinity of the groundwater in the southeast was attributed mostly to reverse ion exchange and somewhat to evaporation.     

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.     

Constraints on the flux of meteoritic and cometary water on the Moon from volatile element (N–Ar) analyses of single lunar soil grains,Luna 24 core

We report new nitrogen and argon isotope and abundance results for single breccia clasts and agglutinates from four different sections of the Luna 24 drill core in order to re-evaluate the provenance of N trapped in lunar regolith, and to place limits on the flux of planetary material to the Moon’s surface. Single Luna 24 grains with ⁴⁰Ar/³⁶Ar ratios <1 show δ¹⁵N values between ?54.5‰ and +123.3‰ relative to terrestrial atmosphere. Thus, low-antiquity lunar soils record both positive and negative δ¹⁵N signatures, and the secular increase of the δ¹⁵N value previously postulated by Kerridge (Kerridge, J.F. [1975]. Science 188(4184), 162–164. doi:10.1126/science.188.4184.162) is no longer apparent when the Luna and Apollo data are combined. Instead, the N isotope signatures, corrected for cosmogenic ¹⁵N, are consistent with binary mixing between isotopically light solar wind (SW) N and a planetary N component with a δ¹⁵N value of +100‰ to +160‰. The lower δ¹⁵N values of Luna 24 grains compared to Apollo samples reflect a higher relative proportion of solar N, resulting from the higher SW fluence in the region of Mare Crisium compared to the central near side of the Moon. Carbonaceous chondrite-like micro-impactors match well the required isotope characteristics of the non-solar N component trapped in low-antiquity lunar regolith. In contrast, a possible cometary contribution to the non-solar N flux is constrained to be ?3–13%. Based on the mixing ratio of SW to planetary N obtained for recently exposed lunar soils, we estimate the flux of micro-impactors to be (2.2–5.7) × 10³ tons yr^?1 at the surface of the Moon. Our estimate for Luna 24 agrees well with that for young Apollo regolith, indicating that the supply of planetary material does not depend on lunar location. Thus, the continuous influx of water-bearing cosmic dust may have represented an important source of water for the lunar surface over the past ～1 Ga, provided that water removal rates (i.e., by meteorite impacts, photodissociation, and sputtering) do not exceed accumulation rates.     

Fractional Crystallization and Liquid Immiscibility Processes in the Alkaline-Carbonatite Complex of Juqui (So Paulo, Brazil)

The Juqui circular intrusion, which is Cretaceous in age (130–135Ma), crops out in the Precambrian gneissic basement in Brazilover an area of 14 km2. It consists of olivine clinopyroxen-itecumulates (with minor olivine gabbros) in the northeastern sector(74 vol.%), whereas ijolites-melteigites-urtites (4%) and nephelinesyenites with minor essexites and syenodiorites (21%) outlinesubannular concentric patterns with an Mg-carbonatite core (1%), in the southwestern part of the complex. Petrographical, bulk rock, and mineral compositional trendsindicate that the origin of the complex can be largely accountedfor by shallow-level fractional crystallization of a carbonatedbasanitic parental magma. Such a magma was generated deep inthe subcontinental lithosphere by low-degree partial meltingof a garnet-phlogopite peridotite source. Mass-balance calculations in agreement with field volume estimatespermit definition of several fractionation stages of the magmaticevolution under nearly closed-system conditions, with inwarddevelopment of zonally arranged side-wall cumulates. These stagesinvolved: (1) fractionation from basanite to essexite magma(liquid fraction F = 33–5%) by crystallization of olivineclinopyroxenite plus minor olivine alkali gabbro cumulates;(2) derivation of the least differentiated mafic nepheline syenite(F = 5–5 %) from essexitic magma by subtraction of a syenodioriteassemblage; (3) exsolution of a carbonatite liquid (5%) froma CO₂-enriched mafic nepheline syenite magma, which also underwentcontinuous fractionation giving rise to ijolite-melteigite-urtitecumulates. The proportion of cumulus clinopyroxene and biotiteand intercumulus nepheline and alkali feldspar in these lastrocks, as well as the absence of alkalis in carbonatite, maybe attributed, at least in part, to loss of alkali-rich hydrousfluids released during and after the unmixing formation of thetwo conjugate liquids. The KD values determined for Mg-carbonatite/nepheline syeniteare lower (1–4–2–9) for light rare earth elements(LREE) than for REE from Eu to Yb (4–6–7–8),in contrast to recent experimental results (Hamilton et al.,1989). A possible explanation is that Juquia Mg-carbonatiterepresents an alreadydifferentiated magma, which underwent extensivefractionation of LREE-enriched calcite. In this way, the highvariability of K₀ REE patterns observed in several alkaline-carbonatitecomplexes can also be accounted for. The remarkably constant initial ⁸⁷Sr/⁸⁶Sr ratios (mostly between0–7052 and 0–7057) support the interpretation ofthe intrusion as having been generated by fractrional crystallizationand liquid immiscibility from a common parental magma. Iligherisotopic ratios (0–7060–0–7078), found mainlyin dykes and in the border facies of the intrusion, may be dueto contamination by the gecissic basement.     

A multi-equation spatial econometric model,with application to EU manufacturing productivity growth

A multi-equation spatial econometric model is used to explain variations across EU regions in manufacturing productivity growth based on recent theoretical developments in urban economics and economic geography. The paper shows that temporal and spatial parameter homogeneity is an unrealistic assumption, contrary to what is typically assumed in the literature. Constraints are imposed on parameters across time periods and between core and peripheral regions of the EU, with the significant loss of fit providing overwhelming evidence of parameter heterogeneity, although the final model does highlight increasing returns to scale, which is a central feature of contemporary theory.