THE COSEISMIC SOURCE SLIP AND COULOMB STRESS TRIGGERING OF 2015 NEPAL GORKHA MW7.9 AND KODARIMW7.3 EARTHQUAKE BASED ON InSAR MEASUREMENTS

According to the structure of the Himalayan orogenic belt, a low-angle antilistric thrust-slip fault model is used to simulate the ramp on the rupture portion of the Main Himalayan Fault. Based on descending Alos -2 and Sentinal -1 data, we invert for the coseismic slip models of the Gorkha earthquake and its largest aftershock, Kodari earthquake. In contrast to the inversion using Alos -2 or Sentinal -1 separately, the joint inversion of both data sets has stronger constraint for the deep slip and can obtain more details in Gorkha earthquake. The rupture depth obtained by joint inversion can be as deep as 24km underground, cutting across the locking line to the transition of locked and the creeping zone. The largest slip is as large as 4.5m appearing 17km underground and the dip angle is between 3°and 10°. Gorkha and Kodari earthquakes have the similar focal mechanisms, both of which are mainly thrusting, and yet some right-lateral slip component in Gorkha earthquake. The inversion results reveal that slip models of the Nepal mainshock and its largest aftershock are complementary in space and the Kodari earthquake occurs in the gaps of slip in Gorkha earthquake. The epicenter of the Kodari earthquake is just right in the transitive zone of the positive and negative Coulomb stress change, where the Coulomb stress change can reach 0.4MPa. We thus argue that Kodari earthquake has been triggered by the Gorkha earthquake.     

THE STATIC STRESS TRIGGERING INFLUENCES OF THE 2015 MW6.4 PISHAN,XINJIANG EARTHQUAKE ON THE NEIGHBORING AREAS

Based on the rupture models of the 2015 Pishan M_W6.4 earthquake and half space homogeneous elastic model, the Coulomb stress changes generated by the earthquake are calculated on the active faults near the earthquake region. The horizontal stress changes and the displacement field are estimated on the area around the epicenter. Results show that:(1)The Coulomb stress is decreased in the west of the western Kunlun frontal thrust fault(9.5×10³Pa), and increased in the east of the western Kunlun frontal thrust fault and the middle of the Kangxiwa faults. More attention should be taken to the seismic rick of the east of the western Kunlun frontal thrust fault; (2)Based on the analysis on the location of the aftershocks, it is found that most of the aftershocks are triggered by the earthquake. In the region of increased Coulomb attraction, the aftershock distribution is more intensive, and in the area of the Coulomb stress reduction, the distribution of aftershocks is relatively sparse; (3)The horizontal area stress increases in the north and south of the earthquake(most part of the Qaidam Basin and the northwest of the Qinghai-Tibet plateau), and decreases in the east and west of the earthquake(northern part of the Qinghai-Tibet plateau and eastern part of the Pamir Mountains). In the epicenter area, the principal compressive stress presents nearly NS direction and the principal extensional stress presents nearly EW direction. The principal compressive stress shows an outward radiation pattern centered on the epicenter with the principal extensional stress along the direction of concentric circles. The principal compressive stress presents NW direction to the west of the epicenter, and NE to the east of the epicenter. With the increase of radius, the stress level gradually decays with 10⁷Pa in the epicenter and hundreds Pa in the Maidan Fault which is in the north of the Qaidam Basin.     

2010年1月12日海地 Mw7.0地震InSAR同震形变观测及同震滑动分布反演

2010年1月12日GMT时间21时53分,在海地境内(72.57°W,18.44°N)发生了Mw7.0地震.文中利用干涉合成孔径雷达(InSAR)方法获得了覆盖整个震区的高精度形变观测资料,用以研究该地震的发震机理.采用ALOs PALSAR数据,分析了轨道、大气等误差源对干涉信号的影响,最终获得了雷达视线向(LOS...     

COSEISMIC DISPLACEMENT AND FAULT SLIP OF THE MW6.1 NAPA EARTHQUAKE IN AMERICA REVEALED BY SENTINEL-1A INSAR DATA

We achieved the coseismic displacements of the Napa M_W6.1 earthquake located in California US occurring on 24 August 2014 by using InSAR data from the newly launched ESA's Sentinel-1A satellite. The 30m×30m ASTER GDEM was used to remove the terrain effect, and phase unwrapping method of branch-cut algorithm was adopted. In order to obtain a better coseismic displacement field, we also tested 90m×90m SRTM data to remove the terrain effect and Minimum Cost Flow algorithm to unwrap the phase. Results showed that the earthquake caused a significant ground displacement with maximum uplift and subsidence of 0.1m and -0.09m in the satellite light of sight(LOS). Based on the Sentinel-1A dataset and sensitivity based iterative fitting(SBIF) method of restrictive least-squares algorithm, we obtained coseismic fault slip distribution and part of the earthquake source parameters. Inversion results show that the strike angle is 341.3°, the dip angle is 80°, rupture is given right-lateral fault, average rake angle is -176.38°, and the maximum slip is ~0.8m at a depth of 4.43km. The accumulative seismic moment is up to 1.6×10¹⁸N·m, equivalent to a magnitude of M_W6.14.     

DATABASE OF LANDSLIDES TRIGGERED BY 2015 GORKHA(NEPAL) MW7.8 EARTHQUAKE

In this study, a detailed database of landslides triggered by the 25 April 2015 Gorkha (Nepal)M_W7.8 earthquake is constructed based on visual interpretation of pre- and post-earthquake high-resolution satellite images and field reconnaissance. Results show the earthquake triggered at least 47 200 landslides, which have a NWW direction spatial distribution, similar with the location and strike of the seismogenic fault. The landslides are of a total area about 110km² and an oval distribution area about 35 700km². On the basis of a scale relationship between landslide area (A)and volume (V), V=1.314 7×A^{1.208 5}, the total volume of the coseismic landslides is estimated to be about 9.64×10⁸m³. In the oval landslide distribution area, the landslide number density, area density, and volume density were calculated and the results are 1.32km^-2, 0.31%, and 0.027m, respectively. This study provides a detailed and objective inventory of landslides triggered by the Gorkha earthquake, which provides very important and essential basic data for study of mechanics of coseismic landslides, spatial pattern, distribution law, and hazard assessment. In addition, the landslide database related to an individual earthquake also provides an important earthquake case in a subduction zone for studying landslides related to multiple earthquakes from a global perspective.     

GPS AND LEVELING CONSTRAINED CO-SEISMIC SOURCE AND SLIP DISTRIBUTION OF THE LUSHAN MS7.0 EARTHQUAKE ON 20 APRIL 2013

Co-seismic deformation results calculated by the observations of GPS continuous sites and precise leveling are used as three-dimensional constraints for the deformation field of Lushan M_S7.0 earthquake on 20 April 2013.The inversion of seismic source parameters are calculated based on the consideration of different value-taking schemes and fault models. Seven data type selection schemes, three fault models and two data coverage selection schemes are designed in order to discuss the effect of data selection and fault model selection on the inversion results. The results show that the fault model using blind reverse fault for the inversion is superior to the model using the fault that ruptures from its upper boundary to the earth'surface, which may indicate that the Lushan earthquake fault is most likely a blind reverse fault; there are no obvious differences in the inversion results between the blind listric reverse fault models and the rectangle blind reverse fault models. The best inversion result of Lushan earthquake seismic moment is M_W6.7.And it also shows that the distribution of dislocations on the fault plane is concentrated in the range of 30km×30km, the northern flank of the seismogenic fault of Lushun earthquake is of dextral faulting and the southern flank of the fault is of sinistral faulting, the sinistral component is larger that dextral component, showing a wedge deformation mode.     

COSEISMIC DISPLACEMENT FIELD OF THE WENCHUAN EARTHQUAKE DERIVED FROM STRONG MOTION RECORDS AND APPLICATION IN SLIP INVERSION

The development of high-rate GNSS seismology and seismic observation methods has provided technical support for acquiring the near-field real-time displacement time series during earthquake. But in practice, the limited number of GNSS continuous stations hardly meets the requirement of near-field quasi-real-time coseismic displacement observation, while the macroseismographs could be an important complement. Compared with high-rate GNSS, macroseismograph has better sensitivity, higher resolution(100~200Hz)and larger dynamic range, and the most importantly, lower cost. However, baseline drift exists in strong-motion data, which limits its widespread use. This paper aims to prove the feasibility and reliability of strong motion data in acquiring seismic displacement sequences, as a supplement to high-rate GNSS. In this study, we have analyzed the strong-motion data of Wenchuan M_S8.0 earthquake in Longmenshan fault zone, based on the automatic scheme for empirical baseline correction proposed by Wang et al., which fits the uncorrected displacement by polynomial to obtain the fitting parameters, and then the baseline correction is completed in the velocity sequence. Through correction processing and quadratic integration, the static coseismic displacement field and displacement time series are obtained. Comparison of the displacement time series from the strong motions with the result of high-rate GPS shows a good coincidence. We have worked out the coseismic displacement field in the large area of Wenchuan earthquake using GPS data and strong motion data. The coseismic displacement fields calculated from GPS and strong motions are consistent with each other in terms of magnitude, direction and distribution patterns. High-precision coseismic deformation can provide better data constraint for fault slip inversion. To verify the influence of strong-motion data on slip distribution in Wenchuan earthquake, we used strong motion, GPS and InSAR data to estimate the stress drop, moment magnitude and coseismic slip model, and our results agreed with those of the previous studies. In addition, the inversion results of different data are different and complementary to some extent. The use of strong-motion data supplements the slip of the fault in the 180km segment and the 270~300km segment, thus making the inversion results of fault slip more comprehensive. From this result, we can draw the following conclusions:1)Based on the robust baseline correction method, the use of strong motion data, as an important complement to high-rate GNSS, can obtain reliable surface displacement after the earthquake. 2)The strong motion data provide an effective method to study the coseismic displacement sequence, the surface rupture process and quick seismogenic parameters acquisition. 3)The combination of multiple data can significantly improve the data coverage and give play to the advantages of different data. Therefore, it is suggested to combine multiple data(GPS, strong motion, InSAR, etc.)for joint inversion to improve the stability of fault slip model.     

USE OF SEISMIC WAVEFORMS AND INSAR DATA FOR DETERMINATION OF THE SEISMOTECTONICS OF THE MAINLING MS6.9 EARTHQUAKE ON NOV.18, 2017

On November 18, 2017, a M_S6.9 earthquake struck Mainling County, Tibet, with a depth of 10km. The earthquake occurred at the eastern Himalaya syntaxis. The Namche Barwan moved northward relative to the Himalayan terrane and was subducted deeply beneath the Lhasa terrane, forming the eastern syntaxis after the collision of the Indian plate and Asian plates. Firstly, this paper uses the far and near field broadband seismic waveform for joint inversion (CAPJoint method)of the earthquake focal mechanism. Two groups of nodal planes are obtained after 1000 times Bootstrap test. The strike, dip and rake of the best solution are calculated to be 302°, 76° and 84° (the nodal plane Ⅰ)and 138°, 27° and 104° (the nodal plane Ⅱ), respectively. This event was captured by interferometric synthetic aperture radar (InSAR)measurements from the Sentinel-1A radar satellite, which provide the opportunity to determine the fault plane, as well as the co-seismic slip distribution, and assess the seismic hazards. The overall trend of the deformation field revealed by InSAR is consistent with the GPS displacement field released by the Gan Wei-Jun's team. Geodesy (InSAR and GPS)observation of the earthquake deformation field shows the northeastern side of the epicenter uplifting and the southwestern side sinking. According to geodetic measurements and the thrust characteristics of fault deformation field, we speculate that the nodal plane Ⅰ is the true rupture plane. Secondly, based on the focal mechanism, we use InSAR data as the constraint to invert for the fine slip distribution on the fault plane. Our best model suggests that the seismogenic fault is a NW-SE striking thrust fault with a high angle. Combined with the slip distribution and aftershocks, we suggest that the earthquake is a high-angle thrust event, which is caused by the NE-dipping thrust beneath the Namche Barwa syntaxis subducted deeply beneath the Lhasa terrane.     

COSEISMIC SURFACE RUPTURES AND SEISMOGENIC MUJI FAULT OF THE 25 NOVEMBER 2016 ARKETAO MW6.6 EARTHQUAKE IN NORTHERN PAMIR

The M_W6.6 Arketao earthquake,which occurred at 14:24:30 UTC 25 November 2016 was the largest earthquake to strike the sparsely inhabited Muji Basin of the Kongur extension system in the eastern Pamir since the M 7 1895 Tashkurgan earthquake.The preliminary field work,sentinel-1A radar interferometry,and relocated hypocenters of earthquake sequences show that the earthquake consists of at least two sub-events and ruptured at least 77km long of the active Muji dextral-slip fault,and the rupture from this right-lateral earthquake propagated mostly unilaterally to the east and up-dip.Tectonic surface rupture with dextral slip of up to 20cm was observed on two tens-meter long segments near the CENC epicenter and 32.6km to the east along the Muji Fault,the later was along a previously existing strand of the Holocene Muji fault scarps.Focal mechanisms are consistent with right-lateral motion along a plane striking 107°,dipping 76° to the south,with a rake of 174°.This plane is compatible with the observed tectonic surface rupture.More than 388 aftershocks were detected and located using a double-difference technique.The mainshock is relocated at the Muji Fault with a depth of 9.3km.The relocated hypocenters of the 2016 Arketao earthquake sequence showed a more than 85km long,less than 8km wide,and 5~13km deep,NWW trending streak of seismicity to the south of the Muji Fault.The focal mechanism and mapping of the surface rupture helped to document the south-dipping fault plane of the mainshock.The listric Muji Fault is outlined by the well-resolved south-dipping streak of seismicity.The 2016 Arketao M_W6.6 and 2015 Murghob M_W7.2 earthquakes highlight the importance role of strike-slip faulting in accommodating both east-west extensional and north-south compressional forces in the Pamir interior,and demonstrate that the present-day stress and deformation patterns in the northern Pamir plateau are dominant by east-west extension in the shallow upper crust.     

2001年昆仑山口西MS8.1地震地表同震位移分布特征

沿长约 4 2 6km的 2 0 0 1年昆仑山口西MS8 1地震地表破裂带共获得 2 91个点的地表同震水平左旋位移数据 ,并在其中 1 1 1个点获得了垂直位移数据。该地震总体以左旋水平位移为主 ,兼具一定的垂直位移。最大地表左旋水平位移值可达 6 4m ,平均水平位移约为 2 7m ,绝大多数测点的垂直位移均 <1m。地表水平位移沿主破裂带走向位移梯度变化于 1 0 - 1～ 1 0 - 4之间 ,这一起伏变化可能起因于野外测量误差、沿主破裂带岩性或松散沉积物厚度的变化、地表破裂带几何结构的不均匀性、地表破裂走向的变化、不同破裂段在昆仑山口西 8 1级地震之前的地震中滑动量的起伏变化 ,以及大量非脆性变形、次级破裂的存在等。水平位移沿主破裂带的长波长 (数十公里至数百公里 )起伏变化较有规律 ,在布喀达坂峰以东表现为分别以 5个水平位移峰值为中心而有规律地起伏变化。这5个位移峰值分别对应于不同的次级地震地表破裂段。各破裂段水平位移峰值均向阶区或拐点逐渐衰减 ,不同地表破裂段位移峰值向两侧衰减的速率是不同的 ,这种位移梯度的不对称分布可能指示了地震破裂的扩展方向。上述位移分布特征真实地反映了地表可见脆     

用连续GPS与远震体波联合反演2015年尼泊尔中部M_S8.1地震破裂过程

由于印度-欧亚板块碰撞,位于板块边界带的喜马拉雅地区大震频繁,但对其活动性的认识仍十分有限.2015年4月25日尼泊尔中东部地区时隔80年再次发生8级地震,为研究板缘地震提供了一次难得机遇.本文用西藏和尼泊尔的GPS连续观测数据和全球分布的远震地震波记录联合反演此次特大地震的破裂过程,结果显示此次地震发生在印度板块与青藏高原接触边界面——喜马拉雅主滑脱断层上.北倾11°、近东西(295°)走向的断层面破裂约100km长(博卡拉到加德满都),130km宽(从加德满都深入我国西藏吉隆县),破裂以逆冲滑动为主,平均幅度达到2.4m,释放的地震矩高达9.4×1020 N·m.反演结果还显示,震源体主要破裂分布深度范围为5~25km,应无地表破裂,属于一次盲地震.基于GPS资料推测的地壳现今运动速率及1833年地震的震源位置,我们推测地震在此次地震破裂区域复发的周期可能为150~200a,而极震区以南的深部滑脱断层仍保持闭锁,未来仍有导致灾害性大震的可能性.     

GPS和InSAR同震形变约束的尼泊尔M_W7.9和M_W7.3地震破裂滑动分布

2015年4月25日尼泊尔爆发MW7.9地震,继而引发5月12日MW7.3级余震,GPS、InSAR监测到震源区及周边大范围同震形变.本文以国内外的GPS和InSAR同震形变为约束,考虑喜马拉雅断裂带岩石圈垂向分层和横向差异的影响,反演主喜马拉雅逆冲断裂在这次主震和余震中破裂面形状和滑动分布.结果显示,主震从USGS确定的震中位置向东偏南延伸100km以上,破裂地面迹线与主前缘逆冲断裂迹线基本一致.破裂面倾角约7°~11°,大部分破裂集中在深度8~20km,同余震分布深度一致.主震最大滑动量约6.0~6.6m,位于14km深处.余震破裂集中在震中附近30km范围内,填补了主震东部破裂空区,最大滑动约3.6~4.6 m,位于13km深.深度20km以下基本没有破裂.地壳介质不均匀性对破裂滑动分布的影响较大,介质不均匀模型的观测值不符值比各向同性弹性半空间模型降低10%以上.本文地震破裂模型特征与地震反射剖面、以及根据震间期大地测量数据反演的喜马拉雅深部蠕滑剖面极其相似.跨喜马拉雅断裂剖面的震间形变量与地震破裂滑移量直接相关.以此推算,尼泊尔中部大震原地复发周期在300年以上.     

包头6.4级地震前后小震活动空间分布特征及其破裂方式研究

对1992 年至1996 年4 月发生在40°00′～41°10′N,108°20′110°30′E 范围内的51 次ML≥2 .0 地震空间位置与包头6 .4 级主震及其余震空间位置的相互关系研究表明,前期地震的活动展布方向与包头6 .4 级主震震源机制解主压应力方向一致,震源深度为9 ～34 km ,78% 的地震活动在20 km 以浅的地层内,前期小震活动使包头—呼和浩特地区形成了空区.包头地震发生在此空区内.包头6 .4 级主震震源较深,余震的震源都比主震的浅.在包头6 .4 级地震发生的前1 ～2 a 时间内,围绕主震空间位置发生了一系列小震,这些小震活动可能对包头地震的成因有着较大影响作用.笔者认为包头6 .4 级地震是脆性破裂,属左旋走滑兼具正断错动.     

THE SOURCE PARAMETERS AND RUPTURE PROCESS OF THE MW7.0 EARTHQUAKE IN ALASKA,USA ON DECEMBER 1, 2018

A magnitude M_W7.0 earthquake struck north of Anchorage, Alaska, USA on 1 December 2018. This earthquake occurred in the Alaska-Aleutian subduction zone, on a fault within the subducting Pacific slab rather than on the shallower boundary between the Pacific and North American plates. In order to better understand the earthquake source characteristics and slip distribution of source rupture process as well as to explore the effect of tectonic environment on dynamic triggering of earthquake, the faulting geometry, slip distribution, seismic moment, source time function are estimated from broadband waveforms downloaded from IRIS Data Management Center. We use the regional broadband waveforms to infer the source parameters with ISOLA package and the teleseismic body wave recorded by stations of the Global Seismic Network is employed to conduct slip distribution inversion with iterative deconvolution method. The focal mechanism solution indicates that the Alaska earthquake occurred as the result of tensile-type normal faulting, the estimated centroid depth from waveform inversion shows that the earthquake occurred at the depth of 56.5km, and the centroid location is 10km far away in northeast direction relative to the location of initial epicenter. We use the aftershock distribution to constrain the fault-plane strike of a normal fault to set up the finite fault model, the finite fault inversion shows that the earthquake slip distribution is concentrated mainly on a rectangular area with 30km×20km, and the maximum slip is up to 3.6m. In addition, the slip distribution shows an asymmetrical distribution and the range of possible rupture direction, the direction of rupture extends to the northeast direction, which is same as that of aftershock distribution for a period of ten days after the mainshock. It is interesting to note that a seismic gap appears in the southwest of the seismogenic fault, we initially determined that the earthquake was a typical normal fault-type earthquake that occurred in the back-arc extensional environment of the subduction collision zone between the Pacific plate and the North American plate, this earthquake was not related to tectonic movement of faults near the Earth's surface. Due to the influence of high temperature and pressure during the subduction of the Pacific plate toward to the north, the subduction angle of the Pacific plate becomes steep, causing consequently the backward bending deformation, thus forming to a tensile environment at the trailing edge of the collision zone and generating the M_W7.0 earthquake in Alaska.     

2013年芦山MS7.0地震和2008年汶川MS8.0地震井水位同震变化的比较分析

以芦山、汶川地震震中为中心,选择了8个井点,对2个地震引起的井水位同震变化进行了比较,重点分析了3口同震变化反向的井点,即重庆的荣昌华江井、北碚柳荫井和四川的泸沽湖井。为使对比研究结果可靠,进一步收集了2011年3月11日日本9.0级地震、2012年4月11日苏门答腊北部海域8.6级地震引起的井水位同震变化资料。对汶川、芦山地震引起的同震体应变的计算结果显示：汶川地震时,荣昌华江井、北碚柳荫井和泸沽湖井水位同震变化与井点位置处的同震体应变一致,即位于同震体应变压缩区的井水位上升,位于体应变膨胀区的井水位下降;芦山地震时,3口井的水位同震变化与同震体应变不一致,表现出与日本、苏门答腊地震时相似的同震变化特征。     

2012年6月30日新疆新源、和静交界M_S6.6地震的破裂过程

利用IRIS全球台网记录到的远场波形对2012年新源、和静交界MS6.6地震的破裂过程进行了反演,同时,结合发震构造喀什河断裂东段的构造特征和MS6.6地震前天山中段区域应力场状态,对MS6.6地震的可能发生机理进行了讨论。结果表明:1)MS6.6地震破裂持续时间约35s,地震能量主要在前16s内释放。整个破裂过程由2次子事件组成,第1次破裂强度大于第2次。此次地震破裂过程相对简单,具有双侧破裂特征,最大滑移量为45cm。MS6.6地震初始破裂点位于高滑动量区域的边缘,深部主体破裂区以逆冲兼右旋走滑错动为主,浅部以右旋走滑错动为主;2)MS6.6地震的3级以上余震主要分布在主震破裂大滑动量区域的外围或滑动量变化梯度较大的区域,主震发生后短期内余震的震源机制特征与主震破裂面在断层面上产生的滑动矢量分布较为一致;3)新源、和静交界MS6.6地震前,新疆地区逆冲类型的中强地震明显增多,尤其是天山中段,且这些地震的震源机制与区域应力场的状态较为一致。显示了天山中段受构造应力场控制作用增强,进而有利于具有右旋扭错性质的喀什河断裂东段发生右旋走滑运动,这可能是6.6级地震发生的直接原因。     

利用InSAR研究乌恰地震同震形变

2008年10月5日在新疆乌恰地区39.50°N、73.64°E处发生MW6.7地震,用ROI_PAC软件处理4景日本ALOS卫星数据得到乌恰地震的三维同震形变场。距离向,断层北侧最大位移达39 cm,南侧最大位移达36 cm;方位向,北侧和南侧的最大位移分别为1.5 m和2 m。采用Okada弹性位错模型分析该地震的滑动分布,模拟得出的断层走滑角为48°,倾滑角为53°,深度为10.5 km,最大形变量2.5 m,分析得出该地震为向北逆冲兼左旋走滑。     

1986年门源地震(MS 6.4)过程地形变演化特征及块体模型解析

分析了1986年门源6．4级地震震源周边地区不同尺度跨断层形变在地震过程中的演化特征,依据非连续变形数值分析方法(DDA)模拟了震源破裂过程引起的地表垂直位移速率变化．分析表明,地表跨越断层的剖面形变变化是垂直和水平栽荷共同作用的结果,而垂直力源的影响可能是引起本次地震形变变化的主要因素．这与地质构造分析和震源机制结果是吻合的．     

PRIMARILY STUDY ON FEATURES OF SURFACE RUPTURES INDUCED BY THE 2016 MW7.8 KARKOURA EARTHQUAKE,NEW ZEALAND

The surface ruptures produced by the 2016 M_W7.8 Karkoura earthquake, New Zealand are distributed in a belt with~170km long and~35km wide, trending generally in the NE-SW direction. There are at least 12 faults on which meter-scale displacements are identified and they were formed across two distinct seismotectonic provinces with fundamental different characteristics(Hamling et al., 2017; Litchfield et al., 2017). Although the trending directions of the seismic surface ruptures vary greatly at different locations, the ruptured faults can be generally divided into two groups with the NE to NEE direction and the NNW to N direction, respectively. The faults in the NNW-near NS direction are nearly parallel with 40~50km apart and featured by reverse movement with the maximum displacement of 5~6m. The faults in the NE-NNE direction, with the maximum of 25~30km apart are not continuous and featured by the dextral strike slip with the largest displacement of 10~12m. Even if some faults along the NE-NEE direction are end to end connected, their strikes differ by about 30°. The combination styles of the strike-slip fault surface ruptures along the NE-NEE direction can be merged into 3 categories, including en-echelon, bifurcation and parallel patterns. The scales of the fault surface ruptures with the same structural style could be obviously different in different areas, which results in significant changes in the widths of deformation zone, from tens of meters to hundreds of meters. En-echelon distributed surface rupture(section)can appear as a combination belt of meter-scale to dozens of meter-scale shear fracture with bulge and compressional shear fractures, and also can be characterized by the combination of the left-step en-echelon tensile shear fractures with a length of more than one hundred meters. The step-overs between surface rupture sections are clearly different in sizes, which can be dozens of meters, hundreds of meters to several kilometers. The spacing between parallel surface ruptures can be several meters, dozens of meters to several kilometers. Besides, as one of the prominent characteristics, the seismic surface ruptures caused by the Karkoura earthquake broke through the known distribution pattern of active faults. The surface ruptures can occur either on the previously thought inactive or unmapped faults, or break through the distribution range of previously realized active faults in the striking or lateral direction. The basic features about the distribution and widths of the surface ruptures induced by the 2016 M_W7.8 Karkoura earthquake, New Zealand presented in this paper might be helpful for understanding some seismic problems such as complex corresponding relationship between the active faults and the deep seismogenic structure, and the necessary measurements for engineering crossing active faults.