Forward Modeling and 3D Inversion of EM Data Collected Over the McArthur River Uranium Deposit in the Athabasca Basin, Canada

Detection and assessment of the deeply buried high-grade uranium deposits in the Athabasca Basin rely on geophysical methods to map conductive rocks. Variable Quaternary surface cover can mask the anomalous signals from depth and affect interpretation of inverted conductivity models. We present the analysis of a number of EM modeling studies and two field data sets, to demonstrate the effects of varying Quaternary cover resistivity and thickness, on the ability to resolve the parameters of underlying sandstone, alteration, and basement conductors. Synthetic data, assuming a typical shallow EM sounding system and realistic resistivities found in the Athabasca Basin, show that resistivity and thickness parameters of the Quaternary cover can be separately recovered in cases where this cover is more conductive than the underlying sandstone, but not when the cover is significantly more resistive. A 3D modeling study shows that using airborne EM data, it is possible to detect a basement conductor of 20 S at a depth of at least 600 m below surface, even in the presence of Quaternary cover thickness variations of the up to 20% (40 m to 60 m). Furthermore, while Quaternary cover variations and deeper sandstone alteration can produce comparable anomalous signal amplitudes in a time-domain EM response, their effects are most visible in distinctly separate time windows. Analysis of a GPR field data set to map the thickness of Quaternary cover in the McArthur River area, indicates that this cover consists mostly of sandy tills and ranges in thickness from 0 to 117 m. Constrained 3D inversion of an airborne EM data set from the same area shows basement conductors consistent with the depth and location of a known fault. Elevated conductivity in the sandstone by up to a factor of two over the background values could indicate possible alteration.

Mineralogy and petrogenesis of fracture coatings in Athabasca Group sandstones from the McArthur River uranium deposit

ABSTRACT The McArthur River unconformity-related uranium deposit, located in the Athabasca Basin of Saskatchewan, Canada, is structurally hosted near the unconformity between Archean to Paleoproterozoic metasedimentary basement and the Proterozoic Athabasca Group sandstones. In this study, the mineralogy and geochemistry of fracture materials within the entire ca. 550 m thickness of the Athabasca Group sandstones and the metasedimentary (host) rocks from the McArthur River area were used to determine the paragenetic sequence and origin of minerals in and near the fractures. Our work sought to determine if the host minerals record elements associated with the uranium deposit at depth and if they could be used to guide exploration (vectoring). Fracture orientations indicate that most are moderately dipping (<50°) and provided permeable pathways for fluid movement within the basin, from below, and through the overlying sedimentary rocks. Many of the fractures and adjacent wall rocks record evidence of multiple distinct fluid events. Seven types of fracture fillings were identified from drill core intersecting the Athabasca Basin and present distinct colors, mineralogy, and chemical features. Brown (Type 1) and pink (Type 7) fractures host paragenetically late botryoidal goethite, Mn oxide minerals, and poorly crystallized kaolinite that formed from relatively recent low-temperature meteoric fluids, as indicated by poor crystallinity and low δ2H values of –198 to –115‰. These minerals variably replaced higher temperature minerals that are rarely preserved on the fractures or in wall rock near the fractures. Hydrothermal alteration associated with the mineralizing system at ca. 200 °C is recorded in assemblages of dickite, well-crystallized kaolinite, and spherulitic dravite in some white and yellow (Type 2) and white (Type 3) fractures, as reflected by the crystal habits and variable δ2H values of –85 to –44‰. Fibrous goethite in white and yellow (Type 2) and black and orange (Type 5) fractures and microfibrous Mn oxy-hydroxide minerals in black (Type 4) fractures also crystallized from hydrothermal fluids, but at temperatures less than 200 °C. White and yellow fractures (Type 2) containing fibrous goethite reflect fracture networks indicative of hydrothermal fluids associated with the mineralizing system during primary dispersion of pathfinder elements and therefore extend the deposit footprint. Brown (Type 1) and pink (Type 7) fractures have low δ2H values in botryoidal goethite and poorly crystallized kaolinite and are indicative of the movement of meteoric waters. Secondary dispersion of elements from the deposit to the surface on some fractures is evidence that fractures are pathways for element migration from the deposit to the surface, over distances exceeding ∼500 m.

Geometallurgical ore characterization of the high-grade polymetallic unconformity-related uranium deposit

ABSTRACT Cigar Lake is a polymetallic, unconformity-related uranium deposit with complex geochemistry and mineralogy located in the eastern Athabasca Basin of northern Saskatchewan, Canada. Variable concentrations and spatial distributions of elements of concern, such as As, Mo, Ni, Co, Se, and Zr, associated with the high-grade tetravalent uranium ores [UO2+x; U(SiO4)1–x(OH)4x] present unique mining, metallurgical, and environmental challenges. Sulfide and arsenide minerals have significant control over As, Mo, Ni, Co, and Se abundances and have properties that affect element of concern mobility, thus requiring consideration during mineral processing, mine-effluent water treatment, and long-term tailings management. The U-bearing (uraninite, coffinite) and metallic arsenide (nickeline, often called “niccolite” in the past), sulfarsenide (gersdorffite, cobaltite), and sulfide (chalcopyrite, pyrite, galena, bornite, chalcocite, sphalerite, pyrrhotite) minerals provide the main controls on the distributions of the elements of concern. Arsenic, Ni, and Co occur primarily in a reduced state as 1:1 molar ratio, Ni-Co:As, arsenide, and sulfarsenide minerals such as gersdorffite, nickeline, and cobaltite. Molybdenum occurs within molybdenite and uraninite. Selenium occurs within coffinite, sulfide, and sulfarsenide minerals. Zirconium is found within detrital zircon and coffinite. The spatial distribution and paragenesis of U-, As-, and S-bearing minerals are a result of the elemental composition, pH, and redox conditions of early formational and later meteoric fluids that formed and have modified the deposit through access along lithostratigraphic permeability and tectonic structures. Using the holistic geometallurgical paradigm presented here, the geochemistry and mineral chemistry at Cigar Lake can be used to optimize and reduce risk during long-term mine and mill planning.

Reactivation mechanisms and related-porosity enhancement of shear zones in context of basement-hosted uranium mineralization: case of the Spitfire discovery in the Patterson Lake corridor, Canada

Uranium mineralization in the Patterson Lake corridor (southwestern Churchill province, Canada) is hosted in the metamorphosed Paleoproterozoic basement covered to the North by the flat-lying sandstone formations of the Athabasca Basin. The mineralization is exclusively contained within inherited ductile structures that were reactivated under a brittle regime. Petrographic and micro-structural studies of drill core samples from the Spitfire discovery (Hook Lake Project) reveal the linkages between structural evolution of the basement, alteration and mineralization. During basement exhumation, localization of non-coaxial deformation led to the formation of a large anastomosing shear zone system made of mylonitic rocks. Strain localization associated with fluid circulation induced strong mineralogical and rheological changes, forming discontinuities in mechanical anisotropy. During and post-deposition of the Athabasca Basin after 1.80 Ga, these zones of anisotropy localized brittle reactivation, expressed by a network of micro-fractures later amplified by dissolution processes which enhanced porosity later filled with phyllosilicates and uranium oxides. Crosscutting relationships between alteration minerals and structures indicate that fluid circulation was active after the basement exhumation. Uranium-bearing fluids moved through the network of micro-fractures. As shown for the Spitfire prospect, fertile structures in the basement below the Athabasca Basin have a combined poly-phase structural and alteration history during which development of ductile shear zones followed by brittle reactivation and dissolution processes led to the formation of superimposed shear and damaged zones in which uranium orebodies are located.Thematic collection: This article is part of the Uranium Fluid Pathways collection available at: https://www.lyellcollection.org/cc/uranium-fluid-pathways

3D electromagnetic modeling of graphitic faults in the Athabasca Basin using a finite-volume time-domain approach with unstructured grids

Uranium exploration in the Athabasca Basin, Canada, relies heavily on ground-based transient electromagnetic (TEM) surveys to target thin, steeply dipping graphitic conductors that are often closely related to the uranium ore deposits. The interpretation of TEM data is important in identifying the locations and trends of conductors in order to guide subsequent drilling campaigns. We present a trial-and-error modeling approach and its application to the interpretation of a data set acquired at Close Lake in the Athabasca Basin. The modeling process has two key tasks: building geo-electric models and computing their TEM responses. The modeling process is repeated with the geo-electric model being iteratively refined based on the match between three-component calculated and measured data from early to late times. To create geo-electric models, we first build a realistic geological model and discretize it using an unstructured tetrahedral mesh, with each mesh cell populated with appropriate resistivities. To calculate the TEM responses of the geo-electric model, we use a 3D finite-volume time-domain (FVTD) algorithm. We construct our initial model based on existing geologic information and drilling data. We show that this modeling process is flexible and can easily handle thin, steeply dipping conductive graphitic fault models with variable resistivities in the fault and background, and with topography. Our interpretation of the Close Lake data matches well with the trend and location of the main conductor as revealed by drilling data, and also confirms the existence of a smaller conductor which only caused noticeable anomalous responses in early-time horizontal-component data. The smaller conductor was suggested by previous electromagnetic data but was missed in a recent interpretation based on the modeling of only late-time vertical component data with plate-based approximate modeling methods.

Fluid evolution along the Patterson Lake corridor in the southwestern Athabasca Basin: constraints from fluid inclusions and implications for unconformity-related uranium mineralization

The Patterson Lake corridor (PLC) in the southwestern margin of the Athabasca Basin hosts several high-grade uranium deposits. These deposits are located in the basement up to 900 m below the unconformity surface, raising questions about their affiliation with typical unconformity-related uranium (URU) deposits elsewhere in the basin. Based on cross-cutting relationships four pre- and three syn- to post-mineralization quartz generations were identified. Fluid inclusion analyses indicate that pre-mineralization fluids have salinities ranging from 0.2 to 27.2 Wt% NaCl equiv. (avg. 9.0 Wt%), whereas syn-mineralization fluids have salinities ranging from 8.8 to 33.8 Wt% NaCl + CaCl2 (avg. 25.4 Wt%), with NaCl- and CaCl2-rich varieties. The homogenization temperatures (Th) of fluid inclusions from pre-mineralization quartz range from 80 ° to 244 ℃ (avg. 147 ℃), and from syn-mineralization quartz range from 64 ° to 248 ℃ (avg. 128 ℃). Fluid boiling is indicated by the co-development of liquid-dominated and vapor-dominated fluid inclusions within individual fluid inclusion assemblages (FIA) from the syn-mineralization quartz and is related to episodic fluid pressure drops caused by reactivation of basement faults. Our results indicate that composition and P-T conditions of the ore fluids in the PLC are comparable to those of typical URU deposits in the Athabasca Basin, indicating that the uranium deposits in the PLC formed under similar hydrothermal conditions. Episodic reactivation of basement faults was an important driving force to draw uraniferous fluids from the basin and reducing fluids from the basement to the mineralization sites, forming deep basement-hosted deposits.Thematic collection: This article is part of the Uranium Fluid Pathways collection available at: https://www.lyellcollection.org/cc/uranium-fluid-pathwaysSupplementary material:https://doi.org/10.6084/m9.figshare.c.5510179

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

Basement rocks below the Athabasca Basin, Saskatchewan, have been intensely altered through paleoweathering and multiple hydrothermal events, including the formation of world-class unconformity-type uranium deposits. Here, we demonstrate the utility of Ti-oxide thermochronology for identifying thermotectonic events in these altered rocks leading to uranium mineralization along basement structures. Rutile grains along the P2 fault, a major fault in the eastern Athabasca Basin, exhibit 207Pb/206Pb ages of ca. 1850−1700 Ma, with a weighted mean of 1757 ± 6 Ma (mean square of weighted deviation [MSWD] = 1.4, n = 116). The older ages (>1770 Ma) record regional metamorphism reaching a temperature of 875 °C during the Trans-Hudson orogeny. Pb diffusion modeling indicates that metamorphic rutile should exhibit cooling ages of 1760−1750 Ma. Rutile grains showing young ages, <1750 Ma, reflect isotopic resetting during regional asthenospheric upwelling between 1770 and 1730 Ma related to the emplacement of the Kivalliq igneous suite to the north. This thermotectonic event (temperature > 550 °C) promoted hydrothermal activity to produce silicified rocks, i.e., “quartzite,” along the P2 fault, which later focused mineralizing fluids for unconformity-type uranium deposits. The young rutile ages also indicate that the basement rocks remained hot until 1700 Ma, providing the maximum age for the deposition of the Athabasca sediments. Anatase yields a concordia age of 1569 ± 31 Ma (MSWD = 0.30, n = 5), which is within uncertainty of the oldest ages for uraninite of the McArthur River deposit. This age corresponds to the incursion of basinal fluids in the basement along the P2 fault during uranium mineralization.

Supplemental Material: 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

Appendix figures and tables; U-Pb data set for rutile and anatase from the P2 fault.