# Montagne Pelée volcano (Martinique, in the French Lesser Antilles) hydrogeological system revealed by high-resolution helicopter-borne electromagnetic imagery

Benoit Vittecoq, Pierre-Alexandre Reninger, Vincent Bellier, Anne-Lise Taïlamé, Laureen Nascimento, Emile Gros, Sophie Violette

<figure><img src="image_1.png" /></figure>

## To cite this version:

Benoit Vittecoq, Pierre-Alexandre Reninger, Vincent Bellier, Anne-Lise Taïlamé, Laureen Nascimento, et al.. Montagne Pelée volcano (Martinique, in the French Lesser Antilles) hydrogeological system revealed by high-resolution helicopter-borne electromagnetic imagery. Hydrogeology Journal, 2023, pp.1-23. 10.1007/s10040-023-02642-5. hal-04155223

HAL Id: hal-04155223

https://brgm.hal.science/hal-04155223

Submitted on 7 Jul 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Montagne Pelée volcano (Martinique, in the French Lesser Antilles) hydrogeological system revealed by high-resolution helicopter-borne electromagnetic imagery

Benoit Vittecoq<sup>1,2 ID</sup> · Pierre-Alexandre Reninger<sup>3</sup> · Vincent Bellier<sup>1</sup> · Anne-Lise Taïlamé<sup>1</sup> · Laureen Nascimento<sup>1</sup> · Emile Gros<sup>1</sup> · Sophie Violette<sup>2,4</sup>

Received: 10 October 2022 / Accepted: 9 May 2023
© The Author(s) 2023

## Abstract

Montagne Pelée, on the French island of Martinique, eastern Caribbean Sea, has been one of the deadliest volcanoes in the world, with 30,000 victims following the 1902 eruption. Thousands of people still live nearby, and this volcano is a strategic “water tank” for Martinique Island, providing 40% of the island’s water supply. This research aimed to better understand its hydro-geological functioning and the relationship with its complex volcanological evolution, taking advantage of a high-resolution helicopter-borne geophysical survey correlated with hydrogeological data from the boreholes and springs databases. Electromagnetic data, correlated with hydrogeological data, allowed for the identification of unsaturated zones, aquifers, and seawater intrusions, as well as the main geological units. In addition, data synthesised from pumping tests revealed that the older the unconsolidated pyroclastic deposits, the lower their hydraulic conductivity. The structural asymmetry between the northeastern and southwestern volcano flanks impacts its hydrogeological functioning. Consequently, the Montagne Pelée hydrogeological conceptual model is marked by several distinguishable aquifers. The upper perched aquifer within recent lava domes is directly involved in, and impacted by, phreatic eruptions, and it supports low flowrate springs. The remaining effective rainfall infiltrates to depth and recharges the hydrothermal system through vertical fractures. The other aquifers are categorized into three groups: northeastern, southeastern and southwestern flank aquifers. This research is a new step toward a better understanding of the Lesser Antilles volcanoes and more broadly of the central and proximal parts of the andesitic active volcanoes.

**Keywords** Volcanic aquifer · Electromagnetic geophysical methods · Conceptual model · Island hydrogeology · Martinique (France)

## Introduction

Over 29 million people live on the flanks of active volcanoes (Brown et al. 2017). Among the essential needs of a population, access to water is crucial for both drinking purposes and irrigation, especially on volcanic islands. Furthermore, surface water and groundwater are often affected

by volcanic activities or are involved in the eruptive processes—for example, groundwater plays a key role in phreatic eruptions (steam-driven explosions). Regarding all these concerns, it is necessary to gain a deep understanding of the geological structure and hydrogeological functioning of active volcanoes.

As shown by Cabrera and Custodio (2019), each island or each volcano has its own hydrogeological characteristics inherited from its geological history, structure and erosion processes. Nonetheless, they also share some characteristics. Thus, for each volcanic island, or for each active volcano, it is necessary to define adequate hydrogeological conceptual models at appropriate scales in order to allow an adaptive management strategy for water resources and to better constrain the interaction between volcanic activities and groundwater.

Historically, hydrogeological functioning of volcanic islands has been synthetized using two conceptual models,

<figure><img src="image_3.png" /></figure>

<footnote>Benoit Vittecoq
b.vittecoq@brgm.fr</footnote>

<footnote n="1">1 BRGM, 97200 Fort de France, Martinique</footnote>

<footnote n="2">2 ENS-PSL Research University and CNRS, UMR.8538 – Laboratoire de Géologie, 24 rue Lhomond, 75231 Paris, France</footnote>

<footnote n="3">3 BRGM, F-45060 Orléans, France</footnote>

<footnote n="4">4 Sorbonne Université, UFR.918, F75005, Paris, France</footnote>
the Hawaiian model (Meinzer 1930; Peterson 1972; Macdonald et al. 1983); and the Canary Islands model (Custodio 1975; Custodio et al. 1988), both for basaltic islands. More recently, a conceptual model of old basaltic islands has also been proposed by Vittecoq et al. (2014), and a relationship between these models, with an evolution from recent Hawaiian type islands (<1 Ma) toward Canarian-type model islands due to weathering processes, has been proposed by Violette et al. (2014). None of these models are appropriate for describing the water resources of andesitic volcanoes.

Vessel and Davis (1981) and Bogie and Mackenzie (1998) have described the geological structure of andesitic volcanoes. Their conceptual model, based on the repartition of volcanoclastic deposits as a function of the distance to the event (central, proximal, median and distal) with heterogeneous facies (lava, pyroclastic deposits, tuff, tephra, air fall, breccia, debris avalanches, lahars and alluvial deposits), allows a better understanding of the hydrogeological functioning (Selles 2014; Selles et al. 2015). For the Lesser Antilles in the Caribbean Sea, Vittecoq et al. (2015), Hemmings et al. (2015) and Vittecoq et al. (2019) have described the hydrogeological functioning of central and proximal parts of mainly fissured and fractured andesitic lava aquifers dating from 1 to 15 Ma. Nevertheless, Lesser Antilles volcanoes are also characterized by young active pyroclastic edifices such as Soufriere Hills Volcano, Montserrat (Wadge et al. 2014), Soufriere, St Vincent (Robertson 1995; Pyle et al. 2018), Soufriere, Guadeloupe (Legendre 2012) and Montagne Pelée, Martinique (Westercamp and Traineau 1983; Westercamp et al. 1990; Michaud-Dubuy 2019).

Montagne Pelée is an iconic volcano, cradle of modern volcanology and one of the deadliest volcanoes of the world with more than 30,000 victims following the 1902–1905 eruption. This volcano is located in the north of Martinique Island (Fig. 1), an andesitic volcanic island of the Lesser Antilles volcanic arc, induced by the subduction zone between the Atlantic Plate and the Caribbean Plate (Westercamp et al. 1989). Thousands of people still live nearby and this volcano is a strategic “water tank” for Martinique Island. River intakes, springs and boreholes provide 40% of the water demands of the island (around 100,000 m³/day). One of its particularities is the presence of porous and permeable pyroclastic deposits superimposed over several tens to hundreds of meters.

Helicopter-borne electromagnetic methods have demonstrated their applicability with respect to basaltic islands (d'Ozouville et al. 2008; Auken et al. 2009; Pryet 2011; Pryet et al. 2011, 2012), active volcanos (Dumont et al. 2019), and andesitic islands (Vittecoq et al. 2015, 2019), especially because of their ability to measure resistivity variations within the first few hundred meters of depth. Nevertheless, to the authors' knowledge, there has been no published research about pyroclastic volcanos and associated hydrogeological

functioning. For Montagne Pelée, being an active volcano, the relationship with the hydrothermal system is also of interest and a better description of groundwater circulations should assist towards a better understanding of the eruptive processes.

This research takes advantage of a unique high-resolution helicopter-borne geophysical survey (Deparis et al. 2014) conducted over the Montagne Pelée (Fig. 2) correlated with geological and hydrogeological data from a database of boreholes and springs (Figs. 2 and 3) with the aim of deepening the comprehension of the structure and hydrogeological functioning of this active volcano and to define its hydrogeological conceptual model.

# Study area

## Geology

Located in the central part of the Lesser Antilles Archipelago (Fig. 1a), Martinique Island is a volcanic island particularly notorious for its active volcano, the Montagne Pelée. This volcano is particularly known for its series of violent, superficial and laterally directed explosions during the 1902–1905 dome-forming eruptions and especially the 1902 May 8th and August 30th events which killed about 30,000 people (Lacroix 1904).

Located in the northern part of the island (Fig. 1b), Montagne Pelée (1,397 m amsl) occupies an area of about 150 km². As described by Vincent et al. (1989), the Montagne Pelée is not a regular andesitic cone. Volcanic deposits are mainly found on its southwestern and northeastern flanks. To the southeast, Montagne Pelée products overlie the older Morne Jacob volcano (5.5–1.5 Ma, Germa et al. 2010). The history of the Montagne Pelée is a good example for demonstrating that eruptive-style changes occur during a volcano's lifetime within the same edifice. The first stage is marked by an effusive style with the dominant production of lava flows and domes, whereas the second and third stages are mainly explosive style, represented by pyroclastic products (Boudon and Balcone-Boissard 2021). The pyroclastic deposits of the Montagne Pelée mostly result from ash and pumice fallouts and diverse pyroclastic density currents (PDC, which are flows of hot gas, ashes and debris from the collapse of a lava dome or an eruptive column) produced during the successive dome-forming, subplinian and plinian eruptions.

Boudon and Balcone-Boissard (2021) have provided a new volcanic activity evolution description of the Montagne Pelée, defined in three stages (Fig. 3). The first stage, from 550–127 ka, is called the primitive Montagne Pelée. The dominant geological facies are andesitic lava flows and lava domes. This primitive volcano activity stopped at 127 ka (Germa et al. 2011a) following a supposed huge flank collapse (The Prêcheur event: 25 km³) that partly destroyed the southwestern flank of the volcano (Vincent et al. 1989; Le
<figure><img src="image_4.png" /><figcaption>**Fig. 1** a Location of the Lesser Antilles Archipelago and of Martinique Island (see arrow) and b location of Montagne Pelée on Martinique Island</figcaption></figure>

Friant et al. 2003). The remains of this volcanic edifice are observable in the northwestern flank and described as the Mont Conil geological unit in previous studies (Westercamp et al. 1989). From 127 to 36 ka, during the second stage, the dominant volcanic activity was andesitic lava dome-forming and consolidated pyroclastic deposits. The third stage (36 ka to present day) is divided in two parts and begins following a supposed flank collapse (36 ka: The Rivière Sèche event, 1.8 km³). The first part (36–25 ka) is dominated by numerous explosive eruptive events, mainly low silica subplinian and plinian deposits, filling the horseshoe-shaped structure induced by the flank collapse and then covering all the flanks of the volcano. The second part (25 ka to the present) diverges from the last one with an increase of silica content leading to volcanic products with an andesitic composition. The deposits result from alternating felsic plinian and subplinian events on one hand and dome-forming eruptions (recent lava domes and PDCs) on the other. Recently, Villemant et al. 2022 consider that the existence of two flank collapses (“Saint Pierre” and “Rivière Sèche”; Le Friant et al. 2003, Germa et al. 2011b) is questionable, as they are not recorded in their offshore cores. Pictures of geological outcrops of each stage are shown in Figures SI_1 and SI_2 in the electronic supplementary material (ESM).

## Fumaroles, warm springs, and the hydrothermal system

Fumaroles activity, summarized by Barat (1986), was observed as being associated with the 1792 and 1851 phreatic eruptions. After each crisis, fumaroles activity rapidly decreases. In 1889, fumaroles restarted at the summit until the 1902 eruptions. Between 1903 and 1929 fumaroles activity

did not stop and increased in 1929 before the eruption. After the 1929 eruption, the temperature of the fumaroles was 350 °C. Temperature rapidly decreased to stabilize at 100 °C, with chemical composition almost exclusively composed of CO₂. Then the fumaroles temperature continued to decrease until 1970 with the disappearance of the fumaroles of the Montagne Pelée. The last analysis (Fabre and Chaigneau 1960) showed a temperature of 80 °C and chemical composition close to the air (79% N₂, 20% O₂ and 0.5% CO₂).

The scarcity of thermal manifestations on the Montagne Pelée is supposed to be the result of permeable aquifers masking ascending fluids (Traineau et al. 1989). Seven groups of warm or thermal springs are observable on the flanks of the Montagne Pelée (Figs. 2 and 3): (1) the Chaude River springs (which appeared after the 1902 eruption according to Lacroix 1904) with current temperatures (recorded in 2022) between 24 and 32 °C (temperature declining over the past 50 years: 60–80 °C in the 1960s, 40–50 °C in the 1990s and the 2000s, and since then lower than 36 °C); (2) the Claire River springs (known before the 1792 phreatic eruption) with temperatures between 25 and 30 °C; (3) the Mitan-Picodo springs (also known before 1792), with temperature between 25 and 36 °C; (4) the Blanche River underwater springs (30–31 °C); (5) an underwater gas bubble emission zone between 2 and 8 m depth associated with several diffuse thermal water springs with temperature around 30 °C, located 1.5 km south from the Blanche River underwater springs; (6) the Grand Rivière spring, at an elevation of 865 m with 27 °C; and (7) the Des Pères River springs (20–23 °C). A small spring would also have been observed in the Prêcheur River but no more precise information is available. Finally, only one borehole has encountered thermal water, a 10-m-deep well
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Locations of flight lines, boreholes, freshwater springs, warm springs and rivers on the digital elevation model (DEM). Red lines correspond to the cross-sections (C1–C8) mentioned in section “Resistivity ranges of pyroclastic aquifers and the main geological units”</figcaption></figure>

approximatively located 100 m upstream of the Blanche River underwater springs, with groundwater temperature around 30–32 °C and a mean water electrical conductivity of 1,267 µS/cm, monitored by the volcanic and seismic observatory of Martinique (Fontaine et al. 2022).

Traineau et al. (1989), then Gadalia et al. (2014), proposed several insights of the hydrothermal system functioning. The heat source, estimated around 800–900 °C, could be a magma chamber located at 5–8 km depth beneath the summit. At least two geochemical zones of the hydrothermal reservoir have been discerned. The first, feeding the Rivière Chaude springs, has bicarbonate-sodium composition (HCO₃⁻-Na). The equilibrium temperature at depth, according to geo-thermometers (Na-K-Ca, Sr-K, Na-Li and Ca-K, Gadalia et al. 2014), would be between 180 and 200 °C. The second,

supplying the Mitan and the Picodo springs, has sodium bicarbonate-chlorinated composition (HCO₃⁻-Cl-Na), with a temperature of 155–180 °C at depth. The Blanche River underwater springs are expected to be the outlets of a lateral flow from the Chaude River springs through pyroclastic deposits (Traineau et al. 1989; Zlotnicki et al. 1998). Based on seismic monitoring, the hydrothermal system is supposed to be between 1 km above the sea level and 3 km below sea level (Boudon and Balcone-Boissard 2021).

## Climate and surface hydrology

The Montagne Pelée climate is a typical humid tropical climate exposed to trade winds, with a rainfall rate increasing with elevation and heavier rainfall on the east flank because

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Locations of boreholes, freshwater springs, warm springs and rivers on the geological map (Westercamp et al. 1990) modified thanks to the new volcanological evolution description by Boudon and Balcone-Boissard 2021</figcaption></figure>

of the orographic effects of the volcano. Average annual precipitation (Fig. 4a) ranges from 2,000 mm at sea level on the west coast and reaching 4,600–6,900 mm at the summit (Guiscafre et al. 1976; Vittecoq et al. 2010).

Around 20 rivers flow on the Montagne Pelée volcano flanks. Their specific low water flowrate ranges between 18 and 35 L/s/km² with a median value of 23 L/s/km² (min 10 L/s/km² and max 60 L/s/km²), the median value being similar on the two flanks. On the southwestern flank the notable difference is that the two rivers, the Sèche River and the Claire River, are mostly dry. The specificity of the Claire River is to be perennial in its upper part, between 800 and 480 m above mean sea level (amsl) elevation, and it fully infiltrates the pyroclastic deposits in an almost endorheic depression at 480 m amsl.

The large amounts of effective rainfall and the low river flows suggest that the amount of groundwater resource may be significant and should be stored in the pyroclastic deposits over thicknesses of a few tens to hundreds of meters. Nevertheless, knowledge regarding the extension and thickness of these deposits is not sufficiently detailed at the volcano scale. Furthermore, the complex geological evolution of the volcano and the impact of its asymmetry between the northeastern and southwestern flanks on groundwater circulation need to be characterized. These pyroclastic deposits also have interstitial porosity and unconfined aquifers, whereas aquifers over the rest of Martinique are mainly confined within fissured and fractured volcanic rocks (Vittecoq et al. 2015, 2019).
Thus, there is a necessity to more deeply characterize
the hydrodynamic characteristics of the pyroclastic depos-
its and the interaction between rivers and groundwater.
Given that the Montagne Pelée is an active volcano with
a hydrothermal system at depth, a better understanding
of the hydrogeological functioning is welcome to better
understand the role of groundwater in the eruptive pro-
cesses. There is then a clear interest to describe a hydro-
geological conceptual model at the volcano scale.

**Materials and methods**

The methodology is based on coupling resistivity data from
the heliborne electromagnetic survey conducted over the
volcano in March 2013 (Fig. 2) with hydrogeological and
geological field data and observations, especially thanks to
58 boreholes and 119 springs located on the volcano flanks,
and to an existing hydrological water budget model.

Borehole and spring databases

A database of the 58 boreholes (Appendix 1) located on the
volcano flanks has been assembled thanks to the BRGM
archives collected since the 1960s. The database contains
the main hydrogeological and geological available data:
location, elevation, borehole depth, water level, aquifer
geology, transmissivity, and hydraulic conductivity). Geo-
logical drilling descriptions were harmonized in order to
take into account the most recent geological knowledge
(Boudon and Balcone-Boissard 2021). Within this data-
base, transmissivity data are available for 21 boreholes.

The springs database (Appendix 2 and Appendix 3) contains 119 springs, compiled from field reconnaissance carried out between 2005 and 2021. The database contains location and spring elevations. Freshwater springs (96, Appendix 2) and thermal springs (23, Appendix 3) are also discerned. This inventory could not be exhaustive but is nevertheless representative of the main springs and of the different hydrogeological contexts.

Helicopter-borne time-domain electromagnetic
method

Martinique Island was covered by airborne electromagnetics (AEM) from February to March 2013 to address geological and hydrogeological purposes (Fig. 2). This survey, supervised by BRGM (French Geological Survey) and presented in Deparis et al. (2014), represented 4,233 line-kilometers for the whole island. The SkyTEM 304 system was chosen to image the shallow subsurface of this volcanic island

(Sørensen and Auken 2004). Developed by the HydroGeo-
physics Group of the University of Aarhus (Denmark), this
system operates in dual transmitter mode. The low moment
(2,826 A.m²: gates from 11 µs to 1 ms) provides information
on the very near surface (the first 50–100 m), while the high
moment (144.440 A.m²: gates from ~70 µs to 8.9 ms) gives
deeper information (~200 m).

The survey was mainly carried out along the N–S direc-
tion with a 400-m spacing and along the W–E direction with
a 4,000-m line spacing; the spacing was locally refined to
200 m over areas of interest such as the Montagne Pelée.
Along each flight line, AEM data were acquired every ~30
m, with an average ground clearance of 64 m due to the sharp
topography of the island. Figure 3 (and Figure SI_3 in the
ESM) shows the locations of the data over the studied area.

The AEM method allows one to image the conductivity/
resistivity contrasts of the subsurface (Ward and Hohmann
1988). Its depth of investigation (DOI), around 200 m for
this research, depends on the emitted magnetic moment, the
bandwidth used, the conductivity of the subsurface and the
signal/noise ratio (Spies 1989). AEM data were processed by
following the procedure described in Reninger et al. (2020),
which is based on the use of the singular value decomposition
(Reninger et al. 2011). The aim of the applied processing was
to keep as much resolution as possible (Reninger et al. 2020).

Finally, a manual editing procedure was performed,
mainly to remove remaining inductive/galvanic coupling
noises. In order to improve the coverage of the dataset, good
quality portions of ferry lines (helicopter routes between
each flight line departure and between the airport and each
day first flight line departure) were also considered dur-
ing the processing (Reninger et al. 2020). Data were then
inverted in “smoothed” one-dimensional (1D) models using
a quasi-3D spatially constrained inversion (SCI) algorithm
(Viezzoli et al. 2008). Each 1D model displays the resis-
tivity variations with depth. More 2D information can be
obtained, interpolating the resistivity layers falling into a
depth or elevation range. This step is generally repeated over
the entire range of investigation. The obtained slices can
then be merged to build a 3D resistivity model allowing
drawn sections in any direction and extraction of interfaces
and/or volumes.

Water budget calculation

A hydrological water budget model of Martinique has been
implemented (Vittecoq et al. 2007; Arnaud and Lanini 2014;
last update: Taïlamé and Lanini 2020) for water stakeholders
in order to assess water resources and water withdrawals.
This model allows daily to inter-annual water budget calcu-
lations at the watershed or subwatershed scale. This lumped
type hydrologic model requires meteorological data series:
rainfall rate and potential evapotranspiration, spatialized at
<figure><img src="image_8.png" /><figcaption>**Fig. 4** Water budget calculation results at the square kilometer scale: a Rainfall rate, b Potential evapotranspiration, c Effective rainfall, and d Infiltration ratio. Blue dots correspond to ground-water level monitoring points and green triangles to river flow-rate monitoring stations, whose data have been used for model calibration (Taïlamé and Lanini, 2020).</figcaption></figure>

NW, NE, SE, SW and the gridded summit zone correspond to the areas reported in section “Water budget”

a 1-km² grid, as well as time series related to water surface and groundwater withdrawals. A cultural map, soil-water-capacity map and runoff coefficient map at the kilometer scale are also needed as input data to calculate the effective rainfall, real evapotranspiration, infiltration/runoff ratio, flow rate at the outlet of rivers and the groundwater level in the underlying aquifers. In this model, Martinique Island has been divided into 1,206 cells, of which 227 concern Montagne Pelée and connected catchments. In each square kilometer cell, the exponential drainage of the soil reservoir is divided between runoff (fast surface and subsurface flow)

and infiltration (slower flow to the underlying first aquifer) considering two parameters—the runoff coefficient and the transit time. Runoff and infiltration flow from each square kilometer cell should then be aggregated at the catchment or subcatchment scale. The last model calibration was carried out (Taïlamé and Lanini 2020) by comparing the results of the model (flowrate at the outlet of rivers and groundwater level) with data from the river flowrate monitoring network and piezometric level monitoring network over the period 1991–2017. The calibration parameters at the square-kilometer-cell scale (soil reservoir drainage time and runoff
coefficient) and at the catchment scale (underground reservoir drainage time and porosity) are adjusted to optimize the model results.

In this research, Montagne Pelée has been divided into five areas: (1) a north-western zone corresponding to watersheds located on the first-stage geological unit, (2) a north-eastern zone and, (3) a southwestern zone, both corresponding to the opposite volcano flanks dominated by the 2nd and the 3rd stage geological units, (4) a southeastern zone corresponding to the Capot River watershed whose particularity is to be half on Montagne Pelée deposits and half on the Morne Jacob andesitic rocks (2.2–5.5 Ma) and, (5) a small zone corresponding to the summit zone (Fig. 4). Rainfall, real evapotranspiration, effective rainfall and infiltration have then been extracted for each area. A schematic view of the model is shown in Fig. SI_4 in the ESM.

# Results

## Resistivity ranges of pyroclastic aquifers and the main geological units

AEM results correlated with boreholes and springs data allow one to identify unsaturated zones, aquifers, and seawater intrusions as well as the main geological units. Eight characteristic cross-sections are presented in Figs. 5 and 6 and their locations are shown in Fig. 2 (and Figure SI 3 in the ESM). The main units highlighted are the following:

* Ranging between 50 and 1000 ohm.m, pyroclastic deposits are particularly well imaged. Within this resistivity range, three subsets can be distinguished:

- The higher resistivity values (>300 ohm.m) correspond to unsaturated and unweathered to slightly weathered pyroclastic deposits and are highlighted by “1” on Figs. 5 and 6.

- Resistivities between 70 and 300 ohm.m correspond to aquifers, confirmed by water level data in boreholes or the position of springs and are highlighted by “2” on Figs. 5 and 6. The higher aquifers values, ranging between 200–300 ohm.m, are observed for the “Pecoul” area and are highlighted by “3” on the C2 cross-sections (Fig. 5) and C6 (Fig. 6). Surficial aquifers of limited extent within 3rd stage pyroclastic deposits, with low flowrate springs, are highlighted by “4” in Figs. 5 and 6.

- Lower values, ranging between 50 and 70 ohm.m, are observed on the southwestern volcano flank between the Claire River and the Sèche River (highlighted by “5” on the C1 cross section on Fig. 5). This 50–100-m-thick layer is interpreted as mineralized (1,100–1,400 µS/cm) and warm (30–32 °C) groundwater.

* Below 50 ohm.m, two geological units, located below the pyroclastic deposits, are imaged:

- On the northwestern side of the volcano, with resistivity around 50 ohm.m, andesitic formations from the first stage (highlighted by “6” on Figs. 5 and 6),

- On the southeastern side, with resistivity below 30 ohm.m, andesitic formations from the Morne Jacob volcano (highlighted by “7” on Figs. 5 and 6).

* Resistivities lower than 15 ohm.m are observed around the coastal boundary of the volcano and correspond to seawater intrusions (highlighted by “8” on C1, C3 and C7 cross-sections on Figs. 5 and 6).

## Hydrogeological characteristics

Transmissivity values, calculated thanks to pumping tests conducted in 21 boreholes, varies between 9 $10^{-5}$ and $8 \times 10^{-2} \text{ m}^2/\text{s}$, with $10^{-3}$ and $6 \times 10^{-3} \text{ m}^2/\text{s}$ as the 1st and 3rd quartiles and with $2 \times 10^{-3} \text{ m}^2/\text{s}$ as the mean value. Hydraulic conductivity values, calculated by dividing transmissivity by the thickness of aquifer crossed by each borehole (Fig. 7), varies between $2 \times 10^{-6}$ and $3 \times 10^{-3} \text{ m/s}$, with $6 \times 10^{-5}$ and $3 \times 10^{-4} \text{ m/s}$ as the 1st and 3rd quartiles and $2 \times 10^{-4} \text{ m/s}$ as the mean value. The harmonization of geological descriptions and ages of aquifers crossed by each borehole, thanks to the new geological history classification (Boudon and Balcone-Boissard 2021), allows a clear correlation between hydraulic conductivities and ages (Fig. 7): the older the unconsolidated pyroclastic deposits are, the lower their hydraulic conductivity.

Only three storage coefficients have been reported. The higher value is $4 \times 10^{-1}$ in the younger pyroclastic aquifer on the southwestern flank (1,902 PDC deposits, 3rd stage, 2nd part). The lower is $3 \times 10^{-4}$ in clayey PDC, 4,400 years before on the northeastern flank (3rd stage, 2nd part) and should be the result of the rapid weathering of pyroclastic deposits due to the higher rainfall exposition of this flank. A middle value of $2 \times 10^{-3}$ in low silica subplinian–plinian deposits is also reported on the northeastern flank (3rd stage, 1st part).

Piezometric levels, measured in the boreholes intersecting water inflows, show two groups. The first one is composed of 16 boreholes located on the southwestern volcano flank (Fig. 8), between 8 and 80 m amsl, close to the sea with a maximum distance of 1,200 m from the coast and a hydraulic gradient of 2.2%. Water levels measured in these boreholes, generally below the level of the rivers, suggest that the rivers are infiltrating towards the aquifer. The second set is composed of 18 boreholes located on the northeastern volcano flank (Fig. 9), between 16 m and 614 m amsl, with water
<figure><img src="image_9.png" /><figcaption>**Fig. 5** Internal resistivity and hydrogeological structures along cross sections C1, C2, C3 and C4. The locations of the cross sections are shown on Fig. 2. Legend: (1) unsaturated and unweathered to slightly weathered pyroclastic deposits (>300 ohm.m), (2) aquifers (70–300 ohm.m), (3) ‘Pecoul’ aquifer (200–300 ohm.m), (4) limited extension</figcaption></figure>

superficial aquifers within 3rd stage pyroclastic deposits, (5) mineralized (1,100–1,400 µS/cm) and warm (30–32 °C) aquifer (50–70 ohm.m), (6) andesitic formations from the first stage (<50 ohm.m), (7) andesitic formations from Morne Jacob Volcano (<30 ohm.m), (8) seawater intrusions (<15 ohm.m)

level depth a few tens meters below ground level (mean: 41 m, min: 3 m, max: 97 m,) and a mean hydraulic gradient of 5.7%. Half of the boreholes have water level above the river level and the other half below, allowing distinguishing sectors with perched aquifers, sectors where groundwater flows towards rivers and sectors where rivers flow to the aquifer.

Freshwater springs are mainly located between 110 and 440 m (1st and 3rd quartile) with 280 m as the mean elevation. The two higher springs are at 900 and 950 m amsl. Most springs have low flowrates (1–10 m³/h) except one with higher flowrate (Morestin spring ≈ 200 m³/h). Almost all springs (90%) are on the northeastern volcano flank (Fig. 9). On this northeastern flank they can be classified into three equivalent distribution groups

(Fig. 9b). The first group corresponds to springs situated in the river's bed, intersecting the basal aquifer water table (depression spring type 1, Fig. 9). The second group corresponds to springs located in the middle part of the steep slopes of the river incisions, corresponding also to perched aquifers of limited extent, emerging thanks to paleosols or geological discontinuities (contact spring type 2, Fig. 9). The third group corresponds to springs located on the gentle slopes of the volcano flanks, corresponding to small-extension perched aquifers with subsurface flows and superficial water tables (depression spring type 3, Fig. 9). On the southwestern flank of the volcano (Fig. 8), the dozen springs are located in the upper part of watersheds (>4,000 m from the sea) and emerge from
<figure><img src="image_10.png" /><figcaption>**Fig. 6** Internal resistivity and hydrogeological structures along cross sections C5, C6, C7 and C8. The locations of the cross sections are shown on Fig. 2. Legend: (1) unsaturated and unweathered to slightly weathered pyroclastic deposits (>300 ohm.m), (2) aquifers (70–300 ohm.m), (3) ‘Pecoul’ aquifer (200–300 ohm.m), (4) limited extension</figcaption></figure>

superficial aquifers within 3rd stage pyroclastic deposits, (5) mineralized (1,100–1,400 µS/cm) and warm (30–32 °C) aquifer (50–70 ohm.m), (6) andesitic formations from the first stage (<50 ohm.m), (7) andesitic formations from Morne Jacob Volcano (<30 ohm.m), (8) seawater intrusions (<15 ohm.m)

perched aquifers. Three of these springs are located near the summit, the Samperre spring is located at the bottom of a large landslide (Peruzzetto et al. 2022), and the eight other springs are located east from the Roxelane River, outside of the flank collapses structures. Lastly, all (except one) of the thermal springs of the Montagne Pelée are also located on this southwestern flank.

## Water budget

Results of the water budget calculation (rainfall, evapotranspiration, effective rainfall, runoff and infiltration) for the five zones and the Montagne Pelée Volcano (Fig. 4) are given in Table 1. Water balance calculations show that, at the volcano scale, mean annual rainfall is 3,639 mm, mean annual potential evapotranspiration is 1,282 mm, mean annual effective rainfall is 2,456 mm and mean annual infiltration is 1,099 mm, corresponding to an annual groundwater recharge volume of 209 Mm³ (almost five times the annual drinking-water need for the island). Depending on the intensity of the rainy season, the interannual variability can be high with a maximum annual value reaching 6,861

mm for rainfall and 6,103 mm for effective rainfall. At the volcano scale, mean infiltration is 45% of effective rainfall and mean runoff is 55%. Depending on the watershed location, mean infiltration varies between 35 and 52% of the effective rainfall (and runoff between 48 and 65%).

Water budget results calculated on the northeastern flank with a repartition between 51% of infiltration and 49% of runoff are in agreement with Barat (1986) who calculated an infiltration rate of 50% (varying between 40 and 60% depending on the watershed). On the southwestern flank, with a repartition between 35% of infiltration and 65% of runoff, the results here are also in agreement with Barat (1986) who calculated an infiltration rate of 38 and 62% of runoff for the Roxelane River. The annual groundwater recharge volume is then greater on the southwestern flank (49 Mm³) than on the northeastern flank (42 Mm³).

## Seawater intrusions

The coastal fringe of the aquifers of this volcanic edifice is hardly affected by seawater intrusions. C1 and C3 cross-sections (Fig. 5) show very low resistivity values (<3 ohm.m)
<figure><img src="image_11.png" /><figcaption>**Fig. 7** Comparison between hydraulic conductivity and pyroclastic deposit ages for 21 boreholes in which pumping tests have been conducted. The older the pyroclastic deposits, the lower their hydraulic conductivity</figcaption></figure>

in the sea, corresponding to marine saline water and resistivity values <10–15 ohm.m in the continuity inland. These cross-sections are representative of the slope of the seawater intrusions (SWI) on most of the others around the volcanic edifice. With a slope of about –40%, SWI rapidly reaches

depths of up to 100 m within a few hundred meters from the coast, which is deeper than the DOI via the SkyTEM 304 system. This strong SWI slope is consistent, considering the strong groundwater recharge rate and the aquifer’s water-table slopes higher than 2%. Electrical conductivity

<figure><img src="image_12.png" /><figcaption>**Fig. 8** Topographic profiles of the Montagne Pelée southwestern flank rivers with piezometric levels in the boreholes and elevations of the fresh-water springs and thermal springs</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Fig. 9** a Topographic profiles of the Montagne Pelée northeastern flank rivers, with piezometric levels in the boreholes and elevations of the freshwater springs. Inset b shows the vertical arrangement of spring outflows</figcaption></figure>

measurements in the three closest boreholes to the coast are also in agreement (153, 280 and 418 µS/cm for distances between 150 and 300 m from the coast).

## Montagne Pelée hydrogeological conceptual model

Correlation between the 3D resistivity model and hydrogeological data from boreholes and springs allows for a description of the hydrogeological conceptual model of the Montagne Pelée volcano. Figure 10 presents a schematic view of the hydrogeological conceptual model, along a theoretical southwest-northeast axis. In addition to the hydrogeological structure of the first hundred meters, this model also shows, thanks to magneto telluric data (Gadalla et al. 2014), the cap rock, the supposed hydrothermal system and the upper crust.

## The upper perched aquifer

The upper aquifer is a perched aquifer (Fig. 10 and the zoom labelled A) located within the recent lava domes (3rd stage, 2nd part). These andesitic domes are deeply fissured and fractured, and the upper part of the volcano between the domes is a chaos of blocs allowing fast infiltration of the huge amount of effective rainfall on the top of the volcano (mean annual amount: 5,813 mm, corresponding to about 7.5–12 Mm³, depending on the surface area considered). Furthermore, an endorheic topography is observed between the bottom of the domes and the crater rim (elevation 1,200 m), also allowing a fast infiltration of rainfall. The area of this aquifer is difficult to determine precisely but is about 1–2 km². The thickness of the saturated zone is unknown, but it is supposed that the thickness is small, as spring flow-rates are low and water is expected to infiltrate at depth toward the deep hydrothermal system.

<table><caption>**Table 1** Water budget calculation results for rainfall, evapotranspiration, effective rainfall, runoff and infiltration for all five zones and the Montagne Pelée Volcano</caption><thead><tr><th></th><th>Area (km²)</th><th>Rainfall (mm)</th><th>PET (mm)</th><th>Effective rainfall (mm)</th><th>Infiltration (mm)</th><th>Runoff (mm)</th><th>% Infiltration</th><th>Infiltration (Mm³)</th></tr></thead><tbody><tr><td>NW</td><td>26.2</td><td>3,586</td><td>1,478</td><td>2,246</td><td>820</td><td>1,420</td><td>37%</td><td>22</td></tr><tr><td>NE</td><td>39.8</td><td>3,417</td><td>1,470</td><td>2,050</td><td>1,052</td><td>994</td><td>51%</td><td>42</td></tr><tr><td>SW</td><td>65.9</td><td>3,227</td><td>1,271</td><td>2,109</td><td>745</td><td>1,358</td><td>35%</td><td>49</td></tr><tr><td>SE</td><td>58.7</td><td>4,275</td><td>1,079</td><td>3,216</td><td>1,652</td><td>1,556</td><td>52%</td><td>97</td></tr><tr><td>Summit</td><td>1.3</td><td>6,571</td><td>754</td><td>5,813</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Volcano scale</td><td>190.6</td><td>3,639</td><td>1,282</td><td>2,456</td><td>1,099</td><td>1,351</td><td>45%</td><td>209</td></tr></tbody></table>

<figure><img src="image_14.png" /></figure>
The other aquifers can be categorized into three groups— the northeastern flank aquifers, the southwestern flank aquifer and the southeastern flank aquifers. The northwestern quarter of the volcano, characterized by first-stage primitive Montagne Pelée lavas flows, remains little known, as only one hiking trail exists near the coast, across the tropical forest, and wind conditions and rough topography did not allow the helicopter-borne survey.

## The northeastern flank aquifers

The northeastern flank is characterized by two types of aquifers (Fig. 10 and the zoom labelled B). In the first 10 m below ground level, small perched aquifers, within pyroclastic deposits of the third stage, give rise to low flowrate springs at various elevations (mainly between 100 and 400 m amsl). Beneath is a second aquifer (second-stage pyroclastic deposits and third stage, with the first part being pyroclastic deposits), corresponding to the basal aquifer flowing to the sea. This basal aquifer is also found at greater elevations (200–400 m amsl) at greater depths (200 m) depending on the morphology and the elevation of the first part lava and breccia located below.

## The southeastern flank aquifer

The geomorphology of the southeastern flank differs from the northeastern flank. At similar elevations to the northeastern flank, pyroclastic deposits, mainly from the third stage, only have a maximum of 200 m thicknesses. These deposits lie on the first stage and the Morne Jacob andesitic lavas which outcrop in some river beds. In the first 10 m, small perched aquifers give rise to low flowrate springs. Under these small aquifers is found the main aquifer of the flank, in which groundwater flows are controlled by (1) the topography of the first stage and Morne Jacob andesitic rocks, and (2) the incision of the Capot River and its tributaries within the pyroclastic deposits (C4 cross-section on Fig. 5). Due to this geomorphology, important groundwater drainage by the river gives the latter the higher flowrate of the island.

## The southwestern flank aquifers

The southwestern flank has several distinctive characteristics. The first particularity is its geomorphology, marked by horseshoe shape relief resulting from two flank collapses (127 ka and 36 ka), as exposed by Le Friant et al. (2003) and Boudon and Balcone-Boissard (2021). The 36-ka Rivière Sèche event is then particularly visible on the C6 cross-section (Fig. 6): several tens of meters of pyroclastics deposits are missing and the basis of the collapse should have been

controlled by the lava from the first stage. The basis of the 127-ka event is less obvious, but well visible on its southern rim (C6 cross-section, Fig. 6). It is reasonable to think that this event may not have happened all at once, but may be the result of successive collapses. Around 20 boreholes have been drilled inside these two structures and crossed mainly the third-stage deposits, enabling identification of a continuous basal aquifer. A deepening of the andesitic lavas (below 50 ohm.m) and consequently an increasing thickness of 200 m of saturated pyroclastic deposits is visible in the axis of the Pecoul sugarcane field (C2 and C6 cross-sections). The extension of this deepening does not correspond to the actual identified flank collapse horseshoe shape's structure and could then correspond to (1) an unknown masked structure or (2) the result of river erosion over thousands of years before being filled by pyroclastic deposits.

The second hydrogeological particularity of the southwestern flank is the quasi absence of small-perched aquifers at various elevation, as observed on the eastern flank, that can be linked to the effect of the successive flank collapses which have not preserved the structures that allow the existence of these small aquifers. The Morestin spring, the higher flowrate spring of the volcano and seven other springs (Fig. 8), located between 4,000 and 6,000 m from the sea and at an elevation between 300 and 600 m, are outside the flank collapse structures and should be associated with a perched aquifer. The third hydrogeological particularity of the southwestern flank is the presence of warm springs, testifying to the existence of a hydrothermal system at depth.

The fourth hydrogeological particularity is the complete infiltration of warm water from the Claire River around 500 m elevation and between 400 and 500 m elevation for the Chaude River, and the existence of mineralized groundwater downward flows particularly visible on the geophysical data (cf. C1 cross-section, Fig. 5). The hydrogeological functioning of this area is synthetized in Fig. 11, with the representation of the lateral extension of the 50 ohm.m layer and color scale according to the elevation of the top of the geophysical layer. Thermal water from the Claire and Chaude rivers infiltrates in highly permeable pyroclastic deposits (3rd part, 2nd stage: 1902–1905 and 1929–1932 PDCs deposits, also known as the “Coulee Rivière Blanche”) and flows toward the sea following the morphology of the flank before the 1902 eruptions. Thermal springs at shallow depth near the shore and the thermal well near the coast confirm this hypothesis.

## The aquitard bedrock

Finally, the hydrogeological aquitard bedrocks of the basal aquifers are the first-stage primitive Montagne Pelée (Mont Conil) lava flows of approximately two thirds in
<figure><img src="image_15.png" /><figcaption>**Fig. 10** a Hydrogeological conceptual model of the Montagne Pelée volcano as a SW-NE cross section, b enlargements for A (summit) and B (NE slope) and c hydrogeological log with geological units and associated stage of building, hydraulic conductivity and resistivity</figcaption></figure>
the northwest direction and the Morne Jacob andesitic lavas of the remaining third in the southeastern direction. Even if the hydraulic conductivity of these lavas is unknown, one can consider that the hydraulic conductivity contrast between the highly permeable pyroclastic deposits and the weathered lavas flows (considering their low electromagnetic resistivity) is enough to play the role of an aquitard bedrock.

Discussion

The hydrogeological conceptual model of the Montagne Pelée volcano has been defined thanks to the correlation of rarely available datasets on the same edifice—helicopter-borne electromagnetic data with a high resolution, a relatively high borehole density database, with water level and hydraulic conductivity data, and a spring database. The results, highlighting thick permeable pyroclastic deposits along the volcano flanks, should help researchers to better understand the hydrogeological functioning on the proximal (2–10 km) part of an andesitic stratovolcano, following the description of Vessell and Davies (1981), Bogie and Mackenzie (1998) and Selles et al. (2015). This model, nevertheless, could not be extrapolated to the rest of Martinique Island, as pyroclastic deposits only concern the Montagne Pelée. In the central part of Martinique Island, two other conceptual models have been developed at a watershed scale (Vittecoq et al. 2019) and an aquifer

scale (Vittecoq et al. 2015). At the Lesser Antilles scale, the classification of high-rise volcanic islands into three categories by Robins et al. (1990) should be updated with the detailed description given here.

The southwestern flank collapses have affected the geological evolution of the Montagne Pelée and have induced a structural difference between the northeastern and the southwestern flanks (Vincent et al. 1989; Le Friant et al. 2003; Boudon and Balcone-Boissard, 2021). These flank collapses are due to an eastern-western topographic asymmetry of the volcano with steeper on-land and underwater slopes in the western flank. This structural asymmetry is reflected in its hydrogeological functioning with perched aquifers and springs on one side and none on the other. Similar observations are expected on other volcanos as flank collapses due to the asymmetry of the volcano have been noticed on some others, whether in the West Indies (Soufriere of Saint Lucia and Soufriere of Saint Vincent (Boudon et al. 2007) or at the smaller-scale Soufriere of Guadeloupe (Rosas Barbajal et al. 2016) in USA (Mount Saint Helen, Lipman and Mullineaux 1981) or in Indonesia (Mount Merapi, Selles 2014). The Merapi volcano is for instance marked by a stable east flank with springs and perched aquifers, whereas the western flank is destabilized by recent eruptions and without springs (Selles 2014). These flank collapses induced a rejuvenation of the geomorphology by thick pyroclastic deposits and a weaker incision. Thus, associated aquifers remain hidden and do not allow springs outflow. As demonstrated

<figure><img src="image_16.png" /><figcaption>**Fig. 11** Interpreted pseudo-3D view of the mineralized (1,100–1,400 µS/cm) and warm (30–32 °C) aquifer (50–70 ohm.m), with representation of the lateral extension of the 50–70 ohm.m layer, with color scale according to the elevation of the top of the geophysical layer</figcaption></figure>
with the hydraulic conductivity dataset for Montagne Pelée, the younger the pyroclastic deposits are, the higher their hydraulic conductivity. The decrease of hydraulic conductivity as a function of time should be the result of weathering processes. If there is enough rainfall, significant aquifers are expected within flank collapses, as supported by self-potential measurements (Zlotnicki et al. 1998) showing the floor of the horseshoe-shaped structure. Conversely, the opposite flank should present a superposition of increasingly old deposits with decreasing hydraulic conductivity and lower storage coefficient. Paleo-soil, erosion and weathering processes (Rad et al. 2013) between each eruptive sequence also lead to horizontal discontinuities with lower permeability, allowing perched aquifers and associated springs. Lastly, helicopter-borne electromagnetic data probably show another horseshoe structure below the Pecoul sugarcane field masked by pyroclastic deposits filling.

Hydraulic conductivity values of the Montagne Pelée pyroclastic deposits, with 2 10⁻⁴ m/s as mean value, are on the same order of magnitude as Mount Mazama in the Oregon Cascade Range, USA (Aldous and Gurrieri 2011), with measured average hydraulic conductivity of plinian pumice layers at 1.4 10⁻⁴ m/s. These values also match with fine sand or gravel hydraulic conductivity, according to Domenico and Schwartz (1990), facilitating the similarities in order to better explain groundwater flows to the public and stakeholders.

Groundwater also plays a key role in phreatic eruptions. The previous eruption cycles at Montagne Pelée (1902–1905 and 1929–1932) started with phreatic eruptions. In 1792 and 1851, phreatic eruptions were also reported, without having been followed by magmatic eruptions. Each time, fumaroles activity is noticed before phreatic eruptions with a rapid decrease after each eruption cycle (Barat 1986). This fumaroles activity should be evidence that the upper-perched aquifer is alternatively heated by conduction from deeper magma and by magmatic gas when the magma rises (going until boiling then vaporized when phreatic eruption occurred) and acting as a buffer against rising gas when the eruption stops. The rapid decrease of fumaroles activity is also evidence that the aquifer rapidly “refills” thanks to the high quantity of effective rainfall thus cooling and diluting the effect of the magmatic gas.

It is assumed that the recharge of the hydrothermal system, located between 1.5 and 3 km depth, is via deep vertical groundwater flow restricted to the central crater system (Traineau et al. 1989; Zlotnicky 1998; Gadalia

et al. 2014). The results here show that the upper perched aquifer gathers rainfall from the endorheic system within the crater. The low real evapotranspiration rate, together with a weak runoff, favours a high infiltration rate. The springs' low flowrates from this upper aquifer and the interpreted structure of the first stage of edification confirmed this hypothesis, but the “connection” with the hydrothermal system remains unclear as the upper magma pathway is supposed to be a maximum of a few tens of meters (corresponding to the diameter of the lava spine) and sealed. Furthermore, the cap rock identified by magneto telluric data (Gadalia et al. 2014) pressurized the hydrothermal system, with few leaks (along potential fractures) leading to the existing thermal springs and past fumaroles. Deep infiltrated water should then refill the hydrothermal system only if the pressure of the water height in the upper aquifer and the supposed fractures near the plumbing system is higher than the pressure of the hydrothermal system. Better knowledge of the hydrothermal system extension and pressure is needed to more deeply characterize the relationship with the shallow aquifer.

## Conclusion

The hydrogeological conceptual model of the Montagne Pelée Volcano has been defined thanks to the correlation of rare data on the same volcanic massif—helicopter-borne electromagnetic data with a high resolution, a relatively high density of boreholes with water level and hydraulic conductivity data, and a spring database. The study has also demonstrated that the older the pyroclastic deposits, the lower their hydraulic conductivity. The on-land and underwater asymmetric topography of the volcano has induced flank collapses leading to a structural difference between the northeastern and southwestern flanks, as reflected in its hydrogeological functioning. Finally, this paper discusses relationships between past phreatic eruptions and the hydrothermal system.

The way forward now is to describe the geochemical characteristics of these multiple aquifers, to more deeply characterize the structure of the volcano below the DOI of the geophysical survey (≈200 m), in order to build a 3D geological and hydrogeological model at the volcano scale and to improve understanding of groundwater flows at the scale of each flank.

<figure><img src="image_17.png" /></figure>
Appendix 1

<table><caption>Table 2 Borehole database (WGS84 coordinate system)</caption><thead><tr><th>National number</th><th>City</th><th>Boreholes name or locality</th><th>X (m)</th><th>Y (m)</th><th>Borehole elevation (m)</th><th>Borehole depth (m)</th><th>Water level elevation (in m amsl)</th><th>T (m<sup>2</sup>/s)</th><th>K (m/s)</th><th>Aquifer geology</th></tr></thead><tbody><tr><td>1165Z0001</td><td>Grand Rivière</td><td>F1</td><td>696,316</td><td>1,644,519</td><td>30</td><td>36</td><td>26</td><td>1.0E-03</td><td>5E-05</td><td>3rd ST. 1st P.</td></tr><tr><td>1165Z0012</td><td>Grand Rivière</td><td>F4</td><td>696,276</td><td>1,644,549</td><td>30</td><td>15</td><td>26</td><td>1.0E-03</td><td>1E-04</td><td>3rd ST. 1st P.</td></tr><tr><td>1165Z0013</td><td>Grand Rivière</td><td>F5</td><td>696,285</td><td>1,644,546</td><td>30</td><td>15</td><td>26</td><td>2.0E-03</td><td>2E-04</td><td>3rd ST. 1st P.</td></tr><tr><td>1165Z0014</td><td>Macouba</td><td>Beauséjour BR1</td><td>696,717</td><td>1,644,902</td><td>136</td><td>160</td><td>39</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0004</td><td>Basse Pointe</td><td>Vivé P3</td><td>705,586</td><td>1,642,059</td><td>16.5</td><td>30</td><td>14.5</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0026</td><td>Basse Pointe</td><td>Chalvet</td><td>705,086</td><td>1,643,059</td><td>30</td><td>29</td><td>12</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0027</td><td>Grand Rivière</td><td>F2</td><td>696,272</td><td>1,644,570</td><td>30</td><td>40</td><td>26</td><td>2.0E-03</td><td>1E-04</td><td>3rd ST. 1st P.</td></tr><tr><td>1166Z0028</td><td>Grand Rivière</td><td>F3</td><td>696,265</td><td>1,644,576</td><td>30</td><td>15</td><td>26</td><td>7.0E-04</td><td>7E-05</td><td>3rd ST. 1st P.</td></tr><tr><td>1166Z0055</td><td>Basse Pointe</td><td>Leyritz</td><td>702,503</td><td>1,642,166</td><td>165</td><td>28</td><td>149</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0056</td><td>Basse Pointe</td><td>Leyritz</td><td>702,652</td><td>1,641,980</td><td>163</td><td>41</td><td>134</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0057</td><td>Basse Pointe</td><td>Leyritz</td><td>701,871</td><td>1,641,560</td><td>230</td><td>50</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1166Z0059</td><td>Macouba</td><td>Loucou CHE1</td><td>698,380</td><td>1,643,494</td><td>259</td><td>160</td><td>147</td><td>2.4E-04</td><td>5E-06</td><td>2nd ST.</td></tr><tr><td>1166Z0060</td><td>Macouba</td><td>Perpigna PE2</td><td>698,450</td><td>1,644,928</td><td>125</td><td>150</td><td>31</td><td>1.2E-03</td><td>2E-05</td><td>2nd ST.</td></tr><tr><td>1166Z0061</td><td>Macouba</td><td>Potiche PR1</td><td>697,705</td><td>1,645,101</td><td>120</td><td>145</td><td>30</td><td>1.9E-04</td><td>4E-06</td><td>2nd ST.</td></tr><tr><td>1167Z0008</td><td>Saint Pierre</td><td>Fond Canonville</td><td>693,856</td><td>1,635,229</td><td>63</td><td>50</td><td>28</td><td>6.0E-03</td><td>3E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0021</td><td>Saint Pierre</td><td>Fond Corre 1</td><td>695,536</td><td>1,632,379</td><td>29</td><td>50</td><td>5</td><td>3.1E-02</td><td>1E-03</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0022</td><td>Saint Pierre</td><td>Fond Corre 2</td><td>695,726</td><td>1,632,589</td><td>55</td><td>60</td><td>13</td><td>1.1E-02</td><td>6E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0023</td><td>Saint Pierre</td><td>Coulée Rivière Blanche</td><td>694,666</td><td>1,633,319</td><td>24.23</td><td>60</td><td>0.23</td><td>4.4E-02</td><td>1E-03</td><td>3rd ST. 2nd P.</td></tr><tr><td>BSS003QEGW</td><td>Saint Pierre</td><td>Coulée Rivière Blanche CB1</td><td>694,574</td><td>1,633,484</td><td>57.7</td><td>82</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0024</td><td>Precheur</td><td>Ptészomètre BRGM</td><td>691,446</td><td>1,637,549</td><td>49</td><td>50</td><td>13</td><td>3.0E-03</td><td>2E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0027</td><td>Saint Pierre</td><td>Fond Canonville</td><td>693,436</td><td>1,634,839</td><td>22</td><td>20</td><td>7</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0028</td><td>Saint Pierre</td><td>Fond Canonville</td><td>693,636</td><td>1,635,049</td><td>48</td><td>40</td><td>21</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0029</td><td>Saint Pierre</td><td>Amont CDST</td><td>696,266</td><td>1,632,119</td><td>57</td><td>63</td><td>15</td><td>1.5E-03</td><td>7E-05</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0031</td><td>Saint Pierre</td><td>Pecoul</td><td>696,881</td><td>1,631,966</td><td>73</td><td>73</td><td>30</td><td>3.0E-02</td><td>1E-03</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0035</td><td>Saint Pierre</td><td>Tombeau des Caraîbes SC4</td><td>694,196</td><td>1,634,499</td><td>55</td><td>20</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0036</td><td>Saint Pierre</td><td>Tombeau des Caraîbes SC5</td><td>693,926</td><td>1,634,279</td><td>100</td><td>20</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0037</td><td>Saint Pierre</td><td>CET Fond Canonville</td><td>693,133</td><td>1,634,897</td><td>8.03</td><td>20.5</td><td>0.03</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0038</td><td>Saint Pierre</td><td>CET Fond Canonville</td><td>693,221</td><td>1,634,855</td><td>12.2</td><td>20.5</td><td>0.2</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0040</td><td>Saint Pierre</td><td>Pecoul FR3</td><td>696,801</td><td>1,632,144</td><td>80.6</td><td>77.5</td><td>28.6</td><td>7.7E-02</td><td>3E-03</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0041</td><td>Saint Pierre</td><td>Pecoul</td><td>696,824</td><td>1,632,078</td><td>79.4</td><td>66</td><td>31.4</td><td>4.7E-03</td><td>3E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0042</td><td>Saint Pierre</td><td>Pecoul FR1 bis</td><td>696,857</td><td>1,632,023</td><td>78.1</td><td>71.9</td><td>30.1</td><td>4.7E-03</td><td>2E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1167Z0045</td><td>Saint Pierre</td><td>Ptészomètre CDST</td><td>695,886</td><td>1,631,929</td><td>30.27</td><td>38.5</td><td>5.27</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1167Z0058</td><td>Saint Pierre</td><td>Pecoul FR7</td><td>696,885</td><td>1,631,971</td><td>77</td><td>72</td><td>28</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>
<figure><img src="image_18.png" /><figcaption>Table 2 (continued)</figcaption></figure>

<table><thead><tr><th>National number</th><th>City</th><th>Boreholes name or locality</th><th>X (m)</th><th>Y (m)</th><th>Borehole elevation (m)</th><th>Borehole depth (m)</th><th>Water level elevation (in m amsl)</th><th>T (m²/s)</th><th>K (m/s)</th><th>Aquifer geology</th></tr></thead><tbody><tr><td>1168ZZ0007</td><td>Ajoupa Bouillon</td><td>Saut Babin S1bis</td><td>704,056</td><td>1,638,469</td><td>158.4</td><td>35</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0008</td><td>Ajoupa Bouillon</td><td>Fond Labour S4</td><td>704,386</td><td>1,638,009</td><td>161.6</td><td>50</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0012</td><td>Ajoupa Bouillon</td><td>Ravine Bœuf S8</td><td>704,916</td><td>1,638,079</td><td>153.9</td><td>30</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0013</td><td>Ajoupa Bouillon</td><td>L'école S9</td><td>705,366</td><td>1,639,269</td><td>169.6</td><td>45</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0015</td><td>Ajoupa Bouillon</td><td>Dufailly S1</td><td>703,906</td><td>1,639,009</td><td>174.4</td><td>50</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0016</td><td>Ajoupa Bouillon</td><td>Vivies S2</td><td>703,966</td><td>1,639,829</td><td>178.5</td><td>50</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0017</td><td>Ajoupa Bouillon</td><td>Les ombrages S3</td><td>704,126</td><td>1,640,089</td><td>145.8</td><td>20</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0018</td><td>Ajoupa Bouillon</td><td>Trou congo S4</td><td>704,631</td><td>1,640,499</td><td>66.3</td><td>10</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0042</td><td>Ajoupa Bouillon</td><td>Grande Savane</td><td>703,326</td><td>1,638,649</td><td>241</td><td>45</td><td>212</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0043</td><td>Basse Pointe</td><td>Morne Balai</td><td>700,846</td><td>1,639,949</td><td>420</td><td>61</td><td>403</td><td>9.0E-05</td><td>2E-06</td><td>3rd ST. 1st P.</td></tr><tr><td>1168ZZ0044</td><td>Ajoupa Bouillon</td><td>Grande Savane</td><td>703,319</td><td>1,638,599</td><td>241</td><td>46.5</td><td>211</td><td>3.0E-03</td><td>2E-04</td><td>3rd ST. 2nd P.</td></tr><tr><td>1168ZZ0045</td><td>Ajoupa Bouillon</td><td>Grande Savane</td><td>703,326</td><td>1,638,649</td><td>230</td><td>45</td><td>201</td><td>1.0E-03</td><td>6E-05</td><td>3rd ST. 2nd P.</td></tr><tr><td>1168ZZ0048</td><td>Basse Pointe</td><td>Démare</td><td>701,586</td><td>1,640,419</td><td>290</td><td>41</td><td>273</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0054</td><td>Basse Pointe</td><td>Chez Lélène</td><td>704,241</td><td>1,640,671</td><td>84</td><td>28</td><td>81</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0061</td><td>Ajoupa Bouillon</td><td>Savane Pascal</td><td>701,340</td><td>1,638,360</td><td>426</td><td>70</td><td>376</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0062</td><td>Ajoupa Bouillon</td><td>Savane Pascal</td><td>701,076</td><td>1,638,480</td><td>440</td><td>70</td><td>386</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0069</td><td>Ajoupa Bouillon</td><td>Fr2</td><td>700,165</td><td>1,637,602</td><td>614</td><td>74</td><td>562</td><td>-</td><td>-</td><td>-</td></tr><tr><td>1168ZZ0070</td><td>Ajoupa Bouillon</td><td>Fr3</td><td>700,349</td><td>1,637,637</td><td>592</td><td>70</td><td>535</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003GXLM</td><td>Girand Rivière</td><td>Beauséjour E1SD1b</td><td>697,168</td><td>1,644,504</td><td>177.69</td><td>25</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003GYRU</td><td>Grand Rivière</td><td>Beauséjour E2SD4</td><td>697,201</td><td>1,644,302</td><td>200.82</td><td>30</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003GYVQ</td><td>Grand Rivière</td><td>Beauséjour E3SP1</td><td>697,021</td><td>1,643,996</td><td>239.25</td><td>35</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003GZUG</td><td>Grand Rivière</td><td>Beauséjour E4SD10</td><td>696,763</td><td>1,643,861</td><td>240.89</td><td>25</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003HALY</td><td>Grand Rivière</td><td>Beauséjour E5SD13</td><td>697,019</td><td>1,643,667</td><td>273.48</td><td>25</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003HBXQ</td><td>Grand Rivière</td><td>Beauséjour E6SP2</td><td>697,091</td><td>1,643,355</td><td>301.48</td><td>35</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>BSS003HBYK</td><td>Grand Rivière</td><td>Beauséjour E7SP3</td><td>697,091</td><td>1,643,104</td><td>323.83</td><td>41</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>
Appendix 2

<table><caption>Table 3 Freshwater springs database (WGS84 coordinate system)</caption><thead><tr><th>National No.</th><th>City name</th><th>X (m)</th><th>Y (m)</th><th>Z (m)</th></tr></thead><tbody><tr><td>BSS002NMEU</td><td>GRAND'RIVIERE</td><td>696,868</td><td>1,644,875</td><td>117</td></tr><tr><td>BSS002NMEV</td><td>GRAND'RIVIERE</td><td>696,425</td><td>1,644,771</td><td>95</td></tr><tr><td>BSS002NMEW</td><td>MACOUBA</td><td>697,327</td><td>1,643,256</td><td>308</td></tr><tr><td>BSS002NMEX</td><td>MACOUBA</td><td>697,410</td><td>1,644,310</td><td>121</td></tr><tr><td>BSS002NMEY</td><td>MACOUBA</td><td>697,391</td><td>1,644,243</td><td>130</td></tr><tr><td>BSS002NMFX</td><td>BASSE-POINTE</td><td>700,006</td><td>1,641,830</td><td>323</td></tr><tr><td>BSS002NMFY</td><td>BASSE-POINTE</td><td>702,186</td><td>1,642,920</td><td>133</td></tr><tr><td>BSS002NMFZ</td><td>MACOUBA</td><td>699,286</td><td>1,641,170</td><td>468</td></tr><tr><td>BSS002NMGA</td><td>MACOUBA</td><td>699,900</td><td>1,643,075</td><td>182</td></tr><tr><td>BSS002NMGB</td><td>MACOUBA</td><td>699,886</td><td>1,643,020</td><td>183</td></tr><tr><td>BSS002NMGC</td><td>MACOUBA</td><td>699,206</td><td>1,645,579</td><td>29</td></tr><tr><td>BSS002NMGD</td><td>BASSE-POINTE</td><td>702,599</td><td>1,644,706</td><td>8</td></tr><tr><td>BSS002NMGE</td><td>MACOUBA</td><td>700,499</td><td>1,644,830</td><td>70</td></tr><tr><td>BSS002NMGG</td><td>BASSE-POINTE</td><td>700,508</td><td>1,642,183</td><td>252</td></tr><tr><td>BSS002NMGH</td><td>BASSE-POINTE</td><td>703,606</td><td>1,643,699</td><td>56</td></tr><tr><td>BSS002NMGJ</td><td>BASSE-POINTE</td><td>705,842</td><td>1,642,545</td><td>13</td></tr><tr><td>BSS002NMGK</td><td>BASSE-POINTE</td><td>702,274</td><td>1,641,452</td><td>181</td></tr><tr><td>BSS002NMGL</td><td>BASSE-POINTE</td><td>701,866</td><td>1,642,373</td><td>154</td></tr><tr><td>BSS002NMGM</td><td>BASSE-POINTE</td><td>703,905</td><td>1,641,285</td><td>124</td></tr><tr><td>BSS002NMGN</td><td>BASSE-POINTE</td><td>702,592</td><td>1,641,441</td><td>180</td></tr><tr><td>BSS002NMGR</td><td>BASSE-POINTE</td><td>703,835</td><td>1,642,325</td><td>79</td></tr><tr><td>BSS002NMGS</td><td>BASSE-POINTE</td><td>703,643</td><td>1,641,552</td><td>137</td></tr><tr><td>BSS002NMGT</td><td>LORRAIN</td><td>705,962</td><td>1,641,365</td><td>37</td></tr><tr><td>BSS002NMGV</td><td>MACOUBA</td><td>698,237</td><td>1,642,492</td><td>368</td></tr><tr><td>BSS002NMGW</td><td>MACOUBA</td><td>699,467</td><td>1,645,187</td><td>87</td></tr><tr><td>BSS002NMGX</td><td>MACOUBA</td><td>700,287</td><td>1,644,296</td><td>108</td></tr><tr><td>BSS002NMGY</td><td>MACOUBA</td><td>699,345</td><td>1,641,899</td><td>325</td></tr><tr><td>BSS002NMGZ</td><td>MACOUBA</td><td>699,013</td><td>1,642,520</td><td>326</td></tr><tr><td>BSS002NMHA</td><td>MACOUBA</td><td>698,849</td><td>1,643,932</td><td>115</td></tr><tr><td>BSS002NMHB</td><td>MACOUBA</td><td>698,633</td><td>1,643,539</td><td>203</td></tr><tr><td>BSS002NMHC</td><td>BASSE-POINTE</td><td>705,368</td><td>1,642,936</td><td>14</td></tr><tr><td>BSS002NMHD</td><td>BASSE-POINTE</td><td>700,282</td><td>1,641,938</td><td>294</td></tr><tr><td>BSS002NMHE</td><td>BASSE-POINTE</td><td>702,252</td><td>1,642,764</td><td>127</td></tr><tr><td>BSS002NMHF</td><td>MACOUBA</td><td>698,849</td><td>1,643,650</td><td>163</td></tr><tr><td>BSS002NMHG</td><td>MACOUBA</td><td>698,830</td><td>1,642,811</td><td>206</td></tr><tr><td>BSS002NMHH</td><td>MACOUBA</td><td>699,461</td><td>1,645,374</td><td>75</td></tr><tr><td>BSS002NMHJ</td><td>MACOUBA</td><td>700,673</td><td>1,644,466</td><td>108</td></tr><tr><td>BSS002NMNJ</td><td>AJOUPA-BOUILLON</td><td>700,457</td><td>1,637,684</td><td>577</td></tr><tr><td>BSS002NMNK</td><td>AJOUPA-BOUILLON</td><td>701,218</td><td>1,637,271</td><td>471</td></tr><tr><td>BSS002NMNM</td><td>BASSE-POINTE</td><td>700,186</td><td>1,639,570</td><td>516</td></tr><tr><td>BSS002NMNN</td><td>SAINT-PIERRE</td><td>699,613</td><td>1,635,940</td><td>509</td></tr><tr><td>BSS002NMNP</td><td>AJOUPA-BOUILLON</td><td>698,353</td><td>1,638,131</td><td>909</td></tr><tr><td>BSS002NMNQ</td><td>SAINT-PIERRE</td><td>698,132</td><td>1,637,694</td><td>942</td></tr><tr><td>BSS002NMNR</td><td>MORNE-ROUGE</td><td>699,413</td><td>1,633,735</td><td>319</td></tr><tr><td>BSS002NMNT</td><td>MORNE-ROUGE</td><td>700,518</td><td>1,636,384</td><td>516</td></tr><tr><td>BSS002NMNU</td><td>BASSE-POINTE</td><td>701,009</td><td>1,640,118</td><td>380</td></tr><tr><td>BSS002NMNW</td><td>MORNE-ROUGE</td><td>700,648</td><td>1,635,733</td><td>509</td></tr></tbody></table>

<table><caption>Table 3 (continued)</caption><thead><tr><th>National No.</th><th>City name</th><th>X (m)</th><th>Y (m)</th><th>Z (m)</th></tr></thead><tbody><tr><td>BSS002NMNX</td><td>AJOUPA-BOUILLON</td><td>701,797</td><td>1,637,786</td><td>338</td></tr><tr><td>BSS002NMNY</td><td>AJOUPA-BOUILLON</td><td>704,116</td><td>1,640,286</td><td>112</td></tr><tr><td>BSS002NMPH</td><td>BASSE-POINTE</td><td>701,579</td><td>1,640,429</td><td>283</td></tr><tr><td>BSS002NMPJ</td><td>AJOUPA-BOUILLON</td><td>701,871</td><td>1,637,816</td><td>357</td></tr><tr><td>BSS002NMPK</td><td>AJOUPA-BOUILLON</td><td>701,888</td><td>1,637,710</td><td>361</td></tr><tr><td>BSS002NMPL</td><td>AJOUPA-BOUILLON</td><td>703,123</td><td>1,638,689</td><td>243</td></tr><tr><td>BSS002NMPM</td><td>AJOUPA-BOUILLON</td><td>701,899</td><td>1,638,287</td><td>353</td></tr><tr><td>BSS002NMPN</td><td>BASSE-POINTE</td><td>701,373</td><td>1,640,649</td><td>284</td></tr><tr><td>BSS002NMPP</td><td>BASSE-POINTE</td><td>701,186</td><td>1,640,475</td><td>306</td></tr><tr><td>BSS002NMPQ</td><td>BASSE-POINTE</td><td>700,871</td><td>1,639,950</td><td>413</td></tr><tr><td>BSS002NMPR</td><td>BASSE-POINTE</td><td>700,145</td><td>1,639,714</td><td>463</td></tr><tr><td>BSS002NMPW</td><td>LORRAIN</td><td>705,210</td><td>1,640,645</td><td>62</td></tr><tr><td>BSS002NMPY</td><td>LORRAIN</td><td>706,501</td><td>1,640,487</td><td>89</td></tr><tr><td>BSS002NMPZ</td><td>LORRAIN</td><td>706,157</td><td>1,640,525</td><td>98</td></tr><tr><td>BSS002NMQA</td><td>LORRAIN</td><td>704,188</td><td>1,637,731</td><td>169</td></tr><tr><td>BSS002NMQB</td><td>LORRAIN</td><td>703,737</td><td>1,637,752</td><td>255</td></tr><tr><td>BSS002NMQC</td><td>LORRAIN</td><td>704,911</td><td>1,639,548</td><td>158</td></tr><tr><td>BSS002NMQD</td><td>LORRAIN</td><td>704,943</td><td>1,638,054</td><td>168</td></tr><tr><td>BSS002NMQE</td><td>LORRAIN</td><td>705,228</td><td>1,637,875</td><td>193</td></tr><tr><td>BSS002NMQF</td><td>LORRAIN</td><td>704,345</td><td>1,637,461</td><td>182</td></tr><tr><td>BSS002NMQN</td><td>MORNE-ROUGE</td><td>700,463</td><td>1,635,844</td><td>521</td></tr><tr><td>BSS002NMQP</td><td>MORNE-ROUGE</td><td>700,078</td><td>1,633,848</td><td>363</td></tr><tr><td>BSS002NMQU</td><td>MORNE-ROUGE</td><td>700,449</td><td>1,634,790</td><td>463</td></tr><tr><td>BSS002NMQV</td><td>MORNE-ROUGE</td><td>700,094</td><td>1,633,080</td><td>340</td></tr><tr><td>BSS002NMQW</td><td>MORNE-ROUGE</td><td>700,532</td><td>1,636,359</td><td>526</td></tr><tr><td>BSS002NMQX</td><td>MORNE-ROUGE</td><td>700,220</td><td>1,636,200</td><td>539</td></tr><tr><td>BSS002NMQY</td><td>MORNE-ROUGE</td><td>699,928</td><td>1,636,314</td><td>600</td></tr><tr><td>BSS002NMQZ</td><td>MORNE-ROUGE</td><td>700,705</td><td>1,636,525</td><td>524</td></tr><tr><td>BSS002NMRE</td><td>AJOUPA-BOUILLON</td><td>701,820</td><td>1,637,709</td><td>366</td></tr><tr><td>BSS002NMRF</td><td>BASSE-POINTE</td><td>704,173</td><td>1,640,755</td><td>98</td></tr><tr><td>BSS002NMRG</td><td>AJOUPA-BOUILLON</td><td>700,961</td><td>1,638,730</td><td>460</td></tr><tr><td>BSS004HMLQ</td><td>AJOUPA-BOUILLON</td><td>699,228</td><td>1,637,997</td><td>674</td></tr><tr><td>BSS004HMLR</td><td>AJOUPA-BOUILLON</td><td>699,283</td><td>1,638,068</td><td>668</td></tr><tr><td>BSS004HMLS</td><td>SAINT-PIERRE</td><td>696,413</td><td>1,637,834</td><td>889</td></tr><tr><td>BSS004HMLT</td><td>MORNE-ROUGE</td><td>699,864</td><td>1,635,313</td><td>445</td></tr><tr><td>BSS004HMLU</td><td>MORNE-ROUGE</td><td>700,391</td><td>1,634,530</td><td>424</td></tr><tr><td>BSS004HMLV</td><td>PRECHEUR</td><td>695,979</td><td>1,638,964</td><td>600</td></tr><tr><td>BSS004HMLW</td><td>PRECHEUR</td><td>696,536</td><td>1,638,484</td><td>895</td></tr><tr><td>BSS004HMLX</td><td>MORNE-ROUGE</td><td>700,551</td><td>1,636,441</td><td>519</td></tr><tr><td>BSS004HMLZ</td><td>AJOUPA-BOUILLON</td><td>701,039</td><td>1,637,100</td><td>419</td></tr><tr><td>BSS004HMMA</td><td>AJOUPA-BOUILLON</td><td>701,052</td><td>1,637,111</td><td>419</td></tr><tr><td>BSS004HMMB</td><td>AJOUPA-BOUILLON</td><td>701,916</td><td>1,638,324</td><td>351</td></tr><tr><td>BSS004HMMC</td><td>AJOUPA-BOUILLON</td><td>704,090</td><td>1,640,130</td><td>132</td></tr><tr><td>BSS004HMMD</td><td>SAINT-PIERRE</td><td>699,182</td><td>1,636,330</td><td>563</td></tr><tr><td>BSS004HMME</td><td>BASSE-POINTE</td><td>699,089</td><td>1,639,074</td><td>677</td></tr><tr><td>BSS004HMMF</td><td>MACOUBA</td><td>697,749</td><td>1,639,784</td><td>689</td></tr><tr><td>BSS004HMMG</td><td>BASSE-POINTE</td><td>698,511</td><td>1,639,771</td><td>723</td></tr><tr><td>BSS002NMGQ</td><td>BASSE-POINTE</td><td>703,986</td><td>1,642,363</td><td>92</td></tr><tr><td>BSS002NMGP</td><td>BASSE-POINTE</td><td>703,950</td><td>1,642,446</td><td>84</td></tr></tbody></table>
# Appendix 3

<table><caption>**Table 4** Thermal springs database (WGS84 coordinate system)</caption><thead><tr><th>X (m)</th><th>Y (m)</th><th>Z (m)</th><th>Group Name</th></tr></thead><tbody><tr><td>696,319</td><td>1,637,000</td><td>541</td><td>Chaude River springs</td></tr><tr><td>696,399</td><td>1,637,060</td><td>583</td><td>Chaude River springs</td></tr><tr><td>696,462</td><td>1,637,050</td><td>604</td><td>Chaude River springs</td></tr><tr><td>696,511</td><td>1,637,050</td><td>604</td><td>Chaude River springs</td></tr><tr><td>696,308</td><td>1,637,660</td><td>768</td><td>Claire River springs</td></tr><tr><td>696,284</td><td>1,637,660</td><td>768</td><td>Claire River springs</td></tr><tr><td>696,306</td><td>1,637,680</td><td>784</td><td>Claire River springs</td></tr><tr><td>696,154</td><td>1,637,530</td><td>727</td><td>Claire River springs</td></tr><tr><td>696,191</td><td>1,637,560</td><td>730</td><td>Claire River springs</td></tr><tr><td>696,001</td><td>1,637,320</td><td>648</td><td>Claire River springs</td></tr><tr><td>695,954</td><td>1,637,130</td><td>576</td><td>Claire River springs</td></tr><tr><td>694,163</td><td>1,636,700</td><td>233</td><td>Mitan and Picodo springs</td></tr><tr><td>694,027</td><td>1,636,650</td><td>219</td><td>Mitan and Picodo springs</td></tr><tr><td>693,663</td><td>1,636,750</td><td>222</td><td>Mitan and Picodo springs</td></tr><tr><td>693,523</td><td>1,636,510</td><td>158</td><td>Mitan and Picodo springs</td></tr><tr><td>694,021</td><td>1,633,250</td><td>0</td><td>Underwater springs</td></tr><tr><td>693,717</td><td>1,633,800</td><td>0</td><td>Underwater springs</td></tr><tr><td>693,832</td><td>1,633,470</td><td>0</td><td>Underwater springs</td></tr><tr><td>693,766</td><td>1,633,680</td><td>0</td><td>Underwater springs</td></tr><tr><td>693,644</td><td>1,633,910</td><td>0</td><td>Underwater springs</td></tr><tr><td>695,233</td><td>1,632,340</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>695,283</td><td>1,632,270</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>695,376</td><td>1,632,140</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>695,386</td><td>1,632,080</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>695,408</td><td>1,632,020</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>695,433</td><td>1,631,990</td><td>0</td><td>Underwater degassing zone</td></tr><tr><td>697,140</td><td>1,639,530</td><td>865</td><td>Grande Rivière spring</td></tr><tr><td>698,039</td><td>1,637,220</td><td>790</td><td>Des Pères River springs</td></tr><tr><td>698,132</td><td>1,637,180</td><td>753</td><td>Des Pères River springs</td></tr></tbody></table>

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s10040-023-02642-5.

**Declarations** Conflict of interest On behalf of all authors, the corresponding author states that there is no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

# References

Aldous A, Gurrieri J (2011) Environmental flows and levels for groundwater-dependent fens of the antelope grazing allotment, Fremont-Winema National Forest, Oregon. The Nature Conservancy, Arlington, VA

Arnaud L, Leclerc B (2010) Forages 1165ZZ0012/F4 et 1165ZZ0013/F5: travaux et pompages d'essais, commune de Grand-Rivière [Boreholes 1165ZZ0012/F4 and 1165ZZ0013/F5: drilling and pumping tests results, municipality of Grand-Rivière]. BRGM open file report RP-58633-FR, BRGM, Orléans, France

Arnaud L, Lanini S (2014) Impact du changement climatique sur les ressources en eau de Martinique [Impact of climate change on Martinique water ressources]. Rapport BRGM/RP-62676-FR, pp 91

Auken E, Christiansen AV, Westergaard JH, Kirkegaard C, Foged N, Viezzoli A (2009) An integrated processing scheme for high-resolution airborne electromagnetic surveys, the SkyTEM system. Explor Geophys 40:184-192

Barat A (1986) Étude du rôle des eaux souterraines dans le mécanisme des éruptions phréatiques. Application à la Montagne Pelée de Martinique et à la Soufrière de Guadeloupe [Study of the role of groundwater in the mechanism of phreatic eruptions: application to the Montagne Pelée of Martinique and the Soufrière de Guadeloupe]. BRGM report no. 115, BRGM, Orléans, France, 205 pp

Bogie I, Mackenzie K (1998) The application of a volcanic facies model to an andesitic stratovolcano hosted geothermal system at Wayang Windu, Java, Indonesia. In: Proceedings 20th NZ Geothermal Workshop, Aukland, 1998, pp 265-270

Boudon G, Balcone-Boissard H (2021) Volcanological evolution of Montagne Pelée (Martinique): a textbook case of alternating plinian and dome-forming eruptions, earth-science reviews, volume 221. https://doi.org/10.1016/j.earscirev.2021.103754

Boudon G, Le Friant A, Komorowski JC, Deplus C, Semet MP (2007) Volcano flank instability in the Lesser Antilles arc: diversity of scale, processes, and temporal recurrence. J Geophys Res: Solid Earth 112(B8)

Brown SK, Jenkins SF, Sparks RSJ, Odbert H, Auker MR (2017) Volcanic fatalities database: analysis of volcanic threat with distance and victim classification. J Appl Volcanol 6:15. https://doi.org/10.1186/s13617-017-0067-4

Cabrera MC, Custodio E (2019) Conceptual hydrogeological models of oceanic intraplate volcanic small. Islands proceedings of IAH2019, the 46th annual congress of the International Association of Hydrogeologists, Málaga, Spain, 22-27 September 2019

Custodio E (1975) Hydrogeologia de las rocas volcanicas [Hydrogeology of volcanic rocks]. 3rd UNESCO-ESA-IHA symposium on groundwater, UNESCO, Paris, pp 23-69

Custodio E, Lopez Garcia L, Amigo E (1988) Simulation par modèle mathématique de l'île volcanique de Ténériffe (Canaries, Espagne) [Mathematical model simulation of the volcanic island of Tenerife (Canary Islands, Spain)]. Hydrogéologie 1988(2):153-167

Deparis J, Reninger PA, Perrin J, Martelet G, Audru JC (2014) Acquisition géophysique héliportée de la Martinique [Helicopter-borne geophysical acquisition of Martinique]. Open file BRGM Report RP-62428-FR. http://infoterre.brgm.fr/rapports/RP-62428-FR.pdf. Accessed May 2023

Dumont M, Peltier A, Roblin E, Reninger PA, Barde-Cabusson S, Fininzola A, Ferrazzini V (2019) Imagery of internal structure and destabilization features of active volcano by 3D high resolution airborne electromagnetism. Sci Rep 9:18280. https://doi.org/10.1038/s41598-019-54415-4

Domenico PA, Schwartz FW (1990) Physical and chemical hydrogeology. Wiley, New York, 824 pp

d'Ozouville N, Auken E, Sorensen KI, Violette S, de Marsily G (2008) Extensive perched aquifer and structural implications
revealed by 3D resistivity mapping in a Galapagos volcano. Earth Planet Sci Lett 269:517–521

Fabre R, Chaigneau M (1960) Technique de prélèvement et d'analyse de gaz de fumerolles de volcans. Exemple d'analyses de gaz prélevés à la montagne Pelée (Martinique) et à la Grande Soufrière (Guadeloupe) [Technique of sampling and analysis of fumarole gases from volcanoes: example of gas analysis from mount Pelee (Martinique) and Grande Soufrière (Guadeloupe)]. Bull Volcanol 23:21–30. https://doi.org/10.1007/BF02596625

Fontaine FR, Corbeau J, Vaton D, the OVSM-IPGP team (2022) Rapport d'activité 2021 de l'Observatoire Volcanologique et Sismologique de Martinique [Activity report 2021 of the Volcanological and Seismological Observatory of Martinique). Institut de Physique du Globe de Paris, Paris

Gadalia A, Baltassat JM, Bouchot V, Caritg S, Coppo N, Gal F, Girard JF, Gutierrez A, Jacob T, Martelet G, Rad S, Tailame AL, Traineau H, Vittecoq B, Wawrzyniak P, Zammit C (2014) Compléments d'exploration géothermique en Martinique: conclusions et recommandations pour les zones de la Montagne Pelée, des Anses d'Arlet, des Pitons du Carbet et du Lamentin [Complementary geothermal exploration in Martinique: conclusions and recommendations for the Montagne Pelée, Anses d'Arlet, Pitons du Carbet and Lamentin areas]. BRGM report RP-63019-FR, BRGM, Orléans, France

Germa A, Quidelleur X, Labanieh S, Lahitte P, Chauvel C (2010) The eruptive history of Morne Jacob volcano (Martinique Island, French West Indies): geochronology, geomorphology and geochemistry of the earliest volcanism in the recent Lesser Antilles Arc. J Volcanol Geotherm Res 198:297–310

Germa A, Quidelleur X, Labanieh S, Chauvel C, Lahitte P (2011a) The volcanic evolution of Martinique Island: insights from K-Ar dating into the Lesser Antilles arc migration since the Oligocene. J Volcanol Geotherm Res 208(2011):122–135

Germa A, Quidelleur X, Lahitte P, Labanieh S, Chauvel C (2011b) The K-Ar Cassignol-Gillot technique applied to Western Martinique lavas: a record of Lesser Antilles arc activity from 2 ma to Mount Pelée volcanism. Quat Geochronol 6:341–355. https://doi.org/10.1016/j.quageo.2011.02.001

Guiscafe J, Klein JC, Moniod F (1976) Les ressources en eau de surface de la Martinique [The surface water resources of Martinique]. Monographies hydrologiques ORSTOM, Noumea, France

Hemmings B, Whitaker F, Gottsmann J, Hughes A (2015) Hydrogeology of Montserrat, review and new insights. J Hydrol Reg Stud 3:1–30. https://doi.org/10.1016/j.ejrh.2014.08.008

Lacroix A (1904) La Montagne Pelée et ses eruptions [Mount Pelee and its eruptions]. Masson, Paris, 662 pp

Le Friant A, Boudon G, Deplus C, Villemant B (2003) Large scale flank-collapse during the recent activity of Montagne Pelee, Martinique. J Geophys Res 108(B1):2055

Legendre Y (2012) Reconstruction fine de l'histoire éruptive et scenarii éruptifs à la soufrière de Guadeloupe: vers un modèle intégré de fonctionnement du volcan [Fine reconstruction of the eruptive history and eruptive scenarios at the Guadeloupe Soufrière: towards an integrated model of the volcano]. Thèse, Université Paris Cité, 675 pp

Lipman PW, Mullineaux DR (1981) The 1980 eruptions of Mount St. Helens. US Geol Surv Prof Pap 1250, 844 pp. https://pubs.usgs.gov/pp/1250/report.pdf. Accessed May 2023

Macdonald GA, Abbott AT, Peterson FL (1983) Volcanoes in the sea: the geology of Hawaii. University of Hawaii Press, Honolulu, 571 pp

Meinzer OE (1930) Ground water in the Hawaiian Islands, in geology and water resources of the Kau District, Hawaii. US Geol Surv Prof Pap 616:1–28

Michaud-Dubuy A (2019) Dynamique des eruptions Pliniennes: réévaluation de L'aléa Volcanique en Martinique [Dynamics of Plinian eruptions: reassessment of the volcanic Hazard in Martinique]. Thèse, Université de Paris, Paris, 200 pp

Peruzzetto M, Legendre Y, Nachbaur A, Dewez TJB, Thiery Y, Levy C, Vittecoq B (2022) How volcanic stratigraphy constrains headscarp collapse scenarios: the Samperre Cliff case study (Martinique island, Lesser Antilles). Nat Hazards Earth Syst Sci 22:3973–3992. https://doi.org/10.5194/nhess-22-3973-2022

Peterson FL (1972) Water development on tropic volcanic islands: type example—Hawaii. Ground Water 10(5):18–23

Pyle DM, Barclay J, Armijos MT (2018) The 1902–3 eruptions of the Soufrière, St Vincent: impacts, relief and response. J Volcanol Geothermal Res 356:183–199. https://doi.org/10.1016/j.jvolgfores.2018.03.005

Pryet A (2011) Hydrogeology of volcanic islands: a case-study in the Galapagos archipelago (Ecuador). PhD Thesis, Université Paris 6 Pierre et Marie Curie, Paris

Pryet A, Ramm J, Chiles JP, Auken E, Deffontaines B, Violette S (2011) 3D resistivity gridding of large AEM datasets: a step toward enhanced geological interpretation. J Appl Geophys 75:277–283

Pryet A, D'ozouville N, Violette S, Deffontaines B, Auken E (2012) Hydrogeological settings of a volcanic island (San Cristobal, Galapagos) from joint interpretation of airborne electromagnetics and geomorphological observations. Hydrol Earth Syst Sci 16(12):4571–4579

Rad S, Rivé K, Vittecoq B, Cerdan O, Allègre JC (2013) Chemical weathering and erosion rates in the Lesser Antilles: an overview in Guadeloupe, Martinique and Dominica. J South Am Earth Sci 45:331–344. https://doi.org/10.1016/j.jsames.2013.03.004

Reninger PA, Martelet G, Deparis J, Perrin J, Chen Y (2011) Singular value decomposition as a denoising tool for airborne time domain electromagnetic data. J Appl Geophys 75:264–276. https://doi.org/10.1016/j.jappgeo.2011.06.034

Reninger PA, Martelet G, Perrin J, Dumont M (2020) Processing methodology for regional AEM surveys and local implications. Explor Geophys 51(1):143–154. https://doi.org/10.1080/08123985.2019.1680249

Robertson REA (1995) An assessment of the risk from future eruptions of the Soufrière volcano of St Vincent, West Indies. Nat Hazards 11:163–191

Robins NS, Lawrence AR, Cripps AC (1990) Problems of groundwater development in small volcanic islands in the eastern Caribbean. In: Krishna JH, Quiñones-Aponte V, Gómez-Gómez F, Morris G (eds) Tropical hydrology and Caribbean water resources. American Water Resources Association, Bethesda, MD, pp 257–267

Rosas-Carbajal M, Komorowski JC, Nicollin F, Gibert D (2016) Volcano electrical tomography unveils edifice collapse hazard linked to hydrothermal system structure and dynamics. Sci Rep 6:29899. https://doi.org/10.1038/srep29899

Selles A (2014) Multidisciplinary study about the hydrogeological behavior of the eastern flank of Merapi volcano, Central Java, Indonesia. PhD Thesis, Université Paris 6 Pierre et Marie Curie, Paris

Selles A, Deffontaines B, Hendrayana H, Violette S (2015) The eastern flank of the Merapi volcano (Central Java, Indonesia): architecture and implications of volcaniclastic deposits. J Asian Earth Sci 108(2015):33–47

Sørensen KI, Auken E (2004) SkyTEM: a new high-resolution helicopter transient electromagnetic system. Explor Geophys 35:191–199

Spies BR (1989) Depth of investigation in electromagnetic sounding methods. Geophysics 54:872–888. https://doi.org/10.1190/1.1442716
Taïlamé AL, Lanini SS (2020) – Définition des volumes prélevables en Martinique. Rapport final [Definition of withdrawable water volumes in Martinique. Final Report]. BRGM/RP-68883-FR, pp 49

Traineau H, Westercamp D, Benderitter Y (1989) Case study of a volcanic geothermal system, Mount Pelée, Martinique. J Volcanol Geotherm Res 38:49–66

Vessell RK, Davies DK (1981) Nonmarine sedimentation in an active fore arc basin. In: Recent and ancient nonmarine depositional environments: models for exploration: Society of Economic Paleon Tologists and Mineralogists Special Paper 31. SEPM Society for Sedimentary Geology, Broken Arrow, OK, pp 31–45. https://doi.org/10.2110/pec.81.31.0031

Viezzoli A, Christiansen AV, Auken E, Sørensen K (2008) Quasi-3D modeling of airborne TEM data by spatially constrained inversion. Geophysics 73:F105–F113

Villemant B, Le Friant A, Caron B, Del Manzo G, Lafuerza S, Emmanuel L, Ishizuka O, Guyard H, Labourdette N, Hidalgo MA (2022) A 1.5 ma marine record of volcanic activity and associated landslides offshore Martinique (Lesser Antilles): sites U1397 and U1399 of IODP 340 expedition. Front Earth Sci 10:767485. https://doi.org/10.3389/feart.2022.767485

Vincent PM, Bourdier JL, Boudon G (1989) The primitive volcano of Montagne Pelée: its construction and partial destruction by flank collapse. In: Boudon G, Bourgaud A (eds) Montagne Pelée. J Volcanol Geotherm Res 38:1–15

Violette S, d'Ozouville N, Pryet A, Deffontaines B, Fortin J, Adelinet M (2014) Hydrogeology of the Galapagos archipelago: an integrated and comparative approach between islands. In: Harpp KS, Mittelstaedt EL, D'Ozouville N, Graham DW (eds) The Galápagos as a laboratory for the earth sciences, vol 9. American Geophysical Union, Washington, DC, pp 167–183

Vittecoq B, Lachassagne P, Lanini S, Ladouche B, Marechal JC, Petit V (2007) Elaboration d'un système d'information sur les eaux souterraines de la Martinique: identification et caractérisations quantitative [elaboration of an information system on groundwater in Martinique: identification and quantitative characterizations]. Open file BRGM report RP-55099-FR. http://infoterre.brgm.fr/rapports/RP-55099-FR.pdf. Accessed May 2023

Vittecoq B, Lachassagne P, Lanini S, Maréchal JC (2010) Evaluation des ressources en eau de la Martinique : calcul spatialisé de la pluie efficace et validation à l'échelle du bassin-versant [Evaluation of water resources in Martinique: spatialized calculation of effective rainfall and validation at the watershed scale]. Rev Sci Eau 23(4):361–373

Vittecoq B, Deparis J, Violette S, Jaouen T, Lacquement F (2014) Influence of successive phases of volcanic construction and erosion on Mayotte Island's hydrogeological functioning as determined from a helicopter-borne resistivity survey correlated with borehole geological and permeability data. J Hydrol 509:519–538. https://doi.org/10.1016/j.jhydrol.2013.11.062

Vittecoq B, Reninger PA, Violette S, Martele G, Dewandel B, Audru JC (2015) Heterogeneity of hydrodynamic properties and groundwater circulation of a coastal andesitic volcanic aquifer controlled by tectonic induced faults and rock fracturing: Martinique Island (Lesser Antilles - FWI). J Hydrol 529:1041–1059. https://doi.org/10.1016/j.jhydrol.2015.09.022

Vittecoq B, Reninger PA, Lacquement F, Martelet G, Violette S (2019) Hydrogeological conceptual model of andesitic watersheds revealed by high-resolution heliborne geophysics. Hydrol Earth Syst Sci 23:2321. https://doi.org/10.5194/hess-23-2321-2019

Wadge G, Voight B, Sparks RSJ, Cole PD, Loughlin SC, Robertson REA (2014) An overview of the eruption of Soufrière Hills volcano, Montserrat from 2000 to 2010. Memoir 39, Geological Society, London, 40 pp

Ward SH, Hohmann GW (1988) Electromagnetic theory for geophysical applications: electromagnetic methods. Appl Geophys 1:130–311

Westercamp D, Traineau H (1983) - Carte géologique de la Montagne Pelée à 1/20 000 (département de la Martinique) et notice explicative simplifiée [Montagne Pelée 1/20 000 geological map (Martinique) with simplified notice]. Éd. BRGM

Westercamp D, Andreieff P, Bouysse P, Cottez S, Battistini R (1989) Notice explicative, carte géol. France (1/50 000), feuille Martinique (explanatory note, Geol. Map France (1/50 000), sheet Martinique. BRGM, Orléans, France, 246 pp

Westercamp D, Pelletier B, Thibaut PM, Traineau H, Andreieff P (1990) Carte géologique de la France (1/50.000), feuille Martinique [Geological map of France (1/50,000), Martinique sheet]. BRGM, Orléans, France. http://infoterre.brgm.fr/. Accessed May 2023

Zlotnicki J, Boudon G, Viodé JP, Delarue JF, Mille A, Bruère F (1998) Hydrothermal circulations beneath Montagne Pelée inferred by self-potential surveying: structural and tectonic implications. J Volcanol Geotherm Res 84:73–91

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.