Le couplage entre la tectonique des plaques et la convection mantellique est généralement étudié avec des modèles numériques dans lesquels le nombre de plaques est fixé a priori. Nous présentons ici un nouvel outil : le modèle MACMA (MAnteau Convectif Multi-Agents) simule une tectonique variable dans le temps, avec des frontières de plaques mobiles dans une géométrie cylindrique 2D. Les vitesses de plaques sont calculées à partir d'un bilan de forces individuel et des mécanismes explicites sont utilisés pour simuler la dynamique des frontières de plaques, décrivant en particulier l'initiation des subductions, la migration des fosses de subduction, la rupture de lithosphère continentale (ouverture d'un océan) et la suture de plaques. La géométrie et l'état thermique du système sont ainsi mis à jour à chaque pas de temps. Dans notre modèle, le nombre de plaques tectoniques émerge donc naturellement de leur interaction. Cette approche est basée sur des systèmes multi-agents en interaction thermique et mécanique. Notre modèle met en jeu quatre types d'agents : cellules de convection, plaques lithosphériques, continents et frontières de plaques. Ces agents collectent des informations dans leur environnement pour prendre des décisions en fonction de règles de comportement incorporées dans le modèle géophysique. Le développement de cette méthode présente deux objectifs : (1) examiner dans quelle mesure des lois de comportement analytiques et empiriques peuvent affecter la dynamique des plaques, et (2) évaluer l'impact des mécanismes de surface de la tectonique des plaques sur l'évolution thermique du manteau. Nous obtenons des prédictions des forces motrices et des vitesses de plaques de la tectonique qui sont en bon accord avec les observations. De plus, notre modèle met en évidence le rôle clé joué par la distribution des âges du plancher océanique sur les fluctuations de flux de chaleur à la surface du manteau, soulignant ainsi l'importance des changements structurels à l'échelle locale (par exemple, la création d'une nouvelle dorsale océanique). D'autre part, en utilisant dans nos simulations les estimations courantes de la température et de la viscosité du manteau, la configuration actuelle des plaques correspond à une phase de décroissance du flux de chaleur, corroborant ainsi le comportement thermique déduit des reconstructions tectoniques récentes. Sur le long terme, l'histoire thermique terrestre est généralement étudiée à partir de lois d'échelle reliant le flux de chaleur surfacique à la température et à la viscosité du manteau convectif. Cependant, une telle paramétrisation du flux conduit à une divergence de la température dans le passé, quand elle est utilisée dans le bilan thermique global de la planète pour calculer l'état thermique de la Terre en remontant le temps à partir de la valeur actuelle de flux de surface. Afin d'étudier le taux de refroidissement de la Terre en utilisant des processus tectoniques explicites, nous avons mis à profit le faible temps de calcul du modèle MACMA pour étudier l'effet d'un large panel de paramètres physiques sur l'évolution thermique à long terme du système. Pour des paramètres terrestres, un taux de refroidissement moyen de 60 K par milliard d'années est obtenu sur 3 Ga, ce qui respecte les contraintes pétrologiques et rhéologiques dont disposent les géologues, sans pour autant invoquer une tectonique plus lente dans le passé, comme cela est souvent fait pour éviter la catastrophe thermique obtenue par les méthodes classiques de paramétrisation. Deux échelles de temps apparaissent donc dans l'évolution du flux de chaleur : une décroissance monotone sur le long terme (plusieurs milliards d'années) et des fluctuations de forte amplitude à court terme dues aux réarrangements structurels de la tectonique des plaques. Nous montrons que la viscosité du manteau n'est pas un paramètre clé de l'évolution thermique du système. Le taux de refroidissement de la planète dépend principalement de sa capacité à remplacer du plancher océanique vieux et isolant par une lithosphère jeune donc fine. Ainsi, les principaux paramètres de contrôle sont liés aux processus de surface, tels que le seuil de rupture continentale et l'âge critique d'initiation des subductions. Nous déduisons de cette étude que la convection mantellique seule ne peut rendre compte de la complexité du refroidissement du manteau à différentes échelles de temps, et que ce sont les processus tectoniques de surface qui contrôlent l'évolution thermique de la Terre.

Les régions boréales seront les plus affectées par le réchauffement climatique, c'est pourquoi cette thèse s'est intéressée à l'étude du cycle hydrologique de ces régions. Une nouvelle méthodologie d'extraction du volume de neige à partir de données radiométriques sur l'ensemble des régions boréales a été validée et a permis de montrer une différence de comportement entre 1988 et 2006 sur les variations du volume de neige de l'Eurasie et celui de l'Amérique du Nord. L'étude des variations de volume des eaux de surface des régions arctiques est par contre plus difficile à estimer avec les données satellitaires actuelles. C'est pourquoi un nouveau projet de satellite, la mission SWOT (Surface Water and Ocean Topography), a été proposé qui vise à fournir des cartes de hauteurs d'eau sur l'ensemble du globe. L'apport de cette mission à l'étude des régions arctiques a été estimé en implémentant une modélisation de l'Ob inférieur, fleuve de l'Ouest sibérien, en couplant un modèle hydrologique à grande échelle et un modèle hydrodynamique d'inondations. En ajustant certains paramètres de ces modèles, il a été possible d'obtenir une modélisation réaliste du débit et des hauteurs d'eau du fleuve. L'utilisation d'un lisseur de Kalman d'ensemble local a permis de montrer que les données SWOT devraient permettre de réduire significativement (de plus de 50%) les erreurs de modélisation. L'intérêt de la mission pour l'observation du débit de l'ensemble des fleuves a aussi été estimé en se basant sur un bilan d'erreur préliminaire. L'utilisation de courbes de tarage, ainsi que la prise en compte des erreurs de mesure SWOT ont permis de montrer que ces nouvelles données devraient permettre d'estimer un débit moyen avec une erreur inférieure à 30% pour tous les fleuves ayant une profondeur de plus de 1 m. D'autre part, il a été montré que l'erreur sur l'estimation du débit mensuel due seulement à l'échantillonnage temporel de SWOT diminue avec l'aire drainée et que, pour une aire drainée supérieure à 6 900 km2, cette erreur devrait être inférieure à 20%. Enfin, une méthodologie simple a permis de calculer que la variation annuelle totale du volume de l'ensemble des lacs est de l'ordre de 9 000 km3. Les données spatiales actuelles ne peuvent pas en observer plus de 15%. Selon notre estimation, SWOT quant à lui devrait pouvoir mesurer entre 50% et 65% de cette variation de volume.

Ce travail de thèse est consacré à l'étude de la structure électrique de la croûte terrestre dans deux régions du Vietnam, au Nord dans le delta du Fleuve Rouge et au Sud dans la région de la faille de la Rivière Saigon, et en Afrique de l'Ouest, où les données utilisées sont celles acquises durant l'Année Internationale de l'Electrojet Equatorial (AIEE). La technique utilisée est celle du sondage magnétotellurique. Deux des régions étudiées (le sud du Vietnam et l'Afrique de l'Ouest) sont situées au voisinage de l'équateur magnétique. Du fait de la présence de l'électrojet équatorial, l'hypothèse de l'onde plane, sur laquelle est basée la magnétotellurique, n'y est vérifiée que pour une gamme de périodes limitée. La modélisation de l'induction par un électrojet gaussien nous a permis de caractériser l'effet de source associé à l'électrojet équatorial, et de déterminer pour chacune des zones étudiées les gammes de périodes pour lesquelles cet effet de source est négligeable. Nous avons par ailleurs montré que l'effet de source associé à l'électrojet équatorial dépend des principales caractéristiques de la structure électrique du sous-sol, et de la distance entre le centre de l'électrojet équatorial et la position examinée. Dans le delta du fleuve Rouge, le profil magnétotellurique est orienté SO-NE. Long d'environ 32 km, il recoupe les failles majeures du système de failles du Fleuve Rouge. Les données magnétotelluriques ont été inversées pour obtenir une coupe transversale de résistivité électrique pour les trois premiers kilomètres. Cette coupe de résistivité électrique et une coupe de densité obtenue par modélisation des données gravimétriques acquises le long du même profil, ainsi que des informations déduites de forages effectués au voisinage du profil ont été utilisées pour proposer une coupe géologique de la région étudiée. Dans la région de la faille de la Rivière Saigon, le profil magnétotellurique, long d'environ 15 km, recoupe le trajet supposé de la faille. La structure électrique des premiers kilomètres de la croûte a été déterminée par inversion 2-D des données pour les gammes de périodes qui ne sont pas affectées par l'électrojet équatorial (< 0,7 s). La coupe transversale de résistivité électrique ainsi obtenue nous a permis de proposer une nouvelle localisation de la faille de la Rivière Saigon, à la limite entre un domaine résistant, le bloc de Da Lat au nord-est, et un domaine moins résistant, le bloc de Can Tho au sud-ouest. Cette nouvelle localisation est à environ 2,2 km au sud-ouest de la localisation généralement supposée pour cette faille. En Afrique de l'Ouest, nous avons montré que les données de jour sont affectées par l'électrojet équatorial pour toute la gamme de période étudiée (>120 s); ces données n'ont donc pas été utilisées pour déterminer la structure électrique de la croûte et du manteau supérieur. L'inversion 2-D des données de nuit, montre la présence d'une structure subverticale, relativement conductrice, qui sépare le Craton Ouest-africain très résistant des terrains situés sous le bassin du fleuve Sénégal.

Le suivi temporel est un processus d'acquisition et d'analyse d'acquisitions multiples répétées au même endroit sur la même cible à différentes périodes de temps. Cela s'applique bien à l'exploration sismique quand les propriétés de la cible varient au cours du temps comme pour les réservoirs pétroliers. Cette technique de sismique, dite 4D en raison de l'intégration du temps dans la construction des images, permet une détection et une estimation des variations du sous-sol survenues lors de l'évolution en temps du milieu. En particulier, dans l'industrie, le suivi et la surveillance peuvent améliorer notre compréhension d'un réservoir de pétrole/gaz ou d'un site de stockage de CO2. Analyser la sismique 4D peut aider à mieux gérer les programmes de production des réservoirs. Ainsi, des acquisitions répétées permettent de suivre l'évolutiondes fronts de fluide injectés: on peut optimiser les programmes d'injection de fluides pour une récupération améliorée des hydrocarbures (enhanced oil recovery). Plusieurs méthodes ont été développées pour l'imagerie variable dans le temps en utilisant les informations des ondes sismiques. Dans ma thèse, je montre que l'inversion de forme d'onde complété (FWI) peut être utilisée pour cette imagerie. Cette m'méthode offre des images sismiques quantitatives haute résolution. Elle est une technique prometteuse pour reconstruire les petites variations de propriétés physiques macro-échelle du sous-sol. Sur une cible identifiée pour ces imageries 4D, plusieurs informations a priori sont souvent disponibles et peuvent être utilisées pour augmenter la résolution de l'image. J'ai introduit ces informations grâce à la définition d'un modèle a priori dans une approche classique FWI en l'accompagnant de la construction d'un modèle d'incertitudes a priori. On peut réaliser deux reconstructions indépendantes et faire la différence les reconstruits: on parle de différence parallèle. On peut aussi effectuer une différence séquentielle o'u l'inversion de l'ensemble de données de la second acquisition, dite moniteur, se fait 'a partir du modèle de base et non plus à partir du modèle utilisé initialement. Enfin, l'approche double-différence conduit à l'inversion des différences entre les deux jeux de données que l'on rajoute aux données synthétiques du modèle de base reconstruit. J'étudie quelle stratégie est à adopter pour obtenir des changements vitesse plus précis et plus robustes. En plus, je propose une imagerie 4D ciblée en construisant un modèle d'incertitude a priori grâce 'a une information (si elle existe) sur la localisation potentielle des variations attendues. Il est démontré que l'inversion 4D ciblée empêche l'apparition d'artéfacts en dehors des zones cibles: on évite la contamination des zones extérieures qui pourrait compromettre la reconstruction des changements 4D réels. Une étude de sensibilité, concernant l'échantillonnage en fréquence pour cette imagerie 4D, montre qu'il est nécessaire de faire agir simultanément un grand nombre de fréquences au cours d'un cycle d'inversion. Ce faisant, l'inversion fournit un modèle de base plus précis que l'approche temporelle, ainsi qu'un modèle des variations 4D plus robuste avec moins d'artéfacts. Toutefois, la FWI effectuée dans le domaine temporel semble être une approche plus intéressante pour l'imagerie 4D. Enfin, l'approche d'inversion 4D régularisée avec un modèle a priori est appliquée sur des ensembles de données réelles d'acquisitions sismiques répétées fournis par TOTAL. Cette reconstruction des variations locales s'inscrit dans un projet d'injection de vapeur pour améliorer la récupération des hydro-carbures: Il est possible de reconstituer des variations de vitesse fines causées par la vapeur injectée.

La composition chimique de la croûte continentale a significativement évolué à la transition Archéen-Protérozoïque (3000-2500 Ma), témoignant de changements géodynamiques majeurs à cette époque. Afin d'étudier l'expression et les origines de ces changements, qui sont encore mal contraints, j'ai étudié une diversité de granitoïdes qui se sont mis en place dans cette gamme d'âges à la marge Nord du craton du Kaapvaal, en Afrique du Sud. Ce travail a permis de préciser la typologie et l'origine des granitoïdes tardi-archéens ; ceux-ci peuvent être classés dans trois grands groupes : (1) Les sanukitoïdes, représentés en Afrique du Sud par le pluton de Bulai, sont des magmas dérivant de l'interaction entre une péridotite mantellique et un composant riche en éléments incompatibles (TTG, liquide issu de la fusion de sédiments, et, plus rarement, fluide aqueux). Les sanukitoïdes peuvent être classés en deux groupes distincts, selon les mécanismes de cette hybridation : les low-Ti sanukitoids proviennent d'une simple hybridation du liquide silicaté avec la péridotite, alors que les high-Ti sanukitoids sont issus de la fusion d'un assemblage métasomatique à amphibole et phlogopite, résultant de ces interactions. Enfin, les mécanismes de différenciation des suites sanukitoïdes au niveau de la croûte sont contrôlées par des mécanismes de cristallisation fractionnée ou (moins vraisemblablement) de fusion partielle. (2) Les sanukitoïdes " marginaux ", représentés dans le craton du Kaapvaal par les plutons de Mashashane, Matlala, Matok et Moletsi, sont des granitoïdes résultant de l'interaction entre des sanukitoïdes et des magmas provenant de la fusion de croûte préexistante. Etant donné la large gamme de sources possibles (TTG, métasédiments, roches mafiques) d'un craton à l'autre, ce groupe est extrêmement diversifié. Leurs mécanismes de différenciation sont contrôlés par la cristallisation fractionnée. (3) Certains granites, tels que le batholite de Turfloop en Afrique du Sud, sont directement issus de la fusion de lithologies crustales (TTG, métasédiments et amphibolites). Au sein du craton du Kaapvaal, l'évolution spatio-temporelle du magmatisme tardi-archéen suit un schéma très caractéristique : les TTG se mettent en place entre ~3300 et ~2800 Ma, puis laissent la place à la genèse de l'ensemble des granitoïdes présentés ci-dessus, qui se déroule entre 2780 et 2590 Ma. Cette séquence d'évènements est reproduite au sein de tous les cratons du monde à la fin de l'Archéen. Elle témoigne de l'avènement des processus de recyclage crustal, puisque, par opposition aux TTG archéennes qui dérivent de métabasaltes juvéniles, les magmas tardi-archéens sont issus à la fois de la différenciation intracrustale et de l'interaction entre une péridotite et du matériel continental introduit dans le manteau. Cette dualité de processus pétrogénétiques est aussi très typique des épisodes magmatiques qui ont lieu à la fin des cycles de subduction-collision post-archéens. Ainsi, l'évolution de la composition des granitoïdes entre 3000 et 2500 Ma traduit vraisemblablement l'initiation d'une forme de tectonique des plaques proche du régime actuel. Celle-ci serait liée au refroidissement planétaire global, qui a probablement entraîné un " effet de seuil " dans l'évolution de l'épaisseur de la croûte océanique ainsi que la rhéologie et le volume de la croûte continentale, permettant ainsi à la subduction et à la collision de ne devenir thermo-mécaniquement stables qu'à partir de la fin de l'Archéen.

Compréhension et quantification du phénomène de l'électrofiltration dans les roches, signaux précurseurs de la rupture, applications aux zones volcaniques et aux zones de subduction - Conversions sismo-électromagnétiques.<br /><br />Mes thèmes de recherche concernent donc principalement les effets électriques liés aux circulations d'eau dans les roches. Les signaux électriques induits par les circulations d'eau sont appelés effets d'électrofiltration ou électrocinétiques et sont dus à la présence d'ions dans l'eau pouvant créer un courant électrique. Ce courant électrique engendre lui-même un champ magnétique. Ces effets reflètent principalement l'interaction électrique entre la roche et l'eau (appelée double- ou triple-couche électrique). En présence d'eau, des réactions chimiques et physiques se produisent à l'interface entre la roche et l'eau, dont il résulte un excédent d'ions négatifs liés à la roche, et par suite, un excédent d'ions positifs sur une couche de faible épaisseur au voisinage de l'interface, le reste du fluide restant neutre [voir Fig. B1 de Guichet et al., GJI, 2006]. Lorsque le fluide est en mouvement, ces ions positifs sont entraînés et engendrent ainsi un courant électrique, dit courant électrique de convection. Ce courant est contre-balancé par un courant de conduction (loi d'Ohm), et il en résulte une différence de potentiel électrique, on parle de potentiel d'électrofiltration ou potentiel électrocinétique.<br />La différence de potentiel électrique ainsi générée est proportionnelle au gradient de pression d'eau, ce qui en fait une méthode de choix pour la détection des circulations d'eau en subsurface. Les mesures en laboratoire sont effectuées de l'échelle centimétrique à l'échelle métrique et sont nécessaires pour contraindre les modèles de calcul des champs électriques et magnétiques induits par les circulations de fluides dans la croûte. Il est à noter que sur le terrain, les Potentiels Spontanés (PS), c'est-à-dire les potentiels électriques mesurés passivement entre deux électrodes, dont les mesures se sont largement développées ces dix dernières années, peuvent être induits par une circulation d'eau (soumise à une force liée à un gradient de pression), mais aussi par un flux de chaleur (soumis à une force liée à un gradient de température), ou encore par un gradient de concentration chimique. Ceci peut se voir dans les relations des flux couplés [Onsager, 1931 ; Sill, 1983] où chaque flux Ji (flux électrique, flux hydrique, flux de chaleur, flux de matière) est lié aux forces Xj (gradient de potentiel électrique, de pression hydrique, de température, de concentration) à travers les coefficients de couplage Lij. Pour ma part, j'ai concentré mes recherches sur une meilleure quantification du champ électrique lié aux circulations d'eau dans les roches. Le coefficient de couplage électrocinétique ou d'électrofiltration, , peut-être mesuré en laboratoire. Il est à noter que le coefficient d'électrofiltration ne dépend pas en principe de propriétés telles que porosité ou perméabilité, ce qui amène le gradient de pression hydrique à jouer le premier rôle dans la génération des potentiels d'électrofiltration [pour plus de détails on peut consulter le paragraphe " Electrokinetic Phenomena " de Jouniaux and Pozzi, JGR, 1997, et l'appendice B de Guichet et al., GJI, 2006].<br />Cette thématique de recherche a donc un intérêt potentiel, et déjà en partie éprouvé, pour la détection à distance et la caractérisation des fluides dans le sous-sol. La communauté scientifique essaie de développer son potentiel en tant qu'outil d'imagerie, d'une part dans le domaine des ressources en eau : profondeur de nappe aquifère, suivi de transferts hydrauliques, caractérisation des flux contaminés, et d'autre part dans le domaine de la géophysique des réservoirs : détection d'hydrocarbones grâce aux conversions sismo-électromagnétiques.

Le fer est nécessaire à la vie. Alors qu'il se trouve en abondance à la surface terrestre, il est extrêmement rare à la surface de l'océan (quelques nanogrammes par litre d'eau de mer), de sorte qu'il y limite la production primaire dans de nombreuses régions et donc le pompage du CO2 atmosphérique par l'océan. Le cycle océanique du fer est donc un maillon essentiel du système climatique. Alors que depuis une trentaine d'années les biogéochimistes marins ont concentré leurs efforts sur l'étude des concentrations de fer dans l'océan, il est devenu possible, depuis quelques années, d'étudier les isotopes du fer dans l'océan. Nous verrons en quoi ces isotopes nous apportent de nouvelles informations sur les sources de fer à l'océan et sur son cycle au sein de l'océan.

Le travail de recherche présenté dans ce mémoire de thèse porte sur l'étude de la conductivité électrique et de la permittivité diélectrique du sol dans le domaine des fréquences moyennes (100 kHz-10 MHz), l'objectif étant de développer une nouvelle méthode de cartographie simultanée de ces paramètres électriques. La conductivité électrique et la permittivité diélectrique étant sensibles aux paramètres d'état du sol, et notamment à sa teneur en eau, le couplage des informations apportées par ces deux paramètres doit permettre d'estimer la valeur de la teneur en eau en chaque point de mesure. Des mesures en laboratoire ont été réalisées afin de connaître l'évolution de ces deux paramètres en fonction de la fréquence. La permittivité diélectrique et la conductivité électrique d'échantillons de sable, de mélanges artificiels de sable et d'argiles (kaolinite et bentonite) en proportions variées et de sol limoneux ont été mesurées pour différentes teneurs en eau. Nous avons étudié l'effet de la teneur en argile, de la minéralogie des argiles impliquées ainsi que de la salinité de l'électrolyte. Dans le cas d'échantillons propres (sans argile), les lois de comportement existantes prédisent correctement la réponse à haute fréquence des matériaux étudiés en fonction de la teneur en eau. Des effets de polarisation d'interface s'ajoutent à la polarisation dipolaire des molécules d'eau pour des fréquences inférieures à 2 MHz. Ces effets sont d'autant plus importants que la teneur en argile augmente. Ils sont liés aux propriétés surfaciques de l'argile présente (capacité d'échange cationique, surface spécifique). Pour comprendre ces relations, des modèles complexes doivent être élaborés et testés. Nous avons privilégié une première approche qui consiste à établir une relation empirique afin de retrouver la valeur de la permittivité diélectrique liée à la polarisation dipolaire (permittivité HF due à la rotation des molécules d'eau). En parallèle a été développé un prototype d'appareil électromagnétique de mesure in situ, de type Slingram, composé de deux bobines perpendiculaires (PERP), qui fonctionne à 1,56 MHz. Après avoir étalonné l'appareil et testé sa réponse à des objets de propriétés connues, plusieurs sites d'expérimentation ont été choisis : des alluvions sableuses à l'INRA d'Orléans (Centre), des luvisols limoneux de l'INRA d'Estrée-Mons (Picardie) et des calcaires fracturés du site de Grand (Vosges). Ils présentent tous des natures de sols différents ; et des contrastes de teneurs en eau, naturels ou artificiels y ont été mesurés. Si, dans le cas des alluvions sableuses, les permittivités diélectriques sont du même ordre de grandeur que les permittivités HF et permettent de retrouver la teneur en eau du sol à partir des modèles théoriques existant, il n'en est pas de même pour des sols plus riches en argile. Les relations empiriques établies à l'aide des expérimentations en laboratoire permettent alors une première estimation de la teneur en eau, conforme aux teneurs en eau massiques des échantillons prélevés.

Cette Habilitation à Diriger des Recherches fait le point sur mes activités passées dans les domaines de la tectonique active, géomorphologie et datations 10Be. Ces activités disparates sont centrées autour d'un même thème qui est la signification du relief au dessus des zones de subduction. Au cours de ces années je me suis progressivement intéressé aux mouvements verticaux. Ceux-ci sont particulièrement bien enregistrés par les morphologies côtières et en particulier par les terrasses marines. Ces dernières se sont formées au cours des précédentes périodes interglaciaires, lorsque le niveau de la mer était haut, à peu près au même niveau que l'actuel. Ainsi, ces terrasses contiennent une information de qualité sur le soulèvement à l'échelle de quelques centaines de milliers d'années. Dans le détail, les séquences de terrasses montrent des variations morphologiques qui renferment probablement plus d'informations que leur altitude seule. C'est pourquoi je me suis intéressé à des cas de terrasses en formation : les plateformes côtières, via l'étude de terrain mais aussi la modélisation. J'ai par ailleurs travaillé sur la signification des concentrations en 10Be, et en particulier mis au point une méthode de quantification long-terme du recul des falaises. Ce mémoire se termine sur les perspectives d'avenir en continuité avec ces travaux sur la formation des terrasses, la vitesse d'évolution des plateformes ainsi que sur leurs caractéristiques en fonction des contextes géodynamiques.

Le marché de la pierre de construction représente une activité économique mondiale importante. Les échanges internationaux favorisent de potentielles substitutions, parfois frauduleuses, entre des pierres aux faciès comparables mais de provenances diverses. Il devient donc important de développer un outil de certification de la provenance des pierres de construction. Cette étude développe différentes méthodes d'identification (pétro-géochimie et/ou magnétisme) de pierres de construction de types granitique (Silverstar du Tarn, rose de La Clarté, gris bleu de Louvigné et bleu de Lanhélin) et sédimentaire (calcaire de Comblanchien) afin de pouvoir en certifier la provenance. De manière univoque, elle individualise ces pierres de construction françaises de potentiels concurrents commerciaux. Notre analyse utilise aussi des outils statistiques pour rendre compte de la variabilité intrinsèque des roches étudiées, quantifier la dissimilarité de roches de provenances différentes ou en établir, le cas échéant, la similarité. Au préalable l'échelle spatiale à laquelle une pierre de construction peut être considérée comme homogène est vérifiée. Des échantillons issus de carrières différentes sont comparés. Nos résultats et l'approche générale illustrée dans ce travail ne sont pas seulement d'intérêts pour les professionnels mais pour toute personne utilisant les pierres de construction à des fins industrielles. Cette étude montre qu'une caractérisation approfondie des pierres de construction peut en permettre la distinction en termes de provenance, et donc offre la possibilité de leur attribuer une carte d'identité, et si nécessaire un label commercial.