Co-seismic slip distribution of the 2011 Tohoku(MW 9.0)earthquake inverted from GPS and space-borne gravimetric data

Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.

Determination of Smoothing Factor for the Inversion of Co-seismic Slip Distribution

For the determination of the smoothing factor (also known as the regularization parameter) in the co-seismic slip distribution inversion, the compromise curve between the model roughness and the data fitting residual is generally used to determine (in order to distinguish the method proposed in this paper, the method is called “L curve” according to its shape). Based on the L-curve, the Eclectic Intersection curve as a new method is proposed to determine the smoothing factor in this paper. The results of the simulated experiment show that the inversion accuracy of the parameters of the seismic slip distribution with the smoothing factor determined by the Eclectic Intersection curve method is better than that of the L curve method. Moreover, the Eclectic Intersection curve method and the L curve method are used to determine the smoothing factor of L’Aquila earthquake and the Taiwan Meinong earthquake slip distribution inversion respectively, and the inversion results are compared and analyzed. The analysis results show that the L’Aquila and the Taiwan Meinong actual earthquake slip distribution results are in the range of other scholars at home and abroad, and compared with the L curve method, the Eclectic Intersection curve method has advantages of high computation efficiency, no need to depend on data fitting degree and more appropriate of smoothing factor and so on.

Co-seismic Slip Distribution of the 2001 Kokoxili Earthquake Inverted by the GPS Data

The authors analyze co_seismic displacement field derived from the Global Position System (GPS) observations collected before and after the 2001 Kokoxili earthquake, western China. Using the co_seismic displacement data, and constrained with surface rupture data, they invert co_seismic slip distribution along the seismic fault. Their result shows that the earthquake ruptured the upper crust down to a depth of 13.1～22km (at 70% certainty), with its optimal estimate at 16.5km. A 2～3m left_lateral strike slip is resolved between the Sun Lake segment and the west end of the main rupture zone, although surface rupture has not been observed there. The surface rupture of this earthquake is ended at the Sun Lake to the west, but left_lateral slip of 1.5～2.0m seems to exist beyond the east end of surface rupture observed from field geology. Seismic moment release estimated using GPS and surface rupture measurement is 6.0×10 20 N·m, which is in good agreement with the result obtained from seismic wave inversion.

Co-seismic surface displacement of the June 21, 2022 M_(W)6 Khōst,Afghanistan earthquake from InSAR observations

A robust estimation of the earthquake location, seismic moment, and fault geometry is essential for objective seismic hazard assessment. Seismic events in a remote location, specifically in the absence of seismic and GNSS networks, can be investigated effectively using the In SAR-based technique. This study adopts the Differential Interferometric SAR(DIn SAR) technique to quantify the co-seismic surface displacement caused by the June 21, 2022, Khōst M_(W)6 earthquake that occurred along the western plate boundary between the Indian and Eurasian plate. The interferograms show that the maximum surface deformation occurred on the northwest and southwest of the fault line. From coherence, the Line of Sight(LOS) displacement, and the co-seismic surface displacement analysis, it has been observed that surface deformation was most pronounced in the southwest region of the fault line, and the surface has moved to the opposite direction along the fault line, which indicates a sinistral slightly oblique strike-slip movement. This In SAR-based observation appears consistent with the seismic waveforms derived from co-seismic surface displacements. Further, it has been argued that the slip deficit accumulated during the period of the last about 48 years along the frontal region of the northward extension of the Suleiman range and associated fault zone is qualitatively estimated at about 1.5 m, which is consistent with the seismic waveforms derived finite slip model.

Co-seismic fault geometry and slip distribution of the 26 December 2004, giant Sumatra–Andaman earthquake constrained by GPS, coral reef, and remote sensing data

We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured ＊200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude（4.0 9 1022 N m）, and smaller than estimation from normal-mode oscillation data modeling（1.0 9 1023 N m）.

足式机器人SLIP模型向上跳跃台阶的运动控制

由于存在地势起伏,台阶对足式机器人运动稳定性会带来较大挑战.弹簧负载倒立摆模型(SLIP)作为研究足式机器人的优良模板,能否完成向上跳跃台阶的动作与其腿部摆角,起跳位置和跳跃高度都有密切的关系.由于调整模型腿部摆角规律容易引发运动失效,故本文在算法中引入虚拟弹簧腿,根据虚拟弹簧腿的运行规律确定合理起跳位置,根据起跳位置来控制系统跳跃高度进而完成跳跃台阶的动作.最后利用仿真软件进行多组仿真,结果表明本文算法对起跳区间划分合理,对起跳高度控制精准,能够实现SLIP模型跳跃台阶前后的稳定运动.

Co-seismic deformation for the 2015 M_(W)7.8 Gorkha earthquake(Nepal)using near-field GPS data

Seasonal variations and common mode errors affect the precision of the Global Positioning System(GPS)time series.In this paper,we explore to improve the precision of coordinate time series,thereby providing a better detection of weak or transient deformation signals,particularly co-seismic signals.Based on 97 GPS stations,including the campaign and continuous GPS stations in Nepal and southern Tibet,we first consider seasonal variations and common errors,then obtain co-seismic deformation of the 2015 Gorkha earthquake in Nepal and southern Tibet.Our co-seismic rupture model is characterized by a shallow ramp and a deeper detachment fault,in agreement with the relocated aftershock sequence.Our results indicate that the earthquake rupture is mainly distributed in the upper-crustal fault,and the maximum slip is up to 8.0 m at~15.0 km depth located in the approximate-80 km east of the epicenter.The average slip is more than 5 m,and the total modelled magnitude is M_(W)7.84,consistent with the observed seismic moment.Our rupture model for the 2015 Gorkha earthquake suggests that the rupture zone is not only in the upper crustal Main Himalayan Thrust(MHT),but also spreads to the northern segment of the MHT.

Revealing the role of pyramidal slip in the high ductility of Mg-Li alloy

The high ductility of Mg-Li alloy has been mainly ascribed to a high activity of pyramidal<c+a>slip to accommodate plastic strain.In the present study,however,a quantitative analysis reveals that Li-addition can only slightly stimulate the activation of pyramidal<c+a>slip under compression along the normal direction of a hot-rolled Mg-4.5 wt.%Li plate,with a relative activity of approximately 18%.Although the limited activity of pyramidal<c+a>slip alone cannot accommodate a large plastic strain,it effectively reduces the number of{10.11}−{10.12}double twins,which are believed to be favorable sites for crack initiation.The evidently reduced activity of double twins leads to a lower cracking tendency,and therefore improves ductility.

Effect of boundary slip on electroosmotic flow in a curved rectangular microchannel

The aim of this study is to numerically investigate the impact of boundary slip on electroosmotic flow(EOF) in curved rectangular microchannels. Navier slip boundary conditions were employed at the curved microchannel walls. The electric potential distribution was governed by the Poisson–Boltzmann equation, whereas the velocity distribution was determined by the Navier–Stokes equation. The finite-difference method was employed to solve these two equations. The detailed discussion focuses on the impact of the curvature ratio, electrokinetic width, aspect ratio and slip length on the velocity. The results indicate that the present problem is strongly dependent on these parameters. The results demonstrate that by varying the dimensionless slip length from 0.001 to 0.01 while maintaining a curvature ratio of 0.5 there is a twofold increase in the maximum velocity. Moreover, this increase becomes more pronounced at higher curvature ratios. In addition, the velocity difference between the inner and outer radial regions increases with increasing slip length. Therefore, the incorporation of the slip boundary condition results in an augmented velocity and a more non-uniform velocity distribution. The findings presented here offer valuable insights into the design and optimization of EOF performance in curved hydrophobic microchannels featuring rectangular cross-sections.

Coseismic deformation and fault slip distribution of the 2023 M_(W)7.8 and M_(W)7.6 earthquakes in Türkiye

On February 6,2023,a devastating earthquake with a moment magnitude of M_(W)7.8 struck the town of Pazarcik in south-central Türkiye,followed by another powerful earthquake with a moment magnitude of M_(W)7.6 that struck the nearby city of Elbistan 9 h later.To study the characteristics of surface deformation caused by this event and the influence of fault rupture,this study calculated the static coseismic deformation of 56 stations and dynamic displacement waveforms of 15 stations using data from the Turkish national fixed global navigation satellite system(GNSS)network.A maximum static coseismic displacement of 0.38 m for the M_(W)7.8 Kahramanmaras earthquake was observed at station ANTE,36 km from the epicenter,and a maximum dynamic coseismic displacement of 4.4 m for the M_(W)7.6 Elbistan earthquake was observed at station EKZ1,5 km from the epicenter.The rupture-slip distributions of the two earthquakes were inverted using GNSS coseismic deformation as a constraint.The results showed that the Kahramanmaras earthquake rupture segment was distinct and exposed on the ground,resulting in significant rupture slip along the Amanos and Pazarcik fault segments of the East Anatolian Fault.The maximum slip in the Pazarcik fault segment was 10.7 m,and rupture occurred at depths of 0–15 km.In the Cardak fault region,the Elbistan earthquake caused significant ruptures at depths of 0–12 km,with the largest amount of slip reaching 11.6 m.The Coulomb stress change caused by the Kahramanmaras earthquake rupture along the Cardak fault segment was approximately 2 bars,and the area of increased Coulomb stress corresponded to the subsequent rupture region of the M_(W)7.6 earthquake.Thus,it is likely that the M_(W)7.8 earthquake triggered or promoted the M_(W)7.6 earthquake.Based on the cumulative stress impact of the M_(W)7.8 and M_(W)7.6 events,the southwestern segment of the East Anatolian Fault,specifically the Amanos fault segment,experienced a Coulomb rupture stress change exceeding 2 bars,warranting further attention to assess its future seismic hazard risk.

Constraints on rupture speed of the 2001 MS 8.1 West Kunlun Mountain Pass earthquake by co-seismic surface rupture slip displacements based on the slip-weakening mechanism with frictional undershoot involved

With co-seismic surface rupture slip displacements provided by the field observation for the 2001 MS8.1 West Kunlun Mountain Pass earthquake, this paper estimates the rupture speed on the main faulting segment with a long straight fault trace on the surface based on a simple slip-weakening rupture model, in which the frictional overshoot or undershoot are involved in consideration of energy partition during the earthquake faulting. In contrast to the study of Bouchon and Vallée, in which the rupture propagation along the main fault could exceed the local shear-wave speed, perhaps reach the P-wave speed on a certain section of fault, our results show that, under a slip-weakening assumption combined with a frictional undershoot (partial stress drop model), average rupture speed should be equal to or less than the Rayleigh wave speed with a high seismic radiation efficiency, which is consistent with the result derived by waveform inversion and the result estimated from source stress field. Associated with the surface rupture mechanism, such as partial stress drop (frictional undershoot) associated with the apparent stress, an alternative rupture mechanism based on the slip-weakening model has also been discussed.

Deformation of the Most Recent Co-seismic Surface Ruptures Along the Garzê–Yushu Fault Zone(Dangjiang Segment)and Tectonic Implications For the Tibetan Plateau

The Garzê–Yushu strike-slip fault in central Tibet is the locus of strong earthquakes（M 〉 7）. The deformation and geometry of the co-seismic surface ruptures are reflected in the surface morphology of the fault and depend on the structure of the upper crust as well as the pre-existing tectonics. Therefore, the most recent co-seismic surface ruptures along the Garzê–Yushu fault zone（Dangjiang segment） reveal the surface deformation of the central Tibetan Plateau. Remote sensing images and field investigations suggest a 85 km long surface rupture zone（striking NW-NWW）, less than 50 m wide, defined by discontinuous fault scarps, right-stepping en echelon tensional cracks and left-stepping mole tracks that point to a left-lateral strike-slip fault. The gullies that cross fault scarps record systematic left-lateral offsets of 1.8 m to 5.0 m owing to the most recent earthquake, with moment magnitude of about M 7.5, in the Dangjiang segment. Geological and geomorphological features suggest that the spatial distribution of the 1738 co-seismic surface rupture zone was controlled by the pre-existing active Garzê–Yushu fault zone（Dangjiang segment）. We confirm that the Garzê–Yushu fault zone, a boundary between the Bayan Har Block to the north and the Qiangtang Block to the south, accommodates the eastward extrusion of the Tibetan Plateau and generates strong earthquakes that release the strain energy owing to the relative motion between the Bayan Har and Qiangtang Blocks.

Surface Rupture and Co-seismic Displacement Produced by the Ms 8.0 Wenchuan Earthquake of May ^(12)th,2008,Sichuan,China:Eastwards Growth of the Qinghai-Tibet Plateau

An earthquake of Ms 8 struck Wenchuan County, western Sichuan, China, on May 12^th, 2008 and resulted in long surface ruptures （〉300 km）. The first-hand observations about the surface ruptures produced by the earthquake in the worst-hit areas of Yingxiu, Beichuan and Qingchuan, ascertained that the causative structure of the earthquake was in the central fault zones of the Longmenshan tectonic belt. Average co-seismic vertical displacements along the individual fault of the Yingxiu-Beiehuan rupture zone reach 2.514 m and the cumulative vertical displacements across the central and frontal Longmenshan fault belt is about 5-6 m. The surface rupture strength was reduced from north of Beichuan to Qingchuan County and shows 2-3 m dextral strike-slip component. The Wenchuan thrust-faulting earthquake is a manifestation of eastward growth of the Tibetan Plateau under the action of continuous convergence of the Indian and Eurasian continents.

Analysis of the Wenchuan Ms8.0 Earthquake's Co-seismic Stress and Displacement Change by Using the Finite Element Method

The variation of in situ stress before and after earthquakes is an issue studied by geologists. In this paper, on the basis of the fault slip dislocation model of Wenchuan Ms8.0 earthquake, the changes of co-seismic displacement and the distribution functions of stress tensor around the Longmen Shan fault zone are calculated. The results show that the co-seismic maximum surface displacement is 4.9 m in the horizontal direction and 6.5 m in the vertical direction, which is almost consistent with the on-site survey and GPS observations. The co-seismic maximum horizontal stress in the hanging wall and footwall decreased sharply as the distance from the Longmen Shan fault zone increased. However, the vertical stress and minimum horizontal stress increased in the footwall and in some areas of the hanging wall. The study of the co-seismic displacement and stress was mainly focused on the long and narrow region along the Longmen Shan fault zone, which coincides with the distribution of the earthquake aftershocks. Therefore, the co-seismic stress only affects the aftershocks, and does not affect distant faults and seismic activities. The results are almost consistent with in situ stress measurements at the two sites before and after Wenchuan Ms8.0 earthquake. Along the fault plane, the co-seismic shear stress in the dip direction is larger than that in the strike direction, which indicates that the faulting mechanism of the Longmen Shan fault zone is a dominant thrust with minor strike-slipping. The results can be used as a reference value for future studies of earthquake mechanisms.

Co-seismic Faults and Geological Hazards and Incidence of Active Fault of Wenchuan Ms 8.0 Earthquake,Sichuan,China

There are two co-seismic faults which developed when the Wenchuan earthquake happened. One occurred along the active fault zone in the central Longmen Mts. and the other in the front of Longmen Mts. The length of which is more than 270 kin and about 80 km respectively. The co-seismic fault shows a reverse flexure belt with strike of N45°-60°E in the ground, which caused uplift at its northwest side and subsidence at the southeast. The fault face dips to the northwest with a dip angle ranging from 50° to 60°. The vertical offset of the co-seismic fault ranges 2.5-3.0 m along the Yingxiu- Beichuan co-seismic fault, and 1.5-1.1 m along the Doujiangyan-Hanwang fault. Movement of the coseismic fault presents obvious segmented features along the active fault zone in central Longmen Mts. For instance, in the section from Yingxiu to Leigu town, thrust without evident slip occurred; while from Beichuan to Qingchuan, thrust and dextral strike-slip take place. Main movement along the front Longmen Mts. shows thrust without slip and segmented features. The area of earthquake intensity more than IX degree and the distribution of secondary geological hazards occurred along the hanging wall of co-seismic faults, and were consistent with the area of aftershock, and its width is less than 40km from co-seismic faults in the hanging wall. The secondary geological hazards, collapses, landslides, debris flows et al., concentrated in the hanging wall of co-seismic fault within 0-20 km from co-seismic fault.

Co-seismic displacements of 2011 Japan Mw9.0 earthquake recorded by far-field GPS stations

Co-seismic displacements of the 2011 Mw9.0 Japan earthquake recorded by GPS stations in China and surrounding areas showed a movement toward the epicenter. The horizontal displacements were up to 1 - 3 cm in northeastern China, 3 -8 mm in the North China, and 2 cm in the Korean peninsula. The vertical movements in China were small uplifts.

Hazard Assessment of Co-seismic Landslides Based on Information Value Method:A Case in 2018 MW6.6 Hokkaido Earthquake,Japan

An MW6.6 earthquake occurred in eastern Hokkaido,Japan on September 6th,2018.Based on the pre-earthquake image from Google Earth and the post-earthquake image from high resolution(3 m)planet satellite,we manually interpret 9293 coseismic landslides and select 7 influencing factors of seismic landslide,such as elevation,slope,slope direction,road distance,flow distance,peak ground acceleration(PGA)and lithology.Then,9293 landslide points are randomly divided into training samples and validation samples with a proportion of 7:3.In detail,the training sample has 6505 landslide points and the validation sample has 2788 landslide points.The hazard risk assessment of seismic landslide is conducted by using the information value method and the study area is further divided into five risk grades,including very low risk area,low risk area,moderate risk area high risk area and very high risk area.The results show that there are 7576 landslides in high risk area and very high risk area,accounting for81.52%of the total landslide number,and the landslide area is 22.93 km^2,accounting for 74.35%of the total area.The hazard zoning is in high accordance with the actual situation.The evaluation results are tested by using the curve of cumulative percentage of hazardous area and cumulative percentage of landslides number.The results show that the success rate of the information value method is 78.50%and the prediction rate is 78.43%.The evaluation results are satisfactory,indicating that the hazard risk assessment results based on information value method may provide scientific reference for landslide hazard risk assessment as well as the disaster prevention and mitigation in the study area.

Multi-scale Decomposition of Co-seismic Deformation from High Resolution DEMs:a Case Study of the 2004 Mid-Niigata Earthquake

Decomposing co-seismic deformation is an immediate need for researchers who are interested in earthquake inversion analysis and geo-hazard mapping. However, conventional InSAR or digital elevation models （DEMs） imagery analyses only provide the displacement in the Line-of-Sight （LOS） direction or elevation changes. The 2004 Mid-Niigata earthquake in Japan provides lessons on how to decompose co-seismic deformation from two sets of DEMs. If three adjacent points undergo a rigid-body-translation movement, their co-seismic deformation can be decomposed by solving simultaneous equations. Although this method has been successfully used to discuss tectonic deformations, the algorithm needed improvement and a more rigorous algorithm, including a new definition of nominal plane, DEMs comparability improvement and matrix condition check is provided. Even with these procedures, the obtained decomposed displacement often showed remarkable scatter prompting the use of the moving average method, which was used to determine both tectonic and localized displacement characteristics. A cut-off window and a pair of band-pass windows were selected according to the regional geology and construction activities to ease the tectonic and localized displacement calculations, respectively. The displacement field of the tectonic scale shows two major clusters of large lateral components, and coincidently major visible landslides were found mostly within them. The localized displacement helps to reveal hidden landslides in the target area. As far as the Kizawa hamlet is concerned, the obtained vectors show down-slope movements, which are consistent with the observed traces of dislocations that were found in the Kizawa tunnel and irrigation wells. The method proposed has great potential to be applied to understanding post-earthquake rehabilitation in other areas.

Numerical simulation of co-seismic deformation of 2011 Japan Mw9.0 earthquake

Co-seismic displacements associated with the Mw9.0 earthquake on March 11, 2011 in Japan are numerically simulated on the basis of a finite-fault dislocation model with PSGRN/PSCMP software. Compared with the inland GPS observation, 90% of the computed eastward, northward and vertical displacements have residuals less than 0.10 m, suggesting that the simulated results can be, to certain extent, used to demon- strate the co-seismic deformation in the near field. In this model, the maximum eastward displacement increa- ses from 6 m along the coast to 30 m near the epicenter, where the maximum southward displacement is 13 m. The three-dimensional display shows that the vertical displacement reaches a maximum uplift of 14.3 m, which is comparable to the tsunami height in the near-trench region. The maximum subsidence is 5.3 m.

Co-seismic deformation derived from GPS observations during April 20th,2013 Lushan Earthquake,Sichuan,China

We process the standard 30 s, static GPS data and the 1 s, high-rate GPS （HRGPS） data provided by the Crustal Movement Observation Network of China with GAMIT/GLOBK software package, and obtain the co- seismic displacements of near field and far field, and the epoch-by-epoch time series of HRGPS during Lushan earthquake. GPS data from about 20 sites in Sichuan province, which located between 40 and 450 km from the epicenter, are analyzed so as to study the characteristics of the static displacements and the dynamic crustal defor- mations, with periods ranging from several minutes to over a month. The result shows that： the static displacements caused by Lushan earthquake are limited to several centi- meters; the nearest station SCTQ at 43 km from the epi- center has the largest static displacement of about 2 cm, while the other stations generally have insignificant dis- placements of less than 5 mm. the stations in the east ofSichuan-Yunnan region shifts 5-10 mm toward the southwest, and the stations in the middle-west of Sichuan Basin moves indistinctively 1-2 mm toward the northwest; station SCTQ has the largest kinematic displacement of about 4 and 3 cm peak-to-peak on the north and east component, respectively, and is much greater than the static permanent displacement; for the stations located at a distance greater than 150 km from the epicenter, the kinematic motions are generally insignificant; exception- ally, station SCNC and station SCSN in central Sichuan Basin have significant kinematic motions although they are more than 200 km away from the epicenter.