Search

In this study, we analyzed low-frequency drone-borne ground-penetrating radar (GPR) and full-wave inversion for soil electrical conductivity mapping. Indeed, in the lowest GPR frequency ranges, the soil surface reflexion coefficient depends more on the soil electrical conductivity than on its permittivity. Numerical experiments were conducted within the frequency range 15–45 MHz to analyze parameter sensitivities, the well-posedness of the inverse problem as well as the depth of sensitivity. The results show that the soil surface reflexion is significantly more sensitive to the soil electrical conductivity than the soil permittivity. Therefore, the conductivity can be retrieved using full-wave inversion within this frequency range, with a characterization depth varying from 0.5 to 1 m, depending on the soil properties. Yet, the permittivity also affects the results and should be accounted for in the inversion strategy. Field measurements were performed using low-frequency drone-borne radar with a 5-m half-wave dipole antenna, and electromagnetic induction (EMI) measurements with different depth sensitivities were conducted for comparison. Kriging interpolation was used to get maps from measurement points. The soil conductivity maps obtained by the proposed GPR and EMI are compliant in terms of absolute values and spatial patterns. This study demonstrated the capacity of low-frequency drone-borne GPR for fast, field-scale soil electrical conductivity mapping. [ABSTRACT FROM AUTHOR]

Drone-borne ground-penetrating radar (GPR) has proven to be efficient and promising for high-resolution topsoil moisture mapping at the field scale. However, in practice, the radar incident angle may change as a function of the terrain slope and unstable flying conditions. In this respect, we analyzed the effect of radar incident angle on full-wave inversion for soil permittivity characterization. Specifically, using the frequency-domain radar equation, antenna characteristic functions with different incident angles were determined for a horn antenna as an example. Then, numerical analyses and field measurements were performed to quantify the errors on the retrieved permittivity resulting from incident angle drifts. In agreement with the radar equation concept, we observed that the angle- and frequency-dependent global transmission functions of the antenna correspond well to its radiation pattern. Numerical and field results show that errors in the permittivity estimation can be very significant as a function of the antenna radiation pattern and the incident angle. Nevertheless, if the incident angle is known and the characteristic antenna functions are determined as a function of the incident angle, full-wave modeling and inversion remain accurate, provided that the signal-to-noise ratio remains sufficient. [ABSTRACT FROM AUTHOR]

We present a closed-form equation for intrinsic modeling of near-field electromagnetic induction (EMI) antennas for planar layered media characterization. Resorting to a decomposition of the backscattered EM field into elementary distributions over the antenna aperture, the EMI transmitting and receiving antennas are modeled using infinitesimal magnetic dipoles and field points, and characteristic frequency-dependent global reflection and transmission coefficients. Low-frequency propagation of the EM fields in the medium is described using 3-D planar layered media Green's functions. We performed measurements with a loop antenna situated at different heights, ranging from near-field to far-field conditions, above water of known electrical conductivity to determine its intrinsic properties, and a range of salinity conditions was applied to subsequently validate the proposed model. The EMI system was set up using a vector network analyzer equipped with a prototype EMI antenna specifically designed for this application. The model showed good accuracy for reproducing the observed data, and model inversion provided good estimates of the medium electrical conductivity. Yet, insensitivity of the EMI signal to water electrical conductivity was encountered for low salinity due to the presence of a copper sheet as the bottom boundary condition of the experimental setup. Moreover, the efficiency of the antenna decreased rapidly as antenna height above water surface increases, leading to increasing discrepancies between estimated and measured water electrical conductivity values as the antenna moves away from the water surface. Although some technical improvements are still needed, the proposed approach is promising for quantitative estimation of soil electrical conductivity from EMI data. [ABSTRACT FROM PUBLISHER]

A new near-field radar modeling approach for wave propagation in planar layered media is presented. The radar antennas are intrinsically modeled using an equivalent set of infinitesimal electric dipoles and characteristic, frequency-dependent, global reflection, and transmission coefficients. These coefficients determine through a plane wave decomposition wave propagation between the radar reference plane, point sources, and field points. The interactions between the antenna and layered medium are thereby inherently accounted for. The fields are calculated using 3-D Green's functions. We validated the model using an ultrawideband frequency-domain radar with a transmitting and receiving Vivaldi antenna operating in the range 0.8-4 GHz. The antenna characteristic coefficients are obtained from near- and far-field measurements over a copper plane. The proposed model provides unprecedented accuracy for describing near-field radar measurements collected over a water layer, the frequency-dependent electrical properties of which were described using the Debye model. Layer thicknesses could be retrieved through full-wave inversion. The proposed approach demonstrated great promise for nondestructive testing of planar materials and digital soil mapping using ground-penetrating radar. [ABSTRACT FROM AUTHOR]

Ground-penetrating radars operating at a stand-off distance from the probed domain are attractive diagnostic tools in several contexts. Despite notable advantages offered in terms of noninvasiveness and flexibility, the usage of stand-off configurations makes the imaging process more difficult. Indeed, an accurate reconstruction of the subsurface/hidden scene requires accounting for the presence of the air–soil interface into the scattering model. This article deals with a comparison of three different linear models for describing the scattering phenomenon under a stand-off configuration, where the presence of the layered medium is explicitly accounted for in the scattering model. The first model is rigorous and based on the spectral domain representation of the incident field and Green’s function for the half-space scenario. The second model exploits a ray representation of the wave propagation from the radar to the target and vice versa by accounting for the refraction phenomenon at the interface. The third model is an approximated one, where the layered nature of the scenario is taken into account by means of an equivalent wavenumber, which allows a straightforward evaluation of the incident field and Green’s function. The resolution performance of each inversion approach is assessed by means of a numerical analysis. Moreover, data referred to an experiment carried out in controlled conditions are processed and reconstructions are provided to support the comparative analysis among the three imaging approaches. [ABSTRACT FROM AUTHOR]

Multifrequency and multioffset ground-penetrating radar data acquisition modes are used to maximize the information content and parameter retrieval capabilities. However, they also increase the computational cost dedicated to the inversion procedure. In this paper, the impact of the number of frequencies and the multistatic configurations on the information retrieval capabilities is investigated through the response surface topographies of the objective functions. We resort to a full-wave-inversion procedure and a recently developed electromagnetic model which takes advantage of a closed-form solution of Maxwell's equations to describe the antenna–medium system. We show with numerical and laboratory experiments the possibility of reducing the number of frequencies from several hundreds to one or several tens of components without affecting the information retrieval capabilities. We also show through several scenarios that the presence of a perfect electrical conductor increases the number of frequencies required to ensure an acceptable retrieval of the subsurface properties whereas the conductivity of the first layer and the relative permittivity of the second layer do not affect it. The results highlight that information content analyses are important in order to study and optimize data acquisition and inversion procedures, and thereby the computation time. [ABSTRACT FROM AUTHOR]

Litter is acknowledged to affect microwave remote sensing data over forests, and accurate quantification of litter radiative properties is essential for proper retrieval of forest soil water content and forest canopy biophysical characteristics from satellite or airborne sensors. Furthermore, detailed characterization of forest soil organic horizons is of paramount importance for ecological studies and for carbon cycle and global changerelated studies. Ultrawideband ground-penetrating radar (GPR) data acquired above beech forest litter were used to analyze the effects of litter on the backscattering and investigate the ability of GPR to quantitatively characterize the properties of the forest floor. Radar measurements were performed over different litter configurations (i.e., single layers of recently fallen litter (OL) or of fragmented litter (OF), and combined OL-OF layers) and thicknesses. Radar data were processed through full-wave inversion. Over the considered frequency range (0.8-4.0 GHz), attenuation of the radar signal increased as both frequency and litter thickness increase, due to scattering and dielectric losses. These effects were accounted for in the radar model through frequency dependence of litter apparent electrical conductivity. Good agreement (RMSE = 9.5 × 10-3 m) was observed between estimated and measured litter thicknesses. Radar data inversions provided reliable estimates of litter electromagnetic properties, with relative dielectric permittivity values around 1.2 and 3.9 for OL and OF litters, respectively, and corresponding values of 0.032 and 0.080 Sm-1 for litter apparent electrical conductivity at 4.0 GHz. These results show great promise for the use of GPR for noninvasive characterization of forest litter. [ABSTRACT FROM PUBLISHER]

We present validation results of a new ground-penetrating radar (GPR) near-field model for determining the electrical properties and correlated water content of a sand using both frequency- and time-domain radars. The radar antennas are intrinsically characterized using an equivalent set of infinitesimal source/field points and characteristic functions of antennas, which were determined using measurements with the antenna at different distances from a copper plane. The antenna radiation was modeled using six source and field points, which was found to be a good compromise between high modeling accuracy and computing efficiency. We validated our model by inverting GPR data to predict the water content of a sand layer subject to seven levels of saturation. A soil dielectric mixing model was integrated into the full-wave GPR inverse modeling to directly estimate the water content and to account for the frequency dependence of the electrical properties. Although the quality of the fit slightly decreased as the antenna approached the sand surface, the results showed a close agreement between measured and modeled data, resulting in accurate estimation of the water content. The average errors of all water content estimates were 0.012 cm3/cm3 for the frequency domain and 0.016 cm3/cm3 for the time-domain GPR. However, the accuracy reduced when the sand became wet. By performing numerical simulations, we found that it is due to the vertical heterogeneity of soil moisture under the effect of the hydrostatic pressure. We also showed that the GPR inversion with the multilayered soil model could account for this heterogeneity and improved the agreement between the modeled and measured GPR data as well as the accuracy of soil moisture estimation. As for the frequency dependence of the electrical properties, in the frequency ranges of both GPR systems, while the dielectric permittivity was approximately constant, the apparent conductivity exponentially increased with increasing frequency. The success of the calibration and validation in laboratory conditions demonstrates a great potential for practical applications of the radar model, notably for the digital soil mapping and nondestructive testing of materials. [ABSTRACT FROM PUBLISHER]

Horizontally stratified media are commonly used to represent naturally occurring and man-made structures, such as soils, roads, and pavements, when probed by ground-penetrating radar (GPR). Electromagnetic (EM) wave scattering from such multilayered media is dependent on the roughness of the interfaces. In this paper, we developed a closed-form asymptotic EM model considering random rough layers based on the scalar Kirchhoff-tangent plane approximation (SKA) model that we combined with planar multilayered media Green’s functions. In order to validate our extended SKA model, we conducted simulations using a numerical EM solver based on the finite-difference time-domain (FDTD) method. We modeled a medium with three layers—a base layer of perfect electric conductor (PEC) overlaid by two layers of different materials with rough interfaces. The reflections at the first and at the second interface were both well reproduced by the SKA model for each roughness condition. For the reflection at the PEC surface, the extended SKA model slightly overestimated the reflection, and this overestimation increased with the roughness amplitude. Good agreement was also obtained between the FDTD simulation input values and the inverted root mean square (rms) height estimates of the top interface, while the inverted rms heights of the second interface were slightly overestimated. The accuracy and the performances of our asymptotic forward model demonstrate the promising perspectives for simulating rough multilayered media and, hence, for the full waveform inversion of GPR data to noninvasively characterize soils and materials. [ABSTRACT FROM AUTHOR]

In this paper, we propose an automated detection algorithm for well- and ill-shaped ground-penetrating radar reflection hyperbolas for complex media, calibrated with human recognition principles. The algorithm detects the apex of the hyperbolas by fitting an analytical function of a hyperbola to the profile edge dots detected with a Canny filter. The existence of a hyperbola is determined using a set of carefully chosen criteria calibrated in order to fit the ambiguity zone for the human brain. The inherent misshapedness of field hyperbolas is further considered by defining a buffer zone around the theoretical hyperbola. First, the method was tested in the laboratory over tree roots and PVC pipes and on field images over tree root systems. Both time- and frequency-domain radars were used on-ground. After around 1–3 min of computation time for 10 000 edge dots in a MATLAB environment (single 1.96-GHz processor), the results showed rates of false alarm and nondetection of maximum 20% and 28%, respectively. In comparison with the semiautomated hyperbola detection provided by a commercial software, these rates were lower. Second, we conducted a sensitivity analysis to estimate the validity of the fitting of a hyperbola equation neglecting the object radius. The fitting was close, but the derivation of the relative permittivity from the analytical equation neglecting the radius led to high errors. In conclusion, owing to the low computational time and its good performances, the proposed algorithm is suitable for complex environments. [ABSTRACT FROM AUTHOR]

The characterization of the subsurface can be performed by full-waveform inversion of electromagnetic data. The modeling process relies on the ability to retrieve the scattered field Green's function from the measured data using sets of antenna characteristic global reflection and transmission coefficients. Crucial for a successful implementation of this technique is the understanding of uncertainties involved in the acquisition of the antenna calibration and survey measurements, and how these propagate in the parameter estimation results. We find that averaging a number of possible Green's functions obtained from one measurement with several antenna characteristic coefficients sets works remarkably well in reducing the uncertainties. The accuracy of the inversions improves using characteristic coefficients acquired as close as possible to the measurement conditions. A clear relation between dynamic range and system resolution is highlighted, based on the number of effective bits contained in the data. [ABSTRACT FROM PUBLISHER]

A combined study of radar profiles and thin section analysis supported by modelling of synthetic radar traces reveals that ground-penetrating radar (GPR) reflections generated in diagenetically altered sandstones cannot always be interpreted unequivocally. This is illustrated in the Luxembourg Sandstone Formation, which has been altered diagenetically by selective carbonate cementation and fracturing. Cemented lenses and concretions developed along the bedding planes, especially at places with high primary carbonate content. Cementation resulted in the alternation of cemented carbonate-rich sandy layers (thickness 30–50 cm and variable length) and uncemented carbonate-poor sandstone layers. The ability of GPR to detect the geometry of these lenses and vertical fractures with centimetre apertures has been tested at several antenna frequencies (100, 200, 250 and 500 MHz). Relative dielectric permittivity calculations were carried out to assess variations of this electric property within the cemented and uncemented layers as a function of porosity, calcite and water content in the pores. Two-dimensional full waveform modelling was also carried out to study the effect of conductivity in the sandstones and the effect of interlayer clay seams. At the penetration depth of the radar (7 m with 250 MHz), cemented lenses and concretions could only be detected with GPR when the porosity contrast was sufficiently high, which is not always the case. This conclusion is supported by the modelling. The data also proved the ability of radar to detect large open vertical fractures along which sandstones are weathered. The study has implications for investigations which will use GPR to detect three-dimensional distribution of diagenetic pore filling precipitates as well as secondary porosity development along fractures. [ABSTRACT FROM AUTHOR]