The Enguri dam and water reservoir, nested in the southwestern Caucasus
(Republic of Georgia), are surrounded by steep mountain slopes. At a distance
of 2.5 km from the dam, a mountain ridge along the reservoir is affected by
active deformations with a double vergence. The western slope, directly
facing the reservoir, has deformations that affect a subaerial area of
1.2 km
GPS data and plate tectonic models indicate that the Greater and Lesser
Caucasus are tectonically very active, with ongoing mountain building
processes comprising complex deformation with vertical and horizontal strain
partitioning (Rebai et al., 1993; Koçyiğit et al., 2001; Reilinger et
al., 1997, 2006; Tan and Taymaz, 2006; Pasquaré et al., 2011). These
processes result from the convergence and continent–continent collision that is still developing between the Eurasian and African–Arabian plates
(Avagyan et al., 2010; Adamia et al., 2017) (Fig. 1). The active deformation
is accompanied by diffuse seismicity that reaches
In the southwestern part of the Georgian Greater Caucasus, the Enguri dam (the world's sixth highest dam) (Figs. 2 and 3a) and hydroelectric plant presently provide a major part of the energy to the country. The Enguri dam is located near a series of geological features that suggest Quaternary uplift, comprising deeply entrenched rivers, deformed river terraces, recent faults and folds and very steep slopes (Tibaldi and Tsereteli, 2017; Tibaldi et al., 2017a). The strategic importance of the Enguri hydroelectric plant and the fact that it is nested within one of the most tectonically active regions of the Caucasus suggest the necessity of carrying out modern studies on the stability of the slopes surrounding the water reservoir, 19 km long, which are completely missing in the scientific literature (only one technical consultant report exists, which is unpublished – CGS, 2015). Mountain slopes, in fact, are subject to gravity effects that can be enhanced by earthquake shaking (Gutierrez-Santolalla et al., 2005); it has been widely demonstrated that this shaking can cause the generation of shallow landslides and deep-seated gravity slope deformations (DSGSDs) in active seismic areas (Beck, 1968; Solonenko, 1977; McCleary et al., 1978; Radbruch-Hall, 1978; Tibaldi et al., 1995; McCalpin, 1999). The dynamic loading produced by seismicity may trigger discrete episodes of fast slope deformation (e.g., Beck, 1968) or may accelerate already existing movements along DSGSD structures (e.g., Pasuto and Soldati, 1996). Whether the movements are continuous (Varnes et al., 1989, 1990, 2000) or intermittent (Beget, 1985), the slopes under active deformation pose a threat to infrastructures (Mahr, 1977; Radbruch-Hall, 1978; McCalpin and Irvine, 1995). Near to the Enguri reservoir is the main road leading to the famous town of Mestia, hosting several ancient monuments and a ski resort. This road is affected by important slope deformations in correspondence with the same slope that faces the Enguri water reservoir. Moreover, in general, slope movements can turn into dangerous catastrophic collapses, examples of which have been reported in Japan (Chigira and Kiho, 1994), Italy (Semenza and Ghirotti, 2000) and Canada (Evans and Couture, 2002).
Tectonic map of the Arabia–Eurasia collision zone (modified from Sosson et al., 2010).
The present paper aims to describe, for the first time, the main active slope deformations affecting a mountain ridge that runs along the eastern side of the Enguri water reservoir. The approach is multidisciplinary since it comprises geological, geomorphological and structural surveys, completed with local seismic peak ground acceleration calculation and analyses of slope instability. The usefulness of similar multidisciplinary approaches has also been demonstrated, for example, by Sajinkumar et al. (2017). The focus of our work is to give the main piece of initial evidence that allows the boundaries of the unstable slopes, their geometries, kinematics, slip rate and volume to be defined. Field observations are completed by preliminary 2-D numerical static, pseudostatic and dynamic analyses. 3-D modeling requires further extensive work that will be carried out in the future and deserves proper publication. The results show the presence of wide active slope deformations that affect the two opposite slopes of this mountain ridge.
The paper has a series of main international impacts: (i) it describes an outstanding example of intense gravity instability with double vergence, (ii) it contributes to answers to the challenging scientific questions on the causes of the increase in the instability of a slope facing a reservoir and (iii) it elucidates the regional processes that can contribute to the induction of intense slope instability. The paper also (iv) represents a contribution to assess the hydrogeologic hazard of the largest hydroelectrical facility of Georgia.
Structural section passing through the study area (box), obtained by integrating our field surveys with geophysical and geological data from Banks et al. (1997). Note the presence of a main basal thrust dipping north, which shows a ramp–flat geometry that produced a frontal asymmetric ramp anticline.
The Greater and Lesser Caucasus are two fold-and-thrust belts separated by the Rioni Basin in the Black Sea and the South Caspian–Kura Intermontane Depression (Fig. 1) (Adamia et al., 1977, 2010; Banks et al., 1997; Mosar et al., 2010; Sosson et al., 2010). The Rioni and Kura regions first developed as foreland basins, in Oligocene–early Miocene times, and were successively affected by the fold-and-thrust belts, representing an example of intra-plate, fast-growing mountain building processes in Neogene times (Forte et al., 2010; Adamia et al., 2011; Alania et al., 2016). Fault-bend folds and fault-propagation folds are widespread, with evidence of thin-skinned tectonics in both the Rioni and Kura fold-and-thrust belts (Adamia et al., 2010; Forte et al., 2010; Alania et al., 2016). These belts were subject to intense deformation during the last 14–15 Ma, characterized by an important fault slip along the main thrust systems that reached its greatest exhumation rate at the end of the Miocene (Adamia et al., 2017; Alania et al., 2016). These structural data are coherent with the results of apatite fission-track studies that revealed that the greatest exhumation rate occurred in the Miocene–Pliocene in the Greater Caucasus, although it started here in the Oligocene (Avdeev and Niemi, 2011; Vincent et al., 2007, 2011).
Successively, the mountain building processes continued up to the Quaternary as evidenced by geodetic and seismic data (Tsereteli et al., 2016). More detailed historical and instrumental catalogues, seismic reflection sections, focal mechanism solutions and field geological–structural surveys show the presence of active compressional tectonics, both in the core of the Greater Caucasus and along the southern border and at the fold-and-thrust belts (Tsereteli et al., 2016; Tibaldi et al., 2017a, b). These authors showed that in the core of the Greater Caucasus there is a series of active reverse faults that are parallel to the mountain range (i.e., WNW–ESE). They dip mostly towards NNE and display pure dip slips. Along the southern front of the Greater Caucasus, where the Enguri dam is located, there is a zone of active thrusting along planes dipping towards NNE, as can be seen in the structural section portrayed in Fig. 2. This section has been obtained by integrating geological–structural field surveys from Tibaldi et al. (2017a) with geophysical and geological data from Banks et al. (1997). The section shows the presence of a main basal thrust dipping north. The basal thrust has a ramp–flat geometry that produced a frontal asymmetric ramp anticline corresponding to the area studied here.
The Transcaucasus depression, located between the Greater and Lesser
Caucasus, shows active inversion tectonics in part of the Rioni Basin, with
uplift, folding and faults mostly hidden under the youngest sedimentary
cover (Tibaldi et al., 2017a, b). The largest instrumental and historical
earthquakes, with
From a stratigraphic and lithological point of view, pre-Mesozoic basement and Jurassic sedimentary rocks characterize the axial zone of the Greater Caucasus, whereas Cretaceous and Cenozoic sedimentary rocks are present in the more external zones (Adamia et al., 2011; Mosar et al., 2010). At the foot of the southwestern Greater Caucasus, where the Enguri dam is located, carbonatic deposits dominate together with rarer terrigenous and tuffaceous rocks of volcanic origin. More detailed data on the local geology are given below.
The studied area shows the presence of Jurassic volcanic and terrigenous
rocks and Cretaceous carbonatic deposits generally dipping to the south
(Fig. 4). Cretaceous strata crop out around the Enguri dam where they strike
around E–W and dip mostly at 60–70
Geological map of the area surrounding the western landslide, based on our new field surveys in the onshore zone and data from Zolotarev et al. (1968) for the zone now submerged by the Enguri water reservoir.
Geomorphological map of the area surrounding the western landslide, based on our new field surveys in the onshore zone and data from Zolotarev et al. (1968) for the zone now submerged by the Enguri water reservoir. The main head scarp of the other landslide facing east is also represented.
From a geomorphological point of view, the studied area is characterized by
the presence of two parallel river valleys trending about N–S (Fig. 3a). The
western valley hosts the Enguri dam, whereas the eastern valley is deeply
entrenched, with differences of altitudes between the valley bottom and the
mountain ridge of up to 700 m. These two valleys are separated by a mountain ridge about
9 km long, whose western slope constitutes the side of the
artificial water reservoir. The two opposite slopes of the mountain ridge
are very steep, locally reaching up to an inclination of 70
This road is also affected by deformations related to the presence of the active head scarps of the two landslides that characterize the opposite slopes of the mountain ridge. The main head scarp of the landslide facing west, towards the Enguri reservoir, is represented in the upper part of Fig. 5. The head scarp of the other landslide is represented in the same figure on the slope facing east.
The western slope of the mountain ridge faces the Enguri water
reservoir directly. At a distance of about 2.5 km from the Enguri dam, this slope
shows a series of landforms typical of recent and active deformation: along the
Jvari–Khaishi–Mestia road, at an altitude of 720–740 m, there is a series of
scarps facing westward (Figs. 3c and 5). They have a sinuous shape in plan
view, given by anastomosed single scarps with a westward concave side. This
suggests that these features resulted from the joining of a series of discrete head
scarps. Each scarp is from a minimum of 20 up to 70 m high. At the scarp
foot, the slope is from subhorizontal to gentle dipping westward (average
inclination of 17
In regard to river streams, they have been outlined based on the present network and topographic maps surveyed before the formation of the water reservoir. Rivers mostly run along the average slope dip down to about half of the original slope, with a dendritic pattern (Fig. 5). At the lower half of the slope, now completely covered by the lake, one single river was draining the landslide area. A few river anomalies are present: at the foot of the slope, now under the lake, the lowermost segment of the aforementioned single river was running parallel to the main Enguri River but with a northward flow. In the northern, upper part of the slope affected by the landslide, one short river segment runs perpendicular to the average slope dip. These river diversions may correspond to anomalies in the average slope topography.
In regards to fractures and fissures, these affect the whole area close to
the head scarps (Fig. 6a). The asphalted surface of the Jvari–Khaishi–Mestia
road is affected by several offsets: deformation is represented by fissures,
up to a few centimeters wide, and scarps facing westwards (Fig. 6b and c). These
structures are parallel to subparallel to the morphological head scarps.
The scarps offsetting the road were repeatedly measured from November 2015
to May 2017, showing that the vertical component of the offset
developed up to 14 cm during this time frame, giving an average slip rate of 9.3 cm yr
In the southern segment of the head scarps, a detailed survey in the forest showed the presence of fissures, tens of meters long and up to 3.8 m wide (e.g., Fig. 6d). All the fissures are vegetated and mostly filled, with the local presence of trees grown inside the fissure, with trunks of about 20 cm in diameter. The bottom of the fissures is filled by coarse debris and soil, while the empty upper part of the fissures has a depth in the range of 1–3 m. This suggests that these fissures have a long history, at least of some decades. The area at the foot of this fissure swarm is locally known as the Khoko landslide. Presently, gypsum is excavated here at a small open mine for economic reasons (Fig. 6a). At the foot of this landslide, near the coast of the artificial water reservoir, we recognized the presence of intensely deformed gypsum rocks (Fig. 6e).
We opened two trenches across some of the main active scarps affecting the
Jvari–Khaishi–Mestia road (Fig. 7). The trenches show evidence of repeated
downthrow of the western block. The western part of Trench 1, in fact,
contains two buried old road surfaces: the oldest is lowered by 1.3 m
with respect to the present road level, whereas the second one is lowered by 20 cm
(Fig. 7c). Interviews with local inhabitants indicated that the oldest
road was built during the Soviet era, about 1960 CE
The total subaerial area affected by slope deformation is 1.2 km
The east-facing slope of the same mountain ridge described above is also affected by active deformations (Figs. 9 and 10). Here, the Jvari–Khaishi–Mestia road, at an altitude of 715–720 m, is characterized by a series of fractures with two main orientations: one set is parallel to the road; i.e., it follows the topographic contour lines with a N–S strike. The other fracture set strikes about E–W, perpendicular to the slope.
Most N–S-striking fractures show evidence of extension, with horizontal and
vertical offsets of several centimeters of the road surface. The N–S fractures follow
an escarpment that runs along part of this east-facing slope for a total
length of about 450 m. The height of this scarp is on the order of 3–20 m.
At the southern end of the escarpment, the road has deformations with lower
amounts of offset. At the foot of the head scarp there is a lake, measuring
43 m
Features observed of the eastern landslide, located in the box shown in
The fractures striking E–W show dominant left-lateral strike-slip motions on the order of several centimeters of the road surface (e.g., Fig. 9b). A secondary,
local extensional (vertical) component is also present. The resulting net
slip indicates left-lateral transtensional motions. An old, N–S-trending
concrete wall, located along the head scarp, is tilted 5
The total area affected by slope deformation is 0.37 km
Geo-mechanical properties of the deposits. The final material properties used for the modeling are presented in bold.
To assess the stability of the slope, the depth of the possible slip surfaces and the contribution of the water table and its variation, comprising different lake levels, and to gauge the impact of possible seismic shaking, we developed a series of analyses of the western landslide; we used the program Rocscience 7.024, with the approaches of Janbu and Morgenstern–Price, for static conditions, whereas for seismic stability evaluation, we employed pseudostatic analysis and Newmark modeling.
The slope stability analyses have been carried out based on a geological–geotechnical model that takes data from our field surveys, geotechnical tests and logs drilled during the Soviet era (1966) and by CGS (2015) into account. The logs of 1966 are located in the lower slope, between 343 and 546 m a.s.l., in part corresponding to the area now under the lake level; the more recent logs are located between 566 and 727 m a.s.l. (Fig. 8b). The reconstruction of the stratigraphy of the landslide zone has been done by merging our field surveys with data from the logs no. 3291 (done in 1966), located at 403.62 m a.s.l., no. 3297 (1966) at 418.87 m a.s.l., no. 3261 (1966) at 435.62 m a.s.l. and no. BH1, BH3, BH4, BH5 and BH6 (done in 2015) located at 570–693.7 m a.s.l., and surveys conducted by Soviet researchers previous to the dam construction.
The present outcropping substrate consists of Jurassic volcano-sedimentary rocks, with interlayering of tuff deposits and sandstones, mostly dipping E to NE (Fig. 8). The same rocks and bedding attitude were surveyed at the foot of the slope before the dam construction. These rocks are locally covered by Quaternary landslide and slope debris deposits. Based on their geotechnical characteristics, the succession has been divided into two geotechnical members: the deepest unit (in brown in Fig. 11) is given by the more intact tuffs and sandstones located below the surface located at a depth between 30 and 80 m. Above this surface, we considered a second geotechnical unit (in green in Fig. 11) made of clay, gypsum, debris and old landslide deposits. Debris contains dominant sandstones and tuff fragmented rocks, with a loam matrix, more rarely a clay matrix, rich in carbonatic and gypsum fragments. Locally, carbonatic fragments dominate. The thickness of this unit is from a few meters to a maximum of 80 m. These two units are covered by late Quaternary deposits that comprise colluvium, eluvium and shallow debris, with a thickness from 1 to 5 m. These shallow deposits have been incorporated in the second geotechnical unit (green unit), due to their low thickness and spatial variability. The geomechanical–geotechnical parameters obtained for the materials affected are presented in Table 1. Considering the high variability of deposits and thus of material properties, we performed a back analysis of the landslide to determine the appropriate material properties.
Location of measured piezometers and water table depth; b.g.s. refers to depths below ground surface.
Sections of the static analysis developed to assess the effect of
variation of the water level in the main landslide body.
Figure 6a locates the piezometric wells that have been used to locate the depth of the water table. Our survey of piezometers showed the presence of water that reached the topographic surface during May 2017 (Table 2). In one well (BH9), the water table showed a temporary artesian behavior. The variation of the Enguri reservoir water level was established based on the information received from the Enguri Dam Company, which indicated the range 430–510 m a.s.l. Taking into account the low hydraulic conductivity of part of the deposits composing the studied slope, the large variations of lake level and the intense rains that characterize this region, we focused on slope stability analyses that consider the stage of saturation of the deposits affected and did not carry out a transient slope stability analysis. The final material properties used for the modeling were derived based on back analysis and field observations and are presented in bold in Table 1. In particular, the values of cohesion used take into consideration the widespread presence of clastic debris deposits locally with abundant silty clays.
In the first static model, the reservoir water level is taken fixed at the
maximum level (510 m) (Fig. 11a). The table water level inside the landslide
body is assumed to be characterized by complete saturation as indicated by piezometers during
the rainy season. The model indicates that upper parts of the slope are
unstable along curvilinear slip surfaces (orange lines) that correspond to
the factor of safety (FS)
Seismic acceleration record used for the pseudostatic and Newmark slope stability analysis. See text for details.
In the second static model, the reservoir water level has been considered at
a minimum (430 m a.s.l.), resulting in more diffuse instability of the
slope, including the part normally below the lake. The relative slip surfaces
(orange in Fig. 11b) correspond to FS
Calculated values of PSH.
The pseudostatic numerical model was developed using a peak ground
acceleration (PGA) that has been selected according to the procedure briefly
outlined here below. Since local strong motion records for Georgia do not
exist, we selected the seismograph record of the northwestern Balkan
Peninsula earthquake (epicentral distance 19.7 km; depth
3.75 km;
Examples of some of the Vp and Vs30 velocity profiles measured along the slope affected by the western landslide. All sites of MASW and seismic noise measurement are shown in the inset.
In a further step, we prepared a digital elevation model (DEM) of the slope
using topographic maps at a 1 : 2500 scale. The DEM was converted to a
slope map and the correlation between the slope and Vs30 30 was adopted for
the whole area of the landslide. Finally, taking into account the inputs of
source models (the active fault source model and the area source models) and
the ground motion prediction equation models, we performed a probabilistic
seismic hazard (PSH) assessment in order to individuate the characteristic
earthquakes for PGA calculation. PSH was calculated for each point taking
into account the site characteristic Vs30 for a probability of exceedance of
10 %, 5 %, 2 % and 1 % in 50 years (Table 3). We used the
EZFRISK software for earthquake ground motion estimation. Results show that
the expected possible scenario is represented by a PGA of 0.46 g for site
class B. Finally, the PGA was multiplied by a reference acceleration
coefficient of 0.17 (Kavazanjian et al., 1997) to obtain a seismic load
coefficient of 0.07 (kh
The pseudostatic analysis is presented in Fig. 14a, with the orange and red
slip surfaces corresponding to FS
Results from seismic slope stability analysis.
Due to the evidence of active slope deformation, we installed two digital
extensometers (Wire Linear Potentiometric Transducer, SF500) that record the deformation across the upper scarp of the western landslide continuously.
One extensometer was installed on 1 November 2016, and the other one a few
months later. In Fig. 15 we show the data collected at the older instrument
(location at Trench 2; Fig. 7a), together with data regarding lake level
variation, amount of rain and internal and external (wire) temperatures,
from 1 November 2016 to 15 May 2018. The total amount of deformation here is
67 mm in 18.5 months, which gives an average extension rate of
4.3 cm yr
Data acquired from one of the two digital extensometers located at the upper scarp of the western landslide (location Trench 2; Fig. 7a). The blue line shows the lake level, the red line the extension, the yellow line the external (wire) temperature, the green line the internal temperature and the light blue lines the amount of rain.
In the study area, we found much evidence of active slope deformation that
affects the same mountain ridge on both the western and eastern flanks.
Although part of the western landslide was already locally known as Khoko,
this is the first time that its characteristics are reported in the
scientific literature, whereas the eastern landslide has never been
recognized before. For example, in the Geoportal of Natural Hazards and Risks
in Georgia (
Regarding the western landslide, the data presented here allowed us to define its real total dimension, suggesting that the unstable slope area is larger than the previously known Khoko landslide. We recognized the presence of a main head scarp that runs upward in respect to the Jvari–Khaishi–Mestia road (i.e., east of the road; Fig. 6a) in the northern-central part, then crosses the road and runs west of it. This scarp is made up of three main segments, with the concave side (in plan view) facing west. The segments are interconnected, and thus the total head scarp is continuous, suggesting that the process of slope instability has already led to the linking together of what originally might have been a series of discrete landslides. Historical and present-day activity along this landslide is manifested by the presence of fissures and slip planes at the foot of the head scarp, subsidiary scarps in the slope and systematic tilted trees all along the slope. Where the Jvari–Khaishi–Mestia road crosses the active structures, it requires continuous maintenance in order to fill the opening fractures and smooth the scarps created along the slip planes.
The eastern landslide shows the main head scarp that runs upward respect to the Jvari–Khaishi–Mestia road (i.e., west of the road; Fig. 10). It is made of a single structure with the concave side (in plan view) facing east. Historical and present-day activity along this landslide is manifested by the presence of fissures along the head scarp, some of which are located in the forest and some of which cross the Jvari–Khaishi–Mestia road that also requires continuous maintenance. Moreover, the 1.9 m left-lateral strike-slip offset of the old concrete wall is compatible with landslide movement along the lateral slip plane.
We found evidence that the studied western slope has a long history of landslides. The geological map surveyed during Soviet times, before the infilling of the water reservoir, shows the presence of prehistoric landslide deposits, with a thickness of up to 35–40 m, at the lower part of the slope, now submerged (see Fig. 4 and section in Fig. 8b). Moreover, at the foot of the present onshore slope, there are gypsum deposits that are extremely deformed; this might result from a long history of tectonic movements or from localized landslide movement, or from a combination of both. Whatever the origin, the presence of gypsum deposits favors slope instability. There is also much evidence suggesting that both the western and eastern landslides started their movement at least several decades ago. The western landslide was already moving during the construction of the pavement of the old Jvari–Khaishi–Mestia road. This road was improved during the construction of the Enguri dam, in the Soviet era, by laying a very thick concrete flooring, as showed in Trench 1, which we opened across one of the active slip planes located at the foot of the head scarp (Fig. 7c). The thickness of this concrete layer of 0.5 m and its presence suggest that it was required possibly due to the recognition at that time of active deformation along this road segment. In fact, during the Soviet era, the general thickness of concrete roads was much less. Moreover, the lowered concrete pavement (1.3 m) with respect to the present road level confirms an ongoing process of subsidence here. Then it was followed by the construction of a further road, found 20 cm below the present road level in Trench 1 and 2, which also confirms continuous downsagging of the area. Here, the present-day road surface suffers continuous offsetting, with a consequent necessity to accumulate further tens of centimeters of asphalt to maintain the road level.
The eastern landslide also has evidence of an old history; the small lake at the foot of the head scarp is interpreted here as being a result of upward tilting of the uppermost part of the landslide body. The thick clay deposits associated with this lake suggest that it is relatively old. The dimension of the trees grown within the fissures located uphill of the head scarp also suggests that these fissures opened at least some decades ago. Moreover, several man-made structures show incremental offset with age (the concrete wall of the 1970s, the water channel, the white strips of the asphalt, etc.), confirming the gradual continuous movements of this landslide.
The height of the head scarp of the western landslide (20–70 m) indicates that there was an important vertical slip in the head zone. The presence of a wide subhorizontal to gentle slope at the scarp foot, which has been exploited for the construction of the Jvari–Khaishi–Mestia road, associated with the presence of the head scarp, suggests the development of rotational movements of the upper slope. Ongoing development of vertical scarps in the asphalt of the road and the presence of gentle scarps in the middle slope indicate that rotation movements are still present in different zones of the landslide area. The orientation of the active fissures and slip planes surveyed on the road, which are parallel to subparallel to the main head scarp, and the presence of dip slips along the vertical scarps, indicate a general movement of this western landslide towards the west. In detail, the movement is more complex since the northern part shows motions towards SSW, the central part toward W and the southern part toward NW. These vectors, together with the general shape of the landslide in plan view, suggest that there is a concentric movement of the landslide body towards the toe zone, where bulging and decrease of motion are expected. Since the tow zone is under the Enguri lake, it is impossible to verify what occurs in the lowermost part of the landslide.
In regard to the eastern landslide, in the field there is clear indication of several meters of vertical slip along the head scarp, at the foot of which there is a flat area that contains the small lake. A portion of this subhorizontal area has also been exploited here to host the Jvari–Khaishi–Mestia road. These data indicate the presence of rotational movements of the upper slope. The northern lateral boundary of the landslide is characterized by a prolongation of the main head scarp that crosses the road. At this point, the concrete wall built in 1974 CE shows tilting northward, compatible with the downthrow of the southern block. The same wall also shows the 1.9 left-lateral strike-slip offset that is compatible with the general eastward motion of the landslide body. The northern boundary of the landslide is thus characterized by transtensional movements that also comply with rotational deformation of the upper part of the unstable slope.
For the western landslide, Trench 1 allows us to estimate a long-term (57 years)
slip rate of the vertical component of motion on the order of
2.3 cm yr
For the eastern landslide, the 190 cm offset of the concrete wall
(November 2015) and its age which dates back to 1974 CE result in a slip
rate of 4.6 cm yr
As regards the dimension and type of the western landslide, for which more
data are available, it is necessary to stress that this slope is
characterized by the presence of a thick shallow succession of debris and
clay–gypsum deposits on the order of at least 30–80 m. This indicates the
presence of a thick sequence of deposits that have poor mechanical properties
and that are widespread all over the investigated slope. In fact, both static
and dynamic analyses indicate the possible development of several slip planes
(characterized by FS
The similar position and features of the western and eastern landslides
cannot be a mere series of coincidences, and we cannot exclude the
possibility that they have a common origin. The toes of the two unstable
slopes are located at 400–300 m a.s.l., which corresponds to the level of
erosion of the present river located at the foot of the eastern landslide,
and to the original river level at the base of the western landslide, now
hidden below the Enguri water reservoir. The two river valleys have a
The above-mentioned valleys have been deeply incised by river activity,
consistent with the recent mountain uplift of the region. Thus, the
formation of the two river valleys and the consequent development of the
mountain ridge that gave rise to the two landslides may also be linked with
tectonic uplift. An increase in Plio-Pleistocene exhumation rates in
northwest Georgia is supported by thermochronometric studies (Vincent et al.,
2011) that show that the southern flank of the western Greater Caucasus range has
undergone rapid exhumation of
A gross estimation of the possible age of the studied landslides can be done
if we assume a constant deformation rate; for the western landslide, with the
average height of the head scarp (45 m) and average slip rate
(3 cm yr
Although this is a preliminary description of the evidence and
characteristics of the landslides located near the Enguri dam, the data
presented here are enough to pinpoint a series of relevant results on the
hydrogeological hazard of the area. First of all, the western landslide faces
the Enguri reservoir directly, with the toe zone of the unstable slope
located below the water level, as shown by the pseudostatic and dynamic
analyses that are justified by the fact that the studied area is located in
an active tectonic region where earthquakes with
Our preliminary static analysis suggests that slope movements can occur both when the lake is at the highest level and when it is lowered. This result does not necessarily suggest a possible catastrophic failure of the slope in static conditions, but is instead coherent with the field evidence of continuous deformation along the slope. We interpret this as a consequence of the characteristics of the deposits along the slope. Slip surfaces can activate especially in the part of the slope above the lake level because here deposits rich in clays, and very altered tuffs, dominate. When the water table lowers from 510 to 430 m, although a rapid drawdown condition is created, the shallow deposits of the upper slope remain saturated due to their low hydraulic conductivity. Moreover, an increase of movements can occur during phases of saturation of the onshore landslide body that may take place particularly during the rainy season. At the lowest lake level, slip surfaces can activate in the lower part of the slope due to debuttressing of the slope toe, although deposits with a relatively high hydraulic conductivity dominate here.
The two studied landslides also pose a threat to other infrastructures like the Jvari–Khaishi–Mestia road, which in fact requires continuous maintenance to avoid the gradual formation of dangerous steps in the asphalt that are tens of decimeters high. We also highlight the presence of an inhabited house resting in the middle of the western landslide, and another house recently built close to the northern boundary slip zone of the eastern landslide.
In order to contribute to a better assessment of the hydrogeological hazard of the area, we installed two digital extensometers in correspondence with the active head structures of the western landslide, together with the benchmarks for GPS measurements that have already been mentioned. This monitoring effort will help to individuate possible alert situations and to better constrain the behavior of this unstable slope.
Here we focus on the western landslide, its relations with the artificial lake and similar settings reported by literature. As shown by the Quaternary geological deposits and by the presence of the high head scarp, the studied area was already affected by the emplacement of landslide deposits in prehistoric times. The already destabilized slope was thus interested by the formation of the Enguri reservoir and its level variations. The creation of artificial lakes can have different effects: on one side, it triggers a possible seepage process, resulting in an increase in porewater pressure within the slope deposits, thus reducing their shear strength. On the other side, the presence of the lake water body induces a stabilizing load at the toe of the submerged part of the slope (Paronuzzi et al., 2013). In transient conditions, the rate of lake filling or drawdown, combined with the permeability of the bank-forming materials, produce different effects of reservoir level changes. Several landslides, in fact, have been triggered by filling–drawdown operations (Schuster, 1979; Kenney, 1992; Zhu et al., 2011); pre-existing, ancient landslides have been reactivated during water reservoir formation as well (e.g., Kaczmarek et al., 2015).
In the case of low permeability of slope materials, rapid drawdown of the lake level can produce a decrease in the factor of safety, possibly leading to slope failure (Kenney, 1992). A rapid decrease in lake level, in fact, results in a short-term increase in hydraulic gradient in adjacent slopes (Jones et al., 1961; Schuster and Wieczorek, 2002; Deying et al., 2010; Pinyol et al., 2012). For example, in Japan about 60 % of reservoir landslides occurred during a sudden decrease of the water level (Nakamura, 1990). In a DSGSD located above the Gepatsch Dam reservoir (Austrian Alps), slope deformation rates correlate with reservoir levels and drawdown conditions (Zangerl et al., 2010). In the case, instead, of high permeability of bank-forming materials, the reservoir level increase can produce a decrease in the factor of safety. For example, for the Vajont (Italy) catastrophic landslide, detailed numerical analyses by Paronuzzi et al. (2013) show the predominant role played by reservoir levels in determining slope instability, with respect to the rain-induced water table in the upper slope. In particular, the initial large slope deformation took place in concomitance with the lake level rise. Consistent with this trend, the final collapse of 9 October 1963 coincided with the maximum reservoir level being reached. Paronuzzi et al. (2013) also show that the decisive geological factor contributing to the Vajont collapse was the presence of a landslide that already existed: the prehistoric rockslide was characterized by materials with high permeability and a thick shear zone at the base, including montmorillonitic clay lenses. The combination of poor mechanical properties of clay beds with the high permeability of the angular limestone gravel determined a rapid reservoir-induced inflow that reduced strength and the factor of safety. A further example is given by the Byford Creek slide, located above the Clyde Dam reservoir (New Zealand). Here, long-term creep movements show a clear reaction to lake filling, with an initial large increase in the deformation rate and long-term slowing of movement (Macfarlane, 2009).
The landslide studied here, more than 40 years since the construction of the Enguri reservoir, still shows a high sensitivity to water infilling operations. The unstable slope behavior, during large water level variations, is similar to other landslides with bank-slope material characterized by high permeability. In particular, the deposits found in the logs of the lower part of the slope (below about 500 m a.s.l.), being characterized by the dominance of angular carbonatic clasts with interspersed clays, are similar to those present at other landslides, for which the factor of safety decreases with increase of lake level. With such characteristics of the materials affected, during reservoir increase the porewater pressure effect on the shear strength prevails over the stabilizing buttressing of the lake water body, inducing an acceleration in slope movements. This is evidenced by the acceleration in slope deformation following the Enguri reservoir infilling (Fig. 15).
We presented, for the first time, the evidence and characteristics of two large landslides located near the Enguri dam and water reservoir, in the southwestern part of the Caucasus (Republic of Georgia).
The landslides affect the two opposite slopes of a mountain ridge that runs
parallel to the Enguri reservoir. The slope directly facing the reservoir
shows active deformations that involve a subaerial area of 1.2 km
The head scarp zones of both landslides affect the main
Jvari–Khaishi–Mestia road with offsets of man-made features that indicate
short-term (last 2 years) slip rates of up to 9.3 cm yr
The possible causes for the past inception of the two landslides are multiple. First of all, we propose the increase above a threshold value of
the slope dip and of the topographic altitude difference created by river
erosion. The toe zones of both landslides, in fact, are located in correspondence with deeply entrenched rivers; the river's location at the foot
of the western landslide is now below the Enguri water reservoir. The
excavation of the two parallel river valleys, in turn, can be linked with the
increase in Plio-Pleistocene exhumation rates of this mountain area of
Extensometer measurements across the head scarp of the western landslide
indicate a present-day variable deformation (average 4.3 cm yr
The data are available from the authors upon request.
AT coordinated the work, contributed to the collection of the field data and wrote most of the manuscript. JG and TO carried out the pseudostatic and dynamic analysis of slope stability, PO carried out the static analysis and they helped to write the related chapters. PO also collected the extensometer data. NT and ZG calculated PSH and selected PGA, carried out the Vp and Vs30 velocity profiles and helped to write the related text. The discussion was agreed by all authors.
The authors declare that they have no conflict of interest.
We acknowledge two anonymous referees for their useful suggestions on a previous version of the paper. This study has been conducted in the framework of the NATO project SfP G4934 “Georgia Hydropower Security”, of the International Lithosphere Program – Task Force II, of the European Space Agency project no. 32309 “Active tectonics and seismic hazard of southwest Caucasus by remotely sensed and seismological data” (Leader Alessandro Tibaldi) and of project 216758 of the Shota Rustaveli National Science Foundation. Edited by: Paola Reichenbach Reviewed by: two anonymous referees