American Mineralogist,Volume76, pages 1931-1939, 1991 Tremolite: New enthalpy and entropy data from a phaseequilibrium study of the reaction tremolite : 2 diopside + 1.5 orthoenstatite+ 0-qurtz * HrO Mlm D. Wnr-cn,* Ausoll R. Plwr,Bv** Department of Geology and Geophysics,University of Edinburgh, West Mains Road, Edinburgh EH9 3Jw, United Kingdom Ansrucr New enthalpy and entropy data have been obtained for a natural end-member tremolite from phaseequilibrium experimentson the high-temperaturedehydration reaction tremolite : 2 diopside + 1.5 orthoenstatite * B-quartz + HrO. Six reversed experimental bracketswereobtainedat0.25kbar (996-1033K), 0.6 kbar (1038-1058K), I kbar (10751 1 0 4K ) , 1 . 5 k b a r ( l l 0 l - 1 1 2 5 K ) , 2 k b a r ( l l l 5 - 1 1 2 8 K ) , a n d 8 k b a r ( 1 1 7 9 - 1 2 0 9 K ) , from which the enthalpy of formation and the third-law entropy of tremolite have been calculated using the internally consistent thermodynamic datasetThermocalc of Holland and Powell (1990).The valuesobtainedare Af{(Tr) : -12299.48 + 14.84kJlmol and S{Tr): 0.5525 + 0.0100 kJ/mol.K, which are in excellentagreementwith the values given in Thermocalcand agreewell with the valuesof Berman et al. (1985).We use the new thermodynamic data to calculate revised limits of tremolite stability in the system CaO-MgO-SiOr-H,O (CMSH). We also discussthe apparent major discrepancybetween values for AIIi(Tr) obtained from solution calorimetry (Kiseleva and Ogorodova, 1984) and the phase equilibrium values. The difference between enthalpies for end-member tremolite determined from phase equilibrium and calorimetric studies (36 kJ) is comparable to that observedfor talc (Kiseleva and Ogorodova, 1984),perhapssuggestinganomalous dissolution behavior of the talclike elementsof tremolite. ria have been made in the systemsCMSH and CMSHTremolitic amphiboles are important constituents of CO, under geologically relevant conditions. These are hydrated basic and ultrabasic rocks and siliceousmarbles summarized in Table l, from which it is evident that in the greenschistand albite-amphibolite facies. These most of the amphiboles studied do not approach ideal amphiboles are the first to grow in these rocks, and as tremolite composition closely.The resultsof thesestudies such they are substratesfor a range of continuous reac- have been used as sourcesfor deriving thermodynamic tions (most importantly coupledexchangereactions);these data for tremolite (Table l). Not surprisingly, given the reactions can lead to zoned amphiboles of use in eluci- chemical diversity of the amphiboles used, there is poor dating the pressure-temperature-timeevolution of meta- agreement between thermodynamic datasets. Furthermorphic rocks (cf. Holland and Richardson, 1979). The more, the microstructural statesofthese amphiboleswere procurement of well constrainedthermodynamic data for not determined.Mareschand Czank (1988) have shown end-member tremolite CarMgrSirOrr(OH),is thereforean that the intensity shifts in X-ray diffraction patterns used important objective. However, despite the dozen or so to infer reaction direction in reversal experiments can experimental studiesof tremolite stability that have been also arise from the preferential elimination of chain-mulmade over the last 30 years, good agreementon experi- tiplicity faults (Czank and Liebau, 1980) in polysomatmentally derived thermodynamic data is lacking. Some ically disordered reactant amphibole, leading to the spurious location of amphibole equilibria in pressuresourcesofdiscrepancy are discussedbelow. space(seealso Graham et al., 1989, p. 70). temperature Boyd (1959) made the first experimental study of tremThe on synthetic tremolite deservefurther experiments olite stability in the syntheticsystemCaO-MgO-SiOr-HrO (CMSH) at H,O pressuresup to 2 kbar. He found that comment. End-member tremolite is notoriously difficult to synthesize:assemblagessuch as diopside + talc + amamphibole was stable at -800 t/0.6 kbar and -870 phibole + quartz, amphibole + diopside * enstatite,and "C/2kbar. A number of later studies of tremolite equilibamphibole + diopside + qvartz are common synthesis products. Complete yields of amphibole from tremolite * Present address:Department of Earth Sciences,University syntheseshave never been achieved, and yet end-memof Cambridge, Downing Street, Cambridge CB2 3EQ, United tremolite, although very rare, does occur in natureber Kingdom. ** Present address:Department of Chemistry, Arizona State the one usedin this study is one example,that of Kiseleva University, Tempe, Arizona 85287,U.S.A. and Ogorodova (1984) is another. Jenkins (1987) adINrnooucrroN 0003-o04x/91l1 I I 2-l 93l$02.00 l93l 1932 WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE Trele 1, Tremolitestarting materialsused in experimentalstudiesof tremolitestabilityand thermochemistry Authors Reaction/technioue Weeks(1956) Hf calorimetry Boyd (19s9) Tr : Di + En + Qz + HrO Robie& Stout (1963) Low-T calorimetry(Cp) Metz (1967,1976) Tr + Dol : Fo + Cc + CO, + H.O Phl + Cc + Oz : Tr + Ksp + CO, + H,O Phl + Cc + Qz : Tr + Ksp + Hoschek(1973) Hewitt (1975) co, + Hro slaughtheret al. (1975) and Eggert & Kerrick ( 19 8 1) Tr+Cc:Dol+Di+CO,+ H.O Dol+Qz+H,O:Tr+Cc+ CO, Tc+Cc+Qz:Tr+CO,+ H.o Dol+Tc+Oz:Tr+CO, Tr+Cc+Qz:Di+CO,+ Jenkins(1983) HrO Tr+Fo:Di+En+HrO Yin & Greenwood (1983) Tr:Di+En+Qz+HrO Skippen& McKinstry (1985) Kiseleva& Ogorodova (1984) Welch(1987) Tr+Fo:Di+En+HrO high-f lead borate calorimetry Tr:Di +En+Qz+HrO Tremolite composition Data set Ga/Mg (Ca,aMno@[email protected] Aloos)Orrlo(OH)rc <0.40 Synthetic tremolite (oxide mix, glass) (Nao11Gorxca1$Mgo,oXMguoo-0.37 Feo01XSi7 * seAlo o8)O*(OH), Nao,1(Ca,7sMgo27xMg4sAlo,.)-0.33 (Si7e3Aloo7)Oz(OH)1$ <0.40 Synthetic tremolite (gel) see Table3 FallsVillage,New York HP, BBG Naturaltremolite,Campo Lungo Naoos(Cale4Mgoo5NMgn."Feoor)0.385 HP Gouverneur,New York BBG, HP St. Gotthard, Switzerland <0.40 B B G ,H P =90% synthetic amphibole yield <0.40 BBG (si, &Alo,o)or2(oH),,4Fo06 (Na,K)o s5(Ca1?3Mgoro[Mgr zoFeoo5Alo,sXSiT osAlo,r)Orr(oH)osFo6 Synthetic tremolite (oxide mix) Synthetictremolite(gel,glass, oxide mix) Synthetic tremolile (gel, oxide mix) Ca, 16(Mg4 e4Feo o3XSi, erAloo,)o,,(oH)roFo03 Caro,Mgo..SirorOrr(OH),"u Naooloa, oo(MglesFeo oJsisooorr(oH)ro Skarn termolite,Balmat, New York H, HP HP, BBG 0.346 <0.40 Ca,16(Mg4e4Feo@) (siTerAlool)or2(oH)roFoo3 R 0.437 0.437 0.404 0.400 St. Gotthard, Switzerland Marble, Karelia Locality unknown. HRTEM (structurally ordered) lR, Raman Note.' From Graham et al. (1989). Abbreviations: Tr : tremolite, Di : diopside, En : enstatite, Qz : quartz, Dol : dolomite, Cc : calcite, Fo : forsterite,Phl : phlogopite,Ksp : potassiumfeldspar,Tc: talc, R : Robieet al. (1978),H : Helgesonet al. (1978),HP: Hollandand Powell(1985), BBG : Bermanet al. (1985). dressedthe problem of tremolite synthesisand presented a cogent case for incongruent dissolution of tremolite (specificallythe SiO, component) in supercritical HrO as a major control on experimental products. With increasing initial HrO contents,the incongruent dissolulion process leads to the following sequenceof assemblagesin which the fluids are SiO, undersaturated: that the stable amphibole was limited to within a few mole percent of TrroMc,o;Tr,oogave Tr-Mc amphibole * clinopyroxene. He also demonstratedfrom dissolution and regrowth experiments that the yield of amphibole could be increased by adding enough quartz to give a SiOr-saturated fluid. Following from this, Jenkins attempted to synthesizetremolite at 6 kbar/850 'C (well inside the nominal stability field calculated from therTr * H,O : Tr"" * Di + En + fluid (l) modynamic data sets)and obtained Tr-Mc amphibole * Tr * HrO: Di + En + fluid (2) clinopyroxene + qvartz, a high-temperatureassemblage Tr * H,O: Di + Fo + En + fluid (3) frequently encountered in previous studies. This result is Tr + HrO: Di + Fo + fluid. (4) strongly suggestedthat end-member tremolite unstable under theseconditions and starts to break down at a lowReaction I (seeAppendix I for abbreviations) involves er temperature by continuous reaction: Tr : Tr., * Cpx a change in amphibole chemistry from tremolite to a + Qz + fluid. The end-member natural tremolite samtremolite-magnesiocummingtonite(Tr-Mc) amphibole ples are then inferred to be stable at temperaturessignif(Jenkins, 1987).Hence,changesin modal proportions and icantly below the univariant dehydration of the limiting amphibole chemistry can occur. The quartz solubility data Mg-saturated Tr-Mc amphiboles: Tr." : Cpx + Opx + of Andersonand Burnham (1965)imply that assemblages Qz + fluid. Amphibole composition is inferred to change 3 and 4 are restrictedto bulk compositions with very high continuously along the univariant dehydration boundary H,O/Tr ratios (i.e., >90 wto/oHrO) at P < l0 kbar and (cf. Graham et al., 1989,their Fig. 9). An important corf < 900'C. Jenkins(1987) attemptedto synthesizeTr- ollary of the incongruent dissolution phenomenonis that Mc amphiboles at 5 molo/ointervals from tremolite to tremolite phaseequilibria in SiOr-saturatedsystemsmay Tr,,Mcr. at 6 kbar and l3 kbar (750-852 "C) and found be thermodynamically incompatible with those deter- WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE l 933 mined in SiOr-undersaturated (e.g., forsterite-bearing) TraLE2. Averagemicroprobeanalysisof naturaltremolite(n: systems. 31) A further complication with tremolite synthesesis the 2o Atoms rapid grofih of diopside (+ talc) in the first hour or so (sample) (O :24) Wt"/" of an experiment(Welch, 1987).When gelsare used,this (0.52) 7.95 sio, 58.61 problem may be even worse becauseof the presenceof (0.01) Tio, 0.o2 0.00 (0.02) Al,o3 0.03 0.01 diopside crystallites produced during firing of the gel at (0.01) CrrO" 0.01 0.00 850-900 "C. Similar crystallization behavior has been (0.02) FeO' 0.43 0.05 observed in some other calcic amphiboles such as eden(0.01) MnO 0.03 0.00 (0.30) Mgo 25.O1 5.06 ite, magnesio-hornblende,and richterite. In some cases (0.1 13.65 8) 1.98 CaO (e.g.,richterite), diopside can be reacted out by crushing (0.01) Naro o.02 0.00 (0.01) KrO o.02 0.00 and re-reacting,suggestingthat the diopside in synthesis (0.00) F.. 0.05 0.00 products is being preservedby an amphibole armor and (0.22) 2.19 1.98 H.ot thereby isolated from the reacting gel or oxide mix. (1.30) 100.07 17.03 This study was undertaken to derive well-constrained (+0.95)+ enthalpy and entropy data for end-member tremolite from * FerO"was not detected in wet-chemicalanalysis(M.J. Saunders,anphase-equilibrium experimentson the reaction Tr : 2Di alyst). ". Wet-chemicalanalysis (M.J. Saunders, analyst). + l.5En + B-Qz + HrO using a natural end-member HrO content determined by H extraction method (see Holdaway et tremolite sample and, if possible,to clarify the nature of a l .f, 1 9 8 6 ) . the breakdown. f This error is the total of analyticaluncertainties(> 2o). ExpnnrlrnnrAr, METHoDS Starting materials Diopside and orthorhombic enstatite were synthesized hydrothermally from stoichiometric gels at 4 kbar/900 1J and 2kbar/900 "C, respectively.Pure natural a-quartz that had been acid washed was used. The specimen of tremolite came from the University of Edinburgh collection, but its locality is unknown. A chemical analysis(the average of 3l microprobe analysesfor the cations and also F and HrO analysesdetermined by other methods) is given in Table 2, from which it can be seen that the specimen has the ideal tremolite composition and is F free. High-resolution transmission electron microscopy revealedthat the amphibole is exceptionally free of polysomatic (chain-width) and polytypic (stacking)faults and mineral inclusions. As such, this tremolite sample is the most suitable specimen for thermodynamic analysis of any tremolite yet studied (seeTable l). Cleavage fragments of tremolite were removed from the hand specimen(a tremolite-calcite-quartzrock) using a razor blade; then they were crushed and left to stand overnight in 30 volo/oethanoic acid to remove any carbonate residuesadhering to crystals.Examination by optical microscopy and powder X-ray diffraction showed that calcite was absent,the sample being amphibole with very minor amounts of quarlz. Two complimentary starting mixes were prepared: 850/oTr:150/oDi-En-Qz and 150/oTr : 850/oDi-En-Qz,the diopside, enstatite,and quartz being mixed in tremolite proportions (4:3:2). All phases were cmshed and sievedto < l5 prm.Approximately 0.02 g of mix was loaded into a Pt capsule(already welded at one end) along with 4-10 wto/odistilled HrO, and the capsulewas welded shut. Eachexperiment was conducted with pairs of mixes. At the conditions of our cold-seal experiments(0.25-2.0kbail700-930 oC),rhe data of Anderson and Burnham (1965) indicate that the hydrous fluids have SiO, molalities of -0. I mol SiO,/kg HrO, which correspondsto <0.01 molo/oSiO, in the fluid. With regard to our experimentsat 8 kbar/900-940'C, we have interpolateda SiO. molality of -1.7 at 900 "C from the data ofAnderson and Burnham (1965,their Fig. 4), corresponding to -3 molo/oSiO, in the fluid. Therefore, we consider the fluids in all our experimentsto be essentially pure HrO, while acknowledgingthat high quartz solubility (and by implication incongruent dissolution) is likely to exert an important control on tremolite stability at higher pressuresand temperatures(> l0 kbar/>800 'C). Apparatus Experimentsat 700-930'C and 0.25-2.0 kbar were performed using vertical cold-sealvessels(Nimonic 105), and those at 8 kbar used an internally heated gas vessel (Holloway, l97l). Ar was used as the pressuremedium. Temperature and pressurewere recorded at least once a day. Temperaturesof the experimentswere measuredusing an external PtlPtl3Rh thermocouple and are accurateto +3 'C (Ford, 1972). Experimental temperatures in experiments were extremely stable: diurnal and long-term variations (2o) were +2 "C. The temperature profiles of the furnaces used in this study were measuredat I atm and 600, 800, and 900 "C using an internal axial Ptl Ptl3Rh thermocoupleThe isothermalhot-spot zone (35 cm long) of each furnace was located and the furnace positioned relative to the vessel so that the charge zone (2 cm long) coincided with the hot spot. Vessel-furnace pairs were kept the same throughout the study. The internally heated pressure vessel was calibrated for temperature at 8 kbar using NaCl melting (detection of latent heat) and is accurateto + 5 'C. In all experiments,sampleswereheatedisobaricallyto temperature.Sampleswere 1934 WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE .t Fig. l. Scanningelectron micrographs of experimental products. (A) Tremolite growth (TD29: 0.6 kbar, 765 "C) in which most tremolite grains have retained their original prismatic facesand sharp outlines. Fretted terminations with protrusions can be seen 'C). on the large tremolite grain (arrowed and inset). (B) Typical product in which dehydration has occurred (TD13: I kbar, 831 quartz grains seen. can be also The view is mostly of pyroxene bladeswith well developed sharp faces;some small equant quenchedin a jet of compressedair. All quencheswere isobaric. Pressuresof the experimentswere measuredby a Heise gaugeand are accurateto +30 bars (+ 10 guaranteedaccuracy, +20 hysteresis),although they could be read to within 5 bars. Cold-seal pressurevesselswere closed off from the gauge.As a check on the extent to which the gauge + line volume buffered pressure drops, dummy experimentswere performedat 0.5-2 kbar and 650-900 'C in which the pressurevesselwas either (1) open continuously to the gaugeor (2) shut offfrom the gauge.By using segments of hollow filler rod it was possible to achieve pressuredrops that were very close to those observedupon opening a vesselto the gauge.Pressuredrops could thus be standardizedso that they correspondedto the true vessel pressures,i.e., gauge + line volume : vesselvolume. The internally heatedgasvesselwas calibrated for pressure (Manganin gauge)using the freezing point of mercury at 298 K, and our quoted pressuresare accurateto +20 bars. Pressurewas recorded by a Manganin gauge open continuously to the vessel. Experimental product examination electron imaging facility; element standardswere wollastonite and periclase;a current of 20 nA and accelerating potential of 20 kV were used. A number of products (including all those within and defining the experimental brackets) were examined by SEM to see if reaction textures could be recognizedand correlated with the X-ray results. SEM could also show the presenceof impurities below the detection limit of XRD (<50/o)and indicate whether dissolution had occurred(Jenkins, I 987). An analytical facility was used that enabled phasesto be identified and distinguished readily. Some representative electron micrographs are shown in Figure l. Figure lA is from an experiment at 0.6 kbar/165 "C (TD29) in which tremolite grew (X-rays). This view is dominated by tremolite grains that have retained their original prismatic faces and sharp outlines. The relatively large grain (arrowed and inset) showsclear ledgesand irregular sharp protruding terminations, the latter giving a "fretted" appearance. This texture is the same as that interpreted by Jenkins (1987, his Fie. 3b) to indicate tremolite erowth. The fretted habit is typical of experiments in which tremolite grew. Other tremolite grains commonly show distinct en-echelon growth ledges.In all casesthe SEM observationsare consistentwith the X-ray results.Evidently, the 250loXRD criterion is good enough to give unequivocal microtextures. We did not observetremolite with smooth rounded edges,such as are characteristicofpartially dissolved and recrystallized cleavagefragments (seeJenkins, 1987, his Fig. 3a), nor did we observe any impurities in grain mounts (RI oils) or SEM samples. Products were examined using the petrological microscope,powder X-ray diffraction (Ni-filtered CuKa, radiation), scanningelectron microscopy (SEM), and in some casesthe electron microprobe. Reaction direction was inferred from changesin the intensities of certain strong diagnostic X-ray reflections. A changein intensity of at least 25o/orelative to starting material diffractogramswas ANALYSTS TrrnnuolvNAMrc taken as evidence of reaction. Multiple repeatedsmears We now give details of our thermodynamic treatment were made to check for any orientational biases,but the of the phase equilibrium data. For internal consistency intensities of all reflections were reproducible to within 100/orelative. Microprobe analysis was performed using we have chosen a single thermodynamic dataset Thera Camebax microprobe equipped with a backscattered mocalc (Holland and Powell. 1990). This datasetwas WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE l 935 TAaLe3, Experimentalresultsfor the reactionfi : 2Di + 1.5En + B-Qz + H,O Experiment P(kbar)' TD1 TD2 TD3$ TD7 TD4 TD8 TDs$ TD6 TD28 TD29$ TD34 TD3O$ TD31 TD32 TD33 TD9 TD1O TD11$ TD15 TD12 TDl3$ TD14 TD22$ TD23 rD24 rD27 TD25$ TD26 TD16 TDl7 TD18$ TD19$ TD2O TD21 TD35$ TD36$ 0.241 0 245 0.244 o.237 o.240 0.222 o.225 0.231 0.616 0.609 0.598 0.604 0.611 0 627 0.596 0.998 1.O14 1.000 1.007 0.989 0.996 1. 0 1 0 1.487 1.472 1.476 1.419 1.477 1.473 1.974 1.986 1.990 1.982 1.985 1.956 8.020-8.060 7.960-7.990 rfc) Duration (h) 700 713 723 606 606 /JC 746 556 746 740 748 760 780 758 /oc 778 785 800 820 839 780 802 812 818 824 831 850 828 834 843 845 852 880 820 835 842 855 880 901 911 931 5b/r c5b 555 475 474 668 474 439 440 437 500 501 503 690 503 RNE 506 413 389 389 629 388 386 454 455 433 433 430 428 120 143 85Tr.. 85Prod'- Tr{ Tr Tr Tr+Pr Tr+Pr Tr+Pr Tr+Pr Pr Tr Tr Tr+Pr Tr+Pr Tr+Pr Pr Pr Tr Tr Tr+Pr Tr+Pr Tr+Pr Tr+Pr Tr+Pr Pt + ?Tl Pr ft+?Pr Tr+Pr Tr+Pr Pr + ?Tl Directionf H H H N N N D D H H N D D D D H H H N N Pl Pr Pr Tr ft+?Pr Tr+Pr Tr+Pr Tr+Pr Pr Pr Tr Tr+Pr Tr+Pr Tr+Pr Pr Pr Tr Tr Tr + ?Pr Pr + ?Tr Pr Pr Tr Pr Tr Tr Tr+Pr Tr+Pr Pr + ?Tr Pr Tr Tt + ?Pl Tr+Pr Tr+Pr Pr Pr Tr Tr Tr Tr+Pr Pr Pr tr Pr D H H N N D D H H H D D D H D - The given pressures are the mean for each experiment(see text). '. Starting materials(mineral mixes) t Interpretationof reaction direction: H : hydration (tremolitegrowth), D : dehydration,N : no reaction. + Experimentalproducts: Tr: tremolite, Pr : dehydrationproducts. $ Experimentdefining a bracket limit. chosenover that ofBerman et al. (1985) becauseit enables the derivation of uncertainties of enthalpies to be extracted. Six reversed experimental brackets were obtained for the reactionTr: 2Di + l.5En + P-Qz + HrO. Data for all experiments are given in Table 3. Experiments lasted 120-746 h. These six brackets were used to derive enthalpy and entropy data for tremolite as follows. At each point along the univariant reaction boundary, the following equilibrium relation holds: AG.:0:^H?- rasg* -T fT(.P * nRTlnfr" (6) + nRZln/",o + Rl"ln K.. RearrangingEquation 5 gives another having a more obvious linear form: + f|*,wdr* av,dp r". + RZln K". P'A(pY) + PAV\ + PL(aY)"(T-298)- - ^He: zA^te I)"or,dr r f'"or.,, u, f'"a,c,ar I A c- p / T d r +Jr| Jzge * [)"or,dr - r f'"or.,, o, o : Ane- z^,se (5) When expansivities and compressibilities are available and P >> I bar, then Equation 5 can be expandedthus: + RI ln K.. nRTlnf^,o (7) The right hand side of Equation 7 is combined into a single term here referred to as G'. Note that this is not the same as the G' of Jenkins(1983, 1984).By plotting WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE l 936 4.2 85 80 Iotercept = -AH"r = -114.23 ! Slope = ASor = 0 1640 t 11.24 kt 00100 8.O kJ K-' 7.4 70 36s ()60 2.O o 1100 1200 Fig.2. Plorof G' againstabsolutetemperature. Linesof maximum and minimum slopethroughall quadrilaterals definethe minimum-maximumuncertainties on 4119and AS9.Seetext. ! t.s o a o a o L o 1.o G' against absolute temperature, values for - AHI and A,$ are obtained as the intercept and slope, respectively. o.5 Water fugacitieswere taken from Thermocalc and are correct to within +0.1 kJ. This uncertainty amounts to < I 'C in the location of the tremolite dehydration reaction. The uncertainty on our 8 kbar data arising from the Thermocalc fit is <10. I kJ. With regard to the equilibrium-constant term, K., electron microprobe data for am900 1000 800 700 phiboles and pyroxenes from the highest- and lowest('C) Temperature temperature experimental brackets indicate that, within analytical error, these have end-member compositions. Fig. 3. Experimentalbracketsfor the reactionTr : 2Di + Hence, K. : I in all our calculations. l.5En + p-Qz+ HrO. Alsoshownis the minimum-maximum The Thermocalc data usedin our calculationsare given envelopecalculated from the minimum-maximumslopedataof in Table 4. Values for fnro and G' are given with the the G' against?Fplot. pressuresand temperaturesof the six experimentalbrackets (expandedby uncertaintiesin pressureof +0.02-0.03 mean and is kbar and temperature of t 5'C) in Table 5. A plot of G' The uncertainty on AIl(Tr) is 2o of the the enthalpies offoron from the uncertainties calculated against absolute temperature is shown in Figure 2. The p-quartz, and HrO orthoenstatite, mation of diopside, limits of the bracket data give rise to six quadrilaterals. A,F1!above. As Lines of maximum and minimum slope passingthrough given in Table 4 and the uncertainty on give uncertainties on third-law ennot Thermocalc does quadrilaterals all define an envelope within which the tropies, the uncertainty on S0(Tr) is taken as the uncerunivariant curve must lie. Maximal uncertainties (>2o) phaseequilibrium data : 114.23 tainty on AS! above. These new are derived from the two limiting slopes: MI? + 11.24kJ, ASP: 0.16400 + 0.0100kJ/K.The corre- for tremolite are discussedin the next section. lation coefficient relating A11! and AS! is 0.985. These Coup-lnrsoNs wITH PUBLISHED data were used to calculate the position of the tremolite DATA THERMODYNAMIC dehydration reaction: Figure 3 showsthe position defined A compilation of enthalpy and entropy data for tremby the above uncertainties, along with our experimental including our new values,is given in Table 6. There olite, (expanded With brackets by experimental uncertainties). is excellent agreement between our phase equilibrium respectto a mean line (not shown in Fig. 3), the reaction of Thermocalc, and good agteement with values and those is located to within + l0 "C by the experimental data. Berman et al. (1985). The value for AI1!(Tr) data of rhe Enthalpy and entropy data for tremolite were calculated of Robie et al. (1978), obtained by integratingthe C" data from Afl! and AS! as follows: ofRobie and Stout (1963),is clearly discrepant.The caAil!(Tr) : )taaKzoi, l.5En,B-Qz,H,o) - aF19 lorimetrically derived value for AH!(Tr) obtained by Kiselevaand Ogorodova(1984) deservessome discussion : -12299.48 + 14.84kJ/mol (8) as there is a major difference between our value for AHf(Tr) and theirs for a natural tremolite that is chemH,o) - a,sg s{Tr) : ) s.1zoi,l.5En,p-Qz, ically identical to ours in terms of cation proportions but : 0.5525+ 0.0010kJ/mol'K. (9) of unknown F content. To facilitate a direct, internally WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE t937 TABLe4. Thermodynamicdata used in calculatingtremoliteequilibria QO Phase LH? oAHl (x10.) Tremolite Diopside Orthoenstatite d-quanz B-quartz Forsterite Talc H"O -12299.48-3200.15 -3089.38 -910.80 -909.07 -2171.87 -5895.23 -241.81 7.42 1.90 2.08 0.74 0.74 1.62 3.70 0.03 552.50' 142.70 132.50 41.50 43.54 94.10 260.80 188.80 aV BV (x10+) (x10{) b Vo 27.2456.619 6.262 2.269 2.367 4.366 13.625 0.000 a 1.2144 0.3145 0.3562 0.0979 0.0979 0.2349 0.5343 0.0401 1x10 5) 2.6528 0.0041 -0.2990 -0.3350 -0.3350 0.1069 3.7416 0.8656 - 12363.0 -2745.9 596.9 -636.2 -636.2 542.9 -8805.2 487.5 -7.3885 -2.0201 -3.1853 -0.7740 -0.7740 - 1.9064 -2.1532 -0.2512 84.5 22.0 18.0 8.0 0.0 16.0 39.0 0.0 36.0 5.5 46 5.9 2.6 3.2 23.O 0.0 /Vote.'Unitsare in kJ, K, and kbar. Numbers in parenthesesare those bv which the values in the column below must be multiplied . This study. consistent comparison between the two studies,we have : recalculatedtheir AIIP(Tr) using their AI1"",,,,"".nrr"(Tr) 388.65 kJ/mol and a calorimetric cycle basedupon the reactiontremolite:2diopside + 2periclase+ 4 B-quartz + brucite. Although various other reactions could have beenchosen,this one resultsin the smallesterror on the derived AIl*,o,.,(Tr). Solution calorimetric data for diopside (Weill et al., 1980), periclase (Davies and Navrotsky, l98l), B-quartz (Akaogi and Navrotsky, 1984), and brucite(Clemenset al., 1987)wereusedwith the heat capacitiesfrom Thermocalc to calculate A11g(Tr)as fotlows: AffiTr) -- ) MI,1ZO|, 2Per, 4B-Qz, Bru) * Morurion.gzo*(Tr) + re16 AG | df. (10) J 29a From this we obtain a value for Afl!(Tr) of -12336-63 + 13.94kJ/mol. When uncertaintiesare considered,there is a minimum differenceof 8.37 kJ. Interestingly, an analTraue5. Valuesof G and f""ofor the experimental brackets P (kbao rfc) 0.202 0.286 0.194 o.276 0.569 0.649 0.560 0.648 0.949 1.051 0.939 1.053 1.426 1.548 1.423 1.531 1.940 2.040 1.9',t7 2.O47 8.000 8.080 7.940 8.010 718 718 765 765 760 760 790 790 807 807 836 836 823 823 857 857 837 837 860 860 906 906 936 936 fH.o(kbar) 0.151 0.261 0.133 o.245 0.532 0.594 0.534 0.596 0.838 0.912 o.847 0.935 1.204 1.296 1.236 1.321 1.623 1.705 1.639 1.749 11.067 11.270 11.'102 11.284 . Valuelying on the lineot minimumslope in d plot. '- Value lying on the line maximum of slope in G plot. G (kJ) 45.34 49.66. 46.85 51.71 57.61 58.34 59.58.60.48 63.82 64.34 65.93 oo.c/ 67.10 67.49 69.80 70.18 69.67 69.89 71.53' 71.83 79.77 79.77* 82.65 82.6s ogous calculation for talc using a value Af1,o,o,,on,nrr*(Tc) : 207.57 + 3.49 kJ/mol taken from Kiselevaand Ogorodova (1984) gives AH!(Tc) : -5923.57 + 10.50 kJl mol comparedwith a value of - 5895.23 + 7.40 kJ/mol derived from phaseequilibria experiments(Thermocalc). The difference is comparable to that observed for tremolite. The crystal-chemicalsimilarities betweentremolite and talc (e.g.,Thompson, 198l) suggesta genuine,possibly anomalous, calorimetric effect associatedwith dissolution of the talc-like elements of the tremolite structure in lead-borate flux. However, we should also point out that as Kiseleva and Ogorodova did not analyze for F, their value for A.F1"o,,,,o,(Tr) may result from a significantly higher F content than ours. In fact, their value is slightly more endothermic than that obtained for synthetic end-member fluortremolite by Graham and Navrotsky(1986). C.ql,cur,arroN oF THE Lrvrrrs oF TREMOLITESTABILITYIN CMSH Our values for the enthalpy and entropy of tremolite derived from phase equilibrium experiments have been used to calculate revised limits of tremolite stability in the CMSH system to 30 kbar. These are shown in Figure 4. In the absenceof any published uncertainties on the third-law entropies of diopside, orthoenstatite,forsterite, quartz, and HrO and the lack of a systematic rationale for combining uncertainties in enthalpy and entropy, we TaBLE of enthalpyandentropydatafor tremolite 6. Compilation Source Robieet al. (1978) Helgesonet al. (1978) Bermanet al. (1985) Holland and Powell (1990) Kiselevaand Ogorodova (1984)f This study AHP (kJ/mol) So (kJ/mol'K) -12355.08 (17.32)' 0.54890'," (0.001 3) -12319.70 0.54890*(0.0013) - 1231 1.30 0.5s115 -1230250 (14.121 0.55000 - 12347.s3(10.88)+ -12336.63(13.94)$ -123OO.47 (14.84) 0.55250 (0.0100) ' Uncertaintiesin brackets are 2r. '- From Robie and Stout (1963); calorimetricallydeterminedvalue. t High temperature solution calorimetry. Their sample was a natural end-membertremolite (possible F) from a Karelianmarble. + Value given by Kiselevaand Ogorodova. $ Valuecalculatedusing internallyconsistentdata from Thermocalcand the heats of solutionfor diopside,periclase,brucite,and B-quartzspecified in the text. WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE 19 3 8 for tremolite give the following thermodynamic quantities:A11! : 4.0 + 0.7 kJ; A.sg: -0.0063 + 0.0010kJ/ -0.382 kJ/kbar. Hence,there will be very small K AZ3 : free-energydifferentials at experimental temperaturesto drive this reaction: AG! is only about 6.5 kJ. The limits set upon the location ofthis reaction by the uncertainty in the enthalpy of reaction are shown in Figure 4 and amount to + 1.5 kbar. Combining this with our envelope for the SiOr-saturated tremolite dehydration reaction gives the location of the invariant point as 27-31 kbar/730800 "C. In passingwe note that our +2d envelopepasses through the bracketsdefined by some recent experiments on this reaction by Jenkins et al. (unpublished data) using synthetic tremolite (T. Holland, personal communication). 20 16 1, 111bar) l2 CoNcr-usroN 700 800 r (t) Fig. 4. Stabilityrelationsof tremolitein the CMSH system calculatedfrom the thermodynamicdata for the naturalstoichiometrictremoliteof this study.J: Jenkins(1983)for the reactionTr * Fo : 2Di + 2.5En+ HrO, usingsynthetictremolite. Seetext for discussion. have chosento use the method of calculating correlated uncertaintiesadvocatedby Powell and Holland (1985), in which only uncertaintiesin enthalpies are considered. The correlation matrices used were taken from Table 3a of Holland and Powell (1985). 2Di + 1.5En + Qz + H,O The minimum-maximum derivation of the median reaction curve from the G' vs. T plot implies that at 853 "C the uncertainty in its location is zero. Clearly, this is unrealistic, and so we have calculatedthe uncertainty on the reaction in the same way as we have for Tr : 2Di + Tc and Tr * Fo : zDi + 2.5En + HrO, using only the uncertainties on the enthalpies of formation given in Thermocalc and our uncertainty on AIIf(Tr) (o : +7.4 kI. This reactionis locatedto t6-12'C and reachesa thermal maximum of 925 "C at about 8 kbar. Tr: Tr * Fo: 2Di + 2.5En + H,O The location ofthis reaction calculated from our data is in good agreementwith the experimental determination of Jenkins (1983) for synthetic"tremolite" (Tr-Mc amphibole?).We obtain an uncertainty in the position of this reactionof + l0-13 oCover the pressurerange l-16 kbar. The position of the reaction is shown as an envelope in Figure 4. Tr: 2Di + Tc This is a "biopyribole" reaction(seeThompson 1981, his Table 3) in which amphibole disproportionates into plroxene and sheet silicate. Our phase equilibrium data We have presentednew enthalpy and entropy data for a natural end-member tremolite determined from phase equilibrium experiments.The data are in excellentagreement with those of Thermocalc and agreewell with those of Berman et al. (1985). The uncertainty on AHKTr) (2o : + 14.84 kJ) is comparable to that of Thermocalc (+14.12 kJ) and is necessarilyquite largebecausetremolite is a "non-anchor phase" (terminology of Holland and Powell, 1990) in the dataset and therefore incorporates uncertainties on substituent phases. Calculations made using the phase equilibrium values give plausible positions for some important tremolite equilibria in CMSH. The calculated position of the reaction Tr * Fo : 2Di + 2.5En + HrO agreeswell with the experimental data ofJenkins (1983). However, at P > l0 kbar and 7" > 800'C it is likely that incongruent dissolution of tremolite will occur, stabilising Tr-Mc amphiboles and thereby displacing the univariant dehydration reaction Tr.. : Di + En + HrO to temperaturesabove the end-member reaction. In this respect the calculated curve for endmember tremolite servesas a referencedatum. Acxlrowr,rncMENTS M.D.W. and A.R.P. acknowledgefinancial support from NERC research studentships. The high-pressure and electron microprobe facilities at University of Edinburgh are supported by NERC. We thank Peter Hill for advice on microprobe analysis, and Colin Graham and Tim Holland for helpful discussions. Frank Spear and David Jenkins are thanked for very thorough and helpful revrews. RnrnnnNcBscrrED transAkaogi, M., and Nawotsky, A. (1984)The quartz-coesite-stishovite formations and calculation ofphase diagrams- Physics ofthe Earth and PlanetaryInteriors, 36, 124-134. Anderson, G.M., and Burnham, C.W. (1965) The solubility of quartz in supercriticalwater. American Joumal of Science,263, 499-511. Berman, R.G., Brown, T.H., and Greenwood, H.J. (1985) An internally consistentthermodynamic databasefor minerals in the systemNa'OK,O-CaO-MgO-FeO-Fe'Or-AlrO3-SiOr-TiOr-HrO-COr.Technical Record TR377. Scientific Document Distribution Office Atomic Energy Authority of CanadaLtd., ResearchCompany, Chalk River, Ontario. l 939 WELCH AND PAWLEY: ENTHALPY AND ENTROPY OF TREMOLITE Boyd, F.R. (1959) Hydrothermal investigations of the amphiboles. In P.H. Abelson, Ed., Researchesin geochemistry,p. 377-396. Wiley, New York. Clemens,J.D., Circone, S., Nawotsky, A., McMillan, p.F., and Wall, y.J. ( I 987) Phlogopite: High-temperature solution calorirnetry, thermodynamic properties, Al-Si and stacking disorder, and phase equilibria. Geochimica et Cosmochimica Acta. 51. 2569-2578 Czank, M., and Liebau, F. (1980) Periodicity faults in chain silicates:A new type ofplanar lattice fault observedwith high-resolutionelectron microscopy. Physicsand Chemistry of Minerals, 6, 85-93 Davies, P.K., and Nawotsky, A. (1981) Thermodynamics of solid solution formation in NiO-MgO and NiO-ZnO. Joumal of Solid Srate Chemistry,38,264-276. Eggert, R.G., and Kerrick, D.M. (1981) Metamorphic equilibria in the siliceousdolomite system: 6 kbar experirnentaldata and geologicimplications. Geochimica et CosmochimicaActa. 45. 1039-1049. Ford, C.E. (1972) Temperature distribution in "Tuttle" hydrothermal pressurevessels.Progressin experimentalpetrology. Natural Environmental ResearchCouncil Publication SeriesD. 2. 87-88. Graham, C.M., and Navrotsky, A. (1986) Thermochemistry of rhe tremolite-edeniteamphiboles using fluorine analogues,and applications to amphibole-plagioclase-quartzequilibria. Contributions to Mineralogy and Petrology, 93, 18-32Graham, C.M., Maresch,W.V., Welch, M.D., and Pawley,A.R (1989) Experimental studies on amphiboles: A review with thermodynamic perspectives.EuropeanJournal ofMineralogy, l, 63-83. Helgeson, H.C., Delaney, J.M., Nesbitt, H.W., and Bird, D.K. (1978) Summary and critique of the thermodynamic proprertiesof minerals. American Journal of Science,278A. l-229. Hewitt, D.A. (1975) Stability of the assemblagephlogopite-calcite-quartz. American Mineralogist 60, 39 l-397. Holdaway, M.J., Dutrow, B.L., Borthwick, J., Shore, P., Harmon, R.S., and Hinton, R.W. (1986) H content of staurolite as determined by H extraction line and ion microprobe. American Mineralogist, 71, I 135I t4t. Holland, T.J.B., and Powell, R. (1985) An internally consistenrthermodynamic dataset with uncertainties and correlations: 2. Data and results. Journal of Metamorphic Geology, 3, 232-370. -(1990) An enlargedand updated internally consistentthermodynamic dataset with uncertainties and correlations: The system KrONarO-CaO-MgO-MnO-FeO-FerO3-Al,Or-TiOr-SiO,-C-H.-Or.Journal of Metamorphic Geology, 8, 89-124. Holland, T.J.B., and Richardson, S.W. (1979) Amphibole zonation in metabasitesas a guide to the evolution of metamorphic conditions. Contributions to Mineralogy and Perrology,70,143-148. Holloway, J.R. (1971) Internally-heatedgas vessels.In G.C. Ulmer, Ed., Research techniques for high pressuresand high temperatures. Springer-Verlag, New York. Hoschek,G. (1973) Die reaktion Phlogopir + Calcir + Quartz: Tremolit + Alkalifeldspat + H,O + COr. Contributions to Mineralogy and petrology,39, 231-237. Jenkins, D.M. (1983) Stability and composition relations of calcic amphiboles in ultramafic rocks. Contributions to Mineralogy and Petrology, 83, 375-384. (1984) Upper-presswestability of synthetic margariteplus quartz. Contributions to Mineralogy and Petrology,88, 332-339. -(1987) Synthesisand characterizationof tremolite in the system HrO-CaO-MgO-SiOr. American Mineralogist, 72, 707-7 | 5. Kiseleva, I.A., and Ogorodova, L.P. (1984) High-temperature solution calorimetry for determining the enthalpies of formation for hydroxylcontaining minerals such as talc and tremolite. Geochemistry International. 12.36-46. Maresch,W.V., and Czank, M. (1988) Crystal chemibtry,growth kinetics and phase relationships of structurally disordered [Mnr-, Mg] amphiboles Fortschritte der Mineralogie, 66, 69-121. Metz, P. (l 967) ExperimentelleBildung von Forsterit und Calcit aus Tremolit und Dolomit. Geochimica et Cosmochimica Acta, 31, l5l7- r532. -(1976) Experimental investigarion of the metamorphism of siliceousdolomites Ill-equilibrium data for the reaction: I Tremolite + ll Dolomite : 8 Forsterite + 13 Calcite + 9 CO, + 1l H"O for the total pressureof 3000 and 5000 bars. Contributions to Mineralogy and Petrology,58, 137-148. Powell, R., and Holland, T.J.B. (1985) An intemally-consistentrhermodynamic datasetwith uncertaintiesand correlations: l. Methods and a worked example.Journal of Metamorphic Geology, 3, 327-342. Robie, R.A., and Stout, J.W. (1963) Heat capacity from 12-305 K, and entropy of talc and tremolite. Journal of PhysicalCh emistry, 67, 22522256. Robie, R.A., Hemingway, B.S., and Fisher, J.R. (1978) Thermodynamic properties of minerals and related substancesat 298.15 K and one atmospherepressureand at higher temp€ratures.U.S. GeologicalSurvey Bulletin 1259,456. Skippen, G., and McKinstry, B.W- (1985) Syntheticand natural tremolite in equilibrium with forsterite, enstatite,diopside and fluid. Contributions to Mineralogy and Petrology, 89, 256-262. Slaughter, J., Kerrick, D.M., and Wall, V.J. (1975) Experimental and thermodynamic study of equilibria in the systemCaO-MgO-SiO,-H,OCO,. American Joumal of Science,275, 143-165. Thompson, J.B ( I 98 l) An introduction to the mineralogy and petrology of the bioppiboles. In Mineralogical Society of America Reviews in Mineralogy,9A, 141-188. Weeks,W.F. (1956) Heats of formation of metamorphic minerals in the systemCaO-MgO-SiOlHrO and their petrologicalsignificance.Journal of Geology, 64, 456-472. Weill, D.F., Hon, R., and Navrotsky, A. (1980) The igneous system CaMgSirOu.CaAlrSirOrNaAlSi,O,:Variations on a classic theme by Bowen. In R.B. Hargreaves,Ed., Physicsofmagmatic processes,p. 4992. Princeton University Press,Princeton, New Jersey. Welch, M.D. (1987) Experimental studies of amphiboles in the system Na,O-CaO-MgO-AlrOj-SiOr-HrO-Si4 and its subsystems,p. I 0 1-l 36. Unpublished Ph.D. thesis,University of Edinburgh,Edinburgh,Unired Kingdom. Yin, H.A., and Greenwood, H.J. (1983) Displacement of equilibria of OH-tremolite and F-tremolite solid solution. I Determination of the equilibrium P-T curve of OH-tremolite. Ens, 64, 347. MaNuscnrrr REcETVED Aucusr 23, 1990 MeNuscnrsr AccEmEDJuNe 26, l99l ApprNorx l. Glossary of thermodynamicterms used in this paper P Total pressure(kbar) T Temperature ('C, K) VP Volume of a solid phase(U/kbar) R Universal gas constant (0.008314 kJlK) 4113 Molar enthalpy of formation of a phase from the elements at 298.15 K and I bar (kJ/mol) So Molar thirdJaw entropy of a phase at 298.15 K and I bar (kJlmol.K) AHP Enthalpy of reaction at 298.15 K and I bar (kJ) A,SP Entropy of reaction at 298.15 K and I bar (kJlK) C, Molar heat capacityof a phaseat constant P (kJlmol'K), formulated as a polynomial in temperatureiCp: a + bT' + CT-2+ dTt/z aV Product of volume and isobaric expansivity for a solid phaseftJ/K.kbar) Product of volume and isothermal compressibility for a 0V solid phase(kJ/kbar) List of formulae and abbreviationsusedin this paper tremolite (Tr) CarMgrSirOrr(OH), Magnesiocummingtonite(Mc) Mg,Si,O,,(OH), Diopside (Di) CaMgSi,Ou Orthoenstatite (En) MgrSirOu sio, Quartz (a, B) (Qz) Forsterite (Fo) MgrSiOn Talc (Tc) MgrSioO,o(OH),
© Copyright 2024 Paperzz