Preliminary results pertaining to coseismic displacement and preseismic strain accumulation of the Lushan M S7.0 earthquake, as reflected by GPS surveying

This paper presents the coseismic displacement and preseismic deformation fields of the Lushan M _S7.0 earthquake that occurred on April 20, 2013. The results are based on GPS observations along the Longmenshan fault and within its vicinity. The coseismic displacement and preseismic GPS results indicate that in the strain release of this earthquake, the thrust rupture is dominant and the laevorotation movement is secondary. Furthermore, we infer that any possible the rupture does not reach the earth’s surface, and the seismogenic fault is most likely one fault to the east of the Guanxian-Anxian fault. Some detailed results are obtainable. (1) The southern segment of the Longmenshan fault is locked preceding the Lushan earthquake. After the Wenchuan earthquake, the strain accumulation rate in the southeast direction accelerates in the epicenter of the Lushan earthquake, and the angle between the principal compressional strain and the seismogenic fault indicates that a sinistral deformation background in the direction of the seismogenic fault precedes the Lushan earthquake. Therefore, it is evident that the Wenchuan M _S8.0 earthquake accelerated the pregnancy of the Lushan earthquake. (2) The coseismic displacements reflected by GPS data are mainly located in a region that is 230 km (NW direction) × 100 km (SW direction), and coseismic displacements larger than 10 mm lie predominantly in a 100-km region (NW direction). (3) On a large scale, the coseismic displacement shows thrust characteristics, but the associated values are remarkably small in the near field (within 70 km) of the earthquake fault. Meanwhile, the thrust movement in this 70-km region does not correspond with the attenuation characteristics of the strain release, indicating that the rupture of this earthquake does not reach the earth’s surface. (4) The laevorotation movements are remarkable in the 50-km region, which is located in the hanging wall that is close to the earthquake fault, and the corresponding values in this case correlate with the attenuation characteristics of the strain release.     

Stress changes from the 2008 Wenchuan earthquake and increased hazard in the Sichuan basin

On 12 May 2008, the devastating magnitude 7.9 (Wenchuan) earthquake struck the eastern edge of the Tibetan plateau, collapsing buildings and killing thousands in major cities aligned along the western Sichuan basin in China. After such a large-magnitude earthquake, rearrangement of stresses in the crust commonly leads to subsequent damaging earthquakes. The mainshock of the 12 May earthquake ruptured with as much as 9 m of slip along the boundary between the Longmen Shan and Sichuan basin, and demonstrated the complex strike-slip and thrust motion that characterizes the region. The Sichuan basin and surroundings are also crossed by other active strike-slip and thrust faults. Here we present calculations of the coseismic stress changes that resulted from the 12 May event using models of those faults, and show that many indicate significant stress increases. Rapid mapping of such stress changes can help to locate fault sections with relatively higher odds of producing large aftershocks.     

Interseismic and coseismic displacements of the Lushan Ms 7.0 earthquake inferred from leveling measurements

By using precise leveling data observed between 1985 and 2010 across the south section of the Longmenshan fault zone, and eliminating the coseismic displacements caused by the Wenchuan Ms 8.0 earthquake, the interseismic vertical deformation field was obtained. The result shows that the Lushan region, located between the Shuangshi-Dachuan fault （front range of the Long- menshan fault） and the Xinkaidian fault （south section of the Dayi fault）, is situated in the intersection zone of positive and negative vertical deformation gradient zones, indicating that this zone was locked within 25 years before the Lushan earthquake. Based on leveling data across the rupture zone surveyed between 2010 and 2013, and by eliminating the vertical deformation within 3 years before the earthquake, the coseismic vertical displacement was derived. The coseismic vertical displacement for the benchmark DD35, which is closest to the epicenter, is up to 198.4 mm （with respect to MY165A）. The coseismic dis- placement field revealed that the northwest region （hanging wall） moved upwards in comparison with the southeastern region （foot wall）, suggesting that the seismogenic fault mainly underwent thrust faulting. By comparing the coseismic and interseismic vertical deformation fields, it was found that the mechanisms of this earthquake are consistent with the elastic rebound theory; the elastic strain energy （displacement deficit） accumulated before the Lu- shan earthquake was released during this quake.     

The spatial distribution pattern of landslides triggered by the 20 April 2013 Lushan earthquake of China and its implication to identification of the seismogenic fault

After the 20 April 2013 Lushan MS6.6 earthquake occurred,investigation and identification of the seismogenic fault for this event have become a focused and debatable issue.This work prepared an initial landslide inventory map related to the Lushan earthquake based on field investigations and visual interpretation of high-resolution aerial photographs and provided evidence for solving the issue aforementioned.The analysis of three landslide-density profiles perpendicular to strike direction of the probable seismogenic fault shows that many landslides occurred on the footwall of the Shuangshi–Dachuan fault(SDF),without sudden change of landslide density near the fault.Very few landslides were detected near the Dayi fault(DF)and also no change of landslide density there.While obvious sudden change of landslide density appeared about 1–2 km from the northwest to the western Shangli fault(WSF),and the landslide density on the hanging wall of the fault is obviously higher than that of on the footwall.Therefore,we infer that the seismogenic fault for the Lushan earthquake is neither the SDF nor the DF,rather probably the WSF located between these two faults,which is an evident linear trace on the earth surface.Meanwhile,the coseismic slip did not propagate upward to the ground,implying the Lushan earthquake was spawned by a blind-thrust-fault beneath the WSF.     

Relocation of the mainshock and aftershock sequences of M S7.0 Sichuan Lushan earthquake

The mainshock of April 20, 2013 Sichuan Lushan M _S7.0 earthquake was relocated using a 3-D velocity model. Double difference algorithm was applied to relocate aftershock sequences of Lushan earthquake. The locations of 2405 aftershocks were determined. The location errors in E-W, N-S and U-D direction were 0.30, 0.29 and 0.59 km on average, respectively. The location of the mainshock is 102.983°E, 30.291°N and the focal depth is 17.6 km. The relocation results show that the aftershocks spread approximately 35 km in length and 16 km in width. The dominant distribution of the focal depth ranges from 10 to 20 km. A few earthquakes occurred in the shallow crust. Focal depth profiles show fault planes dip to the northwest, manifested itself as a listric thrust fault. The dip angle is steep in the shallow crust and gentle in the deep crust. Although the epicenters of aftershocks distributed mainly along both sides of the Shuangshi-Dachuan fault, the seismogenic fault may be a blind thrust fault on the eastern side of the Shuangshi-Dachuan fault. Earthquake relocation results reveal that there is a southeastward tilt aftershock belt intersecting with the seismogenic fault with y-shape. We speculate it is a back thrust fault that often appears in a thrust fault system. Lushan earthquake triggered the seismic activity of the back thrust fault.     

Slip model for the 2011 M_w 9.0 Sendai （Japan） earthquake and its M_w 7.9 aftershock derived from GPS data

Using GPS-measured coseismic and post-seismic displacements for the 8 h following the M w 9.0 Sendai earthquake of March 11, 2011, coseismic and post-seismic fault slip models were developed based on a layered crustal model. The geodetic moment magnitude of the main shock was measured as approximately M w 8.98. The slip exhibits clear reverse characteristics, with a maximum near the hypocenter, and a magnitude of about 23.3 m. Some strike-slip behavior may occur on the two sides of the peak rupture zone. Almost 90% of the seismic moments released by the main shock occurred at depths less than 40 km. The energy released by the fault slip in the 8 h following the main shock is approximately equal to an earthquake of M w 8.13. With a maximum of ~1.5 m, the post-seismic slip was concentrated in the southwestern part of the coseismic rupture fault, which agrees well with the location and behavior of the M w 7.9 aftershock. This implies that the post-seismic deformation in the 8 h after the main shock was mainly induced by the M w 7.9 aftershock. In addition, a post-seismic slip of 0.2-0.4 m was observed at the down-dip extension of the coseismic rupture, which may have been caused by the effect of after-slip during this period.     

The maximum coseismic vertical surface displacement and surface deformation pattern accompanying the Ms 8.0 Wenchuan earthquake

The amount of coseismic deformation and its distribution of the Wenchuan earthquake provide important scientific bases for revealing the mechanisms of earthquake preparation and characterizing the rupture propagation of the Wenchuan earthquake. The previous studies have indicated that the earthquake ruptured the middle-to-north segment of the Longmenshan central fault and the middle segment of the Longmenshan range-front fault, which are characterized by two surface rupture zones of 240 km and 90 km in length, respectively. Based on the pre-earthquake information and photos of landforms and buildings obtained through ge-ologic and geomorphic survey of the area around Shaba Village of Beichuan County, Sichuan Province and the extensive interview with local villagers, we measured the displacements of the major terrain features and the dislocated buildings by total station instruments and differential GPS and obtained the maximum vertical displacement of 9±0.5 m and right-lateral displacement of 2±0.5 m around the Zou’s house in Shaba Village. Though the near-surface deformation exhibits a normal faulting around Shaba Village, the dynamic environment has not changed on the whole. The NW wall of the fault uplifted but without gravity gliding as normally occurring on the hanging wall of a normal fault, which proves that the 9±0.5 m displacement should be the maximum coseismic vertical displacement of the May 12, 2008 Ms 8.0 Wenchuan earthquake.     

Source rupture process of Lushan M S7.0 earthquake, Sichuan, China and its tectonic implications

The source rupture process of the M _S7.0 Lushan earthquake was here evaluated using 40 long-period P waveforms with even azimuth coverage of stations. Results reveal that the rupture process of the Lushan M _S7.0 event to be simpler than that of the Wenchuan earthquake and also showed significant differences between the two rupture processes. The whole rupture process lasted 36 s and most of the moment was released within the first 13 s. The total released moment is 1.9×10¹⁹N m with M _W=6.8. Rupture propagated upwards and bilaterally to both sides from the initial point, resulting in a large slip region of 40 km×30 km, with the maximum slip of 1.8 m, located above the initial point. No surface displacement was estimated around the epicenter, but displacement was observed about 20 km NE and SW directions of the epicenter. Both showed slips of less than 40 cm. The rupture suddenly stopped at 20 km NE of the initial point. This was consistent with the aftershock activity. This phenomenon indicates the existence of significant variation of the medium or tectonic structure, which may prevent the propagation of the rupture and aftershock activity. The earthquake risk of the left segment of Qianshan fault is worthy of attention.     

The Lushan M S7.0 earthquake and activity of the southern segment of the Longmenshan fault zone

Following the Lushan M _S7.0 earthquake on 20 April 2013, a topic of much concern is whether events of M _S7 or greater could occur again on the southern segment of the Longmenshan fault zone. In providing evidence to answer this question, this work analyzes the tectonic relationship between the Lushan event and the 2008 Wenchuan earthquake and the rupture history of the southern segment of the Longmenshan fault zone, through field investigations of active tectonics and paleoearthquake research, and our preliminary conclusions are as follows. The activity of the southern segment of the Longmenshan fault zone is much different to that of its central section, and the late Quaternary activity has propagated forward to the basin in the east. The seismogenic structure of the 2008 Wenchuan earthquake is the central-fore-range fault system, whereas that of the 2013 Lushan event is attributed to the fore-range-range-front fault system, rather than the central fault. The southern segment of the Longmenshan fault zone becomes wider towards the south with an increasing number of secondary faults, of which the individual faults exhibit much weaker surface activity. Therefore, this section is not as capable of generating a major earthquake as is the central segment. It is most likely that the 2013 earthquake fills the seismic gap around Lushan on the southern segment of the Longmenshan fault zone.     

Rupture of the 2004 Sumatra-Andaman earthquake inferred from direct P-wave imaging

The Sumatra-Andaman earthquake on December 26, 2004 is the first well recorded gigantic earthquake (moment magnitude MW 9.3) by modern broadband seismic and Global Positioning System networks. The rich seismic and geodetic recordings have documented unprecedented details about the earthquake rupture, coseismic and postseismic deformations. This is a report of detailed images of the rupture process using the first-arriving compressional waves recorded by the China National Digital Seismic Network (CNDSN). An improved imaging condition was employed to account for the sparse distribution of the CNDSN stations. The resulting images are consistent with the major rupture features reported by previous seismic and geodetic studies. It is found that the earthquake rupture initiated at offshore of northwestern Sumatra and propagated in the north northwest direction at a speed of 2.7 ± 0.2 km/s. The rupture continued for at least 420 s and extended about 1200-1300 km along the Andaman trough with two bursts of seismic energy.     

Interaction between adjacent left-lateral strike-slip faults and thrust faults: the 1976 Songpan earthquake sequence

Based on the published focal mechanisms we have built the fault model of the main shocks of the 1976 Songpan earthquake sequence and calculated the coseismic Coulomb stress changes in the region. The results show that most of the aftershocks had occurred in the region where the Coulomb stresses had been increased, indicating a triggering relationship between the main shocks and the aftershocks. We also show that the first main shock （Ms = 7.2）, which is a left-lateral slip event, had increased the Coulomb stresses by 5×10^5 Pa at the second main shock （a thrust event with Ms = 6.7）. Therefore, we conclude that the first main shock had triggered the second main shock. The third main shock is also a left-lateral event, however, the triggering relationship between the third main shock and the previous two events is less obvious. General model calculations show that there is a good triggering relationship between adjacent left-lateral slip fault and thrust fault, but triggering between parallel slip faults is rather weak.     

Predicting the endpoints of earthquake ruptures

The active fault traces on which earthquakes occur are generally not continuous, and are commonly composed of segments that are separated by discontinuities that appear as steps in map-view. Stress concentrations resulting from slip at such discontinuities may slow or stop rupture propagation and hence play a controlling role in limiting the length of earthquake rupture. Here I examine the mapped surface rupture traces of 22 historical strike-slip earthquakes with rupture lengths ranging between 10 and 420 km. I show that about two-thirds of the endpoints of strike-slip earthquake ruptures are associated with fault steps or the termini of active fault traces, and that there exists a limiting dimension of fault step (3-4 km) above which earthquake ruptures do not propagate and below which rupture propagation ceases only about 40 per cent of the time. The results are of practical importance to seismic hazard analysis where effort is spent attempting to place limits on the probable length of future earthquakes on mapped active faults. Physical insight to the dynamics of the earthquake rupture process is further gained with the observation that the limiting dimension appears to be largely independent of the earthquake rupture length. It follows that the magnitude of stress changes and the volume affected by those stress changes at the driving edge of laterally propagating ruptures are largely similar and invariable during the rupture process regardless of the distance an event has propagated or will propagate.     

Analysing the 1811-1812 New Madrid earthquakes with recent instrumentally recorded aftershocks

Although dynamic stress changes associated with the passage of seismic waves are thought to trigger earthquakes at great distances, more than 60 per cent of all aftershocks appear to be triggered by static stress changes within two rupture lengths of a mainshock. The observed distribution of aftershocks may thus be used to infer details of mainshock rupture geometry. Aftershocks following large mid-continental earthquakes, where background stressing rates are low, are known to persist for centuries, and models based on rate-and-state friction laws provide theoretical support for this inference. Most past studies of the New Madrid earthquake sequence have indeed assumed ongoing microseismicity to be a continuing aftershock sequence. Here we use instrumentally recorded aftershock locations and models of elastic stress change to develop a kinematically consistent rupture scenario for three of the four largest earthquakes of the 1811-1812 New Madrid sequence. Our results suggest that these three events occurred on two contiguous faults, producing lobes of increased stress near fault intersections and end points, in areas where present-day microearthquakes have been hitherto interpreted as evidence of primary mainshock rupture. We infer that the remaining New Madrid mainshock may have occurred more than 200 km north of this region in the Wabash Valley of southern Indiana and Illinois--an area that contains abundant modern microseismicity, and where substantial liquefaction was documented by historic accounts. Our results suggest that future large mid-plate earthquake sequences may extend over a much broader region than previously suspected.     

Earthquakes as beacons of stress change

Aftershocks occurring on faults in the far-field of a large earthquake rupture can generally be accounted for by changes in static stress on these faults caused by the rupture. This implies that faults interact, and that the timing of an earthquake can be affected by previous nearby ruptures. Here we explore the potential of small earthquakes to act as 'beacons' for the mechanical state of the crust. We investigate the static-stress changes resulting from the 1992 Landers earthquake in southern California which occurred in an area of high seismic activity stemming from many faults. We first gauge the response of the regional seismicity to the Landers event with a new technique, and then apply the same method to the inverse problem of determining the slip distribution on the main rupture from the seismicity. Assuming justifiable parameters, we derive credible matches to slip profiles obtained directly from the Landers mainshock. Our results provide a way to monitor mechanical conditions in the upper crust, and to investigate processes leading to fault failure.     

The deterministic nature of earthquake rupture

Understanding the earthquake rupture process is central to our understanding of fault systems and earthquake hazards. Multiple hypotheses concerning the nature of fault rupture have been proposed but no unifying theory has emerged. The conceptual hypothesis most commonly cited is the cascade model for fault rupture. In the cascade model, slip initiates on a small fault patch and continues to rupture further across a fault plane as long as the conditions are favourable. Two fundamental implications of this domino-like theory are that small earthquakes begin in the same manner as large earthquakes and that the rupture process is not deterministic--that is, the size of the earthquake cannot be determined until the cessation of rupture. Here we show that the frequency content of radiated seismic energy within the first few seconds of rupture scales with the final magnitude of the event. We infer that the magnitude of an earthquake can therefore be estimated before the rupture is complete. This finding implies that the rupture process is to some degree deterministic and has implications for the physics of the rupture process.     

Seismology: earthquake risk on the Sunda trench

On 28 March 2005 the Sunda megathrust in Indonesia ruptured again, producing another great earthquake three months after the previous one. The rupture was contiguous with that of the December 2004 Sumatra-Andaman earthquake, and is likely to have been sparked by local stress, although the triggering stresses at its hypocentre were very small - of the order of just 0.1 bar. Calculations show that stresses imposed by the second rupture have brought closer to failure the megathrust immediately to the south, under the Batu and Mentawai islands, and have expanded the area of increased stress on the Sumatra fault. Palaeoseismologic studies show that the Mentawai segment of the Sunda megathrust is well advanced in its seismic cycle and is therefore a good candidate for triggered failure.     

四川芦山Ms 7.0级地震导致周边断层的应力变化

 2013年4月20日,四川雅安市芦山县发生M_s 7.0级强烈地震,造成重大的人员伤亡和经济损失。地震造成的区域库仑应力变化、对周围断层的影响及后续地震带发展趋势是应该关注的问题。利用USGS震源机制解,根据地震静态触发原理,基于弹性位错理论和分层地壳模型,计算得出芦山地震引起同震库仑应力变化从断层的1.0MPa量级减小到200km外的0.1kPa,研究认为在地震之后大部分区域的应力得到释放,鲜水河断裂道孚—康定段和玉龙希断裂南段危险性增加。     

Near-field postseismic deformation along the rupture of 2008 Wenchuan earthquake and its implications

At first time, we observed the postseismic deformation, which actually is a kind of creep after slip, along the intra-plate active fault associated with the Wenchuan earthquake. To understand the near-field postseismic deformation following 2008 Wenchuan earthquake in Sichuan Province of China, we compared the fault scarp or flexure scarp profiles measured in different campaigns. Our result shows that among total 19 observation sites, on 13 sites (68% of 19) fault scarps fall back with an average about 9.7% decrease; on 5 sites (26% of 19) fault scarps present no change; and on one site (6% of 19) fault scarp continues to uplift with 12.8% increase. The variety of fault scarp we observed results mainly from near-field postseismic deformation, after slip occurred in shallow. Based on our observations, the following are demonstrated: except for the southwestern end near Yingxiu Town where coseismic slip deficit and some elastic energy residue exist there, fall back (68%) or non-changing (26%) of fault scarp shows energy balance or energy deficit due to overthrust, implying that the likelihood of occurrence of strong aftershocks of ≥ M7 becomes very small in these energy-released areas. Moreover, we suggest that a minimum of 10% error due to near-field postseismic deformation should be considered when evaluating the magnitude of historic and paleoearthquake or slip rate based on the fault scarp dis-placement, even though the error caused by erosion has been accounted already.     

Plate-boundary deformation associated with the great Sumatra-Andaman earthquake

The Sumatra-Andaman earthquake of 26 December 2004 is the first giant earthquake (moment magnitude M(w) > 9.0) to have occurred since the advent of modern space-based geodesy and broadband seismology. It therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events. Here we report estimates of the ground displacement associated with this event, using near-field Global Positioning System (GPS) surveys in northwestern Sumatra combined with in situ and remote observations of the vertical motion of coral reefs. These data show that the earthquake was generated by rupture of the Sunda subduction megathrust over a distance of >1,500 kilometres and a width of <150 kilometres. Megathrust slip exceeded 20 metres offshore northern Sumatra, mostly at depths shallower than 30 kilometres. Comparison of the geodetically and seismically inferred slip distribution indicates that approximately 30 per cent additional fault slip accrued in the 1.5 months following the 500-second-long seismic rupture. Both seismic and aseismic slip before our re-occupation of GPS sites occurred on the shallow portion of the megathrust, where the large Aceh tsunami originated. Slip tapers off abruptly along strike beneath Simeulue Island at the southeastern edge of the rupture, where the earthquake nucleated and where an M(w) = 7.2 earthquake occurred in late 2002. This edge also abuts the northern limit of slip in the 28 March 2005 M(w) = 8.7 Nias-Simeulue earthquake.     

Lushan M S7.0 earthquake: A blind reserve-fault event

In the epicenter of the Lushan M _S7.0 earthquake there are several imbricate active reverse faults lying from northwest to southeast, namely the Gengda-Longdong, Yanjing-Wulong, Shuangshi-Dachuan and Dayi faults. Emergency field investigations have indicated that no apparent earthquake surface rupture zones were located along these active faults or their adjacent areas. Only brittle compressive ruptures in the cement-covered pavements can be seen in Shuangshi, Taiping, Longxing and Longmen Townships, and these ruptures show that a local crustal shortening occurred in the region during the earthquake. Combining spatial distribution of the relocated aftershocks and focal mechanism solutions, it is inferred that the Lushan earthquake is classified as a typical blind reverse-fault earthquake, and it is advised that the relevant departments should pay great attention to other historically un-ruptured segments along the Longmenshan thrust belt and throughout its adjacent areas.