Vol 52, No 1 (2016)

The P-wave travel time data from the earthquakes offshore and onshore around the Black Sea are used for the tomographic reconstruction of the three-dimensional (3D) velocity distribution in the lithosphere of the region. The preliminary refinement of the foci parameters (the coordinates and origin time) has reduced the random errors in the travel-time data. The earthquake data were supplemented by the previous deep seismic sounding (DSS) data on the profiles in Crimea and offshore off the Black Sea. The dataset included more than 4000 travel times overall. In order to eliminate the crustal effect, the travel times were reduced to a surface at a depth of 35 km corresponding to the mean Moho depth in the region. The improved crustal model was used for removing the contribution of the crust from the initial data. The new tomography method, which was recently developed by one of the authors and which relies on the assumption of smoothness of the lateral velocity variations, was applied for reconstructing the velocity structure of the upper mantle beneath the Black Sea up to a depth of 95 km. The lateral velocity variation maps at different depths and the vertical velocity distributions along the meridional and sublatitudinal cross sections across the Black Sea were constructed. High velocities were revealed in the subcrustal lithosphere, and the structural difference below two subbasins—the West Black Sea (WBS) and the East Black Sea (EBS) ones—was established. It shows that the high-velocity body below the WBS is located deeper than below the EBS and is distinguished by higher velocities. Based on these results, it is concluded that the lithosphere beneath the Black Sea has a continental origin.

The spatial models describing the physico-mechanical properties of the rocks within the seismically active segment of the Altai–Sayan region are constructed from the ground-based geophysical data. The comparative analysis of their probability density functions at the nodes of the spatial grid covering the studied domain and at the hypocenters of the previous earthquakes shows that the bulk and shear moduli are the most efficient markers of the weakened crustal zones, which are prone to earthquakes. The algorithm for constraining the potential earthquake sources based on the spatial distributions of the elastic moduli of the rocks is suggested. When applied to the Taellin segment of the Altai–Sayan region, this algorithm has localized four echelons of the potential source zones of earthquakes at different depths in the crust. The horizontal positions of two of these zones were delineated by the previous studies based on the ground observations of seismicity, whereas the other two zones have been identified by the suggested algorithm for the first time.

An advanced method for estimating the earthquake grouping parameters R_cr and T_cr is proposed in order to identify interrelated seismic events. The method pursues continuity with the previous algorithm suggested in (Mirzoev, 1980; 1988a; 1988b; 1992; Mirzoev and Azizova, 1983; 1984) but uses a more realistic spatial model of the background seismicity. All the calculations in the method can be fully formalized and a preliminary expert estimation of the parameters is not required. The method provides stable estimates of the critical radius R_cr and time T_cr of grouping. Group earthquakes make up 50 to 75% of their total number.

The Mikhnevo seismic station was founded in accordance with directive no. 1134 RS of the Council of Ministers of the Soviet Union of February 6, 1954. The station, installed south of Moscow, began its operations on monitoring nuclear tests in the United States and England in 1954. For dozens of years this station was the leading experimental base for elaborating new technical solutions and methods for monitoring nuclear explosions, equipped with modern seismological instruments. At present, the focus of activities has been moved from military applications to fundamental geophysical research. The station preserves its leading position in seismological observations due to the development of national high-performance digital instruments and creation of the small-aperture seismic array, the only one in the central part of European Russia, which is capable of recording weak seismic events with M_L ≥ 1.5 within a distance of 100 km.

The mantle convection model with phase transitions, non-Newtonian viscosity, and internal heat sources is calculated for two-dimensional (2D) Cartesian geometry. The temperature dependence of viscosity is described by the Arrhenius law with a viscosity step of 50 at the boundary between the upper and lower mantle. The viscosity in the model ranges within 4.5 orders of magnitude. The use of the non-Newtonian rheology enabled us to model the processes of softening in the zone of bending and subduction of the oceanic plates. The yield stress in the model is assumed to be 50 MPa. Based on the obtained model, the structure of the mantle flows and the spatial fields of the stresses σ_xz and σ_xx in the Earth’s mantle are studied. The model demonstrates a stepwise migration of the subduction zones and reveals the sharp changes in the stress fields depending on the stage of the slab detachment. In contrast to the previous model (Bobrov and Baranov, 2014), the self-consistent appearance of the rigid moving lithospheric plates on the surface is observed. Here, the intense flows in the upper mantle cause the drift and bending of the top segments of the slabs and the displacement of the plumes. It is established that when the upwelling plume intersects the boundary between the lower and upper mantle, it assumes a characteristic two-level structure: in the upper mantle, the ascending jet of the mantle material gets thinner, whereas its velocity increases. This effect is caused by the jump in the viscosity at the boundary and is enhanced by the effect of the endothermic phase boundary which impedes the penetration of the plume material from the lower mantle to the upper mantle. The values and distribution of the shear stresses σ_xz and superlithostatic horizontal stresses σ_xx are calculated. In the model area of the subducting slabs the stresses are 60–80 MPa, which is by about an order of magnitude higher than in the other mantle regions. The character of the stress fields in the transition region of the phase boundaries and viscosity step by the plumes and slabs is analyzed. It is established that the viscosity step and endothermic phase boundary at a depth of 660 km induce heterogeneities in the stress fields at the upper/lower mantle boundary. With the assumed model parameters, the exothermic phase transition at 410 km barely affects the stress fields. The slab regions manifest themselves in the stress fields much stronger than the plume regions. This numerically demonstrates that it is the slabs, not the plumes that are the main drivers of the convection. The plumes partly drive the convection and are partly passively involved into the convection stirred by the sinking slabs.