Publications last year
2025
Efficiency and heat transport processes of low-temperature aquifer thermal energy storage systems : new insights from global sensitivity analyses
Authors: Luka Tas, Niels Hartog, Martin Bloemendal, David Simpson, Tanguy Robert, Robin Thibaut, Le Zhang and Thomas Hermans
Type: A1 Journal Article in Geothermal Energy
Abstract:
Aquifer thermal energy storage (ATES) has great potential to mitigate CO2 emissions associated with the heating and cooling of buildings and offers wide applicability. Thick productive aquifer layers have been targeted first, as these are the most promising hydrogeological context for ATES. Regardless, there is currently an increasing trend to target more complex aquifers such as low-transmissivity and alluvial aquifers or fractured rock formations. There, the uncertainty of subsurface characteristics and, with that, the risk of poorly performing systems is considerably higher. Commonly applied strategies to decide upon the ATES feasibility and well design standards for optimization need to be adapted. To further promote the use of ATES in such less favorable aquifers an efficient and systematic methodology evaluating the optimal conditions, while not neglecting uncertainty, is crucial. In this context, the distance-based global sensitivity analysis (DGSA) method is proposed. The analysis focuses on one promising thick productive aquifer, first used to validate the methodology, as well as a complex shallow alluvial aquifer. Through this method, multiple random model realizations are generated by sampling each parameter from a predetermined range of uncertainty. The DGSA methodology validates that the hydraulic conductivity, the natural hydraulic gradient and the annual storage volume dominate the functioning of an ATES system in both hydrogeological settings. The method also advances the state of the art in both settings. It efficiently identifies most informative field data ahead of carrying out the field work itself. In the studied settings, Darcy flux measurements can provide a first estimate of the relative ATES efficiency. It further offers a substantiated basis to streamline models in the future. Insensitive parameters can be fixed to average values without compromising on prediction accuracy. It also demonstrates the insignificance of seasonal soil temperature fluctuations on storage in unconfined shallow aquifers and it clarifies the thermal energy exchange dynamics directly above the storage volume. Finally, it creates the opportunity to explore different storage conditions in a particular setting, allowing to propose cutoff criteria for the investment in ATES. The nuanced understanding gained with this study offers practical guidance for enhanced efficiency of feasibility studies. It proves that the DGSA methodology can significantly speed up the development of ATES in more complex hydrogeological settings
Management strategies for improving maize yield and water productivity under water deficit and soil acidity
Authors: Desale Kidane Asmamaw, Kristine Walraevens, Habtamu Assaye, Fenta Nigate, Enyew Adgo and Wim Cornelis
Type: A2 Journal Article in Experimental Agriculture
Abstract:
Maize holds a key role in ensuring food security in Ethiopia, yet its productivity faces challenges due to water scarcity and soil acidity. Minimizing these problems is crucial to enhance maize yield and maintain food security. This research explored the effects of deficit irrigation (DI) combined with lime, manure, and inorganic fertilizer application on maize yield and water productivity (WP) in Koga, Ethiopia. Three levels of DI, namely 80%, 60%, and 50% of crop evapotranspiration (ETc), alongside 100% ETc (full irrigation) as a reference, were implemented for two consecutive seasons. Five integrated soil fertility management (ISFM) treatments were evaluated over two successive seasons: (i) combining 1.43 Mg ha-1 of lime with 3 Mg ha-1 of manure and full doses of urea + (NPSB, containing 18.9% Nitrogen, 37.7% Phosphorus, 6.95% Sulphur, and 0.1% Boron), referred to as inorganic fertilizer (L1); (ii) combining 1.15 Mg ha-1 of lime with 3 Mg ha-1 of manure and full doses of inorganic fertilizer (L2); (iii) combining 0.86 Mg ha-1 of lime with 3 Mg ha-1 of manure and full doses of inorganic fertilizer (L3); (iv) applying 3 Mg ha-1 of manure and full doses of inorganic fertilizer (M); and (v) using only full doses of inorganic fertilizer (C). Grain yield and biomass were measured at harvest from a 9 m2 sample area in each plot, with three replicates. The combined effects of DI, liming and manuring significantly influenced average grain yield and biomass. Across all irrigation scenarios, higher grain yield and biomass production were found with treatments L1, L2, L3, and M compared to treatment C. The highest WP was found with 50% ETc under all ISFM treatments. The lowest maize yield and WP were recorded with treatment C across all irrigation levels. Manuring combined with reduced irrigation increased grain yield, biomass, and WP compared to the use of inorganic fertilizer alone at 100% ETc. The combined use of lime and manure could mitigate the negative impact of DI on yield.
Microseismic/acoustic emission sources localization method using improved search algorithm for multiple media and irregular complex structures
Authors: Yihan Zhang Zhang, Longbin Yang, Longjun Dong, Jiachuang Wang, Qingchun Hu, Zhongwei Pei and Thomas Hermans
Type: A1 Journal Article in Measurement
Abstract:
Materials and engineering geological environments often present complex media and unique structures, complicating source localization due to diverse and nonlinear velocity models. This paper proposes a microseismic/acoustic emission (MS/AE) source localization method, termed ISACE-GPU, which utilizes an improved A* algorithm and GPU acceleration. The method first enables effective ray path searching for multi-media and irregular structures through the improved A* algorithm. Next, an optimized GPU parallel acceleration framework boosts localization efficiency. Finally, equidistant grid modeling meets diverse environmental needs. The results indicate that this method achieves relative localization errors of 2.5 % and 3.58 % of the monitoring network accuracy in simulations and mining applications, respectively, and improves localization accuracy by up to 92.5% and 28.11 % compared to the microseismic source location method without the pre-measuredwave velocity (MSLM-WV). Under appropriate grid node calculations, it reduces memory usage and increases computational speed by 27 times. The proposed method is suitable for seismic source localization in complex, large-scale environments, ensuring high accuracy and efficiency.
Thermo-hydro-mechanical modeling of geothermal energy systems in deep mines : uncertainty quantification and design optimization
Authors: Le Zhang, Anne-Catherine Dieudonné, Alexandros Daniilidis, Longjun Dong, Wenzhuo Cao, Robin Thibaut, Luka Tas and Thomas Hermans
Type: A1 Journal Article in Applied Energy
Abstract:
We explore the potential of using geothermal energy systems in deep mines to address two critical needs: sustainable energy production and effective mine cooling. This approach not only provides a green energy source but also helps manage high temperatures in deep mine environments. By injecting cold water into the subsurface, we enable heat exchange to cool the mines while producing geothermal energy. Our analysis focuses on how this process influences both the temperature and stability of the mine. Using advanced modeling, we examine the interactions between temperature, pressure, and mechanical stresses caused by geothermal operations. Key parameters, such as the distance between wells and mine workings, the rate of water injection, and material properties of the rock, were identified as crucial to achieving system efficiency and stability. This research offers a practical framework for optimizing geothermal energy systems in deep mines, balancing energy production with the safe and stable operation of mining activities.
2024
Assessing shallow groundwater quality around the Sheba leather tannery area, Wikro, North Ethiopia : a geophysical and hydrochemical study
Authors: Kaleab Adhena Abera, Berhane Abrha Asfaw, Yonatan Garkebo Doyoro, Tesfamichael Gebreyohanes, Abdelwassie Hussien, Gebremedhin Berhane, Miruts Hagos, Abadi Romha and Kristine Walraevens
Type: A2 Journal Article in Geosciences
Abstract:
This study aimed to investigate the shallow groundwater status around the Sheba Leather Tannery area, Wikro, North Ethiopia, through geophysical and hydrochemical methods. Seventeen Vertical Electrical Soundings (VESs) acquisitions, 4 upstream and 13 downstream, of the leather tannery area were conducted. Using the data, four geoelectric profiles were generated. The aquifers’ geoelectrical layers, depth, and lateral extent were delineated. The VES curves depicted three to four resistivity layers. These alternating layers of low, moderate, and high resistivity values, traced at different VES points, were attributed to the formations’ composition and the groundwater quality status. Besides the geophysical survey, 32 water samples were collected from the area. Parameters such as electrical conductivity (EC), total dissolved solids (TDSs), pH, major ions, and heavy metals were analyzed. Moreover, PHREEQC was used to determine the groundwater mineral saturation indices where most minerals, except halite, were found supersaturated. The quality status for drinking purposes was also evaluated using the water quality index (WQI), and the water was classified as good (56.3%), poor (37.5%), and very poor (6.2%). The sodium adsorption ratio (SAR) and the percentage of sodium (Na+%) were calculated, and the results indicated that the water is suitable for direct use in irrigation.
Estimation of groundwater recharge from groundwater level fluctuations and baseflow rates around Mount Meru, Tanzania
Authors: George Bennett, Ceven Shemsanga, Matthieu Kervyn de Meerendre and Kristine Walraevens
Type: A2 Journal Article in Groundwater for Sustainable Development
Abstract:
Estimating groundwater recharge, direct runoff and baseflow is essential for understanding groundwater resource availability and managing groundwater systems. This study estimates groundwater recharge, direct runoff and baseflow on two slopes of Mount Meru: the northern and southern slopes using the water -table fluctuation (WTF) method and baseflow separation technique. High -frequency groundwater level measurements in five shallow wells over three hydrological years from 2018 to 2021 were analysed, while streamflow data in four gauging stations over nine hydrological years from 2010 to 2019 were used. The results of the WTF method show that the aquifer undergoes an average recharge of 544 mm/year and 90 mm/year on the southwestern and north-eastern slopes, respectively. On average, this recharge is about 53% and 13% of the annual rainfall on each slope. The baseflow results show that the aquifer on the south-eastern and north-western slopes recharges an average of 88 mm/year and 54 mm/year, respectively, which is on average about 12% and 7% of annual rainfall, respectively. In general, the high recharge on the south-western slope is attributed to the high rainfall, and the high hydraulic conductivity and high hydraulic diffusivity of the pyroclastic deposits compared to the debris avalanche deposits on the north-eastern slope. In addition, debris avalanche deposits show homogeneous recharge conditions, while pyroclastic deposits show heterogeneous recharge conditions. The WTF method can be useful to identify areas of preferential recharge so that preferential groundwater flow paths can be mapped for focused recharge of surface runoff during the rainy season.
Estimation of groundwater recharge in a volcanic aquifer system using soil moisture balance and baseflow separation methods : the case of Gilgel Gibe catchment, Ethiopia
Authors: Fayera Gudu Tufa, Fekadu Fufa Feyissa, Adisu Befekadu Kebede, Beekan Gurmessa Gudeta, Wagari Mosisa Kitessa, Seifu Kebede Debela, Bekan Chelkeba Tumsa, Alemu Yenehun Beyene, Marc Van Camp and Kristine Walraevens
Type: A2 Journal Article in Hydrology
Abstract:
Understanding the recharge-discharge system of a catchment is key to the efficient use and effective management of groundwater resources. The present study focused on the estimation of groundwater recharge using Soil Moisture Balance (SMB) and Baseflow Separation (BFS) methods in the Gilgel Gibe catchment where water demand for irrigation, domestic, and industrial purposes is dramatically increasing. The demand for groundwater and the existing ambitious plans to respond to this demand will put a strain on the groundwater resource in the catchment unless prompt intervention is undertaken to ensure its sustainability. Ground-based hydrometeorological 36-years data (1985 to 2020) from 17 stations and satellite products from CHIRPS and NASA/POWER were used for the SMB method. Six BFS methods were applied through the Web-based Hydrograph Analysis Tool (WHAT), SepHydro, BFLOW, and Automated Computer Programming (PART) to sub-catchments and the main catchment to estimate the groundwater recharge. The streamflow data (discharge) obtained from the Ministry of Water and Energy were the main input data for the BFS methods. The average annual recharge of groundwater was estimated to be 313 mm using SMB for the years 1985 to 2020 and 314 mm using BFS for the years 1986 to 2003. The results from the SMB method revealed geographical heterogeneity in annual groundwater recharge, varying from 209 to 442 mm. Significant spatial variation is also observed in the estimated annual groundwater recharge using the BFS methods, which varies from 181 to 411 mm for sub-catchments. Hydrogeological conditions of the catchment were observed, and the yielding capacity of existing wells was assessed to evaluate the validity of the results. The recharge values estimated using SMB and BFS methods are comparable and hydrologically reasonable. The findings remarkably provide insightful information for decision-makers to develop effective groundwater management strategies and to prioritize the sub-catchments for immediate intervention to ensure the sustainability of groundwater.
Estimation of spatially distributed groundwater recharge in data-scarce regions
Authors: Ashebir Sewale Belay, Alemu Yenehun Beyene, Fenta Nigate, Seifu A. Tilahun, Mekete Dessie, Michael M. Moges, Margaret Chen, Derbew Fentie, Enyew Adgo, Jan Nyssen, et al.
Type: A1 Journal Article in Journal of Hydrology-regional studies
Abstract:
Study region: Upper Beles Basin, Ethiopia Focus: Data limitations significantly challenge accurate groundwater recharge estimation, especially in regions with sparse gauging stations and highly variable hydro-meteorological conditions, such as the Upper Beles Basin, Ethiopia. Using interpolated data often yields unreliable results. This study evaluates the use of remote sensing-based hydro-meteorological data to estimate spatially distributed groundwater recharge in such data-scarce regions, comparing it with point estimates derived from primary field data, using a case study of Upper Beles Basin. New hydrological insights for the region: The study employed TerraClimate and CHIRPS datasets to estimate spatially distributed recharge using the WetSpass model. Groundwater recharge was also estimated using the Water Table Fluctuation (WTF) and Chloride Mass Balance (CMB) methods, with data from 21 monitoring wells and 45 water samples. The estimated average annual recharge was 420 mm (WTF), 308 mm (CMB), and 365 mm (WetSpass). The substantial variability in point estimates reflects the basin’s recharge heterogeneity, indicating the risks of relying solely on point data. A strong correlation (72 %) was found between the WTF-derived point estimates and WetSpass-generated values. Recharge variability is influenced by land use in the lowlands and by slope, soil, and rainfall in the highlands. This research demonstrates that remote sensing-based data can yield more reliable recharge estimates in regions with sparse gauging stations and highly variable hydro-meteorological conditions.
Hydrogeochemical characterization and processes controlling groundwater chemistry of complex volcanic rock of Jimma area, Ethiopia
Authors: Wagari Mosisa Kitessa, Adisu Befekadu Kebede, Fayera Gudu Tufa, Beekan Gurmessa Gudeta, Alemu Yenehun Beyene, Bekan Chelkeba Tumsa, Seifu Kebede Debela, Fekadu Fufa Feyessa and Kristine Walraevens
Type: A1 Journal Article in Water
Abstract:
The sustainable management of groundwater in the Jimma area is complicated by a lack of comprehensive studies on its chemical makeup and the geochemical processes influencing its hydrochemistry. This research aims to fill that gap by examining 51 groundwater samples from various sources, including deep groundwaters, shallow groundwaters, hand-dug well groundwaters, surface waters, and springs within the area primarily consisting of complex volcanic rocks. The goal is to describe the hydrogeochemical characteristics and determine the key processes affecting groundwater composition in this volcanic area. The study identifies clear patterns in cation and anion concentrations. For deep groundwaters, the average cation concentration is ranked as Na+ > Ca2+ > Mg2+ > K+, while shallow groundwaters, hand-dug well groundwaters, surface waters, and springs show a ranking of Ca2+ > Na+ > Mg2+ > K+. The major anions are typically ordered as HCO3- > NO3- > Cl- > SO42-. The quantitative hydrogeochemical analysis indicates that the freshwater types in the region are primarily Ca-HCO3 and Ca-Mg-HCO3, with some highly mineralized Na-HCO3 waters also detected. The weathering of silicate minerals mainly drives the geochemical processes affecting groundwater chemistry. An increase in mineralization, suggested by saturation indices, points to a longer residence time underground, with deep groundwaters exhibiting the highest saturation levels and springs the lowest. This mineralization is especially significant for Mg-silicates and carbonates. Stability diagrams for feldspar minerals further demonstrate groundwater evolution along flow paths, revealing that shallow systems are in equilibrium with minerals like gibbsite, whereas deeper systems achieve stability with albite, Ca-montmorillonite, and microcline. Higher CO2 levels (10(-1.5) to 10(0.5) atm), likely from mantle-magma degassing, add more HCO3- to the deeper aquifers. This study offers the first thorough characterization of the groundwater composition in the Jimma area and provides important insights into the Jimma area’s hydrogeochemical development, establishing a basis for enhanced groundwater management within this intricate volcanic aquifer system.
Quantifying water loss in leaky micro-dam reservoir through water balance analysis and high-resolution water level data modeling
Authors: Gebremedhin Berhane and Kristine Walraevens
Type: A2 Journal Article in Water practice and Technology
Abstract:
This study aimed to assess the runoff, recharge, and response of a shallow aquifer to leakage from the Arato micro-dam reservoir (MDR). The assessment was conducted using the Soil Conservation Service Curve Number (SCS-CN), soil moisture balance (SMB), and diver (automatic data logger) measurements in both the MDR and a shallow hand-dug well. Recharge was estimated using the chloride mass balance (CMB) and water table fluctuation (WTF) methods. The results revealed that the annual runoff from the catchment was 48.8 mm, which accounted for approximately 0.71 million m(3). The yearly groundwater recharge was estimated to be 104, 92.8, and 100 mm using the SMB, CMB, and WTF methods, respectively. Furthermore, the water balance model of the Arato MDR indicated a leakage rate of 13.2 mm/day. It is noteworthy that the estimated leakage exceeded the seepage initially anticipated during the project’s design phase (9,965 m(3)/year). This research project highlights the significance of utilizing local climatic and physical data from the specific watershed under investigation when planning reservoirs and other water resources. It also underscores the importance of conducting thorough site investigations to accurately quantify hydraulic conductivity for leakage estimation purposes.
A multidimensional AI-trained correction to the 1D approximate model for Airborne TDEM sensing
Authors: Wouter Deleersnyder, David Dudal and Thomas Hermans
Type: A1 Journal Article in Computers & Geosciences
Abstract:
The computational resources required to solve the full 3D inversion of time -domain electromagnetic data are immense. To overcome the time-consuming 3D simulations, we construct a surrogate model, more precisely, a data -driven statistical model to replace the 3D simulations. It is trained on 3D data and predicts the approximate output much faster. We construct a surrogate model that predicts the discrepancy between a 1D subsurface model and a deviation of the 1D assumption. The latter response is fastly computable with a semi -analytical 1D forward model. We exemplify the approach on a two -layered case. The results are encouraging even with few training samples. Given the computational cost related to the 3D simulations, there are limitations in the number of training samples that can be generated. In addition, certain applications require a wide range of parameters to be sampled, such as the electrical conductivity parameters in a saltwater intrusion case. The challenge of this work is achieving the best possible accuracy with only a few thousand samples. We propose to view the performance in terms of learning gain, representing the gain from the surrogate model whilst still acknowledging a residual discrepancy. Our works open new avenues for effectively simulating 3D TDEM data.
Quantifying climate change variability for the better management of water resources: the case of Kobo Valley, Danakil Basin, Ethiopia
Authors: Mengesha Tesfaw, Mekete Dessie, Kristine Walraevens , Thomas Hermans, Fenta Nigate Abera, Tewodros Assefa and Kasye Shitu
Type: A2 Journal Article in Climate
Abstract:
Alterations in the hydrological cycle due to climate change are one of the key threats to the future accessibility of natural resources. This study used 12 GCM climate models from CMIP6 to evaluate future climate change scenarios by applying model performance measures and trend analysis in Kobo Valley, Ethiopia. The models were ranked based on their ability to analyze the historical datasets. The result of this study showed that the outputs of the FIO-ESM-2-0 CIMP6 model had a good overall ranking for both precipitation and temperature. After bias correction of the model-based projections with the observed data, the average annual precipitation in the average scenario (SSP2-4.5) decreased by 4.4% and 13% in 2054 and 2084, respectively. Similarly, in the worst-case scenario (SSP5-8.5), by the end of 2054 and 2084, decreases of 4% and 12.8%, respectively, were predicted. The average annual maximum temperature under the SSP2-4.5 scenario increased by 1.5 °C in 2054 and by 2.1 °C in 2084. The average annual maximum temperature under the worst-case (SSP5-8.5) scenario increased by 1.7 °C in 2054 and by 3.2 °C in 2084. In the middle scenario (SSP4.5), the average annual minimum temperature increased by 2.2 °C in 2054 and by 3 °C in 2084. The average annual minimum temperature under the worst-case (SSP5-8.5) scenario increased by 2.6 °C in 2054 and by 4.3 °C in 2084. The seasonal variability in precipitation in the studied valley will decrease in the winter and increase in the summer. A decrease in precipitation combined with an increase in temperature will strengthen the risk of drought events in the future.
Assessing and Improving the Robustness of Bayesian Evidential Learning in One dimension for inverting Time-Domain main Electromagnetic Data: Introducing a new Threshold Procedure.
Authors: Arsalan Ahmed, Lukas Aigner, Hadrien Michel, Wouter Deleersnyder, David Dudal, Adrian Flores Orozco, Thomas Hermans
Type: A1 Journal Article in Water Journal
Abstract:
Imagine you’re trying to figure out what’s hidden underneath a giant blanket without lifting it. Instead of guessing based on a few peeks under the blanket (like drilling holes in the ground), scientists have ways to ‘feel’ what’s beneath from the surface using geophysical methods like “Time domain electromagnetic data”, which are much cheaper and less harmful to the environment. But the catch is, these ‘feels can be interpreted in many ways, making it tricky to know exactly what’s under there.
This study talks about a smart method, kind of like using an insightful guide, to make better guesses about what’s hidden under our ‘blanket’—the Earth’s surface. This guide, called Bayesian evidential learning, helps us to understand the uncertainties in our guesses, especially when we’re looking at how things might change over time and space underground, like tracking unwanted saltwater getting into freshwater areas in Vietnam.
The cool part? The study shows that with some clever tricks (“the introduction of a new threshold”), we can avoid super expensive or complex computations and still get trustworthy results. This means we can take better care of our environment and resources without breaking the bank or causing unnecessary damage.
Quantifying salinity in heterogeneous coastal aquifers through ERT and IP: Insights from laboratory and field investigations
Authors: Diep Cong-Thi, Linh Pham Dieu, David Caterina, Xavier De Pauw, Huyen Dang Thi, Huu Hieu Ho, Frédéric Nguyen, Thomas Hermans
Type: A1 Journal Article in Journal of Contaminant Hydrology
Abstract:
ERT and IP are considered to be sensitive for classifying saltwater-bearing and clay-rich sediments in the coastal aquifers. To interpret the ERT data, we rely on the petrophysical relationship with mention of the surface conductivity contribution to the total bulk conductivity.
Firstly, the sedimentary samples collected from the co-located boreholes were classified according to their particle size distribution and analyzed in the lab using a spectral IP in controlled salinity conditions to derive their formation factors, surface conductivity, and normalized chargeability. Secondly, the deduced thresholds are applied on the field to distinguish clay-bearing sediments from brackish sandy sediments.
Application of this approach in the Luy River catchment revealed variations in formation factors in various unconsolidated sedimentary classifications, ranging from 4.0 to 8.9 for coarse-grained sand and clay-bearing mixtures, respectively. The presence of clay content is detected in comparison with the normalized chargeability conditions over 1.0 mS.m-1. Clay-bearing sediments were found to be predominantly distributed in discontinuous small lenses, intercalating with sandy layers. This highlights the heterogeneous nature of these coastal aquifers (Diep Cong-Thi et al, 2024).
We hope that these methods can be widely applied to the heterogeneous clay-bearing aquifers in the various coastal regions.
This research is funded by VLIR-UOS and the Belgian development cooperation through the grant VN2019TEA494A103, and theSpecial Research Fund (BOF), Ghent University.
Response to soil compaction of the electrical resistivity tomography, induced polarisation, and electromagnetic induction methods: a case study in Belgium
Authors: Danial Mansourian, Adriaan Vanderhasselt, Wim Cornelis and Thomas Hermans
Type: A1 Journal Article in Soil Research
Abstract:
Context: soil Compaction acts at different scales and is challenging to measure on field scales.
Aims: to evaluate soil compaction under a controlled traffic experiment, using three different geophysical methods.
Methods: Electrical Resistivity Tomography (ERT), Electromagnetic Induction (EMI), and Induced Polarization (IP) were selected to map soil compaction. Two different ERT arrays and EMI geometries were selected with different spacings. The influences of configuration, electrode spacing, and the depth of investigation Index (DOI) were evaluated. Soil physical properties were measured in the Laboratory and in the field. Error models were developed to assess the accuracy of the ERT profiles and later correlated with EMI and soil physical results.
Key results: Penetration resistance measurements identified a compacted layer at 25 to 35 cm depth with a maximum value of 5 MPa under fixed tracks and bulk density of 1.52 Mg m-3, while lowest values were 1.4 MPa and 1.36 Mg m-3. The Dipole-Dipole 10 cm array was more accurate towards both soil properties and locating the zones of high resistivity. The IP method identified chargeability anomalies at the same depth as the resistivity anomalies, possibly indicating a similar origin. The EMI test was less successful in accurately determining the locations of the conductive areas.
Conclusions: A clear relationship between the absolute value of the resistivity/conductivity signals with the level of compaction was not found, yet patterns of lateral variations in resistivity were identified.
Implications: Further studies are needed to establish the concrete relationship between soil compaction and geophysical signals.
Integrated methodology to link geochemical and geophysical-lab data in a geophysical investigation of a slag heap for resource quantification
Authors: Itzel Isunza Manrique, Thomas Hermans , David Caterina, Damien Jougnot, Benoît Mignon, Antoine Masse and Frédéric Nguyen
Type: A1 Journal Article in Journal of Environmental Management
Abstract:
The increasing need to find alternative stocks of critical raw materials drives to revisit the residues generated during the former production of mineral and metallic raw materials. Geophysical methods contribute to the sustainable characterization of metallurgical residues inferring on their composition, zonation and volume(s) estimation. Nevertheless, more quantitative approaches are needed to link geochemical or mineralogical analyses with the geophysical data. In this contribution, we describe a methodology that integrates geochemical and geophysical laboratory measurements to interpret geophysical field data solving a classification problem. The final aim is to estimate volume(s) of different types of materials to assess the potential resource recovery. We illustrate this methodology with a slag heap composed of residues from a former iron and steel factory. First, we carried out a 3D field acquisition using electrical resistivity tomography (ERT) and induced polarization (IP), based on which, a sampling survey was designed. We conducted laboratory measurements of ERT, IP, spectral induced polarization (SIP), and X-ray fluorescence analysis, based on which, 4 groups of different chemical composition were identified. Then we carried out a 3D probabilistic classification of the field data, based on 2D kernel density estimators (for each group) fitted to the inverted data collocated with the samples. The estimated volumes based on the classification model were: 4.17 × 103 m3 ± 12 %, 1.888 × 105 m3 ± 12 %, 59.4 × 103 m3 ± 19 %, and 2.30 × 104 m3 ± 21% for the groups ordered with an increasing metallic content. The uncertainty ranges were derived from comparing the volumes with and without considering the probabilities associated to the classification. We found that a representative sampling and the definition of the KDE bandwidths are defining elements in the classification and ultimately the estimation of volumes. This methodology is suitable to quantitatively interpret geophysical data in terms of the geochemical composition of the materials, integrating uncertainties both in the classification and the estimation of volumes. Furthermore, several crucial elements in the investigation of metallurgical residues could be applied in a real case study, e.g., geophysical field acquisition, sampling and lab measurements.
2023
A new Bayesian framework for the interpretation of geophysical data
Authors: Hadrien Michel, Frédéric Nguyen, Thomas Hermans
Type: Dissertation
Abstract: Providing images of the subsurface from ground-based datasets is at the heart of the geophysicist’s work. Multiple approaches have been applied to tackle this task. Most of the time, this task is performed in a deterministic framework, meaning that for a given dataset, a single model is provided to explain the data. However, those deterministic approaches lack the ability to provide reasonable uncertainty estimations, that take into account the non-unicity of the solution, noise in the data and modelling error. To provide precise and accurate models of the subsurface along with uncertainty, geophysicists use probabilistic approaches. Those approaches are able to sample the ensemble of a priori possible models (the prior) in order to extract models that can reasonably explain the datasets (the posterior). Such approaches, even though superior in terms of the reliability of their results, are rarely applied in practice due to their significant computational requirements. In this manuscript, the aim is to propose a new Bayesian framework to interpret those geophysical datasets. This new framework, called Bayesian Evidential Learning, promises to enable a fast, precise and accurate estimation of the uncertainty. This framework is applied and adapted for 1D geophysical datasets (BEL1D). The new and adapted framework presents several advantages when compared to classical probabilistic approaches: from fast computations due to the limited number of forward runs needed, to providing insight about the experiment sensitivity and the validity of the prior. Moreover, it benefits from its construction as a Machine Learning algorithm, leading to quasi-instantaneous models of uncertainty.
Advancing measurements and representations of subsurface heterogeneity and dynamic processes: towards 4D hydrogeology
Authors: Thomas Hermans (UGent) , Pascal Goderniaux, Damien Jougnot, Jan H. Fleckenstein, Philip Brunner, Frédéric Nguyen, Niklas Linde, Johan Alexander Huisman, Olivier Bour, Jorge Lopez Alvis, et al.
Type: A1 Journal Article in Hydrology and Earth System Sciences
Abstract: Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability generally leads to effective behaviors and emerging phenomena that cannot be predicted from conventional approaches based on homogeneous assumptions and models. However, it is not always clear when, why, how, and at what scale the 4D (3D + time) nature of the subsurface needs to be considered in hydrogeological monitoring, modeling, and applications. In this paper, we discuss the interest and potential for the monitoring and characterization of spatial and temporal variability, including 4D imaging, in a series of hydrogeological processes: (1) groundwater fluxes, (2) solute transport and reaction, (3) vadose zone dynamics, and (4) surface-subsurface water interactions. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. We then highlight recent innovations that have led to significant breakthroughs in high-resolution space-time imaging and modeling the characterization, monitoring, and modeling of these spatial and temporal fluctuations. We finally propose a classification of processes and applications at different scales according to their need and potential for high-resolution space-time imaging. We thus advocate a more systematic characterization of the dynamic and 3D nature of the subsurface for a series of critical processes and emerging applications. This calls for the validation of 4D imaging techniques at highly instrumented observatories and the harmonization of open databases to share hydrogeological data sets in their 4D components.
An integrated geoarchaeological approach to the investigation of multi-period prehistoric settlements – the case of Neolithic Drenovac
Authors: Philippe De Smedt (UGent) , Charles French, Timothy Kinnaird, Tonko Rajkovača, Aleksandar Milekić, Petros Chatzimpaloglou, Jeroen Verhegge (UGent) , Thomas Hermans (UGent) , Gaston Mendoza Veirana (UGent) , Djurdja Obradović, et al.
Type: U Conference
Abstract: A multi-method geoarchaeological investigation was performed to reconstruct multi-phase Neolithic settlement. Invasive and non-invasive surveys showed potential for providing archaeological and environmental landscape data in this complex setting. Large-area geophysical surveys showed potential for deriving stratigraphic information.
Assessing the potential of low‑transmissivity aquifers for aquifer thermal energy storage systems: a case study in Flanders (Belgium)
Authors: Luka Tas (UGent), David Simpson and Thomas Hermans (UGent)
Type: A1 Journal Article in Hydrogeology Journal
Abstract: The Member States of the European Union pledged to reduce greenhouse gas emissions by 80–95% by 2050. Shallow geothermal systems might substantially contribute by providing heating and cooling in a sustainable way through seasonally storing heat and cold in the shallow ground (<200 m). When the minimum yield associated with the installation of a cost-effective aquifer thermal energy storage (ATES) system cannot be met, borehole thermal energy storage, relying mostly on the thermal conductivity of the ground, is proposed. However, for large-scale applications, this requires the installation of hundreds of boreholes, which entails a large cost and high disturbance of the underground. In such cases, ATES systems can nevertheless become interesting. This paper presents a case study performed on a Ghent University campus (Belgium), where the feasibility of ATES in an area with a low transmissivity was determined. The maximum yield of the aquifer was estimated at 5 m3/h through pumping tests. Although this low yield was attributed to the fine grain size of the aquifer, membrane filtering index tests and long-term injection tests revealed that the clogging risk was limited. A groundwater model was used to optimize the well placement. It was shown that a well arrangement in a checkerboard pattern was most effective to optimize the hydraulic efficiency while maintaining the thermal recovery efficiency of the ATES system. Hence, for large-scale projects, efficient thermal energy storage can also be achieved using a (more cost-effective) ATES system even in low-permeability sediments.
Comparison of soft indicator and poisson kriging for the noise-filtering and downscaling of areal data : application to daily COVID-19 incidence rates
Authors: Pierre Goovaerts, Thomas Hermans (UGent) , Peter Goossens and Ellen Van De Vijver (UGent)
Type: A1 Journal Article in ISPRS International Journal of Geo-Information
Abstract: This paper addresses two common challenges in analyzing spatial epidemiological data, specifically disease incidence rates recorded over small areas: filtering noise caused by small local population sizes and deriving estimates at different spatial scales. Geostatistical techniques, including Poisson kriging (PK), have been used to address these issues by accounting for spatial correlation patterns and neighboring observations in smoothing and changing spatial support. However, PK has a limitation in that it can generate unrealistic rates that are either negative or greater than 100%. To overcome this limitation, an alternative method that relies on soft indicator kriging (IK) is presented. The performance of this method is compared to PK using daily COVID-19 incidence rates recorded in 2020–2021 for each of the 581 municipalities in Belgium. Both approaches are used to derive noise-filtered incidence rates for four different dates of the pandemic at the municipality level and at the nodes of a 1 km spacing grid covering the country. The IK approach has several attractive features: (1) the lack of negative kriging estimates, (2) the smaller smoothing effect, and (3) the better agreement with observed municipality-level rates after aggregation, in particular when the original rate was zero.
Editorial : geomechanics and induced seismicity for underground energy and resources exploitation
Authors: Longjun Dong, Wenzhuo Cao and Thomas Hermans (UGent)
Type: Editorial in Frontiers in Earth Science
Keywords: General Earth and Planetary Sciences, geomechanics, microseismic monitoring, underground engineering, induced seismicity geo-energy, mineral resources
Flexible quasi-2D inversion of time-domain AEM data, using a wavelet-based complexity measure
Authors: Wouter Deleersnyder (UGent) , Benjamin Maveau, Thomas Hermans (UGent) and David Dudal (UGent)
Type: A1 Journal Article in Geophysical Journal International
Abstract: Regularization methods improve the stability of ill-posed inverse problems by introducing some a priori characteristics for the solution such as smoothness or sharpness. In this contribution, we propose a multidimensional scale-dependent wavelet-based l(1)-regularization term to cure the ill-posedness of the airborne (time-domain) electromagnetic induction inverse problem. The regularization term is flexible, as it can recover blocky, smooth and tunable in-between inversion models, based on a suitable wavelet basis function. For each orientation, a different wavelet basis function can be used, introducing an additional relative regularization parameter. We propose a calibration method to determine (an educated initial guess for) this relative regularization parameter, which reduces the need to optimize for this parameter and, consequently, the overall computation time is under control. We apply our novel scheme to a time-domain airborne electromagnetic data set in Belgian saltwater intrusion context, but the scheme could equally apply to any other 2D or 3D geophysical inverse problem.
Footprint of fresh submarine groundwater discharge in the Belgian coastal zone : an overview study
Authors: Marieke Paepen (UGent) , Kristine Walraevens (UGent) and Thomas Hermans (UGent)
Type: A2 Journal Article in Leading Edge
Abstract: The presence of fresh groundwater is not limited to land; it also extends offshore and discharges as submarine groundwater discharge (SGD). The freshwater component of SGD (fresh submarine groundwater discharge [FSGD]) can be detected using geophysical techniques that are sensitive to salinity such as resistivity measurements. However, these measurements are often limited to either the land or marine realm, neglecting the land-marine interface. In this study, we focus on this gap by combining onshore and offshore techniques to assess variability in the FSGD footprint near the Belgian coastline through electrical resistivity tomography and continuous resistivity profiling. The difficult working conditions of the highly dynamic North Sea make this offshore survey one of the first of its kind. The footprint varies from limited outflow on the upper beach (e.g., Wenduine) to discharge around and below the low water line (e.g., De Panne, Oostduinkerke, and Knokke-Heist) in the studied areas. The occurrence, footprint, and quantity of SGD seem to be controlled by the presence and size of dune formations that constitute freshwater resources along the shore. Heterogeneity can also play a determining factor in FSGD location.
Ground validation of satellite-based precipitation estimates over poorly gauged catchment: the case of the Drâa basin in Central-East Morocco
Authors: Athmane Khettouch, Mohammed Hssaisoune, Thomas Hermans (UGent) , Aziz Aouijil and Lhoussaine Bouchaou
Type: U Journal Article in Mediterranean Geoscience Reviews
Abstract: In the most ungauged areas, lack of precipitation information limits the accuracy of water balance approaches. The ungauged Drâa river basin (DRB) in center eastern of Morocco is one of the ten driest basins worldwide with lack of adequate rainfall dataset for water resources management. This study assess five satellite precipitation datasets (P-datasets) with high spatial resolution (0.0375°–0.1°) and long period of record (> 40 years), namely CHIRPS V2.0, MSWEP V2.8, PERSIANN-CCS-CDR, TAMSAT V3.1, and ERA5-Land with reference to ground rain observations, based on continuous, categorical, and volumetric indices, and at various elevations, rainfall intensities, and temporal scales (i.e., monthly, seasonal, and sub-seasonal). Moreover, the ability to detect extreme precipitation event and the suitability of the conventional rain gauge to increase the magnitude of error was also quantified. ERA5-Land followed by MSWEP V2.8 have shown the best statistical scores at different times, intensity and elevation scales, while CHIRPS V2.0, PERSIANN-CCS-CDR and TAMSAT V3.1 provide poor estimations of rainfall with high sensitivity to the complexity of terrains. CHIRPS V2.0 was more efficient in detecting wet months, but a large event detected were not confirmed by ground observation. The magnitude of error tends to decrease during summer and for precipitation within the range of 1–12.5 mm for all P-datasets. However, all the five products underestimate the frequency of dry months and overestimate high precipitation intensity. Our findings not only recommend ERA5-Land and MSWEP V2.8 datasets as alternative to rain gauge but also describe the limitation of CHIRPS V2.0, TAMSAT V3.1 and PERSIANN-CSS-CDR for such hydrological studies.
Machine learning for Bayesian experimental design in the subsurface
Authors: Robin Thibaut (UGent)
Type: Dissertation
Abstract: Accurate modeling of the subsurface, a complex and heterogeneous environment that plays a crucial role in the Earth’s water cycle, is challenging due to sparse and incom- plete data. We can reduce the uncertainty associated with subsurface predictions, such as groundwater flow and contaminant transport, by conducting additional observations and measurements in the subsurface. However, practical and economic considerations frequently limit the number of measurements and their locations, such as land occupa- tion, which may limit the number of wells that can be drilled. In this dissertation, we propose simulation-driven methods to reduce uncertainty in subsurface predictions by identifying the most informative data sets to gather. Our method, which is based on Bayesian optimal experimental design and machine learning, determines the nature and location of these data sets, which can include measurements of groundwater levels, tem- perature, and other parameters collected through active or passive sensing methods such as pumping tests, tracer tests, and geophysical surveys. This dissertation is the first to use Bayesian Evidential Learning (BEL) for optimal experimental design, allowing for the optimization of data source locations and the comparison of the utility of different data sources. BEL is a framework for prediction that combines Monte Carlo sampling and machine learning in order to learn a direct relationship between predictor and target variables generated by a simulation model. We demonstrate the efficacy of our methods in three groundwater modeling case studies: (i) wellhead protection area delineation, (ii) an aquifer thermal energy storage monitoring system, and (iii) groundwater-surface water interaction. The case studies show that our approach can significantly reduce the uncertainty in subsurface predictions and guide further subsurface exploration. The first case study, in particular, uses the Traveling Salesman Problem to introduce a novel ap- proach to wellhead protection area delineation. The second case study, which compares well and geophysical data for temperature monitoring, introduces a new method for com- bining observations from multiple data sources in a latent space of the original data. The third case study introduces the Probabilistic Bayesian neural network (PBNN) method to BEL and transitions from a static experimental design framework to a sequential experimental design framework, which estimates groundwater-surface water interaction fluxes from temperature data. We have also developed a Python package, SKBEL, that implements our methods and can be used for a variety of Earth Science applications. Overall, this dissertation demonstrates the utility of BEL for optimal experimental de- sign in groundwater modeling, highlights the potential of BEL for predictive modeling in Earth Sciences, and opens up new avenues for data and simulation-driven subsurface modeling.
Quantitative interpretation of geoelectric inverted data with a robust probabilistic approach
Authors: Itzel Isunza Manrique, David Caterina, Frederic Nguyen and Thomas Hermans (UGent)
Type: A1 Journal Article in Geophysics
Abstract: The nonuniqueness of the solution to the geophysical inverse problem can lead to misinterpretation while characterizing the subsurface. To tackle this situation, ground-truth information from excavations and wells can be used to improve, calibrate, and interpret inverted models. We refer to quantitative interpre-tation as the decision analysis based on probability theory, which is focused on solving a classification problem. First, we present a probabilistic approach to classify the different types of materials or “categories” observed in borehole logs using multiple data sources: inverted 2D electrical resistivity tomography and in-duced polarization data and the positions (x, z) of these boreholes. Then, using the Bayes’ rule and permanence of ratios, we com-pute the joint conditional probabilities of each category, given all data sources in the whole inverted model domain. We validate this approach with synthetic data modeling for a complex anthropogenic-geologic scenario and using real data from an old landfill. Afterward, we assess the performance of the probabilistic approach for classification and compare it with the machine learn-ing algorithm of multilayer perceptron (MLP). In addition, we analyze the effect that the different data sources and the number of boreholes (and their distribution) have on both approaches with the synthetic case. Our results indicate that the MLP performance is better for delineating the different categories where the lateral contrasts in the synthetic resistivity model are small. Never-theless, the classification obtained with the probabilistic approach using real data seems to provide a more geologically realistic distribution. We conclude that the probabilistic approach is robust for classifying categories when high spatial heterogeneity is ex-pected and when ground-truth data are limited or not sparsely distributed. Finally, this approach can be easily extended to inte-grate multiple geophysical methods and does not require the optimization of hyperparameters as for MLP.
The quantitative meaning of resistivity data in a coastal setting : a Belgian case study
Authors: Marieke Paepen (UGent) , Wouter Deleersnyder (UGent) , Kristine Walraevens (UGent) and Thomas Hermans (UGent)
Type: C3 Conference
Abstract: In coastal areas, the natural groundwater flow is affected by human activities, such as managed aquifer recharge (MAR) and groundwater extraction. They can induce saltwater intrusion and impact the fresh submarine groundwater discharge (FSGD). Resistivity methods, such as electrical resistivity tomography (ERT) and continuous resistivity profiling (CRP) are easy to use and very effective to assess the distribution of salt and freshwater in coastal environments. The Western Belgian coast, De Panne and Koksijde, was already investigated with ERT and CRP by Paepen et al. (2022; 2020). In this area, the source of FSGD is a sandy dune ridge of around 2.5 km wide. Now, we compare the FSGD footprint in front of De Panne and Koksijde to other Belgian coastal sites (Raversijde, Wenduine, Knokke-Heist, and Zwin), which have a different structure of the phreatic aquifer and a much smaller dune belt.
The quantitative interpretation of ERT and CRP is not straightforward, but image appraisal tools – such as the model resolution matrix (R), cumulative sensitivity matrix (S), and depth of investigation index (DOI) – can aid (Caterina et al., 2013). To be able to quantitatively assess the resistivity inversion models, five synthetic models were created (Paepen et al., 2022). These models reflect the present situation of the Western Belgian coast, where we find freshwater outflow on the lower beach or below the low water line. Based on the inversion models of the synthetic cases, the model resolution matrix, cumulative sensitivity matrix, and DOI (Oldenburg & Li, 1999) were calculated. The image appraisal tools were then compared to the error on the salinity for each cell in the inversion model (we find an error below 0.05 acceptable). This allows to define a threshold of the different image appraisal tools for which the model can be quantitatively assessed and to apply them to the field data. The thresholds reveal that no quantitative interpretation is possible for the zones of FSGD and that the FSGD resistivity is underestimated by the inversion process, so the salinity of the outflow is overestimated. Nevertheless, lateral qualitative changes can be deduced from the inversion models.