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.

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.

Documentation

In order to become better prepared to support Research Data Management (RDM) practices in sciences and engineering, Queen’s University Library, together with the University Research Services, conducted a research study of all ranks of faculty members, as well as postdoctoral fellows and graduate students at the Faculty of Engineering & Applied Science, Departments of Chemistry, Computer Science, Geological Sciences and Geological Engineering, Mathematics and Statistics, Physics, Engineering Physics & Astronomy, School of Environmental Studies, and Geography & Planning in the Faculty of Arts and Science.

Paleoecology, Diagenesis and Biological Affinity of Tetradiids in the Upper Ordovician (Sandbian) Black River Group of the Kingston Area, Ontario, Canada

Tetradiids are a group of colonial, tubular fossils that occur globally in Middle to Upper Ordovician strata. Tetradiids were first described as a type of tabulate coral; however, based on their four-fold symmetry, division, and presence of a central-sparry canal, they were recently reinterpreted as a florideophyte rhodophyte algae, a reinterpretation that is tested in this thesis. This study focused on understanding the affinity and taphonomy of this order of fossil. Research was conducted by stratigraphic and petrographic analyses of the Black River Group in the Kingston, Ontario region. Tetradiid occurrences were divided into fragment or colonial, with three morphologies of tetradiids described (Tetradium, Phytopsis and Paratetradium). Morphology is specific to depositional environment, with compact Tetradium consistently within ooid grainstones and open branching Phytopsis and chained Paratetradium consistently within mudstones. Two types of patch reefs were recognized: a Paratetradium bioherm, and a Paratetradium, Phytopsis, stromatolite bioherm. The presence of bioherms implies that tetradiids were capable of hypercalcifying. Preservation styles of tetradiids were investigated, and were compared to brachiopods, echinoderms, mollusks, and ooids. Tetradiids were preferentially preserved as molds and demonstrated complete dissolution of skeletal material. Rare specimens, however, demonstrated preserved horizontal partitions, central plates, and a double wall. Skeletal molds were filled with either calcite spar, mud or encrusted by a cryptomicrobial colony. Both calcitic and aragonitic ooids were discovered. The co-occurrence of aragonitic ooids, aragonitic crytodontids, and the evolution of aragonitic, hypercalcifying tetradiids is interpreted as representing the geochemical favoring of aragonite and HMC in a time of global calcite seas. The geochemical favoring of aragonite is interpreted to be independent to global Mg: Ca ratios, but was the result of increased saturation levels and temperature driven by high atmospheric pCO2. Based on the presence of epitheca, tabulae, septa, and the commonality of growth forms, tetradiids are interpreted as an order of Cnidaria. The evolution of an aragonitic skeleton in tetradiids is interpreted to be the result of de novo acquisition of a skeleton from an unmineralized clade.

Understanding arsenic mobility and speciation in lake sediments impacted by ore roasting near Giant mine, NWT

The roasting of gold-bearing arsenopyrite at Giant mine (Northwest Territories) between 1949 and 1999 released approximately 20,000 tonnes of toxic arsenic-bearing aerosols in the local aerial environment. Detailed examination of lake sediments, sediment porewaters, surface waters and lake hydrology sampled from three lakes of differing limnological characteristics was conducted in summer and winter conditions. Samples were analyzed for solid and dissolved elemental concentrations, speciation and mineralogy. The three lakes are located less than 5km from the mine roaster, and downwind, based on predominant wind direction. The objective of the study was to assess the controls on the mobility and fate of arsenic in these roaster-impacted subarctic lacustrine environments. Results show that the occurrence of arsenic trioxide in lake sediments coincides with the regional onset of industrial activities. The bulk of arsenic in sediments is contained in the form of secondary sulphide precipitates, with iron oxides hosting a minimal amount of arsenic near the surface-water interface. The presence of geogenic arsenic is likely contained as dilute impurities in common rock-forming minerals, and is not believed to be a significant source of arsenic to sediments, porewaters or lake waters. Furthermore, the well correlated depth-profiles of arsenic, antimony and gold in sediments may help reveal roaster impact. The soluble arsenic trioxide particles contained in sediments act as the primary source of arsenic into porewaters. Dissolved arsenic in reducing porewaters both precipitate as secondary sulphides in situ, and diffuse upwards into the overlying lake waters. Arsenic diffusion out of porewaters, combined with watercourse-driven residence time, are estimated to be the predominant mechanisms controlling arsenic concentrations in overlying lake waters. The sequestration of arsenic from porewaters as sulphide precipitates, in the study lakes, is not an effective process in keeping lake-water arsenic concentrations below guidelines for the protection of the freshwater environment and drinking water. Seasonal impacts on lake geochemistry derive from ice covering lake waters, cutting them off from of atmospheric oxygen, along with the exclusion of solutes from the ice. Such effects are limited in deep lakes but are can be an important factor controlling arsenic precipitation and mobility in ponds.

Modeling the impact of large-scale solid mass transfers on the gravity field

The focus of this thesis is to explore and quantify the response of large-scale solid mass transfer events on satellite-based gravity observations. The gravity signature of large-scale solid mass transfers has not been deeply explored yet; mainly due to the lack of significant events during dedicated satellite gravity missions‘ lifespans. In light of the next generation of gravity missions, the feasibility of employing satellite gravity observations to detect submarine and surface mass transfers is of importance for geoscience (improves the understanding of geodynamic processes) and for geodesy (improves the understanding of the dynamic gravity field). The aim of this thesis is twofold and focuses on assessing the feasibility of using satellite gravity observations for detecting large-scale solid mass transfers and on modeling the impact on the gravity field caused by these events. A methodology that employs 3D forward modeling simulations and 2D wavelet multiresolution analysis is suggested to estimate the impact of solid mass transfers on satellite gravity observations. The gravity signature of various submarine and subaerial events that occurred in the past was estimated. Case studies were conducted to assess the sensitivity and resolvability required in order to observe gravity differences caused by solid mass transfers. Simulation studies were also employed in order to assess the expected contribution of the Next Generation of Gravity Missions for this application.

Time-dependency of Rocks and Implications Associated with Tunnelling

The formulation of a geotechnical model and the associated prediction of the mechanical behaviour is a challenge engineers need to overcome in order to optimize tunnel design and meet project requirements. Special challenges arise in cases where rocks and rockmasses are susceptible to time-effects and time-dependent processes govern. Progressive rockmass deformation and instability, time-dependent overloading of support and delayed failures are commonly the result of time-dependent phenomena. The research work presented in this thesis serves as an attempt to provide more insight into the time-dependent behaviour of rocks. Emphasis is given on investigating and analyzing creep deformation and time-dependent stress relaxation phenomenon at the laboratory scale and in-depth analyses are presented. This thesis further develops the understanding of these phenomena and practical yet scientific tools for estimating and predicting the long-term strength and the maximum stress relaxation of rock materials are proposed. The identification of the existence of three distinct behavioural stages during stress relaxation is presented and discussed. The main observations associated with time-dependent behaviour are employed in numerical analyses and applied at the tunnel scale. A new approach for simulating and capturing the time-dependent behaviour coupled with the tunnel advancement effect is also developed and analyzed. Guidance is provided to increase the understanding of the support-rockmass interaction and the main implications and significance of time-dependent behaviour associated with rock tunnelling are discussed. The work presented in this thesis advances the scientific understanding of time-dependent rock and rockmass behaviour, increases the awareness of how such phenomena are captured numerically, and lays out a framework for dealing with such deformations when predicting tunnel deformations. Practical aspects of this thesis are also presented, which will increase their usage in the associated industries and close the gap between the scientific and industry communities.

On the feasibility of time-lapse superconducting gravimetry for reservoir monitoring

The feasibility of monitoring fluid flow subsurface processes that result in density changes, using the iGrav superconducting gravimeter, is investigated. Practical targets include steam-assisted gravity drainage (SAGD) bitumen depletion and water pumping from aquifers, for which there is currently a void in low-impact, inexpensive monitoring techniques. This study demonstrates that the iGrav has the potential to be applied to multi-scale and diverse reservoirs. Gravity and gravity gradient signals are forward modeled for a real SAGD reservoir at two time steps, and for surface-fed and groundwater-fed aquifer pumping models, to estimate signal strength and directional dependency of water flow. Time-lapse gravimetry on small-scale reservoirs exhibits two obstacles, namely, a µgal sensitivity requirement and high noise levels in the vicinity of the reservoir. In this study, both limitations are overcome by proposing (i) a portable superconducting gravimeter, and (ii) a pair of instruments under various baseline geometries. This results in improved spatial resolution for locating depletion zones, as well as the cancellation of noise common in both instruments. Results indicate that a pair of iGrav superconducting gravimeters meet the sensitivity requirements and the spatial focusing desired to monitor SAGD bitumen migration at the reservoir scales. For SAGD reservoirs, the well pair separation, reservoir depth, and survey sampling determine the resolvability of individual well pair depletion patterns during the steam chamber rising phase, and general reservoir depletion patterns during the steam chamber spreading phase. Results show that monitoring water table elevation changes due to pumping and tracking whether groundwater or surface water is being extracted are feasible.

Building a Rockfall Database Using Remote Sensing: Techniques for Hazard Management in Canadian Rail Corridors

The section of CN railway between Vancouver and Kamloops runs along the base of many hazardous slopes, including the White Canyon, which is located just outside the town of Lytton, BC. The slope has a history of frequent rockfall activity, which presents a hazard to the railway below. Rockfall inventories can be used to understand the frequency-magnitude relationship of events on hazardous slopes, however it can be difficult to consistently and accurately identify rockfall source zones and volumes on large slopes with frequent activity, leaving many inventories incomplete. We have studied this slope as a part of the Canadian Railway Ground Hazard Research Program and have collected remote sensing data, including terrestrial laser scanning (TLS), photographs, and photogrammetry data since 2012, and used change detection to identify rockfalls on the slope. The objective of this thesis is to use a subset of this data to understand how rockfalls identified from TLS data could be used to understand the frequency-magnitude relationship of rockfalls on the slope. This includes incorporating both new and existing methods to develop a semi-automated workflow to extract rockfall events from the TLS data. We show that these methods can be used to identify events as small as 0.01 m3 and that the duration between scans can have an effect on the frequency-magnitude relationship of the rockfalls. We also show that by incorporating photogrammetry data into our analysis, we can create a 3D geological model of the slope and use this to classify rockfalls by lithology, to further understand the rockfall failure patterns. When relating the rockfall activity to triggering factors, we found that the amount of precipitation occurring over the winter has an effect on the overall rockfall frequency for the remainder of the year. These results can provide the railways with a more complete inventory of events compared to records created through track inspection, or rockfall monitoring systems that are installed on the slope. In addition, we can use the database to understand the spatial and temporal distribution of events. The results can also be used as an input to rockfall modelling programs.

Data Files

In order to become better prepared to support Research Data Management (RDM) practices in sciences and engineering, Queen’s University Library, together with the University Research Services, conducted a research study of all ranks of faculty members, as well as postdoctoral fellows and graduate students at the Faculty of Engineering & Applied Science, Departments of Chemistry, Computer Science, Geological Sciences and Geological Engineering, Mathematics and Statistics, Physics, Engineering Physics & Astronomy, School of Environmental Studies, and Geography & Planning in the Faculty of Arts and Science.

The Provenance and Geochemical Alteration of Bermudan Eolianites

Large calcareous eolianites cover the remote island of Bermuda, accounting for more than 90% of the limestone bedrock. This study examines the sedimentology and geochemistry of these eolianites to better understand Pleistocene oceanography and the meteoric alteration of subtropical carbonate sediments. Cluster analyses reveal that the eolian carbonate sediments fall into two natural groups that represent lagoonal and reefal end members of marine sediment production. Coral fragments are uncharacteristically absent, possibly destroyed prior to their incorporation into eolian deposits by endolithic microboring organisms or broken up during transport. Sediment assemblages lead to the following interpretations of the Bermudan offshore environment: (1) the Ledge Flats reef system along the southwestern coast has been active since MIS 11, contributing coralline algal-rich sediment to the northern beaches of Sandy’s Parish and acting as an energy barrier in the south, allowing for low energy sedimentation in the quiet back- reef region; (2) on the northeastern coast, the low energy back-reef region landward of the Ledge Flats has thrived since MIS 11; (3) during MIS 5e, slightly warmer water temperatures led to the hindrance of coralline algal growth along the southern coast and in the North Lagoon. These are the first interpretations of Pleistocene marine assemblages on Bermuda. Meteoric fluids progressively transformed the pristine carbonate sediments into hardened limestones in a predictable solubility-dependent manner. The progressive alteration is coincident with: (1) divergence of δ18O and δ13C values from those similar to unaltered sediment towards those of calcrete, due to interaction with CO2-charged meteoric fluids; (2) depletion of elements with low partitioning coefficients and low meteoric concentrations, such as barium, boron, magnesium, potassium, sodium, strontium, and uranium; (3) enrichment of iron from Terra Rossa-hosted iron oxides; (4) enrichment of aluminum via detrital minerals sourced from protosol horizons; and (5) manganese concentrations that remain uncharacteristically low, owing to the lack of a consistent manganese source. Elemental correlations are useful for characterizing meteoric diagenesis, assuming the primary mineralogy is recognized, all components have been fully altered, and inter-particle cements are ubiquitous.

Karst hazard analysis towards the development of predictive karst models for southern Ontario

An investigation into karst hazard in southern Ontario has been undertaken with the intention of leading to the development of predictive karst models for this region. The reason these are not currently feasible is a lack of sufficient karst data, though this is not entirely due to the lack of karst features. Geophysical data was collected at Lake on the Mountain, Ontario as part of this karst investigation. This data was collected in order to validate the long-standing hypothesis that Lake on the Mountain was formed from a sinkhole collapse. Sub-bottom acoustic profiling data was collected in order to image the lake bottom sediments and bedrock. Vertical bedrock features interpreted as solutionally enlarged fractures were taken as evidence for karst processes on the lake bottom. Additionally, the bedrock topography shows a narrower and more elongated basin than was previously identified, and this also lies parallel to a mapped fault system in the area. This suggests that Lake on the Mountain was formed over a fault zone which also supports the sinkhole hypothesis as it would provide groundwater pathways for karst dissolution to occur. Previous sediment cores suggest that Lake on the Mountain would have formed at some point during the Wisconsinan glaciation with glacial meltwater and glacial loading as potential contributing factors to sinkhole development. A probabilistic karst model for the state of Kentucky, USA, has been generated using the Weights of Evidence method. This model is presented as an example of the predictive capabilities of these kind of data-driven modelling techniques and to show how such models could be applied to karst in Ontario. The model was able to classify 70% of the validation dataset correctly while minimizing false positive identifications. This is moderately successful and could stand to be improved. Finally, suggestions to improving the current karst model of southern Ontario are suggested with the goal of increasing investigation into karst in Ontario and streamlining the reporting system for sinkholes, caves, and other karst features so as to improve the current Ontario karst database.

Provenance of Holocene Calcareous Aeolianites, Eastern South Australia

The southeastern coast of South Australia contains a spectacular and world-renown suite of Quaternary calcareous aeolianites. This study is focused on the provenance of components in the Holocene sector of this carbonate breach-dune succession. Research was carried out along seven transects from ~30 meters water depth offshore across the beach and into the dunes. Offshore sediments were acquired via grab sampling and SCUBA. Results indicate that dunes of the southern Lacepede and Otway coasts in particular are mostly composed of modern invertebrate and calcareous algal allochems. The most numerous grains are from molluscs, benthic foraminifera, coralline algae, echinoids, and bryozoans. These particles originate in carbonate factories such as macroalgal forests, rocky reefs, seagrass meadows, and low-relief seafloor rockgrounds. The incorporation of carbonate skeletons into coastal dunes, however, depends on a combination of; 1) the infauna within intertidal and nearshore environments, 2) the physical characteristics of different allochems and their ability to withstand fragmentation and abrasion, 3) the wave and swell climate, and 4) the nature of aeolian transport. Most aeolian dune sediment is derived from nearshore and intertidal carbonate factories. This is particularly well illustrated by the abundance of robust infaunal bivalves that inhabit the nearshore sands and virtual absence of bryozoans that are common as sediment particles in water depths > 10mwd. Thus, the calcareous aeolianites in this cool-water carbonate region are not a reflection of the offshore marine shelf factories, but more a product of shallow nearshore-intertidal biomes.

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.

Heavy Minerals in Soils from the Athabasca Basin and the Implications for Exploration Geochemistry of Uranium Deposits at Depth

The Centennial deposit is a high grade (~8% U3O8), deeply buried (~950m), unconformity-related U deposit located in the south-central region of the Athabasca Basin in northern Saskatchewan, Canada. The mineral chemistry of fine fractions (<63 μm) of soils from grids above the Centennial deposit were examined to understand possible controls on the geochemistry and radiogenic 207Pb/206Pb ratios measured in the clay-size (<2 μm) fractions used for exploration. Soil samples distal and proximal to the deposit projection to the surface and geophysically defined structures were selected. Mineral abundances were determined using the scanning electron microscope and Mineral Liberation Analysis. Zircon was the only U-rich mineral identified with modal abundances >0.02% by weight. Monazite, which can be U-rich, was identified, but not in significant abundances. The source of the zircon and other heavy minerals is interpreted to be from sub-cropping sources that are >100 km up-ice from Centennial. Trace element analysis using laser ablation inductively coupled plasma mass spectrometry of hydroseparated zircon grains indicate that zircon abundances and zircon Pb concentrations in surficial samples have minimal effect on the radiogenic 207Pb/206Pb ratios in the clay-fraction of the samples, with the dominant source of radiogenic Pb being clay mineral surfaces that trapped Pb during secondary dispersion from the Centennial uranium deposit through faults and fractures to the surface. The REE patterns indicate HREE enrichment in the clay-fractions of samples that have higher abundances of zircon in the <20 μm fraction. Immobile elements such as HREE that are concentrated in zircon can be used as indicators of radiogenic Pb being sourced from minerals at the surface rather than being sourced from secondary dispersion from deeply buried U deposits.