With the background of global climate change on earth, natural hazards such as landslides, debris/mudflows, floods, erosion, mudflows, landslides are happening more frequent [1] - [9] . There are many such sensitive regions and areas in Georgia. One of the most important of these is riv. Tetri Aragvi gorge and the various dangers emanating from its slopes. River has many large and small tributaries. Under the hazard risk are population and road section of Natakhtari-Mleta of Mtskheta-Stepantsminda-Larsi road [10] - [16] .

These valleys are periodically characterized by strong mudflows. Statistics shows that number of debris/mudflow and landslide events increases annually. Therefore, conducting engineering and geological surveys and identifying geodynamic hazards on the Natakhtari-Mleta section of Mtskheta-Stepantsminda-Larsi road is of great importance for the security and sustainable economic development of our country.

Aim of the research was to study geodynamic hazards activated on the right slope of riv. Tetri Aragvi, in the territory of the district Pasanauri, which may pose a danger to traffic on the Natakhtari-Mleta section of the Georgian Military Road (Mtskheta-Stepantsminda-Larsi) ( Figure 1 ). In particular, we are talking about mudflows that periodically emerge from the aforementioned ravine. If we consider the problem raised in a historical aspect, we will see that these events are not peaceful and harmless in nature. The overall geomorphological picture of the ravine clearly shows a fairly large extraction cone formed at its base ( Figure 2 ), which, at first glance, is much larger than the scale of the ravine, which indicates its very active past ( Figure 3 ). This does not give reason for complacency, since in the recent past, several disasters occurred in Georgia. This is also due to the fact that a fairly large landslide body has formed in the headwaters of the ravine we selected for our research, which includes a proluvial massif formed in the ravine bed at an early stage (pQ_IV). This massif is experiencing severe erosion and feeds the ravine with solid material to produce mudflows. It is built of gravel and crushed stone material, bound together with clay-loam filler. Its thickness reaches 22.0 meters. It is bordered on three sides by flat forms of vertical arrangement and forms an accumulation-erosion form. It is easy to see that the washing away and collapse of such a formation represents a real threat, being a direct triggering factor for mudflows.

As our study of the data from field-scientific research and familiarization with field data from early years has shown, the factors contributing to the occurrence of mudflows in the above-mentioned valleys, in addition to those mentioned above, are the following: Climate; geological structure and geomorphological features, to which, in this case, consideration of seismic factors must also be added [17] . Study area is characterized by diverse climate conditions, main triggering factor for debris/mudflow activation is daily precipitation [18] .

Among the geological factors, the lithological structure of the study area is important. Study area is located in the Mestia-Tianeti tectonic zone and is represented

by carbonate-terrigenous flysch of the Berrian-Valanginian and Hotrivian stages. These floors are built up of clastic limestones and sandstones, turbidites, pelagic marls, limestones, mudstones, and shales. Rocks are easily subject to weathering processes [19] which creates an important prerequisite for the occurrence of modern gravitational phenomena. These phenomena are actively developing in the headwaters of these valleys, where the main mudflows are formed [20] ( Figure 4 ).

Among the geomorphological features of the study area, the steepness of the slopes and the high degree of relief fragmentation should be highlighted [19] [20] ( Figures 5-7 ).

And finally, the most important thing—all the above-mentioned factors unambiguously serve to generate landslide-gravity events within the mentioned valleys, which serve as the main supplier of solid material for mudflows [19] [20] ( Figures 5-7 ). This is clearly illustrated in Table 1 below.

In this regard, rock samples were taken and their physical characteristics were determined. Their mechanical characteristics were selected based on existing normative documents. On the basis of topographic and geological data landslide slide surface was identified and schematic figure was created. Later slope stability K coefficient was calculated, considering various natural factors (seismicity and rock moisture) [19] [20] .

Based on the analyzed data, it can also be noted that the formation of debris/mudflows of various intensity and volume is expected within the study area, which will lead to an increase in ecological risks and, accordingly, the zones of mudflow impact and the vulnerability of various infrastructure facilities (residential buildings, highways, etc.) located in these zones. In the next stage of the research, a digital elevation model (DEM) (Digital Elevation Model—2 × 2 m resolution) will be created using modern technologies of remote sensing and geoinformation systems, which will be used to determine the slope and slope aspect of the study ravine and to construct profiles ( Figure 8 ). The data will be processed using modern software to model landslide and mudflow events and identify their characteristics [12] [21] .

In Table 2 are the results of the calculations performed [22] - [24] :

It should be noted that in future it is planned to use these data for mudflow simulation using software RAMMS ( Figure 9 , Figure 10 ). Which was developed to simulate the runout of mudflow and debrisflows in complex terrain [21] [25] [26] .