The Ca-Eskola component in eclogitic clinopyroxene as a function of pressure, temperature and bulk composition: an experimental study to 15 GPa with possible implications for the formation of oriented SiO2-inclusions in omphacite

Experiments have been conducted in the P-T range 2.5–15 GPa and 850–1,500°C using bulk compositions in the systems SiO₂–TiO₂–Al₂O₃–Fe₂O₃–FeO–MnO–MgO–CaO–Na₂O–K₂O–P₂O₅ and SiO₂–TiO₂–Al₂O₃–MgO–CaO–Na₂O to investigate the Ca-Eskola (CaEs Ca_0.5□_0.5AlSi₂O₆) content of clinopyroxene in eclogitic assemblages containing garnet + clinopyroxene + SiO₂ ± TiO₂ ± kyanite as a function of P, T, and bulk composition. The results show that CaEs_ss in clinopyroxene increases with increasing T and is strongly bulk composition dependent whereby high CaEs-contents are favoured by bulk compositions with high normative anorthite and low diopside contents. In this study, a maximum of 18 mol% CaEs_ss was found at 6 GPa and 1,350°C in a kyanite-eclogite assemblage garnet + clinopyroxene + kyanite + rutile + coesite. By comparison, no significant increase in CaEs_ss with increasing P could be observed. If the formation of oriented SiO₂-rods frequently observed in eclogititc clinopyroxenes is due to the retrogressive breakdown of a CaEs-component then these textures are a cooling rather than a decompression phenomenon and are most likely to be found in kyanite-bearing eclogites cooled from temperatures ≥750°C. The presence of clinopyroxene with approx. 4 mol% CaEs_ss in an experiment conducted at 2.5 GPa/850°C confirms earlier suggestions based on field data that vacancy-rich clinopyroxenes are not necessarily restricted to ultrahigh pressure metamorphic conditions. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

Atomistic simulation of the mechanisms of noble gas incorporation in minerals

Atomistic simulations have been carried out to investigate the mechanisms of noble gas incorporation in minerals using both the traditional two-region approach and the “supercell” method. The traditional two-region approach has been used to calculate defect energies for Ne, Ar, Kr and Xe incorporation in MgO, CaO, diopside and forsterite in the static limit and at one atmosphere pressure. The possibilities of noble gas incorporation via both substitution and interstitial mechanisms are studied. The favored mechanism varies from mineral to mineral and from noble gas to noble gas. In all minerals studied, the variation of the solution energies of noble gas substitution with atomic radius appears approximately parabolic, analogous to those for 1⁺, 2⁺, 3⁺ and 4⁺ trace element incorporation on crystal lattice sites. Noble gas solution energies thus also fall on a curve, similar to those previously observed for cations with different charges, but with much lower curvature.The “supercell” method has been used to investigate the pressure dependence of noble gas incorporation in the same systems. Results indicate a large variation of the solubility of the larger noble gases, Kr and Xe with pressure. In addition, explicit simulation of incorporation at the (0 0 1) surface of MgO shows that the solubility of the heavier noble gases may be considerably enhanced by the presence of interfaces.     

Surface evolution of dissolving minerals investigated with a kinetic Ising model

In natural weathering systems, both the chemistry and the topography of mineral surfaces change as rocks and minerals equilibrate to surface conditions. Most geochemical research has focused on changes in solution chemistry over time; however, temporal changes in surface topography may also yield information about rates and mechanisms of dissolution. We use stochastic dissolution simulations of a regular 2-D lattice with reaction mechanisms defined in terms of nearest neighbor interactions to elucidate how the surface area and reactivity of a crystal evolve during dissolution. Despite the simplicity of the model, it reproduces key features observed or inferred for mineral dissolution. Our model results indicate that: (i) dissolving surfaces reach a steady-state conformation after sufficient dissolution time, (ii) linear defects cause surface area and dissolution rate to vary in concert with one another, (iii) sigmoidal and non-sigmoidal rate vs. free-energy of reaction (ΔG_rxn) behavior can be rationalized in terms of the multiple steps occurring during dissolution, and (iv) surface roughness as a function of ΔG_rxn is highly sensitive to the reaction mechanism. When simulated times to reach steady-state are compared to published time series rate data using suitable scaling, good agreement is found for silicate minerals while the model significantly over-predicts the duration of the transient for Fe and Al oxides. The implication of our simple model is that many aspects of mineral dissolution behavior, including approach to steady-state, sigmoidal vs. non-sigmoidal rate vs. ΔG_rxn behavior, and development of rougher surfaces in conditions further from equilibrium can be explained by nearest neighbor interactions and simple Kossel-type models where reactivity of a surface is defined in terms of perfect surface, step, and kink sites.     

Comparison of arsenic concentrations in simultaneously-collected groundwater and aquifer particles from Bangladesh,India, Vietnam,and Nepal

One of the reasons the processes resulting in As release to groundwater in southern Asia remain poorly understood is the high degree of spatial variability of physical and chemical properties in shallow aquifers. In an attempt to overcome this difficulty, a simple device that collects groundwater and sediment as a slurry from precisely the same interval was developed in Bangladesh. Recently published results from Bangladesh and India relying on the needle-sampler are augmented here with new data from 37 intervals of grey aquifer material of likely Holocene age in Vietnam and Nepal. A total of 145 samples of filtered groundwater ranging in depth from 3 to 36 m that were analyzed for As (1–1000 μg/L), Fe (0.01–40 mg/L), Mn (0.2–4 mg/L) and S (0.04–14 mg/L) are compared. The P-extractable (0.01–36 mg/kg) and HCl-extractable As (0.04–36 mg/kg) content of the particulate phase was determined in the same suite of samples, in addition to Fe(II)/Fe ratios (0.2–1.0) in the acid-leachable fraction of the particulate phase. Needle-sampler data from Bangladesh indicated a relationship between dissolved As in groundwater and P-extractable As in the particulate phase that was interpreted as an indication of adsorptive equilibrium, under sufficiently reducing conditions, across 3 orders of magnitude in concentrations according to a distribution coefficient of 4 mL/g. The more recent observations from India, Vietnam and Nepal show groundwater As concentrations that are often an order of magnitude lower at a given level of P-extractable As compared to Bangladesh, even if only the subset of particularly reducing intervals characterized by leachable Fe(II)/Fe >0.5 and dissolved Fe >0.2 mg/L are considered. Without attempting to explain why As appears to be particularly mobile in reducing aquifers of Bangladesh compared to the other regions, the consequences of increasing the distribution coefficient for As between the particulate and dissolved phase to 40 mL/g for the flushing of shallow aquifers of their initial As content are explored.     

Geochronology of Pliocene volcanism in the Dzhavakheti Highland (the Lesser Caucasus). Part 1: Western part of the Dzhavakheti Highland

Isotopic-geochronological study of the Pliocene magmatic activity in western part of the Dzhavakheti Highland (northwestern region of the Lesser Caucasus) is carried out. The results obtained imply that the Pliocene magmatic activity lasted in this part of the highland approximately 2 million years from 3.75 to 1.75–1.55 Ma. As is established, the studied volcanic rocks correspond in composition mostly to K-Na subalkaline and more abundant normal basalts. Time constraints of main phases in development of basic volcanism within the study region are figured out. We assume that individual pulses of silicic to moderately silicic volcanism presumably took place in the Dzhavakheti Highland about 3.2 and 2.5 Ma ago.     

National flood discharge mapping in Austria

This article presents the approach and the results of a study in which 30, 100 and 200 year return period flood discharges were estimated for 26,000 km of Austrian streams. Three guiding principles were adopted: combination of automatic methods and manual assessments by hydrologists to allow speedy processing and account for the local hydrological situation; combination of various sources of information including flood peak samples, rainfall data, runoff coefficients and historical flood data; and involvement of the Hydrographic Services to increase the accuracy and enhance the acceptance of results. The flood discharges for ungauged catchments were estimated by the Top-kriging approach with manual adjustment to the local flood characteristics. The adopted combination approach proved to be very efficient both in terms of the project time required and in terms of the accuracy and acceptability of the estimated flood discharges of given return periods.     

The role of the volatile phase for REE and Y fractionation in low-silica carbonate magmas: implications from natural carbonatites, Namibia

Summary Two bimodal carbonatite complexes in Namibia of Cretaceous age are explored as to the presence and composition of a coexisting carbonatitic fluid. The Kalkfeld and Ondurakorume complexes contain both Ca- and Mg/Fe-carbonatites, composed of calcite alone or calcite with ferroan dolomite, fluorapatite and strontianite. The major element evolution in the bulk rocks from s?vites to beforsites is due to crystallization of calcite and fluorapatite. All carbonatites show a negative Y anomaly in normalised REE plots. Fractionation is accompanied by successively lower HREE contents between Tb and Yb, expressed by the ratios Nd/Ho and Ho/Lu. The evolution of this downward-facing hump goes along with decreasing Y contents in bulk rocks and minerals. All this requires an additional phase coexisting with the carbonate liquid during fractionation. Comparison between the bulk rocks and the expelled fluid shows that the latter had preferentially accumulated the HREE and Y. Further evidence for this process are hydrothermal, HREE, Y-rich fluorites in other carbonatite complexes which reflect the composition of the expelled fluid. The high strength of fluoride complexes suggests that fluoride complexing in the carbonatitic fluid is the process responsible for extracting HREE and Y from the carbonatite magma, leaving fractionated carbonatite rocks depleted in these elements. The geochemical evolution of carbonatite magmas along fractionation has therefore to be considered in a melt-mineral-fluid system. Correspondence: B. Bühn, Instituto de Geociências, Universidade de Brasilia, Campus Universitário Darcy Ribeiro, Brasilia 70910-900, Brazil