Vol 25, No 3 (2025)
Articles
METHODOLOGICAL ASPECTS OF PREDICTIVE MINERAGENIC STUDIES USING EARTH REMOTE SENSING DATA
Abstract
The article considers methodological aspects of allocation and substantiation of exploration areas for scarce types of ore minerals taking into account the concept of of mineral systems and using Earth remote sensing data with the application of geoinformation and neural network technologies using the example of the Argun metallogenic zone in South-Eastern Transbaikalia. Of the entire range of areas of fundamental and exploratory scientific research, the main attention within the framework of predictive and mineragenic studies is paid to solving the following problems: 1) allocation of lineaments (fault zones) based on processing of digital elevation models; 2) determination of hydraulically active fault structures for the period of ore formation based on tectonophysical reconstructions; 3) analysis of multispectral characteristics of pre-ore, ore-accompanying and post-ore metasomatites based on statistical processing of Landsat-8 satellite data; 4) assessment of fluid-dynamic settings of deposit formation based on data on the composition, properties and genesis of mineral-forming fluids. 5) creation of weight of evidence models based on statistical algorithms for processing data on the dynamics of ore-genetic processes. The feasibility of using such an approach for setting up predictive mineragenic studies on scarce types of strategic mineral raw materials in areas with complex climatic and landscape conditions is shown.
Russian Journal of Earth Sciences. 2025;25(3):ES3001
ES3001
COMBINATORIC-VARIATIONAL APPROACH TO SOLVING INVERSE PROBLEMS IN GEOPHYSICS WITHIN THE FRAMEWORK OF THE DISCRETE POTENTIAL THEORY. COMPLEX GEOLOGICAL ENVIRONMENT
Abstract
Variational-combinatoric approach is used when solving linear and nonlinear inverse problems of geophysics within the framework of discrete potential. The acquired data about physical fields and their sources are of discrete character. Combinatoric methods of discrete mathematics allow us to organize some selective choice of the information on the physical field elements. Analytical approximations are built only by means of verified data which can be recognized as useful signal, not as an arbitrary or systematic noise.
Russian Journal of Earth Sciences. 2025;25(3):ES3002
ES3002
USE OF VARIABLE MATRICES OF LAND COVER CHANGES FOR THE ASSESSMENT OF LAND DYNAMICS BASED ON THE INDICATORS OF LAND DEGRADATION NEUTRALITY
Abstract
The article presents the results of land degradation assessment for the Ryazan and Tula regions based on the concept of land degradation neutrality (LDN). For the first time the approaches were demonstrated for developing a methodology for assessing LDN sub-indicators for humid areas (using the moving average method and modified matrices of changes in land cover type). For forest lands, a matrix modified to take into account the dynamics of woody vegetation, including mutual transitions of coniferous, broad-leaved and small-leaved forests, is proposed. New data on land degradation trends for the studied territories were obtained. The possibility of using the LDN methodology to update official statistics, detail land productivity trends and land use changes is confirmed. Despite the decrease in the degradation rate, in the considered territories, degraded and moderately degraded lands in total still significantly prevail over improved ones, which requires active intervention in regulating land use conditions and the use of new methods. At the same time, the situation in the regions is somewhat different: if in the Tula region degraded lands are mainly confined to the forest-steppe subzone, then in the Ryazan region, foci of degraded lands are equally often found in forest areas. A high dependence of land degradation trends on even short-term (within 5–10 years) climate changes has been revealed. For the Ryazan and Tula regions, such reactions (improvement or deterioration) in individual areas and in the region as a whole can be very contrasting, which requires special consideration when developing measures to adapt to climate change.
Russian Journal of Earth Sciences. 2025;25(3):ES3003
ES3003
FRACTAL CHARACTERISTICS OF THE AURORAL OVAL STRUCTURE ACCORDING TO THE ALL-SKY CAMERA DATA IN APATITY FOR 2013–2020
Abstract
The spatial structure of polar auroras is described by the fractal dimension of glow fluctuations and its anisotropy from direction. The fractal dimension is estimated from the slope in the logarithmic axes of the spectrum in the range of 1.5–50 km, obtained through discrete wavelet transformation of the intensity fluctuation of the glow using Daubechies wavelets of order 5. The variability of the structures is characterized by the slope of the anisotropy variation spectrum over time. The statistics of these characteristics are presented according to the data of the ground-based all-sky camera of the Polar Geophysical Institute in Apatity for 2013–2020 and referenced to the position inside the auroral oval and the values of the geomagnetic field at Lovozero observatory. An algorithm for modeling the structure of polar auroras based on these characteristics is discussed.
Russian Journal of Earth Sciences. 2025;25(3):ES3004
ES3004
SM-ND SYSTEM OF RARE-METAL PEGMATITES OF THE WORLD-CLASS KOLMOZERO LITHIUM DEPOSIT AND SHONGUI BERYLLIUM DEPOSIT, NW RUSSIA: GEOCHEMICAL CAUSES OF DISTURBANCE AND ND MOBILITY
Abstract
The Sm-Nd geochronological study was performed to investigate rare-metal pegmatites from the two unique deposits, i.e., the Kolmozero lithium deposit and Shongui deposit with beryllium mineralization (Kola Peninsula, Russia). For Kolmozero lithium deposit was obtained the Sm-Nd isochrone, corresponding to an age of 1705 ± 60 Ma with high εNd(T ) = +9.1. For Shongui beryllium deposit was obtained the Sm-Nd age, corresponding to an age of 1747 ± 33 Ma with high εNd(T ) = +9.7. This age values are close to the ages of metamorphism obtained earlier by Rb-Sr, K-Ca, and K-Ar methods. Pegmatites are characterized by a wide range of εNd(T ) values from +2 to +16 and 147Sm/144Nd ratios up to 0.3. Possible reasons for disturbance of the Sm-Nd isotope system of pegmatites are analyzed, including multicomponent mixing, fluid influence and metamorphic overprinting. The highly radiogenic signatures of rare metal pegmatites of the Shongui and Kolmozero deposits were found to appear by fractionation of Nd and Sm and their different redistribution with the change of the Sm/Nd ratio. High εNd(T ) values and changes in Sm/Nd ratios indicate the role of REE (rare-earth element) fractionation, while narrow εNd(T ) ranges suggest interaction with fluids during pegmatite formation. These findings emphasize the need for further research into the composition of fluids and their influence on isotope systems.
Russian Journal of Earth Sciences. 2025;25(3):ES3005
ES3005
THE PARADIGM OF A DIGITAL IMAGE OF A GEOTECHNICAL OBJECT AS THE BASIS OF AN AUTOMATED WORKSTATION IN THE TASK OF NON-DESTRUCTIVE TESTING (USING THE EXAMPLE OF A MAIN PIPELINE)
Abstract
The relevance of the work is determined by the automation of monitoring processes of natural and technical systems, not so much in terms of their measurement component, but in terms of interpreting monitoring results in terms of the properties of the desired object. In this work, these are submerged pipeline systems – main oil and gas pipelines, as a special case of geotechnical objects. The goal of developing an automated workstation (AWS) for a specific object is the ultimate parameterization of the interpretation procedure and the unification of the final graphical constructions within the framework of a single digital image of main pipelines. Taking into account the deterministic type of geotechnical object and the objectives of the survey, research methods include, in addition to non-destructive magnetometric testing as a method of obtaining primary data: methods of qualitative and quantitative interpretation of multi-channel magnetometry materials; methods of systematics and graphical display of primary data and the results of their interpretation. Within the framework of approaches generally accepted in the oil and gas industry, the digital image of damaged pipeline sections is based on ideas about the thinning of pipe walls and the stress states associated with thinning zones. The stressed states of the ferromagnetic walls of immersed pipelines form specific responses in an external magnetic field, while the actual thinning of the pipeline walls is detected by contact methods. The results of the surveys are reduced to a system parametrization of magnetometric non-destructive testing data within the framework of the developed automated workplace structure. The content of the conclusions is in the development of the principle for referenceless recognition of risky sections of the pipeline as a fully automated component of diagnostician’s workplace.
Russian Journal of Earth Sciences. 2025;25(3):ES3006
ES3006
NUMERICAL SIMULATION OF RADAR SIGNAL REFLECTED BY SEA SURFACE WITH DIFFERENT SEA ICE CONCENTRATION
Abstract
Currently, methods of radar remote sensing at small incidence angles (from the vertical to 15∘) are actively developed. An important application of these methods is the determination of the presence and sea ice concentration. This paper presents an approach to numerical simulation of an experiment in which a reflecting surface with different sea ice concentrations is modeled and then the characteristics of the radar signal reflected by this surface for a given measurement geometry are modeled. Without loss of generality, in this paper we will consider a specific geometry of the DPR (Dual-frequency Precipitation Radar) radar on the GPM (Global Precipitation Measurement) mission satellite and only the Ku-band of this radar. The signal reflected by sea waves will be calculated within the Kirchhoff approximation. Since there is no generally accepted model for the signal scattered by the sea ice at small incidence angles, an empirical formula obtained from DPR data will be used as a model. The paper discusses a method for determining sea ice concentration using radar sensing data at low incidence angles.
Russian Journal of Earth Sciences. 2025;25(3):ES3007
ES3007
A NEW METHOD FOR RETRIEVING REMOTE SENSING REFLECTANCE FROM FIRST-LEVEL OLCI SATELLITE DATA
Abstract
The paper proposes an alternative method of atmospheric correction using the OLCI satellite data for the Black Sea as an example. Currently, for remote sensing problems, the standard Gordon and Wang atmospheric correction algorithm is used in most cases (GW94). Unfortunately, its operation is often accompanied by the appearance of negative values of the spectral radiance coefficient of the sea (remote sensing reflectance) rs( ) in the shortwave region, which means a sufficient number of physically incorrect values and subsequent incorrect calculation of the concentration of chlorophyll-a and yellow matter. In this paper, a simple algorithm is proposed, built exclusively on analytical formulas, where two procedures of interpolation and extrapolation are conceptually implemented simultaneously, extrapolation - via two channels, interpolation based on the constancy of the color index ratio (CI = rs(412)/ rs(443) = 0.8). Using individual examples of OLCI scanner data, the performance GW94 of the new algorithm was tested for different states of the atmosphere and sea surface by comparing the results with in-kind measurements of the AERONET-OC platforms, with Level-2 data and with the operation of the regional method of additional correction. The new algorithm was tested under the following conditions: clear atmosphere (presence of background aerosol), presence of dust aerosol, cloud boundaries, presence of sun glare, coccolithophore bloom. When analyzing a number of Sentinel 3A/3B satellite images, it was found that the new simple algorithm was, on average, better than the standard one, which means that there is a prospect for its improvement. The advantage of this approach is its universality and the possibility of its implementation for other water areas, if there are patterns in the variability of the "blue" color index.
Russian Journal of Earth Sciences. 2025;25(3):ES3008
ES3008
ON THE PROBLEM OF RADAR MONITORING OF ICE COVER FORMATION AND DESTRUCTION ON INLAND WATERS: FIRST ASSESSMENTS
Abstract
Remote sensing methods are widely used to monitor ongoing climate change, and the area of sea ice in the Arctic and Antarctic is used as one of the criteria. The duration of ice cover on inland waters can be used to assess processes within continents. The small size of internal waters does not allow the use of traditional methods that have proven themselves well in marine conditions. This study considers the possibility of using dual-frequency precipitation radar data to detect ice formation and destruction on small inland waters. Due to the features of backscatter at small incidence angles (< 18∘), inland waters with sizes smaller than the radar resolution (5 km) are “visible” in radar images. Using the Volga River as an example, it is shown that water-ice and ice-water transitions can be detected when analyzing radar images, and thus, possible to estimate the duration of ice cover on inland waters.
Russian Journal of Earth Sciences. 2025;25(3):ES3009
ES3009
THE FUNDAMENTALS OF A TWO-STAGE APPROACH TO SYSTEMATIC EARTHQUAKE PREDICTION
Abstract
A systematic earthquake prediction is performed regularly at fixed intervals within a preselected seismically homogeneous zone. The result of each prediction iteration is a map highlighting the alarm zones, where the epicenters of target earthquakes are expected. The proposed methodology introduces the following innovations: 1 – A prediction is considered successful if all epicenters of the target earthquakes during the forecast interval fall within the alarm zone. 2 – The methodology optimizes both the probability of successfully detecting earthquake epicenters across a series of forecasts and the success rate of predictions in each individual iteration. 3 – The methodology enables the estimation of the probability of success for the next forecast interval. Examples of the method's application are demonstrated for predicting earthquakes in Kamchatka, California, and the island region of Japan.
Russian Journal of Earth Sciences. 2025;25(3):ES3010
ES3010
DETERMINATION OF 3D GRAVITY SOURCE AND ITS DEPTH BENEATH CAMEROON VOLCANIC LINE (CVL) USING DEXP TRANSFORMATION
Abstract
The Cameroon Volcanic Line (CVL), which is around 1600 km long, the Adamawa Plateau, the northern sedimentary basins, the Central African Shear Zone, and the northern boundary between the Pan-African Mobile Belt and Congo Craton are the primary geological features of Cameroon. A good number of authors have attempted to comprehend the geology and gravity field along the CVL by using gravitational data from the EGM2008 model to analyze the gravity effects in areas around Cameroon and the CVL with a focus on its structures and subsurface characteristics. Despite the fact that many authors have written on the subject matter, more emphasis has to be placed on the determination of the gravity source and depth beneath the CVL. Experimental gravity field model XGM2019e_2159 and DTU21 data were used in this research to estimate the depth of the gravity source. Both DEXP (Depth from Extreme Points) and spectral analysis were carried out to complement the results and accuracy of the techniques. The main focus of this research is to investigate the gravity source depth of CVL using DEXP as the main approach to illustrate its application in solving geophysical and geologic problems and reveal details of volcanic structures beneath the CVL. In this work, we describe the steps taken to calculate the anomalous gravity field and regional and residual gravitational effects. We further performed application of the DEXP transformation of 3D gravity field distribution to produce a 3D model for the depth of gravity sources.
Russian Journal of Earth Sciences. 2025;25(3):ES3011
ES3011
ANALYSIS OF WIND WAVES IN ANIVA BAY BASED ON MODEL RESULTS
Abstract
In the Bay of Aniva (Sea of Okhotsk) there is one of the carbon polygons where greenhouse gas balance studies are conducted. Wind waves play an important role in energy, material and gas exchange in the ocean-atmosphere and ocean-land systems, ensuring the transfer of matter and gases between all biotic and abiotic components. This paper presents an analysis of the wave height, period and wavelength based on the database of wave modeling from 1979 to 2019. The maps of the main parameters and their statistical distributions are analyzed. Wind and wave roses are presented for several points in different parts of the bay. The average long-term of significant wave height is up to ∼ 0.5 m in the northern part of the bay and up to ∼ 1 m in the southern part. The maximum significant wave height reaches 8.5 m in the southern part of the bay and about 4–5 m in the north. The most dangerous areas for waves are the southwest, south and southeast.
Russian Journal of Earth Sciences. 2025;25(3):ES3012
ES3012
INTERPRETATION OF A DIGITAL MODEL OF THE ANOMALOUS MAGNETIC FIELD OF THE CRIMEAN REGION AND ADJACENT WATERS, BASED ON RETROSPECTIVE DATA
Abstract
The article presents the results of work on the interpretation of a digital model of the anomalous magnetic field (AMF) of the Crimean region and adjacent waters, based on retrospective data. For interpretation, a priori data in various fields of geology and geophysics were studied in the public domain. The main interpretation tool was the comparison of the calculated high-frequency and low-frequency magnetic field transformants with a priori geological and geophysical information. Based on the results, a scheme of the tectonic structures of the Crimean region reflected in its magnetic field was constructed.
Russian Journal of Earth Sciences. 2025;25(3):ES3013
ES3013
2D TOMOGRAPHIC MAPS OF RAYLEIGH AND LOVE WAVES FOR ARMENIA AND ITS SURROUNDING AREAS
Abstract
This study aims to provide 2D tomography images of surface wave group velocity in Armenia and its neighboring regions within the Eurasian-Arabic plates, aiming to enhance the understanding of the shear velocity structure in the area. In this context, ∼ 516 earthquakes (M ≥ 3.5) recorded by 20 stations between 1999–2018 were analyzed, and the surface wave group velocity dispersion curves for each record (source-station path) were estimated. Subsequently, taking advantage of a 2D-linear inversion procedure, 2D tomography maps for periods of 5–80 s (Depth ≈ 180 km) were computed with a grid spacing of 0.2◦ × 0.5◦. The results of this research indicate that short-periods (5 ≤ T ≤ 25 s) are more influenced by shallow, ever-evolving deformations within various geological units, such as sedimentary basins. A minor low-velocity anomaly (Depth ≈ 27 km) is observed on the northwest slope of the Aragats volcano, which contrasts with the findings of earlier studies. For Rayleigh and Love waves, the North Armenia Block is covered by a high-velocity anomaly. The results for medium periods reveal the presence of very-high-velocity anomalies in some geological units (e.g., Lesser Caucasus), aligning with the ongoing subduction processes. In contrast, very lowvelocity anomalies reflect the uppermost mantle information, revealing an extremely thin lithosphere accompanied with a hot asthenosphere. The findings for long-periods of Love and Rayleigh waves in Armenia reveal an almost uniform velocity distribution pattern, along with very-low-velocity anomalies in the uppermost mantle, attributed to a thin lithosphere or the lack of lithospheric mantle in most units of the study region. The results demonstrate an acceptable accordance with known geological features in the Eurasian-Arabic ongoing collision zone. Overall, the main strengths of this paper lie in the application of a tomographic technique utilizing an important data set. The findings have the potential to provide new insights into the Armenian region and its surrounding areas.
Russian Journal of Earth Sciences. 2025;25(3):ES3014
ES3014


