Kinetics of zircon dissolution and zirconium diffusion in granitic melts of variable water content

The experimental dissolution of zircon into a zircon-undersaturated felsic melt of variable water content at high pressure in the temperature range 1,020° to 1,500° C provides information related to 1) the solubility of zircon, 2) the diffusion kinetics of Zr in an obsidian melt, and 3) the rate of zircon dissolution. Zirconium concentration profiles observed by electron microprobe in the obsidian glass adjacent to a large, polished zircon face provide sufficient information to calculate model diffusion coefficients. Results of dissolution experiments conducted in the virtual absence of water (<0.2% H₂O) yield an activation energy (E) for Zr transport in a melt ofM=1.3 [whereM is the cation ratio (Na+K+2Ca)/(Al·Si)] of 97.7±2.8 kcal-mol^?1, and a frequency factor (D ₀) of 980 _?580 ^+1,390 cm²-sec^?1. Hydrothermal experiments provide an E=47.3±1.9 kcal-mol^?1 andD ₀=0.030 _?0.015 ^+0.030 cm²-sec^?1. Both of these results plot close to a previously defined diffusion compensation line for cations in obsidian. The diffusivity of Zr at 1,200° C increases by a factor of 100 over the first 2% of water introduced into the melt, but subsequently rises by only a factor of five to an apparent plateau value of ～2×10^?9 cm²-sec^?1 by ～6% total water content. The remarkable contrast between the wet and dry diffusivities, which limits the rate of zircon dissolution into granitic melt, indicates that a 50 μm diameter zircon crystal would dissolve in a 3 to 6% water-bearing melt at 750° C in about 100 years, but would require in excess of 200 Ma to dissolve in an equivalent dry system. From this calculation we conclude that zircon dissolution proceeds geologically instantaneously in an undersaturated, water-bearing granite. Estimates of zircon solubility in the obsidian melt in the temperature range of 1,020° C to 1,500° C confirm and extend an existing model of zircon solubility to these higher temperatures in hydrous melts. However, this model does not well describe zircon saturation behavior in systems with less than about 2% water.     

Dislocation-assisted diffusion of oxygen in albite

Oxygen diffusion in albite has been determined by the integrating (bulk ¹⁸O) method between 750° and 450° C, for a P _H2O of 2 kb. The original material has a low dislocation density (<10⁶ cm^?2), and its lattice diffusion coefficient (D ₁), given below, agrees well with previous determinations. A sample was deformed at high temperature and pressure to produce a uniform dislocation density of 5 × 10⁹ cm^?2. The diffusion coefficient (D _a) for this deformed material, given below, is about 0.5 and 0.7 orders of magnitude larger than D ₁ at 700° and 450° C, respectively. This enhancement is believed due to faster diffusion along the cores of dislocations. Assuming a dislocation core radius of 4 Å, the calculated pipe diffusion coefficient (D _p), given below, is about 5 orders of magnitude larger than D ₁. These results suggest that volume diffusion at metamorphic conditions may be only slightly enhanced by the presence of dislocations. $$\begin{gathered} D_1 = 9.8 \pm 6.9 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 33.4 \pm 0.6(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_a = 7.6 \pm 4.0 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 30.9 \pm 1.1(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_p \approx 1.2 \times 10^{ - 1} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 29.8(kcal/mole)/RT]. \hfill \\ \end{gathered} $$      

Diffusion compensation in natural silicates

The temperature dependence of diffusion is usually found to follow the Arrhenius law: D = D₀e^?E/RT Winchell (1969) showed that there is commonly an inter-dependence between D₀ and E (for diffusion in silicate glasses), such that diffusion of different species show a positive correlation on a log D₀ vs E plot. A similar effect was noted by Hofmann (1980) for cation diffusion in basalt. This implies that diffusion rates of different species tend to converge at a particular temperature; this effect is known as the ‘compensation effect’. I will show that this effect is also present for diffusion in feldspars and olivines. The equations for the compensation lines (with E given in kcal/mol) are: basalt—E = 50 + 7.5 log D₀ feldspar—E = 50.7 + 3.4 log D₀ olivine—E = 78.0 + 7.5 log D₀ The convergence, or crossover, temperatures for diffusion in various materials are: obsidian—3400°C basalt—1370°C olivine—1360°C feldspar—460°C Compensation plots are useful for evaluating and comparing experimental diffusion data (though of limited usefulness in a predictive sense) and for understanding ‘closure temperatures’ for diffusion in petrogenetic processes (since closure temperature, the temperature at which natural diffusion processes are frozen in, is dependent on E, log d₀, and cooling rate). I show that most diffusing species in feldspar have a closure-temperature close to the crossover or convergence temperature, implying that all species in feldspars can be expected to ‘freeze-in’ simultaneously at temperatures in the range 400–600°C (for cooling rates in the range 10¹–10⁵°C/myr). Closure temperatures of various species in olivine, on the other hand, span a much larger range (800°C) for a similar range in cooling rates, implying that different elements in olivine will record different time-temperature stages in petrogenetic processes.     

Surface chemistry and dissolution kinetics of glassy rocks at 25°C

The weathering rates and mechanisms of three types of glassy rocks were investigated experimentally at 25 °C, pH 1.0 to 6.2, and reaction times as much as to 3 months. Changes in major element chemistry were monitored concurrently as a function of time in the aqueous solution and within the near surface region of the glass. Leach profiles, obtained by a HF leaching technique, displayed near-surface zones depleted in major cations. These zones increased in depth with increasing time and decreasing pH of reactions. Release rates into the aqueous solution were parabolic for Na and K and linear for Si and Al. A coupled weathering model, involving surface dissolution with concurrent diffusion of Na, K, and Al, produced a mass balance between the aqueous and glass phases. Steady state conditions are reached at pH 1.0 after approximately 3 weeks of reaction. Steady-state is not reached even after 3 months at pH 6.2.An interdiffusion model describes observed changes in Na diffusion profiles for perlite at pH 1.0. The calculated Na self-diffusion coefficient of 5 × 10^?19 cm²·s^?1 at 25°C approximates coefficients extrapolated from previously reported high temperature data for obsidian. The self-diffusion coefficient for H₃O⁺, 1.2 × 10^?20 cm²·s^?1, is similar to measured rates of water diffusion during hydration of obsidian to form perlite.     

An XPS study of the dissolution kinetics of chrysotile in 0.1 N oxalic acid at different temperatures

The dissolution of chrysotile is studied in regard to the surfaces analysis by photoelectron spectrometry. After leaching of chrysotile (Provenance: Thetford; about 200 mg of fibers of 1 cm length) in nonstirred 0.1 N oxalic conditions, the composition of the mineral surfaces is determined by XPS; kinetic curves of dissolution are given in the range 22–80°C. Two conditions for the rate-limiting step are involved for the explanation of the dissolution: diffusion of Mg²⁺ through the fibrous gel or dissociation of chrysotile. By the former, some values of the diffusion coefficient are proposed: D varies from 5·10^?19 cm²s^?1 to 5·10^?16 cm²s^?1, in the range 22–80°C. By the second model, the leaching rate is estimated from 3 Å (22°C) per h to 250 Å (80°C) per h. For the 2 models, the activation heat energy is in the range 15–20 Kcal.     

Oxygen diffusion in amphiboles

The kinetics of oxygen isotope self-diffusion in natural samples of hornblende, tremolite, and richterite have been measured. Samples were run under hydrothermal conditions using ¹⁸O enriched water. Profiles of ${}^{18}\text{O}\text{(}{}^{16}\text{O +}{}^{18}\text{O)}$vs depth into the crystal were obtained using an ion microprobe; the depths of sputtered holes were measured using an optical interferometer. At 1000 bars (100 MPa) water pressure, the activation energies (Q) and pre-exponential factors (D₀) for diffusion parallel to c are: D₀(cm²/sec) Q (kcal/gm-atom) T (°C) Hornblende 1⁺²⁰_?1 × 10^?741 ± 6 650–800 Tremolite 2⁺³⁰_?2× 10^?8 39 ± 5 650–800 Richterite 3⁺⁵_?2 × 10^?4 57 ± 2 650–800The diffusion coefficient (D) for hornblende at 800°C and 1000 bars water pressure measured parallel to the c crystallographic direction is at least ten times greater than that parallel to the a or b directions. An increase in water pressure from 200 to 2000 bars increases D by a factor of 2.7 for diffusion parallel to c at 800°C. The D value for hornblende at 800°C is about 0.01 that for quartz and 0.001 that for anorthite. As a result, closure temperatures for oxygen exchange in natural primary amphiboles are significantly higher than for quartz or feldspars. It is unlikely that amphiboles will exchange oxygen isotopes by diffusion under most crustal conditions.     

Helium measurements of pore fluids obtained from the San Andreas Fault Observatory at Depth (SAFOD, USA) drill cores

⁴He accumulated in fluids is a well established geochemical tracer used to study crustal fluid dynamics. Direct fluid samples are not always collectable; therefore, a method to extract rare gases from matrix fluids of whole rocks by diffusion has been adapted. Helium was measured on matrix fluids extracted from sandstones and mudstones recovered during the San Andreas Fault Observatory at Depth (SAFOD) drilling in California, USA. Samples were typically collected as subcores or from drillcore fragments. Helium concentration and isotope ratios were measured 4?C6 times on each sample, and indicate a bulk ⁴He diffusion coefficient of 3.5?±?1.3?×?10^?C8 cm²?s^?C1 at 21°C, compared to previously published diffusion coefficients of 1.2?×?10^?C18 cm²?s^?C1 (21°C) to 3.0?×?10^?C15 cm²?s^?C1 (150°C) in the sands and clays. Correcting the diffusion coefficient of ⁴He_water for matrix porosity (??3%) and tortuosity (??6?C13) produces effective diffusion coefficients of 1?×?10^?C8 cm^2?s^?C1 (21°C) and 1?×?10^?C7 (120°C), effectively isolating pore fluid ⁴He from the ⁴He contained in the rock matrix. Model calculations indicate that <6% of helium initially dissolved in pore fluids was lost during the sampling process. Complete and quantitative extraction of the pore fluids provide minimum in situ porosity values for sandstones 2.8?±?0.4% (SD, n?=?4) and mudstones 3.1?±?0.8% (SD, n?=?4).     

Inter-diffusion coefficients parallel to the C-Axis in iron-rich clinopyroxenes calculated from microstructures

Subsolidus marginal zoning in calcium-poor clinopyroxenes and intermediate zoning in discontinuously zoned subcalcic- to calcium-rich clinopyroxenes from ironrich igneous rocks is used to calculate the interdiffusion coefficient, D_Ca?(Fe,Mg), parallel to the crystallographic caxis. Wagner's mathematical models describing the displacement of interfaces in solids as the result of isothermal diffusion are adopted. The steady-state heat flow equation is used to approximate the diffusion times. The calculated interdiffusion coefficients are of a reasonable order of magnitude, viz. 6.0×10^?20?2.0×10^?17cm²· sec^?1 at about 900° C.     

Argon retention properties of silicate glasses and implications for <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar age and noble gas diffusion studies

Sized aggregates of glasses (47–84 wt% SiO₂) were fused from igneous-derived cohesive fault rock and igneous rock, and step-heated from ~400 to >1,200 °C to obtain their ³⁹Ar diffusion properties (average E=33,400 cal mol^?1; D _o=4.63×10^?3 cm² s^?1). At T<~1,000 °C, glasses containing <~69 wt% SiO₂ and abundant network-forming cations (Ca, Fe, Mg) reveal moderate to strong non-linear increases in D and E, reflecting structural modifications as the solid transitions to melt. Extrapolation of these Arrhenius properties down to typical geologic T-t conditions could result in a 1.5 log₁₀ unit underestimation in the diffusion rate of Ar in similar materials. Numerical simulations based upon the diffusion results caution that some common geologic glasses will likely yield ⁴⁰Ar/³⁹Ar cooling ages rather than formation ages. However, if cooling rates are sufficiently high, ambient temperatures are sufficiently low (e.g., <65–175 °C), and coarse particles (e.g., radius (r) >~1 mm) are analyzed, glasses with compositions similar to ours may preserve their formation ages.     

Compositional zoning in garnet and kinetics of metamorphic crystallization

The kinetics of zoned garnet porphyroblast growth is exemplified in a sample of garnet-staurolite-biotite schist from the northern Ladoga region. The diffusion-controlled porphyroblast growth was accompanied by a decrease in the kinetic coefficient during phase reactions. Even at insignificant (1–2°C) thermal overstepping, the leading role of diffusion as a factor that controls kinetics of porphyroblast growth in medium-grade metapelites is consistent with the parameters of metamorphic crystallization: T = 500–650°C, t = 1 Ma; D _A1 ^app = 10^?14 cm²/s, L = 0.2–0.6 cm, r = 1–3 mm, ΔC _Al = 1.5 × 10^?4–1.5 × 10^?3 mol/cm³.     

The behavior of apatite during crustal anatexis: Equilibrium and kinetic considerations

The solubility and dissolution kinetics of apatite in felsic melts at 850°–1500°C have been examined experimentally by allowing apatite crystals to partially dissolve into apatite-undersaturated melts containing 0–10 wt% water. Analysis of P and Ca gradients in the crystal/melt interfacial region enables determination of both the diffusivities and the saturation levels of these components in the melt. Phosphorus diffusion was identified as the rate-limiting factor in apatite dissolution. Results of four experiments at 8 kbar run in the virtual absence of water yield an activation energy (E) for P diffusion of 143.6 ± 2.8 kcal-mol^?1 and frequency factor (D₀) of 2.23^+2.88_?1.26 × 10⁹cm²-sec^?1. The addition of water causes dramatic and systematic reduction of both E and D₀ such that at 6 wt% H₂O the values are ~25 kcal-mol^?1 and 10^?5 cm²-sec^?1, respectively. At 1300°C, the diffusivity of P increases by a factor of 50 over the first 2% of water added to the melt, but rises by a factor of only two between 2 and 6%, perhaps reflecting the effect of a concentration-dependent mechanism of H₂O solution. Calcium diffusion gradients do not conform well to simple diffusion theory because the release of calcium at the dissolving crystal surface is linked to the transport rate of phosphorus in the melt, which is typically two orders of magnitude slower than Ca. Calcium chemical diffusion rates calculated from the observed gradients are about 50 times slower than calcium tracer diffusion.Apatite solubilities obtained from these experiments, together with previous results, can be described as a function of absolute temperature (T) and melt composition by the expression: In D^apatite/melt_P = [(8400 + ((SiO₂ ? 0.5)2.64 × 10⁴))/T] ? [3.1 + (12.4(SiO₂ ? 0.5))] where SiO₂ is the weight fraction of silica in the melt. This model appears to be valid between 45% and 75% SiO₂, 0 and 10% water, and for the range of pressures expected in the crust.The diffusivity information extracted from the experiments can be directly applied to several problems of geochemical interest, including I) dissolution times for apatite during crustal anatexis, and 2) pileup of P, and consequent local saturation in apatite, at the surfaces of growing major-mineral phases.     

The effect of dislocations on bulk diffusion in feldspars during metamorphism

ABSTRACT There is no significant difference in the diffusion profiles across albite-adularia bicrystals that were simultaneously deformed at a strain rate of 10^-6S^-1 and those from hydrostatic experiments at the same conditions (1500 MPa and 1000°C for 156 h). This indicates that the bulk alkali diffusion rate, which is the sum of lattice diffusion (D, ₁) and dislocation pipe diffusion (Dp), is not significantly enhanced by dislocations at these conditions, and that the maximum value for the ratio of D_p/D₁ is about 10⁵. This is equal to the value previously reported for‘oxygen’diffusion in albite. If this ratio is independent of temperature, the contribution of either static (pre-deformed) or moving (syn-deformed) dislocations to the bulk diffusion rate of alkalis is probably minor at all metamorphic conditions. For Al and Si diffusion the ratio of D_p/D₁ may be larger if D₁ is lower. Thus a significant contribution from dislocations to bulk diffusion cannot be ruled out, especially during simultaneous deformation.     

Concentrations and behavior of oxygen and oxide ion in melts of composition CaO · MgO · xSiO2

The concentrations and behavior of oxygen and oxide ion were studied in silicate melts of composition CaO · MgO · xSiO₂ (1.25 ≤ x ≤ 3) in the temperature range 1425 to 1575°C by cyclic voltammetry and chronopotentiometry. Electroreduction of oxygen is a reversible, 2 electron process involving dissociated oxygen atoms. The Henry's Law constant for O₂ in molten diopside (CaO · MgO · 2SiO₂) is 0.023 ± 0.004 mole/l atm at 1450°C. The diffusion coefficient for molecular oxygen in diopside melt is 4.5 ± .5 × 10^?6 cm²/sec at 1450°C and the activation energy of diffusion is 80 ± 2 kcal/mole. Oxide ions produced by electroreduction of oxygen, rapidly dissociate silicate polymers, causing the concentration of free oxide ions in diopside melt to be buffered at a low level (4.7 ± .8 × 10^?5 mole/l). The concentration of free oxide ion increases at higher proportions of metal oxides but remains at this value in more silicic melts. The rate of formation of oxide ions by polymerization in diopside melt is 0.021 ± .007 mole/l sec. Thermodynamic parameters (the standard free energy, enthalpy and entropy) for the oxidation of Ni, Co, and Zn in diopside melt in equilibrium with gaseous oxygen agree with those for solid oxide systems. The platinum reference electrode in molten diopside is a reversible, oxygen electrode.     

Co-precipitation of Sr2+ and Ba2+ with aragonite by membrane diffusion of CO2 between 10 and 50 °C

Aragonite is precipitated by a new CO₂-diffusion technique from a Ca²⁺–Mg²⁺–Cl⁻ solution between 10 and 50 °C. Crystallisation of aragonite instead of calcite occurs by maintaining a [Mg²⁺]/[Ca²⁺] ratio of 2 in the fluid. The dissolved inorganic carbon (DIC) is received by diffusion of CO₂ through a polyethylene membrane (diffusion coefficient: D_CO2=10^−6.4 cm² s⁻¹ at 19 °C). It is suggested that significant amounts of DIC may be transferred by diffusion of CO₂ in natural systems if the CO₂ gradient is high. The CO₂-diffusion technique is used as a kind of simple mixed flow reactor for the co-precipitation of barium and strontium with aragonite. The distribution coefficients of Ba²⁺ and Sr²⁺ decrease from 10 to 50 °C according to D_Ba,a*=2.42−0.03595T (°C) and D_Sr,a*=1.32−0.005091T (°C). At 25 °C, the distribution coefficients are D_Ba,a*=1.5±0.1 and D_Sr,a*=1.19±0.03. The effect of temperature on D_Ba,a* is about one order of magnitude higher versus that on D_Sr,a*. Thus, Ba²⁺ may be a potential paleotemperature indicator if the composition of the solution is known.     

Crystal structure, Raman and FTIR spectroscopy, and equations of state of OH-bearing MgSiO3 akimotoite

MgSiO₃ akimotoite is stable relative to majorite-garnet under low-temperature geotherms within steeply or rapidly subducting slabs. Two compositions of Mg–akimotoite were synthesized under similar conditions: Z674 (containing about 550 ppm wt H₂O) was synthesized at 22 GPa and 1,500 °C and SH1101 (nominally anhydrous) was synthesized at 22 GPa and 1,250 °C. Crystal structures of both samples differ significantly from previous studies to give slightly smaller Si sites and larger Mg sites. The bulk thermal expansion coefficients of Z674 are (153–839 K) of a ₁ = 20(3) × 10^?9 K^?2 and a ₀ = 17(2) × 10^?6 K^?1, with an average of α ₀ = 27.1(6) × 10^?6 K^?1. Compressibility at ambient temperature of Z674 was measured up to 34.6 GPa at Sector 13 (GSECARS) at Advanced Photon Source Argonne National Laboratory. The second-order Birch–Murnaghan equation of state (BM2 EoS) fitting yields: V ₀ = 263.7(2) Å³, K _T0 = 217(3) GPa (K′ fixed at 4). The anisotropies of axial thermal expansivities and compressibilities are similar: α _a  = 8.2(3) and α _c  = 10.68(9) (10^?6 K^?1); β _a  = 11.4(3) and β _c  = 15.9(3) (10^?4 GPa). Hydration increases both the bulk thermal expansivity and compressibility, but decreases the anisotropy of structural expansion and compression. Complementary Raman and Fourier transform infrared (FTIR) spectroscopy shows multiple structural hydration sites. Low-temperature and high-pressure FTIR spectroscopy (15–300 K and 0–28 GPa) confirms that the multiple sites are structurally unique, with zero-pressure intrinsic anharmonic mode parameters between ?1.02 × 10^?5 and +1.7 × 10^?5 K^?1, indicating both weak hydrogen bonds (O–H···O) and strong OH bonding due to long O···O distances.     

Interdiffusion of Mg–Fe in olivine at 1,400–1,600 °C and 1 atm total pressure

The interdiffusion coefficient of Mg–Fe in olivine (D _Mg–Fe) was obtained at 1,400–1,600 °C at the atmospheric pressure with the oxygen fugacity of 10^?3.5–10^?2 Pa using a diffusion couple technique. The D _Mg–Fe shows the anisotropy (largest along the [001] direction and smallest along the [100] direction), and its activation energy (280–320 kJ/mol) is ~80–120 kJ/mol higher than that estimated at lower temperatures. The D _Mg–Fe at temperatures of >1,400 °C can be explained by the cation-vacancy chemistry determined both by the Fe³⁺/Fe²⁺ equilibrium and by the intrinsic point defect formation with the formation enthalpy of 220–270 kJ/mol depending on the thermodynamical model for the Fe³⁺/Fe²⁺ equilibrium in olivine. The formation enthalpy of 220–270 kJ/mol for the point defect (cation vacancy) in olivine is consistent with that estimated from the Mg self-diffusion in Fe-free forsterite. The increase in the activation energy of D _Mg–Fe at >1,400 °C is thus interpreted as the result of the transition of diffusion mechanism from the transition metal extrinsic domain to the intrinsic domain at the atmospheric pressure.     

Experimental deformation of sintered albite above and below the order-disorder transition

Abstract

Purified albite powder (44-53 μm) has been sintered to form an albite polycrystal suitable for deformation studies close to the melting temperature. Experiments have been carried out in Griggs solid medium deformation apparatus at 800, 1020 °C and 700 MPa pressure in a dehydrating pyrophyllite confining medium at constant strain rates of 10^?4, 10^?5, 10^?6 and 10^?7/s. At 800 °C the samples were brittle-ductile whereas at 1020 °C they were ductile with a rheology well described by a power law with a stress exponent of 3. The transition from brittle-ductile to ductile also coincided with the order-disorder or low-high albite transition, as indicated by the marked increase in mechanical twinning on the albite law at high temperature. At 1020 °C high dislocation densities (10^?10 - 10^?11/cm² and mechanical twinning characterised the original high albite grains, whereas fine recrystallised grains ( < 5 цm) had low dislocation densities (10⁷ - 10⁸/cm²) and often contained polysynthetic albite and pericline (M-type) twins. It is suggested that the recrystallized grains were monalbite (monoclinic) under test conditions which have inverted to high albite (triclinic) and in so doing produced M-twins, and that the recrystallization mechanism involved grain-boundary bulging to nucleate new high angle boundaries. The implications of the order-disorder transition for twinning and grain boundary migration are discussed and it is suggested that the data cannot be simply extrapolated to natural deformation in the low albite field.     

4He diffusion in specular hematite

In order to test the chronometer qualities of speculante for the (U + Th)/He dating method, ⁴He release experiments by stepwise heating of two specularites from the Rimbach mineralization locality in the southern Vosgues (France) have been carried out. The diffusion coefficients define linear Arrhenius plots within a temperature interval of 250 to 830 °C, which is suggestive of volume diffusion. Extrapolation of the diffusion behavior to 20° C yields diffusion coefficients (D₂₀ values) smaller than 10^?26 [cm² s^?1] for both hematites with activation energies at 116 [kJ/mole]. The results of our study suggest that specularite is a very helium retentive hematite variety which is capable of quantitatively retaining radiogenic helium over geologic periods of time.     

Carbon in olivine single crystals analyzed by the 12C(d,p)13C method and by photoelectron spectroscopy

Carbon subsurface concentration profiles in olivine single crystals from San Carlos, Arizona, and the Sergebet Island. Red Sea, containing total carbon between 60–180 wt.-ppm, were analyzed by means of the ¹²C(d. p)¹³C nuclear reaction and by x-ray induced photoelectron spectroscopy (XPS) in combination with acid etching and with Ar⁺ ion sputtering respectively, between 200–930 K. The (d, p) analysis reveals equilibrium subsurface C profiles extending 1–2 μm or more into the bulk. Their steepness is a function of temperature. Typical mean C concentrations at 300 K in the resolvable layers, 0–0.6, 0.6–1.2, and 1.2–1.8 μm. are 1.8, and 0.6 wt.-%, corresponding to enrichment factors over the mean bulk C concentration of the order of 100, 40 and 30 respectively. In the topmost atomic layers analyzed by XPS the carbon is enriched by a factor of the order of 1000, decreasing with increasing temperature. The results suggest that the carbon is in a truly dissolved state and highly mobile, subject to a reversible subsurface segregation. Most probably local lattice strain associated with the solute C species provide the driving force for this diffusional process. The C diffusion coefficient was determined from the (d, p) data below 300 K: D= 10^?13 exp(?7.8/RT) [m²· sec^?1; KJ · mole^?1] and from XPS data between 450–925 K: D = 10^?14 exp(-6/RT) [m² · sec^?1; KJ · mole^?1] The estimated error of the preexponential factors is ± one order of magnitude, that of the activation energies ±3.5 and ±2 KJ mole^?1 respectively.     

DIffusion of Eu and Gd in basalt and obsidian

Tracer diffusion coefficients of ¹⁵³Gd and ¹⁵²Eu in olivine tholeiite have been determined at temperatures between 1150 and 1440°C. The results are identical for both tracers within experimental error. Between 1440 and 1320°C the diffusion coefficients are given by D(Eu, Gd) = 0.058 exp(?40,600/ RT). Between 1320 and 1210°C, the diffusion coefficients are constant at D = (1.4 ± 0.4) × 10^?7 cm²s^?1 and between 1210 and 1150°C, the D values drop irregularly to 4 × 10^?9 cm²s^?1. The liquidus temperature (1270°C) lies within the region of constant D. Such anomalous behavior has not been encountered in previous studies of Ca, Sr, Ba and Co diffusion in basalt. To explain the constant D value near the liquidus, we speculate that the structure of the melt changes as a function of temperature in such a way that the normal temperature dependence of the diffusivity is compensated. For example, the rare earth ions may be displaced from their (high temperature) octahedral coordination sites to other sites where they are more readily dissociated and therefore become progressively more mobile. The behavior below 1210°C may be the result of relatively stable complexes or molecules in the melt or of the formation of a REE bearing crystalline phase that has so far escaped detection. Preliminary results for Eu diffusion in obsidian are D (Eu, 800°C) = 5 × 10^?13 cm² s^?1 and D (Eu, 950°C) = 1.5 × 10^?11 cm² s^?1. These data are consistent with an activation energy of 59 Kcal mole^?1. These low diffusivities indicate that the partitioning of REE in crystallizing intermediate and acidic melts may be controlled by diffusion in the melt rather than equilibrium between the crystal surface and the bulk melt.The diffusion data are applied to partial melting in the mantle, in an attempt to explain how LREE enriched tholeiites may be derived from a LREE depleted mantle source. In this model LREE diffuse from garnet bearing regions that have small melt fractions into garnet free regions that have relatively large melt fractions. REE diffusion is so slow that this process is quantitatively significant only in small partially molten bodies (diameter ～1 km or less) or in larger, but strongly flattened bodies. Internal convective motion during diapiric rise would also increase the efficiency of the process.