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     

SPRING ICE-JAM FLOODING OF THE PEACE-ATHABASCA DELTA: EVIDENCE OF A CLIMATIC OSCILLATION

A historic record of spring ice-jam floods of the Peace-Athabasca Delta was analyzed for the years 1826–1995. The temporal pattern of flooding is non-random. The likelihood of a flood following a flood, or a non-flood following a non-flood, is greater than expected by chance. Probability analysis of flood occurrence reveals that the period 1860–1880 was a time of unusually few floods, and the period 1915–1950 was a time of unusually frequent floods. The long-term flood frequency is 1 flood in 6.25 years. Changes in flood frequency over the record reveal a pattern of oscillation described by a sine-based model that is correlated with the long-term (Gleissberg) cycle of solar activity. Monte Carlo simulation was used to test a Bennett Dam Model and a Cyclic Model. The Bennett Dam Model is unlikely to have generated the observed flood history (p=0.04). The observed flood history shows a better fit to the Cyclic Model (p=0.65). No correlations between floods and ENSO cold or warm events was detected. The most recent wet period began about 1900 and ended in the early 1960's prior to completion of the W. A. C. Bennett Dam in British Columbia. As independent corroboration of climatically-driven changes in flood frequency we present three additional lines of evidence. The pattern of annual muskrat returns (95 year record) reveals both 10 year cycles and long-term patterns that agree well with the observed flood cycle. The annual area burned in Wood Buffalo National Park is inversely related to flood occurrence. Incised channels and dendritic drainage patterns in the bed of Lake Mamawi provide probable evidence of a previous dry period in the delta. Climatic change or oscillation likely underlies the drying trend observed in recent decades in the Peace-Athabasca Delta.     

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.     

Revised activity–composition models for clinopyroxene and amphibole

Recent activity–composition models for clinopyroxene and amphibole are revised to provide better consistency with observed phase relations in natural rocks. For clinopyroxene, the calibration in NCFMAS is retained, but the incorporation of acmite is revised to improve the partitioning of ferric iron between coexisting clinopyroxenes. For amphibole, the NCFMASH calibration is retained, but the addition of ferric iron is changed to provide consistency with the clinopyroxenes. The thermodynamics of orthoamphibole (gedrite) is also adjusted to resolve an unrelated inconsistency. The effects of these improvements are illustrated through comparison of calculated pseudosections produced with the existing and new models with natural data from lawsonite eclogites.     

Low-pressure granulite facies metamorphism in the Larsemann Hills area, East Antarctica; petrology and tectonic implications for the evolution of the Prydz Bay area

ABSTRACT Thermobarometric studies on various granulite facies areas along the Prydz Bay coast, East Antarctica (73°-79°E, 68°-70°S), show that, at around 1100 Ma, during a late Proterozoic orogeny, the rocks of the Larsemann Hills suffered a lower pressure metamorphic peak than the surrounding areas. Along the Prydz Bay coast, the rocks affected by this event include parts of the Vestfold Hills block plus all of the Rauer Group, the Larsemann Hills and the Munro Kerr Mountains. The dykes in the south-west corner of the Vestfold Hills were recrystallized during this event with little deformation at temperatures not quite as high as in the areas further south-west (650°C, 6.5 kbar) (Collerson et al., 1983), the Rauer Group was metamorphosed at 800°C and 7.5 kbar (Harley, 1987a), the Larsemann Hills at 750°C and 4.5 kbar, and the Munro Kerr Mountains probably at around 850°C and 5 kbar. Retrograde equilibration in the different areas occurred during decompression to about 10 km depth in all areas, followed by isobaric cooling at this depth. This paper shows that the peak metamorphism in the Larsemann Hills occurred at a pressure which is too low to have been the consequence of thermal relaxation of overthickened crust with normal mantle heat flow. Although other areas in Prydz Bay were metamorphosed at sufficiently high pressures so that their decompression paths are not inconsistent with a continental collision model, the inferred pre-metamorphic peak histories and the requirement of consistency with the Larsemann Hills, make it unlikely that collision followed by erosion-driven decompression is an appropriate model. We suggest that the thermal regime of the crust in the Larsemann Hills region was controlled by a perturbation in the asthenosphere, with magma invasion of the crust. We suggest that the 500 Ma event, represented in Prydz Bay by granitic outcrops at Landing Bluff and by several K/Ar ages from the Larsemann Hills area, was responsible for the final excavation of the terrane.     

Decompressional coronas and symplectites in granulites of the Musgrave Complex, central Australia

Abstract In granulite facies metapelitic rocks in the Musgrave Complex, central Australia, reaction between S1 garnet and sillimanite involves the development in S2 of both garnet + cordierite + hercynitic spinel + biotite and hercynitic spinel + cordierite + sillimanite + biotite. The S2 assemblages occur either in coronas and symplectites, mainly around garnet, or, in rocks in which S2 is more strongly developed, as recrystallized assemblages. Ignoring the presence of biotite and ilmenite, the mineral textures can be accounted for qualitatively by a consideration of the model system FeO-MgO-Al₂O₃-SiO₂ (FMAS); the textural relationships accord with decompression accompanying the change from S1 to S2. However, since biotite and ilmenite are involved in the assemblages, the parageneses are better accounted for in terms of equilibria in the expanded model system K₂O-FeO-MgO-Al₂O₃-SiO₂-H₂-TiO₂-Fe₂O₃ (KFMASHTO), i.e. AFM + TiO₂+ Fe₂O₃. The coronas reflect the tectonic unroofing of at least part of the Musgrave Complex from peak S1 conditions of about 8 kbar to S2 conditions of about 4 kbar.     

Calculating phase diagrams involving solid solutions via non-linear equations, with examples using THERMOCALC

Phase diagrams involving solid solutions are calculated by solving sets of non-linear equations. In calculating P–T  projections and compatibility diagrams, the equations used for each equilibrium are the equilibrium relationships for an independent set of reactions between the end-members of the phases in the equilibrium. Invariant points and univariant lines in P–T  projections can be calculated directly, as can coordinates in compatibility diagrams. In calculating P–T  and T–x / P–x pseudosections – diagrams drawn for particular bulk compositions – the equilibrium relationship equations are augmented by mass balance equations. Lines in pseudosections, where the mode of one phase in the lower variance equilibrium is zero, and points, where the modes of two phases are zero, can then be calculated directly. The software, THERMOCALC, allows the calculation of these and a range of other types of phase diagram. Examples of phase diagrams and phase diagram movies, with instructions for their production, along with the THERMOCALC input and output files, and the Mathematica^TM functions for assembling them, are presented in this paper, partly in hard copy and partly on the JMG web sites (http://www.gly.bris.ac.uk/www/jmg/jmg.html, or equivalent Australian or USA sites).