Glacial Geomorphology and Preliminary Glacier Reconstruction in the Jablanica Mountain , Macedonia , Central Balkan Peninsula

Although glacial landforms on the Balkan Peninsula have been studied since the 19th century, only scarce data are available about the extent of the former glaciations in the Central Balkan Peninsula, the transition zone between the Mediterranean and Central Europe. Glacial features of the Jablanica Mt. were mapped, described and classified into morphostratigraphic units. A revised glacio-geomorphological map was produced and glacial landforms were assigned to six morphostratigraphic units. Ten primary and two secondary cirques were identified in the upper parts of the studied valleys, while downstream the valleys were steep and glacially shaped with several glacial steps and thresholds. Cirque and valley morphology indicate that subglacial deepening was limited within the cirques and was more intensive in the valley sections during more extensive glacial phases. The largest reconstructed glaciers were 4.6–7 km long, while the last cirque glaciers were only a few hundred meters long. Using morphostratigraphic data, a glacier reconstruction was carried out for the largest mapped glacial extent. On the basis of glacial geomorphology, a former equilibrium-line altitude (ELA) of ~1800 m and glacier cover of 22.6 km2 were estimated during this stage. The local ELA values were compared to the regional ELA record and enabled to tentatively attribute a MIS 6 age for the reconstructed maximum ice extent in the study area.


Introduction
Glacier reconstruction enables the calculation of palaeoglaciological and palaeoclimatic variables [1].Such data are important for our understanding of past atmosphere-cryosphere interactions [2] and can serve as constraining boundaries for palaeoclimate models.However, in order to reach these goals, reconstruction of palaeoglaciers on the basis of the glacial landforms is essential [3].
Glacial landforms of the Balkan Peninsula have been studied since the late 19th century [4][5][6][7].Several lines of evidence suggest extensive past glaciations of the coastal ranges of the Adriatic Sea in the form of ice caps and valley glaciers [8].However, scarce data are available about the extent of the glaciations in the Central Balkan Peninsula.
The description of the glacier extent during past glaciations was completed by age determination of glacial landforms in Montenegro and northern Greece using U-series techniques [9][10][11][12][13], with the results suggesting repeated glaciations during at least three separate cold stages, the oldest of which occurred earlier than 350 ka.These glaciations were correlated with the MIS 12, MIS 6 and MIS 5d-2 cold stages, with the most extended glaciation during the oldest phase.At other locations surface exposure dating using cosmogenic radionuclides (SED) of moraines was used to estimate the age of glacial stages.In the Šar and Pelister Mountains of Macedonia and in the Rila Mts.(Bulgaria) preliminary 10 Be SED ages were obtained [14][15][16]; and 36 Cl SED dating was recently applied in the Greece [17,18] and southern Macedonia (Galičica Mountain) [19].Other than Mt.Chelmos in Southern Greece [18], no evidence of a glaciation older than the Last Glacial Maximum (LGM, MIS 2) was described by these studies.Extended ice caps were also reported in the coastal ranges of northern and central Croatia, with glaciers extending down to or possibly the present-day sea level [20,21], however no numerical age control is available.
During the Holocene most of the mountains on the Balkan Peninsula were unglaciated.Today only some small (and rapidly shrinking) glaciers or permanent snow patches above 2000 m elevation a.s.l.exist in the region [22].Their mass balance heavily depends on topoclimatic factors, such as receiving extra snow by avalanching and wind drift, and being shaded by steep cirque-walls [23].
The reconstructed equilibrium-line altitudes (ELAs) for the most extended glacial phase in the coastal range are under 1300 m (the lowest in the Mediterranean), and its elevation manifests a rapid increase inland with values around 1800-2200 m in the Central Balkan Peninsula [2,8,14].However, this pattern is ambiguous because glaciers reached their maximum extent asynchronously between MIS 12 and MIS 2. The apparently extensive ice caps in the coastal ranges along the Adriatic Sea and more limited extent of inland glaciations suggest that the coastal ranges with large icefields blocked the inland penetration of humid air masses from south to north, resulting in much drier conditions in the Balkan interior and in the Pannonian Basin compared to the coastal margin of the Balkan Peninsula [10].
Despite the large number of studies devoted to the glaciations of the mountains of the Balkan Peninsula, information on the glacial history of the transitional areas between the coastal ranges and the inland areas is limited.There are many unresolved or contradictory issues regarding the glacial chronology of the Central Balkan Peninsula due to the scarcity of chronological data.This study focused on glacier reconstruction-based on detailed glacial geomorphology of the area-and assessment of the former ELAs of the Jablanica Mt.The mapped glacial landforms will provide the necessary stepping stone towards a future geochronological study [24].

Location and Geological Setting of the Jablanica Mountain
The Jablanica Mountain is a 45 km long, 25 km wide and NNW-SSE oriented mountain range, part of the Dinaride-Hellenide mountain belt, located in the central part of the Balkan Peninsula, along the border between Albania and Macedonia [25] (Figure 1A).
Morphologically, the range is divided in three parts (Figure 1B).The southern part (between the Mokra Mt. (Albania) to the south and the Belica Valley), is mostly plateau-like with karstic landscapes culminating at Čumin Vrv (2125 m).
The middle part (between the Belica and Piskupština Valleys) is the highest part of the mountain (highest peaks: Crn Kamen, 2257; Strižak, 2233) with a glacially sculptured plateau-like crest [25].This study focused on this part of the mountain, comprising the valleys of Vevčani, Podgorci and Labuništa-Pupuljak (Figure 1B).These valleys strike from west to east draining to the Black Drin (Crn Drim) river.
The northern part (from the Piskupština Valley to the Debar Basin) has generally lower elevation (highest peaks: Raduč, 2083; Strmnica, 1904) (Figure 1B).The folded structure of the Jablanica Mt. (mostly N-S-to NW-SE-oriented larger synclines and anticlines) has been cut by several normal faults during the Neogene extension affecting this part of the Balkan Peninsula [26].On the east it is bounded by Neogene sedimentary basins and by the Shebenik Mt. on the west.[25].The basins were initiated during the Late Miocene along a NNWoriented fault line, as part of the South Balkan Extensional System [26,27], and were filled with lacustrine deposits of Upper Miocene, Pliocene and Quaternary age.The Ohrid graben is still actively developing, hosting the 285 m deep Ohrid Lake [28].The Jablanica Mt. has uplifted at the footwall of the normal faults bounding the subsiding grabens to the east, with most likely faster uplift during the Pleistocene [29].
The Jablanica Mt. is mostly composed of Triassic sediments overlying Paleozoic (Devonian) metamorphic rocks and is partly covered by Jurassic sediments.These rocks were intruded by Jurassic mafic-ultramafic magmatic rocks.The eastern and northeastern foothills are covered by Neogene lacustrine sediments [30] (Figure 2).The folded structure of the Jablanica Mt. (mostly N-S-to NW-SE-oriented larger synclines and anticlines) has been cut by several normal faults during the Neogene extension affecting this part of the Balkan Peninsula [26].On the east it is bounded by Neogene sedimentary basins and by the Shebenik Mt. on the west.[25].The basins were initiated during the Late Miocene along a NNW-oriented fault line, as part of the South Balkan Extensional System [26,27], and were filled with lacustrine deposits of Upper Miocene, Pliocene and Quaternary age.The Ohrid graben is still actively developing, hosting the 285 m deep Ohrid Lake [28].The Jablanica Mt. has uplifted at the footwall of the normal faults bounding the subsiding grabens to the east, with most likely faster uplift during the Pleistocene [29].
The Jablanica Mt. is mostly composed of Triassic sediments overlying Paleozoic (Devonian) metamorphic rocks and is partly covered by Jurassic sediments.These rocks were intruded by Jurassic mafic-ultramafic magmatic rocks.The eastern and northeastern foothills are covered by Neogene lacustrine sediments [30] (Figure 2).
The study area, the middle part of the Jablanica Mt., is almost entirely composed of Triassic rocks, with the underlying Paleozoic schist only partly exposed at lower elevations.The stratigraphically lower part of the Triassic section is built of clastic sediments (sandstone and shale), covered by a thick carbonate section, composed of massive to platy and cherty limestone and dolomite.The mountain crest in the study area is almost entirely built of Jurassic conglomerate.This formation is composed of fragments of the Jurassic magmatic and the Triassic clastic and carbonate formations and it is cut by a number of quartz veins [30].The study area, the middle part of the Jablanica Mt., is almost entirely composed of Triassic rocks, with the underlying Paleozoic schist only partly exposed at lower elevations.The stratigraphically lower part of the Triassic section is built of clastic sediments (sandstone and shale), covered by a thick carbonate section, composed of massive to platy and cherty limestone and dolomite.The mountain crest in the study area is almost entirely built of Jurassic conglomerate.This formation is composed of fragments of the Jurassic magmatic and the Triassic clastic and carbonate formations and it is cut by a number of quartz veins [30].

Introduction to the Glacial Geomorphology of the Study Area and Previous Results
Currently, neither glaciers nor permanent firn patches can be found on the Jablanica Mt., however geomorphological evidence of past glaciation is abundant [25,31,32].
The first evidence of former glaciation on the Jablanica Mt. was presented by Dedijer [6].He visited the area north of the Belica Valley, and described the largest cirques with a system of glacial valleys and the development of stepped cirques.He noted the lack of cirques in the Belica Valley (Figure 1B), although it has a clear glacial morphology.He noted a number of cirques in the upper parts of Labuništa, Podgorci and Vevčani valleys, as well as the extensive moraine deposits with large boulders between the Labuništa and Podgorci rivers.In the northern area, a less-developed glacial morphology was described.He reported the lowest moraines at an elevation of 750 m a.s.l. in the north (near Modrič village; Figure 1B), and near or below the villages of Vevčani, Podgorci and Labuništa (850-900 m) in the central part of the range, and suggested two glaciations in the area.In addition, Cvijić [5] described a possible terminal moraine even lower, at 700 m elevation, to the north of Labuništa village, at the confluence of Borovec River and the Black Drin (Figure 1B).
Kolčakovski [25] focused on the glacial and periglacial processes and forms of the Jablanica Mt.He worked generally in the upper parts (above 1500 m).He considered the sediments below ~1700 m as fluvioglacial and, thus did not discuss the glacial deposits previously reported at lower elevations.According to his observations the largest glacier that occupied the Labuništa-Pupuljak valley system descended towards a terminal moraine at 1732 m (Branica locality; Figure 2).Therefore,

Introduction to the Glacial Geomorphology of the Study Area and Previous Results
Currently, neither glaciers nor permanent firn patches can be found on the Jablanica Mt., however geomorphological evidence of past glaciation is abundant [25,31,32].
The first evidence of former glaciation on the Jablanica Mt. was presented by Dedijer [6].He visited the area north of the Belica Valley, and described the largest cirques with a system of glacial valleys and the development of stepped cirques.He noted the lack of cirques in the Belica Valley (Figure 1B), although it has a clear glacial morphology.He noted a number of cirques in the upper parts of Labuništa, Podgorci and Vevčani valleys, as well as the extensive moraine deposits with large boulders between the Labuništa and Podgorci rivers.In the northern area, a less-developed glacial morphology was described.He reported the lowest moraines at an elevation of 750 m a.s.l. in the north (near Modrič village; Figure 1B), and near or below the villages of Vevčani, Podgorci and Labuništa (850-900 m) in the central part of the range, and suggested two glaciations in the area.In addition, Cvijić [5] described a possible terminal moraine even lower, at 700 m elevation, to the north of Labuništa village, at the confluence of Borovec River and the Black Drin (Figure 1B).
Kolčakovski [25] focused on the glacial and periglacial processes and forms of the Jablanica Mt.He worked generally in the upper parts (above 1500 m).He considered the sediments below ~1700 m as fluvioglacial and, thus did not discuss the glacial deposits previously reported at lower elevations.According to his observations the largest glacier that occupied the Labuništa-Pupuljak valley system descended towards a terminal moraine at 1732 m (Branica locality; Figure 2).Therefore, the age of this glacial extent was assigned to the LGM with an estimated ELA of 1964 m.For Podgorci and Vevčani valleys he suggested that the most extensive glaciers reached an elevation of 1600-1700 m a.s.l.

The Climate of the Jablanica Mountain
The climate of the Jablanica Mt. is strongly influenced by the nearby Ohrid Lake and the Adriatic Sea, which is located only 80 km to the west (Figure 1A).The Ohrid Lake has more influence on the temperature regime than the Adriatic Sea, with the latter having higher effect on the precipitation of the study area.According to the WorldClim v2 model [33] and MODIS data [34], the highest parts of the mountain (above 1800 m) have mean annual temperatures between 3-4 • C, while toward the foothills, mean annual temperatures increase to 10-11 • C. The mean summer temperatures (July, August, September) are between 12 • C at the highest elevations, to 22 • C at the mountain foot.The mean January temperatures fall in the range of −5 • C at the highest elevations, to 1.5 • C in the foothills.Based on this temperature range, the local lapse rate is calculated to 0.58 • C/100 m.
According to the WorldClim v2 model, the mean annual precipitation increases from 800 mm at the mountain foot to 1150 mm at the highest parts (above 1800 m).The precipitation regime is Mediterranean, with the greatest amount of precipitation falling during the winter, and lowest during the summer.The mean monthly precipitation in January is about 100-110 mm in the highest part to 80 mm in the bottom, whereas mean August precipitation is only half of that: 50-60 mm at the top, to 40-45 mm at the foothill of the mountain.Winter precipitation is mostly snow, especially in the higher parts of the mountain.Thus, according to the MODIS data, average snow cover in the highest part is about 5 months, i.e., from November to April and the snow depth can be significant-up to 2 m.

Geomorphological Mapping
Geomorphological mapping was conducted on the basis of field investigations and combination of available data sources: topographic maps at 1:25,000 scale [35], a basic geological map in 1:100,000 scale [30], GoogleEarth images at varying scales, and a 9-m resolution digital elevation model (DEM).The DEM was generated from the 5-m (2009) and from the 20-m (2005) DEM model of the Agency for Real Estate and Cadastre of the Republic of Macedonia [36].
For cirque and cirque-floor delineation a combination of field identification, topographic maps, and DEM-derived contour and slope maps was utilized [37].In the field a number of ground control points were collected using a hand-held GPS as a reference for the mapped geomorphological features.Digital cartographic work, morphometric analyses and final map production were done using Global Mapper (v17, Blue Marble Geographics, Hallowell, ME, USA) and ArcGIS (v10.1,ESRI, Redlands, CA, USA) software.
Moraines were classified into morphostratigraphic units [38] based on their elevation, relative position and field observations on their preservation state (e.g., degree of erosion).
Because the study area comprises nearby valleys of similar aspect and elevation range, they are expected to have shared a similar ELA during the same glacial phase.Therefore, ELA estimates of the reconstructed glaciers (see Section 3.3) were used to verify the preliminary morphostratigraphic classification of the studied landforms.
Basic cirque morphometric parameters were derived from the 9-m DEM using the previously delineated cirque boundaries in ArcGIS v10.1 and Global Mapper v17.Cirque aspect is given as median axis aspect, cirque length (L) as the length of the median axis, cirque width (W) as the longest line perpendicular to the median axis, cirque height (H) as the vertical distance between the lowest and highest point within the cirque boundaries, cirque size as the cube root of volume ( 3 √ L × W × H), and cirque floor altitude as the modal floor altitude [37].We also measured cirque profile closure (taken as maximum slope of the headwall) and plan closure [39] and calculated the exponents for logarithmic (power) regressions of the size parameters (length, width, height) on overall cirque size to check for allometric development of the cirques [40].Based on their morphological and morphometric characteristics we attempted to classify the identified glacial cirques according to the classification system of Evans and Cox [41].

Glacier Reconstruction
Glacier reconstruction was based on the position of the landforms of glacial accumulation (lateral and terminal moraines) and erosion (U-shaped valleys, trimlines, glacially sculpted bedrock, glacial steps and over deepened valley sections).Furthermore, where the position of former moraines was inferred from a combination of topographic features (like change of valley slope, narrowing valley section with incised creek, scattered erratic boulders) allowing glacier terminations to be tentatively mapped in places where terminal moraines have been removed by erosion.The reconstructed glaciers were created using the GlaRe (Glacier Reconstruction) ArcGIS toolbox [42] using the 9-m DEM.The GlaRe toolbox offers a semi-automated method of glacier reconstruction.It generates the ice thickness from bed topography along a user defined flowline.It applies the physical laws for ice flow to generate a 3D palaeoglacier surface using several interpolation methods.As the modern land surface differs from the concurrent land surface below the glacier, the postglacial incision or sedimentation will make the reconstructed glacier thicker or thinner, respectively.
The ice thickness was adjusted to the mapped glacial landforms by increasing or decreasing the basal shear stress value (τb) as this is a first-order control on the output glacier surface.The ice thickness is calculated by the software along the flowline using a user defined basal shear stress value (Flowline ice thickness tool).Then the glacier surface was created by the Glacier surface interpolation tool, applying the "Topo to raster" setting.
The most accepted values of the basal shear stress vary within the ~50-150 kPa range with possible increase to 190 kPa [43][44][45].The τb was iteratively tuned within this realistic range until the virtual palaeoglacier surface best matched the mapped glacial features of the topography.Although basal shear stress may vary under a real glacier, we used a uniform value of τb for the preliminary modelling tests.The tests suggested no major change in the shape of the reconstructed glacier when this parameter was set to be variable along the glacier.
Where topography suggested glacier confluence, the present day topographic divides were set as ice divide.The elevation of the neighbouring glaciers could be adjusted by fine-tuning their thickness using the τb value.

Calculation of the Paleo-ELAs
The equilibrium-line altitude (ELA) is the spatially averaged elevation of the points on the surface of the glacier where the climatic mass balance is zero [46].It represents the lower boundary of annual snow accumulation.Accordingly, the ELA is characteristic of the prevailing climate conditions [47].
The palaeo-ELAs were estimated using the ELA Calculation toolbox [48].For the calculations the Accumulation Area Ratio (AAR) method was used.The main input needed for calculation of the AAR is the reconstructed 3D glacier-surface, the creation of which was previously described (Section 3.2).
To consider the glacier geometry and hypsometry, the AAR values were selected depending on the glacier size.The relationship between the glacier area and the balanced-budget AAR can be optimally described by a logarithmic regression equation [49]: where ssAAR 0 is the size-specific balanced-budget AAR, and S is the planimetric area of the reconstructed glacier (km 2 ) The palaeo-ELAs were also estimated by the Area-Altitude Balance Ratio (AABR) method.This method is frequently recognised as more robust (compared to the simple AAR) as it is based on the hypsometry of the former glacier and the ratio between the accumulation and ablation gradients [50,51].Balance ratios vary with regional climate (~1 in polar regions and ~2 in tropical regions).The regional AABR of 1.6, calculated from the contemporary glaciers of the European Alps [52], was assumed as potentially representative for the glaciated Jablanica Mt. and used in the current study.The palaeo-ELAs estimated by both methods enable regional comparisons with other studies using the AAR or the AABR method for ELA calculations.Besides, it gives the opportunity to compare the results provided by the two methods, to assess uncertainties of ELA estimates.

Results
The studied valleys have well developed glacial cirques at their heads along the mountain crest.Cirque floors of all valleys, except for the Pupuljak Valley, host glacial lakes, such as Vevčani Lake located at 1955 m, Podgorci Lake at 1894 m, and Upper and Lower Labuništa Lakes located at 1918 and 1913 m, respectively.Glacially sculptured landscapes are well preserved in the upper parts of the valleys with glacial sediments (i.e., till) found down to elevations of approx.1000 m (Figure 3).
regional AABR of 1.6, calculated from the contemporary glaciers of the European Alps [52], was assumed as potentially representative for the glaciated Jablanica Mt. and used in the current study.The palaeo-ELAs estimated by both methods enable regional comparisons with other studies using the AAR or the AABR method for ELA calculations.Besides, it gives the opportunity to compare the results provided by the two methods, to assess uncertainties of ELA estimates.

Results
The studied valleys have well developed glacial cirques at their heads along the mountain crest.Cirque floors of all valleys, except for the Pupuljak Valley, host glacial lakes, such as Vevčani Lake located at 1955 m, Podgorci Lake at 1894 m, and Upper and Lower Labuništa Lakes located at 1918 and 1913 m, respectively.Glacially sculptured landscapes are well preserved in the upper parts of the valleys with glacial sediments (i.e., till) found down to elevations of approx.1000 m (Figure 3).

Glacial Morphology and Deposits
The valley heads of the studied glacially sculptured valleys (Vevčani, Podgorci, and Labuništa valleys) form simple and compound glacial cirques along the mountain crest (Figures 3 and 5; see more in Section 4.2).Although affected by later slope erosion and post-glacial valley incision, these valleys exhibit a wide, U-shaped morphology narrowing downstream.It is an open question whether the glaciers extended to the foothills during their most extended phases, as it was suggested by Cvijić [5], or terminated as valley glaciers at or above the mountain front (Figure 3).Presently, no evidence of glaciers reaching 700-750 m elevation along the foothills could be found, but the possibility of former piedmont glaciers cannot be excluded.

Glacial Morphology and Deposits
The valley heads of the studied glacially sculptured valleys (Vevčani, Podgorci, and Labuništa valleys) form simple and compound glacial cirques along the mountain crest (Figures 3 and 5; see more in Section 4.2).Although affected by later slope erosion and post-glacial valley incision, these valleys exhibit a wide, U-shaped morphology narrowing downstream.It is an open question whether the glaciers extended to the foothills during their most extended phases, as it was suggested by Cvijić [5], or terminated as valley glaciers at or above the mountain front (Figure 3).Presently, no evidence of glaciers reaching 700-750 m elevation along the foothills could be found, but the possibility of former piedmont glaciers cannot be excluded.The main ridge reaches an elevation between 2000 and 2200 m, and likely preserves a remnant of a Paleogene planation surface, which can be found on other mountains in Macedonia [29].Several sections of the main ridge are smoothed, as well as multiple segments of the secondary ridges separating the cirques (Figure 4A).This morphology may reflect the former existence of a major The main ridge reaches an elevation between 2000 and 2200 m, and likely preserves a remnant of a Paleogene planation surface, which can be found on other mountains in Macedonia [29].Several sections of the main ridge are smoothed, as well as multiple segments of the secondary ridges separating the cirques (Figure 4A).This morphology may reflect the former existence of a major interconnected icefield with limited number of unglaciated peaks (e.g., Crn Kamen, Pupuljak, Strižak).
Most of the valley floors are covered by glacial deposits, with creeks generally incised along the border of the glacial deposits and former glacial valley sides in the lower valleys (Figures 2 and 3).In the upper valley sections morphologically well expressed glacial steps have developed frequently at the contact of lithological boundaries and/or along tectonic structures (Figures 2, 3 and 5C).Most of the valley floors are covered by glacial deposits, with creeks generally incised along the border of the glacial deposits and former glacial valley sides in the lower valleys (Figures 2 and 3).In the upper valley sections morphologically well expressed glacial steps have developed frequently at the contact of lithological boundaries and/or along tectonic structures (Figures 2, 3 and 5C).Glacial deposits (Jablanica Glacial Formation) are found throughout the study area, extending down to elevations of ~1100 m at Vevčani and Podgorci valleys (Figure 4B,C), and possibly to ~950 m at Labuništa Valley.Boulders of 0.5 to 10 m in diameter are the most prominent features of the glacial sediments (Figures 4 and 5).The moraine boulders are mostly (90%) composed of Jurassic conglomerate with smaller amount of Triassic sandstone, shale and limestone.The boulders composed of Jurassic conglomerate are best preserved, and because their source bedrock can only be found along the mountain crest (Figure 2), their presence downvalley was useful for the clear identification of glacial and fluvioglacial sediments.Based on the position of the glacial sediments the Jablanica Glacial Formation was classified in several members.These members are considered to belong to the same chronostratigraphic unit, based on their position, morphology and sediment properties [38].The landforms are best preserved in the Vevčani Valley, therefore most of the type localities and stratigraphic names come from this valley.
The studied glacial valleys are open towards their foothills, with no characteristic glacier terminations.In the Vevčani and Podgorci valleys the presence of large subrounded boulders of Jurassic conglomerate, the Vevčani Member, marks the extent of the largest mapped glaciers.Identification of moraines in the middle and lower part of all valleys was challenging due their poor preservation and dense forest vegetation (Figure 4B-D).However, ridges containing large boulders of Jurassic conglomerate and characteristic changes of valley morphology served as the basis for the identification of the Kutel and Leništa members.
Glacial landforms were better visible on the mountain pastures, although postglacial (e.g., periglacial) modification hindered the identification of primary glacial landforms.Nonetheless, the Golina Member is typically represented by lateral moraines below and occasionally also above the glacial steps, around ~1900-1700 m a.s.l.(Figures 3 and 5B,E).
The upper parts of all of the studied valleys are characterised by wide glacial cirques mostly with an upper and a lower moraine generation (Lokva and Lincura Members; Figure 5B-E).

Vevčani Member
Very large boulders of Jurassic conglomerate can be found in the lowest part of the Vevčani Valley, at the elevation of ~1100 m.The moraine ridge was crosscut by the creek, therefore a latero-terminal moraine was mapped based on the presence of large boulders (Figures 3 and 4B).
The lowest part of the Podgorci Valley, at the elevation of ~1200 m, has a large number of sizable Jurassic conglomerate boulders, but with no clear morphological expression of a moraine (Figure 4C).A blockfield with sub-rounded boulders may evidence fluvioglacial modification of a former moraine deposit, and finer (from cobbles to sand) fluvioglacial sediments were mapped further down (Figure 4E).The former glacier termination was therefore mapped at the down-valley end of the deposit of large sub-rounded boulders.
In the Labuništa Valley, again no clear expression of a terminal moraine could be found.The presumed termination of the glacier tongue was inferred from a slight narrowing and (~10 m high) bulging of the valley floor.In the lower part, at the elevation of ~1325 m, there is a remnant of a left lateral, or mid moraine between the main valley and its northern tributary Sušica Valley, composed of huge (~8 m high) boulders of Jurassic conglomerate, which can be attributed to this phase.
In the middle part of Podgorci and Labuništa valleys an extensive area of glacial deposits developed during this phase, with the current watershed divide located in the area covered by the glacial deposits.This area appears to represent confluence of glaciers in the middle parts of both valleys, which were separated by a large hill of Palaeozoic schists in their lower reach (Figures 2 and 3).

Kutel Member
In the Vevčani Valley two moraine crests with abundant large (up to 3 m high) boulders of Jurassic conglomerates were interpreted as lateral moraines (Figures 3 and 4D).In the Labuništa Valley this member is represented by the termination of a flat valley section followed downstream by a steeper and incised reach, interpreted as a moraine.In the Podgorci Valley glacial sediments are widespread in the valley but with no clear indication of a moraine.The position of the former glacier terminus was inferred by a slight (~10 m high) bulging of the valley-floor covered by glacial sediments.

Leništa Member
The Leništa Member is situated in the middle reach of each valley.In the Vevčani Valley, upstream from the steeper, forested segment, a lateral moraine with an uncertain termination was mapped on the northern side of a gently dipping, wide valley section covered by glacial sediments (Golina; Figure 3).
Similarly, a lateral moraine on the northern side of the Podgorci Valley could be followed below a mid-valley glacial step, down to ~1650 m, where the former glacier terminus is denoted by the termination of the lateral moraine.Glacial sediments re-shaped by periglacial rock-glacier activity above the valley scarp at 1840-1830 m a.s.l.elevation may belong to this phase as well.An inner and an outer lateral moraine ridge were also described to the northeast of the Podgorci Lake, which may belong to this phase.
In the Labuništa Valley, the position of the Leništa Member is represented by the presence of large boulders of Jurassic conglomerate where a break in the valley slope and incision of the creek in glacial sediment was observed (Leništa location; Figures 3 and 4G).This landform has been identified as a poorly preserved terminal moraine.
After the confluence of the Pupuljak and Labuništa valleys, Jurassic conglomerate boulders carried from the highest parts of Labuništa Valley were deposited as erratics on the limestone bedrock sill above the scarp (Branica locality, Figure 4F).Kolčakovski [25] described a terminal moraine at Branica.Most probably the widespread erratic boulders of Jurassic conglomerate occurring on the glacial threshold were mistaken for a moraine (Figure 4F).

Golina Member
On the above described gently dipping valley section (Golina) of the Vevčani Valley, downstream from the glacial step, a well preserved double crested latero-terminal moraine complex was mapped (Figures 3 and 5A).Along the southern valley side, the lateral moraines have been modified by a series of rock glaciers.
Cirque floors of the Podgorci and Labuništa-Pupuljak valleys are at lower elevations compared to the Vevčani Valley.Therefore, in these valleys glaciers had retreated to their cirques by this phase.
In the Podgorci Valley the sediments belonging to the Golina Member are present around Podgorci Lake.The end moraine to the east of the lake is not well expressed due to subsequent periglacial modification.
In the Labuništa-Pupuljak system, sediments of the Golina Member are present around Labuništa Lake, similarly to the Podgorci Valley.
The Pupuljak Valley, a northern tributary of the Labuništa Valley, is the only glacial valley that lacks glacial deposits, although it preserves a well-expressed U-shape morphology.It is almost entirely developed in highly fractured carbonate rocks, which is probably the reason for the poor preservation of glacial deposits.In the Sušica Valley the termination of this phase was outlined based on a faint morphological indication of a small (2 m high) arcuate ridge at the elevation of 1810-1900 m a.s.l.

Lokva Member
In the Vevčani Valley a rock glacier extends down to the northern side of Vevčani Lake.The lake is enclosed on the eastern side by a partly preserved double crested moraine at an elevation of ~1975 m, below which a glacial step can be seen.This moraine is attributed to the Lokva Member (Figures 3 and 5C).
In the Podgorci and Labuništa-Pupuljak valleys the presence of sediments belonging to the Lokva Member is uncertain.In the Podgorci Valley the rock glacier activity may have overwritten the previous glacier deposits (Figures 3 and 5B).Moreover, this section of the valleys is relatively steep, providing no space for the development of well-preserved moraines.Another possibility is that the shrinking glaciers were confined to the cirques in these valleys for longer time, while in the Vevčani Valley the glacier was forced back to the top of the scarp.The latero-terminal moraine to the north of the Lower Labuništa Lake may belong to this member.

Lincura Member
In the Vevčani Valley the highest moraines are found at the elevation of 2040 m, with the best expressed moraine to the east of Crn Kamen (2257 m) (Figures 3 and 5C).
In the Podgorci Valley the highest moraines are at the elevation of 2020-2070 m (Figures 3 and  5B).They are least preserved in the southern cirque, where glacial deposits have been reworked by (up to 0.4 km 2 ) large rock glaciers.
The upper part of the cirques in Labuništa Valley are developed in Jurassic conglomerate, with the highest moraine at the elevation of 2030 m, attributed to the Lincura Member.Most of the Pupuljak Valley has been carved in carbonate rocks and glacial deposits are missing, except in the cirque to the west of the Pupuljak peak, where the bedrock is Jurassic conglomerate, and a moraine was deposited at ~2035 m elevation (Figures 3 and 5D,E).
The elevation of the higher cirque in the Sušica Valley may suggest the survival of a small glacier until this phase, but without clear morphological or sedimentological indications of its termination.

Cirque Morphometry
Ten primary cirques with two secondary cirques were identified in the upper parts of the studied valleys, forming compound cirques (Figure 6).The average elevation of the main cirques ranges from 1929 m to 2078 m, with their total area between 0.3 km 2 and 0.9 km 2 .Cirque floor elevations are between 1838 m and 1995 m, and cirque floor areas between 0.03 km 2 and 0.36 km 2 .The cirques are from 218 m to 399 m high, and from 598 m to 1029 m long.Most of them have an eastern and north-eastern (median axis) aspect, but two have southern and south-eastern aspect (Figure 7A).The two secondary cirques are smaller in all sizes, with cirque floors at similar elevations (2035-2051 m), cirque floor areas of 0.02 km 2 and 0.09 km 2 , and are oriented to the north and east.
the previous glacier deposits (Figures 3 and 5B).Moreover, this section of the valleys is relatively steep, providing no space for the development of well-preserved moraines.Another possibility is that the shrinking glaciers were confined to the cirques in these valleys for longer time, while in the Vevčani Valley the glacier was forced back to the top of the scarp.The latero-terminal moraine to the north of the Lower Labuništa Lake may belong to this member.

Lincura Member
In the Vevčani Valley the highest moraines are found at the elevation of 2040 m, with the best expressed moraine to the east of Crn Kamen (2257 m) (Figures 3 and 5C).
In the Podgorci Valley the highest moraines are at the elevation of 2020-2070 m (Figures 3 and  5B).They are least preserved in the southern cirque, where glacial deposits have been reworked by (up to 0.4 km 2 ) large rock glaciers.
The upper part of the cirques in Labuništa Valley are developed in Jurassic conglomerate, with the highest moraine at the elevation of 2030 m, attributed to the Lincura Member.Most of the Pupuljak Valley has been carved in carbonate rocks and glacial deposits are missing, except in the cirque to the west of the Pupuljak peak, where the bedrock is Jurassic conglomerate, and a moraine was deposited at ~2035 m elevation (Figures 3 and 5D,E).
The elevation of the higher cirque in the Sušica Valley may suggest the survival of a small glacier until this phase, but without clear morphological or sedimentological indications of its termination.

Cirque Morphometry
Ten primary cirques with two secondary cirques were identified in the upper parts of the studied valleys, forming compound cirques (Figure 6).The average elevation of the main cirques ranges from 1929 m to 2078 m, with their total area between 0.3 km 2 and 0.9 km 2 .Cirque floor elevations are between 1838 m and 1995 m, and cirque floor areas between 0.03 km 2 and 0.36 km 2 .The cirques are from 218 m to 399 m high, and from 598 m to 1029 m long.Most of them have an eastern and northeastern (median axis) aspect, but two have southern and south-eastern aspect (Figure 7A).The two secondary cirques are smaller in all sizes, with cirque floors at similar elevations (2035-2051 m), cirque floor areas of 0.02 km 2 and 0.09 km 2 , and are oriented to the north and east.Table 1.General morphometric characteristics of the glacial cirques.Grade system according to [41], with grade 3 representing a "definite" cirque and grade 4 a "poor" cirque.The morphometric values of the main cirques on Jablanica Mt. (Table 1) are comparable to other regions of the world [37], and the main dimensional parameters (area, L/W, L/H, W/H) are very similar to the cirques from other mountains in the region (e.g., Northern Greece [53,54]; Retezat Mts., Romania [55]).Cirque size, as expected, generally decreases with increase of cirque floor elevation (Figure 7B).
All of the cirques have relatively high plan closure values (>174 • ), and based on both plan and profile closure, three groups can be delineated: the first group (cirques 1, 2, 8, 8a, 9 and 10) has high plan closure values, but relatively low profile closure values; the second group (cirques 3, 5, 7) has high profile closure values but relatively lower plan closure values; and the third group (cirques 4, 6 and 6a) has highest values for both plan and profile closure which might reflect higher degree of development (Figure 7C).
The length (L) of the main cirques averages at 838 m, the width (W) averages 806 m and the height (H) averages 288 m.Power exponents of logarithmic regressions of cirque length, width and height against cirques size indicate that cirque growth was allometric (Figure 7D), with cirques growing faster in horizontal than vertical dimensions.Similar findings have been reported from many regions [37] and have been attributed to periglacial processes and/or headwall recession having higher influence over subglacial deepening [56].The morphometric values of the main cirques on Jablanica Mt. (Table 1) are comparable to other regions of the world [37], and the main dimensional parameters (area, L/W, L/H, W/H) are very similar to the cirques from other mountains in the region (e.g., Northern Greece [53,54]; Retezat Mts., Romania [55]).Cirque size, as expected, generally decreases with increase of cirque floor elevation (Figure 7B).
All of the cirques have relatively high plan closure values (>174°), and based on both plan and profile closure, three groups can be delineated: the first group (cirques 1, 2, 8, 8a, 9 and 10) has high plan closure values, but relatively low profile closure values; the second group (cirques 3, 5, 7) has high profile closure values but relatively lower plan closure values; and the third group (cirques 4, 6 and 6a) has highest values for both plan and profile closure which might reflect higher degree of development (Figure 7C).
The length (L) of the main cirques averages at 838 m, the width (W) averages 806 m and the height (H) averages 288 m.Power exponents of logarithmic regressions of cirque length, width and height against cirques size indicate that cirque growth was allometric (Figure 7D), with cirques growing faster in horizontal than vertical dimensions.Similar findings have been reported from many regions [37] and have been attributed to periglacial processes and/or headwall recession having higher influence over subglacial deepening [56].Although all of the cirques have at least some of the textbook attributes of a glacial cirque (e.g., generally high plan and profile closure), most of them have some attributes that are not so well developed or are lacking (e.g., threshold and poorly developed floor).Six primary (3, 4, 5, 6, 7 and 8) and the two secondary (6a and 8a) cirques are classified as "definite" cirques, three of which (5,6,8) also have glacial lakes.They have the morphology of a well-developed cirque, supported by higher Although all of the cirques have at least some of the textbook attributes of a glacial cirque (e.g., generally high plan and profile closure), most of them have some attributes that are not so well developed or are lacking (e.g., threshold and poorly developed floor).Six primary (3, 4, 5, 6, 7 and 8) and the two secondary (6a and 8a) cirques are classified as "definite" cirques, three of which (5,6,8) also have glacial lakes.They have the morphology of a well-developed cirque, supported by higher values for plan and profile closure, but some are lacking threshold (e.g., cirques 3, 7) or have somewhat lower values for profile (e.g., cirques 8, 8a) or plan closure (e.g., cirques 3, 5, 7).The other four primary cirques (1, 2, 9, 10) are classified as "poor" cirques, because they generally have lower profile (although high plan) closure values, poorly developed floor and are lacking a well-developed threshold.

Glacier Reconstruction
The cirque boundaries are well defined, and the valleys are steep and glacially shaped with several glacial steps, and thresholds.Most of the ridges above 1900 m a.s.l.exhibit a smoothed morphology (Figure 4A), probably reflecting the presence of a former icefield.Accordingly, we suggest that glaciers of the reconstructed maximum ice extent (rMIE) likely originated in a large icefield covering most of the area above 1900 m a.s.l.elevation (Figure 8).The largest glaciers reached the foothills but did not extend to the piedmont.During our research no evidence was found of the existence of a previously suggested [5,6] larger glacier phase, with ice lobes extending down to 700-750 m a.s.l.
Glacier reconstruction of the rMIE glaciers (Vevčani Member/phase) suggests that the glaciers occupied an area of 22.6 km 2 within the studied middle part of Jablanica Mt.The longest glacier was in the Labuništa system with a length of 7.2 km (Figure 8).

Glacier Reconstruction
The cirque boundaries are well defined, and the valleys are steep and glacially shaped with several glacial steps, and thresholds.Most of the ridges above 1900 m a.s.l.exhibit a smoothed morphology (Figure 4A), probably reflecting the presence of a former icefield.Accordingly, we suggest that glaciers of the reconstructed maximum ice extent (rMIE) likely originated in a large icefield covering most of the area above 1900 m a.s.l.elevation (Figure 8).The largest glaciers reached the foothills but did not extend to the piedmont.During our research no evidence was found of the existence of a previously suggested [5,6] larger glacier phase, with ice lobes extending down to 700-750 m a.s.l.
Glacier reconstruction of the rMIE glaciers (Vevčani Member/phase) suggests that the glaciers occupied an area of 22.6 km 2 within the studied middle part of Jablanica Mt.The longest glacier was in the Labuništa system with a length of 7.2 km (Figure 8).The ice thickness reached its maximum values (200 m) in the cirque and middle part of the Podgorci and Labuništa valleys (Figure 9).In the lower reach of the Vevčani Valley the thickness of the reconstructed glacier is small (down to zero).This may be an artefact of the reconstruction, as the lateral moraines outlined by the spots of thinned glaciers were not yet filling the valley during this phase.The ice thickness reached its maximum values (200 m) in the cirque and middle part of the Podgorci and Labuništa valleys (Figure 9).In the lower reach of the Vevčani Valley the thickness of the reconstructed glacier is small (down to zero).This may be an artefact of the reconstruction, as the lateral moraines outlined by the spots of thinned glaciers were not yet filling the valley during this phase.The estimated ELAs during this phase varied between 1786 m and 1819 m a.s.l.if calculated by the ssAAR method and between 1756 m and 1809 m by the AABR method.The average values are 1809 m and 1777 m, respectively (Table 2; Figure 8).

Discussion
The glacial landforms mapped on the eastern side of the Jablanica Mt. supported the distinction of six deglaciation phases.Cirque morphometry parameters of the studied ten primary and two secondary cirques suggest that majority of the cirques are not "well-defined" but fall in the "definite" or "poor" classes [41].
The mapped glacier morphology suggest that subglacial deepening was limited within the cirques and was more intensive in the valley sections during more extensive glacial phases.The lack of developed cirque thresholds and poorly developed cirque floors suggest that cirque glaciers existed for a relatively short period or they were not erosive enough to create well-developed cirques [37].
For a tentative estimate of the timing of the maximum glacial advance, the ELAs of the reconstructed glaciers in the Jablanica Mt. were compared to ELAs calculated for surrounding The estimated ELAs during this phase varied between 1786 m and 1819 m a.s.l.if calculated by the ssAAR method and between 1756 m and 1809 m by the AABR method.The average values are 1809 m and 1777 m, respectively (Table 2; Figure 8).

Discussion
The glacial landforms mapped on the eastern side of the Jablanica Mt. supported the distinction of six deglaciation phases.Cirque morphometry parameters of the studied ten primary and two secondary cirques suggest that majority of the cirques are not "well-defined" but fall in the "definite" or "poor" classes [41].
The mapped glacier morphology suggest that subglacial deepening was limited within the cirques and was more intensive in the valley sections during more extensive glacial phases.The lack of developed cirque thresholds and poorly developed cirque floors suggest that cirque glaciers existed for a relatively short period or they were not erosive enough to create well-developed cirques [37].
For a tentative estimate of the timing of the maximum glacial advance, the ELAs of the reconstructed glaciers in the Jablanica Mt. were compared to ELAs calculated for surrounding mountain ranges of the Balkan Peninsula where geochronological data (or at least age estimates) are available (Figure 10).The timing of the maximum ice extent varied between the LGM (Šar, Rila, Retezat Mts.), the MIS 6 (Prokletije Mt.) and MIS 12 (Orjen, Durmitor, Pindus Mts.) phases in the

Conclusions
Mapping of glacial landforms and deposits offered the opportunity for the first reconstruction of the Pleistocene glaciers of the Jablanica Mt. using state-of-art GIS supported tools.A revised glacio-geomorphological map was produced with six members of the Jablanica Glacial Formation identified based on morphostratigraphical criteria.
The morphometric analysis of the cirques showed that most of the cirques on Jablanica Mt., although having at least some of the textbook attributes of a glacial cirque, are not well developed.
The wide and open cirques together with deep and stepped glacial valleys with several thresholds suggest that glacial erosion during the most extended glacial phases was focused along the valley axis to a position closer to the former ELA.The development stage of the cirques may also reflect more dynamic glaciers in the deeper, better developed cirques and less dynamic and/or shorter-lived glaciers in the poorly developed cirques.
The glacier reconstruction suggests that glaciers during the maximum ice extent formed an icefield covering most of the range above 1900 m a.s.l.elevation.The glacier cover for this largest stage could be estimated to 22.6 km 2 with an ELA around 1800 m a.s.l. in the studied valleys of central Jablanica Mt.The existence of a former icefield with several glacier tongues in the valleys is in agreement with results of the morphometric analysis suggesting enhanced glacial erosion along valley axes.
The timing of these former glacial stages remains an open question.However, the revised glacio-geomorphological map and the comparison of the calculated ELA with published ELA data of the region, tentatively assigns the maximum ice extent of the Jablanica Mt. to MIS 6.The presented geomorphological map and glacier reconstruction will serve as a valuable foundation for future efforts to establish a glacial chronological framework on Jablanica Mt.

Figure 2 .
Figure 2. General geological map of the study area on Jablanica Mountain.

Figure 2 .
Figure 2. General geological map of the study area on Jablanica Mountain.

Figure 3 .
Figure 3. Glacial geomorphological map of the middle part of Jablanica Mountain.Blue numbers indicate the locations of the photos on Figures 4 and 5.
Figure 3. Glacial geomorphological map of the middle part of Jablanica Mountain.Blue numbers indicate the locations of the photos on Figures 4 and 5.

Figure 3 .
Figure 3. Glacial geomorphological map of the middle part of Jablanica Mountain.Blue numbers indicate the locations of the photos on Figures 4 and 5.
Figure 3. Glacial geomorphological map of the middle part of Jablanica Mountain.Blue numbers indicate the locations of the photos on Figures 4 and 5.

Figure 4 .
Figure 4. Glacial morphology and deposits on Jablanica Mountain.(A) The uppermost part of Podgorci Valley having ridges with glacially smoothed morphology; Large glacially transported boulders of Jurassic conglomerate (Vevčani Member) in the lower parts of Vevčani (B) and Podgorci (C) valleys; (D) Large boulders of Jurassic conglomerate in the middle part of Vevčani Valley (Kutel Member); (E) Fluvioglacial deposits in the lower parts of Podgorci Valley; (F) Erratic boulder of Jurassic conglomerate (1.5 m high) laying on carbonate bedrock at Branica locality, with view of Pupuljak Valley in the back; (G) Large Jurassic conglomerate boulders of Leništa Member in Labuništa-Pupuljak Valley at Leništa locality.For location of the photos refer to Figure 3.

Figure 4 .
Figure 4. Glacial morphology and deposits on Jablanica Mountain.(A) The uppermost part of Podgorci Valley having ridges with glacially smoothed morphology; Large glacially transported boulders of Jurassic conglomerate (Vevčani Member) in the lower parts of Vevčani (B) and Podgorci (C) valleys; (D) Large boulders of Jurassic conglomerate in the middle part of Vevčani Valley (Kutel Member); (E) Fluvioglacial deposits in the lower parts of Podgorci Valley; (F) Erratic boulder of Jurassic conglomerate (1.5 m high) laying on carbonate bedrock at Branica locality, with view of Pupuljak Valley in the back; (G) Large Jurassic conglomerate boulders of Leništa Member in Labuništa-Pupuljak Valley at Leništa locality.For location of the photos refer to Figure 3.

Figure 5 .
Figure 5. Characteristic glacial and periglacial features of the study area.(A) Moraines of Jurassic conglomerate (Golina Member) in the middle part of Vevčani Valley at Golina locality; (B) Several moraine generations (yellow dashed lines, Lincura (1), Golina (2) and Leništa (3) Members) and rock glaciers (reworked glacial sediments; green dashed lines) in the Upper Podgorci Valley; (C) Moraines from Lincura (1) and Lokva (2) Members and rock glaciers in the Upper Vevčani Valley; (D) Moraine of the smallest glacier (Lincura Member) in the Upper Pupuljak Valley; (E) Moraine of the Lincura Member in the Upper Labuništa Valley (1) and Lateral moraine belonging to the Golina Member in the Upper Podgorci Valley (2).For location of the photos refer to Figure 3.

Figure 5 .
Figure 5. Characteristic glacial and periglacial features of the study area.(A) Moraines of Jurassic conglomerate (Golina Member) in the middle part of Vevčani Valley at Golina locality; (B) Several moraine generations (yellow dashed lines, Lincura (1), Golina (2) and Leništa (3) Members) and rock glaciers (reworked glacial sediments; green dashed lines) in the Upper Podgorci Valley; (C) Moraines from Lincura (1) and Lokva (2) Members and rock glaciers in the Upper Vevčani Valley; (D) Moraine of the smallest glacier (Lincura Member) in the Upper Pupuljak Valley; (E) Moraine of the Lincura Member in the Upper Labuništa Valley (1) and Lateral moraine belonging to the Golina Member in the Upper Podgorci Valley (2).For location of the photos refer to Figure 3.

Figure 6 .
Figure 6.Cirques boundaries and cirque floors mapped in the studied area.

Figure 6 .
Figure 6.Cirques boundaries and cirque floors mapped in the studied area.

Figure 8 .
Figure 8. Reconstructed glaciers during Vevčani phase in the Jablanica Mt.Solid red line indicates the size-specific Accumulation Area Ratio (ssAAR) equilibrium-line altitude (ELA), dashed red line indicates the Area-Altitude Balance Ratio (AABR) ELA.

Figure 8 .
Figure 8. Reconstructed glaciers during Vevčani phase in the Jablanica Mt.Solid red line indicates the size-specific Accumulation Area Ratio (ssAAR) equilibrium-line altitude (ELA), dashed red line indicates the Area-Altitude Balance Ratio (AABR) ELA.

Figure 9 .
Figure 9. Ice thickness of the reconstructed glaciers during Vevčani phase in the Jablanica Mt.

Figure 9 .
Figure 9. Ice thickness of the reconstructed glaciers during Vevčani phase in the Jablanica Mt.

Table 2 .
Results from the glacial reconstruction of the Vevčani phase (reconstructed maximum ice extent (rMIE)) on Jablanica Mt.

Table 2 .
Results from the glacial reconstruction of the Vevčani phase (reconstructed maximum ice extent (rMIE)) on Jablanica Mt.