Slope stability analysis of the Pacaya Volcano, Guatemala, using limit equilibrium and finite element method

The Pacaya volcanic complex is part of the Central American volcanic arc, which is associated with the subduction of the Cocos tectonic plate under the Caribbean plate. Located 30 km south of Guatemala City, Pacaya is situated on the southern rim of the Amatitlan Caldera. It is the largest post-caldera volcano, and has been one of Central America’s most active volcanoes over the last 500 years. Between 400 and 2000 years B.P, the Pacaya volcano had experienced a huge collapse, which resulted in the formation of horseshoe-shaped scarp that is still visible. In the recent years, several smaller collapses have been associated with the activity of the volcano (in 1961 and 2010) affecting its northwestern flanks, which are likely to be induced by the local and regional stress changes. The similar orientation of dry and volcanic fissures and the distribution of new vents would likely explain the reactivation of the pre-existing stress configuration responsible for the old-collapse. This paper presents the first stability analysis of the Pacaya volcanic flank. The inputs for the geological and geotechnical models were defined based on the stratigraphical, lithological, structural data, and material properties obtained from field survey and lab tests. According to the mechanical characteristics, three lithotechnical units were defined: Lava, Lava-Breccia and Breccia-Lava. The Hoek and Brown’s failure criterion was applied for each lithotechnical unit and the rock mass friction angle, apparent cohesion, and strength and deformation characteristics were computed in a specified stress range. Further, the stability of the volcano was evaluated by two-dimensional analysis performed by Limit Equilibrium (LEM, ROCSCIENCE) and Finite Element Method (FEM, PHASE 2 7.0). The stability analysis mainly focused on the modern Pacaya volcano built inside the collapse amphitheatre of “Old Pacaya”. The volcanic instability was assessed based on the variability of safety factor using deterministic, sensitivity, and probabilistic analysis considering the gravitational instability and the effects of external forces such as magma pressure and seismicity as potential triggering mechanisms of lateral collapse. The preliminary results from the analysis provide two insights: first, the least stable sector is on the south-western flank of the volcano; second, the lowest safety factor value suggests that the edifice is stable under gravity alone, and the external triggering mechanism can represent a likely destabilizing factor.

Corrección del deslizamiento de una ladera en la autovia A-4. PK. 340 en el término municipal de Marmolejo, Jaén mediante pantalla de pilotes y sistemas de drenaje

En este proyecto se describe el cálculo y diseño de una serie de medidas correctoras propuestas a raíz de un deslizamiento de ladera que tuvo lugar en noviembre de 2010 y afectó a parte de la calzada de la autovía A-4 en el P.K. 340 entre Córdoba y Jaén; principalmente en la ejecución de una pantalla de pilotes como elemento estructural de estabilización y contención del terreno, por medio del método de equilibrio límite y cálculo del módulo de reacción o coeficiente de balasto horizontal, a través del empleo de los programas Slope/W y RIDO respectivamente. ABSTRACT This project describes the calculation and design of a series of proposed corrective measures following a landslide of which took place in November 2010 and hit part of the road A-4 on the P.K. 340 between Cordoba and Jaen, mainly in the execution of a pile wall as a structural element of stabilization and containment of the terrain, through the limit equilibrium method and calculating the reaction module or horizontal ballast coefficient, through employment of the Slope/W and RIDO programs respectively.

Inception of debris avalanches: remarks on geomechanical modelling

Debris avalanches are complex phenomena due to the variety of mechanisms that control the failure stage and the avalanche formation. Regarding these issues, in the literature, either field evidence or qualitative interpretations can be found while few experimental laboratory tests and rare examples of geomechanical modelling are available for technical and/or scientific purposes. As a contribution to the topic, the paper firstly highlights as the problem can be analysed referring to a unique mathematical framework from which different modelling approaches can be derived based on limit equilibrium method (LEM), finite element method (FEM), or smooth particle hydrodynamics (SPH). Potentialities and limitations of these approaches are then tested for a large study area where huge debris avalanches affected shallow deposits of pyroclastic soils (Sarno-Quindici, Southern Italy). The numerical results show that LEM as well as uncoupled and coupled stress–strain FEM analyses are able to individuate the major triggering mechanisms. On the other hand, coupled SPH analyses outline the relevance of erosion phenomena, which can modify the kinematic features of debris avalanches in their source areas, i.e. velocity, propagation patterns and later spreading of the unstable mass. As a whole, the obtained results encourage the application of the introduced approaches to further analyse real cases in order to enhance the current capability to forecast the inception of these dangerous phenomena.

Geotechnical Data and Numerical Analysis of Edifice Collapse and Related Hazards at Pacaya Volcano, Guatemala

The continual eruptive activity, occurrence of an ancestral catastrophic collapse, and inherent geologic features of Pacaya volcano (Guatemala) demands an evaluation of potential collapse hazards. This thesis merges techniques in the field and laboratory for a better rock mass characterization of volcanic slopes and slope stability evaluation. New field geological, structural, rock mechanical and geotechnical data on Pacaya is reported and is integrated with laboratory tests to better define the physical-mechanical rock mass properties. Additionally, this data is used in numerical models for the quantitative evaluation of lateral instability of large sector collapses and shallow landslides. Regional tectonics and local structures indicate that the local stress regime is transtensional, with an ENE-WSW sigma 3 stress component. Aligned features trending NNW-SSE can be considered as an expression of this weakness zone that favors magma upwelling to the surface. Numerical modeling suggests that a large-scale collapse could be triggered by reasonable ranges of magma pressure (greater than or equal to 7.7 MPa if constant along a central dyke) and seismic acceleration (greater than or equal to 460 cm/s2), and that a layer of pyroclastic deposits beneath the edifice could have been a factor which controlled the ancestral collapse. Finally, the formation of shear cracks within zones of maximum shear strain could provide conduits for lateral flow, which would account for long lava flows erupted at lower elevations.

Evaluation of slope stability under water and seismic load through the Minimum Lithostatic Deviation method

The main objective of this thesis is to obtain a better understanding of the methods to assess the stability of a slope. We have illustrated the principal variants of the Limit Equilibrium (LE) method found in literature, focalizing our attention on the Minimum Lithostatic Deviation (MLD) method, developed by Prof. Tinti and his collaborators (e.g. Tinti and Manucci, 2006, 2008). We had two main goals: the first was to test the MLD method on some real cases. We have selected the case of the Vajont landslide with the objective to reconstruct the conditions that caused the destabilization of Mount Toc, and two sites in the Norwegian margin, where failures has not occurred recently, with the aim to evaluate the present stability state and to assess under which conditions they might be mobilized. The second goal was to study the stability charts by Taylor and by Michalowski, and to use the MLD method to investigate the correctness and adequacy of this engineering tool.

Boundary integral equations on unbounded rough surfaces: Fredholmness and the finite section method

We consider a class of boundary integral equations that arise in the study of strongly elliptic BVPs in unbounded domains of the form $D = \{(x, z)\in \mathbb{R}^{n+1} : x\in \mathbb{R}^n, z > f(x)\}$ where $f : \mathbb{R}^n \to\mathbb{R}$ is a sufficiently smooth bounded and continuous function. A number of specific problems of this type, for example acoustic scattering problems, problems involving elastic waves, and problems in potential theory, have been reformulated as second kind integral equations $u+Ku = v$ in the space $BC$ of bounded, continuous functions. Having recourse to the so-called limit operator method, we address two questions for the operator $A = I + K$ under consideration, with an emphasis on the function space setting $BC$. Firstly, under which conditions is $A$ a Fredholm operator, and, secondly, when is the finite section method applicable to $A$?

Prediction of Soil Properties from the Concepts of Energy Transfer in Dynamic Penetration Tests

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Determination of critical slip surfaces using mutative scale chaos optimization

Limit equilibrium is a common method used to analyze the stability of a slope, and minimization of the factor of safety or identification of critical slip surfaces is a classical geotechnical problem in the context of limit equilibrium methods for slope stability analyses. A mutative scale chaos optimization algorithm is employed in this study to locate the noncircular critical slip surface with Spencer’s method being employed to compute the factor of safety. Four examples from the literature—one homogeneous slope and three layered slopes—are employed to identify the efficiency and accuracy of this approach. Results indicate that the algorithm is flexible and that although it does not generally provide the minimum FS, it provides results that are close to the minimum, an improvement over other solutions proposed in the literature and with small relative errors with respect to other minimum factor of safety (FS) values reported in the literature.

Evaluación de la seguridad al deslizamiento de una presa de fábrica en un cimiento con una familia de discontinuidades y con un criterio de rotura con ley de fluencia no asociada

En España existen del orden de 1,300 grandes presas, de las cuales un 20% fueron construidas antes de los años 60. El hecho de que existan actualmente una gran cantidad de presas antiguas aún en operación, ha producido un creciente interés en reevaluar su seguridad empleando herramientas nuevas o modificadas que incorporan modelos de fallo teóricos más completos, conceptos geotécnicos más complejos y nuevas técnicas de evaluación de la seguridad. Una manera muy común de abordar el análisis de estabilidad de presas de gravedad es, por ejemplo, considerar el deslizamiento a través de la interfase presa-cimiento empleando el criterio de rotura lineal de Mohr-Coulomb, en donde la cohesión y el ángulo de rozamiento son los parámetros que definen la resistencia al corte de la superficie de contacto. Sin embargo la influencia de aspectos como la presencia de planos de debilidad en el macizo rocoso de cimentación; la influencia de otros criterios de rotura para la junta y para el macizo rocoso (ej. el criterio de rotura de Hoek-Brown); las deformaciones volumétricas que ocurren durante la deformación plástica en el fallo del macizo rocoso (i.e., influencia de la dilatancia) no son usualmente consideradas durante el diseño original de la presa. En este contexto, en la presente tesis doctoral se propone una metodología analítica para el análisis de la estabilidad al deslizamiento de presas de hormigón, considerando un mecanismo de fallo en la cimentación caracterizado por la presencia de una familia de discontinuidades. En particular, se considera la posibilidad de que exista una junta sub-horizontal, preexistente y persistente en el macizo rocoso de la cimentación, con una superficie potencial de fallo que se extiende a través del macizo rocoso. El coeficiente de seguridad es entonces estimado usando una combinación de las resistencias a lo largo de los planos de rotura, cuyas resistencias son evaluadas empleando los criterios de rotura no lineales de Barton y Choubey (1977) y Barton y Bandis (1990), a lo largo del plano de deslizamiento de la junta; y el criterio de rotura de Hoek y Brown (1980) en su versión generalizada (Hoek et al. 2002), a lo largo del macizo rocoso. La metodología propuesta también considera la influencia del comportamiento del macizo rocoso cuando este sigue una ley de flujo no asociada con ángulo de dilatancia constante (Hoek y Brown 1997). La nueva metodología analítica propuesta es usada para evaluar las condiciones de estabilidad empleando dos modelos: un modelo determinista y un modelo probabilista, cuyos resultados son el valor del coeficiente de seguridad y la probabilidad de fallo al deslizamiento, respectivamente. El modelo determinista, implementado en MATLAB, es validado usando soluciones numéricas calculadas mediante el método de las diferencias finitas, empleando el código FLAC 6.0. El modelo propuesto proporciona resultados que son bastante similares a aquellos calculados con FLAC; sin embargo, los costos computacionales de la formulación propuesta son significativamente menores, facilitando el análisis de sensibilidad de la influencia de los diferentes parámetros de entrada sobre la seguridad de la presa, de cuyos resultados se obtienen los parámetros que más peso tienen en la estabilidad al deslizamiento de la estructura, manifestándose además la influencia de la ley de flujo en la rotura del macizo rocoso. La probabilidad de fallo es obtenida empleando el método de fiabilidad de primer orden (First Order Reliability Method; FORM), y los resultados de FORM son posteriormente validados mediante simulaciones de Monte Carlo. Los resultados obtenidos mediante ambas metodologías demuestran que, para el caso no asociado, los valores de probabilidad de fallo se ajustan de manera satisfactoria a los obtenidos mediante las simulaciones de Monte Carlo. Los resultados del caso asociado no son tan buenos, ya que producen resultados con errores del 0.7% al 66%, en los que no obstante se obtiene una buena concordancia cuando los casos se encuentran en, o cerca de, la situación de equilibrio límite. La eficiencia computacional es la principal ventaja que ofrece el método FORM para el análisis de la estabilidad de presas de hormigón, a diferencia de las simulaciones de Monte Carlo (que requiere de al menos 4 horas por cada ejecución) FORM requiere tan solo de 1 a 3 minutos en cada ejecución. There are 1,300 large dams in Spain, 20% of which were built before 1960. The fact that there are still many old dams in operation has produced an interest of reevaluate their safety using new or updated tools that incorporate state-of-the-art failure modes, geotechnical concepts and new safety assessment techniques. For instance, for gravity dams one common design approach considers the sliding through the dam-foundation interface, using a simple linear Mohr-Coulomb failure criterion with constant friction angle and cohesion parameters. But the influence of aspects such as the persistence of joint sets in the rock mass below the dam foundation; of the influence of others failure criteria proposed for rock joint and rock masses (e.g. the Hoek-Brown criterion); or the volumetric strains that occur during plastic failure of rock masses (i.e., the influence of dilatancy) are often no considered during the original dam design. In this context, an analytical methodology is proposed herein to assess the sliding stability of concrete dams, considering an extended failure mechanism in its rock foundation, which is characterized by the presence of an inclined, and impersistent joint set. In particular, the possibility of a preexisting sub-horizontal and impersistent joint set is considered, with a potential failure surface that could extend through the rock mass; the safety factor is therefore computed using a combination of strength along the rock joint (using the nonlinear Barton and Choubey (1977) and Barton and Bandis (1990) failure criteria) and along the rock mass (using the nonlinear failure criterion of Hoek and Brown (1980) in its generalized expression from Hoek et al. (2002)). The proposed methodology also considers the influence of a non-associative flow rule that has been incorporated using a (constant) dilation angle (Hoek and Brown 1997). The newly proposed analytical methodology is used to assess the dam stability conditions, employing for this purpose the deterministic and probabilistic models, resulting in the sliding safety factor and the probability of failure respectively. The deterministic model, implemented in MATLAB, is validated using numerical solution computed with the finite difference code FLAC 6.0. The proposed deterministic model provides results that are very similar to those computed with FLAC; however, since the new formulation can be implemented in a spreadsheet, the computational cost of the proposed model is significantly smaller, hence allowing to more easily conduct parametric analyses of the influence of the different input parameters on the dam’s safety. Once the model is validated, parametric analyses are conducting using the main parameters that describe the dam’s foundation. From this study, the impact of the more influential parameters on the sliding stability analysis is obtained and the error of considering the flow rule is assessed. The probability of failure is obtained employing the First Order Reliability Method (FORM). The probabilistic model is then validated using the Monte Carlo simulation method. Results obtained using both methodologies show good agreement for cases in which the rock mass has a nonassociate flow rule. For cases with an associated flow rule errors between 0.70% and 66% are obtained, so that the better adjustments are obtained for cases with, or close to, limit equilibrium conditions. The main advantage of FORM on sliding stability analyses of gravity dams is its computational efficiency, so that Monte Carlo simulations require at least 4 hours on each execution, whereas FORM requires only 1 to 3 minutes on each execution.

Caracterización estructural de los rellenos situados en el trasdós de bóvedas de edificios históricos

El papel esencial de los rellenos en el comportamiento y la estabilidad de las construcciones de fábrica ha sido sorprendentemente ignorado en la bibliografía técnica, contribuyendo este trabajo a cubrir esta laguna del conocimiento. Se presenta en primer lugar la configuración morfológica de los rellenos del trasdós de bóvedas en la Península Ibérica en tres periodos histórico-técnicos: la construcción romana, la románica y la gótica. La clasificación se ha realizado en base a inspecciones visuales de diferentes edificios. En sendos anexos se presentan análisis comparados de la base de datos que se ha reunido y las fichas de las construcciones consideradas. En segundo lugar se ha caracterizado el papel estructural de los rellenos, identificándose las características intrínsecas que se estiman más significativas de dichos rellenos en relación tanto con su propia configuración como con su función estructural. Muy en particular se proponen unas leyes de empujes coherentes con la naturaleza y con las condiciones de contorno de los rellenos y las bóvedas. En esta tesis se ha desarrollado una sencilla, potente y suficiente herramienta de análisis que considera las principales variables que intervienen en el comportamiento de una estructura de fábrica, incluyendo muy particularmente los rellenos. Para ello se ha planteado una herramienta basada en el análisis límite que incluye la comprobación tensional de los diferentes esfuerzos, normales o tangenciales, así como la interfaz entre bóvedas y rellenos, dotándose a los técnicos de una herramienta que permite entender estas estructuras, de complejo funcionamiento. La herramienta propuesta incluye consideraciones energéticas para encontrar la solución más cercana al funcionamiento real. La herramienta se ha calibrado con ejemplos analizados por otros autores con diferentes procedimientos, poniéndose de manifiesto su suficiente precisión y eficacia, especialmente si se tiene en cuenta que requiere sólo de unos datos mínimos. La calibración se ha extendido al ensayo en modelo reducido de este tipo de estructuras en el Laboratorio de Estructuras de la Escuela de Ingenieros de Caminos, Canales y Puertos de Madrid. Fruto de esa calibración, se ha realizado un amplio estudio paramétrico que permite cuantificar la influencia que tienen en el comportamiento estructural variables como el apuntamiento de las bóvedas (relación luz/flecha), su esbeltez (relación luz/canto en clave), la altura de los rellenos rígidos, la altura de los rellenos granulares, etc. Ese estudio ha permitido identificar los puntos críticos, las ratios características que condicionan el comportamiento y, en fin, los aspectos clave del análisis estructural de estas construcciones desde el punto de vista de la interrelación, hasta ahora desconocida, de estructura y relleno, que también es estructura. Ello permitirá orientar de mejor manera las inspecciones de estas estructuras y sus análisis estructurales. Finalmente, la tesis concluye con el desarrollo de tres ejemplos, correspondientes cada uno de ellos a los periodos históricos referidos, en los que se pone de manifiesto el papel crucial de los rellenos en la estabilidad no sólo de las bóvedas sino del resto de elementos estructurales que componen un edificio (pilares, muros, botareles, etc.). Al final se presentan las conclusiones del trabajo y se enumeran las muchas líneas de trabajo que quedan abiertas tras esta tesis. It is surprising that the essential role of the backfill at the extrados of masonry vaults and domes has been ignored or misunderstood in the technical literature. The main purpose of this work is to contribute to overcome such a lacuna. Firstly, the morphological configuration of backfill that can be found in the Iberian Peninsula in Roman, Romanesque and Gothic buildings is also presented. Two annexes are devoted to show the collected database and the result of its analysis. Then, the structural role of backfill is characterised, as well as the intrinsic and most significant properties related both to its configuration and structural function. Thus, specific thrust values are proposed to analyse such structures, taking into account the nature and boundary conditions of walls and vaults. The simple, powerful and sufficient tool of analysis that has been developed is also presented, including the main variables involved, especially the backfill. The limit analysis method, followed by a stress verification, including the extrados of vaults and backfill interface, and an enhanced algorithm based on energetic considerations, are used to provide technicians a tool to understand this type of structures efficiently, in a balanced situation between the rather scarce amount of information and sufficient accuracy and liability. The tool was calibrated by comparing the experimental and theoretical results provided by different authors regarding several examples, as well as our own experimental campaign carried out at the Civil Engineering School of Madrid. A parametric study is also developed, providing original information about the structural importance of several parameters, as the span-to-depth or span-to-rise ratios, relative height of backfill, etc. As a result of this analysis, critical points and tendencies of structural behaviour were detected, helping technicians to inspect and analyse these constructions, also revealing some particularities related to their construction. Furthermore, a set of three examples corresponding to the three structural typologies analysed (Roman, Romanesque and Gothic) is studied. The crucial importance of backfill is then highlighted, not only for the stability of vaults or domes, but also for the equilibrium of walls, pillars and buttresses. Finally, the conclusions of the study and the future research fields are also presented.

Unsaturated Infinite Slope Stability Considering Surface Flux Conditions

A slope stability model is derived for an infinite slope subjected to unsaturated infiltration flow above a phreatic surface. Closed form steady state solutions are derived for the matric suction and degree of saturation profiles. Soil unit weight, consistent with the degree of saturation profile, is also directly calculated and introduced into the analyzes, resulting in closed-form solutions for typical soil parameters and an infinite series solution for arbitrary soil parameters. The solutions are coupled with the infinite slope stability equations to establish a fully realized safety factor function. In general, consideration of soil suction results in higher factor of safety. The increase in shear strength due to the inclusion of soil suction is analogous to making an addition to the cohesion, which, of course, increases the factor of safety against sliding. However, for cohesive soils, the results show lower safety factors for slip surfaces approaching the phreatic surface compared to those produced by common safety factor calculations. The lower factor of safety is due to the increased soil unit weight considered in the matric suction model but not usually accounted for in practice wherein the soil is treated as dry above the phreatic surface. The developed model is verified with a published case study, correctly predicting stability under dry conditions and correctly predicting failure for a particular storm.

The 434(G > C) polymorphism in the eosinophil cationic protein gene and its association with tissue eosinophilia in oral squamous cell carcinomas

Objective: The aim of this study was to investigate the prevalence of the Eosinophil cationic protein (ECP)-gene polymorphism 434(G > C) in oral squamous cell carcinoma (OSCC) patients and its association with tumor-associated tissue eosinophilia (TATE), demographic, clinical, and microscopic variables. Methods: The ECP genotypes of 165 healthy individuals and 157 OSCC patients were detected by PCR-RFLP analysis after cleavage of the amplified DNA sequence with enzyme PstI. TATE was obtained by morphometric analysis. Chi-square test or Fisher`s exact test was used to analyze the association of ECP-gene polymorphism 434(G > C) with TATE, demographic, clinical, and microscopic variables in OSCC patients. Disease-free survival and overall survival were calculated by the Kaplan-Meier product-limit actuarial method and the comparison of the survival curves were performed using log rank test. Results: Most of healthy individuals (53.33%) and OSCC patients (57.97%) were heterozygous for the ECP 434(G > C) polymorphism. Based on numerical differences, our results showed that OSCC patients with intense TATE and at least one C allele had a higher frequency of bilateral neck dissection, local recurrence, vascular embolization, involved resection margins, and postoperative radiotherapy. No statistically significant differences on survival rates were found in OSCC patients presenting different ECP 434(G > C) genotypes. Conclusions: These results suggest a tendency towards a poor clinical outcome in OSCC patients with intense TATE and 434GC/CC genotypes, probably due to an ECP genetic variant with altered cytotoxic activity.

Three-dimensional slope stability analysis of South Peak, Crowsnest Pass, Alberta, Canada

South Peak is a 7-Mm3 potentially unstable rock mass located adjacent to the 1903 Frank Slide on Turtle Mountain, Alberta. This paper presents three-dimensional numerical rock slope stability models and compares them with a previous conceptual slope instability model based on discontinuity surfaces identified using an airborne LiDAR digital elevation model (DEM). Rock mass conditions at South Peak are described using the Geological Strength Index and point load tests, whilst the mean discontinuity set orientations and characteristics are based on approximately 500 field measurements. A kinematic analysis was first conducted to evaluate probable simple discontinuity-controlled failure modes. The potential for wedge failure was further assessed by considering the orientation of wedge intersections over the airborne LiDAR DEM and through a limit equilibrium combination analysis. Block theory was used to evaluate the finiteness and removability of blocks in the rock mass. Finally, the complex interaction between discontinuity sets and the topography within South Peak was investigated through three-dimensional distinct element models using the code 3DEC. The influence of individual discontinuity sets, scale effects, friction angle and the persistence along the discontinuity surfaces on the slope stability conditions were all investigated using this code.

Stability analysis of the 2007 Chehalis lake landslide based on long-range terrestrial photogrammetry and airborne LIDAR data

On December 4th 2007, a 3-Mm3 landslide occurred along the northwestern shore of Chehalis Lake. The initiation zone is located at the intersection of the main valley slope and the northern sidewall of a prominent gully. The slope failure caused a displacement wave that ran up to 38 m on the opposite shore of the lake. The landslide is temporally associated with a rain-on-snow meteorological event which is thought to have triggered it. This paper describes the Chehalis Lake landslide and presents a comparison of discontinuity orientation datasets obtained using three techniques: field measurements, terrestrial photogrammetric 3D models and an airborne LiDAR digital elevation model to describe the orientation and characteristics of the five discontinuity sets present. The discontinuity orientation data are used to perform kinematic, surface wedge limit equilibrium and three-dimensional distinct element analyses. The kinematic and surface wedge analyses suggest that the location of the slope failure (intersection of the valley slope and a gully wall) has facilitated the development of the unstable rock mass which initiated as a planar sliding failure. Results from the three-dimensional distinct element analyses suggest that the presence, orientation and high persistence of a discontinuity set dipping obliquely to the slope were critical to the development of the landslide and led to a failure mechanism dominated by planar sliding. The three-dimensional distinct element modelling also suggests that the presence of a steeply dipping discontinuity set striking perpendicular to the slope and associated with a fault exerted a significant control on the volume and extent of the failed rock mass but not on the overall stability of the slope.