Resolving the location of small intracontinental earthquakes using Open Access seismic and geodetic data: lessons from the 18 January 2017 mb 4.3, Niger, earthquake

A low-magnitude earthquake was recorded on January 18, 2017, in the T\'{e}n\'{e}r\'{e} desert in Niger. This intraplate region is exceptionally sparsely covered with seismic stations and the closest open seismic station, G.TAM in Algeria at a distance of approximately 600 km, was unusually and unfortunately not operational at the time of the event. Body-wave magnitude estimates range from $m_b 4.2$ to $m_b 4.6$ and both seismic location and magnitude constraints are dominated by stations at teleseismic distances. The seismic constraints are strengthened considerably by array stations of the International Monitoring System for verifying compliance with the Comprehensive Nuclear Test-Ban-Treaty. This event, with magnitude relevant to low-yield nuclear tests, provides a valuable validation of the detection and location procedure for small land-based seismic disturbances at significant distances. For seismologists not in the CTBT system, the event is problematic as data from many of the key stations are not openly available. We examine the uncertainty in published routinely-determined epicenters by performing multiple Bayesloc location estimates with published arrival times considering both all published arrival times and those from open stations only. This location exercise confirms lateral uncertainties in seismologically-derived location no smaller than 10 km. Coherence for InSAR in this region is exceptionally high, and allows us to confidently detect a displacement of the order 6 mm in the time-frame containing the earthquake, consistent with the seismic location estimates, and with a lateral length scale consistent with an earthquake of this size, allowing location constraint to within one rupture length ($\leq 5$ km) -- significantly reducing the lateral uncertainty compared with relying on seismological data only. Combining Open Access-only seismological and geodetic data, we precisely constrain the source location, and conclude that this earthquake likely had a shallow source. We then discuss potential ways to continue the integration of geodetic data in the calibration of seismological earthquake location.

SEISMICITY of the ARCTIC in 2015

The article provides an overview and analysis of seismicity within the boundaries of the Arctic region for 2015, a description of seismic station networks, and processing methods. The catalog of earthquakes in the Arctic region was compiled on the basis of catalogs of several organizations and seismological centers. In total, 334 earthquakes are included in the earthquake catalog. Most of the earthquakes that occurred in 2015, including all the strongest earthquakes, were located within the mid-ocean ridges of Mon, Knipovich and Gakkel. In the offshore territories, most of the earthquakes were confined to the Svalbard archipelago, in particular, to the seismically active zone in the Sturfjord strait. The renewal of instrumental seismological observations in 2011 (station ZFI) on Alexandra Land Island in the Franz Josef Land archipelago made it possible to record weak earthquakes in the north of the shelf of the Barents and Kara Seas. For twelve earthquakes, the focal mechanism parameters are presented according to the Global CMT catalog.

SEISMICITY of NORTHERN EURASIA in 2015

The review of the Northern Eurasia seismicity for 2015 includes a description of seismic networks, the results of analysis of the seismic regime and individual noticeable earthquakes in 16 regions of Russia and neighbouring countries. Seismic monitoring was carried out by the networks of seismic station of Russia, Azerbaijan, Armenia, Belarus, Kazakhstan, Kyrgyzstan, Latvia, Moldova, Turkmenistan, Tajikistan, Uzbekistan, Ukraine, including 599 digital, 7 analogue stations and eight seismic groups. In 2015, these networks registered about 27 thousand tectonic earthquakes, over 6 thousand volcanic earthquakes, 599 explosions, 23 mountain-tectonic shocks and induced earthquakes. Focal mechanisms of 592 earthquakes were determined, the information on manifestations of 449 perceptible earthquakes was collected. 26 shocks were felt in settlements of Northern Eurasia with an intensity Ii≥5. According to estimates of the annual number and released seismic energy in 2015 in comparison with the long-term characteristics of the seismic regime, the seismic process in most regions of Northern Eurasia proceeded in the “background” regime. An exception is Tajikistan and adjacent territories, where two strong earthquakes occurred – the Hindu Kush earthquake on October 26 with Mw=7.5, h=230 km in northern Afghanistan, near the border with Tajikistan, and the Sarez earthquake on December 7 with Mw=7.2, Ms=7.6, h=20 km in Tajikistan. Both earthquakes were accompanied by numerous aftershocks and were felt in Tajikistan with intensities Imax=7 and Imax=7–8 respectively, on the MSK-64 scale. Notable event on the territory of Northern Eurasia in 2015 is the emergence of the Muyakan sequence of earthquakes, the largest for the period of instrumental observations in the region "Baikal and Transbaikalia", as a result of which the number of recorded earthquakes in the region quadrupled concerning 2014. The other interesting fact is occurrence of tangible earthquakes in the regions, traditionally considered weakly seismic – near the Semipalatinsk test area in Eastern Kazakhstan (Chingiz earthquake on January 20, Ms=4.1, I0=5–6), in the Middle Urals (Middle Ural earthquake on October 18 with ML=4.7, I0=6) and in the southwest of East -European platform (Poltava earthquake on February 3 with KR=10.7, I0=6).

MUYAKAN EARTHQUAKE SEQUENCE in 2015 (Northern Baikal Region)

We present the results of studies of the seismic regime, focal mechanisms, and macroseismic data in the area of the largest Muyakan activation in 2015 (northern Baikal region). Due to the deployment of a network of temporary seismic stations, the number of registered earthquakes (KR≥3) increased significantly in 2015 and reached  30 thousand. Spatio-temporal development of the considered activation is characterized by stable low values of earthquake hypocenters and dividing the epicentral field into two clusters – eastern and north-western ones. Both clusters are connected with local stress-strain field (rift type and strike-slip, respectively), while the general regime of seismotectonic deformations of the crust in the activation area, calculated from the statistical analysis of 77 focal mechanisms of Muyakan earthquakes (KR≥9.2), demonstrates the predominance of submeri-dional horizontal extension. Macroseismic effects from the largest earthquakes of the Muyakan sequence were felt, predominantly, in Severomuysk settlement (=10–15 km, I=5). New data on strong motions, obtained from the records of the seismic station with the same name, have significantly complemented the database for the territory of the north-eastern flank of the Baikal rift zone. In general, the obtained results could be used to clarify the seismic hazard of the considered area and to improve the instrumental part of seismic scales.

Analysis of the possibility for seismic early warning in the Balkan Peninsula

This study analyses the possibility for seismic early warning (EW) in the Balkan Peninsula. A number of characteristics of seismic record were evaluated for EW utility. Some tests checked the possibility to locate events reliably by Golitsyn’s method, using one seismic station (SS). The distance and relevant travel time from some crustal faults to the nearest SS and big towns were appraised. EW procedures for most of the seismic zones in the Balkan Peninsula are not reliable, excepting the Vrancea zone, because of the small density of the SS, crustal depth of the earthquakes and fault crowdedness of the region.

The mud volcanoes at Santa Barbara and Aragona (Sicily, Italy): a contribution to risk assessment

The Santa Barbara and Aragona areas are affected by mud volcanism (MV) phenomena, consisting of continuous or intermittent emission of mud, water, and gases. This activity could be interrupted by paroxysmal events, with an eruptive column composed mainly of clay material, water, and gases. They are the most hazardous phenomena, and today it is impossible to define the potential parameters for modelling the phenomenon. In 2017, two digital surface models (DSMs) were performed by drone in both areas, thus allowing the mapping of the emission zones and the covered areas by the previous events. Detailed information about past paroxysms was obtained from historical sources, and, with the analysis of the 2017 DSMs, a preliminary hazard assessment was carried out for the first time at two sites. Two potentially hazardous paroxysm surfaces of 0.12 and 0.20 km2 for Santa Barbara and Aragona respectively were defined. In May 2020, at Aragona, a new paroxysm covered a surface of 8721 m2. After this, a new detailed DSM was collected with the aim to make a comparison with the 2017 one. Since 2017, a seismic station was installed in Santa Barbara. From preliminary results, both seismic events and ambient noise showed a frequency of 5–10 Hz.

An Evaluation of the Timing Accuracy of Global and Regional Seismic Stations and Networks

Clock accuracy is a basic parameter of any seismic station and has become increasingly important for seismology as the community seeks to refine structures and dynamic processes of the Earth. In this study, we measure the arrival time differences of moderate repeating earthquakes with magnitude 5.0–5.9 in the time range of 1991–2017 at the same seismic stations by cross-correlating their highly similar waveforms and thereby identify potential timing errors from the outliers of the measurements. The method has very high precision of about 10 ms and shows great potential to be used for routine inspection of the timing accuracy of historical and future digital seismic data. Here, we report 5131 probable cases of timing errors from 451 global and regional stations available from the Incorporated Research Institutions for Seismology Data Management Center, ranging from several tens of milliseconds to over 10 s. Clock accuracy seems to be a prevailing problem in permanent stations with long-running histories. Although most of the timing errors have already been tagged with low timing quality, there are quite a few exceptions, which call for greater attention from network operators and the seismological community. In addition, seismic studies, especially those on temporal changes of the Earth’s media from absolute arrival times, should be careful to avoid misinterpreting timing errors as temporal changes, which is indeed a problem in some previous studies of the Earth’s inner core boundary.

Applicability of On-Site P-Wave Earthquake Early Warning to Seismic Data Observed During the 2011 Off the Pacific Coast of Tohoku Earthquake, Japan

In this study, the on-site P-wave earthquake early warning (EEW) based on the site-specific spectral ratio of S-wave to P-wave to efficiently incorporate the site characteristics, which can potentially issue the earthquake warning by the time of Ts-p, was developed. The spectral ratio of S-wave to P-wave that are related to the source effects, the path effects, and the site effects are significantly affected by the site effects contrast to the source effects and the path effects in practical. At first, the on-site P-wave EEW method which multiplies a site-specific spectral ratio of S-wave to P-wave prepared in advance by P-wave observed in the real-time at seismic stations is applied to seismic data for moderate-sized earthquakes with a magnitude (Mj) of 5.0–6.0, occurred in the eastern Japan, observed at both the sedimentary basin site and the rock site. As a result, this method predicted well the observed S-wave in the single indicator of SI within the logarithmic standard deviation of 0.25 as well as in the frequency of more than 0.5 Hz. It is, also, confirmed that the site-specific spectral ratio of S-wave to P-wave at a seismic station was stably retrieved from 20 data samples at least. To investigate the applicability of this method to earthquake ground motions induced by a large-scaled earthquake, finally, this method is applied to seismic data during the 2011 off the Pacific coast of Tohoku earthquake, Japan (Mw 9.0). The prediction of S-wave using a time-window of 10 s after P-wave arrived, could not reproduce the observation with the underestimation; however, the prediction of S-wave using a time-window of more than 20 s containing P-wave propagated from an area generating strong motions in the fault, could reproduce the observation. Even in the case of the large-scaled earthquake, the on-site P-wave EEW method based on the site-specific spectral ratio of S-wave to P-wave at a seismic station availably works by using the gradually increasing time-windows after P-wave arrived in the single indicator of SI as well as in the frequency content, avoiding the mixture of S-wave into a part of P-wave.

Wind-induced seismic noise at the Princess Elisabeth Antarctica Station

Icequakes are the result of processes occurring within the ice mass or between the ice and its environment. Studying icequakes provides a unique view on ice dynamics, specifically on the basal conditions. Changes in conditions due to environmental or climate changes are reflected in icequakes. Counting and characterizing icequakes is thus essential to monitor them. Most of the icequakes recorded by the seismic station at the Belgian Princess Elisabeth Antarctica Station (PE) have small amplitudes corresponding to maximal displacements of a few nanometres. Their detection threshold is highly variable because of the rapid and strong changes in the local seismic noise level. Therefore, we evaluated the influence of katabatic winds on the noise measured by the well-protected PE surface seismometer. Our purpose is to identify whether the lack of icequake detection during some periods could be associated with variations in the processes generating them or simply with a stronger seismic noise linked to stronger wind conditions. We observed that the wind mainly influences seismic noise at frequencies greater than 1 Hz. The seismic noise power exhibits a bilinear correlation with the wind velocity, with two different slopes at a wind velocity lower and greater than 6 m s−1 and with, for example at a period of 0.26 s, a respective variation of 0.4 dB (m −1 s) and 1.4 dB (m −1 s). These results allowed a synthetic frequency and wind-speed-dependent noise model to be presented that explains the behaviour of the wind-induced seismic noise at PE, which shows that seismic noise amplitude increases exponentially with increasing wind speed. This model enables us to study the influence of the wind on the original seismic dataset, which improves the observation of cryoseismic activity near the PE station.

Unexpected Shallow Earthquake of August 1st, 2020 in the North of Indramayu, West Java, Indonesia

During the last one hundred years, there are no shallow seismicity in the north of Java. This area is dominated by intermediate and deep focus earthquakes due to the subducted Indo-Australian slab. An earthquake with magnitude ML 4.5 struck Indramayu, north of West Java on August 1, 2020. According to the Agency for Meteorology, Climatology, and Geophysics (BMKG), the earthquake was felt III MMI scale in Indramayu and its vicinity. We used waveform data from BMKG seismic station in West Java, then we picked P-and S-waves arrival times from each station and hypocenter location was determined by Geiger method. We have detected Pn before Pg phase on four BMKG seismic stations, indicating a shallow crustal earthquake. Our inversion show that the earthquake occurred in 6.1805° S, 108.2612° E with 5 km focus depth at 16:24:38 GMT+7. Our focal mechanism solution was determined by using moment tensor inversion shows a strike-slip faulting, which corresponds to the active fault in the north of Indramayu.