Late Alpine multistage exhumation of the northwestern Rhodope Metamorphic Complex (northern Rila Mountains, Bulgaria)

The Rhodope crystalline massif is an Alpine metamorphic complex exposed across several mountain ranges in southern Bulgaria and northern Greece which has experienced a complex history including continental collision, partial subduction and syn-metamorphic nappe stacking followed by syn- to post-contractional extension. We present new 40Ar/39Ar and fission-track data from samples taken from both sides of the North Rhodopean Detachment that were combined with detailed structural studies to investigate the tectonothermal evolution of the northern Rila Mountains. A migmatite from the hanging wall of the North Rhodopean Detachment yields a 40Ar/39Ar muscovite age of 100.79 ± 0.55 Ma, a zircon fission-track age of 38.6 ± 1.9 Ma, and an apatite fission-track age of 21.4 ± 1.5 Ma. A biotite schist from the footwall of the detachment yields 40Ar/39Ar biotite age of 34.90 ± 0.15 Ma, and zircon and apatite fission-track ages of 35.6 ± 5.6 and 13.3 ± 1.1 Ma, respectively. Our new data give evidence of a multistage exhumation of the study area. Late Early Cretaceous (~ 101 ± 0.6 Ma) cooling of the Variscan high-grade metamorphic basement through 440–400 °C was caused by either erosion of the emplacing thrust sheet, or post-contractional denudation. Fast exhumation along the North Rhodopean Extensional System drove a pulse of increased tectonic denudation and cooling during the Eocene (39–35 Ma). Exhumation of the rocks in the northern part of the Rila Mountains below 110 ± 10 °C during the middle–late Miocene was associated with displacement along a system of normal faults.


Introduction
The Rhodope Metamorphic Complex (RMC) is exposed in several mountain ranges in southern Bulgaria and northern Greece, including the Rhodope, Rila, Pirin, Vrontous and Pangaion mountains (Fig. 1). In general, the geology of the Rhodope Mountains has attracted far more attention than some other parts of the complex, e.g. the Pirin and Rila Mountains. One probable reason for this is simply accessibility. The Pirin and Rila Mountains, culminating at 2925 m, are the highest and least accessible parts of the Rhodope Metamorphic Complex (Fig. 2). Extensive research in the Rhodope Mountains during the previous 30 years reveals a complex Alpine age (Late Jurassic to Quaternary) tectonic evolution that includes continental collision, partial subduction and syn-to post-contractional crustal extension (e.g. Burg et al. 1995Burg et al. , 1996Dinter and Royden 1993;Dinter 1998;Ricou et al. 1998;Burg 2012;Froitzheim et al. 2014). In comparison, geochronological and structural studies of the rocks exposed in the Rila and Pirin mountains are relatively scarce and only partly elucidates the tectonic evolution of the northwestern Rhodope Metamorphic Complex (Dimov and Damianova 1996;Shipkova and Ivanov 2000;Tueckmantel et al. 2008;Gunnell et al. 2017;Gorinova et al. 2019). The results of this limited research indicate at least two major Alpine contractional events (Gorinova et al. 2019), followed by a period of extension, which drove exhumation along a detachment fault (Shipkova and Ivanov 2000;Tueckmantel et al. 2008;Gorinova et al. 2019). In the Rila Mountains, rocks of the Lower Terrane are not exhumed, unlike the central, east and southern parts of the Rhodope Metamorphic Complex (Fig. 1, Gorinova et al. 2019). Furthermore, the rocks exposed in the Rila Mountains exhibit no evidence of syn-extensional migmatisation, characteristic for other parts of the metamorphic complex      (Gorinova et al. 2019), suggesting the area probably escaped these high-temperature/high-pressure Alpine events. This paper presents new structural and thermochronological (fission-track and 40 Ar/ 39 Ar) data that constrain the late Alpine tectonic evolution of the northwestern Rhodope Metamorphic Complex exposed in the northern parts of the Rila Mountains. Integration of the new data with other previously published data reveal a multistage extensional evolution of the region comprising at least two major phases of exhumation. An initial early phase of very fast exhumation along a ductile to brittle crustal-scale detachment, the North Rhodopean Detachment (Fig. 1), during the late Eocene . A second, later phase of exhumation, beginning in the middle to late Miocene and continuing even in the Quaternary, is accommodated by displacement along brittle, high-angle normal faults that are responsible for the final uplift and creation of the present-day Rila and Pirin highmountain ranges. These two episodes of major extension are separated by an intervening subordinate phase of strike-slip to thrust tectonics, which is probably responsible for some reheating that is observed in all the new rock samples analysed herein as well as those analysed in previous studies.

Kalin Pluton Kalin Pluton Kalin Pluton D o d o v V r a h S Z D o d o v V r a h S Z D o d o v V r a h S Z B is t r it s a S Z B is t r it s a S Z B is t r it s a S Z
The multidisciplinary, structural and thermochronological, analytical methodology this paper employs provide additional degrees of constraint for interpretation of the Cenozoic tectonic evolution of the area, which facilitates determination of the age, distinct structural style of deformation and likely plate tectonic setting of each discrete episode of exhumation and uplift or burial. Analyses of the first phase of Late Eocene age exhumation of middle to lower crustal rocks associated with the formation of crustal-scale, low-angle, ductile to brittle detachment faults accommodating regional extension, suggested deformation in a subduction-related plate tectonic setting. Whilst the second phase of exhumation and uplift, creating the present-day high-mountain topographic relief, associated with the formation of brittle, high-angle normal faults and continental basins is most likely associated with a change in tectonic setting from orogenic subduction to a post-orogenic crustal thinning, inferring a reorganisation of the plate boundaries at this time.

Rhodope metamorphic complex
The Rhodope high-grade metamorphic complex, exposed in southern Bulgaria and northern Greece, represents an Alpine syn-metamorphic nappe stack that experienced massive crustal extension and core complex formation during the Cenozoic (Fig. 1). Although several different tectonometamorphic subdivisions were suggested (e.g. Burg et al. 1996;Turpaud and Reischmann 2010;Janák et al. 2011;Burg 2012;Kydonakis et al. 2015;Schmid et al. 2020), there is general consensus that the metamorphic complex is composed of three main terranes that were initially juxtaposed during a Cretaceous (e.g. Ricou et al. 1998;Burg 2012;Bonev et al. 2015) or Cenozoic collision (e.g. Dinter 1998;Nagel et al. 2011;Froitzheim et al. 2014;Gorinova et al. 2019;Fig. 1). A Lower and Upper terrane represent continental fragments separated by an Intermediate terrane, which has a Jurassic ophiolitic and magmatic arc protolith (e.g. Burg 2012). All of these Rhodopean terranes contain ultra-to high-grade/high-pressure metamorphic rocks, some of which are reported to be of Alpine ages (Liati 2005;Liati et al. 2002;Bosse et al. 2010;Georgieva et al. 2011;Kirchenbaur et al. 2012;Kydonakis et al. 2015;Petrik et al. 2016;Collings et al. 2016). Two differing models were proposed concerning the time of the major orogeny leading to the tectonic juxtaposition of the Rhodopean terranes (e.g. Ricou et al. 1998;Dinter 1998;Nagel et al. 2011;Burg 2012;Froitzheim et al. 2014;Bonev et al. 2015;Gorinova et al. 2019). The first one suggests a single, Late Jurassic-Early Cretaceous, major orogenic episode of northward subduction (e.g. Burg 2012) whereas, the second envisages a two-step scenario, comprising an initial north-vergent Late Jurassic arc-continent collision, followed by Eocene-Oligocene south-vergent continent-continent collision (e.g. Froitzheim et al. 2014;Petrik et al. 2016). The model presented by Gautier et al. (2017) also distinguishes two orogenic cycles in the evolution of the Rhodope metamorphic complex, one of Late Jurassic-Early Cretaceous age and the other of Late Cretaceous-Cenozoic time age, which are separated by ∼45 Myr (117-70 Ma) period of relative tectonic quiescence.
Eocene-Oligocene syn-to post-contractional extension followed by Miocene post-contractional extension generated core complexes that are bound by gently dipping detachment faults (e.g. Dinter and Royden 1993;Krohe and Mposkos 2002;Bonev et al. 2006;Brun and Sokoutis 2007;Georgiev et al. 2010Georgiev et al. , 2016Kounov et al. 2015Kounov et al. , 2020 with cores mainly comprising rocks of the Lower Terrane (Fig. 1). Different names were attributed to these extensional core complexes, whereas we group them here in three major metamorphic complexes-Central, Southern and Eastern (Fig. 1, Kounov et al. 2020).
The extension was also associated with formation of grabens that are filled with Eocene to Quaternary marine to continental deposits (e.g. Cernjavska 1977;Ivanov et al. 2000), accompanied by volcanic and intrusive magmatism (e.g. Harkovska et al. 1998;Marchev et al. 2005). Ductile southward syn-extensional thrusting along the Chepelare and Nestos shear zones (Fig. 1) as young as 36 Ma was also reported, suggesting that the deeper crustal levels of the Rhodope Metamorphic Complex were still accommodating shortening until at least this time (Bosse et al. 2009;Jahn-Awe et al. 2012;Nagel et al. 2011;Gautier et al. 2017).

Rila Mountains
For the Rila Mountains, which occupy the northwestern corner of the Rhodope Metamorphic Complex (Figs. 1 and 2), several slightly diverging tectonometamorphic subdivisions were suggested (Dimov and Damianova 1996;Sarov et al. 2011a;Gorinova et al. 2019). However, there is common agreement that all the units have experienced a major tectonometamorphic event during Alpine time, except the uppermost Variscan high-grade metamorphic rocks (Fig. 2, Sarov et al. 2011a). The main units exposed in the northwestern part of Rila Mountains are Kabul, Malyovitsa and Preka Reka units, whilst the multiphase Rila-West Rhodope Batholith occupies the central and eastern parts (71-39 Ma, von Quadt and Peytcheva 2005;Peytcheva et al. 2007;Figs. 1 and 2). The Chepino Unit, which is previously considered to be part of the Intermediate (Sarov et al. 2008a, Fig. 2) or the Lower Rhodopean terranes (Burg 2012), crops out in the easternmost margin of the Rila Mountains where it is overlain by the North Rhodopean Unit of the Upper Rhodopean Terrane (Burg 2012, Fig. 2) above the North Rhodopean Detachment . Below, the units exposed in the study area and the composition of the North Rhodopean Extensional System ) are described in more detail.

Chepino Unit
The exact structural position of the unit in the Rhodopean crustal stack is unclear. However, Maastrichtian-early Paleocene (63-59 Ma) monazite ages from a leucosome , giving evidence for Alpine migmatisation, suggests a lowermost tectonic position in the Rila Mountains. The unit is composed of migmatized biotite and amphibole-biotite gneisses, thin marble packages and lenses of amphibolites, all of which are affected by upper amphibolite-facies metamorphism at temperatures of 600-650 °C associated with partial melting (Sarov et al. 2008a).
The metamorphic foliation in the study area dips shallowly to the northeast and east-northeast, with stretching lineations that plunge towards the north-northeast (Sarov et al. 2008a). Sense-of-shear criteria record both top-to-the-south and southeast, as well as top-to-the-north and northwest syn-metamorphic tectonic transport. Interpretation of these opposite directed senses-of-shear infer a complex tectonic evolution, which includes initial south-vergent thrusting followed by northward directed extension and the formation of the North Rhodopean Extensional System during the late Eocene .

Preka Reka and Malyovitsa Units
The relationship between the Malyovitsa and Preka Reka units with the Chepino Unit is unclear, because they are separated by the Rila-West Rhodope Batholith (RWRB, Fig. 2). The Preka Reka Unit (Gorinova et al. 2019) mainly consists of orthogneisses and the overlying Malyovitsa Unit comprise paragneisses, orthogneisses, amphibolites, marbles, schists and metaultrabasic rocks (e.g. Dimov and Damianova 1996) (Fig. 2). Dimov and Damianova (1996) reported top-to-thesoutheast displacement along the ductile Malyovitsa Shear Zone that separates both units. Gorinova et al. (2019) consider both units are part of the same litho-tectonic body (the "Malyovitsa Unit"), even though they are separated by a ductile shear zone, because granites of the Preka Reka Unit intrude the Malyovitsa Unit para-series. Gorinova et al. (2019) also interpret the Malyovitsa Unit is regionally very similar to Sarnitsa Unit, which outcrops roughly 100 km south-southeast, to the south of the RWRB (Sarov et al. 2010, Fig. 1 Fig. 2). The rocks of the Malyovitsa and Preka Reka units underwent moderate amphibolitefacies metamorphism with no evidence of a high-pressure event (Gorinova et al. 2019). Gorinova et al. (2019) described an older metamorphic fabric, related to an early Late Cretaceous tectonometamorphic event, which was isoclinally folded during the development of the axial planar foliation that is the main penetrative fabric in both units. A U-Pb concordia age of ca. 82 Ma obtained from recrystallised rims of zircon crystals from the Jurassic orthogneiss body (Gorinova et al. 2019) is interpreted to date the earlier amphibolite-facies metamorphic event. Based on syn-kinematic granitic dykes and veins, the formation of the second main fabric in the Malyovitsa Unit is fixed between 51 and 48 Ma, contemporaneous to the movement along the Dodov Vrah Shear Zone (Fig. 2, Gorinova et al. 2019). This foliation dips moderately to the southwest or west and bears a subhorizontal northwest-southeast-trending stretching lineation that is parallel to the fold axes of the isoclinal folds. Kinematic indicators reveal a syn-metamorphic, top-to-the-southeast sense-of-shear (Dimov and Damianova 1996;Shipkova 1999;Gorinova et al. 2019).
The main penetrative foliation in the Kabul Unit dips to the southwest or northeast. Fold axes plunge at low angles to the northwest or southeast. The foliation bears mostly a northwest-southeast-trending stretching lineation. Kinematic indicators in the non-folded domains show a top-tothe-north-northwest sense-of-shear (Gorinova et al. 2019). This prevalent deformation formed under amphibolite-facies conditions during earliest Late Cretaceous (100-90 Ma), syn-metamorphic, northwestward thrusting of the Variscan high-grade metamorphic rocks over the Kabul Unit (Gorinova 2020). The fabric is sealed by the 60-58 Ma old Kalin pluton (Zagorchev et al. 2014;Gorinova et al. 2019;Fig. 2). Later, between 58 and 48 Ma, the Kabul Unit was thrust over the Malyovitsa Unit along the Dodov Vrah Shear Zone (Gorinova 2020, Fig. 2). This shear zone hosts deformed diorites to gabbro-diorite bodies with a protolith age of ca. 76 Ma (U-Pb zircon ages, Gorinova et al. 2019).
The North Rhodopean Unit consists of biotite, muscovite and two-mica gneisses, marbles and amphibolites, all affected by upper greenschist to low-grade amphibolitefacies metamorphism related to the formation of generally subhorizontal northwest-southeast striking foliation bearing a subhorizontal lineation (Naydenov et al. 2008). Naydenov et al. (2009) report U-Pb zircon mean ages of between 162.9 ± 3.6 and 155.9 ± 4.4 Ma for the youngest deformed magmatic rocks, suggesting an Alpine age for the main metamorphic event in this unit. A Late Jurassic-Early Cretaceous age of the metamorphism is further constrained by Late Cretaceous age granitoids that intrude into the gneisses (Fig. 2, Naydenov et al. 2009). The North Rhodopean Unit is separated from the Chepino Unit by the North Rhodopean Detachment (Fig. 2).

Variscan high-grade metamorphic units
Variscan high-grade metamorphic rocks are reported in the study area both above and below the North Rhodopean Detachment. The rocks west of the Kabul Unit, which were previously attributed to the Ograzhden metamorphic complex (Sarov et al. 2011b, Fig. 2), were later described as the Polich Unit (Gorinova et al. 2019) of the Upper Allochthon of the Rhodope Metamorphic Complex (e.g. Froitzheim et al. 2014). The Variscan rocks below the North Rhodopean Detachment consist of migmatized biotite and twomica gneisses that host bodies of metagabbros (Sarov et al. 2011b). Their lower boundary is the amphibolite-facies Bistritsa Shear Zone, along which the Polich Unit was thrust over the Kabul Unit (Fig. 2). Sense-of-shear criteria indicate top-to-the-northeast (Sarov et al. 2011b) or northwestdirected tectonic transport (Gorinova et al. 2019).
The high-grade metamorphic rocks situated north of the Dzherman and Cherni Iskar rivers, above the North Rhodopean Detachment were equally attributed to the Ograzhden (Sarov et al. 2011a) and later to the Sredna Gora units (Gorinova et al. 2019, Fig. 2). They are made up of migmatized biotite and two-mica gneisses hosting bodies of metagabbros and lenses of serpentinized ultrabasites. Late Ordovician protolith ages are reported for the migmatized gneisses (U-Pb zircon ages between 454 and 446 Ma;Peytcheva et al. 2009). There are no studies for the metamorphic evolution of these rocks. However, similar metamorphic rocks farther north in the basement of the Sredna Gora Zone are intruded by undeformed Carboniferous to Permian granites, suggesting a Variscan age of the regional migmatisation and metamorphism (Carrigan et al. 2005(Carrigan et al. , 2006Gerdjikov et al. 2010). Shipkova and Ivanov (2000) were the first to describe a Cenozoic detachment fault (the Dzherman Detachment) that borders the Rila Mountains to the west and northwest and underlies the late Eocene-Oligocene Padala basin (Fig. 2, Cernjavska 2000). Later, Sarov et al. (2011a, b) describe another detachment located between Sapareva Banya and Govedartsi ( Fig. 2) named the Lakatitsa Detachment. Further east, in the area of Sestrimo (Fig. 2), the brittle-ductile shear zone separating the Chepino Unit from North Rhodopean Unit is termed the North Rhodopean Detachment (Naydenov et al. 2008). More recently, Gorinova et al. (2019) describe the whole extensional structure that borders the northwestern Rila Mountains as the North Rila Detachment, up to where it is cut by the brittle, high-angle, Saparevo normal fault along the northwestern slope of the Rila Mountains (Fig. 2).

North Rhodopean Extensional System
In the study area, Gorinova et al. (2019) describe a lower grade, north-dipping mylonitic foliation overprinting the older fabric along the North Rila Detachment. The generally shallow-dipping fabric of this zone outcrops in a large area along the western bank of the Dzherman River and also on the mountain slopes south of Govedartsi (Fig. 2). The kinematic indicators show a consistent top-to-the-north-northwest sense-of-shear and the age of displacement is latest Eocene to early Oligocene (Gorinova et al. 2019). Large parts of the detachment zone were offset by normal faults formed during later extensional tectonics.
In this study, the extensional shear zones together with the brittle to brittle-ductile North Rhodopean Detachment, developed along the northern slopes of the Rila Mountains, are referred to as the North Rhodopean Extensional System, following .

Rila-West Rhodope Batholith (RWRB)
The batholith is built of a succession of at least four granitic magmatic pulses intruded during the Maastrichtian to the late Eocene time (Valkov et al. 1989;von Quadt and Peytcheva 2005;Peytcheva et al. 2007). The first phase consists of granodioritic and quartz-dioritic rocks that form the largest body in the area, called the Belmeken body emplaced at 70.77 ± 0.09 Ma (U-Pb zircon age) (Valkov et al. 1989;Peytcheva et al. 2007; Fig. 2). These granodiorites, which have a calc-alkaline affinity with a mixed (mantle-crustal) magma source, were formed in a subduction-related setting (Kamenov et al. 1999).
The most widespread lithologies within the batholith are the medium-to coarse-grained biotite granites intruded during the second phase. In the study area, this phase is represented by the Musala body and a small part of the Western Rhodope body in the eastern part of the area, as well as the muscovite-biotite Shpanyovitsa body (Fig. 2). There are no geochronological data for the Musala body, but the Shpanyovitsa granite has a U-Pb zircon age of 39.39 ± 0.21 Ma (Peytcheva et al. 2007). The third phase consists of several small, fine-grained granite to plagiogranite bodies, such as the Chavcha (Fig. 2). The granites of the second and third phase have a high K calc-alkaline affinity with isotopic compositions diagnostic of a mixed mantle-crustal origin. However, their crustal component is considerably higher in comparison to the Late Cretaceous granodiorites of the first phase (Kamenov et al. 1999). Therefore, the granites of the second and third phase are considered to have formed in late-to post-collisional tectonic settings (Kamenov et al. 1999). Stocks and veins of isotropic aplitoid and pegmatoid granites form the fourth phase (Fig. 2). The emplacement of the RWRB led to the formation of the north-northwest to south-southeast elongated Rila-West Rhodope dome (e.g. Ivanov 2017), where the granitoids form the core of the dome, with a rim of high-grade metamorphic rocks (Fig. 1).

Structural studies along the North Rhodopean Extensional System
As defined in this paper, the North Rhodopean Extensional System (NRES) comprises the brittle-ductile to brittle low-angle North Rhodopean Detachment that is underlain by a relatively narrow ductile mylonitic shear zone and a wide protomylonitic zone. Below, is a description of the different elements of the extensional system that outcrop along the northern slopes of the Rila Mountains between the towns of Dupnitsa and Sestrimo (Fig. 2). This system was previously traced as far east as the area of Krichim   (Fig. 1), but its relationship with the detachment system of the Central Rhodope Metamorphic Complex remains unclear ).

Brittle North Rhodopean Detachment
Unfortunately, exposures of this detachment zone in the study area are limited, because it is either buried below lower Oligocene sedimentary rocks and the North Rhodopean Thrust or offset by Miocene-Quaternary normal faults (Fig. 2). Sometimes its existence is only deduced by the presence of brecciated and silicified zones in the footwall.

Sestrimo area
South to southeast of Sestrimo the brittle detachment can be followed (Fig. 3), along the contact between the Chepino and North Rhodopean units. In the rest of the area, it is either buried below the North Rhodopean Thrust or is partially offset by the Miocene-Quaternary Sestrimo normal fault (Valkova and Spiridonov 1979) (Fig. 3). The detachment zone is best studied along the Sestrimska River valley, where the uppermost part of the footwall and tectonites of the detachment zone crop out (Figs. 3, 5a). The fault zone consists of 3 to 5 m thick, black ultra-cataclasites that overlay ultramylonitic gneisses, marbles and pegmatites of the Chepino Unit (Fig. 5a). The direct contact with the hanging wall, comprising strongly faulted gneisses, amphibolites, schists and marbles of the North Rhodopean Unit, is covered by colluvium.

Govedartsi area
In the Govedartsi area, the detachment zone is partly reactivated and/or offset by the Govedartsi normal fault (Fig. 4). The zone is easily mappable by following the lithologically contrasting boundary between the leucocratic orthogneisses of the Preka Reka Unit in the footwall, and the dark Variscan high-grade gneisses and amphibolites in the hanging 1 3 wall (Fig. 4). Brecciated and silicified rocks from the uppermost structural parts of the footwall mark the zone at several localities in the vicinity of Govedartsi. Along the road to Mechit hut (Fig. 4), black ultra-cataclasites crop out above the brecciated and silicified orthogneisses (Fig. 5b).

Dzherman River valley
A wide brittle fault zone that separates the Malyovitsa Unit from the Variscan high-grade metamorphic rocks along the Dzherman River valley has previously been described as a thrust (Marinova and Zagorchev 1991) and a detachment zone (Lakatitsa Detachment of Sarov et al. 2011a, or North Rila Detachment of Gorinova et al. (2019). The brittle zone often has a steep dip due to its reactivation and tilting along the Dzherman normal fault (Figs. 4,5c). To the north, the fault zone is cut by the Miocene-Quaternary Saparevo normal fault, whereas to the south, it is covered by Oligocene and Quaternary sedimentary rocks (Marinova and Zagorchev 1991, Fig. 4  In the area of the Dzherman River, the Malyovitsa gneisses in the footwall of the detachment are intruded by voluminous leucocratic granitic dykes and pegmatites. All these rocks are strongly cataclastically deformed along a ~ 50 m thick zone (Fig. 5c). The zone is capped by several metres thick grey gauge and/or tectonic breccias. Quaternary slope deposits cover the contact with the Variscan migmatites in the hanging wall of the detachment. An almost complete section across the detachment zone was observed to the northeast of the Gradishte Peak, where the deeply incised Dzherman River cuts through the zone (Figs. 4 and 5c). Here, the footwall comprises alternating biotite gneisses, various schists, amphibolites and rare marble layers intruded by pegmatitic dykes. These rocks display only the older penetrative schistosity S 1 that is conformable to the detachment zone (Fig. 4). There is no evidence of intensive ductile deformation that could be related to the extension. The uppermost ca. 10 m of the footwall are affected by strong brecciation,

2.18°D
z h e r m a n r .    silicification and propylitization, and cut by numerous millimetres to centimetres thick ultra-cataclasite zones and faults. The ultra-cataclasite zones generally dip steeply (50-60°) to the east-northeast and often form an anastomosing network.

Mechit Hut Mechit Hut Mechit Hut Rila Lakes Rila Lakes Hut Hut Rila Lakes Hut
With the exception of some mechanically competent, lowstrain lenses (represented mainly by pegmatites and coarse gneisses), the brecciated rocks display a weakly developed S 2 cataclastic foliation that strikes parallel to the contact. Extensional, top-to-the-northeast kinematics can be deduced by the geometric relationship between the shallower S 2 cataclastic foliation and a steeper set of northeast-dipping brittle shear zones (Fig. 6a).

Mylonitic shear zone
A high-strain zone of upper greenschist-to lower amphibolite-facies metamorphic (medium-grade) mylonites (450-600 °C, Passchier and Trouw 2005; Trouw et al. 2009), exposed below the detachment, is best observed at two particular localities in the Sestrimo and Govedartsi areas (Fig. 2). In places, the medium-grade mylonites are overprinted by greenschist-facies metamorphic (lowgrade) mylonites and ultramylonites (300-400 °C, Passchier and Trouw 2005;Trouw et al. 2009). Our observations show that, although differing in relative temperature, both the medium-grade mylonites and the low-grade ultramylonites show a top-to-the-north sense-of-shear and are therefore probably a result of the same deformation event. The relative temperature difference between the mylonites occurs because, as the footwall of the detachment moved upward during crustal extension, the mediumgrade mylonites that formed at depth were displaced to progressively shallower levels, undergoing retrograde deformation and strain localisation causing cataclasis and brecciation (e.g. Lister and Davis 1989). Therefore, we attribute the penetrative ductile fabric associated with the   top-to-the-north sense-of-shear to crustal extension. The absence of a medium-grade mylonitic zone in other parts of the study area may be due to offset by a later brittle detachment and subsequent normal faulting, or because it was never exhumed (e.g. the area of the Dzherman River valley and southwest of Sapareva Banya, Fig. 2).

Sestrimo area
In the Sestrimo area, the mylonitic shear zone is about 100 to 250 m wide (Fig. 5a). It developed on the dark biotite gneisses, amphibolites and migmatites of the Chepino Unit, which are intruded by a variety of dykes, including fine-grained biotite granites, leucogranites and pegmatites, most probably related to the Shpanyovitsa granites (Peytcheva et al. 2007). Most of the granitic dykes are concordant to the mylonitic foliation S 2 in the host rocks and have a pervasive solid-state fabric that in places grades from C/S mylonites to ultramylonites. Lower strains observed in the granitic dykes in comparison to the metamorphic rocks of the Chepino Unit could be due to their higher mechanical competency and/or syn-tectonic emplacement. South of Sestrimo, Von Quadt and Peytcheva (2005) report zircon and monazite single grain U-Pb ages of between 43 and 41 Ma from a deformed granitic vein intruded in the Chepino Unit. The lower contact of the shear zone outcrops along the road from Sestrimo to Chaira (Fig. 3), where it exhibits a sharp shear gradient between the mylonitic host rocks of the Chepino Unit and protomylonitic granodiorites of the RWRB. Previously this shear zone has been described as a part of the Yadenitsa-Grashovo brittle-ductile extensional shear zone, which follows, in places, the contact between the Belmeken granodiorites and the West Rhodopean granites (Naydenov et al. 2008, Fig. 3).
The mylonitic foliation S 2 is principally defined by biotite, muscovite, quartz ribbons and stretched leucosomes and has a generally shallow dip (10-25°) towards the northeast (Fig. 3). However, the orientation of the foliation does vary due to post-mylonitic faulting. Newly formed black hornblende and biotite along the foliation planes, as well as feldspar core and mantle structures and quartz recrystallisation by a subgrain rotation and grain boundary migration mechanisms indicate upper greenschist-to lower amphibolite-facies deformation conditions (ca. 450-600 °C, Stipp et al. 2002;Passchier and Trouw 2005, Fig . 7b, c). Grain-size reduction and the formation of lowgrade ultramylonitic levels are observed within the shear zone as well as densely spaced C' shear bands, which indicate strain localization during the exhumation (Dimov and Georgiev 2000, Figs. 6b and 7a). Isoclinally folded, centimetre-thick quartz-feldspar leucosome bands are often reduced to rootless hinges. Their axes are subparallel to the stretching lineation and axial planes are concordant to the mylonitic foliation S 2 . These folds formed during the folding of an earlier migmatitic fabric S 1 . The mylonites display a prominent north-south trending stretching lineation L 2 that is commonly observed on S 2 foliation and the shear bands (Fig. 3). The lineation is defined by elongated quartz and quartz-feldspar mineral aggregates, the alignment of hornblende crystals and clots of white mica. Kinematic indicators within the medium-grade mylonites and the low-grade ultramylonites include asymmetric, sigma and delta-type K-feldspar porphyroclasts, asymmetric boudinage and C' shear bands (Fig. 6b and c). All these features indicate a general top-to-the-north sense-of-shear, which corroborates previously published data (Dimov and Georgiev 2000;. Within the foliated granitic dykes, a pervasive C/S fabric also shows a top-tothe-north sense-of-shear.  Dam (42.192639 N, 23.885038 E) at the footwall of the North Rhodopean Detachment (Fig. 3). a Hand-sized specimen of mylonitic and ultramylonitic granodiorite with polished surface. C/S fabrics together with C'-type shear bands indicating top-to-the-north sense-of-shear. Due to the high-strain rate the ultramylonitic foliation S um is parallelized to the C-type shear bands. b, c Partially recrystallised K-feldspar porphyroclasts (KFs) within a matrix of dynamically recrystallised by a subgrain rotation and grain boundary migration mechanism quartz indicating upper greenschist to lower amphibolite-facies (medium-grade) deformation conditions (ca. 450-600 °C, Stipp et al. 2002;Passchier and Trouw 2005). d Ultramylonites formed mostly by bulging recrystallization indicating lower greenschistfacies temperature conditions between 300 and 400 °C (Stipp et al. 2002;Passchier and Trouw 2005) overprinting an older upper greenschist to lower amphibolite-facies medium-grade mylonitic C/S fabric both indicating top-to-the-north sense-of-shear. Delta-and sigmatype K-feldspar porphyroclasts (KFs) indicate top-to-the-north senseof-shear. g Mica fish indicating top-to-the-north sense-of-shear

Govedartsi area
The mylonitic shear zone in the Preka Reka Unit orthogneisses is ~ 50 m thick (Figs. 4, 5b). Most outcrops are situated along steep river valleys, because much of the area is covered by Quaternary glacial and alluvial-colluvial deposits. An almost complete section of the shear zone outcrops along the Golyama Lopushnitsa River (Figs. 4). Here, 2-3 m thick, gently north-northwest-dipping ultramylonites are visible (Fig. 5b) immediately below the brittle detachment. The ultramylonitic rocks are massive, homogenous and lack a macroscopically observable foliation fabric. They consist of a black-green, chlorite-rich matrix that hosts K-feldspar porphyroclasts with diameters of ~ 1 cm and rarely, larger pegmatitic clasts. Mylonitic orthogneisses are present structurally below the ultramylonitic zone (Fig. 5b). In fact, the low-grade ultramylonitic zone described above develops on this earlier fabric. The most common lithology is white to pink biotite orthogneiss with up to 2 cm large K-feldspar porphyroclasts (Fig. 6d). More rarely, banded biotite orthogneisses are seen, probably formed at the expense of equigranular granite, which consist of quartz ribbons alternating with a fine-grained matrix containing feldspar porphyroclasts and biotite. Leucogranite and pegmatite dykes concordant to the foliation are also common (Fig. 6d). Most of them are intensely deformed, although some still exhibit a weakly developed foliation. The foliation S 2 in the mylonites dips shallowly (5-20°) to the north-northwest (Fig. 4) and is defined by polycrystalline quartz ribbons recrystallized tails around K-feldspar porphyroclasts and biotite flakes. The stretching lineation L 2 generally plunges towards the north and is defined by aligned mineral aggregates, elongated K-feldspar porphyroclasts and quartz rods (Fig. 4). The foliation planes are often cut by steeper shear bands that are up to 1-2 cm wide and characterised by a pronounced grainsize reduction. The shear bands and asymmetric K-feldspar porphyroclasts indicate a top-to-the-north sense-of-shear (Fig. 6d).

Deformation in the footwall of the NRES below the mylonitic shear zone
Intense ductile deformation related to extensional tectonics is observed in large parts of the footwall below the mylonitic shear zone. The main interpretation challenge is distinguishing younger extension-related fabric from older one formed during Mesozoic to Paleocene crustal shortening. Direct relationships between both fabrics are not observed because they are subparallel. This parallelism could be due to relatively coaxial principal strain axes for both tectonic events or rotation of the older contractional structures. Consequently, the actual extent and character of the extensional deformation can only be evaluated by observations on the Late Cretaceous to late Eocene granites of the RWRB (von Quadt and Peytcheva 2005; Peytcheva et al. 2007;Figs. 3 and 4). Since these granites postdate the main Mesozoic metamorphism in the Rhodopean complex, their deformation can be exclusively attributed to the Cenozoic extension.
Two domains with a differing character of ductile deformation are distinguished in the footwall of the NRES. In the eastern and central parts of the Rila Mountains, east of the Dzherman River valley (Fig. 2), the NRES is better developed and an almost full section of the system is observed, including the brittle low-angle detachment fault and the mylonitic shear zone below it (Fig. 2). The omission of the detachment fault and reduction of the mylonitic shear zones at certain places is due to later thrusting and normal faulting (Figs. 2, 3 and 4). On the other hand, the deformation in the footwall is less intense along the northwest slope of Rila Mountains, west of the Dzherman River valley (Fig. 2) where the upper greenschist-to lower amphibolitefacies mylonitic shear zone, if formed at depth, was never exhumed to the surface. Accordingly, the extension-related deformation in the footwall of the system will be discussed in two separate domains-Eastern and Western.

Eastern domain
In the Eastern domain, to the east of Dzherman River valley, the footwall of the NRES mainly comprises composite RWRB (Fig. 2). Most of the Belmeken granodiorites display a moderate to intensely developed solid-state foliation S 2 , whilst lineations L 2 are less frequent (Fig. 3). The strain intensity varies and protomylonites predominate. The foliation below the mylonitic zone has a shallow dip to the northwest, north or east, whilst in the southwestern part of the study area it has a shallow dip to the southwest (Fig. 3). The Shpanyovitsa and Chavcha bodies, which crop out along the slope of the Rila Mountains to the west of Sestrimo (Figs. 2 and 3), show solid-state deformation. Protomylonites are developed in the northwestern part of the granites, as well as along its northern margin. Locally they grade into mylonitic and ultramylonitic zones that have a thickness of tens of centimetres. In the western part of the granite, the foliation S 2 has a shallow dip to the north-northeast (Fig. 3) and is marked by aligned quartz aggregates and micas. The stretching lineation L 2 plunges towards the north and north-northeast (Fig. 3). C' shear bands indicate a top-to-the-north sense-of-shear (Fig. 6e).
The contact between the Shpanyovitsa granite and the Belmeken granodiorites is well exposed along the road between Sestrimo and Stankovi Baraki Dam (Fig. 3). Here, the granites of the Shpanyovitsa body display a weakly developed southeast-dipping foliation S 2 , whereas a prominent strain gradient is observed towards the contact with the Belmeken granodiorites. The gradient marks a transition from weakly foliated granite to protomylonites and C/S mylonites . Locally the C/S mylonites are overprinted by low-grade ultramylonites. Ultramylonites formed mostly by bulging recrystallization, indicating lower greenschist-facies temperature conditions between 300 and 400 °C (Stipp et al. 2002;Passchier and Trouw 2005; Fig. 7 d). The Belmeken granodiorites are intruded by a large number of leucocratic dykes that often display a solid-state deformational fabric. Numerous shear sense criteria, such as C/S fabric, mica fishes, sigma-and delta-type K-feldspar porphyroclasts indicate a top-to-thenorthwest to north-northwest sense-of-shear (Fig. 7 d, e). The Belmeken granodiorites are cross-cut by a younger and voluminous pulse of coarse-grained biotite granites (Musala body, Valkov et al. 1989, Fig. 2). Field data from the section to the south of Mala Tsarkva (Fig. 4) document a zone of intense solid-state deformation that is several tens of meters thick. The foliation dipping to the north-northwest bears a down-dip stretching lineation (Fig. 4). Shear bands indicate a top-to-the-north-northwest sense-of-shear. Mesocratic, ultramylonitic zones up to 20 cm thick are also observed.
In the westernmost part of the Eastern domain, the Govedartsi area, below the mylonitic zone of the NRES, the foliation in the Jurassic orthogneisses of the Preka Reka Unit dips to the northwest to west (Fig. 4). The direct relationship between this foliation and the extensional mylonitic zone is not observed in this area, and thus it remains unclear if the deformation is related to the extension or the earlier contraction. The shear sense indicators are also not observed. Therefore, only steeply dipping shear bands that intersect the foliation planes and indicate a general top-to-the-north sense-of-shear, are considered to be related to Cenozoic extensional phase.

Western domain
The Western domain extends from the Dzherman River valley to the town of Dupnitsa, along the northwestern slope of the Rila Mountains. Here, the footwall of the NRES consists of amphibolite-facies metamorphic rocks of the Malyovitsa and Kabul units, as well as Variscan high-grade metamorphic rocks (Sarov et al. 2011a;Gorinova et al. 2019;Figs. 2 and 4). The medium-grade mylonitic shear zone is missing in this domain, and the extensional deformation is expressed by low-grade mylonites and extensional shear bands.
In the Sapareva Banya-Panichishte area, the older contraction-related foliation S 1 between the Kalin pluton and the Dzherman River valley generally trends north-northwest to south-southeast and dips in both directions (Fig. 4). First, Dimov and Damianova (1996) and later Gorinova et al. (2019) describe a major top-southeast ductile shear zone, which accommodates thrusting of the Kabul Unit onto the Malyovitsa Unit (Fig. 4). In both units, the foliation planes have a north-northwest to south-southeast-trending stretching lineation L 1 (Fig. 4). Shear sense criteria show a top-tothe-southeast sense-of-shear (Gorinova et al. 2019).
The highest confidence observations of extensional deformation in this area are made between Sapareva Banya and Panichishte (Fig. 4). Here, voluminous granitic melts were intruded into the metamorphic rocks, in a range from fine-grained biotite granites to leucocratic granites and pegmatites (Fig. 6f). Granitic dykes and veins are concordant with, or cross-cut the penetrative foliation S 1 that dips to north-northwest (Fig. 4). A large number of cross-cutting dykes of mostly biotite granites are devoid of any apparent solid-state overprint, whereas some of the concordant dykes display a fine foliation marked by the alignment of micas and stretched quartz. Northeast-dipping extensional C' shear bands that cut a shallow dipping, S1 foliation and granitic dykes, indicate a top-to-the-north sense-of-shear (Fig. 6f). A generally north-dipping stretching lineation, formed by stretched quartz and white mica flakes, is observed on the shear planes.
In the area west of Kalin pluton, Shipkova and Ivanov (2000) describe a brittle detachment fault (Dzherman Detachment) underlain by a ductile mylonitic zone along the northwest slope of the Rila Mountains. Later, Gorinova et al. (2019) confirm the existence of greenschist-facies metamorphic extensional mylonites in this area, where the brittle deformation superimposed on the mylonites is attributed to the movements along the younger Saparevo normal fault (Fig. 2). However, new field observations made for this study do not identify a brittle detachment in the area. Instead, the old contraction-related foliation is cut by steeper shear bands and centimetre-thick, low-grade shear zones indicating a top-to-the-north sense-of-shear, similar to the fabric observed in the Sapareva Banya-Panichishte area. A noteworthy difference in the area west of the Кalin pluton is the lack of any voluminous intrusions of granitic and pegmatitic dykes and veins.

Thermochronology
The 40 Ar/ 39 Ar and fission-track analytical procedures and methodologies are described in Online Appendix A. The geographic location of the analysed samples and the results are presented in Fig. 4 and Tables 1 and 2, whilst relevant 40 Ar/ 39 Ar age spectra are shown in Fig. 8. Use the term plateau age is according to the recommendation of Schaen et al. (2021), such that (i) a plateau must have at least 5 consecutive steps that constitute ≥ 50% of the total 39 Ar released by step-heating, (ii) the plateau steps must not form a positive or negative slope, and (iii) the 40 Ar/ 36 Ar atmospheric value must be 298.56 ± 0.31 (Lee et al. 2006). The 40 Ar/ 39 Ar date  is only considered to be robust if the plateau age overlaps with the inverse isochron date.

Ar/ 39 Ar data
A migmatite sample (SM284) from Variscan high-grade metamorphic rocks is taken from a pristine roadside outcrop east of Sapareva Banya (Fig. 4). The foliation, defined by flattened quartz-feldspar aggregates, white mica and biotite, dips westward at shallow angle. An aliquot of muscovite grains between 300 and 250 µm in size were hand-picked for 40 Ar/ 39 Ar analysis. The initial heating step yields a relatively young date of 76.82 ± 6.30 Ma, although the subsequent 5 heating steps yield a plateau date of 100.79 ± 0.55 Ma, and an inverse isochron date of 100.31 ± 2.98 Ma (MSWD = 2.39) (Fig. 8a). A biotite schist sample (SM286) from the Malyovitsa Unit is taken from a roadside outcrop southeast of Sapareva Banya (Fig. 4). Foliation in the schists dips to the north-northwest. A biotite concentrate was hand-picked from the fraction coarser than 160 µm. The initial (lower temperature) heating steps yield dates that span between 27.39 ± 1.29 and 33.89 ± 0.38 Ma, whilst seven subsequent heating steps, corresponding to ~ 80% of the 39 Ar released, define a plateau date of 34.90 ± 0.15 Ma, and an inverse isochron date of 34.961 ± 0.36 Ma (MSWD = 1.37) (Fig. 8b).

Fission-track data
The low number of apatite grains, combined with their low uranium content and high number of defects, did not allow measurement of sufficient horizontal confined track lengths to obtain useful track length statistics, inhibiting numerical thermal modelling (Rahn and Seward 2000). The migmatite sample (SM284) from the Variscan highgrade metamorphic rocks yielded a zircon FT Central age of 38.6 ± 1.9 Ma with P(χ 2 ) of 62 and an apatite FT Central age of 21.4 ± 1.5 Ma with P(χ 2 ) of 11 ( Fig. 4 and Table 2). Three measured confined tracks in apatite give a mean track length of 14.99 ± 0.95 μm, and a standard deviation of 1.65 μm (Table 2). Biotite schist sample (SM286) from the Malyovitsa Unit yielded a zircon FT Central age of 35.6 ± 5.6 Ma with P(χ 2 ) of 80 and an apatite Central age of 13.3 ± 1.1 Ma with P(χ 2 ) of 83 ( Fig. 4 and Table 2). Two measured confined tracks in apatite give a mean track length of 14.58 ± 0.55 μm, and a standard deviation of 0.78 μm (Table 2).

Interpretation and discussion
The low number of analysed samples and the lack of numerical thermal modelling of the apatite FT data limit constraints on the cooling histories of the rocks exposed in the study area. On the other hand, the combination of the 40 Ar/ 39 Ar and fission-track data from the same rocks in combination with previously published ages (Figs. 3 and 4) permits some conclusions to be made regarding the time and rates of exhumation. The detailed structural analysis and data in this study provide a link between the exhumation histories and structural evolution of the area. Finally, temporal correlations with tectonic phases identified in neighbouring areas strengthen the thermal history and tectonic conclusions.

Late Early to early Late Cretaceous evolution
Assuming the Ar isotope compositions of the muscovite are not modified by recrystallisation/fluid flow, the muscovite 40 Ar/ 39 Ar plateau date of the migmatised orthogneiss sample (SM284) from the Variscan high-grade metamorphic rocks suggest it reached a temperature between 440 and 400 °C (Harrison et al. 2009) at 101 ± 0.6 Ma (Fig. 4). This age corresponds to a late Early Cretaceous shortening event that resulted in north-westward thrusting of the Variscan high-grade metamorphic rocks onto the Kabul Unit along the Bistritsa Shear Zone (Gorinova et al. 2019, Fig. 2). If the sampled Variscan metamorphic rocks to the northeast of Dzherman valley is part of this thrust sheet, the burial depth of these rocks was at least between 10 and 14 km at the end of the Early Cretaceous (assuming a typical for orogenic belts geothermal gradient of 30 and 40 °C/km). Consequently, this overburden has been eroded during the mid-Cretaceous age shortening and/or a subsequent extension. Indeed, it is difficult to evaluate whether the muscovite plateau date of 101 ± 0.6 Ma of the migmatite (SM284) is related to cooling due to erosion synchronous with thrust sheet emplacement or post-contractional cooling. Certainly, these events are not related to the formation of a new ductile fabric in the studied part of the Variscan metamorphic rocks, because the only observed fabric is the migmatitic layering overprinted by high-grade foliation S VR (Fig. 4), which is cut by undeformed Carboniferous to Permian granites farther north (e.g. Carrigan et al. 2005).
North of the study area, Velichkova et al. (2004) report muscovite and biotite 40 Ar/ 39 Ar ages between 105 and 99 Ma related to an Alpine metamorphic overprint associated with ductile deformation along discrete strike-slip shear zones in the Variscan high-grade metamorphic rocks of the Sredna Gora Zone (Fig. 1). On the other hand, thermochronological studies from the neighbouring Kraishte Zone (Kounov et al. 2010) and farther west, in central parts of the Serbo-Macedonian Massif (Antić et al. 2016, Fig. 1), suggest the Albian to Cenomanian (112-90 Ma) is a period of regional extension. In Sredna Gora Zone to the north (Fig. 1), calcalkaline magmatic activity is coeval with this extension and the formation of intra-arc volcano-sedimentary basins, or possibly one large basin with several depocenters filled with Cenomanian to Maastrichtian deposits (e.g. Dabovski et al. 2009;Ivanov 2017). Therefore, it is possible the northwest Rhodope was also under post-contractional extension at that time, like some other parts of the Balkanides.

Time and rate of cooling and exhumation
The zircon FT age of 39 ± 2 Ma from the migmatite sample (SM284) in the hanging wall of the North Rhodopean Extensional System (Ograzhden Unit, Fig. 4) corresponds to the time of emplacement of the Shpanyovitsa granitic phase of the RWRB (U-Pb zircon age of 39 ± 0.2 Ma; Fig. 2) (Peytcheva et al. 2007) and, therefore, could be related to an increase in the thermal gradient in the area. On the other hand, the biotite schist sample (SM286) from the footwall of the extensional system (Malyovitsa Unit, Fig. 4) Schaen et al. (2021) 1 3 from each other. However, no further interpretation of the FT zircon age is attempted because the date is obtained from track density measurements only based on three grains (Table 2) and is thus not considered a statistically robust result.
Elsewhere around the study area, Gunnell et al. (2017) report an apatite FT age of 35 ± 1.8 Ma for a diorite located close to the Rila Lakes hut, at an elevation c. 1000 m above the biotite schist sample (SM286) (Fig. 4). Considering the altitude difference and closure temperatures of both mineral systems, a geothermal gradient of ~ 200 °C/km at ~ 35 Ma is deduced (Ar-in-biotite and FT-apatite; mean effective closure temperature of ~ 310 °C for Ar-biotite, Harrison et al. 1985; 110 ± 10 °C for apatite FT, Gleadow and Duddy 1981). Such a high geothermal gradient seems unrealistic and would be regionally anomalous, because there are no reports of magmatic intrusions of similar age in the area having created such a temperature difference. Consequently, a more likely interpretation is a high exhumation rate of ~ 5 km/Ma during the cooling of the whole ~ 1 km thickness of rocks between the two samples (using a geothermal gradient of 40 °C/km suggested by Gunnell et al. 2017 for eastern Rila during the Paleogene). An even higher estimate of exhumation rate in the area is derived from an 40 Ar/ 39 Ar phlogopite plateau age of 34.2 ± 0.4 Ma obtained from a granitic pegmatite vein located ~ 1 km southeast of Damga Peak (Alexandrov et al. 2001) (Fig. 4). The dated sample is taken from a similar elevation in the vicinity of the diorite sample Gunnell et al. (2017) report (Fig. 4). Considering the overlapping ages of both samples and the closure temperatures of Ar in white mica (~ 440-400 °C, Harrison et al. 2009) and the FT-apatite mineral systems (cited above), a cooling rate of ~ 330-290 °C/Ma are estimated. Again, using the assumption of a geothermal gradient of 40 °C/km (Gunnell et al. 2017), this corresponds to more than 7 km/Ma of exhumation at ~ 35-34 Ma. It is important to note that such high cooling rates are incompatible with simple erosiondriven, subaerial denudation of the Rila Mountains during the Cenozoic, as suggested by Gunnell et al. (2017). In this case, tectonically driven denudation in this area seems a more probable explanation. The possibility that the phlogopite 40 Ar/ 39 Ar age corresponds to the time a pegmatite vein forms, as suggested by Alexandrov et al. (2001), or consequent fluid flow may also be excluded. Late, syn-kinematic granitic pegmatitic and aplitic veins that intrude the Malyovitsa Unit have ages between 47 and 51 Ma (U-Pb zircon ages, Gorinova et al. 2019), whereas the youngest granitic pulse in the RWRB occurs at 39.39 ± 0.21 Ma (U-Pb zircon age, Peytcheva et al. 2007). Therefore, we attribute the very fast exhumation of these rocks at about 35-34 Ma to fault activity along the North Rhodopean Extensional System (Fig. 9a). Fast cooling between 38 and 35 Ma are also recorded in the eastern parts of the Rila Mountains, based on the results of thermal modelling of apatite FT data (Gunnell et al. 2017). A muscovite 40 Ar/ 39 Ar plateau age of 34.8 ± 0.2 Ma was also reported from the deformed Eocene granite in the eastern Rila Mountains near the Stankovi Baraki Dam's wall (Gerdjikov et al. 2006, Fig . 3).
Recently, some authors advocate that fast cooling episodes constrained by thermochronological data could be related to thermal relaxation, postdating a major orogenic event rather than to substantial direct exhumation (e.g. Braun 2016;Wolff et al. 2020). However, this model is discarded as a possibility that explains the exhumation in the study area because the time of the cooling stage estimated by the new data is coeval with both massive erosion and the deposition of a substantial amount of late Eocene age (Priabonian, 37.7-33.9 Ma, Cohen et al. 2022) sediments (Fig. 1, e.g. Zagorchev et al. 1989, see also the discussion below). Significant exhumation is also corroborated by the fact the 39-Ma-old Shpanyovitsa granite is exposed by early Oligocene (33.9-27.8 Ma, Cohen et al. 2022) and reburied under sediments of the Lakatitsa and Kostenets basins (Naydenov et al. 2008, Figs. 2 and 9b). The sedimentary rocks in the Padala Basin are also of early Oligocene age (Figs. 2 and 9b, Cernjavska 2000). In addition, the substantial amount of exhumation is evidenced by the fact that the upper greenschist-to lower amphibolite-facies metamorphic (mediumgrade) mylonites were overprinted by greenschist-facies metamorphic (low-grade) ultramylonites in the footwall of the North Rhodopean Detachment (Fig. 7d). This retrograde overprint and strain localisation in these rocks took place during ascent from deeper to progressively shallower crustal levels as exhumation progressed.

Differential exhumation along the NRES
Structural observations along the NRES suggest the amount of exhumation varies along the zone. The footwall in the eastern part of the system, where the structurally deeper part of the Rhodope Metamorphic Complex is exposed, the Chepino Unit (e.g. Burg 2012), experienced the maximum exhumation (Fig. 2). Logically, in this part of the system (Sestrimo area, Fig. 3), the thickest mylonitic zone is developed below the brittle detachment, where early medium-grade mylonites are later overprinted by low-grade mylonites and ultramylonites. Unlike the eastern domain, in the western domain (west of the Dzherman River valley, Fig. 2), there are no exposures of extensional medium-grade mylonites and deformation along the NRES is accommodated by low-grade mylonites and ductile-brittle to brittle deformation. If medium-grade mylonites did form at depth, they have never been exhumed to the surface. Accordingly, in this part of the NRES, only the uppermost parts of the  Fig. 2). In the area of the Padala basin, the relationship between the North Rhodopean Detachment and the ductile mylonitic shear zone (Shipkova and Ivanov 2000, Fig. 2) is still uncertain. The thermochronological methods in use in this study do not record the differential exhumation along the NRES because the footwall rocks are exhumed from depths with temperatures in excess of ~ 300-400 °C, which is higher than the closure temperature of the applied thermochronological methods.

Syn-extensional sedimentation
In the study area, Eocene age extension was either not accompanied by deposition of syn-tectonic sediments or the sedimentary rocks of this age are not preserved due to subsequent erosion. However, sediments derived from the substantial extension and denudation in the north Rila could possibly be deposited farther east in the Upper Thrace basin (Fig. 1), where Bartonian to Priabonian sedimentary rocks are penetrated in boreholes (e.g. Popov et al. 2015). Upper Eocene-lower Oligocene sediments are deposited in the hanging wall of the North Rhodopean Detachment farther east in the Izvor-Peshtera Basin (Fig. 1, Sarov et al. 2008a, b). Later, the detachment is sealed by the Oligocene volcanic rocks, the Bratsigovo-Dospat ignimbrites (U-Pb age of 31.63 ± 0.40 Ma, Marchev et al. 2022), which erupted into this basin. Late Eocene (Priabonian) sedimentary rocks associated with extension also outcrop in many other parts of the Rhodope Metamorphic Complex, south and southeast of the study area (Fig. 1, e.g. Cernjavska 1977;Boyanov and Goranov 2001), as well as in Kraishte Zone, adjacent to the northwest Rila Mountains (Fig. 1, Zagorchev et al. 1989). This widespread sedimentation is cited by many authors as evidence for the time of onset of regional extensional tectonics (e.g. Pleuger et al. 2011;Kounov et al. 2004Kounov et al. , 2020. The sedimentation in most areas continues into the early Oligocene, which suggests the cessation of activity along the North Rhodopean Detachment (Fig. 9b) is not necessarily coincident with the end of crustal extension in the area.

Correlations with neighbouring areas
The end of the Eocene and the beginning of the Oligocene is a time of substantial crustal extension, affecting the whole southern part of the Balkan Peninsula (e.g. Burchfiel et al. 2008;Kounov et al. 2018), which is characterised by core complexes formation and low-angle normal faulting across all parts of the Rhodope Metamorphic Complex (Dinter and Royden 1993;Krohe and Mposkos 2002;Bonev et al. 2006;Brun and Sokoutis 2007;Kounov et al. 2015Kounov et al. , 2020Georgiev et al. 2010Georgiev et al. , 2016, and the Kraishte area (e.g. Kounov et al. 2004;Antić et al. 2016). For example, the Central Rhodope Metamorphic Complex experiences rapid cooling between 37 and 33 Ma, from temperatures > 650 °C down to ~ 60 °C, associated with partial melting and formation of multiple generations of detachment faults (e.g. Kounov et al. 2020). The Eocene extensional event is also identified in the neighbouring Balkan Fold-Thrust Belt (Fig. 1), where it started in the middle Eocene (~ 44 Ma, Kounov et al. 2018). Recently, Balkanska et al. (2022) reported two stages of post-orogenic middle Eocene to Oligocene extension (45-42 and 34-28 Ma) from the central parts of the Sredna Gora Zone (Fig. 1) crystalline basement, which is ~ 50 km northeast of the studied area.

Late Oligocene-early Miocene heating event
Consideration of Miocene apatite FT ages from both sides of the North Rhodopean Detachment (13.3 ± 1.1 and 21.4 ± 1.5 Ma, Fig. 4) suggests reheating reset the FT ages of these rocks after being at, or close to, the surface in the early Oligocene, as suggested by their proximity to the transgressive contact with sedimentary rocks of the same age (Naydenov et al. 2008, Fig. 4). Based on their apatite (U-Th)/He and FT data, Danišík et al. (2014) define a late Oligocene-early Miocene heating event related to top-tothe-south thrusting of the North Rhodopean Unit over the metamorphic rocks and granites cropping out today in the Rila Mountains. The present-day structure attributed to this event is the North Rhodopean Thrust, which extends along the northern boundary of the Rila-Rhodope Massif, and along which the rocks of the North Rhodopean Unit and Late Cretaceous granites are thrust over lower Oligocene sedimentary rocks (e.g. Naydenov et al. 2008, Fig. 2). However, the tectonic settings and amount of displacement on this thrust are still disputed. Zagorchev (1992) suggests ~ 1 to 5 km displacement along this structure, whereas Naydenov et al. (2008) suggest it is a transpressional domain of the dextral strike-slip Maritsa Fault Zone, which formed during the latest brittle stage of the fault zone tectonic evolution ( Figs. 1 and 9c). Localised occurrences of late Oligocene to early Miocene contractional structures related to a shortlived period of regional shortening are reported from Thrace Basin in Greece (e.g. Kilias et al. 2013), as well as several other places across the Balkan Peninsula (e.g. Dumurdzanov et al. 2005;Burchfiel et al. 2008). This tectonic event is associated with an early Miocene depositional hiatus observed in most basins of the southern Balkan Peninsula (e.g. Zagorchev 1996;Dumurdzanov et al. 2005

Middle-late Miocene extension
The apatite FT ages of the samples from the Malyovitsa Unit (13.3 ± 1.1 Ma; SM286) and the Variscan high-grade metamorphic rocks (21.4 ± 1.5, Ma; SM284) suggest cooling of the analysed rocks in the study area to below 110 ± 10 °C must have taken place in the Miocene (Fig. 4). The brittle-ductile low-angle North Rhodopean Detachment was no longer active at this time as it was already sutured by lower Oligocene sedimentary rocks. However, some differential movement along the contact between these two units must still have occurred, given their different apatite FT ages. Today, normal faults border most of the northern slope of the Rila Mountains and the Saparevo and Sestrimo faults are the largest amongst them (Fig. 2). Displacement along these faults is related to the development of the Miocene to Quaternary Dzherman and Kostenets-Sestrimo basins (Sarov et al. 2011a;Naydenov et al. 2008;Figs. 2, 9d). Some displacement during this time may also occur along the Govedartsi normal fault, which is probably connected to the Saparevo normal fault along the Dzherman fault, which acts as a relay ramp breaching fault (e.g. Fossen and Rotevatn 2016, Fig. 9d). The apatite FT age from the uplifted footwall block of this fault system, suggests the onset of movement along the faults, and the subsequent growth of the basins probably started close to 13 Ma, or after. The diorite sample taken close to the Rila Lakes hut by Gunnell et al. (2017) yields a (U-Th)/He apatite age of 15 ± 1 Ma (Fig. 4), which further suggests some cooling and exhumation occurred during 15-13 Ma. An early to middle Miocene (22-13 Ma) tectonically induced pulse of denudation is suggested by Gunnell et al. (2017) based on apatite (U-Th)/He ages from the Rila Mountains. Middle-late Miocene cooling of the Rila basement complex was also reported by Danišík et al. (2014), based on combined apatite (U-Th)/He and FT data. The middle to late Miocene (~ 12-6 Ma) is also characterised by exhumation along brittle normal faults cutting the Strymon Valley Detachment and the formation of syntectonic sedimentary basins west of the Pirin Mountains, and farther south in the Southern Rhodope Metamorphic Complex (Kounov et al. 2015;Stübner et al. 2016; Fig. 1). These steeply dipping brittle normal faults led to rift-flank uplift and the formation of high topography subject to rapid erosion (Stübner et al. 2016). A similar scenario could be suggested also for our study area where possible high relief formed during the Miocene has partially persisted until today.

Multistage extension evolutionary model
The long and complicated evolution of the syn-to post-orogenic extension and exhumation along the orogenic belts can only be thoroughly revealed by a combination of appropriate thermochronological methods and detailed structural analysis. The coupling of new and previously published 40 Ar/ 39 Ar, fission-track and (U-Th)/He data, together with structural studies of the ductile to brittle deformation of rocks from the northwest part of the Rhodope Metamorphic Complex, allow reconstruction of the exhumation history of these rocks from 600 to less than 110 °C during the Cenozoic extension. The first phase of exhumation in the RMC during the Eocene (39-35 Ma), accommodated by the North Rhodopean Detachment, is related to back-arc extension triggered by slab rollback and Aegean trench retreat (e.g. Brun and Faccenna 2008;Forster and Lister 2009;Kounov et al. 2020). This forms classical extensional detachment zone along which amphibolite-faces metamorphic rocks have been continuously exhumed through greenschist-facies to finally brittle deformation conditions (e.g. Lister and Davis 1989). A decreasing amount of exhumation along the more than 150 km east-west trending detachment zone (from Krichim to Dupnitsa areas, Fig. 1) is observed from west to the east. It may tentatively be suggested that farther west the zone is connected with the Middle Rhodopean Detachment (Fig. 1), which has exhumed even deeper crustal rocks that have experienced anatexis during the extension (e.g. Kounov et al. 2020). After a phase of mostly strike-slip and thrust tectonics, the second middle-late Miocene extensional stage  is associated with the formation of high-angle brittle normal faults. It is difficult to say if these faults are listric and connect with a sole detachment at depth. If yes, such a structure is not exhumed, although the extension continues even during the Quaternary, indicated by the continuous sedimentation in the extensional basins surrounding the Rila Mountains (Fig. 2, Zagorchev 1992). Unlike the first stage, related to the Aegean back-arc extension, the second might just result from a westward extrusion of Anatolian plate along the North Anatolian Fault (NAF, Fig. 1), which pulls the southern Balkan lithosphere slowly to the south, causing regional north-south extension (e.g. Burchfiel et al. 2008). A lack of trench rollback related hot mantle material underlying the Rhodope Metamorphic Complex did not allow melt to spread and trigger stretching in the lower crust which was characteristic of the Eocene extension. Consequently, the middle to late Miocene extension is accommodated only by movement along high-angle brittle normal faults whereas middle to lower crustal rocks are not exhumed. Crustal extension continues even today in the area of Rila Mountains, as suggested by present-day seismicity data indicating active displacement is ongoing along some faults that bound the mountain chain (e.g. Kotzev et al. 2006). The highest present-day topography is in the Rila and Pirin Mountains in comparison to the rest of the Rhodope massif, which could be related to some post-glacial rebound