Fluid-mineral equilibria in a hydrothermal system,Roosevelt hot springs,Utah

The availability of fluids and drill cuttings from the active hydrothermal system at Roosevelt Hot Springs allows a quantitative comparison between the observed and predicted alteration mineralogy, calculated from fluid-mineral equilibria relationships. Comparison of all wells and springs in the thermal area indicates a common reservoir source, and geothermometer calculations predict its temperature to be higher ($\text{288°C ± 10°}$) than the maximum measured temperature of 268°C.The composition of the deep reservoir fluid was estimated from surface well samples, allowing for steam loss, gas release, mineral precipitation and ground-water mixing in the well bore. This deep fluid is sodium chloride in character, with approximately 9700 ppm dissolved solids, a pH of 6.0, and gas partial pressures of O₂ ranging from 10^?32 to 10^?35 atm, CO₂ of 11 atm, H₂S of 0.020 atm and CH₄ of 0.001 atm.Comparison of the alteration mineralogy from producing and nonproducing wells allowed delineation of an alteration pattern characteristic of the reservoir rock. Theoretical alteration mineral assemblages in equilibrium with the deep reservoir fluid, between 150° and 300°C, in the system Na₂O-K₂O-CaO-MgO-FeO-Fe₂O₃-Al₂O₃-H₄SiO₄-H₂O-H₂S-CO₂-HCl, were calculated. Minerals theoretically in equilibrium with the calculated reservoir fluid at $\text{>240°C}$ include sericite, K-feldspar, quartz, chalcedony, hematite, magnetite and pyrite. This assemblage corresponds with observed higher-temperature ($\text{>210°C}$) alteration assemblage in the deeper parts of the producing wells. The presence of montmorillonite and mixed-layer clays with the above assemblage observed at temperatures <210°C corresponds with minerals predicted to be in equilibrium with the fluid below 240°C.Alteration minerals present in the reservoir rock that do not exhibit equilibrium with respect to the reservoir fluid include epidote, anhydrite, calcite and chlorite. These may be products of an earlier hydrothermal event, or processes such as boiling and mixing, or a result of errors in the equilibrium calculations as a result of inadequate thermochemical data.     

The mineral chemistry of new titanates from the jagersfontein kimberlite,South Africa: Implications for metasomatism in the upper mantle

New members of the crichtonite mineral series are described in which K, Ba, Ca and REE are in significant concentrations (5 wt% oxides) filling the A formula position in AM₂₁O₃₈. These phases are chromium (16 wt% Cr₂O₃) titanates (58 wt% TiO₂) enriched in ZrO₂ (5 wt%) and constitute a mineral repository for refractory and large ion lithophile elements in the upper mantle. The mineral senes coexists with Mg-Cr-ilmenite, Nb-Cr-rutile, and Ca-Cr (NbZr) armalcolite that have equally unusual chemistries. Kimberlitic crichtonites are depleted in the intermediate lanthanides but highly enriched in LREE and HREE with chondrite normalized abundances of 10³ to 10⁵. Crichtonite, armalcolite, and Nb-Cr-rutile occupy a compositional range in TiO₂ contents bridging the gap between ilmenite and rutile, two minerals having a widespread distribution in kimberlites and mantle-derived nodule suites.In common with other associations, and based on similarities in mineral chemistry, it is concluded that these minerals formed at P = 20–30 kb, 900–1100°C by reaction of peridotite with metasomatizing fluids. Kimberlitic crichtonite may be expressed as spinel + Cr-ferropseudobrookite, and armalcolite is equivalent to Cr-geikielite + rutile in the system (FeMg)-TiO₂-Cr₂O₃. This system contains a number of Cr-Ti compounds not found as minerals but it is proposed that the ubiquitous occurrence of ilmenite intergrowths in kimberlitic rutile results from decomposition of high pressure αPbO₂-type crystallographic shear structures. The new minerals have exotic chemistries and the high K-affinities broaden the scope for the origin of alkalic rocks, the generation of highly potassic magmas in the upper mantle, and suggest that alkali metasomatism may be pervasive.     

Crystallization of titaniferous chromite,magnesian ilmenite and armalcolite in tholeiitic suites in the Karoo Igneous Province

Titaniferous chromite (up to 8 wt% TiO₂) and magnesian ilmenite (up to 10 wt% MgO) coexist at the base of the differentiated tholeiitic Mount Ayliff Intrusion in the Karoo Province of southern Africa, suggesting that the original magma was TiO₂-rich. Picritic lavas with 3% TiO₂ from the Lebombo monocline of the Karoo Province also contain microphenocrysts of magnesian ilmenite (up to 6 wt% mgO) and armalcolite (up to 7 wt% MgO). These oxide mineral associations and compositions are atypical of tholeiitic magmas, in which chromite usually has less than 1 wt% TiO₂, ilmenite less than 3 wt% MgO and armalcolite is rarely a primary mineral. Experiments have been conducted at one atmosphere pressure on a range of compositions to determine the effect of TiO₂ on the crystallization and composition of chromite, ilmenite and armalcolite. The results indicate that increasing the TiO₂ content of picritic magmas increases the TiO₂ content of the spinel, mainly at the expense of Al₂O₃, whereas Cr₂O₃ is not affected. Spinel compositions in the Mount Ayliff Intrusion (with over 45 wt% Cr₂O₃, less than 10 wt% Al₂O₃ and 8 wt% TiO₂) were duplicated in experiments on a picrite at temperatures of about 1,200°C at the Ni/NiO buffer. Increasing fO₂ from fayalite-magnetite-quartz to Ni/NiO buffer is shown to increase the crystallization temperature of armalcolite and to decrease that of ilmenite. The total FeO content of the liquid has little influence on the crystallization temperature of these phases. The TiO₂ content of the liquid, when either ilmenite or armalcolite crystallizes, varies inversely with SiO₂ content. The MgO content of the liquid at which ilmenite or armalcolite crystallizes depends upon the TiO₂ content of the starting composition, with naturally occurring and experimetally determined saturation being demonstrated for liquids with 5 wt% MgO and 5.5 wt% TiO₂. The partition coefficent for MgO between armalcolite or ilmenite and liquid is about 1.5. Observed magnesian armalcolite and ilmenite compositions in picrite lavas (both minerals) and in the Mount Ayliff Intrusion (ilmenite only) are consistent with crystallization from a TiO₂-rich magma with approximately 5 wt% MgO. The Fe ₂ ³⁺ TiO₅ component of armalcolite in the picrite lavas matches those formed experimentally at temperatures of 1,150–1,110°C and fO₂ of the Ni/NiO to Ni/NiO+1 log unit. Similarities also exist between the compositions of chromite, ilmenite and armalcolite and liquid fraction-ation trends of some Hawaiian high-TiO₂ lavas and the experimental studies presented here.     

Rutile solubility and titanium coordination in silicate melts

The solubility of rutile has been determined in a series of compositions in the K₂O-Al₂O₃-SiO₂ system ($\text{K}{}^1\text{=}\text{K}{}_2\text{O}\text{(K}{}_2\text{O}$ + Al₂O₃) = 0.38–0.90), and the CaO-Al₂O₃-SiO₂ system ($\text{C}{}^1\text{=}\text{CaO}\text{(CaO + Al}{}_2\text{O}{}_3\text{)}\text{= 0.47–0.59}$). Isothermal results in the KAS system at 1325°C, 1400°C, and 1475°C show rutile solubility to be a strong function of the $\text{K}{}^1$ ratio. For example, at 1475°C the amount of TiO₂ required for rutile saturation varies from 9.5 wt% ($\text{K}{}^1\text{= 0.38}$) to 11.5 wt% ($\text{K}{}^1\text{= 0.48}$) to 41.2 wt% ($\text{K}{}^1\text{= 0.90}$). In the CAS system at 1475°C, rutile solubility is not a strong function of $\text{C}{}^1$. The amount of TiO₂ required for saturation varies from 14 wt% ($\text{C}{}^1\text{= 0.48}$) to 16.2 wt% ($\text{C}{}^1\text{= 0.59}$).The solubility changes in KAS melts are interpreted to be due to the formation of strong complexes between Ti and K⁺ in excess of that needed to charge balance Al³⁺. The suggested stoichiometry of this complex is K₂Ti₂O₅ or K₂Ti₃O₇. In CAS melts, the data suggest that Ca²⁺ in excess of A1³⁺ is not as effective at complexing with Ti as is K⁺. The greater solubility of rutile in CAS melts when $\text{C}{}^1$ is less than 0.54 compared to KAS melts of equal $\text{K}{}^1$ ratio results primarily from competition between Ti and Al for complexing cations (Ca vs. K).TiK_β x-ray emission spectra of KAS glasses ($\text{K}{}^1\text{= 0.43–0.60}$) with 7 mole% added TiO₂, rutile, and Ba₂TiO₄, demonstrate that the average Ti-O bond length in these glasses is equal to that of rutile rather than Ba₂TiO₄, implying that Ti in these compositions is 6-fold rather than 4-fold coordinated. Re-examination of published spectroscopic data in light of these results and the solubility data, suggests that the 6-fold coordination polyhedron of Ti is highly distorted, with at least one Ti-O bond grossly undersatisfied in terms of Pauling's rules.     

The system pargasite-H2O-CO2: a model for melting of a hydrous mineral with a mixed-volatile fluid—I. Experimental results to 8 kbar

The stability of the amphibole pargasite [NaCa₂Mg₄Al(Al₂Si₆))O₂₂(OH)₂] in the melting range has been determined at total pressures (P) of 1.2 to 8 kbar. The activity of H₂O was controlled independently of P by using mixtures of H₂O + CO₂ in the fluid phase. The mole fraction of H₂O in the fluid ($\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$) ranged from 1.0 to 0.2.At P < 4 kbar the stability temperature (T) of pargasite decreases with decreasing $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at constant P. Above P ? 4 kbar stability T increases as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ is decreased below one, passes through a T maximum and then decreases with a further decrease in $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$. This behavior is due to a decrease in the H₂O content of the silicate liquid as $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ decreases. The magnitude of the T maximum increases from about 10°C (relative to the stability T for $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{= 1}$) at P = 5 kbar to about 30°C at P = 8 kbar, and the position of the maximum shifts from $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.6}$ at P = 5 kbar to $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}\text{? 0.4}$ at P = 8 kbar.The H₂O content of liquid coexisting with pargasite has been estimated as a function of at 5 and 8 kbar P, and can be used to estimate the H₂O content of magmas. Because pargasite is stable at low values of $\text{X}{}_{\mathrm{H}}{}_2\text{O}{}^{\mathrm{1fl}}$ at high P and T, hornblende can be an important phase in igneous processes even at relatively low H₂O fugacities.     

Sorption of noble gases by solids,with reference to meteorites. I. Magnetite and carbon

To simulate trapping of meteoritic noble gases by solids, 18 samples of Fe₃O₄ were synthesized in a noble gas atmosphere at 350–720 K by the reactions: 3Fe + 4H₂O → Fe₃O₄ + 4H₂ (Ne, Ar, Kr, Xe) 3Fe + 4CO → Fe₃O₃ + 4C + carbides (Xe only) Phases were separated by selective solvents (HgCl₂, HCl). Noble gas contents were analyzed by mass spectrometry, or, in runs where 36 d Xe¹²⁷ tracer was used, by γ-counting. Surface areas, as measured by the BET method, ranged from 1 to 400 m²/g. Isotopic fractionations were below the detection limit of 0.5%/m.u.Sorption of Xe on Fe₃O₄ and C obeys Henry's Law between $\text{1 × 10}{}^{\mathrm{?8}}$ and $\text{4 × 10}{}^{\mathrm{?5}}$ atm, but shows only a slight temperature dependence between 650 and 720 K ($\text{ΔH}{}_{\mathrm{sol}}\text{= ?4 ± 2}\text{kcal/mole}$). The mean distribution coefficient K_Xe is $\text{0.28 ± 0.09}$ cc STP/g atm for Fe₃O₄ and only a factor of $\text{1.2 ± 0.4}$ greater for C; such similarity for two cogenetic phases was predicted by Lewis et al. (1977). Stepped heating and etching experiments show that 20–50% of the total Xe is physically adsorbed and about 20% is trapped in the solid. The rest is chemisorbed with $\text{ΔH}{}_{\mathrm{s}}\text{? ?13}\text{kcal/mole}$. The desorption or exchange half-time for the last two components is >10² yr at room temperature.Etching experiments showed a possible analogy to “Phase $\text{Q}$” in meteorites. A typical carbon + carbide sample, when etched with HNO₃, lost 47% of its Xe but only 0.9% of its mass, corresponding to a ~0.6 Å layer. Though this etchable, surficial gas component was more thermolabile than $\text{Q}$ (release $\text{T}$ below 1000°C, compared to 1200–1600°C), another experiment shows that the proportion of chemisorbed Xe increases upon moderate heating (1 hr at 450°C). Apparently adsorbed gases can become “fixed” to the crystal, by processes not involving volume diffusion (recrystallization, chemical reaction, migration to traps, etc.). Such mechanisms may have acted in the solar nebula, to strengthen the binding of adsorbed gases.Adsorbed atmospheric noble gases are present in all samples, and dominate whenever the noble gas partial pressure in the atmosphere is greater than that in the synthesis. Many of the results of Lancet and Anders (1973) seem to have been dominated by such an atmospheric component; others are suspect for other reasons, whereas still others seem reliable. When the doubtful samples of Lancet and Anders are eliminated or corrected, the fractionation pattern—as in our samples—no longer peaks at Ar, but rises monotonically from Ne to Xe. No clear evidence remains for the strong temperature dependence claimed by these authors.     

Separate or shared metamorphic histories of eclogites and surrounding rocks? An example from the Bohemian Massif

Eclogite boudins occur within an orthogneiss sheet enclosed in a Barrovian metapelite‐dominated volcano‐sedimentary sequence within the Velké Vrbno unit, NE Bohemian Massif. A metamorphic and lithological break defines the base of the eclogite‐bearing orthogneiss nappe, with a structurally lower sequence without eclogite exposed in a tectonic window. The typical assemblage of the structurally upper metapelites is garnet–staurolite–kyanite–biotite–plagioclase–muscovite–quartz–ilmenite ± rutile ± silli‐manite and prograde‐zoned garnet includes chloritoid–chlorite–paragonite–margarite, staurolite–chlorite–paragonite–margarite and kyanite–chlorite–rutile. In pseudosection modelling in the system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O (NCKFMASH) using THERMOCALC, the prograde path crosses the discontinuous reaction chloritoid + margarite = chlorite + garnet + staurolite + paragonite (with muscovite + quartz + H₂O) at 9.5 kbar and 570 °C and the metamorphic peak is reached at 11 kbar and 640 °C. Decompression through about 7 kbar is indicated by sillimanite and biotite growing at the expense of garnet. In the tectonic window, the structurally lower metapelites (garnet–staurolite–biotite–muscovite–quartz ± plagioclase ± sillimanite ± kyanite) and amphibolites (garnet–amphibole–plagioclase ± epidote) indicate a metamorphic peak of 10 kbar at 620 °C and 11 kbar and 610–660 °C, respectively, that is consistent with the other metapelites. The eclogites are composed of garnet, omphacite relicts (jadeite = 33%) within plagioclase–clinopyroxene symplectites, epidote and late amphibole–plagioclase domains. Garnet commonly includes rutile–quartz–epidote ± clinopyroxene (jadeite = 43%) ± magnetite ± amphibole and its growth zoning is compatible in the pseudosection with burial under H₂O‐undersaturated conditions to 18 kbar and 680 °C. Plagioclase + amphibole replaces garnet within foliated boudin margins and results in the assemblage epidote–amphibole–plagioclase indicating that decompression occurred under decreasing temperature into garnet‐free epidote–amphibolite facies conditions. The prograde path of eclogites and metapelites up to the metamorphic peak cannot be shared, being along different geothermal gradients, of about 11 and 17 °C km^?1, respectively, to metamorphic pressure peaks that are 6–7 kbar apart. The eclogite–orthogneiss sheet docked with metapelites at about 11 kbar and 650 °C, and from this depth the exhumation of the pile is shared.     

Effect of pressure on the diffusivity of network-forming cations in melts of jadeitic compositions

The diffusivities of network-forming cations (Si⁴⁺, Al³⁺, Ge⁴⁺ and Ga³⁺) in melts of the jadeitic composition NaAl(Si, Ge)₂O₆ and Na(Al, Ga)Si₂O₆ have been measured at pressures between 6 and 20 kbar at 1400°C. The rates of interdiffusion of Si⁴⁺-Ge⁴⁺ and Al³⁺-Ge³⁺ increase with increasing pressure at constant temperature. The results are consistent with the ion-dynamics computer simulations of Jadeite melt by Angellet al. (1982, 1983). The coefficient measured for the Si⁴⁺-Ge⁴⁺ interdiffusion is between 8 × 10^?10 and 2.5 × 10^?8$\text{cm}{}^2\text{sec}$ at 6 kbar, depending on the composition of the melt, whereas at 20 kbar it is between 7 × 10^?9 and 2 × 10^?7$\text{cm}{}^2\text{sec}$. The effect of pressure is greater for more Si-rich compositions (i.e., closer to NaAlSi₂O₆ composition). The coefficient measured for the Al³⁺-Ga³⁺ inter- diffusion is between 9 × 10^?10 and 3 × 10^?9 cm²/sec at 6 kbar and between 3 × 10^?9 and 1 × 10^?8$\text{cm}{}^2\text{sec}$ at 20 kbar. The rate of increase in diffusivity with pressure of Al³⁺-Ga³⁺ (a factor of 3–4) is smaller than that of Si⁴⁺-Ge⁴⁺ (a factor of 7–17).The Si⁴⁺-Ge⁴⁺ interdiffusion in melts of Na₂O · 4(Si, Ge)O₂ composition has also been measured at 8 and 15 kbar for comparison. The effect of pressure on the diffusivity in this melt is significantly smaller than that for the jadeitic melts. The increase in diffusivity of the network-forming cations in jadeitic melts with increasing pressure may be related to the decrease in viscosity of the same melt. The present results, as well as the ion-dynamics simulations, suggest that the homogenization of partial melts and mixing of magmas would be more efficient at greater depths.     

Experimental determination of the effects of pressure and temperature on the stoichiometry and phase relations of wüstite

When quenched metastable wüstite (Fe_.924O and Fe_.947O) is held at 300°C at pressures up to 200 kbar in a diamond anvil cell, a mixture of magnetite, metallic iron and wüstite is found. We interpret this to indicate that magnetite plus metallic iron constitute the stable phase assemblage at pressures and temperatures below this boundary is stoichiometric FeO ($\text{a}{}_0\text{= 4.332 ± 0.001}\text{A}\text{?}$) at pressures below 110 kbar at 300°C. However, just below the boundary in the pressure range 110 kbar to 200 kbar at 300°C, the residuál wüstite is non-stoichiometric ($\text{a}{}_0\text{< 4.332}\text{A}\text{?}$). Data collected at pressures and temperatures above the boundary indicate that non-stoichiometric wüstite (Fe_xO) plus metallic iron constitute the stable phase assemblage and that the value of x in Fe_xO increases as pressure is increased isothermally to 100 kbar and then decreases as pressure is increased above 100 kbar.     

Subsolidus and melting relations for the join CaCO3-MgCO3 at 10 kbar

Mixtures of pure dry CaCO₃ and MgCO₃ were reacted at 10 kbar in a piston-cylinder apparatus. Solidus and liquidus boundaries were delineated by interpretation of quenched textures. X-ray determined compositions of quenched carbonates are not a reliable guide to the phase relations. The binary melting loop for CaCO₃-MgCO₃ extends from CaCO₃ at 1460°C through a liquidus minimum near 30 wt% MgCO₃ and 1075°C, and it is terminated at the incongruent melting reaction for dolomite solid solution at 1125° C (liquid with 32 wt% MgCO₃) Magnesite solid solution dissociates at 1090°C to produce dolomite + periclase + CO₂, truncating the dolomite-magnesite solvus. The 10 kb liquidus minimum at 1075°C and 30 wt% MgCO₃ occurs at lower temperature and higher $\text{Ca}\text{Mg}$ ratio than the 27 kbar liquidus minimum at 1290°C and 38 wt% MgCO₃. This relationship suggests that the first liquid produced by melting of a carbonate-bearing peridotite has increasing $\text{Mg}\text{Ca}$ ratio with increasing pressure. These phase relations provide part of the framework required to trace paths of crystallization of kimberlite and carbonatite magmas.     

Solubility of sulfur in some magmas at 1 atmosphere

In order to understand the distribution of sulfur in igneous rooks, we determined the solubility of sulfur in volcanic rock melts (tholeiite basalt, hawaiite and rhyodacite from Hawaii) at various gas compositions and at 1250° and 1300°C and 1 atm total pressure. The solubility of sulfur in the melt passes through a minimum with change in oxygen partial pressure, if other factors are held constant. For the basaltic liquid at 1200°C, most sulfur in the melt is as dissolved sulfide (S^?2) at oxygen partial pressures below 10^?8 atm and as dissolved sulfate at oxygen partial pressures above 10^?8 atm. Based on the present solubility data, 5 per cent is inferred for volcanic gas at 1 atm total pressure in equilibrium with subaerially extruded Hawaiian tholeiite basalt (Pele's hair with 180 ppm S) at 1200°C and 10^?8 atm P_O2.     

Volatilization of sodium from silicate melt spheres and its application to the formation of chondrules

The rates of volatilization of Na from liquid spheres of chondrule compositions have been determined as functions of time, temperature, partial pressure of oxygen, and sizes of the spheres. The Na₂O content in the sphere is uniform in each run. but it decreases with time of the run, indicating that the rate of diffusion of Na in the liquid is greater than that of volatilization, and that the latter is the rate-controlling process. The rate of sodium volatilization becomes greater with increasing temperature and with decreasing P_O2 and size of the spheres. The relation of the Na₂O content in the liquid sphere with time and its size indicate that the amount of Na₂O volatilized from the liquid spheres within unit time is proportional to the surface area of the spheres and the concentration of Na₂O in the liquid. From these relations, the rate of volatilization of sodium can be obtained at constant temperature and P_o2. The rate of volatilization of sodium satisfies the Arrhenius relation within the temperature range from about 1450–1600 C at 10^?9,2 atm p_O2; the activation energy for the sodium volatilization is approximately 100 kcal-mole^?1. The rate is also approximately proportional to within the range of p_O2 from 10^?10.2 to 10^?5.0 atm at about 1500° C. Based on the present results and the Na₂O contents in chondrules. it is suggested that they experienced an instant heating with maximum temperature of 1400–2200° C followed by an immediate cooling.     

Estimation of the dielectric constant of H2O from experimental solubilities of quartz,and calculation of the thermodynamic properties of aqueous species to 900°C at 2 kb

Experimental quartz solubilities in H₂O (Anderson and Burnham, 1965, 1967) were used together with equations of state for quartz and aqueous species (Helgesonet al., 1978; Walther and Helgeson, 1977) to calculate the dielectric constant of H₂O () at pressures and temperatures greater than those for which experimental measurements (Heger, 1969; Lukashovet al., 1975) are available ($\text{0.001 ? P ? 5}\text{kb}$ and $\text{0 ? T ? 600°C}$). Estimates of computed in this way for 2 kb (which are the most reliable) range from 9.6 at 600°C to 5.6 at 800°C. These values are 0.5 and 0.8 units greater, respectively, than corresponding values estimated by Quist and Marshall (1965), but they differ by <0.3 units from extrapolated values computed from Pitzer's (1983) adaptation of the Kirkwood (1939) equation. The estimates of generated from quartz solubilities at 2 kb were fit with a power function of temperature, which was then used together with equations and data given by Helgeson and Kirkham (1974a,b, 1976) Helgesonet al. (1981), and Helgeson (1982b, 1984) to calculate Born functions, Debye Hückel parameters, and the thermodynamic properties of Na⁺, K⁺, Mg⁺⁺, Ca⁺⁺, and other aqueous species of geologic interest at temperatures to 900°C.     

Pro- and retrograde P–T evolution of granulites of the Beit Bridge Complex (Limpopo Belt, South Africa): constraints from quantitative phase diagrams and geotectonic implications

Interpretations based on quantitative phase diagrams in the system CaO–Na₂O–K₂O–TiO₂–MnO–FeO–MgO–Al₂O₃–SiO₂–H₂O indicate that mineral assemblages, zonations and microstructures observed in migmatitic rocks from the Beit Bridge Complex (Messina area, Limpopo Belt) formed along a clockwise P–T path. That path displays a prograde P–T increase from 600 °C/7.0 kbar to 780 °C/9–10 kbar (pressure peak) and 820 °C/8 kbar (thermal peak), followed by a P–T decrease to 600 °C/4 kbar. The data used to construct the P–T path were derived from three samples of migmatitic gneiss from a restricted area, each of which has a distinct bulk composition: (1) a K, Al‐rich garnet–biotite–cordierite–sillimanite–K‐feldspar–plagioclase–quartz–graphite gneiss (2) a K‐poor, Al‐rich garnet–biotite–staurolite–cordierite–kyanite–sillimanite–plagioclase–quartz–rutile gneiss, and (3) a K, Al‐poor, Fe‐rich garnet–orthopyroxene–biotite–chlorite–plagioclase–quartz–rutile–ilmenite gneiss. Preservation of continuous prograde garnet growth zonation demonstrates that the pro‐ and retrograde P–T evolution of the gneisses must have been rapid, occurring during a single orogenic cycle. These petrological findings in combination with existing geochronological and structural data show that granulite facies metamorphism of the Beit Bridge metasedimentary rocks resulted from an orogenic event during the Palaeoproterozoic (c. 2.0 Ga), caused by oblique collision between the Kaapvaal and Zimbabwe Cratons. Abbreviations follow Kretz (1983 ).     

An empirical solvus for CO2-H2O-2.6 wt% salt

The solvus in the system CO₂-H₂O-2.6 wt% NaCl-equivalent was determined by measuring temperature of homogenization in fluid inclusions which contained variable $\text{CO}{}_2\text{H}{}_2\text{O}$ but the same amount of salt dissolved in the aqueous phase at room temperature. The critical point of the solvus is at 340 ± 5°C, at pressures between 1 and 2 kbar; this is about 65°C higher than for the pure CO₂-H₂O system. The solvus is assymetrical, with a steeper H₂O-rich limb and with the critical point at mole fraction of water between 0.65 and 0.8.     

Thermochemistry of the high structural state plagioclases

The enthalpies of solution of a suite of 19 high-structural state synthetic plagioclases were measured in a Pb₂B₂O₅ melt at 970 K. The samples were crystallized from analyzed glasses at 1200°C and 20 kbar pressure in a piston-cylinder apparatus. A number of runs were also made on Amelia albite and Amelia albite synthetically disordered at 1050–1080°C and one bar for one month and at 1200°C and 20 kbar for 10 hr. The component oxides of anorthite, CaO, Al₂O₃ and SiO₂, were remeasured.The ΔH of disorder of albite inferred in the present study from albite crystallized from glass is 3.23 kcal, which agrees with the 3.4 found by Holm and Kleppa (1968). It is not certain whether this value includes the ΔH of a reversible displacive transition to monoclinic symmetry, as suggested by Helgesonet al. (1978) for the Holm-Kleppa results. The enthalpy of solution value for albite accepted for the solid solution series is based on the heat-treated Amelia albite and is 2.86 kcal less than for untreated Amelia albite.The enthalpy of formation from the oxides at 970 K of synthetic anorthite is ?24.06 ± 0.31 kcal, significantly higher than the ?23.16 kcal found by Charluet al. (1978), and in good agreement with the value of ?23.89 ± 0.82 given by Robieet al. (1979), based on acid calorimetry.The excess enthalpy of mixing in high plagioclase can be represented by the expression, valid at 970 K: ΔH^ex(±0.16 kcal) = 6.7461 X_abX²_An + 2.0247 X_AnX²_Ab where X_Ab and X_An are, respectively, the mole fractions of NaAlSi₃O₈ and CaAl₂Si₂O₈. This ΔH^ex, together with the mixing entropy of Kerrick and Darken's (1975) Al-avoidance model, reproduces almost perfectly the free energy of mixing found by Orville (1972) in aqueous cation-exchange experiments at 700°C. It is likely that Al-avoidance is the significant stabilizing factor in the high plagioclase series, at least for X_An≥ 0.3. At high temperatures the plagioclases have nearly the free energies of ideal one-site solid solutions. The Al-avoidance model leads to the following Gibbs energy of mixing for the high plagioclase series: $\text{ΔG}{}^{\mathrm{mix}}\text{= ΔH}{}^{\mathrm{ex}}\text{+ RT X}{}_{\mathrm{Ab}}\text{ln[X}{}^2{}_{\mathrm{Ab}}\text{(2 ? X}{}_{\mathrm{Ab}}\text{)]+ X}{}_{\mathrm{An}}\text{ln}\text{[X}{}_{\mathrm{An}}\text{(1+X}{}_{\mathrm{An}}\text{)}{}^2\text{]}\text{4}$. The entropy and enthalpy of mixing should be very nearly independent of temperature because of the unlikelihood of excess heat capacity in the albite-anorthite join.     

Pressure dependence of structural parameters of paragonite

The compressibility and structure of a 2M₁ paragonite with composition [Na_0.88K_0.10Ca_0.01Ba_0.01] [Al_1.97Ti_0.007Fe_0.01Mn_0.002Mg_0.006]Si_3.01Al_0.99O₁₀OH₂ were determined at pressures between 1 bar and 41 kbar, by single crystal X-ray diffraction using a Merrill-Bassett diamond anvil cell. Compressibility turned out to be largely anisotropic, linear compressibility coefficients parallel to the unit cell edges being β_a=3.5(1)·10^?4, β_b=3.6(1)·10^?4, β_c=8.3(3)·10^?4 kbar^?1 (β_a:β_b·β_c=1:1028:2.371). The isothermal bulk modulus, calculated as the reciprocal of the mean compressibility of the cell volume, was 650(20) kbar. The main features of the deformation mechanism resulting from structural refinements at pressures of 0.5, 25.4, 40.5 kbar were: –?variation in sheet thickness, showing that compression of the c parameter was mainly due to the interlayer thickness reduction from 3.07 Å at 0.5 kbar to 2.81 Å at 40.5 kbar; –?the compressibility of octahedra was greater than that of tetrahedra, the dimensional misfit between tetrahedral and octahedral sheets increased with P, so that tetrahedral rotation angel α increased from 15° at 0.5 kbar to 21.6° at 40.5 kbar; –?the basal surface corrugation (Δz) of the tetrahedral layer, due to the different dimensions of M1 and M2 octahedra and to the octahedral distortion, decreased with P (Δz=0.19 and 0.12 Å at 0.5 and 40.5 kbar respectively). Comparison of the new data on paragonite with those of a K-muscovite and a Na-rich muscovite (Comodi and Zanazzi 1995) revealed a clear trend toward decreasing of compressibility when Na substitutes for K atoms in the interlayer sites.     

The Thermodynamics of Ordered Perovskites on the CaTiO3-FeTiO3 Join

The enthalpies of formation from ilmenite, FeTiO₃, and perovskite, CaTiO₃, of two intermediate ordered perovskite phases, CaFeTi₂O₆ and CaFe₃Ti₄O₁₂, have been measured at 801°C using oxide melt solution calorimetry. These data, in combination with experiments at high pressure and temperature, indicate that below 1518±50°C CaFe₃Ti₄O₁₂ is stable at lower pressures (∼7 GPa at 1200°C) than CaFeTi₂O₆ (∼13 GPa at 1200°C). This relationship should be reversed, and CaFeTi₂O₆ should become stable at lower pressures than CaFe₃Ti₄O₁₂, at temperatures above 1518±50°C. These intermediate phases are of petrological interest because they form as a reaction between two minerals, ilmenite and perovskite, which are commonly associated in kimberlites, and because their pressure-temperature range of formation overlaps that of origin of kimberlites. Received: 10 November 1997 / Revised; accepted: 15 January 1998