A gravimetric assessment of the Gotthard Base Tunnel geological model: insights from a novel gravity terrain-adaptation correction and rock physics data

Abstract

The Gotthard Base Tunnel (GBT) is a 57 km long railway tunnel, constructed in the Central Alps in Switzerland and extending mainly North–South across numerous geological units. We acquired 80 new gravity data points at the surface along the GBT profile and used 77 gravity measurements in the tunnel to test and constrain the shallow crustal, km-scale geological model established during the tunnel construction. To this end, we developed a novel processing scheme, which computes a fully 3D, density-dependent gravity terrain-adaptation correction (TAC), to consistently compare the gravity observations with the 2D geological model structure; the latter converted into a density model. This approach allowed to explore and quantify candidate rock density distributions along the GBT modelled profile in a computationally-efficient manner, and to test whether a reasonable fit can be found without structural modification of the geological model. The tested density data for the various lithologies were compiled from the SAPHYR rock physical property database. The tested models were evaluated both in terms of misfit between observed and synthetic gravity data, and also in terms of correlation between misfit trend and topography of the target profile. The results indicate that the locally sampled densities provide a better fit to the data for the considered lithologies, rather than density data averaged over a wider set of Alpine rock samples for the same lithology. Furthermore, using one homogeneous and constant density value for all the topographic corrections does not provide an optimal fit to the data, which instead confirms density variations along the profile. Structurally, a satisfactory fit could be found without modifying the 2D geological model, which thus can be considered gravimetry-proof. From a more general perspective, the gravity data processing routines and the density-dependent corrections developed in this case study represent a remarkable potential for further high-resolution gravity investigations of geological structures.