Abstract
The frequency-magnitude distribution follows the Gutenberg-Richter empirical law, in which the scaling between small and large earthquakes is represented by the b-value. Laboratory experiments have shown that the b-value is related to fault mechanics with an inverse dependency to the differential stress, as was also inferred from observational datasets through relations with earthquake depth and style of faulting. In this study, we aim to obtain a better understanding of the geological structure and tectonics along the Dead Sea transform (DST), by examining relations of the b-value to three source parameters: the earthquake depth, the seismic moment release, and the predominant style of faulting. We analyse a regional earthquake catalogue of ∼20,300 earthquakes that were recorded between 1983 and 2020 in a regional rectangle between latitudes 27.5°N−35.5°N and longitudes 32°E−38°E. We convert the duration magnitudes, Md, to moment magnitudes, Mw, applying a new regional empirical relation, by that achieving a consistent magnitude type for the entire catalogue. Exploring the variations in the b-value for several regions along and near the DST, we find that the b-value increases from 0.93 to 1.19 as the dominant style of faulting changes from almost pure strike-slip, along the DST, to normal faulting at the Galilee, northern Israel. Focusing on the DST, our temporal analysis shows an inverse correlation between the b-value and the seismic moment release, whereas the spatial variations are more complex, showing combined dependencies on seismogenic depth and seismic moment release. We also identify seismic gaps that might be related to locking or creeping of sections along the DST and should be considered for hazard assessment. Furthermore, we observe a northward decreasing trend of the b-value along the DST, which we associate to an increase of the differential stress due to structural variations, from more extensional deformation in the south to more compressional deformation in the north.
1 Introduction
One of the most fundamental observations in earthquake seismology is that the frequency-magnitude relation follows the Gutenberg-Richter empirical law (Gutenberg and Richter, 1944):where N is the cumulative number of earthquakes of at least a magnitude M; and a and b are seismicity parameters.
Whilst the a-value reflects the seismic activity level and unless normalised, varies in different time windows, the b-value reflects the proportion between small and large earthquakes, and thus has a significant impact for hazard evaluations (Frankel, 1995; Petersen et al., 2011; Marzocchi and Taroni, 2014; Magrin et al., 2017; Sokolov et al., 2017). Yet, some observations suggest that the frequency-magnitude relation (Eq. 1) might deviate at the upper bound of the magnitude distribution, since some fault zones are characterised with repeated large earthquakes of approximately the same size (i.e., “characteristic” behaviour; Wesnousky et al., 1983; Schwartz and Coppersmith, 1984).
Although a b-value of a unity has been previously suggested on a global scale based on real-data analyses (Frohlich and Davis, 1993; Felzer et al., 2004; El-Isa and Eaton, 2014) and theoretical considerations (King, 1983), other studies have been aimed to understand its variations with regards to earthquake mechanics (e.g. Scholz, 1968, 2015; Main et al., 1989; Henderson and Main, 1992; Lei et al., 2000; Tan et al., 2019). It has been previously demonstrated that the b-value varies with the style of faulting, with typical values of 0.7–0.8, 0.9–1.0, and 1.1–1.2 for thrust, strike-slip and normal faulting, respectively, with values in-between for oblique faulting (Schorlemmer et al., 2005; Petruccelli et al., 2019a). This dependency has been shown by other studies as well (Gulia and Wiemer, 2010; Scholz, 2015; Bora et al., 2018; ), suggesting an inverse correlation of the b-value with the differential stress, considering the Anderson theory of faulting (). The inverse dependency of the b-value with the differential stress has been also shown through rock fracture experiments (Scholz, 1968; ; Rivière et al., 2018), and dependency with the earthquake depth considering rheological models (Mori and Abercrombie, 1997; Gerstenberger et al., 2001; Spada et al., 2013; Scholz, 2015; Rigo et al., 2018). In such cases, the b-value has been observed to decrease with depth, within the brittle part of the earth’s crust.
Previous estimations of the frequency-magnitude relation in the region of Israel (Ben-Menahem, 1981, ; Salamon et al., 1996; Shapira and Hofstetter, 2002; Hamiel et al., 2009) were mainly focused on hazard perspective, and had limited data for investigating variations of the b-value. In this study we systematically calculate the b-value, focusing on ∼500-km of the southern part of the Dead Sea transform (DST), a ∼1,000-km long continental transform plate boundary that links between the Red Sea spreading centre and the convergence zone in southern Turkey (e.g., Garfunkel, 2014). Our purpose is to investigate variations of the b-value and their relation to the earthquake depth, the style of faulting, and the seismic moment release; hence examining whether and to what extent does the b-value relate to mechanical properties of the fault zone. We examine spatial variations of the b-value with the seismogenic depth, which we estimate according to the 75th and 95th percentiles of the seismicity depth distribution, because they reflect the alteration in the seismogenic depth along the DST (e.g., Shalev et al., 2013; Wetzler and Kurzon, 2016). Possible correlations with the seismic moment release are examined in both space (e.g., Bora et al., 2018) and time (e.g. Cao and Gao, 2002) analyses, for receiving insights about the seismicity of the region and for further understanding variations of the b-value.
Determination of the b-value is often limited within spatial zones, based on seismological or tectonic considerations, which in many cases reflect the dominant style of faulting (e.g. Radulian et al., 2000, 2018; ; Bus et al., 2009) and may also be implemented for hazard evaluation purposes (Frankel, 1995; Helmstetter et al., 2006; Yadav et al., 2012; ; Maiti and Kamai, 2020; Mandal et al., 2021; Yagoda-Biran et al., 2021). A preliminary seismogneic zonation in the region (Shamir et al., 2001) and its following analysis of the frequency-magnitude relation (Shapira and Hofstetter, 2002) were based on a rather sparse seismological dataset. A more recent seismogenic zonation (Sharon, 2020) was based on the density distribution of epicentres and seismic moment release, and their spatial relation to the main seismic sources (Sharon et al., 2020). However, this zonation is not adequate to show systematic relation between the b-value and the faulting style because many of these previous zones are too small to achieve well-determined b-value with our current data. For this purpose, we introduce here a rather simplified and broad new tectonic zonation in the DST and its periphery, according to the seismic network capability and the local tectonics.
In this study we first relocate the earthquakes recorded between January 1983 and September 2020 (Wetzler and Kurzon, 2016). Then, due to inconsistent magnitude type within the original catalogue, we generated a catalogue, choosing Mw as the preferred homogenous magnitude type, since it is physical-based and is not saturated at high magnitudes (Kanamori, 1977; Hanks and Kanamori, 1979). Subsequently, we estimate the completeness magnitude, , and calculate the frequency-magnitude parameters with respect to the seismogenic depth, the seismic moment release and the tectonic regime associated with the prevailing style of faulting. Specifically, we perform a systematic and thorough spatio-temporal investigation of the b-value at the DST and its periphery.
2 Tectonic settings
The DST was formed during the Miocene, as the African-Arabian plate broke, generating the Suez rift and the DST. While the Suez rift has shown minor signs of post-Miocene deformation, the DST is considered to be the main source of post-Miocene deformation in the region (Garfunkel and Bartov, 1977; Joffe and Garfunkel, 1987; Steckler et al., 1988). It consists of a ∼1000-km long ∼N-S orientated fault system, which is the largest in the Levant (Figure 1). Evaluation from geologic and geodetic sources indicate Quaternary slip rates of 4–5 mm/yr (Garfunkel, 2010; Sadeh et al., 2012; Marco and Klinger, 2014; Hamiel et al., 2018). Our study focuses on the southern section of the DST (Figure 1), dominated by a left-lateral strike-slip overall displacement of ∼105-km accumulated over the past ∼16–20 million years (Quennell, 1959; Garfunkel, 1981, 2014; Nuriel et al., 2017).
FIGURE 1
The lateral motion on the DST occurs on left-stepping strike-slip and oblique-slip fault segments that delimit a string of en-echelon arranged pull-apart basins (Garfunkel, 1981; Zak and Freund, 1981; Garfunkel and Ben-Avraham, 2001). The DST is topographically expressed by a pronounced 5–25 km wide valley, bordered by normal faults that extend along the valley margins. The north-eastern edge of the study area comprises the Lebanon restraining bend (LRB; Figure 1), where the DST is branched into several segments, transferring the strike-slip motion into the Lebanon area (Gomez et al., 2003, 2007). The northern section of the DST crosses northwest Syria in a N-S orientation, most of it outside our study area.
South of Lebanon, the Sinai sub-plate has several fault systems, associated with Quaternary internal-deformation: the Carmel-Tirza fault zone (CTF; Figure 1) divides the Israel-Sinai sub-plate into two tectonic domains (Neev et al., 1976; ; Sadeh et al., 2012) where the southern part is more rigid, while the northern consists of a set of graben-and-horst structures with E-W-striking normal faults associated with S-N extension (Ron and Eyal, 1985). The CTF consists of SE-NW orientated fault segments, with normal and oblique motions (Freund, 1970). It is associated with coeval motion of ∼0.7 mm/yr left-lateral slip and ∼0.6 mm/yr extension rates (Sadeh et al., 2012), with recent seismicity in its eastern side that sprawl over several parallel fault segments (Hofstetter et al., 1996; Sharon et al., 2020).
To the south of the CTF, several ∼E-W striking faults are associated with mainly dextral slip and some normal faulting (Bentor and Vroman, 1954; ; Zilberman et al., 1996), which occurred mainly during the Neogene or in earlier periods (Weinberger et al., 2020). An additional fault system of ∼NNE striking faults in southern Israel is associated with normal faulting and minor extension component along the transform system, and was more active during Quaternary times (; , ; Calvo and Bartov, 2001; Marco, 2007).
The potential for strong earthquakes along the DST is demonstrated by the earthquakes of the 1927 ML 6.2 near Jericho, north of the Dead Sea (Shapira et al., 1993), and the 1995 MW 7.2 Nuweiba (Hofstetter et al., 2003) in the Gulf of Elat (Aqaba), along with pre-instrumental records, leading to estimations of up to MW ∼7.5 (; ; Marco and Klinger, 2014; Zohar et al., 2016, 2017; Lu et al., 2020). Deep-crust seismicity, correlated with low heat flow areas, particularly in the Dead Sea basin, probably indicates a cold crust, with deep brittle to ductile transition zone (; Shalev et al., 2007, 2013; ; Wetzler and Kurzon, 2016).
3 Seismological dataset
We analyse an earthquake catalogue from 1 January 1983 until 28 September 2020, recorded by ∼210 stations. The majority of the data originates in the stations of the Israel Seismic Network (ISN), including many new stations deployed under the framework of the TRUAA network, established for the purpose of Earthquake Early Warning System for the state of Israel (Kurzon et al., 2020). In addition, some of the data comes from local CTBT and CNF stations (Comprehensive Nuclear Test-Ban Treaty, and Cooperating National Facility, respectively), and a minority originates in stations of other networks: the GEOFON global network of GFZ, the JSO seismic observatory of Jordan and the CQ seismic network of Cyprus.
In the original catalogue, documented by the ISN, there are 23,316 earthquakes between the latitudes and longitudes , referred here as the Rectangular Area (RA; Figure 1). Following Wetzler and Kurzon (2016), we applied similar procedure and methods for earthquake relocation: applying the regional velocity model of Gitterman et al. (2005), earthquakes are located by the genloc library (Pavlis et al., 2004) of the Antelope seismic software package (<www.brtt.com>). Horizontal median errors are 420 m in longitude, 510 m in latitude, and the vertical median error is 750 m, calculated according to Pavlis (1986). The relocated catalogue has ∼23,200 events (Figure 1), with improved locations, and loosing less than 0.5% of the events due to large location errors. The magnitude range of the catalogue is , from which ∼2,900 are with unknown magnitudes (Table 1).
TABLE 1
| None | ||||
|---|---|---|---|---|
| RA | 16,546 | 5,537 | 1,756 | 2,895 |
| NCA | 3,965 | 2,895 | 810 | 938 |
The number of and magnitudes, determined for the investigated catalogue. The right-most field is for events that contain no magnitude.
The catalogue includes two magnitude types: duration magnitude (Md) and the moment magnitude (Mw). Table 1 summarises the two types of magnitude that were determined for the relocated catalogue. The magnitudes of the Mw 7.2 1995 Nuweiba earthquake and the Mw 5.1 2004 Dead Sea earthquake were fixed according to Hofstetter et al. (2003) and Hofstetter et al. (2008), respectively.
The detection sensitivity of the seismic network is primarily dependent on the background noise and the distribution of seismic stations. The seismic network coverage area (NCA; Figure 2) was recently determined by Sharon et al. (2020) according to the total threshold number of wave arrivals at seismic stations, accounting for their spatial distribution, hence, including hypocentres that are relatively well-constrained, with relatively low magnitude of completeness (. Therefore, determination of the frequency-magnitude parameters within the NCA can be based on more data, consisting of a larger range of magnitudes.
FIGURE 2
4 The frequency-magnitude relation
4.1 Magnitude conversion
Earthquake magnitude can be estimated by a wide range of methods and parameters, depending on the spectral properties of the source, the seismic phases, the instrumentation capabilities, and the consideration of path and site effects. Therefore, magnitude estimation does not behave uniformly for all magnitude ranges, and also saturates at different levels (Kanamori, 1977, 1983; Utsu, 2002). Thus, magnitude conversion to a single type of magnitude is vital for seismological analyses that require homogenous earthquake catalogue with consistent magnitude type, such as analysis of the frequency-magnitude relation.
The moment magnitude (Mw) is based on the physical dimensions of the seismic source, does not saturate at extreme rupture size, and hence is more adequate to represent a wide range of magnitudes (Kanamori, 1977; Hanks and Kanamori, 1979; Choy and Boatwright, 1995). Therefore, a common practice is converting other magnitude types to Mw (Scordilis, 2006; Yadav et al., 2009, 2012; Ross et al., 2016; Kumar et al., 2020), which is widely used as a unified magnitude for seismological and hazard-related applications (e.g. Bormann and Di Giacomo, 2011). In the original Israel catalogue, about 70% of the events are assigned only with Md, without Mw, and less than 10% are estimated by both magnitude types. Therefore, by converting Md to Mw, we intent to achieve a catalogue of a uniform magnitude type, Mw, that can be the basis for the current seismological analysis, as well as for future investigations.
Considering that many of the magnitude types are determined by a set of parameters, varying according to geological, seismological and instrumentational settings, there is no global formula that converts between and ; hence, unique formulas are deduced for different regions. A few regression techniques have been suggested and employed for converting magnitudes, and the justification for applying a single method or the other, are based either on theory or empirical observations (e.g. Yadav et al., 2012; Kadirioğlu and Kartal, 2016; Das et al., 2018). In this study we examine Md to Mw conversion, applying both, ordinary least squares (OLS) regressions, (e.g. Scordilis, 2006; Castello et al., 2007; ; Ross et al., 2016), and orthogonal regression (OR; e.g., Glaister, 2005; Castellaro et al., 2006; Kane and Mroch, 2020).
Although it is more common to convert magnitudes through linear regressions, added a quadratic regression and obtained better fit to observations in the same region. We examine their approach for Md to Mw conversion, fitting both linear and quadratic regressions. For comparison, we also examine part of their obtained regressions on the much larger data utilised here.
The regression coefficients published by were deduced in the forms of:andfor the linear and quadratic regressions, respectively. We convert M0 to Mw following the relationship (Hanks and Kanamori, 1979; ):referring to units of N-m (e.g. Shearer, 2009).
In this manner we achieve the standard coefficient correlations summarised in Table 2 along with our current results. Examining the performance of the linear and quadratic fits using regression parameters, in comparison to the fits obtained in the current study, we note a clear improvement in the goodness of fit, observed by the lower values of the Root Mean Square error (RMSE) of the current study (Table 2). This result is not surprising, as results are based only on ∼100 earthquakes, compared with ∼1,800 earthquakes for the current results (Table 1).
TABLE 2
| a | |||||
|---|---|---|---|---|---|
| Regression type and data source | a | b | c | RMSE | AICc |
| Linear () | 0.90 | 0.11 | - | 0.287 | - |
| Quadratic () | 0.13 | −0.13 | 2.11 | 0.400 | - |
| Linear (RA) | 0.81 | 0.52 | - | 0.170 | −4.46e + 03 |
| Quadratic (RA) | 0.03 | 0.65 | 0.69 | 0.167 | −4.52e + 03 |
| b | |||||
|---|---|---|---|---|---|
| Regression type | a | b | c | RMSE | AICc |
| Linear Orthogonal (RA) | 0.77 | 0.64 | - | 0.139 | −6.71e + 03 |
| Quadratic (RA) | 0.03 | 0.65 | 0.69 | 0.132 | −6.90e + 03 |
Coefficient correlation parameters (a, b, c) that correspond to linear and quadratic regressions of the forms and , respectively; RMSE and AICc of these regressions are in respect to earthquakes located within RA; a) The regression coefficients of are from S-wave analysis and the magnitude range of , and were achieved through Eqs 1–3 (see Methods section); b) Comparison between orthogonal regression and quadratic OLS regresion.
In Figure 3, we compare between two OLS regressions, linear and quadratic, finding a relatively lower RMSE values obtained by the OLS quadratic fit (Figure 3C) compared with the linear OLS fit (Figure 3A). The parameters retain stable residuals around zero throughout the entire magnitude range for both regression types (Figure 3B, D). However, both the lowest and highest magnitude ranges suggest that the quadratic fit (Figure 3D) shows more constrained convergence around the line y=0 in comparison to some bias in the linear fit (Figure 3B). In addition, the total RMSE of the quadratic fit is slightly smaller (Table 2), and also the AICc method (Hurvich and Tsai, 1989), based on the Akaike Information Criteria (Akaike, 1974), shows lower values for the quadratic fit (Table 2). We also examined a linear orthogonal regression (OR), comparing it to the quadratic OLS fit. The quadratic OLS has better fit, reflected in its lower RMSE and AICc values (Table 2).
FIGURE 3
The stability of the regressions and corresponding residuals with detailed RMSE were also tested for earthquakes within NCA (Supplementary Figure S1), showing similar results (Table 2; Supplementary Table S1). For a more robust conversion that is based on a larger number of earthquakes, we choose to use the parameters obtained by the entire research area (RA); hence, our Md to Mw conversion formula (Table 2) is:
4.2 The magnitude of completeness and the b-value
The magnitude of completeness (Mc) is the minimum magnitude of which the seismic network detects all the events. Therefore, it marks the point from which the frequency-magnitude distribution is linear. We apply a few algorithms, by Goebel et al. (2017), following and Clauset et al. (2009), based on Kolmogorov-Smirnov test, and by Mizrahi et al. (2021), is also based on Clauset et al. (2009), deducing Mc with high statistical robustness. Both algorithms suggest a completeness magnitude of 2.1 for the NCA. Similarly, in a previous work (Sharon, 2020), Mc was estimated as 2.0 for the NCA, for the years 1983–2017. This is in consent with prior estimations of Shapira (1992), obtaining Mc = 2.0 for a region, approximately overlapping the NCA, between the years 1984–1991. The predominant (or even the only) magnitude type in the analyses of Sharon (2020) and Shapira (1992) was Md. Our conversion formulation (Eq. 5) indicates that Md=2.0 is approximately equivalent to Mw of 2.1. Thus, we conclude that the magnitude of completeness for our homogenised catalogue is 2.1. As the NCA dataset is the one used later for more detailed analysis, by plotting the moment magnitude as a function of the sequential number of events (e.g., Zhuang et al., 2017; Bustos et al., 2022), we demonstrate that there are no clear short-term periods of magnitude incompleteness, within the NCA dataset (Figure 4).
FIGURE 4
In general, higher Mc is expected for remote seismicity for which the seismic sensitivity decreases. This is demonstrated for earthquakes recorded in the entire research area (RA), where both algorithms of Goebel et al. (2017) and Mizrahi et al. (2021) suggest that Mc = 3.8.
We calculate the b-value through the Maximum Likelihood Estimation by the following equation (Marzocchi and Sandri, 2003; following Fisher, 1950; Utsu, 1966; and Bender, 1983):where is the sampling average of the magnitudes, Mc is the completeness magnitude, and is the magnitude bin interval in which the data is examined. The standard error of the b-value is obtained from the following formula (Shi and Bolt, 1982):where n is the number of earthquakes.
We calculate the frequency-magnitude parameters for several zones (Figures 5, 6): first, for regional zones: RA and NCA, and then for zones that differ in their local tectonics (Figure 5; Table 3). The zone of the DST is defined here by a ∼25-km width polygon representing the deformation zone of the southern section of the DST (excluding the Gulf of Elat due to the lack of seismic coverage). The width of the polygon increases up to 28-km in pull-apart basins, and decreases to 20-km at more localised sections that consist of long straight fault segments. For the rest of the NCA, the polygon of the CTF seismogenic zone (Sharon, 2020) is adopted to divide the off-fault seismicity into southern and northern Off-Fault Zones, OFZ(S) and OFZ(N), respectively (Figure 5). These two zones are also bounded by the DST polygon from the east and hence reflect two local tectonic provinces (Neev et al., 1976; ; Sadeh et al., 2012) that are separated from the main active fault zones. The frequency-magnitude parameters are presented in Table 3 for each of these tectonic zones (Figure 5), and their detailed plots are provided in Figure 6. In addition, the frequency-magnitude parameters of the seismogenic zones (Sharon, 2020) are provided in Supplementary Table S2.
FIGURE 5
FIGURE 6
TABLE 3
| Zone | a-value | b-value | Events | Magnitude range (Mc in bold) | M0 release(N*m) | 75th depth percentile (km) | 95th depth percentile (km) |
|---|---|---|---|---|---|---|---|
| RA | 660 | 9.90e + 18 | 34.0 | 45.3 | |||
| NCA | 1561 | 3.16e + 17 | 16.8 | 25.3 | |||
| DST | 1301 | 3.00e + 17 | 17.1 | 25.2 | |||
| CTF | 368 | 1.27e + 17 | 14.7 | 21.9 | |||
| OFZ(S) | 154 | 2.05e + 16 | 19.2 | 28.4 | |||
| OFZ(N) | 80 | 2.74e + 15 | 7.9 | 18.1 |
5 Variations of the b-value
5.1 Spatial variations of the b-value
We first examine whether the differences in the b-values between the tectonic zones (Table 3; Figure 5) reflect different style of faulting. We examine only the tectonic zones, within the NCA, all affected by a similar network coverage, reflected also by similar magnitude of completeness values of Mc=2.1 (Figure 5). The strike-slip dominated DST (Garfunkel, 1981, 2014; Hofstetter et al., 2007; Marco and Klinger, 2014) shows a b-value of 0.93. The CTF zone, which accommodates an extensional and strike-slip associated deformation (Freund, 1970; ; Rotstein et al., 1993; Hofstetter et al., 1996; Sadeh et al., 2012), has a b-value of 0.99. The Galilee area, marked by OFZ(N), is mostly extension-dominated normal faulting (Freund, 1970; Ron et al., 1984; Hofstetter et al., 2007), and shows a b-value of 1.19. The OFZ(S), an area of sparse seismicity associated with strike-slip and extensional structures (Bentor and Vroman, 1954; ; Zilberman et al., 1996; ; Ginat et al., 2000, 2002), shows a b-value of 0.96. These b-values, including their errors (Shi and Bolt, 1982), show a subtle preference of higher b-values for extensional faulting. This trend is most prominent for the OFZ(N) (Figure 5) with a b-value of 1.19. We further explore the significance of this trend through two statistical methods: the Utsu’s method (Utsu, 1999; e.g., Xie et al., 2019) and the Student’s t-test (Gosset, 1908; e.g., Petruccelli et al., 2018). In both methods, we test the statistical significance for ten pairs, comprising the relations between the five tectonic zones (including the “parent” NCA zone), and present the results by significance level matrices (Table 4), in which the significance level () is between 0 and 1; is highly significant and indicates identical datasets (e.g., Utsu, 1999).
TABLE 4
| NCA | DST | CTF | OFZ(S) | OFZ(N) | |
|---|---|---|---|---|---|
| A. Utsu’s test significance level matrix | |||||
| NCA | 1 | 0.257 | 0.319 | 0.367 | 0.070 |
| DST | 0.257 | 1 | 0.211 | 0.343 | 0.043 |
| CTF | 0.319 | 0.211 | 1 | 0.349 | 0.126 |
| OFZ(S) | 0.367 | 0.343 | 0.349 | 1 | 0.112 |
| OFZ(N) | 0.070 | 0.043 | 0.126 | 0.112 | 1 |
| B. Student’s t-test significance level matrix | |||||
| NCA | 1 | <0.01 | 0.014 | 0.328 | <0.01 |
| DST | <0.01 | 1 | <0.01 | 0.085 | <0.01 |
| CTF | 0.014 | <0.01 | 1 | 0.120 | <0.01 |
| OFZ(S) | 0.328 | 0.085 | 0.120 | 1 | <0.01 |
| OFZ(N) | <0.01 | <0.01 | <0.01 | <0.01 | 1 |
| C. Color scheme | |||||
| s.l<0.01 | |||||
| 0.01<s.l<0.05 | |||||
| 0.05<s.l<0.2 | |||||
| 0.2<s.l<1 | |||||
| s.l=1 | |||||
Analysis of statistical significance for the Tectonic Zones, showing the significance level () matrices according to a) Utsu’s test; and b) the Student’s t-test. The statistical significance follows a color scheme (Table 4c) defined by the range of , in which the green shades show threshold required for obtaining statistical significance. While the Utsu’s test show high similarity between the datasets, and therefore except for one pair (DST ←→OFZ(N)) does not indicate statistical significance, the Student’s t-test shows that seven out of 10 pairs are significantly different from each other, with . The stages of calculation for the t-test are provided in the Tables A3 and A4 of the Supplementary material.
In Table 4 we present the results obtained by Utsu’s method. In this method, we first compute the AIC (Akaike Information Criteria; Akaike, 1974) for each examined dataset. Then the difference between each dataset-pair () is calculated for each dataset, and the significance level is given by the probability (Utsu, 1999); when p the two examined datasets are significantly different. According to Utsu’s method only one pair shows clear statistical significance (DST←→OFZ(N); Table 4). In Table 4 we present the results obtained by the Student’s t-test method. The type of t-test is chosen according to the ratio between each pair’s b-value standard deviations, (Rice, 2006) defining whether the variances of the b-values are similar or not between two tectonic zones (Supplementary Table S3). By setting the statistical significance requirement to (e.g., Utsu, 1999) and by using the Student’s t-distribution table (Supplementary Table S4), we determine whether the observed differences between the b-values of two tectonic zones are statistically significant or a matter of chance, hence insignificant (Table 4). In addition to Utsu’s method, the Student’s t-test method suggests that except for the OFZ(S) related pairs, all the other pairs are significantly different from each other, meaning that the differences between the b-values are not randomly distributed. Hence, they may be attributed to physical causes, such as the predominant faulting style (Schorlemmer et al., 2005). Such differences between the two methods have been observed, for example, by Petruccelli et al. (2018), as the Student’s t-test method has shown more dataset-pairs with statistical significance than the Utsu’s method, for all the intermediate cases (those that do not reflect a pure faulting style; see Tables 3, 4 in Petruccelli et al., 2018).
5.2 Spatial variations of the b-value along the DST
We further explore in detail the spatial variations of the b-value along the main tectonic feature, within the DST polygon (Figure 5), by systematically scanning it in a 60-km long latitudinal bands, at 1-km increments, using the DST’s polygon completeness level of Mc=2.1 (Table 3; Supplementary Figure S2). These spatial bands represent the approximate minimum length to capture at least 100 events at each sample, with only one gap in between the bands (Figure 7). The b-value and its error (Shi and Bolt 1982) are calculated at each sample, where the mean number of events in these spatial bands is 185. An additional sparse-seismicity profile (50–99 events per sample; Figure 7) aligns with the main profile, showing a consistent behaviour of the b-value profile along the DST. In addition, the error bar band around the profile shows that the changes in the b-value are significant, substantially larger than the errors. The same spatial profiling is generated for the 75th and 95th hypocentre depth percentiles, and for the accumulated seismic moment converted from the magnitudes (Figure 7). The depth percentiles mark the depth profiles for which 75% and 95% of the events within the spatial bands are shallower, respectively. The accumulated seismic moment, which is the sum of seismic moments of all events within the spatial bands, is calculated at each sample.
FIGURE 7
The profiles presented in Figure 7 show rather complex correlations between the b-value and both the seismogenic depth and seismic moment release. We observe an inverse correlation of the b-value with the seismogenic depth along the DST, with lower b-values for deep seismogenic zones (i.e., latitudes 30.3°–30.7°, 30.9°–31.3°). A similar trend is presented in Supplementary Figure S3) for the tectonic (Figure 5) and seismogenic zones (Sharon, 2020; provided also in the Supplementary Figure S4). In addition, the seismic moment release shows at the southern section an inverse correlation with the b-value (latitudes 29.3°–31.0°; Figure 7); for example, this is emphasized at the “quasi-spike” of the b-value in latitude ∼30.8. In contrast, the b-value increase at ∼31.7, cannot be clearly correlated either to the seismogenic depth or to the seismic moment release, and seems to correspond to the Mw 5.1 mainshock (Hofstetter et al., 2008) with its aftershock sequence (see Figure 7A). This increase may reflect delocalization at the northern Dead Sea basin associated with its structural complexity (e.g., Garfunkel and Ben-Avraham, 1996). Yet, with all the complexity of the b-value profile, the overall trend shows that the b-value decreases northward along the DST (dashed blue linear trend line in Figure 7).
5.3 Temporal variations of the b-value along the DST
Since the relationship between the b-value and the seismic moment release has not been comprehensively studied, it is not clear whether the inverse spatial correlation (southern section, Figure 7) is a typical relationship, or whether it is biased by specific significant seismic activity (e.g., mainshocks, aftershocks, swarms). Therefore, some insights may be obtained by examining the temporal relations between the b-value and the seismic moment release. Statistical methods have been developed for examining high-resolution temporal variations of the b-value (e.g. Wiemer et al., 1998; Mousavi et al., 2017). However, applying them requires focusing on specific sub-zones, reducing to too small datasets in the case of the DST, due to the relatively low-intensity seismicity. Therefore, we investigate the temporal correlation of the b-value with the seismic moment release (Figure 8), by applying a rather simple technique, calculating the b-value in a moving time window, with a fixed number of 100 events per window, and a time interval of five events; also here, the DST’s polygon completeness level of Mc2.1 is applied (Table 3; Supplementary Figure S5). Similar temporal analysis of the b-value, using a fixed number of events, has been shown and discussed in other studies (e.g., Nuannin et al., 2004; Tormann et al., 2013). In addition, as the catalogue consists of earthquakes recorded continuously by the ISN, we assume the data is inherently homogenous, with no temporal bias that can be caused by data comprising a combination of different catalogues (e.g. Zúñiga and Wiemer, 1999; Tormann et al., 2010).
FIGURE 8
The temporal analysis (Figure 8) mainly shows an inverse correlation between the b-value and the seismic moment release (Figure 8A). We also examine the b-value in specific time windows, with a minimum of 450 events per time window, achieving higher statistical confidence. Their frequency-magnitude plots are presented in Figures 8B–E. These time-windows deduce b-values that fit the fluctuating b-value profile obtained by the fixed-amount-of-events method.
6 Discussion
6.1 Relation to the differential stress and faulting mechanism
The seismogenic depth along the DST fluctuates correlatively with the thermal profile (Shalev et al., 2013; Wetzler and Kurzon, 2016). The >25-km deep seismicity in some parts (Figure 7), particularly in the Dead Sea basin, indicates a cold crust with a deep brittle-to-ductile transition zone (; Shalev et al., 2013; ; Wetzler and Kurzon, 2016), whereas shallow brittle-to-ductile transition occurs in zones of shallower seismogenic depth. The differential stress within the brittle crust increases with depth (e.g. Sibson, 1974, 1984; Llana-Fúnez and Lõpez-Fernández, 2015) as the confining pressure grows (Byerlee, 1968). Hence, deeper seismogenic zones are associated with higher differential stress because they include deeper portion of earth’s (brittle-) crust. We therefore interpret the inverse correlation between the b-value and the seismogenic depth (Figure 7) as related to the differential stress within the brittle crust. This interpretation is consistent with previous laboratory () and seismic observations (Gerstenberger et al., 2001; Spada et al., 2013; Scholz, 2015; Petruccelli et al., 2019b), which showed inverse correlation between the b-value and the depth of earthquakes, and hence with their corresponding differential stress.
The increase of the b-value in tectonic zones, as the predominant faulting style changes from strike-slip dominated DST to normal faulting dominated OFZ(N) (Figure 5), is interpreted here as another expression of the inverse dependency of b-value with differential stress environments (Schorlemmer et al., 2005; Llana-Fúnez and Lõpez-Fernández, 2015; Petruccelli et al., 2019b). This relation is shown by normal faults that are associated with relatively lower differential stresses and higher b-values, while strike-slip faults that are associated with higher differential stresses, and lower b-values.
We also observe an overall trend of the b-value decreasing northwards, within the DST polygon (Figure 7). This trend may be caused by variations of mechanical parameters such as the internal friction angle and faulting geometry (; Petruccelli et al., 2019a). Such variations can be linked to geometrical changes of the plate boundary (Garfunkel, 1981; Joffe and Garfunkel, 1987) and perhaps a transition from the extensional pull-apart basins of the Gulf of Elat (Reilinger et al., 2006) in the southern end of the DST polygon, closer to the Red Sea rift, to the transpressional LRB in the north (Gomez et al., 2007; Weinberger et al., 2009). Differences in the differential stress associated with these structural variations may also be related to post-Neogene widening of the fault zone along the southern section of the DST in the study area (latitudes<31°; Figure 7), while its northern section (latitudes>33°; Figure 7) underwent convergence and shortening (; Marco, 2007; Weinberger et al., 2009; Garfunkel, 2010).
6.2 The b-value and the seismic moment release
The b-value shows an inverse correlation with the seismic moment release, at least in the southern section of the DST (latitudes<31°; Figure 7), and a clearer inverse correlation is seen in the temporal analysis (Figure 8). The inverse correlation aligns with the theoretical study of Wyss (1973) and with observations from other tectonic domains (e.g., Cao and Gao, 2002; Bora et al., 2018).
Seismic moment release is affected by two main factors: the accumulation of seismic moment of small events, and the addition of seismic moment from more significant events. The latter, when exists, tends to dominate over the former. In many aftershock sequences it has been observed that the b-value increases (Suyehiro, 1966; Gibowicz, 1974; King, 1983; Wiemer and Katsumata, 1999; Gulia et al., 2018), due to the stress relaxation after the mainshock and activation of fault branches at the periphery of the rupture, while some have shown a decrease in the b-value around the mainshock (e.g., Mukhopadhyay et al., 2013). A decrease in b-value prior to the main shock has been observed in other studies, and was associated to an increase of the differential stress (e.g. Schorlemmer and Wiemer, 2005; Rivière et al., 2018).
In the case of the DST, the seismicity observed here is within magnitudes of up to Mw 5.1, consisting of: one significant mainshock of Mw 5.1 in 2004 at the Dead Sea basin, two swarms at the Sea of Galilee of Mw ∼3.5 and ∼4.5, in 2013 and 2018, respectively, and sporadic seismicity of up to Mw 4.4, throughout the DST. Except for the swarms, which show a clear decrease of the b-values (Figure 8), all the main shocks, including the Mw 5.1, are characterised by a relatively short aftershock activity (Hofstetter et al., 2008); for the Mw 5.1 event, the complete aftershock sequence includes 23 events above the Mc of 2.1. The b-value spatial peak, corresponding to this sequence in the spatial b-value profile (Lat ∼31.7; Figure 7), is not aligned with the spatial seismogenic depth profile and the spatial seismic moment release profiles. Therefore, it seems that the temporal variation of the b-value within the DST mainly reflects the effect of moderate seismicity (Mw 3–5) on the whole size distribution, a phenomena that has been observed particularly in small datasets (Marzocchi et al., 2020), and in this case - with low aftershock activity. For the swarms (Figure 8), the reduction in the b-value reflects the repetition of similar size mainshocks, typical for swarms (e.g., Goebel et al., 2016, Wetzler et al., 2019).
6.3 Seismic gaps along the DST
Some of the sections along the DST, which show extremely low seismic activity, were documented by Sharon et al. (2020). The most prominent ones (Figure 7) are the Hazeva Seismic Gap (HSG) and the Beit She’an Seismic Gap (BSG). The HSG contains a “quasi-spike” b-value rise, which correlates to a rather similar “quasi-spike” of low seismic moment release (Figure 7). Since high b-values also relate to weak zones that may be associated with creeping (Schorlemmer et al., 2004; Wyss et al., 2004; Tormann et al., 2013, 2014), the combination of a high b-value with low seismic activity implies that local creeping occurs in this zone. Geodetic data from the BSG has been interpreted to show shallow crustal creep (1.5 ± 1.0 km) in its northernmost part (Hamiel et al., 2016). Sharon et al. (2020) suggested that the relative absence of seismicity in this gap is caused by slip partitioning of the DST activity to the CTF (Sadeh et al., 2012; Hamiel et al., 2018), and to some extent, shallow creeping. Assuming the low b-value in the southern part of the BSG is not an artefact of the small amount of data, it may indicate a particularly locked section (e.g., asperity; Wiemer and Wyss, 1997, 2002; Wyss, 2001; Spada et al., 2013; Tormann et al., 2014). The sparse seismicity observed here and the relative quiescence along the segment in the past centuries, as well as the associated large displacement deficit (Marco and Klinger, 2014; Lefevre et al., 2018), should be accounted for in seismic hazard assessment.
6.4 Comparison to other continental transform plate boundaries
Comparison of b-values between different tectonic domains may be problematic, since some factors can dramatically affect the b-value: the magnitude type (Wyss, 2020) and its calibration (Tormann et al., 2010), changes in the network operating procedures (Zúñiga and Wiemer, 1999), the spatial frame (e. g., Page and Felzer, 2015), and the time window in respect of the seismic cycle (e.g., Raub et al., 2017); that is, assuming a correct choice of the completeness magnitude and well-calculated b-value (i.e., the method used). Thus, although the b-value is expected to vary between regions according to the faulting style (Schorlemmer et al., 2005), such interpretations between different faults of different regions should be made carefully. Here, we compare our results, regarding the DST, with three other major continental transform fault systems (e.g., Gauriau and Dolan, 2021).
A b-value of 1.03 was deduced for a 20-km width spatial zone surrounding the southern San Andreas fault (Page and Felzer, 2015). Rather similar b-values were achieved in a wider zone (0.99–1.01; Hutton et al., 2010). Two zones, comprised of widths of 10 and 20 km off the northern part of the San Andreas fault, California, yielded b-values of 0.99 and 0.93, respectively (Wyss, 2020). In New Zealand, a b-value of 0.98 was obtained for a 20-km wide zone along the central Alpine fault (Wyss, 2020), which accommodates reverse slip rates of approximately 20–40% of its horizontal motion (Norris and Cooper, 2001). A larger zone surrounding the central Alpine fault deduced a b-value of 0.85 (Michailos et al., 2019). Seismicity from the area of the North Anatolian fault, Turkey, divided into three parts, yielded b-values of 1.01–1.02 in its eastern and central parts, and 1.13 in the western part (Öztürk, 2011).
Our results indicate a b-value of 0.93 for the DST with a gentle spatial trend associated with tensional and compressional components in the edges of the study area. Considering that the northern part of the San Andreas fault is somewhat affected by a compression regime (e.g., Schwartz et al., 1990; Williams et al., 2006), and that the western part of the North Anatolian fault accommodates normal faulting components (Reilinger et al., 1997), their corresponding b-values imply also a gentle faulting-style dependency of the b-value, similar to our DST observations. In addition, the b-value of the DST is somewhat lower compared with these other major strike-slip faults, despite reverse component in parts of the San Andreas fault and particularly in the Alpine fault. This might be due to relatively deeper seismicity of the DST (up to 21–27 km), in comparison to the other examined faults, possibly associated with increased differential stress along significant parts of the DST.
7 Conclusions
We analyse ∼20,300 earthquakes along the DST and its periphery and provide a new conversion formula from to . After converting magnitudes, we determine the frequency-magnitude parameters for the complete catalogue, and to five significant subsets. Our b-value results align with previous studies in the region (Salamon et al., 1996; Shapira and Hofstetter, 2002; Hamiel et al., 2009; Lu et al., 2020), providing a high resolution characterisation of the frequency-magnitude statistical parameters, including within well-defined tectonic zones, consisting the network coverage area.
Variations of the b-value within the tectonic zones correspond to changes in the tectonic regime and presumably to the associated style of faulting, with values that vary from 0.93 in the strike-slip dominated DST to 1.19 in the normal faulting dominated OFZ(N)). In addition, high-resolution spatial profile along the DST reveals a decreasing trend of the b-value towards the north. This trend corresponds to a possible increased differential stress associated with a gradual change in the faulting regime, from an extension component at the south, where the DST shifts from the extensional pull-apart basins of the Gulf of Elat, closer to the Red Sea rift, to a compression component at the north, approaching the LRB. Hence, we propose that these variations reflect the inverse dependency of the b-value with the differential stress as revealed from the theory of faulting.
The DST spatial profile reveals an inverse dependency between the b-value and the seismogenic depth as the b-value decreases with the deepening of the seismogenic depth. Considering that deeper seismogenic zones are distributed over depths of higher differential stress, we suggest that this inverse correlation results from variations in the differential stress, which increases with depth within the brittle part of the crust, showing rheological dependency of the b-value.
Our analyses show that within the DST zone, the b-value inversely correlates with seismic moment release, reflecting a weak role of aftershocks within the DST. While the temporal fluctuations of the b-value better reflect this correlation, its spatial variations are linked to both the seismic moment release and the seismogenic depth.
The observations from this study are in line with previously-suggested dependency between the b-value and the differential stress (Scholz, 1968, 2015; ; Goebel et al., 2013). Considering this dependency and associated previous observations, we also suggest that anomalous b-values may indicate creeping and locked sections of the fault in two seismic gaps defined here—the Hazeva and Beit She’an seismic gaps, respectively. Our results contribute to the seismo-tectonic understanding of the DST, and should also be considered for seismic hazard evaluation in the region. In addition, as the observed spatio-temporal relations between the b-value and the seismic moment release are not well understood, this relation should be further investigated in higher resolution, and when more detailed data will be available for the region.
7.1 Data and resources
The seismic catalogue used in this study is an updated relocated catalogue, based on what was presented in Wetzler and Kurzon, 2016, and can be downloaded at: <https://www.gov.il/en/Departments/General/seismic-catalogs-files>. Maps were generated using QGIS open source mapping software [(<https://www.qgis.org/>), last accessed November 2022], and the global high resolution population GIS layers were downloaded from <https://www.worldpop.org/>. Many of the figures were generated by Matlab (<https://www.mathworks.com/products/matlab.html>), last accessed November 2022, and all figures were finalized using Microsoft PowerPoint. Naturally, we cannot cover all aspects in the main manuscript, and accordingly, we have added a Supplemental Material section, providing more details and additional supporting analyses.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://www.gov.il/en/Departments/General/seismic-catalogs-files.
Author contributions
MS and IK contributed to the conception and the design of the study. MS and IK performed the majority of the data analysis. MS and NW performed additional data analysis. MS wrote the first draft of the manuscript. IK and NW wrote sections of the manuscript. AS and SM provided essential geological and geophysical insights. All authors contributed to manuscript revision, read, and approved the submitted version.
Funding
This study was supported by the Ministry of Energy grants: 214-11-10, 219-11-054 and 222-11-001.
Acknowledgments
We thank Andrey Polozov, Ran Nof, Marina Gorstein, Yuval Tal, Yariv Hamiel and Marcelo Rosensaft for their assistance and their useful comments. We thank the guest associate editor, Guido Adinolfi, and two reviewers for their helpful comments, assisting us in improving the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2022.1074729/full#supplementary-material
References
1
AchmonM. (1986). The Carmel border fault between yoqneam and nesher. M.Sc. thesis. Jerusalem: Hebrew University of Jerusalem.
2
AgnonA. (2014). “Pre-instrumental earthquakes along the Dead Sea rift,”in Dead Sea transform fault system: Reviews. Editors GarfunkelZ.Ben-AvrahamZ.KaganE. (Dordrecht: Springer Netherlands), 207–262.
3
AkiK. (1965). Maximum likelihood estimate of b in the formula log N = a - bM and its confidence limits. Bull. Earthq. Res. Inst. Tokyo Univ.43, 237–239.
4
AkiK.RichardsP. G. (2002). Quantitative seismology. Sausalito, CA: University Science Books.
5
AkaikeH. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control19 (6), 716–723. 10.1109/TAC.1974.1100705
6
AldersonsF.Ben-AvrahamZ.HofstetterA.KisslingE.Al-YazjeenT. (2003). Lower-crustal strength under the Dead Sea basin from local earthquake data and rheological modeling. Earth Planet. Sci. Lett.214 (1–2), 129–142. 10.1016/S0012-821X(03)00381-9
7
AldersonsF.Ben-AvrahamZ. (2014). “The seismogenic thickness in the Dead Sea area,”in Dead Sea transform fault system: Reviews. Editors GarfunkelZ.Ben-AvrahamZ.KaganE. J. (Dordrecht: Springer Netherlands), 53–89. 10.1007/978-94-017-8872-4
8
AmbraseysN. (2009). Earthquakes in the mediterranean and Middle East. Cambridge: Cambridge University Press. 10.1017/CBO9781139195430
9
AmitranoD. (2003). Brittle-ductile transition and associated seismicity: Experimental and numerical studies and relationship with the b value. J. Geophys. Res.108 (B1), 1–15. 10.1029/2001jb000680
10
AndersonE. M. (1905). The dynamics of faulting.Trans. Edinb. Geol. Soc.8 (3), 387–402. 10.1144/transed.8.3.387
11
AshishC.LindholmI. A.ParvezKühnD.KuhnD. (2016). Probabilistic earthquake hazard assessment for Peninsular India. J. Seismol.20 (2), 629–653. 10.1007/s10950-015-9548-2
12
AtaevaG.ShapiraA.HofstetterA. (2015). Determination of source parameters for local and regional earthquakes in Israel. J. Seismol.19 (2), 389–401. 10.1007/s10950-014-9472-x
13
AvniY.BartovY.GarfunkelZ.GinatH. (2000). Evolution of the Paran drainage basin and its relation to the Plio-Pleistocene history of the Arava Rift Western margin, Israel. Israel J. Earth Sci.494, 215–238. 10.1560/W8WL-JU3Y-KM7W-8LX4
14
AvniY.BartovY.GarfunkelZ.GinatH. (2001). The arava formation-A pliocene sequence in the arava valley and its Western margin, southern Israel. Isr. J. Earth Sci.50 (2–4), 101–120. 10.1092/5U6A-RM5E-M8E3-QXM7
15
BalaA.RadulianM.PopescuE. (2003). Earthquakes distribution and their focal mechanism in correlation with the active tectonic zones of Romania. J. Geodyn.36 (1–2), 129–145. 10.1016/S0264-3707(03)00044-9
16
BartovY. (1974). A Structural and paleogeographical study of the central Sinai faults and domes. Hebrew, Englilsh abstract. PhD thesis (Jerusalem: Hebrew University of Jerusalem).
17
BartovY.AvniY.CalvoR.FrieslanderU. (1998). The Zofar fault - a major intra-rift feature in the Arava rift valley,. Geol. Soc. Isr. Curr. Res.11, 27–32.
18
BeallA.van den EndeM.AmpueroJ. P.CapitanioF. A.Åfagereng (2022). Linking earthquake magnitude-frequency statistics and stress in visco-frictional fault zone models. Geophys. Res. Lett.49 (20), e2022GL099247. 10.1029/2022gl099247
19
Ben-AvrahamZ.GinzburgA. (1990). Displaced terranes and crustal evolution of the Levant and the eastern Mediterranean. Tectonics9 (4), 613–622. 10.1029/TC009i004p00613
20
Ben-MenahemA. (1991). Four thousand years of seismicity along the Dead Sea rift. J. Geophys. Res.96 (B12), 20195–20216. 10.1029/91jb01936
21
Ben-MenahemA. (1981). Variation of slip and creep along the levant rift over the past 4500 years. Tectonophysics80 (1–4), 183–197. 10.1016/0040-1951(81)90149-9
22
BentorY. K.VromanA. (1954). A Structural contour map of Israel (1:250, 000) with remarks on its dynamic interpretation. Bull. Res. Counc. Isr.4 (2), 125–135.
23
BenderB. (1983). Maximum likelihood estimation of b values for magnitude grouped data. Bull. Seismol. Soc. Am. 73 (3), 831–851.
24
BoraD. K.BorahK.MahantaR.MadhabJ. (2018). Seismic b - values and its correlation with seismic moment and Bouguer gravity anomaly over Indo-Burma ranges of northeast India : Tectonic implications. Tectonophysics728–729, 130–141. 10.1016/j.tecto.2018.01.001
25
BormannP.Di GiacomoD. (2011). The moment magnitude Mw and the energy magnitude Me: Common roots and differences,. J. Seismol.15 (2), 411–427. 10.1007/s10950-010-9219-2
26
BusZ.GrenerczyG.TóthL.MónusP. (2009). Active crustal deformation in two seismogenic zones of the Pannonian region - GPS versus seismological observations. Tectonophysics474 (1–2), 343–352. 10.1016/j.tecto.2009.02.045
27
BustosJ. P. A.TaroniM.AdamL. (2022). A Robust Statistical Framework to Properly Test the Spatiotemporal Variations of the b-Value: An Application to the Geothermal and Volcanic Zones of the Nevado del Ruiz Volcano. Seismol. Res. Lett.93, 2793–2803. 10.1785/0220220004
28
ByerleeJ. D. (1968). Brittle-ductile transition in rocks. J. Geophys. Res.73 (14), 4741–4750. 10.1029/JB073i014p04741
29
CalvoR.BartovY. (2001). Hazeva Group, southern Israel: New observations, and their implications for its stratigraphy, paleogeography, and tectono-sedimentary regime. Isr. J. Earth Sci.50 (2–4), 71–100. 10.1092/B02L-6K04-UFQL-KUE3
30
CaoA.GaoS. S. (2002). Temporal variation of seismic b-value beneath northeastern Japan Island Arc. Geophys. Res. Lett.29 (9), 48-1–48-3. 10.1029/2001GL013775
31
CastellaroS.MulargiaF. S.KaganY. Y. (2006). Regression problems for magnitudes. Geophys. J. Int.165 (3), 913–930. 10.1111/j.1365-246X.2006.02955.x
32
CastelloB.OlivieriM.SelvaggiG. (2007). Local and duration magnitude determination for the Italian earthquake catalog, 1981-2002. Bull. Seismol. Soc. Am.97 (1), 128–139. 10.1785/0120050258
33
ChoyG. L.BoatwrightJ. L. (1995). Global patterns of radiated seismic energy and apparent stress. J. Geophys. Res.100 (B9), 18205–18228. 10.1029/95jb01969
34
ClausetA.ShaliziC. R.NewmanM. E. J. (2009). Power-law distributions in empirical data. SIAM Rev.51 (4), 661–703. 10.1137/070710111
35
DasR.WasonH. R.GonzalezG.SharmaM. L.ChoudhuryD.LindholmC.et al (2018). Earthquake magnitude conversion problem. Bull. Seismol. Soc. Am.108 (4), 1995–2007. 10.1785/0120170157
36
El-IsaZ. H.EatonD. W. (2014). Spatiotemporal variations in the b-value of earthquake magnitude-frequency distributions: Classification and causes. Tectonophysics615–616, 1–11. 10.1016/j.tecto.2013.12.001
37
FelzerK. R.AbercrombieR. E.EkströmG. (2004). A common origin for aftershocks, foreshocks, and multiplets. Bull. Seismol. Soc. Am.94 (1), 88–98. 10.1785/0120030069
38
FisherR. A. (1950). Contribution to Mathematical Statistics. New York, NY: Wiley.
39
FrankelA. (1995). Mapping seismic hazard in the central and eastern United States. Seismol. Res. Lett.66 (4), 8–21. 10.1785/gssrl.66.4.8
40
FreundR. (1970). The geometry of faulting in the Galilee. Isr. J. Earth Sci.19 (1943), 117–140.
41
FrohlichC.DavisS. D. (1993). Teleseismic b values; or, much ado about 1.0. J. Geophys. Res.98 (B1), 631–644. 10.1029/92jb01891
42
GarfunkelZ. (1981). Internal structure of the Dead Sea leaky transform (rift) in relation to plate kinematics. Tectonophysics80 (1–4), 81–108. 10.1016/0040-1951(81)90143-8
43
GarfunkelZ. (2014).“Lateral motion and deformation along the Dead Sea transform,” in Dead Sea transform fault system: Reviews. Editors GarfunkelZ.Ben-AvrahamZ.KaganE. (Dordrecht: Springer Netherlands), 109–150.
44
GarfunkelZ. (2010). The long- and short-term lateral slip and seismicity along the Dead Sea Transform: An interim evaluation. Israel J. Earth Sci.58 (3–4), 217–235. 10.1560/IJES.58.3-4.217
45
GarfunkelZ.BartovY. (1977). The tectonics of the Suez rift,. Geol. Surv. Isr. Bull.71, 1–44.
46
GarfunkelZ.Ben-AvrahamZ. (2001). Basins along the Dead Sea transform, mémoires du Muséum natl. d’histoire Nat.186, 607–627.
47
GarfunkelZ.Ben-AvrahamZ. (1996). The structure of the Dead Sea basin,. Tectonophysics266 (1-4), 155–176. 10.1016/s0040-1951(96)00188-6
48
GauriauJ.DolanJ. F. (2021). Relative structural complexity of plate-boundary fault systems controls incremental slip-rate behavior of major strike-slip faults. Geochem. Geophys. Geosyst.22, e2021GC009938. 10.1029/2021GC009938
49
GerstenbergerM.WiemerS.GiardiniD. (2001). A systematic test of the hypothesis that the b value varies with depth in California. Geophys. Res. Lett.28 (1), 57–60. 10.1029/2000GL012026
50
GibowiczS. J. (1974). Frequency-magnitude, depth, and time relations for earthquakes in an island arc: North Island, New Zealand. Tectonophysics23 (3), 283–297. 10.1016/0040-1951(74)90028-6
51
GinatH.ZilbermanE.AmitR. (2002). Red sedimentary units as indicators of Early Pleistocene tectonic activity in the southern Negev desert, Israel. Geomorphology45 (1–2), 127–146. 10.1016/S0169-555X(01)00193-3
52
GinatH.ZilbermanE.AvniY. (2000). Tectonic and paleogeographic significance of the edom river, a pliocene stream that crossed the Dead Sea rift valley,. Israel J. Earth Sci.49 (3), 159–177. 10.1560/N2P9-YBJ0-Q44Y-GWYN
53
GittermanY.PinskyV.ShapiraA.ErginM.KalafatD.GurbuzG.et al (2005). Improvement in detection, location, and identification of small events through joint data analysis by seismic stations in the Middle East/Eastern Mediterranean region. Final Report DTRA01-00-C-0119.
54
GlaisterP. (2005). The use of orthogonal distances in generating the total least squares estimate. Math. Comput. Educ.39 (1), 21–30.
55
GoebelT. H. W.HosseiniS. M.CappaF.HaukssonE.AmpueroJ. P.AminzadehF.et al (2016). Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California. Geophys. Res. Lett.43, 1092–1099. 10.1002/2015GL066948
56
GoebelT. H. W.KwiatekG.BeckerT. W.BrodskyE. E.DresenG. (2017). What allows seismic events to grow big? Insights from b-value and fault roughness analysis in laboratory stick-slip experiments. Geology45 (9), 815–818. 10.1130/G39147.1
57
GoebelT. H. W.SchorlemmerD.BeckerT. W.DresenG.SammisC. G. (2013). Acoustic emissions document stress changes over many seismic cycles in stick-slip experiments. Geophys. Res. Lett.40 (10), 2049–2054. 10.1002/grl.50507
58
GomezF.KaramG.KhawlieM.McCluskyS.VernantP.ReilingerR.et al (2007). Global Positioning System measurements of strain accumulation and slip transfer through the restraining bend along the Dead Sea fault system in Lebanon. Geophys. J. Int.168 (3), 1021–1028. 10.1111/j.1365-246X.2006.03328.x
59
GomezF.MeghraouiM.DarkalA. N.HijaziF.MoutyM.SuleimanY.et al (2003). Holocene faulting and earthquake recurrence along the Serghaya branch of the Dead Sea fault system in Syria and Lebanon. Geophys. J. Int.153 (3), 658–674. 10.1046/j.1365-246X.2003.01933.x
60
GossetW. S. (1908). The probable error of a mean,. Biometrika6 (1), 1–25. 10.2307/2331554
61
GuliaL.RinaldiA. P.TormannT.VannucciG.EnescuB.WiemerS. (2018). The effect of a mainshock on the size distribution of the aftershocks. Geophys. Res. Lett.45, 13,277–13,287. 10.1029/2018GL080619
62
GuliaL.WiemerS. (2010). The influence of tectonic regimes on the earthquake size distribution: A case study for Italy. Geophys. Res. Lett.37 (10), 1–6. 10.1029/2010GL043066
63
GutenbergB.RichterC. F. (1944). Frequency of earthquakes in California. Bull. Seismol. Soc. Am.34, 185–188. 10.1785/bssa0340040185
64
HamielY.AmitR.BeginZ. B.MarcoS.KatzO.SalamonA.et al (2009). The seismicity along the Dead Sea fault during the last 60, 000 years. Bull. Seismol. Soc. Am.99 (3), 2020–2026. 10.1785/0120080218
65
HamielY.MassonF.PiatibratovaO.MizrahiY. (2018). GPS measurements of crustal deformation across the southern Arava Valley section of the Dead Sea Fault and implications to regional seismic hazard assessment. Tectonophysics724–725, 171–178. 10.1016/j.tecto.2018.01.016
66
HamielY.PiatibratovaO.MizrahiY. (2016). Creep along the northern Jordan valley section of the Dead Sea fault. Geophys. Res. Lett.43 (6), 2494–2501. 10.1002/2016GL067913
67
HanksT. C.KanamoriH. (1979). A moment magnitude scale,. J. Geophys. Res.84 (B5), 2348–2350. 10.1029/JB084iB05p02348
68
HelmstetterA.KaganY. Y.JacksonD. D. (2006). Comparison of short-term and time-independent earthquake forecast models for southern California. Bull. Seismol. Soc. Am.96 (1), 90–106. 10.1785/0120050067
69
HendersonJ.MainI. (1992). A simple fracture-mechanical model for the evolution of seismicity. Geophys. Res. Lett.19 (4), 365–368. 10.1029/92gl00274
70
HofstetterA.ThioH. K.ShamirG. (2003). Source mechanism of the 22/11/1995 gulf of aqaba earthquake and its. J. Seismol.7, 99–114. 10.1023/a:1021206930730
71
HofstetterA.Van EckT.ShapiraA. (1996). Seismic activity along fault branches of the Dead Sea-Jordan transform system: The Carmel-Tirtza fault system. Tectonophysics267 (1–4), 317–330. 10.1016/S0040-1951(96)00108-4
72
HofstetterR.GittermanY.PinksyV.KraevaN.FeldmanL. (2008). Seismological observations of the northern Dead Sea basin earthquake on 11 February 2004 and its associated activity. Israel J. Earth Sci.57 (2), 101–124. 10.1560/IJES.57.2.101
73
HofstetterR.KlingerY.AmratA. Q.RiveraL.DorbathL. (2007). Stress tensor and focal mechanisms along the Dead Sea fault and related structural elements based on seismological data. Tectonophysics429 (3–4), 165–181. 10.1016/j.tecto.2006.03.010
74
HurvichC. M.TsaiC.-L. (1989). Regression and time series model selection in small samples. Biometrika76, 297–307. 10.1093/biomet/76.2.297
75
HuttonK.WoessnerJ.HaukssonE. (2010). Earthquake monitoring in southern California for seventy-seven years (1932–2008). Bull. Seismol. Soc. Am.100 (2), 423–446. 10.1785/0120090130
76
JoffeS.GarfunkelZ. (1987). Plate kinematics of the circum Red Sea-a re-evaluation,. Tectonophysics141 (1–3), 5–22. 10.1016/0040-1951(87)90171-5
77
KadirioğluF. T.KartalR. F. (2016). The new empirical magnitude conversion relations using an improved earthquake catalogue for Turkey and its near vicinity (1900–2012). Turk. J. Earth Sci.25 (4), 300–310. 10.3906/yer-1511-7
78
KanamoriH. (1983). Magnitude scale and quantification of earthquakes. Tectonophysics93 (3–4), 185–199. 10.1016/0040-1951(83)90273-1
79
KanamoriH. (1977). The energy release in great earthquakes. J. Geophys. Res.82 (20), 2981–2987. 10.1029/jb082i020p02981
80
KaneM. T.MrochA. A. (2020). Orthogonal regression, the cleary criterion, and lord's paradox: Asking the right questions. ETS Res. Rep. Ser.2020, 1–24. 10.1002/ets2.12298
81
KingG. (1983). The accommodation of large strains in the upper lithosphere of the Earth and other solids by self-similar fault systems: The geometrical origin of b-value. Pure Appl. Geophys.121 (5–6), 761–815. 10.1007/BF02590182
82
KumarR.YadavR. B. S.CastellaroS. (2020). Regional earthquake magnitude conversion relations for the himalayan seismic belt. Seismol. Res. Lett.91 (6), 3195–3207. 10.1785/0220200204
83
KurzonI.NofR. N.LaporteM.LutzkyH.PolozovA.ZakoskyD.et al (2020). The “TRUAA” seismic network: Upgrading the Israel seismic network—toward national earthquake early warning system. Seismol. Res. Lett.91 (6), 3236–3255. 10.1785/0220200169
84
LefevreM.KlingerY.Al-QaryoutiM.Le BéonM.MoumaniK. (2018). Slip deficit and temporal clustering along the Dead Sea fault from paleoseismological investigations. Sci. Rep.8 (1), 4511–4519. 10.1038/s41598-018-22627-9
85
LeiX.KusunoseK.RaoM. V. M. S.NishizawaO.SatohT. (2000). Quasi-static fault growth and cracking in homogeneous brittle rock under triaxial compression using acoustic emission monitoring. J. Geophys. Res.105 (3), 6127–6139. 10.1029/1999jb900385
86
Llana-FúnezS.Lõpez-FernándezC. (2015). The seismogenic zone of the continental crust in Northwest Iberia and its relation to crustal structure. Tectonics34 (8), 1751–1767. 10.1002/2015TC003877
87
LuY.WetzlerN.WaldmannN.AgnonA.BiasiG. P.MarcoS. (2020). A 220, 000-year-long continuous large earthquake record on a slow-slipping plate boundary. Sci. Adv.6 (48), eaba4170. 10.1126/sciadv.aba4170
88
MagrinA.PeresanA.KronrodT.VaccariF.PanzaG. F. (2017). Neo-deterministic seismic hazard assessment and earthquake occurrence rate,. Eng. Geol.229, 95–109. 10.1016/j.enggeo.2017.09.004
89
MainI. G.MeredithP. G.JonesC. (1989). A reinterpretation of the precursory seismic b-value anomaly from fracture mechanics. Geophys. J. Int.96 (1), 131–138. 10.1111/j.1365-246X.1989.tb05255.x
90
MaitiS. K.KamaiR. (2020). Interaction between fault and off-fault seismic sources in hazard analysis – a case study from Israel. Eng. Geol.274, 105723. 10.1016/j.enggeo.2020.105723
91
MandalP.SrinageshD.SureshG.NareshB.NaiduM.SinghD. K.et al (2021). Characterization of earthquake hazard at the Palghar and Pulichintala swarm activity regions (India) through three-dimensional modelling of b-value and fractal (correlation) dimensions. Nat. Hazards108, 1183–1196. 10.1007/s11069-021-04726-5
92
MarcoS. (2007). Temporal variation in the geometry of a strike-slip fault zone: Examples from the Dead Sea Transform. Tectonophysics445 (3–4), 186–199. 10.1016/j.tecto.2007.08.014
93
MarcoS.KlingerY. (2014). “Review of on-fault palaeoseismic studies along the Dead Sea fault,”in Dead Sea transform fault system: Reviews. Editors GarfunkelZ.Ben-AvrahamZ.KaganE. (Dordrecht: Springer Netherlands), 183–205. 10.1007/978-94-017-8872-4
94
MarzocchiW.TaroniM. (2014). Some thoughts on declustering in probabilistic seismic-hazard analysis. Bull. Seismol. Soc. Am.104 (4), 1838–1845. 10.1785/0120130300
95
MarzocchiW.SandriL. (2003). A review and new insights on the estimation of the b value and its uncertainty. Ann. Geophys. 46 (6), 1271–1282.
96
MarzocchiW.SpassianiI.StalloneA.TaroniM. (2020). How to be fooled searching for significant variations of the b-value. Geophys. J. Int.220 (3), 1845–1856. 10.1093/gji/ggz541
97
MichailosK.SmithE. G. C.ChamberlainC. J.SavageM. K.TownendJ. (2019). Variations in seismogenic thickness along the central Alpine Fault, New Zealand, revealed by a decade’s relocated microseismicity. Geochem. Geophys. Geosyst.20, 470–486. 10.1029/2018GC007743
98
MizrahiL.NandanS.WiemerS. (2021). The effect of declustering on the size distribution of mainshocks. Seismol. Res. Lett.92 (4), 2333–2342. 10.1785/0220200231
99
MoriJ.AbercrombieR. E. (1997). Depth dependence of earthquake frequency-magnitude distributions in California: Implications for rupture initiation. J. Geophys. Res.102 (B7), 15081–15090. 10.1029/97jb01356
100
MousaviS. M.OgwariP. O.HortonS. P.LangstonC. A. (2017). Spatio-temporal evolution of frequency-magnitude distribution and seismogenic index during initiation of induced seismicity at Guy-Greenbrier, Arkansas. Phys. Earth Planet. Interiors267, 53–66. 10.1016/j.pepi.2017.04.005
101
MukhopadhyayB.MogrenS.MukhopadhyayM.DasguptaS. (2013). Incipient status of dyke intrusion in top crust - evidences from the Al-Ays 2009 earthquake swarm, Harrat Lunayyir, SW Saudi Arabia. Geomat. Nat. Hazards Risk4 (1), 30–48. 10.1080/19475705.2012.663794
102
NeevD.AlmagorG.AradA.GinzburgA.HallJ. K. (1976). The geology of the southeastern Mediterranean Sea,. Geol. Surv. Isr. Bull.68, 1–51.
103
NorrisR. J.CooperA. F. (2001). Late Quaternary slip rates and slip partitioning on the Alpine Fault, New Zealand. J. Struct. Geol.23 (2), 507–520. 10.1016/s0191-8141(00)00122-x
104
NuanninP.KulhanekO.PerssonL. (2005). Spatial and temporal b value anomalies preceding the devastating off coast of NW Sumatra earthquake of December 26, 2004. Geophys. Res. Lett.32 (11), 1–4. 10.1029/2005GL022679
105
NurielP.WeinbergerR.Kylander-ClarkA. R. C.HackerB. R.CraddockJ. P. (2017). The onset of the Dead Sea transform based on calcite age-strain analyses,. Geology45 (7), 587–590. 10.1130/G38903.1
106
ÖztürkS. (2011). Characteristics of seismic activity in the western, central and eastern parts of the North Anatolian fault zone, Turkey: Temporal and spatial analysis. Acta Geophys.59 (2), 209–238. 10.2478/s11600-010-0050-5
107
PageM.FelzerK. (2015). Southern San Andreas fault seismicity is consistent with the Gutenberg–Richter magnitude–frequency distribution. Bull. Seismol. Soc. Am.105 (4), 2070–2080. 10.1785/0120140340
108
PavlisG. L. (1986). Appraising earthquake hypocenter location errors: A complete, practical approach for single-event locations. Bull. Seismol. Soc. Am.76 (6), 1699–1717.
109
PavlisG. L.VernonF. L.HarveyD.QuinlanD. (2004). The generalized earthquake-location (GENLOC) package: An earthquake-location library. Comput. Geosci.30 (9/10), 1079–1091. 10.1016/j.cageo.2004.06.010
110
PetersenM. D.DawsonT. E.ChenR.CaoT.WillsC. J.SchwartzD. P.et al (2011). Fault displacement hazard for strike-slip faults. Bull. Seismol. Soc. Am.101 (2), 805–825. 10.1785/0120100035
111
PetruccelliA.SchorlemmerD.TormannT.RinaldiA. P.WiemerS.GasperiniP.et al (2019b). The influence of faulting style on the size-distribution of global earthquakes,. Earth Planet. Sci. Lett.527, 115791. 10.1016/j.epsl.2019.115791
112
PetruccelliA.GasperiniP.TormannT.SchorlemmerD.RinaldiA. P.VannucciG.et al (2019a). Simultaneous dependence of the earthquake-size distribution on faulting style and depth,. Geophys. Res. Lett.46 (20), 11044–11053. 10.1029/2019GL083997
113
PetruccelliA.VannucciG.LolliB.GasperiniP. (2018). Harmonic fluctuation of the slope of the frequency–magnitude distribution (b-value) as a function of the angle of rake. Bull. Seismol. Soc. Am.108 (4), 1864–1876. 10.1785/0120170328
114
QuennellA. M. (1959). “Tectonics of the Dead Sea rift,” in Proc. 20th int. Geol. Congr. Mex., 385–403.
115
RadulianM.BălăA.PopescuE.Toma-DănilăD. (2018). Earthquake mechanism and characterization of seismogenic zones in south-eastern part of Romania. Ann. Geophys.61 (1). 10.4401/ag-7443
116
RadulianM.MândrescuN.PanzaG. F.PopescuE.UtaleA. (2000). Characterization of seismogenic zones of Romania. Pure Appl. Geophys.157 (1–2), 57–77. 10.1007/978-3-0348-8415-0_4
117
RaubC.Martínez-GarzónP.KwiatekG.BohnhoffM.DresenG. (2017). Variations of seismic b-value at different stages of the seismic cycle along the North Anatolian Fault Zone in northwestern Turkey. Tectonophysics712, 232–248. 10.1016/j.tecto.2017.05.028
118
ReilingerR.McCluskyS.VernantP.LawrenceS.ErgintavS.CakmakR.et al (2006). GPS constraints on continental deformation in the Africa-Arabia-Eurasia continental collision zone and implications for the dynamics of plate interactions. J. Geophys. Res.1115, 1–26. 10.1029/2005JB004051
119
ReilingerR.E.McCluskyS. C.OralM. B.KingR. W.ToksozM. N.BarkaA. A.et al (1997). Global Positioning System measurements of present-day crustal movements in the Arabia-Africa-Eurasia plate collision zone. J. Geophys. Res. Solid Earth102 (B5), 9983–9999.
120
RiceJohn A. (2006). Mathematical statistics and data analysis. 3rd ed. Duxbury Advanced.
121
RigoA.SouriauA.SylvanderM. (2018). Spatial variations of b-value and crustal stress in the Pyrenees. J. Seismol.22, 337–352. 10.1007/s10950-017-9709-6
122
RivièreJ.LvZ.JohnsonP. A.MaroneC. (2018). Evolution of b-value during the seismic cycle: Insights from laboratory experiments on simulated faults. Earth Planet. Sci. Lett.482, 407–413. 10.1016/j.epsl.2017.11.036
123
RonH.EyalY. (1985). Intraplate deformation by block rotation and mesostructures along the Dead Sea transform, northern Israel. Tectonics4 (1), 85–105. 10.1029/tc004i001p00085
124
RonH.FreundR.GarfunkelZ.NurA. (1984). Block rotation by strike-slip faulting: Structural and paleomagnetic evidence. J. Geophys. Res.89B7, 6256–6270. 10.1029/JB089iB07p06256
125
RossZ. E.Ben-ZionY.WhiteM. C.VernonF. L. (2016). Analysis of earthquake body wave spectra for potency and magnitude values: Implications for magnitude scaling relations. Geophys. J. Int.207 (2), 1158–1164. 10.1093/gji/ggw327
126
RotsteinY.BrunerI.KafriU. (1993). High-resolution seismic imaging of the Carmel Fault and its implications for the structure of Mt Carmel. Isr. J. Earth Sci.42 (2), 55–69.
127
SadehM.HamielY.ZivA.BockY.FangP.WdowinskiS. (2012). Crustal deformation along the Dead Sea transform and the Carmel fault inferred from 12 years of GPS measurements. J. Geophys. Res.117 (8), 1–14. 10.1029/2012JB009241
128
SalamonA.HofstetterA.GarfunkelZ.RonH. (1996). Seismicity of the eastern mediterranean region: Perspective from the Sinai subplate. Tectonophysics263 (1–4), 293–305. 10.1016/S0040-1951(96)00030-3
129
ScholzC. H. (2015). On the stress dependence of the earthquake b value. Geophys. Res. Lett.42 (5), 1399–1402. 10.1002/2014GL062863
130
ScholzC. H. (1968). The frequency - magnitude relation of microfracturing in rock and its relation to earthquakes. Bull. Seismol. Soc. Am.58 (1), 399–415. 10.1785/bssa0580010399
131
SchorlemmerD.WiemerS.KleiberP.MaunderM. N.HarleyS. J. (2005). Microseismicity data forecast rupture area. Nature434 (7037), E1–E2. 10.1038/nature03581
132
SchorlemmerD.WiemerS.WyssM. (2004). Earthquake statistics at parkfield: 1. Stationarity of b values. J. Geophys. Res.109 (12), B12307–B12317. 10.1029/2004JB003234
133
SchorlemmerD.WiemerS.WyssM. (2005). Variations in earthquake-size distribution across different stress regimes. Nature437 (7058), 539–542. 10.1038/nature04094
134
SchwartzD. P.CoppersmithK. J. (1984). Fault behavior and characteristic earthquakes: Examples from the wasatch and san Andreas Fault Zones. J. Geophys. Res.89 (B7), 5681–5698. 10.1029/JB089iB07p05681
135
SchwartzS. Y.OrangeD. L.AndersonR. S. (1990). Complex fault interactions in a restraining bend on the San Andreas fault, southern Santa Cruz Mountains, California. Geophys. Res. Lett.17 (8), 1207–1210. 10.1029/gl017i008p01207
136
ScordilisE. M. (2006). Empirical global relations converting MS and mb to moment magnitude. J. Seismol.10 (2), 225–236. 10.1007/s10950-006-9012-4
137
ShalevE.LyakhovskyV.WeinsteinY.Ben-AvrahamZ. (2013). The thermal structure of Israel and the Dead Sea fault,. Tectonophysics602, 69–77. 10.1016/j.tecto.2012.09.011
138
ShalevE.LyakhovskyV.YechieliY. (2007). Is advective heat transport significant at the Dead Sea basin?Geofluids7 (3), 292–300. 10.1111/j.1468-8123.2007.00190.x
139
ShamirG.BartovY.SnehA.FleisherL.AradV.RosensaftM. (2001). Preliminary seismic zonation in Israel.GII report No. 550/95/01(1) and report No. GSI12/2000, 01, 28. no. 550.
140
ShapiraA. (1992). Detectability of regional seismic network: Analysis of the Israel seismic networks. Isr. J. Earth-Sciences41 (1), 21–25.
141
ShapiraA.AvniR.NorA. (1993). A new estimate for the epicenter of the Jericho earthquake of 11 July 1927. Isr. J. Earth-Sciences42 (2), 93–96.
142
ShapiraA.HofstetterA. (2002). Seismicity parameters of seismogenic zones. RELEMR-MERC Rep. 2001, no. 1980. Available at: http://www.relemr–merc.org/Other_information/seis.
143
SharonM. (2020). “Mapping and characterising active tectonic sources in Israel and adjacent areas,” in Geological survey of Israel report No. GSI/19/20M.Sc. thesis. Jerusalem (Israel): Tel Aviv University.
144
SharonM.SagyA.KurzonI.MarcoS.RosensaftM. (2020). Assessment of seismic sources and capable faults through hierarchic tectonic criteria: Implications for seismic hazard in the Levant. Nat. Hazards Earth Syst. Sci.20 (1), 125–148. 10.5194/nhess-20-125-2020
145
ShearerP. M. (2009). Introduction to seismology. New York: Cambridge University Press.
146
ShiY.BoltB. A. (1982). The standard error of the magnitude-frequency b value,. Bull. Seismol. Soc. Am.72 (5), 1677–1687. 10.1785/bssa0720051677
147
SibsonR. H. (1974). Frictional constraints on thrust, wrench and normal faults. Nature249 (5457), 542–544. 10.1038/249542a0
148
SibsonR. H. (1984). Roughness at the base of the seismogenic zone: Contributing factors. J. Geophys. Res.89 (B7), 5791–5799. 10.1029/JB089iB07p05791
149
SokolovV.ZahranH. M.YoussefS. E. H.El-HadidyM.AlraddadiW. W. (2017). Probabilistic seismic hazard assessment for Saudi Arabia using spatially smoothed seismicity and analysis of hazard uncertainty. Bull. Earthq. Eng.15 (7), 2695–2735. 10.1007/s10518-016-0075-5
150
SpadaM.TormannT.WiemerS.EnescuB. (2013). Generic dependence of the frequency-size distribution of earthquakes on depth and its relation to the strength profile of the crust. Geophys. Res. Lett.40 (4), 709–714. 10.1029/2012GL054198
151
StecklerM. S.BerthelotF.LyberisN.Le PichonX. (1988). Subsidence in the gulf of suez: Implications for rifting and plate kinematics. Tectonophysics153 (1–4), 249–270. 10.1016/0040-1951(88)90019-4
152
SuyehiroB. Y. S. (1966). Difference between aftershocks and foreshocks in the relationship of magnitude to frequency of occurrence for the great Chilean earthquake of 1960. Bull. Seismol. Soc. Am.56 (1), 185–200. 10.1785/bssa0560010185
153
TanJ. Y.WaldhauserF.TolstoyM.WilcockW. S. D. (2019). Axial Seamount: Periodic tidal loading reveals stress dependence of the earthquake size distribution (b value). Earth Planet. Sci. Lett.512, 39–45. 10.1016/j.epsl.2019.01.047
154
TormannT.WiemerS.HaukssonE. (2010). Changes of reporting rates in the southern California earthquake catalog, introduced by a new definition of ML. Bull. Seismol. Soc. Am.100 (4), 1733–1742. 10.1785/0120090124
155
TormannT.WiemerS.MetzgerS.MichaelA.HardebeckJ. L. (2013). Size distribution of parkfield’s microearthquakes reflects changes in surface creep rate. Geophys. J. Int.193 (3), 1474–1478. 10.1093/gji/ggt093
156
TormannT.WiemerS.MignanA. (2014). Systematic survey of high-resolution b value imaging along Californian faults : Inference on asperities. J. Geophys. Res. Solid Earth119, 2029–2054. 10.1002/2013JB010867
157
UtsuT. (2002). 44 Relationships between magnitude scales. Int. Geophys.81, 733–746. 10.1016/S0074-6142(02)80247-9
158
UtsuT. (1999). “Representation and analysis of the earthquake size distribution: A historical review and some new approaches,” in Seismicity patterns, their statistical significance and physical meaning (Basel: Birkhäuser), 509–535.
159
UtsuT. (1966). A statistical significance test of the difference in b-value between two earthquake groups. J. Phys. Earth14 (2), 37–40.
160
WeinbergerR.GrossM. R.SnehA. (2009). Evolving deformation along a transform plate boundary: Example from the Dead Sea Fault in northern Israel. Tectonics28 (5), 1–19. 10.1029/2008TC002316
161
WeinbergerR.NurielP.Kylander-clarkA. R. C.CraddockJ. P. (2020). Temporal and spatial relations between large-scale fault systems: Evidence from the Sinai-Negev shear zone and the Dead Sea Fault. Earth. Sci. Rev.211, 103377. 10.1016/j.earscirev.2020.103377
162
WesnouskyS. G.ScholzC. H.ShimazakiK.MatsudaT. (1983). Earthquake frequency distribution and the mechanics of faulting. J. Geophys. Res.B11, 9331–9340. 10.1029/jb088ib11p09331
163
WetzlerN.KurzonI. (2016). The earthquake activity of Israel: Revisiting 30 years of local and regional seismic records along the Dead Sea transform. Seismol. Res. Lett.87 (1), 47–58. 10.1785/0220150157
164
WetzlerN.ShalevE.GoebelT.AmelungF.KurzonI.LyakhovskyV.et al (2019). Earthquake swarms triggered by groundwater extraction near the Dead Sea Fault. Geophys. Res. Lett.46 (14), 8056–8063. 10.1029/2019GL083491
165
WiemerS.KatsumataK. (1999). Spatial variability of seismicity parameters in aftershock zones. J. Geophys. Res.104 (B6), 13135–13151. 10.1029/1999jb900032
166
WiemerS.McNuttS. R.WyssM. (1998). Temporal and three-dimensional spatial analyses of the frequency-magnitude distribution near Long Valley Caldera, California. Geophys. J. Int.134 (2), 409–421. 10.1046/j.1365-246X.1998.00561.x
167
WiemerS.WyssM. (2002). Mapping spatial variability of the frequency-magnitude distribution of earthquakes. Adv. Geophys.45, 259. 10.1016/S0065-2687(02)80007-3
168
WiemerS.WyssM. (1997). Mapping the frequency-magnitude distribution in asperities: An improved technique to calculate recurrence times?J. Geophys. Res.102 (B7), 15115–15128. 10.1029/97jb00726
169
WilliamsT. B.KelseyH. M.FreymuellerJ. T. (2006). GPS-derived strain in northwestern California: Termination of the San Andreas fault system and convergence of the Sierra Nevada–Great Valley block contribute to southern Cascadia forearc contraction. Tectonophysics413 ((3–4)), 171–184. 10.1016/j.tecto.2005.10.047
170
WyssM. (2001). Locked and creeping patches along the Hayward fault, California. Geophys. Res. Lett.28 (18), 3537–3540. 10.1029/2001gl013499
171
WyssM. (2020). Return times of large earthquakes cannot Be estimated correctly from seismicity rates: 1906 san francisco and 1717 Alpine fault ruptures. Seismol. Res. Lett.91 (4), 2163–2169. 10.1785/0220200008
172
WyssM. (1973). Towards a physical understanding of the earthquake frequency distribution, geophys. Geophys. J. Int.31 (4), 341–359. 10.1111/j.1365-246X.1973.tb06506.x
173
WyssM.SammisC. G.NadeauR. M.WiemerS. (2004). Fractal dimension and b-value on creeping and locked patches of the San Andreas fault near Parkfield, California. Bull. Seismol. Soc. Am.94 (2), 410–421. 10.1785/0120030054
174
XieW.HattoriK.HanP. (2019). Temporal variation and statistical assessment of the b value off the pacific coast of tokachi, hokkaido, Japan. Entropy21 (3), 249. 10.3390/e21030249
175
YadavR. B. S.BayrakY.TripathiJ. N.ChopraS.SinghA. P.BayrakE. (2012). A probabilistic assessment of earthquake hazard parameters in NW himalaya and the adjoining regions. Pure Appl. Geophys.169 (9), 1619–1639. 10.1007/s00024-011-0434-8
176
YadavR. B. S.BormannP.RastogiB. K.DasM. C.ChopraS. (2009). A Homogeneous and complete earthquake catalog for northeast India and the adjoining region. Seismol. Res. Lett.80 (4), 609–627. 10.1785/gssrl.80.4.609
177
Yagoda-BiranG.MaitiS. K.WetzlerN.NofR. N.PashcurY.KamaiR. (2021). A ground-motion database for Israel with its corresponding point-source parameters, for engineering seismology applications. Seismol. Res. Lett.92 (4), 2679–2690. 10.1785/0220200477
178
ZakI.FreundR. (1981). Asymmetry and basin migration in the Dead Sea rift. Tectonophysics80 (1–4), 27–38. 10.1016/0040-1951(81)90140-2
179
ZhuangJ.OgataY.WangT. (2017). Data completeness of the Kumamoto earthquake sequence in the JMA catalog and its influence on the estimation of the ETAS parameters. Earth Planets Space69, 36. 10.1186/s40623-017-0614-6
180
ZilbermanE.BaerG.AvniY.FeiginD. (1996). Pliocene fluvial systems and tectonics in the central Negev, southern Israel. Isr. J. Earth Sci.45 (3), 113–126.
181
ZoharM.SalamonA.RubinR. (2017). Earthquake damage history in Israel and its close surrounding - evaluation of spatial and temporal patterns. Tectonophysics696–697, 1–13. 10.1016/j.tecto.2016.12.015
182
ZoharM.SalamonA.RubinR. (2016). Reappraised list of historical earthquakes that affected Israel and its close surroundings. J. Seismol.20 (3), 971–985. 10.1007/s10950-016-9575-7
183
ZúñigaF. R.WiemerS. (1999). Seismicity patterns : Are they always related to natural causes ?Pure Appl. Geophys.155, 713–726. 10.1007/s000240050285
Summary
Keywords
seismic moment release, seismogenic zones, seismic gap, spatial variation analysis, temporal variation
Citation
Sharon M, Kurzon I, Wetzler N, Sagy A, Marco S and Ben-Avraham Z (2022) Variations of the seismic b-value along the Dead Sea transform. Front. Earth Sci. 10:1074729. doi: 10.3389/feart.2022.1074729
Received
19 October 2022
Accepted
30 November 2022
Published
21 December 2022
Volume
10 - 2022
Edited by
Guido Maria Adinolfi, University of Sannio, Italy
Reviewed by
Matteo Taroni, Istituto Nazionale di Geofisica e Vulcanologia (INGV), Italy
Alessandro Vuan, Istituto Nazionale di Oceanografia e di Geofisica Sperimentale (Italy), Italy
Updates
Copyright
© 2022 Sharon, Kurzon, Wetzler, Sagy, Marco and Ben-Avraham.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ittai Kurzon, ittaik@gsi.gov.il
This article was submitted to Structural Geology and Tectonics, a section of the journal Frontiers in Earth Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.