In this study, we investigated temporal variations in stress drop and b‐value in the earthquake swarm that occurred at the Yamagata‐Fukushima border, NE Japan, after the 2011 Tohoku‐Oki earthquake. In this swarm, frictional strengths were estimated to have changed with time due to fluid diffusion. We first estimated the source spectra for 1,800 earthquakes with 2.0 ≤ MJMA < 3.0, by correcting the site‐amplification and attenuation effects determined using both S waves and coda waves. We then determined corner frequency assuming the omega‐square model and estimated stress drop for 1,693 earthquakes. We found that the estimated stress drops tended to have values of 1–4 MPa and that stress drops significantly changed with time. In particular, the estimated stress drops were very small at the beginning, and increased with time for ~50 days. Similar temporal changes were obtained for b‐value; the b‐value was very high (b ~ 2) at the beginning, and decreased with time, becoming approximately constant (b ~ 1) after ~50 days. Patterns of temporal changes in stress drop and b‐value were similar to the patterns for frictional strength and earthquake occurrence rate, suggesting that the change in frictional strength due to migrating fluid not only triggered the swarm activity but also affected earthquake and seismicity characteristics. The estimated high Q−1 value, as well as the hypocenter migration, supports the presence of fluid, and its role in the generation and physical characteristics of the swarm.

A shallow M7.3 event with a M6.5 foreshock occurred along the Futagawa‐Hinagu fault zone in Kyushu, SW Japan. We investigated the spatiotemporal variation of the stress orientations in and around the source area of this 2016 Kumamoto earthquake sequence by inverting 1218 focal mechanisms. The results show that the σ3 axis in the vicinity of the fault plane significantly rotated counterclockwise after the M6.5 foreshock and rotated clockwise after the M7.3 main shock in the Hinagu fault segment. This observation indicates that a significant portion of the shear stress was released both by the M6.5 foreshock and M7.3 main shock. It is estimated that the stress release by the M6.5 foreshock occurred in the shallower part of the Hinagu fault segment, which brought the stress concentration in its deeper part. This might have caused the M7.3 main shock rupture mainly along the deeper part of the Hinagu fault segment after 28 h.

To understand role of fluid on earthquake rupture processes, we investigated effects of thermal pressurization on spatial variation of dynamic rupture by computing spontaneous rupture propagation on a rectangular fault. We found thermal pressurization can cause heterogeneity of rupture even on a fault of uniform properties. On drained faults, tractions drop linearly with increasing slip in the same way everywhere. However, by changing the drained condition to an undrained one, the slip‐weakening curves become non‐linear and depend on locations on faults with small shear zone thickness w, and the dynamic frictional stresses vary spatially and temporally. Consequently, the super‐shear transition fault length decreases for small w, and the final slip distribution can have some peaks regardless of w, especially on undrained faults. These effects should be taken into account of determining dynamic rupture parameters and modeling earthquake cycles when the presence of fluid is suggested in the source regions.

This paper reports stick-slip behaviors of Indian gabbro as studied using a new large-scale biaxial friction apparatus, built in the National Research Institute for Earth Science and Disaster Prevention (NIED), Tsukuba, Japan. The apparatus consists of the existing shaking table as the shear-loading device up to 3,600 kN, the main frame for holding two large rectangular prismatic specimens with a sliding area of 0.75 m 2 and for applying normal stresses r n up to 1.33 MPa, and a reaction force unit holding the stationary specimen to the ground. The shaking table can produce loading rates v up to 1.0 m/s, accelerations up to 9.4 m/s 2 , and displacements d up to 0.44 m, using four servocontrolled actuators. We report results from eight preliminary experiments conducted with room humidity on the same gabbro specimens at v = 0.1-100 mm/s and r n = 0.66-1.33 MPa, and with d of about 0.39 m. The peak and steady-state friction coefficients were about 0.8 and 0.6, respectively, consistent with the Byerlee friction. The axial force drop or shear stress drop during an abrupt slip is linearly proportional to the amount of displacement, and the slope of this relationship determines the stiffness of the apparatus as 1.15 9 10 8 N/m or 153 MPa/m for the specimens we used. This low stiffness makes fault motion very unstable and the overshooting of shear stress to a negative value was recognized in some violent stick-slip events. An abrupt slip occurred in a constant rise time of 16-18 ms despite wide variation of the stress drop, and an average velocity during an abrupt slip is linearly proportional to the stress drop. The use of a large-scale shaking table has a great potential in increasing the slip rate and total displacement in biaxial friction experiments with large specimens.

No abstract