Meaning of the λ parameter of skew–normal and log–skew normal distributions in fluid geochemistry

The literature related to skew–normal distributions has grown rapidly in recent years but at the moment few applications concern the description of natural phenomena with this type of probability models, as well as the interpretation of their parameters. The skew–normal distributions family represents an extension of the normal family to which a parameter (λ) has been added to regulate the skewness. The development of this theoretical field has followed the general tendency in Statistics towards more flexible methods to represent features of the data, as adequately as possible, and to reduce unrealistic assumptions as the normality that underlies most methods of univariate and multivariate analysis. In this paper an investigation on the shape of the frequency distribution of the logratio ln(Cl−/Na+) whose components are related to waters composition for 26 wells, has been performed. Samples have been collected around the active center of Vulcano island (Aeolian archipelago, southern Italy) from 1977 up to now at time intervals of about six months. Data of the logratio have been tentatively modeled by evaluating the performance of the skew–normal model for each well. Values of the λ parameter have been compared by considering temperature and spatial position of the sampling points. Preliminary results indicate that changes in λ values can be related to the nature of environmental processes affecting the data

Thermodynamics and log–contrast analysis in fluid geochemistry

There are two principal chemical concepts that are important for studying the natural environment. The first one is thermodynamics, which describes whether a system is at equilibrium or can spontaneously change by chemical reactions. The second main concept is how fast chemical reactions (kinetics or rate of chemical change) take place whenever they start. In this work we examine a natural system in which both thermodynamics and kinetic factors are important in determining the abundance of NH+4 , NO−2 and NO−3 in superficial waters. Samples were collected in the Arno Basin (Tuscany, Italy), a system in which natural and antrophic effects both contribute to highly modify the chemical composition of water. Thermodynamical modelling based on the reduction-oxidation reactions involving the passage NH+4 -> NO−2 -> NO−3 in equilibrium conditions has allowed to determine the Eh redox potential values able to characterise the state of each sample and, consequently, of the fluid environment from which it was drawn. Just as pH expresses the concentration of H+ in solution, redox potential is used to express the tendency of an environment to receive or supply electrons. In this context, oxic environments, as those of river systems, are said to have a high redox potential because O2 is available as an electron acceptor. Principles of thermodynamics and chemical kinetics allow to obtain a model that often does not completely describe the reality of natural systems. Chemical reactions may indeed fail to achieve equilibrium because the products escape from the site of the rection or because reactions involving the trasformation are very slow, so that non-equilibrium conditions exist for long periods. Moreover, reaction rates can be sensitive to poorly understood catalytic effects or to surface effects, while variables as concentration (a large number of chemical species can coexist and interact concurrently), temperature and pressure can have large gradients in natural systems. By taking into account this, data of 91 water samples have been modelled by using statistical methodologies for compositional data. The application of log–contrast analysis has allowed to obtain statistical parameters to be correlated with the calculated Eh values. In this way, natural conditions in which chemical equilibrium is hypothesised, as well as underlying fast reactions, are compared with those described by a stochastic approach

A Log-linear model of grain size influence on the geochemistry of sediments

Sediment composition is mainly controlled by the nature of the source rock(s), and chemical (weathering) and physical processes (mechanical crushing, abrasion, hydrodynamic sorting) during alteration and transport. Although the factors controlling these processes are conceptually well understood, detailed quantification of compositional changes induced by a single process are rare, as are examples where the effects of several processes can be distinguished. The present study was designed to characterize the role of mechanical crushing and sorting in the absence of chemical weathering. Twenty sediment samples were taken from Alpine glaciers that erode almost pure granitoid lithologies. For each sample, 11 grain-size fractions from granules to clay (ø grades <-1 to >9) were separated, and each fraction was analysed for its chemical composition. The presence of clear steps in the box-plots of all parts (in adequate ilr and clr scales) against ø is assumed to be explained by typical crystal size ranges for the relevant mineral phases. These scatter plots and the biplot suggest a splitting of the full grain size range into three groups: coarser than ø=4 (comparatively rich in SiO2, Na2O, K2O, Al2O3, and dominated by “felsic” minerals like quartz and feldspar), finer than ø=8 (comparatively rich in TiO2, MnO, MgO, Fe2O3, mostly related to “mafic” sheet silicates like biotite and chlorite), and intermediate grains sizes (4≤ø <8; comparatively rich in P2O5 and CaO, related to apatite, some feldspar). To further test the absence of chemical weathering, the observed compositions were regressed against three explanatory variables: a trend on grain size in ø scale, a step function for ø≥4, and another for ø≥8. The original hypothesis was that the trend could be identified with weathering effects, whereas each step function would highlight those minerals with biggest characteristic size at its lower end. Results suggest that this assumption is reasonable for the step function, but that besides weathering some other factors (different mechanical behavior of minerals) have also an important contribution to the trend. Key words: sediment, geochemistry, grain size, regression, step function

Aqueous geochemistry and oxygen isotope compositions of acid mine drainage from the Río Tinto, SW Spain, highlight inconsistencies in current models

The Rio Tinto river in SW Spain is a classic example of acid mine drainage and the focus of an increasing amount of research including environmental geochemistry, extremophile microbiology and Mars-analogue studies. Its 5000-year mining legacy has resulted in a wide range of point inputs including spoil heaps and tunnels draining underground workings. The variety of inputs and importance of the river as a research site make it an ideal location for investigating sulphide oxidation mechanisms at the field scale. Mass balance calculations showed that pyrite oxidation accounts for over 93% of the dissolved sulphate derived from sulphide oxidation in the Rio Tinto point inputs. Oxygen isotopes in water and sulphate were analysed from a variety of drainage sources and displayed delta O-18((SO4-H2O)) values from 3.9 to 13.6 parts per thousand, indicating that different oxidation pathways occurred at different sites within the catchment. The most commonly used approach to interpreting field oxygen isotope data applies water and oxygen fractionation factors derived from laboratory experiments. We demonstrate that this approach cannot explain high delta O-18((SO4-H2O)) values in a manner that is consistent with recent models of pyrite and sulphoxyanion oxidation. In the Rio Tinto, high delta O-18((SO4-H2O)) values (11.2-13.6 parts per thousand) occur in concentrated (Fe = 172-829 mM), low pH (0.88-1.4), ferrous iron (68-91% of total Fe) waters and are most simply explained by a mechanism involving a dissolved sulphite intermediate, sulphite-water oxygen equilibrium exchange and finally sulphite oxidation to sulphate with O-2. In contrast, drainage from large waste blocks of acid volcanic tuff with pyritiferous veins also had low pH (1.7). but had a low delta O-18((SO4-H2O)) value of 4.0 parts per thousand and high concentrations of ferric iron (Fe(III) = 185 mM, total Fe = 186 mM), suggesting a pathway where ferric iron is the primary oxidant, water is the primary source of oxygen in the sulphate and where sulphate is released directly from the pyrite surface. However, problems remain with the sulphite-water oxygen exchange model and recommendations are therefore made for future experiments to refine our understanding of oxygen isotopes in pyrite oxidation. (C) 2009 Elsevier B.V. All rights reserved.

Geology and geochronology of the Betara region in south-southeastern Brazil: Evidence for possible Statherian (1.80-1.75 Ga) and Calymmian (1.50-1.45 Ga) extension events

The results of geological mapping, chemical analysis and radiometric dating of metabasic rocks of Betara Formation, and mapping and dating of those present in the Betara basement nucleus together with mylonitic granodiorite and syenogranite are reported here. U-Pb analysis of bulk zircon fractions from the metabasic rocks of the basement nucleus yielded a Statherian age of 1790 +/- 22 Ma, while the metabasic rocks from the upper part of the Betara Formation yielded a Calymmian age between 1500 and 1450 Ma. This age is a minimum for the deposition of the Betara Formation. The older metabasic rocks are associated with post-tectonic, possibly anorogenic syenogranite, while the younger ones are gabbro or very porphyritic ankaramite whose REE patterns are consistent with crystallization from an N-MORB parent magma. The observations and data point to the probable events associated with extensional processes of the end of Paleoproterozoic and early Mesoproterozoic. Similar registers of Statherian (1.80-1.75 Ga) and Calymmian (1.50-1.45 Ga) extensional events are recorded in other parts of the South American and African continents. The Neoproterozoic witnessed the formation and junction of the tectonic slices which formed the Apiai domain during the assemblage of western Gondwana. (C) 2010 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

Petrography, magnetic susceptibility and geochemistry of the Rio Branco Granite, Carajás Province, southeast of Pará, Brazil

Petrografia, suscetibilidade magnética e geoquímica do Granito Rio Branco, Província Carajás, sudeste do Pará, Brazil. O Granito Rio Branco é um stock paleoproterozoico intrusivo no biotita-monzogranito arqueano Cruzadão. Ocorre a oeste da cidade de Canaã dos Carajás, nas proximidades da mina de cobre do Sossego na Província Carajás. É constituído por sienogranitos não deformados e isotrópicos, hololeucocráticos, em geral de granulação média. A mineralogia é formada por feldspato alcalino pertítico, quartzo e plagioclásio. A biotita, intensamente cloritizada, é a principal fase máfica, acompanhada por flluorita, allanita, zircão, pirita e calcopirita como minerais acessórios. Albitização e, com menor intensidade greisenização, afetaram o granito, sendo a mineralogia secundária albita, fluorita, topázio, clorita, muscovita, siderofilita e óxidos e/ou hidróxidos de ferro. O Granito Rio Branco apresenta valores sistematicamente baixos de suscetibilidade magnética (SM) variando de 1,3 x 10^-5 a 6,96 x 10^-4 (SI). Geoquimicamente, é metaluminoso a peraluminoso, possui altas razões FeOt/(FeOt + MgO) e mostra afinidades com granitos ferrosos, tipo-A do subtipo A2. Os padrões dos ETR revelam um ligeiro enriquecimento de ETR leves em relação ao ETR pesados e anomalia negativa acentuada de Eu (Eu/Eu* = 0,08 - 0,13), resultando feição em "gaivota", característica de granitos evoluídos. O conjunto de dados obtidos demonstra o caráter evoluído do Granito Rio Branco e sua derivação a partir de líquidos reduzidos e enriquecidos em voláteis, causadores das transformações hidrotermais tardias. O estudo comparativo deste corpo com aqueles das suítes anorogênicas da Província Carajás sugere que o Granito Rio Branco possui maior afinidade com os granitos das suítes Velho Guilherme e, em menor grau, Serra dos Carajás. Por outro lado, é claramente distinto da Suíte Jamon. Embora apresente características similares às dos granitos especializados em estanho, não há mineralizações desta natureza associadas ao corpo.

Provenance discrimination of Lower Cretaceous synrift sandstones (eastern Iberian Chain Spain) Constraints from detrital modes heavy minerals and geochemistry

Petrography, geochemical whole-rock composition, and chemical analyses of tourmaline were performed in order to determine the source areas of Lower Cretaceous Mora, El Castellar, and uppermost Camarillas Formation sandstones from the Iberian Chain, Spain. Sandstones were deposited in intraplate subbasins, which are bound by plutonic and volcanic rocks of Permian, Triassic, and Jurassic age, Paleozoic metamorphic rocks, and Triassic sedimentary rocks. Modal analyses together with petrographic and cathodoluminescence observations allowed us to define three quartz-feldspathic petrofacies and recognize diagenetic processes that modified the original framework composition. Results from average restored petrofacies are: Mora petrofacies = P/F >1 and Q(r)70 F(r)22 R(r)9; El Castellar petrofacies = P/F >1 and Q(r)57 F(r)25 R(r)18; and Camarillas petrofacies = P/F ∼ zero and Q(r)64 F(r)28 R(r)7 (P—plagioclase; F—feldspar; Q—quartz; R—rock fragments; r—restored composition). Trace-element and rare earth element abundances of whole-rock analyses discriminate well between the three petrofacies based on: (1) the Rb concentration, which is indicative of the K content and reflects the amount of K-feldspar modal abundance, and (2) the relative modal abundance of heavy minerals (tourmaline, zircon, titanite, and apatite), which is reproduced by the elements hosted in the observed heavy mineral assemblage (i.e., B and Li for tourmaline; Zr, Hf, and Ta for zircon; Ti, Ta, Nb, and their rare earth elements for titanite; and P, Y, and their rare earth elements for apatite). Tourmaline chemical composition for the three petrofacies ranges from Fe-tourmaline of granitic to Mg-tourmaline of metamorphic origin. The three defined petrofacies suggest a mixed provenance from plutonic and metamorphic source rocks. However, a progressively major influence of granitic source rocks was detected from the lowermost Mora petrofacies toward the uppermost Camarillas petrofacies. This provenance trend is consistent with the uplift and erosion of the Iberian Massif, which coincided with the development of the latest Berriasian synrift regional unconformity and affected all of the Iberian intraplate basins. The uplifting stage of Iberian Massif pluton caused a significant dilution of Paleozoic metamorphic source areas, which were dominant during the sedimentation of the lowermost Mora and El Castellar petrofacies. The association of petrographic data with whole-rock geochemical compositions and tourmaline chemical analysis has proved to be useful for determining source area characteristics, their predominance, and the evolution of source rock types during the deposition of quartz-feldspathic sandstones in intraplate basins. This approach ensures that provenance interpretation is consistent with the geological context.

Geochemistry and Stable Isotopes of Surface Water and Groundwater in the Continental Pit in Butte, Montana, USA

The Continental porphyry Cu‐Mo mine, located 2 km east of the famous Berkeley Pit lake of Butte, Montana, contains two small lakes that vary in size depending on mining activity. In contrast to the acidic Berkeley Pit lake, the Continental Pit waters have near-neutral pH and relatively low metal concentrations. The main reason is geological: whereas the Berkeley Pit mined highly‐altered granite rich in pyrite with no neutralizing potential, the Continental Pit is mining weakly‐altered granite with lower pyrite concentrations and up to 1‐2% hydrothermal calcite. The purpose of this study was to gather and interpret information that bears on the chemistry of surface water and groundwater in the active Continental Pit. Pre‐existing chemistry data from sampling of the Continental Pit were compiled from the Montana Bureau of Mines and Geology and Montana Department of Environmental Quality records. In addition, in March of 2013, new water samples were collected from the mine’s main dewatering well, the Sarsfield well, and a nearby acidic seep (Pavilion Seep) and analyzed for trace metals and several stable isotopes, including dD and d18O of water, d13C of dissolved inorganic carbon, and d34S of dissolved sulfate. In December 2013, several soil samples were collected from the shore of the frozen pit lake and surrounding area. The soil samples were analyzed using X‐ray diffraction to determine mineral content. Based on Visual Minteq modeling, water in the Continental Pit lake is near equilibrium with a number of carbonate, sulfate, and molybdate minerals, including calcite, dolomite, rhodochrosite (MnCO3), brochantite (CuSO4·3Cu(OH)2), malachite (Cu2CO3(OH)2), hydrozincite (Zn5(CO3)2(OH)6), gypsum, and powellite (CaMoO4). The fact that these minerals are close to equilibrium suggests that they are present on the weathered mine walls and/or in the sediment of the surface water ponds. X‐Ray Diffraction (XRD) analysis of the pond “beach” sample failed to show any discrete metal‐bearing phases. One of the soil samples collected higher in the mine, near an area of active weathering of chalcocite‐rich ore, contained over 50% chalcanthite (CuSO4·5H2O). This water‐soluble copper salt is easily dissolved in water, and is probably a major source of copper to the pond and underlying groundwater system. However, concentrations of copper in the latter are probably controlled by other, less‐soluble minerals, such as brochantite or malachite. Although the acidity of the Pavilion Seep is high (~ 11 meq/L), the flow is much less than the Sarsfield Well at the current time. Thus, the pH, major and minor element chemistry in the Continental Pit lakes are buffered by calcite and other carbonate minerals. For the Continental Pit waters to become acidic, the influx of acidic seepage (e.g., Pavilion Seep) would need to increase substantially over its present volume.

Geology of the Butte Ore Deposit

Dr. John Dilles discusses the geochemistry of the porphyry Cu-Mo resource found at at Butte, Montana. The porphyry formed from dilute magmatic fluids that contained 1,000s of ppm Cu between 66 and 64 Ma, and at depths of ~8 km. The porphory is zoned from innermost Cu (Ag) ore; to mixed intermediate Cu (Zn, Ag) / Zn-MJn-Ag (Cu, Pb, Au) ore; to an outer Mn-Ag (Pb) ore that grades to barren quartz.

Weathering geochemistry and Sr-Nd fingerprints of equatorial upper Nile and Congo muds

This study investigates processes of sediment generation in equatorial central Africa. An original, complete and integrated mineralogical-geochemical database on silt-sized sediments derived from different parent rocks (basalt, granite, gneiss, metapsammite, sandstone) along the East African Rift from 5°S in Tanzania to 5°N in Sudan is presented and used to assess the incidence of diverse factors controlling sediment composition (source-rock lithology, geomorphology, hydraulic sorting, grain size, recycling), with particular emphasis on chemical weathering.

Geochronology and isotope geochemistry of Eocene dykes intruding the Ladakh Batholith

In order to determine the extent and timing of dyke formation in the Ladakh Batholith we examined about 30 mostly andesitic dykes intruding the Ladakh batholith in a ca. 50 km wide area to the west of Leh (NW India). The dykes in the east of the area trend E-NE and those in the west trend N-NW. The difference in orientation is also evident in the petrography and isotopic signatures. The eastern dykes contain corroded quartz xenocrysts and show negative ε0(Nd) and positive ε0(Sr) values, where as the western dykes do not contain quartz xenocrysts and exhibit positive ε0(Nd) and near-zero ε0(Sr) values. The variability in Sr-Nd isotopes (ε0(Nd) = 3.6 to −9.6, ε0(Sr) = 0.4 to 143) and the quartz xenocrysts can best be explained by (differing degrees of) crustal assimilation of the parent magma of the dykes. Separated minerals from five dykes were dated by 40Ar-39Ar incremental heating: amphibole ages range between 50 and 54 Ma, and one biotite dated both by Rb-Sr and by 40Ar-39Ar gave an age of 45 Ma. One dated pseudotachylyte sample attests to brittle faulting at ca. 54 Ma. The combination of structural field evidence with petrographic, isotopic and geochronological analyses demonstrates that the dykes did not form from a single, progressively differentiating magma chamber, despite having formed in the same tectonic setting around the same time, and that processes such as crustal assimilation and magma mixing/mingling also played a significant role in magma petrogenesis.