Vol 54, No 4 (2018)

Correcting the effects of the sphericity of the Earth in the results of the interpretation of gravimetric data is a topical issue in modern gravimetry. Estimating the error of the gravity field calculations due to the replacement of the spherical Earth model by the plane model is an important part of this problem. In this paper, a method is proposed for transforming the plane density models into spherical ones and vice versa. Algorithms for calculating the vertical component of gravity field for both model types are presented. For two extensive plane models of the Earth’s density, their transformation into spherical models is carried out and the resulting gravity fields are compared. The relative root mean square residuals between the fields calculated with this replacement are at most 5%.

The results of the laboratory and numerical experiments in circular rotating trays with thin layers of a conductive fluid under the MHD generation of small-scale velocity fields are presented. The configurations of constant magnets for MHD generation were determined based on the numerical calculations with shallow water equations. Both the laboratory and numerical experiments with rotating trays demonstrate the emergence of nonaxisymmetric structures and large-scale near-circular vortices caused by the energy transfer from the system of the externally generated small-scale vortices to the large-scale velocity fields under the action of the Coriolis force. The near-circular vortex has areas with differential rotation when the angular velocity of rotation decreases with the radius. The single large-scale vortices and wide jet flows arise in the regimes of subrotation and superrotation relative to the external rotation depending on its angular velocity. The emergence of the flow structures with the azimuthal wave number m = 2 is demonstrated, and their probable relation to the anomalies of the geomagnetic field observed on the Earth’s surface is considered.

New representations of the elements of the fields of attraction (potential and its first derivatives) are presented for the important approximating models such as polyhedron and polygonal plates with the density varying by the linear law. It is shown that these elements are determined through the elements of the fields from the models with the known analytical representations (polyhedron, polygonal plate, and material segment with constant densities) and the additional integrals for which the explicit analytical expressions exist.

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.

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.

The airborne gravity measurements both in Russia and abroad widely use the GT-2A aerogravimetric complex. This complex has a number of benefits but suffers from a serious limitation: in high latitudes, it is incapable of making measurements on westward courses. In this paper, we consider the causes of this drawback and propose theoretical solutions to fix it. Based on the obtained extensive empirical data, we select the most suitable solution. This enables us to suggest a modification of the GT-2A gravimetric complex for conducting airborne gravity surveys in the regions of the polar caps of the Earth.