3-D Modelling of Campi Flegrei Ground Deformations: Role of Caldera Boundary Discontinuities

Campi Flegrei is a caldera complex located west of Naples, Italy. The last eruption occurred in 1538, although the volcano has produced unrest episodes since then, involving rapid and large ground movements (up to 2 m vertical in two years), accompanied by intense seismic activity. Surface ground displacements detected by various techniques (mainly InSAR and levelling) for the 1970 to 1996 period can be modelled by a shallow point source in an elastic half-space, however the source depth is not compatible with seismic and drill hole observations, which suggest a magma chamber just below 4 km depth. This apparent paradox has been explained by the presence of boundary fractures marking the caldera collapse. We present here the first full 3-D modelling for the unrest of 1982–1985 including the effect of caldera bordering fractures and the topography. To model the presence of topography and of the complex caldera rim discontinuities, we used a mixed boundary elements method. The a priori caldera geometry is determined initially from gravimetric modelling results and refined by inversion. The presence of the caldera discontinuities allows a fit to the 1982–1985 levelling data as good as, or better than, in the continuous half-space case, with quite a different source depth which fits the actual magma chamber position as seen from seismic waves. These results show the importance of volcanic structures, and mainly of caldera collapses, in ground deformation episodes.     

Simultaneous inversion of deformation and gravity changes in a horizontally layered half-space: Evidences for magma intrusion during the 1982–1984 unrest at Campi Flegrei caldera (Italy)

A very large uplift (about 1.8 m) occurred in the period 1982–1984 at Campi Flegrei caldera, Italy, without culminating in an eruption. A still-standing controversy accompanies the interpretation of deformation and gravity changes recorded during the unrest, which were interpreted to result from the sub-surface magmatic reservoir by some authors and from the hydrothermal system or hybrid sources by others. Here for the first time we take into account crustal layering while inverting leveling, EDM, and gravity data using uniformly-pressurized sources, namely small vertical spheroids and finite horizontal penny-shaped sources. The weight of EDM data in the misfit function is chosen from a trade-off curve in order to balance the compromise between fitting the leveling and the EDM data well. Models using a homogeneous medium cannot give a good simultaneous fit to leveling and EDM deformation data of the 1982–1984 unrest, whereas incorporating a layered structure (determined from seismically derived estimates of the P wave speed for the crust, and not adjusted to improve the fit in any of the inversions) allows a significantly better fit. Also, layering affects the sub-surface mass redistribution effects on gravity changes, and we show that the retrieved intrusion density is in full agreement with densities for highly evolved magmas expected at the Campi Flegrei caldera for depths of 3 to 4 km, ruling out hydrothermal fluids as the primary cause of the 1982–1984 unrest. The source of the 1982–1984 CF unrest was probably a shallow (about 3-km deep) penny-shaped magma intrusion fed by a deeper magma chamber; source overpressure was few MPa.     

Evaluating fracture patterns within a resurgent caldera: Campi Flegrei, Italy

Understanding deformation of active calderas allows their dynamics to be defined and their hazard mitigated. The Campi Flegrei resurgent caldera (Italy) is one of the most active and hazardous volcanoes in the world, characterized by post-collapse resurgence, eruptions, ground deformation, and seismicity. An original structural analysis provides an overview of the main fracture zones. NW-SE and NE-SW fractures (normal or transtensive faults and extensional fractures) predominate along the rim and within the caldera, suggesting a regional control, both during and after the collapses. While the NE-SW fractures are ubiquitous in the deposits of the last ∼37 ka, NW-SE fractures predominate in the last 4.5 ka, during resurgence. The most recently (<4.5 ka) strained area lies in the caldera center (Solfatara area), where the faults, with an overall ∼ENE-WSW extension direction, appear to be associated with the bending due to resurgence. Solfatara lies immediately to the east of the most uplifted part of the caldera (Pozzuoli area), where domes form and culminate both on the long-term (resurgence, accompanied by volcanic activity) and short-term deformation (1982–1984 bradyseism, accompanied by seismic and hydrothermal activity). Similar volcano-tectonic behavior characterizes the short- and long-term uplifts, and only the intensity of the tectonic and volcanic activity varies, being related to varying amounts of uplift. Seismicity and hydrothermal manifestations occur during the bradyseisms, with moderate uplift, while surface faulting and eruptions occur during resurgence, with higher uplift. The features observed at Campi Flegrei are found at other major calderas, suggesting consistent behavior of large magmatic systems.     

Vertical ground movements in the Campi Flegrei caldera as a chaotic dynamic phenomenon

The ground level in the Campi Flegrei caldera has never been stationary in the last 2,000 years. Historical data, and a nearly continuous tide-gauge record 20 years long, show that uplift and sinking have taken place on a variety of different time scales. In addition, the Campi Flegrei volcanic system appears to be sensitive to weak external forces such as tidal forces. We infer from these elements that the Campi Flegrei system is far from thermodynamic equilibrium, and suggest that its dynamics may be chaotic. We analyze the short-term variations of the ground level, and find that they can be described in a low-dimensional phase space. The dynamics of the Campi Flegrei system seems to have been phase-locked with tidal forces in the period following the 1970–1972 climax, and to have undergone a transition to chaos in some moment that preceded the presently continuing sinking phase.     

A mechanical fluid-dynamical model for ground movements at Campi Flegrei caldera

We present here a consistent model, which explains the mechanisms of unrest phenomena at Campi Flegrei (Italy), both at short-term (years) and at secular scales. The model consists basically of two effects: the first one is related to the elastic response of the shallow crust to increasing pressure within a shallow magma chamber; the second involves the fluid-dynamics of shallow aquifers in response to increasing pressure and/or temperature at depth. The most important roles in the proposed model are played by the effect of lateral stress–strain discontinuities marking the inner caldera borders and by the response of the geothermal system to the increase of pressure and/or temperature coming from the magma chamber. The model takes into account most results of recent research, and it is able to explain in a unitary way most of geophysical observations recently collected at Campi Flegrei. Some new results, presented in this paper, explain the mechanism of generation of seismicity in terms of Coulomb stress changes, and the difference between secular and short-term deformations, in terms of the non linear effect of the bordering discontinuities. The coupled effects of magma chamber stresses with the dynamics of shallow fluids provide results for a semi-quantitative interpretation of observed ground deformation, seismicity and of their time evolution. The model has very important implications for the hazard assessment during unrest phenomena. It represents, for the part involving the fluid-dynamical response, the ideal continuation of the basic, fundamental models and intuitions of Oliveri del Castillo and coworkers (Oliveri del Castillo, A., Quagliariello, M.T., 1969. Sulla genesi del bradisimo flegro. Atti Associazione Geofisica Italiana, 18th Congress, Napoli, pp. 557–594; Casterano, L., Oliveri del Castillo, A., Quagliariello, M.T., 1976. Hydromdynamics and geodynamics in the Phlegraean Fields area of Italy. Nature, 264, 161–154).     

Structure and evolution of the Bay of Pozzuoli (Italy) using marine seismic reflection data: implications for collapse of the Campi Flegrei caldera

Digital marine seismic reflection data acquired in 1973 in the Bay of Pozzuoli, and recently reprocessed, were used to study the volcanological evolution of the marine sector of Campi Flegrei Caldera during the last 37 ka. In order to gain more information, interpretation also involved estimation of the "pseudo-velocity" and the "pseudo-density" from the resistivity logs of two onshore deep exploration wells. The main results are: (1) discovery of ancient pre-18 ka and post-37 ka submarine and mainly effusive volcanic activity, along coeval emission centers located at the edges of Campi Flegrei Caldera; (2) confirmation that the caldera collapse in the marine sector of Campi Flegrei seems strongly controlled by regional NE–SW and NW–SE structural discontinuities; (3) the finding of at least two episodes of collapse in the bay; and (4) identification of a post-18 ka volcanic deflation phase that has caused about 150–200 m of subsidence in the central sector of the Bay of Pozzuoli in the last 18 ka.Editorial responsibilty: T. Druitt     

Modelling gravity variations consistent with ground deformation in the Campi Flegrei caldera (Italy)

In recent years (1970–72 and 1982–84) two inflation episodes took place in the Campi Flegrei caldera (Italy), characterized by significant ground uplift and gravity variations. An elastic half-space model with vertical density stratification is employed to compute the displacement field and the gravity variations produced by the deformation of buried layers, following the inflation of a spherically symmetric deformation source. Contributions to gravity variations are produced by dilation/contraction of the medium, by the displacements of density interfaces (the free surface and subsurface layers) and of source boundaries and, possibly, by new mass input from remote distances into the source volume. Three cases were examined in detail: In case I, the magma chamber is identified as the deformation source and volume and pressure increase in the magma chamber is due to input of new magma from remote distances; in case II deformation is due to magma differentiation within the magma chamber (deformation source with constant mass); in case III the geothermal system is identified as the deformation source and a pressure increase, possibly driven by the exsolution of high temperature and high pressure volatiles in the magma chamber, is assumed to play a dominant role. From the comparison between measured and computed gravity residuals (free-air-corrected gravity variations) we can assess that, in case I, an inflation source with constant density would predict gravity residuals compatible with observations, whereas an expansion at constant mass (case II) would predict gravity residuals much lower than observed. The resolving power of gravity data however prevents accurate assessment of the density of the emplaced material. In case III, the pervasive density increase of the geothermal fluids induced by pressure increase is assumed to be the main source of gravity variations. The average porosity value required for this model to match both the ground deformation and the gravity residuals is found to be ˜10%, a value which is compatible with measured porosity values at Campi Flegrei in deep wells. The subsidence phases following both inflation episodes and the gravity residuals during subsidence lead us to consider case III as more plausible, even if a suitable combination of cases I and III cannot be discarded.     

Hazard assessment at volcanic fields: the Campi Flegrei case history

Volcanological analysis of the 10 000 yr –1538 explosive activity at Campi Flegrei shows that the most common explosive eruptions are characterized by the emplacement of flow or surge deposits, originating from the interaction between magma and shallow and/or sea water. The minimum volumes of pyroclastic products range between 0.04 and 0.7 km³; the proximal areas covered by these products range from 3–4 to 40–50 km². The pyroclastic flow and surge deposits occurring inside the caldera have been strongly controlled by pre-existent morphology; because of this, the area of present Napoli city was blanketed by approximately 5 m of pyroclastic deposits, during the last 5000 yr.Previous analysis suggests that the presence of even very low topographic obstacles may influence pyroclastic density current run out such that future eruptive deposits would mainly be confined inside the caldera rim. We suggest that a future eruption at Campi Flegrei would not seriously involve the urbanized area of Napoli city located on the hills. On the contrary, the plains located on the eastern side of the caldera (Fuorigrotta, Bagnoli) would be the most damaged area.     

Unrest episodes at Campi Flegrei: A reconstruction of vertical ground movements during 1905–2009

Geodetic observations at Campi Flegrei caldera were initiated in 1905. Historical observations and the few measurements made before 1970 suggested a deflationary trend. Since 1969, the ground started to inflate during two major uplift episodes in 1969–72 and 1982–1985. We collected and reanalyzed all available punctual observations of vertical ground displacement taken in the period 1905–2009 with special attention to the period before 1969, to reconstruct in greater detail the deformation history of the caldera. We make use of the many photographs of the sea level in a roman ruin (the Serapeum Market) taken during the period between 1905 and 1969 to infer with more accuracy its relative height with respect to the sea level. We identify a previously disregarded major episode of ground uplift occurred between 1950 and 1952 with a maximum uplift of about 73 cm. This finding suggests that Campi Flegrei is currently experiencing a prolonged period of unrest longer than previously thought. The higher seismicity associated with the later episodes of unrest suggests that the volcano has approached an instability threshold, which may eventually result in a volcanic eruption.     

Probability hazard map for future vent opening at the Campi Flegrei caldera,Italy

The Campi Flegrei caldera is a restless structure affected by general subsidence and ongoing resurgence of its central part. The persistent activity of the system and the explosive character of the volcanism lead to a very high volcanic hazard that, combined with intense urbanization, corresponds to a very high volcanic risk. One of the largest sources of uncertainty in volcanic hazard/risk assessment for Campi Flegrei is the spatial location of the future volcanic activity. This paper presents and discusses a long-term probability hazard map for vent opening in case of renewal of volcanism at the Campi Flegrei caldera, which shows the spatial conditional probability for the next vent opening, given that an eruption occurs. The map has been constructed by building a Bayesian inference scheme merging prior information and past data. The method allows both aleatory and epistemic uncertainties to be evaluated. The probability map of vent opening shows that two areas of relatively high probability are present within the active portion of the caldera, with a probability approximately doubled with respect to the rest of the caldera. The map has an immediate use in evaluating the areas of the caldera prone to the highest volcanic hazard. Furthermore, it represents an important ingredient in addressing the more general problem of quantitative volcanic hazards assessment at the Campi Flegrei caldera.     

The Campania-Campi Flegrei area: a contribution to discern the best structural model from gravity interpretation

Structural interpretations of negative gravity anomalies at Campi Flegrei and Campania are reviewed. The interpretations are not uniquely determined, as different models correspond to equally good anomaly fits.These models are examined using a new method based on the interpretation of the apparent density map computed from the gravity anomaly map.As regards the Campi Flegrei anomaly, a funnel-shaped caldera 3D model is readily obtained. This is more appropriate than previous interpretations, as not only anomaly but also apparent densities are fitted.As about the regional structures, the used approach suggests that the more pertinent model is a large magmatic intrusion, 9 to 14 km deep.A final scheme is proposed linking this intrusion to the funnel-shaped caldera as the main regional reservoir.     

Seismic and ground deformation crises at Rabaul Caldera: Prelude to an eruption?

A dramatic short-term increase in seismicity and ground deformation took place at Rabaul Caldera on 19 September 1983, and marked the start of a period of frequent episodes of high seismic energy release and concurrent rapid ground deformation. Together with increased background levels of seismicity and ground deformation, these phenomena are interpreted as indications of higher rates of magma injection at shallow depths within the caldera, which greatly increases the likelihod of an eruption at Rabaul in the near future. A modest volume of magma, about 100 million m³, could be available for eruption from two shallow reservoirs, but a somewhat deeper and much larger magma body — residual from the latest major eruption about 1400 yr BP — may also exist beneath the caldera.     

The Neapolitan Yellow Tuff — A large volume multiphase eruption from Campi Flegrei,Southern Italy

the Neapolitan Yellow Tuff (NYT) (12 ka BP) is considered to be the product of a single eruption. Two different members (A and B) have been identified and can be correlated around the whole of Campi Flegrei. Member A is made up of at least 6 fall units including both ash and lapilli horizons. The basal stratified ash unit (A1) is interpreted to be a phreatoplinian fall deposit, since it shows a widespread dispersal (>1000 km²) and a constant thickness over considerable topography. The absence of many lapilli fall units in proximal and medial areas testifies to the erosive power of the intervening pyroclastic surges. The overlying member B was formed by many pyroclastic flows, radially distributed around Campi Flegrei, that varied widely in their eruptive and emplacement mechanisms. In some of the most proximal exposures coarse scoria and lithic-rich deposits, sometimes welded, have been identified at the base of member B. Isopach and isopleth maps of fall-units, combined with the distribution of the coarse proximal facies, indicate that the eruptive vent was located in the NE area of Campi Flegrei. It is considered that the NYT eruption produced collapse of a caldera approximately 10 km diameter within Campi Flegrei. The caldera rim, located by geological and borehole evidence, is now largely buried by the products of more recent eruptions. Initiation of caldera collapse may have been contemporaneous with the start of the second phase (member B). It is suggested that there was a single vent throughout the eruption rather than the development of multiple or ring vents. Chemical data indicate that different levels of a zoned trachyte-phonolite magma chamber were tapped during the eruption. The minimum volume of the NYT is calculated to be about 50 km³ (DRE), of which 35 km³ (70%) occurs within the caldera.     

The structure of the Campanian Plain and the activity of the Neapolitan volcanoes (Italy)

The central Campanian Plain is dominated by the structural depression of Acerra whose origin is tectonic, but may have been enlarged and further depressed after the eruption of the Campanian Ignimbrite (42-25 ka). The deposits of the Campanian Ignimbrite are possibly the results of multiple eruptions with huge pyroclastic deposits that covered all the Campanian Plain.The more recent activity of Vesuvius, Campi Flegrei and Procida occurred on the borders of Acerra depression and resulted from a reactivation of regional faults after the Campanian Ignimbrite cycle. The activity of Vesuvius produced the building of a stratovolcano mostly by effusive and plinian explosive eruptions. The Campi Flegrei area, on the contrary, was dominated by the eruption of the Neapolitan Yellow Tuff at 12 ka that produced a caldera collapse of the Gulf of Pozzuoli. The caldera formation controlled the emplacement of the recent activity of Campi Flegrei and the new volcanoes were formed only within the caldera or along its rim.     

Eruptive versus non-eruptive behaviour of large calderas: the example of Campi Flegrei caldera (southern Italy)

Caldera eruptions are among the most hazardous of natural phenomena. Many calderas around the world are active and are characterised by recurrent uplift and subsidence periods due to the dynamics of their magma reservoirs. These periods of unrest are, in some cases, accompanied by eruptions. At Campi Flegrei caldera (CFc), which is an area characterised by very high volcanic risk, the recurrence of this behaviour has stimulated the study of the rock rheology around the magma chamber, in order to estimate the likelihood of an eruption. This study considers different scenarios of shallow crustal behaviour, taking into account the earlier models of CFc ground deformation and caldera eruptions, and including recent geophysical investigations of the area. A semi-quantitative evaluation of the different factors that lead to magma storage or to its eruption (such as magma chamber size, wall-rock viscosity, temperature, and regional tectonic strain rate) is reported here for elastic and viscoelastic conditions. Considering the large magmatic sources of the CFc ignimbrite eruptions (400–2,000 km³) and a wall-rock viscosity between 10¹⁸ and 10²⁰ Pa s, the conditions for eruptive failure are difficult to attain. Smaller source dimensions (a few cubic kilometres) promote the condition for fracture (eruption) rather than for the flow of wall rock. We also analyse the influence of the regional extensional stress regime on magma storage and eruptions, and the thermal stress as a possible source of caldera uplift. The present study also emphasises the difficulty of distinguishing eruption and non-eruption scenarios at CFc, since an unambiguous model that accounts for the rock rheology, magma-source dimensions and locations and regional stress field influences is still lacking.     

Hot fluid migration: an efficient source of ground deformation: application to the 1982–1985 crisis at Campi Flegrei-Italy

Some solutions of the forced heat advection problem in compressible media are worked out employing a perturbative approach and their implication for thermoelastic deformation are discussed. A sharp temperature front, which migrates at a speed in the order of Darcy flow rate, develops in the medium, giving rise to significant deformation via thermal expansion. A thermally induced pressure source accompanies the temperature front, which may be significant only in very high temperature cases. Results are applied to interpreting the uplift episode of 1982–1984 at Campi Flegrei (C.F.), near Naples, Italy. A mechanism is envisaged for uplift at Campi Flegrei in which a sudden connection is established between a deep, hot, high-pressure fluid reservoir and a shallow, relatively cold, low-pressure aquifer. The inclusion of fluid migration in the deformation model allows simple explanations of several geochemical and geophysical observations made during the bradyseismic crisis. It appears that the proposed mechanism may explain the large observed uplift, without requiring unreasonable pressure increase within the magma chamber. Furthermore, the deformation source may be allowed to be shallower than the magma chamber, as required by any reasonable deformation model at C.F.     

Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.     

Understanding the Seismic Velocity Structure of Campi Flegrei Caldera (Italy): From the Laboratory to the Field Scale

We report laboratory measurements of P- and S-wave velocities on samples of tuff from Campi Flegrei (Italy), and a new tomographic velocity map of the Campi Flegrei caldera. Laboratory measurements were made in a hydrostatic pressure vessel during both increasing and decreasing effective pressure cycles. Selected samples were also thermally stressed at temperatures up to 600°C to induce thermal crack damage. Acoustic emission output was recorded throughout each thermal stressing experiment, and velocities were measured after thermal stressing. Laboratory P- and S-wave velocities are initially low for the tuff, which has an initial porosity of ~45%, but both increase by between 25 and 50% over the effective pressure range of 5 to 80 MPa, corresponding to a decrease of porosity of ~70%. Marked velocity hysteresis, due to inelastic damage processes, is also observed in samples subjected to a pressurization-depressurization cycle. Tomographic seismic velocity distributions obtained from field recordings are in general agreement with the laboratory measurements. Integration of the laboratory ultrasonic and seismic tomography data indicates that the tuffs of the Campi Flegrei caldera can be water or gas saturated, and shows that inelastic pore collapse and cracking produced by mechanical and thermal stress can significantly change the velocity properties of Campi Flegrei tuffs at depth. These changes need to be taken into account in accurately interpreting the crustal structure from tomographic data.     

The Christmas Mountains caldera complex,Trans-Pecos Texas

The Christmas Mountains caldera complex developed approximately 42 Ma ago over an elliptical (8×5 km) laccolithic dome that formed during emplacement of the caldera magma body. Rocks of the caldera complex consist of tuffs, lavas, and volcaniclastic deposits, divided into five sequences. Three of the sequences contain major ash-flow tuffs whose eruption led to collapse of four calderas, all 1–1.5 km in diameter, over the dome. The oldest caldera-related rocks are sparsely porphyritic, rhyolitic, air-fall and ash-flow tuffs that record formation and collapse of a Plinian-type eruption column. Eruption of these tuffs induced collapse of a wedge along the western margin of the dome. A second, more abundantly porphyritic tuff led to collapse of a second caldera that partly overlapped the first. The last major eruptions were abundantly porphyritic, peralkaline quartz-trachyte ash-flow tuffs that ponded within two calderas over the crest of the dome. The tuffs are interbedded with coarse breccias that resulted from failure of the caldera walls. The Christmas Mountains caldera complex and two similar structures in Trans-Pecos Texas constitute a newly recognized caldera type, here termed a laccocaldera. They differ from more conventional calderas by having developed over thin laccolithic magma chambers rather than more deep-seated bodies, by their extreme precaldera doming and by their small size. However, they are similar to other calderas in having initial Plinian-type air-fall eruption followed by column collapse and ash-flow generation, multiple cycles of eruption, contemporaneous eruption and collapse, apparent pistonlike subsidence of the calderas, and compositional zoning within the magma chamber. Laccocalderas could occur else-where, particularly in alkalic magma belts in areas of undeformed sedimentary rocks.     

Hydrothermal circulation beneath Mount Pelée inferred by self potential surveying. Structural and tectonic implications

Self-potential (SP) surveys were made on Mount Pelée volcano (Martinique Island, French West Indies) in 1991 and 1992 in order to recognize its hydrothermal system, the associated groundwater channeling and the main superficial structures of the massif. Almost 70 km of profiles were carried out with an average sample spacing of 50 m. Measurements essentially reveal negative SP anomalies, down to −1700 mV, with high gradients (−1.83 mV/m) due to the infiltration of meteoric water into the massif. Rims of summit calderas Morne Macouba and Etang-Sec present sharp negative SP anomalies on the western, northern, and eastern flanks. Negative SP anomalies indicate no upward water flow beneath Mount Pelée summit. On the southwestern volcano flank, a 3.5×6 km horseshoe-shaped structure corresponding to a southwest flank collapse event, older than 25,000 years BP, is clearly identified by the SP mapping. High gradients border the inner southern rim from Morne Calebasse to St Pierre town and the Caribbean Sea. Along the northern rim of the horseshoe-shaped structure the negative SP anomalies give place to a positive SP anomaly, up to 200 mV, of SW–NE trend. This zone covers the area of two active hot springs (Sources Chaudes and Puits Chaud: 40–65°C). Marine magnetic surveys and bathymetry show that the horseshoe-shaped structure spreads into the Caribbean Sea up to about 10 km from the coast. Buried structural discontinuities are evidenced inside the flank collapse structure. The upper one deviates the groundwater flow coming from the summit toward the south flank where the flow finds an indentation to expand again downwards. This discontinuity is either an old hypothetical caldera rim partly destroyed by the collapse of the south–southwestern flank and covered by recent pyroclastic deposits, or more probably the trace of a bulge landslide. A circulation model of the hydrothermal waters is proposed. Rainfall (5–6 m/year) is partly drained inside the summital calderas and the flank collapse zone through pyroclastic flows down to an impermeable basement. There the groundwater constitutes perched aquifers at the contact of the bulge landslide, or of the hypothetical old caldera rim. Along the inner northern border of the flank collapse structure the phreatic water is reheated. Warm groundwater flows along the northern avalanche structure rim and discharges near the coast in ground and marine outcrops, of medium temperature. Finally, the main part of the meteoric water is channeled along the old caldera rim, or along the bulge landslide towards the south flank of Mount Pelée, where some gaps in the rim exist. There the groundwater finds again a subhorizontal gravitational circulation along Mount Pelée slopes into the Caribbean Sea.