Les aquifères karstiques constituent l'essentiel des ressources en eau du pourtour Méditerranéen. Au delà de la zone saturée de ces systèmes, la zone non saturée constitue une entité importante quant au transfert et au stockage temporaire de l'eau. La structure et le fonctionnement de cette entité du karst sont complexes et mal connus à cause du processus de karstification qui hiérarchise la circulation hydrodynamique en son sein. Dans cette étude nous utilisons des méthodes géodésiques et hydrogéologiques afin d'observer de façon directe et de quantifier les processus de transfert et de stockage de la zone non saturée. Les mesures gravimétriques en surface et profondeur sur différents systèmes karstiques ont permis de quantifier les variations du stockage en eau saisonnier dans la zone non saturée. La gravimétrie différentielle montre que ces variations sont réparties dans les quelques premiers mètres de la zone non saturée. Il semble, de plus, que la capacité de stockage de la zone non saturée soit fonction de la lithologie de l'encaissant (calcaire et dolomie). Enfin, des observations directes du flux d'eau en zone non saturée ont été effectuées par le biais de mesures dans différentes cavités. Nous montrons que certaines caractéristiques du flux d'eau en zone non saturée, comme la présence d'un flux bipolaire lent-rapide, peuvent être généralisées quelques soit la profondeur de mesure et la surface d'impluvium. La modélisation de ce flux, à partir d'un schéma simple, a montré la complexité des processus de transfert engendrés au sein de la zone non saturée. / Karst systems are the most important groundwater reservoirs in Mediterranean area. Beyond saturated zone, unsaturated zone constitute an entity where water storage and transfer play a major role in a global behavior of spring water. This structure and functioning are complex and not well known because of the high heterogeneity created and organized by groundwater flow. In this study, we used gravimetric and hydrogeologic method to directly observe and quantify water storage and transfer processes in unsaturated zone of karst system. Surface to depth gravity measurements on several karst system allow quantifying seasonal water storage changes in unsaturated zone. We show that these variations occur in the first tens meters of unsaturated zone. Furthermore, water storage capacity of unsaturated zone seems to be influenced by lithology (limestone or dolomite). Finally, direct observations of water flow in unsaturated zone are done using caves access. We show that some characteristics of unsaturated water flow such as bipolar slow-fast flow can be widespread for whatever depth of measurement and recharge area. Model of this flow, using a simple modeling scheme, show the complexity of water transfer processes in unsaturated zone of karst system.

Les techniques actuelles de mesure de pesanteur terrestre, depuis les observations spatiales jusqu'aux mesures ponctuelles en surface, couvrent un large spectre de longueurs d'onde du champ de gravité de la Terre. Ces données acquises permettent non seulement la réalisation de modèles étendus de géoïde, mais aussi une meilleure compréhension des transferts de masse affectant la Terre.<br />Cependant, les couvertures spatiale et spectrale des mesures gravimétriques ne sont pas homogène à la surface de la Terre et certaines régions difficiles d'accès sont quasiment vierges de toute mesure. De plus, la gamme des longueurs d'onde intermédiaires (10/100 km) est mal couverte par la gravimétrie terrestre ou spatiale. La gravimétrie mobile depuis un véhicule terrestre, un bateau ou un avion apparaît donc comme la technique complémentaire capable de pallier les insuffisances des techniques gravimétriques actuelles.<br />C'est dans l'objectif de développer un système de gravimétrie mobile autonome que ce travail a été entrepris. Le système «Limo-G» (LIght MOving Gravimetry system) est un gravimètre vectoriel absolu, composé de trois accéléromètres et d'un système GPS à quatre antennes, embarquable dans toutes sortes de véhicule, fiable, ergonomique et peu onéreux.<br />Un travail métrologique effectué au laboratoire et sur le terrain a permis de concevoir et de tester une méthode d'étalonnage de l'instrument. Puis, une méthode originale de traitement conjoint des données accélérométriques et des données GPS a été développée et testée à partir de simulations issues de données réelles acquises lors d'un levé expérimental en mer. Les résultats de ces simulations indiquent que les performances du système Limo-G permettent de mesurer la pesanteur au milligal. Les améliorations nécessaires pour accroître la précision ont aussi été mises en évidence et validées par simulation.

Iran couvre une grande superficie en longitude entre les méridiens 44°E et 64°E et en latitude entre les parallèles 25°N et 40°N. la cartographie de champ de pesanteur sur l'Iran est de première importance pour des considérations d'ordre géodésique, géophysique et géodynamique. Dans cette thèse, les mesures de pesanteur sont utilisées pour déterminer le géoïde gravimétrique sur l'Iran. Ce géoïde est couplé à la mesure de hauteur GPS et d'altitude (nivellement) pour réaliser une surface verticale opérationnelle sur le territoire Iranien. La contribution aux principaux travaux géodésiques et gravimétriques réalisés ces dernières années porte sur l'établissement: 1) du réseau national de gravimétrique absolu de l'Iran (NAGNI09), 2) de la ligne nationale d'étalonnage de gravimétrique de l'Iran (NGCLI10), 3) du réseau multi- observations géodésiques et gravimétriques de l'Iran (MPGGNI10). Le réseau gravimétrique absolu, comporte 24 stations où la mesure de la pesanteur a été réalisée à l'aide de gravimètres FG5, avec une précision meilleure que 5 Gal. La répétition des observations sur deux sites entre 2000 et 2007 met en évidence des variations inter-annuelles de la pesanteur en relation avec l'évolution du contenu en eau du sous-sol et (ou) la déformation tectonique. Le réseau a servi de point d'appui pour la réalisation du réseau géodésique et gravimétrique MPGGNI10 de maille 55 km sur lequel a été mesuré la pesanteur à l'aide des gravimètres relatifs CG5 et CG-3/M, la hauteur de GPS et l'altitude avec une précision respectivement de 0.010 mGal, 0.03 m et . La technique de retrait- restauration couplée à la méthode de condensation de Helmert a permis de calculer un nouveau modèle de géoïde gravimétrique, IRGeoid10, avec une précision absolue et relatif respectivement de l'ordre de 0.26 m et 2.8 ppm. Ce géoïde est ajusté aux points GPS nivelés pour définir un nouveau référentiel des altitudes sur l'Iran. / Iran covers a large area limited in longitude by the meridians 44°E and 64°E and in latitude by the parallels 25°N and 40°N. Mapping a new gravity field over Iran is the first important data for geodetic, geophysical and geodynamical considerations. In this thesis, the gravity measurements are used to determine the gravimetric geoid over Iran. This geoid is coupled with the GPS height and altitude (levelling) to realize an operational vertical surface at the territory of Iran. The contribution of the principal geodetic and gravimetric works realized in recent years are the establishment of: 1) the national absolute gravity network of Iran (NAGNI09), 2) the national gravity calibration line of Iran (NGCLI10) and 3) the multi-observations geodetic and gravimetric network of Iran (MPGGNI10). The absolute gravity network consists in 24 stations where the gravity measurement has been realized with the help of gravimeters FG5, with a precision better than 5 Gal. The repetition of the observations at two stations between 2000 and 2007 makes obvious the inter-annual variations of gravity in relation of the amount of underground water changes and (or) tectonic deformation. The absolute gravity network has served the base stations for the realization of the MPGGNI10 geodetic and gravimetric network with a mesh of 55 km, at which the gravity is measured with the help of relative gravimeters CG-5 and CG-3/M, the GPS height and the altitude with a precision of 0.010 mGal, 0.03 m, and respectively. The remove-restore technique coupled with the Helmert's condensation method is chousen to compute a new gravimetric geoid model, IRGeoid10, with a absolute and relative precision of the order of 0.26 m and 2.8 ppm respectively. The gravimetric geoid is adjusted at the GPS/levelling points to define new vertical reference surface over Iran

Nous avons mené cet exercice au travers d’une étude pluridisciplinaire combinant des couvertures géophysiques (gravimétrie et aéromagnétisme) réactualisées ou nouvelles, des informations géologiques acquises sur le terrain ou sur carottes de forages pétroliers et un catalogage de propriétés pétrophysiques des roches. Sur une large moitié sud du Bassin parisien, du Massif armoricain aux Vosges, nous avons : i) réalisé un schéma structural du substratum du Bassin parisien, ii) proposé une carte géologique détaillée, sous couverture, iii) réactualisé la carte de zonation de la chaîne varisque en France, intégrant de nouvelles datations sur carottes de forages ayant atteint le substratum. D’après ces résultats, la suture éo-varisque se poursuit d’ouest en est, de la partie sud-armoricaine jusqu’au nord des Vosges, le domaine nord-armoricain se prolonge jusqu’à la faille de Bray et la zone saxothuringienne est continue d’est (Vosges) en ouest (Manche). Egalement, à l’aide de modélisations conjointes de gravimétrie et magnétisme, intégrant l’information recueillie dans le forage de Couy, la nature du substratum environnant désormais mieux connu, ainsi que les contraintes pétrophysiques nouvelles, nous proposons que la source de l’AMBP soit un corps plutonique intermédiaire à acide. Pour finir, en compilant la bibliographie et nos résultats, nous proposons que ce corps se soit mis en place entre le Viséen inférieur et le Viséen supérieur, dans une ouverture NNE-SSW selon un raccourcissement maximal NNW-SSE. / Mapping the cadomian/variscan basement buried beneath the Paris basin is a key objective to improve the geological knowledge of the Variscan belt in West Europe. We address this objective using a combination of new and updated potential field data (aeromagnetic and gravity), as well as information provided by field studies and drill cores, and a library of petrophysical characteristics measured on field samples. On a large southern half of the Paris basin, from the Armorican Massif to the Vosges: i) we made a new interpretative structural sketch map under the sedimentary cover, ii) we proposed a detailed geological map of the Paris basin substratum, iii) we updated the tectonic sketch of the Variscan belt in France, supported by new geochronological determinations on deep drill cores. Altogether, these results show that the eo-Variscan suture is extended from the southern part of the Armorican Massif to the northern Vosges. Eastwards, the North-Armorican domain reaches the Bray fault; and the Saxothuringian zone is extended from the Vosges to the English Channel. In addition, based on combined aeromagnetic and gravity 2D modelling constrained by Couy deep borehole data, the PBMA source is interpreted as a deep intermediate-alkaline plutonic body. Finally, compiling bibliographic information together with our new results, we hypothesize that this magmatic body could have been emplaced in the early to late Visean period, during a NNE-SSW opening under a NNW-SSE shortening.

Au cours du dernier siècle, plusieurs séismes majeurs ont affecté l'Himalaya. Cependant, la taille maximale de ces événements et la probabilité d'occurrence de méga-séismes avec des magnitudes proches de 9 sont toujours matière à débat. L'étude de la segmentation de l'arc Himalayen est donc primordiale afin de comprendre les mécanismes qui contrôlent ces séismes ainsi que leur extension spatiale. La compréhension du cycle sismique en Himalaya est aujourd'hui essentiellement basée sur des études menées au Népal central, ce qui limite notre connaissance de son fonctionnement tri-dimensionnel. Ce travail de thèse permet d'apporter de nouvelles contraintes sur les variations latérales de la déformation crustale dans la zone comprise entre l'ouest Népal (~80°E) et l'est Bhoutan (~92°E). La première partie de cette thèse est consacrée à l'étude des variations latérales de la structure lithosphérique. Quatre campagnes gravimétriques ont été réalisées entre 2010 et 2012. Elles ont permis d'étendre le réseau gravimétrique Népalais jusqu'à l'ouest Népal et d'établir le premier réseau gravimétrique Bhoutanais. Ce nouveau jeu de données, combiné avec les données sismologiques disponibles, permet de contraindre des modèles thermo-mécaniques de la flexure de la plaque Indienne sous l'Himalaya. Les résultats suggèrent qu'il n'existe pas de variations latérales majeures du comportement mécanique de la lithosphère Indienne entre le centre et l'est Népal contrairement au Bhoutan où une rigidité flexurale plus faible est nécessaire pour expliquer les données.La deuxième partie de cette thèse est consacrée à l'étude des variations latérales des déformations récentes dans le prisme Himalayen. Les études paléosismologiques menées depuis 15 ans en Himalaya ont permis d'étendre le catalogue des événements majeurs sur le dernier millénaire. En intégrant à la fois les contraintes disponibles sur ces séismes majeurs et le chargement séculaire, nous étudions les variations spatiales et temporelles de l'état de contrainte du chevauchemlent Himalayen principal le long de l'arc. Nos calculs montrent que même dans la cas où un méga-séisme se serait produit en ~1100 en Himalaya, le chargement séculaire a quasiment compensé la chute de contrainte associée. Les modèles montrent aussi que le séisme du Shillong ne joue pas un rôle majeur sur les contraintes accumulées sur le chevauchement Himalayen depuis 1897 au niveau du Bhoutan. Enfin, nous présentons les premières contraintes sur la tectonique active au Bhoutan. L'étude morphotectonique réalisée au centre du Sud-Bhoutan a permis de montrer que la majeure partie de la déformation Holocène est accommodée au niveau du front Himalayen, comme au Népal. Nous montrons aussi des évidences de ruptures de surface avec des décalages verticaux de plusieurs mètres associés à deux séismes majeurs sur le dernier millénaire. Ces évidences remettent en cause l'interprétation d'un méga-séisme en ~1100. / Several major earthquakes have struck the 2500 km long Himalayan range during the past century. However, both the maximum size of such earthquakes and the probability of occurrence of a magnitude 9 megaquake are still a mater of debate. In this framework it is of key importance to investigate the lateral segmentation of the Himalayan arc in order to deepen our understanding of the mechanisms that control such events as well as their extents. Besides the remarkably uniform tectonic set-up and the fact that most studies on the Himalayas have focused on Central Nepal, several indications of along strike changes can be gathered.This thesis brings new constraints on lateral variations of crustal structure and of deformation between western Nepal (~80°E) and eastern Bhutan (~92°E).The first part of this thesis focuses on the lateral variations of crustal structure. In the past three years we have carried out four land gravity campaigns. We have established 366 new gravity measurement points to fill the data gaps in central and western Nepal as well as in Bhutan. This new dataset, together with available seismological data, is used to constrain thermo-mechanical models of the bending of the India plate underneath the Himalayas. While the inferred crustal geometry does not exhibit major along-strike variations over central and eastern Nepal, the shorter wavelength flexure of the lithosphere in Bhutan is associated with a weaker upper mantle rheology.The second part focuses on the lateral variations of recent deformation in the Himalayas. Fifteen years of paleoseismological investigations have extended the catalogue of major Himalayan events over the last millennium. Combining these information with secular loading, we assess the spatial and temporal stress changes on the Main Himalayan Thrust along the orogen over the last nine centuries. Our calculations indicate that inter-seismic loading has now nearly overcome the Coulomb stress decrease caused by the great ~1100 medieval earthquake. Our results also point out that the 1897 Shillong plateau earthquake does not have a major influence on the stress accumulated on the Main Himalayan Thrust since then. In order to better characterize active tectonics in the Bhutan Himalayas where no studies were done so far, we carried out a morphotectonic analysis in the south-central part of the kingdom. We show that the same amount of Holocene deformation is accommodated on the frontal thrust in Bhutan as on the neighbouring portions of the Himalayan arc. We also find evidences for two major (M>8) earthquakes on this thrust in the Bhutan Himalayas during the last millennium. Our results therefore show that Bhutan cannot be considered as a seismic gap. They also challenge the interpretation of one single ~1100 medieval mega-event that would have ruptured the Main Frontal Thrust from central Nepal to eastern Bhutan.

Les aquifères karstiques constituent l'essentiel des ressources en eau souterraine sur le pourtour de la Méditerranée. Leur structure et leur fonctionnement sont complexes et mal connus à cause du processus de karstification qui hiérarchise la circulation hydrodynamique en leur sein. La localisation du stockage de l'eau et de ses variations au sein de tels aquifères est difficile à mesurer et à modéliser. Nous utilisons les outils de la géodésie pour apporter de nouvelles contraintes sur le fonctionnement hydrologique et les variations de stock d'eau en milieu karstique. Ainsi, le système karstique du Durzon, (Causse du Larzac, France) fait l'objet d'un suivi par plusieurs méthodes géodésiques incluant: - Des mesures de gravimétrie absolue en trois sites répétées de manière mensuelle ; - Des mesures répétées de gravimétrie relative entre la surface et la profondeur du karst ; - Des mesures répétées de microgravimétrie ; - Des mesures inclinométriques longue base en deux sites. Nous démontrons que le couplage entre la gravimétrie absolue et relative permet d'appréhender les variations de stock d'eau sur le bassin d'alimentation étudié. Les variations inter-anuelles observées sont de l'ordre de 400 mm, et nous montrons que celles-ci s'opèrent dans certaines zones au sein des horizons supérieurs du karst : épikarst et zone d'infiltration. Une forte hétérogénéité spatiale des variations de stockage est aussi mise en évidence à l'échelle du bassin d'alimentation. Les mesures inclinométriques permettent de suivre les phénomènes d'infiltration rapide et lente dans la masse karstique fracturée. Le suivi géodésique de ce système karstique permet d'apprécier quantitativement et de façon intégrée son fonctionnement hydrologique et pourra être transposé pour l'étude d'autres aquifère karstiques.

La mesure des dérivées spatiales du champ de pesanteur, ou gradiométrie en gravimétrie, est née au début du XXe siècle avec la balance de torsion d'Eötvös, premier gradiomètre. Utilisée avec succès pour la prospection géophysique, cette technique fut cependant délaissée pour un temps au profit de la gravimétrie, plus facile et moins coûteuse. Cependant, les développements instrumentaux en gradiométrie ont continué, et les gradiomètres actuels permettent la mesure simultanée des gradients de pesanteur dans trois directions indépendantes de l'espace. La gradiométrie en gravimétrie connaît ainsi une renaissance, et l'Agence Spatiale Européenne lancera en 2008 le satellite GOCE, avec à son bord le premier gradiomètre spatial, permettant la cartographie globale de la pesanteur avec une résolution spatiale sans précédent. Cette étude est consacrée au signal de gradiométrie en gravimétrie, de l'acquisition des données à leur interprétation. Plus spécifiquement, dérivant d'un même potentiel, l'accélération de la pesanteur et ses dérivées peuvent être considérées comme des mesures interdépendantes d'une même quantité. Nous avons élaboré une méthode permettant, en exploitant cette redondance, de réduire le bruit dans les données de gradiométrie et, plus efficacement, dans celles de gravimétrie. Nous utilisons également une propriété spécififique des gradients de pesanteur, l'existence d'invariants scalaires combinant les différents gradients, que nous relions à la géométrie des sources à l'origine des anomalies de pesanteur. Nous avons ainsi développé une méthode, complémentaire à la déconvolution d'Euler des données de gravimétrie, qui améliore la localisation de ces sources.

Les techniques de télédétection spatiale constituent un apport majeur pour l'étude des variations de masses d'eau dans les grands bassins fluviaux, en permettant un suivi homogène de ces fluctuations dans l'espace et dans le temps. <br />L'objectif de cette thèse a été de développer de nouvelles applications hydrologiques au moyen des mesures spatiales acquises par différents types de mission satellitaire : altimétrie radar, imagerie satellitaire, gravimétrie spatiale. L'altimétrie spatiale offre la possibilité d'étudier les variations de niveau d'eau des grands fleuves, des lacs et des zones d'inondation, garantissant ainsi une surveillance continue et globale des eaux de surface. Elle donne aussi accès à des produits hydrologiques nouveaux comme le profil hydrologique ou la pente des fleuves. Elle permet en outre de définir des réseaux limnimétriques nivelés, dont les stations peuvent être définis tant sur les fleuves que sur les zones d'inondation, complémentaires des réseaux in-situ. Combinée à l'imagerie spatiale, l'altimétrie satellitaire a été utilisée pour déterminer les variations de volume d'eau dans les grands bassins fluviaux. Ces paramètres revêtent, en effet, une importance fondamentale pour les hydrologues car le premier est à la base des études hydrodynamiques et le second apporte des contraintes sur la répartition des masses d'eau entre zones inondées et réseau hydrographique, avec des applications au transport des sédiments et à la disponibilité des ressources en eau à l'échelle régionale. Des exemples d'utilisation de ces techniques sont présentés pour les bassins amazoniens et du Mékong. La mission de gravimétrie spatiale GRACE, lancée en mars 2002, fournit, quant à elle, les variations spatio-temporelles des stocks d'eaux continentales (eau des sols et manteau neigeux) et de paramètres hydrologiques dérivés comme l'évapotranspiration. Une analyse de l'évolution des stocks d'eau et de neige est présentée à partir des premiers géoïdes mensuels issus de la mission GRACE, aux échelles globale et régionale, ainsi que le calcul du paramètre d'évapotranspiration, à l'échelle du bassin versant. Ces résultats sont comparés à la variation des volumes d'eau de surface obtenue précédemment pour le bassin du Mékong.

L'Iran couvre une grande superficie en longitude entre les méridiens 44°E et 64°E et en latitude entre les parallèles 25°N et 40°N. la cartographie de champ de pesanteur sur l'Iran est de premiére importance pour des considérations d'ordre géodésique, géophysique et geodynamique. Dans cette thèse, les mesures de pesanteur sont utilisées pour déterminer le géoïde gravimétrique sur l'Iran. Ce géoïde est couplé à la mesure de hauteur GPS et d'altitude (nivellement) pour réaliser une surface verticale opérationnelle sur le territoire Iranien. La contribution aux principaux travaux géodésiques et gravimétriques réalisés ces dernières années porte sur l'établissement: 1) du réseau national de gravimétrique absolu de l'Iran (NAGNI09), 2) de la ligne nationale d'étalonnage de gravimétrique de l'Iran (NGCLI10), 3) du réseau multi- observations géodésiques et gravimétriques de l'Iran (MPGGNI10). Le réseau gravimétrique absolu, comporte 24 stations où la mesure de la pesanteur a été réalisée à l'aide de gravimètres FG5, avec une précision meilleure que 5μ Gal. La répétition des observations sur deux sites entre 2000 et 2007 met en évidence des variations inter-annuelles de la pesanteur en relation avec l'évolution du contenu en eau du sous-sol et (ou) la déformation tectonique. Le réseau a servi de point d'appui pour la réalisation du réseau géodésique et gravimétrique MPGGNI10 de maille 55 km sur lequel a été mesuré la pesanteur à l'aide des gravimètres relatifs CG5 et CG-3/M, la hauteur de GPS et l'altitude avec une précision respectivement de 0.010 mGal, 0.03 m et 3mm km . La technique de retrait- restoration couplée à la méthode de condensation de Helmert a permis de calculer un nouveau modèle de géoïde gravimétrique, IRGeoid10, avec une précision absolue et relatif respectivement de l'ordre de 0.26m et 2.8 ppm. Ce géoïde est ajusté aux points GPS nivelés pour définir un nouveau référentiel des altitudes sur l'Iran.

Les volcans sont des édifices complexes et hétérogènes. Imager leur structure demande des méthodes adaptées au contexte ainsi qu'à la nature, à la qualité et à la quantité des données disponibles.Ce travail porte sur l'imagerie géophysique de l'île volcanique de Basse-Terre en Guadeloupe, par inversion de données gravimétriques et magnétiques, par tomographie d'ondes de surface extraites de corrélations de bruit sismique et par tomographie de temps d'arrivée des séismes de l'arc des Petites Antilles. Nous obtenons ainsi les premiers modèles tridimensionnels de l'île de Basse-Terre en termes de densité, aimantation et vitesses sismiques. La gravimétrie et le magnétisme conduisent à des modèles de densité et d'aimantation avec une très bonne résolution latérale dans les premiers kilomètres, mais avec une résolution verticale très faible. Les ondes de surface issues des corrélations de bruit présentent une bonne résolution verticale dans les premiers kilomètres, mais résolvent mal les variations latérales de vitesses. La tomographie de temps d'arrivée permet d'imager les vitesses sismiques jusqu'à des profondeurs élevées, avec une bonne résolution latérale, mais une mauvaise résolution verticale. Nous proposons une méthode pour inverser conjointement les données sismologiques et gravimétriques pour produire un modèle de vitesses et de densité qui bénéficie des résolutions complémentaires des différents jeux de données.Nos résultats soulignent une tendance NNO-SSE, avec des vitesses et des densités qui diminuent vers le sud en moyenne, conformément à la migration du volcanisme du nord vers le sud au cours de l'histoire géologique de l'île. Nous observons de faibles densités, vitesses sismiques et aimantations au niveau du dôme du volcan de la Soufrière ainsi que dans la région de Bouillante sur la côte ouest, reflétant la forte altération hydrothermale de ces zones. Deux anciens centres éruptifs majeurs, caractérisés par de fortes densités et vitesses sismiques, sont identifiés : l'un au niveau de la vallée de Beaugendre, l'autre au nord au niveau de la vallée de Petite Plaine. / Volcanoes are complex and heterogeneous edifices. Imaging their structures requires methods that are adapted to the context and to the nature, quality and quantity of available data. This work focuses on the geophysical imaging of the volcanic island of Basse-Terre in Guadeloupe, by performing inversions of gravimetric and magnetic data, ambient noise based surface wave tomography, and travel time tomography from earthquakes of the Lesser Antilles arc. We therefore obtain the first three-dimensional models of Basse-Terre island in terms of density, magnetisation and seismic velocities. Gravimetry and magnetism lead to density and magnetisation models with a very good lateral resoltuion within the first kilometres but with very low vertical resolution. Surface wave tomography from noise correlations shows a good vertical resolution within the first kilometres, but poorly resolves the lateral velocity variations. Travel time tomography allows imaging seismic velocities up to great depths with a good lateral resolution but a poor vertical resolution. We propose a method to invert jointly the seismological and gravimetric data in order to produce a velocity and density model that takes advantage of the complementary resolutions associated with the different datasets.Our results outline a NNW-SSE trend, with velocities and densities roughtly decreasing southwards, in accordance with the southward migration of volcanism across the geological history of the island. We observe low densities, seismic velocities and magnetisations in the area of the dome of the Soufrière volcano as well as in the Bouillante area along the western coast, reflecting the high hydrothermal alteration of these areas. Two major old eruptive centres are identified in the areas of Beaugendre Valley and Petite Plaine Valley.