A Brief Introduction to Ice-Penetrating Radar
lindzey.github.io
lindzey.github.io
I wonder how big the reflection coefficient is for the air/ice interface. It seems like it would be huge. So maybe ground based techniques offer better coupling at the cost of not being able to survey as much area?
We transmit 8kW, and the air/ice surface reflection coefficient is ~0.08 (~-11dB). Flying at ~600m above the surface, spreading loss actually contributes more to signal attenuation (1/(2*h)^2 ~= -62dB).
Our instrument is optimized for seeing through the entire ice sheet, mapping deep layers and the bed. Other (also airborne) instruments operate at higher frequencies, trading higher resolution for less penetration. I'm not super familiar with groups using ground-based ice-penetrating radar, but one big tradeoff is $$$. The airplane is hugely expensive to operate, whereas ground-based just needs a snowmobile.
This is trickier than normal imaging, because if you ignore the multiple reflections you're basically just inverting a matrix (a fourier transform in the case of single line imaging). With multipath the measurement becomes nonlinear in your scene and inversion isn't as trivial.
Unfortunately, this only improves the resolution along the flight track; features parallel to the flight path and offset to the side are the hardest to filter out.
There has been some fun work along the lines you describe that uses maps of the surface shape (generated from camera imagery or scanning laser data) to discriminate which apparently-subsurface echoes are most likely due to the surface topography. So far as I know, this hasn't been automated - it's more an aid to human interpretation (we have an army of undergrads that "picks" the most likely bed location).
Do you have papers describing y'all's processing? I'd love to learn more!
http://euler.msi.ucsb.edu/papers/2014_emery_apm_from_soo.pdf
It has some of the basic processing techniques.
We do use a chirped signal in order to improve vertical resolution - convolving the outgoing pulse with the returned trace gives us the equivalent resolution of a ~80ns pulse, but with more power.