Large subduction-zone earthquakes generate long-lasting and wide-spread aftershock sequences. The physical and statistical patterns of these aftershock sequences are of considerable importance for better understanding earthquake dynamics and for seismic hazard assessments and earthquake risk mitigation. In this work, we analyzed the statistical properties of 42 aftershock sequences in terms of their temporal evolution. These aftershock sequences followed recent large subduction-zone earthquakes of M ≥ 7.0 with focal depths less than 70 km that have occurred worldwide since 1976. Their temporal properties were analyzed by investigating the probability distribution of the interevent times between successive aftershocks in terms of non-extensive statistical physics (NESP). We demonstrate the presence of a crossover behavior from power-law (q ≠ 1) to exponential (q = 1) scaling for greater interevent times. The estimated entropic q-values characterizing the observed distributions range from 1.67 to 1.83. The q-exponential behavior, along with the crossover behavior observed for greater interevent times, are further discussed in terms of superstatistics and in view of a stochastic mechanism with memory effects, which could generate the observed scaling patterns of the interevent time evolution in earthquake aftershock sequences.
Keywords: Tsallis entropy; aftershock sequences; interevent times; megathrust earthquakes; subduction zones; superstatistics.