@article{SharmaHainzlZoelleretal.2020, author = {Sharma, Shubham and Hainzl, Sebastian and Z{\"o}ller, Gert and Holschneider, Matthias}, title = {Is Coulomb stress the best choice for aftershock forecasting?}, series = {Journal of geophysical research : Solid earth}, volume = {125}, journal = {Journal of geophysical research : Solid earth}, number = {9}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1029/2020JB019553}, pages = {12}, year = {2020}, abstract = {The Coulomb failure stress (CFS) criterion is the most commonly used method for predicting spatial distributions of aftershocks following large earthquakes. However, large uncertainties are always associated with the calculation of Coulomb stress change. The uncertainties mainly arise due to nonunique slip inversions and unknown receiver faults; especially for the latter, results are highly dependent on the choice of the assumed receiver mechanism. Based on binary tests (aftershocks yes/no), recent studies suggest that alternative stress quantities, a distance-slip probabilistic model as well as deep neural network (DNN) approaches, all are superior to CFS with predefined receiver mechanism. To challenge this conclusion, which might have large implications, we use 289 slip inversions from SRCMOD database to calculate more realistic CFS values for a layered half-space and variable receiver mechanisms. We also analyze the effect of the magnitude cutoff, grid size variation, and aftershock duration to verify the use of receiver operating characteristic (ROC) analysis for the ranking of stress metrics. The observations suggest that introducing a layered half-space does not improve the stress maps and ROC curves. However, results significantly improve for larger aftershocks and shorter time periods but without changing the ranking. We also go beyond binary testing and apply alternative statistics to test the ability to estimate aftershock numbers, which confirm that simple stress metrics perform better than the classic Coulomb failure stress calculations and are also better than the distance-slip probabilistic model.}, language = {en} } @article{PowaliSharmaMandaletal.2020, author = {Powali, Debarchan and Sharma, Shubham and Mandal, Riddhi and Mitra, Supriyo}, title = {A reappraisal of the 2005 Kashmir (M-w 7.6) earthquake and its aftershocks}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {789}, 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.2020.228501}, pages = {11}, year = {2020}, abstract = {We study the source properties of the 2005 Kashmir earthquake and its aftershocks to unravel the seismotectonics of the NW Himalayan syntaxis. The mainshock and larger aftershocks have been simultaneously relocated using phase data. We use back-projection of high-frequency energy from multiple teleseismic arrays to model the spatio-temporal evolution of the mainshock rupture. Our analysis reveal a bilateral rupture, which initially propagated SE and then NW of the epicenter, with an average rupture velocity of similar to 2 km s(-1). The area of maximum energy release is parallel to and bound by the surface rupture. Incorporating rupture propagation and velocity, we model the mainshock as a line source using P- and SH-waveform inversion. Our result confirms that the mainshock occurred on a NE dipping (similar to 35 degrees) fault plane, with centroid depth of similar to 10 km. Integrated source time function show that majority of the energy was released in the first similar to 20 s, and was confined above the hypocenter. From waveform inverted fault dimension and seismic moment, we argue that the mainshock had an additional similar to 25 km blind rupture beyond the NW Himalayan syntaxis. Combining this with findings from previous studies, we conjecture that the blind rupture propagated NW of the syntaxis underneath a weak detachment overlain by infra-Cambrian salt layer, and terminated in a wedge thrust. All moderate-to-large aftershocks, NW of the mainshock rupture, are concentrated at the edge of the blind rupture termination. Source modeling of these aftershocks reveal thrust mechanism with centroid depths of 2-10 km, and fault planes oriented subparallel to the mainshock rupture. To study the influence of mainshock rupture on aftershock occurrence, we compute Coulomb failure stress on aftershock faults. All these aftershocks lie in the positive Coulomb stress change region. This suggest that the aftershocks have been triggered by either co-seismic or post-seismic slip on the mainshock fault.}, language = {en} } @phdthesis{Sharma2024, author = {Sharma, Shubham}, title = {Integrated approaches to earthquake forecasting}, doi = {10.25932/publishup-63612}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-636125}, school = {Universit{\"a}t Potsdam}, pages = {xvi, 76}, year = {2024}, abstract = {A comprehensive study on seismic hazard and earthquake triggering is crucial for effective mitigation of earthquake risks. The destructive nature of earthquakes motivates researchers to work on forecasting despite the apparent randomness of the earthquake occurrences. Understanding their underlying mechanisms and patterns is vital, given their potential for widespread devastation and loss of life. This thesis combines methodologies, including Coulomb stress calculations and aftershock analysis, to shed light on earthquake complexities, ultimately enhancing seismic hazard assessment. The Coulomb failure stress (CFS) criterion is widely used to predict the spatial distributions of aftershocks following large earthquakes. However, uncertainties associated with CFS calculations arise from non-unique slip inversions and unknown fault networks, particularly due to the choice of the assumed aftershocks (receiver) mechanisms. Recent studies have proposed alternative stress quantities and deep neural network approaches as superior to CFS with predefined receiver mechanisms. To challenge these propositions, I utilized 289 slip inversions from the SRCMOD database to calculate more realistic CFS values for a layered-half space and variable receiver mechanisms. The analysis also investigates the impact of magnitude cutoff, grid size variation, and aftershock duration on the ranking of stress metrics using receiver operating characteristic (ROC) analysis. Results reveal the performance of stress metrics significantly improves after accounting for receiver variability and for larger aftershocks and shorter time periods, without altering the relative ranking of the different stress metrics. To corroborate Coulomb stress calculations with the findings of earthquake source studies in more detail, I studied the source properties of the 2005 Kashmir earthquake and its aftershocks, aiming to unravel the seismotectonics of the NW Himalayan syntaxis. I simultaneously relocated the mainshock and its largest aftershocks using phase data, followed by a comprehensive analysis of Coulomb stress changes on the aftershock planes. By computing the Coulomb failure stress changes on the aftershock faults, I found that all large aftershocks lie in regions of positive stress change, indicating triggering by either co-seismic or post-seismic slip on the mainshock fault. Finally, I investigated the relationship between mainshock-induced stress changes and associated seismicity parameters, in particular those of the frequency-magnitude (Gutenberg-Richter) distribution and the temporal aftershock decay (Omori-Utsu law). For that purpose, I used my global data set of 127 mainshock-aftershock sequences with the calculated Coulomb Stress (ΔCFS) and the alternative receiver-independent stress metrics in the vicinity of the mainshocks and analyzed the aftershocks properties depend on the stress values. Surprisingly, the results show a clear positive correlation between the Gutenberg-Richter b-value and induced stress, contrary to expectations from laboratory experiments. This observation highlights the significance of structural heterogeneity and strength variations in seismicity patterns. Furthermore, the study demonstrates that aftershock productivity increases nonlinearly with stress, while the Omori-Utsu parameters c and p systematically decrease with increasing stress changes. These partly unexpected findings have significant implications for future estimations of aftershock hazard. The findings in this thesis provides valuable insights into earthquake triggering mechanisms by examining the relationship between stress changes and aftershock occurrence. The results contribute to improved understanding of earthquake behavior and can aid in the development of more accurate probabilistic-seismic hazard forecasts and risk reduction strategies.}, language = {en} }