@article{DurandBentzKwiateketal.2020,
  author    = {Durand, Virginie and Bentz, Stephan and Kwiatek, Grzegorz and Dresen, Georg and Wollin, Christopher and Heidbach, Oliver and Martinez-Garzon, Patricia and Cotton, Fabrice Pierre and Nurlu, Murat and Bohnhoff, Marco},
  title     = {A two-scale preparation phase preceded an M-w 5.8 earthquake in the sea of marmara offshore Istanbul, Turkey},
  series = {Seismological research letters},
  volume    = {91},
  journal   = {Seismological research letters},
  number    = {6},
  address   = {Boulder},
  issn      = {0895-0695},
  doi       = {10.1785/0220200110},
  pages     = {3139 -- 3147},
  year      = {2020},
  abstract  = {We analyze the spatiotemporal evolution of seismicity during a sequence of moderate (an M-w 4.7 foreshock and M-w 5.8 mainshock) earthquakes occurring in September 2019 at the transition between a creeping and a locked segment of the North Anatolian fault in the central Sea of Marmara, northwest Turkey. To investigate in detail the seismicity evolution, we apply a matched-filter technique to continuous waveforms, thus reducing the magnitude threshold for detection. Sequences of foreshocks preceding the two largest events are clearly seen, exhibiting two different behaviors: a long-term activation of the seismicity along the entire fault segment and a short-term concentration around the epicenters of the large events. We suggest a two-scale preparation phase, with aseismic slip preparing the mainshock final rupture a few days before, and a cascade mechanism leading to the nucleation of the mainshock. Thus, our study shows a combination of seismic and aseismic slip during the foreshock sequence changing the strength of the fault, bringing it closer to failure.},
  language  = {en}
}
@article{ZiebarthvonSpechtHeidbachetal.2020,
  author    = {Ziebarth, Malte J. and von Specht, Sebastian and Heidbach, Oliver and Cotton, Fabrice Pierre and Anderson, John G.},
  title     = {Applying conservation of energy to estimate earthquake frequencies from strain rates and stresses},
  series = {Journal of geophysical research : Solid earth},
  volume    = {125},
  journal   = {Journal of geophysical research : Solid earth},
  number    = {8},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1029/2020JB020186},
  pages     = {25},
  year      = {2020},
  abstract  = {Estimating earthquake occurrence rates from the accumulation rate of seismic moment is an established tool of seismic hazard analysis. We propose an alternative, fault-agnostic approach based on the conservation of energy: the Energy-Conserving Seismicity Framework (ENCOS). Working in energy space has the advantage that the radiated energy is a better predictor of the damage potential of earthquake waves than the seismic moment release. In a region, ENCOS balances the stationary power available to cause earthquakes with the long-term seismic energy release represented by the energy-frequency distribution's first moment. Accumulation and release are connected through the average seismic efficiency, by which we mean the fraction of released energy that is converted into seismic waves. Besides measuring earthquakes in energy, ENCOS differs from moment balance essentially in that the energy accumulation rate depends on the total stress in addition to the strain rate tensor. To validate ENCOS, we exemplarily model the energy-frequency distribution around Southern California. We estimate the energy accumulation rate due to tectonic loading assuming poroelasticity and hydrostasis. Using data from the World Stress Map and assuming the frictional limit to estimate the stress tensor, we obtain a power of 0.8 GW. The uncertainty range, 0.3-2.0GW, originates mainly from the thickness of the seismogenic crust, the friction coefficient on preexisting faults, and models of Global Positioning System (GPS) derived strain rates. Based on a Gutenberg-Richter magnitude-frequency distribution, this power can be distributed over a range of energies consistent with historical earthquake rates and reasonable bounds on the seismic efficiency.},
  language  = {en}
}
@article{PenaHeidbachMorenoetal.2019,
  author    = {Pena, Carlos and Heidbach, Oliver and Moreno, Marcos and Bedford, Jonathan and Ziegler, Moritz 0. and Tassara, Andres Ollero and Oncken, Onno},
  title     = {Role of Lower Crust in the Postseismic Deformation of the 2010 Maule Earthquake: Insights from a Model with Power-Law Rheology},
  series = {Pure and applied geophysics},
  volume    = {176},
  journal   = {Pure and applied geophysics},
  number    = {9},
  publisher = {Springer},
  address   = {Basel},
  issn      = {0033-4553},
  doi       = {10.1007/s00024-018-02090-3},
  pages     = {3913 -- 3928},
  year      = {2019},
  abstract  = {The surface deformation associated with the 2010 M-w 8.8 Maule earthquake in Chile was recorded in great detail before, during and after the event. The high data quality of the continuous GPS (cGPS) observations has facilitated a number of studies that model the postseismic deformation signal with a combination of relocking, afterslip and viscoelastic relaxation using linear rheology for the upper mantle. Here, we investigate the impact of using linear Maxwell or power-law rheology with a 2D geomechanical-numerical model to better understand the relative importance of the different processes that control the postseismic deformation signal. Our model results reveal that, in particular, the modeled cumulative vertical postseismic deformation pattern in the near field (< 300 km from the trench) is very sensitive to the location of maximum afterslip and choice of rheology. In the model with power-law rheology, the afterslip maximum is located at 20-35 km rather than > 50 km depth as suggested in previous studies. The explanation for this difference is that in the model with power-law rheology the relaxation of coseismically imposed differential stresses occurs mainly in the lower crust. However, even though the model with power-law rheology probably has more potential to explain the vertical postseismic signal in the near field, the uncertainty of the applied temperature field is substantial, and this needs further investigations and improvements.},
  language  = {en}
}
@misc{GholamrezaieScheckWenderothBottetal.2019,
  author    = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Bott, Judith and Heidbach, Oliver and Strecker, Manfred},
  title     = {3-D crustal density model of the Sea of Marmara},
  series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe},
  number    = {737},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-43466},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-434661},
  pages     = {785 -- 807},
  year      = {2019},
  abstract  = {Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible "end-member" solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.},
  language  = {en}
}
@article{GholamrezaieScheckWenderothBottetal.2019,
  author    = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Bott, Judith and Heidbach, Oliver and Strecker, Manfred},
  title     = {3-D crustal density model of the Sea of Marmara},
  series = {Solid Earth},
  volume    = {10},
  journal   = {Solid Earth},
  publisher = {Copernicus Publ.},
  address   = {G{\"o}ttingen},
  issn      = {1869-9510},
  doi       = {10.5194/se-10-785-2019},
  pages     = {785 -- 807},
  year      = {2019},
  abstract  = {Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible "end-member" solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.},
  language  = {en}
}
@misc{HeidbachRajabiCuietal.2018,
  author    = {Heidbach, Oliver and Rajabi, Mojtaba and Cui, Xiaofeng and Fuchs, Karl and Mueller, Birgit and Reinecker, John and Reiter, Karsten and Tingay, Mark and Wenzel, Friedemann and Xie, Furen and Ziegler, Moritz O. and Zoback, Mary-Lou and Zoback, Mark},
  title     = {The World Stress Map database release 2016},
  series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  volume    = {744},
  journal   = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0040-1951},
  doi       = {10.1016/j.tecto.2018.07.007},
  pages     = {484 -- 498},
  year      = {2018},
  abstract  = {Knowledge of the present-day crustal in-situ stress field is a key for the understanding of geodynamic processes such as global plate tectonics and earthquakes. It is also essential for the management of geo-reservoirs and underground storage sites for energy and waste. Since 1986, the World Stress Map (WSM) project has systematically compiled the orientation of maximum horizontal stress (S-Hmax). For the 30th anniversary of the project, the WSM database has been updated significantly with 42,870 data records which is double the amount of data in comparison to the database release in 2008. The update focuses on areas with previously sparse data coverage to resolve the stress pattern on different spatial scales. In this paper, we present details of the new WSM database release 2016 and an analysis of global and regional stress pattern. With the higher data density, we can now resolve stress pattern heterogeneities from plate-wide to local scales. In particular, we show two examples of 40 degrees-60 degrees S-Hmax rotations within 70 km. These rotations can be used as proxies to better understand the relative importance of plate boundary forces that control the long wave-length pattern in comparison to regional and local controls of the crustal stress state. In the new WSM project phase IV that started in 2017, we will continue to further refine the information on the S-Hmax orientation and the stress regime. However, we will also focus on the compilation of stress magnitude data as this information is essential for the calibration of geomechanical-numerical models. This enables us to derive a 3-D continuous description of the stress tensor from point-wise and incomplete stress tensor information provided with the WSM database. Such forward models are required for safety aspects of anthropogenic activities in the underground and for a better understanding of tectonic processes such as the earthquake cycle.},
  language  = {en}
}
@article{vonSpechtHeidbachCottonetal.2018,
  author    = {von Specht, Sebastian and Heidbach, Oliver and Cotton, Fabrice Pierre and Zang, Arno},
  title     = {Uncertainty reduction of stress tensor inversion with data-driven catalogue selection},
  series = {Geophysical journal international},
  volume    = {214},
  journal   = {Geophysical journal international},
  number    = {3},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggy240},
  pages     = {2250 -- 2263},
  year      = {2018},
  abstract  = {The selection of earthquake focal mechanisms (FMs) for stress tensor inversion (STI) is commonly done on a spatial basis, that is, hypocentres. However, this selection approach may include data that are undesired, for example, by mixing events that are caused by different stress tensors when for the STI a single stress tensor is assumed. Due to the significant increase of FM data in the past decades, objective data-driven data selection is feasible, allowing more refined FM catalogues that avoid these issues and provide data weights for the STI routines. We present the application of angular classification with expectation-maximization (ACE) as a tool for data selection. ACE identifies clusters of FM without a priori information. The identified clusters can be used for the classification of the style-of-faulting and as weights of the FM data. We demonstrate that ACE effectively selects data that can be associated with a single stress tensor. Two application examples are given for weighted STI from South America. We use the resulting clusters and weights as a priori information for an STI for these regions and show that uncertainties of the stress tensor estimates are reduced significantly.},
  language  = {en}
}
@article{LiMorenoSwittBedfordetal.2017,
  author    = {Li, Shaoyang and Moreno Switt, Marcos and Bedford, Jonathan and Rosenau, Matthias and Heidbach, Oliver and Melnick, Daniel and Oncken, Onno},
  title     = {Postseismic uplift of the Andes following the 2010 Maule earthquake},
  series = {Geophysical research letters},
  volume    = {44},
  journal   = {Geophysical research letters},
  number    = {4},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1002/2016GL071995},
  pages     = {1768 -- 1776},
  year      = {2017},
  abstract  = {Postseismic surface deformation associated with great subduction earthquakes is controlled by asthenosphere rheology, frictional properties of the fault, and structural complexity. Here by modeling GPS displacements in the 6 years following the 2010 M-w 8.8 Maule earthquake in Chile, we investigate the impact of heterogeneous viscosity distribution in the South American subcontinental asthenosphere on the 3-D postseismic deformation pattern. The observed postseismic deformation is characterized by flexure of the South America plate with peak uplift in the Andean mountain range and subsidence in the hinterland. We find that, at the time scale of observation, over 2 orders of magnitude gradual increase in asthenosphere viscosity from the arc area toward the cratonic hinterland is needed to jointly explain horizontal and vertical displacements. Our findings present an efficient method to estimate spatial variations of viscosity, which clearly improves the fitting to the vertical signal of deformation. Lateral changes in asthenosphere viscosity can be correlated with the thermomechanical transition from weak subvolcanic arc mantle to strong subcratonic mantle, thus suggesting a stationary heterogeneous viscosity structure. However, we cannot rule out a transient viscosity structure (e.g., power law rheology) with the short time span of observation.},
  language  = {en}
}
@article{ZieglerHeidbachZangetal.2017,
  author    = {Ziegler, Moritz O. and Heidbach, Oliver and Zang, Arno and Martinez-Garzon, Patricia and Bohnhoff, Marco},
  title     = {Estimation of the differential stress from the stress rotation angle in low permeable rock},
  series = {Geophysical research letters},
  volume    = {44},
  journal   = {Geophysical research letters},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0094-8276},
  doi       = {10.1002/2017GL073598},
  pages     = {6761 -- 6770},
  year      = {2017},
  abstract  = {Rotations of the principal stress axes are observed as a result of fluid injection into reservoirs. We use a generic, fully coupled 3-D thermo-hydro-mechanical model to investigate systematically the dependence of this stress rotation on different reservoir properties and injection scenarios. We find that permeability, injection rate, and initial differential stress are the key factors, while other reservoir properties only play a negligible role. In particular, we find that thermal effects do not significantly contribute to stress rotations. For reservoir types with usual differential stress and reservoir treatment the occurrence of significant stress rotations is limited to reservoirs with a permeability of less than approximately 10(-12)m(2). Higher permeability effectively prevents stress rotations to occur. Thus, according to these general findings, the observed principal stress axes rotation can be used as a proxy of the initial differential stress provided that rock permeability and fluid injection rate are known a priori.},
  language  = {en}
}
@article{RajabiZieglerTingayetal.2016,
  author    = {Rajabi, Mojtaba and Ziegler, Moritz O. and Tingay, Mark and Heidbach, Oliver and Reynolds, Scott},
  title     = {Contemporary tectonic stress pattern of the Taranaki Basin, New Zealand},
  series = {Journal of geophysical research : Solid earth},
  volume    = {121},
  journal   = {Journal of geophysical research : Solid earth},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {2169-9313},
  doi       = {10.1002/2016JB013178},
  pages     = {6053 -- 6070},
  year      = {2016},
  abstract  = {The present-day stress state is a key parameter in numerous geoscientific research fields including geodynamics, seismic hazard assessment, and geomechanics of georeservoirs. The Taranaki Basin of New Zealand is located on the Australian Plate and forms the western boundary of tectonic deformation due to Pacific Plate subduction along the Hikurangi margin. This paper presents the first comprehensive wellbore-derived basin-scale in situ stress analysis in New Zealand. We analyze borehole image and oriented caliper data from 129 petroleum wells in the Taranaki Basin to interpret the shape of boreholes and determine the orientation of maximum horizontal stress (S-Hmax). We combine these data (151 S-Hmax data records) with 40 stress data records derived from individual earthquake focal mechanism solutions, 6 from stress inversions of focal mechanisms, and 1 data record using the average of several focal mechanism solutions. The resulting data set has 198 data records for the Taranaki Basin and suggests a regional S-Hmax orientation of N068 degrees E (22 degrees), which is in agreement with NW-SE extension suggested by geological data. Furthermore, this ENE-WSW average S-Hmax orientation is subparallel to the subduction trench and strike of the subducting slab (N50 degrees E) beneath the central western North Island. Hence, we suggest that the slab geometry and the associated forces due to slab rollback are the key control of crustal stress in the Taranaki Basin. In addition, we find stress perturbations with depth in the vicinity of faults in some of the studied wells, which highlight the impact of local stress sources on the present-day stress rotation.},
  language  = {en}
}