The AlpArray Gravity Research Group (AAGRG), as part of the European AlpArray program, focuses on the compilation of a homogeneous surface-based gravity data set across the Alpine area. In 2017 10 European countries in the Alpine realm agreed to contribute with gravity data for a new compilation of the Alpine gravity field in an area spanning from 2 to 23^∘ E and from 41 to 51^∘ N. This compilation relies on existing national gravity databases and, for the Ligurian and the Adriatic seas, on shipborne data of the Service Hydrographique et Océanographique de la Marine and of the Bureau Gravimétrique International. Furthermore, for the Ivrea zone in the Western Alps, recently acquired data were added to the database. This first pan-Alpine gravity data map is homogeneous regarding input data sets, applied methods and all corrections, as well as reference frames. Here, the AAGRG presents the data set of the recalculated gravity fields on a 4 km × 4 km grid for public release and a 2 km × 2 km grid for special request. The final products also include calculated values for mass and bathymetry corrections of the measured gravity at each grid point, as well as height. This allows users to use later customized densities for their own calculations of mass corrections. Correction densities used are 2670 kg m⁻³ for landmasses, 1030 kg m⁻³ for water masses above the ellipsoid and −1640 kg m⁻³ for those below the ellipsoid and 1000 kg m⁻³ for lake water masses. The correction radius was set to the Hayford zone O₂ (167 km). The new Bouguer anomaly is station completed (CBA) and compiled according to the most modern criteria and reference frames (both positioning and gravity), including atmospheric corrections. Special emphasis was put on the gravity effect of the numerous lakes in the study area, which can have an effect of up to 5 mGal for gravity stations located at shorelines with steep slopes, e.g., for the rather deep reservoirs in the Alps. The results of an error statistic based on cross validations and/or “interpolation residuals” are provided for the entire database. As an example, the interpolation residuals of the Austrian data set range between about −8 and +8 mGal and the cross-validation residuals between −14 and +10 mGal; standard deviations are well below 1 mGal. The accuracy of the newly compiled gravity database is close to ±5 mGal for most areas. A first interpretation of the new map shows that the resolution of the gravity anomalies is suited for applications ranging from intra-crustal- to crustal-scale modeling to interdisciplinary studies on the regional and continental scales, as well as applications as joint inversion with other data sets. The data are published with the DOI https://doi.org/10.5880/fidgeo.2020.045 (Zahorec et al., 2021) via GFZ Data Services.

We investigated the mechanisms leading to the formation of tunnel valleys in the Swiss foreland near Bern. We proceeded through producing 3D maps of the bedrock topography based on drillhole information and a new gravimetric survey combined with modelling. In this context, the combination of information about the densities of the sedimentary fill and of the bedrock, together with published borehole data and the results of gravity surveys along 11 profiles across the valleys, served as input for the application of our 3D gravity modelling software referred to as PRISMA. This ultimately allowed us to model the gravity effect of the Quaternary fill of the overdeepenings and to produce cross-sectional geometries of these troughs. The results show that 2–3 km upstream of the city of Bern, the overdeepenings are approximately 3 km wide. They are characterized by steep to oversteepened lateral flanks and a wide flat base, which we consider as a U-shaped cross-sectional geometry. There, the maximum residual gravity anomaly ranges between − 3 to − 4 mGal for the Aare valley, which is the main overdeepening of the region. Modelling shows that this corresponds to a depression, which reaches a depth of c. 300 m a.s.l. Farther downstream approaching Bern, the erosional trough narrows by c. 1 km, and the base gets shallower by c. 100 m as revealed by drillings. This is supported by the results of our gravity surveys, which disclose a lower maximum gravity effect of c. − 0.8 to − 1.3 mGal. Interestingly, in the Bern city area, these shallow troughs with maximum gravity anomalies ranging from − 1.4 to − 1.8 mGal are underlain by one or multiple inner gorges, which are at least 100 m deep (based on drilling information) and only a few tens of meters wide (disclosed by gravity modelling). At the downstream end of the Bern area, we observe that the trough widens from 2 km at the northern border of Bern to c. 4 km approximately 2 km farther downstream, while the bottom still reaches c. 300 to 200 m a.s.l. Our gravity survey implies that this change is associated with an increase in the maximum residual anomaly, reaching values of − 2.5 mGal. Interestingly, the overdeepening’s cross-sectional geometry in this area has steeply dipping flanks converging to a narrow base, which we consider as V-shaped. We attribute this shape to erosion by water either underneath or at the snout of a glacier, forming a gorge. This narrow bedrock depression was subsequently widened by glacial carving. In this context, strong glacial erosion upstream of the Bern area appears to have overprinted these traces. In contrast, beneath the city of Bern and farther downstream these V-shaped features have been preserved. Available chronological data suggest that the formation of this gorge occurred prior to MIS 8 and possibly during the aftermath of one of the largest glaciations when large fluxes of meltwater resulted in the fluvial carving into the bedrock.

Astrogeodetic deflections of the vertical (DoVs) are close indicators of the slope of the geoid. Thus, DoVs observed along horizontal profiles may be integrated to create geoid undulation profiles. In this study, we collected DoV data in the Eastern Swiss Alps using a Swiss Digital Zenith Camera, the COmpact DIgital Astrometric Camera (CODIAC), and two total station-based QDaedalus systems. In the mountainous terrain of the Eastern Swiss Alps, the geoid profile was established at 15 benchmarks over a two-week period in June 2021. The elevation along the profile ranges from 1185 to 1800 m, with benchmark spacing ranging from 0.55 km to 2.10 km. The DoV, gravity, GNSS, and levelling measurements were conducted on these 15 benchmarks. The collected gravity data were primarily used for corrections of the DoV-based geoid profiles, accounting for variations in station height and the geoid-quasigeoid separation. The GNSS/levelling and DoV data were both used to compute geoid heights. These geoid heights are compared with the Swiss Geoid Model 2004 (CHGeo2004) and two global gravity field models (EGM2008 and XGM2019e). Our study demonstrates that absolute geoid heights derived from GNSS/leveling data achieve centimeter-level accuracy, underscoring the precision of this method. Comparisons with CHGeo2004 predictions reveal a strong correlation, closely aligning with both GNSS/leveling and DoV-derived results. Additionally, the differential geoid height analysis highlights localized variations in the geoid surface, further validating the robustness of CHGeo2004 in capturing fine-scale geoid heights. These findings confirm the reliability of both absolute and differential geoid height calculations for precise geoid modeling in complex mountainous terrains.