Statistical Methods and Modeling of Seismogenesis. Eleftheria Papadimitriou
the cells of each segment from A to E in a time span of 2,000 years. Vertical bars indicate the occurrence time of earthquakes whose ruptures significantly modified the stress of more than one segment. The increase of stress is mostly due to tectonic loading, but is also due to stress transfer from the cells of a segment to the others. The average stress drops at the time of each earthquake according to the size of the rupture on the specific segment participating in the earthquake.
Figure 2.18(b) shows the standard deviation (again in a normalized scale) of stress on the cells of each segment. Figure 2.18(c) shows the ratio between the values of the average stress and its standard deviation. This parameter increases always and only before a strong earthquake. Assuming that the trend shown in Figure 2.17 before large earthquakes (i.e. stress becomes higher and more uniform on a relevant patch of the fault system) is a phenomenon that happens in reality, it could be regarded as a promising tool for forecasting strong events. Unfortunately, it is not feasible to map the stress on the fault surface in the real environment with the necessary resolution. However, some seismicity patterns observable by means of modern high-quality seismic networks, such as spatio-temporal changes of b-value, could be associated with stress variation in the crust and applied in an operational forecasting system (Montuori et al. 2016; Gulia et al. 2016; Gulia and Wiemer 2019).
Figure 2.17. Stress time history on the Nankai seismogenic structure for the first 76 years. For a color version of this figure, see www.iste.co.uk/limnios/statistical.zip
2.5. Appendix 1: Relations among source parameters adopted in the simulation model
Figure 2.18. (a) The time history (in a normalized scale) of the average stress computed on all the cells of each segment from A to E in a time span of 2,000 years. Vertical bars indicate the occurrence time of earthquakes whose ruptures significantly modified the stress of more than one segment. The increase in stress is mostly due to tectonic loading, but is also due to stress transfer from the cells of a segment to the others. The average stress drops at the time of each earthquake according to the size of the rupture on the specific segment participating in the earthquake. (b) The same as in the top panel for the standard deviation of stress on the cells of each segment. The standard deviation decreases (i.e. the stress becomes more uniform) when the occurrence time of earthquakes approaches. (c) Ratio between the values of the average stress and its standard deviation. This parameter increases always and only before a strong earthquake, as a possible precursor of large size events. For a color version of this figure, see www.iste.co.uk/limnios/statistical.zip
This appendix provides the theoretical framework of the model adopted in the simulation algorithms described in section 2.3. Because the magnitudes of the earthquakes of the synthetic catalog generated by the algorithm are obtained from their seismic moments, and the number of these earthquakes per unit time must fit the constraint of the slip rate assigned to any fault segment of the model, it is necessary to adopt relationships between the magnitude and both the rupture area and the average slip of each event.
The scalar seismic moment of an earthquake is defined as:
where μ is the shear modulus of the elastic medium,
where Δσ is the static stress drop of the earthquake. From [2.5] and [2.6], we obtain:
Two widely used formulas link the seismic moment and the seismic energy of an earthquake with magnitude:
Substituting [2.8] and [2.9] into [2.7], we obtain:
[2.10]
and adopting a value μ = 3.3 • 1010Pa for rocks in the Earth’s crust:
[2.11]
For a rectangular rupture, in analogy with the theoretical formulas introduced for a circular fault by Keilis Borok (1959), the following relations were proposed by Console and Catalli (2006):
[2.12]
[2.13]
where x and y are the coordinates on the fault plane, counted from the center of the rectangular source, D(x, y) is the displacement on the point x and y coordinates and Dmax is the displacement on the fault center. From the definition of seismic moment [2.5]:
Combining [2.5],