Search

On February 12, 2014, an ML 4.2 earthquake (named 0212 EQ) occurred at the Tatun Volcano Group (TVG) in northern Taiwan. Seismic data from both the regional broadband array in Taiwan and a local dense seismic network in the TVG indicate that the 0212 EQ was located at a depth of 4.3 km near Shamoshan, one of the latest eruptive units in the TVG. Referring to the aftershocks surrounding the 0212 EQ, we further found that (1) they formed a linear seismic zone dipping ∼50° towards the southeast to reach the Shanchiao Fault that cuts through the TVG, (2) focal mechanism of more than half of the larger aftershocks (MD ≥ 1.5) matched this linear seismic zone for a normal faulting, and (3) the seismic b-value was calculated to be 3.55. They are used to delineate the geometry of the Shanchiao Fault at the upper crust and to infer that localized geothermal fluid might have affected the activity of aftershocks as manifested by the very high seismic b-value. The strongly compensated linear vector dipole of the centroid moment tensor suggests the 0212 EQ had a non-double-couple source resulted from the gravitational collapse. Thus, we conclude that the 0212 EQ was a combined synergistic effect of a collapsed reservoir filled with volcanic fluids and faulting activities along the Shanchiao Fault.

In this study, an empirical transfer function (ETF), which is the spectrum difference in Fourier amplitude spectra between observed strong ground motion and synthetic motion obtained by a stochastic point-source simulation technique, is constructed for the Taipei Basin, Taiwan. The basis stochastic point-source simulations can be treated as reference rock site conditions in order to consider site effects. The parameters of the stochastic point-source approach related to source and path effects are collected from previous well-verified studies. A database of shallow, small-magnitude earthquakes is selected to construct the ETFs so that the point-source approach for synthetic motions might be more widely applicable. The high-frequency synthetic motion obtained from the ETF procedure is site-corrected in the strong site-response area of the Taipei Basin. The site-response characteristics of the ETF show similar responses as in previous studies, which indicates that the base synthetic model is suitable for the reference rock conditions in the Taipei Basin. The dominant frequency contour corresponds to the shape of the bottom of the geological basement (the top of the Tertiary period), which is the Sungshan formation. Two clear high-amplification areas are identified in the deepest region of the Sungshan formation, as shown by an amplification contour of 0.5 Hz. Meanwhile, a high-amplification area was shifted to the basin’s edge, as shown by an amplification contour of 2.0 Hz. Three target earthquakes with different kinds of source conditions, including shallow small-magnitude events, shallow and relatively large-magnitude events, and deep small-magnitude events relative to the ETF database, are tested to verify site correction. The results indicate that ETF-based site correction is effective for shallow earthquakes, even those with higher magnitudes, but is not suitable for deep earthquakes. Finally, one of the most significant shallow large-magnitude earthquakes (the 1999 Chi-Chi earthquake in Taiwan) is verified in this study. A finite fault stochastic simulation technique is applied, owing to the complexity of the fault rupture process for the Chi-Chi earthquake, and the ETF-based site-correction function is multiplied to obtain a precise simulation of high-frequency (up to 10 Hz) strong motions. The high-frequency prediction has good agreement in both time and frequency domain in this study, and the prediction level is the same as that predicted by the site-corrected ground motion prediction equation.

In this study, microtremor array measurements were conducted at 45 sites in the Western Plain of Taiwan. The arrays were approximately 30 m or 60 m in radius, depending on the site. The maximum-likelihood frequency–wavenumber method was adopted to obtain the phase velocities of Rayleigh waves, and then a genetic algorithm technique based on an inversion scheme of the fundamental mode of the Rayleigh waves’ dispersion curves was applied to calculate a preliminary S-wave velocity (Vs) profile at each site. Because a layer of thick sediment covers the bedrock in the Western Plain of Taiwan, microtremor arrays in this size range cannot estimate the structure of the entire sediment. Therefore, this study implemented further inversion of the horizontal-to-vertical spectral ratios of the microtremors to estimate the deeper structures up to the bedrock of Vs greater than 1000 m/s. Previously logged velocity profiles for different depths at or near our study sites were collected and compared with the Vs profiles derived from our microtremor array measurements; the results were found to be highly comparable. Therefore, we could delineate the depth distributions for the layer depths for Vs = 600 m/s and 1000 m/s in this region. The depth for Vs = 600 m/s is approximately 50 m in the piedmont area and approximately 300 m at the coastline; moreover, the depths for Vs = 1000 m/s increase from 200 m in the piedmont area to approximately 1000 m at the coastline. The depths for Vs = 1.0 km/s (Z1.0), which is an important parameter that accounts for the basin effect in recent ground motion prediction equations, are consequently available at the study sites. The distribution of Z1.0 as a function of Vs30 indicates higher similarity to that in Japan than in the San Francisco Bay area.