Comparative Study of Ground Penetrating Radar Systems with Time Delay Methods for Accurate Pavement Thickness and Dielectric Constant Estimation

Abstract

The applicability of using ground penetrating radar (GPR) as a non-destructive technology for estimating the thickness and dielectric constant of pavement layers has been extensively documented in the literature. Time delay methods (TDMs) have been the most commonly used by researchers. However, various studies have reported different results with respect to the suitability of different GPR systems and TDMs for dielectric constant and thickness estimation. Thus, this study aims to compare the use of ground-coupled GPR and air-coupled GPR (AC-GPR) for estimating the dielectric constant of pavement materials. In the context of TDMs, the study examines three methodologies for assessing the traveled distance in the investigated layer. The study involves laboratory experiments using asphalt and dry sand. Also, it comprises a sensitivity analysis to examine the effect of different experimentation setups. The study employed two GSSI antennas: 51600S (GC-51600S) and GSSI 42000S (AC-42000S). Using the manufacturer’s recommended operation setups, the results generally indicated the possibility of simplifying the TDM calculation by considering the 42000S and 51600S antennas as monostatic for pavement materials with different dielectric constants (≥3) with an error rate within 5%. Moreover, the study indicates that using AC-GPR may be inappropriate for common midpoint analysis using the default mounting setups. These findings provide valuable insights into the suitability and limitations of different GPR systems and calculation methods for accurately estimating pavement’s dielectric constant and thickness.

Keywords

infrastructure highway maintenance pavements testing and evaluation of transportation structures ground penetrating radar testing

Ground penetrating radar (GPR) is a widely used non-destructive technology in various civil engineering applications. Initially, GPR was mainly used to locate buried objects in soil ( 1 ), for example, water, sewer, and gas pipelines. However, over the past three decades, significant improvements in GPR theory and technology have assisted in extending the application of GPR in civil engineering. Such applications include condition assessment and monitoring of concrete bridges ( 2 – 4 ), tunnels (5, 6), and pavement ( 7 – 9 ).

In particular, GPR is praised as a cost- and time-effective technology for various applications in pavement engineering. These include determining layer thicknesses, detecting pavement distresses, locating utilities, estimating density and compaction, and assessing air content. Notably, GPR has found widespread use in evaluating the layer thickness and density of in situ pavement materials, which usually necessitates the estimation of the dielectric constant ( 7 – 9 ).

The application principle of GPR is based on transmitting electromagnetic waves (EMs) and detecting the reflected parts at the interfaces of dielectric contrast. The information obtained from the received GPR signal, including the two-way travel time (TWTT) and amplitude, is utilized to characterize the properties of pavement materials. These properties are influenced by the material composition and affect the speed and deterioration of the wave. Among these properties, the dielectric constant is deemed the most crucial material property for evaluating layer thickness of in situ pavement materials (7, 8).

Several methods have been developed to estimate the dielectric constant of pavement materials. One approach involves determining the dielectric constant using numerical methods based on the composition of the pavement materials, including binding materials, aggregates, air, and water content ( 8 ). The surface reflection method (SRM) is used to estimate the dielectric constant for pavement layers using information on the amplitude of the reflected wave from pavement layer interfaces ( 10 ). However, researchers have revealed that the measured dielectric constant values are likely to experience significant uncertainties because of various factors ( 11 –13). Primarily, the SRM assumes homogeneity of the pavement layer through its depth. Consequently, the accuracy of this method can be compromised if the properties of the pavement material are inconsistent throughout its depth. To address this limitation of the SRM, time delay methods (TDMs) are commonly used to estimate the dielectric constant. The application of TDMs is discussed in the following section.

Dielectric Constant and Thickness Calculation Using Time Delay Methods

TDMs are commonly used to estimate the dielectric constant based on the TWTT of the GPR wave in the surveyed pavement material. TDMs are based on the relationship established between the GPR wave’s velocity and the medium’s dielectric constant ( 9 ). The related calculation is shown in Equations 1–3 ( 9 ). Equation 1 establishes the relationship among the velocity of the GPR wave in layer n ( $v_{n}$ ), the traveled distance across the layer ( $z_{n}$ ), and the TWTT $t_{n}$ in the same layer. Equation 2 defines the velocity $v_{n}$ as a ratio between the velocity of the EM wave in air (c) to the square root of the dielectric constant ( $ε_{n}$ ) of layer n. By substituting Equation 1 into Equation 2 and rearranging the terms, Equation 3 can be derived for calculating the dielectric constant of layer n:

v_{n} = \frac{2 z_{n}}{t_{n}}

(1)

v_{n} = \frac{c}{\sqrt{ε_{n}}}

(2)

ε_{n} = {(\frac{t_{n} \times c}{2 \times z_{n}})}^{2}

(3)

The classical form of the TDM involves the use of sampled cores to estimate the traveled distance $z_{n}$ for the estimation of the dielectric constant ( 12 ). The estimated dielectric constant, determined at discrete locations, is then extrapolated for the entire study area. Then, Equation 3 can be rearranged and utilized to calculate the thickness of the whole layer.

This TDM is more representative of the dielectric constant of the whole layer thickness compared to the SRM. Loizos and Plati ( 12 ) compared the SRM and TDMs to calculate the dielectric constant in 16 test locations selected within nine pavement sections. The TDM achieved a median error of 5%, while the SRM had a notably higher median error of 9.6%. However, obtaining the TWTT information for thin layers can be challenging because of signal overlapping. Also, the availability of layer thickness information necessitates the extraction of core samples, which can be time-consuming and expensive. Moreover, the extracted core samples may not fully represent the entire study area because of the presence of variations in compaction, void content, water content, and other influential factors.

Alternatively, the common midpoint (CMP) method is a widely employed technique with GPR to estimate the thickness and the dielectric constant of pavement layers without the need for core extraction ( 9 ). The application of the CMP involves the separation of the transmitter and the receiver of the GPR system at different distances. Figure 1 illustrates the application of the CMP in the first layer. Subsequently, triangulation principles are leveraged to calculate the thickness of the layers and dielectric constant ( 14 ).

Figure 1.

Application of the common midpoint method.

The CMP method offers several advantages, including applicability to estimating the dielectric constant and pavement thickness without the need for core extraction. In addition, it relies on time delay, which makes it less prone to variations throughout the thickness of the layer. However, it should be noted that the CMP may not be effective for thin layers because of the challenges posed by signal overlapping. Furthermore, the performance of the method could be hindered when using ground-coupled GPR (GC-GPR) because of ground–antenna coupling challenges and limited resolution. On the other hand, the use of air-coupled GPR (AC-GPR) may pose limitations on the applicability of the CMP considering the conventional mounting setups, as demonstrated in the subsequent sensitivity analysis conducted in this study (the recommended mounting height for widely used AC-GPR systems from GSSI is typically between 46 and 51 cm [ 15 ]).

Proper consideration of the surveyed case, particularly the thickness range and dielectric constant of the surveyed layer, is crucial when selecting the appropriate GPR system and calculation method for dielectric constant estimation in the pavement. In the context of the TDM, assessing the traveled distance by the EM wave may become a critical factor, and researchers have approached this aspect differently in various applications. Some studies opt for a simplified approach by treating commercially available bistatic GPR systems as monostatic, which can ease the estimation process. However, other researchers explore triangulations to simulate the traveled path of the EM waves, aiming to enhance the accuracy of results or facilitate the utilization of AC-GPR for the CMP.

Monostatic Simplification

Most commercially available GPR systems are bistatic, featuring fixed transmitters and receivers within a single unit. However, in practice, the distance is often calculated by treating the system as monostatic. Figure 2 illustrates the traveled distance considered for the monostatic simplification (MS) method when using TDMs. In such cases, the traveled distance across the first layer ( $z_{1}$ ) is assumed to be equal to thickness of the first layer ( $h_{1}$ ).

Figure 2.

Determination of the traveled distance using the monostatic simplification method.

Bistatic Simplification

In the bistatic simplification (BS) method, the bistatic nature of the GPR system is recognized. However, Snell’s law is ignored when estimating the traveling distance of the signal. Snell’s law explains the behavior of the EM waves when traveling across two materials of different dielectric constants ( 9 ). Figure 3 illustrates the travel path of the GPR signal reflected from the bottom of the first layer. As shown in Figure 3, the traveled distance can be linked to the thickness of the layers using the Pythagoras theorem.

Figure 3.

Determination of traveled distance using the bistatic simplification method.

For estimating the traveled distance by the GPR wave in the first layer, the triangulation relationships can be applied, as shown in Equations 4 and 5. Equation 6 can be obtained by substituting Equation 4 into Equation 5 ( 5 ). Using Equation 6, $z_{1}$ can be calculated as a function of the transmitter–receiver offset distance ( $d_{TR}$ ), thickness of the first layer ( $h_{1}$ ), and the elevation of the transmitter/receiver from the bottom of the investigated layer ( $h$ ):

\frac{h_{1}}{h} = \frac{w}{\frac{d_{T R}}{2}}

(4)

z_{1} = \sqrt{w^{2} + {(h_{1})}^{2}}

(5)

z_{1} = \sqrt{{(\frac{h_{1} \times (d_{TR / 2})}{h})}^{2} + {(h_{1})}^{2}}

(6)

where $w$ is the offset distance between the incident point on the surface layer (q) and the midpoint distance between the transmitter and the receiver (m), as shown in Figure 3.

The BS method was utilized in the work of Ahmed et al. ( 16 ), where AC-GPR was used for pavement thickness estimation. The study claimed a significant improvement in the accuracy of thickness estimation. Hamrouche and Saarenketo ( 17 ) proposed a method for evaluating the dielectric constant by employing a horizontally arranged single transmitter and multiple receivers of the AC-GPR type. The study utilized BS to estimate the travel path of the GPR waves and apply CMP calculations.

Bistatic Correction

Some studies have recognized the transmitter–receiver offset distance as well as the effect of signal refraction, as outlined in Snell’s law to estimate the path of the GPR waves. Figure 4 provides an illustration of the travel path of the GPR signal reflected from the bottom of the first layer according to the bistatic correction (BC) method.

Figure 4.

Determination of layer thickness using the bistatic correction method.

Equation 7 defines the relationship between the traveled distance in the air with the elevation of the transmitter/receiver from the top surface of the investigated layer and the distance $x_{1} .$ Equation 8 establishes the relationship between the traveled distance in the upper pavement layer ( $z_{1})$ with the thickness of the layer $(h_{1})$ and the distance $x_{0}$ . Here, $x_{0}$ is the offset distance between the incident point on the surface layer (q) and the midpoint distance between the transmitter and the receiver (m), as shown in Figure 4. Equation 9 relates $x_{0}$ , $x_{1}$ and the transmitter–receiver offset distance ( $d_{TR})$ . Equation 10 demonstrates the application of Snell’s law by relating the incident angle $(α)$ and the refraction angel ( $θ)$ with the velocity of the GPR wave in the air ( $v_{0})$ and the upper pavement layer $(v_{1})$ ( 14 ). The sine of the incident and the refracted angles can be estimated using Equations 11 and 12, respectively. Equation 13 is obtained by combining Equations 2, 7, 8, 9, 10, and 11. Likewise, Equation 14 is derived by substituting Equation 8 into Equation 2. Subsequently, $x_{0}$ can be estimated by solving Equations 13 and 14. The traveled distance ( $z_{1}$ ) and dielectric constant ( $ε_{1})$ of the upper pavement layer can be estimated using Equations 9 and 14 ( 14 ), respectively:

y = \sqrt{{(x_{1})}^{2} + {(h_{air})}^{2}}

(7)

z_{1} = \sqrt{{(x_{0})}^{2} + {(h_{1})}^{2}}

(8)

d_{TR / 2} = x_{0} + x_{1}

(9)

\frac{v_{0}}{v_{1}} = \frac{\sin (α)}{\sin (θ)}

(10)

\sin (α) = \frac{x_{1}}{y}

(11)

\sin (θ) = \frac{x_{0}}{z_{1}}

(12)

ε_{1} = {(\frac{d_{TR} - x_{0}}{\sqrt{{(d_{TR} - x_{0})}^{2} + {(h_{air})}^{2}}} \times \frac{\sqrt{{(x_{0})}^{2} + {(h_{1})}^{2}}}{x_{0}})}^{2}

(13)

ε_{1} = {(\frac{t_{1} \times c}{2 \times \sqrt{{(x_{0})}^{2} + {(h_{1})}^{2}}})}^{2}

(14)

Lahouar et al. ( 10 ) developed a modified CMP method using a monostatic GC-GPR and bistatic GC-GPR. This method accounts for signal refraction through pavement layers by incorporating Snell’s law. Leng and Al-Qadi ( 9 ) proposed a modified common midpoint (XCMP) technique employing two air-coupled bistatic systems. In the XCMP method, the two systems are aligned along the surveying direction, and each receiver is assumed to collect signals from both transmitters. The XCMP method incorporates Snell’s law to address signal refraction ( 9 ).

Research Motivation and Objectives

The applicability of the GPR to measure the dielectric constant and layer thickness of pavement has been extensively documented. However, researchers have reported various performance results based on the surveyed site, the used GPR system, and the employed estimation method ( 9 ). This study aims to provide valuable insights to aid researchers and practitioners in selecting the most suitable method and GPR system for their specific study, considering the type and thickness of the investigated material.

The objectives of this study include comparing the applicability of AC-GPR and GC-GPR to evaluate the dielectric constant and thickness of pavement materials using TDMs. By examining the performance of these two GPR systems, the study aims to identify their respective strengths and limitations. In addition, this study aims to compare the different methods for estimating the distance traveled by GPR waves in the context of TDMs. Understanding the impact of these distance estimation techniques on the applicability of AC-GPR for the CMP and the overall accuracy of dielectric constant and thickness estimation is a key aspect of this investigation.

Methodology

The investigation was conducted in the lab using two GPR systems (AC-51600S and GC-51600S) and two pavement materials: asphalt and dry sand. The laboratory setup for the experiment is shown in Figure 5. Furthermore, three methods for evaluating the traveled distance in the pavement layer (MS, BS, and BC) were compared in the context of TDM calculations. In addition, sensitivity analysis was conducted to investigate the effect of GPR system characteristics and mounting setups on the accuracy of dielectric estimation and their applicability for CMP calculations. This section outlines the calculation methods for the dielectric constant and the description of the experimental setups.

Figure 5.

Lab setup showing the tested the materials (A50, A100, and DS75) and the AC-51600S antenna.

Materials

Two types of pavement materials, namely asphalt and dry sand, were used to compare the different GPR systems and dielectric constant estimation methods. The type and thickness of the pavement layers investigated in this study are presented in Table 1. Specimen A50 represents an asphalt layer with a thickness of 50 mm, whereas A100 represents an asphalt layer with a thickness of 100 mm. Both A50 and A100 were prepared using the same mix design to ensure consistent properties for comparative analysis. In addition, the study included DS75, which represents a dry sand layer with a thickness of 75 mm. The dimensions of the specimens are 30 cm × 60 cm.

Table 1.

Pavement Materials Used in the Experiment

Material	Thickness (mm)	Specimen name
Asphalt	50	A50
Asphalt	100	A100
Dry sand	75	DS75

GPR System

This study utilized two GPR systems for data collection: an AC-51600S system and a GC-51600S system, both manufactured by GSSI (Geophysical Survey Systems, Inc.). The AC-51600S system comprised a SIR 30 control unit and a GSSI 42000S high-frequency air launched antenna. The antenna used in this system had a central frequency of 2.0 GHz and dimensions of 21 cm × 55.6 cm × 49.5 cm (width, length, and height). The distance between the centroid of the transmitter and the receiver is 37.5 cm ( 16 ). The elevation of the transmitter and receiver from the soffit of the antenna box was estimated at 13 cm. The ground-coupled system consisted of the SIR-4000 control unit and the GSSI ground-coupled antenna model 51600S. The antenna used in this system had a central frequency of 1.6 GHz and dimensions of 3.8 cm × 10 cm × 16.5 cm (thickness, width, and length). The distance between the centroid of the transmitter and the receiver is around 6.2 cm. During regular use of GC-51600S, the estimated distance between the transmitter and receiver from the surface of the examined material was assumed to be within 2 cm. For system operation, RADAN software (version 7.5.18.02270) was employed. The collected data were exported in CSV format for further processing and analysis. Figure 6, a–e, visually represents the SIR-30 unit, 42000S antenna, SIR-4000 system, 51600S unit, and RADAN software interface, respectively.

Figure 6.

Ground penetrating radar systems: (a) GSSI 51600S antenna, (b) SIR 30 control unit, (c) GSSI 51600S antenna, (d) SIR 4000 unit, and (e) RADAN interface.

The setup of the GPR systems used in this study is presented in Table 2. The number of samples per scan was set as 512. The time range was set as 12 ns for the air-coupled system and 7.65 ns for the ground-coupled system. The data collection mode was configured as time. Moreover, infinite impulse response (IIR) and finite impulse response (FIR) filters were applied to improve the quality of the collected data. The IIR filter of a high-pass cutoff frequency of 100 MHz was applied to both systems. A custom FIR filter, GSSI 42000S, was applied for the AC-51600S. A FIR filter with a low-pass cutoff frequency of 3000 MHz and a high-pass cutoff frequency of 470 MHz was applied to the GC-51600S.

Table 2.

Data Collection Setting for the Air- and Ground-Coupled Systems

Attribute	AC-GPR	GC-GPR
Samples/scan	512	512
Time range (ns)	12	7.65
Scan mode	Time	Time
IIR	High-pass 100 (MHz)	High-pass 100 (MHz)
FIR	GSSI 42000S custom filter	High-pass 470 (MHz) Low-pass 3000 (MHz)

Note: AC-GPR = air-coupled ground penetrating radar; GC-GPR = ground-coupled ground penetrating radar; IIR = infinite impulse response; FIR = finite impulse response.

Experimental Setup

The experiment involves collecting GPR data for the three samples after placing them on the top of a metal plate (MP). The MP was approximately the same size as the specimens. Figures 7 and 8 depict the experimental setup and GPR reading using the AC-42000S and GC-51600S for A50, respectively. Furthermore, Figures 7 and 8 present the A-scan and B-scan acquired during the experiment on specimen A50 using the AC-51600S and GC-51600S, respectively.

Figure 7.

Illustration diagram of the experimental setup using the AC-51600S for specimen A50.

Figure 8.

Illustration diagram of the experimental setup using the GC-51600S for specimen A50.

Sensitivity Analysis

This study involved laboratory experimental work with a limited number of samples and setups. Thus, a sensitivity analysis was conducted to understand the impact of the investigated pavement material’s thickness and dielectric constant, antenna parameters, and the experiment setup on the criticality of addressing the transmitter–receiver offset. In addition, the analysis was used to investigate the applicability of using the AC-51600S for the CMP. The sensitivity analysis was conducted using a theoretical framework. By systematically examining these factors and their effects on GPR measurements, the sensitivity analysis offered valuable information to guide researchers and practitioners in optimizing GPR methodologies for accurate dielectric constant estimation in pavement engineering. The details of the parameters considered in the sensitivity analysis are shown in Table 3.

Table 3.

Details of the Sensitivity Analysis

S.N.	GPR system	Variables	Values
1	AC-GPR GC-GPR	$h_{1}$ (cm)	10, 15, 20, 30, 40, 50, 60
		Dielectric constant	1–20
		Calculation method of $z_{0}$	MS, BS, BC
2	AC-GPR	$h_{1}$ (cm)	10, 20, 30, 40, 50
2	AC-GPR	$d_{TR}$ (cm)	20–200

Note: GPR = ground penetrating radar; AC-GPR = air-coupled ground penetrating radar; GC-GPR = ground-coupled ground penetrating radar; MS = monostatic simplification; BS = bistatic simplification; BC = bistatic correction.

Results and Discussion

Experimental Results

Table 4 presents the calculated dielectric constant values of the tested specimens using different TDMs with the AC-42000S and GC-51600S. The dielectric constant values were estimated by employing MS, BS, and BC methods to assess the GPR wave’s traveled distance. Generally, dielectric constant values measured with the AC-42000S exhibited greater consistency compared to those measured with the GC-51600S, particularly for specimens A50 and DS75. This can be attributed to the shorter wavelength and higher resolution of the AC-42000S in contrast to the GC-51600S. In addition, dielectric constant values obtained using individual TDMs for the asphalt specimens (A50 and A100) were notably consistent, especially with the AC-42000S. For instance, dielectric constant values measured using the different methods for the two asphalt specimens exhibited only 0.1 difference with the AC-42000S. Similarly, the application of MS on the asphalt specimens yielded a highly consistent dielectric constant calculation using the GC-51600S, showing only a 0.1 difference. The relatively high divergence between the dielectric constant values calculated for the two asphalt specimens (A50 and A100) using the GC-51600S while employing the BS and BC methods can be attributed to the limited thickness of the A50 specimens compared to the wavelength of the GPR wave.

Table 4.

Dielectric Constant of the Tested Specimens

Setup no.	Specimen	Dielectric constant
Setup no.	Specimen	GC-MS	GC-BS	GC-BC	AC-MS	AC-BS	AC-BC
1	A50	3.7^a	3.1^a	3.4^a	4.2	3.8	4.1
2	A100	3.8	3.6	3.7	4.3	3.9	4.2
3	DS75	1.8^b	1.6^b	1.7^b	2.7	2.4	2.6

Note: GC = ground-coupled; AC = air-coupled; MS = monostatic simplification; BS = bistatic simplification; BC = bistatic correction.

Suspected corrupted value caused by wave overlapping.

Severely corrupted value caused by wave overlapping.

However, a different observation arises when comparing the dielectric constant values obtained from the two employed GPR systems using the same TDM for the same material. Figure 9, a and b, illustrates the variation between the dielectric constant values calculated using the three TDMs and the two GPR systems for specimens A50 and A100, respectively. The figures illustrate that the dielectric constant values obtained using the AC-42000S were larger than those obtained using the AC-51600S by about 18% and 12% for specimens A50 and A100, respectively. This suggests that the antenna type and associated settings have a more pronounced effect compared to the effect of the TDM, particularly for MS and BC.

Figure 9.

Comparison of the dielectric constant values calculated using the different combinations of the employed time delay method and ground penetrating radar systems for samples (a) A50 and (b) A100.

The dielectric constant measured for DS75 using the AC-42000S was approximately 2.6, while it was significantly lower at 1.7 when measured with the GC-51600S. This difference can be attributed to the overlapping of waves when using the GC-51600S. Figure 10a shows the combined wave obtained in RADAN, representing a superposition of the direct and the reflected waves. The reflected wave was obtained by reducing the direct wave, as shown in Figure 10a. Subsequently, the wave representing the air/DS75 interface was reduced, as shown in Figure 10b. As presented in Figure 10b, the first cycle of the reflection from the DS75/MP interfaces with the second cycle of the reflection wave from the air/DS75 interface. This overlap resulted in a shift in the peak of the first cycle of the wave reflection from the DS75/MP by three samples. After removing the overlapping effect, the corrected value of the dielectric constant for DS75 was obtained as 2.0. However, a detailed discussion on how to obtain the corrected peak by subsequent reduction of the waves is beyond the scope of this study. Thus, the TDM can be characterized as unsuitable for evaluating the dielectric constant of DS75 in the context of this study. The overlapping of waves in the GC-51600S led to inaccuracies in the dielectric constant estimation, highlighting the limitation of TDMs in the evaluation of the dielectric constant of materials of low depth and dielectric constant. The lower dielectric constant values measured for A50 compared to A100 using the GC-51600S can be primarily attributed to wave interference effects. The relatively long wavelength of the GC-51600S, when investigating A50, led to increased interference between the reflected waves. Nevertheless, the dielectric constant value for A50 is believed to be minimally affected by wave overlapping, as it is relatively close to that measured for the A100 specimen.

Figure 10.

Subsequent decomposition of the wave obtained for testing DS75 using the GC-51600S: (a) decomposition of the combined wave into direct and reflected waves and (b) decomposition of the reflected wave into DS75 reflection and resultant waves.

Indeed, overlapping of reflected waves can pose challenges when measuring the thickness or evaluating the dielectric constant of thin pavement layers with a low dielectric constant. Thin layers can be defined as layers where the GPR reflections at the surface and bottom overlap. In other words, thin layers refer to those with a thickness smaller than the wavelength of the GPR waves when passing through them. Considering that the current wave of GPR consists of two cycles, the shift of the peak in the first cycle of the wave from the bottom of one layer may happen when the thickness of the layer is smaller than three-quarters of the wavelength of the GPR waves when passing through the same layer. The wavelength of the predominant frequency of a GPR of the central frequency (f) in a layer can be calculated using the frequency formula, as shown in Equation 15:

λ = \frac{c}{f \times \sqrt{ε}}

(15)

Table 5 presents the wavelength of the GC-51600S and AC-42000S waves in asphalt and dry sand as calculated using Equation 15. As presented in Table 5, the layer thicknesses of specimens A50 and DS75 are less than three-quarters of the wavelength of the GC-51600S, indicating overlapping that may result in discrepancies in the measured dielectric constant values.

Table 5.

Wavelength for the Ground Penetrating Radar (GPR) Signals in the Tested Materials

Material	h (cm)	$ε$	v (cm/ns)	1.6 GHz>		2.0 GHz
Material	h (cm)	$ε$	v (cm/ns)	$λ$ (cm)	$0.75 λ$ (cm)	$λ$ (cm)	$0.75 λ$ (cm)
Asphalt (A50)	5	4.2	14.8	9.1	6.9	7.3	5.5
Asphalt (A100)	10	4.2	14.8	9.1	6.9	7.3	5.5
Dry sand (DS75)	7.5	2.6	18.6	11.6	8.7	9.3	7.0

Note: h = specimen thickness; $ε$ = dielectric constant; $λ$ = GPR wavelength.

Considering the setting of the GPR systems, as shown in Table 2, the resolution effect on the obtained results is investigated. As the GPR readings are collected in steps called samples, the start and end of the layer can cumulatively involve uncertainty of up to almost two scans. In such a way, the uncertainty range can be within ±2 samples. The dielectric constant uncertainty was calculated using the AC-42000S and the GC-51600S for the different specimens.

Figure 11 presents the uncertainty calculated for the dielectric constant as a percentage of the estimated dielectric constant. As shown in Figure 11, when comparing the results individually for each GPR system, the uncertainty range is more extensive for A50 compared to that of A100 for both the AC-42000S and GC-51600S. In general, it is expected that the larger the thickness and the dielectric constant, the lower the effect of the resolution on the uncertainty in the calculated value of the dielectric constant. Also, Figure 11 demonstrates that the uncertainty for the GC-51600S is lower compared to that of the AC-42000S when comparing the values for the same specimens. This can be attributed to the higher resolution of the GC-51600S because of the setup used. The time range was set at 7.65 ns for the GC-51600S compared to 12 ns for the AC-42000S, while the number of samples per scan (A-scan) remained the same for both systems.

Figure 11.

Error range for the calculated dielectric constant caused by resolution.

The resolution of the GPR system can be adjusted by manipulating the number of samples per scan and the targeted time range. Increasing the targeted time range and reducing the number of samples can have a negative impact on the vertical resolution. As a consequence, uncertainties with respect to the calculation of the dielectric constant and pavement thickness increase. Moreover, when expressed as a ratio, the uncertainties can be influenced by other factors, including the thickness of the surveyed layer and the dielectric constant. Thus, to enhance the accuracy of GPR-based measurements, careful consideration of the resolution settings and a thorough understanding of the factors affecting uncertainties is crucial.

Sensitivity Analysis

A comprehensive sensitivity analysis was conducted to understand the impact of the dielectric constant and thickness of the surveyed material on the error introduced by employing the simplification methods (MS and BC). The analysis considered the same mounting and data collection setup, calculating the error across a range of seven different thicknesses and dielectric constant values, spanning from 1 to 20.

Figure 12, a and b, plots the relative error introduced in the calculated dielectric constant caused by the adoption of the simplification methods of MS and BS, respectively. As shown in Figure 12a, the use of the MS method resulted in an overestimation of the dielectric constant values. Even for materials with a dielectric constant equal to two, the error is below 5%, which is much smaller than the error that might originate from the resolution problems of the GPR in the settings used in the current study. Thus, the MS method can be applied with a negligible error rate in such scenarios. In contrast, Figure 12b shows that the adoption of the BS method caused an underestimation of the dielectric constant values for the analyzed scenarios. Moreover, the absolute relative error value increased with the increase in the dielectric constant and layer thickness. If the dielectric constant is larger than three or the thickness is lower than 20 cm, the error introduced using the MS is limited and lower than that introduced using BS. Therefore, in such scenarios, it is reasonable to consider the GPR system as monostatic.

Figure 12.

Introduced errors in the calculation of the dielectric constant using the AC-42000S based on the (a) monostatic simplification (MS) and (b) bistatic simplification (BS) methods.

Figure 13, a and b, shows the relative error introduced using the MS and BS methods with the GC-51600S. Similar trends to those observed for the AC-42000S are seen for the GC-51600S, with the relative error being lower for the GC-51600S. Except for small layer thicknesses (≤15 cm), the relative error is lower than 2% when using the MS method. In addition, the error introduced by using MS or BS is negligible for higher thicknesses, indicating the possibility of using the GC-51600S as a monostatic system.

Figure 13.

Introduced errors in the calculation of the dielectric constant using the GC-51600S based on the (a) monostatic simplification (MS) and (b) bistatic simplification (BS) methods.

As previously discussed, when the receiver and transmitter are elevated at approximately 40 cm with a transmitter–receiver offset of 37.5 cm (the setup used in the current study), they can be treated as monostatic with negligible error, falling within the resolution error uncertainty. Therefore, we analyzed the effect of antenna elevation and transmitter–receiver offset to investigate the applicability of using AC-GPR for the CPM calculation. The analysis considered a pavement material with a dielectric constant equal to 4.2.

Figure 14, a and b, illustrates the relationship between the distance $2 x_{0}$ and the transmitter–receiver offset at two elevations of 50 and 20 cm, respectively. Here, $2 x_{0}$ represents the distance between the incident and reflection points on the surface of the tested layer.

Figure 14.

Relationship between the transmitter–receiver offset of air-coupled ground penetrating radar and the distance $2 x_{0}$ when the antenna elevation equals (a) 50 cm and (b) 20 cm.

Figure 14a shows that at an elevation of 50 cm and $d_{TR}$ less than 70 cm, the value of $2 x_{0}$ is expected to be less than 20 cm. Particularly for lower thicknesses (10 and 20 cm), the value of $2 x_{0}$ remains less than 20 cm, even after spreading the transmitter and the receiver by 2 m. On the other hand, Figure 14b indicates that lowering the elevation of the antenna to 20 cm can boost the value of $2 x_{0}$ . However, for low thicknesses (10 and 20 cm), $2 x_{0}$ does not exceed 20 cm, even after spreading the transmitter and the receiver by about 2 m. Only spreading the transmitter and the receiver beyond 80 and 100 cm at an elevation of 20 cm can result in $2 x_{0}$ values larger than 30 cm when investigating layers of thickness 40 and 50 cm, respectively.

This analysis indicates that the practicality of using AC-GPR for the CMP method should be carefully evaluated based on the thickness of the layer and the expected dielectric constant. The difference in $2 x_{0}$ values needs to be large enough between the two readings (corresponding to the first and second spread distance) to exceed the uncertainty range introduced by the resolution error. For lower thicknesses (less than 20 cm) and relatively low dielectric constants (less than 4), it may not be practical to use AC-GPR, even at an elevation of 20 cm and a transmitter–receiver offset of 2 m.

Conclusions

The extensive literature on GPR has established its applicability as a non-destructive technology for estimating the thickness and dielectric constant of pavement layers. Researchers have widely used TDMs for this purpose. However, the variability in results across different studies with respect to the suitability of different GPR systems and TDM calculation methods for dielectric constant and thickness estimation calls for a comprehensive investigation.

This study involved a comparative analysis of GC-GPR and AC-GPR for estimating the dielectric constant of pavement materials. Within the framework of TDMs, the study examined three distinct methodologies for estimating the traveled distance in the investigated layer. The investigations encompassed laboratory experiments using asphalt and dry sand and simulations with various experimentation setups. The GSSI 51600S and GSSI 42000S antennas were utilized to acquire data for our study while adopting the default setups recommended by the vendor.

The results indicated that the AC-42000S and GC-51600S yielded dielectric constants that differed by about 12% for A100. Notably, the differences were more significant for the other two specimens with lower thickness, which points to a limitation of TDMs when dealing with thin layers. Also, the effect of resolution on the certainty of the results is more prominent for lower thicknesses and dielectric constants. The results revealed that 42000S and 51600S antennas could be dealt with as monostatic antennas for pavement materials with different dielectric constants (≥3) with an error rate within 5%. However, this study highlights the need for caution when employing AC-GPR in common CMP analysis with default mounting setups, as it may yield unreliable results.

By providing a comprehensive comparison of GC-GPR and AC-GPR, together with insights into optimizing TDM calculations, this study provides practical guidance for researchers and practitioners in making informed decisions with respect to GPR system selection and calculation methods. This study adds to the existing body of knowledge about applications of GPR in pavement engineering and emphasizes the importance of customized methods in specific survey conditions. However, this study was limited to laboratory experiments employing only two GPR systems and three specimens, supplemented with a sensitivity analysis conducted using a theoretical framework. The impact of the physical sample size, particularly for the application of AC-GPR, was not investigated. Therefore, further research is needed to complement this study by incorporating more GPR systems, specimens of different types, and varying physical dimensions.

Footnotes

Author Contributions

The authors confirm contribution to the paper as follows: study conception and design: A. Fares, T. Zayed; data collection: A. Fares, N. Faris; analysis and interpretation of results: A. Fares, N. Faris; draft manuscript preparation: A. Fares, T. Zayed, N. Faris, S. Abdelkhalek, E.M. Abdelkader. All authors reviewed the results and approved the final version of the manuscript.

Declaration of Conflicting Interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: Ali I.F. Fares would like to express his acknowledgment and thanks to the Research Grants Council (RGC) of Hong Kong for supporting his PhD study through the Hong Kong Ph.D. Fellowship Scheme (HKPFS).

ORCID iDs

Ali Fares

Tarek Zayed

Nour Faris

Sherif Abdelkhalek

References

Alsharqawi

Dawood

Abdelkhalek

Abouhamad

Zayed

Condition Assessment of Concrete-Made Structures Using Ground Penetrating Radar. Automation in Construction, Vol. 144, 2022, p. 104627. https://doi.org/10.1016/j.autcon.2022.104627.

Abdelkader

E. M.

Zayed

Faris

Synthesized Evaluation of Reinforced Concrete Bridge Defects, Their Non-Destructive Inspection and Analysis Methods: A Systematic Review and Bibliometric Analysis of the Past Three Decades. Buildings, Vol. 13, No. 3, 2023, p. 800. https://doi.org/10.3390/buildings13030800.

Abdelkader

Zayed

Comprehensive Inspection System for Concrete Bridge Deck Application: Current Situation and Future Needs. Journal of Performance of Constructed Facilities, Vol. 34, No. 5, 2020, p. 03120001. https://doi.org/10.1061/(asce)cf.1943-5509.0001484.

Faris

Zayed

Abdelkader

E. M.

Fares

Corrosion Assessment Using Ground Penetrating Radar in Reinforced Concrete Structures: Influential Factors and Analysis Methods. Automation in Construction, Vol. 156, 2023, p. 105130. https://doi.org/10.1016/j.autcon.2023.105130.

Qin

Zhang

Tang

Wang

Automatic Recognition of Tunnel Lining Elements from GPR Images Using Deep Convolutional Networks with Data Augmentation. Automation in Construction, Vol. 130, 2021, p. 103830. https://doi.org/10.1016/j.autcon.2021.103830.

Dawood

Zhu

Zayed

Deterioration Mapping in Subway Infrastructure Using Sensory Data of GPR. Tunnelling and Underground Space Technology, Vol. 103, 2020, p. 103487. https://doi.org/10.1016/j.tust.2020.103487.

Evans

Frost

Stonecliffe-Jones

Dixon

Assessment of In Situ Dielectric Constant of Pavement Materials. Transportation Research Record: Journal of the Transportation Research Board, 2007. 2037: 128–135.

Cao

Al-Qadi

I. L.

Development of a Numerical Model to Predict the Dielectric Properties of Heterogeneous Asphalt Concrete. Sensors, Vol. 21, No. 8, 2021, p. 2643. https://doi.org/10.3390/s21082643.

Leng

Al-Qadi

I. L.

An Innovative Method for Measuring Pavement Dielectric Constant Using the Extended CMP Method with Two Air-Coupled GPR Systems. NDT & E International, Vol. 66, 2014, pp. 90–98. https://doi.org/10.1016/j.ndteint.2014.05.002.

10.

Lahouar

Al-Qadi

I. L.

Loulizi

Clark

T. M.

Lee

D. T.

Approach to Determining In Situ Dielectric Constant of Pavements: Development and Implementation at Interstate 81 in Virginia. Transportation Research Record: Journal of the Transportation Research Board, 2002. 1806: 81–87.

11.

Porubiaková

Komačka

A Comparison of Dielectric Constants of Various Asphalts Calculated from Time Intervals and Amplitudes. Procedia Engineering, Vol. 111, 2015, pp. 660–665. https://doi.org/10.1016/j.proeng.2015.07.129.

12.

Loizos

Plati

Accuracy of Pavement Thicknesses Estimation Using Different Ground Penetrating Radar Analysis Approaches. NDT & E International, Vol. 40, No. 2, 2007, pp. 147–157. https://doi.org/10.1016/j.ndteint.2006.09.001.

13.

Marecos

Fontul

Solla

de Lurdes Antunes

Evaluation of the Feasibility of Common Mid-Point Approach for Air-Coupled GPR Applied to Road Pavement Assessment. Measurement, Vol. 128, 2018, pp. 295–305. https://doi.org/10.1016/j.measurement.2018.06.062.

14.

Al-Qadi

I. L.

Lahouar

Use of GPR for Thickness Measurement and Quality Control of Flexible Pavements. Journal of the Association of Asphalt Paving Technologists, Vol. 73, 2004, pp. 501–528.

15.

Geophysical Survey Systems, Inc. GSSI Antenna Manual. Geophysical Survey Systems, Inc., Nashua, NH, 2017.

16.

Ahmed

M. U.

Tarefder

R. A.

Maji

A. K.

Incorporating Transmitter–Receiver Offset to Interpret Pavement Layer Thicknesses by GPR. Case Studies in Nondestructive Testing and Evaluation, Vol. 6, 2016, pp. 94–104. https://doi.org/10.1016/j.csndt.2016.11.002.

17.

Hamrouche

Saarenketo

Improvement of a Coreless Method to Calculate the Average Dielectric Value of the Whole Asphalt Layer of a Road Pavement. Proceedings of the 15th International Conference on Ground Penetrating Radar, GPR 2014, Brussels, Belgium, IEEE, New York, 2014, pp. 899–902. https://doi.org/10.1109/ICGPR.2014.6970556.