skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Joint inversion of full-waveform ground-penetrating radar and electrical resistivity data — Part 2: Enhancing low frequencies with the envelope transform and cross gradients
Title: Joint inversion of full-waveform ground-penetrating radar and electrical resistivity data — Part 2: Enhancing low frequencies with the envelope transform and cross gradients
Recovering material properties of the subsurface using ground-penetrating radar (GPR) data of finite bandwidth with missing low frequencies and in the presence of strong attenuation is a challenging problem. We have adopted three nonlinear inverse methods for recovering electrical conductivity and permittivity of the subsurface by joining GPR multioffset and electrical resistivity (ER) data acquired at the surface. All of the methods use ER data to constrain the low spatial frequency of the conductivity solution. The first method uses the envelope of the GPR data to exploit low-frequency content in full-waveform inversion and does not assume structural similarities of the material properties. The second method uses cross gradients to manage weak amplitudes in the GPR data by assuming structural similarities between permittivity and conductivity. The third method uses the envelope of the GPR data and the cross gradient of the model parameters. By joining ER and GPR data, exploiting low-frequency content in the GPR data, and assuming structural similarities between the electrical permittivity and conductivity, we are able to recover subsurface parameters in regions where the GPR data have a signal-to-noise ratio close to one.  more » « less
Award ID(s):
1720472
PAR ID:
10309733
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
GEOPHYSICS
Volume:
85
Issue:
6
ISSN:
0016-8033
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, GEOPHYSICS)
    We have developed an algorithm for joint inversion of full-waveform ground-penetrating radar (GPR) and electrical resistivity (ER) data. The GPR data are sensitive to electrical permittivity through reflectivity and velocity, and electrical conductivity through reflectivity and attenuation. The ER data are directly sensitive to the electrical conductivity. The two types of data are inherently linked through Maxwell’s equations, and we jointly invert them. Our results show that the two types of data work cooperatively to effectively regularize each other while honoring the physics of the geophysical methods. We first compute sensitivity updates separately for the GPR and ER data using the adjoint method, and then we sum these updates to account for both types of sensitivities. The sensitivities are added with the paradigm of letting both data types always contribute to our inversion in proportion to how well their respective objective functions are being resolved in each iteration. Our algorithm makes no assumptions of the subsurface geometry nor the structural similarities between the parameters with the caveat of needing a good initial model. We find that our joint inversion outperforms the GPR and ER separate inversions, and we determine that GPR effectively supports ER in regions of low conductivity, whereas ER supports GPR in regions with strong attenuation. 
    more » « less
  2. ; ; (, GEOPHYSICS)
    We exploit the different but complementary data sensitivities of ground-penetrating radar (GPR) and electrical resistivity (ER) by applying a multiphysics, multiparameter, simultaneous 2.5D joint inversion without invoking petrophysical relationships. Our method joins full-waveform inversion (FWI) GPR with adjoint derived ER sensitivities on the same computational domain. We incorporate a stable source estimation routine into the FWI-GPR. We apply our method in a controlled alluvial aquifer using only surface-acquired data. The site exhibits a shallow groundwater boundary and unconsolidated heterogeneous alluvial deposits. We compare our recovered parameters to individual FWI-GPR and ER results, and we compare them to log measurements of capacitive conductivity and neutron-derived porosity. Our joint inversion provides a more representative depiction of subsurface structures because it incorporates multiple intrinsic parameters, and it is therefore superior to an interpretation based on log data, FWI-GPR, or ER alone. 
    more » « less
  3. ; ; (, GEOPHYSICS)
    We have developed a memory and operation-count efficient 2.5D inversion algorithm of electrical resistivity (ER) data that can handle fine discretization domains imposed by other geophysical (e.g, ground penetrating radar or seismic) data. Due to numerical stability criteria and available computational memory, joint inversion of different types of geophysical data can impose different grid discretization constraints on the model parameters. Our algorithm enables the ER data sensitivities to be directly joined with other geophysical data without the need of interpolating or coarsening the discretization. We have used the adjoint method directly in the discretized Maxwell’s steady state equation to compute the data sensitivity to the conductivity. In doing so, we make no finite-difference approximation on the Jacobian of the data and avoid the need to store large and dense matrices. Rather, we exploit matrix-vector multiplication of sparse matrices and find successful convergence using gradient descent for our inversion routine without having to resort to the Hessian of the objective function. By assuming a 2.5D subsurface, we are able to linearly reduce memory requirements when compared to a 3D gradient descent inversion, and by a power of two when compared to storing a 2D Hessian. Moreover, our method linearly outperforms operation counts when compared with 3D Gauss-Newton conjugate-gradient schemes, which scales cubically in our favor with respect to the thickness of the 3D domain. We physically appraise the domain of the recovered conductivity using a cutoff of the electric current density present in our survey. We evaluate two case studies to assess the validity of our algorithm. First, on a 2.5D synthetic example, and then on field data acquired in a controlled alluvial aquifer, where we were able to match the recovered conductivity to borehole observations. 
    more » « less
  4. ; (, The Cryosphere)
    null (Ed.)
    Abstract. We have examined a general expression giving the specularreflection coefficient for a radar wave approaching a reflecting interfacewith normal incidence. The reflecting interface separates two homogeneousisotropic media, the properties of which are fully described by three scalarquantities: dielectric permittivity, magnetic permeability, and electricalconductivity. The derived relationship indicates that electricalconductivity should not be neglected a priori in glaciological investigations ofsubglacial materials and in ground-penetrating radar (GPR) studies of saturated sediments and bedrock,even at the high end of typical linear radar frequencies used in suchinvestigations (e.g., 100–400 MHz). Our own experience in resistivitysurveying in Antarctica, combined with a literature review, suggests that awide range of geologic materials can have electrical conductivity that ishigh enough to significantly impact the value of radar reflectivity.Furthermore, we have given two examples of prior studies in which inclusionof electrical conductivity in calculation of the radar bed reflectivity mayprovide an explanation for results that may be considered surprising if theimpact of electrical conductivity on radar reflection is neglected. Thecommonly made assumption that only dielectric permittivity of the two medianeeds to be considered in interpretation of radar reflectivity can lead toerroneous conclusions. 
    more » « less
  5. (, Environmental Data Initiative)
    30-minute Soil Moisture data at 10 cm at CSIS Block-13. Collection of soil volumetric water content data at Jornada LTER Cross-scale Interactions Study (CSIS) blocks, New Mexico, supports the environmental monitoring objectives of the Jornada LTER monitoring program that look at plant-soil water dynamics. Soil volumetric water content is measured at 10 cm depth at 3 locations at each of the 4 plots per block. The locations in each plot are under bare soil, under perennial grass, and under shrub. Volumetric water content, bulk electrical conductivity, soil temperature, and bulk dielectric permittivity are recorded every 30 minutes at an automated meteorological station installed at Jornada LTER CSIS Block-13 site. This is an ONGOING dataset. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email