Subglacial tunnel valleys in the Alpine foreland: an example from Bern, Switzerland

The morphology of the Alpine and adjacent landscapes is directly related to glacial erosion and associated sediment transport. Here we report the effects of glacio-hydrologic erosion on bedrock topography in the Swiss Plateau. Specifically, we identify the presence of subsurface valleys beneath the city of Bern and discuss their genesis. Stratigraphic investigations of more than 4,000 borehole data within a 430km2-large area reveal the presence of a network of >200m-deep and 1,000m-wide valleys. They are flat floored with steep sided walls and are filled by Quaternary glacial deposits. The central valley beneath Bern is straight and oriented towards the NNW, with valley flanks more than 20° steep. The valley bottom has an irregular undulating profile along the thalweg, with differences between sills and hollows higher than 50-100m over a reach of 4km length. Approximately 500m high bedrock highlands flank the valley network. The highlands are dissected by up to 80m-deep and 500m-broad hanging valleys that currently drain away from the axis of the main valley. We interpret the valleys beneath the city of Bern to be a tunnel valley network which originated from subglacial erosion by melt water. The highland valleys served as proglacial meltwater paths and are hanging with respect to the trunk system, indicating that these incipient highland systems as well as the main gorge beneath Bern formed by glacial melt water under pressur

Overdeepenings in the Swiss plateau: U-shaped geometries underlain by inner gorges

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

LGM glaciations in the northeastern Anatolian mountains: New insights

Barhal Valley belongs to the Çoruh Valley System in the Kaçkar Mountains of northeastern Anatolia. This 13 km long valley is located to the south of the main weather divide and to the east of Mt. Kaçkar, with the highest peak of the mountain range being 3932 m. Today, source of an average yearly precipitation of 2000 mm of moisture is the Black Sea, situated approximately 40 km to the north of the study site. Glaciers of the Last Glacial Maximum (LGM) descended directly from Mt. Kaçkar and reached an altitude of ca. 1850 m a.s.l. (above sea level). In this study, we are exploring whether the position of Barhal Valley to the south of the main weather divide and its east–west orientation have an influence on the existence and expansion of paleoglaciers. Here, we present 32 new cosmogenic 36Cl dates on erratic boulders from the Çoruh Valley System. We reconstructed three geomorphologically well-contained glacier advances in the Barhal Valley, namely at 34.0 ± 2.3 ka, 22.2 ± 2.6 ka, and 18.3 ± 1.7 ka within the time window of the global LGM. Field evidence shows that the glacier of the 18.3 ± 1.7 ka advance disappeared rapidly and that by the latest time, at 15.6 ± 1.8 ka, the upper cirques were ice-free. No evidence for Lateglacial glacier fluctuations was found, and the Neoglacial activity is restricted to the cirques with rock glaciers. A range of 2700 to 3000 m for the Equilibrium Line Altitude (ELA) at the LGM was reported based on modeling of the glacial morphology. We determined that the most likely position of the LGM ELA in the Çoruh Valley System was at 2900 m a.s.l. We suggest an alternative moisture source to the direct transport from the Black Sea for the ice accumulation in the Eastern Black See Mountains. The shift of the Polar Front and of the Siberian High Pressure System to the south during the LGM resulted in the domination of easterly airflow to the Caucasus and Kaçkar Mountains with moisture from expanded lakes in central–western Siberia and from the enlarged Aral- and Caspian Seas

Timing of retreat of the Reuss Glacier (Switzerland) at the end of the Last Glacial Maximum

We used cosmogenic 10Be and 36Cl to establish the timing of the onset of deglaciation after the Last Glacial Maximum of the Reuss Glacier, one of the piedmont lobes of the Alpine ice cap that reached the northern Alpine foreland in Switzerland. In this study, we sampled erratic boulders both at the frontal position in the foreland (Lenzburg and Wohlen, canton Aargau) and at the lateral Alpine border position (Seeboden moraine, Rigi, canton Schwyz). The minimum age for the beginning of retreat is 22.2±1.0ka at the frontal (terminal) position and 20.4±1.0ka at the lateral position. These ages are directly comparable with exposure ages from the other piedmont lobes in the northern Alpine foreland. Our data from the mountain called Rigi, do not support the hypothesis that boulders located external to the Seeboden moraine were deposited prior to the last glacial cycle. We present a first exposure age from an erratic boulder in a retreat position in the Alpine foreland. The Reuss Glacier was approximately 12km behind the maximal extent no later than at 18.6±0.9ka

Unravelling the moisture sources of the Alpine Glaciers using tunnel valleys as constraints

A lack of archives has impeded reconstructions of moisture pathways for past glaciations in the European Alps. Here, we focus on the confluence area of two palaeoglaciers in the Swiss Plateau that were sourced on the northern (Aare glacier) and southern sides (Valais glacier) of the European Alps. We mapped tunnel valleys in the region using a drilling database, based on which we inferred the relative extent of each glacier c. 270 ka ago when the valleys were formed. We then compared this situation with that of the LGM. We found that, while the Valais glacier expanded farther into the foreland than the Aare glacier during the LGM, the opposite was the case c. 270 ka ago. We also found that LGM glaciers were non-erosive in the distal foreland. These contrasts in extents and erosional efficiencies imply differences in moisture pathways between the LGM and the time when the tunnel valleys were formed

Overdeepenings in the Swiss plateau: U-shaped geometries underlain by inner gorges

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.ISSN:1661-8734ISSN:1661-872

Environmental controls on 10 Be-based catchment-averaged denudation rates along the western margin of the Peruvian Andes

We present 10Be-based basin-averaged denudation rates for the entire western margin of the Peruvian Andes. Denudation rates range from c. 9 mm/ka to 190 mm/ka and are related neither to the subduction of the Nazca plate, nor to the current seismicity along the Pacific coast and the occurrence of raised Quaternary marine terraces. Therefore, we exclude a tectonic control on denudation on a millennial timescale. Instead, we explain >60% of the observed denudation rates with a model where erosion rates increase either with mean basin slope angles or with mean annual water discharge. These relationships suggest a strong environmental control on denudation