# δ13C, CO2/3He and 3He/4He ratios reveal the presence of mantle gas in the CO2-rich groundwaters of the Ardennes massif (Spa, Belgium)

Agathe Defourny, Pierre-Henri Blard, Laurent Zimmermann, Patrick Jobé, Arnaud Collignon, Frédéric Nguyen, Alain Dassargues

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

## To cite this version:

Agathe Defourny, Pierre-Henri Blard, Laurent Zimmermann, Patrick Jobé, Arnaud Collignon, et al. δ13C, CO2/3He and 3He/4He ratios reveal the presence of mantle gas in the CO2-rich groundwaters of the Ardennes massif (Spa, Belgium). Hydrology and Earth System Sciences, 2022, 26, pp.2637-2648. 10.5194/hess-26-2637-2022. insu-03712868

HAL Id: insu-03712868

https://insu.hal.science/insu-03712868

Submitted on 4 Jul 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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
<figure><img src="image_3.png" /></figure>

# $\delta^{13}$C, CO$_2$ / $^3$He and $^3$He / $^4$He ratios reveal the presence of mantle gas in the CO$_2$-rich groundwaters of the Ardennes massif (Spa, Belgium)

Agathe Defourny$^{1,2}$, Pierre-Henri Blard$^{3,4}$, Laurent Zimmermann$^3$, Patrick Jobé$^2$, Arnaud Collignon$^2$, Frédéric Nguyen$^1$, and Alain Dassargues$^1$

**Correspondence:** Agathe Defourny (adefourny@uliege.be)

Received: 3 December 2021 – Discussion started: 17 January 2022

Revised: 8 April 2022 – Accepted: 13 April 2022 – Published: 19 May 2022

**Abstract.** Although natural CO$_2$-rich groundwaters of eastern Belgium have been known for centuries, the exact origin of their gas is still unclear. This paper presents the results of a sampling campaign in Belgium (Spa, Stoumont, Malmedy): 30 samples of both carbogaseous and non-carbogaseous groundwaters were analyzed for major elements, CO$_2$ content and carbon isotopic composition. Among them, 13 samples were also analyzed for $^3$He/$^4$He and $^4$He/$^{20}$Ne ratios. The combination of $\delta^{13}$C (between ca. -9‰ VPDB¹ and -2‰ VPDB), CO$_2$/$^3$He ratio (between $1.9 \times 10^8$ and $2.9 \times 10^9$) and $^3$He/$^4$He (between 0.92 and 2.70 Ra) shows with a high level of confidence that the CO$_2$ in the carbogaseous groundwater of Spa and Bru has a mantle origin. It can likely be attributed to the degassing of mantle from the neighboring Eifel volcanic fields, located at a distance of 100 km eastwards.

The identity and nature of the deep-rooted fractures that act as CO$_2$ transport pathways to the surface are still to be clarified, but several major thrust faults exist in the Rhenish Massif and could connect the Eifel volcanic fields with the studied area.

## 1 Introduction

CO$_2$-rich groundwaters have always been a very fascinating geомanifestation, and their presence in a region is often the trigger of a strong economic and touristic activity. Lately, a better assessment of CO$_2$ circulation modes in the subsurface has gained interest, as it is important to finely document their contribution to the natural budget of atmospheric CO$_2$. Understanding the saturation state of CO$_2$ in groundwaters and in geological reservoirs is also important for CO$_2$ storage projects. Moreover, in areas where CO$_2$-rich groundwater is bottled as mineral water, it is essential to have a complete understanding of the whole system, to ensure sustainable exploitation of the resource. The presence of CO$_2$ in groundwater – in excess compared to the atmospheric equilibrium – can result from different phenomena, with the most common ones being a direct contribution from the mantle, the dissolution of carbonate rocks or an organic origin (Agnew, 2018). The discharge of CO$_2$ from deep geological structures to the surface is always the result of a specific geological context which involves a source of CO$_2$ at depth and an intricate system of faults acting as transport pathways to the local groundwaters, which is considered the final receptor.

This study focuses on the dozens of CO$_2$-rich groundwater springs that exist in the Ardennes massif in eastern Belgium. The most famous ones are located in the small city of Spa, whose springs have been known since the Roman Empire.

<footnote>$^1$Urban and Environmental Engineering, University of Liège, Liège, Belgium</footnote>

<footnote>$^2$Water Resource Department, Spadel S.A., Spa, Belgium</footnote>

<footnote>$^3$Centre de Recherches Pétrographiques et Géochimiques (CRPG), Université de Lorraine, CNRS, UMR7358, Nancy, France</footnote>

<footnote>$^4$Laboratoire de Glaciologie, Université libre de Bruxelles (ULB), Brussels, Belgium</footnote>

<footnote>$^1$Vienna Pee Dee Belemnite</footnote>
The name of the town became famous thanks to the devel-
opment of thermalism in the 19th century and is now used
worldwide to refer to wellness and bathing activities. These
groundwaters yield natural springs at the surface, but they
are also exploited as mineral water from boreholes. Their
CO₂ content (2 ± 0.5 g L⁻¹) makes them slightly acidic (pH
around 5.7). They have a high content of iron (17±10 mg L⁻¹
on average but up to 50 mg L⁻¹). They, however, bear a
relatively low TDS (total dissolved solid) value with a dry
residue ranging between 80 and 160 mg L⁻¹ in compari-
son with other naturally sparkling mineral waters bottled in
other European geological contexts (e.g., 940 mg L⁻¹ for San
Pellegrino, 1100 mg L⁻¹ for Badoit, 3325 mg L⁻¹ for Vichy
Célestins, as stated on bottle labels).

Although these springs have been bottled for centuries and
studied for many years, the origin of their high CO2 content
has not been established already. Helium isotopes, elemental
⁴He/²⁰Ne and CO₂ isotopic composition of dissolved gas are
powerful tools to identify the sources of these gas (Sano and
Marty, 1995; Karolyte et al., 2019; Gilfillan et al., 2019). In
this paper, we present the ³He, ⁴He, ²⁰Ne and CO₂ concen-
trations measured in 13 groundwater samples of the Spa and
Bru areas (Ardennes, Belgium), together with hydrochemi-
cal analysis on major elements for 30 samples, to identify
the origin of dissolved CO₂ in groundwater and to explore
the potential hydro-connection with the Eifel volcanic fields
(western Germany), where similar CO₂-rich groundwaters
are found.

**2 Geological context**

**2.1 Regional geology**

CO₂-rich mineral waters from eastern Belgium are located
in the Rhenish Massif, which is part of the Rhenohercy-
nian fold belt (Vanbrabant et al., 2002). This massif extends
through western Germany, eastern Belgium, Luxembourg,
and a part of France, as shown in Fig. 1. The Rhenish Mas-
sif is dominated by Paleozoic rocks and is separated into
two parts by the Rhine Graben. The western part, the Ar-
dennes (eastern Belgium), is bordered to the north by the
Midi-Eifel thrust fault. In the Ardennes region, the Rhenish
Massif is dominated by the Ardennes Allochthon, consisting
of a Cambro-Ordovician basement unconformably overlain
by Devono-Carboniferous sandstones and limestones (Bar-
ros et al., 2021) (Fig. 1).

The present regional geology is the result of several stages. The oldest rocks that can be observed in Belgium are found in the Ardennes region; they consist of Cambrian to Ordovician sediments deposited in deep-platform marine environments. They mainly consist of fine clays in alternation with sandstones. During the Late Silurian, the Caledonian orogeny faulted and folded these layers and induced a strong metamorphism. Compression and shearing of clay-

stone produced a well-expressed schistosity, while sandstone
evolved into quartzites. The metasediments observed today
are therefore an alternation of clays and sandstones meta-
morphosed into slates ("phyllades") and quartzites, which
are called "quartzophyllades" in the region. Then, after an
emersion and erosion period, sedimentary deposition started
again during the Lower Devonian, in unconformity over the
folded rocks. During the Lower Devonian period, limestones
were deposited in vast carbonate platforms. Then, sandstones
accumulated in coastal detrital environments during the De-
vonian.

At the end of the Devonian, the Variscan orogeny took
place, inducing another phase of tectonic deformation. Dur-
ing this orogeny, the Dinant Synclinorium and the Ardennes
massif were displaced several kilometers northward. This
(latter) unit constitutes the Ardennes Allochthon, a great an-
ticlinal ensemble that is limited to the north by the Midi-Eifel
thrust fault. Three Cambrian massifs are identified within
the Ardennes Allochthon: the Rocroi, the Givonne and the
Stavelot massifs. The studied springs are located within the
Stavelot Massif and at its border with the Dinant Synclino-
rium. The Malmedy Graben, filled with Permian conglom-
erates, developed in the center of the Stavelot Massif with
an SW-NE orientation is separating the Stavelot Massif into
two parts (Barros et al., 2021).

**2.2 CO₂-rich groundwater springs in the Rhenish Massif and the Eifel volcanic fields**

Numerous occurrences of CO₂-rich groundwater springs are recorded in the Rhenish Massif, as shown in Fig. 2. CO₂-rich groundwaters that are present in Belgium are cold waters (12°C on average). They are fed by the recharge from local or regional precipitation, as confirmed by ¹⁸O and ²H isotopic measurements (Barros et al., 2021). The system is dominated by metamorphized sedimentary rocks, and the aquifer zones lie in the first hundreds of meters of these deposits, in the fractured and weathered parts.

Naturally sparkling groundwaters are bottled by three
different companies in the area (Bru-Chevron and Spa
Monopole in Belgium and Gerolsteiner Brunnen in Ger-
many). While it has been proven several years ago that the
dissolved CO₂ present in the springs of western Germany
was the result of mantle degassing, this was not confirmed for
the Belgian springs yet (May et al., 1996; Aeschbach-Hertig
et al., 1996; Barros et al., 2021).

However, the proximity between the two gas-rich ground-
water areas, their common geological context and the scale
of the Eifel volcanism suggest that the gas they contain might
have a common origin. In the Eifel area, volcanic activ-
ity has occurred during the Quaternary. It is thought that
the lower-mantle upwelling under central Europe may feed
smaller upper-mantle plumes (Goes et al., 1999). Indeed,
small plume structures have been identified below the Eifel
area by teleseismic tomography (see Ritter et al., 2001), and
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Simplified geology of the Ardennes Allochthon and its localization in western Europe and in the Rhenish Massif, modified after Fielitz and Mansy (1999).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Occurrences of CO₂-rich groundwater springs in the Rhenish Massif, together with the main bottling sites. The spring locations were compiled from May et al. (1996), Bräur et al. (2013) and internal data from Spadel. The projection of the Eifel volcanic area is depicted after Bräur et al. (2013).</figcaption></figure>
a recent study by Kreemer et al. (2020) has shown that this plume was still buoyant (Kreemer et al., 2020). However, recent studies on the isotopic composition of noble gas in the Eifel area and in Massif Central (France) tend to confirm a direct contribution from the upper mantle (mid-ocean-ridge basalt, MORB). The exact origin of the Eifel volcanism is thus still a matter of debate (Bekaert et al., 2019; Moreira et al., 2018).

## 3 Sampling and analysis

For this study, water samples from both natural springs and wells were collected at 30 different locations distributed within an area (30 km × 20 km) belonging to the Spa, Bru or Malmedy regions (Table 1, Fig. 3). Each of them has been analyzed for major and trace elements, physicochemical parameters, dissolved gases (O₂, CO₂, He, Ne), and carbon isotopes (δ¹³C). During a second campaign, 13 of these 30 sites were resampled, with the specific purpose of analyzing their ³He/⁴He and ⁴He/²⁰Ne isotopic ratios. Samples were stored in copper tubes clamped on both sides to prevent any degassing or air contamination.

Major and trace elements, together with physicochemical parameters and dissolved O₂ and CO₂, were measured at the Spadel hydrochemistry laboratory, following the specific certified procedures: ISO 10523 (pH), ISO 7888 (electrical conductivity), ISO 10304-1 (chlorides, sulfates and nitrates), ISO 9963 (bicarbonates), ISO 17289 (dissolved oxygen) and ISO 17294 (ICP-MS<fnref n="2" /> for the rest of the elements). Dissolved CO₂ was quantified by mineral sequestration via a barium sulfate saturated solution, followed by an inverse titration. Carbon isotopic ratios were measured by a private lab using isotope-ratio mass spectrometry (IRMS) and gas chromatography isotope-ratio mass spectrometry (GC-IRMS). Finally, the dissolved concentrations of helium isotopes were analyzed at the CRPG of Nancy (CNRS-UMR 7358) by static vacuum mass spectrometry after vacuum extraction and purification according to Zimmermann and Bekaert (2020) and Zimmermann et al. (2015). The ⁴He/²⁰Ne ratio was also measured in the same water aliquots, with a quadrupole installed on the extraction line. Instrument sensitivity was determined against a gas standard having an atmospheric composition. Tap water that had been placed in equilibrium with the atmosphere was also analyzed for comparison (Table 1).

Figure 3 shows the location of each sample in the local geology.

## 4 Results

Table 1 presents the gas results measured in our 30 new samples. Data from Victoriaquelle (VQ) and Schwefelquelle (SQ) are from Marty et al. (2020). Data from the Laacher See

Mofetta (LaS) and Wehr 10 well (W10) (Bräur et al., 2013) are also included (i.e., a *mofetta* is a fumarole discharging mostly carbon dioxide). All these values can be compared to the mid-ocean-ridge basalt (MORB) endmember, which is commonly accepted to represent the signature of the upper mantle and has been studied in detail by Graham (2002).

All the data obtained from these groundwater samples of Spa, Bru and Malmedy (cations, anions, dissolved CO₂, ³He/⁴He and ⁴He/²⁰Ne ratios) are available in Table 1.

## 5 Discussion

### 5.1 Chemical composition of groundwater and dissolved gases

The Piper diagram presented in Fig. 4 shows the relative proportion of cations and anions for each sample. The figure shows a calcium–magnesium–bicarbonate type for most samples. Regarding the cations (Ca, Mg, Na + K), the majority of the samples have rather closed compositions and display a globally balanced composition, with the less abundant cation representing in any case at least 10 % of the total. CGB samples are composed of a smaller proportion of Na + K with a bigger proportion of Mg compared to CGS samples. One exception to this balanced composition in cations is the *PdI* sample, which is computed on a near-pure calcium pole, probably because it was sampled in a zone influenced by the carbonate-rich conglomerate of the Malmedy Graben (shown in Fig. 3). The distinction between CG and NCG groundwater samples is much clearer looking for their anion compositions. There is a significant relative enrichment of CG samples in bicarbonates, due to the presence of dissolved CO₂. Only the samples *Art* and *PlG* are out of the range of this carbonate-rich cluster. However, this lower carbonate content is only relative as the sample *Art* has a much higher concentration in sulfates (16.7 mg L⁻¹, whereas most samples contain < 5 mg L⁻¹), presumably from a local geological origin, and the sample *PlG* has a higher chloride content (58.2 mg L⁻¹, whereas the average composition of the other samples is around 5 mg L⁻¹). The latter is presumably from an anthropic origin.

This figure indicates that both CO₂-rich and non-CO₂-rich groundwaters have initially similar compositions, mainly driven by the local lithologies. The enrichment of CO₂ and bicarbonates in groundwater leads to an acidification of groundwater. Hence, CO₂-rich groundwaters are generally more corrosive and more mineralized. However, the waters of this dataset keep the same relative proportions in cation composition, whatever their carbonate contents (Fig. 4). This observation is not really in line with the hypothesis that carbonate dissolution explains the main origin of this dissolved CO₂. As discussed in Barros et al. (2021), this hypothesis was until now the most commonly accepted one for the CO₂ origin. Indeed, if the presence of dissolved CO₂ in carbo-

<footnote n="2">²Inductively coupled plasma mass spectrometry</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Location of each sample with regard to the main geological features.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Piper diagram representing the relative proportion of cations and anions for each sample. CGB = carbogaseous water from the Bru area, CGS = carbogaseous water from the Spa area, NCGB = non-carbogaseous water from the Bru area, NCGS = non-carbogaseous water from the Spa area, and CGM = carbogaseous water from the Malmedy area.</figcaption></figure>

gaseous waters was due to carbonate dissolution, it would be expected that this groundwater would also be enriched in Ca²⁺ and possibly Mg²⁺ ions in comparison to non-carbogaseous waters, but this is not the case in these Spa-Bru waters.

## 5.2 Constraining gas origin with helium and carbon isotopes

Inorganic carbon isotopes have proven to be a very powerful tool to make the distinction between different carbon sources. This isotopic proxy is particularly adapted to sparkling mineral waters (Fillimonova et al., 2020; Carreira et al., 2014; Redondo and Yelamos, 2004). However, the contribution of each carbon source is sometimes difficult to deconvolve. For example, a bulk C composition resulting from the mixing between marine limestones (δ¹³C ≈ 0‰) and organic sediments (δ¹³C ≈ -20‰) may have a δ¹³C similar to that of mid-ocean-ridge basalt (MORB) (δ¹³C ≈ -6.5 ± 2.5‰), as shown by Sano and Marty (1995). Atmospheric δ¹³C is estimated to be around -8‰ (Karolyte et al., 2019).

The δ¹³C values of inorganic carbon measured in our samples range between -8.8‰ VPDB and -1.8‰ VPDB for CG groundwater samples and between -27‰ VPDB and -15.2‰ VPDB for NCG groundwater samples. The distinction between CG and NCG groundwaters is clear: NCG groundwaters have much lighter carbon isotopic ratios. The majority of δ¹³C values in CO₂-rich groundwater are clus-
<figure><img src="image_8.png" /></figure>

**Figure 5.** CO₂/³He ratios vs. δ¹³C values for groundwater samples in relation to mixing between the mantle, carbonate and organic CO₂ endmembers based on Sano and Marty (1995). CGB = carbogaseous water from the Bru area, CGS = carbogaseous water from the Spa area, CGE = carbogaseous water from the Eifel area, NCGB = non-carbogaseous water from the Bru area, and NCGS = non-carbogaseous water from the Spa area. Gas Eifel composition from Marty et al. (2020). VQ, SQ and MORB are depicted from the values presented in Marty et al. (2020). All CO₂-rich groundwater samples fall in or close to the mantle range. Non-CO₂-rich groundwater samples present much lighter C isotopic composition. The impact of calcite precipitation at three temperature (50, 100 and 150°C) is modeled assuming an open system Rayleigh distillation, following Ray (2009) and Barry et al. (2020).

tered around the mantle endmember, but some of them may
also be compatible with a limestone origin (Fig. 5).

The combination of δ¹³C with the CO₂/³He ratio permits
us to make the distinction between the three main sources of
dissolved carbon (mantle, limestone, organic carbon; Sano
and Marty, 1995), as shown in Fig. 5: with the production of
³He in the crust being negligible (Andrews and Kay, 1982),
a high ³He content (and hence a low CO₂/³He ratio) is the
signature of a mantle carbon input. It can be observed in
Fig. 5 that CO₂-rich groundwater samples have CO₂/³He
and δ¹³C values that are compatible with (or very close to)
the mantle MORB endmember. These compositions are quite
close to those of the volcanic gas of the Eifel area (Marty
et al., 2020) (SQ, VQ samples). Although the Spa–Bru CG
groundwaters have δ¹³C values quite close to a limestone
source, their CO₂/³He ratio is 2 to 4 orders of magnitude
lower than the limestone endmember. Their clear enrichment
in ³He demonstrates that the gas dissolved in the Spa–Bru
groundwaters has a mantle origin, probably from the nearby
Eifel volcanic fields.

Although most of the CG waters are close to the mantle
endmember in this CO₂/³He vs. δ¹³C signature, some data
points do not perfectly match the mixing lines between the

three endmembers. Other processes that may have fraction-
ated the initial CO₂/³He and δ¹³C signatures thus need to be
considered (e.g., Ray, 2009; Barry et al., 2020). With CO₂
and helium having different solubilities, partial degassing
may fractionate the CO₂/³He ratio. However, with the sol-
ubility of CO₂ being larger than the one of helium, this pro-
cess should lead to an increase in the CO₂/³He ratio of wa-
ters affected by degassing ratio, which is contrary to what
we observe here. Moreover, there is almost no correlation
between the CO₂/³He ratios and the dissolved helium con-
centrations (after correction for atmospheric helium), which
is an observation that makes this process unlikely (Table 1).
Another physical process that has the ability to modify the
initial CO₂/³He-δ¹³C signature is the precipitation of calcite,
leading to lower the δ¹³C and the CO₂/³He by a Rayleigh
distillation (in open system). In Fig. 5, following Ray (2009)
and Barry et al. (2020), we modeled the effect of calcite pre-
cipitation at various temperatures (50, 100 and 150 °C), as-
suming an initial gas composition similar to the Eifel end-
member (Marty et al., 2020). Although this process may, in
theory, explain part of the scatter observed within the carbo-
gaseous sources of Spa and Bru, we should, however, be cau-
tious, because calcite precipitation has never been observed
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Proportion of crustal and mantle He in the samples, based on R/Ra and ⁴He/²⁰Ne values. CGB = carbogaseous water from the Bru area, CGS = carbogaseous water from the Spa area, CGE = carbogaseous water from the Eifel area, NCGB = non-carbogaseous water from the Bru area, and NCGS = non-carbogaseous water from the Spa area. Mixing lines are computed from Eqs. (3) and (2). Values from Eifel and MORB are depicted from Marty et al. (2020). He atmospheric contribution is negligible. He crustal contribution appears to be more important for CG samples from the Spa area, as the local lithology is richer in uranium minerals.</figcaption></figure>

in any of the many boreholes that have been installed in the area for water production, i.e., in the underlying rocks where these fluids are supposed to have transited. It is nevertheless important to stress that this fractionation does not hamper the identification of a clear mantle signature.

The low δ¹³C and CO₂/³He values of the non-carbogaseous sources stand below the pure mixing curve between the mantle and organic carbon endmembers (Fig. 5). Our geochemical dataset (Table 1) shows that these gas-poor fluids are non-carbogaseous groundwaters that have never been enriched with gas, explaining their very different composition and depletion of δ¹³C and CO₂.

## 5.3 Discriminating He origin with He and Ne isotopes

He isotopic ratios are normalized against the atmospheric isotopic composition and are expressed in R_a (considering atmospheric ³He/⁴He = R_a = 1.382 × 10⁻⁶ (Sano and Fischer, 2013)). In these waters, R_a ranges between 0.92 and 2.70 (±0.02). ⁴He/²⁰Ne ratios are very variable with values ranging from 1.7 to 348.7 (±4%). Based on the measured ⁴He/²⁰Ne ratios and assuming all ²⁰Ne of atmospheric origin, the contribution of atmospheric He can be computed using mixing Eq. (1), where ⁴He/²⁰Ne_a is considered equal to 0.267 according to Holocher et al. (2001), and

⁴He/²⁰Ne_mantle is equal to 1000 according to Dunai and Baur (1995).

$$ \%He_{atm} = \frac{1 - \frac{^{4}He/^{20}Ne_{mantle}}{^{4}He/^{20}Ne_{sample}}}{1 - \frac{^{4}He/^{20}Ne_{mantle}}{^{4}He/^{20}Ne_{air}}} \quad (1) $$

This shows that atmospheric helium is negligible for all CG samples (less than 0.5%) but more important for NCG samples (between 5% and 15%). The measured ³He/⁴He ratios allow for the computation of the origin of helium, which appears to be a mixture of crustal and mantle helium. The proportion of each source (crust and mantle) may be computed according to mixing Eqs. (2) and (3), where ³He/⁴He_crust and ³He/⁴He_mantle are taken to be equal to 0.02 and 6.5 Ra, respectively, according to Dunai and Baur (1995). We chose this 6.5 Ra “mantle” endmember as representative of the subcontinental lithospheric mantle (Dunai and Baur, 1995).

$$ \%He_{mantle} = \frac{(1 - \%He_{atm})(^{3}He/^{4}He_{crust} - ^{3}He/^{4}He_{sample})}{^{3}He/^{4}He_{crust} - ^{3}He/^{4}He_{mantle}} \quad (2) $$

$$ \%He_{crust} = 1 - \%He_{mantle}. \quad (3) $$

Crust-mantle mixing lines are also displayed in Fig. 6.

Helium present in groundwater samples appears to be between 50% and 80% from crustal origin. A distinction can
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Updated conceptual model for the existence of CG and NCG groundwaters. Aquifer zones are located in the weathered parts of the bedrock. The aquifer in fractured, and slate beds isolate aquifer zones from each other. Degassing CO₂ from the mantle is brought to the system through regional thrust faults, and dissolved in groundwater.</figcaption></figure>

be made between CG samples from Spa and Bru area, samples from Spa displaying more crustal helium (from 66% to 81%) than those of Bru (ranging from 54% to 57%). This may result from the local lithologies: Cambrian and Ordovician rocks from the Stavelot Massif are known for their high uranium content, often leading to high radon concentrations in cellars and underground buildings and also ⁴He production through alpha decay (Depret et al., 2021; Vanderschueren, 2011). CG samples of Spa also differ from Eifel samples since they bear a larger proportion of crustal helium than the groundwaters of this volcanic region, where more than 80% of He is from the mantle.

## 5.4 Conceptual model of groundwater and gas circulation

CO₂/³He and δ¹³C measured in the samples from eastern Belgian springs show that the dissolved CO₂ present in the springs from eastern Belgium originates from a mantle contribution. The similar cation proportions for both CG and NCG groundwaters support the hypothesis that both groundwater types are initially the same (i.e., having the same origin). This is at odds with the commonly accepted hypothesis according to which CG groundwater would travel at several kilometers depth to meet and dissolve carbonated layers,

where it could have acquired its dissolved CO₂ content before flowing upward rapidly to the surface. This hypothesis is also whittled by the fact that CG groundwaters are cold groundwaters, having similar temperatures to NCG groundwaters (between 12 and 13 °C). Assuming a vertical temperature gradient of 30 °C km⁻¹, this suggests that water circulation does not occur deeper than a few hundred meters below the surface. ³He/⁴He ratios indicate that crustal ⁴He enrichment occurred during the water circulation in the aquifer, but such addition of crustal fluid is impossible in less than 1 km flow water paths.

The conceptual circulation model of water and gas at depth has then to be updated and is shown in Fig. 7. The conceptual model is very similar to one of the models described by Pisolkar (2017) in his work aiming at developing integrated hydrogeological models representing the classical properties of CO₂-rich mineral water systems in different contexts.

Considering all our data with those by Barros et al. (2021) (notably the isotopic δ¹⁸O and D/H composition of the water), we propose here a revised model (Fig. 7) to explain the origin of the gas-rich groundwaters of the Spa-Bru massif, involving the input of CO₂-³He-rich gas from the nearby Eifel massif through deep crustal faults.

The geomorphology of the system is mostly controlled by basin structures and anticlines. Slate beds act as low per-
meability barriers partitioning the aquifer and isolating the different compartments. Small faults and surface weathering enable the infiltration and storage of water from the surface to the underground, whereas major deep-rooted faults act as CO₂ transport pathway from degassing mantle that can be located several tens of kilometers away from the local aquifers. Where major faults reaching the surface do not encounter a sufficiently permeable and water-saturated zone (i.e., in the slates covered by a clayey colluvium as a result of weathering processes), CO₂-³He-rich dry gases are discharged to the surface and mofettes are observed. The discharge points are mainly springs (at low topographic points or at geological low permeability/high permeability contacts) in the hill slopes or water abstraction wells.

The CO₂/³He and δ¹³C composition of the samples is very close to samples taken in the Eifel volcanic fields (Fig. 5, Table 1); considering the proximity of the two sites (distance less than 100 km), it seems very likely that the gas found in the Belgian springs comes from the degassing of the mantle in the Eifel volcanic fields. In the future, the link between the two regions needs to be further explained in terms of structural geology.

The main structural features existing in the area are the European thrust fault, the Xhoris thrust fault and the normal faults linked to the opening of the Malmedy Graben (see Fig. 1). All these faults could be connected at several kilometers depth to the Midi-Eifel thrust fault, which is a major thrust fault corresponding to the northernmost front of the Variscian orogeny and connecting both regions. This fault acts as an important seismic reflector and could thus be observed below the Stavelot Massif thanks to seismic measurement surveys lead by the DEKORP research group in the early 90s (Stiller et al., 1987).

The exact origin of the Eifel volcanism, between a plume or an upper-mantle contribution, is still a matter of debate. Unfortunately, although the results confirm the magmatic origin of the gas dissolved in groundwater, they do not provide new information to support one or the other hypothesis.

## 6 Conclusions and future research

This paper answers a long-standing question regarding the origin of the dissolved CO₂ in the naturally sparkling mineral waters of eastern Belgium. The combination of δ¹³C and ³He isotopes have shown with a high level of confidence that the dissolved CO₂ in groundwater from the springs and in boreholes was from mantle origin, and it can be very likely attributed to the degassing of the mantle in the neighboring Eifel volcanic fields, located at a distance less than 100 km eastward. The role of the deep-rooted fractures that act as CO₂ transport pathway to the surface are still to be clarified, but several major thrust faults exist in the Rhenish Massif and could have connected the Eifel area with the studied area.

**Code availability.** No specific code was used for this publication. All data and equations were implemented in Microsoft Excel.

**Data availability.** All the data that made it possible to draw the conclusions of this article are presented in Table 1. The company's internal hydrochemical analysis data on mineral waters cannot be shown for reasons of confidentiality. When data from other publications are used, they are accessible in the publication mentioned and included in Table 1.

**Author contributions.** AD is the main contributor of this paper and worked on the conceptualization, the methodology, the field investigation, the writing of the paper and the visualization of the data. P-HB worked consequently on the validation of the methodology and data, the field investigation, the writing, editing and proof reading of the work, and the funding acquisition. LZ was in charge of the noble gas analysis. PJ, AC and FN supervised the work and helped with access to internal data for wells and springs where samples were taken. AD is the main supervisor of the work and consequently also worked on the project administration, paper proofreading and editing, and funding acquisition.

**Competing interests.** The contact author has declared that neither they nor their co-authors have any competing interests.

**Disclaimer.** Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Acknowledgements.** This work is part of the ROSEAU project, as part of the Walloon program "Doctorat en Entreprise", co-funded by the SPW Région Wallonne of Belgium and the company Bru-Chevron S.A. (Spadel S.A.), under grant number 7984. The collaboration with the CRPG of Nancy was made possible thanks to the Europlanet Transnational Access Program, as this project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement no. 871149. The authors acknowledge these two institutions for their support.

**Financial support.** This research has been supported by the Service Public de Wallonie (grant no. 7984) and the Horizon 2020 (EPN-2024-RI (grant no. 871149)).

**Review statement.** This paper was edited by Brian Berkowitz and reviewed by two anonymous referees.

## References

* Aeschbach-Hertig, W., Kipfer, R., Hofer, M., Imboden, D., Wieler, R., and Signer, P.: Quantification of gas fluxes from sub-
continental mantle: the example of Laacher See, a maar lake in Germany, Geochim. Cosmochim. A., 60, 31–41, https://doi.org/10.1016/0016-7037(95)00370-3, 1996.

Agnew, R. J.: Why springs bubble: A framework for gas discharge in groundwater, Groundwater, 56, 859–870, https://doi.org/10.1111/gwat.12789, 2018.

Andrews, J. and Kay, R.: Natural production of tritium in permeable rocks, Nature, 298, 361–363, https://doi.org/10.1038/298361a0, 1982.

Barros, R., Defourny, A., Collignon, A., Jobé, P., Dassargues, A., Piessens, K., and Welkenhuysen, K.: A review of the geology and origin of CO₂ in mineral water springs in east Belgium, Geol. Belg., 24, 17–31, https://doi.org/10.20341/gb.2020.023, 2021.

Barry, P. H., Negrete-Aranda, R., Spelz, R. M., Seltzer, A. M., Bekaert, D. V., Virrueta, C., and Kulongoski, J. T.: Volatile sources, sinks and pathways: A helium-carbon isotope study of Baja California fluids and gases, Chem. Geol., 550, 119722, https://doi.org/10.1016/j.chemgeo.2020.119722, 2020.

Bekaert, D., Broadley, M., Caracausi, A., and Marty, B.: Novel insights into the degassing history of the Earth's mantle from high precision noble gas analysis of magmatic gas, Earth Planet. Sci. Lett., 525, 115766, https://doi.org/10.1016/j.epsl.2019.115766, 2019.

Bräur, K., Kämpf, H., Niedermann, S., and Strauch, G.: Inferences for the existence of different magmatic reservoirs beneath the Eifel area (Germany): A multi-isotope (C,N,He,Ne,Ar) approach, Chem. Geol., 56, 193–208, https://doi.org/10.1016/j.chemgeo.2013.08.013, 2013.

Carreira, P., Marques, J., Carvalho, M., Nunes, D., and da Silva, M. A.: Carbon isotopes and geochemical processes in CO₂-rich cold mineral water, N_Portugal, Environ. Earth Sci., 71, 2941–2953, 2014.

Depret, M., Bruni, Y., Dassargues, A., Defourny, A., Marion, J.-M., Vanderschueren, W., and Hatert, F.: Mineralogical and hydrogeological study of “pouhons” in the lower Paleozoic formations of the Stavelot-Venn Massif, Geologica Belg., 24, 109–124, 2021.

Dunai, T. and Baur, H.: Helium, neon, and argon systematics of the European subcontinental mantle: Implications for its geochemical evolution, Chechim. Cosmochim. A., 59, 2767–2783, https://doi.org/10.1016/0016-7037(95)00172-V, 1995.

Fielitz, W. and Mansy, J.-L.: Pre- and synorogenic burial metamorphism in the Ardenne and neighbouring areas (Rhenohercynian zone, central European Variscides, Tectonophysics, 309, 227–256, https://doi.org/10.1016/S0040-1951(99)00141-9, 1999.

Fillimonova, E., Lavrushin, V., Kharitonova, N., Sartykov, A., Maximova, E., Baranovskaya, E., Korzun, A., Maslov, A., and Baidariko, E.: Hydrogeology and hydrochemistry of mineral sparkling groundwater within Essentuki area (Caucasian mineral water region), Environ. Earth Sci., 79, 15, https://doi.org/10.1007/s12665-019-8721-2, 2020.

Gilfillan, S., Györe, D., Flude, S., Johnson, G., Bond, C., Hicks, N., Lister, R., jones, D., Kremer, Y., Haszeldine, R., and Stuart, F.: Noble gases confirm plume-related mantle degassing beneath Southern Africa, Nat. Commun., 10, 5028, https://doi.org/10.1038/s41467-019-12944-6, 2019.

Goes, S., Spakman, W., and Bijwaard, H.: A lower mantle source for central European volcanism, Science, 286, 1928–1931, 1999.

Graham, D. W.: Noble gas isotope geochemistry of mid-ocean ridge and ocean island basalts: Characterization of mantle source reservoirs, 2002.

Holocher, J., Matta, V., Aeschbach-Hertig, W., Beyerle, U., and Hofer, M.: Noble gas and major element constraints on the water dynamics in an alpine floodplain, Groundwater, 39, 841–852, 2001.

Vanderschueren, H. W.: Le radon dans l'air, dans l'eau et dans les roches. Mesure dynamique de son exhalation et contribution à la caractérisation géologique des matériaux, Geological survey of Belgium professional paper, 308, 96 pp., http://biblio.naturalsciences.be/rbins-publications/professional-papers-of-the-geological-survey-of-belgium/bibliographic-references/vanderschueren_radon_2011 (last access: 18 May 2022), 2011.

Karolyte, R., Johnson, G., Györe, D., Serno, S., Flude, S., Stuart, F., Chivas, A., Boyce, A., and Gilfillan, S.: Tracing the migration of mantle CO₂ in gas fields and mineral water springs in south-east Australia using noble gas and stable isotopes, Geochim. Cosmochim. A., 259, 109–128, https://doi.org/10.1016/J.gca.2019.06.002, 2019.

Kreemer, C., Blewitt, G., and Davis, P. M.: Geodetic evidence for a buoyant mantle plume beneath the Eifel volcanic area, NW Europe, Geophys. J. Int., 222, 1316–1322, https://doi.org/10.1093/gji/ggaa227, 2020.

Marty, B., Almayrac, M., Barry, P., Bekaert, D., Broadley, M., Byrne, D., Ballentine, C., and Caracausi, A.: An evaluation of the C / N ratio of the mantle from natural CO₂-rich gas analysis : Geochemical and cosmochemical implications, Earth Planet. Sci. Lett., 551, 116574, https://doi.org/10.1016/J.epsl.2020.116574, 2020.

May, F., Hoernes, S., and Neugebauer, H.: Genesis and sitribution of mineral waters as a consequence of recent lithospheric dynamics: the Rhenish Massif, Central Europe, Geol. Rundsch., 85, 782–799, https://doi.org/10.1007/s005310050112, 1996.

Moreira, M., Rouchon, V., Muller, E., and Noirez, S.: The xenon isotopic signature of the mantle beneath Massif Central, Gechem. Perspect., 6, 28–32, 2018.

Pisolkar, M.: Why CO₂-rich mineral waters present peculiar signatures, when compared with other types of mineral waters? Vidago-Pedras Salgadas case study, Master Thesis, Tecnico Lisboa, Lisboa, Portugal, https://fenix.tecnico.ulisboa.pt/downloadFile/1689244997257994/dissertacao.pdf (last access: 9 May 2022), 2017.

Ray, J. S.: Carbon isotopic variations in fluid-deposited graphite: evidence for multicomponent Rayleigh isotopic fractionation, Int. Geol. Rev., 51, 45–57, https://doi.org/10.1080/00206810802625057, 2009.

Redondo, R. and Yelamos, J. G.: Determination of CO₂ origin (natural or industrial) in sparkling bottled waters by ¹³C/¹²C isotope ratio analysis, Food Chem., 92, 507–514, 2004.

Ritter, J., Jordan, M., Christensen, U., and Achauer, U.: A mantle plume below the Eifel volcanic fields, Germany, Earth Planet. Sci. Lett., 186, 7–14, https://doi.org/10.1016/S0012-821X(01)00226-6, 2001.

Sano, Y. and Fischer, T.: The analysis and interpretation of noble gases in modern hydrothermal systems, in: The Noble Gases as Geochemical Tracers, edited by: Burnard, P., 249–317, Springer, https://doi.org/10.1007/978-3-642-28836-4_10, 2013.
Sano, Y. and Marty, B.: Origin of carbon in fumarolic gas from island arcs, Chem. Geol., 119, 265–274, 1995.

Stiller, M., Kaerger, L., Agafonova, T., Krawczyk, C., Oncken, O., Weber, M., Former DEKORP Project Leaders, Former DEKORP/BELCORP Research Group, Former DEKORP Processing Centre: Deep seismic reflection profile DEKORP 1987-1A across the western Rhenish Massif, West Germany/ East Belgium, GFZ Data Services, https://doi.org/10.5880/GFZ.DEKORP-1A.001, 2020.

Vanbrabant, Y., Braun, J., and Jongmans, D.: Models of passive margin inversion: implications for the Rhenohercynian fold-and-thrust belt, Belgium and Germany, Earth Planet. Sci. Lett., 202, 15–29, https://doi.org/10.1016/S0012-821X(02)00751-3, 2002.

Zimmermann, L. and Bekaert, D.: Analyse des gaz rares par spectrométrie de masse statique – Mesures et applications.J6637 V1(Techniques de l’ingénieur), 2020.

Zimmermann, L., Füri, E., and Burnard, P.: Purification des gaz rares sous ultravide – Méthodes de purification. Techniques de l’Ingénieur, J6 635, 2015.