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.

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Depths of Partial Crystallization of H2O-bearing MORB: Phase Equilibria Simulations of Basalts at the MAR near Ascension Island (7 11{degrees}S)

Phase equilibria simulations were performed on naturally quenchedbasaltic glasses to determine crystallization conditions priorto eruption of magmas at the Mid-Atlantic Ridge (MAR) east ofAscension Island (7–11°S). The results indicate thatmid-ocean ridge basalt (MORB) magmas beneath different segmentsof the MAR have crystallized over a wide range of pressures(100–900 MPa). However, each segment seems to have a specificcrystallization history. Nearly isobaric crystallization conditions(100–300 MPa) were obtained for the geochemically enrichedMORB magmas of the central segments, whereas normal (N)-MORBmagmas of the bounding segments are characterized by polybariccrystallization conditions (200–900 MPa). In addition,our results demonstrate close to anhydrous crystallization conditionsof N-MORBs, whereas geochemically enriched MORBs were successfullymodeled in the presence of 0·4–1 wt% H₂O in theparental melts. These estimates are in agreement with direct(Fourier transform IR) measurements of H₂O abundances in basalticglasses and melt inclusions for selected samples. Water contentsdetermined in the parental melts are in the range 0·04–0·09and 0·30–0·55 wt% H₂O for depleted and enrichedMORBs, respectively. Our results are in general agreement (within±200 MPa) with previous approaches used to evaluate pressureestimates in MORB. However, the determination of pre-eruptiveconditions of MORBs, including temperature and water contentin addition to pressure, requires the improvement of magma crystallizationmodels to simulate liquid lines of descent in the presence ofsmall amounts of water. KEY WORDS: MORB; Mid-Atlantic Ridge; depth of crystallization; water abundances; phase equilibria calculations; cotectic crystallization; pressure estimates; polybaric fractionation     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

陆面过程中稳定水同位素逐日变化的数学模拟——引入稳定同位素效应的CLM实例分析

将稳定同位素效应引入CLM(Community Land Model),并对巴西马瑙斯站在平衡年的稳定水同位素的逐日变化进行模拟和分析.结果表明: 降水、水汽和地表径流中δ18O存在明显的季节变化,并与相应的水量存在显著的负相关关系,但凝结物中δ18O与地面凝结量存在显著的正相关关系,蒸发水汽中δ18O与蒸发量之间无显著的相关关系.受土壤贮水削峰功能的影响,表层土壤和根区水中δ18O的季节变化全无.植被层蒸发水汽中稳定同位素的丰度与大气的干湿程度存在密切联系: 当降水量少时,大气干燥,植被层的蒸发较少,植被蒸发中δ18O较高;当降水量较大时,空气湿润,植被层的蒸发量较大,蒸发中δ18O则较低.植被蒸腾中δ18O的变化与源区水体中δ18O的变化保持一致,尤其是与根区水中的δ18O.由于地下径流直接源自根区水的补充,因此,地下径流中δ18O等于根区水中的δ18O.模拟结果还显示,降水MWL (大气水线)的梯度项和常数项均比全球平均MWL略偏小.尽管主要来自降水的贡献,但地表径流和植被层水体的MWLs与降水MWL存在较大的差异,这一方面与两类水体在蒸发过程中的稳定同位素的富集作用有关,另一方面与CLM模拟的水量有关.大气水汽线与降水的MWL的梯度值相近,说明大气水汽与降水近似处于稳定同位素平衡状态.另外,模拟的地面的凝结线与植被层的凝结线均与全球大气水线相近,且具有非常高相关程度,说明CLM的模拟是合理的.     

绿片岩-低角闪岩相变质条件下磁铁铁矿与黄铁矿间的Fe同位素分馏

对辽宁省鞍山一本溪地区经历了绿片岩一低角闪岩相变质的新太古代条带状铁建造中磁铁矿和黄铁矿矿物对的Fe同位素分析结果显示:相对于标准IRMM-014,所有样品的磁铁矿和黄铁矿均显示Fe的重同位素富集;且黄铁矿的Fe同位素比值均大于磁铁矿的Fe同位素比值(ε57Fe黄铁矿ε57Fe磁铁矿),两种矿物的Fe同位素比值之差为△57Fe黄铁矿-磁铁矿=2.23～5.13.黄铁矿富集铁的重同位素表明矿物的Fe同位素组成并不代表其原始沉积的特征,而是在区域变质作用过程中Fe同位素发生了交换的结果.由同位素平衡判别图解可知,在绿片岩一低角闪岩相变质作用中,磁铁矿-黄铁矿间的Fe同位素基本达到了平衡,且在平衡条件下黄铁矿比磁铁矿更富集Fe的重同位素,二者之间的Fe同位素平衡分馏系数α黄铁矿-磁铁矿≈1.000 4‰±0.06‰(2σ).这一研究成果是对变质作用过程中Fe同位素的地球化学行为认识的重要进展.     

草莓红绿柱石(Pezzottaite)的晶体结构测定及其意义

利用CCD单晶衍射仪对草莓红绿柱石进行晶体结构测定,测定其空间群为R3c,并讨论了结构中Li 有序地置换Be2 造成c/a值的变化特征。晶体中发现复三方偏三角面体的负晶包裹体,其宏观对称属m,它与空间群Rc相对应,但与测出晶体结构的空间群有矛盾。X射线结构测定已否定了其空间群为Rm的看法。要解决草莓红绿柱石晶体空间群是否有对称中心(Rc或Rc)的问题,还需进一步将X射线结构的数据及其有关的电性实验结合分析。     

九龙江河口区表层水铀同位素的粒级分布

利用微滤和超滤技术研究了九龙江河口区表层水中铀及其同位素组成的粒级分布和地球化学行为.结果表明,溶解态(<0.4μm)中低分子量组分(<10 000 u)占主要份额,胶体态(10 000 u~0.4μm)238U所占比例不足1%,且随盐度的增加其所占份额逐渐降低.溶解态、低分子量组分和胶体态238U的比活度与盐度之间存在良好的线性正相关关系,证实它们在九龙江河口区呈现保守行为.在颗粒态(>0.4μm)中,各粒级组分238U所占份额主要受控于相应颗粒物的浓度,在盐度小于20的区域,各粒级颗粒组分238U占颗粒态的份额有如下变化次序:10~53μm>2~10μm>0.4~2μm>大于53μm,而在盐度大于30的近外海站位,该次序发生一些变化:0.4~2μm>10~53μm>2~10μm>大于53μm,最小粒级颗粒组分238U的贡献有所增加,反映了自生铀贡献的加强.九龙江河口区表层水中溶解态(包括低分子量组分和胶体态)的234U/238U)_A.R.均大于1,显示出234U过剩的特征,而各粒级颗粒组分中的234U/238U)_A.R.则接近于平衡值(1.0).这一现象与陆地岩石风化过程中水体对铀的淋滤释出量及234U的优先浸出有关.对232Th/238U质量比的研究显示,溶解态及其所包括的低分子量组分和胶体态的232Th/238U质量比均小于1,而颗粒态及其所包括的4个粒级组分中的232Th/238U质量比均大于1,反映了向外海输送过程中铀、钍地球化学行为的差异.     

蒸发皿中水面蒸发氢氧同位素分馏的实验研究

气象要素与蒸发密切相关,通过室内外不同气象条件下的器皿水蒸发实验,获得了水面蒸发氢氧稳定同位素分馏因子与气象要素的关系。实验结果表明,随着蒸发的进行,剩余水体中逐渐富集重同位素;自由水体蒸发同位素分馏在垂线上有分层现象,表层水体同位素值比垂线平均的同位素值略富集;不同温度条件下的室内蒸发实验中,温度越高,液-气间分馏系数越小,相应于同一剩余水体体积比,剩余水体稳定同位素值则越低。室外器皿水自由蒸发实验中得出的蒸发线方程斜率较大地偏离了当地降水线,表明实验期间水体蒸发分馏作用较明显。该研究为进一步揭示水体蒸发分馏规律提供了可靠的实验依据。