Wenchuan Earthquake Deformation 3D Modelling based on ALOS/PALSAR Data

A devastating earthquake of magnitude Mw 7.9 occurred in Wenchuan area of Sichuan Province, China on 12th May 2008 and caused great casualties and economic damage. This study is aiming to investigate the faulting geometry and motion of the major seismic faults in Longmenshan fault thrust belt that caused this earthquake, based on the surface rupture displacement data measured using differential interferometric synthetic aperture radar (DInSAR) and SAR amplitude pixel-offset techniques. The cross-event Japanese ALOS PALSAR data have been used for this study. First, the methodology for recovering the missing data in the decoherence zone of the DInSAR line-of-sight (LOS) surface motion maps was developed. In the area along the seismic fault zone, the coherence between pre- and post-event SAR images is completely lost because of the earthquake induced violent and chaotic destruction on the land surface and as the result, no surface displacement can be measured using the DInSAR technique. An Adaptive Local Kriging (ALK) technique has then been developed to retrieve the interferometric fringe patterns in the decoherence zone. The novel ALK operating in a multi-step approach enables to retrieve and interpolate the values with high fidelity to the original dataset. Thus a map of continuous radar LOS displacement was generated. Then, the horizontal displacement motion maps in ground range and azimuth direction were derived from cross-event SAR amplitude image pairs using advanced sub-pixel offset technique, Phase Correlation based Image Analysis System (PCIAS). Though the ground range pixel-offset is proportional to the LOS displacement, the azimuth pixel-offset data provide extra information of the coseismic motion. Thus the horizontal displacement vector field can be obtained in order to constrain the faulting motions in key areas. Finally, with the constraints of the ALK refined DInSAR data and the horizontal displacement data together with the published seismic focal mechanism solutions, seismic reflection profiles and field observations, forward modelling was proceeded using the Poly3D software to decide the most likely faulting geometry based on the optimal matching between the simulated and the measured surface displacement. In the much disputed Beichuan – Pengguan area, the best fit is achieved only when the Pengguan fault is set as the main fault that intercept the Yingxiu-Beichuan fault at a depth of about 13 kilometres. This geometric relationship between the two faults and the distribution of slip is compatible with them being two adjacent splay faults on a propagating thrust

Patterns and mechanisms of coseismic and postseismic slips of the 2011 M W 7.1 Van (Turkey) earthquake revealed by multi-platform synthetic aperture radar interferometry

On 23rd October 2011, a MW 7.1 reverse slip earthquake occurred in the Bardakçı-Saray thrust fault zone in the Van region, Eastern Turkey. Earlier geodetic studies have found different slip distributions in terms of both magnitude and pattern. In this paper, we present several COSMO-SkyMED (CSK), Envisat ASAR and RADARSAT-2 interferograms spanning different time intervals, showing that significant postseismic signals can be observed in the first three days after the mainshock. Using observations that combine coseismic and postseismic signals is shown to significantly underestimate coseismic slip. We hence employed the CSK pair with the minimum postseismic signals to generate one conventional interferogram and one along-track interferogram for further coseismic modelling. Our best-fit coseismic slip model suggests that: (1) this event is associated with a buried NNW dipping fault with a preferable dip angle of 49° and a maximum slip of 6.5 m at a depth of 12 km; and (2) two unequal asperities can be observed, consistent with previous seismic solutions. Significant oblique aseismic slip with predominant left-lateral slip components above the coseismic rupture zone within the first 3 days after the mainshock is also revealed by a postseismic CSK interferogram, indicating that the greatest principal stress axis might have rotated due to a significant stress drop during the coseismic rupture

Modelling co- and post-seismic displacements revealed by InSAR, and their implications for fault behaviour

The ultimate goal of seismology is to estimate the timing, magnitude and potential spatial extent of future seismic events along pre-existing faults. Based on the rate-state friction law, several theoretical physical earthquake models have been proposed towards this goal. Tectonic loading rate and frictional properties of faults are required in these models. Modern geodetic observations, e.g. GPS and InSAR, have provided unprecedented near-field observations following large earthquakes. In theory, according to the frictional rate and state asperity earthquake model, velocity-weakening regions holding seismic motions on faults should be separated with velocity-strengthening regions within which faults slip only aseismically. However, early afterslip following the 2011 MW 9.1 Tohoku-Oki earthquake revealed from GPS measurements was largely overlaid on the historical rupture zones, which challenged the velocity weakening asperity model. Therefore, the performance of the laboratory based friction law in the natural events needs further investigation, and the factors that may affect the estimates of slip models through geodetic modelling should also be discussed systematically. In this thesis, several moderate-strong events were investigated in order to address this important issue. The best-fit co- and post-seismic slip models following the 2009 MW 6.3 Haixi, Qinghai thrust-slip earthquake determined by InSAR deformation time-series suggest that the maximum afterslip is concentrated in the same area as the coseismic slip model, which is similar to the patterns observed in the 2011 Japan earthquake. In this case, complex geometric asperity may play a vital role in the coseismic nucleation and postseismic faulting. The major early afterslip after the 2011 MW 7.1 Van mainshock, which was revealed by one COSMO-SkyMed postseismic interferogram, is found just above the coseismic slip pattern. In this event, a postseismic modelling that did not allow slip across the coseismic asperity was also tested, suggesting that the slip model without slip in the asperities can explain the postseismic observations as well as the afterslip model without constraints on slip in the asperities. In the 2011 MW 9.1 Tohoku-Oki earthquake, a joint inversion with the GRACE coseismic gravity changes and inland coseismic GPS observations was conducted to re-investigate the coseismic slip model of the mainshock. A comparison of slip models from these different datasets suggests that significant variations of slip models can be observed, particularly the locations of the maximum slips. The joint slip model shows that the maximum slip of ~42 m appears near the seafloor surface close to the Japan Trench. Meanwhile, the accumulative afterslip patterns (slip >2 m) determined in previous studies appear in spatial correlation with the Coulomb stress changes generated using the joint slip model. As a strike-slip faulting event, the 2011 MW 6.8 Yushu earthquake was also investigated through co- and post-seismic modelling with more SAR data than was used in previous study. Best slip models suggest that the major afterslip is concentrated in shallow parts of the faults and between the two major coseismic slip patterns, suggesting that the performance of the rate and state frictional asperity model is appropriate in this event. Other postseismic physical mechanisms, pore-elastic rebound and viscoelastic relaxation have also been examined, which cannot significantly affect the estimate of the shallow afterslip model in this study. It is believed that the shallow afterslip predominantly controlled the postseismic behaviour after the mainshock in this case. In comparison to another 21 earthquakes investigated using geodetic data from other studies, complementary spatial extents between co- and post-seismic slip models can be identified. The 2009 MW 6.3 Qinghai earthquake is an exceptional case, in which the faulting behaviours might be dominated by the fault structure (e.g. fault bending). In conclusion, the major contributions from this thesis include: 1) the friction law gives a first order fit in most of natural events examined in this thesis; 2) geometric asperities may play an important role in faulting during earthquake cycles; 3) significant uncertainties in co- and post-seismic slip models can appreciably bias the estimation of fault frictional properties; 4) new insights derived from each earthquake regarding their fault structures and complex faulting behaviours have been observed in this thesis; and (5) a novel package for geodetic earthquake modelling has been developed, which can handle multiple datasets including InSAR, GPS and land/space based gravity changes

Crustal structure and the seismogenic environment in Yunnan imaged by double-difference tomography

The large-scale faulting and earthquake activities that developed extensively in the Yunnan area are associated with the collision of India and Eurasia. The fine crustal structure can provide a better understanding of the crustal deformation, seismogenic environment, and rupture processes. We performed a new 3-dimensional (3D) P wave velocity structure and seismic relocation using double-difference tomography based on seismic observations. The tomography images show that large-scale low-velocity anomalies spread around the margin of the south Chuan–Dian Block, Xiaojiang fault (XJF), and the Lijiang–Xiaojinhe fault (LJ-XJHF) in the middle and lower crust. There is an obvious high-speed anomaly in the Emeishan large igneous province (ELIP). We infer that the low-velocity anomaly under the LJ-XJHF zone may be derived from the lower crustal flow extruded from the central Tibetan plateau and obstructed by the ELIP, while the velocity anomalies around the XJF might be caused by shear heating, which is associated with the large-deep strike–slip fault and the transmission of stress in the southeast direction. The inversion results also show that the Yangbi earthquake occurred at the NW–SE boundary of high and low velocity from the upper crust to the lower crust, which coincides well with the location of the Yangbi earthquake sequence and the Weixi–Qiaohou fault. Meanwhile, the earthquake relocations show that the aftershocks are mainly distributed at low velocities. All the aforementioned research results indicate that the Yangbi earthquake might be attributed to the intrusion of the soft material flow along the Weixi–Qiaohou fault in the NW–SE direction. These low-viscosity crustal materials would cause brittle fractures and result in NW–SE sinistral strike–slip faults

Measuring Coseismic Deformation With Spaceborne Synthetic Aperture Radar: A Review

In the past 25 years, space-borne Synthetic Aperture Radar imagery has become an increasingly available data source for the study of crustal deformation associated with moderate to large earthquakes (M > 4.0). Coseismic surface deformation can be measured with several well-established techniques, the applicability of which depends on the ground displacement pattern, on several radar parameters, and on the surface properties at the time of the radar acquisitions. The state-of-the-art concerning the measurement techniques is reviewed, and their application to over 100 case-studies since the launch of the Sentinel-1a satellite is discussed, including the performance of the different methods and the data processing aspects, which still constitute topics of ongoing research

Advances in Geotechnical Earthquake Engineering

This book sheds lights on recent advances in Geotechnical Earthquake Engineering with special emphasis on soil liquefaction, soil-structure interaction, seismic safety of dams and underground monuments, mitigation strategies against landslide and fire whirlwind resulting from earthquakes and vibration of a layered rotating plant and Bryan's effect. The book contains sixteen chapters covering several interesting research topics written by researchers and experts from several countries. The research reported in this book is useful to graduate students and researchers working in the fields of structural and earthquake engineering. The book will also be of considerable help to civil engineers working on construction and repair of engineering structures, such as buildings, roads, dams and monuments

The fate of landslide debris after large earthquakes, Wenchuan, China

Large earthquakes are significant geomorphic events in mountain ranges;they can cause uplift while also triggering huge volumes of bedrock landsliding. Constraining theerosion of coseismic landslide deposits is important forunderstandingthe long-term impact of earthquakes on landscape evolution and natural hazards. The 2008 Mw7.9 Wenchuan earthquake provides an excellent opportunity to constrain the processes and rates that erode landslide deposits within and out of mountain catchments. I constructed a large multitemporal landslide inventory of the epicentral regionof the earthquaketo map the erosion of landslide deposits in the 10 years after the earthquake. In these 10 years less than 13% of the sediment hadbeen deposited into the stores associated with the main trunk of the orogen draining river leaving up to 88% stored on the hillslope. Debris flows were the key remobilising process, transporting over 50% of all the sediment that is mobilised. Using slope stability modelling and analysis of the landslide inventory I identified most of the debris flows were triggered on the hill slope and their frequency decreased rapidly through time. Initially debris flows were triggered in small deposits with low drainage areasbut through time the landslides closer to channel network become more likely to be the source of debris flow activity, a pattern that highlights a change in the physical properties of the deposits. I evaluate several theories that may explain this pattern and highlight strategies to identify the stabilising processes. Finally, I present a 0-Dimensional model investigatingthe impact of storage of coseismic sediment on landscape evolution. This model suggests coseismic sediment can remain in the mountain range for timescales greater than the recurrence time of large earthquakes causing a decrease in exhumation ratesas sediment is remobilised. The long storage time of sediment in orogens reduces the impact of earthquakes on long term (>1000 years) exhumation records so that single earthquakes are only recorded in the most local or highest resolution records

Analysis of earthquake signals by spaceborne gravimetry

The Gravity Recovery And Climate Experiment (GRACE) mission was launched on Mar. 17, 2002 and has provided the scientists with the gravity data for nearly ten years. The time variable gravity field provided by the GRACE has improved our knowledge of the earth in many fields such as hydrology, oceanography and glaciology. But compared to those “hot” fields, the publications of GRACE in seismology is considerably less. However, GRACE can provide scientists with an independent observation of the earthquake process. Coincidentally, some of the largest earthquakes are within GRACE’s life span - Sumatra-Andaman Earthquake (Indonesia) 2004, Maule Earthquake (Chile) 2010 and Tohoku Earthquake (Japan) 2011. Furthermore, a smaller earthquake - Sichuan Earthquake (China) 2008 has also been examined to test whether the GRACE can detect earthquakes smaller than Mw = 8.0. Different from the traditional methods of the earthquake researches, the gravity method has its advantages: 1. Massive: global scale; 2. Insight: gravity changes can reveal the underground mass changes which do not cause so much motion on the earth surface; 3. Convenient: superior to the traditional methods, the spaceborne gravimetry can get the data from the ocean and glacier parts. The conditions of the data are different among these four earthquakes. The procedures to eliminate the GRACE observation errors and unwanted geophysical data are necessary. First, the C20 term should be replaced by the Satellite Laser Ranging (SLR) data. Second, the hydrology signal especially in the regions of Chile and Sichuan should be eliminated by the Global Land Data Assimilation System (GLDAS) model. Third, Fan filter or Gauss filter 350 km should be applied. Time series analysis by the two-phase changepoint detection and hypothesis testing are applied for each earthquake which is a point-wise analysis. Least squares adjustment is performed on each point to display the coseismic and postseismic signals. Meanwhile, the surface analysis is done by the Empirical Orthogonal Functions (EOF) as it has a flexible base which can suit the data automatically. Although the observation errors have been removed as much as possible, the limited spatial and time resolutions of the GRACE satellite and to retrieve relatively weak earthquake signal among the strong hydrological signals are still problems in the analysis. GRACE can detect some of the large earthquakes, but it depends on the earthquake type, area and the length of the time-series before and after the earthquake. Both coseismic signal and postseismic signal are detected in Sumatra-Andaman Earthquake. Meanwhile, there is no significant coseismic signal in the time series of Sichuan Earthquake, but the EOF detects suspicious earthquake signal in mode 2 with the magnitude less than 1 µGal. For Maule Earthquake, only the coseismic signal is detected. Due to the limited dataset, the detection of the coseismic signal is successful but the postseismic signal is not long enough to be detected in Tohoku Earthquake. However, the different filters will affect the magnitude of the gravity change, so the real gravity changes of those four areas are still under debate. Last, EOF can be used for the separation of the earthquake signals. Compared to other geodetic technics the gravity method can detect the signals underground and in the ocean areas. The coseismic and postseismic signals detected by GRACE show underground processes of the earthquakes which can help scientists better understand the earthquake mechanism and will contribute to the earthquake prediction in the future