Effects of strong earthquakes in variations of electrical and meteorological parameters of the near-surface atmosphere in Kamchatka region

The diurnal variations in electrical (quasistatic electric field and electrical conductivity) and meteorological (temperature, pressure, relative humidity of the atmosphere, and wind speed) parameters, measured simultaneously before strong earthquakes in Kamchatka region (November 15, 2006, М = 8.3; January 13, 2007, М = 8.1; January 30, 2016, М = 7.2), are studied for the first time in detail. It is found that a successively anomalous increase in temperature, despite the negative regular trend in these winter months, was observed in the period of six–seven days before the occurrences of earthquakes. An anomalous temperature increase led to the formation of “winter thunderstorm” conditions in the near-surface atmosphere of Kamchatka region, which was manifested in the appearance of an anomalous, type 2 electrical signal, the amplification of and intensive variations in electrical conductivity, heavy precipitation (snow showers), high relative humidity of air, storm winds, and pressure changes. With the weak flow of natural heat radiation in this season, the observed dynamics of electric and meteorological processes can likely be explained by the appearance of an additional heat source of seismic nature.     

Spectra of atmospheric waves in the dynamo region of the ionosphere at Kamchatka

A spectral analysis of the diurnal variations in the geomagnetic field horizontal component, observed at Kamchatka and Barrow polar observatory in September–October 1999, has been performed. The complete set of oscillations of thermal tidal atmospheric waves with T = 24, 12, 8, and 4 h has been detected in the variation spectral power (Sq) at Kamchatka, and only the fundamental harmonic with T = 24 h has been distinguished at Barrow. The above periods vary in both directions relative to stable maximums during strong geomagnetic disturbances. The relative spectral intensity at subharmonics also vary toward the fundamental harmonic with a period of 24 h. In the frequency band 0.5–3 h (IGW periods), the maximal intensity in the background spectra is observed at T ～ 2 h and increases by an order of magnitude with increasing geomagnetic activity at both Kamchatka and Barrow. A day before earthquakes, the intensity of this maximum is below the rms background values, and the spectra widen toward the region of periods shorter than 2 h. A similar effect was previously observed in the power spectra of the diurnal variations in the quasistatic electric field and VLF noise, simultaneously measured in September–October 1999.     

Experience in registration of variations caused by strong earthquakes in the level and physicochemical parameters of ground waters in the piezometric wells: the case of Kamchatka

The hardware complex that was produced by OOO Polynom, Khabarovsk, for registration of the level, temperature, and electrical conductivity of ground water in wells and meteorological parameters (atmospheric pressure, air temperature) at a measurement frequency from 5 min to 1 h is described. The equipment is installed in the wells of Kamchatka and has been used for several years to register variations caused by earthquakes in the parameters of ground waters. Different variations in measured parameters of ground waters due to strong earthquakes of February 28, 2013, M _W = 6.8 and May 24, 2013, M _W = 8.3 are registered with this equipment in wells YuZ-5 and E-1. The registered variations and their systematization are described taking into account the mechanisms of a seismic impact on the state of the well–water-saturated rock system.     

Spatial relation between the tidal component of seismic noise and development zones of strong earthquakes (from long-term regular observations on the Kamchatka Peninsula)

The results of long-term studies of seismic noise before strong regional earthquakes on the Kamchatka Peninsula in 1992–2006 are presented. These results show that parameters of seismic noise variations caused by the tidal effect depend on the source position of the forthcoming earthquake. The reproducibility of the observed effects is demonstrated by the example of two strong deep (～200 km) subduction earthquakes with similar parameters that occurred on June 16, 2003 (M = 6.9), and on June 10, 2004 (M = 6.8). The physical mechanism of the synchronization of the tidal component extracted from high-frequency seismic noise with the wave of the gravitational tidal potential can be related to the possible development of near-surface dilatancy zones.     

Problem of the nature of the sunrise effect in diurnal variations in the electric field in Kamchatka: 1. Time variations in the electric field

The effect of sunrise in time variations in the electric field in the near-Earth atmosphere at the Kamchatka Paratunka observatory has been studied. Twenty-nine records under fair-weather conditions have been selected. It has been indicated that the estimated effect parameters—the times of the effect’s onset and field strength maximum relative to the sunrise time, as well as the ratio of the strength maximum to its value before sunrise and the effect duration—coincide with the previously published data. Thereby, the conclusion is confirmed that the sunrise effect in diurnal variations in the electric field in the near-Earth atmosphere is related to the turbulence and convection processes in the atmospheric boundary layer at a change in atmospheric temperature.     

Electric and electromagnetic processes in the near-Earth atmosphere before earthquakes on Kamchatka

Simultaneous records of the diurnal variations in the quasistatic electric field in the near-Earth atmosphere, fluxes of discrete electromagnetic pulses in the VLF band, source azimuths, narrowband filter output emission intensity at frequencies of 4.65 and 5.3 kHz, and time forms and spectra of VLF pulses have been analyzed. The anomalous behavior of these parameters in October 2002 and August 2004 with different time delay was accompanied by earthquakes near the southeastern coast of Kamchatka at distances of 250–400 km from the registration points. Based on the results of a fine frequency-time analysis of the broadband records of VLF signals, it has been indicated for the first time that discrete electromagnetic pulses observed in anomalous fluxes before earthquakes were signals of local thunderstorm processes.     

Anomalous solar-diurnal variations in cosmic rays related to crossing of the IMF sector boundaries by the Earth and the problem of earthquakes

Based on a comparison of the cases of a decrease in the ratio of A _n/A _μ (where A _n and A _μ are the amplitudes of the diurnal variations of the neutron and hard cosmic ray components) to the instants of the Earth crossing the neutral IMF, it has been indicated that the process of such crossing is most effective for stimulating large destructive earthquakes with a magnitude of M ≥ 6. The 11-year period in the cyclicity of the occurrence probability of the above earthquakes has been revealed.     

Rationale for applying short-period sensors with extended frequency response for recording strong earthquakes

The article proves the application of short-period sensors with extended frequency response as a reasonable alternative to broadband seismometers. We assessed uncertainties between magnitudes determined by short-period and broadband sensor data for earthquakes with M > 6 recorded at the Mikhnevo geophysical observatory in 2014. Data analysis consisted of body wave magnitude and surface wave magnitude estimates by standard seismological methods. Magnitude errors corresponded to the magnitude uncertainties in seismological catalogs.     

Determination of source parameters of strong earthquakes of Kamchatka,the Kurile Islands,and Japan from aftershock sequences

Aftershock sequences of some strong earthquakes of Kamchatka, the Kurile Islands, and Japan are examined. Such source parameters as the length L, along-dip width W, motion on fault D, and stress drop Δσ are determined from the aftershock sequences considered. The values of these parameters were obtained by the formal estimation of linear source parameters (lower bound estimates) and visually (upper bound estimates). The correlation dependences of the obtained parameters on the surface wave (M _S ) and seismic moment (M _W ) magnitudes are calculated.     

On the question of influence of remote earthquakes on seismicity

The responses of the Kamchatka earthquakes with a minimum energy class of completeness K ≥ 8.5 to 214 strong worldwide earthquakes with magnitudes M ≥ 7.5 and to 40 earthquakes with M ≥ 8 are studied. The analysis covers the time interval of 1963–2012. The distances from the sources of the strongest earthquakes to the center of the seismically active Kamchatka zone range from 600 to 16000 km. It is established that the remote earthquakes enhanced seismic activity in Kamchatka, at least in the cases when the dynamic strain was above 10^?6, which corresponds to the additional stresses of 10^?2 MPa, accelerations above 0.1 cm/s², and the periods of the surface waves of ～20 s. The response to the remote events gradually increased within a few days. The sensitivity of the response to the remote earthquakes varied in the course of time, which is identified on the intervals of a few dozens of years.     

Detailed Study of Time Variations in the Gutenberg–Richter <Emphasis Type="Italic">b</Emphasis>-value Based on Highly Accurate Seismic Observations at the Garm Prognostic Site,Tajikistan

The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log₁₀E^2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log₁₀E^2/3, and log₁₀N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log₁₀E^2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log₁₀N_i and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log₁₀N_i time series in the different magnitude ranges (M ≥ Mc_i) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔM_i) in the range of low magnitudes are known. It is concluded that using the parameter log₁₀N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.     

Discussion on the feature of strong earthquake: Orderly distribution in time,space and intensity before the Western Kunlun Mountain Pass <Emphasis Type="Italic">M</Emphasis>=8.1 earthquake

In the paper, the feature of strong earthquake orderly distribution in time, space and intensity before the Western Kunlun Mountain Pass M=8.1 earthquake is preliminarily studied. The modulation and triggering factors such as the earth rotation, earth tides are analyzed. The results show that: the giant earthquakes with the magnitude more than 8 occurred about every 24 years and the earthquakes with the magnitude more than 7 about every 7 years in Chinese mainland. The Western Kunlun Mountain M=8.1 earthquake exactly occurred at the expected time; The spatial distance show approximately the same distances between each two swarms. The earth rotation, earth tide, sun tide and sun magnetic field have played a role of modulation and triggering in the intensity. At last, the conditions for earthquake generation and occurrence are also discussed.     

Dynamics of Induced Seismicity during the Filling of the Nurek Reservoir

The seismological data in the area of induced seismicity in the region of the Nurek reservoir are analyzed. The analysis is based on the developed database for the earthquakes that occurred from 1955 to 1989 and is aimed at finding the regularities in the variations of the parameters of the transitional seismic regime caused by filling a reservoir. These parameters include the b-value—the slope of the graph of the Gutenberg–Richter frequency–magnitude relationship, the fractal dimension d of the set of the epicenters, and fracture cycle parameter q = αb ? d, where coefficient α determines the ratio between the magnitude and source size M = α log l + β. It is shown that during the filling of a reservoir, these parameters undergo statistically reliable variations: at the initial stages, the b-value increases, the fractal dimension of the set of epicenters decreases, and the fracture cycle parameter q grows and becomes positive in the middle of the time interval of reservoir filling. After a reservoir is filled, these parameters recover their background values. The aftershock sequences of the three strongest earthquakes—before, in the beginning, and in the middle of the reservoir filling period—are studied. It is confirmed that the Omori parameter p for the aftershock sequences during filling is smaller than for the earthquake before filling. Based on the dynamics of the studied parameters, it is conjectured that the relaxation time of the transitional seismic regime after the emergence of induced seismicity is about 10 years.     

Temporary seismological observations in the area of the 2012–2013 Tolbachik Fissure Eruption: Results

The seismicity that accompanied the Tolbachik Fissure Eruption was recorded by additional seismic stations that were installed in the southern Klyuchevskoi Volcanic Cluster area in January to October 2013. We used broadband (0.033–50 Hz) three-component digital Guralp CMG-6TD seismometers. This temporary network provided seismicity data at a lower energy level than can be done using the regional seismograph network of Kamchatka. The processing of the resulting digital records supplied data for compiling a catalog of over 700 M _L = 0–3.5 (K _S = 1.5–8.5) earthquakes, which is an order of magnitude greater than the number of events located by the regional network for the same period of time. The seismicity in the area of Ploskii Tolbachik Volcano was found to concentrate mostly in spatially isolated areas during the eruption. The main isolated clusters of earthquakes were identified both in the eruption area itself and along the periphery of Ploskii Tolbachik Volcano, in the area of the Zimina volcanic massif, and in the Tolud epicenter zone; the eruption zone was not dominant in the seismicity. The region of a shallow seismicity increase beneath Ploskii Tolbachik before the eruption was not found to exhibit any increased activity during the time the temporary seismograph network was operated, which means that a seismicity inversion took place at the beginning of the eruption. We discuss the question of what the earthquake-generating features are that we have identified.     

Parameter Variations in the Subsoil Radon Field at the Paratunka Station of the Petropavlovsk-Kamchatsky Geodynamic Test Site in 2011–2016

A network of stations for subsoil radon monitoring is in operation at the Petropavlovsk-Kamchatsky geodynamic testing area and is aimed at detection of strong earthquake precursors. At all stations, measurements are carried out using gas-discharge counters located at different depths within aeration zones of soft sediments. The volume activity of radon (VA Rn) is monitored at the most equipped station Paratunka (PRT) at three measurement sites located across the regional fault. The radon flux density (RFD) is measured from the surface. The article reviews responses in the dynamics of VA Rn and RFD from the surface at the PRT station prior to the Kamchatka earthquakes with magnitudes М _W > 5 that occurred over the period of 2011–2016. The revealed RFD seasonal cycle is likely related to the seasonal variations in air temperature. The postseismic effect caused by the strongest deep Okhotsk earthquake (May 24, 2013, М _W = 8.3) is detected in the RFD data. The behavior of VA Rn dynamics during time periods of the strong earthquakes is different. The results confirm the existing opinion on the formation of narrowly localized zones of Rn runoff to the atmosphere owing to both vertical and horizontal irregularities in the top layer of soil, which can react differently to changes in the geoenvironment stress–strain. On the basis of the real-time radon monitoring data, the authors have issued partially successful short-term prediction for several earthquakes. The results of this work confirm the opinion of many researchers that radon monitoring can be used in the short-term prediction of strong earthquakes.     

Features of Radiation and Propagation of Seismic Waves in the Northern Caucasus: Manifestations in Regional Coda

Based on the Anapa (ANN) seismic station records of ~40 earthquakes (M_W > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.     

Creepex of Kamchatka shallow earthquakes

Results of investigation of the lithosphere in the Kamchatka seismic focal zone from dynamic characteristics of earthquake records obtained at regional stations are presented. It is assumed that the specificity of the source zone can be estimated by the relation Cr = K _P ? bK _S ? c characterizing relative energies (energy classes, according to [Fedotov, 1972]) of short period transverse and longitudinal waves in the source. Azimuthal, spatial, and temporal variations in Cr and their relation to focal mechanisms are examined. Spatiotemporal variations in this parameter are shown to be caused by the influence of variations in the conditions in the source zone (its substance or process) on the radiation of P and S waves.     

On the possibility of using the tidal modulation of seismic waves for forecasting earthquakes

The earthquake forecasting technique based on the effect of the modulation of high frequency seismic noise (HFSN) by the Earth’s tides is discussed and the experience of its practical application in Kamchatka is described based on the extensive measurement data for 1996–2013. The empirical development of the lower magnitude threshold on the epicentral distance is refined. The reliability and validity of the precursor are assessed. The efficiency of the precursor is estimated in two ways. It is shown that the results of applying the prognostic procedure statistically significantly differ from a random guess. The presented estimates are based on the catalog containing 68 earthquakes of 1996–2013 in the predicted magnitude interval M ≥ 5.0 at the epicentral distances that are linked to the magnitude by functional dependence.     

Earthquake precursory signatures in electromagnetic radiation measurements in terms of day-to-day fractal spectral exponent variation: analysis of the eastern Aegean 13/04/2017–20/07/2017 seismic activity

The recent seismic activity in the eastern Aegean Sea, foregrounded by two major (M_L?=?6.1 and M_L?=?6.2) earthquakes occurring near Lesvos and Kos islands, respectively, is investigated in this work. Electromagnetic radiation measurements across different frequency bands and antenna orientations from a monitoring station in Agia Paraskevi, Lesvos Island, are analysed in order to reveal earthquake precursory signatures hidden in the electromagnetic data sequence. A straightforward, data-driven approach is employed in which day-to-day variations of the fractal characteristics of the measurements are adaptively monitored via a fractal spectral exponent similarity measure. The evolution of the fractal day-to-day variation in a 99-day period, which includes the two major earthquakes, shows a sustained, sudden increase lasting 1 to 4 days before every earthquake cluster of M_L?=?4.0 and above. Most importantly, this day-to-day variation subsides shortly after and remains relatively low during the absence of earthquakes, thus alleviating the emergence of false alarms.     

Propagation of internal gravity waves to the altitudes of the ionospheric <Emphasis Type="Italic">D</Emphasis> and dynamo regions in the seismic region (Kamchatka): Preliminary results

A spectral analysis of simultaneous diurnal variations in the E _z component of the quasi-static electric field in the near-Earth atmosphere, VLF radio noise, and the horizontal component of the geomagnetic field, observed at Kamchatka in September 1999, has been performed. These geophysical parameters are indirectly used to study wave processes in the near-Earth atmosphere and in the ionospheric D and dynamo regions within the band of periods of internal gravity waves (T = 0.5?3.5 h). The correlation method in the frequency region is used to analyze the interrelation between the wave processes in these atmospheric regions. The power cross-spectra of various pairs of geophysical parameters have been studied depending on meteorological, seismic, and geomagnetic activities. It is shown that the oscillations in the power spectra in the T ～ 1–1.5 h band of periods are caused by the sources of internal gravity waves in the near-Earth atmosphere and by the remote sources above the dynamo region of the ionosphere within the T ～ 1.5–3 h band of periods.