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.     

Diurnal Variation of Gravity Wave Activity at Midlatitudes in the Ionospheric F Region

We investigate the short-term fluctuations in the period range from 15 to 180 minutes in the electron density variations of the F region ionosphere. Electron density profiles obtained at the ionospheric stations of Pruhonice (49.9° N, 14.5° E) and Ebro (40.8° N, 0.5° E) at five minute time sampling have been used for this analysis. The diurnal changes of the activity of the acoustic gravity wave fluctuations (AGW) show a clear enhancement during and several hours after sunrise. The periods of such AGW's are about 60 to 75 minutes and these waves propagates vertically through the ionosphere from a source located at an altitude of 180-220 km. The most likely source for these events seems to be passage of the Solar terminator.     

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.     

Martian water vapor: Mars Express PFS/LW observations

We present the seasonal and geographical variations of the martian water vapor monitored from the Planetary Fourier Spectrometer Long Wavelength Channel aboard the Mars Express spacecraft. Our dataset covers one martian year (end of Mars Year 26, Mars Year 27), but the seasonal coverage is far from complete. The seasonal and latitudinal behavior of the water vapor is globally consistent with previous datasets, Viking Orbiter Mars Atmospheric Water Detectors (MAWD) and Mars Global Surveyor Thermal Emission Spectrometer (MGS/TES), and with simultaneous results obtained from other Mars Express instruments, OMEGA and SPICAM. However, our absolute water columns are lower and higher by a factor of 1.5 than the values obtained by TES and SPICAM, respectively. In particular, we retrieve a Northern midsummer maximum of 60 pr-μm, lower than the 100-pr-μm observed by TES. The geographical distribution of water exhibits two local maxima at low latitudes, located over Tharsis and Arabia. Global Climate Model (GCM) simulations suggest that these local enhancements are controlled by atmospheric dynamics. During Northern spring, we observe a bulge of water vapor over the seasonal polar cap edge, consistent with the northward transport of water from the retreating seasonal cap to the permanent polar cap. In terms of vertical distribution, we find that the water volume mixing ratio over the large volcanos remains constant with the surface altitude within a factor of two. However, on the whole dataset we find that the water column, normalized to a fixed pressure, is anti-correlated with the surface pressure, indicating a vertical distribution intermediate between control by atmospheric saturation and confinement to a surface layer. This anti-correlation is not reproduced by GCM simulations of the water cycle, which do not include exchange between atmospheric and subsurface water. This situation suggests a possible role for regolith-atmosphere exchange in the martian water cycle.     

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.     

On the sampling of airborne measurements of radiative properties of cloud fields: Results from International Cirrus Experiment (ICE), 1987

Summary During the first pilot phase of the International Cirrus Experiment (ICE) in September/October 1987 radiation measurements were carried out above and below broken cloud fields over the German Bight. Analyses of the measurements of 25.9.1987 (broken cumulus cloud field) showed an unexpected high absorption of solar radiation but not those of 28.9.1987 (stratocumulus/stratus cloud field). This might be explained by a partial negligence of cloud side factor under the viewing geometry of the pyranometer and other reasons. The light guide effect (leakage) within holes can account for the depletion of cloud reflectance and the enhancement of cloud transmittance at large cloud cover, as observed during the measurements on 28.9.1987. The fluctuation of the irradiances, which can be characterized by the intermittency parameter, has been related to the corresponding cloud cover by a regression equation. On the condition of similarity, this scaling parameter derived from sample measurements could be extrapolated to that of the whole cloud field. Thus, the estimation of the irradiance spectrum for the whole cloud field may be possible, from samples taken with aircraft and imaging data.With 5 Figures     

Cosmic‐ray rejection for single spectroscopic CCD images by means of template matching

This article describes a novel algorithm for cosmic‐ray rejection in single spectroscopic CCD images. This algorithm is based on a variation of template matching. It focuses on identifying those pixels belong to spectra, while other conventional algorithms tried to locate the cosmic‐ray hits directly. The main principle is applying template matching to find suspicious blocks, which is followed by surface patching to locate the legitimate pixels accurately. Therefore, the rest pixels are the ones corrupted by cosmic‐ray hits. Meanwhile, almost all the parameters are automatically extracted from the images. Examples of its performance are given for both simulated and observed images. It shows an advantage of significantly low false alarm rate with relatively high detection rate (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.