CoulombAnalysis: Program for Flexibly Calculating Coulomb Stress Changes on Receiver Faults with Any Orientations and Non-vertical Profiles with Any Trends

Over the past more than thirty years, the Coulomb failure model has been employed to compute Coulomb stress changes imparted by mainshocks, magmatic motions or slow slip events for the sake of deciphering potential physical causal linkages among such events. For the computation of Coulomb stress changes, two special programs of Coulomb3.4 and PSGRN-PSCMP have been popularly used. Yet they cannot be directly used to work out Coulomb stress changes on various profiles and on diverse optimally oriented failure planes (OOPs). We here present a program that can be utilized to wholesale and automatically calculate Coulomb stress changes on nodal planes of focal mechanisms, profiles, OOPs and major faults. In addition, the explicit mathematical framework underlying the program is offered for the first time to our knowledge. Both the program and the mathematical framework associated with the Coulomb failure model should be of great benefit to deepening insights on earthquake stress triggering and fault-network interaction in the near future.