Transport properties of olivine grain boundaries from electrical conductivity experiments

Grain boundary processes contribute significantly to electronic and ionic transports in materials within Earth’s interior. We report a novel experimental study of grain boundary conductivity in highly strained olivine aggregates that demonstrates the importance of misorientation angle between adjacent grains on aggregate transport properties. We performed electrical conductivity measurements of melt-free polycrystalline olivine (Fo₉₀) samples that had been previously deformed at 1200 °C and 0.3 GPa to shear strains up to γ?=?7.3. The electrical conductivity and anisotropy were measured at 2.8 GPa over the temperature range 700–1400 °C. We observed that (1) the electrical conductivity of samples with a small grain size (3–6 µm) and strong crystallographic preferred orientation produced by dynamic recrystallization during large-strain shear deformation is a factor of 10 or more larger than that measured on coarse-grained samples, (2) the sample deformed to the highest strain is the most conductive even though it does not have the smallest grain size, and (3) conductivity is up to a factor of ~?4 larger in the direction of shear than normal to the shear plane. Based on these results combined with electrical conductivity data for coarse-grained, polycrystalline olivine and for single crystals, we propose that the electrical conductivity of our fine-grained samples is dominated by grain boundary paths. In addition, the electrical anisotropy results from preferential alignment of higher-conductivity grain boundaries associated with the development of a strong crystallographic preferred orientation of the grains.     

Atomic structures and energies of grain boundaries in Mg2SiO4 forsterite from atomistic modeling

Grain boundaries influence many physical and chemical properties of crystalline materials. Here, we perform molecular dynamics simulations to study the structure of a series of [100] symmetric tilt grain boundaries in Mg₂SiO₄ forsterite. The present results show that grain boundary energies depend significantly on misorientation angle. For small misorientation angles (up to 22°), grain boundary structures consist of an array of partial edge dislocations with Burgers vector $\frac{1}{2}[001]$ associated with stacking faults and their energies can be readily fit with a model which adds the Peach-Koehler equation to the Read-Shockley dislocation model for grain boundaries. The core radius of partial dislocations and the spacing between the partials derived from grain boundary energies show that the transition from low- to high-angle grain boundaries occurs for a misorientation angle between 22° and 32°. For high misorientation angles (32.1° and 60.8°), the cores of dislocations overlap and form repeated structural units. Finally, we use a low energy atomic configuration obtained by molecular dynamics for the misorientation of 12.18° as input to simulate a high-resolution transmission electron microscopy (HRTEM) image. The simulated image is in good agreement with an observed HRTEM image, which indicates the power of the present approach to predict realistic atomic structures of grain boundaries in complex silicates.     

The rheology of olivine and spinel magnesium germanate (Mg2GeO4): TEM study of the defect microstructures

Synthetic polycrystals of α-Mg₂GeO₄ (with the olivine structure) and γ-Mg₂GeO₄ (with the spinel structure) deformed at high temperature and pressure in their respective stability fields were investigated by analytical transmission electron microscopy. Specimens with a mean grain size of 20–30 µm deform by dislocation glide and/or climb. The predominance of glide versus climb depends on stress and grain orientation. The defect microstructures of both polymorphs are very similar to those observed in their respective silicate analogues, α- and γ-(Mg,Fe)₂SiO₄, and, in the case of the spinel phase, very similar to those observed in magnesium aluminate spinels. These observations suggest that Mg₂GeO₄ is a good rheological analogue for the Earth’s upper mantle. A spinel specimen deformed under the same conditions of temperature and strain rate as an olivine specimen was approximately three times stronger than olivine. In specimens of both phases deformed at or above 1400 K, a thin amorphous film composed of Mg, Ge, and O was detected along some grain boundaries. Grains ≤10 µm diameter surrounded by a film of amorphous phase (>10 nm thick) exhibited low dislocation densities, and deformation appeared to have occurred by grain boundary sliding.     

Quantitative characterization of plastic deformation of zircon and geological implications

The deformation-related microstructure of an Indian Ocean zircon hosted in a gabbro deformed at amphibolite grade has been quantified by electron backscatter diffraction. Orientation mapping reveals progressive variations in intragrain crystallographic orientations that accommodate 20° of misorientation in the zircon crystal. These variations are manifested by discrete low-angle (<4°) boundaries that separate domains recording no resolvable orientation variation. The progressive nature of orientation change is documented by crystallographic pole figures which show systematic small circle distributions, and disorientation axes associated with 0.5–4° disorientation angles, which lie parallel to rational low index crystallographic axes. In the most distorted part of the grain (area A), this is the [100] crystal direction. A quaternion analysis of orientation correlations confirms the [100] rotation axis inferred by stereographic inspection, and reveals subtle orientation variations related to the local boundary structure. Microstructural characteristics and orientation data are consistent with the low-angle boundaries having a tilt boundary geometry with dislocation line [100]. This tilt boundary is most likely to have formed by accumulation of edge dislocations associated with a 〈001〉{100} slip system. Analysis of the energy associated with these dislocations suggest they are energetically more favorable than TEM verified 〈010〉{100} slip. Analysis of minor boundaries in area A indicates deformation by either (001) edge, or [100](100) and [001](100) screw dislocations. In other parts of the grain, cross slip on (111), and (112) planes seems likely. These data provide the first detailed microstructural analysis of naturally deformed zircon and indicate ductile crystal-plastic deformation of zircon by the formation and migration of dislocations into low-angle boundaries. Minimum estimates of dislocation density in the low-angle boundaries are of the order of ∼3.10¹⁰ cm⁻². This value is sufficiently high to have a marked effect on the geochemical behavior of zircon, via enhanced bulk diffusion and increased dissolution rates. Therefore, crystal plasticity in zircon may have significant implications for the interpretation of radiometric ages, isotopic discordance and trace element mobility during high-grade metamorphism and melting of the crust.     

Frequency dependent electrical properties of dunite as functions of temperature and oxygen fugacity

Using impedance spectroscopy, we have measured the electrical properties of two dunites and a single crystal olivine sample from 1000 to 1200° C as a function of oxygen fugacity (f _{o 2}). Two conduction mechanisms with resistances that add in series are observed for the dunites corresponding to grain interior and grain boundary conduction mechanisms. The conductivities for each mechanism were determined by analyzing the data using a complex nonlinear least squares fitting routine and the equivalent circuit approach. The grain interiors display a conductivity dependent on f _{o 2} to the 1/5.5–1/7 power, consistent with other determinations, and interpreted as indicating small polaron transport (Fe _Mg ^· ). The grain boundaries demonstrate a weaker f _{o 2} dependence that is dependent on temperature and material. Under certain conditions the f _{o 2} dependence of the grain boundary conductivity is negative. This result indicates that oxygen ion transport is probably not the dominant grain boundary charge transport mechanism; however, an unequivocal determination of the grain boundary mechanism has not been achieved. In some dunites the grain boundaries are more conductive than the grain interiors; in other dunites they are more resistive than grain interiors. The grain boundaries do not enhance the total conductivity of any of the materials of this study but are the controlling mechanism in some instances. Measurement of the complex electrical response at frequencies as low as 10^-4 Hz is required to determine the role of grain boundaries on the overall electrical properties of polycrystalline dunite.     

Fast diffusion along mobile grain boundaries in calcite

Experimental measurements of grain boundary diffusion are usually conducted on static boundaries, despite the fact that grain boundaries deep in the Earth are frequently mobile. In order to explore the possible effect of boundary mobility on grain boundary diffusion rates we have measured the uptake of ⁴⁴Ca from a layer of ⁴⁴Ca-enriched calcite powder during the static recrystallization of a single crystal of calcite at 900°C. A region about 500 μm wide adjacent to the powder layer is heterogeneously enriched in ⁴⁴Ca, and complex zoning patterns, including sharp steps in composition and continuous increases and decreases in ⁴⁴Ca content, are developed. In metamorphic rocks, these would normally be interpreted in terms of changes in pressure or temperature, Rayleigh fractionation, or episodic fluid infiltration. These explanations cannot apply to our experiments, and instead the zoning patterns are interpreted as being due to variations in grain boundary migration rate. We have applied an analytical model which allows the product of grain boundary diffusion coefficient and grain boundary width (D _GB δ) to be calculated from the grain boundary migration rate and the compositional gradient away from the powder layer. The value of D _GB δ in the mobile grain boundaries is at least five orders of magnitude greater than the published value for static boundaries under the same conditions. In order to allow the scale of chemical equilibrium (and hence textural evolution) to be predicted under both experimental and geological conditions, we present quantitative diffusion-regime maps for static and mobile boundaries in calcite, using both published values and our new values for grain boundary diffusion in mobile boundaries. Enhanced diffusion in mobile boundaries has wide implications for the high temperature rheology of Earth materials, for geochronology, and for interpretations of the length- and time-scales of chemical mass-transport. Moreover, zones of anomalously high electrical conductivity in the crust and mantle could be regions undergoing recrystallization such as active shear zones, rather than regions of anomalous mineralogy, water- or melt-content as is generally suggested.     

Hydrogen diffusion in spinel grain boundaries and consequences for chemical homogenization in hydrous peridotite

Hydrogen can be stored in the structure of nominally anhydrous minerals as point defects, and these impurities substantially modify many physical properties of Earth’s mantle minerals. However, mantle rocks are composed of mineral grains separated by grain boundaries and interphase grains boundaries. Therefore, as a potential hydrogen reservoir, grain boundaries should be given proper attention. Here, I report an experimental investigation into hydrogen diffusion through grain boundaries in polycrystalline aggregates. Sintering and diffusion experiments were performed using a gas-medium high-pressure vessel at under pressure of 300 MPa and over a temperature range of 900–1,250°C. The diffusion assembly consisted of a polycrystalline cylinder of aluminous spinel + olivine crystals with a talc cylinder as the main hydrogen source. A Ni capsule was used to buffer the oxygen fugacity at Ni–NiO. Experimental durations varied from 5 min to 5 h. The presence of hydrogen in the crystals was measured by Fourier-transform infrared spectroscopy. The calculation of the diffusion coefficients was based on the estimation of the characteristic distance. The absence or presence of hydrogen recorded by the ‘hydrogen sensor’ olivines embedded in the aggregate allows the estimation of bounds on this characteristic distance. Results presented here suggest that hydrogen effective diffusion coefficients are only one order of magnitude faster (~10⁻⁹ m²s⁻¹ at 1,000°C) than in an olivine single crystal along the [100] axis. Resulting diffusion coefficients for hydrogen in grain boundary are four orders of magnitude faster than in a single crystal, but this diffusivity is not fast enough to affect hydrogen mobility in mantle rocks with grain sizes greater than ~1 mm. Thus, very limited chemical homogenization would occur using grain boundaries diffusion in mantle hydrous peridotite for incompatible and volatile element, such as hydrogen.     

Grain growth in porous olivine aggregates

Grain growth experiments have been performed at 1 atm on fine grain size (<10 μm) synthetic olivine (Fo₉₁) aggregates at various temperatures (1200° to 1400° C), oxygen fugacities (10^-4 to 10^-11 atm) and total anneal times (10, 30, 60, 100 and 200 h). The rate of grain growth increased with increasing temperature and with increasing oxygen fugacity. The presence of a second phase (residual porosity), introduced during sample fabrication, has a significant effect on grain growth, with evolution in grain size paralleled by changes in the size and frequency of the pores. When the grain growth data were fit to a growth law G ⁿ ?G _O ⁿ =κ₀ tf ₀ ^m ₂e^?Q/RT, the growth exponents fall in the range of n=4 to 5, suggesting that grain growth may be controlled by the coalescence of the second phase. The evolution in pore size and frequency may occur either by the transport of the ionic species constituting olivine between the pores or by the movement of the pores themselves along the grain boundaries and edges. Thus, the rate of growth of the pores and grains is probably limited by diffusion of the slowest ionic species constituting olivine (magnesium, iron, silicon, or oxygen) moving along the fastest path for that species (through the lattice, along the grain boundaries, around the surface of the moving pores, or through the vapor phase in the pores). Activation energies for grain growth of Q=290 ± 20 kJ/mol and 345 ± 25 kJ/mol were calculated from our results for n=4 and 5, respectively. These activation energies preclude vapor-phase transport and iron diffusion along grain boundaries but do not otherwise permit a discrimination between the rate limiting species or path. The oxygen fugacity exponent of m ≈0.12 suggests that lattice diffusion does not control the grain growth. However, the lack of data for magnesium, iron, silicon and oxygen surface and grain boundary diffusion in olivine makes definitive determination of the mechanism controlling grain growth difficult.     

The legacy of crystal-plastic deformation in olivine: high-diffusivity pathways during serpentinization

Crystal-plastic olivine deformation to produce subgrain boundaries composed of edge dislocations is an inevitable consequence of asthenospheric mantle flow. Although crystal-plastic deformation and serpentinization are spatio-temporally decoupled, we identified compositional readjustments expressed on the micrometric level as a striped Fe-enriched ( [`(X)]<sub>\textFe</sub> \bar{X}_{\text{Fe}}  = 0.24 ± 0.02 (zones); 0.12 ± 0.02 (bulk)) or Fe-depleted ( [`(X)]_\textFe \bar{X}_{\text{Fe}}  = 0.10 ± 0.01 (zones); 0.13 ± 0.01 (bulk)) zoning in partly serpentinized olivine grains from two upper mantle sections in Norway. Focused ion beam sample preparation combined with transmission electron microscopy (TEM) and aberration-corrected scanning TEM, enabling atomic-level resolved electron energy-loss spectroscopic line profiling, reveals that every zone is immediately associated with a subgrain boundary. We infer that the zonings are a result of the environmental Fe²⁺Mg₋₁ exchange potential during antigorite serpentinization of olivine and the drive toward element exchange equilibrium. This is facilitated by enhanced solid-state diffusion along subgrain boundaries in a system, which otherwise re-equilibrates via dissolution-reprecipitation. Fe enrichment or depletion is controlled by the silica activity imposed on the system by the local olivine/orthopyroxene mass ratio, temperature and the effect of magnetite stability. The Fe-Mg exchange coefficients K_\textD^{\textAtg/\textOl} K_{\text{D}}^{{{\text{Atg}}/{\text{Ol}}}} between both types of zoning and antigorite display coalescence toward exchange equilibrium. With both types of zoning, Mn is enriched and Ni depleted compared with the unaffected bulk composition. Nanometer-sized, heterogeneously distributed antigorite precipitates along olivine subgrain boundaries suggest that water was able to ingress along them. Crystallographic orientation relationships gained via electron backscatter diffraction between olivine grain domains and different serpentine vein generations support the hypothesis that serpentinization was initiated along olivine subgrain boundaries.     

Thin amorphous films (1–2 nm) at olivine grain boundaries in mantle xenoliths from San Carlos,Arizona

 Olivine grain boundaries and phase boundaries in xenoliths from San Carlos have been investigated by high-resolution transmission electron microscopy (HREM) and analytical electron microscopy (AEM). Thin amorphous intergranular layers with variable width (1–2 nm) were detected along olivine grain boundaries. The Al₂O₃, TiO₂ and CaO concentrations of the amorphous layers increase with increasing width of the layer. The composition of the amorphous intergranular layers depends on the interface type – grain or phase boundary. Morphology, amorphous state and chemical composition of the intergranular layer suggest the presence of a melt film at olivine grain boundaries. Since the composition of the amorphous phase strongly depends on the type of interface, the melt must have been generated at the grain boundary. Also, the melt chemistry is different from the composition of partial melts produced from possible hydrous phases, such as phlogopite or amphibole, and from the host basanite. The mobility of very thin melt films is assumed to be very limited due to the strong interface forces between the melt and the grain boundary. It is concluded that grain boundary melting occurred at the interfaces due to decompression during uplift. The melt wetted olivine grain boundaries as well as olivine-opx phase boundaries. The thin amorphous layers formed melt microsystems. Mixing of melts from different microsystems is suggested to occur in wider melt films, melt veins or melt pockets thus creating a magmatic melt that could be extracted from its source. Received: 6 November 1995 / Accepted: 24 January 1996     

Flow behavior and microstructures of hydrous olivine aggregates at upper mantle pressures and temperatures

Deformation experiments on olivine aggregates were performed under hydrous conditions using a deformation-DIA apparatus combined with synchrotron in situ X-ray observations at pressures of 1.5–9.8 GPa, temperatures of 1223–1800 K, and strain rates ranging from 0.8 × 10^?5 to 7.5 × 10^?5 s^?1. The pressure and strain rate dependencies of the plasticity of hydrous olivine may be described by an activation volume of 17 ± 6 cm³ mol^?1 and a stress exponent of 3.2 ± 0.6 at temperatures of 1323–1423 K. A comparison between previous data sets and our results at a normalized temperature and a strain rate showed that the creep strength of hydrous olivine deformed at 1323–1423 K is much weaker than that for the dislocation creep of water-saturated olivine and is similar to that for diffusional creep and dislocation-accommodated grain boundary sliding, while dislocation microstructures showing the [001] slip or the [001](100) slip system were developed. At temperatures of 1633–1800 K, a much stronger pressure effect on creep strength was observed for olivine with an activation volume of 27 ± 7 cm³ mol^?1 assuming a stress exponent of 3.5, water fugacity exponent of 1.2, and activation energy of 520 kJ mol^?1 (i.e., power-law dislocation creep of hydrous olivine). Because of the weak pressure dependence of the rheology of hydrous olivine at lower temperatures, water weakening of olivine could be effective in the deeper and colder part of Earth’s upper mantle.     

Grain boundary diffusion in enstatite

We have investigated grain boundary diffusion rates in enstatite by heating single crystals of quartz packed in powdered San Carlos olivine (Mg_0.90Fe_0.10)₂SiO₄ at controlled oxygen fugacities in the range 10^?5.7 to 10^?8.7?atm and temperatures from 1350° to 1450?°C for times from 5 to 100?h at 1?atm total pressure. Following the experiments, the thickness of the coherent polycrystalline reaction rim of pyroxene that had formed between the quartz and olivine was measured using backscatter scanning imaging in the electron microprobe. Quantitative microprobe analysis indicated that the composition of this reaction phase is (Mg_0.92Fe_0.08)₂Si₂O₆. The rate of growth of the pyroxene increases with increasing temperature, is independent of the oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by ionic diffusion through the pyroxene rim. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the olivine-pyroxene interface, and is therefore controlled by the diffusion of silicon and oxygen. The parabolic rate constants determined from the experiments were analyzed in terms of the oxide activity gradient across the rim to yield mean effective diffusivities for the rate-limiting ionic species, assuming bulk transport through the pyroxene layer. These effective diffusivities are faster than the lattice diffusivities for the slowest species (silicon) calculated from creep experiments, but slower than measured lattice diffusivities for oxygen in enstatite. Thus, silicon grain boundary diffusion is most likely to be the rate-limiting process in the growth of the pyroxene rims. Also, as oxygen transport through the pyroxene rims must be faster than silicon transport, diffusion of oxygen along the grain boundaries must be faster than through the lattice. The grain boundary diffusivity for silicon in orthopyroxenite is then given by D¯gbSiδ=(3.3±3.0)×10^?9f0.0O₂e^?400±65/RT?m³s^?1, where the activation energy for diffusion is in kJ/mol, and δ is the grain boundary width in m. Calculated growth rates for enstatite under these conditions are significantly slower than predicted by an extrapolation from similar experiments performed at 1000?°C under high pressure (hydrous) conditions by Yund and Tullis (1992), perhaps due to water-enhancement of diffusion in their experiments.     

Thin intergranular melt films and melt pockets in spinel peridotite xenoliths from the Rhön area (Germany): early stage of melt generation by grain boundary melting

Optical microscopy and transmission electron microscopy (TEM) on a porphyroclastic high temperature spinel peridotite from the Rhön area reveal fine, irregular glass layers and pockets along mineral interfaces, cracks in olivine, inside olivine crystals and in spongy rims of clinopyroxene. The chemical composition of the glass deviates significantly from the composition of the host basanite. Electron diffraction technique confirms the amorphous nature of the glass, thus classifying it as a former melt. Every grain or phase boundary shows amorphous intergranular glass layers of variable thickness and characteristic chemical composition with distinct chemical inhomogeneities. Olivine grain boundaries, as the most common type of interfaces, exhibit two different types of melt glasses: (1) Type I melt at olivine grain boundaries, which is characterized by low contents of SiO₂ (～37?wt%) and Al₂O₃ (～5?wt%) and elevated contents of MgO (～31?wt%) and FeO (～22?wt%), is supposed to have formed prior to or during the thermal overprint and the dynamic recrystallisation of the xenolith in the mantle. Melt inclusions inside olivine grains with an average composition of type I melt are suggested to be earlier melt droplets at olivine interfaces, overgrown by migrating olivine grain boundaries during recrystallization in the mantle prior to the uplift of the xenolith. (2) Type II melt, the most common type of melt in the xenolith, shows higher contents of SiO₂ (～48?wt%) and Al₂O₃ (～17?wt%) but lower contents of MgO (～20?wt%) and FeO (～11?wt%). The observation of different types of glass within a single xenolith indicates the development of different chemical melt equilibria at interfaces or triple junctions in the xenolith. The absence of geochemical trends in bivariate plots excludes a unifying process for the genesis of these glasses. Melt inclusions in the spongy rims of clinopyroxene are interpreted to be the product of a potassium-rich metasomatism. The formation of most amorphous intergranular melt layers and pockets at the mineral interfaces including type II melt at olivine grain boundaries is suggested to result from decompression melting during the uplift with the basalt magma. We suggest that these glasses were produced by grain boundary melting due to lattice mismatch and impurity segregation. The observed intergranular amorphous layers or melts represent the very beginning of mineral melting by grain boundary melting.     

Synthetic near Σ5 (210)/[100] grain boundary in YAG fabricated by direct bonding: structure and stability

Several macroscopic physical and chemical properties, such as rheology, elasticity, or transport properties are governed by grain boundary processes. An improved understanding of the structure and evolution of grain boundaries has thus become a key challenge in geosciences and material sciences. Here, we report the structure of near Σ5 (210)/[100] grain boundaries in Y₃Al₅O₁₂ (YAG), which were synthesised by the wafer direct bonding method. The produced grain boundaries were annealed at different temperatures, ranging from 673 to 1,873 K. The grain boundaries annealed at different temperatures are not distinguishable based on their flatness and apparent cohesiveness in high resolution TEM (HRTEM) micrographs, but show a considerable step in their mechanical stability at around 1,273 K, a temperature that corresponds to roughly half the melting temperature of YAG. This study further focuses on the effect of a slight misorientation of the two crystals on the grain boundary structure and we discuss if the boundary can reach a state of minimum energy configuration during annealing. Along the grain boundaries, we observed a long-range strain contrast with a periodicity of 40 nm, which has not been reported for high-angle grain boundaries so far. We conclude that this contrast is caused by faceting along the grain boundary plane, which is needed to achieve minimum energy configuration of the grain boundary plane.     

Kinetics of diffusion-controlled growth of fayalite

The rate of growth of fayalite (Fe₂SiO₄) has been measured at one atmosphere total pressure, temperatures from 1000° to 1120° C, and oxygen fugacities controlled by CO/CO₂ gas-mixing from 10^-9.9 to 10^-13.0atm, chosen to span the fayalite stability field. The fine-grained polycrystalline fayalite layer was formed by reacting the oxides FeO or Fe₃O₄ with a thin slice of single-crystal quartz. The rate of growth of the fayalite increases with increasing temperature and decreasing oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by diffusion through the fayalite. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the quartz-fayalite interface, and is therefore controlled by the diffusion of iron and oxygen. The parabolic rate constant was analyzed in terms of the oxide activity gradient to yield mean chemical diffusivities for the rate-limiting ionic species, assuming bulk transport through the fayalite layer. Given that iron diffusion in olivine polycrystals occurs either by lattice diffusion, which shows a positive dependence on oxygen activity, or by grain boundary diffusion, which would result in growth rates significantly faster than we observe, we conclude that the diffusivities derived in this study represent oxygen diffusion. However, since oxygen lattice diffusion in fayalite has been established to be much slower than our measurements, it is likely that the transport path for oxygen is along the grain boundaries. Thus, the mean grain boundary diffusivity of oxygen in fayalite $\bar D$ _O ^gb (m² s^-1), using the measured grain size of 0.25 μm, is then given by $$\bar D_O^{gb} {\mathbf{ }}\delta = 1.28 \times 10^{ - 3} f_{O_2 }^{ - 0.17} {\mathbf{ }}e^{ - 540/RT} $$ , where δ is the grain boundary width (in m), and the activation energy is in kJ/mol. Assuming δ=10^-9 m (Ricoult and Kohlstedt 1983), the oxygen grain boundary diffusivities are about a factor of 30 × slower than those reported by Watson (1986) for Fo₉₀ olivine.     

Application of the cBΩ model to the calculation of diffusion parameters of He in olivine

The validity of the thermodynamic cBΩ model is tested in terms of the experimentally determined diffusion coefficients of He in a natural Fe-bearing olivine (Fo₉₀) and a synthetic end-member forsterite (Mg₂SiO₄) over a broad temperature range (250–950 °C), as reported recently by Cherniak and Watson (Geochem Cosmochim Acta 84:269–279, 2012). The calculated activation enthalpies for each of the three crystallographic axes were found to be (134 ± 5), (137 ± 13) and (158 ± 4) kJ mol^?1 for the [100], [010] and [001] directions in forsterite, and (141 ± 9) kJ mol^?1 for the [010] direction in olivine, exhibiting a deviation of <1 % with the corresponding reported experimental values. Additional point defect parameters such as activation volume, activation entropy and activation Gibbs free energy were calculated as a function of temperature. The estimated activation volumes (3.2–3.9 ± 0.3 cm³ mol^?1) of He diffusion in olivine are comparable with other reported results for hydrogen and tracer diffusion of Mg cations in olivine. The pressure dependence of He diffusion coefficients was also determined, based on single experimental diffusion measurements at 2.6 and 2.7 GPa along the [001] direction in forsterite at 400 and 650 °C.     

Grain boundary dissolution porosity in quartzofeldspathic ultramylonites: Implications for permeability enhancement and weakening of mid-crustal shear zones

Quartzofeldspathic ultramylonites from the Alpine Fault Zone, one of the world's major, active plate boundary-scale fault zones have quartz crystallographic preferred orientations (CPO) and abundant low-angle (<10° misorientation) boundaries, typical microstructures for dislocation creep-dominated deformation. Geometrically necessary dislocation density estimates indicate mean dislocation densities of ∼10⁹ cm⁻². A significant proportion (∼30%) of grain boundaries (>10° misorientation) are decorated by faceted pores, commonly with uniformly-oriented pyramidal shapes. Only grain boundaries with >10° misorientation angles in polymineralic aggregates are decorated by pores. Mean grain boundary pore densities are ∼5 × 10⁸ cm⁻². Grain boundary pores are dissolution pits generated during syn-deformational transient grain boundary permeability, nucleating on dislocation traces at dilatant grain boundary interfaces. They have not been removed by subsequent grain boundary closure or annealing. Pore decoration could have led to grain boundary pinning, triggering a switch in the dominant deformation mechanism to grain boundary sliding, which is supported by evidence of CPO destruction in matrix quartz. Pore-decorated grain boundaries have significantly reduced surface area available for adhesion and cohesion, which would reduce the tensile and shear strength of grain boundaries, and hence, the bulk rock. Grain boundary decoration also significantly decreased the mean distance between pores, potentially facilitating dynamic permeability. Consequently, these microstructures provide a new explanation for strain weakening and evidence of fluid flow along grain boundaries in mylonites at mid-crustal conditions.     

Diffusion along interphase boundaries and its effect on retrograde zoning patterns of metamorphic minerals

In this study we use two dimensional chemical patterns and numerical modeling to estimate the relative rates of chemical transport along interphase boundaries (ib) and through grain (s) interiors during retrograde Fe–Mg exchange between garnet and biotite at a garnet–biotite–quartz triple junction. We demonstrate that systematic variations in garnet–rim compositions, which are frequently observed along garnet–quartz interfaces, and deviations from concentric retrograde zoning patterns start to develop when chemical transport along the interphase boundaries becomes slow during cooling. The capacities for chemical transport along an interphase boundary depend on the product D ^ib K ^ib/s a, where D ^ib is the diffusion coefficient of the exchangeable components within the interphase boundary medium, K ^ib/s is the equilibrium partitioning coefficient between the cation exchange partners and the interphase boundary medium and a is the interphase boundary width. The model is applied to analyze the retrograde zoning patterns in garnets from the Mozambique belt (SE-Kenya), which cooled from 820°C at a rate of ca. 2°C/my. It is found that non-equilibrated compositions in garnet along garnet/quartz interphase boundaries started to develop below 700°C due to insufficient rates of chemical transport along these boundaries. The transport capacities of garnet/quartz interphase boundaries was estimated to have been between about 1E-23 cm³/s (575°C) and 1E-20 cm³/s (700°C) from modeling the observed X _Fe pattern in garnet close to a garnet–quartz–biotite triple junction and relying on published data on the diffusivity of Fe²⁺ in garnet. Similar transport capacities are obtained; when the interphase boundary is assumed to be filled with a material that has the transport properties and chemical composition of a free melt in equilibrium with garnet, biotite and quartz at the respective conditions. In contrast, if the transport properties of the interphase boundary medium are related to the diffusivities and solubility of Fe²⁺ and FeOH⁺ within a free aqueous solution, chemical transport along the interphase boundaries would be much more efficient, and exchange equilibrium would have been maintained during the entire cooling history of the rocks. The observation of systematic deviations from local equilibrium along the garnet–quartz interphase boundaries leads us to exclude the presence of an aqueous fluid along the interphase boundary at any time during cooling.     

Sintering of olivine and olivine-basalt aggregates

The sintering behavior of olivine and olivine-basalt aggregates has been examined at temperatures near 1,300° C. Experimental factors contributing to rapid sintering kinetics and high-density, fine-grained specimens include: (i) the uniform dispersion of basalt throughout the specimen, (ii) a very fine, uniform particle size for the olivine powder, (iii) oxygen fugacities near the high P _O2 end of the olivine stability field, and (iv) rapid heating to the sintering temperature. Olivine-basalt specimens prepared from olivine particles coated with a synthetic basalt achieve chemical and microstructural equilibrium more rapidly, as well as produce higher density and finer grain-sized aggregates, than do specimens prepared by mechanical mixing of olivine and natural basalt powders. The grain boundary mobility for olivine, measured for olivine-basalt aggregates which have undergone secondary recrystallization, is on the order of 2×10^?15 (m/s)/(N/m²) in the temperature range 1,300–1,400° C. Solution-precipitation (pressure-solution) processes make an important contribution to the development of the microstructure in olivine-basalt aggregates.     

Olivine [100] normal to foliation: lattice preferred orientation in prograde garnet peridotite formed at high H<Subscript>2</Subscript>O activity,Cima di Gagnone (Central Alps)

Automated electron backscattered diffraction (EBSD) was applied using a scanning electron microscope to obtain lattice preferred orientation (LPO) data for olivine in garnet peridotites of the Central Alps. As a reference frame, the LPOs of enstatite were also investigated. In the garnet peridotite at Cima di Gagnone (CDG), a weak foliation carrying a distinct lineation is present. The lineation is characterized by elongated enstatite, olivine and poikiloblastic garnet. Olivine shows a very unusual LPO with [100] normal to foliation and [001] parallel to lineation. Achsenverteilungsanalyse (AVA) maps demonstrate that [001] of olivine grains corresponds quite well to their maximum length axes which are preferentially parallel to the lineation. Numerous planar hydrous defects within (001) planes of olivine are marked by palisades of ilmenite rods and show a preferred orientation normal to lineation. Calculated P-wave velocities for CDG are fastest (8.32 km s^у) normal to foliation with a relatively low anisotropy (2.9%). Compared to mantle peridotites with the usual (010)[100] LPO where the fastest Vp direction is towards the lineation, the relationship between flow geometry and seismic anisotropy is significantly different at CDG. Several mechanisms for the formation of the LPO type at CDG are considered, with glide possible on (100)[001] of olivine. On the basis of field data as well as petrographic and petrologic evidence, it has been demonstrated that the CDG garnet peridotite formed by prograde metamorphism from a hydrous protolith at pressures and temperatures of about 3.0 GPa and 750 °C, respectively. The CDG LPO is interpreted to have formed during hydrous subduction zone metamorphism. The same interpretation may hold for the previously investigated olivine LPO at Alpe Arami, which is similar to that at the nearby CDG. The observed anomalous LPO is no proof for ultradeep (>3.0 GPa) conditions.