Wen-Nan Wu

Could a Sumatra-like megathrust earthquake occur in the south Ryukyu subduction zone?

A comparison of the geological and geophysical environments between the Himalaya-Sumatra and Taiwan-Ryukyu collision-subduction systems revealed close tectonic similarities. Both regions are characterized by strongly oblique convergent processes and dominated by similar tectonic stress regimes. In the two areas, the intersections of the oceanic fracture zones with the subduction systems are characterized by trench-parallel high free-air gravity anomaly features in the fore-arcs and the epicenters of large earthquakes were located at the boundary between the positive and negative gravity anomalies. These event distributions and high-gravity anomalies indicate a strong coupling degree of the intersection area, which was probably induced by a strong resistance of the fracture features during the subduction. Moreover, the seismicity distribution in the Ryukyu area was very similar to the pre-seismic activity pattern of the 2004 Sumatra event. That is, thrust-type earthquakes with a trench-normal P-axis occurred frequently along the oceanward side of the mainshock, whereas only a few thrust earthquakes occurred along the continentward side. Therefore, the aseismic area located west of 128°E in the western Ryukyu subduction zone could have resulted from the strong plate locking effect beneath the high gravity anomaly zone. By analogy with the tectonic environment of the Sumatra subduction zone, the occurrence of a potential Sumatra-like earthquake in the south Ryukyu arc is highly likely and the rupture will mainly propagate continentward to fulfill the region of low seismicity (approximately 125° E to 129° E; 23° N to 26.5° N), which may generate a hazardous tsunami.

Possible northward extension of the Philippine Fault Zone offshore Luzon Island (Philippines)

The Philippine Fault Zone, a system of left-lateral strike-slip faults traversing the length of the Philippine Islands, is associated with the oblique convergence between the Philippine Sea Plate (PSP) and the Eurasian Plate (EP). Although it is a major deformational structure within the diffuse PSP–EP convergent boundary, some of its segments, particularly its marine extensions, are not well studied. To investigate the crustal deformation in the marine prolongation of the Philippine Fault Zone offshore Luzon Island, multi-channel seismic (MCS) data, gravity data and centroid moment tensor solutions were used in this study. Focal mechanism solutions from the Global CMT catalog were inverted to determine the average principal stress directions and consequently understand the prevailing stress regime in the study area. The stress inversion results indicate that the direction of maximum compression (σ1) is 321°N, which coincides with the PSP–EP convergence direction. From the MCS profiles, the study area was subdivided into deformation zone and a relatively stable zone. Thrust faulting, folding and general uplift are observed in the deformation zone. This zone is further subdivided into the active and inactive segments. In the active segment, uplift is occurring in the submarine ridge. This deformation pattern can be related to the ongoing uplift in some regions bisected by the PFZ. The inactive segment, characterized by intense folding of the sequences and faulting of the basement and overlying sequences, is suggested as the precursor of the Philippine Fault Zone. Deformation appears to be recently shifted to the east as delineated by an uplifted N-NW trending submarine ridge offshore NW Luzon Island.

Spatio-temporal distribution of seismic moment release near the source area of the 2011 Tohoku-Oki earthquake

To understand the generation mechanism of the 2011 Tohoku-Oki earthquake, we analyzed the spatial and temporal distribution of the cumulative seismic moment (ΣM0) for all earthquakes along the Japan Trench listed in the global centroid moment tensor catalog between January 1976 and November 2011. Three areas with distinct characteristics of ΣM0 are identified: (1) in the 2011 Tohoku-Oki source area, the ΣM0 released in the down-dip portion (≥30 km) was greater than that in the up-dip portion (<30 km) before the mainshock; (2) the Σ M0 of the up-dip portion in the region where slow slip activity prevails is greater than that of the down-dip portion throughout the study period; (3) in the surroundings of the source area, where interplate coupling is relatively low and the largest foreshock occurred, values of ΣM0 for the down-dip and up-dip portions are both intermediate. Our results show that a megathrust rupture could be generated by large accumulations of strain energy at the subduction interface, resulting from the differential strain energy released in the down-dip and up-dip portions during the interseismic interval. We propose that the variation pattern of ΣM0 may offer spatial constraints for seismic hazard assessment in the future.

Delineating seismogenic structures by a revised best estimate method: application to the Taiwan orogenic belt

The best estimate method (BEM), which shifts an object event towards the centroid of the events within its location uncertainty, was originally proposed to clarify the pattern of seismicity distribution. It was later widely applied to constrain the location and geometry of seismogenic structures. The goal of this study is to improve the efficiency, robustness and effectiveness of the BEM. We perform a series of synthetic experiments by generating a three-dimensional (3D) event dataset representing complex tectonic features and perturbing the presumed hypocenter locations. The synthetic experiments show that appropriate expression for the location uncertainty of the object event and restrictive thresholds of the event selection can significantly improve the collapsing seismicity image. We also demonstrate that the performance of our revised BEM (RBEM) is better than that of the original BEM. We then apply our RBEM to the seismicity catalog for the Taiwan region, collapsing the diffusive catalog hypocenters into sharp images of seismicity. Our RBEM result delineates a clearer low-seismicity zone in central Taiwan and a wider separation between the Wadati–Benioff double seismic layers in northeastern Taiwan than the corresponding features in either the catalog seismicity or the 3D relocation seismicity. Because of the simplicity and efficiency of the RBEM, it is perfectly suited to delineate the location and geometry of active faults in routine operations.