Geophysical investigation of landslides : a review

International audienceIn the last two decades, shallow geophysics has considerably evolved with the emergence of 2D spatial imaging, then 3D spatial imaging and now 4D time and space imaging. These techniques allow the study of the spatial and temporal variations of geological structures. This paper aims at presenting a current state-of-the-art on the application of surface geophysical methods to landslide characterization and focuses on recent papers (after 1990) published in peer-reviewed International Journals. Until recently, geophysical techniques have been relatively little used for the reconnaissance of landslides for at least two main reasons. The first one is that geophysical methods provide images in terms of physical parameters which are not directly linked to the geological and mechanical properties required by geologists and engineers. The second reason shown through this study probably comes from a tendency among a part of the geophysicists to overestimate the quality and reliability of the results. This paper gave the opportunity to review recent applications of the main geophysical techniques to landslide characterisation, showing both their interest and their limits. We also emphasized the geophysical image characteristics (resolution, penetration depth) which have to be provided for assessing their reliability, as well as the absolute requirements to combine geophysical methods and to calibrate them with existing geological and geotechnical data. We hope that this paper will contribute to fill the gaps between communities and to strength of using appropriate geophysical methods for landslide investigation

A 2D numerical study of the effect of particle shape and orientation on resistivity in shallow formations

International audienceSurficial heterogeneous soils such as till, alluvial fans or slope deposits are difficult to characterize by geotechnical tests due to the presence of decimeter to meter sized pebbles or rocks. The effective resistivity of such two-component medium composed of a percentage of resistive particles embedded in a conductive matrix is given by the Bussian's equation. The application of this equation allows the concentration of resistive particles to be determined if the resistivity values of each component and of the mixture, as well as the cementation exponent m, are known. However, previous theoretical and experimental studies have shown that the effective resistivity is affected by the shape of the particles. The objective of this study is to numerically determine the 2D effects of particle shape and orientation on the resistivity. Two configurations have been considered in the Finite Element modeling: laboratory like measurements and field layout. For circular particles, the numerical results fit the Bussian's equation with an exponent m of 2. Aligned elongated particles induce an anisotropy which can raise or diminish the exponent m, depending on the particle orientation and on the tortuosity of the current paths. Field experiment simulations showed that m varies 2 between 2.5 and 3.1 for an aspect ratio of 5 and that anisotropy resulting from the particle shape has little effect (m close to 2) when this ratio is lower than 2.5. This increase of m with the aspect ratio is in agreement with the theoretical model of Mendelson and Cohen and is consistent with the results of experimental studies. For laboratory measurement simulations, m values vary between 1.3 and 4 for a particle aspect ratio of 5, whatever the resistivity contrast between the particles and the matrix. The difference of results between the two configurations is explained by the paradox of anisotropy

Are Cities of Northern Europe at Risk ?

A benchmark study on seimic risk has been realised on Liege, Belgium, at the request of the Regional Authority. Its main interest is that it deals with the seismic risk of a city in a low seismicity region. The work involves a hazard study based on the recently defined map of seismicity of Belgium, the definition of the individual vulnerability of buildings, the combination of hazard and vulnerability to define risk and static evaluations of connecting details in non engineered structures. For the evaluation of vulnerability, a simplified screening method has been defined. The main conclusion is that in regions where the Peak Ground Acceleration is higher than 0,1g and the building stock does not possess good structural quality, the seismic risk may be considered high.Peer reviewe

Site Amplification Effects in the Ubaye Valley (France): Measurements and Modeling

Local soil conditions can yield to significant variations in ground motion generated by earthquakes. During the last two decades, these amplification effects have been observed by numerous authors for well-documented earthquakes. On the other hand, theoretical models and numerical techniques have been proposed to physically understand site effects. However, few comparisons have been made between observations and theoretical results for well-known 2D structures (sediment-filled valleys). The aim of this paper is to present a study of the response of a valley in the French Alps. The structure (Ubaye valley) was chosen for its moderate dimensions (500 meters wide and 65 meters thick) which allow an accurate determination of the deposit characteristics, and for the relatively high seismicity of the region. The study has included the set-up of a temporary array of five seismological stations, a geophysical survey of the valley to determine the dynamic properties and the geometry of the soft deposits, and numerical modeling (1D and 2D cases) of the response. Comparisons between observed and computed amplifications show a good agreement for particular input motions (corresponding to the SH case). In the other cases, the spectral ratios exhibit a great variability between different groups of similar earthquakes and more developed simulations should be used (2D P-SV, 3D)

Propositions pour le développement d'une méthodologie d'évaluation quantitative de l'aléa rocheux A new methodology for quantitative rock fall hazard assessment

Document de travail (3 pages) Working paperAn innovative methodology is proposed for assessing rock fall hazard. Diffuse hazard can be assessed from rock fall frequency obtained by laser scanner or photogrammetry. For located hazards, the probability of occurrence is difficult to estimate quantitatively, due to two key problems: (a) the poor knowledge of the internal rock mass structure (or geometry); (b) the poor knowledge of the processes which lead to failure. A new approach is proposed to constraint the geometrical and mechanical model of the rock wall, which is based on the dynamic response of the potentially unstable rock compartments to seismic noise. To describe the temporal evolution of the stability and predict the time to failure, a simple rheological law is suggested. It is proposed to fit the rheological parameter of this law, by taking into account the expected number of rock falls in the studied rock wall, derived from the rock fall frequency. The study of the dynamic response is also suggested as an innovative method for monitoring of unstable rock masses

Ambient Seismic Noise and Microseismicity Monitoring of a Prone-To-Fall Quartzite Tower (Ormea, NW Italy)

Remote sensing techniques are leading methodologies for landslide characterization and monitoring. However, they may be limited in highly vegetated areas and do not allow for continuously tracking the evolution to failure in an early warning perspective. Alternative or complementary methods should be designed for potentially unstable sites in these environments. The results of a six-month passive seismic monitoring experiment on a prone-to-fall quartzite tower are here presented. Ambient seismic noise and microseismicity analyses were carried out on the continuously recorded seismic traces to characterize site stability and monitor its possible irreversible and reversible modifications driven by meteorological factors, in comparison with displacement measured on site. No irreversible modifications in the measured seismic parameters (i.e., natural resonance frequencies of the tower, seismic velocity changes, rupture-related microseismic signals) were detected in the monitored period, and no permanent displacement was observed at the tower top. Results highlighted, however, a strong temperature control on these parameters and unusual preferential vibration directions with respect to the literature case studies on nearly 2D rock columns, likely due the tower geometric constraints, as confirmed by 3D numerical modeling. A clear correlation with the tower displacement rate was found in the results, supporting the suitability of passive seismic monitoring systems for site characterization and early waning purposes

Multiconfiguration GPR measurements for geometric fracture characterization in limestone cliffs (Alps)

Until now, geophysical methods have been rarely used to investigate vertical limestone cliffs, mainly due to the extreme conditions for data acquisition. Nevertheless, these techniques are the only available methods which could provide information on the internal state or a rock mass in terms of discontinuities, which play a major role in rock-fall hazards. In this case study, detailed GPR measurements were carried out on a test site with different acquisition configurations deployed on vertical cliff faces. Conventional 2D profiles, common midpoints (CMP) and transmission data were acquired to evaluate the potential of radar waves to improve the characterization of the geometry and properties of the main discontinuities (fractures) within the massif. The results show that the 3D geometry of fractures, which is a crucial parameter for stability assessment, can be retrieved by combining vertical and horizontal profiles performed along the cliff. CMP profiles acquired along the cliff allow a velocity profile to be obtained as a function of depth. Finally, transmission experiments, which generate complex radargrams, have provided valuable and quantitative information on the rock mass, through the modelling of the waves generated. On the other hand, a velocity tomography obtained from the first arrivals travelling through the rock mass from the transmitters to the receivers, shows an image of the investigated zone with a poor resolution

Characterisation of soils with stony inclusions using geoelectrical measurements

Characterisation and sampling of coarse heterogeneous soils is often impossible using common geotechnical in-situ tests once the soil contains particles with a diameter larger than a few decimetres. In this situation geophysical techniques - and particularly electrical measurements - can act as an alternative method for obtaining information about the ground characteristics. This paper deals with the use of electrical tomography on heterogeneous diphasic media consisting of resistive inclusions embedded in a conductive matrix. The adopted approach articulates in three steps: numerical modelling, measurements on a small-scale physical model, and field measurements. Electrical measurements were simulated using finite element analyses, on a numerical model containing a random concentration of inclusions varying from 0 to 40 %. It is shown that for electrode spacing 8 times greater than the radius of inclusions, the equivalent homogeneous resistivity is obtained. In this condition, average measured resistivity is a function of the concentration of inclusions, in agreement with the theoretical laws. To apply these results on real data, a small-scale physical model has been built, where electrical measurements were conducted both on the model and on each phase. From these laboratory measurements, a very satisfying estimation of the percentage of inclusions has been obtained. Finally, the methodology applied to a real experimental site composed of alluvial fan deposits made of limestone rocks embedded in a clayey matrix. The estimated percentage of rock particles obtained via electrical measurements was in accordance with the real grain size distribution

Modal and thermal analysis of Les Arches unstable rock column (Vercors massif, French Alps)

A potentially unstable limestone column (∼1000 m3, Vercors, French Alps) delineated by an open rear fracture was continuously instrumented with two three-component seismic sensors from mid-May 2009 to mid-October 2011. Spectral analysis of seismic noise allowed several resonance frequencies to be determined, ranging from 6 to 21 Hz. The frequency domain decomposition (FDD) technique was applied to the ambient vibrations recorded on the top of the rock column. Three vibration modes were identified at 6, 7.5 and 9 Hz, describing the upper part of corresponding modal shapes. Finite element numerical modelling of the column dynamic response confirmed that the first two modes are bending modes perpendicular and parallel to the fracture, respectively, while the third one corresponds to torsion. Seismic noise monitoring also pointed out that resonance frequencies fluctuate with time, under thermomechanical control. For seasonal cycles, changes in frequency are due to the variations of the bulk elastic properties with temperature. At daily scale, increase in fundamental frequency with temperature has been interpreted as resulting from the rock expansion inducing a closure of the rear fracture rock bridges, hence stiffening the contact between the column and the rock mass. Conversely, the rock contraction induces a fracture opening and a decrease in resonance frequency. In winter, when the temperature drops below 0 ◦C, a dramatic increase in fundamental frequency is observed from 6 Hz to more than 25 Hz, resulting from ice formation in the fracture. During spring, the resonance frequency gradually diminishes with ice melting to reach the value measured before winter

Permafrost extension modeling in rock slope since the Last Glacial Maximum: application to the large Séchilienne landslide (French Alps).

12 pagesInternational audienceRecent dating performed on large landslides in the Alps reveal that the initiation of instability did not immediately follow deglaciation but occurred several thousand years after ice down-wastage in the valleys. This result indicates that debuttressing is not the immediate cause of landslide initiation. The period of slope destabilization appears to coincide with the wetter and warmer Holocene Climatic Optimum, indicating a climatic cause of landslide triggering, although the role of seismic activity cannot be ruled out. A phenomenon which may partly explain the delay between valley deglaciation and gravitational instability is the temporal persistence of thick permafrost layers developed in the Alps since the Last Glacial Maximum (LGM). This hypothesis was tested through 2D thermal numerical modeling of the large Séchilienne landslide (Romanche valley, French Alps) using plausible input parameter values. Simulation results suggest that permafrost vanished in the Séchilienne slope at 10 to 11 ka, 3,000 to 4,000 years following the total ice down-wastage of the Romanche valley at 14.3 ka. Permafrost persistence could have contributed to the failure delay by temporally strengthening the slope. Numerical simulations also show that the permafrost depth expansion approximately fits the thickness of ground affected by gravitational destabilization, as deduced from geophysical investigations. These results further suggest that permafrost development, associated with an ice segregation mechanism, damaged the rock slope and influenced the resulting landslide geometry