Abridged version
1 Introduction
The ‘Seuil du Poitou’ area (Fig. 1) was selected by the French National radioactive waste management agency (ANDRA) as a potential site for building an underground laboratory in granites. Seventeen cored boreholes, completed by petrographical, geochemical, geophysical and structural studies [9,26], allowed us to characterise the Charroux–Civray plutonic complex, concealed under Jurassic sediments. Three main types of granitoids were distinguished: medium and high-K calk-alkaline, and peraluminous ones. The aim of this work is to constrain the timing of peraluminous magmatism using muscovite 40Ar/39Ar dating. The closure temperatures of the 40Ar/39Ar isotopic systems have been assumed to be 300±30 °C for biotites, 325±30 for muscovites and 500±50 °C for amphiboles [16,25].
2 The granitoids of the Charroux–Civray plutonic complex
Geophysical, geological, and geochemical studies have shown that the Charroux–Civray complex is composite and dominated by deeply rooted medium- to high-K calk-alkaline granitoids. Geophysical data indicated the occurrence of a low density and magnetic susceptibility body, below these granitoids. This body was interpreted as corresponding to a large leucogranite pluton (about 30 km long, 10 to 15 km large and 5 to 10 km thick) [26].
The medium- and high-K calk-alkaline granitoids belong to the ‘Limousin Tonalitic line’ [18]. Their emplacement ages range from 356±5 Ma to 349±6 Ma (U/Pb dating on zircon [2]), at an estimated depth of 14±2 km [13]. 40Ar/39Ar dating on amphibole and biotite of these granitoids [1] indicates a fast global cooling (>100 °C Ma−1) between 350 and 340 Ma, which is interpreted as corresponding to the thermal equilibration of the last intrusions (around 350 Ma) with surrounding metamorphic rocks.
Peraluminous magmatism has been encountered in two boreholes: (i) the CHA109 borehole, as a muscovite ± biotite, highly fractionated leucogranite, to which wolframite–cassiterite veinlets and greisens are associated; (ii) the CHA111 borehole, as a muscovite–garnet leucogranite, belonging to a set of decimetric to metric leucogranitic veins, aplites and pegmatites cross-cutting a monzogranite. However, the main leucogranitic body (hereafter referred in French as CLP, ‘corps leucogranitique principal’), indicated by geophysical studies, has not been reached by any of the boreholes.
3 Samples and results
40Ar/39Ar dating was performed on three samples: a muscovite from the CHA111 borehole (sample 111), and two muscovites from the CHA109 borehole, one from the main facies (sample 109), another one from a wolframite–cassiterite greisen (sample 109g).
The three muscovite single grains analysed gave plateau ages (Table 1): 311.6±0.3 and 309.3±0.4 Ma respectively for samples 109 and 109g, 338.8±0.4 Ma for sample 111.
Résultats des datations 40Ar/39Ar des muscovites des leucogranites et greisen du complexe de Charroux–Civray. Les résultats sont présentés avec des marges d'erreur à .
40Ar/39Ar data on muscovites of the leucogranites and the greisen of the Charroux–Civray complex. The results are presented with a error.
Échantillon | Profondeur | Âge plateau | % de 39Ar sur plateau | Âge isochrone | 40Ar/36Ar à l'origine | Âge intégré |
109 | 171 m | 50 | 313,0±2,9 | 289 | 314,9±1,1 | |
109g | 171 m | 100 | 308,9±2,9 | 289 | 308,8±0,6 | |
111 | 162 m | 93 | 338,2±1,6 | 279 | 330,2±0,2 |
4 Discussion and conclusion
Assuming a conductive regime, a leucogranitic intrusion, whose size ranges from some hundreds of metres to some kilometres, will reach thermal equilibrium with enclosing rocks respectively in few tens of kiloyears to few megayears. 40Ar/39Ar dating of Late Devonian granitoids [1] indicates that a regional isotherm of 300±30 °C was reached at ca 345–340 Ma. Therefore, taking into account the small size of both leucogranitic intrusions, the obtained 40Ar/39Ar ages should be close to the emplacement ages.
The obtained data highlight the occurrence of two episodes of peraluminous leucogranitic magmatic activity. The first one, at ca 339 Ma, corresponding to the age of the muscovite garnet vein, is synchronous with several leucogranitic intrusions known in the Limousin area (100 km to the southeast): a vein in the Tulle antiform (338±4 Ma by 40Ar/39Ar, [21]) and two granite intrusions belonging to the Millevaches plutonic complex (ca 335 Ma, [17]). Systematic 40Ar/39Ar dating of metamorphic rocks along more than 100 km in the northern Limousin [1] indicates that this Mid-Carboniferous leucogranite magmatism could be synchronous with a regional exhumation stage.
A second one, at ca 312 Ma, corresponds to the age of the muscovite ± biotite leucogranite. Interestingly, this age is similar to the age of all Rare Metal Granite intrusions dated in the northern French Massif Central [1,6,10]. The onset of this specialised magmatism has been attributed to the thermal anomaly created by the delamination of the lower lithosphere, occurring by the end of the Variscan orogeny [7,19].
The interpretation of geophysical data [26] indicating the existence of the ‘CLP’ may be re-visited using the new 40Ar/39Ar dating on the leucogranites and those obtained recently [1] on the Late Devonian granitoids (amphibole and biotite ages ranging from 348 to 343 Ma, Fig. 1). Several scenarios may be suggested.
(1) The ‘CLP’ emplacement occurred between 350 and 340 Ma. In this case, 40Ar/39Ar dating obtained on Late Devonian granitoids would correspond to a nearly synchronous cooling of these granitoids and the ‘CLP’. Several peraluminous intrusions of similar age are known in the Millevaches massif (French Massif Central, [17]) and Vendée [8]. However, a leucogranitic intrusion with a thickness similar to the one proposed for the ‘CLP’ (5 to 10 km thick) is unknown in the French Massif Central during the Early Carboniferous. Thus, the CLP emplacement may have occurred at 340 Ma, at the same age as the muscovite–garnet leucogranite vein from borehole CHA111, which could represent a distal injection from the CLP.
(2) The emplacement of the ‘CLP’ occurred after 340 Ma, and thus may be synchronous with the main peraluminous magmatism episode in Limousin. It could have occurred at ca 312 Ma, in relation with the muscovite ± biotite leucogranite of borehole CHA109. In this case, a simplified numerical model of the ‘CLP’ cooling (Annex) shows that Late Devonian granitoids should have been reheated by the thermal aureole of the ‘CLP’ at ca 300–400 °C (Fig. 4). The lack of biotite 40Ar/39Ar age disturbance in the Late Devonian granitoids [1] may imply, either (i) that the amplitude and time duration of the thermal perturbation indicated by the model was not high enough to re-open significantly the 40Ar/39Ar isotopic system of biotites, or (ii) that the estimated thermal perturbation was indeed overestimated for several reasons, such as the depth location of the ‘CLP’, the span of time during which the leucogranitic magmas composing the ‘CLP’ have been emplaced. However, the new 40Ar/39Ar dating results do not allow us to give more precise information concerning the age and the geometry at the depth of the large leucogranitic body shown by geophysics.
1 Introduction
Le complexe plutonique de Charroux–Civray (Vienne) (Fig. 1) appartient au socle cristallin sous couverture mésozoı̈que du seuil du Poitou. Les 17 sondages carottés réalisés par l'Agence nationale de gestion de déchets radioactifs (Andra) ont fait l'objet de nombreuses études scientifiques, afin de reconnaı̂tre les caractéristiques de ce complexe. Trois types de granitoı̈des (calco-alcalin moyennement et fortement potassique, ou encore peralumineux) ont été identifiés [9]. Ce travail apporte de nouvelles données thermochronométriques sur les intrusions peralumineuses de ce complexe par datation 40Ar/39Ar sur muscovite.
En se basant sur les paramètres thermochronologiques de Dodson [11], des températures de fermeture isotopique de 300±30 °C pour les biotites, 325±30 °C pour les muscovites et 500±50 °C pour les amphiboles ont été utilisées [16,25].
2 Les complexe magmatique de Charroux–Civray
2.1 Les granitoı̈des fini-dévoniens
Les granitoı̈des calco-alcalins moyennement à fortement potassiques du complexe magmatique de Charroux–Civray appartiennent à la « ligne tonalitique limousine » [18]. La mise en place de ces granitoı̈des s'étale de 356±5 Ma à 349±6 Ma (UPb sur zircon [2]), à une profondeur qui a été estimée à 14±2 km (barométrie sur amphibole [13]). Les datations 40Ar/39Ar sur amphibole et biotite des granitoı̈des calco-alcalins [1] suggèrent un refroidissement global rapide (>100 °C Ma−1) entre 350 et 345 Ma. Ceci pourrait traduire la remise à l'équilibre thermique des dernières intrusions avec leur encaissant, à la fin de l'activité magmatique vers 350 Ma.
2.2 Le magmatisme peralumineux
Des leucogranites peralumineux ont été recoupés : (i) dans le sondage CHA109 sous la forme d'un leucogranite à muscovite dominante et biotite accessoire très fractionné (riche en Sn, Rb, Ta, U, W), associé à des filons à wolframite et cassitérite à épontes greisenisées ; (ii) dans le sondage CHA111, sous forme d'un ensemble de filons décimétriques à métriques recoupant un monzogranite. Cet ensemble est constitué de leucogranite à muscovite et grenat, leucogranites à deux micas, aplites et pegmatites, qui ne présentent, en revanche, aucun enrichissement significatif en éléments incompatibles.
Par ailleurs, l'inversion couplée des données gravimétriques et magnétiques [26] suggère la présence d'un corps léger et peu magnétique, de 5 à 10 km d'épaisseur, 10 à 15 km de large et 30 km de long, sous le complexe de Charroux–Civray (Fig. 2). Il a été interprété comme correspondant à un massif leucogranitique peralumineux, dont le toit serait à une profondeur minimale de l'ordre du kilomètre [26], car il n'a pas été atteint par les forages les plus profonds. Ce corps leucogranitique principal sera appelé CLP. Il est à noter que l'épaisseur suggérée pour le CLP est très importante pour un leucogranite. Le corps de faible densité indiqué par les données gravimétriques pourrait aussi correspondre à des formations de type « Unité inférieure des gneiss », mais les caractéristiques très magnétiques de cette dernière (orthogneiss à magnétite, amphibolites) ne permettent pas de rendre compte de la faiblesse de l'anomalie magnétique couplée à l'anomalie gravimétrique négative. Seul un corps de type leucogranitique peut rendre compte des anomalies géophysiques observées. Il pourrait s'agir là de l'un des plus volumineux massifs leucogranitiques de l'Ouest du Massif central, qui serait toutefois de puissance équivalente au massif leucogranitique de Mortagne, en Vendée.
Les caractéristiques minéralogiques et géochimiques des filons à grenat–muscovite du sondage CHA111 ne permettent donc pas de les considérer comme des apophyses du leucogranite à deux micas du CHA109, malgré sa proximité géographique. Il pourrait, en revanche, représenter une émission distale du CLP [9]. Afin de préciser les relations spatio-temporelles de ces deux unités entre elles et avec le CLP, des datations 40Ar/39Ar ont été effectuées sur une muscovite d'un filon de leucogranite du sondage CHA111 (échantillon 111) et sur deux muscovites du CHA109, l'une provenant du faciès principal (échantillon 109) et l'autre d'un greisen (échantillon 109g).
3 Datation 40Ar/39Ar
3.1 Méthodologie
Après broyage et tamisage, les grains à analyser ont été prélevés dans la fraction 0,3–0,8 μm par sélection sous la loupe binoculaire. Les analyses isotopiques ont été réalisées sur grains individuels. La méthode analytique est décrite dans [1] et [22]. Les échantillons ont été irradiés à Mc Master (Hamilton, Canada, flux total reçu : 9×1018 n cm−2) et le moniteur de flux utilisé est l'amphibole Hb3gr (âge : 1072 Ma ; [20]).
Les rapports isotopiques mesurés sont corrigés du bruit de fond, des valeurs du blanc, des interférences dues à K, Ca et Cl, de la contamination atmosphérique et de la discrimination de masse, les erreurs sur ces différents paramètres étant propagées dans le calcul d'erreur sur l'âge de l'échantillon. Le calcul d'erreur inclut l'erreur sur le rapport 40Ar/39Ar du moniteur, mais non l'erreur sur son âge. Toutes les erreurs reportées sont à .
3.2 Résultats
L'échantillon 111 présente un spectre d'âge très plat (Fig. 3), définissant un âge plateau de 338,8±0,4 Ma (93 % de l'39Ar, Tableau 1). Les deux premières étapes sont légèrement plus jeunes, mais cette perte limitée d'argon radiogénique reste négligeable vis-à-vis de la qualité de l'âge plateau défini sur la majeure partie du spectre.
Pour l'échantillon 109, deux étapes présentent des marges d'erreur plus élevées, qui correspondent à des analyses sur la cage de Faraday. Ces deux étapes ont été exclues pour préserver l'homogénéité du calcul de l'âge plateau, et seules les étapes sur photomultiplicateur ont été utilisées. Ceci a permis de définir un âge plateau à 311,6±0,3 Ma (sur 54 % de l'39Ar, Tableau 1). L'échantillon 109g permet de définir un excellent âge plateau (100 % de l'39Ar mesuré) (Fig. 3, Tableau 1) à 309,3±0,4 Ma. Si l'on exclut les deux premières étapes (plateau sur 82 % des gaz), un âge à 309,5±0,3 Ma peut être obtenu. Ces deux âges étant très proches ; c'est l'âge présentant la plus grande marge d'erreur (309,3±0,4 Ma) qui a été choisi pour l'échantillon 109g.
4 Discussion
4.1 Interprétation géodynamique des âges 40Ar/39Ar
Une intrusion leucogranitique, dont la taille pourrait varier de celle des filons du sondage CHA111 à celle du corps leucogranitique principal (CLP), s'équilibre thermiquement avec l'encaissant entre quelques dizaines de ka à environ 5 Ma selon sa taille et le mode de dissipation de la chaleur [5,14]. Or, les âges 40Ar/39Ar des biotites du complexe de Charroux–Civray [1] suggèrent que la température régionale était inférieure à 300±30 °C après 340 Ma. De ce fait, l'âge 40Ar/39Ar des micas des leucogranites devrait donc être proche, à quelques millions d'années près, au plus, de l'âge de leur mise en place.
L'âge de mise en place du filon leucogranitique du sondage CHA111, à 339 Ma, est à rapprocher de celui d'intrusions leucogranitiques connues 100 km au sud-est. Dans l'antiforme de Tulle [21], un âge 40Ar/39Ar sur muscovite à 337±4 Ma a été obtenu sur un filon leucogranitique, interprété comme étant proche de l'âge de mise en place. Dans le complexe granitique du Millevaches, les leucogranites de Saint Julien (15×7 km) et de Goulles (20×8 km) ont été datés respectivement à 332±6,5 Ma et 336±7 Ma, par la méthode Rb/Sr sur roche totale [17]. De plus, les âges 40Ar/39Ar des micas des métamorphites du Limousin sont constants à ca 340 Ma sur plus de 100 km et interprétés comme traduisant une exhumation du socle au Carbonifère moyen [1]. Ainsi, un épisode magmatique leucogranitique, datant du Carbonifère moyen à l'échelle du Limousin, synchrone d'une phase d'exhumation, est mis en évidence.
Les deux muscovites du sondage 109 (leucogranite et greisen) ont fourni des âges voisins (Tableau 1), de l'ordre de 312 et 309 Ma. Comme pour l'échantillon 111, l'âge de la muscovite du leucogranite peut être interprété comme étant très proche de l'âge de mise en place. L'âge du greisen pourrait indiquer des circulations hydrothermales (1–2 kbar, 300–400 °C [4,12]) légèrement plus tardives, associées aux intrusions leucogranitiques lors d'un second épisode d'exhumation du socle, ayant conduit à des minéralisations à Sn et W [10].
Il est intéressant de noter que cet âge à environ 310 Ma est commun aux intrusions leucogranitiques à métaux rares du Nord du Massif central français (Chédeville et Beauvoir [6], Richemont [1], Montebras [10]) et au gisement de tungstène de Vaulry [10]. En outre, l'étude des inclusions fluides [4] postérieures à la formation du greisen indique une décompression, liée à la remontée rapide du bâti hercynien, similaire à celle proposée au Stéphanien dans le Limousin [12,24]. Le complexe de Charroux–Civray semble donc avoir eu une évolution semblable à celle du Nord-Ouest du Massif central entre 310 et 300 Ma, avec la mise en place de magmas leucogranitiques spécialisés, lors du passage à un régime d'extension généralisé, en relation avec une possible délamination de la lithosphère inférieure [7,19].
En résumé, dans le complexe de Charroux–Civray comme dans le Limousin, la période 340–310 Ma semble marquer les limites d'un ensemble d'épisodes leucogranitiques, de plus en plus fractionnés lorsque l'on se rapproche de 310 Ma. Cette évolution magmatique pourrait traduire les changements tectoniques survenus dans l'Ouest du Massif central, lors d'une phase d'extension synconvergence [3], transition entre la fin de la l'épaississement crustal vers 350–360 Ma et l'effondrement généralisé de la chaı̂ne vers 300 Ma [3,15].
4.2 Relations entre les intrusions leucogranitiques et le corps leucogranitique principal (CLP)
L'interprétation des données géophysiques [26] suggérant l'existence d'un CLP sous les granitoı̈des dévoniens peut être rediscutée à la lumière des âges 40Ar/39Ar sur les leucogranites et sur les granitoı̈des du complexe de Charroux–Civray (donnant des âges de fermeture sur biotite de 348 à 343 Ma [1], Fig. 1). Plusieurs scénarios peuvent être envisagés.
(1) La mise en place du CLP se situe entre 350 et 340 Ma. Dans ce cas, les âges 40Ar/39Ar obtenus sur les granitoı̈des dévoniens [1] traduiraient un refroidissement de ces derniers de manière sub-synchrone à celui du CLP. Des intrusions peralumineuses de taille pluri-kilométrique et d'âge voisin sont connues à 100 km au sud-est, dans le complexe granitique du Millevaches (cf. ci-dessus) et en Vendée (âge UPb sur zircon à 340±4 Ma pour le granite de Largeasse, de type Guéret [8]). De ce fait, la mise en place du CLP pourrait avoir eu lieu vers 340 Ma, âge du filon de leucogranite à muscovite–grenat du sondage CHA111, qui en serait ainsi un émissaire. Cependant, une intrusion leucogranitique d'une épaisseur semblable à celle proposée pour le CLP (5 à 10 km) n'a jamais été signalée dans le Massif central au Carbonifère Inférieur.
(2) La mise en place du CLP est postérieure à 340 Ma, par exemple vers 312 Ma, synchrone du leucogranite à muscovite ± biotite du sondage CHA109, en conformité avec l'âge de la fin de l'épisode majeur de magmatisme leucogranitique en Limousin. Dans ce cas, un modèle simplifié de l'influence thermique du CLP (annexe) montre que des points situés dans les granitoı̈des calco-alcalins dévoniens à environ 1000 à 2000 m du contact avec le CLP, initialement à une température de 200 °C environ (Fig. 4), devraient avoir été réchauffés à des températures allant de 300 à 400 °C (Fig. 4). Ces températures sont conformes à celles dérivant de l'étude des paragenèses d'altération et des inclusions fluides dans ces granitoı̈des, qui indiquent une phase de circulations fluides microfissurales, à des températures minimales de 130 à 350 °C, attribuées à la mise en place de magmas leucogranitiques [4]. Or, les systèmes isotopiques 40Ar/39Ar des minéraux des granitoı̈des calco-alcalins sus-jacents au CLP (notamment la biotite) ne semblent pas avoir été perturbés [1]. Ceci pourrait signifier, soit (i) que l'amplitude et la durée de la perturbation thermique indiquée par le modèle (température maximale de 300 à 400 °C, pour une durée inférieure à 1,5 Ma, quelle que soit la distance retenue, Fig. 4), n'ont pas été suffisantes pour ré-ouvrir de façon significative le système isotopique 40Ar/39Ar des biotites, soit (ii) que la perturbation thermique a été plus faible que celle indiquée par le modèle (par exemple, à cause d'une mise en place du CLP sous forme d'intrusions multiples étagées dans le temps, ou bien à une profondeur légèrement plus grande, hypothèses qui restent compatibles avec les incertitudes des données géophysiques). Ainsi, les données disponibles ne permettent pas de trancher entre ces différents scénarios.
5 Conclusion
Les données obtenues dans cette étude permettent de mettre en évidence deux épisodes magmatiques peralumineux de type leucogranitique dans le complexe de Charroux–Civray. Un premier épisode s'est exprimé à 340 Ma, sous la forme d'un filon de leucogranite à muscovite–grenat (sondage CHA111). Ce magmatisme, que l'on retrouve à la même époque dans le Limousin, pourrait être synchrone d'une première phase d'exhumation du socle varisque. Un second épisode (sondage CHA109) s'est exprimé vers 310 Ma sous la forme d'un leucogranite à muscovite ± biotite spécialisé, associé à des minéralisations W ± Sn. Cet âge est caractéristique de toutes les intrusions de granite à métaux rares du Nord du Massif central et semble marquer le passage à un régime d'extension généralisé. En revanche, l'ensemble des données disponibles ne permet pas de préciser l'âge et la géométrie du complexe leucogranitique profond intrusif dans les granitoı̈des dévoniens, mis en évidence par la géophysique.
Remerciements
Les auteurs remercient les deux rapporteurs (M. Pichavant et P. Monié) pour leurs remarques constructives. Ce travail a bénéficié d'un financement de l'Andra. Les mesures 40Ar/39Ar ont été effectuées au laboratoire du CNRS 6526 Geosciences Azur de l'université de Nice–Sophia Antipolis. Contribution CRPG n∘ 1613.
Appendix A
Pour étudier l'influence de la mise en place du CLP sur son encaissant, un modèle au premier ordre des transferts de chaleur en régime conductif transitoire a été réalisé par une méthode d'éléments finis 2D. Le calcul est effectué sur une grille non uniforme, de dimensions globales 40 km × 20 km de profondeur, représentant une coupe simplifiée du complexe de Charroux–Civray selon un axe SW–NE. Il est supposé que le CLP puisse être représenté en moyenne par un corps de 10 km de largeur sur 8 km d'épaisseur (Fig. 4, [26]), se mettant en place en une seule fois. De ce fait, ce modèle ne pourra donner qu'une estimation de la durée maximale de la perturbation thermique induite par la mise en place du CLP.
Les paramètres physiques des ensembles lithologiques simplifiés (encaissant indifférencié, d'une part, et CLP, d'autre part) sont donnés dans le Tableau 2, d'après [14].
Paramètres physiques du modèle 2D des transferts de chaleur (voir annexe).
Physical parameters used for the 2D model of heat transfer (see appendix).
Lithologie | Conductivité thermique | Production de chaleur | Chaleur spécifique |
(W m−1 K−1) | (μW m−3) | (J m−3 K−1) | |
Encaissant | 2,5 | 1 | 3 × 106 |
Leucogranites | 3 | 3,5 | 3,5 × 106 |
Les conditions aux limites sont les suivantes : en surface, une température imposée à 10 °C ; sur les frontières latérales, conditions adiabatiques (flux latéral nul) ; à la base, à 20 km, un flux de 45 mW m−2 a été imposé. Le champ de température initial utilise un gradient moyen au niveau de l'intrusion de 30 °C km−1, en accord avec les données thermométriques [13]. La nature supposée du CLP suggère des températures du magma de l'ordre de 700 °C [23]. Celle-ci a été augmentée de 5 à 10 %, pour simuler l'effet de la chaleur latente de cristallisation [14], ce qui conduit à une température du magma de l'ordre de 750 °C.
Pour évaluer l'influence des incertitudes géométriques (position des points de données par rapport au CLP) sur l'évolution du champ de température, quatre points repères dans l'encaissant ont été choisis (Fig. 4b). Les deux premiers (points A et B) sont à 1000 m au-dessus du CLP, respectivement sur son axe central et à sa périphérie. Les deux derniers points (C et D) ont le même positionnement horizontal, mais sont situés 2000 m au-dessus du CLP. Les profondeurs de ces quatre points sont respectivement de 5,5 km et 6,5 km, correspondant à une température avant intrusion de 200±15 °C, en accord avec les données thermiques [4,13].