Rockfall susceptibility assessment and remote geological mapping with LiDAR point clouds

Characterizing the geological features and structures in three dimensions over inaccessible rock cliffs is needed to assess natural hazards such as rockfalls and rockslides and also to perform investigations aimed at mapping geological contacts and building stratigraphy and fold models. Indeed, the detailed 3D data, such as LiDAR point clouds, allow to study accurately the hazard processes and the structure of geologic features, in particular in vertical and overhanging rock slopes. Thus, 3D geological models have a great potential of being applied to a wide range of geological investigations both in research and applied geology projects, such as mines, tunnels and reservoirs. Recent development of ground-based remote sensing techniques (LiDAR, photogrammetry and multispectral / hyperspectral images) are revolutionizing the acquisition of morphological and geological information. As a consequence, there is a great potential for improving the modeling of geological bodies as well as failure mechanisms and stability conditions by integrating detailed remote data. During the past ten years several large rockfall events occurred along important transportation corridors where millions of people travel every year (Switzerland: Gotthard motorway and railway; Canada: Sea to sky highway between Vancouver and Whistler). These events show that there is still a lack of knowledge concerning the detection of potential rockfalls, making mountain residential settlements and roads highly risky. It is necessary to understand the main factors that destabilize rocky outcrops even if inventories are lacking and if no clear morphological evidences of rockfall activity are observed. In order to increase the possibilities of forecasting potential future landslides, it is crucial to understand the evolution of rock slope stability. Defining the areas theoretically most prone to rockfalls can be particularly useful to simulate trajectory profiles and to generate hazard maps, which are the basis for land use planning in mountainous regions. The most important questions to address in order to assess rockfall hazard are: Where are the most probable sources for future rockfalls located? What are the frequencies of occurrence of these rockfalls? I characterized the fracturing patterns in the field and with LiDAR point clouds. Afterwards, I developed a model to compute the failure mechanisms on terrestrial point clouds in order to assess the susceptibility to rockfalls at the cliff scale. Similar procedures were already available to evaluate the susceptibility to rockfalls based on aerial digital elevation models. This new model gives the possibility to detect the most susceptible rockfall sources with unprecented detail in the vertical and overhanging areas. The results of the computation of the most probable rockfall source areas in granitic cliffs of Yosemite Valley and Mont-Blanc massif were then compared to the inventoried rockfall events to validate the calculation methods. Yosemite Valley was chosen as a test area because it has a particularly strong rockfall activity (about one rockfall every week) which leads to a high rockfall hazard. The west face of the Dru was also chosen for the relevant rockfall activity and especially because it was affected by some of the largest rockfalls that occurred in the Alps during the last 10 years. Moreover, both areas were suitable because of their huge vertical and overhanging cliffs that are difficult to study with classical methods. Limit equilibrium models have been applied to several case studies to evaluate the effects of different parameters on the stability of rockslope areas. The impact of the degradation of rockbridges on the stability of large compartments in the west face of the Dru was assessed using finite element modeling. In particular I conducted a back-analysis of the large rockfall event of 2005 (265'000 m3) by integrating field observations of joint conditions, characteristics of fracturing pattern and results of geomechanical tests on the intact rock. These analyses improved our understanding of the factors that influence the stability of rock compartments and were used to define the most probable future rockfall volumes at the Dru. Terrestrial laser scanning point clouds were also successfully employed to perform geological mapping in 3D, using the intensity of the backscattered signal. Another technique to obtain vertical geological maps is combining triangulated TLS mesh with 2D geological maps. At El Capitan (Yosemite Valley) we built a georeferenced vertical map of the main plutonio rocks that was used to investigate the reasons for preferential rockwall retreat rate. Additional efforts to characterize the erosion rate were made at Monte Generoso (Ticino, southern Switzerland) where I attempted to improve the estimation of long term erosion by taking into account also the volumes of the unstable rock compartments. Eventually, the following points summarize the main out puts of my research: The new model to compute the failure mechanisms and the rockfall susceptibility with 3D point clouds allows to define accurately the most probable rockfall source areas at the cliff scale. The analysis of the rockbridges at the Dru shows the potential of integrating detailed measurements of the fractures in geomechanical models of rockmass stability. The correction of the LiDAR intensity signal gives the possibility to classify a point cloud according to the rock type and then use this information to model complex geologic structures. The integration of these results, on rockmass fracturing and composition, with existing methods can improve rockfall hazard assessments and enhance the interpretation of the evolution of steep rockslopes. -- La caractérisation de la géologie en 3D pour des parois rocheuses inaccessibles est une étape nécessaire pour évaluer les dangers naturels tels que chutes de blocs et glissements rocheux, mais aussi pour réaliser des modèles stratigraphiques ou de structures plissées. Les modèles géologiques 3D ont un grand potentiel pour être appliqués dans une vaste gamme de travaux géologiques dans le domaine de la recherche, mais aussi dans des projets appliqués comme les mines, les tunnels ou les réservoirs. Les développements récents des outils de télédétection terrestre (LiDAR, photogrammétrie et imagerie multispectrale / hyperspectrale) sont en train de révolutionner l'acquisition d'informations géomorphologiques et géologiques. Par conséquence, il y a un grand potentiel d'amélioration pour la modélisation d'objets géologiques, ainsi que des mécanismes de rupture et des conditions de stabilité, en intégrant des données détaillées acquises à distance. Pour augmenter les possibilités de prévoir les éboulements futurs, il est fondamental de comprendre l'évolution actuelle de la stabilité des parois rocheuses. Définir les zones qui sont théoriquement plus propices aux chutes de blocs peut être très utile pour simuler les trajectoires de propagation des blocs et pour réaliser des cartes de danger, qui constituent la base de l'aménagement du territoire dans les régions de montagne. Les questions plus importantes à résoudre pour estimer le danger de chutes de blocs sont : Où se situent les sources plus probables pour les chutes de blocs et éboulement futurs ? Avec quelle fréquence vont se produire ces événements ? Donc, j'ai caractérisé les réseaux de fractures sur le terrain et avec des nuages de points LiDAR. Ensuite, j'ai développé un modèle pour calculer les mécanismes de rupture directement sur les nuages de points pour pouvoir évaluer la susceptibilité au déclenchement de chutes de blocs à l'échelle de la paroi. Les zones sources de chutes de blocs les plus probables dans les parois granitiques de la vallée de Yosemite et du massif du Mont-Blanc ont été calculées et ensuite comparés aux inventaires des événements pour vérifier les méthodes. Des modèles d'équilibre limite ont été appliqués à plusieurs cas d'études pour évaluer les effets de différents paramètres sur la stabilité des parois. L'impact de la dégradation des ponts rocheux sur la stabilité de grands compartiments de roche dans la paroi ouest du Petit Dru a été évalué en utilisant la modélisation par éléments finis. En particulier j'ai analysé le grand éboulement de 2005 (265'000 m3), qui a emporté l'entier du pilier sud-ouest. Dans le modèle j'ai intégré des observations des conditions des joints, les caractéristiques du réseau de fractures et les résultats de tests géoméchaniques sur la roche intacte. Ces analyses ont amélioré l'estimation des paramètres qui influencent la stabilité des compartiments rocheux et ont servi pour définir des volumes probables pour des éboulements futurs. Les nuages de points obtenus avec le scanner laser terrestre ont été utilisés avec succès aussi pour produire des cartes géologiques en 3D, en utilisant l'intensité du signal réfléchi. Une autre technique pour obtenir des cartes géologiques des zones verticales consiste à combiner un maillage LiDAR avec une carte géologique en 2D. A El Capitan (Yosemite Valley) nous avons pu géoréferencer une carte verticale des principales roches plutoniques que j'ai utilisé ensuite pour étudier les raisons d'une érosion préférentielle de certaines zones de la paroi. D'autres efforts pour quantifier le taux d'érosion ont été effectués au Monte Generoso (Ticino, Suisse) où j'ai essayé d'améliorer l'estimation de l'érosion au long terme en prenant en compte les volumes des compartiments rocheux instables. L'intégration de ces résultats, sur la fracturation et la composition de l'amas rocheux, avec les méthodes existantes permet d'améliorer la prise en compte de l'aléa chute de pierres et éboulements et augmente les possibilités d'interprétation de l'évolution des parois rocheuses.

Geological mapping of the Heber, Rhüdener Sattel, Northern Germany

This master thesis describes the geological mapping of an 8 km**2 area of the 300 m high elevation HEBER in Northern Germany which is part of the Rhüdener Sattel (Harzvorland). The geology consists of mesozoic rocks from Buntsandstein to Keuper; major parts are Muschelkalk, partly covered with a thin layer of Quaternary loess. The map includes measurements of strike and dip in various outcrops, drafting a scetch of the probable tectonic structure of Rühdener Sattel.

Geological Evolution of Coastal and Marine Environments off the Campania Continental Shelf Through Marine Geological Mapping - The Example of the Cilento Promontory

The geological evolution of coastal and marine environments offshore the Cilento Promontory through marine geological mapping is discussed here. The marine geological map n. 502 “Agropoli,” located offshore the Cilento Promontory (southern Italy), is described and put in regional geologic setting. The study area covers water depths ranging between 30 and 200 m isobaths. The geologic map has been constructed in the frame of a research program financed by the National Geological Survey of Italy (CARG Project), finalized to the construction of an up-to-date cartography of the Campania region. Geological and geophysical data on the continental shelf and slope offshore the southern Campania region have been acquired in an area bounded northward by the Gulf of Salerno and southward by the Gulf of Policastro. A high-resolution multibeam bathymetry has permitted the construction of a digital elevation model (DEM). Sidescan sonar profiles have also been collected and interpreted, and their merging with bathymetric data has allowed for the realization of the base for the marine geologic cartography. The calibration of geophysical data has been attempted through sea-bottom samples. The morpho-structures and the seismic sequences overlying the outcrops of acoustic basement reported in the cartographic representation have been studied in detail using single-channel seismics. The interpretation of seismic profiles has been a support for the reconstruction of the stratigraphic and structural setting of the Quaternary continental shelf successions and the outcrops of rocky acoustic basement in correspondence to the Licosa Cape morphostructural high. These areas result from the seaward prolongation of the stratigraphic and structural units, widely cropping out in the surrounding emerged sector of the Cilento Promontory. The cartographic approach is based on the recognition of laterally coeval depositional systems, interpreted in the frame of system tracts of the Late Quaternary depositional sequence.

Geological transect across the Northwestern Himalaya in eastern Ladakh and Lahul (A model for the continental collision of India and Asia)

The detailed geological mapping and structural study of a complete transect across the northwestern Himalaya allow to describe the tectonic evolution of the north Indian continental margin during the Tethys ocean opening and the Himalayan Orogeny. The Late Paleozoic Tethys rifting is associated with several tectonomagmatic events. In Upper Lahul and SE Zanskar, this extensional phase is recorded by Lower Carboniferous synsedimentary transtensional faults, a Lower Permian stratigraphic unconformity, a Lower Permian granitic intrusion and middle Permian basaltic extrusions (Panjal Traps). In eastern Ladakh, a Permian listric normal fault is also related to this phase. The scarcity of synsedimentary faults and the gradual increase of the Permian syn-rift sediment thickness towards the NE suggest a flexural type margin. The collision of India and Asia is characterized by a succession of contrasting orogenic phases. South of the Suture Zone, the initiation of the SW vergent Nyimaling-Tsarap Nappe corresponds to an early phase of continental underthrusting. To the S, in Lahul, an opposite underthrusting within the Indian plate is recorded by the NE vergent Tandi Syncline. This structure is associated with the newly defined Shikar Beh Nappe, now partly eroded, which is responsible for the high grade (amphibolite facies) regional metamorphism of South Lahul. The main thrusting of the Nyimaling-Tsarap Nappe followed the formation of the Shikar Beh Nappe. The Nyimaling-Tsarap Nappe developed by ductile shear of the upper part of the subducted Indian continental margin and is responsible for the progressive regional metamorphism of SE Zanskar, reaching amphibolite facies below the frontal part of the nappe, near Sarchu. In Upper Lahul, the frontal parts of the Nyimaling-Tsarap and Shikar Beh nappes are separated by a zone of low grade metamorphic rocks (pumpellyite-actinolite facies to lower greenschist facies). At high structural level, the Nyimaling-Tsarap Nappe is characterized by imbricate structures, which grade into a large ductile shear zone with depth. The related crustal shortening is about 87 km. The root zone and the frontal part of this nappe have been subsequently affected by two zones of dextral transpression and underthrusting: the Nyimaling Shear Zone and the Sarchu Shear Zone. These shear zones are interpreted as consequences of the counterclockwise rotation of the continental underthrusting direction of India relative to Asia, which occurred some 45 and 36 Ma ago, according to plate tectonic models. Later, a phase of NE vergent `'backfolding'' developed on these two zones of dextral transpression, creating isoclinal folds in SE Zanskar and more open folds in the Nyimaling Dome and in the Indus Molasse sediments. During a late stage of the Himalayan Orogeny, the frontal part of the Nyimaling-Tsarap Nappe underwent an extension of about 15 km. This phase is represented by two types of structures, responsible for the tectonic unroofing of the amphibolite facies rocks of the Sarchu area: the Sarchu high angle Normal Fault, cutting a first set of low angle normal faults, which have been created by reactivation of older thrust planes related to the Nyimaling-Tsarap Nappe.

Correlating soil porosity and respective geological unit in Paraíba do sul Valley, Brazil: a geostatistical methodology proposal

This manuscript aims proposing a methodology for correlating soil porosity to the respective geological units using geostatistical analysis techniques, including interpolation data by kriging. The site studied was in Lorena municipality, Paraíba do Sul Valley, southeastern Brazil. Specifically all studies were carried out within an area of 12 km2 located at Santa Edwirges farm. The database comprehended 41 soil samples taken at different geological and geomorphologic units at three different depths: surface, 50 cm and 100 cm depth. The geostatistical analyses results were correlated to a geological mapping specifically elaborated for the site. This mapping accounts for two different geological formations and a geological contact characterized by a shearing zone. The results indicate the existence of a significant relationship between the soil porosity and the respective geological units. The studies revealed that the residual soils from weathered granitic rocks tend to have higher porosities than the residual soils from weathered biotite gneiss rocks, while the soil porosity within the shearing zone is relatively un-sensitive to the respective geological formation. The spatial patterns observed were efficient to evaluate the relationship between the soil porosity, geology unit and the and geomorphology showing a good potential for correlating with others soil properties such as hydraulic conductivity, soil water retention curves and erosion potentials.

Engineering geological data in support of municipal land use planning-a case study in Analandia, southeast Brazil

This paper presents a method for transforming the information of an engineering geological map into useful information for non-specialists involved in land-use planning. The method consists of classifying the engineering geological units in terms of land use capability and identifying the legal and the geologic restrictions that apply in the study area. Both informations are then superimposed over the land use and a conflict areas map is created. The analysis of these data leads to the identification of existing and forthcoming land use conflicts and enables the proposal of planning measures on a regional and local scale. The map for the regional planning was compiled at a 1:50,000 scale and encompasses the whole municipal land area where uses are mainly rural. The map for the local planning was compiled at a 1:10,000 scale and encompasses the urban area. Most of the classification and operations on maps used spatial analyst tools available in the Geographical Information System. The regional studies showed that the greater part of Analandia's territory presents appropriate land uses. The local-scale studies indicate that the majority of the densely occupied urban areas are in suitable land. Although the situation is in general positive, municipal policies should address the identified and expected land use conflicts, so that it can be further improved.

Geological map of HANSSONHORNA, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (hight differences) are accurate to approximately ±10 m.

Geological map of Geltinger Birk, western Baltic Sea, Germany

The geological structure of a Holocene sand spit system and the adjacent Weichselian glacial deposits in the northeastern part of Schleswig-Holstein have been investigated and presented in a geological map. Thin meltwater deposits overlie the glacial tills in the area of the former Beverö lsland in the west. To its north and northeast, the modern Sand spit system is present. Its basal transgression horizon is composed mainly of gravels and boulders, and directly overlie the Pleistocene deposits. Further up the succession, fine graind sands are present, in turn overlain by the coarser grained sands of the barrier bar. To the east, under the protection of the sand spit, gyttyas and peats which sometimes attain large thicknesses have been deposited under lacustrinellagoonal conditions. Closer to the shore, these sediments are covered by marine sands.

Geological map of RYGHNUTEN, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (hight differences) are accurate to approximately ±10 m.

Geological map of WORSFOLDFJELLET, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (height differences) are accurate to approximately ±10 m.

Geological map of JUCKESKAMMEN, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (hight differences) are accurate to approximately ±10 m.

Geological map of GRAMKROKEN, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (hight differences) are accurate to approximately ±10 m.

Geological map of KVITEBOTNEN, Heimefrontfjella, Antarctica (Scale 1:25,000)

Topographic data of this geological map were obtained through stereoscopic aerial photo interpretation. The photogrammetric photo flights were undertaken in 1986 by the Institut für Angewandte Geodäsie, Frankfurt. Horizontal ground control points required for aerial photo interpretation were determined by means of Doppler satellite observation during the 2nd German Neuschwabenland Expedition 1985/86. Vertical ground control points were taken from unpublished map drafts at 1:100 000 scale by Norsk Polarinstitutt, Oslo. The elevation above mean sea level was transferred to Heimefrontfjella barometrically. For this reason assertions concerning the absolute elevation (referred to sea level) are uncertain. Contours and spot heights presented on the map were obtained from the photogrammetric evaluation of the photography taken in 1986; relative elevation data (hight differences) are accurate to approximately ±10 m.