Search

Dense strong motion observation networks provided us with valuable data for studying strong motion generation from large earthquakes. From kinematic waveform inversion of seismic data, the slip distribution on the fault surface of large earthquakes is known to be spatially heterogeneous. Because heterogeneities in the slip and stress drop distributions control the generation of near-source ground motion, it is important to characterize these heterogeneities for past earthquakes in constructing a source model for reliable prediction of strong ground motion. The stress changes during large earthquakes on the faults recently occurring in Japan are estimated from the detailed slip models obtained by the kinematic waveform inversion. The stress drops on and off asperities are summarized on the basis of the stress change distributions obtained here. In this paper, we define the asperity to be a rectangular area whose slip is 1.5 or more times larger than the average slip over the fault according to the previous study for inland crustal earthquakes. The average static stress drops on the asperities of the earthquakes studied here are in the range 6-23 MPa, whereas those off the asperities are below 3 MPa. We compiled the stress drop on the asperities together with a data set from previous studies of other inland earthquakes in Japan and elsewhere. The static stress drop on the asperity depends on its depth, and we obtained an empirical relationship between the static stress drop and the asperity's depth. Moreover, surface-breaking asperities seemed to have smaller stress drops than buried asperities. Simple ground motion simulations using the characterized asperity source models reveal that deep asperities generate larger ground motion than shallow asperities, because of the different stress drops of the asperities. These characteristics can be used for advanced source modeling in strong ground motion prediction for inland crustal earthquakes. [ABSTRACT FROM AUTHOR]

The paper addresses the interpretation of the location, type, and size of the source for the earth-quake of March 11, 2011. The source-a subvertical reverse fault trending in the azimuth of ∼25° along the island arc-is located in the middle part of the Pacific slope of Honshu Island, between 38°-38.5°N and 35.5°N. The length of the source, about 350 km, approximately corresponds to a magnitude ∼8.7 earthquake. In the north, the source is bounded by a sublatitudinal reverse fault, which generated an earthquake with magnitude 7.2-7.5 in 1978. On this segment of the Pacific slope of Honshu Island, there are probably another one or a few other large seismic sources, which are still latent. They are longitudinal reverse faults, which are comparable in scale with the source of the March, 2011 earthquake. The recurrence period of the maximal earthquakes in such sources is more than 1000 years. [ABSTRACT FROM AUTHOR]

The 1703 Genroku Kanto earthquake and the resulting tsunami caused catastrophic damage in the Kanto region of Japan. Previous modeling of the 1703 earthquake applied inversion analyses of the observed terrestrial crustal deformations along the coast of the southern Boso Peninsula and revealed that the tsunami was generated along the Sagami Trough. Although these models readily explained the observed crustal deformation, they were unable to model an offshore fault along the Sagami Trough because of difficulties related to the distance of the offshore fault from the shoreline. In addition, information regarding the terrestrial crustal deformation is insufficient to constrain such inverted models. To model an offshore fault and investigate the triggering of large tsunamis off the Pacific coast of the Boso Peninsula, we studied historical documents related to the 1703 tsunami from Choshi City. Based on these historical documents, we estimated tsunami heights of ≥5.9, 11.4-11.7, ≥7.7, 10.8 and ≥4.8 m for the Choshi City regions of Isejiga- ura, Kobatake- ike, Nagasaki, Tokawa and Na'arai, respectively. Although previous studies assumed that the tsunami heights ranged from 3.0 to 4.0 m in Choshi City, we revealed that the tsunami reached heights exceeded 11 m in the city. We further studied the fault model of the 1703 Genroku Kanto earthquake numerically using the newly obtained tsunami height data. Consequently, we determined that the source of the 1703 earthquake was a 120-km-long offshore fault along the Sagami Trough, which is in close proximity to the Japan Trench. Our results suggest that earthquake energies resulting in magnitudes greater than Mw 8.32 along the entire length of the Sagami Trough could have been released during the 1703 Genroku Kanto earthquake. [ABSTRACT FROM AUTHOR]

Understanding the effect of the complex fault geometry on the dynamic rupture process and discriminating it from the complexity originating from the rheological properties of faults, is an essential subject in earthquake science. The 2014 Northern Nagano earthquake, which occurred near the end zone of a major active fault system, provided unique observations that enabled us to investigate both the detailed geometrical fault structure and surface deformation patterns as well as the temporal sequence led up from a prominent foreshock activity. We first develop a geometrical fault model with a substantial variation along strike, and a model for the regional stress field, which is well constrained by the observations. This significant along-strike variation in fault geometry probably reflects the difference of fault maturity at the end zone of the complex fault system. We used this model in order to conduct a set of dynamic rupture simulations using the highly efficient spatiotemporal boundary integral equation method. Based on our simulations, we show that the observed surface deformation can be reasonably explained as the effect of the non-planar fault geometry with a number of branch faults and bends. [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]

Large earthquakes occurred in the central part of Kumamoto Prefecture on April 14-16, 2016, causing severe damage to the northern segment of the Hinagu faults and the eastern segment of the Futagawa faults. Earthquake surface ruptures appeared along these faults and on the Aso volcanic edifice, which in turn generated landslides. We conducted landform change analysis of the central cones of Aso volcano by using satellite and aerial photographs. First, we categorized the topographical changes as surface scarps, arc-shaped cracks, and linear cracks. Field survey indicated that landslides caused the scarps and arc-shaped cracks, whereas faulting caused the linear cracks. We discovered a surface rupture concentration zone (RCZ) formed three ruptures bands with many surface ruptures and landslides extending from the west foot to the center of the Aso volcanic edifice. The magmatic volcanic vents that formed during the past 10,000 years are located along the north margin of the RCZ. Moreover, the distribution and dip of the core of rupture concentration zone correspond with the Nakadake craters. We conclude that a strong relationship exists between the volcanic vents and fault structures in the central cones of Aso volcano. [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]

The 2016 Kumamoto earthquake sequence started with a M (Japan Meteorological Agency magnitude) 6.5 event on April 14, and culminated in a M 7.3 event on April 16. Associated with the sequence, approximately 34-km-long surface ruptures appeared along the eastern part of the Futagawa fault zone and the northernmost part of the Hinagu fault zone. We carried out an urgent field investigation soon after the earthquake to map the extent and displacement of surface ruptures with the following results. (1) The rupture zone generally consisted of a series of left-stepping en echelon arrays of discontinuous fault traces of various lengths. (2) Slip exceeding 100 cm occurred on previously unrecognized fault traces in the alluvial lowland of the Kiyama plain and on the western rim of the Aso volcano caldera. (3) Large slip with maximum dextral slip of 220 cm was measured throughout the central section of the rupture zone along the Futagawa segment, and the slip gradually decreased bilaterally on the adjoining northeastern and southwestern sections. (4) The surface rupture mostly occurred along fault traces mapped in previous active fault investigations. (5) Most of the surface ruptures were produced by the mainshock, and significant postseismic slip occurred after the mainshock. [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]

In 1995, the Hyogoken-Nanbu earthquake ( M 7.3) ruptured the Nojima fault, Awaji Island, central Japan. To investigate the recovery process of a fault zone after a large earthquake, repeated water injection experiments have been conducted every few years in an 1800-m-long borehole near the Nojima fault since 1997. In addition, the groundwater discharge rate and pressure have been observed in an 800-m borehole. From the resulting data, the macroscopic permeability of the fault fracture zone was estimated to range roughly from 1 × 10 to 2 × 10 m/s. The macroscopic permeability of the fault fracture zone decreased until 2003, and then, it stabilized or increased slightly through 2006. These changes in permeability indicate that the fault fracture zone stabilized within 8 years after the occurrence of the earthquake. [ABSTRACT FROM AUTHOR]

The 2016 Kumamoto earthquakes, including an Mw-7 right-lateral earthquake on April 15 (UTC), occurred along faults within the Beppu−Shimabara graben in central Kyushu, Japan. Previous studies showed that the graben area was under heterogeneous stress conditions with north-south T-axes and spreading in a north-south direction. Here, we construct a detailed crustal deformation model using three-component Global Navigation Satellite System data in 2000-2010 and considering the distribution of geological fault traces in this area. Our inversion analysis suggests that the strain accumulation rate for the right-lateral seismic slip segment (corresponding to the Futagawa fault), where the largest of the 2016 Kumamoto earthquakes ruptured, was several times smaller than the other segments in the Beppu−Shimabara graben. Furthermore, we observe distinct subsidence along the Beppu−Shimabara graben. Our base model attributes the subsidence to deflation of magma reservoirs beneath volcanoes, but the observed vertical velocities are poorly fit. In order to improve the fitting results for the vertical deformation, we need more sophisticated volcano-deformation model (such as a sill-like deformation source for Mt. Aso) or graben model. [ABSTRACT FROM AUTHOR]

The M7.3 Kumamoto earthquake occurred on April 16, 2016, in the western part of Kyushu, at a depth of 12 km, on an active strike-slip fault. Here, we report on a relatively widespread activation of small remote earthquakes, which occurred as far as Hokkaido, detected by analyzing the continuous waveform data recorded at seismic stations all over Japan. Such relatively widespread remote seismicity activation, following a large inland earthquake, has not been reported before for Japan. Our analysis demonstrates that the remote events were triggered dynamically, by the passage of the surface waves from the Kumamoto earthquake. Most of the remotely triggered events in the Tohoku and Hokkaido regions, as well as close to Izu Peninsula, occur at or close to volcanoes, which suggests that the excitation of crustal fluids, by the passage of Rayleigh waves, played an important triggering role. Nevertheless, remote activation in other regions, like Noto Peninsula, occurred away from volcanoes. The relatively large-amplitude Love waves, enhanced by a source directivity effect during the Kumamoto earthquake, may have triggered seismicity on local active faults. The dynamic stresses in the areas where remote activation has been observed range from several kPa to tens of kPa, the thresholds being lower than in previous dynamic triggering cases for Japan; this might relate to a change in the crustal conditions following the 2011 M9.0 Tohoku-oki earthquake, in particular at volcanoes in NE Japan. [ABSTRACT FROM AUTHOR]

We estimate corner frequencies and stress drops for 298 events ranging from M 3.2-7.0 in 17 inland crustal earthquake sequences in Japan to investigate the source scaling and variation in stress drops. We obtain the source spectral ratio from observed records by the S-wave coda spectral ratio method. The advantage of using the S-wave coda is in obtaining much more stable source spectral ratios than using direct S-waves. We carefully examine the common shape of the decay of coda envelopes between event pair records. The corner frequency and stress drop are estimated by modeling the observed source spectral ratio with the omega-square source spectral model. We investigate the dependences of stress drops on some tectonic effects such as regionality, focal mechanism, and source depth. The principal findings are as follows: (1) a break in self-similar source scaling is found in our dataset. Events larger than M 4.5 show larger stress drops than those of smaller events. (2) Stress drops of aftershocks are mostly smaller than those of mainshocks in each sequence. (3) There are no systematic differences between stress drops of events occurring inside and outside the Niigata-Kobe Tectonic Zone in Japan. (4) Clear dependence of the faulting type on stress drops cannot be seen. (5) Stress drops of aftershocks depend on their source depth. (6) The crack size obtained from the corner frequency corresponds to the total rupture area of heterogeneous slip models for large events. [ABSTRACT FROM AUTHOR]

The 1995 Kobe (Hyogo-ken Nanbu) earthquake, M = 7.2, ruptured the Nojima fault in southwest Japan. We have studied core samples taken from two scientific drillholes that crossed the fault zone SW of the epicentral region on Awaji Island. The shallower hole, drilled by the Geological Survey of Japan (GSJ), was started 75 m to the SE of the surface trace of the Nojima fault and crossed the fault at a depth of 624 m. A deeper hole, drilled by the National Research Institute for Earth Science and Disaster Prevention (NIED) was started 302 m to the SE of the fault and crossed fault strands below a depth of 1140 m. We have measured strength and matrix permeability of core samples taken from these two drillholes. We find a strong correlation between permeability and proximity to the fault zone shear axes. The half-width of the high permeability zone (approximately 15 to 25 m) is in good agreement with the fault zone width inferred from trapped seismic wave analysis and other evidence. The fault zone core or shear axis contains clays with permeabilities of approximately 0.1 to 1 microdarcy at 50 MPa effective confining pressure (10 to 30 microdarcy at in situ pressures). Within a few meters of the fault zone core, the rock is highly fractured but has sustained little net shear. Matrix permeability of this zone is approximately 30 to 60 microdarcy at 50 MPa effective confining pressure (300 to 1000 microdarcy at in situ pressures). Outside this damage zone, matrix permeability drops below 0.01 microdarcy. The clay-rich core material has the lowest strength with a coefficient of friction of approximately 0.55. Shear strength increases with distance from the shear axis. These permeability and strength observations reveal a simple fault zone structure with a relatively weak fine-grained core surrounded by a damage zone of fractured rock. In this case, the damage zone will act as a high-permeability conduit for vertical and horizontal flow in the plane of the fault. The fine-grained core region, however, will impede fluid flow across the fault. [ABSTRACT FROM AUTHOR]

The 2004 Mid Niigata Prefecture earthquake ( M JMA 6.8) and its aftershock sequences generated complicated, i.e., several conjugate fault planes in their source region. In order to understand the generating process of these earthquakes, we estimated a 3-D distribution of relative scattering coefficients in the source region. The large slip area during the main shock rupture seems to be bounded by strong heterogeneous zones with larger scattering coefficients. Hypocenters of the main shock and major large aftershocks with M 5-6 classes tend to be located close to stronger scattering areas. We found that one of these strong heterogeneities already existed before the occurrence of the M 5.9 aftershock on November 8. We suppose that heterogeneous structures in the source region of this earthquake sequence affected the initiation and growth of ruptures of the main shock and major large aftershocks. [ABSTRACT FROM AUTHOR]

Non-volcanic tremor is a weak, extended duration seismic signal observed episodically on some major faults, often in conjunction with slow slip events. Such tremor may hold the key to understanding fundamental processes at the deep roots of faults, and could signal times of accelerated slip and hence increased seismic hazard. The mechanism underlying the generation of tremor and its relationship to aseismic slip are, however, as yet unresolved. Here we demonstrate that tremor beneath Shikoku, Japan, can be explained as a swarm of small, low-frequency earthquakes, each of which occurs as shear faulting on the subduction-zone plate interface. This suggests that tremor and slow slip are different manifestations of a single process. [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]

The 2011 Tohoku-oki earthquake unexpectedly ruptured to the shallowest portion of the plate boundary fault and triggered a large tsunami. The shallow portion had generally been regarded as a seismically stable zone until this event, but its significance has now been dramatically revealed for future disaster mitigation. This research approaches the shallow portion, especially the formation process of its structure and plate boundary faults. Scientific drilling conducted near the Japan Trench after the earthquake reported a thin plate boundary fault (~ 7 m) and thick deformation zone (~ 100 m). This thin fault would be expected given the relatively small displacement near the trench (~ 3.2 km), but the deformation zone thickness is anomalously wide given this small magnitude of slip. To understand the dynamic deformation processes that lead to the development of a thick deformation zone surrounding a thin fault core, we conducted forward modeling of an analog experiment with the technique to visualize fault activity. Sandbox experiments are effective for the approximation of the geological phenomenon and structure. The seismic profile of the largest slip region in the 2011 earthquake shows that a half-graben structure has been subducted underneath the frontal wedge, thus we focused on this structural architecture. As a result, we found a new fault formation pattern, i.e., the frontal thrust (the most frontal part of décollement) periodically partitioned into pieces, which connect again to form a large-continuous fault. The fault also oscillates up and down during this process, which we call “dancing,” and a thick shear zone is formed in a relatively short time where this occurs, even though the fault only has a small displacement. By analogy, the thick deformation zone observed at the Japan Trench could be formed by such fault dancing. The energy of the fault activity is commonly estimated from the fault displacement derived from the thickness of the shear zones. Applying the thickness-displacement law without considering the effect of the dancing may cause overestimation. The architecture of the shear zone formed is similar to that of a mélange, and the origin of tectonic mélanges may be explained by this mechanism. [ABSTRACT FROM AUTHOR]