Explore the vast range of titles available.

The 2018, Hokkaido Eastern Iburi, Japan, earthquake is an event characterized by complexity of the rupture process and slip pattern, which may involve both reverse and strike-slip motion depending on the locations on the fault surface. We perform dynamic rupture simulations based on simple physical laws, conditions for stressing and fault friction, and the non-planar fault geometry constrained by the aftershock observation. The complex fault geometry is numerically treated by the boundary integral equation method accelerated by the fast domain portioning method. The fault geometry is characterized primarily by the combination of six fault planes. As a result, we are able to explain several observed features of the event, including the spatial variation of the final fault slip and rupture velocity, which are inferred from the kinematic slip inversion. We also succeed in refining the constraint of the regional stress field in the focal area based on the simulation. Our results show that the overall patterns of the complex rupture event can be reproduced by a relatively simple model of the regional stress and the fault friction, if the geometrical complexity of the fault is properly taken into account.

The Hokkaido Eastern Iburi Earthquake (M = 6.7) occurred on Sep. 6, 2018 in the southern part of Central Hokkaido, Japan. Since Paleogene, this region has experienced the dextral oblique transpression between the Eurasia and North American (Okhotsk) Plates and the subsequent collision between the Northeast Japan Arc and the Kuril Arc due to the oblique subduction of the Pacific Plate. This earthquake occurred beneath the foreland fold-and-thrust belt of the Hidaka Collision zone developed by the collision process, and is characterized by its deep focal depth (~ 37 km) and complicated rupture process. The reanalyses of controlled source seismic data collected in the 1998–2000 Hokkaido Transect Project revealed the detailed structure beneath the fold-and-thrust belt, and its relationship with the aftershock activity of this earthquake. Our reflection processing using the CRS/MDRS stacking method imaged for the first time the lower crust and uppermost mantle structures of the Northeast Japan Arc underthrust beneath a thick (~ 5–10 km) sedimentary package of the fold-and-thrust belt. Based on the analysis of the refraction/wide-angle reflection data, the total thickness of this Northeast Japan Arc crust is only 16–22 km. The Moho is at depths of 26–28 km in the source region of the Hokkaido Eastern Iburi Earthquake. Our hypocenter determination using a 3D structure model shows that most of the aftershocks are distributed in a depth range of 7–45 km with steep geometry facing to the east. The seismic activity is quite low within the thick sediments of the fold–thrust belt, from which we find no indication on the relationship of this event with the shallow (< 10–15 km) and rather flat active faults developed in the fold-and-thrust belt. On the other hand, a number of aftershocks are distributed below the Moho. This high activity may be caused by the cold crust delaminated from the Kuril Arc side by the arc–arc collision, which prevents the thermal circulation and cools the forearc uppermost mantle to generate an environment more favorable for brittle fracture.

The Noto Peninsula, extending northward into the Sea of Japan, features a narrow, elongated shape, complex coastal topography, and numerous active faults along its coastline. Since December 2020, intense earthquake swarms accompanied by crustal deformation have occurred in the northeastern peninsula, likely caused by fluid upwelling from deep underground. The largest event, a Magnitude 7.6 earthquake, struck on January 1, 2024, with aftershock distributions indicating multiple faults ruptured over approximately 150 km. This study aimed to clarify the temporal and spatial variation in seismic wave radiation and investigate the source process of the M7.6 event using the back-projection method. This method estimates the origin of wave packets recorded by a seismic array. In Japan, seismic networks operated by local governments often include densely distributed stations to evaluate seismic intensity. We used these dense sites as a seismic array complemented by strong ground motion data from NIED K-NET and KiK-net. The analysis assumed three fault planes, based on previous studies. Velocity waveforms in two frequency bands (0.05–2.0 Hz and 0.5–5.0 Hz) were used to estimate areas of strong radiation intensity, representing the sources of seismic waves. In the low-frequency band, strong radiation intensity was observed near the rupture initiation point and in shallow regions of the northern Noto Peninsula, corresponding to large fault slips that caused the uplift of the coastline. In contrast, no strong radiation intensity was detected off the northeast coast of the Noto Peninsula in the low-frequency band, suggesting the absence of a significant slip. High-frequency analysis revealed distributions of strong radiation intensities complementary to those in the low-frequency band. A subevent occurring around 20 s after the rupture initiation was found to originate near the northern coast of the Noto Peninsula. Graphical Abstract

We investigate spatiotemporal changes of intermediate‐depth earthquakes in the double seismic zone beneath Central and Northeastern Japan before and after the 2011 magnitude 9 Tohoku earthquake. We build a template‐matching catalog 1 year before and 1 year after the Tohoku earthquake using Hi‐net recordings. The new catalog has a six‐fold increase in earthquakes compared to the Japan Meteorological Agency catalog. Our results show no significant change in the intermediate‐depth earthquake rate prior to the Tohoku earthquake, but a clear increase in both planes following the Tohoku earthquake. The regions with increased intermediate‐depth earthquake activity and the post‐seismic slips following the Tohoku earthquake are spatially separate and complementary with each other. Aftershock productivity of intermediate‐depth earthquakes increased in both planes following the Tohoku earthquake. Overall, aftershock productivity of the upper plane is higher than the lower plane, likely indicating that stress environments and physical mechanisms of intermediate‐depth earthquakes in the two planes are distinct. Plain Language Summary Intermediate‐depth earthquakes occur at 70–350 km depth below the Earth's surface. Because of the high pressure and temperature at such depths, the physical mechanism of intermediate‐depth earthquakes is still poorly understood. In this study, we try to obtain insights into this problem by investigating the behaviors of intermediate‐depth earthquakes before and after the 2011 magnitude 9 Tohoku earthquake. We build an intermediate‐depth earthquake catalog in Central and Northeastern Japan around the Tohoku earthquake. This new catalog contains six times more earthquakes and is more complete than the standard catalog used in previous studies. After analyzing the earthquakes from the new catalog, we find no significant precursory acceleration of intermediate‐depth earthquake activities prior to the Tohoku earthquake. However, we observe a clear increase in intermediate‐depth earthquake rate following the Tohoku earthquake. We also observe that the aftershock productivity of intermediate‐depth earthquakes increased following the Tohoku earthquake. The subduction slab beneath our studied area is characterized by a double seismic zone, and we find that aftershock productivity in the upper plane of seismicity is higher than that in the lower plane. This phenomenon may be due to differences in the environments or physical mechanisms in the two planes of intermediate‐depth earthquakes. Key Points We build a new catalog of intermediate‐depth earthquakes in Japan 1 year before and 1 year after the 2011 M9 Tohoku earthquake There is no significant increase in intermediate‐depth earthquake activities prior to the Tohoku earthquake but a clear increase after it Aftershock productivities of double seismic zones increase after the Tohoku earthquake, and it in upper plane is higher than in lower plane

In the original publication of this article (Iwasaki et al. 2019), the author name ‘Matsubara Makoto’ in author list is not correct. The correct name should be ‘Makoto Matsubara’.

We conducted seismic tomography for entire Japanese Islands including the Sea of Japan and the Pacific Ocean using arrival times from reflection survey as well as the routine seismic network. We successfully imaged the shallow zone along the Sea of Japan from offshore Yamagata to the Noto Peninsula by using air gun data. An extremely low-V shallow zone is imaged between Sado Island and Noto Peninsula. We also obtained detailed seismic velocity structure beneath the Pacific Ocean at depths of 20–50 km using S-net data. The 2007 Noto Peninsula, the 2007 offshore Chuetsu, and the 2019 offshore Yamagata earthquakes occurred at the boundary between high-Vp and low-Vp zones. The west side of the hypocenter of the 2019 offshore Yamagata earthquake at depths of 10–30 km has high-V corresponding to the Mogami Trough. This high-V zone passes through Awa Island and reaches Sado Basin between Sado Island and Honshu. A major rift zone in the Tohoku Arc extending from the Akita region to the Niigata region along the coast of Sea of Japan corresponds to high-V lower crust and a shallow Moho.

We present the results of relocating 17,544 hypocenters determined from data recorded during the 2016 Kumamoto Earthquake sequence, during the interval between April 14, 2016, and August 31, 2016. For this, we used a double-difference relocation method to constrain high-resolution hypocenter locations by cross-correlation differential times as well as the NIED Hi-net catalog differential times. The sequence included two large events (on 14 April: M JMA 6.5 and on 16 April: M JMA 7.3) that occurred in a complicated region where the Hinagu and Futagawa faults meet. By comparing these high-resolution earthquake locations in three different periods [(P1) between 2001 and 2012; (P2) between M JMA 6.5 and M JMA 7.3; and (P3) between M JMA 7.3 and August 31, 2016], we present the significant seismicity after the mainshock relative to the background seismicity. Events during the Kumamoto Earthquake sequence occurred generally within the same sites of known faults and background seismicity. For an example, the seismicity during period P2 formed a sharp linear shape along the northern part of the Hinagu fault for about 20 km. A series of linear seismicity events occurred during period P3 along the Futagawa fault to the east (for about 28 km), in the northern part of the Aso caldera, and in the Oita region around the Beppu–Haneyama fault zone. These events also extended to the mid- and southern parts of the Hinagu fault zone and were shaped only after the M 7.3 event. Moreover, high-resolution hypocenter locations also allowed us to identify some clusters of events that occurred in regions where background seismicity has not been confirmed. For instance, activity on the northwestern edge of the Aso caldera and in small areas within the Beppu–Haneyama fault zone became apparent with new seismic activity. We also demonstrate herein the absence of seismicity between the northeast extension of the Futagawa fault zone and the Aso caldera region, which became clearly shown after the M 7.3 event. This low-seismicity region is located at the boundary of the low- and high-velocity structures and different focal mechanisms, but is also close to the maximum slip area of the M 7.3 event.

We conduct seismic tomography using arrival time data picked by NIED Hi-net, including earthquakes off the coast, outside the seismic network. For these offshore events, we use the NIED F-net focal depth. We detect two low- V . zones in the uppermost subducting oceanic crust. The landward low- V zone with a large anomaly corresponds to the western edge of the coseismic slip zone of the 2011 off the Pacific coast of Tohoku Earthquake. The asperities of the previously known Off-Miyagi and Off-Fukushima earthquakes with magnitudes around 7.0 are also located at the boundary of the low- V and the eastern high- V zones. The initial break point (hypocenter) is associated with the edge of a slightly low- V and low- V p / V s zone. The trenchward low- V and low- V p / V s zone extending southwestward from the hypocenter may indicate the existence of a subducted seamount. The high- V zone and low- V p / V s zone might have accumulated the strain and resulted in the huge coseismic slip zone of the 2011 Tohoku Earthquake. The low- V and low- V p / V s zone is a slight fluctuation within the high- V zone and might have acted as the initial break point of the 2011 Tohoku Earthquake.

High-quality seismic data with broad frequency band are essential for the study of seismic energy, Es, and its scaling with seismic moment, Mo. The 2000 Western Tottori earthquake (Mw 6.6) and its aftershocks as recorded by NIED Hi-net including undistributed high-sampling data provide an excellent data set for this purpose. In this study we use: 1) regular data sampled at 100 sps of small and intermediate (M 2–4) aftershocks just after the mainshock, and, 2) 100 sps data and high sampling 1000 sps data of small events (M 0–3) about two years after the mainshock. Spectral ratios are calculated between all combinations of events that both occurred close to one another and had similar mechanisms. We calculated seismic energies of P and S waves for each event by fitting omega-square spectral models to the spectral ratios. Analysis of both the high and lower sampling rate data results in statistically significant size dependence of Es/Mo; however, none of these trends can explain overall scaling when all events, including the mainshock, are considered. Artificial size dependence due to band limitation and omega-square assumption may be responsible for the apparent trends.

We have elucidated depth variations in the stress field associated with the 2007 Noto Hanto, Japan, earthquake by stress tensor inversion using high‐quality aftershock data obtained by a dense seismic network. Aftershocks that occurred above 4 km in depth indicated a strike‐slip stress regime. By contrast, aftershocks in deeper parts indicated a thrust faulting stress regime. This depth variation in the stress regime correlates well with that in the slip direction derived from a finite source model using geodetic data. Furthermore, the maximum principal stress (σ1) axis was stably oriented approximately W20°N down to the depth of the mainshock hypocenter, largely in agreement with the regional stress field, but, below that depth, the σ1 axis had no definite orientation, indicating horizontally isotropic stress. One likely cause of these drastic changes in the stress regime with depth is the buoyant force of a fluid reservoir localized beneath the seismogenic zone.