Hostname: page-component-76fb5796d-45l2p Total loading time: 0 Render date: 2024-04-25T10:11:52.588Z Has data issue: false hasContentIssue false
  • English
  • Français

A Bayesian Approach for Timing the Neolithization in Mediterranean Iberia

Published online by Cambridge University Press:  10 August 2017

Oreto García-Puchol Open the ORCID record for Oreto García-Puchol [Opens in a new window] ,
Joan Bernabeu-Aubán Open the ORCID record for Joan Bernabeu-Aubán [Opens in a new window] ,
C Michael Barton ,
Salvador Pardo-Gordó Open the ORCID record for Salvador Pardo-Gordó [Opens in a new window] ,
Sarah B McClure  and
Agustín Diez-Castillo Open the ORCID record for Agustín Diez-Castillo [Opens in a new window]
Oreto García-Puchol*
Affiliation:
PREMEDOC Research Group. Departament de Prehistòria, Arqueologia i Història Antiga, Universitat de València, Av. Blasco Ibáñez 28, 46010 Valencia, Spain
Joan Bernabeu-Aubán
Affiliation:
PREMEDOC Research Group. Departament de Prehistòria, Arqueologia i Història Antiga, Universitat de València, Av. Blasco Ibáñez 28, 46010 Valencia, Spain
C Michael Barton
Affiliation:
Center for Social Dynamics and Complexity, Arizona State University, Tempe, AZ, USA
Salvador Pardo-Gordó
Affiliation:
PREMEDOC Research Group. Departament de Prehistòria, Arqueologia i Història Antiga, Universitat de València, Av. Blasco Ibáñez 28, 46010 Valencia, Spain
Sarah B McClure
Affiliation:
Department of Anthropology, Pennsylvania State University, 409 Carpenter Building, University Park, PA, 16802, USA
Agustín Diez-Castillo
Affiliation:
Research Group GRAM, Departament de Prehistòria, Arqueologia i Història Antiga, Universitat de València, Spain
*
*Corresponding author. Email: oreto.garcia@uv.es.
Rights & Permissions [Opens in a new window]

Abstract

In this paper, we compile recent 14C dates related to the Neolithic transition in Mediterranean Iberia and present a Bayesian chronological approach for testing the dual model, a mixed model proposed to explain the spread of farming and husbandry processes in eastern Iberia. The dual model postulates the coexistence of agricultural pioneers and indigenous Mesolithic foraging groups in the Middle Holocene. We test this general model with more regional models of four geographical areas (Northeast, Upper, and Middle Ebro Valley, and Eastern and South/Southeastern regions) and present a filtered summed probability of all 14C dates known in the region in order to compare socioecological dynamics over a long period. Finally, we discuss the results and analyze how certain specific characteristics of sites and their chronologies can serve for timing the Neolithic expansion in Mediterranean Iberia.

Keywords

Bayesian chronologyMediterranean IberiaNeolithization processSCDPDsocioecological dynamics
Type
Research Article
Information
Radiocarbon , Volume 60 , Issue 1 , February 2018 , pp. 181 - 205
Copyright
© 2017 by the Arizona Board of Regents on behalf of the University of Arizona 

INTRODUCTION

The study of the Neolithic transition around the Mediterranean has given rise to a long-standing debate about the key drivers of change, involving explanations largely rooted in cultural diffusion theory and differentially emphasizing demic and/or indigenous components. Building on the pioneer work of Anmerman and Cavalli Sforza (Reference Ammerman and Cavalli-Sforza1984), current literature presents research on the Neolithic spread at both pan-Mediterranean and regional scales, assessed against various compilations (filtered or unfiltered) of available radiocarbon datasets (Zilhao Reference Zilhão2001, Reference Zilhão2011; Manen and Sabatier Reference Manen and Sabatier2003; Fort Reference Fort2012; Bernabeu and Martí Reference Bernabeu Aubán and Martí2014; Isern et al. Reference Isern, Fort, Carvalho, Gibaja and Ibañez2014).

The great increase in the number of radiocarbon (14C) dates in recent years, together with the more accurate selection of samples—thanks to advances in accelerator mass spectrometry (AMS) techniques and the introduction of filters to control the quality of the samples—is providing a more robust chronological framework to support this work. This is the case in Iberia, where a notable growth of 14C dates supplies the basis for exploring alternative hypotheses proposed to account for the spread of farming using Bayesian models for understanding the timing of the Neolithic transition.

In this paper, we use Bayesian chronology models to evaluate the dual model suggested to explain the appearance of the Neolithic in the Valencia region—the Eastern Mediterranean area of the Iberian Peninsula (Figure 1; Bernabeu Reference Bernabeu Aubán2002; Juan Cabanilles and Martí Reference Juan Cabanilles and Martí2002; Bernabeu and Martí Reference Bernabeu Aubán and Martí2014). The dual model combines demic spread of Neolithic farmers with the possibility of contact between Neolithic settlers and Mesolithic groups in a way similar to what has been proposed in northern Europe (Zvelevil and Rowley Conwy Reference Zvelebil and Rowley Conwy1986). The scenario postulates the coexistence of some pioneer Neolithic settlers with the Final Mesolithic groups in different territories during an indeterminate temporal interval.

Figure 1 Map of Iberia showing the regions discussed.

To test this model, we built a Bayesian framework with accurate 14C dates taken from a carefully assessed archaeological context. We use the Bayesian chronological models that have been created by OxCal developers (Bronk Ramsey Reference Bronk Ramsey2000, Bronk Ramsey Reference Bronk Ramsey2009). As a complementary approach, we have used summed 14C probability distributions (SCDPD) to observe fluctuations in population dynamics following assumptions well explained in several key contributions (Shennan and Edinborough Reference Shennan and Edinborough2007; Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013). We filtered the 14C dataset in order to produce a more accurate chronology that allows us to distinguish economic patterns (foragers or producers) within a well-defined chrono-cultural context.

As a first step, we compiled all the available 14C data, and then selected only single fragment samples from short-life plant and animal taxa to include in the analysis. Within the general area of Mediterranean Spain, we examine the evidence from several sub-regions in order to compare chronological models at different geographical scales.

Archaeological Background

Mediterranean Iberia is a wide region with a varied terrain of littoral planes and mountainous regions, transected by two main natural corridors for human travel: the Mediterranean coast (N/S) and the Ebro valley (E/W) (Figure 1). With generally temperate weather conditions (though varying by latitude and altitude) charcoal and pollen analyses reveal the progressive increase of forest cover since the beginning of the Holocene period, reflected by the expansion of oak and the upward migration of Pinus sp and Juniperus sp (Carrión et al. Reference Carrión, Fernández, González Sampériz, Gil Romera, Badal, Carrión Marco, López-Merino, López Sáez, Fierro and Burjachs2010). Holocene sea-level rise was a major factor in modifying the littoral landscape. Early Mesolithic settlements are well known, although the evidence is not uniformly distributed, concentrated in some areas and absent from others. The Early Mesolithic lithic industry is based on flakes, with notches and denticulates, and extends over a long period (ca. 8100–6800 cal BC) in the territory between Catalonia, the Ebro valley, and the Valencian region (Alday 2006). Evidence from the Andalusian region is not currently included in this techno-typological complex.

Mesolithic subsistence included a wide spectrum of medium and small mammals, wild plants, and terrestrial and marine mollusks. Although cave and rock-shelter sites are common, an increasing number of open-air sites are being identified, sometimes including necropoli such as El Collado (Aparicio Reference Aparicio2008; Gibaja et al. Reference Gibaja, Subirà, Terradas, Santos, Agulló, Gómez Martínez, Allièse and Fernández-López de Pablo2015). Coincident with the arrival of the Atlantic Period, from the middle of the 7th millennium cal BC, an abrupt change in lithic industries can be observed. They include standardized bladelets employed to make geometric microliths using the microburin technique. Regionally, this is termed the Late/Final Mesolithic or Geometric Mesolithic. Within this tradition, geographic variation is observed, such as the absence of Geometric Mesolithic assemblages in Catalonia (Vaquero and García Argüelles Reference Vaquero Rodríguez and García-Argüelles i Andreu2009) and the scarcity of evidence in Andalusia (Aura et al. Reference Aura Tortosa, Jordá Pardo, García Borja, García Puchol, Badal García, Pérez Ripoll, Pérez Jordá, Pascual Benito, Carrión Marco and Morales Pérez2013). On the other hand, the Ebro valley and the central (Valencian) portion of Mediterranean Spain show an increased number of sites in caves and rock shelters, and open air sites, although not all have been studied nor are all well preserved. Coastal-interior seasonal movement seems to characterize the subsistence strategy in the Valencian region (Martí et al. Reference Martí, JE, Juan-Cabanilles, García Puchol and Fernández2009), while Mediterranean links are evident in the presence of marine shells (Columbella rustica) at several sites located in the upper Ebro valley.

Two main phases can be recognized with some distinctive regional features. In the first phase (second half of the 7th millennium cal BC) trapeze forms predominate among geometrics, whereas triangular forms (Cocina type triangles) are more common in the second phase (first half of the 6th millennium cal BC). Around 5650 cal BC some archaeological deposits first show evidence of domesticates, together with some new cultural material such as pottery, within the Mediterranean Impressed Ware tradition (see Bernabeu and Martí Reference Bernabeu Aubán and Martí2014). Recently, new discoveries indicate the possibility of at least two independent trajectories of Neolithic colonization or influence to the Mediterranean Iberian Peninsula (Bernabeu and Martí Reference Bernabeu Aubán and Martí2014). This new information may also reveal different temporal rhythms in Neolithic expansion, as well as possible different routes (Manen et al. Reference Guilaine and Manen2007, García Borja et al. Reference García Borja, Aura Tortosa, Jordá Pardo and Salazar García2014). As we show below, several chronological patterns can be observed in relation to the transition to farming throughout this large area, and we explore their possibilities by applying SCDPD and Bayesian chronological models.

MATERIAL AND METHODS

Assessment of 14C Dataset

As stated above, we compiled all available 14C data for Iberia between 7650 and 6000 uncalibrated BP (ca. 6600/5000 cal BC). This time span covers the Late/Final Mesolithic and Early Neolithic periods in Mediterranean Iberia. The resulting database was designed to also collect relevant information regarding the spatial and cultural context of each published 14C date. For constructing Bayesian chronological models, we work only with dates taken from short-lived, single-fragment samples, with SD <100, in order to establish a more accurate approximation of the chronology of archaeological deposits (Table 1). Several studies have pointed out the potentially problematic effects of using 14C dates from all contexts for estimating the timing of the Neolithic expansion in the Western Mediterranean (Zilhao Reference Zilhão2001; Bernabeu Reference Bernabeu Aubán, Molina, Díez Castillo and Orozco Köhler2006), particularly the old-wood effect that can produce anomalously older dates. In addition, problems related to the absence of an accurate sample selection protocol have been observed in multi-component sites (Bernabeu et al. Reference Bernabeu Aubán, Barton and Perez Ripoll2001; Bernabeu Reference Bernabeu Aubán, Molina, Díez Castillo and Orozco Köhler2006). We also discarded samples of marine origin because the reservoir effect can be problematic, and fluctuates in relation to space and time (e.g., Ascough et al. Reference Ascough, Cook and Dugmore2005). Given these considerations, we use only short-lived, single-fragment samples (animal bones, human bones, seeds and fruits, bone artifacts, and charcoal identified as coming from short-lived plants), to help control for the possibility of post-depositional disturbance to archaeological deposits. The presence/absence of domestic species at a site constitutes the premise for determining Neolithic or Mesolithic status.

Table 1 Radiocarbon dates (short-life, singular samples) used to build the chronological models. A: area, PH: phase, GM: Geometric Mesolithic; EN: Early Neolithic.

Table 1 shows the detail of the total number of dates (2016) used here. The start of the Final Mesolithic is taken from the most ancient published date clearly associated with human activity (Cueva de la Cocina Beta-267440, 7610±40 uncalibrated BP: wild goat bone with processing marks). Several dates relating to the end of Mesolithic occupations and the first Neolithic evidence merit special consideration:

Mendandia Ua-34366 (7265±60 BP), level III sup., is from a sample identified originally as aurochs, in a level with pottery but without other domestic elements. Although recent mtDNA analysis of this bone suggests the possibility that it is domestic (Alday 2012), we remain cautious due the absence of any other domesticates at the site. We group this with Mesolithic sites.

Forcas II Beta-247404 (6750±40 BP) and Beta-247405 (6740±40 BP) are from two levels with pottery but without domestic elements, and are therefore also considered as Mesolithic.

Carigüela Col-1565 (6749±39 BP), level XVI—a sample identified as domestic cattle—provides the most ancient Neolithic date in Iberia. This date is part of a series of dates from Neolithic levels recovered from older excavations at this site, which is located in the southern region of Iberia (Mevdev 2013). This series also produced an even older radiocarbon date Coll-1563 (6950, 40) identified as a sheep/goat, but located in a typologically Middle Neolithic level. OxA-1131 (7010±40 BP) is another anomalous date from Carigüela taken from a single sample of a wild horse, from level XV of the older excavations (Carrasco et al. Reference Carrasco, Gámiz, Pachón and Martínez-Sevilla2010). These two earlier dates do not agree with the general chronology of Neolithic expansion elsewhere in the Western Mediterranean. On the other hand, no Mesolithic elements have been identified at the site so far. This suggests a need for a thorough reanalysis of the stratigraphic and archaeological context of the early excavations at the site. For these reasons, we only include Col-1565 in our models.

Peña Larga Beta-242783 (6720±40 BP) comes from a bone identified as domestic sheep/goat from a typologically Neolithic level (together with other domestic animals and pottery) of a site located in the Upper Ebro Valley (Fernández Eraso Reference Fernández Eraso2011). The result is surprising if we consider the oldest expected distribution of early radiocarbon dates from the Mediterranean coast.

Botiqueria GrA-13268 (6040±50 BP) level 6 and GrA-13270 (6240±50 BP) level 8, is associated with two stratigraphical units with impressed pottery separated from a Final Mesolithic level GrA-13267 (6830±50 BP) by a sterile layer (Barandiarán and Cava Reference Barandiarán and Cava2002). No domestic taxa have been recognized, but the assemblage is of small size and poor quality. The presence of triangles and segments with bifacial retouches in units 6 and 8 also suggest late age, in the second half of the 6th millennium cal BC.

A similar case is observed for Abrigo de Ángel 2, Beta-254048 (6390±40 BP) where only two fragments of pottery have been recovered (Domingo et al. 2010), in a similar context with segments and triangles with bifacial retouches characteristic of Neolithic sites. No faunal remains have been found except a few unidentifiable fragments of animal bones.

Valmayor II Beta-341168 (6570±30 BP) and III Beta-341167 (6090±30 BP) date two levels with a significant number of Early Neolithic sherds, and lithic remains similar to those in Botiqueria and Angel 2 (Rojo et al. 2015). But in this case, the faunal remains have been identified as wild species. No domestic seeds have been recovered, despite a well-controlled excavation. Consequently, the authors consider this site as evidence for the acculturation of Mesolithic groups in the area (Rojo et al. 2015).

Botiqueria GrA-13268, GrA-13270, Ángel 2 Beta-254048 and Valmayor Beta-341168, Beta-341167 have been tested alternatively as Mesolithic or Neolithic in the models.

In this paper, we consider the Late/Final Mesolithic as a unit, although archaeological research distinguishes between an early phase where trapezes predominate and a final phase characterised by triangular armatures (Cocina type, similar to Muge type). The same consideration is applied to the Early Neolithic.

We present a dual scale of analysis: macro-regional and regional levels. The macro-regional level comprises a wide coastal region between the Pyrenees and Andalusia, encompassing Catalonia, Valencia, Murcia, and eastern Andalusia; along with the upper and middle Ebro Valley, including the southern Basque Country, Navarra, and Aragon (Figure 1). For the regional level analysis, we have used the major natural corridors mentioned above (Ebro River and the border of the coastal littoral) to divide the region into four main sub-areas. We have designated the areas as northeast (north of the Ebro River in Catalonia), upper Ebro Valley (north of the Aragon Ebro valley, Navarra region, and South Basque country), east (Catalonia south of the Ebro, lower Aragon, and the Valencian Autonomous region), and south/southeast (Murcia and Eastern Andalusia).

Calibration and Phase Modeling

In our analysis, we use Bayesian models as they have been implemented in Oxcal 4.2.4 (Bronk Ramsey Reference Bronk Ramsey2009) to identify the beginning and end of two independent but overlapping phases within the neolithization process: Late/Final Mesolithic and Early Neolithic. The Boundary command identifies an event that has not been directly dated and estimates a probability distribution for its occurrence based on the known dates included in the phases preceeding and following the event (Bronk Ramsey Reference Bronk Ramsey2009). The 14C dates are not ordered stratigraphically, because they derive from different sites. The resulting start and end dates for each phase are presented with 1σ and 2σ ranges similar to calibrated 14C dates. An agreement of 60.0 or higher is considered an acceptable result. We tested two models, using alternative interpretations of the recent 14C dates from Botiquería levels 6 (GrA-13268) and 8 (GrA-13270), Abrigo de Ángel 2 (Beta-254048), Valmayor II and III (Beta-341168 and Beta-341167): as Neolithic (model A) or Mesolithic (model B).

We also examined regional variability in the modeled Mesolithic/Neolithic boundary for the upper/middle Ebro Valley and eastern region. Below, we discuss the chronological distinctions between the two phases (Mesolithic and Early Neolithic) and some problems linked with the differences detected in the regional archaeological contexts (i.e., information gaps).

Building Summed Probabilities Distributions of 14C dates and Sites

With the idea of comparing trends in settlement patterns during the neolithization process, we use a method similar to common summed probabilities that includes a correction for the impact of different numbers of dates for different sites and levels. Alternative proposals for mitigate the impact of sites with multiple radiocarbon dating have been used in other works (Edinborough Reference Edinborough2005; Shennan and Edinborough Reference Shennan and Edinborough2007; Edinborough Reference Edinborough2008; see also Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013 for a more refined method). Although an increasing number of archaeological studies use summed radiocarbon probability distributions (SCDPD) as a proxy for studying long-term prehistoric demography (Gamble et al. Reference Gamble, Davies, Pettitt, Hazelwood and Richards2005; Collard et al. Reference Collard, Edinborough, Shennan and Thomas2010; Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013), indiscriminate use can present several problems (Bamforth and Grund Reference Bamforth and Grund2012; Williams Reference Williams2012; Combre and Robinson 2014). Nevertheless, some authors have developed robust protocols that allow supporting the use of SCDPD as an indicator for fluctuations in human population levels (Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013; Timpson et al. Reference Timpson, Colledge, Crema, Edinborough, Kering, Manning, Thomas and Shennan2014). The assumption of SCDPD analysis is that multiple small, unsystematic samples from a large sample of sites constitute a quasi-random sample of regional occupation trends. In this paper, we consider 14C dates from short-lived and long-lived samples with a standard deviation of less than or equal to 100. We do not demand the same level of reliability as for the chronological model, because to observe trends we assume that an increased number of dates provides more robust results over a wider perspective. We include 14C dates from the beginning of the Holocene until 6000 uncalibrated BP. In all we use a total of 402 dates. We calibrated all the dates using Oxcal Program 4.2.4 and IntCal13 curve for the Northern Hemisphere (Reimer et al. Reference Robb2013) with 1σ (68.2). Then we distributed all dates within intervals of 200 yr. Finally, we selected and summed only one date per site per interval. In the discussion that follows, this is a salient point to add to the explanation over the chronological models constructed.

RESULTS

To quote Box (Reference Box1979): “All models are wrong but some are useful.” With that in mind, we try to interpret the results obtained, and to point out several aspects relevant to socioecological dynamics which can be explored through these Bayesian analyses (Whittle et al. Reference Whittle, Healy and Bayliss2011). We test the dual model building different scenarios in several scales of analysis (see Riede and Edinborough Reference Riede and Edinborough2012; Burley et al. Reference Burley, Edinborough, Weisler and Zhao2015). At the macro scale of the analysis, Model A has an acceptable agreement index of 73.7, while Model B has an agreement index of 59 (Table 2, and supplementary material) and consequently it cannot be considered acceptable. Model A considers Botiquería level 6 (GrA-13268) and 8 (GrA-13270), Abrigo de Ángel 2 (Beta-254048) and Valmayor II and III Beta-341168 and Beta-341167 as Neolithic samples. Therefore, in this model we do not consider the absence of domestic elements as a definitive indicator of a Mesolithic occupation.

Table 2 Bayesian sequence modeled dates from Oxcal 4.2.3 (Bronk Ramsey Reference Cámalich Massieu and Martín Socas2009). Calibration conducted using IntCal13 curve: Northern Hemisphere (Reimer et al. Reference Robb2013).

If this seems to contradict our prior assumptions regarding Neolithic and Mesolithic status (mentioned above), in reality it serves to highlight the temporal distance between the latest Mesolithic 14C date in Botiqueria (GrA-13267, 6830±50) and the first Neolithic dates (GrA13270, 6240±50), and reflects a possible gap in human occupation indicated by the excavators (Barandiarán and Cava Reference Barandiarán and Cava2001). In our opinion these arguments reinforce the results of Model A versus model B. On the other hand, both models present problems with the agreement index for the following dates: Col-1565 (6749±39) and Beta-242783 (6720±40) with an index of 6.2 and 27.4 for model A, and 7 and 30.4 for model B. Treating the date Beta-341167 and Botiquería date GrA-13268 (6040±50) as Mesolithic also fails, with an index of 51.6 and 33.4. Both, Col-1565 (6749±39) and Beta-242783 (6720±40) 14C dates, then, are probably best considered the most ancient dates for Neolithic contexts in Iberia. Col-l565 corresponds to a cattle bone located at the bottom of the Neolithic levels in Carigüela, of the southern region (Mevdev 2013). Unfortunately, no detailed publication is yet available for this interesting cave, which presents one of the most complete assemblage series of impressed pottery in southern Iberia (Navarrete Reference Navarrete Enciso1976; Mevdev 2013). Currently, it is not possible to evaluate the evidence from this site. Recent 14C analyses suggest caution in interpreting certain 14C dates, such as Col-1563 (6950±40) and OxA-1131 (7010±40), due to problems linked with the sample identification (sheep/goat) and/or with the cultural associations of the samples (Mevdev 2013). Alternatively, some problems linked with bone dating protocols could be considered (Wood 2015). While this sample could be associated with pre-Neolithic occupations, we have not yet discarded Col-1565 (6749±39), because it is close in age to other early Neolithic contexts in Iberia.

Beta-242783 (6720±40) corresponds to a sheep/goat sample date found at the Peña Larga rock shelter in the upper Ebro valley (Figure 2). This merits a special mention due to its location in an inland area distant from the Mediterranean Sea and along a main corridor for Neolithic expansion. Accepting this date implies a rapid neolithization process via different inland routes and not only along the Mediterranean coast as suggested by Zilhao (Reference Zilhão2001, Reference Zilhão2011). Alternatively, a combination of taphonomic and/or technical problems affecting the 14C sample could be responsible for this result. For the time being, we consider the date valid, while awaiting further clarifying information, as there are no immediately apparent problems with the calibration curve.

Figure 2 Mesolithic and Neolithic sites in Mediterranean Iberia region used for chronological modeling.

We applied the overlapping model in two regional sub-areas: the upper-middle Ebro valley and the eastern coast of Iberia. Regional models have only been constructed for these two sub-areas, due to the absence of late Mesolithic deposits in the northern (Catalonia) sub-region and the Murcia/eastern Andalusia sub-region. The overlapping model in the upper-middle Ebro valley has an overall agreement index of 91.2, which constitutes a good result and indicates the consistency of the model (Figure 3). The date of Peña Larga (Beta-242783, 6720±40) shows a low agreement (28.1), which again highlights its distinct context discussed above. As we did for the entire area, we consider two hypotheses (A and B) for the Eastern region. In this case Model A considers Botiquería levels 6 (GrA-13268) and 8 (GrA-13270), Abrigo de Ángel 2 (Beta-254048) and Valmayor II and III (Beta-341168 and Beta-341167) as Neolithic, while Model B considers the same dates as Mesolithic dates. Model A in the Eastern region generates the larger agreement index value (Amodel 99), although both present acceptable index values (Figure 4). The most ancient Neolithic date in this region is En Pardo (Beta-231880, 6660±40) (Soler et al. Reference Soler Diaz, García Atienzar, Duque Espino, Ferrer García, Gómez Pérez, Guillem Calatayud, Iborra Eres, Martínez Valle, Pérez Jordà, Roca de Togores Muñoz and Ximénez de Embún Sánchez2013). Model B has an agreement index of 90.4 (Figure 5).

Figure 3 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Ebro valley region.

Figure 4 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Eastern area (model A).

Figure 5 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Eastern area (model B).

Overall, the general models that reflect the validity of hypothesis A (Botiquería, Ángel 2 and Valmayor as Neolithic occupations) are more robust than those of hypothesis B (Table 2 and Figure 6). Accepting hypothesis A implies the coexistence of Mesolithic and Neolithic societies. But thereafter (Table 3 and Figure 6) it is not possible to follow their immediate evolution from the current archaeological record. Model B reflects an even longer time span of coexistence, although the resolution of the stratigraphic context remains a challenge. On a regional level, the models offer more robust results. The Ebro model presents a consistent model where the coexistence between Mesolithic and Neolithic groups is more visible from 14C (DistMin/Max=6/144) (Table 3 and Figure 6). In the Eastern area, the possibility of recognizing an acculturation process through 14C dates is not well reflected by Model A (more robust Amodel: 99).

Figure 6 Phase modeling results showing the end of the Mesolithic and the beginning of the Neolithic according with the different models proposed (general A model, and regional models: Eastern A and B models, and Ebro model).

Table 3 Bayesian sequence modeled dates from Oxcal 4.2.3 (Bronk Ramsey Reference Bronk Ramsey2009). Interval span and distance between the end of the Mesolithic and the start of the Neolithic. Calibration conducted using IntCal13 curve: Northern Hemisphere (Reimer et al. Reference Reimer, Bard, Bayliss, Beck, Blackwell, Bronk Ramsey, Buck, Cheng, Edwards, Friedrich, Grootes, Guilderson, Haflidason, Hajdas, Hatté, Heaton, Hoffmann, Hogg, Hughen, Kaiser, Kromer, Manning, Niu, Reimer, Richards, Scott, Southon, Staff, Turney and van der Plicht2013).

Consequently, we feel this suggests that coexistence of Mesolithic hunter/gatherers and Neolithic farmers is possible even though not well reflected in the archaeological record. The regional variability reinforces this conclusion, and highlights regional differences between the Mediterranean corridor (a quick Neolithic expansion) and the Ebro valley (where acculturation processes had a greater opportunity).

As a complementary analysis to better understand trends in socioecological dynamics among the last hunter gatherers and first farmers in Mediterranean Iberia, we combined all filtered radiocarbon dates (single fragments, short-lived taxa, SC<100) for the Mediterranean region into a SCDPD curve (Figure 7). We compare this curve with a curve of the number of distinct occupations dated to each 200-yr interval throughout the same region. Although there are relatively fewer sites than dates for the later Neolithic, both datasets display similar patterns. An interesting consideration is the increase of dates coincident with the beginning of the Late Mesolithic (Figure 7) associated with the appearance of standardized blade technology and trapezes (Martí et al. Reference Martí, JE, Juan-Cabanilles, García Puchol and Fernández2009; Binder et al. Reference Binder, Collina, Guilbert, Perrin and García Puchol2012). The peaks in both curves show when the Neolithic is well in evidence over much of the Iberian Peninsula (ca. 5300 cal BC) (Bernabeu and Martí Reference Bernabeu Aubán and Martí2014; Bernabeu et al. 2014). The 14C curve begins a steep climb around 5700 cal BC, with the arrival of the Neolithic to the region. That is, a substantial change in 14C densities is coincident with the introduction and expansion of agricultural practices associated with the initial appearance of the Neolithic.

Figure 7 Comparative sumprob radiocarbon dates/sites distributed within intervals of 200 years. We calibrated all the dates using Oxcal Program 4.2 and IntCal13 curve for the Northern Hemisphere (Reimer et al. Reference Reimer, Bard, Bayliss, Beck, Blackwell, Bronk Ramsey, Buck, Cheng, Edwards, Friedrich, Grootes, Guilderson, Haflidason, Hajdas, Hatté, Heaton, Hoffmann, Hogg, Hughen, Kaiser, Kromer, Manning, Niu, Reimer, Richards, Scott, Southon, Staff, Turney and van der Plicht2013) with 1σ (68.2).

DISCUSSION

As pointed out by Bocquet-Appel (Reference Bocquet-Appel2008), indicators of demographic transition have been detected around the world at the beginning of the Neolithic, linking the increase of population to the new modes of subsistence and residential patterns. A growing number of researchers have used the density curves of large series of radiocarbon data as a viable proxy for these signals (Gamble et al. Reference Gamble, Davies, Pettitt, Hazelwood and Richards2005; Shennan and Edinborough Reference Shennan and Edinborough2007; Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013).

Recent studies use this information as demographic proxies at different geographical levels in Europe (Shennan and Ediborough Reference Shennan and Edinborough2007) as well as to show trends of settlement densities (Bernabeu et al. 2014; Bernabeu et al. Reference Bernabeu Aubán, García Puchol, Barton, McClure and Pardo Gordó2015b). The results presented here indicate the clarity of this signal for the Mediterranean Iberian region, where the arrival of agriculture and husbandry coincides with a rapid increase in radiocarbon densities—a pattern that is also observed in other European regions (Shennan et al. Reference Shennan, Downey, Timpson, Edinborough, Colledge, Kerig, Manning and Thomas2013; Timpson et al. Reference Timpson, Colledge, Crema, Edinborough, Kering, Manning, Thomas and Shennan2014).

The general debate in Europe about the social processes responsible for the Neolithic transition centers on the importance of demic expansion (see Ammerman and Cavalli Sforza Reference Ammerman and Cavalli-Sforza1984; Perles Reference Perlès2001; Zilhao Reference Zilhão2001; Guilaine Reference Guilaine2001; Ammerman and Biagi Reference Ammerman and Biagi2003; Forenbaher and Miracle Reference Forenbaher and Miracle2006; Guilaine and Manen et al. Reference Guilaine and Manen2007; Robb Reference Robb2007; Bernabeu et al. Reference Bernabeu Aubán, Molina, Esquembe, Ramón and Boronat2009; Colledge et al. Reference Colledge, Connolly, Dobney, Manning and Shennan2013; Ozdogan Reference Özdoğan2011; Rowley-Conwy Reference Rowley-Conwy2011; Juan-Cabanilles and Garcia Puchol Reference Juan-Cabanilles and García Puchol2013; McClure Reference McClure2013; McClure and Podrug Reference McClure and Podrug2016) and the pace of Neolithic expansion (Bocquet-Appel et al. Reference Bocquet-Appel, Naji, Vander Linden and Kozlowski2012, Fort Reference Fort2012). Despite the discussion about the interpretation of particular sites and dates (see Zilhao Reference Zilhão2011; Alday 2012), or the respective impact of demic diffusion versus indigenous components (Bernabeu Reference Bernabeu Aubán2002; Zilhao Reference Zilhão2003; García Puchol et al. Reference García Puchol, Molina, Aura and Bernabeu Aubán2009; Diaz del Río Reference Diaz del Rio2011; Cruz Berrocal Reference Cruz Berrocal2012; Bernabeu and Martí Reference Bernabeu Aubán and Martí2014; Isern et al. Reference Isern, Fort, Carvalho, Gibaja and Ibañez2014; Bernabeu et al. Reference Bernabeu Aubán, Barton, Pardo Gordó and Bergin2015a, Reference Bernabeu Aubán, García Puchol, Barton, McClure and Pardo Gordó2015b), general agreement exists among researchers in the Iberian Peninsula on the existence of an initial impact of foreign components and the hypothesis of dual drivers of demic expansion and acculturation. In accordance with the dual model proposal we assume that an initial foreign expansion constituted the driving force of the subsequent spread in Iberia.

Chronological analysis using a Bayesian approach allows the exploration of alternative hypotheses about dual demic/acculturation models. In this respect, several points are worth noting. One concerns the latest Mesolithic context considered. We tested two possibilities for the cultural affiliations of several ambiguous dates from sites in the eastern region: Botiqueria, Angel 2, Valmayor II and III. The Bayesian models indicate a more probable Neolithic affiliation for these dated contexts. The existence of a chronological gap between Mesolithic and Neolithic levels in Botiqueria reinforce this interpretation. This would make Aizpea GrA-779 (6600, 50) (Cava 1997) the last authenticated Mesolithic date currently known. Other dates position the initial Neolithic expansion in Iberia at around 5650 cal BC (Col-1565, 6749±39: 5725–5575 cal BC; Beta-242783, 6720±40: 5715–5561 cal BC; OxA-26068, 6655±45: 5644–5491 cal BC, 2σ range; Beta-131577, 6590±40: 5616–5481), all based on bones of possibly domestic animals, although some problems can be highlighted about their domestic status only confirmed for Beta-131577 (Cueva de Nerja, in southern Iberia) (Martins et al. Reference Martins, Oms, Pereira, Pike, Rowsell and Zilhão2015). The oldest domestic seed corresponds to Mas d’Is (Beta-162092, 6600±40, 5617–5485 cal BC), an open-air site located in Eastern Iberia region (Bernabeu et al. 2014). Their wide distribution around the region fuels debate about the existence of different entry routes of Neolithic farmers coming from the north and from the south along the Iberian coast (García Borja et al. Reference García Borja, Aura Tortosa, Jordá Pardo and Salazar García2014).

CONCLUSION

We have tested alternative variants of the dual model, a mixed model that proposes the coexistence between migrant farmers and Mesolithic groups during the Neolithic arrival in Mediterranean Iberia. To do that we have used a large and carefully filtered radiocarbon dataset to compare a SCDPD method as a proxy for diachronic fluctuations in population levels of and Bayesian boundary modeling to examing the potential overlap in final Mesolithic and initial Neolithic occupations. We found variable support for Mesolithic acculturation at the scales of sub-regions and over all of Mediterranean Iberia. The most robust models do not support the existence of a lengthy acculturation process across the entire region, nor in the eastern sub-area. In the Ebro area, there is stronger support for a longer period of acculturation between Mesolithic and Neolithic groups. Although we were unable to construct sub-regional models for the northern (Catalonia) sub-region and the Murcia/eastern Andalusia sub-region, current archaeological data supports a model of Neolithic colonization in those areas.

The combined SCDPD curve built for the entire region of Mediterranean Iberia also sheds new light on discussions of neolithization. An initial increase in 14C densities is linked with the appearance of blade technology and trapezes in the Mesolithic sequence in the region; a much more marked peak in the SCDPD curve corresponds with the first appearance of domestic remains. Consequently, we note that agro-pastoral practices significantly increased the densities—and possibly increased populations across the peninsula (see Figure 7). This rapid increase does not seem to reflect the kind of long-term acculturation process in Mediterranean Iberia as a whole observed in the archaeological record of other European regions (Collard et al. Reference Collard, Edinborough, Shennan and Thomas2010; Rowley Conwy Reference Rowley-Conwy2011).

The detection of necropoli from the Early Mesolithic onward in the area, indicates differences in settlement patterns that could reflect more residential stability, an aspect which also coincides with the increase in radiocarbon date densities. At the same time, the rise in the diversity of resources used is noted from the beginning of the Holocene in hunter-gatherer subsistence patterns (Martí et al. Reference Martí, JE, Juan-Cabanilles, García Puchol and Fernández2009), including sites located along coastal and wetland areas. Despite these data, current information does not seem to indicate high population levels. On the other hand, the signal detected in SCDPD at the beginning of the Neolithic could be linked with the hypothesis of a foreign introduction related to farming and herding practices by migrant farmers. The results of Bayesian boundary modeling offer support for some degree of regional variability in farmer-hunter interaction. A rapid assimilation seems to be the case in the east, and a longer coexistence in the Ebro valley. If this were the case, we would expect to find sites and layers with Mesolithic industries and some Neolithic artifacts, such as pottery. However, a more refined chronology is needed if we want to gain further insight into the problems relating to Mesolithic/Neolithic relationships.

Acknowledgments

This research was funded by the Spanish Ministry of Economy, Industry and Competitiveness grants HAR2012-33111 “MesoCocina: Los últimos caza-recolectores y el paradigma de la neolitización en el mediterráneo peninsular” and HAR2015-68962: EVOLPAST: Dinámicas evolutivas y patrones de variabilidad cultural de los últimos caza-recolectores y el primer neolítico en el este peninsular (circa 7000–4500 cal BC). We would like to acknowledge the anonymous reviewers who helped to improve the earlier version of this paper. SPG is supported by University of Valencia (subprograma Atracció de Talent - Contractes Postdoctorals).

SUPPLEMENTARY MATERIAL

To view supplementary material for this article, please visit https://doi.org/10.1017/RDC.2017.61

References

Alday Ruiz, A. 2009. El final del Mesolítico y los inicios del Neolítico en la península Ibérica: Cronología y Fases. Munibe Antropologia – Arkeologia 60:157–173.Google Scholar
Ammerman, AJ, Biagi, P, editors. 2003. The Widening Harvest: The Neolithic Transition in Europe: Looking Back, Looking Forward. Boston: Archaeological Institute of America. 343 p.Google Scholar
Ammerman, AJ, Cavalli-Sforza, LL. 1984. Neolithic Transition and the Genetics of Populations in Europe. Princeton: Princeton University Press.CrossRefGoogle Scholar
Aparicio, Pérez J. 2008. La necrópolis mesolítica de El Collado (Oliva, Valencia). Varia VIII. Serie Arqueológica. Diputación Provincial, Sección de Estudios Arqueológicos Valencianos.Google Scholar
Ascough, P, Cook, G, Dugmore, A. 2005. Methodological approaches to determining the marine radiocarbon reservoir effect. Progress in Physical Geography 29(4):532–547.CrossRefGoogle Scholar
Aura Tortosa, JE, Jordá Pardo, JF, García Borja, P, García Puchol, O, Badal García, E, Pérez Ripoll, M, Pérez Jordá, G, Pascual Benito, JLl, Carrión Marco, Y, Morales Pérez, JV. 2013. Una perspectiva mediterránea sobre el proceso de neolitización. Los datos de la cueva de Nerja en el contexto de Andalucía (España). Revista de Prehistoria de Andalucía 4:53–77.Google Scholar
Baldellou, V. 2012. La Cueva de Chaves (Bastarás-Casbas, Huesca). Saguntum: Papeles del Laboratorio de Arqueología de Valencia, Extra 12:141–144.Google Scholar
Bamforth, DB, Grund, B. 2012. Radiocarbon calibration curves, summed probability distributions, and Early Paleoindian population trends in North America. Journal of Archaeological Science 39(6):1768–1774.CrossRefGoogle Scholar
Barandiarán, I, Cava, A. 2001. Cazadores-recolectores en el Pirineo navarro: el sitio de Aizpea entre 8.000 y 6.000 años antes de ahora. Veleia, Anejo, series maior, 10: Universidad del País Vasco, Vitoria.Google Scholar
Barandiarán, I, Cava, A. 2002. A propósito de unas fechas del Bajo Aragón: reflexiones sobre el Mesolítico y el Neolítico en la cuenca del Ebro. Spal, (2000)9. Sevilla 293–326.Google Scholar
Bernabeu Aubán, J. 2002. The social and the symbolic context of Neolithization. Saguntum: Papeles del Laboratorio de Arqueología de Valencia, Extra 5:209–233.Google Scholar
Bernabeu Auban, J. 2006. Una visión actual sobre el origen y difusión del Neolítico. En La Península Ibérica. In: García O, Aura JE, editors. 8.000 Años de Ocupación Humana en la Cabecera del río de Alcoi . Alcoi: Museu d’Alcoi. p 189–211.Google Scholar
Bernabeu Aubán, J, Martí, B. 2014. The first agricultural groups in the Iberian Peninsula. In: Manen C, Perrin T, Guilaine J, editors. Transitions in the Mediterranean: How Hunters Became Farmers. Errance. p 419–438.Google Scholar
Bernabeu Aubán, J, Molina, L, editors. 2009. La Cova de Les Cendres, Serie Mayor. Alicante: Fundación MARQ.Google Scholar
Bernabeu Aubán, J, Barton, CM, Perez Ripoll, M. 2001. A taphonomic perspective on Neolithic beginnings: theory, interpretation, and empirical data in the Western Mediterranean. Journal of Archaeological Science 28(6):597–612.CrossRefGoogle Scholar
Bernabeu Aubán, J, Molina, LL, Díez Castillo, A, Orozco Köhler, T. 2006. Inequalities and power. Three millennia of Prehistory in Mediterranean Spain (5.600-2.000 cal BC). In: Diaz-del-Río P, García Sanjuán L, editors. Social Inequality in Iberian Late Prehistory. British Archaeological Reports, International Series 1525:97–116.Google Scholar
Bernabeu Aubán, J, Molina, LL, Esquembe, MA, Ramón, J, Boronat, JD. 2009. La cerámica impresa mediterrána en el origen del Neolítico de la península Ibérica. In: De Méditerranée et d’ailleurs...Mélanges offerts à Jean Guilaine. Archives d’Écologie Préhistorique. p 83–95.Google Scholar
Bernabeu Aubán, J, Díez, A, Orozco, T. 2014. Campañas de excavación recientes en el yacimiento del Mas d’Is (Penàguila, Alacant). In: Olcina Doménech MH, Soler Díaz JA, editors. Arqueología en Alicante en la primera década del siglo XXI. II Jornadas de Arqueología y Patriomonio Alicantino, MARQ. Arqueología y Museos, Extra 01: Alicante. p 183–8.Google Scholar
Bernabeu Aubán, J, García Puchol, O, Pardo, S, Barton, M, McClure, SB. 2014. Socioecological dynamics at the time of Neolithic transition in Iberia. Environmental Archaeology. DOI: 10.1179/1749631414Y.0000000032.CrossRefGoogle Scholar
Bernabeu Aubán, J, Barton, CM, Pardo Gordó, S, Bergin, SM. 2015a. Modeling initial Neolithic Dispersal. The first agricultural groups on West Mediterranean. Ecological Modeling 307:22–31.CrossRefGoogle Scholar
Bernabeu Aubán, J, García Puchol, O, Barton, CM, McClure, S, Pardo Gordó, S. 2015b. Radiocarbon dates, climatic events, and social dynamics during the Early Neolithic in Mediterranean Iberia. Quaternary International. http://dx.doi.org/10.1016/j.quaint.2015.09.020.CrossRefGoogle Scholar
Binder, D, Collina, C, Guilbert, R, Perrin, T, García Puchol, O. 2012. Pressure-knapping blade production in the north-western Mediterranean region during the seventh millennium cal B.C. In: Desrosiers PM, editor. The Emergence of Pressure Blade Making: from Origin to Modern Experimentation. Springer. p 199–217.CrossRefGoogle Scholar
Bocquet-Appel, J-P. 2008. Explaining the Neolithic demographic transition. In Bocquet-Appel J-P, Bar-Yosef O, editors. The Neolithic Demographic Transition and its Consequences. New York: Springer. p 35–56.CrossRefGoogle Scholar
Bocquet-Appel, JP, Naji, S, Vander Linden, M, Kozlowski, J. 2012. Understanding the rates of expansion of the farming system in Europe. Journal of Archaeological Science 39:531–546.CrossRefGoogle Scholar
Box, GEP. 1979. Robustness in the strategy of scientific model building. In: Launer RL, Wilkinson GN, editors. Robustness in Statistics. New York: Academic Press. p 201–236.CrossRefGoogle Scholar
Bronk Ramsey, C. 2000. Comment on the use of Bayesian statistics for 14C dates of chronological ordered samples: a critical analysis. Radiocarbon 42(2):199–202.CrossRefGoogle Scholar
Bronk Ramsey, C. 2009. Bayesian analysis of radiocarbon dates. Radiocarbon 51(1):337–360.CrossRefGoogle Scholar
Burley, DV, Edinborough, K, Weisler, M, Zhao, J. 2015. Bayesian modelling and chronological precision for Polinesian settlement of Tonga. PLoS One: 23 March 2015.CrossRefGoogle Scholar
Cámalich Massieu, MD, Martín Socas, D. 2013. Los Inicios Del Neolítico En Andalucía. Entre La Tradición Y La Innovación. Menga 4:103–132.Google Scholar
Carrasco, J, Gámiz, J, Pachón, JA, Martínez-Sevilla, F. 2010. El poblamiento neolítico en los dominios penibéticos del poniente granadino. Cuadernos de Prehistoria de la Universidad de Granada 20:223–294.Google Scholar
Carrión, JS, Fernández, S, González Sampériz, P, Gil Romera, G, Badal, E, Carrión Marco, Y, López-Merino, L, López Sáez, JA, Fierro, E, Burjachs, F. 2010. Expected trends and surprises in the Lateglacial and Holocene vegetation history of the Iberian Peninsula and Balearic Islands. Review of Palaeobotany and Palynology 162:458–475.CrossRefGoogle Scholar
Cebrià, A, Fontanals, M, Martín, P, Morales, JI, Oms, JF, Rodríguez Hidalgo, A, Soto, M, Vergès, JM. 2014. Nuevos datos para el Neolítico antiguo en el nordeste de la Península Ibérica procedentes de la Cova del Toll (Moià, Barcelona) y de la Cueva de la Font Major (L´Espluga de Francolí, Tarragona). Trabajos de Prehistoria 71(1):134–145.CrossRefGoogle Scholar
Collard, M, Edinborough, K, Shennan, S, Thomas, MG. 2010. Radiocarbon evidence indicates that migrants introduced farming to Britain. Journal of Archaeological Science 37:866–870.CrossRefGoogle Scholar
Colledge, S, Connolly, J, Dobney, K, Manning, K, Shennan, S, editors. 2013. The Origins and Spread of Domestic animals in southwest Asia and Europe. Walnut Creek, CA: Left Coast Press.Google Scholar
Cortés, M, Jiménez, FJ, Simón, MD, Gibaja, JF, Carvalho, AF, Martínez-Ruiz, F, Rodrigo, M, Flores, JA, Paytan, A, López, JA, Peña-Chocarro, L, Carrión, JS, Morales, A, Roselló, E, Riquelme, JA, Dean, RM, Salgueiro, E, Martínez, RM, De la Rubia, JJ, Lozano, MC, Vera, JL, Llorente, L, Bicho, NF. 2012. The Mesolithic-Neolithic transition in southern Iberia. Quaternary Research 77(2):221–234.CrossRefGoogle Scholar
Cruz Berrocal, M. 2012. The Early Neolithic in the Iberian Peninsula and the Western Mediterranean: a review of the evidence on migration. Journal of World Prehistory 25:123–156.CrossRefGoogle Scholar
Diaz del Rio, P. 2011. The Neolithic Argonauts of the Western Mediterranean and other underdetermined hypotheses of colonial encounters. In: Bolender D, editor. Eventful Archaeologies New Approaches to Social Transformation in the Archaeological Record. New York: SUNY Press. p 88–99.Google Scholar
Edinborough, K. 2005. Evolution of bow-arrow technology. Doctoral thesis. University of London.Google Scholar
Edinborough, K. 2008. Population history and the evolution of Mesolithic arrowhead technology in south Scandinavia. In: Shennan SJ, editor. Pattern and Process in Cultural Evolution. Berkeley: University of California Press. p 191–202.Google Scholar
Fernández Eraso, J. 2011. Las cerámicas neolíticas de La Rioja alavesa en su contexto. Los casos de Peña Larga y Los Husos I y II. In: Bernabeu Aubán J, Rojo Guerra MA, Molina Balaguer Ll, editors. Las primeras producciones cerámicas: el VI milenio cal AC en la Península Ibérica. Saguntum, PLAV, Extra 12: Universitat de Valencia. p 117–129.Google Scholar
Fernández López de Pablo, J, Gómez Puche, M. 2009. Climate change and population dynamics during the Late Mesolithic and the Neolithic transition in Iberia. Documenta Praehistorica XXXVI:67–96.CrossRefGoogle Scholar
Fernández López de Pablo, J, Salazar-García, DC, Subirà-Galdacano, ME, Roca de Togores, C, Gómez Puche, M, Richards, MP, Esquembre, MA. 2013. Late Mesolithic burials at Casa Corona (Villena, Spain): direct radiocarbon and paleodietary evidence of the last forager populations in Eastern Iberia. Journal of Archaeological Science 40:671–680.CrossRefGoogle Scholar
Forenbaher, S, Miracle, PT. 2006. The spread of farming in the eastern Adriatic. Antiquity 79:514–528.CrossRefGoogle Scholar
Fort, J. 2012. Synthesis between demic and cultural diffusion in the Neolithic transition in Europe. PNAS 109(46):18669–18673.CrossRefGoogle ScholarPubMed
Gamble, C, Davies, W, Pettitt, P, Hazelwood, L, Richards, M. 2005. The archaeological and genetic foundations of the European population during the Late Glacial. Cambridge Archaeological Journal 15:193–223.CrossRefGoogle Scholar
García Borja, P, Salazar García, DC, Martins, H, Pérez Jordà, G, Sanchís Serra, A. 2012. Dataciones radiocarbónicas de la Cova de la Sarsa (Bocairent, Valencia). Recerques del Museu d´Alcoi 21:19–24.Google Scholar
García Borja, P, Aura Tortosa, JE, Jordá Pardo, JE, Salazar García, D. 2014. La cerámica neolítica de la Cueva de Nerja (Málaga, España): salas del Vestíbulo y la Mina. Archivo de Prehistoria Levantina XXX:81–131.Google Scholar
García Puchol, O, Molina, LL, Aura, JE, Bernabeu Aubán, J. 2009. From the Mesolithic to the Neolithic on the Mediterranean coast of the Iberian Peninsula: global action in human context. adapting to the Holocene in Iberia. Journal of Anthropological Research 65(2):237–251.CrossRefGoogle Scholar
Gibaja, JF, Subirà, ME, Terradas, X, Santos, FJ, Agulló, L, Gómez Martínez, I, Allièse, F, Fernández-López de Pablo, J. 2015. The emergence of Mesolithic cemeteries in SW Europe: insights from the El Collado (Oliva, Valencia, Spain) radiocarbon record. PlosOne: 7 January 2015. DOI: 10.1371/journal.pone.0115505.CrossRefGoogle Scholar
Guilaine, J. 2001. La diffusion de l’agriculture en Europe: une hypothèse arythmique. Zephyrus 53–54:267–272.Google Scholar
Guilaine, J, Manen, C. 2007. From Mesolithic to Early Neolithic in the Western Mediterranean. In: Whittle A, Cummings V, editors. Going Over: The Mesolithic-Neolithic Transition in North-West Europe. Oxford: Proceedings of British Academy. p 21–51.Google Scholar
Isern, N, Fort, J, Carvalho, A, Gibaja, J, Ibañez, J. 2014. The Neolithic Transition in the Iberian Peninsula: data analysis and modeling. Journal of Archaeological Method and Theory 21(2):447–460.CrossRefGoogle Scholar
Juan Cabanilles, J, Martí, B. 2002. Poblamiento y procesos culturales en la península Ibérica del VII al V milenio ac. Saguntum: Papeles del Laboratorio de Arqueología de Valencia, Extra 5:45–87.Google Scholar
Juan-Cabanilles, J, García Puchol, O. 2013. Rupture et continuité dans la néolithisation du versant méditerranéen de la péninsule Ibérique: mise à l’épreuve du modèle de dualité culturelle. In: Jaubert, et al. editors. Transitions, ruptures et continuité en Préhistoire. Vol. 1. Évolution des techniques - Comportements funéraires - Néolithique ancien (XXVIIe Congrès Préhistorique de France, Bordeaux - Les Eyzies, 2010). Paris: Société Préhistorique Française. p 405–417.Google Scholar
Manen, C, Sabatier, P. 2003. Chronique radiocarbone de la néolithisation en Méditerranée occidental. Bulletin de la Société préhistorique française 100(3):479–504.CrossRefGoogle Scholar
Manen, C, Marchand, G, Carvalho, AF. 2007. Le Néolithique ancient de la peninsula Ibérique: vers une nouvelle evaluation du mirage africain? Actes du Congress du Centenaire: Un siècle de construction du discours scientifique en Préhistoire 3:133–151.Google Scholar
Martí, B. 2011. La Cova de l’Or (Beniarrés, Alacant). Saguntum. Papeles del Laboratorio de Arqueología de Valencia, nº extra 12:183–186.Google Scholar
Martí, B, JE, Aura, Juan-Cabanilles, J, García Puchol, O, Fernández, J. 2009. El mesolítico Geométrico de tipo «Cocina» en el País Valenciano. In: Utrilla P, Montes L, editors. El Mesolítico Geométrico en la Península Ibérica. Zaragoza: Universidad de Zaragoza. Monografías Arqueológicas 44:205–258.Google Scholar
Martins, H, Oms, FX, Pereira, L, Pike, AW, Rowsell, K, Zilhão, J. 2015. Radiocarbon dating the beginning of the Neolithic in Iberia: new results, new problems. Journal of Mediterranean Archaeology 28(1):105–131.CrossRefGoogle Scholar
McClure, SB. 2013. Domesticated animals and biodiversity: early agriculture at the gates of Europe and long-term ecological consequences. Anthropocene 4:57–68.CrossRefGoogle Scholar
McClure, SB, Podrug, E. 2016. Villages, landscapes and early farming in northern Dalmatia. In: Lillios KT, Chazan M, editors. Fresh Fields and Pastures New: Papers Presented in Honor of Andrew M.T. Moore. Leiden: Sidestone Press. p 117–144.Google Scholar
Medved, I. 2013. Continuty vs. discontinuity. Epipaleolithic and Early Neolithic in the southeast of the Iberian Peninsula. PhD dissertation. Cologne. Available at kups.ub.uni-koeln.de/5217/1/Medved2013_Continuity.pdf.Google Scholar
Montes, L, Domingo, R, González-Sampériz, P, Sebastián, M, Aranbarri, J, Castaños, P, García-Simón, LM, Alcolea, M, Laborda, R. 2016. Landscape, resources and people during the Mesolithic and Neolithic times in NE Iberia: the Arba de Biel Basin. Quaternary International 403:133–150. Available at http://dx.doi.org/10.1016/j.quaint.2015.05.041.CrossRefGoogle Scholar
Navarrete Enciso, MS. 1976. La Cultura de las Cuevas con cerámica decorada en Andalucía Oriental. Granada: Universidad de Granada.Google Scholar
Olalde, I, Schroeder, H, Sandoval-Velasco, M, Vinner, L, Lobón, I, Ramirez, O, Civit, S, García Borja, P, Salazar-García, DC, Talamo, S, Fullola, JM, Oms, FX, Pedro, M, Martínez, P, Sanz, M, Daura, J, Zilhão, J, Marquès-Bonet, T, Thomas, M, Gilbert, P, Lalueza-Fox, C. 2015. A common genetic origin for early farmers from Mediterranean Cardial and Central European LBK cultures. Molecular Biology and Evolution: September 2015. DOI: 10.1093/molbev/msv181.CrossRefGoogle Scholar
Oms, FX, Esteve, X, Mestres, J, Martín, P, Martins, YH. 2014. La neolitización del nordeste de la Península Ibérica: datos radiocarbónicos y culturales de los asentamientos al aire libre del Penedès. Trabajos de Prehistoria 71(1):42–55 DOI: 10.3989/tp.2014.12123.CrossRefGoogle Scholar
Özdoğan, M. 2011. Archaeological evidence on the westward expansion of farming communities from Earstern Anatolia to the Aegean and the Balkans. Current Anthropology 54(Supplement 4):S415–S430.CrossRefGoogle Scholar
Peña Chocarro, L, Pérez Jordà, G, Morales Mateos, J, Vera Rodríguez, JC. 2013. Y llegaron los agricultores: agrigultura y recolección en el occidente del mediterráneo. Menga. Revista de Prehistoria de Andalucía 4:15–33.Google Scholar
Perlès, C. 2001. The Early Neolithic in Greece: the First Farming Communities in Europe. Cambridge: Cambridge University Press.CrossRefGoogle Scholar
Reimer, PJ, Bard, E, Bayliss, A, Beck, JW, Blackwell, PG, Bronk Ramsey, C, Buck, CE, Cheng, H, Edwards, RL, Friedrich, M, Grootes, PM, Guilderson, TP, Haflidason, H, Hajdas, I, Hatté, C, Heaton, TJ, Hoffmann, DL, Hogg, AG, Hughen, KA, Kaiser, KF, Kromer, B, Manning, SW, Niu, M, Reimer, RW, Richards, DA, Scott, EM, Southon, JR, Staff, RA, Turney, CSM, van der Plicht, J. 2013. IntCal13 and Marine13 radiocarbon age calibration curves 0–50,000 years cal BP. Radiocarbon 55(4):1869–1887.CrossRefGoogle Scholar
Riede, F, Edinborough, K. 2012. Bayesian radiocarbon models for the cultural transition during the Allerød in southern Scandinavia. Journal of Archaeological Science 3:744–756.CrossRefGoogle Scholar
Robb, J. 2007. The Early Mediterranean Village: Agency, Material Culture, and Social Change in Neolithic Italy. Cambridge University Press.CrossRefGoogle Scholar
Rojo, MA, Garrido, R, García-Martínez de Lagrán, I, editors. 2012. El neolítico en la Península Ibérica y su contexto europeo. Madrid: Cátedra.Google Scholar
Rojo Guerra, MA, Peña Chocarro, L, Royo Guillen, JI, Tejedor Rodríguez, C, García Martínez de Lagrán, I, Arcusa Magallón, H, Garrido Peña, R, Moreno García, M, Mazzuco, N, Gibaja Bao, JF, Ortega, D, Kromer, B, Kurt, WA. 2013. Pastores trashumantes del neolítico antiguo en un entorno de alta montaña: secuencia crono-cultural de la Cova de els Trocs (San Feliú de Veri, Huesca). BSAA Arqueología LXXIX:9–55.Google Scholar
Rowley-Conwy, P. 2011. Westward ho! The spread of agriculture from central Europe to the Atlantic. Current Anthropology 52(4):431–451.CrossRefGoogle Scholar
Salazar García, DC, Aura, JE, Olària, CR, Talamo, S, Morales, JV, Richards, MP. 2014. Isotope evidence for the use of marine resources in the Eastern Iberia Mesolithic. Journal of Archaeological Science 42:231–240.CrossRefGoogle Scholar
Shennan, S, Edinborough, K. 2007. Prehistoric population history: from the Late Glacial to the Late Neolithic in central and northern Europe. Journal of Archaeological Science 34:1339–1345.CrossRefGoogle Scholar
Shennan, S, Downey, SS, Timpson, A, Edinborough, K, Colledge, S, Kerig, T, Manning, K, Thomas, MG. 2013. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nature Communications 4:2486.CrossRefGoogle ScholarPubMed
Soler Diaz, JA, García Atienzar, G, Duque Espino, D, Ferrer García, C, Gómez Pérez, O, Guillem Calatayud, PM, Iborra Eres, P, Martínez Valle, R, Pérez Jordà, G, Roca de Togores Muñoz, C, Ximénez de Embún Sánchez, T. 2013. La primera ocupación neolítica de la Cova d´En Pardo (Planes, Alicante). Avance de estudio pluridisciplinar de los niveles VIII y VIIIb. Universitat de València. Saguntum PLAV 45(9):24.Google Scholar
Timpson, A, Colledge, S, Crema, E, Edinborough, K, Kering, T, Manning, K, Thomas, MG, Shennan, S. 2014. Reconstructing regional population fluctuations in the European Neolithic using radiocarbon dates: a new case-study using an improved method. Journal of Archaeological Science 52:549–557.CrossRefGoogle Scholar
Utrilla, P, Montes, L, Mazo, C, Martínez Bea, M, Domingo, R. 2009. El Mesolítico geométrico en Aragón. In: Utrilla P, Montes L, editors. El Mesolítico Geométrico en la Península Ibérica. Zaragoza: Universidad de Zaragoza. Monografías Arqueológicas 44 131–190.Google Scholar
Utrilla, P, Berdejo, A, Obón, A. 2012. El Esplugón: un gran abrigo mesolítico en el valle del Guarga (Huesca). In: Muñiz Álvarez JR, coordinator. AD ORIENTEM. Del final del Paleolítico en el norte de España a las primeras civilizaciones del Oriente Próximo. Universidad de Oviedo: Ménsula Ediciones. p 235–251.Google Scholar
Vaquero Rodríguez, M, García-Argüelles i Andreu, P. 2009. Algunas reflexiones sobre la ausencia de mesolítico geométrico en Catalunañ. In: Utrilla P, Montes L, editors. El Mesolítico Geométrico en la Península Ibérica. Zaragoza: Universidad de Zaragoza. Monografías Arqueológicas 44 191–203.Google Scholar
Whittle, AWR, Healy, FMA, Bayliss, A. 2011. Gatherting Time: Dating the Early Neolithic Enclosures of Southern Britain and Ireland. Oxford: Oxbow Books.CrossRefGoogle Scholar
Williams, AN. 2012. The use of summed radiocarbon probability distributions in archaeology: a review of methods. Journal of Archaeological Science 39(6):578–589.CrossRefGoogle Scholar
Zilhão, J. 2001. Radiocarbon evidence for maritime pioneer colonization at the origins of farming in west Mediterranean Europe. Proceedings of the National Academy of Sciences 98:14180–14185.CrossRefGoogle ScholarPubMed
Zilhão, J. 2003. The Neolithic Transition in Portugal and the role of demic diffusion in the spread of agriculture across West Mediterranean Europe. In: Ammerman AJ, Biagi P, editors. The Widening Harvest. The Neolithic Transition in Europe: Looking Back, Looking Forward. Boston: Archaeological Institute of America. p 207–223.Google Scholar
Zilhão, J. 2011. Time is on my side. In: Hadjikoumis A, Robinson E, Viner S, editors. The Dynamics of Neolithisation in Europe: Studies in Honour of Andrew Sherratt. Oxford: Oxbow Books. p 46–65.CrossRefGoogle Scholar
Zvelebil, M, Rowley Conwy, P. 1986. Foragers and farmers in Atlantic Europe. In: Zvelebil M, editor. Hunters in Transition: Mesolithic Societies of Temperate Eurasia and Their Transition to Farming. Cambridge University Press. p 67–93.Google Scholar
Figure 0

Figure 1 Map of Iberia showing the regions discussed.

Figure 1

Table 1 Radiocarbon dates (short-life, singular samples) used to build the chronological models. A: area, PH: phase, GM: Geometric Mesolithic; EN: Early Neolithic.

Figure 2

Table 2 Bayesian sequence modeled dates from Oxcal 4.2.3 (Bronk Ramsey 2009). Calibration conducted using IntCal13 curve: Northern Hemisphere (Reimer et al. 2013).

Figure 3

Figure 2 Mesolithic and Neolithic sites in Mediterranean Iberia region used for chronological modeling.

Figure 4

Figure 3 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Ebro valley region.

Figure 5

Figure 4 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Eastern area (model A).

Figure 6

Figure 5 Phase modeling results showing the probability distributions of radiocarbon dates for Mesolithic and Neolithic overlapping phases in the Eastern area (model B).

Figure 7

Figure 6 Phase modeling results showing the end of the Mesolithic and the beginning of the Neolithic according with the different models proposed (general A model, and regional models: Eastern A and B models, and Ebro model).

Figure 8

Table 3 Bayesian sequence modeled dates from Oxcal 4.2.3 (Bronk Ramsey 2009). Interval span and distance between the end of the Mesolithic and the start of the Neolithic. Calibration conducted using IntCal13 curve: Northern Hemisphere (Reimer et al. 2013).

Figure 9

Figure 7 Comparative sumprob radiocarbon dates/sites distributed within intervals of 200 years. We calibrated all the dates using Oxcal Program 4.2 and IntCal13 curve for the Northern Hemisphere (Reimer et al. 2013) with 1σ (68.2).

Supplementary material: PDF

García-Puchol supplementary material

García-Puchol supplementary material

Download García-Puchol supplementary material(PDF)
PDF 928.5 KB