A shape derivative algorithm for reconstructing elastic dislocations in geophysics

We consider the inverse problem of determining an elastic dislocation that models a seismic fault in the quasi-static regime of aseismic, creeping faults, from displacement measurements made at the surface of Earth. We derive both a distributed as well as a boundary shape derivative that encodes the change in a misfit functional between the measured and the computed surface displacement under infinitesimal movements of the dislocation and infinitesimal changes in the slip vector, which gives the displacement jump across the dislocation. We employ the shape derivative in an iterative reconstruction algorithm. We present some numerical test of the reconstruction algorithm in a simplified 2D setting.