The Niagara Escarpment is a fractured Palaeozoic sedimentary cuesta, subject to year-round weathering in a temperate climate. We examined the temperature of the rock surface and fractures at three in situ sites with varying aspect and lithology, as well as the surface and interior of three control blocks maintained in outdoor conditions between December 2020 and March 2021. The objectives were to examine the interplay between freeze–thaw and thermal weathering in the winter months and to identify potential factors influencing these processes. Both diurnal-scale and prolonged freeze–thaw cycles differing in spatial and temporal extent were identified, coincident with periods of high moisture. We frequently observed rapid temperature changes (>1 °C min−1) at sites with strong insolation, which implies that the temperature regime is suitable for thermal shock and fatigue to occur. Site-specific factors, such as the aspect of the escarpment face and lithology, impact the mechanism and extent of weathering. Southeast-facing sites with high insolation are dominated by diurnal-scale freeze–thaw; west- and east-facing sites with lower insolation experience a more prominent prolonged freeze–thaw cycle. Across all sites there is a gradient between surface and fracture temperature that follows diurnal trends in air temperature and insolation. Variability in the surface-fracture gradient may enhance weathering processes by shifting the orientation and magnitude of stress, and by changing the spatial distribution of freezing and thawing. Our research indicates that site-specific factors and pre-existing fractures moderate the influence of air temperature and insolation on thermal gradients, and ultimately the weathering regime.

Thrombolite reefs with archaeocyaths are common in the subtidal limestones of the lower Cambrian in the western Anti-Atlas of Morocco. The Igoudine Formation of the Tata Group recorded the first replacement of the microbial consortium (stromatolite-dominated) by thrombolite reefs with archaeocyaths and shelly metazoans. In order to better understand the role of the microbial community in the formation of thrombolite reefs with archaeocyaths across this critical transition, the macro-, micro- and ultra-fabric of thrombolites have been studied in detail. Three major components are identified within the first thrombolytic reef: archaeocyaths, calcimicrobes and micritic matrix. The studied thrombolites are typically dominated by the calcimicrobe Renalcis with subordinate Epiphyton and Girvanella. Scanning electron microscopy of the dark micrite of the Renalcis chambers showed amorphous translucent sheet-like structures interpreted as extracellular polymeric substances, closely associated with organominerals including nanoglobules and polyhedrons. Exceptionally well-preserved Renalcis chambers contain bacterial fossils similar to those described in modern microbialites, including microspherical coccoid fossils and filamentous bacteria that are either spaced or in close associations forming colonies. These organomineralization-related features suggest a bacterial origin for the Renalcis calcimicrobe. The matrices between the Renalcis chambers consist predominantly of clotted peloidal micrite. Mineralization of Renalcis microframes may involve two major biomineralization processes: (1) replacement of organic matter by organominerals resulting from anaerobic degradation of extracellular polymeric substances and bacterial sheaths and (2) encrustation of bacterial sheaths and extracellular polymeric substances due to increasing alkalinity of the microenvironment. These mechanisms played a crucial role in the early diagenetic cementation and preservation of the studied reefs.

Jurassic carbonate strata in the eastern Paris Basin exhibit several generations of faults, tension gashes and stylolites. Although their relative chronology can sometimes be determined according to cross-cutting relationships, the duration of major deformation phases and their influence on fluid flow and carbonate cementation are still uncertain. This contribution aims to clarify the timing of brittle deformation and associated calcite cementation. Tension gashes filled by calcite in Jurassic carbonates were sampled in outcrops and boreholes and dated through U–Pb geochronology. Almost all the sampled fractures were cemented during the Cenozoic period. Continuous deformation spread from c. 50 to 30 Ma. Tension gashes oriented N10° to N20° dated at 48–43 Ma show the main Pyrenean contractional stage. A second set of calcites were dated at c. 35–33 Ma and document a Late Eocene – Oligocene extension. A transition from the compressional to the extensional regime is expressed by tension gashes dated between 43 and 35 Ma. Finally, tension gashes oriented N150° to N175°, dated between 32 and 18 Ma, may result from the propagation of the horizontal stress generated by the Alpine orogen or by late Pyrenean deformation. Clumped isotope thermometry on five samples revealed both low crystallization temperatures (from 27 to 53 °C) and the meteoric origin of calcite-precipitating fluids. Our research therefore documents a continuous fracturing from Ypresian to Rupelian times, and less expressed brittle deformation during the Miocene period.

To study the ore mineralization at the outcrop scale we merge an advection–diffusion simulation with the geochemical software iphreeqc to model the mixing of two realistic fluids. We simulate the infiltration of a metal-rich fluid into a rock that is saturated with pore fluid. We test the feedback effects with a number of scenarios based on an outcrop-scale 5 × 5 m model consisting of two high-permeable vertical faults within a low-permeable host rock that lead into a permeable layer. The hot metal-rich fluid enters the model through the faults from below. We solve the advection–diffusion equation for 12 chemical species and temperature, and use iphreeqc to determine the resulting properties of local fluid domains as well as related saturation indices for minerals. The faults in the model act as pathways for the metal-rich fluid, with the infiltrating fluid displacing the pore fluid. Mixing in the model takes place as a function of advection along permeable faults coupled with diffusion of chemical species at the interface between two fluids, while heat diffusion is fast enough (103 times faster) to equilibrate temperature. Simulations show a high saturation index of mixing-derived minerals such as barite at the interface between the two fluids as a result of fluid mixing. Fast fluid pathways (i.e. faults) show travelling waves of high saturation indices of barite, while low-permeability zones such as fault walls and areas below less permeable layers experience stationary waves of high saturation indices. Our results show that, depending on the dominant transport process (advection or diffusion), mineralization will localize next to permeability contrasts in zones where local diffusion dominates.

Fractures that form when fluid pressure ruptures the rock are referred to as fluid-driven fractures or hydrofractures. These include most dykes, inclined sheets and sills, but also many mineral veins and joints, as well as human-made hydraulic fractures. While considerable field and theoretical work has focused on the geometry and arrest of hydrofractures, how they select their propagation paths, particularly in layered and faulted rocks, has received less attention. Here I propose that of all the possible paths that a given hydrofracture may follow, it selects the path of least (minimum) action as determined by Hamilton’s principle. This means that the selected path is that along which the energy transformed (released) multiplied by the time taken for the propagation is a minimum. Hydrofractures advance their tips/fronts in steps, with a time lag between the fracture front and the fluid front. In the present framework, each step is then controlled by Hamilton’s principle. The results suggest that when the hosting rock body is regarded as homogeneous, isotropic and non-fractured, hydrofracture paths are everywhere perpendicular to the trajectories of the minimum compressive (maximum tensile) principal stress σ3 and follow the trajectories of the maximum principal compressive stress σ1. When applied to layered and faulted rock body, the results indicate that hydrofracture paths may follow existing faults for a while, depending primarily on (1) the dip of the fault (steep faults are the most likely to be used by vertically propagating hydrofractures), and (2) the tensile strength across the fault compared with the tensile strength of the host rock along a path following the direction of σ 1. The results suggest that hydrofractures may use faults as parts of their paths primarily if the fault is steeply dipping and with close to zero tensile strength.

During the early Bajocian, a conspicuous coal-bearing siliciclastic succession was deposited in the northern Tabas Bock, which is important for understanding the regional geodynamics of the Central-East Iranian Microcontinent (CEIM) as well as for the Jurassic coal genesis in this part of Laurasia. Sedimentary facies analysis in a well-exposed section of the lower Bajocian Hojedk Formation (Kalshaneh area, northern Tabas Block) led to the recognition of ten characteristic sedimentary facies and three facies associations, representing channels with point bars and floodplains of a Bajocian meandering river system. Modal analysis indicates that the mature quartz arenites and quartzo-lithic sandstones of the Hojedk Formation originated from the erosion and recycling of older, supracrustal sedimentary rocks on the Yazd Block to the west. The coal petrography and maturity show an advanced maturation stage, whereas the great thickness of these continental strata points to a pronounced extension-related subsidence in the northern Tabas Block. The rapid rate of differential subsidence can be explained by accelerated normal block-faulting in the back-arc extensional basin of the CEIM, facing the Neotethys to the south. Compared to the thick Jurassic, the post-Jurassic strata are relatively thin and played a limited role in the thermal history of the coal in the northern Tabas Block. A relatively high geothermal gradient in the tectonically highly mobile area of the northern Tabas Block and/or heating by regionally widespread Palaeogene intrusions were most probably the key drivers of the thermal maturation of the Middle Jurassic coals.

Horseshoe crabs as a group are renowned for their morphological conservatism punctuated by marked shifts in morphology associated with the occupation of non-marine environments and have been suggested to exhibit a consistent developmental trajectory throughout their evolutionary history. Here, we report a new species of horseshoe crab from the Ordovician (Late Sandbian) of Kingston, Ontario, Canada, from juvenile and adult material. This new species provides critical insight into the ontogeny and morphology of the earliest horseshoe crabs, indicating that at least some Palaeozoic forms had freely articulating tergites anterior to the fused thoracetron and an opisthosoma comprising 13 segments.

Rounded to angular granoblastic textured mafic xenoliths, ranging from ∼1 to 6 cm in dimension, occur together with mantle peridotite xenoliths in a small Neogene/Quaternary alkali basalt cone in northeastern Iran. These crustal xenoliths consist of plagioclase feldspar, clinopyroxene, orthopyroxene and minor olivine, spinel, titanomagnetite and apatite. Their bulk compositions are similar to tholeiitic basalts and they are interpreted as meta-gabbroic rocks derived from mid- to lower crustal depths of 18 to 30 km. Rb–Sr dating suggests an age of c. 457 ± 95 Ma for these crustal xenoliths, and their geochemistry shows some similarities to Ordovician gabbros that crop out ∼20 km to the west. The data suggest that the gabbroic proto-lithologies of the xenoliths formed by intrusion of mafic magmas into the mid- to lower crust, possibly during extension and magmatism related to the opening of the Hercynian Palaeotethys ocean that separated central and eastern Iran from the Eurasian plate during the Late Palaeozoic.

Turbidity currents commonly bypass sediment in submarine channels on the continental slope, and deposit sediment lobes farther down-dip on the flat and unconfined abyssal plain. Seafloor and outcrop data have shown that the transition from bypass to deposition usually occurs over complex zones referred to as channel–lobe transition zones (CLTZs). Recognition of these zones in cores and outcrop remains challenging due to a lack of characteristic sedimentary facies and structures. This paper focuses on Unit E of the Permian Fort Brown Formation in the Karoo Basin, South Africa, in the Slagtersfontein outcrop complex, which has previously been interpreted as a CLTZ. This study integrates thin-section micrographs, sedimentary facies, bed-set and stratigraphic architecture, and palaeoflow directions to achieve a multiscale analysis of CLTZ features. A novel process-based facies scheme is developed to evaluate deposits in terms of the depositional or erosional tendencies of the flows that formed them. This scheme allows bypass to be distinguished from depositional zones by the spatial distribution of certain sediment facies. Areas of net sediment bypass were predominantly marked by erosive sediment facies and a larger variability in palaeoflow direction while depositional areas showed a lower variability in palaeoflow directions. Metre-scale structures in the bypass-dominated area reveal seafloor erosion and scour formation. Field relations suggest the presence of a ∼500 m long mega-scour in the CLTZ. The characteristic structures documented here are applicable for identifying CLTZs in sparse datasets such as outcrops with limited palaeogeographical context and sediment cores obtained from subsurface systems.

The late Palaeozoic continental-arc magmatic rocks in the Gongzhuling area are located in the Liaoyuan Accretionary Belt. Here we present new zircon U–Pb ages, whole-rock major- and trace-element compositions, Li and zircon Hf isotopic compositions and oxygen fugacity of these rocks with an aim to constrain the lithium isotopic composition of the source region and origin of the magmas. These rocks were formed during 269–258 Ma in middle–late Permian time. The dioritic rocks were formed through mixing processes, with the mafic melts originating from a metasomatized mantle wedge and the felsic melts from the lower crust of a Neoproterozoic arc. The mantle wedge has been metasomatized by Li-rich fluids derived from subducted oceanic crust, as indicated by the δ7Li values of +0.4 ‰ to +3.5 ‰ and positive ϵHf(t) values (+0.7 to +13.1). Redox-sensitive Ce in the zircons indicates the fO2 of the magmas to be low to intermediate (FMQ−2.2 to FMQ+2.6; FMQ is the fayalite–magnetite–quartz redox buffer), precluding large-scale porphyry Cu–Mo mineralization. The middle–late Permian magmatic rocks represent the terminal magmatic record of the subduction of the Palaeo-Asian oceanic crust, meaning that the final closure of the Palaeo-Asian Ocean in the eastern Central Asian Orogenic Belt occurred at the end of the Permian Period. Recent identification of Mesoproterozoic (c. 1400 Ma) granites suggests some Palaeoproterozoic crustal fragments still exist in the Liaoyuan Accretionary Belt, but only in a small amount; therefore, it is concluded that the crustal growth of the Liaoyuan Accretionary Belt occurred mainly during the Neoproterozoic period.

This study presents zircon U–Pb geochronology and Hf and O isotope data for granitic rocks in the Zijinshan ore field, southwestern Fujian Province, China. The intrusive rocks comprise monzogranite, granite and granodiorite. The magmatic zircon U–Pb ages from these granitic rocks can be divided into four episodes: episodes 1 (157.9–159.9 Ma) and 2 (141 Ma) in the Late Jurassic and episodes 3 (108.1–103.2 Ma) and 4 (97.5–99.7 Ma) in the Middle and Late Cretaceous, respectively. Patterns of rare-earth elements (REE) show enrichment in light REE and obvious negative Eu anomalies. These rocks are also enriched in Rb, Th, U, La, Ce, Nd and Hf, and depleted in Ba, Nb, Sr, P and Ti. The increasing La/Yb ratio and decreasing heavy REE content with decreasing age may imply an increasing contribution of mantle-derived materials from the Late Jurassic to Cretaceous. The zircon ϵ Hf(t) and δ18O values, ranging from −37.7 to −2.8 and 12.0 ‰ to 6.3 ‰, respectively, indicate that the lower crust is an important source of granitic rocks. There was a significant increase in ϵ Hf(t) values and a decrease in δ18O values in the younger magmatic episodes (3 and 4), which ranged from −11.4 to −0.6 and 10.7 to 6.3 ‰, respectively. This suggests an increasing contribution of mantle-derived magma to the crustal melts from the Late Jurassic to late Early Cretaceous in response to the changing regional tectonic setting from compression to extension and an increasing interaction between the crust and mantle.

We investigate the shape and strength of the magnetic fabrics (anisotropy of magnetic susceptibility (AMS) data) of various massive granitic plutons from different parts of India, using the eigenvalue method. The study aims to analyse eigenvalues and establish their relationship with various deformational attributes. It involves: (1) calculating eigenvectors and their corresponding eigenvalues from magnetic fabric datasets; (2) finding a link between the geometrical appearance of eigenvectors and the mechanistic issues involved with a specific deformation scenario; and (3) determining shape and strength parameters from the magnetic foliation data distribution.

The statistical analysis for the unimodal magnetic fabric dataset of orthorhombic symmetry class implies that the plane, consisting of intermediate (V2) and minimum (V3) eigenvectors with pole V1, accurately traces the instantaneous stretching axis (ISAmax) of a particular material flow system under a pure shear regime. Moreover, for the distributions of similar symmetry and modality, we infer that the rotational characteristics of eigenvectors with respect to a fixed coordinate cause a distinct shift of such planes (V2–V3) from the ISAmax of a steady-state flow system under simple shear, where a substantial amount of rotational strain is involved. However, our findings also suggest that variation in symmetry and modality of magnetic fabric data distribution of different studied granitoids can directly influence the relative disposition of V2–V3 with respect to the direction of ISAmax. We conclude that eigenvalue analysis of magnetic fabrics is a powerful approach, which can be utilized while studying the salient deformational aspects of any syntectonic massive granitic body.

Scapolite occurrences are widely observed in the metasedimentary rocks exposed around the Khetri Copper Belt and adjoining Nim ka Thana copper mineralized area in western India. Amoeboidal to well-developed and rounded/elliptical-shaped marialitic scapolite (Na-rich end-member) rich zones with variable Cl contents ranging from 1.0 wt % to 2.9 wt % have been identified in proximity to the ore-bearing hydrothermal fluid activity zones. Although scapolite is formed as a product of regional metamorphism in many places, in this study, we propose a strong possibility that scapolite was formed by hydrothermal ore-bearing fluid interaction with metasediments. The evidence of hydrothermal activity and Cl sourcing is attributed to (i) the absence of evaporite beds in the area and no Na-rich plagioclase as inclusions within the scapolite suggesting the formation of marialitic scapolite from sodic plagioclase in the metasediments with the interacting hydrothermal fluid; (ii) an epithermal to mesothermal hydrothermal fluid with moderate salinity responsible for the Cu mineralization that is ascribed to be the source of Cl for the formation of marialitic scapolite; (iii) diffusion of SO2 in the scapolite in close association with the sulfide mineral phase (chalcopyrite) supporting the involvement of ore-bearing fluid in the development of scapolite; (iv) the absence of zoned scapolite, the spatial distribution of scapolite in a particular lithology, the occasional incorporation of sulfur into marialitic scapolite and the texture/geometry in the scapolite suggesting a broad hydrothermal linkage instead of a pure metamorphic origin.

I welcome the comment from Dewey and Ryan on my paper (Searle, 2021), and the chance to clarify some conflicting Caledonian tectonic models and new ideas. As outlined in the paper, my review was based on the geology of Scotland, not Ireland, Newfoundland or the rest of the Caledonian orogeny beyond, mainly because of length constraints, and the fact that my personal field experience lies mainly in Scotland. The Great Glen fault and the Moine Supergroup are not exposed at all in Ireland, so critical parts, and half of the Caledonian orogeny are missing in Ireland. I argue that the Caledonian geology of Scotland can be explained by a Laurentia – Midland valley arc-microcontinent collision, progressing from NW-directed ophiolite obduction in the Late Cambrian to crustal thickening and regional metamorphism during the Ordovician and Silurian, across the Dalradian and Moine metamorphic rocks. Caledonian structures do continue north into Greenland, as summarized by Leslie et al. (2008), but the geological evolution of the Scottish Caledonides does not require any collision with Baltica. The Dewey and Ryan comment raises major tectonic questions, which this reply addresses.

Searle (2022) argues, based on data from Scotland alone, that the deformation and metamorphism of the Scottish Highlands is continuously diachronous from the early Ordovician (Grampian) in the Dalradian of the Grampian Highlands to the mid-Silurian (Scandian) in the Moine of the Northwest Highlands. Necessarily, he disputes substantial offset along the Great Glen Fault (GGF). He compares the Scottish Caledonides with the Himalayas. We take the different view, which also considers evidence from along-strike in Ireland and Newfoundland, that the Grampian and Scandian orogenies were discrete events. The Grampian Orogeny involved major deformation and metamorphism across the Grampian and Northwest Highlands to the Moine Thrust Zone, whereas the Scandian Orogeny affected only the Northwest Highlands. The early Ordovician Grampian Orogeny was generated by the collision and obduction of an oceanic arc-ophiolite with and onto the Laurentian continent. The mid-Silurian Scandian Orogeny was caused by continental collision between Laurentia and Baltica. Late Silurian – early Devonian sinistral slip along the GGF emplaced the the Moine of the Northwest Highlands from a distant location to the north in the Scandian collision zone, against the Dalradian, which had not been involved in the Scandian collision. We outline several (eight) arguments that we believe contradict evidence for progressive northwestward propagation of a Himalaya-style thrust wedge from the early Ordovician to the mid-Silurian.

The coastal region of Guangdong Province, China, is characterized by well-developed dyke swarms. The dykes with widths of decimetres to metres and lengths of tens of metres occur along straight and planar fractures cutting granites or volcanic rocks with Jurassic to Early Cretaceous ages. They show steep attitudes with strikes varying from NNW to NNE and from NE to SE, consistent with a stress regime transition from E–W to N–S extension. Major-element analysis on representative dyke samples reveals a composition range from basaltic to andesitic with a few dacitic outliers. Trace elements of most samples show notable Nb–Ta negative anomalies and Pb positive anomalies on the primitive mantle normalized spidergrams, characteristic of arc-related rocks. Among these, several dacite samples show notable fractionated medium rare earth elements and heavy rare earth elements and high Sr/Y (47–74) and La/Yb (15–21) ratios indicative of adakitic affinity. A few samples lack anomalies in Nb and Ta and have primitive trace-element ratios (e.g. Nb/La and Sr/Nd) or initial ϵNd values resembling ocean island basalt-like rocks. Rare earth element distribution patterns of all samples are right inclined and generally absent of Eu anomalies, which together with other trace-element indexes indicate an origin from heterogeneous mantle sources with depths below the stability field of plagioclase. Radiogenic isotopes, especially those of Nd and Sr (or Nd–Pb), define an array extended from the depleted to enriched mantle (EMII) provenance. LA-ICP-MS dating of zircon obtains ages between 110 and 70 Ma indicating emplacement of these dykes in latest Early to Late Cretaceous time, a time just before initiation of the South China Sea basin. In the context of regional geodynamics, it is proposed that these dykes were formed in a tectonic setting transiting from a Pacific-related back-arc to a passive continental margin pertaining to the development of the South China Sea basin.

During late Carboniferous time, the residual ocean basin gradually closed in West Junggar and only a small amount of seawater remained in the Hala’alat Mountain area, where discussions of provenance and tectonics are limited. In this study, LA-ICP-MS U–Pb dating and heavy mineral identification are conducted on the upper Carboniferous tuffaceous sandstones from the Hala’alat and Aladeyikesai formations in the Hala’alat Mountain area. The results reveal the low maturity of the clastic sediments, indicating proximal deposition. The Hala’alat Formation detrital zircons present a single peak (c. 330 Ma). Speculatively, the primary provenance is the Boshchekul–Chingiz Arc, and the secondary sources are the Darbut Tectono-Magmatic Belt and island arcs in the basin. The main peak and provenance of the Aladeyikesai Formation are similar to those of the Hala’alat Formation. Moreover, several age groups, namely, 370–344 Ma, 427–404 Ma and 478–476 Ma, potentially correspond to provenances of the Darbut Tectono-Magmatic Belt, the Boshchekul–Chingiz Arc and the Kujibai–Hongguleleng Ophiolitic Mélange Belt. In addition, the maximum depositional ages of the Hala’alat and Aladeyikesai formations calculated are 314.6 ± 1.54 Ma and 330.8 ± 0.61 Ma, respectively. Comprehensive analysis shows a relatively singular provenance of the Hala’alat Formation. While the provenance of the Aladeyikesai Formation shows inheritance, the provenance area extends northwards to the Kujibai–Hongguleleng Ophiolitic Mélange Belt. Furthermore, the closure of the Junggar Ocean during Carboniferous time caused the potential source region of the Hala’alat Mountain area to migrate northeastwards from Barleik Mountain to Xiemisitai Mountain. This study provides a basis for the analysis of regional geological evolution.

In volcano-tectonic regions, dyke propagation from shallow magmatic chambers is often controlled by the interaction of the local and regional stress fields. The variations of the stress fields result from a combination of factors including the regional tectonic stress, the geometry of pressurized magma chambers, the layering and the pre-existing discontinuities (e.g. fractures). In this contribution, we describe and apply a new multiparametric inversion technique based on geomechanics that can invert for both the far field stress attributes and the internal pressure of magma chambers or stocks, constrained by observed dyke or eruptive fissure orientations. This technique is based on the superposition principle and uses linear elastic models that can be solved using many types of numerical methods. For practical reasons, we chose a 3D boundary element method (BEM) for a heterogeneous elastic half-space, where magma chambers are modelled as pressurized cavities. To verify this approach, the BEM solution has been validated against the known 3D analytical solution of a pressurized cylindrical cavity. Then the effectiveness of this technique and its practical use is demonstrated through application to natural examples of dyke network development around two different volcanic systems, the Spanish Peaks (USA) and the Galapagos Islands (Ecuador). Results demonstrate that regional stress characteristics as well as the internal pressure of magma chambers can be estimated from observed radial and circumferential dyke patterns and some knowledge of magma chamber geometry.