Amphibole composition in tonalite as a function of pressure: an experimental calibration of the Al-in-hornblende barometer

The Al-in-hornblende barometer, which correlates Al^tot content of magmatic hornblende linearly with crystallization pressure of intrusion (Hammarstrom and Zen 1986), has been calibrated experimentally under water-saturated conditions at pressures of 2.5–13 kbar and temperatures of 700–655°C. Equilibration of the assemblage hornlende-biotite-plagioclase-orthoclasequartz-sphene-Fe-Ti-oxide-melt-vapor from a natural tonalite 15–20° above its wet solidus results in hornblende compositions which can be fit by the equation: P(±0.6 kbar) = –3.01 + 4.76 Al _hbl ^tot r ²=0.99, where Al^tot is the total Al content of hornblende in atoms per formula unit (apfu). Al^tot increase with pressure can be ascribed mainly to a tschermak-exchange ( ) accompanied by minor plagioclase-substitution ( ). This experimental calibration agrees well with empirical field calibrations, wherein pressures are estimated by contact-aureole barometry, confirming that contact-aureole pressures and pressures calculated by the Al-in-hornblende barometer are essentially identical. This calibration is also consistent with the previous experimental calibration by Johnson and Rutherford (1989b) which was accomplished at higher temperatures, stabilizing the required buffer assemblage by use of mixed H₂O-CO₂ fluids. The latter calibration yields higher Al^tot content in hornblendes at corresponding pressures, this can be ascribed to increased edenite-exchange ( ) at elevated temperatures. The comparison of both experimental calibrations shows the important influence of the fluid composition, which affects the solidus temperature, on equilibration of hornblende in the buffering phase assemblage.     

Reaction spaces and P-T paths: from amphibole eclogite to greenschist facies in the Austroalpine domain (Oetztal Complex)

The transition from feldspar amphibolite to eclogite is a very wide P-T field that extends from some-where close to 5 kbar where the garnet-amphibole pair starts to appear, to 10–20 kbar at albite-out reaction, then up to 25–30 kbar where an hydrated phase such as amphibole can be stable with pyroxene and garnet. Thus the assemblage garnet (py)+ amphibole (tr)+epidote (cz)±plagioclase (ab)±clinopyroxene (di)±quartz (qz)±fluid is commonly reported in a large number of metamorphic terrains. These mineral phases are complex solid-solutions which adapt to variations in environmental conditions mainly by means of continuous reactions. The reaction space, introduced by. Thompson in 1982a, provides a very elegant and powerful tool to approach these high-variance assemblages. The reactions:

Experimental calibration of Forsterite–Anorthite–Ca-Tschermak–Enstatite (FACE) geobarometer for mantle peridotites

The crystallization of plagioclase-bearing assemblages in mantle rocks is witness of mantle exhumation at shallow depth. Previous experimental works on peridotites have found systematic compositional variations in coexisting minerals at decreasing pressure within the plagioclase stability field. In this experimental study we present new constraints on the stability of plagioclase as a function of different Na₂O/CaO bulk ratios, and we present a new geobarometer for mantle rocks. Experiments have been performed in a single-stage piston cylinder at 5–10 kbar, 1050–1150?°C at nominally anhydrous conditions using seeded gels of peridotite compositions (Na₂O/CaO?=?0.08–0.13; X _Cr = Cr/(Cr?+?Al)?=?0.07–0.10) as starting materials. As expected, the increase of the bulk Na₂O/CaO ratio extends the plagioclase stability to higher pressure; in the studied high-Na fertile lherzolite (HNa-FLZ), the plagioclase-spinel transition occurs at 1100?°C between 9 and 10 kbar; in a fertile lherzolite (FLZ) with Na₂O/CaO?=?0.08, it occurs between 8 and 9 kbar at 1100?°C. This study provides, together with previous experimental results, a consistent database, covering a wide range of P–T conditions (3–9 kbar, 1000–1150?°C) and variable bulk compositions to be used to define and calibrate a geobarometer for plagioclase-bearing mantle rocks. The pressure sensitive equilibrium:

$$\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{4}}}^{{\text{Ol}}}}\limits_{{\text{Forsterite}}} +\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{8}}}^{{\text{Pl}}}}\limits_{{\text{Anorthite}}~} =\mathop {{\text{CaA}}{{\text{l}}_{\text{2}}}{\text{Si}}{{\text{O}}_{\text{6}}}^{{\text{Cpx}}}}\limits_{{\text{Ca-Tschermak}}} +{\text{ }}\mathop {{\text{M}}{{\text{g}}_{\text{2}}}{\text{S}}{{\text{i}}_{\text{2}}}{{\text{O}}_{\text{6}}}^{{\text{Opx}}}}\limits_{{\text{Enstatite}}} ,$$

has been empirically calibrated by least squares regression analysis of experimental data combined with Monte Carlo simulation. The result of the fit gives the following equation:

$$P=7.2( \pm 2.9)+0.0078( \pm 0.0021)T{\text{ }}+0.0022( \pm 0.0001)T{\text{ }}\ln K,$$

$${R^2}=0.93,$$

where P is expressed in kbar and T in kelvin. K is the equilibrium constant K?=?a _CaTs × a _en/a _an × a _fo, where a _CaTs, a _en, a _an and a _fo are the activities of Ca-Tschermak in clinopyroxene, enstatite in orthopyroxene, anorthite in plagioclase and forsterite in olivine. The proposed geobarometer for plagioclase peridotites, coupled to detailed microstructural and mineral chemistry investigations, represents a valuable tool to track the exhumation of the lithospheric mantle at extensional environments.

     

Stabilitätsbeziehungen zwischen Chlorit,Cordierit und Almandin bei der Metamorphose

The stability relations between cordierite and almandite in rocks, having a composition of CaO poor argillaceous rocks, were experimentally investigated. The starting material consisted of a mixture of chlorite, muscovite, and quartz. Systems with widely varying Fe²⁺/Fe²⁺+Mg ratios were investigated by using two different chlorites, thuringite or ripidolite, in the starting mixture. Cordierite is formed according to the following reaction: $${\text{Chlorite + muscovite + quartz}} \rightleftharpoons {\text{cordierite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}}$$ . At low pressures this reaction characterizes the facies boundary between the albite-epidotehornfels facies and the hornblende-hornfels facies, at medium pressures the beginning of the cordierite-amphibolite facies. Experiments were carried out reversibly and gave the following equilibrium data: 505±10°C at 500 bars H₂O pressure, 513±10°C at 1000 bars H₂O pressure, 527±10°C at 2000 bars H₂O pressure, and 557±10°C at 4000 bars H₂O pressure. These equilibrium data are valid for the Fe-rich starting material, using thuringite as the chlorite, as well as for the Mg-rich starting mixture with ripidolite. At 6000 bars the equilibrium temperature for the Mg-rich mixture is 587±10°C. In the Fe-rich mixture almandite was formed instead of cordierite at 6000 bars. The following reaction was observed: $${\text{Thuringite + muscovite + quartz}} \rightleftharpoons {\text{almandite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + H}}_{\text{2}} {\text{O}}$$ . Experiments with the Fe-rich mixture, containing Fe²⁺/Fe²⁺+Mg in the ratio 8∶10, yielded three stability fields in a P,T-diagram (Fig.1):

1. Above 600°C/5.25 kb and 700°C/6.5 kb almandite+biotite+Al₂SiO₅ coexist stably, cordierite being unstable.
2. The field, in which almandite, biotite and Al₂SiO₅ are stable together with cordierite, is restricted by two curves, passing through the following points:
   1. 625°C/5.5 kb and 700°C/6.5 kb,
   2. 625°C/5.5 kb and 700°C/4.0 kb.
3. At conditions below curves 1 and 2b, cordierite, biotite, and Al₂SiO₅ are formed, but no garnet.

An appreciable MnO-content in the system lowers the pressures needed for the formation of almandite garnet, but the quantitative influence of the spessartite-component on the formation of almandite could not yet be determined. the Mg-rich system with Fe²⁺/Fe²⁺+Mg=0.4 garnet did not form at pressures up to 7 kb in the temperature range investigated. Experiments at unspecified higher pressures (in a simple squeezer-type apparatus) yielded the reaction: $${\text{Ripidolite + muscovite + quartz}} \rightleftharpoons {\text{almandite + biotite + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + H}}_{\text{2}} {\text{O}}$$ . Further experiments are needed to determine the equilibrium data. The occurence of garnet in metamorphic rocks is discussed in the light of the experimental results.     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

A note on the significance of the assemblage calcite-quartz-plagioclase-paragonite-graphite

The biotite zone assemblage: calcite-quartz-plagioclase (An₂₅)-phengite-paragonite-chlorite-graphite, is developed at the contact between a carbonate and a pelite from British Columbia. Thermochemical data for the equilibrium paragonite+calcite+2 quartz=albite+ anorthite+CO₂+H₂O yields: $$\log f{\text{H}}_{\text{2}} {\text{O}} + \log f{\text{CO}}_{\text{2}} = 5.76 + 0.117 \times 10^{ - 3} (P - 1)$$ for a temperature of 700°K and a plagioclase composition of An₂₅. By combining this equation with equations describing equilibria between graphite and gas species in the system C-H-O, the following partial pressures: \(P{\text{H}}_2 {\text{O}} = 2572{\text{b, }}P{\text{CO}}_2 = 3162{\text{b, }}P{\text{H}}_2 = 2.5{\text{b, }}P{\text{CH}}_4 = 52.5{\text{b, }}P{\text{CO}} = 11.0{\text{b}}\) are obtained for \(f{\text{O}}_2 = 10^{ - 26}\) . If total pressure equals fluid pressure, then the total pressure during metamorphism was approximately 6 kb. The total fluid pressure calculated is extremely sensitive to the value of \(f{\text{O}}_2\) chosen.     

Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

Water diffusion in Mount Changbai peralkaline rhyolitic melt

Diffusion couple experiments with wet half (up to 4.6 wt%) and dry half were carried out at 789–1,516 K and 0.47–1.42 GPa to investigate water diffusion in a peralkaline rhyolitic melt with major oxide concentrations matching Mount Changbai rhyolite. Combining data from this work and a related study, total water diffusivity in peralkaline rhyolitic melt can be expressed as:

$ D_{{{\text{H}}_{ 2} {\text{O}}_{\text{t}} }} = D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} \left( {1 - \frac{0.5 - X}{{\sqrt {[4\exp (3110/T - 1.876) - 1](X - X^{2} ) + 0.25} }}} \right), $

$ {\text{with}}\;D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} = \exp \left[ { - 1 2. 7 8 9- \frac{13939}{T} - 1229.6\frac{P}{T} + ( - 27.867 + \frac{60559}{T})X} \right], $

where D is in m² s^?1, T is the temperature in K, P is the pressure in GPa, and X is the mole fraction of water and calculated as X = (C/18.015)/(C/18.015 + (100 ? C)/33.14), where C is water content in wt%. We recommend this equation in modeling bubble growth and volcanic eruption dynamics in peralkaline rhyolitic eruptions, such as the ~1,000-ad eruption of Mount Changbai in North East China. Water diffusivities in peralkaline and metaluminous rhyolitic melts are comparable within a factor of 2, in contrast with the 1.0–2.6 orders of magnitude difference in viscosities. The decoupling of diffusivity of neutral molecular species from melt viscosity, i.e., the deviation from the inversely proportional relationship predicted by the Stokes–Einstein equation, might be attributed to the small size of H₂O molecules. With distinct viscosities but similar diffusivity, bubble growth controlled by diffusion in peralkaline and metaluminous rhyolitic melts follows similar parabolic curves. However, at low confining pressure or low water content, viscosity plays a larger role and bubble growth rate in peralkaline rhyolitic melt is much faster than that in metaluminous rhyolite.

     

The effect of liquid composition on the partitioning of Ni between olivine and silicate melt

We report the results of experiments designed to separate the effects of temperature and pressure from liquid composition on the partitioning of Ni between olivine and liquid, \(D_{\text{Ni}}^{\text{ol/liq}}\). Experiments were performed from 1300 to 1600 °C and 1 atm to 3.0 GPa, using mid-ocean ridge basalt (MORB) glass surrounded by powdered olivine in graphite–Pt double capsules at high pressure and powdered MORB in crucibles fabricated from single crystals of San Carlos olivine at one atmosphere. In these experiments, pressure and temperature were varied in such a way that we produced a series of liquids, each with an approximately constant composition (~12, ~15, and ~21 wt% MgO). Previously, we used a similar approach to show that \(D_{\text{Ni}}^{\text{ol/liq}}\) for a liquid with ~18 wt% MgO is a strong function of temperature. Combining the new data presented here with our previous results allows us to separate the effects of temperature from composition. We fit our data based on a Ni–Mg exchange reaction, which yields \(\ln \left( {D_{\text{Ni}}^{\text{molar}} } \right) = \frac{{ -\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{RT} + \frac{{\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{R} - \ln \left( {\frac{{X_{\text{MgO}}^{\text{liq}} }}{{X_{{{\text{MgSi}}_{ 0. 5} {\text{O}}_{ 2} }}^{\text{ol}} }}} \right).\) Each subset of constant composition experiments displays roughly the same temperature dependence of \(D_{\text{Ni}}^{\text{ol/liq}}\) (i.e.,\(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\)) as previously reported for liquids with ~18 wt% MgO. Fitting new data presented here (15 experiments) in conjunction with our 13 previously published experiments (those with ~18 wt% MgO in the silicate liquid) to the above expression gives \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 3641 ± 396 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 1.597 ± 0.229. Adding data from the literature yields \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 4505 ± 196 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 2.075 ± 0.120, a set of coefficients that leads to a predictive equation for \(D_{\text{Ni}}^{\text{ol/liq}}\) applicable to a wide range of melt compositions. We use the results of our work to model the melting of peridotite beneath lithosphere of varying thickness and show that: (1) a positive correlation between NiO in magnesian olivine phenocrysts and lithospheric thickness is expected given a temperature-dependent \(D_{\text{Ni}}^{\text{ol/liq}} ,\) and (2) the magnitude of the slope for natural samples is consistent with our experimentally determined temperature dependence. Alternative processes to generate the positive correlation between NiO in magnesian olivines and lithospheric thickness, such as the melting of olivine-free pyroxenite, are possible, but they are not required to explain the observed correlation of NiO concentration in initially crystallizing olivine with lithospheric thickness.     

Some geobarometers involving cordierite in the FeO-Al2O3-SiO2(±H2O) system: refinements,thermodynamic calibration,and applicability in granulite facies rocks

Five geobarometers involving cordierite have been formulated for quantitative pressure sensing in high grade metapelites. The relevant reactions in the FeO-Al₂O₃-SiO₂ (±H₂O) system are based on the assemblages (A) cordierite-garnet-sillimanite-quartz, (B) cordierite-spinel-quartz, (C) cordierite-garnet-spinel-sillimanite, (D) cordierite-garnet-orthopyroxene-quartz and (E) cordierite-orthopyroxene-sillimanite-quartz. Application of the barometric formulations to a large number of granulite grade rocks indicates that the cordierite-garnet-sillimanite-quartz equilibrium is widely applicable and registers pressures which are in good agreement with the “consensus” pressure estimates. The dispersion in the computed P values, expressed as one standard deviation, is within ±1.2 kbar. The geobarometers (B) and (C) also yield pressures which are reasonable and compare well with those computed from equilibrium (A). The estimated pressures from (D) and (E), both involving orthopyroxene, are at variance with these estimates. It has been argued that the discrepancy in pressures obtained from these geobarometers stems from an inadequate knowledge of activity-composition relations and/or errors in input thermodynamic data of aluminous orthopyroxene. The convergence of pressure values estimated from the barometric formulations, especially (A), (B) and (C), implies that the present formulations are more dependable than the existing formulations and are also capable of setting limits on P values in response to varying $$\begin{gathered} {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ {\text{ = 1/3Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2/3Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 5/6SiO}}_{\text{2}} {\text{. (A)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ = FeAl}}_{\text{2}} {\text{O}}_{\text{4}} {\text{ + 5/2SiO}}_{\text{2}} {\text{. (B)}} \hfill \\ {\text{Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + FeAl}}_{\text{2}} {\text{O}}_{\text{4}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{. (C)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 3/2SiO}}_{\text{2}} .{\text{ (D)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}{}_{\text{4}}{\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ = 1/2{\text{Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} {\text{ + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 1/2SiO}}_{\text{2}} .{\text{ (E)}} \hfill \\ \end{gathered}$$ . The present communication addresses the calibration, applicability and reliability of these barometers with reference to granulite facies metapelites.     

Phase relationships in granulitic metapelites from the Ivrea-Verbano zone (Northern Italy)

Paragneisses of the Ivrea-Verbano zone exhibit over a horizontal distance of 5 km mineralogical changes indicative of the transition from amphibolite to granulite facies metamorphism. The most obvious change is the progressive replacement of biotite by garnet via the reaction: a $${\text{Biotite + sillimanite + quartz }} \to {\text{ Garnet + K - feldspar + H}}_{\text{2}} {\text{O}}$$ which results in a systematic increase in the modal ratio g = (garnet)/(garnet + biotite) with increasing grade. The systematic variations in garnet and biotite contents of metapelites are also reflected by the compositions of these phases, both of which become more magnesian with increasing metamorphic grade. The pressure of metamorphism has been estimated from the Ca₃Al₂Si₃O₁₂ contents of garnets coexisting with plagioclase, sillimanite and quartz. These phases are related by the equilibrium: b $$\begin{gathered} 3 CaAl_2 {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} \rightleftharpoons Ca_3 Al_2 {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} + 2 Al_2 {\text{SiO}}_{\text{5}} + {\text{SiO}}_{\text{2}} \hfill \\ plagioclase garnet sillimanite quartz \hfill \\ \end{gathered} $$ which has been applied to these rocks using the available data on the mixing properties of plagioclase and garnet solid solutions. Temperature and f _{H 2O} estimates have been made in a similar way using thermodynamic data on the biotite-garnet reaction (a) and the approximate solidus temperatures of paragneisses. Amphibolite to granulite facies metamorphism in the Ivrea-Verbano zone took place in the P-T ranges 9–11 kb and 700–820 °C. The differences in temperature and pressure of metamorphism between g= 0 and g = 1 (5 kms horizontal distance) were less than 50° C and approximately 1 kb. Retrogression and re-equilibration of garnets and biotites in the metapelites extended to temperatures more than 50° C below and pressures more than 1.5 kb below the peak of metamorphism, the degree of retrogression increasing with decreasing grade of the metamorphic “peak”. The pressure and temperature of the peak of metamorphism are not inconsistent with the hypothesis that the Ivrea-Verbano zone is a slice of upthrusted lower crust from the crust-mantle transition region, although it appears that the thermal gradient was too low for the zone to represent a near-vertical section through the crust. The most reasonable explanation of the granulite facies metamorphism is that it arose through intrusion of mafic rocks into a region already undergoing recrystallisation under amphibolite facies conditions.     

The origin of garnet in the anorthosite-charnockite suite of the Adirondacks

Detailed analysis of textural and chemical criteria in rocks of the anorthosite-charnockite suite of the Adirondack Highlands suggests that development of garnet in silica-saturated rocks of the suite occurs according to the reaction: $$\begin{gathered} {\text{Anorthite}} {\text{Orthopyroxene}} {\text{Quartz}} \hfill \\ {\text{2CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} + (6 - \alpha )({\text{Fe,Mg}}){\text{SiO}}_{\text{3}} + \alpha {\text{Fe - Oxide + (}}\alpha {\text{ - 2)SiO}}_{\text{2}} \hfill \\ {\text{Garnet}} {\text{Clinopyroxene}} \hfill \\ = {\text{Ca(Fe,Mg)}}_{\text{5}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{6}} {\text{O}}_{{\text{24}}} + {\text{Ca(Fe,Mg)Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ \end{gathered} $$ , where α is a function of the distribution of Fe and Mg between the several coexisting ferromagnesian phases. Depending upon the relative amounts of Fe and Mg present, quartz may be either a reactant or a product. Using an aluminum-fixed reference frame, this reaction can be restated in terms of a set of balanced partial reactions describing the processes occurring in spatially separated domains within the rock. The fact that garnet invariably replaces plagioclase as opposed to the other reactant phases indicates that the aluminum-fixed model is valid as a first approximation. This reaction is univariant and produces unzoned garnet. It differs from a similar equation proposed by de Waard (1965) for the origin of garnet in Adirondack metabasic rocks, i.e. 6 Orthopyroxene+2 Anorthite = Clinopyroxene+Garnet+2 Quartz, the principle difference being that iron oxides (ilmenite and/or magnetite) are essential reactant phases in the present reactions. The product assemblage (garnet+clinopyroxene+plagioclase ± orthopyroxene ± quartz) is characteristic of the clinopyroxene-almandine subfacies of the granulite facies.     

Raman spectra of gypsum under pressure

Raman sprectra of a gypsum crystal were made at pressures between 0.001 and 7 kbar using He gas as the pressure medium. \(\frac{{{\text{d}}v}}{{dP}}\) values for bands in the range 3,600–100 cm^?1 were obtained. Comparison of results with \(\frac{{{\text{d}}v}}{{{\text{d}}T}}\) from the literature for temperatures of 77 and 300° K. shows that the internal modes of the SO₄ units are more sensitive to pressure than to temperature. The effect is small. Coupled H₂O-SO₄ translational modes are greatly affected by both pressure and temperature while coupled Ca-SO₄ mode are less so. It was found that stretching vibrations of water molecules were affected differently under pressure. The band at 3,500 cm^?1 is more greatly displaced by pressure \(\left( {\frac{{{\text{d}}v}}{{{\text{d}}P}} = {\text{2}}{\text{.11cm}}^{{\text{ - 1}}} /{\text{kbar}}} \right)\) than the band at 3,400 cm^?1 \(\left( {\frac{{{\text{d}}v}}{{{\text{d}}P}} \simeq {\text{2}}{\text{.11cm}}^{{\text{ - 1}}} /{\text{kbar}}} \right)\) . Assuming two different hydrogen bond intensities for the water molecules, one can attribute this difference in behavior of stretching modes to and increase in hydrogen bonding of one of the hydrogens which is exterior to the double H₂O planes in the gypsum structure. The great variety of pressure derivatives for the different types of vibrational modes observed indicates that each molecular unit readjusts internally to pressure induced volume changes and the some of the chemical bonds between the units are significantly affected.     

The thermal expansion and crystal structure of mirabilite (Na<Subscript>2</Subscript>SO<Subscript>4</Subscript>·10D<Subscript>2</Subscript>O) from 4.2 to 300 K,determined by time-of-flight neutron powder diffraction

We have collected high resolution neutron powder diffraction patterns from Na₂SO₄·10D₂O over the temperature range 4.2–300 K following rapid quenching in liquid nitrogen, and over a series of slow warming and cooling cycles. The crystal is monoclinic, space-group P2₁/c (Z = 4) with a = 11.44214(4) Å, b = 10.34276(4) Å, c = 12.75486(6) Å, β = 107.847(1)°, and V = 1436.794(8) Å³ at 4.2 K (slowly cooled), and a = 11.51472(6) Å, b = 10.36495(6) Å, c = 12.84651(7) Å, β = 107.7543(1)°, V = 1460.20(1) Å³ at 300 K. Structures were refined to R _P (Rietveld powder residual, \( R_{P} = {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } \mathord{\left/ {\vphantom {{\sum {\left| {I_{\text{obs}} - I_{\text{calc}} } \right|} } {\sum {I_{\text{obs}} } }}} \right. \kern-\nulldelimiterspace} {\sum {I_{\text{obs}} } }} \)) better than 2.5% at 4.2 K (quenched and slow cooled), 150 and 300 K. The sulfate disorder observed previously by Levy and Lisensky (Acta Cryst B34:3502–3510, 1978) was not present in our specimen, but we did observe changes with temperature in deuteron occupancies of the orientationally disordered water molecules coordinated to Na. The temperature dependence of the unit-cell volume from 4.2 to 300 K is well represented by a simple polynomial of the form V = ? 4.143(1) × 10^?7 T ³ + 0.00047(2) T² ? 0.027(2) T + 1437.0(1) Å³ (R ² = 99.98%). The coefficient of volume thermal expansion, α _V , is positive above 40 K, and displays a similar magnitude and temperature dependence to α _V in deuterated epsomite and meridianiite. The relationship between the magnitude and orientation of the principal axes of the thermal expansion tensor and the main structural elements are discussed; freezing in of deuteron disorder in the quenched specimen affects the thermal expansion, manifested most obviously as a change in the behaviour of the unit-cell parameter β.     

Thermodynamic properties of almandine-grossular garnet solid solutions

The mixing properties of Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnet solid solutions have been studied in the temperature range 850–1100° C. The experimental method involves measuring the composition of garnet in equilibrium with an assemblage in which the activity of the Ca₃Al₂Si₃O₁₂ component is fixed. Experiments on the assemblage garnet solid solution, anorthite, Al₂SiO₅ polymorph and quartz at known pressure and temperature fix the activity of the Ca₃Al₂Si₃O₁₂ component through the equilibrium: 1 $$\begin{gathered} {\text{3CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} \rightleftarrows {\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} \hfill \\ {\text{Anorthite garnet}} \hfill \\ {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + SiO}}_{\text{2}} \hfill \\ {\text{ sillimanite/kyanite quartz}}{\text{.}} \hfill \\ \end{gathered}$$ This equilibrium, with either sillimanite or kyanite as the aluminosilicate mineral, was used to control \({\text{a}}_{{\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} }^{{\text{gt}}} \) . The compositions of the garnet solutions produced were determined by measurement of their unit cell edges. At 1 bar Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnets exhibit negative deviations from ideality at the Fe-rich end of the series and positive deviations at the calcium end. With increasing pressure the activity coefficients for the Ca₃Al₂Si₃O₁₂ component increase because the partial molar volume of this component is greater than the molar volume of pure grossular. Previous studies indicate that the activity coefficients for the Ca₃Al₂Si₃O₁₂ component also increase with increasing (Mg/Mg+Fe) ratio of the garnet. The region of negative deviation from ideality implies a tendency towards formation of a stable Fe-Ca garnet component. Evidence in support of this conclusion has been found in a natural Fe-rich garnet which was found to contain two different garnet phases of distinctly different compositions.     

The high-pressure stability of chlorite and other hydrates in subduction mélanges: experiments in the system Cr<Subscript>2</Subscript>O<Subscript>3</Subscript>–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>–H<Subscript>2</Subscript>O

The solubility of chromium in chlorite as a function of pressure, temperature, and bulk composition was investigated in the system Cr₂O₃–MgO–Al₂O₃–SiO₂–H₂O, and its effect on phase relations evaluated. Three different compositions with X _Cr = Cr/(Cr + Al) = 0.075, 0.25, and 0.5 respectively, were investigated at 1.5–6.5 GPa, 650–900 °C. Cr-chlorite only occurs in the bulk composition with X _Cr = 0.075; otherwise, spinel and garnet are the major aluminous phases. In the experiments, Cr-chlorite coexists with enstatite up to 3.5 GPa, 800–850 °C, and with forsterite, pyrope, and spinel at higher pressure. At P > 5 GPa other hydrates occur: a Cr-bearing phase-HAPY (Mg_2.2Al_1.5Cr_0.1Si_1.1O₆(OH)₂) is stable in assemblage with pyrope, forsterite, and spinel; Mg-sursassite coexists at 6.0 GPa, 650 °C with forsterite and spinel and a new Cr-bearing phase, named 11.5 Å phase (Mg:Al:Si = 6.3:1.2:2.4) after the first diffraction peak observed in high-resolution X-ray diffraction pattern. Cr affects the stability of chlorite by shifting its breakdown reactions toward higher temperature, but Cr solubility at high pressure is reduced compared with the solubility observed in low-pressure occurrences in hydrothermal environments. Chromium partitions generally according to \(X_{\text{Cr}}^{\text{spinel}}\) ? \(X_{\text{Cr}}^{\text{opx}}\) > \(X_{\text{Cr}}^{\text{chlorite}}\) ≥ \(X_{\text{Cr}}^{\text{HAPY}}\) > \(X_{\text{Cr}}^{\text{garnet}}\). At 5 GPa, 750 °C (bulk with X _Cr = 0.075) equilibrium values are \(X_{\text{Cr}}^{\text{spinel}}\) = 0.27, \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.08, \(X_{\text{Cr}}^{\text{garnet}}\) = 0.05; at 5.4 GPa, 720 °C \(X_{\text{Cr}}^{\text{spinel}}\) = 0.33, \(X_{\text{Cr}}^{\text{HAPY}}\) = 0.06, and \(X_{\text{Cr}}^{\text{garnet}}\) = 0.04; and at 3.5 GPa, 850 °C \(X_{\text{Cr}}^{\text{opx}}\) = 0.12 and \(X_{\text{Cr}}^{\text{chlorite}}\) = 0.07. Results on Cr–Al partitioning between spinel and garnet suggest that at low temperature the spinel- to garnet-peridotite transition has a negative slope of 0.5 GPa/100 °C. The formation of phase-HAPY, in assemblage with garnet and spinel, at pressures above chlorite breakdown, provides a viable mechanism to promote H₂O transport in metasomatized ultramafic mélanges of subduction channels.     

H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].     

Majoritic garnet grains within shock-induced melt veins in amphibolites from the Ries impact crater suggest ultrahigh crystallization pressures between 18 and 9 GPa

Shock-induced melt veins in amphibolites from the Nördlinger Ries often have chemical compositions that are similar to bulk rock (i.e., basaltic), but there are other veins that are confined to chlorite-rich cracks that formed before the impact and these are poor in Ca and Na. Majoritic garnets within the shock veins show a broad chemical variation between three endmembers: (1) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}{\text{Al}}_{2} ({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (normal garnet, Grt), (2) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}[{\text{M}}^{2 + } ({\text{Si,Ti}})]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\)  (majorite, Maj), and (3) \({}^{\text{VIII}}({{\text {Na} {\text M}^{2+}}_2}) {}^{\text{VI}}[ ({\text{Si,Ti}}){\text {Al}}]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (Na-majorite₅₀Grt₅₀), whereby M²⁺ = Mg²⁺, Fe²⁺, Mn²⁺, Ca²⁺. In particular, we observed a broad variation in ^VI(Si,Ti) which ranges from 0.12 to 0.58 cations per formula unit (cpfu). All these majoritic garnets crystallized during shock pressure release at different ultrahigh pressures. Those with high ^VI(Si,Ti) (0.36–0.58 cpfu) formed at high pressures and temperatures from amphibole-rich melts, while majoritic garnets with lower ^VI(Si,Ti) of 0.12–0.27 cpfu formed at lower pressures and temperatures from chlorite-rich melts. Furthermore, majoritic garnets with intermediate values of ^VI(Si,Ti) (0.24–0.39) crystallized from melts with intermediate contents of Ca and Na. To the best of our knowledge the ‘MORB-type’ Ca–Na-rich majoritic garnets with maximum contents of 2.99 wt% Na₂O and calculated crystallisation pressures of 16–18 GPa are the most extreme representatives ever found in terrestrial shocked materials. At the Ries, the duration of the initial contact and compression stage at the central location of impact is estimated to only ~ 0.1 s. We used a ~ 200-µm-thick shock-induced vein in a moderately shocked amphibolite to model its pressure–temperature–time (P–T–t) path. The graphic model manifests a peak temperature of ~ 2600 °C for the vein, continuum pressure lasting for ~ 0.02 s, a quench duration of ~ 0.02 s and a shock pulse of ~ 0.038 s. The small difference between the continuum pressure and the pressure of majoritic garnet crystallization underlines the usefulness of applying crystallisation pressures of majoritic garnets from metabasites for calculation of dynamic shock pressures of host rocks. Majoritic garnets of chlorite provenance, however, are not suitable for the determination of continuum pressure since they crystallized relatively late during shock release. An extraordinary glass- and majorite-bearing amphibole fragment in a shock-vein of one amphibolite documents the whole unloading path.     

Experimental determination of the reaction paragonite=jadeite+kyanite+H2O,and internally consistent thermodynamic data for part of the system Na2O−Al2O3−SiO2−H2O,with applications to eclogites and blueschists

Hydrothermal reversal experiments have been performed on the upper pressure stability of paragonite in the temperature range 550–740 ° C. The reaction $$\begin{gathered} {\text{NaAl}}_{\text{3}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{1 0}}} ({\text{OH)}}_{\text{2}} \hfill \\ {\text{ paragonite}} \hfill \\ {\text{ = NaAlSi}}_{\text{2}} {\text{O}}_{\text{6}} + {\text{Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}} \hfill \\ {\text{ jadeite kyanite vapour}} \hfill \\ \end{gathered}$$ has been bracketed at 550 ° C, 600 ° C, 650 ° C, and 700 ° C, at pressures 24–26 kb, 24–25.5 kb, 24–25 kb, and 23–24.5 kb respectively. The reaction has a shallow negative slope (? 10 bar °C^?1) and is of geobarometric significance to the stability of the eclogite assemblage, omphacite+kyanite. The experimental brackets are thermodynamically consistent with the lower pressure reversals of Chatterjee (1970, 1972), and a set of thermodynamic data is presented which satisfies all the reversal brackets for six reactions in the system Na₂O-Al₂O₃-SiO₂-H₂O. The Modified Redlich Kwong equation for H₂O (Holloway, 1977) predicts fugacities which are too high to satisfy the reversals of this study. The P-T stabilities of important eclogite and blueschist assemblages involving omphacite, kyanite, lawsonite, Jadeite, albite, chloritoid, and almandine with paragonite have been calculated using thermodynamic data derived from this study.     

Weathering Processes in the Min Jiang: Major Elements, {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} }

We investigated the dissolved major elements, $ {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ composition of the Min Jiang, a headwater tributary of the Chang Jiang (Yangtze River). A forward calculation method was applied to quantify the relative contribution to the dissolved load from rain, evaporite, carbonate, and silicate reservoirs. Input from carbonate weathering dominated the major element composition (58–93%) and that from silicate weathering ranged from 2 to 18% in unperturbed Min Jiang watersheds. Most samples were supersaturated with respect to calcite, and the CO₂ partial pressures were similar to or up to ～5 times higher than atmospheric levels. The Sr concentrations in our samples were low (1.3–2.5 μM) with isotopic composition ranging from 0.7108 to 0.7127, suggesting some contribution from felsic silicates. The Si/(Na^* + K) ratios ranged from 0.5 to 2.5, which indicate low to moderate silicate weathering intensity. The $ \delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } \;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ for five select samples showed that the source of dissolved sulfate was combustion of locally consumed coal. The silicate weathering rates were 23–181 × 10³ mol/km²/year, and the CO₂ consumption rates were 31–246 × 10³ mol/km²/year, which are moderate on a global basis. Upon testing various climatic and geomorphic factors for correlation with the CO₂ consumption rate, the best correlation coefficients found were with water temperature (r ² = 0.284, p = 0.009), water discharge (r ² = 0.253, p = 0.014), and relief (r ² = 0.230, p = 0.019).