Clinopyroxene dissolution in basaltic melt

The history of magmatic systems may be inferred from reactions between mantle xenoliths and host basalt if the thermodynamics and kinetics of the reactions are quantified. To study diffusive and convective clinopyroxene dissolution in silicate melts, diffusive clinopyroxene dissolution experiments were conducted at 0.47–1.90 GPa and 1509–1790 K in a piston-cylinder apparatus. Clinopyroxene saturation is found to be roughly determined by MgO and CaO content. The effective binary diffusivities, D_MgO and D_CaO, and the interface melt saturation condition, , are extracted from the experiments. D_MgO and D_CaO show Arrhenian dependence on temperature. The pressure dependence is small and not resolved within 0.47–1.90 GPa. in the interface melt increases with increasing temperature, but decreases with increasing pressure. Convective clinopyroxene dissolution, where the convection is driven by the density difference between the crystal and melt, is modeled using the diffusivities and interface melt saturation condition. Previous studies showed that the convective dissolution rate depends on the thermodynamics, kinetics and fluid dynamics of the system. Comparing our results for clinopyroxene dissolution to results from a previous study on convective olivine dissolution shows that the kinetic and fluid dynamic aspects of the two minerals are quite similar. However, the thermodynamics of clinopyroxene dissolution depends more strongly on the degree of superheating and composition of the host melt than that of olivine dissolution. The models for clinopyroxene and olivine dissolution are tested against literature experiments on mineral–melt interaction. They are then applied to previously proposed reactions between Hawaii basalts and mantle minerals, mid-ocean ridge basalts and mantle minerals, and xenoliths digestion in a basalt at Kuandian, Northeast China.     

Diffusive dissolution in ternary systems: analysis with applications to quartz and quartzite dissolution in molten silicates

Exact solutions to equations governing isothermal diffusive dissolution of a crystalline slab in a ternary liquid were obtained to include the effect of coupled chemical diffusion in the liquid. These analytical results, supplemented by approximate solutions valid for slow dissolving, provide new insights into the characteristics of diffusive dissolution in ternary systems. Dissolution rate is proportional to square root of time in diffusive dissolution. The coefficient of proportionality is a function of diffusion coefficients, liquidus relation, melt composition at the crystal–melt interface, and compositions of the dissolving crystal and starting melt. In the limit of slow dissolving, the dissolution rate can be written in terms of three dimensionless parameters that are functions of the aforementioned parameters. Dissolution rate is proportional to the diffusion rate of the slow eigen component in the melt when the diffusion rate of the minor eigen component is much slower than the diffusion rate of the major eigen component.Laboratory experiments of diffusive dissolution of single crystals and polycrystalline aggregates of quartz in a haplodacitic melt (25 wt.% CaO, 15 wt.% Al₂O₃, and 60 wt.% SiO₂) were conducted at 1500°C and 1 GPa. Measured dissolution distances (X_b, in microns) are proportional to the square root of experimental run time (t, in seconds), X_b = −0.620 (±0.019) √t. Chemical concentration profiles measured from quenched melts are invariant with time when displayed against the distance (measured from the crystal–melt interface) normalized by the square root of time. The melt compositions at the crystal–melt interface, extrapolated from the measured diffusion profiles in the quenched melts, are within 0.2 wt.% of the independently measured quartz liquidus in the ternary CaO–Al₂O₃–SiO₂ at 1500°C and 1 GPa. These results suggest that crystal and melt are in chemical equilibrium at their interface shortly after the onset of dissolution. Diffusive dissolution of quartz and quartzite is characterized by slow dissolving. Using quartz liquidus as one of the boundary conditions, it has been shown that the calculated dissolution distances and concentration profiles are in good agreement with the experimentally measured ones. Coupled diffusion played an essential role in quartz and quartzite dissolution in haplodacitic to haplobasaltic melts, and is likely to play an important role in diffusion-limited kinetic processes such as crystal growth and dissolution in natural melts of basaltic–rhyolitic compositions.     

Diffusive crystal dissolution

Crystal dissolution may include three component processes: interface reaction, diffusion and complications due to convection. We report here a theoretical and experimental study of crystal dissolution in silicate melt without convection. A reaction-diffusion equation is developed and numerically solved. The results show that during non-convective crystal dissolution in silicate melt, the interface melt composition reaches a constant or stationary saturation composition in less than a second, hence interface reaction is not the rate-determining step and crystal dissolution in silicate melt is usually diffusion-controlled. Crystal dissolution experiments (designed to suppress convection) show that the concentration profiles of all components propagate into the melt according to the square root of run duration, and that the dissolution distance is also proportional to the square root of run duration. Thus our experiments confirm that the dissolution is diffusion controlled, which is consistent with our numerical calculations. For some principal equilibrium-determining components, concentration profiles conform approximately to the analytical solution of the diffusion equation with a constant effective binary diffusion coefficient. Diffusive dissolution rates (which are inversely proportional to square root of time) can thus be predicted from the phase equilibria and the effective binary diffusion coefficients. To predict steady-state convective dissolution rates, the thickness of the boundary layer must be known. If the convective compositional boundary layer thickness around a dissolving crystal aggregate or near the wall of a magma chamber during convection is about 2 cm or larger, then convective dissolution would rarely result in any significant alteration of original melt. Our dissolution experiments also illustrate the complexity of the diffusion process. Uphill diffusion is common, especially during olivine dissolution into andesitic melt where a majority of the components show the effect of diffusion up their own concentration gradients. Uphill diffusion has implications to the understanding of crystal zoning, and suggests caution is required in applying least squares mass balance analysis to magmatic rocks affected by processes involving diffusion.     

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.     

Mg-Fe partitioning between olivine and ultramafic melts at high pressures

Precise determination of the partitioning of Mg and Fe²⁺ between olivine and ultramafic melt has been made at pressures from 5 to 13 GPa using a MA-8 type multi-anvil high-pressure apparatus (PREM) installed at Earthquake Research Institute, University of Tokyo. A very short rhenium capsule (<100 μm sample thickness) was adopted to minimize temperature variation within the sample container. Synthetic gels with the composition of the upper mantle peridotite were used as starting materials to promote the homogeneity. Analyses of quenched melts and coexisting olivines were made with an electron probe microanalyzer. The obtained partition coefficient, K_D [=(FeO/MgO)^ol/(FeO/MgO)^melt], decreases from 0.35 to 0.25 with increasing pressure from 5 to 13 GPa, suggesting a negative correlation between pressure and K_D above 5 GPa. Our result is consistent with a parabolic relationship between K_D and degree of polymerization (NBO/T) of melts reported by previous studies at lower pressures. The negative correlation between pressure and K_D suggests that olivine crystallizing in a magma ocean becomes more Mg-rich with depth and that primary magmas generated in the upper mantle become more Fe-rich with depth than previously estimated.     

Temperature dependence of sulfide and sulfate solubility in olivine-saturated basaltic magmas

The sulfur concentration at pyrrhotite- and anhydrite-saturation in primitive hydrous basaltic melt of the 2001-2002 eruption of Mt. Etna was determined at 200 MPa, T = 1050-1250 °C and at log fO₂ from FMQ to FMQ+2.2 (FMQ is Fayalite-Magnetite-Quartz oxygen buffer). At 1050 °C Au sample containers were used. A double-capsule technique, using a single crystal olivine sample container closed with an olivine piston, embedded in a sealed Au₈₀Pd₂₀ capsule, was developed to perform experiments in S-bearing hydrous basaltic systems at T > 1050 °C. Pyrrhotite is found to be a stable phase coexisting with melt at FMQ-FMQ+0.3, whereas anhydrite is stable at FMQ+1.4-FMQ+2.2. The S concentration in the melt increases almost linearly from 0.12 ± 0.01 to 0.39 ± 0.02 wt.% S at FeS-saturation and from 0.74 ± 0.01 to 1.08 ± 0.04 wt.% S at anhydrite-saturation with T ranging from 1050-1250 °C. The relationships between S concentration at pyrrhotite and/or anhydrite saturation, MgO content of the olivine-saturated melt, T, and log fO₂ observed in this study and from previous data are used to develop an empirical model for estimating the magmatic T and fO₂ from the S and MgO concentrations of H₂O-bearing olivine-saturated basaltic melts. The model can also be used to determine maximum S concentrations, if fO₂ and MgO content of the melt are known. The application of the model to compositions of melt inclusions in olivines from Mt. Etna indicates that the most primitive magmas trapped in inclusions might have been stored at log fO₂ slightly higher than FMQ+1 and at T = 1100-1150 °C, whereas more evolved melts could have been trapped at T ? 1100 °C. These values are in a good agreement with the estimates obtained by other independent methods reported in the literature.     

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)     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

The effect of experiment geometry on the mechanism and rate of dissolution of quartz in basanite at 0.5 GPa and 1350 °C

Mineral dissolution is an important factor in many magmatic processes such as melting, assimilation and magma mixing. Since it is not possible to determine dissolution rates or mechanisms from natural samples, experimental measurements are very useful. However, the geometry of the crystal–melt system can have a large effect on the measured rate, depending on whether the contaminated melt formed during dissolution is gravitationally stable or unstable. This study examines the effects of the crystal–melt geometry on the dissolution rate and mechanism. The experiments were performed using basanite melt and cylinders and spheres prepared from a single crystal of natural quartz. All of the experiments were performed in the piston cylinder apparatus at 0.5 GPa and 1350 °C. Four crystal–melt geometries were used: (1) quartz cylinders on top of a column of melt; (2) quartz cylinders beneath a column of basanite melt; (3) quartz cylinders in the middle of column of melt; (4) quartz spheres on top of a column of basanite melt. These geometries allow an examination of non-convective, convective and mixed non-convective/convective dissolution. Sphere experiments were included, as this has been the most commonly used geometry in previous experimental studies. In all of the experiments quartz dissolves directly into the basanite without formation of cristobalite or tridymite. Quartz on top of a column of melt dissolves at a rate almost proportional to the square root of time and forms a silica-rich compositional boundary layer that is gravitationally stable. All of the samples show well-defined compositional gradients in the boundary layer; however, the melt at the interface varies in composition with time and plots of concentration as a function of distance normalized to time show that the diffusion rate of SiO₂ increases with time. These data suggest that the rate-controlling step during quartz dissolution is interface reaction rather than cation diffusion. Quartz on the bottom of a column of basanite dissolves much more quickly than in the quartz-on-top experiments and the dissolution rate is linear, due to the periodic gravitational instability and resultant convection of the boundary layer. Even though interface kinetics are the rate-controlling step in quartz dissolution, convection causes an increase in dissolution rate because it replenishes the boundary layer with new, silica-undersaturated melt, which dissolves the quartz more quickly than the contaminated melt. These data suggest that the interface reaction rate is controlled by the degree of undersaturation of the solvent melt in the dissolving component. Both quartz-in-middle and quartz sphere experiments dissolve at a rate intermediate between the two extremes and both show a power law rate. Both dissolve by a combination of convective and non-convective dissolution but the sphere experiments are affected by an additional factor. During the experiment the sphere can sink through the capsule causing forced convection which adds another complication to the interpretation of the dissolution rate data. The results of this study indicate that the choice of experiment geometry plays a major role in determining the observed dissolution rate. Mineral spheres, which have been widely used in the past, are not ideal for dissolution studies. Instead, dissolution rates and mechanisms are best determined in the absence of convection. These experiments have an additional advantage in that for diffusion-controlled dissolution, they allow determination of cation diffusivity. Received: 2 March 2000 / Accepted: 11 April 2000     

An experimental study of anorthosite dissolution in lunar picritic magmas: Implications for crustal assimilation processes

Experiments characterizing the kinetics of anorthosite dissolution in lunar picritic magmas (very low-Ti, low-Ti, and high-Ti picritic glasses) were conducted at 0.6 GPa and 1250-1400 °C using the dissolution couple method. Reaction between the anorthosite and lunar picritic magmas at 1250-1300 °C produced a spinel + melt layer. Reaction between the anorthosite and an olivine-saturated low-Ti magma at 1250-1300 °C produced a crystal-free region between the spinel + melt layer and the olivine-saturated magma. The anorthosite dissolution experiments conducted at 1400 °C simply dissolved anorthosite and did not result in a crystal-bearing region. The rate of anorthosite dissolution strongly depends on temperature and composition of the reacting melt. Concentration profiles that develop during anorthosite dissolution are nonlinear and extend from the picritic glass compositions to anorthite. These profiles feature a large and continuous variation in melt density and viscosity from the anorthosite-melt interface to the initial picritic magmas. In both the low-Ti and high-Ti magmas the diffusive fluxes of TiO₂, Al₂O₃, and SiO₂ are strongly coupled to the concentration gradients of CaO and FeO. Anorthosite dissolution may play an important role in producing the chemical variability of the lunar picritic magmas, the origin of spinel in the lunar basalts and picritic glasses, and the petrogenesis of the high-Al basalts.     

Contrasting types of metasomatism in dunite, wehrlite and websterite xenoliths from Kimberley, South Africa

Dunite, wehrlite and websterite are rare members of the mantle xenolith suite in the Kimberley kimberlites of the Kaapvaal Craton in southern Africa. All three types were originally residues of extensive melt extraction and experienced varying amounts and types of melt re-enrichment. The melt depletion event, dated by Re-Os isotope systematics at 2.9 Ga or older, is evidenced by the high Mg# (Mg/(Mg + Fe)) of silicate minerals (olivine (0.89-0.93); pyroxene (0.88-0.93); garnet (0.72-0.85)), high Cr# (Cr/(Cr + Al)) of spinel (0.53-0.84) and mostly low whole-rock SiO₂, CaO and Al₂O₃ contents. Shortly after melt depletion, websterites were formed by reaction between depleted peridotites and silica-rich melt (>60 wt% SiO₂) derived by partial melting of eclogite before or during cratonization. The melt-peridotite interaction converted olivine into orthopyroxene.All three xenolith types have secondary metasomatic clinopyroxene and garnet, which occur along olivine grain boundaries and have an amoeboid texture. As indicated by the preservation of oxygen isotope disequilibrium in the minerals and trace-element concentrations in clinopyroxene and garnet, this metasomatic event is probably of Mesozoic age and was caused by percolating alkaline basaltic melts. This melt metasomatism enriched the xenoliths in CaO, Al₂O₃, FeO and high-field-strength-elements, and might correspond to the Karoo magmatism at 200 Ma. The websterite xenoliths experienced both the orthoyproxene-enrichment and clinopyroxene-garnet metasomatic events, whereas dunite and wehrlite xenoliths only saw the later basaltic melt event, and may have been situated further away from the source of melt migration channels.     

Pressure effect on the divalent cation distribution in nonideal solid solution of forsterite and fayalite

The divalent cation distribution in olivine (Mg, Fe)₂SiO₄ under high pressure and temperature was studied to clarify the detailed state of olivine in the mantle. Single crystal samples were heated for a sufficient period of time for the cations to migrate and quenched fast enough to preserve the equilibrated state under high pressures, and the crystal structure was determined with X-ray method. The pressure effect on the distribution coefficient K _D[= (Fe/Mg) _M1/(Fe/Mg) _M2] was determined for the first time; dK _D/dP?0.02 GPa^?1. A set of five thermodynamic parameters required to describe the regular solution model was determined from data concerning the pressure dependence and the known temperature and compositional effects. As a result we have shown how K _D depends on pressure, temperature, and composition. The notable feature clarified is the very large contribution of nonideality in the olivine solid solution. The K _D of olivine crystals in the mantle is predicted; K _D increases to ～ 2.2 at the depth of 400 km, in contrast to 0.9 ～ 1.2 of natural samples available at the surface of the Earth.     

Experimental determination of the role of diffusion on Li isotope fractionation during basaltic glass weathering

In order to use lithium isotopes as tracers of silicate weathering, it is of primary importance to determine the processes responsible for Li isotope fractionation and to constrain the isotope fractionation factors caused by each process as a function of environmental parameters (e.g. temperature, pH). The aim of this study is to assess Li isotope fractionation during the dissolution of basalt and particularly during leaching of Li into solution by diffusion or ion exchange. To this end, we performed dissolution experiments on a Li-enriched synthetic basaltic glass at low ratios of mineral surface area/volume of solution (S/V), over short timescales, at various temperatures (50 and 90 °C) and pH (3, 7, and 10). Analyses of the Li isotope composition of the resulting solutions show that the leachates are enriched in ⁶Li (δ⁷Li = +4.9 to +10.5‰) compared to the fresh basaltic glass (δ⁷Li = +10.3 ± 0.4‰). The δ⁷Li value of the leachate is lower during the early stages of the leaching process, increasing to values close to the fresh basaltic glass as leaching progresses. These low δ⁷Li values can be explained in terms of diffusion-driven isotope fractionation. In order to quantify the fractionation caused by diffusion, we have developed a model that couples Li diffusion with dissolution of the glassy silicate network. This model calculates the ratio of the diffusion coefficients of both isotopes (a = D₇/D₆), as well as its dependence on temperature, pH, and S/V. a is mainly dependent on temperature, which can be explained by a small difference in activation energy (0.10 ± 0.02 kJ/mol) between ⁶Li⁺ and ⁷Li⁺. This temperature dependence reveals that Li isotope fractionation during diffusion is low at low temperatures (T < 20 °C), but can be significant at high temperatures. However, concerning hydrothermal fluids (T > 120 °C), the dissolution rate of basaltic glass is also high and masks the effects of diffusion. These results indicate that the high δ⁷Li values of river waters, in particular in basaltic catchments, and the fractionated values of hydrothermal fluids are mainly controlled by precipitation of secondary phases.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

An experimental study of the kinetics of lherzolite reactive dissolution with applications to melt channel formation

The kinetics of lherzolite dissolution in an alkali basalt and a basaltic andesite was examined experimentally at 1,300°C and 1 GPa using the dissolution couple method. Dissolution of lherzolite in basaltic liquids produces either the melt-bearing dunite–harzburgite–lherzolite (DHL) sequence or the melt-bearing harzburgite–lherzolite sequence depending on whether the reacting melt is or close to olivine saturation (alkali basalt) or olivine + orthopyroxene saturation (basaltic andesite). The dunite in the DHL sequence is pyroxene-free and the harzburgites in both sequences are clinopyroxene-free. The melt fraction and olivine grain size in the dunite are larger than those in the harzburgite. The olivine grain size in the dunite and harzburgite in the DHL sequence also increases as a function experimental run time. Across the sharp dunite–harzburgite and harzburgite–lherzolite interfaces, systematic compositional variations are observed in the reacting melt, interstitial melt, olivine, and to a lesser extent, pyroxenes as functions of distance and time. The systematic variations in lithology, grain size, mineral chemistry, and melt compositions are broadly similar to those observed in the mantle sections of ophiolites. The processes of lherzolite dissolution in basaltic liquids involve dissolution, precipitation, reprecipitation, and diffusive transport in the interstitial melts and surrounding minerals. Preferential dissolution of olivine and clinopyroxene and precipitation of orthopyroxene in the basaltic andesite produces the melt-bearing harzburgite–lherzolite sequence. Preferential dissolution of clinopyroxene and orthopyroxene and precipitation of olivine results in the melt-bearing DHL sequence. Preferential mineral dissolution can also affect the composition of the through-going melt in a dunite channel or harzburgite matrix. Systematic variations in melt fraction and mineral grain size in the peridotite sequences are likely to play an important role in the development of channelized or diffuse porous melt flow in the mantle.An erratum to this article can be found at     

Evidence for heterogeneous enriched shergottite mantle sources in Mars from olivine-hosted melt inclusions in Larkman Nunatak 06319

Larkman Nunatak (LAR) 06319 is an olivine-phyric shergottite whose olivine crystals contain abundant crystallized melt inclusions. In this study, three types of melt inclusion were distinguished, based on their occurrence and the composition of their olivine host: Type-I inclusions occur in phenocryst cores (Fo_77-73); Type-II inclusions occur in phenocryst mantles (Fo_71-66); Type-III inclusions occur in phenocryst rims (Fo_61-51) and within groundmass olivine. The sizes of the melt inclusions decrease significantly from Type-I (∼150-250 μm diameter) to Type-II (∼100 μm diameter) to Type-III (∼25-75 μm diameter). Present bulk compositions (PBC) of the crystallized melt inclusions were calculated for each of the three melt inclusion types based on average modal abundances and analyzed compositions of constituent phases. Primary trapped liquid compositions were then reconstructed by addition of olivine and adjustment of the Fe/Mg ratio to equilibrium with the host olivine (to account for crystallization of wall olivine and the effects of Fe/Mg re-equilibration). The present bulk composition of Type-I inclusions (PBC1) plots on a tie-line that passes through olivine and the LAR 06319 whole-rock composition. The parent magma composition can be reconstructed by addition of 29 mol% olivine to PBC1, and adjustment of Fe/Mg for equilibrium with olivine of Fo₇₇ composition. The resulting parent magma composition has a predicted crystallization sequence that is consistent with that determined from petrographic observations, and differs significantly from the whole-rock only in an accumulated olivine component (∼10 wt%). This is consistent with a calculation indicating that ∼10 wt% magnesian (Fo_77-73) olivine must be subtracted from the whole-rock to yield a melt in equilibrium with Fo₇₇. Thus, two independent estimates indicate that LAR 06319 contains ∼10 wt% cumulate olivine.The rare earth element (REE) patterns of Type-I melt inclusions are similar to that of the LAR 06319 whole-rock. The REE patterns of Type-II and Type-III melt inclusions are also broadly parallel to that of the whole-rock, but at higher absolute abundances. These results are consistent with an LAR 06319 parent magma that crystallized as a closed-system, with its incompatible-element enrichment being inherited from its mantle source region. However, fractional crystallization of the reconstructed LAR 06319 parent magma cannot reproduce the major and trace element characteristics of all enriched basaltic shergottites, indicating local-to-large scale major- and trace-element variations in the mantle source of enriched shergottites. Therefore, LAR 06319 cannot be parental to the enriched basaltic shergottites.     

Water diffusion in dacitic melt

H₂O diffusion in dacitic melt was investigated at 0.48-0.95 GPa and 786-893 K in a piston-cylinder apparatus. The diffusion couple design was used, in which a nominally dry dacitic glass makes one half and is juxtaposed with a hydrous dacitic glass containing up to ∼8 wt.% total water (H₂O_t). H₂O concentration profiles were measured on quenched glasses with infrared microspectroscopy. The H₂O diffusivity in dacite increases rapidly with water content under experimental conditions, similar to previous measurements at the same temperature but at pressure <0.15 GPa. However, compared with the low-pressure data, H₂O diffusion at high pressure is systematically slower. H₂O diffusion profiles in dacite can be modeled by assuming molecular H₂O (H₂O_m) is the diffusing species. Total H₂O diffusivity D_H2Ot within 786-1798 K, 0-1 GPa, and 0-8 wt.% H₂O_t can be expressed as: where D_H2Ot is in m²/s, T is temperature in K, P is pressure in GPa, K = exp(1.49 − 2634/T) is the equilibrium constant of speciation reaction (H₂O_m+O?2OH) in the melt, X = C/18.015/[C/18.015 + (100 − C)/33.82], C is wt.% of H₂O_t, and 18.015 and 33.82 g/mol correspond to the molar masses of H₂O and anhydrous dacite on a single oxygen basis. Compared to H₂O diffusion in rhyolite, diffusivity in dacite is lower at intermediate temperatures but higher at superliquidus temperatures. This general H₂O diffusivity expression can be applied to a broad range of geological conditions, including both magma chamber processes and volcanic eruption dynamics from conduit to the surface.     

Partitioning of Ni2+, Co2+, Fe2+, Mn2+ and Mg2+ between olivine and silicate melts: compositional dependence of partition coefficient

Partitioning of Ni²⁺, Co²⁺, Fe²⁺, Mn²⁺ and Mg²⁺ between olivine and silicate melts has been determined near the join (Mg_0.5·-Fe_0.5)₂SiO₄-K₂O 4SiO₂ and for seven different basaltic compositions. The experiments were made at 1 atm total pressure, 1500-1100°C, and under moderate to reducing oxygen fugacities. The concentration factor, defined as K^MO = (MO)^ol/(MO)^liq (molar ratio), increases markedly for all the cations studied as the olivine component of the liquid decreases. Much of the increase in K^MO is considered as due to the compositional effect of the coexisting liquid: the temperature effect on K^MO is probably opposite to the compositional effect (K^MO decreases as temperature decreases).The partition coefficient K^MO-MgO = (MO/MgO)^ol/(MO/MgO)^liq for the reaction, M_ol²⁺ + Mg_liq²⁺ = M_liq²⁺ + Mg_ol²⁺. is relatively constant over a wide range of SiO₂ content of the liquid, except in the case of Ni²⁺. The partition coefficients have similar ranges both in synthetic and natural rock systems: K^NiO-MgO = 1.8–3.0, K^CoO-MgO = 0.6–0.8, K^FeO-MgO = 0.27–0.38, and K^MnO-MgO = 0.23–0.32. There is a systematic variation in the partition coefficient K^MO-MgO with the composition of liquid; K^MO-MgO increases with increasing SiO₂ content of melt. The partition coefficients for Co²⁺, Fe²⁺ and Mn²⁺ are useful to test the equilibration of olivine with magma of a wide compositional range.     

Component mobility at 900 °C and 18 kbar from experimentally grown coronas in a natural gabbro

Several approximately 100-μm-wide reaction zones were grown under experimental conditions of 900 °C and 18 kbar along former olivine-plagioclase contacts in a natural gabbro. The reaction zone comprises two distinct domains: (i) an irregularly bounded zone with idiomorphic grains of zoisite and minor corundum and kyanite immersed in a melt developed at the plagioclase side and (ii) a well-defined reaction band comprising a succession of mineral layers forming a corona structure around olivine. Between the olivine and the plagioclase reactant phases we observe the following layer sequence: olivine|pyroxene|garnet|partially molten domain|plagioclase. Within the pyroxene layer two micro-structurally distinct layers comprising enstatite and clinopyroxene can be discerned. Chemical potential gradients persisted for the CaO, Al₂O₃, SiO₂, MgO and FeO components, which drove diffusion of Ca, Al and Si bearing species from the garnet-matrix interface to the pyroxene-olivine interface and diffusion of Mg- and Fe-bearing species in the opposite direction. The systematic mineralogical organization and chemical zoning across the corona suggest that the olivine corona was formed by a “diffusion-controlled” reaction. We estimate a set of diffusion coefficients and conclude that L_AlAl < L_CaCa < (L_SiSi, L_FeFe) < L_MgMg during reaction rim growth.     

Experimental study of polybaric REE partitioning between olivine, pyroxene and melt of the Yamato 980459 composition: Insights into the petrogenesis of depleted shergottites

A synthetic composition representing the Yamato 980459 martian basalt (shergottite) has been used to carry out phase relation, and rare earth element (REE) olivine and pyroxene partitioning experiments. Yamato 980459 is a sample of primitive basalt derived from a reduced end-member among martian mantle sources. Experiments carried out between 1-2 GPa and 1350-1650 °C simulate the estimated pressure-temperature conditions of basaltic melt generation in the martian mantle. Olivine-melt and orthopyroxene-melt partition coefficients for La, Nd, Sm, Eu, Gd and Yb (D_REE values) were determined by LA-ICPMS, and are similar to the published values for terrestrial basaltic systems. We have not detected significant variation in D-values with pressure over the range investigated, and by comparison with previous studies carried out at lower pressure.We apply the experimentally obtained olivine-melt and orthopyroxene-melt D_REE values to fractional crystallization and partial melting models to develop a three-stage geochemical model for the evolution of martian meteorites. In our model we propose two ancient (∼4.535 Ga) sources: the Nakhlite Source, located in the shallow mantle, and the Deep Mantle Source, located close to the martian core-mantle boundary. These two sources evolved distinctly on the ε¹⁴³Nd evolution curve due to their different Sm/Nd ratios. By partially melting the Nakhlite Source at ∼1.3 Ga, we are able to produce a slightly depleted residue (Nakhlite Residue). The Nakhlite Residue is left undisturbed until ∼500 Ma, at which point the depleted Deep Mantle Source is brought up by a plume mechanism carrying with it high heat flow, melts and isotopic signatures of the deep mantle (e.g., ε¹⁸²W, ε¹⁴²Nd, etc.). The plume-derived Deep Mantle Source combines with the Nakhlite Residue producing a mixture that becomes a mantle source (herein referred to as “the Y98 source”) for Yamato 980459 and the other depleted shergottites with the characteristic range of Sm/Nd ratios of these meteorites. The same hot plume provides a heat source for the formation of enriched and intermediate shergottites. Our model reproduces the REE patterns of nakhlites and depleted shergottites and can explain high ε¹⁴³Nd in depleted shergottites. Furthermore, the model results can be used to interpret whole rock Rb-Sr and Sm-Nd ages of shergottites.