Earthquake hazard and risk assessment based on Unified Scaling Law for Earthquakes: Greater Caucasus and Crimea

We continue applying the general concept of seismic risk analysis in a number of seismic regions worldwide by constructing regional seismic hazard maps based on morphostructural analysis, pattern recognition, and the Unified Scaling Law for Earthquakes (USLE), which generalizes the Gutenberg-Richter relationship making use of naturally fractal distribution of earthquake sources of different size in a seismic region. The USLE stands for an empirical relationship log₁₀N(M, L)?=?A?+?B·(5 – M)?+?C·log₁₀L, where N(M, L) is the expected annual number of earthquakes of a certain magnitude M within a seismically prone area of linear dimension L. We use parameters A, B, and C of USLE to estimate, first, the expected maximum magnitude in a time interval at seismically prone nodes of the morphostructural scheme of the region under study, then map the corresponding expected ground shaking parameters (e.g., peak ground acceleration, PGA, or macro-seismic intensity). After a rigorous verification against the available seismic evidences in the past (usually, the observed instrumental PGA or the historically reported macro-seismic intensity), such a seismic hazard map is used to generate maps of specific earthquake risks for population, cities, and infrastructures (e.g., those based on census of population, buildings inventory). The methodology of seismic hazard and risk assessment is illustrated by application to the territory of Greater Caucasus and Crimea.     

Some Probabilistic and Statistical Properties of the Seismic Regime of Zemmouri (Algeria) Seismoactive Zone

Statistical tests have been used to adjust the Zemmouri seismic data using a distribution function. The Pareto law has been used and the probabilities of various expected earthquakes were computed. A mathematical expression giving the quantiles was established. The extreme values limiting law confirmed the accuracy of the adjustment method. Using the moment magnitude scale, a probabilistic model was made to predict the occurrences of strong earthquakes. The seismic structure has been characterized by the slope of the recurrence plot γ, fractal dimension D, concentration parameter K_sr, Hurst exponents H_r and H_t. The values of D, γ, K_sr, H_r, and H_t diminished many months before the principal seismic shock (M = 6.9) of the studied seismoactive zone has occurred. Three stages of the deformation of the geophysical medium are manifested in the variation of the coefficient G% of the clustering of minor seismic events.     

Detailed Study of Time Variations in the Gutenberg–Richter <Emphasis Type="Italic">b</Emphasis>-value Based on Highly Accurate Seismic Observations at the Garm Prognostic Site,Tajikistan

The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log₁₀E^2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log₁₀E^2/3, and log₁₀N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log₁₀E^2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log₁₀N_i and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log₁₀N_i time series in the different magnitude ranges (M ≥ Mc_i) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔM_i) in the range of low magnitudes are known. It is concluded that using the parameter log₁₀N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.     

Geodynamic and Seismotectonic Activity in Eastern Tibet in the 21st Century

At the beginning of the 21st century, a series of great earthquakes were recorded in northeastern Tibet, along the periphery of the Bayan Hara lithospheric block. An earthquake with M_S = 8.1 occurred within the East Kunlun fault zone in the Kunlun Mountains, which caused an extended surface rupture with left-lateral strike slip. An earthquake with M_S = 8 occurred in Wenchuan (China) on May 12, 2008, giving rise to an extended overthrust along the Lunmanshan fault zone. An earthquake with M_S = 7.1 occurred in Yushu (China) on April 14, 2010; its epicenter was on the Grazze–Yushu–Funchuoshan fault; a left-lateral strikeslip offset was observed on the surface. An earthquake with M_S = 7 occurred in the vicinity of Lushan on April 20, 2013; its epicenter was within the Lunmanshan fault zone, 103 km southwest of the zone of the catastrophic Wenchuan earthquake. An earthquake with M_S = 8.2 occurred in Nepal on April 25, 2015. Based on the CSN seismic catalog, the energy of all earthquakes in eastern Tibet at the end of the 20th and beginning of the 21st centuries was estimated. It was found that Tibet was seismically quiet from 1980 to 2000. The beginning of the 21st century has been marked by seismic activation with earthquake sources migrating southward to surround the Bayan Hara lithospheric block from every quarter. Therefore, this block can be regarded as one of the most seismically active regions of China.     

Features of Radiation and Propagation of Seismic Waves in the Northern Caucasus: Manifestations in Regional Coda

Based on the Anapa (ANN) seismic station records of ~40 earthquakes (M_W > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.     

Precursory seismic quiescence along the Sumatra-Andaman subduction zone: past and present

In this study, the seismic quiescence prior to hazardous earthquakes was analyzed along the Sumatra-Andaman subduction zone (SASZ). The seismicity data were screened statistically with mainshock earthquakes of M _w?≥?4.4 reported during 1980–2015 being defined as the completeness database. In order to examine the possibility of using the seismic quiescence stage as a marker of subsequent earthquakes, the seismicity data reported prior to the eight major earthquakes along the SASZ were analyzed for changes in their seismicity rate using the statistical Z test. Iterative tests revealed that Z factors of N?=?50 events and T?=?2?years were optimal for detecting sudden rate changes such as quiescence and to map these spatially. The observed quiescence periods conformed to the subsequent major earthquake occurrences both spatially and temporally. Using suitable conditions obtained from successive retrospective tests, the seismicity rate changes were then mapped from the most up-to-date seismicity data available. This revealed three areas along the SASZ that might generate a major earthquake in the future: (i) Nicobar Islands (Z?=?6.7), (ii) the western offshore side of Sumatra Island (Z?=?7.1), and (iii) western Myanmar (Z?=?6.7). The performance of a stochastic test using a number of synthetic randomized catalogues indicated these levels of anomalous Z value showed the above anomaly is unlikely due to chance or random fluctuations of the earthquake. Thus, these three areas have a high possibility of generating a strong-to-major earthquake in the future.     

Site Specific Ground Motion Modeling and Seismic Response Analysis for Microzonation of Baku,Azerbaijan

We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V _P values measured for the soils. The downtown part of Baku city is characterized by low V_S30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

The <Emphasis Type="Italic">Q</Emphasis>-factor estimates for the crust and upper mantle in the vicinity of Sochi and Anapa (North Caucasus)

The regularities in the radiation and propagation of seismic waves in the regions of the North Caucasus are analyzed for estimating the ground motion parameters during the probable future strong earthquakes. Based on the records of the regional earthquakes with magnitudes M_W ~ 3.9–5.6 within epicentral distances up to ~300 km obtained during the period of digital measurements at the Sochi and Anapa seismic stations, the Q-factors in the vicinities of these sites are estimated at ~55 f^0.9 and ~90f^0.7, respectively. The estimates were obtained by the coda normalization method developed by Aki, Rautian, and other authors. This method is based on the phenomenon of suppression of the earthquake (source) effects and local (site) responses by coda waves in the S-wave spectra. The obtained Q-factor estimates can be used for forecasting the ground shaking parameters for the future probable strong earthquakes in the North Caucasus in the vicinities of Sochi and Anapa.     

The ratio between corner frequencies of source spectra of <Emphasis Type="Italic">P</Emphasis>- and <Emphasis Type="Italic">S</Emphasis>-waves—a new discriminant between earthquakes and quarry blasts

This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f ₀(P)/f ₀(S). The ratios were computed for each event with corner frequencies f ₀ of P- and S-wave, which were obtained from the measured f ₀ ^I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f ₀(P)/f ₀(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M _w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f ₀(P)/f ₀(S), with the empirical threshold value of the ratio for Israel as 1.48.     

Temporal and spatial distributions of precursory seismicity rate changes in the Thailand-Laos-Myanmar border region: implication for upcoming hazardous earthquakes

To study the prospective areas of upcoming strong-to-major earthquakes, i.e., M _w  ≥ 6.0, a catalog of seismicity in the vicinity of the Thailand-Laos-Myanmar border region was generated and then investigated statistically. Based on the successful investigations of previous works, the seismicity rate change (Z value) technique was applied in this study. According to the completeness earthquake dataset, eight available case studies of strong-to-major earthquakes were investigated retrospectively. After iterative tests of the characteristic parameters concerning the number of earthquakes (N) and time window (T _w ), the values of 50 and 1.2 years, respectively, were found to reveal an anomalous high Z-value peak (seismic quiescence) prior to the occurrence of six out of the eight major earthquake events studied. In addition, the location of the Z-value anomalies conformed fairly well to the epicenters of those earthquakes. Based on the investigation of correlation coefficient and the stochastic test of the Z values, the parameters used here (N = 50 events and T _w  = 1.2 years) were suitable to determine the precursory Z value and not random phenomena. The Z values of this study and the frequency-magnitude distribution b values of a previous work both highlighted the same prospective areas that might generate an upcoming major earthquake: (i) some areas in the northern part of Laos and (ii) the eastern part of Myanmar.     

<Emphasis Type="Italic">M</Emphasis>9 Tohoku Earthquake Hydro- and Seismic Response in the Caucasus and North Turkey

Presently, there are a lot of observations on the significant impact of strong remote earthquakes on underground water and local seismicity. Teleseismic wave trains of strong earthquakes give rise to several hydraulic effects in boreholes, namely permanent water level changes and water level oscillations, which closely mimic the seismograms (hydro-seismograms). Clear identical anomalies in the deep borehole water levels have been observed on a large part of the territory of Georgia during passing of the S and Love–Rayleigh teleseismic waves (including also multiple surface Rayleigh waves) of the 2011 Tohoku M9 earthquake. The analysis carried out in order to find dynamically triggered events (non-volcanic tremors) of the Tohoku earthquake by the accepted methodology has not revealed a clear tremor signature in the test area: the Caucasus and North Turkey. The possible mechanisms of some seismic signals of unknown origin observed during passage of teleseismic waves of Tohoku earthquake are discussed.     

Possible Short-Term Precursors of Strong Crustal Earthquakes in Japan based on Data from the Ground Stations of Vertical Ionospheric Sounding

We have studied changes in the ionosphere prior to strong crustal earthquakes with magnitudes of М ≥ 6.5 based on the data from the ground-based stations of vertical ionospheric sounding Kokobunji, Akita, and Wakkanai for the period 1968–2004. The data are analyzed based on hourly measurements of the virtual height and frequency parameters of the sporadic E layer and critical frequency of the regular F2 layer over the course of three days prior to the earthquakes. In the studied intervals of time before all earthquakes, anomalous changes were discovered both in the frequency parameters of the Es and F2 ionospheric layers and in the virtual height of the sporadic E layer; the changes were observed on the same day at stations spaced apart by several hundred kilometers. A high degree of correlation is found between the lead-time of these ionospheric anomalies preceding the seismic impact and the magnitude of the subsequent earthquakes. It is concluded that such ionospheric disturbances can be short-term ionospheric precursors of earthquakes.     

On the Problem of Deep Earthquakes in Crimea–Black Sea Region

It is a common opinion that only crustal earthquakes can occur in the Crimea–Black Sea region. Since the existence of deep earthquakes in the Crimea–Black Sea region is extremely important for the construction of a geodynamic model for this region, an attempt is made to verify the validity of this widespread view. To do this, the coordinates of all earthquakes recorded by the stations of the Crimean seismological network are reinterpreted with an algorithm developed by one of the authors. The data published in the seismological catalogs and bulletins of the Crimea–Black Sea region for 1970–2012 are used for the analysis. To refine the coordinates of hypocenters of earthquakes in the Crimea–Black Sea region, in addition to the data from stations of the Crimean seismological network, information from seismic stations located around the Black Sea coast are used. In total, the data from 61 seismic stations were used to determine the hypocenter coordinates. The used earthquake catalogs for 1970–2012 contain information on ~2140 events with magnitudes from–1.5 to 5.5. The bulletins provide information on the arrival times of P- and S-waves at seismic stations for 1988 events recorded by three or more stations. The principal innovation of this study is the use of the original author’s hypocenter determination algorithm, which minimizes the functional of distances between the points (X, Y, H) and (x, y, h) corresponding to the theoretical and observed seismic wave travel times from the earthquake source to the recording stations. The determination of the coordinates of earthquake hypocenters is much more stable in this case than the usual minimization of the residual functional for the arrival time of an earthquake wave at a station (the difference between the theoretical and observed values). Since determination of the hypocenter coordinates can be influenced by the chosen velocity column beneath each station, special attention is focused on collecting information on velocity profiles. To evaluate the influence of the upper mantle on the results of calculating the velocity model, two different low-velocity and high-velocity models are used; the results are compared with each other. Both velocity models are set to a depth of 640 km, which is fundamentally important in determining hypocenters for deep earthquakes. Studies of the Crimea–Black Sea region have revealed more than 70 earthquakes with a source depth of more than 60 km. The adequacy of the obtained depth values is confirmed by the results of comparing the initial experimental data from the bulletins with the theoretical travel-time curves for earthquake sources with depths of 50 and 200 km. The sources of deep earthquakes found in the Crimea–Black Sea region significantly change our understanding of the structure and geotectonics of this region.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

Anomalous variations in the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript>-layer structure at geomagnetic midlatitudes of the Southern and Northern hemispheres at the transition from summer to winter conditions under low solar activity

The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N₂ and O₂ molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O⁺ on N _e ions under winter and summer conditions forms not more than 11% of the N _e winter anomaly event in the F ₂ layer and topside ionosphere. The model without electronically excited O⁺ ions reduces the duration of the N _e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.     

Seismic network detection probability assessment using waveforms and accounting to event association logic

The geographical area where a seismic event of magnitude M?≥?M _t is detected by a seismic station network, for a defined probability is derived from a station probability of detection estimated as a function of epicentral distance. The latter is determined from both the bulletin data and the waveforms recorded by the station during the occurrence of the event with and without band-pass filtering. For simulating the real detection process, the waveforms are processed using the conventional Carl Johnson detection and association algorithm. The attempt is presented to account for the association time criterion in addition to the conventional approach adopted by the known PMC method.     

Propagation of internal gravity waves to the altitudes of the ionospheric <Emphasis Type="Italic">D</Emphasis> and dynamo regions in the seismic region (Kamchatka): Preliminary results

A spectral analysis of simultaneous diurnal variations in the E _z component of the quasi-static electric field in the near-Earth atmosphere, VLF radio noise, and the horizontal component of the geomagnetic field, observed at Kamchatka in September 1999, has been performed. These geophysical parameters are indirectly used to study wave processes in the near-Earth atmosphere and in the ionospheric D and dynamo regions within the band of periods of internal gravity waves (T = 0.5?3.5 h). The correlation method in the frequency region is used to analyze the interrelation between the wave processes in these atmospheric regions. The power cross-spectra of various pairs of geophysical parameters have been studied depending on meteorological, seismic, and geomagnetic activities. It is shown that the oscillations in the power spectra in the T ～ 1–1.5 h band of periods are caused by the sources of internal gravity waves in the near-Earth atmosphere and by the remote sources above the dynamo region of the ionosphere within the T ～ 1.5–3 h band of periods.     

Source parameters for the 2013–2015 earthquake sequence in Nógrád county,Hungary

Between 2013 June and 2015 January, 35 earthquakes with local magnitude M _L ranging from 1.1 to 4.2 occurred in Nógrád county, Hungary. This earthquake sequence represents above average seismic activity in the region and is the first one that was recorded by a significant number of three-component digital seismographs in the county. Using a Bayesian multiple-event location algorithm, we have estimated the hypocenters of 30 earthquakes with M _L ≥1.5. The events occurred in two small regions of a few squared kilometers: one to the east of Érsekvadkert and the other at Iliny. The uncertainty of the epicenters is about 1.5–1.7 km in the E-W direction and 1.8–2.1 km in the N-S direction at the 95 % confidence level. The estimated event depths are confined to the upper 3 km of the crust. We have successfully estimated the full moment tensors of 4 M _w ≥3.6 earthquakes using a probabilistic waveform inversion procedure. The non-double-couple components of the retrieved moment tensor solutions are statistically insignificant. The negligible amount of the isotropic component implies the tectonic nature of the investigated events. All of the analyzed earthquakes have strike-slip mechanism with either right-lateral slip on an approximately N-S striking or left-lateral movement on a roughly E-W striking nodal plane. The orientations of the obtained focal mechanisms are in good agreement with the main stress pattern published for the epicentral region. Both the P and T principal axes are horizontal, and the P axis is oriented along a NE-SW direction.