Nearby faults interact through stress changes induced by fault slip and viscoelastic flow. The process is, however, often elusive and can be geometry-dependent and time-variant. Here, we combine geodetic and field observations to characterize the interaction of two head-to-head, conjugate faults in eastern Taiwan during the 2022 Chihshang earthquake sequence. We map the coseismic slip on the Central Range fault and dynamically-triggered shallow slip on the Longitudinal Valley fault, which has been creeping interseismically. Overlapping of seismic and aseismic slip suggests that the Longitudinal Valley fault is capable of hosting a variety of distinct slip behaviors. Moreover, substantial slip on the Central Range fault suppresses Coulomb stress on the Longitudinal Valley fault, and vice versa, resulting in seismic bursts in an out-of-phase pattern on the two faults as seen in the hundred-year historical records. Such fault interaction implies the need for time-dependent seismic hazard reassessment for the complex fault system.

The ascending and descending InSAR deformations derived from ALOS-2 and Sentinel-1 satellite SAR images and GPS displacements are used to estimate the fault model of the 2018 Mw 6.4 Hualien earthquake. The sinistral strike-slip fault dipping to the west with a high dip angle of 89.4° and a rake angle of 201.7° is considered as the seismogenic fault of this event. This seismogenic fault also triggered the ruptures of the Milun fault, which dips to the east with a dip angle of ~72°, and an unknown west-dipping fault with a dip angle of 85.2°. Two predicted faulting models indicate that the InSAR deformation fields include more postseismic slip than those of the GPS data. The north segment of the Milun fault and west-dipping fault have been triggered by the rupture of the seismogenic fault, but the postseismic slip occurred only in the south segment of the Milun fault. The InSAR-derived co-seismic and postseismic faulting model suggests that the significant slip concentrates at depths of 2.4–15.0 km of the main fault, 0.0–14.0 km of the Milun fault. Only minor slip is detected on the west-dipping fault. The maximum fault slip of ca. 2.1 m is located at the depth of ca. 2.4 km under the Meilun Tableland. The Coulomb failure stress (CFS) change calculated by the co-seismic and postseismic faulting model shows that there is a significant CFS increase in the east of the south segment of the Milun fault, but few of the aftershocks occur in this area, which indicates a high risk of future seismic hazard.

The earthquake doublet of Mw 6.5 and 6.9 occurred along the west-dipping Central Range fault (CRF) adjacent to the east-dipping creep segment of the Longitudinal Valley fault in eastern Taiwan on 17–18 September 2022. The faulting model derived from the Interferometric Synthetic Aperture Radar and Global Positioning System observation suggests that the west-dipping CRF rupture is responsible for the two mainshocks. Meanwhile, the two major earthquakes resulted in ∼100 km of fault slip along the double-convergent plate boundary faults near the Longitudinal Valley in eastern Taiwan. The 2022 Chihshang earthquake sequence filled the seismic gap of the CRF located between the 2006 ML 6.1 Taitung earthquake and the 2013 ML 6.1 Ruisui earthquake. Finally, the 2022 Chihshang earthquake sequence increased the Coulomb failure stress in the southernmost segment of the CRF, which may cause a higher seismic potential in the future.

Applications of low-cost single-frequency continuous GPS receivers for monitoring volcano and landslide activities as well as to complement dual-frequency receivers have been demonstrated to produce stable and accurate positioning. In studies of crustal deformation, the relative distance between monitoring stations may vary from several kilometers to tens of kilometers, hence the differential single-frequency observations cannot model the ionospheric delay or other distance dependent errors. The 55 low-cost single-frequency continuous stations have been deployed together with 52 continuous dual-frequency stations in southeastern Taiwan since 2008. All of the single-frequency stations have applied corrections using dual-frequency stations to eliminate the distance dependent errors. Comparing velocity estimates from 8 co-located, the differences in horizontal and vertical components are less than 3 mm/yr and 6 mm/yr, respectively. Our study shows that the combination of single- and dual-frequency GPS data can provide robust results to study the fault slip behavior on the Longitudinal Valley fault.

Two subsequent strong earthquakes—the ML 6.6 Guanshan and the ML 6.8 Chihshang earthquakes—struck southeast Taiwan on 17 and 18 September 2022, causing extensive damage associated with multiple fault systems. Mapping their coseismic deformation is crucial for constraining the slip behavior and possible fault interaction. To enhance the coverage and resolution of the coseismic deformation features for the two events, data from two dense strong-motion networks are utilized to retrieve coseismic displacements using a two-stage baseline correction method. The reliability of derived displacements is validated with cosite Global Navigation Satellite Systems (GNSS) data. The denser distribution of the seismic stations supplements the GNSS and seismic data previously reported and reveals a more detailed view of regional and near-fault deformation. The integrated results confirm that the surface ruptures caused by the Chihshang earthquake are located to the east of the responsible Central Range fault, north of the epicenter. The complex coseismic deformation pattern of the Guanshan earthquake indicates the possible involvement of local structures especially around Taitung. More critically, localized subsidence created by the southward rupture propagation of the Guanshan earthquake is depicted between the epicenters of the two earthquakes. This subsidence may impart large stress on the northern segment of the fault and trigger the Chihshang earthquake. The enriched near-fault deformation patterns revealed by our newly compiled dataset offer valuable constraints for understanding shallow slip distribution, shallow versus deep slip budget, and possible multifault interaction in future investigations.