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.     

ANALYTICAL HAYAMI SOLUTION FOR THE DIFFUSIVE WAVE FLOOD ROUTING PROBLEM WITH LATERAL INFLOW

The diffusive wave equation is generally used in flood routing in rivers. The two parameters of the equation, celerity and diffusivity, are usually taken as functions of the discharge. If these two parameters can be assumed to be constant without lateral inflow, the diffusive wave equation may have an analytical solution: the Hayami model. A general analytical method, based on ‘Hayami’s hypothesis, is developed here which resolves the diffusive wave flood routing equation with lateral inflow or outflow uniformly distributed over a channel reach. Flood routing parameters are then identified using observed inflow and outflow and the Hayami model used to simulate outflow. Two examples are discussed. Firstly, the prediction of the hydrograph at a downstream section on the basis of a knowledge of the hydrograph at an upstream section and the lateral inflow. The second example concerns lateral inflow identification between an upstream and a downstream section on the basis of a knowledge of hydrographs at the upstream and downstream sections. The new general Hayami model was applied to flood routing simulation and for lateral inflow identification of the River Allier in France. The major advantages of the method relate to computer simulation, real-time forecasting and control applications in examples where numerical instabilities, in the solution of the partial differential equations must be avoided.     

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.