3D geological modelling of the western Aar Massif (external Central Alps, Switzerland)
Authors: Musso Piantelli Ferdinando , Herwegh Marco , Berger Alfons , Wiederkehr Michael , Kurmann Eva , Möri Andreas , Baumberger Roland 
Published: March 2020

&lt;p&gt;3D modelling of complex and irregular geological bodies is an expanding discipline that combines two-dimensional cartographic and structural data managed with GIS technology. This study presents a complete workflow developed to process geological information to build a 3D model of major stratigraphic, structural and tectonic boundaries. The investigated area is located in the western part of the Aar Massif (external Central Alps, Switzerland) characterized by pronounced topographic (600&amp;#8211;&lt;4000&amp;#160;m) relief, making it prone for surface based 3D depth constructions. The workflow comprises four major steps:&lt;/p&gt;&lt;p&gt;(1)&amp;#160; &lt;strong&gt;Generation of 2D polylines in a map view&lt;/strong&gt;: a two-dimensional dataset of sequences of polylines has been generated in ArcGIS (10.3.1) defining the starting dataset for the major stratigraphic and tectonic boundaries of the bedrock units. This dataset has been compiled and integrated by using: (i) GeoCover vector datasets 1:25&amp;#160;000 of the Swiss Geological Survey; (ii) The Geological Special Map 1:100&amp;#160;000 of the Aar Massif and the Tavetsch and Gotthard Nappes of the Swiss Geological Survey; (iii) data from literature; and (iv) additional field work conducted for this study in key-locations.&lt;/p&gt;&lt;p&gt;(2) &lt;strong&gt;Projection of 2D information onto 3D digital elevation model&lt;/strong&gt;: with the 3D structural modelling software Move (Petex/Midland Valley; 2019.1) the boundaries have then been projected on a digital elevation model (swissALTI3D) with 2&amp;#160;m resolution.&lt;/p&gt;&lt;p&gt;(3) &lt;strong&gt;Construction of tectonic cross sections&lt;/strong&gt;: the use of geometric arguments as well as structural measurements allows for projection of these boundaries into a dense regularly spaced network of 2D cross-sections.&lt;/p&gt;&lt;p&gt;(4) &lt;strong&gt;Interpolation of 3D surfaces&lt;/strong&gt;: the surface and cross-sections boundaries can be interpolated by applying 3D projection and meshing techniques resulting in a final 3D structural model.&lt;/p&gt;&lt;p&gt;Generally, steps (2&amp;#8211;4) require iterative adaptations particularly in the case of surface areas being covered by glaciers or unconsolidated Quaternary sediments. In the model, special emphasis is given to visualize the current structural disposition of the western Aar Massif as well as the relative geometric and overprinting relationships of the deformation sequence that shaped the investigated area throughout the Alpine deformation. Finally, since in the investigated area underground data are scarce, an assessment of the relative uncertainties related to input data and is intended to be performed following the approach proposed by Baumberger (2014) and Fer&amp;#324;andez (2005). The workflow presented here offers the chance to gain validation approaches for domains only weakly constrained or with no surface data available, by generating a 3D model that integrates all accessible geological information and background knowledge.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;REFERENCES&lt;/p&gt;&lt;p&gt;Baumberger, R. (2014): Quantification of Lineaments: Link between internal 3D structure and surface evolution 328 of the Hasli valley (Aar massif, central alps, Switzerland), University of Bern, PhD Thesis, unpublished.&lt;/p&gt;&lt;p&gt;Fer&amp;#324;andez, O. (2005): Obtaining a best fitting plane through 3D georeferenced data, Journal of Structural Geology 27, pp. 855&amp;#8211;858&lt;/p&gt;