Next Article in Journal
The Effect of Operating Variables on the Performance of Column Flotation of Silica Sand
Previous Article in Journal
Fluorite and Gibbsite Solubility Controls the Vertical Transport of Fluoride and Aluminum during Rainwater Percolation through Ashfall Deposits in La Palma (Canary Islands, Spain)
Previous Article in Special Issue
Petrology and Geochemistry of Highly Differentiated Tholeiitic Magmas: Granophyres in the Messejana–Plasencia Great Dyke (Central Iberia)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Method Approach to Geophysical Imaging of a Composite Pluton in North Portugal

Departamento de Geociências Ambiente e Ordenamento do Território, Instituto de Ciências da Terra—Pólo Porto, Faculdade de Ciências, Universidade do Porto, Rua do Campo Alegre 687, 4169-007 Porto, Portugal
*
Author to whom correspondence should be addressed.
Minerals 2024, 14(4), 342; https://doi.org/10.3390/min14040342
Submission received: 19 February 2024 / Revised: 18 March 2024 / Accepted: 20 March 2024 / Published: 26 March 2024

Abstract

:
Potassium (K), thorium (Th), and uranium (U) are good markers of magmatic or alteration processes and the surface concentrations of these radioelements can be mapped at the regional or local scale through radiometric (gamma) surveys. In this study, a radiometric survey was performed in a post-orogenic pluton located in North Portugal, namely the Lamas de Olo Pluton, composed by three granitic facies. This pluton has already been intensively studied, including magnetic susceptibility, gravimetric, geochemical, and petrographic studies. The main objective of this work is to evaluate the radiometric data and combine them with other characteristics, such as magnetic susceptibility, and gravimetry, as well as to elucidate structures such as faults and fractures, outline geological boundaries, and identify alteration zones within various granites of the pluton. The radiometric approach reveals the spatial distribution of radioelements, offering a more distinct portrayal of the geology in the studied area. The radioactive heat production rate was calculated for the studied pluton, showing that the mean value is 4.09 µW m−3, surpassing the known mean values for granites. Our study highlights that radiometric measurements unveil compositional variations within granitic pluton and aid in identifying feeder zones. Furthermore, these measurements can be correlated with each type of granites, demonstrating associations with surface concentrations of K-Th-U. Our findings indicate a spatial alignment between the NE feeder root and a U-rich granite (Barragem granite), as evidenced by its elevated concentration of this radioelement. Conversely, the other root displays a notable relative concentration of Th, consistent with the Th-rich characteristics observed of the two other granites (Lamas de Olo and Alto dos Cabeços granites).

Graphical Abstract

1. Introduction

Geophysical techniques form an integral part of most geological mapping works and have a wide application in reconnaissance surveys. In recent years, there has been an increasing use of geophysical methods, owing to their ability to map and delineate subtle physical variations in the geology most often missed during fieldwork (e.g., [1,2,3,4,5,6]).
Geophysical methods used in this work include radiometry, magnetic susceptibility, and gravimetry. These methods, like other geophysical techniques, use specific physical properties such as radioactivity, magnetic susceptibility, and density to map the characteristics of the outcropping and non-outcropping rocks. Both magnetic susceptibility and gravity surveys are useful for delineating geological boundaries. However, gravimetric surveys have a greater depth of investigation and are commonly used to map geological features at a deeper level, allowing us to infer the geometry of a pluton. Concerning radiometric surveys, they are useful in geological mapping, since the different types of rock can be recognized from their distinctive radioactive signature [7,8].
Radiometric surveys measure variations in three radioelements, namely potassium (K%), thorium (Th ppm), and uranium (U ppm) in the soil and exposed rocks. It is important to understand that the radiometric method has several characteristics that make it unique amongst the geophysical methods, namely: (i) the measured radioactivity emanates only from the Earth’s surface’s top few centimeters, and (ii) the identification of elemental sources through the energy of emitted γ-rays allows radiometric data to map variations in the chemical, rather than physical, characteristics of the studied area [2]. The total gamma radiation constitutes a “frozen image” of the last stage of the magmatic building of plutons and is statistically reliable for distinguishing between rock units (e.g., [1]). Therefore, variations in radiation concentration are related to the presence of lithologies with different chemical and mineralogical composition.
During radioactive decay, mass is converted to energy [9]. The radioactive decay of K, Th, and U is recognized as an important heat source in continental crust with typical heat-production values of 0.1–3.0 μW m−3 (e.g., [10,11]).
The generation of heat within a particular rock is influenced by the concentrations of K, Th, and U, which exhibit significant variations depending on the rock type. However, certain patterns emerge due to the consistent geochemical behavior of these elements during processes such as magmatic differentiation, sedimentation, and metamorphism, which contribute to the distribution of these radioelements in rocks. In the case of granites, the heat production varies between 1.25 ± 0.83 μW m−3 and 4.36 ± 2.17 μW m−3 [12], and the average value is reported to be ca. 2.00 μW m−3 in [12] studies, and 2.45 μW m−3 according to [13].
The radioactive heat production rate in a rock can be calculated using the following equation [9]:
H = 10 5 ρ ( 9.52 c U + 2.56 c T h + 3.48   c K ,   in   µ W   m 3
where H represents the radioactive heat production rate, ρ is density (in kg m−3), 9.52 × 10−5, 2.56 × 10−5 and 3.48 × 10−9 are the heat production constants for U, Th, and K, respectively (in W kg−1), and the c U , c T h , and c K are the measured content of U (ppm), Th (ppm), and K (%).
In the present work, the geophysical methods were applied to post-orogenic Variscan granites from the Lamas de Olo pluton. The magnetic susceptibility and gravimetric studies were previously performed by Cruz et al. [5,14,15], and the radiometric approach is now overlapped and analyzed together with the data previously obtained.
The main objective of this work is to map the radiometric response of the Lamas de Olo pluton, emphasizing the correlation between new radiometric data and the magnetic characteristics of the pluton’s facies, its geometry in depth, and the positioning of its feeding zones. The overarching goal is to enhance our understanding of this pluton.

2. Geological Setting

The Iberian Variscan belt, formed by the collision of Laurussia and Gondwana during the Devonian–Carboniferous times, is a significant curved section within the broader European Variscan belt [16,17,18]. In Portugal, three main ductile deformation phases (D1, D2, and D3) have been identified [19]. In the Central Iberian Zone (CIZ), only D1 (ca. 360–337 Ma; e.g., [20,21,22,23,24]), and D3 (ca. 321–300 Ma; e.g., [24,25]) were distinguished. D1 resulted in isoclinal folds with a sub-horizontal axial plane, NW-SE strike, and D3 characterized by wide folding with small amplitude, featuring a subvertical axial plane and NW-SE sub-horizontal axis. Simultaneously, during D3, vertical ductile shear zones developed in lower structural levels, while upper levels experienced brittle deformation, leading to fracture-conjugate systems with NNE-SSW and NNW-SSE orientations [26]. The late-Variscan period is recognized as a time of widespread fracturing in the Variscan belt of southwest Europe [26,27,28].
Northwestern Portugal is occupied by a large volume of granites emplaced during the post-collisional stage of the Variscan orogeny (312–290 Ma; e.g., [24,25,29]).
The target of this study is a pluton located in the northern part of the Central Iberian Zone, the Lamas de Olo pluton, in northern Portugal (LOP, Figure 1).
The LOP is a composite pluton constituted by three outcropping biotite granites, similar in mineralogical composition but with different textures and different biotite contents: (i) Lamas de Olo (LO), (ii) Alto dos Cabeços (AC), and (iii) Barragem (BA) granites. LOP granites are mostly composed of quartz, plagioclase (oligoclase, oligoclase–andesine, and albite), K-feldspars (perthitic orthoclase, and microcline), and biotite. Muscovite I (Ms I), muscovite II (Ms II; post-magmatic muscovite), zircon, sphene, allanite, fluorite, hematite, magnetite, ilmenite, chlorite, rutile, apatite, goethite, epidote, and tourmaline are present as accessory minerals (e.g., [14]). The LO granite is the most abundant in outcrop and is characterized by a medium-to-coarse-grained porphyritic granite (Biotite (Bt) > Muscovite (Ms)). The AC corresponds to a fine-to-medium-grained porphyritic granite (Bt > Ms). The younger granite, which cut the other two granites, is BA which outcrops in the center of the pluton and near a dam and is classified as leucocratic fine-to-medium-grained and slightly porphyritic granite (Bt = Ms I + Ms II) [14,34,35].
The LOP intrudes metasediments from the Douro Group (DG) of the Schist–Graywacke Complex supergroup (CXG), and syntectonic (syn-D3) peraluminous granites from Vila Real Massif (VRM) (Figure 1; [34]). The VRM is a composite massif constituted by two-mica granites, with the monazite and zircon analysis yielding a weighted average 207Pb/235U age of ca. 311 ± 1 Ma [36]. In the field, the contacts between LO and AC granites are generally diffuse; however, field observations show that BA granite crosscuts the LO and the AC granites. 207Pb/235U dating yields an age of 297.19 ± 0.73 Ma to LO granite [35].
Previous studies (e.g., [5]) show that this pluton has a rhombus shape controlled by an NNW-SSE fault system conjugated with both the NNE-SSW and WNW-ESE system faults.
The NNE-SSW system is parallel to the Penacova-Régua-Verin fault, still active and marked by hot CO2-rich mineral water springs [37,38]. The gas-carbonic springs are reported in the localities of Chaves, Vidago, and Pedras Salgadas [37,38]. However, there is no evidence of water springs in the studied area. Regarding the WNW-ESE system fault, commonly affected by late episyenitization processes, this occurs mostly in the ESE border of the pluton. Beyond that, W-Mo(-Sn) mineralizations occurred in association with some of the WNW-ESE quartz veins [5,31].
Recently, and based on the new classification of Hildenbrand et al. [39], LOP is classified as a late-C3 pluton. Expanding on this, Sant’Ovaia et al. [40] delve into the magnetic fabric and show that LOP is an example of a pluton emplaced during the final stage of collision, controlled by regional late C3, NNW-SSE strike-slip faults within a wrench regime.

3. Materials and Methods

3.1. Radiometric Imaging

Two types of data were analyzed in this work: (i) available radiometric data, and (ii) data collected in the field with specific equipment.
(i)
Available radiometric data
The radiometric data were provided by “Unidade de Recursos Minerais e Geofísica” from “Laboratório de Geologia e Minas—Laboratório Nacional de Energia e Geologia” (LGM-LNEG; [41]). The data correspond to an area that covers the “Folha 101—Lordelo-Vila Real“ from the Military Map of Portugal at a scale of 1:25,000. LNEG provided 4118 radiometric measurements of total gamma radiation (nGy h−1).
Different methodologies were used to obtain these radiometric measurements, namely using a scintillator SAPHYMO SRAT, model SPP2. Complementary surveys were made in regional scaling, using spectrometers such as Exploranium GR-256 hand held, and Exploranium GR-650 4 L volume of crystal coupled to a car [41].
(ii)
Radiometric survey in the field
The instrument employed for radiometric surveying is referred to as a gamma-ray spectrometer, and in this study, the portable GR-320 Exploranium model (Manufactured by SAIC Exploranium, 6108 Edwards Blvd, Mississauga, ON, Canada) was used. The portable equipment provides measurements for individual isotopes of K (in %), Th (in ppm), and U (in ppm).
The radiometric survey comprised 57 measuring stations (Figure 1), with three measurements conducted at each station and subsequent calculation of the averages. In the area studied, 48 measurements are performed in LOP granites, and only 9 measurements corresponded to the lithologies of the LOP surrounding rocks (7 in VRM, and 2 in CXG).
The terrain in the study area is characterized by hills, making it challenging to reach numerous outcrops. Accordingly, data collection primarily focused on considering local roads and access points along two main directions, N-S and W-E. Therefore, while the analysis includes the entire pluton, it should be noted that the NW and SW areas of the pluton, as well as the central area of the pluton, composed by BA and AC granite, pose greater challenges in terms of data assessment. Consequently, any conclusions drawn should not be considered definitive but rather justify detailed examination.

3.2. Magnetic Susceptibility and Gravimetric Approach

These techniques were previously explained, and their results were discussed in Cruz et al. [5,14,15].
The magnetic susceptibility measurements were performed using the KLY-4S Kappabridge susceptometer Agico model (Manufactured by Advanced Geoscience Instruments Company, Agico, situated in Brno, Czech Republic) from the Instituto de Ciências da Terra—Polo Porto. At the lab, the data were gathered in a reversible magnetization field, at 0.3 mT and room temperature. The magnetic susceptibility was measured in 350 oriented cores from 48 sampling sites roughly distributed across Lamas de Olo granites [5,14,15].
The gravity measurements were performed using a LaCoste & Romberg G1054 gravimeter (Manufacturered by LaCoste & Romberg, Austin, TX, USA) from the Instituto de Ciências da Terra—Polo do Porto. In total, 271 gravimetric stations were measured along 3 gravimetric profiles both on the LOP granites and host rocks [5]. However, the comparison between magnetic susceptibility, gravimetric response, and radiometric imaging has not been performed yet.

3.3. Coordinate System and Map Specification

The data are in the PT-TM06/ETRS89 coordinate system and the map interpolation was made using Oasis Montaj 2023.2 and ArcGIS PRO 3.2.2 software. The minimum curvature interpolation method was used to create all the interpolated maps in this work, namely: magnetic susceptibility and residual anomaly maps, as well as all the maps related to the radiometric data. The white lines present in Figures 2–7 represent the isolines of the residual anomaly (3rd order) and are superimposed to help the interpretation of the radiometric data.

4. Results and Discussion

4.1. Radiometric Data Interpretation

As a starting approach, the gamma radiation map was interpolated using data provided by LNEG (Figure 2).
The radiometric map (Figure 2) shows the higher values of gamma radiation in the SE zone of the pluton.
Since the provided data only present the total radiation without the individual spectra of K, Th, and U, a more detailed study of the LOP granites has been carried out.
A radiometric survey was conducted, mostly in the LO granite, and the data were analyzed. Several interpolation maps were created using the obtained data: (i) a map displaying the distribution of total radiation values (Figure 3a); (ii) a map showing the distribution of K isotope values (%; Figure 3b); (iii) a map corresponding to the interpolation of the Th isotope values (ppm; Figure 3c); (iv) a map depicting the distribution of U isotope values (ppm; Figure 3d), and (v) a red–green–blue (RGB) ternary map (Figure 4).
As an alternative to the four concentration maps (Figure 3), Figure 4 shows a commonly used approach to visualize gamma-ray spectrometry data, allowing the simultaneous representation of three isotopes. The ternary K–Th–U radiometric map was created through a red–green–blue (RGB) additive chromatic combination applied point by point for the three radiometric measurements. The ternary map was created through the “Grid and Image” tool of the Oasis Montaj 2023.2 software. Despite the impact of weathering and soil development on the area, the ternary K–Th–U radiometric maps can effectively display geological features [1].
In the following discussion concerning the ternary K–Th–U radiometric map, terms such as “K-relative rich” or “K-relative poor” do not depict absolute variations in K content but rather indicate high or low K/(K + Th + U) ratios. This principle similarly applies to the analysis of Th and U content. In the NE section of the ternary map, a blueish zone is evident, reflecting a U-relative high concentration compared to K and Th isotopes. This zone corresponds to the northern outcrop of Barragem granite. Another light blue zone at the south is identified, corresponding to a relatively high U and Th concentration, and depleted in K. At the north, a greenish zone appears, corresponding to a Th-relative high concentration. The reddish zone located at NW in the map suggests the presence of an area richer in K. However, this area is difficult to assess, and no data have been collected on site. Accordingly, interpretations rely on the interpolation of available data, emphasizing the need to refrain from drawing direct conclusions. Much of the eastern part of the ternary map comprises white tones, which are characteristic of a relative high radioelement concentration and represent an area intensely fractured in the Lamas de Olo granite. Yellowish tones appear in the center of LOP, which means that this zone has a relatively high K-Th concentration. It is important to note that this region is predominantly occupied by the BA and AC granites. However, no specific measurements were obtained on these granites. The map is generated through data interpolation based on measurements from the LO granite. Therefore, the observed enrichment in K-Th and relative depletion in U are characteristic features of the LO granite that outcrops in this area. Otherwise, in the western area of LOP, the green dark color is dominant, which corresponds to a relatively low radioelement concentration (Figure 4).
Table 1 displays the descriptive statistics for values obtained across the study area, as well as specifically for the LOP granites and surrounding rocks (mostly VRM). The total radiation data (the sum of K-Th-U isotopes concentration) range between 27.7 and 59.9 ppm (Table 1). Comparing the interpolation maps of the total gamma radiation data (Figure 3a) with the map of total gamma radiation interpolated with the LNEG data, some similarities are noticeable (Figure 2).
Looking to the values of LOP granites, for the K isotope, the minimum value is 2.8% and the maximum is 6.0%; for Th, the minimum is 13.1 ppm and the maximum is 30.3 ppm; for U, the minimum is 5.1 ppm and the maximum is 18.5 ppm. The Th is the isotope that has the higher variation in the LOP granites (the standard deviation is 3.7 ppm), and K is the one that has the lower variability (the standard deviation is 0.6%; Table 1).
A detailed examination of the radiometric data interpolation maps and the geological map reveals that:
  • Both total radiation maps, the one using LNEG data interpolation (Figure 2) and the other using data obtained with portable device (Figure 3a), reveal significant variability across the pluton. The high concentration is located at the east border of LOPl;
  • The higher concentration of Th is spatially related to areas with intense fracturing, coinciding with locations where NNE-SSW and NE-SW fractures intersect with WNW-ESE structures. This correlation is particularly notable in the middle of pluton, and at NE and ESE borders of LOP (Figure 1 and Figure 3c);
  • The U concentration is higher in the eastern area of LOP, exhibiting a NNE-SSW trend. A closer examination shows that the U concentration is highest in the ESE zone of the pluton, where there is intense WNW-ESE fracturing. Additionally, the high U concentration is spatially related to the NNE-SSW fault, parallel to the regional Penacova-Régua-Verin fault (Figure 1 and Figure 3d).
While radiometric data offer insights into isotopic concentration and whole-rock geochemical analysis provides information on the chemical and mineralogical composition of rocks, it is possible to establish a connection between the two.
Previous geochemical studies [15] reveal that LOP granites exhibit characteristics of being calc-alkaline to alkali-calcic as well as peraluminous in nature. Regarding to the rare earth element (REE) spectra, all granites display a flat pattern, with a slight pronounced Eu anomaly in LO and AC granites and a more pronounced Eu anomaly in BA granite. The findings imply the presence of two distinct magma sources: LO and AC granites have geochemical similarities, suggesting a common magma origin, while BA granite originates from a different magma source and represents the most evolved granite. The BA granite, characterized as an alkaline granite, is the youngest and most evolved, displaying REE spectra with a similar pattern to anorogenic granites. These observations suggest that BA granite was the last to ascend and be emplaced in the crust.
The geochemical data (Table 2) show that the major element of K2O (%) has similar values in all the LOP granites ( σ = 0.24) which is in accordance with the variability of K isotope (Figure 3b and Table 1). Regarding the Th average content, is quite similar in the Lamas de Olo ( x ¯ = 15.2 ppm) and Alto dos Cabeços ( x ¯ = 16.5 ppm) granites, having the lower values in the Barragem ( x ¯ = 9.4 ppm) granite. Otherwise, the U element mean values are much higher in the Barragem ( x ¯ = 28.0 ppm) granite, reaching values of 38.9 ppm, while in other granites the values are lower (LO: x ¯ = 10.0 ppm; AC: x ¯ = 7.4 ppm) (Table 2; [10]).
It is noteworthy to highlight that, by comparing with the mean values of the composition of continental crust, namely 5.6 ppm and 1.3 ppm for Th and U values [42], the LOP granites exhibit significantly higher Th and U values, given the mean values of 14.1 ppm and 13.9 ppm, respectively. Works performed by Teixeira et al. [43]), show that Th and U are mostly concentrated in zircon, monazite, allanite, and xenotime.
When comparing geochemical data, particularly focusing on Th and U content, with the interpolation maps of Th and U isotopes, a strong positive correlation is observed. There is an enrichment of all elements in the ESE border of LOP, aligning with a highly fractured area. Additionally, looking at Figure 4, a distinct blue spot is evident near the northernmost outcrop of the Barragem granite, suggesting a higher concentration of U in this region, aligning with the occurrence of the granite with the highest U content (Table 2). As stated earlier, only the northern outcrop of the BA granite was surveyed, which could account for why the blue spot is only depicted in this outcrop and not in the most representative outcrop of the BA, located in the central zone of the pluton.
To better understand the elemental distribution in LOP granites, predictive maps of eU/eTh, eU/K, and eTh/K ratios were made (Figure 5a–c).
The predictive maps for eU/eTh, eU/eK, and eTh/eK confirm the high concentration of U and Th radioelements in LOP, as well as highlight the heterogeneous distribution of radioelements in LOP granites.
Given that radioactive decay generates heat, the radioactive heat production rate was calculated for LOP. Considering the mean density of LOP as 2557 kg m−3 [5] and using mean values of 9.1 ppm, 22.5 ppm and 4.5% for the U, Th, and K content, respectively, the heat production (H) of LOP is determined to be 4.09 µW m−3.
Comparing the H of LOP to the mean value observed in granites worldwide, the LOP exhibits higher mean values. This suggests that the intrusion of this pluton has resulted in moderated heat flow in the crust; however, this is not enough to develop significant geothermal energy resources potential.

4.2. Magnetic Susceptibility

The magnetic susceptibility values show a heterogeneous behavior over the pluton (21 µSI < Km < 44,382 µSI), with both high (≥1000 μSI) and low (<1000 μSI) values of magnetic susceptibility (Km). The Lamas de Olo granite shows a higher variability in the Km values, and the Barragem granite has the one that has the lowest values (e.g., [5,14,15]). The Km higher values (≥1000 μSI) are located at the east border of LOP, as well as in the central zone of the pluton (red and pink areas; Figure 6).
According to the classification of Ishihara [44], LOP belongs to the group of magnetite-type granites (e.g., [5,14,15]). The Km data, combined with petrographic, micro-Raman spectra, and other magnetic mineralogy studies (e.g., thermomagnetic curves and isothermal remanence magnetization), show the presence of pure magnetite, Ti-poor magnetite, and hematite in granites with high and low Km values, albeit in different proportions [14]. Beyond that, the petrographic studies show martitization processes which allow us to infer that LOP granites have originally a magnetite-type behavior, but later processes lead to a partial or total alteration of magnetite into hematite [5,14].
Comparing the spatial distribution of magnetic susceptibility (Figure 6) with the diverse radiometric maps (Figure 2, Figure 3, Figure 4 and Figure 5), no straightforward correlation is evident. These data highlight, once again, the heterogeneous nature of LOP granites, characterized by a wide range of magnetic susceptibility values and a varied distribution of radioelements, as well as geochemical concentrations (Table 2).

4.3. Gravimetric Survey

The gravimetric survey was an efficient way to identify discontinuities, establish the geometry and depth of LOP granites, and infer about the pluton feeder zones. The Bouguer anomaly for each gravimetric station was calculated (using a reduction density of 2.66 g/cm3), and the Bouguer anomaly map was obtained. The Bouguer anomaly in the studied area ranges between −37.53 mGal and −16.22 mGal, having −31.27 mGal as the mean value [5]. The regional and residual anomalies were calculated using 1st, 2nd, and 3rd order polynomial surfaces. Combining the residual maps with geological information (such as geological contacts, LOP outcrop shape, and existing faults), the model that best fits the surface of the Bouguer anomaly was identified as the residual anomaly map obtained through a 3rd order polynomial surface (Figure 7; [5]).
The low negative anomalies are spatially correlated with the presence of less dense rocks, specifically granites (LOP granites and Vila Real two-mica granites). Conversely, the high values, predominantly in the southwest area, correspond to the presence of rocks with high-density values, such as those belonging to the Douro Group Metasedimentary rocks (e.g., schists; Figure 7).
The lower values in the central zone of the residual map are interpreted as the location of the LOP’s feeding zone. Two deeper zones are identified: one in the NE zone of the LOP, and another one on the WSW border of the pluton, slightly displaced to the west of the current contours of the pluton (Figure 7). The location of the feeder zone was deduced based on the lowest residual anomaly values and local geological knowledge, specifically the intersection of late-Variscan fault systems (ENE-WSW with both NNE-SSW and NNW-SSE; [5]).
The correlation, between the radiometric values and the residual anomaly values, was carried out by XY-plots (Figure 8). These plots represent data collected at sampling sites where both radiometric and gravimetric values were acquired, accounting for 44 sites. Figure 8 provides a clear depiction that K isotopes exhibit the least variability across the study area, while Th and U values display significant variations throughout the pluton. However, no clear relationship between radioelements distribution and residual anomaly values was evident, emphasizing the heterogeneous distribution of radioelements across the pluton.
To help our interpretations, the isolines of the residual anomaly (Figure 7) were superimposed on the various radiometric maps. Therefore, it can be concluded that:
  • There is a positive spatial correlation between the northern feeder root (at the NE zone) and the higher concentration in the U isotope (Figure 4 and Figure 5a,b);
  • A positive spatial correlation between the western feeder root and the relative concentration of the Th radioelement, when comparing with K and U (Figure 4 and Figure 5a,c).
Although no direct numerical relation between the K-Th-U and the residual anomaly values (Figure 8) was observed, there is a clear spatial relation between the higher concentration of U and Th in the feeder zones of the BA and LO + AC granites, respectively.
The distribution of U and Th isotopes in the feeder zones corresponds with the results outlined by Cruz et al. [5]. Their study indicates that the Barragem granite, characterized by a higher uranium content, primarily ascends through the NE root, using the intersection of NNE-SSW and WNW-ESE structures. Conversely, the ascent of Lamas de Olo and Alto dos Cabeços granites, richer in Th, predominantly occurred from the west zone root, guided by the intersection of W-E and NW-SE faults.

5. Conclusions

From this study, the following conclusions have been reached:
  • Radiometric measurements obtained through portable gamma-ray spectrometers effectively reveal the compositional heterogeneities within granitic plutons;
  • The LOP’s radioactive heat production rate is 4.09 µW m−3, which is higher than the global average for granites;
  • Feeder zones, identified in the residual anomaly map, can be correlated with the surface location of high-U and high-Th zones, corresponding to the northern and western roots, respectively. This correlation is in line with the geochemical character of the Barragem granite (rich in U), as well as the Lamas de Olo and Alto dos Cabeços granites (rich in Th);
  • Insights into granite petrogenesis and emplacement conditions can be gained through the examination of variations in isotopic distribution illustrated in ternary radiometric maps. These maps provided a “frozen image” capturing the final stages of magmatic construction.

Author Contributions

Conceptualization, C.C., H.S. and F.N.; methodology, C.C.; software, C.C.; validation, C.C., H.S. and F.N.; formal analysis, C.C., H.S. and F.N.; investigation, C.C., H.S. and F.N.; resources, C.C.; data curation, C.C.; writing—original draft preparation, C.C.; writing—review and editing, C.C., H.S. and F.N.; visualization, C.C., H.S. and F.N.; supervision, H.S. and F.N.; project administration, H.S.; funding acquisition, C.C. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is a contracted researcher under the UIDP/04683/2020 project (DOI: 10.54499/UIDP/04683/2020) award by FCT—Foundation for Science and Technology, I.P. This work is supported by national funding awarded by FCT—Foundation for Science and Technology, I.P., projects UIDB/04683/2020 (DOI: 10.54499/UIDB/04683/2020) and UIDP/04683/2020 (DOI: 10.54499/UIDP/04683/2020).

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors thank the Unidade de Recursos Minerais e Geofísica from the Laboratório de Geologia e Minas—Laboratório Nacional de Energia e Geologia for providing radiometric measurements to the authors. The authors also thank Ricardo Ribeiro and António Oliveira for their help in the acquisition of radiometric data in the field. The authors express their gratitude to the anonymous reviewers for dedicating their valuable time, expertise, and providing constructive feedback that contributed significantly to the improvement of this article. The authors also thank Chuthathip Mangkonsu and Jingjing Yang for their guidance and insightful recommendations during the review process.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tartèse, R.; Boulvais, P.; Poujol, M.; Vigneresse, J.L. Granite petrogenesis revealed by combined gravimetric and radiometric imaging. Tectonophysics 2011, 501, 98–103. [Google Scholar] [CrossRef]
  2. Dentith, M.; Mudge, S.T. Geophysics for the Mineral Exploration Geoscientist; Cambridge University Press: Cambridge, UK, 2014; 516p. [Google Scholar]
  3. Fall, M.; Baratoux, D.; Ndiaye, P.M.; Jessell, M.; Baratoux, L. Multi-scale distribution of Potassium. Thorium and Uranium in Paleoproterozoic granites from eastern Senegal. J. Afr. Earth Sci. 2018, 148, 30–51. [Google Scholar] [CrossRef]
  4. Kumwenda, J.; Lackie, M. Geophysical interpretation of the geology of the Stanthorpe region using aeromagnetic, gravity and radiometric data. Explor. Geophys. 2019, 50, 653–666. [Google Scholar] [CrossRef]
  5. Cruz, C.; Sant’ovaia, H.; Raposo, M.I.B.; Lourenço, J.M.; Almeida, F.; Noronha, F. Unraveling the emplacement history of a Po tuguese post-tectonic Variscan pluton using fabrics and gravimetry. J. Struct. Geol. 2021, 153, 104470. [Google Scholar] [CrossRef]
  6. Cruz, C.; Nogueira, P.; Máximo, J.; Noronha, F.; Sant’Ovaia, H. Granite emplacement within a suture shear zone in the Southwestern Iberian Variscan Belt: The Santa Eulália Plutonic Complex. Geol. Soc. 2023, 180, 13. [Google Scholar] [CrossRef]
  7. Moxham, R.M. Natural radioactivity in Washington County, Maryland. Geophysics 1963, 28, 262–272. [Google Scholar] [CrossRef]
  8. Pires, A.C.B.; Harthill, N. Statistical analysis of airborne gamma-ray data for geologic mapping purposes: Crixas-Itapaci area, Goias, Brazil. Geophysics 1989, 54, 1326–1332. [Google Scholar] [CrossRef]
  9. Rybach, L. Determination of heat production rate. In Handbook of Terrestrial Heat Flow Density Determination; Hänel, R., Rybach, L.L., Stegena, Eds.; Kluwer: Dordrecht, The Netherlands, 1988; pp. 125–142. [Google Scholar]
  10. Vilà, M.; Fernández, M.; Jiménez-Munt, I. Radiogenic heat production variability of some common lithological groups and its significance to lithospheric thermal modeling. Tectonophysics 2010, 490, 152–164. [Google Scholar] [CrossRef]
  11. Clark, C.; Fitzsimons, I.C.; Healy, D.; Harley, S.L. How does the continental crust get really hot? Elements 2011, 7, 235–240. [Google Scholar] [CrossRef]
  12. Artemieva, I.M.; Thybo, H.; Jakobsen, K.; Sørensen, N.K.; Nielsen, L.S.K. Heat production in granitic rocks: Global analysis based on a new data compilation GRANITE2017. Earth-Sci. Rev. 2017, 172, 1–26. [Google Scholar] [CrossRef]
  13. Rybach, L. Radioactive Heat Production: A Physical Property Determined by the Chemistry of Rocks. In The Physics and Chemistry of Minerals and Rocks; Strens, R.G.J., Ed.; Wiley & Sons: London, UK, 1976; pp. 309–318. [Google Scholar]
  14. Cruz, C.; Sant’Ovaia, H.; Noronha, F. Magnetic mineralogy of Variscan granites from northern Portugal: An approach to their petrogenesis and metallogenic potential. Geol. Acta 2020, 18, 1–20. [Google Scholar] [CrossRef]
  15. Cruz, C.; Góiis, J.; Sant’Ovaia, H.; Noronha, F. Geostatistical approach in the study of the magnetic susceptibility variation: Lamas de Olo Pluton case study. J. Iber. Geol. 2020, 46, 279–289. [Google Scholar] [CrossRef]
  16. Franke, W. Variscan plate tectonics in Central Europe–current ideas and open questions. Tectonophysics 1989, 169, 221–228. [Google Scholar] [CrossRef]
  17. Ribeiro, A.; Pereira, E.; Dias, R. Structure in the Northwest of the Iberian Peninsula. In Pre-Mesozoic Geology of Iberia; Dallmeyer, R.D., Martínez Garcia, E., Eds.; Springer: Berlin/Heidelberg, Germany, 1990; pp. 220–236. [Google Scholar]
  18. Kroner, U.; Romer, R.L. Two plates-Many subduction zones: The Variscan orogeny reconsidered. Gondwana Res. 2013, 24, 298–329. [Google Scholar] [CrossRef]
  19. Noronha, F.; Ramos, J.M.F.; Rebelo, J.; Ribeiro, A.; Ribeiro, M.L. Essai de corr’elation des phases de d’eformation hercyniennes dans le NW de la p’eninsule Ibérique. Leidse Geol. Meded. 1981, 52, 87–91. [Google Scholar]
  20. Dallmeyer, R.D.; Martinez Catalán, J.R.; Arenas, R.; Gil Ibarguchi, J.I.; Gutiérrez-Alonso, G.; Farias, P.; Bastida, F.; Aller, J. Diachronous Variscan tectonothermal activity in the NW Iberian Massif: Evidence from 40Ar/39Ar dating of regional fabrics. Tectonophysics 1997, 277, 307–337. [Google Scholar] [CrossRef]
  21. Castiñeiras, P.; Villaseca, C.; Barbero, L.; Martín Romera, C. SHRIMP U-Pb zircon dating of anatexis in high-grade migmatite complexes of Central Spain: Implications in the Hercynian evolution of Central Iberia. Int. J. Earth Sci. 2008, 97, 35–50. [Google Scholar] [CrossRef]
  22. Martínez Catalán, J.R.; Rubio Pascual, F.J.; Díez Montes, A.; Díez Fernández, R.; G’omez Barreiro, J.; Dias da Silva, I.; González Clavijo, E.; Ayarza, P.; Alcock, J.E. The Late Variscan HT/LP Metamorphic Event in NW and Central Iberia: Relationships to Crustal Thickening, Extension, Orocline Development and Crustal Evolution. Geol. Soc. 2014, 405, 225–247. [Google Scholar] [CrossRef]
  23. Díez Fernández, R.; Arenas, R.; Pereira, M.F.; Sánchez Martínez, S.; Albert, R.; Martín Parra, L.M.; Rubio Pascual, F.J.; Matas, J. Tectonic evolution of variscan Iberia: Gondwana–Laurussia collision revisited. Earth Sci. Rev. 2016, 162, 269–292. [Google Scholar] [CrossRef]
  24. Pereira, M.F.; Díez Fernández, R.; Gama, C.; Hofmann, M.; Gärtner, A.; Linnemann, U. S-type granite generation and emplacement during a regional switch from extensional to contractional deformation (Central Iberian Zone, Iberian autochthonous domain, Variscan Orogeny). Int. J. Earth Sci. 2018, 107, 251–267. [Google Scholar] [CrossRef]
  25. Ferreira, N.; Iglésias, M.; Noronha, F.; Pereira, E.; Ribeiro, A.; Ribeiro, M.L. Granitóides da Zona Centro Ibérica e seu enquadramento geodinâmico. In Geología de los Granitoides y Rocas Asociadas del Macizo Hesperico. Libro de Homenaje a L.C. García de Figuerola; Bea, F., Carnicero, A., Gonzalo, J., Lopez Plaza, M., Rodriguez Alonso, M., Eds.; Editorial Rueda: Madrid, Spain, 1987; pp. 37–51. [Google Scholar]
  26. Marques, F.O.; Mateus, A.; Tassinari, C. The Late-Variscan fault network in central–northern Portugal (NW Iberia): A re-evaluation. Tectonophysics 2002, 359, 255–270. [Google Scholar] [CrossRef]
  27. Ribeiro, A. Contribution à l’étude de Trás-os-Montes Oriental. Serviços Geológicos de Portugal. Ph.D. Thesis, Memórias Serviços Geológicos de Portugal, Lisboa, Portugal, 1974; 168p. [Google Scholar]
  28. Arthaud, F.; Matte, P. Les décrochements Tardi-Hercyniens du Sud-ouest de l’Europe. Géometrie et essai de reconstitution des conditions de la déformation. Tectonophysics 1975, 25, 139–171. [Google Scholar] [CrossRef]
  29. Dias, G.; Noronha, F.; Almeida, A.; Simões, P.P.; Martins, H.C.B.; Ferreira, N. Geocronologia e petrogénese do plutonismo tardi-Varisco (NW de Portugal): Síntese e inferências sobre o processo de acreção e reciclagem crustal na Zona Centro-Ibérica. In Geologia Clássica, Volume I—Ciências Geológicas: Ensino, Investigação e sua História; Cotelo Neiva, J.M., Ribeiro, A., Victor, M., Noronha, F., Ramalho, M., Eds.; Geologia Clássica, Lisboa, Escolar Editora: Lisboa, Portugal, 2010; pp. 143–160. [Google Scholar]
  30. Pereira, E.; Silva, N.; Moreira, A.; Ribeiro, A. Carta Geológica de Portugal à Escala 1:50,000. Folha 10-A (Celorico de Basto); Serviços Geológicos de Portugal: Lisboa, Portugal, 1987. [Google Scholar]
  31. Helal, B. Granitoïdes, Granites à Métaux Rares et Hydrothermalisme Associe: Géologie, Minéralogie et Géochimie de Plusieurs suites Tardi-Hercyniennes (Nord du Portugal). Ph.D. Thesis, Ecole Nacionale Superieure des Mines de Saint-Etienne, Saint-Étienne, France, 1992. [Google Scholar]
  32. Julivert, M.; Fontboté, J.M.; Ribeiro, A.; Conde, L. Mapa Tectónico de la Península Ibérica y Baleares a Escala 1:1,000,000 y Memoria Explicativa; Instituto Geologico y Mineiro de España: Madrid, Spain, 1974. [Google Scholar]
  33. Farias, P.; Gallastegui, G.; González Lodeiro, L.; Marquínez, J.; Martín Parra, L.M.; Martínez Catalán, J.R.; Pablo Maciá, J.G.; Rodríguez Fernández, L.R. Aportaciones al conocimiento de la litoestratigrafia y estructura de Galicia Central. In Memórias n◦1, IX Reunião Sobre a Geologia do Oeste Peninsular, Museu Laboratário Mineralogia e Geologia-Faculdade de Ciências da Universidade do Porto; Universidad de Oporto: Porto, Portugal, 1987; pp. 411–431. [Google Scholar]
  34. Pereira, E. Notícia Explicativa da Folha 10-A (Celorico de Basto) da Carta Geológica de Portugal: 1:50,000; Serviços Geológicos de Portugal: Lisboa, Portugal, 1989; p. 53. [Google Scholar]
  35. Fernandes, S.; Gomes, M.; Teixeira, R.; Corfu, F. Geochemistry of biotite granites from the Lamas de Olo Pluton, northern Portugal. In Proceedings of the Geophysical Research Abstracts, EGU General Assembly, Vienna, Austria, 7–12 April 2013. [Google Scholar]
  36. Almeida, A.; Leterrier, J.; Noronha, F.; Bertrand, J.M. U-Pb zircon and monazite geochronology of the Hercynian two mica granite composite pluton of Cabeceiras de Basto (Northern Portugal). C. R. L’acad. Sci. 1998, 326, 779–785. [Google Scholar]
  37. Baptista, J.C.; Coke, C.; Dias, R.; Ribeiro, A. Tectónica e Geomorfologia da Região de Pedras Salgadas/Vidago e as Nascentes Minerais Associadas; Comunicações da XII Reunião de Geologia do Oeste Peninsular: Évora, Portugal, 1993; pp. 125–139. [Google Scholar]
  38. Baptista, J.C. Estudo Neotectónico da Zona de Falha Penacova-Régua-Verin. Ph.D. Thesis, UTAD, Vila Real, Portugal, 1998; 344p. [Google Scholar]
  39. Hildenbrand, A.; Marques, F.O.; Quidelleur, X.; Noronha, F. Exhumation history of the Variscan orogen in western Iberia as inferred from new K-Ar and 40Ar/39Ar data on granites from Portugal. Tectonophysics 2021, 812, 228863. [Google Scholar] [CrossRef]
  40. Sant’Ovaia, H.; Cruz, C.; Gonçalves, A.; Nogueira, P.; Noronha, F. Deciphering Iberian Variscan Orogen Magmatism Using the Anisotropy of Magnetic Susceptibility from Granites. Minerals 2024, 14, 309. [Google Scholar] [CrossRef]
  41. Batista, M.J.; Leote Torres, L.; Leote, J.; Prazeres, C.; Saraiva, J.; Carvalho, J. Carta Radiométrica de Portugal (1:500,000); Laboratório de Energia e Geologia: Amadora, Portugal, 2013; ISBN 978-989-675-027-5. [Google Scholar]
  42. Rudnick, R.L.; Gao, S. The Composition of the Continental Crust. In Treatise on Geochemistry; Holland, H.D., Turekian, K.K., Eds.; The Crust, Elsevier-Pergamon: Oxford, UK, 2003; Volume 3, pp. 1–64. [Google Scholar] [CrossRef]
  43. Teixeira, R.J.S.; Gomes, M.E.P.; Martins, L.M.O.; Pereira, A.J.C.S.; Neves, L.J.P.F. Natural Radiation and Geochemistry of the Lamas de Olo Biotite Granite, Northern Portugal. In Proceedings of the Goldschmidt, Sacramento, CA, USA, 8–13 June 2014; p. 2458. [Google Scholar]
  44. Ishihara, S. The Magnetite-series and Ilmenite-series Granitic Rocks. Min. Geol. 1977, 27, 292–305. [Google Scholar]
Figure 1. (a) Simplified geologic map of northern and central Portugal with the plutons and principal faults and shear zones (adapted from Ferreira et al. [25]); the red square indicates the pluton under study: the Lamas de Olo pluton. (b) Simplified geological map of the Lamas de Olo Pluton (modified: [30,31]) with radiometric measurement points; in the Iberian Variscan belt the location of the studied pluton was identified by a red star. LOP—Lamas de Olo Pluton. CZ—Cantabric Zone, WALZ—West Asturian Leonese Zone, GTMZ—Galicia Trás-os-Montes Zone, CIZ—Central Iberian Zone, OMZ—Ossa Morena Zone, SPZ—Sul Portuguesa Zone (based on: [32,33]).
Figure 1. (a) Simplified geologic map of northern and central Portugal with the plutons and principal faults and shear zones (adapted from Ferreira et al. [25]); the red square indicates the pluton under study: the Lamas de Olo pluton. (b) Simplified geological map of the Lamas de Olo Pluton (modified: [30,31]) with radiometric measurement points; in the Iberian Variscan belt the location of the studied pluton was identified by a red star. LOP—Lamas de Olo Pluton. CZ—Cantabric Zone, WALZ—West Asturian Leonese Zone, GTMZ—Galicia Trás-os-Montes Zone, CIZ—Central Iberian Zone, OMZ—Ossa Morena Zone, SPZ—Sul Portuguesa Zone (based on: [32,33]).
Minerals 14 00342 g001
Figure 2. Interpolation of radiometric data in nGy h−1 (values provided by LNEG) in the LOP (Radiometric data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software: modified from: [41]).
Figure 2. Interpolation of radiometric data in nGy h−1 (values provided by LNEG) in the LOP (Radiometric data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software: modified from: [41]).
Minerals 14 00342 g002
Figure 3. Radiometric data interpolated maps of the distribution of (a) total radiation values (ppm); (b) K values (%); (c) Th values (ppm); and (d) U values (ppm; radiometric data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software).
Figure 3. Radiometric data interpolated maps of the distribution of (a) total radiation values (ppm); (b) K values (%); (c) Th values (ppm); and (d) U values (ppm; radiometric data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software).
Minerals 14 00342 g003
Figure 4. Composite ternary RGB radiometric map, with K, Th, and U in the red (R), green (G), and blue (B) channels, respectively (map made in the Oasis Montaj 2023.2 software).
Figure 4. Composite ternary RGB radiometric map, with K, Th, and U in the red (R), green (G), and blue (B) channels, respectively (map made in the Oasis Montaj 2023.2 software).
Minerals 14 00342 g004
Figure 5. Predictions maps of (a) eU/eTh, (b) eU/eK, and (c) eTh/eK ratios (maps made in Oasis Montaj 2023.2 software).
Figure 5. Predictions maps of (a) eU/eTh, (b) eU/eK, and (c) eTh/eK ratios (maps made in Oasis Montaj 2023.2 software).
Minerals 14 00342 g005
Figure 6. Magnetic susceptibly map of the Lamas de Olo Pluton (Magnetic susceptibility data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software; modified from Cruz et al. [5,14,15]).
Figure 6. Magnetic susceptibly map of the Lamas de Olo Pluton (Magnetic susceptibility data interpolation made by minimum curvature in the Oasis Montaj 2023.2 software; modified from Cruz et al. [5,14,15]).
Minerals 14 00342 g006
Figure 7. Residual anomaly map (mGal) of Lamas de Olo pluton and host rocks (residual map interpolated through the Oasis Montaj 2023.2 software and modified from Cruz et al. [5]). The white lines represent the isolines of the residual anomaly (3rd order) and are superimposed on the maps in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 to help the interpretation of the radiometric data.
Figure 7. Residual anomaly map (mGal) of Lamas de Olo pluton and host rocks (residual map interpolated through the Oasis Montaj 2023.2 software and modified from Cruz et al. [5]). The white lines represent the isolines of the residual anomaly (3rd order) and are superimposed on the maps in Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 to help the interpretation of the radiometric data.
Minerals 14 00342 g007
Figure 8. XY-plots of radiometric values versus (vs.) residual anomaly values (3rd order in mGal) of Lamas de Olo pluton and host rocks: (a) TOT mean (total radiation in ppm) vs. residual anomaly; (b) eK (values of K isotope in %) vs. residual anomaly; (c) eTh (values of Th isotope in ppm) vs. residual anomaly; (d) eU (values of U isotope in ppm) vs. residual anomaly.
Figure 8. XY-plots of radiometric values versus (vs.) residual anomaly values (3rd order in mGal) of Lamas de Olo pluton and host rocks: (a) TOT mean (total radiation in ppm) vs. residual anomaly; (b) eK (values of K isotope in %) vs. residual anomaly; (c) eTh (values of Th isotope in ppm) vs. residual anomaly; (d) eU (values of U isotope in ppm) vs. residual anomaly.
Minerals 14 00342 g008
Table 1. Descriptive statistics of radiometric data in the study area ( x ¯ : mean value; σ : standard deviation; min.: minimum value; max.: maximum value).
Table 1. Descriptive statistics of radiometric data in the study area ( x ¯ : mean value; σ : standard deviation; min.: minimum value; max.: maximum value).
Total (ppm)K (%)Th (ppm)U (ppm)eU/eTheTh/KeU/K
LOP + surrounding rocks x ¯ 41.14.422.18.80.45.02.0
σ 6.20.63.92.70.10.50.5
Min.27.72.812.93.30.23.30.9
Max.59.96.030.318.51.26.63.9
LOP x ¯ 42.04.522.59.10.45.02.0
σ 6.00.63.72.70.10.60.5
Min.27.72.813.15.10.23.31.3
Max.59.96.030.318.51.26.63.9
Surrounding rocks x ¯ 36.14.019.76.80.44.91.7
σ 5.00.54.11.60.10.80.3
Min.27.83.212.93.30.23.90.9
Max.43.84.724.68.90.56.01.9
Table 2. Descriptive statistics of major and trace elements in LOP granites based on geochemical data ( x ¯ : mean value; σ : standard deviation; Min.: minimum value; Max.: maximum value; adapted from: Cruz et al. [15]).
Table 2. Descriptive statistics of major and trace elements in LOP granites based on geochemical data ( x ¯ : mean value; σ : standard deviation; Min.: minimum value; Max.: maximum value; adapted from: Cruz et al. [15]).
Lamas de Olo GraniteAlto dos Cabeços GraniteBarragem GraniteLOP
x ¯ σ Min.Max. x ¯ σ Min.Max. x ¯ σ Min.Max. x ¯ σ Min.Max.
SiO2 (%)73.631.4171.6975.0472.501.2771.1173.6076.840.7376.2377.6574.152.0471.1177.65
TiO2 (%)0.1850.020.1680.2060.2410.030.2150.2790.0910.0070.0840.0980.1750.0590.0840.279
Al2O3 (%)13.720.2913.3014.0413.500.3313.2113.8612.340.0312.3112.3713.320.6412.3114.04
Fe2O3 (%)0.950.180.601.130.770.180.640.970.300.100.210.410.740.320.211.13
FeO (%)1.00.20.71.21.20.11.11.30.70.20.50.90.90.20.51.3
Fe2O3 (T%)2.010.121.912.232.070.311.872.421.080.300.831.411.790.480.832.42
MnO (%)0.0590.010.0390.0670.0600.020.0430.0750.0310.0180.0200.0520.0520.0180.0200.075
MgO (%)0.330.050.290.410.500.080.460.590.160.010.160.170.330.130.160.59
CaO (%)0.980.140.731.161.260.260.961.420.390.140.290.550.900.370.291.42
Na2O (%)3.490.133.393.743.000.382.563.233.100.092.993.153.270.302.563.74
K2O (%)4.690.214.394.984.260.094.184.364.520.144.414.674.540.244.184.98
P2O5 (%)0.060.010.050.070.070.010.070.080.040.010.030.040.060.020.030.08
LOI (%)0.540.150.390.731.280.810.662.190.810.260.510.990.790.490.392.19
Total (%)99.580.8598.34100.598.630.3598.3199.0099.320.5798.9499.9799.280.7698.31100.5
Sc (ppm)5055505551455045
Be (ppm)52410735114873613217364132
V (ppm)1511417223202582710155725
Cr (ppm)117288015011017100130127117402601185440260
Co (ppm)3033303311123113
Ni (ppm)138103033630403738108024211080
Cu (ppm)50552010103060692014023385140
Zn (ppm)341015404764050554915110432415110
Ga (ppm)1811718181171817117181811718
Ge (ppm)2023202231232023
As (ppm)3033303352373237
Rb (ppm)24215223265253292212773041529132126032221321
Sr (ppm)98138111414751421522792137934621152
Y (ppm)33329372571933371328523281952
Zr (ppm)1077100118123141131395074558973045139
Nb (ppm)181162018216203711274823101648
Mo (ppm)3124202242263126
Ag (ppm)0.30.10.30.50.30.00.30.30.40.20.30.60.30.10.30.6
In (ppm)0.10.00.10.10.10.00.10.10.10.00.10.10.10.00.10.1
Sn (ppm)72310946139471483314
Sb (ppm)0.30.00.30.30.30.00.30.30.30.00.30.30.30.00.30.3
Cs (ppm)11.13.87.618.513.94.110.318.411.82.39.413.912.03.57.618.5
Ba (ppm)388782824835315348158781455013234718250587
Hf (ppm)3.30.33.03.83.70.23.53.92.60.12.52.73.20.52.53.9
Ta (ppm)4.30.63.75.23.71.22.85.010.92.97.813.65.83.42.813.6
W (ppm)4228303373495329
Tl (ppm)1.30.11.21.41.40.11.31.51.50.21.31.61.40.11.21.6
Pb (ppm)3112932291283041635473362847
Bi (ppm)0.40.30.20.80.20.00.20.26.08.60.716.01.74.50.216.0
Th (ppm)15.21.813.218.316.51.215.617.99.40.19.49.514.13.19.418.3
U (ppm)10.05.65.620.97.41.85.49.028.014.112.138.913.911.25.438.9
La (ppm)20.62.918.526.029.53.527.233.512.10.212.012.320.76.912.033.5
Ce (ppm)42.35.136.450.655.45.150.260.322.95.817.428.940.713.017.460.3
Pr (ppm)4.820.574.285.866.510.526.197.113.120.372.823.534.821.342.827.11
Nd (ppm)18.52.316.622.823.41.122.624.711.51.210.612.918.04.810.624.7
Sm (ppm)4.50.63.95.64.80.34.65.13.60.53.03.94.30.73.05.6
Eu (ppm)0.580.040.540.640.710.040.680.750.220.040.180.250.520.190.180.75
Gd (ppm)4.70.44.35.54.40.83.75.34.21.03.35.24.50.73.35.5
Tb (ppm)0.90.10.71.00.70.20.60.90.90.40.61.30.80.20.61.3
Dy (ppm)5.20.54.55.84.41.33.35.86.22.44.38.95.31.43.38.9
Ho (ppm)1.10.10.91.20.90.30.71.21.30.50.91.81.10.30.71.8
Er (ppm)3.30.42.73.82.70.72.13.54.11.62.85.83.30.92.15.8
Tm (ppm)0.530.060.450.610.420.120.330.560.710.250.510.990.550.160.330.99
Yb (ppm)3.60.33.14.13.00.92.44.05.52.03.87.73.91.42.47.7
Lu (ppm)0.560.060.470.640.450.090.380.550.890.320.621.240.620.230.381.24
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cruz, C.; Noronha, F.; Sant’Ovaia, H. A Multi-Method Approach to Geophysical Imaging of a Composite Pluton in North Portugal. Minerals 2024, 14, 342. https://doi.org/10.3390/min14040342

AMA Style

Cruz C, Noronha F, Sant’Ovaia H. A Multi-Method Approach to Geophysical Imaging of a Composite Pluton in North Portugal. Minerals. 2024; 14(4):342. https://doi.org/10.3390/min14040342

Chicago/Turabian Style

Cruz, Cláudia, Fernando Noronha, and Helena Sant’Ovaia. 2024. "A Multi-Method Approach to Geophysical Imaging of a Composite Pluton in North Portugal" Minerals 14, no. 4: 342. https://doi.org/10.3390/min14040342

APA Style

Cruz, C., Noronha, F., & Sant’Ovaia, H. (2024). A Multi-Method Approach to Geophysical Imaging of a Composite Pluton in North Portugal. Minerals, 14(4), 342. https://doi.org/10.3390/min14040342

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop