Three-dimensional finite-element modelling of coseismic Coulomb stress changes on intra-continental dip-slip faults

• Systematic study about coseismic Coulomb stress change on dip-slip fault planes
• 3D finite-element models with normal fault and thrust fault arrays
• Fault dip, dip direction and position influence the distribution of stress changes.
• Coseismic slip and friction coefficient affect the magnitude of the stress changes.
• Results reveal correlation between the coseismic stress changes and strain field.

Investigating fault interaction plays a crucial role in seismic hazard assessment. The calculation of Coulomb stress changes allows quantifying the stress changes on so-called receiver faults in the surrounding of the fault that experienced the earthquake. A positive stress change implies that the earthquake brought the receiver fault closer to failure while a negative value indicates a delay of the next earthquake. So far, most studies focussed on stress changes for particular faults and earthquakes. Here we present a general analysis of the Coulomb stress changes on intra-continental dip-slip faults using finite-element models with normal and thrust faults arrays, respectively. Our models allow calculating coseismic (“static”) stress changes on pre-defined fault planes, whose dip and position can be varied. Gravity and ongoing regional deformation (i.e. shortening or extension) are included. The results for thrust and normal faults show that synthetic receiver faults located in the hanging wall and footwall of the source fault exhibit a symmetric stress distribution, with large areas of negative and small areas of positive Coulomb stress changes. In contrast, faults positioned in along-strike prolongation of the source fault and outside of its immediate hanging wall and footwall undergo mostly positive stress changes. The stress changes are largest at the fault tip that is closer to the source fault. Our results show that the stress change distribution depends on the fault dip while the magnitude depends on the friction coefficient and the amount of coseismic slip. The Coulomb stress changes can be explained by the spatial distribution of the coseismic strain, which shows domains of horizontal extension and shortening that alternate both at the surface and with depth. Our models allow identifying the general patterns of Coulomb stress changes on dip-slip faults, which are often concealed by the peculiarity of the specific fault or earthquake in nature.