We analyse seismicity during the 6-year period 2012–2017 in the new time domain termed natural time in the Chiapas region where the M8.2 earthquake occurred, Mexico’s largest earthquake in more than a century, in order to study the complexity measures associated with fluctuations of entropy as well as with entropy change under time reversal. We find that almost three months before the M8.2 earthquake, i.e., on 14 June 2017, the complexity measure associated with the fluctuations of entropy change under time reversal shows an abrupt increase, which, however, does not hold for the complexity measure associated with the fluctuations of entropy in forward time. On the same date, the entropy change under time reversal has been previously found to exhibit a minimum [Physica A 506, 625–634 (2018)]; we thus find here that this minimum is also accompanied by increased fluctuations of the entropy change under time reversal. In addition, we find a simultaneous increase of the Tsallis entropic index q.

This article is focused on a new procedure concerning a more accurate identification of the occurrence time of an impending major earthquake (EQ). Specifically, we first recapitulate that, as was recently shown [P. Varotsos et al., Communications in Nonlinear Science and Numerical Simulation 125 (2023) 107370], natural time analysis of seismicity supplemented with the non-additive Tsallis entropy $S_q$ leads to a shortening of the time window of an impending major EQ. This has been shown for the Tohoku mega-EQ of magnitude M9 that occurred in Japan on 11 March 2011, which is the largest event ever recorded in Japan. Here, we also show that such a shortening of the time window of an impending mainshock can be achieved for major, but smaller EQs, of the order of M8 and M7. In particular, the following EQs are treated: the Chiapas M8.2 EQ, which is Mexico’s largest EQ for more than a century that took place on 7 September 2017 near the coast of Chiapas state in Mexico, the 19 September 2017 M7.1 EQ that occurred within the Mexican flat slab, and the M7.1 Ridgecrest EQ on 6 July 2019 in California.

The changes in the seismic signals generated by avalanches recorded at three sites along a path at the Vallée de la Sionne (VdlS) experimental site are presented. We discuss and correlate the differences in the duration, signal amplitudes, and frequency content of the sections (Signal ONset (ON), Signal Body (SBO), and Signal TAil and Signal ENd STA-SEN) of the spectrograms with the evolution of the powder, transitional and wet snow avalanches along a path. The development of the avalanche front was quantified using the exponential function in time F (t) = K’ exp (β t) fitted to the shape of the signal ONset (SON section of the spectrogram. The speed of the avalanche front is contained in β. To this end, a new method was developed. The three seismic components were converted into one seismic component (FS), when expressing the vector in polar coordinates. We linked the theoretical function of the shape of the FS-SON section of the spectrogram to the numerical coefficients of its shape after considering the spectrogram as an image. This allowed us to obtain the coefficients K’ and β. For this purpose, the Hough Transform (HT) was applied to the image. The values of the resulting coefficients K’ and β are included in different ranges in accordance with the three types of avalanche. Curves created with these coefficients enable us to estimate the development of the different avalanche types along the path. Our results show the feasibility of classifying the type of avalanche through these coefficients. Average speeds of the avalanches approaching the recording sites were estimated. The speed values of wet and transitional avalanches are consistent with those derived from GEODAR (GEOphysical Doppler radAR) measurements, when available. The absence of agreement in the speed values obtained from seismic signals and GEODAR measurements for powder snow avalanches indicates, for this type of avalanche, a different source of the measured signal. Hence, the use of the two measuring systems proves to be complementary.

One of the most important subduction zones in the world is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate. One part of it is located in the Mexican Pacific Coast, where the Cocos plate inserts beneath the North American plate with different dip angles, showing important seismicity. Under the central Mexican area, such a dip angle becomes practically horizontal and such an area is known as flat slab. An earthquake of magnitude M7.1 occurred on 19 September 2017, the epicenter of which was located in this flat slab. It caused important human and material losses of urban communities including a large area of Mexico City. The seismicity recorded in the flat slab region is analyzed here in natural time from 1995 until the occurrence of this M7.1 earthquake in 2017 by studying the entropy change under time reversal and the variability β of the order parameter of seismicity as well as characterize the risk of an impending earthquake by applying the nowcasting method. The entropy change ΔS under time reversal minimizes on 21 June 2017 that is almost one week after the observation of such a minimum in the Chiapas region where a magnitude M8.2 earthquake took place on 7 September 2017 being Mexico’s largest quake in more than a century. A minimum of β was also observed during the period February–March 2017. Moreover, we show that, after the minimum of ΔS, the order parameter of seismicity starts diminishing, thus approaching gradually the critical value 0.070 around the end of August and the beginning of September 2017, which signals that a strong earthquake is anticipated shortly in the flat slab.

It has been shown that some dynamic features hidden in the time series of complex systems can be unveiled if we analyze them in a time domain termed natural time. In this analysis, we can identify when a system approaches a critical point (dynamic phase transition). Here, based on natural time analysis, which enables the introduction of an order parameter for seismicity, we discuss a procedure through which we could achieve the identification of the occurrence time of the M8.2 earthquake that occurred on 7 September 2017 in Mexico in Chiapas region, which is the largest magnitude event recorded in Mexico in more than a century. In particular, we first investigated the order parameter fluctuations of seismicity in the entire Mexico and found that, during an almost 30-year period, i.e., from 1 January 1988 until the M8.2 earthquake occurrence, they were minimized around 27 July 2017. From this date, we started computing the variance of seismicity in Chiapas region and found that it approached the critical value 0.070 on 6 September 2017, almost one day before this M8.2 earthquake occurrence.