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.     

The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

GPS signal acquisition using the repeatability of successive code phase measurements

This paper introduces a new method for GPS signal acquisition, which is based on the repeatability of successive code phase measurements and the M-of-N search algorithm. The performance of the proposed method in terms of probability of signal detection is similar to that of traditional methods, except that the calculation of the probability of detection does not rely on the noise distribution or the Carrier-to-Noise ratio (C/N₀). The code phase repeatability-based method is presented along with equations for probability of detection and probability of false detection. If the distribution of the noise is known, it also provides an estimate of the C/N₀. The proposed method is illustrated for coherent and non-coherent acquisition and C/N₀ estimation.     

Preservation/exhumation of ultrahigh-pressure subduction complexes

Ultrahigh-pressure (UHP) metamorphic terranes reflect subduction of continental crust to depths of 90–140 km in Phanerozoic contractional orogens. Rocks are intensely overprinted by lower pressure mineral assemblages; traces of relict UHP phases are preserved only under kinetically inhibiting circumstances. Most UHP complexes present in the upper crust are thin, imbricate sheets consisting chiefly of felsic units ± serpentinites; dense mafic and peridotitic rocks make up less than  10% of each exhumed subduction complex. Roundtrip prograde–retrograde P–T paths are completed in 10–20 Myr, and rates of ascent to mid-crustal levels approximate descent velocities. Late-stage domical uplifts typify many UHP complexes.

Sialic crust may be deeply subducted, reflecting profound underflow of an oceanic plate prior to collisional suturing. Exhumation involves decompression through the P–T stability fields of lower pressure metamorphic facies. Scattered UHP relics are retained in strong, refractory, watertight host minerals (e.g., zircon, pyroxene, garnet) typified by low rates of intracrystalline diffusion. Isolation of such inclusions from the recrystallizing rock matrix impedes back reaction. Thin-aspect ratio, ductile-deformed nappes are formed in the subduction zone; heat is conducted away from UHP complexes as they rise along the subduction channel. The low aggregate density of continental crust is much less than that of the mantle it displaces during underflow; its rapid ascent to mid-crustal levels is driven by buoyancy. Return to shallow levels does not require removal of the overlying mantle wedge. Late-stage underplating, structural contraction, tectonic aneurysms and/or plate shallowing convey mid-crustal UHP décollements surfaceward in domical uplifts where they are exposed by erosion. Unless these situations are mutually satisfied, UHP complexes are completely transformed to low-pressure assemblages, obliterating all evidence of profound subduction.     

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     

An automated strategy for calculation of phase diagram sections and retrieval of rock properties as a function of physical conditions

We formulate an algorithm for the calculation of stable phase relations of a system with constrained bulk composition as a function of its environmental variables. The basis of this algorithm is the approximate representation of the free energy composition surfaces of solution phases by inscribed polyhedra. This representation leads to discretization of high variance phase fields into a continuous mesh of smaller polygonal fields within which the composition and physical properties of the phases are uniquely determined. The resulting phase diagram sections are useful for understanding the phase relations of complex metamorphic systems and for applications in which it is necessary to establish the variations in rock properties such as density, seismic velocities and volatile‐content through a metamorphic cycle. The algorithm has been implemented within a computer program that is general with respect to both the choice of variables and the number of components and phases possible in a system, and is independent of the structure of the equations of state used to describe the phases of the system.     

The geological and geodynamic evolution of the eastern Black Sea basin: insights from 2-D and 3-D tectonic modelling

Subsidence mechanisms that may have controlled the evolution of the eastern Black Sea have been studied and simulated using a numerical model that integrates structural, thermal, isostatic and surface processes in both two- (2-D) and three-dimensions (3-D). The model enables the forward modelling of extensional basin evolution followed by deformation due to subsequent extensional and compressional events. Seismic data show that the eastern Black Sea has evolved via a sequence of interrelated tectonic events that began with early Tertiary rifting followed by several phases of compression, mainly confined to the edges of the basin. A large magnitude (approximately 12 km) of regional subsidence also occurred in the central basin throughout the Tertiary. Models that simulate the magnitude of observed fault controlled extension (β=1.13) do not reproduce the total depth of the basin. Similarly, the modelling of compressional deformation around the edges of the basin does little to enhance subsidence in the central basin. A modelling approach that quantifies lithosphere extension according to the amount of observed crustal thinning and thickening across the basin provides the closest match to overall subsidence. The modelling also shows that deep crustal and mantle–lithosphere processes can significantly influence the rate and magnitude of syn- to post-rift subsidence and shows that such mechanisms may have played an important role in forming the anomalously thin syn-rift and thick Miocene–Quaternary sequences observed in the basin. It is also suggested that extension of a 40–45 km thick pre-rift crust is required to generate the observed magnitude of total subsidence when considering a realistic bathymetry.     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.     

Structure sismique du socle paléozoı̈que du bassin des Doukkala,Môle côtier,Maroc occidental. Indication en faveur de l'existence d'une phase éo-varisqueSeismic structure of the Doukkala basin,Palaeozoic basement,western Morocco: a hint for an Eovariscan fold-and-thrust belt

Seismic profiles and well data from the Doukkala basin unravel the structure of the Palaeozoic basement and suggest that this coastal zone of western Morocco was affected by a compressive phase during the Frasnian. This resulted in the formation of upright, plurikilometric folds associated with reverse faults (North Doukkala), and of asymmetrical folds associated with mostly west verging ramps (South Doukkala). Folding involved all pre-Upper Frasnian formations and caused partial or total hiatus of Upper Frasnian–Strunian strata. This event can be correlated with the orogenic phase reported from more internal domains of the Morocco Hercynian belt, where it is referred to as the ‘Bretonne’ or ‘Eovariscan’ phase. To cite this article: H. Echarfaoui et al., C. R. Geoscience 334 (2002) 13–20