GEOLOGIJA 50/1, 121–129, Ljubljana 2007 Rheological investigation for the landslide Slano Blato near Ajdov{~ina (Slovenia) Reolo{ke raziskave za plaz Slano Blato pri Ajdov{~ini Karmen FIFER BIZJAK1 & Andreja ZUPAN^I^2 1Slovenian National Building and Civil Engineering Institute, Dimi~eva 12, 1000 Ljubljana, Slovenia 2Faculty of Chemistry and Chemical Technology, Department of Chemical Engineering, University of Ljubljana, A{ker~eva 5, 1001, Ljubljana, Slovenia Klju~ne besede: zemljinski tok, plaz, geotehnika, reolo{ke lastnosti materiala, numeri~no modeliranje, Burgerjev model, Slano Blato, Slovenija Key words: earth flow, landslide, geotechnical engineering, rheological properties, numerical modelling, Burger model, Slano Blato, Slovenia Abstract The landslide Slano Blato, is situated above the village Lokavec near Ajdov{~ina in the weat of Slovenia. It has a relatively long history and was first mentioned in a document in 1887. At that time it destroyed a part of a main road and reconstruction works took 17 years. In the last decade, movement of the landslide was observed in November 2000, when it reached distances of 60–100 m/day. By means of geotechnical research work on the landslide in the year 2004, several rheological tests were also carried out, which is not usual for geotechnical research work. A stability analysis was carried out numerically by applying the Burger elasto-plastic model. The model took into account geomechanical and rheologi-cal characteristics of the landslide. Izvle~ek Plaz Slano Blato leži nad vasjo Lokavec pri Ajdov{~ini v zahodnem delu Slovenije. Ima že dokaj dolgo zgodovino, saj je bil v dokumentih omenjen že leta 1887. V tem ~asu je uni~il glavno cesto in njegova sanacija je trajala 17 let. V zadnjem desetletju so bili prvi ve~ji premiki na plazu v novembru 2000. Najve~ji izmerjeni pomiki splazele mase so bili 60–100 m/dan. V letu 2004 so se izvajale geotehni~ne raziskave za namen pridobitve podatkov za sanacijo plazu. V sklopu teh raziskav so bile izvedene tudi obsežne reolo{ke raziskave. Stabilnostne analize so bile izvedene z Burgerjevim elasto-plasti~nim modelom. Model pri izra~unu upo{teva geomehanske in reolo{ke karakteristike splazele mase. Introduction Landslide is an important erosion process in Slovenia, affecting 8000 km2 of labile or potentially unstable slopes (around 40 % of the country’s area) mainly composed of unconsolidated or partially consolidated fine – grained soils (Miko{ et al., 2004). The Slano Blato landslide is situated in the west of Slovenia (Figure 1), above Lokavec village (Figure 2), on the border between the Alps and the Mediterranean region. Within 100 km of Slano Blato there are three large landslides, each with a potential sliding mass of over 500,000 m3. Also some other relatively large landslides are located in the Slovenian Alps. 122 Karmen Fifer Bizjak & Andreja Zupan~i~ The Slano Blato landslide has a relatively long history and having been first mentioned in a document in 1887. At that time it destroyed a part of a main road after which reconstruction works took 17 years. More recently, movement of the landslide was observed in November 2000, when it reached peak displacement rates of 60–100 m/day. The landslide Slano Blato was probably activated, as an earth flow, by a combination of several events. In 1998, a strong earthquake occurred in the Upper So~a River valley (Vidrih, 2001). And next, the year 2000 was very wet and the old drainage system was not maintained any more. Intense rainfalls which cause large landslides are not rare in the Alpine regions of Italy, (Guzzeti, 2004), Switzerland and Austria (Moser, 2002). Detailed investigation of the landslide was undertaken in the year 2004, to determine the depth of the landslide, its geotech-nical and rheological parameters. The models, which are usually used for the analysis of landslides are the Mohr-Coulomb, the Drucker Prager or the von Mises model. A numerical stability analysis was performed by applying the finite difference method with Burger visco-plastic model. Figure 2. The Landslide Slano Blato, above the village Lokavec Slika 2. Plaz Slano Blato, nad vasjo Lokavec Rheological investigation for the landslide Slano Blato near Ajdov{~ina (Slovenia) 123 The model takes into account the geome-chanical and rheological characteristics of the landslide. The model has not been used in geotechnical practice very often. Probably is the reason in the relevant rheological soil parameters for the calculation. For that aim exacting rheological tests have to be done that are not usual tests in the geotechnical investigation work. The finite element model has already been used for the Val Pola landslide, using the classical elasto-plastic law and quasi static time stepping (Crosta, 2001). The Burger visco-plastic model, which was used in our case, has not been used for the landslide problems yet. Geology The Slano Blato landslide is situated at the contact between Triassic limestone and Eocene flysch formation (Figure 3). The Eocene flysch consists of marl and layers of sandstone with thickness of centimetres or several meters. The rock is highly tectoni-cally deformed. The limestone was overthrusted on to the flysch over a very large distance along the Trnovski overthrust. In consequence the region consists of large synclines and anticlines. The upper part of the limestone is at 670 m above sea level. Limestone is fissured into blocks with dimensions of several cubic meters. The limestone is overthrusted on to the massive sandstone, which belongs to the flysch series. The contact dips at approximately 10° to the NW. Based on the sediment texture, inverse positions of layers were established in the upper part of the landslide. Inverse layers are observed even up to 550 m. Throughout the landslide several faults with dinaritic dip direction 330–345/ 55–75° were observed. The flysch in the region of the landslide can be divided into three parts: – layers of sandstone of thickness between one and several meters – a region with alternation of marl and 10 cm thick layers of sandstone – a region where layers of marl prevail under the layers of sandstone. The dip direction of layers is WN, which is favourable for slope stability. The flysch is covered with clayey gravel, that forms the landslide mass. This is very Figure 3. Geological map and profile of the landslide Slano Blato Slika 3. Geolo{ka karta in profil za plaz Slano Blato 124 Karmen Fifer Bizjak & Andreja Zupan~i~ heterogeneous, with blocks of limestone, sandstone and clayey silty gravel. The thickness of the gravel is 3 to 10 m, as ascertained by borehole logging and geophysical investigations. ples that were taken during a rainy period of drilling, on the other hand, had water content of 16–17 % and lower strength. The index of plasticity, Ip, was between 19 and 32 % in keeping with results obtained for Austrian earth flows (Moser, 2002). Field and laboratory measurements Geomechanical characterisation of the landslide Geomechanic characterisation of the gravel material and the flysch have been undertaken as follows. Flysch layers Logging of the boreholes determined that the RQD index of the flysch is 0. The marl is highly tectonized, in places even into silty grey clay. The intermediate layers of sandstone are also highly fissured. In the laboratory tests the water content (w), liquid limit (wl), plastic limit (wp), dry unit weight (?dry), uniaxial strength (Qu), angle of friction (?) and cohesion (c) were determined (table 1). The uniaxial strength Qu proved to be very low. Shear vane test, yielded on intact friction angle of 8.4°, and a residual friction angle of, ?rez = 4.4°. These values are very low. Table 1. Geomechanical properties of the flysch layers Tabela 1.Geomehanske lastnosti fli{nih plasti w % 10-14 y dry Qu KN/m3 kN/m2 c KN/m2 18-20.6 155-323 22-27 0-37 Landslide gravel material The results of geotechnical testing of the landslide gravel material proved to be highly depended upon the condition at the take of sampling. Samples collected in dry weather had a much lower water content (around 12 %) and higher strength (Table 2). The sam- Rheological investigations For rheological tests samples were prepared from sieved landslide materials with particles below 63 ΅m. This particle fraction represented 36 wt % of the total sieved material, which was proved as the critical part for material sliding. Samples of different water contents, from 35 to 60 %, were prepared by adding water to dry powders. Samples were taken from the surface and at a depth of 8 m at different locations of the upper part of the landslide. In order to allow complete wetting of the particles, rheologi-cal tests were performed two days after the preparation. The rheological tests were carried out using the controlled stress rheometer, Haake RheoStress RS 150, equipped with a parallel plates sensor geometry (with serrated surface), PP 25, of 1.2 mm gap. The measurements were carried out under destructive and non destructive shear conditions at 23 °C. A small amount of water, placed around the measuring sensor and covered by a solvent trap, was necessary to prevent water evaporation from the samples. Due to the peculiar behaviour of the investigated semisolid samples, measuring protocols were predetermined in order to obtain repeat-able experimental data. Shear stress was increased for 3 minutes under destructive shear conditions. The shear stress ranges, which depended on the amount of water in the sample, were selected in the range where the transition from creep flow to shear flow was expected. The viscoelastic properties of samples were examined under non-destructive conditions of oscillatory shear and by creep-recovery tests. The upper limit of the linear viscoelastic response (LVR) was determined Table 2. Main geomechanical properties of sliding mass Tabela 2. Geomehanske lastnosti splazele mase Properties w % Boring in dry period Boring in wet period 12 16-25.6 % 21 20.4 wl % 52 50 Ydry KN/m3 19.6-20.3 15.2-15.6 Qu kN/m2 64-125 E MPa 29 24 c KN/m2 15 18 Rheological investigation for the landslide Slano Blato near Ajdov{~ina (Slovenia) 125 w = 50% w = 40% GM – Maxwell, (Pa) 1.92E + 05 1.11E + 06 ?M Maxwell, (Pa.s) 1.95E + 07 1.30E + 09 ?K (Kelvin Voight), (Pa.s) 2.14E + 06 5.02E + 06 GK (Kelvin Voight) (Pa) 5.92E + 04 1.53E + 05 Table 3. The material properties, the elastic modulus GM,K and the viscosities, ?M,K, evaluated as parameters of the Burger model for sample 405 at water content of 40 % and 50 % Tabela 3. Reolo{ke lastnosti, elasti~ni modul GM,K viskoznost, ?M,K, kot parametri Burgerjevega modela za vzorec 405 pri vsebnosti vlage 40 % in 50% from stress sweep tests under oscillatory shear conditions. The mechanical spectra of the examined suspensions were evaluated by applying frequency sweep experiments in the LVR. Creep-recovery tests were performed by applying a constant shear stress for 5 min and measuring the increased shear deformation in the sample during creep and recovery. By increasing the shear stress during the rheological tests under destructive shear conditions the material deformed (creep flow) and in narrow shear stress range the viscosity drastically dropped, as shown in Figure 4. Due to high solid loadings and the nature of the samples, a homogeneous shear flow field was not achieved. Critical shear stress, at which the transition from creep flow to shear flow occurred, was determined for samples with different contents of water. Viscosities of the samples at different water contents were taken from the plateau region of the creep flow. As shown in Figure 5, the critical shear stress and the viscosity strongly decreased with increased water content in the investigated samples, independently of i.oe+08 -i *s US> 1 8 1 1.0E+06 • - fN W-3JW ¦f? ¦» ^ r * 1 l.OE+04- *» ¦ W-50% W-4QM [".!¦:,; . Si LI 1 : J A :unH!r]iii1i]tl]ic bubble 1.0E-HI2 1.0E+« 1.0E444 1.0E44K Slii-Hr fin.»-. (Pill Figure 4. Determination of critical shear stress from stress sweep tests for the sample from middle part of the landslide contained different water content, from 35 % to 60 % Slika 4. Dolo~itev kriti~ne strižne trdnosti iz reolo{ke preiskave za vzorec iz osrednjega dela plazu, pri vsebnosti vlage med 30 in 60 % sample location in the landslide. From rheo-logical tests and from the geomechanical characterisation of the examined landslide materials it was concluded that the critical water content for the formation of earth flow could be about 40 %. Figure 5. Influence of water content on critical shear stress and viscosity of examined samples taken from different regions of the landslide: A: middle part, B: from borehole of the upper part, C: upper part – at the same location as B, D: upper part Slika 5. Vpliv vlage na kriti~no strižno trdnost preiskanih vzorcev vzetih iz A: osrednjega dela plazu, B: iz vrtine zgornjega dela plazu, C: povr{ina, ista lokacija kot B, D: spodnji del plazu 126 In order to determine the parameters required for numerical simulations of the earth flow with the FLAC program, the rheological tests of the samples having water content of 40 % and 50 % were performed under nondestructive shear conditions. Measurements showed that the investigated samples with water content of 40 % exhibited viscoelastic behaviour with predominant elastic contribution to the viscoelastic response. From the experimental data measured in creep and recovery tests, the parameters of the Burger model were evaluated. The Burger model (Barnes, 2000) describes the response of many real viscoelastic materials on the applied constant shear stress (?c) in the range of linear response. It consists of four simple mechanical elements, two springs (G – elastic modulus) and two dash-pots (? – viscosity) and represents a combination of the Maxwell (describes viscoelastic-liquid response) and the Kelvin-Voight (describes viscoelastic-solid response) mechanical model in a series (Figure 6). Figure 6. Burghers mechanical model - a combination of Kelvin-Voight mechanical model and Maxwell model Slika 6. Burghersov model sestavljata zaporedno vezana Kelvin-Voightov in Maxwellow model The Burger model describes shear deformation (y) of the material during creep tests as: y(t)=^+^+^ • (l-exp(-tAK)} t\m gm gk and y(t) = x(t)/G(t) = xJG(t) where XK represents the Kelvin-Voight relaxation time, XK = T|K/GK and xc the applied shear stress and G(t) the time dependence of the shear modulus. The Burger model is often written in terms of time dependence of compliance J(t), which is defined as [J(t) = 1/G(t)]. This means that time dependence of shear deformation can be written as: y(t) = x(t) · J(t) = xc · J(t) Karmen Fifer Bizjak & Andreja Zupan~i~ Then the Burger model can be expressed: J(t)=—+JM+JK -{l-exp(-tAK)} For the evaluation of the model parameters it is necessary to determine the dependence of y(t) experimentally. Viscous (t|m and r|K) and elastic (GM and GK) contributions to the viscoelastic response of the investigated sample, taken from the upper part of the landslide, at water contents of 40 % and 50 %, were calculated from the creep tests. For calculation the examined range of water content was selected in order to compare the model parameters evaluated for the sample with solid-like viscoelastic response (w = 40 %) with the sample with fluid-like response (w = 50 %). As reported in Table 3, the decrease of humidity in the examined range increased the values of material properties by at least ten times. The experimental data and the corresponding curves calculated by using the Burger model are shown in Figure 7. It is evident that the model correlated with the experimental data in a satisfactory way. For calculation of all Burger’s parameters the SOLVER protocol in Excel was used. Experimental data from creep and recovery curve were fitted with Burger’s parameters. Figure 7. Creep and recovery tests for the sample taken from the upper part of the landslide at water content of 40% and 50%. Curves passing the experimental data are correlated by the Burger model Slika 7. Test lezenja za vzorce vzete iz zgornjega dela plazu, z vsebnostjo vode 40 in 50 %. Krivulja je korelirana z Burgerjevim modelom Geomechanical and rheological investigations of landslide materials showed that the material properties of the samples are strongly influenced by water content, as B-B Rheological investigation for the landslide Slano Blato near Ajdov{~ina (Slovenia) 127 well as by the time in which the materials were exposed to wetting, therefore by the precipitation conditions on the landslide. During drilling, influx of underground water was observed in layers at different depths. Numerical analyses The aim of numerical modelling was to establish the most critical parts of the landslide and the mechanism of failure. For the calculation, the FLAC program was used. This is an explicit finite difference program that performs Lagrangian analyses. Because of the large dimension of the landslide, it was divided for the numerical calculations into 5 regions. Each region begins and finishes with stable layers of sandstone. In that way the influence of the border effects was reduced. Burger - visco-plastic model The Burger visco-plastic model considers visco-elasto-plastic deviatoric behaviour and elasto-plastic volumetric behaviour. The visco-elastic and plastic strain-rate component act in a series. The visco-elastic part corresponds to the Burger model and the plastic part to the Mohr-Coulomb model. The deviatoric component was described with the relation: e.. = eK. + e™ + ep. V 9 8 V where eθ - Kelvin strain component, e? -Maxwell strain component, e?- plastic strain component, "J where they are; g - potential function, A, plastic flow parameter, e^ - plastic volume strain rate. The volumetric behaviour is determined with: where e, is volume strain rate. Results of the Burger visco-plastic model The Burger visco-plastic model was used only for the gravel landslide mass. The aim of this calculation was to simulate the earth flow of the landslide gravel mass at the time of the largest movements. For the flysch layers, the Mohr-Coulomb model was used. Input geomechanical and rheological data are given in Table 4. Table 4. Geomechanical parameters for the Burgher visco-plastic model of gravel and slide mass Tabela 4. Geomehanski parametri za Burgerjev visko-plasti~en model za splazelo maso y KN/ /m3 E MPa ° c KN/ /m2 V MPa ?M Pas GK MPa ?K Pas 21 3 24 1 0,3 1,11e6 1,3e9 1,53e5 5,02e6 Considering the material properties in the Burger visco-elasto-plastic model, the maximal deformations developed in one day were calculated with numerical simulations. At the time of maximal landslide movement the observed deformations were in the range of 60–100 m/day. The best results of numerical simulations were obtained when the geomechanical parameters were taken into account for the material water content ranging between 35 and 40 %. With these material properties, taken at different parts of the landslide, the simulated deformations were in the range 70–80 m/day. The results of the numerical analyses are presented in Figure 8, from the upper part of the landslide to the bottom. At the end of the landslide, several meters thick layer of sandstone represents a natural barrier before the village. There was an important question whether that thick layer of sandstone was strong enough or whether it could collapse under additional gravel landslide mass from the upper part of the landslide. Geomechanical investigation proved, with inclinometer measurements, that the natural barrier is stable. Also numerical analyses confirm that it is stable even under additional landslide mass. In the case of landslide moving, the landslide mass will pour over the sandstone layer, but the natural barrier of sandstone will remain stable. 128 Karmen Fifer Bizjak & Andreja Zupan~i~ Figure 8. Results of Burger elasto-plastic mode; displacements per day Slika 8. Rezultat izra~una z Burgerjevim visko-plasti~nim modelom; premiki na dan Rheological investigation for the landslide Slano Blato near Ajdov{~ina (Slovenia) 129 Conclusions With the visco-plastic Burger model the rheological characteristics of gravel landslide mass were considered. Taking into account also the rheological properties of the materials at different landslide locations, the model allowed us to decribe the actual situation on the landslide. In this way, we were able to simulate even the largest deformations when the smallest particles fraction had a water content between 35 and 40 %. The Burger model is not often used for this type of geotechnical problem. Our experience is that together with good rheological tests it could present quite a reliable prediction of landslide movement with time. Geomechanical laboratory tests showed that shear properties depend on the percentage of moisture. The samples taken from boreholes during a wet period had poorer geomechanical properties than those taken during a dry period, but the evaluated results did not differ drastically. Greater differences were observed from rheological characterisations. For example, the viscosity changes were in the range between 103 and 108 Pas, depending on the content of water. Five reinforced concrete wells – dowels were constructed in the upper part of the landslide in year 2005. Wells were used for dewatering and as a retaining structure. Until that larger movements on landslide have not observed. Acknowledgement This work was financialy supported by Ministry of Environment, Spatial planing and Energy. References Alleoti, P. 2004: A warning system for rainfall-induced shallow failures. – Engineering geology, 73, 247–265. D'Amato, A. G., Gianecchini, R. & Puc-cinelli, A. 2004: The influence of the geological and geomorphological settings on shallow landslides. An example in a temperate climate environment: the June 19, 1996 event in northwestern Tuscany (Italy). – Engineering geology, 73, 215–228. Duan, J. 1996: A coupled hydrologic-geomor-phic model for evaluating effects of vegetation change on watersheds. – Ph.D. thesis, Dept. of Forest Engineering, Oregon State University. Barnes, H. A. 2000: A handbook of elementary rheology. – University of Wales, Institute of non Newtonian fluid mechanics, Aberystwth SY23 3BZ. Buma, J. & Dehn, M. 1998: A method for predicting the impact of climate change on slope stability. – Envir. Geol., 35, 190–196. Crosta, G., Calvetti, F. , Imposimato, S., Roddeman, D., Frattini, O. & Agliardi, F. 2001: Granular flows and numerical modelling of landslides. – Damocles, Universita di Milano. Guzzeti, F. , Cardinali, M., Reichenbach, P. , Cipolla, F. , Sebastiani, C., Galli, M. & Salvati, P. 2004. Landslides triggered by 23 November 200 rainfall event in the Imperia province, Western Liguira, Italy. – Engineering geology, 73, 229–245. Hutchinson, J.N. 1988: General Report: morphological and Geotechnical parameters of landslides in relation of geology and hydrogeo-logy. – In. Bonnard, C, (Editor), Proceeding, Fifth International Symposium of Landslides, Christ-church, N.Z., 3, 1805–1841. Iverson, R.M. & Valance, J.M. 2001: New views of granular mass flows. – Geology, 29, 115– 118. Miko{, M., ^etina, M. & Brilly, M. 2004: Hydrologic conditions responsible for triggering the Stože landslide, Slovenia. – Eng. geol., 73, 3/4, 193–213. Ko~evar, M. 2001: Poro~ilo o geolo{kih in hidrogeolo{kih geomehanskih raziskavah na plazu Slano Blato. – Internal report, archives MOP, Ljubljana. Moser, M. 2000: Geotechnical aspects of landslide in the Alps, – In: Rybar, J., Stember, J., Wagner, P. (Eds), Landslides, 23–43. Marinos, P. & Hoek, E. 2000: GSI: A geological friendly tool for rock mass strength estimation. – Proceedings of the International Conference on Geotechnical & Geological Engineering (GeoEng, 2000), Technomic Publishing Co. Inc., Melbourne, Australia, 1422–1440. Petkov{ek, A. 2002: Poro~ilo o geomehanskih laboratorijskih raziskavah zemljine iz plazu Slano Blato nad Lokavcem. – Internal report, archives ZAG, Ljubljana. Sasaki, Y. , Fujii, A. & Asai, K. 2000: Soil creep process and its role in debris slide generation – filed measurements on the north side of Tsukuba Mountain in Japan. – Engineering Geology, 56, 163–183. Vidrih, R., Ribi~i~, M. & Suhadolc, P. 2001: Seismological effects on rocks during the 12 April 1998 upper So~a territory earthquake (NW Slovenia). – Tectonophysics, 330 3/4, 153–175.