Partitioning behavior of chlorine and fluorine in the system apatite-melt-fluid. II: Felsic silicate systems at 200 MPa

Hydrothermal experiments were conducted to determine the partitioning of Cl between rhyolitic to rhyodacitic melts, apatite, and aqueous fluid(s) and the partitioning of F between apatite and these melts at ca. 200 MPa and 900-924 °C. The number of fluid phases in our experiments is unknown; they may have involved a single fluid or vapor plus saline liquid. The partitioning behavior of Cl between apatite and melt is non-Nernstian and is a complex function of melt composition and the Cl concentration of the system. Values of D_Cl^apat/melt (wt. fraction of: Cl in apatite/Cl in melt) vary from 1 to 4.5 and are largest when the Cl concentrations of the melt are at or near the Cl-saturation value of the melt. The Cl-saturation concentrations of silicate melts are lowest in evolved, silica-rich melts, so with elevated Cl concentrations in a system and with all else equal, the maximum values of D_Cl^apat/melt occur with the most felsic melt. In contrast, values of D_F^apat/melt range from 11 to 40 for these felsic melts, and many of these are an order of magnitude greater than those applying to basaltic melts at 200 MPa and 1066-1150 °C. The Cl concentration of apatite is a simple and linear function of the concentration of Cl in fluid. Values of D_Cl^fluid/apat for these experiments range from 9 to 43, and some values are an order of magnitude greater than those determined in 200-MPa experiments involving basaltic melts at 1066-1150 °C.In order to determine the concentrations and interpret the behavior of volatile components in magmas, the experimental data have been applied to the halogen concentrations of apatite grains from chemically evolved rocks of Augustine volcano, Alaska; Krakatau volcano, Indonesia; Mt. Pinatubo, Philippines; Mt. St. Helens, Washington; Mt. Mazama, Oregon; Lascar volcano, Chile; Santorini volcano, Greece, and the Bishop Tuff, California. The F concentrations of these magmas estimated from apatite-melt equilibria range from 0.06 to 0.12 wt% and are generally equivalent to the concentrations of F determined in the melt inclusions. In contrast, the Cl concentrations of the magmas estimated from apatite-melt equilibria (e.g., ca. 0.3-0.9 wt%) greatly exceed those determined in the melt inclusions from all of these volcanic systems except for the Bishop Tuff where the agreement is good. This discrepancy in estimated Cl concentrations of melt could result from several processes, including the hypothesis that the composition of apatite represents a comparatively Cl-enriched stage of magma evolution that precedes melt inclusion entrapment prior to the sequestration of Cl by coexisting magmatic aqueous and/or saline fluid(s).     

Low-temperature anisotropic diffusion of helium in zircon: Implications for zircon (U-Th)/He thermochronometry

In the last decade the zircon (U-Th)/He (ZHe) thermochronometer has been applied to a variety of geologic problems. Although bulk diffusion coefficients for He in zircon are available from laboratory step-heating experiments, little is known about the diffusion mechanism(s) and their dependence on the crystallographic structure of zircon. Here, we investigate the diffusion of He in perfectly crystalline zircon using atomistic simulation methods that provide insights into the structural pathways of He migration in zircon. Empirical force fields and quantum-mechanical calculations reveal that the energy barriers for He diffusion are strongly dependent on structure. The most favorable pathway for He diffusion is the [0 0 1] direction through the open channels parallel to the c-axis (, activation energy for tracer diffusion of a He atom along [0 0 1]). In contrast, energy barriers are higher in other directions where narrower channels for He diffusion are identified, such as [1 0 0], [1 0 1], and [1 1 0] (ΔE^∗ of 44.8, 101.7, and 421.3 kJ mol⁻¹, respectively). Molecular dynamics simulations are in agreement with these results and provide additional insight in the diffusion mechanisms along different crystallographic directions, as well as the temperature dependence. Below the closure temperature of He in zircon [T_c ∼ 180 °C, Reiners P. W., Spell T. L., Nicolescu S., and Zanetti K. A. (2004) Zircon (U-Th)/He thermochronometry: He diffusion and comparisons with Ar-40/Ar-39 dating. Geochim. Cosmochim. Acta68, 1857-1887], diffusion is anisotropic as He moves preferentially along the [0 0 1] direction, and calculated tracer diffusivities along the two most favorable directions differ by approximately five orders of magnitude (D_[001]/D_[100] ∼ 10⁵, at T = 25 °C). Above this temperature, He atoms start to hop between adjacent [0 0 1] channels, along [1 0 0] and [0 1 0] directions (perpendicular to the c-axis). The diffusion along [1 0 0] and [0 1 0] is thermally activated, such that at higher temperatures, He diffusion in zircon becomes nearly isotropic (D_[001]/D_[100] ∼ 10, at T = 580 °C). These results suggest that the anisotropic nature of He diffusion at temperatures near the closure temperature should be considered in future diffusivity experiments. Furthermore, care should be taken when making geologic interpretations (e.g., exhumation rates, timing of cooling, etc.) from this thermochronometer until the effects of anisotropic diffusion on bulk ages and closure temperature estimates are better quantified.     

The record of temperature, wind velocity and air humidity in the δD and δO of water inclusions in synthetic and Messinian halites

Deuterium and oxygen isotope fractionations between liquid and vapor water were experimentally-determined during evaporation of a NaCl solution (35 g L⁻¹) as a function of water temperature and wind velocity. In the case of a null wind velocity, slopes of δD-δ¹⁸O trajectories of residual waters hyperbolically decrease with increasing water temperatures in the range 23-47 °C. For wind velocities ranging from 0.8 to 2.2 m s⁻¹, slopes of the δD-δ¹⁸O trajectories linearly increase with increasing wind velocity at a given water temperature. These experimental results can be modeled by using Rayleigh distillation equations taking into account wind-related kinetics effects. Deuterium and oxygen isotope compositions of water inclusions trapped by the precipitated halite crystals were determined by micro-equilibration techniques.These isotopic compositions accurately reflect those of the surrounding residual waters during halite growth. Isotopic compositions of water inclusions in twenty natural halites from the Messinian Realmonte mine in Sicily suggest precipitation temperatures of that match the homogenization temperatures obtained by microthermometry (median = 34 ± 5 °C). The similarity between the measured and experimental slopes of the δD-δ¹⁸O evaporation trajectories suggests that the effect of wind was negligible during the genesis of these halite deposits. Hydrogen and oxygen isotope compositions of water inclusions from Realmonte halite also define a linear trend whose extrapolation until intersection with the Mediterranean Meteoric Water Line allows the characterization of the water source with δD and δ¹⁸O values of −70 ± 10‰ and −11.5 ± 1.5‰, respectively. These results reveal that the huge amounts of salts deposited in Sicily result from the evaporation of seawater mixed with a dominant fraction (?50%) of meteoric waters most likely deriving from alpine fluvial discharge.     

Partitioning of Nb and Ta between rutile and felsic melt and the fractionation of Nb/Ta during partial melting of hydrous metabasalt

In order to fully assess the role of rutile in fractionation of Nb/Ta during partial melting of hydrous metabasalt, we have measured rutile - felsic melt partition coefficients (D values) for Nb and Ta with tonalitic to trondhjemitic compositions at 1.5-3.5 GPa, 900-1350 °C and ∼5.0-20 wt% H₂O. D_Nb, D_Ta and D_Nb/D_Ta range from 17 ± 1 to 246 ± 13, 34 ± 2 to 232 ± 25 and 0.51 ± 0.04 to 1.06 ± 0.13, respectively. For the compositions investigated, melt composition appears to have no observable effect on the partitioning; the effect of pressure is also slight; whereas temperature and H₂O have marked effects. D_Nb, D_Ta and D_Nb/D_Ta increase with decreasing temperature and H₂O content, showing a reversal of D_Nb/D_Ta from <1.0 to >1.0. Using the data that approached equilibrium and obeyed Henry’s law, expressions describing the dependences of D_Nb, D_Ta and D_Nb/D_Ta on temperature, pressure and melt H₂O content were obtained:

(1)     

The distribution of sulfur between haplogranitic melts and aqueous fluids

The distribution of sulfur between haplogranitic melt and aqueous fluid has been measured as a function of oxygen fugacity (Co-CoO-buffer to hematite-magnetite buffer), pressure (0.5-3 kbar), and temperature (750-850 °C). Sulfur always strongly partitions into the fluid. At a given oxygen fugacity, pressure and temperature, the distribution of sulfur between melt and fluid can be described by one constant partition coefficient over a wide range of sulfur concentrations. Oxygen fugacity is the most important parameter controlling sulfur partitioning. While the fluid/melt partition coefficient of sulfur is 468 ± 32 under Co-CoO buffer conditions at 2 kbar and 850 °C, it decreases to 47 ± 4 at an oxygen fugacity 0.5-1 log unit above Ni-NiO at the same pressure and temperature. A further increase in oxygen fugacity to the hematite-magnetite buffer has virtually no effect on the partition coefficient (D^fluid/melt = 49 ± 2). The dependence of D^fluid/melt on temperature and pressure was systematically explored at an oxygen fugacity 0.5-1 log units above Ni-NiO. At 850 °C, the effect of pressure on the partition coefficient is small (D^fluid/melt = 58 ± 3 at 0.5 kbar; 94 ± 9 at 1 kbar; 47 ± 4 at 2 kbar and 68 ± 5 at 3 kbar) and temperature also has only a minor effect on partitioning.The data show the “sulfur excess” observed in many explosive volcanic eruptions can easily be explained by the presence of a small fraction of hydrous fluid in the magma chamber before the eruption. The sulfur excess can be calculated as the product of the fluid/melt partition coefficient of sulfur and the mass ratio of fluid over melt in the erupted material. For a plausible fluid/melt partition coefficient of 47 under oxidizing conditions, a 10-fold sulfur excess corresponds to a 17.6 wt.% of fluid in the erupted material. Large sulfur excesses (10-fold or higher) are only to be expected if only a small fraction of the magma residing in the magma chamber is erupted.The behavior of sulfur, which seems to be largely independent of pressure and temperature under oxidizing conditions is very different from chlorine, where the fluid/melt partition coefficient strongly increases with pressure. Variations in the SO₂/HCl ratio of volcanic gases, if they reflect primary processes in the magma chamber, therefore provide an indicator of pressure variations in a magma. In particular, major increases in the S/Cl ratio of an aqueous fluid coexisting with a felsic magma suggest a pressure reduction in the magma chamber and/or magma rising to the surface.     

Extreme oxygen isotope signature of meteoric water in magmatic zircon from metagranite in the Sulu orogen, China: Implications for Neoproterozoic rift magmatism

Unusual ¹⁸O depletion, with δ¹⁸O values as negative as −10‰ to −4‰ relative to VSMOW, was reported in zircons from ultrahigh-pressure eclogite-facies metamorphic rocks in the Dabie-Sulu orogenic belt, China. But it is critical for the negative δ¹⁸O zircons to be distinguished between magmatic and metamorphic origins, because the ¹⁸O depletion can be acquired by high-T eclogite-facies metamorphism of meteoric-hydrothermally altered low δ¹⁸O rocks. While zircon O diffusion kinetics has placed a reasonable constraint on this, zircon trace element compositions can provide a straightforward distinction between the magmatic and metamorphic origins. This paper reports our finding of unusual ¹⁸O depletion in zircon from granitic gneiss in the northeastern end of the Sulu orogen. Zircon δ¹⁸O values vary from −7.8‰ to −3.1‰ along a profile of 50 m length at Zaobuzhen. They are close to extremely low δ¹⁸O values of −9.0‰ to −5.9‰ for metagranite at Qinglongshan and adjacent areas in the southwestern end of the Sulu orogen. CL imaging suggests that the low δ¹⁸O zircons at Zaobuzhen are primarily of magmatic origin, but underwent different degrees of metamorphic modification. Zircon U-Pb dating yields middle Neoproterozoic ages of 751 ± 27 to 779 ± 25 Ma for protolith crystallization and Triassic ages of 214 ± 10 to 241 ± 33 Ma for metamorphic resetting. However, no metamorphic modification occurs in zircon REE patterns that only indicate magmatic recrystallization and hydrothermal alteration, respectively. Thus, the negative δ¹⁸O zircons are interpreted as crystallizing from negative δ¹⁸O magmas due to melting of meteoric-hydrothermally altered negative δ¹⁸O rocks in an active rift setting at about 780 Ma. The variation in zircon δ¹⁸O values indicates considerable O isotope heterogeneity in its granitic protolith. Zircon Lu-Hf isotope analyses give positive ε_Hf(t) values of 1.6-4.1 and Hf model ages of 1.18-1.30 Ga. This suggests that the granitic protolith was derived from the mid-Neoproterozoic reworking of late Mesoproterozoic juvenile crust. The metagranites at Zaobuzhen and Qinglongshan, about 450 km apart, are two known occurrences of the unusually low δ¹⁸O zircons below −6‰ so far reported in the Sulu orogen. They are similar to each other in both protolith and metamorphic ages, so that they share the same nature of both Neoproterozoic protolith and Triassic metamorphism. Therefore, the locally negative δ¹⁸O zircons may register centers of low δ¹⁸O magmatism during the supercontinental rifting.     

Bromine cycle in subduction zones through in situ Br monitoring in diamond anvil cells

The geochemical partitioning of bromine between hydrous haplogranitic melts, initially enriched with respect to Br and aqueous fluids, has been continuously monitored in situ during decompression. Experiments were carried out in diamond anvil cells from 890 °C to room temperature and from 1.7 GPa to room pressure, typically from high P, T conditions corresponding to total miscibility (presence of a supercritical fluid). Br contents were measured in aqueous fluids, hydrous melts and supercritical fluids. Partition coefficients of bromine were characterized at pressure and temperature between fluids, hydrous melts and/or glasses, as appropriate: D^Br_fluid/melt = (Br)_fluid/(Br)_melt, ranges from 2.18 to 9.2 ± 0.5 for conditions within the ranges 0.66-1.7 GPa, 590-890 °C; and D^Br_fluid/glass = (Br)_fluid/(Br)_glass ranges from 60 to 375 at room conditions. The results suggest that because high pressure melts and fluids are capable of accepting high concentrations of bromine, this element may be efficiently removed from the slab to the mantle source of arc magmas. We show that Br may be highly concentrated in subduction zone magmas and strongly enriched in subduction-related volcanic gases, because its mobility is strongly correlated with that of water during magma degassing. Furthermore, our experimental results suggest that a non negligible part of Br present in the subducted slab may remain in the down-going slab, being transported toward the transition zone. This indicates that the Br cycle in subduction zones is in fact divided in two related but independent parts: (1) a shallower one where recycled Br may leave the slab with a water and silica-bearing “fluid” leading to enriched arc magmas that return Br to the atmosphere. (2) A deeper cycle where Br may be recycled back to the mantle maybe to the transition zone, where it may be present in high pressure water-rich metasomatic fluids.     

Partitioning of Na between olivine and melt: An experimental study with application to the formation of meteoritic Na2O-rich chondrule glass and refractory forsterite grains

We report experimentally determined 1 atm olivine/melt D_Na partitioning data for low f_O2, a variety of melt compositions and a temperature range of 1325-1522 °C. We demonstrated that high-current electron microprobe analyses (EPMA, I = 500 nA, 600 s on the peak) allow quantitative determination of Na₂O in olivine down to ∼10 μg/g. The mean olivine/melt D_Na from 12 experimental runs is 0.0031 ± 0.0007 (1σ). This is the recommended value for low pressures and a wide range of natural compositions.This result is applied to the problem of the origin of alkalis in chondrules and the formation of chondritic refractory forsterite grains. The data on Semarkona (LL3.0) chondrules show that Na₂O is primordial and was present during olivine crystallization. For refractory forsterite grains from Murchison (CM2), we demonstrate that high CaO contents are not a result of equilibration with Na₂O-rich melts, but require high activities of CaO during their formation.     

He diffusion systematics in minerals: Evidence from synthetic monazite and zircon structure phosphates

Helium diffusivity was measured in synthetic rare-earth-element orthophosphates with systematically varying properties to evaluate potential controls on He transport in minerals. In the zircon structure phosphates (in this study, the phosphates of Tb, Dy, Ho, Er, Tm, Yb, and Lu as well as synthetic xenotime, YPO₄), He diffusion is strongly anisotropic. Transport apparently proceeds preferentially through channels aligned with the c-axis. The activation energy for diffusion is almost the same (122 ± 6 kJ/mol) in all members of this family, but there is a monotonic decrease in D_o with atomic number from TbPO₄ (∼10⁵ cm²/s) to LuPO₄ (∼10 cm²/s). The c-parallel channels become increasingly constricted in the same sequence, likely accounting for the systematically decreasing diffusivity. The He closure temperature (r = 1 cm, dT/dt = 10 °C/Myr) increases with atomic number from 44 °C for TbPO₄ to 88 °C for LuPO₄. Diffusion of radiogenic helium from natural zircon and xenotime is much slower than these synthetic analogs predict, suggesting that coupled substitution of REE and P for Zr and Si and/or radiation damage profoundly modify the energetics of interstitial He diffusion. In particular, α-recoil may play a key role by damaging the continuity and integrity of the channels.Monazite structure phosphates (here La, Ce, Pr, Nd, Sm, and Gd phosphate) are far more He retentive than those of the zircon structure. Activation energies increase smoothly with atomic number from LaPO₄ (183 kJ/mol) to NdPO₄ (224 kJ/mol) then decrease again to GdPO₄ (198 kJ/mol). D_o values mimic this pattern, spanning a range from ∼10⁻¹ cm²/s (GdPO₄) to 10⁴ cm²/s (NdPO₄). Nevertheless, He closure temperatures increase monotonically with atomic number, from 300 °C in LaPO₄ to 410 °C in GdPO₄. No evidence was obtained bearing on diffusion anisotropy, but the monazite structure lacks through-going channels so it is not expected. Diffusion parameters for radiogenic helium in natural monazite are similar to those obtained on the synthetic analogs.Ionic porosity is not the primary control on He diffusion in the orthophosphates. Within a given structure and with limited elemental substitution, ionic porosity and He closure temperature are negatively correlated, as predicted. However, differences between crystal structures are far more important than ion packing density: at comparable ionic porosity the monazite structure phosphates have He closure temperatures ∼300 °C higher than the xenotime structure phosphates. Modifications to the structures by radiation damage likely play a similarly significant role in controlling He diffusion.     

High-pressure phase equilibria and element partitioning experiments on Apollo 15 green C picritic glass: Implications for the role of garnet in the deep lunar interior

We report results of nominally anhydrous near-liquidus experiments on a synthetic analog to very low-titanium Apollo 15 green C lunar picritic glass from ∼2 to 5 GPa. Apollo 15 green C glass (A15C) is saturated with garnet and pyroxene on the liquidus at ∼3 GPa. However, such an assemblage is unlikely to represent the lunar-mantle source region for this glass, and instead an olivine + orthopyroxene-dominated source is favored, in accord with earlier lower-pressure experiments on A15C. Near-liquidus garnet has a slight but significant majorite component at ∼5 GPa in this iron-rich bulk composition, as expected from our previous work in ordinary-chondritic bulk compositions. Ion microprobe measurements of partitioning of Sr, Ba, Sc, Nd, Sm, Dy, Yb, Y, Zr, Hf, and Th between garnet and coexisting melt in these experiments are the first garnet partition coefficients (D values) available that are directly relevant to lunar compositions. D values for these garnets differ significantly compared to D values for garnets grown in more magnesian, terrestrial bulk compositions, which until now are all that have been available in modeling the possible role of garnet in the lunar interior. For example, D values for heavy rare earth elements are lower than are those from terrestrial basaltic systems. These partitioning values are well-described by the lattice-strain partitioning model, but predictive relationships for garnet partitioning using that model fail to match the measured values, as was the case in our earlier work on chondritic compositions. Using our new D values in place of the “terrestrial” values in a variety of models of lunar petrogenesis, we suggest that garnet is unlikely to be present in the source regions for very titanium-poor lunar liquids despite its appearance on the liquidus of A15C.     

In situ U-Pb, O and Hf isotopic compositions of zircon and olivine from Eoarchaean rocks, West Greenland: New insights to making old crust

The sources and petrogenetic processes that generated some of the Earth’s oldest continental crust have been more tightly constrained via an integrated, in situ (U-Pb, O and Hf) isotopic approach. The minerals analysed were representative zircon from four Eoarchaean TTG tonalites and two felsic volcanic rocks, and olivine from one harzburgite/dunite of the Itsaq Gneiss Complex (IGC), southern West Greenland. The samples were carefully chosen from localities with least migmatisation, metasomatism and strain. Zircon was thoroughly characterized prior to analysis using cathodoluminescence, scanning electron, reflected and transmitted light imaging. The zircon from all but one sample showed only minor post-magmatic recrystallisation. ²⁰⁷Pb/²⁰⁶Pb dating of oscillatory-zoned zircon using SHRIMP RG (n = 142) indicates derivation of the felsic igneous rocks from different batches of magma at 3.88, 3.85, 3.81, 3.80 and 3.69 Ga.Analyses of ¹⁸O/¹⁶O compositions of olivine from a harzburgite/dunite (n = 8) using SHRIMP II in multi-collector mode, indicate that the oxygen isotopic composition of this sample of Eoarchaean mantle (δ¹⁸O_Ol = 6.0 ± 0.4‰) was slightly enriched in ¹⁸O, but not significantly different from that of the modern mantle. Zircon δ¹⁸O measurements from the six felsic rocks (n = 93) record mean or weighted mean compositions ranging from 4.9 ± 0.7‰ to 5.1 ± 0.4‰, with recrystallised domains showing no indication of oxygen isotopic exchange during younger tectonothermal events. δ¹⁸O_Zr compositions indicate that the primary magmas were largely in equilibrium with the mantle or mantle-derived melts generated at similar high temperatures, while calculated tonalite δ¹⁸O_WR compositions (6.7-6.9‰) resemble those of modern adakites.LA-MC-ICPMS zircon ¹⁷⁶Hf/¹⁷⁷Hf analyses were obtained from six samples (n = 122). Five samples record weighted mean initial ε_Hf compositions ranging from to 0.5 ± 0.6 to −0.1 ± 0.7 (calculated using λ¹⁷⁶Lu = 1.867 × 10⁻¹¹ yr⁻¹), while one sample records a composition of 1.3 ± 0.7, indicating the magmas were generated from a reservoir with a time averaged, near chondritic Lu/Hf. The derivation of TTG magmas from a chondritic Lu/Hf source implies either that there was not voluminous continental crustal growth nor major mantle differentiation leading to Lu/Hf fractionation during the Hadean or Eoarchaean, or alternatively that rapid recycling of an early formed crust allowed the early mantle to maintain a chondritic Lu/Hf.Previous studies have demonstrated that ancient TTG rocks were mostly produced by dehydration melting of mafic rocks within the stability field of garnet, probably in flatly-subducted or buried oceanic crust. The oxygen isotopic signatures measured here at high spatial resolution allow the source materials to be better defined. Melting of a mixed mafic source consisting of ∼80% unaltered gabbro (δ¹⁸O_WR = 5.5‰) with ∼20% hydrothermally altered gabbro/basalt (δ¹⁸O_WR = 4.0‰) would produce tonalite magmas within the average compositional range observed. ¹⁸O-enriched components such as altered shallow basaltic oceanic crust and pelagic or continental sediments were not present in the sources of these TTG melts. The absence of high ¹⁸O signatures may indicate either the rarity of low temperature altered sediments, or their effective removal from the down-going slab.     

Ti site occupancy in zircon

Ti site occupancy in zircon (ZrSiO₄) is fundamental to thermobarometry because substitution mechanisms control Ti content-temperature relations. Here we describe the results of three independent methods used to demonstrate that Ti substitutes for Si and not Zr in zircon. Zircon grains were synthesized from oxide powders held in a Na₂WO₄ flux at 1 bar and 1300 °C. Zircon grains equilibrated with rutile + cristobalite show Ti contents (1201 ppm) nearly half that of zircon grains equilibrated with srilankite ((Ti,Zr)O₂) + tetragonal zirconia (2640 ppm). The lower Ti content of zircon grains produced at silica-saturated conditions indicates that Ti substitution predominately occurs on the Si site. Moreover, the higher Ti contents of silica-saturated experiments at 1 bar (1201 ppm), relative to those at 1 GPa (457 ppm, Ferry and Watson, 2007), indicates a substantial pressure effect on Ti solubility in zircon. Measured Ti K-α edge X-ray Absorption Near Edge Structure (XANES) spectra of synthetic zircon grains show energies and normalized intensities akin to those seen among tetrahedrally coordinated Ti-bearing standard minerals, strongly suggesting that Ti occupies the Si site. Density functional theory (DFT) calculations confirm that Ti substitution is most likely to occur on the Si site and predict a Ti-O bond length of 1.797 Å (compared to an average of 2.160 Å for substitution on the Zr site), in excellent agreement with X-ray Absorption Fine Structure (EXAFS) spectra of experimentally grown zircon grains which indicate a value of 1.76(1) Å. The software FEFF 8.4 was used to simulate XANES spectra from the defect structures determined by DFT for Ti substituting on both the Si and Zr sites. The predicted spectrum for Ti on the Si site reproduces all the key features of the experimental zircon spectra, whereas Ti on the Zr site is markedly different. All applied methods confirm that Ti substitutes for Si in zircon. Consequently, the Ti content of zircon at a given pressure is not only a function of temperature, but will increase with decreasing silica activity. Because elements that activate or quench cathodoluminescence (CL) in zircon are incorporated into the Zr site, a decoupling of CL from Ti contents - incorporated on the Si site in zircon is expected. This hypothesis has been verified by a systematic CL-trace element study of natural and experimental zircon.     

O/O and D/H ratios of pedogenic kaolinite in a North American Cenomanian laterite: Paleoclimatic implications

Kaolinite, gibbsite and quartz are the dominant minerals in samples collected from two outcrops of a Cenomanian (∼95 Ma) laterite in southwestern Minnesota. A combination of measured yields and isotope ratios permitted mass balance calculations of the δD and δ¹⁸O values of the kaolinite in these samples. These calculations yielded kaolinite δD values of about −73‰ and δ¹⁸O values of about +18.7‰. The δD and δ¹⁸O values appear to preserve information on the ancient weathering system.If formed in hydrogen and oxygen isotope equilibrium with water characterized by the global meteoric water line (GMWL), the kaolinite δD and δ¹⁸O values indicate a crystallization temperature of 22 (±5) °C. A nominal paleotemperature of 22 °C implies a δ¹⁸O value for the corresponding water of −6.3‰. The combination of temperature and meteoric water δ¹⁸O values is consistent with relatively intense rainfall at that mid-paleolatitude location (∼40°N) on the eastern shore of the North American Western Interior Seaway. The inferred Cenomanian paleosol temperature of ∼22 °C is in general accord with published mid-Cretaceous continental mean annual temperatures (MAT) estimated from leaf margin analyses of fossil plants.When compared with results from a published GCM-based Cenomanian climate simulation which specifies a latitudinal sea surface temperature (SST) gradient that was either near modern or smaller-than-modern, the kaolinite paleotemperature of 22 °C is closer to the GCM-predicted MAT for a smaller equator-to-pole temperature difference in the mid-Cretaceous. Moreover, the warm, kaolinite-derived, mid-paleolatitude temperature of 22 °C is associated with proxy estimates of high concentrations of atmospheric CO₂ in the Cenomanian. The overall similarity of proxy and model results suggests that the general features of Cenomanian continental climate in that North American locale are probably being revealed.     

Zircons from kimberlite: New insights from oxygen isotopes, trace elements, and Ti in zircon thermometry

Zircons found in mantle-sourced kimberlite provide probes into the isotopic chemistry of the asthenosphere and subcontinental lithospheric mantle. However, little is known about the conditions of formation of these zircons. A suite of 88 zircons found in kimberlites from Africa, Siberia, Brazil, and the United States have been analyzed for their Ti concentration and selected zircons were analyzed for their Rare Earth Element (REE) concentrations by ion microprobe. In addition, precise and accurate laser-fluorination oxygen isotope data were obtained for zircons from Brazil (5.1 ± 0.3‰, 1SD) and the Midwest United States (5.3 ± 0.3‰), yielding mantle-like δ¹⁸O values similar to published data for Africa (5.2 ± 0.3‰) and Siberia (5.3 ± 0.2‰). Most megacrysts in this study preserve fine-scale, oscillatory zoning in CL and are generally homogenous in oxygen isotopic composition, consistent with preservation of primary compositions. A few zircons from Brazil show some evidence of chemical zoning due to recrystallization. The Ti content of mantle zircon is in general low with average compositions from each locality of 13 ± 8.4 ppm (1SD, Kaapvaal craton), 12 ± 8.7 ppm (Siberian platform), 18 ± 11 ppm (Brazil), and 4.8 ±4.3 ppm (United States). The recently calibrated Ti in zircon thermometer yields an average temperature of 744 ±62 °C (1SD) for the average of 13 ± 9 ppm Ti, with no correction for pressure, a_TiO2, or a_SiO2. The Ti content of zircons found within rutile nodules from the Orapa kimberlite (Kaapvaal craton) is almost indistinguishable from those with no constraint on a_TiO2, suggesting that reduced a_TiO2 is not responsible for lower than expected mantle temperatures. The average temperature in this study corresponds to ∼3 GPa on a 40 mW/m² cratonic geotherm. If correct, this would suggest that zircon megacrysts from all four cratons formed in the shallow lithospheric mantle. However, there are several possibly confounding effects to this thermometer, including: a pressure correction and disequilibrium zircon growth. Zircons from rutile nodules have REE contents that span the range of mantle zircon REE and are similar to both zircon megacrysts and zircons from metasomatic assemblages.     

Trace element partitioning between apatite and silicate melts

We present new experimental apatite/melt trace element partition coefficients for a large number of trace elements (Cs, Rb, Ba, La, Ce, Pr, Sm, Gd, Lu, Y, Sr, Zr, Hf, Nb, Ta, U, Pb, and Th). The experiments were conducted at pressures of 1.0 GPa and temperatures of 1250 °C. The rare earth elements (La, Ce, Pr, Sm, Gd, and Lu), Y, and Sr are compatible in apatite, whereas the larger lithophile elements (Cs, Rb, and Ba) are strongly incompatible. Other trace elements such as U, Th, and Pb have partition coefficients close to unity. In all experiments we found D_Hf > D_Zr, D_Ta ≈ D_Nb, and D_Ba > D_Rb > D_Cs. The experiments reveal a strong influence of melt composition on REE partition coefficients. With increasing polymerisation of the melt, apatite/melt partition coefficients for the rare earth elements increase for about an order of magnitude. We also present some results in fluorine-rich and water-rich systems, respectively, but no significant influence of either H₂O or F on the partitioning was found. Furthermore, we also present experimentally determined partition coefficients in close-to natural compositions which should be directly applicable to magmatic processes.     

Crystal-melt partitioning of noble gases (helium, neon, argon, krypton, and xenon) for olivine and clinopyroxene

Mineral-melt partition coefficients of all noble gases (^min/meltD_i) have been obtained for olivine (ol) and clinopyroxene (cpx) by UV laser ablation (213 nm) of individual crystals grown from melts at 0.1 GPa mixed noble gas pressure. Experimental techniques were developed to grow crystals virtually free of melt and fluid inclusions since both have been found to cause profound problems in previous work. This is a particularly important issue for the analysis of noble gases in crystals that have very low partition coefficients relative to coexisting melt and fluid phases. The preferred partitioning values obtained for the ol-melt system for He, Ne, Ar, Kr, and Xe are 0.00017(13), 0.00007(7), 0.0011(6), 0.00026(16), and , respectively. The respective cpx-melt partition coefficients are 0.0002(2), 0.00041(35), 0.0011(7), 0.0002(2), and . The data confirm the incompatible behaviour of noble gases for both olivine and clinopyroxene but unlike other trace elements these values show little variation for a wide range of atomic radius. The lack of dependence of partitioning on atomic radius is, however, consistent with the partitioning behaviour of other trace elements which have been found to exhibit progressively lower dependence of ^min/meltD_i on radius as the charge decreases. As all noble gases appear to exhibit similar ^min/meltD_i values we deduce that noble gases are not significantly fractionated from each other by olivine and clinopyroxene during melting and fractional crystallisation. Although incompatible, the partitioning values for noble gases also suggest that significant amounts of primordial noble gases may well have been retained in the mantle despite intensive melting processes. The implication of our data is that high primordial/radiogenic noble gas ratios (³He/⁴He, ²²Ne/²¹Ne, and ³⁶Ar/⁴⁰Ar) characteristic of plume basalt sources can be achieved by recycling a previously melted (depleted) mantle source rather than reflecting an isolated, non-degassed primordial mantle region.     

Response of Precambrian zircon to the chemical abrasion (CA-TIMS) method and implications for improvement of age determinations

Electron backscatter imaging, Raman spectroscopy and U-Pb geochronology have been applied to Precambrian zircon grains that were annealed at 1000 and 1450 °C for various times, then leached with HF to constrain the conditions for healing radiation damage and attaining primary U-Pb zircon ages using the chemical abrasion (CA-TIMS) method. SEM images reveal a variety of textures for ZrO₂ overgrowths on 1450 °C annealed and leached zircon surfaces that depend on the degree of radiation damage and annealing history. Highly damaged zircon produces finer textures than zircon with less damage.Raman spectroscopy indicates that crystals with different levels of radiation damage are only partially restored by annealing at 1000 °C for 2-3 days. Longer annealing periods of 20 days are not noticeably more effective. Annealing at 1450 °C for 1 h results in partial breakdown of zircon but restores Raman peak widths and wave numbers to values characteristic of undamaged zircon after ZrO₂ overgrowths are removed by HF. Raman spectra are much less sensitive to polarization angle for annealed highly damaged grains than for weakly damaged zircon.U-Pb isotopic analyses of low to moderately damaged zircon (alpha fluence ranging up to 10¹⁹/g corresponding to an amorphization volume fraction of 80% or more) yield almost concordant data (0.3-0.5% discordance) after high-temperature annealing at 1450 °C followed by HF leaching at 195 °C. Analyses of cracked zircon annealed at 1450 °C and leached may remain discordant but those of uncracked grains are concordant. Most analyses show primary ²⁰⁷Pb/²⁰⁶Pb ages although cracked grains annealed at 1450 °C may produce discordant data with ²⁰⁷Pb/²⁰⁶Pb ages that are too young after leaching. The solubility of highly damaged, very disordered zircon (amorphization level of 99%) is only slightly reduced by annealing, and analyses of leach residues are strongly discordant although primary ²⁰⁷Pb/²⁰⁶Pb ages are obtained.Annealing of highly damaged zircon under any conditions apparently results in a mass of randomly oriented micro-crystals that pseudomorph the original grain. This could explain the fine-scale pattern observed on etched crystal surfaces, reduced anisotropy at the 5 μm scale of the Raman laser beam and high solubility in HF. It may be impossible to restore primary U-Pb isotopic ages in such cases but precise ages can still potentially be determined from ²⁰⁷Pb/²⁰⁶Pb ratios or by application of the air abrasion method.     

Experimentally determined biomediated Sr partition coefficient for dolomite: Significance and implication for natural dolomite

Two strains of moderately halophilic bacteria were grown in aerobic culture experiments containing gel medium to determine the Sr partition coefficient between dolomite and the medium from which it precipitates at 15 to 45 °C. The results demonstrate that Sr incorporation in dolomite does occur not by the substitution of Ca, but rather by Mg. They also suggest that Sr partitioning between the culture medium and the minerals is better described by the Nernst equation (D_Sr^dol = Sr_dol/Sr_bmi), instead of the Henderson and Kracek equation (D_Sr^dol = (Sr/Ca)_dol/(Sr/Ca)_solution. The maximum value for D_Sr^dol occurs at 15 °C in cultures with and without sulfate, while the minimum values occur at 35 °C, where the bacteria exhibit optimal growth. For experiments at 25, 35 and 45 °C, we observed that D_Sr^dol values are greater in cultures with sulfate than in cultures without sulfate, whereas D_Sr^dol values are smaller in cultures with sulfate than in cultures without sulfate at 15 °C.Together, our observations suggest that D_Sr^dol is apparently related to microbial activity, temperature and sulfate concentration, regardless of the convention used to assess the D_Sr^dol. These results have implications for the interpretation of depositional environments of ancient dolomite. The results of our culture experiments show that higher Sr concentrations in ancient dolomite could reflect microbial mediated primary precipitation. In contrast, previous interpretations concluded that high Sr concentrations in ancient dolomites are an indication of secondary replacement of aragonite, which incorporates high Sr concentrations in its crystal lattice, reflecting a diagenetic process.     

O and H diffusion in uraninite: Implications for fluid-uraninite interactions, nuclear waste disposal, and nuclear forensics

Diffusion coefficients for oxygen and hydrogen were determined from a series of natural uraninite-H₂O experiments between 50 and 700 °C. Under hydrous conditions there are two diffusion mechanisms: (1) an initial extremely fast-path diffusion mechanism that overprinted the oxygen isotopic composition of the entire crystals regardless of temperature and (2) a slower volume-diffusive mechanism dominated by defect clusters that displace or eject nearest neighbor oxygen atoms to form two interstitial sites and two partial vacancies, and by vacancy migration. Using the volume diffusion coefficients in the temperature range of 400-600 °C, diffusion coefficients for oxygen can be represented by D = 1.90e⁻⁵ exp (−123,382 J/RT) cm²/s and for temperatures between 100 and 300 °C the diffusion coefficients can be represented by D = 1.95e⁻¹⁰ exp (−62484 J/RT) cm²/s, where the activation energies for uraninite are 123.4 and 62.5 kJ/mol, respectively. Hydrogen diffusion in uraninite appears to be controlled by similar mechanisms as oxygen. Using the volume diffusion coefficients for temperatures between 50 and 700 °C, diffusion coefficients for hydrogen can be represented by D = 9.28e⁻⁶ exp (−156,528 J/RT) cm²/s for temperatures between 450 and 700 °C and D = 1.39e⁻¹⁴ exp (−34518 J/RT) cm²/s for temperatures between 50 and 400 °C, where the activation energies for uraninite are 156.5 and 34.5 kJ/mol, respectively.Results from these new experiments have implications for isotopic exchange during natural UO₂-water interactions. The exceptionally low δ¹⁸O values of natural uraninites (i.e. −32‰ to −19.5‰) from unconformity-type uranium deposits in Saskatchewan, in conjunction with theoretical and experimental uraninite-water and UO₃-water fractionation factors, suggest that primary uranium mineralization is not in oxygen isotopic equilibrium with coeval clay and silicate minerals. The low δ¹⁸O values have been interpreted as resulting from the low temperature overprinting of primary uranium mineralization in the presence of relatively modern meteoric fluids having δ¹⁸O values of ca. −18‰, despite petrographic and U-Pb isotope data that indicate limited alteration. Our data show that the anomalously low oxygen isotopic composition of the uraninite from the Athabasca Basin can be due to meteoric water overprinting under reducing conditions, and meteoric water or groundwater can significantly affect the oxygen isotopic composition of spent nuclear fuel in a geologic repository, with minimal change to the chemical composition or texture. Moreover, the rather fast oxygen and hydrogen diffusion coefficients for uraninite, especially at low temperatures, suggest that oxygen and hydrogen diffusion may impart characteristic isotopic signals that can be used to track the route of fissile material.     

Solubility behavior of alkali aluminosilicate components in aqueous fluids and silicate melts at high pressure and temperature

The solubility behavior of K₂O, Na₂O, Al₂O₃, and SiO₂ in silicate-saturated aqueous fluid and coexisting H₂O-saturated silicate melts in the systems K₂O-Al₂O₃-SiO₂-H₂O and Na₂O-Al₂O₃-SiO₂-H₂O has been examined in the 1- to 2-GPa pressure range at 1100°C. Glasses of Na- and K-tetrasilicate compositions with 0, 3, and 6 mol% Al₂O₃ were used as starting materials. In both systems, the oxides dissolve incongruently in aqueous fluid and silicate melt. When recalculated to an anhydrous basis, the aqueous fluids are enriched in alkalis and depleted in silica and alumina relative to their proportions in the starting materials. The extent of incongruency is more pronounced in the Na₂O-Al₂O₃-SiO₂-H₂O system than in the K₂O-Al₂O₃-SiO₂-H₂O system.The partition coefficients of the oxides, D_oxide^fluid/melt, are linear and positive functions of the oxide concentration in the fluid for each composition. There is a slight dependence of the partition coefficients on bulk composition. No effect of pressure could be discerned. For alkali metals, the fluid/melt partition coefficients range from 0.06 to 0.8. For Al₂O₃ this range is 0.01 to 0.2, and for SiO₂, it is 0.01 to 0.32. For all compositions, D_K2O^fluid/melt∼D_Na2O^fluid/melt>D_SiO2^fluid/melt>D_Al2O3^fluid/melt for the same oxide concentration in the fluid. D_K2O^fluid/melt, D_Na2O^fluid/melt, and D_SiO2^fluid/melt correlate negatively with the Al₂O₃ content of the systems. This correlation is consistent with a solubility model of alkalis that involve associated KOH°, NaOH°, silicate, and aluminate complexes.