Search

An Mw 7.5 earthquake struck Palu in the northern coast of Sulawesi island, Indonesia, on 28 September 2018. Its focal mechanism was determined to be a left‐lateral strike‐slip fault, which is generally expected to not produce a tsunami. However, a large tsunami with runup heights of more than 6 m was observed along the coast of Palu city. Here, we show a complex triggering supershear source model as determined by teleseismic waveform inversion. Three asperities with different slip characteristics were found on the 120‐kilometer‐long rupture zone. Significant triggering rupture with a supershear speed was observed south of the epicenter, which was just beneath Palu city. This special rupture process can cause a strong directivity effect that produced anomalously large ground shaking with nonlinear effects in Palu area. The coseismic deformation determined from the inverted source model showed large horizontal displacements. These horizontal movements combined with complex bathymetry and topography could have pushed seawater to generate a tsunami even though the Palu earthquake was a strike‐slip event.

The 5 February 2016, Meinong, Taiwan, earthquake brought extensive damage to nearby cities with significant pulse‐like velocity ground motions. In addition to the spatial slip distribution determination using filtered strong‐motion data, we show that, with the advantage of the densely distributed seismic network as a seismic array, we can project the earthquake sources (asperities) directly using nearly unfiltered data, which is crucial to the understanding of the generation of the pulse‐like velocity ground motions. We recognize that the moderate but damaging M_L 6.6 Meinong earthquake was a composite of an M_w 5.5 foreshock and an M_w 6.18 mainshock with a 1.8–5.0 s time delay. The foreshock occurred at the hypocenter reported by the official agency, followed by the mainshock with a centroid located at 12.3 km to the north‐northwest of the hypocenter and at a depth of 15 km. This foreshock–mainshock composition is not distinguishable in the finite‐fault inversion because it filtered the seismic data to low frequencies. Our results show that the pulse‐like velocity ground motions are mainly attributed to the source of mainshock with its directivity and site effects, resulting in the disastrous damages in the city of Tainan. Although finite‐fault inversion using filtered seismic data for spatial slip distribution on the fault has been a classic procedure in understanding earthquake rupture processes, using a dense seismic network as a seismic array for unfiltered records helps us delineate the earthquake sources directly and provide more delicate information for future understanding of earthquake source complexity.

Despite a moderate magnitude, M_w = 6.4, the 5 February 2016 Meinong, Taiwan, earthquake caused significant damage in Tainan City and the surrounding areas. Several seismograms display an impulsive S-wave velocity pulse with an amplitude of about 1 m s-1, which is similar to large S-wave pulses recorded for the past several larger damaging earthquakes, such as the 1995 Kobe, Japan, earthquake (M_w = 6.9) and the 1994 Northridge, California, earthquake (M_w = 6.7). The observed PGV in the Tainan area is about 10 times larger than the median PGV of M_w = 6.4 crustal earthquakes in Taiwan. We investigate the cause of the localized strong ground motions. The peak-to-peak ground-motion displacement at the basin sites near Tainan is about 35 times larger than that at a mountain site with a similar epicentral distance. At some frequency bands (0.9 - 1.1 Hz), the amplitude ratio is as large as 200. Using the focal mechanism of this earthquake, typical “soft” and “hard” crustal structures, and directivity inferred from the observed waveforms and the slip distribution, we show that the combined effect yields an amplitude ratio of 17 to 34. The larger amplitude ratios at higher frequency bands can be probably due to the effects of complex 3-D basin structures. The result indicates that even from a moderate event, if these effects simultaneously work together toward amplifying ground motions, the extremely large ground motions as observed in Tainan can occur. Such occurrences should be taken into consideration in hazard mitigation measures in the place with frequent moderate earthquakes.

Although surface geology, eruption information and clustering seismicity all suggest Turtle Island (Kueishantao) of northern Taiwan is an active volcano, there was no direct evidence to conclude that magma reservoirs exist beneath it. Even less evidence is available to determine their spatial configuration. If the magma reservoirs are filled by liquids and melt, S-waves are totally reflected and leave behind a shadow, like when passing through the Earth’s outer core. We detect both these S-wave shadows and strong reflections from the surface using earthquakes at different depths and azimuths. These observations identify a km-scale molten-filled volume located beneath Turtle Island. The magmatic nature of the reservoir is supported by the onset of non-double-couple earthquakes with strong CLVD (Compensated Linear Vector Dipole) and ISO (Isotropic) components, which show a tensor crack compatible with some volume changes within the reservoir. Combining these results with two independent 3-D velocity models and aeromagnetic anomalies recorded in Taiwan, a partially-molten ~19% low-velocity volume is estimated in the mid-crust (13–23 km), with spatial uncertainties of ~3 km. The elongated direction approximately follows the strike of the Okinawa trough, indicating that the source of the magma reservoir might be a back-arc opening.

On 6 February 2016, an M L 6.6 earthquake occurred in the Meinong area of southern Taiwan, causing anomalously large ground shaking. In this study, a joint source inversion was performed to understand the rupture process of this event and to identify how the anomalous strong motion was produced. The results show that the rupture process of this event was complex; at least two asperities were identified on the fault plane. The rupture mainly developed in the down‐dip direction and propagated toward the northwest. Results from a 3D wave propagation simulation further indicate that the strong ground shaking observed in southwest Taiwan was caused by a combination of three effects: (1) rupture directivity, (2) source radiation pattern, and (3) sedimentary amplification. The results demonstrate that constructive source rupture and wave propagation effects due to a moderate earthquake can cause anomalously strong ground motion and serious seismic hazards in nearby densely populated areas.

The 1909 Taipei earthquake (M 7.3) occurred beneath the Taipei metropolitan area (TMA) causing substantial damage according to the historical literature. According to the hypocenter relocation and tectonic implications provided in a previous study, we simulated ground motions within the TMA using a hybrid simulation method involving the spectral-element method (SEM) and the empirical Green's function method (EGFM). We used the SEM for simulating low-frequency components and the EGFM for simulating high-frequency components. These high and low frequency components were subsequently combined. For the EGFM we used the records from a recent M_L 4.9 earthquake (11 October 2013, depth = 143.8 km) in the Taipei area as the empirical Green's function. According to the historical literature, the observed PGA (peak ground acceleration) values are 59.2 and 67.0 gal at ancient stations TAP and KEE, with periods of 1.21 and 1.34 s, respectively. By comparing the simulated PGA values at modern stations TAPB and WFSB to the historical documented ones for 12 different models, our result suggests that the 1909 Taipei earthquake was an event with a magnitude of about M_w 7.3 and stress drop of approximately 30 bars, or a smaller equivalent magnitude between M_w 6.8 - 7.3 but with much higher average stress drop of more than 100 bars. For a deep event beneath TMA a larger vertical P-wave motion and longer period shaking wave, as addressed in the historical literature, might be expected with prolonged shaking as found in the simulation. A seismic hazard assessment is necessary for metropolitan Taipei to better understand the long period shaking from deep subduction zone intra plate events.

The Mw8.4 Illapel earthquake occurred on 16 September was the largest global event in 2015. This earthquake was not unexpected because the hypocenter was located in a seismic gap of the Peru-Chile subduction zone. However, the source model derived from 3-D spectral-element inversion of teleseismic waves reveals a distinct two-stage rupture process with completely different slip characteristics as a composite megathrust event. The two stages were temporally separated. Rupture in the first stage, with a moment magnitude of Mw8.32, built up energetically from the deeper locked zone and propagated in the updip direction toward the trench. Subsequently, the rupture of the second stage, with a magnitude of Mw8.08, mainly occurred in the shallow subduction zone with atypical repeating slip behavior. The unique spatial-temporal rupture evolution presented in this source model is key to further in-depth studies of earthquake physics and source dynamics in subduction systems.

An Mw 6.2 earthquake struck Hualien on the eastern coast of Taiwan on April 18, 2019. Its focal mechanism was determined to be a thrust fault with a depth of 18.8 km. A wide range of strong ground shaking with an intensity value from 4 (25–80 Gal) to 7 (greater than 400 Gal) occurred in northern Taiwan, which is unusual for earthquakes of this size. Here we show a source model as determined by local seismic waveform inversion. Two asperities were found, the smaller one occurred at the hypocenter and the larger one was found 10 km north of the hypocenter with a depth of 18–25 km. The initial rupture started from the hypocenter and then rapidly propagated all the way to the large asperity in the north. A near supershear speed of close to 4.0 km/s was found to have occurred during this northward rupture process. The inverted source model was evaluated using forward three-dimensional ground motion simulations. The results show strong agreement with island-wide observed displacement waveforms and peak ground acceleration for periods from 3.0 to 20.0 s. The rapid northward rupture caused a strong directivity effect coupled with the specific source radiation pattern, resulting in a large area of strong ground shaking in northern Taiwan, even though the Hualien earthquake was a moderately-sized event.

The southern Ryukyus represents an area where different tectonic stress regimes result in high seismicity and increased seismic hazard for nearby areas such as Taiwan. On 18 December 2001 at 04:03 (GMT) a strong earthquake (Mw 6.8) occurred in the forearc area of the southern Ryukyu subduction zone. Revised moment tensor solutions published by GCMT and BATS groups show a normal faulting mechanism with some strike-slip component and also point to a shallow focal depth (∼12 km). We use arrival times picked at both Taiwanese and Japanese stations along with a 3D geo-realistic a priori velocity model in order to obtain accurate absolute locations for the mainshock and 153 of its aftershocks. Locations are derived by using the Maximum intersection (MAXI) algorithm which has been used in many previous seismicity studies in the southern Ryukyus. These improved locations indicate that the mainshock was caused by the failure of a NE–SW oriented fault that extends from the edge of the Nanao forearc sedimentary basin to the Ryukyu arc basement. Far-field P and SH waveforms of the mainshock recorded at stations surrounding the source and at distances 30–100°, were inverted for the purpose of investigating its rupture process. A non-negative least-squares inversion technique utilizing multiple time windows was used to derive the spatio-temporal slip distribution. The preferred slip distribution model shows that there is one large area of high slip (∼0.9 m) at 5–15 km depth that essentially represents the crystalline rocks of the Ryukyu arc basement. Another smaller area with lower slip (∼0.4 m) extends at 10–15 km depth beneath the Nanao basin. Most aftershocks are located in areas of low slip (

Topography influences ground motion and, in general, increases the amplitude of shaking at mountain tops and ridges, whereas valleys have reduced ground motions, as is observed from data recorded during and after real earthquakes and from numerical simulations. However, recent publications have focused mainly on the implications for ground motion in the mountainous regions themselves, whereas the impact on surrounding low-lying areas has received less attention. Here, we develop a new spectral-element mesh implementation to accommodate realistic topography as well as the complex shape of the Taipei sedimentary basin, which is located close to the Central Mountain Range in northern Taiwan. Spectral-element numerical simulations indicate that high-resolution topography can change peak ground velocity (PGV) values in mountainous areas by ±50% compared to a half-space response. We further demonstrate that large-scale topography can affect the propagation of seismic waves in nearby areas. For example, if a shallow earthquake occurs in the I-Lan region of Taiwan, the Central Mountain Range will significantly scatter the surface waves and will in turn reduce the amplitude of ground motion in the Taipei basin. However, as the hypocenter moves deeper, topography scatters body waves, which subsequently propagate as surface waves into the basin. These waves continue to interact with the basin and the surrounding mountains, finally resulting in complex amplification patterns in Taipei City, with an overall PGV increase of more than 50%. For realistic subduction zone earthquake scenarios off the northeast coast of Taiwan, the effects of topography on ground motion in both the mountains and the Taipei basin vary and depend on the rupture process. The complex interactions that can occur between mountains and surrounding areas, especially sedimentary basins, illustrate the fact that topography should be taken into account when assessing seismic hazard.

International audience; The Taiwanese orogen is the result of the active and vigorous collision between the Eurasian and the Philippine Sea plates. In this strong convergence context, a one month-long earthquake sequence composed of two mainshocks of magnitude 6 occurred in April 2006 in the Taitung area. The first mainshock (T1) mobilized a structure just west to the plate boundary, while the second one (T2) is located on the other side. In order to retrieve the exact fault geometry of T2, we inverted waveforms from stations located at local and teleseismic distances. We performed a grid search on strike and dip (each strike-dip couple defining a geometry) considering the source as a point source then as an extended source. In the last case, a simple average over stations misfit failed to isolate a precise geometry. To overcome this problem, we used a more statistical approach on the distribution of misfits to define the best geometry. Comparing the best model to local structures, it appears the generative fault was the plate boundary that rotates from a strike pointing at N20°E north of the event to N0°E in T2 area with an identical eastward dip of 35°. For this model, the fault slip inversion provides a critical slip of 3.5 cm above which slip, rake and rupture time are constrained with uncertainties of 29%, 14° and 0.47 s respectively. The average slip along the rupture was 20 cm with a maximum of 46 ± 13 cm. The movement was inverse with a minor left-lateral component similar to the faulting behavior of the plate boundary. In addition, the slip pattern of T2 is contained within the southern portion of the deepest segment of the plate boundary and at the edge of the rupture area of the 2003 Chengkung earthquake (MW 6.8), a large event also generated by the plate boundary but 2.5 years earlier. After 1 s of aseismic spreading, the rupture propagated seismically and circularly outward before being stopped by the fault bending.

The 2013 M_L 6.4 Ruisui earthquake struck a seismic gap around the northern part of the Longitudinal Valley where rare earthquakes have occurred for at least two decades. Seismic data with different frequency bands connote a diverse scale of source characteristics. The slip model derived by a previous study using long-period seismic data provided us with the preliminary earthquake source properties in low-frequency energy radiation. However, the high-frequency part of the seismic energy needs to be considered for the strong motions in the comprehensive frequency band. As a case study, we carried out broadband strong motion simulations through the hybrid scheme, which simulated waveforms combining the low-frequency motions (frequency-wavenumber integration method) and high-frequency motions (empirical Green's function method) to reveal the strong ground motion and source characteristics of the 2013 Ruisui earthquake. The results show that even though the local structure is complicated this hybrid simulation can effectively rebuild the overall broadband strong ground motion characteristics for the 2013 Ruisui earthquake. Additional study issues are required to better understand the nature of localized high peak ground acceleration and the related seismic hazards for future earthquakes.

We use the spectral-element method to simulate strong ground motion throughout the Taipei metropolitan area. Mesh generation for the Taipei basin poses two main challenges: (1) the basin is surrounded by steep mountains, and (2) the city is located on top of a shallow, low-wave-speed sedimentary basin. To accommodate the steep and rapidly varying topography, we introduce a thin high-resolution mesh layer near the surface. The mesh for the shallow sedimentary basin is adjusted to honor its complex geometry and sharp lateral wave-speed contrasts. Variations in Moho thickness beneath Northern Taiwan are also incorporated in the mesh. Spectral-element simulations show that ground motion in the Taipei metropolitan region is strongly affected by the geometry of the basin and the surrounding mountains. The amplification of ground motion is mainly controlled by basin depth and shallow shear-wave speeds, although surface topography also serves to amplify and prolong seismic shaking.