


Vol 53, No 3 (2017)
- Year: 2017
- Articles: 14
- URL: https://journals.rcsi.science/1069-3513/issue/view/13775
Article
Modeling the dynamics of the block structure and seismicity of the Caucasus
Abstract
Based on the morphostructural zoning scheme of the Caucasus, the block structure reflecting the real fault geometry and the block formation of the region is constructed. Several dozens of numerical experiments are conducted for simulating the dynamics of the block structure and the arising seismicity. The modeling relies on the following principles. It is assumed that the structure is composed of perfectly rigid blocks separated by infinitely thin fault planes. On the fault planes and on the blocks' bottoms, the blocks viscoelastically interact with each other and with the underlying medium. At each time instant, the translational displacements and rotations of the blocks are calculated from the condition of the quasi-static equilibrium of the entire block structure. The earthquakes occur in accordance with the dry friction model at the time instants when within a certain segment of the fault the stress-to-pressure ratio exceeds the given threshold. The modeling yields the synthetic catalog of the Caucasian earthquakes the spatial distribution of which reflects a set of characteristic features of the real seismicity. The similarity is observed in the magnitude–frequency diagrams of the synthetic and real seismicity. The comparison of the positions of the epicenters of the strong synthetic earthquakes with the results of recognizing the highly seismically active areas in the Caucasus demonstrates the presence of such epicenters in a few highly active areas where, according to the observations, strong earthquakes have not occurred to date.



Siberian traps: Hypotheses and seismology data
Abstract
Siberian traps are the result of huge basalt eruptions which took place about 250 Ma ago over a vast territory of Siberia. The genesis of Siberian traps is attributed to a mantle plume with a center in the region of Iceland or beneath the central Urals in terms of their present coordinates. The eruption mechanism is associated with delamination—replacement of the mantle lithosphere by the deep magma material. The receiver function analysis of the records from the Norilsk seismic station (NRIL) allows comparing these hypotheses with the factual data on the depth structure of the region of Siberian traps. The S-wave velocity section place the seismic lithosphere/asthenosphere boundary (LAB) at a depth of 155–190 km, commensurate with the data for the other cratons. The mantle lithosphere has a high S-wave velocity characteristic of cratons (4.6–4.8 km/s instead of the typical value 4.5 km/s). The seismic boundary, which is located at a depth around 410 km beneath the continents is depressed by ~10 km in the region of the NRIL station. The phase diagram of olivine/wadsleyite transformation accounts for this depression by a 50–100°С increase in temperature. At the depths of 350–400 km, the S-wave velocity drops due to partial melting. A new reduction in the S-wave velocities is observed at a depth of 460 km. The similar anomalies (deepening of the 410-km seismic boundary and low shear wave velocity at depths of 350–400 and 460–500 km, respectively) were previously revealed in the other regions of the Meso-Cenozoic volcanism. In the case of a differently directed drift of the Siberian lithosphere and underlying mantle at depths down to 500 km, these anomalies are barely accountable. In particular, if the mantle at a depth ranging from 200 to 500 km is fixed, the anomalies should be observed at the original locations where they emerged 250 Ma ago, i.e. thousands of km from the Siberian traps. Our seismic data suggest that despite the low viscosity of the asthenosphere, the mantle drift at depths ranging from 200 to 500 km is correlated with the drift of the Siberian lithospheric plate. Furthermore, the position of the mantle plume beneath the Urals is easier to reconcile with the seismic data than its position beneath Iceland because of the Siberian traps being less remote from the Urals.



S-wave velocity model for several regions of the Kamchatka Peninsula from the cross correlations of ambient seismic noise
Abstract
The data from the seismic networks of the Kamchatka Branch of the Geophysical Survey of the Russian Academy of Sciences are used for calculating the cross correlations of seismic noise for the stationary digital stations over 2013 and for radio telemetric stations (RTS) in the region of the Klyuchevskoy volcano over the period from January 1, 2009 to May 31, 2013. Four hundred and two correlations overall are calculated. The fundamental-mode group velocities of the Rayleigh waves are calculated in the periods ranging from 5 to 50 s. The calculations for the region of the Klyuchevskaya group of volcanoes are based on the RTS data and cover the periods from 2 to 8 s. The two-dimensional (2D) maps of group velocity distributions in different periods are constructed with the use of the algorithm of surface wave tomography (Barmin, 2001). The velocity sections for the selected Kamchatka regions are reconstructed by the dispersion curve inversion technique (Mordret, 2014). For each region, the structure of the Earth’s crust and upper mantle down to a depth of 50 km was obtained.



Formalized clustering and significant earthquake-prone areas in the Crimean Peninsula and Northwest Caucasus
Abstract
With the use of the modified version of the original algorithmic Formalized Clustering and Zoning (FCAZ) system, the areas prone to the probable emergence of the epicenters of significant earthquakes are recognized in the joint region of the Crimea and western part of the Northern Caucasus. The selection of this region is justified by the tectonic structure and the presence of the active junction zone of the meganticlinoria. The reliability of the obtained recognition is substantiated by the comparative analysis of the actual and random FCAZ-recognition. For the first time, the problem of recognizing the locations of the probable emergence of the earthquakes' epicenters is solved for two different magnitude thresholds. This allows us to interpret the areas prone to the probable emergence of the epicenters of significant earthquakes as fuzzy sets.



Forecasting aftershock activity: 2. Estimating the area prone to strong aftershocks
Abstract
The technique for forecasting the spatial domain where fairly intense aftershocks should be expected after a strong earthquake is considered. The paper presents the task of estimating the area prone to the strong future aftershocks using the data for the first 12 h after the main shock. The existing aftershock identification techniques are inapplicable to this task because they either analyze the distributions of the epicenters of the aftershock process that has been already completed or only consider the parameters of the main shock and only provide rough estimates. Using the developed criteria of estimating the quality of the prediction, we quantitatively compared quite a few different candidates. The latter included the main known techniques and their modifications suggested by us. In these modifications, we took into account the results of the recent studies on the dynamics of the aftershock process. This enabled us to select the optimal procedure which demonstrated the best results of the quantitative tests for more than 120 aftershock sequences with the magnitudes starting from 6.5 all over the world. This procedure can be used in the seismological monitoring centers for forecasting the area prone to the aftershock activity after a strong earthquake based on the data of operative processing.



A study of the energy released by great (M ≥ 7) deep focus seismic events with allowance for the Mw 8.3 earthquake of May 24, 2013 in the Sea of Okhotsk, Russia
Abstract
The distribution of the hypocenters of great seismic events with M ≥ 7.0 and, consequently, the seismic energy released in their sources is asymmetric along the Earth radius. According to our estimates, 90% of the energy is released relatively close to the Earth’s surface, at the average depth of 50 km. The bulk of the remaining 10% is associated with seismic events that take place at large depths, on average 630–640 km above the boundary between the transition zone and lower mantle. These very deep earthquakes (depth ≥350 km) significantly differ from the shallow events. Their sources, in contrast to the shallow focus events, are located inside the plate. The examination of seven seismic zones described in the present study, except for the Honsu-Kamchatka zone which accommodates both the shallow and deep M ≥ 7.0 earthquakes, shows that the linear distribution of the hypocenters of deep earthquakes is considerably shorter than that of the shallow earthquakes, which, in turn, determines the length of the seismic zones at depth. In the zones of the Solomon Islands and the Izu–Bonin–Mariana arc, there are no seismic events with M ≥ 7.0 deeper than 450 km. In the zones of Indonesia, Philippines, Tonga–Kermadec, and Chile–Peru, the mentioned length’s shortening at the top of the lower mantle (660-km discontinuity) relative to the length of the zone observed close to the surface is unequivocal. The relationship between the lithospheric plates is supported by the spatial distribution of the hypocenters. The position of the foci of very deep (≥500 km) earthquakes indicates where the descending lithospheric plates conflict with the upper boundary of the lower mantle, and where in some cases they cross it. This passage generates the compression and elongation inside the slab. A comparison of the time distribution of the shallow and deep seismic events suggests the absence of direct relationship between these two different earthquake activities. For studying the fairly uncommon deep earthquakes, important additional information was provided by the largest of the deep earthquakes, the May 24, 2013 M 8.3 event beneath the Sea of Okhotsk, in an area where significant deep earthquakes have already occurred. Based on our studies of the records provided by the Russian and Hungarian national seismological networks, we concluded that this seismic event was preceded by an earthquake swarm, which consisted of 58 M ≥ 5 events between May 15 and 24, 2013 in the higher part of the sinking slab east of Kamchatka within the segment of increased historical seismicity. Most probably, the interaction of two distinct active source zones took place. The aftershock activity beneath the Sea of Okhotsk was moderate: thirteen events with magnitudes above М ≥ 4.0 were observed by June 27, 2013. Nevertheless, the unusually small number of aftershocks determined a fault area (2.64 × 104 km2), generally similar to that in the case of the assumed shallow M 8.3 event.



Paleomagnetism of early paleozoic geological complexes of Mongolia
Abstract
Most of the studied Early Phanerozoic rocks of West Mongolia have undergone repeated remagnetization. Secondary magnetization components with normal and reversed polarity are isolated. The magnetization components with normal polarity are associated with the Mesozoic remagnetization of the rocks. The components with reversed polarity were probably formed during the Carboniferous–Permian superchron of reversed polarity. The analysis of the distribution of the reversed-polarity magnetization component in the structure of Mongolia permits some zonation. Within Mongolia, the regions with insignificant post- Permian deformations and complicated post-Permian deformations are identified; also the area of rotations of large geological blocks about the horizontal axis (Khan-Khukhei Ridge) is distinguished. It is hypothesized that in the Ordovician rocks of West Mongolia, the magnetization component that is close to primary was identified. If this is the case, the paleolatitude calculated from this magnetization direction corresponds to the interval 14°–17°–20° (minimal–mean–maximal) of probably northern latitude



Deep electrical conductivity features in the transition zone from the Pacific to Eurasia
Abstract
Understanding the processes that occur in the transition from the Pacific Ocean to Eurasia is key to constructing the tectonic models of the Earth’s shells and the convection models of the upper mantle. The electromagnetic methods permit estimating the temperature and fluid content (and/or carbon (graphite) content) in the Earth’s interior. These estimates are independent of the traditionally used estimates based on seismic methods because the dependence of electrical conductivity on the physical properties of the rock is based on different principles than the behavior of the elastic waves. The region is characterized by a complicated geological structure with intense three-dimensional (3D) surface heterogeneities, which significantly aggravate the retrieval of the information about the deep horizons in the structure of the Earth’s mantle from the observed electromagnetic (EM) fields. The detailed analysis of the nature of the deep electrical conductivity and structural features of the transition from the Pacific to Eurasia included numerical modeling of the typical two- and three-dimensional models has been carried out. Based on this analysis, the approaches that increase the reliability of the interpretation of the results of the EM studies are suggested.



On the origin of background fluctuations in electric field measurements on the seafloor
Abstract
The origin of background fluctuations in the natural electric field observed in the measurements on the seafloor of shallow seas and oceans, as well as in the lakes and estuaries of large rivers is considered. It is shown that the presence of these fluctuations, which have a quasi-harmonic behavior, is due to the surface waves and can be accounted for by the action of the (1) magnetohydrodynamic and/or (2) electrokinetic effects, with the predominance of the latter in the freshwater basins. The excitation models of the background fluctuations in the electric field by both the long surface waves and interacting short waves are analyzed. The estimates of the electric field of background fluctuations are obtained.



On the solution stability in the problem of determining the time variations of the tidal responses of a medium in the vicinity of the sources of severe earthquakes
Abstract
In (Molodenskii M.S. et al., 2016), the data from horizontal pendulums recording the tilts in the closest vicinity of the Great Tohoku earthquake of March 11, 2011 in Japan were analyzed. A significantly improved method for statistical analysis of the observational data enabled the authors to reveal a slow growth in tidal tilts during a period of six years before the earthquake, which was superseded by an instantaneous drop in the amplitudes at the time of the earthquake. After this, during the subsequent four years, the tidal amplitudes have remained at a significantly lower level than their average values before the earthquake. These changes in tidal amplitudes testify to the nonlinear character of the tidal response of the medium in the presence of large tectonic stresses: as is well known, the linear relationship between stresses and strains in a real medium is only the case for stresses that are far below the yield stress. When the stresses approach the failure limit, two counteracting effects come into play: (1) the shear moduli in some areas decrease as a result of the avalanche growth of the crack formation processes, and (2) the moduli increase due to the compression in the other areas. Irrespective of which particular effect of these two is predominant, in either case the linearity of the relationship between the stresses and strains should be violated. This violation cannot but affect the amplitudes of the tidal tilts and strains characterizing this relationship in the presence of fairly low additional tidal stresses (i.e., the derivative of the off-diagonal stress tensor components with respect to the same components of the strain tensor). Since there is presently a sufficiently dense network of the horizontal pendulums recording the tilts (the global IRIS network and the particularly dense F-NET network in Japan), monitoring the changes in the amplitudes of tidal tilts can be considered as a key instrument for capturing the signs of the approach of tectonic stresses to their critical values. The increase in tidal amplitudes before the Tohoku earthquake and their drop at the moment of the earthquake, which were revealed by us, as well as the constancy of the amplitudes during four years after the event, unambiguously indicate that the accumulation of tectonic stresses caused the growth in tidal amplitudes, whereas the stress release by the earthquake caused their diminution. This does not however mean that stress accumulation is accompanied by a decrease in the elastic moduli and that the release of stresses is accompanied by the growth of elastic moduli all over the source area. As was shown in (Molodenskii M.S. et al., 2012), even in the simplest model of spatially homogeneous variations of elastic modules, the variations in tidal tilts are an odd function of the distance from the epicenter. Therefore, irrespective of whether the elastic moduli decrease or increase, the amplitudes of tidal tilts should decrease in some areas and increase in other areas. Hence, the very fact of the growth of tidal tilt amplitudes with time cannot be considered as a sign of the growth of tectonic stresses. To be positive about the latter, one should make sure that the consistent (unidirectional) changes have been observed during a sufficiently long time interval and that their magnitudes were significantly larger than the measurement errors. Hence, it is important to reliably estimate the errors of the observational data.



Tectonic stress accumulation in the region of the Kuril–Kamchatka Island Arc system
Abstract
The time variations in the tidal response of the medium in 2011–2015 according to the measurements at the global seismographic network (GSN) in Kamchatka are considered. Based on the data from the horizontal pendulums recording the eastward tilts at the station, it is established that there was a linear growth in the tidal tilt amplitudes up to May 24, 2013 due to the changes in the elastic moduli caused by tectonic stress accumulation. The growth phase was followed by the decline in the tidal tilt amplitudes induced by the release of stresses after the Sea-of-Okhotsk earthquake.



Paleomagnetism of carboniferous formations of the Russian platform: Reinvestigation of old collections
Abstract
The collections of Carboniferous rocks from sections of the Russian Platform (Gzhelian, Moscovian, Bashkirian, and Visean stages) are studied. The new mean paleomagnetic poles are obtained from the Gzhelian, Moscovian, and Visean layers of the Carboniferous of the Russian Platform. In the redbed Gzhelian and Moscovian rocks, the natural remanent magnetization (NRM) components with the inclination shallowing are revealed, which is due to the presence of the large hematite particles or particle aggregates associated with the interaction between the magnetic and clay particles. Based on the obtained determinations and the results contained in the World paleomagnetic database, the trajectory of the apparent polar wander path (APWP) for the East European Platform is constructed in the interval from the Devonian to Early Permian. The Carboniferous kinematics of the East European Platform is estimated.



The T-matrix approach for the mathematical modeling of the effective elastic properties of hydrocarbon reservoirs
Abstract
Various methods of the effective medium theory (EMT)—the T-matrix approach (optical potential approximation and coherent potential approximation), Mori–Tanaka method, generalized singular approximation, etc.—for modeling the effective elastic properties of hydrocarbon reservoir rocks are considered. The relationship of the different approaches (perturbation theory, self-consistent methods, and variational principles) is demonstrated. The classification of the methods by the degree of complexity of solving the inverse problem is suggested. An example of the theoretical modeling of an oil-, gas-, and water-saturated rock with the oriented axis of the fractures is presented. Conclusions about the applicability of different methods for rock modeling are made. In this study we review and compare different methods of the EMT and present the guidelines for determining the effective properties of the medium in rock physics modeling.



Discussion
On the role of magnetostrophic waves in geodynamo
Abstract
It is shown that magnetostrophic waves which are generated in the equatorial plane of the Earth’s core due to the instability of the equatorial jet and which propagate almost transversely to the rotational axis off the tangent cylinder, have a negative helicity in the northern hemisphere and positive helicity in the southern hemisphere. When the wave trains propagate through the regions with a constant azimuthal magnetic field caused by the Ω-effect, this helicity distribution induces an electromotive force (emf) (due to the α-effect), which may lead to the maintenance of the initial dipole field by the scenario of the α-Ω dynamo.


