The stratigraphic architecture of the Swiss Molasse basin, situated on the northern side of the evolving Alps, reveals crucial information about the basin's geometry, its evolution, and the processes leading to the deposition of the siliciclastic sediments. Nevertheless, the formation of the Upper Marine Molasse (OMM) and the controls on the related Burdigalian transgression have still been a matter of scientific debate. During the time period from ca. 20 to 17 Ma, the Swiss Molasse basin was partly flooded by a shallow marine sea striking SW–NE. Previous studies have proposed that the transgression occurred in response to a rise in global sea level, a reduction of sediment flux, or an increase in tectonically controlled accommodation space. Here, we readdress this problem and extract stratigraphic signals from the Burdigalian molasse deposits that can be related to changes in sediment supply rate, variations in the eustatic sea level, and subduction tectonics. To achieve this goal, we conducted sedimentological and stratigraphic analyses of several sites across the entire Swiss Molasse basin.
Field investigations show that the transgression and the subsequent evolution of the Burdigalian seaway was characterized by (i) a deepening and widening of the basin, (ii) phases of erosion and non-deposition during Lower Freshwater Molasse (USM), OMM, and Upper Freshwater Molasse (OSM) times, and (iii) changes in along-strike drainage reversals. We use these changes in the stratigraphic record to disentangle tectonic and surface controls on the facies evolution at various scales. As the most important mechanism, rollback subduction of the European mantle lithosphere most likely caused a further downwarping of the foreland plate, which we use to explain the deepening and widening of the Molasse basin, particularly at distal sites. In addition, subduction tectonics also caused the uplift of the Aar massif. This process was likely to have shifted the patterns of surface loads, thereby resulting in a buckling of the foreland plate and influencing the water depths in the basin. We use this mechanism to explain the establishment of distinct depositional settings, particularly the formation of subtidal shoals wherein a bulge in relation to this buckling is expected. The rise of the Aar massif also resulted in a reorganization of the drainage network in the Alpine hinterland, with the consequence that the sediment flux to the basin decreased. We consider this reduction in sediment supply to have amplified the tectonically controlled deepening of the Molasse basin. Because the marine conditions were generally very shallow, subtle changes in eustatic sea level contributed to the formation of several hiatuses that chronicle periods of erosion and non-sedimentation. These processes also amplified the tectonically induced increase in accommodation space during times of global sea level highstands. Whereas these mechanisms are capable of explaining the establishment of the Burdigalian seaway and the formation of distinct sedimentological niches in the Swiss Molasse basin, the drainage reversal during OMM times possibly requires a change in tectonic processes at the slab scale, most likely including the entire Alpine range between the Eastern and Central Alps.
In conclusion, we consider rollback tectonics to be the main driving force controlling the transgression of the OMM in Switzerland, with contributions by the uplift of individual crustal blocks (here the Aar massif) and by a reduction of sediment supply. This reduction of sediment flux was likely to have been controlled by tectonic processes as well when basement blocks became uplifted, thereby modifying the catchment geometries. Eustatic changes in sea level explain the various hiatuses and amplified the deepening of the basin during eustatic highstand conditions.
Foreland basins and their deposits have often been used to explore the tectonic evolution of their hinterlands, mainly because these basins are mechanically coupled with the adjacent mountain belts (Beaumont, 1981; Jordan, 1981; DeCelles, 2004). The formation of these foreland basins occurs through the flexural downwarping of the underlying lithosphere in response to loading, which results in the formation of a wedge-shaped trough where sediment accumulates (DeCelles and Giles, 1996; Allen and Allen, 2005). The shape of the foreland trough depends on the mechanical properties of the foreland plate (Sinclair et al., 1991; Flemings and Jordan, 1990; Jordan and Flemings, 1991), the load of the sedimentary fill itself (Flemings and Jordan, 1990; Jordan and Flemings, 1991), and predominantly the tectonic and geodynamic processes leading to changes in plate loading (Beaumont, 1981; Jordan, 1981; Allen et al., 1991; Sinclair et al., 1991; DeCelles and Gilles, 1996). Additionally, a foreland basin can either be occupied by a peripheral sea when sediment flux is lower than the formation rate of accommodation space or by a fluvial system if the opposite is the case (Sinclair and Allen, 1992). A shift from terrestrial to marine conditions, for instance, can occur through a reduction in sediment flux, an increase in tectonically controlled subsidence rate (Sinclair et al., 1991), or a rise in the eustatic sea level (Reichenbacher et al., 2013). This has particularly been inferred for the North Alpine Foreland Basin, or the Molasse basin, situated on the northern side of the European Alps (Fig. 1a), which experienced a change from terrestrial to marine conditions during Burdigalian times ca. 20 Myr ago (e.g. Matter et al., 1980; Pfiffner, 1986; Schlunegger et al., 1997a; Kempf et al., 1999; Kuhlemann and Kempf, 2002; Ortner et al., 2011). In the Swiss part, this transgression resulted in the establishment of a shallow marine seaway, linking the Paratethys in the NE with the Tethys in the SW (Allen et al., 1985), and it is recorded by the deposition of the Upper Marine Molasse group (OMM; Matter et al., 1980; P. A. Allen, 1984). Although the history and geometry of the Burdigalian seaway and the related sedimentary processes are well known through detailed sedimentological and chronological investigations (e.g. Lemcke et al., 1953; P. A. Allen, 1984; Allen et al., 1985; Homewood et al., 1986; Doppler, 1989; Keller, 1989; Jin et al., 1995; Salvermoser, 1999; Strunck and Matter, 2002; Kuhlemann and Kempf, 2002; Reichenbacher et al., 2013), the controls on this transgression have still been a matter of ongoing scientific debate. Previous authors (e.g. Allen et al., 1985; Homewood et al., 1986; Keller, 1989; Strunck and Matter, 2002; Reichenbacher et al., 2013) proposed a combination of a reduced sediment flux and a rise in global sea level as possible mechanisms. However, thermo-chronometric data from the core of the Alps (Lepontine dome, Fig. 1; Boston et al., 2017) and structural work in the external Aar massif (Fig. 1; Herwegh et al., 2017) revealed that the Burdigalian was also the time of major tectonic events, including thrusting in the Aar massif and tectonic exhumation of the Lepontine dome (Schmid et al., 1996), leading to changes in surface loads in the Alps. Accordingly, it was proposed that tectonic processes could also have controlled the transgression of the peripheral sea through downwarping of the foreland plate (Sinclair et al., 1991). It is possible that the Swiss part of the Molasse basin bears key information to differentiate between these underlying mechanisms (eustasy, sediment flux, tectonics) because it is situated to the north of the Aar massif and the Lepontine dome (Fig. 1).
The aim of this paper is to disentangle tectonic processes, a reduction of sediment flux, and changes in the eustatic sea level as controls on the transgression of the OMM in the Swiss Molasse. We analysed OMM outcrops and sections along the entire Swiss Molasse basin at both proximal and distal positions relative to the Alpine front (Fig. 1a). To this extent, we measured palaeo-flow directions, and we explored the OMM deposits for their sedimentary facies and related depositional settings. We also determined the thickness and grain size of sedimentary bedforms, and we applied hydrological concepts to calculate palaeo-water depths based on the measured parameters. However, a chronological framework is absolutely required for correlating sections across a basin where facies relationships are strongly heterochronous, as is the case for the Swiss Molasse basin (Matter et al., 1980). Therefore, we also reassessed the temporal framework of the analysed sections through a compilation of previously published magnetostratigraphic and biostratigraphic data, and we correlated the individual sections across the basin.
The doubly vergent Alpine orogen (Fig. 1) is the consequence of the late Cretaceous to present continent–continent collision between the European and Adriatic plates (Schmid et al., 1996; Handy et al., 2010). It is comprised of a double vergent nappe stack wherein the material on the northern side was derived from the European and Adriatic plates, while the rocks on the southern side is only of Adriatic provenance. In the centre, the Alps expose a crystalline core of European origin referred to as the Lepontine dome and the external massifs (Fig. 1; e.g. Aar massif; Spicher, 1980). At deeper crustal levels, the Alpine orogen is underlain by a thick crustal root made up of a stack of lower crustal material derived from the European continental plate (Fry et al., 2010; Fig. 1b). Beneath the core of the orogen, the ca. 160 km long (Lippitsch et al., 2003) lithospheric mantle slab of the European continental plate bends and thus downwarps the foreland plate towards the SE (Fig. 1b). This bending was mainly driven by slab load due to the relatively large density of the subducted lithospheric mantle in comparison with the surrounding asthenosphere as seismo-tomography data reveal (Lippitsch et al., 2003). On the northern side of the Alps, the structurally highest unit is made up of Austroalpine nappes that structurally overlie the Penninicum, which in turn are underlain by the Helvetic thrust nappes (Fig. 1a). The front of the Helvetic and Penninicum units is referred to as the basal Alpine thrust (Fig. 1b). On the southern side, the Alps are made up of the south alpine thrust sheets that consist of crystalline basement rocks and sedimentary units of African origin. This fold-and-thrust belt is bordered to the south by the Po basin (Fig. 1a). The northern side of the Alps are separated from the southern side by the north-dipping Periadriatic line that accommodated most of the shortening during the Oligocene and the early Miocene by backthrusting and right-lateral slip (Schmid et al., 1996).
Recently, slab load has been considered the major driving force of the subduction history of the European plate and for the exhumation of crystalline rocks (Kissling and Schlunegger, 2018). This also concerns the exhumation of the Lepontine dome (Fig. 1), where normal faulting along the Simplon detachment fault resulted in rapid tectonic exhumation of the dome between late Oligocene and early Miocene times with a peak recorded by thermo-chronometric data at 20 Ma (Fig. 1b; Hurford, 1986; Mancktelow, 1985; Mancktelow and Grasemann, 1997; Schlunegger and Willett, 1999; Boston et al., 2017). This was also the time when the Aar massif, situated on the European continental plate (Fig. 1b), experienced a period of rapid vertical extrusion (Herwegh et al., 2017). Herwegh et al. (2017) and Kissling and Schlunegger (2018) proposed a mechanism referred to as rollback subduction to explain these observations. According to these authors, delamination of crustal material from the European mantle lithosphere along the Moho resulted in a stacking of buoyant lower crustal rocks beneath the Lepontine dome and the Aar massif, forming the crustal root (Fry et al., 2010; Fig. 1b). These processes are considered to have maintained isostatic equilibrium between the subducted lithospheric mantle and the crust and thus the elevated topography (Schlunegger and Kissling, 2015). Additionally, they most likely balanced, through the stacking of the crustal root (Fry et al., 2010, Fig. 1b), the rapid removal of upper crustal material in the Lepontine dome at ca. 20 Ma (Schlunegger and Willet, 1999; Boston et al., 2017). Delamination of crustal material has also been invoked to explain the contemporaneous rapid exhumation and northward thrusting of the Aar massif along steeply dipping thrusts (Herwegh et al., 2017). These processes were contemporaneous with (i) the reorganization of the drainage network of the Central Alps (Kuhlemann et al., 2001a; Schlunegger et al., 1998), (ii) the decrease in sediment flux to the basin, as revealed by sediment budgets (Kuhlemann, 2000; Kuhlemann et al., 2001a, b), and (iii) the Burdigalian transgression in the Swiss part of the Molasse basin. We will thus refer to these processes when discussing the controls on the Burdigalian transgression within a geodynamic framework.
The Molasse basin is approximately 700 km long and striking ENE–WSW from
France to Austria, where it broadens from
Reconstructions of the evolution of the Molasse basin (Fig. 2a) have been the focus of many research articles over the past years (e.g. Matter et al., 1980; Homewood and Allen, 1981; P. A. Allen, 1984; Keller, 1989; Schlunegger et al., 1996; Sinclair, 1997; Kempf et al., 1999; Kuhlemann and Kempf, 2002; Ortner et al., 2011; Reichenbacher et al., 2013). This has resulted in the general notion that the large-scale subsidence history of the Molasse basin was closely linked with the geodynamic evolution of the Alps (Sinclair et al., 1991; Kuhlemann and Kempf, 2002; Pfiffner et al., 2002; Ortner et al., 2011; Schlunegger and Kissling, 2015). The development of this basin as a foreland trough has been considered to commence with the closure of the Alpine Tethys in late Cretaceous times (Lihou and Allen, 1996; Schmid et al., 1996). This was the time when subduction of the European oceanic lithosphere with a large density beneath the Adriatic continental plate started. Large slab load forces resulted in a downwarping of the European foreland plate and the formation of a deep marine trough (Schmid et al., 1996), where sedimentation occurred by turbidites (Sinclair, 1997) on submarine fans (Allen et al., 1991; Sinclair, 1997; Lu et al., 2018; Reichenwallner, 2019). This first period of basin evolution has been referred to as the flysch stage in the literature (Fig. 2b; Sinclair and Allen, 1992). The situation changed at 35–32 Ma when the buoyant continental lithosphere of the European plate started to enter the subduction channel. Strong tension forces started to operate at the stretched margin of the European continental crust, particularly beneath the Central Alps (Schmid et al., 1996), with the result that the subducted oceanic lithosphere of the European plate broke off (Davies and von Blanckenburg, 1995). The consequence was a rebound of the European plate, a rise of the Central Swiss Alps, and an increase in sediment flux to the Swiss Molasse basin (Sinclair, 1997; Kuhlemann et al., 2001a, b; Willett, 2010; Garefalakis and Schlunegger, 2018), which became overfilled at ca. 30 Ma (Sinclair and Allen, 1992; Sinclair, 1997). The subsequent post-30 Ma stage of basin evolution has been referred to as the Molasse stage (Sinclair and Allen, 1992).
Molasse sedimentation occurred from ca. 30 Ma onward and was recorded by two large-scale, transgressive–regressive mega-sequences (Fig. 2b; e.g. Sinclair, 1997; Kempf et al., 1997, 1999; Kuhlemann and Kempf, 2002; Cederbom et al., 2004, 2011). These two mega-sequences consist of four lithostratigraphic groups. The first mega-sequence comprises the Lower Marine Molasse (UMM; Diem, 1986) and the Lower Freshwater Molasse (USM; Platt and Keller, 1992), and the second mega-sequence consists of the Upper Marine Molasse (OMM; Homewood et al., 1986) and the Upper Freshwater Molasse (OSM; Matter et al., 1980). Sedimentation in the Molasse basin continued up to ca. 10–5 Ma, when a phase of uplift during Pliocene times resulted in erosion and recycling of the previously deposited Molasse units (Mazurek et al., 2006; Cederbom et al., 2004, 2011). This erosion reached deeper stratigraphic levels in the western part of the Molasse basin than in the eastern segment (Baran et al., 2014) with the consequence that the OMM deposits are only fragmentarily preserved in the west (Fig. 2a).
Sediment dispersal changed during the Molasse stage of basin evolution. Prior to 20 Ma, during USM times, measurements of sediment transport directions (Kempf, 1998; Kempf et al., 1999) and sediment provenance analysis (Füchtbauer, 1964) revealed that the sedimentary material was transported to the east by braided to meandering streams. At that time, a coastline was situated near Munich, separating a deep marine trough farther east from a terrestrial environment to the west of Munich (Kuhlemann and Kempf, 2002). During OMM times, heavy mineral assemblages reveal that the Swiss Molasse basin operated as a closed sedimentary trough, where all supplied material was locally stored (Allen et al., 1985). During OSM times, from ca. 16.5 to 5 Ma, heavy mineral data imply that material with sources in the Hercynian basement north of Munich or the Bohemian massif was supplied to the Swiss Molasse (“Graupensandrinne”; Allen et al., 1985; Berger, 1996), suggesting that material transport occurred towards the west (Kuhlemann and Kempf, 2002). The details of the reversal of the drainage direction have not yet been elaborated, and related scenarios lack a database with palaeo-flow information, particularly for the OMM. The establishment of such a database, particularly the assignment of a more precise age to the drainage reversal, will be part of the scope of this article.
The Upper Marine Molasse (OMM) deposits, which are the focus of this study, mainly consist of a suite of shallow marine sandstones and mudstones that were deposited between ca. 20 and 17 Ma (Fig. 2b) in a ca. 70–80 km wide seaway (Allen and Homewood, 1984; Allen et al., 1985; Keller, 1989; Strunck and Matter, 2002). Close to the Alpine thrust front, the OMM successions are up to 900 m thick and thin to a few tens of metres towards the distal basin margin farther northwest.
Large streams with sources in the Central Alps supplied their material to
the Molasse basin, thereby forming megafans and conglomerate deposits with
diameters
West of the Napf, the OMM at the proximal basin border has also been categorized in two units (Fig. 3a) but with a different scheme and different names. These are the Sense Formation at the base (suite of sandstones with mudstone interbeds) and the Kalchstätten Formation at the top (alternation of sandstones and mudstones). Further up-section, these marine deposits (Strunck and Matter, 2002) grade into the fluvial conglomerates of the Guggershorn Formation and thus into the OSM (Strunck and Matter, 2002) (Fig. 3a). Magnetostratigraphic dating showed that the accumulation of the Guggershorn Formation occurred contemporaneously with marine (OMM) sedimentation at more distal sites. These conglomerates thus represent the deposits of a braided stream that shed its material to the OMM sea, similar to the Napf conglomerates.
In the central basin near Fribourg (Fig. 2a), also on the western side of the Napf, sedimentological investigations of sand waves (Homewood and Allen, 1981; Homewood et al., 1986) disclosed the occurrence of tidal bundles and a bidirectional dispersal of sedimentary material. These deposits have been assigned to a subtidal environment (Homewood and Allen, 1981) where material with sources in the Central Alps was redistributed in the basin by strong tidal currents that entered the basin from the Tethys in the south and the Paratethys in the northeast (Allen et al., 1985; Kuhlemann and Kempf, 2002; Bieg et al., 2008). Lithostratigraphic correlations suggest that the deposits at Fribourg most likely correspond to the lower Sense Fm (Python, 1996) and thus to the OMM-Ia (see the Discussion section and Fig. 3a).
In the basin axis, between the Lake Neuchâtel and Wohlen areas (Fig. 2a), coarse-grained sandstones with large-scale cross-beds in which individual grains are larger than 2 mm have been interpreted as subtidal sand waves (Allen and Homewood, 1984; Allen et al., 1985). These deposits are calcareous sandstones with shelly fragments, referred to as the “Muschelsandstein” (Allen and Homewood, 1984; Allen et al., 1985, Fig. 3a). Alternatively, they occur as coarse-grained cross-bedded sandstones with large lithoclasts, also called the “Grobsandstein” (Jost et al., 2016).
Ages for the OMM deposits have been established by multiple authors through
palaeontological analyses of mammalian fragments and teeth (Keller, 1989;
Schlunegger et al., 1996; Kempf et al., 1999) as well as
Following the scope of the paper, we established the facies relationships
and sediment transport patterns during the transgressive phase of the OMM
and thus mainly focused on the OMM-I. We proceeded through sedimentological
investigations of key sites (Fig. 2a), which expose the related succession
in proximal (Entlen section east of the Napf units; Sense section west of
the Napf units), central (St. Magdalena site, Gurten drill core), and
distal positions (Lake Neuchâtel and Wohlen areas). The lithofacies in
the Entlen and Sense sections was investigated in the field at the scale of
Sedimentological logs of
Lithofacies encountered in the Entlen section.
Lithofacies encountered in the Napf units.
The lithofacies have been identified (e.g. Schaad et al., 1992; Keller, 1989) based on the assemblage of sedimentary characteristics, including grain size, thickness, lateral extent if applicable, sedimentary structures, basal contact, colour, and fossil content (Tables 1 to 5). The lithofacies types correspond to individual bedforms (see Tables 1 to 5 for references), which bear information on flow strengths, flow directions, sediment supply, and water depths (e.g. Keller, 1989). The combination of these parameters, usually recorded by distinct assemblages of lithofacies types, can be used to identify distinct sedimentary settings. Related concepts of facies analysis have been documented for fluvial deposits (Miall, 1978, 1985, 1996; Platt and Keller, 1992) but are less standardized for shallow marine deposits. Here, we followed Keller (1989) and Schaad et al. (1992), who developed a concept for shallow marine deposits whereby lithofacies types are grouped into facies assemblages in a hierarchic order, based on which distinct shallow marine settings can be interpreted. We followed these authors and assembled the various lithofacies types into five depositional settings, which are from land to sea: terrestrial, backshore, foreshore, nearshore, and offshore. We then mapped the depositional settings at the scale of 1 : 25 000 at various sites across the Swiss Molasse basin where suitable outcrops were present (see Fig. 2a for visited sites).
Sediment transport directions were determined from orientations of clast imbrications, gutter casts, and dip directions of cross-beds. In addition, the orientation of the coastline can be inferred from sediment transport within the surf-and-swash zone at the wet beach where rolling grains carve millimetre-thin rills in the beach deposits, which are oriented perpendicular to the coast. These rills are recorded by linear grooves, or parting lineations, on the surface of sandstones (J. R. L. Allen, 1982; Hammer, 1984). We thus measured the orientation of these features where visible. We also determined the strike direction of oscillation ripple marks to infer the orientation of waves and thus of the coastline.
We analysed the sediments in the key sections according to their palaeo-water depths. For oscillation ripple marks, the ripple metrics (spacing between ripple crests and ripple heights) together with the grain size can be used to infer water depths at the time the oscillation ripples were formed (Diem, 1985; Allen et al., 1985; see the Supplement). We thus measured the ripple metrics with a metre stick together with grain sizes in the field and calculated the water depths following Allen (1997). Likewise, minimum water depths can be inferred from the heights of cross-beds as examples from modern streams have shown (Bridge and Tye, 2000; Leclair and Bridge, 2001). Please refer to the Supplement for the deviation of the related equations. Published information from deep drillings (Boswil 1; Hünenberg 1; Schlunegger et al., 1997a) and seismostratigraphic data (line 8307, Schlunegger et al., 1997a; line BEAGBE.N780025; Fig. S2 in the Supplement) completed the available database.
The sedimentary suite of the ca. 370 metre thick OMM-Ia at Entlen (Figs. 4a and
S4a in the Supplement) records a large diversity of lithofacies types
(Table 1). Parallel-laminated (Sp), fine- to medium-grained sandstone
packages are centimetres to decimetres thick and normally graded. These deposits alternate
with decimetre-thick low-angle cross-bedded units with tangential lower boundaries
(Sc, Sct
The overlying ca. 430 m thick OMM-Ib (Figs. 4a and S4a in the Supplement)
comprises fine- to medium-grained sandstone packages with mudstone interbeds
(Table 1). Low-angle trough (Sct
Estimates of palaeo-water depths (see the Supplement) reveal that the OMM-Ia
sedimentary rocks were deposited in shallow conditions
Measurements of the bedform orientations of the OMM-Ia deposits reveal sediment
transport directions between 315
The Napf units, which are a terrestrial interval of the OMM and the OSM (Fig. 3; Schlunegger et al., 1996), are ca. 1550 m thick and include a succession of conglomerates, sandstones, and mudstone interbeds (Matter, 1964), which we categorize into five lithofacies types (Table 2, Fig. S4d in the Supplement). Individual conglomerate beds are up to 10 m thick and display stacks of 2–3 m thick beds with massive- (Gm) to cross-bedded (Gc) geometries. The sandstone beds occur as massive-bedded (Sm) and cross-bedded (Sc) units. Interbedded mudstones are horizontally bedded (Mp) and have a yellowish–reddish mottling, caliche nodules, and root casts. Palaeo-flow measurements imply a change from a NE-directed transport during USM times to a NW-directed sediment transport between OMM-I and OSM times.
The OMM at Sense (Figs. 4b and S4b) starts with a ca. 200 m thick
succession of predominantly sandstones with some mudstone interbeds. Medium-
to coarse-grained sandstone beds, up to 2–3 m thick, are massive-bedded
(Sm), parallel-laminated (Sp), and trough cross-bedded (Sct
Lithofacies encountered in the Sense section.
The first occurrence of 5–10 m thick sandstone beds at the 200 m
stratigraphic level (Fig. 4b) marks a distinct shift in the stratigraphic
record where several metre-thick cross-bedded sandstone beds dominate the
sedimentary succession. At this level, (i) 5–10 m thick normally graded
sandstone beds overlie an erosive base and display epsilon cross-beds (Sce);
(ii) cross-bedded sandstones (Sct
Estimates of palaeo-water depths range between 5 and 10 m (Fig. 4b and
Table S4 in the Supplement) during the deposition of the lowermost 200 m.
Conditions were deepest at the 200 m stratigraphic level, reaching water
depths in the range of up to 30 m (see the Supplement). Measurements of sediment
transport directions cover the range between ca. 0
The sandstones within a cave system near Fribourg (St. Magdalena; Figs. 2a
and S4c in the Supplement) are medium- to coarse-grained and display an
amalgamation of up to 1–3 m wide cross-bedded troughs (Sct
In the nearby ca. 260 m deep Gurten (Fig. 2a) drill core (Fig. S3 in
the Supplement), OMM-Ia deposits occur as cross-bedded sandstones (Sc) topped
with mudstone drapes (Md). These lithofacies associations (Table 4) are most
abundant within the drill core and make up ca. 200 m of the log. However,
because drill cores offer limited information about the dimensions of the
encountered sediments, we were not able to determine if cross-beds can be
assigned to tabular beds (Sct
Lithofacies encountered at the St. Magdalena site and Gurten drill core.
Calcareous, shelly sandstones (Scc) occur in the basin axis and are an assemblage of various lithofacies. This Scc facies association is made up of 5–10 m thick, coarse-grained sandstone beds with low-angle cross-beds (Sc) that contain coquinas, shell fragments (Shf), and pebbles (Sg) in places. Interbedded fine-grained sandstones contain current ripple marks (Scr) recording an opposite flow direction relative to the cross-beds (Sc).
In the west (sites at Lake Neuchâtel area; Figs. 2a and S4c in
the Supplement), mapping shows that Scc Muschelsandstein deposits are ca. 5 m
thick and record NNE- to NE-directed sediment transport. Foreset thicknesses
of these deposits thin to
The OMM-Ia sedimentary rocks of the Entlen section are assigned to a backshore to upper nearshore realm within a wave-dominated environment (Fig. 4a, and please see Table 1 for references). Records of waves are inferred from (i) tabular, parallel-laminated and normally graded (Sp) sandstones, which are interpreted to represent sediments of the surf-and-swash zone near the wet beach where sedimentation occurs in the upper flow regime; (ii) low-angle cross-beds (Sc) with pebbly lags and shell fragments (Shf), which could reflect sand reefs (or shoals), rip channel fills, or storm layers; and (iii) oscillation (Sos) as well as branching ripple marks (Sbr) pointing to wave activity (Fig. 4a). Gravels and pebbly lags (Sg) are either evidence of high-energy storm events or river inflow from the backshore. Finer-grained lithofacies, which are either indicative of rapid sedimentation (Sv, Mcl; Table 1) or incipient pedogenesis (Mp, Mm; Table 1), are consistent with a shallow marine, wave-dominated environment.
The basal part of the OMM-Ib suite is assigned to a foreshore to lower
nearshore setting shaped by the combined effect of wave and tidal activity
(Fig. 4a). This is inferred by the observation that current ripple marks
(Scr), which are situated on top of lamina sets of tabular (Sct
The change from a nearshore, wave-dominated environment (OMM-Ia) to an
environment with tidal records (OMM-Ib) was additionally associated with a
deepening from
We interpret the association of massive- (Gm) to cross-bedded (Gc) conglomerates and massive- to cross-bedded sandstones (Sm, Sc) as deposits within a braided river system (Table 2; please see references there). In such an environment, conglomerates are common records of active channels. Massive- (Sm) to cross-bedded (Sc) sandstones alternating with mottled mudstones were most likely formed on the floodplains bordering the network of braided channels when bursts resulted in the accumulation of crevasse splay deposits (Platt and Keller, 1992). Mudstone interbeds (Mp) with evidence of palaeosol genesis formed when channel belts shifted away from the axis of the section (Platt and Keller, 1992). This facies association was mapped over tens of kilometres, both across and along strike of the basin orientation. It is thus assigned to an alluvial megafan (Schlunegger and Kissling, 2015), which was deposited by braided streams.
We assign the OMM of the Sense section to a tidal-dominated environment
in which deltaic estuaries dominated the sedimentary facies (Fig. 4b, and please
see Table 3 for references). This is inferred from (i) sigmoidal
cross-bedded sandstones (Sc) with distinct top sets, foresets, bottom sets, and
pebbly lags (Sg) that are indicative of a delta and (ii) trough
cross-bedded (Sct
The first occurrence of 5–10 m thick sandstone beds at the 200 m
stratigraphic level (Fig. 4b) records a remarkable increase in the water
depth when 5–10 m deep tidal channels (Sce sandstone beds with epsilon
cross-beds) grade into several metre-thick subtidal sand waves (Sct
The several metre-thick outcrops near Fribourg are interpreted as subtidal shoal
deposits, which accumulated within a tidal-dominated environment (Table 4;
please see references there). This is inferred from metre-scale
cross-bedded troughs (Sct
We interpret the Scc Muschelsandstein (Table 5) sediments (containing
coquinas, shell fragments (Shf), and pebbles – Sg) to have been deposited
within the topographic axis of the Burdigalian seaway (see also Allen et
al., 1985, and Jost et al., 2016, for a similar interpretation). Metre-scale
cross-beds (Sct
Lithofacies encountered in the Lake Neuchâtel and Wohlen areas.
No magnetostratigraphic data are available for the Entlen section, but a temporal calibration of the deposits can be achieved through indirect lines of evidence. This particularly concerns the reconstruction of an age constraint for the basal transgression, which is accomplished using two lines of evidence: first, the transgression post-dates the deposition of the USM, which terminated at C6An1 at the Fischenbach section (Fig. 3b; Schlunegger et al., 1996). Second, based on stratigraphic interpretations of palaeo-flow direction data, Strunck and Matter (2002) considered the transgression of the OMM to have progressed from the east towards the west, where the first marine sediments have been dated with C6r in the Sense section (Fig. 3b; see next section). An E–W transgression of the OMM is also seen in seismic line BEAGBE.N780025 (Fig. S2 in the Supplement), where OMM deposits onlap the USM in a westward direction. Accordingly, the onset of the OMM at Entlen on the eastern side of the Napf predates the transgression at Sense farther west. Based on these arguments, we set an age of ca. 20 Ma for the base of the OMM-I in the eastern Swiss Molasse basin (Figs. 3b and 5), which is consistent with Kälin and Kempf (2009). For the top of the OMM-I, we determine an age using the magnetostratigraphy of the Napf section (Fig. 3b) ca. 10 km to the west of Entlen. This section includes an alternation of six reversed- and five normal-polarized magnetozones (Schlunegger et al., 1996). The lowermost, very long normally polarized interval (N1, Fig. 3b) includes the mammalian fossil site Hasenbach 1, recording an MN3b age (Schlunegger et al., 1996) or a lower MN3b age, as a revision of the mammalian material has shown (Kälin and Kempf, 2009). This allows for a correlation of the normally polarized interval N1 with chron 5En of the MPTS (Cande and Kent, 1992, 1995) and the ATNTS (Lourens et al., 2004), respectively (Fig. 3b). Since the third reversed magnetozone of the Napf section is very short (R2), and since the ATNTS chron 5D spans several hundred thousand years and is thus quite long, it is most likely that a hiatus encloses C5Dr2 to C5Dr1 (Fig. 3b). In addition, because (i) the change from MN3b to MN4a has been calibrated with C5Dr2 (Jost et al., 2016), and since (ii) the base of the overlying OMM-II (Figs. 3b and 5) has been dated with MN4 (Keller, 1989), we suggest that the inferred hiatus coincides with the boundary between the OMM-I and the OMM-II (Figs. 3b and 5). This age assignment is consistent with magneto-polarity stratigraphy in the Molasse basin ca. 70 km farther to the east (Kempf and Matter, 1999). It is also consistent with micro-mammalian investigations in the distal Molasse basin ca. 30 km farther north where Jost et al. (2016) found that deposits spanning sometime between late MN3b and MN4a are missing. Based on these constraints, we suggest that the top of the OMM-I correlates with C5Dr of the MPTS and C5Dr2 of the ATNTS, respectively, followed by a ca. 0.5 Myr long hiatus (Figs. 3b and 5). According to this correlation, the sediments recording the maximum flooding conditions in the Entlen section are ca. 19 Myr old.
Correlations of the OMM deposits from the Entlen section with the Wohlen area were accomplished by Schlunegger et al. (1997a) through a seismostratigraphic analysis of the seismic line 8307 (please see Fig. 2a for a trace of the line). The seismic data show that the OMM-I deposits onlap USM strata and then overlap this unit (Fig. 5b). Schlunegger et al. (1997a) correlated the OMM-I sequence with their Unit B in the Entlen section, which corresponds to the top of the OMM-Ia in our stratigraphic scheme. In addition, our field investigations and micro-mammalian data by Jost et al. (2016) revealed that the Muschelsandstein follows on top of the OMM-Ia and most likely corresponds in age to the OMM-Ib. Based on these arguments, we constrain the deposition in the distal basin to the time interval between ca. 19 and 18 Ma (Fig. 5b). However, based on seismostratigraphic investigations of line 8307, Schlunegger et al. (1997a) proposed that sedimentation was interrupted at ca. 18 Ma by a ca. 0.5 Myr long or possibly longer hiatus. This time span was later specified through new micro-mammalian discoveries by Jost et al. (2016), who noted that a record of MN4a is missing in the Wohlen area and that the base of the OMM-II hosts mammalian fragments that correspond to MN4b. The interpretation of an inferred unconformity is additionally supported through observations of vadose cements (Allen et al., 1985) within the Muschelsandstein and through evidence of a thick palaeosol separating OMM-I from OMM-II in the Entlen section (see Sects. 4 and 5). We use the occurrence of vadose cements and the palaeosol to propose that the uppermost beds of the OMM-Ib (including the Muschelsandstein) were exposed to erosion, or non-sedimentation, after deposition. Furthermore, because the Muschelsandstein unit records the deepest water depth during OMM-I times at distal sites, we tentatively suggest that the deposition of these mega-sand waves started at the same time when the deepest conditions (MFS) were recorded within the Entlen section (Fig. 5b).
Magnetostratigraphic data for the Sense section were presented by Strunck and Matter (2002). These authors placed the USM–OMM boundary at this site within C6r of the MPTS (Cande and Kent, 1992, 1995) or alternatively of the ATNTS (Fig. 3b; Lourens et al., 2004). The subsequent alternation of normal and reverse magnetozones was correlated by these authors with chrons 6r through 5Dn of Cande and Kent's MPTS (1992, 1995), the latter of which corresponds to C5Dn1 of the ATNTS (Lourens et al., 2004). Following Strunck and Matter (2002), a possible hiatus prior to ca. 17.7 Ma is also likely to be recorded within the Sense section (Figs. 3b and 5a). This correlation implies that the lower Sense Fm corresponds to the OMM-Ia, whereas the upper Sense Fm and the Kalchstätten Fm are time-equivalent units of the OMM-Ib (Fig. 3a and b). The topmost 50 m of the Kalchstätten Fm follows upon this hiatus and corresponds to the OMM-II in our scheme (Figs. 3 and 5a). This further implies that a hiatus separates OMM-I from OMM-II across the entire basin (Fig. 5a). In addition, subsequent sedimentation (base of OMM-II) progressed from the west to the east (Fig. 5a). In the same sense, the sediment packages recording the maximum flooding conditions (MFS) most likely have the same age across the entire basin (Fig. 5a).
No micro-mammalian sites have been reported for the OMM deposits in the
distal western Molasse basin. Therefore, we cannot provide further
constraints on the history of sedimentation. However, our field inspections
in the area of Lake Neuchâtel (Fig. 2a) show a sedimentary succession
similar to that in the east, where amalgamated sandstone beds are overlain
by the Muschelsandstein unit. Our field inspections also show that these
calcareous, shelly sandstones thin from ca.
The chronostratigraphic framework together with the sedimentological data and palaeo-flow directions are used to propose a scenario of how the basin evolved through time. During USM times (Fig. 6a), prior to the Burdigalian transgression, the basin was occupied by alluvial megafans at the proximal basin border, which gave way to an axially directed channel-belt system in the distal basin (Fig. 6a; Kuhlemann and Kempf, 2002). Analysis of heavy mineral assemblages by Füchtbauer (1964) and measurements of palaeo-flow directions in our study area (sole casts and cross-beds; this paper) and in eastern Switzerland (Kempf et al., 1999) revealed a NE-directed material transport towards the Munich region (Fig. 6a), which is consistent with the results of previous syntheses (e.g. Pfiffner et al., 2002; Kuhlemann and Kempf, 2002). In this area, the Molasse streams ended in a peripheral sea where neritic to open marine conditions prevailed (Kuhlemann and Kempf, 2002). Within the basin, a possible divide for sediment transport was situated somewhere SW of Geneva. We infer such a separation of sediment dispersal based on published sediment transport directions. In particular, south of Geneva, tidal cross-beds imply a sediment dispersal towards the south and thus towards the Tethys (Allen et al., 1991; Allen and Baas, 1993). To the northeast of Geneva, however, our own measurements and data from Kempf et al. (1999) reveal a NE-directed sediment transport to the Paratethys.
Palaeogeographical reconstructions of the Molasse basin
at different stages:
The palaeogeographic reconstruction based on our and published data for the period between ca. 20 and 19 Ma is shown in Fig. 6b. It illustrates that the central part of the Molasse basin changed to a shallow marine sea, which was ca. 40 km wide at that time. Our estimates of palaeo-bathymetric conditions and sedimentological data (Fig. 4a and b) reveal that the water depths corresponded to a subtidal and nearshore setting. Nearshore to possibly offshore conditions (30–50 m) are recorded by subtidal mega-sand waves (Allen and Bass, 1993) south of Geneva and by the predominance of sandstone–mudstone alternations within the Boswil and Hünenberg drill cores in the NE (Wohlen area, Fig. 2a; Schlunegger et al., 1997a). Subtidal shoals, in up to 5 m deep water, occupied the western part of the central Swiss Molasse (Fig. 6b). This was already proposed by Homewood and Allen (1981), and it is confirmed here by our sedimentological data and estimates of palaeo-water depths (see the Supplement). Measurements of sediment transport directions from the shoal deposits reveal bimodal, SW–NE-directed transport with a dominant NE orientation. This is particularly the case at the proximal basin border near the Sense section (Fig. 2a) where deltaic foresets accumulated within an estuarine setting (Fig. 4b). Mapping of depositional settings allowed us to trace the shoal deposits towards the northern tip of the Napf megafan, from which the shoals narrow from ca. 20 to ca. 10 km over a 70 km long distance along strike. It thus appears that the shoals were deflected towards the topographic axis through a dominant NE-directed material transport (Fig. 6b). This interpretation is additionally supported by measurements of the transport directions of the Napf megafan (i.e. clast imbrications) and the coastal deposits at the Entlen section (i.e. parting lineation, cross-beds) pointing material transport towards the NE (Figs. 4a and 5a). At the distal margin of the basin, field inspections show that beach sandstones gave way to subtidal shoal deposits up-section. It thus appears that the Molasse basin between the Lake Neuchâtel and Wohlen areas (Fig. 2a) was a region of sediment export to the NE and to the SW. Material transport was most likely accomplished through strong tidal currents that entered the Swiss Molasse as two major tidal waves from the Tethys in the south and the Paratethys in the northeast (Bieg et al., 2008).
The situation during ca. 19–18 Ma (Fig. 6c) started with the time when the
maximum flooding surface was formed in the depositional record (MFS;
Fig. 5a and b). The sedimentological data reveal that this time was
characterized by a widening of the basin to widths up to 80 km, and it was
dominated by offshore conditions in the topographic axis with water depths
The palaeogeographic situation shown in Fig. 6d comprises the time span between ca. 18 and ca. 14 Ma and displays the evolution from the OMM-II to the OSM. The OMM-II period followed a phase of non-sedimentation and possibly erosion across the entire Swiss Molasse basin, as our reassessment of the chronological framework of the OMM reveals (Fig. 5a and b). In addition, measurements of sediment transport directions reveal a SW-oriented sediment transport at proximal positions (Fig. 5a), which is consistent with the results of previous syntheses (e.g. Pfiffner et al., 2002; Kuhlemann and Kempf, 2002). This also implies that a possible E–W divide for sediment transport shifted towards the region near Munich or even farther east. Similar to Kuhlemann and Kempf (2002), we infer such a scenario from the supply of material with sources in the Hercynian basement north of Munich (Fig. 1a) or the Bohemian massif (Graupensandrinne, Fig. 6d; Allen et al., 1985; Berger, 1996; see also Sect. 2), which implies a westward tilt of the basin axis. This period ended with the progradation of the alluvial megafans during the time of the OSM.
In conclusion, this study confirms the results and syntheses of previous authors on the general sedimentation and material transport pattern during the deposition of the OMM. Nevertheless, our refinement of the chronological framework in combination with additional sediment transport data allow us to specify some further details on the development of the transgression of the OMM. These include (i) the establishment that the Burdigalian seaway was accompanied by both a deepening and widening of the basin. These mechanisms occurred contemporaneously and were associated with a northward shift of the topographic axis to the distal basin margin, where offshore and thus the deepest marine conditions established at 19 Ma. (ii) The reversal of the sediment transport direction from an originally NE-oriented sediment dispersal to a SW-oriented sediment transport started sometime after 20 Ma and was completed at 18 Ma at the latest. (iii) A wave-dominated coastline (with some tidal records) established on the eastern side of the Napf megafan (Fig. 4a), whereas a tidal-dominated estuarine environment characterized the proximal coastal margin on the western side of the Napf (Fig. 4b). These variations in sedimentation pattern appear to explain why the lithostratigraphic framework differs between the two regions (see Fig. 3a).
We relate the reversal of the drainage direction between the OMM-I and the OMM-II to tectonic processes operating at deeper crustal levels beneath the Alps. This interpretation is guided by Pfiffner et al. (2002), who related changes in sediment dispersal within the Swiss part of the basin to a possible tilt of the foreland plate caused by the westward shift of the Ivrea body. This tectonic unit is comprised of mantle rocks with a high density (Fig. 1a) and thus could have influenced the deflection of the foreland plate (Pfiffner et al., 2002). While this mechanism is a viable explanation for the westward-directed tilt of the basin axis, we argue that a complementary driving force beneath the Eastern Alps is required to explain the drainage reversal across the entire basin, at least between Germany and Switzerland. We thus present a hypothesis of a possible geodynamic scenario to explain the 18 Myr old change in the drainage direction in the next section, but we also acknowledge that this interpretation is speculative at this stage and warrants further investigations. Such an exploration, however, requires the geodynamic processes between 33 and 30 Ma to also be considered. At that time, the subducted European oceanic lithosphere was considered to have broken off beneath the Central Alps (Davies and von Blanckenburg, 1995; Schmid et al., 1996). However, underneath the Eastern Alps the European oceanic lithosphere remained attached to the continental plate as palinspastic restorations revealed (Handy et al., 2015). The consequence of slab break-off beneath the Central Alps was the rise of the Alpine topography and a large sediment flux to the Swiss Molasse basin (Sinclair, 1997; Kuhlemann et al., 2001a, b; Willett, 2010; Garefalakis and Schlunegger, 2018), which became overfilled at ca. 30 Ma (Sinclair and Allen, 1992; Sinclair, 1997; see also Sect. 2.1). East of Munich, however, the basin still remained underfilled until ca. 20 Ma as evidenced by deep marine sedimentation, whereby debris flows and proximal turbidites accumulated within the basin axis (Fertig et al., 1991; Malzer et al., 1993). We use these observations to propose that vertically directed slab load forces were still downwarping the foreland plate beneath the Eastern Alps to allow such a deep trough to form. In contrast, slab break-off beneath the Central Alps most likely caused a rebound of the foreland plate in Switzerland (Schmid et al., 1996; Schlunegger and Castelltort, 2016). We interpret the consequence to be a stronger downward deflection of the European foreland plate beneath the Eastern Alps compared to the Central Alps, which could explain the east-directed sediment transport prior to ca. 20 Ma (Fig. 6a).
Between ca. 20–17 Ma, i.e. during OMM times, a remarkable change was recorded in the Molasse basin. The eastern Molasse basin experienced a change from deep to shallow marine conditions (Kuhlemann and Kempf, 2002), and the entire basin recorded a reversal of the drainage direction from the E to the W (see above). We relate these shifts to a change in the pattern of slab load forces underneath the Eastern and Central Alps. Particularly in the eastern Molasse basin, the change from deep to shallow marine conditions could reflect a response to slab unloading through delamination, or break-off, of the subducted European lithosphere underneath the Eastern Alps (Ustaszewski et al., 2008), while rollback subduction of the European plate beneath the Central Alps of Switzerland continued, as Kissling and Schlunegger (2018) proposed. This could have resulted in a rebound of the European plate beneath the Eastern Alps, whereas plate downwarping continued beneath the Central Alps. We interpret these along-strike differences in the plate deflection to have caused a westward tilt of the foreland plate, which in turn could have controlled the drainage reversal. The reasons for the along-strike differences in the subduction mechanisms are not clear at this stage and could either be related to (i) inheritance related to the Mesozoic phase of rifting (Schmid et al., 2004; Handy et al., 2010) or (ii) differences in the mechanical strengths and rheological conditions of the foreland plate between the Swiss and the German–Austrian Molasse basins (Tesauro et al., 2009, 2013).
The drainage reversal occurred simultaneously with the widening of the basin in central Switzerland, as our chronological refinement shows. Because both changes occurred at the same time, one could infer a causality of the underlying controls. We refrain from such a view at this stage since we lack a detailed 3-D restoration of the tectonic and geodynamic situation of the Alps–Molasse basin system for that time. Nevertheless, a possible control on the widening of the basin can be identified based on a cross-sectional view across the Central Alps and the related geodynamic processes during the Burdigalian ca. 20 Myr ago (Fig. 7). At that time, we propose that the velocity of rollback subduction beneath the Central Alps of Switzerland was likely to have accelerated. We justify this interpretation through the observation that (i) tectonic exhumation of the Lepontine dome (Fig. 1), accomplished through slip along the Simplon detachment fault (Mancktelow and Grasemann, 1997), occurred at the highest rates at that time (Boston et al., 2017; Schlunegger and Willet, 1999) and that (ii) the rapid rise of the Aar massif (Fig. 1) also commenced at 20 Ma (Herwegh et al., 2017). Following the concepts of Kissling and Schlunegger (2018), these processes require a mechanism whereby buoyant crustal material several tens of kilometres thick was delaminated from the subducting European continental plate and accumulated within the crustal root (Fig. 7) within a short time period. In agreement with Kissling and Schlunegger (2018), we interpret a rapid phase of rollback subduction as also being able to shift the basin axis to more distal sites. We use these mechanisms to explain the cross-sectional widening of the Molasse basin at 19 Ma and the northward shift of the basin axis (Fig. 7), thereby giving way to the deposition of the offshore Muschelsandstein (Fig. 6b and c). This lithofacies association has been mapped along the distal basin border adjacent to the external massifs only.
Simplified geological–geophysical model of the Alpine orogen for the time between 20 and 18 Ma, showing the most important geodynamic forces that might have shaped the Molasse basin and induced the transgression. Modified after Schlunegger and Kissling (2015) and Kissling and Schlunegger (2018).
We interpret the different coastal morphologies between the eastern and
western sides of the Napf megafan (see Sect. 5) to have rather been controlled by
tectonic processes than by contrasts between the tidal waves from the
Paratethys and the Tethys. We tentatively exclude a surface control at this
stage because the pattern of tidal sand waves in the basin axis is not
reflected by a corresponding coastal morphology at the proximal basin
margin. In particular, the Muschelsandstein cross-beds are thicker in
the Wohlen area (up to 10 m) than in the Neuchâtel region (
Structural mapping in the Aar massif (Wehrens, 2015; Wehrens et al., 2017) has revealed that crustal blocks were rising along steeply SE-dipping thrust faults (Fig. 1b). This process was related to rollback subduction of the European mantle lithosphere and the related delamination of crustal material, which resulted in the rise of the Aar massif (Herwegh et al., 2017; Fig. 1b). This mechanism also lifted the topography surrounding the massif to higher elevations (Fig. 7), thereby forming a positive anomaly in the topographic load in the region. Sinclair et al. (1991) and Sinclair (1996) explored a possible stratigraphic response to topographic loading associated with the uplift of the Aar massif (Sinclair et al., 1991) through the application of a linear elastic plate model where thrusting and erosion are dynamically coupled. In their model, the distance between the location of thrusting (Aar massif) and the site in the basin wherein a signal is expected depends primarily on the flexural rigidity (or alternatively the elastic thickness or the Te value) of the crustal rocks underlying the foreland basin (Sinclair, 1996). The flexural rigidity of the rocks beneath the Swiss Molasse basin has been quantified with an elastic thickness of ca. 10 km using stratigraphic constraints (Sinclair et al., 1991). This estimate is particularly based on thickness gradients of accumulated Molasse deposits across a section from the distal basin border to the Alpine thrust front. This pattern, however, could have been influenced by upper crustal inhomogeneities (Waschbusch and Royden, 1992) such as pre-existing faults (Pfiffner, 1986). This could explain why estimates of Te values that are based on stratigraphic data are lower (Sinclair et al., 1991; Schlunegger et al., 1997b) than estimates that are based on the curvature of the entire European foreland plate from the distal Molasse border to the core of the Alps and even deeper (Pfiffner et al., 2002; Schlunegger and Kissling, 2015), where Te values up to 50 km have been proposed (see the discussion in Pfiffner et al., 2002). However, because the flexural response of the Molasse basin to local topographic loads was likely to have been influenced by inherited faults in the basement (Pfiffner, 1986), lower Te values appear more appropriate (Waschbusch and Royden, 1992). Accordingly, if we consider a local and thus upper crustal response to loading characterized by a Te value of ca. 10 km (Sinclair et al., 1991; Schlunegger et al., 1997b), then shifts in surface loads through the kilometre-thick stacking of additional material in the Aar massif are likely to have resulted in the formation of several tens of metres of supplementary accommodation space at the proximal basin border, as the models of Sinclair et al. (1991) predict. As a result, depocentres in the Molasse basin are predicted to step back to proximal positions (Sinclair et al., 1991), which is consistent with our interpretation of the flow directions and the inferred backstepping of the Napf megafan (see Sect. 6.1 and Fig. 6). In addition, according to Sinclair (1996), upward-directed bulging of a few tens of metres is expected at the distal (forebulge) and at the lateral margins of the load (lateral bulge). The spacing between an expected lateral bulge and the location of the surface forcing ranges between 50 and 100 km (using a Te value of 10 km), which is consistent with the distance between the Aar massif and the inferred subtidal shoals in the western Swiss Molasse basin (near Fribourg; Figs. 2a and 6b). Accordingly, we suggest that the establishment of subtidal shoals at the northern tip and on the western side of the Napf is the consequence of this bulging (Fig. 6b and c). Because the plate had an eastward tilt at that time, as inferred from palaeo-flow directions, such a flexural signal could possibly not be recorded on the eastern side of the Aar massif where the marine conditions were too deep (Fig. 6b and c).
We additionally use these mechanisms to explain the development of different
depositional settings at the proximal basin border of the Swiss Molasse
basin. East of the Napf megafan, a relatively high subsidence (rate of
ca. 340 m Ma
The time around 20 Ma was also characterized by a continuous reduction in
sediment flux from originally 25 000 km
Molasse stages (USM, OMM, and OSM) with mean palaeo-transport directions (black arrows), hiatuses plotted against stable oxygen isotope stages (Miller et al., 1996), and sediment flux (Kuhlemann, 2000; Willet, 2010). The red triangle demarcates the onset of delamination and rapid exhumation of the Aar massif (Herwegh et al., 2017). Grey arrows demarcate decreases in sediment flux and falls in the eustatic sea level, possibly contributing to the related hiatuses. MFS: maximum flooding surface.
Whereas tectonic processes are recorded in the arrangement of depositional
settings in the entire Swiss Molasse basin, signals related to the eustatic
changes in sea level are possibly recorded by several hiatuses. This
particularly concerns the times of non-sedimentation between OMM-I and
OMM-II and between the OMM and the OSM phase, which we have elaborated on in
Sect. 5.1 (Figs. 5a and 8). In this context,
In summary, we suggest that the Burdigalian transgression was related to a combination of a deepening and widening of the basin and a reduction of sediment supply rates, which we ultimately relate to tectonic processes in the Alpine hinterland. In this context, we consider rollback subduction to have most likely been responsible for the widening of the basin in the foreland and for the shift of the basin axis to distal positions. In addition, rollback subduction of the European mantle lithosphere and delamination of crustal material most likely resulted in the rapid exhumation of the Lepontine dome (Boston et al., 2017) and the associated rise of the Aar massif (Herwegh et al., 2017). These processes are interpreted to have triggered the change in the configuration of the drainage network (Schlunegger et al., 2001; Kühni and Pfiffner, 2001), with the consequence that the sediment flux to the basin decreased. This reduction in sediment flux, together with the tectonic widening of the basin, was thus likely to have shifted the basin to underfilled conditions, which could have allowed the transgression of the peripheral sea in Switzerland (Fig. 8). In addition, shifts in surface loads, caused by the rise of the Aar massif, resulted in flexural adjustments in the Molasse basin through the buckling of the foreland plate at a smaller scale. We suggest that this influenced the water depths within the basin, which could explain the development of distinct depositional settings and the formation of subtidal shoals wherein a lateral bulge is expected. Because of the formation of shallow marine conditions, subtle changes in eustatic sea level contributed to the occurrence of several hiatuses (Sant et al., 2017). Whereas these mechanisms are capable of explaining the establishment of the Burdigalian seaway and the formation of distinct sedimentological niches in Switzerland, the drainage reversal during OMM times possibly requires a change in the tectonic processes at a scale that includes the subduction history of the entire mountain range, at least between the Eastern and Central Alps. Current explanations are still speculative and await the results of ongoing research in the framework of the AlpArray initiative. At this stage, we conclude that the geodynamic processes in the Alps include subduction mechanisms, delamination of crustal material, the uplift of the Aar massif, reorganization of the drainage network, and lower sediment fluxes, which are reflected in the Swiss Molasse basin through the establishment of shallow marine conditions and a shift of the topographic axis towards more distal sites at 19 Ma. Accordingly, the Burdigalian transgression in Switzerland most likely had a tectonic driving force but with amplifications through responses occurring on the surface of the Alps and the Molasse basin.
The dataset is available from the senior author upon request.
The supplement related to this article is available online at:
FS designed the study. PG carried out the experiments and collected and interpreted the data with support by FS. The figures and photos were created by PG with support by FS. FS and PG wrote the text.
The authors declare that they have no conflict of interest.
Special thanks go to Kenneth Eriksson, Kei Ogata, and an anonymous reviewer for their constructive and insightful comments. The supportive editorial handling by Elias Samankassou is kindly acknowledged. We would also like to thank SEAG (Aktiengesellschaft für Schweizerisches Erdöl), who provided us with a scan of the seismic line BEAGBE.N780025, and Kellerhals and Haefeli AG for permission to use the photos of the Gurten drill core.
This research has been supported by the Swiss National Science Foundation (grant no. 154198).
This paper was edited by Elias Samankassou and reviewed by Kei Ogata, Kenneth Eriksson, and one anonymous referee.