Seismic Refraction Tomography Research Articles

Prospection of faults on the Southern Erftscholle (Germany) with individually and jointly inverted refraction seismics and electrical resistivity tomography

As part of the Lower Rhein Embayment (LRE), the Southern Erft block is characterized by a complex tectonic setting that influences hydrological and geological conditions on a local as well as regional level. The study area is located in the south of North Rhine-Westphalia and traversed by several NW-SE-oriented fault structures. Since the tectonic structures were located by past studies based on a sparse foundation of geological data, the positions include considerable uncertainties. Therefore, it was decided to re-evaluate and refine the assumed fault locations by conducting geophysical measurements. Seismic Refraction Tomography (SRT) as well as Electrical Resistivity Tomography (ERT) was performed along seven measurement profiles with a length of up to 1.1 km. In addition to compiling individual resistivity and velocity models for all deduced measurements, ERT and SRT datasets were cooperatively inverted using the Structurally Coupled Cooperative Inversion (SCCI). This algorithm strengthens structural similarities between velocity and resistivity by adapting the individual regularizations after each model iteration. Previously assumed locations of the tectonic structures diverge from the new evidence based on ERT and SRT surveys. Especially in the western and eastern parts of the research area, differences between the survey results and formerly assumed locations are in the order of 100 m. Seismic and geoelectric measurements further indicate a fault structure in the southern part of the area, which remained undetected by past studies. The cooperative inversions do not improve the geophysical models qualitatively, since the individually inverted datasets already provide results of good quality and resolution.

Application of Experimental Configurations of Seismic and Electric Tomographic Techniques to the Investigation of Complex Geological Structures

The main purpose of this study is the subsurface investigation of two complex geological environments focusing on the improvement of data acquisition and processing parameters regarding electric and seismic tomographic techniques. Two different study areas, in central–east Peloponnese and SE Attica, were selected, where detailed geological mapping and surface geophysical survey were carried out. The applied geophysical survey included the application of electrical resistivity tomography (ERT), seismic refraction tomography (SRT) and ground penetrating radar (GPR). The geoelectrical measurements were acquired with different arrays and electrode configurations. Moreover, various types of seismic sources were used at seventeen shot locations along the seismic arrays. For the processing of geoelectrical data, clustered datasets were created, increasing the depth of investigation and discriminatory capability. The seismic data processing included the following: (a) the creation of synthetic models and seismic records to determine the effectiveness and capabilities of the technique, (b) spectral analysis of the seismic records to determine the optimal seismic source type and (c) inversion of the field data to create representative subsurface velocity models. The results of the two techniques successfully delineated the complex subsurface structure that characterizes these two geological environments. The application of the ERT combined with the SRT are the two dominant, high-resolution techniques for the elucidation of complex subsurface structures.

Engineering geological investigation of Gololcha dam for evaluation of leakage and abutment slope stability, Eastern Ethiopia

The Ethiopian economy is primarily dependent on agriculture, making the construction of water harvesting facilities, such as dams, crucial for improving the productivity of this sector. The ongoing construction of the Gololcha dam on the Kurkura River located in Eastern Ethiopia aims to enhance irrigation schemes in the region. However, the dam site's complex geological and structural conditions pose challenges related to leakage and slope instability. Hence, this study focuses on addressing the abutment slope stability and leakage condition of this dam. This study employed kinematic analysis and the Limit Equilibrium Method (LEM) to assess slope stability. Additionally, engineering geological mapping, discontinuity surveys, seismic refraction tomography (SRT), and in-situ permeability testing were used to evaluate the leakage condition of the dam site. Notably, the permeability and SRT survey results identified potential leakage zones to the depth of 35, 30, and 35 m at the left, right, and central foundations of the dam, respectively. The kinematic method revealed one planar and two wedge modes of failure in the slope section covered by slightly weathered and fractured basalt rock at the right abutment. Further stability analysis of these two modes of failures via LEM analysis indicated slope instability under saturated conditions, emphasizing the role of pore water pressure. Furthermore, LEM modeling was directly utilized using the Slide 6.0 software to analyze the slope stability condition of the left abutment of the dam. This modeling also uncovered instability under saturated conditions. Based on the study findings, this study recommended curtain grouting to address potential leakage, as well as slope flattening and removing unstable rock wedges and loose material to stabilize unstable slope sections.

Cross-gradient joint inversion and clustering of ERT and SRT data on structured meshes incorporating topography

Summary The combination of electrical resistivity and seismic refraction tomography is a common practice for the characterization of subsurface features. Presently, the cross-gradient inversion scheme stands out as one of the most robust joint approaches, and some authors modified it to manage complex topographies on unstructured meshes even if at the expense of introducing additional parameters in the inversion process. We propose in this work a cross-gradient algorithm for jointly inverting electrical and seismic tomographic data on structured meshes in cases with non-flat topography. The proposed approach preserves the benefit of the classical cross-gradient approach without the need to impose physical length scales, as for irregular meshes. The quality of the results is evaluated in comparison with independent inversion through a new standardized cross-gradient index and a fuzzy c-means analysis that provides an assessment of the reconstruction accuracy through the membership function. The proposed method was applied to both synthetic models and field-scale examples located in Central Italy, where an accurate geophysical reconstruction is needed for the rehabilitation of existing dams. For all cases, joint inversion yielded superior results compared to independent inversion, demonstrating better agreement with available borehole data. The effectiveness of the joint approach was also demonstrated by the post-inversion tools, where the new cross-gradient index highlighted changes in structural similarity whilst fuzzy c-means clustering allowed for a quantitative reconstruction (position and shape) of the main units at the sites, facilitating the detection of site layering modifications.

A comparative study of using geophysical methods for imaging subsurface voids of various sizes and at different depths

Subsurface voids pose significant geohazards, underscoring the need for their timely detection in order to mitigate the associated hazard. We report on a field study aimed at the comparative assessment of electrical resistivity tomography (ERT), seismic refraction tomography (SRT), and the multichannel analysis of surface waves (MASW) for void mapping in karstic regions. The field surveys were conducted at a site in central Texas, of typical karstic geomorphology, and involved the co-located deployment of ERT, SRT and MASW arrays. Post-surveying, boreholes were drilled at select locations for verification purposes. It is shown that MASW demonstrated limited ability to resolve voids due to its inherent theoretical limitations. In contrast, ERT revealed high-resistivity air-filled zones, and low-resistivity soil-filled regions, which aligned well with post-survey borehole logs, although deeper voids remained largely undetected. SRT clearly delineated voids through velocity reductions, but smoothing effects overestimated void velocities. Using ERT and SRT jointly provided improved void characterization compared to single-method-based interpretations, with ERT determining void type and SRT delineating boundaries. Despite the relative success of the joint ERT-SRT application, it is evident that without the corroboration provided by invasive testing, definitive void localization and characterization under arbitrary site conditions remains elusive.

Minimum entropy constrained cooperative inversion of electrical resistivity, seismic and magnetic data

Geophysical methods are widely used to gather information about the subsurface as they are non-intrusive and comparably cheap. However, the solution to the geophysical inverse problem is inherently non-unique, which introduces considerable uncertainties. As a partial remedy to this problem, independently acquired geophysical data sets can be jointly inverted to reduce ambiguities in the resulting multi-physical subsurface images. A novel cooperative inversion approach with joint minimum entropy constraints is used to create more consistent multi-physical images with sharper boundaries with respect to the single-method inversions. Here, this approach is implemented in an open-source software and its applicability on electrical resistivity tomography (ERT), seismic refraction tomography (SRT), and magnetic data is investigated. A synthetic 2D ERT and SRT data study is used to demonstrate the approach and to investigate the influence of the governing parameters. The findings showcase the advantage of the joint minimum entropy (JME) stabilizer over separate, conventional smoothness-constrained inversions. The method is then used to analyze field data from Rockeskyller Kopf, Germany. 3D ERT and magnetic data are combined and the results confirm the expected volcanic diatreme structure with improved details. The multi-physical images of both methods are consistent in some regions, as similar boundaries are produced in the resulting models. Because of its sensitivity to hydrologic conditions in the subsurface, observations suggest that the ERT method senses different structures than the magnetic method. These structures in the ERT result do not seem to be enforced on the magnetic susceptibility distribution, showcasing the flexibility of the approach. Both investigations outline the importance of a suitable parameter and reference model selection for the performance of the approach and suggest careful parameter tests prior to the joint inversion. With proper settings, the JME inversion is a promising tool for geophysical imaging, however, this work also identifies some objectives for future studies and additional research to explore and optimize the method.

Evaluation of watertightness and slope stability analysis of Upper Guder dam, West Showa, Central Ethiopia

Ethiopia is currently emphasizing the construction of micro dams to enhance agricultural productivity and assure food security. The Upper Guder Dam, situated in the West Shoa Zone, Central Ethiopia, is one such project. However, the dam faces challenges due to complex geological and structural conditions, including leakage and slope instability. Hence, this study focuses on addressing the abutment slope stability and watertightness condition of this dam. Accordingly, the study employed discontinuity surveying, Seismic Refraction Tomography (SRT), Rock Quality Designation (RQD), and Lugeon testing to assess the water tightness. Additionally, kinematic analysis and the Limit Equilibrium Method (LEM) were used to evaluate abutment slope stability. Results from surface geological mapping and core drilling showed that the dam site is constituted by quaternary soil, tuff, and basalt. The Lugeon test results indicate that the left abutment of the dam is susceptible to leakage down to a depth of 40 m, which suggests that the permeable zone extends deeper than the results revealed by the SRT. Furthermore, this test demonstrated that the dam's right and central foundations are susceptible to leakage to the drilled depths. The SRT revealed that the central foundation, right, and left abutments of the dam are susceptible to leakage to the depth of 35 m, 30 m, and 34 m, respectively. Moreover, the kinematic analysis revealed that a section of the left abutment of the dam is susceptible to wedge mode of failure due to the intersection of JS1 and JS2. The LEM modeling of the right abutment of the dam also depicted that this section of the dam is unstable under saturated conditions which illustrates the importance of precipitation as the major slope destabilizing factor in the study area. Based on the study findings, this study recommended the use of curtain grouting to address the water tightness issue and slope angle reduction to mitigate the slope instability problem.

Brief communication: On the potential of seismic polarity reversal to identify a thin low-velocity layer above a high-velocity layer in ice-rich rock glaciers

Abstract. Seismic refraction tomography is a commonly used technique to characterise rock glaciers, as the boundary between unfrozen and ice-bearing layers represents a strong impedance contrast. In several rock glaciers, we observed a reversed polarity of the waves refracted by an extended ice-bearing layer compared to direct-wave arrivals. This phase change may be related to the presence of a thin low-velocity layer (LVL), such as fine- to coarse-grained sediments, above a thicker ice-rich layer. Our results are confirmed by the modelling and analysis of synthetic seismograms to demonstrate that the presence of a low-velocity layer can produce a polarity reversal on the seismic gather.

Geophysical and numerical approaches to solving the mechanisms of landslides triggered by earthquakes: A case study of Kahramanmaraş (6 February, 2023)

Many landslides occurred in the 11 Turkish provinces affected by the earthquakes centered in Kahramanmaraş in February 2023, and many people lost their lives due to those landslides. To reduce the risks and hazards brought on by landslides induced by earthquakes, it is essential to identify the mechanisms causing landslides to arise in such situations. In this study, the geological strength index and both seismic refraction tomography and electrical resistivity tomography (ERT) geophysical methods were used to determine the failure mechanism of layered sedimentary units defined as flysch (claystone-siltstone-sandstone-marl). The shear strength reduction finite element method (SSR-FEM) was used with the RS2 computer program for model solutions. One of the geophysical methods, ERT, facilitated an interpretation of the cyclic failure mechanism that occurs in flysch, and it was found to be congruent with numerical models. Landslide slip circle depth, water-saturated layers, and resistivity values were determined with the ERT method. Seismic methods were employed to determine the velocity values (Vp, Vs) of the landslip mass. These velocity data, along with the elasticity modulus and Poisson ratios derived from them, were utilized in our numerical model. In the maximum shear strength and horizontal displacement analyses performed on the section line determined before the landslide, the natural state and the situation under the influence of groundwater were examined and the strength reduction factor (SRF) was determined as 1.35 and 1.24, respectively. In an evaluation conducted under the influence of the maximum horizontal ground acceleration (0.22 g) determined by taking into account the seismicity of the study area, it was concluded that instability occurred in the flysches and the SRF value was 0.72. An important finding of this study is that the results of geophysical methods used to determine the lithological properties of sedimentary rock masses and understand the causes of landslides align closely with those of numerical analysis methods, providing high accuracy. These methods also support new approaches, particularly in identifying landslip boundaries and detecting areas with high water saturation. It has been discovered that earthquakes actively contribute to the instability of these kinds of soils.

Deciphering of karst geomorphology and sinkhole (doline) structures using multiple geophysical and geological methods (Trabzon, NE Türkiye)

To date, in Türkiye only a limited number and volume of combined geophysical and geological studies about karst have been performed. In this study, karstification and geomorphological features were examined with geophysical and geological methods together and initial results were obtained for Türkiye. Although the geology of the limestone forming the Şahinkaya Member, which contains Çayırbağı, Çalköy, and Çal Cave, near the Düzköy district of Trabzon/Türkiye province was studied by many researchers to date, there is no geophysical study to determine the internal structural features, groundwater, dolines, and karstic voids. The aim of this study was to identify karst formations and their structural extensions in Şahinkaya Member with geophysical methods. Therefore, three different study locations with a total surface area of approximately 3.2 km2 were examined with electrical resistivity tomography, self-potential, seismic refraction tomography, multichannel analysis of surface waves, and ground penetrating radar. These geophysical applications in limestone helped to identify karst cavities, water-saturated zones and dolines. Finally, the order of priority and efficiency of the five applied geophysical methods was compared, and the stages of the applications were outlined. In addition, the origin of karstification in the area investigated in this study was supported by petrographic, petrophysical and rock mechanic data.

A Comparison of Three Geophysical Methods for Determining the Depth to Bedrock

The location of the study site for the presented research here is at Illinois in USA. Geophysical investigations were conducted to characterize the subsurface. At the site, borehole control, downhole seismic (DHS), seismic refraction tomography (SRT) and multichannel analysis of surface waves (MASW) data were acquired for the purpose of seismic site characterization. Shear wave and compressional wave velocities were used to estimate depth to bedrock. The data analyses showed the depth to bedrock determined by DHS, MASW and SRT is in good agreement with borehole data.

Integration of ground geophysical methods to characterize near‐surface aquifer zones within an active mine

AbstractUnderstanding near‐surface groundwater storage, flow patterns, surface and groundwater interactions in mining areas can assist in making mining more efficient and profitable. This is especially important in opencast mines affected by water inflows that may negatively affect production and increase mining costs. We map and characterize the near‐surface aquifer zones at the opencast site of Tharisa Minerals, located in the southwestern region of the Bushveld Complex (South Africa). The main goal is to infer pit water inflow at the mine site and determine how it may be better controlled. The Bushveld Complex hosts partially connected and unconfined alluvial, shallow‐weathered and crystalline bedrock aquifers, which are often connected by small‐scale permeable zones. Seismic refraction tomography, multichannel analysis of surface waves, electrical resistivity tomography and borehole data are used to map and understand the different aquifer zones in the vicinity of the mine, as well as infer their relation to water inflow in the mine pits. The geophysical surveys map the overburden, weathered bedrock aquifer zone, and the top of the crystalline aquifer rock zone reasonably well. They reveal extensive and deep weathering, and possible high hydraulic conductivity in the vicinity of the mine. The results provide a better understanding of the mine's near‐surface environment, which could be used to implement effective and targeted dewatering techniques, thus enabling better pit inflow water control to improve mine working conditions and production.

Integrated of Seismic Refraction Tomography (SRT), and Multichannel Analysis of Surface Waves (MASW) for Verifying and Mapping a postulated Fault

The location of the study site for the presented paper here is at Arkansas in the United States of America. Geophysical investigations using seismic refraction tomography (SRT) and multichannel analysis of surface waves (MASW) surveys have been jointly carried out on the basis of an experimental investigation in the field. The motive of this study is to validate and map a postulated fault, in the eastern part of Arkansas in the United State of America. The SRT techniques have shown their potentiality in successfully identifying the fault and generated remarkably similar images of the fault, while the MASW tool generated a slightly different image of the fault due to condition that the earth under survey made up of layers of material that increase in seismic velocity which consider the perfect case for SRT techniques. The results show that the integrated use of seismic refraction tomography and multichannel analysis of surface waves is an effective tool with great potential. They provide very useful data for shallow mapping applications.

Assessment of different agricultural soil compaction levels using shallow seismic geophysical methods

A well-developed soil structure is crucial for a fertile soil and rules its influence on the ecosystem. Slight variations in the way soil components are organized and connected can greatly impact the soil mechanical and hydraulic characteristics. However, these features are challenging to measure at a scale that is relevant for agricultural management. In this study, we assess the capability of two seismic geophysical methods, i.e. Seismic Refraction Tomography and Multichannel Analysis of Surface Waves, to monitor the effects of compaction on agricultural bare soil. The purpose is to highlight the different mechanical response caused both by soil plastic deformation and soil water distribution due to increasing compaction. Results demonstrate that both geophysical techniques provide sufficient information to capture the effects of soil compaction and distinguish its significance, while traditional direct measurements, being punctual, lack sufficient spatial coverage. P-wave velocities carry a strong imprint of soil compaction, provided by seismic refraction, incorporating the information given by the multiphase medium in its entirety. On the other hand, S-wave velocity derived from Surface Waves discriminates the effect of solid matrix structure. This work, moreover, aims to pave the way for seismic methods to spatially characterize compaction at field scale, illustrating its potential and suggesting possible developments.

Assessment of Site Effects and Numerical Modeling of Seismic Ground Motion to Support Seismic Microzonation of Dushanbe City, Tajikistan

In the territory of Dushanbe city, the capital of Tajikistan, detailed geological and geophysical data were collected during geophysical surveys in 2019–2020. The data comprise 5 microtremor array measurements, 9 seismic refraction tomography profiles, seismological data from 5 temporary seismic stations for standard spectral ratio calculations, 60 borehole datasets, and 175 ambient noise measurements. The complete dataset for Dushanbe was used to build a consistent 3D geologic model of the city with a size of 12 × 12 km2. The results of the seismological and geophysical surveys were compared and calibrated with borehole data to define the boundaries of each layer in the study area. The Leapfrog Works software was utilized to create a 3D geomodel. From the 3D geomodel, we extracted six 12 km long 2D geological cross-sections. These 2D geological cross-sections were used for 2D dynamic numerical modeling with the Universal Distinct Element Code software to calculate the local seismic response. Finally, the dynamic numerical modeling results were compared with the amplification functions obtained from the seismological and ambient noise data analysis. The 2D dynamic numerical modeling results allowed a better assessment of the site effects in the study area to support seismic microzonation and the determination of local peak ground acceleration changes in combination with regional seismic hazard maps. In addition, our results confirm the strong seismic amplification effects noted in some previous studies, which are attributed to the influence of local topographic and subsurface characteristics on seismic ground motions.

Seismic characterisation of the subsoil under a historic building: Cathedral Church of Saint Mary in Murcia case study

This research focuses on the Cathedral Church of Saint Mary in Murcia, Spain, which is located in a moderate-to-high seismic risk zone in the Spanish context. The study uses geophysical techniques and geotechnical investigation to characterise seismic ground models at the building's scale, aiming to understand the real amplification of the ground under seismic effect in historic buildings. Three soil layers (silt, sand with gravel, and gravel) were identified through core borings. Multichannel Analysis of Surface Waves (MASW) profiles and Mini-array profiles revealed Vs values of 305 ± 32 m/s, 296 ± 62 m/s, and 440 ± 38.5 m/s, respectively, for these materials. Seismic Refraction Tomography (SRT) showed Vp values of 586 ± 73 m/s, 700 m/s, and 1466 ± 489 m/s for the corresponding layers. The horizontal-to-vertical spectral ratio (HVSR) approach identified a ground predominant period ranging from 0.37 to 0.38 s. Another significant peak at 2.3 s is observed, probably associated to the Triassic basement. Three seismic events with magnitudes Mw between 4.9 and 5.1 were used as inputs to determine the amplification factor (AF). The results indicate a PGA amplification factor between 1.7 and 2.1. These results contribute to the conservation and mitigation of seismic risk of this cultural heritage generating input data that enables precise computation of the soil-structure interaction.

Assess the Effectiveness of Multichannel Analysis of Surface Waves Method in Mapping Ancient Structures in Ultrashallow Aquatic Environments: The Case of Agioi Theodoroi, Greece

ABSTRACTThe present geophysical research aims to evaluate the applicability of multichannel analysis of surface waves (MASW) on mapping ultrashallow underwater ancient masonry remnants. The work presents the analysis from a single seismic line using MASW and seismic refraction tomography (SRT) methods and its corresponding electrical resistivity tomography (ERT) section surveyed at the submerged prehistoric site of Agioi Theodoroi area located 10 km eastern of Heraklion, Crete, Greece. The 2D MASW velocity model exhibits significant correspondence with the resistivity structure extracted from the ERT data, showing lateral S‐wave velocity (Vs) variations at the positions where the high resistivity anomalies exist. The analysis of synthetic seismic data calculated from a respective model reproduced a comparable S‐wave velocity pseudo‐section with the real data. However, the investigated targets (submerged buried masonry) appear shallower and wider in MASW sections than in the real world and the corresponding synthetic models, due to insufficient vertical and horizontal resolution of this method. Surface waves travelling through the seafloor sediments (Scholte‐waves) demonstrate very low velocity values. This makes them suitable for the detection of shallow and relatively large (> 0.5 m) underwater manmade structures, providing the enhancement of MASW method resolution, by utilizing a high frequency (> 100 Hz) seismic source, recording short Scholte wavelengths (≤ 1 m) and using shorter (≤ 0.5 m) receiver spacing and array length. Consequently, the results of this work demonstrate the potential in employing conventional seismic techniques in the delineation of underwater antiquities and the revealing of the cultural dynamics in very shallow off‐shore archaeological sites.

Factors Affecting the Stability of Loess Landslides: A Review

The stability of loess landslides affects the production and livelihood of the people in its vicinity. The stability of loess landslides is influenced by various factors, including internal structure, collapsibility, water content, and shear strength. The landslide stability of loesses can be analyzed by several geophysical methods, such as seismic refraction tomography (SRT), electrical resistivity tomography (ERT), micro-seismic technology, and ground penetrating radar (GPR). Geotechnical tests (compression and shear tests) and remote sensing techniques (Global Navigation Satellite System (GNSS), Interferometric Synthetic Aperture Radar (InSAR) and airborne 3D laser technology) are used for studying the landslide stability of loesses as well. Some of the methods above can measure parameters (e.g., fractures, water content, shear strength, creep) which influence the stability of loess landslides, while other methods qualitatively indicate the influencing factors. Integrating parameters measured by different methods, minimizing disturbances to landslides, and assessing landslide stability are important steps in studying landslide hazards. This paper comprehensively introduces the methods used in recent studies on the landslide stability of loesses and summarizes the factors which affect the landslide stability. Furthermore, the relationships between different parameters and methods are examined. This paper enhances comprehension of the underlying mechanisms of the stability of loess landslides to diminish disastrous consequences.

Computation and Modelling of Seismic Refraction Tomography for Anisotropic Medium Based on Pseudo Bending Method

Seismic refraction tomography is one of the imaging techniques in geophysical methods used to remodel the near-surface velocity layer structure of the Earth. In this study, we carried out a new computational approach and modelling of seismic refraction tomography using the pseudo-bending method. The true model of the near-surface is designed to be anisotropic medium which is having a low velocity anomaly distribution. This anomaly is constructed in such a way as to be similar to a model of liquid waste away spreading, which exhibits seismic velocities ranging from 1600 m/s to 1800 m/s. Based on our computations and numerical modelling results, it was found that the ray tracing path using pseudo-bending method displays an asymmetrical trajectory when the positions of the source and geophone are exchanged. Altering the shooting configuration from direct shoot (DS) to reversed shoot (RS) also reveals a significant difference in travel time values. The results of delay time tomography inversion, which represents the difference between travel times in the true model and the initial model using the SIRT method, indicate the presence of a low velocity anomaly that can be interpreted as the distribution of liquid waste.

Impacts of hydrofacies geometry designed from seismic refraction tomography on estimated hydrogeophysical variables

Abstract. Understanding the critical zone processes related to groundwater flows relies on subsurface structure knowledge and its associated parameters. We propose a methodology to draw the patterns of the subsurface critical zone at the catchment scale from seismic refraction data and show its interest for hydrological modelling. The designed patterns define the structure of a physically based distributed hydrological model applied to a mountainous catchment. With that goal, we acquired 10 seismic profiles covering the different geomorphology zones of the studied catchment. We develop a methodology to analyse the geostatistical characteristics of the seismic data and interpolate them over the whole catchment. The applied geostatistical model considers the scale variability of the subsurface structures observed from the seismic data analysis. We use compressional seismic wave velocity thresholds to identify the depth of the soil and saprolite bottom boundaries. Assuming that such porous compartments host the main part of the active aquifer, their patterns are embedded in a distributed hydrological model. We examine the sensitivity of classical hydrological data (piezometric heads) and geophysical data (magnetic resonance soundings) to the applied velocity thresholds used to define the soil and saprolite boundaries. Different sets of hydrogeological parameters are used in order to distinguish general trends or specificities related to the choice of parameter values. The application of the methodology to an actual catchment illustrates the interest of seismic refraction in constraining the structure of the critical zone subsurface compartments. The sensitivity tests highlight the complementarity of the analysed hydrogeophysical data sets.