Application of lithogeochemical and pyrite trace element data for the determination of vectors to ore in the Raja Au-Co prospect, northern Finland

. Discovering ore deposits is becoming increasingly difficult, and this is particularly true in areas of glaciated terrains. As a potential exploration tool for such terrains, we test the vectoring capacities of trace element and sulfur isotope characteristics of pyrite, combined with quantitative statistical methods of whole-rock geochemical datasets. Our target is the Rajapalot gold-cobalt project in northern Finland, where metamorphosed Paleoproterozoic volcanic and-sedimentary rocks of the Peräpohja belt host recently discovered gold prospects, which also have significant cobalt enrichment. The focus is 15 particularly put on a single gold-cobalt prospect, known as Raja, an excellent example of this unusual cobalt-enriched gold deposit, common in the metamorphosed terranes of northern Finland. The major lithologies at Rajapalot comprise amphibolite facies metamorphosed and polydeformed calcsilicate rocks that alternate with albitic units, mafic volcanic rocks, mica schist and quartzite. Mineralization at Rajapalot prospects is characterized by an older Co-mineralizing event and a younger high-grade Au-mineralization with re-mobilization and re-deposition of Co. Detailed in situ laser ablation inductively coupled 20 plasma-mass spectrometry (LA-ICP-MS) is a powerful technique that produces robust trace element and sulfur isotope databases from paragenetically and texturally well-characterized pyrite from the Raja prospect. The results are treated with appropriate log-ratio transformations and used for multivariate statistical data analysis, such as the computation of principal components. Application of these methods revealed that elements such as Co, Ni, Cu, Au, As, Ag, Mo, Bi, Te, Se, Sn, U, Tl and W have high vectoring capacities to discriminate between Co-only and Au-Co zones, as well as between mineralization 25 stages. The systematic pyrite study suggests that homogenous sulfur isotopic characteristics (+1.3 ‰ to +5.9 ‰.) and positive loadings of Co, Se, As, Te, Bi and Au onto PC1 are reflective of an early stage of Co-mineralization, while the opposing negative loadings of Mo, Ni, W, Tl, Cu and Ag along PC1 are associated with pyrites from the Au-mineralizing event. The sulfur isotopic signature of the latter pyrite type is between -1.2 ‰ and +7.4 ‰. Subtle patterns recognized from the whole-rock geochemistry favor an As-Au-Se-Te-W-U signature along the positive axis of PC1 for the localization of high-grade Au-30 Co-zones, whereas the element group Ni, Cu, Co, Te, Se and As, which has negative loadings onto PC2, will predict Co-only zones.

This study shows the efficiency of trace element geochemistry in mineral exploration targeting, which has the capacity to define future targets by characterizing the metallogenic potential of a host rock, as well as distinguishing various stages of mineralization.

Introduction
Orogenic gold deposits in Paleoproterozoic belts in northern Finland are attractive exploration targets.Some have been explored not only for gold but because of the occurrences of economically important concentrations of base metals, especially cobalt.The latter style of epigenetic-hydrothermal gold deposits of metamorphic terranes are classified as "orogenic gold deposits with atypical metal associations" (Goldfarb et al., 2001;Eilu, 2015).An example of this kind of gold mineralization has recently been discovered in the Paleoproterozoic Peräpohja belt, within the Rajapalot project area owned by the Mawson Gold Ltd (centered at coordinates 408750.412294Eby 7369184.18986N of the EUREF_FIN_TM35FIN coordinate system).
Seven prospects make up the broader Rajapalot Au-Co project, and this study focuses on a single prospect, Raja.Similar Corich orogenic gold deposits also occur in the Central Lapland belt and in the Kuusamo belt in northern Finland, and some of them are also undergoing active mineral exploration.Therefore, better knowledge of formation conditions of this specific type of deposits, and, especially, establishment of geochemical vectors towards ore, supports the recognition of occurrences of these ores in greenfield exploration projects, as well as the recognition of new mineralized zones in already known deposits.
The orogenic gold deposits in northern Finland are characterized by intense hydrothermal alteration that is associated with the occurrence of significant amounts of sulfides.The main sulfides are pyrrhotite, pyrite, chalcopyrite and pentlandite, along with other minerals including cobaltite, arsenopyrite and ± uraninite.All of these minerals can concentrate trace elements in their crystal structure, which can be used to understand the controls on the partitioning behavior of elements, the stages of mineralization, as well as overprinting features and the characteristics of ore-forming fluids (Godel et al., 2007;Dare et al., 2012;Piña et al., 2013;King et al., 2014;Duran et al., 2015;Cook et al., 2016;Keith et al., 2018).Trace element distribution patterns of those sulfide minerals not only provide a proxy for ore-forming conditions but are becoming increasingly important as vectors to ore (Keith et al., 2018).Significant improvements in trace element analytical techniques by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) allow time and cost-effective in situ analysis of precious and base metals with low detection limits in various minerals including sulfides.Therefore, this method can be a routine tool even in the early stages of mineral exploration projects.In that regard, pyrite is particularly suitable for this type of analysis, considering its capability of hosting a variety of elements of economic interest (Au, Ag, Bi, Co, Cu, Mo, Ni, Sb, Se and Te) and its compositional adaption to changes in the physicochemical fluid parameters (Liu et al., 2018;Large et al., 2009;Ulrich et al., 2011;Deditius et al., 2014;Keith et al., 2018;Liu et al., 2018).In addition to trace element analysis, sulfur isotope data are of significant assistance during exploratory work.The well-established sulfur isotope fractionation in pyrite improves our understanding of ore deposits and serves as a valuable and supplemental vectoring technique to the trace elemental and lithogeochemical approach (Ohomoto and Rye, 1979;Barker et al., 2009;Voute et al., 2019).Several authors have shown that a combination of aforementioned analytical techniques also highly support refining of genetic models of ore deposits (Barker et al., 2009;Hodkiewicz et al., 2009;Ulrich et al., 2011;Liu et al., 2018;Mukherjee et al., 2019;Voute et al., 2019;Meng et al., 2020;Vasilopoulos et al., 2021).
In this study we show how lithogeochemical compositions and a robust database of paired multi-element and sulfur isotope in situ analyses in pyrite by LA-ICP-MS, can be used for the identification of ore-related alteration zones and for the development of vectors towards the ore.By applying computational intelligence on these compositional datasets, we were able to develop multi-element associations as new efficient variables for specific metal associations.In order to investigate and identify the hidden trends and clusters in our geochemical data, we used multivariate statistical techniques that aim for dimension reduction, but still capture the essential variability in the original dataset (Aitchison, 1984).Following the suggestions of Aitchison (1986), Reinmann et al. (2012Reinmann et al. ( , 2017) ) and Hron et al. (2020) we discuss the problems of a proper pre-processing of compositional data for appropriate multivariate statistical analysis, and show how to avoid severe misinterpretations when performing principal component analysis.

Regional geology
The bedrock of Finland is defined by an Archean basement (3.5-2.5 Ga), its Paleoproterozoic sedimentary-volcanic cover (2.5-1.9Ga) and the Svecofennian orogenic domain (1.93-1.8Ga; Hanski, 2015; see Fig. 1a).Archean crustal segments are attributed to the Kola and Karelian cratons and are separated by the Lapland Granulite Belt and the Belomorian terrain (Fig. 1a; Lahtinen et al., 2005).Throughout the Paleoproterozoic tectonic evolution (2.5-1.9Ga), the Archean cratons underwent several stages of intracontinental extension and rifting, as well as continental margin rifting, which resulted in the formation of deep-scale structures and shallow water-basins (Lahtinen et al., 2005).The depositional history of these Karelian basins coincides with the Great Oxygenation Event, which created favorable conditions for a pre-concentrations of metals (e.g.Co, Cu, Ni, Mo and Au) in sulfidic sediments, as well as the deposition of carbonaceous material (Melezhik et al. 2013;Lyons et al., 2014;Large et al., 2014Large et al., , 2015Large et al., , 2017;;Gregory et al., 2017;Johnson et al., 2017;Armstrong et al., 2018.).Köykkä et al. (2019) proposed five basin evolution stages (see Fig. 1c) and summarized the volcano-sedimentary successions and intrusions with the following generalized lithostratigraphy (from bottom to top): (i) mafic volcanics and minor conglomerates (ii) clastic sedimentary rocks, (iii) subaerial mafic volcanics, komatiites and carbonate rocks; (iv) greywackes, carbonaceous-bearing rocks and sulfur-rich pelitic rocks and (v) phyllites and greywackes.All these units were then deformed and metamorphosed during the Svecofennian Orogeny (Lahtinen et al., 2005;Fig. 1c).Such a depositional and orogenic evolution of Archean to Paleoproterozoic settings is worldwide referred to as Paleoproterozoic greenstone belts.Within the northern part of the Karelian domain, the Svecofennian orogenic gold deposits were formed in the early stages of the accretion of microcontinents between 1.92 and 1.86 Ga, which resulted in the formation of the Fennoscandia Plate, and the far field effect of the collision of Fennsocandia and Sarmatia in the SE (Svecobaltic orogeny) and Amazonia in the west (Nordic orogeny) between 1.85 and 1.79 Ga (Weihed et al., 2005;Lahtinen et al., 2005;Molnár et al., 2017Molnár et al., , 2018) ) ).The gold mineralization is mainly located in complexly folded thrust zones and hosted near the boundary between the oxidized and reduced strata, which have experienced several stages of deformation and alteration (Hanski, 2015).One of these Paloproterozoic rift-related basins is the Peräpohja belt (PB) in northern Finland, located between the Archean granitoid Pudasjärvi complex to the southeast and the Central Lapland granitoid complex (CLGC) to the north (Figs. 1b and 2a;Vanhanen et al., 2015;Nironen, 2017).The maximum depositional age for the PB is defined by the NE trending 2.44 Ga layered intrusions of the Torino-Näränkävaara belt (Iljina and Hanski et al., 2005), scattered along the northern boundaries of the Pudasjärvi and Lentua complexes (Fig. 1b; Ranta et al., 2015;Nironen, 2017).After emplacement, normal faulting of these intrusions led to partial erosion of igneous layers, onto which the lowermost and oldest units of the PB were deposited (Sompujärvi conglomerates and Runkaus volcanic sequence, Fig. 1c; Hanski et al, 2005;Nironen, 2017).In the western part of the belt, the youngest supracrustal metasediments are cut by 1879 ± 3 Ma monzonite intrusions, which constrain the minimum age of the PB at 1.88 Ga (Lehtinen et al., 2005;Hanski et al., 2005;Ranta et al., 2015;Nironen, 2017).
The tectonic evolution of the PB is characterized by a polyphase deformation history (between ~1.9 and 1.8 Ga) and increasing metamorphic conditions towards the northern parts from lower-greenschist to upper amphibolite facies and local migmatization along the northeastern marginal zone (Fig. 1c; Hanski et al., 2005;Lahtinen et al., 2005;Laajoki, 2005;Hölttä and Heilimo, 2017.).A detailed description of the Svecofennian tectonic evolution of the PB and related emplacement periods of granitic intrusions is provided by Lahtinen et al. (2015) and Nironen (2017), where the authors discuss the five deformation stages that affected the PB in great detail (D1-D5; see Fig. 1c).

Geology, stratigraphy and mineralization of the Rajapalot Au-Co project
Several occurrences of gold mineralization are located within the northern domain of the PB near the contact with the Central Lapland Granitoid Complex (CLGC) and comprises a group of prospects, all owned by Mawson Gold Ltd. (Fig. 2a;Cook and Hudson, 2018).Intense exploration in this area started in the Rompas prospect where native gold bearing pockets and veinlets were recognized in uraninite-bearing calc-silicate veins (Vanhanen et al., 2015;Molnár et al., 2016a).Geochemical anomalies in boulders and till also indicate the presence of gold-bearing mineralization in the Rajapalot area, about 8 kilometers east of Rompas.At the current stage of the intense drilling program, six major mineralized zones have been recognized in this latter area and one of them is the Raja prospect (Fig. 2b).Inferred mineral resource-estimate at Raja only are at 2.97 million tonnes @ 2.9 g/t gold and covers average 383 ppm cobalt, which form 46 % of the gold-cobalt resources at Rajapalot (Mawson Gold Ltd., 2021).
Rocks in the area of the Raja prospect are multiply deformed and metamorphosed to amphibolite-facies, which makes the determination of protolithic features particularly difficult.According to Cook and Hudson (2018) the Au-Co mineralization is 145 hosted by two isoclinally folded sequences (Fig. 2c): (i) Sequence 1 comprises siliciclastic, albitized and carbonatized, largely oxidized metasedimentary rocks from a continental margin setting; (ii) Sequence 2 represents a metamorphosed sedimentary sequence formed under reduced conditions consisting of pelitic turbidites, sandstones, carbonates and sulfidic carbonaceous rocks.Mafic rocks (lava flows, dykes and volcanoclastic sediments) are common within both sequences.In the mineralized zones, domains of retrograde alteration to chlorite or epidote occur.150 The sampled drill cores are from two profiles (profile 1 and profile 2; Fig. 2b-c) from Sequence 2 of the Raja prospect.Major lithologies are mafic metavolcanic rocks (e.g.pillow basalt to amphibolite), albite-calcsilicate rocks, biotite-calcsilicate rocks, albitite, mica schist and muscovite-bearing quartzite (Fig. 2c).
The zones of high-grade Au-Co mineralization are characterized by sulfide disseminations adjacent to linear, or sub-linear near-vertical structures (Cook and Hudson, 2018;Cook et al., 2020).The Raja Au-Co resource extends 240 meters parallel to strike, 950 meters down plunge reaching a vertical depth of 560 meters (Mawson Gold Ltd., 2020).The most significant host lithologies are (brecciated) albitite and foliated muscovite-biotite schist within the lower drill core intercepts (sulfidic K-Mg-Fe-rich rocks, see Fig. 2c).The best grades are detected at the boundary between these two units, indicating that rheological differences played an important role in focused fluid flow and metal precipitation within the system (Farajewicz & Cook, 2021).
Given the low K-concentrations, this type of mineralization is also referred to as the sodic-type (Farajewicz & Cook, 2021).
Mica schist-hosted mineralization in contrast, is of the K±Fe-type, and occurs in strongly deformed, medium-to coarse-grained metapelitic rocks, which consist of quartz and albite porphyroclasts enveloped by mica-rich bands (muscovite/sericite and biotite) that include chlorite as a replacement product of biotite (Farajewicz & Cook, 2021).The most common ore minerals are pyrrhotite, pyrite, chalcopyrite, Bi-Te-rich phases, molybdenite, native gold, cobaltite and scheelite.
The structurally-controlled gold mineralization formed subsequent to the peak of amphibolite facies metamorphism (at around 1.78 Ga; Molnár et al.;2017) coinciding with the D5-deformation event of the Peräpohja belt (see Fig. 1c).This age was also confirmed by Ranta et al. (2017Ranta et al. ( , 2018)), who conducted a trace element, boron isotope and fluid inclusion study on tourmaline from quartz-tourmaline-sulfide-native gold veins at Rajapalot (Palokas prospect).By using tourmaline, the authors could evaluate the physico-chemical characteristics of the mineralizing fluid, which provide a genetic link to the 1.78 Ga granitoid magmatism and the fracture-related gold mineralization (Ranta et. al, 2017).This is consistent with the results from Re-Os isotope dating on molybdenite (1.78 ± 0.01 Ga) associated with gold in tourmaline-rich zones at Palokas prospect (Molnár et al., 2017).
3 Analytical methods and data processing

Sampling
For the purpose of this study, 83 samples were collected from two profile sections of drill cores within the Raja Au-Co prospect.
Sampling was oriented according to drill core logs and geochemical data (e.g.trace element concentrations) available in the Mawson Gold Ltd. database at the time of sampling.Sampled drill cores include PAL0104, PAL0092, PAL0093, PAL0097 and PAL0100, as well as PAL0088 and PAL0116 which intersected the zones of major mineralization and areas farther away from them (see Fig. 2c).The samples collected represents the major stratigraphic units present within the Raja prospect, as well as highly mineralized Co-Au zones and barren domains.

Analytical methods
Whole-rock geochemical analysis was performed by MSALABS on behalf of Mawson Gold Ltd using a 4-acid digestion and ICP-AES finish.The methodology used is according to the official procedure published by the MSALABS analysis sheet: Drill core samples are ground to meet the passing criteria of 85% -75μm and are then weighed and digested using hydrofluoric acid, hydrochloric acid, nitric acid and perchloric acid.After the digestion process, samples were treated with deionized water and analyzed by Inductively Coupled Plasma-Atomic Emission Spectroscopy and Inductively Coupled Plasma Mass Spectrometry.The reporting units are in ppm.To ensure quality, reference materials, blanks and duplicates were used during analysis; spectral inter-element interferences were corrected.The dataset comprises major and minor trace elements.
Mineralogical and textural observations were performed on polished thin sections, using transmitted and reflected light polarizing microscopes, as well as back-scattered electron imaging with a high-resolution scanning electron microscope (JEOL JSM-5900LV) at the Geological Survey of Finland (Espoo, Finland).This process built the foundation for understanding the mineralization processes and was significant in establishing a systematic, time-efficient plan for trace element and sulfur isotopic analysis.
Trace element analyses in pyrite were performed by a laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) at the Geological Survey of Finland (GTK), using a Nu AttoM single collector (SC)-ICPMS (Nu Instruments Ltd., Wrexham, UK) and an Analyte Excite 193 ArF laser-ablation system (Photon Machines, San Diego, USA).The laser was run at a pulse frequency of 10 Hz and a pulse energy of 5 mJ at 30% attenuation to produce an energy flux of 2.5 J/cm 2 on the sample surface with a 50 µm spot size.Each analysis was initiated with a 20 second baseline measurement followed by switching on the laser for 40 seconds for signal acquisition.Analyses were made using time resolved analysis (TRA) with continuous acquisition of data for each set of points (2 standards, 15 unknowns, 1 quality control standard).Synthetic pressed nanopellets of MASS-1 and UQAC standards were used for external standardization and the isotope 57 Fe as the internal standard.Measurements were performed over 68 isotopes and 36 elements at low resolution (∆M/M = 300) using the fast scanning mode.Data reduction was handled using the software GLITTER TM (Van Achterbergh et al, 2001) which allows baseline subtraction, integration of the signal over a selected time window and quantification using known concentrations of the external and internal standards.Analyses with large inclusions were discarded, but in the cases where inclusions were ablated towards the end of the spot analysis, only mixed signal portions were excluded.For a closer estimation of the actual mean, analyses with certain elements below the minimum detection limit (MDL) were treated as MDL/2 (Winderbaum et al., 2012).The following suite of isotopes was measured during spot analysis: 27 Al, 29 Si, 33 S, 34 S, 45 Sc, 48 Ti, 49 Ti, 51 V, 52 Cr, 53  Sulfur isotope analyses of pyrite were performed using a Nu Plasma HR multicollector (MC)-ICPMS at the Geological Survey of Finland in Espoo together with a Photon Machine Analyte G2 laser system.Samples were ablated in He gas (gas flows = 0.4 and 0.1 l/min) within a HelEx ablation cell (Müller et al., 2009).S isotopes were analyzed at medium resolution.During the ablation the data were collected in static mode ( 32 S, 34 S).Single spot samples were ablated at a spatial resolution of 40 micrometers, using a fluence of 3.5 J/cm2 and at 3 Hz on thin sections.The total S signal obtained was between 0.5 and 3.0 V.Under these conditions, after a 20 s baseline, 50-60 s of ablation is needed to obtain an internal precision of 34 S/ 32 S ≤ ± 0.000005 (1 SE).One pyrite standard was used for external standard bracketing (PPP-1; Gilbert et al, JAAS, 2014) and the inhouse standards Py1 and Py2 were used for quality control, both previously measured by gas mass spectrometry.For Py1, with a δ 34 SCDT(‰) reference value of -0.6 ± 0.3 ‰ (1s) we measured an average value of -0.6 ± 0.5 ‰ (1s, n=26).For Py 2 with a δ 34 SCDT(‰) reference value of -0.3 ± 0.3 ‰ (1s), we measured an average value of -0.2 ± 0.2 ‰ (1s, n=26).One in-house chalcopyrite standard (Cpy1) was used for external bracketing, while another was used for quality control.The in-house chalcopyrite standard Cpy 2, previously measured by gas mass spectrometry, has a δ 34 SCDT(‰) reference value of -0.7 ± 0.5 ‰ (1s) for which we measured an average value of -0.9 ± 0.3 ‰ (1s, n=18).

Data processing
Whole-rock geochemical and mineral trace element data are classical examples of compositional data, which consist of random vectors with components that are parts or percentages of a whole.These vectors only provide relative information and their components are constrained to sum to unity (Aitchison, 1983(Aitchison, , 1986)).Thus, when using geochemical data for multivariate analysis, one must be aware of the radical difference between the real Euclidean space associated with unconstrained data and the sample space for compositional data, as well as of the severe misinterpretations that can result from improper data preprocessing (Aitchison andGreenacre, 2002, Filzmoser et al., 2010).In order to filter out the relative and relevant information of analyzed elements, and to change the geometry of original data to real coordinates in the Euclidean space, compositional data are treated with the centered log-ratio (clr) methodology, which calculates the logarithm of the ratio between the element and the geometric mean of all other parts of the composition for each single sample (Aitchison and Greenacre, 2002;Reimann et al., 2008;Reimann et al., 2017, Hron et al., 2020).Following this approach, we avoid using the log-normal transformation and are preferring the clr transformation as a data pre-processing method for statistical data analysis, such as principal component analysis.

Principal component analysis
The concept of principal component analysis (PCA) is based on the dimension reduction of a large number of possibly correlated variables to a significant set of linearly uncorrelated variables, which are used to obtain element associations of ores (Macheyeki, et al., 2020).This correlation-based multivariate analysis technique models a complex set of variables in simplified patterns, but still accounts for most of the data variance (Reimann et al, 2008).Choosing the number of variables from a log-ratio transformed compositional dataset is crucial when expecting stable results from a PCA (see Le Maitre, 1982).
In general terms, the computation of principal component analysis (PCA) consists of the calculation of eigenvectors and the corresponding eigenvalues with respect to a covariance matrix (see Eq.1): with entries (see Eq. 2): for a given  ×  data matrix  = ( i,j ) i=1,…,n;j=1,…,p where  equals the number of observations and  the number of variables.Mathematically speaking, each row of the data matrix can be interpreted as the realization of a  dimensional random vector ( 1 , … ,   ), while the space that contains all possible realizations of this random vector is referred to as feature space.This  dimensional space will be denoted as   .
By selecting a subset of the eigenvectors   ,  = 1, … , , which result from the PCA procedure, one can define the axes of a subspace   ⊆   , where  < .In other words, a PCA enables us to express the variability of our data in terms of a lower dimensional subspace   .The axes of   are obtained by selecting a subset of eigenvectors from  1 , …,   .In order to determine this subset, each eigenvalue  1 ,…,   of the corresponding eigenvectors need to be taken into consideration.The objective is to achieve the maximum possible variance with respect to   .This is achieved by the eigenvectors which correspond to large eigenvalues.
In order to prevent a dominance of variables with a large absolute magnitude, a correlation matrix PCA, which is a PCA on the standardized data, was computed (Jolliffe & Cadima, 2016).In this context, the initial matrix  = ( , ), is replaced with the standardized data matrix  = ( , ), where  = 1, … ,  and  = 1, … ,  and (see Eq. 3): Each data value  , is centered and then divided by the standard deviation   of the  observations of variable  (Joliffe & Cadima, 2006).After this standardization, the data matrix  is used for the computation of the correlation matrix PCA.In order to assure that most of the data variability is preserved, only principal components with large eigenvalues (> 1) are considered (see Kaiser, 1960;Joliffe, 2002).
Biplots were produced as a graphical representation of the correlation matrix PCA.They describe the elements of the eigenvectors   (PC loadings) and the values that each observation would score on a particular principal component (PC scores; see Jolliffe & Cadima, 2016).Computation of PC was performed using the software ioGAS™ 7.2.Before entering the PCA, the datasets were treated with the appropriate clr transformation.Out of the 50 elements analyzed within the whole-rock geochemical dataset, a subset of 12 elements was used (As, Au, Co, Cu, Mo, Ni, Se, Sn, Te, Tl, U and W).From the 61 elements analyzed in pyrites, a total of 12 were used to produce geochemical groupings with the PCA (Ag, As, Au, Bi, Co, Cu, Mo, Ni, Se, Te, Tl and W).Resulting geochemical patterns provide further insights into mineralization processes and have the capacity to differentiate between mineralized and barren zones.

Molar element ratio analysis
The geochemical data were further used to determine the effects of hydrothermal alteration on the rock types and the mineralization (Madeisky and Stanley, 1993;Stanley and Madeisky, 1996;Stanley, 1998Stanley, , 2017)).Molar element ratio (MER) analysis and corresponding plots of 3K/(Al + Mg + (Fe − S) versus Al/(Al + Mg + (Fe − S), and (Na + K)/Al versus Na/Al, are the variables and ratios that best describe the observed alteration assemblages.For this purpose, the reported concentration of major elements that is stated in weight percent must be recalculated to a 100% volatile-free composition and converted to moles, in order to be used as ratios for the above-mentioned MER plots.

Description of major rock units
Given the intense metamorphism and multi-stage hydrothermal alteration history at Raja, determination of primary textures and mineral phases is quite complicated.The most common alteration types include albitization (pre-post-syn-orogenic), retrograde alteration assemblages of chloritized biotite, chlorite and Fe-Mg amphiboles (syn-orogenic), sulfidation-induced formation of massive pyrrhotite and subordinate pyrite and chalcopyrite (syn-to post-orogenic), as well as sections where plagioclase is affected by the sericitic alteration characterize the latest alteration event (post-orogenic; see Table 1).Based on the abundance of rock-forming minerals, major rock-types can be classified as albitite, albite-calcsilicate rocks, biotitecalcsilicate rock, metavolcanic rocks, mica schists and quartzites (see Fig. 3).The calcsilicate rocks define fine-to mediumgrained, folded and crenulated lithologies that display a color range from grey to dark grey (biotite-calcsilicate rocks; Fig. 3a), and light grey to light pink (albite-calcsilicate rocks; Fig. 3b).They consist of albite, quartz, amphibole, biotite and pyroxene (± carbonates, ± chlorite; ± hematite, ± magnetite) and can exhibit alternating bands of felsic and mafic mineral phases.
Disseminated pyrrhotite is the most common sulfide mineral in these rocks, followed by pyrite and chalcopyrite (Figs.3a-b).
Fine-to medium-grained albitite is composed of an albite-quartz-rich matrix with a light-greyish appearance.The color, however, displays wide variations raging from nearly white to pink.Sheet silicates such as biotite, muscovite/sericite and chlorite (retrograde phase) occur as interstitial phases or in veinlets and fractures (Fig. 3e).Brecciated and fractured albitite hosts the mineralization, commonly along the boundary with mica schist at lower stratigraphic positions.Pyrrhotite, pyrite and chalcopyrite are the most common sulfides; euhedral cobaltite is the dominant Co-bearing phase.
The fine-to medium-grained mica schist is a light grey-to grey-colored, foliated to crenulated muscovite and biotite-rich metapelitic rock (sulfidic K-Fe-Mg-rich rock) that also hosts the Raja mineralization, along with albitite (Figs.3e-f).This rock commonly contains quartz and albite porphyroclasts with enveloping muscovite and biotite, which often exhibit replacement of biotite by chlorite in sheared zones.Accessory phases are represented by magnetite, apatite, rutile ± tourmaline.Sulfide mineralization is dominated by foliated pyrrhotite ± pyrite, while the most common ore minerals are cobaltite, native gold, molybdenite, scheelite and Bi-Te-rich phases.
Light-grey quartzites are usually fine to medium-grained and display weak bedding textures and contain muscovite bands, as well as sericitized domains (Fig. 3g).

Ore textures and mineralogy
Ore minerals either occur in foliation-related dissemination or in veins that cut across foliation.The most common sulfide minerals are pyrrhotite, pyrite, chalcopyrite, pentlandite, molybdenite, galena and sphalerite; oxides are mainly ilmenite, magnetite and rutile.As the most abundant sulfide, pyrrhotite is accumulated in biotite-rich bands that follow the rock fabric.
In this textural setting pyrrhotite is associated with chalcopyrite, cobaltite and rare scheelite.Pyrrhotite also surrounds foliationrelated pyrite crystals in a fringe-like habit.The second most abundant sulfide at Raja is pyrite.It displays three distinct textural types, all of which are locally present in some samples (Fig. 4).They show no restriction to certain lithologies or specific stratigraphic levels.Pyrite A occurs as irregularly deformed or elongated grains following the foliation of the host rock (Fig. 4a) and is often associated with pyrrhotite that grows in antitaxial fringe-like shapes around the crystals of pyrite (Figs.4b-c).
Pyrite B is also observed along the foliation fabric, however, in a porphyroclastic texture with up to 1 cm large grains.These grains often have magnetite-amphibole-chlorite fringes (Fig. 4a and Fig. 4d-e).Less commonly, pyrite B also forms idiomorphic, coarse-grained crystals (up to 1.5 cm) without visible pressure shadows in the plane of thin sections.350 Occasionally, pyrite B hosts rare remnants of colloform pyrite, which shows concentric banding with uni-directional growth as a gradual filling of cavities (Fig. 4d).The anhedral pyrite C is associated with pyrrhotite and chalcopyrite (Fig. 4f) within chlorite and quartz-carbonate veinlets which cut across the foliation.Pentlandite and molybdenite are accessory sulfide minerals.Pentlandite occurs as flame-like exsolution lamellae (<20µm) or as anhedral inclusions in pyrrhotite.Molybdenite forms acicular crystals (50 µm to 250 µm) that are locally interlayered with Bi-Te-rich phases along rims of pyrrhotite (Fig. 5a).The Bi-Te-rich phases also occur along grain contacts with cobaltite (Fig. 5b).Given the small grain size of the Bi-Te-rich phases, no useful analyses could be produced to define their respective mineral names.The major Co-bearing phases at Raja are cobaltite, linnaeite and (±) Co-pentlandite.While linnaeite and Co-pentlandite mostly form as exsolutions in pyrrhotite (Farajewicz and Cook, 2021), the more abundant cobaltite is observed in two principal textural settings: (i) cobaltite I occurs in biotite-amphibole-rich domains with retrograde chlorite alteration.The euhedral (Fig. 5b) to anhedral (Figs.5c-d) crystals of cobaltite I cover a size range between 30 and 300 µm and are either enclosed in or occur along rims of rock-foliation-related pyrrhotite, suggesting cobaltite formation prior to deformation.Adding to this indication, is the presence of foliation-related biotite, partially replaced by chlorite.(ii) Anhedral cobaltite II grains (with up to 200 µm sizes) are embedded in biotite or muscovite and hosted by mica schist (Fig. 5e).Frequently, native gold and Bi-Te-rich phases are observed along the rims of cobaltite II (Fig. 5e).Apart from the close association with cobaltite II (Fig. 5e), native gold is also found as an interstitial phase in mica-rich domains (Fig. 5f).
Considering that both cobaltite types follow the rock-foliation within an anhedral shape, a formation prior to the Au-bearing stage is suggested.Whether the euhedral cobaltite I crystals (Figs.5b) and the anhedral cobaltite II (Fig. 5e) are both 385 representatives of the same Co-mineralization event or are attributed to two different Co-stages is difficult to distinguish.
Nevertheless, the overall textural relationships clearly suggest that native gold, Bi-Te-rich phases and molybdenite have grown after the formation of the Co-enriched phases.Table 1 puts the geological processes, alteration, as well as mineralization into context and provides a summary of the paragenetic succession of ore mineral precipitation at Raja.

Detection of ore zones
According to the whole-rock geochemical data from the sampled drill cores two significant types of ore zones can be distinguished at the Raja prospect: (i) Co-only zones are characterized by horizons with Co-enrichments of up to 1085 ppm and gold concentrations below the detection limit (b.d.l.); these horizons are primarily hosted by (biotite) calcsilicate rocks and amphibole-rich rocks in the stratigraphically high position in the profile defined by the drill cores (Fig. 6a and Fig. 7); (ii) Au-Co-zones are preferentially located at greater depths and hosted by albite-rich rocks and mica schist(Fig.6b and Fig. 7).
The Au-Co zones have up to 115 ppm gold and 2340 ppm cobalt (Fig. 6 and Fig. 7).Above the main mineralized position, a second type of Au-Co zone occurs (up to 2.67 ppm Au and 160 ppm Co; Fig. 6b) within metavolcanic rocks and albitecalcsilicate rocks (Fig. 6b and Fig. 7).The remaining barren areas have Co-concentrations below the estimated rock-average for bedrock units of the Karelian domain (see Rasilainen et al., 2008) and gold below the limit of detection.

Principal component analysis of whole-rock trace element datarecognition of ore-related and barren element associations
In order to substantiate meaningful geochemical patterns within mineralized and barren zones, multivariate statistical analyses, such as principal component analysis are applied to our large multi-element data set (see next section).To achieve stable results, it is recommended to enter the principial component analysis with variables selected based on geochemical reasoning, and not with the full set of analyzed elements (see Le Maitre, 1982;Reimann et al., 2008).For this purpose, metals, pathfinder elements and elements of economic interest that cover broad compositional fields within Co-Only, Au-Co and barren zones are particularly suitable.The boxplot in Figure 8 shows the visual distributions for the log-transformed concentrations of the key elements (As, Au, Co, Cu, Mo, Ni, Se, Sn, Te, Tl, U and W), which will enter the principal component analysis as clrtransformed variables (see next section).Within this set of variables, elements of economic interest (such as Au) and elements with a long-term recovery potential (such as Co, Se and W) are included to better understand their behavior and establish element associations.The selected variables for the  = 1084 dataset are displayed in biplots that illustrate the new coordinate system of principal components (see Fig. 9).
How much each observation scores on a particular principal component is described by PC scores, while the elements of an eigenvector define the PC loadings (or weights; Reimann et al., 2008;Jolliffe & Cadima, 2016;Holland, 2019).Positive loadings indicate a positive correlation with the principal component, while negative loadings point to a negative correlation.
The vector-length is proportional to the variability of the two presented principal components, while the direction and angle of the vectors defines the correlation between the variables.
Table 2 summarizes the contribution of each variable to the first four principal components, whose explained variation is above 1 (see eigenvalues), which accounts for almost 70 % of the total variance.(see % explained variance).This is sufficient to extract the following relevant information contained in the correlation matrix PCA (see Table 2): The largest positive loadings onto PC1 has As, followed by Mo, W, Se, Au and Te.Elements such as Sn, Tl and Ni on the contrary, display the largest negative loadings on PC1, indicating lower concentrations of these elements on the positive axis 1. Uranium contributes to the PC2 with large positive loadings, along with Sn, W, Tl, Mo and Au, whereas Co and Cu display large negative loadings, which corresponds well with the negative to almost no correlation between the Co-Cu-group and the Au-Mo-W-U-group and the Sn-Tl-association, respectively (see Fig. 9a).PC3 is dominated by large positive loadings for Au, W and Cu, as well as large negative loadings for Se and U, as well as negative loadings for Mo, Co and Ni (Fig. 9b).Variables that strongly contribute to PC4 with positive loadings are Te, Au, Ni and Se; the largest negative loadings are displayed by Mo (Fig. 9c).The biplots in Figure 9 capture the loadings of the twelve variables in relation to the first four principal components and shows the separation of samples in Au-Co zones, Co-only zones and barren zones within the new coordinate system.In the PC1-PC2 projection (46.9% of the total variance; Table 2) the majority of datapoints from Au-Co zones is concentrated in the positive domain of PC1 and distributed between the loadings of U, W, Mo, Au, As, Te and Se, where it shows overlapping features with datapoints from barren zones and Co-only zones within the positive and negative PC2 space, respectively (Fig. 9a).Along the negative axis of PC1, between the negative loadings of Ni, Sn and Tl, overlapping features are observed for samples from Co-only and barren zones (Fig. 9a).A separation between these two zones, however, is better visualized with the PC2, considering that almost all samples from barren zones are distributed along the positive axis of PC2, where the large positive loading for U in the barren zones is opposing the large negative loadings for the Co-Cu-group in the negative PC2 space of the Co-only zones.Given this observation, PC2 appears to be the principial component that is particularly suitable to distinguish between mineralized and barren zones.
Although the first two principal components account for most of the data variability (46.9%) and provide by far the largest capacities to differentiate between Co-only zones, Au-Co zones and barren sections, it is worth discussing the geochemical associations within the PC1-PC3 projection, which represents a total of 37.9% of the variance (Table 2; Fig. 9b).The dominant characteristic in this projection is that the PC3 is also capable to separate the Au-Co zones, which are distributed along the positive axis of PC3 (and PC1), in the domain of the positive loadings of Au, W, Cu and As, from the Co-only zones, which dominate the negative axis of both PC3 and PC1. (Fig. 9b).470 In the PC1-PC4 projection (36.8% of the variance; Table 2) the dominant positive loadings on PC4 are Te and Au, which are opposing the negative loadings of Mo, As, Cu and Tl along PC4 (Fig. 9c).While the majority of observations from Au-Co zones are constrained within the positive PC1 and along the positive and negative axis of PC4, samples from Co-only and barren zones show overlapping features along the positive and negative axis of both principal components (see Fig. 9c.).Considering that the PC1-PC2 projection has the largest capacity to distinguish between mineralized and barren zones, it is used to unravel underlying geochemical patterns by correlating the concentrations of elements with the distribution of samples in the loadings plot (see Figs. 9d-o).By doing so, the following observations are made: (i) Along the positive axis of PC1, a group of seven elements (As, Au, Mo, Se, Te, U and W; Figs.9d-j) describes the high-grade fracture-related Au-Co zones from albite-rich rocks and mica schists, while the few datapoints along the negative axes of PC1 (and PC2), between the loadings of Ni, Cu and Co, are reflective of the stratigraphically higher positioned Au-Co zones from metavolcanic units with lower Au and Co contents (Figs.9k-m).
(ii) The majority of Co-only observations is concentrated within the negative PC2 space and along the negative and positive axis of PC1, where two groups can be distinguished: (a) high-grade Co-only systems, positioned stratigraphically above the high-grade Au-Co zones, are distributed between the positive loadings of As, Te and Se on PC1 (Fig. 9d; Figs9g-h; Fig. 9l) and (b) Co-only zone from stratigraphically higher positioned metavolcanic rocks with lower Co-concentrations within the negative PC1-PC2 space (Fig. 9k-m).Lower concentrations of Co are also detected in the few samples from Co-only zones that plot along the positive axes of PC2, hosted by metavolcanic rocks well above and well below the high-grade Au mineralization (Fig. 9l).
(iii) The positive association of Sn and Tl (Fig. 9n-o), as well as Sn and U (Fig. 9i) in the positive PC2 space best applies to barren rocks, which usually show an inverse positioning of most base metals.

Relationships between whole-rock major element geochemistry and mineralization
The evaluation of whole rock major oxide data together with the trace element data reveals useful and important information about the lithological control on the mineralization, the mechanisms of ore minerals precipitation, as well as the effect of hydrothermal alteration on the rock composition (Madeisky andStanley, 1993, Stanley andMadeisky, 1996;Stanley, 1998Stanley, , 2017;;Stanley et al., 1996).For this purpose, we are using molar element ratio analysis and corresponding plots to graphically represent mass transfer processes responsible for hydrothermal alteration.The data presented in Figure 10 show the major element geochemistry of collected drill core samples, plotted in molar element ratio diagrams of 3K/(Al + Mg + (Fe − S) versus Al/(Al + Mg + (Fe − S), and (Na + K)/Al versus Na/Al, which correspond to mineral compositions of the observed alteration minerals.The illustration of multiple hydrothermal alteration events is facilitated by labeling specific mineralogical assemblages according to different molar values from the whole-rock chemistry (Fig. 10).Geochemical gradients recognized from the correlation of alteration mineralogy with Au-and Co-grades in MER dimensions exhibit distinct trends (Figs.10a-b; Figs.10d-e): mica schist samples record the highest Au-grades (up to 33 ppm Au) and are connected to sericitization by showing a clear trend towards the muscovite node.Significant concentrations of Au are also detected in albitites (up to 4.96 ppm) and albite-calcsilicate rocks (up to 3.22 ppm) that plot close to the albite node or between the albite-muscovite nodes.
The biotite-calcsilicates that plot between the albite-chlorite domain (Fig. 10a), in addition to a pronounced retrograde chlorite alteration (see petrography section), exhibit a trend towards the muscovite node with up to 1.91 ppm Au.Samples that are trending towards the albite and sericite end member generally display higher Au concentrations (Figs.10a-b; Figs.10d-e).
The lowest Au-grades are detected in quartzite (up to 0.57 ppm) and mafic metavolcanic rocks (up to 0.6 ppm).Petrographic analysis suggests, the latter underwent retrograde chlorite alteration, consistent with the concentration of data towards the chlorite node (Figs.10a & 10d).The highest concentrations of Co are detected in albite-calcsilicate rocks (up to 2190 ppm) and albitite (up to 1445 ppm) between the albite-chlorite end members, as well as in biotite-calcsilicate rocks (up to 2150 ppm) and mineralized mica schists (up to 1880 ppm) when close to the muscovite node (Figs.10c &10f).When considering that albitization is a pre-mineralization event, lithologies that experienced this type of alteration, most likely generated a nano-scale porosity that provided potential pathways for the early Co-and late Au-transporting fluids.This scenario could explain the enrichments of Co, as well as Au in the competent albite-bearing rocks (Figs.10b-c; Figs.10e-f).Important to note is that the rocks that experienced the late sericitization (e.g., mica schists) do not only record the highest Au, but also have high Co-concentrations.

Trace elements
The three categories of texturally distinct pyrite types (A, B and C) are further distinguished by their respective trace element compositions.The descriptive statistics of trace element concentrations from each pyrite type are reported in Table 3 and provide an overview of minimum, median, geometric mean and maximum values of elements recorded within the pyrite structure.
Accordingly, rock-foliation-related pyrite A is enriched in Co (up to 3.5 %) and contains up to 4142 ppm As (Fig. 11).Further compositional characteristics are significant concentrations of Au (up to 4.38 ppm), Bi (up to 94.45 ppm), Se (up to 425.4 ppm) and Te (up to110.4ppm), but low concentrations of Ni (up to 138 ppm), Cu (up to 158 ppm), Mo (up to 0.94 ppm) and W (up to 0.01 ppm; see Table 3 and Fig. 11).Trace element concentrations of recrystallized pyrite B are irregular and cover a broad compositional field.The majority of cobalt data are distributed between 20.6 ppm and 3.5 %.As an effect of core to rim zoning, wide variations of cobalt, ranging between 2 ppm and 5900 ppm, are predominately detected in pyrite porphyroclasts.
No compositional zoning is reported from the late fracture-related pyrite C, which is also distinct from other types of pyrite due to elevated Ni-contents (up to 3800 ppm).Apart from nickel, this pyrite type is characterized by enrichments in Co (up to 8920 ppm) and the highest concentrations of Mo (up to 713 ppm), Au (up to 2.27 ppm), Cu (up to 1.14 %), Pb (up to 8179 ppm), Tl (up to 123 ppm) and W (up to 242 ppm; Fig. 11).Regarding Se/S ratios, pyrite A covers a relatively narrow range with the lowest mean value of 2.93 × 10 −4 , while pyrite B displays a much wider scattering of Se/S (from 4.00 × 10 −6 to 0.003) and a mean value of 4.12 × 10 −4 ; the Se/S ratio in pyrite C has an overall mean value of 3.60 × 10 −4 (see Table 3).In order to avoid unstable results or misinterpretations for the subsequent multivariate analysis, a set of variables that includes metals and pathfinder elements with distinct compositional differences between all three pyrite types is used.The boxplot in

Table 3: Summary of the trace element composition of pyrite A, B and C (in ppm), as well as results of the sulfur isotopic analyses (in ‰).
Figure 11 shows the 12 elements that will be entered into the PCA as clr-transformed variables (As, Co, Cu, Ni, Mo, Tl and W, Ag, Au, Se, Bi and Te).After performing the PCA from the  = 260 dataset, certain element groupings are observed in the biplots of Figure 12, which only include principal components with eigenvalues > 1 (see Table 4).The foliation-related generations of pyrite A and B are chiefly concentrated in the right part of the PC1-PC2 projection, which represents a total of 49.44% of the data variability.This domain is characterized by positive loadings of Te, Bi, Se, As, Au and Co onto PC1 (Fig. 12a; Table 4).In contrast, the late fracture-related pyrite C are scattered along the negative axis of PC1 on the left part of the plot, which highlights the loadings of Tl, W, Ag, Cu and Mo (Fig. 12a; Table 4).The disassociation between these two element groups builds the basis to distinguish between the early (A and B) and late generation pyrites (C) and is expressed as the new variable PC1.Regarding PC2, it is particularly noticeable that the observations from pyrite A form two distinct groups along the positive and negative axes, which correspond to the positive loadings of As, Se and Co, and to negative loadings of Au, Bi and ± Te, respectively (Fig. 12a).Additional elements with positive loadings onto PC2 are Mo, Ni and W, while Ag, Tl and Cu are restricted to the negative space of PC2.Although most of the variance is captured by the first two principal components, the PC1-PC3 projection still accounts for a total of 41.55% of the data variability and carries information which is worth being presented (Fig. 12b).Particularly notable in this plot is the clear distinction between pyrite  b).While the PC2-PC3 projection only accounts for 32.35% of the data variability, it is still capable to differentiate between pyrite A and B, as well as pyrite B and C, along the PC3 axis, with the large positive loadings of Co and W that are essentially opposing the negative loadings of the remaining elements (Fig. 12c; Table 4).
Focusing on the PC1-PC2 projection, which accounts for most of the data variability, it appears that the chemistry of pyrite A reveals much more than just being enriched in Co. Pyrite A is distributed within the Co-As-Se group and the Au-Bi group, which are separated by the loading of Te (Fig. 12d-l).By correlating the concentrations of variables with projected element associations in the loadings plots, it shown that Co and Se concentrations are evenly distributed within the two aforementioned domains (Figs. 12d-e), whereas Au and Bi, and are only enriched in grains within the Au-Bi group (Fig. 12f-g).Concentrations of As and Te are higher within the positive and negative PC2 space, respectively (Figs. 12h-i).Additional characteristics of pyrite A are the negative correlations between Co-Ni, Co-W, and Au-Tl and almost no correlation between Co-Mo.
In the case of pyrite B, where trace element concentrations fall within a broad range, an interpretation of results may be a difficult task.The PCA, however, helps in this regard and shows the real data structure and clustering.Pyrite B is a major host for Se (Fig. 12e) and contains significant concentrations of As (Fig12h), Te (Fig. 12i), Ni (Fig. 12j) and Cu (Fig. 12k).Another noteworthy observation is that high concentrations of Au, Bi, Te and Cu are only valid for observations that plot along the negative axis of PC2.
The fracture-related pyrite C essentially occupies the entire left domain of the PC1-PC2 projection, where observations are scattered between the loadings of Mo, Ni, W, Tl, Cu and Ag.Pyrite C is further enriched in all these redox sensitive elements (Figs.12j-o) and has a significant Au-concentrations when close to the loadings of Cu, Ag, as well as Tl and W (Fig. 12 f).

Sulfur isotope data
For most of the pyrite grains both trace element and sulfur isotopes were analyzed by placing the analysis spots next to each other, separated by <30µm, in order to determine possible correlations.Unfortunately, no significant trends are observed between the trace element and sulfur isotope data.The range of δ 34 SVCDTdata from pyrite is between -1.8 ‰ and +7.4 ‰ (Fig. 13a; Table 3).Pyrite A, which has the highest Co/Ni ratios, displays positive only δ 34 S signatures (+1.3 ‰ to +5.9 ‰), contrasting with pyrite B (-1.8 ‰ to +7.3 ‰) and pyrite C (-1.2 ‰ to +7.4 ‰; Fig. 13b; Table 3).The standard deviation of δ 34 SVCDTdata for pyrite A is generally smaller (1.3‰), compared to pyrite B (2.0‰), and pyrite C (2.3‰).Pyrite B records diffuse within-grain and within-sample variations of the sulfur isotope compositions, while the heterogeneity of δ 34 S values from pyrite C seems to be more lithology-dependent.Mafic metavolcanic and biotite-calcsilicate units, cover a range of -1.2 ‰ and +5.4 ‰, whereas pyrite C from albitite, and albite-calcsilicate rocks records measured δ 34 SVCDTdata from +5.3 to 7.4 ‰.

Relationship between alteration and concentration of elements in the host rocks
The area of the Rajapalot Au-Co project hosts an epigenetic hydrothermal Au-Co-system that has properties comparable to an orogenic gold deposit with atypical metal association.Considering that the supracrustal sequences of the PB formed in an old intracratonic basin, multiple hydrothermal events could have invoked metal precipitation during pre-, syn-and post-orogenic processes (Lahtinen et al., 2012;Eilu, 2015;Molnár et al., 2016Molnár et al., , 2017;;Ranta et al., 2018).There is even the possibility that prior to regional metamorphism the episodic intrusion of felsic and mafic magmas locally altered the supracrustal sequences and thus enabled a pre-metamorphic (uneconomic) enrichment of base metals (Eilu, 2015;Molnár et al, 2017;Ranta et al., 2020).Referring to the PB, rocks at Rajapalot were affected by extensive alteration summarized in three events: (i) Regional pre-orogenic or early-metamorphic albitization, which predates the peak deformation and gold mineralization, was caused by greenstone belt-scale hydrothermal circulation of magmatic and/or basinal brines (see Eilu et al., 2003), which prepared the country rocks for a possible preenrichment of metals and for later Au-mineralizing processes (Pankka and Vanhanen, 1992;Vanhanen, 2001).In order to discuss the pre-enrichment possibility of metals, we should consider the MER diagrams (Fig. 10), which correlate mineralization with different styles of hydrothermal alteration that have affected the rocks.What stands out immediately is the positive correlation of high Co-grades with albitization.This observation either indicates that these rocks were already enriched in cobalt during the basin stage (pre-Svecofennian orogeny) of structural evolution or suggests that extensive albitization had a major control on the Co-mineralization (e.g., metal transportation by saline basinal fluids), or else points to an overprint of albitized rocks by the later Co-mineralization, unrelated to the albitization event.However, it should be considered that albitization is a process that generates nano-scale porosity, which can provide pathways for early Co-transporting fluids.A further alteration that could have played a role in the Co-enrichment is the retrograde chloritization, which affected biotite-rich rocks, such as mafic metavolcanic or biotite-calcsilicate rocks.However, if retrograde chlorite alteration was the event the contributed to Co-enrichment, significant amounts of cobalt should be recorded by these rocks, which is not observed from the MER diagrams (<206 ppm). (ii) The late Au-mineralization is observed in albite-rich rocks, as well as mica schists, some of which are also enriched in Co.These rocks underwent several stages of alteration, such as albitization, chloritization and sericitization producing widely varying rheological hosts to mineralization (Farajawicz & Cook, 2021).These contrasting host rheologies will play a significant part in the design of suitable crushing circuits capable of releasing gold and cobalt into chemical or flotation treatment cells.More specifically, the albite-rich strata are more enriched in Co, while the muscoviterich units have higher Au concentrations (see Fig. 10).To clarify the succession of mineralizing events, the results from the PCA should be considered, which show the negative correlation of Au and Co, as well as their opposing loadings onto the first three PC (see Fig. 9 and Table 2).These features suggest that different processes could have generated the Co-only mineralization, which was later overprinted by the Au-mineralization, resulting in Co-only and Au-Co-zones, respectively.
Nevertheless, the crucial questions in the viewpoint of mineral exploration are whether there is an explanation for the

Sources and precipitation of cobalt and gold
Here we discuss the possible sources of metals enriched in the two types of ore zones that have been discovered at Raja: (i) The Co-only zones hosted by albitite, amphibolite and biotite-calcsilicate-rich rocks above the main ore body, and (ii) the Au-Co zones hosted by fractured albitite and foliated mica schist.
The presence of Co-only zones, as well as the abundance of Co-rich but Ni-poor pyrite (A) and cobaltite, both following the rock-foliation, further supports the model of an early cobalt mineralization that was later overprinted by the orogenic gold mineralization.In order to mobilize cobalt, a hydrothermal tetrahedral chloride-complex is required, since tetrahedral complexes display a higher affinity for chlorine than octahedral complexes (Brugger et al., 2016).Contrasting to NiCl2(aq) species, which favor higher temperature and higher salinity when transitioning from the octahedral to the tetrahedral forms, tetrahedral CoCl2(aq) species are the overall dominant species at high temperature in solutions with low to moderate salinity (Brugger at al., 2016). Brugger et al. (2016) thus conclude that the solubility of cobalt is more sensitive to salinity than nickel, which is why the mobility of cobalt in chlorine-rich hydrothermal fluids is higher than that of nickel.Considering the intratectonic basinal setting at Rajapalot, the source rocks for the Co-transporting high-salinity fluids, could have been the evaporitic Petäjäskoski Formation (see Fig. 2; Ranta et al., 2018).This is in good agreement with the extensive pre-orogenic albitization event, which is widespread in the PB and is typically initiated during sedimentation and diagenetic processes in evaporitic environments (Eilu et al., 2003;2015;Ranta et al., 2018).This pre-orogenic Na-metasomatism mechanically prepared the rocks for the Co-rich mineralization within the Co-Only zones.Postdating the albitization, retrograde chloritization resulted in the replacement of biotite, which is why chlorite is locally accompanying the earlier precipitated, foliation-related cobaltite, as well as Co-rich pyrite A. In addition, chlorite-amphibole-magnetite-fringes commonly surround the porphyroclastic pyrite B, suggesting that chloritization postdated the Co-mineralization, rather than heavily influencing it.
The gold-cobalt mineralization at Raja is hosted by mafic metavolcanic units and underlying albite-calcsilicate rocks and by the stratigraphically lower positioned fractured albitite and foliated mica schists.Deformation promoted the re-activation of major shear zones, along which circulation of metal-rich and Au-bearing hydrothermal fluids was possible (1.84-1.77Ga).
The potential sources for a large-scale fluid mobilization, capable of transporting metals to form an orogenic gold deposit, however, are the subject of discussion.In general, orogenic gold deposits, as defined by Groves et al. (1998), are vertically extensive deposits at crustal depths from 2 to 20 km in accretionary or collisional fore-arc or back-arc settings.The Autransporting agent is a low salinity, near neutral H2O-CO2-(± CH4 ± N2 ± H2S) fluid that possibly formed from metamorphic devolatilization of supracrustal rocks at mid-upper greenschist-to amphibolite-facies conditions (Powell et al., 1991;Phillips and Groves, 1983;Goldfarb et al., 1986Goldfarb et al., , 1988Goldfarb et al., , 1989;;Groves et al., 1987;Bierlein and Crowe, 2000;Goldfarb et al., 2001;2005;Eilu, 2015;Goldfarb and Groves, 2015;Groves et al., 2019;Patten et al., 2020).According to Phillips and Powell (2010), the volume of fluids released from these rocks is enough for the formation of giant deposits.The genetic concept of the supracrustal metamorphic model supported by aforementioned authors, emphasizes upwards advection of mineralizing low salinity fluids along fault networks.The migrating reductive fluids that carry gold as thiobisulfide or other sulfuric complexes, then deposit the gold into zones of low tensile strength, where it is precipitated from solution by reducing fluid-rock interaction processes, as well as sulfidation reactions (Phillips and Powell, 2010;Garofalo et al., 2014;Goldfarb et al., 2019).
Another process that may have influenced the hydrothermal activity that led to a high-grade gold accumulation at Raja, is the emplacement of post-orogenic granitic intrusions.However, genetically linking the post-tectonic intrusions with fracturerelated gold does not fit the supracrustal metamorphic model.In association with the tectonic evolution of the PB, three emplacement periods of granitic intrusions are distinguished: (i) pre-orogenic Kierovaara granite (1.99 Ga; Ranta et al, 2018), (ii) the syn-orogenic Haparanda series monzonitic intrusions (1.89-1.86Ga; Perttunen and Vaasjoki, 2001) and (iii) the emplacement of late-to post-orogenic pegmatitic (tourmaline-rich) granites between the 1.79 and 1.77 Ga (Ranta et al., 2015(Ranta et al., , 2018)).The age of the latest, post-tectonic intrusive event falls within the age of gold-mineralization, which was estimated by Molnár et al. (2017) to be around 1.78 Ga.Other studies from the Rompas-Rajapalot area (Molnár et al., 2016a;Ranta et al., 2017Ranta et al., , 2018)), also emphasize that the circulation of gold-bearing fluids along fracture systems could be connected with the post-Svecofennian thermal event (1.79-1.75Ga) indicated by the post-orogenic granitoid magmatism (1.79-1.77Ga).
At Raja, the PC1 carries enough information about the pyrite trace element chemistry, to distinguish between the early Comineralization (pyrite A and B) with positive loadings of Se, As, Co, Te, Bi and Au and the late Au-mineralization (pyrite C) with negative loadings of Mo, Ni, W, Tl, Cu and Ag (see Fig. 12, and Table 4).
The dominant characteristics between Co and Ni in the PC1-PC2 projection are their opposing loadings onto PC1 and their negative correlation (see Fig. 12 and Table4).The pyrite data cluster very well in this plot and concentrations of both elements are always present above their respective detection limits.Such a behavior is not observed from the whole-rock data, in which   1).
Another important feature of Raja pyrites is the behavior of arsenic, which shows a positive correlation with Co, Te and Se.
The Co-As association is in good agreement with the presence of foliation-related cobaltite and Co-enriched pyrite A, while the As, Te and Se association can be interpreted as an anion substitution for S, considering that all three elements can play the role of anions in pyrites (George et al., 2008).Despite the high concentrations of As (up to 7715 ppm), only a minority of pyrite grains at Raja contain Au concentrations above the detection limit with a maximum of 4.38 ppm.Deditius et al. (2014) have demonstrated that invisible Au and As generally correlate positively in pyrite when As contents are above 0.01-0.1 wt%, and that both elements favor precipitation at decreasing temperatures (between 500 and 200°C) in various type of gold deposits.
According to the trace element compositions of pyrite in our study, however, Au-concentrations are independent of the Ascontent: this is in agreement with presence of free gold in the Au-Co ore at Raja.In addition, the PC1-PC2 projection displays opposing As-and Au-loadings along the axis of PC2 (see Fig. 12a and Table 4).The PC3, however, points out that some grains could show an association between Au and As, when considering their negative loadings (see Fig. 12b).Particularly noteworthy is the positive correlation of Au and the Bi-Te association in all PC projections.Referring to petrography, this observation is consistent with the overall association of native Au with Bi-Te-rich phases in mineralized sections.The Au and Bi signature could further indicate the coupled substitution 2Fe 2+ ↔ Au + + Bi 3+ (George et al., 2008).A geochemical affinity between Au and Bi is common and has been reported from different deposits (Nurmi and Sorjonen-Ward, 1993;Skirrow and Walshe, 2002;George et al., 2008).↔ (Cu, Ag) + + Sb 3+ .Considering the remarkable positive correlations among these elements, their incorporation in the pyrite structure via substitution for Fe 2+ is likely.Alternatively, Tl, Ag and Cu could be incorporated into nanoparticles within pyrite, commonly observed in hydrothermal ore systems that host pyrites with a wide range of trace elements (Deditius et al., 2011).
Scenarios that could facilitate the formation of nanoparticles include a direct precipitation from hydrothermal fluids, or exsolution from the pyrite matrix (Palenik et al., 2004;Deditius et al., 2011).Although a careful data checking of LA-ICP-MS ablation profiles for possible ablated inclusions was carried out, the presence of nanoparticles containing portions of trace elements should not be excluded (Cook et al., 2009).In order to determine whether the relatively high concentration of all analyzed trace elements in Raja pyrites are due to nanoparticles would require detailed high-resolution transmission electron microscopic measurements, which however is beyond the scope of this study.

Signatures of sulfur isotopes in pyrite and ore-forming processes
The overall δ 34 SVCDT data of pyrite grains at Raja are between -1.8 ‰ and +7.4 ‰, covering a common range of sulfur isotope composition for orogenic gold deposits (Molnár et al., 2016a).Considering the similarities of δ 34 SVCDT values in each pyrite type, these signatures are not diagnostic, but rather indicative of anomalous Co-or Au-zones.However, the combined use of trace element geochemistry and S-isotopic signatures in pyrite can help to interpret the behavior of S, hydrothermal processes, fluid/rock interactions and accumulation of metals during the formation of the Raja Au-Co prospect.
Formation of foliation-related pyrite A was facilitated by Co-transporting hydrothermal fluids in a reduced environment.More specifically, the δ 34 SVCDTvalues of the early Co-enriched pyrite A type only record positive results between +1.3 ‰ to +5.9 ‰.This range indicates a largely homogenous reduced source of sulfur, which is either related to a reduced magmatic and sedimentary environment, or to a reducing fluid that was released during prograde metamorphic devolatilization of metasedimentary rocks (Seal, 2006;Molnár et al., 2016a).This is probably also reflected in the relatively homogenous trace element characteristics of the majority of pyrite A, as well as in the homogenous Se/S ratios (see Table 3).Eventually these fluids could have originated from the evaporitic Petäsjäskoski Formation and caused extensive albitization, as well as a significant pre-enrichment of cobalt.Considering that significant concentrations of Au are detected in this pyrite type, we further propose an early, local deposition of gold.
The significant changes in Co/Ni ratios in pyrite B, and the wider range in sulfur isotopic values (-1.8 ‰ to +7.3 ‰), as well as the erratic within-grain variations of δ 34 SVCDT (up to 4‰.), may be caused by the variation in the redox conditions in the fluid system (Ohmoto, 1972;Hodkiewicz et al. 2009;Molnár et al, 2016b).Such a variation is not uncommon in magmatic hydrothermal systems with predominantly magmatic origin of sulfur.However, mixing of sulfur from different sources (e.g., sulfur isotopic data could provide a link to the latest, post-tectonic intrusive event at PB, that not only falls within the age of gold-mineralization (1.78 Ga; Molnár et al., 2017;Ranta et al., 2018), but also could have induced magmatic fluid circulation that was potentially involved in the formation of some of orogenic gold deposits (Dubé and Gosselin, 2007;Duuring et al., 2007;Kitney et al., 2011;Molnár et al., 2016Molnár et al., , 2017)).

Conclusion and implications for mineral exploration
The purpose of this study was to establish new synthetic geochemical variables by combining the whole-rock geochemistry with a robust database on trace element measurements of different pyrite types in an orogenic gold-base metal setting.We conclude that these results not only have the capacity to define the metallogenic potential of the host rock but also are capable of distinguishing various stages of mineralization.Even though this method requires a sophisticated laboratory background with a LA-ICP-MS analytical infrastructure, sampling, as well as petrographic observation prior to analysis, acquiring the laser based analytical data is a relatively quick process that generates data significantly faster than the whole rock geochemical approach.By compiling a large geochemical database for pyrite and other sulfides from different orogenic Au-systems in northern Finland (Raič et al., in prep), discrimination patterns can be established, which will help to define the metallogenic context for target areas with more accuracy.Pairing of these results with the sulfur isotopic signatures of measured pyrites allows detection of varying redox conditions, which is crucial when trying to delineate future target areas.This study further shows that with appropriate log-ratio transformations, principal components can be computed, which provide an insight into the true structure and unravel subtle geochemical patterns associated with ore-forming processes.Results from these multivariate statistical data analyses may help to define future targets in the Rajapalot project, and potentially elsewhere in hydrothermal zones with Au-Co-(base metals) in northern Finland, and are summarized as follows: (i) Geochemical variables identified from the whole-rock geochemical dataset suggest that host lithologies of Au-Co zones are in favor of an As-Au-Se-Te-W-U signature along the positive axis of PC1.The negative PC2 represents the group of elements Ni, Cu, Co, Te, Se and As, which are a characteristic association in Co-only zones, whereas the Sn-Tl group along the positive axis of PC2 will point to barren zones.

Figure 1 :
Figure 1: (a) Geologic map of the main bedrock units of the eastern part of the Fennoscandian Shield and the Peräpohja belt as highlighted area (black rectangle).Simplified geological map is based on Koistinen et al. (2001) (b) Geologic map of northern Finland showing Archean and Proterozoic rocks.The Rajapalot Au-Co project is located within the shaded rectangle.Simplified map is modified from the Bedrock of Finland -DigiKP (2021; digital map database of the Geological Survey of Finland).(c) Lithostratigraphy and structural evolution of the Peräpohja belt modified after Köykkä et al., (2019) and Lahtinen et al., (2015), respectively.

Figure 2 :
Figure 2: (a) Detailed geologic map showing the location of the Rajapalot Au-Co area.(b) Drill hole location map of Rajapalot and its major prospects and location of investigated drill holes within the Raja area (indicated as white circles).Maps (a) and (b) are based on DigiKP, the digital map database of the Geological Survey of Finland.A-B cross section shown in figure 2c.(c) Block model of the corresponding A-B cross section from figure 2b showing the stratigraphy of the Raja prospect area and location of sampled drill holes included in profile 1 and profile 2. Corresponding NE-SW profiles from profile 1 and location of sampled drill holes: PAL0104, PAL0092, PAL0093, PAL0097 and PAL0100, as well as drill holes from profile 2 (PAL0088 and PAL0116).Block model and profile sections are drawn by Jackson van den Hove (PGN Consulting) are modified from Cook and Hudson (2018) and Cook et al. (2020).

Figure 4 :
Figure 4: Textural setting of three distinct pyrite types.(a) Circled figures on the left show their occurrence in hand specimen.Reflected light photomicrographs show (b) and (c) pyrite A surrounded by fringe-shaped pyrrhotite following the foliation of the host rock; (d) and (e) show pyrite B porphyroclasts with magnetite (mag)-amphibole-chlorite fringes and (f) fracture-filling pyrite C associated with chalcopyrite (cp).

Figure 5 :
Figure 5: Photomicrographs showing the textural occurrence of pyrrhotite, cobaltite, molybdenite (mo), Bi-Te-rich phases and native gold.(a) Sandwich-textured molybdenite-Bi-Te-rich phase assemblage at pyrrhotite rim; (b) Cobaltite I accumulations with Bi-Te-rich phases between grain contacts.This whole assemblage is located along a pyrrhotite rim.(c) Cobaltite I inclusion in pyrrhotite that follows the foliation of the host rock.(d) Cobaltite I at pyrrhotite rim.Assemblage is associated with biotite following rock-foliation.Pentlandite (Pn) forming flame-like inclusions in pyrrhotite.(e) Cobaltite II, native Au and Bi-Te-rich phases in micarich domains.(f) Void-filling native Au textures in domains of sericitization.

Figure 6 :
Figure 6: 3D-model of sampled drill cores showing the distribution of gold and cobalt.Ore zones are characterized as (a) Co-only zones (orange), where Au is b.d.l. and as (b) Au-Co-zones (purple), with Au up to 115 ppm and Co up to 2341 ppm.Barren zones are indicated in grey colors.

Figure 7 :
Figure 7: Downhole plots illustrating the distribution of Au (in log-scale) and Co (in ppm) along the stratigraphic column of sampled drill cores.

Figure 8 :
Figure 8: Logarithmic-scaled boxplot illustrating the range in concentration of metals and pathfinder elements in Co-only zones, Au-Co-zones and barren zones.

Figure 11 :
Figure 11: Logarithmic-scaled boxplot illustrating the range in concentration of metals and pathfinder elements in analyzed pyrite A, B and C by LA-ICP-MS.

A
and B along the positive and negative axis of PC3, respectively.The high positive loadings of Co and the large negative loading of Ni onto PC3, correspond well with the chemistry of Co-rich, but Ni-poor pyrite A. Most of the observation from pyrite B are essentially distributed within the positive PC1 and negative PC3 space, corresponding to the loadings of As, Bi, Au and Se; pyrite C observations are scattered between the loadings of Tl, Ag and W along the negative axis of PC1.The two opposing groups of Au-Bi-Te-Co-As-Se and Ag-Cu-Tl-W-Ni-Mo remain consistent in this projection (Figs 12a-

Figure 13 :
Figure 13: Frequency δ 34 S histograms for (a) analyzed pyrite.(b) Comparison of analyzed δ 34 S-values from pyrite A, B and C.
lower Co-concentrations in Co-only zones (up to 1085 ppm) compared to the Au-Co-zones (up to 2340 ppm),and what could be the reason for the compositional differences within the two distinct Au-Co domains that either give up to 115 ppm Au and 2340 ppm Co when hosted by albitite and mica schists and their contact zones at lower stratigraphic levels, or noticeably lower concentrations in gold (up to 2.67 ppm) and cobalt (up to Au and 160 ppm), when distributed within mafic metavolcanic units and underlying by albite-calcsilicate rocks at higher stratigraphic levels.The following scenarios could provide clarification: (i) During pre-metamorphic albitization the mobilization of basinal fluids that were moderately enriched in cobalt, resulted in the formation of the initial Co-only zones.At a later Svecofennian deformation stage, Au-transporting fluids loaded with Co were deposited during fracturing of the most competent units that now host the highest Au-and Co-grades (fractured albitite and mica schist), compared to the less competent units (mafic metavolcanics) with Au-Co-concentrations. (ii) (ii) Albitization and the resulting generation of varying porosity properties in host rocks, could enable the circulation and precipitation of Co-rich fluids along potential pathways.These units had varying concentrations of Co and were overprinted by late-orogenic Au-only fluids, which preferred the most reactive rock for precipitation.
Co and Ni have a positive correlation in Au-Co-zones, as well as in Co-only zones.At the scale of pyrite crystals however, a somewhat positive correlation of Ni and Co is chiefly observed in pyrite C, while a negative correlation characterizes the pyrite A type (high Co/Ni ratio) and no correlation in the pyrite B type ( Additional significant geochemical patterns observed from the pyrite chemistries are the negative loadings of Mo, Ni, W, Tl, Cu and Ag that are opposing the positive loadings of Co, As, Se, Te, Bi and Au along the axis of PC1.All elements within the negative space of PC1 are enriched in the fracture-related pyrite C, which is part of the late molybdenite-scheelite-native gold-Bi-Te-rich phases and fracture-related chalcopyrite assemblage.In addition, Mo, Ni and W have positive loadings onto PC2, opposing the negative loadings of Tl, Ag and Cu.According toDmitrijeva et al. (2020), the significant W and Mo signals in pyrite represent newly precipitated sub-micron-scale scheelite and molybdenite phases that are the products of a coupled dissolution-reprecipitation reaction of a parent W-Mo-bearing hematite.Regarding the divalent cations in pyrite C, their presence is likely due to the replacement of Fe 2+ ↔ Ni 2+ (or Cu 2+ ).Additionally, the noteworthy concentrations of Cu, Tl, and Ag could be explained by the two following coupled substitutions: (i) 2Fe 2+ ↔Tl + +Sb 3+(D'Orazio et al., 2017) and (ii) 2Fe 2+ from magmatic fluids and from metasediment-related fluids) cannot be excluded and in this latter case the elevated δ 34 SVCDT values may indicate influx of sedimentary sulfur produced by sulphate-reducing bacteria.It appears that pyrite B characterizes a stage between the older Co-rich pyrite A and the younger Ni-rich pyrite C endmember, during which the wide compositional ranges of trace elements and sulfur isotopes either indicate the changing redox conditions of the system, or a mixing of reduced and oxidized fluids from the pyrite A and pyrite C stages, respectively.The documented wider range δ 34 SVCDT values for pyrite C (-1.2 ‰ to +7.4 ‰), as well as the elevated concentrations of redox sensitive elements (Ni, Cu, W Ag Hg, Mo and Zn) may also indicate the occasional shift of conditions to a more oxidizing state.The Ni-rich nature of the fracture-related pyrite C, its close relation to the gold mineralization, and the fluctuations in Cr,55Mn, 57 Fe,59Co,60Ni,61Ni, 62 Ni, 65 Cu, 67 Zn, 69 Ga, 70 Zn, 71 Ga, 72 Ge, 73 Ge, 75 As, 77 Se, 78 Se, 93 Nb, 95 Mo, 97 Mo, 99 Ru, 101 Ru, 103 Rh, 105 Pd, 106 Pd, 107 Ag, 108 Pd, 109 Ag, 111 Cd, 115 In, 118 Sn, 119 Sn, 120 Sn, 121 Sb, 123 Sb, 124 Te, 125 Te, 126 Te, 137 Ba, 181 Ta, 182 W, 183 W, 194 Pt, 195 Pt, 197 Au, 202 Hg, 203 Tl, 205 Tl, 208 Pb, 209 Pb, 232 Th and 238 U.

Table 4 : Results from principal component analysis from the 𝒏 = 𝟐𝟔𝟎 trace element data from pyrite, including the 12 variables in relation to the first three principal components, which contain most of the data variance (see corresponding eigenvalues). The proportion of the variance in each eigenvector is displayed as the percentage of the explained variance. The first three elements with the largest positive and negative loadings along the positive and negative axis of each PC are marked in bold print.
Table3).This wide range of Co/Ni is further reflected in the PC1-PC3 projection, where PC3 is dominated by the Co and Ni loadings, which however have an inverse positioning.Such observations indicate a non-uniform evolution of Co and Ni in pyrite during pre-, syn-, and late-orogenic processes (Table