利用宏观烈度资料估计历史地震震源参数的研究概述

历史地震宏观烈度资料指史料中对地震破坏程度的记录,是研究历史地震最基本也是最重要的数据来源.由于有仪器记录地震波形的时间相对短暂,对历史地震震源参数的估算成为地震危险性分析的重要组成部分,前人尝试使用不同方法解决这一问题.本文以使用的工具和计算手段为依据,分类介绍了基于烈度资料估计震源参数的主要方法,对这些方法的优势、不足和应用范围做了概述,并展望了利用烈度资料估计震源参数研究的发展方向.     

Transfrontier macroseismic observations of the Ml = 5.4 earthquake of February 22, 2003 at Rambervillers,France

The 2003 Ml = 5.4 Rambervillers earthquake, north-east of France, is the largest seismic event recorded north of the Alps since the 1992 Ms = 5.3, I₀ = VII, Roermond earthquake, Netherlands. With a maximum macroseismic intensity of VI-VII EMS-98, the 2003 event was broadly felt to a distance of 300 km from the epicentre. It provides a unique opportunity to test and compare the different procedures used in France, Germany and Switzerland when evaluating macroseismic intensities. The main purpose of this paper is to present a common transfrontier macroseismic map based on the EMS-98 intensity scale. Maximum horizontal accelerations recorded in the area are compared to the intensity values, and we propose to use a differential technique to re-estimate the magnitude of the 1682 Remiremont, I₀ = VIII, earthquake, which occurred 40 km south of Rambervillers.     

Source mechanism and source parameters of May 28, 1998 earthquake,Egypt

On May 28, 1998, a moderate size earthquake of mb 5.5 occurred offshore the northwestern part of Egypt (latitude 31.45°N and longitude 27.64°E). It was widely felt in the northern part of Egypt. Being the largest well-recorded event in the area for which seismic data from the global digital network are available, it provides an excellent opportunity to study the tectonic process and present day stress field occurring along the offshore Egyptian coast. The source parameters of this event are determined using three different techniques: modeling of surface wave spectral amplitudes, regional waveform inversion, and teleseismic body waveform inversion. The results show a high-angle reverse fault mechanism generally trending NNW–SSE. The P-axis trends ENE–WSW consistently with the prevailed compression stress along the southeastern Hellenic arc and southwestern part of the Cyprean arc. This unexpected mechanism is most probably related to a positive inversion of the NW trending offshore normal faults and confirms an extension of the back thrusting effects towards the African margin. The estimated focal depth ranges from 22 to 25 km, indicating a lower crustal origin earthquake owing to deep-seated tectonics. The source time function indicates a single source with rise time and total rupture duration of 2 and 5 s, respectively. The seismic moment (M _o) and the moment magnitude (M _w) determined by the three techniques are 1.03 × 10¹⁷ Nm, 5.28; 1.24 × 10¹⁷ Nm, 5.33; and 1.68 × 10¹⁷ Nm, 5.42; respectively. The calculated fault radius, stress drop, and the average dislocation assuming a circular fault model are 7.2 km, 0.63 Mpa, and 0.11 m, respectively.     

Source model and the macroseismic effect of the 1991 Racha earthquake

Coherency of the source model of the 1991 Racha earthquake in the Greater Caucasus with different data types is analyzed. Authors, when interpreting macroseismic data, accept complex nature of macroseismic effects generation but, nevertheless, consider that its spatial distribution follows certain regularities. First time in the practice, method of evaluation of the mecroseismic material completeness is proposed based on the intensity attenuation along with distance. It is demonstrated the character of macroseismic intensity attenuation can be used for verification of the source model elements constructed based on other seismological data. Dependence of the macroseismic effect distribution on azimuth in near field of the 1991 Racha earthquake is recognized.     

Influence of fluids on deep crustal Jabalpur earthquake of 21, May 1997: Geophysical evidences

We present the geophysical evidences on the role of fluids for generation of the lower crustal Jabalpur earthquake (21 May 1997, mb 6.0, Mw 5.8), in the mid-continental fracture zone of the Indian Peninsular Shield. With a focal depth of 35 km, it indicates a high angled (< 62^∘ enclosed with maximum principal stress direction) reverse fault with small component of left-lateral strike slip in the lower crust. The Son-Narmada-Tapti (SONATA) magalineament, during the past two centuries, has experienced about 25 moderate to strong earthquakes; two of which namely the Son Valley (1927, M 6.5) and Jabalpur (21 May 1997) were disastrous. Historical earthquakes and recent earthquake swarms indicate a moderate to high seismicity in SONATA belt that is due to high strain accumulation, flexuring of the crust and neotectonic movements of the faults in the rift zones. By analyzing geophysical parameters such as Zero-Free air-based (ZFb) gravity anomalies (∼ −10 to –30 mGals), heat flow values (45–47 mWm⁻²), magneto-telluric values (1- Ohm m), strain rate (1.5 × 10⁻⁸) and failure stress conditions, we identify plausible causative factors for the occurrence of lower crustal earthquake in this region Fluids, due to dehydration of serpentinite in the lower crust, are suggested to be present in the earthquake source zone. The estimated pore-fluid factor for the Jabalpur earthquake (λ _v ) is 0.95. The diffusion of pore-pressure relaxation, represented as pressure perturbation generated by coseismic stress change was seen in the form of swarm activity two years prior to the Jabalpur earthquake. We suggest the existence of a deep pre-fractured zone with low shear stress (τ = 15–18 MPa) that indicates the presence of fluid filled fractured mafic material in the felsic crust, in critical state of unstable failure condition, and also fluid driven migration of swarm activity before the Jabalpur earthquake.     

Geological and macroseismic effects of the Muya, 1957 earthquake and palaeoearthquakes in Baikal region

Intensity of the Muya, 1957 earthquake is assessed in localities based on macroseismic data and in epicentral area based on effects in natural environment; it is analyzed how these assessments correspond to each other and to instrumental location of epicenter, hypocentral depth, and magnitude; it is evaluated, how seismodislocations of the Muya earthquake could serve as control of palaeoseismostructure parameters in this region. Spatial distribution of macroseismic effect confirms relatively deep source (20–22 km). Deep source agrees with anomalously short surface rupture length (not more than 25 km); only a part of the source exposed on the surface. Comparison with length of palaeoseismostructures shows that it is a regional feature. Epicentral intensity based on surface ruptures is assed X degrees in ESI2007 scale. Ignoring geological effects will underestimate epicentral intensity up to two degrees. Source mechanism with three sub-sources is in agreement with segmentation of surface ruptures. Sub-sources are of strike-slip type with small normal component; essential normal slip at surface is probably not representative for the source and is due to accommodation of strike-slip movement along with a system of sub-parallel en echelon ruptures under tension.     

Structural heterogeneity of the media and macroseismic field formation: Racha,April 29, 1991 earthquake

The points with normal, anomalously low, and anomalously high shaking intensities are recognized in the spatial distribution of macroseismic effects from the 1991 Racha earthquake, Greater Caucasus. Distribution of these points in the epicentral area is not random. Comparison between this distribution and the results of local tomography reveals that seismic wave velocities do not increase in the upper layers (from 0 to 3 km) beneath the points with anomalously high intensity, while a sharp increase in velocity is observed in the depth interval from 6 to 9 km. An original method of b-value mapping is suggested. Application of the method demonstrates that anomalously low intensities correlate to high b-values. This likely reflects higher intensity attenuation associated with higher b-value.     

Combined ESR and EISCAT observations of the dayside polar cap and auroral oval during the May 15, 1997 storm

The high-latitude ionospheric response to a major magnetic storm on May 15, 1997 is studied and different responses in the polar cap and the auroral oval are highlighted. Depletion of the F₂ region electron density occurred in both the polar cap and the auroral zone, but due to different physical processes. The increased recombination rate of O⁺ ions caused by a strong electric field played a crucial role in the auroral zone. The transport effect, however, especially the strong upward ion flow was also of great importance in the dayside polar cap. During the main phase and the beginning of the recovery phase soft particle precipitation in the polar cap showed a clear relation to the dynamic pressure of the solar wind, with a maximum cross-correlation coefficient of 0.63 at a time lag of 5 min.     

ERNE observations of energetic particles associated with Earth-directed coronal mass ejections in April and May, 1997

Two Earth-directed coronal mass ejections (CMEs), which were most effective in energetic (1–50 MeV) particle acceleration during the first 18 months since the Solar and Heliospheric Observatory (SOHO) launch, occurred on April 7 and May 12, 1997. In the analysis of these events we have deconvoluted the injection spectrum of energetic protons by using the method described by Anttila et al. In order to apply the method developed earlier for data of a rotating satellite (Geostationary Operational Environmental Satellites, GOES), we first had to develop a method to calculate the omnidirectional energetic particle intensities from the observations of Energetic and Relativistic Nuclei and Electrons (ERNE), which is an energetic particle detector onboard the three-axis stabilized SOHO spacecraft. The omnidirectional intensities are calculated by fitting an exponential pitch angle distribution from directional information of energetic protons observed by ERNE. The results of the analysis show that, compared to a much faster and more intensive CMEs observed during the previous solar maximum, the acceleration efficiency decreases fast when the shock propagates outward from the Sun. The particles injected at distances <0.5 AU from the Sun dominate the particle flux during the whole period, when the shock propagates to the site of the spacecraft. The main portion of particles injected by the shock during its propagation further outward from the Sun are trapped around the shock, and are seen as an intensity increase at the time of the shock passage.     

Assessing the intensity of the September 17, 2003, Khoito-Gol earthquake from instrumental and macroseismic data

The present study is devoted to the assessment of shaking intensities due to the September 17, 2003, earthquake occurring in the Khoito-Gol basin (southwestern flank of the Baikal Rift System). The instrumental and macroseismic data used here were acquired during an investigation into the impact of this seismic event. The highest intensity of shaking, VI, was recorded at Khoito-Gol, the human settlement that was the nearest to the instrumental epicenter. A nonuniform distribution of intensity was found for different directions from the epicenter. A scatter of as much as two intensity units was recorded at one settlement.     

GNSS观测的2021年5月22日玛多MS7.4地震同震位移及其约束反演的滑动破裂分布

2021年5月22日青海省果洛州玛多县发生M_W7.4地震,此次地震产生的地表破裂在空间上表现出明显的分段特征.本文基于不同来源的GNSS连续观测网数据获取了此次地震的精细三维同震形变场,结果显示:观测到的最大水平位移量达到280 mm,最大垂直形变量仅为25 mm,暗示此次地震的逆冲分量较小;此次地震具有较为明显的左旋走滑特征,同震形变基本对称,在NW-SE向的影响范围更广,该方向上水平同震形变大于3 mm的震中距范围超过500 km.进而,本文以余震精定位结果和GNSS观测到的三维同震形变场为约束,构建了地表破裂线为折线、倾角为85°、倾向西南的断层模型,反演了滑动破裂分布.结果显示:滑动破裂分布在震中两侧不均匀,均破裂到地表,破裂深度达到15 km左右,最大滑移量为4.73 m,计算的矩震级为M_W7.37.该结果与余震精定位结果具有很好的一致性,破裂的极值区正好位于早期余震空区,推测该余震空区未来的发震风险性较低.最后基于反演结果模拟计算了震中区域形变和应变场,结合应变值在断层地表迹线东南侧呈现挤压特征和已有的研究成果,推测此次地震增强了巴颜喀拉块体在东部地区挤压应力的积累特征,导致东部地区发震危险性增强,值得后续跟踪研究.     

Linking ground motion measurements and macroseismic observations in France: a case study based on accelerometric and macroseismic databases

Comparison between accelerometric and macroseismic observations is made for three M _w?=?4.5 earthquakes, which occurred in north-eastern France and south-western Germany in 2003 and 2004. Scalar and spectral instrumental parameters are processed from the accelerometric data recorded by nine accelerometric stations located between 29 and 180 km from the epicentres. Macroseismic data are based on French Internet reports. In addition to the single questionnaire intensity, analysis of the internal correlation between the encoded answers highlights four predominant fields of questions bearing different physical meanings: (1) “vibratory motions of small objects”, (2) “displacement and fall of objects”, (3) “acoustic noise” and (4) “personal feelings”. Best correlations between macroseismic and instrumental observations are obtained when the macroseismic parameters are averaged over 10-km-radius circles around each station. Macroseismic intensities predicted by published peak ground velocity (PGV)–intensity relationships agree with our observed intensities, contrary to those based on peak ground acceleration (PGA). Correlation between the macroseismic and instrumental data for intensities between II and V (EMS-98) is better for PGV than for PGA. Correlation with the response spectra exhibits clear frequency dependence for all macroseismic parameters. Horizontal and vertical components are significantly correlated with the macroseismic parameters between 1 and 10 Hz, a range corresponding to both natural frequencies of most buildings and high energy content in the seismic ground motion. Between 10 and 25 Hz, a clear lack of correlation between macroseismic and instrumental observations exists. It could be due to a combination of the decrease in the energy signal above 10 Hz, a high level of anthropogenic noise and an increase in variability in soil conditions. Above 25 Hz, the correlation coefficients between the acceleration response spectra and the macroseismic parameters are close to the PGA correlation level.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

On the influence of earthquake energy range on the estimated seismicity parameters before the magnitude 7.9 Kronotskii earthquake of 1997

Results are reported from a detailed study of central Kamchatka seismicity for the period 1962–1997 based on a modification of the traditional approach. The approach involves (a) a detailed structure of the seismic region that recognizes the Kronotskii and Shipunskii geoblocks and two further blocks, the continental slope, and the offshore portion, (b) a study of variations in the rate of M = 3.0–7.2 earthquakes and the amount of seismic energy released at depths of 0–50 and 51–100 km, (c) a study of seismicity variability, and (d) separate estimates of the recurrence of crust-mantle earthquakes (depths 0–50 km) and mantle events (51–100 km). As a result, apart from corroborating the fact of a quiescence preceding the December 5, 1997 Kronotskii earthquake (M 7.9), we also found that a relationship exists between its beginning and the position of the earthquake-generating region relative to the mainshock epicenter. The quiescence dominates the seismic process during the pre-mainshock period and is characterized by a decreased rate of earthquakes (the first feature) and a decreased amount of seismic energy release (the second feature). Based on the first feature, we found that the quiescence started in 1987 throughout the entire depth range (0–100 km) in both parts of the Kronotskii geoblock close to the rupture zone of the eponymous earthquake. As to the Shipunskii geoblock, which is farther from the rupture zone, the quiescence began in the mantle of the inner area first (1988) and somewhat later at depths of 0–50 km within the continental slope (1989). By the second feature, the quiescence began at shallower depths in the inner area of the Kronotskii geoblock at the same time and later on (a year later) in the mantle (1988). Under the continental slope of the trench in the Shipunskii geoblock the shallower quiescence also began in 1987, while it was 3 years late in the inner zone (1990) and involved the earthquake-generating earth volume at depths of 0–100 km. These data are identical with or sufficiently close to the estimate for the beginning of this quiescence using a circular area of radius 150 km that combines the Kronotskii and Shipunskii geoblocks by the RTL method (1990).     

2017年精河MS6.6地震的断层参数和破裂过程

利用远震资料、近场强震资料和合成孔径雷达干涉同震形变资料确定了2017年8月9日精河M_S6.6地震的断层面参数及震源破裂细节。为得到可靠的断层几何参数,发展了一套基于InSAR数据滑动分布反演的三维格点搜索流程,对本次地震断层面的走向、倾角和震源深度进行了格点搜索。结果显示,地震断层面走向为95°,倾角为47°,震源深度为14 km。基于搜索得到的断层模型进行破裂过程联合反演的结果显示:精河M_S6.6地震为一次单侧破裂事件,最大滑动量约为0.8 m,滑动区域集中在断层面上震源以西5—15 km,沿倾向15—25 km,破裂主要发生在10 km深度以下区域。断层面上的平均滑动角为106°。整个破裂过程释放的标量地震矩为3.6×1018 N·m,对应矩震级为M_W6.3。破裂过程持续约9 s,期间的破裂速度约为2.1—2.6 km/s。由于地震破裂主要集中在10 km以下,未来可能需要关注该区域0—10 km发生潜在地震的可能性。      

Pyroclastic surges and flows from the 8–10 May 1997 explosive eruption of Bezymianny volcano, Kamchatka, Russia

The 8-10 May 1997 eruption of Bezymianny volcano began with extrusion of a crystallized plug from the vent in the upper part of the dome. Progressive gravitational collapses of the plug caused decompression of highly crystalline magma in the upper conduit, leading at 13:12 local time on 9 May to a powerful, vertical Vulcanian explosion. The dense pyroclastic mixture collapsed in boil-over style to generate a pyroclastic surge which was focused toward the southeast by the steep-walled, 1956 horseshoe-shaped crater. This surge, with a temperature <200 °C, covered an elliptical area >30 km² with deposits as much as 30 cm thick and extending 7 km from the vent. The surge deposits comprised massive to vaguely laminated, gravelly sand (Md -1.2 to 3.7J sorting 1.2 to 3J) of poorly vesiculated andesite (mean density 1.82 g cm^-3; vesicularity 30 vol%; SiO₂ content ~58.0 wt%). The deposits, with a volume of 5-15᎒⁶ m³, became finer grained and better sorted with distance; the maximal diameter of juvenile clasts decreased from 46 to 4 cm. The transport and deposition of the surge over a snowy landscape generated extensive lahars which traveled >30 km. Immediately following the surge, semi-vesiculated block-and-ash flows were emplaced as far as 4.7 km from the vent. Over time the juvenile lava in clasts of these flows became progressively less crystallized, apparently more silicic (59.0 to 59.9 wt% SiO₂) and more vesiculated (density 1.64 to 1.12 g cm^-3; vesicularity 37 to 57 vol%). At this stage the eruption showed transitional behavior, with mass divided between collapsing fountain and buoyant column. The youngest pumice-and-ash flows were accompanied by a sustained sub-Plinian eruption column ~14 km high, from which platy fallout clasts were deposited (~59.7% SiO₂; density 1.09 g cm^-3; vesicularity 58 vol%). The explosive activity lasted about 37 min and produced a total of ~0.026 km³ dense rock equivalent of magma, with an average discharge of ~1.2᎒⁴ m³ s^-1. A lava flow ~200 m long terminated the eruption. The evolutionary succession of different eruptive styles during the explosive eruption was caused by vertical gradients in crystallization and volatile content of the conduit magma, which produced significant changes in the properties of the erupting mixture.     

Geological and macroseismic manifestations of the earthquake on October 16, 2011 in the Skovorodino region, Amur oblast

The epicentral zone and settlements that suffered from the M _S = 6.1 earthquake in the northwest Amur oblast are examined. Only secondary seismic dislocations were revealed and mapped in detail. The inspection of settlements and inhabitansts inquiry allowed estimation of the intensity of the macroseismic effect based on the MSK-64 scale. These forthwith primary factual data give an idea on the shaking intensity at the distant and nearest zones and precise location of the earthquake focus. The map of isoseists of the highest (7–8) intensity is drawn.     

Source parameters of bohai earthquake, July 18, 1969 and yongshan earthquake, May 11, 1974 determined by synthetic seismograms of teleseismic P waves

The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M _o=0.9×10¹⁹Nm,L=59.9km.V _R=3.5km/s, ∧ _R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M _o=1.3×10¹⁸Nm,L=48.8km,V _R=3km/s, ∧ _R =−10°, where φ is strike, δ dip angle, λ slip angle,M _o seismic moment,L rupture length,V _R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ _R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.