Water release patterns of heated speleothem calcite and hydrogen isotope composition of fluid inclusions

Speleothem fluid inclusions are a potential paleo-precipitation proxy to reconstruct past rainwater isotopic composition (δ¹⁸O, δD). To get a better insight in the extraction of inclusion water from heated speleothem calcite, we monitored the water released from crushed and uncrushed speleothem calcite, heated to 900 °C at a rate of 300 °C/h, with a quadrupole mass spectrometer. Crushed calcite released water in three not well individualised peaks between 25 and 360 °C, 360 and 650 °C and between 650 and 800 °C while uncrushed calcite released water in two distinct temperature intervals: between 25 and 550 °C and between 550 and 900 °C.Water from two speleothems from the Han-sur-Lesse cave was recovered using three different techniques: i) the crushing and heating to 360 °C technique, ii) the decrepitation by heating to 550 °C and iii) the decomposition by heating to 900 °C technique. Measurements of the δD of water recovered by the decomposition of Han-sur-Lesse calcite heated to 900 °C did not show a 20 to 30‰ offset as found by previous authors. However a difference of 7‰ was observed between water released before and after decomposition of the calcite. Water recovery from the Han-sur-Lesse samples suggests that a simple heating technique (up to 550 °C) without crushing could both (a) recover water with δD representative of that of the drip water and (b) double the water yield as compared to the crushing and heating method.Our study warns for possible contamination of the recovered inclusion water with hydration water of lime, responsible for the recovery of water with very negative δD values.     

Chemical and boron-isotope variations in tourmalines from an S-type granite and its source rocks: the Erongo granite and tourmalinites in the Damara Belt,Namibia

Tourmaline is widespread in metapelites and pegmatites from the Neoproterozoic Damara Belt, which form the basement and potential source rocks of the Cretaceous Erongo granite. This study traces the B-isotope variations in tourmalines from the basement, from the Erongo granite and from its hydrothermal stage. Tourmalines from the basement are alkali-deficient schorl-dravites, with B-isotope ratios typical for continental crust (δ¹¹B average −8.4‰ ± 1.4, n = 11; one sample at −13‰, n = 2). Virtually all tourmaline in the Erongo granite occurs in distinctive tourmaline-quartz orbicules. This “main-stage” tourmaline is alkali-deficient schorl (20–30% X-site vacancy, Fe/(Fe + Mg) 0.8–1), with uniform B-isotope compositions (δ¹¹B −8.7‰ ± 1.5, n = 49) that are indistinguishable from the basement average, suggesting that boron was derived from anatexis of the local basement rocks with no significant shift in isotopic composition. Secondary, hydrothermal tourmaline in the granite has a bimodal B-isotope distribution with one peak at about −9‰, like the main-stage tourmaline, and a second at −2‰. We propose that the tourmaline-rich orbicules formed late in the crystallization history from an immiscible Na–B–Fe-rich hydrous melt. The massive precipitation of orbicular tourmaline nearly exhausted the melt in boron and the shift of δ¹¹B to −2‰ in secondary tourmaline can be explained by Rayleigh fractionation after about 90% B-depletion in the residual fluid. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Enhancing understanding of breakdown and collapse in the Yorkshire Dales using ground penetrating radar on cave sediments

Historical evidence shows block breakdown and collapse are actively occurring in large fault aligned caverns in the Yorkshire Dales karst. Deployment of ground penetrating radar at two such sites has provided detailed images of the sedimentary sequences below the present day cavern floor but no large blocks are imaged within the sediments. Solutional processes must be removing limestone from the sediment to allow continued cavern growth. Possible mechanisms to account for the lack of large blocks within the sediment fill are discussed.     

Assessing soil erosion in Mediterranean karst landscapes of Lebanon using remote sensing and GIS

Maps showing the potential for soil erosion at 1:100,000 scale are produced in a study area within Lebanon that can be used for evaluating erosion of Mediterranean karstic terrain with two different sets of impact factors built into an erosion model. The first set of factors is: soil erodibility, morphology, land cover/use and rainfall erosivity. The second is obtained by the first adding a fifth factor, rock infiltration. High infiltration can reflect high recharge, therefore decreasing the potential of surface runoff and hence the quantity of transported materials. Infiltration is derived as a function of lithology, lineament density, karstification and drainage density, all of which can be easily extracted from satellite imagery. The influence of these factors is assessed by a weight/rate approach sharing similarities between quantitative and qualitative methods and depending on pair-wise comparison matrix.The main outcome was the production of factorial maps and erosion susceptibility maps (scale 1:100,000). Spatial and attribute comparison of erosion maps indicates that the model that includes a measure of rock infiltration better represents erosion potential. Field investigation of rills and gullies shows 87.5% precision of the model with rock infiltration. This is 17.5% greater than the precision of the model without rock infiltration. These results indicate the necessity and importance of integrating information on infiltration of rock outcrops to assess soil erosion in Mediterranean karst landscapes.     

Spectroscopic investigation of uranyl(VI) and citrate coadsorption to Al2O3

ATR-FTIR spectroscopy is used to understand the adsorption of uranyl-citrate complexes to Al₂O₃. Spectral data indicate that uranyl-citrate complexes partially dissociate upon adsorption, allowing full or partial hydrolysis of the uranyl ion. K_ads values determined for free citrate adsorption are similar to those for citrate in uranyl-citrate complexes, indicating that the complexation of uranyl by citrate does not significantly affect the ability of citrate to bond with the surface. The isotherm data also indicate enhanced citrate adsorption to Al₂O₃ in the presence of uranyl, suggesting that uranyl may be the central link between two citrate ligands, and that uranyl is associated with the surface through a bridging citrate ligand. Finally, uranyl-citrate complexes interact with citrate adsorbed to Al₂O₃ through outer sphere interactions.     

An absolute paleotemperature record from 10 to 6 Ka inferred from fluid inclusion D/H ratios of a stalagmite from Vancouver Island, British Columbia, Canada

By using continuous helium flow during the crushing of calcite speleothem samples, we are able to recover liberated inclusion waters without isotopic fractionation. A paleotemperature record for the Jacklah Jill Cave locality, Vancouver Island, BC, was obtained from a 30-cm tall stalagmite that grew 10.3-6.3 Ka ago, using δ¹⁸O values of the crushed calcite and of the inclusion water as inferred from its δD. It is found that the locality experienced mean annual temperature variations up to 11 °C over a 4-Ka period in the early Holocene. At the beginning of the period, local temperature quickly increased from a minimum of ∼1 °C to around 10 °C, but this early climate optimum, about 3 °C warmer than today, only lasted for ∼1200 years. About 8.6 Ka ago, temperature had declined to ∼7 °C, approximately the same as the modern cave temperature. Since then, the study area has experienced only minor temperature fluctuations, but there was a brief fall to ∼4 °C at around 7 Ka ago, which might be caused by a short lived expansion of local alpine glaciers. The long-term T-dependence of δD was 1.47‰/°C, identical to the value in modern precipitation.     

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.     

Evidence for crustal degassing of CF4 and SF6 in Mojave Desert groundwaters

Dissolved tetrafluoromethane (CF₄) and sulfur hexafluoride (SF₆) concentrations were measured in groundwater samples from the Eastern Morongo Basin (EMB) and Mojave River Basin (MRB) located in the southern Mojave Desert, California. Both CF₄ and SF₆ are supersaturated with respect to equilibrium with the preindustrial atmosphere at the recharge temperatures and elevations of the Mojave Desert. These observations provide the first in situ evidence for a flux of CF₄ from the lithosphere. A gradual basin-wide enhancement in dissolved CF₄ and SF₆ concentrations with groundwater age is consistent with release of these gases during weathering of the surrounding granitic alluvium. Dissolved CF₄ and SF₆ concentrations in these groundwaters also contain a deeper crustal component associated with a lithospheric flux entering the EMB and MRB through the underlying basement. The crustal flux of CF₄, but not of SF₆, is enhanced in the vicinity of local active fault systems due to release of crustal fluids during episodic fracture events driven by local tectonic activity. When fluxes of CF₄ and SF₆ into Mojave Desert groundwaters are extrapolated to the global scale they are consistent, within large uncertainties, with the fluxes required to sustain the preindustrial atmospheric abundances of CF₄ and SF₆.