Vol 54, No 2 (2018)

Highly efficient approximate ray tracing techniques which can be used in seismic emission tomography and in other methods requiring a large number of raypaths are described. The techniques are applicable for the gradient and plane-layered velocity sections of the medium and for the models with a complicated geometry of contrasting boundaries. The empirical results obtained with the use of the discussed ray tracing technologies and seismic emission tomography results, as well as the results of numerical modeling, are presented.

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.

The study and radiocarbon dating of the low alluvial terraces of the Chon-Aksuu River, in the Northern Issyk-Kul region, which were broken by the Kebin (Kemin) earthquake of 1911 (Ms = 8.2, Io = 10 to 11), are carried out. The obtained radiocarbon dated ages refer to the second half of the Holocene. Since that time, at least eight strong earthquakes took place along this (Chon-Aksuu) segment of the Aksuu border fault. Three seismic events, including the earthquake of 1911 occurred in the second millennium A.D. This outburst of seismic energy was preceded by two millennia of seismic quiescence, which set in after another pulse of seismic activation. The latter lasted for 1.5 millennia and included five strong earthquakes. The recurrence period of seismic events during the activations is 300–600 years. Hence, the seismic regime along the Chon–Aksuu segment of the Aksuu border fault in the second half of the Holocene was a succession of two seismic activations, each with a duration of 1.0–1.5 ka, which were separated by a 2-ka interval of seismic quiescence. Therefore, the absolute datings of the river terraces of different ages which have been broken by a seismogenic rupture can serve as a reliable source of information about the age of the strong earthquakes that occurred along the seismogenic fault.

The paper continues the series of our works on recognizing the areas prone to the strongest, strong, and significant earthquakes with the use of the Formalized Clustering And Zoning (FCAZ) intellectual clustering system. We recognized the zones prone to the probable emergence of epicenters of the strongest (M ≥ 74/3) earthquakes on the Pacific Coast of Kamchatka. The FCAZ-zones are compared to the zones that were recognized in 1984 by the classical recognition method for Earthquake-Prone Areas (EPA) by transferring the criteria of high seismicity from the Andes mountain belt to the territory of Kamchatka. The FCAZ recognition was carried out with two-dimensional and three-dimensional objects of recognition.

Schemes and criteria are developed for using the measured and modeled geotemperatures for studying the thermal regime of the source rock formations, as well as the tectonic and sedimentary history of sedimentary basins, by the example of the oil fields of the Yamal Peninsula. The method of paleotemperature modeling based on the numerical solution of the heat conduction equation for a horizontally layered solid with a movable upper boundary is used. The mathematical model directly includes the climatic secular trend of the Earth’s surface temperature as the boundary condition and the paleotemperatures determined from the vitrinite reflectance as the measurement data. The method does not require a priori information about the nature and intensities of the heat flow from the Earth’s interior; the flow is determined by solving the inverse problem of geothermy with a parametric description of the of the sedimentation history and the history of the thermophysical properties of the sedimentary stratum. The rate of sedimentation is allowed to be zero and negative which provides the possibility to take into account the gaps in sedimentation and denudation. The formation, existence, and degradation of the permafrost stratum and ice cover are taken into account as dynamical lithological–stratigraphic complexes with anomalously high thermal conductivity. It is established that disregarding the paleoclimatic factors precludes an adequate reconstruction of thermal history of the source-rock deposits. Revealing and taking into account the Late Eocene regression provided the computationally optimal and richest thermal history of the source-rock Bazhenov Formation, which led to more correct volumetric–genetic estimates of the reserves. For estimating the hydrocarbon reserves in the land territories of the Arctic region of West Siberia by the volumetric–genetic technique, it is recommended to use the Arctic secular trend of temperatures and take into account the dynamics of the Neoplesitocene permafrost layers 300–600 m thick. Otherwise, the calculated hydrocarbon reserves could be underestimated by up to 40%.

Reliable use of strain data in geophysical studies requires their preliminary correction for ocean loading and various local distortions. These effects, in turn, can be estimated from the tidal records which are contributed by solid and oceanic loading. In this work, we estimate the oceanic tidal loading at two European strain stations (Baksan, Russia, and Gran Sasso, Italy) by analyzing the results obtained with the different Earth and ocean models. The influence of local distortions on the strain measurements at the two stations is estimated.

For theoretically studying the intensity of the influence exerted by the polarization of the rocks on the results of direct current (DC) well logging, a solution is suggested for the direct inner problem of the DC electric logging in the polarizable model of plane-layered medium containing a heterogeneity by the example of the three-layer model of the hosting medium. Initially, the solution is presented in the form of a traditional vector volume-integral equation of the second kind (IE2) for the electric current density vector. The vector IE2 is solved by the modified iteration–dissipation method. By the transformations, the initial IE2 is reduced to the equation with the contraction integral operator for an axisymmetric model of electrical well-logging of the three-layer polarizable medium intersected by an infinitely long circular cylinder. The latter simulates the borehole with a zone of penetration where the sought vector consists of the radial J_r and J_z axial (relative to the cylinder’s axis) components. The decomposition of the obtained vector IE2 into scalar components and the discretization in the coordinates r and z lead to a heterogeneous system of linear algebraic equations with a block matrix of the coefficients representing 2x2 matrices whose elements are the triple integrals of the mixed derivatives of the second-order Green’s function with respect to the parameters r, z, r', and z'. With the use of the analytical transformations and standard integrals, the integrals over the areas of the partition cells and azimuthal coordinate are reduced to single integrals (with respect to the variable t = cos ϕ on the interval [−1, 1]) calculated by the Gauss method for numerical integration. For estimating the effective coefficient of polarization of the complex medium, it is suggested to use the Siegel–Komarov formula.

The paradox of the long-term existence of ferromanganese concretions on the ocean floor is considered. It is known that nodules do not sink in the less dense medium of liquefied sediments and have been located on the bottom surface since their genesis over many thousands or even millions of years. At the same time, river and sea sludge easily sucks up other objects (of the same density), which in some way fall to the bottom. We propose an original mechanism explaining this paradox, which is based on a combination of three physicochemical phenomena.