Hydrogeophysics and remote sensing for the design of hydrogeological conceptual models in hard rocks – Sardón catchment (Spain)

Abstract

Abstract Hard rock aquifers are highly heterogeneous and hydrogeologically complex. To contribute to the design of hydrogeological conceptual models of hard rock aquifers, we propose a multi-techniques methodology based on a downward approach that combines remote sensing (RS), non-invasive hydrogeophysics and hydrogeological field data acquisition. The proposed methodology is particularly suitable for data scarce areas. It was applied in the pilot research area of Sardon catchment (80 km2) located west of Salamanca (Spain). The area was selected because of hard-rock hydrogeology, semi-arid climate and scarcity of groundwater resources. The proposed methodology consisted of three main steps. First, we detected the main hydrogeological features at the catchment scale by processing: (i) a high resolution digital terrain model to map lineaments and to outline fault zones; and (ii) high-resolution, multispectral satellite QuickBird and WorldView-2 images to map the outcropping granite. Second, we characterized at the local scale the hydrogeological features identified at step one with: i) ground penetrating radar (GPR) to assess groundwater table depth complementing the available monitoring network data; ii) 2D electric resistivity tomography (ERT) and frequency domain electromagnetic (FDEM) to retrieve the hydrostratigraphy along selected survey transects; iii) magnetic resonance soundings (MRS) to retrieve the hydrostratigraphy and aquifer parameters at the selected survey sites. In the third step, we drilled 5 boreholes (25 to 48 m deep) and performed slug tests to verify the hydrogeophysical interpretation and to calibrate the MRS parameters. Finally, we compiled and integrated all acquired data to define the geometry and parameters of the Sardon aquifer at the catchment scale. In line with a general conceptual model of hard rock aquifers, we identified two main hydrostratigraphic layers: a saprolite layer and a fissured layer. Both layers were intersected and drained by fault zones that control the hydrogeology of the catchment. The spatial discontinuities of the saprolite layer were well defined by RS techniques while subsurface geometry and aquifer parameters by hydrogeophysics. The GPR method was able to detect shallow water table at depth between 1 and 3 m b.g.s. The hydrostratigraphy and parameterization of the fissured layer remained uncertain because ERT and FDEM geophysical methods were quantitatively not conclusive while MRS detectability was restricted by low volumetric water content. The proposed multi-technique methodology integrating cost efficient RS, hydrogeophysics and hydrogeological field investigations allowed us to characterize geometrically and parametrically the Sardon hard rock aquifer system, facilitating the design of hydrogeological conceptual model of the area.