Laser-scan and gravity joint investigation for subsurface cavity exploration — The Grotta Gigante benchmark

Abstract

We have studied a big karstic cave (Grotta Gigante) in northern Italy using an innovative combination of laser-scan and gravity data. We aimed to forward model the gravity anomaly due to the cavity, verify its compatibility with the Bouguer field, and identify the eventual presence of other sources of gravity anomalies. A sensitivity study was performed preliminarily to assess the minimum size of bodies that could be detected by the gravity surveys. The 3D density model of the Grotta Gigante was constructed using as a geometric constraint the laser-scan data set, which mapped the internal morphologies of the cave, and density measurements on collected rock samples. The laser point cloud was reduced in data density, filtered from the outliers, and subdivided into two surfaces representing the vault and the floor of the cave, to correctly define the prism model. Then, a mean density value, obtained from laboratory measurements, was assigned to the prisms. We computed the gravity effect of the model in the same points at which the gravity field had been measured. Excellent correlation was found for the cavity; some gravity anomalies were revealed in the surrounding area of the Grotta Gigante that could be effected by other underground karstic morphologies. We attempted to estimate the probable size and depth of the causative bodies, compatible with the geologic environment. This site testified to the goodness of gravity methods for the exploration of such structures, that is, particularly important for risk assessment in a karstic area. The cave itself, the biggest tourist cave worldwide, represents an upper limit for expected gravity signals. The combination of exact knowledge of the causative body and the related gravity anomalies composed a unique data set (that we released to the public, as a benchmark), useful for testing inversion and forward model gravity algorithms.

Gravity as a tool to improve the hydrologic mass budget in karstic areas

“…We obtained these two surfaces by firstly dividing the point cloud between ceiling and bottom exploiting a surface that locally follows the median plane of the cave. A similar approach was used by Pivetta and Braitenberg (2015) to process the point cloud of laser-scan data from the Grotta Gigante. Then we interpolated the scattered points separately into two regular grids with spatial resolution of 2 m x 2 m. The two surfaces are plotted in Figure B2a and B2b; the color code is proportional to the height of the surface from sea level.…”

Gravity as a tool to improve the hydrologic mass budget in karstic areas

“…We obtained these two surfaces by firstly dividing the point cloud between ceiling and bottom, exploiting a surface that locally follows the median plane of the cave. A similar approach was used by Pivetta and Braitenberg (2015) to process the point cloud of laser-scan data from the Grotta Gigante. Then we interpolated the scattered points separately into two regular grids with spatial resolution of 2 m × 2 m. The two surfaces are plotted in Fig.…”

GTE: a new software for gravitational terrain effect computation: theory and performances

“…Some new developments (SZWILLUS et al 2012), based on adaptive changing the topography resolution with Tesseroids, need to be mentioned here, since they can probably considerably reduce the computational time. Tesseroids is mainly used for geophysical studies at different scales from the very local one, such as the reconstruction and analysis of the Grotta Gigante cave signal (a Karstic cave in the Northern part of Italy PIVETTA and BRAITENBERG 2015), to the regional ones such as the study of the crustal structure in the Andean region (ALVAREZ et al 2014) or the study of the European Alps orogenetic belt (BRAITENBERG et al 2013). All the tests have been performed on a single node of a supercomputer equipped with two 8-cores Intel Haswell 2.40 GHz processors (for a total of 16 cores) with 128 GB RAM.…”