SEISMICITY of KAMCHATKA and COMMANDER ISLANDS in 2015

The seismicity review of Kamchatka and surrounding territories for 2015 is given. In the Kamchatka earthquake catalogue, the minimum local magnitude of completeness is MLmin=3.5, and for earthquakes with h≥350 km under the Okhotsk sea MLmin=3.6. The Kamchatka earthquake catalogue for 2015 with ML3.5, published in the Appendix to this issue, includes 1213 events. 92 earthquakes of the catalogue with ML=3.0–6.5 were felt in Kamchatka and surrounding areas with seismic intensity I=2–6 according to the MSK-64 scale. For all events with ML5.0 that occurred in 2015 in the KB GS RAS area of responsibility, an attempt to calculate the seismic moment tensor (SMT) was made. There are 32 such events in the regional catalogue. For 28 earthquakes the SMT and depth h of the equivalent point source were calculated successfully. The calculations were performed for the SMT double-couple model using a nonlinear algorithm. In 2015, a typical location of the earthquake epicenters was observed in the Kamchatka zone. In 2015, the seismicity level in all selected zones and in the region as a whole correspond to the background one according to the “SESL’09” scale. The number of recorded events with ML3.5 and strong earthquakes with ML5.0 is close to the average annual value. Anomalous and significant events were not recorded.

Supplemental Material: Earthquake weather and climate change: Should we stress about the forecast?

Table S1. Onshore UK earthquake catalogue for 1980–2012 with local magnitudes ≥1.5<br>

Supplemental Material: Earthquake weather and climate change: Should we stress about the forecast?

Table S1. Onshore UK earthquake catalogue for 1980–2012 with local magnitudes ≥1.5<br>

Characteristics of Foreshocks Occurrence of Onshore Earthquakes In Japan For Mj3.0 To 7.2 Mainshocks

We use the Japan Meteorological Agency (JMA) earthquake catalogue from 2001 to 2021 to investigate the spatiotemporal distribution of foreshocks for shallow mainshocks (Mj3.0–7.2) that are located onshore of Japan. We find clear peaks for the earlier small earthquakes within 10 days and 3 km prior to the mainshocks, which are considered as our definition of foreshocks. After removing the aftershocks, earthquake swarms and possible earthquakes triggered by the 2011 Mw9.0 Tohoku-oki earthquake, we find that for the 2,066 independent earthquakes, 783 (37.9%) have one or more foreshocks. There is a decreasing trend of foreshock occurrence with mainshock depth. Also, normal faulting earthquakes have higher foreshock occurrence than reverse faulting earthquakes. We calculate the rates of foreshock occurrence as a function of the magnitudes of foreshocks and mainshocks, and we have found no clear trend between the magnitudes of foreshocks and mainshocks.

A homogeneous earthquake catalogue for Turkey

A new homogenized earthquake catalogue for Turkey is compiled for the period 1900–2018. The earthquake parameters are obtained from the Bulletin of International Seismological Centre that was fully updated in 2020. New conversion equations between moment magnitude and the other scales (md, ML, mb, Ms, and M) are determined using the general orthogonal regression method to build up a homogeneous catalogue, which is the essential database for seismic hazard studies. The 95 % confidence intervals are estimated using the bootstrap method with 1000 samples. The equivalent moment magnitudes (Mw*) for the entire catalogue are calculated using the magnitude relations to homogenize the catalogue. The magnitude of completeness is 2.7 Mw*. The final catalogue is not declustered or truncated using a threshold magnitude in order to be a widely usable catalogue. It contains not only Mw* but also the average and median of the observed magnitudes for each event. Contrary to the limited earthquake parameters in the previous catalogues for Turkey, the 45 parameters of ∼378 000 events are presented in this study.

A pragmatic approach to adjusting early instrumental local magnitudes for seismic hazard assessments in Australia

Prior to the development of Australian-specific magnitude formulae, the 1935 magnitude correction factors by Charles Richter—originally developed for southern California—were almost exclusively used to calculate earthquake magnitudes throughout Australia prior to the 1990s. Due to the difference in ground-motion attenuation between southern California and much of the Australian continent, many earthquake magnitudes from the early instrumental era are likely to have been overestimated in the Australian earthquake catalogue. A method is developed that adjusts local magnitudes (ML) using the difference between the original (inappropriate) magnitude formulae (or look-up tables) and the Australian-specific formulae at a distance determined by the nearest recording station likely to have recorded the earthquake. Nationally, these adjustments have reduced the number of earthquakes of ML ≥ 4.5 in the early instrumental catalogue by approximately 25% since 1900, while the number of ML ≥ 5.0 earthquakes has reduced by approximately 32% over the same time period. The reduction in the number of moderate-to-large-magnitude earthquakes over the instrumental period yields long-term earthquake rates that are more consistent with present-day rates, since the development of Australian-specific magnitude formulae (approximately 1990). The adjustment of early instrumental magnitudes to obtain consistently derived earthquake catalogue is important for seismic hazard assessments.

The Seismicity of Ischia Island, Italy: An Integrated Earthquake Catalogue From 8th Century BC to 2019 and Its Statistical Properties

Ischia is a densely inhabited and touristic volcanic island located in the northern sector of the Gulf of Naples (Italy). In 2017, the Mw 3.9 Casamicciola earthquake occurred after more than one century of seismic quiescence characterized only by minor seismicity, which followed a century with three destructive earthquakes (in 1828, 1881, and 1883). These events, despite their moderate magnitude (Mw &lt; 5.5), lead to dreadful effects on buildings and population. However, an integrated catalogue systematically covering historical and instrumental seismicity of Ischia has been still lacking since many years. Here, we review and systematically re-analyse all the available data on the historical and instrumental seismicity, to build an integrated earthquake catalogue for Ischia with a robust characterization of existing uncertainties. Supported by new or updated macroseismic datasets, we significantly enriched existing catalogues, as the Italian Parametric Earthquake Catalogue (CPTI15) that, with this analysis, passed from 12 to 57 earthquakes with macroseismic parametrization. We also extended back by 6 years the coverage of the instrumental catalogue, homogenizing the estimated seismic parameters. The obtained catalogue will not only represent a solid base for future local hazard quantifications, but also it provides the unique opportunity of characterizing the evolution of the Ischia seismicity over centuries. To this end, we analyse the spatial, temporal, and magnitude distributions of Ischia seismicity, revealing for example that, also in the present long-lasting period of volcanic quiescence, is significantly non-stationary and characterized by a b-value larger than 1.

Do spatial and/or temporal b-value variations exist within the Groningen induced earthquake catalogue?

&lt;p&gt;We present a data-driven analysis to derive whether statistically significant spatial and/or temporal Gutenberg-Richter b-value variations exist within the induced earthquake catalogue of the Groningen gas field. We utilize the method developed by Kamer and Hiemer (2015; J. Geophys. Res. Solid Earth, 120, doi:10.1002/2014JB011510 ) which is based on optimal partitioning using Voronoi tessellation, penalized likelihood, and wisdom of the crowd philosophy. Our implementation derives both the magnitude of completeness and the b-values simultaneously. The magnitude of completeness is computed with the maximum curvature method with a correction applied to avoid bias due to catalogue incompleteness. Finally, following Marzocchi et al. (2020; Geophys. J. Int. 220, doi: 10.1093/gji/ggz541) the b-values computed are corrected for bin size and small sample sizes.&lt;/p&gt;&lt;p&gt;In a first step we have limited the analysis to spatial variations in the b-values. A significant advantage of the approach taken is that it is feasible to also derive b-values in regions of very low data density. We will show that a statistically significant variation in b-values is obtained. Very low b-values (b&lt;0.8) are observed in the central-northern part of the gas field. However, in the west near the production cluster Eemskanaal (EKL) and in the east near the city of Delfzijl significantly higher b-values (b&gt;1.1) are observed. A Kolmogorov-Smirnov test of frequency-magnitude distributions for the two areas obtains a p-value of 1.5 10-13 and 2.3 10-12 for the EKL region and Delfzijl regions, respectively, rendering the difference more than statistically significant at the 99% confidence level.&lt;/p&gt;&lt;p&gt;In a second step we extended the spatial analysis to a spatial-temporal analysis. The results of the analysis show that the Groningen earthquake database is too small to derive meaningful spatial results for the full Groningen gas field based on multiple random temporal nodes.&amp;#160; We divided the dataset in two almost equal datasets: both containing roughly 50% of the data and of comparable spatial resolution. Spatial analysis of these two subsets of the catalogue shows a significant decrease of the b-values in the central and southern regions. Particularly in the western EKL region the b-value decreases from 1.2 to 0.92. The decrease is close to significant at the 90% confidence level. The northern region exhibits comparable low b-values in both periods. As the data in the first decade is primarily concentrated in the northern region, we have attempted to assess the spatial b-value here in the period prior to 2005. We find the high b-value area is significantly smaller and the minimum value is higher (b = 0.96 pre-2005 versus b = 0.88 post-2012). The difference is significant only at the interquartile level, but the model resolution is low.&lt;/p&gt;&lt;p&gt;Based on our results, we could conclude a spatial and temporal variation in b-value is observed. However, despite our efforts to limit bias in the derivation, variations could still result from the presence of a truncation. Hence, we will extend the current analysis by a comparable analysis assuming a constant b-value and estimating the corner magnitude of a taper truncation.&lt;/p&gt;