Single Fault Model Research Articles

A New Method to Invert for Interseismic Deep Slip Along Closely Spaced Faults Using Surface Velocities and Subsurface Stressing‐Rate Tensors

AbstractInversions of interseismic geodetic surface velocities often cannot uniquely resolve the three‐dimensional slip‐rate distribution along closely spaced faults. Microseismic focal mechanisms reveal stress information at depth and may provide additional constraints for inversions that estimate slip rates. Here, we present a new inverse approach that utilizes both surface velocities and subsurface stressing‐rate tensors to constrain interseismic slip rates and activity of closely spaced faults. We assess the ability of the inverse approach to recover slip rate distributions from stressing‐rate tensors and surface velocities generated by two forward models: (a) a single strike‐slip fault model and (b) a complex southern San Andreas fault system (SAFS) model. The single fault model inversions reveal that a sparse array of regularly spaced stressing‐rate tensors can recover the forward model slip distribution better than surface velocity inversions alone. Because focal mechanism inversions currently provide normalized deviatoric stress tensors, we perform inversions for slip rate using full, deviatoric or normalized deviatoric forward‐model‐generated stressing‐rate tensors to assess the impact of removing stress magnitude from the constraining data. All the inversions, except for those that use normalized deviatoric stressing‐rate tensors, recover the forward model slip‐rate distribution well, even for the SAFS model. Jointly inverting stressing rate and velocity data best recovers the forward model slip‐rate distribution and may improve estimates of interseismic deep slip rates in regions of complex faulting, such as the southern SAFS; however, successful inversions of crustal data will require methods to estimate stressing‐rate magnitudes.

Study on Evolution Characteristics of Stress Field of Stepped Fault

With the development of mining to the deep, the ground stress increases correspondingly, and the geological conditions become more and more complicated. Mine dynamic disasters mainly occur near faults. Numerous studies have been conducted on single faults, but few studies have been done on stepped faults formed by double faults under complex geological conditions. Based on this, with the 3110-working face of Zhong Tai Mining in He bi as the engineering background, Single fault model a and different step fault models b, c, d, e and f were established by Flac3D numerical simulation software. The characteristics of stress field evolution of single fault and step fault under the influence of deep mining are studied. The results show that the stress distribution in the upper wall of the fault forms a stress concentration zone, and in the lower wall forms a stress weakening zone. Under the influence of excavation disturbance, the vertical fault distance and the horizontal distance between faults will affect the stress evolution characteristics of the stepped fault.

Interseismic slip rate and fault geometry along the northwest Himalaya

SUMMARY Geodetic networks enable us to investigate interseismic crustal deformation along the northwest Himalaya. Using 144 GNSS surface velocities and a Bayesian inversion model, we estimate the slip rate and fault geometry of the Main Himalayan Thrust (MHT) along six arc-normal transects in the northwest Himalaya. We consider that the fault plane consists of three sections along the décollement, namely the locking zone (0−12 km), the transition zone (10−22 km) and the creeping zone (≥22 km). The MHT is found to be completely locked from the surface down to an average depth of 6 ± 2 km. The locking-to-creeping transition zone along the décollement extends from the edge of the fully locked area to a deeper depth (14 ± 3 km) to the tip of the creeping zone of the MHT (17 ± 2 km) with a slip rate of 1.6 ± 0.9 to 3.7 ± 1.1 mm yr−1. Considering the range of uncertainties between 1−2 mm yr−1 for the GNSS velocities, the inverted slip rate along the transition zone of MHT turns out to be insignificant. Thus, the locking zone along the northwest Himalaya extends from the MFT to ∼111 ± 6 km in the north with a locking depth of ∼17 ± 2 km. The deeper part of the MHT is inferred to be creeping with an average slip rate of ∼19.1 ± 1.9 mm yr−1 along the northwest Himalaya. In addition, we have also illustrated a splay-fault model to account for the fault kinematics along the splay faults and the main décollement. The splay-fault model indicates a distributed slip rate at the locking-to-creeping transition zone and about ∼15 per cent smaller slip rate of the MHT than that of the single-fault model. Further, the checkerboard test and the uniform slip model exhibit the reliability of the current GNSS network and the inversion model (single- and splay-fault models). Overall, the updated fault kinematics inevitably contribute to the improvement of seismic hazard evaluation along the northwest Himalaya.

Exploring the Impact of Soft Errors on the Reliability of Real-Time Embedded Operating Systems

The continuous scaling of electronic components has led to the development of high-performance microprocessors that are suitable even for safety-critical applications where radiation-induced errors such as Single Event Effects (SEEs) can have a significant impact on the performance and reliability of the system. This work is dedicated to investigating the reliability of systems based on programmable hardware and Real-time operating Systems (RTOS) in the presence of architectural faults induced by soft errors in the configuration memory of the programmable hardware. We performed a proton radiation test campaigned at PSI radiation facility to identify the fault model affecting the configuration memory of Xilinx Zynq-7020 reconfigurable AP-Soc Device. The identified fault model in terms of SEU and MBU clusters has been used to evaluate the impact of proton-induced faults on applications running within FreeRTOS on a Microblaze soft processor. A Single Event Multiple Upset fault model resulting from a proton test is presented, focusing on characteristics such as shape, size, and frequency of observed cluster of errors. We conduct two fault injection campaigns and analyze the results to assess the effect of cluster size on system reliability. Moreover, we discuss software exceptions caused by faults that can affect the hardware structure of the soft processor.

Fine divided method of active tectonic blocks based on adaptive density peak clustering and improved Euler model, revealed by GNSS velocities in Southwest Japan

The rational division of active tectonic blocks is an important premise for distinguishing active faults, discovering hidden faults, and an in-depth understanding of regional crustal deformation characteristics, seismicity, and tectonic dynamics. In recent years, clustering algorithms have been widely used in active block division owing to their high efficiency and reliability. However, the classical clustering algorithm sensitivity to the block boundary under different conditions is inconsistent. Additionally, the Euler model constraint (EMC) process involved in the classical clustering algorithm ignores the influence of the boundary points and gross error points of the active block. Therefore, in this study, a new active block division method based on adaptive density peak clustering (ADPC) and an improved Euler model constraint (IEMC) is proposed. The simulation test results show that compared with the traditional clustering algorithm, the mean error rate of the new method can be reduced by up to 40% in the simple single-fault model and by 10% in the complex three-segment model. The new method was further applied to an actual case of an active block division in southwest Japan. The results further show that the new method is less affected by the regional environmental conditions. Compared with previous results in southwest Japan, the proposed method has higher accuracy and higher reliability of the active block division results. The new method for active block division proposed in this study can be used to accurately explore the location and quantity of block boundaries.

Rapid and quantitative uncertainty estimation of coseismic slip distribution for large interplate earthquakes using real-time GNSS data and its application to tsunami inundation prediction

This study proposes a new method for the uncertainty estimation of coseismic slip distribution on the plate interface deduced from real-time global navigation satellite system (GNSS) data and explores its application for tsunami inundation prediction. Jointly developed by the Geospatial Information Authority of Japan and Tohoku University, REGARD (REal-time GEONET Analysis system for Rapid Deformation monitoring) estimates coseismic fault models (a single rectangular fault model and slip distribution model) in real time to support tsunami prediction. The estimated results are adopted as part of the Disaster Information System, which is used by the Cabinet Office of the Government of Japan to assess tsunami inundation and damage. However, the REGARD system currently struggles to estimate the quantitative uncertainty of the estimated result, although the obtained result should contain both observation and modeling errors caused by the model settings. Understanding such quantitative uncertainties based on the input data is essential for utilizing this resource for disaster response. We developed an algorithm that estimates the coseismic slip distribution and its uncertainties using Markov chain Monte Carlo methods. We focused on the Nankai Trough of southwest Japan, where megathrust earthquakes have repeatedly occurred, and used simulation data to assume a Hoei-type earthquake. We divided the 2951 rectangular subfaults on the plate interface and designed a multistage sampling flow with stepwise perturbation groups. As a result, we successfully estimated the slip distribution and its uncertainty at the 95% confidence interval of the posterior probability density function. Furthermore, we developed a new visualization procedure that shows the risk of tsunami inundation and the probability on a map. Under the algorithm, we regarded the Markov chain Monte Carlo samples as individual fault models and clustered them using the k-means approach to obtain different tsunami source scenarios. We then calculated the parallel tsunami inundations and integrated the results on the map. This map, which expresses the uncertainties of tsunami inundation caused by uncertainties in the coseismic fault estimation, offers quantitative and real time insights into possible worst-case scenarios.Graphical

Real-time automatic uncertainty estimation of coseismic single rectangular fault model using GNSS data

Rapid estimation of the coseismic fault model for medium-to-large-sized earthquakes is key for disaster response. To estimate the coseismic fault model for large earthquakes, the Geospatial Information Authority of Japan and Tohoku University have jointly developed a real-time GEONET analysis system for rapid deformation monitoring (REGARD). REGARD can estimate the single rectangular fault model and slip distribution along the assumed plate interface. The single rectangular fault model is useful as a first-order approximation of a medium-to-large earthquake. However, in its estimation, it is difficult to obtain accurate results for model parameters due to the strong effect of initial values. To solve this problem, this study proposes a new method to estimate the coseismic fault model and model uncertainties in real time based on the Bayesian inversion approach using the Markov Chain Monte Carlo (MCMC) method. The MCMC approach is computationally expensive and hyperparameters should be defined in advance via trial and error. The sampling efficiency was improved using a parallel tempering method, and an automatic definition method for hyperparameters was developed for real-time use. The calculation time was within 30 s for 1 × 106 samples using a typical single LINUX server, which can implement real-time analysis, similar to REGARD. The reliability of the developed method was evaluated using data from recent earthquakes (2016 Kumamoto and 2019 Yamagata-Oki earthquakes). Simulations of the earthquakes in the Sea of Japan were also conducted exhaustively. The results showed an advantage over the maximum likelihood approach with a priori information, which has initial value dependence in nonlinear problems. In terms of application to data with a small signal-to-noise ratio, the results suggest the possibility of using several conjugate fault models. There is a tradeoff between the fault area and slip amount, especially for offshore earthquakes, which means that quantification of the uncertainty enables us to evaluate the reliability of the fault model estimation results in real time.

Seismogenic potential of the High Durance Fault constrained by 20 yr of GNSS measurements in the Western European Alps

SUMMARYDue to the steady moderate seismicity observed along the Briançon seismic arc, in the south-western French Alps, three temporary GNSS (Global Navigation Satellite System) surveys took place in 1996, 2006 and 2011, across a ∼50 × 60 km² wide area, to investigate the surface deformation field. The horizontal velocity field computed from these three surveys showed an east–west extension in the network. A fourth campaign was led in 2016, creating a 20 yr observation span, resulting in measurements which reach a sufficient accuracy to assess whether extension found within the Briançon network is localized onto any particular tectonic feature. Several faults in this area are known to be active normal faults. Assessing the localization of the deformation may lead to a better understanding of the active tectonics of the Alpine belt. To address this issue, a robust velocity field was computed from the combination of the different campaign and permanent GNSS data. Strain rate tensors were derived for the first time in this area on a 0.1 × 0.1 deg grid to assess the distribution of the deformation. The regional deformation appears localized in the Briançon area and reaches up to 20 ± 5 nanostrain yr−1 in the centre of the network. The observed velocities were projected on a profile across the network and compared with modelled interseismic deformation to characterize the behaviour of the major active faults known in the study zone. While a two-fault model provides the best fit to the data, a single fault model has only marginally higher residuals, with parameters which are more consistent with the seismotectonics of the region. The localization of the single modelled fault is consistent with the location of the High Durance Fault (HDF). Therefore, we used the known geological location of this structure as a priori information in a block model to compute a fault slip rate at the interface between the two blocks. The velocities on the interface indicate 0.4–0.5 mm yr−1 of extension, and therefore strain accumulates along the HDF throughout the seismic cycle. The geodetically derived fault slip rate is converted into an equivalent seismic moment release rate, which is consistent within its uncertainty bounds with the known historical and instrumental seismicity of the Briançon area.

Investigating thrust-fault growth and segment linkage using displacement distribution analysis in the active Duzhanzi thrust fault zone, Northern Tian Shan of China

Few studies investigate the detailed character of thrust-fault population or their along-strike fault-displacement patterns. Our study presents new displacement-length data collected from well-exposed thrust faults within the Dushanzi Thrust Fault Zone (DTFZ), an active thrust fault-fold system in the Northern Tianshan of China. Numerous surface-rupturing thrust faults cut through upper Quaternary glacial outwash terraces and younger alluvial fans. Based on 0.5 m-resolution LiDAR-derived DEMs, we measured along-strike displacement and lengths of ~30 fault strands to construct a three-dimensional deformation field recorded by these deformed marker surfaces. Three fault segments constitute the DTFZ, and each presents different geometric structures and displacement patterns. Through detailed analysis of the displacement distribution and geometrical attributes of thrust faults, we constrain the growth and propagation history of faults and of the systems that they form. The scaling relationship between maximum displacement (Dmax) and fault length (L) measured from thrust faults within the DTFZ displays similar linear scaling to other thrust faults, but different Dmax/L ratios due to underestimated displacement of the DTFZ. Two conceptual models: a single-fault model and segment-linkage model are explored to explain the growth of faults or fault systems and the changes of Dmax/L ratios during fault evolution.

Performance of uniform and heterogeneous slip distributions for the modeling of the November 2016 off Fukushima earthquake and tsunami, Japan

The Mw 6.9 earthquake off Fukushima Prefecture, Japan, of 22 November 2016 was followed by a tsunami that struck the Japanese coast from Hokkaido in northern Japan to Wakayama Prefecture in western Japan. We compared the performance of a seismologically deduced single-fault model, a seismologically deduced finite fault slip model (FFM), an optimized single-fault model based on tsunami data, the FFM with horizontal shift, and the tsunami waveform inversion models of the previous studies considered for this earthquake regarding reproduction of tsunami waves by tsunami computations. It is important to discuss how these models work well because it is sometimes desirable to obtain an earthquake source model to estimate tsunami waves with a simple process obtained with limited data from the viewpoint of tsunami prediction. The seismologically deduced FFM has an advantage in terms of the information of slip regions of fault plane and was superior to the seismologically deduced single-fault model, especially in predicting amplitudes of tsunami waves. This means that when only with seismic data, the FFM could narrow the range of forecast of tsunami amplitude. In the comparison of models optimized with tsunami data, the single-fault model showed the almost equivalent performance of the tsunami waveform inversion models of previous studies regarding the waveform coincidence with observations and the horizontal location at the negative peak of the initial sea surface displacement. In case the main generation region of the tsunami is concentrated in one place, the tsunamis can be expressed by a single-fault model by conducting the detailed grid search. We also confirmed that the centroid location of centroid moment tensor (CMT) solution and the absolute location of the FFM were not necessarily suitable to express tsunamis, while the moment magnitude, the focal mechanism, the centroid depth of CMT solution, and the relative slip distribution of the FFM were effective to represent tsunamis. Since this event occurred at the shallow depth, the speed of tsunami wave is particularly slow. Therefore, it would be advisable to pay attention to the horizontal uncertainty to apply seismologically obtained solution to tsunami forecast, especially when a tsunami occurs in shallow water.

Testing Different Tectonic Models for the Source of the Mw 6.5, 30 October 2016, Norcia Earthquake (Central Italy): A Youthful Normal Fault, or Negative Inversion of an Old Thrust?

AbstractWe adopted a multidisciplinary approach to investigate the seismotectonic scenario of the 30 October 2016, Mw 6.5, Norcia earthquake, the largest shock of the 2016–2017 central Italy earthquake sequence. First, we used seismological and geodetic data to infer the dip of the main slip patch of the seismogenic fault that turned out to be rather low‐angle (~37°). To evaluate whether this is an acceptable dip for the main seismogenic source, we modeled earthquake deformation using single‐ and multiple‐fault models deduced from aftershock pattern analyses. These models show that the coseismic deformation generated by the Norcia earthquake is coherent with slip along a rather shallow‐dipping plane. To understand the geological significance of this solution, we reconstructed the subsurface architecture of the epicentral area. As the available data are not robust enough to converge on a single fault model, we built three different models encompassing all major geological evidence and the associated uncertainties, including the tectonic style and the location of major décollement levels. In all models the structures derived from the contractional phase play a significant role: from controlling segmentation to partially reusing inherited faults, to fully reactivating in extension a regional thrust, geometrically compatible with the source of the Norcia earthquake. Based on our conclusions, some additional seismogenic sources falling in the eastern, external portions of the Apennines may coincide with inherited structures. This may be a common occurrence in this region of the chain, where the inception of extension is as recent as Middle‐Upper Pleistocene.

Modeling High Stress Drops, Scaling, Interaction, and Irregularity of Repeating Earthquake Sequences Near Parkfield

AbstractRepeating earthquake sequences have been actively investigated to clarify many aspects of earthquake physics. The two particularly well‐studied sequences, known as the Los Angeles and San Francisco repeaters, have several intriguing observations, including their long (for the seismic moment) recurrence times that would suggest stress drops of 300 MPa based on typical assumptions, near‐syncronized timing prior to 2004, and higher than typical inferred stress drops (of 25 to 65 MPa, up to 90 MPa locally), but not as high as the recurrence times suggest. Here we show that all these observations are self‐consistent, in the sense that they can be reproduced in a single fault model. The suitable models build on the standard rate‐and‐state fault models, with velocity‐weakening patches imbedded into a velocity‐strengthening region, by adding either enhanced dynamic weakening during seismic slip or elevated normal stress on the patches, or both, to allow for the higher stress drops. Such models are able to match the observed average properties of the San Francisco and Los Angeles repeaters, as well as the overall nontrivial scaling between the recurrence time and seismic moment exhibited by many repeating sequences as a whole, for reasonable parameter choices based on experiments and theoretical studies. These models are characterized by the occurrence of substantial and variable aseismic slip at the locations of the repeating sources, which explains their atypical relation between recurrence interval and seismic moment, induces variability in the repeating source properties as observed, and results in their neither slip‐ nor time‐predictable behavior.

Study of characteristic of tsunami base on the coastal morphology in north Donggala, Central Sulawesi

The northern arm of Sulawesi potentials to generate earthquake and Tsunami due to the existence of subduction zone in sulawesi sea. It makes the North Donggala as an area with active seismicity. One of the earthquake and Tsunami events occurred is the earthquake and tsunami of Toli-Toli 1996 (M 7.9) causing 9 people are killed and severe damage in Tonggolobibi, Siboang, and Balukang. This earthquake induced tsunami runup of 3.4 m and inundated as far as 400 meters. The aims of this study is to predict runup and inundation area using numerical model and to find out the characteristics of Tsunami wave on straight, bay and cape shape coastal morphology and slopes of coastal. The data in this research consist of are the Etopo2 bathymetry data in data obtained from NOAA (National Oceanic and Atmospheric Administration), Toli-toli’s main earthquakes focal mechanism data 1st January1996 from GCMT (Global Centroid Moment Tensor), the data gained from the SRTM (Shuttle Radar Topography Mission) data 30 m and land cover data in 1996 from Ministry of environment and forestry . Single fault model is used to predict the high of tsunami run-up and to inundation area along Donggala coastal area. Its reviewed by morphology of coastal area that higher run up shows occurs at coastline type like bay have higher run up compare to area with cape and straight coastline. The result shows that the slopes have negative or contras correlation with Tsunami runup and its inundation area.

Trial of Quantitative Uncertainty Estimation for the Single Rectangular Coseismic Fault Model Deduced from Real-time GNSS Data Using the MCMC Approach : A Case Study of the 2011 Tohoku-Oki Earthquake

Trial of Quantitative Uncertainty Estimation for the Single Rectangular Coseismic Fault Model Deduced from Real-time GNSS Data Using the MCMC Approach : A Case Study of the 2011 Tohoku-Oki Earthquake

Tsunami Source Inversion Using Tide Gauge and DART Tsunami Waveforms of the 2017 Mw8.2 Mexico Earthquake

On September 8, 2017 (UTC), a normal-fault earthquake occurred 87 km off the southeast coast of Mexico. This earthquake generated a tsunami that was recorded at coastal tide gauge and offshore buoy stations. First, we conducted a numerical tsunami simulation using a single-fault model to understand the tsunami characteristics near the rupture area, focusing on the nearby tide gauge stations. Second, the tsunami source of this event was estimated from inversion of tsunami waveforms recorded at six coastal stations and three buoys located in the deep ocean. Using the aftershock distribution within 1 day following the main shock, the fault plane orientation had a northeast dip direction (strike = 320^{circ }, dip = 77^{circ }, and rake =-92^{circ }). The results of the tsunami waveform inversion revealed that the fault area was 240 km × 90 km in size with most of the largest slip occurring on the middle and deepest segments of the fault. The maximum slip was 6.03 m from a 30 × 30 km2 segment that was 64.82 km deep at the center of the fault area. The estimated slip distribution showed that the main asperity was at the center of the fault area. The second asperity with an average slip of 5.5 m was found on the northwest-most segments. The estimated slip distribution yielded a seismic moment of 2.9 times 10^{21} Nm (Mw = 8.24), which was calculated assuming an average rigidity of 7times 10^{10} N/m2.

Hypocenter relocations and tsunami simulation for the 15 November 2014 Northern Molucca Sea earthquake in Indonesia

A reverse fault earthquake (Mw 7.1) occurred in the Northern Molucca Sea, Indonesia, on 15 November 2014 at 2:31:40 UTC. The earthquake produced small tsunami waves that are recorded at Jailolo (9 cm), Tobelo (1 cm), and Menado (3 cm) tide gauges. The Indonesian Agency for Climatology, Meteorology, and Geophysics (BMKG) issued a timely (5 minutes after the earthquake) tsunami warning for the event. We used the teleseismic double‐difference seismic tomography method (teletomoDD) to relocate the hypocenters of the mainshock and the aftershocks. The relocated hypocenter of the mainshock for the 2014 Northern Molucca Sea earthquake is located at 1.923°N, 126.539°E, and depth of 48.87 km. In general, the relocated aftershock hypocenters are shallower than those from the BMKG catalog. The relocated hypocenters are distributed within a depth range of 6 to 64 km. The aftershock area from the relocated hypocenters is 80 km long and 55 km wide. The estimated seismic moment from the Global CMT solution (GCMT) was 4.75 × 1019 Nm. We simulated the tsunami from fault model of each GCMT nodal plane to find a fault model that can best explain the observed tsunami heights at Jailolo, Tobelo, and Menado tide gauges. The best single fault model for this event is dipping to the west, has fault length, width, and slip amount of 47 km, 25 km, and 1.16 m, respectively. The K value calculated using the observed and simulated tsunami heights for this best model is 1.026, suggests a very good fit to tsunami observations.

Multiplicity of solutions to geophysical inversion reflected by rupture slip distribution of the 2015 Nepal earthquake

The equivalence of geophysical fields, the finiteness of measurements and the measurement errors make the result of geophysical inversion non-unique. For example, the measurements and inversion method used, the priori rupture model determined and the slip distribution smoothing factor selected will have significant influences on the earthquake rupture slip distribution. Using different data and methods, different authors have given different rupture slip distribution models of the 2015 Mw7.9 Nepal earthquake, with the maximum slip ranging from 3.0 m to 6.8 m. In this paper, geometry parameters of the single rectangular fault model in elastic half-space were inferred constraining with the Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar (InSAR) coseismic deformations and bounding the slip with approximate average value; and then, the single rectangular fault was divided into multiple sub-faults, and the final slip smoothing factor, the final slip distribution and the maximum slip were determined with the misfit–roughness tradeoff curve, the cross-validation sum of squares (CVSS) and the third-party observation data or indexes being comprehensively taken into account. The results show that, the rupture of the Nepal earthquake extended by over 100 km east by south. The maximum slip of the earthquake was about 6.5–6.7 m, and most of the slip is confined at depths of 8–20 km, consistent with the depth distribution of aftershocks. The method for reducing the multiplicity of solutions to rupture slip distribution in this paper was ever used in inversion of rupture slip distribution for the 2008 Wenchuan and 2013 Lushan earthquakes, and the third-party measurement – surface dislocation has very large effect on reducing the multiplicity of solutions to inversion of the Wenchuan earthquake. Other priori information or indicators, such as fault strike, dip, earthquake magnitude, seismic activity, Coulomb stress, and seismic period, can be used for beneficial validation of and comparison with inversion results.

COSEISMIC DEFORMATION FIELD AND FAULT SLIP DISTRIBUTION OF THE 2015 CHILE Mw8.3 EARTHQUAKE

On September 16, 2015, a magnitude 8.3 earthquake struck west of Illapel, Chile. We analyzed Sentinel-1A/IW InSAR data on the descending track acquired before and after the Chile Mw8.3 earthquake of 16 September 2015. We found that the coseismic deformation field of this event consists of many semi circular fringes protruding to east in an approximately 300km long and 190km wide region. The maximum coseismic displacement is about 1.33m in LOS direction corresponding to subsidence or westward shift of the ground. We inverted the coseismic fault slip based on a small-dip single plane fault model in a homogeneous elastic half space. The inverted coseismic slip mainly concentrates at shallow depth above the hypocenter with a symmetry shape. The rupture length along strike is about 340 km with maximum slip of about 8.16m near the trench. The estimated moment is 3.126×1021 N.m (Mw8.27)，the maximum depth of coseismic slip near zero appears to 50km. We also analyzed the postseismic deformation fields using four interferograms with different time intervals. The results show that postseismic deformation occurred in a narrow area of approximately 65km wide with maximum slip 11cm, and its predominant motion changes from uplift to subsidence with time. that is to say, at first, the postseismic deformation direction is opposite to that of coseismic deformation, then it tends to be consistent with coseismic deformation.It maybe indicates the differences and changes in the velocity between the Nazca oceanic plate and the South American continental plate.

COSEISMIC DEFORMATION FIELD AND FAULT SLIP DISTRIBUTION OF THE 2015 CHILE Mw8.3 EARTHQUAKE

Abstract. On September 16, 2015, a magnitude 8.3 earthquake struck west of Illapel, Chile. We analyzed Sentinel-1A/IW InSAR data on the descending track acquired before and after the Chile Mw8.3 earthquake of 16 September 2015. We found that the coseismic deformation field of this event consists of many semi circular fringes protruding to east in an approximately 300km long and 190km wide region. The maximum coseismic displacement is about 1.33m in LOS direction corresponding to subsidence or westward shift of the ground. We inverted the coseismic fault slip based on a small-dip single plane fault model in a homogeneous elastic half space. The inverted coseismic slip mainly concentrates at shallow depth above the hypocenter with a symmetry shape. The rupture length along strike is about 340 km with maximum slip of about 8.16m near the trench. The estimated moment is 3.126×1021 N.m (Mw8.27)，the maximum depth of coseismic slip near zero appears to 50km. We also analyzed the postseismic deformation fields using four interferograms with different time intervals. The results show that postseismic deformation occurred in a narrow area of approximately 65km wide with maximum slip 11cm, and its predominant motion changes from uplift to subsidence with time. that is to say, at first, the postseismic deformation direction is opposite to that of coseismic deformation, then it tends to be consistent with coseismic deformation.It maybe indicates the differences and changes in the velocity between the Nazca oceanic plate and the South American continental plate.

Source Model from ALOS-2 ScanSAR of the 2015 Nepal Earthquakes

Abstract The 2015 Gorkha Nepal Earthquake sequence started with a magnitude Mw 7.8 main shock and continued with several large aftershocks, particularly the second major shock of Mw 7.3. Both earthquake events were captured using ALOS-2 ScanSAR images to determine the coseismic surface deformation and the source models. In this paper, the displacement maps were produced and the corresponding modelling results were discussed. The single fault model of the main shock suggests that there was nearly 6 m of right-lateral oblique slip motion with fault struck of 292° and dipped gently Northeast at 7°, indicating that the main shock was on a thrust fault. Moreover, a single fault model for the Mw 7.3 quake with striking of 312° and dipping of 11° was derived from observed result. Both results showed the fault planes struck generally to South and dipped northeast, which depicted the risks since the main shock occurred.