1. Introduction
As part of the quest for a carbon-free energy mix, or at least one that is less CO
2-intensive, H
2, once a chemical product, is now being considered as a means of storing intermittent electricity and decarbonizing certain industrial processes. H
2 is still manufactured, but its presence underground, long known but poorly understood, has been reassessed [
1]. Indeed, H
2 is not only an energy carrier but also a source of energy, and continuous production from a field, discovered by chance in Mali, has confirmed this for over 10 years [
2,
3]. Many countries have adapted their mining legislation to allow prospecting for this new resource, and companies have embarked on exploration [
4,
5].
The understanding of hydrogen systems is still in its infancy although the main generating rocks have begun to be inventoried and mapped [
5,
6]. The main reactions are also fairly well known (see [
6,
7] and reference therein), even though their kinetics remain to be better defined. The main reaction is the oxidoreduction (reduction of water and oxidation of an Fe rich mineral), which may be active in different geological contexts. Radiolysis also leads to the generation of H
2 from water [
8]. Additional H
2 generating rocks are the organic rich facies, such as coal, as the late maturation of the organic matter generates H
2, which may remain as a free gas at high temperatures (above 200 °C) [
9,
10]. When subsurface data exist, the presence and current depth of these generating rocks can be investigated. For hydrocarbon exploration, the presence of reservoirs, traps and seals is verified before drilling. For the H
2 system, the characteristics of the seal and even the existence of perfect seals are still debated. At Bourakebougou, in Mali, the only field currently in production, the seals of the various reservoirs are dolerites and shale [
11,
12]. However, some authors doubt the sealing capacity of regular poor permeability rocks against H
2 dismigration [
13], while others note that salt could be an excellent cap rock [
14]. Thus, the presence of H
2 leaks is being investigated for the time being. These leakages are interpreted as proof of an existing and active H
2 system. The fault zones, as with all fluids, are a preferential pathway for H
2 migration [
15,
16], where the fault valve effect could result in pulses [
17]. In these areas, leakages are usually greatest, and more continuous.
Various tools enable the H
2 explorationist to map H
2 presence in the soil. As H
2 is virtually absent from air (0.5 ppm), the presence of H
2 in the soil is interpreted as evidence of a leakage from a deeper system, generating rock, degassing aquifer, accumulation or migration pathway. As H
2 presence affects vegetation and sometimes soil geometry, remote sensing can be used. This includes satellite images (such as Google Earth), Digital elevation models (DEM), and vegetation indexes often based on infrared acquisition. Sub Circular Depression (SCD) with abnormally high H
2 content in the soil, also called fairy circle, where vegetation disappears or is different, has been described first in Russia [
18] and, then, in the USA [
19], Brazil [
20], Australia [
21,
22] and Namibia [
23]. The total absence of vegetation observed in Brazil [
20] is not universal and just different; less tall, vegetation can also be observed as in Colombia [
24], where sugarcanes continue to grow but plants are smaller in the H
2 leakage zones. Usually, field acquisitions are prepared on the basis of the existence of potential H
2-generating rocks in the areas where vegetation anomalies have been noted [
6]. This methodology is basic but appears to be effective and, very often, the soil in these areas of missing or anomalous vegetation contains an unusual H
2 percentage.
Since [
18,
25], the vast majority of soil gas measurements have been carried out in the same way, with just a few variations [
6,
26]. This method is known to have some limitations [
27,
28], but it is the only one that is easy to implement. In this paper, after describing the usual way of measuring H
2 content in the soil, we present a new tool that has been developed to overcome these limitations. We will discuss what we have learned thanks to it and, especially, the depth-dependent measurement of H
2 content that the tool allows. The developed tool can also be used to sample gases present in the soil. Different containers could be used and the comparison between them will be presented in another publication.
The measurements published in this article all come from the Colombian Andean foreland; see the location in
Figure 1, called the Llanos Basin for the northern part and the Putumayo Basin for the South [
29]. The Llanos Basin is a flat area of about 200,000 km
2 that extends eastward to the Guiana craton. It corresponds to the foreland of the today Eastern Cordillera, a Mesozoic depocenter inverted during the Tertiary [
29,
30]. This basin is rich in hydrocarbons and, therefore, the subsurface geometry is rather well known [
31]. Between the basement and the tertiary deposits, the Paleozoic and the Mesozoic are still present in the western part but pinch out eastward. The Putumayo Basin southward extends to the Ecuadorian border; it is less deep, the oil and gas exploration is, however, active [
32]. The exact H
2 potential of these two basins, and repairing their more prospective areas, is not the subject of this article. For reasons of confidentiality, we will not provide further details about the location of the studied areas or on the generating rocks we are considering.
2. Existing Soil Gas Measurement Methods
Soil H
2 content is measured using sensors that analyze the gas contained in the soil porosity. This gas is largely air—around 20% of O
2 and 80% of N
2—but it can vary if deep sources bring other gases to the near surface, or if microorganisms generate gases in the soil. The typical additional gases are CH
4 and CO
2. The sensors used analyze all these gases. Specifically, they analyze O
2, CH
4 and CO
2, in % and calculate a “balance” that is interpreted as N
2. The sensors chosen for H
2 exploration also measure H
2, in ppm, as this gas is in smaller quantities in the soil. The GA5000 or BIOGAS5000 (manufactured by QED Geotech) used by many authors [
19,
33,
34] may also be equipped to measure H
2S and CO, which is mainly a safety feature when working in volcanic zones, as these two gases are toxic. Other commercial gas analyzers are also performing the same kind of measurements; in our case, we used an Optimax (manufactured by MRU Instruments Inc.), and laboratories such as the CSIRO in Australia are developing their own tools [
27,
33]. Permanent H
2 sensors, called Parhys, have been also developed by Engie [
20].
Microorganisms are always present in the soil and many of them are known to consume H
2 [
35,
36,
37,
38]. These authors have shown that correcting this consumption to calculate a flow rate at depth from surface measurements is not easy, as it depends on many factors, including, for instance, permeability and water content. The consumption has been measured in the laboratory in volume per time. As a result, the speed of the transport of the H
2 is a key parameter; if it is too slow, the microorganisms have the time to consume all the H
2 [
36]. Inverting the signal to deduce from surface flow the existence or depth of a leaking reservoir or degassing aquifer is highly speculative [
6,
35]. To compute a flow from punctual data, even numerous, 50 sensors on a single SCD recorded H
2 for 8 months in Brazil, also required hypotheses [
20]. Even if these H
2 near-surface data sets (80–100 cm depth) cannot be directly inverted, they are key to eliminating unrealistic hypotheses on the deep H
2 flow. For instance, in Albania, the H
2 flow measured in a mine is too high to be directly related to a H
2 generation zone and is interpreted in terms of leakage from a deeper reservoir [
39].
In the field, all the teams who performed soil gas measurements have noticed that H
2 contents are low when measurements are taken very close to the surface (10 to 20 cm), and all measurements are generally taken at around 80/100 cm from the surface (see additional material in [
6]). This value has no scientific basis, but is the result of various constraints, such as the length of the drill bit that you can buy and take with you in your luggage. Furthermore, to penetrate the ground with a drill bit, we have to put weight on the tool, which is impossible with a tube that is too long. In some cases, where the ground is too hard, soil may be missing in some basins; we are unable to drill as deeply and measurements are taken at a shallower depth.
2.1. Tube Installation
How the hole should be made is debated, as H
2 can be generated by perforation, as it is well known on a larger scale in the oil industry [
40]. Some authors suggest avoiding the rotary mode with the drill [
28], while others note that results are generally insensitive to this factor [
6]. From the experience of one of the authors, the key factors are the water content in the soil (after or during a rain for instance) and the temperature of the drill bit, usually related to the drilling time or rotation speed; the H
2, in this case, comes from the oxidation of the drill bit.
Once the hole is made (
Figure 2), it normally remains open if the soil is not too loose, and a cannula is inserted into the ground. This tube, made in copper or stainless steel, is pierced laterally at the end which is in the soil. The drill has usually a diameter of 14 or 16 mm; the tube, which is slightly thinner, is inserted into the vacuum as soon as possible after drilling and the measurement is taken: the sensor pump pumps air which is into the soil and analyzes it.
There are two stages in the procedure: first, drilling, and then, inserting the tube as quickly as possible. This process is clearly handcrafted, but it was the best that could be performed with the available equipment. When permanent monitoring is not possible, an alternative to obtain complementary data is to leave the tube in the ground, close its upper end, and measure the H
2 content several times over a day or two, to ensure that the initial drilling no longer affected the measurement. Field constrains do not always allow this but, usually, a recharge of at least 1 or 2 h long could be tested. By experience, recharge could be visible after just 1 h and permits optimism regarding the fact that the H
2 leakage is active. Long-term monitoring has the advantage of ensuring that the flow is continuously active, as opposed to measuring residual H
2 adsorbed, as in shale minerals in the soil, but such data sets are still rare [
20]. Conversely, you can be unlucky and measure zero at a location where a few hours or days later, the H
2 content of the soil could have been high, as all monitoring data have shown [
20].
The fact that the measurement changes near the surface has been noticed by many authors, but as they also noted that, in SCDs, soil H2 content varies throughout the day and spatially across the structure; they did not attempt to quantify these variations.
2.2. Gas Measurement and Sampling
In terms of commercial gas analyzers, many field campaigns have been conducted with the GA2000, GA5000 or BIOGAS5000 equipped with various sensors. The difference between the three tools is mainly their age and the presence, or absence, of an integrated GPS. They are robust devices; the laboratory tests that many teams have carried out with known blends of gases are good, and in the field, they do not seem to be affected by weather conditions. They need, however, to be protected from the dust and the water and filters have to be installed. H
2S has also to be avoided, but this kind of gas is rather scarce outside of volcanic contexts. When we want to sample soil gas, taking advantage of the sensor pump is useful, even if vacuum tubes have a slight gas-sucking effect. The gas can be sampled either at the inlet or outlet of the sensor, but it has been found that there is a bias at the outlet, it is better to carry this out before the gas circulate in the sensor. So, the various researchers have either fitted a T-junction on the pipe between the rod and the sensor, or inserted a spike directly into the silicone pipe, which means changing it from time to time [
6]. The correlations between the field measurements and the measurements that could be performed in laboratory are usually good [
6,
19] but may slightly differ. The reasons for this shift are, essentially, the gap between the timing of the measurements and the filling of the exetainer. The maximum of H
2 content is usually reached after 40 s, which is mainly due to the sensor response time, although the speed of the gas transport in the soil and up to the analyzer may also play a role. As a result, the best time to sample the gas is unclear. If we wait until the maximum is recorded on the sensor, the H
2 content of the circulating gas is lower. Here too, everyone’s technique for obtaining a representative sample is based on trial and error. Exetainers are rather small and this method of sampling also precludes some analyses.
In our case, in Colombia, we worked with an Optimax, a multi-gas sensor sold by MRU Instruments Inc. The analyzed gases are the same as with the GA, four, by default, with others optional, and the computation of a balance is interpreted as N
2; in addition, the temperature is displayed. The application provides a time-indexed record, allowing users to visualize changes in gas concentrations through graphical representations that could be interesting for the sampling.
Figure 3 shows the available display for two measurements. Measuring may start before actioning the pump (
Figure 3b). As already stated, N
2 is not measured but deduced from the other gases, implying that its variations follow mainly O
2 peaks. With this sensor, the O
2 measurement is fine after 10 to 15 s and the H
2 measurement needs about 40 s. The shift between the O
2 and H
2 peaks is likely due to the sensors. With this real-time measurement, one may select the point and depth to take samples for laboratory analysis. These laboratory analyses have a cost and carrying out many field measurements without additional cost is key to selecting the right samples. One may note that without pumping, the signal reflects only air composition, which is coherent with the fact that one needs to extract the gas from the porosity of the soil. Since the soil permeability is not large, the pressure gap between the surface and 1 or a few meters is not large enough to create an immediate flow. This observation raises doubts about the ability of certain laboratories to carry out good analyses without a pump. Based on the pump power and the soil characteristics, the authors computed the volume of the soil that is drained around the tube, which is a cylinder about 8 cm in diameter [
20]. The height of the drained cylinder corresponds to the perforated part of the stem. This very small volume is reduced to nearly zero without pumping. Before the pump is started (
Figure 3b), the analyzer only measures the air in the stem.
In conclusion, it appears that the four main limitations of the current protocol are that (1) it separates drilling from setting up the measurement tube, (2) it limits the measurement depth, (3) it does not include monitoring and (4) the quantity and container of the gas samples are limited. The proposed solution addresses points 1, 2 and 4.
4. Results: H2 Content in the Soil of the Colombian Forelands
Four field campaigns have been organized between April 2022 and October 2023. More than 530 in situ soil gas measurements have been performed and the values are presented in
Figure 6. The data show the excellent potential of these two basins, with more than 290 values above 100 ppm and 170 above 200 ppm. There is, as usual, a large number of near-zero values. Even a few hundred ppm may look like a small quantity; H
2 emanations in the soil are not random nor ubiquitous. The average value of 158 ppm and the maximum value of 1629 ppm indicate that hydrogen generation is active in these basins. However, the goal of this article is the understanding of H
2 displacement in the soils, so we will only focus on the presentation of the results and discussion of those points.
In terms of the variability of the results, the fact that we have values at various depths does not change the general feature noted by previous authors: H2 content in the soil is highly variable, even in the areas, or especially in the areas, where leakage is ongoing. In other areas, the data are uniformly zero or near zero. This means that many measurements must always be made when evaluating an area.
Figure 6b shows this variability versus depth on a given zone. In the studied areas, there are no fairy circles but some vegetation anomalies, and the measures have been planned in cluster around selected locations. The distances between the measured points in a given site are, therefore, rather small, ranging from a few tens of a meter to less than 300 m. We can reasonably consider that the subsoil is fairly homogeneous at this scale, and the absence of any noticeable variation with depth suggests that the soil’s characteristics do not allow H
2 to accumulate on a specific layer. We can, therefore, interpret the vertical profiles as due to variations in H
2 flow.
4.1. Maximum H2 Content Versus Depth
The data set presented here includes values measured from near surface to 21 feet. As explained previously, when we open and close the system too many times to pump the gas, due to the characteristics of the soil in the studied areas, the grid that protects the tube often becomes plugged by a sort of mudcake and the gas cannot get through. Except for a few trials, the data have been recorded with a maximum of five extensions (i.e., 15 feet—4.6 m). The data for O
2 and H
2 contents in the soil are presented in
Figure 7. The first important information is that there is no general trend versus depth. In neither of the two basins shown here does the H
2 content of the gas clearly increase with depth. Conversely, the H
2 content is variable whatever the depth. It should also be noted that there are no major differences between the areas studied. For comparison, the O
2 content has been displayed; it is also variable and not depth-dependent. In these basins, the CH
4 is almost always close to zero and it is the CO
2 that balances the O
2.
4.2. Relation to Other Gases
As stated before, various gases are measured, and one may study the relationship between the H
2 and other gases.
Figure 7 shows the O
2 values versus depth and
Figure 8 shows the CH
4 and CO
2 versus H
2. Surprisingly, since Llanos and Putumayo Basins are petroleum provinces, the CH
4 content is always close to zero. The CO
2 content also remains fairly low in the Llanos Basin, at less than 6%, with an average of less than 0.5%. In the Putumayo Basin, the maximum is 13% and the average is 1.7%. Regarding the methane content, except for one point at 10% (located at a depth of 4 feet), all the values are close to zero, and the average is 0.09%. Interestingly,
Figure 8 shows that there is no relation between the CH
4 and H
2 content.
4.3. Gradient of the H2 Content Versus Depth
With the new tool developed for this study, the H
2 content in the soil has been measured in 180 points with an average of three measurements at different depths by points.
Figure 9 shows the vertical profiles of a couple of these points. The dots correspond to the data and the curve is an extrapolation to allow visualizing the H
2 gradient versus depth. Clearly, there is no constant and no typical profile. The gradients are sometimes negative and sometimes positive, and they change according to the depth. In the vast majority of cases, the H
2 content profiles display maximum and minimum values even along this rather short distance. Only in three cases is the gradient constantly positive; it means that the H
2 content is increasing between three and 12 feet. Due to the spacing of the measurements, one every 3 feet, the width of the peaks is, obviously, of little significance, but the presence of maximum and minimum values, close to zero, over similar vertical distances (P4, P5, L7) draws attention. Variation at wavelengths shorter than 3 feet cannot be quantified with this approach, but at one location, we performed a higher-frequency measurement to ensure that the values were representative (
Figure 10).
Leakage, and, therefore, H
2 emanations, are a dynamic system. For the large part, the H
2 is moving in the soil. To avoid the overinterpretation of these data, we have to keep in mind that the curves in
Figure 9 do not exactly represent vertical profiles at a given time. The time required to drill and measure four points is around ½ h. The duration of measurement was, for instance, 35 min for P1, 20 min for P2 and 36 min for P6, which was the deepest point, with five measurements. Interpreting this vertical profile as a picture at a given time is, therefore, an approximation, but it is the best we can do without having all the soil characteristics to recompute the signal. As an indication of the time required for the H
2 pulse to be transported upwards by 30 m, in the modeling carried out previously, it was around 6 h [
35]. It means that, for the 5 m displayed here, the gas will need 1 h in sandy soil; in shaly soil, the gas transport velocity is lower. The curves are, therefore, squeezed in comparison to the vertical profile at a given instant.
In the Putumayo Basin, the wavelength is between 5 and 9 feet, about 2.5 m ± 20%, whereas, in the Llanos, it is larger, about 12 feet. The difference could be due to different pulsing rates or different velocities in the upward transport. These values cannot be considered as quantitative, but they show the signal shape. The peak width is related to data spacing, but
Figure 10 gives a better definition of the signal; it corresponds to an acquisition almost every foot (the spacing has not been regular due to acquisition constrains). The H
2 peaks are much narrower, since, in
Figure 9, they were only defined by one value. The oxygen and CO
2 curves are more linear, with just three or four anomalous values, in comparison to the average value. A more global trend seems to start below 15 feet, with a decrease in O
2 content from 20 to 15%, but we do not have enough data to be definitive on that conclusion. The non-decrease in oxygen content for the first meters is, however, confirmed. It could be also noticed that the H
2 ambient noise, the minimum between the peaks, is around 100 ppm, which is a rather high value in comparison to that of many basins [
6,
26].
The H
2 content in the soil is known to be variable in time and space [
41]. All the published data confirm this observation, and various explanations have been proposed, which are all related to the way the H
2 is transported in the soils. This transport is influenced by the atmospheric pressure [
42], by the soil permeability and the microorganism consumption [
35,
36] and by the water content [
27]. Some scientists also tested the Earth and Moon tide effects [
43]; their conclusions were negative but other investigators think it is important. The 24 h wavelength highlighted by the Brazilian data set [
20] is correctly expdlained by the atmospheric pressure and the spatial variation by changes in soil characteristics (mineralogy, permeability, water content and microorganism consumption). For instance, lower permeability or an increase in water content will slow down the upward H
2 transport to the surface, resulting in a larger consumption by the microorganism, as quantified by Myagkiy and coauthors [
35,
36]. However, the charge of H
2, which rises from the rocks below the soil, can also be sporadic, if, for instance, if it is linked to a deeper fault valve effect [
17].
4.4. Model of H2 Rising Through a Soil
The two signals, soil breathing and deep flow, are distinct. In the context of native H
2 exploration, only the deep flow is of interest to us, whether sporadic or not. Spending too much time discussing what is happening at the near surface is irrelevant, but it must be understood in order to abstract from it. The data set of the long-term monitoring recorded in Brazil is compatible with the breathing of the soil, influenced by the air pressure with a signal on the order of 100 ppm and subject to the 24 h cycle and sporadic pulses, with much higher values and often less regular ones (see [
6,
20] for the details of the discussion and the data).
The origin of the deep H
2 flow may be diverse, directly from the rock that generates the H
2, but more likely, from a migration pathway or from a reservoir. The solubility of H
2 in water is high at large depths but decreases when the pressure, and so the depth, decreases [
44,
45]. In the case of H
2 dissolved in an aquifer, degassing will take place when the H
2 dissolved amount overpasses the solubility. The free gas thus released rises, like bubbles in an aqueous medium. The measured signal strongly suggests that this is what we are seeing in Colombia.
To our knowledge, data on H
2 content versus depth have only been published once, in the Carolina Bay in the US [
19]. The numerous data of H
2 versus depth recorded in Colombia are, therefore, a new and complementary element to the spatial variability already described many times over, and to the variability along the day and according to the dryness of the season, which are also well known [
27].
A linear increase in H
2 content with depth would have been compatible with a constant, or at least permanent, deep upward H
2 flow attenuated by near-surface consumption by microorganisms. This is not what has been observed. Gradient variations with depth strongly suggest a sporadic source. Such a deep signal results in variable gradient, both positive and negative, as modelling has already shown [
35].
In terms of pulse frequency, observing two maxima over 5 m, or at least one minimum and one maximum over 2 m, gives an indication of this frequency. It has to be higher than the 24 h or 6 h as in the simple model published by these authors [
35]. Reversing this signal would require in-depth knowledge of soil permeability, hydrometry (both in 3D) and H
2 transport mode (advection vs. diffusion). It could help to be more quantitative to pass from the data set of H
2 content in the soil to an H
2 flow rate, but acquiring so many data to characterize the soil will not help to decrease the uncertainties about the H
2 resources. Similarly, good quantification of oil or gas seeps does not help to prove reserves. We interpret these depth variations as proof that the signal is dynamic and, therefore, as proof that leakage is active today.
The hypotheses in
Figure 11 are that microorganism consumption decreases with depth, and that a clay layer close to the surface slows down the H
2 flow and causes it to be consumed. The zero near-surface signal in the center comes from these hypotheses. The signal decay from 100% at 30 m to 35% 8 h later at 7 m is due to consumption, but also to the spreading of the H
2 flow. This decrease in the maximum does not fit at all with the Colombian data set. Of course, there is no guarantee that the pulses are of the same intensity every time, but the fact that, on average, the maxima do not increase with depth (
Figure 8) also indicates that the consumption is not as assumed in this previous study [
35,
36].
When bubbling is observed in lakes or at sea, we see gas bubbles, and it is easy for anyone to realize that the gas content is not constant and cannot be. This has been pictured in New Caledonia [
46] and, more recently, in Indonesia [
47]. In soil, diffusion and variations in the permeability and/or absorption capacity of certain elements will make gas transport more complicated, but the source is probably just as sporadic.
Figure 12 is a simplified diagram, unlike
Figure 11, which is the result of a calculation, but it illustrates what is observed with these vertical profiles with rising pulse of H
2.
5. Conclusions
One of the main advantages of this new tool is the continuous measurement capability. There is no opening of the system and, therefore, no gas escape during perforation. The stem is pressed in without rotation and the measurements could be performed at various depths without reopening the hole. This method, therefore, enables accurate measurement of H2 concentrations at depths of up to 15 feet. It also allows us to sample the gas at various depths and the methods developed and presented in this paper offer enhanced reliability.
In return, the equipment is heavier, its use requires more time for a single measurement performed with the classical method (as presented in
Figure 2) and, since it is purely “hammer like”, in the case of outcropping rocks, the penetration will be limited. One good policy could be to use both methods in parallel, dividing up the team in the field. Quick and numerous measurements to make a map and deeper measurements could be performed on certain points to see the profile in depth.
Our dataset (180 different locations, 535 measurements) does not show, on average, higher values at 4 m than at 1 m. The H2 content values at 1 m look representative and it seems reasonable to continue to perform H2 concentration mapping with the classical method. However, the data clearly show the pulses of the H2 leakages. This confirms, once again, that you need to take different measurements at different times to know whether there is a flow of H2. The wavelength of the signal versus depth may be interpreted in terms of flow but only with a numerical tool.
The expected increase in H
2 with depth is not observed between 1 and 5 m. This may mean that microorganism consumption is only very active in the near surface and/or that, in that case, the upward H
2 flow is fast enough to not be affected. As a range of values, in the modeling performed previously, the upward H
2 near-surface velocity was in the range of 1 m per hours (strongly dependent on the water content and diffusion/advection ratio). The data do not show any correlation between the CO
2, CH
4 and H
2 content (
Figure 8 and
Figure 10), which is also not very compatible with a signal strongly influenced by microorganisms.
It was already known that the H2 signal within the soil is variable in 3D, (time and two spatial directions). We show, with these new data, that it is variable in 4D. This stands as another good reason to be careful not to interpret this type of data too quickly and, above all, not to consider scattered and too-sparse data as quantitative.
In addition, it could be noted that the H2 potential of these two basins, where H2 emanations are numerous and high, is very good.