Skip to main content

ORIGINAL RESEARCH article

Front. Earth Sci., 21 December 2022
Sec. Structural Geology and Tectonics
This article is part of the Research Topic Tectonic Faults: Types, Characteristics and Behavior View all 7 articles

Variations of the seismic b-value along the Dead Sea transform

  • 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):

logN=abM(1)

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, Mc, 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
www.frontiersin.org

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 27.5°35.5°N and longitudes 32°38°E, 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 0.1M7.2, from which ∼2,900 are with unknown magnitudes (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. The number of Md and Mw 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 (Mc. 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
www.frontiersin.org

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 Md and Mw; 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; Ataeva et al., 2015; 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, 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:

logM0=aMd+b(2)

and

logM0=aMd2+bMd+c(3)

for the linear and quadratic regressions, respectively. We convert M0 to Mw following the relationship (Hanks and Kanamori, 1979; Aki and Richards, 2002):

Mw=23logM09.1(4)

referring to M0 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 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
www.frontiersin.org

TABLE 2. Coefficient correlation parameters (a, b, c) that correspond to linear and quadratic regressions of the forms Mw=a*Md+b and Mw=a*Md2+b*Md+c, respectively; RMSE and AICc of these regressions are in respect to earthquakes located within RA; a) The regression coefficients of Ataeva et al. (2015) are from S-wave analysis and the magnitude range of 2.7Md5.6, and were achieved through Eqs 13 (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
www.frontiersin.org

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:

Mw=0.03Md2+0.65Md+0.69(5)

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
www.frontiersin.org

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):

b=1ln10M¯McM/2(6)

where M¯ is the sampling average of the magnitudes, Mc is the completeness magnitude, and ΔM 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):

σb=2.30b2i=1nMiM¯2/nn10.5(7)

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
www.frontiersin.org

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
www.frontiersin.org

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.

TABLE 3
www.frontiersin.org

TABLE 3. Seismic characteristics of tectonic zones (Figures 4, 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 (s.l.) is between 0 and 1; s.l.<0.01 is highly significant and s.l.=1 indicates identical datasets (e.g., Utsu, 1999).

TABLE 4
www.frontiersin.org

TABLE 4. Analysis of statistical significance for the Tectonic Zones, showing the significance level (s.l.) 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 s.l., in which the green shades show s.l.0.05 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 s.l.0.05. 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 (ΔAIC) is calculated for each dataset, and the significance level is given by the probability P=eΔAIC22 (Utsu, 1999); when p 0.05 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, σ1/σ2 (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 s.l.<=0.05 (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
www.frontiersin.org

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°. 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 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
www.frontiersin.org

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 Md to Mw. 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 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.

CrossRef Full Text | Google Scholar

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.

Google Scholar

Aki, K., and Richards, P. G. (2002). Quantitative seismology. Sausalito, CA: University Science Books.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Ambraseys, N. (2009). Earthquakes in the mediterranean and Middle East. Cambridge: Cambridge University Press. doi:10.1017/CBO9781139195430

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Anderson, E. M. (1905). The dynamics of faulting. Trans. Edinb. Geol. Soc. 8 (3), 387–402. doi:10.1144/transed.8.3.387

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

Bender, B. (1983). Maximum likelihood estimation of b values for magnitude grouped data. Bull. Seismol. Soc. Am. 73 (3), 831–851.

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Byerlee, J. D. (1968). Brittle-ductile transition in rocks. J. Geophys. Res. 73 (14), 4741–4750. doi:10.1029/JB073i014p04741

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Fisher, R. A. (1950). Contribution to Mathematical Statistics. New York, NY: Wiley.

Google Scholar

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

CrossRef Full Text | Google Scholar

Freund, R. (1970). The geometry of faulting in the Galilee. Isr. J. Earth Sci. 19 (1943), 117–140.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Garfunkel, Z., and Bartov, Y. (1977). The tectonics of the Suez rift,. Geol. Surv. Isr. Bull. 71, 1–44.

Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

Glaister, P. (2005). The use of orthogonal distances in generating the total least squares estimate. Math. Comput. Educ. 39 (1), 21–30.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Gosset, W. S. (1908). The probable error of a mean,. Biometrika 6 (1), 1–25. doi:10.2307/2331554

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Gutenberg, B., and Richter, C. F. (1944). Frequency of earthquakes in California. Bull. Seismol. Soc. Am. 34, 185–188. doi:10.1785/bssa0340040185

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Hanks, T. C., and Kanamori, H. (1979). A moment magnitude scale,. J. Geophys. Res. 84 (B5), 2348–2350. doi:10.1029/JB084iB05p02348

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Kanamori, H. (1983). Magnitude scale and quantification of earthquakes. Tectonophysics 93 (3–4), 185–199. doi:10.1016/0040-1951(83)90273-1

CrossRef Full Text | Google Scholar

Kanamori, H. (1977). The energy release in great earthquakes. J. Geophys. Res. 82 (20), 2981–2987. doi:10.1029/jb082i020p02981

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Ö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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Quennell, A. M. (1959). “Tectonics of the Dead Sea rift,” in Proc. 20th int. Geol. Congr. Mex., 385–403.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

CrossRef Full Text | Google Scholar

Rice, John A. (2006). Mathematical statistics and data analysis. 3rd ed. Duxbury Advanced.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Scholz, C. H. (2015). On the stress dependence of the earthquake b value. Geophys. Res. Lett. 42 (5), 1399–1402. doi:10.1002/2014GL062863

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

Shapira, A. (1992). Detectability of regional seismic network: Analysis of the Israel seismic networks. Isr. J. Earth-Sciences 41 (1), 21–25.

Google Scholar

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.

Google Scholar

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.

Google Scholar

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.

Google Scholar

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

CrossRef Full Text | Google Scholar

Shearer, P. M. (2009). Introduction to seismology. New York: Cambridge University Press.

Google Scholar

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

CrossRef Full Text | Google Scholar

Sibson, R. H. (1974). Frictional constraints on thrust, wrench and normal faults. Nature 249 (5457), 542–544. doi:10.1038/249542a0

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Utsu, T. (2002). 44 Relationships between magnitude scales. Int. Geophys. 81, 733–746. doi:10.1016/S0074-6142(02)80247-9

CrossRef Full Text | Google Scholar

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.

CrossRef Full Text | Google Scholar

Utsu, T. (1966). A statistical significance test of the difference in b-value between two earthquake groups. J. Phys. Earth 14 (2), 37–40.

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Wyss, M. (2001). Locked and creeping patches along the Hayward fault, California. Geophys. Res. Lett. 28 (18), 3537–3540. doi:10.1029/2001gl013499

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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.

Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Zúñiga, F. R., and Wiemer, S. (1999). Seismicity patterns : Are they always related to natural causes ? Pure Appl. Geophys. 155, 713–726. doi:10.1007/s000240050285

CrossRef Full Text | Google Scholar

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, 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

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=

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.