Search

We propose a statistical analysis method to identify common features of seismic activity that are indistinguishable from most other seismicity, and to find anomalous activity that differs from these common features. Using the hypocenter catalog of earthquakes which occurs in and around Japan during the past 20 years, we apply this method for the parameters of the frequency-magnitude distribution and a parameter that expresses the correlation of seismicity with tides as indices, with a focus on objectively understanding the regional characteristics of seismicity. As a result, we extracted a "typical" probability density distribution of each index value common to the most analysis regions and "anomalous" regions with index-value distributions that differ significantly from the typical distributions. In terms of the frequency-magnitude distribution, most estimated values of indices in the anomalous activity areas can be explained as variations corresponding to the effects of fluids, interplate coupling, and stress fields that control faulting styles that have been pointed out in previous studies. By extracting typical index values for the frequency-magnitude distribution, common features of the frequency-magnitude distribution that depend on the earthquake occurrence interval were identified. Although seismicity showed no clear correlation with tides, the index value for tidal correlation changes to reflect the proportion of earthquakes occurring in a series of periods shorter than the tidal period; it is therefore useful as an index to capture the characteristics of such earthquake occurrence intervals. We also show that the typical probability density distribution of these index values can be represented by existing models or their extensions. By using the proposed models as a reference, it is possible to quantify the degree of anomaly using the same concept as that of the method presented here; hence, such a method should be applicable to monitor seismic activity. [ABSTRACT FROM AUTHOR]

Investigation of the characteristic behavior of successive earthquakes that closely occur in space and time is important to understand the generation mechanism of earthquakes and useful to assess a triggered earthquake, especially around the area, where a first large earthquake took place. Here, we analyzed the Global Centroid Moment Tensor catalog from 1976 to 2016 for shallow earthquakes with a moment magnitude, M w , of at least 5.5, and the F-net catalog, Japan, for 4 ≤ M w < 5.5 , to clarify the spatio-temporal characteristics of the successive earthquakes. We first sorted all of the earthquakes in time and removed the aftershocks that occurred in and around the faults of earthquakes with M w larger than the target magnitude range we investigated. Then, we selected source events from the beginning and searched for earthquakes that occurred within a horizontal distance (D) and a lapsed time ( T a ) from the source event to group them in clusters. Then, the source event was selected from the catalog in order, and the same procedure was repeated. We counted the number of clusters, each of which consisted of successive earthquakes, for different D and T a. To examine whether successive earthquakes were explained by random occurrences, we compared the results with simulations in which earthquakes occurred randomly in time but at the same locations matching the centroids in the real data. The comparison showed that the number of clusters for the simulation rapidly increased with D and merged with that for real data at a short distance, which is defined here as the triggering distance. We find that triggering distance is proportional to about 1/5 to 1/4 of the seismic moment (M 0) of the source event, and exponentially decreases with increasing T a . Relating the derived empirical scaling relations between M 0 and triggering distance from the equations in the ETAS model, we show that the observed exponents of 1/5 to 1/4 were well predicted from the estimated ETAS parameters in various regions around the world. These consistencies first show that successive occurrence of earthquakes is well explained by the ETAS model. [ABSTRACT FROM AUTHOR]

Investigation of the characteristic behavior of successive earthquakes that closely occur in space and time is important to understand the generation mechanism of earthquakes and useful to assess a triggered earthquake, especially around the area, where a first large earthquake took place. Here, we analyzed the Global Centroid Moment Tensor catalog from 1976 to 2016 for shallow earthquakes with a moment magnitude, M w , of at least 5.5, and the F-net catalog, Japan, for 4 ≤ M w < 5.5 , to clarify the spatio-temporal characteristics of the successive earthquakes. We first sorted all of the earthquakes in time and removed the aftershocks that occurred in and around the faults of earthquakes with M w larger than the target magnitude range we investigated. Then, we selected source events from the beginning and searched for earthquakes that occurred within a horizontal distance (D) and a lapsed time ( T a ) from the source event to group them in clusters. Then, the source event was selected from the catalog in order, and the same procedure was repeated. We counted the number of clusters, each of which consisted of successive earthquakes, for different D and T a. To examine whether successive earthquakes were explained by random occurrences, we compared the results with simulations in which earthquakes occurred randomly in time but at the same locations matching the centroids in the real data. The comparison showed that the number of clusters for the simulation rapidly increased with D and merged with that for real data at a short distance, which is defined here as the triggering distance. We find that triggering distance is proportional to about 1/5 to 1/4 of the seismic moment (M 0) of the source event, and exponentially decreases with increasing T a . Relating the derived empirical scaling relations between M 0 and triggering distance from the equations in the ETAS model, we show that the observed exponents of 1/5 to 1/4 were well predicted from the estimated ETAS parameters in various regions around the world. These consistencies first show that successive occurrence of earthquakes is well explained by the ETAS model. [ABSTRACT FROM AUTHOR]

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

By interpreting the 2D/3D seismic survey data acquired in the surrounding ocean areas of the Northeast (NE) Japan Arc, we clarified the detailed geological structure and demonstrated that the basic structure in the hanging-wall plate of the subduction system consists of many structural blocks (segments) separated by NW–SE trending large transcurrent faults (strike-slip faults). This structural configuration showed a close relationship with the distribution of foreshocks, mainshock, and aftershocks, coseismic slip models of the 2011 M9.0 Tohoku-Oki megathrust earthquake, coseismic slip area of M-7 class earthquakes, quasi-static slip rates, back slip rate, and seismic tomography images. In addition, the coseismic slip models revealed that the trenchward forearc of the structural blocks between the Offshore Hidaka tectonic line and the Honjo-Sendai tectonic line fitted well with the coseismic slip area of the 2011 Tohoku-Oki earthquake. These findings suggest that the structural blocks bounded by these two tectonic lines slipped rapidly trenchward when the mainshock occurred. The M7 earthquakes were also concentrated along these two tectonic lines, thereby suggesting a close relationship between seismic activity and the inherited geological structure of the overriding plate in the NE Japan forearc. [ABSTRACT FROM AUTHOR]

To explore a local relationship between geological structures and the occurrence of very-low-frequency earthquakes (VLFEs), a particular class of slow earthquakes with characteristic periods of 10–100 s, we investigated three-dimensional (3D) structural features using reprocessed 3D seismic data from the Nankai Trough off Kumano, southwestern Japan. In this region, VLFEs have been observed along the subducting Philippine Sea Plate. Although the detailed source distribution of VLFEs was estimated by means of recent land-based and offshore seismic networks, the relation with geological features is not well understood. First, we reprocessed the 3D seismic data with advanced techniques and reinterpreted the fault distribution in the sediment layer of the accretionary prism and tracked two key horizons: a décollement and the oceanic crust surface. In the accretionary prism sediments, multiple continuous reflectors of basal detachments in the underthrust sequence and conjugate faults cutting the shallow imbricated thrust sequence were identified. In contrast to the gentle variation in the décollement surface, the topographic relief of the oceanic crust was prominent, with ridges and surface displacement due to faults in the oceanic crust. Then, we compared the structural features with the VLFE source locations. Most VLFEs were located deep in the underthrust sediments where the sediments may consist of underconsolidated muds. Furthermore, a high spatial correlation was observed between the VLFE distribution and the oceanic crust topographic relief. The maximum stress direction, which was inferred from the conjugate faults in the imbricated thrust zone, was consistent with the spatial relation between the VLFE localization and the oceanic crust central ridge. Oceanic crust ridges may cause strain accumulation in the underthrust sediments on the landward sides of the ridges, and low-angle slow thrust movements might be caused using weak slip planes in the underthrust muddy sediments. That is, the topographic relief of the oceanic crust may control the occurrence of shallow VLFEs in the Nankai Trough. [ABSTRACT FROM AUTHOR]

Stress drops of small earthquakes have been estimated under the assumption that the rupture propagates symmetrically within a circular fault. However, recent studies have observed directivity effects on seismic waves even for small earthquakes. In this study, rupture directivity was investigated systematically for small-to-moderate-sized earthquakes (M 3.5–5.5) that occurred beneath inland Japan from 2004 to April 2019. Apparent moment rate functions were determined for 1463 earthquakes, and their corner frequencies were inverted for their rupture parameters. The results indicate that most of the analyzed earthquakes (1217 of 1463) are characterized by significantly asymmetrical rupture propagation. The stress drops of the earthquakes estimated by considering asymmetrical rupture propagation were 16.8 MPa, which are almost twice the estimates based on the commonly used symmetrical rupture model. This shows the importance of recognizing the diversity of ruptures, even for small earthquakes, for extracting information about earthquake sources and the Earth's structure. The prevailing rupture directivity can provide useful information on source parameters including the fault size, fault geometry, and rupture velocity of small earthquakes. [ABSTRACT FROM AUTHOR]

In order to clarify the origin of aftershocks, we precisely analyze the hypocenters and focal mechanisms of the aftershocks following the 2000 Western Tottori Earthquake, which occurred in the western part of Japan, using data from dense seismic observations. We investigate whether aftershocks occur on the mainshock fault plane on which coseismic slip occurred or they represent the rupture of fractures surrounding the mainshock fault plane. Based on the hypocenter distribution of the aftershocks, the subsurface fault structure of the mainshock is estimated using principal component analysis. As a result, we can obtain the detail fault structure composed of 8 best-fit planes. We demonstrate that the aftershocks around the mainshock fault are distributed within zones of 1.0-1.5 km in thicknesses, and their focal mechanisms are significantly diverse. This result suggests that most of the aftershocks represent the rupture of fractures surrounding the mainshock fault rather than the rerupture of the mainshock fault. The aftershocks have a much wider zone compared with the exhumed fault zone in field observations, suggesting that many aftershocks occur outside the fault damage zone. We find that most aftershocks except in and around the large-slip region are well explained by coseismic stress changes. These results suggest that the thickness of the aftershock distribution may be controlled by the stress changes caused by the heterogeneous slip distribution during the mainshock. The aftershock is also distributed within a much wider zone than the hypocenter distribution observed in swarm activity in the geothermal region, which is thought to be caused by the migration of hydrothermal fluid. This result implies a difference in generation processes: Stress changes due to the mainshock contribute primarily to the occurrence of aftershocks, whereas earthquake swarms in the geothermal region are caused by fluid migration within the localized zone. [ABSTRACT FROM AUTHOR]

To evaluate local site effects in Kumamoto City, we installed six temporary seismic stations along a 6-km north-south survey line in the city immediately after the 2016 Kumamoto earthquake foreshock (Mj 6.4), which occurred on April 14, 2016. Seismic data from the 2016 Kumamoto earthquake (Mj 7.3), which occurred on April 16, 2016, were successfully recorded at two sites and indicated large amplitudes in the frequency range of 0.5-3 Hz. Site amplifications estimated from weak ground motion data, with a station at Mt. Kinbo used as a reference, are relatively variable along this survey line; however, site amplification factors in the frequency range of 0.5-3 Hz are not large enough to explain the amplitudes produced by the main shock. Nevertheless, site amplifications estimated from strong ground motion data recorded at the two sites during the main shock are large in the frequency range of 1-3 Hz. These findings reveal that the strong ground motions in the frequency range of 1-3 Hz were enhanced by nonlinear behavior of the subsurface soil in Kumamoto City. Moreover, it is observed that the frequency contents of the main shock data in the frequency range of 0.7-3 Hz differ significantly between the two sites, despite the proximity of these sites (600-m interval). Therefore, we also performed single-station microtremor measurements with an interval distance of approximately 100 m between these two sites. We confirmed that the peak frequencies of the horizontal-to-vertical spectral ratios of microtremors have trends that are similar to those of the site amplification factors between the two sites. However, these results could not explain the differences in strong ground motions observed at the two sites during the 2016 Kumamoto earthquake. [ABSTRACT FROM AUTHOR]

Shallow S-wave velocity V profiles were estimated for 26 temporary strong motion observation sites surrounding the epicenters of a sequence of the 2016 Kumamoto earthquake. The microtremor array method was used to gather the dispersion characteristics of Rayleigh waves. V profiles were obtained by inverting the dispersion curves for each site and those of three permanent strong motion stations that recorded the sequence of seismic events. The shallow V profiles near two of the permanent strong motion stations in the town of Mashiki were almost identical. However, the V profiles at other stations varied. The V profiles were found to have the common feature of the uppermost low-velocity layer being widely distributed from Mashiki to the village of Minami-Aso, and it was especially thick in the areas that suffered heavy damage. This low-velocity layer was a major contributor to the site amplification. The horizontal-to-vertical spectral ratios of the microtremors indicate that both the shallow soil and deep sedimentary layers may control the site response characteristics over a broad frequency range. [ABSTRACT FROM AUTHOR]

At the 2016 Kumamoto earthquake, surface ruptures were observed not only along the Futagawa fault, where main ruptures occurred, but also along the Hinagu fault. To estimate the slip distribution on these faults, we extend a method of nonlinear inversion analysis (Fukahata and Wright in Geophys J Int 173:353-364, 2008) to a two-fault system. With the method of Fukahata and Wright (2008), we can simultaneously determine the optimal dip angle of a fault and the slip distribution on it, based on Akaike's Bayesian information criterion by regarding the dip angle as an hyperparameter. By inverting the InSAR data with the developed method, we obtain the dip angles of the Futagawa and Hinagu faults as 61° ± 6° and 74° ± 12°, respectively. The slip on the Futagawa fault is mainly strike slip. The largest slip on it is over 5 m around the center of the model fault (130.9° in longitude) with a significant normal slip component. The slip on the Futagawa fault quickly decreases to zero beyond the intersection with the Hinagu fault. On the other hand, the slip has a local peak just inside Aso caldera, which would be a cause of severe damage in this area. A relatively larger reverse fault slip component on a deeper part around the intersection with Aso caldera suggests that something complicated happened there. The slip on the Hinagu fault is almost a pure strike slip with a peak of about 2.4 m. The developed method is useful in clarifying the slip distribution, when a complicated rupture like the Kumamoto earthquake happens in a remote area. [ABSTRACT FROM AUTHOR]

We report precursory seismic patterns prior to the 2016 Kumamoto earthquakes, as measured by four different methods based on changes in seismicity that can be used for earthquake forecasting: the b-value method, two methods of seismic quiescence evaluation, and an analysis of seismicity density in space and time. The spatial extent of precursory patterns differs from one method to the other and ranges from local scales (typically, asperity size) to regional scales (e.g., 2° × 3° around the source zone). The earthquakes were preceded by periods of pronounced anomalies, which lasted in yearly scales (1.5 years), or longer (>3 years). We demonstrate that a combination of multiple methods detected different signals prior to the Kumamoto earthquakes. This indicates great potential to reduce the hazard at possible future sites of earthquakes relative to long-term seismic hazard assessment. We also found that the seismic quiescence in a regional-scale area, detected by using the two methods of seismic quiescence evaluation, was a common precursor to the 2016 Kumamoto earthquakes and 2015 Off Satsuma Peninsula earthquake. The result allows us to interpret both events as the onset that occurred at a section along the tectonic line from the Okinawa Trough through the Beppu-Shimabara graben. [ABSTRACT FROM AUTHOR]

An ENE-trending ~30-km-long surface rupture emerged during the Mw = 7.0 16 April 2016 Kumamoto earthquake along the previously mapped Futagawa and northern Hinagu faults. This included a previously unknown 5-km-long fault within the Aso Caldera, central Kyushu. The rupture zone is mostly composed of right-lateral slip sections, with a maximum of 2-m coseismic slip. One of the noteworthy features we observed in the field are ~10-km-long segmented normal fault scarps, dipping to the northwest, along the previously mapped Idenokuchi fault, 1.2-2.0 km south of and subparallel to the Futagawa fault. The maximum amount of coseismic throw on the Idenokuchi fault is ~2 m, which is nearly equivalent to the maximum slip on the strike-slip rupture. The locations and slip motions of the 2016 rupture are also manifested as interferogram fringe offsets in InSAR images. Together with geodetic and seismic inversions of subsurface fault slip, we present a schematic structural model where oblique motion occurred on a northwest-dipping subsurface fault and the slip is partitioned at the surface into strike-slip and normal fault scarps. Our simple dislocation model demonstrates that this bifurcation into pure strike-slip and normal faults likely occurs for optimally oriented failure near the surface. The Kumamoto case, with detailed geological observations and geophysical models, would be the second significant slip-partitioned earthquake around the globe. It provides an important insight into scale- and depth-dependent stress heterogeneity and an implication to a proper estimate of seismic hazard in complex and broad multiple fault strands. [ABSTRACT FROM AUTHOR]

We present the initial results of rapid fault estimations for the 2016 Kumamoto earthquake on April 16 ( M 7.3), and coseismic displacements caused by the two large foreshocks that occurred on April 14 ( M 6.5) and April 15 ( M 6.4) from the GEONET real-time analysis system (REGARD), which is based on a Global Navigation Satellite System (GNSS) kinematic positioning technique. The real-time finite-fault estimate ( M 6.85) was obtained within 1 min and converged to M 6.96 within 5 min of the origin time of the mainshock ( M 7.3). The finite-fault estimate shows right-lateral strike-slip fault along the Futagawa fault segment, which is consistent with the finite-fault model inferred from post-processed GNSS and InSAR analysis. Furthermore, significant coseismic displacements were observed due to the April 14 and April 15 foreshocks at nearby sites, though these earthquakes were smaller than the pre-assigned system threshold. Our results also demonstrate the potential for the GNSS-based earthquake early warning system for inland earthquakes. [ABSTRACT FROM AUTHOR]

We estimate high-frequency (4-20 Hz) energy release due to the 2016 Kumamoto, Japan, earthquake sequence, within a time period from April 14 to 26 through envelope inversion analysis applied to the Hi-net continuous seismograms. We especially focus on energy releases after each of the April 14 M 6.5 and the April 16 M 7.3 earthquakes. The cumulative energy release from aftershocks of the April 14 event reaches 60% of that from the April 14 event itself by the lapse time of 27 h (pre-April 16 period). On the other hand, the cumulative energy release from aftershocks of the April 16 event reaches only 11 and 13% of that from the April 16 event itself by the lapse times of 27 h and 10 days (post-April 16 period), respectively. This discrepancy in the normalized cumulative energy release (NCER) indicates that the April 14 event was followed by much larger relative aftershock productivity than the April 16 event. Thus, NCER would provide information that reflects relative aftershock productivity and ongoing seismicity pattern after a large earthquake. We also find that the temporal decay of the energy release rate obeys the power law. The exponent p of the power-law decay is estimated to be 1.7-2.1, which is much larger than the typical p value of the Omori-Utsu law: slightly larger than 1. We propose a simple relationship given by p = βp/ b, where p value, b value of the Gutenberg-Richter law, and β value of the magnitude-energy release relationship are combined. [ABSTRACT FROM AUTHOR]

An M6.5 earthquake (foreshock) and M7.3 earthquake (mainshock) struck Kumamoto Prefecture on April 14, 2016, and April 16, 2016. To evaluate the effect of crustal deformation due to the earthquake on the Aso magma system, we detected crustal deformation using InSAR and GNSS. From InSAR analysis, we detected large crustal deformations along the Hinagu Fault, the Futagawa Fault, and the northeast extension of the latter fault. It extended to more than 50 km, and the maximum slant-range change exceeded 1 m. Although the obtained crustal deformation was approximately explained by the right-lateral strike-slip on the fault, its details could not be explained by such simple faulting. Additionally, we found complex surface deformation west of the Aso caldera rim, suggesting that shallow fault slips occurred in many known and unknown faults associated with the earthquake. Most of the crustal deformation could be reasonably explained by four rectangle faults located along the Futagawa Fault, in the northeast extension of the Futagawa Fault, alongside the Hinagu Fault, and in the eastern part of the Futagawa Fault. The first three of faults have high dip angles and right-lateral slip. The other was a fault with a low dip angle that branched from the shallow depth of the fault along the Futagawa Fault. The normal-dip right-lateral slip was estimated for this segment. Based on the estimated fault model, we calculated the displacement and stress field around the Aso volcano by the finite-element method (FEM) to evaluate the effects on the Aso magma system. In this calculation, we assumed a spherical soft medium located at a 6-km depth beneath the area south of the Kusasenri region as the magma system and considered only static effects. The result shows complex distributions of displacements and stresses, but we can notice the following significant points. (1) The spherical magma system deformed to an ellipsoid, and the total volume was slightly increased, less than 1%. (2) The differential stress around the upper portion of the magma system was as large as 3.5 MPa. This is strong enough to open pre-existing cracks and can cause the migration of magma. [ABSTRACT FROM AUTHOR]

The 2016 Kumamoto earthquake sequence was a rare event worldwide in that the surface ruptures associated with the largest foreshock (Mj 6.5) of 21:26 (JST), 14 April ruptured again during the mainshock (Mj 7.3) of 01:25 (JST), 16 April. The 14 April Mj 6.5 earthquake produced 6-km-long surface ruptures along the central portion of the Futagawa-Hinagu fault zone (FHFZ). The mainshock produced 31-km-long surface ruptures along the central to northeastern part of the FHFZ. Field observations and eyewitness accounts documented that the offsets of the ruptures associated with the 14 April foreshock became larger after the 16 April mainshock, suggesting that the same portion of the fault ruptured to the surface twice in the Kumamoto earthquake sequence. The 6-km-long surface ruptures associated with the largest foreshock are located near a geometric bend of the FHFZ characterized by ~50° change in strike. The epicenter of the mainshock is also located near the bend. These observations imply that the Kumamoto earthquake sequence was initiated due to a stress concentration on the bend of the FHFZ, and the mainshock was initiated approximately at the same place about 28 h after the largest foreshock. This foreshock/mainshock sequence of the Kumamoto earthquake is not successive events on the adjacent different fault zones, because the 6-km-long surface ruptures of the largest foreshock are part of the 31-km-long surface ruptures of the mainshock. [ABSTRACT FROM AUTHOR]

A driving stress of the Mw5.8 reverse-faulting Awaji Island earthquake (2013), southwest Japan, was investigated using focal mechanism solutions of earthquakes before and after the mainshock. The seismic records from regional high-sensitivity seismic stations were used. Further, the stress tensor inversion method was applied to infer the stress fields in the source region. The results of the stress tensor inversion and the slip tendency analysis revealed that the stress field within the source region deviates from the surrounding area, in which the stress field locally contains a reverse-faulting component with ENE–WSW compression. This local fluctuation in the stress field is key to producing reverse-faulting earthquakes. The existing knowledge on regional-scale stress (tens to hundreds of km) cannot predict the occurrence of the Awaji Island earthquake, emphasizing the importance of estimating local-scale (< tens of km) stress information. It is possible that the local-scale stress heterogeneity has been formed by local tectonic movement, i.e., the formation of flexures in combination with recurring deep aseismic slips. The coseismic Coulomb stress change, induced by the disastrous 1995 Mw6.9 Kobe earthquake, increased along the fault plane of the Awaji Island earthquake; however, the postseismic stress change was negative. We concluded that the gradual stress build-up, due to the interseismic plate locking along the Nankai trough, overcame the postseismic stress reduction in a few years, pushing the Awaji Island earthquake fault over its failure threshold in 2013. The observation that the earthquake occurred in response to the interseismic plate locking has an important implication in terms of seismotectonics in southwest Japan, facilitating further research on the causal relationship between the inland earthquake activity and the Nankai trough earthquake. Furthermore, this study highlighted that the dataset before the mainshock may not have sufficient information to reflect the stress field in the source region due to the lack of earthquakes in that region. This is because the earthquake fault is generally locked prior to the mainshock. Further research is needed for estimating the stress field in the vicinity of an earthquake fault via seismicity before the mainshock alone. [ABSTRACT FROM AUTHOR]

We investigate S-wave attenuation (Qs) structure in and around Hokkaido, Japan, including the source area of the 2018 Hokkaido Eastern Iburi Earthquake (M 6.7) and its aftershocks. From the strong-motion seismograms recorded in the nationwide seismograph network in Japan, we compute spectrum amplitude at 1–10 Hz with every 1 Hz and then estimate Qs structure by adopting tomographic inversion technique. The obtained Qs structure suggests that lateral variations of Qs are remarkably formed in the study areas. The low-Qs anomalies distributing north to south are imaged in both sides of the Hidaka mountain range which is built in the central Hokkaido. These characteristics structures of Qs are extended from the ground surface to depths by about 50 km, indicating that heterogeneous structures resulted from the collision of the Kuril and northeastern Japan arcs develop at the depth ranges. Around the source area of the 2018 Hokkaido Eastern Iburi Earthquake, the abrupt changes of Qs values are identified. The 2018 Hokkaido Eastern Iburi Earthquake and its aftershocks seem to be arraigned on the high-Qs zones, whereas the low-Qs anomalies are also imaged next to the high-Qs zones. The boundary of the low- and high-Qs zones lies with near vertical, corresponding to aftershock distributions of the 2018 Hokkaido Eastern Iburi Earthquake. This agreement would propose that the heterogeneities formed due to the arc–arc collision would characterize faulting processes in the main shock of the Hokkaido Eastern Iburi Earthquake and the activities of its aftershocks. [ABSTRACT FROM AUTHOR]

In this study, we investigated the influence of local structural heterogeneities on aftershocks of the 2008 Iwate–Miyagi Nairiku earthquake (M 7.2) that occurred in the northeastern Japan arc. Although this area is characterized by an E–W compressional thrust faulting stress regime, many N–S compressional thrust-type and strike-slip-type earthquakes were also observed in the aftershock activity of the 2008 event. The occurrence of such N–S compressional aftershocks indicates that the differential stress was very low before the main shock. We investigated Vp/Vs structure in the aftershock area in detail to find most of the area showed Vp/Vs as low as less than 1.70, which is consistent with the previous studies. Our result shows, however, the aftershock area was dotted with high Vp/Vs small areas, which suggests the locations of abundant fluid-filled cracks. The comparison of the distributions of Vp/Vs, the stress changes due to the mainshock, and focal mechanisms of the aftershocks indicate that the N–S compressional aftershocks are concentrated in regions with high Vp/Vs and large N–S compressional stress change. From these results, we concluded that the N–S compressional aftershocks occurred in regions with low strength due to the high pore pressure and with N–S compressional stress caused by the stress change due to the mainshock in a low differential stress regime. [ABSTRACT FROM AUTHOR]

We reconstructed the spatiotemporal evolution of seismicity associated with the 2018 Mw 5.6 northern Osaka earthquake, Japan, to discuss the source fault model of the mainshock rupture, the possible link between this rupture and known active faults, and subsequent crustal deformation. We first relocated the hypocenters listed in the earthquake catalog determined by the Japan Meteorological Agency using a double-difference relocation algorithm. We then searched for the earthquake waveforms that closely resembled those of the relocated hypocenters by applying a matched filtering technique to the continuous waveform data. The relocated hypocenters revealed two distinct planar alignments with different fault geometries. A combination of the relocated hypocenters and focal mechanisms suggests that the mainshock rupture initiated on a NNW-SSE-striking thrust fault, dipping ~ 45° to the east, with the rupture propagating to an adjacent sub-vertical ENE-WSW-striking strike-slip fault ~ 0.3 s after the initial mainshock rupture, resulting in the simultaneous propagation of dynamic rupture along the two faults. The strike-slip fault is oblique to the strike of the Arima-Takatsuki Fault, indicating that blind strike-slip faulting occurred. While the east-dipping thrust fault is located deeper than the modeled extent of the Uemachi Fault, a simple extrapolation of the near-surface geometry of the Uemachi fault partially overlaps the mainshock rupture area. Although it is unclear as to whether a blind thrust fault or a deep portion of the Uemachi Fault ruptured during this mainshock-aftershock sequence, a mainshock rupture would have transferred a static stress change of > 0.1 MPa to a portion of the east-dipping thrust fault system. Intensive aftershocks have persisted along the northern and southern edges of the source area, including moderate-magnitude events, whereas the seismicity in the central part of the source area has shown a rapid decay over time. Delayed triggered aftershocks were clearly identified along the northern extension of the rupture area. Because the background seismicity is predominant in this northern area, we interpret that aseismic deformation, such as cataclastic flow lubricated by crustal fluids, triggered this off-fault seismicity. [ABSTRACT FROM AUTHOR]