Late Quaternary glacial phases in the Iberian Peninsula

Manuscript Details Manuscript number EARTH_2018_636 Title Late Quaternary glacial phases in the Iberian Peninsula Article type Review Article Abstract The only glaciers existing now in the Iberian Peninsula are small features located in the Pyrenees, though their number and extension has undergone significant changes over the Late Quaternary. The wide range of glacial landforms and deposits distributed across different Iberian ranges suggests the occurrence of several past stages with larger glacial systems. The objective of this research is to summarize the current knowledge on the spatial and temporal patterns of glacial activity in the Iberian mountains during the Late Quaternary. To this purpose, the chronological framework was divided in six periods: glaciations prior to the Last Glacial Cycle (Middle Pleistocene), Last Glacial Cycle (Late Pleistocene), Termination-1, Holocene, Little Ice Age (LIA) and present-day. The data were geographically divided considering the mountain systems where glacial evidence exists: Pyrenees, Cantabrian Range, NW ranges, Central Range, Iberian Range and Sierra Nevada. During cold Quaternary stages, ice accumulated in the head valleys of these mountain ranges and glaciers flowed down-valleys. In all cases, glaciers remained confined within the mountain systems and did not reach the surrounding lowlands. Depending on the combination of temperatures and moisture conditions, more or less ice was stored. In some ranges, there is evidence of Middle Pleistocene glaciations, one potentially correlating with marine isotope stage (MIS) 12 and another correlating with MIS 6 with glaciation dated to ca. 130-170 ka. However, most of the glacial records correspond to the Last Glacial Cycle and subsequent Termination. The maximum glacial expansion of this last Pleistocene glaciation stage occurred well before the global Last Glacial Maximum (LGM) between 30 and 60 ka in the Cantabrian Mountains and Pyrenees, at ca. 30 ka in Sierra Nevada and NW ranges, and (almost) synchronously to the LGM in the Central Range and Iberian Range. A massive glacial retreat occurred in all ranges at 19-20 ka, but the long-term deglaciation process was interrupted by cold stages, such as the Oldest and Younger Dryas, which favoured glacial expansion in the highest mountains. Temperature increase recorded during the Holocene conditioned the melting of glaciers, which only reappeared in the highest massifs during the coldest stages, such as the LIA. However, post-LIA warming led to glacier disappearance in the Cantabrian Mountains, Sierra Nevada and most massifs of the Pyrenees, together with an accelerated shrinking of the small glaciers still existing in this range at elevations near 3000 m. Keywords Iberian Peninsula, glaciation, Last Glacial Maximum, Termination-1, Holocene, Little Ice Age. Corresponding Author Marc Oliva Corresponding Author's Institution Department of Geography, University of Barcelona, Spain Order of Authors Marc Oliva, palacios david, José M. Fernández-Fernández, Laura RodríguezRodríguez, Jose-Maria Garcia Ruiz, Nuria Andrés, Rosa M. Carrasco, Javier Pedraza, Augusto Perez Alberti, Marcos Valcarcel, Philip. Hughes Suggested reviewers Jesus Ruiz-Fernández, Susan Ivy-Ochs, Jamie Woodward, Mauro Guglielmin Submission Files Included in this PDF File Name [File Type] Cover letter.pdf [Cover Letter] Abstract.docx [Abstract] Manuscript.docx [Manuscript File] Figure 1.tif [Figure] Figure 2.tif [Figure] Figure 3.tif [Figure] Figure 4.tif [Figure] Figure 5.tif [Figure] Figure 6.tif [Figure] Figure 7.tif [Figure] Figure 8.jpg [Figure] Figure 9.tif [Figure] Figure 10.tif [Figure] Figure 11.tif [Figure] Figure 12.tif [Figure] Figure 13.tif [Figure] Figure 14.tif [Figure] Tables.docx [Table] Submission Files Not Included in this PDF File Name [File Type] Supplementary material.kmz [Interactive Map Data (.kml, .kmz)] To view all the submission files, including those not included in the PDF, click on the manuscript title on your EVISE Homepage, then click 'Download zip file'. This manuscript contains content innovation file(s). General instructions for reviewing content innovation files can be found here. Abstract The only glaciers existing now in the Iberian Peninsula are small features located in the Pyrenees, though their number and extension has undergone significant changes over the Late Quaternary. The wide range of glacial landforms and deposits distributed across different Iberian ranges suggests the occurrence of several past stages with larger glacial systems. The objective of this research is to summarize the current knowledge on the spatial and temporal patterns of glacial activity in the Iberian mountains during the Late Quaternary. To this purpose, the chronological framework was divided in six periods: glaciations prior to the Last Glacial Cycle (Middle Pleistocene), Last Glacial Cycle (Late Pleistocene), Termination-1, Holocene, Little Ice Age (LIA) and present-day. The data were geographically divided considering the mountain systems where glacial evidence exists: Pyrenees, Cantabrian Range, NW ranges, Central Range, Iberian Range and Sierra Nevada. During cold Quaternary stages, ice accumulated in the head valleys of these mountain ranges and glaciers flowed downvalleys. In all cases, glaciers remained confined within the mountain systems and did not reach the surrounding lowlands. Depending on the combination of temperatures and moisture conditions, more or less ice was stored. In some ranges, there is evidence of Middle Pleistocene glaciations, one potentially correlating with marine isotope stage (MIS) 12 and another correlating with MIS 6 with glaciation dated to ca. 130-170 ka. However, most of the glacial records correspond to the Last Glacial Cycle and subsequent Termination. The maximum glacial expansion of this last Pleistocene glaciation stage occurred well before the global Last Glacial Maximum (LGM) between 30 and 60 ka in the Cantabrian Mountains and Pyrenees, at ca. 30 ka in Sierra Nevada and NW ranges, and (almost) synchronously to the LGM in the Central Range and Iberian Range. A massive glacial retreat occurred in all ranges at 19-20 ka, but the long-term deglaciation process was interrupted by cold stages, such as the Oldest and Younger Dryas, which favoured glacial expansion in the highest mountains. Temperature increase recorded during the Holocene conditioned the melting of glaciers, which only reappeared in the highest massifs during the coldest stages, such as the LIA. However, post-LIA warming led to glacier disappearance in the Cantabrian Mountains, Sierra Nevada and most massifs of the Pyrenees, together with an accelerated shrinking of the small glaciers still existing in this range at elevations near 3000 m. Key words: Iberian Peninsula, glaciation, Last Glacial Maximum, Termination-1, Holocene, Little Ice Age. Abstract The only glaciers existing now in the Iberian Peninsula are small features located in the Pyrenees, though their number and extension has undergone significant changes over the Late Quaternary. The wide range of glacial landforms and deposits distributed across different Iberian ranges suggests the occurrence of several past stages with larger glacial systems. The objective of this research is to summarize the current knowledge on the spatial and temporal patterns of glacial activity in the Iberian mountains during the Late Quaternary. To this purpose, the chronological framework was divided in six periods: glaciations prior to the Last Glacial Cycle (Middle Pleistocene), Last Glacial Cycle (Late Pleistocene), Termination-1, Holocene, Little Ice Age (LIA) and present-day. The data were geographically divided considering the mountain systems where glacial evidence exists: Pyrenees, Cantabrian Range, NW ranges, Central Range, Iberian Range and Sierra Nevada. During cold Quaternary stages, ice accumulated in the head valleys of these mountain ranges and glaciers flowed downvalleys. In all cases, glaciers remained confined within the mountain systems and did not reach the surrounding lowlands. Depending on the combination of temperatures and moisture conditions, more or less ice was stored. In some ranges, there is evidence of Middle Pleistocene glaciations, one potentially correlating with marine isotope stage (MIS) 12 and another correlating with MIS 6 with glaciation dated to ca. 130-170 ka. However, most of the glacial records correspond to the Last Glacial Cycle and subsequent Termination. The maximum glacial expansion of this last Pleistocene glaciation stage occurred well before the global Last Glacial Maximum (LGM) between 30 and 60 ka in the Cantabrian Mountains and Pyrenees, at ca. 30 ka in Sierra Nevada and NW ranges, and (almost) synchronously to the LGM in the Central Range and Iberian Range. A massive glacial retreat occurred in all ranges at 19-20 ka, but the long-term deglaciation process was interrupted by cold stages, such as the Oldest and Younger Dryas, which favoured glacial expansion in the highest mountains. Temperature increase recorded during the Holocene conditioned the melting of glaciers, which only reappeared in the highest massifs during the coldest stages, such as the LIA. However, post-LIA warming led to glacier disappearance in the Cantabrian Mountains, Sierra Nevada and most massifs of the Pyrenees, together with an accelerated shrinking of the small glaciers still existing in this range at elevations near 3000 m. Key words: Iberian Peninsula, glaciation, Last Glacial Maximum, Termination-1, Holocene, Little Ice Age. 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 1Introduction The landscape of the highest Iberian mountain ranges has been intensely shaped by both Quaternary glaciations and postglacial periglacial processes prevailing during interglacial periods (Oliva et al., 2016). Until the mid19th century, the biblical Great Flood theory was widely accepted to explain the existence of glacial deposits and landforms in mountain areas around the world (Ehlers et al., 2016). The first notes including a scientific perspective were carried out in the Alps by Agassiz (1840), who described the impact on Alpine landscapes of large glaciers flowing down-valleys during ancient ice ages. Taking into account these observations, early scientists from Central Europe (biologists, geologists, geographers, etc.) – who had previously conducted research in the Alps – exported their knowledge to southern European ranges, such as Iberian mountains (i.e. Boissier, 1839), describing similar geomorphological features and processes to those observed in their home countries (Gómez-Ortiz et al., 2018). The visits of Albrecht Penck – a geomorphologist with a long previous field experience on the study of glacial features in the Alps –, to various Iberian ranges allowed the identification of the size and distribution of past glaciers, for the very first time, as well as the existence of geomorphic features of several glacial cycles (see, for example, Penck, 1883). In a parallel way, the Quaternary scientist Hugo Obermaier (1877-1946) promoted the study of Iberian glaciations from the beginning of the 20th century. In the Iberian Peninsula, early reports on the role of glaciers shaping the high lands focused on the still glaciated massifs at the end of the Little Ice Age (LIA), namely the Pyrenees (see González-Trueba et al., 2008), Cantabrian Mountains (see González-Trueba, 2006, 2007) and Sierra Nevada (see Gómez-Ortiz et al., 2006, 2009). The first texts were accompanied by geographical sketches, paintings and photographs of those glaciers (e.g. Bide, 1893; Schrader, 1895; Briet, 1902), together with accurate descriptions of their topographical characteristics, elevations and dimensions. Although some researchers described cold-climate geomorphological features in the main massifs over the first half of the 20th century (Panzer, 1926; GarcíaSainz, 1935; Dresch, 1937; Nussbaum, 1956), the turning point for glacial and periglacial research in the Iberian Peninsula was the organization of the International Union for Quaternary Research (INQUA) meeting in Barcelona-Madrid in 1957 (Gómez-Ortiz and Palacios, 1995; Gómez-Ortiz and Vieira, 2006). This conference was a major boost for research on past and present glacial and periglacial processes, favouring networking and promoting internationalization, which resulted in the publication of key studies over the next decades (Barrère, 1963; Messerli, 1967; Serrat, 1977; García-Ruiz, 1979; Pérez-Alberti, 1979; Gómez-Ortiz, 1980; Vilaplana, 1983; Bru, 1985; Ortigosa, 1986). Most of these works included new advances on the monitoring of geomorphic processes, sedimentological analysis of glacial and periglacial deposits as well as new observations in unexplored areas that opened new perspectives on the impact of Quaternary climate variability on Iberian mountain landscapes. For some ranges, some authors even proposed the existence several glacial stages based on the existence of landforms and processes left by past glaciers, such as in the Pyrenees where three glacial cycles were described (Penck, 1883). The relative chronological reconstruction of glacial stages proposed in some of these studies was progressively complemented with the implementation of absolute dating techniques – namely radiocarbon dating – that allowed placing in time environmental events, and therefore providing ages for past glacial activity. The first works including radiocarbon ages were based on organic remnants trapped within glacial sediments that provided evidence on glacial advances occurred during the Last Glacial Cycle (Late Pleistocene) until ca. 40 ka cal BP. This technique was subsequently complemented with the use of Optically Stimulated Luminescence (OSL) dating on fluvio-glacial sediments in the late 1990s and early 2000s, which allowed expanding the chronology until ca. 80-90 ka BP (García-Ruiz et al., 2010). However, the use of these two dating methods generated a deep discussion on the chronology of the local Maximum Ice Extent (MIE) of the Last Glacial Cycle in the Iberian mountains (Pérez-Alberti et al., 2004; Hughes et al., 2006a; Hughes and Woodward, 2008; García-Ruiz et al., 2010), which was even more intense when the use of surface exposure dating using CosmicRay Exposure (CRE) dating became widespread (Pallàs et al., 2006; García-Ruiz et al., 2010; Palacios et al., 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 2011, 2018). The global Last Glacial Maximum is defined as the interval 26.5-20/19 ka (Clark et al., 2009) or 27.5-23.3 ka (Hughes and Gibbard, 2015), with both these overlapping definitions falling within marine isoptope stage (MIS) 2. Whereas chronologies based on 14C and OSL suggested a local MIE predating the global LGM in most ranges (García-Ruiz et al., 2003, 2012; Lewis et al., 2009; Ruiz-Fernández et al., 2016; Serrano et al., 2013, 2015, 2017), CRE dates approached the age of the MIE in the Pyrenees and Central Range mountains closer to the LGM timing (Pallàs et al., 2006; Palacios et al., 2011, 2012a,b, 2015; DomínguezVillar et al., 2013) but confirmed a MIE older than the LGM in the NW ranges (Rodríguez-Rodríguez et al., 2011, 2014), Cantabrian Mountains (Rodríguez-Rodríguez et al., 2015, 2016) and Sierra Nevada (GómezOrtiz et al., 2012a, 2015; Palacios et al., 2016). Recent advances on surface exposure dating have also favoured the chronological reconstruction of older and younger glaciations in some mountain ranges (García-Ruiz et al., 2014, 2016; Palacios et al., 2017a, b, 2018), which is crucial to better understand Quaternary climate variability in southern Europe. With this large increase in the number of studies focusing on past glaciations in the Iberian Peninsula, the objective of this work is to review the different glacial stages occurred in the Iberian mountains from a spatiotemporal perspective. We have reviewed all available dates (14C, OSL, U-Th series, 210Pb, and historical sources), as well as unified the criteria used from several authors in different massifs including the last advances on the production rate of CRE (namely 10Be, 36Cl and 21Ne) with the purpose of giving answer to the following questions: - What sort of geomorphic evidence is there from glaciations occurred prior to the Last Glacial Cycle in the Iberian Peninsula? - What is the exact timing of the local MIE in the different Iberian mountain ranges? - What do we know about the chronology and the impact of glacial stages following the post-LGM massive deglaciation on the current landscape of these massifs? - Are these glacial advances and retreats synchronous to patterns occurred in other high mountain ranges from Europe and northern Africa? - What is the main factor (temperature vs moisture) controlling the major glacial expansion and subsequent advances and retreats in the Iberian Peninsula? - Are results from the different dating methods used until now comparable? - What are the temporal and spatial gaps on our current understanding of Late Quaternary glacial processes in the Iberian Peninsula? 2Study area Extending over an area of 582,925 km², the Iberian Peninsula is located in the SW corner of Europe between latitude 43° 47′ N and 36° 01′ N and longitude 9º 30’ W and 3º 19’. Mountain ranges in the Iberian Peninsula are generally aligned W-E and distributed in the periphery, separating the relatively flat areas of the central part of the peninsula from the surrounding coastal fringes. This rough topography, together with the differing influences affecting the Iberian Peninsula, such as the maritime (Atlantic/Mediterranean), climatic (subtropical high pressure belt/mid-latitude westerlies) and biomes (Europe/Africa), result in the variety of landscapes existing across Iberia (Oliva et al., 2018). The Iberian Peninsula includes six mountain ranges with peaks above 2000 m a.s.l.: the Pyrenees, the Cantabrian Mountains, the NW ranges, the Central Range, the Iberian Range and the Betic Range. Glacial landscape features are widespread in these mountains, though evidence of Quaternary glacial activity is also found in other mountains at lower altitudes, particularly in the NW corner (Figure 1). Figure 1 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 Iberian climate is highly variable and affected by continental and maritime air masses from distinct origins. The complex terrain determines the existence of a wide range of microclimatic regimes across the Iberian Peninsula and within the individual mountain systems. Whereas precipitation is mostly concentrated between October and May brought by the mid-latitude cyclones associated to the prevailing zonal circulation, the summer is generally a dry season due to the influence of the Azores anticyclone (Trigo et al., 2004). Consequently, precipitation in the Iberian Peninsula decreases generally from N to S and from W to E, with values over 2000-2500 mm in Atlantic-influenced ranges and minima of 600-900 mm in the Sierra Nevada. Mean annual air temperatures (MAAT) follow an opposite pattern, increasing towards the S and the E. The regional 0 °C isotherm is placed at ca. 2400-2500 m in the Cantabrian Mountains (Muñoz, 1982), ca. 2950 m in the Pyrenees (Chueca et al., 2005) and at 3400 m in the Sierra Nevada (Oliva et al., 2016b). The highest lands of these mountain ranges have been intensely shaped by glacial processes during the Quaternary, as well as by postglacial environmental dynamics, namely by periglacial, slope and alluvial processes and shallowand deep-seated landslides (Oliva et al., 2016a). Quaternary cold stages (glacial) favoured the development of glaciers filling the valley heads and flowing down-slope. The combination of temperatures and moisture conditions controlled the elevation shifts of the Equilibrium Line Altitude (ELA), and therefore the ice volume stored in the Iberian mountains and the length of the glaciers. Apart from climate conditions prevailing in the North Atlantic region, the latitude as well as the geographical influence of sea surface water temperatures (cool Atlantic Ocean vs warm Mediterranean Sea) determined the spatial distribution of the glaciated domain in the Iberian ranges during Quaternary glacial stages, generally increasing towards the S and E (Pérez-Alberti et al., 2004). Ice-free areas located above the ELA and below the glaciated environments were affected by very active periglacial processes, with the formation of permafrost landforms and seasonal frost features that are inactive under present-day climate regime (Oliva et al., 2016a). Quaternary warm periods (interglacial) promoted the migration of the periglacial belt to higher areas, and cryonival dynamics reshaped the formerly glaciated environments. This is what is occurred in the present-day interglacial, the Holocene, and very small remnants of Quaternary glaciers currently exist only in the Pyrenees and the highest massifs are affected by periglacial activity. 3Methodology This paper presents a thorough review of all existing scientific literature on glacial processes in the Iberian mountains, including research papers published in international peer-reviewed journals, book chapters and conference proceedings, theses and books as well as other regional publications published in local languages. With the purpose of better understanding the spatial and temporal patterns of glacial activity in the Iberian Peninsula, chronological and geomorphological data were compiled and divided considering the different mountain ranges: Pyrenees, Cantabrian Mountains, NW ranges, Central Range, Iberian Range and Betic Range (i.e. Sierra Nevada). In addition, for each study area, data are organized based on six main periods: glaciations prior to the Last Glacial Cycle (Middle Pleistocene), Last Glacial Cycle (Late Pleistocene), Termination-1, Holocene, LIA and present-day. In this paper these intervals are informally defined and are not intended to indicate or replace any existing formal stratigraphical basis, as the intervals sometimes represent transitions between major climatic changes. Several definitions can be used to define the Last Glacial Cycle (Hughes and Gibbard, 2018). In this paper, the Last Glacial Cycle is defined by the marine isotopic stage boundaries and includes MIS 5d-2 starting at the end of Eemian period, at 115 ka (Dahl-Jensen et al., 2013). The onset of global climate changes leading to Termination-1 (which marks the boundary between the Late Pleistocene and the Holocene; Hughes and Gibbard, 2018) started at 19-20 ka when a widespread glacial retreat is detected across the Northern Hemisphere (Clark et al., 2009). This includes several cold and warm stages until the Holocene, namely: the Oldest Dryas (OD; 17.5-14.7 ka, stadial GS-2.1a), Bølling-Allerød (BA; 14.7-12.9 ka, interstadial GI-1) and 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 Younger Dryas (YD; 12.9-11.7 ka, stadial GS1). All of these intervals correspond to the classical Late-glacial in northwest Europe (Mangerud, 1974). Following Walker et al. (2012), the Holocene is subdivided in the Early (11.7-8.2 ka BP), Middle (8.2-4.2 ka) and Late Holocene (since 4.2 ka BP), which includes the LIA that in the Iberian Peninsula extends from 1300 to 1850 (Oliva et al., 2018). Finally, we also include a description on present-day geomorphological processes to better frame post-LIA environmental dynamics. Table 1 Figure 2 The different glacial phases have been inferred using different dating methods, including 351 (36Cl, 10Be and 21Ne) CRE dates, 68 14C, 43 OSL and 9 other (such as U-Th series, 210Pb, and historical sources). Their distribution considering the diverse time periods and mountain ranges is summarized in Table 1 and Figure 2. All 14C radiocarbon dates have been calibrated using the CALIB 7.1 program. Also, in this study, CRE ages have been re-calculated for those samples from which enough information is available so that a new reassessment can be done. In the cases where required data is not available for recalculation, the original published age is indicated (Supplementary material). The 36Cl CRE ages have been recalculated applying the same parameters aiming to achieve comparable results, so the following 36Cl production rates have been implemented: 42.2 ± 4.8 atoms 36Cl (g Ca)-1 yr-1 from Ca spallation (Schimmelpfennig et al., 2011), 148.1 ± 7.8 atoms 36Cl (g K)-1 yr-1 from K spallation (Schimmelpfennig et al., 2014), 13 ± 3 atoms 36Cl (g Ti)-1 yr-1 from Ti spallation (Fink et al., 2002) and 1.9 atoms 36Cl (g Fe)-1 yr-1 from Fe spallation (Stone et al., 2005). A value of 696 ± 185 neutrons (g air)-1 yr-1 was applied as the production rate of the epithermal neutrons from fast neutrons in the land/atmosphere interface (Marrero et al., 2016). Scaling factors for nucleonic and mounic production were recalculated following the formulae in Stone (2000). On the other hand, 10Be CRE ages have been recalculated by using the “CREp” (Cosmic Ray Exposure Program) online calculator (Martin et al., 2017; http://crep.crpg.cnrs-nancy.fr/#/). With the aim of unify the exposure age calculation for all samples, we have applied the LSD (Lifton-Dunai-Sato) scaling scheme (Lifton et al., 2014), the ERA40 atmospheric model (Uppala et al., 2005) and the geomagnetic database based on the LSD Framework (Lifton et al., 2014). Applying the aforementioned parameters implies, in turn, a SLHL (Sea-Level High-Latitude) 10Be production rate from Be spallation of 3.99 ± 0.22 atoms g-1 yr-1. Topographic shielding factor of each 10Be/36Cl sampling site has been re-calculated through the Topographic Shielding Factor through the “Topographic Shielding Calculator v.2” belonging to the “CRONUSCalc” Program (Marrero et al., 2016). However, for those samples whose field measurements for topographic shielding factor calculation are unreliable or unavailable, it has been obtained from the “Point-based Shielding Model” GIS-tool devised by Li (2018), which implements the method proposed by Balco et al. (2008). It only requires a Digital Elevation Model, a point shapefile with the location of the samples and the dipping and strike data stored in two separate fields. For each mountain range, all the data were summarized in a table including available chronology (if existing), geomorphic evidence (environments, landforms and processes), main references as well as a figure with features representative of each of the stages. The different time periods were mapped in GIS environment to better represent the spatio-temporal patterns of glacial activity in each mountain range. 4Results The study of glacial processes in the Iberian Peninsula has been solely focused on mountain ranges, where the impact of cold-climate geomorphological processes on the landscape is widespread as shown by the existence of a wide range of landforms and deposits described below. 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 4.1. The Pyrenees Quaternary glaciations left a deep imprint on the relief of the Pyrenees because of the large extent of glaciers, with ice tongues exceeding 30–50 km long and 400 m thick, and even reaching up to 800 m in the overdeepened basin of Benasque in the Ésera valley (Bordonau, 1992) as well as in the upper Garonne valley (Fernandes et al., 2017). Altitude and latitude played a major role to explain the size and geomorphological features of glaciers in the Pyrenees: Thus, altitude exceeds 2500 m in the main divides from the headwater of the Ansó valley eastwards, and frequently surpasses 3000 m between the headwaters of the Gállego valley and the Noguera Ribagorçana valley (Aneto, 3404 m; Posets, 3375 m; Monte Perdido: 3355 m). Besides, The Pyrenees are, together with the Cantabrian Mountains and the NW ranges, those located at higher latitudes in the Iberian Peninsula, resulting in a relatively low position of the 0 ºC isotherm and the ELA during the coldest periods. Here, more than in any other Mediterranean mountains, a large variety of well-developed landforms and deposits can be found: glacial cirques, U-shaped valleys (like, for instance, those of the Aragón Subordán, Estarrún, Aragón, Gállego, Ara, Cinca, Cinqueta, Ésera, Noguera Ribagorçana and Noguera Pallaresa valleys), verrous or rocky thresholds, hanging-tributary valleys (like the Ip valley, a tributary of the Aragón valley; the Arazas valley, a tributary of the Ara valley; and most of glacial tributaries of the Ésera paleoglacier), glaciolacustrine deposits in ice/moraine-dammed lakes (many examples, particularly the Linás de Broto paleolake dammed by a lateral moraine of the Ara glacier; Sancho et al., 2018), glacial transfluences, erratic blocks, glacial-origin lakes, roches moutonnées and proglacial plains (e.g. that located at the front of the Senegüé moraine, in the Gállego valley), as well as lateral and frontal moraines (like those of the Villanúa-Castiello de Jaca basin in the Aragón valley, the Senegüé moraine in the Gállego valley, the Sant Antoni moraine in the Noguera Ribagorçana valley, and the Puigcerdá complex at the Cerdanya plain), and flute moraines caused by surging glaciers (e.g. in the Marboré cirque, Monte Perdido massif: Serrano and Martín-Moreno, 2018). Most of the main glaciers started from cirques shaped in the paleozoic axis of the Pyrenees (granite, quartzite and limestone), although relatively big glaciers also started in the Inner Sierras, composed of Mesozoic and Cenozoic limestone and sandstone (García-Ruiz et al., 2000), affected by the Alpine tectonics. Their thickness decreased relatively rapid as clearly represented by the declining height of the lateral moraines near their terminal basins. The outermost moraines are located at approximately 800–900 m, and even less in valleys descending to the northern slope of the range (i.e. Garone valley glacier) and in some southern valleys (i.e. Gállego valley glacier). The Pyrenees include deposits from various glacial cycles. Studies on Pyrenean Quaternary glaciers account for a long tradition since the last decades of the 19th century, thanks to prestigious foreigner geomorphologists. This was the case of Penck (1883) and Panzer (1926). The first one visited the Aragón, Gállego and Ara valleys and raised the question of how many glaciations are represented in the morainic deposits of the Pyrenees; the second one recognized two glacial cycles in the terminal glacial basin of Castiello de Jaca-Villanúa, Aragón River valley, which were then attributed to Riss and Würm – following the Alpine terminology proposed by Penck (1883) – according to the connections between the main moraines and the 60 and 20 m fluvial terraces. The pioneer works of both geomorphologists were the basis for posterior studies focused for decades on: (i) examining the maximum extent recorded by the main Pyrenean glaciers; and (ii) discussing, in absence of direct dates, if the morainic deposits corresponded to one or more glacial cycles. Among such studies, the works of Obermaier (1921), Solé Sabarís (1941, 1951), Llopis-Lladó (1947), Nussbaum (1949), Barrère (1963), Gómez-Ortiz (1987), Bordonau (1992), Serrat et al. (1994), Martí-Bono and García-Ruiz (1994), Chueca et al. (1998) and Serrano (1998) described the position of the main glacial deposits and tried to establish a sequence of glacial periods based on the fabric characteristics and location of the deposits, e.g. Glacial Maximum, Intermediate Stable Stage, Disjunctive Stage and Finiglacial Stage (Martínez de Pisón, 1989). Recent studies carried out during the last two decades have decisively contributed to answer such questions and to pose new ones related with the age of the maximum extent of Pyrenean glaciers and of the distinct stages during deglaciation processes. Currently, the main features of the Late Pleistocene Pyrenean glacial cycle are relatively well known, although new chronological data are needed to establish better correlations between Pyrenean valleys with the glacial sequences identified in other Iberian and European ranges. 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 ages from erratic and moraine boulders from the terminal zone of the Porma paleoglacier provided results of 173–131 ka that suggest glacial activity during MIS 6 (Rodríguez-Rodríguez et al., 2016). Last Glacial Cycle During the MIE of the last glacial cycle, glaciers covered a total surface extent of about 2240 km2, showing remarkable asymmetry on glacier distribution between both sides of the range (Rodríguez-Rodríguez et al., 2015). The combination of pre-glacial topography, distance to the sea, and climate most likely controlled the distribution of glaciers, whose longest glacier tongues were shorter in the steeper northern slope (<15 km) compared to the gentler southern one of the range (15 to 51 km) (Serrano et al., 2017). In the central Cantabrian Mountains, where broad land areas are above 2000 m asl, some extensive ice fields formed close to the N-S divide showing asymmetric development towards the southern slope of the range, with the exception of Picos de Europa and Pas Mountains. Occasionally, outlet glaciers were thick enough to fulfill the valleys and overflowed towards adjacent tributaries (Figure 5), forming interconnected systems of valley glaciers or transection glaciers (e.g. the Porma paleoglacier; Rodríguez-Rodríguez et al., 2016). The front of the longest glaciers reached minimum elevations that were generally lower in the northern slope of the range (most frequently between 800 and 1200 m altitude, locally down to 400–500 m) compared to the southern one (1100– 1250 m, locally down to 725–905 m). Regional paleoELA varied from ca. 1000 to 2000 m during the MIE, showing a distribution pattern similar to the present-day winter precipitation and summer temperature (SantosGonzález et al., 2013). The lowest ELA values were recorded towards the eastern and western ends of the range (ca. 1000 m), while it remained at ca. 1500 and 1600 m in the Central Cantabrian Mountains divide. The highest ELA values (> 1750 m) were recorded further inland, in Peña Prieta Massif (Pellitero, 2013). Regarding the timing of glacial stages, the application of radiocarbon in glacio-lacustrine environments, deposited synchronously or subsequently to the MIE moraines, provided minimum ages in the range 45–35 14C cal ka BP for the Western and Central massifs of Picos de Europa (Jiménez-Sánchez and Farias, 2002; Moreno et al., 2010; Serrano et al., 2012; Nieuwendam et al., 2015; Ruíz-Fernández et al., 2016), > 44–35 14C cal ka BP in Somiedo-Laciana (Jalut et al., 2010), > 33.5 14C cal ka BP in the Redes Natural Park (JiménezSánchez and Farias, 2002), and > 30 cal ka BP in the Pas Mountains (Serrano et al., 2013). In the Cares valley, at least two fluvial terraces were deposited at 41.1–42.4 14C cal ka BP (+8–10 m) and >48 ka (+20–22 m) coevally with glacier occupation in the Picos de Europa (Ruíz-Fernández and Poblete-Piedrabuena, 2011). The OSL analysis performed in supraglacial tills from Vega del Naranco moraine complex, in Fuentes Carrionas Massif, provides a reference age of ca. 36 ka for the local MIE and involves younger ages for the six moraine ridges preserved upvalleys (Serrano et al., 2013). Similarly, the second outermost ridge of the Brañagallones lateral moraine complex, in the Redes Natural Park, yielded an OSL age of ca. 24 ka for a glacier re-advance of similar extent than the previous MIE (Jiménez-Sánchez et al., 2013). Numerical ages obtained in the Porma valley using a combination of 10Be CRE, OSL and 14C have provided the most complete record up to date of past glaciations in the Cantabrian Mountains (Rodríguez-Rodríguez et al., 2016, 2018a). Results suggest continuous glacial occupation of the valley during the Last Glacial Cycle, recording at least three major glacial stages at ca. 110 ka (MIS-5d), ca. 56 ka (MIS-3); and 33–24 ka (MIS-2) for which the glacier front reached minimum elevations of 1110–1130 m. The OSL analysis performed on till samples in the northern slope of Pas Mountains also suggest glacial occupation from 78.5 to 40 ka in the Ason valley and 45–41.6 ka in the Gandara valley (Frochoso et al., 2013). The glacial re-advances recorded in some Cantabrian valleys such as the Porma, Monasterio or Vega del Naranco are synchronous with the global LGM. Termination-1 The minimum 10Be CRE ages obtained from lateral moraines in the Porma catchment place the onset of the last Termination at ca. 22–20 ka (Rodríguez-Rodríguez et al., 2018a), and both 10Be CRE and 14C dating suggest considerably glacier thinning by ca. 18–17.5 ka (Rodríguez-Rodríguez et al., 2016, 2018a). RuizFernández et al (2016) showed evidence of intense periglacial conditions with thinner glaciers between 22.5 and 18.7 14C cal ka BP in the Western Massif of Picos de Europa. The sequence of recessional moraines dated 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 in the Monasterio valley provides consistent evidence for deglaciation in the opposite mountain slope (Rodríguez-Rodríguez et al., 2017), suggesting glacier front stagnations and/or minor re-advances at ca. 19, 17.5 and 14.6 10Be ka, for which the glacier front was located at altitudes between 1150 and 1540 m. In the northern slope of Fuentes Carrionas massif, proglacial lacustrine sedimentation occurred at 18.9–18.8 cal ka BP in Vega del Naranco (1540–1550 m altitude) due to valley impoundment by the moraine complex (Serrano et al., 2013). Westward, the minimum radiocarbon age obtained at the basis of Lago de Valle (also known as Lago de Ajo, placed at 1566 m altitude) suggests remarkable retreat of NW-oriented valley glaciers in Somiedo by ca. 17.6–17.1 14C cal ka BP (Allen et al., 1996). Thus, broad areas of the Central Cantabrian Mountains with top elevations between 1900 and 2100 m, like Monasterio valley in Redes Natural Park or Lago de Valle in Somiedo, were possibly fully deglaciated after the OD. A phase of glacier tongue individualization or disjunction with the fronts reaching 1600–1700 m is identified in some massifs, suggesting an ELA pattern similar to the local MIE, recording minima values of 1310 and 1500 m in the eastern and western ends of the range and increasing towards the middle part (Serrano et al., 2017). Besides an advanced retreat of valley glaciers, the Late-glacial was characterized by a progressive replacement from glacial to periglacial processes with development of large rock glaciers that could possibly be active intermittently during the Holocene and the LIA (Alonso and Trombotto-Liaudat, 2009). The distribution of relict rock glaciers shows initiation lines between 1600 and 2000 m and toes reaching minimum elevations of 1500 m (Gómez-Villar et al., 2011). Up to date, the 10Be CRE dating of two rock glacier toes located at altitudes of 1540 and 1650 m provided reference ages of ca. 16.2–13.6 ka for the stabilization of their lowest ridge (Rodríguez-Rodríguez et al., 2016, 2017). Paraglacial rock slope instabilities started during the OD, soon after glacier retreat (ca. 16–15 ka) as indicates preliminary numerical ages from Braña Creek in the San Isidro mountain pass area, and continued during the Holocene (Rodríguez-Rodríguez et al., 2018b). Although the occurrence of a YD re-advance of glaciers has been hypothesized for the highest massifs of the Cantabrian Mountains, such as the Picos de Europa and Fuentes Carrionas as well as local settings like Peña Ubiña, no dating evidence is yet available (Serrano et al., 2013, 2017). In the Central Massif of Picos de Europa, small glaciers reaching up to 1 km in length developed, leaving a thick debris cover on moraines at altitudes of 1800–2000 m in shelter areas of glacier cirques, suggesting that the local ELA was at about 2200 m altitude (Serrano et al., 2017). Similarly, moraines preserved in Fuentes Carrionas at altitudes between 1900– 2000 m and 2100–2200 m suggests two stages of glacier stabilization that could correlate with the YD for which the ELA was set at 2300 and 2350 m, respectively (Pellitero, 2013). However, broad areas of the Central Cantabrian Mountains where the highest peaks do not exceed 2000–2100 m remained glacier-free during the YD, possibly because the regional ELA was set at about 2500 m altitude (Serrano et al., 2017). Little Ice Age Only six glaciers existed in the highest cirques of Picos de Europa at the end of the LIA, extending over a surface of 25.5 ha (15.5 ha in the Central Massif and 10 ha in the Western Massif), and always distributed at the foot of northern walls in sheltered environments snow-fed by avalanches and wind action (Miotke, 1968; González-Trueba et al., 2008). During the mid-19th century these glaciers reached their maximum extent and built frontal moraines at 2200–2320 m, and the ELA was located slightly lower in the Western Massif (2252 m) than in the Central Massif (2341 m; González-Trueba, 2005; González-Trueba et al., 2008). Present-day Four ice-patches remains in the Western and Central massifs of Picos de Europa, nested within the LIA moraines but showing remarkably reduced surface extent compared to the former LIA glaciers (52 to 76% total surface reduction; González-Trueba et al., 2007). The most controverted one is the Jou Negro ice-patch, the single one that is not fully covered by debris. Although some authors described it as a glacier (GonzálezSuárez and Alonso, 1994; Alonso and González, 1998), most authors agree that it cannot be considered a 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 glacier because it lacks movement (Frochoso and Castañon, 1995; González-Trueba et al., 2005). The Forcadone buried ice-patch is another remnant of a LIA glacier with permanent negative temperature values at deeper layers (Ruiz-Fernández et al., 2017). 4.3 NW ranges The NW ranges comprise a series of mountain ranges that continue the relief of the Cantabrian Mountains towards the NW corner of the Iberian Peninsula (longitude between 6ºW and 8ºW), acquiring W-SW distribution trends. Bedrock geology mostly corresponds to Palaeozoic metamorphic and igneous rocks. Like in the Cantabrian Mountains, the distribution of Pleistocene glaciations in the NW ranges of the Iberian Peninsula was strongly conditioned by both the preglacial landscape topography and the proximity to the Atlantic Ocean. The NW ranges are generally arranged as smooth, near horizontal topographic highs set at increasing altitudes inland. Glacial evidence has been extensively documented since the 1910s in the highest massifs of Trevinca (2128 m), Serra dos Ancares (1998 m) and Manzaneda (1781 m) (Taboada, 1913), but also in other relatively low mountain settings like Serra do Courel (1654 m), Serra de Xurés-Gerêz (1548 m), Serra do Oribio (1484 m), Montes do Cebreiro (1474 m), and Serra do Xistral (1031 m) (Pérez-Alberti et al., 2004, 2018). Glaciations prior to the Last Glacial Cycle The idea of mountain glaciations occurring prior to the Last Glacial Cycle was initially proposed by Hernández-Pacheco (1949, 1957) to explain the glacial record existing in the Manzaneda massif. He ascribed the lowest moraines to the Riss glaciation – following the Alpine terminology – and estimated the altitude of the paleoELA in 1428 m. Similarly, Llopis (1957) interpreted the Sanabria Lake moraine complex preserved in the eastern slope of the Trevinca massif at ca. 1000 m as the result of three superimposed glaciations (Mindel, Riss and Würm). However, absolute ages based on 10Be CRE dating indicate that the Sanabria Lake moraine complex was deposited during the Last Glacial Cycle (Rodríguez-Rodríguez et al., 2014), showing time consistence with OSL analysis of glacio-fluvial sediments from the proglacial site of Pias in the opposite slope (Pérez-Alberti et al., 2011). Preliminary CRE dating carried out in the Manzaneda and Xurés-Gerêz massifs provided non-conclusive results to support the occurrence of glaciations before the Last Glacial Cycle (Table 4). Glacial evidence in the Manzaneda massif is indicative of a mountain ice cap covering a total surface area of 130 km2 (Vidal-Romaní and Santos-Fidalgo 1994) drained by glacier outlets up to 7-10 km in length and glacier thickness up to 150200 m (Pérez-Alberti et al., 1993). The 21Ne dating of a single boulder from the Castañeda moraine (1210 m) yielded an age older than 120 ka, whereas other three landforms sampled just 0.5 km upstream in the Cenza valley revealed 10Be and 21Ne ages ascribable to the Last Glacial Cycle (Vidal-Romaní et al., 1999, 2015). Similarly, two 21Ne analysis from two glacial polished surfaces sampled 4 km downslope in the Vilameá valley and at the mountain divide in the Serra de Xurés-Gerêz yielded surprisingly old age results of 238 ka and 130 ka for the maximum extent of an ice cap of ca. 64 km2 surface that was installed on the Xurés-Gerêz massif (Vidal-Romaní et al., 1999) and drained to the N by outlet glaciers up to 5-6 km long (Vidal-Romaní et al., 1990). If correct, these 21Ne ages would imply that glaciers did not grow in mountain environments below 1500 m altitude during the Last Glacial Cycle, which mismatches the chronological results obtained in the nearby Trevinca massif. Although preliminary numerical ages are still scarce to sustain that glaciations older than the last glacial Pleistocene cycle occurred, the presence of poorly preserved evidence like big erratic boulders distributed outside the limits of the lowest moraines in some moraine complexes (e.g. the Sanabria moraine complex) might probably represent the last evidence remaining from old glaciations (Figure 6). Table 4 Last Glacial Cycle 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 Most authors have interpreted the glacial record preserved in the NW ranges as formed during multiple stages of the Last Glacial Cycle (Pérez-Alberti et al., 2004). Glacial evidence suggests that mountain glaciers developed extensively in the highest massifs of the NW Iberian corner (e.g. Pérez-Alberti et al., 1992, 1993, 2002; Vidal-Romaní and Santos-Fidalgo, 1994), but also in relatively low mountainous areas (e.g. Valcárcel, 1998, Valcárcel et al., 2002a), being possible to distinguish up to four glacier typologies (Pérez-Alberti et al., 1993; Pérez-Alberti and Valcárcel, 1998): (i) mountain ice caps drained by outlet glacier; (ii) Alpine glaciers with coalescent tributaries; (iii) Alpine glaciers composed by a single ice tongue; and (iv) cirque glaciers displaying incipient ice tongue development (Figure 6). Up to date, the best studied glacial record in the NW ranges is that located in the Trevinca massif. Glacial evidence suggests that a mountain ice cap developed across an area of 475 km2 drained by several outlet glaciers radially disposed along the pre-existing fluvial valleys (Rodríguez-Rodríguez et al., 2014). Numerical reconstructions indicate that the Trevinca ice cap reached an ice thickness up to 200–300 m on the high plateau (Cowton et al., 2009; Rodríguez-Rodríguez et al., 2014), while the Cepedelo moraines suggest that the Bibei outlet glacier was up to 500 m thick, being by far the thickest one in the entire NW Iberian corner (PérezAlberti et al., 1993, 2002). The Bibei paleoglacier was also the longest outlet glacier (30 km; SW aspect), followed by the Tera ice tongue (24 km; SE aspect), which had their glacier fronts at minimum altitudes of 900 to 950 m, respectively (Pérez-Alberti et al., 1993; Rodríguez-Rodríguez et al., 2011). Glacio-fluvial sediments from the Pias proglacial site, deposited within the margins of the lowest front in the Bibei valley, have provided minimum ages in the range between 27 and 31 OSL ka for the local MIE of the Last Glacial Cycle (Pérez-Alberti et al., 2011). Similarly, radiocarbon analysis of pollen concentrates retrieved from the base of the Sanabria Lake sedimentary record also revealed a minimum age of 25.7–25.9 14C ka cal BP for the local MIE (Rodríguez-Rodríguez et al., 2011). Few chronological data are available for other massifs of the NW ranges, though it is likely to consider that they followed a similar temporal pattern than that documented in Trevinca. The preglacial topography in other mountain areas like Serra dos Ancares or Serra do Courel did not favour the formation of ice caps, and Alpinetype glaciers of different complexity developed instead (Pérez-Alberti et al., 1993; Pérez-Alberti and Valcárcel, 1998). Valley glaciers featured asymmetric distribution patterns between opposite mountain slopes, most probably influenced by a combination of controlling factors such as the preglacial landscape topography, variable topoclimatic conditions and the snowdrift effect of wind (Hernández-Pachecho, 1957; Pérez-Alberti et al., 1993; Vidal-Romaní et al., 1994; Valcárcel et al., 2009). In Serra dos Ancares, coalescent alpine glaciers up to 7-13 km long flowed along the SE slope whereas alpine glaciers did not overcome 10 km in length in the NW slope. However, the glacier fronts reached, in general, lower elevations in the NW slope (700–1000 m) compared to the SE one (825–920 m), most likely due to differences in the regional slope angle and relief between both mountain sides (Pérez-Alberti et al., 1992, 1993; Valcárcel et al., 2002b). Glaciers were slightly shorter in lower ranges like Serra do Courel (Pérez-Alberti, 2018), forming coalescent Alpine glaciers up to 5.5 to 7.8 km long along the Seara (150 m thick; front at 950 m) and Vilarbacú valleys (100 m thick; front at 900 m), while other valleys recorded simple Alpine glacier tongues significantly shorter (1.5 to 3 km; front at 650–950 m). The shortest glaciers of the NW ranges have been documented in Montes do Cebreiro and Serra do Oribio, where cirque glaciers with incipient glacier tongues up to 0.7–2.5 km long and less than 80 m thick developed, reaching minimum elevations of 900–1000 m (Pérez-Alberti et al., 1993; Valcárcel et al., 2002a). The fact that these two massifs do not reach 1500 m altitude, suggests that climate conditions during the Last Glacial Cycle in the NW ranges allowed sustaining glaciers at extremely low elevations possibly favoured by remarkable moisture supply from the Atlantic Ocean. Moreover, glacial evidence described in the coastal ranges of Serra do Xistral and Faro do Avión, distributed at elevations as low as 600–700 m indicating extremely low ELA conditions (ca. 900 m altitude) during the MIE of the Last Glacial Cycle (Schmizt, 1969; Pérez-Alberti et al., 1993). An extreme case is detected in the Capelada massif, where moraines formed during the MIE of the Last Glacial Cycle are located just above present-day sea level (Pérez-Alberti, 2014; Figure 6c). In contrast, the paleoELA rose progressively inland to 1250–1350 m in Serra do Courel, 1350–1450 m in 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 Serra dos Ancares and up to 1500–1600 m altitude in the Trevinca massif (Valcárcel et al., 2002b; PérezAlberti et al., 2004). Figure 6 Thus, the well-constrained glacial chronology of the Trevinca massif indicates that the MIE of Last Glacial Cycle probably took place prior to 26–31 ka. Glacier fronts remained stable until 22 ka, as suggests the oldest 10Be ages obtained in the outermost lateral moraine of the Sanabria Lake moraine complex (RodríguezRodríguez et al., 2014). Radiocarbon analysis of the lowest unit retrieved from the San Martin de Castañeda kame terrace – deposited outside this lateral moraine – also provided a minimum age of 21.8–22.1 14C ka cal BP for the impoundment of the lateral tributary due to moraine build-up (Rodríguez-Rodríguez et al., 2011). Non-glaciated mountain areas and lower coastal ranges about 600 m were affected by periglaciation, with sedimentary sequences that alternate colluvium deposits with organic-rich paleosoil intervals indicating periglacial activity between 33 and 48.7 14C ka cal BP (Butzer, 1967; Brosche, 1982; Costa-Casais et al., 1994, 1996; Cano et al., 1997; Costa-Casais, 2001; Blanco-Chao et al., 2007; Pérez-Alberti et al., 2009; Oliva et al., 2016). Therefore, the periglacial belt was close to the coastline during the Last Glacial Cycle, supporting the idea that the glacial record preserved in most mountain areas of the NW ranges corresponds to the Last Glacial Cycle. Termination-1 The moraine complex preserved around the Sanabria Lake (1000 m), in the Trevinca massif, is the only one that has provided detailed numerical age to constrain the timing of the onset of Termination-1 based on 10Be CRE dating of moraine boulders (Rodríguez-Rodríguez et al., 2014). Recalibrated results support that glacier recession started after 22 10Be ka and recorded multiple glacier front stagnations with construction of recessional moraines enclosing the lake until 18.5 10Be ka (Figure 6b). The basal age of lacustrine sequences like Laguna de las Sanguijuelas, deposited between two recessional moraines, provided basal radiocarbon ages of 17.8–18.2 14C ka cal BP (Muñoz-Sobrino et al., 2004) that are time consistent with moraine ages. Radiocarbon ages from the base of Laguna La Roya (15.3–15.6 14C ka cal BP; Allen et al., 1996), located at ca. 1630 m on top of the Trevinca high plateau, suggests that the lowest sectors of the Segundera highlands were deglaciated by the time of the OD. Similar conclusions are obtained from Lagoa Grande o das Lamas Lake (1364 m) in the Manzaneda massif, which basal age of 15.0–15.6 14C ka cal BP indicates that glaciers also disappeared from the lowest sector of this massif during the OD (Maldonado, 1994), consistently with the 21Ne age of 15.4 ka obtained in a drumlin of the Cenza valley (Vidal-Romaní et al., 1999); however, the existence of a drumlin in the area is highly controversial and is not supported by most of the scientific community. The radiocarbon date obtained from a pollen concentrate at the top of the Sanabria Lake basal unit, characterized by inorganic massive sand, suggests that the influence of ablation waters discharge from the glacier disappeared from the eastern sub-basin of the lake around 14.1–14.5 14C ka cal BP (RodríguezRodríguez et al., 2011; Jambrina-Enríquez et al., 2014). The basal age of the Lleguna sequence, deposited in a glacial over-deepening depression located to the South of the Sanabria Lake, also provided a minimum age of 13.8–14.1 14C ka cal BP for the time of deglaciation (Muñoz-Sobrino et al., 2004). Thus, numerical ages from the Sanabria Lake and Lleguna sequence suggests that the recession of the Tera glacier front accelerated at ca. 14.5 14C ka cal BP, at the beginning of the BO. A thin inorganic interval inter-bedded in the organic-rich sequence of Sanabria Lake record has been interpreted as a possible short-lived glacier readvance between 12.1 and 13.1 14C ka cal BP, before the YD (Jambrina-Enríquez et al., 2014). Most moraines preserved at higher elevations than the Sanabria moraine complex in the Tera valley (1600–1880 m) have not been dated yet, nor those preserved in other massifs (e.g. at 1200–1800 m in Ancares, 1200–1500 m in Courel and Manzaneda). In any case, glacier retreat was probably faster in mountain settings with top elevations under 1700 m altitude, as suggests the basal radiocarbon age obtained in Lucenza Lake (1374 m), in Serra do Courel (20.6–21.4 14C ka cal BP; Muñoz-Sobrino et al., 2001). 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 Present day Since 12.1 14C ka cal BP (Jambrina-Enríquez et al., 2014), there is neither geomorphic nor sedimentological evidence of recent glacial activity in the NW ranges. Currently, there is only evidence of active nivation processes related to long-lying snow patches in some mountain environments, with formation of protalus ramparts, intense abbrasion of rock surfaces – even forming striations – and mobilization of a large amount of material down-slope (Carrera-Gómez and Valcárcel, 2018). Based on lichenometric measurements, these authors also suggest the occurrence of different stages with more intense nivation processes during the LIA in some NW ranges. 4.4 Central Range The Central Range, running SW-NE between parallels 40–41ºN, divides in part the Iberian Peninsula approximately in the center. This mountain range is composed of a series of sierras, which are in fact tectonic blocks, formed mostly by crystalline Paleozoic and Precambrian rocks, arranged longitudinally along the axis of the range and separated by large fractures and grabens. At present, there are no active glaciers in the Central Range, but most of its sierras conserve glacial landforms on their summits, which become more significant depending on how far west they are located, and therefore how close they are to the influence of the humid Atlantic air masses. At the western end of the range, the Serra da Estrela (Pico Torre, 1993 m) has glacial valleys more than 13 km long. At the central part, the Sierra de Gredos (Pico Almanzor, 2591 m) and Sierra de Béjar (Pico Calvitero, 2399 m) include also glacial valleys more than 14 km long. The Sierra de Guadarrama (Pico Peñalara, 2428 m) and Sierra de Ayllón (Pico Lobo, 2274 m) at the eastern end of the range, have only glacial cirques less than 2 km long. Glaciations prior to the Last Glacial Cycle In the Central Range, the existence of glacial landforms predating the Last Glacial Cycle was proposed for the Laguna Cirque, Sierra de Guadarrama, under the eastern face of Pico Peñalara (Table 5), where there are small maximum advance moraine ridges formed in the Riss glaciation, following the Alpine terminology (Penck, 1883). This relative chronology was accepted later by several other researchers (Fernández Navarro, 1915; Obermaier and Carandell 1917, and others) until the mid-20th century, when these moraines were assigned to the Würm, although no dates were provided (Butzer and Fränzle, 1959). In the early 21st century, CRE dating methods were applied to these moraines (Palacios et al., 2012a), obtaining results that show that they formed during the last glacial cycle, synchronous with or slightly earlier than the LGM. This chronology was confirmed in subsequent studies also in the Sierra de Guadarrama that also used CRE dates, both in the moraines of the Laguna Cirque (Domínguez-Villar et al., 2013) and in the Pelados-El Nevero massif (Carrasco et al., 2016). By contrast, in the Sierra de Gredos and nearby mountains similar evidence of glacial landforms that may have originated during previous glaciations has never been found (Schmieder, 1915; Huget del Villar, 1915, 1917; Obermaier and Carandell, 1916; Vidal Box, 1936; Martínez de Pisón and Muñoz, 1972; Pedraza and Fernández, 1981), and this has been confirmed by more recent studies (Palacios et al., 2011, 2012b; Pedraza et al., 2013; Carrasco et al., 2013, 2015; Domínguez-Villar et al., 2013). The same also occurs in the Serra da Estrela, from the earliest studies of relative chronology (Lautensach, 1929; 1932; Daveau, 1971; Daveau et al., 1997) to the most recent using TL absolute dating techniques (Vieira et al., 2001; Vieira, 2004, 2008); however, ongoing studies have revealed MIS 6 ages for the most external lateral moraines using CRE techniques (Vieira, Palacios and Lorenzo, unpublished data). Table 5 Last Glacial Cycle Our current knowledge of the glacial evolution of the Central Range limits the chronology of the glacial landforms to the last Quaternary glacial cycle, and specifically, to its last phases, with the MIE occurring 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 slightly prior to the LGM and important advance and retreat stages during Termination I, but with no geomorphic evidence of Neoglacial ice expansion. As mentioned above, the distribution and typology of glacial landforms in the Central Range change significantly from the W to the E. In the eastern Sierras de Ayllón and Guadarrama only small cirques are preserved, mostly to leeward of winds bringing the westerly storms (Hernández-Pacheco, 1930) and in most cases there is only a single moraine ridge, or at most two. These small glaciers formed under summits at altitudes higher than 2100 m and their fronts reach elevations of ca. 1800 m. The paleoELA for the Sierra de Guadarrama was estimated at 1900-2000 m (Sanz-Herráiz, 1988). The largest of these cirques is the Laguna cirque in the Peñalara massif, which is approximately 1.8 km long. In this cirque, small moraine ridges – which were considered older than the last glacial cycle (Fernández Navarro, 1915; Obermaier and Carandell 1917) – are located in front of a large complex frontal moraine. CRE dating using 36Cl revealed the ages of these small moraines of the maximum advance ranging from 29–31 ka for a boulder on a northern side of the ridge to 16–22 ka for various boulders on a southern slope (Palacios et al., 2012a). In parallel, the main moraine system composed of several ridges superimposed on the front and aggregated into only one on the laterals yielded ages ranging from 16 to 27 ka. Consequently, these data suggest that the local MIE in this massif preceded the LGM by several thousand years, when the front was relatively stable, with subsequent minor advances and retreats resulting in a large polygenic moraine. The magnitude of this moraine was such that large snow patches developed on its fronts to leeward of the westerly winds, displacing boulders which then formed what seem to be small maximum advance ridges. This process may explain why the same ages are found in the main moraine as well as in the small ridges outside the glaciated environment (Palacios et al., 2012a). In the same Peñalara massif, CRE dating also confirmed a LGM age of 26 36Cl ka for the main moraine of the Pepe Hernando cirque (Palacios et al., 2012a). The chronological framework for Peñalara massif coincides with CRE 10Be dates in other small Guadarrama cirques, such as in the Los Pelados-El Nevero massif, 16 km NE of the Peñalara massif (Domínguez-Villar et al., 2013; Carrasco et al., 2016). Other glacial landforms existing in Guadarrama cirques are possibly synchronous to those of Peñalara massif based on their similar morphology and degree of preservation, though no chronological data are available to date (Fernández Navarro, 1915; Obermaier and Carandell, 1917, Fränzle, 1959; Sanz-Herráiz, 1988; among others). A block field located on the Peñalara summit has been CRE dated at 80 36Cl ka, which confirms that in most of the eastern sierras of the Central Range the summits were ice-free during the local MIE (Palacios et al., 2012a). A small plateau glacier with some ice-free summits has only been described in the Pelados-el Nevero massif, a very small area in these mountains. The age of these nunataks protruding the glacial ice ranged from 33.5 to 59.7 10Be ka, which is slightly younger than those obtained on the summit of Peñalara, but they confirm that these surfaces were ice-free during the MIE (Carrasco et al., 2016). On the north face of the Pico Mujer Muerta, also in the Sierra de Guadarrama, some sedimentary formations have been classified as alluvial fans formed during the paraglacial stage and dated at 30-35 OSL ka (Bullón, 2016). Contrary to what occurs in the eastern sierras of the Central Range, the central Sierras de Gredos and Béjar were covered by ice caps with a larger extension to the N and E, and with a few short tongues descending to the S (Pedraza et al., 2013, Carrasco et al., 2013). These authors synthesized the geomorphological evolution that is common to most glaciated valleys in these central sierras of the Central Range, where a large moraine was formed. But in many cases, in sectors further from the sides or front of the main moraine, there are isolated glacial boulders or small moraine ridges (Figure 7a). These deposits were dated by CRE methods using 10Be isotope, and as in Peñalara, the ages obtained resulted in a glacial stage that occurred several thousand years before the LGM (26–31 ka) (Domínguez-Villar et al., 2013, Carrasco et al., 2013, 2015). This period of maximum ice extent would have coincided with a wet, cold phase, as shown by the U-Th speleotheme dating in a cave in the southern foothills of the Sierra de Gredos (Domínguez-Villar et al., 2013). 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 Figure 7 There have been several atempts to constrain the age of the main moraine in some valleys in Gredos and Béjar, which must have formed during several stages of advance and stabilization. The largest glaciers (ca. 14 km long) formed in the highest northern valleys of the Sierra de Gredos and reached elevations down to 1400 m, with the ELA estimated at 1830 m (Palacios et al., 2011). In the Gredos valley, under the north face of the Pico Almanzor, some boulders of the main moraine system were dated in areas just where the valley widens and the slope gradient decreases. This moraine, which was a single lateral polygenic ridge, divides in this sector to form several very stable ridges. CRE dating using 36Cl isotope reported ages of 20–26 ka for these boulders, suggesting a LGM age (Palacios et al., 2011). In the nearby Pinar valley, the same dating method applied to several boulders of the main moraine system showed very similar results, with ages ranging from 19 to 24 ka (Palacios et al., 2012b). The Sierra de Béjar was covered by an ice cap of 57 km2, with the glacier tongues (max. length 8 km) descending to around 1200 m, and a mean ELA of 2010 m (Carrasco et al., 2013). Boulders from moraine formations in the Duque, Trampal, Endrinal and Cuerpo de Hombre valleys dated by CRE using isotope 10Be yielded ages from 19 to 27 ka (Carrasco et al., 2013, 2015), very similar to those of Gredos. The largest ice cap of the Central Range was located in its western fringe, the Serra da Estrela, encompassing a surface of 66 km2 and a thickness of 340 m in some valleys. At that time, the ELA was placed at 1640 m (Vieira, 2004, 2008). Only few data on glacial landforms are available; some fluvio-glacial deposits were dated through OSL at 30–35 ka (Vieira et al., 2001) and a preliminary approach using CRE dating by the 36Cl isotope on boulders of the large moraine complexes yielded ages around 21–26 ka, very similar to the rest of the Central Range (Palacios et al., 2012c). In the Serra da Estrela, the largest moraines were dated at 31–26 36Cl ka, with a sequence of internal moraine ridges indicative of several minor advances and retreats until 18–16 36Cl ka (Vieira and Palacios, 2010; Vieira, Palacios and Lorenzo, unpublished data). In summary, the available data suggest that the maximum ice extent occurred slightly before the LGM, at 27– 35 ka, although it was a short glacial advance that left only minor landforms. The LGM advance was very significant and lasted for thousands of years, with small advances and retreats, which led to the development of a large polygenic moraine in each valley. Termination-1 In the Peñalara massif, moraines that may have formed during the OD are spatially close to those of the LGM, with ages of 16–17 36Cl ka. Some of these areas were affected by intense paraglacial dynamics that triggered the development of rock glaciers, though they stabilized shortly after their formation at around 15 36Cl ka (Palacios et al., 2012a). Something similar occurred in the Gredos valley where polished surfaces located behind the main moraines have been dated at around 16 ka 36Cl. In the Pinar valley and in Peñalara, moraines of 16–17 36Cl ka are located next to those of the LGM features. Similarly, in the Cuerpo de Hombre valley, moraines dated at 15–17 10Be ka are very close to the LGM. This evidence highlights the occurrence of important glacial advances in the Central Range during the OD, with slightly smaller glaciers than during the LGM and even, in come cases, next to them (Palacios et al., 2017). In the Serra da Estrela, a polished threshold near the summit shows that the ice cap disappeared at 15 36Cl ka (Vieira and Palacios, 2010; Palacios et al., 2012c). Only recently, some data are indicative of small glacial advances during the YD in the Cuerpo de Hombre valley, with moraines dated at 11–13 10Be ka (Carrasco et al., 2015). In conclusion, the deglaciation of the Central Range seems to have followed a similar pattern to that detected in other European mountains, with a major advance during the OD that favoured the expansion of glaciers 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 until almost the LGM moraines or until a few hundred meters from them (Palacios et al., 2017). The recent dating of a moraine of the YD in the Sierra de Béjar may suggest that the YD also promoted glacial expansion in the Central Range, though few dates are already available for this period in contrast to in other Iberian mountains where YD moraines are widespread (García-Ruiz et al., 2016). Holocene Three CRE datings reveal the definitive disappearance of glaciers in the Central Range during the Early Holocene. A polished bedrock surface below the summit of Peñalara indicates an approximate age of 11 36Cl ka (Palacios et al., 2012a). At the head of the Pinar valley, the dating of a abraded surface also yielded an age between 10 and 11 36Cl ka (Palacios et al., 2012b). In the Cuerpo de Hombre valley, the moraine located under the headwall yielded a minimum age of 11 10Be ka (Carrasco et al., 2015). Although there is still much to research, the available data seem to indicate that the glaciers of the Central Range disappeared definitively in the transition between the Pleistocene and the Holocene, or in the onset of the Holocene. LIA There is no geomorphic evidence of subsequent formation of glaciers, neither during the Neoglacial period nor during the LIA. There is only one protalus rampart located in the Gredos cirque that was dated by lichenometry and showed its development during the LIA (García-Sancho et al., 2001), although there might be other similar landforms in other Central Range cirques. Present day Current cold-climate geomorphological processes in the Central Range are only limited to the erosive activity of long-term snow patches, particularly during colder periods, such as the late 1960s and early 1970s (Palacios et al., 2003). 4.5 Iberian Range The Iberian Range limits the Central Iberian plateau by its N and NE sides. This range, extended in NW-SE disposition along 500 km, is composed both by basement fragments and sectors of sedimentary cover with variable thickness. From N to S and W to E the main ranges are: Sierra de la Demanda (San Lorenzo, 2271 m), Sierra de Neila (Campiña, 2049 m), Picos de Urbión (Urbión, 2228 m), Sierra Cebollera (La Mesa, 2168 m) and Sierra del Moncayo (Moncayo Peak, 2314 m). The first approach to the Iberian Range glaciers corresponded to Carandell and Gómez de Llarena (1918), mainly focused on the Sierra de Urbión, which was also studied in detail by Thornes (1968). The Sierra de la Demanda was analyzed by García-Ruiz (1979), Antón-Burgos (1985), Arnáez-Vadillo (1987), Arnáez-Vadillo and García-Ruiz (1990) and recently, Fernández-Fernández et al. (2017), whereas the Sierra Cebollera was studied by Ortigosa (1985, 1986), Sanz Pérez and Pellicer (1994) and García-Ruiz et al. (2007), who published a geomorphologic map of the Sierra Cebollera. A brief description of the geomorphology of the Sierra de Neila was performed by Ortega and Centeno (1987). The Sierra del Moncayo has been the object of detailed studies by Martínez de Pisón and Arenillas (1976) and Pellicer (1984). Besides, García-Ruiz et al. (1998) synthesized the available information on glacial landforms and deposits in the Iberian Range. Table 6 Glacier snouts reached 1650 m during the MIE in the Sierra de la Demanda, similarly to the Sierra del Moncayo, where glaciers remained confined in three cirques of steep slopes. By contrast, glaciers in the Neila, Urbión and Cebollera massifs were hosted in larger cirques favoured by the structural setting shaped by the main divide formed by a cuesta scarp with low-dip strata (Ortigosa, 1986). The best examples of glacial valleys in the Iberian Range are found in these ranges, with well-developed glacial troughs that can reach 6-km long 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 in the Urbión valley, and over 3-km long in the NE cirque of La Mesa de Cebollera Peak. The Urbión valley encompasses the glacial deposits at elevations of ca. 1270 m, and even lateral obturation deposits at the confluence with tributary ravines (García-Ruiz et al., 1998). CRE dating techniques, through 10Be, have only been applied in the Sierra de la Demanda, in the westernmost cirque of the Mencilla Peak and in the SE cirque of the San Lorenzo peak (Fernández-Fernández et al., 2017). Figure 8 Glaciations prior to the Last Glacial Cycle Lotze (1962) attributed the glacial morphologies existing in the Sierra de la Demanda to the Riss glaciation but did not provide any absolute date. Nonetheless, the hypothesis of glacial landforms being originated in a previous glaciation in Sierra de la Demanda was rejected by García-Ruiz (1979) in favor of the Würm glaciation as the only glacial stage including different climatic oscillations. The author defended this origin based on their fresh morphology, with scarcely developed landforms and absence of postglacial regressive erosion, together with the similar grain size and composition of the moraines as well as the short distance between them. The best example of this fact was found in the westernmost cirque of the Mencilla Peak (GarcíaRuiz, 1979), although the author did not dismiss the possibility that the outermost moraine in the SE San Lorenzo cirque could have formed during a previous glaciation as suggested by the existence of undifferentiated till distributed across the slope as well as the presence of a soil layer on the surface. Based on the similar gran size of the frontal moraines, Thornes (1968) also attributed the glacial features existing in Urbión to the Würm glaciation. In fact, the existence of older moraines in the SE cirque of the San Lorenzo peak cannot be discarded due to the presence of erratic boulders outside the outermost moraine (FernándezFernández et al., 2017). Last Glacial Cycle The outermost moraine of the Mencilla cirque, Sierra de la Demanda, could not be dated due to the fact that intense weathering has resulted in the inexistence of suitable boulders for CRE dating. In the case of the SE San Lorenzo cirque, the outermost moraine yielded a minimum age of 18.1 ± 2.3 10Be ka (Table 6). Although the intermediate moraine could not be dated, considering the Late Pleistocene chronology of neighbouring glaciated areas, the relative position of the moraines within the cirque and the partial overlapping of the age of the outermost and the innermost (16.8 ± 1.4 10Be ka) moraines, it could be argued that the intermediate moraine could have been formed during the LGM (Figure 8). Thus, the outermost moraine might have been deposited during a previous glacial expansion (Fernández-Fernández et al., 2017). Glaciers hardly reached 1 km long during their maximum extent, and their ELAs descended to 1480 m in the Mencilla peak cirques and to 1671 m in the SE cirque of the San Lorenzo peak. The different massifs of the Iberian Range include abundant examples of moraine systems, sometimes up to three frontal moraines suggesting different stages of advance and retreat. Some of these moraines can be of great size, with large boulder accumulations resulting from the intense network of fractures affecting the cirque headwalls. Apart from those existing in the Sierra de la Demanda, there are also several good moraine complexes, such as in the Laguna Negra valley (Sierra de Urbión), in the NE valley of La Mesa de Cebollera (Sierra de Cebollera) and in the Morca and San Miguel cirques in the Sierra del Moncayo (Martínez de Pisón and Arenillas Parra, 1976; Pellicer, 1984). A well-developed moraine arc has been recently identified close to Orihuela del Tremedal, Albarracín Massif, in the SE limit of the Iberian Range, dominated by a quartzitic crest at 1920 m. This deposit has not been yet dated, although it most likely represents the maximum extent of a small, isolated glacier in a relatively low and southern massif (González-Sampériz et al., 2006-2007). Termination-1 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 fragmentary. For example, in the Italian Appennines evidence for glaciations have been identified in a proglacial lake basin dating to MIS 14, 12, 10, 8 and 6, in addition to moraines dating to the Last Glacial Cycle (Giraudi and Giaccio, 2017). 5.2 Chronology of the maximum ice extent of the Last Glacial Cycle The Last Glacial Cycle corresponds to the period MIS 5d-2 following the Eemian interglacial (MIS 5e) that ended at 115 ka (Dahl-Jensen et al., 2013). It was a period with highly variable temperatures alternating extreme cold stages (stadials; MIS-2, 4) with warmer stages with temperatures almost as high as present-day (interstadial; MIS-3). As a result of the combined effect of decreased northern summer insolation, lower tropical Pacific sea surface temperatures as well as low atmospheric CO2, ice accumulated in large ice sheets of the Northern Hemisphere between 33.0 and 26.5 ka, reaching their maximum position between 26.5 and 19 ka (Clark et al., 2009; Hughes et al., 2016a). This time period constitutes the so-called LGM, though some authors suggest a slightly older age of 27.5–23.3 ka (Hughes and Gibbard, 2015) based on the global dust record in polar ice cores. This is the coldest stage of MIS-2 and corresponds to the global LGM, when the greatest ice volume was stored in land masses, and thus sea level was lower. In the case of Iberian mountains, radiocarbon dating was the most widespread technique to determine the age of glacial activity during the Last Glacial Cycle until the implementation of other techniques, such as OSL and CRE dating. In the late 1990 and early 2000s several authors proposed a MIE of the Last Glacial Cycle occurring during MIS-4. These studies were based on OSL dating on fluvio-glacial sediments located near the outermost moraine complexes and suggested that the largest glacial advance had occurred between 50 and 70 ka. This is the case of several valleys in the Central Pyrenees (García-Ruiz et al., 2013) and is time consistent with the coldest global temperatures of MIS-4 recorded at 65–70 ka (Kindler et al., 2014). Some authors also suggested that older dates than those reported by 14C dating could reveal periods of glacial expansion at that stage, but the chronological limitations of radiocarbon dating impeded confirming this fact. Consequently, radiocarbon dating constrained the age of the MIE in most ranges at ca. 40–50 14C cal ka BP or older than this age, but without providing approximate dates. It must be taken into account that 14C requires the dating of organic fragments, but biological particles near glaciated environments must have been scarce during the deposition of the glacial features where these particles were trapped, namely moraines and till, due to the very cold climate regime; in addition, the long time passed since their deposition and often the poor degree of preservation of these deposits makes it difficult to find appropriate samples to be dated. Despite these limitations, in some mountain environments, such as in most NW massifs, this has been an extensively used dating approach; here, the available data indicate that very intense periglacial conditions probably associated to the MIE occurred within the range 33–48 14C cal ka BP, or even older (Table 4). Similarly, in the Cantabrian Mountains 14C dates also suggest a MIE of the Last Glacial Cycle occurring between 35 and 45 14C cal ka BP (Figure 12). Figure 12 The accuracy on the reconstruction of glacial oscillations has improved substantially the last decade with the use of CRE dating. This technique allowed inferring the age directly from glacial records (moraines, erratic boulders and polished bedrock) and expanding the number of sites from which chronological data are available. Consequently, we can infer the sequence of environmental events occurred in the Iberian massifs since the Last Glacial Cycle, which has major implications for paleoclimate reconstructions. In addition, the use of CRE allows the validation of 14C and OSL dates but has also introduced new uncertainties on the chronology of past glaciations. In most areas where CRE dating has been applied, age results suggest a MIE occurring during MIS-2 but show also contrasting temporal and spatial patterns. In some areas, such as in some valleys of the Eastern Pyrenees (Delmas et al., 2008; Palacios et al., 2015a; Andrés et al, 2018) and Central Range (Palacios et al., 2012a; Domínguez-Villar et al., 2013; Carrasco et al., 2016), CRE yielded ages for the MIE almost synchronous to the LGM (Figure 10). In other massifs, glaciers reached their maximum extent several 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 millennia before the LGM, at 30–35 ka, such as in the Sierra Nevada at ca. 30 36Cl ka (Gómez-Ortiz et al., 2012a; Palacios et al., 2016a) and in the Trevinca massif, NW ranges, at 33 10Be ka (Rodríguez-Rodríguez et al., 2011). In other mid-latitude mountain environments of the European continent, the LGM promoted the most extensive glaciers of the Last Glacial Cycle, such as in the Alps (Ivy-Ochs et al., 2008; Scapozza et al, 2015; Federici et al., 2016; Wirsig et al., 2016), Apennine (Giraudi, 2011, 2015), Tatra Mountains (Makos et al., 2014), Krkonoše Mountains (Engel et al., 2014) or Balkan mountain ranges (Hughes et al., 2010; 2011; Kuhlemann et al., 2013). The LGM coincided with a minima in solar radiation in the northern hemisphere at 24 ka (Alley et al., 2002). After this, during a period of increasing northern hemisphere insolation a massive and rapid deglaciation was underway at 19-20 ka both in the large northern ice sheets (Hughes et al., 2016a) and in midlatitude mountain environments due to changing orbital forcing patterns (Alley et al., 2002; Clark et al., 2009; Shakun et al., 2015). 5.3 Glacial advances and retreats during Termination-1 Orbital forcing variations induced changes on the northern summer insolation that resulted in the onset of Termination-1 at high latitudes and mountain regions at 19–20 ka, as well as an abrupt rise in sea level (Clark et al., 2009; Shakun et al., 2015). Changing orbital patterns together with large-scale events in the North Atlantic region induced a number of centennial to millennial-scale climatic fluctuations during the long-term deglaciation that favoured glacial expansion (Clark et al., 2012). For example, the OD occurred close in time to a major phase of iceberg rafting during Heinrich Event I (Palacios et al., 2016a) whilst the YD was associated with major meltwater discharge into the North Atlantic from the Laurentide ice sheet (Renssen et al., 2018). Greenland ice cores and Laurentine initial retreating point to a temperature increase at ca. 19 ka following the LGM (Lambeck, et al, 2014), which induced a major deglaciation of the large ice sheets of the Northern Hemisphere as well as of mountain glaciers that dramatically reduced their dimensions by 17–18 ka (Marks, 2015; Vasskog et al., 2015; Stroeven et al., 2016; Wirsig et al., 2016). The deglaciation was a rapid process in most Iberian mountains and in only few millennia glaciers were probably confined within the high lands of the highest massifs. This pattern also occurred in other European mountains, such as in the Alps, where by 18 ka the Rhine Valley had lost 80% of its ice mass (Ivy-Ochs et al., 2004). Subsequently, at ca. 17 ka there was a period of significant cooling across the Northern Hemisphere, the OD, with severe winters and mild summers that resulted in glacial expansion (Denton et al., 2005; Williams et al., 2012). In the Iberian mountains glaciers expanded again across the valley floors and deposited moraines at ca. 16.8–16.5 ka relatively close to the LGM advance, with a second readvance at 15.5 ka followed by a major retreat (Figure 13): 1 ka later glaciers were again confined within the limits of the cirques (Palacios et al., 2017a). This is the case of the Central (Palacios et al., 2015b; 2016b) and Eastern Pyrenees (Palacios et al., 2015a), Central Range (Palacios et al., 2012c, 2017b) and Sierra Nevada (Palacios et al., 2016a). Evidence of the OD is widespread in the Alps between 17 and 16 ka (Ivy-Ochs, 2015; Wirsig et al., 2016) as well as in several other mountain ranges of the Mediterranean region, where glaciers expanded significantly between 17 and 16 ka and formed moraines at only a few km from LGM moraines, such as in the Tatra Mountains (Makos et al., 2013, 2014), Carpathian Mountains (Reuther et al., 2007) and mountains of Anatolia (Sarikaya et al., 2008; Zahno et al., 2010; Akçar et al., 2014). Figure 13 However, the OD was followed by the BO, a period of abrupt warming of ca. 9 ºC inferred from Greenland ice cores (Clark et al., 2012; Buizert et al., 2014), slightly less in western Europe of 3–5 ºC (Clark et al., 2012). This warming was detected in the entire North Atlantic region and led to a major shrinking of Northern Hemisphere ice sheets (Briner at al., 2014; Marks, 2015; Stroeven et al., 2016) as well as to the disappearance 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 of glaciers from Alpine environments or their spatial confinement in the highest cirques from the highest massifs. In the case of Iberia, during the BA glaciers probably melted away completely or only persisted as very small features in the Central Pyrenees (Palacios et al., 2017a). Outside the Alps, where glaciers persisted during the OD due to their much higher altitudes (Ivy-Ochs, 2015), ice volumes shrunk in most Mediterranean massifs (Palacios et al., 2017a). Greenland ice cores suggest that the YD was ca. 4.5° ± 2°C warmer than the OD (Buizert et al., 2014), with temperatures in western Europe significantly lower (5–10 ºC) than those prevailing during the BA (Clark et al., 2012). The cooling across Europe was more remarkable in northern (4–5 ºC) than in southern latitudes (2–3 ºC) (Moreno et al, 2014; Heiri et al., 2015). This cold stage has been reported in many Iberian records associated with cold and dry conditions that resulted in hydrological changes in lakes and rivers with the development of fluvial terraces, reduction of speleothems growth, soil development in high cirques, forest decline, etc. (Fletcher et al. 2010; Moreno, 2014; García-Ruiz et al., 2014; Oliva et al., 2014a). Some terrestrial records in the Iberian Peninsula suggest a temperature decrease of 2.5 ºC with respect to the BA period (IriarteChiapusso et al, 2017), a trend that was even higher than 4–5 ºC in sea surface temperatures in SE Iberia (Cacho et al., 2002) and Portuguese coast (Rodrigues et al., 2010), with values similar or lower than during the LGM. Despite the prevailing aridity in the Iberian Peninsula (Naughton et al., 2015), the intense cold allowed the development of small glaciers in the Pyrenees (García-Ruiz et al., 2014), Central Range (Carrasco et al., 2015) and Sierra Nevada (Gómez-Ortiz et al., 2012a; Palacios et al., 2016a), and probably in the Cantabrian Mountains (Figure 14). Ice tongues during this stage were short and exceeded only the limits of the glacial cirques in a few km. YD glaciers are recorded elsewhere in European mountains, such as the Alps where two moraine systems developed between 13.5 and 12 ka (Ivy-Ochs, 2015) or Tatra Mountains where moraines date to 12.5 ka (Makos et al., 2013). In other mountains of the Mediterranean region moraines of similar ages to those dated in the Iberian Peninsula have been reported; this is the case in the Italian Apennines where radiocarbon ages from an ice-marginal lake basin suggests a YD age (Giraudi and Frezzotti, 1992; Giraudi, 2015), the Moroccan High Atlas where moraines have yielded an average CRE age of 12.2 ± 1.6 ka (Hughes et al., 2011) and Anatolia where moraines date to between 11.5 ± 0.8 and 13.0 ± 0.8 10Be ka in the Kaçkar Mountain range (Akçar et al., 2007) and around 11.5 ± 1.0 and 13.3 ± 1.1 10Be ka in the Uludağ Mountain (Zahno et al., 2010). Whilst glaciers were present in many of the mountains of Iberia and the wider Mediterranean during the YD, it is unclear whether the YD was a phase of significant glacier expansion or simply a phase of glacier stabilisation in a period of overall glacier retreat during Termination I. Evidence from further north along the NE Atlantic margin in the British Isles now suggests that the YD glaciers may have been smaller than those during the preceding BA (Bromley et al., 2018) indicating a period of overall glacier recession from the OD through the BA to the YD. As noted above, in Iberia the climate conditions were less favourable for glaciers in the YD compared with earlier OD and glaciers are likely to have stabilised as a result of the large depression in temperatures between the BA and YD despite the latter being a very dry. How widely this climatic scenario can be applied is open to debate and further research is needed to understand the glacier fluctuations across the Mediterranean mountains close to Termination I. Figure 14 5.4 Holocene glacial stages in the Iberian highest ranges Following the YD, temperatures increased substantially on the order of ca. 10 ºC in Greenland and 4 ºC in western Europe (Clark et al., 2012), which led to a migration of cold-climate processes to higher elevations in the entire North Atlantic region. The Holocene Thermal Optimum was a period of higher temperatures associated with the local orbitally forced summer insolation maximum in the Northern Hemisphere lasting from 11 to 5 ka (Renssen et al., 2018). Consequently, glaciers developed during the YD gradually shrunk and many finally disappeared during the Early Holocene, such as the large Pleistocene ice sheets existing in North America (Stokes et al., 2017) and Fennoscandia (Stroeven et al., 2016). In the case of Iberian glaciers, they only persisted during the Early Holocene as features of reduced dimensions in the highest mountain ranges, 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 namely in the Central Pyrenees (Palacios et al., 2017a; Tomkins et al., 2018), Central Range (Palacios et al., 2012a, 201b; Carrasco et al., 2015) and Sierra Nevada (Palacios et al., 2016a). As temperatures kept rising, their disappearance favoured paraglacial activity and the development of permafrost-related features, such as rock glaciers, which became inactive soon after their formation (Oliva et al., 2016a; Palacios et al., 2016b; Andrés et al., 2018). The temperature increase recorded during the Early Holocene in other Alpine regions of the European continent also led to significant glacial retreat, such as in the Alps (Ivy-Ochs et al., 2009; Scapozza et al, 2015; Federici et al., 2016). Climate variability increased in the North Atlantic region during the Mid and Late Holocene, with temperature oscillations on the order of ± 2 ºC (Mayewski et al., 2004). However, there is little evidence of glacial activity during the Mid and Late Holocene in the Iberian mountains, with glacier occurrence only identified in the Central Pyrenees and Sierra Nevada. It is likely that some of the moraines systems existing in other massifs of the Central Pyrenees and Cantabrian Mountains formed during Holocene glacial advances, but no absolute dates confirm that fact yet. García-Ruiz et al. (2014) dated a glacial advance at 5.1 ± 0.1 36Cl ka in the Marboré Cirque as well as another Neoglacial expansion that was followed by a period of glacial shrinking between ca. 3.4 and 2.5 36Cl ka and another glacial stage at 1.4–1.2 36Cl ka. Oliva and Gomez Ortiz (2012) found evidence of the development of glaciers in the Sierra Nevada during Bond events at 2-8–2.7 and 1.4–1.2 14C cal ka BP. In both cases, the last stage of glacier expansion during the Holocene recorded in both cirques occurred during the Dark Ages and the LIA. These two last glacial advances have been also recorded in many mountain regions and polar environments (Solomina et al., 2016). Yet, as highlighted by Matthews (2013), the chronology of Neoglacial glacial advances in European mountains is still not properly understood and needs to be reassessed. However, in some cases, there are no moraine systems between the YD and LIA moraines, which may suggest that the terminal moraines correspond to polygenic systems that make impossible to differentiate previous Holocene glacial stages from LIA moraines (Crest et al., 2017), such as in the Alps and Norway (Matthews, 2013) and other mountain regions across the Northern Hemisphere (Solomina et al., 2015). 5.5 The LIA, background of present-day glacial retreat Climate variability during the Late Holocene includes several warm (e.g. Medieval Climate Anomaly, recent warming) and cold stages (e.g. Dark Ages, LIA). Precisely, the LIA has been defined as one of the coldest stages in the Northern Hemisphere of the last 10000 years (Bradley and Jones, 1992). In the Iberian Peninsula, a recent multi-proxy reconstruction for that period indicates that the onset of colder conditions following the medieval warm epoch started around 1300 and ended by 1850 (Oliva et al., 2018). The progressive decline of temperatures during the early 14th century in the Iberian Peninsula favoured the gradual development of glaciers at the foot of the highest peaks, namely in the Central Pyrenees, Picos de Europa and Sierra Nevada. The combination of precipitation and temperatures conditioned their dimensions, enlarging during cold and wet stages and shrinking (and possibly disappearing) during warm and dry phases within the LIA. During the LIA, the Sierra Nevada and Picos de Europa only encompassed 2 and 6 glaciers, respectively, whereas 111 glaciers existed in 15 different massifs in the Central Pyrenees. In all cases, glaciers were of reduced dimensions (several ha) – Aneto glacier was the largest, with a glaciated surface of 236 ha (GonzálezTrueba et al., 2005), and 610 ha for the entire Aneto-Maladeta massif (Rico et al., 2017) –, north exposed at the foot of steep cirque walls protecting ice and snow melting and favouring its accumulation by wind effect. Glaciers in the Sierra Nevada and Picos de Europa massifs remained confined within the cirques, whereas in the Pyrenees there were some small alpine glaciers flowing down-valleys as well as some glaciers in southern cirques of the highest massifs of this mountain range. The elevation difference between the ELA and the highest summits as well as the steepness of the cirques conditioned the amount of ice stored in these mountains. The ELA in the Picos de Europa was located at 2250–2340 m, in the Pyrenees was placed at 2620–2945 m, and in the Sierra Nevada around 3000–3100 m (Oliva et al., 2018). However, little is known about when these LIA glaciers started to form, except for Sierra Nevada. There is a well-established age of formation of the ice mass in the Mulhacén cirque, which was reconstructed based on 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 the sediments of La Mosca Lake that revealed its appearance by 1440 and disappearance by 1710, at the end of the Maunder Minimum cold stage (1645-1715). In the nearby Veleta cirque, a glacier formed during the early 14th century and readvanced during the 17th century (Palacios et al., 2019). During this time, solar radiation minimum conditioned MAAT ca. 2 °C below present-day values (Oliva et al., 2018), with glaciers reaching their maximum volume of the LIA in the Sierra Nevada and the Pyrenees. A significant expansion was also recorded between 1805 and 1830 following the Dalton Minimum (1790–1820) (García-Ruiz et al., 2014). A similar timing was also recorded in the Alps, when glaciers expanded during these stages (Zemp et al., 2015; Zumbuhl et al., 2018), as well as in several other mid-latitude mountain regions across the world (Solomina et al., 2016). In other Mediterranean massifs glaciers generally reached their maximum advance of the LIA – and, in some cases, of the Holocene – during the 17th century with a significant expansion also during the mid-19th century (Hughes, 2014). 5.6 The Pyrenean, last Iberian current glaciers MAAT have increased by approximately 1 ºC since the end of the LIA by 1850, which has conditioned a shift of cold geomorphological processes to higher elevations, with the complete deglaciation of the Picos de Europa and Sierra Nevada as well as of most massifs in the Pyrenees (Oliva et al., 2018). Warmer temperatures determine an ELA above (or very close to) the summits in the Picos de Europa and Sierra Nevada, therefore impeding glacial development in these massifs. In LIA glaciated environments, paraglacial processes have supplied abundant debris covering the remnants of those glaciers. In Picos de Europa, two of the six LIA glaciers are (semi-)permanent snow fields and the other four constitute (exposed or buried) ice patches (González-Trueba, 2005, 2006, 2007; González-Trueba et al., 2008; Serrano et al., 2011; Ruiz-Fernández et al., 2016, 2017). In the Sierra Nevada, the debris cover distributed inside the LIA moraine complexes shows evidence of subsidence and collapse, which is related to the degradation of the permafrost and buried ice masses located beneath them (Gómez-Ortiz et al., 2012b). In the same line, the kinematic monitoring of a rock glacier existing in the Veleta cirque revealed higher vertical than horizontal displacement rates, which is also indicative of the degradation of the frozen body existing in its interior (Gómez-Ortiz et al., 2014, 2018). In the case of the Pyrenees, the decrease of the glaciated surface from 1850 to present-day almost accounts for ca. 90% of the LIA area, decreasing from 2060 ha to 321 ha in 2008 (René, 2013), and to 244.6 ha in 2016 (Rico et al., 2017). This decrease has not been continuous, with several periods of glacial advance, stabilization and retreat, as also recorded in several other mid-latitude environments such as the Alps (Zemp et al., 2015). However, the strong disequilibrium between glaciers and present-day climate conditions has favoured an accelerated shrinking of the ice volume stored in the still glaciated massifs over the last two decades (Chueca et al., 2007; López-Moreno et al., 2016). The rate of shrinkage depends on the local topographical setting and microclimatic conditions, with solar radiation playing the major role on present-day glacier mass balances (López-Moreno et al., 2006). Some LIA glaciers have originated semi-permanent snow fields or ice patches whereas in other glaciated cirques during the LIA there is only a thick debris cover resulting from intense postLIA paraglacial dynamics (Serrano et al., 2018). In the Pyrenees, as in other mountains of the Mediterranean region still encompassing small glaciers, only avalanche-fed glaciers – which are the most resilient to climate change (Hughes, 2018) –, are those more likely to survive over the next few decades in an area where temperatures will increase more than global average and precipitation is also expected to decrease (IPCC, 2013). 6Conclusions The location of the Iberian Peninsula between the Atlantic Ocean and the Mediterranean Sea, between polar and subtropical air masses and between the shifting major high and low pressure systems affecting the peninsula has led to significant changes in climate regimes during the Quaternary period. This strong climate 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 sensitivity depicted in the rough Iberian terrain has determined the environmental dynamics prevailing in the highest mountain ranges, resulting in changing glacial to periglacial regimes following climate oscillations. The impact of Quaternary glaciers on mid-latitude alpine landscapes has been profusely examined from different perspectives over the last two centuries. In the Iberian Peninsula, a deep knowledge on the spatial distribution of glacial features is generally known since the 1970s and 1980s, although the temporal framework of glacial activity has been a matter of debate until nowadays, when only exist few and very small remnants of the large glaciers extending across the main mountain systems during the Late Quaternary. The use of different dating techniques over time has conditioned our understanding of glacial stages in the Iberian mountains. Despite the need to improve the accuracy of these techniques, as well as to extend it to other unexplored areas and different glacial records (moraines, polished bedrock and erratic boulders), results show evidence of different contrasting spatio-temporal patterns for past glaciations resulting from the high sensitivity of Iberian climate. Today, the Iberian mountains are one of the best examined mid-mountain environments with regards to its glacial evolution. Whereas 14C and OSL usually reported older ages for the MIE of the Last Glacial Cycle, CRE data suggest that it took place almost synchronously to the LGM or slightly before it in most mountain regions. There is evidence in some ranges of the occurrence of past glaciations prior to the Last Glacial Cycle between ca. 130 and 170 ka, as well as of several phases of advance and retreat within this last glacial Pleistocene cycle. Following the massive deglaciation of most Iberian ranges around 19–20 ka, the coldest stages recorded until the Holocene, namely the OD and YD, also promoted the expansion of glaciers. Finally, during the present interglacial, glaciers have been almost inexistent in the Iberian mountains, only reappearing as small spots during the coldest periods, such as the LIA. Post-LIA warming favoured the gradual disappearance of those small glaciers, which only persist nowadays in the highest massifs of the Central Pyrenees undergoing an accelerated melting process. However, from a critical perspective, our current understanding of glacial processes in the Iberian Peninsula during the Quaternary still presents several gaps that need to be addressed by the scientific community in the near future: - Glacial chronology for certain periods needs to be improved (e.g. Holocene). Did the 8.2 ka event promoted glacial expansion in the Iberian Peninsula? - There is increasing evidence of the occurrence glaciations predating the MIE of the Last Glacial Cycle, though they have been only detected in the Pyrenees, NW ranges and Sierra Nevada. Is there geomorphic evidence from these periods also in other ranges? - Chronological data in some areas (e.g. NW ranges, Iberian Range) are scarce and disconnected both in space and time. - There is the need to validate results with other complementary dating techniques. In some areas (e.g. Iberian Range) the chronology is (almost uniquely) based on one dating method (CRE) and this should be contrasted with other dating techniques (OSL and 14C). - The impact of past glaciations has been mostly studied in mountain environments but the connection with other processes in the lowlands has not been addressed. What is the impact of glacial stages on other records located down-valley of the glaciated environments (lacustrine, fluvial)? - Environmental dynamics following the deglaciation have not been examined in detail. What processes control the geography, intensity and timing of paraglacial response and postglacial periglacial dynamics in the different Iberian mountain ranges formerly glaciated to a greater or lesser extent? Chronological data of glacial oscillations are a consequence of climate fluctuations, and therefore should be incorporated in other paleoclimate data sources to better understand climate variability during the Late Quaternary in the North Atlantic region, and particularly in southern Europe and the Mediterranean basin, one of the regions of the global climate system with the greatest climate sensitivity. 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 Acknowledgements Marc Oliva is supported by the Ramón y Cajal Program (RYC-2015-17597) and the Research Group ANTALP (Antarctic, Arctic, Alpine Environments; 2017-SGR-1102) funded by the Government of Catalonia through the AGAUR agency. Laura Rodríguez is supported by the Marie Curie‐ Clarín COFUND program financed jointly by Gobierno del Principado de Asturias and the 7th WP of the European Union / Marie Curie Actions (reference ACA‐17‐19). The work complements the research topics examined in the PALEOGREEN project (CTM2017-87976-P), ESPAS project (CGL2015-655698-R), CGL2016-78380-P and CGL2015-65813-R project of the Ministry of Economy, Industry and Competitiveness, Spain, and the NUNANTAR project (02/SAICT/2017 - 32002) of the Fundação para a Ciência e a Tecnologia, Portugal. Research conducted in the Central Range is also supported by the National Parks Autonomous Agency (OAPN, MAGRAMA, project 1092/2014). 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 References Agassiz, L., 1840. Etudes sur les glaciers. Neuchatel. Akçar, N., Yavuz, V., Ivy-Ochs, S., Kubik, P.W., Vardar, M., Schlüchter, C., 2007. Paleoglacial records from Kavron Valley, NE Turkey: field and cosmogenic exposure dating evidence. Quat. Int. 164-165, 170-183. Akçar, N., Yavuz, V., Ivy-Ochs, S., Reber, R., Kubik, P.W., Zahno, C., Schlüchter, C., 2014. Glacier response to the change in atmospheric circulation in the eastern Mediterranean during the Last Glacial Maximum. Quat Geochronol. 19, 17-41. Alonso, V., González-Suárez, J.J., 1998. Presencia de hielo glaciar en los Picos de Europa (Cordillera Cantábrica). El helero del Jou Negro. Cuatern. Geomorfol. 12(1-2), 35-44. Alonso, V., Trombotto-Liaudat, T., 2009. Periglacial geomorphology of El Miro area, Cantabrian Mountains, NW Spain. Z. Geomorphol. 53(3), 335-357. Alonso, J.L., Pulgar, J.A., García-Ramos, J.C., Barba, P., 1996. Tertiary basins and Alpine tectonics in the Cantabrian Mountains (NW Spain), in: Friend, P., Dabrios, C. (Eds.), Tertiary basins of Spain: the stratigraphic record of crustal kinematics. Cambridge University Press, Cambridge, pp. 214-227. Allen, J.R.M., Huntley, B., Watts W.A., 1996. The vegetation and climate of northwest Iberia over the last 14000 yr. J. Quat. Sci. 11, 125-147. Alley, R.B., Brook, E.J., Anandakrishnan, S., 2002. A northern lead in orbital band: north-south phasing of Ice-Age events. Quat. Sci. Rev. 21, 431-441. Andrés, N., Gómez-Ortiz, A., Fernández-Fernández, J.M., Tanarro, L.M., Salvador-Franch, F., Oliva, M., Palacios, D., 2018. Deglaciation timing and rock glacier origin in La Cerdanya Mountains and Puigmal massif, Southeastern Pyrenees. Boreas 47(4), 1050-1071. Andrieu, V., 1991. Dynamique du paléoenvironnement de la vallée montagnarde de la Garonne (Pyrénées centrales, France) de la fin des temps glaciaires à l'actuel. PhD Thesis, Toulouse 2. Andrieu, V., Hubschman, J., Jalut, G., Hérail, G., 1988. Chronologie de la déglaciation des Pyrénées françaises, dynamique de sedimentation et contenu pollinique des paleolacs: Application à l’interprétation du retrait glaciaire. Quaternaire 2-3, 55-67. Arenillas-Parra, M., Cobos-Campos, G., Navarro-Carballo, J., 2008. Datos sobre la nieve y los glaciares en las cordilleras españolas. El Programa ERHIN (1984-2008). Ministerio de Medio Ambiente y Medio Rural y Marino, Madrid. Arnaez-Vadillo, J., 1987. Formas y procesos en la evolución de vertientes de la Sierra de la Demanda (Sistema Ibérico). Cuad. Invest. Geogr. 13, 9-153. Arnáez-Vadillo, J., García-Ruiz, J.M., 1990. Mapa geomorfológico de España. Hoja 240 (Ezcaray). Geoforma Ediciones, Logroño. Balco, G., Stone, J.O., Lifton, N.A., Dunai, T.J., 2008. A complete and easily accessible means of calculating surface exposure ages or erosion rates from 10Be and 26Al measurements. Quat. Geochronol. 3, 174-195. Barrère, P., 1963. La période glaciaire dand l’Ouest des Pyrénées centrales Franco-Espagnoles. Bull. Soc. Geol. Fr. 7(5), 516-525. Bickel, L., Lüthgens, C., Lomax, J., Fiebig, M., 2015. The timing of the penultimate glaciation in the northern Alpine Foreland: new insights from luminescence dating. P. Geol. A. 126(4-5), 536-550. Bide, J.B., 1893. Deuxième excursion dans la Sierra Nevada. Annuaire du Club Alpin Français 20, 276-305. Blanco-Chao, R., Pérez-Alberti, A., Trenhaile, A.S., Costa-Casais, M., Valcárcel, M., 2007. Shore platform abrasion in a para-periglacial environment, northwestern Spain. Gemorphology 83, 136-151. Boissier, E., 1839. Voyage botanique dans le midi de l´Espagne pendant l´année 1837. Fundación Caja de Granada y Universidad de Málaga. Granada, 1995. Bordonau, J., 1992. Els complexos glacio-lacustres relacionats amb el darer cicle glacial als Pirineus. Geoforma Ediciones, Logroño. Böse, M., Lüthgens, C., Lee, J.R., Rose, J., 2012. Quaternary Glaciations of northern Europe. Qua. Sci. Reviews 44, 1-25. Bradley, R.S., Jones, P.D., 1992. Climate Since 1500 AD. Routeledge, London. Briet, L., 1902. Le glacé de Mont-Perdú. Revue Club Alpine Francaise, secc. Sud Ouest. 1500, 183-186. Briner, J.P., Svendsen, J.I., Mangerud, J., Lohne, Ø.S., Young, N.E., 2014. A 10Be chronology of southwestern Scandinavian Ice Sheet history during the Lateglacial period. J. Quat. Sci. 29, 370-380. Bromley, G., Putnam, A., Borns, H., Lowell, T., Sandford, T., Barrell, D., 2018. Interstadial Rise and Younger Dryas Demise of Scotland’s Last Ice Fields. Paleoceanography and Paleoclimatology 33, 412-429. Brosche K.U., 1982. Studien zu jungpleistozänen un holozänen sedimenten und fossilen böden in Küstengebiet von West-Galizien (NW Spanien). Eiszeitalter und Gegenwart 32, 63-80. 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 Bru, M.J., 1985. Estudi geomorfologic: el modelatge glacial d'un sector del Pirineu Central (Valls Ferrera i de Cardos). PhD Thesis, University of Barcelona. Burgos, F.J.A., 1985. Características morfológicas de la Sierra de la Demanda. Actas del I Coloquio sobre Geografía de La Rioja: Geografía Física 31-42. Buizert, C., Gkinis, V., Severinghaus, J.P., He, F., Lecavalier, B.S., Kindler, P., Leuenberger, M., Carlson, A.E., Vinther, B.M., Masson-Delmotte, V., White, J.W.C., Liu, Z.Y., Otte-Bliesner, B., Brook, E.J., 2014. Greenland temperature response to climate forcing during the last deglaciation. Science 345(6201), 11771180. Bullón, T., 2016. The upper Pleistocene on the northern face of the Guadarrama Mountains (central Spain): Palaeoclimatic phases and glacial activity. Geomorphology 268, 233-245. Butzer K.W., 1967. Geomorphology and stratigraphy of the Paleolithic site os Budiño (Pontevedra, Spain). Eiszeitalter und Gegenwart 33, 95-117. Butzer, K.W., Franzle, O., 1959. Observations on pre-Würm glaciations of the Iberian Peninsula. Z. Geomorphol. 3(1), 85-87. Cacho, I., Grimalt, O., Canals, M., 2002. Response of the Western Mediterranean Sea to rapid climatic variability during the last 50,000 years: a molecular biomarker approach. J. Mar. Syst. 33-34, 253-272. Cano J., Fumanal M. P., Ferrer C., Usera J., Blázquez A. M., Olmo, J., 1997. Evolución de la costa meridional de Galicia durante el Cuaternario superior, in Rodríguez-Vidal, J., (Ed.), Cuaternario Ibérico, AEQUA, Huelva, pp. 33-46. Carandell, J., de Llarena, J.G., 1918. El glaciarismo cuaternario en los Montes Ibéricos. Trabajos del Museo Nacional de ciencias Naturales. Serie Geología 22, 104 pp. Carrasco, R.M., Pedraza, J., Domínguez-Villar, D., Villa, J., Willenbring, J.K. 2013. The plateau glacier in the Sierra de Béjar (Iberian Central System) during its maximum extent. Reconstruction and chronology. Geomorphology 196, 83-93. Carrasco, R.M., Pedraza, J., Domínguez-Villar, D., Willenbring, J.K., Villa, J., 2015. Sequence and chronology of the Cuerpo de Hombre paleoglacier (Iberian Central System) during the last glacial cycle. Quat. Sci. Rev. 129, 163-177. Carrasco, R.M., de Pedraza, J., Willenbring, J.K., Karampaglidis, T., Soteres, R.L., Martín-Duque, J.F., 2016. Morfología glaciar del Macizo de Los Pelados-El Nevero (Parque Nacional de la Sierra de Guadarrama). Nueva interpretación y cronología. Bol. R. Soc. Esp. Hist. Nat. Sec. Geol. 110, 49-66. Castañón, J.C., Frochoso, M., 1992. La glaciación Würm en las montañas cantábricas, in: Cearreta, A., Ugarte, F.M. (Eds.), The late quaternary in the Western Pyrenean region. Universidad del País Vasco. Bilbao, pp. 319-332. Carrera-Gómez, P., Valcárcel, M., 2018. The geomorphological role of snow since the Little Ice Age in the Sierra de Ancares (NW Spain). Cuad. Invest. Geogr. 44(1), 171-185. Chueca, J., Peña, J.L., Lampre, F., García-Ruiz, J.M., Martí-Bono, C., 1998. Los glaciares del Pirineo aragonés: Estudio de su evolución y extensión actual. Universidad de Zaragoza, Zaragoza. Chueca, J., Julián, A., Saz-Sánchez, M.A., Creus, J., López-Moreno, J.I., 2005. Responses to climatic changes since the ‘Little Ice Age’ on Maladeta Glacier (Central Pyrenees). Geomorphology 68, 167-182. Chueca, J., Julián, A., López-Moreno, J.I., 2007. Recent evolution (1981–2005) of the Maladeta glaciers, Pyrenees, Spain: extent and volume losses and their relation with climatic and topographic factors. J. Glaciol. 53, 547-557. Clark, P.U., Dyke, A.S., Shakun, J.D., Carlson, A.E., Clark, J., Wohlfarth, B., Mitrovica, J.X., Hostetler, S.W., McCabe, A.M., 2009. The last glacial maximum. Science 325, 710-714. Clark, P.U., Shakun, J.D., Baker, P.A., Bartlein, P.J., Brewer, S., Brook, E., Carlson, A.E., Cheng, H., Kaufman, D.S., Liu, Z., Marchitto, T.M., Mix, A.C., Morrill, C., Otto-Bliesner, B.L., Pahnke, K., Russell, J.M., Whitlock, C., Adkins, J.F., Blois, J.L., Clark, J., Colman, S.M., Curry, W.B., Flower, B.P., He, F., Johnson, T.C., Lynch-Stieglitz, J., Markgraf, V., McManus, J., Mitrovica, J.X., Moreno, P.I., Williams, J.W., 2012. Global climate evolution during the last deglaciation. Proc. Natl. Acad. Sci. U. S. 109, E1134E1142. Crest, Y., Delmas, M., Braucher, R., Gunnell, Y., Calvet, M., ASTER Team, 2017. Cirques have growth spurts during deglacial and interglacial periods: Evidence from 10Be and 26Al nuclide inventories in the central and eastern Pyrenees. Geomorphology 278, 60-77. Costa-Casais M., 2001. Análise sedimentaria e reconstrucción da costa Atlántica de Galicia. PhD. Thesis, Universidade de Santiago de Compostela. Costa-Casais, M., Martínez-Cortizas, A., Pérez-Alberti, A., 1994. Caracterización de un depósito costero de la Ría de Muros-Noia (A Coruña, Galicia), in: Arnáez-Vadillo, J., García-Ruíz J.M., Gómez-Villar, A., 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182 2183 (Eds.), Geomorfología en España. III Reunión de Geomorfología. Sociedad Española de Geomorfología. Universidad de La Rioja. Logroño, pp. 356-368. Costa-Casais, M., Martínez-Cortizas, A., Pérez-Alberti, A., 1996. Tipo de depósitos costeiros antigos entre Cabo de Fisterra e o Cabo de Corrubedo, in: Pérez-Alberti, A., Martini, P., Chesworth, W., MartínezCortizas, A. (Eds.). Dinámica y Evolución de Medios Cuaternarios. Xunta de Galicia. Santiago de Compostela, pp. 417-430. Cowton, T., Hughes, P.D., Gibbard, P.L., 2009. Palaeoglaciation of Parque Natural Lago de Sanabria, northwest Spain. Geomorphology 108(3-4), 282-291. Dahl-Jensen, D., NEEM community, 2013. Eemian interglacial reconstructed from a Greenland folded ice core. Nature. 493(7433): 489-494. Daveau, S., 1971. La Glaciation de la Serra de Estrela. Finisterra 6(11), 5-40. Daveau, S., Ferreira, A. B., Ferreira, N., Vieira, G., 1997. Novas observações sobre a glaciação da Serra da Estrela. Estud. Quat. 1, 41-51. Dehnert, A., Preusser, F., Kramers, J.D., Akcar, N., Kubik, P.W., Reber, R., Schlüchter, C., 2010. A multidating approach applied to proglacial sediments attributed to the Most Extensive Glaciation of the Swiss Alps. Boreas 39, 620-632. Delmas, M., 2015. The last maximum ice extent and subsequent deglaciation of the Pyrenees: an overview of recent research. Cuad. Invest. Geogr. 41, 359-387. Delmas, M., Gunnell, Y., Braucher, R., Calvet, M., Bourlès, D., 2008. Exposure age chronology of the last glaciation in the eastern Pyrenees. Quat. Res. 69, 231-241. Delmas, M., Calvet, M., Gunnell, Y., Braucher, R., Bourlès, D., 2011. Palaeogeography and 10Be exposure age chronology of Middle and Late Pleistocene glacier systems in the northern Pyrenees: implications for reconstructing regional palaeoclimates. Paleogeogr. Paleoclimatol. Paleoecol. 305. 109-122. Delmas, M., Gunnell, Y., Calvet, M., 2014. Environmental controls on alpine cirque size. Geomorphology 206, 318-329. Denton, G.H., Alley, R.B., Comer, G.C., Broecker, W.S., 2005. The role of seasonality in abrupt climate change. Quat. Sci. Rev. 24, 1159-1182. Domínguez-Villar, D., Carrasco, R.M., Pedraza, J., Cheng, H., Edwards, R.L., Willenbring, J.K., 2013. Early maximum extent of paleoglaciers from Mediterranean mountains during the last glaciation. Sci Rep-UK. 3, 20-34. Dresch, J., 1937. De la Sierra Nevada au Grand Atlas, formes glaciaires et de nivation. Mélanges de géographie et d´orientalisme offerts à E.F. Gautiers. Tours. 194-212. Ehlers, J., Gibbard, P.L., Hughes, P.D., 2011. Introduction, in: Ehlers, J., Gibbard, P.L., Hughes, P.D., (Eds.) Quaternary Glaciations-Extent and Chronology. A Closer Look. Elsevier, Amsterdam, pp. 1-14. Ehlers, J., Hughes, P.D., Gibbard, P.L., 2016. The Ice Age. John Wiley & Sons, New Jersey. Emiliani, C., 1955. Pleistocene temperatures. J Geol. 63, 538-578. Emiliani, C., 1966. Paleotemperature analysis of Caribbean cores P6304-8 and P6304-9 and a generalized temperature curve for the past 425,000 years. J Geol. 74, 109-126. Engel, Z., Braucher, R., Traczyk, A., Laetitia, L., ASTER team, 2014. 10Be exposure age chronology of the last glaciation in the Krkonoše Mountains, Central Europe. Geomorphology 206, 107-121. Evans, I.S., 2006. Local aspect asymmetry of mountain glaciation: A global survey of consistency of favored directions for glacier number and altitudes. Geomorphology 73, 166-184. Federici, P.R., Ribolini, A., Spagnolo, M., 2016. The glacial history of the Maritime Alps from the Last Glacial Maximum to the Little Ice Age. Geol. Soc. Spec. Publ. 433, 433-439. Felicísimo, A.M., 1992. El clima de Asturias. Geografía de Asturias, vol 1. Editorial Prensa Asturiana, pp. 1732. Fernandes, M., Oliva, M., Palma, P., Ruiz-Fernández, J., Lopes, L., 2017. Glacial stages and post-glacial environmental evolution in the Upper Garonne valley, Central Pyrenees. Sci. Total Environ. 584, 12821299. Fernandes, M.; Palma, P.; Lopes, L.; Ruiz-Fernández, J.; Pereira, P., Oliva, M., 2018. Spatial distribution and morphometry of permafrost-related landforms in the Central Pyrenees and associated paleoclimatic implications. Quat. Int. 470(A), 96-108. Fernández-Fernández, J.M., Palacios, D., García-Ruiz, J.M., Andrés, N., Schimmelpfennig, I., Gómez-Villar, A., Santos-González, J., Álvarez-Martínez, J., Arnáez, J., Úbeda, J., Léanni, L., ASTER Team, 2017. Chronological and geomorphological investigation of fossil debris-covered glaciers in relation to deglaciation processes: A case study in the Sierra de la Demanda, northern Spain. Quat. Sci. Rev. 170, 232-249. 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 2226 2227 2228 2229 2230 2231 2232 2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 Nieuwendam, A., Ruiz-Fernández, J., Oliva, M., Freitas, C., Paio, V., Cruces, A., 2015. Postglacial landscape changes and cryogenic processes in the Picos de Europa (Northern Spain) reconstructed from geomorphological mapping and microstructures on quartz grains. Permafrost Periglacial. Process. 27, 96108. Nussbaum, F., 1949. Sur les traces des glaciers quaternaires dans la région de l’Aragón. Pirineos 13-14, 497518. Nussbaum, F., 1956. Observations morphologiques dans la région de la Noguera Pallaresa. Pirineos 39-42, 5797. Obermaier, H., 1914. Estudio de los glaciares de los Picos de Europa. Trabajos del Museo Nacional de Ciencias Naturales: Serie Geológica 9, 1-42. Obermaier, H., 1921. El glaciarismo cuaternario en el valle del río Ara y en el Parque Nacional de Ordesa (Pirineos). Boletín de la Sociedad Ibérica de Ciencias Naturales 20, 106-110. Obermaier, H., Carandell-Pericay, J., 1916. Contribución al estudio del glaciarismo cuaternario de la Sierra de Gredos, Trabajos del Museo Nacional de Ciencias Naturales. Serie Geológica, 14, pp. 54. Obermaier, H., Carandell-Pericay, J., 1917. Nuevos datos acerca de la extensión del glaciarismo cuaternario en la Cordillera Central. Extracto del Bol. Soc. Esp. Hist. Nat., XVII. Madrid. Oliva, M., 2009. Holocene alpine environments in Sierra Nevada (Southern Spain). Ph.D. Thesis. University of Barcelona. Oliva, M., 2018. The Little Ice Age, the climatic background of present-day warming in Europe. Cuad. Invest. Geogr. 44(1), 7-13. Oliva, M., Gómez-Ortiz, A., 2012. Late Holocene environmental dynamics and climate variability in a Mediterranean high mountain environment (Sierra Nevada, Spain) inferred from lake sediments and historical sources. Holocene 22(8), 915-927. Oliva M, Gómez-Ortiz A, Schulte L., 2010. Tendencia a la aridez en Sierra Nevada desde el Holoceno Medio inferida a partir de sedimentos lacustres. Bol. Asoc. Geogr. Esp. 52, 27-42. Oliva, M., Schulte, L., Gómez-Ortiz, A., 2011. The role of aridification in constraining the elevation range of Holocene solifluction processes and associated landforms in the periglacial belt of the Sierra Nevada (Southern Spain). Earth Surf. Process. Landf. 36(10), 1279-1291. Oliva, M., Gómez-Ortiz, A., Palacios, D., Salvador-Franch, F., Salvà-Catarineu, M., 2014a. Environmental evolution in Sierra Nevada (South Spain) since the Last Glaciation based on multi-proxy records. Quat. Int. 353, 195-209. Oliva, M., Gómez-Ortiz, A., Palacios, D., Salvador-Franch, F., Salvà-Catarineu, M., 2014b. El Cuaternario en el macizo de Sierra Nevada. Evolución paleoambiental y paisaje a partir de la interpretación de registros naturales y documentos de época. Scr. Nova, XVIII (472), 1-16. Oliva, M., Gómez-Ortiz, A., Palacios, D., Franch, F. S., Catarineu, M. S., 2015. El Cuaternario en el macizo de Sierra Nevada. Estado de la cuestión, in: Gómez-Ortiz, A., Salvador-Franch, F., Oliva, M., SavàCaterineu, M., (Eds.), Avances, métodos y técnicas en el estudio del periglaciarismo. Universidad de Barcelona, Barcelona, pp. 65.-78. Oliva, M., Serrano, E., Gómez-Ortiz, A., González-Amuchastegui, M.J., Nieuwendam, A., Palacios, D., Pellitero-Ondicol, R., Pérez-Alberti, A., Ruiz-Fernández, J., Valcárcel, M., Vieira, G., Antoniades, D., 2016a. Spatial and temporal variability of periglaciation of the Iberian Peninsula. Quat. Sci. Rev. 137, 176-199. Oliva, M., Gómez-Ortiz, A., Salvador-Franch, F., Salvà-Catarineu, M., Ramos, M., Palacios, D., Tanarro, L., Pereira, P., Ruiz-Fernández, J., 2016b. Inexistence of permafrost at the top of Veleta peak (Sierra Nevada, Spain). Sci. Total Environ. 550, 484-494. Oliva, M.; Ruiz-Fernández, J., Barriendos, M., Benito, G., Cuadrat, J.M., García-Ruiz, J.M., Giralt, S., GómezOrtiz, A.; Hernández, A., López-Costas, O., López-Moreno, J.I., López-Sáez, J.A., Martínez-Cortizas, A., Moreno, A., Prohom, M., Saz, M.A., Serrano, E., Tejedor, E., Trigo, R., Valero-Garcés, B.L., VicenteSerrano, S., 2018. The Little Ice Age in Iberian mountains. Earth-Sci. Rev. 177, 175-208. Ortega, L.I., Centeno, J.D., 1987. Nota sobre nuevos focos glaciares detectados en la Sierra de Neila (NW del Sitema Ibérico, España). Bol. R. Soc. Esp. Hist. Nat. Sec. Geol. 82(1-4), 147-150. Ortigosa, L.M., 1985. Glaciares rocosos en la Sierra de Cebollera (Sistema Ibérico). Actas del I Coloquio sobre Geografía de La Rioja. 55-67. Ortigosa, L.M., 1986. Geomorfología glaciar de la Sierra Cebollera (Sistema Ibérico). Instituto de Estudios Riojanos, Logroño. Palacios, D., Andres, N., Luengo, E., 2003. Distribution and effectiveness of nivaltion in Mediterranean mountains: Peñalara (Spain). Geomorphology 54, 157-178. 2597 2598 2599 2600 2601 2602 2603 2604 2605 2606 2607 2608 2609 2610 2611 2612 2613 2614 2615 2616 2617 2618 2619 2620 2621 2622 2623 2624 2625 2626 2627 2628 2629 2630 2631 2632 2633 2634 2635 2636 2637 2638 2639 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650 2651 2652 2653 2654 2655 Palacios, D., Marcos, J., Vázquez-Selem, L., 2011. Last Glacial Maximum and deglaciation of Sierra de Gredos, central Iberian Peninsula. Quat. Int. 233, 16-26. Palacios, D., Andrés, N., Marcos, J., Vázquez-Selem, L., 2012a. Glacial landforms and their paleoclimatic significance in Sierra de Guadarrama, Central Iberian Península. Geomorphology 139-140, 67-78. Palacios, D., Andrés, N., Marcos, J., Vázquez-Selem, L., 2012b. Maximum glacial advance and deglaciation of the Pinar Valley (Sierra de Gredos, Central Spain) and its significance in the Mediterranean context. Geomorphology 177-178, 51-61. Palacios, D., Andrés, N., Vieira, G., de Marcos, J., Vázquez-Selem, L., 2012c. Last Glacial Maximum and deglaciation of the Iberian Central System. Geophysical Research Abstracts 14, EGU20123738. Palacios, D., Gómez-Ortiz, A., de Andrés, N., Vázquez-Selem, L., Salvador-Franch, F., Oliva, M., 2015a. Maximum extent of Late Pleistocene glaciers and last deglaciation of La Cerdanya mountains, Southeastern Pyrenees. Geomorphology 231, 116-129. Palacios, D., Andrés, N., López-Moreno, J.I., García-Ruiz, J.M., 2015b. Late Pleistocene deglaciation in the upper Gállego Valley, Central Pyrenees. Quat. Res. 83, 397-414. Palacios, D., Gómez-Ortiz, A., Andrés, N., Salvador, F., Oliva, M., 2016a. Timing and new geomorphologic evidence of the Last Deglaciation stages in Sierra Nevada (southern Spain). Quat. Sci. Rev. 150, 110-129. Palacios, D., Andrés, N., Gómez-Ortiz, A., García-Ruiz, J.M., 2016b. Evidence of glacial activity during the Oldest Dryas in the mountains of Spain, in: Hughes, P.D., Woodward J.C (Eds.), Quaternary glaciation in the Mediterranean mountains, Geol. Soc. Spec. Publ. 433, 87-110. Palacios, D., García-Ruiz, J.M., Andrés, N., Schimmelpfennig, I., Campos, N. Leanni, L., ASTER Team, 2017a. Deglaciation in the central Pyrenees during the Pleistocene-Holocene transition: Timing and geomorphological significance. Quat. Sci. Rev. 162, 111-127. Palacios, D., Andrés, N., Gómez-Ortiz, A., García-Ruiz, J. M., 2017b. Evidence of glacial activity during the Oldest Dryas in the mountains of Spain. Geol. Soc. Spec. Publ. 433(1), 87-110. Palacios, D., Gómez-Ortiz, A., Alcalá-Reygosa, J., Andrés, N., Oliva, M., Tanarro, L. M., Salvador-Franch, F., Schimmelpfennig I, Fernández-Fernández J. M., Lénni L., ASTER Team. 2019. The challenging application of cosmogenic dating methods in residual glacial landforms: The case of Sierra Nevada (Spain). Geomorphology 325, 103-118. Palade, B., Palacios, D., Gómez-Ortiz, A., 2011. Los glaciares rocosos de Sierra Nevada y su significado paleoclimático. Una primera aproximación. Cuad. Invest. Geogr. 37(2), 95-118. Pallàs, R., Rodés, A., Braucher, R., Carcaillet, J., Ortuño, M., Bordonau, J., Bourles, D., Vilaplana, J.M., Masana, E., Santanach, P., 2006. Late Pleistocene and Holocene glaciation in the Pyrenees: a critical review and new evidence from 10Be exposure ages, south-central Pyrenees. Quat. Sci. Rev. 25, 29372963. Pallàs, R., Rodés, A., Braucher, R., Bourlès, D., Delmas, M., Calvet, M., Gunnell, Y., 2010. Small, isolated glacial catchments as priority targets for cosmogenic exposure dating of Pleistocene climate fluctuations, southeastern Pyrenees. Geology. 38, 891-894. Palma, P., Oliva, M., García-Hernández, C., Gómez-Ortiz, A., Ruiz-Fernández, J., Salvador-Franch, F., SalvàCatarineu, M., 2017. Spatial characterization of glacial and periglacial landforms in the highlands of Sierra Nevada (Spain). Sci. Total Environ. 584, 1256-1267. Panzer, W., 1926. Talentwicklung und Eiszeitklima in nordöstlichen Spanien. Abhandlungen der Senckenbergischen Naturforschenden Gesellschaft. 386. Pedraza, J., Fernández, P., 1981. Cuaternario del Mapa Geológico de de Bohoyo. Mapa Geológico de Bohoyo (577). Madrid. Pedraza, J., Carrasco, R.M., Domínguez-Villar, D., Villa, J., 2013. Late Pleistocene glacial evolutionary stages in the Gredos mountains (Iberian Central System). Quat. Int. 302, 88-100. Penck, A., 1883. La période glaciaire dans les Pyrénées. Bulletin de la Société d’Histoire Naturelle de Toulouse 19, 105-200. Penck, A., Brückner, E., 1909. Die Alpen im Eiszeitalter. 3, Taunitz, Leipzig. Pellicer, F., 1984. Geomorfología de las Cadenas Ibéricas entre el Jalón y el Moncayo. Cuad. Estud. Borj. 11, 9-378 Pellitero-Ondicol, R., 2013. Evolución finicuaternaria del glaciarismo en el macizo de Fuentes Carrionas (Cordillera Cantábrica), propuesta cronológica y paleoambiental. Cuatern. Geomorfol. 27(1-2), 71-90. Pérez-Alberti, A., 1979. Nuevas observaciones sobre glaciarismo y periglaciarismo en el NW de la Península Ibérica. La Galicia sud oriental. Acta Geol. Hisp. 14, 441-444. Pérez-Alberti, A., 2014. Geomorfología, in: Ferreiro, D.A., Macías, F., Calvo de Anta, R., Pérez-Alberti, A., Otero-Pérez, X.L., Verde-Vilanova, J.R., Pérez-Llaguno, C., Saiz-Rubio, R., Díez, E., Bolaños-Guerrón, 2656 2657 2658 2659 2660 2661 2662 2663 2664 2665 2666 2667 2668 2669 2670 2671 2672 2673 2674 2675 2676 2677 2678 2679 2680 2681 2682 2683 2684 2685 2686 2687 2688 2689 2690 2691 2692 2693 2694 2695 2696 2697 2698 2699 2700 2701 2702 2703 2704 2705 2706 2707 2708 2709 2710 2711 2712 2713 2714 D., (Eds.), Complejo Básico-Ultrabásico de Capelada-Cabo Ortegal. Andavira Editorial, Santiago de Compostela, pp. 24-30. Pérez-Alberti, A., 2018. Xeomorfoloxia das Montañas do Courel and Mapa Xeomorfolóxico do Courel. Grupo de Desenvolvemento Rural Ribeira Sacra-Courel, 70 pp. Pérez-Alberti, A., Blanco-Chao, R., 1995. Importancia paleoclimática de las formas y depósitos de origen glaciar y periglaciar en la costa noroccidental de la Península Ibérica. (Serra da Capelada, Cedeira, A Coruña). Simposio Internacional Paleoambiente Cuaternario en la Península Ibérica. Resúmenes. Santiago de Compostela. Pérez-Alberti, A., Valcárcel, M., 1998. Caracterización y distribución espacial del glaciarismo en el Noroeste de la Península Ibérica, in: Gómez-Ortiz, A., Pérez-Alberti, A. (Eds.), Las huellas glaciares de las montañas españolas. Universidade de Santiago de Compostela. Santiago de Compostela, pp. 17-62. Pérez-Alberti, A., Rodríguez-Guitián, M.A, Valcárcel, M., 1992. El modelado glaciar en la vertiente oriental de la Sierra de Ancares (Noroeste de la Península Ibérica). Pap. geogr. 18, 39-51. Pérez-Alberti, A., Rodríguez-Guitián, M.A., Valcárcel, M., 1993. Las formas y depósitos glaciares en las Sierras Orientales y Septentrionales de Galicia (NW Península Ibérica), in: Pérez Alberti, A., Ramil-Rego, P., Guitián, L., (Eds.), La evolución del paisaje en las montañas del entorno de los Caminos Jacobeos. Xunta de Galicia. Santiago de Compostela, pp. 61-90. Pérez-Alberti, A., Valcárcel, M., Martini, P.I., 2002. La deglaciación en la vertiente suroccidental del Macizo de Peña Trevinca (Montañas de Trevinca, NW Ibérico), in: Serrano, E., García, A., Morales, J.C., Ortega, M.Y. (Eds.), Estudios recientes (2000-2002) en Geomorfología, Patrimonio, Montaña, Dinámica Territorial, pp. 205-216. Pérez-Alberti, A., Valcárcel, M., Blanco, R., 2004. Pleistocene glaciation in Spain, in: Ehlers, J., Gibbard, P.L. (Eds.), Quaternary GlaciationsExtent and Chronology: Part I: Europe (2). Elsevier, London, pp. 389394. Pérez-Alberti, A., López-Bedoya, J., Cunha, P.P., 2009. Sedimentological analysis of cold-climate stratified slope deposits of Galicia, NW Iberia. In IAS 2009 Meeting: Book of Abstracts vol. 6 (nº 6) Alghero, Italy. Pérez-Alberti, A., Valcárcel, M., Martini, P.I., Pascucci, V., Andreucci, S., 2011. Upper Pleistocene glacial valley-junction sediments at Pias, Trevinca Mountains, NW Spain. Geol. Soc. Spec. Publ. 354(1), 93-110. Preusser, F., Schlüchter, C., 2004. Dates from an important early Late Pleistocene ice advance in the Aare Valley, Switzerland. Eclogae Geol. Helv. 97(2), 245-253. René, P., 2013. Le réchauffement climatique en images. Cairn, Pau. Renssen, H., Seppä, H., Heiri, O., Roche, D.M., Goosse, H., Fichefet, T., 2009. The spatial and temporal complexity of the Holocene thermal maximum. Nat. Geosci. 2, 411-414. Renssen, H., Goosse, H., Roche, D., Sepp, H., 2018. The global hydroclimate response during the Younger Dryas event. Quat. Sci. Rev. 193, 84-97. Reuther, A.U., Urdea, P., Geiger, C., Ivy-Ochs, S., Niller, H.P., Kubik, P.W., Heine, K., 2007. Late Pleistocene glacial chronology of the Pietrele Valley, Retezat Mountains, Southern Carpathians constrained by 10Be exposure ages and pedological investigations. Quat. Int. 164-165, 151-169. Rico, I., Izagirre, E., Serrano, E., López-Moreno, J.I., 2017. Current glacier area in the Pyrenees: An updated assessment 2016. Pirineos 172, 10-3989. Rodés, A., Pallàs, R., Braucher, R., Bourlès, D., 2008. La última deglaciación en los Pirineos a partir de la datación de superficies de exposición mediante 10Be. Geo-Temas 10, 755-758. Rodrigues, T., Grimalt, J.O., Abrantes, F., Naugton, F., Flores, J.A., 2010. The last glacial-interglacial transition (LGIT) in the western mid-latitudes of the North Atlantic: abrupt sea surface temperature changes and sea level implications. Quat. Sci. Rev. 29(15-16), 1853-1862. Rodríguez-Rodríguez, L., Jiménez-Sánchez, M., Domínguez-Cuesta, M.J. Rico, M.T., Valero-Garcés, B.L., 2011. Last deglaciation in northwestern Spain: New chronological and geomorphologic evidence from the Sanabria region. Geomorphology 135, 48-65. Rodríguez-Rodríguez, L., Jiménez-Sánchez, M., Domínguez-Cuesta, M.J., Rinterknecht, V., Pallàs, R., Bourlès, D., Valero-Garcés, B.L., 2014. A multiple dating-method approach applied to the Sanabria Lake moraine complex (NW Iberian Peninsula, SW Europe). Quat. Sci. Rev. 83, 1-10. Rodríguez-Rodríguez, L., Jiménez-Sánchez, M., Domínguez-Cuesta, M.J., Aranburu, A., 2015. Research history on glacial geomorphology and geochronology of the Cantabrian Mountains, north Iberia (4342ºN/7-2ºW). Quat. Int. 364, 6-21. Rodríguez-Rodríguez, L., Jiménez-Sánchez, M., Domínguez-Cuesta, M.J., Rinterknecht, V., Pallàs, R., Bourlès, D., 2016. Chronology of glaciations in the Cantabrian Mountains (NW Iberia) during the Last Glacial Cycle based on in situ-produced 10Be. Quat. Sci. Rev.138, 31-48. 2715 2716 2717 2718 2719 2720 2721 2722 2723 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734 2735 2736 2737 2738 2739 2740 2741 2742 2743 2744 2745 2746 2747 2748 2749 2750 2751 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761 2762 2763 2764 2765 2766 2767 2768 2769 2770 2771 2772 2773 Rodríguez-Rodríguez, L., Jiménez-Sánchez, M., Domínguez-Cuesta, M.J., Rinterknecht, V., Pallàs, R., ASTER Team, 2017. Timing of last deglaciation in the Cantabrian Mountains (Iberian Peninsula; North Atlantic Region) based on in situ-produced 10Be exposure dating. Quat. Sci. Rev. 171, 166-181. Rodríguez-Rodríguez, L., Domínguez-Cuesta, M.J., Rinterknecht, V., Jiménez-Sánchez, M., GonzálezLemos, S., Léanni, L., Sanjurjo, J., Ballesteros, D., Valenzuela, P., Llana-Fúnez, S., ASTER Team, 2018a. Constraining the age of superimposed glacial records in mountain environments with multiple dating methods (Cantabrian Mountains, Iberian Peninsula). Quat. Sci. Rev. 195, 215-231. Rodríguez-Rodríguez, L., González-Lemos, S., Ballesteros, D., Valenzuela, P., Domínguez-Cuesta, M.J., Llana-Fúnez, S., Jiménez-Sánchez, M., 2018b. Timing of paraglacial rock-slope failures and denudation signatures in the Cantabrian Mountains (North Iberian Peninsula). Land Degrad. Dev. 29(9), 3159-3173. Ruddiman, W.F., McIntyre, A., 1981a. The Mode and Mechanism of the Last Deglaciation: Oceanic Evidence. Quat. Res. 16, 125-134. Ruddiman, W.F., McIntyre, A., 1981b. The North Atlantic Ocean during the Last Deglaciation. Paleogeogr. Paleoclimatol. Paleoecol. 35, 145-214. Ruiz-Fernández J, Poblete-Piedrabuena M.A, Serrano-Muela M.P., Martí-Bono, C., García-Ruiz, J.M., 2009. Morphometry of glacial cirques in the Cantabrian Range (Northwest Spain). Z. Geomorphol. 53(1), 4768. Ruiz Fernández, J., Poblete Piedrabuena, M.A., 2011. Las terrazas fluviales del río Cares: aportaciones sedimentológicas y cronológicas (Picos de Europa, Asturias). The fluvial terraces of the Cares river: sedimentological and chronological contributions (Picos de Europa, Asturias). Estud. Geogr. 72(270), 173-202. Ruiz-Fernández, J., Oliva, M., Cruces, A., Lopes, V., Freitas, C., García-García, C., Geraldes, M., López-Sáez, J.A., 2016. Environmental evolution in the Picos de Europa (Cantabrian Mountains, SW Europe) since the Last Glaciation. Quat. Sci. Rev. 138, 87-104. Ruiz-Fernández, J., Oliva, M., Hrbáček, F., Vieira, G., García-Hernández, C., 2017. Soil temperatures in an Atlantic high mountain environment: The Forcadona buried ice patch (Picos de Europa, NW Spain). Catena 149(2), 637-647. Sanchez-Gómez, S., 1990. Aplicación del estudio de suelos a la dinámica de la cuenca del río Lanjarón: relación suelos-geomorfología. PhD Thesis. University of Granada. Sancho, C., Arenas, C., Pardo, G., Peña-Monné, J.L., Rodés, E.J., Bartolomé, M., García-Ruiz, J.M., MartíBono, C., 2018. Glaciolacustrine deposits formed in an ice-dammed tributary valley in the south-central Pyrenees: New evidence for late Pleistocene climate. Sediment. Geol. 366, 47-66. Sanz Pérez, E.S, Pellicer-Corellano, F., 1994. El periglaciarismo de la Sierra de Cebollera (Soria y La Rioja). Cuatern. Geomorfol. 8(1-2), 123-138. Sarikaya, M.A., Zreda, M., Ciner, A., Zweck, C., 2008. Cold and wet Last Glacial Maximum on Mount Sandiras, SW Turkey, inferred from cosmogenic dating and glacier modelling. Quat. Sci. Rev. 27, 769789. Santos-González, J., Redondo-Vega, J.M., González-Gutiérrez, R.B., Gómez-Villar, A., 2013. Applying the AABR method to reconstruct equilibrium-line altitudes from the last glacial maximum in the Cantabrian Mountains (SW Europe). Paleogeogr. Paleoclimatol. Paleoecol. 387, 185-199. Sanz-Herraiz, C., 1988. El relieve del Guadarrama Oriental. Consejería de Política Territorial, Comunidad de Madrid. Madrid. Scapozza, C., Castelletti, C., Soma, L., Dall'Agnolo, S., Ambrosi, C., 2015. Timing of LGM and deglaciation in the southern Swiss Alps. Geomorphologie 4, 307-322. Schmieder O., 1915. Die Sierra de Gredos, Mitteilunger der Geographischen Gesellschaft in München, Zehnter Band, 1, 1-55. Schimmelpfennig, I., Benedetti, L., Garreta, V., Pik, R., Blard, P-H., Burnard, P., Bourles, D., Finkel, R., Ammon, K., Dunai, T., 2011. Calibration of cosmogenic 36Cl production rates from Ca and K-spallation in lava flows from Mt. Etna (38ºN, Italy) and Payun Matru (36ºS, Argentina), Geochim. Cosmochim. Acta. 75, 2611-2632. Schimmelpfennig, I., Schaefer, J. M., Putnam, A. E., Koffman, T., Benedetti, L., Ivy-Ochs, S., Aster Team and Schlüchter, Ch., 2014. 36Cl production rate from K-spallation in the European Alps (Chironico landslide, Switzerland). J. Quat. Sci. 29. 407-413. Schmitz, H., 1969. Glazialmorphologische Untersuchungen im Bergland Nordwestspaniens: (Galicien/León). Geographischen Inst. der Univ zu Köln im Selbstverl. Wiesbaden. Schrader, F., 1895. Sur l’étendue des glaciers des Pyrénées. Annuaire du Club Alpin Français 1894, 403-423. 2774 2775 2776 2777 2778 2779 2780 2781 2782 2783 2784 2785 2786 2787 2788 2789 2790 2791 2792 2793 2794 2795 2796 2797 2798 2799 2800 2801 2802 2803 2804 2805 2806 2807 2808 2809 2810 2811 2812 2813 2814 2815 2816 2817 2818 2819 2820 2821 2822 2823 2824 2825 2826 2827 2828 2829 2830 2831 2832 Schulte, L., 2002. Climatic and human influence on river systems and glacier fluctuations in southeast Spain since the Last Glacial Maximum. Quat. Int. 93, 85-100. Sermet, J., 1942. Sierra Nevada. Estudios Geográficos 9, 727-749. Serrano, E., 1998. Geomorfología del Alto Gállego (Pirineo Aragonés). Institución Fernando el Católico. Diputación de Zaragoza. Zaragoza. Serrano, E., 1991. Glacial evolution of the Upper Gallego Valley (Panticosa Mountains and Ribera de Biescas, Aragonese Pyrenees, Spain). Pirineos 138, 83-78. Serrano, E., 1992. Huellas de una glaciación pre-pleniglaciar en la Ribera de Biescas (Cuenca del Gállego, Pirineo aragonés), in: Libro Homenaje a Félix Ugarte, Eusko Ikaskuntza, 213-227, San Sebastián. Serrano, E., Agudo, D., 1988. La deglaciación del valle de los ibones Azules (Panticosa). Estudio glaciomorfológico. Cuatern. Geomorfol. 2, 115-123. Serrano, E., Martín-Moreno, R., 2018. Surge glaciers during the Little Ice Age in the Pyrenees. Cuad. Invest. Geogr. 44 (1), 213-244. Serrano, E., González-Trueba, J.J., Sanjosé, J.J., 2011. Dinámica, evolución y estructura de los glaciares rocosos de los Pirineos. Cuadernos de Investigación Geográfica 37(2), 145-170. Serrano, E., González-Trueba, J.J., González-García, M., 2012. Mountain glaciation and paleoclimate reconstruction in the Picos de Europa (Iberian Peninsula, SW Europe). Quat. Res. 78(2), 303-31. Serrano, E., González-Trueba, J.J., Pellitero, R., González-García, M., Gómez-Lende, M., 2013. Quaternary glacial evolution in the Central Cantabrian Mountains. Geomorphology 196, 65-82. Serrano, E., Gómez-Lende, M., González-Amuchástegui, M.J., González-García, M., González-Trueba, J.J., Pellitero, R., Rico, I., 2015. Glacial chronology, environmental changes and implications for human occupation during the upper Pleistocene in the eastern Cantabrian Mountains. Quat. Int. 364, 22-34. Serrano, E., González-Trueba, J.J., Pellitero, R., Gómez-Lende, M., 2017. Quaternary glacial history of the Cantabrian Mountains of northern Spain: a new synthesis. Geol. Soc. Spec. Publ. 433(1), 55-85. Serrano, E., Oliva, M., González-García, M., López-Moreno, J.I., González-Trueba, J.J, Martín-Moreno, R., Martín-Díaz, J., Gómez-Lende, M., Nofre, J., Palma, P., 2018. Post-Little Ice Age paraglacial processes and landforms in the high mountains of the Iberian Peninsula. Land Degrad. Dev. 29, 4186-4208. Serrat, D., 1977. L 'estudi geomorfològic del Pirineu oriental (Puigmal i Costabona). PhD Thesis, University of Barcelona. Serrat, D., Vilaplana, J.M., Martí-Bono, C.E., 1983. Some depositional models of glaciolacustrine environments in Southern Pyrenees, in: Evenson, E., Schluchter, Ch., Rabassa, J.P., (Eds.), Tills and related deposits. Balkema, Rotterdam, pp. 231-244. Serrat, D., Bordonau, J., Bru, J., Furdada, G., Gómez-Ortiz, A., Martí, J., Martí, M., Salvador-Franch, F., Ventura, J., Vilaplana, J.M., 1994. Síntesis cartográfica del glaciarismo surpirenaico oriental, in: Martí Bono, C., García Ruiz, J.M., (Eds.), El glaciarismo surpirenaico: nuevas aportaciones. Geoforma Ediciones, Logroño, pp. 9-15. Sermet, J., 1942. Sierra Nevada. Estud. Geogr. 3(9), 727. Shackleton, N.J., 1967. Oxygen isotope analyses and Pleistocene temperatures re-assessed. Nature 215, 15-17. Shakun, J.D., Clark, P.U., He, F., Lifton, N.A., Liu, Z., Otto-Bliesner, B.L., 2015. Regional and global forcing of glacier retreat during the last deglaciation. Nat. Commun. 6, 8059. Solé-Sabarís, L., 1941. La Canal de Berdún. Primera Reunión de Estudios Geográficos. Instituto Juan Sebastián Elcano, Madrid, pp. 113-160. Solé Sabarís, L., 1951. Los Pirineos. El medio y el hombre. Editorial Martín, Barcelona. Solomina, O.N., Bradley, R.S., Hodgson, D.A., Ivy-Ochs, S., Jomelli, V., Mackintosh, A.N., Nesje, A., Owen, L.A., Wanner, H., Wiles, G.C., Young, N.E., 2015. Holocene glacier fluctuations. Quat. Sci. Rev. 111, 934. Solomina, O.N., Bradley, R.S., Jomelli, V., Geirsdottir, A., Darrell, S K., Koch, J., McKay, N.P., Masiokas M., Miller, G., Nesje, A., Nicolussi, K., Owen, L.A., Putnam, A.E., Wanner, H., Wiles, G., Yan, B., 2016. Glacier fluctuations during the last 2000 years. Quat. Sci. Rev. 149, 61-90. Stone, J.O., 2000. Air pressure and cosmogenic isotope production. J. Geophys. Res. 105, 23753-23759. Stone, J.O., Fifield, L.K., Vasconcelos, P., 2005. Terrestrial Chlorine-36 Production from Spallation of Iron, in Abstracts of 10th International Conference on Accelerator Mass Spectrometry, September 5-10, Berkeley. Stokes, C.R., 2017. Deglaciation of the Laurentide Ice Sheet from the Last Glacial Maximum. Cuad. Invest. Geogr. 43(2), 377-428. 2833 2834 2835 2836 2837 2838 2839 2840 2841 2842 2843 2844 2845 2846 2847 2848 2849 2850 2851 2852 2853 2854 2855 2856 2857 2858 2859 2860 2861 2862 2863 2864 2865 2866 2867 2868 2869 2870 2871 2872 2873 2874 2875 2876 2877 2878 2879 2880 2881 2882 2883 2884 2885 2886 2887 2888 2889 2890 2891 Stroeven, A.P., Hättestrand, C., Kleman, J., Heyman, J., Fabel, D., Fredin, O., Goodfellow, B.W., Harbor, J.M., Jansen, J.D., Olsen, L., Caffee, M.W., Fink, D., Lundqvist, J., Rosqvist, G.C., Strömberg, B., Jansson, K.N., 2016. Deglaciation of Fennoscandia. Quat. Sci. Rev. 147, 91-121. Taboada, J., 1913. El lago de San Martín de Castañeda. Bol. Soc. Esp. Hist. Nat. 13, 860-883. Thornes, J. B., 1968. Glacial and periglacial features in the Urbion Mountains, Spain. Estud. geol. (Madr.). 24, 249-258. Tomkins, M.D., Dortch, J.M., Hughes, P.D., Huck, J.J., Stimson, A.G., Delmas, M., Calvet, M., Pallàs, R., 2018. Rapid age assessment of glacial landforms in the Pyrenees using Schmidt hammer exposure dating (SHED). Quat. Res. 90(1), 26-37. Trigo, R.M., Pozo-Vázquez, D., Osborn, T.J., Castro-Diez, Y., Gámis-Fortis S., Esteban-Parra, M.J., 2004. North Atlantic Oscillation influence on precipitation, river flow and water resources in the Iberian Peninsula. Int. J. Climatol. 24, 925-944. Turu, V., Calvet, M., Bordonau, J., Gunnell, Y., Delmas, M., Vilaplana, J.M., Jalut, G., 2016. Did Pyrenean glaciers dance to the beat of global climate events? Evidence of the Würmian sequence stratigraphy of an ice-dammed palaeolake depocentre in Andorra, in: Hughes, P.D., Woodward, J.C. (Eds.), Quaternary glaciation in the Mediterranean mountains, Geol. Soc. Spec. Publ. 433, 111-136. Uppala, S.M., Kållberg, P., Simmons, A., Andrae, U., Bechtold, V., Fiorino, M., Gibson, J., Woollen, J., 2005. The ERA-40 reanalysis. Q.J.R. Meteorological Soc. 131, 2961-3012. Valcárcel, M., 1998. Evolución geomorfológica y dinámica de las vertientes en el NE de Galicia. Importancia de los procesos de origen frío en un sector de las montañas lucenses. PhD Thesis. Universidade de Santiago de Compostela. Valcárcel, M., Pérez-Alberti, A., 2002 a. Un ejemplo de glaciarismo de baja cota en el NW de la Península Ibérica: el valle de Queixadoiro, in: Serrano, E., García, A., Morales, J.C., Ortega, M.Y. (Eds.), Estudios recientes en Geomorfología. Patrimonio, montaña y dinámica territorial. Universidad de Valladolid, Valladolid, 205-215. Valcárcel, M., Pérez-Alberti, A. 2002b. La glaciación finipleistocena en el sector noroccidental de las montañas leonesas: la Sierra de los Ancares, in: Redondo Vega, J.M., Gómez Villar, A., González Gutiérrez, R.B., Carrera Gómez, P., (Eds.), El modelado de origen glaciar en las montañas leonesas. Universidad de León, pp. 67-102. Valcárcel, M., Blanco-Chao, R., Feal-Pérez, A., 2009. Fases criológicas reconocidas en el NW de la Península Ibérica durante el Pleistoceno Final a partir de datos geomorfológicos: implicaciones en la modelización de las condiciones paleoambientales, in: Flor-Rodríguez, G., Gallastegui, J., Flor-Blanco, G., MartínLlaneza, J. (Eds.), 6º Simposio sobre el Margen Ibérico Atlántico MIA09. Universidad de Oviedo, pp. 165-168. Vasskog, K., Langebroek, P., Andrews, J.T., Nilsen, J.E.Ø., Nesje, A., 2015. The Greenland Ice Sheet during the last glacial cycle: Current ice loss contribution to sea-level rise from a palaeoclimatic perspective. Earth-Sci. Rev. 150, 45-67. Vegas, J., 2006. Los sistemas lacustres de las sierras de Neila y Urbión. Análisis sedimentológico y climático del Pleistoceno superior y Holoceno. PhD Thesis, Complutense University of Madrid. Vegas, J., 2007a. Análisis digital de imágenes de microscopía electrónica de ritmias siliclásticas. Interpretación limnogeológica de los registros lacustres de origen glaciar. (No sist. Ibérico). Cuatern. Geomorfol. 21, 49-75. Vegas, J., 2007b. Caracterización de eventos climáticos del Pleistoceno Superior-Holoceno mediante el estudio sedimentológico de la Laguna Grande (Sierra de Neila, NO Sistema Ibérico). Revista de la Sociedad Geológica de España 20(1-2), 53-70. Vieira, G., 2004. Geomorfologia dos planaltos e altos vales da Serra da Estrela. Ambientes frios do Plistocénico Superior e dinâmica actual. PhD Thesis, Universidade de Lisboa. Vieira, G., 2008. Combined numerical and geomorphological reconstruction of the Serra da Estrela plateau icefield, Portugal. Geomorphology 97(1-2), 190-207. Vieira, G., Ferreira, A.B. Mycielska-Dowgiallo, E. Woronko, B. Olszak, I., 2001.Thermoluminiscence dating of fluvioglacial sediments (Serra da Estrela, Portugal). V REQUI/I, Lisboa, pp. 85-92. Vieira, G., Palacios, D., 2010. New cosmogenic exposure dates for the Serra da Estrela glaciation. Preliminary results, in: Bateira, C., Soares, L., Gomes, A., Chaminé, H.I., (Eds.), Proceedings do V Congresso Nacional de Geomorfologia, pp. 14-25. Vidal-Romaní, J.R., Vilaplana, J.M., de Brum-Ferreira, A., Zezere, J.L., Rodrigues, L., Monge, С., 1990. Los tills de la Serra de Gerês-Xurés y la glaciación pleistocena (Minho, Portugal-Ourense, Galicia). Cuatern. Geomorfol. 4, 13-25. 2892 2893 2894 2895 2896 2897 2898 2899 2900 2901 2902 2903 2904 2905 2906 2907 2908 2909 2910 2911 2912 2913 2914 2915 2916 2917 2918 2919 2920 2921 2922 2923 2924 2925 2926 2927 2928 2929 2930 2931 2932 2933 2934 2935 2936 2937 2938 2939 2940 2941 2942 2943 2944 2945 2946 2947 2948 2949 2950 Vidal-Romaní, J.R., Santos-Fidalgo, L., 1994. La deglaciación finicuaternaria en el noroeste peninsular (Serra de Queixa-Invernadoiro, Ourense, Galicia): datos geomorfológicos y paleobotánicos. Cuatern. Geomorfol. 8(1-2), 33-44. Vidal-Romaní, J.R., Fernández Mosquera, D., Marti, K., Ferreira, A.B., 1999. Nuevos datos para la cronología glaciar pleistocena en el NW de la Península Ibérica. Cad. Lab. Xeol. Laxe. 24, 7-29. Vidal-Romaní, J.R., Fernández-Mosquera D., Marti, K., 2015. The glaciation of Serra de Queixa-Invernadoiro and Serra do Gerês-Xurés, NW Iberia. A critical review and a cosmogenic nuclide (10Be and 21 Ne) chronology. Cad. Lab. Xeol. Laxe. 38, 27-45. Vidal-Box, C., 1936. Contribución al conocimiento morfológico del segmento occidental de la Sierra de Gredos (Bohoyo). Bol. Soc. Esp. Hist. Nat. 36, 17-31. Vilaplana, J.M., 1983. Estudi del glacialisme de les valls de la Valira d'Ordino i d'Arinsal. PhD Thesis, University of Barcelona. Villa, E., Stoll, H., Farias, P., Adrados, L., Edwards, R.L., Cheng, H., 2013. Age and significance of the Quaternary cemented deposits of the Duje Valley (Picos de Europa, Northern Spain). Quat. Res. 79(1), 1-5. Walker, M.J.C., Berkelhammer, M., Björck, S., Cwynar, L.C., Fisher, D.A., Long, A.J., Lowe, J.J., Newnham, R.M., Rasmussen, S.O., Weiss, H., 2012. Formal subdivision of the Holocene Series/Epoch: a Discussion Paper by a Working Group of INTIMATE (Integration of ice-core, marine and terrestrial records) and the Subcommission on Quaternary Stratigraphy (International Commission on Stratigraphy). J. Quat. Sci. 27, 649-659. Williams, C., Flower, B.P., Hastings, D.W., 2012. Seasonal Laurentide Ice Sheet melting during the ‘Mystery Interval’ (15.5–14.5 ka). Geology 40, 955-958. Wirsig, C., Zasadni, J., Christl, M., Akçar, N., Ivy-Ochs, S., 2016. Dating the onset of LGM ice surface lowering in the High Alps. Quat. Sci. Rev. 143, 37-50. Zahno, C., Akçar, N., Yavuz, V., Kubik, P.W., Schlüchter, C., 2010. Chronology of Late Pleistocene glacier variations at the Uludag Mountain, NW Turkey. Quat. Sci. Rev. 29, 1173-1187. Zemp, M., Frey, H., Gärtner-Roer, I., Nussbaumer, S.U., Hoelzle, M., Paul, F., Haeberli, W., Denzinger, F., Ahlstrøm, A.P., Anderson, B., Bajracharya, S., Baroni, C., Braun, L.N., Cáceres, B.E., Casassa, G., Cobos, G., Dávila, L.R., Delgado Granados, H., Demuth, M.N., Espizua, L., Fischer, A., Fujita, K., Gadek, B., Ghazanfar, A., Hagen, J.O., Holmlund, P., Karimi, N., Li, Z., Pelto, M., Pitte, P., Popovnin, V.V., Portocarrero, C.A., Prinz, R., Sangewar, C.V., Severskiy, I., Sigursson, O., Soruco, A., Usubaliev, R., Vincent, C., 2015. Historically unprecedented global glacier decline in the early 21st century. J. Glaciol. 61(228), 745-762. 2951 2952 2953 2954 2955 2956 2957 2958 2959 2960 2961 2962 2963 2964 2965 2966 2967 2968 2969 2970 2971 2972 2973 2974 2975 2976 2977 2978 2979 2980 2981 2982 2983 2984 2985 2986 2987 2988 2989 2990 2991 2992 2993 2994 2995 2996 2997 2998 2999 3000 3001 3002 3003 3004 3005 3006 3007 3008 3009 Figure captions Figure 1. Location of the study sites in the different Iberian mountain ranges examined in this research. Figure 2. Number and type of absolute dating methods used in each mountain range. Figure 3. Glacial deposits and landforms from distinct glacial stages in the Pyrenees: (A) Main moraines and the 60 m fluvioglacial terrace corresponding to stages prior to LGM in the terminal basin of the Aragón valley glacier (M1 was dated at 171 OSL ka, and M2 at 68 OSL ka, whereas the fluvioglacial terrace was dated at 263 OSL ka); (B) LGM moraine at the front of the Malniu-Guils cirque complex, between the Duran valley to the west and the Querol valley to the East, Cerdanya, Eastern Spanish Pyrenees (the oldest boulders were dated at 23 10Be ka). The flat surface at the background culminates in the Puigpedrós Peak (2914 m a.s.l.); (C) The OD moraine, in a tributary of the Acherito valley, headwater of the Aragón Subordán valley glacier; (D) The YD moraine near the cirque headwall of the Escarra valley, a tributary of the Gállego valley glacier; (E) Mid-Holocene and LIA moraines in the Marboré cirque, at the foot of the Monte Perdido peak (3355 m); (F) LIA moraines above the YD complexes in the Maladeta massif. Figure 4. (A) Geomorphological map with the location of the main fluvial terrace sand moraine systems in the terminal basin of the Aragón River glacier; (B) Cross profile with the location of the main fluvial and glacial features. Figure 5. Examples of glacial landforms deposited during the Last Glaciation in the Cantabrian Mountains: (A) Local MIE central moraine developed between the Enol and Ercina lakes in Lagos de Covadonga, Western Massif of Picos de Europa; (B) Lateral moraine complex and related kame terrace deposit marking the local MIE and LGM glacial stages in the Monasterio valley, Redes Natural Park; (C) Erratic boulders deposited during the last deglaciation cold stages on the Fronfría ice-moulded surface in the Porma valley; (D) Recessional moraine and foremost ridge of a rock glacier deposited during the OD at the Valdevezón cirque in the Monasterio valley, Redes Natural Park; (E) Moraine possibly formed during either the YD or the Holocene at the base of Peña Agujas peak, Sierra de Sentiles; (F) Moraine formed during the LIA inside the Jou Negro cirque. Figure 6. Examples of landforms deposited during thein the NW Ranges: (A) Erratic boulders preserved outside the limits of the moraines of the Last Glaciation in the Tera valley, Trevinca massif; (B) Recessional moraine complex of Sanabria Lake, deposited during the last deglaciation in the Tera valley, Trevinca massif; (C) Moraine system developed in the NW slope of Capelada massif and ascribed to periglacial conditions in the coastal areas of Galicia during the Last Glaciation; (D) As Lamas Lake dammed by recessional moraines dated at ca. 15 ka in Requixo valley, Manzaneda massif; (F) Protalus ramparts formed during the LIA near Pico Cuiña (1998 m), Ancares massif. Figure 7. Examples of moraine complexes in the Central Range: (A) Principal moraine and peripheral ridges in the Laguna cirque, eastern face of Peñalara peak (Sierra de Guadarrama), both from LGM age; (B) LGM, OD and YD moraines in Cuerpo de Hombre, Sierra de Béjar; (C) Pinar valley in Sierra de Gredos, with peripheral ridges of unknown ages, principal moraine from the LGM and several ridges from the OD; (D) Deglaciation moraine sequences in La Serra paleoglacier, La Nava Massif, Sierra de Gredos; (E) LGM, OD; YD and Holocene glacial features in La Galana cirque, Gredos. Figure 8. Examples of moraines generated during the different glacial stages in the Iberian range: (A) Overhead view of the complete moraine sequence in the SE San Lorenzo cirque; (B) Oblique field view of the intermediate moraine in such cirque (LGM?) enclosing the fossil debris-covered glacier (DCG), and an older moraine corresponding to an earlier stage (pre-LGM?); (C) Overhead view of the complete moraine sequence in the westernmost Mencilla cirque; (D) Field view of the loose block accumulations interpreted as a fossil debris-covered glacier in such cirque, and the innermost moraine (OD?). Aerial orthophotos were obtained from the CNIG/IGN (http://www.ign.es/web/ign/portal). Figure 9. Examples of moraines generated during the different glacial stages in Sierra Nevada: (A) Moraines predating the Last Glaciation in Naute valley (dated at 135-140 ka); (B Moraines of the LGM and subsequent glacial advances in Hoya de la Mora cirque; (C) OD and YD moraines in San Juan valley; (D) Holocene moraine sequence at the foot of the highest Iberian peak (Mulhacén, 3478 m); LIA moraine ridge including two glacial advances at ca. 1350 and 1700 AD as revealed by surface exposure dating and historical sources. 3010 3011 3012 3013 3014 3015 3016 3017 3018 3019 3020 3021 3022 3023 3024 3025 3026 3027 3028 3029 3030 3031 3032 3033 3034 3035 3036 3037 3038 3039 3040 3041 3042 3043 3044 3045 3046 3047 3048 3049 3050 3051 3052 3053 3054 3055 3056 3057 3058 3059 3060 3061 3062 3063 3064 3065 3066 3067 3068 Figure 10. Summary of the chronological framework of glacial phases of the Last Glacial Cycle in Iberian mountains. Figure 11. Distribution of geomorphic evidence of glacial activity during MIS 6 in Iberian mountain ranges. Figure 12. For each of the mountain ranges, the figures include the dates of (a) the MIE of the Last Glaciation, (b) the MIS 4, (c) the MIS 2, and (d) the onset of Termination-1. Figure 13. Summary of the chronological framework of glacial stages during the Late Pleistocene-Holocene transition in the Iberian mountains. Figure 14. For each of the mountain ranges, the figures include the dates of glacial evidence during (a) the OD, (b) the BO, (c) the YD, (d) the Holocene, and (d) LIA. 3069 3070 3071 3072 3073 3074 3075 3076 3077 3078 3079 3080 3081 3082 3083 3084 3085 3086 3087 3088 3089 3090 3091 3092 3093 3094 3095 3096 3097 3098 3099 3100 3101 3102 3103 3104 3105 3106 3107 3108 3109 3110 3111 3112 3113 3114 3115 3116 3117 3118 3119 3120 3121 3122 3123 3124 3125 3126 3127 Table 1. Distribution of the different dating methods used to infer the deglaciation in each mountain range. Phase Pyrenees Cantabrian Mountains NW ranges Central Range Iberian Range Sierra Nevada Total Previous glaciations 0 / 0 / 4 / 04 / 0 / 0 / 23 / 0 / 0 / 0--3 / 0 / 0 / 09 / 0 / 4 / 2 Last Glacial Cycle 20 / 7 / 20 / 029 / 32 / 16 / 019 / 5 / 3 / 042 / 0 / 0 / 01 / 2 / 0 / 019 / 5 / 3 / 0118 / 46 / 39 / 0 Termination-1 72 / 4 / 1 / 027 / 6 / 0 / 19 / 10 / 0 / 025 / 0 / 0 / 05 / 2 / 0 / 09 / 10 / 0 / 0161 / 23 / 1 / 1 Holocene 30 / 0 / 0 / 01 / 1 / 0 / 3-5 / 0 / 0 / 011 / 0 / 0 / 00 / 0 / 0 / 060 / 4 / 0 / 3 LIA - - - - -3 / 0 / 0 / 33 / 0 / 3 / 0 Total 122 / 11 / 25 / 061 / 39 / 16 / 631 / 15 / 3 / 072 / 0 / 0 / 017 / 4 / 0 / 048 / 4 / 0 / 3351 / 73 / 47 / 6 CRE / C14 / OSL / other Table 2. Glacial activity in the Pyrenees during the Late Quaternary. Phase Chronology Environments, landforms and processes References Previous glaciations ~ 263 ka ~ 170-180 ka Some glacial deposits in the Gállego valley are located several km away from the main moraine system. Previous glaciations are well represented in the main and outermost moraine of the Aragón valley, as well as in fluvial deposits of the Aragón and Gállego valleys Martí-Bono (1978), Vilaplana (1983), Serrano (1992), Lewis et al. (2009), García-Ruiz et al. (2013) Last Glacial Cycle MIE occurred at 55-65 ka Some glaciolacustrine deposits show evidence of glacial expansion between 30 and 45 ka The LGM is well represented in the Eastern Pyrenees at ~21 ka Occurrence of a clear asynchroneity between the global LGM (MIS-2) and the regional MIE (MIS-4), particularly in the Central Pyrenees. The Gállego, Ara and Cinca Valleys show several examples of MIS-4 moraines. In the Eastern Pyrenees, there is increasing evidence of the presence of moraines of polygenic nature, as a result of the superposition of pulsations along the Last Glacial Cycle. García-Ruiz et al. (2003), GonzálezSampériz et al. (2006), Lewis et al. (2009), Pallàs et al. (2010), Turu et al. (2016), Andrés et al. (2018), Sancho et al. (2018), Guerrero et al. (2018) Termination1 Intense deglaciation following the LGM, at ca. 19 ka. Most of OD moraines were deposited between 16 and 17 ka. YD moraines between 11.5 and 13 ka Deglaciation in the Pyrenees is very well represented by a range of deposits whose distance to the headwater is directly related to the altitude of mountain massifs. The OD glacial expansion represented the growth of several km long glaciers leaving a number of morainic deposits. Development of rock glaciers at the end of the OD. A new, short glacial advance occurred during the YD, with deposition of big moraines close to the cirque limits. The longest glaciers were ca. 4 km long. Pallàs et al. (2006), Palacios et al. (2015a, 2015b, 2016a, 2017), GarcíaRuiz et al. (2016), Crest et al. (2017) Holocene Early to Mid-Holocene Most of rock glaciers that developed at the end of OD or YD melted out. Mid-Holocene: Reactivation of glaciers in the highest massifs, particularly Monte Perdido and Maladeta. A long and large moraine corresponding to Mid-Holocene has been identified in the Marboré cirque, northern face of the Monte Perdido massif. Probable existence of polygenic moraines formed during successive glacial advances in the Maladeta massif. García-Ruiz et al. (2014), Crest et al. (2017) LIA 15th-19th centuries: more than 100 small glaciers controlled by aspect and altitude Development of large, fresh moraines at the front of Quaternary cirques, mainly in massifs peaking at more than 3000. González-Trueba et al. (2008), GarcíaRuiz et al. (2014), Serrano and MartínMoreno (2018) Present-day Approx. 19 small glaciers remaining in north facing, well-protected slopes in the highest massifs. By 2016, glaciers occupy an area of 242 ha. Glaciers are clearly declining, with spatial and volumetric contraction. Some glaciers appear the more and more covered with debris. René (2013), López-Moreno et al. (2016), Rico et al. (2017) Table 3. Glacial activity in the Cantabrian Mountains during the Late Quaternary. Phase Chronology Environments, landforms and processes References Previous glaciations >394 ka (MIS 12?) ~173–131 ka (MIS 6) Precipitation of carbonate cements in old breccia deposits resting on a possible glacial-polished surface. Some erratic boulders in the terminal zone of the former Porma glacier suggest glaciation during MIS 6. Villa et al. (2013), Rodríguez-Rodríguez et al. (2016), Rodríguez-Rodríguez et al., (2018a) Last Glacial Cycle Pre-MIE at 110 ka (MIS 5d) MIE prior to 45–35 ka (MIS 3) LGM between 33–24 ka (MIS 2) The front of the Porma glacier reached an elevation of ~1110 m, during MIS 5d and MIS 3, leaving moraines and erratic boulders. Proglacial glacio-lacustrine sequences and kame terraces deposited outside the MIE lateral moraines provide minimum ages of ~45–33 cal ka BP in Picos de Europa, Redes, Somiedo-Laciana and Pas Mountains. Inherited boulders in the Porma lateral moraines in combination with OSL analysis of till and kame terraces suggest continuous glacial occupation from MIS 3 to the LGM of MIS 2. Jiménez-Sánchez and Farias (2002), Jalut et al. (2010), Moreno et al. (2010), Serrano et al. (2012), Frochoso et al. (2013), JiménezSánchez et al. (2013), Serrano et al. (2013), Ruíz-Fernández et al. (2016), RodríguezRodríguez et al. (2016, 2018a) Termination1 Onset of deglaciation between ~22–20 ka Glacier front stagnations until 14 ka Rock glacier stabilization at 15.7–13 ka Onset of paraglacial instabilities Abandonment of lateral moraines in the Porma valley and its main tributaries after ~22 ka. Glaciers were remarkably thinned and retreated from their MIE position by ~18–17 ka in the different massifs of the Cantabrian Mountains. Recessional moraines deposited in the Monasterio valley at ~19, 17.5, and 14.6 ka. Preliminar 10Be CRE dating of relict rock glaciers toes lying at 1600– 1550 m altitude suggest their stabilization between ~16.2–13.6 ka. Start of paraglacial rock slope failures in San Isidro valley at ~16–15 ka. Rodríguez-Rodríguez et al. (2016, 2017, 2018a, 2018b), Serrano et al. (2013), Allen et al. (1996) LIA Mid-19th century Six glaciers developed in the highest areas of the Central and Western massifs of Picos de Europa, leaving frontal moraines at elevations of 2200–2320 m Miotke (1968), González-Trueba (2005), González-Trueba et al. (2008), Present-day 2007 Only a few permanent ice-patches remain in the Central and Western massifs of Picos de Europa, fully or partly cover by debris mounds. González-Suárez and Alonso (1994), Frochoso and Castañón (1995), Alonso and González (1998), González-Trueba (2005), González-Trueba et al. (2007)