###### Abstract  
The introduction of modern methods for the mathematical processing of geological data is one of the promising areas of study and development in the field of geosciences. For example, today mathematical geology makes it possible to reliably identify astronomical cycles by measuring the scalar magnetic parameters of rocks (magnetic susceptibility). The main aim of this study is to develop a mathematical tool for identifying stable oscillation cycles (periods) in the dataset of the magnetic susceptibility of rocks in a geological section. The author's method (algorithm) is based on the concept of discrete mathematical analysis--an innovative mathematical approach to the analysis of discrete geological and geophysical data. Its reliability is also demonstrated, by comparison with the results obtained by classical methods: Fourier analysis, Lomb periodogram, and REDFIT. The proposed algorithm was applied by the authors to analyze the material of field geological studies of the Zhelezny Rog section (Taman Peninsula). As a result, stable cycles were determined for the Pontian and Lower Maeotian sedimentary strata of the Black Sea Basin (Paratethys). s 2022 2  
stratigraphic sequence. In particular, cyclostratigraphy involves applying this knowledge to geochronology by increasing the accuracy and resolution of chronostratigraphic units. It uses the astronomical cycles of the currently known periodicity and an interpretation of the sedimentation conditions. According to the astronomical theory of paleoclimates, largely developed by [PERSON], cyclical changes in the Earth's orbit (orbital cycles) are the main causes of climate change [2]. These cycles translate into climatic, oceanographic, sedimentary, and biological changes recorded in sedimentary deposits over geological time. Numerous case studies have shown that detailed analysis of sedimentary rocks allows these cycles to be identified with a high degree of confidence. An unprecedented high temporal resolution is available once the relationship between the sedimentary record and orbital forcing is established. This relationship provides a basis for the timing of the processes occurring in the Earth system [3]. The main task of cyclostratigraphy is to analyze the structure of sedimentary deposits and identify stable signals reflecting the orbital influence. An essential aspect of time series analysis in cyclostratigraphy is transforming raw data into the time domain using calibration points. Once the time series is established, mathematical methods can be used to detect orbital periodicities. Determination of astronomical cyclicity using geochronological methods, including paleomagnetic reconstruction and lithological analysis, is possible by defining the magnetic susceptibility of sediments. Magnetic susceptibility depends on the amount of solar radiation reaching the Earth's surface and reflects climatic fluctuations. Spectral analysis of scalar magnetic parameters allows the detection of long-term period oscillations of the Earth's axis, the angle of inclination of the Earth's axis to the plane of the ecliptic, and eccentricity. This method gives a possibility to define absolute ages of sediments with accuracy in the order of 20,000 to 400,000 years. In this paper, the application of mathematical methods to solve some critical problems of cyclostratigraphy and paleogeography of Eastern Paratethys Maeotian sediments is proposed. It should be noted that Eastern Paratethys has, for many years, been the research subject of one of the authors of the present paper. The study is aimed at the analysis of sedimentation processes of the Eastern Paratethys by the identification of stable oscillation cycles in the magnetic susceptibility data. To identify such cycles (periods), the authors developed a mathematical algorithm that is based on the concept of discrete mathematical analysis (DMA).