榴辉岩中单斜辉石-石榴子石镁同位素地质温度计评述

榴辉岩中单斜辉石和石榴子石之间显著的镁同位素平衡分馏,使其成为一种具有潜力的高精度地质温度计。为此本文选取文献中已报道的来自西南天山洋壳冷俯冲造山带、大别-苏鲁陆壳碰撞造山带和南非卡普瓦尔克拉通三种构造环境中的64对单斜辉石—石榴子石矿物对镁同位素数据利用δ~(26)Mg_(Cpx)-δ~(26)Mg_(Grt)图解筛选出50对达到镁同位素平衡分馏的数据,采用Huang et al.(2013)通过第一性原理计算和Wang et al.(2012)、Li et al.(2016)分别通过经验估计得出的镁同位素温度计计算榴辉岩的峰期温度并与前人研究结果进行对比。通过分析计算结果发现对于造山带榴辉岩,Huang et al.(2013)的温度计计算结果与前人通过传统温度计和相平衡模拟研究结果较一致,较好地重现了榴辉岩的峰期温度,而Wang et al.(2012)和Li et al.(2016)的温度计计算结果明显低于前人获得的峰期温度;对于克拉通榴辉岩,三种镁同位素温度计的计算结果与前人通过传统温度计获得的研究结果大多数相差在50℃以上,很可能是早期退变质过程中镁同位素在高温条件下再平衡导致的这表明三种镁同位素温度计均不适用于克拉通榴辉岩。同时,基于这些榴辉岩样品数据,本文利用经验估计的方法进行校准,得到了新的单斜辉石-石榴子石镁同位素公式△~(26)Mg_(Cpx-Grt)=1.11×10~6/[T(K)]~2(R~2=0.92)。此外,本文也对单斜辉石-石榴子石镁同位素地质温度计的应用前景及应用时应当注意的问题进行了简单的探讨。     

榴辉岩富铝单斜辉石中出溶石榴子石现象研究

本文研究了江苏孟中榴辉岩富铝单斜辉石中出溶石榴子石现象，建立了出溶母体辉石成分，提出了出溶石榴子石的矿物反应。结果表明，富铝单斜辉石首先出溶富镁石榴子石，并形成绿辉石。后二者正是榴辉岩的典型矿物组合，因此，母体辉石可能正代表了榴辉岩的母体。本文中变质成因榴辉宕中出溶现象的研究尚属少见。     

榴辉岩富铝单斜辉石出溶石榴子石现象研究

本文研究了江苏孟中榴辉岩富铝单斜辉石中出溶石榴子石现象，建立了出溶母体辉石成分，提出了出溶石榴子石的矿物反应。     

孟中榴辉岩形成和演化和P-T轨迹

本文研究了江苏东海县孟中榴辉岩富铝单斜辉石中石榴子石出溶作用,利用温压计对榴辉岩中各阶段矿物组合进行了温压估算,由此建立了孟中榴辉岩形成和演化的P-T轨迹。结果表明,该榴辉岩经历了地幔条件下由于降温和可能的升压作用,富铝单斜辉石中出溶石榴子石并形成榴辉岩;榴辉岩构造侵位到高角闪岩相片麻岩基底中;退化变质这3个阶段的演化过程。      

北京密云变质杂岩矿物学和结晶的P-T条件

密云变质杂岩分布于北京密云县半城子以南到沙厂一带,由上壳岩和侵入岩组成。对麻粒岩、透辉斜长角闪岩中的石榴子石、斜方辉石,单斜辉石、角闪石和斜长石进行了电子探针分析;研究了镁铁矿物中的Fe、Mg、Ca和Na的分配系数,认为达到了化学平衡。根据麻粒岩中变质基性岩脉片麻理方向与麻粒岩的不同,本区可区分出两次变质作用。 本文应用4种二辉石地质温度计和5种石榴子石-单斜辉石地质温度计公式以及一些可靠的地质压力计公式,计算了变质作用的P-T条件。第一次变质作用为760—870℃和1.1—1.3GPa;第二次为740—850℃和0.8—0.84GPa。     

榴辉岩常用温压计在应用中应注意的问题

本文通过再现相平衡实验数据和检查热力学活度模型两种手段,对榴辉岩中几种常用的温压计进行了检验,发现榴辉岩中某些常用温压计存在以下问题:(1)在 Eliis and Green(1979)、Powell(1985)、Krogh(1988)和 Ravna(2000)四种石榴石-单斜辉石温度计中只有 Ravna(2000)的版本能较好的再现相平衡实验数据。(2)将石榴石-单斜辉石温度计应用于含 X_(Jd)>0.55绿辉石的榴辉岩中会出现很大的误差。(3)Green and Hellman(1982)的石榴石-多硅白云母温度计计算的高压含多硅白云母榴辉岩变质温度普遍偏高,但是计算超高压榴辉岩的结果能较好的与 Ravna(2000)的石榴石-单斜辉石温度计计算结果保持一致。(4)Waters and Martin(1993)的石榴石-单斜辉石-多硅白云母压力计、Ravna and Terry(2004)的石榴石-单斜辉石-多硅白云母-蓝晶石-柯石英/石英温压计的精度都受到了 Holland(1990)的单斜辉石活度模型的限制,它们不能适用于绿辉石 X_(Jd)>0.55的榴辉岩,而 Waters and Manin(1996)对 Waters and Martin(1993)的版本做了一个经验校正,弥补了单斜辉石活度模型的缺陷,因此可以适用于绿辉石 X_(Jd)>0.55的榴辉岩。(5)Nakamura and Banno(1997)的石榴石-绿辉石-蓝晶石-柯石英温压计因运用了不恰当的石榴石和铁钙辉石的活度模型,从而使得计算结果与岩相学观察结果不一致。因此,我们建议:对于绿辉石 X_(Jd)<0.55的多硅白云母榴辉岩,可以运用 Waters and Martin(1993)压力计和 Ravna(2000)的温度计联合求解温压;对于含高硬玉组分(X_(Jd)>0.55)绿辉石的超高压多硅白云母榴辉岩,可选用 Waters and Martin(1996)压力计和 Green andHellman(1982)的温度计联合求解温压;对于含有石榴石 绿辉石 蓝晶石 柯石英矿物组合的榴辉岩,在 X_(Jd)<0.55的情况下。可选用 Ravna and Terry(2004)的温压计求解温压。在应用这些温压计时,应注意各温压计的适用温压范围和成分范围,尤其是石榴石 X_(Jd)、Mg~#和绿辉石 X_(Jd)的范围。另外,由于矿物中 Fe~(3 )的含量对温度计算结果影响很大,所以还必须合理地校正所选矿物的 Fe~(3 )。     

青藏高原拉萨地块松多榴辉岩的岩相学特征和变质演化过程

松多榴辉岩出露于拉萨地块的石英片岩中,主要由较为基性的金红石榴辉岩和较为酸性的石英榴辉岩组成。榴辉岩相矿物组合为石榴子石 绿辉石 绿帘石±多硅白云母±石英±金红石。岩石发生了较强烈的退变质作用,退变质矿物有角闪石、绿帘石、石英、钠长石及绿泥石。石榴子石变斑晶具有生长环带结构,变斑晶和基质石榴子石主要落入C类榴辉岩区,少数石榴子石变斑晶边部和基质石榴子石落入B类榴辉岩区;单斜辉石主要为绿辉石,少数Ⅰ世代和Ⅲ世代为普通辉石;角闪石均为钙质角闪石。根据石榴子石-绿辉石-多硅白云母矿物温压计计算,获得的温压范围为630～777℃和2.58～2.70GPa,峰期变质条件接近于石英-柯石英转变线。榴辉岩的原岩经历了从高绿片岩相、角闪岩相、榴辉岩相、角闪岩相到高绿片岩相的变质过程,这反映了与古特提斯洋闭合有关的俯冲进变质作用和随后的折返退变质作用。     

大别山碧溪岭榴辉岩峰期变质条件地质温度计适应性探讨

榴辉岩石榴石-单斜辉石(Gt-Cpx)地质温度计常用于研究超高压变质作用p-t条件。业已存在多种Gt-Cpx温度计算公式和绿辉石Fe2+校正方法,本文通过碧溪岭地区浅色及深色2种榴辉岩p-t条件计算,对几种常见温度计算公式和绿辉石Fe2+校正方法进行比较,认为Droop(1987)电价平衡法和Ravna(2000)温度计算公式最适合榴辉岩的温度计算。计算结果表明,碧溪岭地区浅色榴辉岩峰期变质p-t条件为3.0～3.4GPa和660～731℃。     

南乌拉尔玛克苏托夫榴辉岩蓝片岩杂岩体的岩石学

玛克苏托夫 (Maksyutov)榴辉岩是乌拉尔—昆仑—祁连—秦岭—大别山元古代褶皱带之一部分。榴辉岩呈团块或布丁夹于片麻岩内 ,由于显生宙与花岗岩侵入有关的退变质叠加影响 ,时代变化自元古代至下古生代。榴辉岩矿物组分恒定但不同程度受绿帘蓝片岩相退变质叠加 ,榴辉岩中石榴石以富铁为主 (f=50 %～ 75% ) ;单斜辉石为绿辉石 Jd40 .3～ 52 .2。据 Krogh石榴子石—单斜辉石温度计 ,榴辉岩形成温度为 1 0 60～ 686℃ ,压力为 2 2 .3× 1 0 5～ 1 4.1 6× 1 0 5k Pa;PT轨迹具岩浆 (榴辉岩浆 )自上地幔向下地壳上侵的特点。     

北大别的高温超高压变质作用与多阶段折返

近期的变质岩石学、地球化学及同位素年代学研究表明,北大别整体经历了高温超高压变质作用和多阶段折返历史,因而表现为广泛发育的多期减压结构和极少保留早期的超高压变质记录。北大别榴辉岩以高温变质作用以及折返期间因麻粒岩相和角闪岩相退变质变质作用而形成的多期后成合晶为显著特征。石榴子石中伴有放射状胀裂纹的单晶和多晶石英包体指示早期柯石英的转变结果,这已被锆石中发现的柯石英残晶所证实。结合北大别北部榴辉岩和片麻岩中发现的金刚石等超高压证据以及三叠纪变质记录,由此证明北大别整体经历过深俯冲和印支期超高压变质作用。北大别榴辉岩的多阶段高温条件主要来自石榴子石-绿辉石矿物对温度计、单斜辉石中紫苏辉石+石英针状矿物出熔体以及金红石中较高的Zr含量和变质锆石中较高的Ti含量等得出的温度证据。此外,多期后成合晶以及石榴子石和单斜辉石等矿物中成分分带的存在,证明该区榴辉岩经历了一个多阶段、快速折返过程;而不同变质阶段的温度、压力和形成时代,却反映该区榴辉岩经历了长时间的高温变质演化与缓慢冷却过程。长时间的高温变质作用与缓慢冷却过程也许正是北大别长期难以发现柯石英和有关超高压记录的重要原因。因此,这些成果为大别山三个不同超高压岩片的差异折返模型的建立提供了新的证据。     

Regression diagnostics and robust regression in geothermometer/geobarometer calibration: the garnet-clinopyroxene geothermometer revisited

Abstract The calibration of geothermometers and geobarometers should involve not only the determination of the parameters in the equation used, but also the uncertainties on, and the correlations between, these parameters. This necessitates the use of a technique such as least squares. Given the poor performance of least squares in the presence of outliers in the data, techniques for identifying outliers for exclusion—regression diagnostics, and techniques for handling data which include outliers—robust regression and jackknifing, are essential. These techniques are summarized and their importance is emphasized, and they are applied to the calibration of the garnet-clinopyroxene Fe-Mg exchange geothermometer.
The experimental data of Raheim & Green (1974) and Ellis & Green (1979) are explored using regression diagnostics to discover outliers in the data. After exclusion of the two influential outliers found, a new geothermometer equation for garnet-clinopyroxene Fe-Mg exchange is derived using robust regression and based on all the data: thus, T (K) = 2790 + 10 P + 3140x_ca,g/1.735 + In K _D where T is in Kelvin and P is in kbar. This equation, as might be hoped, is essentially identical to that of Ellis & Green (1979). Equations for calculating the uncertainty in a calculated temperature, contributed by uncertainties in the calibration, are also derived.     

A study on metamorphic temperature and pressure of Archean Miyun Group,Beijing

Metamorphic temperatures and pressures of the Archean Miyun Group were determined from orthopyroxene-clinopyroxene, garnet-clinopyroxene, garnet-biotite and δO _Q ¹⁸ -δO _Mt ¹⁸ geothermometers and orthopyroxene barometer. The results show that the temperature in the first metamorphic stage of the Miyun Group is 820°+50°C and the pressure about 10 kb, which suggests that the granulite facies occurs under moderate pressures with a geothermal gradient of 22°–25°C/km. The corresponding burial depth is about 35 km. The temperature prevailing during the second metamorphic stage is in the range 650°–700°C, indicating a moderate condition between granulite and high grade amphibolite facies. Oxygen isotope data also show that the temperature of later superimposed regression metamorphism of high green schist facies in this area may be within the range of 470°–560°C.     

Mixing behavior in quaternary garnet solid solution and an extended Ellis and Green garnet-clinopyroxene geothermometer

Applications of the different formulations of garnet-clinopyroxene geothermometers to crustal granulites and particularly to metamorphosed manganese formations show a wide scatter of estimates. This is primarily related to the uncertainty associated with the formulation of the activity-composition relationship in garnet solid solution besides that in the clinopyroxene solid solution. The largest uncertainty appears to have resulted from the adopted Mg-Mn interaction parameter in garnet which is yet to be experimentally determined. A fresh statistical regression of eleven pairs of garnet and clinopyroxene equilibrated at identical P-T conditions from a manganese formation in India yields W _Mg-Mn ^Gt (=W_Mn-Mg) value of 1600±500 cal/mols. An extension of the Ellis and Green (1979) geothermometer, known to work well in normal granulites, incorporating terms associated with the activity coefficients of the exchangeable cations, now yields consistent temperature estimation from normal granulites as well as from metamorphosed manganese formations. The geothermometer, however, will not work well for garnet having X_Fe0.65 and clinopyroxenes containing appreciable Na, Al and/or Fe³⁺.     

Garnet-pyroxene and lawsonite-bearing rocks of the Maksyutov Complex (Southern Urals)

The new petrological and geochrological data are used to constrain the nature of garnet-clinopyroxene and lawsonite-bearing rocks, which contain a rare variety of Fe-Ca-rich garnets. These rocks associated with antigorite serpentinites have no equivalents in the other units of the Maksyutov eclogite-glaucophane schist complex and, thus, can be regarded as a separate “mafic-ultramafic” unit. Based on their mineral and chemical composition, the garnet-clinopyroxene and lawsonite-bearing rocks can be interpreted as HP associations formed within a deep continental rifting setting. They experienced a series of metasomatic alterations during decompression exhumation and were accreted to the Maksyutov Complex as a result of the arc-continent collision. The U–Pb zircon data indicate that a Late Riphean–Lower Paleozoic age (824 and 440–470 Ma) was a crystallization age of garnet-clinopyroxene rocks and Ar–Ar white mica age (341 ± 2 Ma) represents the timing of the final accretion of all structural unit to the Maksyutov Complex.     

P–T conditions for the Arendal granulites, southern Norway: implications for the roles of P, T and CO2 in deep crustal LILE-depletion

Abstract The orthopyroxene-clinopyroxene, garnet-orthopyroxene and garnet-clinopyroxene geothermometers, and the garnet-orthopyroxene-plagioclase, garnet-clinopyroxene-plagioclase and anorthite-ferrosilite-grossular-almandine-quartz geobarometers are applied to metabasites and the garnetplagioclase-sillimanite-quartz geobarometer is applied to a metapelite from the Proterozoic Arendal granulite terrain, Bamble sector, Norway. P–T conditions of metamorphism were 7.3 ± 0.5 kbar and 800 ± 60°C.
This terrain shows a regional gradation from the amphibolite facies, into normal LILE content granulite facies rocks and finally strongly LILE deficient granulite facies gneisses. Neither P nor T vary significantly across the entire transition zone. The change in 'grade'parallels the increasing dominance of CO₂ over H₂O in the fluid phase.
LILE-depletion is not a pre-condition of granulite facies metamorphism: granulites may have either 'depleted'or 'normal'chemistries. The results presented herein show that LILE-deficiency in granulite facies orthogneisses is not necessarily related to variations in either P or T . The important mechanisms in the Arendal terrain were (a) direct synmetamorphic crystallization from magma, with primary LILE-poor mineralogies imposed by the prevailing fluid regime, and (b) metamorphic depletion, involving scavenging of LILEs during flushing by mantle-derived CO₂-rich fluids. The latter process is constrained by U–Pb and Rb–Sr isotopic work to have occurred no later than 50 Ma after intrusion of the acid-intermediate gneisses, and was probably associated with contemporary basic magmatism in a tectonic environment similar to a present day cordilleran continental margin.     

Garnet peridotite from Colorado Plateau ultramafic diatremes: Hydrates,carbonates, and comparative geothermometry

Crystal fragments of pyrope from diatremes of ultramafic microbreccia in the Navajo Province of the Colorado Plateau contain inclusions of olivine, pyroxene, spinel, chlorite, amphibole, chlorapatite, and dolomite. The included suite supports earlier hypotheses that hydrous phases and carbonates were primary parts of some garnet peridotite assemblages in the Plateau lithosphere. Garnets with spinel and orthopyroxene inclusions likely all were sampled at pressures less than 36 kb and perhaps as low as 15–20 kb; no evidence was found for inclusions from greater depths. Temperature estimates are 800°–900° C for garnet-clinopyroxene equilibration, but only 500°–700° C for garnetolivine equilibration. Inherent differences between geothermometry methods account for only part of the discrepancy. Pronounced Fe-Mg zoning in garnet at olivine contacts and the lack of such zoning at clinopyroxene contacts are evidence that the difference in part relates to relative reequilibration rates with cooling. The garnet-olivine temperature estimates may be the best approximations to mantle temperatures before eruption. Our data are compatible both with the hypothesis that the garnet peridotite was emplaced in the mantle by large-scale, horizontal transport in the lithosphere and with the hypothesis that rocks were sampled from Precambrian lithosphere cooled to temperatures like those along a low heat flow geotherm. Discordances between the geothermometers here and in other lherzolite localities may be keys to evaluating tectonic histories of lherzolite masses.     

高温地热系统中粘土矿物形成对Na-K和K-Mg地球化学温标准确性的影响

传统地球化学温标在估算高温地热系统内浅层热储温度(一般为100～200℃)时存在局限性,其中应用广泛的Na-K温标和K-Mg温标出现误差的原因仍不清楚.在收集了全球代表性热田内采自地热井的201个流体样品的水文地球化学数据后,利用软件WATCH将井口流体地球化学数据还原为热储条件下的对应值;在此基础上,对Na-K温标和K-Mg温标进行了评价,发现钾长石和常见富钾双八面体粘土矿物均可能对浅层热储内地热流体中的钾含量产生影响,富镁双八面体粘土矿物也可达到与地热流体的平衡,而地热流体中钠含量则受水-岩相互作用的影响很小.因此,浅层地热流体的Na-K比值与热储温度不具有对应关系,而K-Mg温标在计算浅层热储温度时虽然具有一定指示意义,但仍无法得到足够准确的结果.      

Solute and gas geothermometry of geothermal wells: a geochemometrics study for evaluating the effectiveness of geothermometers to predict deep reservoir temperatures

Deep reservoir temperatures of 10 important geothermal systems of the world were estimated by applying 13 solute (Na/K) and 21 gas geothermometers. The predicted temperatures were comprehensively evaluated and compared with measured bottom-hole temperatures using geochemometric techniques. The present study reveals (1) high prediction performances in most of the Na/K geothermometers for the majority of the geothermal fields with liquid-dominated reservoirs, whereas low prediction performances were indicated for the geothermal fields with vapour-dominated and high-temperature reservoirs; (2) the gas geothermometers, in comparison to Na/K, are more successful in predicting the subsurface temperatures in high-temperature geothermal systems; (3) the geothermal systems for which Na/K geothermometers have indicated a high prediction performance, the gas geothermometers have specified a low prediction performances, and vice versa; (4) both Na/K and gas geothermometers, generally, overestimated the reservoir temperatures for the majority of the low-enthalpy geothermal fields and underestimated for the majority of the high-enthalpy geothermal fields; (5) the reservoir temperature predictions of gas geothermometers have more scatter than those temperatures inferred from Na/K geothermometers; and (6) in general, Na/K geothermometers seem to be a more successful geochemical tool in predicting reliable reservoir temperatures than gas geothermometers.     

An evaluation of garnet-clinopyroxene geothermometry in granulites

Four new formulations of the garnet-clinopyroxene geothermometer (Ellis and Green 1979; Ganguly 1979; Saxena 1979; Dahl 1980) have been evaluated in the Adirondacks and five other granulite terranes using results from 94 mineral pairs. The Saxena and Ganguly formulations give temperatures that are generally 100–150° C above those constrained by phase equilibria and other independent thermometry while the empirical calibration of Dahl gives widely scattered, erratic results. Despite some scatter in the data, the Ellis and Green calibration appears to be more accurate and precise than the others and is the most useful garnet-clinopyroxene thermometer currently available for quantitative thermometry in granulites. All four formulations are sensitive to large variations (>70–80° C) in temperature suggesting that problems with accuracy and precision can be improved with further refinement of model-based aspects of the thermometers.     

Solute geothermometry of springs and wells of the Los Azufres and Las Tres Vírgenes geothermal fields,Mexico

Subsurface reservoir temperatures of two important Mexican geothermal systems (Los Azufres and Las Tres Vírgenes) were estimated by applying all available solute geothermometers for 88 and 56 chemical data measurements of the spring waters and fluids of the deep geothermal wells, respectively. Most of the chemical data for spring water of these two geothermal fields are for HCO₃ water, followed by SO₄ and Cl types. For the Los Azufres geothermal field (LAGF), the reservoir temperatures estimated by Na-K geothermometers for springs of HCO₃ and SO₄ waters, and by Na-Li and Li-Mg geothermometers for Cl water, are close to the average bottom-hole temperature (BHT) of the geothermal wells. However, all reservoir temperatures for spring waters from the Las Tres Vírgenes geothermal field (LTVGF) estimated by all solute geothermometers indicated significantly large differences (low temperatures) compared to the BHT. Evaluation of inferred reservoir temperatures for spring waters of the LAGF and LTVGF suggests that not all springs nor all solute geothermometers provide reliable estimation of the reservoir temperatures. Even though chemical equilibrium probably was not achieved in the water–rock system, Na-K geothermometers for HCO₃ water (peripheral water mainly of meteoric origin with little geothermal component) and SO₄ water (geothermal steam heated) and Na-Li and Li-Mg geothermometers for Cl-rich spring water (fully mature geothermal water) of the LAGF indicated reservoir temperatures close to the BHT. However, in comparison with the geothermometry of spring water of the LAGF and LTVGF, fluid measurements from geothermal wells of these two fields indicated reservoir temperatures in close agreement with their respective BHTs. For the best use of the solute geothermometry for spring water, it is advisable to: (1) chemically classify the springs based on water types; (2) identify and eliminate the discordant outlier observations by considering each water type as a separate sampled population; (3) apply all available solute geothermometers employing a suitable computer program such as SolGeo instead of using some specific, arbitrarily chosen geothermometers; and (4) evaluate the temperatures obtained for each solute geothermometer by considering the subsurface lithology, hydrological conditions, and BHTs or static formation temperatures whenever available.