# Mammal extinction facilitated biome shift and human population change during the last glacial termination in East-Central Europe

Enikő Katalin Magyari, Mihály Gasparik, István Major, György Lengyel, Ilona Pál, Attila Virág, János Korponai, Aritina Haliuc, Zoltán Szabó, Piroska Pazonyi

► To cite this version:

Enikő Katalin Magyari, Mihály Gasparik, István Major, György Lengyel, Ilona Pál, et al. Mammal extinction facilitated biome shift and human population change during the last glacial termination in East-Central Europe. Scientific Reports, 2022, 12, 10.1038/s41598-022-10714-x. insu-03868226

HAL Id: insu-03868226

https://insu.hal.science/insu-03868226

Submitted on 23 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Mammal extinction facilitated biome shift and human population change during the last glacial termination in East-Central Europe

Enikő Katalin Magyari<fnref n="1" />,²,³✉, Mihály Gasparik⁴, István Major³, György Lengyel⁵, Ilona Pál³, Attila Virág<fnref n="1" />,⁶, János Korponai⁷,⁸,⁹, Aritina Haliuc¹⁰,¹¹, Zoltán Szabó² & Piroska Pazonyi<fnref n="1" />

The study of local extinction times, together with the associated environmental and human population changes in the last glacial termination, provides insights into the causes of mega- and microfauna extinctions. In East-Central (EC) Europe, groups of Palaeolithic humans were present throughout the last glacial maximum, but disappeared suddenly around 15,200 cal BP. In this study cave sediment profiles dated using radiocarbon techniques and a large set of mammal bones dated directly by AMS ¹⁴C were used to determine local extinction times. These were, in turn, compared to changes in the total megafauna population of EC Europe derived from coprophilous fungi, the Epigravettian population decline, quantitative climate models, pollen and plant macrofossil inferred climate, as well as to biome reconstructions. The results suggest that the population size of large herbivores decreased in the area after 17,700 cal BP, when temperate tree abundance and warm continental steppe cover both increased in the lowlands. Boreal forest expansion started around 16,200 cal BP. Cave sediments show the decline of narrow-headed vole and arctic lemming populations specifically associated with a tundra environment at the same time and the expansion of the common vole, an inhabitant of steppes. The last dated appearance of arctic lemming was at ~16,640 cal BP, while that of the narrow-headed vole at ~13,340, and the estimated extinction time of woolly mammoth was either at 13,830 (GRIWM) or 15,210 (PHASE), and reindeer at 11,860 (GRIWM) or 12,550 cal BP (PHASE). The population decline of the large herbivore fauna slightly preceded changes in terrestrial vegetation, and likely facilitated it via a reduction in the intensity of grazing and the concomitant accumulation of plant biomass. Furthermore, it is possible to conclude that the Late Epigravettian population had high degree of quarry-fidelity; they left the basin when these mammals vanished.

Despite the richness of vertebrate records from cave deposits and open air sites from the Late Pleistocene<fnref n="1" />, twentieth century paleontological research did not take the opportunity to date rigorously the latest appearance of large mammals and rodents during the last glacial termination in EC Europe (the eastern part of Central Europe). It also failed to date systematically cave mammal bone stratigraphies for the last 40,000 years, a time period for which AMS ¹⁴C dating would have been a perfect choice²,³. Only in recent years has a resurgence of

<footnote n="1">¹MTA-MTM-ELTE Research Group for Palaeontology, Budapest, Ludovika tér 2, 1083, Hungary. ²Department of Environmental and Landscape Geography, Eötvös Loránd University, Budapest, Pázmány Péter stny 1/c, 1117, Hungary. ³Isotope Climatology and Environmental Research Centre (ICER), Institute for Nuclear Research, Hungarian Academy of Science, Debrecen, Bem tér 18/c, 4026, Hungary. ⁴Department of Palaeontology and Geology, Hungarian Natural History Museum, Budapest, Ludovika tér 2, 1083, Hungary. ⁵Department of Prehistory and Archaeology, University of Miskolc, Miskolc-Egyetemváros 3515, Hungary. ⁶Department of Mineralogy and Geology, University of Debrecen, Debrecen, Egyetem tér 1, 4032, Hungary. ⁷Department of Water Supply and Sewerage, University of Public Service, Baja, Bajcsy-Zsilinszky utca 12-14, 6500, Hungary. ⁸Department of Environmental Sciences, Sapientia Hungarian University of Transylvania, Calea Turzii str. 4, 400193 Cluj-Napoca, Romania. ⁹Department of Danube's Diversity, Institute of Aquatic Ecology, Centre for Ecological Research, Budapest, Karolina út 29, 1113, Hungary. ¹⁰Institute of Speleology, Romanian Academy, 5 Clinicilor str, 400006 Cluj-Napoca, Romania. ¹¹EPOC, UMR 5805, Université de Bordeaux, Allée Geoffroy St Hilaire Bat B18N, CS 50023, 33615 Pessac Cedex, France. ✉email: eniko.magyari@ttk.elte.hu</footnote>
<figure><img src="image_3.png" /></figure>

**Figure 1.** Location of the studied cave sites (red, 1–2), pollen records (green, 3–7), dated reindeer (*Rangifer tarandus*) (orange, 1, 8–11) and dated woolly mammoth (*Mammuthus primigenius*) bones (yellow, 12–32) in the Carpathian Basin (a) and in East-Central Europe (b). 1: Jankovich Cave; 2: Rejtek I Rock Shelter.; 3: Taul dintre Brazi¹⁶; 4: Lake St Anne¹³; 5: Kelemér Nagymohos¹⁷; 6: Kokad Mire¹⁸; 7: Kardoskút Fehér Lake¹⁹; 8: Arka; 9: Ságvár; 10: Jászfelsőszentgyörgy; 11: Peskő Cave; 12: Csajág; 13: Pilismarót; 14: Feldebrő; 15: Szurdokpüspöki; 16: Budapest-Csillaghegy; 17: Szeged-Öthalom; 18: Esztergom-Gyurgyalag; 19: Tápiósüly; 20: Zók; 21: Mátraderecske; 22: Tarcal; 23: Bodrogkeresztúr; 24: Mende; 25: Tata-Porhanyó quarry; 26: Kiskunlacháza; 27: Ócsa-Felsőbabád; 28: Lágymányos; 29: Gyoma, River Tisza; 30: River Tisza; 31: Törökbecse; 32: Zebegény. The software used to create the maps was ArcGIS 10.2.2 for desktop, version 10.2.2.3552. Software url: https://suppo rt.esri.com/en/products/desktop/arcgis-desktop/arcmap/10-2-2.; the photo of Rejtek I Rock Shelter was taken by Tivadar Czina under licence CC BY-SA 3.0, source: https://hu.wikipedia.org/wiki/Rejteki_1._sz._kőfülke; the photo of Jankovich Cave was taken by Mihály Gasparik.

Upper Weichselian megafauna research been seen, with researchers revisiting key localities, often associated with human occupation, and reassessing their dating⁴⁻⁶. As our understanding of Late Epigravettian human lifestyle, economy and hunting strategy increases, the clearer it becomes that the disappearance of mobile Epigravettian groups from the Carpathian Basin (CB) around 15,200 cal BP is related to changing environmental conditions, such as the rapid disappearance of their main quarry, the reindeer (*Rangifer tarandus*), the decimation of secondary quarries such as the wild horse (*Equus ferus*), not to mention the worsening visibility for hunting in a landscape subject to an increase in forest cover⁷,⁸.

Here, the focus is on one species of megafauna, the woolly mammoth (*Mammuthus primigenius*), and several ungulate and rodent species that co-inhabited the CB in a fully developed glacial ecological setting, and what is more, one that is relatively well understood⁹⁻¹⁴.

Using a large dataset comprising dated mammoth bones, two profiles of caves with known human occupation and five multi-proxy paleoecological records from lake and mire sediments covering the last glacial maximum (LGM) and last glacial termination (Fig. 1), a specific paleoecological hypothesis relating rapid climate change to population dynamics is put to the test, namely, that transitions from cold to warm intervals were briefly optimal for grazing megafauna, but these brief optima were followed by rapid regional extinctions¹⁵. The question of the order of faunistic and vegetation biome changes and its casual linkage is also examined.

**Studied cave sites and the mammoth bone collections.** The caves under consideration are in Northern Hungary, in the Pre-Carpathian hill region (Fig. 1). The Jankovich Cave is an archaeological site in the Gerecse Hills. In the present study the Late Upper Palaeolithic layer is in focus²⁰,²¹ that likely belongs to the Late
Epigravettian. Situated at 330 m a.s.l., its 88 m long chamber has been excavated repeatedly since 1913; Block II, excavated by László Vértes in 1956, was selected for examination (Suppl. Fig. 1). This revealed that microfauna adapted to a cold climate were replaced by temperate species²¹, and was therefore assigned to the Pleistocene/Holocene boundary. 21 bone samples from 13 layers were selected for AMS ¹⁴C dating (Suppl. Tables 1 and 2).

The Rejtek I Rock Shelter is situated in the Bükk Mts at 534 m a.s.l. (Fig. 1, Suppl. Fig. 2). Here, the AMS ¹⁴C measurements were carried out on Profiles II (upper 4 layers) and III (layers 5–8)²². These two profiles had earlier been combined on the basis of depth²³, and faunistic and floristic changes were published accordingly (Suppl. Figs. 2 and 3; Suppl. Tables 1 and 3). Knapped lithics were found²⁴, charcoal, mollusc and bone assemblages were analysed⁴,²⁵,²⁶. 24 bone samples from 8 layers were selected for dating (Suppl. Tables 1 and 3). Priority was given to mammal species with missing or uncertain data on their last occurrence: e.g. steppe pika (*Ochotona pusilla*), arctic lemming (*Dicrostonyx torquatus*), narrow-headed vole (*Lasiopodomys (S.) gregalis*), reindeer (*R. tarandus*), and willow grouse (*Lagopus lagopus*).

The woolly mammoth was common in the CB during the Weichselian glaciation; *ca.* 400 specimens and 6 skeletons had been found in Hungary, but ¹⁴C dates were only available for 8 specimens²⁷. Within the scope of this study, 21 localities (Suppl. Table 4, Fig. 1) were selected for AMS ¹⁴C dating.

The cave localities, lake and mire sites used for the pollen-based climate and biome reconstructions in this study are described in detail in the Supplementary Material.

## Results

**Last glacial termination faunal changes and last detection times.** Stratigraphically consistent dates were obtained in both caves below a depth of 1 m of sediment (layers 8–12 in Rejtek, layers 5–10 in Jankovich), while in the top meter, mixed Holocene and Pleistocene ages indicated significant disturbance (Suppl. Tables 1, 2 and 3). The age range of the bone assemblages in the two caves is different: Jankovich layers 5–10 date between 17,550 and 15,300 cal BP, while the Rejtek layers 8–12 span 13,450–9950 cal BP. The Pleistocene-Holocene transition could only be traced at Rejtek. The age-depth models (Suppl. Fig. 4) suggest that sediment accumulation rates were relatively even in both caves, *ca.* 0.5 mm yr⁻¹.

Major vole faunal changes in the Jankovich profile suggest that the common vole (*Microtus arvalis*), narrow-headed vole (*L. (S.) gregalis*) and arctic lemming (*D. torquatus*) were dominant after the last glacial maximum (LGM: 26–21 ka cal BP) up until 16,360 cal BP. This assemblage was replaced by common vole (*M. arvalis*) and bank vole (*Clethrionomys glareolus*), which achieved dominance by 15,540 cal BP, suggesting significant warming and the disappearance of steppe-tundra habitats.

Radiocarbon dating confirms that this change started after Heinrich event 1 (HE-1) and had ended by the onset of the Bølling/Allerød Interstadial (GI-1e: ~ 14,700 cal BP) (Fig. 2; Suppl. Fig. 5). The onset of the Holocene is characterised by an increase in the forest dwelling bank vole (*C. glareolus*). A transition from wooded steppe tundra to continental steppe may be inferred as taking place between 15,540 and 16,000 cal BP (Fig. 2). The Jankovich vole record implies that the dominance of the arctic lemming and tundra vole (*D. torquatus* & *M. oeconomus*) between 17,000 and 16,300 cal BP may be connected to the HE-1 cooling (Suppl. Fig. 5). The steppe pika (*O. pusilla*) retreated after this interval (Suppl. Table 1).

The large mammalian fauna of Jankovich also indicates this climatic change: the cave bear (*Ursus spelaeus*), reindeer (*R. tarandus*), cave lion (*Panthera (Leo) spelea*), arctic fox (*Vulpes lagopus*) and horse (*Equus sp.*) that characterize the lower layers gradually disappear upwards in the stratigraphy. It should, however, be noted that the large mammalian bones were not directly dated using ¹⁴C.

Macrofauna, vole and rodent relative frequency from Rejtek I (Fig. 2; Suppl. Fig. 6) demonstrate that accumulation started in the Allerød (~GI-1b) and a major change is the shift in dominance from the common vole to the bank vole (*M. arvalis* - > *C. glareolus*) at 12,200–11,700 cal BP. This correlates with the Holocene transition. The relative frequencies of the European pine vole (*M. subterraneus*) and wood mouse (*Apodemus sylvaticus*) also increase indicating warming and forest cover increase. In the Allerød warm period, willow grouse (*L. lagopus*; dated to 14,610–15,190 cal BP; Suppl. Table 1), rock ptarmigan (*L. mutus*), steppe pika (*O. pusilla*), reindeer (*R. tarandus*) and steppe bison (*Bison priscus*; bone with low collagen) were still present, but they disappear in the Holocene layers (Fig. 2; Suppl. Fig. 6). This suggests a Late Glacial (LG) persistence of some species adapted to cold taiga, forest tundra and steppe habitats.

This cave sequence is exceptional in terms of the detailed charcoal record present in the same sediment layers as those in which the mammal bones were found. As shown in Suppl. Fig. 3, the dominant trees in the LG landscape of Rejtek were needle leaved: the charred remains of Norway spruce (*Picea abies*) and Arolla pine (*Pinus cembra*) were found together with Scots pine (*P. sylvestris*). From *ca.* 13,000 cal BP onwards, a dominance of mixed deciduous, maple (*Acer platanoides*), field elm (*Ulmus campestris*), hornbeam (*Carpinus betulus*) and beech (*Fagus sylvatica*) suggested a LG expansion of several early and late successional deciduous trees⁴,²⁵. Particularly striking was the appearance of beech and hornbeam. These spread later at 300 m a.s.l. in the Bükk Mts, around 6500 cal BP²⁸,²⁹, and their LG presence supports isoensyme and chloroplast (cp) microsatellite DNA studies inferring possible cryptic populations in the N Hungarian Hills³⁰. The mollusc fauna of Rejtek contained over 65% forest dwellers (*Clausiliidae*, *Cochlodina cerata*) in the LG and Early Holocene layers, indicating a forested landscape⁴,³¹.

¹⁴C dates suggest that the arctic lemming's (*D. torquatus*) population significantly decreased in the Gerecse Hills around 16,640 cal BP, followed by the narrow-headed vole (*L. (S.) gregalis*) around 13,340 cal BP, while the steppe pika (*O. pusilla*) survived into the Early Holocene; the latest detection time is around 11,015 cal BP at Rejtek (Suppl. Table 1; Fig. 3). Several reindeer (*R. tarandus*) bones were also dated; the latest date was from a stratigraphically consistent, undisturbed layer: ~ 15,195 cal BP (Suppl. Table 5; Suppl. Fig. 7). Taking into account other published data on reindeer ¹⁴C dates in the CB (15 sites), the Phase model in OxCal³² and the
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Vole community changes through time in two Carpathian Basin cave sediment profiles AMS ¹⁴C dated within the scope of this study (Jankovich Cave, Rejtek Rock Shelter); vegetation inferred from the vole community; reconstruction of the biome on the basis of five pollen records from the Carpathian Region (for locations see Fig. 1); last detection times of micro- and macrofauna elements in this study (from left: *Ochotona pusilla*, *Lasiopodomys* (S.) *gregalis*, *Rangifer tarandus*, *Mammuthus primigenius*, *Dicrostonyx torquatus*), pollen inferred summer mean temperature reconstructions (WAPLS, MAT) on the basis of the Kokad Mire pollen records (redrawn and modified from Magyari et al.¹⁸) and CCSM3 climate model based summer mean temperature simulation for the Kokad grid cell from PaleoView; animal drawings on this figure were made in Corel Draw version 23.0.0.363 by the first author.</figcaption></figure>

GRIWM model³³ were employed to estimate the last appearance of reindeer. The end boundary in the Phase model provides an estimate between 8300 and 14,480 cal BP, with a median at 12,550, while the GRIWM model places the extinction between 11,470 and 12,815 with an inverse weighted terminal age estimate of 11,860 cal BP (see Supplementary Material).

In the radiocarbon dating of the woolly mammoth (*M. primigenius*) bone assemblages from the CB (Suppl. Table 4; Suppl. Fig. 8), 12 of 40 bones had an insufficient collagen content, 16 samples were dated beyond 40,000 ¹⁴C BP, and 12 provided dates suitable for calibration.

Combined with previously published dates from the CB (Suppl. Fig. 8), the 2σ calibrated age ranges using Intcal20³⁵ showed that the youngest woolly mammoth was from Csajág³⁶, dated to 16,180 (15,830–16,010) cal BP. There are eight bones dated to <20 ka, nine between 20 and 30 ka, and twenty beyond 30 ka cal BP. Although the number of dated bones is a fraction of the bones available in museum collections, an absence between 22–24.5 ka and 27.4–32.5 ka cal BP is apparent. These data further suggest that woolly mammoth was common in the CB after the LGM, until 18 ka cal BP, though population decline was strong by 16.2 ka cal BP, and the species likely persisted locally during the coldest LGM period. The Phase model (n = 26) suggests a last appearance between 12,690 and 16,150 cal BP, with a median at 15,210 cal BP. The GRIWM model narrows down the confidence interval, but also estimates a considerably younger extinction time (13,717–13,960, terminal age 13,830 cal BP). If the problematic Zók site and burnt bones are removed from the Phase modelling (n=22), the estimated last appearance interval is still very similar (12,050–16,110 cal BP; median: 15,270).

**Last glacial termination vegetation changes in pollen and plant macrofossil records.** Pollen based biome reconstruction of four EC European localities (Fig. 2) suggests that the vegetation was diverse during and after the LGM, depending on elevation, water availability and soil type. Alluvial plains and mid eleva-
<figure><img src="image_5.png" /><figcaption>**Figure 3.** INTIMATE event stratigraphy plotted along the NGRIP and GISP2 δ¹⁸O and Ca²⁺ ice core records from Greenland³⁴; selected dung fungal spores (*Sporomiella*, HdV55a) indicating the presence of large bodied grazers in the Kokad pollen profile¹⁸; calibrated BP intervals of presence for selected mega- and microfauna elements in the Carpathian Basin dated within the scope of this study; summary pollen diagrams from two sites (Lake St Anne—950 m a. s. l.¹³ and Kokad Mire—112 m a. s. l.¹⁸) showing the cumulative percentage of total tree, shrub and herb pollen types and elm (*Ulmus*) pollen, and separate non-cumulative relative frequency curves for major herb pollen types (*Chenopodium*-type, *Artemisia*, *Poaceae*, *Cyperaceae*), palynological richness, microcharcoal influx, macrocharcoal inferred fire frequency curves, peak fire magnitudes and fire events (red crosses) from the same pollen sites. The methodology of macrocharcoal analyses is described in the Supplementary Material.</figcaption></figure>
tion mountain sites were covered by cool coniferous forest steppes (e.g. around Fehér Lake and Lake St Anne), locally tundra vegetation was dominant in cold wet basins (Nagymohos Peat Bog, Kokad Mire), while xerophytic woodland and scrub and cold steppe biomes were also present (Fig. 2).

Biome shifts begin to prevail earliest between 17,000 and 16,000 cal BP. In the mountains, cold deciduous forest, cool coniferous forest and forest steppe replaced the dominant steppe and taiga biome of the preceding HE-1 cooling. A much more abrupt biome shift from tundra to cool coniferous forest took place at Kokad around 16,200 cal BP (Fig. 2)¹⁸; this was followed by a second biome shift to cool mixed forest at 14,740 cal BP, a biome that then persisted into the Early Holocene here. A similar biome shift towards mixed leaved (boreo-nemoral) forest was also detected in the Eastern Carpathian site at 14,500 cal BP, suggesting the expansion of temperate deciduous trees during the LG interstadial (GI-1).

At the lowland Kokad site, grazing indicator fungi (Sporomiella, Sordariales and Sordaria-type (HdV55a)) were present on pollen slides. The disappearance of these fungi took place at 16,780 cal BP (Fig. 3)¹⁸. Around the same time, the micro- and macro-charcoal accumulation rate increased in the sediment and stayed high until ca. 15,000 cal BP, suggesting recurrent wildfires in the region between 17,700 and 15,000 cal BP (Fig. 3).

Comparing selected pollen types, palynological richness (representing floristic diversity here), evenness (a measure of abundance distribution within the pollen assemblages), microcharcoal accumulation rate (indicative of regional fires) and macrocharcoal-based fire frequency and peak magnitude records (indicative of local fires), several differences become apparent between the lowland and mid mountain localities (Fig. 3). First, warming at ~16,200 cal BP and even at 14,700 cal BP did not lead to massive afforestation in the lowland; instead, a modest increase in forest cover was accompanied by the expansion of grass-steppes. A second striking difference is in fire history. The lowland site was characterised by increased fire activity from 17,800 cal BP until the LG onset (14,700 cal BP), when broad-leaved trees became dominant in this region. The onset of fire increase pre-dates the estimated disappearance of the megafauna, but appears at the time of the drastic population decrease of these megaherbivores.

It also seems that floristic diversity correlates well with vegetation openness and vegetation type; in case of the lowland site, the sequence tundra → cool steppe → cool coniferous forest steppe resulted in an increase in palynological diversity and evenness; at the mid mountain site, however, palynological diversity decreased when closed mixed leaved and taiga forests started to dominate from the LG onset (~14,700 cal BP). Taxa accumulation curves (Suppl. Figs. 9 and 10) also indicate a change in the local species pool at the LG onset. This has implications for the grazing megafauna as well: first, the least diverse floristically was the steppe tundra (E(T) < 20), while the forest steppe biomes were more diverse (E(T) ~ 30). On the basis of the radiocarbon dates from woolly mammoth and reindeer (Suppl. Figs. 7 and 8), a decline in these species is observed when locally tundra (steppe tundra) disappeared.

Figure 2 shows pollen-inferred July mean temperatures at Kokad using two different transfer functions¹⁸. While the WA-PLS reconstruction indicates that between 19,440 and 16,700 cal BP pollen-inferred summer mean temperatures were stable and fluctuated between 14.5 and 15.7 °C, the MAT based reconstruction yields lower values and a high amplitude fluctuation. A distinct cooling is present between 17,900 and 16,280 cal BP, with inferred summer mean temperatures between 7.5 and 12.7 °C. This interval coincides with HE-1 (17,850–16,200 cal BP), and the reconstructed temperature values follow the summer mean temperature values indicated by the vole fauna at the Jankovich Cave²¹ (Fig. 2). This suggests that the MAT based reconstruction is likely to be realistic. Following this event, July mean temperatures increased rapidly to 16–18 °C. Both the vole and pollen based summer temperature reconstruction demonstrate that the highest amplitude warming took place before the LG onset, around 16.2 ka cal BP in the CB and summer mean temperatures decreased only modestly during the Younger Dryas climatic reversal (Fig. 2), corroborating earlier proxy based reconstructions from the CB lowlands and SE Carpathian mountain sites¹⁶,³⁷⁻³⁹.

## Discussion

**The order of events: ecosystem changes in the Carpathian Basin and Europe during the last glacial termination.** The multi-proxy data (Fig. 3) suggest that the collapse of the large herbivore fauna slightly preceded or coincided with the terrestrial vegetation change in the CB. In addition to the radiocarbon dated bones of woolly mammoth (Suppl. Fig. 8), fungal spores (Sporomiella, Sordaria-type (HdV55a), 1.4–230 spores cm⁻² yr⁻¹) living in the dung of large grazing mammals confirm a population decrease from 17.2 ka cal BP that likely culminated around 16.2 ka cal BP, contemporaneous with the last dated mammoth, even though the final disappearance is modelled to 13.8–15.2 ka cal BP. As demonstrated by the Kokad micro- and macro-charcoal records, regional and local fires appeared from 17.8 ka cal BP in the eastern lowland of the CB up until 14.3 ka cal BP, when the spread of temperate woody species started. A correlation is also apparent between the last ¹⁴C date for a mammoth and the transformation of the herbaceous vegetation at the Kokad locality. The picture is therefore less definite here than in North America, where fungal spore, charcoal and pollen studies have demonstrated that the population decrease of megaherbivores preceded the vegetation change⁴⁰. In the CB grazing indicator fungal spores are less frequent in the LP and LG deposits, rendering this less conclusive; in addition, radiocarbon dates suggest coincidence of the change in vegetation and megaherbivore latest radiocarbon dates in case of the woolly mammoth. Overall, the data presented here suggest that a previously large herbivore population probably decreased earlier than the biome shifts took place, and thus the limiting factor for large grazing mammals was not the decreasing availability of food resources, but their intolerance of increased warmth. Furthermore, regional wildfire histories were divergent in the mountain and lowland localities, depending largely on tree cover (Fig. 3). From regional summaries it is known that wildfires were generally more frequent during the Early Holocene in EC Europe, and in conifer forest dominated regions during the LG⁹,⁴¹. This increase in fires may be explained by orbital forcing (warmer than present summers, cold winters,
seasonal drought stress); however, the same studies also concluded that once temperate broadleaved forest had become established, biomass burning was high at ~45% tree cover and decreased to a minimum at between 60 and 70% tree cover. These data corroborate the determination of biomass burning by tree compositional change and tree cover. It was the highest at medium needle-leaved tree cover in the region.

The microfaunal changes and the vegetation changes inferred from the vole community from Jankovich and Rejtek I are the first rigorously dated cave sequences that yield information about ecosystem changes during the last glacial termination in Northern Hungary. Their merit is that they demonstrate that the point at which faunal turnover took place was not the LG onset, but the post HE-1 warming occurring at ~16.2 ka cal BP (see Jankovich Cave record on Fig. 2 and Suppl. Figs. 5 and 6).

Such high-resolution microfaunal records are rare in Europe. In Hungary, two cave faunal assemblages were revisited recently⁴,⁴², including Rejtek I and Petényi Cave in the Bükk Mts, where mollusc shells were used for dating the earlier sequences, which also contained undiagnostic Late Upper Palaeolithic artefacts, such as retouched bladelets. The dates based on mollusc shells are close to the oldest estimates given here (Suppl. Table 3), in the LG interstadial layers (11–12) of Rejtek I, and the mollusc $^{14}$C ages are also older, while in the Early Holocene layers there is a good agreement between the bone collagen and mollusc carbonate based $^{14}$C dates.

Overall, this discrepancy, plus the deposition hiatus that demonstrably took place during the LG interstadial (between ~14,530 and 9270 cal BP) in the Petényi Cave, and the older dates for some large grazer or predator bones in Rejtek I and Jankovich Cave to warn researchers that the interpretation of various cave sediment fauna has to be treated with caution, and the best dating results and faunal based inferences can be obtained from micromammal assemblages. Furthermore, the bone and charcoal assemblages in the bottom layer of Petényi Cave (dated between 15,180 and 14,530 cal BP) indicate the development of a transitional flora and fauna (boreo-nemoral forest dominated by spruce) in the Bükk Plateau⁴,²³,⁴³ supporting the pollen records from the CB that afforestation and warming started directly after HE-1.

In Europe systematically analysed and radiocarbon dated Late Pleistocene—Early Holocene cave sediment sequences with rich bone assemblages are rare. The few that cover a similar time period are in Western France⁴⁴,⁴⁵, Spain⁴⁶ and Ukraine⁴⁷. In addition, a recent summary work⁴⁸ compares ecosystem changes at a regional scale over the last 50,000 years.

The faunal turnover at 16.2 ka cal BP identified in Jankovich Cave has also been detected in two French cave sequences (Peyrazet and Coulet des Roches). In both cases, changes in the small mammal communities between the Pleistocene and Holocene were the result of a succession of climatic events starting at the end of HE-1. Several rodents that occur in temperate and forested habitats today (e.g. the garden dormouse (*Eliomys quercinus*), wood mouse (*A. sylvaticus*), bank vole (*C. glareolus*), Mediterranean pine vole (*Microtus* (Terricola) duodecimcostatus*)) appeared in the middle or at the end of the LG, while species adapted to cold (e.g. arctic lemming, narrow-headed vole) gradually disappeared⁴⁴,⁴⁵. The three radiocarbon dated reindeer bones from the Peyrazet Cave gave an age range 13,835–15,410 cal BP⁴⁴, i.e. slightly more recent than the results from Jankovich Cave, falling mainly within the Bølling interstadial. Similar faunal turnover and climatic changes were observed in the El Mirón Cave sequence⁴⁶, where an increase in forest species diversity was detected between 18 and 11 ka cal BP, while *Pliomys lenki* (an extinct Pleistocene vole) and some other species mainly disappeared at the end of the YD.

According to Puzachenko and Markova⁴⁸ the CB belongs to the Central European South region, where after the significant decrease in species richness of the LGM, the number of species was gradually restored to the value characteristic of the late MIS-3 between ~17.5 and 13 ka BP. If local last detection times from the Rejtek I and Jankovich are compared with the regional disappearance times reported, it is possible to conclude that all Late Pleistocene species disappeared earlier from the CB than from the rest of the region. Even though the number of systematically dated cave sequences is still low in the CB, and thus later local extinction times are plausible, these relatively early local last detection times may, in all likelihood, be connected to the southern geographical position of the basin, and also to the biome mosaic that characterised this region during the LPG and LG.

In comparison with the aforementioned areas, the Grot Skeliastyi Rock Shelter in south-western Crimea indicates a different Pleistocene-Holocene faunal transition⁴⁷. Only the large herbivore species became extinct from this assemblage, while most other taxa persisted from the Pleistocene into the Holocene without losses, a phenomenon likely to have been influenced by the Crimea's geographic position and milder climate. Its relevance to the CB lies in the vegetation and fauna of the south-eastern lowland areas that via the Iron Gate are directly connected to the Pontic Crimean territories and contain several common taxa⁴⁹. Many of the Pontic species originated from this climatically relatively stable area, and as the differences of the faunal records demonstrate, it is likely that Pontic species migrated into the CB during the last glacial termination period⁵⁰.

**Plaid ecosystems reverting to equilibrium ecosystem mosaics: key to steppe fauna survival.** Sommer and Nadachowski⁵¹ have demonstrated that faunal communities during the LGM contained a combination of cold and temperate faunal elements in the Balkans (except Greece), in SW France and in the Carpathian Region. Consequently, in these regions delayed expansion of new faunal communities in response to climate change were not influenced by delayed immigration. When the order of changes in different ecosystem components is examined¹²,¹⁷,⁵⁰, a similar situation applies to certain woody and herbaceous elements of the temperate forest and forest steppe flora.

According to the Plaids and Stripes Hypothesis⁵², the main cause of the Late Pleistocene megafauna extinction was the cessation of short-term climate fluctuation during the last ice age that added a dimension of temporal complexity, which is now missing from many modern ecosystems⁵³. Millennial- and centennial-scale high amplitude climate fluctuations kept ecosystems out of balance, as plant and animal species struggled to keep up with repeated shifts in their environments. Advantageous eco-physiological attributes, such as greater mobility, lower cost of locomotion, greater dietary breadth and higher metabolic efficiency allowed the Late Pleistocene
fauna to flourish in such plaid settings characterised by disequilibrium. Due to frequent climatic disturbance, early successional plant communities dominated by forbs and gaminoids were the key elements in supporting dense populations of megafaunal herbivores⁵⁴,⁵⁵ on immature, and thus more productive, soils. In support of this hypothesis, the transformation of the plaid ecosystems towards a striped structure during the LG and Early Holocene can be easily traced on Eurasian biome simulation maps⁵⁶.

The pollen based biome reconstruction, microfaunal change and climate model simulations presented here all suggest that the rapid transformation of the plaid landscape took place from 16.2 ka cal BP in the CB, and the major element was the expansion of temperate and boreal woodland species and the overall gradual increase in woodland cover that did not favour mammoth and reindeer. On the other hand, if current and simulated plant biomes in the CB are examined, it is possible to see that the striped boreal and cold temperate biomes break up in the lowlands of the basin, where the so called equilibrium mosaic ecosystem (with edaphic steppes and temperate forest steppe) persisted throughout the Holocene for edaphic and hydrological reasons⁴⁹.

This deviation from the regional trend probably had an overarching consequence during the Holocene climate stability, allowing the longer subsistence of mega-herbivore mammal species in the lowlands, as demonstrated by several studies³,⁵⁷,⁵⁸. Although climate change led to the replacement of the ungulate species spectrum due to partial afforestation, early warming around 16.2 ka cal BP was detrimental to the mammoth adapted to cold-steppe, tundra-steppe environments. As a very distant and indirect parallel, the lowlands of the CB are somewhat similar to the African savannah, where dryland ecosystems are particularly susceptible to millennial-scale boom-and-bust cycles in primary productivity⁵⁹, and therefore plaid ecosystems are pertinent.

**Local herbivore extinctions in context of European extinction history.** Reindeer. It is known from range modelling that the global potential range of reindeer declined by 84% between 21 and 6 kyr BP⁶⁰. This is explicable in terms of rapid climate change, particularly after HE-1. Starting from their modern July mean temperature tolerance of <12–13 °C and a metabolic adaptation to <15 °C⁶¹, their distribution dynamics in the CB suggest that reindeer were common in the CB during the LGM (from ca. 23 ka cal BP) and the population had declined steeply by ca. 15.2 ka cal BP (Suppl. Table 5), with a modelled last appearance at ~11,860–12,550 cal BP. The range dynamics of reindeer in Europe summarized recently⁶² showed only four context-dated reindeer findings from the CB, all of which were dated between 18 and 25 ka cal BP. Recently, the Zöld Cave from the central CB⁸ yielded a reindeer bone dated to 15.4–16 ka cal BP, and further reindeer remains were associated with a charcoal date 14.9–15.3 ka cal BP. In this context, the dating results of several reindeer bones from Jankovich Cave presented here with a calibrated (2σ) age range of 15,085–20,540 cal BP (Suppl. Table 5) suggest that reindeer persisted in the basin after the LGM, and their local extinction probably occurred later than assumed by Sommer et al.⁶².

The present data suggest that the reindeer population declined considerably around 15.2 ka cal BP. Comparing this timing with the climate reconstruction inferred from pollen and chironomid, together with biome records (Figs. 2 and 3)³⁸, it is possible to conclude that reindeer persisted in the cool conifer forest steppe environment of the basin for about 1000 years, and their population decline predated the emergence of cool mixed (coniferous-deciduous) forests around 14.7 ka cal BP (Fig. 2; Suppl. Fig. 11). The species was abundant during the LGM and persisted during the subsequent period, as demonstrated by the relatively large number of Late Epigravettian sites in Hungary where reindeer bones are present.

In Southern Sweden, reindeer extinction took place at the transition from open pine-birch forest to pine-deciduous dominated forest transition. In the CB this coincidence of deciduous tree expansion and decline in the reindeer population cannot be demonstrated, although both elm (Ulmus) and hazel (Corylus) were already expanding regionally in the Great Hungarian Plain at 15.2 ka cal BP without a biome shift (Fig. 2). It is likely that rapidly increasing summer temperatures had a direct effect on the local reindeer population, and that this probably left the basin as a result of its metabolic adaptation to <15 °C. As the climate reconstructions used herein demonstrate, not only were the lowlands certainly too warm for reindeer by 15.2 ka cal BP, but also the mid mountain regions (Fig. 2). Another striking feature of the plain is that during the abrupt biome shift from a tundra to cool coniferous forest biome, the lowlands and low hills of the CB remained partially steppe covered, and this characteristic of the landscape must have helped the survival of grazers if their metabolic/physiological adaptation made it possible. The timing of the reindeer's withdrawal from the CB agrees well with the dates of the youngest/latest Epigravettian campsites⁸,⁶³. Even though the Epigravettian population also hunted horses by this time⁸, the coincidence of these two events suggests quarry fidelity and environmental determinism.

**Woolly mammoth.** Available summaries on European proboscidean extinction times suggest that woolly mammoth (*Mammuthus primigenius*) was present in the ice-free parts of Europe during the Weichselian Glacial until 14 ka cal BP, when its populations collapsed due to warming⁶⁴. Moreover, the endemic European mammoth population became extinct after 24 ka cal BP and was replaced by members of a Siberian mammoth genetic clade which had been colonizing Europe since 34 ka cal BP⁶⁵.

In the Austrian Alps, an area occupied by an extensive ice-stream network during the LGM, these animals migrated several tens of kilometres into alpine valleys during the first half of MIS 3⁶⁶, when ice-free conditions prevailed in the major valleys. Over 230 bones have been examined in Austria, and their distribution suggests that mammoth were present in river valleys and adjacent loess covered forelands of the Alps.

Considering the two periods of apparent absence of woolly mammoth in the CB (32.5–27.4 and 24.5–22 cal BP; Suppl. Fig. 8.), the onset of the latter period coincides with a massive dust accumulation period above Greenland (see Ca²⁺ on Fig. 3) followed by two short interstadials (GI-2.1 & GI-2.2). Regarding the earlier time interval, four Greenland interstadials fall within this (GI-3-4-5.1-5.2; Suppl. Fig. 8)³⁴, and it is known, mainly from loess mollusc studies, that the lowlands of the CB were covered by boreal parkland forests at these milder
time intervals⁶⁷. Furthermore, variations in glacial dust deposition on centennial–millennial timescales in the CB and Greenland were synchronous⁶⁸. Even though the number of ¹⁴C dated mammoth bones is still low, if the apparent low figures for the mammoth population or indeed its absence are valid, then it is likely that the warmer and more forested periods were disadvantageous to its populations in the CB. Moreover, the European mammoth population extinction after 24 ka cal BP likely also affected the CB, where on the basis of ¹⁴C measurements that form part of this study, the Siberian clade's expansion is probable after 22 ka cal BP.

**Epigravettian hunters and megafauna extinction in the CB: the relationship between human population and faunal change.** Even though human activity as a result of hunting and habitat modification are often cited as the principal driving force in megafauna extinction⁶⁹, the diversity of extinction patterns observed on different continents has led to an increasing recognition of the potential synergistic role of climate change⁶⁰. As demonstrated by Cooper et al.⁷⁰, in many cases the extinction of genetic clades coincided with rapid warming events in the absence of a human presence in North America. These events involved the rapid replacement of one species or population by a conspecific or congeneric one across a broad area. It appears that cold conditions were not an important driver of extinctions even in the presence of anatomically modern humans in Europe.

In the CB, Early Epigravettian groups were present during the LGM (GS-3–2.1c) and Late Epigravettian in the GS-2.1a-b and early GI-1 periods⁷. According to Lengyel et al.⁷¹ Early Epigravettian (26–20 ka cal BP) hunter-gatherers subsisted on reindeer and the wild horse, with reindeer predominating. A marked change was detected in the dominant quarry at the Late Epigravettian sites (20–15.2 ka cal BP), when reindeer decreased, while the wild horse became dominant, and mammoth was present again. These changes suggest that the LGM reindeer population thinned out in the CB during the last glacial termination. It has also been demonstrated that a decrease in human population of the CB took place at the end of the Late Epigravettian⁸. So far only Lovas (14–13 ka cal BP)⁷², and a stray find from Mezőlak (13.7–13.46 ka cal BP)⁷³,⁷⁴ are known from Transdanubia that are contemporaneous with the Epi-Magdalenian Culture of Czechia⁷⁵ and the Arched Backed Point techno complex of Poland⁷⁶.

These findings attest to the fact that human groups with new persistence strategies appeared in the western CB about 1000 years after the mammoth, reindeer and wild horse hunters left. Lovas provided evidence for the hunting of the Eurasian elk (A. alces) and Red deer (C. elaphus) and the use of their bones for ochre mining⁷⁷,⁷⁸. From these still fragmentary data it is possible to conclude that the Late Epigravettian population that seems to have left the CB around 15 ka cal BP had strong quarry fidelity and hunting habits, with the result that it followed the megafauna elements. Supporting this argument are the recently dated Late Palaeolithic camp sites (13–11.7 ka cal BP) further north in Slovakia, and a few in Southern Poland⁷⁹, where the fauna is too poorly preserved to determine hunted species.

## Conclusions

In this paper the hypothesis that transitions from cold to warm intervals were briefly optimal for grazing megafauna, followed by rapid extinctions has been tested¹⁵. The results support this hypothesis in that the dated bone assemblages of both woolly mammoth and reindeer attest to relatively large and increasing population sizes after the LGM until c. 16.2 ka cal BP. The paleoclimate and biome reconstructions suggest that a major warming at 16.2 ka cal BP was, however, detrimental to these populations, probably due to the intolerance to increased warmth in the case of reindeer. It has also been demonstrated that in the CB, vegetation regime shifts inferred from pollen data coincided with or slightly postdated the local population decline of woolly mammoth and reindeer. Furthermore, wildfires played a role in the transformation of the vegetation at the last glacial termination, when the array of biome changes was divergent. Afforestation by cold deciduous and cool coniferous trees was rapid in the mid mountains from 16.2 kyr cal BP, while steppe tundra biomes transformed into boreo-nemoral forest steppe in the lowland localities.

When the increasing number of well dated Late Epigravettian sites and their associated fauna in the CB are compared with the biome and climate reconstruction on the basis of pollen, it is possible to infer that the reason for the disappearance of the Epigravettian hunters was the diminishing population size of their dominant quarries, while vegetation change was more gradual. Placing this in a European context, the CB behaved similarly to some south and south-west European areas, and underwent earlier (post HE-1) and more drastic faunal and human population changes in response to early warming after the LGM than more north-westerly and northerly locations in Europe.

It was found that cave faunas in general can be used to estimate local extinction times. Direct AMS ¹⁴C measurements on the investigated species are required; large mammal bones from the same sediment layer often give deviating (older) ages to the associated micro-mammal fauna. Since Late Pleistocene cave faunas have extensively been studied and preserved in museum collections in the CB, there is still ample opportunity to resolve the many open questions of rapid faunal changes at the last glacial termination in this region of Europe.

## Methods

**Bone sampling from the museum collections.** The primary consideration in selecting remains was to choose large or medium sized mammal species or larger sized micromammals that were members of the vertebrate fauna in the Carpathian Basin during the Late Pleistocene but which became extinct in the latest Pleistocene or in the Holocene. Since the aim was to determine the probable age of disappearance, preference was given to bones of the same species from all layers in the sequence, where possible.
**Radiocarbon dating, age-depth modelling.** The selected bones were first subjected to physical preparation, including surface removal using a Dremel device, followed by grinding and sieving. Then 600 mg of the 0.5–1 mm fraction was placed into a test tube to perform acid–base-acid treatment (ABA), using 0.5 M HCl and 0.1 M NaOH reagent, rinsing the sample with ultrapure water in between. The pH of the samples was then adjusted to 3 and the test tubes were put in a block heater at 75 °C for 24 h to gelatinize the collagen. The liquid gelatine samples were then filtered using a 2 µm glass fibre filter and freeze dried for two days. For combustion, a ~4 mg gelatine sample together with MnO₂ reagent was placed to a glass combustion tube, which was, after sealing, heated to 550 °C⁸⁰,⁸¹. The CO₂ gas thus obtained was then purified and reduced to solid graphite, applying a customized sealed tube graphitization method⁸². The ¹⁴C measurements were performed using the EnvironMICADAS AMS instrument at the ICER laboratory⁸³.

¹⁴C ages were calibrated into calendar years using LibCal Rev. v. 8.2html software and the IntCal20 calibration curve³⁵. Age-depth modelling for the cave sediment sequences was performed with the Bacon package in R using Bayesian probability statistics⁸⁴. The woolly mammoth and reindeer radiocarbon dates presented in Supplementary Tables 4 and 5 were calibrated against the IntCal20 dataset using OxCal version 4.3 and incorporated within a single Phase model in OxCal version 4.3 in order to provide an estimate of the last appearance dates³²,⁸⁵. The 'Phase' command is a grouping model. It assumes no geographic relationship between samples, and that the ages represent a uniform distribution between a start and end boundary. The posterior distributions allowed the determination of probability distribution functions (PDFs) for the beginning and end of the phase. Modelled ages are reported here at the 95% probability range in thousands of calendar years BP (years; relative to AD 1950). An alternative model was also used to estimate extinction times of woolly mammoth and reindeer. The GRIWM model³³ inversely weights the contribution of each dated record to determine the terminal date. Equations and further details of the method are described in the Supplementary Material.

**Re-analysis and re-plotting of the mammal bone assemblages.** Bone assemblages were re-analysed for taxonomic inconsistencies by Piroska Pazonyi. Stratigraphic diagrams of the vole and other small mammal faunas were plotted in Psimpoll 2.27 and edited in Corel Draw 16.

**Quantitative climate reconstruction (pollen & vole thermometers).** The method of the pollen-based summer mean temperature ($T_{JJA}$) reconstruction from Kokad Mire was described in detail in Magyari et al.¹⁸ In brief, the European Modern Pollen Database (EMPD) was used as training set⁸⁶. 209 terrestrial pollen types were used for the summer mean temperature reconstruction, and this was further reduced to the dominant 43 taxa for the modern analogue technique reconstructions. The training set has 2687 modern pollen samples, of which 1240 are located below 600 m above sea level. Modern $T_{JJA}$ (mean air temperature for June, July and August) for each site was included in the EMPD climate data file. Quantitative pollen-based transfer functions were developed for $T_{JJA}$ using weighted averaging-partial least squares (WA-PLS) regression⁸⁷ with five components. Leave-one-out cross-validation⁸⁸ was employed to evaluate the model performance and estimated performance error statistics such as coefficient of determination ($r²$) between measured and predicted values, root mean square error of prediction (RMSEP) and maximum bias for each WA-PLS transfer function. The surface pollen data was transfomed to square-roots to reduce the noise of the data⁸⁹, a randomization t-test⁹⁰ was applied to select the most appropriate WA-PLS component for TJJA reconstructions, and sample-specific standard errors were calculated for the reconstructions, using a bootstrapping procedure with 1000 cycles⁹¹. The WA-PLS transfer functions and associated $T_{JJA}$ reconstructions were performed using the RIOJA package in R⁹². The statistical significance of the WA-PLS based $T_{JJA}$ reconstructions was assessed using an approach involving 999 randomizations developed by Telford and Birks⁹³, and these significance tests were conducted using the R package PALAEOSIG⁹³. An alternative transfer function, the modern analogue technique (MAT) was also used to derive $T_{JJA}$. This method does not require real calibration and is based on a comparison of past pollen assemblages to modern assemblages. The similarity between each fossil and modern pollen assemblage is evaluated using the chord distance metric⁹⁴,⁹⁵. In this study, the six modern pollen spectra that had the smallest distance were considered the best modern analogues of the given pollen spectrum, and were subsequently used for the reconstruction. If the chord distance was above a threshold defined by the Monte-Carlo method⁹⁴, the modern sample was considered a bad analogue, and was not taken into account in the reconstruction. Estimates of climatic parameters were obtained by taking a weighted average of the values for all selected best modern analogues, where the weights used were the inverse of the chord distance. MAT based $T_{JJA}$ reconstruction and statistical tests were run in R using the RIOJA package. In this study the focus was on the reconstruction of a single climatic parameter, $T_{JJA}$. The aim of the reconstruction is to compare the pollen-based $T_{JJA}$ record with the extinction times of both mega and microfauna elements and the pollen based TJJA reconstruction were also compared with the vole-thermometer based $T_{July}$ reconstruction from the Jankovich Cave and Rejtek Rock Shelter²¹,²³. The method is a paleozoogeographical calculation based on the principle of actualism; it is used to determine the July (summer) temperature of the accumulation period of different samples²¹,⁹⁶. The present-day optimum temperature for the distribution of certain vole-species is given in Jánossy and Kordos²³. These are 15 °C for *Clethrionomys*, 17.5 °C for *Arvicola*, 21 °C for *Microtus arvalis*, 19 °C for *Microtus agrestis*, 10 °C for *Lasiopodomys* (S.) *gregalis*, 12.5 °C for *Microtus oeconomus* and 7.5 °C for *Dicrostonyx*. After the multiplication of the values of mean July temperature by the percentage of the species in question compared to other vole-species, these values were summed up and divided by 100. *Microtus nivalis* was not taken into account because it is a mountain species, so its distribution is not zonal²¹,²³.

**Biome reconstruction.** For the pollen assemblage based biome reconstructions the technique of Prentice et al.⁹⁷ was used, as improved further by Tarasov et al.⁹⁸. The first step was the assignment of the pollen taxa to
plant functional types (PFT). The biome-PFT-taxon matrix published in Allen et al.⁹⁹ was used here. Plant functional types occupy specific bioclimatic spaces that can be assigned to one or several biomes. The biome-PFT matrix is a list of biomes, indicating which PFTs are characteristic of each biome. The next step is the calculation of biome affinity scores⁹⁷. This equation sums up the square roots of pollen percentages within a PFT and sums the affinity scores within a biome characterised by a set of PFTs. Eventually, the pollen sample is assigned to the biome with which it has the maximum affinity. The threshold value above which a taxon is considered is 0.5%. If the arboreal pollen sum is < 70%, the temperate forest biome is replaced by wooded steppe⁹⁹. In this paper, the stratigraphic plot of the main biome affinity scores and the assigned biomes is presented.

## Data availability

The data that support the findings of this study are openly available in Mendeley Data at http://dx.doi.org/10.17632/9kwbj5y54j.1.

Received: 3 August 2021; Accepted: 1 April 2022

Published online: 26 April 2022

## References

1. Vörös, I. Large mammal remains from the Upper Palaeolithic site at Esztergom-Gyurgyalag. *Acta Archaeol. Hung.* **43**, 261–263 (1991).

2. Jánossy, D. *Pleistocene Vertebrate Faunas of Hungary. Journal of Chemical Information and Modeling* (Akadémiai Kiadó, 1986).

3. Kordos, L. A sketch of the vertebrata biostratigraphy of the Hungarian Holocene. *Földrajzi Közlemények* **101**, 144–160 (1978).

4. Sümegi, P., Rudner, E. & Töröcsik, T. Environmental and chronological reconstruction problems during the Pleistocene/Holocene transition in Hungary (Magyarország pleisztocén végi és kora holocén környezeti változások kronológiai, tér és időbeli rekonstrukciós problémái). In *Öskoros Kutatók IV. Összejövetelének Konferenciakötete* (ed. Kolozsi, B.) 279–298 (Hajdú-Bihar Megyei Múzeumok Igazgatósága, 2012).

5. Bősken, J. et al. Investigating the last glacial Gravettian site 'Ságvár Lyukas Hill' (Hungary) and its paleoenvironmental and geo-chronological context using a multi-proxy approach. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **509**, 77–90 (2018).

6. Wilczyński, J. et al. Mammoth hunting strategies during the Late Gravettian in Central Europe as determined from case studies of Milovice I (Czech Republic) and Kraków Spadzista (Poland). *Quat. Sci. Rev.* **223**, 105919 (2019).

7. Lengyel, G. Reassessing the middle and late upper palaeolithic in Hungary. *Acta Archaeol. Carpathica* **51**, 47–66 (2016).

8. Béres, S. et al. Zöld cave and the late epigravettian in eastern central Europe. *Quat. Int.* **587–588**, 158–171 (2021).

9. Feurdean, A. et al. Trends in biomass burning in the Carpathian region over the last 15,000 years. *Quat. Sci. Rev.* **45**, 111–125 (2012).

10. Kuneš, P. et al. Interpretation of the last-glacial vegetation of eastern-central Europe using modern analogues from southern Siberia. *J. Biogeogr.* **35**, 2223–2236 (2008).

11. Pazonyi, P. Mammalian ecosystem dynamics in the Carpathian Basin during the last 27,000 years. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **212**, 295–314 (2004).

12. Sümegi, P. et al. Climatic fluctuations inferred for the middle and late plenigracial (MIS 2) based on high-resolution (~ca. 20 y) preliminary environmental magnetic investigation of the loess section of the Madaras brickyard (Hungary). *Cent. Eur. Geol.* **55**, 329–345 (2012).

13. Magyari, E. K. et al. Vegetation and environmental responses to climate forcing during the last glacial maximum and deglaciation in the East Carpathians: attenuated response to maximum cooling and increased biomass burning. *Quat. Sci. Rev.* **106**, 278–298 (2014).

14. Feurdean, A. et al. Climate variability and associated vegetation response throughout central and eastern Europe (CEE) between 60 and 8 ka. *Quat. Sci. Rev.* **106**, 206–224 (2014).

15. Mann, D. H. et al. Life and extinction of megafauna in the ice-age Arctic. *Proc. Natl. Acad. Sci. U. S. A.* **112**, 14301–14306 (2015).

16. Magyari, E. K. et al. Rapid vegetation response to Lateglacial and early Holocene climatic fluctuation in the South Carpathian Mountains (Romania). *Quat. Sci. Rev.* **35**, 116–130 (2012).

17. Magyari, E. K. et al. Late Plenigracial vegetation in eastern-central Europe: are there modern analogues in Siberia? *Quat. Sci. Rev.* **95**, 60–79 (2014).

18. Magyari, E. K. et al. Warm Younger Dryas summers and early late glacial spread of temperate deciduous trees in the Pannonian Basin during the last glacial termination (20–9 kyr cal BP). *Quat. Sci. Rev.* **225**, 105980 (2019).

19. Sümegi, P., Magyari, E., Dániel, P., Molnár, M. & Töröcsik, T. Responses of terrestrial ecosystems to Dansgaard-Oeschger cycles and Heinrich-events: a 28,000-year record of environmental changes from SE Hungary. *Quat. Int.* **293**, 34–50 (2013).

20. Hillebrand, J. *Paleolithic History (Az őskör története)* (Magyar Szemle Társaság, 1934).

21. Vértes, L., Kretzoi, M. & Herrmann, M. Neuere Forschungen in der Jankovich-Höhle. *Folia Archaeol.* **9**, 3–23 (1957).

22. Jánossy, D. Preliminary results of the paleontological investigations of a yet unknown rock shelter in the Bükk Mountains (A Bükk-hegység eddig ismeretlen kőfülkéjében végzett őslénytani ásatás előzetes eredménye, Répáshuta, Rejtek). *Karszt- és Barlkut. Tájékoztató* 72 (1963).

23. Jánossy, D. & Kordos, L. Pleistocene-Holocene Mollusc and Vertebrate Fauna of two caves in Hungary. *Ann. Hist. Musei Natl. Hungarici* **68**, 5–29 (1976).

24. Vértes, L. *Paleolithic and Mesolithic Remains in Hungary (Az őskör és az Átmeneti Kökor Emlékei Magyarországon)* (Akadémiai Kiadó, 1965).

25. Stieber, J. Oberpleistozäne Vegetationsgeschichte Ungarns im Spiegel anthrakotomischer Ergebnisse (bis 1957) (A magyarországi felsőpleisztocén vegetáció-története az anthrakotómiai eredmények (1957-ig) tükrében). *Földtani Közlöny* **97**, 305–317 (1967).

26. Jánossy, D. Vorläufige Ergebnisse der Ausgrabungen in der Felsnische Rejtek I. (Bükkgebirge, Gem. Répáshuta). *Karszt- és Barlangkutatás* 3, 49–58 (1961).

27. Kovács, J. Radiocarbon chronology of late Pleistocene large mammal faunas from the Pannonian basin (Hungary). *Bull. Geosci.* **87**, 13–19 (2012).

28. Willis, K. J., Braun, M., Sümegi, P. & Tóth, A. Does soil change cause vegetation change or vice versa? A temporal perspective from Hungary. *Ecology* **78**, 740–750 (1997).

29. Magyari, E. Holocene biogeography of *Fagus sylvatica* L. and *Carpinus betulus* L. in the Carpathian-Alpine Region. *Folia Hist. Musei Matra.* **26**, 15–35 (2002).

30. Magri, D. Persistence of tree taxa in Europe and quaternary climate changes. *Quat. Int.* **219**, 145–151 (2010).

31. Fükőh, L. Biostratigraphical investigation of the mollusc fauna of Rejtek I. rock-niche and Petényi Cave: Bükk Mountains, Hungary (Rejtek kőfülke és a Petényi-barlang (Bükk-hegység) Mollusca faunájának malakosztratigráfiai vizsgálata). *Folia Hist. Musei Matra.* **12**, 9–13 (1987).
32. Ramsey, C. B. & Lee, S. Recent and planned developments of the program OxCal. *Radiocarbon* **55**, 720–730 (2013).

33. Bradshaw, C. J. A., Cooper, A., Turney, C. S. M. & Brook, B. W. Robust estimates of extinction time in the geological record. *Quat. Sci. Rev.* **33**, 14–19 (2012).

34. Rasmussen, S. O. *et al*. A stratigraphic framework for abrupt climatic changes during the last glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. *Quat. Sci. Rev.* **106**, 14–28 (2014).

35. Reimer, P. J. *et al*. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). *Radiocarbon* **62**, 725–757 (2020).

36. Katona, L., Kovács, J., Kordos, L., Szappanos, B. & Linkai, I. The Csajág mammoths (*Mammuthus primigenius*): late Pleniglacial finds from Hungary and their chronological significance. *Quat. Int.* **255**, 130–138 (2012).

37. Buczkó, K. *et al*. Responses of diatoms to the Younger Dryas climatic reversal in a South Carpathian mountain lake (Romania). *J. Paleolimnol.* **48**, 417–431 (2012).

38. Tóth, M. *et al*. A chironomid-based reconstruction of late glacial summer temperatures in the southern Carpathians (Romania). *Quat. Res.* **77**, 122–131 (2012).

39. Sümegi, P. *et al*. Radiocarbon-dated paleoenvironmental changes on a lake and peat sediment sequence from the central Great Hungarian Plain (Central Europe) during the last 25,000 years. *Radiocarbon* **53**, 85–97 (2011).

40. Gill, J. L., Williams, J. W., Jackson, S. T., Lininger, K. B. & Robinson, G. S. Pleistocene megafaunal collapse, novel plant communities, and enhanced fire regimes in North America. *Science* **326**, 1100–1103 (2009).

41. Feurdean, A. *et al*. Fire hazard modulation by long-term dynamics in land cover and dominant forest type in eastern and central Europe. *Biogeosciences* **17**, 1213–1230 (2020).

42. Sümegi, P. *et al*. Radiocarbon dated complex paleoecological and geoarcheological analyses at the Bodrogkeresztúr–Henye Gravettian site (Ne Hungary). *Archeometriai Műhely* **13**, 31–41 (2016).

43. Herrmann, M., Jánossy, D., Steiber, J. & Vértes, L. Ausgrabungen in der Petényi- und Pesko-Höhle (Bükk-Gebirge). *Folia Archaeol.* **8**, 3–22 (1956).

44. Royer, A. How complex is the evolution of small mammal communities during the Late Glacial in southwest France?. *Quat. Int.* **414**, 23–33 (2016).

45. Crégut-Bonnouire, E. *et al*. The karst of the Vaucluse, an exceptional record for the last glacial maximum (LGM) and the Late-glacial period palaeoenvironment of southeastern France. *Quat. Int.* **339–340**, 41–61 (2014).

46. Cuenca-Bescós, G., Straus, L. G., González Morales, M. R. & García Pimienta, J. C. The reconstruction of past environments through small mammals: from the Mousterian to the Bronze Age in El Mirón Cave (Cantabria, Spain). *J. Archaeol. Sci.* **36**, 947–955 (2009).

47. Kovalchuk, O. *et al*. Living in a time of change: late Pleistocene/Holocene transitional vertebrate fauna of Grot Skeliastiy (Crimea, Ukraine). *Hist. Biol.* https://doi.org/10.1080/08912963.2020.1769094 (2020).

48. Puzachenko, A. Y. & Markova, A. K. Evolution of mammal species composition and species richness during the Late Pleistocene–Holocene transition in Europe: a general view at the regional scale. *Quat. Int.* **530–531**, 88–106 (2019).

49. Varga, Z. Extra-Mediterranean refugia, post-glacial vegetation history and area dynamics in Eastern Central Europe. In *Relict Species: Phylogeography and Conservation Biology* (eds Habel, J. C. & Assmann, T.) 57–87 (Springer Berlin Heidelberg, 2010).

50. Magyari, E. K. *et al*. Holocene persistence of wooded steppe in the Great Hungarian Plain. *J. Biogeogr.* **37**, 915–935 (2010).

51. Sommer, R. S. & Nadachowski, A. Glacial refugia of mammals in Europe: evidence from fossil records. *Mamm. Rev.* **36**, 251–265 (2006).

52. Mann, D. H., Groves, P., Gaglioti, B. V. & Shapiro, B. A. Climate-driven ecological stability as a globally shared cause of Late Quaternary megafaunal extinctions: the Plaids and Stripes Hypothesis. *Biol. Rev.* **94**, 328–352 (2019).

53. Lister, A. M. & Sher, A. V. Ice cores and mammoth extinction. *Nature* **378**, 23–24 (1995).

54. Owen-Smith, N. R. Megaherbivores: The Influence of Very Large Body Size on Ecology (Cambridge University Press, 1988).

55. Guthrie, R. D. *Frozen Fauna of the Mammoth Steppe: The story of Blue Babe* (The University of Chicago Press, 1990).

56. Huntley, B. *et al*. Millennial climatic fluctuations are key to the structure of last glacial ecosystems. *PLoS One* **8**, e61963 (2013).

57. Vörös, I. Large mammalian faunal changes during the Late Upper Pleistocene and Early Holocene times in the Carpathian Basin. In *Pleistocene Environment in Hungary* (ed. Pécsi, M.) 81–102 (Geographical Research Institute HAS, 1987).

58. Németh, A. *et al*. Holocene mammal extinctions in the Carpathian Basin: a review. *Mamm. Rev.* **47**, 38–52 (2017).

59. Merchant, R., Brewer, S., Webb, T. I. & Turvey, S. T. Holocenedeforestation: a history of human-environmental interactions, climate change, and extinction. In *Holocene Extinctions* (ed. Turvey, S. T.) 213–234 (Oxford University Press, 2009).

60. Lorenzen, E. D. *et al*. Species-specific responses of Late Quaternary megafauna to climate and humans. *Nature* **479**, 359–364 (2011).

61. Herre, W. Rangifer tarandus—Ren, Rentier. In *Handbuch der Saugetiere Europas 2/II Paarhufer—Artiodactyla* (eds Niethammer, J. & Krapp, F.) 198–216 (Aula Publisher, 1986).

62. Sommer, R. S., Kalbe, J., Ekström, J., Benecke, N. & Liljegren, R. Range dynamics of the reindeer in Europe during the last 25,000 years. *J. Biogeogr.* **41**, 298–306 (2014).

63. Lengyel, G. & Wilczyński, J. (2018) The Gravettian and the Epigravettian chronology in eastern central Europe: a comment on Bösken et al. (2017). *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **506**, 265–269 (2018).

64. Sommer, R. S. Late Pleistocene and Holocene history of mammals in Europe. *Handb. Mamm. Eur.* https://doi.org/10.1007/978-3-319-65038-8_3-1 (2020).

65. Palkopoulou, E. *et al*. Holarctic genetic structure and range dynamics in the woolly mammoth. *Proc. R. Soc. B Biol. Sci.* **280**, 20131910 (2013).

66. Spötl, C., Reimer, P. J. & Göhlich, U. B. Mammoths inside the Alps during the last glacial period: radiocarbon constraints from Austria and palaeoenvironmental implications. *Quat. Sci. Rev.* **190**, 11–19 (2018).

67. Sümegi, P. *Loess and Upper Paleolithic Environment in Hungary: An Introduction to the Environmental History of Hungary* (Aurea, 2005).

68. Újvári, G. *et al*. Coupled European and Greenland last glacial dust activity driven by North Atlantic climate. *Proc. Natl. Acad. Sci. U. S. A.* **114**, E10632–E10638 (2017).

69. Haynes, G. Extinctions in North America's late glacial landscapes. *Quat. Int.* **285**, 89–98 (2013).

70. Cooper, A. *et al*. Abrupt warming events drove Late Pleistocene Holarctic megafaunal turnover. *Science* **349**, 602–606 (2015).

71. Lengyel, G. *et al*. The Epigravettian chronology and the human population of eastern Central Europe during MIS2. *Quat. Sci. Rev.* **271**, 1071–1078 (2012).

72. Sajó, I. E. *et al*. Core-shell processing of natural pigment: upper Palaeolithic red ochre from Lovas, Hungary. *PLoS One* **10**, e1–18 (2015).

73. Horváth, T. & Ilon, G. Mezőlak-Szélmező-peatbog: an unusual prehistoric site (Mezőlak-szélmező-tözegtelep: egy nem hétköznapi őskori lelőhely). *Archeometriai Műhely* **14**, 143–183 (2017).

74. Zalai-Gaál, I. Possibilites of the social-archaeological studies of the Neolithic. *Antaeus* **27**, 449–471 (2004).

75. Reade, H. *et al*. Magdalenian and Epimagdalenian chronology and palaeoenvironments at Kúlna Cave, Moravia, Czech Republic. *Archaeol. Anthropol. Sci.* https://doi.org/10.1007/s12520-020-01254-4 (2021).

76. Łanczont, M. *et al*. Late Glacial environment and human settlement of the Central Western Carpathians: a case study of the Nowa Biała 1 open-air site (Podhale Region, southern Poland). *Quat. Int.* **512**, 113–132 (2019).

77. Mészáros, G. & Vértes, L. A paint mine from the early Upper Palaeolithic age near Lovas (Hungary, county Veszprém). *Acta Archaeol. Acad. Sci. Hung.* **5**, 5–34 (1955).
78. Pathou-Mathis, M. Nouvelle analyse du métérial osseux du site de Lovas. *Praehistoria* **3**, 161–175 (2002).

79. Sobkowiak-Tabaka, I. & Diachenko, A. Approaching daily life at Late Palaeolithic camps: the case of Lubrza 10, Western Poland. *Prahistorische Z.* **95**, 311–333 (2020).

80. Molnár, M. et al. EnvironMICADAS: a mini 14C AMS with enhanced gas ion source. *Radiocarbon* **55**, 338–344 (2013).

81. Major, I. et al. Assessment and development of bone preparation for radiocarbon dating at HEKAL. *Radiocarbon* **61**, 1551–1561 (2019).

82. Rinyu, L. et al. Optimization of sealed tube graphitization method for environmental C-14 studies using MICADAS. *Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater.* **294**, 270–275 (2013).

83. Molnár, M. et al. Status report of the new AMS 14 C sample preparation lab of the Hertelendi laboratory of environmental studies (Debrecen, Hungary). *Radiocarbon* **55**, 665–676 (2013).

84. Blaauw, M. & Christeny, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. *Bayesian Anal.* **6**, 457–474 (2011).

85. Kosintsev, P. et al. Evolution and extinction of the giant rhinoceros Elasmotherium sibiricum sheds light on late Quaternary megafaunal extinctions. *Nat. Ecol. Evol.* **3**, 31–38 (2019).

86. Davis, B. A. S. et al. The European modern pollen database (EMPD) project. *Veg. Hist. Archaeobot.* **22**, 521–530 (2013).

87. ter Braak, C. J. F. & Juggins, S. Weighted averaging partial least squares regression (WA-PLS): an improved method for reconstructing environmental variables from species assemblages. *Hydrobiologia* **269–270**, 485–502 (1993).

88. Birks, H. J. B., Line, J. M., Juggings, S., Stevenson, A. C. & ter Braak, C. J. F. Diatoms and pH reconstruction. *Philos. Trans. R. Soc. B* **327**, 263–278 (1990).

89. Prentice, I. C. Multidimensional scaling as a research tool in quaternary palnology: a review of theory and methods. *Rev. Palaeobot. Palynol.* **31**, 71–104 (1980).

90. van der Voet, H. Comparing the predictive accuracy of models using a simple randomization test. *Chemom. Intell. Lab. Syst.* **25**, 313–323 (1994).

91. Birks, H. J. B. Quantitative palaeoenvironmental reconstructions from holocene biological data. *Glob. Change Holocene* https://doi.org/10.4324/9780203785027 (2003).

92. Rioja, J. S. Analysis of Quaternary Science Data, R package version (0.8-5). (2012).

93. Telford, R. J. & Birks, H. J. B. A novel method for assessing the statistical significance of quantitative reconstructions inferred from biotic assemblages. *Quat. Sci. Rev.* **30**, 1272–1278 (2011).

94. Guiot, J. Methodology of the last climatic cycle reconstruction in France from pollen data. *Palaeogeogr. Palaeoclimatol. Palaeoecol.* **80**, 49–69 (1990).

95. Birks, H. J. B. Ecological palaeoecology and conservation biology: controversies, challenges, and compromises. *Int. J. Biodivers. Sci. Ecosyst. Serv. Manag.* **8**, 292–304 (2012).

96. Kordos, L. Climatostratigraphy of Upper Pleistocene vertebrates and the conditions of loess formation in Hungary. *GeoJournal* **15**, 163–166 (1987).

97. Prentice, I. C., Guiot, J., Huntley, B., Jolly, D. & Cheddadi, R. Reconstructing biomes from palaeoecological data: a general method and its application to European pollen data at 0 and 6 ka. *Clim. Dyn.* **12**, 185–194 (1996).

98. Tarasov, P. E. et al. Present-day and mid-Holocene biomes reconstructed from pollen and plant macrofossil data from the former Soviet Union and Mongolia. *J. Biogeogr.* **25**, 1029–1053 (1998).

99. Allen, J. R. M., Watts, W. A. & Huntley, B. Weichselian palnostratigraphy, palaeovegetation and palaeoenvironment; the record from Lago Grande di Monticchio, southern Italy. *Quat. Int.* **73–74**, 91–110 (2000).

## Acknowledgements

This research was supported by the European Union and the State of Hungary, co-financed by the European Regional Development Fund in the project of GINOP-2.3.2.-15-2016-00009 'ICER' (Dating megafauna extinction in the Carpathian basin and examining its relation to paleoenvironmental and paleoclimate changes). EKM acknowledges the support of OTKA Research Fund (K129167). This is MTA-MTM Paleo Contribution No. 349. Hereby we declare that we obtained permission from the Hungarian Natural History Museum and the Mining and Geological Survey of Hungary for the subsampling of bone specimens for ¹⁴C measurements.

## Author contributions

E.K.M. wrote the manuscript, was responsible for the pollen analyses, for the conceptualization, methodology, supervision, project administration, funding acquisition, review & editing; Z.S.Z. and J.K. took part in the visualisation and statistical analyses of the data; P.P. and M.G. analysed the bone assemblages, wrote part of the manuscript, took part in the conceptualization; I.M. was responsible for the stable and radioactive isotope analyses, wrote part of the original draft, reviewed & edited the manuscript; I.P. analysed pollen samples and took part in data visualisation; G.Y.L. was responsible for archaeological data acquisition, reviewed & edited the manuscript; A.V. took part in the sample collection, reviewed & edited the manuscript; A.H. did the statistical analysis of the macrocharcoal records and plotted diagrams.

## Funding

Open access funding provided by Eötvös Loránd University.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-10714-x.

**Correspondence** and requests for materials should be addressed to E.K.M.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022