Equations of state of Plagioclase Feldspars

The volume variation with pressure of seven intermediate plagioclase feldspars has been determined by high-pressure single-crystal X-ray diffraction. The bulk moduli of plagioclases for a 3rd-order Birch-Murnaghan EoS can be described by the following pair of equations:

The magnesiowüstite: iron equilibrium and its implications for the activity-composition relations of (Mg,Fe)<Subscript>2</Subscript>SiO<Subscript>4</Subscript> olivine solid solutions

The chemical potential of oxygen (µO₂) in equilibrium with magnesiowüstite solid solution (Mg, Fe)O and metallic Fe has been determined by gas-mixing experiments at 1,473 K supplemented by solid-cell EMF experiments at lower temperatures. The results give:

Crystal chemical controls on rare earth element partitioning between epidote-group minerals and melts: an experimental and theoretical study

We have experimentally determined the partitioning of REE (rare earth elements) between zoisite and hydrous silicate melt at 1,100 °C and 3 GPa. All REE behave moderately compatible in zoisite with respect to the melt and all show a smooth parabolic dependence on ionic radius. The partitioning parabola peaks at Nd , and the compatibility slightly decreases towards La and decreases by half an order of magnitude towards Yb . Application of the elastic strain model of Blundy and Wood (1994) to the available zoisite and allanite REE mineral/melt partitioning data and comparison with partitioning pattern calculated from a combination of structural and physical data (taken from the literature) with the elastic strain model suggest that in zoisite REE prefer the A1-site and that only La and Ce are incorporated into the A2-site in significant amounts. In contrast, in allanite, all REE are preferentially incorporated into the large and highly co-ordinated A2 site. As a result, zoisite fractionates the MREE effectively from the HREE and moderately from the LREE, while allanite fractionates the LREE very effectively from the MREE and HREE. Consequently, the presence of either zoisite or allanite during slab melting will lead to quite different REE pattern in the produced melt.Editorial responsibility: J. Hoefs     

Impact of land use change on groundwater quality in a typical karst watershed of southwest China: a case study of the Xiaojiang watershed,Yunnan Province

The impact of land-use change on the quality of groundwater in the Xiaotjiang watershed, China was assessed for the period 1982–2004. Groundwater samples were collected from 30 monitoring points across the watershed, and were representative of the various changes, determined by remote sensing and geographical information systems. The results indicate that 610 km² (60% of the total watershed area) were subject to land-use change during the period. The most important changes were the conversion of 135 km² of forested land to cultivated land, and 211 km² of unused land to cultivated land. The main impact was ascribed to diffuse pollution from fertilizers applied to newly cultivated land, and from building development. Overall the groundwater pH value was significantly increased, as were the concentrations of ions , , , , and Cl⁻ in groundwater whilst the concentrations of Ca²⁺ and declined. More precisely, in the regions where forested land and unused land were converted into cultivated land, the pH value and the concentrations of Mg²⁺, , , , , Cl⁻ increased whilst the concentrations of Ca²⁺ and declined. However in the region where cultivated land was converted into construction land, the pH value and the concentrations of Ca²⁺, Mg²⁺, , , , , , Cl⁻ increased.

---

Résumé  L’impact des changements de l’utilisation du territoire sur la qualité de l’eau souterraine dans le bassin versant de Xiaojiang, en Chine, a été évalué de 1982 à 2004. Des échantillons d’eau souterraine ont été récoltés à partir de 30 points d’observation éparpillés sur le bassin, représentant les divers changements déterminés par télédétection et système d’information géographique. Les résultats indiquent que 610 km² (soit 60% de la surface du bassin) ont été sujets à des modifications de l’utilisation du territoire sur cette période. Les changements les plus importants furent la conversion de 135 km² de forêt et 211 km² de terres inutilisées en terres cultivées. Le principal impact est attribué à la pollution diffuse des engrais utilisés en agriculture et pour les batiments. De manière générale le pH de l’eau souterraine a augmenté significativement, ainsi que les concentrations des ions , , , , et Cl⁻, tandis que les concentration en Ca²⁺ et ont diminué. Plus précisément dans les régions transformées en terres cultivées, la valeur du pH et les concentrations en Mg²⁺, , , , , Cl⁻ ont augmenté tandis que les concentrations en Ca²⁺ et ont diminué. Toutefois dans les régions cultivées converties en zones de construction, le pH et les concentrations en Ca²⁺, Mg²⁺, , , , , , Cl⁻ ont augmenté.

---

Resumen  El impacto del cambio en uso de la tierra en la calidad del agua en la cuenca Xiaojiang, China fue evaluado para el periodo 1982–2004. Muestras de agua subterránea fueron tomadas de 30 puntos de monitoreo a través de la cuenca, y fueron representativas de los múltiples cambios, determinados por sensores remotos y sistemas de información geográfica. Los resultados indican que 610 km² (60% del área total de la cuenca) estaban sujetos a cambios de uso de la tierra durante el periodo estudiado. Los cambios más importantes fueron la conversión de 135 km² de bosques a tierra cultivada, y 211 km² de tierra sin uso (ociosa) a tierra cultivada. El impacto principal fue causado por contaminación difusa de fertilizantes aplicados a la tierra recientemente cultivada, y a desarrollo de construcciones. En general el pH en agua subterránea creció significantemente, al igual que las concentraciones de los iones , , , , y Cl⁻ en agua subterránea mientras que las concentraciones de Ca²⁺ y decrecieron. Mas precisamente, en las regiones donde bosque y tierra ociosa fueron convertidas en tierra cultivada, el valor de pH y las concentraciones de Mg²⁺, , , , , Cl⁻ crecieron mientras las concentraciones de Ca²⁺ y decrecieron. Sin embargo en la región donde tierra cultivada fue convertida en construcciones, el valor de pH y las concentraciones de Ca²⁺, Mg²⁺, , , , , , Cl⁻ crecieron.

     

Nature and profile distribution of iron and aluminum in relation to pedogenic processes in some soils developed under tropical environment in India

Different extractants were used to determine various forms of Fe and Al in three soil profiles developed under a tropical environment in India. The average contents of Fe and Al, extracted by different extracting reagents were found to be in descending order, as follows:

Systematic study of liquidas phase relations in olivine melilitite +H2O +CO2 at high pressures and petrogenesis of an olivine melilitite magma

Near-liquidus phase relationships of a spinel lherzolite-bearing olivine melilitite from Tasmania were investigated over a P, T range with varying , , and . At 30 kb under MH-buffered conditions, systematic changes of liquidus phases occur with increasing ( = CO₂/CO₂ +H₂O+olivine melilitite). Olivine is the liquidus phase in the presence of H₂O alone and is joined by clinopyroxene at low . Increasing eliminates olivine and clinopyroxene becomes the only liquidus phase. Further addition of CO₂ brings garnet+orthopyroxene onto the liquidus together with clinopyroxene, which disappears with even higher CO₂. The same systematic changes appear to hold at higher and lower pressures also, only that the phase boundaries are shifted to different . The field with olivine- +clinopyroxene becomes stable to higher with lower pressure and approaches most closely the field with garnet+orthopyroxene+clinopyroxene at about 27 kb, 1160 °C, 0.08 and 0.2 (i.e., 6–7% CO₂+ 7–8% H₂O). Olivine does not coexist with garnet+orthopyroxene+clinopyroxene under these MH-buffered conditions. Lower oxygen fugacities do not increase the stability of olivine to higher and do not change the phase relationships and liquidus temperatures drastically. Thus, it is inferred that olivine melilitite 2927 originates as a 5% melt (inferred from K_{2 O} and P₂O₅ content) from a pyrolite source at about 27kb, 1160 dg with about 6–7% CO₂ and 7–8% H₂O dissolved in the melt. The highly undersaturated character of the melt and the inability to find olivine together with garnet and orthopyroxene on the liquidus (in spite of the close approach of the respective liquidus fields) can be explained by reaction relationships of olivine and clinopyroxene with orthopyroxene, garnet and melt in the presence of CO₂.     

The Skaergaard liquid line of descent revisited

There is a fundamental conflict between the suggestion that the iron content of Skaergaard liquids increases during Fe–Ti oxide fractionation and the observation that at the same time oxygen fugacity () drops by two log-units below the fayalite-magnetite-quartz oxygen buffer (FMQ). A new petrographic study of average Skaergaard gabbros shows that the total modal content of Fe–Ti oxides is about 22% in the early LZc and markedly decreases to below 5% in the UZc. Forward modeling based on these modal constraints, as well as experimental results on Skaergaard-related dikes, predicts that fractionation of troctolitic LZa gabbros drives the derivative liquid towards a high-iron content. Strong iron enrichment continues, together with a small decline in silica, during LZb crystallization due to the appearance of augite as a fractionating phase. The fractionation of Fe–Ti oxides in the LZc initially suppresses iron enrichment and reverses the silica trend to one of slight enrichment. However, continued evolution into the UZ produces liquids with maximum UZc FeO* content of 23–25 wt.% and SiO₂ content of 53 wt.% (FeO* is total iron as FeO). The maximum in FeO* is dependent on several factors of which the Fe–Ti oxide mode has the strongest effect. The during crystallization of the LZc is widely thought to have been at, or slightly below, the fayalite-magnetite-quartz oxygen buffer (FMQ). Under closed system evolution, incorporation of ferric iron into augite during formation of the LZb restricts the increase in to about 0.1 log-units above FMQ (=0.1 ΔFMQ). Likewise, crystallization of the LZc through the UZa, involving Fe–Ti oxide minerals, leads to a decline in of less than 0.1 ΔFMQ. Crystallization of the UZb-c gabbros results in oxidation to a maximum of 0.5 ΔFMQ. This behavior can account for the iron-rich character of the UZ gabbros, as well as, the low modal content of Fe–Ti oxides. Thus, evolved Skaergaard liquids are high in iron and contain a modest amount of SiO₂. Our modeling result do not account for a strong drop in through the layered series. Such a drop would require an unacceptably high proportion of Fe–Ti oxides and high-magnetite content in the fractionating assemblage. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

New thermodynamic models and revised calibrations for the Ti-in-zircon and Zr-in-rutile thermometers

The models recognize that ZrSiO₄, ZrTiO₄, and TiSiO₄, but not ZrO₂ or TiO₂, are independently variable phase components in zircon. Accordingly, the equilibrium controlling the Zr content of rutile coexisting with zircon is ZrSiO₄ = ZrO₂ (in rutile) + SiO₂. The equilibrium controlling the Ti content of zircon is either ZrSiO₄ + TiO₂ = ZrTiO₄ + SiO₂ or TiO₂ + SiO₂ = TiSiO₄, depending whether Ti substitutes for Si or Zr. The Zr content of rutile thus depends on the activity of SiO₂ as well as T, and the Ti content of zircon depends on and as well as T. New and published experimental data confirm the predicted increase in the Zr content of rutile with decreasing and unequivocally demonstrate that the Ti content of zircon increases with decreasing . The substitution of Ti in zircon therefore is primarily for Si. Assuming a constant effect of P, unit and that and are proportional to ppm Zr in rutile and ppm Ti in zircon, [log(ppm Zr-in-rutile) + log] = A₁ + B₁/T(K) and [log(ppm Ti-in-zircon) + log − log] = A₂ + B₂/T, where the A and B are constants. The constants were derived from published and new data from experiments with buffered by either quartz or zircon + zirconia, from experiments with defined by the Zr content of rutile, and from well-characterized natural samples. Results are A₁ = 7.420 ± 0.105; B₁ = −4,530 ± 111; A₂ = 5.711 ± 0.072; B₂ = −4,800 ± 86 with activity referenced to α-quartz and rutile at P and T of interest. The zircon thermometer may now be applied to rocks without quartz and/or rutile, and the rutile thermometer applied to rocks without quartz, provided that and are estimated. Maximum uncertainties introduced to zircon and rutile thermometry by unconstrained and can be quantitatively assessed and are ≈60 to 70°C at 750°C. A preliminary assessment of the dependence of the two thermometers on P predicts that an uncertainty of ±1 GPa introduces an additional uncertainty at 750°C of ≈50°C for the Ti-in-zircon thermometer and of ≈70 to 80°C for the Zr-in-rutile thermometer.     

Solubility of Ca and Al in orthopyroxene from spinel peridotite: an improved version of an empirical geothermometer

Geothermometric equations for spinel peridotites by Fujii (1976), Gasparik and Newton (1984), and Chatterjee, and Terhart (1985) based on the reaction enstatite (en)+spinel (sp)Mg–Tschermaks (mats)+forsterite (fo) were tested using a nearly isothermal suite of mantle xenoliths from the Eifel, West Germany. In spite of using activities of MgAl₂O₄, en, and mats to allow for the non-ideal solution behaviour of the constituent phases, temperatures calculated from these equations systematically change as a function of Cr/(Cr+AL+Fe³⁺) in spinel. We propose an improved version of the empirical geothermometer for spinel peridotites of Sachtleben and Seck (1981) derived from the evaluation of the solubilities of Ca and Al in orthopyroxene from more than 100 spinel peridotites from the Rhenish Volcanic Province. A least squares regression yielded a smooth correlation between

Fluorite stability in silicic magmas

Recent experimental evidence is used to assess the conditions under which fluorite forms an early crystallising phase in silicic magmas. Fluorite solubility primarily depends on the (Na + K)/Al balance in the coexisting silicic melt, reaching a minimum in metaluminous melts. It can display reaction relationships with topaz and titanite, depending on changes in melt composition during crystallisation. An empirical model of fluorite stability in Ca-poor peralkaline rhyolite melts is derived and applied to selected rocks:

Water in microcrystalline quartz of volcanic origin: Agates

Agates of volcanic origin, containing the different quartz species, fibrous, length-fast chalcedony (CH), granular fine quartz (FQ), and fibrous, length-slow, to lepidospheric quartzine (QN), have been investigated to evaluate possible relations between microstructure, i.e. crystallite size and texture, refractive indices, densities, contents of trace elements and of water, as well as dehydration behaviour. By means of near infrared spectroscopy, total water contents , could be differentiated quantitatively into contents of molecular water, , and silanole-group water, . Despite the low total water contents of the agates studied ( between 1 and 2 wt.%), near infrared spectroscopy results in reliable data on and .Wall-layering CH consists of fibrous quartz crystals and exhibits higher C-ratios, , than horizontally layered FQ which consists predominantly of granular quartz crystals (C _CH=0.45±0.11 (N=6), C _FQ=0.36±0.10 (N=4). This result is interpreted to be due to analogy with the behaviour of C-ratios in fluid phase-deposited opals-AN (hyalithe) and liquid phase-deposited opals-AG (non-crystalline opal) or -CT (common opal) (Langer and Flörke 1974).Translucent layers of CH show mostly lower refractive indices, when measured parallel than when measured perpendicular to the axes of the quartz fibers. The same is true for milky layers of CH. Crystallite sizes are smaller in the latter than in the former.For all samples studied, exists a positive correlation between at% (1/2Ca+1/2Mg+Na+K+Li) and at% (Al³⁺+Fe³⁺). This indicates that at least parts of (A1³⁺+ Fe³⁺) substitute for Si in the quartz structure. The charge is balanced by incorporation of di- and mono-valent cations in structural interstices. When the quantity at % H⁺, as obtained from , is included into the sum at% (1/2 Me²⁺+Me⁺), the above correlation is destroyed. This result could be indicative for a strong concentration of the Si-OH groups in the surface of the quartz microcrystallites.     

An experimental study of the diffusion of oxygen in quartz and albite using an overgrowth technique

Diffusion rates of¹⁸O tracer in quartz ( c, 1 Kb H₂O) and Amelia albite ( 001, 2 Kb H₂O) have been measured, using Secondary Ion Mass Spectrometry (SIMS). A new technique involving hydrothermal deposition of labelled materials has removed the possibility of pressure solution-reprecipitation processes adversely affecting the experiments. Reported diffusion constants are:-quartz ( c), ,Q=98±7 KJ mol^–1 (600–825° C, 1 Kb); Amelia albite ( 001), ,Q=85±7 KJ mol^–1, (400–600° C, 2 Kb). Measured quartz¹⁸O diffusivities decrease discontinuously at the- transition, reflecting strong structural influences. The reported albite data agree with previously recorded studies, but-quartz data indicate significantly lower activation energies. Possible causes of this discrepancy, and some geological consequences, are noted.     

The Influence of Redox Reactions on the Uptake of Dissolved Ce by Suspended Fe and Mn Oxide Particles

Laboratory experiments were conducted to evaluate the partitioning ofrare earth elements (REE) between solution and suspended particles. Becauseof their strong tendency to complex, the REE can be used to study a varietyof marine processes and in particular particle scavenging. In this study, anemphasis was placed on examining abiotic redox processes that influence theuptake of dissolved Ce by particles. Batch sorption experiments wereconducted with REE and synthetic mineral phases over the range of pH4–9. The solutions varied in ionic strength between 0 and 0.7 M andconsisted of individual solutes (NaNO₃, NaCl, andNa₂SO₄), ionic mixtures that duplicate theseawater composition, and natural seawater. The uptake of REE from solutionwas also studied at a Pt electrode coated with using cyclic voltametry. Experimental results are consistent with uptake of dissolved Ce onto occurring by a combination of oxidativescavenging and surface complexation. The contribution of oxidativescavenging to the removal of Ce from solution is most pronounced at acidicpH, where the strictly trivalent REE exhibit little propensity for sorptiononto . Sorption of dissolved Ce onto FeOOH occursin a manner analogous to that of the other strictly trivalent REE and nocontribution from oxidative scavenging is observed on this mineral phase atlow pH. Our work also substantiates the hypothesis that anions in solution,particularly and Cl^-, aswell as those adsorbed on the surface of the particles, influence the extentof Ce uptake by . This revised version was published online in July 2006 with corrections to the Cover Date.     

Fe2⇄Mg2 and TiAl2⇄MgSi2 exchange reactions between clinopyroxenes and silicate melts

The interdependence of the Fe(Mg)_–1 (e.g., FeO-MgO in silicate melt; CaFeSi₂O₆-CaMgSi₂O₆ in pyroxene) and TiAl₂(MgSi₂)_–1 exchange reactions between silicate melts and coexisting Ca-pyroxene has been examined. High-calcium clinopyroxenes were grown in 1 atmosphere melting and crystallization experiments on rock powders spanning the composition range tholeiite to melilitite (1,0922+Mg²⁺ exchange and suggest that at given values of extent of Fe(Mg)_–1 substitution is strongly coupled with the TiAl₂(MgSi₂)_–1 substitution in pyroxenes near the five-component space CaMg(Si₂O₆-CaFe(Si)₂O₆-CaTi(Al)₂O₆-CaFe(Al,Si)₂O₆-CaAl(Al,Si)₂O₆. The inferred stabilization of Ti in iron-rich relative to magnesium pyroxene is consistent with the operation of Fe²⁺Ti⁴⁺ intervalence charge transfer interactions (e.g., Rossman 1980) and observations on zoning in natural titanaugites (e.g., Tracy and Robinson 1977). Although the rims of some pyroxenes grown in some melting experiments exhibit prominent zoning in TiAl₂(MgSi₂)_–1, the average values of inferred from the compositions of these pyroxenes, together with those of the relatively homogeneous pyroxenes produced in crystallization experiments, exhibit a 11 correlation with values of derived from the solution model of Ghiorso et al. (1983) with a standard error of 750 calories. The Ti contents of Ca-rich pyroxenes crystallizing from a wide range of natural silicate liquids can therefore be predicted.     

The generation and crystallization conditions of the Proterozoic Harney Peak Leucogranite,Black Hills,South Dakota,USA: Petrologic and geochemical constraints

The mineralogy, petrology and geochemistry of the Proterozoic Harney Peak Granite, Black Hills, South Dakota, were examined in view of experimentally determined phase equilibria applicable to granitic systems in order to place constraints on the progenesis of peraluminous leucogranites and commonly associated rare-element pegmatites. The granite was emplaced at 3–4 kbar as multiple sills and dikes into quartz-mica schists at the culmination of a regional high-temperature, low-pressure metamorphic event. Principally along the periphery of the main pluton and in satellite intrusions, the sills segregated into granite-pegmatite couplets. The major minerals include quartz, K-feldspar, sodic plagioclase and muscovite. Biotite-{Mg No. [Molar MgO/(MgO+FeO)]=0.32-0.38} is the predominant ferromagnesian mineral in the granite's core, whereas at the periphery of the main pluton and in the satellite intrusions tourmaline (Mg No.=0.18–0.48) is the dominant ferromagnesian phase. Almandine-spessartine garnet is also found in the outer intrusions. There is virtually a complete overlap in the wide concentration ranges of SiO₂, CaO, MgO, FeO, Sr, Zr, W of the biotite- and tourmaline-bearing granite suites with no discernable differentiation trends on Harker diagrams, precluding the derivation of one suite from the other by differentiation following emplacement. This is consistent with the oxygen isotope compositions which are 11.5 ± 0.6 for the biotite granites and 13.2 ± 0.8 for the tourmaline granites, suggesting derivation from different sources. The concentrations of TiO₂ and possibly Ba are higher and of MnO and B are lower in the biotite granites. The normative Orthoclase/Albite ratio is extremely variable ranging from 0.26 to 1.65 in the biotite granites to 0.01–1.75 in the tourmaline granites. Very few sample compositions fall near the high-pressure, watersaturated haplogranite minima-eutectic trend, indicating that the granites for the most part are not minimum melts generated under conditions with =1. Instead, most biotite granites are more potassic than the water-saturated minima and eutectics and in analogy with experimentally produced granitic melts, they are best explained by melting at 6 kbar, <1 and temperatures 800°C. Such high temperatures are also indicated by oxygen isotope equilibration among the constituent minerals (Nabelek et al. 1992). Several of the tourmaline granite samples contain virtually no K-feldspar and have oxygen isotope equilibration temperatures 716–775°C. Therefore, they must represent high-temperature accumulations of liquidus minerals crystallized under equilibrium conditions from melts more sodic than the water-saturated haplogranite minima or during fractionation of intruded melts into granite-pegmatite couplets accompanied by volatile-aided differentiation of the alkali elements. The indicated high temperatures, <1, the relatively high TiO₂ and Ba concentrations and the relatively low values of the biotite granites suggest that they were generated by high-extent, biotite-dehydration melting of an immature Archean metasedimentary source. The ascent of the hot melts may have triggered low-extent, muscovite-dehydration melting of schists higher in the crust producing the high-B, low-Ti melts comprising the periphery of the main pluton and the satellite intrusions. Alternatively, the different granite types may be the result of melting of a vertical section of the crust in response to the ascent of a thermal pulse, with the low- biotite granites generated at a deeper, hotter region and the high- tourmaline granites at a higher, cooler region of the crust. The low-Ti and high-B concentrations in the high- melts resulted in the crystallization of tourmaline rather than biotite, which promoted the observed differentiation of the melts into the granitic and pegmatitic layers found along the periphery of the main pluton and the satellite intrusions.     

A regional gradient in the composition of metamorphic fluids in pelitic schist,Pecos Baldy,New Mexico

Two metamorphic isograds cut across graphitic schist near Pecos Baldy, New Mexico. The southern isograd marks the first coexistence of staurolite with biotite, whereas the northern isograd marks the first coexistence of andalusite with biotite. The isograds do not record changes in temperature or pressure. Instead, they record a regional gradient in the composition of the metamorphic fluid phase. Ortega Quartzite, which contains primary hematite, lies immediately north of the graphitic schist. Mineral compositions within the schist change gradually toward the quartzite, reflecting gradients in and . The chemical potential gradients, locally as high as 72 cal/m in and 9 cal/m in , controlled the positions of the two mapped isograds. The staurolite-biotite isograd records where fell below 0.80, at near 10^–23 bars; the andalusite-biotite isograd records where fell below 0.25, at near 10^–22 bars. Dehydration and oxidation were coupled by graphite-fluid equilibrium.The chemical potential gradients apparently formed during metamorphism, as graphite in schist reacted with hematite in quartzite. Local oxidation of graphite formed CO₂ which triggered dehydration reactions along the schistquartzite contact. This process created a C-O-H fluid which infiltrated into overlying rocks. Upward infiltration, local fluid-rock equilibration and additional infiltration proceeded until the composition of the infiltrating fluid evolved to that in equilibrium with the infiltrated rock. This point occurs very close to the staurolite-biotite isograd. Pelitic rocks structurally above this isograd show no petrographic evidence of infiltration, even though calculations indicate that volumetric fluid/rock ratios may have exceeded 15 and thin, rare calc-silicate beds show extensive K-metasomatism and quartz veining.     

H<Subscript>2</Subscript>O diffusion in peralkaline to peraluminous rhyolitic melts

The diffusion of water in a peralkaline and a peraluminous rhyolitic melt was investigated at temperatures of 714–1,493 K and pressures of 100 and 500 MPa. At temperatures below 923 K dehydration experiments were performed on glasses containing about 2 wt% H₂O _t in cold seal pressure vessels. At high temperatures diffusion couples of water-poor (<0.5 wt% H₂O _t ) and water-rich (~2 wt% H₂O _t ) melts were run in an internally heated gas pressure vessel. Argon was the pressure medium in both cases. Concentration profiles of hydrous species (OH groups and H₂O molecules) were measured along the diffusion direction using near-infrared (NIR) microspectroscopy. The bulk water diffusivity () was derived from profiles of total water () using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between and Both methods consistently indicate that is proportional to in this range of water contents for both bulk compositions, in agreement with previous work on metaluminous rhyolite. The water diffusivity in the peraluminous melts agrees very well with data for metaluminous rhyolites implying that an excess of Al₂O₃ with respect to alkalis does not affect water diffusion. On the other hand, water diffusion is faster by roughly a factor of two in the peralkaline melt compared to the metaluminous melt. The following expression for the water diffusivity in the peralkaline rhyolite as a function of temperature and pressure was obtained by least-squares fitting:

Mechanisms of hydrogen incorporation and diffusion in iron-bearing olivine

The incorporation and diffusion of hydrogen in San Carlos olivine (Fo₉₀) single crystals were studied by performing experiments under hydrothermal conditions. The experiments were carried out either at 1.5 GPa, 1,000°C for 1.5 h in a piston cylinder apparatus or at 0.2 GPa, 900°C for 1 or 20 h in a cold-seal vessel. The oxygen fugacity was buffered using Ni–NiO, and the silica activity was buffered by adding San Carlos orthopyroxene powders. Polarized Fourier transform infrared (FTIR) spectroscopy was utilized to quantify the hydroxyl distributions in the samples after the experiments. The resulting infrared spectra reproduce the features of FTIR spectra that are observed in olivine from common mantle peridotite xenoliths. The hydrogen concentration at the edges of the hydrogenated olivine crystals corresponds to concentration levels calculated from published water solubility laws. Hydrogen diffusivities were determined for the three crystallographic axes from profiles of water content as a function of position. The chemical diffusion coefficients are comparable to those previously reported for natural iron-bearing olivine. At high temperature, hydrogenation is dominated by coupled diffusion of protons and octahedrally coordinated metal vacancies where the vacancy diffusion rate limits the process. From the experimental data, we determined the following diffusion laws (diffusivity in m² s⁻¹, activation energies in kJ mol⁻¹): for diffusion along [100] and [010]; for diffusion along [001]. These diffusion rates are fast enough to modify significantly water contents within olivine grains in xenoliths ascending from the mantle.     

APL computer programs for thermodynamic calculations of equilibria in P-T-X CO 2 space

APL computer programs for the thermodynamic calculation of devolatilization and solid-solid equilibria operate using stored values for the molar volume and entropy of solids, the free energies of H₂O and CO₂, and the free energies of formation for 110 geologically-important phases. P-T-X _{CO 2} calculations of devolatilization equilibria can be made at pressures from 0.2 through 10 kb, and temperatures from 200 through 1,000° C. P-T-X calculations of solid-solid equilibria may be accomplished at pressures to 30 kb and temperatures to 1,000° C. Calculations can be extrapolations from experimental points, or direct calculations from thermochemical data alone. Options are available in these programs to consider effects of: real vs. ideal gas mixing, thermal expansion and compressibility, solid solution, fluid pressure differing from solid pressure, and uncertainties in high-temperature entropies.A collection of thermodynamic data programs accompanies the programs for calculating P-T-X _{CO 2} equilibria. Over a wide range of physical conditions, the data functions report free energies, entropies, fugacities of H₂O and CO₂, high temperature entropies of solids, and activities of components in H₂O-CO₂ mixtures.List of Symbols Activity of H₂O and CO₂ - G_f Free energy of formation of a phase from elements - G_r Free energy change of reaction - G _r ^o Standard state free energy change of a reaction - Free energies of pure H₂O and CO₂ - H _r ^o Standard state enthalpy change for a reaction - K Equilibrium constant - R Gas constant - S _r ^o Standard state entropy change of reaction - S _s ^o Standard state entropy change of solids in a reaction - V_s ^o Standard state volume change of a reaction - V_s ^o Standard state volume change of solids in a reaction - Mole fraction of H₂O and CO₂ - Activity coefficient of H₂O and CO₂