- 1Geological Survey of Israel, Jerusalem, Israel
- 2Department of Geophysics, Tel Aviv University, Tel Aviv, Israel
- 3Institute of Earth Sciences, Hebrew University of Jerusalem, Jerusalem, Israel
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; Beall et al., 2022), suggesting an inverse correlation of the b-value with the differential stress, considering the Anderson theory of faulting (Anderson, 1905). The inverse dependency of the b-value with the differential stress has been also shown through rock fracture experiments (Scholz, 1968; Amitrano, 2003; 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, 1991; 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; Bala et al., 2003; Bus et al., 2009) and may also be implemented for hazard evaluation purposes (Frankel, 1995; Helmstetter et al., 2006; Yadav et al., 2012; Ashish et al., 2016; 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,
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. Seismic map of the study area showing the main fault segments of the Dead Sea transform (DST) and of the northwest orientated Carmel-Tirza fault system (CTF) system (Sharon et al., 2020); recorded seismicity of 1983–2020 (expansion of Wetzler and Kurzon 2016 catalogue), highlighting the most significant events in the past century: the 1927 M six at the northern Dead Sea, and the MW 7.2 at the Gulf of Elat (Aqaba). Black arrows denote the relative plate motion along the DST.
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; Ben-Avraham and Ginzburg, 1990; 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; Bartov, 1974; 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 (Bartov et al., 1998; Avni et al., 2000, 2001; 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 (Ambraseys, 2009; Agnon, 2014; 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 (Aldersons et al., 2003; Shalev et al., 2007, 2013; Aldersons and Ben-Avraham, 2014; 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
TABLE 1. The number of
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 (
FIGURE 2. Seismic stations utilised for recording the earthquakes of the examined catalogue, and the ensuing seismic network coverage area (NCA) delimited by the green polygon.
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
Although it is more common to convert magnitudes through linear regressions, Ataeva et al. (2015) 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 Ataeva et al. (2015) were deduced in the forms of:
and
for the linear and quadratic regressions, respectively. We convert M0 to Mw following the relationship (Hanks and Kanamori, 1979; Aki and Richards, 2002):
referring to
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 Ataeva et al. (2015) 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 Ataeva et al. (2015) results are based only on ∼100 earthquakes, compared with ∼1,800 earthquakes for the current results (Table 1).
TABLE 2. Coefficient correlation parameters (a, b, c) that correspond to linear and quadratic regressions of the forms
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. Linear and Quadratic fits of the current study, within the RA. Black lines in (A,C) are linear and quadratic fits, respectively, within RA, superimposed on the associated earthquake data. Dashed yellow lines represent 1:1 ratio. Residuals from these fits are scattered for the linear (B) and quadratic (D) fits; where the black line is the linear fit of these residuals, dashed yellow line represents 1:1 ratio; Green and Red lines show the detailed RMSE in intervals of 25 events in an ascending order of magnitudes, and in intervals of a single magnitude, respectively, calculated separately for the negative and positive residuals.
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 Aki (1965) 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. Moment magnitudes arranged by the sequential number of events. Mw is set to the obtained magnitude of completeness, Mc=2.1, demonstrating that there is no clear short-term magnitude incompleteness in the data used for the following analysis, within the network coverage area (NCA).
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
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; Ben-Avraham and Ginzburg, 1990; 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. Tectonic zones and their b-values, within the NCA. DST is in yellow, CTF is in purple; Off Fault Zones (OFZ) have two sections: south [OFZ (S), in orange], and north [OFZ (N), in green] of the CTF, differing in their prevailing style of faulting. Note the partial overlap between DST and both CTF and OFZ (N). Bold and thin black lines are the main seismic sources and Quaternary faults, respectively (Sharon et al., 2020). Also shown are focal mechanism diagrams, demonstrating the dominant style of faulting in each tectonic zone, inferred from geological, seismological and geodetic observations (see text for further details).
FIGURE 6. The frequency-magnitude relation for all earthquakes recorded between 1983 and 2020, within (A) RA; (B) NCA; (C) DST; (D) CTF; (E) OFZ (S); (F) OFZ (N). Dashed parts are extrapolations. Note that the y-axis marks the annual rate of seismicity, and not the total number of events. The tectonic zones of sub-plots c–f are presented in Figure 5.
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; Achmon, 1986; 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; Bartov, 1974; Zilberman et al., 1996; Avni et al., 2000; 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 (
TABLE 4. Analysis of statistical significance for the Tectonic Zones, showing the significance level (
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 (
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. Top part: Map showing the DST polygon (yellow), defining the data for the analysis. The black lines show the main sources (see Sharon et al., 2020). Bottom part: Spatial variations of the b-value (in blue), of the 75th and 95th hypocentre depth percentiles (dashed yellow and red lines, respectively) and of the cumulative seismic moment release (in Green). Note the left axes are directed downwards. Solid and dashed lines of the b-value indicate at least 100 events, and 50–99 events per sample, respectively, and are shown with their light blue error bar band surrounding them; the blue linear dashed curve marks the spatial trend of the b-value. Within the Top part, the brown polygon marks zones with high values of seismic density and seismic moment density (see Sharon et al., 2020), overlapping the DST polygon. These zones emphasize two seismic gaps: 1) HSG—Hazeva Seismic Gap (latitudes ∼30.7°–30.8°N), and 2) BSG—Beit She’an Seismic Gap (latitudes ∼32.4°–32.6°N); the latter is prominent enough to also leave a gap within the b-value profile (Bottom Part).
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
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 Mc 2.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. (A) Temporal variations of the b-value (in blue) with their light blue error bar band around them, and of the cumulative seismic moment release (in green; note the axis on the left is directed downwards) in the DST (Figures 5, 7), based on 100-event time windows in 5-event intervals. Also marked are the b-values calculated for specific time-windows [(B–E); dashed grey lines], verifying the fluctuations seen in the b-value profile; their frequency magnitude curves are presented in panel 8 (B–E). Two reductions of the b-value can be related to swarm activities at the Sea of Galilee, during 2013 and 2018 (dashed orange rectangles; Wetzler et al., 2019).
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 (Aldersons et al., 2003; Shalev et al., 2013; Aldersons and Ben-Avraham, 2014; 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 (Amitrano, 2003) 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 (Amitrano, 2003; 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 (Avni et al., 2000; 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
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 Anderson (1905) 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; Amitrano, 2003; 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.
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
Achmon, M. (1986). The Carmel border fault between yoqneam and nesher. M.Sc. thesis. Jerusalem: Hebrew University of Jerusalem.
Agnon, A. (2014). “Pre-instrumental earthquakes along the Dead Sea rift,”in Dead Sea transform fault system: Reviews. Editors Z. Garfunkel, Z. Ben-Avraham, and E. Kagan (Dordrecht: Springer Netherlands), 207–262.
Aki, K. (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.
Aki, K., and Richards, P. G. (2002). Quantitative seismology. Sausalito, CA: University Science Books.
Akaike, H. (1974). A new look at the statistical model identification. IEEE Transactions on Automatic Control 19 (6), 716–723. doi:10.1109/TAC.1974.1100705
Aldersons, F., Ben-Avraham, Z., Hofstetter, A., Kissling, E., and Al-Yazjeen, T. (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. doi:10.1016/S0012-821X(03)00381-9
Aldersons, F., and Ben-Avraham, Z. (2014). “The seismogenic thickness in the Dead Sea area,”in Dead Sea transform fault system: Reviews. Editors Z. Garfunkel, Z. Ben-Avraham, and E. J. Kagan (Dordrecht: Springer Netherlands), 53–89. doi:10.1007/978-94-017-8872-4
Ambraseys, N. (2009). Earthquakes in the mediterranean and Middle East. Cambridge: Cambridge University Press. doi:10.1017/CBO9781139195430
Amitrano, D. (2003). Brittle-ductile transition and associated seismicity: Experimental and numerical studies and relationship with the b value. J. Geophys. Res. 108 (B1), 1–15. doi:10.1029/2001jb000680
Anderson, E. M. (1905). The dynamics of faulting. Trans. Edinb. Geol. Soc. 8 (3), 387–402. doi:10.1144/transed.8.3.387
Ashish, C., Lindholm, I. A., ParvezKühn, D., and Kuhn, D. (2016). Probabilistic earthquake hazard assessment for Peninsular India. J. Seismol. 20 (2), 629–653. doi:10.1007/s10950-015-9548-2
Ataeva, G., Shapira, A., and Hofstetter, A. (2015). Determination of source parameters for local and regional earthquakes in Israel. J. Seismol. 19 (2), 389–401. doi:10.1007/s10950-014-9472-x
Avni, Y., Bartov, Y., Garfunkel, Z., and Ginat, H. (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. doi:10.1560/W8WL-JU3Y-KM7W-8LX4
Avni, Y., Bartov, Y., Garfunkel, Z., and Ginat, H. (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. doi:10.1092/5U6A-RM5E-M8E3-QXM7
Bala, A., Radulian, M., and Popescu, E. (2003). Earthquakes distribution and their focal mechanism in correlation with the active tectonic zones of Romania. J. Geodyn. 36 (1–2), 129–145. doi:10.1016/S0264-3707(03)00044-9
Bartov, Y. (1974). A Structural and paleogeographical study of the central Sinai faults and domes. Hebrew, Englilsh abstract. PhD thesis (Jerusalem: Hebrew University of Jerusalem).
Bartov, Y., Avni, Y., Calvo, R., and Frieslander, U. (1998). The Zofar fault - a major intra-rift feature in the Arava rift valley,. Geol. Soc. Isr. Curr. Res. 11, 27–32.
Beall, A., van den Ende, M., Ampuero, J. P., Capitanio, F. A., and Åfagereng, (2022). Linking earthquake magnitude-frequency statistics and stress in visco-frictional fault zone models. Geophys. Res. Lett. 49 (20), e2022GL099247. doi:10.1029/2022gl099247
Ben-Avraham, Z., and Ginzburg, A. (1990). Displaced terranes and crustal evolution of the Levant and the eastern Mediterranean. Tectonics 9 (4), 613–622. doi:10.1029/TC009i004p00613
Ben-Menahem, A. (1991). Four thousand years of seismicity along the Dead Sea rift. J. Geophys. Res. 96 (B12), 20195–20216. doi:10.1029/91jb01936
Ben-Menahem, A. (1981). Variation of slip and creep along the levant rift over the past 4500 years. Tectonophysics 80 (1–4), 183–197. doi:10.1016/0040-1951(81)90149-9
Bentor, Y. K., and Vroman, A. (1954). A Structural contour map of Israel (1:250, 000) with remarks on its dynamic interpretation. Bull. Res. Counc. Isr. 4 (2), 125–135.
Bender, B. (1983). Maximum likelihood estimation of b values for magnitude grouped data. Bull. Seismol. Soc. Am. 73 (3), 831–851.
Bora, D. K., Borah, K., Mahanta, R., and Madhab, J. (2018). Seismic b - values and its correlation with seismic moment and Bouguer gravity anomaly over Indo-Burma ranges of northeast India : Tectonic implications. Tectonophysics 728–729, 130–141. doi:10.1016/j.tecto.2018.01.001
Bormann, P., and Di Giacomo, D. (2011). The moment magnitude Mw and the energy magnitude Me: Common roots and differences,. J. Seismol. 15 (2), 411–427. doi:10.1007/s10950-010-9219-2
Bus, Z., Grenerczy, G., Tóth, L., and Mónus, P. (2009). Active crustal deformation in two seismogenic zones of the Pannonian region - GPS versus seismological observations. Tectonophysics 474 (1–2), 343–352. doi:10.1016/j.tecto.2009.02.045
Bustos, J. P. A., Taroni, M., and Adam, L. (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. doi:10.1785/0220220004
Byerlee, J. D. (1968). Brittle-ductile transition in rocks. J. Geophys. Res. 73 (14), 4741–4750. doi:10.1029/JB073i014p04741
Calvo, R., and Bartov, Y. (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. doi:10.1092/B02L-6K04-UFQL-KUE3
Cao, A., and Gao, S. S. (2002). Temporal variation of seismic b-value beneath northeastern Japan Island Arc. Geophys. Res. Lett. 29 (9), 48-1–48-3. doi:10.1029/2001GL013775
Castellaro, S., Mulargia, F. S., and Kagan, Y. Y. (2006). Regression problems for magnitudes. Geophys. J. Int. 165 (3), 913–930. doi:10.1111/j.1365-246X.2006.02955.x
Castello, B., Olivieri, M., and Selvaggi, G. (2007). Local and duration magnitude determination for the Italian earthquake catalog, 1981-2002. Bull. Seismol. Soc. Am. 97 (1), 128–139. doi:10.1785/0120050258
Choy, G. L., and Boatwright, J. L. (1995). Global patterns of radiated seismic energy and apparent stress. J. Geophys. Res. 100 (B9), 18205–18228. doi:10.1029/95jb01969
Clauset, A., Shalizi, C. R., and Newman, M. E. J. (2009). Power-law distributions in empirical data. SIAM Rev. 51 (4), 661–703. doi:10.1137/070710111
Das, R., Wason, H. R., Gonzalez, G., Sharma, M. L., Choudhury, D., Lindholm, C., et al. (2018). Earthquake magnitude conversion problem. Bull. Seismol. Soc. Am. 108 (4), 1995–2007. doi:10.1785/0120170157
El-Isa, Z. H., and Eaton, D. W. (2014). Spatiotemporal variations in the b-value of earthquake magnitude-frequency distributions: Classification and causes. Tectonophysics 615–616, 1–11. doi:10.1016/j.tecto.2013.12.001
Felzer, K. R., Abercrombie, R. E., and Ekström, G. (2004). A common origin for aftershocks, foreshocks, and multiplets. Bull. Seismol. Soc. Am. 94 (1), 88–98. doi:10.1785/0120030069
Frankel, A. (1995). Mapping seismic hazard in the central and eastern United States. Seismol. Res. Lett. 66 (4), 8–21. doi:10.1785/gssrl.66.4.8
Frohlich, C., and Davis, S. D. (1993). Teleseismic b values; or, much ado about 1.0. J. Geophys. Res. 98 (B1), 631–644. doi:10.1029/92jb01891
Garfunkel, Z. (1981). Internal structure of the Dead Sea leaky transform (rift) in relation to plate kinematics. Tectonophysics 80 (1–4), 81–108. doi:10.1016/0040-1951(81)90143-8
Garfunkel, Z. (2014).“Lateral motion and deformation along the Dead Sea transform,” in Dead Sea transform fault system: Reviews. Editors Z. Garfunkel, Z. Ben-Avraham, and E. Kagan (Dordrecht: Springer Netherlands), 109–150.
Garfunkel, Z. (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. doi:10.1560/IJES.58.3-4.217
Garfunkel, Z., and Bartov, Y. (1977). The tectonics of the Suez rift,. Geol. Surv. Isr. Bull. 71, 1–44.
Garfunkel, Z., and Ben-Avraham, Z. (2001). Basins along the Dead Sea transform, mémoires du Muséum natl. d’histoire Nat. 186, 607–627.
Garfunkel, Z., and Ben-Avraham, Z. (1996). The structure of the Dead Sea basin,. Tectonophysics 266 (1-4), 155–176. doi:10.1016/s0040-1951(96)00188-6
Gauriau, J., and Dolan, J. 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. doi:10.1029/2021GC009938
Gerstenberger, M., Wiemer, S., and Giardini, D. (2001). A systematic test of the hypothesis that the b value varies with depth in California. Geophys. Res. Lett. 28 (1), 57–60. doi:10.1029/2000GL012026
Gibowicz, S. J. (1974). Frequency-magnitude, depth, and time relations for earthquakes in an island arc: North Island, New Zealand. Tectonophysics 23 (3), 283–297. doi:10.1016/0040-1951(74)90028-6
Ginat, H., Zilberman, E., and Amit, R. (2002). Red sedimentary units as indicators of Early Pleistocene tectonic activity in the southern Negev desert, Israel. Geomorphology 45 (1–2), 127–146. doi:10.1016/S0169-555X(01)00193-3
Ginat, H., Zilberman, E., and Avni, Y. (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. doi:10.1560/N2P9-YBJ0-Q44Y-GWYN
Gitterman, Y., Pinsky, V., Shapira, A., Ergin, M., Kalafat, D., Gurbuz, G., 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.
Glaister, P. (2005). The use of orthogonal distances in generating the total least squares estimate. Math. Comput. Educ. 39 (1), 21–30.
Goebel, T. H. W., Hosseini, S. M., Cappa, F., Hauksson, E., Ampuero, J. P., Aminzadeh, F., et al. (2016). Wastewater disposal and earthquake swarm activity at the southern end of the Central Valley, California. Geophys. Res. Lett. 43, 1092–1099. doi:10.1002/2015GL066948
Goebel, T. H. W., Kwiatek, G., Becker, T. W., Brodsky, E. E., and Dresen, G. (2017). What allows seismic events to grow big? Insights from b-value and fault roughness analysis in laboratory stick-slip experiments. Geology 45 (9), 815–818. doi:10.1130/G39147.1
Goebel, T. H. W., Schorlemmer, D., Becker, T. W., Dresen, G., and Sammis, C. G. (2013). Acoustic emissions document stress changes over many seismic cycles in stick-slip experiments. Geophys. Res. Lett. 40 (10), 2049–2054. doi:10.1002/grl.50507
Gomez, F., Karam, G., Khawlie, M., McClusky, S., Vernant, P., Reilinger, R., 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. doi:10.1111/j.1365-246X.2006.03328.x
Gomez, F., Meghraoui, M., Darkal, A. N., Hijazi, F., Mouty, M., Suleiman, Y., 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. doi:10.1046/j.1365-246X.2003.01933.x
Gulia, L., Rinaldi, A. P., Tormann, T., Vannucci, G., Enescu, B., and Wiemer, S. (2018). The effect of a mainshock on the size distribution of the aftershocks. Geophys. Res. Lett. 45, 13,277–13,287. doi:10.1029/2018GL080619
Gulia, L., and Wiemer, S. (2010). The influence of tectonic regimes on the earthquake size distribution: A case study for Italy. Geophys. Res. Lett. 37 (10), 1–6. doi:10.1029/2010GL043066
Gutenberg, B., and Richter, C. F. (1944). Frequency of earthquakes in California. Bull. Seismol. Soc. Am. 34, 185–188. doi:10.1785/bssa0340040185
Hamiel, Y., Amit, R., Begin, Z. B., Marco, S., Katz, O., Salamon, A., et al. (2009). The seismicity along the Dead Sea fault during the last 60, 000 years. Bull. Seismol. Soc. Am. 99 (3), 2020–2026. doi:10.1785/0120080218
Hamiel, Y., Masson, F., Piatibratova, O., and Mizrahi, Y. (2018). GPS measurements of crustal deformation across the southern Arava Valley section of the Dead Sea Fault and implications to regional seismic hazard assessment. Tectonophysics 724–725, 171–178. doi:10.1016/j.tecto.2018.01.016
Hamiel, Y., Piatibratova, O., and Mizrahi, Y. (2016). Creep along the northern Jordan valley section of the Dead Sea fault. Geophys. Res. Lett. 43 (6), 2494–2501. doi:10.1002/2016GL067913
Hanks, T. C., and Kanamori, H. (1979). A moment magnitude scale,. J. Geophys. Res. 84 (B5), 2348–2350. doi:10.1029/JB084iB05p02348
Helmstetter, A., Kagan, Y. Y., and Jackson, D. D. (2006). Comparison of short-term and time-independent earthquake forecast models for southern California. Bull. Seismol. Soc. Am. 96 (1), 90–106. doi:10.1785/0120050067
Henderson, J., and Main, I. (1992). A simple fracture-mechanical model for the evolution of seismicity. Geophys. Res. Lett. 19 (4), 365–368. doi:10.1029/92gl00274
Hofstetter, A., Thio, H. K., and Shamir, G. (2003). Source mechanism of the 22/11/1995 gulf of aqaba earthquake and its. J. Seismol. 7, 99–114. doi:10.1023/a:1021206930730
Hofstetter, A., Van Eck, T., and Shapira, A. (1996). Seismic activity along fault branches of the Dead Sea-Jordan transform system: The Carmel-Tirtza fault system. Tectonophysics 267 (1–4), 317–330. doi:10.1016/S0040-1951(96)00108-4
Hofstetter, R., Gitterman, Y., Pinksy, V., Kraeva, N., and Feldman, L. (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. doi:10.1560/IJES.57.2.101
Hofstetter, R., Klinger, Y., Amrat, A. Q., Rivera, L., and Dorbath, L. (2007). Stress tensor and focal mechanisms along the Dead Sea fault and related structural elements based on seismological data. Tectonophysics 429 (3–4), 165–181. doi:10.1016/j.tecto.2006.03.010
Hurvich, C. M., and Tsai, C.-L. (1989). Regression and time series model selection in small samples. Biometrika 76, 297–307. doi:10.1093/biomet/76.2.297
Hutton, K., Woessner, J., and Hauksson, E. (2010). Earthquake monitoring in southern California for seventy-seven years (1932–2008). Bull. Seismol. Soc. Am. 100 (2), 423–446. doi:10.1785/0120090130
Joffe, S., and Garfunkel, Z. (1987). Plate kinematics of the circum Red Sea-a re-evaluation,. Tectonophysics 141 (1–3), 5–22. doi:10.1016/0040-1951(87)90171-5
Kadirioğlu, F. T., and Kartal, R. 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. doi:10.3906/yer-1511-7
Kanamori, H. (1983). Magnitude scale and quantification of earthquakes. Tectonophysics 93 (3–4), 185–199. doi:10.1016/0040-1951(83)90273-1
Kanamori, H. (1977). The energy release in great earthquakes. J. Geophys. Res. 82 (20), 2981–2987. doi:10.1029/jb082i020p02981
Kane, M. T., and Mroch, A. A. (2020). Orthogonal regression, the cleary criterion, and lord's paradox: Asking the right questions. ETS Res. Rep. Ser. 2020, 1–24. doi:10.1002/ets2.12298
King, G. (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. doi:10.1007/BF02590182
Kumar, R., Yadav, R. B. S., and Castellaro, S. (2020). Regional earthquake magnitude conversion relations for the himalayan seismic belt. Seismol. Res. Lett. 91 (6), 3195–3207. doi:10.1785/0220200204
Kurzon, I., Nof, R. N., Laporte, M., Lutzky, H., Polozov, A., Zakosky, D., 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. doi:10.1785/0220200169
Lefevre, M., Klinger, Y., Al-Qaryouti, M., Le Béon, M., and Moumani, K. (2018). Slip deficit and temporal clustering along the Dead Sea fault from paleoseismological investigations. Sci. Rep. 8 (1), 4511–4519. doi:10.1038/s41598-018-22627-9
Lei, X., Kusunose, K., Rao, M. V. M. S., Nishizawa, O., and Satoh, T. (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. doi:10.1029/1999jb900385
Llana-Fúnez, S., and Lõpez-Fernández, C. (2015). The seismogenic zone of the continental crust in Northwest Iberia and its relation to crustal structure. Tectonics 34 (8), 1751–1767. doi:10.1002/2015TC003877
Lu, Y., Wetzler, N., Waldmann, N., Agnon, A., Biasi, G. P., and Marco, S. (2020). A 220, 000-year-long continuous large earthquake record on a slow-slipping plate boundary. Sci. Adv. 6 (48), eaba4170. doi:10.1126/sciadv.aba4170
Magrin, A., Peresan, A., Kronrod, T., Vaccari, F., and Panza, G. F. (2017). Neo-deterministic seismic hazard assessment and earthquake occurrence rate,. Eng. Geol. 229, 95–109. doi:10.1016/j.enggeo.2017.09.004
Main, I. G., Meredith, P. G., and Jones, C. (1989). A reinterpretation of the precursory seismic b-value anomaly from fracture mechanics. Geophys. J. Int. 96 (1), 131–138. doi:10.1111/j.1365-246X.1989.tb05255.x
Maiti, S. K., and Kamai, R. (2020). Interaction between fault and off-fault seismic sources in hazard analysis – a case study from Israel. Eng. Geol. 274, 105723. doi:10.1016/j.enggeo.2020.105723
Mandal, P., Srinagesh, D., Suresh, G., Naresh, B., Naidu, M., Singh, D. 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. Hazards 108, 1183–1196. doi:10.1007/s11069-021-04726-5
Marco, S. (2007). Temporal variation in the geometry of a strike-slip fault zone: Examples from the Dead Sea Transform. Tectonophysics 445 (3–4), 186–199. doi:10.1016/j.tecto.2007.08.014
Marco, S., and Klinger, Y. (2014). “Review of on-fault palaeoseismic studies along the Dead Sea fault,”in Dead Sea transform fault system: Reviews. Editors Z. Garfunkel, Z. Ben-Avraham, and E. Kagan (Dordrecht: Springer Netherlands), 183–205. doi:10.1007/978-94-017-8872-4
Marzocchi, W., and Taroni, M. (2014). Some thoughts on declustering in probabilistic seismic-hazard analysis. Bull. Seismol. Soc. Am. 104 (4), 1838–1845. doi:10.1785/0120130300
Marzocchi, W., and Sandri, L. (2003). A review and new insights on the estimation of the b value and its uncertainty. Ann. Geophys. 46 (6), 1271–1282.
Marzocchi, W., Spassiani, I., Stallone, A., and Taroni, M. (2020). How to be fooled searching for significant variations of the b-value. Geophys. J. Int. 220 (3), 1845–1856. doi:10.1093/gji/ggz541
Michailos, K., Smith, E. G. C., Chamberlain, C. J., Savage, M. K., and Townend, J. (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. doi:10.1029/2018GC007743
Mizrahi, L., Nandan, S., and Wiemer, S. (2021). The effect of declustering on the size distribution of mainshocks. Seismol. Res. Lett. 92 (4), 2333–2342. doi:10.1785/0220200231
Mori, J., and Abercrombie, R. E. (1997). Depth dependence of earthquake frequency-magnitude distributions in California: Implications for rupture initiation. J. Geophys. Res. 102 (B7), 15081–15090. doi:10.1029/97jb01356
Mousavi, S. M., Ogwari, P. O., Horton, S. P., and Langston, C. A. (2017). Spatio-temporal evolution of frequency-magnitude distribution and seismogenic index during initiation of induced seismicity at Guy-Greenbrier, Arkansas. Phys. Earth Planet. Interiors 267, 53–66. doi:10.1016/j.pepi.2017.04.005
Mukhopadhyay, B., Mogren, S., Mukhopadhyay, M., and Dasgupta, S. (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 Risk 4 (1), 30–48. doi:10.1080/19475705.2012.663794
Neev, D., Almagor, G., Arad, A., Ginzburg, A., and Hall, J. K. (1976). The geology of the southeastern Mediterranean Sea,. Geol. Surv. Isr. Bull. 68, 1–51.
Norris, R. J., and Cooper, A. F. (2001). Late Quaternary slip rates and slip partitioning on the Alpine Fault, New Zealand. J. Struct. Geol. 23 (2), 507–520. doi:10.1016/s0191-8141(00)00122-x
Nuannin, P., Kulhanek, O., and Persson, L. (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. doi:10.1029/2005GL022679
Nuriel, P., Weinberger, R., Kylander-Clark, A. R. C., Hacker, B. R., and Craddock, J. P. (2017). The onset of the Dead Sea transform based on calcite age-strain analyses,. Geology 45 (7), 587–590. doi:10.1130/G38903.1
Öztürk, S. (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. doi:10.2478/s11600-010-0050-5
Page, M., and Felzer, K. (2015). Southern San Andreas fault seismicity is consistent with the Gutenberg–Richter magnitude–frequency distribution. Bull. Seismol. Soc. Am. 105 (4), 2070–2080. doi:10.1785/0120140340
Pavlis, G. L. (1986). Appraising earthquake hypocenter location errors: A complete, practical approach for single-event locations. Bull. Seismol. Soc. Am. 76 (6), 1699–1717.
Pavlis, G. L., Vernon, F. L., Harvey, D., and Quinlan, D. (2004). The generalized earthquake-location (GENLOC) package: An earthquake-location library. Comput. Geosci. 30 (9/10), 1079–1091. doi:10.1016/j.cageo.2004.06.010
Petersen, M. D., Dawson, T. E., Chen, R., Cao, T., Wills, C. J., Schwartz, D. P., et al. (2011). Fault displacement hazard for strike-slip faults. Bull. Seismol. Soc. Am. 101 (2), 805–825. doi:10.1785/0120100035
Petruccelli, A., Schorlemmer, D., Tormann, T., Rinaldi, A. P., Wiemer, S., Gasperini, P., et al. (2019b). The influence of faulting style on the size-distribution of global earthquakes,. Earth Planet. Sci. Lett. 527, 115791. doi:10.1016/j.epsl.2019.115791
Petruccelli, A., Gasperini, P., Tormann, T., Schorlemmer, D., Rinaldi, A. P., Vannucci, G., et al. (2019a). Simultaneous dependence of the earthquake-size distribution on faulting style and depth,. Geophys. Res. Lett. 46 (20), 11044–11053. doi:10.1029/2019GL083997
Petruccelli, A., Vannucci, G., Lolli, B., and Gasperini, P. (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. doi:10.1785/0120170328
Quennell, A. M. (1959). “Tectonics of the Dead Sea rift,” in Proc. 20th int. Geol. Congr. Mex., 385–403.
Radulian, M., Bălă, A., Popescu, E., and Toma-Dănilă, D. (2018). Earthquake mechanism and characterization of seismogenic zones in south-eastern part of Romania. Ann. Geophys. 61 (1). doi:10.4401/ag-7443
Radulian, M., Mândrescu, N., Panza, G. F., Popescu, E., and Utale, A. (2000). Characterization of seismogenic zones of Romania. Pure Appl. Geophys. 157 (1–2), 57–77. doi:10.1007/978-3-0348-8415-0_4
Raub, C., Martínez-Garzón, P., Kwiatek, G., Bohnhoff, M., and Dresen, G. (2017). Variations of seismic b-value at different stages of the seismic cycle along the North Anatolian Fault Zone in northwestern Turkey. Tectonophysics 712, 232–248. doi:10.1016/j.tecto.2017.05.028
Reilinger, R., McClusky, S., Vernant, P., Lawrence, S., Ergintav, S., Cakmak, R., 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. doi:10.1029/2005JB004051
Reilinger, R.E., McClusky, S. C., Oral, M. B., King, R. W., Toksoz, M. N., Barka, A. 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 Earth 102 (B5), 9983–9999.
Rigo, A., Souriau, A., and Sylvander, M. (2018). Spatial variations of b-value and crustal stress in the Pyrenees. J. Seismol. 22, 337–352. doi:10.1007/s10950-017-9709-6
Rivière, J., Lv, Z., Johnson, P. A., and Marone, C. (2018). Evolution of b-value during the seismic cycle: Insights from laboratory experiments on simulated faults. Earth Planet. Sci. Lett. 482, 407–413. doi:10.1016/j.epsl.2017.11.036
Ron, H., and Eyal, Y. (1985). Intraplate deformation by block rotation and mesostructures along the Dead Sea transform, northern Israel. Tectonics 4 (1), 85–105. doi:10.1029/tc004i001p00085
Ron, H., Freund, R., Garfunkel, Z., and Nur, A. (1984). Block rotation by strike-slip faulting: Structural and paleomagnetic evidence. J. Geophys. Res. 89B7, 6256–6270. doi:10.1029/JB089iB07p06256
Ross, Z. E., Ben-Zion, Y., White, M. C., and Vernon, F. 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. doi:10.1093/gji/ggw327
Rotstein, Y., Bruner, I., and Kafri, U. (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.
Sadeh, M., Hamiel, Y., Ziv, A., Bock, Y., Fang, P., and Wdowinski, S. (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. doi:10.1029/2012JB009241
Salamon, A., Hofstetter, A., Garfunkel, Z., and Ron, H. (1996). Seismicity of the eastern mediterranean region: Perspective from the Sinai subplate. Tectonophysics 263 (1–4), 293–305. doi:10.1016/S0040-1951(96)00030-3
Scholz, C. H. (2015). On the stress dependence of the earthquake b value. Geophys. Res. Lett. 42 (5), 1399–1402. doi:10.1002/2014GL062863
Scholz, C. H. (1968). The frequency - magnitude relation of microfracturing in rock and its relation to earthquakes. Bull. Seismol. Soc. Am. 58 (1), 399–415. doi:10.1785/bssa0580010399
Schorlemmer, D., Wiemer, S., Kleiber, P., Maunder, M. N., and Harley, S. J. (2005). Microseismicity data forecast rupture area. Nature 434 (7037), E1–E2. doi:10.1038/nature03581
Schorlemmer, D., Wiemer, S., and Wyss, M. (2004). Earthquake statistics at parkfield: 1. Stationarity of b values. J. Geophys. Res. 109 (12), B12307–B12317. doi:10.1029/2004JB003234
Schorlemmer, D., Wiemer, S., and Wyss, M. (2005). Variations in earthquake-size distribution across different stress regimes. Nature 437 (7058), 539–542. doi:10.1038/nature04094
Schwartz, D. P., and Coppersmith, K. J. (1984). Fault behavior and characteristic earthquakes: Examples from the wasatch and san Andreas Fault Zones. J. Geophys. Res. 89 (B7), 5681–5698. doi:10.1029/JB089iB07p05681
Schwartz, S. Y., Orange, D. L., and Anderson, R. 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. doi:10.1029/gl017i008p01207
Scordilis, E. M. (2006). Empirical global relations converting MS and mb to moment magnitude. J. Seismol. 10 (2), 225–236. doi:10.1007/s10950-006-9012-4
Shalev, E., Lyakhovsky, V., Weinstein, Y., and Ben-Avraham, Z. (2013). The thermal structure of Israel and the Dead Sea fault,. Tectonophysics 602, 69–77. doi:10.1016/j.tecto.2012.09.011
Shalev, E., Lyakhovsky, V., and Yechieli, Y. (2007). Is advective heat transport significant at the Dead Sea basin? Geofluids 7 (3), 292–300. doi:10.1111/j.1468-8123.2007.00190.x
Shamir, G., Bartov, Y., Sneh, A., Fleisher, L., Arad, V., and Rosensaft, M. (2001). Preliminary seismic zonation in Israel.GII report No. 550/95/01(1) and report No. GSI12/2000, 01, 28. no. 550.
Shapira, A. (1992). Detectability of regional seismic network: Analysis of the Israel seismic networks. Isr. J. Earth-Sciences 41 (1), 21–25.
Shapira, A., Avni, R., and Nor, A. (1993). A new estimate for the epicenter of the Jericho earthquake of 11 July 1927. Isr. J. Earth-Sciences 42 (2), 93–96.
Shapira, A., and Hofstetter, A. (2002). Seismicity parameters of seismogenic zones. RELEMR-MERC Rep. 2001, no. 1980. Available at: http://www.relemr–merc.org/Other_information/seis.
Sharon, M. (2020). “Mapping and characterising active tectonic sources in Israel and adjacent areas,” in Geological survey of Israel report No. GSI/19/20 M.Sc. thesis. Jerusalem (Israel): Tel Aviv University.
Sharon, M., Sagy, A., Kurzon, I., Marco, S., and Rosensaft, M. (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. doi:10.5194/nhess-20-125-2020
Shi, Y., and Bolt, B. A. (1982). The standard error of the magnitude-frequency b value,. Bull. Seismol. Soc. Am. 72 (5), 1677–1687. doi:10.1785/bssa0720051677
Sibson, R. H. (1974). Frictional constraints on thrust, wrench and normal faults. Nature 249 (5457), 542–544. doi:10.1038/249542a0
Sibson, R. H. (1984). Roughness at the base of the seismogenic zone: Contributing factors. J. Geophys. Res. 89 (B7), 5791–5799. doi:10.1029/JB089iB07p05791
Sokolov, V., Zahran, H. M., Youssef, S. E. H., El-Hadidy, M., and Alraddadi, W. 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. doi:10.1007/s10518-016-0075-5
Spada, M., Tormann, T., Wiemer, S., and Enescu, B. (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. doi:10.1029/2012GL054198
Steckler, M. S., Berthelot, F., Lyberis, N., and Le Pichon, X. (1988). Subsidence in the gulf of suez: Implications for rifting and plate kinematics. Tectonophysics 153 (1–4), 249–270. doi:10.1016/0040-1951(88)90019-4
Suyehiro, B. 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. doi:10.1785/bssa0560010185
Tan, J. Y., Waldhauser, F., Tolstoy, M., and Wilcock, W. 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. doi:10.1016/j.epsl.2019.01.047
Tormann, T., Wiemer, S., and Hauksson, E. (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. doi:10.1785/0120090124
Tormann, T., Wiemer, S., Metzger, S., Michael, A., and Hardebeck, J. L. (2013). Size distribution of parkfield’s microearthquakes reflects changes in surface creep rate. Geophys. J. Int. 193 (3), 1474–1478. doi:10.1093/gji/ggt093
Tormann, T., Wiemer, S., and Mignan, A. (2014). Systematic survey of high-resolution b value imaging along Californian faults : Inference on asperities. J. Geophys. Res. Solid Earth 119, 2029–2054. doi:10.1002/2013JB010867
Utsu, T. (2002). 44 Relationships between magnitude scales. Int. Geophys. 81, 733–746. doi:10.1016/S0074-6142(02)80247-9
Utsu, T. (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.
Utsu, T. (1966). A statistical significance test of the difference in b-value between two earthquake groups. J. Phys. Earth 14 (2), 37–40.
Weinberger, R., Gross, M. R., and Sneh, A. (2009). Evolving deformation along a transform plate boundary: Example from the Dead Sea Fault in northern Israel. Tectonics 28 (5), 1–19. doi:10.1029/2008TC002316
Weinberger, R., Nuriel, P., Kylander-clark, A. R. C., and Craddock, J. 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. doi:10.1016/j.earscirev.2020.103377
Wesnousky, S. G., Scholz, C. H., Shimazaki, K., and Matsuda, T. (1983). Earthquake frequency distribution and the mechanics of faulting. J. Geophys. Res. B11, 9331–9340. doi:10.1029/jb088ib11p09331
Wetzler, N., and Kurzon, I. (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. doi:10.1785/0220150157
Wetzler, N., Shalev, E., Goebel, T., Amelung, F., Kurzon, I., Lyakhovsky, V., et al. (2019). Earthquake swarms triggered by groundwater extraction near the Dead Sea Fault. Geophys. Res. Lett. 46 (14), 8056–8063. doi:10.1029/2019GL083491
Wiemer, S., and Katsumata, K. (1999). Spatial variability of seismicity parameters in aftershock zones. J. Geophys. Res. 104 (B6), 13135–13151. doi:10.1029/1999jb900032
Wiemer, S., McNutt, S. R., and Wyss, M. (1998). Temporal and three-dimensional spatial analyses of the frequency-magnitude distribution near Long Valley Caldera, California. Geophys. J. Int. 134 (2), 409–421. doi:10.1046/j.1365-246X.1998.00561.x
Wiemer, S., and Wyss, M. (2002). Mapping spatial variability of the frequency-magnitude distribution of earthquakes. Adv. Geophys. 45, 259. doi:10.1016/S0065-2687(02)80007-3
Wiemer, S., and Wyss, M. (1997). Mapping the frequency-magnitude distribution in asperities: An improved technique to calculate recurrence times? J. Geophys. Res. 102 (B7), 15115–15128. doi:10.1029/97jb00726
Williams, T. B., Kelsey, H. M., and Freymueller, J. 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. Tectonophysics 413 ((3–4)), 171–184. doi:10.1016/j.tecto.2005.10.047
Wyss, M. (2001). Locked and creeping patches along the Hayward fault, California. Geophys. Res. Lett. 28 (18), 3537–3540. doi:10.1029/2001gl013499
Wyss, M. (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. doi:10.1785/0220200008
Wyss, M. (1973). Towards a physical understanding of the earthquake frequency distribution, geophys. Geophys. J. Int. 31 (4), 341–359. doi:10.1111/j.1365-246X.1973.tb06506.x
Wyss, M., Sammis, C. G., Nadeau, R. M., and Wiemer, S. (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. doi:10.1785/0120030054
Xie, W., Hattori, K., and Han, P. (2019). Temporal variation and statistical assessment of the b value off the pacific coast of tokachi, hokkaido, Japan. Entropy 21 (3), 249. doi:10.3390/e21030249
Yadav, R. B. S., Bayrak, Y., Tripathi, J. N., Chopra, S., Singh, A. P., and Bayrak, E. (2012). A probabilistic assessment of earthquake hazard parameters in NW himalaya and the adjoining regions. Pure Appl. Geophys. 169 (9), 1619–1639. doi:10.1007/s00024-011-0434-8
Yadav, R. B. S., Bormann, P., Rastogi, B. K., Das, M. C., and Chopra, S. (2009). A Homogeneous and complete earthquake catalog for northeast India and the adjoining region. Seismol. Res. Lett. 80 (4), 609–627. doi:10.1785/gssrl.80.4.609
Yagoda-Biran, G., Maiti, S. K., Wetzler, N., Nof, R. N., Pashcur, Y., and Kamai, R. (2021). A ground-motion database for Israel with its corresponding point-source parameters, for engineering seismology applications. Seismol. Res. Lett. 92 (4), 2679–2690. doi:10.1785/0220200477
Zak, I., and Freund, R. (1981). Asymmetry and basin migration in the Dead Sea rift. Tectonophysics 80 (1–4), 27–38. doi:10.1016/0040-1951(81)90140-2
Zhuang, J., Ogata, Y., and Wang, T. (2017). Data completeness of the Kumamoto earthquake sequence in the JMA catalog and its influence on the estimation of the ETAS parameters. Earth Planets Space 69, 36. doi:10.1186/s40623-017-0614-6
Zilberman, E., Baer, G., Avni, Y., and Feigin, D. (1996). Pliocene fluvial systems and tectonics in the central Negev, southern Israel. Isr. J. Earth Sci. 45 (3), 113–126.
Zohar, M., Salamon, A., and Rubin, R. (2017). Earthquake damage history in Israel and its close surrounding - evaluation of spatial and temporal patterns. Tectonophysics 696–697, 1–13. doi:10.1016/j.tecto.2016.12.015
Zohar, M., Salamon, A., and Rubin, R. (2016). Reappraised list of historical earthquakes that affected Israel and its close surroundings. J. Seismol. 20 (3), 971–985. doi:10.1007/s10950-016-9575-7
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.
Edited by:
Guido Maria Adinolfi, University of Sannio, ItalyReviewed by:
Matteo Taroni, Istituto Nazionale di Geofisica e Vulcanologia (INGV), ItalyAlessandro Vuan, Istituto Nazionale di Oceanografia e di Geofisica Sperimentale (Italy), Italy
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, aXR0YWlrQGdzaS5nb3YuaWw=