An inverse diffusivity problem for the helium production–diffusion equation

Abstract

Thermochronology is a technique for the extraction of information about the thermal history of rocks. Such information is crucial for determining the depth below the surface at which rocks were located at a given time (Bao G et al 2011 Commun. Comput. Phys. 9 129). Mathematically, extracting the time–temperature path can be formulated as an inverse diffusivity problem for the helium production–diffusion equation which is the underlying process of thermochronology. In this paper, to reconstruct the diffusivity which depends on space only and accounts for the structural information of rocks, a local Hölder conditional stability is obtained by a Carleman estimate. A uniqueness result is also proven for extracting the thermal history, i.e. identifying the time-dependant part of the diffusion coefficient, provided that it is analytical with respect to time. Numerical examples are presented to illustrate the validity and effectiveness of the proposed regularization scheme.