Changes in the stress state of faults and their surroundings is a highly plausible mechanism explaining earthquake interaction. These stress changes can impact the seismicity rate and the size distribution of earthquakes. However, the effect of large earthquakes on the earthquake size distribution along the Longmenshan fault has not been quantified. We evaluated the levels of the b value for the stable state before and after the large earthquakes on 12 May 2008 (Wenchuan, M_S 8.0) and 20 April 2013 (Lushan, M_S 7.0) along the Longmenshan fault. We found that after the mainshocks, the size distribution of the subsequent earthquakes shifted toward relatively larger events in the Wenchuan aftershock zone (b value decreased from 1.21 to 0.84), and generally remained invariable in the Lushan aftershock zone (b value remained at 0.76). The time required for the b value to return to stable states after both mainshocks was entirely consistent with the time needed by the aftershock depth images to stop visibly changing. The result of the temporal variation of b values shows decreasing trends for the b value before both large earthquakes. Our results are available for assessing the potential seismic risk of the Longmenshan fault as a reference.

Seismically active fault zones are complex natural systems and they exhibit multifractal correlation between earthquakes in space and time. In this paper, the seismicity of the Taiwan seismic region was studied through the multifractal characteristics of the spatial-temporal distribution of earthquakes from 1st January 1995 to 1st January 2019. We quantified the multifractal characteristics of Taiwan at different scales and defined them as ΔD values. Furthermore, we studied the relationship between the ΔD and b values, which signifies the average size distribution of those earthquakes. The results are as follows. (1) The temporal multifractal curve changes substantially before and after the strong earthquakes. (2) The maximum ΔD value of the seismic region in Taiwan occurs at depths of 0~9 km, indicating that geological structures and focal mechanisms is the most complex at these depths compared with other depths. (3) ΔD values for different regions range from 0.2~1.5, and b values range from 0.65~1.3, with a significant positive correlation between them (ΔD = 1.5 × b − 0.68). For this purpose, a statistical relationship is developed between b and ΔD values, and regional and temporal changes of these parameters are analyzed in order to reveal the potential of future earthquakes in the study region.