Non-invasive Geophysical Surveys in Search of the Roman Temple of Augustus Under the Cathedral of Tarragona (Catalonia,Spain): A Case Study

An integrated geophysical survey has been conducted at the Tarragona’s Cathedral (Catalonia, NE Spain) with the aim to confirm the potential occurrence of archaeological remains of the Roman Temple dedicated to the Emperor Augustus. Many hypotheses have been proposed about its possible location, the last ones regarding the inner part of the Cathedral, which is one of the most renowned temples of Spain (twelfth century) evolving from Romanesque to Gothic styles. A geophysical project including electrical resistivity tomography (ERT) and ground probing radar (GPR) was planned over 1 year considering the administrative and logistic difficulties of such a project inside a cathedral of religious veneration. Finally, both ERT and GPR have been conducted during a week of intensive overnight surveys that provided detailed information on subsurface existing structures. The ERT method has been applied using different techniques and arrays, ranging from standard Wenner–Schlumberger 2D sections to full 3D electrical imaging with the advanced Maximum Yield Grid array. Electrical resistivity data were recorded extensively, making available many thousands of apparent resistivity data to obtain a complete 3D image after a full inversion. In conclusion, some significant buried structures have been revealed providing conclusive information for archaeologists. GPR results provided additional information about shallowest structures. The geophysical results were clear enough to persuade religious authorities and archaeologists to conduct selected excavations in the most promising areas that confirmed the interpretation of geophysical data. In conclusion, the significant buried structures revealed by geophysical methods under the cathedral were confirmed by archaeological digging as the basement of the impressive Roman Temple that headed the Provincial Forum of Tarraco, seat of the Concilium of Hispania Citerior Province.     

Buried Quaternary valleys in western Denmark—occurrence and inferred implications for groundwater resources and vulnerability

Numerous geophysical investigations in the western part of onshore Denmark constitute the basis for a delineation of buried Quaternary valleys. The geophysical methods comprise primarily Transient ElectroMagnetic (TEM) and reflection seismic surveys, and the geophysical data have been combined with lithological data from boreholes. Buried valleys appear both as single valleys and in dense networks. The internal structure of the valleys is typically complex due to repeated erosional and depositional events. Buried valleys are common geological structures in the region and they influence the distribution of Tertiary and Quaternary sediments greatly. A large number of buried valleys in the region contain important aquifers, whose natural protection varies depending on thickness and character of overlying clay layers. Many of these aquifers are deep-seated and well protected, but because of the prevailing heterogeneity of the valley infill and the erosional incisions created by different valley generations, preferential flow paths for downward transport of contaminated water from shallow aquifers may occur.     

Geophysical study of the Peinan Archaeological Site,Taiwan

The Peinan archaeological site is the most intact Neolithic village with slate coffin burial complexes in Taiwan. However, the area that potentially contains significant ancient remains is covered by dense vegetation. No reliable data show the distribution of the ancient village, and no geophysical investigation has been performed at this site. To evaluate various geophysical methods under the geological setting and surface condition of the site, the physical properties of the remains were measured and four geophysical methods involving magnetic, electromagnetic (EM), electrical resistivity tomography (ERT), and ground-penetrating radar (GPR) were tested along three parallel profiles. The results imply that the EM and magnetic methods are much cost-effective and suitable for investigating the entire area. GPR and ERT methods can provide high resolution subsurface image, which are much suitable for subsequently detail investigation.The EM and magnetic surveys were thus conducted over the entire Peinan Cultural Park to understand the distribution of the ancient building remains at the Peinan site. The results of this study were verified by subsequent excavations, which indicate that the EM survey was successful in delineating the majority of the ancient village because the basements of building are highly resistive in comparison to the background sediment. The results of this investigation suggest that the ancient village was broadly distributed over the eastern part of the Peinan Culture Park and extended to the southeast.     

Assessment of an ancient bridge combining geophysical and advanced photogrammetric methods: Application to the Pont De Coq,France

A high resolution geophysical survey was carried out on the Pont De Coq, a medieval stone arch bridge located in Normandy (France) in 2011 and 2012. Two complementary methods are used: Electrical Resistivity Tomography (ERT) and Ground PenetratingRadar (GPR). They allow to evaluate the structural state of the bridge and to characterize the subsurface around and beneath the bridge. An excellent correlation is obtained between the geophysical methods and the geological data obtained around the bridge. In order to improve the restitution of the geophysical data, an advanced photogrammetric method is performed, providing a high resolution 3D Digital Terrain Model (DTM) of the Pont de Coq. The advanced photogrammetry enhances the presentation of the GPR and ERT data. This approach is an easy-to-use, rapid and cost-effective tool for stakeholders. Finally, it is a promising and original method for improved interpretations of future geophysical surveys.     

Imaging the Mariánské Lázně Fault (Czech Republic) by 3-D ground-penetrating radar and electric resistivity tomography

Geodynamic activity in the area of West Bohemia is typified by the occurrence of earthquake swarms, Quaternary volcanism and high flux of mantle-derived CO₂. The highest swarm activity occurs beneath the eastern edge of the Cheb basin, which is delineated by the NW-SE trending morphologically pronounced Mariánské Lázn?? Fault (MLF) controlling the formation of the basin. The previous trenching survey across the MLF zone has identified several fault strands with possible Quaternary activity. In this paper we present the results of the geophysical survey focused to trace the faults signatures in geophysical sections and to build an image of near surface tectonics. The method of electric resistivity tomography (ERT) along two profiles parallel to the trench identified a strong resistivity contrast between the bodies of sandy gravels in the middle and conductive clayey sands to the west and weathered crystalline basement to the east. The 2-D ground penetration radar (GPR) sections show direct correlation of reflections with lithological boundaries identified in the trench. As expected, the GPR signal amplitudes increase with the resistivities found in the ERT sections. Two of the four faults identified in the trench are indicated in the resistivity and GPR sections. A 3-D GPR measurement has identified a spot of high amplitudes elongated parallel to the MLF trend, which coincides with the high resistivity body. To improve the signal-to-noise ratio of the time slices we stacked the GPR time slices within vertically homogeneous blocks. This provided a contrast image of the sand-gravel body including its boundaries in three dimensions. The detailed analysis of the 3-D GPR cube revealed additional fault that limits the highly reflective sands and appears to be offset by another younger fault. Our results suggest a complex fault pattern in the studied area, which deserves a further study.     

Using geophysical surveys to test tracer‐based storage estimates in headwater catchments

Hydrogeophysical surveys were carried out in a 3.2 km² Scottish catchment where previous isotope studies inferred significant groundwater storage that makes important contributions to streamflow. We used electrical resistivity tomography (ERT) to characterize the architecture of glacial drifts and make an approximation of catchment‐scale storage. Four ERT lines (360–535 m in length) revealed extensive 5–10 m deep drift cover on steeper slopes, which extends up to 20–40 m in valley bottom areas. Assuming low clay fractions, we interpret variable resistivity as correlating with variations in porosity and water content. Using Archie's Law as a first approximation, we compute likely bounds for storage along the ERT transects. Areas of highest groundwater storage occur in valley bottom peat soils (up to 4 m deep) and underlying drift where up to 10 000 mm of precipitation equivalent may be stored. This is consistent with groundwater levels which indicate saturation to within 0.2 m of the surface. However, significant slow groundwater flow paths occur in the shallower drifts on steeper hillslopes, where point storage varies between ~1000 mm–5000 mm. These fluxes maintain saturated conditions in the valley bottom and are recharged from drift‐free areas on the catchment interfluves. The surveys indicate that catchment scale storage is >2000 mm which is consistent with tracer‐based estimates. Copyright © 2016 John Wiley & Sons, Ltd.     

Combining Electrical Resistivity Tomography and Ground Penetrating Radar to study geological structuring of karst Unsaturated Zone

This paper highlights the efficiency and complementarity of a light package of geophysical techniques to study the structure of karst Unsaturated Zone (UZ) in typical Mediterranean environment where soil cover is thin or absent. Both selected techniques, 2D Ground Penetrating Radar (GPR) and Electrical Resistivity Tomography (ERT), are widely used in environmental studies and their application is accessible for a lot of scientists/engineers. However, GPR or ERT alone is not able to provide an enhanced characterization of geological features in karst media. In the present study, GPR results supply a near surface high resolution imaging and thus can provide relevant geological information such as stratifications and fractures. Despite the quality of the results GPR's investigation depth remains limited to around 12 m. Apparent and inverted resistivity provided by ERT surveys shows strong lateral and vertical variations. These variations can inform about general geological structuring and feature orientation. ERT is able to prospect down to 40 m but it's a low resolution integrative technique. In the study area the investigated limestone is a commonly electrical resistive formation (more than 2000 Ω.m). However deeper than 5–7 m, the ERT profiles reveal several zones of moderate resistivity (around 900 Ω.m). In these zones a stratification change corresponding to slanted bedding is clearly identified by GPR results. The combination of both GPR and ERT results can allow a well-established geological interpretation. These moderate resistivity zones with slanted beddings can explain the presence of a perennial water flow point 35 m below the surface of the studied site within the underground gallery of the Low-Noise Underground Laboratory (LSBB).     

Evaluation of trenching,ground penetrating radar (GPR) and electrical resistivity tomography (ERT) for sinkhole characterization

This paper explores the suitability and advantages of combining the trenching technique with geophysical surveys [ground penetrating radar (GPR) and electrical resistivity tomography (ERT)] for sinkhole characterization in a mantled karst area. The approach is applied to two active sinkholes concealed by anthropogenic deposits and formed by contrasting subsidence mechanisms; collapse and sagging. The ERT section acquired across the collapse sinkhole images the clayey fill of the depression as an obvious low resistivity area, showing the approximate location of the sinkhole edges. Spatially dense GPR surveys provide information on the position of the boundaries of the concealed subsidence structures and their three‐dimensional (3D) internal geometry, revealing the dominant subsidence mechanism. We illustrate the impact of several factors on the quality of the GPR data such as sinkhole size, nominal frequency of the antennas, antenna shielding, and the presence of backfilled excavations and above‐surface objects. Trenches provided detailed information on the subsurface structure of the sinkhole, subsidence magnitude, partitioning of the strain, and the position of the sinkhole edges, especially when they are deep enough and excavated across the central sector and perpendicular to the boundaries. The stratigraphic and structural relationships observed in the trench were then used to infer the spatial evolution of the sinkholes (e.g. enlargement), their kinematic behavior (episodic versus progressive), and to differentiate discrete subsidence events and their associated magnitude. Numerical dates were used to estimate average subsidence rates and the recurrence of subsidence events. Such integrated data sets may be used as an objective basis to forecast the future behavior of potentially damaging sinkholes and to assess the associated hazard and risk. Copyright © 2013 John Wiley & Sons, Ltd.     

Geophysical investigation of the Sandalp rock avalanche deposits

In the study of rock avalanche phenomena, numerical modelling makes use of back analyses of the rock avalanche propagation for calibration of the modelling assumptions and parameters. The back analyses require knowledge of the run-out area boundaries and the thickness distribution of the deposit. Geophysical methods can be applied to retrieve the thickness distribution, but, due to strong heterogeneities and logistic problems they are seldom applied. The aim of this work is to assess the potential of integrated geophysical methods to recognise and characterise a deposit created by two rock avalanches which occurred in the Sandalp valley (Switzerland) in 1996. The topography of the site before and after the rock avalanche is known and can be used as a benchmark. Resistivity tomography, seismic P-wave tomography, and active and passive surface wave analysis have been applied on several profiles deployed both on the rock avalanche deposit and in the surrounding area. Innovative approaches for surface wave analysis based on laterally constrained inversion and multimodal inversion have been applied to the data. A comparison of the results of the geophysical investigations with the topographic benchmark has shown the capability of the geophysical methods to locate the bottom of the deposit in the areas where the contrast with the host sediments properties is significant. In these areas, the deposit has higher resistivities and lower seismic velocities than the underlying materials. In the areas where the deposit is thicker and richer in fine-grained materials the geophysical parameters are not able to discriminate between the rock avalanche deposit and the underlying sediments. As a secondary task, the geophysical methods also allowed the bedrock pattern to be outlined.     

Integrated Electrical Resistivity Tomography and Ground Penetrating Radar Measurements Applied to Tomb Detection

The integrated use of electrical resistivity tomography (ERT) and ground penetrating radar (GPR) measurements, and in particular the joint analysis of 2D and 3D data, can represent a valid solution for target identification at complex archaeological sites. A good example, in this respect, is given by the case study of a Phoenician–Punic necropolis in the archaeological site of Nora, in southern Sardinia (Italy), where GPR and ERT measurements were collected before site excavation. In this specific case, the mix of soil and air in the buried chambers, as well as the orientation and the complex spatial distribution of these structures into the sandstone bedrock, generated a number of anomalies difficult to interpret only using 2D results. Only the integration of all GPR and ERT data in a 3D view, and the comparison with archaeological evidence after the excavation, allowed a solid interpretation of geophysical anomalies visible in the 2D sections. Overall, this case study demonstrates the efficiency of the combined use of GPR and ERT acquisitions and shows how, in general, only the joint analysis of 2D data and in a 3D view can help the interpretation of the real distribution of the buried archaeological remains at similar archaeological complex sites.     

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework.The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy–sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements.The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux.As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers.     

EM and GPR investigations of contaminant spread around the Hoc Mon waste site,Vietnam

This paper is presenting the results from near-surface geophysical surveys near the waste site of Hoc Mon in southern Vietnam where leachate contamination has been recognized at the surface. Using EM and GPR surveys, we were able to determine the lateral extent of a contaminated area of high electrical conductivity and have identified channels that concentrate the contaminant flow. The simple relationship between the electrical resistivity and the leachate concentration is suggested and estimated the in situ leachate concentration from the inversion of the EM data; values are as high as 40%. Thanks to a permeability barrier leachate flow is confined to the shallow subsurface, making it easier to apply possible site remediation projects.     

Characterizing hydrological processes in a semiarid rangeland watershed: A hydrogeophysical approach

The complex ecohydrological processes of rangelands can be studied through the framework of ecological sites (ESs) or hillslope‐scale soil–vegetation complexes. High‐quality hydrologic field investigations are needed to quantitatively link ES characteristics to hydrologic function. Geophysical tools are useful in this context because they provide valuable information about the subsurface at appropriate spatial scales. We conducted 20 field experiments in which we deployed time‐lapse electrical resistivity tomography (ERT), variable intensity rainfall simulation, ground‐penetrating radar (GPR), and seismic refraction, on hillslope plots at five different ESs within the Upper Crow Creek Watershed in south‐east Wyoming. Surface runoff was measured using a precalibrated flume. Infiltration data from the rainfall simulations, coupled with site‐specific resistivity–water content relationships and ERT datasets, were used to spatially and temporally track the progression of the wetting front. First‐order constraints on subsurface structure were made at each ES using the geophysical methods. Sites ranged from infiltrating 100% of applied rainfall to infiltrating less than 60%. Analysis of covariance results indicated significant differences in the rate of wetting front progression, ranging from 0.346 m min^?1/2 for sites with a subsurface dominated by saprolitic material to 0.156 m min^?1/2 for sites with a well‐developed soil profile. There was broad agreement in subsurface structure between the geophysical methods with GPR typically providing the most detail. Joint interpretation of the geophysics showed that subsurface features such as soil layer thickness and the location of subsurface obstructions such as granite corestones and material boundaries had a large effect on the rate of infiltration and subsurface flow processes. These features identified through the geophysics varied significantly by ES. By linking surface hydrologic information from the rainfall simulations with subsurface information provided by the geophysics, we can characterize the ES‐specific hydrologic response. Both surface and subsurface flow processes differed among sites and are directly linked to measured characteristics.     

Geophysical Discovery of a New LNAPL Plume at the Former Wurtsmith AFB, Oscoda, Michigan

Alight nonaqueous phase liquid (LNAPL) ground water contaminant plume has been discovered by purely geophysical means at the former Wurtsmith Air Force Base (AFB) near Oscoda, Michigan. It is located near another plume called FT-02, which is a well-studied area undergoing natural bioremediation. The plume was discovered by ground penetrating radar (GPR) profiling while extending a long line from FT-02 to establish background variability around that plume. The new plume was apparent because of a high-conductivity "shadow' or GPR reflection attenuation observed below the conductive zone at the top of the aquifer, identical to the pattern observed at the FT-02 plume. Further GPR surveys were conducted by students of a Western Michigan University geophysics field course to outline the proximal part of the plume. The GPR survey was supplemented by an electromagnetic induction (EM) survey which showed a group of four cables crossing the area. Finally, a magnetometer survey was conducted to search for any buried steel objects which might have been missed by the EM survey. The results of the three geophysical surveys were then used by students of a University of Michigan field course to guide subsurface soil and fluid sampling, which verified the presence of residual LNAPL product and ground water with conductivities 2.5 to 3.3 times above background. The plume source is in the vicinity of a vaulted underground storage tank (UST) formerly used for the collection of waste solvents and fuels for subsequent use in the fire training exercises at FT-02. This newly discovered LNAPL plume, along with other "mature' plumes, fits the electrical model which predicts conductive ground water below the decomposing but electrically resistive LNAPLs. Finally, this is a fine example of the cooperative use of a dedicated research site for training by students of two different universities.     

Postglacial sediment budget of Chilliwack Valley,British Columbia

The paraglacial reworking of glacial sediments by rivers and mass wasting is an important conditioning factor for modern sediment yields in mountainous catchments in formerly glaciated regions. Catchment scale and patterns of sediment storage are important influences in the rate of postglacial adjustment. We develop a quantitative framework to estimate the volume, sediment type, and fractional size distribution of legacy glacial materials in a large (1230 km²) watershed in the North Cascade Mountains in south‐western British Columbia, Canada. Chilliwack Valley is exceptional because of the well‐dated bounds of deglaciation. Interpolation of paleo‐surfaces from partially eroded deposits in the valley allows us to estimate the total evacuated sediment volume. We present a chronology of sediment evacuation from the valley and deposition in the outlet fan, based on infrared stimulated luminescence (IRSL) and ¹⁴ C dating of river terraces and fan strata, respectively. The effects of paraglacial sedimentation in Chilliwack Valley were intensified through a major fall in valley base‐level following ice retreat. The steepened mainstem valley gradient led to deep incision of valley fills and fan deposits in the lower valley network. The results of this integrated study provide a postglacial chronology and detailed sediment budget, accounting for long‐term sorting of the original sediments, lag deposit formation in the mainstem, deposition in the outlet fan, and approximate downstream losses of suspended sediment and wash load. The mass balance indicates that a bulk volume of approximately 3.2 km³ of glacial material has been evacuated from the valley. Copyright © 2012 John Wiley & Sons, Ltd.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

Geoelectrical and hydrochemical investigations for characterizing the salt water intrusion in the Khanasser valley, northern Syria

An integrated approach of geoelectrical and hydrochemical investigation surveys was proposed for indicating contact regions between saline and fresh groundwater in the Khanasser valley region, northern Syria. The qualitative and quantitative interpretations of 34 vertical electrical soundings (VES) enable to characterize the salt water intrusion laterally and vertically. The established iso-apparent resistivity maps for different AB/2 spacings obviously indicate the presence of a lowresistivity (less than 4 Ohm·m) zone related to the salt water intrusion in the Quaternary and Paleogene deposits. The different hydrochemical and geophysical parameters, such as electrical resistivity, total dissolved solids (TDS) and major ions concentrations used to characterize the salt water intrusion gave almost similar results in locating and mapping the different boundaries of the groundwater salinity. The proposed approach is useful for mapping the interface between different groundwater qualities, and can be therefore used to successfully characterize the salt water intrusion phenomenon in other semi-arid regions. The application of such an approach is a powerful tool and can be used for water resource management in the water scarce areas.     

Time‐Lapse Electrical Geophysical Monitoring of Amendment‐Based Biostimulation

Biostimulation is increasingly used to accelerate microbial remediation of recalcitrant groundwater contaminants. Effective application of biostimulation requires successful emplacement of amendment in the contaminant target zone. Verification of remediation performance requires postemplacement assessment and contaminant monitoring. Sampling‐based approaches are expensive and provide low‐density spatial and temporal information. Time‐lapse electrical resistivity tomography (ERT) is an effective geophysical method for determining temporal changes in subsurface electrical conductivity. Because remedial amendments and biostimulation‐related biogeochemical processes often change subsurface electrical conductivity, ERT can complement and enhance sampling‐based approaches for assessing emplacement and monitoring biostimulation‐based remediation. Field studies demonstrating the ability of time‐lapse ERT to monitor amendment emplacement and behavior were performed during a biostimulation remediation effort conducted at the Department of Defense Reutilization and Marketing Office (DRMO) Yard, in Brandywine, Maryland, United States. Geochemical fluid sampling was used to calibrate a petrophysical relation in order to predict groundwater indicators of amendment distribution. The petrophysical relations were field validated by comparing predictions to sequestered fluid sample results, thus demonstrating the potential of electrical geophysics for quantitative assessment of amendment‐related geochemical properties. Crosshole radar zero‐offset profile and borehole geophysical logging were also performed to augment the data set and validate interpretation. In addition to delineating amendment transport in the first 10 months after emplacement, the time‐lapse ERT results show later changes in bulk electrical properties interpreted as mineral precipitation. Results support the use of more cost‐effective surface‐based ERT in conjunction with limited field sampling to improve spatial and temporal monitoring of amendment emplacement and remediation performance.