The Evolution of the Los Negritos Porphyry Copper Deposit, Coquimbo, Chile: Geological, Geochronological, Mineralogical and Geochemical Evidence

The Los Negritos porphyry copper deposit is located ~ 4 km to the northeast of Carmen de Andacollo Mine in the Chilean Cretaceous metallogenic belt. The mineralization is hosted in andesite of the Quebrada Marquesa Formation and a series of at least four early to intramineral porphyry intrusive rock types: plagioclase quartz biotite porphyry (P1b and P1a dated at 109.60± 0.75 Ma and 107.22± 0.40 Ma); plagioclase biotite porphyry (P2: 106.30 ± 0.47 Ma); and quartz plagioclase biotite porphyry (P3: 106.19 ± 0.42 Ma). These units are cut by late‐ to post‐mineral plagioclase‐hornblende porphyritic rocks (P4b: 106.20 ± 0.69 Ma and P4a: 106.50 ± 0.68 Ma). The earliest intrusive units (P1) were affected by an initial stage of K‐feldspar‐biotite alteration, with chalcopyrite, molybdenite (date at 108.5 ± 0.5 Ma) and gold (up to 0.11 ppm), and the surrounding volcanic host rock was overprinted by chlorite‐epidote dominated (propylitic) alteration. Subsequent to the P2 and P3 intrusion, these rocks were affected by albite and then a second stage of potassic alteration. The Ti and Ba contents in hydrothermal biotite are notably lower (typically Ti = 0.100‐0.144 a.p.f.u. and Ba = 0.001‐0.005 a.p.f.u) than in magmatic ones (generally Ti = 0.186‐0.222 a.p.f.u. and Ba = 0.014‐0.023 a.p.f.u.), and constitute an excellent discriminant of the nature of biotite. These early stages of alteration were overprinted by copper‐molybdenum bearing chlorite‐sericite alteration at 106.60 ± 0.5 Ma (Re‐Os age in molybdenite) and by quartz‐sericite‐pyrite veins (phyllic), respectively in the southwest and northeast areas. The average temperature associated with these two alteration facies is estimated around 305 °C. Weak albite‐calcite alteration, spatially associated with sulfosalts and distributed along the margins of P3, overprinted the phyllic facies. The intrusive rock units at the Los Negritos and Carmen de Andacollo deposits are geochemically classified as diorite to granodiorite with a calc‐alkaline magmatic affinity, and formed in a volcanic arc setting from partial melting of a metasomatized mantle wedge. They are interpreted to be cogenetic, and related to a common long‐lived magma chamber that emplaced during a period of tectonic inversion known as the Subhercynian, Peruvian or Pacific event.

CHARACTERIZATION AND THERMAL-HYDRAULIC MODELING OF ENGINEERED BARRIERS FOR A CANADIAN DEEP GEOLOGICAL REPOSITORY FOR USED NUCLEAR FUEL

Within Canada there are more than 2.5 million bundles of spent nuclear fuel with another approximately 2 million bundles to be generated in the future. Canada, and every country around the world that has taken a decision on management of spent nuclear fuel, has decided on long-term containment and isolation of the fuel within a deep geological repository. At depth, a deep geological repository consists of a network of placement rooms where the bundles will be located within a multi-layered system that incorporates engineered and natural barriers. The barriers will be placed in a complex thermal-hydraulic-mechanical-chemical-biological (THMCB) environment. A large database of material properties for all components in the repository are required to construct representative models. Within the repository, the sealing materials will experience elevated temperatures due to the thermal gradient produced by radioactive decay heat from the waste inside the container. Furthermore, high porewater pressure due to the depth of repository along with possibility of elevated salinity of groundwater would cause the bentonite-based materials to be under transient hydraulic conditions. Therefore it is crucial to characterize the sealing materials over a wide range of thermal-hydraulic conditions. A comprehensive experimental program has been conducted to measure properties (mainly focused on thermal properties) of all sealing materials involved in Mark II concept at plausible thermal-hydraulic conditions. The thermal response of Canada’s concept for a deep geological repository has been modelled using experimentally measured thermal properties. Plausible scenarios are defined and the effects of these scenarios are examined on the container surface temperature as well as the surrounding geosphere to assess whether they meet design criteria for the cases studied. The thermal response shows that if all the materials even being at dried condition, repository still performs acceptably as long as sealing materials remain in contact.

Multi-Scale Visualization and Interpretation of Geological Relationships

Visualization and interpretation of geological observations into a cohesive geological model are essential to Earth sciences and related fields. Various emerging technologies offer approaches to multi-scale visualization of heterogeneous data, providing new opportunities that facilitate model development and interpretation processes. These include increased accessibility to 3D scanning technology, global connectivity, and Web-based interactive platforms. The geological sciences and geological engineering disciplines are adopting these technologies as volumes of data and physical samples greatly increase. However, a standardized and universally agreed upon workflow and approach have yet to properly be developed. In this thesis, the 3D scanning workflow is presented as a foundation for a virtual geological database. This database provides augmented levels of tangibility to students and researchers who have little to no access to locations that are remote or inaccessible. A Web-GIS platform was utilized jointly with customized widgets developed throughout the course of this research to aid in visualizing hand-sized/meso-scale geological samples within a geologic and geospatial context. This context is provided as a macro-scale GIS interface, where geophysical and geodetic images and data are visualized. Specifically, an interactive interface is developed that allows for simultaneous visualization to improve the understanding of geological trends and relationships. These developed tools will allow for rapid data access and global sharing, and will facilitate comprehension of geological models using multi-scale heterogeneous observations.

Axis Mapping: The Estimation of Surface Orientations and its Applications in Vehicle Localization and Structural Geology

The map representation of an environment should be selected based on its intended application. For example, a geometrically accurate map describing the Euclidean space of an environment is not necessarily the best choice if only a small subset its features are required. One possible subset is the orientations of the flat surfaces in the environment, represented by a special parameterization of normal vectors called axes. Devoid of positional information, the entries of an axis map form a non-injective relationship with the flat surfaces in the environment, which results in physically distinct flat surfaces being represented by a single axis. This drastically reduces the complexity of the map, but retains important information about the environment that can be used in meaningful applications in both two and three dimensions. This thesis presents axis mapping, which is an algorithm that accurately and automatically estimates an axis map of an environment based on sensor measurements collected by a mobile platform. Furthermore, two major applications of axis maps are developed and implemented. First, the LiDAR compass is a heading estimation algorithm that compares measurements of axes with an axis map of the environment. Pairing the LiDAR compass with simple translation measurements forms the basis for an accurate two-dimensional localization algorithm. It is shown that this algorithm eliminates the growth of heading error in both indoor and outdoor environments, resulting in accurate localization over long distances. Second, in the context of geotechnical engineering, a three-dimensional axis map is called a stereonet, which is used as a tool to examine the strength and stability of a rock face. Axis mapping provides a novel approach to create accurate stereonets safely, rapidly, and inexpensively compared to established methods. The non-injective property of axis maps is leveraged to probabilistically describe the relationships between non-sequential measurements of the rock face. The automatic estimation of stereonets was tested in three separate outdoor environments. It is shown that axis mapping can accurately estimate stereonets while improving safety, requiring significantly less time and effort, and lowering costs compared to traditional and current state-of-the-art approaches.

Towards Mapping of Rock Walls Using a UAV-Mounted 2D Laser Scanner in GPS Denied Environments

In geotechnical engineering, the stability of rock excavations and walls is estimated by using tools that include a map of the orientations of exposed rock faces. However, measuring these orientations by using conventional methods can be time consuming, sometimes dangerous, and is limited to regions of the exposed rock that are reachable by a human. This thesis introduces a 2D, simulated, quadcopter-based rock wall mapping algorithm for GPS denied environments such as underground mines or near high walls on surface. The proposed algorithm employs techniques from the field of robotics known as simultaneous localization and mapping (SLAM) and is a step towards 3D rock wall mapping. Not only are quadcopters agile, but they can hover. This is very useful for confined spaces such as underground or near rock walls. The quadcopter requires sensors to enable self localization and mapping in dark, confined and GPS denied environments. However, these sensors are limited by the quadcopter payload and power restrictions. Because of these restrictions, a light weight 2D laser scanner is proposed. As a first step towards a 3D mapping algorithm, this thesis proposes a simplified scenario in which a simulated 1D laser range finder and 2D IMU are mounted on a quadcopter that is moving on a plane. Because the 1D laser does not provide enough information to map the 2D world from a single measurement, many measurements are combined over the trajectory of the quadcopter. Least Squares Optimization (LSO) is used to optimize the estimated trajectory and rock face for all data collected over the length of a light. Simulation results show that the mapping algorithm developed is a good first step. It shows that by combining measurements over a trajectory, the scanned rock face can be estimated using a lower-dimensional range sensor. A swathing manoeuvre is introduced as a way to promote loop closures within a short time period, thus reducing accumulated error. Some suggestions on how to improve the algorithm are also provided.

Alteration mineralogy and pathfinder element inventory of the McArthur River unconformity-related uranium deposit, Canada

The chemical compositions, modal mineralogy, and textural variability of interstitial minerals in sandstones of the Athabasca Group strata in the vicinity of the McArthur River unconformity-related uranium deposit were characterized using a combination of short wave infrared spectroscopy (SWIR), lithogeochemistry, scanning electron microscopy (SEM), electron probe microanalysis (EPMA) and laser ablation mass spectrometry (LA-ICP-MS) to determine the residence sites of pathfinder trace elements. The importance of integrating in-situ mineral chemistry with whole-rock analyses resides in the possibility to establish the mineralogical and paragenetic context of geochemical signatures in defining the footprint of the deposit. Located in the Athabasca Basin, Saskatchewan, Canada, the deposit is situated below ~550 m of quartz arenitic sandstones that are strongly silicified between depths of approximately 200-400 m. The silicified layer exhibits significant control on the distribution of alteration minerals, and appears to have restricted both the primary and secondary dispersion of pathfinder trace elements, which include U, radiogenic Pb isotopes, V, Ni, Co, Cu, Mo, As, Zn, and REEs. Diagenetic background sandstones contain assemblages of illite, dickite, aluminum-phosphate-sulfate (APS) minerals, apatite, and Fe-Ti oxide minerals. Altered sandstones contain assemblages of Al-Mg chlorite (sudoite), alkali-deficient dravite, APS minerals, kaolinite, illite, and oxide minerals. Throughout the sandstones, APS minerals account for the majority of the Sr and LREE concentrations, whereas late pre-ore chlorite, containing up to 0.1 wt.% Ni, accounts for the majority of Ni concentrations. Cobalt, Cu, Mo, and Zn occur predominantly in cryptic sub-micron sulfide and sulfarsenide inclusions in clay mineral aggregates and in association with paragenetically-late Fe-Ti oxides. Uranium occurs predominantly in cryptic micro-inclusions associated with pyrite in late-stage quartz overgrowths, and with paragenetically late Fe-Ti oxide micro-inclusions in kaolinite. Additionally, up to 0.2 wt.% U is cryptically distributed in post-ore Fe-oxide veins. Early diagenetic apatite, monazite and apatite inclusions in detrital quartz, and detrital zircon also contribute significant U and HREE to samples analyzed with an aggressive leach such as Aqua Regia. Detailed LA-ICP-MS chemical mapping of interstitial assemblages, detrital grains, and cements provides new insights into the distribution and inventory of pathfinder elements in the footprint of the McArthur River uranium deposit.