Next Article in Journal
Spatio-Temporal Variability of the Aerosol Optical Depth over the Gorky and Cheboksary Reservoirs in 2022–2023
Previous Article in Journal
The Detection of Desert Aerosol Incorporating Coherent Doppler Wind Lidar and Rayleigh–Mie–Raman Lidar
Previous Article in Special Issue
Common-Mode Clutter Filtering for the Problem of Sounding Multilayer Media Using Ground-Penetrating Radar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Antenna Radiation Pattern and Penetration Depth in Ground Penetrating Radar Field Missions

1
Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (IZMIRAN), Troitsk, 108840 Moscow, Russia
2
JSC Company VNIISMI, 127566 Moscow, Russia
3
Faculty of Fundamental Physics, Moscow Pedagogical State University (MPGU), 118881 Moscow, Russia
4
Faculty of Managemrnt, Bauman Moscow State Technical University (MSTU), 105005 Moscow, Russia
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(23), 5452; https://doi.org/10.3390/rs15235452
Submission received: 10 October 2023 / Revised: 8 November 2023 / Accepted: 15 November 2023 / Published: 22 November 2023

Abstract

:
This article discusses the methods and results of assessing the angular resolution and sounding depth of enhanced-power ground penetration radars obtained during archaeological and geographical expeditionary works in various natural areas. Elongated local objects were used as test objects to evaluate the horizontal radiation pattern of the Loza–V georadar in the upper- and lower-half spaces. The depth of operation of the Loza–N low-frequency radar was estimated during a geophysical study of a unique natural object in the Siberian taiga. The variability of the GPR antenna radiation patterns in different materials (air, dry, or wet soils) confirms the necessity of quantitative measurements with controlled electrophysical parameters.

1. Introduction

Ground penetrating radar (GPR, georadar) [1,2,3] is widely used in archeology, geology, civil engineering, and security issues [4,5,6]. The problem of GPR radiation directivity plays an important role in planning the measurement paths in rural and urban conditions. On the one hand, a wide radiation pattern increases the field of view of the device and makes it possible to estimate the permittivity of the subsurface medium from the attenuation of the signal by separating the transmitter and receiver antennas [1] (other proposed approaches require an improvised test platform [7,8] or double passing the same path with different antennas [9]). On the other hand, especially when working in an urban environment, side reflections make it difficult to detect and identify objects of interest. Of special interest are the long resistively loaded antennas of low-frequency georadars used in geophysics and geological exploration. While the high-frequency antennas can be tested in an electromagnetic anechoic chamber, the performance of long low-frequency GPR antennas can be studied only during practical fieldwork.
In this article, we describe the directional and penetration characteristics of powerful GPR of the Loza series (“deep penetration radar”, DPR) [9,10,11]. All of them have antennas made as dipoles of different lengths, with the resistive loading roughly following Wu–King’s law [12] and minimizing “ringing” at the antenna’s central frequency. Although the radiating and receiving properties of GPR antennas have been studied in a number of theoretical and experimental works (see [13,14,15,16,17,18,19,20]), georadar missions provide a good opportunity to assess their performance under realistic operating conditions.
In the first part of this paper, we present the results of theoretical and experimental studies of the GPR antenna radiation patterns in the subsurface medium and in the upper hemisphere (the former characteristic is necessary for the interpretation of georadar B-scans [1], whereas the latter is useful when working in urban conditions, with multiple interfering aerial reflections). These methodological works have been performed by the research team of the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radio Wave Propagation (IZMIRAN) and JSC “Company VNIISMI” during a joint expedition organized by the Institute of Archeology of Crimea of the Russian Academy of Sciences in July, 2022. Among the studied archeological objects, the ancient Jewish cemetery in the vicinity of Kerch city, the remains of an ancient water pipe near the village Geroevka (El Tiygen), and the destroyed Venetian fortress in Tikhaya Balka should be mentioned. The results of the field campaign were reported at a scientific-practical conference [21] in December 2022. Along with archaeological research, the features of the local relief and subsurface artifacts were used in order to evaluate the technical characteristics of the Loza–V GPR [10]. First of all, we were interested in estimating the radiation pattern of the resistively loaded dipole antennas used in this enhanced georadar power. This important parameter has been studied in a number of theoretical and experimental works, e.g., see [14,22,23,24,25]. However, taking into account the construction, size, and potential of our DPR antennas, designed for geophysical and industrial applications, full-scale field measurements were invaluable.
The second part of the presented experimental material was obtained during the expedition organized in September 2022 by the Russian Geographic Society, the newspaper “Komsomolskaya Pravda” and the TV channel “Russia1”. The object of the study was the Patomsky crater, an unusual geological structure in the Siberian taiga discovered in 1949 in the Irkutsk region [26]. This rare object (a stone ring 80 m in diameter with a central cone composed of large limestone fragments) attracted the attention of researchers, and various hypotheses of its origin were put forward. Our GPR cross sections, performed with Loza–N GPR [11], were aimed at choosing between the most probable volcanic model [27] and the early meteorite hypothesis [28]. From a technical point of view, this experiment confirmed the possibility of GPR sounding the first hundreds of meters of a subsurface medium. Geophysics revealed a fundamental difference between deep GPR echoes from a massive stone ground and from a probable water-filled volcanic crater.

2. Materials and Methods

Experimental studies of the directional pattern of dipole antennas were carried out using a GPR of the Loza series [10,11]. The main design features of the device were outlined in [29]. In this “deep penetration radar” (DPR) construction, everything was performed to increase its depth sounding capabilities (high voltage EM pulse and long resistively loaded antennas). The design features made it possible to use DPR in the study of subsurface structures at depths of up to 100–150 m in “heavy” low-resistivity soils and up to 200–300 m in high-resistivity rocks.
The fundamental features of the Loza GPR series design are follows:
-
Ultra-high power: The peak power of the transmitted EM pulse was brought up to a practical limit, determined only by the insulating properties of the environment, by means of a high-voltage discharger supplying a probing pulse of voltage from 5 to 21 kV to the antenna.
-
Concentration of signal energy in the low-frequency part of the spectrum: In order to achieve maximum depths, the maximum energy of the probing pulse is shifted to the lowest frequencies, within the 1–50 MHz frequency band of the receiver determined by the length of the transmitter antenna. The medium-frequency version Loza–V, operating in the frequency band of 50–300 MHz, is equipped with 100, 200, and 300 MHz dipole antennas (Figure 1). The Loza–N DPR is supplied with 50 MHz (3 m long), 25 MHz (6 m), 15 MHz (10 m), and 10 MHz (15 m) transmitter and receiver half-wave dipole antennas (Figure 2). All antennas are designed with a resistive load that gradually grow towards the ends of the dipole [12].
-
Large dynamic range of the reflected signal registration: The peak power of the transmitter is brought to a physical limit and is limited only by the electrical breakdown of the environment (ground and air). The use of ultra-powerful transmitters and low-frequency antennas provides a dynamic range of reflected signals of more than 120 dB, which allows one to work in environments with high conductivity, such as loam or wet clay. The registration system allows one to digitize the signal in the entire dynamic range without changing the GPR settings. The resistive damping makes it possible to obtain probing pulses of a practically non-oscillating nature [14,22,23,24], see Figure 3.
The efficiency of GPR sounding (first of all, its depth) depends on the parameters of the environment. The radiation patterns of the antennas are largely determined by the electrical properties of the ground on which the antennas are placed. When the GPR antenna (linear electric dipole) is placed tightly at the interface between two media (ground surface), a radiation pattern is formed, essentially oriented towards the medium with higher dielectric permittivity ε r = n 2 (or soil refractive index) n = ε r [30,31]. Approximately ε r = n 2 times more radiation is emitted towards the lower hemisphere (into the ground) than upwards (into the air).
Figure 4a,b illustrate the effect of the interface on the harmonic and pulsed EM radiation of an idealized infinite-line antenna into the subsurface medium. In both cases, a noticeable maximum of the radiation pattern in the far zone is observed in the directions corresponding to the angle of total internal reflection Θ = arcsin 1 n , cf. [13,18,32]. These estimates allow one to assess the GPR radiation pattern in the transversal “H plane”. A similar but smoother pattern, without “horns”, demonstrates the popular GSSI bow-tie antenna [20]. In both cases, directivity is rather poor.
The problem of characterizing the directional properties of ground penetration radar has been discussed in a number of classical works, e.g., see [13,14,15,32]. The simplest approach, in the far zone of harmonic radiation, is just the superposition of the Tx and Rx radiation patterns. The situation becomes more complicated for non-harmonic radiation, even in the model case of an infinitely long GPR line antenna [30] excited by a single current step [31]. This can be illustrated using the following set of color plots (Figure 5). The left frame (a) shows the spatial distribution of the pulsed radiation of an infinite linear vibrator placed on the ground at a certain moment ( c t = 1 m). A similar spatio-temporal radiation pattern has a receiver antenna (b). The radar spatial resolution is determined by the product (c) showing a bright spot at the intersection of (a) and (b) patterns—cf. [15].
These analytical results are presented here in order to provide a rough qualitative picture of realistic GPR radiation characteristics and illustrate the problem of their quantitative estimates. First of all, the Engheta–Papas–Elachi model of an infinitely long line antenna [30] works only at distances small or similar, compared with the antenna length. Moreover, numerical calculations show that in time domain, the peak amplitude values tend to a far-field pattern with sharp edges at the critical total reflection angles |θ| = η quite gradually, at distances of hundreds of meters, unrealistic for most applications [31], cf. Figure 4b. Finally, Figure 5 highlights the complexity and ambiguity of the term “pulsed GPR radiation pattern”, cf. [13,14,15,16,17].
The experimental evaluation of the radiation pattern of subsurface radar requires either the creation of an expensive measuring stand or the ingenious use of local objects as natural test objects [29]. The latter technique was used by participants of the Crimean expedition in 2022 in the spare intervals between archaeological works. In our measurements, natural test objects were at hand (power line wires, as shown in Figure 6, and non-metallic pipes of an ancient aqueduct were discovered during archaeological excavation, as shown in Figure 7).
Before proceeding to the presentation of the experimental results, the following remarks should be made. We did not seek to express the characteristics of GPR emitting and receiving systems in the exact formulas of the antenna theory [33]. Parameters such as gain and directivity cannot be accurately measured in the field, and also do not have a strict definition for monopulse GPR antennas due to the variability of the pulse waveform. In addition, the resulting radiation pattern depends not only on the design and size of the Tx and Rx antennas but also on the properties of the underlying surface [13,18,19,20], which may change considerably along the measurement path. Nevertheless, a practical assessment of the directivity characteristics is possible and very useful, whereas the calculation of the antenna gain for an enhanced power GPR, such as Loza–V, is of no practical importance.
In order to confirm the need to take into account the radiation pattern of the GPR antennas in the upper hemisphere, we provide an example of working with the Loza–V GPR in a rural area. The task was to localize an underground cable along a high-voltage power line. The primary B-scan across the route is shown in Figure 8a. First of all, the two systems of hyperbolas strike the eye in the left and right parts of the figure, which could be interpreted as an image of localized objects at a depth of 10 and 3 m, if it were not for the excessively wide angular opening of the hyperbola wings, not leading to their focusing by standard wave migration methods [34,35], as shown in Figure 8b. The soil dielectric constant ε ~ 4 was assessed when using the standard CMP method [1] (other practical methods for estimating the electrical parameters of the soil are discussed, for example, in Refs. [7,9]). Visual control of the situation prompted the operator to perform the wave migration procedure with a value ε = 1 corresponding to the dielectric constant of free space, as shown in Figure 8c. As a result, we obtained an impressive picture of five wires of a high-voltage power line at a height of 19–25 m, compared with Figure 8d, and a single overhead cable at a height of about 6 m in the right corner of the picture. An image of the desired underground cable is visible as a narrow hyperbola in the upper-right corner of Figure 8a and as a faint spot in Figure 8b.
The use of natural linear objects to evaluate the Loza–V GPR directivity pattern in the upper and lower hemispheres is described below in Section 3.1 and Section 3.2. Section 3.3 and Section 3.4 are devoted to the phenomenon of the Patom crater, where we discuss the results of the GPR cross-section of this natural object and apply the solution of the 1D inverse problem to explain the peculiarities of the deep subsurface sounding.

3. Experimental Results

3.1. GPR Directivity Pattern in Upper Hemisphere

In the spring of 2020, scientists from the Institute of Archeology of the Crimea and the Institute of Archeology of the Russian Academy of Sciences identified one of the earliest Jewish necropolises in the vicinity of the modern city of Kerch. The site of the necropolis belonged to the community of Panticapaeum, mentioned by ancient authors and located at Cape Ak-Burun. It was here that the archeologists collected more than three dozen Jewish tombstones and fragments containing inscriptions and images (Figure 8a). The revealed names of the deceased members of the Jewish community add to our knowledge of the onomastics of the “Bosporan” civilization.
In the framework of the 2022 expedition, the territory of the necropolis was studied via Loza–V scanning, with 100 MHz antennas and a 5 kV transmitter. A dense accumulation of anomalies, interpreted as graves, was found in a relatively small area. Five previously unknown objects on the radar image were identified as crypts, and up to 40 tombstones with Jewish symbols were found on the site (Figure 9). Further study of the necropolis during the 2023 summer field campaign allowed the archeologists to trace the path of Jewish migration to the Bosporus in the 1st century BC and the first centuries of our era.
The first purely radiophysical task in this expedition was to evaluate the GPR radiation patterns in the upper hemisphere. This was performed according to the following method: the Loza–V GPR antenna set was installed in the center of a circle divided into 30-degree sectors (Figure 6). In each of the positions, four measurements were taken, after which the antennas were rotated by 30 degrees, with the radar handle remaining above the center of the circle. The elongated scattering test object was a power line wire at a height of 15 m (one of the pillars is visible in Figure 6). The operator moved under the power line in a circle synchronously with the rotation of the antennas.
As a result of the GPR signal measurements, the horizontal radiation pattern of the 200 MHz dipole Loza–V antenna in the upper hemisphere (in air) under real field conditions (on wet ground with ε r 8 ) was obtained. From visual observation, it was understood that in the radargram, the subsurface reflections could be mixed with the echo of the power line above the GPR set at a height of about 15 m (round-trip delay of about 100 ns) (see Figure 6). In order to obtain a realistic radar image, clearly seen on the radargram and B-scan (Figure 10a), it was converted using “radar velocity” v r a d a r = 0.5 c = 15 10 7 m/s corresponding to ε r = 1 . The resulting radiation pattern, calculated as the angular dependence of the relative maximum of the radar pulse scattered by the power line cable, is shown in Figure 10b. The difference between the maximum and minimum of the directivity pattern is more than 40 dB (it is interesting to note that the radiation pattern has some asymmetry; the closer the receiver antenna is to the object, the greater the amplitude of the reflected signal). The obtained estimate of the radiation pattern is of practical importance when planning a GPR survey in urban conditions with unremovable objects of “aerial” reflection.

3.2. GPR Directivity Pattern in Subsurface Medium

The evaluation of the Loza–V GPR directivity pattern in the subsurface medium was performed during the Crimean summer expedition of 2022 at two archeological sites. One of them was found in the area of the Hebrew necropolis mentioned in Section 3.1. It was a crypt examined at the request of the Institute of Archeology of Crimea. The crypt was opened earlier by robbers and had no archeological value. However, it was useful for us as a reference for identifying the radar image of a known object and estimating the GPR radiation pattern. The radargram in Figure 11a shows both the measured radar delay time t in nanoseconds (right scale) and the object depth d in meters, converted using the estimated “radar velocity” V r = c 2 ε 5 10 7 m/s (left scale). The measured depth at the top of the crypt (3.5 m) was taken into account in the interpretation of the radargram.
The measured radiation pattern, as shown in Figure 11b, reflects the good focusing of the GPR signal and the high level of the reflected pulse, which is due to the high dielectric permittivity and low loss rate of the propagation medium.
Another test object for GPR directivity measurements was found in the course of archeological work near the village of Geroevka (El Tiygen). In its vicinity, about 2.5 km from the village, there is an ancient well (Figure 12). According to archaeologists, this well and other objects unearthed as a result of landslides or during the construction of a gas pipeline (Figure 13a), are related to the ancient settlement of Nymphaeum.
The objective of our GPR survey was to obtain, in a non-disruptive manner, information confirming the existence of an ancient underground water supply in the area between the identified objects. It was carried out, following the instructions of archaeologists, using Loza–V GPR, 200 MHz antennas, and a 5 kV transmitter (attenuated by 30 dB).
A structure similar to a cross-section of two parallel non-metallic pipes was identified in all the marked areas, as shown in Figure 13b and A1 in Figure 14a,b. The radar image can be interpreted as the two ceramic pipes of an ancient aqueduct. Apart from the pure archeological interest, they served as a convenient elongated test object for evaluating the GPR antenna directivity under typical field conditions. Figure 14c illustrates the practical method of directivity pattern assessment.
By moving horizontally at a distance x from the top of the hyperbola (image of a compact object at a depth z = h ), we increased the distance according to the law h 2 + x 2 = h / cos α . The amplitude of the reflected signal at the observation point x is taken from the waveform in Figure 14b at the corresponding time delay t = t 0 / cos α , where t 0 = 2 h c ε r (see Figure 14c). In this way, the effective radiation pattern F ( h , α ) of the transmit-receiving GPR antenna was built. Although it depends not only on the soil parameters but also on the depth and the scattering pattern of the elongated buried test object, these estimates are qualitatively true and very useful when working in typical field conditions. The results of the GPR radiation pattern measurements via radar pulse scattering by the buried pipe are shown in Figure 15.
Another assessment of GPR antenna directivity was carried out during a survey of the Venetian settlement territory (Tikhaya Balka, XIV-XV centuries AD) at the order of the Archaeological Expedition of the State Hermitage Museum. As a test object, we used the foundation of the destroyed fortress wall at a depth of 1.5 m. This experiment gave a larger width of the main lobe—about 45 degrees at the level of 0.7 (see Figure 16). This difference might be due to the lower refractive index of dry sandy soil.
In both cases, one can notice a rather smooth directivity pattern, which seemingly contradicts the classical Engetta–Papas–Elachi results [30] predicting “horns” of the GPR antenna radiation pattern at the angle of total internal reflection, as shown in Figure 4a. However, as thorough analysis, laboratory experiment, and numerical calculation confirm, realistic resistively loaded GPR antennas have smooth radiation patterns [19,20,32].

3.3. Deep Radar Probing of Patomsky Crater

In what follows, we discuss the experimental material obtained during the research mission organized in September 2022 by the Russian Geographic Society, the newspaper “Komsomolskaya Pravda” and the TV channel “Russia1”. The object of the study was the Patomsky crater, an unusual geological structure in the Siberian taiga discovered in the Bodaibo district of the Irkutsk region [26]. This rare object (a stone ring 80 m in diameter, with a central cone composed of large fragments of limestone, as shown in Figure 17) attracted the attention of geo-scientists. The organizers invited experienced geophysicists V. L Sulyandikov and F. P. Morozov for the scientific support of the project.
The origin of the mysterious Patomsky crater, located in the north of the Irkutsk region, attracted the attention of many researchers. V.V. Kolpakov, who discovered the Patomsky crater in 1949 during geological survey work, put forward a hypothesis about its formation as a result of a meteorite fall [26]. The famous volcanologist S.V. Obruchev argued with this hypothesis, leaning toward the volcanic origin of the crater. This point of view was seriously supported by geological and dendrological studies [27]. Nevertheless, some experts defend the meteorite hypothesis, assuming the presence of a foreign body in the Patomsky crater at a depth of 180–200 m [28].
Our GPR cross-sections (Figure 18) performed by F. P. Morozov confirm the most probable volcanic model [27]. The survey was carried out from the Earth’s surface using a JSC VNIISMI low-frequency Loza–N DPR [9,26] equipped with a 21 kV pulse transmitter and antennas with a central frequency of 25 MHz (antenna length of 6 m). Unlike high-frequency 200–400 MHz GPR [36,37] designed for probing the first subsurface meters, the time scale for this measurement was set to 4096 ns (maximum depth of about 200–250 m). Similarly to Figure 10, both radar delay time (right scale) and object depth (left scale) are shown in Figure 19b and below, with an estimated “radar velocity” of 5.5 cm/ns.
GPR B–scan No 6 recorded along the outer ring of the crater (Figure 19a) reflects the geological structure outside the crater cone. According to the interpretation of experienced geologists, it is composed of the following elements. Horizon 0–60 m is represented by limestone with a homogeneous undisturbed structure. Deeper than 60 m, undisturbed rocks are recorded, markedly differing in physical properties from the overlying horizon. The main result of the GPR data analysis along this profile is that there are no vertical anomalous structures outside the crater.
The next B–scan, No 3, was recorded along the latitudinal diameter of the crater (Figure 20a). Here, the geological structure is represented by the following features: horizon 0–100 m is filled with destroyed limestone; undisturbed rocks lie deeper. The most interesting anomalous structure of vertical development is recorded on the central hill of the crater. Two more similar smaller-scale anomalies are noted in the ring moat on the northern and southern sides of the crater—see Figure 20b.
A series of horizontal cross-sections of the obtained GPR data (Figure 21a) made it possible to construct an impressive 3D electrodynamic model of the Patomsky crater (Figure 21b). The depths of the selected cross-sections were estimated using a hypothetical radar velocity of v r = 5.5 cm/ns.
Of special interest are the deepest sections—from 180 to 280 m. All of them clearly show the contour of the central vertical anomalous structure, located under the central hill of the Patomsky crater. At depths of less than 200 m, anomalous structures of a smaller scale appear. It is noteworthy to mention that the well-distinguishable ring structures around the main anomaly are registered in all sections. Such ring structures (concentric zones of compression and expansion) are typical for the areas of explosive impact on the rock.
Additional information about the nature of the Patomsky crater is provided by the analysis of the reflected signal waveforms—see Figure 22. The signals (A) and (C), measured above the rock mass, are characterized by the amplitude deviation towards negative values (−53 and −51 dB, respectively). Such properties of the reflected signal, according to the results of theoretical analysis, indicate low values of the propagation medium permittivity [11]. In practice, this means that the probing signal propagates in a high-resistance medium (in this case, monolithic limestone).
Waveform (B) shows a significant deviation in the signal amplitude in the region of positive values.
The deviations in the reflected signal amplitude indicate that the EM wave propagates in a medium with high values of permittivity or conductivity [38,39,40]. Such changes in the properties in a limited vertical portion of the limestone mass can only be associated with a local influx of moisture, the appearance of which can substantially increase the dielectric constant. A local fracture, which can arise due to tectonic processes or from a local impact, can serve as a collector for moisture ingress into the rock mass.
The characteristic shape of the Patomsky crater surface, which consists of a ring structure with a bulk cone, an annular moat, and a central hill, can also, with some probability, indicate that the crater originally arose due to the impact of a meteorite, which would support the original hypothesis of A. M. Portnov [28]. The existence of such a local channel of limited depth, which ensures the transport of surface waters in the cold Siberian climate, would lead to the second mechanism—ice diapirus (frost heaving). The processes of frost heaving every winter straighten and restore the shape of the Patomsky crater with a bulk cone, annular ditch, and central hill. The latter version also seems quite probable. The only indisputable fact is the existence of a water-saturated channel in a homogeneous limestone mass under the central hill of the Patomsky crater, which was convincingly confirmed by the georadar.

3.4. Qualitative Solution of Inverse Problem

Many analytical and numerical approaches have been developed to solve the problem of GPR pulse propagation in a smoothly layered subsurface medium—e.g., see [41]. For a qualitative analysis of the aforementioned GPR survey results, we used our approximate solution of the one-dimensional inverse problem of electromagnetic sounding [42,43]—a time-domain version of the classical coupled wave method [44,45]. It provides a general picture of the reflected pulse formation on the subsurface medium gradients, consistent with numerous experimental results, and yields a closed-form solution of the simplified inverse problem. We represent the electromagnetic signal received by the GPR antenna as E ( 0 , s ) = f ( s ) + g ( s ) . Here, s = c t is the normalized propagation time of the probing pulse ( c is the speed of light), f ( s ) is the signal waveform, and the integral
g ( s ) = 1 4 0 ε ( z ) ε ( z ) f s 2 0 z ε ( z ) d ζ d z
describes the cumulative effect of partial reflections of the incident wave from dielectric permittivity gradients in the subsurface medium. This approximation, obtained using the coupled wave method [31], makes it possible to find an explicit solution ε ( z ) to the inverse problem. Its parametric representation has the form
ε z ( s ) = ε 0 exp 4 0 s Q ( s ) d s z ( s ) = 1 2 ε 0 0 s exp 2 0 s Q ( r ) d r d s ,
where the function Q ( s ) = ε [ z ( s ) ] 8 ε 3 / 2 [ z ( s ) ] is determined by the inverse transform of the ratio
Q ( s ) = 1 2 π + i δ + i δ Q ˜ ( k ) e i k s d k , Q ˜ ( k ) = g ˜ ( k ) f ˜ ( k )
of the incident and reflected pulse Fourier transforms: f ˜ ( k ) = 0 f ( s ) e i k s d s and g ˜ ( k ) = 0 g ( s ) e i k s d s . Absorption can be taken into account by introducing the complex permittivity, ε ˜ = ε + 4 π i σ ω , which leads to a significant complication of the analysis [31]. Under the condition 4 π σ < < ω , the perturbation method can be used. The calculations are simplified and lead to a slight modification of the main formula, Formula (1):
g ( s ) = 1 4 0 ε ( z ) ε ( z ) exp 4 π c 0 z σ ( ζ ) ε ( ζ ) d ζ f s 2 0 z ε ( ζ ) d ζ d z
Unfortunately, one measurement of the reflected signal waveform is not enough to determine the two unknown functions: ε ( z ) and σ ( z ) . To determine them, one can use an a priori electrodynamic model of the subsurface medium (for example, by assuming a constant loss angle: q = arctan 4 π σ ( z ) c ε ( z ) = C o n s t ) or perform two measurements of the reflected signal with the antennas having different frequency responses [42]. In the first case, we obtain an integral equation
g ( s ) = 1 8 0 s ε [ z ( r ) ] ε 3 / 2 [ z ( r ) ] exp q 2 r f s r d r = 0 s P ( r ) f s r d r ,
with P ( s ) = 1 2 π + i δ + i δ g ˜ ( k ) f ˜ ( k ) e i k s d k = ε [ z ( s ) ] 8 ε 3 / 2 [ z ( s ) ] exp q 0 z ε ( ζ ) d ζ and some generalization of the parametric solution (2)–(3):
ε z ( s ) = ε 0 0 s exp 4 0 s P ( r ) exp ( q 2 r ) d r z ( s ) = 1 2 ε 0 0 s exp 2 0 s P ( r ) exp ( q 2 r ) d r d s ,
The second approach requires some complications of the experimental technique.
In this work, to obtain a qualitative estimate, we neglect ohmic absorption, which is apparently not very significant under experimental conditions, and consider the level of the reflected signal as a measure of the vertical gradients of the subsurface medium permittivity. For an asymptotic estimate, we transform Equation (5) via integration by parts, defining f 0 = f = 0 , f s = d d s h ( s ) , where the function h ( s ) has a single maximum at s = l . Let, for example, h ( s ) = 1 cos ( 2 π e a s ) , where l = log 2 a is the characteristic pulse length. For a short pulse, integral (5) is determined using a narrow neighborhood of the function h ( s ) maximum. By removing slowly varying functions from under the integral sign (5), we obtain the following:
g ( s ) = 1 8 0 s d d s ε [ z ( r ) ] ε 3 / 2 [ z ( r ) ] exp q 2 r h s r d r 1 16 ε 1 / 2 ε ε 3 / 2 e q 2 r 0 s h s r d r .
The integral H ( s ) = 0 s h s r d r = s + 1 a c i ( 2 π e a s ) c i ( 2 π ) for small values of s can be approximated using a linear function H ( s ) 2 s 0.86 l , which rapidly tends to the constant H ( s ) 3.5 l depending on the antenna half-length l (see Figure 23).
Following the method developed in our works [42,43], we consider Formula (7) as a differential equation for determining the ground permittivity profile ε ( z ) . Indeed, for a given model of the GPR pulse f ( s ) and the measured waveform of the reflected signal waveform g ( s ) in the differential equation ε 1 / 2 ( z ) ε ( z ) ε 3 / 2 ( z ) q ε ( z ) ε 3 / 2 ( z ) = 4 Φ ( s ) , the right-hand side Φ ( s ) = 4 g ( s ) H ( s ) is a known function that can be integrated by substitution ε ( z ) = ε 0 exp V ( s ) , s = 2 0 z ε ( r ) d r : ε ( z ) ε 3 / 2 ( z ) = 2 V ˙ ( s ) , ε ( z ) ε 3 / 2 ( z ) = 4 ε 1 / 2 V ¨ ( s ) ; V ¨ ( s ) q 2 V ˙ ( s ) = Φ ( s ) . Taking into account the initial condition ε ( 0 ) = ε 0 its explicit solution is given by the integral
V ( s ) = 2 q 0 exp ( q 2 r ) Φ ( r ) Φ ( r + s ) d r 0 s Φ ( r ) d r
Equation (9) is simplified in the absence of losses ( q = 0 ):
ε ( z ) = ε 0 exp V ( s ) , V ( s ) 0 Φ ( r + s ) Φ ( r ) r d r = 0 Φ ( r ) min ( r , s ) d r
In combination with the integral z = Z ( s ) = 1 2 ε 0 0 s exp 1 2 V ( r ) d r , Equation (9) yields an explicit solution to the subsurface sounding problem. Three issues remain as follows:
(1)
There is no registration of the emitted pulse amplitude in the current GPR models (some approaches to the problem are outlined in [46]);
(2)
Since drilling near this rare natural object is excluded, the numerical value of the permittivity at the depth of far reflections can be estimated only from model calculations (we tried several likely values for different moisture saturation).
Moreover, our 1D model does not describe the divergence of the probing signal in the subsurface medium, which leads to some errors in the return pulse amplitude. If necessary, this effect can be taken into account by introducing an appropriate divergence factor.
With these considerations, an approximate solution of the 1D inverse problem was used to process the GPR sounding data of the Patomsky crater, as performed by F. P. Morozov, as a part of the Komsomolskaya Pravda expedition. Three A-scans were selected in Figure 24a: in the center of the stone “dome” (A 38) and at two points on the annular shaft—A 21 and A 48. For convenience, the color amplitude scale used throughout the paper is repeated in Figure 24b. Attention was paid to the difference in the sign of deep reflections when sounding the central dome and its periphery—see Figure 22. To construct the dielectric permittivity profile, taking into account the above considerations, our analytical solution was written in a form convenient for practical assessments, with explicitly indicated hypothetical permittivity values at the reflection depths of the initial and tail parts of the GPR pulse ε 0 = ε ( 0 ) and ε 1 = ε ( z 1 ) :
ε = ε [ Z ( s ) ] = ε 0 ε 1 ε 0 R ( s ) , R ( s ) = 0 g ( r ) H ( r ) min ( r , s ) d r / 0 g ( r ) H ( r ) r d r .
Figure 25 shows plots of the soil dielectric permittivity reconstructed using Equation (10) from the initial value consistent with geological data to the deep ground values, corresponding to the hypothesis of meltwater saturation of the crater, for three hypothetical values of ε 1 = ε ( z 1 ) = 8 , 12 , 16 .
The maximum dielectric permittivity at the depths of about 100 m may indicate increased moisture saturation of the soil filling the well.
These estimates are rather speculative and preliminary. A detailed study of the Patomsky crater structure requires additional georadar and magnetometric sounding. The only conclusion that can be drawn from the DPR measurements is the absence of a compact foreign body in the crater, which would result in a noticeable backscattering.

4. Conclusions

The results of our Crimean campaign have cultural, historical, and methodological value. The usefulness of georadar surveys in archeological research can hardly be overestimated; the number of new artefacts found during a short GPR inspection of an ancient Jewish cemetery in the vicinity of Kerch City is comparable with the total volume of previous findings. Archeological work near the village of Geroevka not only helped to outline the territory of a rare ancient water pipeline but may also guide the direction of further construction works. The objects found in the territory of the ancient Venetian settlement of the XIV-XV centuries in Tikhaya Balka will also contribute to the material history of the Middle Ages.
From a technical point of view, the Crimean expedition provided a good chance to estimate the GPR antenna radiation patterns in different materials (air, dry, or wet soils) and emphasized the necessity to organize quantitative measurements in typical media with controlled electrophysical parameters.
Our Siberian mission must be evaluated from different points of view. First of all, it was the first attempt at the objective electrophysical inspection of this rare natural object. The technical parameters of the enhanced-power GPR allowed us to obtain a radar image of the Patomsky crater up at a depth of about 200 m. The results of the deep georadar sounding did not support the hypothesis of a massive foreign body present in the crater, which makes the less exotic volcanic origin of the Patomsky phenomenon more probable. Other non-destructive techniques, such as multi-frequency radar probing combined with magnetic measurements, can provide new information on the origin and structure of the crater. In this regard, as well as to popularize the tourist route and organize the protection of this unique natural site, the “Komsomolskaya Pravda” expedition can play a good role.

Author Contributions

Conceptualization, P.M. and A.P.; methodology, P.M.; software, A.P. and M.L.; investigation, F.M., M.L. and L.B.; formal analysis, A.P.; resources, P.M.; data curation, F.M., M.L. and L.B.; writing and funding acquisition, P.M. and A.P.; English editing, M.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation in the framework of the research project “Deep Penetration Radar: Theory, Methods, Experiment”, grant № 22-1200083. Scientific management and technical support have been provided by IZMIRAN and JSC Company VNIISMI.

Data Availability Statement

The data presented in this study are available upon request from IZMIRAN and JSC VNIISMI.

Acknowledgments

The first part of the study was conducted with legal support from the Institute of Archeology of Crimea of the Russian Academy of Sciences, the State Hermitage Museum, and the Research Center for the History and Archeology of Crimea. The Siberian expedition was organized and supported by the editors of the newspaper “Komsomolskaya Pravda”, and the scientific equipment was provided by VNIISMI and JSC Lenzoloto.

Conflicts of Interest

The Authors Pavel Morozov and Fedor Morozov were partly employed by the JSC Company VNIISMI, with no conflict of interest. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Daniels, D.J. Ground Penetrating Radar, 2nd ed.; IEE Radar, Sonar and Navigation Series 15; The Institution of Electrical Engineers: London, UK, 2004. [Google Scholar]
  2. Benedetto, A.; Pajewski, L. (Eds.) Civil Engineering Applications of Ground Penetrating Radar; Book Series: Transactions in Civil and Environmental Engineering; Springer: New Delhi, India, 2015; p. 385. [Google Scholar]
  3. Persico, R. Introduction to Ground Penetrating Radar: Inverse Scattering and Data Processing; Wiley-IEEE Press: New York, NY, USA, 2014; p. 392. [Google Scholar]
  4. Hou, F.; Rui, X.; Fan, X.; Zhang, H. Review of GPR Activities in Civil Infrastructures: Data Analysis and Applications. Remote Sens. 2022, 14, 5972. [Google Scholar] [CrossRef]
  5. Srimuk, P.; Boonpoonga, A.; Kaemarungsi, K.; Athikulwongse, K.; Dentri, S. Implementation of and experimentation with ground-penetrating radar for real-time automatic detection of buried improvised explosive devices. Sensors 2022, 22, 8710. [Google Scholar] [CrossRef] [PubMed]
  6. Barone, P.M. Understanding Buried Anomalies: A Practical Guide to GPR; Lambert Academic Publishing: Saarbrücken, Germany, 2016; p. 92. [Google Scholar]
  7. Álvarez López, Y.; Fernández, M.G.; Arboleya Arboleya, A.; González-Valdés, B.; Rodríguez Vaqueiro, Y.; Las-Heras Andrés, F.; Pino García, A. SAR-based technique for soil permittivity estimation. Int. J. Remote Sens. 2017, 38, 5168–5185. [Google Scholar] [CrossRef]
  8. Lauro, S.E.; Mattei, E.; Barone, P.M.; Pettinelli, E.; Vannaroni, E.G.; Valerio, G.; Comite, D.; Galli, D.A. Estimation of subsurface dielectric target depth for GPR planetary exploration: Laboratory measurements and modeling. J. Appl. Geophys. 2013, 93, 93–100. [Google Scholar] [CrossRef]
  9. Popov, A.; Morozov, P.; Edemsky, D.; Edemsky, F.; Pavlovski, B.; Zapunidi, S. Expedient GPR survey schemes. In Proceedings of the 11th International Radar Symposium, Vilnius, Lithuania, 16–18 June 2010. [Google Scholar]
  10. Kopeikin, V.V.; Krasheninnikov, I.V.; Morozov, P.A.; Popov, A.V.; Fang, G.; Liu, X.; Zhou, B. Experimental verification of Loza–V GPR penetration depth and signal quality. In Proceedings of the 4th International Workshop on Advanced GPR, Naples, Italy, 27–29 June 2007; pp. 230–233. [Google Scholar]
  11. Berkut, A.I.; Edemsky, D.E.; Kopeikin, V.V.; Morozov, P.A.; Prokopovich, I.V.; Popov, A.V. Deep penetration subsurface radar: Hardware, results, interpretation. In Proceedings of the 9th International Workshop on Advanced Ground Penetrating Radar (IWAGPR), Edinburgh, UK, 28–30 June 2017; pp. 1–6. [Google Scholar]
  12. Wu, T.T.; King, R.W.P. The cylindrical antenna with nonreflecting resistive loading. IEEE Trans. Antennas Propag. 1965, 13, 369–373. [Google Scholar] [CrossRef]
  13. Smith, G. Directive properties of antennas for transmission into a material half-space. IEEE Trans. Antennas Propag. 1984, 32, 232–246. [Google Scholar] [CrossRef]
  14. Arcone, S. Numerical studies of the radiation patterns of resistively loaded dipoles. J. Appl. Geophys. 1995, 33, 39–52. [Google Scholar] [CrossRef]
  15. Millard, S.; Shaari, A.; Bungey, J. Field pattern characteristics of GPR antennas. NDT E Int. 2002, 35, 473–482. [Google Scholar] [CrossRef]
  16. Poljak, D.; Dori, V.; Birki, M.; El Khamlichi Drissi, K.; Lallechere, S.; Pajewski, L. A simple analysis of dipole antenna radiation above a multilayered medium. In Proceedings of the 9th Interational Workshop IWAGPR, Edinburgh, UK, 28–30 June 2017; pp. 1–6. [Google Scholar]
  17. Perez-Gracia, V.; González-Drigo, R.; Di Capua, D. Horizontal resolution in a non-destructive shallow GPR survey: An experimental evaluation. NDT E Int. 2008, 41, 611–620. [Google Scholar] [CrossRef]
  18. Warren, C.; Giannopoulos, A. Creating FDTD models of commercial GPR antennas using Taguchi’s optimisation method. Geophysics 2011, 76, G37–G47. [Google Scholar] [CrossRef]
  19. Diamanti, N.; Annan, A.P. Characterizing the energy distribution around GPR antennas. J. Appl. Geophys. 2013, 99, 83–90. [Google Scholar] [CrossRef]
  20. Warren, C.; Giannopoulos, A. Characterisation of a ground penetrating radar antenna in lossless homogeneous and lossy heterogeneous environments. Signal Process. 2017, 132, 221–226. [Google Scholar] [CrossRef]
  21. Gukin, V.D.; Morozov, P.A.; Andriushchenko, V.V. Prospecning of region “Calm bay” Nº2 on the 2022 south-eastern Crimean archaeological expedition of the State Hermitage Museum. In Proceedings of the VII Scientific-Practical Conference “Crimea: Science, Culture, Politics; MBUK FMD: Feodosia, Ukraine, 2023; pp. 149–163. ISBN 987-5-6045030-7-2. [Google Scholar]
  22. Wagner, S.; Pham, A.-V. The ultrawideband elliptical resistively loaded vee dipole. IEEE Trans. Antennas Propag. 2020, 68, 2523–2530. [Google Scholar] [CrossRef]
  23. Esselle, K.P.; Stuchly, S.S. Pulse-receiving characteristics of resistively loaded dipole antennas. IEEE Trans. Antennas Propag. 1990, 38, 1677–1683. [Google Scholar] [CrossRef]
  24. Samaddar, S.N.; Mokole, E.L. Transient behavior of radiated and received fields associated with a resistively loaded dipole. In Proceedings of the Ultra-Wideband Short-Pulse Electromagnetics 4, Tel-Aviv, Israel, 14–19 June 1999; pp. 165–180. [Google Scholar]
  25. Montoya, T.P.; Smith, G.S. Land mine detection using a ground-penetrating radar based on resistively loaded Vee dipoles. IEEE Trans. Antennas Propag. 1999, 47, 1795–1806. [Google Scholar] [CrossRef]
  26. Kolpakov, V.V. Misterious crater in Patom Highlands. Priroda 1951, 2, 58–61. (In Russian) [Google Scholar]
  27. Antipin, V.S.; Fedorov, A.M. The origin of Patom Crater, East Siberia, from geological and geochemical data. Dokl. Earth Sci. 2008, 423, 1335–1339. [Google Scholar] [CrossRef]
  28. Portnov, A.M. Crater on the Patom Highlands. Priroda 1962, 11, 102–103. (In Russian) [Google Scholar]
  29. Berkut, A.I.; Popov, A.V.; Kopeikin, V.V.; Morozov, P.A. Theprinciples of increasing the penetration depth and application experience of “Loza” GPR. In International School «Registration of Subsurface Targets by Short-Range Radar», Bauman MSTU; Bauman MSTU: Moscow, Russia, 2016. [Google Scholar]
  30. Engheta, N.; Papas, C.H.; Elachi, C. Interface extinction and subsurface peaking of the radiation pattern of a line source. Appl. Phys. B 1981, 26, 231–238. [Google Scholar] [CrossRef]
  31. Edemsky, F.; Popov, A.; Zapunidi, S. A time domain model of GPR antenna radiation pattern. Int. J. Electron. Telecommun. 2011, 57, 407–411. [Google Scholar] [CrossRef]
  32. Lampe, B.; Holliger, K. Resistively loaded antennas for ground-penetrating radar: A modeling approach. Geophysics 2005, 70, K23–K32. [Google Scholar] [CrossRef]
  33. Balanis, A. Antenna Theory: Analysis and Design, 4th ed.; John Wiley & Sons: Hoboken, NJ, USA, 2016; 1095p. [Google Scholar]
  34. Johansson, E.M.; Mast, J.E. Three-dimensional ground penetrating radar imaging using synthetic aperture time-domain focusing. In Advanced Microwave and Millimeter-Wave Detectors; LLNL: Livermore, CA, USA, 1994. [Google Scholar] [CrossRef]
  35. Martinez-Lorenzo, J.A.; Rappaport, C.M.; Quivira, F. Physical Limitations on detecting tunnels using underground-focusing spotlight synthetic aperture radar. IEEE Trans. Geosci. Remote Sens. 2011, 49, 65–70. [Google Scholar] [CrossRef]
  36. Poluha, B.; Porsani, J.L.; Almeida, E.R.; dos Santos, V.R.N.; Allen, S.J. Depth estimates of buried utility systems using the GPR method: Studies at the IAG/USP Geophysics test site. Int. J. Geosci. 2017, 8, 726–742. [Google Scholar] [CrossRef]
  37. Fediuk, A.; Wunderlich, T.; Wilken, D.; Rabbel, W. Ground penetrating radar measurements in shallow water environments—A case study. Remote Sens. 2022, 14, 3659. [Google Scholar] [CrossRef]
  38. Zeng, Q.; Delisle, G.Y. Transient analysis of electromagnetic wave reflection from a stratified medium. In Proceedings of the 2010 Asia-Pacific International Symposium on Electromagnetic Compatibility, Beijing, China, 12–16 April 2010; pp. 881–884. [Google Scholar]
  39. Kaganovsky, Y.; Heyman, E. Pulsed beam propagation in plane stratified media: Asymptotically exact solutions. In Proceedings of the 2010 URSI International Symposium on Electromagnetic Theory, Berlin, Germany, 16–19 August 2010; pp. 829–832. [Google Scholar]
  40. Epstein, P.S. Reflection of waves in an inhomogeneous absorbing medium. Proc. Natl. Acad. Sci. USA 1930, 16, 627–637. [Google Scholar] [CrossRef] [PubMed]
  41. Lambot, S.; Slob, E.; Vereecken, H. Fast evaluation of zero-offset Green’s function for layered media with application to ground-penetrating radar. Geophys. Res. Lett. 2007, 34, L21405. [Google Scholar] [CrossRef]
  42. Vinogradov, V.A.; Kopeikin, V.V.; Popov, A.V. An approximate solution of 1D inverse problem. In Proceedings of the 10th International Conference on Ground Penetrating Radar, Delft, The Netherlands, 21–24 June 2004; pp. 95–98. [Google Scholar]
  43. Prokopovich, I.; Popov, A.; Pajewski, L.; Marciniak, M. Application of Coupled-Wave Wentzel-Kramers-Brillouin Approximation to Ground Penetrating Radar. Remote Sens. 2018, 10, 22. [Google Scholar] [CrossRef]
  44. Brekhovskikh, L.M. Waves in Layered Media, 2nd ed.; Academic Press: New York, NY, USA, 1980; 520p, ISBN 9780323161626. [Google Scholar]
  45. Bremmer, H. Propagation of Electromagnetic Waves. In Handbuch der Physik/Encyclopedia of Physics; Flugge, S., Ed.; Springer: Berlin/Goettingen/Heidelberg, Germany, 1958; Volume 4/16, pp. 423–639. [Google Scholar]
  46. Popov, A.; Prokopovich, I.; Edemskii, F. GPR image deconvolution from antenna current waveform. In Proceedings of the Days on Diffraction (DD), St. Petersburg, Russia, 4–8 June 2018; pp. 231–236. [Google Scholar]
Figure 1. Loza–V GPR with 200 MHz (0.75 m long) antennas.
Figure 1. Loza–V GPR with 200 MHz (0.75 m long) antennas.
Remotesensing 15 05452 g001
Figure 2. Loza–N GPR with 50 MHz (3 m long) antennas.
Figure 2. Loza–N GPR with 50 MHz (3 m long) antennas.
Remotesensing 15 05452 g002
Figure 3. (a) Measured pulse in a resistively loaded dipole antenna (negative of the original screen [14], a.u.), and (b) GPR pulse model [11].
Figure 3. (a) Measured pulse in a resistively loaded dipole antenna (negative of the original screen [14], a.u.), and (b) GPR pulse model [11].
Remotesensing 15 05452 g003
Figure 4. (a) Radiation pattern of an infinite line vibrator on the ground–air boundary [30]. Harmonic radiation and different values of the ground permittivity ε r = 2 (∙∙∙); 4 (−−−); 9 (–––). (b) Peak radiation pattern of an impulse line antenna [31] for ε r = 4 and a one-period sinusoidal current pulse form. Different distances between the transmitter dipole and point receiver: (1) 200 m, (2) 12 m, (3) 4 m, and (4) 2 m. (c) Total energy radiation patterns of a 1.5 GHz GSSI bow-tie antenna over a lossless half-space of permittivity ε r = 5, at distances from 0.1 to 0.57 m (H-plane, colors mark amplitude levels [20]).
Figure 4. (a) Radiation pattern of an infinite line vibrator on the ground–air boundary [30]. Harmonic radiation and different values of the ground permittivity ε r = 2 (∙∙∙); 4 (−−−); 9 (–––). (b) Peak radiation pattern of an impulse line antenna [31] for ε r = 4 and a one-period sinusoidal current pulse form. Different distances between the transmitter dipole and point receiver: (1) 200 m, (2) 12 m, (3) 4 m, and (4) 2 m. (c) Total energy radiation patterns of a 1.5 GHz GSSI bow-tie antenna over a lossless half-space of permittivity ε r = 5, at distances from 0.1 to 0.57 m (H-plane, colors mark amplitude levels [20]).
Remotesensing 15 05452 g004
Figure 5. The spatio−temporal radiation pattern of an infinite line GPR, with 1 m separation between the antennas. (a) Transmitter pulse propagating into the ground ( n = 2 ) for the distance c t = 1 m; (b) receiver directivity pattern; (c) Resulting GPR pattern.
Figure 5. The spatio−temporal radiation pattern of an infinite line GPR, with 1 m separation between the antennas. (a) Transmitter pulse propagating into the ground ( n = 2 ) for the distance c t = 1 m; (b) receiver directivity pattern; (c) Resulting GPR pattern.
Remotesensing 15 05452 g005
Figure 6. Measurements of GPR radiation pattern in the upper hemisphere under an electric power line.
Figure 6. Measurements of GPR radiation pattern in the upper hemisphere under an electric power line.
Remotesensing 15 05452 g006
Figure 7. Ancient well and pipe outlets.
Figure 7. Ancient well and pipe outlets.
Remotesensing 15 05452 g007
Figure 8. GPR search of a buried cable in a rural area. (a) Primary B–scan; (b) wave focusing with soil dielectric constant ε ~ 4 ; (c) wave focusing with vacuum permittivity value ε = 1 ; (d) overhead electric power line.
Figure 8. GPR search of a buried cable in a rural area. (a) Primary B–scan; (b) wave focusing with soil dielectric constant ε ~ 4 ; (c) wave focusing with vacuum permittivity value ε = 1 ; (d) overhead electric power line.
Remotesensing 15 05452 g008
Figure 9. Fragments of Jewish tombstones.
Figure 9. Fragments of Jewish tombstones.
Remotesensing 15 05452 g009
Figure 10. (a) GPR B–scan of the electric power line. Conversion of the radar delay time in the apparent “object depth” is performed with ε r = 1 ( v r a d a r = 0.5 c ) ; (b) Horizontal directivity pattern. Amplitude color scale is shown in the bottom of (a) plot (relative units).
Figure 10. (a) GPR B–scan of the electric power line. Conversion of the radar delay time in the apparent “object depth” is performed with ε r = 1 ( v r a d a r = 0.5 c ) ; (b) Horizontal directivity pattern. Amplitude color scale is shown in the bottom of (a) plot (relative units).
Remotesensing 15 05452 g010
Figure 11. (a) Subsurface radar B–scan of the ancient Jewish cemetery (Kerch area). The image of the ancient crypt top is marked with a white arrow. The same color scale as in Figure 10; (b) Measured directivity pattern of Loza–V GPR calculated for wet soil with ε r 8 (orange color intensity reflects field strength).
Figure 11. (a) Subsurface radar B–scan of the ancient Jewish cemetery (Kerch area). The image of the ancient crypt top is marked with a white arrow. The same color scale as in Figure 10; (b) Measured directivity pattern of Loza–V GPR calculated for wet soil with ε r 8 (orange color intensity reflects field strength).
Remotesensing 15 05452 g011
Figure 12. Ancient well.
Figure 12. Ancient well.
Remotesensing 15 05452 g012
Figure 13. The photo of antique water supply ceramic pipes (a). Three B–scans revealing cross-sections of two ceramic pipes, with enlarged square inserts (b). The same color scale as in Figure 10 and Figure 11.
Figure 13. The photo of antique water supply ceramic pipes (a). Three B–scans revealing cross-sections of two ceramic pipes, with enlarged square inserts (b). The same color scale as in Figure 10 and Figure 11.
Remotesensing 15 05452 g013
Figure 14. GPR B–scan with marked cross-sections A1 and A2 (a); two A–scans in the selected cross-sections, with logarithmic color scale. Scan A1 reveals the subsurface pipe at about 2 m depth (b); practical method of GPR directivity pattern assessment (c). The amplitude color scale is shown in the bottom of (b).
Figure 14. GPR B–scan with marked cross-sections A1 and A2 (a); two A–scans in the selected cross-sections, with logarithmic color scale. Scan A1 reveals the subsurface pipe at about 2 m depth (b); practical method of GPR directivity pattern assessment (c). The amplitude color scale is shown in the bottom of (b).
Remotesensing 15 05452 g014
Figure 15. Loza–V GPR directivity estimated on the subsurface object shown in Figure 13a (yellow color intensity reflects field strength).
Figure 15. Loza–V GPR directivity estimated on the subsurface object shown in Figure 13a (yellow color intensity reflects field strength).
Remotesensing 15 05452 g015
Figure 16. Loza–V GPR directivity on dry sand ground (orange color intensity reflects field strength).
Figure 16. Loza–V GPR directivity on dry sand ground (orange color intensity reflects field strength).
Remotesensing 15 05452 g016
Figure 17. Patomsky crater—general view.
Figure 17. Patomsky crater—general view.
Remotesensing 15 05452 g017
Figure 18. GPR survey paths over the Patomsky crater (numbered red lines, passed along the arrows). Isohypses with heights are marked in meters (yellow lines).
Figure 18. GPR survey paths over the Patomsky crater (numbered red lines, passed along the arrows). Isohypses with heights are marked in meters (yellow lines).
Remotesensing 15 05452 g018
Figure 19. (a) GPR probing of the crater outer ring; (b) Radar B–scan No 6. Here and below a conditional amplitude scale is used.
Figure 19. (a) GPR probing of the crater outer ring; (b) Radar B–scan No 6. Here and below a conditional amplitude scale is used.
Remotesensing 15 05452 g019
Figure 20. GPR probing across the latitudinal diameter of the crater. (a) Aerial photo; (b) GPR B–scan No 3 with marked A–scans (A, B, C ) at 12, 40 and 72 meters of the profile.
Figure 20. GPR probing across the latitudinal diameter of the crater. (a) Aerial photo; (b) GPR B–scan No 3 with marked A–scans (A, B, C ) at 12, 40 and 72 meters of the profile.
Remotesensing 15 05452 g020
Figure 21. Loza–N GPR qualitative image of the Patomsky crater. (a) Horizontal sections; (b) 3D electrodynamic model.
Figure 21. Loza–N GPR qualitative image of the Patomsky crater. (a) Horizontal sections; (b) 3D electrodynamic model.
Remotesensing 15 05452 g021
Figure 22. Three GPR signal waveforms from B–scan shown in Figure 20b: x = 12 , 40 , 72 m (AC) and the amplitude color scale.
Figure 22. Three GPR signal waveforms from B–scan shown in Figure 20b: x = 12 , 40 , 72 m (AC) and the amplitude color scale.
Remotesensing 15 05452 g022
Figure 23. Model source function H ( s ) for l = 3 , 6 , 10 m. Explanation of colors is given in the box; the dashed lines mark the common tangent.
Figure 23. Model source function H ( s ) for l = 3 , 6 , 10 m. Explanation of colors is given in the box; the dashed lines mark the common tangent.
Remotesensing 15 05452 g023
Figure 24. (a) Radargram of deep subsurface probing of the Patomsky crater. B-scan with marked profiles: A 21 m, A 38 m, and A 48 m; (b) Logarithmic amplitude color scale is shown in the bottom of the linear plot.
Figure 24. (a) Radargram of deep subsurface probing of the Patomsky crater. B-scan with marked profiles: A 21 m, A 38 m, and A 48 m; (b) Logarithmic amplitude color scale is shown in the bottom of the linear plot.
Remotesensing 15 05452 g024
Figure 25. Depth permittivity profile in the center (A 38 m) and on the periphery of the crater (A 21 and A 48 m) reconstructed from the GPR data.
Figure 25. Depth permittivity profile in the center (A 38 m) and on the periphery of the crater (A 21 and A 48 m) reconstructed from the GPR data.
Remotesensing 15 05452 g025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Morozov, P.; Morozov, F.; Lazarev, M.; Bogolyubov, L.; Popov, A. Characterization of Antenna Radiation Pattern and Penetration Depth in Ground Penetrating Radar Field Missions. Remote Sens. 2023, 15, 5452. https://doi.org/10.3390/rs15235452

AMA Style

Morozov P, Morozov F, Lazarev M, Bogolyubov L, Popov A. Characterization of Antenna Radiation Pattern and Penetration Depth in Ground Penetrating Radar Field Missions. Remote Sensing. 2023; 15(23):5452. https://doi.org/10.3390/rs15235452

Chicago/Turabian Style

Morozov, Pavel, Fedor Morozov, Maxim Lazarev, Leonid Bogolyubov, and Alexei Popov. 2023. "Characterization of Antenna Radiation Pattern and Penetration Depth in Ground Penetrating Radar Field Missions" Remote Sensing 15, no. 23: 5452. https://doi.org/10.3390/rs15235452

APA Style

Morozov, P., Morozov, F., Lazarev, M., Bogolyubov, L., & Popov, A. (2023). Characterization of Antenna Radiation Pattern and Penetration Depth in Ground Penetrating Radar Field Missions. Remote Sensing, 15(23), 5452. https://doi.org/10.3390/rs15235452

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop