1. Introduction
Precision agriculture (PA), supported by various technologies, is a rapidly emerging field for managing the agricultural landscape on a large scale by considering the field variability for increasing agricultural productivity while minimizing negative environmental impacts and the production cost
[1][2][3]. Agricultural practices such as irrigation, land preparation, and fertilization are typically applied uniformly across the entire field, treating the heterogeneous fields as homogeneous. However, the spatial and temporal heterogeneity of the agricultural landscape can be effectively monitored through implementing different technologies to support PA, including (1) geographic information systems (GIS), (2) satellite-based global positioning systems (GPS), (3) remote sensing, (4) drones, (5) the Internet of Things, (6) artificial intelligence, and (7) different proximal sensors, for example, geophysical techniques
[1][4][5]. These digital technologies are applied to collect, process, monitor, and map the spatiotemporal variability of the agricultural landscape. Their purpose is to improve agronomic performance, enhance crop productivity and develop decision support tools
[1][3][6][7][8][9][10][11][12][13]. By utilizing the latest technologies, such as unmanned aerial vehicles, agricultural machinery, and robotic technologies, the required amounts of water, nutrients, and agrochemicals for plant growth and development can be accurately applied to specific locations within the agricultural landscape and within the appropriate timeframe
[1][4][6][8][9][14][15][16]. PA offers several benefits, including improved soil fertility and health, increased water productivity and food security, minimized soil and water pollution, and reduced labor force requirements, as well as the overuse of resources such as water, fertilizer, seeds, and energy, leading to lower production costs
[4][6][14][17][18]. The incorporation of spatiotemporal variability mapping with geo-statistics and GIS marked a major advance in PA to a new level to identify management zones
[7][8][19].
The spatiotemporal variability of crop factors and subsurface physical, chemical, and hydrological properties, and processes are crucial in PA
[1]. Hydro-geophysics is an efficient approach that includes multi-scale probing and high-resolution imaging techniques for accurately obtaining the spatiotemporal variability of the subsurface hydrological processes and soil properties
[20][21][22][23]. Traditional methods such as soil sampling and laboratory analysis are destructive, labor-intensive, costly in large scales, time-consuming, and mainly provide point scale measurements only
[1][20][24][25]. Commonly used geophysical techniques in hydro-geophysics are electrical resistivity tomography, electromagnetic induction (EMI), self-potential, ground-penetrating radar (GPR), induced polarization, surface nuclear magnetic resonance, gravity, magnetics, and seismic methods
[23][24][25][26].
2. Basic Operating Principles of Ground-Penetrating Radar
GPR is a near-surface electromagnetic proximal sensor commonly used in agriculture and environmental applications
[27][28][29]. Commercially available GPR systems use unguided electromagnetic waves with frequencies ranging from 10 MHz to 3600 MHz
[28][30]. Resolution and depth of penetration are vital factors of GPR applications. Radar wave propagation velocity and wave attenuation depend on primary electromagnetic properties such as the relative dielectric permittivity (
εr), electric conductivity σ, and magnetic permeability (µ) of the soil/media
[28][29].
εr is mainly controlled by the water content of the subsurface as the permittivity of liquid water overwhelms those of other soil constituents
[28][31]. Equation (1)
shows the relationship between velocity (
v) and
εr, where
c is the propagation velocity of electromagnetic waves in free space, which is equal to the speed of light (0.3 m/ns)
[28]. The penetration depth of GPR is mainly determined by the σ of subsurface materials and the operation frequency
[27][28][32][33]. When EM waves travel through conductive materials, EM energy is lost as heat through the electrical current. In addition, energy loss in EM waves is due to increasing operating frequency and scattering. This energy loss is called attenuation and, therefore, reduces the penetration depth of GPR waves
[27][28][34][35]. GPR wave attenuation also results from geometrical spreading in 3D wave propagation. The resolution of GPR increases with increasing frequency, increasing bandwidth, and decreasing wave velocity, since the wavelength (λ) is inversely proportional to the frequency (f) and directly proportional to the velocity (λ = v/f)
[35][36][37]. The radiation pattern of a GPR antenna is the primary determinant of its footprint, which is the subsurface area illuminated by the radiation emitted from the antenna. The shape and extent of the radiation pattern depend on several factors, including the antenna design, frequency, size, shape, and dielectric properties of the subsurface. The beam width of the radiation pattern is also affected by the dielectric permittivity of the subsurface, with higher permittivity resulting in a narrower beam.
The most common GPR system consists of an impulse generator that repeatedly sends a source of a particular voltage and frequency to a transmitting antenna. It is worth noting that frequency-domain radars are also available in addition to impulse radars and are also referred to as time-domain radars
[38]. They operate by successively transmitting continuous waves over a specific frequency range. Ground-coupled and air-coupled horn antennas are used for on-ground and air-borne GPR surveys, respectively (
Figure 1).
Figure 1. Ground-penetrating radar (GPR) surveys with (a) ground-coupled antennas (GPR instrument by S. Pathirana), (b) air-coupled antennas (GPR prototype by S. Lambot).
A commonly used GPR system consists of two antennas: a transmitter (T
x) and a receiver (R
x) or multiple receivers. The T
X generates electromagnetic waves, following the frequency-dependent antenna radiation pattern, which travels through the air, air-surface interface, and subsurface. As shown in
Figure 2, the R
x receives airwave, direct ground wave (DGW), reflected wave, and refracted waves
[27][28][39]. EM wave reflection or refraction at a boundary depends on different electrical properties (mainly ε
r in most GPR applications) of the layer above and below the boundary
[40].
Figure 2. Ray paths of ground-penetrating radar (GPR) wave propagation in a two-layer soil which has different dielectric permittivity values (modified from Huisman et al.
[32]).
GPR applications have three main data acquisition modes: (1) reflection profiling, (2) velocity-sounding, and (3) transillumination
[27][28][35]. In the reflection profiling method, the T
x and R
x antennas are kept in a fixed antenna separation, orientation, and station interval, and both antennas are moved along the survey direction
[27]. This method is called the common offset method or fixed offset method (FOM) and produces a vertical 2D image of the subsurface reflections (
Figure 3a)
[35][41]. The GPR velocity-sounding method can be performed using the common mid-point (CMP) method or the wide-angle reflection and refraction (WARR) method. In the CMP method, T
x and R
x are moved apart from each other by keeping the midpoint between the two antennas fixed (
Figure 3b). However, in the WARR method, one antenna (T
x, for example) is kept at a fixed location, and the other antenna (R
x) is moved away by increasing the antenna offset (
Figure 3c). The objectives of the velocity-sounding methods are to estimate the velocity by measuring the two-way travel times for different antenna offsets and then finding the slope of the relationship between the squared travel time and antenna offset
[27][35]. The transillumination survey method places the Tx and Rx on opposite sides of the medium. Under transillumination, zero-offset profiling (ZOP) is a quick method to find anomalous zones by moving T
x and R
x from one station to another at a predetermined step size (
Figure 3d). In the multi-offset gathering (MOG) mode, one antenna (T
x) is kept stationary while the other antenna (R
x) is moved to multiple locations to produce tomographic imaging (
Figure 3e). In vertical reflection profiling (VRP), T
x is placed on the surface, and R
x is placed in the borehole since this method has many advantages over ZOP and MOG (
Figure 3f). The data acquisition method and frequency selected in each survey are based on the application requirements and field conditions
[27][28][42][43].
Figure 3. Data acquisition methods in ground penetration radar (GPR) applications, (
a) fixed offset, (
b) common mid-point, (
c) wide-angle reflection and refraction, (
d) zero offset profiling, (
e) multiple offsets gathering, and (
f) vertical reflection profiling methods (modified from Liu et al.
[44]).
3. Applications of Ground-Penetrating Radar in Soil Studies
This section discusses the applications of GPR to estimate soil properties; states such as SWC, porosity, compaction, salinity, texture, SOM, and other applications such as depth to the groundwater table and capillarity.
3.1. Soil Water Content
SWC estimation is the most extensively used and well-developed GPR application in soil studies. The large-scale estimation and mapping of the spatiotemporal variability of SWC are critical across the agricultural landscape, but these are difficult tasks to complete with traditional methods. For example, mapping high resolution SWC variability using GPR provides the information necessary to optimize the amount of water required for and the timing of crop irrigation. According to previous studies, volumetric soil water content (SWC
V) can be estimated from GPR using various methods such as the reflected wave velocity method, DGW velocity method, transmitted wave velocity method (in boreholes), surface reflection coefficient method, average envelope amplitude (AEA) method, and full-waveform inversion (FWI) method. Subsurface SWC influences the propagation velocity of GPR waves since ε
r varies with SWC and velocity varies with ε
r [35][45]. The well-established and widely used Topp’s equation
[31] estimates SWC
V from ε
r. Capacitance probes, neutron probes, the gravimetric method, and time domain reflectometry (TDR) are used for validating GPR-estimated SWC
[41][43][46][47][48][49][50][51][52].
The reflected wave velocity method is used to acquire SWC
V in deeper soil layers. However, reflectors such as isolated objects or interfaces are required to obtain the travel time of the reflected radar wave. Natural reflectors such as rocks and roots and artificial reflectors such as pipes are point reflectors, while lithologic layers and the water table are interfaces
[33][47][48][53][54][55]. The two-way travel time of the GPR wave above the reflector is measured and converted into radar wave velocity, then to ε
r, and is then used to calculate the SWC
V [52][56][57]. The velocity of the reflected wave is obtained using an FOM or multiple offset methods (MOM): CMP or WARR. When the depth of the reflector is known, the reflected wave velocity is estimated in the FOM by dividing the depth of the reflector by the two-way travel time. Lunt et al.
[56] applied the reflected wave velocity method under natural conditions, while Stoffregen et al.
[58], Loeffler and Bano
[59], Ercoli et al.
[57], and Zhou et al.
[60] applied it under controlled conditions to obtain SWC
V. With advanced analysis techniques, the average velocity of a reflector can be obtained by fitting a hyperbola
[52][61][62][63]. On the other hand, CMP and WARR methods are applied as multiple offset methods when the depth of the reflector is unknown. However, data collection and processing in this method are time-consuming and labor-intensive
[30][44]. Steelman and Endress
[64] applied the multi-frequency CMP method in three different sites to estimate the vertical variation of SWC
V with the reflected wave method during a complete annual cycle, including wetting/drying and freezing/thawing. The reflection method for determining SWC
V has some limitations, as it depends on the continuous availability of reflectors, especially in shallow soils
[52].
- b.
-
Direct ground wave velocity method
On the other hand, the DGW velocity method of GPR can be used to estimate the SWC
V of shallow (uppermost centimeters) soils without employing a reflector
[41][44][47][50][65]. Both the FOM and MOM can be used to estimate SWC
V from the direct ground wave, though the MOM is time-consuming and laborious compared to the FOM
[41][65]. To address this, Sperl
[37] established a method to estimate SWC
V over a large area by combining the MOM and FOM. In this method, the MOM only decides the most suitable antenna separation to distinguish between the DGW and the direct airwave. Then, the FOM is applied with the decided antenna separation to collect data effectively over a large area
[44][66]. Huisman et al.
[65] assessed the accuracy of DGW to estimate and map SWC
V over a large area with WARR and the FOM (single trace analysis). The accuracy of SWC estimated based on WARR measurements was ±0.030 m
3/m
3, while it was ±0.037 m
3/m
3 based on the FOM when compared with the TDR method. Therefore, the authors suggested that the available TDR calibrations, such as Topp’s equation
[31], can also be applied to GPR
[65]. Furthermore, the authors suggested that the most appropriate assessment between electromagnetic methods (TDR and GPR) is ε
r, rather than SWC
V. Galagedara et al.
[41] and Huisman and Bouten
[49] discussed the importance of accurate time-zero calibration to estimate SWC
V using the DGW method. Huisman and Bouten
[49] conducted a sensitivity analysis to study time-zero error at zero antennae offset with WARR measurements. The authors found that GPR simulations could not account for the meantime shift at zero offset due to the time-picking error and SWC
V heterogeneity in the sensitivity analysis. Galagedara et al.
[41] suggested an accurate and stable time-zero picking methodology to estimate SWC
V with an error rate of less than 1% under variable water contents.
Huisman et al.
[67] found that GPR is an efficient technique to map SWC
V over an agricultural landscape, and GPR-estimated water content matched well with TDR results. The accuracy and spatial and temporal variability of the SWC under different water contents, such as irrigation and drainage conditions, have been studied by several researchers
[41][48][50][53][67][68][69][70]. In addition, SWC variation under dry and rainy seasons was studied by Thitimakorn et al.
[71], while Cao et al.
[72] studied three-dimensional soil water dynamics before and after heavy rainfall. These studies compared SWC estimated from the DGW and found good agreements with TDR or gravimetric methods as standard methods. Weihermüller et al.
[73] mapped the spatial variability of SWC with the DGW method at a silty loam site with 450 MHz center frequency antennas. However, Weihermüller et al.
[73] compared GPR (450 MHz)-estimated SWC
V with TDR and volumetric samples and found that the results of the GPR did not agree well with that of the TDR and volumetric samples. The reason behind this was reported to be the high signal attenuation of the GPR due to the relatively high clay content present
[73].
However, previous studies have also identified some issues with the DGW method. It is challenging to distinguish the direct airwave and DGW under dry conditions since the velocity of GPR is high. Therefore, the DGW also reaches the R
x rapidly and interferes with the direct airwave
[41][53][70]. The effective penetration depth of the GPR DGW varies with the soil texture, antenna frequency, and wet/dry conditions present
[41][46][50][66][74]. The heterogeneity of the upper soil layer can produce reflections which interfere with DGW. Radar wave attenuates with higher conductivity, and therefore, the penetration depth decreases with increases in high-conductive materials such as clay. Galagedara et al.
[75] applied a numerical model to investigate the effective sampling depth of the GPR DGW in terms of antenna frequency and dry/wet two-layer conditions. The authors found that the penetration depth of the DGW decreases with increasing frequency and wetness since the sampling depth is a function of the wavelength
[37][66][75].
Table 2 summarizes the previous studies related to the penetration depth of the GPR DGW when estimating SWC.
Table 2. Summary of the effective depth of the ground-penetrating radar direct ground wave studies.
- c.
-
Transmitted wave velocity method
In the transmitted wave velocity method, both the T
x and R
x antennas are placed in boreholes or in surface drains, and the direct wave passing through the media is used to estimate the SWC
V [42][43][77][78][79][80]. In the early stage, ZOP and MOG measurements were widely applied. However, using ZOP and MOG, direct waves interfere with reflected and critically refracted waves in low-velocity zones and underestimate the SWC. Therefore, VRP was introduced
[81][82][83]. VRP requires only one borehole; thus, the ground disturbance and cost are relatively low. Across several studies, the transmitted wave velocity method was applied to estimate SWC
V [42][77][79][80], even though the method is not widely used to estimate SWC
V at the root zone in agriculture
[42][78][80]. Nevertheless, Wijewardana and Galagedara
[43] applied the transmitted direct wave method to estimate SWC
V in raised bed agricultural fields, where T
x and R
x were moved along the surface drains of raised beds.
- d.
-
Surface reflection coefficient method
The surface reflection coefficient method is an off-ground GPR technique; both antennas are moved above ground, and the SWC
V is estimated based on the amplitude of the reflected wave at the soil surface. The underlying equations and modeling hypotheses are detailed and discussed in particular detail in Lambot et al.
[84]. This method determines the reflection coefficient (R) using the amplitude of the reflections from the soil surface (A). The amplitude of the reflections from the perfect electric conductor (A
PEC) is positioned at the same distance as the soil (Equations (2) and (3)). Adekani.
[70], al Hagrey and Müller
[85], Redman et al.
[86], and Redman et al.
[87] applied this method to estimate SWC. This method is more suitable for agricultural applications such as optimizing seeding depth and irrigation management in very thin upper soil layers [0–0.10 m]
[70]. However, the practical applicability of the GPR surface reflection coefficient method for estimating soil moisture is constrained by its sensitivity to calibration height, which is difficult to maintain consistently in real-world field applications. Moreover, the method relies on a simplified assumption of 1D propagation, which does not account for the antenna, further limiting its accuracy.
- e.
-
Average envelope amplitude method
When estimating the SWC using the DGW with the FOM, differentiating the DGW and the direct air wave is often challenging. To avoid this problem, Pettinelli et al.
[88] proposed a method to analyze early time signal amplitude (first arrival direct wave—a combined direct airwave and DGW) without considering the separate travel times of the direct air and ground waves. In the AEA method, SWC is assessed by correlating ε
r variation with attributes of the early time signal of GPR
[88][89]. Furthermore, the AEA is sensitive to both ε
r and σ
a and changes in waveform attributes such as shape, amplitude, and duration with changing ε
r and σ
a [88][90]. Pettinelli et al.
[88] applied the AEA of the early time signal to estimate the SWC in a controlled field condition, Ferrara et al.
[91] applied it under natural field conditions, and Ferrara et al.
[92] applied it under laboratory conditions. Algeo et al.
[93] compared two early time signal analysis methods—the AEA method and the carrier frequency amplitude method—to map SWC and found that both methods have strengths and weaknesses. Another study assessed the applicability of AEA in early time signals to study SWC during irrigation and showed the possibility of this method in estimating SWC in clay-rich soils
[90]. However, this AEA method is still in development, and further research studies are needed to estimate soil properties under different field conditions.
- f.
-
Full waveform inversion method
FWI is a numerical modeling method that retrieves the unknown ε
r and/or σ distribution from a known EM field by fitting a full-wave EM model to the radar measurements
[38]. The EM model tries to describe the radar signal, including the radar source, antenna(s), antenna(s)—medium interactions, and the medium, as accurately as possible. Ernst et al.
[94], Klotzsche et al.
[95], Meles et al.
[96], Klotzsche et al.
[97], Gueting et al.
[98], and Yu et al.
[99] applied FWI for borehole GPR experiments. Lambot et al.
[38][100][101], Jonard et al.
[102][103], Minet et al.
[69], and de Mahieu et al.
[104] applied FWI for air-coupled GPR configurations to estimate SWC. The method was recently used with a drone-borne GPR (
Figure 4) for high-resolution soil moisture mapping in agricultural fields
[105]. The radar equation used in these studies was also generalized to near-field or on-ground GPR conditions
[101][106][107], thereby opening new avenues for agricultural applications. FWI has proven to be a powerful tool for extracting the maximum information from GPR data and facilitating automated data processing. However, the application of FWI has been limited by the inherent complexity of the underlying electromagnetic model and the associated data processing requirements. In that respect, within the framework of the EU’s agROBOfood project, called MIRAGE (grant agreement Nº 825395, 2021–2023), a specifically dedicated radar and software for soil moisture mapping, namely, gprSense, has been developed (
https://www.gprsense.com, accessed on 26 December 2022, Sensar Consulting, Belgium)
[108]. gprSense implements the full-wave radar equation introduced by Lambot et al.
[109] in a user-friendly software platform with an intuitive interface. This tool enables real-time automated FMI and the streaming of soil moisture data, making it accessible to both basic users and advanced scientists.
Figure 4. High-resolution soil moisture map obtained with drone-borne GPR and full-wave inversion (FWI) in an agricultural field in Belgium (GPR prototype by S. Lambot).
Traditionally, the reflection coefficient at the soil surface is assumed to depend solely on the contrast in dielectric permittivity (
εr) between the soil and the air. However, for frequencies below 100 MHz, the soil conductivity (
σ) can also significantly influence the reflection coefficient
[84]. At these low frequencies, the sensitivity of the reflection to
σ becomes higher than that of
εr. Wu and Lambot
[105] proposed a new method for mapping soil
σ using relatively low-frequency drone-borne GPR and full-wave inversion techniques. Their method works best when the dielectric permittivity and conductivity are not too large. Although the sensitivity to
εr is not negligible in this frequency range, their method demonstrated good agreement with the
σ values obtained from EMI surveys.
3.2. Soil Porosity and Soil Compaction
Soil porosity indirectly influences GPR propagation wave velocity and amplitude. When pore spaces are filled with water, it changes the ε
r, which changes the v. Therefore, ε
r at saturation can be used as an indicator of total porosity
[110][111][112][113][114]. Macro-pores can trap more water during saturation (irrigation and excess rainfall) since their infiltration rate is high, but water will drain quickly due to gravity and become dry
[115]. This phenomenon will affect the velocity of GPR waves, where velocity decreases during infiltration and increases during drainage. This velocity variation is rapid in coarse-textured soils, and wetting and drying are faster than in fine-textured soils. Micropores can hold more water against gravity through capillary action, and micropores are more common in clay soil, even though the GPR waves attenuate due to the high σ of clay
[112][115].
Different researchers initially assessed soil porosity using GPR wave velocity
[110][111][113][116]. However, a reliable method to estimate porosity, bulk density, soil compaction, or soil penetration resistance using GPR has yet to be developed and tested. In unsaturated soils, pore spaces are filled with both water and air. The relationship between soil porosity, soil water saturation, and SWC
v is crucial when finding the porosity with GPR. If the soil water saturation is known, porosity can be estimated since SWC
v can be obtained from GPR
[114]. Additionally, under fully saturated conditions, SWC
V is equal to the soil porosity
[115][116]; thus, GPR can be used to estimate soil porosity by measuring the SWC
V at saturation.
Due to the difficulty of estimating porosity alone with GPR, researchers have integrated other geophysical techniques and different theoretical and empirical approaches with GPR. Laboratory-scale experiments were conducted to estimate soil porosity using GPR
[111][113][117] and field-scale experiments
[110][114]. However, the differences between the experimental scales and laboratory-scale measurements must be tested before their application to field conditions
[113][117].
Turesson
[110] estimated the porosity and soil water saturation of sand, and Khalil et al.
[114] estimated the porosity and soil water saturation of a sandstone aquifer using GPR and resistivity techniques using both Topp’s equation (Equation (1) and Archie’s equation (Equation (2)). Ghose and Slob
[113] assessed the porosity and soil water saturation of shallow subsoil by developing integrated GPR-seismic techniques through numerical modeling. Meanwhile, Lai et al.
[111] proposed a new method to determine porosity using GPR, namely, the cyclic moisture variation technique. In this method, the authors determined the variation of ε
r and soil water saturation from a partially saturated state to a fully saturated state of the soil and modified the CRIM equation (Equation (3)) to obtain the porosity.
Researchers have also applied GPR to assess the effect of heavy machinery on soil compaction, a serious problem in agricultural fields
[118][119][120]. With compaction, porosity decreases, consequently increasing the bulk density and penetration resistance, thereby reducing water infiltration
[118][121]. Previous studies showed that compaction changes GPR attributes, i.e., propagation wave velocity and amplitude
[118][119][120][121][122][123]. Moreover, researchers found a negative correlation between GPR wave amplitude and compaction/penetration resistance
[118][119][120][121]. This negative correlation may be due to reduced free water and increased bound water in the soil structure, resulting in increased water-soluble salts in soil pore water, thereby increasing the soil’s bulk σ and attenuating GPR waves and decaying the wave amplitudes
[121]. Another reason could be that increasing the soil density increases the EM wave reflection more than the transmission; hence, the EM energy, wave amplitude, and penetration depth decrease
[121]. On the other hand, Akinsunmade et al.
[118] and Akinsunmade
[119] found that GPR signals penetrate deeper depths in compacted areas than in uncompacted areas. With soil compaction, porosity, and SWC, ε
r decreases, while GPR wave velocity and wavelength increase, increasing the penetration depth
[118][119][120].
3.3. Soil Salinity
Soil σ is the best parameter to use to estimate soil salinity
[27][124][125]. In agricultural fields, soil σ is temporally unstable since it frequently changes, mainly with SWC, due to, for example, irrigation, drainage, leaching, evapotranspiration, fertilizer application, and other soil amendments
[125][126]. Soil σ is mainly measured using: (1) saturated paste extract electrical conductivity (σ
e); (2) apparent electrical conductivity (σ
a); and (3) soil water electrical conductivity (σ
w). The estimation of the σ
e of the soil using soil samples in the laboratory is the standard method
[34][127] for soil salinity measurements. This standard method is time-consuming, laborious, and costly for large-scale applications. Therefore, electrical resistivity, EMI, and TDR techniques are widely used as alternative methods. However, these alternative methods provide the σ
a of the subsurface but not the σ
w or σ
e [125][128][129][130][131][132][133]. The three main current flow pathways contributing to the σ are the liquid, solid–liquid, and solid phases
[129][134]. However, σ
w is the most appropriate measurement since it is the salinity experienced by plant roots. σ
w is impossible to obtain directly using alternative methods in the field, and it is difficult to measure in the laboratory
[128][129]. In the literature, soil salinity is expressed in different ways using terms such as σ
a, σ
w, σ
e, and σ at different soil: water (σ
1:1, σ
1:2, and σ
1:5) ratios. The σ of the solid phase (σ
s) is an important property in agricultural soils and a key variable in PA. It strongly correlates to clay content, a textural property that strongly influences soil water storage, dynamics, and plant growth. Together with SWC, σ
s significantly impacts the measured σa, making it an excellent surrogate for mapping clay content, usually with EMI
[128][129].
Researchers have identified that the influence of σ on GPR wave propagation limits GPR applications, including in agricultural soils. In most applications, such as locating buried objects or utilities and stratigraphic studies, a high σ restricts the GPR signal penetration due to attenuation
[28][34][36][45]. Soils with high clay content (illite and montmorillonite) have a high CEC and high σ, leading to higher GPR signal attenuation in clay-rich soils
[27][35][36]. However, this limitation (i.e., the influence of σ on GPR wave propagation) can be used as an opportunity in other applications, such as contamination mapping, identifying highly saline areas in PA, and seawater intrusion.
Mimrose et al.
[135] studied the influence of irrigation water salinity on the GPR signals using different salt-water concentrations. The authors showed that the amplitude of the GPR reflected wave is inversely proportional to irrigation water salinity. Ferrara et al.
[92] applied the early time GPR amplitude analysis method to find the influence of the σ on GPR wave amplitude under uniform ε
r conditions using salt-water. In this study, the authors found a high correlation between σ and GPR reflected wave amplitude. Alsharahi et al.
[136] used a numerical modeling approach to estimate reflections from iron bars and plastic bottles to evaluate the effect of ε
r and σ on GPR waves and found that the reflected wave amplitude decreases as σ increases. Wu et al.
[133] assessed σ variations by applying the waveform comparison method under different conductivities by implementing the GPR reflection coefficient method.
In previous studies, GPR was applied to identify and map soil contaminants such as excess fertilizers, soil amendments, leachates from dump sites, and hydrocarbons. Wijewardana et al.
[137][138][139] studied the effect of inorganic contaminants produced from landfill leachate on the GPR responses through field studies, controlled lysimeters, and a modeling approach. The authors found that GPR signal strength decreases with increasing contaminant concentrations due to increasing σ in contaminant plumes. The reflected wave disappeared completely at high σ levels due to attenuation
[139]. Reflected wave amplitude decreased with increasing σ, and further research in this field was suggested, since these methods could potentially evaluate σ variation in the subsurface and contamination mapping using the reflected wave amplitude
[135][139][140]. The R at an interface increases as the σ contrast increases. Nevertheless, as σ attenuates GPR waves in a given layer, it lowers the amplitude of the reflection amplitude for the lower interfaces. Under laboratory conditions, full-wave inversion was successfully used to retrieve soil conductivity (σ), as demonstrated by Lambot et al.
[38]. However, in field conditions with unknown subsurface layering, the inverse problem becomes ill-posed, and the retrieval of σ becomes challenging.
3.4. Soil Hydraulic Properties
As described by the water retention curve and hydraulic conductivity function, unsaturated soil hydraulic properties govern subsurface water dynamics
[115]. Hence, as GPR permits the characterization of SWC, time-lapse GPR offers the possibility to characterize these properties through the monitoring of SWC and its dynamics
[141][142]. This requires the coupling of the GPR derived-soil moisture or a GPR data processing algorithm with a soil hydrodynamic model, for example, one based on Richards’ equation
[143]. For instance, Binley et al.
[77], Rucker and Ferré
[82], Cassiani and Binley
[144], and Kowalsky et al.
[145] applied borehole GPR and tomographic inversion to monitor the distribution of water between boreholes and infer soil hydraulic properties. Lambot et al.
[38] remotely characterized the hydraulic properties of a laboratory soil column using full-wave GPR data inversion and subsequent soil hydrodynamic inversion. Lambot et al.
[84] introduced an integrated 3D full-wave electromagnetic and 1-D hydrodynamic inverse modeling procedure to estimate the soil hydraulic properties from far-field GPR measurements. The method was further studied and applied in the field by Jadoon et al.
[146][147]. Tran et al.
[106] used data assimilation techniques based on a maximum likelihood ensemble filter to estimate the soil hydraulic properties and reconstruct continuous soil moisture profiles. Despite the promising perspectives for environmental engineering applications shown in these studies, the utilization of joint GPR and hydrodynamic modeling approaches has proven to be challenging, particularly in the agricultural context, because the parameterization of these models is complex and requires a detailed understanding of the soil stratigraphy, petrophysical relationships, and boundary conditions for the hydrodynamic model.
3.5. Groundwater Table and Capillary Fringe Reflection
Determining the depth to the groundwater table (DGWT) is crucial in water management because DGWT affects groundwater recharge, water supply to plants, and contaminant accumulation and transport (especially agrochemicals)
[148][149]. Capillary fringe and the groundwater table fluctuate with seasonal variations, affecting agricultural water management, especially during the growing season. Indirect geophysical techniques such as GPR and seismic and resistivity techniques have been employed as alternative methods to traditional destructive piezometer installation when estimating DGWT. GPR is appropriate to estimate the DGWT of shallow aquifers (0–30 m) non-destructively on a large-scale
[145][150][151][152][153]. Information on the groundwater table fluctuation during the growing season is vital to understand the water availability for crops through capillarity and groundwater contamination potentials due to agricultural inputs.
The soil above (i.e., unsaturated) and the soil below (i.e., saturated) the water table have different SWCs and thus have different ε
r values
[31][148][151][154][155][156]. Therefore, due to the contrast in ε
r at the interface, the water table can be identified in the radargrams
[151][154][157]. Nevertheless, due to the capillary rise, the transition from the saturated zone to the unsaturated zone (capillary fringe) is not sharp, especially in fine-textured soils. Indeed, the observed reflection occurs some distance above the water table, depending on the shape of the capillary fringe. Under hydrostatic conditions, the shape corresponds to the soil’s water retention curve
[38]; otherwise, it can be relatively variable depending on the hydrodynamic conditions of the soil
[148][155][157][158][159][160]. The top of the capillary fringe is partly saturated and the bottom is fully saturated (the bottom is the water table). Thus, there is an SWC variation through the capillary fringe, and GPR wave velocity decreases from top to bottom (ε
r increases from top to bottom)
[115][156][160]. Because of this heterogeneity, GPR wave reflection varies, along with the capillary fringe
[148][151][156]. The height of the capillary fringe varies with the texture, pore size, and pore size distribution
[148][151][154][156][161]. In coarse grain soils (e.g., sand), the capillary height is less, and the contrast in ε
r between dry and saturated sand is sharp; consequently, the water table can be distinguished from easily GPR reflections
[153][161][162]. Conversely, the capillarity is high in fine-grain soils. Therefore, it is difficult to distinguish the actual water table in clay soils due to low contrast at the interface, decreasing the accuracy of DTWT estimation using GPR
[151][157][158][159][160][162][163].
3.6. Other Soil Properties
Other than SWC, soil compaction, and soil salinity, GPR has been applied to estimate other soil properties, such as soil texture and clay content
[164][165][166][167] and SOM/soil organic carbon (SOC)
[168][169][170][171][172]. Soil profile stratigraphy studies focusing on other soil properties were carried out by Doolittle and Collins
[173], Stroh et al.
[174], Meadows et al.
[164], André et al.
[175], and Nováková et al.
[176], and soil organic horizons were studied in particular by Winkelbauer et al.
[177].
In agricultural soil studies, GPR is extensively used for SWC estimation, followed by the estimation of soil salinity, porosity, and bulk density, while other properties are currently being researched. In addition to soil properties and states, GPR estimates soil horizons, stratigraphy, and water table mapping. The estimation of these properties in the agricultural landscape using GPR will provide essential information needed for farmland management to support PA.
This entry is adapted from the peer-reviewed paper 10.3390/rs15112932