Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

Documentation of soil conditions at liquefaction and non-liquefaction sites from 1999 Chi–Chi (Taiwan) earthquake

The 1999 Chi–Chi, Taiwan, earthquake provides case histories of ground failure and non-ground failure that are valuable to the ongoing development of liquefaction susceptibility, triggering and surface manifestation models because the data occupy sparsely populated parameter spaces (i.e. high cyclic stress ratio and high fines content with low to moderate soil plasticity). In this paper, we document results from several large site investigation programs conducted in Nantou, Wufeng and Yuanlin, Taiwan. The seismic performance of the investigated sites include non-ground failure building and free-field sites, building sites with partial foundation bearing failures, free-field lateral spread sites, and free-field level ground sites with sediment boils. Field and laboratory investigation protocols for the sites are described, including cone penetration testing (some with pore pressure and shear wave velocity measurements) and rotary wash borings with standard penetration testing (including energy measurements). Implications of the SPT energy measurements with respect to established guidelines for the estimation of SPT energy ratio (including short rod corrections) are presented. Finally, data for three example sites are shown that illustrate potential applications of the data set, and which also demonstrate a condition where existing liquefaction analysis procedures fail to predict the observed field performance.     

Evaluation of hard rock spectral models for the Taiwan region on the basis of the 1999 Chi–Chi earthquake data

We analyze the ability of different spectral models to describe the frequency content of ground motion during the 1999 Chi–Chi earthquake (M_W=7.6, Taiwan) and two large (M_L=6.8) aftershocks. The spectral models evaluated include the one-corner model of Brune applied with various key parameters (seismic moment and stress drop), and the two-corner-frequency models proposed for eastern North America [Bull. Seismol. Soc. Am. 83 (1993) 1778] and California [Bull. Seismol. Soc. Am. 90 (2000) 255]. The ground-motion spectra predicted by these spectral models for hypothetical very hard rock site were compared with the Chi–Chi earthquake data obtained on rock (class B) and soft rock or very dense soil (class C) sites. The approach also allows us evaluating the generalized empirical amplification function for class B and C sites in the region.

It has been found that, the amplitude spectra of recorded ground acceleration (the mainshock and aftershocks) for frequencies larger than 0.3–0.4 Hz agree with the modelled two-corner-frequency spectra calculated using the model proposed for California. The single-corner-frequency model also provides a good agreement with the observations when using so-called ‘short-period seismic moment’ [Phys. Earth Planet. Interiors 37 (1985) 108] instead of the reported values obtained from long-period waves. The key parameters used in the single-corner model coincide with parameters of subsources evaluated for the complicated mainshock source. Therefore, it is possible to confirm the suggestion that the short-period seismic waves, at least for the thrust earthquakes, are generated mainly from the fracture of small-scale heterogeneities. The use of two-corner-frequency source model for earthquake spectrum that is based on long-period seismic moment value is equivalent, for frequencies larger than 0.3–0.4 Hz, to the use of single-corner-frequency model that is based on the parameters of major subsource.     

Nonlinear identification of the soil response at Dahan downhole array site during the 1999 Chi–Chi earthquake

Records of the 1999 Chi–Chi earthquake, provided by the Dahan downhole array were analyzed, and nonlinear hysteretic stress–strain relations in the soil layers were estimated by the method previously used for studying the response of soils during the 1995 Kobe and the 2000 Tottori earthquakes. The obtained models of the soil behavior were applied for evaluating changes of the shear moduli in the soil layers and for the nonlinear identification of the soil behavior at Dahan site during the Chi–Chi earthquake. We found that reduction of the shear moduli in the soil layers did not exceed 5%, and the soil response was virtually linear. The content of the nonlinear components in the soil response was about 5% of the intensity of the response, and it was mostly due to the odd-order nonlinearities. A similarity in the stress–strain relations describing the behavior of soils during the 1995 Kobe earthquake, the 2000 Tottori earthquake, and the 1999 Chi–Chi earthquake was found, which indicates the possibility to describe the behavior of similar types of soils at different sites by similar stress–strain relations and predict soil behavior in future earthquakes.     

Investigating the impact of the Chi‐Chi earthquake on the occurrence of debris flows using artificial neural networks

Debris flows have caused enormous losses of property and human life in Taiwan during the last two decades. An efficient and reliable method for predicting the occurrence of debris flows is required. The major goal of this study is to explore the impact of the Chi‐Chi earthquake on the occurrence of debris flows by applying the artificial neural network (ANN) that takes both hydrological and geomorphologic influences into account. The Chen‐Yu‐Lan River watershed, which is located in central Taiwan, is chosen for evaluating the critical rainfall triggering debris flows. A total of 1151 data sets were collected for calibrating model parameters with two training strategies. Significant differences before and after the earthquake have been found: (1) The size of landslide area is proportioned to the occurrence of debris flows; (2) the amount of critical rainfall required for triggering debris flows has reduced significantly, about half of the original critical rainfall in the study case; and (3) the frequency of the occurrence of debris flows is largely increased. The overall accuracy of model prediction in testing phase has reached 96·5%; moreover, the accuracy of occurrence prediction is largely increased from 24 to 80% as the network trained with data from before the Chi‐Chi earthquake sets and with data from the lumped before and after the earthquake sets. The results demonstrated that the ANN is capable of learning the complex mechanism of debris flows and producing satisfactory predictions. Copyright © 2009 John Wiley & Sons, Ltd.     

Fault rock analysis of the northern part of the Chelungpu Fault and its relation to earthquake faulting of the 1999 Chi-Chi earthquake, Taiwan

Abstract   The 1999 Chi-Chi earthquake in Taiwan ( M _w = 7.6) produced a surface rupture along the north–south-striking Chelungpu thrust fault with pure dip-slip (east side up) and left lateral strike-slip displacements. Near-field strong-motion data for the northern part of the fault illustrate a distinct lack of the high-frequency seismic radiation associated with a large slip (10–15 m) and a rapid slip velocity (2–4 m/s), suggesting a smooth seismic slip associated with low dynamic frictional resistance on the fault. A drillhole was constructed at shallow depths in the possible fault zones of the northern part of the Chelungpu Fault, which may have slipped during the 1999 earthquake. One of the zones consists of a 20-cm-thick, unconsolidated fault breccia with a chaotic texture lacking both discrete slip surfaces (e.g. Riedel shears) and grain crushing. Other possible fault zones are marked by the narrow (less than a few centimeters) gouge zone in which clayey material intrudes into the damaged zone outside of the gouge zone. These characteristic fault rock textures suggest that the slip mechanisms at shallow levels during the earthquake involved either granular flow of initially unconsolidated material or slip localization under elevated pore pressure along the narrow clayey gouge zone. Because both mechanisms lead to low dynamic frictional resistance on the fault, the rapid seismic slip in the deep portions of the fault (i.e. the source region of strong-motion radiation) could have been accommodated by frictionless slip on the shallow portions of the fault. The combination of strong-motion data and fault rock analysis suggests that smooth slip associated with low dynamic friction occurred on both the deep and shallow portions of the fault, resulting in a large slip between the source region and the surface in the northern region.     

2007年6月3日宁洱6.4级地震地表变形的构造分析和解释

以2007年6月3日云南省宁洱6.4级地震地表裂缝、喷砂冒水、地震滑坡、地震崩塌等资料为基础,结合震区的地质构造、震源机制解、余震分布等资料,研究了本次地震的发震断层及其强震频繁发生的动力学机制. 地震地表裂缝资料表明,北西向的宁洱断裂班海段具有右旋走滑的特征,北北东向的断层具有左旋走滑的性质. 地震地质灾害集中分布在330deg;方向上的长13.5 km、宽4 km的范围内. 等震线的长轴方向亦为330deg;,Ⅷ度区长轴长度为17 km. 震源机制解资料显示,宁洱6.4级地震的北西向节面为右旋错动,北东向节面为左旋错动. 大于等于2.0级的余震分布优势方向为330deg;,深度为3——12 km,优势深度为3——10 km,余震分布与地震地质灾害集中分布带一致. 以上资料说明,宁洱断裂班海段为这次地震的发震断层. 最后以活动地块理论为基础,讨论了宁洱地震的动力来源为印度板块的向北推挤使青藏高原向东滑移,在滇西南地区转化为向南南东方向的挤出,使宁洱附近网络状的北西向断裂发生右旋活动,北东向断裂发生左旋活动. 这种构造格局可能是该区频繁发生6.0——6.9级地震的原因.      

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Compositional characteristics and geodynamic significance of late Miocene volcanic rocks associated with the Chah Zard epithermal gold–silver deposit,southwest Yazd,Iran

Whole‐rock geochemical and Sr–Nd isotopic data are presented for late Miocene volcanic rocks associated with the Chah Zard epithermal Au–Ag deposit in the Urumieh‐Dokhtar Magmatic Arc (UDMA), Iran, to investigate the magma source, petrogenesis and the geodynamic evolution of the study area. The Chah Zard andesitic to rhyolitic volcanic rocks are characterized by significant Large Ion Lithophile Element (LILE) and Light Rare Earth Element (LREE) enrichment coupled with High Field Strength Element (HFSE) depletion. Our geochemical data indicate an adakitic‐like signature for the volcanic rocks (e.g. SiO₂ > 62 wt%, Al₂O₃ > 15 wt%, MgO < 1.5 wt%, Sr/Y > 70, La/Yb > 35, Yb < 1 ppm, and Y < 18 ppm, and no significant Eu anomalies), distinguishing them from the other volcanic rocks of the UDMA. The Chah Zard volcanic rocks have similar Sr and Nd isotopic compositions; the ⁸⁷Sr/⁸⁶Sr(i) ratios range from 0.704 902 to 0.705 093 and the ε_Nd(i) values are from +2.33 to +2.70. However, the rhyolite porphyry represents the final stage of magmatism in the area and has a relatively high ⁸⁷Sr/⁸⁶Sr ratio (0.705 811). Our data suggest that the andesitic magmas are from a heterogeneous source and likely to result from partial melting of a metasomatized mantle wedge associated with a mixture of subducted oceanic crust and sediment. These melts subsequently underwent fractional crystallization along with minor amounts of crustal assimilation. Our study is consistent with the model that the volcanic host rocks to epithermal gold mineralization in the UDMA are genetically related to late Miocene Neo‐Tethyan slab break‐off beneath Central Iran.