{'fraction_non_alphanumeric': 0.03965914915427986, 'fraction_numerical': 0.022648641722193746, 'mean_word_length': 4.674786402472278, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mineralogical and hydrogeological study of “pouhons” in the lower Palaeozoic formations of the Stavelot-Venn Massif, Belgium

Numerous naturally CO2-rich mineral water springs, locally called pouhons, occur in the Stavelot-Venn Massif. These water springs show a particular composition with a high content of iron, manganese and lithium, and are characterised by a redorange colour resulting from iron hydroxide precipitation near the land surface. Radon measurements have shown that these ferruginous deposits are weakly radioactive. The Upper Cambrian black shales of the La Gleize Formation are also known to display radioactive anomalies. These rocks show enrichment in HFSE (Pb, U, Y, Ce, Zr, Ti, Nb) and are depleted in transition metals (Co, Ni, Cu, Zn). Specific minerals such as florencite-(Ce), monazite-(Ce), xenotime-(Y) and zircon have been identified and are probably at the origin of the radioactive anomalies. Uranium was gradually leached from these minerals, transported in solution, and finally concentrated in ferruginous muds. These muds are mainly composed of goethite (most often amorphous), residual quartz and calcite in some samples. The most probable hypothesis is that uranium is adsorbed in small concentrations on the goethite surface. On the other hand, the Ottré Formation (Ordovician) appears to be the main source of lithium, iron and manganese. Pouhon waters have therefore probably leached rocks of various mineralogy and chemical composition during their sub-surface circulation.


Introduction
Numerous naturally ferruginous and CO 2 -rich mineral water springs, locally called pouhons, occur in southeast Belgium. These water sources are mainly concentrated in the Cambro-Ordovician formations of the Stavelot-Venn Massif (SVM). Some are exploited, such as the mineral and thermal waters of Spa (Marie-Henriette) or the mineral waters of Bru-Chevron. These sources are also known for their curative properties and were first mentioned in journals (Naturalis Historiae) of Pliny the Elder (77-79 AD).
Although these water springs have been known and exploited for a long time, many questions remain as their age, their underground circulation pathway, the origins of the CO 2 and their particular composition. Indeed, hydrogeological interpretations are relatively difficult in the SVM due to its complex geological setting. Pouhons are characterised by a redorange colour resulting from iron oxide and hydroxide precipitation linked to changing redox conditions near the surface. Radon measurements have shown that these ferruginous deposits are weakly radioactive (Vanderschueren, 2011). Uranium, at the origin of radon activity, was probably gradually leached from rocks, transported in solution, and finally concentrated in ferruginous muds. In this paper, we present a combined geological, mineralogical, geochemical, and hydrogeological study of pouhons; these data will improve our global understanding of these ferruginous springs.

Geological setting
The SVM, mainly located in Belgium (Fig. 1), is the largest and the best-studied Cambro-Ordovician inlier of the Ardenne Allochthon. It is part of the Rhenohercynian fold-and-thrust belt (Fielitz & Mansy, 1999;Verniers et al., 2001) and is situated close to the Variscan front complex. The SVM has undergone a complex geodynamic history. It was affected by the Caledonian orogeny during the Ordovician and the Silurian, and by the Variscan orogeny at the end of the Carboniferous (Ziegler et al., 1979;Cocks & Fortey, 2009;Torsvik et al., 2012). The SVM was therefore severely deformed and fractured.
Finally, the Salm Group (lower to middle Ordovician) consists of three formations subdivided into eight members. (1) The Jalhay Formation (Sm1) comprises the Solwaster, Spa, and Lierneux members. This formation is formed of greenish to bluish slates, quartzites, and silty slates (locally called quartzophyllade).

Uranium prospecting in Wallonia
In 1852, André Dumont discovered some crystals of torbernite (Cu(UO 2 ) 2 (PO 4 ) 2 ·12H 2 O) during an excursion in the Vielsalm region and was the first to record the occurrence of uraniumbearing minerals in the Ardennes (Dumont, 1852). However, it was not until after the Second World War that interest emerged in the radioactivity of the Belgian subsoil.
From 1955 to 1962, the Geological Survey of Belgium carried out a preliminary radiometric prospection using a portable scintillometer. Some radioactive anomalies have been identified in the SVM during this campaign. Most of them are located in the black shales of the La Gleize Formation (top of the Revin Group) with a radioactivity 4 to 5 times higher than the regional background (Charlet et al., 1983). Additional analyses by gamma spectrometry carried out by Charlet et al. (1977) indicated a relatively high uranium content (about 20 ppm) and a low Th/U ratio (~0.4) in these black shales.
Then, the Geological Survey of Belgium carried out a carborne prospection to measure uranium, thorium, and potassium contents in Wallonia subsoil. Again, they observed a series of anomalies whose spatial distribution seems to confirm the existence of a more radioactive level at the top of the Revin Group (Charlet et al., 1983).
Finally, in 1994 and 1995, the Geological Surveys of Belgium and Grand Duchy of Luxembourg carried out an airborne campaign. The aircraft followed a north-south trajectory spaced of one kilometre with measurements taken every second, i.e. approximately every 70 metres. Data were interpolated using a statistical kriging method to obtain a continuous spatial visualization (Vanderschueren et al., 2007). Results obtained show that the most uranium-rich areas are located along the Cambro-Ordovician limit, in the upper part of the La Gleize Formation (RV5) and the lower part of the Jalhay Formation (Solwaster Member). Both are composed of black shales which are probably at the origin of the high uranium contents. However, these formations have never been studied from a mineralogical point of view.

Hydrogeological setting
Two types of aquifer are known in the SVM: (1) Shallow aquifers related to differential alteration of the bedrock. The most interesting aquifer areas correspond to the alteration of quartzite whereas clay products resulting from slate alteration form zones of low permeability. (2) Deep aquifer zones distributed in the more quartzitic and fractured/faulted zones of the bedrock which probably fed CO 2 -rich mineral water springs (Monjoie 1995;Jobé, 2014).
Pouhons result from complex underground processes involving sequentially water infiltration, interactions with rocks, CO 2 assimilation, upwelling and finally water emergence (Barros et al., 2021). However, their spatial distribution shows that they are not randomly located but that their location seems to be controlled by geological and geomorphological factors. First of all, pouhons are mainly located in the immediate vicinity of the tectonic structures mapped by F. Geukens in 1999(Monjoie, 1995 Fig. 1). Geophysical profiles combining ERT (electrical resistivity tomography) and IP (induced polarization) methods were carried out above known uplift zones in the Ardennes region (Defourny et al., 2020). These profiles always show important contrasts in both electrical resistivity and chargeability relatively close to the vertical axis of pouhons (Defourny et al., 2020). Resistivity anomalies are interpreted as saturated fractures while chargeability anomalies are probably linked to iron oxides and hydroxides precipitation in the upper part of the aquifer (Defourny et al., 2020). Faults and fractures therefore probably act as preferential drains enabling the uplift of CO 2 -rich groundwater to the surface. Then, pouhons are generally located in valleys (Monjoie, 1995;Fig. 2A-B). Nearly 80% of pouhons are situated closer than 100 m from a river and at an altitude between 300 and 400 m (Fig.  2B). This is probably due to the hydrostatic pressure decrease in the valleys which would favour the release of CO 2 and then the upwelling of water (Monjoie, 1995). Finally, we notice that pouhons located in the SVM are mainly concentrated in rocks of the Revin Group ( Fig. 1).

Hydrochemistry
CO 2 -rich mineral water springs in the SVM generally show calco-magnesian bicarbonate hydrochemical facies (Fig. 3). The average chemical composition of pouhons is given in Table 1. They are naturally oversaturated in CO 2 (up to 4 g/l). These high carbon dioxide contents imply a slightly to moderately acidic pH inducing a higher aggressiveness. However, waters are moderately mineralised as shown by the electrical conductivity. Pouhons show a particular composition with a high content of iron and manganese (Jobé, 2014; Table 1). They also frequently contain unusually high concentrations of lithium (Table 1). On the other hand, they are characterised by low content of nitrate indicating the absence of anthropogenic pollution. Relatively close springs can show significant hydrochemical variations due to the different pathways of waters through the complex network of fractures or due to the aquifer from which they originally came (Monjoie, 1995).
The 18 O-2 H analyses show that these waters are situated on the meteoric waters line, excluding a volcanic origin ( Fig. 4A-B). Isotopic signature suggests that pouhons are primarily fed by meteoric water (Barros et al., 2021). The origin of CO 2 in pouhons is still unknown and seems relatively complex. According to Jobé (2014) and Barros et al. (2021), CO 2 could result from a combination of various processes. Two hypotheses are generally advanced: (1) CO 2 would be produced by dissolution of marine limestones and/or carbonate nodules according to this reaction: CaCO 3 + 2H + → Ca 2+ + CO 2 (g). However, no carbonate rocks were found in the SVM. They may occur at depth below the Midi-Eifel overthrust suggesting that meteoric water should infiltrate at least 2 km (Barros et al., 2021). (2) CO 2 would be of mantellic origin and produced by volcanic degassing related to the neighbouring Eifel area. The extension of the Eifelian mantle plume at depth or the presence of an unknown shallow magmatic reservoir under the SVM may then both constitute a potential source of CO 2 (Barros et al., 2021).
Recent isotope analyses indicate that the ẟ 13 C of CO 2 ranges from -9 to -2‰ and up to ~15% of He from magmatic origin has been identified in dissolved gases (Barros et al., 2021;Fig. 5). Therefore, it seems that CO 2 shows a mixed origin, i.e. both carbonate and magmatic (Barros et al., 2021). The low temperatures of springs (~10 °C) and the presence of CO 2 gas emanation on the surface (locally called "mofettes") suggest that CO 2 probably rises from its source and is added to groundwaters at shallower depths (Barros et al., 2021).

Radioactivity of pouhons
The natural radioactivity of waters is directly related to the nature of rocks that they interact with during their underground circulation pathway. Most natural groundwaters contain radon produced by the disintegration of uranium. Measurements of radon activity in pouhon waters were carried out by Vanderschueren (2011). Most of them show moderate activity (10 to 100 Bq/l) with some exceptions which show higher radiation levels (1220 Bq/l).
As mentioned before, pouhons show an orange-red colour linked to iron oxide and hydroxide precipitation. Iron is initially reduced to Fe 2+ in conditions of pH (between 4 and 6) and Eh (between -0.3 and +0.2 V) prevailing at depth (Defourny et al., 2020). However, conditions become more oxidizing during groundwater upflows inducing rapid oxidation of Fe 2+ to Fe 3+ according to this reaction: 4Fe 2+ + 3O 2 + 2H 2 O = 4FeO(OH). In addition, the release of CO 2 close to the surface induces an increase in pH which can also lead to iron hydroxide formation. Radon measurements have shown that these ferruginous deposits are also weakly radioactive (Vanderschueren, 2011). Uranium at the origin of radon activity was therefore solubilised  from an initial source, transported in solution probably under the form of UO 2 2+ and finally concentrated in ferruginous muds. It is necessary to characterise the mineralogy of these muds to determine which mineral phase carries the high uranium (or thorium) contents.

Data collection and analyses
Various rocks collected in the SVM and almost 200 samples of ferruginous muds analysed by H.W. Vanderschueren to measure their radon activity were available to characterise their mineralogy ( Fig. 1; Fig. 6A-H). All these samples were milled to obtain a homogeneous powder and analysed from a mineralogical and geochemical point of view, with particular attention to black shales of the La Gleize Formation recognised as presenting radioactive anomalies.
Major oxides (SiO 2 , Al 2 O 3 , TiO 2 , Fe 2 O 3tot , MnO, MgO, CaO, Na 2 O, K 2 O, P 2 O 5 , LOI) and some trace elements (Nb, Y, Zr, Sr, Rb, U, Th, Pb, Zn, Ni, Co, Ga, Ba, Cr, Ce) were obtained by X-ray fluorescence (XRF) with an ARL 9400XP spectrometer. Lithium borate fused glass discs were made to measure major elements while trace elements were analysed on pressed powder pellets. Matrix corrections were realised, and accuracy is estimated as better than 1% for major elements and 5% for trace elements.
Mineralogical analyses were performed by powder X-ray diffraction (XRD) using a Bruker D8-ECO diffractometer (CuKα radiation, λ 1.5418 Å). X-rays are filtered through a Ni filter to limit the Kβ contribution. The goniometer covers a 2θ angle between 2 and 70° with a step of 0.02° and a counting time of 1 second per step. X-ray powder patterns were interpreted using the Bruker EVA 3.2 programme with the ICDD database (PDF-2 version). Quantification of mineral phases was performed by Rietveld refinement of the powder diffraction patterns, by using the Bruker TOPAS 4.2 programme. Chemical analyses have not been used to constrain the Rietveld refinements.
Qualitative chemical analyses were performed with a JCM-6000 NeoscopeTM electron microprobe (Mineralogy Laboratory, Natural History Museum, Luxembourg) using an accelerating voltage of 15 kV on polished section and ferruginous muds. EDX detector is an SDD of JEOL technology.
Infrared spectra were performed on ferruginous muds by a Thermo Nicolet NEXUS 470 FTIR spectrometer with a 1 cm -1 resolution, over the 400-4000 cm -1 range. The samples were prepared by mixing 2 mg of a sample with KBr and compressed to form a pellet which was dried for a few hours at 110 °C.  (Barros et al., 2021).

Florencite-(Ce), CeAl 3 (PO 4 ) 2 (OH) 6
Florencite-(Ce) has been identified in various samples, as irregular grains with diameters between 5 and 50 µm. They are generally associated with pyrite although isolated grains have also been observed (Fig. 7C). Semi-quantitative EDX analyses indicate that florencite mainly contains Ce, Al and P, with smaller amounts of La and Nd.

Monazite-(Ce), Ce(PO 4 )
Monazite-(Ce) is relatively abundant and was identified in most samples, as polygranular masses with a size ranging from 50 to 200 µm (Fig. 7D-E). The EDX analyses show that these grains mainly contain Ce and P, with smaller amounts of La and Nd.

Sphalerite, (Zn,Fe)S
Sphalerite grains are very small (<10 µm) and always occur as inclusions in pyrite (Fig. 7F). Their EDX analyses show the presence of Zn and S, with smaller amounts of Fe and Cd.

Xénotime-(Y), YPO 4
One grain of xenotime-(Y) has been identified in the Viel sample, as an isolated grain of very small size (1-2 µm). The EDX analysis shows Y and P as main constituents, thus confirming the identification.

Zircon, Zr(SiO 4 )
Zircon has been frequently identified as isolated crystals of variable size (10 to 80 µm). Some crystals are euhedral and well zoned (Fig. 7G-H), and some of these zonations are enriched in calcium, uranium, thorium, scandium and cerium.
Geochemical data are given in Table 2 and Table 3. The black shales of the La Gleize Formation mainly contain SiO 2 and Al 2 O 3 , as well as K 2 O, MgO and Fe 2 O 3 in lower concentrations. The proportions of TiO 2 , MnO, CaO, Na 2 O and P 2 O 5 generally do not exceed 1 wt.%. Trace element contents are normalised to the upper continental crust (UCC, McLennan, 2001) and plotted in a semi-logarithmic diagram (Fig. 8A). These rocks are mainly enriched in V and high field strength elements (HFSE; Pb, U, Y, Ce, Zr, Ti, Nb) except in Th and Cr, and depleted in transition metals (Co, Ni, Cu, Zn) and Sr. The U, V, Ce and Ti enrichments are the most pronounced, particularly in black graphitic shales (Viel).
A thin level strongly affected by weathering was observed and sampled (Spox) within the black shales of Spa (SP2 and SR2; N50.48802, E5.84391). It is mainly composed of goethite, quartz, muscovite, and chlorite. For this sample, major and trace elements are normalised to the protolith to determine the behaviour of these elements during weathering (Fig. 8B). It shows a significant enrichment in P 2 O 5 and transition metals as well as a depletion in HFSE except for Zr. Copper shows the highest enrichment, with an average content of 20.5 ppm in slates and 366 ppm in this sample (Table 3). Uranium previously enriched in the protolith was not detected in the Spox sample.

Other formations of the SVM
Rock samples collected in the other formations have only been studied by X-ray powder diffraction and XRF (Tables 2 and 3). The Hourt, La Venne and Jalhay formations show similar mineralogy and are mainly composed of quartz, muscovite, chlorite, feldspar and rutile ( Table 2). The Ottré Formation shows more various mineralogy ( Table  2). The reddish to purplish colour of these rocks results from the occurrence of hematite (Fe 2 O 3 ). This mineral has been detected in all samples except in the coticule (Cot2). The Lienne valley ore (UXO) displays the highest proportion in hematite (33.45 wt.%), whereas the average percentage of this mineral in the host red slates is around 8 wt.%. Spessartine (Mn 3 Al 2 (SiO 4 ) 3 ) has been detected in the coticule (Cot2), the pseudocoticule (Pseu), the ore from the Lienne valley but also in the slates of Thier del Preu (Pcot1, Pcot2) and Salmchâteau (FCU).  Table 2). However, the Ottré Formation stands out and shows significant enrichments in Fe 2 O 3 and MnO (Fig. 9A). The highest manganese contents are observed in the pseudocoticule (Pseu = 15.65 wt.% MnO), the Lienne valley manganese ore (UXO = 10.11 wt.% MnO) and in the coticule (Cot2 = 5.61 wt.% MnO). The CaO and P 2 O 5 contents also tend to increase in the Ottré Formation.
Trace element contents (Table 3) were once again normalised to the upper continental crust (UCC, McLennan, 2001), and plotted in a semi-logarithmic graph (Fig. 9B). The Hourt Formation is principally depleted in trace elements. The Jalhay and La Venne formations show a geochemical signature similar to that of the La Gleize Formation but are generally less enriched in trace elements, except in Zn, Ce, Co and Ni. The Ottré Formation is mainly enriched in Ba, Cr and Co, and depleted in U and Th. Finally, the Lienne ores, the coticule and the pseudocoticule show particularly high chromium contents of 676 ppm, 1067 ppm and 940 ppm, respectively.
A correlation matrix grouping majors and trace elements has been calculated (Fig. 10). Pearson coefficients show that TiO 2 is relatively well correlated with Al 2 O 3 and K 2 O. CaO seems to be correlated with MnO and more particularly with P 2 O 5 . Strong correlations have been observed between Ga, Nb, Rb, TiO 2 , Al 2 O 3 , MgO and K 2 O. Co and Cr show a high correlation with MnO. Uranium does not seem correlated with other elements, except vanadium (0.82).

Ferruginous muds
Mineralogical analyses have been performed on ferruginous mud samples. X-ray powder patterns are generally marked by the presence of two very large peaks, located around 2θ values of 35° and 63°. This shows the extremely poor crystallinity and probably the amorphous nature of muds. However, some X-ray powder patterns show better crystallinity. Goethite is the most frequently identified mineral, but quartz, calcite and dolomite have also been detected in some samples. Infrared spectra were also performed, showing absorption bands that coincide fairly well with those of goethite (3300, 1800, 1000, 890, 795 and 595, 470 cm -1 ), indicating once again that these muds are mainly composed of this mineral (Fig. 11). Other well-marked bands have been observed, showing that muds contain other phases. The 1495 cm -1 band is certainly produced by vibrations of carbonate groups, while the band at 1630 cm -1 corresponds to the deformation vibrations of the water molecule (Kodama, 1985). The 1365 cm -1 peak is probably associated with the vibrations of nitrate groups (Kodama, 1985).
From a geochemical point of view (Table 4), these muds are mainly composed of Fe 2 O 3 (± 65 wt.%), as well as smaller proportions of SiO 2 (± 3.86 wt.%) and P 2 O 5 (± 3.78 wt.%). Loss on ignition (LOI) is relatively high indicating a significant loss in volatiles. Sampling conditions did not systematically allow the collection of pure ferruginous deposits. Contaminations are noticeable in GER, PIL and SAUV samples with a significant increase of the silica contents compared to Fe 2 O 3 . Regarding trace elements, muds are mainly enriched in Sr, Ba, Ce, Zr, Cu and Zn ( Fig. 12; Table 4). The RUY sample shows the highest zirconium content (1843 ppm), and the BAR sample shows a relatively high zinc level (847 ppm). Uranium and thorium were not detected.
Images obtained with a scanning electron microscope show that grains in the muds are formed by an agglomeration of spheroids (1 µm diameter) comparable to spherulites. An EDX elemental mapping has also been performed, which clearly shows that the muds are mainly composed of iron with some areas enriched in Si, Al and Ca. Some tiny spots of uranium have also been detected.

Origin of uranium in the black shales of the La Gleize Formation
Black shales of the La Gleize Formation show classical mineralogy since they are mainly composed of quartz, muscovite, chlorite, feldspars, with variable quantities of pyrite and rutile (Table 2). More specific minerals, such as florencite-(Ce), monazite-(Ce), xenotime-(Y), and zircon, have also been identified by SEM analyses. It is interesting to note that alluvial monazite-(Ce) was previously observed by Burnotte et al. (1989), in rivers cross-cutting the slates of the Revin Group in the SVM.
The ionic radius of Ce 3+ (1.01 Å; Shannon, 1976) and Y 3+ (0.90 Å) are relatively close to that of U 4+ (0.89 Å). Uranium can consequently be incorporated into the crystallographic site of cerium or yttrium by coupled heterovalent substitutions with calcium at a single site, without affecting the mineral structure (Ce 3+ → 0.5 Ca 2+ + 0.5 U 4+ ). This substitution explains the occurrence of Th and U in many rare-earth minerals such as florencite-(Ce), monazite-(Ce) or xenotime-(Y). U 4+ can also be incorporated in the crystal structure of zircon by homovalent substitution with Zr 4+ at a single site (Meis & Gale, 1998), and rutile may also insert some uranium, probably by homovalent substitutions with Ti 4+ (Mezger et al., 1989). All these minerals may incorporate small quantities of uranium into their structure and are therefore probably at the origin of the radioactive anomalies observed in these black shales. Uranium was likely gradually leached from these minerals, transported in solution, and finally concentrated in ferruginous muds.
From a geochemical point of view, these black shales are mainly enriched in V and HFSE (Pb, U, Y, Ce, Zr, Ti, Nb; Fig.  8A). Uranium, cerium and vanadium display the highest enrichments observed in the SVM, although uranium contents are relatively low and do not exceed 20 ppm. It is noteworthy that Ce and Y contents, reaching around 200 and 40 ppm respectively (Table 3), are similar to those reported by Burnotte et al. (1989) for the black slates of the Revin Group from La Vecquée and Grand-Halleux. On the other hand, these rocks are depleted in transition metals (Co, Ni, Cu, Zn) and in Sr (Fig. 8A).
Such low trace element contents are unusual in these types of rocks. Black shales are generally enriched in numerous redox -sensitive metals (V, Cr, Ni, U, Co Cu, Zn) due to the conditions prevailing in their formation environment (Vine & Tourtelot, 1970). The enrichment was probably inhibited in our samples by the unavailability of metals in the water column or by a high sedimentation rate inducing a high clastic dilution. The black graphitic shales of Vielsalm (Viel) show the most significant enrichments. Trace element contents in these rocks are therefore probably controlled by the organic matter and/or by the clay material contents.

Geochemistry of the Ottré Formation and correlations
The Ottré Formation is characterised by a significant enrichment in Fe 2 O 3 , MnO, Cr and Co, compared to the other rocks of the SVM ( Fig. 9 A-B). The Fe 2 O 3 enrichment results from the occurrence of hematite, whereas the MnO enrichment is directly related to the occurrence of spessartine and rhodochrosite. These geochemical features are well known in the SVM, and were already described in several papers dealing with the geochemistry of coticules and pseudocoticules (Krosse & Schreyer, 1993;Baijot et al. 2011;Herbosch et al., 2016). Uranium was not detected by XRF in these rocks.
Some correlations have been observed between major and trace elements (Fig. 10). Titanium and aluminium are both highly immobile, and the strong correlation between these elements probably results from their similar geochemical behaviour. The correlation between calcium and phosphorus is certainly due to the presence of accessory fluorapatite in these rocks (Hatert et al., 2002), as well as to the incorporation of calcium in the structure of cerium phosphate minerals. It is noteworthy that monazite-(Ce), florencite-(Ce), and xenotime-(Y) were yet known to occur in the slates of the Ottré Formation at Salmchâteau (Hatert et al., 2002;Blondieau et al., 2017). The correlation between TiO 2 and Nb can be explained by the ability of rutile to incorporate niobium into its crystal structure. Indeed, Nb 5+ may substitute to Ti 4+ in rutile, the excess charge is then compensated by metal vacancies or by a complementary amount of Fe 2+ (Deer et al., 1962;Graham & Morris, 1973). Gallium and rubidium are particularly well correlated with Al 2 O 3 and K 2 O; these elements are therefore probably associated with muscovite. The correlation between chromium, cobalt and manganese is not well understood. Cr could substitute to Al in spessartine thus explaining the strong correlation between Cr and Mn.

Composition and origin of ferruginous muds
Ferruginous muds are mainly composed of goethite (most often amorphous), residual quartz, calcite, and dolomite in some samples. This assemblage suggests that iron initially dissolved in water precipitates essentially as goethite during groundwater upflow. Calcite and dolomite could precipitate due to a rapid drop in CO 2 partial pressure, a rise in ambient temperature and an increase in oxidation potential. Infrared spectra confirm that muds are mainly composed of goethite; however, other peaks, probably corresponding to carbonate, phosphate and nitrate groups have also been identified.
Chemical XRF analyses show that the muds are mainly composed of Fe 2 O 3 (~65 wt.%), with smaller proportions of SiO 2 (± 3.86 wt.%) and P 2 O 5 (± 3.78 wt.%; Table 4). The high LOI observed are probably related to hydroxide and carbonate groups. Trace element analyses show relatively high content of Ba, Ce, Zr, Cu and Zn (with some Sr, V, Co, Y; Table 4; Fig.  12), but no correlations have been established between the geochemical composition of ferruginous muds and the hydrochemical composition of the waters. This suggests that the composition of muds is not directly related to the initial composition of groundwaters. The high zirconium contents observed in muds certainly result from the gradual leaching of zircon releasing small amounts of uranium at the same time. RUY sample shows the highest radon activity and also the highest zirconium content (1843 ppm). There is a potential relationship between the zirconium content and the uranium content in muds although the radon activity does not appear to effectively correlate with the zirconium content in other samples.
Uranium and thorium have not been detected by XRF spectrometry in the mud samples (Table 4). This indicates that the U and Th contents are relatively low. An elementary EDX mapping has made it possible to detect some tiny spots of uranium, which do not seem correlated with other elements. The most likely hypothesis is that uranium is adsorbed in small concentrations on goethite surface. Indeed, goethite is characterised by a surface charge and has a high adsorption capacity for dissolved uranium (Hsi & Langmuir, 1985). These authors demonstrated that the speciation of uranium in solution (dependent on pH) significantly controls its degree of adsorption. The adsorption of U 6+ is optimal in a pH range between 5 and 8.5.

Hydrochemistry of pouhon waters
Pouhon waters are characterised by a particular hydrochemical composition. They are significantly enriched in iron and manganese and have frequently abnormally high lithium concentrations. In the same way, the Ottré Formation is enriched in Fe 2 O 3 and MnO compared to other formations of the SVM. Dekoninck et al. (2019) carried out geochemical analyses on slates of the Ottré Formation. Data indicate that these rocks are also relatively rich in lithium with contents ranging from 50 to 328 ppm. Lithiophorite ((Al,Li)MnO 2 (OH) 2 ) is a weathering mineral abundantly found in this formation (Fransolet, 1979). The absence of this mineral in the other formations of SVM shows that lithium is mainly concentrated in the slates of the Ottré Formation.  Although the composition of groundwaters appears to be closely related to the slates of the Ottré Formation, the extent of these rocks is relatively limited and pouhons are mainly located in rocks of the Revin Group. Waters have therefore interacted and potentially leached these rocks during their upwelling. The black shales of the La Gleize Formation appear to be the main source of uranium at the origin of radon activity observed in ferruginous muds. Moreover, the loss of uranium in the alteration products (Spox) compared to the protolith demonstrate the relatively high mobility of this element during weathering processes. These results suggest that the naturally CO 2 -rich waters occurring in the pouhons of the SVM have probably leached rocks of various mineralogy and chemical composition during their circulation in the sub-surface.

Acknowledgements
The authors thank N. Delmelle, who performed the XRF analyses, as well as J. Vander Auwera who gave access to the XRF instrument. Many thanks are due to S. Philippo for the Table 4. Major oxides (wt.% ) and trace elements (ppm) geochemical XRF analyses of ferruginous mud samples. LOI = Loss on ignition; ND = not detected.