Trace Modelling: A Quantitative Approach to the Interpretation of Ground-Penetrating Radar Profiles

The analysis of ground-penetrating radar data generally relies on the visual identification of structures on selected profiles and their interpretation in terms of buried features. In simple cases, inverse modelling of the acquired data set can facilitate interpretation and reduce subjectivity. These methods suffer from severe restrictions due to antenna resolution limits, which prevent the identification of tiny structures, particularly in forensic, stratigraphic, and engineering applications. Here, we describe a technique to obtain a high-resolution characterization of the underground, based on the forward modelling of individual traces (A-scans) of selected radar profiles. The model traces are built by superposition of Ricker wavelets with different polarities, amplitudes, and arrival times and are used to create reflectivity diagrams that plot reflection amplitudes and polarities versus depth. A thin bed is defined as a layer of higher or lower permittivity relative to the surrounding material, such that the top and bottom reflections are subject to constructive interference, determining the formation of an anomalous peak in the trace (tuning effect). The proposed method allows the detection of ultra-thin layers, well beyond the Rayleigh vertical resolution of GPR antennas. This approach requires a preliminary estimation of the instrumental uncertainty of common monostatic antennas and takes into account the frequency-dependent attenuation, which causes a spectral shift of the dominant frequency acquired by the receiver antenna. Such a quantitative approach to analyzing radar data can be used in several applications, notably in stratigraphic, forensic, paleontological, civil engineering, heritage protection, and soil stratigraphy applications.