The Niagara Escarpment is a geological feature located in southern Ontario, Canada, and the northeastern United States, comprising highly fractured sandstone, shale and carbonates deposited during the Ordovician and Silurian periods. Differential erosion of the strata has generated a steep cliff face which bisects the city of Hamilton, Ontario. Geological fractures are widespread in the escarpment and result in the formation of unstable blocks of rock subject to erosion through rockfall. This presents structural stability issues of concern due to the proximity of the escarpment to urban infrastructure. We quantify and analyse fracture networks in the escarpment using a combined field- and numerical-modelling-based approach. The location, orientation and aperture of fractures were documented from local outcrops using scanline and area survey methods. Clusters of poles describing the orientation of geological discontinuities were identified in spherical projections, defining three sets: (1) a sub-vertical stratabound set striking N–S, (2) a sub-vertical stratabound set striking E–W, and (3) a set parallel to horizontal sedimentary bedding planes which we infer controlled sub-vertical fracture geometry. Discrete fracture network modelling of fracture sets highlights their high degree of connectivity, and contribution to local geohazards, and quantifies their role in controlling fluid flow through escarpment strata, which is dependent on fracture aperture. Additionally, bedding planes have the potential to act as free surfaces, facilitating stress conditions in which approximately cuboid blocks are produced, and increasing the risk of rockfalls. We conclude that fractures present a first-order control on the fluid-flow properties and stability of the Niagara Escarpment.

U–Pb geochronology, Hf isotopes and trace-element chemistry of zircon grains from migmatite of the upper Sutlej valley (Leo Pargil), Northwest Himalaya, reveal a protracted geological evolution and constrain anatexis and tectonothermal processes in response to Himalayan orogenesis. U–Pb geochronology and ϵHf record separate clusters of ages on the concordia plots in the migmatite (1050–950 Ma, 850–790 Ma and 650–500 Ma). The 1050–950 Ma zircon population supports a provenance from magmatic units related to the assembly of Rodinia. A minor amount of Palaeoproterozoic grains were likely derived from the Indian craton. The potential source rock of the 930–800 Ma detrital zircons may be granitoid present in Greater Himalayan rocks themselves and the Aravalli Range, which has 870–800 Ma granitic rocks. The arc-type basement within the Himalayan–Tibet orogen recorded (900–600 Ma) igneous activity, which may depict a northeasterly extension of juvenile terranes in the Arabian–Nubian Shield. The granitoid of 800 Ma may be a potential source for 790 Ma detrital zircons owing to scatter in 206/238 dates. The 650–500 Ma zircon population suggests their derivation from the East African Orogen and Ross–Delamerian Orogen of Gondwana. The Cambrian–Ordovician magmatism during the Bhimphedian Orogeny and observed late Neoproterozoic to Ordovician detrital zircons have been derived to some extent from Greater Himalayan magmatic sources. We found no detrital zircon grains that cannot be explained as coming from local sources. One sample yielded a discordia lower intercept age of 15.6 ± 2.2 Ma, the age of melt crystallization.

Temporal constraints on the closure of the eastern segment of the Palaeo-Asian Ocean along the northern margin of the North China Craton (NCC) remain unclear. As a part of the NCC, the sedimentation and tectonic evolution of the Late Palaeozoic Ordos Basin were closely related to the opening and closing of the Palaeo-Asian Ocean. We use petrology, quantitative mineralogical analysis, U–Pb geochronology and trace element signatures of detrital zircons of the Lower Shihezi Formation from two sections in the eastern north Ordos Basin and two sections in the western north Ordos Basin to reconstruct the sedimentary provenance and tectonic background of the northern Ordos Basin. The results show that the sediments of the western sections were mainly derived from the Yinshan orogenic belt and Alxa block, and that those in the eastern sections only came from the Yinshan orogenic belt. The trace element ratios in detrital zircons from the Late Palaeozoic sandstones indicate that the source areas were mainly subduction-related continental arcs, closely related to the continued subduction of the Palaeo-Asian Ocean in the Late Palaeozoic. Since the main Late Palaeozoic magmatic periods vary on the east and west sides of the northern margin of the Ordos Basin, two main collisions related to Palaeo-Asian Ocean closure are recorded. The collision on the west side occurred significantly earlier than that in the east. This study implies that the Palaeo-Asian Ocean began to subduct beneath the NCC in the Carboniferous and gradually closed from west to east thereafter.

In the Singhbhum Craton of the Indian shield, the Remal granite-gneiss preserves felsic magmatic fabrics onto which a low-temperature segregation layering has been superposed. Planar, sub-horizontal to gently dipping layers (Sign1) comprise K-feldspar megacrysts, plagioclase and quartz, with the base of each layer defined by segregations of biotite. Sign2 consists of trough cross-bedded layers composed of K-feldspar phenocrysts, plagioclase and quartz with biotite schlieren defining the base of each layer. Microstructural features such as concentrically arranged mineral inclusions in K-feldspar phenocrysts and graphic intergrowth textures testify to the magmatic origin of these fabrics, with insignificant subsequent metamorphic reconstitution. The tectonic fabric S1 has developed sub-parallel to localized greenschist-facies mylonite bands, and is defined by weakly aligned flakes of biotite. Crystallographic preferred orientations away from the mylonitized domains show a strong alignment of K-feldspar, quartz and biotite parallel to the magmatic fabric due to efficient segregation during magmatic flow. Quartz crystallographic preferred orientations within the mylonitized domains show a strong preferred orientation and dextral asymmetry. Temperature constraints from synkinematic chlorites along with estimates of deformation temperature from quartz crystallographic preferred orientations indicate that mylonitization occurred at the lower limits of quartz crystal plasticity. The results of combined thermodynamic and multiphysics modelling studies show that felsic magmas can undergo significant convective motion for a wide range of crystallinities and water contents before solidification. Additionally, segregation layering resembling a gneissosity can develop at low temperatures owing to localized mylonitization and concomitant dissolution–precipitation of biotite.

We conducted a U–Pb–Hf isotope study on zircon crystals from ignimbrites of the Changhsingian to Olenekian (253–248 Ma) Los Menucos Basin in the North Patagonian Massif (NPM), Argentina, in order to evaluate the age and petrogenesis of the magmas. Additionally, a compilation of whole-rock geochemistry and U–Pb–Hf in zircon isotope data for the Permian to Middle Triassic rocks of the NPM, for comparison with our data, was made to assess whether Patagonia would have been an exotic terrane accreted to SW Gondwana during the late Palaeozoic. We interpret the available U–Pb–Hf data to suggest that northern Patagonia experienced eastward arc expansion from the early Permian, about 273 Ma ago. This ∼820 km arc expansion event involved crustal shortening and magmatism with high-silica adakitic affinity, resulting in Hf-isotopic pull-down. At 253 Ma, slab steepening became associated with the coeval emplacement of ignimbrites of the Los Menucos Basin, which involved post-orogenic to intraplate magmatism. During the Middle Triassic, a slab break-off triggered uplift and basaltic underplating, promoting the emplacement of dike swarms with C-type adakitic signature at 246–244 Ma. The Hf isotope data for SW Gondwana for the same period indicate distinct trends that are explained here by differential slab roll-back since the Guadalupian, in a slab-tearing setting. Therefore, Permian to Middle Triassic magmatism is interpreted as having been associated with an eastward-directed proto-Pacific subduction system, which ultimately supports an autochthonous origin for Patagonia.

Detailed structural investigations across the Achankovil terrane boundary shear zone (AKSZ) system show distinct differences in the geometry of superposed fold structures between the Trivandrum (TB) and Southern Madurai (SMB) blocks separated by the AKSZ. The metamorphic temperatures and pressures estimated from the SMB and TB have a similar range, that is, 600–880°C and 5–8 kbar. The similar clockwise P-T paths retrieved by phase equilibrium modelling from both the blocks represent the last deformation and metamorphism shared between them during their accretion along the AKSZ. The distinct evolutionary history of the SMB and TB prior to their amalgamation is supported by the contrasting structural, metamorphic and chronological patterns, particularly the lack of prominent middle Neoproterozoic ages and the presence of Palaeoproterozoic ages in the TB, and vice versa in the SMB. The prominent 600–500 Ma monazite ages in the TB, SMB and AKSZ attest to their timing of accretion along the AKSZ. The study corroborates the S-directed subduction model proposed for terrane accretion along the AKSZ, and provides further insight into the subduction–accretion–collision tectonics associated with the late Neoproterozoic – Cambrian evolutionary history of this region.

The tectonothermal histories and geochronology of the SMB and TB are compared to the once adjoined crustal domains of Madagascar within Gondwanaland. It is suggested that the TB is equivalent to the Androyan and Anosyan domains of southern Madagascar, the SMB is equivalent to the Antananarivo and Itremo–Ikalamavony blocks of the central Madagascar, and the sinistral AKSZ is contiguous with the Ranotsara Shear Zone in southern Madagascar.

Field investigations into the sedimentary evolution of the Mulargia–Escalaplano late- to post-Variscan basin in Sardinia have been carried out. This basin is the only Variscan molassic basin in the Southern Variscan Realm where relationships between all the Late Pennsylvanian to Middle Triassic sedimentary cycles have been detailed for the first time. Here the interplay between tectonics and sedimentation can also be inferred. The investigations evidence that the sedimentary response recorded in the basin matches the progress and final collapse of the southern Variscan chain branch. The basin preserves continental deposits of decreasing energy organized in cycles, separated by weak unconformities related to abrupt volcano-tectonic climaxes and influenced by progressively drying climates. The sedimentary facies and stratigraphy suggest a progressive widening of the basin and decreasing depositional energy to the southeast as well as the migration of the basin depocentre in the same direction. Stratigraphic, environmental and evolutionary correlations between the presently separated NW and SE parts of the basin are reconstructed by considering the role of tectonics. As a result, a detailed lithostratigraphic framework for the Mulargia–Escalaplano basin is proposed. The findings in the sedimentary evolution match the evolution of similar coeval basins of France and Germany, following a pull-apart to extensional model. Based on comparisons with the other Sardinian and SW Europe coeval basins, relationships between ‘limnic’ and red bed successions in SW Europe during the collapse of the Variscan chain are discussed. The renaming of the European tectonomagmatic ‘Mid-Permian Episode’ to the ‘Lower Permian Episode’ is proposed.

Palaeopascichnida is a problematic group of extinct organisms that is globally distributed in Ediacaran sequences of Avalonia, Baltica, Siberia, South China and Australia. The fossils related to Palaeopascichnida consist of serially or cluster-like arranged, millimetre- to centimetre-scale globular or allantoid chambers, which are characterized by substantial differences in preservation, leading to no consistent diagnosis for these organisms. Here we integrate morphometric variation, stratigraphic distribution and habitat settings of more than 1200 specimens from all known fossil localities. The results of the morphological analysis demonstrate variation in chamber shape and size, and allow us to recognize six valid species within the group. Statistical analysis of the specimen distribution with respect to sedimentary environments indicates a significant difference in palaeoecological settings between species, making a significant contribution to the evolution and systematic palaeontology of these problematic organisms and perspective on their use in Neoproterozoic biostratigraphy. Our revision and systematic study sheds new light on one of the least studied groups of the late Ediacaran biota.

Cheloniellida is a rare but significant clade of artiopod euarthropods, and a major component of Vicissicaudata alongside aglaspidids. The first cheloniellid to be described was Triopus draboviensis, based on a single specimen from the Upper Ordovician Letná Formation in Bohemia. Because this specimen was believed to be missing its anterior, Triopus draboviensis has been the subject of considerable speculation regarding its affinities, with several authors combining various isolated cephalic shields with its trunk. Notwithstanding, Triopus draboviensis has eluded detailed restudy in the 150 years since its first description. Here we provide a redescription of the holotype and hitherto only known specimen of Triopus draboviensis, and describe a recently discovered second partial specimen as Triopus sp. It is shown that the holotype of Triopus draboviensis preserves the complete cephalic shield, while the trunk reveals a likely articulating device, which may represent an apomorphy for Cheloniellida. Further information on the trunk morphology provided by the new specimen has allowed a complete reconstruction of the dorsal exoskeleton of Triopus. An assessment is made of fossils previously suggested to represent the cephalic shield of Triopus draboviensis; while it is shown that none of those fossils can be attributed to Triopus, it is reaffirmed that Drabovaspis complexa likely does belong to Cheloniellida, making Drabovaspis the third cheloniellid genus known from the Letná Formation, after Triopus and Duslia. A revised diagnosis for Cheloniellida is provided, and the recently described problematic euarthropod Parioscorpio venator from the Silurian of Wisconsin is firmly rejected from this clade.

The subduction model of the Neo-Tethys during the Early Cretaceous has always been a controversial topic, and the scarcity of Early Cretaceous magmatic rocks in the southern part of the Gangdese batholith is the main cause of this debate. To address this issue, this article presents new zircon U–Pb chronology, zircon Hf isotope, whole-rock geochemistry and Sr–Nd isotope data for the Early Cretaceous quartz diorite dykes with adakite affinity in Liuqiong, Gongga. Zircon U–Pb dating of three samples yielded ages of c. 141–137 Ma, indicating that the Liuqiong quartz diorite was emplaced in the Early Cretaceous. The whole-rock geochemical analysis shows that the Liuqiong quartz diorite is enriched in large-ion lithophile elements (LILEs) and light rare-earth elements (LREEs) and is depleted in high-field-strength elements (HFSEs), which are related to slab subduction. Additionally, the Liuqiong quartz diorite has high SiO2, Al2O3 and Sr contents, high Sr/Y ratios and low heavy rare-earth element (HREE) and Y contents, which are compatible with typical adakite signatures. The initial 87Sr/86Sr values of the Liuqiong adakite range from 0.705617 to 0.705853, and the whole-rock ϵNd(t) values vary between +5.78 and +6.24. The zircon ϵHf(t) values vary from +11.5 to +16.4. Our results show that the Liuqiong adakite magma was derived from partial melting of the Neo-Tethyan oceanic plate (mid-ocean ridge basalt (MORB) + sediment + fluid), with some degree of subsequent peridotite interaction within the overlying mantle wedge. Combining regional data, we favour the interpretation that the Neo-Tethyan oceanic crust was subducted at a low angle beneath the Gangdese during the Early Cretaceous.

Ophiolitic mélanges in Anatolia represent Mesozoic subduction-accretion complexes, which are unusually poor in land-derived coarse-clastic rocks. A segment of the ophiolitic mélange in the Beynam region south of Ankara was studied. The ophiolitic mélange consists of three accretionary units (AUs), which are distinguished by lithology, structure, age and geochemistry. At the base there is a serpentinite mélange, which is overlain by a semi-intact Upper Jurassic ophiolite with boninite geochemistry. The topmost AU consists of ocean-island-like alkali basalts with seamount-derived Triassic shallow-marine limestones and Jurassic radiolarian cherts, which are stratigraphically overlain by Upper Cretaceous fore-arc turbidites. The base of the fore-arc sequence is palaeontologically and isotopically dated to the early to middle Campanian (c. 81 Ma). Detrital zircons from the fore-arc sequence indicate a Late Cretaceous (87–81 Ma) magmatic arc as a source. The formation of the subduction-accretion complex was a two-stage process. The first stage took place during the Late Jurassic – Early Cretaceous, when supra-subduction type oceanic crust was generated, and subduction accretion was intra-oceanic. In the second stage during the Late Cretaceous the subduction jumped inboard, creating an Andean-type convergent margin, and the Jurassic oceanic crust was incorporated in the subduction-accretion complex. The lack of land-derived sandstones in the ophiolitic mélange can be attributed to the intra-oceanic subduction and to the limestone deposition in the upper plate during the main phase of subduction accretion in the Late Jurassic – Early Cretaceous.

Marble in the supracrustal rocks of the Lewisian Complex, Tiree, includes chlorine-bearing amphiboles, chlorine-rich apatite, sulphur-rich scapolite, albite and phlogopite, all of which are regarded as evidence for evaporites in other metamorphosed sequences. Titanite yields U–Pb ages of ∼1.6 Ga, i.e. late Laxfordian, which excludes a younger imprint of sodium metasomatism. Traces of anhydrite, and isotopically heavy pyrite, also indicate deposition from seawater. Elsewhere in the Hebrides, tourmaline in Lewisian Complex marbles may represent seafloor exhalative deposits. Combined, the evidence suggests Lewisian Complex supracrustal marbles formed in an evaporative environment, like other Palaeoproterozoic successions across the North Atlantic region.

Understanding the distribution and geometry of faults and fractures is critical for predicting both subsurface permeability architecture and the integrity of geological fluid barriers, particularly in rocks with low primary porosity and permeability. While fracture patterns in relatively competent, weathering-resistant (therefore often well-exposed) rocks are generally well studied in outcrop, the role of mechanically weak layers in defining fracture patterns is frequently overlooked or under-represented. Here we show that rock composition, specifically clay and silicate minerals versus carbonate content, exerts a strong control on fault and fracture propagation and bed-containment within a mechanically layered, Cretaceous carbonate sequence at Canyon Lake Gorge, Texas. We find that relatively incompetent, clay-rich layers limit fault and fracture propagation, and cause bed-containment of fractures in more competent beds. In our results, no clear relationships exist between mechanical layer thickness and fracture abundance. These results are important for understanding the relative importance of composition versus bed thickness on fracture abundance in the subsurface, and for predicting fracture-controlled fluid flow pathways, seals and fracture connectivity across beds with variable compositions, thicknesses and competences.

Geothermal systems in terrains affected by polyphase deformation have reservoirs with a 3D geometry that is always difficult to predict. In this paper we describe a fossil exhumed geothermal system exposed in eastern Elba Island that developed in a polyphase folded and faulted setting, which can help us to understand how geothermal fluids circulate in geological bodies with inherited structures. Geothermal circulation at Elba allowed the deposition of Fe-ore deposits (haematite/magnetite and pyrite) and altered rock volumes, which represent tracers of the palaeo-fluid flow. Normal and oblique-slip faults dissected a polyphase folded metasiliciclastic succession and produced a secondary permeability in the range of 5 × 10−13 to 5 × 10−16 m2. From the permeable fault zones acting as feeder conduits, geothermal fluids permeated the hydraulically connected metasiliciclastic rock bodies previously deformed by two generations of folds. Geothermal fluids followed the already defined geometry, thus giving rise to apparent folded mineralized levels. Fluid migration into the metasiliciclastic rocks was possible due to their chemical aggression, which favoured the dissolution and reprecipitation of quartz, and Fe-oxide and sulphide deposition. Renewed fluids maintained their chemical properties (pH value and temperature, mostly). This conclusion provides inputs for reconstructing geothermal conceptual models and evaluating the geothermal potentiality of exploitable areas developing in similar geological settings.

Calcite veins are common in organic-rich mudrocks, but their genesis and ability to transmit fluids are debated. A combined microstructural and isotopic investigation of an array of calcite veins recovered in core from the Marcellus Formation reveals that the veins grew via a combination of continuous fibrous growth and punctuated fracture-opening increments. Continuous opening is the result of pressure-solution creep and involves no mechanical fracturing, but rather the growth of a pressure fringe around a pre-existing, sealed fracture. In contrast, incremental opening is accomplished by overpressured, mineral-saturated fluid, which repeatedly ruptures the rock at the cement / host-rock interface. Punctuated growth increments occurred repeatedly throughout an otherwise protracted, continuous growth history, indicating that the present structures preserve hybrid deformation conditions between brittle, fluid-assisted cracking and plastic strain. Stable isotopic signatures match those of a regional opening-mode fracture set that formed in response to catagenetic fluid overpressures amid a tectonically imposed (Alleghanian) stress field. It is concluded that calcite veins form as opening-mode hydraulic fractures and are susceptible to increments of brittle reactivation, even while inelastic growth processes widen and fill the veins with fibrous cement.

Faulted and fractured systems form a critical component of fluid flow, especially within low-permeable reservoirs. Therefore, developing suitable methodologies for acquiring structural data and simulating flow through fractured media is vital to improve efficiency and reduce uncertainties in modelling the subsurface. Outcrop analogues provide excellent areas for the analysis and characterization of fractures within the reservoir rocks where subsurface data are limited. Variation in fracture arrangement, distribution and connectivity can be obtained from 2D fractured cliff sections and pavements. These sections can then be used for efficient discretization and homogenization techniques to obtain reliable predictions on permeability distributions in the geothermal reservoirs. Fracture network anisotropy in the Malm reservoir unit is assessed using detailed structural analysis and numerical homogenization of outcrop analogues from an open pit quarry within the Franconian Basin, Germany. Several events are recorded in the fracture networks from the Late Jurassic the Alpine Orogeny and are observed to be influenced by the Kulmbach Fault nearby with a reverse throw of 800 m. The fractured outcrops are digitized for fluid flow simulations and homogenization to determine the permeability tensors of the networks. The tensors show differences in fluid transport direction where fracture permeability is controlled by orientation compared to a constant value. As a result, it is observed that the orientation of the tensor is influenced by the Kulmbach Fault, and therefore faults within the reservoirs at depth should be considered as important controls on the fracture flow of the geothermal system.

Widely distributed Mid-Neoproterozoic mafic rocks of the Qilian – Qaidam – East Kunlun region record the tectonic evolution of the northeastern Tibetan Plateau. This study presents whole-rock geochemistry, zircon U–Pb geochronology and Hf isotopes for the Xialanuoer gabbros of the central South Qilian Belt (SQB). Zircon laser ablation – inductively coupled plasma – mass spectrometry (LA-ICP-MS) U–Pb dating indicates that the gabbros were emplaced at ca. 738 Ma, indicating they are contemporaneous with mafic magmatism elsewhere in the northeastern Tibetan Plateau. The gabbros have low SiO2, Cr and Ni contents and Mg# values, are relatively enriched in light rare-earth elements (LREEs) and depleted in high-field-strength elements (HFSEs; e.g. Nb and Ta), have no positive Zr or Hf anomalies and have relatively high Nb/Ta but low Nb/La ratios. These data indicate that the Xialanuoer gabbros formed from calc-alkaline basaltic magmas that were originally generated by the partial melting of an enriched mantle of type-I (EMI-type) enriched region of the lithospheric mantle, underwent little to no crustal contamination prior to their emplacement, and have within-plate basalt geochemical affinities. Combining these data with the presence of widespread contemporaneous continental rift-related magmatism and sedimentation in the North Qilian, Central Qilian, South Qilian, Quanji, North Qaidam and East Kunlun regions suggests that the northeastern Tibetan Plateau underwent Mid-Neoproterozoic continental rifting, which also affected other Rodinian blocks (e.g. Tarim, South China, Australia, North America and Southern Africa).

New palynological results from the Gafo Formation (southern sector of the Pulo do Lobo Domain, South Portuguese Zone) are integrated with recently studied sections and drillholes from the Portuguese and Spanish sectors. A total of 44 samples were studied, 27 of which were positive for palynomorph taxonomy. This research revealed well preserved palynological assemblages, including 73 spore species allocated to 28 genera, four acritarch genera, three prasinophyte algae genera plus common chitinozoan remains. Some additional forms were retained under open nomenclature. From this, the first complete age determination for the Gafo Formation in Portugal and Spain was achieved, indicating a middle Famennian (Grandispora gracilis–Grandispora famenensis, GF Biozone) to a late Famennian (Grandispora echinata, VH Biozone) age. A greywacke sample from the same Gafo Formation was dated by U–Pb zircon geochronology and a maximum depositional age of 369 ± 2.5 Ma was obtained. A correlation between these palynological and U–Pb zircon data and the palynoflora assemblages of previous authors was made, as well as with the ages of felsic volcanic rocks found intercalated in the Gafo Formation, confirming the complex stratigraphy of Pulo do Lobo Domain. The results are consistent with stratigraphic mapping and structural interpretations, allowing a middle–late Famennian age (GF and VH Biozone) to be assigned to the Gafo Formation sedimentary rocks. This work has also contributed to a reinterpretation of Gafo Formation depositional facies correlatives (e.g. the Santa Barbara Group in Spain) as being the same lithological unit.

In the Palghat–Cauvery Shear/Suture Zone of the Southern Granulite Terrane, the Sittampundi Layered Complex occurs as a mappable unit. The Sittampundi Layered Complex consists of mafic, ultramafic and anorthositic rocks with chromitite layers and has been interpreted as an arc/ophiolite complex that formed in an Archaean suprasubduction zone arc setting. In the Sittampundi Layered Complex, reddish-black metabasites occur as layers or boudins with a rim of amphibolite within anorthosite. The peak metamorphic assemblage of the metabasites is garnet + clinopyroxene + quartz + rutile ± plagioclase ± orthopyroxene, and symplectite consisting of amphibole and plagioclase formed around the garnet during retrograde metamorphism. The protolith of metabasite may have intruded in a suprasubduction zone arc setting during the late Neoarchaean (c. 2540–2520 Ma), and then underwent high-pressure granulite-facies peak metamorphism (900–800 °C and 11–14 kbar) in the early Palaeoproterozoic (c. 2460–2440 Ma), followed by amphibolite-facies metamorphism (550–480 °C and 5.5–4.5 kbar) in the middle Palaeoproterozoic (c. 1900–1850 Ma). The results obtained in this study, together with previous studies, indicate that: (1) the high-pressure granulite-facies metamorphism in the study area indicates that subduction occurred during the Archaean–Proterozoic boundary with a higher apparent average geothermal gradient (∼20–16 °C km−1) than the modern-day Earth, and (2) the apparent average geothermal gradient of the subduction zone was ∼29–14 °C km−1 during the Archaean–Proterozoic boundary, which was still too high to enter the realm of eclogite-facies metamorphism.