3D Ground Penetrating Radar Survey Research Articles

Ground penetrating radar investigations of buried remnants at ancient capital cities of Panchala and Vatsa kingdoms spread along Ganga-Yamuna doab of India from 600 BCE to 1100 CE

Ahichhatra and Kaushambi, ancient settlements located in Uttar Pradesh, comprise differential cultural deposits and architectural remnants from 2500 BCE to 1200 CE. These settlements flourished as the capital cities of once-prosperous urban kingdoms known as Panchala and Vatsa Mahajanapadas but remain abandoned in the present day. To provide ancillary data for the simplification of traditional archaeological explorational practices, a few portions located in and around the fortified areas of both settlements were considered for subsurface investigations using Ground Penetrating Radar (GPR). In Ahichhatra and Kaushambi, 2D and 3D GPR survey was carried out in the open areas devoid of vegetation for the understanding of the spatial alignment of the buried features within the fortified walls. The subsurface investigations at Ahichhatra have shown high amplitude reflections referring to buried wall structures, residential complexes, raised platforms, past-used surfaces and structural continuity. Also, the partial excavation carried out at Ahichhatra has revealed walls and compacted floors validating our findings. From Kaushambi, the findings include wall structures and interconnected chambered structures. The anomalies obtained in the data collected from considered survey sites located in Ahichhatra are 1–2 m wide and 2–4 m deep whereas in Kaushambi they are 0.5–1 m wide and 2–4 m deep. These findings helped us to determine the progressive urban relics of the ancient kingdoms which were known for their second phase of urbanism in Ganga Valley that prevailed until their abandonment. For a better understanding of the cultural history, detailed mapping of the contemporary sites can be carried out using integrated survey methods.

Integration of Terrestrial Laser Scanner (TLS) and Ground Penetrating Radar (GPR) to Characterize the Three-Dimensional (3D) Geometry of the Maoyaba Segment of the Litang Fault, Southeastern Tibetan Plateau

High-resolution topographic and stratigraphic datasets have been increasing applied in active fault investigation and seismic hazard assessment. There is a need for the comprehensive analysis of active faults on the basis of the correlating geomorphologic features and stratigraphic data. The integration of TLS and GPR was adopted to characterize the 3D geometry of the fault on the Maoyaba segment of Litang fault. The TLS was used to obtain the high-resolution topographic data for establishing the 3D surficial model of the fault. The 2D 250 MHz and 500 MHz GPR profiles were carried out to image the shallow geometry of the fault along four survey lines. In addition, the 3D GPR survey was performed by ten 2D 500 MHz GPR profiles with 1 m spacing. From the 2D and 3D GPR results, a wedge-shaped deformation zone of the electromagnetic wave was clearly found on the GPR profiles, and it was considered to be the main fault zone with a small graben structure. Three faults were identified on the main fault zone, and fault F1 and F3 were the boundary faults, while the fault F2 was the secondary fault. The subsurface geometry of the fault on the GPR interpreted results is consistent with the geomorphologic features of the TLS-derived data, and it indicates that the Maoyaba fault is a typical, normal fault. For reducing the environmental disruption and economic losses, GPR was the most optimal method for detecting the subsurface structures of active faults in the Litang fault with a non-destructive and cost-effective fashion. The 3D surface and subsurface geometry of the fault was interpreted from the integrated data of TLS and GPR. The fusion data also offers the chance for the subsurface structures of active faults on the GPR profiles to be better understood with its corresponding superficial features. The study results demonstrate that the integration of TLS and GPR has the capability to obtain the high-resolution micro geomorphology and shallow geometry of active faults on the Maoyaba segment of the Litang fault, and it also provides a future prospect for the integration of TLS and GPR, and is valuable for active fault investigation and seismic hazard assessment, especially in the Qinghai-Tibet Plateau area.

Automatic modelling of urban subsurface with ground-penetrating radar using multi-agent classification method

ABSTRACT The subsurface of urban cities is becoming increasingly congested. In-time records of subsurface structures are of vital importance for the maintenance and management of urban infrastructure beneath or above the ground. Ground-penetrating radar (GPR) is a nondestructive testing method that can survey and image the subsurface without excavation. However, the interpretation of GPR relies on the operator’s experience. An automatic workflow was proposed for recognizing and classifying subsurface structures with GPR using computer vision and machine learning techniques. The workflow comprises three stages: first, full-cover GPR measurements are processed to form the C-scans; second, the abnormal areas are extracted from the full-cover C-scans with coefficient of variation-active contour model (CV-ACM); finally, the extracted segments are recognized and classified from the corresponding B-scans with aggregate channel feature (ACF) to produce a semantic map. The selected computer vision methods were validated by a controlled test in the laboratory, and the entire workflow was evaluated with a real, on-site case study. The results of the controlled and on-site case were both promising. This study establishes the necessity of a full-cover 3D GPR survey, illustrating the feasibility of integrating advanced computer vision techniques to analyze a large amount of 3D GPR survey data, and paves the way for automating subsurface modeling with GPR.

Three-Dimensional Reverse Time Migration of Ground-Penetrating Radar Signals

Three-dimensional (3D) ground-penetrating radar (GPR) systems and 3D seismic imaging techniques have been developing fast and evolving rapidly in the last decade. Ray-based migration methods have been successfully applied to processing 3D GPR signals based on the similarity between electromagnetic and seismic waves. However, reverse time migration (RTM) of 3D GPR signals has not been well studied in the past. In this paper, we present a 3D RTM based on Maxwell’s equations for 3D GPR surveys. Migration recovers the true subsurface structure from a distorted and unfocused time profile and suppresses common electromagnetic clutter and noise. RTM based on Maxwell’s equations can consider conductivity directly and compensate for the attenuation within a high-conductivity zone. Compared with 2D RTM, 3D RTM back-propagates both in-line and cross-line signals simultaneously and can include complex 3D permittivity and conductivity models. We have integrated a parallel finite-difference time-domain (FDTD) algorithm based on a hybrid MPI and OpenMP scheme to reduce the computational cost of 3D problems. The 3D RTM experiments on an anomaly of “EM” shape and a realistic sand dune model demonstrate the effective recovery of 3D subsurface structures.

Locating and characterizing burials using 3D ground-penetrating radar (GPR) and terrestrial laser scanning (TLS) at the historic Mueschke Cemetery, Houston, Texas

We use ground-penetrating radar (GPR), terrestrial laser scanning (TLS), cemetery records, and oral histories to locate and characterize unmarked burials at the Mueschke Cemetery in Houston, Texas. We employ TLS surveying for the reconnaissance of the various types of burials in the cemetery. Concrete burials were used after 1940 while wooden burials were common before 1940. Strong correlations of TLS shallow surface mounds and depressions with known burials provide substantial motivation to follow-up with a detailed investigation with GPR. To calibrate the GPR survey, we undertook excavation and lab experiments to measure relative dielectric constant (3−10), moisture content (14%–21%), and electrical conductivity (0.36–0.38S/m) of the location's soil. Grain-size experiments indicated that the soil at the survey location is reasonable for GPR, which was subsequently reinforced by the high quality GPR images with a depth of penetration of about 2.5m and a vertical resolution of 6cm. Numerical modeling provided a signature of a concrete burial with a flat rectangular top while a wooden burial resembles a typical hyperbola shape. 2D GPR template surveys over concrete and wooden burials show strong response and follow the numerical modeling findings. Three methods were used to estimate the soil radar velocity: CMP, two-way travel time (TWTT), and hyperbola fitting. These techniques yielded a soil velocity of about 0.055–0.067m/ns. 3D GPR surveys produced anomalies that were consistent with headstones, cemetery records, and oral histories. Taking GPR depth slices overlain on the TLS elevation profile and cemetery records strongly indicated the location of two suspected unmarked graves including the oldest burial (1875) of James West.

Imaging of an active fault: Comparison between 3D GPR data and outcrops at the Castrovillari fault, Calabria, Italy

We have integrated and analyzed a 3D ground-penetrating radar (GPR) volume with a trenching exposure data set to evaluate the potential of these methods individually and combined for study of a fault zone. We chose a test site across a branch of the active Castrovillari fault in the Northern Calabria (Southern Italy). This tectonic structure is one of the most active in the area, and it has generated strong earthquakes in the past. Based on analysis of previously collected data, a 3D GPR survey was carried out 1.2 m from a fault outcrop. The goal was to use the GPR volume to guide and optimize the excavation of a trench and then to use the trenching data to validate the GPR volume interpretation. We used seismic interpretation software to display vertical and horizontal sections and for horizon tracking and attribute analyses. We obtained quantitative information on the geometry of structural and geologic features, such as fault strike and dip angle, defining the boundaries of different stratigraphic units. We validated our GPR data interpretation with the outcrop section and trench wall demonstrating the benefits of GPR in extensional tectonics environments and the great potential of the combined geologic and geophysical approach.

4D ground-penetrating radar during a plot scale dye tracer experiment

Flow phenomena in the unsaturated zone are highly variable in time and space. Thus, it is challenging to measure and monitor such processes under field conditions. Here, we present a new setup and interpretation approach for combining a dye tracer experiment with a 4D ground-penetrating radar (GPR) survey. Therefore, we designed a rainfall experiment during which we measured three surface-based 3D GPR surveys using a pair of 500MHz antennas. Such a survey setup requires accurate acquisition and processing techniques to extract time-lapse information supporting the interpretation of selected cross-sections photographed after excavating the site. Our results reveal patterns of traveltime changes in the measured GPR data, which are associated with soil moisture changes. As distinct horizons are present at our site, such changes can be quantified and transferred into changes in total soil moisture content. Our soil moisture estimates are similar to the amount of infiltrated water, which confirms our experimental approach and makes us confident for further developing this strategy, especially, with respect to improving the temporal and spatial resolution.

Ground penetrating radar and microwave tomography 3D applications for the deck evaluation of the Musmeci bridge in Potenza, Italy

An extensive experimental and numerical investigation has been carried out to assess the status of the 'Ponte sul Basento' (1967–1976), in the town of Potenza (Basilicata region, southern Italy), better known as the Musmeci bridge. Architecturally, the bridge is a considerable reinforced 20th century concrete structure that was designed and built by the Italian architect Sergio Musmeci (1926–1981). Moreover, the bridge represents an important element of the infrastructural network, linking the city centre to the Potenza-Sicignano highway, crossing the Basento river and the railway close to the main train station of the city. Recently, due to ageing and continuous and significant traffic, the bridge started to be affected by several problems such as water infiltration. Within the presented study, a widespread ground penetrating radar (GPR) survey has been designed to investigate the geometrical characteristics of the bridge deck (Gerber saddles, internal stiffening walls, pillar supports) and detect the presence of defects or damage due to water infiltration and traffic fatigue. Concerning this, a 900 MHz 3D GPR survey has been performed along a zone of one of the lanes on the road surface. Moreover, a second 1500 MHz 3D survey has been carried out at the bottom of the bridge deck in order to gain detailed information about an important structural element of the bridge, the Gerber saddle. Both results have been processed following two approaches: the first a classical time-domain processing session based on commercial software and the use of migration; the second in microwave tomography, an advanced frequency domain automatic PC-based inversion algorithm. In this paper, we present a comparative interpretation of both kinds of processed results, and provide considerations about the investigated structures.

A 3D ground penetrating radar imaging of the heavy rainfall-induced deformation around a river levee: a case study of Ara River, Saitama, Japan

This paper describes a three-dimensional ground penetrating radar (GPR) survey carried out around a levee of the Ara River in Saitama, Japan, where deformation of the ground was observed after heavy rainfall associated with the typhoon of September 2007. The high-density 3D GPR survey was conducted as a series of closely adjacent four directional sets of 2D surveys at an area surrounding vertical cracks on the paved road caused by deformations induced by heavy rain. The survey directions of the 2D surveys were 0, 90, 45, and –45 degrees with respect to the paved road and the intervals between lines were less than 0.5 m. The 3D subsurface structure was accurately imaged by the result of data processing using Kirchhoff-type 3D migration. As a result, locations and vertical continuities of the heavy rainfall induced cracks in the paved road were clearly imaged. This will be a great help in considering the generation mechanisms of the cracks. Moreover, the current risk of a secondary disaster was found to be low, as no air-filled cavities were detected by the 3D GPR survey.

Sinkhole structure imaging in covered Karst terrain

Ground penetrating radar (GPR) and resistivity techniques have been widely used to map the locations of sinkholes in covered karst terrain. To determine whether a sinkhole is a likely preferential conduit for groundwater flow, however, requires higher‐resolution imaging than that used in conventional sinkhole mapping surveys. Field observations combined with simulated surveys for a 15‐m diameter 3‐m deep sinkhole in west‐central Florida are used to assess the resolution of GPR and resistivity surveys targeting the semiconfining unit that floors the sinkhole depression. 2D resistivity surveys clearly show the central depression as well as resistivity contrasts between the cover sediments within and outside of the sinkhole, but are inadequate for resolving breaches in the semiconfining unit or underlying conduits. A 3D GPR survey resolves vertical structure on the order of tens of centimeters within the semiconfining unit, as well as indicators of conduits that extend several meters beneath the central depression. 3D GPR thus holds promise for imaging hydrologically significant features of sinkholes.

Full-resolution 3D GPR imaging

Noninvasive 3D ground-penetrating radar (GPR) imaging with submeter resolution in all directions delineates the internal architecture and processes of the shallow subsurface. Full-resolution imaging requires unaliased recording of reflections and diffractions coupled with 3D migration processing. The GPR practitioner can easily determine necessary acquisition trace spacing on a frequency-wavenumber (f-k) plot of a representative 2D GPR test profile. Quarter-wavelength spatial sampling is a minimum requirement for full-resolution GPR recording. An intensely fractured limestone quarry serves as a test site for a 100-MHz 3D GPR survey with 0.1 m × 0.2 m trace spacing. This example clearly defines the geometry of fractures in four different orientations, including vertical dips to a depth of 20 m. Decimation to commonly used half-wavelength spatial sampling or only 2D migration processing makes most fractures invisible. The extra data-acquisition effort results in image volumes with submeter resolution, both in the vertical and horizontal directions. Such 3D data sets accurately image fractured rock, sedimentary structures, and archeological remains in previously unseen detail. This makes full-resolution 3D GPR imaging a valuable tool for integrated studies of the shallow subsurface.

Hydraulic effects of shales in fluvial-deltaic deposits: Ground-penetrating radar, outcrop observations, geostatistics, and three-dimensional flow modeling for the Ferron Sandstone, Utah

Ground-penetrating radar (GPR) surveys, outcrop measurements, and cores provide a high-resolution 3D geologic model to investigate the hydraulic effects of shales in marine-influenced lower delta-plain distributary channel deposits within the Cretaceous-age Ferron Sandstone at Corbula Gulch in central Utah, USA. Shale statistics are computed from outcrop observations. Although slight anisotropy was observed in mean length and variogram ranges parallel and perpendicular to pale of low \((\overline L _ \bot \approx 1.2\overline L _ )\), the anisotropy is not statistically significant and the estimated mean length is 5.4 m. Truncated Gaussian simulation was used to create maps of shales that are placed on variably dipping stratigraphic surfaces interpreted from high-resolution 3D GPR surveys, outcrop interpretations, and boreholes. Sandstone permeability is estimated from radar responses calibrated to permeability measurements from core samples. Experimentally designed flow simulations examine the effects of variogram range, shale coverage fraction, and trends in shale coverage on predicted upscaled permeability, breakthrough time, and sweep efficiency. Approximately 1500 flow simulations examine three different geologic models, flow in the 3 coordinate directions, 16 geostatistical parameter combinations, and 10 realizations for each model. ANOVA and response models computed from the flow simulations demonstrate that shales decrease sweep, recovery, and permeability, especially in the vertical direction. The effect on horizontal flow is smaller. Flow predictions for ideal tracer displacements at Corbula Gulch are sensitive to shale-coverage fraction, but are relatively insensitive to twofold variations in variogram range or to vertical trends in shale coverage. Although the hydraulic effects of shale are statistically significant, the changes in flow responses rarely exceed 20%. As a result, it may be reasonable to use simple models when incorporating analogous shales into models of reservoirs or aquifers.