Discrimination Between Small Earthquakes and Quarry Blasts in Egypt Using Spectral Source Characteristics

Source parameters calculated from displacement spectra of both P and S waves are used to discriminate between earthquakes and quarry blasts in three regions of Egypt during the 2009–2015 period. We use vertical component seismograms from 440 earthquakes and 450 quarry explosions with MD 1.5 to 3.3 to calculate source parameters, including scalar moments and corner frequencies. The Mo(P,S) vs. fc(P,S) and P- to S-wave corner frequency fc(P)/fc(S) ratios are used to distinguish quarry blasts from earthquakes. A comparison of Mo(P,S) vs. fc(P,S) for both earthquakes and explosions in Egypt demonstrates that explosions had significantly lower corner frequencies than earthquakes, particularly for S-wave displacement spectra. In contrast to the Northern and Central regions, the Southern Egyptian region provides a perfect separation of corner frequencies of earthquakes and explosions for both P- and S-waves. The empirically derived average ratio of fc(P)/fc(S) for earthquakes is 1.28, 1.26 and 1.26 in the Northern, Central and Southern Egyptian regions, respectively. For explosions, average fc(P)/fc(S) ratios are 1.89, 1.86 and 2.0 in the three Egyptian regions, respectively. According to these findings, the average ratio of fc(P) to fc(S) for explosions is higher than those for earthquakes, implying that the differences in ratios enhance the ability of the fc(P) vs. fc(S) approach to discriminate between earthquakes and explosions. Based on the average fc(P)/fc(S) ratios vs. Mw in the whole of Egypt, the observed fc(P)/fc(S) discrimination threshold value for separating quarry explosions from earthquakes is 1.51–1.52.


Introduction
The Egyptian National Seismological Network (ENSN), which is spread across the country, plays a key role in monitoring seismic events, both natural earthquakes and artificial explosions. Several quarries in Egypt use blasting to handle the mobility of rock masses. Trustable discrimination of human-made explosions is presently becoming a significant issue in seismic observations by both local and regional seismic networks. Distinguishing between smallmagnitude events remains a challenging task.
Among the techniques used for discrimination are ratios of amplitudes of various seismic phases (e.g., Anderson et al., 2009;McLaughlin et al., 2004) and velocity spectra assessment (e.g., Gitterman et al., 1998;Walter et al., 1995). Other investigation methods include spectral recognition of certain wave groups, spectral ratios of the seismic phases and spectrograms (e.g., Allmann et al., 2008;Dahy & Hassib, 2010;Hammer et al., 2013;Koch & Fäh, 2002). A significant number of studies have used short-period records for discrimination in specific regions (e.g., Horasan et al., 2009;Yilmaz et al., 2013). Yavuz et al. (2019) performed the discrimination between earthquakes and explosions that occurred in the Armutlu Peninsula area, northwestern Turkey, using the power spectrum density (PSD) discrimination method. They noted that the power spectra of quarry blasts decreased more sharply with increasing frequencies than those of earthquakes. Sertçeliket al. (2020) conducted a study in the Eastern Marmara region, Turkey, using the amplitude peak ratio, complexity, spectral ratio and continuous wavelet transform (CWT) discrimination methods. They found that 182 of the recorded seismic events were caused by artificial explosions, while 395 events were classified as earthquakes.
Displacement spectra of the P-and S-waves recorded by short-period stations for both local earthquakes and quarry blasts of comparable magnitudes in Israel and Lebanon have been inspected and evaluated by Ataeva et al. (2017). According to their findings, the calculated ratios of P-to S-wave corner frequency f c (P)/f c (S) computed for each event have demonstrated significant segregation between earthquake and quarry blast groups. Explosions also have lower corner frequencies (f c ) than earthquakes, particularly for S-wave displacement spectra, as shown by the M o vs. f c relationship for earthquakes and explosions.
This study applies the spectral analysis technique for discrimination between earthquakes and quarry blasts (explosions). Following the approach of Ataeva et al. (2017), the average values of the corner frequency ratio f c (P)/f c (S) and M o vs. f c of multiple displacement spectra stations for both earthquakes and quarry explosions in three Egyptian regions (Northern, Central and Southern) are investigated to establish statistical discrimination criteria for segregation. Each region is checked individually to identify variations that can distinguish any of them. We also investigate the relationship between the f c (P)/f c (S) ratio and the moment magnitude M w of both P-and S-waves for the whole of Egypt.

Data Source
The waveform data used in this study consist of 3416 vertical component seismograms from 440 earthquakes and 450 quarry explosions with M D 1.5-3.3, which occurred at epicentral distances up to 200 km during the 2009-2015 periods. The broadband and short-period stations data of ENSN are extracted from the database. Earthquakes of depths \ 25 km were used. Most quarry blasts are identified from monthly reports of quarries received by the ENSN main center independently of seismic means. In this study, we rely on solutions under the condition that the seismic event was recorded by at least three seismic stations located at different epicentral distances with a reasonable signal-to-noise ratio (SNR C 10). We used the software code offered by Kumar et al. (2014) to determine the SNR of P-and S-waves. SNR is calculated by dividing the Fourier amplitude spectrum estimated from a P or S window of the signal by the Fourier amplitude spectrum calculated from a noisy window of the same duration, which is frequently selected before the event. The data processing software EQK_SRC_PARA (Kumar et al., 2012) was used in the analysis. The hypocentral parameters and duration magnitudes of the studied events were determined using ATLAS program, provided by Nanometrics Inc., Canada, for seismic data analysis using the 1D velocity-depth model of EI-Hadidy (1995).
ENSN stations involve different seismometers: short period SS-1 Ranger with 1-Hz natural frequency, L4C type, broadband 120 s Streckeisen STS-2 seismometers, Trillium 40 and Trillium 120 s (Fig. 1). The sampling rate of the recorded data is 100 samples per second. According to the distribution of both earthquakes and quarry explosions, Egypt is divided into three main regions (Fig. 2): Northern (including subzones 1, 2 and 3); Central (including subzones 4 and 5) and Southern Egypt (including subzone 6). For reliable and robust discrimination analysis, comparable magnitudes and close spatial distribution for earthquakes and explosions are maintained (Figs. 2, 3). Figures 4, 5 demonstrate a time-domain comparison between earthquake and explosion, using waveforms for two events having similar magnitudes and were recorded by the same stations and located almost at the same epicentral distances. It is evident that the displayed waveforms reveal the difference between an earthquake and an explosion, clarifying the quality of S-waves produced by explosions. In some cases, S-waves for quarry blasts are poorly created and S-wave windows could be possibly contaminated by long period surface waves (Rg-type energy) (Fig. 4)

Data Processing
Spectral analysis is one of the most valuable screening methods to differentiate between earthquakes and explosions. The source parameters for all events were investigated using the P-and S-wave displacement amplitude spectra obtained from the vertical component seismograms according to the omega-square source model of Brune (1970Brune ( , 1971, considering attenuation effects to geometrical spreading and anelastic attenuation. The main reason for using the vertical component is that most ENSN stations during the study period have a single vertical component. Moreover, the vertical component averts any site effects or spectral complexity. EQK_SRC_-PARA software of Kumar et al. (2012) was used to calculate the displacement amplitude spectra. The first step in the signal pre-processing is to remove the mean and linear trends. The second step is to perform deconvolution to remove the instrumental effect from each station seismogram, using a designed filter with the coefficient of poles, zeros and combined sensitivity of both sensor and digitizer. Once the signal was corrected, the seismograms were transformed from counts to m/s. Consequently, we manually picked the appropriate range of P-and S-wave time windows in seconds for each seismogram, always starting at the first arrival times of P-and S-waves, which were picked manually on the vertical component seismogram. Several checks with different durations of P-and S-waves were performed to select the optimum signal window length that would avoid interference from other phases and ensure the resolution and stabilization of amplitude spectra. The selected P-wave window length covers the first P-phase group; it ranges from 1.5 to 4 s. The S-wave's selected time window usually ranges from 10 to 20 s, depending on distances. This window includes the most significant part of the signal energy, as generally accepted for estimating the seismic moment (Baumbach & Borman, 2012;Brune, 1970Brune, , 1971Hanks & Wyss, 1972). At short distances, most quarry blast S-wave time windows also include the surface Rg phase (actually very weak S-wave ? Rg-wave windows) (Fig. 3). The amplitude Rg wave is strongly depth-dependent; it is often observed in near-surface events such as quarry blasts with large amplitudes (e.g., McLaughlin et al., 2004) and not usually observed in earthquakes with depths greater than a  2017), the long S-wave window offers the optimum results for discrimination between earthquakes and explosions. After selecting the windows, integration must be done to obtain the displacement waveforms using the cumulative trapezoidal numerical integration technique. To obtain the displacement spectra of the specific time window of P-and S-waves for each vertical velocity record, the fast Fourier transform (FFT) is carried out. The obtained displacement spectrum was corrected for both geometric spreading and anelastic attenuation. If the attenuation correction was not neglected, two effects were produced; the spectrum form changes, which affects the corner frequency assessment and influences the low-frequency level, which consequently affects the seismic moment calculation. Therefore, the correction for attenuation is particularly important to obtain a reliable corner frequency. It also contributes to correcting the flat part of the spectrum in the case of small earthquakes with high-frequency content (Süle & Wébe, 2013). In this study, based on the equation G(R) = 1/R, geometrical spreading is assessed for body waves Observed P-wave (a) and S-wave (b) displacement spectra of an earthquake (28 October 2012, M d = 2.5) before (green) and after (red) correction for attenuation effects together with the best-fitting theoretical Brune's spectra (blue). This is calculated from seismograms of the HAG station located 50 km from the earthquake epicenter. c Observed S-wave displacement spectra of an earthquake (31 January 2012, M d = 1.6) before (green) and after (red) correction for attenuation effects together with the best-fitting theoretical Brune's spectra (blue). This is calculated from seismograms of the HAG station located 26 km from the earthquake epicenter. Also, the measured M o , f c , absolute value of spectral decay slope at high frequencies (P) and moment magnitude M w are marked. The value of corner frequency is marked by red and green lines Vol. 179, (2022) Discrimination Between Small Earthquakes and Quarry Blasts in Egypt (Havskov & Ottemöller, 2010). This research focuses on short distances; S-waves can be handled as P-waves. Correction of the anelastic attenuation that includes attenuation due to path effect exp[-pft/ Q(f)], with the frequency-dependent quality factors (Q P and Q S ) of P-and S-waves and the attenuation near the site exp [-pf j], with (j) parameters (Singh et al., 1982), was performed. The following rela- The arrivals of P-wave, S-wave and epicentral distance values D of the recorded station (HAG) are also shown. The corrected P-and S-wave displacement spectra after data processing were calculated from seismograms of the HAG station located 35 km from the earthquake and 22 km from the quarry blast epicenter. The analyzed time windows of P-and S-waves are marked by highlighted green rectangles. The best-fitting theoretical Brune's source model is marked in blue. Also, the measured M o , f c , absolute value of spectral decay slope at high frequencies (P) and moment magnitude M w are shown. The red line shows the value of the corner frequency f c relationships for Egypt. Therefore, the findings of Abercrombie and Leary (1993) will be considered with Q P = 2Q S near the surface. Anderson et al. (1965) obtained similar results in a study of Q in the upper mantle where Q P = 2.25 Q S . Furthermore, Barton (2007) calculated Q P independently from the empirical equation Q P = (V P /V S )ÁQ S . Adopting the above results, we assume the relation Q P (f) = 2QS(f) for the frequency-dependent P-wave quality factor in this study. In Egypt, the kappa (j) factor estimated by previous studies varied from 0.01 to 0.03 s (Abou Elenean et al., 2010). For the small earthquakes with very high-frequency content, the near surface correction (j) enabled recognizing the masked corner frequencies in the high-frequency spectral decay. Figure 6 illustrates the effect of attenuation correction on displacement spectra together with the theoretical Brune's model. Because the exponent of the frequency f in Q(f) relations used in this study is close to unity, the path-dependent attenuation Flowchart of procedure adopted in EQK_SRC_PARA software for estimation of displacement spectra and source parameters for studied events depends slightly on frequency. In this case, the Q correction practically has a slight change in the shape of the P-and S-wave spectra. It affects only the low-frequency spectral level, affecting the scalar moment (M o and M w ) estimation. On the other hand, the attenuation correction for the small earthquakes with significant high-frequency content causes slight changes in the low spectral level. This affects the spectrum's shape, causing a shift in the corner frequency f c . Moreover, the near-surface attenuation correction is also essential for obtaining the true values of corner frequency f c for small earthquakes.
Once the corrected displacement spectra were obtained, the nonlinear least-square inversion technique was performed to fit Brune's source model using EQK_SRC_PARA software to estimate the spectral parameters, including the low-frequency spectral amplitude (X 0 ) and the corner frequency (f c ) at the source according to Eq. (1). Consequently, EQK_SRC_PARA software automatically calculates the source parameters for both studied earthquakes and explosions, including scalar moments (M o and M w ) based on the Eqs. (2) and (3). The fitted theoretical Brune's curve of the computed displacement spectra is defined by two lines, which refer to the low-frequency spectral plateau and high-frequency decay of f -2 . The calculated corner frequency is a boundary point in the frequency response at which spectral energy flowing through the signal begins to be reduced (attenuated) rather than passing through (at which the spectrum begins to decay at a higher frequency). It is also defined as the intersection point of the low-and high-frequency asymptotes (e.g., Baumbach & Borman, 2012;Brune, 1970;Havskov & Ottemöller, 2010;Kiratzi & Louvari, 2001;Motlagh & Mostafazadeh, 2008). Figure 7 shows an example of the corrected displacement spectra of an earthquake and an explosion in the Northern region. Additional examples are available in the supplementary materials.
A flow chart of the procedure used to estimate the displacement spectra and source parameters for the studied events is given in Fig. 8.
The seismic moments (M o (P) and M o (S)) of the investigated earthquakes and explosions for each event-station pair are estimated (Eq. 2) according to Brune's source model for both P-and S-waves separately (Brune, 1970;Hanks & Wyss, 1972): where V is the velocity (V P or V S ) extracted from Table 1, according to the source depth, X 0 is the lowfrequency spectral level, R is the epicentral distance, q is the density,<hu is the radiation pattern, and F is the free surface effect (assumed to be 2) calculated for the individual stations based on V P /V S and the emergence angles of P-and S-waves using the free surf program of Snoke (2003). The average values of radiation pattern coefficient < p ðh; uÞ = 0.64 (Baumbach & Borman, 2012) for P-waves and < s ðh; uÞ = 0.55 for S-waves (Boore & Boatwright, 1984) are used in this study. According to the concept of Hanks and Kanamori (1979), the moment magnitude M w of the examined earthquakes and explosions  at the source is calculated using the following equation: The average value for corner frequencies f c , seismic moments M o and moment magnitudes M w for each event, and their standard deviation, is estimated from the spectra for P-and S-waves recorded by different stations (i). Thus, a more accurate evaluation of these parameters is obtained. Averaging the source parameters reduces the scattering induced by different station azimuths, particularly for P-waves. Following the equations suggested by Archuleta et al. (1982), average parameters are estimated for both P-and S-waves separately: SD log x avg where x points to f c (P, S) or M o (P, S) or M w (P, S), and x avg indicates the average value of x. N is the number of the stations used; SD log x avg Â Ã represents the standard deviation of the logarithm.

Seismic Moment Vs. Corner Frequency
The average f c , M o and M w values estimated from the multiple available stations for every earthquake and explosion are used in the comparative discrimination analysis for P-and S-waves separately. The standard deviations of source parameters (f c , M o and M w ) are estimated for each event to obtain a reliable discrimination analysis. The maximum number of stations used in these estimations is 8 and 6 for each earthquake and explosion, respectively. Table 2 shows the calculated values of standard deviations of the estimated corner frequencies f c , seismic moment M o and moment magnitudes M w from Pand S-waves separately in logarithmic units for all studied earthquakes and explosions in Egypt. The average standard deviations values r of corner frequency estimations were acceptable and practically validated for both earthquakes and explosions, as required for the comparative discrimination analysis confidence (Fig. 9). For both earthquakes and explosions in Northern, Central and Southern Egypt, seismic moments M o vs. corner frequencies f c were compared individually for P-and S-waves (Fig. 10). This figure demonstrates that the seismic moments M o of earthquakes and explosions are almost located within the same range. The corner frequencies of the S-wave displacement spectra show a better separation between earthquakes and explosions than those of the P-wave spectra. The relation between f c and M o also shows a more significant separation between earthquakes and explosions for P-and S-waves in Southern Egypt than in Northern and Central Egypt regions. Comparing the corner frequencies of earthquakes and explosions in Northern, Central and Southern Egypt for P-and S-waves indicates that explosions show markedly lower corner frequencies f c than earthquakes, particularly for S-wave displacement spectra (Fig. 11). This figure displays a minor overlap in the corner frequencies of earthquakes and explosions for both Northern and Central Egypt. In Northern Egypt, this overlap is smaller for S-waves (1.25%) than for P-waves (11.7%), while Central Egypt shows an overlap of 18% and 0.06% for P-and S-waves, respectively.
In contrast, the Southern Egypt region observed a strong complete distinction between corner frequency values of earthquakes and explosions. In all Egyptian zones for earthquakes, the corner frequency values f c (P) and f c (S) manifest wide frequency ranges while explosions display narrow frequency bands. The shift of f c (P) relative to f c (S) is smaller for earthquakes than for explosions in Egypt. Table 3 demonstrates the estimated corner frequency values for earthquakes b Figure 10 Seismic moment M o vs. corner frequency f c for studied earthquakes (red) and explosions (blue) in Northern, Central and Southern Egypt obtained from P-waves (a, c, e) and S-waves (b, d, f and explosions for P-and S-waves in Northern, Central and Southern Egypt, respectively.

Corner Frequency of P-Waves Vs. Corner Frequency of S-Waves
Ataeva et al. (2017) showed that the average ratio f c (P)/f c (S) is higher for explosions than earthquakes. This has highlighted the ability to use this approach as a discriminating method. Figure (Table 4). Based on our findings presented in Fig. 12, the average f c (P) is larger than the average f c (S) for both events, but the average ratio f c (P)/f c (S) is larger for explosions. Therefore, this result proved the ability of using f c (P)/f c (S) ratio as a discriminant tool between earthquakes and explosions. A relation between f c (P)/ f c (S) ratio and moment magnitude M w estimated from P-and S-waves are constructed for all three regions b Figure 11 Histograms showing distributions of the corner frequency for studied earthquake (red) and explosions (blue cross) in Northern, Central and Southern Egypt obtained from P-wave (a, c, e) and S-wave spectra (b, d, f) Table 3 Values of corner frequencies f c (Hz) of the studied earthquakes and explosions in Northern, Central and Southern Egypt for P-and S-waves, separately  in Egypt (Fig. 13). The ratios f c (P)/f c (S) are independent of estimated moment magnitude M w and range from 1.04 to 1.51 (average 1.28) for earthquakes, while explosions display a range between 1.52 and 2.84 (average 1.97). The f c (P)/f c (S) discrimination threshold value for segregating of quarry explosions from earthquakes in Egypt is estimated to be around * 1.51-1.52 (Fig. 13). The empirical threshold value is derived by averaging the highest earthquake ratio f c (P)/f c (S) and the lowest explosion ratio f c (P)/f c (S). Furthermore, the analysis shows a broad dispersion in the corner frequency ratios of explosions.

Discussion and Conclusions
In this study, the spectral analysis technique is used to differentiate between earthquakes and quarry blasts in three regions (Northern, Central and Southern) Egypt. The most common approach in spectral analysis is to analyze and compare the corrected P-and S-waves displacement spectra for earthquakes and explosions with comparable magnitudes in each Egyptian region. This study manifested a significant difference between earthquakes and explosions in their spectral properties, focusing on the displacement spectra corner frequency that provided better information about the type of seismic source. Table 4 Obtained average f c (P)/f c (S) ratios with their standard deviations r and determination coefficient R 2 for studied earthquakes and explosions in Northern, Central and Southern Egypt  The waveform data processing involved instrumental correction, windowing of selected P and S waves of time histories, integration, calculating the displacement Fourier spectrum and correcting the path effect. Two principal spectral parameters are then extracted from the corrected spectra: the lowfrequency spectral level X 0 and the corner frequency f c . Based on Brune's source model, these two parameters are identified using nonlinear least square inversion techniques. The source parameters (f c , M o and M w ) are estimated with their standard deviations averaged over several recording stations for the same event. The obtained relatively small values of standard deviations of the average source parameters for earthquakes and explosions suggested that our approach is valid. Then, a quantitative and comparative discrimination analysis was performed for Pand S-waves separately to establish some criteria for distinguishing between earthquakes and explosions in each of the three Egyptian regions.
Our findings showed that the displacement spectra of the quarry blasts demonstrated the lack of highfrequency energy, especially in S-waves. Moreover, we observed a steeper fall off (decrease more sharply with increasing frequency) of the quarry blasts spectra than those of earthquakes spectra that showed a more stable pattern in the three regions. Therefore, the apparent corner frequencies of the P-wave and especially the quarry blasts' S-wave are lower than those for earthquakes of comparable magnitudes. Many researchers also observed similar findings in quarry blast spectra (Bennett & Murphy, 1986;Gitterman, et al., 1998;Shearer & Allmann, 2007;Taylor et al., 1988). In addition, Su et al. (1991), Aki (1995) and Yavuz et al. (2019) also observed that the quarry blast spectra decreased more sharply with increasing frequency than earthquake spectra, leading to a lower corner frequency estimate. Allmann et al. (2008) analyzed and compared the P-wave displacement spectra of different events in Southern California. They obtained a lower corner frequency of quarry blast than earthquake spectra. A steeper fall of the spectra was also observed for quarry blasts.
In this study, our findings showed that the highfrequency P/S spectral amplitude ratio resulting from the difference in corner frequencies for both earthquakes and explosions acted as an effective discriminator between earthquakes and quarry blasts. Ataeva et al. (2017) compared the displacement spectra of the P/S ratio for earthquakes and explosions and obtained the same findings. However, similar findings were also observed for nuclear explosion monitoring; the spectral ratio of regional phases (i.e., Pn/Lg, Pn/Sn, Pg/Lg, Pg/Sn) was an effective discriminant at high frequencies. Many studies have shown that regional spectral amplitude ratios effectively discriminate between earthquakes and nuclear explosions, typically at higher frequencies above about 3 to 4 Hz for test sites in China and northern Russia (e.g., Fisk, 2006;Xie & Patton, 1999). Results obtained from these studies indicated that the spectral amplitude ratios Pn/Lg, Pn/Sn, Pg/Lg and Pg/Sn manifested a good separation between earthquakes and explosions at higher frequencies at all of the hard-rock test sites, which was explained in terms of the difference in the corner frequencies of P phases and Sn-or Lg-waves for nuclear explosions.
For each region in Egypt, we constructed discrimination relations M o (P,S) vs. f c (P,S) and f c (P) vs. f c (P,S). The first relation manifested that earthquakes and explosions are clearly isolated by the S-wave compared with the P-wave, and the calculated corner frequencies of explosions are lower than those of earthquakes for both types of waves. As a result, M o vs. f c can be used as a discriminant graph, depending on the corner frequency. Furthermore, comparing the three Egyptian regions, the Southern region gives a strong and clear separation of corner frequency values of earthquakes and explosions for both waves, in contrast to the other two regions, which showed a slight interference in the corner frequency distributions. However, the variation in M o vs. f c among the three Egyptian regions could be attributed to various near-surface geologies, particularly in explosion-prone areas. The near-source (explosion) geological properties caused a change in the spectral characteristics (spectral slope and corner frequency) (Ford et al., 2011;Walter et al., 1995). The different blasting techniques and practices used in the various regions, including delay in ripple-firing (delay time) and blast hole arrangement, also caused changes in the spectral characteristics, particularly the corner frequency. Vol. 179, (2022) Discrimination Between Small Earthquakes and Quarry Blasts in Egypt 615