SAR差分干涉测量技术在地震区域形变测量中的应用研究

合成孔径雷达(SAR)差分干涉测量是一种最新的大地形变测量遥感技术方法.利用欧洲空间局1996年4月15日获取的地球资源卫星ERS-1数据和4月16日获取的ERS-2数据,以及1997年12月2日获取的ERS-2数据,应用差分干涉测量技术对发生于1997年11月8日的西藏玛尼地震进行了提取区域形变场的应用研究.从得到的变化检测条纹图中可识别出地表破裂带,并可定量推算震中周围和两条断裂带附近的变形情况.分析结果表明:断裂为左行平移性质,断裂带南侧水平错动的最大变形量至少超过2.8 m,北侧水平错动的最大变形量至少超过1.75 m.差分干涉测量分析结果与地面调查及已有的相关研究资料符合得很好.     

活动断层分段研究

地震发生在活动断层上,而断层活动不一定都发震,为进行地震危险性分析,确定未来可能发生地震的地段和强度。因此,把活动断层分成地震破裂段进行研究,分析断层活动在时间上、空间上的不均匀性,了解破裂段特征和分段标志,确定破坏性地震发震的构造条件是极为重要的。     

利用相对反投影方法对2010年Mw 8.8级智利地震破裂过程成像

基于美国流动台网的三分量远震P波资料,利用相对反投影方法,对2010年2月27日Mw 8.8智利地震的破裂过程成像。根据美国地质调查局(USGS)给出的震源参数及源区地质资料,设定一个小倾角断层面。将远震P波反投影至该断层面上,获取震源破裂过程。结果表明,智利地震破裂由3个子事件组成,破裂长度至少为513 km,宽度至少为100 km。破裂从震中开始向北方和南方两个方向双侧扩展,但北侧的破裂明显强于南侧。北侧的破裂长度约为340 km,破裂持续时间约125 s,破裂速度约为2.87 km/s;而南侧破裂长度约为173 km,破裂持续时间约99 s,破裂速度约为1.84 km/s。另外,北侧两个子事件的能量释放高峰分别出现在破裂开始之后16 s和79 s,与初始破裂点的距离分别为75 km和230 km。南侧子事件的能量释放高峰出现时刻为80 s,与初始破裂点的距离为145 km。智利地震在破裂过程中可能受到了破裂区域的障碍体影响。     

Extent, duration and speed of the 2004 Sumatra-Andaman earthquake imaged by the Hi-Net array

The disastrous Sumatra-Andaman earthquake of 26 December 2004 was one of the largest ever recorded. The damage potential of such earthquakes depends on the extent and magnitude of fault slip. The first reliable moment magnitude estimate of 9.0 was obtained several hours after the Sumatra-Andaman earthquake, but more recent, longer-period, normal-mode analyses have indicated that it had a moment magnitude of 9.3, about 2.5 times larger. Here we introduce a method for directly imaging earthquake rupture that uses the first-arriving compressional wave and is potentially able to produce detailed images within 30 min of rupture initiation. We used the Hi-Net seismic array in Japan as an antenna to map the progression of slip by monitoring the direction of high-frequency radiation. We find that the rupture spread over the entire 1,300-km-long aftershock zone by propagating northward at roughly 2.8 km s(-1) for approximately 8 minutes. Comparisons with the aftershock areas of other great earthquakes indicate that the Sumatra-Andaman earthquake did indeed have a moment magnitude of approximately 9.3. Its rupture, in both duration and extent, is the longest ever recorded.     

联合slipBERI与D-InSAR技术的青海省玛多县地震形变场提取与模拟

针对2021年5月22日青海玛多县Ms7.4地震震区形变信息、形变特征以及滑动断裂特性的提取与模拟等问题,获取了玛多县地震区的Sentinel-1A影像,采用双轨差分干涉法,并优化各项参数,提取出同震形变场,利用slipBERI(slip from Bayesian Regularized Inversion)方法对断层的几何参数及形变场进行反演和模拟。结果表明：玛多地震同震形变场的形状近似于一个椭圆,断层整体呈西北-东南走向,其上部为沉降区,下部为隆升区,最大LOS(line of sight,视线向)形变分别为0.65m和0.81m。地震形变场的运动主要以东西方向的水平运动为主,并伴有明显的左旋走滑,断层上下方相对视线向运动可达1.50m,表明此次地震的地表破裂有明显的错位移动。通过分析形变信息和地表破裂特征,可以判断该破裂带位于巴颜喀拉块体,为昆仑山口-江口断裂,反演结果与观测结果相符,这表明观测结果较可靠。     

汶川地震地表破裂带宽度与断层上盘效应

以汶川地震产生的地表破裂宽度为研究对象, 以地表微变形、裂隙和鼓包、叠瓦状断层及反冲断层等为标志界定破裂的宽度, 查明中央破裂带和前山破裂带宽度的分布。结果表明, 中央断裂带与汉旺?白鹿断裂带的地表破裂宽度不同, 中央断裂带不同段落的地表破裂宽度也不相同。断层还表现出明显的断层上盘效应, 北川?映秀断裂上盘最小避让距离为70 m, 下盘最小避让距离为23.5 m, 断层两侧避让总宽度不应小于108 m; 汉旺?白鹿断裂上盘最小避让距离为35 m, 下盘最小避让距离20 m, 总避让宽度不应小于58 m。地表破裂宽度及断层上盘效应与破裂带陡坎高度、地震烈度分布、强震记录具有较好的一致性。     

巨形大震与全球变化

所指的巨形大震是指发震断层极长、断层面积甚大和矩震级Mw为9级和9级以上的地震。这种大震的发生不但与北半球和中国的气候降温有一定相关性,而且与中国的特大洪灾也有一定关系。2004年印尼巨形Mw9.0级大震已给我们提出了警告。     

The M S7.1 Yushu earthquake surface rupture and large historical earthquakes on the Garzê-Yushu Fault

As revealed by field investigations, the co-seismic surface rupture zone of the 2010 M_S7.1 Yushu earthquake, Qinghai is a char-acteristic sinistral strike-slip feature consisting of three distinct sinistral primary ruptures, with an overall strike of 310°–320° and a total length of 31 km. In addition, an approximately 2-km-long en-echelon tensile fissure zone was found east of Longbao Town; if this site is taken as the north end of the rupture zone, then the rupture had a total length of ~51 km. The surface rupture zone is composed of a series of fissures arranged in an en-echelon or alternating relationship between compressive bulges and tensile fissures, with a measured maximum horizontal displacement of 1.8 m. The surface rupture zone extends along the mapped Garzê-Yushu Fault, which implicates it as the seismogenic fault for this earthquake. Historically, a few earthquakes with a magnitude of about 7 have occurred along the fault, and additionally traces of paleoearthquakes are evident that characterize the short-period recurrence interval of large earthquakes here. Similar to the seismogenic process of the 2008 Wenchuan earthquake, the Yushu earthquake is also due to the stress accumulation and release on the block boundaries resulting from the eastward expansion of Qinghai-Tibet Plateau. However, in contrast with the Wenchuan earthquake, the Yushu earthquake had a sinistral strike-slip mechanism resulting from the uneven eastward extrusion of the Baryan Har and Sichuan-Yunnan fault blocks.     

Particle size and energetics of gouge from earthquake rupture zones

Grain size reduction and gouge formation are found to be ubiquitous in brittle faults at all scales, and most slip along mature faults is observed to have been localized within gouge zones. This fine-grain gouge is thought to control earthquake instability, and thus understanding its properties is central to an understanding of the earthquake process. Here we show that gouge from the San Andreas fault, California, with approximately 160 km slip, and the rupture zone of a recent earthquake in a South African mine with only approximately 0.4 m slip, display similar characteristics, in that ultrafine grains approach the nanometre scale, gouge surface areas approach 80 m2 g(-1), and grain size distribution is non-fractal. These observations challenge the common perception that gouge texture is fractal and that gouge surface energy is a negligible contributor to the earthquake energy budget. We propose that the observed fine-grain gouge is not related to quasi-static cumulative slip, but is instead formed by dynamic rock pulverization during the propagation of a single earthquake.     

Seismology: earthquake risk on the Sunda trench

On 28 March 2005 the Sunda megathrust in Indonesia ruptured again, producing another great earthquake three months after the previous one. The rupture was contiguous with that of the December 2004 Sumatra-Andaman earthquake, and is likely to have been sparked by local stress, although the triggering stresses at its hypocentre were very small - of the order of just 0.1 bar. Calculations show that stresses imposed by the second rupture have brought closer to failure the megathrust immediately to the south, under the Batu and Mentawai islands, and have expanded the area of increased stress on the Sumatra fault. Palaeoseismologic studies show that the Mentawai segment of the Sunda megathrust is well advanced in its seismic cycle and is therefore a good candidate for triggered failure.     

Predicting the endpoints of earthquake ruptures

The active fault traces on which earthquakes occur are generally not continuous, and are commonly composed of segments that are separated by discontinuities that appear as steps in map-view. Stress concentrations resulting from slip at such discontinuities may slow or stop rupture propagation and hence play a controlling role in limiting the length of earthquake rupture. Here I examine the mapped surface rupture traces of 22 historical strike-slip earthquakes with rupture lengths ranging between 10 and 420 km. I show that about two-thirds of the endpoints of strike-slip earthquake ruptures are associated with fault steps or the termini of active fault traces, and that there exists a limiting dimension of fault step (3-4 km) above which earthquake ruptures do not propagate and below which rupture propagation ceases only about 40 per cent of the time. The results are of practical importance to seismic hazard analysis where effort is spent attempting to place limits on the probable length of future earthquakes on mapped active faults. Physical insight to the dynamics of the earthquake rupture process is further gained with the observation that the limiting dimension appears to be largely independent of the earthquake rupture length. It follows that the magnitude of stress changes and the volume affected by those stress changes at the driving edge of laterally propagating ruptures are largely similar and invariable during the rupture process regardless of the distance an event has propagated or will propagate.     

Static slip model of the M w 9.0 Tohoku (Japan) earthquake: Results from joint inversion of terrestrial GPS data and seafloor GPS/acoustic data

Based on co-seismic displacements recorded by terrestrial GPS stations and seafloor GPS/acoustic stations, the static slip model of the 2011 Mw 9.0 Tohoku earthquake was determined by inverting the data using a layered earth model. According to a priori information, the rupture surface was modeled with a geometry that is close to the actual rupture, in which the fault dip angle increases with depth and the fault strike varies with the trend of the trench. As shown by the results inferred from the joint inversion, the "geodetic" moment is 3.68 × 10 22 Nm, corresponding to Mw 9.01, and the maximum slip is positioned at a depth of 13.5 km with a slip magnitude of 45.8 m. Rupture asperities with slip exceeding 10 m are mainly distributed from 39.6 to 36.97°N, over a length of almost 240 km along the trench. The slip was mostly concentrated at depths shallower than 40 km, up-dip of the hypocenter. "Checkerboard" tests reveal that a joint inversion of multiple datasets can resolve the slip distribution better than an inversion with terrestrial GPS data only-especially when aiming to resolve slip at shallow depths. Thus, the joint inversion results obtained by this work may provide a more reliable slip model than the results of other studies that are only derived from terrestrial GPS data or seismic waveform data.     

Damage to the shallow Landers fault from the nearby Hector Mine earthquake

Crustal faults have long been identified as sites where localized sliding motion occurs during earthquakes, which allows for the relative motion between adjacent crustal blocks. Although there is a growing awareness that we must understand the evolution of fault systems on many timescales to relate present-day crustal stresses and fault motions to geological structures formed in the past, fault-zone damage and healing have been documented quantitatively in only a few cases. We have been monitoring the healing of damage on the shallow Johnson Valley fault after its rupture in the 1992 magnitude-7.3 Landers earthquake, and here we report that this healing was interrupted in 1999 by the magnitude-7.1 Hector Mine earthquake rupture, which occurred 20-30 km away. The Hector Mine earthquake both strongly shook and permanently strained the Johnson Valley fault, adding damage discernible as a temporary reversal of the healing process. The fault has since resumed the trend of strength recovery that it showed after the Landers earthquake. These observations lead us to speculate that fault damage caused by strong seismic waves may help to explain earthquake clustering and seismicity triggering by shaking, and may be involved in friction reduction during faulting.     

Fracture surface energy of the Punchbowl fault, San Andreas system

Fracture energy is a form of latent heat required to create an earthquake rupture surface and is related to parameters governing rupture propagation and processes of slip weakening. Fracture energy has been estimated from seismological and experimental rock deformation data, yet its magnitude, mechanisms of rupture surface formation and processes leading to slip weakening are not well defined. Here we quantify structural observations of the Punchbowl fault, a large-displacement exhumed fault in the San Andreas fault system, and show that the energy required to create the fracture surface area in the fault is about 300 times greater than seismological estimates would predict for a single large earthquake. If fracture energy is attributed entirely to the production of fracture surfaces, then all of the fracture surface area in the Punchbowl fault could have been produced by earthquake displacements totalling <1 km. But this would only account for a small fraction of the total energy budget, and therefore additional processes probably contributed to slip weakening during earthquake rupture.     

Plate-boundary deformation associated with the great Sumatra-Andaman earthquake

The Sumatra-Andaman earthquake of 26 December 2004 is the first giant earthquake (moment magnitude M(w) > 9.0) to have occurred since the advent of modern space-based geodesy and broadband seismology. It therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events. Here we report estimates of the ground displacement associated with this event, using near-field Global Positioning System (GPS) surveys in northwestern Sumatra combined with in situ and remote observations of the vertical motion of coral reefs. These data show that the earthquake was generated by rupture of the Sunda subduction megathrust over a distance of >1,500 kilometres and a width of <150 kilometres. Megathrust slip exceeded 20 metres offshore northern Sumatra, mostly at depths shallower than 30 kilometres. Comparison of the geodetically and seismically inferred slip distribution indicates that approximately 30 per cent additional fault slip accrued in the 1.5 months following the 500-second-long seismic rupture. Both seismic and aseismic slip before our re-occupation of GPS sites occurred on the shallow portion of the megathrust, where the large Aceh tsunami originated. Slip tapers off abruptly along strike beneath Simeulue Island at the southeastern edge of the rupture, where the earthquake nucleated and where an M(w) = 7.2 earthquake occurred in late 2002. This edge also abuts the northern limit of slip in the 28 March 2005 M(w) = 8.7 Nias-Simeulue earthquake.     

The deterministic nature of earthquake rupture

Understanding the earthquake rupture process is central to our understanding of fault systems and earthquake hazards. Multiple hypotheses concerning the nature of fault rupture have been proposed but no unifying theory has emerged. The conceptual hypothesis most commonly cited is the cascade model for fault rupture. In the cascade model, slip initiates on a small fault patch and continues to rupture further across a fault plane as long as the conditions are favourable. Two fundamental implications of this domino-like theory are that small earthquakes begin in the same manner as large earthquakes and that the rupture process is not deterministic--that is, the size of the earthquake cannot be determined until the cessation of rupture. Here we show that the frequency content of radiated seismic energy within the first few seconds of rupture scales with the final magnitude of the event. We infer that the magnitude of an earthquake can therefore be estimated before the rupture is complete. This finding implies that the rupture process is to some degree deterministic and has implications for the physics of the rupture process.     

Preliminary results pertaining to coseismic displacement and preseismic strain accumulation of the Lushan M S7.0 earthquake, as reflected by GPS surveying

This paper presents the coseismic displacement and preseismic deformation fields of the Lushan M _S7.0 earthquake that occurred on April 20, 2013. The results are based on GPS observations along the Longmenshan fault and within its vicinity. The coseismic displacement and preseismic GPS results indicate that in the strain release of this earthquake, the thrust rupture is dominant and the laevorotation movement is secondary. Furthermore, we infer that any possible the rupture does not reach the earth’s surface, and the seismogenic fault is most likely one fault to the east of the Guanxian-Anxian fault. Some detailed results are obtainable. (1) The southern segment of the Longmenshan fault is locked preceding the Lushan earthquake. After the Wenchuan earthquake, the strain accumulation rate in the southeast direction accelerates in the epicenter of the Lushan earthquake, and the angle between the principal compressional strain and the seismogenic fault indicates that a sinistral deformation background in the direction of the seismogenic fault precedes the Lushan earthquake. Therefore, it is evident that the Wenchuan M _S8.0 earthquake accelerated the pregnancy of the Lushan earthquake. (2) The coseismic displacements reflected by GPS data are mainly located in a region that is 230 km (NW direction) × 100 km (SW direction), and coseismic displacements larger than 10 mm lie predominantly in a 100-km region (NW direction). (3) On a large scale, the coseismic displacement shows thrust characteristics, but the associated values are remarkably small in the near field (within 70 km) of the earthquake fault. Meanwhile, the thrust movement in this 70-km region does not correspond with the attenuation characteristics of the strain release, indicating that the rupture of this earthquake does not reach the earth’s surface. (4) The laevorotation movements are remarkable in the 50-km region, which is located in the hanging wall that is close to the earthquake fault, and the corresponding values in this case correlate with the attenuation characteristics of the strain release.     

Source rupture process of Lushan M S7.0 earthquake, Sichuan, China and its tectonic implications

The source rupture process of the M _S7.0 Lushan earthquake was here evaluated using 40 long-period P waveforms with even azimuth coverage of stations. Results reveal that the rupture process of the Lushan M _S7.0 event to be simpler than that of the Wenchuan earthquake and also showed significant differences between the two rupture processes. The whole rupture process lasted 36 s and most of the moment was released within the first 13 s. The total released moment is 1.9×10¹⁹N m with M _W=6.8. Rupture propagated upwards and bilaterally to both sides from the initial point, resulting in a large slip region of 40 km×30 km, with the maximum slip of 1.8 m, located above the initial point. No surface displacement was estimated around the epicenter, but displacement was observed about 20 km NE and SW directions of the epicenter. Both showed slips of less than 40 cm. The rupture suddenly stopped at 20 km NE of the initial point. This was consistent with the aftershock activity. This phenomenon indicates the existence of significant variation of the medium or tectonic structure, which may prevent the propagation of the rupture and aftershock activity. The earthquake risk of the left segment of Qianshan fault is worthy of attention.     

九寨沟Ms 7.0级地震的左旋走滑作用与动力机制

2017年8月8日四川九寨沟发生的Ms 7.0级地震是继2008年汶川Ms 8.0地震、2013年芦山Ms 7.0级地震后在青藏高原东缘发生的又一次强震。本文通过综合分析九寨沟Ms 7.0级地震及历史地震的震源机制解、余震和历史地震分布、区域应力场、活动断层等资料,来揭示九寨沟地震的发震构造与动力机制。初步研究结果表明:(1)此次地震的震中位于塔藏断裂、岷江断裂和虎牙断裂之间的交汇区,显示活动断裂的交汇区对此次地震的发生具有控制作用;(2)发震断裂为虎牙断裂,断裂走向为北西西向,倾向南西,倾角较陡,属于高倾角左旋走滑型地震;(3)震中位于虎牙断裂北段的北部地震空区,充填了1973年和1976年4次大于Mw6.0级地震空区;(4)此次地震位于2008年汶川Ms 8.0级地震的库仑应力增加区,应是汶川地震的应力传递和触发的结果;(5)此次地震位于巴颜喀拉块体的东北部顶角区,青藏高原东缘下地壳流向北东方向的挤出是驱动此次地震的动力机制。     

Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system

The San Andreas fault in California is a mature continental transform fault that accommodates a significant fraction of motion between the North American and Pacific plates. The two most recent great earthquakes on this fault ruptured its northern and central sections in 1906 and 1857, respectively. The southern section of the fault, however, has not produced a great earthquake in historic times (for at least 250 years). Assuming the average slip rate of a few centimetres per year, typical of the rest of the San Andreas fault, the minimum amount of slip deficit accrued on the southern section is of the order of 7-10 metres, comparable to the maximum co-seismic offset ever documented on the fault. Here I present high-resolution measurements of interseismic deformation across the southern San Andreas fault system using a well-populated catalogue of space-borne synthetic aperture radar data. The data reveal a nearly equal partitioning of deformation between the southern San Andreas and San Jacinto faults, with a pronounced asymmetry in strain accumulation with respect to the geologically mapped fault traces. The observed strain rates confirm that the southern section of the San Andreas fault may be approaching the end of the interseismic phase of the earthquake cycle.