Lithium isotope constraints on crust–mantle interactions and surface processes on Mars

Lithium abundances and isotope compositions are reported for a suite of martian meteorites that span the range of petrological and geochemical types recognized to date for Mars. Samples include twenty-one bulk-rock enriched, intermediate and depleted shergottites, six nakhlites, two chassignites, the orthopyroxenite Allan Hills (ALH) 84001 and the polymict breccia Northwest Africa (NWA) 7034. Shergottites unaffected by terrestrial weathering exhibit a range in δ7Li from 2.1 to 6.2‰, similar to that reported for pristine terrestrial peridotites and unaltered mid-ocean ridge and ocean island basalts. Two chassignites have δ7Li values (4.0‰) intermediate to the shergottite range, and combined, these meteorites provide the most robust current constraints on δ7Li of the martian mantle. The polymict breccia NWA 7034 has the lowest δ7Li (−0.2‰) of all terrestrially unaltered martian meteorites measured to date and may represent an isotopically light surface end-member. The new data for NWA 7034 imply that martian crustal surface materials had both a lighter Li isotope composition and elevated Li abundance compared with their associated mantle. These findings are supported by Li data for olivine-phyric shergotitte NWA 1068, a black glass phase isolated from the Tissint meteorite fall, and some nakhlites, which all show evidence for assimilation of a low-δ7Li crustal component. The range in δ7Li for nakhlites (1.8 to 5.2‰), and co-variations with chlorine abundance, suggests crustal contamination by Cl-rich brines. The differences in Li isotope composition and abundance between the martian mantle and estimated crust are not as large as the fractionations observed for terrestrial continental crust and mantle, suggesting a difference in the styles of alteration and weathering between water-dominated processes on Earth versus possibly Cl–S-rich brines on Mars. Using high-MgO shergottites (>15 wt.% MgO) it is possible to estimate the δ7Li of Bulk Silicate Mars (BSM) to be 4.2 ± 0.9‰ (2σ). This value is at the higher end of estimates for the Bulk Silicate Earth (BSE; 3.5 ± 1.0‰, 2σ), but overlaps within uncertainty.

Continuous eclogite melting and variable refertilisation in upwelling heterogeneous mantle

Large-scale tectonic processes introduce a range of crustal lithologies into the Earth's mantle. These lithologies have been implicated as sources of compositional heterogeneity in mantle-derived magmas. The model being explored here assumes the presence of widely dispersed fragments of residual eclogite (derived from recycled oceanic crust), stretched and stirred by convection in the mantle. Here we show with an experimental study that these residual eclogites continuously melt during upwelling of such heterogeneous mantle and we characterize the melting reactions and compositional changes in the residue minerals. The chemical exchange between these partial melts and more refractory peridotite leads to a variably metasomatised mantle. Re-melting of these metasomatised peridotite lithologies at given pressures and temperatures results in diverse melt compositions, which may contribute to the observed heterogeneity of oceanic basalt suites. We also show that heterogeneous upwelling mantle is subject to diverse local freezing, hybridization and carbonate-carbon-silicate redox reactions along a mantle adiabat.

Benchmark problems for reactive transport modeling of the generation and attenuation of acid rock drainage

Acid rock drainage (ARD) is a problem of international relevance with substantial environmental and economic implications. Reactive transport modeling has proven a powerful tool for the process-based assessment of metal release and attenuation at ARD sites. Although a variety of models has been used to investigate ARD, a systematic model intercomparison has not been conducted to date. This contribution presents such a model intercomparison involving three synthetic benchmark problems designed to evaluate model results for the most relevant processes at ARD sites. The first benchmark (ARD-B1) focuses on the oxidation of sulfide minerals in an unsaturated tailing impoundment, affected by the ingress of atmospheric oxygen. ARD-B2 extends the first problem to include pH buffering by primary mineral dissolution and secondary mineral precipitation. The third problem (ARD-B3) in addition considers the kinetic and pH-dependent dissolution of silicate minerals under low pH conditions. The set of benchmarks was solved by four reactive transport codes, namely CrunchFlow, Flotran, HP1, and MIN3P. The results comparison focused on spatial profiles of dissolved concentrations, pH and pE, pore gas composition, and mineral assemblages. In addition, results of transient profiles for selected elements and cumulative mass loadings were considered in the intercomparison. Despite substantial differences in model formulations, very good agreement was obtained between the various codes. Residual deviations between the results are analyzed and discussed in terms of their implications for capturing system evolution and long-term mass loading predictions.

Experimental Characterization and Quantification of Cement-Bentonite Interaction using Core Infiltration Techniques Coupled with X-ray Tomography

Deep geological storage of radioactive waste foresees cementitious materials as reinforcement of tunnels and as backfill. Bentonite is proposed to enclose spent fuel canisters and as drift seals. Sand/bentonite (s/b) is foreseen as backfill material of access galleries or as drift seals. The emplacement of cementitious material next to clay material generates an enormous chemical gradient in pore-water composition that drives diffusive solute transport. Laboratory studies and reactive transport modeling predicted significant mineral alteration at and near interfaces, mainly resulting in a decrease of porosity in bentonite. The goal of this thesis was to characterize and quantify the cement/bentonite interactions both spatially and temporally in laboratory experiments. A newly developed mobile X-ray transparent core infiltration device was used to perform X-ray computed tomography (CT) scans without interruption of running experiments. CT scans allowed tracking the evolution of the reaction plume and changes in core volume/diameter/density during the experiments. In total 4 core infiltration experiments were carried out for this study with the compacted and saturated cores consisting of MX-80 bentonite and sand/MX-80 bentonite mixture (s/b; 65/35%). Two different high-pH cementitious pore-fluids were infiltrated: a young (early) ordinary Portland cement pore-fluid (APWOPC; K+–Na+–OH-; pH 13.4; ionic strength 0.28 mol/kg) and a young ‘low-pH’ ESDRED shotcrete pore-fluid (APWESDRED; Ca2+–Na+–K+–formate; pH 11.4; ionic strength 0.11 mol/kg). The experiments lasted between 1 and 2 years. In both bentonite experiments, the hydraulic conductivity was strongly reduced after switching to high-pH fluids, changing eventually from an advective to a diffusion-dominated transport regime. The reduction was mainly induced by mineral precipitation and possibly partly also by high ionic strength pore-fluids. Both bentonite cores showed a volume reduction and a resulting transient flow in which pore-water was squeezed out during high-pH infiltration. The outflow chemistry was characterized by a high ionic strength, while chloride in the initial pore water got replaced as main anionic charge carrier by sulfate, originating from gypsum dissolution. The chemistry of the high-pH fluids got strongly buffered by the bentonite, consuming hydroxide and in case of APWESDRED also formate. Hydroxide got consumed by mineral reactions (saponite and possibly talc and brucite precipitation), while formate being affected by bacterial degradation. Post-mortem analysis showed reaction zones near the inlet of the bentonite core, characterized by calcium and magnesium enrichment, consisting predominately of calcite and saponite, respectively. Silica got enriched in the outflow, indicating dissolution of silicate-minerals, identified as preferentially cristobalite. In s/b, infiltration of APWOPC reduced the hydraulic conductivity strongly, while APWESDRED infiltration had no effect. The reduction was mainly induced by mineral precipitation and probably partly also by high ionic strength pore-fluids. Not clear is why the observed mineral precipitates in the APWESDRED experiment had no effect on the fluid flow. Both s/b cores showed a volume expansion along with decreasing ionic strengths of the outflow, due to mineral reactions or in case of APWESDRED infiltration also mediated by microbiological activity, consuming hydroxide and formate, respectively. The chemistry of the high-pH fluids got strongly buffered by the s/b. In the case of APWESDRED infiltration, formate reached the outflow only for a short time, followed by enrichment in acetate, indicating most likely biological activity. This was in agreement to post-mortem analysis of the core, observing black spots on the inflow surface, while the sample had a rotten-egg smell indicative of some sulfate reduction. Post-mortem analysis showed further in both cores a Ca-enrichment in the first 10 mm of the core due to calcite precipitation. Mg-enrichment was only observed in the APWOPC experiment, originating from newly formed saponite. Silica got enriched in the outflow of both experiments, indicating dissolution of silicate-minerals, identified in the OPC experiment as cristobalite. The experiments attested an effective buffering capacity for bentonite and s/b, a progressing coupled hydraulic-chemical sealing process and also the preservation of the physical integrity of the interface region in this setup with a total pressure boundary condition on the core sample. No complete pore-clogging was observed but the hydraulic conductivity got rather strongly reduced in 3 experiments, explained by clogging of the intergranular porosity (macroporosity). Such a drop in hydraulic conductivity may impact the saturation time of the buffer in a nuclear waste repository, although the processes and geometry will be more complex in repository situation.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

The supergene thorium and rare-earth element deposit at Morro do Ferro, Poços de Caldas, Minas Gerais, Brazil

The thorium and rare-earth element (Th-REE) deposit at Morro do Ferro formed under supergene lateritic weathering conditions. The ore body consists of shallow NW-SE elongated argillaceous lenses that extend from the top of the hill downwards along its south-eastern slope. The deposit is capped by a network of magnetite layers which protected the underlying highly weathered, argillaceous host rock from excessive erosion. The surrounding country rocks comprise a sequence of subvolcanic phonolite intrusions that have been strongly altered by hydrothermal and supergene processes. From petrological, mineralogical and geochemical studies, and mass balance calculations, it is inferred that the highly weathered host rock was originally carbonatitic in composition, initially enriched in Th and REEs compared to the surrounding silicate rocks. The intrusion of the carbonatite caused fenitic alteration in the surrounding phonolites, consisting of early potassic alteration followed by a vein-type Th-REE mineralization with associated fluorite, carbonate, pyrite and zircon. Subsequent weathering has completely decomposed the carbonatite forming a residual supergene enrichment of Th and REEs. Initial weathering of the carbonatite has created a chemical environment that might have been conductive to carbonate and phosphate complexing of the REEs in groundwaters. This may have appreciably restricted the dissolution of primary REE phases. Strongly oxidic weathering has resulted in a fractionation between Ce and the other light rare earth elements (LREEs). Ce3+ is oxidized to Ce4+ and retained together with Th by secondary mineral formation (cerianite, thorianite), and by adsorption on poorly crystalline iron- and aluminium-hydroxides. In contrast, the trivalent LREEs are retained to a lesser degree and are thus more available for secondary mineral formation (Nd-lanthanite) and adsorption at greater depths down the weathering column. Seasonally controlled fluctuations of recharge waters into the weathering column may help to explain the observed repetition of Th-Ce enriched zones underlain by trivalent LREE enriched zones.

Cosmic-ray exposure ages of pallasites

We analyzed cosmogenic nuclides in metal and/or silicate (primarily olivine) separated from the main-group pallasites Admire, Ahumada, Albin, Brahin, Brenham, Esquel, Finmarken, Glorieta Mountain, Huckitta, Imilac, Krasnojarsk, Marjalahti, Molong, Seymchan, South Bend, Springwater, and Thiel Mountains and from Eagle Station. The metal separates contained an olivine fraction which although small, <1 wt% in most cases, nonetheless contributes significantly to the budgets of some nuclides (e.g., up to 35% for Ne-21 and Al-26). A correction for olivine is therefore essential and was made using model calculations and/or empirical relations for the production rates of cosmogenic nuclides in iron meteoroids and/or measured elemental concentrations. Cosmic-ray exposure (CRE) ages for the metal phases of the main-group pallasites range from 7 to 180 Ma, but many of the ages cluster around a central peak near 100 Ma. These CRE ages suggest that the parent body of the main-group pallasites underwent a major break-up that produced most of the meteorites analyzed. The CRE age distribution for the pallasites overlaps only a small fraction of the distribution for the IIIAB iron meteorites. Most pallasites and IIIAB irons originated in different collisions, probably on different parent bodies; a few IIIABs and pallasites may have come out of the same collision but a firm conclusion requires further study. CRE ages calculated from noble gas and radionuclide data of the metal fraction are higher on average than the Ne-21 exposure ages obtained for the olivine samples. As the metal and olivine fractions were taken in most cases from different specimens, the depth-dependency of the production rate ratio Be-10/Ne-21 in metal, not accounted for in our calculations, may explain the difference.

Correlated cosmogenic W and Os isotopic variations in Carbo and implications for Hf-W chronology

An obstacle for establishing the chronology of iron meteorite formation using 182Hf-182W systematics (t1/2 = 8.9 Myr) is to find proper neutron fluence monitors to correct for cosmic ray modification of W isotopic composition. Recent studies showed that siderophile elements such as Pt and Os could serve such a purpose. To test and calibrate these neutron dosimeters, the isotopic compositions of W and Os were measured in a slab of the IID iron meteorite Carbo. This slab has a well-characterized noble gas depth profile reflecting different degrees of shielding to cosmic rays. The results show that W and Os isotopic ratios correlate with distance from the pre-atmospheric center. Negative correlations, barely resolved within error, were found between epsilo190Os-epsilo189Os and epsilo186Os-epsilo189Os with slopes of -0.64 ± 0.45 and -1.8(+1.9/-2.1), respectively. These Os isotope correlations broadly agree with model predictions for capture of secondary neutrons produced by cosmic ray irradiation and results reported previously for other groups of iron meteorites. Correlations were also found between epsilo182W-epsilo189Os (slope = 1.02 ± 0.37) and epsilo182W-epsilo190Os (slope = -1.38 ± 0.58). Intercepts of these two correlations yield pre-exposure epsilo182W values of -3.32 ± 0.51 and -3.62 ± 0.23, respectively (weighted average epsilo182W = -3.57 ± 0.21). This value relies on a large extrapolation leading to a large uncertainty but gives a metal-silicate segregation age of -0.5 ± 2.4 Myr after formation of the solar system. Combining the iron meteorite measurements with simulations of cosmogenic effects in iron meteorites, equations are presented to calculate and correct for cosmogenic effects on 182W using Os isotopes.

The unstable CO2 feedback cycle on ocean planets

Ocean planets are volatile-rich planets, not present in our Solar system, which are thought to be dominated by deep, global oceans. This results in the formation of high-pressure water ice, separating the planetary crust from the liquid ocean and, thus, also from the atmosphere. Therefore, instead of a carbonate-silicate cycle like on the Earth, the atmospheric carbon dioxide concentration is governed by the capability of the ocean to dissolve carbon dioxide (CO2). In our study, we focus on the CO2 cycle between the atmosphere and the ocean which determines the atmospheric CO2 content. The atmospheric amount of CO2 is a fundamental quantity for assessing the potential habitability of the planet's surface because of its strong greenhouse effect, which determines the planetary surface temperature to a large degree. In contrast to the stabilizing carbonate-silicate cycle regulating the long-term CO2 inventory of the Earth atmosphere, we find that the CO2 cycle feedback on ocean planets is negative and has strong destabilizing effects on the planetary climate. By using a chemistry model for oceanic CO2 dissolution and an atmospheric model for exoplanets, we show that the CO2 feedback cycle can severely limit the extension of the habitable zone for ocean planets.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.