Advances in research on earthquake fluids hydrogeology in China:a review

Monitoring of subsurface fluid （underground fluid） is an important part of efforts for earthquake prediction in China. The nationwide network, which monitors groundwater level, water temperature, and radon and mercury in groundwater, has been constructed in the last decades. Large amounts of abnormal fluid changes before and after major earthquakes have been recorded, providing precious data for research in earthquake sciences. Many studies have been done in earthquake fluid hydrogeology in order to probe the nature of the earthquake. Much progress in earthquake fluid hydrogeology has been made in the last decades. The paper provides a review of the advances in research on earthquake fluid hydrogeology over the last 40 years in China. It deals with the following five aspects： （1） an introduction to the development history of monitoring networks construction; （2） cases of different subsurface fluid changes recorded before some major earthquakes which occurred in the last decades; （3） characteristics of subsurface fluid changes following major earthquakes; （4） mechanism of subsurface fluid changes before and following earthquakes; （5） application of earthquake fluids in the hydrogeology field.     

山东及邻区地震重新定位及其与活动构造的关系

Earthquakes taking place from 1975 to 2010 in and around Shandong Province are relocated using double-difference （HypoDD） and Hypoinvers 2000 （Hypo2000） methods, after correction of the onset times of seismic phases. The results show that the relocated seismicity is clearly associated with regional tectonics in space, and is also in agreement with the existence of deep faults imaged by wide-angle and deep seismic reflection profiling ; most of the focal depths are in the range of 5 - 25km, and there are clearly two predominant depths： 10km and 16km, which are inferred to be on the bottom of the upper crust and in the middle crust, respectively. The pattern of seismic activity indicates that moderate and strong earthquakes are likely to occur in the brittle-ductile transition zone between the upper and the lower crust, as the outcome of the deep tectonic dynamic process and the movement and deformation of faults in the upper and shallow crust under the regional stress field.     

玉树MS 7. 1 地震同震效应模拟与绝对重力检验

Using the results of aftershock relocation, inversion on seismic waves and InSAR results, and surface rupture displacements obtained by geological survey after the earthquake, this paper constructs a fault model of the Yushu Ms7. 1 earthquake. Based on rectangular dislocation theory in an elastic-viscoelastic layered medium, we have simulated the co- seismic deformation and gravity change with gravitational effect considered. The pictures show that the absolute gravity measuring point is beside the extremum of coseismic gravity change, and the numerical value reaches 25.02 x 10-Sm. s-2. After a discussion about the gravity changes before the earthquake and the coherence consistency between two FG-5 absolute gravimeters, we think that the measured value 27.2 × 10^-8 m· s^-2 at Yushu station is coseismic gravity change. It＇s coincident with the simulation results based on dislocation theory. Therefore it is a good tool to test the near-field changes found by dislocation theory.     

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

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

叠溪地震和地震地质概说

1调查经过1933年8月25日大约在下午2时半,四川茂县之北叠溪,忽然发生大地震,即时间附近群山崩倒,叠溪城全部毁灭,岷江断流。叠溪周围30里内,属于强烈地震区,在这一范围内,任何房屋庐舍,全部破坏,群山崖壁,全部崩坍。在叠溪圆径百里之内,所有道路都被破坏,不能通行。邮电线路,全部不通。当时只知茂汶     

Field survey around strong motion stations and its implications on the seismic intensity in the Lushan earthquake on April 20, 2013

The Ms7.0 Lushan earthquake on April 20, 2013 is another destructive event in China since the Ms8.0 Wenchuan earthquake in 2008 and Ms7.1 Yushu earth- quake in 2010. A large number of strong motion recordings were accumulated by the National Strong Motion Obser- vation Network System of China. The maximum peak ground acceleration （PGA） at Station 51BXD in Baoxing Country is recorded as -1,005.3 cm/s2, which is even larger than the maximum one in the Wenchuan earthquake. A field survey around three typical strong motion stations confirms that the earthquake damage is consistent with the issued map of macroseismic intensity. For the oscillation period 0.3-1.0 s which is the common natural period range of the Chinese civil building, a comparison shows that the observed response spectrums are considerably smaller than the designed values in the Chinese code and this could be one of the reasons that the macroseismic intensity is lower than what we expected despite the high amplitude of PGAs. The Housner spectral intensities from 16 stations are also basically correlated with their macroseismic intensities, and the empirical distribution of spectral intensities from Lushan and Wenchuan Earthquakes under the Chinese scale is almost identical with those under the European scale.     

A theoretical study of correlation and earthquake magnitude based models between scaled energy on two source displacement

The correlation of the scaled energy,e = Es/ Mo, versus earthquake magnitude, Ms, is studied based on two models： （1） Model 1 based on the use of the time function of the average displacements, with a ω^-2 source spectrum, across a fault plane; and （2） Model 2 based on the use of the time function of the average displacements, with a ω^-3 source spectrum, across a fault plane. For the second model, there are two cases： （a） As ζ= T, where r is the rise time and T the rupture time, lg（e） - -Ms; and （b） As ζ 〈〈 T, lg（e）- -（1/2）Ms. The second model leads to a negative value of e. This means that Model 2 cannot work for studying the present problem. The results obtained from Model 1 suggest that the source model is a factor, yet not a unique one, in controlling the correlation of e versus Ms.     

Focal Mechanism Solutions of the 2008 Wenchuan earthquake sequence from P-wave polarities and SH/P amplitude ratios:new results and implications

Abstract The 2008 Wenchuan earthquake, a major intraplate earthquake with Mw 7.9, occurred on the slowly deforming Longmenshan fault. To better understand the causes of this devastating earthquake, we need knowledge of the regional stress field and the underlying geodynamic processes. Here, we determine focal mechanism solutions （FMSs） of the 2008 Wenchuan earthquake sequence （WES） using both P-wave first-motion polarity data and SH/P amplitude ratio （AR） data. As P-wave polarities are more reliable information, they are given priority over SH/ PAR, the latter of which are used only when the former has loose constraint on the FMSs. We collect data from three categories： （1） permanent stations deployed by the China Earthquake Administration （CEA）; （2） the Western Sichuan Passive Seismic Array （WSPSA） deployed by Institute of Geology, CEA; （3） global stations from Incorporated Research Institutions for Seismology. Finally, 129 events with magnitude over Ms 4.0 in the 2008 WES are identified to have well-constrained FMSs. Among them, 83 are well constrained by P-wave polarities only as shown by Cai et al. （Earthq Sci 24（1）：115-125,2011）, and the rest of which are newly constrained by incorporating SH/P AR. Based on the spatial distribution and FMSs of the WES, we draw following conclusions： （1） the principle compressional directions of most FMSs of the WES are subhorizontal, generally in agreement with the conclusion given by Cai et al. （2011） but with a few modifications that the compressional directions are WNW-ESE around Wenchuan and ENE-WSW around Qingchuan, respectively. The subhorizontal compressional direction along the Longmenshan fault from SW to NE seems to have a leftlateral rotation, which agrees well with regional stress field inverted by former researchers （e.g., Xu et al., Acta Seismol Sin 30（5）, 1987; Acta Geophys Sin 32（6）, 1989; Cui et al., Seismol Geol 27（2）：234-242, 2005）; （2） the FMSs of the events not only reflected the regional stress state of the Longmenshan region, but also were obviously controlled by the faults to some extent, which was pointed out by Cai et al. （2011） and Yi et al. （Chin J Geophys 55（4）：1213-1227, 2012）; （3） while the 2008 Wenchuan earthquake and some of its strong aftershocks released most of the elastic energy accumulated on the Longmen- shan fault, some other aftershocks seem to occur just for releasing the elastic energy promptly created by the 2008 Wenchuan earthquake and some of its strong aftershocks. （4） Our results further suggest that the Longmenshan fault from Wenchuan to Beichuan was nearly fully destroyed by the 2008 Wenchuan earthquake and accordingly propose that there is less probability for great earthquakes in the middle part of the Longmenshan fault in the near future, although there might be a barrier to the southwest of Wenchuan and it is needed to pay some attention on it in the near future.     

汶川8．0级特大地震汉源县震害特点与烈度异常成因探讨

本文在汉源县地震灾害调查的基础上,通过收集前人在该区已取得的成果资料,采用综合分析、类比等方法,对汉源县地震烈度异常区震害及造成异常的原因进行了分析探讨。汉源县建筑物破坏形式以水平地震作用造成的“X”型剪切破坏为主,竖向地震造成的破坏较少。地震波传播路径、场地土层条件和地形地貌条件三者的耦合是造成汉源县烈度异常的主要原因,其中场地土层中占优势厚度的砾石层是造成场地地震波放大的主要内因。本文取得的初步结论为更加深入地科学研究本区地震烈度异常的原因提供了新思路。     

新搜集到的《嘉陵江日报》关于叠溪地震的报道(根据原件电子扫描稿校对)

民国期间的《嘉陵江日报》对1933年发生在四川省茂县叠溪的地震进行过系列报道,其中根据地质学家常隆庆先生的信函对于该地震考察情况的报道尤为详尽。过去,这些资料从未被地震历史资料的文献收录,鲜为地震工作者所知,具有较高的史料价值。整理转载这6000余字的资料,并作了注释和说明,相信能够为今后的地震史料研究提供参考。     

Processes of the displacement field change of the 2009 April 6 M_W6.3 L’Aquila earthquake using persistent scatterer and small baseline methods

Using a time series method that combines both the persistent scatterer and small baseline approaches, we analyzed 9 scenes Envisat ASAR data over the L＇Aquila earthquake, and obtained a Shocke＇s displacement field and its evolution processes. The results show that： （1） Envisat ASAR clearly detected the whole processes of displacement field of the L＇Aquila earthquake, and distinct variations at different stages of the displacement field. （2） Preseismic creep displacement → displacement mutation when faulting → constantly slowed down after the earthquake. （3） The area of the strongest deformation and ground rupture was a low-lying oval depression region to the southeast. Surface faulting within a zone of about 22 km× 14 km, with an orientation of 135°, occurred along the NW-striking and SW-dipping Paganica-S. Demetrio normal fault. （4） In analyzing an area of about 54 km x 59 km, bounded by north-south axis to the epicenter, the displacement field has significant characteristics of a watershed： westward of the epicenter shows uplift with maximum of 130 mm in line-of-sight （LOS）, and east of the epicenter was a region with 220 mm of maximum subsidence in the LOS, concentrating on the rupture zone, the majority of which formed in the course of faulting and subsequence.     

天津近海海域隐伏断裂地震危险性评价1

Based on faults surveying and research data in the Tianjin offshore areas, through studying tectonic structure, Quaternary activity, deep structure, stress and strain fields and seismicity in the Tianjin offshore areas, the activity and tectonic features of the faults are determined synthetically. Using seismo-geological data, and the historical and modern seismicity data, the frequency-magnitude relationship model normalized by 500a is established and based on the relationship between the upper limit of maximum magnitude Mu and at/b, the maximum magnitudes of the sea section of the Haihe river fault and the Haiyi fault are calculated. Then Poisson probability model is adopted and the quantitative parameters, such as the maximum magnitude, occurrence probability, recurrence cycle of the faults in the south Tianjin offshore areas in the coming 50 - 200a, are calculated.     

断层活动协调比计算方法的改进及震例分析

According to the concept of the fault motion coordination ratio （FCR） , this paper discusses the effect of the starting point on the result of FCR calculation and puts forward the calculation method for FCR using the sliding window and the index for evaluating the dispersion. Earthquake cases analysis shows that at the Lijiang site across the fault： the FCR value varied greatly and its distribution was dispersive before the Ms7.0 Lijiang earthquake, while the value was stable and less dispersive after the earthquake, which reflects the strain accumulation of the fault during the seismogenic process and the poor movement coordination between the motion of the three components. After the earthquake, the fault was in a free activity state, the accumulated strain energy released, and the movement of the three components was coordinated mutually. At present, FCR dispersion of Lijiang is at a low value, and fault strain accumulation is at a low level.     

“霍山窗”地区现代构造应力场研究

The focal mechanisms of 62 moderate-small earthquakes since 1980 in the ＂ Huoshan seismic window＂ region are calculated with the method developed recently by Snoke, combining the use of the first motion of P, SV and SH waves with their amplitude ratios. Based on these abundant focal mechanisms, the mean tectonic stress field in the ＂Huoshan seismic window＂ region is inverted with the average stress tensor method, and the result shows that the ＂Huoshan seismic window＂ region is horizontally compressed in the near EW direction and horizontally dilated in the near NS direction, which is in accord with statistical results of focal mechanism parameters. We estimate the difference （also referred to as consistency parameter 0） between the force axis direction of the focal mechanism solution and the mean stress tensor, then further analyze the variation characteristics of 0 versus time, and the relationship with moderately strong earthquakes in the east China region. The result indicates that 0 in the ＂ Huoshan seismic window＂ region is in good correspondence with moderately strong earthquakes in the East China region. When 0 is lower than the mean value, corresponding moderately strong earthquakes may occur in the East China region.     

Gravity observations along the eastern margin of the Tibetan Plateau and an application to the Lushan Ms7.0 earthquake

This paper introduces relative and absolute gravity change observations in the eastern portion of the Tibetan Plateau. We analyze and discuss a change that occurred in 2010 in the gravity along the eastern margin of the plateau and the relationship between this change and the 2013 Lushan Ms7.0 earthquake. Our results show that： （1） before the Lushan MsT.0 earthquake, gravity anomalies along the eastern margin of the Tibetan Plateau changed drastically. The Lushan earthquake occurred at the bend of the high gradient zone of gravity var- iation along the southern edge of the Longmenshan fault zone. （2） The 2013 Lushan earthquake occurred less than 100 km away from the epicenter of the 2008 Wenchuan earthquake. Lushan and Wenchuan are located at the center of a four- quadrant section with different gravity anomalies, which may suggest that restoration after the Wenchuan earthquake may have played a role in causing the Lushan earthquake. （3） A medium-term prediction based on changes in gravity anoma- lies was made before the Lushan Ms7.0 earthquake, in par- ticular, a prediction of epicenter location.     

An explanation of the large PGA value of 2013 Ms7.0 Lushan earthquake at 51BXD station through topographic analysis

In Ms7.0 Lushan earthquake, a large amount of strong ground motion recordings were collected. In this paper, we analyze the recordings carefully. The abnor- mality of ground motion recordings is identified through a log linear regression. In the station of 51BXD, the PGA value has exceeded 1 g, which is the biggest peak ground acceleration （PGA） value obtained from all recordings in this earthquake. The log linear relation shows the PGA value in this station is abnormally large. As this station is located on the footage of a hill, the topographic amplifi- cation factor is explored in order to explain this abnor- mality. Through 3D numerical modeling using spectral element method with transmitting boundary conditions, the amplification factor is quantized. In this station, the topo- graphic amplification is highly polarized in the direction of East-West which agrees with the empirical recordings. This research result suggests us in future directionality of topographic amplification should be considered in the aseismic design.     

Failure time variation derived from R-S relation： the role of the static stress perturbation

In general, earthquake cycle related to earthquake faulting could include four major processes which could be described by(1) fault locking,(2) self-acceleration or nucleation(possible foreshocks),(3) coseismic slip, and(4) post-stress relaxation and afterslip. A sudden static stress change/perturbation in the surrounding crust can advance/ delay the fault instability or failure time and modify earthquake rates. Based on a simple one-dimensional spring-slider block model with the combination of rate-and-statedependent friction relation, in this study, we have approximately derived the simple analytical solutions of clock advance/delay of fault failures caused by a sudden static Coulomb stress change applied in the different temporal evolution periods during an earthquake faulting. The results have been used in the physics-based explanation of delayed characteristic earthquake in Parkfield region, California, in which the next characteristic earthquake of M 6.0 after 1966 occurred in 2004 instead of around 1988 according to its characteristic return period of 22 years. At the same time, the analytical solutions also indicate that the time advance/ delay in Coulomb stress change derived by the dislocation model has a certain limitation and fundamental flaw. Furthermore, we discussed the essential difference between rate- and state-variable constitutive(R–S) model and Coulomb stress model used commonly in current earthquake triggering study, and demonstrated that, in fact, the Coulomb stress model could be involved in the R–S model. The results, we have obtained in this study, could be used in the development of time-dependent fault interaction model and the probability calculation related to the time-dependent and renewal earthquake prediction model.     

1992～2011年四川省地震局学术论文统计分析

以中国知网期刊全文数据库（CNKI）和维普科技期刊全文数据库为基础,统计分析了四川省地震局1992～2011年期间发表的地震期刊中文学术论文的情况。四川省地震局20年间共计发表论文1171篇,平均每年产出地震期刊学术论文58．5篇,以四川省地震局为第一署名单位论文925篇,平均每年46篇,论文数量逐年分布呈上升趋势,且与强震活动相关;论文刊载于74种中文科技期刊,总频次大于20次的高频期刊有8个;论文基金资助率逐年增加。论文统计结果显示四川地震科技研究工作一直在持续进步,在科技发展的主动性和创新性以及基层台站科研方面需进一步加强。     

汶川地震对水电工程的影响

2008年5月12日四川省汶川县发生8.0级地震,这是自1976年7月28日唐山7.8级地震后,中国发生的最大地震之一,给地震灾区人们带来深重的灾难。灾区大中型水电工程经受了远超工程设防标准的强震考验,大坝无一溃决,结构整体安全、稳定。部分建筑设施受到地震滑坡、崩塌等次生灾害影响,经及时抢修维护,均在震后迅速恢复发电及生产。灾区外大中型水电工程经受了强震考验,大坝完好无损,工程建设及运行正常有序。"5.12"汶川地震不是水库蓄水造成的,水库蓄水也不可能诱发"5.12"汶川地震。必须加强对地震机理、机制、地球板块运动、震源活动规律及趋势的探讨、认识与科学研究,逐步提高对地震的预警与预报,预防地震及地质灾害的水平,贯彻执行《防震减灾法》及有关律条例,复核、修订城镇工业、民用建筑和水利、水电工程抗震设防烈度及规范;结合生态环境修复,切实做好灾后重建工作。保护人民生命和财产安全,保障社会主义建设顺利进行,坚持可持续发展,构建和谐社会,确保长治久安。     

A study on the uncertainties of the centroid depth of the 2013 Lushan earthquake from teleseimic body wave data

Centroid depth of earthquakes is essential for seismic hazard mitigation. But, various studies provided different solutions for the centroid depth of the damaging 2013 Lushan earthquake, thus hindering further studies of the earthquake processes. To resolve its centroid depth and assess the uncertainties, we apply the teleseismic cut and paste method to invert for centroid depth with teleseismic body waves in the epicentral distance of 300-90~. We performed the inversion for P waves only as well the case of both P and SH waves and found that both cases lead to depth solutions with difference less than 0.5 km. We also investigated the effects on depth inversion from azimuth gap of seismic stations, source duration, and comer fre- quency of filter. These various tests show that even azi- muthal distribution of seismic stations is helpful for accurate depth inversion. It is also found that estimate of centroid depth is sensitive to source duration. Moreover, the depth is biased to larger values when corner frequency of low-pass filter is very low. The uncertainty in the velocity model can also generate some error in the depth estimation （- 1.0 km）.With all the above factors consid- ered, the centroid depth of Lushan earthquake is proposed to be around 12 km, with uncertainty about 2 km.