Gediminas Hill Slopes Behavior in 3D Finite Element Model

: In this manuscript, we present the results of three-dimensional ﬁnite element analyses performed in the cloud of a large-scale model of the Gediminas Hill buildings and the construction remains of the Upper and Lower castles in Vilnius city. The greatest challenge associated with the simulated numerical model is the difﬁcult geological layer surface inclinations and soil–structure interaction behavior prognosis, which require signiﬁcant computational resources. The purpose of this research work is to present current and possible worst-case scenarios for Gediminas Hill, considering its buildings and construction remains, regarding the stability of its slopes through a safety analysis. The construction of a numerical three-dimensional model of Gediminas Hill allows for us to assess the soil–structure interaction behavior. The results of non-linear analysis on the created model are in agreement with the tendencies observed in direct geodetic measurements and the relevant landslide history.


Introduction
Numerical modeling techniques have evolved with the improved use of computational resources. At first, two-dimensional (2D) numerical modeling was created, while three-dimensional (3D) numerical modeling was introduced later. Initially, when solving geotechnical problems by applying 3D modeling [1], only small-scale volumes with very limited finite element quantities or large-scale volumes with enormous finite elements could be used; in this regard, Kausel [2] has summarized the history of soil-structure interaction modeling very well. At present, 3D models are larger, more complex, and more computationally reliable [3][4][5][6]. The progress of high-performance calculation methods has allowed for the performance of stochastic finite element analyses using large-scale meshes or a large amount of Monte Carlo simulations [7,8].
Three-dimensional geotechnical models are more suitable for large projects [9]. When such projects are analyzed, it is easier to spend more time on 3D numerical model preparation and obtain the results at any necessary area or point than to conduct many 2D calculations in different cross-sections and obtain results that are only valid in the chosen 2D simulation cross-section. Furthermore, in 2D slope stability analysis, the landslide width cannot be obtained. This is why 3D finite element modeling is very important for slope stability analysis [10]. Additionally, 3D geotechnical modelling can help to define representative borehole logs indicating the stratification of soil [11], the geometry of different sedimentary layers [12], and whether it is necessary to add the superstructure to the model [13]. The 3D geotechnical model can help to determine the weakest areas in the model by analyzing the reduction in safety [14].
One of the most interesting and actual objects-from a geotechnical and geological engineering point of view-in Lithuania, to which a 3D model can be applied, is Gediminas Hill and the remains of its Upper castle in Vilnius city. This object is very important to the local community, as it is included in the Lithuanian culture heritage register. Various studies on Gediminas Hill have already been carried out [15][16][17][18][19][20][21][22][23]. The purpose of this research is to present the current and possible worst-case scenarios of Gediminas Hill and the remains of the Upper castle, regarding the stability of slopes, by conducting a factor of safety analysis. The reliability of the 3D geotechnical model is ensured by validating obtained results with respect to the existing monitoring data.

Input Data
This section introduces the initial data entered into the 3D model, which consisted of geological engineering and geotechnical conditions; the geodetic surface; building remains and structure plans; pedestrian pathways; and tunnels. All these data were prepared for entry into 3D model, which was calculated using a VILNIUS TECH virtual computer.

Computational Resources
The numerical modeling was simulated in a private VILNIUS TECH cloud. Cloud services were chosen to quickly respond to the changing needs of the computing resource. This was the first time that such large-scale calculations were performed using the university's private computing resources, making it difficult to determine the associated resource needs. A physical server was prepared in the VILNIUS TECH data center for these calculations, which was connected to the VILNIUS TECH private cloud. The cloud service uses opensource OpenStack software that allows for the construction of private clouds. The physical server has 2 Intel Xeon Gold 5220R 2.2 GHz processor units, 1536 GB RAM, and uses 6 TB space, which was combined from 8 SSD drives. For the first calculations with the virtual server, 20 core CPU resources and 320 GB of RAM were assigned; however, the calculation capacity was later increased to 1015 GB RAM. To ensure successful numerical simulations with the DIANA finite element analysis (FEA) release 10.3 [24], the virtual server was connected to the department of Information systems resource monitoring system, which allowed for the server loads to be monitored and faults in the numerical simulation to be quickly identified.

Geodetic Data
The investigated surface of the Gediminas Hill area was scanned and a 3D surface model was created from simple photographs using the Acute3D Viewer [25] software. To ensure the high resolution of the 3D surface model, photogrammetry, computer vision, and computational geometry algorithms were used, fulfilling industrial quality requirements in terms of precision, scalability, efficiency, usage, robustness, and interoperability. A general view of the Gediminas Hill 3D surface is presented in Figure 1. One of the most interesting and actual objects-from a geotechnical and geological engineering point of view-in Lithuania, to which a 3D model can be applied, is Gediminas Hill and the remains of its Upper castle in Vilnius city. This object is very important to the local community, as it is included in the Lithuanian culture heritage register. Various studies on Gediminas Hill have already been carried out [15][16][17][18][19][20][21][22][23]. The purpose of this research is to present the current and possible worst-case scenarios of Gediminas Hill and the remains of the Upper castle, regarding the stability of slopes, by conducting a factor of safety analysis. The reliability of the 3D geotechnical model is ensured by validating obtained results with respect to the existing monitoring data.

Input Data
This section introduces the initial data entered into the 3D model, which consisted of geological engineering and geotechnical conditions; the geodetic surface; building remains and structure plans; pedestrian pathways; and tunnels. All these data were prepared for entry into 3D model, which was calculated using a VILNIUS TECH virtual computer.

Computational Resources
The numerical modeling was simulated in a private VILNIUS TECH cloud. Cloud services were chosen to quickly respond to the changing needs of the computing resource. This was the first time that such large-scale calculations were performed using the university's private computing resources, making it difficult to determine the associated resource needs. A physical server was prepared in the VILNIUS TECH data center for these calculations, which was connected to the VILNIUS TECH private cloud. The cloud service uses open-source OpenStack software that allows for the construction of private clouds. The physical server has 2 Intel Xeon Gold 5220R 2.2 GHz processor units, 1536 GB RAM, and uses 6 TB space, which was combined from 8 SSD drives. For the first calculations with the virtual server, 20 core CPU resources and 320 GB of RAM were assigned; however, the calculation capacity was later increased to 1015 GB RAM. To ensure successful numerical simulations with the DIANA finite element analysis (FEA) release 10.3 [24], the virtual server was connected to the department of Information systems resource monitoring system, which allowed for the server loads to be monitored and faults in the numerical simulation to be quickly identified.

Geodetic Data
The investigated surface of the Gediminas Hill area was scanned and a 3D surface model was created from simple photographs using the Acute3D Viewer [25] software. To ensure the high resolution of the 3D surface model, photogrammetry, computer vision, and computational geometry algorithms were used, fulfilling industrial quality requirements in terms of precision, scalability, efficiency, usage, robustness, and interoperability. A general view of the Gediminas Hill 3D surface is presented in Figure 1.  The generated geodetic model point cloud was very dense and not suitable for 3D modelling. It was necessary to prepare simplifications of the geodetic data, as follows: for slopes and pedestrian pathways, the surface point density must be no less than 0.5 m, and the construction surface density must be no less than 0.2 m. Furthermore, it was very important to know where the construction edge ends and the soil starts; therefore, the contact points of different objects also were included into the geodetic surface model. According to accepted simplifications of the geodetic data, a three-dimensional surface of Gediminas Hill was created ( Figure 2). The generated geodetic model point cloud was very dense and not suitable for 3D modelling. It was necessary to prepare simplifications of the geodetic data, as follows: for slopes and pedestrian pathways, the surface point density must be no less than 0.5 m, and the construction surface density must be no less than 0.2 m. Furthermore, it was very important to know where the construction edge ends and the soil starts; therefore, the contact points of different objects also were included into the geodetic surface model. According to accepted simplifications of the geodetic data, a three-dimensional surface of Gediminas Hill was created ( Figure 2).

Engineering Geological and Geotechnical Conditions
From a geomorphological perspective, Gediminas Hill is located on the north-western edge of the Medininkai Highland, part of the Ašmena Upland, formed during the Middle Pleistocene Glaciation [26]. According to [27], the geological structure of the Gediminas Hill is represented by deposits from a few Pleistocene glaciations and, according to the Quaternary period stratigraphic scheme currently approved for use in the Lithuanian Geological Survey, the deposits depend on Dainava (Elsterian 2) and Žeimena formations of the Middle Pleistocene [28]. The Žeimena Formation can be subdivided into deposits of the Žemaitija (Saalian) and Medininkai (Warthian) sub-formations. The Middle Pleistocene deposits are affected by deluvial processes and human activities on the surface and slopes of the hill [27].
The last major engineering geological and geotechnical investigations of Gediminas Hill were carried out in 2019 and 2020 [29]. Detailed hydrogeological conditions of the site have been presented by [27]. During these investigations, geological engineering layers, along with their physical and mechanical properties, were identified (Table 1).

Engineering Geological and Geotechnical Conditions
From a geomorphological perspective, Gediminas Hill is located on the north-western edge of the Medininkai Highland, part of the Ašmena Upland, formed during the Middle Pleistocene Glaciation [26]. According to [27], the geological structure of the Gediminas Hill is represented by deposits from a few Pleistocene glaciations and, according to the Quaternary period stratigraphic scheme currently approved for use in the Lithuanian Geological Survey, the deposits depend on Dainava (Elsterian 2) and Žeimena formations of the Middle Pleistocene [28]. The Žeimena Formation can be subdivided into deposits of the Žemaitija (Saalian) and Medininkai (Warthian) sub-formations. The Middle Pleistocene deposits are affected by deluvial processes and human activities on the surface and slopes of the hill [27].
The last major engineering geological and geotechnical investigations of Gediminas Hill were carried out in 2019 and 2020 [29]. Detailed hydrogeological conditions of the site have been presented by [27]. During these investigations, geological engineering layers, along with their physical and mechanical properties, were identified (Table 1).
Seeking to simplify the 3D geotechnical model, the geological layers given in Table 1 were formed by joining some layers with the same stratigraphy index; that is, some geological layers were averaged in height and their properties were determined based on existing local geotechnical experience. The simplified geological layer properties used in the 3D geotechnical model are given in Table 2. The dilatancy angle ( Table 2) was also accepted, according to [30], where the dilatancy of clay (except for over-consolidated clay) is mostly ψ ≈ 0 • , while the dilatancy of sand depends on both the density and the angle of internal friction, and can be found using the equation

Buildings and Structures Remains, Path, Funicular and Tunnels
Ten construction remains objects, a pedestrian path, two tunnels, and one funicular are situated in the area of Gediminas Hill ( Figure 3). All the buildings and construction remains were included in the numerical 3D model. The top of Gediminas Hill can be reached by a pedestrian pathway or funicular (Figure 3). At present, visitors are not allowed to see the tunnels, which were installed during World War I and World War II. One tunnel was excavated vertically from the top of Gediminas Hill (depth~37-40 m) and has two exit ways (to the south-western and north-western slopes). The south-western tunnel length is~88 m and the north-western tunnel length is~79 m. This tunnel, due to its age and the rotten timber construction damages that occurred in 1960 [33], was filled with a liquid sand mixture to avoid tunnel construction failure. Another tunnel, which is close to the castle-keeper's building [34], is quite short (~13 m), but has not been filled with liquid sand mixture. It is worth noting that this tunnel was only found at the end of 2019, at which point its timber construction was in a very bad condition; therefore, during 2020, bearing construction reconstruction of this tunnel was carried out. reached by a pedestrian pathway or funicular (Figure 3). At present, visitors are not allowed to see the tunnels, which were installed during World War I and World War II. One tunnel was excavated vertically from the top of Gediminas Hill (depth ~37-40 m) and has two exit ways (to the south-western and north-western slopes). The south-western tunnel length is ~88 m and the north-western tunnel length is ~79 m. This tunnel, due to its age and the rotten timber construction damages that occurred in 1960 [33], was filled with a liquid sand mixture to avoid tunnel construction failure. Another tunnel, which is close to the castle-keeper's building [34], is quite short (~13 m), but has not been filled with liquid sand mixture. It is worth noting that this tunnel was only found at the end of 2019, at which point its timber construction was in a very bad condition; therefore, during 2020, bearing construction reconstruction of this tunnel was carried out.

Development of a Numerical Three-Dimensional Model of Gediminas Hill
According to the existing investigations of Gediminas Hill buildings and structure remains, as well as slope-monitoring, it was observed that shallow landslides form on the slopes. For this reason, in 3D numerical simulations, the Mohr-Coulomb material model was accepted for soil, while linear constitutive material models were used for buildings and remains. After analysis, the relevant literature, existing drawings, and information sources, and the geometry and locations of certain objects, were determined, and a three-

Development of a Numerical Three-Dimensional Model of Gediminas Hill
According to the existing investigations of Gediminas Hill buildings and structure remains, as well as slope-monitoring, it was observed that shallow landslides form on the slopes. For this reason, in 3D numerical simulations, the Mohr-Coulomb material model was accepted for soil, while linear constitutive material models were used for buildings and remains. After analysis, the relevant literature, existing drawings, and information sources, and the geometry and locations of certain objects, were determined, and a threedimensional numerical model of Gediminas Hill was created. The dimensions of the numerical model were 270 m in the east-west direction and 228 m in the south-north direction. The height of the model-from the horizontal supporting plane to the top of the west tower-was 76.47 m. The distance from the ground surface to the horizontal plane at different angles of the model varied from 7.9 m to 15.4 m. The lower horizontal plane of the numerical model (as well as the lateral surrounding surfaces) was constrained in all directions. In total, 1,847,241 finite elements were used, the size of which was not constant.
A numerical model of Gediminas Hill, comprising the soil layers forming the mountain, as well as the buildings and structures placed on it (Figure 3), was created in the finite element program Diana 10.3. The hill array consisted of eight major soil layers (Table 2), classified as solid volumetric bodies. Surfaces were used to separate the soil layers, which were composed of 3D Face-type elements interconnected in three-dimensional (x, y, z) space (as in Figure 2). The approximate dimension of one 3D Face element was 3.0 × 3.0 m. The three-dimensional face surfaces and their nodes, created in the drawing program, were localized using the Lithuanian coordinate system 94 (LKS94) coordinate system, and the heights were in correspondence with the Lithuanian elevation system (LAS) system.
The behavior of the numerical model and slope stability were determined through a four-stage calculation analysis. All four stages were assigned to three phases. In the first stage, the self-weight of the natural soil layers and tunnel backfill were evaluated (Figure 4, Phase 1). In the second stage ( Figure 4, Phase 2), the self-weight of the technogenic soil layer (tIV-dIV) and gabions in the northern slope part (Table 3) were evaluated. The gabions were installed on natural soil layers to replace technogenic soil due to the occurrence of a shallow landslide in 2016 [18]. The third stage (Figure 4, Phase 3) involved evaluating the construction and self-weight of buildings, as well as additional loads on the pedestrian pathway ( Figure 3) and the remains of the Upper castle pitch (Figure 3). In the fourth stage ( Figure 4, Phase 3), slope stability calculations assessing the reduction in the strength of the technogenic soil layer saturation were initiated [35].  (Table 2), classified as solid volumetric bodies. Surfaces were used to separate the soil layers, which were composed of 3D Face-type elements interconnected in three-dimensional (x, y, z) space (as in Figure 2). The approximate dimension of one 3D Face element was 3.0 × 3.0 m. The three-dimensional face surfaces and their nodes, created in the drawing program, were localized using the Lithuanian coordinate system 94 (LKS94) coordinate system, and the heights were in correspondence with the Lithuanian elevation system (LAS) system.
The behavior of the numerical model and slope stability were determined through a four-stage calculation analysis. All four stages were assigned to three phases. In the first stage, the self-weight of the natural soil layers and tunnel backfill were evaluated ( Figure  4, Phase 1). In the second stage ( Figure 4, Phase 2), the self-weight of the technogenic soil layer (tIV-dIV) and gabions in the northern slope part (Table 3) were evaluated. The gabions were installed on natural soil layers to replace technogenic soil due to the occurrence of a shallow landslide in 2016 [18]. The third stage (Figure 4, Phase 3) involved evaluating the construction and self-weight of buildings, as well as additional loads on the pedestrian pathway ( Figure 3) and the remains of the Upper castle pitch (Figure 3). In the fourth stage ( Figure 4, Phase 3), slope stability calculations assessing the reduction in the strength of the technogenic soil layer saturation were initiated [35].   A non-linear numerical analysis of slopes displacements was conducted using different material models for soils, structures, and building remains. For the soil layers, the Mohr-Coulomb plasticity material model was applied, the properties of which are given in Tables 2 and 3. Additional properties are given in Table 3, detailing the properties of the accepted gabions and gabions surrounding soil in the northern slope. For the simulated constructions and building remains, elastic body mechanical properties were assigned, which represent the behavior of stone masonry and bricks up to their plastic limits. Structures and buildings transferred loads to the soil layers. The mechanical properties of the evaluated structures and building remains are presented in Table 4. The remains of the Old Arsenal building were simulated with low density (Table 4), as the construction of this building was not investigated by the authors; therefore, the volume of internal spaces was evaluated to reduce the total masonry density.  Before conducting the numerical simulation stages and phases (Figure 4), the numerical model must be filled with volumetric finite elements. All the finite element volumes were divided into hexa/quad-type elements, where, for total element volume coverage, tetrahedral/triangle-type elements can additionally be used. The size of a finite element depends on the meshed volumetric body edge length, as the length is divided into two parts. If the division number is increased, the required computational resource quantity automatically increases; thus, it is necessary to evaluate the capacity of computational resources based on the difficulty associated with the simulated geometry. The calculation time also depends on the quantity of the finite elements.
In the third calculation phase (Figure 4), all volumetric finite elements are active. When the self-weight evaluation reaches 100% in the third phase, strength reduction analysis is only initiated for the tIV-dIV soil layer. For other layers, the reduction in safety is not assessed.
The strength reduction method (SRF), sometimes referred to as the strength reduction finite element method [36], is a technique for calculating the slope factor of safety (FoS) [37], defined as the ratio of current soil shear strength to the minimum shear strength necessary to avoid failure. Alternatively, the FoS is considered as the factor required to reduce the soil strength to produce imminent failure. In SRF, the strength characteristics of the soil materials (cohesion and tangent angle of internal friction) are reduced by a factor until the loss of stability or failure of the structure occurs, as follows [38,39]: If the factor of safety calculated by Formula (1) does not converge, then a new lower increment is evaluated to obtain a new factor of safety, which reaches convergence: If this is not converged after the application of Formula (2), then the results are evaluated using the last convergence factor: At first, the strength reduction analysis factor of safety is equal to 1.0: With each chosen factor, the soil cohesion (5) and tangent angle of internal friction are reduced (6): tan ϕ n+1 = tan ϕ n FS n+1 (6)

Behavior of Numerical Model
The interaction behavior of finite elements in a numerical model depends on their physical and mechanical properties, material models, supports or loads, finite element size and quantity [40], and slope width [41], among other factors. Historically, the behavior of Gediminas Hill has depended on additional factors, including climate change and the history of landslides [33]. All these factors increase the required computational resources. The considered model of Gediminas Hill demanded a calculation capacity of up to 1015 GB RAM. Processing such a numerical model can take from 25 to 26 h, of which from 4 to 5 h are required for the generation of the finite element mesh and from 20 to 22 h for numerical analysis. The finite element generation time depends on the complexity of the model geometry and the set of mesh parameters. The numerical analysis time also depends on the number of finite elements that are generated and the development of displacements. If larger displacements occur, the calculation process requires more iterations and, with each safety analysis step, the number of iterations increases. Therefore, the numerical analysis in the third phase took the longest. The obtained results on slope displacements from the third simulation phase without a strength reduction analysis (Figure 4) are presented in Figure 5.
At first, the strength reduction analysis factor of safety is equal to 1.0: With each chosen factor, the soil cohesion (5) and tangent angle of internal friction are reduced (6): , tan .

Behavior of Numerical Model
The interaction behavior of finite elements in a numerical model depends on their physical and mechanical properties, material models, supports or loads, finite element size and quantity [40], and slope width [41], among other factors. Historically, the behavior of Gediminas Hill has depended on additional factors, including climate change and the history of landslides [33]. All these factors increase the required computational resources. The considered model of Gediminas Hill demanded a calculation capacity of up to 1015 GB RAM. Processing such a numerical model can take from 25 to 26 h, of which from 4 to 5 h are required for the generation of the finite element mesh and from 20 to 22 h for numerical analysis. The finite element generation time depends on the complexity of the model geometry and the set of mesh parameters. The numerical analysis time also depends on the number of finite elements that are generated and the development of displacements. If larger displacements occur, the calculation process requires more iterations and, with each safety analysis step, the number of iterations increases. Therefore, the numerical analysis in the third phase took the longest. The obtained results on slope displacements from the third simulation phase without a strength reduction analysis ( Figure  4) are presented in Figure 5.  The displacements presented in Figure 5 showed the same tendency as the monitoring data (presented in Figure 6). Thus, the compiled 3D numerical model can present the same slope movement tendencies. It is very problematic to obtain authentic results when historic object behavior is simulated, as the buildings in the model were built of bricks and stones in 1419, partially demolished in 1655-1661, restored and conserved in 1816-1817, 1930-1939, and 1956-1963, and conserved and partially restored in 1990-1993 [42]. Thus, the biggest concentration of soil displacements was close to the remains of the Upper castle western tower and the Upper castle remains ( Figure 5). Seeking to develop a prognosis for the possible formation of landslide areas, strength reduction calculations were carried out in the third simulation phase (Figure 4). The strength reduction analysis was applied only to the tIV-dIV soil layer (Table 2), as existing knowledge and monitoring data [43] have shown that landslides only form in the tIV-dIV layer. These landslides are shallow.  The displacements presented in Figure 5 showed the same tendency as the monitoring data (presented in Figure 6). Thus, the compiled 3D numerical model can present the same slope movement tendencies. It is very problematic to obtain authentic results when historic object behavior is simulated, as the buildings in the model were built of bricks and stones in 1419, partially demolished in 1655-1661, restored and conserved in 1816-1817, 1930-1939, and 1956-1963, and conserved and partially restored in 1990-1993 [42]. Thus, the biggest concentration of soil displacements was close to the remains of the Upper castle western tower and the Upper castle remains ( Figure 5). Seeking to develop a prognosis for the possible formation of landslide areas, strength reduction calculations were carried out in the third simulation phase (Figure 4). The strength reduction analysis was applied only to the tIV-dIV soil layer (Table 2), as existing knowledge and monitoring data [43] have shown that landslides only form in the tIV-dIV layer. These landslides are shallow. The displacements presented in Figure 5 showed the same tendency as the monitoring data (presented in Figure 6). Thus, the compiled 3D numerical model can present the same slope movement tendencies. It is very problematic to obtain authentic results when historic object behavior is simulated, as the buildings in the model were built of bricks and stones in 1419, partially demolished in 1655-1661, restored and conserved in 1816-1817, 1930-1939, and 1956-1963, and conserved and partially restored in 1990-1993 [42]. Thus, the biggest concentration of soil displacements was close to the remains of the Upper castle western tower and the Upper castle remains ( Figure 5). Seeking to develop a prognosis for the possible formation of landslide areas, strength reduction calculations were carried out in the third simulation phase (Figure 4). The strength reduction analysis was applied only to the tIV-dIV soil layer (Table 2), as existing knowledge and monitoring data [43] have shown that landslides only form in the tIV-dIV layer. These landslides are shallow.  Strength reduction analysis of the technogenic soil (tIV-dIV) layer ( Figure 7) allowed for us to determine areas with the highest concentrations of displacements. Strength reduction analysis indicated not only displacements of the technogenic soil layer (tIV-dIV), but also displacements of constructions and buildings (in some places, the technogenic soil layer is under the foundations). In each SRF calculation, the depth of the technogenic soil (tIV-dIV) layer is different due to the soil displacements, where the newly formed soil mass is evaluated. When shallow landslides form, they include only technogenic soil (tIV-dIV). The increasing mass and volume of the moving soil causes stability failure, and a local shallow landslide occurs. Strength reduction analysis of the technogenic soil (tIV-dIV) layer ( Figure 7) allowed for us to determine areas with the highest concentrations of displacements. Strength reduction analysis indicated not only displacements of the technogenic soil layer (tIV-dIV), but also displacements of constructions and buildings (in some places, the technogenic soil layer is under the foundations). In each SRF calculation, the depth of the technogenic soil (tIV-dIV) layer is different due to the soil displacements, where the newly formed soil mass is evaluated. When shallow landslides form, they include only technogenic soil (tIV-dIV). The increasing mass and volume of the moving soil causes stability failure, and a local shallow landslide occurs.  According to Figure 7, the highest possibility of shallow landslide formation was observed in the eastern and south-western slopes. The formation of these landslides correlates with tIV-dIV soil layer thickness-a thicker layer increases the landslide possibility and vice versa. A shallow landslide in the eastern slope occurred on 07 February 2022 (Figure 8), which appeared under the most difficult environmental conditions: the top layer of the soil was not frozen, and there was a lot of melting water under snow cover. According to Figure 7, the highest possibility of shallow landslide formation was observed in the eastern and south-western slopes. The formation of these landslides correlates with tIV-dIV soil layer thickness-a thicker layer increases the landslide possibility and vice versa. A shallow landslide in the eastern slope occurred on 07 February 2022 (Figure 8), which appeared under the most difficult environmental conditions: the top layer of the soil was not frozen, and there was a lot of melting water under snow cover.
In Figure 8, the assigned tags correspond to those given in Figure 3. The landslide location and size (width, 6-6.5 m; length, 11-11.5 m) correlated well with the displacement zone obtained in the numerical simulation (Figure 7). In Figure 8, the assigned tags correspond to those given in Figure 3. The landslide location and size (width, 6-6.5 m; length, 11-11.5 m) correlated well with the displacement zone obtained in the numerical simulation (Figure 7).

Conclusions
To create a volumetric numerical model of Gediminas Hill, it was necessary to utilize various recent and historical data. The constructed numerical model allowed for us to evaluate the slope stability behavior of Gediminas Hill, based on the accepted material models, loads, calculation phases, and so on, from which we drew the following conclusions: 1. Constructions, buildings, and additional loads increased soil displacement, mostly in the technogenic (tIV-dIV) soil layer. The analyzed monitoring data, including recent and historical investigation reports on the constructions, as well as the numerical model behavior, allowed us to assume that the deformation of building remains depends on technogenic (tIV-dIV) soil layer displacements. The influence of the climate (i.e., precipitation) on slope stability was also evaluated, considering strength reductions, which demonstrated that the slope stability depends on the technogenic soil layer height and the inclination of undisturbed soil layers. 2. The surface inclination of natural soil layers at the point of contact with the technogenic (tIV-dIV) soil layer in some areas is higher than that of the technogenic soil layers. In these areas, the biggest deformation zones may form. Numerical model analysis allowed for us to make the conclusion that technogenic (tIV-dIV) soil layer shallow landslide occurrence depends on the layer height, changes in mechanical parameters, and the surface inclination of the natural soil layers. 3. The numerical Gediminas Hill 3D model allowed for us to assess the soil-building interaction behaviors. The results of non-linear analysis of the created model corresponded to the tendencies observed in direct geodetic measurements and were in agreement with the relevant landslide history.

Conclusions
To create a volumetric numerical model of Gediminas Hill, it was necessary to utilize various recent and historical data. The constructed numerical model allowed for us to evaluate the slope stability behavior of Gediminas Hill, based on the accepted material models, loads, calculation phases, and so on, from which we drew the following conclusions:

1.
Constructions, buildings, and additional loads increased soil displacement, mostly in the technogenic (tIV-dIV) soil layer. The analyzed monitoring data, including recent and historical investigation reports on the constructions, as well as the numerical model behavior, allowed us to assume that the deformation of building remains depends on technogenic (tIV-dIV) soil layer displacements. The influence of the climate (i.e., precipitation) on slope stability was also evaluated, considering strength reductions, which demonstrated that the slope stability depends on the technogenic soil layer height and the inclination of undisturbed soil layers.

2.
The surface inclination of natural soil layers at the point of contact with the technogenic (tIV-dIV) soil layer in some areas is higher than that of the technogenic soil layers. In these areas, the biggest deformation zones may form. Numerical model analysis allowed for us to make the conclusion that technogenic (tIV-dIV) soil layer shallow landslide occurrence depends on the layer height, changes in mechanical parameters, and the surface inclination of the natural soil layers. 3.
The numerical Gediminas Hill 3D model allowed for us to assess the soil-building interaction behaviors. The results of non-linear analysis of the created model corresponded to the tendencies observed in direct geodetic measurements and were in agreement with the relevant landslide history.