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 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.

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.

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）.

Distribution patterns of rock mass displacement in deeply buried areas induced by active fault creep slip at engineering scale

Active fault creep slip induces deformation of rock mass buried deeply in fault zones that significantly affect the operational safety of long linear projects passing through it.Displacement distribution patterns of rock masses in active fault zones which have been investigated previously are the key design basis for such projects.Therefore,a discrete element numerical model with different fault types,slip time,dip angles,and complex geological features was established,and then the creep slip for normal,reverse,and strike-slip faults were simulated to analyze the displacement distribution in the fault rock mass.A disk rotation test system and the corresponding laboratory test method were developed for simulating rock mass displacement induced by creep slippage of faults.A series of rotation tests for softand hard-layered specimens under combined compression and torsional stress were conducted to verify the numerical results and analyze the factors influencing the displacement distribution.An S-shaped displacement distribution independent of fault dip angle was identified corresponding to reverse,normal,and strike-slip faults.The results indicated that the higher the degree of horizontal extrusion,the softer the rock mass at the fault core,and the higher the degree of displacement concentration in the fault core;about 70%of the creep slip displacement occurs within this zone under 100 years of creep slippage.

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 earth- quake, 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 Posi- tioning 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 kin, 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 distri- bution 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 slip distribution of 2009 L'Aquila earthquake derived from InSAR and GPS data

To better understand the mechanism of the Mw6.3 L＇Aquila （Central Italy） earthquake occurred in 2009, global positioning system （GPS） and interferometric synthetic aperture radar （InSAR） data were used to derive the coseismic slip distribution of the earthquake fault. Firstly, based on the homogeneous elastic half-space model, the fault geometric parameters were solved by the genetic algorithm. The best fitting model shows that the fault is a 13.7 km×14.1 km rectangular fault, in 139.3° strike direction and 50.2° southwest-dipping. Secondly, fixing the optimal fault geometric parameters, the fault plane was extended and discretized into 16× 16 patches, each with a size of 1 kmx 1 krn, and the non-uniform slip distribution of the fault was inverted by the steepest descent method with an appropriate smoothing ratio based on the layered crustal structure model. The preferred solution shows that the fault is mainly a normal fault with slight right-lateral strike slip, the maximum slip of 1.01 m is located in the depth of 8.28 km, the average rake is -100.9°, and the total geodetic moment is about 3.34× 1018 N.m （Mw 6.28）. The results are much closer than previous studies in comparison with the seismological estimation. These demonstrate that the coseismic fault slip distribution of the L＇Aauila earthauake inverted by the crustal model considering layered characters is reliable.

Characteristics of InSAR Coseismic Deformation and Slip Distribution of the Akto M_S6. 7 Earthquake,Xinjiang

The Akto MS6. 7 earthquake occurred near the western end of the Muji fault basin in the top of the Pamir syntaxis. The main shock of this earthquake is complicated and the focal mechanism solutions based on the seismic wave inversions are different. Based on the Sentinel-1 SAR data,the coseismal deformation field of the earthquake is obtained by In SAR technique. Based on the elastic half-space dislocation model,the geometrical parameters and the slip distribution model are determined by nonlinear and linear inversion algorithms. The results show that the distributed slip model can well explain the coseismic deformation field. The earthquake includes at least two rupture events,which are located at 7 km(74. 11°E,39. 25°N)and 33 km(74. 49°E,39. 16°N)east from the epicenter according to the CENC. The deformation field caused by the earthquake shows a symmetry distribution,with the maximum LOS deformation of 20 cm. The main seismic slip is concentrated in the 0-20 km depth,and the maximum slip is 0. 84 m. The seismic fault is the Muji fault,and this earthquake indicates that the northeastward push of the Indian plate is enhanced.

Resolving co- and early post-seismic slip variations of the 2021 MW 7.4 Madoi earthquake in east Bayan Har block with a block-wide distributed deformation mode from satellite synthetic aperture radar data

On 21 May 2021(UTC),an MW 7.4 earthquake jolted the east Bayan Har block in the Tibetan Plateau.The earthquake received widespread attention as it is the largest event in the Tibetan Plateau and its surroundings since the 2008 Wenchuan earthquake,and especially in proximity to the seismic gaps on the east Kunlun fault.Here we use satellite interferometric synthetic aperture radar data and subpixel offset observations along the range directions to characterize the coseismic deformation of the earthquake.Range offset displacements depict clear surface ruptures with a total length of~170 km involving two possible activated fault segments in the earthquake.Coseismic modeling results indicate that the earthquake was dominated by left-lateral strike-slip motions of up to 7 m within the top 12 km of the crust.The well-resolved slip variations are characterized by five major slip patches along strike and 64%of shallow slip deficit,suggesting a young seismogenic structure.Spatial-temporal changes of the postseismic deformation are mapped from early 6-day and 24-day InSAR observations,and are well explained by time-dependent afterslip models.Analysis of Global Navigation Satellite System(GNSS)velocity profiles and strain rates suggests that the eastward extrusion of plateau is diffusely distributed across the east Bayan Har block,but exhibits significant lateral heterogeneities,as evidenced by magnetotelluric observations.The block-wide distributed deformation of the east Bayan Har block along with the significant co-and post-seismic stress loadings from the Madoi earthquake imply high seismic risks along regional faults,especially the Tuosuo Lake and Maqên-Maqu segments of the Kunlun fault that are known as seismic gaps.

Torque Distribution of Electric Vehicle with Four In-Wheel Motors Based on Road Adhesion Margin

With the worsening of energy crisis and environmental pollution,electric vehicles with four in?wheel motors have been paid more and more attention. The main research subject is how to reasonably distribute the driving torque of each wheel. Considering the longitudinal motion,lateral motion,yaw movement and rotation of the four wheels,the tire model and the seven DOF dynamic model of the vehicle are established in this paper. Then,the torque distribution method is proposed based on road adhesion margin,which can be divided into anti ? slip control layer and torque distribution layer. The anti?slip control layer is built based on sliding mode variable structure control,whose main function is to avoid the excessive slip of wheels caused by road conditions. The torque distribution layer is responsible for selecting the torque distribution method based on road adhesion margin. The simulation results show that the proposed torque distribution method can ensure the vehicle quickly adapt to current road adhesion conditions,and improve the handling stability and dynamic performance of the vehicle in the driving process.

Study on Running Performance of Flexible Wheels for Lunar Rovers： Relation between Stress Distribution of Flexible Wheel and Running Performance

On lunar exploration missions, the rovers which can move and explore directly are considered by various agency like NASA （National Aeronautics and Space Administration）, JAXA （Japan Aerospace Exploration Agency）, ESA （European Space Agency）. Lunar rovers are required to move on rough terrains such as craters and rear cliffs where it is scientifically very important to explore. However, there is a problem that the rovers have possibility of stack because of the lunar surface is covered with loose soil named Regolith. Therefore, this paper investigates a mechanism of kinetic behavior between the wheels of the exploration rovers and loose soil. And then, this paper proposed a flexible wheel to solve like that problems. The flexible wheel has the surface which can be changed flexibly toward rough terrain. Running experiments on loose soil which imitated regolith were carried out to observe the traversability of the flexible wheel using slip ratio. Traversality of flexible wheel was better than the circular rigid wheel. The authors believe that stress distribution is important. The stress distribution of the flexible wheels is horizontally long and stress value is small. However, the stress distribution can be changed by loaded more weight. Therefore, the relationship between the stress and the running performance was considered using this differential stress distribution. In experiments, the authors used the flexible wheel with simple structure （3 limbs）. From these considerations, the relationship between the stress of the flexible wheel and the running performance was described.

GNSS约束的2022年泸定M6.8地震滑动分布及同震应力变化

2022年9月5日四川泸定县发生M6.8级强震,震中位于鲜水河断裂东南末端,野外地质调查初步结果显示本次地震并未发现明显的地表破裂迹象.本文基于震后科考GNSS流动观测和震中周边连续站观测资料,解算并提取震中90 km范围内31个测站的静态水平向同震位移.结果显示:GNSS同震形变场空间分布呈现明显的左旋走滑特征,观测到的最大水平向同震位移达23 cm,震中距50 km范围内同震位移量普遍大于1 cm.基于GNSS观测资料反演得到的同震滑动分布显示泸定地震地表破裂主要集中在磨西至田湾之间,主破裂深度2~8 km,最大滑动量~1.96 m,地震矩9.25×10^(18)N·m,对应矩震级M_W6.6.静态库仑应力结果显示本次地震增强了震源破裂区周边活动断层的库仑应力,并触发了大量余震,余震主要发生在库仑应力增强区域.结合震间闭锁分布、历史地震及库仑应力变化,我们认为未来需要密切关注与磨西断裂交接的安宁河、大凉山断裂以及康定—磨西段的地震危险性.

2020年西藏定日M_(W)5.6地震震源参数估计和应力触发研究

2020年3月20日青藏高原西南缘定日县发生M_(W)5.6地震,距2015年尼泊尔M_(W)7.9地震~250 km.尼泊尔地震,尤其是震后余滑是否触发了此次定日地震还有待研究.本文联合合成孔径雷达和区域地震波资料研究定日地震的破裂特征.首先利用近场形变和宽频带地震波资料,通过贝叶斯自举优化算法揭示定日地震的均匀滑动模型;然后在此基础上构建断层几何模型并反演震源滑动分布.研究发现定日地震的发震断层走向~334°,倾角~51°.破裂主要集中在约2.0~5.5 km深度范围内.破裂范围~5.6 km×4.4 km,释放总的地震矩~3.33×10^(17)N·m.最大滑动量~1.27 m,发生在3.786 km深度.破裂以正断滑动为主兼少许右旋走滑分量,同区域历史地震表现出相似的破裂机制,表明印度板块向北东方向挤压欧亚板块,在藏南地区产生了近东西向的张应力.库仑应力变化研究表明,尼泊尔M_(W)7.9地震主余震和定日地区四次历史地震共同触发了2020年定日M_(W)5.6地震,其中尼泊尔地震震后2年的余滑引起的库仑应力变化占库仑应力增加总量的~40%,震后余滑在未来地震危险性评估中发挥的作用不容忽视.

锚杆岩体界面载荷传递规律及锚固长度设计

为揭示巷道围岩锚固界面脱粘失效机理,量化锚杆支护设计参数,采用三线性粘结滑移模型进行理论分析,对轴向载荷作用下锚固脱粘失效全过程中,锚固段界面剪应力、锚杆轴力分布演化规律以及界面极限锚固力进行研究,根据锚固段长度不同,得到两种界面剪应力分布演化类型。研究结果表明:当锚固长度较短时,界面剪应力存在全长软化阶段;当锚固长度较长时,界面剪应力存在弹性-软化-滑移三段共存阶段。锚固粘结界面弹性段、软化段、摩擦段内的剪应力分别呈现双曲余弦函数衰减分布、余弦函数上升分布、均匀分布规律,锚杆轴力随界面剪应力分布演化呈现多种形态的衰减分布规律。根据锚固界面模型解析计算获得极限锚固力,当不考虑脱粘摩擦力时,极限锚固力随锚固长度的增加趋近于某一固定值;当考虑脱粘摩擦力时,增加锚固长度能够持续提高锚固界面安全系数。研究成果可为锚固机制分析、锚杆支护参数设计提供理论参考。

2022年泸定M_(S)6.8地震震前区域变形背景及同震形变特征

文中利用Sentinel-1 SAR影像获取了2022年泸定6.8级地震的震间形变场(2014—2020年)和同震形变场,估算了泸定地震的震间断层滑动速率、闭锁深度和同震滑动分布,分析了此次地震对周边断裂的影响,并进一步探讨泸定地震的发震构造和磨西断裂及其周边断裂的未来地震发生趋势。震间InSAR形变场及InSAR和GNSS融合的三维形变场结果显示了2022年泸定地震发震断层的长期滑动特征,滑动速率为(5.9±1.8)mm/a,主要表现为剪切变形明显且应变集中的特征,震前磨西断裂存在明显闭锁,具有强震发生背景。同震滑动分布反演结果表明,此次泸定地震是以左旋走滑为主的高倾角走滑地震,最大滑动量可达1.71m,最大滑动位于10km深处。基于文中获得的同震滑动分布计算了同震位错效应引起的断层面上的库仑应力变化,结果显示磨西断裂南段的地震破裂段形成应力影区,而在此次地震中未破裂的北段的库仑应力显著升高,同时,折多塘断裂南东部分、安宁河断裂带石棉—冕宁段北西部分、大凉山断裂竹马段南东部分和公益海段南东部分的断层面库仑应力显著增加。2022年泸定地震的发震构造为鲜水河断裂带的磨西断裂,此次泸定地震的发生并没有完全降低鲜水河断裂带的地震风险,仍需重点关注磨西断裂未破裂段及其周边具有大地震发生背景的强震破裂空段。

四轮独立电驱动越野车驱动力优化控制

为了充分利用四轮独立电驱动越野车各轮转矩独立可控的优势,提高越野车的牵引力和越野能力,提出了一种驱动转矩协调控制策略.根据越野车前后轴载荷预分配驱动电动机转矩,利用基于车轮滑转率-路面附着系数的极值寻求算法实现路面最优滑转率估计.采用比例-积分-微分(proportional-integral-differential, PID)控制-滑模控制(sliding mode control, SMC)和基于状态的驱动力再分配方法实现各轮驱动转矩的协调分配,抑制越野车在恶劣路面条件下的车轮打滑.通过CarSim/MATLAB联合仿真以及硬件在环(hardware-in-the-loop, HIL)测试,进行了爬陡坡、单一附着路面、对开路面及连续起伏路面工况的试验验证.结果表明:提出的控制策略能根据实际工况分配驱动轮的转矩,降低驱动轮滑转率,实现整车驱动力的优化.

2022年门源M_(S)6.9地震的同震滑动分布:联合InSAR、GPS和地表位错的贝叶斯建模

2022年1月8日青海省门源M_(S)6.9地震发生于托莱山断裂东段和冷龙岭断裂西段的交汇部位,其运动学特征关系到青藏高原北缘的动力学行为,同时也是区域地震危险性评估中的基础数据.本研究针对该地震,综合InSAR和GPS观测的同震位移场,以及野外地质调查的同震位错数据,采用贝叶斯反演方法,构建同震滑动分布,断层模型采用三角网格更好地刻画断层的不规则几何形态.结果显示,本次左旋走滑型地震的主要滑动发生在冷龙岭断裂西段,深度范围约2~6 km,计算的地震矩为0.95×10^(19)N·m,对应矩震级M_(W)6.65.利用多数据的联合反演解析出,沿走向存在两个滑动集中区,东侧的最大滑动达约4.8 m.震间亏损能量和库仑应力分析表明,托莱山断裂和冷龙岭断裂仍存在未来发震潜能;该地震造成的应力变化,对长期地震空区金强河—毛毛山段落有加载作用,因此,该段未来地震危险性值得关注.

2024年1月23日M_(W)7.0乌什地震InSAR同震形变场和断层滑动分布

2024年1月23日,新疆维吾尔自治区乌什县发生了M_(W)7.0地震,本次地震是南天山断裂带近一个世纪以来的最大地震事件.为确定乌什地震的发震断层位置及其产状,本研究基于Sentinel-1A数据和差分干涉测量、像素偏移技术获得乌什地震升降轨同震形变场,并通过贝叶斯非线性反演方法估计了发震断层的几何参数,获得乌什地震同震滑动分布.研究结果显示,升轨最大视线向抬升形变达80 cm,最大视线向沉降约16 cm.视线向同震形变与像素偏移量结果表明乌什地震显著的垂向运动特征,而升降轨像素偏移方位向观测的反向运动特征表明乌什地震发震断层明显的左旋走滑运动分量.反演结果显示乌什地震发震断层为东北东—西南西走向,北西倾的断裂产状,最优断层倾角约67°,平均滑动角约60°,表明乌什地震同震滑动以逆冲运动为主,兼具少量左旋走滑分量.本文获取的乌什地震变形方式和断层参数与天山南缘盆—山边界的迈丹—沙依拉姆断裂的运动学和几何学相符.乌什地震的高倾角断层走滑-逆冲活动表明南天山与塔里木盆地之间的斜向汇聚作用,天山地区的地壳缩短由山前新生的逆冲推覆带和造山区的高角度逆断层共同吸收.

基于InSAR形变的滑面摩擦残余估计与滑坡危险性评价

传统滑坡灾害危险性评价往往依赖于坡表位移或形变速率等间接特征开展,但较少直接利用滑面自身的物理特征评估滑坡的危险性。在此,本文发展了一种依据滑面摩擦残余的斜坡危险性评价理论与方法,其基本流程为:首先利用InSAR技术提取滑坡时序变形数据,其次在坡表形变的约束下反演获得斜坡沿滑面的滑移时空演化特征;进一步依据速度摩擦衰减定律计算获得坡体的摩擦系数残余及稳定摩擦系数,并据此分析滑面稳定摩擦系数的时空变化,进而定量评估滑坡的危险性。本文选取西藏山南市隆子县庞村滑坡作为研究对象,利用发展理论与方法获得庞村滑坡的多期滑面稳定摩擦系数,结果显示滑面稳定摩擦系数虽然会随季节出现一定的周期性波动,但总体上呈现显著的下降趋势,并由滑坡后缘向滑坡前缘发展,同时初步推测在滑面稳定摩擦系数小于0.4时,滑坡可能会进入临滑失稳状态。本文提出的利用InSAR形变的滑面摩擦残余与滑坡危险性评价方法,可相对便捷、迅速地实现对滑坡状态的定量评价,且成本较低,可为广域滑坡灾害的监测预警及风险性评估提供有意义的科学参考。