Imaging hydrological dynamics in karst unsaturated zones by time-lapse electrical resistivity tomography

Jian Zhang; Colette Sirieix; Dominique Genty; Fabien Salmon; Cécile Verdet; Sylvain Mateo; Shan Xu; Stéphane Bujan; Ludovic Devaux; Marie Larcanché

doi:10.1016/j.scitotenv.2023.168037

Imaging hydrological dynamics in karst unsaturated zones by time-lapse electrical resistivity tomography

Sci Total Environ. 2024 Jan 10:907:168037. doi: 10.1016/j.scitotenv.2023.168037. Epub 2023 Oct 23.

Authors

Affiliations

¹ Université de Bordeaux, CNRS, Bordeaux INP, I2M, UMR 5295, F-33400 Talence, France; Arts et Metiers Institute of Technology, CNRS, Bordeaux INP, INRAE, I2M, UMR 5295, F-33400 Talence, France; Environnements et Paléoenvironnements Océaniques et Continentaux (EPOC), UMR CNRS, 5805, Université de Bordeaux, 33615 Pessac Cedex, France. Electronic address: jian.zhang@u-bordeaux.fr.
² Université de Bordeaux, CNRS, Bordeaux INP, I2M, UMR 5295, F-33400 Talence, France; Arts et Metiers Institute of Technology, CNRS, Bordeaux INP, INRAE, I2M, UMR 5295, F-33400 Talence, France. Electronic address: colette.sirieix@u-bordeaux.fr.
³ Environnements et Paléoenvironnements Océaniques et Continentaux (EPOC), UMR CNRS, 5805, Université de Bordeaux, 33615 Pessac Cedex, France.
⁴ Université de Bordeaux, CNRS, Bordeaux INP, I2M, UMR 5295, F-33400 Talence, France; Arts et Metiers Institute of Technology, CNRS, Bordeaux INP, INRAE, I2M, UMR 5295, F-33400 Talence, France.
⁵ School of Civil Engineering and Mechanics, Yanshan University, Qinhuangdao, PR China.

PMID: 37879471
DOI: 10.1016/j.scitotenv.2023.168037

Abstract

The hydrodynamics of karst terrain are highly complex due to the diverse fractures and reservoirs within limestone formations. The time delay between rainfall events and subsequent flow into reservoirs exhibits significant variability. However, these hydrological processes are not easily visualized in karst topography. Subsurface geophysics, specifically 2D time-lapse electrical resistivity tomography (ERT), provides an effective method for studying the relationships between hydrological and geophysical features. In our research, we adopted ERT in the Karst Critical Zone (KCZ) to visualize specific karstic zones, including cave galleries, water storage reservoirs, wetting fronts, soil layers, and potential preferential flow paths down to a depth of 20 m. To capture spatial and seasonal variations in resistivity, we presented a comprehensive approach by combining sixteen inversion models obtained between February 2020 and September 2022 above the Villars Cave in SW-France-a well-known prehistoric cave. We used a multi-dimensional statistical technique called Hierarchical Agglomerative Clustering (HAC) to create a composite model that divided the synthetic ERT image into eight clusters representing different karst critical zones. The ERT image clearly visualized the cave gallery with high resistivity values that remained consistent throughout the seasons. Our analysis revealed a close seasonal relationship between water excess and resistivity variations in most infiltration zones, with time delays increasing with depth. The karst reservoirs, located at significant depths compared to other clusters, displayed sensitivity to changes in water excess but were primarily affected by fluctuations in water conductivity, particularly during summer or dry periods. These findings have significant implications for predicting rainwater infiltration pathways into caves, thereby assisting in the conservation and preservation of prehistoric caves and their cultural heritage.

Keywords: Drip rate; Electrical resistivity tomography; Hierarchical agglomerative clustering; Infiltration; Karst reservoirs.