Supplementary material for Large-wavelength late Miocene thrusting in the North Alpine foreland : Implications for late orogenic processes

Electrodynamic disaggregation (selFrag) To release the apatite crystals from the rock samples, we used the electrodynamic disaggregation technique (selFrag). This method exposes the rock specimen to a high voltage pulse and fractures it along its grain boundaries. As opposed to separation using a jaw crusher, this method is less time consuming and the rock is disintegrated along the grain boundaries (Giese et al., 2010). This ensures individual grains are less prone to damaging during processing.


Introduction
Deep crustal processes and slab dynamics have been considered to influence the evolution of mountain belts (e.g. Davies and von Blanckenburg, 1995;Molnar et al., 1993;Oncken et al., 2006). However, these deep-seated signals may be masked by tectonic forcing at upper-crustal levels and by enhanced surface erosion related to climate change (e.g. Champagnac et al., 2007;Chemenda et al., 2000;Ganti et al., 2016;Whipple, 2009;Willett et al., 2006). In nearsurface crustal domains, it is thus challenging to isolate the exhumation signal related to slab dynamics. In this context, foreland basins offer suitable archives as they potentially bear information that allows us to resolve the influence of deep-seated processes on mountain building. This is the case because these basins not only record signals that are related to surface dynamics such as changes in sediment fluxes and eustacy (e.g. Pippèrr and Reichenbacher, 2017;Sinclair and Allen, 1992), but they also preserve information on the tectonic processes at the crustal and possibly mantle scales that operate on long timescales and at long spatial wavelengths (e.g. DeCelles and Giles, 1996;Garefalakis and Schlunegger, 2018;Leary et al., 2016). The north Alpine foreland basin, or Molasse Basin, is particularly suited to constrain the geodynamic evolution of the collisional Alpine orogen because the history of this sedimentary trough has been well established through numerous magneto-and tectonostratigraphic (e.g. Burkhard and Sommaruga, 1998;Ganss and Schmidt-Thomé, 1953;Homewood et al., 1986;Kempf et al., 1999;Pfiffner, 1986;Schlunegger et al., 1996;Sinclair et al., 1991), seismic (Hinsch, 2013;Mock and Herwegh, 2017;Ortner et al., 2015;Sommaruga et al., 2012), and low-temperature thermochronological analyses (Cederbom et al., 2004(Cederbom et al., , 2011Gusterhuber et al., 2012;von Hagke et al., 2012von Hagke et al., , 2014bMazurek et al., 2006).
Studies from the forelands of the European Alps have shown that the most external parts of the orogen were incorporated into the orogenic wedge in Miocene times (e.g. Becker, 2000;Burkhard, 1990;von Hagke et al., 2012von Hagke et al., , 2014bHinsch, 2013;Ortner et al., 2015;Pfiffner, 1986;Schmid et al., 1996;Schönborn, 1992). In the case of the Swiss Molasse Basin, late-Miocene deformation has been kinematically and spatially linked to the uplift and exhumation of the external crystalline massifs (ECMs), which represent basement units derived from the subducting European plate. This late-Miocene deformation stage was also linked to the main deformation in the Jura fold-and-thrust belt (FTB) situated at the northern margin of the Molasse Basin (Figs. 1 and 2; e.g. Boyer and Elliot, 1982;Burkhard, 1990;Burkhard and Sommaruga, 1998;von Hagke et al., 2012;Laubscher, 1961Laubscher, , 1992Mosar, 1999;Pfiffner, 1986;Sommaruga, 1999). The inferred linkages between the uplift of the ECMs, the imbrication of the proximal molasse deposits, and the main deformation in the Jura FTB are mainly based on a classical scenario of ongoing continent-continent collision, where compressional wedge tectonics and shortening result in crustal thickening of basement units in the hinterland and a propagation of the orogenic wedge towards the foreland, including imbricate thrusting (Pfiffner, 1986;Rosenberg and Berger, 2009;Schmid et al., 1996).
New studies from the Aar Massif of the Central Alps challenge this view. Based on geometric, kinematic, metamorphic, and geodynamic arguments, Herwegh et al. (2017Herwegh et al. ( , 2020 suggest that in the eastern ECMs (Aar, Mont Blanc, and Aiguilles Rouges massifs) the mode of uplift and consequentially exhumation switches between "vertical" and "horizontal tectonics". Note that the terms "vertical" and "horizontal tectonics" as used by Herwegh et al. (2017Herwegh et al. ( , 2020 are based on geometric and kinematic considerations; i.e. they imply a steeper or less steep orientation of the main faults along which strain is accommodated. Thus, in a compressional framework, these terms describe whether the vertical or the horizontal displacement components play the dominant role at a given point in time. In other words, horizontal shortening also continues during phases of vertical tectonics. The evolution of the Helvetic fold-and-thrust belt with the associated nappe stacking represents an Oligocene stage of horizontal tectonics. A major switch to vertical tectonics occurred in early Miocene times when major parts of the Aar and Mont Blanc-Aiguilles Rouges massifs experienced a strong vertical displacement component mainly along steep to subvertical reverse faults (see Herwegh et al., 2017, 2020, andreferences therein). This resulted in the differential rock uplift and the doming of these eastern ECMs as well as in a passive upward bulging of the entire nappe stack above. From a geodynamic point of view, (i) a retreating European slab, (ii) delamination of lower crust from the mantle, and (iii) a buoyancy-driven rock uplift component within a compressional regime were suggested by the aforementioned authors to be the driver for this vertical-tectonic forcing. It is important to stress that this vertical displacement component became dominant at a time when the un-thinned or only slightly thinned part of the buoyant former European passive continental margin entered the subduction channel. This stage of vertical tectonics in the eastern ECMs was followed by a second, late-Miocene phase of horizontal tectonics and associated en bloc exhumation (Herwegh et al., 2020), which in the Aar Massif was accomplished through slip along shallow south-dipping basal thrust systems Herwegh et al., 2017Herwegh et al., , 2020Mair et al., 2018). A similar scenario of a late-stage Miocene en bloc exhumation due to NWvergent thrusting was proposed for the evolution of the Mont Blanc-Aiguilles Rouges massifs (see Herwegh et al., 2020).
In addition to this ECM-related information, structural, chronostratigraphic, and low-temperature thermochronological data (Burkhard and Sommaruga, 1998;von Hagke et al., 2012von Hagke et al., , 2014bOrtner et al., 2015) suggest that the Alpine deformation front remained stationary within the Subalpine Molasse during the time of vertical tectonics. The deformation then propagated to the Jura FTB during the late Miocene when the vertical tectonics was superseded by the thrust-related exhumation in the Aar Massif (e.g. Burkhard and Sommaruga, 1998;von Hagke et al., 2012;Herwegh et al., 2020). The corresponding along-strike partitioning of strain between the Jura FTB and the Subalpine Molasse has, however, only been incompletely resolved so far (e.g. Burkhard, 1990). Likewise, it has been unclear whether the amount of shortening within the Subalpine Molasse is consistent along-strike, particularly with respect to the highly non-cylindrical architecture of the Alpine hinterland as well as the lentoid-shaped map appearance of the eastern ECMs ( Fig. 2; Burkhard, 1990). However, this information is vital for understanding how deformation in the ECMs is potentially linked to foreland FTB tectonics. Accordingly, as a first and major contribution of our paper, we aim at reconstructing the chronology and amount of shortening of the molasse sequences during the late stage of Alpine orogeny since the mid-Miocene. We present new low-temperature thermochronological data from the Lake Thun region, which we consider a key area for the understanding of Subalpine Molasse tectonics for the following reasons.
i. With the Aar Massif being one of the best-studied ECM in the Alps in terms of its Alpine tectonic evolution and geodynamic significance, information on the history of thrusting of the Subalpine Molasse in front of the Aar Massif is of great importance if the scope is reconstructing the post-collisional tectonic evolution of the Central Alps. ii. The thickness of the basal décollement zone within the Mesozoic cover sediments below the Molasse Basin decreases significantly west of ca. 7 • E (Landesgeologie, 2017; Mock and Herwegh, 2017;Sommaruga et al., 2012;Sommaruga et al., 2017). Hence, the Lake Thun area is located in an ideal place where we can expect the eastward transition from the detached Plateau Molasse to the non-detached Foreland Molasse.
iii. While east of the Lake Thun area, the frontal part of the Subalpine Molasse is often delineated by a backthrust and a triangle zone, this feature is largely missing farther west, where predominantly northwestward imbricate thrusting occurs.
We compare the thermochronological data with previously published data from farther east (von Hagke et al., 2012(von Hagke et al., , 2014b. We combine this information with information from well-established work on the chronology, tectonics, and stratigraphy of the thrust Subalpine Molasse between 6.8 • E (Lake Geneva) and 12.8 • E (near Salzburg; e.g. Burkhard and Sommaruga, 1998;von Hagke et al., 2012von Hagke et al., , 2014bHinsch, 2013;Kempf et al., 1999;Ortner et al., 2015;Pfiffner, 1986;Schlunegger et al., 1997;and many others) in order to better constrain the timing and spatial pattern of thrusting in the Subalpine Molasse. We will relate this to the history of shortening of the Jura FTB and to Alpine tectonic events in an effort to identify possible relationships with the geodynamic forcing at a larger scale.
Geological mapping of the Subalpine Molasse (e.g. Ganss and Schmidt-Thomé, 1953;Haldemann et al., 1980;Weidmann et al., 1993;Zaugg et al., 2011) as well as stratigraphic work (e.g. Bachmann and Müller, 1992;Kempf et al., 1999;Lemcke, 1988;Schlunegger et al., 1993Schlunegger et al., , 1997Schlunegger and Kissling, 2015) has shown that the proximal basin border is characterized by large orogen-parallel lithologic variations where kilometer-thick conglomerate suites with high mechanical strengths alternate with mudstones and sandstones with low at-yield conditions over lateral distances of a few kilometers. As a consequence, the patterns of thrust faults and folds within the Subalpine Molasse change along-strike. In regions where kilometer-thick conglomerate packages are present, the geometry of the Subalpine Molasse is characterized by kilometer-spaced thrust faults with relatively large  Buxtorf, 1916;Herwegh et al., 2017;Mock and Herwegh, 2017;Pfiffner, 2009;Rosenberg and Kissling, 2013). (b) Engadine Window (adapted from Rosenberg et al., 2015). (c) Western Tauern Window (adapted from Rosenberg et al., 2015). displacements on them (e.g. Ortner et al., 2015;von Hagke et al., 2012von Hagke et al., , 2014b. In areas however, where the foreland mechanical stratigraphy mainly consists of sandstone-mudstone alternations, the structural style is characterized by closely spaced folds and thrust faults with possibly smaller displacements (e.g. Kempf et al., 1999;Ortner et al., 2015;Schlunegger et al., 1993Schlunegger et al., , 1997. We expect that these differences in mechanical stratigraphy leave an imprint on where strain was accommodated during late orogenic shortening. Therefore, as a related second contribution, we aim at exploring how lithological variations within the foreland basin, which were controlled by the paleogeographic conditions, contributed to the spatial distribution of late orogenic strain.
2 Geological and tectonic setting

The Central Alps
Classically, the boundary (or transition) between the Central and the Eastern Alps has been located in eastern Switzerland (e.g. Schmid et al., 2004). There, remnants of the Piemont-Ligurian Ocean separate units of European origin (Central and Western Alps) from thrust nappes that are part of the Adriatic continental plate (Figs. 1 and 2; Eastern Alps; Schmid et al., 2004). However, this surface observation is not reflected in the slab geometry at depth. Instead, seismotomography studies have shown that the subducting European slab beneath the Central Alps extends to the east until ca. 11-13 • E, i.e. between the Giudicarie-Brenner fault system and the central Tauern Window, where a segmenta-tion of the slab structure is observed (Hetényi et al., 2018b;Kästle et al., 2020;Lippitsch et al., 2003;Mitterbauer et al., 2011;Zhao et al., 2016; see further details in the "Discussion" section). Similar to Rosenberg et al. (2018), who place the boundary between the Central and Eastern Alps at the Giudicarie-Brenner fault system, we shift this boundary farther east than classically placed (e.g. Schmid et al., 2004). We thus follow the argumentation given above, which is also advocated by Kissling and Schlunegger (2018), and place the Central Eastern Alps transition where the current slab configuration changes at depth. This is particularly important when the influence of subduction mechanisms on the surface geology during the past few millions of years is considered, which is the scope of the discussion of the results at the end of this paper. In this study, we therefore focus on the Central Alps because late-Miocene thrusting in the Subalpine Molasse terminates near Salzburg ( Fig. 1; Beidinger and Decker, 2014;Hinsch, 2013;Ortner et al., 2015), thus spatially coinciding with the aforementioned along-strike changes in the configuration of the lithospheric mantle slabs at depth between ca. 11 and 13 • E. However, we will elaborate on the influences of Eastern Alpine tectonic processes in the discussion part of this paper.
The Central Alps are situated almost entirely on top of the subducted Central European lithospheric slab (Fig. 2a and b;Schmid et al., 2017), which steeply dips into the asthenospheric mantle as imaged by teleseismic tomography (Lippitsch et al., 2003). The orogen is the result of the Late Cretaceous to Eocene subduction of the oceanic part of the European plate under the Adriatic continental plate and subsequent post-35 Ma continent-continent collision (e.g. Schmid et al., 1996). Subduction initially affected the European oceanic lithosphere, which was followed by the subduction of the distal spur of the Iberian plate that became the Penninic thrust nappes. During the Eocene, subduction started to involve the distal and stretched European continental margin (e.g. Cardello et al., 2019;Mosar et al., 1996;Stampfli and Marchant, 1997). The subduction system then became clogged when the thick and buoyant European crust started to enter the subduction channel. As a result, the subducted oceanic slab of the European plate supposedly broke off, which has been inferred from widespread plutonism ca. 32 Myr ago (Davies and von Blanckenburg, 1995). Slab unloading and basal accretion of crustal segments to the upper plate resulted in rock and surface uplift. This gave way to the incipient rise of the Alpine topography and consequently exhumation of uplifted rocks. Uplift was mainly accomplished through backthrusting along the Insubric Line (Hurford, 1986;Schmid et al., 1996), and this was also manifested in the increase in sediment discharge into the Alpine foreland (Schlunegger and Castelltort, 2016;Sinclair, 1997).
During the Oligocene, Europe-derived sedimentary units were sheared off from their substratum and emplaced to the north over a distance of several tens of kilometers, thereby forming the present-day Helvetic cover nappes (Pfiffner, 2011, and references therein). Ongoing convergence resulted in delamination of lower European crustal segments from the lithospheric mantle (Fry et al., 2010). This mechanism eventually resulted in the inferred buoyancy-driven subvertical uplift of the thickened crust (Kissling and Schlunegger, 2018) along steeply dipping shear zones and in the exhumation of the ECMs ca. 20 Myr ago ( Fig. 2a; e.g. Egli et al., 2017;Glotzbach et al., 2011;Herwegh et al., 2017Herwegh et al., , 2020. Farther east, between the Aar Massif and the Brenner Fault (a segment, which we here regard as the eastern Central Alps), such Europe-derived crustal blocks are not exposed at the surface. East of the Brenner Fault ( Fig. 1), northward indentation of the Adriatic plate resulted in eastwarddirected lateral extrusion of crustal blocks, which started in early Miocene times (Frisch et al., 1998;Handy et al., 2015;Ratschbacher et al., 1991). This process was associated with the exhumation of the Tauern Window, which was accomplished by upward folding and erosion of the nappe pile and by normal faulting along its bounding low-angle normal faults (see Rosenberg et al., 2018, and references therein).
Late stages in the evolution of the Central Alps are dominated by orogen-perpendicular growth when the deformation front of the Alps propagated to the external parts, resulting in the development of the Jura FTB and the Subalpine Molasse in the north and in the ongoing growth of the Southern Alps on the southern side of the Alps (Figs. 1 and 2; Burkhard, 1990;Caputo et al., 2010;Castellarin and Cantelli, 2000;Nussbaum, 2000;Pfiffner, 1986;Schmid et al., 1996;Schönborn, 1992).

Structures and tectonic evolution
The Molasse Basin of the north Alpine foreland, which contains deposits derived from the progressive erosion of the evolving Alps since 32 Ma , can tectonically be subdivided into three parts. (i) The Plateau Molasse is the gently folded part of the basin, which evolved into a wedge-top basin (e.g. Willett and Schlunegger, 2010) during the late Miocene when the Jura FTB experienced its major phase of deformation. Folding of the Jura FTB was accomplished as the Mesozoic suite of limestones and marls became detached from the underlying basement along a main décollement zone within the Triassic evaporites. To the east, with the Triassic evaporites diminishing and the Jura FTB tapering off, the basin gradually changes into a non-detached configuration, (ii) the Foreland Molasse (e.g. Berge and Veal, 2005;Ortner et al., 2015;Pfiffner, 1986). At the southern, proximal basin border, (iii) the Subalpine Molasse extends continuously as a narrow band of imbricates from south of Geneva to Salzburg, where it disappears beneath Helvetic and Penninic units before it emerges again in upper Austria. In this contribution we mainly focus on the Subalpine Molasse between Lake Geneva and Salzburg, i.e. the part of the fold-thrust belt where the tectonic processes may be associated with the subduction of the European lithospheric slab of the Central Alps.
The Subalpine Molasse consists of south-dipping imbricated thrust sheets. In large parts and predominantly northeast of the Lake Thun area, the structures of the Subalpine Molasse also include the north-dipping backthrusts that form a triangle zone at the transition to the Plateau and Foreland Molasse ( Fig. 2; Berge and Veal, 2005;Fuchs, 1976;Müller et al., 1988;Ortner et al., 2015;Schuller et al., 2015;Sommaruga et al., 2012). The Subalpine Molasse started to become incorporated into the orogenic wedge shortly after deposition in Oligocene times (Hinsch, 2013;Kempf et al., 1999;Pfiffner, 1986). After ca. 20 Ma, contemporaneously with the development of the frontal triangle zone, the northern Alpine thrust front remained stationary in the area of the Subalpine Molasse (Burkhard and Sommaruga, 1998;von Hagke et al., 2014b;Ortner et al., 2015). It was not until ca. 12 Ma when parts of the deformation in the western Molasse Basin (i.e. west of Lake Constance) propagated along a basal décollement zone within the Triassic evaporites into the thin-skinned Jura FTB (Becker, 2000;Burkhard and Sommaruga, 1998;Laubscher, 1961;Philippe et al., 1996), although some deformation occurred already in the late Oligocene in the area of today's Jura FTB (Aubert, 1958;Liniger, 1967). Hence, during the late Miocene, the Molasse Basin experienced along-strike changes in tectonic style and locus of deformation. While sediment accumulation in the eastern Molasse Basin (i.e. here the Foreland Molasse between Lake Constance and Salzburg) still occurred in a foredeep setting, at this stage the western part evolved into a wedge-top basin (i.e. Plateau Molasse; Willett and Schlunegger, 2010). Both the evaporite basal décollement and the thrusts of the Subalpine Molasse are considered to root below and in the ECMs. Accordingly, they were kinematically linked to the late-Miocene exhumation of the ECMs ( Fig. 2a; e.g. Burkhard, 1990), which was driven at that time by north-directed thrusting along shallow-dipping faults (Herwegh et al., , 2020, thereby causing a phase of accelerated exhumation at ca. 10 Ma (Fox et al., 2016;Glotzbach et al., 2010;Valla et al., 2012;Vernon et al., 2009;Weisenberger et al., 2012). Based on stratigraphic data, Haus (1935) inferred already in 1935 that the Subalpine Molasse of Central Switzerland was subject to major thrusting in late-Miocene times. Recently, studies have revisited this topic and documented that the Subalpine Molasse was subject to break-back thrusting between ca. 13 and 4 Ma (von Hagke et al., 2012(von Hagke et al., , 2014bOrtner et al., 2015;Schuller et al., 2015), which was thus coeval to folding and thrusting in the Jura FTB. We note that late-Miocene thrusting is not recorded in the Subalpine Molasse east of Salzburg (Beidinger and Decker, 2014;Hinsch, 2013;Ortner et al., 2015).

Stratigraphic development
The clastic infill of the Oligocene to Miocene peripheral Molasse Basin consists of eroded sediments of the evolving Alps and in the northern parts partly of material shed from the Black Forest and the Bohemian Massif. Accommodation space was formed through subsidence, classically related to flexural bending of the European plate in response to the combined effect of subsurface slab loading and topographic loading of the advancing Alpine thrust wedge during Paleogene and Neogene times Burkhard and Sommaruga, 1998;Karner and Watts, 1983;Pfiffner, 1986;Zweigel et al., 1998). For the Central Alps, this view has recently been challenged by Schlunegger and Kissling (2015), who favor a slab rollback mechanism to explain foreland plate flexure and accommodation space formation. In this scenario, slab rollback is driven by the interplay between vertically directed slab loads exerted by the subducted European lithospheric mantle and buoyancy-driven crustal delamination (Kissling and Schlunegger, 2018).
In the Swiss part of the basin, the molasse sediments form two regressive and coarsening-upward megasequences Kuhlemann and Kempf, 2002;Schlunegger et al., 2007). The first megasequence describes the transition from Rupelian (ca. pre-30 Ma) sedimentation in underfilled conditions to Chattian-Aquitanian (ca. 28-20 Ma) sedimentation when the basin was overfilled. The Burdigalian (ca. 20-17 Ma) and post-Burdigalian (ca. post-17 Ma) stratigraphic records then chronicle the second megasequence during filled to overfilled conditions (Sinclair and Allen, 1992). West of ca. 11 • E, large alluvial megafans developed at the mountain front during the overfilled stage of the basin (Frisch et al., 1998;Kuhlemann and Kempf, 2002;Ortner et al., 2015). At their cores, close to the apex, large and kilometer-thick conglomerate sequences were deposited, while at the margins, the sedimentation was mainly sandand mudstone dominated. At the proximal basin border, numerous locally derived torrents resulted in the construction of bajada fans, which further contributed to the high alongstrike stratigraphic variability at the proximal basin border (Kempf et al., 1999;Schlunegger et al., 1997;Spiegel et al., 2001). Farther east, the Molasse Basin prevailed in an underfilled stage until at least 17 Ma, when sedimentation of sandstones and marls occurred under brackish to shallow marine conditions (Hinsch, 2013;Kuhlemann and Kempf, 2002;Lemcke, 1988;Ortner et al., 2015). Large alluvial fans are missing due to the channelizing effect of the paleo-Inn river, which transported the erosional detritus effectively farther to the east (Frisch et al., 1998;Kuhlemann and Kempf, 2002). For more details on the laterally varying stratigraphy, including Wheeler diagrams, the reader is referred to Kempf et al. (1999), Schlunegger et al. (1997), and Schlunegger and Norton (2013) as well as Hinsch (2013), Kuhlemann and Kempf (2002), Lemcke (1988), and Ortner et al. (2015) for the Molasse Basin west and east of Lake Constance, respectively.

Litho-tectonic architecture of the Lake Thun area
In the sampling area, mapping shows that the litho-tectonic architecture of the Subalpine Molasse contrasts between both sides of Lake Thun ( Fig. 3; Beck, 1945;Blau, 1966;Haus, 1937;Jordi, 2012;Landesgeologie, 2005;Mock and Herwegh, 2017;Rutsch, 1947;Schlunegger et al., 1993Schlunegger et al., , 1997Vollmayr, 1992). On the eastern side, the basin is generally made up of amalgamated conglomerates, which are cut from west to east by two narrow sand-and mudstone-dominated bands. The south-dipping thrust north of the Falkenfluh anticline as well as the three tightly spaced thrusts further south all emerge along these mechanically weaker bands (Fig. 3a). In between these thrust domains, we observe large tectonic slices (6-8 km wide) made up of mechanically strong conglomerates of Burdigalian (north) and Chattian-Aquitanian age (south). The latter is bordered to the north by three thrusts for which we can infer a break-back (i.e. back-stepping) sequence of thrusting from cross-cutting relationships in map view (Fig. 3a). To the south, across a south-dipping thrust, Rupelian sand-and mudstones crop out. West of the Aar valley, conglomerates are largely absent and the stratigraphy is sand-and mudstone dominated. Compared to the Subalpine Molasse east of the Aar valley, the FTB is much narrower, thrusts are more evenly distributed, and the tectonic slices are narrower (2-3 km; Fig. 3). Across the whole area we can see a clear correlation between the distribution of mechanically weak sand-and mudstones and the location and alignment of thrust faults. The subsurface continuation of the thrusts has been inferred from seismic interpretation (for more details, see Mock and Herwegh, 2017), which suggests that the thrusts merge in a detachment level at a depth of ca. 2 km (Fig. 3b). However, the type of lithology of the detachment remains unknown, since the poor signal-to-noise ratio of the seismic images does not allow a clear distinction in this regard. Due to the lack of subsurface data further south it also remains speculative how the structures extend southward. As indicated in the sections (Fig. 3b), it is, however, likely that the thrusts root within the Rupelian sand-and mudstones at the base of the molasse sequence on top of the Mesozoic cover sediments. In summary, there are pronounced differences both in the mechanical stratigraphy as well as in the structural configuration of the Subalpine Molasse west and east of the Aar valley. We discuss this in detail in Sect. 5.1 below.

Selection of sample sites and litho-tectonic architecture
We collected 12 samples across the Subalpine Molasse east and west of Lake Thun for apatite (U−Th−Sm)/He (AHe) dating (Fig. 3). The northernmost samples represent coarsegrained Burdigalian sandstones of the Plateau Molasse. Further to the southwest, we sampled Chattian-Aquitanian and Rupelian sandstones within the Subalpine Molasse. Samples have been collected in the hangingwalls and footwalls of the individual thrusts (Fig. 3). To control the depositional age of the sediments, samples were taken in the vicinity of sites with known mammal ages and magneto-polaritybased chronologies (Schlunegger et al., 1996;Strunck and Matter, 2002) wherever possible. In order to test the influence of along-strike changes in the mechanical stratigraphy of the Subalpine Molasse on the thrusting pattern, we compiled more details about the geological architecture of the proximal molasse from published geological maps (Landesgeologie, 2005) and use available tectonic sections (Ortner et al., 2015;Sommaruga et al., 2012) to illustrate the related tectonic style.

(U−Th−Sm)/He (AHe) thermochronology and thermal modeling
We determined the most recent exhumation history of the Subalpine Molasse of the Lake Thun region through AHe dating. This method is based on the α decay of 238 U, 235 U, 232 Th, and 147 Sm isotopes and the retention of the radiogenic product 4 He in the crystal lattice below a certain temperature (e.g. Farley, 2002). Diffusive loss of 4 He in the lattice depends on the grain size, shape, chemical composition, distribution of the mother isotopes, radiation damage density, and the time-temperature evolution of the crystal (Farley, 2002;Wolf et al., 1996). Consequently, AHe ages can provide estimates for about the time when the mineral passed through the diffusion-sensitive temperature interval between ca. 80 and 40 • C , which is referred to as the partial retention zone (PRZ). Hence, this method allows constraining the tectono-thermal history of the studied rocks in the uppermost few kilometers of Earth's crust. Detrital apatite grains deposited in sedimentary basins primarily carry a cooling history of the hinterland at the time of erosion. Subsequent burial due to sedimentation or tectonic loading may reheat the detrital grains to temperatures above the closure temperature, thereby resetting the chronometer. Subsequent exhumation will chronicle the basin's exhumation, whereas grains that have not been reset during the basin's burial history still carry a signal of older cooling events. Consequently,  (2017), Rutsch (1947), and Schlunegger et al. (1993Schlunegger et al. ( , 1997. Lithological information is compiled from Landesgeologie (2005) and Schlunegger et al. (1993).
the relation between cooling age and stratigraphic age may provide estimates on the burial as well as the exhumation history (e.g. Reiners and Brandon, 2006). We used a combination of standard techniques for the separation of apatite minerals, which particularly includes electrodynamic disaggregation, and magnetic and heavy liquid separation. Single crystals were handpicked under a binocular and checked for inclusions and imperfections under an optical microscope with cross-polarized light (more information on the mineral separation techniques and picking criteria are given in Sect. S1 in the Supplement). Helium extraction and measurement of parent isotope contents has been conducted at the GÖochron laboratories of the University of Göttingen. Raw ages were corrected for α ejection (Table 1). We measured four to eight single-grain ages per sample and calculated average ages using the unweighted arithmetic mean for completely reset samples. We excluded single-grain ages for subsequent geological interpretation based on the following criteria (Table 1): (i) high analytical errors (> 10 %), (ii) very low U content (< 10 ppm), (iii) a substantial amount of He on the first re-extract (> 4 %), or (iv) erroneous old ages stemming most likely from U−Th-rich mineral inclusions which produce parentless He. For the latter case, we plotted the He content versus the present-day He production rate in order to detect these ages (for details, see Vermeesch, 2008).
We constrained the thermal histories of the sampled sediments (Fig. 4) by modeling the AHe age data with the HeFTy software (Ketcham, 2005). We gave the algorithm a high degree of freedom by using only a few modeling constraints: the age of sedimentation, the paleo-temperature (Mosbrugger et al., 2005), the present-day annual average temperature, and maximum post-depositional heating rates inferred from maximum sedimentation rates (Schlunegger and Norton, 2015) and a paleo-geothermal gradient of 28 • C km −1 (Schegg and Leu, 1998).

AHe age data
Samples SM-7 and SM-16 from the Plateau Molasse show a wide spread in AHe ages, indicating a partially exhumed fossil PRZ. Hence, they did not experience enough postdepositional heating for a full reset of the thermochronometer ( Fig. 3b and Table 1). All other samples show single-grain ages, which plot within error for each sample. Ages are significantly younger than the corresponding depositional ages and are thus considered to represent completely reset ages (fully exhumed fossil PRZ), hence inferring substantial postdepositional burial and heating.
Average ages for completely reset samples range from 6.0 ± 0.4 to 10.8 ± 0.6 Ma ( Fig. 3b and Table 1). All samples from the Subalpine Molasse show post-depositional burial and heating to > 60 • C (Fig. 4). Sample SM-7 from the Plateau Molasse did not experience enough post-depositional heating to fully reset the thermochronometer (Fig. 4). Thermal modeling supports the young exhumation of the Subalpine Molasse between 12 and 5 Ma (Fig. 4). The thermal histories of SM-7 and SM-13 indicate an exhumation signal of the Plateau Molasse at ca. 10 Ma.
Samples SM-8, SM-11, and SM-15 were collected from the same tectonic sliver (Figs. 3a and 5), and corresponding average ages of 6.0 ± 0.4, 6.6 ± 0.4, and 6.1 ± 0.4 Ma group well within their margins of error. A jump in average ages occurs across the thrust to the adjacent tectonic slivers and the Plateau Molasse in the north, where samples SM-12, SM-13, and SM-14 yield average ages of 9.2 ± 1, 8.9±0.7, and 10.8±0.6 Ma, respectively, thus also grouping within their margins of error. This pattern is also reproduced by the modeled thermal histories, where a jump in exhumation ages across tectonic boundaries is recognized (Fig. 4). In the western area, the average ages do not show such a close correlation with the tectonic position (Figs. 3a and 5). However, the thermal models for samples SM-5 and SM-6 indicate younger cooling than for samples SM-4 and SM-7 (Fig. 4), thus suggesting a tectonic control on exhumation.

Late-Miocene shortening estimates
Based on published restored cross sections (Burkhard and Sommaruga, 1998;von Hagke et al., 2014b;Ortner et al., 2015), retro-deformed and balanced palinspastic maps (Philippe et al., 1996), thermochronological data (von Hagke et al., 2012(von Hagke et al., , 2014b, and our own cross section restorations (see Sect. S2 and Fig. S1 in the Supplement), we estimated the amount of late-Miocene horizontal shortening of the Subalpine Molasse and the Jura FTB from Lake Geneva to Salzburg (Fig. 6b). Post-12 Ma horizontal shortening in the Jura FTB decreases from a maximum of ca. 32 km in the west to 0 km at the eastern tip (Philippe et al., 1996). Contrariwise, minimum horizontal shortening in the Subalpine Molasse increases from ca. 10 km in the west (Lake Geneva) to ca. 20 km farther east (Lake Constance) before decreasing again to below 1 km in the area near Salzburg (Beidinger and Decker, 2014;Burkhard and Sommaruga, 1998;von Hagke et al., 2012von Hagke et al., , 2014bHinsch, 2013;Ortner et al., 2015). These values do not account for the shortening taken up by the frontal triangle zone between ca. 20 Ma and 12 Ma (von Hagke et al., 2014b;Kempf et al., 1999;Ortner et al., 2015). We acknowledge that shortening estimates for the Subalpine Molasse are subject to uncertainties (Burkhard and Sommaruga, 1998;von Hagke and Malz, 2018;Ortner et al., 2015) and likely represent minimum estimates due to (i) the unconstrained large parts of proximal Subalpine Molasse, which are hidden below the frontal thrusts of the Helvetic nappes and Penninic Klippen units, and (ii) the nonpreservation of the hangingwall cut-offs of individual thrust sheets. Despite these uncertainties, we observe that the late- The amount of helium is given in nano-cubic-centimeters in standard temperature and pressure. b Ejection correction (Ft): correction factor for α ejection Hourigan et al., 2005). c Uncertainty of the single-grain age is given as 2σ in percent (or in millions of years), and it includes both the analytical uncertainty and the estimated uncertainty of the Ft. d Average ages for totally reset samples were calculated as the unweighted arithmetic mean. g Uncertainty of the sample average age is 2 standard error as (SD)/(n)1/2, where SD is standard deviation of the age replicates and n is the number of age determinations. f Ages with a substantial first He re-extract (> 4 %) and/or a total analytical error of > 10 % have been excluded. Outliers on the [He]−P plot have also been excluded (Vermeesch, 2008).
Miocene cumulative amount of shortening in the Subalpine Molasse and the Jura FTB decreases from the west to the east from ca. 30 km near Geneva to ca. 20 km near Lake Constance before decreasing to < 1 km between ca. 10.5 • E and Salzburg (12.8 • E; Fig. 6b).

Uncertainty of the thermochronological data
Interpreting thermochronological ages from sedimentary basins may be challenging, as they potentially have been influenced by factors that are difficult to quantify. For instance, heat flux within the basin may have been variable through time when hydrothermal fluids may have heated or cooled the grains after deposition (Louis et al., 2019), or individual grains dated have different provenance histories that change their respective temperature sensitivity (Luijendijk, 2019). These uncertainties can result in a spread of single-grain ages within individual samples or along-strike the orogen, which cannot be explained by factors such as effective uranium content (e.g. Flowers et al., 2009) or broken grains (Brown et al., 2013). Provenance age is a particular problem for the interpretation of samples which have not been reset because, commonly, only three to five single-grain ages are reported, which makes statistical analysis of the grain age distributions difficult.
In our data set, we report two samples that show unreset ages (samples SM-7 and SM-16). These samples can only be used to infer a maximum amount of cooling of less than 60 • C. Interpreting the single-grain ages either as detrital ages or partially reset ages and deriving t-T histories for exhumation of the hinterland would require additional data. For completely reset samples, the provenance history of the single-grain ages is erased (Luijendijk, 2019). Single-grain ages of all other samples in our data set  may be influenced by heat flux variations. Paleo-heat flux is notoriously difficult to constrain (Louis et al., 2019). However, as shown from a compilation of geothermal springs within the Alps, heat flux variations are spatially localized (Luijendijk et al., 2020). Consequently, this should be reflected in local outliers within the data set, for instance a suite of very young ages within a region. As our data are consistent with other data along-strike the Subalpine Molasse as far east as Lake Constance (von Hagke et al., 2012(von Hagke et al., , 2014b, we consider local paleo-heat flux variations to be negligible and infer that the data reflect cooling associated with thrusting in the Subalpine Molasse.  (Ketcham, 2005) show the time-temperature history of the samples discussed in the text. Modeling constraints are shown as black boxes. The bold black lines at 10 Ma serve as a visual time reference. The thermal histories for the different samples are aligned from north to south and from east to west according to their sample location.

Local-scale stratigraphic architecture conditioning the pattern of strain accommodation
In the sampling area, mapping discloses along-strike differences in the litho-tectonic architecture of the Subalpine Molasse (Fig. 3). The sediments are thrust northwestward along SW-NE-striking thrusts. East of the Aar valley, a backthrust emerges at the surface and forms a frontal triangle zone (Figs. 3 and 5), a structure which is also known from parts of the Subalpine Molasse farther east ( Fig. 7; e.g. Berge and Veal, 2005;von Hagke and Malz, 2018;Müller et al., 1988;Ortner et al., 2015;Schuller et al., 2015 (Beck, 1945;Haus, 1937;Jordi, 2012;Rutsch, 1947;Schlunegger et al., 1993Schlunegger et al., , 1997 and subsurface (2D seismic interpretation; Mock and Herwegh, 2017) geological information. PRZ, partial retention zone. al., 2012; Stäuble and Pfiffner, 1991). The Aar valley, running across the study area, is characterized by a low relief and is filled by > 100 m thick Quaternary deposits (Fig. 5;Dürst Stucki et al., 2010). Accordingly, the structural configuration of this part of the study area is only poorly resolved. However, the structures of the Subalpine Molasse change abruptly across the Aar valley (Figs. 3 and 5), as has been described by many authors (Beck, 1945;Blau, 1966;Haus, 1937;Pfiffner, 2011;Rutsch, 1947;Vollmayr, 1992). Based on this observation and under the consideration of available interpreted reflection seismic lines, the presence of a possible syn-thrusting strike-slip fault zone running along the valley axis has been proposed (Mock and Herwegh, 2017;Pfiffner, 2011;Vollmayr, 1992). The presence of such a fault is, however, speculative due to the low resolution of the seismic data. The latter stems from (i) the thick Quaternary cover resulting in a very poor signal-to-noise ratio and (ii) the poorly resolved molasse strata as a result of the frequency, which was chosen in order to optimize for the targeted Mesozoic horizons below (Mock and Herwegh, 2017). An alternative explanation for the sudden along-strike change in the tectonic architecture has first been proposed by Rutsch (1947). He reported that the change from a mainly conglomeratic (east) to a sand-and mudstone-dominated (west) lithofacies coincides with an increase in the folding intensity along the frontmost anticline (Falkenfluh anticline; Fig. 3a). Indeed, while conglomerates are the dominant lithofacies in the eastern part of the study area, they are vastly absent west of the Aar valley, where mainly alternating sequences of sandstones and mudstones outcrop (Figs. 3 and 7;Landesgeologie, 2005). This difference in mechanical stratigraphy east and west of the Aar valley is probably a result of an asymmetric dispersal system with a distinct northeastward direction of sediment discharge (Schlunegger and Norton, 2015). The distribution of mechanically different lithologies seems to control the pattern of thrusting. While east of the Aar valley, the mechanically stronger thick conglomeratic sequences deform en bloc (i.e. as one unit) and thrusts are concentrated in narrow bands following mechanically weak zones of sand-and mudstones, strain is accommodated in a much more distributed pattern along more closely spaced thrusts in the western part of the study area (Fig. 3). With deformation localizing along the pronounced mechanical contrast between conglomerates and sand-or mudstones, the large tectonic slice made up of amalgamated conglomerates east of the Aar valley is thrust northwestward en bloc. This is also manifested in well-defined and constant AHe ages of ca. 6 Ma within this tectonic slice (samples SM-8, Figs. 3 and 5).
The observation that along-strike variations in the stratigraphic architecture lead to complex patterns of strain accommodation is not unique to our study area but can be also made at the Mont Pèlerin, the Rigi, or the Hörnli in western, central, and eastern Switzerland, respectively (Fig. 7). Salients seem to occur at the apex of former alluvial megafan depositional systems, while recesses are observed in regions in between these large dispersal systems, where thick conglomerate sequences are missing. The lithological control on the strain accommodation and exhumation pattern is particularly well observed in the Rigi area (profile 3 in Fig. 7; Sommaruga et al., 2012). The thick conglomerate sequence of the Rigi thrust sheet was (re-)activated en bloc at ca. 5 Ma, while the adjacent sand-and mudstone-dominated part to the north experienced a period of thrusting at ca. 9 Ma along evenly spaced faults (von Hagke et al., 2012). Similar dependencies between the style of deformation and lateral changes in lithology have also been described for the Subalpine Molasse in Bavaria and western Austria (for detailed information, see Ortner et al., 2015), as well as for the basal detachments of the Subalpine Molasse and the Jura FTB (von Hagke et al., 2014a). While large alluvial fans formed thick conglomerate-sandstone sequences in western Bavaria, brackish to shallow marine conditions prevailed farther east where marls and sandstones were deposited (Frisch et al., 1998;Kuhlemann and Kempf, 2002). This lateral change in the molasse's mechanical stratigraphy has been described as having a direct influence on the deformation style of the Subalpine Molasse (Ortner et al., 2015). While stacks of tectonic horse structures and a pronounced triangle zone developed in western Bavaria, the deformation style changes to buckle folding farther east, which decreases in amplitude and the triangle zone disappears (a detailed description is provided in Ortner et al., 2015). Numerical models of syntectonic sedimentation support the idea that sediments shed on an evolving FTB strongly control its geometry and may include the formation of backthrusts (e.g. Fillon et al., 2013).

The link to exhumation of the external crystalline massifs
Classically, thrusting in the Jura FTB and the Subalpine Molasse between Lake Geneva and Lake Constance has been kinematically and spatiotemporally related to the uplift and exhumation of the ECMs and to the propagation of the deformation front towards the foreland (Burkhard, 1990;Burkhard and Sommaruga, 1998;Pfiffner et al., 1997). Note that according to Campani et al. (2012), von Hagke et al. (2014b, and Schlunegger and Kissling (2015), the mean elevation of the Alpine topography has not changed since the early Miocene, and we can assume steady-state conditions between rock uplift and surface erosion for this time interval. In this case, we can consider rock uplift and exhumation rates to be similar. The Aar, and possibly also the Mont Blanc, and Aiguilles Rouges massifs were exhumed in early to mid-Miocene times during a stage of differential rock uplift along subvertical reverse faults (vertical tectonics; Herwegh et al., 2017Herwegh et al., , 2020. This stage of vertical tectonics was followed by late-Miocene en bloc rock uplift and exhumation above a series of shallow southwest-dipping basal thrusts at the massif's northwestern front (Egli et al., 2017;Herwegh et al., 2017Herwegh et al., , 2020. It is during this late (ca. post-13 Ma) thrusting-dominated stage (i.e. horizontal tectonics) that deformation was translated northward into the Alpine foreland. In the non-detached part of the north Alpine foreland, thrusting of the proximal molasse occurred. In the western north Alpine foreland, large-scale strain partitioning occurred when deformation propagated 50-90 km to the north forming the Jura FTB, while at the same time imbricates of Subalpine Molasse were thrust northward (e.g. Becker, 2000;Burkhard, 1990;Burkhard and Sommaruga, 1998;von Hagke et al., 2012). However, it has also been reported that the proximal foreland basin east of the easternmost ECM, i.e. the Aar Massif, has been subject to post-12 Ma thrusting and horizontal shortening ( Fig. 6a and b; Ortner et al., 2015). Our AHe age data from the Subalpine Molasse (Figs. 3, 4, and 5) fit with AHe ages farther east (von Hagke et al., 2012(von Hagke et al., , 2014b and chronicle a period of thrusting and exhumation of the Subalpine Molasse between 12 and 4 Ma (Fig. 6a). This occurred coevally with the main deformation phase in the Jura FTB, which lasted from ca. 12 to 4 Ma (e.g. Becker, 2000). Similar ages constrained from geological and seismic interpretation, and observations of growth strata in the youngest preserved sediments have been reported for the Subalpine Molasse between Lake Constance and Salzburg ( Fig. 6a; Ortner et al., 2015). Late-Miocene thrusting in the Subalpine Molasse of Central Switzerland has also been inferred from stratigraphic data (Haus, 1935). Since the youngest AHe ages are recorded from the internal tectonic slices of the Subalpine Molasse (Figs. 3, 4, and 5), we can infer the occurrence of break-back thrusting (i.e. back-stepping sequence of thrusting), a characteristic feature, which has been confirmed so far along the Subalpine Molasse between Bern and Salzburg based on thermochronological data and cross section restorations (von Hagke et al., 2012(von Hagke et al., , 2014bOrtner et al., 2015;Schuller et al., 2015) but has been locally argued for as early as the 1930s (Haus, 1935(Haus, , 1937. The break-back thrusts are supposedly younger than the development of the frontal triangle zone, which formed the active northern deformation front from ca. 20 to 12 Ma (von Hagke et al., 2014b;Ortner et al., 2015). Furthermore, our AHe ages from the Plateau Molasse record a partially exhumed PRZ (Figs. 3,4,and 5) and thus corroborate the occurrence of substantial exhumation of the flat-lying Plateau and Foreland Molasse (Cederbom et al., 2011;Genser et al., 2007;Gusterhuber et al., 2012;von Hagke et al., 2012;Zweigel et al., 1998), indicating a largescale exhumation signal across the entire basin.
Although foreland deformation is kinematically linked to the late stage of thrust-dominated exhumation of the ECMs (horizontal tectonics; Herwegh et al., 2017Herwegh et al., , 2020, the occurrence of late-Miocene thrusting in the foreland is not restricted to the areas in front of the ECMs. This can be observed east of the Aar Massif, where the basement topography in the hinterland decreases rapidly towards the east (Fig. 6c), while horizontal shortening within the foreland remains constant at ca. 20-18 km (Fig. 6b). Furthermore, we observe that short-wavelength along-strike variations in the hinterland's basement topography, i.e. the different levels of exhumed basement blocks, are not reflected in the observed amounts of horizontal shortening within the foreland. However, looking at the pattern over a longer wavelength (ca. 500 km, i.e. the length of the northern Central Alps), the general disappearance of the ECMs towards the east and the decrease in foreland shortening are in good correlation (Fig. 6). Combining these observations with the aforementioned data about the timing and kinematics of ECM exhumation and foreland deformation, we suggest that late-Miocene thrusting is a long-wavelength feature which occurs along the entire northern Central Alps, encompassing both the ECMs and the foreland (i.e. Subalpine Molasse and Jura FTB; Fig. 8c). Along the western portion of the Central Alps, this follows an early to mid-Miocene stage, which was characterized by the subvertical extrusion of the eastern ECMs and a stationary deformation front in the Subalpine Molasse. Hence, along the northern front of the western Central Alps, the mid-late-Miocene boundary marks the transition from an inferred buoyancy-driven regime with the subvertical rise of the eastern ECMs (vertical tectonics; Herwegh et al., 2017Herwegh et al., , 2020 to a large-scale thrust-dominated regime (horizontal tectonics) encompassing also the northern front of the eastern Central Alps (i.e. east of the ECMs) thereafter.

Possible link to geodynamic processes beneath the core of the Central Alps
After the early stages of post-35 Ma continent-continent collision, plate convergence rates seemed to decrease noticeably at ca. 20 Ma ( Fig. 8b; Handy et al., 2010;Schmid et al., 1996). At this stage, pro-wedge widening of the Molasse Basin came to a relative halt (Fig. 8a). In the Central Alps, the proximal part of the basin kept subsiding by an additional 2-3 km (Burkhard and Sommaruga, 1998;Schlunegger and Kissling, 2015), while the distal realm became subject to erosion (Kuhlemann and Kempf, 2002). This phase was also associated with the period of differential vertical uplift of the Aar Massif (and probably also the Mont Blanc and Aiguilles Rouges massifs), i.e. its rise along steeply dipping reverse shear zones between ca. 20 and 12 Ma (Herwegh et al., , 2020Wehrens et al., 2017), while the northern deformation front remained stationary in the Subalpine Molasse ( Fig. 8c; Burkhard and Sommaruga, 1998;von Hagke et al., 2014b;Ortner et al., 2015). East of Munich, however, the Molasse Basin experienced, at this time, a period of uniform subsidence and even a short-lived phase of uplift (ca. 17-16 Ma) as evidenced by horizontal to 5 • westward tilting strata of post 20 Ma molasse sediments (Gusterhuber et al., 2012;Zweigel et al., 1998). The latter was attributed to eastward-directed lateral extrusion of crustal blocks (Gusterhuber et al., 2012). It also correlates with the observed westward increasing gradient of crustal thickening associated with an eastward decrease in the mean elevation of the Eastern Alps (Rosenberg et al., 2018). Alternatively, Handy et al. (2015) explained this signal as a result of deep crustal unloading due to a proposed slab tear and a corresponding subduction polarity reversal beneath the Eastern Alps, which has been proposed to have occurred at ca. 20 Ma (Handy et al., 2015;Schmid et al., 2013;Ustaszewski et al., 2008). These observations imply that important along-strike changes and a large-scale geodynamic (Hetényi et al., 2018b;Kissling et al., 2006;Lippitsch et al., 2003;Mitterbauer et al., 2011) and tectonic (e.g. Handy et al., 2015;Ratschbacher et al., 1991;Rosenberg et al., 2018) decoupling between the Central and the Eastern Alps and their Molasse Basins were established in Miocene times. We suggest that this is also mirrored by the along-strike pattern in late-Miocene cumulative horizontal shortening in the north Alpine foreland, which, without fully considering the possibly large uncertainties of shortening estimation in the Subalpine Molasse, seems to decrease from ca. 30 km near Lake Geneva to ca. 20 km near Lake Constance. A decrease from ca. 18 to < 1 km is finally recorded east of ca. 10.5 • E ( Fig. 6). East of Salzburg at the front of the Eastern Alps, zero late-Miocene shortening is recorded in the proximal molasse (Beidinger and Decker, 2014;Hinsch, 2013;Ortner et al., 2015). This pattern has been attributed to an increasing transfer of shortening towards the internal parts of the orogen (i.e. the Tauern Window) and to the out-of-section removal of crust through lateral extrusion in the Eastern Alps (Ortner et al., 2006(Ortner et al., , 2015. While such upper-crustal processes may certainly mask the signal of late-Miocene foreland deformation, the latter may also reflect the response to the change in a geodynamic driving force operating at a larger scale, situated at deeper crustal levels, and encompassing the entire Central Alps (Fig. 9).
In the context of the aforementioned observations, we propose that the attached European lithospheric slab beneath the Central Alps presents a possible candidate for driving the inferred long-wavelength signal of late-Miocene foreland deformation (Fig. 9). This hypothesis is based on the following main observations and considerations. (i) The alongstrike extent of late-Miocene thrusting in the foreland correlates remarkably well spatially with the proposed extent of the steeply south-dipping Central Alpine European slab imaged by seismic tomography (Fig. 9; Kästle et al., 2020;Lippitsch et al., 2003;Zhao et al., 2016). (ii) The long spatial wavelength of tectonically driven exhumation of the Subalpine Molasse, which has been attributed to distinct thrusting events between ca. 12 and 4 Ma, can be viewed as an upper-crustal expression of a lithospheric-scale tectonic driver acting at that wavelength. (iii) The segmentation of the deep structure of the Central and Eastern Alps (Hetényi et al., 2018b;Kästle et al., 2020;Kissling et al., 2006;Lippitsch et al., 2003;Mitterbauer et al., 2011) has been conjectured to have occurred at ca. 20 Ma and induced a geodynamic and tectonic reorganization along the Alpine chain thereafter (Handy et al., 2015;Schmid et al., 2013;Ustaszewski et al., 2008). Hence, this time constraint allows us to discuss the Central Alpine late-Miocene foreland deformation (which is decoupled from the Eastern Alps) in the context of the ongoing evolution of the segmented Central Alpine slab. The corresponding transition in the macro-tectonic regime of the northern Central Alps from inferred buoyancy-driven vertical to large-scale horizontal tectonics may reflect the response to the decreasing rates of European slab rollback and a late phase of post-collisional indentation of the Adriatic plate as recently proposed by Herwegh et al. (2017Herwegh et al. ( , 2020.

The exceptional position of the Bavarian Subalpine Molasse
The configuration of the lithospheric mantle slabs is inherently different beneath the Central and the Eastern Alps (Kästle et al., 2020;Kissling et al., 2006;Kissling and Schlunegger, 2018;Lippitsch et al., 2003;Mitterbauer et al., 2011;Zhao et al., 2016). At the lithospheric scale, a polarity reversal between the European slab and the Adriatic slab has been  Schlunegger and Kissling, 2015). The location of the section is given in Fig. 1. (b) Rates of plate convergence between the Adriatic plate and Europe (Handy et al., 2010;Schmid et al., 1996). (c) Tectonic evolution of the major tectonic units of the northern Central Alps. FTB, fold-and-thrust belt.
conjectured beneath the western Tauern Window (Lippitsch et al., 2003). However, this model is currently debated as new geophysical data indicate that at depth the southward-dipping European slab extends from the Central Alps until east of the Giudicarie and Brenner Faults (Kästle et al., 2020) and possibly as far east as the central Tauern Window (Figs. 1 and 9; Qorbani et al., 2015). Based on tomographic results from the AlpArray initiative, Hetényi et al. (2018b) propose that at ca. 13 • E a northward-dipping slab of the Adriatic plate is present. All these studies show that the deep structure remains very uncertain and debated since tomographic data have been interpreted in different ways regarding the slab geometries at depth (Hetényi et al., 2018b;Kästle et al., 2020;Lippitsch et al., 2003;Mitterbauer et al., 2011;Zhao et al., 2016). However, they all concur on the observation that between the deep velocity anomalies of the Central and the Eastern Alps a major discontinuity is present east of the Giudicarie and Brenner Faults (see also Handy et al., 2015). This correlates with the eastward termination of the late-Miocene Subalpine Molasse near Salzburg (Fig. 9). We thus propose that the geophysical data place the link between the deep structure and deformation of the Subalpine Molasse between Lake Constance and Salzburg in a Central Alpine rather than an Eastern Alpine context. While at the mantle scale, a segmentation of the slab structure is observed between ca. 11 and 13 • E (Fig. 9), an along-strike segmentation at crustal levels occurs at the Brenner Fault. East of the Brenner Fault (Figs. 1 and 9), the post-collisional evolution of the Eastern Alps is characterized by the northward indentation of the Dolomite indenter and the related eastward lateral extrusion of crustal blocks (Ratschbacher et al., 1991). At the larger scale, this was possibly facilitated by slab retreat beneath the Carpathians and the associated rifting in the Pannonian Basin (Peresson and Decker, 1997). The main phase of lateral extrusion occurred in early and mid-Miocene times (Frisch et al., 1998;Ratschbacher et al., 1991), possibly extending into the late Miocene (Ortner et al., 2015). These processes were additionally associated with the collisional exhumation of the Tauern Window and normal faulting along its bounding lowangle normal faults (see Rosenberg et al., 2018). In this respect, the Bavarian Subalpine Molasse is particularly interesting, since it extends over this transition area (Fig. 9). We suggest that in this transient position, the tectonics of this segment of the Subalpine Molasse between Lake Constance and Salzburg was probably affected by the deep-seated dynamics of the Central Alpine European slab, while at the same time also being masked by the aforementioned uppercrustal processes of Eastern Alpine tectonics (e.g. lateral extrusion).
In order to resolve the influence of deep-seated processes on the tectonics of the Subalpine Molasse and the foreland basin in general, we propose that the present-day slab geometries underneath the entire Alps must be resolved at higher resolution. Furthermore, the time evolution of the slabs must be constrained with a focus on if, when, and how a potential subduction polarity reversal occurred in the Eastern Alps. These studies should be supported by source-to-sink analyses linking the stratigraphy of the foreland to hinterland processes. New thermochronological data from the Bavarian Subalpine and Foreland (i.e. undeformed) Molasse will be essential. Furthermore, we expect that ongoing seismotomographic investigations will disclose further details to constrain the deep-seated driving mechanisms (Hetényi et al., 2018a).

Conclusions
In this paper, we presented new low-temperature thermochronological age data from the Subalpine Molasse of the Central European Alps. By comparing our results to published age and stratigraphic data along-strike the Alps from Lake Geneva to Salzburg we conclude the following.
-(U−Th−Sm)/He ages along the Subalpine Molasse of the Central Alps record thrust-related exhumation between 12 and 4 Ma with at least two thrusting pulses associated with it at ca. 10 and 6 Ma.
-The pattern of strain accommodation is strongly conditioned by the local-scale mechanical stratigraphy, since the locus of deformation depends on the distribution of mechanically weak (sand-and mudstone) and strong (conglomerate) lithologies. Accordingly, a general widening of the foreland thrust belt (i.e. tectonic salients) and en bloc deformation of large conglomerate bodies occurs at the locations where the lithology is conglomerate dominated, i.e. at the site of former megafan dispersal systems depositing large and kilometer-thick conglomerate sequences into the foreland. Conversely, tectonic recesses and a more distributed deformation pattern developed where the mechanical stratigraphy is dominated by sand-and mudstone deposits.
Hence, we observe that the deformation style during late-Miocene thrusting of the Subalpine Molasse is masked by local variations in the stratigraphic architecture. However, in terms of timing and kinematics, thrust-related shortening, though decreasing in intensity from the west to the east, is consistent along-strike in those parts of the Alps which correlate with the lateral (i.e. along-strike) extent of the Central Alpine lithospheric mantle slab at depth. Thus, we may interpret late-Miocene long-wavelength thrusting as an uppercrustal signal resulting from this lithospheric driving force. The latter possibly leads to a change in the macro-tectonic regime in the western Central Alps at ca. 12 Ma, from vertical to thrusting-dominated horizontal tectonics.
Finally, although we lack the required data to precisely determine the geodynamic processes responsible for the late phase of shortening in the molasse, we are able to constrain the timing of this event to 12-4 Ma. In addition, this work shows that low-temperature thermochronological data yield an improved understanding of the chronology of orogenic processes where the late orogenic stage of a mountain belt may be characterized by a complex pattern of strain release conditioned by site-specific stratigraphic and thus lithological conditions at the local scale and by widespread tectonism possibly resulting from changes in deep-seated driving forces.
Data availability. The research data are included in this paper and in the Supplement and can be freely accessed.
Author contributions. SM designed the study with support from MH, FS, and CvH. FS and MH assisted SM during sampling in the field. SM carried out mineral separation and picked the apatite crystals. ID carried out the helium extraction and the ICP-MS measurements at the University of Göttingen. CvH and ID assisted in analyzing and interpreting the apatite (U−Th−Sm)/He ages. SM prepared the paper, with contributions from all co-authors.