1975年海城地震和1976年和林格尔强震与1976年唐山大震的关系讨论

从北纬40°纬向地震带活动的关系探讨了华北3次强震发生的同时性的原因,这条EW向地震带首先发生的地震传递的能量促使唐山积累了发震能量,而另一条由邢台、河间NE走向的下地壳蠕滑断层的能量传至唐山断层后,解锁了该孕震断层,从而发生了唐山大震。地震波的触发作用也是同步的另一个原因,除了震动的断层面的直接触发之外,对组合模式中单地震波的积累与调整也会产生触发作用而改变其性质。另外还从中国8级大地震发生的25年周期讨论了同步性。     

Toward a physical undenstanding of earthquake scaling relations

In seismological literature, there exist two competing theories (the so-calledW model andL model) treating earthquake scaling relations between mean slip and rupture dimension and between seismic moment and rupture dimension. The core of arguments differentiating the two theories is whether the mean slip should scale with the rupture width or with the rupture length for large earthquakes. In this paper, we apply the elastic theory of dislocation to clarify the controversy. Several static dislocation models are used to simulate strike-slip earthquakes. Our results show that the mean slip scales linearly with the rupture width for small earthquakes with a rupture length smaller than the thickness of the seismogenic layer. However, for large earthquakes with a rupture length larger than the thickness of the seismogenic layer, our models show a more complicated scaling relation between mean slip and rupture dimension. When the rupture length is smaller than a cross-over length, the mean slip scales nearly linearly with the rupture length. When the rupture length is larger than a cross-over length, the mean slip approaches asymptotically a constant value and scales approximately with the rupture width. The cross-over length is a function of the rupture width and is about 75 km for earthquakes with a saturated rupture width of 15 km. We compare our theoretical predictions with observed source parameters of some large strike-slip earthquakes, and they match up well. Our results also suggest that when large earthquakes have a fixed aspect ratio of rupture length to rupture width (which seems to be the case for most subduction earthquakes) the mean slip scales with the rupture dimension in the same way as small earthquakes.     

InSAR数据约束的2021年5月21日云南漾濞MS6.4地震发震构造研究

运用Sentinel-1A卫星数据和D-InSAR技术,获取2021-05-21云南漾濞M_S6.4地震的同震形变场。结果显示,漾濞地震同震形变场长轴近NW展布升降轨形变场符号相反,视线向最大沉降量和抬升量为0.1 m。InSAR同震形变场反演的滑动分布主要集中在沿走向2～12 km,倾向1～9 km的范围内,最大滑动量0.35 m,发震断层长9.8 km、宽4 km,滑动量主要集中在地下3～6 km范围内,滑动角-146.7°。同震位移场及滑动分布模型反映本次地震为发震断层的右旋走滑事件,地震破裂未达到地表。断层模型反演结果显示,矩震级为M_W6.1,发震断层以北西走向右旋走滑运动为主,初步认为本次M_W6.1地震发震断裂可能是一条NW向的维西—乔后断裂西侧的隐伏次生断裂。     

青藏高原东部及邻区地震断层面上的潮汐应力触发效应

根据已发布震源机制解目录（哈佛大学CMT）,将青藏高原东部及邻区划分为5个构造应力场分区,并对各分区的地震逐个计算其发震断层面上的固体潮汐正应力、剪应力、库仑破裂应力及相位角,分析潮汐应力分量对不同类型发震断层的作用效果及其随深度变化特征.基于库仑破裂应力判断准则,研究潮汐应力对各种类型地震的触发作用;基于Schuster检验方法,统计分析潮汐应力对各个震级档、不同构造类型地震的影响.综合运用上述两种分析方法,探讨潮汐应力对不同震级地震以及处于不同构造块体、发震断层、震源深度地震的触发机制.结果表明,潮汐应力对印度块体和拉萨块体的正断和逆断型地震,滇缅泰块体、印支块体和松潘-甘孜块体的走滑和斜滑型地震,川滇菱形块体的斜滑型地震均存在不同程度的触发效应,且触发效应的强弱依赖于震源深度、震级大小、发震断层类型及其所处区域构造应力场.     

Retardations in fault creep rates before local moderate earthquakes along the San Andreas fault system,central California

Records of shallow aseismic slip (fault creep) obtained along parts of the San Andreas and Calaveras faults in central California demonstrate that significant changes in creep rates often have been associated with local moderate earthquakes. An immediate postearthquake increase followed by gradual, long-term decay back to a previous background rate is generally the most obvious earthquake effect on fault creep. This phenomenon, identified as aseismic afterslip, usually is characterized by above-average creep rates for several months to a few years. In several cases, minor step-like movements, called coseismic slip events, have occurred at or near the times of mainshocks. One extreme case of coseismic slip, recorded at Cienega Winery on the San Andreas fault 17.5 km southeast of San Juan Bautista, consisted of 11 mm of sudden displacement coincident with earthquakes ofM _L =5.3 andM _L =5.2 that occurred 2.5 minutes apart on 9 April 1961. At least one of these shocks originated on the main fault beneath the winery. Creep activity subsequently stopped at the winery for 19 months, then gradually returned to a nearly steady rate slightly below the previous long-term average.The phenomena mentioned above can be explained in terms of simple models consisting of relatively weak material along shallow reaches of the fault responding to changes in load imposed by sudden slip within the underlying seismogenic zone. In addition to coseismic slip and afterslip phenomena, however, pre-earthquakeretardations in creep rates also have been observed. Onsets of significant, persistent decreases in creep rates have occurred at several sites 12 months or more before the times of moderate earthquakes. A 44-month retardation before the 1979M _L =5.9 Coyote Lake earthquake on the Calaveras fault was recorded at the Shore Road creepmeter site 10 km northwest of Hollister. Creep retardation on the San Andreas fault near San Juan Bautista has been evident in records from one creepmeter site for the past 5 years. Retardations with durations of 21 and 19 months also occurred at Shore Road before the 1974 and 1984 earthquakes ofM _L =5.2 andM _L =6.2, respectively.Although creep retardation remains poorly understood, several possible explanations have been discussed previously. (1) Certain onsets of apparent creep retardation may be explained as abrupt terminations of afterslip generated from previous moderate-mainshock sequences. (2) Retardations may be related to significant decreases in the rate of seismic and/or aseismic slip occurring within or beneath the underlying seismogenic zone. Such decreases may be caused by changes in local conditions related to growth of asperities, strain hardening, or dilatancy, or perhaps by passage of stress-waves or other fluctuations in driving stresses. (3) Finally, creep rates may be lowered (or increased) by stresses imposed on the fault by seismic or aseismic slip on neighboring faults. In addition to causing creep-rate increases or retardations, such fault interactions occasionally may trigger earthquakes.Regardless of the actual mechanisms involved and the current lack of understanding of creep retardation, it appears that shallow fault creep is sensitive to local and regional effects that promote or accompany intermediate-term preparation stages leading to moderate earthquakes. A strategy for more complete monitoring of fault creep, wherever it is known to occur, therefore should be assigned a higher priority in our continuing efforts to test various hypotheses concerning the mechanical relations between seismic and aseismic slip.     

2008年5月12日四川汶川8.0级地震与部分余震的震源机制解

采用区域和远台Pn或Pg初至波初动符号, 利用下半球等面积投影, 求解了2008年5月12日四川汶川8.0级地震和截止到2008年12月10日发生的部分4级以上余震的震源机制解。 汶川8.0级地震的震源机制为: 节面Ⅰ的走向为5°, 倾角为48°, 滑动角为39°; 节面Ⅱ的走向为247°, 倾角为62°, 滑动角为131°。 P轴方位角为309°, 仰角为8°, T轴方位角为208°, 仰角为54 °, B轴方位角为44°, 仰角为35°。 结合地质构造和余震空间分布, 可以确定节面Ⅱ为发震断层面。 根据震源机制解, 引发本次地震的断层活动主要表现为逆冲, 主破裂面为S67°W与该地震所在断层的走向基本一致(断裂总体走向N45°E)^[1]; 主压应力轴P轴为N51 °W, 主压应力轴P轴方位与该区域构造应力场方向基本一致。 根据余震震源机制解结果, 龙门山断裂带南段发生的余震与北段发生的余震的震源机制都具有优势分布, 且两者差异明显。 早期发生在南段的余震的破裂是以逆倾滑动为主, 兼有走向滑动; 而随着时间的推移, 余震向北段迁移, 在龙门山构造的北段地震震源的破裂方式以走向滑动为主, 兼有一定的逆倾滑动; 龙门构造带南段震源应力场受主震应力场的控制, 而龙门构造带北段震源应力场不仅受区域应力场的影响, 还受主震应力场的影响。     

Coseismic deformation of the 2021 MW7.4 Maduo earthquake from joint inversion of InSAR,GPS, and teleseismic data

The M_W7.4 Maduo earthquake occurred on 22 May 2021 at 02:04 CST with a large-expansion surface rupture. This earthquake was located in the Bayan Har block at the eastern Tibetan Plateau, where eight earthquakes of M_S >7.0 have occurred in the past 25 years. Here, we combined interferometric synthetic aperture radar, GPS, and teleseismic data to study the coseismic slip distribution, fault geometry, and dynamic source rupture process of the Maduo earthquake. We found that the overall coseismic deformation field of the Maduo earthquake is distributed in the NWW-SEE direction along 285°. There was slight bending at the western end and two branches at the eastern end. The maximum slip is located near the eastern bending area on the northern branch of the fault system. The rupture nucleated on the Jiangcuo fault and propagated approximately 160 km along-strike in both the NWW and SEE directions. The characteristic source rupture process of the Maduo earthquake is similar to that of the 2010 M_W6.8 Yushu earthquake, indicating that similar earthquakes with large-expansion surface ruptures and small shallow slip deficits can occur on both the internal fault and boundary fault of the Bayan Har block.     

A Probabilistic Assessment of Earthquake Hazard Parameters in NW Himalaya and the Adjoining Regions

The maximum likelihood estimation method is applied to study the geographical distribution of earthquake hazard parameters and seismicity in 28 seismogenic source zones of NW Himalaya and the adjoining regions. For this purpose, we have prepared a reliable, homogeneous and complete earthquake catalogue during the period 1500–2010. The technique used here allows the data to contain either historical or instrumental era or even a combination of the both. In this study, the earthquake hazard parameters, which include maximum regional magnitude (M _max), mean seismic activity rate (λ), the parameter b (or β?=?b/log e) of Gutenberg–Richter (G–R) frequency-magnitude relationship, the return periods of earthquakes with a certain threshold magnitude along with their probabilities of occurrences have been calculated using only instrumental earthquake data during the period 1900–2010. The uncertainties in magnitude have been also taken into consideration during the calculation of hazard parameters. The earthquake hazard in the whole NW Himalaya region has been calculated in 28 seismogenic source zones delineated on the basis of seismicity level, tectonics and focal mechanism. The annual probability of exceedance of earthquake (activity rate) of certain magnitude is also calculated for all seismogenic source zones. The obtained earthquake hazard parameters were geographically distributed in all 28 seismogenic source zones to analyze the spatial variation of localized seismicity parameters. It is observed that seismic hazard level is high in Quetta-Kirthar-Sulaiman region in Pakistan, Hindukush-Pamir Himalaya region and Uttarkashi-Chamoli region in Himalayan Frontal Thrust belt. The source zones that are expected to have maximum regional magnitude (M _max) of more than 8.0 are Quetta, southern Pamir, Caucasus and Kashmir-Himanchal Pradesh which have experienced such magnitude of earthquakes in the past. It is observed that seismic hazard level varies spatially from one zone to another which suggests that the examined regions have high crustal heterogeneity and seismotectonic complexity.     

Interaction between shallow and subcrustal dislocations on a normal fault

We propose a model which may explain seismic sequences which are often observed in seismogenic regions, as for example in the Apenninic chain (Italy). In particular, we consider a normal fault and earthquakes taking place at different depths: a first shock in a shallower layer and a second in a deeper one. The normal fault is embedded in a viscoelastic half-space. As a consequence of the rheology, there are two different brittle layers, a shallower and a deeper one, where earthquakes can nucleate. Between these two layers, the rheological behavior is ductile. The thicknesses of the layers depend on the geothermal profile that is calculated taking into account the profile of the thermal and rheological parameters with depth. The fault plane, crossing layers with different rheological behavior, is heterogeneous in respect to the slip style: seismic in the brittle layers, aseismic in the ductile layer. Dislocations in the shallower layer are assumed to produce aseismic slip in the area of the fault belonging to the ductile layer. The stress concentrated, by the seismic and aseismic dislocations, on the fault plane section in the deeper brittle layer is evaluated. It is compared with the tectonic stress rate in order to calculate how much earlier the second earthquake would occur compared to if just the bare tectonic sstress was acting. It results that such an advance is comparable with typical recurrence times of earthquakes and so a mechanism of interaction between different seismic sources, mediated by aseismic slip, can be supposed. The strains and displacements at the Earth’s surface due to seismic and aseismic slip are calculated. They are large enough to be detected by present geodetic techniques.     

Temporal evolution of the focal mechanism consistency of the 2021 Yangbi MS 6.4 earthquake sequence in Yunnan

Using the Cut And Paste (CAP) method, we invert the focal mechanism of 38 moderate earthquakes (M_S ≥ 3.0) recorded by Yunnan seismic network and analyze the corresponding focal mechanism consistency based on the minimum spatial rotation angle. Our results indicate that the M_S 6.4 mainshock is induced by a lateral strike slip fault (with a rake angle of ～ ?165°) and a little normal-faulting component event along a nearly vertical plane (dipping angle～ 79° and strike ～138°). Combining our results with high resolution catalog, we argue that the seismogenic fault of this earthquake sequence is a secondary fault western to the major Weixi-Qiaohou-Weishan fault. The focal mechanism evolution can be divided into three periods. During the first period, the foreshock sequence, the focal mechanism consistency is the highest (KA<36°); during the second period which is shortly after the mainshock, the focal mechanism shows strong variation with KA ranging from 8° to 110°; during the third period, the seismicity becomes weak and the focal mechanism of the earthquakes becomes more consistent than the second period (18°<KA<73°). We suggest that the KA, to some extent, represents the coherence between local tectonic stress regime and the stress state of each individual earthquake. Furthermore, high focal mechanism consistency and high linearity of seismic distribution may serve as indicators for the identification of foreshock sequence.     

The December 26, 2004, off the west coast of northern Sumatra,Indonesia, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript>=9.0, earthquake and the critical-point-like model of earthquake preparation

Long-term seismic activity prior to the December 26, 2004, off the west coast of northern Sumatra, Indonesia, M _W=9.0 earthquake was investigated using the Harvard CMT catalogue. It is observed that before this great earthquake, there exists an accelerating moment release (AMR) process with the temporal scale of a quarter century and the spatial scale of 1 500 km. Within this spatial range, the M _W=9.0 event falls into the piece-wise power-law-like frequency-magnitude distribution. Therefore, in the perspective of the critical-point-like model of earthquake preparation, the failure to forecast/predict the approaching and/or the size of this earthquake is not due to the physically intrinsic unpredictability of earthquakes.     

Evolution of Stress Fields and Faulting in Seismic Zones

—Measurements indicate that stress magnitudes in the crust are normally limited by the frictional equilibrium on pre-existing, optimally oriented faults. Fault zones where these limitations are frequently reached are referred to as seismic zones. Fault zones in the crust concentrate stresses because their material properties are different from those of the host rock. Most fault zones are spatially relatively stable structures, however the associated seismicity in these zones is quite variable in space and time. Here we propose that this variability is attributable to stress-concentration zones that migrate and expand through the fault zone. We suggest that following a large earthquake and the associated stress relaxation, shear stresses of a magnitude sufficient to produce earthquakes occur only in those small parts of the seismic zone that, because of material properties and boundary conditions, encourage concentration of shear stress. During the earthquake cycle, the conditions for seismogenic fault slip migrate from these stress-concentration regions throughout the entire seismic zone. Thus, while the stress-concentration regions continue to produce small slips and small earthquakes throughout the seismic cycle, the conditions for slip and earthquakes are gradually reached in larger parts of, and eventually the whole, seismogenic layer of the seismic zone. Prior to the propagation of an earthquake fracture that gives rise to a large earthquake, the stress conditions in the zone along the whole potential rupture plane must be essentially similar. This follows because if they were not, then, on entering crustal parts where the state of stress was unfavourable to this type of faulting, the fault propagation would be arrested. The proposed necessary homogenisation of the stress field in a seismic zone as a precursor to large earthquakes implies that by monitoring the state of stress in a seismic zone, its large earthquakes may possibly be forecasted. We test the model on data from Iceland and demonstrate that it broadly explains the historical, as well as the current, patterns of seismogenic faulting in the South Iceland Seismic Zone.     

2013年8月河北蔚县小震群遗漏地震检测与发震构造分析

华北地区近年来小震群活动频繁, 在有数字波形记录的中强地震相对缺乏的背景下, 小震群发震构造精细研究可为华北地区地震危险性分析和地震趋势判断提供重要依据. 本文利用匹配滤波技术对2013年8月22—25日河北蔚县小震群遗漏地震事件进行检测, 并通过地震精定位和震源机制求解分析此次震群的发震构造. 计算结果显示, 通过互相关扫描检测到18次被地震台网常规分析遗漏的地震, 约为地震目录给出的13次地震事件的1.38倍. 该震群发震构造有北东向和北西向两组断裂, 震群活动前期以北东向构造活动为主, 后期地震主要发生在北西向构造, 北西向构造在此次震群活动中地震频度和强度均高于北东向构造. 震源机制计算结果显示北西向构造发震机制以正断拉张为主.      

SIMULATION OF SEISMIC RISK IN THE DALIANGSHAN SUB-BLOCK AND ADJACENT AREAS USING THE NONLINEAR FRICTION FEM METHOD

Most earthquakes result from fault activity under heterogeneous loading and complex physical properties, also affected by fault structure and interaction between faults. Such a complicated mechanism makes often failures of the "seismic gap" theory in the effort of medium-and long-term earthquake prediction. This study attempts to address this issue using the finite element method(FEM).The friction behavior of faults can be used to simulate the non-uniformity of rupture processes of the seismogenic structure. So we use the FEM containing non-linear friction to simulate fault ruptures in the Daliangshan sub-block and adjacent areas, and compare the results with time-space evolution of historical M_S ≥ 7 earthquakes since 1840 in this region. In the simulation, the sequence of large-batch fault contact nodes change from "stick state" to "slip state" in short time, which mimics the sudden fault slip and the occurrence of major earthquakes. The results show that the fault breaking lengths from simulation are largely consistent with the magnitudes of historical earthquakes in the study area, such as the 1850 Puge-Xichang M_S7.5, and 1887 Shiping M_S7.0 earthquakes. The simulation also shows the development of seismic gaps and "gap breaks" by major earthquakes on the Xianshuihe fault, such as 1955 Kangding M_S7.5 earthquake. Especially, the results illustrated the very long time of the seismogenic process of the 2008 Wenchuan M_S8.0 earthquake, and the corresponding sudden big rupture along the Longmenshan Fault, which is very similar to the observed surface rupture and very long incubation time and sudden co-seismic process. Then, this simulation is further applied to long-term earthquake prediction for the study area by calculation on a much longer time. The simulation results suggest that the Xiaojiang fault and the Zemuhe fault have relatively higher seismic risk, while moderate-sized earthquakes might occur on the Daliangshan fault and the Aninghe fault, and major earthquakes might rupture the northern segment of the Xianshuihe fault in a much longer time.     

Studies on the June 23, 2010 north Ottawa M W 5.2 earthquake and vicinity seismicity

The basic parameters for the earthquake with a moment magnitude (M _W) of 5.2 on the 23rd of June 2010 have been investigated. The earthquake occurred on a hidden fault in the northwest direction about 60?km north-northeast of Ottawa in the Western Quebec Seismic Zone (WQSZ) and had a focal depth of about 21?km. The focal mechanism was a thrust type with strike in the northwest direction and dipping in the northeast direction. The relative relocations of seven larger aftershocks show that the source rupture area was about 6?km². The b value of the aftershock sequence was 0.8?C1.0, and the decay rate of the aftershocks was faster than normal cases. The dominant seismogenic depths are about 12 to 22?km in most parts of the WQSZ, while the seismogenic depth along the Ottawa?CBonnechere Graben can be as deep as 28?km. Based on the seismic activity in the WQSZ and vicinity since 1961, it seems that the periods of moderate earthquakes are about 6?C10?years.     

Rupture model of the 2013 M_W 6.6 Lushan (China) earthquake constrained by a new GPS data set and its effects on potential seismic hazard

Vertical records are critically important when determining the rupture model of an earthquake, especially a thrust earthquake. Due to the relatively low fitness level of near-field vertical displacements, the precision of previous rupture models is relatively low, and the seismic hazard evaluated thereafter should be further updated. In this study, we applied three-component displacement records from GPS stations in and around the source region of the 2013 MW6.6 Lushan earthquake to re-investigate the rupture model.To improve the resolution of the rupture model, records from both continuous and campaign GPS stations were gathered, and secular deformations of the GPS movements were removed from the records of the campaign stations to ensure their reliability. The rupture model was derived by the steepest descent method(SDM), which is based on a layered velocity structure. The peak slip value was about 0.75 m, with a seismic moment release of 9.89 × 10~(18) N·m, which was equivalent to an M_W6.6 event. The inferred fault geometry coincided well with the aftershock distribution of the Lushan earthquake. Unlike previous rupture models, a secondary slip asperity existed at a shallow depth and even touched the ground surface. Based on the distribution of the co-seismic ruptures of the Lushan and Wenchuan earthquakes, post-seismic relaxation of the Wenchuan earthquake, and tectonic loading process, we proposed that the seismic hazard is quite high and still needs special attention in the seismic gap between the two earthquakes.     

鲜水河断裂带上特征地震的初步研究

特征地震是大地震原地重复的重要表现形式.现有资料的初步研究表明,鲜水河断裂带上大地震属特征地震模式,其地震破裂长度、同震位错量以及断层错动方式,在原地保持较长时间的一致性.由于大地震屡屡在原地重复发生,沿断裂特定地段累积位错分布与一次地震的位错相一致,从而导致断层滑动速率的同步变化.本文以1973年炉霍地震为例,研究了鲜水河断裂的特征地震现象.该段的地震活动属特征地震模式,不服从古登堡-里克特的线性震级频度关系.特征地震不仅对断错地貌、滑动速率有重要的影响,由于这种地震模式是以特定震级的大地震为主导,几乎没有中等震级地震发生,这对地震活动性研究也具有重要的意义.      

Stress drop assessment of the August 8, 2017, Jiuzhaigou earthquake sequence and its tectonic implications*

By using a broadband Lg attenuation model developed for the Tibetan Plateau, we isolate source terms by removing attenuation and site effects from the observed Lg-wave displacement spectra of the M7.0 earthquake that occurred on August 8, 2017, in Jiuzhaigou, China, and its aftershock sequence. Thus, the source parameters, including the scalar seismic moment, corner frequency and stress drop, of these events can be further estimated. The estimated stress drops vary from 47.1 kPa to 7149.6 kPa, with a median value of 59.4 kPa and most values falling between 50 kPa and 75 kPa. The estimated stress drops show significant spatial variations. Lower stress drops were mainly found close to the mainshock and on the seismogenic fault plane with large coseismic slip. In contrast, the highest stress drop was 7.1 MPa for the mainshock, and relatively large stress drops were also found for aftershocks away from the major seismogenic fault and at depths deeper than the zone with large coseismic slip. By using a statistical method, we found self-similarity among some of the aftershocks with a nearly constant stress drop. In contrast, the stress drop increased with the seismic moment for other aftershocks. The amount of stress released during earthquakes is a fundamental characteristic of the earthquake rupture process. As such, the stress drop represents a key parameter for improving our understanding of earthquake source physics.     

The seismicity and tectonic stress field characteristics of the Longmenshan fault zone before the Wenchuan MS8.0 earthquake

The seismicity of Longmenshan fault zone and its vicinities before the 12 May 2008 Wenchuan M_S8.0 earthquake is studied. Based on the digital seismic waveform data observed from regional seismic networks and mobile stations, the focal mechanism solutions are determined. Our analysis results show that the seismicities of Longmenshan fault zone before the 12 May 2008 Wenchuan earthquake were in stable state. No obvious phenomena of seismic activity intensifying appeared. According to focal mechanism solutions of some small earthquakes before the 12 May 2008 Wenchuan earthquake, the direction of principal compressive stress P-axis is WNW-ESE. The two hypocenter fault planes are NE-striking and NW-striking. The plane of NE direction is among N50°?70°E, the dip angles of fault planes are 60°?70° and it is very steep. The faultings of most earthquakes are dominantly characterized by dip-slip reverse and small part of faultings present strike-slip. The azimuths of principal compressive stress, the strikes of source fault planes and the dislocation types calculated from some small earthquakes before the 12 May 2008 Wenchuan earthquake are in accordance with that of the main shock. The average stress field of micro-rupture along the Longmenshan fault zone before the great earthquake is also consistent with that calculated from main shock. Zipingpu dam is located in the east side 20 km from the initial rupture area of the 12 May 2008 Wenchuan earthquake. The activity increment of small earthquakes in the Zipingpu dam is in the period of water discharging. The source parameter results of the small earthquakes which occurred near the initial rupture area of the 12 May 2008 Wenchuan earthquake indicate that the focal depths are 5 to 14 km and the source parameters are identical with that of earthquake.