United States

Landslides represent a serious threat in hilly and mountainous areas of the world. Usually, shallow landslides are triggered by intense or prolonged rainfalls and frequently may evolve into potentially catastrophic flow-like movements. Because of the capacity of these phenomena to travel long distances, buildings and infrastructures located in areas improperly deemed safe can be affected. Landslides can be destructive events when mobilize large volumes of soil and/or rock and undergo long distances. It is impossible to stop such events and the best way to prevent disasters is to reduce or avoid the related risk (Carrara et al., 1999). Landslide risk assessment is a research topic of growing interest across the scientific community. Such interest stems from the awareness of the still dramatic impact of landslides on the socio-economic development of vast areas and the concurrent boosting request for urbanization across landslide-prone areas (Petley, 2012). Risk mitigation measures must necessarily focus on the vulnerability of the exposed elements and, in particular, of buildings potentially impacted by landslides. The destruction of these structures can provoke enormous economic losses and endanger the people inside them. On this basis, a study concerning landslide impacting settlements in the Cuenca area (Ecuador) has been carried out thanks to the availability of two fundamental datasets: landslide inventory map (LIM) and information on building&#39;s presence. A LIM represents a fundamental tool to display information on landslides activity and their multi-temporal evolution (Corominas et al., 2014). In most cases, inventories are discontinuous over time and so-called event-based landslide inventories are available, which are landslide inventories based on a singular severe event (rainfall- or earthquake-induced). In this case, LIM was derived from interferometric data and the visual interpretation of aerial photos integrated and validated by field investigation. The surveyed slope failures were classified according to Cruden &amp; Varnes (1996). As concerns buildings information, it was possible to obtain multi-temporal data on the evolution of buildings constructed from 2010 to 2020. From this knowledge, it was possible to note that over time, due to the steady population increase, most of the buildings were built in correspondence of the mountainous areas surrounding the Cuenca center. Thanks to the availability of these data, it can be stated that the buildings, over time, have been built in areas that can potentially be affected by landslides.

SGI 2021

Landslides affecting buildings in rapidly growing areas of Cuenca (Ecuador)

**Attenzione, gli orari riportati nel PROGRAMMA (icona calendario) si riferiscono al fuso orario locale dell’utente. Il fuso orario ufficiale del Congresso è Central European Time (CET).**

L'Istituto Nazionale di Oceanografia e Geofisica Sperimentale - OGS e il Dipartimento di Matematica e Geoscienze dell'Università di Trieste sono orgogliosi di ospitare il 90° Congresso della Società Geologica Italiana a Trieste dal 14 al 16 settembre 2021. La modalità del congresso è interamente online.

Con il motto "Geologia senza frontiere" questo evento mira a promuovere le Scienze della Terra come disciplina scientifica proiettata nel futuro, oltre i confini geografici, politici e disciplinari, per comprendere meglio l'interno del Pianeta Terra e la struttura del sottosuolo, difendere la società da pericolosi eventi geologici, comprendere meglio la variabilità climatica, pianificare l'utilizzo delle georisorse in modo etico e nel rispetto dell'ambiente.


**Please note, times listed in the SCHEDULE tab (calendar icon) reflect YOUR LOCAL TIME ZONE. The main timezone for this conference will be Central European Summer Time.**

The National Institute of Oceanography and Experimental Geophysics - OGS and the Department of Mathematics and Geosciences of the University of Trieste are proud to host the 90th Congress of the Italian Geological Society in Trieste on 14-16 September 2021.  The Congress will be held entirely online.

With the motto "Geology Without Borders" this event aims to promote Earth Sciences as a scientific discipline projected into the future, beyond geographical, political and disciplinary boundaries, to better understand the interior of Planet Earth and the structure of the subsurface, defend society from dangerous geological events, better understand climatic variability, plan the use of geo-resources in an ethical way and respecting the environment.


Welcome to the SGI2021 virtual conference.  In order to access the conference and content, [please register for the meeting here](https://www.geoscienze.org/trieste2021/index.php/registrazioni/registrati)

Please register for the SGI conference to join the event.

Earth Sciences should be understood not only as necessary for a better understanding of the interior of the planet Earth, but also as a defense of society from dangerous geological events, for the understanding of climatic variations, for the planning and use of georesources in a ethical and respecting and protecting the environment.

Landslides represent a serious threat in hilly and mountainous areas of the world. Usually, shallow landslides are triggered by intense or prolonged rainfalls and frequently may evolve into potentially catastrophic flow-like movements. Because of the capacity of these phenomena to travel long distances, buildings and infrastructures located in areas improperly deemed safe can be affected. Landslides can be destructive events when mobilize large volumes of soil and/or rock and undergo long distances. It is impossible to stop such events and the best way to prevent disasters is to reduce or avoid the related risk (Carrara et al., 1999). Landslide risk assessment is a research topic of growing interest across the scientific community. Such interest stems from the awareness of the still dramatic impact of landslides on the socio-economic development of vast areas and the concurrent boosting request for urbanization across landslide-prone areas (Petley, 2012). Risk mitigation measures must necessarily focus on the vulnerability of the exposed elements and, in particular, of buildings potentially impacted by landslides. The destruction of these structures can provoke enormous economic losses and endanger the people inside them. On this basis, a study concerning landslide impacting settlements in the Cuenca area (Ecuador) has been carried out thanks to the availability of two fundamental datasets: landslide inventory map (LIM) and information on building's presence. A LIM represents a fundamental tool to display information on landslides activity and their multi-temporal evolution (Corominas et al., 2014). In most cases, inventories are discontinuous over time and so-called event-based landslide inventories are available, which are landslide inventories based on a singular severe event (rainfall- or earthquake-induced). In this case, LIM was derived from interferometric data and the visual interpretation of aerial photos integrated and validated by field investigation. The surveyed slope failures were classified according to Cruden & Varnes (1996). As concerns buildings information, it was possible to obtain multi-temporal data on the evolution of buildings constructed from 2010 to 2020. From this knowledge, it was possible to note that over time, due to the steady population increase, most of the buildings were built in correspondence of the mountainous areas surrounding the Cuenca center. Thanks to the availability of these data, it can be stated that the buildings, over time, have been built in areas that can potentially be affected by landslides.

technical paper

Geological maps have been historically a central part of the geological knowledge acquisition process: they are used to summarize all the relevant information of a site of interest and they provide a subdivision of an area into geological units. In its simplest form, a digital geological map produced in a Geographical Information System (GIS) typically takes the shape of a polygonal layer, which is needed to visualize the map. For this reason, it is quite common that mappers start new maps by populating a polygonal layer. Although this approach is easier to understand, it inherently has several drawbacks. In particular, the mapper will draw a contact twice for each neighboring polygon, resulting in the need of specialized topological editing tools. Furthermore, the contacts should be stylable (e.g., to represent how much certain we are about the contact). This is not possible with a polygonal representation alone and would require duplicating those contacts as an additional layer, which will, in turn, require additional editing whenever the map is updated. The solution to these issues requires adopting a different approach to avoiding data duplication inherent in a polygonal representation. We hereby present a python plugin for the popular QGIS software, which helps mappers adopt a more topologically-safe approach. Mappy expects its input to be constituted by lines, representing the geological contacts. Provided the lines intersect each other, these suffice to define the partitioning of the map in different units. To assign to each region additional information (i.e., unit’s names) a point layer, containing the needed attributes, is used. Mappy then takes care of updating a polygonal layer by generating the polygons and merging the attributes. This approach, similar to the one suggested by Skinner et al. (2018) has several benefits: any number of segments can be used to define a contact, making it possible to assign different attributes to different portions of the same line. The maps are easier to edit because any time a contact is modified the changes will automatically propagate to all affected polygons, without resorting to specialized topological editing tools. Furthermore, the contacts can be reused for map layout and directly styled. In this way, the mapper can draw the contacts working on just two units at a time, and the map can be progressively refined with no significant effort. These ideas can be implemented in any GIS of choice by performing a polygonization of the contacts and a spatial join of the points to assign the attributes, or Mappy can be used. Mappy enables the operator to adopt a polygons-free strategy during the whole creation of the map and provides the algorithms needed to re-generate the polygonal layer at any time. It also comprises additional tools for clipping the contacts at the intersection points and for automatic placement of the points. Once the base layers are crated (points and contacts) Mappy will update the resulting polygonal layer automatically, which style can be defined progressively during the creation of the map. Mappy has been already tested during the PLANMAP-GMAP mapping winter school, simplifying the generation of consistent geological maps to the students, and demonstrating the benefits of this approach, especially for newcomers which were able to produce new maps without the need of understanding complex topological editing tools. This tool will be further developed to support the upcoming mapping efforts during the JUICE mission, ESA ExoMars, and other missions interacting with the GMAP project.

Mappy: a python plugin to ease geological mapping with QGIS

In this work  we propose the geologic mapping of the Copernicus crater and the characterization of a potential landing site for a rover mission. Moreover, we have selected 68 scientific targets for in situ analyses and sampling, which can answer the ESA’s scientific questions about the history and geology of the Moon. Copernicus Crater is a relatively young impact complex crater, located on the near side and of about 100km indiameter. Studying Copernicus means having the possibility to get information about lunar crust, lunar interior, cratering process, and lunar explosive volcanism as well as constraining the age of a reference stratigraphic plane of the whole Moon, with an important impact on the dating system applied to the entire inner Solar System. Olivine spectral signature was detected (Liu F. et al., 2011), and this can provide information about the lunar crust and mantle. Absolute-dating Copernicus would improve the lunar chronological curve, used to date surfaces of all the bodies of the inner Solar System, through the crater size frequency distribution. Improving this curve means improving the age estimations of all these bodies. Moreover, this work can assist the foreseen human assisted robotic missions. As an incremental update of the work performed on Tusberti, (2019) we firstly performed a large-scale geological characterization and mapping of Copernicus crater, using remote sensing techniques, analyzing topographic, imaging and spectral data coming from Lunar Reconnaissance Orbiter (LRO), Selenological and Engineering Explorer (SELENE), and Clementine missions where many interesting features were visible (like cinder cones, vents, melt differentiated, fresh craters and presence of olivine). Therefore, we generated two geological maps, one covering the entire Copernicus crater (at a scale of 1:400000) and one focused on its NW quadrant (scale 1:100000). In fact, by analyzing the topography in terms of slope and roughness we infer that this quadrant could be the best candidate for the landing area. Based on those elaborations and on a preliminary hazard study (following the constraints from Pajola M. et al., 2017), we selected a tentative landing ellipse with in mind the engineering constraints of the now deprecated HERACLES mission and the forthcoming EL3 (European Large Lunar Lander). Moreover, the landing ellipse was selected not only for safety reasons but also because 68 potential sampling sites which are located in reachable areas. These sampling sites are carefully selected in order to fulfill the major ESA questions related to Moon exploration. 

Geologic mapping and landing site characterization in Copernicus Crater (Moon)

Rembrandt, a 715 km-diameter basin on Mercury centered at 32.89°S,  87.86°E, constitutes a remarkable example of the stratigraphic variability on Mercury. The uppermost smooth plains of volcanic origin, which are temporally unrelated to the impact event, have been distinguished by the surrounding smooth plains utilizing spectral information. Moreover, it is crosscut by a prominent thrust system more than 1000 km long with NE-SW orientation. This consists of a set of lobate scarps and wrinkle ridges with inverse to transpressional kinematics (Massironi et al., 2015) and is named Enterprise Rupes (Ferrari et al., 2015; Galluzzi et al., 2015; Semenzato et al., 2020). Its formation is attributed mainly to crustal shortening induced by global contraction (see Semenzato et al., 2020 and references therein). The aim of this work, as part of the PLANMAP Horizon 2020 project, is to use all the available geologic information to provide a comprehensive 3D model of both the structural setting and the thicknesses of the main units of the Rembrandt basin infilling using explicit modeling approach.  While implicit modelling and its effectiveness rely on the amount of available data both on the surface and the subsurface to provide a best-fit solution that satisfies all the constraints, the explicit geological modelling implies a more extended geologic interpretation based on the available data, which is also subject to the operator’s experience. In practical terms, this modeling approach is essentially based on the construction of several interpreted parallel and intersecting cross-sections that are then interpolated generating mesh surfaces. The advantage is that it can be applied in areas with low amount of data, especially with few subsurface constraints. We mainly based our reconstruction on the geo-stratigraphic map and geologic cross-section of Rembrandt basin from Semenzato et al. (2020). Additional constraints, regarding fault kinematics and fault plane geometry were derived by Galluzzi et al. (2015) and Crane (2020). As a basemap, we used the MESSENGER MDIS global mosaic at 166 m/pixel together with the MDIS global DEM (665 m/pixel). The base unit of Rembrandt (named IT - para-autoctonous intracrater plains) is overlaid by two different events of infilling: the YIP (Young Interior smooth Plains) and OIP (Older Interior smooth Plains) units. Their thickness was estimated in Semenzato et al., (2020),and among the methods one implies the presence of impact craters excavating and exposing the inner basin stratigraphy. We also used such craters for the construction of the interpreted cross-sections in complement to the one presented in that work. We built 14 cross sections that were subdivided into two intersecting sub-parallel sets. The interpreted subsurface horizons were used for mesh interpolation together with the geologic contacts between the different units on the surface creating meshes representin ghte top of the YIP, OIP and IT basin infilling units. The volume between them was then calculated by means of a tetrahedral mesh. The structural modelling was performed by means of the structural mapping from Semenzato et al., (2020) and Ferrari et al., (2015). The faults’ vergence and their angles were extracted thanks to the illumination conditions on the image mosaic and topography, and also to the measurements provided by Galluzzi et al., (2015) on displaced impact craters, to the work by Crane (2020) and Massironi et al., (2015). The Enterprise rupes thrust geometry was constructed including its splays and backthrusts, as well as the contractional lobate scarps and extensional radial features within the basin units. Also in planetary geology and with relative scarcity of data, the 3D geomodelling proves to be a useful tool in order to check surface interpretations and mapping, that might improve or change once reconstructed in 3D and helps also to refine volume calculations and the shape of geoogic bodies. Further work will focus on forward and backwards modeling of the Enterprise Rupes system deformation.

3D geo-model of Rembrandt basin on Mercury: structural framework and infilling 

Noctis Labyrinthus (NL) is located in the western part of the large extensional province of Valles Mariners and close to the vast volcanic plateau of Tharsis Montes. Within this area, linear pit-chains commonly occur in close vicinity to faults, inside the graben, and mainly connected to swarms of small cracks, generating a complex pattern of interconnected faults and a widespread grid of scalloped pit-chains. Such a striking spatial correlation between faults and pit-chains makes this area a location of interest to study the possible mechanisms of formation and evolution of these features. Previous work interpreted NL as a karst landscape with associated caves and water-related processes (Baioni et al.,2017; Weitz et al., 2016; Jackson et al.,2011; Thollot et al.,2010; Grindrod et Balme.,2010; Bérezi., 2008; Murchie et al.,2007; Johnston et al.,2006; Bibring et al.,2006; Boston.,2004;), a network of lava tubes (leone, 2014) or the result of volcano-tectonic activity (Mège et al, 2003, Schultz, 1998). In this work, we used the orthoimages H3210_0000 and H3221_0000 acquired by the High-Resolution Stereo Camera (HRSC) onboard Mars Express, bearing a resolution of 12.5 m/pixel, and for the surface topography, we utilized an individual Digital Elevation Model (DEM) from Mars Orbiter Laser Altimeter (MOLA) onboard Mars Global Surveyor, with a resolution of ~460 m/pixel. We have mapped all apparent faults and grabens within the selected area (N=3781). Each fault was traced by a single polyline or by a group of polylines when segmented, and we distinguished their size by assigning different colors, while the pit-chains were delimited by a polygon shapefile we performed a detailed analysis of the faults patterns and geometries over different lengths with the aim to provide a comprehensive model for the formation of the complex structural features in Noctis Labyrinthus. According to the occurrence of faults and pit-chains, a clear subdivision can be distinguished between the northern and the southern sectors of the studied area; The northern sector of Noctis Labyrinthus basically consists of smooth terrains, lacking troughs and pit-chains, while in the southern sector, pit-chains are well developed and mostly parallel to the ENE-WSW fault system. The azimuth distribution of fault populations shows two main trends in the rose diagrams: ENE-WSW and NS. The NS faults systems are apparently cross-cut by the ENE-WSW faults system, which seems to occur before pit chains formation. We infer a development of two fault systems potentially produced in response of bidirectional or radial extension, where distinct faults develop perpendicular to each extensional direction. The extension is driven by invariant regional stress field, in pure shear, where faults grow at different scales (Schultz et al., 2000; Turabi and Berg., 2011). We propose that a Synchronous or 2-Phase bidirectional regional extension can be the prime processes behind the formation of the rift systems in NL and the generation of T, L, and X-Shaped intersections between faults (Collanega et al., 2020). The presence of pit chains remains difficult to explain within this setting: the interaction with magmatic activity from the large nearby Syria Planum magmatic complex might have played a role in shaping these structures. In this scenario, dikes intruded below the graben floors might have exploited the local extensional setting creating dilatant conditions in correspondence of major structures. A subsequent pressure drop in the plumbing system might have led to a surface collapse forming the large pit chains. In conclusion, we hypothesize that the complex graben and pit chains system of NL is the combined result of the activity of a bidirectional extension likely in conjunction with the activity of a magmatic plumbing system.

A Volcano-tectonic activity: possible scenario beyond the formation of the rift systems in Noctis Labyrinthus (Mars)

The spectral characterization of planetary surfaces and other solar system objects allows us to obtain geochemical information useful for geological mapping and for understanding the evolutionary history of planets and other celestial objects. Furthermore, identifying planetary surfaces geochemical properties, allow us to characterize volcanic features in order to investigate new volcanic mechanism anomalous on planet Earth, but common in other places. The goal of this work is to geochemically characterize planetary surfaces, through the investigation of their mid-infrared range spectra. We propose a spectral characterization of geological properties by means of comparison between planetary surfaces emissivity spectra, measured by probes, and analog rock samples emissivity spectra, measured in laboratory, obtained from natural melting samples. Principal component analysis, clustering and categorization machine learning algorithms and statistical methodology are utilized to spectral comparison as well as spectral data analysis. These techniques allow us to analyse a large amount of data. We have analysed about 1000 emissivity spectra of a Martian basaltic lava flow located at west of Arsia Mons (3° S, 221,8° E, IAU 2000), at the T=273-300K; about 500 emissivity spectra of Nili Patera, a Martian caldera situated to Syrtis Major Planum (9° N, 67° E IAU 2000), at T=273-300K; both groups measured by Thermal Emission Spectrometer (TES) onboard NASA’s Mars Global Surveyor. Regarding our sample spectra, we have analysed 11 emissivity spectra of samples synthesized at the UNIPG laboratory, measured in reflectance and emissivity at the temperatures of 450, 600, 750 and 900K; other 80 emissivity spectra of rock samples from dataset of the Arizona State University were also analysed. On Mars, magmatic differentiation is present, as occurred in Nili Patera caldera, where THEMIS multispectral images show an evident spectral contrast due to magmatic differentiation. The same contrast is identifiable in our TES spectral analysis. The preliminary results show similar trend for chemical component that we have been analysed. Once chemical and spectral features of our samples will be better characterized, we will be able to define quantitative geochemical constraints in order to make geological maps based on mid-infrared measures.

Comparison of emissivity spectra from Mars surface and laboratory rock samples: preliminary results

The dark terrain of Ganymede is the oldest geologic unit of the satellite, subsequently crosscut by the youngest light terrain. It shows morphotectonic structures called furrows, which are radial or circumferential troughs bounded by high albedo rims and a lower albedo floor. It is assumed that furrows represent a multi-ring basin formed by an ancient impact crater (McKinnon & Melosh, 1980; Prockter et al., 2000). In fact, Galileo Regio’s furrows show a concentric pattern centered at around 20° S 180° W, which is the assumed location of the crater (Prockter et al., 2002; Hirata et al., 2020). Most of them have been preserved in their pristine form, while some have acted as zones of weakness and have been re-used by later tectonic deformation. These zones are likely considered as precursors of grooved terrain formation. In this contribution, we present the structural mapping at the regional scale of the furrows in Galileo Regio, and we identify their hierarchy and tectonic setting originated by strike-slip and shortening. Riedel shear, synthetic and antithetic structures, and conjugate fractures have been identified. Hence, the performed mapping allows to infer the tectonic setting of the furrows and their associated structures, and in turn the (paleo)stress field that affected the area. Galileo Regio has been deformed by the same regional kinematics consistent with the right-lateral transpression that has affected Uruk Sulcus (Rossi et al., 2018). We infer a ≈ NE-SW trending maximum horizontal stress responsible for a shortening associated with strike-slip in the area extending from 180°-120° W to 0°-60° N. This work confirms the key role that transpression has played in the leading hemisphere of Ganymede and further statistical analysis will be performed to support this reconstruction. The produced tectonic model will be used for the scientific preparation of future high-resolution observation that will be performed by the JANUS instrument onboard JUICE mission (Palumbo et al., 2014).

Strike-slip and shortening in the dark terrain of Ganymede

The surface of Ganymede appears as a deeply tectonized environment mainly composed of a brittle icy crust overlying a subsurface global liquid ocean, whose depth was constrained to be between 100 and 150 km (Lucchetti et al., 2020 and references therein). The surface of Ganymede can be subdivided into dark and bright terrains depending on their crater density, differences in albedo and surface morphology. Indeed, the majority of the terrains can be subdivided according to the presence of furrows and grooves features. These structures represent the brittle deformation of the satellite’s icy shell. The study location is between 150W -180W and 30N-10 S in the region of Uruk Sulcus, a NW-SE terrain ~400 km wide and ~2500 km long strongly deformed by grooves. Several hypotheses account for the presence of such structures, namely tilt-block normal faulting (Pappalardo et al., 1998), crustal necking similar to rifting stage in mid ocean ridges (Pizzi et al., 2017 and references therein) or - in the case of Uruk Sulcus - the effect of a major dextral transpression  (Rossi et al., 2020). Also, it has been demonstrated that the spatial distribution of such structures shows a fractal spatial clustering, that implies that the fracture network below the surface can be interconnected and percolating, reaching up to the ~100-150 km-deep ocean (Lucchetti et al., 2020). To better constrain the transpression hypothesis in Uruk Sulcus, we are testing different 3D geomodelling approaches hence deriving more information on both the volume of the fractured brittle ice crust and the most favored location of actual percolation. At present day, 3D geolocial modelling is approaching the planetary geology, as also testified by the outcomes of the PLANMAP Horizon 2020 project (Massironi et al., 2018, 2021; Pozzobon et al., 2020; Penasa et al., 2020; Pozzobon et a., 2021). With the information either derived from direct surface observations and measurements and/or also constrained by geologic interpretation, it is possible to derive a 3D geomodel of the subsurface that represents specific geologic properties or geometries. Such approach allows to refine or confirm the geologic interpretations or provide new useful qualitative and quantitative insights. A preliminar model of Uruk sulcus was built using the global DEM by Zubarev et al., (2017)  as an initial, low-resolution surface. This served as an upper surface boundary to project the structural mapping from Rossi et al., (2020) and also as an anchor point to reconstruct the major Uruk sulcus structures. It was then textured it with the controlled global mosaic from DLR of Voyager 1, 2 and Galileo images at an average resolution of 359 m/pixel. As a lower model boundary, the icy crust is constrained at depth by the ocean-ice interface. In particular we used the depth value calculated using the fractal clustering analysis in Lucchetti et al., (2020) related specifically to Uruk sulcus. We were able to reconstruct a preliminary qualitative model of the fault geometries for GUS1 (Grooves of Uruk Sulcus 1) system deformation (Rossi et al., 2020) as a large-scale structural framework. We performed a first test on a smaller region in the sigmoidal area named S2 and containing the GUS3 (also in Rossi et al., 2020). We extracted a voxel-based volume of icy crust comprised between the surface, the lower ocean-ice interface bounded by the major fault lines delimiting GUS3. At this location a major set of NE-SW high-angle extensional structures arranged in a right-stepping en-echelon geometry are present together with other two grooves families with NW-SE and NS orientations, although lower in number. By exploiting the voxel-based volume and knowing the scaling laws that rule of the size distributions of the grooves and their azimuths, it is possible to  simulate a DFN (Digital Fracture Network) for each of the fractures’ families and to predict possible high-connectivity locations in the fractured ice volumes and their location at the surface. Further analyses will attempt to simulate the movement along the major structures of Uruk backwards and forward in time, in order and to capture the strain in the ice volume and also test the transpression vs necking approach. The 3D geomodelling approach proves to be a promising tool not only for scientific analysis but also to identify interesting target of future JUICE mission (Grasset et al., 2013) observation through the JANUS instrument (Palumbo et al, 2014), where resurfacing or deep fluid expulsions could be more favored.

3D geo-model of Uruk Sulcus structures (Ganymede): insights and implications

Multidisciplinary efforts are currently boosted by an urgent need to improve our understanding of Earth’s lava systems in sight of their future exploration on other rocky bodies of the Solar System. Among the many structures documented in lava fields, lava tubes are one of the most enigmatic. They constitute a peculiar type of caves dug by molten lava flows. These roofed conducts are very efficient thermal structures enabling channelizing lava transport over long distances. The longest lava tubes are found on volcanic plateau characterised by a gentle slope (<2°) or on volcanic islands (e.g. Hawai’i, Canaries, Iceland, etc.). These structures are easily recognizable from the surface by skylights and collapses on the roof forming pits chains, which allow the reconstruction of the tube path. Due to the similar characteristics of basaltic volcanism on other rocky bodies (e.g. Mars and the Moon), it is expected that lava tubes have similar origins and morphologies to those on Earth. Indeed, analogous conformations of aligned collapses have been seen on the surfaces of Mars and the Moon (Haruyama et al., 2012). This discovery has led to a growing interest as these structures could be suitable sites for future exploration and/or permanent human settlement. Located in the north-eastern part of the island of Lanzarote (Canary Islands), the La Corona lava tube system with its 7.6 km of total length (~8.9 km of cave development) and 10-20 m diameter is one of the terrestrial largest volcanic cave complexes. Accordingly, to the particular geological context in which it arose, La Corona system is one of the most studied. Indeed, the Canarian archipelago represents both long-term and spatially focused volcanic activity over a poorly mobile tectonic plate [less than ~2 cm/yr, during the last 30 Ma (Gaina et al., 2013)]. This environment identifies the Canaries as one of the best analogues of the Martian one-shell plate volcanism on Earth (Meyzen et al., 2015). What makes this inflated lava tube so interesting is a pyroclastic layer, derived by the initial Strombolian activity of La Corona vent (Carracedo et al., 2013) and interleaved within the lava flows crossed by the tube. The layer follows the tube for at least one-third of its extent and we speculate that it could have been pivotal for the inception of the inflation process. By analogy, similar geological settings could be favourable for the formation of lava tubes on rocky bodies like Mars and the Moon, where weak layers of pyroclastic deposits or fine regolith are thought to be common.

Development of inflated lava tubes in analogous planetary environments: the case of La Corona system (Lanzarote, Canary Islands)

Research in urban geomorphology is a relatively recent topic and its theory and practice require continuous updates. Few studies have focussed on issues specifically connected with the cartographic representation of geomorphological processes and landforms in urban environments, with the aim of tentatively reconstruct and clarify the original physical landscape. Among these, recent urban geomorphological surveys conducted in Italian cities proposed new strategies for detecting and mapping human induced modifications to the topography, and anthropogenic landform classification has been recently deepened in the revised Italian national guidelines for geomorphological mapping (Campobasso et al., 2018), proposed by the Working Groups of the Italian Association of Physical Geography and Geomorphology (AIGeo) and the Italian Institute for the Environmental Protection and Research (ISPRA). In this work, natural landforms modified by man and artificial landforms are explored and compared in three Italian cities: Rome, Genoa and Perugia. The comparison is interesting because of the different morphogenetic environment characterizing the 3 cities and also for the differences in the geomorphological survey and mapping. At the same time, in these urban areas, human activities have similarly remodelled the topographic surface and become the most important modelling agent of relief. Many reliefs have been erased, whereas thick layers of anthropogenic deposits covered most of the natural landforms, with the result, for example, of hiding or greatly modifying the natural valley incisions. Geomorphological studies in urban environments allow the reconstruction of the physical landscape before human modifications and the implementation of innovative methodological tools, and an effective legend to represent the urban landforms and their evolution over time.

Geomorphological classification of the landscape in urban areas: anthropogenic landforms in Genoa, Rome and Perugia

Since early 1990s, the advent of satellite-based techniques, such as the Differential Interferometry Synthetic Aperture (DInSAR) (e.g. Curlander & McDonough, 1991) and Global Navigation Satellite System (GNSS), has revolutionized the monitoring of ground level displacements thanks to satellite-based techniques. These observation systems enabled us to control with unprecedented accuracy and resolution deformation phenomena caused by transient events, such as earthquakes, volcanic eruptions and sinkholes, or long-term processes, like subsidence and landslides, thus revealing useful applications in the geological, geotechnical and structural analysis fields (e.g. Lanari et al., 1998; Herrera et al., 2009). Today, the Advanced DInSAR (A-DInSAR) technique is based on high resolution sensors, which monitor the evolution of instability phenomena in any atmospheric and lighting condition, and in both historical centers and archaeological sites. Within this framework, we propose the application of an exhaustive approach for studying both permanent and periodic components of vertical land deformation in the urban area of Naples (southern Italy). Such methodology, successfully applied in the Po Delta area (Vitagliano et al., 2020), is proposed to study the geological risks in urban areas and to evaluate the effects of the anthropic context on the periodic processes related to climate variability. Firstly, a spatial analysis has been carried out to calculate the interferometric velocities of Sentinel-1 and COSMO-SkyMed data in the time intervals January 2016 – December 2019 and March 2017 – December 2020, respectively. The obtained ground displacement and velocity have been used for mapping the stable and unstable urban areas and for evaluating the dependencies between the distribution of the velocity zones and some geological features (e.g. ancient alluvial deposits and paleo-coast lines). Afterwards, a spatial-temporal analysis of continuous GPS data available for a site located in the city center has been performed, which provided a reliable identification of the seasonal trends of the ground displacement. Finally, an automatic cross-correlation procedure developed in Matlab environment was applied to systematically correlate the interferometric trends of each reflection point with the rainfall and sea level oscillations, measured in the same time interval. The correlation process allowed the identification of the most significant dependencies and the triggering mechanisms of the slow soil oscillations, and gave also interesting insights into urban instability processes, such as sinkholes.

Landslides affecting buildings in rapidly growing areas of Cuenca (Ecuador)

Downloads

Next from SGI 2021

Mappy: a python plugin to ease geological mapping with QGIS

Similar lecture

Identifying areas suitable for Sustainable Drainage Systems and Aquifer Storage and Recovery to mitigate stormwater flooding phenomena in Rome (Italy)

Landslides affecting buildings in rapidly growing areas of Cuenca (Ecuador)

.css-70qvj9{display:-webkit-box;display:-webkit-flex;display:-ms-flexbox;display:flex;-webkit-align-items:center;-webkit-box-align:center;-ms-flex-align:center;align-items:center;}Downloads

Next from SGI 2021

Mappy: a python plugin to ease geological mapping with QGIS

Similar lecture

Identifying areas suitable for Sustainable Drainage Systems and Aquifer Storage and Recovery to mitigate stormwater flooding phenomena in Rome (Italy)

Downloads