Estimation of the 2001 Kunlun Earthquake Fault Slip from GPS Coseismic Data Using Hori’S Inverse Method∗

Total Page:16

File Type:pdf, Size:1020Kb

Estimation of the 2001 Kunlun Earthquake Fault Slip from GPS Coseismic Data Using Hori’S Inverse Method∗ Earthq Sci (2009)22: 609−614 609 Doi: 10.1007/s11589-009-0609-x Estimation of the 2001 Kunlun earthquake fault slip from GPS coseismic data using Hori’s ∗ inverse method Honglin Jin and Hui Wang Institute of Earthquake Science, China Earthquake Administration, Beijing 100036, China Abstract The Hori’s inverse method based on spectral decomposition was applied to estimate coseismic slip distribution on the rupture plane of the 14 November 2001 MS8.1 Kunlun earthquake based on GPS survey results. The inversion result shows that the six sliding models can be constrained by the coseismic GPS data. The established slips mainly concentrated along the eastern segment of the fault rupture, and the maximum magnitude is about 7 m. Slip on the eastern segment of the fault rupture represents as purely left-lateral strike-slip. Slip on the western segment of the seismic rupture represents as mainly dip-slip with the maximum dip-slip about 1 m. Total predicted scalar seismic moment is 5.196×1020 N⋅m. Our results constrained by geodetic data are consistent with seismological results. Key words: Kunlun MS8.1 earthquake; coseismic GPS data; fault slip; inversion CLC number: P315.72+5 Document code: A placements at more than 40 sites along the fault rupture, 1 Introduction the epicenter, the distribution of surface rupture and co- seismic GPS data, etc. The research fields include geo- The 14 November 2001 MS8.1 Kunlun earthquake occurred in the west of Kunlun mountain pass in Qing- logical investigation, seismometry, level survey and hai province. It is the largest earthquake that occurred on GPS survey, etc. Field investigation found that the the east Kunlun fault. The surface rupture produced by maximum surface left-lateral coseismic dislocation is − the event was the longest one among historical earth- 6 7 m, and the maximum vertical displacement reaches quakes. It spans from west of the Hoh Sai Lake to east 4 m (Xu et al, 2002; Chen et al, 2003). However, field of the Kunlun mountain pass. The western end of the investigation just provided data on several sites along fault rupture is located near the Buka Daban Peak, and the fault rupture, not giving deep slip on the rupture the eastern end of the fault rupture is located 70 km east plane. Wan et al (2004, 2008) inverted distribution of of the Kunlun mountain pass. The total length of the deep slips on fault plane basing on geodetic data by us- fault rupture is about 426 km. The maximum horizontal ing Okada’s method. In Okada’s inverse method, dis- dislocation along the fault rupture is more than 6 m. The crete observation data were directly used to estimate slip maximum vertical dislocation is less than 1 m (Xu et al, distribution with some prior information. As a result, the data conditions, such as data density, distribution and 2002). The 1973 M7.3 Mayi earthquake, the 1997 MS7.5 Mayi earthquake, the 1937 M7.5 Huashixia earthquake, accuracy, directly affect inversion results. So, their re- sults are rough and cannot constrain detailed fault slip. and the 1963 MS7.1 Dulan earthquake occurred along the east Kunlun fault. The mechanisms of these earth- The representation theorem in elastostatics relates quakes are consistent with local deformation pattern surface displacements to slip distribution on a fault sur- (Shen et al, 2003). face through a linear integral equation that uses Green’s Previous studies have provided plenty of data on function (Maruyama, 1964; Okada, 1985). Therefore we can set up an inverse problem for reconstructing the the MS8.1 Kunlun earthquake. The data include dis- static image of a seismic or aseismic source from ob- ∗ served geodetic data. The Hori’s inverse method that Received 20 July 2009; accepted in revised form 30 October 2009; published 10 December 2009. based on a spectral decomposition of the Green’s func- Corresponding author. e-mail: [email protected] 610 Earthq Sci (2009)22: 609−614 tion, is similar to the SVD (singular-value decomposi- ated with k eigenvalues. tion) analysis method (Menke, 1984; Parker, 1977, 1994; ∞ = αϕ α + Kirsch, 1996). As a continuous inverse method, the Ho- u(x) ∑u (x) noise (3) α =1 ri’s method separates the inverse into two steps: firstly ∞ α u α calculate an inverse operator based on given domains of p( y) = ψ ( y) (4) ∑ α surface and fault structure; and secondly estimate the α=1 λ α α slip distribution from observed data. Moreover, the me- In above equations, u are coefficients of ϕ (x), and k is thod does not require any prior information about the measurable number of the eigenvalues (or modes). The smoothness of fault slip distribution since it assumes second term in equation (3) contains modes given by that the response function and source function are in- α ϕ (x) for α > k, but they cannot be distinguished from trinsically smooth. errors or noises due to the limited accuracy. Therefore, In order to get detailed slip distribution function on even though a response function cannot be measured, it fault plane, we will apply Hori’s method to invert the is possible to accurately estimate a continuous function coseismic slip distribution using observed GPS data. from a set of discrete data by using a set of eigenfunc- First, we introduce the inverse analysis method, then, α tions. The coefficients {u } are estimated by numerical invert the coseismic slip distribution of 2001 Kunlun computation, generally using the least squares method. earthquake using coseismic GPS data. Finally, we dis- Detailed interpretation for the formulation and the pa- cuss the tectonic implication of our results. rameters used refer to Hori (2001) and Jin et al (2007). 2 Inversion method based on a spectral 3 Estimation of back-slip distribution decomposition of Green’s function 3.1 Surface and fault structure We assume the medium to be an elastic half-space. First, we determine the domains of the surface and Let S be the surface plane and F be the fault plane, and the fault structure. Fault structure is determined from denote the points on S and F by vectors x and y, respec- geological results. Referring to previous researches tively. According to the representation theorem in elas- (Chen et al, 2003; Tan et al, 2004; Wan et al, 2004; Xu tostatics (Maruyama, 1964; Yabuki and Matsu’ura, and Chen, 2004), we take the depth range of the source 1992), the static surface displacements u(x) are ex- domain as 0.5−25 km to fully cover the rupture fault pressed as region. The fault rupture length and dip were taken as u(x) = g(x, y) p( y)dy (1) 390 km and 86°, respectively. ∫F The surface S and fault F are decomposed into N in terms of source function p( y) (slip or back slip) and and M segments. The convergence is checked by in- Green’s function. According to function analysis theory, creasing N and M which correspond to the discretization the Green’s function g(x, y) is regarded as a compact for surface and fault. It is important to examine the nu- operator that maps from source to response, and the op- merical validity of spectral decomposition. We note that erator admits the spectral decomposition (Hori, 2001) Green’s function must be discretized using the (N, M) ∞ points by employing a set of smooth functions when α α α g(x, y) = ∑λ ϕ (x)ψ ( y), (2) spectral decomposition is obtained numerically. α =1 A smoothed fault is plotted on a grid with 2 106 α α α where λ is the α-th eigenvalue, ϕ (x) and ψ ( y) are the elements for numerical computations. In this study, we associated eigenfunctions of the surface and the fault, use the point dislocations embedded in the smoothed respectively. The α-th eigenvalue is positive value that curved surface, instead of multi-rectangular fault planes α α converges to zero as α increases, and ϕ (x) and ψ ( y) as used in the previous studies. form orthonormal basis functions on S and F respec- The Earth’s surface is represented by a curved tively. plane of 870 km×610 km along the rupture fault and the u(x) is expressed as a linear combination of the surface was gridded into 3 315 elements for numerical 1 2 k measurable modes {ϕ (x), ϕ (x), ⋅⋅⋅, ϕ (x)}, and p( y) is computations. With the selected set of discrete points expressed as the linear combination of source eigen- representing the surface and fault domains (Figure 1), 1 2 k functions {ψ ( y), ψ ( y), ⋅⋅⋅, ψ ( y)}, which are associ- we calculate the Green’s function, g(x, y), and separate Earthq Sci (2009)22: 609−614 611 it into a set of surface and fault eigenfunctions by spec- α − α ⋅⋅⋅ changes ( ui+1 ui ) with different i (i =1, 2, , 25). The tral decomposition. Note that the surface deformation smaller coefficient of variation tends to more stable. caused by fault coseismic slip can be expressed by a Figure 3 shows stability of coefficients u1 and u2, from linear combination of surface eigenfunctions, and then which we can found the suitable k value with coefficient slip distribution is expressed by a linear combination of uk is 6. If k >6, large error resulted from the instable co- fault eigenfunctions. efficients will be introduced to inversion result. Once k=6 is determined, we can estimate the surface dis- placements function, as shown in Figure 2 (blue). It means the available GPS data can extract slip modes up to 6 for the given fault structure.
Recommended publications
  • Rupture Geometry and Multi-Segment Rupture of the November 2001 Earthquake in the Kunlun Fault System, Northern Tibet, China
    活断層・古地震研究報告,No. 4, p. 243-264, 2004 Rupture geometry and multi-segment rupture of the November 2001 earthquake in the Kunlun fault system, northern Tibet, China Bihong Fu1, Yasuo Awata2, Jianguo Du3 and Wengui He4 1, 2Active Fault Research Center, GSJ/AIST ([email protected], [email protected]) 3Institute of Earthquake Science, CEA ([email protected]) 4Lanzhou Institute of Seismology, CEA ([email protected]) Abstract: We document the spatial distribution and geometry of surface rupture zone produced by the 2001 Mw 7.8 Kunlun earthquake, based on high-resolution satellite images combined with the field measurements. Our results show that the 2001 surface rupture zone can be divided into five segments according to the geometry of surface rupture. They are the Sun Lake, Buka Daban-Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments from west to east. These segments, varying from 55 to 130 km in length, are separated by step-overs or bends. The Sun Lake segment extends about 65 km with a strike of N45º-75ºW (between 90º05'E-90º50'E) along the previously unrecognized West Sun Lake fault. A gap about 30 km long exists between the Sun Lake and Buka Daban Peak where no obvious surface ruptures can be observed either from the satellite images or the field observations. The Buka Daban-Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments run about 365 km striking N75º-85ºW along the southern slope of the Kunlun Mountains (between 91º07'E-94º58'E). This segmentation of the 2001 surface rupture zone is well correlated with the pattern of slip distribution measured in the field; that is, the abrupt changes of slip distribution occur at segment boundaries.
    [Show full text]
  • The Mw 7.8, 2001 Kunlunshan Earthquake: Extreme Rupture Speed Variability and Effect of Fault Geometry D
    JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111, B08303, doi:10.1029/2005JB004137, 2006 Click Here for Full Article The Mw 7.8, 2001 Kunlunshan earthquake: Extreme rupture speed variability and effect of fault geometry D. P. Robinson,1 C. Brough,1 and S. Das1 Received 2 November 2005; revised 19 January 2006; accepted 12 April 2006; published 11 August 2006. [1] By analyzing body wave seismograms, we show that the rupture speed on the Main Kunlun Fault during the Mw 7.8 2001 Kunlunshan, Tibet, earthquake was highly variable and the rupture process consisted of three stages. In the first stage, the rupture accelerated from rest to an average speed of 3.3 km/s over a distance of 120 km. The rupture then propagated for another 150 km at an apparent rupture speed exceeding the P wave speed. In the final stage, the earthquake fault bifurcates, and the rupture front slowed down. The region of highest rupture velocity is found to coincide with the region of highest fault slip, has the longest slip duration, and is where off-fault ground cracking is observed in the field. Stress drops are found to be higher in regions of higher rupture speeds. The greatest concentration of aftershocks is located near the fault bifurcation zone and hence coincides with the region of highest fault slip, highest stress drop and highest rupture velocity. The fault width is no more than 10 km in most places and is about 20 km in the region of highest slip. This narrow fault width is attributed to the fact that crust below this depth is sufficiently warm not to permit brittle failure to occur.
    [Show full text]
  • Systematical Stream Offsets Resulting from Large
    The Open Geology Journal, 2008, 2, 1-9 1 Systematical Stream Offsets Resulting from Large Earthquakes Along the Strike-Slip Kunlun Fault, Northern Tibetan Plateau: Evidence from the 2001 Mw 7.8 Kunlun Earthquake Aiming Lin* Graduate School of Science and Technology, Shizuoka University, Shizuoka 422-8529, Japan Abstract: Field investigations and interpretations of high-resolution remote sensing images reveal geomorphic features of cumulative offsets and deflections of Holocene streams and gullies along the western 450-km segment of the Kunlun fault which triggered the 2001 Mw 7.8 Kunlun earthquake in the northern Tibetan Plateau. The streams and gullies developed on Holocene alluvial fans are sinistrally displaced by up to 115 m, with 3-5 m offsets caused by the 2001 Kunlun earth- quake. Radiocarbon ages show that the alluvial fans formed in the past 6,000-9,000 years, suggesting an averaged slip rate of 16±3 mm/yr over the past 6,000-9,000 years on the western 450-km segment of the Kunlun fault. Geomorphic and geologic evidence confirms that the systematical offsets of streams and gullies are the results of repeated large earth- quakes and these topographic features are reliable indicators of seismic displacements accumulated on active strike-slip faults. Keywords: Stream offset, slip rate, strike-slip Kunlun fault, 2001 Mw7.8 Kunlun earthquake, Tibetan Plateau. INTRODUCTION images including 1m-resolution IKONOS and 0.61m-reso- lution QuickBird images and field investigations to examine During the past two decades, the integration of geologic, systematic offset features of streams, gullies and terrace ris- geomorphic, seismic, and geophysical information has led to ers related to the 2001 co-seismic surface ruptures developed increased recognition and understanding of the tectonic sig- upon Holocene alluvial fans along the western segment of nificance of geomorphic features caused by strike-slip along the Kunlun fault.
    [Show full text]
  • RESEARCH Slip Rate and Recurrence Intervals of The
    RESEARCH Slip rate and recurrence intervals of the east Lenglongling fault constrained by morphotectonics: Tectonic implications for the northeastern Tibetan Plateau Wenliang Jiang1, 2,*, Zhujun Han2, Peng Guo2, Jingfa Zhang1, Qisong Jiao1, Shuai Kang1, and Yunfeng Tian1 1KEY LABORATORY OF CRUSTAL DYNAMICS, INSTITUTE OF CRUSTAL DYNAMICS, CHINA EARTHQUAKE ADMINISTRATION (CEA), BEIJING 100085, CHINA 2INSTITUTE OF GEOLOGY, CHINA EARTHQUAKE ADMINISTRATION (CEA), BEIJING 100029, CHINA ABSTRACT The Lenglongling fault located in the northeast margin of the Tibetan Plateau plays an important role in accommodating the tectonic defor- mation of the Tibetan Plateau relative to the Gobi–Ala Shan platform to the north and the North China craton to the east. However, little is known about the fault due to a lack of previous research. In this study we use terrestrial light detection and ranging (LiDAR) data combined with high-resolution remote sensing images to survey offset landforms in the east part of the Lenglongling fault. Microtopographic analysis of well-preserved offset terraces, gullies, ridges, and pluvial fans in the highland environment allows evaluation of single-event slip and multievent cumulative slip. Our study provides an important assessment of the horizontal offset associated with the latest earthquake and four paleoearthquakes that were identified from a series of offset bedrock terraces by constructing a morphotectonic evolution model. Terrestrial LiDAR data indicate that the east Lenglongling fault follows a characteristic slip model. The single-event slip of this section is ~9.4 m; 7–8 paleoearthquakes are thought to have occurred during the Holocene, and a left-lateral strike-slip rate of 6.6 ± 0.3 mm/yr is estimated.
    [Show full text]
  • China National Report on Geodesy
    2003-2006 China National Report on Geodesy For The 24rd General Assembly of IUGG Perugia, Italy, 2-13 July, 2007 By Chinese National Committee for The International Union of Geodesy and Geophysics Beijing, China June, 2007 PREFACE The Chinese National Committee of IAG is pleased to present the 2003-2006 quadrennial China National Report on Geodesy to the Chinese National Committee for IUGG. During the last four years, significant advances have been made in the study of Geodesy in China. The presentation of the 2003-2006 China National Report on Geodesy is a reflection to these advances, and provides the record of Chinese contributions to geodesy. The report includes the following contents: (1)Crustal movement and astro-geodynamics research. (2)Gravity measurement and partial gravity field refining in China. (4)Geocentric coordinate framework maintenance. (5)Ocean geodesy in China. (6)Advances made in data processing in China. (7)The height determination of Qomolangma Feng (Mt. Everest). The current report is published only in CD. It is hoped that this National Report would be of help for Chinese scientists in exchanging the results and ideas in the research and application of geodesy with colleagues all over the world. Prof. CHENG Pengfei Chairman of Chinese National Committee of IAG Prof. DANG Yamin Secretary of Chinese National Committee of IAG China National Report on Geodesy (2003-2006) Report No.1 THE HEIGHT DETERMINATION OF QOMOLANGMA FENG (MT. EVEREST) IN 2005 CHEN Junyong, ZHANG Yanping, YUAN Janli, GUO Chunxi and ZHANG Peng (State Bureau of Surveying and Mapping, Beijing 100830, China) Qomolangma Feng – Mt.
    [Show full text]
  • (Ms = 8.1), Earthquake from Seismological and Geological
    Geophys. J. Int. (2009) 177, 555–570 doi: 10.1111/j.1365-246X.2008.04063.x Validation of the rupture properties of the 2001 Kunlun, China (M s = 8.1), earthquake from seismological and geological observations ∗ Yi-Ying Wen,1 Kuo-Fong Ma,1 Teh-Ru Alex Song,2 and Walter D. Mooney3 1Graduate Institute of Geophysics, National Central University, Taiwan. E-mail: [email protected] 2Seismological Laboratory, California Institute of Technology, USA 3US Geological Survey, Menlo Park, CA, USA Accepted 2008 October 31. Received 2008 October 31; in original form 2008 April 30 SUMMARY We determine the finite-fault slip distribution of the 2001 Kunlun earthquake (M s = 8.1) by inverting teleseismic waveforms, as constrained by geological and remote sensing field observations. The spatial slip distribution along the 400-km-long fault was divided into five segments in accordance with geological observations. Forward modelling of regional surface waves was performed to estimate the variation of the speed of rupture propagation during faulting. For our modelling, the regional 1-D velocity structure was carefully constructed for each of six regional seismic stations using three events with magnitudes of 5.1–5.4 distributed along the ruptured portion of the Kunlun fault. Our result shows that the average rupture velocity is about 3.6 km s−1, consistent with teleseismic long period wave modelling. The initial rupture was almost purely strike-slip with a rupture velocity of 1.9 km s−1, increasing to 3.5 km s−1 in the second fault segment, and reaching a rupture velocity of about 6 km s−1 in the third segment and the fourth segment, where the maximum surface offset, with a broad fault zone, was observed.
    [Show full text]
  • Supershear Earthquakes
    Copyright © 2020 Harsha S. Bhat SELFPUBLISHED TUFTE-LATEX.GITHUB.IO/TUFTE-LATEX Licensed under the Apache License, Version 2.0 (the “License”); you may not use this file except in compliance with the License. You may obtain a copy of the License at http://www.apache.org/licenses/LICENSE-2.0. Unless required by applicable law or agreed to in writing, software distributed under the License is distributed on an “AS IS” BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the License for the specific language governing permissions and limitations under the License. FOREWORD This manuscript encloses all documents required to obtain a Habilitation à Diriger des Recherches from École Normale Supérieure. I have synthesised my research work on supershear earthquakes, listed below, that were con- ducted over the last 13 years in collaboration with students, postdocs and colleagues at various institutions around the world. 1. Jara, J., L. Bruhat, S. Antoine, K. Okubo, M. Y. Thomas, Y. Klinger, R. Jolivet, and H. S. Bhat (2020). “Signature of supershear transition seen in damage and aftershock pattern”. submitted to Nat. Comm. 2. Amlani, F., H. S. Bhat, W. J. F. Simons, A. Schubnel, C. Vigny, A. J. Rosakis, J. Efendi, A. Elbanna, and H. Z. Abidin (2020). “Supershear Tsunamis and insights from the Mw 7.5 Palu Earthquake”. to appear in Nat. Geosci. 3. Mello, M., H. S. Bhat, and A. J. Rosakis (2016). “Spatiotemporal properties of sub-Rayleigh and supershear rupture velocity fields : Theory and Experiments”. J. Mech. Phys. Solids. DOI: 10.1016/j. jmps.2016.02.031.
    [Show full text]
  • Sagaing Fault
    SOE THURA TUN & WATKINSON Sagaing Fault 1 SOE THURA TUN & WATKINSON, I. M. 2016. The Sagaing Fault. In: BARBER, A. J., RIDD, M. F., KHIN 2 ZAW & RANGIN, C. (eds.) Myanmar: Geology, Resources and Tectonics. Geological Society, London, 3 Memoir. 4 The Sagaing Fault 5 SOE THURA TUN1 & IAN M. WATKINSON2* 6 7 1Myanmar Earthquake Committee, Myanmar Engineering Society Building, 8 Hlaing Universities Campus, Yangon, Myanmar 9 2Department of Earth Sciences, 10 Royal Holloway University of London, Egham, Surrey TW20 0EX, United Kingdom 11 *Corresponding author (e-mail: [email protected]) 12 Text words 15,888 References number (words) 163 (4,831) Tables words 0 Figures number 15 13 Abbreviated title: Sagaing Fault 14 The Sagaing Fault is amongst the longest and most active strike-slip faults in the world (e.g. Molnar 15 & Dayem 2010; Robinson et al. 2010; Searle & Morley 2011). It accommodates more than half of the 16 right-lateral motion between Sundaland and India, within the diffuse plate boundary along the eastern 17 margin of India, which occupies much of Myanmar (e.g. Vigny et al. 2003; Nielsen et al. 2004; 18 Socquet et al. 2006) (Fig. 1). Acting as a ridge-subduction transform (Le Dain et al. 1984; Guzmán- 19 Speziale & Ni 1996; Yeats et al. 1997), the 1500 km long Sagaing Fault links major thrust systems in 20 the north, such as the Naga, Lohit and Main Central thrust zones near the eastern Himalayan syntaxis, 21 to the Andaman Sea spreading centre in the south (e.g. Win Swe 1970; Le Dain et al.
    [Show full text]
  • Tectonophysics Global Catalog of Earthquake Rupture Velocities
    Tectonophysics xxx (xxxx) xxx–xxx Contents lists available at ScienceDirect Tectonophysics journal homepage: www.elsevier.com/locate/tecto Global catalog of earthquake rupture velocities shows anticorrelation between stress drop and rupture velocity Agnès Chouneta,*, Martin Valléea, Mathieu Causseb, Françoise Courboulexc a Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, CNRS, Paris, France b Université Grenoble Alpes, IFSTTAR, CNRS, IRD, ISTerre, Grenoble 38000, France c Université Côte d’Azur, CNRS, Observatoire de la Côte d’Azur, IRD, Géoazur, Valbonne, France ARTICLE INFO ABSTRACT Keywords: Application of the SCARDEC method provides the apparent source time functions together with seismic moment, Rupture velocity depth, and focal mechanism, for most of the recent earthquakes with magnitude larger than 5.6–6. Using this Rupture propagation large dataset, we have developed a method to systematically invert for the rupture direction and average rupture Source time functions velocity V r, when unilateral rupture propagation dominates. The approach is applied to all the shallow Global earthquakes seismology (z < 120 km) earthquakes of the catalog over the 1992–2015 time period. After a careful validation process, Stress drop variability rupture properties for a catalog of 96 earthquakes are obtained. The subsequent analysis of this catalog provides Peak Ground Acceleration variability several insights about the seismic rupture process. We first report that up-dip ruptures are more abundant than down-dip ruptures for shallow subduction interface earthquakes, which can be understood as a consequence of the material contrast between the slab and the overriding crust. Rupture velocities, which are searched without any a-priori up to the maximal P wave velocity (6000–8000 m/s), are found between 1200 m/s and 4500 m/s.
    [Show full text]
  • Supershear Earthquake Ruptures – Theory, Methods, Laboratory Experiments and Fault Superhighways: an Update
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Springer - Publisher Connector Chapter 1 Supershear Earthquake Ruptures – Theory, Methods, Laboratory Experiments and Fault Superhighways: An Update Shamita Das Abstract The occurrence of earthquakes propagating at speeds not only exceeding the shear wave speed of the medium (~3 km/s in the Earth’s crust), but even reaching compressional wave speeds of nearly 6 km/s is now well established. In this paper, the history of development of ideas since the early 1970s is given first. The topic is then discussed from the point of view of theoretical modelling. A brief description of a method for analysing seismic waveform records to obtain earth- quake rupture speed information is given. Examples of earthquakes known to have propagated at supershear speed are listed. Laboratory experiments in which such speeds have been measured, both in rocks as well as on man-made materials, are discussed. Finally, faults worldwide which have the potential to propagate for long distances (> about 100 km) at supershear speeds are identified (“fault superhighways”). 1.1 Introduction Seismologists now know that one of the important parameters controlling earth- quake damage is the fault rupture speed, and changes in this rupture speed (Madariaga 1977, 1983). The changes in rupture speed generate high-frequency damaging waves Thus, the knowledge of how this rupture speed changes during earthquakes and its maximum possible value are essential for reliable earthquake hazard assessment. But how high this rupture speed can be has been understood only relatively recently. In the 1950–1960s, it was believed that earthquake ruptures could only reach the Rayleigh wave speed.
    [Show full text]
  • Millennial Recurrence of Large Earthquakes on the Haiyuan Fault
    Bulletin of the Seismological Society of America, Vol. 97, No. 1B, pp. 14–34, February 2007, doi: 10.1785/0120050118 E Millennial Recurrence of Large Earthquakes on the Haiyuan Fault near Songshan, Gansu Province, China by Jing Liu-Zeng, Yann Klinger, Xiwei Xu, Ce´cile Lasserre, Guihua Chen, Wenbing Chen, Paul Tapponnier, and Biao Zhang Abstract The Haiyuan fault is a major active left-lateral fault along the northeast edge of the Tibet-Qinghai Plateau. Studying this fault is important in understanding current deformation of the plateau and the mechanics of continental deformation in general. Previous studies have mostly focused on the slip rate of the fault. Paleo- seismic investigations on the fault are sparse, and have been targeted mostly at the stretch of the fault that ruptured in the 1920 M ϳ8.6 earthquake in Ningxia Province. To investigate the millennial seismic history of the western Haiyuan fault, we opened two trenches in a small pull-apart basin near Songshan, in Gansu Province. The excavation exposes sedimentary layers of alternating colors: dark brown silty to clayey deposit and light yellowish brown layers of coarser-grained sandy deposit. The main fault zone is readily recognizable by the disruption and tilting of the layers. Six paleoseismic events are identified and named SS1 through SS6, from youngest to oldest. Charcoal is abundant, yet generally tiny in the shallowest parts of the trench exposures. Thirteen samples were dated to constrain the ages of paleoseismic events. All six events have occurred during the past 3500–3900 years. The horizontal offsets associated with these events are poorly known.
    [Show full text]
  • Shallow Seismic Structure of Kunlun Fault Zone in Northern Tibetan Plateau, China: Implications for the 2001 Ms8.1 Kunlun Earthquake
    Geophys. J. Int. (2009) 177, 978–1000 doi: 10.1111/j.1365-246X.2009.04049.x Shallow seismic structure of Kunlun fault zone in northern Tibetan Plateau, China: implications for the 2001 Ms8.1 Kunlun earthquake Chun-Yong Wang,1 Walter D. Mooney,2 Zhifeng Ding,1 Jiansi Yang,1 Zhixiang Yao1 and Hai Lou1 1Institute of Geophysics, China Earthquake Administration, Beijing 100081, China E-mail: [email protected] 2US Geological Survey, Menlo Park, CA 94025, USA Accepted 2008 November 12. Received 2008 November 11; in original form 2007 October 10 SUMMARY The shallow seismic velocity structure of the Kunlun fault zone (KLFZ) was jointly deduced from seismic refraction profiling and the records of trapped waves that were excited by five explosions. The data were collected after the 2001 Kunlun Ms8.1 earthquake in the northern Tibetan Plateau. Seismic phases for the in-line record sections (26 records up to a distance of 15 km) along the fault zone were analysed, and 1-D P- and S-wave velocity models of shallow crust within the fault zone were determined by using the seismic refraction method. Sixteen seismic stations were deployed along the off-line profile perpendicular to the fault zone. Fault-zone trapped waves appear clearly on the record sections, which were simulated with a 3-D finite difference algorithm. Quantitative analysis of the correlation coefficients of the synthetic and observed trapped waveforms indicates that the Kunlun fault-zone width is 300 m, and S-wave quality factor Q within the fault zone is 15. Significantly, S-wave velocities within the fault zone are reduced by 30–45 per cent from surrounding rocks to a depth of at least 1–2 km, while P-wave velocities are reduced by 7–20 per cent.
    [Show full text]