Structural complexity inferred from anisotropic resistivity: Example from airborne EM and compilation of historical resistivity/induced polarization data from the gold-rich Canadian Malartic district, Québec, Canada

This paper is © 2019 Society of Exploration Geophysicists. The posting is available free of charge and its use is subject to the SEG terms and conditions: https://seg.org/Terms-of-UseStructurally complex zones within orogenic terranes typically correspond to areas where there is interference between multiple fold generations and are known to be favorable pathways for fluid flow because of their higher permeability. In the Canadian Malartic district, gold anomalies have been linked with zones of structural complexity that have been quantified by outcrop bedding orientation measurements and calculation of bedding variance maps. In this work, historical apparent resistivity and induced polarization data in the Canadian Malartic district were reprocessed and combined with new surveys to create a compilation of inverted chargeability and resistivity, which were then interpreted together with airborne electromagnetics and outcrop structural data. The results indicate chargeability anomalies, up to five times the background value, associated with the sulfide mineral content in monzodioritic dikes that are thickened in folds and hydrothermally altered. Although the airborne apparent half-space resistivity is mostly sensitive to conductive surficial cover, the inverted ground resistivity method is sensitive to deeper structure and likely represents bedrock signal at depths greater than 25 m. Inverted ground resistivity exhibits strong anisotropy in areas of subvertical bedding, where measured resistivities can vary by up to a factor of two, over the same location, depending on whether the survey lines are perpendicular or parallel to the strike of bedding. This result is observed at scales of 50 cm up to 100 m. Analysis of inverted ground resistivity together with bedding variance indicates a strong correlation between structurally complex zones with high bedding variance and a decrease in resistivity at depths greater than 25 m. This suggests that in places where the presence of disseminated gold cannot be directly detected, or where the outcrop exposure is limited due to overburden cover, geophysical data may still succeed in identifying structural complexity zones that could potentially host mineralization.Natural Sciences and Engineering Research Council of Canada and the Canada Mining Innovation Council (NSERC-CMIC Mineral Exploration Footprints Project Contribution 178

Factors influencing consumer wine choice: The case of wine tourism

This research was undertaken to examine the relationship between the wine tourism activities of consumers and their subsequent wine choices. Exploratory research was undertaken giving specific attention to wine tourism, first, as a factor of influence on decision-making, second, on levels of consumer involvement and third on long-term consumer loyalty. A convenience sample of 12 wine consumers in Dublin participated in the study. Semi-structured in-depth interviews were used to gather the qualitative data used. The findings clearly establish links made by consumers between their wine tourism experiences and their subsequent wine purchase preferences. The findings highlight the importance of facilitating the co-creation of memorable wine tourism experiences as a determinant of improved consumer involvement as well as the development of long-term consumer loyalty. Implications for industry practitioners suggest a continued focus on the provision of unique wine tourism experiences; while the findings support the theoretical links between consumers own past experiences and their use as a valuable information source when making purchase decisions

Palaeotectonic setting of the south-eastern Kédougou-Kéniéba Inlier, West Africa: new insights from igneous trace element geochemistry and U-Pb zircon ages

New U-Pb zircon ages and geochemistry from the eastern Kédougou-Kéniéba Inlier are presented and integrated with published data to generate a revised tectonic framework for the westernmost Birimian terranes. The Falémé Volcanic Belt and Kofi Series are highly prospective, hosting several multi-million ounce gold deposits and a significant iron ore resource, but remain under-researched. It is therefore important to constrain the fundamental geological setting. The igneous rocks of the eastern Kédougou-Kéniéba Inlier are dominantly of high-K calc-alkaline affinity, with fractionated REE patterns and negative Nb-Ta anomalies. The plutonic rocks in the Falémé Belt are dioritic to granodioritic in composition, with moderately fractionated REE patterns and metaluminous A/CNK signatures. Felsic, peraluminous granite stocks, dykes and plutons with fractionated REE patterns and negative Eu, Ti and P anomalies intruded both the Falémé Belt and Kofi Series. Albitisation masks the affinity of some units, although use of the Th-Co diagram shows that prior to albitisation, all igneous units belonged to the high-K calc-alkaline series. New U-Pb age data for the Boboti and Balangouma plutons indicate crystallisation at 2088.5 ± 8.5 Ma and at 2112 ± 13 Ma, respectively. Inherited zircons in the Boboti pluton indicate magmatic activity in the Falémé Belt at 2218 ± 83 Ma coincided with the oldest dated units in the Mako Belt to the West. Systematic changes in Dy/Yb, Sm/La, Nb/Zr, Rb concentration, Eu-anomaly and ɛNdt over ∼200 Ma reveal that the tectonic setting in the KKI evolved from a volcanic island arc environment to an active continental margin. Crustal thickening, as a result of a shift to collisional tectonic setting, combined with magmatic differentiation, led to the generation of peraluminous, granitic melts with a significant crustal component. A small suite of more basic intrusive and extrusive rocks on the eastern margin of the Dialé-Daléma basin are highly metaluminous and display limited LILE enrichment, with normalised HREE values close to unity. The Daléma igneous rocks may have formed in an extensional back arc, related to the arc system

Integrated Multi-Parameter Exploration Footprints of the Canadian Malartic Disseminated Au, McArthur River-Millennium Unconformity U, and Highland Valley Porphyry Cu Deposits: Preliminary Results from the NSERC-CMIC Mineral Exploration Footprints Research Network

Mineral exploration in Canada is increasingly focused on concealed and deeply buried targets, requiring more effective tools to detect large-scale ore-forming systems and to vector from their most distal margins to their high grade cores. A new generation of ore system models is required to achieve this. The Mineral Exploration Footprints Research Network is a consortium of 70 faculty, research associates, and students from 20 Canadian universities working with 30 mining, mineral exploration, and mining service providers to develop new approaches to ore system modelling based on more effective integration and visualization of multi-parameter geological-structural-mineralogical-lithogeochemical-petrophysical-geophysical exploration data. The Network is developing the next generation ore system models and exploration strategies at three sites based on integrated data visualization using self-consistent 3D Common Earth Models and geostatistical/machine learning technologies. Thus far over 60 footprint components and vectors have been identified at the Canadian Malartic stockwork-disseminated Au deposit, 20–30 at the McArthur-Millennium unconformity U deposits, and over 20 in the Highland Valley porphyry Cu system. For the first time, these are being assembled into comprehensive models that will serve as landmark case studies for data integration and analysis in the today’s challenging exploration environment

Development and application of feature engineered geological layers for ranking magmatic, volcanogenic, and orogenic system components in Archean greenstone belts

Geologically representative feature engineering is a crucial component in geoscientific applications of machine learning. Many commonly applied feature engineering techniques used to produce input variables for machine learning apply geological knowledge to generic data science techniques, which can lead to ambiguity, geological oversimplification, and/or compounding subjective bias. Workflows that utilize minimally processed input variables attempt to overcome these issues, but often lead to convoluted and uninterpretable results. To address these challenges, new and enhanced feature engineering methods were developed by combining geological knowledge, understanding of data limitations, and a variety of data science techniques. These include non-Euclidean fluid pre-deformation path distance, rheological and chemical contrast, geologically constrained interpolation of characteristic host rock geochemistry, interpolation of mobile element gain/loss, assemblages, magnetic intensity, structural complexity, host rock physical properties. These methods were applied to compiled open-source and new field observations from Archean greenstone terranes in the Abitibi and western Wabigoon sub-provinces of the Superior Province near Timmins and Dryden, Ontario, respectively. Resulting feature maps represent conceptually significant components in magmatic, volcanogenic, and orogenic mineral systems. A comparison of ranked feature importance from random forests to conceptual mineral system models show that the feature maps adequately represent system components, with a few exceptions attributed to biased training data or limited constraint data. The study also highlights the shared importance of several highly ranked features for the three mineral systems, indicating that spatially related mineral systems exploit the same features when available. Comparing feature importance when classifying orogenic Au mineralization in Timmins and Dryden provides insights into the possible cause of contrasting endowment being related to fluid source. The study demonstrates that integrative studies leveraging multi-disciplinary data and methodology have the potential to advance geological understanding, maximize data utility, and generate robust exploration targets

Revised Eburnean geodynamic evolution of the gold-rich southern Ashanti Belt, Ghana, with new field and geophysical evidence of pre-Tarkwaian deformations

Integration of regional geophysical datasets and detailed field observations provide new insights into the paleoproterozoic structural evolution of southwestern Ghana. The study area is dominated by three metavolcanic and metasedimentary packages known as the Sefwi Group, the Kumasi Group (Birimian) and the Tarkwa Group (Tarkwaian) that were intruded by abundant TTG granitoids during the Eoeburnean and Eburnean phases of an event termed the "Eburnean Orogeny". This study identifies an Eoeburnean (pre-Tarkwaian) deformation event (D1) that produced significant deformation in the Sefwi Group metavolcanics. D1 is associated with N-S shortening manifested as regional scale folding in the southern Ashanti Belt. D1 synorogenic granitoids were intruded between 2187 Ma and 2158 Ma under greenschist metamorphic condition. Syn-D1 gold mineralisation associated with quartz veining could be the original source of Tarkwaian paleo-placers and/or remobilised gold concentrations along major shear zones. D2 represents an extensional phase associated with the Kumasi Group sedimentation (2154-2125 Ma) which could be related to activation of major structures such as the Ashanti Fault as low angle detachments that controlled the deposition of the Kumasi Group and the opening of the Kumasi and Akyem Basin. The Tarkwa Group (2107-2097 Ma) unconformably overlies the Birimian Supergroups and was deposited in response to D3 shortening. D3 resulted in the inversion of syn-D2 detachments faults within the Ashanti Belt. NW-SE D3 shortening produced regional scale folding within the Birimian and the Tarkwaian metasediments. D4 deformation corresponds with sinistral reactivation of D3 thrust faults, and is locally associated with macro-scale folding at Obuasi and Wassa gold mines. By the end of D4, the regional scale architecture was built and was only slightly modified by the two last events. D5 postdates the Eburnean metamorphic peak and corresponds to open recumbent folds associated with a subhorizontal crenulation cleavage. D6 is present as a subvertical crenulation cleavage and reverse faults associated with NE-SW shortening

Geodiversity : exploration of 3D geological model space

The process of building a 3D model necessitates the reconciliation of field observations, geophysical interpretation, geological data uncertainty and the prevailing tectonic evolution hypotheses and interpretations. Uncertainty is compounded when clustered data points collected at local scales are statistically upscaled to one or two points for use in regional models. Interpretation is required to interpolate between sparse field data points using ambiguous geophysical data in covered terranes. It becomes clear that multiple interpretations are possible during model construction. The various interpretations are considered as potential natural representatives, but pragmatism typically dictates that just a single interpretation is offered by the modelling process. Uncertainties are introduced into the 3D model during construction from a variety of sources and through data set optimisation that produces a single model. Practices such as these are likely to result in a model that does not adequately represent the target geology. A set of geometrical 'geodiversity' metrics are used to analyse a 3D model of the Gippsland Basin, southeastern Australia after perturbing geological input data via uncertainty simulation. The resulting sets of perturbed geological observations are used to calculate a suite of geological 3D models that display a range of geological architectures. The concept of biodiversity has been adapted for the geosciences to quantify geometric variability, or geodiversity, between models in order to understand the effect uncertainty has models geometry. Various geometrical relationships (depth, volume, contact surface area, curvature and geological complexity) are used to describe the range of possibilities exhibited throughout the model suite. End-member models geodiversity metrics are classified in a similar manner to taxonomic descriptions. Further analysis of the model suite is performed using principal component analysis (PCA) to determine important geometrical characteristics. The configuration of the model space is determined through identifying 'outlier' model examples, which potentially represent undiscovered model 'species'

Tectonics, mineralisation and regolith evolution of the West African Craton

Studies of mafic dyke swarms may simultaneously provide information on the mechanical, geochemical, geochronological and magnetic environments at the time of their formation. The mafic intrusive history of different cratons can also be potentially used to unravel their assembly into their current configuration. The identification and classification of dykes is a first step to all these studies. Fortunately, even in regions with poor outcrop, we can use the strong magnetic response of mafic dykes to identify and map their extent. In West Africa the first maps of mafic dyke distribution were made over 40 years ago, but there are still large areas where there are almost no published data. In this paper we present a significantly updated map of mafic dykes for the West Africa Craton based in large part on new interpretations of the regional airborne magnetic database. This map includes the locations of over three thousand dykes across the craton, which locally shows several orientation clusters that provide a minimum estimate for the total number of dyke swarms in this region. Whilst we will have to wait until systematic dating of the different swarms is completed, we can demonstrate that there is a long and complex history of mafic magmatism across the craton, with up to 26 distinct dyke swarms mapped based according to their orientation. The mapping and dating of these swarms will provide key constraints on the assembly of the fragments that make up the modern continents

The nature of the southern West African craton lithosphere inferred from its electrical resistivity

The West-African craton is defined by a combination of Archean and Palaeoproterozoic rocks that stabilised at ~2 Ga towards the end of the Paleoproterozoic Eburnean Orogeny, and therefore may reflect the transition from Archean to modern tectonic processes. Exploring its present lithospheric architecture aids further understanding of not only the craton’s stability through its history but also its formation. We investigate the lithospheric structure of the craton through analysing and modelling magnetotelluric (MT) data from a 500-km-long east-west profile in northern Ghana and southern Burkina Faso crossing part of the Baoulé-Mossi Domain and reaching the Volta Basin in the south-eastern part of the craton. Although the MT stations are along a 2D profile, due to the complexity of the structures characterising the area, 3D resistivity modelling of the data is performed to obtain insights on the thermal signature and composition of the subcontinental lithosphere beneath the area. The thermal structure and water content estimates from different resistivity models highlight a strong dependence on the starting model in the 3D inversions, but still enable us to put constraints on the deep structure of the craton. The present‐day thermal lithosphere‐asthenosphere boundary (LAB) depth is estimated to be at least 250 km beneath the Baoulé-Mossi domain. The area likely transitions from a cold and thick lithosphere with relatively low water content into thinner, more fertile lithosphere below the Volta Basin. Although the inferred amount of water could be explained by Paleoproterozoic subduction processes involved in the formation of the Baoulé-Mossi domain, later enrichment of the lithosphere cannot be excluded

Geometry of two glacial valleys in the northern Pyrenees estimated using gravity data

We estimate using gravity data the thickness of post-glacial unconsolidated sediment filling two major glacial valleys in northern Pyrenees: the Gave de Pau valley between Pierrefitte-Nestalas and Lourdes, and the Garonne valley between Saint-Beat and Barbazan. One hundred and eighty-four new gravity data complete 74 measurements obtained from the International Gravimetric Bureau database. Negative residual anomalies resulting from the presence of small-density unconsolidated sediment approach 4 mgal in both the Gave de Pau and the Garonne valleys. Estimating the sediment thickness requires knowing the density contrast between Quaternary sediments and the underlying bedrock. Supposing this density contrast is 600 kg/m(3), the maximum estimated thickness of post-glacial sediment is similar to 230 and 300 meters, and the volume of sediment is 2.1 and 3.2 km(3) in the Gave de Pau and Garonne valleys, respectively. In both valleys, the depth of Quaternary sediment suddenly increases at the confluence between two major glacial valleys (Gave de Pau - Gave de Cauterets, and Garonne - Pique confluences). Overdeepened basins are less deep downstream when approaching terminal moraines (Lourdes and Barbazan area), illustrating that the efficiency of glacial erosion depends on the ice flux flowing through valleys