Next Article in Journal
Evidence of Seismic-Related Liquefaction Processes within the Volcanic Record of the Campi Flegrei Caldera (Italy)
Previous Article in Journal
Improved Planktonic, Benthic Foraminiferal and Nannofossil Biostratigraphy Aids the Interpretation of the Evolution at Hole U1468A: IODP Expedition 359, the Maldives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Petrophysical Characterization of Non-Magnetic Granites; Density and Magnetic Susceptibility Relationships

by
Emilio L. Pueyo
1,*,
Teresa Román-Berdiel
2,
Pablo Calvín
1,
Jean Luc Bouchez
3,
Elisabet Beamud
4,
Conxi Ayala
5,6,
Francesca Loi
2,
Ruth Soto
1,
Pilar Clariana
1,
Aina Margalef
7,
Nuria Bach
4,
Nia Schamuells
6,
Félix M. Rubio
5,
Ana Gimeno
5,
Elena Fernández de Arévalo
5,
Carmen Rey-Moral
5,
Arturo García
5,
Joan Martí
4,
Antonio M. Casas
2 and
José Luis García-Lobón
5
1
IGME (CSIC), Unidad de Zaragoza, C/Manuel Lasala 44, 9 B, 50006 Zaragoza, Spain
2
Geotransfer Research Group, Departamento de Ciencias de la Tierra, Instituto de Investigación en Ciencias Ambientales (IUCA), Universidad de Zaragoza, 50006 Zaragoza, Spain
3
Géosciences Environnement Laboratory, Federal University of Toulouse, Avenue E. Belin, 31400 Toulouse, France
4
Paleomagnetic Laboratory CCiTUB-Geo3Bcn, Geosciences Barcelona–CSIC, C/Lluis Solé i Sabarís s/n, 08028 Barcelona, Spain
5
IGME (CSIC), C/La Calera 1, Tres Cantos, 28760 Madrid, Spain
6
Currently Visiting Researcher at Geosciences Barcelona–CSIC, C/Lluis Solé i Sabarís s/n, 08028 Barcelona, Spain
7
Andorra Research + Innovation (AR + I), Av Rocafort 21-23, 3, AD600 Sant Julià de Lòria, Andorra
*
Author to whom correspondence should be addressed.
Geosciences 2022, 12(6), 240; https://doi.org/10.3390/geosciences12060240
Submission received: 4 May 2022 / Revised: 4 May 2022 / Accepted: 24 May 2022 / Published: 8 June 2022

Abstract

:
In this work we establish reliable correlations between density and magnetic susceptibility in three paramagnetic granites from the Pyrenees. In total, 128 sites (310 density measurements and >2600 susceptibility ones) were studied in the Mont Louis-Andorra, Maladeta and Marimanha granitic plutons covering the main range of variability of magnetic susceptibility. Regressions were calculated for every granitic body and an integrated linear function was obtained for the entire dataset: ρ (kg/m3) = 2566 (kg/m3) + 0.541κ (10−6 S.I.) (R:0.97). This relationship is only valid in the paramagnetic domain, where iron is mostly fractioned in iron-bearing phyllosilicates and the occurrence of magnetite is negligible (or at least its contribution to the bulk susceptibility). This relationship, likely different in other bodies, allows for transforming magnetic susceptibility data into density data, helping to constrain gravity modelling when density data from rock samples are scarce. Given the large amount of AMS studies worldwide, together with the quickness and cost-effectiveness of susceptibility measurements with portable devices, this methodology allows for densifying and homogenizing the petrophysical data when modelling granite rock volumes based on both magnetic and gravimetric signals.

1. Introduction

Beyond classic mineral exploration in hydrothermally altered or contact metamorphism rocks (e.g., Joly et al. [1]; Dentith and Mudge [2]) or as radioactive waste disposals (e.g., Wang [3]), geophysical surveying of granite bodies is increasingly of interest due to, among other aspects, their potential in geothermal energy (Genter et al. [4] Huenges and Ledru [5]; Moore et al. [6]; Zhang and Zhao [7]). Potential fields geophysical techniques, both gravity and magnetism, have been of great help for 2D and 3D modeling of granitic bodies for a long time (Bott et al. [8]; Bott and Smithson [9], Henkel [10] Vigneresse [11] Ameglio et al. [12]; Cruz et al. [13]) inasmuch they are quick, repeatedly resolute and cost-effective methods for an initial characterization of those bodies at depth. As in any potential-field study, petrophysical data are an essential constraint to link the geophysical signal with the geological features, and thus a keystone to reduce the ambiguity and uncertainty in the interpretations (Henkel [10,14,15]; Enkin et al. [16]; Dentith et al. [17]; Pueyo et al. [18]).
Density and susceptibility relationships have been well-known in granitic rocks from a petrological and petrophysical point of view since the pioneering works by Henkel ([10,14]). Some works have demonstrated that δ18O and the wt.% SiO2 display a negative correlation to magnetic susceptibility in granites (Ellwood and Werner [19], Villaseca et al. [20] respectively). Moreover, it is well known that the magnetic properties of granites depend on their chemical and tectonic affinity (Kanaya and Ishihara [21]; Ishihara [22], 1981 [23]) and they are classified as magnetic granites (κ ranging from 10−3 to 10−2 S. I.) and non-magnetic granites (κ from 10−5 to 10−4 S.I.) (Ellwood & Werner [19]). Non-magnetic granites usually correspond to supracrustal sources, S-type (Chappell & White [24]), where iron is fractionated mainly in ilmenite (paramagnetic at room temperature) and biotite (strongly paramagnetic; Martín-Hernández & Hirt [25]). In magnetic granites, usually related to deeper and/or igneous sources (I-type), the iron mainly forms magnetite crystals (ferrimagnetic). Elming [26] in the Caledonides of Jämtland (Sweden), or more recently Terrinha et al. [27] in Sintra granite, showed coarse relations between density and magnetic susceptibility along five orders of magnitude (κ from 10−5 to 10−1 S.I.). Criss and Champion [28] measured the density in paleomagnetic samples of the Idaho granite and recognized that variations in susceptibility were mostly related to the rock density at a regional scale. Bourne [29] stablished that density and magnetic susceptibility (separately) inversely correlate with SiO2 concentrations in ferromagnetic granites (ca. κ > 500 × 10−6 S.I.; Bouchez [30,31]).
A big step further was given by Gleizes [32]; Bouchez et al. [33]; Gleizes et al. [34] and in the Foix and Mont-Louis Andorra granites, respectively (in the Axial Zone of the Pyrenees). Gleizes et al. [34] related the magnetic susceptibility with the iron composition [Fe+2 and Fe+3] and the expected susceptibility derived from theoretical estimates (after Rochette [35] and Rochette et al. [36]), and also proposed a correlation between magnetic susceptibility and magmatic facies in paramagnetic granites (ca. κ < 500 × 10−6 S.I). Their work has had a profound impact in the application of anisotropy of magnetic susceptibility (AMS) techniques to characterize the internal structure and composition of granitic bodies. Currently, hundreds of granitic bodies all over the world have been systematically studied by AMS, helping us to understand their emplacement models and generating vast databases of susceptibility data (Román-Berdiel et al. [37]; Aranguren et al. [38]; Trindade et al. [39]; Sant’Ovaia et al. [40]; Ferré et al. [41]; Kratinová et al. [42]; Joly et al. [43]; Porquet et al. [44] among many others).
Focusing only on paramagnetic granites (in terms of Bouchez [30], or on the paramagnetic trend (in terms of Henkel [15] and Enkin et al. [16]), very little has been done comparing these petrophysical properties. Ameglio et al. [12] and references therein) found fine correlations among density and magnetic susceptibility in a few samples from three calc-alkaline S-type (paramagnetic) granites from French and Spanish Variscan massifs (namely Sidobre, Cabeza de Araya and Mont Louis-Andorra).
In this paper, the main goal is to build reliable correlations between magnetic susceptibility and density in paramagnetic granites. Three plutons from the Pyrenees are the study-cases: Mont Louis-Andorra, Maladeta and Marimanha. Previous AMS studies in these bodies (Gleizes et al. [34]; Leblanc et al. [45]; Antolín-Tomás et al. [46], respectively) allow the conversion of a vast susceptibility database into robust and reliable density data ready to be modelled together with gravity signals.

2. Geological Setting

The granitic massifs of Maladeta (MAL), Marimanha (MAR) and Mount-Louis Andorra (MLA) crop out in the Axial Zone of the Pyrenees and belong to the tardi-tectonic granitic bodies of the Variscan Orogen (Figure 1). All three plutons were intruded in the upper crust during the latest stages of the Variscan orogeny (Gleizes and Bouchez [47]; Gleizes [32]; Leblanc et al. [45]; Antolín-Tomás [48]; Antolín-Tomás et al. [46]).
The Variscan belt of Europe is part of a vast Paleozoic chain built between 500 and 250 Ma, due to the convergence and collision of two large continental masses, Laurentia-Baltica and Africa (Matte. [49]). It is a sinuous orogen that can be followed discontinuously from the south of Spain (Martínez-Catalán [50]; Pastor-Galán et al. [51]) to the Bohemian massif, and which very probably extends under the Carpathians to the Variscan Caucasus (Figure 1A). It is therefore a Paleozoic chain of nearly 5000 km in length that constitutes the southwest edge of stable Europe (Matte [52]); 250 Ma ago, most of the Variscan belt was completely eroded and subsequent extension and formation of basins occurred during the Mesozoic (Pangea breakup). Then, in the Cenozoic, the Alpine orogeny developed a new collisional chain between Iberia and the SW margin of the Eurasian plate between 65 and 20 M.a. (e.g., Muñoz [52]) resulting in the current Pyrenean Mountain range. The Alpine orogeny gave rise to the uplift of Paleozoic units in the core of the Pyrenean range (i.e., the Axial zone) during the collision by developing a central antiformal stack of several southward-facing basement-involved thrust sheets (Muñoz [53]; Martínez-Peña and Casas-Sainz [54]; Casas et al. [55]). Recently, some authors claim that the persistence of a relatively flat envelope for the Paleozoic sedimentary pile and Variscan isograds, and the absence of Alpine crustal-scale faults in the core of the Axial Zone, suggest that the Axial Zone constitutes a large Variscan structural unit preserved during Pyrenean orogeny (Cochelin et al. [56]).
In this context, the Marimanha granite crops out within the Alpine Gavarnie/Nogueras thrust sheet, the uppermost unit of the south-verging piggyback-sequence in this part of the Pyrenees, whereas the Maladeta and Andorra-Mont-Louis granites, two of the largest granitic bodies of the Axial Pyrenees, crop out in the underlying Bielsa/Orri thrust sheet (Figure 1B).
During the Variscan orogeny, low-grade Paleozoic metasedimentary rocks were affected by progressive polyphasic deformation (Druguet [57], García-Sansegundo et al. [58], Casas et al. [59]), which mostly occurred during a compressional tectonic setting in which fold and thrust systems developed, giving rise to the crustal thickening of the Variscan cordillera in this region (Soula et al. [60]; Carreras and Capella [61], Gutierrez-Medina [62]; Clariana and García-Sansegundo [63]). Subsequently, in the last states of the variscan deformation, an extensional deformation event occurred. According to Soula [64]) and Autran [65]) the main deformation event is synchronous with the peak metamorphism, while other authors (Guitard [66], Zwart [67,68]; Liesa [69] and Aguilar [70]) consider the deformation event to have continued during and after the metamorphic peak. Early south-verging thrust sheets involve Silurian to Carboniferous rocks in the hanging wall and Cambro-Ordovician rocks in the footwall (e.g., Majesté-Menjoulas [71]; Raymond [72]; Bodin and Ledru [73]; Losantos et al. [74]). The Silurian slates act commonly as a detachment level between the two units, and the deformation above it is characterized by fault propagation folds affecting Silurian and Devonian rocks (García-Sansegundo [75], Clariana [76], Margalef [77]). The last stages of the variscan deformation have been characterized as dextral shear motion accompanied by granite intrusion (Leblanc et al. [78]; Evans et al. [79]; Gleizes et al. [80,81,82]; Olivier et al. [83]; Román-Berdiel et al. [84,85]; Auréjac et al. [86]). Due to crustal thickening, the Variscan segment of the Pyrenees experienced crustal flow, gneiss dome formation and subsequent granitic massifs intrusion in the upper crust by in- situ ballooning favored by the boundary between the Cambro-Ordovician and Siluro-Devonian rocks (Antolín-Tomás et al. [46]). The U–Pb ages published for the Pyrenean granites indicate that the Variscan plutonism of the Pyrenees is mainly Carboniferous (Romer and Soler [87]; Paquette et al. [88]; Guerrot [89,90]; Roberts et al. [91]; Maurel et al. [92]; Olivier et al. [93]; Gleizes et al. [82]) and Permian in age (Denèle et al. [94]). Recent works claim that partial melting, crustal flow, gneiss dome formation and pluton emplacement occurred over a very short period at the time scale of the Variscan belt formation, on the order of 5 Ma, at ca. 304 Ma (Denèle et al. [95]). In general, these magmatic bodies follow the dominant trend of the Variscan structure.
Figure 1. (A) European Variscan belt (Modified after Franke [96]); CZ-Cantabrian Zone, WALZ-West Asturian-Leonese Zone, GTOMZ-Galicia Tras-os-Montes Zone, CIZ-Central Iberian Zone, OMZ-Ossa Morena Zone, SPZ-Southportuguese Zone, PAZ-Pyrenean Axial Zone, CCR-Catalonian Coastal Ranges, NPM-Nord-Pyrenean Massifs, MB-Massifs Basques. (B) Geological sketch map of the Central part of the Pyrenees, showing the situation of the Pyrenean granites.
Figure 1. (A) European Variscan belt (Modified after Franke [96]); CZ-Cantabrian Zone, WALZ-West Asturian-Leonese Zone, GTOMZ-Galicia Tras-os-Montes Zone, CIZ-Central Iberian Zone, OMZ-Ossa Morena Zone, SPZ-Southportuguese Zone, PAZ-Pyrenean Axial Zone, CCR-Catalonian Coastal Ranges, NPM-Nord-Pyrenean Massifs, MB-Massifs Basques. (B) Geological sketch map of the Central part of the Pyrenees, showing the situation of the Pyrenean granites.
Geosciences 12 00240 g001
The Marimanha granite is a triangular pluton in map view, with an outcrop area of about 45 km2. The available geochronological data for the Marimanha pluton consists of a poorly constrained Rb–Sr age (on whole rock) of 315–290 Ma (Palau i Ramirez [96]). At the petrological level, this massif includes gabbro, diorite, granodiorite and leucogranite facies (Palau i Ramirez [97]). Its intrusion occurred in the western end of the Pallaresa Dome, in Cambro–Ordovician siliciclastic metasediments, cutting the Silurian and Devonian limestones and slates involved in the Variscan Roca Blanca thrust (Losantos et al. [74]; Bodin and Ledru [73]). The metamorphic aureole overprints the Variscan structures, and the whole is cut by Late Variscan and/or Alpine WNW–ESE reverse faults.
The Maladeta granite massif (Evans [98]; Leblanc et al. [45]; Evans et al. [79]) is elongated in the E-W direction with an area of 400 km2. Available geochronological data provide a rather young date of 277 ± 7 Ma (Rb/Sr; whole-rock, Vitrac-Michard et al. [99]). Petrographic zoning is defined in the western Aneto unit, with more basic rock-types surrounding others that appear to be more silicic, and the presence of gabbros is defined in the south-eastern part of the massif (Charlet [100,101]). The massif is penetratively deformed in some places near its southern edge, and major fault zones separate different units. It was emplaced into Cambro-Ordovician to Carboniferous rocks at a shallow level of about 6 km in depth (Pouget et al. [102]). It developed a contact metamorphism aureole in which the sillimanite grade has been reached.
The Mont-Louis-Andorra massif (Autran et al. [103]; Debon et al. [104]; Gleizes & Bouchez [47]; Gleizes et al. [34]; Bouchez and Gleizes [105]) is elongated in the E-W direction over 55 km and covers nearly 600 km2. Recent available geochronological data provide an age of 301–303 Ma (U/Pb, Denèle et al. [95]). Petrographically, this massif consists of leucogranites, monzogranites and granodiorites with dioritic microgranular enclaves. Its intrusion occurred mainly in Ediacaran-Cambrian metasediments (Laumonier et al. [106]) and its top reaches Devonian sediments to the SW (Andorra). It developed a metamorphic contact aureole with isograds cutting across the regional metamorphic isograds (Guitard et al. [107]; Laumonier et al. [108]).
The bulk susceptibility from almost 600 studied AMS sites (Figure 2) displays very similar distributions for the three plutons and non-significant differences in the mean values (Table 1). These available and densely sampled nets of AMS data from the 90s and early 2000s (Figure 2) are revisited with new data acquired during the development of this work in the three granites to establish reliable correlations between density and magnetic susceptibility.

3. Methodology

In total, both properties (density and magnetic susceptibility) were studied in 128 sites from the Mont Louis-Andorra (MLA), Maladeta (MAL) and Marimanha (MAR) plutons (71, 42 and 15 sites, respectively) in the Central Pyrenees (Figure 3).
Four different types of samples were considered (Figure 4): (A) large hand samples of a few kilograms (type 1) were taken in 92 new sites (MLA and MAL plutons), (B) from them, subsampling of smaller blocks (mini-blocks; type 2; 20 to 200 g) were measured in 48 of those sites. (C) several previous samples from AMS studies (paleomagnetic standard specimens, type 3, about 25–30 g) were remeasured in 36 sites from MLA (21 sites; Gleizes et al. [33]) and MAR (15 sites; Antolín et al. [45]). Finally, (D) subsampling of some standard cores with a 0.5 cm Ø non-magnetic drill bit (1 cm in length, type 4-minicores, below 1 g) in the MLA granite allowed for the determination of susceptibility and other magnetic properties using vibrating and superconducting magnetometers. All the studied samples from these 128 sites tried to evenly cover the entire susceptibility spectra of the studied plutons (Figure 2).

3.1. Density Determinations

In type 1 samples (large hand blocks), the apparent density estimation was performed following the European standard UNE EN 1936:2006 (Natural stone test methods-Determination of real density and apparent density, and of total and open porosity [CEN/TC 246-Natural stones; CEN/TC 246/WG 2-Test methods]). More than two density data were taken in average per site (ranging between 1 and 5) at the IGME laboratories (Tres Cantos, Madrid), yielding a total of 95 determinations. Two main initial conditions must be honored: the specimens must be larger than 60 mL in volume, and their surface/volume ratio must range between 0.08 and 0.20 mm−1 (regular cubes with sides between 30 and 75 mm). In the laboratory, the specimens are first dried in an oven at 70 ± 5 °C until a constant mass (mD) is achieved; the difference between two consecutive weighing carried out after 24 h must be lower than 0.1% of the initial mass (this mass will also allow obtaining the dry bulk density). Then, the specimens are subjected to a vacuum of 2.0 (±0.7) kPa for 2 (±0.2) h. Distilled water is slowly added without losing the vacuum, and the samples are submerged more than 15 min. The water temperature should be 20 ± 5 °C. Atmospheric pressure is slowly restored and the samples are left 24 (±2) additional hours in water. Then, we proceed to the weighing of the specimens immersed in a hydrostatic balance (mH) and subsequently dried with a damp cloth, and its saturated mass is determined (mS). See additional details and instrumentation in Rubio et al. [1]).
Then the apparent (saturated) density is (in kg/m3):
ρ = ρW × mD/([mS − mH])
where ρW is the density of distilled water at 20 °C: 988 kg/m3.
Type 2 samples (more than 80 mini-blocks) were measured in the University of Barcelona laboratories applying the Archimedes principle (without a previous drying and without using paraffin); part of this set was previously measured in the IGME laboratories in Madrid (as type 1 samples). In type 3 AMS samples (111 specimens from MAR), density was estimated in three different ways: using the Archimedes principle (samples with and without paraffin) and estimation of the rock volume with a Vernier caliper on cylindrical samples. Using this third method, 21 additional estimates from MLA were also taken. For that purpose, only regular (cylindrical sections) and complete samples (whole, unbroken, etc.) were used. Any broken, incomplete, irregular or cracked specimen was ruled out for this purpose. Apart from the weighting of the sample (m), the maximum and minimum diameters (Ø) and heights of the specimen (H) were measured with a Vernier caliper. Afterwards both measurements were averaged out (Øm and Hm) and the volume was rapidly calculated: V = π (Øm/2)2·Hm, as well as the density ρ = m/(π [Øm/2]2·Hm). Considering all types together, 310 density determinations were obtained; 69 in MAL, 130 in MLA and 111 in MAR plutons (Table 2).

3.2. Susceptibility and Rock Magnetism Measurements

In total, more than 2600 susceptibility determinations were obtained in 128 studied sites (Table 2). More than 80% were taken directly in 44 outcrops accompanying type-1 samples (21 selected sites from MAL and 23 ones from MLA plutons) using SM20 (1687 measurements) and KT20 (428 measurements) portable susceptometers (by GF Instruments and Terraplus, respectively). On average, more than 50 readings were measured in these outcrops (ca. 10–15 m2), and at least one large hand sample (type 1) was taken approximately in the center of the outcrop surface. Type 2 samples (>200 readings) were measured in a KLY-2 kappabridge located in the Paleomagnetic laboratory of the CCiTUB and Geociencias Barcelona (CSIC). Type 3 samples (paleomagnetic standard) were measured in susceptibility bridges available at the laboratories of the universities of Toulouse and Zaragoza (AGICO kappameters models KLY-2 and KLY-3, respectively). Selected specimens from Toulouse were measured again (1 per selected site) and contrasted to previously published works in the Mont Louis-Andorra granite (Gleizes et al. [34]). In the case of Marimanha, all available stored specimens at the Geotransfer Group (University of Zaragoza) were re-measured in some particular sites (Antolín et al. [46]; Rubio et al. [110]; Loi et al. [111]).
Finally, some comparisons between high (κHF) and low field (κLF) susceptibilities were performed using a MPMS superconducting magnetometer (model 5S by Quantum Design Ltd.) of the University of Zaragoza to determine the para- and ferromagnetic s.l. fractions. Selected cores from MLA pluton were re-measured in the KLY-3 and then subsampled. Those minicores (type 4 samples) were measured in the MPMS instrument under the same conditions as the AGICO instruments (external magnetic field of 0.4 mT, 4 Oe in AC at about 900 Hz) to obtain the low-field value (κLF). Eventually, the same minicores were subjected to two high DC fields (0.9 and 2.5 T), after the total magnetic saturation, to estimate the high-field paramagnetic susceptibility (κHF). Since the κHF only reflects the paramagnetic contribution, meaning κLFκHF, then a dominant paramagnetic contribution can be assumed. Otherwise, if κLF > κHF, then the ferromagnetic contribution cannot be neglected. Additionally, the same collection of mini cores was re-measured in a MVSM (model 3900 by Princeton Measurements Corporation) hosted in the CEREGE paleomagnetic laboratory at Aix-en-Provence while trying to characterize the ferromagnetic hysteresis as well as the high field (up to 1 T) response (paramagnetic slope).

4. Results

4.1. Paramagnetic Contribution to Bulk Susceptibility

Linear correlations of density and magnetic susceptibility have been only proven in a few samples from calc-alkaline (paramagnetic) granites from SE European Variscides (Ameglio et al. [12]). Therefore, the dominant paramagnetic contribution to the susceptibility must be firstly evaluated in our dataset. Standard laboratory measurements of magnetic (bulk) susceptibility are often taken with AGICO susceptometers at room temperature and low magnetic fields (0.4 mT or below, about 10 times larger than the Earth’s magnetic field) and depend on the contribution of all kinds of minerals (diamagnetic, paramagnetic and ferromagnetic). In non-magnetic granites (κ from 10−5 to 10−4 S.I. after Ellwood & Werner [19]) bulk susceptibility is reasonably assumed to be dominated by paramagnetic minerals (Bouchez [30]).
However, some previous comparisons of high-field and low-field susceptibilities in the Bielsa and Millares massifs (Román-Berdiel et al. [84,85], respectively) from the Pyrenees, near the studied plutons (see Figure 1B), attested to the occurrence of significant ferromagnetic fractions in the bulk susceptibility value: up to 20% and 40% in Millares and Bielsa, respectively (but never higher). Conversely, Marimanha data by Antolín et al. [46]), displayed a pure paramagnetic signal (Figure 5a).
Therefore, we run high-field and low-field susceptibility measurements on a set of representative samples (type 4) from the MLA pluton to unambiguously prove or disprove this assumption. The MPMS measurements confirm the pure paramagnetic behavior of the LF signal, since not a single analyzed sample from MLA gave any significant proportion of ferromagnetic components (Figure 5a). Besides, the induced magnetization measurements taken in the MVSM show a pure paramagnetic (linear) behavior and no significant ferromagnetic hysteresis (Figure 5b). Taking into account all of these results, together with the similar distribution of the magnetic susceptibility in the three studied granites (κ < 400 × 10−6 S.I. in 93% of the studied sites; Figure 6), we can extend this conclusion to our entire dataset and a generalized paramagnetic behavior can be assumed.

4.2. Data Distribution

The overall magnetic susceptibility distribution obtained in the analyzed samples (all types together) falls within the entire range of variability of the studied granitic bodies obtained in previous works (Gleizes et al. [34]; Leblanc et al. [45]; Antolín et al. [46], see Figure 2), as was originally designed. In Figure 6, only site means from the studied Pyrenean granites were plotted. The averages of the new susceptibility data (Table 3) mimic the mean from previous papers (Table 1). Therefore, our sampling can be considered representative for the studied bodies. Identically, the density distribution (again, only site means) also falls within the expected range of variability of granitic, syenitic, dioritic and gabbroic rocks (e.g., Dortman [112]; Schön [113]); 2500–2900 kg/m3.
The Henkel’s plotting (1976) of density versus magnetic susceptibility data for each pluton (MAL, MLA and MAR) (only site means) evinces, in the first instance, for noisy petrophysical relationships between both variables (Figure 7) with poor correlation coefficients. At this point, linear regressions are only reliable for the Marimanha pluton (R: 0.81), however, a deeper analysis must be performed and discussed, focusing on the sample types and their reliability.

5. Discussion

5.1. Noise in Magnetic Susceptibility and Density Measurements

Simple box-plots (Figure 8) of both variables and sample types separately attest for mean values relatively robust but also witnesses for several outliers in our dataset. These outliers reflect, in part, the extreme members of the observed variables, but not only. For example, the Maladeta susceptibility data seem to concentrate most of the outliers as well as the most extreme values found in the density distribution. It is worth noticing that the Marimanha pluton does not show any outlier.
The natural variability of the petrophysical properties in a given outcrop (site) as well as the type of samples used to characterized them are key variables to understand the noise in the Henkel’s plot (Figure 7). Magnetic susceptibility within an outcrop (ca. 10–15 m2) can easily spread along one order of magnitude or more. This behavior was already described, to a certain extent, in a homogeneity analysis at the meter and decameter scales performed in AMS sites of granites (Olivier et al. [114]). Density, at this scale, is much less variable in these rocks and hardly exceeds 1% in studies with 7 to 9 determinations per site (Loi et al. [111]; Rubio et al. [110]). In this sense, the total rock volume where both properties are measured is also critical. Type 3 samples are much smaller than type 1 (hand blocks), but since 9–15 standard specimens evenly distributed in the outcrop are measured (Enkin [115]; Pueyo et al. [116,117]), the total 90–150 cm3 rock volume measured at the site scale is comparable with type 1 samples but likely more representative since it covers a larger outcrop surface.
On the other hand, the comparison of susceptibility measurements from different instruments is not necessarily straightforward even in laboratory instruments (Sagnotti et al. [118]). Focusing on portable devices, few works have provided comparative analysis so far (Izquierdo-Llavall et al. [119]; Terrinha et al. [27];). Measuring settings (applied magnetic field and frequency) as well as coil diameters are not always the same and, therefore, measurements from different instruments are not necessarily comparable. This is similar to what happens with density determinations from different methods. For example, the geometric method for determining the sample volume in standard AMS or paleomagnetic samples used to underestimate the rock density. Assumable differences between 1–2%, depending upon the lithology, have been observed in Pyrenean rocks (Pueyo et al. [18]).
Therefore, the number and type of data in a given outcrop, as well as the instrumental and methodological approaches, are critical when attempting to ascertain the reliability and representativeness of the data, and thus they may represent significant sources of noise in the dataset. However, these issues are beyond the scope of this paper and should be tackled in future studies. Meanwhile, the experience gained in this study allows us to propose some best practices in the field and in the laboratory to improve the quality of the petrophysical characterization of granites:
Hand samples (type 1) are always desirable. Ideally, two–three samples of some dm3 each will allow obtaining some regular parallelepipeds in the laboratory (edges about 4–8 cm) that will yield a reliable density estimation and its uncertainty. Six susceptibility measurements (one per face) have to be directly taken in the parallelepiped faces of every sample for the same reason.
Susceptometer coil diameter should be optimized in relation to the sample volume, many portable coil models display diameters about 5–6 cm (ideal for type-1 samples), and smaller coils should be avoided.
Additionally, outcrop measurements should be taken (more than 40 readings) to allow the comparison of the variability between the decametric and centimetric scales, and to characterize the representativeness of the measurements.
Small rock blocks (type 2) may provide reliable couplets of data (density and susceptibility) in every sample, but can only be used if a sufficient number of them (ca. 10 samples) are measured in a given outcrop to guarantee the representativeness of the site mean.
Standard AMS or paleomagnetic samples (type 3) provide a reasonable characterization of both variables if nine or more specimens are measured. Geometric determination of volume is not recommended unless laser or photogrametric techniques are used for the volume estimation; Whiting et al. [120]; (Moret-Fernández et al. [121]).

5.2. Susceptibility and Density Regressions in the Paramagnetic Domain

After filtering the outliers, we have estimated the regressions in the three studied plutons, considering the sample types separately (Table 4 and Figure 9a). Marimanha displays the best definition (only type-3 samples) (R: 0.91) and a full characterization of the uncertainty for both variables (the only case). Mont Louis-Andorra shows a good definition, especially in type-2 samples (R: 0.90), although the three analyzed sample types show relatively good resemblance among the regressions (slopes and intercepts) and all R coefficients are above 0.80. Finally, Maladeta displays the larger discrepancies in all coefficients and the worst regressions (R: 0.54 for type-1 samples), although the regression parameters are very similar to those found in the other two granites (Table 4).
Merging all data together (Figure 9d), the regression yields a reasonable fit: slope 0.458, intercept at 2586 (R: 0.74). The combined observation of all these regressions together (Figure 9e) allows us to be confident to establish a reliable fit since all the parameters are very similar (Table 4 and Figure 9); the intercept (b) ranges between 2548 and 2616 with a mean of 2578 (±23 as standard deviation), the slope (a) between 0.37 and 0.57 mean 0.49 (±0.05) and the squared root of the determination coefficients (R) between 0.54 and 0.91 with a mean of 0.77 (±0.12).
Finally, and focusing on obtaining a unique function needed for the conversion of susceptibility into density data, we estimate the averages in distinct increments (Figure 9f) using all data together (those in Figure 9d, without the outliers in red). A susceptibility window of 50 × 10−6 S.I. units was used for this purpose and all means (couplets of density and susceptibility data) were also characterized by their respective standard deviations.
ρ (kg/m3) = 2566 (kg/m3) + 0.541 κ (10−6 S.I.) × (R:0.97)

5.3. Conversion of Magnetic Susceptibility Data into Density Data

This final regression has been applied to the Mont Louis-Andorra, Maladeta and Marimanha AMS databases (Gleizes et al. [34]; Leblanc et al. [45]; Antolín et al. [46] respectively) to obtain the density data distribution shown as histograms (Figure 10). The three bodies display very similar mean densities (Table 5) around 2675 kg/m3, slightly higher than our own estimations at the laboratory based on more than 300 determinations (Table 2). This apparent underestimation is partially caused by the dominance of more basic (dense) facies at the pluton scale. MAL and MLA show unimodal and sharp distributions (each body is characterized by more than 250 AMS sites) while MAR (only based in 62 sites and much smaller outcropping size) seems to display a bimodal distribution. The transformed density data have also allowed us to build the density maps (Figure 11). The derived new density grids could be used for the modeling of the gravimetric data in future studies.

5.4. Relationships between Geochemical, Petrological and Petrophysical Properties

A negative correlation between density and magnetic susceptibility with respect to the % in wt of SiO2 is remarkable (Figure 12). This relationship is well known in magnetic granites (κ > 500 × 10−6 S.I.) (Bourne [29]; Dentith et al. [17]) and, to a lesser extent, in paramagnetic ones (Villaseca et al. [20]). If we attend to the classification of Gleizes et al. [122]), the correlation extends between the felsic facies (leucogranites) up to the more basic ones (quartz diorites and gabbros). Conversely, this correlation is positive if we focus on other oxides (Al2O3, Fe2O3 or Ti02). The linear increase of the magnetic susceptibility (and density) with the iron content will only work in the paramagnetic domain (κ < 500 × 10−6 S.I.).
On the other hand, the density contour map in the studied granites (Figure 11) mimics the susceptibility one since the relation between both properties is linear. Therefore, the concentric zoning recognized mainly in Maladeta (both western and eastern sectors), as well as in Marimanha plutons, can also be recognized in the density map. The resemblance between the petrological features (Charlet [100,101]) and the susceptibility patterns found in map-view was already recognized in the paper by Leblanc et al. [45]) and the one by Antolín et al. [46]) in Marimanha, and we can extend this comparison to the new density maps. Leucogranites display densities below ca. 2600 kg/m3, monzogranites between 2600 and ca 2690 kg/m3, granodiorites up to ca. 2775 kg/m3 and quartz diorites and more basic facies (gabbros) display larger density values (see also Figure 13). These ranges of variability for granitic (s.l.) facies are in total agreement with classic petrophysical classifications (Dortman [111]; Olhoeft and Johnson [123] Schön [113]).

5.5. The Pyrenean Petrophysical Database of Granites

Finally, we can tentatively extend the robust linear regression found in the Maladeta, Marimanha and Mont Louis-Andorra granites (Equation (3)) to the entire AMS Pyrenean database (Porquet et al. [43] and references therein). It is worth mentioning that, ideally, this regression should be established for every single pluton. However, the similarity of the linear regression found among the three studied plutons (Table 4) together with the similarity of the susceptibility distributions of the 22 granitic Pyrenean bodies (Porquet et al. [44] and reprocessed data from Table 6) allow us to be confident about this preliminary exercise that should be improved in future studies.
We first built a combined histogram merging all data together (more than 2200 sites and 10,000 susceptibility determinations, Figure 13 and Figure 14). There, the density and magnetic susceptibility variables are related together. Besides, we also included the petrological facies boundaries defined by Gleizes et al. [34]). Monzodioritic and Granodioritic facies prevail in the distribution, and the overall mean density value of this vast dataset is: 2680.1 kg/m3 (±112.1), a magnitude that could be used in the regional modeling of gravimetric data in the Pyrenean range. Locally, the density of other bodies may change between 2651.3 kg/m3 (minimum value in Eaux Chaudes) and 2742.5 kg/m3 (maximum mean density found in Borderes).
Interestingly, a clear pattern seems to exist between the mean magnetic susceptibility (and the equivalent density) of the granitic bodies and their outcropping surface (Figure 13b). Although an accurate relationship between the outcropping surface and the granitic size (volume) cannot be established, we can cautiously assume this qualitative relationship. After reaching a critical size around 40 km2, the mean value of the petrophysical variables remains relatively constant at 200 × 10−6 S.I. and 2680 kg/m3. Most smaller bodies display larger mean densities, larger magnetic susceptibilities and the proven occurrence of significant ferromagnetic fractions, as in the Bielsa or Millares plutons (Figure 5; Román-Berdiel et al. [84,85], respectively). This critical threshold composition of these plutonic rocks could be related to the source region of magma generation. Large volumes always correspond to very homogeneous granodiorites with very small variations in the content of ferromagnesian minerals and, consequently, of bulk susceptibility of the rock. These rocks generate by partial melting of the lower crust, while more basic facies may derive from deeper mantle and less differentiated melts. This observation should be studied in future works.

6. Conclusions

The comprehensive study of petrophysical properties (density and magnetic susceptibility) in 128 sites from three plutonic bodies of the Central Pyrenees (310 density measurements and more than 2600 susceptibility ones) has allowed for establishing the following conclusions:
A sufficient number of samples are needed to characterize the outcrop variability since magnetic susceptibility may significantly vary (1–2 orders of magnitude) in a given outcrop, in this sense:
2–3 Type-1 samples (hand blocks) yielding 6–7 samples (preferably cube sides between 5–7 cm) are necessary to fully characterize the outcrop variability (both in density and in magnetic susceptibility).
Alternatively, 9 or more mini blocks (type-2), or, preferably, standard paleomagnetic or AMS (type-3) samples, are also able to reliably describe the outcrop variability.
The plotting density and magnetic susceptibility together (Henkel, 1994) allows for observing similar linear relationships between these variables in the three studied Pyrenean granites.
This relationship: ρ (kg/m3) = 2566 (kg/m3) + 0.541 κ (10−6 S.I.) (R:0.97) falls within the paramagnetic domain in non-magnetic granites (κ < 500 × 10−6 S.I.), where the iron is mostly fractioned in iron-bearing phyllosilicates and the occurrence of magnetite is negligible (or at least its contribution to the bulk susceptibility), as was proven in our dataset.
This regression (which may be different in other bodies) allows for transforming magnetic susceptibility data into density data and, thus, can be very helpful to improve the petrophysical knowledge during the modeling of potential field signals (gravity and magnetism).
The tentative conversion of the entire granitic database of the Pyrenees (>10,000 susceptibility determinations from 22 granitic bodies) yields an average density value of: 2680 kg/m3 (+/−112 kg/m3) that could be used in the regional modeling of gravimetric signals in the Pyrenean range.
Considering the large amount of AMS studies of granites plutons worldwide, usually characterized by dense sampling nets, a calibrated and specific function in every case can be useful to derive reliable and robust databases of density and to improve the modelling of potential fields geophysical data (gravimetric and magnetic). Besides, field, quick, and cost-effective susceptibility measurements performed with portable kappameters may help in densifying and homogenizing these datasets.

Author Contributions

Conceptualization, E.L.P., T.R.-B., J.L.B., F.M.R., J.M. and J.L.G.-L.; Data curation, E.B., C.A., F.L., R.S., P.C. (Pilar Clariana), A.M., N.B., N.S., A.G. (Ana Gimeno), E.F.d.A., C.R.-M., A.G. (Arturo García) and J.L.G.-L.; Formal analysis, P.C. (Pablo Calvín), E.B., C.A., F.L., R.S., P.C. (Pilar Clariana), N.B., N.S., A.G. (Ana Gimeno), E.F.d.A., C.R.-M., A.G. (Arturo García) and J.L.G.-L.; Funding acquisition, C.A. and J.M.; Methodology, E.L.P., E.B., N.S. and F.M.R.; Validation, F.L., A.M., F.M.R. and A.M.C.; Visualization, P.C. (Pilar Clariana) and P.C. (Pablo Calvín); Writing—original draft, E.L.P.; Writing—review & editing, E.L.P., T.R.-B. and P.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financed by the projects GeoPiri3D (CGL2017-84901-C2-2-P), UKRIA4D (PID2019-104693GB-I00/CTA) and IMAGYN (PID2020-114273GB-C22) from the Spanish Ministry of Science (funded by MCIN/AEI/10.13039/501100011033 and “ERDF A way of making Europe”). This work is methodologically related also with the 3DGeoEU project financed by the European Commission under the ERANET Cofound action GeoERA (Grant No.: 731166). The GeoAp Research group from the Aragonian Government is also acknowledged. We are also in debt to the staff of the Petrophysical Laboratory (IGME, Tres Cantos) and to the Geophysics technicians (José Mª Llorente and Agustín González). The help of the Rock Magnetism Laboratory of the Centre Européen de Recherche et d’Enseignement de Géosciences de l’Environnement (CEREGE) in Aix en Provence-Marseille is also acknowledged. P.C. acknowledges funding from PTA2017-14779-I and FJC2019-041058-I (AEI-Spain) contracts. E.B. thanks the Geomodels Research Institute of the University of Barcelona and GGAC-2017SGR596 (Generalitat de Catalunya).

Data Availability Statement

Petrophysical data shown in this paper can be consulted in IGME data repositories, in particular in the SIGEOF database (https://info.igme.es/SIGEOF/ (accessed on 5 January 2022)) that follows FAIR principles. There, the petrophysical window tab must be selected.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Joly, A.; Porwal, A.; McCuaig, T.C. Exploration targeting for orogenic gold deposits in the Granites-Tanami Orogen: Mineral system analysis, targeting model and prospectivity analysis. Ore. Geol. Rev. 2012, 48, 349–383. [Google Scholar] [CrossRef]
  2. Dentith, M.; Mudge, S.T. Geophysics for the Mineral Exploration Geoscientist; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  3. Wang, J. High-level radioactive waste disposal in China: Update 2010. J. Rock Mech. Geotech. Eng. 2010, 2, 1–11. [Google Scholar]
  4. Genter, A.; Evans, K.; Cuenot, N.; Fritsch, D.; Sanjuan, B. Contribution of the exploration of deep crystalline fractured reservoir of Soultz to the knowledge of enhanced geothermal systems (EGS). Comptes Rendus Geosci. 2010, 342, 502–516. [Google Scholar] [CrossRef]
  5. Huenges, E.; Ledru, P. (Eds.) Geothermal Energy Systems: Exploration, Development, and Utilization; John Wiley & Sons: New York, NY, USA, 2011. [Google Scholar]
  6. Moore, J.; McLennan, J.; Allis, R.; Pankow, K.; Simmons, S.; Podgorney, R.; Wannamaker, P.; Bartley, J.; Jones, C.; Rickard, W. The Utah Frontier Observatory for Research in Geothermal Energy (FORGE): An International Laboratory for Enhanced Geothermal System Technology Development. In Proceedings of the 44th Workshop on Geothermal Reservoir Engineering, Stanford, CA, USA, 11–13 February 2019. [Google Scholar]
  7. Zhang, Y.; Zhao, G.-F. A global review of deep geothermal energy exploration: From a view of rock mechanics and engineering. Géoméch. Geophys. Geo-Energy Geo-Resour. 2019, 6, 4. [Google Scholar] [CrossRef]
  8. Bott, M.H.P.; Day, A.A.; Masson-Smith, D. The geological interpretation of gravity and magnetic surveys in Devon and Cornwall. Philosophical Transactions of the Royal Society of London. Ser. A Math. Phys. Sci. 1958, 251, 161–191. [Google Scholar]
  9. Bott, M.H.P.; Smithson, S.B. Gravity Investigations of Subsurface Shape and Mass Distributions of Granite Batholiths. GSA Bull. 1967, 78, 859–878. [Google Scholar] [CrossRef]
  10. Henkel, H. Studies of density and magnetic properties of rocks from Northern Sweden. Pure Appl. Geophys. 1976, 114, 235–249. [Google Scholar] [CrossRef]
  11. Vigneresse, J.L. Use and misuse of geophysical data to determine the shape at depth of granitic intrusions. Geol. J. 1990, 25, 249–260. [Google Scholar] [CrossRef]
  12. Améglio, L.; Vigneresse, J.L.; Bouchez, J.L. Granite Pluton Geometry and Emplacement Mode Inferred from Combined Fabric and Gravity Data. In Granite: From Segregation of Melt to Emplacement Fabrics; Springer: Berlin/Heidelberg, Germany, 1997; pp. 199–214. [Google Scholar] [CrossRef]
  13. Cruz, C.; Sant’Ovaia, H.; Raposo, M.I.B.; Lourenço, J.M.; Almeida, F.; Noronha, F. Unraveling the emplacement history of a Portuguese post-tectonic Variscan pluton using magnetic fabrics and gravimetry. J. Struct. Geol. 2021, 153, 104470. [Google Scholar] [CrossRef]
  14. Henkel, H. Petrophysical properties (density and magnetization) of rocks from the northern part of the Baltic Shield. Tectonophysics 1991, 192, 1–19. [Google Scholar] [CrossRef]
  15. Henkel, H. Standard diagrams of magnetic properties and density—A tool for understanding magnetic petrology. J. Appl. Geophys. 1994, 32, 43–53. [Google Scholar] [CrossRef]
  16. Enkin, R.J.; Hamilton, T.S.; Morris, W.A. The Henkel Petrophysical Plot: Mineralogy and Lithology From Physical Properties. Geochem. Geophys. Geosyst. 2020, 21, e2019GC008818. [Google Scholar] [CrossRef]
  17. Dentith, M.; Enkin, R.J.; Morris, W.; Adams, C.; Bourne, B. Petrophysics and mineral exploration: A workflow for data analysis and a new interpretation framework. Geophys. Prospect. 2019, 68, 178–199. [Google Scholar] [CrossRef] [Green Version]
  18. Pueyo, E.L.; Ayala, C.; Izquierdo-Llavall, E.; Rubio, F.M.; Santolaria, P.; Clariana, P.; Soto, R.; Müller, C.O.; Rey-Moral, C.; Zehner, B.; et al. Optimized 3D Reconstruction Workflow Based on Gravimetric, Structural and Petrophysical Data. Deliverable 6.4. Project 3dgeo-EU, Geoera- 3D Geomodeling for Europe, Project Number GeoE.171.005. Report. 2021. 264p. Available online: https://geoera.eu/wp-content/uploads/2022/01/3DGEO-EU_D8.5_Summary-of-project-work-and-results.pdf (accessed on 14 January 2022).
  19. Ellwood, B.B.; Wenner, D.B. Correlation of magnetic susceptibility with 18O/16O data in late orogenic granites of the southern Appalachian Piedmont. Earth Planet. Sci. Lett. 1981, 54, 200–202. [Google Scholar] [CrossRef]
  20. Villaseca, C.; Ruiz-Martínez, V.C.; Pérez-Soba, C. Magnetic susceptibility of Variscan granite-types of the Spanish Central System and the redox state of magma. Geol. Acta 2017, 15, 379–394. [Google Scholar]
  21. Kanaya, H.; Ishihara, S. Regional variation of magnetic susceptibility of the granitic rocks in japan. J. Jpn. Assoc. Miner. Pet. Econ. Geol. 1973, 68, 211–224. [Google Scholar] [CrossRef]
  22. Ishihara, S. The magnetite-series and ilmenite-series granitic rocks. Min. Geol. 1977, 27, 293–305. [Google Scholar]
  23. Ishihara, S. The Granitoid Series and Mineralization. In 75th Anniversary Volume: Economic Geology; Society of Economic Geologists: Littleton, CO, USA, 1981. [Google Scholar]
  24. Chappell, B.W.; White, A.J. Two contrasting granite types: 25 years later. Aust. J. Earth Sci. 2001, 48, 489–499. [Google Scholar] [CrossRef]
  25. Hernandez, F.M.; Hirt, A. The anisotropy of magnetic susceptibility in biotite, muscovite and chlorite single crystals. Tectonophysics 2003, 367, 13–28. [Google Scholar] [CrossRef]
  26. Elming, S.-A. Density and magnetic properties of rocks in the Caledonides of Jämtland, Sweden. Geol. Föreningen I Stockh. Förhandlingar. 1980, 102, 439–453. [Google Scholar] [CrossRef]
  27. Terrinha, P.; Pueyo, E.L.; Aranguren, A.; Kullberg, J.C.; Kullberg, M.C.; Casas-Sainz, A.; Azevedo, M.D.R. Gravimetric and magnetic fabric study of the Sintra Igneous complex: Laccolith-plug emplacement in the Western Iberian passive margin. Geol. Rundsch. 2017, 107, 1807–1833. [Google Scholar] [CrossRef]
  28. Criss, R.E.; Champion, D.E. Magnetic properties of granitic rocks from the southern half of the Idaho Batholith: Influences of hydrothermal alteration and implications for aeromagnetic interpretation. J. Geophys. Res. Earth Surf. 1984, 89, 7061–7076. [Google Scholar] [CrossRef]
  29. Bourne, J.H. Use of magnetic susceptibility, density, and modal mineral data as a guide to the composition of granitic plutons. Math. Geol. 1993, 25, 357–375. [Google Scholar] [CrossRef]
  30. Bouchez, J.L. Granite is Never Isotropic: An Introduction to AMS Studies of Granitic Rocks. In Granite: From Segregation of Melt to Emplacement Fabrics; Springer: Berlin/Heidelberg, Germany, 1997; pp. 95–112. [Google Scholar]
  31. Bouchez, J.L. Anisotropie de susceptibilité magnétique et fabrique des granites. Comptes Rendus De L’académie Des Sci.-Ser. IIA-Earth Planet. Sci. 2000, 330, 1–14. [Google Scholar] [CrossRef]
  32. Gleizes, G. Structure Des Granites Hercyniens Des Pyrénées de Mont Louis-Andorre à la Maladeta. Ph.D. Thesis, Université de Toulouse III, Toulouse, France, 1990. [Google Scholar]
  33. Bouchez, J.L.; Gleizes, G.; Djouadi, T.; Rochette, P. Microstructure and magnetic susceptibility applied to emplacement kinematics of granites: The example of the foix pluton (French pyrenees). Tectonophysics 1990, 184, 157–171. [Google Scholar] [CrossRef]
  34. Gleizes, G.; Nédélec, A.; Bouchez, J.-L.; Autran, A.; Rochette, P. Magnetic susceptibility of the Mont-Louis andorra ilmenite-type granite (Pyrenees): A new tool for the petrographic characterization and regional mapping of zoned granite plutons. J. Geophys. Res. Earth Surf. 1993, 98, 4317–4331. [Google Scholar] [CrossRef]
  35. Rochette, P. Magnetic susceptibility of the rock matrix related to magnetic fabric studies. J. Struct. Geol. 1987, 9, 1015–1020. [Google Scholar] [CrossRef]
  36. Rochette, P.; Jackson, M.; Aubourg, C. Rock magnetism and the interpretation of anisotropy of magnetic susceptibility. Rev. Geophys. 1992, 30, 209–226. [Google Scholar] [CrossRef]
  37. Román-Berdiel, T.; Pueyo, E.L.; Casas-Sainz, A. Granite emplacement during contemporary shortening and normal faulting: Structural and magnetic study of the Veiga Massif (NW Spain). J. Struct. Geol. 1995, 17, 1689–1706. [Google Scholar] [CrossRef]
  38. Aranguren, A.; Tubia, J.; Bouchez, J.L.; Vigneresse, J.-L. The Guitiriz granite, Variscan belt of northern Spain: Extension-controlled emplacement of magma during tectonic escape. Earth Planet. Sci. Lett. 1996, 139, 165–176. [Google Scholar] [CrossRef]
  39. Trindade, R.I.; Raposo, M.B.; Ernesto, M.; Siqueira, R. Magnetic susceptibility and partial anhysteretic remanence anisotropies in the magnetite-bearing granite pluton of Tourão, NE Brazil. Tectonophysics 1999, 314, 443–468. [Google Scholar] [CrossRef]
  40. Sant’Ovaia, H.; Bouchez, J.L.; Noronha, F.; Leblanc, D.; Vigneresse, J.L. Composite-laccolith emplacement of the post-tectonic Vila Pouca de Aguiar granite pluton (northern Portugal): A combined AMS and gravity study. Trans. R. Soc. Edinb. Earth Sci. 2000, 91, 123–137. [Google Scholar]
  41. Ferre, E.; Gleizes, G.; Caby, R. Obliquely convergent tectonics and granite emplacement in the Trans-Saharan belt of Eastern Nigeria: A synthesis. Precambrian Res. 2002, 114, 199–219. [Google Scholar] [CrossRef]
  42. Kratinová, Z.; Schulmann, K.; Edel, J.-B.; Jezek, J.; Schaltegger, U. Model of successive granite sheet emplacement in transtensional setting: Integrated microstructural and anisotropy of magnetic susceptibility study. Tectonics 2007, 26, 1–7. [Google Scholar] [CrossRef]
  43. Joly, A.; Faure, M.; Martelet, G.; Chen, Y. Gravity inversion, AMS and geochronological investigations of syntectonic granitic plutons in the southern part of the Variscan French Massif Central. J. Struct. Geol. 2009, 31, 421–443. [Google Scholar] [CrossRef] [Green Version]
  44. Porquet, M.; Pueyo, E.L.; Román-Berdiel, T.; Olivier, P.; Longares, L.A.; Cuevas, J.; Ramajo, J.; Antolín, B.; Aranguren, A.; Auréjac, J.B.; et al. Anisotropy of magnetic susceptibility of the Pyrenean granites. J. Maps 2017, 13, 438–448. [Google Scholar] [CrossRef]
  45. Leblanc, D.; Gleizes, G.; Lespinasse, P.; Olivier, P.; Bouchez, J.-L. The maladeta granite polydiapir, Spanish Pyrenees: A detailed magnetostructural study. J. Struct. Geol. 1994, 16, 223–235. [Google Scholar] [CrossRef]
  46. Antolín-Tomás, B.; Román-Berdiel, T.; Casas-Sainz, A.; Gil-Peña, I.; Oliva, B.; Soto, R. Structural and magnetic fabric study of the Marimanha granite (Axial Zone of the Pyrenees). Geol. Rundsch. 2007, 98, 427–441. [Google Scholar] [CrossRef]
  47. Gleizes, G.; Bouchez, J.L. Le granite de Mont-Louis (Zone axiale des Pyrénées): Anisotropie magnétique, structures et microstructures. Comptes Rendus De L’académie Des Sci. Paris 1989, 309, 1075–1082. [Google Scholar]
  48. Antolín-Tomás, B. Estudio Estructural Del Plutón De Marimanha (Pirineo Central). Posgrado de Iniciación a La Investigación en Geología; Universidad de Zaragoza: Zaragoza, Spain, 2006; p. 128. [Google Scholar]
  49. Matte, P. Accretionary history and crustal evolution of the Variscan belt in Western Europe. Tectonophysics 1991, 196, 309–337. [Google Scholar] [CrossRef]
  50. Martínez-Catalán, J.R. A non-cylindrical model for the northwest Iberian allochthonous terranes and their equivalents in the Hercynian belt os Western Europe. Tectonophysics 1990, 179, 253–272. [Google Scholar] [CrossRef]
  51. Pastor-Galán, D.; Groenewegen, T.; Brouwer, D.; Krijgsman, W.; Dekkers, M. One or two oroclines in the Variscan orogen of Iberia? Implications for Pangea amalgamation. Geology 2015, 43, 527–530. [Google Scholar] [CrossRef] [Green Version]
  52. Matte, P. La chaîne varisque parmi les chaînes paléozoïques periatlantiques, modèle d’évolution et position des grands blocs continentaux au Permo-Carbonifère. Bull. Soc. Géol. Fr. 1986, 2, 9–24. [Google Scholar] [CrossRef]
  53. Muñoz, J.A. Evolution of a continental collision belt: ECORS-Pyrenees crustal balanced cross-section. In Thrust Tectonics; McClay, K.R., Ed.; Chapman and Hall: London, UK, 1992; pp. 235–246. [Google Scholar]
  54. Martínez-Peña, M.B.; Casas-Sainz, A.M. Cretaceous-tertiary tectonic inversion of the Cotiella Basin (southern Pyrenees, Spain). Int. J. Earth Sci. 2003, 92, 99–113. [Google Scholar] [CrossRef]
  55. Casas, A.M.; Oliva-Urcia, B.; Román-Berdiel, T.; Pueyo, E. Basement deformation: Tertiary folding and fracturing of the Variscan Bielsa granite (Axial zone, central Pyrenees). Geodin. Acta 2003, 16, 99–117. [Google Scholar] [CrossRef]
  56. Cochelin, B.; Lemirre, B.; Denèle, Y.; Blanquat, M.D.S.; Lahfid, A.; Duchene, S. Structural inheritance in the Central Pyrenees: The Variscan to Alpine tectonometamorphic evolution of the Axial Zone. J. Geol. Soc. 2017, 175, 336–351. [Google Scholar] [CrossRef]
  57. Druguet, E. Development of high thermal gradients by coeval transpression and magmatism during the Variscan orogeny: Insights from the Cap de Creus (Eastern Pyrenees). Tectonophysics 2001, 332, 275–293. [Google Scholar] [CrossRef]
  58. García-Sansegundo, J.; Poblet, J.; Alonso, J.L.; Clariana, P. Hinterland–Foreland Zonation of the Variscan Orogen in the Central Pyrenees: Comparison with the Northern Part of the Iberian Variscan Massif. In Kinematic Evolution and Structural Styles of Fold-and-Thrust Belts; Poblet, J., Lisle, R.J., Eds.; Special Publications; Geological Society: London, UK, 2011; Volume 349, pp. 169–184. [Google Scholar]
  59. Casas, J.M.; Clariana, P.; García-Sansegundo, J.; Margalef, A. The Pyrenees. In The Geology of Iberia: A Geodynamic Approach. Regional Geology Reviews; Quesada, C., Oliveira, J., Eds.; Springer: Cham, Switzerland, 2019; pp. 335–337. [Google Scholar] [CrossRef]
  60. Soula, J.; Debat, P.; Deramond, J.; Pouget, P. A dynamic model of the structural evolution of the Hercynian Pyrenees. Tectonophysics 1986, 129, 29–51. [Google Scholar] [CrossRef]
  61. Carreras, J.; Capella, I. Tectonic levels in the Paleozoic basement of the Pyrenees: A review and a new interpretation. J. Struct. Geol. 1994, 16, 1509–1524. [Google Scholar] [CrossRef]
  62. Gutierrez-Medina, M. Nuevas aportaciones al conocimiento de la estructura varisca y alpina de la lámina cabalgante de Bono, Zona Axial, Pirineos Centrales. Trab. De Geol. Univ. De Oviedo 2007, 27, 159–177. [Google Scholar]
  63. Clariana, P.; García-Sansegundo, J. Variscan structure in the Eastern part of the Pallaresa massif Axial Zone of the Pyrenees (NW Andorra). Tectonic implications. Bull. Soc. Géol. Fr. 2009, 180, 501–511. [Google Scholar] [CrossRef]
  64. Soula, J.-C. Characteristics and mode of emplacement of gneiss domes and plutonic domes in central-eastern Pyrenees. J. Struct. Geol. 1982, 4, 313–342. [Google Scholar] [CrossRef]
  65. Autran, A. Relations Des Cristallisations Métamorphiques Avec Les Phases De Déformation Successives. Métamorphisme Hercynien Des Métapélites Dans Le Massif Du Canigou-Caranca. Chapitre 10. Métamorphisme Hercynien. In Synthèse Géologique Et Géophysique Des Pyrénées; Barnolas, A., Chiron, J.C., Guérangé, B., Eds.; Introduction, Géophysique, Cycle Hercynien; BRGM: Orleans, France; ITGE: Madrid, Spain, 1996; Volumn 1, pp. 547–548. [Google Scholar]
  66. Guitard, G. Phases de plissement dans les terrains métamorphiques de la zone axiale pyrénéenne du Canigou durant l’orogénèse hercynienne. CRAS Paris 1965, 265, 1357–1360. [Google Scholar]
  67. Zwart, H.J. The Geology of the Central Pyrenees. Leidse Geol. Meded 1979, 50, 1–74. [Google Scholar]
  68. Zwart, H. The variscan geology of the Pyrenees. Tectonophysics 1986, 129, 9–27. [Google Scholar] [CrossRef]
  69. Liesa, M. Relations of Hercynian metamorphism with magmatism and deformation in the Eastern Pyrenees. Implications for Hercynian evolution. Geol. Mijnb. 1993, 72, 295–304. [Google Scholar]
  70. Aguilar, C. P–T–t–d Constrains on the Late Variscan Evolution of the Eastern Pyrenees. Ph.D. Thesis, Universitat de Barcelona, Barcelona, Spain, 2013. [Google Scholar]
  71. Majesté-Menjoulas, C. L’unité paleozoïque de Bachebirou-Chinipro, temoin d’une tectonique tangentielle varisque dans les Pyrénées Centrales. CRAS 1982, 294, 145–150. [Google Scholar]
  72. Raymond, D. Tectonique tangentielle varisque dans le Paléozoique supérieur de l’Est des Pyrénées françaises; l’exemple du Pays de Sault (Nord du granite de Quérigut, Aude et Ariège) et des régions voisines. Bull. Soc. Géol. Fr. 1986, 2, 479–485. [Google Scholar] [CrossRef]
  73. Bodin, J.; Ledru, P. Nappes hercyniennes précoces à matériel dévonien hétéropique dans les Pyrénées ariégeoises. CRAS Paris 1986, 302, 969–974. [Google Scholar]
  74. Losantos, M.; Palau, J.; Sanz-López, J. Considerations about Hercynian thrusting in the Marimanya massif (central Pyrenees). Tectonophysics 1986, 129, 71–79. [Google Scholar] [CrossRef]
  75. García-Sansegundo, J. Estratigrafia Y Estructura De La Zona Axial Pirenaica En La Transversal Del Valle De Aran Y Alta Ribagorza. Ph.D. Thesis, Universidad de Oviedo, Oviedo, Spain, 1991. [Google Scholar]
  76. Clariana, P. Estratigrafía, Estructura y Su Relación Con El Metamorfismo De La Zona Axial Pirenaica En La Transversal Del Noroeste De Andorra y Comarcas Del Pallars Sobirá Y El Alt Urgell (Lleida). Ph.D. Thesis, Universidad de Zaragoza, Zaragoza, Spain, 2015; 190p. [Google Scholar]
  77. Margalef, A. Estudi Estructural i Estratigràfic Del Sud d’Andorra. Ph.D. Thesis, Universitat de Barcelona, Barcelona, Spain, 2015. [Google Scholar]
  78. Leblanc, D.; Gleizes, G.; Roux, L.; Bouchez, J. Variscan dextral transpression in the French Pyrenees: New data from the Pic des Trois-Seigneurs granodiorite and its country rocks. Tectonophysics 1996, 261, 331–345. [Google Scholar] [CrossRef]
  79. Evans, N.G.; Gleizes, G.; Leblanc, D.; Bouchez, J.L. Syntectonic emplacement of the Maladeta granite (Pyrenees) deduced from relationships between Hercynian deformation and contact metamorphism. J. Geol. Soc. 1998, 155, 209–216. [Google Scholar] [CrossRef]
  80. Gleizes, G.; Leblanc, D.; Bouchez, J.L. The main phase of the Hercynian orogeny in the Pyrenees is a dextral transpression. Geol. Soc. Lond. Spéc. Publ. 1998, 135, 267–273. [Google Scholar] [CrossRef]
  81. Gleizes, G.; Leblanc, D.; Olivier, P.; Bouchez, J. Strain partitioning in a pluton during emplacement in transpressional regime: The example of the Néouvielle granite (Pyrenees). Geol. Rundsch. 2001, 90, 325–340. [Google Scholar] [CrossRef]
  82. Gleizes, G.; Crevon, G.; Asrat, A.; Barbey, P. Structure, age and mode of emplacement of the Hercynian Bordères-Louron pluton (Central Pyrenees, France). Int. J. Earth Sci. 2006, 95, 1039–1052. [Google Scholar] [CrossRef]
  83. Olivier, P.; Améglio, L.; Richen, H.; Vadeboin, F. Emplacement of the Aya Variscan granitic pluton (Basque Pyrenees) in a dextral transcurrent regime inferred from a combined magneto-structural and gravimetric study. J. Geol. Soc. 1999, 156, 991–1002. [Google Scholar] [CrossRef]
  84. Román-Berdiel, T.; Casas, A.; Oliva-Urcia, B.; Pueyo, E.L.; Rillo, C. The main Variscan deformation event in the Pyrenees: New data from the structural study of the Bielsa granite. J. Struct. Geol. 2004, 26, 659–677. [Google Scholar] [CrossRef]
  85. Román-Berdiel, T.; Casas, A.M.; Oliva-Urcia, B.; Pueyo, E.L.; Liesa, C.; Soto, R. The Variscan Millares granite (central Pyrenees): Pluton emplacement in a T fracture of a dextral shear zone. Geodin. Acta 2006, 19, 197–211. [Google Scholar] [CrossRef]
  86. Auréjac, J.B.; Gleizes, G.; Diot, H.; Bouchez, J.L. Le complexe granitique de Quérigut (Pyrénées, France) ré-examiné par la technique de l’ASM: Un pluton syntectonique de la transpression dextre hercynienne. Bull. Soc. Géol. Fr. 2004, 175, 157–174. [Google Scholar]
  87. Romer, R.; Soler, A. U-Pb age and lead isotopic characterization of Au-bearing skarn related to the Andorra granite (central Pyrenees, Spain). Miner. Depos. 1995, 30, 374–383. [Google Scholar] [CrossRef]
  88. Paquette, J.L.; Gleizes, G.; Leblanc, D.; Bouchez, J.L. Le granite de Bassiès (Pyrénées): Un pluton syntectonique d’âge Westphalien. Géochronologie U–Pb sur zircons. CRAS Paris 1997, 324, 387–392. [Google Scholar]
  89. Guerrot, C. Résultats de datation U-Pb par dissolution sur zircons pour le granite de Cauterets, Pyrénées. Rapport inédit, SMN/PEA/ISO 146/98 CG/NB, BRGM. In Notice Explicative, Carte géol. France (1/50 000), Feuille Gavarnie (1082); Majesté-Menjoulas, C., Debon, F., Barrère, P., Eds.; BRGM: Orléans, France, 1998; p. 4. [Google Scholar]
  90. Guerrot, C. Datation du pluton des Eaux-Chaudes. In Notice Explicative, Carte géol. France (1/50,000), Feuille Laruns-Somport (1069); Ternet, Y., Majesté-Menjoulas, C., Canérot, J., Baudin, T., Cocherie, A., Guerrot, C., Rossi, P., Eds.; BRGM: Orléans, France, 2001; pp. 185–187. [Google Scholar]
  91. Roberts, M.P.; Pin, C.; Clemens, J.D.; Paquette, J.L. Petrogenesis of Mafic to Felsic Plutonic Rock Associations: The Calc-alkaline Quérigut Complex, French Pyrenees. J. Pet. 2000, 46, 809–844. [Google Scholar] [CrossRef]
  92. Maurel, O.; Respaut, J.P.; Monié, P.; Arnaud, N.; Brunel, M. U-Pb emplacement and 40Ar/39Ar cooling ages of the Mont-Louis granite massif (eastern Pyrenees, France). C. R. Geosci. 2004, 336, 1091–1098. [Google Scholar] [CrossRef]
  93. Olivier, P.; Gleizes, G.; Paquette, J. Gneiss domes and granite emplacement in an obliquely convergent regime: New interpretation of the Variscan Agly Massif (Eastern Pyrenees, France). Spec. Pap.-Geol. Soc. Am. 2004, 380, 229–242. [Google Scholar]
  94. Denèle, Y.; Paquette, J.-L.; Olivier, P.; Barbey, P. Permian granites in the Pyrenees: The Aya pluton (Basque Country). Terra Nova 2011, 24, 105–113. [Google Scholar] [CrossRef]
  95. Denèle, Y.; Laumonier, B.; Paquette, J.-L.; Olivier, P.; Gleizes, G.; Barbey, P. Timing of Granite Emplacement, Crustal Flow and Gneiss Dome Formation in the Variscan Segment of the Pyrenees. In The Variscan Orogeny: Extent, Timescale and the Formation of the European Crust; Schulmann, K., Martínez Catalán, J.R., Lardeaux, J.M., Janousek, V., Oggiano, G., Eds.; Special Publications; Geological Society: London, UK, 2014; Volume 405, pp. 265–287. [Google Scholar]
  96. Franke, W. Tectonostratigraphic units in the Variscan belt of central Europe. Geol. Soc. Am. Spec. Pap. 1989, 230, 67–90. [Google Scholar]
  97. Palau i Ramirez, J. El magmatisme calcoalcalí del massís de Marimanya i les mineralitzations As–Au–W associades. Institut cartogràfic de Catalunya. Monogr. Tèchniques 1998, 4, 1–340. [Google Scholar]
  98. Evans, N.G. Deformation in the Aureole of the Maladeta Granodiorite, Spanish Pyrenees: A Record of Pluton Emplacement. In Proceedings of the 14e Réunion Sciences Terre, Toulouse, France, 13 April 1992; p. 59. [Google Scholar]
  99. Vitrac-Michard, A.; Albaède, F.; Dupuis, C.; Tailor, H.P. The genesis of Variscan (Hercynian) plutonic rocks: Inferences from Sr, Pb and O studies on the Maladeta igneous complex, Central Pyrenees (Spain). Contr. Miner. Petrol. 1980, 72, 57–72. [Google Scholar] [CrossRef]
  100. Charlet, J.M. Etude préliminaire du massif granitique de la Maladeta (Pyrénées centrales espagnoles). Ann. Soc. Géol. Nord. Lille 1968, 88, 65–75. [Google Scholar]
  101. Charlet, J.M. Le massif granitique de la Maladeta (Pyrénées centrales espagnoles), synthèse des données géologiques. Ann. Soc. Géol. Belg. 1979, 102, 313–323. [Google Scholar]
  102. Pouget, P.; Lamouroux, C.; Dahmani, A.; Debat, P.; Driouch, Y.; Mercier, A.; Soula, J.C.; Vezat, R. Typologie et mode de mise en place des roches magmatiques dans les Pyrénées hercyniennes. Geol. Rdsch. 1989, 78, 537–554. [Google Scholar] [CrossRef]
  103. Autran, A.; Fonteilles, M.; Guitard, G. Relations entre les intrusions de granitoides, l’anatexie et le métamorphisme général: Cas de la chaîne hercynienne des Pyrénées orientales. Bull. Soc. Géol. Fr. Sér. 1970, 7, 673–731. [Google Scholar] [CrossRef]
  104. Debon, F.; Enrique, P.; Autran, A. Magmatisme Hercynien. In Synthèse Géologique et Géophysique des Pyrénées; Barnolas, A., Chiron, J.C., Eds.; BRGM: Orléans, France; ITGE: Madrid, Spain, 1996; pp. 361–499. [Google Scholar]
  105. Bouchez, J.L.; Gleizes, G. Two-stage deformation of the Mont-Louis-Andorra granite pluton (Variscan Pyrenees) inferred from magnetic susceptibility anisotropy. J. Geol. Soc. 1995, 152, 669–679. [Google Scholar] [CrossRef]
  106. Laumonier, B.; Autran, A.; Barbey, P.; Cheilletz, A.; Baudin, T.; Cocherie, A.; Guerrot, C. Conséquences de l’absence de socle cadomien sur l’âge et la signification des séries pré-varisques (anté- Ordovicien supérieur) du sud de la France (Pyrénées, Montagne Noire). Bull. De La Société Géologique De Fr. 2004, 175, 105–117. [Google Scholar]
  107. Guitard, G.; Laumonier, B.; Autran, A.; Bandet, Y.; Berger, G.M. Notice Explicative, Carte Géologique de France (1/50 000), Feuille de Prades (1095); BRGM: Orléans, France, 1998. [Google Scholar]
  108. Laumonier, B.; Marignac, C.; Kister, P. Polymetamorphism and crustal evolution of the eastern Pyrenees during the Late Carboniferous Variscan orogenesis. BSGF-Earth Sci. Bull. 2010, 181, 411–428. [Google Scholar] [CrossRef]
  109. Clariana, P.; Soto, R.; Ayala, C.; Casas-Sainz, A.M.; Román-Berdiel, T.; Oliva-Urcia, B.; Pueyo, E.L.; Beamud, E.; Rey-Moral, C.; Rubio, F.; et al. Basement and cover architecture in the Central Pyrenees constrained by gravity data. Geol. Rundsch. 2021, 111, 641–658. [Google Scholar] [CrossRef]
  110. Rubio Sánchez de Aguililla, F.; Rey Moral, C.; Ayala, C.; Pueyo, E.L. Bases Metodológicas para la Obtención de Cartografía Gravimétrica. Petrofísica y Modelización Gravimétrica 2D y 3D; Instituto Geológico y Minero de España (IGME, CSIC): Madrid, Spain, 2021; 149p. [Google Scholar]
  111. Loi, F.; Román-Berdiel, M.T.; Casas-Sainz, A.M.; Pueyo, E.L. Correlación Entre Propiedades Petrofísicas Del Granito De Marimanha (Pirineo Axial). In MAGIBER XI-Paleomagnetismo en Espanha e Portugal (Livro de Resumos); Font, E., Beamud, E., Oliva-urcia, B., Pueyo, E.L., Lopez, F.C., Eds.; Universidade de Coimbra: Coimbra, Portugal, 2019; pp. 87–90. ISBN 9789899891470. [Google Scholar]
  112. Dortman, N.B. Fizicheskiye Svoystva Gornykh Porod Ipoleznykh Iskopayemykh (Petrofizika); spravochnik geofizika; USSR Izdat Nedra: Moscow, Russia, 1976; p. 527. [Google Scholar]
  113. Schön, J.H. Physical Properties of Rocks: Fundamentals and Principles of Petrophysics; Elsevier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  114. Olivier, P.; Blanquat, M.D.S.; Gleizes, G.; Leblanc, D. Homogeneity of Granite Fabrics at the Metre and Dekametre Scales. In Granite: From Segregation of Melt to Emplacement Fabrics; Springer: Berlin/Heidelberg, Germany, 1997; pp. 113–127. [Google Scholar]
  115. Enkin, R.J. The Canadian Rock Physical Property Database: First public release. GSC Open File 2018, 8460, 60. [Google Scholar] [CrossRef]
  116. Pueyo, E.L.; Román-Berdiel, T.; Bouchez, J.L.; Casas, A.; Larrasoaña, J.C. Statistical significance of magnetic fabric data in studies of paramagnetic granites. Geol. Soc. Lond. Spéc. Publ. 2004, 238, 395–420. [Google Scholar] [CrossRef]
  117. Pueyo, E.L.; Oliva-Urcia, B.; Sussman, A.J.; Cifelli, F. Palaeomagnetism in Fold and Thrust Belts: Use with caution. In Geological Society of London Special Publication on Palaeomagnetism in Fold and Thrust Belts: New Perspectives; Pueyo, E.L., Cifelli, F., Sussman, A., Oliva-Urcia, B., Eds.; Geological Society: London, UK, 2016; Volume 425, pp. 259–276. [Google Scholar]
  118. Sagnotti, L.; Rochette, P.; Jackson, M.; Vadeboin, F.; Dinarès-Turell, J.; Winkler, A. Inter-laboratory calibration of low-field magnetic and anhysteretic susceptibility measurements. Phys. Earth Planet. Inter. 2003, 138, 25–38. [Google Scholar] [CrossRef]
  119. Izquierdo-Llavall, E.; Ayala, C.; Rodríguez Pintó, A.; Pueyo, E.L.; Oliva-Urcia, B.; Casas, A.M.; Rubio, F.M.; García Crespo, J.; Llorente Delgado, J.M.; González Durán, A.; et al. Informe de la investigación del reservorio potencial de CO2 “Zona de Enlace” (Cordilleras Ibérica y Costero Catalanas), Proyecto ALGECO2 (Fase 2). IGME-IRMC (MITT), technical report 229p., Madrid, Spain, 2014.
  120. Whiting, M.; Salley, S.W.; James, D.K.; Karl, J.W.; Brungard, C.W. Rapid bulk density measurement using mobile device photogrammetry. Soil Sci. Soc. Am. J. 2020, 84, 811–817. [Google Scholar] [CrossRef]
  121. Moret-Fernández, D.; Latorre, B.; Peña, C.; González-Cebollada, C.; López, M.V. Applicability of the photogrammetry technique to determine the volume and the bulk density of small soil aggregates. Soil Res. 2016, 54, 354. [Google Scholar] [CrossRef] [Green Version]
  122. Gleizes, G.; Leblanc, D.; Santana, V.; Olivier, P.; Bouchez, J.L. Sigmoidal structures featuring dextral shear during emplacement of the Hercynian granite complex of Cauterets–Panticosa (Pyrenees). J. Struct. Geol. 1998, 20, 1229–1245. [Google Scholar] [CrossRef]
  123. Olhoeft, G.R.; Johnson, G.R. Densities of Rocks and Minerals. In Practical Handbook of Physical Properties of Rocks and Minerals; Carmichael, R.S., Ed.; Section II; CRC Press: Boca Raton, FL, USA, 1989; pp. 1–38. [Google Scholar]
Figure 2. Susceptibility histograms of site means from previously published AMS studies (data by Gleizes et al. [34] for MLA; Leblanc et al. [44] for MAL; Antolín et al. [46], for MAR, partially reprocessed by Porquet et al. [44]).
Figure 2. Susceptibility histograms of site means from previously published AMS studies (data by Gleizes et al. [34] for MLA; Leblanc et al. [44] for MAL; Antolín et al. [46], for MAR, partially reprocessed by Porquet et al. [44]).
Geosciences 12 00240 g002
Figure 3. Location of the studied granitic sites in the Central Pyrenees. Geological map modified from Clariana et al. [109]).
Figure 3. Location of the studied granitic sites in the Central Pyrenees. Geological map modified from Clariana et al. [109]).
Geosciences 12 00240 g003
Figure 4. Types of petrophysical samples. Type 1 hand samples (A), type 2 mini blocks (B) and types 3 standard paleomagnetic-AMS samples (C) and subsampling for magnetometers (type 4). (D) Caliper estimate of diameter of type-3 samples.
Figure 4. Types of petrophysical samples. Type 1 hand samples (A), type 2 mini blocks (B) and types 3 standard paleomagnetic-AMS samples (C) and subsampling for magnetometers (type 4). (D) Caliper estimate of diameter of type-3 samples.
Geosciences 12 00240 g004
Figure 5. (A) HF vs. LF susceptibility data obtained in minicores in the MPMS instrument for selected samples of the MLA granite (this work), the Marimanha pluton (Antolín et al. [46]) and compared with samples from other Pyrenean granites (data from Bielsa and Millares granites by Román-Berdiel et al. [84,85]). (B) Induced magnetic moment (emu) as a function of the external magnetic field (T). All samples from Mont-Louis Andorra granite (MLA).
Figure 5. (A) HF vs. LF susceptibility data obtained in minicores in the MPMS instrument for selected samples of the MLA granite (this work), the Marimanha pluton (Antolín et al. [46]) and compared with samples from other Pyrenean granites (data from Bielsa and Millares granites by Román-Berdiel et al. [84,85]). (B) Induced magnetic moment (emu) as a function of the external magnetic field (T). All samples from Mont-Louis Andorra granite (MLA).
Geosciences 12 00240 g005
Figure 6. Magnetic susceptibility (A) and density histograms (B) of the studied samples (only site mean values were used). The three plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis-Andorra (green). Ordinates (#) represent the number of data.
Figure 6. Magnetic susceptibility (A) and density histograms (B) of the studied samples (only site mean values were used). The three plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis-Andorra (green). Ordinates (#) represent the number of data.
Geosciences 12 00240 g006
Figure 7. Density vs. magnetic susceptibility in Pyrenean granites (only data from this study are shown, see data summary in Table 3). The three studied plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis- Andorra (green).
Figure 7. Density vs. magnetic susceptibility in Pyrenean granites (only data from this study are shown, see data summary in Table 3). The three studied plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis- Andorra (green).
Geosciences 12 00240 g007
Figure 8. Natural variability of density and magnetic susceptibility as a function of sample type and granitic body (only data from this study are shown, see data overview in Table 3). The three plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis- Andorra (green). The red box-plots comprise all data together.
Figure 8. Natural variability of density and magnetic susceptibility as a function of sample type and granitic body (only data from this study are shown, see data overview in Table 3). The three plutons are plotted with different colors Maladeta (orange), Marimanha (blue), Mont-Louis- Andorra (green). The red box-plots comprise all data together.
Geosciences 12 00240 g008
Figure 9. Density versus magnetic susceptibility graphs and regressions. (AC) pluton by pluton (sample types are differentiated by different colors). (D) All data together, outliers (in red) are excluded for the calculation. (E) All regressions together plus the running window averages (F).
Figure 9. Density versus magnetic susceptibility graphs and regressions. (AC) pluton by pluton (sample types are differentiated by different colors). (D) All data together, outliers (in red) are excluded for the calculation. (E) All regressions together plus the running window averages (F).
Geosciences 12 00240 g009
Figure 10. Density histograms for the three studied granitic massifs derived from the conversion of bulk susceptibility (from previous AMS studies) (counting interval 50 kg/m3).
Figure 10. Density histograms for the three studied granitic massifs derived from the conversion of bulk susceptibility (from previous AMS studies) (counting interval 50 kg/m3).
Geosciences 12 00240 g010
Figure 11. Magnetic susceptibility (κ) and density (ρ) mapping in the central Pyrenees. From west to east: Maladeta, Marimanha and Mont Louis-Andorra plutons.
Figure 11. Magnetic susceptibility (κ) and density (ρ) mapping in the central Pyrenees. From west to east: Maladeta, Marimanha and Mont Louis-Andorra plutons.
Geosciences 12 00240 g011
Figure 12. Geochemical data (% in wt of SiO2, Al2O3, Fe2O3 and Ti02) against rock density (ρ) and magnetic susceptibility (κ) as well as the linear regressions stablished among them in some type-2 samples frm Maladeta and Mont Louis Andorra plutons.
Figure 12. Geochemical data (% in wt of SiO2, Al2O3, Fe2O3 and Ti02) against rock density (ρ) and magnetic susceptibility (κ) as well as the linear regressions stablished among them in some type-2 samples frm Maladeta and Mont Louis Andorra plutons.
Geosciences 12 00240 g012
Figure 13. (A) Histogram of the magnetic susceptibility (site means) and density data in the Pyrenees (derived from the AMS database). (B) mean magnetic susceptibility vs. granite body size (km2).
Figure 13. (A) Histogram of the magnetic susceptibility (site means) and density data in the Pyrenees (derived from the AMS database). (B) mean magnetic susceptibility vs. granite body size (km2).
Geosciences 12 00240 g013
Figure 14. Individual box-plots of density and magnetic susceptibility (raw data) after the conversion of susceptibility into density in the Pyrenean granites (from the compilation by Porquet et al. [44]. The granites are sorted by increasing values of the density and susceptibility means. Ansignan (ANS) Aston (AST), Aya (AYA), Bassies (BAS), Bielsa (BIE), Borderes (BOR), Cauterets (CAU), Eaux Chaudes (ECH), Erce (ERC), Foix (FOI), Lacourt (LAC), Lys (LYS), Maladeta (MAL), Marimanha (MAR), Millares (MIL), Mont Luis-Andorra (MLA), Neouville (NEV), Panticosa (PAN), Posets (POS), Querigut (QUE), StArnac (SAR), Trois (TSG). All plutons together (ALL in red).
Figure 14. Individual box-plots of density and magnetic susceptibility (raw data) after the conversion of susceptibility into density in the Pyrenean granites (from the compilation by Porquet et al. [44]. The granites are sorted by increasing values of the density and susceptibility means. Ansignan (ANS) Aston (AST), Aya (AYA), Bassies (BAS), Bielsa (BIE), Borderes (BOR), Cauterets (CAU), Eaux Chaudes (ECH), Erce (ERC), Foix (FOI), Lacourt (LAC), Lys (LYS), Maladeta (MAL), Marimanha (MAR), Millares (MIL), Mont Luis-Andorra (MLA), Neouville (NEV), Panticosa (PAN), Posets (POS), Querigut (QUE), StArnac (SAR), Trois (TSG). All plutons together (ALL in red).
Geosciences 12 00240 g014
Table 1. Main statistical variables from previous studies in Maladeta (MAL), Marimanha (MAR) and the Mont-Louis Andorra (MAL) granites (based on data by Leblanc et al., 1994 [45]; Antolín-Tomás et al. [46] and Gleizes et al. [34] respectively). O.S: Outcropping Surface, S. km2; sites per sq. km/. All susceptibility data are expressed as µ S.I. Note that these colors are kept all along the paper to refer those granites.
Table 1. Main statistical variables from previous studies in Maladeta (MAL), Marimanha (MAR) and the Mont-Louis Andorra (MAL) granites (based on data by Leblanc et al., 1994 [45]; Antolín-Tomás et al. [46] and Gleizes et al. [34] respectively). O.S: Outcropping Surface, S. km2; sites per sq. km/. All susceptibility data are expressed as µ S.I. Note that these colors are kept all along the paper to refer those granites.
SitesO.S. (km2)S./km2minMaxMeanMedianRMSStan DevStand ErrorSkewnessKurtosis
MaladetaMAL2534150.617676213.7196.0227.478.04.92.0828.753
MarimanhaMAR62321.9430476203.2209.0225.598.512.50.2550.032
Mont Louis-AndorraMLA2545500.4614366197.6188.5207.964.54.00.2070.196
AllAll5699970.577676205.4194.0218.775.23.21.1905.222
Table 2. Summary of new density and susceptibility measurements (depending upon the sample type) carried out in this paper, specifying laboratory procedures, instrumentation, mean values and standard deviation per granite pluton (see body text for further explanations). Groups/Labs: MAD: Madrid, BCN: Barcelona, ZAZ: Zaragoza, TLSE: Toulouse. Density measurements in kg/m3 and susceptibility in µ S.I.
Table 2. Summary of new density and susceptibility measurements (depending upon the sample type) carried out in this paper, specifying laboratory procedures, instrumentation, mean values and standard deviation per granite pluton (see body text for further explanations). Groups/Labs: MAD: Madrid, BCN: Barcelona, ZAZ: Zaragoza, TLSE: Toulouse. Density measurements in kg/m3 and susceptibility in µ S.I.
Granite BodyGroupsSample TypeTypeWeight#Sites#MethodMeanStaDevLabs
MaladetaZAZ+MADHand blocks13–5 k2129Archimedes+UNE268563MAD
MaladetaBCN+MADMini blocks2<20 g2140Archimedes265848BCN+ MAD
Mont Louis-AndorraZAZ+MADHand blocks13–5 k2366Archimedes+UNE267047MAD
Mont Louis-AndorraBCN+MADMini blocks2<20 g2743Archimedes266152BCN+ MAD
Mont Louis-AndorraZAZ+TLSEPmag cores3&4≈25 g2121Geometric265754ZAZ
MarimanhaZAZPmag cores3≈25 g15111Archimedes w&w/o paraffine269778ZAZ
Total 128310 267112
Granite body#KLY2/KLY3BartingtonSM20KT20MPMSMVSMMethodMeanStaDevlabs
Maladeta1134 381096 Field and Lab171.9131.9MAD
Maladeta1269036 Labs187.0155.9BCN+ MAD
Mont Louis-Andorra1059 40591428 Field and Lab152.775.2MAD
Mont Louis-Andorra15612036 Labs233.1176.0BCN+ MAD
Mont Louis-Andorra4821 27Labs192.594.4ZAZ+ TLSE
Marimanha106106 Lab232.5133.9ZAZ
Table 3. Synthesis of statistics of petrophysical data (density and susceptibility) obtained in this paper. Dens(ity) and Magnetic Sus(ceptibility) separately for each pluton, as well as the global means.
Table 3. Synthesis of statistics of petrophysical data (density and susceptibility) obtained in this paper. Dens(ity) and Magnetic Sus(ceptibility) separately for each pluton, as well as the global means.
VariableMinMaxPointsMeanMedianRMSStd DevStd ErrorSkewnessKurtosis
Dens-MAL24372821482658266026597811−0.4450.594
Dens-MLA25202797652660265926615264−0.2640.849
Dens-MAR260128651526972691269878200.8880.115
Density All243728651282664266226646659−0.1060.149
Sus-MAL−13.2762.545204.5169.1257.6158.323.619.78141.531
Sus-MLA−5.4688.563176.3154207.8110.81417.05255.973
Sus-MAR31.9454.315232.5234.4266.0133.934.60.126−12.126
Sus-All−13.2762.5123193.5166.8234.6133.11217.9894.64
Table 4. Summary of density and magnetic susceptibility regressions (ρ = a κ + b) found in our dataset as a function of the pluton and the sample type (ρ in kg/m3 and x = κ in 10−6 S.I., n; number of sites considered).
Table 4. Summary of density and magnetic susceptibility regressions (ρ = a κ + b) found in our dataset as a function of the pluton and the sample type (ρ in kg/m3 and x = κ in 10−6 S.I., n; number of sites considered).
GraniteTypenabR
MLA1230.52125900.83
MLA2150.51825650.90
MLA3210.46525680.81
MAL1190.48926160.54
MAL2180.52925480.80
MAR3140.57325530.91
MAL1+2370.37026090.55
MLA1+2+3590.46325820.79
All1+31100.45825860.74
Type 11420.50426030.69
Type 22330.49325650.86
Type 33350.53325580.87
Run Wind1+2+31100.54125660.97
Table 5. Density main statistical variables for the three studied bodies after the conversion of susceptibility data; Maladeta (MAL), Marimanha (MAR) and Mont-Louis Andorra (MAL). Density expressed as kg/m3.
Table 5. Density main statistical variables for the three studied bodies after the conversion of susceptibility data; Maladeta (MAL), Marimanha (MAR) and Mont-Louis Andorra (MAL). Density expressed as kg/m3.
MinMax# SitesMean ρ (kg/m3)MedianRMSStd. Dev.% Std DevStd. ErrorSkewnessKurtosis
MAL25702932253268226722682421.6%32.0848.784
MAR2582282362267626802677532.0%70.2520.025
MLA25742764254267326682673351.3%20.2090.200
Table 6. Overview of deduced density (this work) and magnetic susceptibility statistical parameters (compilation by Porquet et al. [44] found in the 22 granitic bodies from the Pyrenees. Aston (AST), Aya (AYA), Bassies (BAS), Bielsa (BIE), Borderes (BOR), Cauterets (CAU), Eaux Chaudes (ECH), Erce (ERC), Foix (FOI), Lacourt (LAC), Lys (LYS), Maladeta (MAL), Marimanha (MAR), Millares (MIL), Mont Luis-Andorra (MLA), Neouville (NEV), Panticosa (PAN), Posets (POS), Querigut (QUE), StArnac (SAR), Trois (TSG).
Table 6. Overview of deduced density (this work) and magnetic susceptibility statistical parameters (compilation by Porquet et al. [44] found in the 22 granitic bodies from the Pyrenees. Aston (AST), Aya (AYA), Bassies (BAS), Bielsa (BIE), Borderes (BOR), Cauterets (CAU), Eaux Chaudes (ECH), Erce (ERC), Foix (FOI), Lacourt (LAC), Lys (LYS), Maladeta (MAL), Marimanha (MAR), Millares (MIL), Mont Luis-Andorra (MLA), Neouville (NEV), Panticosa (PAN), Posets (POS), Querigut (QUE), StArnac (SAR), Trois (TSG).
PlutonMinMax# SitesMeanMedianRMSStand DevStand ErrorSkewnessKurtosis
ANS6792222305.7230.5376.1224.347.81.4251.095
AST42221247183.2136.0291.8227.614.55.95045.635
AYA1263393160.785.0226.0159.816.61.5140.967
BAS3450888172.0167.5180.454.55.83.13117.679
BIE1782260180.7147.5226.2137.217.72.5137.584
BOR9364861326.4280.0362.4158.720.30.436−1.056
CAU18592199187.7179.0219.7114.48.10.9640.672
ECH1932028157.5151.5171.067.712.80.4140.262
ERC9435546169.2164.5176.751.37.61.6454.072
FOI2051669172.0197.0198.098.811.90.4610.747
LAC15946221274.0264.0281.264.714.11.0861.992
LYS3257599211.1187.0239.1112.911.30.8370.384
MAL7676253213.7196.0227.478.04.92.0828.753
MAR3047662203.2209.0225.598.512.50.2550.032
MIL7068954285.0253.5318.6143.719.61.0140.673
MLA14366254197.7188.5207.964.54.00.2070.196
NEV0541132212.9182.5226.678.06.80.9312.411
PAN63609111255.0246.0278.9113.410.80.340−0.183
POS17941269270.9270.0274.041.65.00.6541.009
QUE9799121227.6199.0263.4133.112.11.1301.744
SAR32465696303.1203.0671.4602.261.56.13838.426
TSG18132934237.7232.5240.436.56.30.7940.139
ALL046562219213.1188.0278.1178.73.812.282258.47
ANS2602306522273126912734121261.4261.101
AST2568376724726652640266812385.95145.650
AYA25732908932653261226548691.5150.969
BAS25842841882659265726593033.12517.710
BIE257530116026642646266574102.5167.614
BOR261629166127432717274486110.436−1.058
CAU257628861992668266326686240.9630.671
ECH25762739282651264826523770.4030.244
ERC26172758462658265526582841.6514.088
FOI25772845692659267326605360.4620.742
LAC26522816212714270927143581.0942.004
LYS25842877992680266726816160.8370.383
MAL257029322532682267226824232.0858.784
MAR25822823622676268026775370.2520.025
MIL260429395427202704272178111.0170.685
MLA257427642542673266826733520.2090.200
NEV256628591322681266526824240.9352.436
PAN260028951112704269927056160.341−0.183
POS26632789692713271227132230.6621.021
QUE257129981212689267426907271.1291.738
SAR2583508496273026762749326336.13838.426
TSG26642744342695269226952030.7910.142
ALL25665084222126802668268311221.729288.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pueyo, E.L.; Román-Berdiel, T.; Calvín, P.; Bouchez, J.L.; Beamud, E.; Ayala, C.; Loi, F.; Soto, R.; Clariana, P.; Margalef, A.; et al. Petrophysical Characterization of Non-Magnetic Granites; Density and Magnetic Susceptibility Relationships. Geosciences 2022, 12, 240. https://doi.org/10.3390/geosciences12060240

AMA Style

Pueyo EL, Román-Berdiel T, Calvín P, Bouchez JL, Beamud E, Ayala C, Loi F, Soto R, Clariana P, Margalef A, et al. Petrophysical Characterization of Non-Magnetic Granites; Density and Magnetic Susceptibility Relationships. Geosciences. 2022; 12(6):240. https://doi.org/10.3390/geosciences12060240

Chicago/Turabian Style

Pueyo, Emilio L., Teresa Román-Berdiel, Pablo Calvín, Jean Luc Bouchez, Elisabet Beamud, Conxi Ayala, Francesca Loi, Ruth Soto, Pilar Clariana, Aina Margalef, and et al. 2022. "Petrophysical Characterization of Non-Magnetic Granites; Density and Magnetic Susceptibility Relationships" Geosciences 12, no. 6: 240. https://doi.org/10.3390/geosciences12060240

APA Style

Pueyo, E. L., Román-Berdiel, T., Calvín, P., Bouchez, J. L., Beamud, E., Ayala, C., Loi, F., Soto, R., Clariana, P., Margalef, A., Bach, N., Schamuells, N., Rubio, F. M., Gimeno, A., de Arévalo, E. F., Rey-Moral, C., García, A., Martí, J., Casas, A. M., & García-Lobón, J. L. (2022). Petrophysical Characterization of Non-Magnetic Granites; Density and Magnetic Susceptibility Relationships. Geosciences, 12(6), 240. https://doi.org/10.3390/geosciences12060240

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