Structural analysis, petrochronology and metamorphic petrology enable identification and bracketing of the timing of a newly mapped high-temperature ductile shear zone (Jagat Shear Zone (JSZ)) in the Himalayan metamorphic core in Central-Western Nepal. In situ U-Th-Pb monazite petrochronology constrains the timing of top-to-the-S/SW shearing between 28–27 Ma and 17 Ma. Burial and prograde metamorphisms in footwall rocks were linked to thrust-sense movement along the JSZ, while the hanging wall rocks were retrogressed and exhumed. The identification and age of the JSZ (as part of a regional system of shear zones: the High Himalayan Discontinuity (HHD)) coupled with the localization and timing of activity of the Main Central Thrust (MCT) (i) fills a gap in tracing the HHD along orogenic strike, (ii) supports the identification of the position and timing of the long-debated MCT and (iii) helps to place the boundaries of the Himalayan metamorphic core and its internal architecture. Thus, our study is a significant step towards a precise identification of the burial, assembly and exhumation mechanisms of the Himalayan metamorphic core.

Late Kungurian (Lower Permian) conodonts are described from the Kanthan Limestone, Perak, Peninsular Malaysia and for the first time from the Sibumasu Terrane of the Malay Peninsula. The co-occurrence of Gullodus duani, Gullodus hemicircularis, Gullodus sicilianus, Mesogondolella lamberti and Mesogondolella siciliensis represent the Mesogondolella lamberti International Conodont Zone and the broadly equivalent Mesogondolella siciliensis Regional Conodont Zone. A small fault-bounded basal Pennsylvanian (basal Bashkirian) conodont fauna including Gnathodus girtyi simplex and Declinognathodus inaequalis is also reported. The late Kungurian conodonts from the Kanthan Limestone were deposited in a relatively deep-water environment on the northern passive margin of the Sibumasu Terrane of the eastern Cimmerian Continent located at c. 35oS latitude. Biogeographically, the fauna represents the southern peri-Gondwana Cool Water Province which is consistent with its palaeogeographic location. A new scheme, utilizing characteristics of P1 elements, including position of the 1st denticle, location of 2nd and 3rd denticles, platform shape, platform cross-section, denticle shape in cross-section, and lateral denticle development is proposed for distinguishing between species of the hindeodid genera Gullodus, Hindeodus and Isarcicella. The late Kungurian fauna from the Kanthan Limestone represents the southern peri-Gondwana Cool Water Province supporting palaeogeographic reconstructions placing the Sibumasu Terrane in moderate southern palaeolatitudes in the Kungurian.

Fossils are more and more used in phylogenetic evolutionary studies either for clade calibration, or as terminals in a dataset including morphological characters. The strength of these methodological advances relies however on the quality and completeness of the fossil record. For crickets (Insecta, Orthoptera, Gryllidea), few ancient (pre-Cenozoic) well-preserved fossils are known, except for isolated wings often classified in purely fossil groups and a few fossils found in Cretaceous amber. Here, we present two remarkable fossils from mid-Cretaceous amber of France, that were imaged using X-ray synchrotron microtomography and exhibit an exquisite preservation allowing description with a precision similar to that of extant taxa. Palaeonemobius occidentalis Laurent and Desutter-Grandcolas, gen. nov., sp. nov. and Picogryllus carentonensis Josse and Desutter-Grandcolas, gen. nov., sp. nov. are the oldest representatives of the Nemobiinae and Podoscirtinae subfamilies of the Trigonidiidae and Oecanthidae families respectively. P. carentonensis Josse and Desutter-Grandcolas, gen. nov., sp. nov. is also the smallest adult male with a full stridulatory apparatus ever documented in crickets (body length 3.3 mm), and the first taxon of the cricket clade for which male genitalia can be partly described. We discuss the significance of Cretaceous fossils of crickets for future evolutionary studies of this clade.

We present U–Pb dates from peridotitic pyrope-rich garnet from four mantle xenoliths entrained in a kimberlite from Bultfontein, South Africa. Garnet dates magmatic emplacement due to the high mantle residence temperatures of the source material prior to eruption, which were most likely above the closure temperature for the pyrope U–Pb system. We determine a U–Pb date of 84.0 ± 8.1 Ma for the emplacement of the Bultfontein kimberlite from garnet in our four xenolith samples. The date reproduces previous dates obtained from other mineral-isotope systems (chiefly Rb–Sr in phlogopite). Garnet can be dated despite extremely low concentrations of U (median ∼0.05 μg/g), because concentrations of common Pb are often low or non-detectable. This means that sub-concordant garnets can be dated with moderate precision using very large laser-ablation spots (130 μm) measured by quadrupole inductively coupled plasma – mass spectrometry (LA-Q-ICP-MS). Our strategy demonstrates successful U–Pb dating of a U-poor mineral due to high initial ratios of U to common Pb in some grains, and the wide spread of isotopic compositions of grains on a concordia diagram. In addition, the analytical protocol is not complex and uses widely available analytical methods and strategies. This new methodology has some advantages and disadvantages for dating kimberlite emplacement versus established methods (U-based decay systems in perovskite and zircon, or Rb- or K-based systems in phlogopite). However, this method has unique promise for its potential application to detrital diamond prospecting and, more speculatively, to the dating of pyrope inclusions in diamond.

Loose limestone blocks of a newly recognized hydrocarbon-seep deposit from the lower Oligocene Jansen Creek Member of the Makah Formation were collected on a beach terrace close to the mouth of Bullman Creek in Washington State, USA. The limestone consists largely of authigenic carbonate phases, including 13C-depleted fibrous cement forming banded and botryoidal crystal aggregates with δ13C values as low as –23.5 ‰. Lipids extracted from the limestone yielded molecular fossils of anaerobic methanotrophic archaea (ANME), dominated by compounds of an ANME-2/DSS consortium with δ13C values as low as −106 ‰, indicating formation at an ancient methane seep. The fossil inventory of the seep deposit is remarkable, consisting almost solely of echinoid remains, whereas typical seep biota are absent. Varying preservation of the echinoid fossils indicates parautochthonous deposition, corroborated by evidence for high fluid flow at the ancient seep, possibly responsible for displacement of echinoids after death. Although a full taxonomic description of the echinoids cannot be given, almost all fossils were assigned to one taxon of irregular spatangoids, except for a single regular echinoid. Abundance and lifestyle of the irregular spatangoids in the Bullman Creek echinoid seep deposit resemble those of the fossil Tithonia oxfordiana from an upper Jurassic seep deposit in France and extant Sarsiaster griegii from modern seeps in the Gulf of Mexico. The Bullman Creek echinoid deposit probably represents a fossil analogue of the Gulf of Mexico Sarsiaster mass occurrence, indicating that the adaptation of spatangoid echinoids to chemosynthesis-based ecosystems ranges back at least to the earliest Oligocene.

Malformed trilobite specimens present important insight into understanding how this extinct arthropod group recovered from developmental or moulting malfunctions, pathologies, and injuries. Previously documented examples of malformed trilobite specimens are often considered in isolation, with few studies reporting on multiple malformations in the same species. Here we report malformed specimens of the ellipsocephaloid trilobite Estaingia bilobata from the Emu Bay Shale Konservat-Lagerstätte (Cambrian Series 2, Stage 4) on Kangaroo Island, South Australia. Ten malformed specimens exhibiting injuries, pathologies, and a range of teratologies are documented. Furthermore, five examples of mangled exoskeletons are presented, indicative of predation on E. bilobata. Considering the position of malformed and normal specimens of E. bilobata in bivariate space, we demonstrate that the majority of malformed specimens cluster among the larger individuals. Such specimens may exemplify larger forms successfully escaping predation attempts, but could equally represent individuals exhibiting old injuries that were made during earlier (smaller) growth stages that have healed through subsequent moulting events. The available evidence from the Emu Bay Shale suggests that this small, extremely abundant trilobite likely played an important role in the structure of the local ecosystem, occupying a low trophic level and being preyed upon by multiple durophagous arthropods. Furthermore, the scarcity of malformed E. bilobata specimens demonstrates how rarely injuries, developmental malfunctions, and pathological infestations occurred within the species.

The Palaeocene–Eocene Thermal Maximum (PETM) is one of the most important Cenozoic climatic events shaping modern biodiversity, yet reconstructions of its palaeobiomes remain controversial. Here we describe Gesomyrmex gallicus sp. nov., a new, extinct species of the ant genus Gesomyrmex Mayr, 1868, based on minor and major workers preserved in the early Eocene amber of Oise, France. Although there are only seven known extant species in the genus, all confined to warm and humid forests of SE Asia, the fossil record of Gesomyrmex indicates that the genus once radiated across Eurasia. The new species therefore provides direct evidence that this radiation likely co-occurred with the PETM. This observation constrains palaeoclimatic reconstructions of the Early Eocene by requiring the presence of an extensive, homogeneous and interconnected rainforest-like biome across palaeo-Eurasia, a scenario otherwise corroborated by pollen assemblages. The phases of regression of warm and humid forests and the widening of dry biomes probably occurred later during the Eocene, in between hyperthermals, and were probably less extensive than some computer models predicted.

The rock record from the late Early Jurassic in southern Africa encompasses the history of voluminous continental flood basalt outpourings associated with the magmatic events in the Karoo–Ferrar Large Igneous Province (LIP) in southern and eastern Gondwana. This multiphase magmatism produced one of Earth’s largest continental flood basalt successions volumetrically and is assumed to have been a main driving mechanism in late Early Jurassic global environmental perturbations, including mass extinctions and changes in climate. In southern Africa, these Lower Jurassic flood basalts are interbedded with fossiliferous sedimentary rocks, which in turn host the last signs of ‘Karoo life’ in the form of fossil plants, invertebrates and vertebrates, including the trackways of hopping mammals and the ultimate Karoo dinosaurs. The sedimentology and palaeontology of the interbeds archived depositional and biotic processes in running water as well as in and around shallow, up to ∼10 m deep freshwater lakes and ponds in the late Early Jurassic. This study explains how a complex freshwater palaeo-habitat prevailed – albeit temporarily – in this extremely stressful environment, which was unlike any modern volcanic system. The evidence collectively points to seasonally wet, warm temperate climatic conditions during the early phases of Karoo volcanism. Moreover, the evidence in the rocks also suggests that the dynamic volcanic conditions resulted in shifting habitats that likely facilitated the migration of the ultimate Karoo biota towards the north and west, away from the main Karoo land of fire, just before Gondwana started to disassemble. This refinement of the environmental dynamics in southern Gondwana presented herein lays the groundwork for future high-resolution volcanological, geochronological and chemostratigraphical studies aimed at the nuanced understanding of the global environmental effect of the Karoo–Ferrar LIP.

Hydrofractures, or hydraulic fractures, are fractures where a significantly elevated fluid pressure played a role in their formation. Natural hydrofractures are abundant in rocks and are often preserved as magmatic dykes or sills, and mineral-filled fractures or mineral veins. However, we focus on the formation and evolution of non-igneous hydrofractures. Here we review the basic theory of the role of fluid pressure in rock failure, showing that both Terzaghi’s and Biot’s theories can be reconciled if the appropriate boundary conditions are considered. We next discuss the propagation of hydrofractures after initial failure, where networks of hydrofractures may form or hydrofractures may ascend through the crust as mobile hydrofractures. As fractures can form as a result of both tectonic stresses and an elevated fluid pressure, we address the question of how to ascertain whether a fracture is a hydrofracture. We argue that extensional or dilational fractures that formed below c. 2–3 km depth are, under normal circumstances, hydrofractures, but at shallower depth they may, but must not be hydrofractures. Since veins and breccias are often the products of hydrofractures that are left in the geological record, we discuss these and critically assess which vein structures can, and which do not necessarily, indicate hydrofracturing. Hydrofracturing can suddenly and locally change the permeability in a rock by providing new fluid pathways. This can lead to highly dynamic self-organization of crustal-scale fluid flow.

The Twyfelskupje carbonatite complex, Southern Namibia, exhibits the typical igneous emplacement structures of carbonatites, including plugs, cone sheets and dyke stockworks. The excellent exposure allows for detailed studies of the high-level geochemical and structural evolution of the carbonatite, and the nature of the concomitant rare earth element mineralization. Radiogenic isotope analyses (Sr, Nd, Pb) reveal that, in common with many other carbonatites, the Twyfelskupje carbonatite complex appears to be predominantly derived from mixing between HIMU and EM1 mantle end-members. Following partial melting of these mantle sources, the geochemical and structural evolution of the Twyfelskupje carbonatite complex proceeded by a staged process involving separate magma pulses, a series of emplacement structures, sub-solidus crystallization, fractionation and low-temperature hydrothermal alteration. The dominant rare earth element minerals in the Twyfelskupje carbonatite complex are fluorcarbonates and monazite, and are characterized by variable Ca, high F and light rare earth elements in the order Ce>La>Nd. Comparison between the rare earth element concentrations of the whole rocks, dominant rare earth element minerals and carbonates suggests that ∼95 % of the total rare earth element abundance of the Twyfelskupje carbonatite complex is contained within fluorcarbonates and monazite. Overall, the early calcio-carbonatite plugs are rare earth element enriched (mean 4.47 wt % rare earth oxides) relative to the magnesio-carbonatite cone sheets (mean 2.51 wt % rare earth oxides).

The Arenig Fawr area of North Wales constitutes the type area for the British Lower to Middle Ordovician Arenig Series and is complemented by sections in the Carmarthen and Whitland areas of South Wales. We describe chitinozoan assemblages from both areas in order to aid correlation of the Arenig Series in its type region with the global Ordovician series and stages. Chitinozoans recorded from Arenig Fawr provide permissive rather than conclusive evidence but suggest that the Henllan Ash Member correlates with the upper Floian Stage Slice Fl3 or lower Dapingian Stage Slice Dp1. Better results were obtained from South Wales where six chitinozoan assemblages are distinguished, ranging in age from late Tremadocian to middle Darriwilian (early Llanvirn). Most species are known from South China, Gondwana and/or Baltica where there are controls on ranges. They show that much of the lower Arenig (Moridunian) succession in South Wales correlates with the upper Floian Stage (Fl3). Correlatives of the lower and middle Floian Stage (Fl1, Fl2), if present, must be represented by the Ogof Hên Formation and lowest Carmarthen Formation. Chitinozoan assemblages from the upper Arenig Series (Fennian Stage) are more readily correlated with Gondwanan biozones and indicate correlation of the Fennian Stage with the Dapingian and lower Darriwilian (Dw1) stages. The middle Arenig Whitlandian Stage is constrained in South Wales to an interval from the uppermost Floian Stage to the basal Dapingian Stage, resulting in an inferred increased rate of sediment accumulation.

Mudstones from the Argiles de Châtillon Formation exposed in the Boulonnais region of Northern France represent a proximal lateral equivalent of the organic-rich Kimmeridge Clay Formation. The Argiles de Châtillon Formation is composed of two subunits that straddle the Kimmeridgian–Tithonian boundary. Each subunit contains an organic-rich interval. The two conspicuous organic-rich intervals have been linked to either periods of high sea level or greenhouse warming. Here, we use palynology to further understand climate and environmental mechanisms that drove organic matter enrichment. We use bulk organic carbon isotope records (δ13Corg) to correlate the Boulonnais sections with those of the Kimmeridge Clay Formation. The palynological results suggest that the stratigraphically lower organic-rich interval (Kimmeridgian) was deposited under suboxic to anoxic stratified conditions. A large-scale climate shift from cooler/humid to warmer/arid conditions marked the Kimmeridgian–Tithonian boundary, influencing organic matter enrichment in the stratigraphically higher organic-rich interval (Tithonian). In contrast with the lower organic-rich interval, there are no indications of stratified conditions for the higher organic-rich interval. Within this thicker organic-rich interval, cyclic variations in amorphous organic matter distribution, total organic carbon and δ13Corg trends on a 2 m scale are observed. They co-occur with fluctuations of the palynological assemblages, indicative of more humid versus arid climate conditions, likely alternating on a ∼100 kyr eccentricity timescale. Our results show that under the most humid phases of these overall arid climate conditions, sulfurization of carbohydrates was the dominant control on organic matter preservation. This climate-controlled process that drives organic matter enrichment in the Tithonian can be recognized on a basin-wide scale.

Planktonic foraminifera were studied at Lydden Spout, near Folkestone (southeast England, UK), the reference section of the middle Cenomanian Event 1 (MCE 1) characterized by a prominent double-peak δ13C excursion of 1 ‰ identified in different ocean basins and considered a global event. Biostratigraphic and quantitative analysis of planktonic foraminifera are correlated to the δ13C perturbation, to the positive δ18O shifts identified within MCE 1 and to the occurrence of Boreal macrofossils (the bivalves Chlamys arlesiensis and Oxytoma seminudum, and the belemnite Praectinocamax primus). Variations in abundance and species richness of planktonic foraminifera and the inferred palaeoecological preferences of taxa permit the identification of distinct palaeoenvironmental settings across MCE 1. The stratigraphic interval corresponding to MCE 1 is characterized by the absence of oligotrophic rotaliporids, and by the evolutionary appearance of meso-eutrophic dicarinellids and of Muricohedbergella portsdownensis, a cold-water species that occurs at the same level as the Boreal macrofossils. These observations indicate a palaeoceanographic scenario characterized by reduced stratification of surface waters and absence/disruption of the thermocline in a dominantly eutrophic regime during MCE 1. Evidence provided by planktonic foraminifera, Boreal macrofossils and δ18O records documented for the late Cenomanian Plenus Cold Event (PCE) at Eastbourne (UK) reveal similarities that confirm the periodic inflow of cold Boreal seawater originating in the Norwegian Sea as previously postulated to explain the occurrence of Boreal fauna in the Anglo-Paris Basin. The southerly extension of this water mass may be related to the reorganization of circulation driven by the long eccentricity cycle.

Palaeontological field surveys in remote regions are a challenge, because of uncertainty in finding new specimens, high transportation costs, risks for the crew and a long time commitment. The effort can be facilitated by using high-resolution satellite imagery. Here we present a new opportunity to investigate remote fossil localities in detail, mapping the optical signature of individual fossils. We explain a practical workflow for detecting fossils using remote-sensing platforms and cluster algorithms. We tested the method within the Petrified Forest National Park, where fossil logs are sparse in a large area with mixed lithologies. We ran both unsupervised and supervised classifications, obtaining the best estimations for the presence of fossil logs using the likelihood and spectral angle mapper algorithms. We recognized general constraints and described logical and physical pros and cons of each estimated map. We also explained how the outcomes should be critically evaluated with consistent accuracy tests. Instead of searching for fossiliferous outcrops, our method targets single fossil specimens (or highly condensed accumulations), obtaining a significant increase in potential efficiency and effectiveness of field surveys. When repeatedly applied to the same region over time, it could also be useful for monitoring palaeontological heritage localities. Most importantly, the method here described is feasible, easily applicable to both fossil logs and bones, and represents a step towards standard best practices for applying remote sensing in the palaeontological field.

This contribution provides insights into ocean-floor hydrothermal metamorphism of the fast-evolving Dinaridic Neotethys. Mineralogical, geochemical and Sr isotope data collected from altered ophiolites and non-ophiolite basalts/andesites and tuffs of the active continental margin are consistent with hydrothermal alteration trajectories that reflect the host-rock composition. This suggests that hydrothermal fluxes were restricted within a simple closed seawater-fed system. Based on the initial isotopic ratios of Sr, two fluid–rock interaction trends are established: (a) low-to-medium degrees of metasomatism in pre-Middle Jurassic anorogenic ophiolites that progressively abated, and (b) increased intensities of metasomatism in post-Middle Jurassic orogenic ophiolites. This agrees with chlorite thermometry and Ca-Al-(Fe)-silicate phase chemistry. The metamorphic assemblages belong to the zeolite, prehnite-pumpellyite, prehnite-actinolite and greenschist facies. The facies is reliant on the temperature of hydrothermal systems and their fluid chemistry. Rare earth element (REE) phase geochemistry shows (a) variable fluid–rock ratios in chlorite and pumpellyite dependent on fluid temperatures, (b) prominent Eu and Ce anomalies that reflect the fluid oxidation state, (c) light REE/heavy REE mobilization attributed to prevalent ligand complexation, and (d) multi-phase fluid percolation across reaction zones of heterogeneous permeability. This study proposes initiation of simple hydrothermal system(s) at or near a spreading centre(s) in the infancy of the Dinaridic Neotethys. Such a system became more complex during Middle Jurassic and Early Cretaceous time with reactive hydrothermal fluids passing the recharge area and reaching the hot reaction zone. An abrupt obliteration of the established high-temperature regime ensued, following the final closure of the Neotethys.

Brittle reactivation of plastic shear zones is frequently observed in geologically old terranes. To better understand such deformation zones, we have studied the >700 Ma long structural history of the Himdalen–Ørje Deformation Zone (HØDZ) in SE Norway by K–Ar and 40Ar–39Ar geochronology, and structural characterization. Several generations of mylonites make up the ductile part of HØDZ, the Ørje Shear Zone. A 40Ar–39Ar white mica plateau age of 908.6 ± 7.0 Ma constrains the timing of extensional reactivation of the Ørje mylonite. The mylonite is extensively reworked during brittle deformation events by the Himdalen Fault. 40Ar–39Ar plateau ages of 375.0 ± 22.7 Ma and 351.7 ± 4.4 Ma from pseudotachylite veins and K–Ar ages of authigenic illite in fault gouge at c. 380 Ma are interpreted to date initial brittle deformation, possibly associated with the Variscan orogeny. Major brittle deformation during the Early–Mid Permian Oslo Rift is documented by a 40Ar–39Ar pseudotachylite plateau age of 294.6 ± 5.2 Ma and a K–Ar fault gouge age of c. 270 Ma. The last datable faulting event is constrained by the finest size fraction in three separate gouges at c. 200 Ma. The study demonstrates that multiple geologically significant K–Ar ages can be constrained from fault gouges within the same fault core by combining careful field sampling, structural characterization, detailed mineralogy and illite crystallinity analysis. We suggest that initial localization of brittle strain along plastic shear zones is controlled by mechanical anisotropy of parallel-oriented, throughgoing phyllosilicate-rich foliation planes within the mylonitic fabric.

Conodonts are an important biostratigraphic tool for many Phanerozoic stages. Along with graptolites, they define all global Ordovician Stage boundaries. Within the Upper Ordovician interval, a known species of Amorphognathus tvaerensis (Bergström) is present in both Sandbian and Katian stratotype sections. Study of changes in the succession of A. tvaerensis revealed that elements in the upper part of its range differ morphologically quite distinctly from those in its lower part. Here, they are described as a new conodont species, A. viirae sp. nov. This new species is recognized in several Estonian and Swedish sections, with apparent occurrence also in Mójcza Quarry, Holy Cross Mountains, Poland and Black Knob Ridge, Oklahoma, USA. Detailed analysis of early Amorphognathus elements from Estonian and Swedish sections revealed the absence of A. inaequalis (Rhodes) in both regions, although a conodont subzone based on this species was identified earlier by some authors. Both the absence of A. inaequalis (Rhodes) and recognition of the new species A. viirae sp. nov. resulted in the revision of the conodont zonation, and a new version of it is proposed for the Sandbian Stage in the Atlantic Realm. The new zonation includes (from below) Pygodus anserinus, Baltoniodus variabilis, A. tvaerensis, B. gerdae, A. viirae and B. alobatus Conodont zones.

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.

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.