Transient electromagnetic calculations using the Gaver-Stehfest inverse Laplace transform method

Calculations for the transient electromagnetic (TEM) method are commonly performed by using a discrete Fourier transform method to invert the appropriate transform of the solution. We derive the Laplace transform of the solution for TEM soundings over an N-layer earth and show how to use the Gaver-Stehfest algorithm to invert it numerically. This is considerably more stable and computationally efficient than inversion using the discrete Fourier transform.