A large underground inhomogeneity, such as a salt dome or cavity, is known to disturb the subsurface temperature field. Such anomalies appear in many geophysical surveys. Detection and knowledge of the magnitude of these disturbances is the objective of both near surface and deep borehole temperature surveys aimed at delineating the inhomogeneities. It also impacts surface temperature history analysis which reconstructs past climate change in an effort to study the recent global warming. This work is aimed at quantifying these effects by solving a problem of heat conduction in Earth's subsurface in the presence of a spherical inhomogeneity. Both the steady state temperature field pertaining to the constant geothermal gradient and the time dependent field caused by a surface jump in temperature are solved. A solution is derived for both cases as an infinite series of spherical harmonics and Bessel functions (in the Laplace domain) for the steady and unsteady problems, respectively. It is shown that an accurate solution can be achieved by a small number of terms. The results are illustrated and analysed for a given accuracy and for a few values of the governing parameters. The general solution can be simplified considerably for asymptotic values of the parameters. Comparison with the exact solution shows that these approximations are accurate for a wide range of parameter values. Some examples of applying the solution to the geophysical methods stated above are discussed. In the case of surface temperature surveys, an example showing good agreement between theoretical and measured heat flux above a salt diapir is given. Though the sphere is an idealized shape, the simplicity of the solution makes possible a general analysis toward gaining a better understanding of the process. Furthermore, it can be employed for preliminary assessment of the impact of a body and may serve as a benchmark for numerical solutions.
- Composition of the continental crust
- Downhole methods
- Heat flow