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.     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

Composition and pre-metamorphic geodynamic setting of the ultrahigh-pressure metabasic rocks from Dabie Shan,E-China

As part of the Yangtze plate, segments of the Dabie Shan terrane of Central China underwent ultra-high pressure metamorphism during Triassic subduction. We studied the geochemistry of the abundant eclogites to evaluate the nature of the protoliths and their geodynamic setting. Although some previous geochemical work exists, the analyses and interpretation herein are based on a new subdivision of the ultra-high pressure sequence into basement and cover units (Changpu and Ganghe Unit), revealing new and important results. In addition, eclogites of the so-called HP Unit south of the UHP units were studied. Whereas the large ion lithophile elements indicate postmagmatic, metasomatic changes of some samples, the high-field strength elements and the rare earth elements display original magmatic trends. The geochemical characteristics of the eclogites of the ultra-high pressure areas display a strong dependence on their “structural” and geographic position. The eclogites of the basement and the Changpu Unit indicate melt intrusion and extrusion in a continental rift system, i.e. during extensional tectonics. In contrast, the Ganghe Unit is characterized by a pronounced chemical homogeneity. The composition of the eclogites indicates generation from a mantle source highly influenced by slab-derived fluids. Those of the HP Unit show similar characteristics. Magmatism of the Ganghe and HP Unit probably occurred in a continental arc setting. A similar age for both units, geographically and/or tectonically separated, is possible. The geodynamic interpretation based on the geochemistry of the four units points to a Neoproterozoic scenario in which the protoliths of the HP and the cover units could have been of similar age and deposited in one evolving geological system. A rift-related larger-scale basin might have formed, e.g. a continental back-arc basin behind a magmatic arc after or simultaneous to sedimentation and magmatism in the magmatic arc. Alternatively, magmatism occurred in independent geodynamic settings, distinct in time and space. The units were juxtaposed during exhumation, after subduction to varying depths.     

黑龙江省鸡西盆地南、北两带城子河组精细对比

鸡西盆地被平麻断裂分割成南、北两个条带,多数研究者将南带1煤之下的海相层叠伏于北带海相层或底砾岩之下,称为“城子河组下段”或“石河北组”。利用事件地层学的原理,分析盆地内城子河组早期盆地裂陷、海侵、火山喷发、聚煤等重大事件在地层精细对比方面的意义,甄别出了盆地中重大事件的地层记录．以等时层为基础建立对比框架。通过大量钻井资料追溯对比,认为盆地中坳陷部位的渴湖海湾沉积是南北两带的沉积过渡类型,根据微体古生物研究,南北两带海相层中的沟鞭藻类化石面貌基本一致,南北带的海相层应是同一层位,初步实现了盆地内南、北两个条带城子河组的精细对比,为恢复盆地古地理环境奠定了可靠的基础。     

地铁施工对济南白泉泉群流量的影响

为研究地铁建设对济南白泉泉群的影响，在综合分析白泉泉域地质、水文地质条件的基础上，假定研究区岩溶强径流带位置及水力性质，利用FEFLOW软件建立地下水流数值模型。以规划地铁M1号线为研究对象，分析了济南东站、梁王站、梁王东站分别施工及3个站同时施工4种情景下，采用施工降水或施工降水+人工回灌两种施工方式对白泉泉群流量的影响。结果表明：单独采用施工降水的施工方式使得白泉泉群流量衰减，其中3个站同时施工对泉流量的影响最大，泉流量最大衰减达5.48%；各站分别施工时，济南东站对泉流量影响最大，泉流量较未施工时减少了0.043×10⁴ m³/d。采用施工降水+人工回灌的施工方式，能够有效缓解泉流量的衰减，各车站施工时的泉流量衰减由仅施工降水时的2.26%~5.48%降低至0.08%~1.21%。岩溶强径流带有利于地下水形成优势径流，促进白泉泉群补给，一定程度上缓解因地铁施工引起的泉流量衰减。     

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.     

Evolution from Oceanic Subduction to Continental Collision: a Case Study from the Northern Tibetan Plateau Based on Geochemical and Geochronological Data

Two apparently distinct, sub-parallel, paleo-subduction zonescan be recognized along the northern margin of the Tibetan Plateau:the North Qilian Suture Zone (oceanic-type) with ophioliticmélanges and high-pressure eclogites and blueschistsin the north, and the North Qaidam Belt (continental-type) inthe south, an ultrahigh-pressure (UHP) metamorphic terrane comprisingpelitic and granitic gneisses, eclogites and garnet peridotites.Eclogites from both belts have protoliths broadly similar tomid-ocean ridge basalts (MORB) or oceanic island basalts (OIB)in composition with overlapping metamorphic ages (480–440Ma, with weighted mean ages of 464 ± 6 Ma for North Qilianand 457 ± 7 Ma for North Qaidam), determined by zirconU–Pb sensitive high-resolution ion microprobe dating.Coesite-bearing zircon grains in pelitic gneisses from the NorthQaidam UHP Belt yield a peak metamorphic age of 423 ±6 Ma, 40 Myr younger than the age of eclogite formation, anda retrograde age of 403 ± 9 Ma. These data, combinedwith regional relationships, allow us to infer that these twoparallel belts may represent an evolutionary sequence from oceanicsubduction to continental collision, and continental underthrusting,to final exhumation. The Qilian–Qaidam Craton was probablya fragment of the Rodinia supercontinent with a passive marginand extended oceanic lithosphere in the north, which was subductedbeneath the North China Craton to depths >100 km at c. 423Ma and exhumed at c. 403 Ma (zircon rim ages in pelitic gneiss). KEY WORDS: HP and UHP rocks; subduction belts; zircon SHRIMP ages; Northern Tibetan Plateau     

Topography as a major factor in the development of arcuate thrust belts: insights from sandbox experiments

We have used sandbox experiments to investigate and to illustrate the effects of topography upon the development of arcuate thrust belts. In experiments where a sand pack shortened and thickened in front of an advancing rectilinear piston, the geometry of the developing thrust wedge was highly sensitive to variations in surface topography. In the absence of erosion and sedimentation, the surface slope tended to become uniform, as predicted by the theory of critical taper. Under these conditions, the wedge propagated by sequential accretion of new thrust slices. In contrast, where erosion or sedimentation caused the topographic profile to become irregular, thrusts developed out of sequence. For example, erosion throughout a hinterland caused underlying thrusts to remain active and inhibited the development of new thrusts in the foreland. Where initial topography was irregular in plan view, accreting thrusts tended to be arcuate. They were convex towards the foreland, around an initially high area; concave towards the foreland, around an initially low area. Initial plateaux tended to behave rigidly, while arcuate thrust slices accreted to them. Thrust motions were radial with respect to each plateau. Within transfer zones to each side, fault blocks rotated about vertical axes and thrust motions were oblique-slip. At late stages of deformation, the surface slope of the thrust wedge tended towards a uniform value. Initial mountains of conical shape (representing volcanoes) also escaped deformation, except at depth, where they detached. Arcuate thrust slices accreted to front and back. Where a developing thrust wedge was subject to local incision, accreting thrust slices dipped towards surrounding areas of high topography, forming Vs across valleys.Arcuate structural patterns are to be found around the three highest plateaux on Earth (Tibet, Pamirs and Altiplano) and around the Tromen volcanic ridge in the Neuquén Basin of northern Patagonia. We infer that these areas behaved in quasi-rigid fashion, protected as they were by their high topography.     

Crustal structure of the Ruby Mountains and southern Carlin Trend region, Nevada, from magnetotelluric data

We have collected about 150 magnetotelluric (MT) soundings in northeastern Nevada in the region of the Ruby Mountains metamorphic core complex uplift and southern Carlin mineral trend, in an effort to illuminate controls on core complex evolution and deposition of world-class gold deposits. The region has experienced a broad range of tectonic events including several periods of compressional and extensional deformation, which have contributed to the total expression of electrical resistivity. Most of the soundings reside in three east–west profiles across increasing degrees of core uplift to the north (Bald Mountain, Harrison Pass, and Secret Pass latitudes). One short cross-line was also taken to assess an east–west structure to the north of the northern profile. Model resistivity cross-sections were derived from the MT data using a 2-D inversion algorithm, which damps departures of model parameters from an a priori structure. Geological interpretation of the resistivity combines previous seismic, potential field and isotope models, structural and petrological models for regional compression and extension, and detailed structural/stratigraphic interpretations incorporating drilling for petroleum and mineral exploration. To first order, the resistivity structure is one of a moderately conductive, Phanerozoic sedimentary section fundamentally disrupted by intrusion and uplift of resistive crystalline rocks. Late Devonian and early Mississippian shales of the Pilot and Chainman Formations together form an important conductive marker sequence in the stratigraphy and show pronounced increases in conductance (conductivity–thickness product) from east to west. These increases are attributed to graphitization caused by Elko–Sevier era compressional shear deformation and possibly by intrusive heating. The resistive crystalline central massifs adjoin the host stratigraphy across crustal-scale, steeply dipping fault zones. The zones provide pathways to the lower crust for heterogeneous, upper crustal induced, electric current flow. Resistive core complex crust appears steeply bounded under the middle of the neighboring grabens and not to deepen at a shallow angle to arbitrary distances to the west. The numerous crustal breaks imaged with MT may contribute to the low effective elastic thickness (T_e) estimated regionally for the Great Basin and exemplify the mid-crustal, steeply dipping slip zones in which major earthquakes nucleate. An east–west oriented conductor in the crystalline upper crust spans the East Humboldt Range and northern Ruby Mountains. The conductor may be related to nearby graphitic metasediments, with possible alteration by middle Tertiary magmatism. Lower crustal resistivity everywhere under the profiles is low and appears quasi one-dimensional. It is consistent with a low rock porosity (<1 vol.%) containing hypersaline brines and possible water-undersaturated crustal melts, residual to the mostly Miocene regional extension. The resistivity expression of the southern Carlin Trend (CT) in the Pinon Range is not a simple lineament but rather a family of structures attributed to Eocene intrusion, stratal deformation, and alteration/graphitization. Substantial reactivation or overprinting by core complex uplift or Basin–Range extensional events seems likely. We concur with others that the Carlin Trend may result in part from overlap of the large Eocene Northeast Nevada Volcanic Field with Precambrian–Paleozoic deep-water clastic source rocks thickening abruptly to the west of the Pinon Range, and projecting to the north–northwest.     

Models of passive margin inversion: implications for the Rhenohercynian fold-and-thrust belt, Belgium and Germany

Positive tectonic inversion is related to the transmission of compressional stresses along a décollement into the foreland of an orogenic zone. This stress and strain concentration in regions remote from the main orogenic front is commonly related to the presence of pre-existing rheological heterogeneities such as normal syn-depositional faults. During inversion, these pre-existing normal faults are reactivated as reverse faults. Tectonic inversion in the Rhenohercynian fold-and-thrust belt during the Variscan Orogeny shows that inversion is likely synchronous with the onset of collision in the hinterland. Here, we present the results of a simplified thermo-mechanical model (STM) which allows one to study strain partitioning between two orogenic zones. We show that, if the two orogenic zones have the same mechanical properties, the viscosity of the décollement, which links them, controls the initial strain partitioning. During subsequent finite shortening, erosional processes determine the partitioning of strain rate. The presence of a weak structure in the inverted zone and of a low-viscosity décollement leads to initial strain concentration in the inverted track rather than in the collision zone and a progressive decrease in strain partitioning between the two orogenic zones. The STM results are in good agreement with results of a 2D finite-element model. We conclude that, in the western part of the Rhenohercynian Massif, simultaneous uplift and deformation within the Mid-German Crystalline Rise (the main collision zone) and the Ardenne Anticlinorium (the inverted zone) lead to interpreting this orogenic event as a case of vice tectonic rather than the propagation of a ‘wave of folding’ towards the Variscan front, as suggested by previous authors.