Silicon isotope fractionation between rice plants and nutrient solution and its significance to the study of the silicon cycle

The silicon isotope fractionation between rice plant and nutrient solution was studied experimentally. Rice plants were grown to maturity with the hydroponic culture in a naturally lit glasshouse. The nutrient solution was sampled for 14 times during the whole rice growth period. The rice plants were collected at various growth stages and different parts of the plants were sampled separately. The silica contents of the samples were determined by the gravimetric method and the silicon isotope compositions were measured using the SiF₄ method.In the growth process, the silicon content in the nutrient solution decreased gradually from 16 mM at starting stage to 0.1-0.2 mM at harvest and the amount of silica in single rice plant increased gradually from 0.00013 g at start to 4.329 g at harvest. Within rice plant the SiO₂ fraction in roots reduced continuously from 0.23 at the seedling stage, through 0.12 at the tiller stage, 0.05 at the jointing stage, 0.023 at the heading stage, to 0.009 at the maturity stage. Accordingly, the fraction of SiO₂ in aerial parts increased from 0.77, through 0.88, 0.95, 0.977, to 0.991 for the same stages. The silicon content in roots decreased from the jointing stage, through the heading stage, to the maturity stage, parallel to the decrease of silicon content in the nutrient solution. At the maturity stage, the silicon content increased from roots, through stem and leaves, to husks, but decreased drastically from husks to grains. These observations show that transpiration and evaporation may play an important role in silica transportation and precipitation within rice plants.It was observed that the δ³⁰Si of the nutrient solution increased gradually from −0.1‰ at start to 1.5‰ at harvest, and the δ³⁰Si of silicon absorbed by bulk rice plant increased gradually from −1.72‰ at start to −0.08‰ at harvest, reflecting the effect of the kinetic silicon isotope fractionation during silicon absorption by rice plants from nutrient solutions. The calculated silicon isotope fractionation factor between the silicon instantaneously absorbed by rice roots and the silicon in nutrient solution vary from 0.9983 at start to 0.9995 at harvest, similar to those reported for bamboo, banana and diatoms in direction and extent. In the maturity stage, the δ³⁰Si value of rice organs decreased from −1.33‰ in roots to −1.98‰ in stem, and then increased through −0.16‰ in leaves and 1.24‰ in husks, to 2.21‰ in grains. This trend is similar to those observed in the field grown rice and bamboo.These quantitative data provide us a solid base for understanding the mechanisms of silicon absorption, transportation and precipitation in rice plants and the role of rice growth in the continental Si cycle.     

植物中硅矿化作用的硅同位素示踪研究

项目首次对田地生长的水稻与竹子和室内栽培水稻中氧化硅的含量、形态、分布及硅同位素组成进行了系统研究。研究发现水稻中的氧化硅含量有由根到茎、叶、稻壳逐渐增高的趋势, 但在米粒中含量急剧降低。竹子中的氧化硅含量也显由杆到枝、叶逐渐增高的趋势。在竹子和水稻的根部, 氧化硅都集中在内皮层;而在其地上部分(竿、枝、叶、壳), 氧化硅主要出现于外皮层。在单株水稻和竹子中都发现不同器官间存在显著的、系统的硅同位素分馏。水稻的? 30Si显示有由根到茎降低, 而后向叶、壳和米逐渐增高的趋势。竹子的? 30Si也显由根到竿降低, 而后向枝、叶增高的趋势。这种硅同位素变化可能是由植株内体液中的溶解硅在竿、枝、叶、壳相继沉淀出氧化硅时, 产生瑞利过程的硅同位素分馏的结果。研究得出竹子和水稻中溶解硅与沉淀硅间的硅同位素分馏系数分别为0.9981和0.9996。研究发现水稻根和竹根从土壤溶液中吸取硅时, 也存在硅同位素动力分馏。竹子与水稻吸收硅与土壤可溶硅之间的硅同位素分馏系数分别为0.9988和0.9989。研究得出：1)水稻与竹子由外界吸收的含硅化合物主要为正硅酸;2)被动吸收是其吸收硅的重要形式;3)蒸发作用是硅在这些植物中迁移和沉淀的主要机制。研究结果为理解植物中硅吸收、搬运和沉淀硅的方式与机制和探讨植物在硅、碳生物地球化学循环方面的作用提供了可靠的证据     

Silicon isotope compositions of dissolved silicon and suspended matter in the Yangtze River, China

Silicon isotope compositions of main channel samples of the Yangtze River were systematically investigated along with their chemical compositions. The concentration of suspended matter in the Yangtze River tends to decrease from the upper reaches to the lower reaches, corresponding to settling of the sediments in the lakes and reservoirs due to reduction of the velocity of water flow. The silica contents of suspended matter vary from 52.1% to 56.9% and their δ³⁰Si values vary from 0 to −0.7‰, both similar to those of shales. From the upper to lower reaches, the silica contents of suspended matter tend to increase, whilst their δ³⁰Si values tend to decrease. Both trends reflect the increase of clay minerals and decrease of carbonates in suspended matter.The concentrations of dissolved silicon vary from 97 to 121 μmol/L and their δ³⁰Si values vary over a wide range from 0.7 to 3.4‰. From the upper to lower reaches, dissolved silica concentrations tend to decrease and their δ³⁰Si values tend to increase. These trends mainly reflect the change of chemical and isotopic characteristics of the tributaries from the upper to lower reaches. The major factors responsible for these changes may be the high meteoric precipitation and significant silicon absorption by grass (in wetlands) and rice (in paddy fields) in drainage areas of the middle and lower reaches.There is no correlation between δ³⁰Si of dissolved silicon and that of suspended matter. The Δ³⁰Si_Diss-SPM values vary over a wide range of 1.0-3.7‰, indicating that (1) they are out of isotopic equilibrium, (2) dissolved silicon and the associated suspended matter do not belong to one physico-chemical system, and (3) isotopic exchange rate between them is very slow.The δ³⁰Si value of dissolved silicon output from the Yangtze River to the East Sea is estimated to be 3.0‰, much higher than the values reported for the Amazon and Congo rivers. This increases the δ³⁰Si range of dissolved silicon in the world’s rivers from 0.4-1.2%; to 0.4-3.4%.     

Silicon isotope and trace element constraints on the origin of ∼3.5 Ga cherts: Implications for Early Archaean marine environments

Silicon (Si) isotope variability in Precambrian chert deposits is significant, but proposed explanations for the observed heterogeneity are incomplete in terms of silica provenance and fractionation mechanisms involved. To address these issues we investigated Si isotope systematics, in conjunction with geochemical and mineralogical data, in three well-characterised and approximately contemporaneous, ∼3.5 Ga chert units from the Pilbara greenstone terrane (Western Australia).We show that Si isotope variation in these cherts is large (−2.4‰ to +1.3‰) and was induced by near-surface processes that were controlled by ambient conditions. Cherts that formed by chemical precipitation of silica show the largest spread in δ³⁰Si (−2.4‰ to +0.6‰) and are characterised by positive Eu, La and Y anomalies and overall depletions in lithophile trace elements. Silicon isotope systematics in these orthochemical deposits are explained by (1) mixing between hydrothermal fluids and seawater, and/or (2) fractionation of hydrothermal fluids by subsurface losses of silica due to conductive cooling. Rayleigh-type fractionation of hydrothermal fluids was largely controlled by temperature differences between these fluids and seawater. Lamina-scale Si isotope heterogeneity within individual chemical chert samples up to 2.2‰ is considered to reflect the dynamic nature of hydrothermal activity. Silicified volcanogenic sediments lack diagnostic REE+Y anomalies, are enriched in lithophile elements, and exhibit a much more restricted range of positive δ³⁰Si (+0.1‰ to +1.1‰), which points to seawater as the dominant source of silica.The proposed model for Si isotope variability in the Early Archaean implies that chemical cherts with the most negative δ³⁰Si formed from pristine hydrothermal fluids, whereas silicified or chemical sediments with positive δ³⁰Si are closest to pure seawater deposits. Taking the most positive value found in this study (+1.3‰), and assuming that the Si isotope composition of seawater is governed by input of fractionated hydrothermal fluids, we infer that the temperature of ∼3.5 Ga seawater was below ∼55 °C.     

Mechanisms controlling silicon isotope distribution in the Eastern Equatorial Pacific

The distribution of silicon isotopes along a meridional transect at 140°W longitude in the Eastern Equatorial Pacific was used to test the hypothesis that δ³⁰Si of silicic acid in surface waters should correlate with net silica production rates (gross silica production minus silica dissolution) rather than rates of gross silica production due to the opposing Si isotope fractionations associated with silica production and silica dissolution. Variations in δ³⁰Si appeared significantly correlated with net silica production rates in equatorial surface waters and not with gross production rates. Around the Equator, values of δ³⁰Si as low as deep water values occurred in the upper mesopelagic in a zone of net silica dissolution and high detrital biogenic silica content, where the release of low δ³⁰Si silicic acid from opal dissolution would be expected to decrease δ³⁰Si. The δ³⁰Si of the deep water at 140°W appears constant for depths >2000 m and is similar to the deep water at 110°W. This study brings to light the importance of considering Si fractionation during diatom silica dissolution, the biological fractionation during silica production and physical factors such as currents and mixing with adjacent water masses when interpreting silicon isotope distributions.     

Equilibrium and kinetic Si isotope fractionation factors and their implications for Si isotope distributions in the Earth’s surface environments

Several important equilibrium Si isotope fractionation factors among minerals, organic molecules and the H₄SiO₄ solution are complemented to facilitate the explanation of the distributions of Si isotopes in Earth’s surface environments. The results reveal that, in comparison to aqueous H₄SiO₄, heavy Si isotopes will be significantly enriched in secondary silicate minerals. On the contrary, quadra-coordinated organosilicon complexes are enriched in light silicon isotope relative to the solution. The extent of ²⁸Si-enrichment in hyper-coordinated organosilicon complexes was found to be the largest. In addition, the large kinetic isotope effect associated with the polymerization of monosilicic acid and dimer was calculated, and the results support the previous statement that highly ²⁸Si-enrichment in the formation of amorphous quartz precursor contributes to the discrepancy between theoretical calculations and field observations. With the equilibrium Si isotope fractionation factors provided here, Si isotope distributions in many of Earth’s surface systems can be explained. For example, the change of bulk soil δ³⁰Si can be predicted as a concave pattern with respect to the weathering degree, with the minimum value where allophane completely dissolves and the total amount of sesqui-oxides and poorly crystalline minerals reaches their maximum. When, under equilibrium conditions, the well-crystallized clays start to precipitate from the pore solutions, the bulk soil δ³⁰Si will increase again and reach a constant value. Similarly, the precipitation of crystalline smectite and the dissolution of poorly crystalline kaolinite may explain the δ³⁰Si variations in the ground water profile. The equilibrium Si isotope fractionations among the quadra-coordinated organosilicon complexes and the H₄SiO₄ solution may also shed light on the Si isotope distributions in the Si-accumulating plants.     

Impact of soil weathering degree on silicon isotopic fractionation during adsorption onto iron oxides in basaltic ash soils, Cameroon

The sequestration of silicon in soil clay-sized iron oxides may affect the terrestrial cycle of Si. Iron oxides indeed specifically adsorb aqueous monosilicic acid (H₄SiO₄⁰), thereby influencing Si concentration in soil solution. Here we study the impact of H₄SiO₄⁰ adsorption on the fractionation of Si isotopes in basaltic ash soils differing in weathering degree (from two weathering sequences, Cameroon), hence in clay and Fe-oxide contents, and evaluate the potential isotopic impact on dissolved Si in surrounding Cameroon rivers. Adsorption was measured in batch experiment series designed as function of time (0-72 h) and initial concentration (ic) of Si in solution (0.61-1.18 mM) at 20 °C, constant pH (5.5) and ionic strength (1 mM). After various soil-solution contact times, the δ³⁰Si vs. NBS28 compositions were determined in selected solutions by MC-ICP-MS (Nu Plasma) in medium resolution, operating in dry plasma with Mg doping with an average precision of ±0.15‰ (±2σ_SEM). The quantitative adsorption of H₄SiO₄⁰ by soil Fe-oxides left a solution depleted in light Si isotopes, which confirms previous study on synthetic Fe-oxides. Measured against its initial composition (δ³⁰Si = +0.02 ± 0.07‰ (±2σ_SD)), the solutions were systematically enriched in ³⁰Si reaching maximum δ³⁰Si values ranging between +0.16‰ and +0.95‰ after 72 h contact time. The enrichment of the solution in heavy isotopes increased with increasing values of three parameters: soil weathering degree, iron oxide content, and proportion of short-range ordered Fe-oxide. The Si-isotopic signature of the solution was partly influenced by Si release, possibly through mineral dissolution and Si desorption from oxide surfaces, depending on soil type, highlighting the complex pattern of natural soils. Surrounding Cameroon rivers displayed a mean Si-isotopic signature of +1.19‰. Our data imply that in natural environments, H₄SiO₄⁰ adsorption by soil clay-sized Fe-oxides at least partly impacts the Si-isotopic signature of the soil solution exported to water streams.     

Sources and biological fractionation of Silicon isotopes in the Eastern Equatorial Pacific

Silicon isotopes in dissolved silicic acid were measured in the upper four kilometers between 4°N and 3°S latitude at 110°W longitude in the eastern Equatorial Pacific. Silicon isotopes became progressively heavier with silicic acid depletion of surface water as expected from biological fractionation. The value of ε estimated by applying a steady-state isotope fractionation model to data from all stations between 4°N and 3°S was −0.77 ± 0.12‰ (std. err.). When the analysis was restricted to those stations whose temperature and salinity profiles indicated that they were directly influenced by upwelling of the Equatorial Undercurrent (EUC), the resulting value of ε was −1.08 ± 0.27‰ (std. err.) similar to the value established in culture studies (−1.1‰). When the non steady state Rayleigh model was applied to the same restricted data set the resulting value of ε was significantly more positive, −0.61 ± 0.16‰ (std. err.). To the extent that the equatorial system approximates a steady state these results support a value of −1.1‰ for the fractionation factor for isotopes of Si in the sea. Without the assumption of steady state the value of ε can only be constrained to be between −0.6 and −1.1‰. Silicic acid in Equatorial Pacific Deep Water below 2000 m had a near constant δ³⁰Si of +1.32 ± 0.05‰. That value is significantly more positive than obtained for North Pacific Deep Water at similar depths at stations to the northwest of our study area (0.9-1.0‰) and it is slightly less positive than new measures of the δ³⁰Si of silicic acid from the silicic acid plume centered over the Cascadia basin in the Northeast Pacific (Si(OH)₄ > 180  μM, δ³⁰Si = +1.46 ± 0.12‰ (SD, n = 4). We show that the data from the equator and Cascadia basin fit a general trend of increasing δ³⁰Si(OH)₄ with increasing silicic acid concentration in the deep sea, but that the isotope values from the Northeast Pacific are anomalously light. The observed level of variation in the silicon isotope composition of deep waters from this single ocean basin is considerably larger than that predicted by current models based on fractionation during opal formation with no isotope effect during dissolution. Confirmation of such high variability in deep water δ³⁰Si(OH)₄ within individual ocean basins will require reassessment of the mechanisms controlling the distribution of isotopes of silicon in the sea.     

SIMS analyses of silicon and oxygen isotope ratios for quartz from Archean and Paleoproterozoic banded iron formations

Banded iron formations (BIFs) are chemical marine sediments dominantly composed of alternating iron-rich (oxide, carbonate, sulfide) and silicon-rich (chert, jasper) layers. Isotope ratios of iron, carbon, and sulfur in BIF iron-bearing minerals are biosignatures that reflect microbial cycling for these elements in BIFs. While much attention has focused on iron, banded iron formations are equally banded silica formations. Thus, silicon isotope ratios for quartz can provide insight on the sources and cycling of silicon in BIFs. BIFs are banded by definition, and microlaminae, or sub-mm banding, are characteristic of many BIFs. In situ microanalysis including secondary ion mass spectrometry is well-suited for analyzing such small features. In this study we used a CAMECA IMS-1280 ion microprobe to obtain highly accurate (±0.3‰) and spatially resolved (∼10 μm spot size) analyses of silicon and oxygen isotope ratios for quartz from several well known BIFs: Isua, southwest Greenland (∼3.8 Ga); Hamersley Group, Western Australia (∼2.5 Ga); Transvaal Group, South Africa (∼2.5 Ga); and Biwabik Iron Formation, Minnesota, USA (∼1.9 Ga). Values of δ¹⁸O range from +7.9‰ to +27.5‰ and include the highest reported δ¹⁸O values for BIF quartz. Values of δ³⁰Si have a range of ∼5‰ from −3.7‰ to +1.2‰ and extend to the lowest δ³⁰Si values for Precambrian cherts. Isua BIF samples are homogeneous in δ¹⁸O to ±0.3‰ at mm- to cm-scale, but are heterogeneous in δ³⁰Si up to 3‰, similar to the range in δ³⁰Si found in BIFs that have not experienced high temperature metamorphism (up to 300 °C). Values of δ³⁰Si for quartz are homogeneous to ±0.3‰ in individual sub-mm laminae, but vary by up to 3‰ between multiple laminae over mm-to-cm of vertical banding. The scale of exchange for Si in quartz in BIFs is thus limited to the size of microlaminae, or less than ∼1 mm. We interpret differences in δ³⁰Si between microlaminae as preserved from primary deposition. Silicon in BIF quartz is mostly of marine hydrothermal origin (δ³⁰Si < −0.5‰) but silicon from continental weathering (δ³⁰Si ∼ 1‰) was an important source as early as 3.8 Ga.     

Experimental study of silicon isotope dynamic fractionation and its application in geology

Silicon shows no variation in its chemical valence in nature and exists mainly in the form of silicon-oxygen tetrahedra, so very small silicon isotope thermodynamic fractionation occurs and the resultant silicon isotope variation is limited. Dynamic fractionation of Si isotopes during precipitation of SiO₂ from a solution is a main factor leading to substantial variations in silicon isotopes in nature. In this experimental study, we determined the dynamic fractionation factorα for silicon isotopes during precipitation of SiO₂ from the solution. And in combination of α, a theoretical explanation is presented of the considerably low δ³⁰Si values of black smokers on modern seafloor, Archean banded magnetite-quartzite and clay minerals of weathering origin, and of clearly high δ³⁰Si values of siliceous rocks in shallow-sea carbonate platforms. This paper won the Paper of Excellence in the Second National Young Scientist Symposium on Geochemistry of Minerals and Rocks.     

The silicon isotopic composition of surface waters in the Atlantic and Indian sectors of the Southern Ocean

We report here the silicon isotopic composition (δ³⁰Si) of dissolved silicon (DSi) from 42 surface water samples from the Drake Passage, the Weddell Gyre, other areas south of the Southern Boundary of the Antarctic Circumpolar Current (ACC), and the ACC near the Kerguelen Plateau, taken between the beginning of February and the end of March 2007. From the beginning to end of the cruise (ANTXXIII/9), DSi diminished in the Antarctic by 50 μmol L⁻¹ while concentrations of nitrate + nitrite and phosphate showed no net decline, indicating that the high seasonal Si/N removal ratios well known for the Southern Ocean may be more related to the strength of the silicate pump in the Southern Ocean than to the instantaneous Si/N uptake ratio of diatoms. The δ³⁰Si of DSi in samples containing more than 20 μM DSi were strongly negatively correlated to DSi concentrations, supporting the use of δ³⁰Si as a proxy for DSi removal. The “open system” fractionation observed, ε = −1.2 ± 0.11‰, agrees well with results from previous work in other areas, and the estimate of the initial δ³⁰Si of DSi of +1.4‰ is not far off observations of the δ³⁰Si of DSi in Winter Water (WW) in this area. Results were used to model DSi draw down in the past from the δ³⁰Si of sediment cores, although isotopic fractionation during silica dissolution appeared to influence the δ³⁰Si of some surface water samples, inviting further study of this phenomenon.     

Natural variations of δSi ratios during progressive basalt weathering, Hawaiian Islands

Silicon stable isotopes can be used to trace the biogeochemical pathways of Si as it moves from its continental sources to its sink in ocean sediments. Along the way, Si is incorporated into clay minerals, taken up by plants where it forms plant opal, and leached into rivers, the major land-to-ocean conduit. Compared to igneous rocks, the waters that drain continents are enriched in heavy Si isotopes, but the mechanisms that control fractionation have not been elucidated. We studied Si isotope fractionation along a 4 million yr basaltic soil chronosequence on the Hawaiian Islands. Using the natural context of these samples in combination with laboratory experiments, we demonstrate that the isotopic composition of dissolved Si in weathering systems is determined by the combined effects of rock disintegration, clay mineral neosynthesis, and Si biocycling. Weathering preferentially releases ²⁸Si into solution, whereas secondary mineral formation preferentially removes ²⁸Si from solution. In humid environments, leached soils have lost large amounts of this soluble Si, thus creating a net loss of ³⁰Si from the entire soil system. As soils develop and greater fractions of Si reside in neoformed clay minerals, δ³⁰Si_{bulk soil} values change progressively toward more negative values; basalt δ³⁰Si values are about −0.5‰, but older soils have δ³⁰Si values up to −2.5‰. The difference between the solid and solution δ³⁰Si values remains more or less constant with progressive weathering, and therefore, soil water from older soils has a more negative δ³⁰Si composition. In the upper horizons of the Hawaiian soils, this weathering-driven δ³⁰Si shift is modified by the addition of unweathered primary minerals via dust, carrying δ³⁰Si values of about −0.5‰, and by biocycling of Si via plants, producing negative δ³⁰Si values in phytoliths and positive δ³⁰Si values in soil solutions derived from upper horizons. Due to the high concentrations of dissolved Si in these near-surface layers, rivers have more positive δ³⁰Si values than predicted based on the weathering status of the lower horizons. When combined with published δ³⁰Si values from large rivers worldwide, we find that the results from Hawaii point to weathering control of Si isotopes delivered to the oceans, and thus, to an important continent-ocean linkage that warrants further investigation.     

First-principles calculations of equilibrium fractionation of O and Si isotopes in quartz,albite, anorthite,and zircon

In this study, we used first-principles calculations based on density functional theory to investigate silicon and oxygen isotope fractionation factors among the most abundant major silicate minerals in granites, i.e., quartz and plagioclase (including albite and anorthite), and an important accessory mineral zircon. Combined with previous results of minerals commonly occurring in the crust and upper mantle (orthoenstatite, clinoenstatite, garnet, and olivine), our study reveals that the Si isotope fractionations in minerals are strongly correlated with SiO₄ tetrahedron volume (or average Si–O bond length). The ³⁰Si enrichment order follows the sequence of quartz > albite > anorthite > olivine ≈ zircon > enstatite > diopside, and the ¹⁸O enrichment follows the order of quartz > albite > anorthite > enstatite > zircon > olivine. Our calculation predicts that measurable fractionation of Si isotopes can occur among crustal silicate minerals during high-temperature geochemical processes. This work also allows us to evaluate Si isotope fractionation between minerals and silicate melts with variable compositions. Trajectory for δ³⁰Si variation during fractional crystallization of silicate minerals was simulated with our calculated Si isotope fractionation factors between minerals and melts, suggesting the important roles of fractional crystallization to cause Si isotopic variations during magmatic differentiation. Our study also predicts that δ³⁰Si data of ferroan anorthosites of the Moon can be explained by crystallization and aggregation of anorthite during lunar magma ocean processes. Finally, O and Si isotope fractionation factors between zircon and melts were estimated based on our calculation, which can be used to quantitatively account for O and Si isotope composition of zircons crystallized during magma differentiation.     

Experimental determination of magnesium isotope fractionation during higher plant growth

Two higher plant species (rye grass and clover) were cultivated under laboratory conditions on two substrates (solution, phlogopite) in order to constrain the corresponding Mg isotope fractionations during plant growth and Mg uptake. We show that bulk plants are systematically enriched in heavy isotopes relative to their nutrient source. The Δ²⁶Mg_plant-source range from 0.72‰ to 0.26‰ for rye grass and from 1.05‰ to 0.41‰ for clover. Plants grown on phlogopite display Mg isotope signatures (relative to the Mg source) ∼0.3‰ lower than hydroponic plants. For a given substrate, rye grass display lower δ²⁶Mg (by ∼0.3‰) relative to clover. Magnesium desorbed from rye grass roots display a δ²⁶Mg greater than the nutrient solution. Adsorption experiments on dead and living rye grass roots also indicate a significant enrichment in heavy isotopes of the Mg adsorbed on the root surface. Our results indicate that the key processes responsible for heavy isotope enrichment in plants are located at the root level. Both species also exhibit an enrichment in light isotopes from roots to shoots (Δ²⁶Mg_leaf-root = −0.65‰ and −0.34‰ for rye grass and clover grown on phlogopite respectively, and Δ²⁶Mg_leaf-root of −0.06‰ and −0.22‰ for the same species grown hydroponically). This heavy isotope depletion in leaves can be explained by biological processes that affect leaves and roots differently: (1) organo-Mg complex (including chlorophyll) formation, and (2) Mg transport within plant. For both species, a positive correlation between δ²⁶Mg and K/Mg was observed among the various organs. This correlation is consistent with the link between K and Mg internal cycles, as well as with formation of organo-magnesium compounds associated with enrichment in heavy isotopes. Considering our results together with the published range for δ²⁶Mg of natural plants and rivers, we estimate that a significant change in continental vegetation would induce a change of the mean river δ²⁶Mg that is comparable to analytical uncertainties.     

Micro‐scale silicon isotope heterogeneity observed in hydrothermal quartz precipitates from the >3.7 Ga Isua Greenstone Belt,SW Greenland

Pillow basalt and chert form integral lithologies comprising many Archean greenstone belt packages. To investigate details of these lithologies in the >3.7 Ga Isua Greenstone Belt, SW Greenland, we measured silicon isotope compositions of quartz crystals, by secondary ion mass spectrometry, from a quartz‐cemented, quartz‐amygdaloidal basaltic pillow breccia, recrystallized chert and chert clasts thought to represent silica precipitation under hydrothermal conditions. The recrystallized chert, chert clasts and quartz cement have overlapping δ³⁰Si values, while the δ³⁰Si values of the quartz amygdules span nearly the entire range of previously published values for quartz precipitates of any age, despite amphibolite facies metamorphism. We suggest that the heterogeneity is derived from kinetic isotope fractionation during quartz precipitation under disequilibrium conditions in a hydrothermal setting, consistent with the pillow breccia origin. On the basis of the present data, we conclude that the geological context of each sample must be carefully evaluated when interpreting δ³⁰Si values of quartz.     

Weathering, dust, and biocycling effects on soil silicon isotope ratios

Silicon isotope ratios (δ³⁰Si) of bulk mineral materials in soil integrate effects from both silicon sources and processing. Here we report δ³⁰Si values from a climate gradient of Hawaiian soils developed on 170 ka basalt and relate them to patterns of soil chemistry and mineralogy. The results demonstrate informative relationships between the mass fraction of soil Si depletion and δ³⁰Si. In upper (<1 m deep) soil horizons along the climate gradient, Si depletion correlates with decreases of residual δ³⁰Si values in low rainfall soils and increases in high rainfall soils. Strong positive correlation between soil δ³⁰Si and dust-derived quartz and mica content show that both trends are largely controlled by the abundance of these weathering-resistant minerals. The data also lend support to the idea that fractionation of Si isotopes in secondary phases is controlled by partitioning of silicon between dissolved and precipitated products during the initial weathering of primary basalt. Secondary mineral δ³⁰Si values from lower (>1 m deep) soil horizons generally correlate with the isotope fractionation predicted by a study of dissolved Si in basalt-watershed rivers and driven by preferential ²⁸Si removal from the dissolved phase during precipitation. In contrast, after correcting for the influence of dust, secondary mineral Si depletion and δ³⁰Si values in shallow (<1 m deep) soil horizons showed evidence of biocycling induced Si redistribution and substantially lower δ³⁰Si values than predicted. Low δ³⁰Si values in shallow soil horizons compared to predictions can be attributed to repeated fractionation as secondary minerals undergo additional cycles of dissolution and precipitation. Primary mineral weathering, secondary mineral weathering, dust accumulation, and biocycling are major processes in terrestrial Si cycling and these results demonstrate that each can be traced by δ³⁰Si values interpreted in conjunction with mineralogy and measures of Si depletion.     

Silicon isotopes in the inner Solar System: Implications for core formation, solar nebular processes and partial melting

We report Si isotopic data on a suite of terrestrial mantle-derived samples, meteorites and a lunar sample. Our data on co-existing mantle minerals, peridotites and basalts demonstrate lack of any resolvable high temperature fractionation during igneous processes. We show that the δ³⁰Si of the bulk silicate Earth (BSE) is identical, within analytical uncertainties, to carbonaceous and ordinary chondrites (CHUR). Based on our data the difference between δ³⁰Si_BSE and δ³⁰Si_CHUR is 0.035 ± 0.035. Whole-rock differentiated meteorites from different parent bodies (Mars, Vesta) and a lunar breccia sample also show similar δ³⁰Si suggesting broad-scale Si isotope homogeneity in the inner Solar System with an average δ²⁹Si = −0.20 ± 0.01 and δ³⁰Si = −0.39 ± 0.02 relative to the NBS28 Si isotope standard.A difference between δ³⁰Si_BSE and δ³⁰Si_CHUR of 0.035, as observed in our study, translates to less than 1.67 wt.% Si in the core considering a continuous accretion model whereas estimates using a batch model are even lower. Within uncertainties (±0.035‰) in the δ³⁰Si difference between the BSE and CHUR, a maximum of 3.84 wt.% Si could be present in the Earth’s core whereas at δ³⁰Si_BSE-δ³⁰Si_CHUR = 0, there is no requirement of Si in the Earth’s core. Such low Si in the core necessitates the presence of other light elements in the core to explain its density deficit. Our data also places constraints on the oxidation state of the Earth’s mantle during core segregation. The uncertainties in estimating the concentration of oxidized Fe in the mantle during the first 90% of accretion arise from uncertainties in the estimates of the equilibrium partition coefficient of silicon between metal and silicate at conditions relevant to core formation. For δ³⁰Si_BSE-δ³⁰Si_CHUR = 0.035 ± 0.035, the concentration of oxidized Fe in the mantle during the first 90% of accretion could be as low as ∼1%. However, at δ³⁰Si_BSE-δ³⁰Si_CHUR = 0, the Si isotope data do not require any change in the mantle concentration of oxidized Fe during accretion from the present day value of 6.26%.     

Silicon isotopic composition of dissolved silicon and suspended particulate matter in the Yellow River, China, with implications for the global silicon cycle

The silicon isotopic composition of dissolved silicon and suspended particulate matter (SPM) were systematically investigated in water samples from the mainstem of the Yellow River and 4 major tributaries. The SPM content of the Yellow River varied from 1.4 to 38,560 mg/L, averaging 3568 mg/L, and the δ³⁰Si of suspended particulate matter (δ³⁰Si_SPM) varied from 0.3‰ to −0.4‰, averaging −0.02‰. The major factors affecting the SPM content and the δ³⁰Si_SPM values in the Yellow River were inferred to be the mineralogical, chemical and isotopic characteristics of the sediments from the Loess Plateau and a combination of the climate and the flow discharge of the river.The major ions in the Yellow River water were Na⁺, Ca²⁺, Mg²⁺, HCO₃⁻, SO₄²⁻ and Cl⁻. High salt concentration was observed in samples from the middle and lower reaches, likely reflecting the effects of evaporation and irrigation because the Na⁺, Mg²⁺, SO₄²⁻ and K⁺ concentrations were correlated with the Cl⁻ concentration. The dissolved Si concentration (D_Si) increased downstream, varying from 0.016 to 0.323 mM. The δ³⁰Si of dissolved Si (δ³⁰Si_Diss) varied from 0.4‰ to 2.5‰, averaging 1.28‰. The major processes controlling the D_Si and δ³⁰Si_Diss of the Yellow River are (a) the weathering of silicate rocks, (b) the formation of phytoliths in plants, (c) the evaporation of water from and the addition of meteoric water to the river system, which only affects concentrations, (d) the adsorption and desorption of aqueous monosilicic acid on iron oxide, and (e) the dissolution of phytoliths in soils.The D_Si and δ³⁰Si_Diss values of global rivers vary spatially and temporally in response to changes in climate, chemical weathering intensity and biological activity. The moderately positive δ³⁰Si_Diss values observed in the Yellow River may be attributed to the higher rates of chemical weathering and biological activities that have been observed in this catchment in comparison with those of other previously studied catchments, excluding the Yangtze River. Human activities may also potentially influence chemical weathering and biological activities and affect the D_Si and δ³⁰Si_Diss values of the major rivers of the world. Further river studies should be performed to gain a better understanding of the global Si isotope budget.     

Silicon isotope fractionation during magmatic differentiation

The Si isotopic composition of Earth’s mantle is thought to be homogeneous (δ³⁰Si = −0.29 ± 0.08‰, 2 s.d.) and not greatly affected by partial melting and recycling. Previous analyses of evolved igneous material indicate that such rocks are isotopically heavy relative to the mantle. To understand this variation, it is necessary to investigate the degree of Si isotopic fractionation that takes place during magmatic differentiation. Here we report Si isotopic compositions of lavas from Hekla volcano, Iceland, which has formed in a region devoid of old, geochemically diverse crust. We show that Si isotopic composition varies linearly as a function of silica content, with more differentiated rocks possessing heavier isotopic compositions. Data for samples from the Afar Rift Zone, as well as various igneous USGS standards are collinear with the Hekla trend, providing evidence of a fundamental relationship between magmatic differentiation and Si isotopes. The effect of fractionation has been tested by studying cumulates from the Skaergaard Complex, which show that olivine and pyroxene are isotopically light, and plagioclase heavy, relative to the Si isotopic composition of the Earth’s mantle. Therefore, Si isotopes can be utilised to model the competing effects of mafic and felsic mineral fractionation in evolving silicate liquids and cumulates.At an average SiO₂ content of ∼60 wt.%, the predicted δ³⁰Si value of the continental crust that should result from magmatic fractionation alone is −0.23 ± 0.05‰ (2 s.e.), barely heavier than the mantle. This is, at most, a maximum estimate, as this does not take into account weathered material whose formation drives the products toward lighter δ³⁰Si values. Mass balance calculations suggest that removal of continental crust of this composition from the upper mantle will not affect the Si isotopic composition of the mantle.     

CaO-MgO-Al2o3-SiO2 liquids: chemical and isotopic effects of Mg and Si evaporation in a closed system of solar composition

A method is shown for calculating vapor pressures over a CMAS droplet in a gas of any composition. It is applied to the problem of the evolution of the chemical and Mg and Si isotopic composition of a completely molten droplet having the composition of a likely refractory inclusion precursor during its evaporation into the complementary, i.e. modified solar, gas from which it originally condensed, a more realistic model than previous calculations in which the ambient gas is pure H_2(g). Because the loss rate of Mg is greater than that of Si, the vapor pressure of Mg_(g) falls and its ambient pressure rises faster than those of SiO_(g) during isothermal evaporation, causing the flux of Mg_(g) to approach zero faster and MgO to approach its equilibrium concentration sooner than SiO₂. As time passes, δ²⁵Mg and δ²⁹Si increase in the droplet and decrease in the ambient gas. The net flux of each isotope crossing the droplet/gas interface is the difference between its outgoing and incoming flux. δ²⁵Mg and δ²⁹Si of this instantaneous gas become higher, first overtaking their values in the ambient gas, causing them to increase with time, and later overtaking their values in the droplet itself, causing them to decrease with time, ultimately reaching their equilibrium values. If the system is cooling during evaporation and if mass transfer ceases at the solidus temperature, 1500 K, final MgO and SiO₂ contents of the droplet are slightly higher in modified solar gas than in pure H_2(g), and the difference increases with decreasing cooling rate and increasing ambient pressure. During cooling under some conditions, net fluxes of evaporating species become negative, causing reversal of the evaporation process into a condensation process, an increase in the MgO and/or SiO₂ content of the droplet with time, and an increase in their final concentrations with increasing ambient pressure and/or dust/gas ratio. At cooling rates <∼3 K/h, closed-system evaporation at P^tot ∼ 10⁻³ bar in a modified solar gas, or at lower pressure in systems with enhanced dust/gas ratio, can yield the same δ²⁵Mg in a residual CMAS droplet for vastly different evaporated fractions of Mg. The δ²⁵Mg of a refractory residue may thus be insufficient to determine the extent of Mg loss from its precursor. Evaporation of Mg into an Mg-bearing ambient gas causes δ²⁶Mg and δ²⁵Mg of the residual droplet to fall below values expected from Rayleigh fractionation for the amount of ²⁴Mg evaporated, with the degree of departure increasing with increasing fraction evaporated and ambient pressure of Mg. δ²⁶Mg and δ²⁵Mg do not depart proportionately from Rayleigh fractionation curves, with δ²⁵Mg being less than expected on the basis of δ²⁶Mg by up to ∼1.2‰. Such departures from Rayleigh fractionation could be used in principle to distinguish heavily from lightly evaporated residues with the same δ²⁵Mg.