Nemuro-Oki earthquake of June 17, 1973: A lithospheric rebound at the upper half of the interface

P-wave first motions, radiation patterns and amplitudes of long-period surface waves, relocated aftershock distributions, leveling and tsunami data indicate that the 1973 Nemuro-Oki earthquake is caused by a low-angle thrust-faulting, representing a rebound at the upper 50 km of the interface between the continental and oceanic lithospheres. Rebound, most likely aseismic, at depths below 50 km, is suggested to take place in the near future from a comparison of recent geologic crustal deformation with pre-seismic and co-seismic data. The estimated seismic moment is about $\text{1}\text{3}\text{–}\text{1}\text{4}$ of that for the neighboring great earthquakes. The macro-seismic data suggest that the 1973 earthquake is smaller than the 1894 Nemuro-Oki earthquake, the last great earthquake in this region.The 1973 earthquake had been predicted on the basis of a seismic gap. Although the prediction was successful as to the location and nature of the faulting and partly as to the occurrence time, it is smaller than the predicted one. A part of the seismic gap may still remain. The difference between the observed seismic slip (1.6 m) and that predicted on the basis of the pre-seismic crustal deformation (3.0 m) indicates either (1) the 1973 earthquake relieved only a part of the strain accumulated in the upper 50 km, or (2) a significant amount of aseismic slip took place on the seismic fault and completely relieved the accumulated strain in the focal region of the 1973 earthquake. If the former is the case, the remaining strain, not only in the focal region, but also in the remaining seismic gap adjoining it, may be relieved in a larger earthquake in the future.The source parameters obtained are as follows: fault plane, dip direction = N40°W, dip angle = 27°; seismic moment = 6.7 · 10²⁷ dyn cm; average slip dislocation, 1.6 m in N63°W direction; stress drop = 35 bars. In these calculations, the fault dimension and the rigidity are assumed to be 100 · 60 km² and 7.0 · 10¹¹ dyn/cm², respectively.     

Ionospheric perturbations before Pu’er earthquake observed on DEMETER

A destructive shock with magnitude of 6.3 occurred on 2 June 2007 at 21h34min56s UT in Pu’er region (23.0°N, 101.1°E), Yunnan Province, China. The data from DEMETER satellite during the period from 23 May to 2 June, i.e., ten days before the earthquake and one day just on the day of earthquake occurrence, were analyzed. Among the 284 orbits of DEMETER during the period, 29 orbits with the trace passing through the region within 1 888 km from the epicenter were selected to be studied. Seven anomalous events were found on the dataset of the seven orbits among the 29 ones. There existed synchronous perturbations on the variations of the spectrogram of the electric field and the variations of the density and temperature of the ions and electron, in contrast with the variations of its surround- ing area. And five events appeared in the space within 1 888 km from the epicenter while the other two were out of the studied area. Electrostatic turbulences were also recorded with the synchronous perturbations with that in the electron density and ions density in plasma in the region near the epicenter in the five events, which seems to sug- gest that there be some physical relation between these events and the preparation processes of Pu’er earthquake.     

Earthquake responses in the high-latitude ionosphere

The data, obtained using the methods of partial reflections and ionosphere vertical sounding on the Kola Peninsula and in Scandinavia, at Tumannyi (69.0° N, 35.7° E) and Sodankyla (67.37°N, 26.63°E) observatories, have been analyzed in order to detect earthquake responses. The strong earthquakes have been considered: one earthquake with a magnitude of 7.7 occurred at 0819:25 UT on July 17, 2006, on the western coast of Indonesia (9.33° S, 107.26° E), and another earthquake with a magnitude of 6.2 occurred 2253:59 UT on May 26, 2006, on Yava (7.94° S, 110.32° E). These earthquakes, the epicenters of which were located in the same region and at identical depths (10 km), were observed under quiet conditions in the geomagnetic field (ΣK _p = 5.7 and 6.3) and during small solar flares. The response of the ionosphere to these flares was mainly observed in the parameters of the lower ionosphere in the D and E regions. It has been found out that the period of variations in the ordinary component of the partially reflected signal at altitudes of the E region increased before the earthquake that occurred on July 17, 2006. The f _min variations at Sodankyla observatory started 20 h before the earthquake. The periods of these variations were 3–6 h. The same periods were found in the variations in other ionospheric parameters (foEs and h’Es). The variations in the ordinary component of partially reflected signals with periods of 2–5 hours were observed on the day of another earthquake (May 26, 2006). Internal gravity waves with periods of several hours, which can be related to the earthquakes, were detected in the amplitude spectra of the ordinary component of partially reflected signals and in other parameters in the lower ionosphere.     

利用DEMETER卫星数据分析汶川地震前的电离层异常

DEMETER卫星由法国于2004年6月发射,高度约为665 km,圆轨道且经过极区,其主要目的是研究与地震有关的电离层异常.2008年5月12日6∶28（世界时）在四川省汶川县（31.02°N, 103.37°E）发生了Ms8.0级的特大地震.本文分析了在震前一周内DEMETER卫星观测到的电离层异常.这些异常主要包括：（1）卫星在震前4天和5天经过震中附近时,记录到电子浓度、电子温度和氧离子浓度出现剧烈变化（变化率均超过20%）;（2）在震前3天和7天,卫星探测到增强的电磁辐射.同时利用ELF频段（0~1250 Hz）的三分量电场和磁场数据,对异常电磁辐射作了进一步的分析（计算了电场功率谱、磁场功率谱、波矢方向和极化特性等）.另外,考查了观测期间的地磁指数,由此得出观测到的异常现象不大可能由地磁活动引起.最后,对震前出现的这些电离层异常的可能的物理机理进行了讨论.     

Dispersive O+ conics observed in the plasma-sheet boundary layer with CRRES/LOMICS during a magnetic storm

We present initial results from the Low-energy magnetospheric ion composition sensor (LOMICS) on the Combined release and radiation effects satellite (CRRES) together with electron, magnetic field, and electric field wave data. LOMICS measures all important magnetospheric ion species (H⁺, He⁺⁺, He⁺, O⁺⁺, O⁺) simultaneously in the energy range 60 eV to 45 keV, as well as their pitch-angle distributions, within the time resolution afforded by the spacecraft spin period of 30 s. During the geomagnetic storm of 9 July 1991, over a period of 42 min (0734 UT to 0816 UT) the LOMICS ion mass spectrometer observed an apparent O⁺ conic flowing away from the southern hemisphere with a bulk velocity that decreased exponentially with time from 300 km/s to 50 km/s, while its temperature also decreased exponentially from 700 to 5 eV. At the onset of the O⁺ conic, intense low-frequency electromagnetic wave activity and strong pitch-angle scattering were also observed. At the time of the observations the CRRES spacecraft was inbound at L\approx7.5 near dusk, magnetic local time (MLT), and at a magnetic latitude of -23°. Our analysis using several CRRES instruments suggests that the spacecraft was skimming along the plasma sheet boundary layer (PSBL) when the upward-flowing ion conic arrived. The conic appears to have evolved in time, both slowing and cooling, due to wave-particle interactions. We are unable to conclude whether the conic was causally associated with spatial structures of the PSBL or the central plasma sheet.     

Focal process of the great Chilean earthquake May 22, 1960

Long-period strain seismogram recorded at Pasadena is used to determine the focal process of the 1960 Chilean earthquake. Synthetic seismograms computed for various fault models are matched with the observed strain seismogram to determine the fault parameters. A low-angle (～ 10°) thrust model with rupture length of 800 km and rupture velocity of 3.5 km/sec is consistent with the observed Rayleigh/Love wave ratio and the radiation asymmetry. A seismic moment of 2.7 · 10³⁰ dyn · cm is obtained for the main shock. This value, together with the estimated fault area of 1.6 · 10⁵ km², gives an average dislocation of 24 m. The strain seismogram clearly shows unusually long-period (300–600 sec) wave arriving at the P time of a large foreshock which occurred about 15 minutes before the main shock, suggesting a large slow deformation in the epicentral area prior to the major failure. A simple dislocation model shows that a dislocation of 30 m, having a time constant of 300–600 sec, over a fault plane of 800 × 200 km² is required to explain this precursory displacement. The entire focal process may be envisaged in terms of a large-scale deformation which started rather gradually and eventually triggered the foreshocks and the “main” shock. This mechanism may explain the large premonitory deformations documented, but not recorded instrumentally, for several Japanese earthquakes. The moments of the main shock and the precursor add to 6 · 10³⁰ dyn · cm which is large enough to affect the earth's polar motion.     

A surface-wave investigation of the rupture mechanism of the Gobi-Altai (December 4, 1957) earthquake

Long-period records of multiple Love waves from the 1957 earthquake in Mongolia (M_S = 8.0) at Pasadena are analysed and compared to synthetic seismograms, generated by the method of Kanamori. A fit in the time domain shows that the records are not consistent with the previous solution, achieved through a frequency-domain analysis of directivity by Ben-Menahem and Toksöz. The solution asks for a shorter rupture of 270 km at a velocity of 3.5 km/s. The focal parameters are constrained by updating all the reported first motion and are found to be: strike = 103°, dip = 53°, slip = 32°. A seismic moment of 1.8 · 10²⁸ dyn · cm is obtained. These figures are also consistent with a time-domain analysis of Love waves at Palisades and Strasbourg, and of Rayleigh waves at Pasadena, with a directivity study of Love waves at Pasadena, and with static deformation and isoseismal data. A discussion is given of the relation between moment, magnitude and rupture area, and a comparison is made with other events in the same region: it is concluded that this earthquake does not exhibit an “intra-plate” behavior, but rather compares better with “inter-plate” events, such as the great Assam earthquake.     

基于Sentinel-1A数据反演漾濞MS6.4地震的同震形变场及断层几何参数

2021年5月21日晚21时48分,云南省大理州漾濞县(震中:25.67°N,99.87°E)发生M_S6.4地震,震源深度8 km。为快速获得此次地震同震形变场及断层几何参数,研究该次地震的发震构造等,文章基于震前、震后的sentinel-1A卫星升降轨SAR数据进行二轨法差分雷达干涉测量(DInSAR),并基于Okada弹性半空间位错模型反演断层几何参数。研究结果如下:(1)此次地震造成的同震形变场长约19 km,宽约20 km;(2)升轨雷达视线向最大形变约为8.2 cm,降轨雷达视线向最大形变约为8.7 cm;(3)地震断层走向为313.7°,倾角为87°,滑动角为175°,为右旋走滑型断层,最大滑动量为0.79 m,反演得出的地震矩为1.48×10~(18) N·m,矩震级为M_W6.1。在川滇块体向南挤出的构造背景下,块体西边界的维西—乔后断裂、红河断裂发生右旋走滑,本次地震便是维西—乔后断裂南段分支断裂右旋走滑活动的体现。     

运用InSAR技术解算四川长宁MS6.0地震三维形变场及时间序列形变分析

运用升、降轨Sentinel-1A 卫星的差分干涉影像,获取2019-06-17四川宜宾长宁 MS6．0地震的三维同震形变场.在此基础上,以升降轨同震形变数据为约束条件,基于 Okada弹性半空间位错模型反演得到发震断层符合走滑和逆冲特征,断层破裂尺度约为15km×20km,断层滑动角为 44．37°,断层倾角为56．42°,震源深度约为10．2km,矩震级为 M W5．8.最后,采用SBAS-InSAR 技术获取该地区2019-04-05至2019-08-03各时间段的累计形变,结果认为该区域在震前近场形变波动较小,震后一段时间累积形变增长,分析原因可能是余震分布使得地表变化处于不稳定状态.通过与已有研究文献的比较和对该区域断层构造的分析,推测此次长宁地震发震断层由反演出的断层滑动引起,滑动面上缘接近地表,主震引起的次级断层活动触发短期内强余震频发.     

Precursory Activation of Seismicity in Advance of the Kobe, 1995, M = 7.2 Earthquake

—A succession of precursory changes of seismicity characteristic to earthquakes of magnitude 7.0–7.5 occurred in advance of the Kobe 1995, M = 7.2, earthquake. Using the Japan Meteorological Agency (JMA) regional catalog of earthquakes, the M8 prediction algorithm (Keilies-Borko and Kossobokov, 1987) recognizes the time of increased probability, TIP, for an earthquake with magnitude 7.0–7.5 from July 1991 through June 1996. The prediction is limited to a circle of 280-km radius centered at 33.5°N, 133.75°E. The broad area of intermediate-term precursory rise of activity encompasses a 175 by 175-km square, where the sequence of earthquakes exhibited a specific intermittent behavior. The square is outlined as the second-approximation reduced area of alarm by the "Mendocino Scenario" algorithm, MSc (Kossobokov et al., 1990). Moreover, since the M8 alarm starts, there were no swarms recorded except the one on 9–26 Nov. 1994, located at 34.9°N, 135.4°E. Time, location, and magnitude of the 1995 Kobe earthquake fulfill the M8-MSc predictions. Its aftershock zone ruptured the 54-km segment of the fault zone marked by the swarm, directly in the corner of the reduced alarm area. The Kobe 1995 epicenter is less than 50 km from the swarm and it coincides with the epicenter of the M 3.5 foreshock which took place 11 hours in advance.     

Source process and tectonic implications of the Spanish deep-focus earthquake of March 29, 1954

The source process of the deep-focus Spanish earthquake of March 29, 1954 (m_b = 7.1, h = 630 km) has been studied by using seismograms recorded at teleseismic distances. Because of its unusual location, this earthquake is considered to be one of the most important earthquakes that merit detailed studies. Long-period body-wave records reveal that the earthquake is a complicated multiple event whose wave form is quite different from that of usual deep earthquakes. The total duration of P phases at teleseismic distances is as long as 40 s. This long duration may explain the considerable property damage in Granada and Malaga, Spain, which is rather rare for deep earthquakes. Using the azimuthal distribution of the differences between the arrival times of the first, the second and later P phases, the hypocenters of the later events are determined with respect to the first event. The focus of the second event is located on the vertical nodal plane of the first shock suggesting that this vertical plane is the fault plane. This fault plane which strikes in N2°E and dips 89.1°E defines a nearly vertical dip-slip fault, the block to the west moving downwards. The time interval and spatial separation between the first and the second events are 4.3 s and 19 km respectively, giving an apparent rupture velocity of 4.3 km/s which is about 74% of the S-wave velocity at the source. A third event occurred about 8.8 s after the first event and about 35.6 km from it. At least six to ten events can be identified during the whole sequence. The mechanism of some of the later events, however, seems to differ from the first two events. Synthetic seismograms are generated by superposition of a number of point sources and are matched with the observed signals to determine the seismic moment. The seismic moments of the later events are comparable to, or even larger than, that of the first. The total seismic moment is determined to be 7 · 10²⁷ dyn cm while the moments of the first and the second shocks are 2.1 · 10²⁶ dyn cm and 5.1 · 10²⁶ dyn cm, respectively. The earthquake may represent a series of fractures in a detached piece of the lithosphere which sank rapidly into the deep mantle preserving the heterogeneity of material property at shallow depths.     

Synthesis of accelerations of the 1976 Tangshan earthquake ( M S 7.8) in near- and far-field by using semi-empirical method

In this paper a heterogeneous fault model of the Tangshan earthquake is suggested, which consists of two southern sub-faults striking N30°E and two northern sub-faults striking N50°E. Total length of the main shock fault is 114 km and seismic moment is about 1.4 × 10²⁰ N·m. The epicentre of the main shock is located at the southern part, near the intersection of the two bands. Accelerations of two aftershocks (M _L 5.5, M _S 6.9) were used as empirical Green’s functions to synthesize the accelerations of the main shock in near and far field. A method that small events and main shock are considered not satisfying the similarity relationship in the improved empirical Green’s function is also applied in this paper. Peak values, duration and response spectra of synthesized accelerations in far field are in agreement with the observed records. The synthesized results in near field are also in agreement with the epicentral intensity distribution of the main shock. The results show that the peak acceleration of Tangshan earthquake in epicentral region exceeds 1.1 g. It is consistent with the peak accelerations recorded in some large earthquakes occurred in recent years. The research is supported by the National Science Foundation of China.     

Advanced sensitivity calibration of the Los Angeles strong motion array

Results are presented of recent sensitivity calibration of 76 accelerographs (SMA-1) of the Los Angeles Strong Motion Array. These have pendulum-like transducers and optical recording system. One characteristic of their design is off-axis sensitivity, which is magnified by transducer misalignment. A new calibration procedure was applied, which considers off-axis sensitivity and measures the angles of misalignment (φ and ψ), as well as the incident angle of the light beam onto the film (θ₀). These are required (1) for accurate estimation of sensitivity, and (2) for proper instrument correction of recorded accelerograms which considers also cross-axis sensitivity and misalignment. These effects are important near large acceleration peaks (approaching and exceeding 1g), e.g. like the ones recorded near the source of the 1994 Northridge earthquake (M_L=6·4). This earthquake was recorded by 65 stations of the Los Angeles Strong Motion Array, at epicentral distances from 2 to 85 km. Histograms showing distribution of the misalignment angles, light beam incidence angle θ₀ (for unloaded position) and the transducer sensitivities are presented. These indicate that the misalignment angles are typically 1–1·5°, but may also be 3–4°. Angle θ₀ (usually neglected), is mostly between ±8°, but may reach ±12°. Assuming θ₀=0 leads to systematically smaller values of the measured sensitivity (e.g. by ∼3% for θ₀=8° and ∼4% for θ₀=12°). Comparison of the newly measured sensitivities with those measured prior to installation (in 1979/1980), s^old, shows that, in general, the new values are systematically smaller. The difference is typically within 5 per cent, but in some cases is as large as 10 per cent. Other principal sources of the observed differences and their mechanisms are discussed. Those include long-term changes in the transducers (e.g. change of stiffness, reflected in changes of the natural frequency) and differences in the calibration procedure (e.g. errors associated with manual reading film records with tilt test data, and with transducer and instrument housing misalignment). The presented results may be considered typical of similar strong motion arrays worldwide. © 1998 John Wiley & Sons, Ltd.     

Seasonal and diurnal variation of electron and iron concentrations at the meteor heights above Arecibo

We report observations of seasonal and local time variation of the averaged electron and iron concentrations, as well as simultaneous measurements of the two species, above the Arecibo Observatory (18.35°N, 66.75°N), Puerto Rico. The average Fe profile between 21:00 and 24:00 LT has a single peak at about 85 km with the exception of the summer when an additional peak exists at about 95 km. The higher Fe peak in the summer is correlated with higher electron concentrations in this season. The three nights of simultaneous measurements of electron and iron concentrations show that narrow layers of Fe and electrons are well correlated. Comparison of the climatological and simultaneous Fe and electron data suggests that recombination of Fe⁺ plays an important role in determining the Fe profile in the upper part of the Fe layer. Above 93 km, the Fe concentration appears to increase after sunset if the electron concentration exceeds about 4000 electrons cm⁻³. The average rate of Fe production is about 0.1 atom cm⁻³ s⁻¹ for all seasons at 100 km in the early evening hours. A chemical model reveals that the concentration of Fe⁺ must be 50–80% of the total ionization over Arecibo for typical equinox conditions to explain the observed rate of Fe production. These high relative Fe⁺ concentrations are consistent with in situ observations that Fe⁺ is usually the dominant ion in sporadic E layers in the nighttime lower E region. This suggests that the source of Fe⁺ is provided by sporadic E layers descending over Arecibo after sunset. The Fe density between 80 and 85 km decreases during the night, for all seasons. This is attributed to the formation of stable molecular Fe species, such as FeOH, due to the increase in O₃ and decrease in atomic O and H during the night at these altitudes.     

GCITEM-IGGCAS: A new global coupled ionosphere–thermosphere-electrodynamics model

The Global Coupled Ionosphere–Thermosphere-Electrodynamics Model developed at Institute of Geology and Geophysics, Chinese Academy of Sciences (GCITEM-IGGCAS), is introduced in this paper. This new model self-consistently calculates the time-dependent three-dimensional (3-D) structures of the main thermospheric and ionospheric parameters in the height range from 90 to 600 km, including neutral number density of major species O₂, N₂, and O and minor species N(²D), N(⁴S), NO, He and Ar; ion number densities of O⁺ ,O₂⁺, N₂⁺, NO⁺, N⁺ and electron; neutral, electron and ion temperature; and neutral wind vectors. The mid- and low-latitude electric fields can also be self-consistently calculated. GCITEM-IGGCAS is a full 3-D code with 5° latitude by 7.5° longitude cells in a spherical geographical coordinate system, which bases on an altitude grid. We show two simulations in this paper: a March Equinox one and a June Solstice one, and compare their simulation results to MSIS00 and IRI2000 empirical model. GCITEM-IGGCAS can reproduce the main features of the thermosphere and ionosphere in both cases.     

The Campania-Lucania (southern Italy) earthquake of 23 November 1980

We present the main seismological results of our study of the Campania-Lucania earthquake of 23 November 1980. A complete set of far field and local data has been analysed. From long-period body waves data we determine the fault plane solution (φ₁ = 140°,δ₁ = 60°,φ₂ = 75°,δ₂ = 54°), a depth of 15 km and calculate a seismic moment of 6 × 10²⁵ dyne cm and a source duration of 6 s. From data of a local network deployed immediately after the event we determine aftershock locations: they are aligned in a direction NW-SE that fit extremely well with the focal solution determined above. We can choose as fault plane the plane striking 140° and dipping at 60° and the event is a normal event with a large component of left-lateral strike slip. The source area evaluated from this aftershock distribution 14 km × 40 km is quite suitable for an earthquake of a seismic moment of 6 × 10²⁵ dyne cm.     

Rupture process of the M L=4.1 earthquake in Huailai Basin on July 20, 1995

On July 20, 1995, an earthquake of M _L=4.1 occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, seismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocentres were precisely located. About 2 hours after the occurrence of the main shock, a smaller event of M _L=2.0 took place at 40.323°N, 115.447°E with a focal depth of 5.0 km, which is very close to the main shock. Using the M _L=2.0 earthquake as an empirical Green’s function, a regularization method was applied to retrieve the far-field source-time function (STF) of the main shock. Considering the records of HDSN are the type of velocity, to depress high frequency noise, we removed instrument response from the records of the two events, then integrated them to get displacement seismogram before applying the regularization method. From the 5 field stations, P phases in vertical direction which mostly are about 0.5 s in length were used. The STFs obtained from each seismic phases are in good agreement, showing that the M _L=4.1 earthquake consisted of two events. STFs from each station demonstrate an obvious “seismic Doppler effect”. Assuming the nodal plane striking 37° and dipping 40°, determined by using P wave first motion data and aftershock distribution, is the fault plane, through a trial and error method, the following results were drawn: Both of the events lasted about 0.1 s, the rupture length of the first one is 0.5 km, longer than the second one which is 0.3 km, and the rupture velocity of the first event is 5.0 km/s, larger than that of the second one which is about 3.0 km/s; the second event took place 0.06 s later than the first one; on the fault plane, the first event ruptured in the direction γ=140° measured clockwise from the strike of the fault, while the second event ruptured at γ=80°, the initial point of the second one locates at γ=−100° and 0.52 km from the beginning point of the first one. Using far-field ground displacement spectrum measurement method, the following source parameters about the M _L=4.1 earthquake were also reached: the scalar earthquake moment is 3.3×10¹³ N·m, stress drop 4.6 MPa, rupture radius 0.16 km. Contribution No. 99FE2022, Institute of Geophysics, China Seismological Bureau. This study is supported by the Chinese Joint Seismological Science Foundation (95-07-411).     

Altitude of the upper boundary of AAR based on observations of ion beams in inverted-V structures: A case study

Outflowing ion beams forming four successive inverted-V structures in the energy-time spectrograms of H+, He+, and O+ were observed at an altitude of 3.4 RE by Cluster satellites travelling above the auroral acceleration region (AAR) in the southern hemisphere on February 14, 2001. Energization by negative U-shaped potential structures in the AAR is believed to be responsible for the formation of these outflowing ion inverted-V structures. Thus, utilizing the different motion properties of the three ion species, the altitude of the upper boundary of the AAR is estimated to be ~11100 km. Moreover, based on multi-satellite observations, each of these U-shaped potential structures involved in this event crosses the latitudinal direction at ~0.4°–1° invariantlatitude (ILAT), moving poleward at an average speed of ~0.2° ILAT per minute, before disappearing at ~71.5° ILAT.     

Response of the ionosphere to the Baltic Sea earthquake of 21 September 2004

GPS observations of the European permanent network were used to identify seismo-ionospheric precursors of Baltic Sea earthquake of 21 September 2004. It is a very rare event for this region of Europe (magnitude of about 5.0). This value is the threshold for the occurrence of seismic effects in the ionosphere. In total electron content (TEC) data over the region of the earthquake, a specific ionospheric anomaly appeared one day before the earthquake was detected. The ionospheric variability had a positive sign with an enhancement of about 4–5 TECU (1 TECU = 10¹⁶ electrons/m²) relative to the non-disturbed state of the ionosphere. The anomaly had a duration of 4–5 hours in the day time. The special size of this anomaly was about 1000 km. The characteristic parameters of the anomaly show that it can be associated with ionospheric precursors of an earthquake.     

Prediction and characteristics of the 1975 eruption of Tolbachik volcano,Kamchatka

According to a long-term prediction. Tolbachik volcano was expected to erupt with a 0.7 probability, some time in the period 1964–1978. An eruption of Tolbachik commenced at 21.45 GMT on July 5, 1975. It took place on the southwestern Hank of the volcano at an altitude of 880 m a.s.l. about 18 km from the central crater. An earthquake swarm preceded it. The place and time of eruption were predicted three days belore it began on the basis of epicenter locations and characteristics of recorded seismic activity. During the period July 5–28 gases and incandescent magma were continuously ejected to a height of 2,000 m above ground level. Ash clouds rose to a height of 6 to 8 km, with a train of ash extending over a horizontal distance of 300 km. The velocity of jets from the crater was about 200 m/see. During the first days of the eruption the quantity of materials erupted and the eruption power amounted to 1.25 · 10⁵ kg/see and 2.1 · 10¹¹W, respectively. The vertical growth of the scoria cone was consistent with the lawH=2.15³√vt, where the time and height are expressed in seconds and meters, respectively. The mouth of the volcano conduit was estimated to be 12 m in diameter. Lava began to erupt at 22^h23^m on July 28. During the period July 5–31 about 3 · 10¹¹ kg of magmatic material, consisting of ash, scoria and lava, was erupted onto the earth’s surface. The energy released over the period of eruption accounted for 5 · 10¹⁷ J.