# Role of Hydrothermal Circulation along and above Inherited Basement Structures Relating to Unconformity-Related Uranium Mineralization

Jonathan Poh, Khalifa Eldursi, Patrick Ledru, Philippe Yamato, Guoxiang

Chi, Antonio Benedicto

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jonathan Poh, Khalifa Eldursi, Patrick Ledru, Philippe Yamato, Guoxiang Chi, et al. Role of Hydrothermal Circulation along and above Inherited Basement Structures Relating to Unconformity-Related Uranium Mineralization. Geofluids, 2022, 2022, pp.Article ID 9131289. 10.1155/2022/9131289. insu-03844761

HAL Id: insu-03844761

https://insu.hal.science/insu-03844761

Submitted on 9 Nov 2022

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research Article

# Role of Hydrothermal Circulation along and above Inherited Basement Structures Relating to Unconformity-Related Uranium Mineralization

Jonathan Poh ID,<fnref n="1" /> Khalifa Eldursi ID,<fnref n="2" />,<fnref n="3" /> Patrick Ledru ID,<fnref n="3" /> Philippe Yamato ID,<fnref n="1" />,<fnref n="4" />
Guoxiang Chi ID,<fnref n="5" /> and Antonio Benedicto ID,⁶

Correspondence should be addressed to Jonathan Poh; jonathanpoh87@gmail.com

Received 18 February 2022; Revised 25 September 2022; Accepted 5 October 2022; Published 7 November 2022

Academic Editor: Dan Ma

Copyright © 2022 Jonathan Poh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

The reactivation of inherited tectonic structures formed during the Paleoproterozoic Trans-Hudson Orogeny (THO) has played a significant role in generating high-grade unconformity-related uranium deposits in the eastern Athabasca Basin. The role of these tectonic structures is now investigated through a series of two-dimensional hydrothermal numerical models. Two modelling scenarios are considered: (1) models during the THO peak of metamorphism and (2) models with a permeable layer mimicking the presence of the Athabasca Basin, deposited unconformably over the THO basement. In the first scenario, general fluid patterns are strongly affected by the applied permeability configurations. Unidirectional high fluid flow zones (from 10⁻⁹ to 10⁻⁸ m·s⁻<fnref n="1" />) and high thermal gradients (up to 65°C·km⁻<fnref n="1" />) can be observed above and within the deep-seated tectonic structures. In the second scenario, well-established fluid convection cells or unidirectional fluid flow zones are observed within the basin layer, with upflow originating from the core of the deep-seated structures, regardless of the applied fluid pressure regime. These results highlight that these deep-seated structures can efficiently transport fluids and heat towards the upper parts of the crust and the basin. In the second scenario, the loci for preore alteration are then evaluated by computing a rock alteration index based on temperature and fluid velocity constraints. These alteration areas reside along and above the deep-seated structures and are potential regions for structural reactivation during mineralization. These results imply that the analysis of the inherited tectonic structures, combined with the alteration regions, can serve as markers for uranium exploration.

## 1. Introduction

The formation of mineral deposits requires multiple processes related to hydrothermal systems. According to the “mineral system” framework, four critical elements must be defined to evaluate the conditions to generate mineral deposits: lithosphere architecture, favourable geodynamics and transient self-organized hydrothermal systems, the fertility of the geological domains, and the preservation of the

depositional environment [1, 2]. This framework indicates that large long-lived, deeply penetrating, and steeply dipping structures that juxtapose distinctly different basement domains are geological features that lead towards giant deposits [1]. Several studies have evaluated how the “mineral system” framework can be applied from the regional to the camp scale and translated in terms of project generation and exploration targeting [1, 3, 4]. This concept is currently guiding geophysical surveys, highlighting the limits of

<footnote n="1">¹Univ Rennes, CNRS, Géosciences Rennes, UMR 6118, 35000 Rennes, France</footnote>

<footnote n="2">²Department of Earth Sciences, University of Bashqi, PO Box 9480-1308, Bashqi, Libya</footnote>

<footnote n="3">³Université de Lorraine, GeoRessources 2, Rue du Doyen Marcel Roubault TSA 70605 54518 Vandoeuvre-lés-Nancy Cedex, France</footnote>

<footnote n="4">⁴Institut Universitaire de France, Paris, France</footnote>

<footnote n="5">⁵Department of Geology, University of Regina, Regina, Saskatchewan, Canada S4S 0A2</footnote>

<footnote>Université Paris Saclay, UMR GEOPS, 91405 Orsay, France</footnote>
crustal and lithospheric blocks [5, 6], the reassessments of the potential of regional metallogenic provinces (e.g., [3, 7]), and the detailed mineralogical and geochemical studies of the successive hydrothermal alteration episodes [8–12]. However, the contribution of numerical modelling to constrain the role of such deeply penetrating and steeply dipping structures as long-lived vectors for fluids and heat transport has not been investigated.

The Athabasca Basin (Saskatchewan, Canada), a Paleo- to Mesoproterozoic siliciclastic basin, hosts giant unconformity-related uranium deposits and a major metallogenic province [13]. At the eastern extent of the Athabasca Basin, multiple high-grade and high tonnage unconformity-related uranium deposits reside along a NE-SW trending basement-hosted structure that extend over several hundred kilometres along strike. This is the Wollaston-Mudjatik Transition Zone (WMTZ) formed during the Trans-Hudson Orogeny (THO; Figure 1(a); [14, 15]). As revealed by exploration drill holes and its magnetic geophysical signature below the Athabasca Basin, this tectonic structure is composed of steeply dipping anastomosed shear zones (Figure 1(b), [16–18]). At its western extent, recent discoveries of basement hosted world class uranium deposits along the Patterson Lake conductive trend [19] and extensive geophysical and geological surveys [6, 20–22] reveal the presence of crustal to lithospheric scale deeply rooted structures. These structural features can be considered as critical elements for the location of major hydrothermal systems and generation of giant uranium deposits.

At the unconformity-related uranium deposits in the eastern Athabasca Basin, uranium mineralization has been related to the steeply dipping inherited basement structures, considering their ability to move large volumes of fluid (e.g., [23–26]). However, the mechanism for fluid transport and metal deposition remains unclear. Li et al. [25] argued that short-term deformation-driven flows might be incapable of transporting sufficient metalliferous-enriched fluid needed to form large unconformity-related uranium deposits. They proposed that thermally driven fluid convection is more capable of transporting large volumes of uranium-bearing fluids to deposition sites and that alternating deformation-driven fluid flow and thermally driven fluid convection are responsible for uranium mineralization. The establishment of convection cells at thermal equilibrium is long-lasting, even though fault reactivations play an important role in enhancing the permeability of the conduits.

This numerical modelling study is aimed at addressing the following questions concerning fluid flow patterns and physical conditions required to form uranium deposits within and above deep-seated structures during periods of quiescence:

(1) What were the fluid flow patterns and their main drivers during the metamorphic retrogression of the Wollaston-Mudjatik Transition Zone?

(2) What was the role of the ancient deep-seated inherited structures and how the pre-Athabasca fluid flow pattern and related alteration impacted or controlled

the location of the hydrothermal activity at the origin of the unconformity-related uranium deposits?

(3) What are the implications of deep-seated structures for future ore exploration?

For this study, we consider two 2-dimensional modelling scenarios that include inherited tectonic structures formed before the formation of the Athabasca Basin. The first modelling scenario mimics the tectonic conditions following the evolution of the WMTZ during the peak of metamorphism of the THO. The performed models allow to test different permeability configurations of increasing complexity before the deposition of the Athabasca Basin. The second modelling scenario involves a permeable layer representing the presence of the Athabasca Basin. A rock alteration index constrained by fluid velocity and temperature distribution [27] is then computed to predict zones where preore chlorite alteration is most likely to occur when the second modelling scenario has reached steady-state conditions.

## 2. The Unconformity-Related Uranium Mineral System of Eastern Athabasca Basin

The crustal architecture of the eastern part of the Athabasca Basin is inherited from the Paleoproterozoic Trans-Hudson Orogeny (THO; Figure 1) which resulted from the collision between the Hearne craton, the Sask craton, and the Superior craton between ca 1.87 and 1.79 Ga [15, 28–34]. The Hearne craton, consisting of Archean granulitic tonalitic-trondhjemitic gneiss domes of the Mudjatik Domain and the Paleoproterozoic metasedimentary gneisses of the Wollaston Domain, is interpreted as a passive margin sequence [31]. These domains were intensively imbricated and deformed during the Trans-Hudson Orogen that led to the formation of the WMTZ [14, 28, 32]. The WMTZ is characterized as subvertical anastomosing fault systems that are highlighted by narrow magnetic low signatures (Figures 1 and 2(a); [16–18, 35]). Mylonite zones are characterized by steeply dipping foliations bearing subvertical stretching lineation [36, 37] and by the presence of Paleoproterozoic sediments that reached amphibolite metamorphic facies conditions (i.e., $P \approx 0.8–1$ GPa, $T \approx 700–800°C$, [14, 15, 38]).

Previous thermomechanical numerical models have verified that such deep-seated major shear zones could be formed from a protracted shortening of the lithosphere [39] and that sediments can reach burial $P-T$ conditions in agreement with the peak of metamorphism known for the THO at ca. 1.84–1.82 Ga ([39] and Figure 2(b)). Thermobarometry studies indicate a decrease in pressure from 0.8–1 GPa to ~0.5 GPa occurring at 1.82 to 1.77 Ga due to exhumation of the Archean basement and changes in the tectonic regime [15, 38]. Subsequently, a protracted cooling period occurred in the region at 1.76–1.71 Ga, marking the late-THO tectonic evolution [28, 38, 40–42].

Tectonic evolution of the THO is finally sealed by the deposition of detrital sediments that form the Paleoproterozoic to Mesoproterozoic basins (including Athabasca) which rest unconformably over the exhumed domains, preserving
<figure><img src="image_3.png" /></figure>

FIGURE 1: (a) Generalised tectonic map of the Rae-Hearne craton, Western Canadian shield, with named unconformity-related uranium deposits and main structural features that has been modified from Annesley et al. [14], Alexandre et al. [28], and Eglington et al. [107]. A composite interpreted aeromagnetic map (provided by Orano Canada Inc) is superimposed on the map as an inset for (a), highlighting the different structural types located in the Mudjatik Domain and Wollaston-Mudjatik Transition Zone. (b) Interpreted structural features of the basement along the Wollaston-Mudjatik Transition Zone with named unconformity-related uranium deposits and occurrences (modified from [16, 24]). Dashed lines indicated main structural features. The solid red line (A-A') within inset represents a schematic cross-section of the Wollaston-Mudjatik Transition Zone in Figure 2(a).

the underlying basement-hosted structures. These steeply dipping tectonic structures within the basement and in particular the graphite-bearing shear zones were reactivated after the deposition of the Athabasca Basin. This reactivation is marked in many fault segments by an offset of the unconformity that reaches about 100 m in the McArthur River uranium deposits [13, 16] and by the development of cataclastic breccia and gouge in the basement faults and brittle faults into the basin (e.g., [13, 43-47]).

This polyphase tectonic evolution and the existence of favourable transient events are also considered as critical elements to characterize the unconformity-related uranium mineral system [16]. The continuity between the reactivated faults of the basement and the brittle faults in the basin and the development of alteration haloes and mineralization within these structures suggest that the reactivation has controlled the circulation of the mineralizing fluids ([13, 16, 48]). Geochemical and fluid inclusion studies indicate that

uranium mineralization is coeval with the tectonic reactivation of deep-seated inherited structures [26, 49-52]. Previous hydrothermal numerical results indicate that mineralization can be associated with tectonic reactivation occurring after the thermal blanketing effect generated by the Athabasca Basin [24, 25, 53-55]. The way how fluids circulate between the basement and basin through the inherited structures leads to two types of unconformity-related uranium deposits characterized by different alteration zones [13, 49, 56-58]. The first type is the sandstone-hosted type associated with the discharge (egress) of the fluids from the basement. The sandstone-hosted type deposits have two end-member alteration patterns that are located within the sandstone overlying the unconformity: (1) zones of quartz dissolution+illite and quartz overgrowth+illite-kaolinite-chlorite+dravite as observed at the Cigar Lake deposit [43] and (2) zones of silicified-dravite-chlorite-kaolinite and illite-dravite as characterized in the McArthur River deposit
<figure><img src="image_4.png" /></figure>

FIGURE 2: (a) Generalised cross-section of the Wollaston-Mudjatik Transition Zone with the location marked in Figure 1(a) inset. Modified from [16] and Benedicto et al. [35]. Thermomechanical model outputs, (b) M1 and (c) M2, from Poh et al. [39] are used as inputs to generate reference models for fluid-thermal simulation shown in Figure 3. Red rectangular boxes in (b) and (c) show the physical dimensions of the hydrothermal models which are conducted in this study, respectively. The geometry and spacing of the anastomosing fault systems in (a) appear to be in good agreement with the vertical structures generated in the thermomechanical models in (b) and (c).
[46]. The second type is the basement-hosted deposits resulting from recharge (ingress) of the fluids from the basin. Narrow, inverted alteration halos of sudoite±illite and chlorite+biotite±sudoite are located at the outer margins of the inherited basement structures [58]. However, the spatial distribution of the alteration halos is extensive, up to 400 m wide at the base of the sandstone, and it may exceed several thousand metres in strike length (e.g., [13]). These unconformity-related uranium deposits and their alteration zones attest to the widespread fluid-rock interactions and fluid circulation between the major fault systems of the basement and the Athabasca Basin.

# 3. Modelling Approach

## 3.1. Governing Equations.
The hydrothermal models are performed with FLAC3D (the Fast Lagrangian Analysis of Continua in 3 Dimensions). This code solves the coupled fluid flow and heat transfer equations on a Cartesian grid via an iterative process that involves sequential stepping of fluid and thermal modules [59]. It has been widely used for modelling fluid flow, heat transfer, and deformation, especially at an outcrop to upper crustal scale [24, 53–55, 60–63].

Heat transfer is solved by the following series of equations that define the energy balance in a convective-diffusive heat setting. Conservation of energy (Equation (1)) is achieved through

$$C_p^r \frac{\partial T}{\partial t} + \nabla \mathbf{q} + \rho_0 C_p^w \mathbf{q}_w \cdot \nabla T = q_v, \quad (1)$$

where T is temperature; q is the thermal flux; q_w is the fluid specific discharge; q_v is the volumetric heat source intensity; ρ_0 and C_p^w are the reference density and specific heat capacity of the fluid, respectively; and C_p^r represents the specific heat capacity of the rocks.

Conductive heat transfer (Equation (2)) is solved by Fourier's law which defines specific heat flow (q) as heat flow normalised by area:

$$\mathbf{q} = -k_r \nabla T, \quad (2)$$

where k_r is the thermal conductivity of rocks and ∇T is the temperature gradient vector.

The temperature values from the previous time step (Equation (3)) are used to calculate the fluid density (ρ_f) in the following time steps, which is expressed as

$$\rho_f = \rho_0 [1 - \alpha_f (T - T_0)], \quad (3)$$

where α_f is the volumetric thermal expansion of the fluid at T = T_0. The calculated fluid density (ρ_f) is used to compute fluid velocity (q_w) according to Darcy's Law (Equation (4)), which is expressed as

$$\mathbf{q}_w = -\frac{k}{\mu} \nabla P + \rho_f g \nabla z, \quad (4)$$

where μ is the fluid dynamic viscosity, g is the gravitational acceleration, k is the assigned permeability of the rock, and ∇P is the fluid gradient pressure.

Boussinesq approximation is applied to the models.

## 3.2. Model Setup

### 3.2.1. Geological Setting for Modelling Sections.
The first main geological scenario, M1 (Figure 2(b)), mimics the evolution of the WMTZ during the peak of metamorphism of the THO. The model consists of two homogenous layers of rocks simulating the presence of metasediments (representing the Wollaston Supergroup) deposited unconformably on top of the Archean granitic basement. The geometry of the metasediments that intersects the basement is determined by the sedimentary cusps formed due to pop-down tectonics [36, 39, 64, 65]. From M1, three modelling subscenarios will be evaluated as M1A, M1B, and M1C (see Section 3.2.3 for more details).

The second main geological scenario, M2 (Figure 2(c)), takes into account a period of tectonic quiescence of 120–160 Ma occurring between the peak of metamorphism of the THO (~1.84 Ga; [14]) and the deposition of the Athabasca Basin (~1.76–1.72 Ga; [38, 42]). During this period, the eroded region undergoes a slow postorogenic cooling [41]. The region underwent 20 km of distributed exhumation leading to an uplift of the basement up to the depth of 7 km. This erosion phase removed most of the complex structural features, leaving behind a remnant of the deep-seated structures starting at the depth of 6 km. The resultant geological model agrees with previous seismic surveys and thermotectonic interpretations [14, 66, 67]. From M2, three modelling subscenarios will be evaluated as M2A, M2B, and M2C (see Section 3.2.3 for more details).

### 3.2.2. Geometry, Boundary, and Gridding of the Models.
The upper and lower boundaries are left open for fluids to leave the system, whereas the side boundaries remain closed and impermeable. Initial fluid pressure at the bottom boundary is set to the calculated fluid pressure (see Section 3.2.3 for more details). Initial temperature at the top of the model is maintained at 20°C. At the bottom boundary, the initial temperature is assigned with a value based on a linear geothermal gradient of 30°C·km⁻¹.

All M1 models have a physical model domain of 200 km length, 0.1 km width, and 15 km depth (see Figures 3(a)–3(c)). For M1A and M1B in Figures 3(a) and 3(b), the structural model was designed using an equidistant rectangular grid with a cell resolution of 1 km by 0.15 km (20000 elements; Figure 4(a), A1 and A2). The width of the vertical structures penetrating the basement is 3 km. For M1C, a Tetrahedral Irregular Network (TIN) was used to generate a more accurate representation of the geometry for the structures and deformed regions (Figures 3(b) and 4(b)). With 17710 elements, high-resolution zones are assigned to where we expected to have high hydrothermal activity (i.e., within the metasediments and tectonic structures) to reduce computation time (Figure 4, B1 and B2). The size of the grid
<figure><img src="image_5.png" /><figcaption>FIGURE 3: Permeability configurations for two types of models: M1 and M2. Model M1 has three permeability configurations denoted as (a) M1A, (b) M1B, and (c) M2C. Permeability configuration for models (c) M2A, M2B, and M2C remains the same. Details regarding permeability setups for models are found in Table 1 and Section 3.2.2.</figcaption></figure>

for low-permeability regions such as the basement is assigned a coarser grid (up to 2 km spacing in Figure 4, B2).

To focus our study on the hydrothermal activity involving the basin and the surviving deep-seated structures, the length of M2 models is reduced from 200 km to 110 km to provide focus in our study. The same TIN grid system of 14514 elements is also applied (Figure 4(c)). Between the depths of 3 to 8 km, the high-permeability zones (i.e., basin and deep-seated structures) are refined to have high-resolution regions (tetrahedral grid smaller than 200 m, seen in Figure 4, C1 and C2). Low-resolution grids (i.e., tetrahedral grid larger than 1.5 km) are employed to the basement and deformed basement rocks as both represent permeability regions.

### 3.2.3. Model Fluid Configuration.
The hydrothermal modeling protocol follows a simple physical procedure that does not consider vapour-fluid multiphase, fluid-fluid, and fluid-rock interaction. Among the fluid properties, only density varies with temperature. Other parameters such as fluid viscosity, heat capacity, and thermal conductivity remain constant (see Table 1). According to Sibson [68], fluid pressure increases with depth following different regimes (see the blue dashed line in Figure 5). The upper part of the profile (from 0 to 4 km depth) reflects the hydrostatic pressure conditions where the pore fluid factor, $\lambda_v = P_f/\sigma_v \approx 0.4$ (where $P_f$ and $\sigma_v$ are the fluid pressure and the vertical stress, respectively). Below, the fluid pressure profile corresponds
<figure><img src="image_6.png" /><figcaption>FIGURE 4: Different types of grids (Equidistant Rectangle and Tetrahedral Irregular grid) are used for different models. While Equidistant Rectangle grid is used for simple models (two rock layer models; (a) M1A and M1B), the Tetrahedral Irregular grid is employed to design complex models ((b) M1C and (c) M2A, M2B, and M2C). High-resolution regions are assigned to rock layers such as the basin and deep-seated structures, where hydrothermal activity is expected. Subsequently, low-permeability regions (i.e., the basement and the impermeable cover) are assigned as low resolution.</figcaption></figure>
<table><caption>TABLE 1: Hydrothermal properties of rocks and fluid used in the models.</caption><thead><tr><th></th><th>Metasediments<br>(sedi)</th><th>Basement<br>(base)</th><th>Deformed metasediments<br>(def_sedi)</th><th>Deformed basement<br>(def_base)</th><th>Pure water</th><th>Lower Athabasca Basin<br>(basin)</th><th>Impermeable cover<br>(cover)</th></tr></thead><tbody><tr><td><strong>Hydrological properties</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>φ (%)</td><td>10.0</td><td>5.0</td><td>20.0</td><td>10.0</td><td>—</td><td>10.0</td><td>5.0</td></tr><tr><td>k (m<sup>2</sup>)</td><td>M1A: 10<sup>-16</sup><br>M1B: 10<sup>-13</sup> to 10<sup>-17</sup><br>M1C: 5.0 × 10<sup>-17</sup></td><td>10<sup>-18</sup></td><td>M1C: 10<sup>-16</sup><br>M2: 5.0 × 10<sup>-15</sup></td><td>M1C: 7.0 × 10<sup>-16</sup><br>M2: 10<sup>-17</sup></td><td>—</td><td>10<sup>-13</sup></td><td>10<sup>-17</sup></td></tr><tr><td>ρ (kg·m<sup>-3</sup>)</td><td>2600.0</td><td>2650.0</td><td>2400.0</td><td>2400.0</td><td>1000.0</td><td>2500.0</td><td>2500.0</td></tr><tr><td>η (Pa·s)</td><td>—</td><td>—</td><td>—</td><td>—</td><td>1.0 × 10<sup>-3</sup></td><td>—</td><td>—</td></tr><tr><td><strong>Thermal properties</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>C<sub>p</sub> (J·kg<sup>-1</sup>·K<sup>-1</sup>)</td><td>803.0</td><td>803.0</td><td>803.0</td><td>803.0</td><td>4185.0</td><td>803.0</td><td>803.0</td></tr><tr><td>k<sub>f</sub> (W·K<sup>-1</sup>·m<sup>-1</sup>)</td><td>2.5</td><td>3.5</td><td>M1: 2.5<br>M2: 5.0</td><td>3.5</td><td>0.6</td><td>2.5</td><td>2.5</td></tr><tr><td>α (K<sup>-1</sup>)</td><td>10<sup>-5</sup></td><td>8.0 × 10<sup>-6</sup></td><td>10<sup>-5</sup></td><td>8.0 × 10<sup>-6</sup></td><td>1.85 × 10<sup>-3</sup></td><td>10<sup>-5</sup></td><td>8.0 × 10<sup>-6</sup></td></tr></tbody></table>
<figure><img src="image_7.png" /></figure>

FIGURE 5: Profiles of fluid pressure versus depth within the crust depending on the pore-fluid factor, $λ_v$. The points on the dotted-dashed line are taken from the M2 reference model. The blue dashed line assumes a smooth progression from hydrostatic towards lithostatic pressure as permeability steadily reduces with increasing depth from Sibson [86]. Rock density values used to compute the curve can be referred in Table 1.

to the transition from a hydrostatic to a lithostatic regime. This regime (where $0.4 < λ_v < 0.9$) can extend up to 15 km depth (see [68]). However, below 10 km depth, fluid pressure is very close to lithostatic pressure ($λ_v > 0.9$) without considering other potential horizontal stresses.

In our models, a simplified dual-pressure regime is considered (black dashed line in Figure 5), as seen in previous published numerical models (e.g., [69]). Hydrostatic pressure is first applied up to 7 km depth ($λ_v < 0.4$) in order to account for the porous sandstone-dominated nature of the Athabasca strata that form a relatively permeable environment [49, 50, 70, 71]. Fluid pressure increases to near lithostatic pressure from 7 to 8 km depth allowing to simulate the transition from hydrostatic to lithostatic pressure.

For the M2 models in Figure 4(c), three variants of fluid pressure regimes are considered. M2A uses the same simplified dual-pressure regime as in M1. The second model, M2B, follows the same simplified dual-pressure system of M2A but fluid pressure evaluated at the bottom of the model is kept constant. It allows to simulate the presence of a constant fluid flow maintained throughout the more permeable regions (i.e., deformed basement and deep-seated structures). M2C uses a simplified hydrostatic pressure regime only (i.e., $λ_v = 0.4$) to evaluate the effect of fluid pressure change in fluid flow patterns.

### 3.2.4. Model Permeability Configurations.
Estimating rock permeability ($k$) in any geologic model is required for fluid flow modelling but assigning permeability values to rocks is challenging, considering their wide range of possible

values and their sensitivity to physical processes (e.g., [72, 73]). Hence, a suitable permeability value is required primarily for large and deep-seated structures (10 km by 15 km) to establish a baseline of the fluid flow patterns. Previously published numerical works simulated large and deep fault zones with a constant permeability of $10^{-15} \text{ m}^2$ that extends to depths of 15 km (e.g., [27, 61, 74]). High fluid flux values in our initial numerical tests of $10^{-15} \text{ m}^2$ for the deep-seated structures cause fluid flow values not to converge, resulting in the model being numerically unstable. The permeability of the deep-seated structures ($k_{\text{def\_base}} = 10^{-16} \text{ m}^2$) is then decreased to allow the fluid flow and temperature equations to converge, thus achieving a thermally steady-state condition.

In M1 subscenarios, permeability of the rocks is first varied to evaluate their effect on the fluid flow patterns. The porosity and permeability assigned to the various rocks remain constant regardless of model scenario (Table 1). Three permeability configurations are considered. The first subscenario (M1A in Figure 3(a)) has a bimodal permeability configuration with a $k_{\text{sed1}}$ of $10^{-16} \text{ m}^2$ for the metasediments and a $k_{\text{base}}$ of $10^{-18} \text{ m}^2$ for the basement. The second subscenario (M1B in Figure 3(b)) utilises a simplified depth-dependent permeability function for metasediments, as indicated by experimental data [75–77]. The depth-dependent function is expressed as $\log k = -14 - 3.2 \log z$, where $k$ ($\text{m}^2$) and $z$ (km) correspond to rock isotropic permeability and depth, respectively. The permeability of the basement is kept constant. High-strain regions found within the metasediments and basement are incorporated into the
third subscenario (M1C in Figure 3(c)) to simulate tectonic inheritance (M1 in Figure 2(b)). The spacing and geometry of the inherited zones are then in good agreement with the anastomosed fault systems observed nature (in Figure 2(a)). These high-strain regions in the models for both metasediments and basement are labelled as deformed metasediments and deformed basement, respectively. The permeability of these deformed metasediments ($k_{\text{def\_sedi}} = 10^{-16} \text{ m}^2$) and basement ($k_{\text{def\_base}} = 7.0 \times 10^{-16} \text{ m}^2$) is assumed to be more permeable than their nondeformed counterparts that mimic the damage zones during fault generation (e.g., [78, 79]). The change of permeability from their nondeformed counterparts is highly uncertain due to the heterogeneity of the fault zone and sampling process [79]. The permeabilities of the deformed metasediments and basement are considered constant.

The Athabasca Basin in M2 is composed of a 2 km permeable homogenous sandstone layer ($k_{\text{basin}} = 10^{-13} \text{ m}^2$) followed by a 4 km low-permeability eroded cover ($k_{\text{cover}} = 10^{-17} \text{ m}^2$) that represents the clay-rich impermeable Wolverine Point Formation overlying the Manitou Falls Formation (Figure 3(d); [42]). The thermal conductivity of the remaining deep-seated structures (in particular, the exhumed deformed metasediments cover located within the basement) was increased (from 2.5 to 5.0 W·K⁻¹·m⁻¹) to reflect the presence of alteration minerals (i.e., graphite, quartz, and chlorite) that have been identified as vectors for uranium uptake and precipitation [49, 55, 80–83].

### 3.3. Thermally Constrained Restricted Rock Alteration Index (RRAI)
Hydrothermal fluids contain dissolved chemical species that flow through rock columns. Tectonic events that lead to strain localisation perturb the environment via permeability enhancement, which allows the transport of reactants in solution towards the reaction site by advection and diffusion (e.g., [73, 84–86]). Potential alterations and mineralization patterns can be deduced by looking at zones where the cooling rate is essential. This cooling rate is referred to as the “rate of mineralisation” [87] that corresponds to the negative values of the rock alteration index (RAI). RAI is evaluated as $\vec{u} \cdot \Delta T$, where u and T represent fluid velocity and temperature, respectively. RAI has been used in previous hydrothermal models to evaluate and predict ore deposition potential (e.g., [63, 74, 88]). RAI values highlight the importance of the deep-seated structures concerning unconformity-related uranium deposits. Negative RAI indicate fluid flows to regions of low temperature. Potential zones for mineralization and alteration can be predicted in a metamorphic setting when $u > 10^{-10} \text{ m} \cdot \text{s}^{-1}$ (see [27], for more details). This fluid velocity value is then applied as a constraint to RAI which we now consider as the restricted rock alteration index (RRAI).

Alteration halos delineate the presence of these unconformity-related deposits following the main deposit types. Preore chlorite alteration is then considered as a pathfinder and a potential chemical trap for uranium mineralization [9, 28, 89]. We propose that RRAI can be constrained further by associating it with the formation temperature of

the preore chlorite alteration. These temperatures are found to vary widely in deposits: from 275 to 300°C in the Millenium deposit [9]; 230–250°C at Dawn Lake, McArthur River, and Rabbit Lake [89]; and 235–245°C in McArthur River [11]. Therefore, in combination with RRAI and the formation temperature of preore alteration, this RRAI certifies that the physical conditions for chlorite alteration are met; RRAI = 1 if fluid temperature is within 230–300°C and RAI is negative and $u > 10^{-10} \text{ m} \cdot \text{s}^{-1}$. This thermally constrained RRAI will only be applied to the M2 models because formation temperatures of preore alterations are found in deposits located during the presence of the Athabasca Basin and the metasediments of the WMTZ, and the earliest uranium mineralization event occurred with the presence of the Athabasca Basin (ca. 1.68 Ga; [28]).

### 3.4. Model Limitations
The primary assumption in our numerical models is that the fluid viscosity is kept constant. In FLAC3D, rock permeability interacts with the rock's isotropic permeability coefficient and is evaluated as $k = k_h/\rho_g g$, where $k_h$ and $g$ are hydraulic conductivity and gravitational acceleration, respectively. Furthermore, intrinsic permeability ($k$) is evaluated to be inversely proportional to dynamic viscosity ($\mu$). Making fluid viscosity variable will cause unwanted permeability enhancement to the initial permeability setup seen in Sheldon [61]. In the context of the unconformity-related uranium deposits at the Athabasca Basin, the source of the basement-associated fluids is highly debatable (e.g., [12]). Considering the complex geological layout from the thermomechanical model, we can only investigate the fluid flow patterns by keeping fluid viscosity constant and preventing any permeability enhancement. Variations in the fluid viscosity can be considered in future model iterations.

## 4. Results

### 4.1. Model M1
In the bimodal permeability configuration for M1A, macroscopic fluid flow patterns indicate unidirectional fluid flow within the metasediments. Fluid is observed to flow along the contact between the impermeable basement and in regions where sedimentary cusping occurred at $x = -40 \text{ km}, -20 \text{ km}, 0 \text{ km}$, and ca. +25 km. No convection cells are observed (Figure 6(a)). Fluid temperature in the metasediments reached at least 500°C where sedimentary cusps are present and 475°C when they are not present at depth. Three regions of increased fluid flow velocity of at least an order of magnitude ($u \approx 10^{-11} \text{ m} \cdot \text{s}^{-1}$ to $10^{-10} \text{ m} \cdot \text{s}^{-1}$) in comparison with the basement are located: (1) within the metasediments where sedimentary cusps are present at depth, (2) below the contact of the metasediments and the basement, and (3) within the basement at ~13 km depth. This region of increased fluid velocity corresponds to the transition from hydrostatic to lithostatic pressure (as indicated in Figure 5) and is present in all M1 models (Figure 6).

The presence of a variable permeability in metasediments of model M1B enhances a more variable and periodic fluid flow. Upward flows occur where sedimentary cusps are present (at $x \approx -40 \text{ km}, -20 \text{ km}, 0 \text{ km}$, and +25 km in
<figure><img src="image_8.png" /><figcaption>FIGURE 6: Thermal and fluid flow results for models (a) M1A, (b) M1B, and (c) M1C. Permeability setups for each model can be referred to Figure 3. Isotherms are represented by solid-coloured lines in increments of 50°C. Colour fill of the plots and isotherms is based on scientific colour maps [108].</figcaption></figure>

Figure 6(b)), and downward flows occur aside them (at $x \approx -50$ km, -30 km, +10 km, and +60 km in Figure 6(b)). The pairing of upward and downward flows coupled with the corresponding deflections of the isotherms indicates the presence of well-established convection cells.

Unidirectional fluid flow vectors in M1C (where $k_{\text{sedi}}$ in M1C is reduced by an order of magnitude) indicate that no convection cells form in these conditions (Figure 6(c)). This, compared with model M1A (Figure 6(a)), prevents the generation of high temperatures at the upper parts of the model. Fluid velocities show regions of high fluid velocities ($u > 10^{-8} \text{ m} \cdot \text{s}^{-1}$) that correspond to the geometries of tectonic inheritance present in the basement and

metasediments. Upward deflections of the isotherms indeed clearly define areas of increased fluid flow and geometries of tectonic inheritance. Additionally, fluid vectors are also observed to be oriented preferentially along these regions. At the transition between hydrostatic to lithostatic pressure regime, the temperature is at least 250°C (ca. 40°C km$^{-1}$) and up to 400°C (ca. 70°C km$^{-1}$) at the undeformed basement and the opening of the deep-seated structures, respectively.

4.2. Model M2. In model M2A, the primary convective cells occur directly above the deep-seated structures, with a wavelength of 5 km and high fluid velocities of $10^{-7} \text{ m} \cdot \text{s}^{-1}$. Temperatures at the contact between deep-seated structures
and the basin peaked at 350°C and 375°C (see Figure 7, A1 and A2, respectively). At the 6 km depth, a temperature of 375°C in the basin is near the contact of the deep-seated structure. Between the 6 and 8 km depths, regions of increased fluid velocity ($u \approx 10^{-9} \text{ m} \cdot \text{s}^{-1}$ in Figure 7(a)) are generated at the contact of the deformed basement and deep-seated structures and the impermeable basement. This depth range also corresponds to the depth of the transition from hydrostatic to lithostatic fluid pressure (in Figure 5), which results to higher fluid pressure within more permeable regions. Fluid flows within the basin layer indicate the presence of well-established convection cells (in Figures 7(a) and 7(b)). Widespread convection cells within the basin follow downward deflections of the 200°C isotherms. The wavelengths of these convection cells within the basin are smaller than the convection cells above the deep-seated structures (from 2.5 km to 1 km). Fluid temperature in these convection cells is between 200 and 225°C. Around the basin, fluid velocity is around $10^{-11} \text{ to } 10^{-10} \text{ m} \cdot \text{s}^{-1}$. When compared to M2A, fluid flow patterns, velocities, and temperature distribution in M2B are similar even though the applied dual fluid pressure regime in M2B has an open boundary (in Figure 7(b), B1 and B2). It is apparent that fluid velocity is controlled by the permeability of the rocks as the permeability profile is similar in the two models.

Regions of increased fluid pressure (i.e., 6–8 km and 12–14 km depth) are not present in model M2C (Figure 7(c)) unlike in M2A and M2B. Well-established convection cells are present at the contact of the deep-seated structures, with a wavelength of 5 km and fluid velocities of $10^{-7} \text{ m} \cdot \text{s}^{-1}$. At the contact between the deep-seated structures and the basin, the elevated temperature peaked at 300°C at 6 km depth. Fluid flows near the deep-seated structures indicate two types of flow within the basin layer (in Figure 7, C1 and C2): (1) a lateral flow at the lower parts of the basin layer and (2) a vertical fluid flow, as seen as seen from the cyclical deflection of the 200°C isotherm at the upper part of the basin layer (Figure 7, C1 and C2). The wavelength and amplitude of the 200°C isotherm deflections are small (<500 m). Furthermore, fluid flow in the basin close to the deformed basement ($x \approx -5$ to +10 km in Figure 7(c)) presents some well-established convective cells with short wavelengths (<500 m).

Loci for chlorite alteration can be found in three areas in M2A and M2B (in Figures 8(a) and 8(b)). Alteration region 1 is located at ca. 7 km depth and corresponds to high fluid velocity area within the basement. It marks the transition from hydrostatic to lithostatic pressures (in Figure 5). This alteration region corresponds to basement-hosted and not confined by the deep-seated structures. Alteration region 2 corresponds to zones associated with the deformed basement, bounded by the 300°C isotherms, and is located at 6 to 8 km depth. Here, fluids have the potential to flow from zones of lower to higher permeability [24]. The vertical geometry of the chlorite alteration aligns with the geometry of the deformed basement. The loci of the chlorite alteration is comparable with observations from Millennium deposit (e.g., Figure 3 in [9]). Alteration region 3 corresponds to the loci of chlorite alteration bounded by the 300°C isotherms within the basin layer and above the contact of the deep-seated structures. This locus of alteration region 3 is comparable to the distribution of the illite-dominated alteration region observed in Cigar Lake [43]. These results show that, at the first order, a dual fluid pressure regime can set the physical conditions allowing to create the alteration distribution for a basin-hosted unconformity-related uranium deposit type.

In M2C, the loci for chlorite alteration are restricted to alteration regions 2 and 3 (in Figure 8(c)). As it is bounded by the 200 and 300°C isotherms, alteration distribution within the deep-seated structures is limited. The loci of alteration in M2C in the basin are more extensive than M2A and M2B. The presence of well-established convective cells in our model may have sufficient fluid-rock interaction to form large alteration haloes as illite in Cigar Lake (e.g., [24, 50]).

# 5. Discussion

## 5.1. General Fluid Patterns and Convection Cells along and above Deep-Seated Structures.
In all M1 models, fluid convection cells occurred only in M1B (Figure 6(b)) in a depth-dependent permeability setting. These convection cells occur due to the presence of more impermeable zones in the metasediments at depth and the cooling effect of the top boundary (20°C). Applying less permeable zones in higher temperature regions towards the bottom boundary reduces fluid velocity, allowing thermal conduction to dominate. The resultant temperature field becomes more homogeneous when a less permeable metasediments zone in M1C is applied ($k_{\text{sedi}} = 5 \times 10^{-17} \text{ m}^2$ in Figure 6(c)). Reducing the rate of fluid-thermal advection also prevents the cooling effect of the top boundary from getting overwhelmed, causing cooler and denser fluids to flow downwards. In M1A, the absence of the downward fluid flow indicates that the cooling effect from the top boundary is insufficient and is overwhelmed by the high fluid advection.

Two main drivers of fluid flow are present in all M2 models regardless of the applied pressure regime (Figure 7): (1) thermal conduction is dominant within an impermeable basement as fluid flow is less active. This leads to unidirectional fluid flow and horizontal isotherms. (2) Thermal convective flow is dominant in permeable rock regions (i.e., the deep-seated structures and the basin). This type of flow is characterized by the deflections of the isotherms, high fluid velocities, and the convective motion of fluids. These convection cells occur in all M2 models. The maximum temperature of the fluids flowing above the deep-seated structures in M2A and M2B (ca. 350°C in Figures 7(a) and 7(b)) is higher compared to M2C (ca. 290°C in Figure 7(c)). This increase in temperature of the fluids is caused by the presence of a fluid pressure gradient from lithostatic (at depth) to hydrostatic (towards the surface) within the deep-seated structures. This fluid pressure gradient leads to an increased fluid flux and, consequently, to a heating controlled by advection.

In all M2 models (in Figure 7), widespread convective cells within the basin layer are consistent with previously published numerical results (e.g., [24, 25, 54]). The main
<figure><img src="image_9.png" /><figcaption>FIGURE 7: Thermal and fluid flow results for M2 and their variations: (a) M2A in the simplified dual-pressure regime, (b) M2B with the simplified dual-pressure regime but with an open bottom boundary, and (c) M2C in a hydrostatic pressure regime. Fluid flow vectors shown are normalised. Permeability setup for the model can be referred to Figure 3. Isotherms are represented by solid coloured lines and are in increments of 30°C. Colour fill of the plots and isotherms is based on scientific colour maps [108].</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>FIGURE 8: RRAI application to delineate loci for potential preore chlorite alteration seen with (a) M2A with the simplified dual-pressure regime, (b) M2B with the simplified dual-pressure regime and an open bottom boundary, and (c) M2C with the hydrostatic pressure regime. Partially shaded ellipses labelled 1, 2, and 3 in M2A correspond to the spatial distribution of the alteration regions and are described in detail in the text. Colour fill for the isotherms is based on scientific colour maps [108].</figcaption></figure>

difference is the scale of the numerical model. Li et al. [25, 54] considered fluid flow activity on a local deposit (i.e., the McArthur River deposit at the 5 km by 7 km depth), whereas this study considers fluid flow activity at a regional scale (i.e., WMTZ at 150 km by 15 km depth). Additionally, the deep and large structures in the M2 models are considered as large homogenous rock units within an impermeable basement (i.e., 3–5 km wide). In contrast, the numerical models conducted by Eldursi et al. [24] and Li et al. [25] considered the effects of one and multiple discrete structures

mimicking superficial fault systems at a deposit scale in a homogenous layout, respectively. Our results in the M2 models show that well-established thermal convections do occur above the deep-seated structures. At the deposit scale, our results indicate preferential fluid flow within these structures in the basin and in the deeper parts of the basement (see Figure 4 in [24]). The depth of this fluid flow within the structures also corresponds to the transition from hydrostatic to lithostatic pressure regime, favouring advective flow [90]. Our results indicate that the presence of these deep-
seated structures is essential to initiate thermal convections. Furthermore, both occurrences of long-lived convection cells above the deep-seated structures and fluid flow occurring preferentially within these structures also explain why the Wollaston-Mudjatik Transition zone is more favourable for mineral deposits.

## 5.2. Role of Inherited Deep-Seated Structures in Preparing the Geological Environment for Mineralization and Significance for Exploration.
The depth where the transition from litho-static to hydrostatic pressure regime occurs fulfils one of the conditions for the permeable structures to operate as potential fluid-pressure-activated valves [68]. These fault systems are indeed favourable loci to form alteration zones and correspond to preferential area where local deformation can occur during a transient tectonic events [14, 29]. Furthermore, the reactivation of these deep-seated structures can cause localised fluid-thermal fluxes, creating microscale networks of permeability swarms and deposit minerals (e.g., [8, 26, 91]). At a local deposit scale, hydrothermal fluids do flow out and into the structures during compres-sional and extensional strain, respectively [24, 25, 53]. The onset of fault activation can also temporary disrupt prevail-ing thermal convection within the basin [25].

Preore chlorite alteration, which acts as a chemical reduc-tant in proximity to graphite, has been identified as a poten-tial chemical trap ([28] and references therein) and as an essential vector for uranium mineralization in the Millen-nium deposit [9] and P-Patch [52] and Kiggavik-Andrew Lake structural trend prospects (e.g., [92]). The deep-seated structures are efficient conduits for transporting fluids and heat from the basement to the basin. Thus, they become the ideal loci for hosting key reductant minerals such as graphite and rehydrate metamorphic minerals such as chlorite (e.g., [13, 52]). Hydrothermal simulation results indicate that the fault systems maintain some localised fluid-thermal fluxes (i.e., natural hydrothermal state). They also show that these structures have the physical conditions for chlorite alteration to occur prior tectonic reactivation regardless of permeability configuration and pressure profile (see Figures 7 and 8, respectively).

Dynamic pressure on these structures may also play a vital role in fault reactivation. Hence, total pressure within the brittle crust is estimated to be up to ~2 times the litho-static pressure under compression and ~2/3 times the litho-static pressure under extension [93–96]. These brittle structures would become focal points for pressure variation, depending on their rheology (e.g., [97, 98]). In the Athabasca Basin, it is postulated that the impermeable seal of the fault and fluids is subjected to this pressure variation due to far-field stresses. Sources for these far-field stresses could be from the deposition of the thick (~7 km) Athabasca Basin itself (e.g., [42]), protracted deep erosion and exhumation of the THO (e.g., [41]), and the formation and breakup of supercontinents Columbia/Nuna and Rodinia (e.g., [99]). The oldest preore alteration and uranium mineralization in the Athabasca Basin are ca. 1.68 Ga and ca. 1.59 Ga, respec-tively [28] which provides about 90 Myrs for tectonic load-ing to occur. At the onset of fault formation, permeability

increases significantly and is followed by high fluid flow from depth to the basin layer. The amount of overpressured fluids ejected from fault failure can be significant, especially during the first few years (e.g., [100]). Hydrothermal precip-itation occurs if the physical and chemical conditions for preore chlorite alteration and uranium mineralization are met (Figure 8). This is also in good agreement with the fault-valve hypothesis proposed by Sibson [86, 90]. How-ever, this appears to be a transient process. After all, hydro-thermal fluids have the potential to precipitate additional minerals such as quartz, uraninite, and chlorite within the newly made permeability swarm and change the physical conditions needed for rock failure [8, 91, 96].

At a regional scale, magnetotelluric profiles have been conducted to characterize the structure of the lithosphere and structures' association to mineral deposits (e.g., [5, 6, 101]). Magnetotelluric profiles conducted towards southern Trans-Hudson Orogeny by Jones et al. [5] and in the Taltson orogeny towards the SW Athabasca Basin show low-resistivity anomalies that have a close spatial relationship with high-strain tectonic domains and to major mineral occurrences [6, 21, 22]. These observed anomalies extend over strike for hundreds of kilometres, (e.g., North American Central Plans, [5]). Moreover, field data, such as the pres-ence of carbonatites within the associated deep-seated shear zones of the Patterson Lake corridor, indicates that these deep-seated structures reach mantle depths [20]. Similar crustal-scale low-resistivity corridor networks are also well documented in the Gawler craton, Australia, and refer these high-strain zones for potential fluid flow [22, 102–104]. Sub-sequently, these tectonic structures are assumed to be strongly modified through long-standing chemical and physical processes. The rocks associated with these deep-seated inherited structures in our M2 models would have undergone multiple syntectonic depositions, tectonic reactivations, and fluid flow activity for at least 120–160 Myrs (i.e., from the peak of metamorphism of the THO to the deposition of the Athabasca Basin). Fluid flows preferentially from the impermeable basement to the permeable basin through the deep-seated structures to provide fluids and heat, potentially fulfilling critical conditions for uranium mineralization [13, 23]. Such interpretations correspond very well to the tectonic structures generated from the thermomechanical models proposed in Poh et al. [39, 65] and the structural cross-sections ([17, 35] and Figure 2(a)).

Therefore, the role of the deep-seated structures can be surmised in two main points: (1) natural candidates for fluid flow channelling and fluid-rock interaction and (2) preferen-tial zones for tectonic reactivation during transient geo-dynamic activity.

## 5.3. Refinement of Exploration Strategies.
The unconformity-related uranium deposits hosted by the Athabasca Basin are spatially associated with three main tectonic zones that still have the potential for further exploration: (1) WMTZ, (2) Snowbird tectonic zone, and (3) Patterson Lake deformation corridor. These deep-seated structures can establish convec-tion cells with a significant temperature increase (up to 350°C) and discharge a significant amount of fluid from
the basement to the basin (Figure 7). Physical conditions within the basin, deep-seated structures, and the deformed basement are indeed suitable for potential preore chlorite alteration in our models and are in good agreement with previous fluid inclusion and clay mineral studies ([9, 57, 82, 83, 105, 106]; Figure 8). The predicted alteration regions also represent a high potential for a chemical trap located within the basin layer, the deep-seated structures, and the deformed basin (Figure 8). Naturally, these structures become potential sites for fluid accumulation and high fluid-rock interaction. Within the WMTZ, the footprint of unconformity-related uranium ore bodies is spatially associated with second-order fault zones or transecting late brittle faults ([16, 17], and Figure 2(a)). It is suggested that the exploration strategy can be expanded to include deep-seated structures at depth and out of the basin, and it will potentially lead to significant discoveries.

## 6. Conclusions

This work is carried out on two tectonic scenarios mimicking the WMTZ before and after the Athabasca Basin formation and using 2-D hydrothermal models. The aim is to establish the role of deep-seated structures on fluid flow and the consequences to mineralization potential. Results show that deep-seated structures (up to 30 km depth) constitute fluid-thermal conduits capable of bringing fluids and heat towards the upper portions of the crust. Additionally, well-established and robust convective cells can occur near the deep-seated structures. Loci for preore chlorite alteration are also determined when applying the RRAI. The distribution of these areas strongly depend on the geometry and permeability of the inherited deep-seated structures. Furthermore, our results also show that these preore chlorite alteration corridors within the deep-seated structures also indicate that mechanically weak zones are prone to fault reactivation. This study may provide insight into the potential hydrodynamic activity in regions presenting a basin cover. An upscaling from camp to regional scale to properly characterize these deep-rooted, ancient structures and paleo-plate boundaries is, therefore, a viable exploration strategy.

## Data Availability

Data can be made available upon reasonable request.

## Disclosure

This work represents the third and final part of the author's Ph.D. thesis.

## Conflicts of Interest

The authors declare that they have no conflicts of interest.

## Acknowledgments

This study was supported financially by Orano Canada Inc (Project number 39M.PSGEPA0425). Numerical simulations were performed at the GeoFluids lab at the University

of Regina. The authors would like to thank Thibault Duretz, Laurent Guillou-Frottier, Remy Chemillac, and the team of Orano for their feedback and discussions on earlier versions of the manuscript.

## References

[1] T. C. McCuaig and J. M. Hronsky, “The mineral system concept: the key to exploration targeting,” *Society of Economic Geologists Special Publication*, vol. 18, pp. 153–175, 2014.

[2] L. Wyborn, C. Heinrich, and A. Jaques, “Australian Proterozoic mineral systems: essential ingredients and mappable criteria,” in *Proceedings of the Australasian Institute of Mining and Metallurgy Annual Conference*, pp. 109–115, Darwin, 1994.

[3] S. Jaireth, I. C. Roach, E. Bastrakov, and S. Liu, “Basin-related uranium mineral systems in Australia: a review of critical features,” *Ore Geology Reviews*, vol. 76, pp. 360–394, 2016.

[4] T. C. McCuaig, S. Beresford, and J. Hronsky, “Translating the mineral systems approach into an effective exploration targeting system,” *Ore Geology Reviews*, vol. 38, no. 3, pp. 128–138, 2010.

[5] A. G. Jones, J. Ledo, and I. J. Ferguson, “Electromagnetic images of the Trans-Hudson orogen: the North American Central Plains anomaly revealed,” *Canadian Journal of Earth Sciences*, vol. 42, no. 4, pp. 457–478, 2005.

[6] V. Tschirhart, S. Pehrsson, C. Card, E. G. Potter, J. Powell, and D. Pană, “Interpretation of buried basement in the southwestern Athabasca Basin, Canada, from integrated geophysical and geological datasets,” *Geochemistry: Exploration, Environment, Analysis*, vol. 21, no. 11, 2020.

[7] M. Bruce, O. Kreuzer, A. Wilde, A. Buckingham, K. Butera, and F. Bierlein, “Unconformity-type uranium systems: a comparative review and predictive modelling of critical genetic factors,” *Minerals*, vol. 10, no. 9, p. 738, 2020.

[8] M. Abdelrazek, A. Benedicto, O. Gerbeaud, and P. Ledru, “Reactivation mechanisms and related-porosity enhancement of shear zones in the context of basement-hosted uranium mineralization: case of the Spitfire discovery in the Patterson Lake corridor, Canada,” *Geochemistry: Exploration, Environment, Analysis*, vol. 21, no. 3, 2021.

[9] J. Cloutier, K. Kyser, G. R. Olivo, P. Alexandre, and J. Halaburda, “The Millennium uranium deposit, Athabasca Basin, Saskatchewan, Canada: an atypical basement-hosted unconformity-related uranium deposit,” *Economic Geology*, vol. 104, no. 6, pp. 815–840, 2009.

[10] D. D. Johnstone, K. M. Bethune, C. D. Card, and V. Tschirhart, “Structural evolution and related implications for uranium mineralization in the Patterson Lake corridor, southwestern Athabasca Basin, Saskatchewan, Canada,” *Geochemistry: Exploration, Environment, Analysis*, vol. 21, no. 1, 2021.

[11] R. Ng, P. Alexandre, and K. Kyser, “Mineralogical and geochemical evolution of the unconformity-related McArthur River Zone 4 orebody in the Athabasca Basin, Canada: implications of a silicified zone,” *Economic Geology*, vol. 108, no. 7, pp. 1657–1689, 2013.

[12] M. Rabiei, G. Chi, E. G. Potter et al., “Fluid evolution along the Patterson Lake corridor in the southwestern Athabasca Basin: constraints from fluid inclusions and implications for unconformity-related uranium mineralization,” *Geochemistry: Exploration, Environment, Analysis*, vol. 21, no. 3, 2021.
[13] C. W. Jefferson, D. J. Thomas, S. S. Gandhi et al., "Unconformity-associated uranium deposits of the Athabasca Basin, Saskatchewan and Alberta," *Bulletin-Geological Survey of Canada*, vol. 588, p. 23, 2007.

[14] I. R. Annesley, C. Madore, and P. Portella, "Geology and thermotectonic evolution of the western margin of the Trans-Hudson Orogen: evidence from the eastern sub-Athabasca basement, Saskatchewan," *Canadian Journal of Earth Sciences*, vol. 42, no. 4, pp. 573-597, 2005.

[15] P. Jeanneret, P. Goncalves, C. Durand et al., "Tectono-metamorphic evolution of the pre-Athabasca basement within the Wollaston-Mudjatik Transition Zone, Saskatchewan," *Canadian Journal of Earth Sciences*, vol. 53, pp. 231-259, 2016.

[16] P. Ledru, A. Benedicto, G. Chi et al., "The unconformity-related uranium mineral system of the Athabasca Basin (Canada)," in *Dans Ressources métalliques 2*, S. Decrée, Ed., ISTE Editions, Londres, 2023.

[17] P. Ledru, A. Benedicto, J. Poh, and J. Robbins, "Implications of a new tectonic model on the role of pre-Athabasca basement-hosted fault systems in the formation of the unconformity related uranium deposits," in *Open House 2019 Abstract Volume, Saskatchewan Ministry of Energy and Resources*, p. 4, Saskatchewan Geological Survey, 2019.

[18] K. Witherly, "Geophysical expressions of ore systems—our current understanding," in *Building Exploration Capability for the 21st Century*, K. D. Kelley and H. C. Golden, Eds., pp. 177-208, Society of Economic Geologists, 2014.

[19] S. Hillacre, K. Ansdell, and B. McEwan, "Geology, Structural Analysis, and Paragenesis of the Arrow Uranium Deposit, Western Athabasca Basin, Saskatchewan, Canada: Implications for the Development of the Patterson Lake Corridor," *Economic Geology*, vol. 116, no. 2, pp. 285-321, 2021.

[20] C. D. Card, "The Patterson Lake corridor of Saskatchewan, Canada: defining crystalline rocks in a deep-seated structure that hosts a giant, high-grade Proterozoic unconformity uranium system," *Geochemistry: Exploration, Environment, Analysis*, vol. 21, no. 1, 2021.

[21] E. G. Potter, V. Tschirhart, J. W. Powell et al., "Targeted geoscience initiative 5: integrated multidisciplinary studies of unconformity-related uranium deposits from the Patterson Lake corridor, northern Saskatchewan," Technical Report, 2020.

[22] V. Tschirhart, E. G. Potter, J. W. Powell, E. A. Roots, and J. A. Craven, "Deep geological controls on formation of the highest-grade uranium deposits in the world: magnetotelluric imaging of unconformity-related systems from the Athabasca Basin, Canada," *Geophysical Research Letters*, vol. 49, no. 15, 2022.

[23] G. Chi, H. Chu, D. Petts, E. Potter, S. Jackson, and A. Williams-Jones, "Uranium-rich diagenetic fluids provide the key to unconformity-related uranium mineralization in the Athabasca Basin," *Scientific Reports*, vol. 9, no. 1, p. 5530, 2019.

[24] K. Eldursi, G. Chi, K. Bethune, Z. Li, P. Ledru, and D. Quirt, "New insights from 2- and 3-d numerical modelling on fluid flow mechanisms and geological factors responsible for the formation of the world-class cigar lake uranium deposit, eastern Athabasca basin, Canada," *Mineralium Deposita*, vol. 56, no. 7, pp. 1365-1388, 2021.

[25] Z. Li, G. Chi, K. M. Bethune et al., "Interplay between thermal convection and compressional fault reactivation in the formation of unconformity-related uranium deposits," *Mineralium Deposita*, vol. 56, no. 7, pp. 1389-1404, 2021.

[26] P. Martz, M. Cathelineau, J. Mercadier et al., "C-O-H-N fluids circulations and graphite precipitation in reactivated Hudsonian shear zones during basement uplift of the Wollaston-Mudjatik Transition Zone: example of the Cigar Lake U deposit," *Lithos*, vol. 294-295, pp. 222-245, 2017.

[27] K. Eldursi, Y. Branquet, L. Guillou-Frottier, and E. Marcoux, "Numerical investigation of transient hydrothermal processes around intrusions: heat-transfer and fluid-circulation controlled mineralization patterns," *Earth and Planetary Science Letters*, vol. 288, no. 1-2, pp. 70-83, 2009.

[28] P. Alexandre, K. Kyser, D. Thomas, P. Polito, and J. Marlat, "Geochronology of unconformity-related uranium deposits in the Athabasca Basin, Saskatchewan, Canada and their integration in the evolution of the basin," *Mineralium Deposita*, vol. 44, no. 1, pp. 41-59, 2009.

[29] P. Audet and R. Bürgmann, "Dominant role of tectonic inheritance in supercontinent cycles," *Nature Geoscience*, vol. 4, no. 3, pp. 184-187, 2011.

[30] C. D. Card, K. M. Bethune, W. J. Davis, N. Rayner, and K. E. Ashton, "The case for a distinct Taltson orogeny: evidence from northwest Saskatchewan, Canada," *Precambrian Research*, vol. 255, pp. 245-265, 2014.

[31] D. Corrigan, Z. Hajnal, B. Németh, and S. B. Lucas, "Tectonic framework of a Paleoproterozoic arc-continent to continent-continent collisional zone, Trans-Hudson Orogen, from geological and seismic reflection studies," *Canadian Journal of Earth Sciences*, vol. 42, no. 4, pp. 421-434, 2005.

[32] D. Corrigan, S. Pehrsson, N. Wodicka, and E. de Kemp, "The Palaeoproterozoic Trans-Hudson Orogen: a prototype of modern accretionary processes," *Geological Society, London, Special Publications*, vol. 327, no. 1, pp. 457-479, 2009.

[33] P. F. Hoffman, "United plates of America, the birth of a craton: Early Proterozoic assembly and growth of Laurentia," *Annual Review of Earth and Planetary Sciences*, vol. 16, no. 1, pp. 543-603, 1988.

[34] S. J. Whitmeyer and K. E. Karlstrom, "Tectonic model for the Proterozoic growth of North America," *Geosphere*, vol. 3, no. 4, pp. 220-259, 2007.

[35] A. Benedicto, M. Abdelrazek, P. Ledru, C. MacKay, and D. Kinar, "Structural controls of uranium mineralization in the basement of the Athabasca Basin, Saskatchewan, Canada," *Geofluids*, vol. 2021, Article ID 3853468, 30 pages, 2021.

[36] D. Gapais, "Tectonics-mineralisation relationships within weak continental lithospheres: a new structural framework for Precambrian cratons," *Bulletin de la Société Géologique de France*, vol. 189, no. 3, p. 14, 2018.

[37] H. Tran, G. Yeo, and K. Bethune, "Geology of the McKenzie Falls area, Haultain River, Wollaston-Mudjatik domains boundary (NTS 74B-7 and-8)," *Summary of Investigations*, vol. 2, p. 99-4, 1999.

[38] P. Jeanneret, P. Goncalves, C. Durand et al., "Geochronological constraints on the Trans-Hudsonian tectono-metamorphic evolution of the pre-Athabasca basement within the Wollaston-Mudjatik Transition Zone, Saskatchewan," *Precambrian Research*, vol. 301, pp. 152-178, 2017.

[39] J. Poh, P. Yamato, T. Duretz, D. Gapais, and P. Ledru, "Precambrian deformation belts in compressive tectonic regimes: a numerical perspective," *Tectonophysics*, vol. 777, article 228350, 2020.
[40] E. Adlakha and K. Hattori, "Thermotectonic events recorded by U-Pb geochronology and Zr-in-rutile thermometry of Ti oxides in basement rocks along the P2 fault, eastern Athabasca Basin, Saskatchewan, Canada," *GSA Bulletin*, vol. 134, no. 3-4, pp. 567-576, 2021.

[41] J. Chakungal, P. H. Reynolds, R. A. Jamieson, and D. Corrigan, "Slow post-orogenic cooling in a deeply eroded orogen: Reindeer Zone, Trans-Hudson Orogen, Saskatchewan," *Canadian Journal of Earth Sciences*, vol. 41, no. 7, pp. 867-880, 2004.

[42] P. Ramaekers, C. Jefferson, G. Yeo et al. et al., "Revised geological map and stratigraphy of the Athabasca Group, Saskatchewan and Alberta," *Bulletin-Geological Survey of Canada*, vol. 588, p. 155, 2007.

[43] C. S. Bishop, A. G. Mainville, and L. D. Yesnik, *Cigar Lake Operation Northern Saskatchewan, Canada*, Technical Report. Cameco Corporation, 2016.

[44] C. D. Card, D. Pan, R. A. Stern, and N. Rayner, "New insights into the geological history of the basement rocks to the southwestern Athabasca Basin, Saskatchewan and Alberta," *EXTECH IV: Geology and Uranium Exploration Technology of the Proterozoic Athabasca Basin, Saskatchewan and Alberta*, C. W. Jefferson and G. D. Delaney, Eds., pp. 119-134, 2007.

[45] C. D. Card, D. Pana, P. Portella, D. J. Thomas, and I. R. Annesley, "Basement rocks to the Athabasca Basin, Saskatchewan and Alberta," *EXTECH IV: Geology and Uranium Exploration Technology of the Proterozoic Athabasca Basin, Saskatchewan and Alberta*, C. W. Jefferson and G. D. Delaney, Eds., pp. 69-87, 2007.

[46] B. McGill, J. Marlatt, R. Matthews, V. Sopuck, L. Homeniuk, and J. Hubregtse, "The P2 North uranium deposit, Saskatchewan, Canada," *Exploration and Mining Geology*, vol. 2, pp. 321-331, 1993.

[47] G. Tourigny, D. H. Quirt, N. S. F. Wilson, S. Wilson, G. Breton, and P. Portella, "Geological and structural features of the Sue C uranium deposit," *EXTECH IV: Geology and Uranium Exploration Technology of the Proterozoic Athabasca Basin, Saskatchewan and Alberta*, C. W. Jefferson and G. D. Delaney, Eds., pp. 229-247, 2007.

[48] Z. Li, K. M. Bethune, G. Chi, S. A. Bosman, C. D. Card, and J. Hanley, "Topographic features of the sub-Athabasca Group unconformity surface in the southeastern Athabasca Basin and their relationship to uranium ore deposits," *Canadian Journal of Earth Sciences*, vol. 52, no. 10, pp. 903-920, 2015.

[49] J. Hoeve and D. Quirt, "A stationary redox front as a critical factor in the formation of high-grade, unconformity-type uranium ores in the Athabasca basin, Saskatchewan, Canada," *Bulletin de Minéralogie*, vol. 110, no. 2, pp. 157-171, 1987.

[50] P. Martz, J. Mercadier, M. Cathelineau et al., "Formation of U-rich mineralizing fluids through basinal brine migration within basement-hosted shear zones: a large-scale study of the fluid chemistry around the unconformity-related Cigar Lake U deposit (Saskatchewan, Canada)," *Chemical Geology*, vol. 508, pp. 116-143, 2019.

[51] J. Mercadier, I. Annesley, C. McKechnie, T. Bogdan, and S. Creighton, "Magmatic and metamorphic uraninite mineralization in the western margin of the Trans-Hudson Orogen (Saskatchewan, Canada): a uranium source for unconformity-related uranium deposits?," *Economic Geology*, vol. 108, no. 5, pp. 1037-1065, 2013.

[52] J. Mercadier, A. Richard, M. C. Boiron, M. Cathelineau, and M. Cuney, "Migration of brines in the basement rocks of the Athabasca Basin through microfracture networks (P-Patch U deposit, Canada)," *Lithos*, vol. 115, no. 1-4, pp. 121-136, 2010.

[53] T. Cui, J. Yang, and I. M. Samson, "Tectonic deformation and fluid flow: implications for the formation of unconformity-related uranium deposits," *Economic Geology*, vol. 107, no. 1, pp. 147-163, 2012.

[54] Z. Li, G. Chi, K. M. Bethune et al., "Synchronous egress and ingress fluid flow related to compressional reactivation of basement faults: the Phoenix and Gryphon uranium deposits, southeastern Athabasca Basin, Saskatchewan, Canada," *Mineralium Deposita*, vol. 53, pp. 277-292, 2017.

[55] Z. Li, G. Chi, and K. Bethune, "The effects of basement faults on thermal convection and implications for the formation of unconformity-related uranium deposits in the Athabasca Basin, Canada," *Geofluids*, vol. 16, no. 4, pp. 729-751, 2016.

[56] M. Fayek and T. K. Kyser, "Characterisation of multiple fluid-flow events and rare-earth-element mobility associated with formation of unconformity-type uranium deposits in the Athabasca Basin, Saskatchewan," *The Canadian Mineralogist*, vol. 3, pp. 627-658, 1997.

[57] T. G. Kotzer and T. K. Kyser, "Petrogenesis of the Proterozoic Athabasca Basin, northern Saskatchewan, Canada, and its relation to diagenesis, hydrothermal uranium mineralization and paleohydrogeology," *Chemical Geology*, vol. 120, no. 1-2, pp. 45-89, 1995.

[58] D. H. Quirt, "Athabasca unconformity-type uranium deposits: one deposit type with many variations," in *Uranium Geochemistry 2003 Proceedings*, pp. 309-312, Nancy, France, 2003.

[59] Itasca, *FLAC3D: Fast Lagrangian Analysis of Continua in 3 Dimensions*, User Manual, Version 5.0, Itasca Consulting Group Inc., Minneapolis, Minnesota, 2012.

[60] N. H. Oliver, J. G. McLellan, B. E. Hobbs, J. S. Cleverley, A. Ord, and L. Feltrin, "100th anniversary special paper: numerical models of extensional deformation, heat transfer, and fluid flow across basement-cover interfaces during basin-related mineralization," *Economic Geology*, vol. 101, no. 1, pp. 1-31, 2006.

[61] H. A. Sheldon, "Simulation of magmatic and metamorphic fluid production coupled with deformation, fluid flow and heat transport," *Computers & Geosciences*, vol. 35, no. 11, pp. 2275-2281, 2009.

[62] P. Upton and D. Craw, "Extension and gold mineralisation in the hanging walls of active convergent continental shear zones," *Journal of Structural Geology*, vol. 64, pp. 135-148, 2014.

[63] C. Zhao, B. E. Hobbs, and H. B. Mühlhaus, "Finite element modelling of temperature gradient driven rock alteration and mineralization in porous rock masses," *Computer Methods in Applied Mechanics and Engineering*, vol. 165, no. 1-4, pp. 175-187, 1998.

[64] D. Gapais, J. Jaguin, F. Cagnard, and P. Boulvais, "Pop-down tectonics, fluid channelling and ore deposits within ancient hot orogens," *Tectonophysics*, vol. 618, pp. 102-106, 2014.

[65] J. Poh, P. Yamato, T. Duretz, D. Gapais, and P. Ledru, "The transition from ancient to modern-style tectonics: insights from lithosphere dynamics modelling in compressional regimes," *Gondwana Research*, vol. 99, pp. 77-92, 2021.
[66] Z. Hajnal, J. Lewry, D. White et al., "The Sask Craton and Hearne Province margin: seismic reflection studies in the western Trans-Hudson Orogen," *Canadian Journal of Earth Sciences*, vol. 42, no. 4, pp. 403-419, 2005.

[67] J. F. Lewry and T. I. Sibbald, "Thermotectonic evolution of the Churchill province in northern Saskatchewan," *Tectono-physics*, vol. 68, no. 1-2, pp. 45-82, 1980.

[68] R. H. Sibson, "Crustal stress, faulting and fluid flow," *Geological Society, London, Special Publications*, vol. 78, no. 1, pp. 69-84, 1994.

[69] H. A. Sheldon and A. Ord, "Evolution of porosity, permeability and fluid pressure in dilatant faults post-failure: implications for fluid flow and mineralization," *Geofluids*, vol. 5, no. 4, pp. 272-288, 2005.

[70] G. Chi, S. Bosman, and C. Card, "Numerical modeling of fluid pressure regime in the Athabasca basin and implications for fluid flow models related to the unconformity-type uranium mineralization," *Journal of Geochemical Exploration, vol. 125*, pp. 8-19, 2013.

[71] M. Pagel, "Détermination des conditions physico-chimiques de la silicification diagénétique des grès Athabasca (Canada) au moyen des inclusions fluids," *Comptes Rendus Série D*, vol. 280, pp. 2301-2304, 1975.

[72] R. A. Freeze and J. A. Cherry, *Groundwater*, Prentice-Hall, Englewoods Cliffs, New Jersey, 1979.

[73] J. Rutqvist and O. Stephansson, "The role of hydromechanical coupling in fractured rock engineering," *Hydrogeology Journal*, vol. 11, no. 1, pp. 7-40, 2003.

[74] V. Harcouët-Menou, L. Guillou-Frottier, A. Bonneville, P. M. Adler, and V. Mourzenko, "Hydrothermal convection in and around mineralized fault zones: insights from two- and three-dimensional numerical modeling applied to the Ashanti belt, Ghana," *Geofluids*, vol. 9, no. 2, pp. 116-137, 2009.

[75] S. Ingebritsen and C. E. Manning, "Geological implications of a permeability-depth curve for the continental crust," *Geology*, vol. 27, no. 12, pp. 1107-1110, 1999.

[76] S. E. Ingebritsen and C. E. Manning, "Permeability of the continental crust: dynamic variations inferred from seismicity and metamorphism," *Geofluids*, vol. 10, no. 1-2, pp. 193-205, 2010.

[77] V. M. Shmonov, V. M. Vitiovtova, A. V. Zharikov, and A. A. Grafchikov, "Permeability of the continental crust: implications of experimental data," *Journal of Geochemical Exploration*, vol. 78-79, pp. 697-699, 2003.

[78] H. K. H. Olierook, N. E. Timms, and P. J. Hamilton, "Mechanisms for permeability modification in the damage zone of a normal fault, northern Perth Basin, Western Australia," *Marine and Petroleum Geology*, vol. 50, pp. 130-147, 2014.

[79] J. Scibek, "Multidisciplinary database of permeability of fault zones and surrounding protolith rocks at world-wide sites," *Scientific Data*, vol. 7, no. 1, p. 95, 2020.

[80] M. C. Boiron, M. Cathelineau, and A. Richard, "Fluid flows and metal deposition near basement/cover unconformity: lessons and analogies from Pb-Zn-F-Ba systems for the understanding of Proterozoic U deposits," *Geofluids*, vol. 10, no. 1-2, pp. 270-292, 2010.

[81] C. Clauser and E. Huenges, "Thermal conductivity of rocks and minerals," *Rock Physics & Phase Relations: a Handbook of Physical Constants*, vol. 3, pp. 105-126, 1995.

[82] J. Hoeve and D. Quirt, *Mineralisation and host rock alteration in relation to clay mineral diagenesis and evolution of the Middle-Proterozoic, Athabasca Basin, northern Saskatchewan, Canada*, Technical Report, Saskatchewan Research Council, 1984.

[83] J. Hoeve, D. Quirt, and D. Alonso, "Clay mineral stratigraphy of the Athabasca Group: correlation inside and outside the Carswell Structure uranium deposit, Saskatchewan," *Geological Association of Canada, Special Paper*, vol. 29, p. 12, 1985.

[84] L. Cathles and J. J. Adams, "Fluid flow and petroleum and mineral resources in the upper (<20-km) continental crust," *Economic Geology* 100th Anniversary, vol. 77-110, 2005.

[85] S. F. Cox, "Coupling between deformation, fluid pressures, and fluid flow in ore-producing hydrothermal systems at depth in the crust," in *One Hundredth Anniversary Volume. Society of Economic Geologists*, pp. 39-75, 2005.

[86] R. H. Sibson, "Preparation zones for large crustal earthquakes consequent on fault-valve action," *Earth, Planets and Space*, vol. 72, no. 1, 2020.

[87] O. M. Phillips, *Geological Fluid Dynamics: Sub-Surface Flow and Reactions*, Cambridge University Press, 2009.

[88] B. E. Hobbs, Y. Zhang, A. Ord, and C. Zhao, "Application of coupled deformation, fluid flow, thermal and chemical modelling to predictive mineral exploration," *Journal of Geochemical Exploration*, vol. 69-70, pp. 505-509, 2000.

[89] P. Alexandre, K. Kyser, P. Polito, and D. Thomas, "Alteration mineralogy and stable isotope geochemistry of Paleoproterozoic basement-hosted unconformity-type uranium deposits in the Athabasca Basin, Canada," *Economic Geology*, vol. 100, no. 8, pp. 1547-1563, 2005.

[90] R. H. Sibson, "Dual-driven fault failure in the lower seismic zone," *Bulletin of the Seismological Society of America*, vol. 110, no. 2, pp. 850-862, 2020.

[91] S. F. Cox, "Injection-driven swarm seismicity and permeability enhancement: Implications for the dynamics of hydrothermal ore systems in high fluid-flux, overpressured faulting regimes - An invited paper," *Economic Geology*, vol. 111, no. 3, pp. 559-587, 2016.

[92] T. Riegler, J. L. Lescuyer, P. Wollenberg, D. Quirt, and D. Beaufort, "Alteration related to uranium deposits in the Kiggavik-Andrew Lake structural trend, Nunavut, Canada: new insights from petrography and clay mineralogy," *The Canadian Mineralogist*, vol. 52, pp. 27-45, 2014.

[93] T. Gerya, "Tectonic overpressure and underpressure in lithospheric tectonics and metamorphism," *Journal of Metamorphic Geology*, vol. 33, no. 8, pp. 785-800, 2015.

[94] K. Petrini and Y. Podladchikov, "Lithospheric pressure-depth relationship in compressive regions of thickened crust," *Journal of Metamorphic Geology*, vol. 18, no. 1, pp. 67-77, 2000.

[95] P. Yamato and J. P. Brun, "Metamorphic record of catastrophic pressure drops in subduction zones," *Nature Geoscience*, vol. 10, pp. 46-50, 2017.

[96] A. Bauville and P. Yamato, "Pressure-to-depth conversion models for metamorphic rocks: derivation and applications," *Geochemistry, Geophysics, Geosystems*, vol. 22, no. 1, article e2020GC009280, 2021.

[97] E. Moulas, J. P. Burg, and Y. Podladchikov, "Stress field associated with elliptical inclusions in a deforming matrix: mathematical model and implications for tectonic overpressure in the lithosphere," *Tectonophysics*, vol. 631, pp. 37-49, 2014.

[98] S. M. Schmalholz and Y. Y. Podladchikov, "Tectonic overpressure in weak crustal-scale shear zones and implications
for the exhumation of high-pressure rocks,” *Geophysical Research Letters*, vol. 40, no. 10, pp. 1984–1988, 2013.

[99] R. M. Palin, M. Santosh, W. Cao, S.-S. Li, D. Hernández-Uribe, and A. Parsons, “Secular change and the onset of plate tectonics on Earth,” *Earth-Science Reviews*, vol. 207, article 103172, 2020.

[100] S. Micklethwaite, A. Ford, W. Witt, and H. A. Sheldon, “The where and how of faults, fluids and permeability – insights from fault stepovers, scaling properties and gold mineralisation,” *Geofluids*, vol. 15, no. 1-2, pp. 240–251, 2015.

[101] A. Jones, C. Lawley, D. Snyder, J. C. Afonso, H. Grutter, and G. Begg, *Geophysical and geochemical imaging of the continental upper mantle for mineral systems regional targeting (short course)*, Prospectors & Developers Association of Canada, 2019.

[102] G. Heinson, Y. Didana, P. Soeffky, S. Thiel, and T. Wise, “The crustal geophysical signature of a world-class magmatic mineral system,” *Scientific Reports*, vol. 8, no. 1, p. 10608, 2018.

[103] J. G. Motta, P. G. Betts, C. R. de Souza Filho, S. Thiel, S. Curtis, and R. J. Armit, “Proxies for basement structure and its implications for Mesoproterozoic metallogenic provinces in the Gawler Craton,” *Journal of Geophysical Research* - *Solid Earth*, vol. 124, no. 3, pp. 3088–3104, 2019.

[104] S. Thiel and G. Heinson, “Crustal imaging of a mobile belt using magnetotellurics: an example of the Fowler Domain in South Australia,” *Journal of Geophysical Research*, vol. 115, no. B6, p. B06102, 2010.

[105] D. Derome, M. Cathelineau, M. Cuney, C. Fabre, T. Lhomme, and D. A. Banks, “Mixing of sodic and calcic brines and uranium deposition at McArthur River, Saskatchewan, Canada: a Raman and laser-induced breakdown spectroscopic study of fluid inclusions,” *Economic Geology*, vol. 100, no. 8, pp. 1529–1545, 2005.

[106] A. Richard, M. Cathelineau, M. C. Boiron, J. Mercadier, D. A. Banks, and M. Cuney, “Metal-rich fluid inclusions provide new insights into unconformity-related U deposits (Athabasca Basin and basement, Canada),” *Mineralium Deposita*, vol. 51, pp. 249–270, 2015.

[107] B. M. Eglington, S. J. Pehrsson, K. M. Ansdell et al., “A domain-based digital summary of the evolution of the Palaeoproterozoic of North America and Greenland and associated unconformity-related uranium mineralization,” *Precambrian Research*, vol. 232, pp. 4–26, 2013.

[108] F. Crameri, *Scientific colour maps*, 2021.