Separate lead isotope analyses of leachate and residue rock fractions: implications for metal source tracing in ore deposit studies

Separate lead isotope analyses of leachate and residue fractions are applied to a broad spectrum of rocks commonly investigated in metallogenic studies. Resulting data highlight a systematic behavior of leachate and residue fractions with respect to lead isotope compositions, which essentially depends on the mineralogical composition of the rock. Granitoid and high-grade metamorphic rocks have residue compositions virtually identical to common lead. In contrast, low-grade metasedimentary rocks may have residue compositions swamped by radiogenic lead of leach-resistant zircons. Mafic magmatic rocks have residues that are often more radiogenic than leachates, depending on the ratio of leach-refractory zircons to common lead in the residual fraction of these rocks. Separate leachate and residue analyses of source rocks provide two lead isotope end members whose mixture may represent lead with the appropriate ore fluid composition. Our leaching experiments indicate that hot acid solutions (and by inference hydrothermal fluids) preferentially leach radiogenic lead from medium- to high-grade metamorphic and granitoid rocks, whereas they preferentially leach common lead from low-grade metasedimentary and mafic magmatic rocks. The method presented in this study provides a reliable alternative to other methods (i.e., age-correction of bulk-rock compositions) for the determination of the common lead signature of felsic to intermediate magmatic rocks. This may be preferable to age-corrected bulk-rock analyses, where ages to apply for corrections of bulk-rock data are not known or where moderately to highly altered rocks must be used. Case studies of orogenic gold and MVT districts of Peru (Pataz and San Vicente, respectively) show that separate leachate and residue lead isotope analyses carried out systematically on whole rocks allow a more thorough evaluation of metal source reservoirs than does the standard method of age-corrected or uncorrected bulk-rock analyses.     

Predictive GIS-Based Model of Rockfall Activity in Mountain Cliffs

Rockfall susceptibility has been analysed in mountain cliffs of the Cantabrian Range, North Spain. The main aim of this analysis has been to build a predictive model of rockfall activity from a low number of environmental and geological variables. The rockfall activity has been quantified in a GIS. The cartographic information used shows the spatial distribution of all the recent talus screes as well as their associated source areas in the rock-slopes. The area relation At/Ar (recent talus scree polygon/source basins) in the rock slopes has been used as the rockfall activity indicator. This relation has been validated in 50 pilot rock-slopes and compared with the relation number of recent rock fragments/source basin, obtained from field work. The environmental factors causing rockfall depend on the rock slope situation, and these are: altitude and sun radiation on the rock cliff. The geological factors considered are: lithology, relative position of the main discontinuities with respect to the topographic surface and two morphologic parameters: the roughness and slope gradient. A logistic regression analysis has been applied to a population of 442 limestone and quartzite rock cliffs. The dependent variable is the rockfall activity indicator, which allows the definition of two classes of rock cliff units: low and high activity. The independent variables are altitude, sun radiation (equinox radiation, summer solstice radiation, winter solstice radiation), slope roughness, slope gradient,anisotropy and lithology. Results suggest that it is possible tobuild a valid cartographic predictive model for rockfall activity in mountain rock cliffs from a limited number of easily obtainable variables. The method is especially applicable in massive rock slopes or in regions with uniform rock mass characteristics.     

In-situ Formation of Light-Absorbing Organic Matter in Cloud Water

Current climate models seem to underestimate the flux of solar energy absorbed by the global troposphere. All of these models are constrained with the assumption that cloud droplets consist of pure water. Here we demonstrate in a simple laboratory experiment that aromatic hydroxy-acids which are found in continental fine aerosol can react with hydroxyl radicals under typical conditions prevalent in cloud water influenced by biomass burning. The reactions yield colored organic species which do absorb solar radiation. We also suggest that the products of such reactions may be humic-like substances whose presence in continental aerosol has been confirmed but their source mechanisms are still much sought after. We also attempt to give a first order estimate of the enhancement of water absorption at a visible wavelength under atmospheric conditions.     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.     

The age and thermal history of Cerro Rico de Potosi, Bolivia

 Cerro Rico de Potosi, Bolivia, is the world’s largest silver deposit and has been mined since the sixteenth century for silver, and for tin and zinc during the twentieth century, together with by-product copper and lead. The deposit consists primarily of veins that cut an altered igneous body that we interpret to be a dacitic volcanic dome and its underlying tuff ring and explosion breccia. The deposit is compositionally and thermally zoned, having a core of cassiterite, wolframite, bismuthinite, and arsenopyrite surrounded by a peripheral, lower-temperature mineral assemblage consisting principally of sphalerite, galena, lead sulfosalt, and silver minerals. The low-temperature assemblage also was superimposed on the high-temperature assemblage in response to cooling of the main hydrothermal system. Both the dacite dome and the ore fluids were derived from a larger magmatic/hydrothermal source at depth. The dome was repeatedly fractured by recurrent movement on the fault system that guided its initial emplacement. The dome was extruded at 13.8±0.2 Ma (2σ), based on U-Th-Pb dating of zircon. Mineralization and alteration occurred within about 0.3 my of dome emplacement, as indicated by a ⁴⁰Ar/³⁹Ar date of 13.76±0.10 Ma (1σ) for sericite from the pervasive quartz-sericite-pyrite alteration associated with the main-stage, high-temperature, mineralization. The last thermal event able to reset zircon fission tracks occurred no later than 12.5±1.1 Ma (1σ), as indicated by fission-tract dating. Minor sericite, and magmatic-steam alunite veins, were episodically formed around 11 Ma and between 8.3 and 5.7 Ma; the younger episodes occurring at the time of extensional fracturing at Cerro Rico and widespread volcanism in the adjacent Los Frailes volcanic field. None of these younger events appear to be significant thermal/mineralizing events; the exceptionally flat thermal release pattern of ³⁹Ar from sericite and the results of the fission-tract dating of zircon show that none of the younger events was hot enough, and lasted long enough, to cause significant loss of Ar or annealing of zircon fission tracks. U-Th-Pb dating of zircon cores indicates a Precambrian progenitor for some zircons, and REE analyses of dated samples of hydrothermally altered dacite show the presence of a prominent positive Eu anomaly, which constrains interpretations of the origin and evolution of the magmatic/hydrothermal system. Received: 14 October 1995/Accepted: 29 January 1996     

The replacement of anthophyllite by jimthompsonite: a model for hydration reactions in biopyriboles

Anthophyllite crystals found in ultramafic lenses of the Lepontine Alps (Switzerland) contain coherent, submicroscopic intergrowths of ordered and disordered biopyribole polysomes. The chain width distributions of disordered polysomes were analyzed using high resolution transmission electron microscopy (HRTEM). Chains wider than triple were interpreted as intermediate products in the transformation of anthophyllite to the triple chain silicate jimthompsonite. The concentration of individual chain types is strongly correlated with the reaction progress. Based on observed zipper terminations and the transformation rules given by Veblen and Buseck (1980) a scheme of possible reaction paths leading from anthophyllite to jimthomp sonite is proposed. The reaction scheme and a simple kinetic model for elementary reactions allow modeling of the observed chain width distributions. The model suggests that the complex reaction paths involving steps with increasing and decreasing chain width are more important in the formation of jimthompsonite than the direct transformation from anthophyllite. The wide chains (>triple) occurring as intermediate products of the multi-step paths are structurally closer to talc than jimthompsonite. The back-transformation of these wide chains to triple chains is, therefore, a strong argument that jimthompsonite is a stable phase and not only a metastable intermediate product in the transformation of anthophyllite to talc. Received: 8 July 1996 / Accepted: 13 December 1996     

Methodical approach to distribution of epilithic and drifting algae communities in a French subalpine river: Inferences on water quality assessment

In a sub-alpine river, the Asse, with an unpredictable discharge regime, chlorophyll, density and taxonomic diversity of both drifting and periphytic communities were studied from a methodical approach. The investigations based on stamping two close cross-sections to know their spatial distributions in a heterogeneous bed substrate showed a heterogeneous pattern of colonization mainly influenced by hydrological variables as the current velocity, the size of the substrata (gravel-pebble-boulder). For the drifting algae, mainly composed with benthic species temporarily unhooked to the substrata, the heterogeneity of the distribution was reduced in fast flow which induced a mixing of the whole mass of water. For the periphyton, some results indicated that the greatest diversity and the highest chlorophyll concentration were found on pebble substratum where the current velocity was moderate. These results have also shown that the classical methods used to day for investigations were inadequate for most Mediterranean rivers. It was really difficult to estimate the water quality from classical index created for use in regular large plain rivers.     

Nuclear Waste Disposal Simulations: Couplex Test Cases

We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods.     

A Finite Volume Scheme for the Transport of Radionucleides in Porous Media

The COUPLEX1 Test case (Bourgeat et al., 2003) is devoted to the comparison of numerical schemes on a convection–diffusion–reaction problem. We first show that the results of the simulation can be mainly predicted by a simple analysis of the data. A finite volume scheme, with three different treatments of the convective term, is then shown to deliver accurate and stable results under a low computational cost.