Bor-Shouh Huang

Underplating in the Himalaya-Tibet Collision Zone Revealed by the Hi-CLIMB Experiment

Himalayan-Tibetan Underplate The Himalayas formed from the collision of India with Eurasia beginning about 50 million years ago, but the fate and position of the subducted Indian crust was not well defined until the Hi-CLIMB seismic experiment was initiated. The centerpiece of the project is an 800-kilometer-long, closely spaced, linear array of broadband seismographs, extending from the Ganges lowland, across the Himalayas, and onto the central Tibetan plateau. Nábělek et al. (p. 1371) present images of the crust and upper mantle of the Southern Tibetan plateau underthrust northward by the Indian plate, in which they trace the base of the Indian plate to 31°N. The character of the crust-mantle interface in this region suggests that the Indian crust is at least partly decoupled from the mantle beneath. A seismic study delineates the position and local thickening of the Indian plate underlying the Himalayas and southern Tibet. We studied the formation of the Himalayan mountain range and the Tibetan Plateau by investigating their lithospheric structure. Using an 800-kilometer-long, densely spaced seismic array, we have constructed an image of the crust and upper mantle beneath the Himalayas and the southern Tibetan Plateau. The image reveals in a continuous fashion the Main Himalayan thrust fault as it extends from a shallow depth under Nepal to the mid-crust under southern Tibet. Indian crust can be traced to 31°N. The crust/mantle interface beneath Tibet is anisotropic, indicating shearing during its formation. The dipping mantle fabric suggests that the Indian mantle is subducting in a diffuse fashion along several evolving subparallel structures.

Moment-tensor inversion for offshore earthquakes east of Taiwan and their implications to regional collision

Reliable determination of source parameters for offshore earthquakes east of Taiwan with mb<5.5 was a difficult task because of the poor azimuthal coverage by local network and the lack of signals at teleseismic distances. We take advantage of the recently established “Broadband Array in Taiwan for Seismology” (BATS) to invert seismic moment tensors for 7 such events occurred in 1996. To cope with different patterns of background noise and unknown structural details, we utilize variable frequency bands in the inversion and adapt a two-step procedure to select best velocity models for individual epicenter-station paths. Our results are consistent with the overall patterns of regional collision and indicate that the resulting compressive stress has caused significant intraplate deformation within the Philippine Sea plate. Simulation of the regions geological evolution and orogenic processes should take this factor into account and allow the Philippine Sea plate to deform internally.

Effects of Topography on Seismic-Wave Propagation: An Example from Northern Taiwan

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.

Effects of Realistic Surface Topography on Seismic Ground Motion in the Yangminshan Region of Taiwan Based Upon the Spectral-Element Method and LiDAR DTM

We combine light detection and ranging (LiDAR) digital terrain model (DTM) data and an improved mesh implementation to investigate the effects of high- resolution surface topography on seismic ground motion based upon the spectral- element method. In general, topography increases the amplitude of shaking at mountain tops and ridges, whereas valleys usually have reduced ground motion, as has been observed in both records from past earthquakes and numerical simulations. However, the effects of realistic topography on ground motion have not often been clearly characterized in numerical simulations, especially the seismic response of the true ground surface. Here, we use LiDAR DTM data, which provide two-meter reso- lution at the free surface, and a spectral-element method to simulate three-dimensional (3D) seismic-wave propagation in the Yangminshan region in Taiwan, incorporating the effects of realistic topography. A smoothed topographic map is employed beneath the model surface in order to decrease mesh distortions due to steep ground surfaces. Numerical simulations show that seismic shaking in mountainous areas is strongly affected by topography and source frequency content. The amplification of ground motion mainly occurs at the tops of hills and ridges whilst the valleys and flat-topped hills experience lower levels of ground shaking. Interaction between small-scale to- pographic features and high-frequency surface waves can produce unusually strong shaking. We demonstrate that topographic variations can change peak ground accel- eration (PGA) values by 50% in mountainous areas, and the relative change in PGA between a valley and a ridge can be as high as a factor of 2 compared to a flat surface response. This suggests that high-resolution, realistic topographic features should be taken into account in seismic hazard analysis, especially for densely populated moun- tainous areas.

Three-dimensional simulations of seismic-wave propagation in the Taipei basin with realistic topography based upon the spectral-element method

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.

Review: Progress in Rotational Ground-Motion Observations from Explosions and Local Earthquakes in Taiwan

Rotational motions generated by large earthquakes in the far field have been successfully measured, and observations agree well with the classical elasticity theory. However, recent rotational measurements in the near field of earthquakes in Japan and in Taiwan indicate that rotational ground motions are 10 to 100 times larger than expected from the classical elasticity theory. The near-field strong-motion records of the 1999 Mw 7:6 Chi-Chi, Taiwan, earthquake suggest that the ground motions along the 100 km rupture are complex. Some rather arbitrary baseline corrections are necessary in order to obtain reasonable displacement values from double integra- tion of the acceleration data. Because rotational motions can contaminate acceleration observations due to the induced perturbation of the Earths gravitational field, we started a modest program to observe rotational ground motions in Taiwan. Three papers have reported the rotational observations in Taiwan: (1) at the HGSD station (Liu et al., 2009), (2) at the N3 site from two TAiwan Integrated GEodynamics Research (TAIGER) explosions (Lin et al., 2009), and (3) at the Taiwan campus of the National Chung-Cheng University (NCCU )( Wuet al., 2009). In addition, Langston et al. (2009) reported the results of analyzing the TAIGER explosion data. As noted by several authors before, we found a linear relationship between peak rotational rate (PRR in mrad=sec) and peak ground acceleration (PGA in m=sec 2 ) from local earthquakes in Taiwan, PRR 0:002 1:301 PGA, with a correlation coefficient of 0.988.

1999 Chi-Chi Earthquake: A Case Study on the Role of Thrust-Ramp Structures for Generating Earthquakes

The 21 September 1999 Chi-Chi earthquake ( M w 7.6) occurred on east-dipping shallow thrust faults that produced a high-relief surface rupture. Extraordinary surface breaks appeared that could be clearly traced for about 100 km across many counties. These thrust faults, the Chelungpu and Shihkang, are part of an active fold-and-thrust belt related to ongoing recent arc-continent collision. Measurement of slip vectors along the earthquake rupture indicates that the orientation of the maximum shear stress changed from a westward direction (N70-90°W) on the Chelungpu fault to a northwestward direction (N30-40°W) on the Shihkang fault. The stress field underwent a clockwise rotation of about 40° during the Chi-Chi earthquake. Near-rupture vertical displacements in the hanging wall of the Shihkang fault have more cumulative displacement than on the Chelungpu fault, which is consistent with Global Positioning System (GPS) measurements. Maximum vertical offset on the rupture was found to be about 10 m by the surficial rupture and GPS measurements. In addition, analysis of crustal deformation by GPS measurements on the hanging wall defines a coseismic uplift related to a fault ramp structure. Our synthesis of geological and geodetic analyses shows the importance of ramp structures associated with thrust faults for generating large earthquakes and provides a general framework for understanding earthquake in fold-and-thrust belts. Large surficial coseismic uplift and strong asperities appear to be a function of fault ramp geometry. Our analysis also indicates that, in general, ramp structures in fold-and-thrust belts may have a high potential in generating large earthquakes. Manuscript received 14 November 2000.

Three-dimensional elastic wave velocity structure of the Hualien region of Taiwan: Evidence of active crustal exhumation

The Hualien region of Taiwan is located at a complex transition of the boundary between the Eurasian and Philippine Sea Plates. To the southwest, the mountains of Taiwan are uplifting rapidly as a consequence of an ongoing arc-continent collision, while to the east the oceanic Philippine Sea Plate is subducting northward beneath Eurasia. We investigated the structure and dynamics of this region by analyzing seismograms of local earthquakes recorded during a deployment of the Portable Array for Numerical Data Acquisition II network. P and S wave velocity structures deduced from travel time tomography analysis show that the collisional suture to the south of Hualien is characterized by a narrow (< 10 km width), near vertically dipping zone of low velocities that extends to depths in excess of 20 km. Velocities in the Eastern Central Range west of the suture zone are significantly higher and define a feature 10–15 km wide that appears to be continuous from the near surface to depths as great as 40 km. Farther to the west beneath the Western Central Range, the velocities again decrease. Focal mechanisms of local earthquakes show that while thrust faulting is the predominate mode of deformation throughout the region, normal faulting occurs as well beneath the Eastern Central Range. Thus the rapid uplift of the mountains of Taiwan may be a result not only of compressional shortening but also of an excess of positive buoyancy. We suggest that the higher velocities and extensional mechanisms in the Eastern Central Range are caused by the ongoing exhumation of previously subducted continental crust, while the lower velocities to the west reflect continued underthrusting of the crust beneath the Eastern Central Range.

Variations of P wave speeds in the mantle transition zone beneath the northern Philippine Sea

Using waveforms and travel times from deep earthquakes, we constructed 16 seismic profiles, each of which constrains the radial variation in Vp over a small area beneath the northern Philippine Sea. Taken together, the azimuthal coverage of these profiles also places tight bounds on the lateral extent of a region of anomalously high Vp (up to 3% faster than average Earth models) originally suggested by travel time tomography. Unlike travel time tomography, which relies heavily on arrival times of the direct P phase, we utilize the waveforms and move-out of later arrivals that mainly sample the mantle transition zone of interest. Our results identify three important characteristics of the northern Philippine Sea anomaly that are distinct from previous results. First, being approximately a subhorizontal, laterally uniform feature, the anomaly is localized beneath the northwestern corner of the Philippine Sea, within a region of approximately 500×500 km2 immediately east of the Ryukyu arc. Second, the anomaly is well constrained to occur in the lower portion of the transition zone, extending all the way down to the 660-km discontinuity. Third, the presence of such a distinct anomaly reduces the contrast in Vp across the 660-km discontinuity from approximately 6% to 3%. Such a configuration is consistent with the interpretation that the anomaly is caused by a remnant of subducted slab, as negative buoyancy should rest the slab just above the 660-km discontinuity where resistance to subduction is expected from a negative Clapeyron slope during the spinel—Mg-Fe-perovskite transition.

An Observation of Rupture Pulses of the 20 September 1999 Chi-Chi, Taiwan, Earthquake from Near-Field Seismograms

The ground-velocity recordings of the 20 September 1999, Chi-Chi, Taiwan earthquake recorded at stations near the ruptured fault trace show a simple, large-amplitude, and long-period pulse following the S wave, which is closely associated with the surface faulting and the rupture process of thrust faulting. The conspicuous pulse on the ground-velocity seismogram following the S -wave arrival, called the S 1 phase, is interpreted as the superposition of the rupture pulses that nucleate at an asperity near and underneath the station and propagate up-dip and laterally along the fault plane toward the surface stations. The arrival times of the S 1 phase and the onsets of the permanent displacement at stations near and along the ruptured fault trace increase with hypocentral distance, suggesting that the rupture of the Chi-Chi earthquake might have initiated at the hypocenter of the mainshock and propagated both upward and laterally from south to north. On the basis of the travel-time differences between the S 1 phase and the direct S wave at the stations near and along the ruptured fault trace, the rupture velocities varied from 2.28 to 2.69 km/sec, with an average rupture velocity of about 2.49 km/sec. The rupture velocities decreased from south to north.