Next Article in Journal
Acknowledgment to Reviewers of Quaternary in 2020
Next Article in Special Issue
Prey Size Decline as a Unifying Ecological Selecting Agent in Pleistocene Human Evolution
Previous Article in Journal
Quaternary Permafrost in China: Framework and Discussions
Previous Article in Special Issue
Magnetostratigraphy and Chronology of the Lower Pleistocene Primate Bearing Dafnero Fossil Site, N. Greece
 
 
Correction published on 8 April 2022, see Quaternary 2022, 5(2), 21.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dating of the Lower Pleistocene Vertebrate Site of Tsiotra Vryssi (Mygdonia Basin, Greece): Biochronology, Magnetostratigraphy, and Cosmogenic Radionuclides

1
Palaeoanthropology, Senckenberg Centre for Human Evolution and Palaeoenvironment, Eberhard Karls University of Tübingen, Rümelinstr. 23, 72070 Tübingen, Germany
2
Laboratory of Geology and Palaeontology, School of Geology, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
3
Dipartimento di Scienze della Terra “Ardito Desio”, University of Milan, via Mangiagalli 34, I-20133 Milan, Italy
4
Department of Geosciences, Eberhard Karls University of Tübingen, Schnarrenbergstr. 94-96, 72076 Tübingen, Germany
5
Department of Geophysics, School of Geology, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Submission received: 11 November 2020 / Revised: 30 December 2020 / Accepted: 5 January 2021 / Published: 8 January 2021 / Corrected: 8 April 2022

Abstract

:
Background and scope: The late Villafranchian large mammal age (~2.0–1.2 Ma) of the Early Pleistocene is a crucial interval of time for mammal/hominin migrations and faunal turnovers in western Eurasia. However, an accurate chronological framework for the Balkans and adjacent territories is still missing, preventing pan-European biogeographic correlations and schemes. In this article, we report the first detailed chronological scheme for the late Villafranchian of southeastern Europe through a comprehensive and multidisciplinary dating approach (biochronology, magnetostratigraphy, and cosmogenic radionuclides) of the recently discovered Lower Pleistocene vertebrate site Tsiotra Vryssi (TSR) in the Mygdonia Basin, Greece. Results: The minimum burial ages (1.88 ± 0.16 Ma, 2.10 ± 0.18 Ma, and 1.98 ± 0.18 Ma) provided by the method of cosmogenic radionuclides indicate that the normal magnetic polarity identified below the fossiliferous layer correlates to the Olduvai subchron (1.95–1.78 Ma; C2n). Therefore, an age younger than 1.78 Ma is indicated for the fossiliferous layer, which was deposited during reverse polarity chron C1r. These results are in agreement with the biochronological data, which further point to an upper age limit at ~1.5 Ma. Overall, an age between 1.78 and ~1.5 Ma (i.e., within the first part of the late Villafranchian) is proposed for the TSR fauna. Conclusions: Our results not only provide age constraints for the local mammal faunal succession, thus allowing for a better understanding of faunal changes within the same sedimentary basin, but also contribute to improving correlations on a broader scale, leading to more accurate biogeographic, palaeoecological, and taphonomic interpretations.

1. Introduction

The fascinating world of the Pleistocene has long attracted the interest of scientists from different disciplines. The rich mammal record, including iconic animals such as mammoths, saber-toothed cats, giant hyenas, and giant deer; the evolution and dispersals of our own genus Homo; and the large-scale global climatic fluctuations are only some of the highly intriguing research areas that characterize the Pleistocene. The late Villafranchian large mammal age, defined in Europe as starting at ~2.0 Ma [1,2,3] and ending at ~1.2 Ma [4,5,6] at the beginning of the subsequent Epivillafranchian, is a crucial period in the evolution of the Early Pleistocene terrestrial ecosystems of Eurasia. The climatic conditions during the late Villafranchian resulted in the (intermittent) expansion of open-forest and savanna-like landscapes from the Transcaucasian to the European peri-Mediterranean regions [7,8]. In turn, this affected the structure and dynamics of the western Eurasian mammal palaeocommunities tremendously, triggering faunal turnovers, extinctions, and dispersal events [2,7,9], including the first arrival of hominins at the gates of Europe [10,11]. Indeed, the fossil record of the northern/eastern peri-Black Sea region and southern Europe documents the arrival during this period of new immigrants of Asian and African origin along with the regional extinction of several, previously well-adapted lineages [1,2,7,12]. The analysis and interpretation of these mammal bioevents and their underlying mechanisms necessitate the establishment of an accurate chronological framework. The discovery of several new fossiliferous localities and the re-assessment of older ones, in combination with improved dating methodologies, have substantially enhanced our knowledge about the Early Pleistocene of Europe. However, there are still several open issues that need to be addressed. For instance, the dating methods (either direct or indirect) that are commonly applied to Quaternary deposits and fossils, have their own assumptions, as well as assets and limitations (see [13] for an overview) that bring uncertainties in age estimations, thus, limiting high-resolution correlations across different sites. Moreover, although Early Pleistocene local European records and faunal successions have been relatively well explored and defined through time, pan-European biogeographic correlations and schemes (large-scale biochronology) are still limited (see also [1]) and poorly supported by an accurate chronological framework, especially for the Balkans and adjacent territories. As a result, our understanding of the palaeoecological conditions, including corridors or barriers that may have enabled or impeded the dispersal of several large mammals in this timeframe, is highly obscured. Additionally, these continental-scale biogeographic extrapolations rely mostly on spot data from sedimentary basins from distant geographic regions, such as Italy, Spain, France, and the Balkans (Figure 1), which, as stated above, are difficult to correlate, resulting in the following uncertainties especially when attempting to resolve: (i) possible discontinuities in the fossil record, both local or composite (faunal units); (ii) the degree of diachrony or synchrony of the various dispersal events in response, for example, to Pleistocene climate variability; and (iii) the different geological, climatic, environmental, and taphonomic conditions at the various sites (see also [14,15]).
The Lower Pleistocene fossiliferous localities of the Mygdonia Basin (central Macedonia, Greece) have been known since the end of the 1970s [16,17,18]. Numerous fieldwork campaigns, carried out for several decades, have led to the discovery of several sites, and to the excavation and study of a great number of fossils. This work has significantly enriched the Pleistocene fossil record of Greece [19,20,21,22]. Currently, there are twelve open-air Lower Pleistocene vertebrate sites known in the Mygdonia Basin (Figure 2), located in at least three distinct stratigraphic horizons and correlated to different parts of the late Villafranchian. Owing to their rich and diversified fauna of different evolutionary stages, the Mygdonia sites (1) permit identification and tracking of faunal and evolutionary changes within the stratigraphic succession of the same sedimentary basin, and thus enhance our knowledge on the evolution and taxonomy of several vertebrate taxa and mammal turnovers; (2) allow a more precise biostratigraphic correlation with other fossil assemblages of Greece and southeastern Europe; and (3) contribute to palaeogeographical interpretations and the better understanding of the palaeoenvironmental conditions, offering overall valuable insights into the Early Pleistocene terrestrial ecosystems of Europe. Moreover, Greece and the Balkans in general, located in the southeastern corner of Europe, occupy a critical position for the dispersals of mammals and hominins to and from Europe, western Asia, and the Levant/Africa [22,23,24,25,26,27], highlighting the need for systematic and deeper investigations in key locations, such as the Mygdonia Basin. Importantly for the Pleistocene continental biostratigraphy of Greece, the late Villafranchian localities of the Mygdonia Basin, together with the older ones from the Aliakmon Basin (e.g., Dafnero) and the younger ones from the Megalopolis Basin (e.g., Marathousa-1 and Kyparissia sites), all of which are characterized by long stratigraphic sequences, altogether cover a rather continuous faunal sequence spanning from the middle Villafranchian to the Galerian [22,28,29,30,31], and permit the detailed study of the Early–Middle Pleistocene fossil and palaeoenvironmental record of Greece.
In this article, we report the first detailed chronological scheme for the late Villafranchian from southeastern Europe. This has been achieved through a comprehensive and multidisciplinary dating approach (biochronology, magnetostratigraphy, and cosmogenic radionuclides) of the recently discovered Lower Pleistocene vertebrate site Tsiotra Vryssi (TSR) in the Mygdonia Basin. More specifically, the aims of the present study are the following: (1) to precisely constrain the age of the site by means of a combination of dating methods; (2) to follow more thoroughly the regional mammal succession through time and to enhance our knowledge about the relationships among the various Mygdonia Basin localities; and (3) to correlate TSR with other Lower Pleistocene localities of western Eurasia, thus, placing TSR within the larger discussions on Early Pleistocene mammal diversity and bioevents. As one of the first such studies from the Early Pleistocene of Greece and Southeastern Europe in general, and the first from the Mygdonia Basin, this study aspires to render TSR and the surrounding fossil sites to be a western Eurasian reference region for the chronological correlation of other localities and for the study and comparison of their vertebrate assemblages.

2. Geological Setting and Fossiliferous Locality

The Mygdonia Basin is located in central Macedonia (Greece), to the northeast of Thessaloniki, and represents an elongated east-west trending tectonic depression, the formation of which started during the early–middle Miocene with the incipient development of the Pre-Mygdonia Basin (Figure 2). During the Neogene–Εarly Pleistocene, the basin was filled with fluvial/fluvio-torrential and lacustrine sediments, while at the beginning of the Μiddle Pleistocene a new tectonic event led to further subdivisions into smaller basins, which were subsequently filled mostly by lacustrine sediments [32]. These Neogene to Quaternary deposits were subdivided into two high-rank lithostratigraphic units, namely the Pre-Mygdonian Group (Neogene to Εarly Pleistocene) and the Mygdonian Group (Middle Pleistocene to Holocene). The Pre-Mygdonian Group can be further subdivided into three successive lithostratigraphic formations (Fm), the Chryssavgi Fm, the Gerakarou Fm and, the Platanochori Fm [19]. The following six Lower Pleistocene fossiliferous sites have been discovered, to date, in the upper part of the Gerakarou Fm: Gerakarou-1 (GER), Vassiloudi (VSL), Krimni-1 and -2 (KRI and KRM), Kalamoto-2 (KLT), and Tsiotra Vryssi (TSR). The study of their assemblages indicated that, although the localities are not isochronous, all of them can be correlated to the first part of the late Villafranchian (Early Pleistocene) [19,20,21]. The overlying Platanochori Fm includes six more fossiliferous sites, i.e., Apollonia-1 (APL), Ravin of Voulgarakis (RVL), Riza-1 (RIZ), Marathoussa (MAR), Kalamoto-1 (KAL), and Platanochori-1 (PLN) (Figure 2), which are correlated to the second part of the late Villafranchian [19,20,21]. More information on the stratigraphy and the fossiliferous sites can be found in [19] and [21].
TSR (GPS coordinates 40°31′23.2″ N, 23°26′07.5″ E and elevation 250 m a.s.l.) is situated in the Gerakarou Fm, the most widely exposed sedimentary unit across the Mygdonia Basin (Figure 2) with a thickness in excess of 100 m. The Gerakarou Fm consists of red-brown sands, gravels, sandy silts, and clays, deposited in a fluvial/fluvio-torrential environment [19]. TSR was discovered in 2014 during a targeted survey by a joint team from the Aristotle University of Thessaloniki and the Paleoanthropology working group at the Eberhard Karls University of Tübingen [21]. It has since yielded a wealth of large mammal fossils (Figure 3), as well as some micromammals, reptiles, and birds. The ~1 m thick fossiliferous unit Geo 2a (Figure 4) consists of pale brown-dark yellowish brown poorly sorted silts, locally intercalated by cm-thick lenses of medium-coarse grained sands and is relatively more clayey in its uppermost part. The spatial taphonomic analysis of the faunal remains and the sedimentological and micromorphological studies indicated multiple dispersion events and recurrent spatial rearrangement of a lag, (peri)autochthonous assemblage, consistent with the cyclical lateral switching of a braided fluvial system [33]. The large mammal assemblage is rich and diversified and, so far, comprises 14 species ([21,34,35,36] and unpublished data; a cranium belonging to the felid Megantereon was discovered in the excavation season of 2018, whose complete study is pending, however, the genus is included in the herein updated faunal list). The fauna is dominated by horses (Equus), but also includes several other herbivore taxa, such as bovids (Leptobos, Bison, Pontoceros), deer (Cervus, Praemegaceros), rhinos (Stephanorhinus), mammoths (Mammuthus), and giraffids (Palaeotragus). Carnivores are well represented and include wolf-like dogs (Canis), giant hyenas (Pachycrocuta), bears (Ursus), and saber-toothed cats (Megantereon) (Table 1). TSR has yielded more than 2000 mammalian bones and teeth to date, which are stored at the Museum of Geology, Palaeontology, and Palaeoanthropology of the Laboratory of Geology and Palaeontology at the Aristotle University of Thessaloniki, Greece (LGPUT).

3. Methodology

Cosmogenic radionuclides: Pairs of in situ-produced cosmogenic radionuclides analyzed in fluvial sediment can be used for burial dating based on different decay rates of the two isotopes (e.g., [37,38]). Age determination of sediment deposition in the study area makes use of simple burial dating with in situ-produced 26Al and 10Be concentrations in quartz (e.g., [39]). Three sand samples from the same stratigraphic depth located at ~1 m below the bottom of the Geo 2b fossiliferous layer (Figure 4a) were sieved to 0.25 to 0.50 mm. Quartz was separated and cleaned to reach Al concentrations of ~100 ppm. Al and Be were extracted according to protocols outlined in [40]. The 26Al/27Al and 10Be/9Be ratios of samples, blanks, and standards were measured at the AMS facility of Cologne (Germany). The 26Al and 10Be concentrations of samples were corrected for chemistry blanks of 1.93 × 105 and 4.97 × 104 atoms, respectively. Standards used for calibration of the 26Al/27Al and 10Be/9Be ratios were KN01-5-3 and KN01-4-3 [41], as well as KN01-6-2 and KN01-5-3 [42], respectively. All sample preparation was conducted at the University of Tübingen.
Simple burial dating was performed with a MatLab script modified from [40] based on [39]. Decay constants used for 26Al and 10Be were (9.830 ± 0.250) × 10−7 [43] and (4.997 ± 0.043) × 10−7 [44,45]. The sea level and high latitude (SLHL) production rates for 26Al were 28.54, 0.84, and 0.081 atoms/(g(qtz) yr) for nucleonic, stopped muonic, and fast muonic production, respectively [46,47]. The SLHL production rates used were 3.92, 0.012, and 0.039 atoms/(g(qtz) yr) for nucleonic, stopped muonic, and fast muonic 10Be production, respectively [46,47]. The 26Al/10Be ratio of production at SLHL was ~7.4. The SLHL production rates were scaled to the sample location with the online calculator of Marrero et al. [48] using the scaling procedure “SA” based on Lifton et al. [49]. Production rates at individual sampling depths were calculated based on nucleonic, stopped muonic, and fast muonic adsorption lengths, which were 157, 1500, and 4320 g/cm2, respectively [46]. Overburden was assumed to be at least 3000 cm with a constant density over time of 2.0 ± 0.1 g/cm3. This overburden results in a calculated simple burial age, which can be considered to be a minimum burial age [39]. The production rate at the sample location was used as the production rate of the sediment source area, because the production rate of the palaeo-sediment source area was unknown. Therefore, this production rate is considered to be a minimum production rate for a source area located at higher elevation, and hence results in a maximum calculated burial age. Thus, based on the above, the calculated simple burial age has the potential to be either younger or older than the actual age, and we interpreted it to be the best estimate of the actual burial age given with available constraints.
Palaeomagnetism: A total of 52 standard (10 cc) samples were recovered from 4.5 stratigraphic meters of the Tsiotra Vryssi section (Figure 4a), either by manually inserting plastic boxes in the outcrop section (32 specimens) or by hand-sampling six oriented blocks, each yielding up to 5 specimens, which were subsequently consolidated and cored in the laboratory. All analyses were conducted at the Alpine Laboratory of Paleomagnetism (ALP) of Peveragno (Italy), with the exception of the thermal demagnetization experiments, which were performed at the Department of Geophysics of the Aristotle University of Thessaloniki.
The low-field magnetic susceptibility (κ) was measured with an AGICO KLY-3 Kappabridge. Then, samples were AF (Alternating Field) demagnetized in steps of 5–10 mT up to 100 mT with a 2G Enterprises AF Degausser. The natural remanent magnetization (NRM) was measured after each demagnetization step with a 2G Enterprises 755 DC-SQUID cryogenic magnetometer located in a magnetically shielded room. The directions of the NRM were plotted on standard vector end-point demagnetization diagrams and the characteristic remanent magnetization (ChRM), where present, was isolated with standard principal component analysis on selected data. Isothermal remanent magnetization (IRM) acquisition experiments were performed on 6 representative samples using an ASC Scientific IM-10-30 impulse magnetizer and an AGICO JR-6 spinner magnetometer. Finally, thermal demagnetization experiments of a three-component IRM adopting 2 T, 0.5 T, and 0.1 T orthogonal fields were performed on another set of 4 samples with a Minispin spinner magnetometer (Molspin) at the Department of Geophysics of the Aristotle University of Thessaloniki.
Faunal similarity analysis: In order to examine the similarities and relationships of TSR with other Early Pleistocene localities, we performed hierarchical cluster analysis at the genus level. Twenty-four faunas/sites were selected ranging geographically from Georgia to Spain and spanning chronologically from ~2.5 to ~0.8 Ma (middle Villafranchian to Epivillafranchian). A total of 56 genera of large mammals (excluding primates) were considered. The analysis was conducted with the PAST software, version 4.04 [50]. The Q-mode dendrogram (with unweighted pair-group average (UPGMA) algorithm) was acquired by applying the similarity index Raup-Crick, which is used to statistically test presence and absence data through a randomization method [51]. The faunal assemblages considered are Fonelas P-1 [52], Untermassfeld [53,54], Gerakarou-1, Krimni-1, Apollonia-1, Libakos [55], Slivnitsa [56], Poggio Rosso [57], Kalamoto-2 [20], Dmanisi [11], Pirro 10 [58], Sima del Elefante (TE-) Lower Red Unit (LRU) TE8–14 [59], Vallparadís Estació Middle Unit, Layers EVT6–VT8 [60,61], Venta Micena [62], Val di Chiana [9], Pietrafitta [9], Trlica 11–10 [63], Vallonnet [64], Fuente Nueva 3 [62], Barranco León [62], Pantalla [65], Dafnero [66] and Taurida [67], with faunal updates and taxonomical reviews by us for the Greek sites.

4. Results

4.1. Cosmogenic Radionuclides

In situ produced cosmogenic radionuclide concentrations of 26Al and 10Be from three samples (Figure 4a) range from (31.88 ± 1.74) to (37.07 ± 1.91) × 104 atoms/g(qtz) and (9.40 ± 0.33) to (10.51 ± 0.42) × 104 atoms/g(qtz). The 26Al/10Be ratios are 3.53 ± 0.23, 3.18 ± 0.20, and 3.39 ± 0.22 for the three samples from the same stratigraphic layer. These ratios result in simple burial ages of 1.88 ± 0.16 Ma, 2.10 ± 0.18 Ma, and 1.98 ± 0.18 Ma, respectively (Table 2).

4.2. Magnetic Properties

Values of κ are relatively low all along the section (Figure 4c), with an average value of 105 × 10−6 SI. Samples from the basal part of the sequence, between 0 and 2 m from the base, show the lowest values (15–60 × 10−6 SI), whereas peak κ values were measured from samples collected at 2 m (ca. 150 × 10−6 SI), 3.5 m (ca. 215 × 10−6 SI), and 4.5 m (ca. 275 × 10−6 SI). The IRM acquisition curves indicate the dominance of low to medium coercivity phases (80 to 400 mT) with saturation reached at 2500 mT, suggesting the presence of a mixture of magnetite/maghemite and hematite (Figure 5A). The three-axis IRM experiments [68] show that the medium (0.5 T) and high-coercivity (2 T) curves have maximum unblocking temperatures of 650–675 °C (Figure 5B), consistent with the presence of hematite. The low-coercivity (0.1 T) curve declines at a rather constant pace from room temperature up to 350–400 °C, before recording a rapid decrease at ~600 °C, followed by a slower decrease up to 650–675 °C (Figure 5B); this behavior could be interpreted to reflect progressive oxidation of an original maghemite/magnetite phase into a secondary hematite.

4.3. Magnetostratigraphy

Bipolar ChRM component directions, oriented north-and-down or south-and-up in in situ coordinates and within a field range of ~10–70 mT, were isolated in 23 of 38 demagnetized samples from the Tsiotra Vryssi section (Figure 6). The mean ChRM direction, calculated by applying standard Fisher statistics on n = 23 ChRM directions (declination = 176.8 °E, inclination = −40.5°, α95 = 10.7°, Table 3 and Figure 7), led to a palaeomagnetic pole position of latitude 72.4 °N and longitude 213.2 °E (A95 = 7.8°). The McFadden and McElhinny [69] reversal test applied to the ChRM component directions in in situ coordinates result was indeterminate (k1 ≠ k2, type I), probably due to the small number of components used.
When plotted versus stratigraphic depth, the latitude (positive for normal polarity, negative for reverse polarity) of virtual geomagnetic poles (VGPs) calculated for each ChRM component direction (Figure 4) define a lower incomplete normal magnetic polarity interval overlain by an incomplete reverse polarity interval (Figure 4).

4.4. Biochronology

The presence of the large-sized hyena Pachycrocuta brevirostris in the TSR fauna indicates a post-middle Villafranchian age, as the migration of this species to Europe, referred to as the “Pachycrocuta brevirostris event”, is tentatively dated at ~2.0 Ma [3] and marks the beginning of the late Villafranchian [2]. Approximately at the same time, Canis sensu stricto also becomes widespread [70,71]. This taxon’s earliest occurrences in the Balkan Peninsula are traced at Slivnitsa (Bulgaria) and Gerakarou 1, correlated to the middle/late Villafranchian boundary [24,72,73]. The bear Ursus etruscus from TSR differs from the middle Villafranchian specimens of Saint-Vallier (France) and is more similar to the typical U. etruscus from Upper Valdarno (Olivola and Tasso Faunal Units (FUs) of Italy); however, it is most similar to the specimens from Pietrafitta (Farneta FU) [36], tentatively dated at ca. 1.7–1.5 Ma [1,3,74]. The large-sized deer Praemegaceros and the spiral-horned antelope Pontoceros first occur at Dmanisi and are present soon after in the northern peri-Black Sea region and the Balkan Peninsula [75,76], while Praemegaceros appears in Italy during the Farneta FU [1] and in Spain after ~1.6 Ma, at the Orce sites [77]. The giraffid Palaeotragus is well recorded during the middle Villafranchian, and is also present in the late Villafranchian of Fonelas P-1, Dmanisi, Libakos, and Denizli-Kocabaş [78], but has not been documented, as yet, in the fossil-rich Pirro Nord, the Orce sites, Apollonia 1, or Untermassfeld, collectively correlated to the second half of the late Villafranchian and Epivillafranchian [1,7]. The strongest biochronological indicator recovered so far from TSR is the co-existence of the large-sized bovids Leptobos and Bison. Leptobos was present in Europe during the early and middle Villafranchian but also survived during the first part of the late Villafranchian (e.g., [65,79]). The earliest record of bisons in western Eurasia is documented at Dmanisi (~1.8 Ma), while they subsequently dispersed into peri-Mediterranean Europe from ~1.7–1.5 Ma onwards (e.g., Farneta FU, Orce sites, Pirro Nord, Apollonia-1, Sima del Elefante-TE9). The Leptobos-Bison co-existence is relatively rare; apart from TSR, it is recorded in Mygdonia Basin at Krimni-1 and Kalamoto-2, as well as further afield at Trlica 11–10 in Montenegro, the Taurida Cave in Russia and Pietrafitta in Italy [35,67,75,79,80].
In order to examine more comprehensively the overall faunal relationships and similarities between TSR and other middle Villafranchian–Epivillafranchian localities of western Eurasia, we conducted hierarchical cluster analysis at the genus level. The resulting Q-mode dendrogram (Figure 8) clusters the faunas chronologically (rather than geographically). Two main clusters emerge, reflecting the significant faunal reorganization (even at genus level) close to ~2.0 Ma at the middle/late Villafranchian boundary, i.e., cluster A, including the middle and middle/late Villafranchian faunas (Dafnero and Slivnitsa) and cluster B encompassing the late Villafranchian to Epivillafranchian ones. Cluster B is further divided into two subclusters, i.e., B1 consisting of the late Villafranchian to Epivillafranchian faunas younger than ~2.0–1.9 Ma and B2 including those faunas that are at or slightly after the middle/late Villafranchian boundary, at ca. 2.0 (Fonelas P-1 and Gerakarou-1). Subcluster B1 is separated into B1.1 of the late Villafranchian to Epivillafranchian faunas of Europe and B1.2 of Dmanisi. The former is further divided into B1.1.1 (~2.0–1.5 Ma) and B1.1.2 (after ~1.6–1.5 Ma). This distinction signifies the small-scale reorganization of the faunas at the genus level during the middle part of the late Villafranchian, leading to the faunas (B1.1.2) of the second part of the late Villafranchian and the Epivillafranchian, which have “lost” middle Villafranchian relics and have mostly “renewed” their faunal elements. B1.1.1 is separated into B1.1.1a, which includes the 2.0–1.8 Ma faunas of Poggio Rosso and Pantalla (both within Olivola/Tasso FUs [1,81]; and B1.1.1b, with the faunas dated between ~1.8 and ~1.5 Ma. B1.1.1b, including TSR, represents the time interval, which documents the co-existence of Leptobos and Bison, the earliest occurrences of Praemegaceros in Europe, and is overall characterized by the mixture of middle (“archaic“) Villafranchian components and “modern” late Villafranchian newcomers.
Considering all biochronological data, TSR is younger than the middle Villafranchian (~2.6–2.0 Ma) faunas and older than the latest Villfranchian and Epivillafranchian ones (~1.5–0.8 Ma, e.g, the Orce sites, Apollonia-1, Pirro Nord, Untermassfeld, and Vallonnet). On the basis of the large mammals, a biochronological age between ~1.8 and 1.6–1.5 Ma is expected for TSR.

5. Discussion

5.1. Age Estimation: Age Constraints

The minimum burial ages (1.88 ± 0.16 Ma, 2.10 ± 0.18 Ma, and 1.98 ± 0.18 Ma) provided by the cosmogenic radionuclide method indicate that the normal magnetic polarity interval encompassing Geo 2b and the bottom part of Geo 2a correlates to the Olduvai subchron (1.95–1.78 Ma, C2n). As a result, the whole magnetic polarity sequence corresponds to the late part of the Olduvai and the early part of the succeeding Matuyama Chron, whose base is currently dated to 1.78 Ma [82]. Therefore, an age younger than 1.78 Ma is indicated for the fossiliferous layer (the bulk of Geo 2a), which was deposited during reverse polarity chron C1r. These results are in very good agreement with the biochronological data, which also indicate an upper age limit at ~1.5 Ma for the fossil accumulation. Therefore, an age between 1.78 and ~1.5 Ma (within the first part of the late Villafranchian) is proposed for the TSR fauna (Figure 9).

5.2. Placement of TSR within the Greek and European Context

The TSR fauna lacks typical middle Villafranchian (~2.6–2.0 Ma) large mammal associations, including the cervids Croizetoceros, Metacervocerus, and Eucladoceros; the bovids Gazella, Gazellospira, and Gallogoral; the canid Nyctereutes; and the hyaenids Chasmaporthetes and Pliocrocuta, present in the Saint-Vallier FU, Varshets (Bulgaria), as well as in the Greek fossil record, for example, in Dafnero, Sesklo, Vatera, and Volax [1,22,56,83]. TSR also appears younger than the “Canis” or “Panthera” sites of the Costa San Giacomo FU, Senéze (France), La Puebla de Valverde (Spain) and Slivnitsa, correlated close to the middle/late Villafranchian boundary, and which preserve most of the genera mentioned above [56,61,81,84]. However, TSR retains some middle Villafranchian taxa, which also survived until the beginning of the late Villafranchian, such as Leptobos and Palaeotragus. On the other hand, the TSR fauna includes several late Villafranchian newcomers, such as Pachycrocuta, Pontoceros, Bison and Praemegaceros. The association Palaeotragus + Praemegaceros (+ Leptobos in [85]) is so far only recorded in Fîntîna lui Mitilan (Romania, [86]), the association Leptobos + Bison + Pontoceros + Praemegaceros + two Equus species is documented to date only in Trlica 11–10 [75], while the association Leptobos + Bison + Pontoceros + two Equus species in Taurida Cave [67]. TSR is also close to Dmanisi, whose faunal record includes Palaeotragus, Bison, Pontoceros, Praemegaceros, and two Equus species [11]. With respect to localities of the central/western sector of the Mediterranean region, TSR is chronologically younger than Fonelas-P1, Poggio Rosso, and Upper Valdarno localities, all correlated to the beginning of the late Villafranchian (ca. 2.0–1.8 Ma), and it is older than Venta Micena, Pirro Nord, Fuente Nueva-3, and Barranco León, collectively dated after ca. 1.6–1.5 Ma, fitting best with the Farneta FU of Italy.
Within the regional faunal succession, TSR is younger than Gerakarou 1 yielding the late Villafranchian “marker” Pachycrocuta brevirostris, but otherwise still retaining a more “middle Villafranchian character” in its fauna including the bovid Gazella (in Italy not later than the Costa San Giacomo FU [79,87]) and the suid Sus strozzi (see [88]). On the other hand, TSR is older than Apollonia 1, where middle Villafranchian relics are absent, and typical late Villafranchian components are present, such as Canis apolloniensis (Canis ex gr. mosbachensis), Bison, Praemegaceros, and Pontoceros, accompanied by various caprines/ovibovines [22], the latter, so far, absent from TSR. The current data indicate that TSR is perhaps chronologically close (although faunal and perhaps chronological differences do exist) with Krimni-1 and Kalamoto-2, where leptobovine bovids also coexist with bisons. Outside the Mygdonia Basin, Libakos (Aliakmon Basin) might be of comparable age to TSR. Libakos has yielded Leptobos, Palaeotragus, Pachycrocuta, Pontoceros, Equus altidens, and two Canis species, in association, however, with the more “modern” Hippopotamus antiquus, and the cervids Praedama savini and Dama vallonnetensis [22,89].

6. Conclusions

Cosmogenic radionuclide dating provides numerical age constraints on the normal polarity detected from the magnetostratigraphic analysis below the fossiliferous layer, correlating it to the Olduvai Subchron. The reverse polarity interval identified in the fossiliferous layer, combined with the biochronological data, indicate that the TSR sediments were deposited during the early part of the Matuyama Chron, between 1.78 and ~1.5 Ma. Overall, TSR presents faunal similarities with Dmanisi, the Farneta FU, Trlica 11–10, Taurida, and the Greek faunas of Libakos, Krimni-1, and Kalamoto-2. Our results not only provide constraints for the age of the local mammal faunal succession, but also help to improve correlations on a broader scale between northern Mediterranean and peri-Black Sea sites, thus, leading to more accurate biogeographic, palaeoecological, and taphonomic interpretations during a crucial interval of time for mammal/hominin migrations and faunal turnovers.
The dating of TSR at ca. 1.78–1.5 Ma places the site close to (or shortly after) the age of Dmanisi, which documents the earliest presence of hominins in western Eurasia. Additionally, this interval records the last occurrences of a number of Villafranchian taxa, as well as the onset of a gradual displacement of other taxa by large mammals that characterize the following Epivillafranchian and Galerian large mammal ages [1,7,9]. Inter-regional differences in the timing and extent of faunal renewals, dispersals, and changes in the structure of mammal palaeocommunities reflect the variability in biotic responses to ecological stimuli according to the geographic setting and prevailing environmental conditions [9,90]. Chronological frameworks serve as the building blocks that allow us to assess potential synchroneity, asynchrony, or diachrony in faunal dynamics. Such detailed studies have been lacking until recently for the Pleistocene of the Greek peninsula. Together with the current investigations for the middle Villafranchian [31,66], our study contributes to filling this gap by providing, for the first time, palaeomagnetic and radiometric constraints to a late Villafranchian locality in the Balkans, and in the Mygdonia Basin in particular, which is one of the most important and richest fossiliferous basins in southeastern Europe. Further refinement of the chronostratigraphic framework of the Mygdonia Basin would provide a better understanding of faunal changes within this sedimentary basin and help calibrate biochronological correlations with other, supra-regional faunal complexes.

Author Contributions

Conceptualization, G.E.K., D.S.K., V.T., G.D.K., and K.H.; methodology, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., T.A.E., E.A., M.M. (Mattia Marini), and G.M.; software, G.E.K., M.M. (Matteo Maron), and M.M. (Mattia Marini); validation, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., T.A.E., E.A., M.M. (Mattia Marini), V.T., G.M., G.D.K., and K.H.; formal analysis, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., T.A.E., E.A., and M.M. (Mattia Marini); data curation, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., T.A.E., E.A., M.M. (Mattia Marini), G.M., and G.D.K.; writing—original draft preparation, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., E.A., M.M. (Mattia Marini), and V.T.; writing—review and editing, G.E.K., D.S.K., M.M. (Matteo Maron), M.S., T.A.E., E.A., M.M. (Mattia Marini), V.T., G.M., G.D.K., and K.H.; visualization, G.E.K., M.M. (Matteo Maron), and M.M. (Mattia Marini); project administration, G.E.K., D.S.K., V.T., G.D.K., and K.H.; funding acquisition, K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the European Research Council STG no. 283503 (PaGE) and the European Research Council CoG 724703 (CROSSROADS), both awarded to K. Harvati (University of Tübingen).

Data Availability Statement

Data is contained within the article.

Acknowledgments

We are grateful to D. Giusti, N. Thompson (University of Tübingen) and all the field team members for their valuable contribution. We thank the four reviewers for their constructive comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rook, L.; Martínez-Navarro, B. Villafranchian: The long story of a Plio-Pleistocene European large mammal biochronologic unit. Quat. Int. 2010, 219, 134–144. [Google Scholar] [CrossRef]
  2. Martínez-Navarro, B. Early Pleistocene faunas of Eurasia and hominin dispersals. In Out of Africa I. The first hominin colonization of Eurasia; Fleagle, J.G., Shea, J.J., Grine, F.E., Baden, A.L., Leakey, R.E., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 207–224. [Google Scholar]
  3. Napoleone, G.; Albianelli, A.; Azzaroli, A.; Bertini, A.; Magi, M.; Mazzini, I. Calibration of the upper Valdarno basin to the Plio-Pleistocene for correlating the Apennine continental sequences. Il Quaternario 2003, 16, 131–166. [Google Scholar]
  4. Kahlke, R.D. Late Early Pleistocene European large mammals and the concept of an Epivillafranchian biochron. In Late Neogene and Quaternary Biodiversity and Evolution: Regional Developments and Interregional Correlations. Vol. II. Proceedings of the 18th International Senckenberg Conference (VI International Palaeontological Colloquium in Weimar); Kahlke, R.D., Maul, L.C., Mazza, P., Eds.; Courier Forschungsinstitut Senckenberg: Frankfurt a. M., Germany, 2007; Volume 259, pp. 265–278. [Google Scholar]
  5. Palombo, M.R. Large mammals faunal dynamics in Southwestern Europe during the late Early Pleistocene: Implications for the biochronological assessment and correlation of mammalian faunas. Alpine Mediterr. Quat. 2016, 29, 143–168. [Google Scholar]
  6. Bellucci, L.; Sardella, R.; Rook, L. Large mammal biochronology framework in Europe at Jaramillo: The Epivillafranchian as a formal biochron. Quat. Int. 2015, 389, 84–89. [Google Scholar] [CrossRef]
  7. Kahlke, R.D.; García, N.; Kostopoulos, D.S.; Lacombat, F.; Lister, A.M.; Mazza, P.P.A.; Spassov, N.; Titov, V.V. Western Palaearctic palaeoenvironmental conditions during the Early and early Middle Pleistocene inferred from large mammal communities, and implications for hominin dispersal in Europe. Quat. Sci. Rev. 2011, 30, 1368–1395. [Google Scholar] [CrossRef]
  8. Leroy, S.A.G.; Arpe, K.; Mikolajewicz, U. Vegetation context and climatic limits of the Early Pleistocene hominin dispersal in Europe. Quat. Sci. Rev. 2011, 30, 1448–1463. [Google Scholar] [CrossRef] [Green Version]
  9. Palombo, M.R. Deconstructing mammal dispersals and faunal dynamics in SW Europe during the Quaternary. Quat. Sci. Rev. 2014, 96, 50–71. [Google Scholar] [CrossRef]
  10. Ferring, R.; Oms, O.; Agustí, J.; Berna, F.; Nioradze, M.; Shelia, T.; Tappen, M.; Vekua, A.; Zhvania, D.; Lordkipanidze, D. Earliest human occupations at Dmanisi (Georgian Caucasus) dated to 1.85–1.78 Ma. Proc. Natl. Acad. Sci. USA 2011, 108, 10432–10436. [Google Scholar] [CrossRef] [Green Version]
  11. Lordkipanidze, D.; Jashashvili, T.; Vekua, A.; de Leon, M.S.P.; Zollikofer, C.P.E.; Rightmire, G.P.; Pontzer, H.; Ferring, R.; Oms, O.; Tappen, M.; et al. Postcranial evidence from early Homo from Dmanisi, Georgia. Nature 2007, 449, 305–310. [Google Scholar] [CrossRef]
  12. Palombo, M.R. Faunal dynamics in SW Europe during the late Early Pleistocene: Palaeobiogeographical insights and biochronological issues. In European Early Pleistocene Biogeography and Ecology Based on the Mammal Record: Case Studies and Preliminary Syntheses; Kostopoulos, D.S., Konidaris, G., Tesakov, A., van den Hoek Ostende, L.W., Rook, L., Eds.; Comptes Rendus Palevol: Amsterdam, The Netherlands, 2018; Volume 17, pp. 247–261. [Google Scholar]
  13. Muttoni, G.; Scardia, G.; Kent, D.V. Early hominins in Europe: The Galerian migration hypothesis. Quat. Sci. Rev. 2018, 180, 1–29. [Google Scholar] [CrossRef] [Green Version]
  14. Palombo, M.R. Biochronology of terrestrial mammals and Quaternary subdivisions: A case study of large mammals from the Italian peninsula. Il Quaternario 2009, 22, 291–306. [Google Scholar]
  15. Palombo, M.R. Discrete dispersal bioevents of large mammals in Southern Europe in the post-Olduvai Early Pleistocene: A critical overview. Quat. Int. 2017, 431, 3–19. [Google Scholar] [CrossRef]
  16. Sakellariou-Mane, E.; Psilovikos, A.; Koufos, G. Contribution to the study of Villafranchian in Northern Chalkidiki (Macedonia, Greece). Sci. Annals, Fac. Phys. Mathem., Univ. Thessaloniki 1979, 19, 279–293. [Google Scholar]
  17. Zamanis, A.; Faugères, L.; Bonis, L.d.; Fountoulis, D.; Simeakis, C.; Panayotis, I.; Dimitrakopoulos, R.; Mercier, J.L. Découverte d’une faune de Mammifères du Quaternaire ancien dans les formations du lac Langhada (Macédoine centrale, Grèce). Implications néotectoniques. C. R. Acad. Sci. Paris 1980, 291, 813–816. [Google Scholar]
  18. Koufos, G.D.; Melentis, J.K. New data from the Villafranchian mammal locality of Gerakarou (Macedonia-Greece). Prak. Akad. Ath. 1983, 58, 185–191. [Google Scholar]
  19. Koufos, G.D.; Syrides, G.E.; Kostopoulos, D.S.; Koliadimou, K.K. Preliminary results about the stratigraphy and the palaeoenvironment of Mygdonia Basin, Macedonia, Greece. Geobios 1995, M.S. 18, 243–249. [Google Scholar] [CrossRef]
  20. Tsoukala, E.; Chatzopoulou, K. A new Early Pleistocene (latest Villafranchian) site with mammals in Kalamotó (Mygdonia Basin, Macedonia, Greece) – Preliminary report. Mitt. Komm. Quartärforsch. Österr. Akad. Wiss. 2005, 14, 213–233. [Google Scholar]
  21. Konidaris, G.E.; Tourloukis, V.; Kostopoulos, D.S.; Thompson, N.; Giusti, D.; Michailidis, D.; Koufos, G.D.; Harvati, K. Two new vertebrate localities from the Early Pleistocene of Mygdonia Basin (Macedonia, Greece): Preliminary results. C. R. Palevol 2015, 14, 353–362. [Google Scholar] [CrossRef]
  22. Koufos, G.D.; Kostopoulos, D.S. The Plio-Pleistocene large mammal record of Greece: Implications for early human dispersals into Europe. In Paleoanthropology of the Balkans and Anatolia; Harvati, K., Roksandic, M., Eds.; Springer: Dordrecht, The Netherlands, 2016; pp. 269–280. [Google Scholar]
  23. Kostopoulos, D.S.; Vassiliadou, K.; Koufos, G.D. The beginning of Pleistocene in the Balkan area according to the mammal record; palaeozoogeographical approach. Ann. Géol. Pays Hellén. 2002, 39, 253–278. [Google Scholar]
  24. Spassov, N. Southeastern Europe as a route for the earliest dispersal of Homo toward Europe: Ecological conditions and the timing of the first human occupation of Europe. In Paleoanthropology of the Balkans and Anatolia; Harvati, K., Roksandic, M., Eds.; Springer: Dordrecht, The Netherlands, 2016; pp. 281–290. [Google Scholar]
  25. Tourloukis, V.; Harvati, K. The Palaeolithic record of Greece: A synthesis of the evidence and a research agenda for the future. Quat. Int. 2018, 66, 48–65. [Google Scholar] [CrossRef]
  26. Harvati, K.; Panagopoulou, E.; Runnels, C. The paleoanthropology of Greece. Evol. Anthropol. 2009, 18, 131–143. [Google Scholar] [CrossRef]
  27. Harvati, K.; Röding, C.; Bosman, A.M.; Karakostis, F.A.; Grün, R.; Stringer, C.; Karkanas, P.; Thompson, N.C.; Koutoulidis, V.; Moulopoulos, L.A.; et al. Apidima Cave fossils provide earliest evidence of Homo sapiens in Eurasia. Nature 2019, 571, 500–504. [Google Scholar] [CrossRef]
  28. Athanassiou, A. Pleistocene vertebrates from the Kyparíssia lignite mine, Megalopolis Basin, S. Greece: Rodentia, Carnivora, Proboscidea, Perissodactyla, Ruminantia. Quat. Int. 2018, 497, 198–221. [Google Scholar] [CrossRef]
  29. Athanassiou, A.; Michailidis, D.; Vlachos, E.; Tourloukis, V.; Thompson, N.; Harvati, K. Pleistocene vertebrates from the Kyparíssia lignite mine, Megalopolis Basin, S. Greece: Testudines, Aves, Suiformes. Quat. Int. 2018, 497, 178–197. [Google Scholar] [CrossRef]
  30. Konidaris, G.E.; Athanassiou, A.; Tourloukis, V.; Thompson, N.; Giusti, D.; Panagopoulou, E.; Harvati, K. The skeleton of a straight-tusked elephant (Palaeoloxodon antiquus) and other large mammals from the Middle Pleistocene butchering locality Marathousa 1 (Megalopolis Basin, Greece): Preliminary results. Quat. Int. 2018, 497, 65–84. [Google Scholar] [CrossRef]
  31. Benammi, M.; Aidona, E.; Merceron, G.; Koufos, G.D.; Kostopoulos, D.S. Magnetostratigraphy and chronology of the Lower Pleistocene primate bearing Dafnero fossil site, N. Greece. Quaternary 2020, 3, 22. [Google Scholar] [CrossRef]
  32. Psilovikos, A.A. Palaeogeographic development of the basin and lake of Mygdonia (Langada-Volvi area), Greece. Ph.D. Thesis, Aristotle University of Thessaloniki, Thessaloniki, Greece, 1977. [Google Scholar]
  33. Giusti, D.; Konidaris, G.E.; Tourloukis, V.; Marini, M.; Maron, M.; Zerboni, A.; Thompson, N.; Koufos, G.D.; Kostopoulos, D.S.; Harvati, K. Recursive anisotropy: A spatial taphonomic study of the Early Pleistocene vertebrate assemblage of Tsiotra Vryssi, Mygdonia Basin, Greece. Boreas 2019, 48, 713–730. [Google Scholar]
  34. Konidaris, G.E.; Kostopoulos, D.S.; Koufos, G.D.; Tourloukis, V.; Harvati, K. Tsiotra Vryssi: A new vertebrate locality from the Early Pleistocene of Mygdonia Basin (Macedonia, Greece). In Proceedings of the XIV Annual Meeting of the European Association of Vertebrate Palaeontologists, Haarlem, The Netherlands, 6–10 July 2016; p. 37. [Google Scholar]
  35. Kostopoulos, D.S.; Maniakas, I.; Tsoukala, E. Early bison remains from Mygdonia Basin (Northern Greece). Geodiversitas 2018, 40, 283–319. [Google Scholar] [CrossRef] [Green Version]
  36. Koufos, G.D.; Konidaris, G.E.; Harvati, K. Revisiting Ursus etruscus (Carnivora, Mammalia) from the Early Pleistocene of Greece with description of new material. Quat. Int. 2018, 497, 222–239. [Google Scholar]
  37. Granger, D.E.; Kirchner, J.W.; Finkel, R.C. Quaternary downcutting rate of the New River, Virginia, measured from differential decay of cosmogenic 26Al and 10Be in cave-deposited alluvium. Geology 1997, 25, 107–110. [Google Scholar] [CrossRef]
  38. Granger, D.E.; Muzikar, P.F. Dating sediment burial with in situ-produced cosmogenic nuclides: Theory, techniques, and limitations. Earth Planet. Sci. Lett. 2001, 188, 269–281. [Google Scholar] [CrossRef] [Green Version]
  39. Granger, D.E. Cosmogenic nuclide burial dating in archaeology and paleoanthropology. In Treatise on Geochemistry, 2nd ed.; Holland, H.D., Turekian, K.K., Eds.; Elsevier: Amsterdam, The Netherlands, 2014; pp. 81–97. [Google Scholar]
  40. Schaller, M.; Ehlers, T.A.; Stor, T.; Torrent, J.; Lobato, L.; Christl, M.; Vockenhuber, C. Timing of European fluvial terrace formation and incision rates constrained by cosmogenic nuclide dating. Earth Planet. Sci. Lett. 2016, 451, 221–231. [Google Scholar] [CrossRef]
  41. Nishiizumi, K. Preparation of 26Al AMS standards. Nucl. Instrum. Methods Phys. Res. B 2004, 223–224, 388–392. [Google Scholar] [CrossRef]
  42. Nishiizumi, K.; Imamura, M.; Caffee, M.W.; Southon, J.R.; Finkel, R.C.; McAninch, J. Absolute calibration of 10Be AMS standards. Nucl. Instrum. Methods Phys. Res. B 2007, 258, 403–413. [Google Scholar] [CrossRef]
  43. Norris, T.L.; Gancarz, A.J.; Rokop, D.J.; Thomas, K.W. Half-life of 26Al. J. Geophys. Res. Solid Earth 1983, 88, B331–B333. [Google Scholar] [CrossRef]
  44. Chmeleff, J.; von Blanckenburg, F.; Kossert, K.; Jakob, D. Determination of the 10Be half-life by multicollector ICP-MS and liquid scintillation counting. Nucl. Instrum. Methods Phys. Res. B 2010, 268, 192–199. [Google Scholar] [CrossRef] [Green Version]
  45. Korschinek, G.; Bergmaier, A.; Faestermann, T.; Gerstmann, U.C.; Knie, K.; Rugel, G.; Wallner, A.; Dillmann, I.; Dollinger, G.; von Gostomski, C.L.; et al. A new value for the half-life of 10Be by Heavy-Ion Elastic Recoil Detection and liquid scintillation counting. Nucl. Instrum. Methods Phys. Res. B 2010, 268, 187–191. [Google Scholar] [CrossRef]
  46. Braucher, R.; Merchel, S.; Borgomano, J.; Bourlès, D.L. Production of cosmogenic radionuclides at great depth: A multi element approach. Earth Planet. Sci. Lett. 2011, 309, 1–9. [Google Scholar] [CrossRef]
  47. Borchers, B.; Marrero, S.; Balco, G.; Caffee, M.; Goehring, B.; Lifton, N.; Nishiizumi, K.; Phillips, F.; Schaefer, J.; Stone, J. Geological calibration of spallation production rates in the CRONUS-Earth project. Quat. Geochronol. 2016, 31, 188–198. [Google Scholar] [CrossRef] [Green Version]
  48. Marrero, S.M.; Phillips, F.M.; Borchers, B.; Lifton, N.; Aumer, R.; Balco, G. Cosmogenic nuclide systematics and the CRONUScalc program. Quat. Geochronol. 2016, 31, 160–187. [Google Scholar] [CrossRef] [Green Version]
  49. Lifton, N.; Sato, T.; Dunai, T.J. Scaling in situ cosmogenic nuclide production rates using analytical approximations to atmospheric cosmic-ray fluxes. Earth Planet. Sci. Lett. 2014, 386, 149–160. [Google Scholar] [CrossRef]
  50. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 2001, 4, 1–9. [Google Scholar]
  51. Hammer, Ø.; Harper, D.A.T. Paleontological data analysis; Blackwell Publishing: Oxford, UK, 2006; p. 351. [Google Scholar]
  52. Arribas, A.; Garrido, G.; Viseras, C.; Soria, J.M.; Pla, S.; Solano, J.G.; Garcés, M.; Beamud, E.; Carrión, J.S. A mammalian lost world in Southwest Europe during the Late Pliocene. PLoS ONE 2009, 4, e7127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kahlke, R.D. Die Unterpleistozäne Komplexfundstelle Untermaßfeld-Zusammenfassung des Kenntnisstandes sowie synthetische Betrachtungen zu Genesemodell, Paläoökologie und Stratigraphie. In Das Pleistozän von Untermaßfeld bei Meiningen (Thüringen), Teil 3; Kahlke, R.D., Ed.; Monographien des Römisch-Germanischen Zentralmuseums Mainz, Dr. Rudolf Habelt: Bonn, Germany, 2001; pp. 931–1030. [Google Scholar]
  54. Kahlke, R.D. (Ed.) The Pleistocene of Untermassfeld near Meiningen (Thüringen, Germany). Part 4; Monographien des Römisch-Germanischen Zentralmuseums Mainz: Mainz, Germany, 2020; Volume 40. [Google Scholar]
  55. Koufos, G.D. The Villafranchian mammalian faunas and biochronology of Greece. Boll. Soc. Paleontol. Ital. 2001, 40, 217–223. [Google Scholar]
  56. Spassov, N. The Plio-Pleistocene vertebrate fauna in South-Eastern Europe and the megafaunal migratory waves from the east to Europe. Rev. de Paléobiologie 2003, 22, 197–229. [Google Scholar]
  57. Mazza, P.P.A.; Bertini, A.; Magi, M. The Late Pliocene site of Poggio Rosso (Central Italy): Taphonomy and paleoenvironment. Palaios 2004, 19, 227–248. [Google Scholar] [CrossRef]
  58. Pavia, M.; Zunino, M.; Coltorti, M.; Angelone, C.; Arzarello, M.; Bagnus, C.; Bellucci, L.; Colombero, S.; Marcolini, F.; Peretto, C.; et al. Stratigraphical and palaeontological data from the Early Pleistocene Pirro 10 site of Pirro Nord (Puglia, south eastern Italy). Quat. Int. 2012, 267, 40–55. [Google Scholar] [CrossRef]
  59. Huguet, R.; Saladié, P.; Cáceres, I.; Díez, C.; Rosell, J.; Bennàsar, M.; Blasco, R.; Esteban-Nadal, M.; Gabucio, M.J.; Rodríguez-Hidalgo, A.; et al. Successful subsistence strategies of the first humans in south-western Europe. Quat. Int. 2013, 295, 168–182. [Google Scholar] [CrossRef]
  60. Madurell-Malapeira, J.; Minwer-Barakat, R.; Alba, D.M.; Garcés, M.; Gómez, M.; Aurell-Garrido, J.; Ros-Montoya, S.; Moyà-Solà, S.; Berástegui, X. The Vallparadís section (Terrassa, Iberian Peninsula) and the latest Villafranchian faunas of Europe. Quat. Sci. Rev. 2010, 29, 3972–3982. [Google Scholar] [CrossRef]
  61. Madurell-Malapeira, J.; Ros-Montoya, S.; Espigares, M.P.; Alba, D.M.; Aurell-Garrido, J. Villafranchian large mammals from the Iberian Peninsula: Paleobiogeography, paleoecology and dispersal events. J. Iber. Geol. 2014, 167–178. [Google Scholar] [CrossRef] [Green Version]
  62. Medin, T.; Martínez-Navarro, B.; Rivals, F.; Madurell-Malapeira, J.; Ros-Montoya, S.; Espigares, M.P.; Figueirido, B.; Rook, L.; Palmqvist, P. Late Villafranchian Ursus etruscus and other large carnivorans from the Orce sites (Guadix-Baza basin, Andalusia, southern Spain): Taxonomy, biochronology, paleobiology, and ecogeographical context. Quat. Int. 2017, 431, 20–41. [Google Scholar] [CrossRef]
  63. Vislobokova, I.A.; Agadjanyan, A.K. New data on age of the Pleistocene fauna from the Trlica locality (Montenegro, Central Balkans) and its correlation with other faunas of Europe. Stratigr. Geol. Correl. 2016, 24, 188–202. [Google Scholar] [CrossRef]
  64. Michel, V.; Shen, C.C.; Woodhead, J.; Hu, H.M.; Wu, C.C.; Moullé, P.É.; Khatib, S.; Cauche, D.; Moncel, M.H.; Valensi, P.; et al. New dating evidence of the early presence of hominins in Southern Europe. Sci. Rep. 2017, 7, 10074. [Google Scholar] [CrossRef]
  65. Cherin, M.; D’Allestro, V.; Masini, F. New bovid remains from the Early Pleistocene of Umbria (Italy) and a reappraisal of Leptobos merlai. J. Mamm. Evol. 2019, 26, 201–224. [Google Scholar] [CrossRef]
  66. Kostopoulos, D.S.; Aidona, E.; Benammi, M.; Gkeme, A.; Grasset, L.; Guy, F.; Koufos, G.D.; Kynigopoulou, Z.; Le Maitre, A.; Novello, A.; et al. The Lower Pleistocene primate-bearing fossil site of Dafnero (W. Macedonia, Greece): New data from classic and innovative approaches. In Proceedings of the 15th International Congress of the Geological Society of Greece, Athens, Greece, 22–24 May 2019. [Google Scholar]
  67. Lopatin, A.V.; Vislobokova, I.A.; Lavrov, A.V.; Startsev, D.B.; Gimranov, D.O.; Zelenkov, N.V.; Maschenko, E.N.; Sotnikova, M.V.; Tarasenko, K.K.; Titov, V.V. The Taurida cave, a new locality of Early Pleistocene vertebrates in Crimea. Dokl. Biol. Sci. 2019, 485, 40–43. [Google Scholar] [CrossRef] [PubMed]
  68. Lowrie, W. Identification of ferromagnetic minerals in a rock by coercivity and unblocking temperature properties. Geophys. Res. Lett. 1990, 17, 159–162. [Google Scholar] [CrossRef]
  69. McFadden, P.L.; McElhinny, M.W. Classification of the reversal test in palaeomagnetism. Geophys. J. Int. 1990, 103, 725–729. [Google Scholar] [CrossRef] [Green Version]
  70. Sotnikova, M.; Rook, L. Dispersal of the Canini (Mammalia, Canidae: Caninae) across Eurasia during the Late Miocene to Early Pleistocene. Quat. Int. 2010, 212, 86–97. [Google Scholar] [CrossRef]
  71. Cherin, M.; Bertè, D.F.; Rook, L.; Sardella, R. Re-defining Canis etruscus (Canidae, Mammalia): A new look into the evolutionary history of early Pleistocene dogs resulting from the outstanding fossil record from Pantalla (Italy). J. Mamm. Evol. 2014, 21, 95–110. [Google Scholar] [CrossRef]
  72. Koufos, G.D. The Pleistocene carnivores of the Mygdonia basin (Macedonia, Greece). Ann. de Paléontol. 1992, 78, 205–257. [Google Scholar]
  73. Spassov, N. Biochronology and zoogeographic affinities of the Villafranchian faunas of Bulgaria and South Europe. Hist. Nat. Bulgarica 2000, 12, 89–128. [Google Scholar]
  74. Masini, F.; Sala, B. Large- and small-mammal distribution patterns and chronostratigraphic boundaries from the Late Pliocene to the Middle Pleistocene of the Italian peninsula. Quat. Int. 2007, 160, 43–56. [Google Scholar] [CrossRef]
  75. Vislobokova, I.A.; Agadjanian, A.K. New data on large mammals of the Pleistocene Trlica fauna, Montenegro, the Central Balkans. Paleontol. J. 2015, 49, 651–667. [Google Scholar] [CrossRef]
  76. Vislobokova, I.A.; Titov, V.V.; Lavrov, A.V.; Gimranov, D.O.; Startsev, D.B.; Tasasenko, K.K. Early Pleistocene spiral-horned antelopes (Artiodactyla, Bovidae) from the Taurida Cave (Crimea, Russia). Paleontol. J. 2020, 54, 81–90. [Google Scholar] [CrossRef]
  77. Abbazzi, L. La fauna de cérvidos de Barranco León y Fuente Nueva-3. In Ocupaciones Humanas en el Pleistoceno Inferior y Medio de la Cuenca de Guadix-Baza; Toro, I., Martínez-Navarro, B., Agustí, J., Eds.; Consejería de Cultura, Junta de Andalucía, Arqueología Monografías: Sevilla, Spain, 2010; pp. 273–290. [Google Scholar]
  78. Athanassiou, A. New giraffid (Artiodactyla) material from the Lower Pleistocene locality of Sésklo (SE Thessaly, Greece): Evidence for an extension of the genus Palaeotragus into the Pleistocene. Zitteliana B 2014, 32, 71–89. [Google Scholar]
  79. Masini, F.; Palombo, M.R.; Rozzi, R. A reappraisal of the Early to Middle Pleistocene Italian Bovidae. Quat. Int. 2013, 288, 45–62. [Google Scholar] [CrossRef]
  80. Sorbelli, L.; Madurell-Malapeira, J.; Prat-Vericat, M.; Azzarà, B.; Romi, P.; Gentili, S.; Cherin, M. Primitive Bison from the late Villafranchian site of Pietrafitta (Italy). In Proceedings of the Invertebrate and Vertebrate Paleontology day (IVPday), Pisa, Italy, 17–20 November 2020; pp. 112–116. [Google Scholar]
  81. Sardella, R.; Bellucci, L.; Bona, F.; Cherin, M.; Iurino, D.A.; Rook, L. Before and after the earliest Homo dispersal in Europe: Evidence from the early Pleistocene sites of the Italian Peninsula. In European Early Pleistocene Biogeography and Ecology Based on the Mammal Record: Case Studies and Preliminary Syntheses; Kostopoulos, D.S., Konidaris, G., Tesakov, A., van den Hoek Ostende, L.W., Rook, L., Eds.; Comptes Rendus Palevol: Amsterdam, The Netherlands, 2018; Volume 17, pp. 287–295. [Google Scholar]
  82. Cohen, K.M.; Gibbard, P.L. Global chronostratigraphical correlation table for the last 2.7 million years, version 2019 QI-500. Quat. Int. 2019, 500, 20–31. [Google Scholar] [CrossRef]
  83. Guérin, C.; Faure, M.; Argant, A.; Argant, J.; Crégut-Bonnoure, É.; Debard, É.; Delson, E.; Eisenmann, V.; Hugueney, M.; Limondin-Lozouet, N.; et al. Le gisement pliocène supérieur de Saint-Vallier (Drôme, France): Synthèse biostratigraphique et paléoécologique. Geobios 2004, 37, S349–S360. [Google Scholar]
  84. Delson, E.; Faure, M.; Guerin, C.; Aprile, L.; Argant, J.; Blackwell, B.A.B.; Debard, E.; Harcourt-Smith, W.; Martin-Suarez, E.; Monguillon, A.; et al. Franco-American renewed research at the Late Villafranchian locality of Senèze (Haute-Loire, France). Cour. Forsch.-Inst. Senckenberg 2006, 256, 275–290. [Google Scholar]
  85. Radulescu, C.; Samson, P.M. Biochronology and evolution of the Early Pliocene to the Early Pleistocene mammalian faunas of Romania. Boll. Soc. Paleontol. Ital. 2001, 40, 285–291. [Google Scholar]
  86. Terhune, C.E.; Curran, S.; Croitor, R.; Drăgușin, V.; Gaudin, T.; Petculescu, A.; Robinson, C.; Robu, M.; Werdelin, L. Early Pleistocene fauna of the Olteţ River Valley of Romania: Biochronological and biogeographic implications. Quat. Int. 2020, 553, 14–33. [Google Scholar] [CrossRef]
  87. Bellucci, L.; Sardella, R. The last Antilopini bovids from the Early Pleistocene of Italy. Quat. Int. 2015, 357, 245–252. [Google Scholar] [CrossRef]
  88. Martínez-Navarro, B.; Madurell-Malapeira, J.; Ros-Montoya, S.; Espigares, M.P.; Medin, T.; Hortolà, P.; Palmqvist, P. The Epivillafranchian and the arrival of pigs into Europe. Quat. Int. 2015, 389, 131–138. [Google Scholar] [CrossRef]
  89. Steensma, K.J. Plio-/Pleistozäne Großsäugetiere (Mammalia) aus dem Beckem von Kastoria/Grevena, südlich von Neapolis - NW-Griechenland. Ph.D. Thesis, Technische Universität Clausthal, Clausthal, Germany, 1988. [Google Scholar]
  90. Palombo, M.R. To what extent could functional diversity be a useful tool in inferring ecosystem responses to past climate changes? Quat. Int. 2016, 413, 15–31. [Google Scholar] [CrossRef]
Figure 1. Selected sites of the Caucasus, the northern peri-Black Sea territory, and the northern peri-Mediterranean region, attributed to the late Villafranchian and Epivillafranchian, ~2.0–0.8 Ma (map from Natural Earth, naturalearthdata.com).
Figure 1. Selected sites of the Caucasus, the northern peri-Black Sea territory, and the northern peri-Mediterranean region, attributed to the late Villafranchian and Epivillafranchian, ~2.0–0.8 Ma (map from Natural Earth, naturalearthdata.com).
Quaternary 04 00001 g001
Figure 2. Geological map of the Neogene and Quaternary lithostratigraphic units of the Mygdonia Basin (left) and simplified composite stratigraphic column of the basin indicating the position of the fossiliferous localities (right). Modified from [21] and data from [19]. Map of Greece from Copernicus Land Monitoring Service (https://land.copernicus.eu/).
Figure 2. Geological map of the Neogene and Quaternary lithostratigraphic units of the Mygdonia Basin (left) and simplified composite stratigraphic column of the basin indicating the position of the fossiliferous localities (right). Modified from [21] and data from [19]. Map of Greece from Copernicus Land Monitoring Service (https://land.copernicus.eu/).
Quaternary 04 00001 g002
Figure 3. Tsiotra Vryssi site and excavation photos. (a) General view of the site; (b) Closer view of the site; (c) Cranium of a bear (Ursus etruscus); (d) Complete front limb of a horse (Equus sp.); (e) Accumulation of isolated and articulated bones of a large bovid, a horse, and a rhino (Stephanorhinus sp.); (f) Skeletal remains of a mammoth (Mammuthus meridionalis).
Figure 3. Tsiotra Vryssi site and excavation photos. (a) General view of the site; (b) Closer view of the site; (c) Cranium of a bear (Ursus etruscus); (d) Complete front limb of a horse (Equus sp.); (e) Accumulation of isolated and articulated bones of a large bovid, a horse, and a rhino (Stephanorhinus sp.); (f) Skeletal remains of a mammoth (Mammuthus meridionalis).
Quaternary 04 00001 g003
Figure 4. (a) Stratigraphic log (modified from [33]) of the Tsiotra Vryssi section with location of the samples for magnetostratigraphy and experiments for rock magnetism, and for cosmogenic radionuclides dating; (b) Latitude of the virtual geomagnetic poles (VGP) with interpreted magnetostratigraphy (black is normal polarity, white reverse polarity); (c) Magnetic susceptibility κ.
Figure 4. (a) Stratigraphic log (modified from [33]) of the Tsiotra Vryssi section with location of the samples for magnetostratigraphy and experiments for rock magnetism, and for cosmogenic radionuclides dating; (b) Latitude of the virtual geomagnetic poles (VGP) with interpreted magnetostratigraphy (black is normal polarity, white reverse polarity); (c) Magnetic susceptibility κ.
Quaternary 04 00001 g004
Figure 5. (A) Isothermal remanent magnetization (IRM) backfield acquisition curves and (B) thermal decay curves of three-axis IRM of selected samples from Tsiotra Vryssi (see Figure 4 for location).
Figure 5. (A) Isothermal remanent magnetization (IRM) backfield acquisition curves and (B) thermal decay curves of three-axis IRM of selected samples from Tsiotra Vryssi (see Figure 4 for location).
Quaternary 04 00001 g005
Figure 6. Vector end-point demagnetization diagrams of key samples from the Tsiotra Vryssi section in in situ coordinates. Closed squares are projections onto the horizontal plane, and open squares are projections onto the vertical plane. Applied AF covers a range 0–100 mT with steps of 5 mT up to 50 mT and of 10 mT up to 100 mT.
Figure 6. Vector end-point demagnetization diagrams of key samples from the Tsiotra Vryssi section in in situ coordinates. Closed squares are projections onto the horizontal plane, and open squares are projections onto the vertical plane. Applied AF covers a range 0–100 mT with steps of 5 mT up to 50 mT and of 10 mT up to 100 mT.
Quaternary 04 00001 g006
Figure 7. Equal-area projections of the ChRM component directions from the Tsiotra Vryssi section (see Table 3 for mean direction).
Figure 7. Equal-area projections of the ChRM component directions from the Tsiotra Vryssi section (see Table 3 for mean direction).
Quaternary 04 00001 g007
Figure 8. Hierarchical cluster analysis at the genus level showing the faunal similarities among selected middle Villafranchian to Epivillafranchian mammal assemblages of western Eurasia.
Figure 8. Hierarchical cluster analysis at the genus level showing the faunal similarities among selected middle Villafranchian to Epivillafranchian mammal assemblages of western Eurasia.
Quaternary 04 00001 g008
Figure 9. Correlation with the geomagnetic polarity time scale (GPTS) [82], simple burial ages provided by cosmogenic radionuclides, and biochronological range of selected large mammals from Tsiotra Vryssi.
Figure 9. Correlation with the geomagnetic polarity time scale (GPTS) [82], simple burial ages provided by cosmogenic radionuclides, and biochronological range of selected large mammals from Tsiotra Vryssi.
Quaternary 04 00001 g009
Table 1. Preliminary faunal list of the large mammals from Tsiotra Vryssi (excavation seasons 2014–2019). Faunal list updated from [21], based on [34,35,36] and unpublished data.
Table 1. Preliminary faunal list of the large mammals from Tsiotra Vryssi (excavation seasons 2014–2019). Faunal list updated from [21], based on [34,35,36] and unpublished data.
OrderFamilyGenusSpecies
ProboscideaElephantidaeMammuthusmeridionalis
CarnivoraCanidaeCanissp.
UrsidaeUrsusetruscus
HyaenidaePachycrocutabrevirostris
FelidaeMegantereonsp.
PerissodactylaEquidaeEquussp. (medium-sized)
Equussp. (large-sized)
RhinocerotidaeStephanorhinussp.
ArtiodactylaGiraffidaePalaeotragussp.
CervidaeCervussp.
Praemegacerossp.
BovidaePontocerossp.
Leptobossp.
Bisoncf. degiulli
Table 2. Cosmogenic radionuclide information for samples from the stratigraphic layer Geo 2b of Tsiotra Vryssi.
Table 2. Cosmogenic radionuclide information for samples from the stratigraphic layer Geo 2b of Tsiotra Vryssi.
Sample IDLab ID Grain Size (μm)Qtz Dissolved
(g)
m(27Al)
(ppm)
26Al/27Al 1σ Error
(%)
26Al Conc.
104 atoms/g(qtz)
m(9Be)
(mg)
10Be/9Be 1σ Error
(%)
10Be Conc.
104 atoms/g(qtz)
26Al/10Be Simple Burial Age
(Myr)
COSMO 1 GL11250–50046.76121 1.39 × 10−134.46 37.07 ± 1.910.3397 2.16 × 10−133.6010.51 ± 0.423.53 ± 0.231.54 ± 0.10
COSMO 2 GL7250–50076.5281 1.80 × 10−135.23 32.48 ± 1.700.2868 4.10 × 10−133.3410.23 ± 0.343.18 ± 0.201.75 ± 0.10
COSMO 3 GL8250–50055.9765 2.20 × 10−135.44 31.88 ± 1.740.2874 2.76 × 10−133.519.40 ± 0.333.39 ± 0.221.62 ± 0.10
Table 3. Palaeomagnetic mean direction and pole form the Tsiotra Vryssi characteristic remanent magnetization (ChRM) component directions. N, number of samples; k, Fisher precision parameter; α95, radius of the cone of 95% confidence; Dec., mean declination of the ChRM component; Inc., mean inclination of the ChRM component; Lat., latitude of palaeomagnetic pole; Long., longitude of palaeomagnetic pole; A95, radius of the cone of 95% confidence around the palaeomagnetic pole.
Table 3. Palaeomagnetic mean direction and pole form the Tsiotra Vryssi characteristic remanent magnetization (ChRM) component directions. N, number of samples; k, Fisher precision parameter; α95, radius of the cone of 95% confidence; Dec., mean declination of the ChRM component; Inc., mean inclination of the ChRM component; Lat., latitude of palaeomagnetic pole; Long., longitude of palaeomagnetic pole; A95, radius of the cone of 95% confidence around the palaeomagnetic pole.
Mean Directions
In Situ
Palaeopole
Nkα95Dec.Inc.Lat.Long.A95
239.010.7°176.8 °E−40.5°72.4 °N213.2 °E7.8°
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Konidaris, G.E.; Kostopoulos, D.S.; Maron, M.; Schaller, M.; Ehlers, T.A.; Aidona, E.; Marini, M.; Tourloukis, V.; Muttoni, G.; Koufos, G.D.; et al. Dating of the Lower Pleistocene Vertebrate Site of Tsiotra Vryssi (Mygdonia Basin, Greece): Biochronology, Magnetostratigraphy, and Cosmogenic Radionuclides. Quaternary 2021, 4, 1. https://0-doi-org.brum.beds.ac.uk/10.3390/quat4010001

AMA Style

Konidaris GE, Kostopoulos DS, Maron M, Schaller M, Ehlers TA, Aidona E, Marini M, Tourloukis V, Muttoni G, Koufos GD, et al. Dating of the Lower Pleistocene Vertebrate Site of Tsiotra Vryssi (Mygdonia Basin, Greece): Biochronology, Magnetostratigraphy, and Cosmogenic Radionuclides. Quaternary. 2021; 4(1):1. https://0-doi-org.brum.beds.ac.uk/10.3390/quat4010001

Chicago/Turabian Style

Konidaris, George E., Dimitris S. Kostopoulos, Matteo Maron, Mirjam Schaller, Todd A. Ehlers, Elina Aidona, Mattia Marini, Vangelis Tourloukis, Giovanni Muttoni, George D. Koufos, and et al. 2021. "Dating of the Lower Pleistocene Vertebrate Site of Tsiotra Vryssi (Mygdonia Basin, Greece): Biochronology, Magnetostratigraphy, and Cosmogenic Radionuclides" Quaternary 4, no. 1: 1. https://0-doi-org.brum.beds.ac.uk/10.3390/quat4010001

Article Metrics

Back to TopTop