Deformation induced by distributions of single forces in a layered half-space: EFGRN/EFCMP

In the present paper we introduce a numerical model for the representation of displacement, strain and stress due to single forces embedded in a layered elastic half-space. The code EFGRN/EFCMP (Elastic Forces GReeN functions/Elastic Forces CoMPutation) is able to represent the mechanical effects due to preassigned distributions of single forces. Even if internal deformation sources can be described by distributions of equivalent body forces with vanishing resultant and moment, single forces are employed in geophysics to represent hydraulic and/or lithostatic loads, effects of internal density anomalies, and even some kind of seismic events. A distribution of single forces is also used to describe the effects of an inelastic inclusion located inside an elastic medium. In fact, the recent literature shows that poro-elastic and thermo-elastic inclusions can be represented using single forces distributed on their boundaries. EFGRN/EFCMP shares the benefits of rapid and semi-analytical calculation offered by the parent code, EDGRN/EDCMP, which is instead suitable for the representation of extended dislocation sources, as seismic faults. The present code also provides an option for computing the effects of a distribution of single forces embedded in a homogeneous half-space, by using the analytical solutions of Mindlin. Accordingly, EFGRN/EFCMP can be a valid support both for the representation of forward models of deformation sources and for the procedures of inversion of geodetic data in a layered medium. We show some applications of the code and we provide several scripts in MATLAB language which help the user to quickly start using EFGRN/EFCMP.

Automated summary

This paper introduces a numerical model for representing displacement, strain, and stress due to single forces in a layered elastic half-space. The model can support forward modeling of deformation sources and inversion of geodetic data in a layered medium.