Metamorphism in Archaean greenstone belts: calculated fluid compositions and implications for gold mineralization

An inescapable consequence of the metamorphism of greenstone belt sequences is the release of a large volume of metamorphic fluid of low salinity with chemical characteristics controlled by the mineral assemblages involved in the devolatilization reactions. For mafic and ultramafic sequences, the composition of fluids released at upper greenschist to lower amphibolite facies conditions for the necessary relatively hot geotherm corresponds to those inferred for greenstone gold deposits (X_CO2= 0.2–0.3). This result follows from the calculation of mineral equilibria in the model system CaO–MgO–FeO–Al₂O₃–SiO₂–H₂O–CO₂, using a new, expanded, internally consistent dataset. Greenstone metamorphism cannot have involved much crustal over-thickening, because very shallow levels of greenstone belts are preserved. Such orogeny can be accounted for if compressive deformation of the crust is accompanied by thinning of the mantle lithosphere. In this case, the observed metamorphism, which was contemporaneous with deformation, is of the low-P high-T type. For this type of metamorphism, the metamorphic peak should have occurred earlier at deeper levels in the crust; i.e. the piezothermal array should be of the ‘deeper-earlier’type. However, at shallow crustal levels, the piezothermal array is likely to have been of ‘deeper-later’type, as a consequence of erosion. Thus, while the lower crust reached maximum temperatures, and partially melted to produce the observed granites, mid-crustal levels were releasing fluids prograde into shallow crustal levels that were already retrograde. We propose that these fluids are responsible for the gold mineralization. Thus, the contemporaneity of igneous activity and gold mineralization is a natural consequence of the thermal evolution, and does not mean that the mineralization has to be a consequence of igneous processes. Upward migration of metamorphic fluid, via appropriate structurally controlled pathways, will bring the fluid into contact with mineral assemblages that have equilibrated with a fluid with significantly lower X_CO2. These assemblages are therefore grossly out of equilibrium with the fluid. In the case of infiltrated metabasic rocks, intense carbonation and sulphidation is predicted. If, as seems reasonable, gold is mobilized by the fluid generated by devolatilization, then the combination of processes proposed, most of which are an inevitable consequence of the metamorphism, leads to the formation of greenstone gold deposits predominantly from metamorphic fluids.     

Calculated Phase Relations in the System Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O with Applications to UHP Eclogites and Whiteschists

Pressure–temperature grids in the system Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O and its subsystems have been calculatedin the range 15–45 kbar and 550–900°C, usingan internally consistent thermodynamic dataset and new thermodynamicmodels for amphibole, white mica, and clinopyroxene, with thesoftware THERMOCALC. Minerals considered for the grids includegarnet, omphacite, diopside, jadeite, hornblende, actinolite,glaucophane, zoisite, lawsonite, kyanite, coesite, quartz, talc,muscovite, paragonite, biotite, chlorite, and plagioclase. Compatibilitydiagrams are used to illustrate the phase relationships in thegrids. Coesite-bearing eclogites and a whiteschist from Chinaare used to demonstrate the ability of pseudosections to modelphase relationships in natural ultrahigh-pressure metamorphicrocks. Under water-saturated conditions, chlorite-bearing assemblagesin Mg- and Al-rich eclogites are stable at lower temperaturesthan in Fe-rich eclogites. The relative temperature stabilityof the three amphiboles is hornblende > actinolite > glaucophane(amphibole names used sensu lato). Talc-bearing assemblagesare stable only at low temperature and high pressure in Mg-and Al-rich eclogites. For most eclogite compositions, talccoexists with lawsonite, but not zoisite, in the stability fieldof coesite. Water content contouring of pressure–temperaturepseudosections, along with appropriate geotherms, provides newconstraints concerning dehydration of such rocks in subductingslabs. Chlorite and lawsonite are two important H₂O-carriersin subducting slabs. Depending on bulk composition and pressure–temperaturepath, amphibole may or may not be a major H₂O-carrier to depth.In most cases, dehydration to make ultrahigh-pressure eclogitestakes place gradually, with H₂O content controlled by divariantor higher variance assemblages. Therefore, fluid fluxes in subductionzones are likely to be continuous, with the rate of dehydrationchanging with changing pressure and temperature. Further, eclogitesof different bulk compositions dehydrate differently. Dehydrationof Fe-rich eclogite is nearly complete at relatively shallowdepth, whereas Mg- and Al-rich eclogites dehydrate continuouslydown to greater depth. KEY WORDS: dehydration; eclogites; phase relations; THERMOCALC; UHP metamorphism; whiteschists     

Calculated Phase Relations in the System NCKFMASH (Na2O-CaO-K2O-FeO-MgO-Al2O3-SiO2-H2O) for High-Pressure Metapelites

Petrogenetic grids in the system NCKFMASH (Na₂O–CaO–K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O)and the subsystems NCKMASH and NCKFASH calculated with the softwareTHERMOCALC 3.1 are presented for the P–T range 7–30kbar and 450–680°C, for assemblages involving garnet,chloritoid, biotite, carpholite, talc, chlorite, kyanite, staurolite,paragonite, glaucophane, jadeite, omphacite, diopsidic pyroxene,plagioclase, zoisite and lawsonite, with phengite, quartz/coesiteand H₂O in excess. These grids, together with calculated compatibilitydiagrams and P–T and T–X_Ca and P–X_Ca pseudosectionsfor different bulk-rock compositions, show that incorporationof Ca into the NKFMASH system leads to many of the NKFMASH invariantequilibria moving to lower pressure and/or lower temperature,which results, in most cases, in the stability of jadeite andgarnet being enlarged, but in the reduction of stability ofglaucophane, plagioclase and AFM phases. The effect of Ca onthe stability of paragonite is dependent on mineral assemblageat different P–T conditions. The calculated NCKFMASH diagramsare powerful in delineating the phase equilibria and P–Tconditions of natural pelitic assemblages. Moreover, contoursof the calculated phengite Si isopleths in P–T and P–X_Capseudosections confirm that phengite barometry in NCKFMASH isstrongly dependent on mineral assemblage. KEY WORDS: phase relations; metapelites; NCKFMASH; THERMOCALC; phengite geobarometry     

Internal architecture and sedimentary evolution of coarse-grained, turbidite channel-levee complexes, Early Eocene Regência Canyon, Espírito Santo Basin, Brazil

Successions of Early Eocene coarse-grained turbidites up to 400 m thick fill fault-controlled canyons along the eastern Brazilian continental margin. They form part of a Late Albian to Early Eocene transgressive succession characterized by onlapping, deepening-upward sedimentation. In the Lagoa Parda oil field (Regência Canyon, Espírito Santo Basin) the turbidite facies consist mostly of unstratified conglomerate and sandstone, with interbedded bioturbated mudstone and thin-bedded, stratified sandstone. Within the main Regência Canyon, the coarser grained facies occur within 38 deeply incised channels. The fills are 9 to >50 m thick, 210 to >1050 m wide and >1 km long. The finer grained facies build asymmetrical levees that are higher and thicker on the left side (looking downstream) of their channels, probably as an effect of the Coriolis force (to the left in the Southern Hemisphere). Nine levee successions up to 50 m thick are associated with the 20 youngest channels. The deposits filling the low-sinuosity Lagoa Parda channels record successive channel abandonment through relatively rapid avulsions. Avulsions of unleveed channels took place randomly, but channels with well-developed levees show preferential avulsion to the right (looking downstream), opposite to the direction of preferential levee growth. Lagoa Parda channels can be grouped into three complexes 20–100 m thick. These complexes have an estimated duration of about 140 000 years. It is suggested that control of the development of individual channel complexes was related to variation in sediment supply, in turn probably related to climatic changes. The deposition of each channel complex would have followed an increase in sediment supply into the Regência Canyon through delta/fan-delta and littoral drift systems, which in turn would have responded to phases of higher denudation rates in the high-relief, ancestral coastal ranges of south-eastern Brazil. Overall, the three Lagoa Parda channel complexes form a turbidite succession characterized by channel fills that become narrower, thinner and finer grained upward. These trends were induced mostly by a longer term (>400 000 years) decrease in sediment supply, which in turn resulted from the combined effects of a long-term (second-order) trend of sea-level rise, and the decreasing fault activity at the basin margin and source area.     

Partial melting during tectonic exhumation of a granulite terrane: an example from the Larsemann Hills, East Antarctica

Anatectic migmatites in medium- to low-pressure granulite facies metasediments exposed in the Larsemann Hills, East Antarctica, contain leucosomes with abundant quartz and plagioclase and minor interstitial K-feldspar, and assemblages of garnet–cordierite–spinel–ilmenite–sillimanite. Qualitative modelling in the system K₂O–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–O₂, in conjunction with various P–T calculations indicate that the high-grade retrograde evolution of the terrane was dominated by decompression from peak conditions of c. 7 kbar at c. 800 °C to 4–5 kbar at c. 750 °C. Extensive partial melting during decompression involved the replacement of biotite by the assemblage cordierite–garnet–spinel within the leucosomes. These leucosomes represent the site of partial melt generation, the cordierite–garnet–spinel–ilmenite assemblage representing the solid products and excess reactants from the melting reaction. The extraction and accumulation of this decompression-generated melt led to the formation of syntectonic pegmatites and extensive granitic plutons. Leucosome development and terrane decompression proceeded during crustal transpression, synchronous with upper crustal extension, during a progressive Early Palaeozoic collisional event. Subsequent retrograde evolution was characterized by cooling, as indicated by the growth of biotite replacing spinel and garnet, thin mantles of cordierite replacing spinel and quartz within metapelites, and garnet replacing orthopyroxene and hornblende within metabasites. P–T calculations on late mylonites indicate lower grade conditions of formation of c. 3.5 kbar at c. 650 °C, consistent with the development of late cooling textures.     

The pressure dependence of the zirconium-in-rutile thermometer

The solubility of ZrO₂ in rutile is strongly temperature-dependent and has been identified as a potentially powerful thermometer when the rutile coexists with an appropriate buffer assemblage, e.g. zircon + quartz. In combination with experimental data at 10 kbar, previous consideration of data on natural rutile has not identified a pressure dependence for the thermometer. However, the expected volume change as a result of substitution of the larger Zr⁴⁺ cation for Ti⁴⁺ suggests that the Zr content of rutile should decrease with increasing pressure. To investigate the pressure dependence of the thermometer, piston cylinder (at 10, 20 & 30 kbar) and 1 atm furnace experiments were performed in the system ZrO₂-TiO₂-SiO₂. The solubility of ZrO₂ in rutile, in the presence of zircon and quartz was reversed at each pressure value. From these experiments, the thermodynamics of the end-member reaction ZrSiO₄ = SiO₂ + ZrO₂ (in rutile) have been determined. There is a secondary pressure effect accompanying the primary temperature dependence of the Zr content of rutile. New thermometer equations are, in the α -quartz field: in the β -quartz field and in the coesite field in which φ is ppm Zr, P is in kbar and R is the gas constant, 0.0083144 kJ K⁻¹. Thermometric results using these equations are shown for a range of geological settings.