eprintid: 1239 rev_number: 19 eprint_status: archive userid: 2 dir: disk0/00/00/12/39 datestamp: 2015-02-17 11:28:05 lastmod: 2018-03-28 09:58:55 status_changed: 2015-04-27 12:10:31 type: article metadata_visibility: show creators_name: Martinez-Garzon, Patricia creators_name: Kwiatek, Grzegorz creators_name: Ickrath, Michèle creators_name: Bohnhoff, Marco creators_id: creators_id: kwiatek@gfz-potsdam.de creators_id: creators_id: bohnhoff@gfz-potsdam.de corp_creators: Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences, Section 3.2: Geomechanics and Rheology, Telegrafenberg, Potsdam, Germany patricia@gfz-potsdam.de title: MSATSI: A MATLAB Package for Stress Inversion Combining Solid Classic Methodology, a New Simplified User-Handling, and a Visualization Tool ispublished: pub subjects: MP6_1 divisions: IP13 full_text_status: none abstract: The estimation of the stress‐field orientation from focal mechanisms of earthquakes is a relevant tool to understand crustal mechanics and the physics of earthquakes. In global seismology, Formal Stress Inversion (FSI) is a well‐established technique to study tectonic processes associated with the occurrence of large earthquakes (e.g., Hardebeck and Michael, 2004; Yoshida et al., 2012). Information on the stress‐field orientation is also relevant for the exploitation of hydrocarbon and geothermal reservoirs. Knowledge of the orientation of the maximum horizontal stress (σHmax) is crucial for reservoir development, such as drilling and leakage risk assessment (Terakawa et al., 2012; Martínez‐Garzón et al., 2013). Additionally, the FSI technique may be useful for understanding the physics of rupture processes at a microscale in the frame of rock deformation experiments in the laboratory. In seismological practice, the stress tensor is obtained either from inverting earthquake focal mechanisms or directly from first‐motion polarities. Most of the developed FSI methods share two main assumptions: The stress field is homogeneous within the considered rock volume. The slip of the fault is parallel to the direction of the tangential traction (Wallace, 1951; Bott, 1959). Estimation of stress‐field orientation is a nonlinear inverse problem that can be solved either directly or linearized. When solving the nonlinear inverse problem, the best‐fitting stress tensor to the group of focal mechanisms is obtained using grid‐search algorithms (Gephart and Forsyth, 1984; Gephart, 1990; Arnold and Townend, 2007) or Monte Carlo sampling‐based optimization methods (Angelier, 1984; Xu, 2004). The linearized inversion scheme is solved by a generally less computationally extensive least‐squares technique (Michael, 1984; Hardebeck and Michael, 2006). date: 2014-07 publication: Seismological Research Letters volume: 85 number: 4 publisher: Seismological Society of America pagerange: 896-904 id_number: doi:10.1785/0220130189 refereed: TRUE issn: 0895-0695 official_url: http://dx.doi.org/10.1785/0220130189 access_IS-EPOS: limited software_references: Stress_Inversion owner: Publisher citation: Martinez-Garzon, Patricia and Kwiatek, Grzegorz and Ickrath, Michèle and Bohnhoff, Marco (2014) MSATSI: A MATLAB Package for Stress Inversion Combining Solid Classic Methodology, a New Simplified User-Handling, and a Visualization Tool. Seismological Research Letters, 85 (4). pp. 896-904. DOI: https://doi.org/10.1785/0220130189