Geochemistry of crustal formation and evolution

Terrestrial planets produce crusts as they differentiate. The Earth’s bi-modal crust, with a high-standing granitic continental crust and a low-standing basaltic oceanic crust, is unique in our solar system and links the evolution of the interior and exterior of this planet. Here I present geochemical observations to constrain processes accompanying crustal formation and evolution. My approach includes geochemical analyses, quantitative modeling, and experimental studies. The Archean crustal evolution project represents my perspective on when Earth’s continental crust began forming. In this project, I utilized critical element ratios in sedimentary records to track the evolution of the MgO content in the upper continental crust as a function time. The early Archean subaerial crust had >11 wt. % MgO, whereas by the end of Archean its composition had evolved to about 4 wt. % MgO, suggesting a transition of the upper crust from a basalt-like to a more granite-like bulk composition. Driving this fundamental change of the upper crustal composition is the widespread operation of subduction processes, suggesting the onset of global plate tectonics at ~ 3 Ga (Abstract figure). Three of the chapters in this dissertation leverage the use of Eu anomalies to track the recycling of crustal materials back into the mantle, where Eu anomaly is a sensitive measure of the element’s behavior relative to neighboring lanthanoids (Sm and Gd) during crustal differentiation. My compilation of Sm-Eu-Gd data for the continental crust shows that the average crust has a net negative Eu anomaly. This result requires recycling of Eu-enriched lower continental crust to the mantle. Mass balance calculations require that about three times the mass of the modern continental crust was returned into the mantle over Earth history, possibly via density-driven recycling. High precision measurements of Eu/Eu* in selected primitive glasses of mid-ocean ridge basalt (MORB) from global MORs, combined with numerical modeling, suggests that the recycled lower crustal materials are not found within the MORB source and may have at least partially sank into the lower mantle where they can be sampled by hot spot volcanoes. The Lesser Antilles Li isotope project provides insights into the Li systematics of this young island arc, a representative section of proto-continental crust. Martinique Island lavas, to my knowledge, represent the only clear case in which crustal Li is recycled back into their mantle source, as documented by the isotopically light Li isotopes in Lesser Antilles sediments that feed into the fore arc subduction trench. By corollary, the mantle-like Li signal in global arc lavas is likely the result of broadly similar Li isotopic compositions between the upper mantle and bulk subducting sediments in most arcs. My PhD project on Li diffusion mechanism in zircon is being carried out in extensive collaboration with multiple institutes and employs analytical, experimental and modeling studies. This ongoing project, finds that REE and Y play an important role in controlling Li diffusion in natural zircons, with Li partially coupling to REE and Y to maintain charge balance. Access to state-of-art instrumentation presented critical opportunities to identify the mechanisms that cause elemental fractionation during laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) analysis. My work here elucidates the elemental fractionation associated with plasma plume condensation during laser ablation and particle-ion conversion in the ICP.

Nickel isotope fractionation during laterite Ni ore smelting and refining: implications for tracing the sources of Ni in smelter-affected soils

Nickel isotope ratios were measured in ores, fly ash, slags and FeNi samples from two metallurgical plants located in the Goiás State, Brazil (Barro Alto, Niquelândia). This allowed investigating the mass-dependent fractionation of Ni isotopes during the Ni-laterite ore smelting and refining. Feeding material exhibits a large range of δ60Ni values (from 0.02 ± 0.10 ‰ to 0.20 ± 0.05 ‰, n=7), explained by the diversity of Ni-bearing phases, and the average of δ60Nifeeding materials was found equal to 0.08 ± 0.08‰ (2SD, n=7). Both δ60Ni values of fly ash (δ60Ni = 0.07 ± 0.07‰, n=10) and final FeNi produced (0.05 ± 0.02 ‰, n=2) were not significantly different from the feeding materials ones. These values are consistent with the very high production yield of the factories. However, smelting slags present the heaviest δ60Ni values of all the smelter samples, with δ60Ni ranging from 0.11 ± 0.05 ‰ to 0.27 ± 0.05 ‰ (n=8). Soils were also collected near and far from the Niquelândia metallurgical plant, to evaluate the potential of Ni isotopes for tracing the natural vs anthropogenic Ni in soils. The Ni isotopic composition of the non-impacted topsoils developed on ultramafic rocks ranges from -0.26 ± 0.09 ‰ to -0.04 ± 0.05 ‰ (n=20). On the contrary, the Ni isotopic composition of the non-ultramafic topsoils, collected close to the plant, exhibit a large variation of δ60Ni, ranging from -0.19 ± 0.13 ‰ up to 0.10 ± 0.05 ‰ (n=4). This slight but significant enrichment in heavy isotopes highlight the potential impact of smelting activity in the surrounding area, as well as the potential of Ni isotopes for discerning anthropogenic samples (heavier δ60Ni values) from natural ones (lighter δ60Ni values). However, given the global range of published δ60Ni values (from -1.03 to 2.5 ‰) and more particularly those associated to natural weathering of ultramafic rocks (from -0.61 to 0.32‰), the use of Ni isotopes for tracing environmental contamination from smelters will remain challenging.

Geochemistry and mineralogy of a silica chimney from an inactive seafloor hydrothermal field (East Pacific Rise, 18°S)

An inactive vent field comprised of dead chimneys was discovered on the ultrafast East Pacific Rise (EPR) at 18°S during the research campaign NAUDUR with the R/V Le Nadir in December 1993. One of these chimneys was sampled, studied and found to be largely composed of silica-mineralized bacterial-like filaments. The filaments are inferred to be the result of microbial activity leading to silica (± Fe-oxyhydroxide) precipitation. The chimney grew from the most external layer (precipitated 226 ± 4 yr. B.P.) towards the central chimney conduit. Hydrothermal activity ceased 154 ± 13 yr. B.P. and the chimney conduit was completely sealed. Mixing between an end-member hydrothermal fluid and seawater explains the Sr–Nd isotopic composition of the chimney. Seawater was the major source of Sr to the chimney, whereas the dominant Nd source was the local mid-ocean ridge basalt (MORB) leached by the hydrothermal fluids. The mixing scenarios point to a dynamic hydrothermal system with fluctuating fluid compositions. The proportion of seawater within the venting fluid responsible for the precipitation of the silica chimney layers varied between 94 and 85%. Pb-isotope data indicates that all of the Pb in the chimney was derived from the underlying MORB. The precipitation temperatures of the chimney layers varied between 55 and 71 °C, and were a function of the seawater/end-member hydrothermal fluid mixing ratio. δ30Si correlates with the temperature of precipitation implying that temperature is one of the major controls of the Si-isotope composition of the chimney. Concentrations of elements across the chimney wall were a function of this mixing ratio and the composition of the end-member hydrothermal fluid. The inward growth of the chimney wall and accompanying decrease in wall permeability resulted in an inward decrease in the seawater/hydrothermal fluid mixing ratio, which in turn exerted a control on the concentrations of the elements supplied mainly by the hydrothermal fluids. The silica chimney is significantly enriched in U, likely a result of bacterial concentration of U from the seawater-dominated vent fluid. The chimney is poor in rare earth elements (REE). It inherited its REE distribution patterns from the parent end-member hydrothermal fluids. The dilution of the hydrothermal fluid with over 85% seawater could not obliterate the particular REE features (positive Eu anomaly) of the hydrothermal fluids.

Multiple sulfur isotope fractionations in inorganic aqueous systems

New constraints on isotope fractionation factors in inorganic aqueous sulfur systems based on theoretical and experimental techniques relevant to studies of the sulfur cycle in modern environments and the geologic rock record are presented in this dissertation. These include theoretical estimations of equilibrium isotope fractionation factors utilizing quantum mechanical software and a water cluster model approach for aqueous sulfur compounds that span the entire range of oxidation state for sulfur. These theoretical calculations generally reproduce the available experimental determinations from the literature and provide new constraints where no others are available. These theoretical calculations illustrate in detail the relationship between sulfur bonding environment and the mass dependence associated with equilibrium isotope exchange reactions involving all four isotopes of sulfur. I additionally highlight the effect of isomers of protonated compounds (compounds with the same chemical formula but different structure, where protons are bound to either sulfur or oxygen atoms) on isotope partitioning in the sulfite (S4+) and sulfoxylate (S2+) systems, both of which are key intermediates in oxidation-reduction processes in the sulfur cycle. I demonstrate that isomers containing the highest degree of coordination around sulfur (where protonation occurs on the sulfur atom) have a strong influence on isotopic fractionation factors, and argue that isomerization phenomenon should be considered in models of the sulfur cycle. Additionally, experimental results of the reaction rates and isotope fractionations associated with the chemical oxidation of aqueous sulfide are presented. Sulfide oxidation is a major process in the global sulfur cycle due largely to the sulfide-producing activity of anaerobic microorganisms in organic-rich marine sediments. These experiments reveal relationships between isotope fractionations and reaction rate as a function of both temperature and trace metal (ferrous iron) catalysis that I interpret in the context of the complex mechanism of sulfide oxidation. I also demonstrate that sulfide oxidation is a process associated with a mass dependence that can be described as not conforming to the mass dependence typically associated with equilibrium isotope exchange. This observation has implications for the inclusion of oxidative processes in environmental- and global-scale models of the sulfur cycle based on the mass balance of all four isotopes of sulfur. The contents of this dissertation provide key reference information on isotopic fractionation factors in aqueous sulfur systems that will have far-reaching applicability to studies of the sulfur cycle in a wide variety of natural settings.

Geochemistry Of Acid Mine Drainage From A Coal Mining Area And Processes Controlling Metal Attenuation In Stream Waters, Southern Brazil.

Acid drainage influence on the water and sediment quality was investigated in a coal mining area (southern Brazil). Mine drainage showed pH between 3.2 and 4.6 and elevated concentrations of sulfate, As and metals, of which, Fe, Mn and Zn exceeded the limits for the emission of effluents stated in the Brazilian legislation. Arsenic also exceeded the limit, but only slightly. Groundwater monitoring wells from active mines and tailings piles showed pH interval and chemical concentrations similar to those of mine drainage. However, the river and ground water samples of municipal public water supplies revealed a pH range from 7.2 to 7.5 and low chemical concentrations, although Cd concentration slightly exceeded the limit adopted by Brazilian legislation for groundwater. In general, surface waters showed large pH range (6 to 10.8), and changes caused by acid drainage in the chemical composition of these waters were not very significant. Locally, acid drainage seemed to have dissolved carbonate rocks present in the local stratigraphic sequence, attenuating the dispersion of metals and As. Stream sediments presented anomalies of these elements, which were strongly dependent on the proximity of tailings piles and abandoned mines. We found that precipitation processes in sediments and the dilution of dissolved phases were responsible for the attenuation of the concentrations of the metals and As in the acid drainage and river water mixing zone. In general, a larger influence of mining activities on the chemical composition of the surface waters and sediments was observed when enrichment factors in relation to regional background levels were used.

Rare earth element and yttrium geochemistry applied to the genetic study of cryolite ore at the Pitinga Mine (Amazon, Brazil)

This work aims at the geochemical study of Pitinga cryolite mineralization through REE and Y analyses in disseminated and massive cryolite ore deposits, as well as in fluorite occurrences. REE signatures in fluorite and cryolite are similar to those in the Madeira albite granite. The highest ΣREE values are found in magmatic cryolite (677 to 1345 ppm); ΣREE is lower in massive cryolite. Average values for the different cryolite types are 10.3 ppm, 6.66 ppm and 8.38 ppm (for nucleated, caramel and white types, respectively). Disseminated fluorite displays higher ΣREE values (1708 and 1526ppm) than fluorite in late veins(34.81ppm). Yttrium concentration is higher in disseminated fluorite and in magmatic cryolite. The evolution of several parameters (REEtotal, LREE/HREE, Y) was followed throughout successive stages of evolution in albite granites and associated mineralization. At the end of the process, late cryolite was formed with low REEtotal content. REE data indicate that the MCD was formed by, and the disseminated ore enriched by (additional formation of hydrothermal disseminated cryolite), hydrothermal fluids, residual from albite granite. The presence of tetrads is poorly defined, although nucleated, caramel and white cryolite types show evidence for tetrad effect.

Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Geodynamic evolution of the SW Variscides: Orogenic collapse shown by new tectonometamorphic and isotopic data from western Ossa-Morena Zone, SW Iberia

The pre-Mesozoic geodynamic evolution of SW Iberia has been investigated on the basis of detailed structural analysis, isotope dating, and petrologic study of high-pressure (HP) rocks, revealing the superposition of several tectonometamorphic events: (1) An HP event older than circa 358 Ma is recorded in basic rocks preserved inside marbles, which suggests subduction of a continental margin. The deformation associated with this stage is recorded by a refractory graphite fabric and noncoaxial mesoscopic structures found within the host metasediments. The sense of shear is top to south, revealing thrusting synthetic with subduction (underthrusting) to the north. (2) Recrystallization before circa 358 Ma is due to a regional-scale thermal episode and magmatism. (3) Noncoaxial deformation with top to north sense of shear in northward dipping large-scale shear zones is associated with pervasive hydration and metamorphic retrogression under mostly greenschist facies. This indicates exhumation by normal faulting in a detachment zone confined to the top to north and north dipping shear zones during postorogenic collapse soon after 358 Ma ago (inversion of earlier top to south thrusts). (4) Static recrystallization at circa 318 Ma is due to regional-scale granitic intrusions. Citation: Rosas, F. M., F. O. Marques, M. Ballevre, and C. Tassinari (2008), Geodynamic evolution of the SW Variscides: Orogenic collapse shown by new tectonometamorphic and isotopic data from western Ossa-Morena Zone, SW Iberia, Tectonics, 27, TC6008, doi:10.1029/2008TC002333.

Determination of the geographical origin of Brazilian wines by isotope and mineral analysis

In the present research, we studied wines from three different south Brazilian winemaking regions with the purpose of differentiating them by geographical origin of the grapes. Brazil`s wide territory and climate diversity allow grape cultivation and winemaking in many regions of different and unique characteristics. The wine grape cultivation for winemaking concentrates in the South Region, mainly in the Serra GaA(0)cha, the mountain area of the state of Rio Grande do Sul, which is responsible for 90% of the domestic wine production. However, in recent years, two new production regions have developed: the Campanha, the plains to the south and the Serra do Sudeste, the hills to the southeast of the state. Analysis of isotopic ratios of (18)O/(16)O of wine water, (13)C/(12)C of ethanol, and of minerals were used to characterize wines from different regions. The isotope analysis of delta(18)O of wine water and minerals Mg and Rb were the most efficient to differentiate the regions. By using isotope and mineral analysis, and discrimination analysis, it was possible to classify the wines from south Brazil.

ACCUMULATION AND TRANSLOCATION OF SILICON IN RICE AND BEAN PLANTS USING THE 30SI STABLE ISOTOPE

The 30Si silicon isotope stable was used for assessing the accumulation and translocation of Si in rice and bean plants grown in labeled nutritive solution. The isotopic silicon composition in plant materials was determined by mass spectrometry (IRMS) using the method based on SiF4 formation. Considering the total-Si added into nutritive solutions, the quantity absorbed by plants was near to 51% for rice and 15% for bean plants. The accumulated amounts of Si per plant were about 150g in rice and 8.6g in bean. Approximately 70% of the total-Si accumulated was found in leaves. At presented experimental conditions, the results confirmed that once Si is accumulated in the old parts of rice and bean plant tissues it is not redistributed to new parts, even when Si is not supplied to plants from nutritive solution.

Facies, delta(13)C, delta(15)N and C/N analyses in a late Quaternary compound estuarine fill, northern Brazil and relation to sea level

The present work integrates sedimentary facies, (14)C dating, delta(13)C, delta(15)N, and C/N with geologic and geomorphologic data available from literature. The aim was to characterize the depositional settings of a late Quaternary estuary in northeastern Marajo Island and analyze its evolution within the context of relative sea level fluctuations. The data derive from four continuous cores along a proximal-to-distal transect of a paleoestuary, previously recognized using remote sensing information. Fifteen sediment samples recorded ages ranging from 42,580 +/- 1430 to 3184 +/- 37 (14)C yr B.P. Fades analysis indicated fine- to coarse-grained sands with parallel lamination or cross stratification, massive or laminated muds and heterolithic deposits. delta(13)C (-28.1 parts per thousand to -19.7 parts per thousand, mean = -23.0 parts per thousand), delta(15)N (+ 14.8 parts per thousand to + 4.7 parts per thousand, mean = + 9.2 parts per thousand) and C/N (14.5 to 1.5, mean = 7.9) indicate mostly marine and freshwater phytoplankton sources for the organic matter. The results confirm a large late Quaternary paleoestuary in northeastern Marajo Island. The distribution of delta(13)C, delta(15)N, and C/N, together with fades associations, led to identify depositional settings related to fluvial channel, floodplain, tidal channel/tidal flat, central basin, tidal delta, and tidal inlet/sand barrier. These deposits are consistent with a wave-dominated estuary. Variations in stratigraphy and geochemistry are controlled by changes in relative sea level, revealing a main transgression from an undetermined time around 42,000 (14)C yr B.P. and 29,340 (+/- 200) (14)C yr B.P., which is synchronous to the overall drop in sea level after the last interglacial. Following this period, and probably until 9110 +/- 37 (14)C yr B.P., i.e., during a time interval encompassing two glacial episodes including the Last Glacial and the Younger Dryas, there was a pronounced drop in sea level, recorded by the development of a major erosional discontinuity due to valley re-incision. Sea level rose again until 5464 +/- 40 (14)C yr B.P, just before the main worldwide mid-Holocene transgressive peak. Mid to late Holocene coastal progradation ended the Marajo paleoestuarine history, and promoted the establishment of continental conditions throughout the island. The divergence comparing the Marajo sea level behavior with the eustatic curve allows hypothesizing that post-rifting tectonics along the Brazilian Equatorial margin influenced the sedimentary evolution of the studied paleoestuary. Considering that sedimentary facies in estuarine settings are highly variable both laterally and vertically, the present integration of facies with isotope and elemental analyses was crucial to provide a more precise interpretation of the Late Pleistocene and Holocene Marajo paleoestuary, and analyze its sea level history within the eustatic and tectonic context. (C) 2010 Elsevier B.V. All rights reserved.

Spatial variation in the strength of mutualism between a jumping spider and a terrestrial bromeliad: Evidence from the stable isotope N-15

Psecas chapoda, a neotropical jumping spider strictly associated with the terrestrial bromeliad Bromelia balansae in cerrados and semi-deciduous forests in South America, effectively contributes to plant nutrition and growth. In this study, our goal was to investigate if spider density caused spatial variations in the strength of this spider-plant mutualism. We found a positive significant relationship between spider density and delta N-15 values for bromeliad leaves in different forest fragments. Open grassland Bromeliads were associated with spiders and had higher delta N-15 values compared to forest bromeliads. Although forest bromeliads had no association with spiders their total N concentrations were higher. These results suggest that bromeliad nutrition is likely more litter-based in forests and more spider-based in open grasslands. This study is one of the few to show nutrient provisioning and conditionality in a spider-plant system. (c) 2008 Elsevier Masson SAS. All rights reserved.

Frozen Chicken for Wild Fish: Nutritional Transition in the Brazilian Amazon Region Determined by Carbon and Nitrogen Stable Isotope Ratios in Fingernails

Objectives: Amazonian populations are experiencing dietary changes characteristic of the nutrition transition. However, the degree of change appears to vary between urban and rural settings. To investigate this process, we determined carbon and nitrogen stable isotope ratios in fingernails and dietary intake of Amazonian populations living along a rural to urban continuum along the Solimoes River in Brazil. Methods: Carbon and nitrogen stable isotope ratios were analyzed from the fingernails of 431 volunteer subjects living in different settings ranging from rural villages, small towns to urban centers along the Solimoes River. Data from 200 dietary intake surveys were also collected using food frequency questionnaires and 24-h recall interviews in an effort to determine qualitative aspects of diet composition. Results: Fingernail delta(13)C values (mean standard deviation) were -23.2 +/- 1.3, 20.2 +/- 1.5, and 17.4 +/- 1.3 parts per thousand and delta(15)N values were 11.8 +/- 0.6, 10.4 +/- 0.8, and 10.8 +/- 0.7 parts per thousand for those living in rural villages, small towns, and major cities, respectively. We found a gradual increase in the number of food items derived from C(4) plant types (meat and sugar) and the replacement of food items derived from C(3) plant types (fish and manioc flour) with increasing size of urban centers. Conclusion: Increasing urbanization in the Brazilian Amazon is associated with a significant change in food habits with processed and industrialized products playing an increasingly important role in the diet and contributing to the nutrition transition in the region. Am. J. Hum. Biol. 23:642-650, 2011. (C) 2011 Wiley-Liss, Inc.

Deciphering the weathering processes using environmental mineralogy and geochemistry: Towards an integrated model of laterite and podzol genesis in the Upper Amazon Basin

This article reports major results from collaborative research between France and Brazil on soil and water systems, carried out in the Upper Amazon Basin. It reveals the weathering processes acting in the partly inundated, low elevation plateaus of the Basin, mostly covered by evergreen forest. Our findings are based on geochemical data and mineral spectroscopy that probe the crystal chemistry of Fe and Al in mineral phases (mainly kaolinite, Al- and Fe-(hydr)oxides) of tropical soils (laterites). These techniques reveal crystal alterations in mineral populations of different ages and changes of metal speciation associated with mineral or organic phases. These results provide an integrated model of soil formation and changes (from laterites to podzols) in distinct hydrological compartments of the Amazon landscapes and under altered water regimes. (C) 2010 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.