Impact statement
Present in Europe for more than 700,000 years, there has long been debate concerning the presence and loss of Acheulean Palaeolithic culture on the continent. Attention has often focused on the role of glaciation and demographic factors in northern and southern regions. Here, the temporal cohesion of the European Acheulean archaeological record is statistically assessed under multiple spatial scenarios. Few breaks in the archaeological record are identified, suggesting Acheulean cultural information to have only occasionally been regionally absent. Four absences appear linked to northern glacial cycles, while a fifth is observed in Iberia soon after the Acheulean’s introduction to Europe c. 880,000 years ago. This study represents the first to assess an exhaustive database of reliably dated European Acheulean sites in the pursuit of identifying cultural and demographic patterns during this pivotal point in the early colonisation of Europe. At a continental level, the Acheulean appears to have constantly been present in Europe after its first introduction, demonstrating the importance of these technologies to hominin populations and its durability as a cultural phenomenon. These data have implications for understanding the persistence of Acheulean culture in Africa and elsewhere in Eurasia across its c. 1.5 million years.
Introduction
Understanding the spatial and temporal character of the archaeological record is a fundamental goal of archaeologists. Diverse processes determine when and where we observe past human material culture, but by gaining an accurate picture of the artefact record it becomes easier to identify these cultural evolutionary, demographic, ecological and geological influences, among others. In the case of prehistoric humans, archaeologists are faced with identifying these varied and dynamic processes using a sparse and highly fragmented archaeological record (Isaac, Reference Isaac1969; Binford, Reference Binford and Soffer1987; Lycett and Eren, Reference Lycett and Eren2013; Kuhn and Clark, Reference Kuhn and Clark2015; French, Reference French2016; Gallotti, Reference Gallotti2016; Pope et al., Reference Pope, Blundell, Scott and Cutler2016; Kuhn, Reference Kuhn2020; Key et al., Reference Key, Roberts and Jarić2021).
Europe displays perhaps the richest record of dated prehistoric sites in the world. As such, it provides an opportunity to gain insight into large, continental-scale cultural processes at resolutions that may be unachievable elsewhere. Nowhere is this more apparent than when considering the Lower Palaeolithic/Early Stone Age (ESA), given more than 160 years of research on the continent (Prestwich, Reference Prestwich1860; Evans, 1872; White, Reference White2022). The Acheulean represents the most heavily studied of European Lower Palaeolithic cultural phenomena and covers a majority of the period from c. 900 to 150 thousand years ago (ka) (Moncel et al., Reference Moncel, Santagata, Pereira, Nomade, Voinchet, Bahain, Daujeard, Curci, Lemorini, Hardy, Eramo, Berto, Raynal, Arzarello, Mecozzi, Iannucci, Sardella, Allegretta, Delluniversita, Terzano, Dugas, Jouanic, Queffelec, d’Andrea, Valentini, Minucci, Carpentiero and Piperno2020c; Ashton and Davis, Reference Ashton and Davis2021; Key et al., Reference Key, Roberts and Jarić2021; Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023). As such, it provides a relatively rich record of sites, many of which have been robustly dated using modern radiometric techniques (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a).
The Acheulean replaces the purely flake-and-core focused technologies observed in earlier, more sporadically evidenced European populations, likely representing a dispersal of new cultural information from western/central Eurasia and, potentially, Africa (Dennell and Roebroeks, Reference Dennell and Roebroeks1996; Sharon, Reference Sharon2011; Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; Gallotti, Reference Gallotti2016; Sharon and Barsky, Reference Sharon and Barsky2016; Arroyo et al., Reference Arroyo, Proffitt and Key2019; Méndez-Quintas et al., Reference Méndez-Quintas, Santonja, Arnold, Cunha-Ribeiro, da Silva, Demuro, Duval, Gomes, Meireles, Monteiri-Rodrigues and Pérez-González2020; Moncel et al., Reference Moncel, Despriee, Courcimaut, Voinchet and Bahain2020b,Reference Moncel, Santagata, Pereira, Nomade, Voinchet, Bahain, Daujeard, Curci, Lemorini, Hardy, Eramo, Berto, Raynal, Arzarello, Mecozzi, Iannucci, Sardella, Allegretta, Delluniversita, Terzano, Dugas, Jouanic, Queffelec, d’Andrea, Valentini, Minucci, Carpentiero and Pipernoc. The tradition persists until Neanderthals and later prepared core technologies appear from c. 400 to 300 ka (Arsuaga et al., Reference Arsuaga, Martinez, Arnold, Aranburu, Gracia-Tellez, Sharp, Quam, Falgueres, Pantoja-Perez, Bischoff, Poza-Rey, Pares, Carretero, Demuro, Lorenzo, Sala, Martinon-Torres, Garcia, Alcazar, Cuenca-Bescos, Gomez-Olivencia, Moreno, Pablos, Shen, Rodriguez, Ortega, Garcia, Bonmati, Bermudez de Castro and Carbonell2014; Ollé et al., Reference Ollé, Verges, Rodriguez-Alvarez, Caceres, Angelucci, Vallverdu, Demuro, Arnold, Falgueres, Bennasar, Lopez-Garcia, Blain, Banuls-Cardona, Burjachs, Exposito, Lopez-Polin and Lopez-Ortega2016b; Moncel et al., Reference Moncel, Ashton, Arzarello, Fontana, Lamotte, Scott, Muttillo, Berruti, Nenzioni, Tuffreau and Peretto2020a; Key et al., Reference Key, Roberts and Jarić2021), with many of the youngest known Acheulean sites being observed in southern France and Iberia (Michel et al., Reference Michel, Shen, Valensi and de Lumley2009; Monteiro-Rodrigues and Cunha-Ribeiro, Reference Monteiro-Rodrigues and Cunha-Ribeiro2014; Méndez-Quintas et al., Reference Méndez-Quintas, Demuro, Arnold, Duval, Perez-Gonzalez and Santonja2019) (Table 1, Figure 1).
Note: The full database of retrieved information for each site is available in Supplementary material S1. This includes citations to the relevant articles from which data were retrieved. Dates are in years before present.
Morphological, technological, temporal and spatial evidence points to a single, but variable, cultural tradition being represented by the Acheulean (Gowlett, Reference Gowlett1979; Lycett and Gowlett, Reference Lycett and Gowlett2008; Shipton, Reference Shipton and Groucutt2020; Key, Reference Key2022), which itself is most often defined by the presence of bifacially flaked core technologies (Sharon, Reference Sharon2010; de la Torre and Mora, Reference de la Torre, Mora, Cole, McNabb, Grove and Hosfield2020; Kuhn, Reference Kuhn2020). Within Europe, two forms of bifaces – handaxes and cleavers – are produced, although each varies within and between assemblages, and variants such as picks and ficrons have been defined (Santonja and Villa, 1990; Wymer, Reference Wymer1999; Vaughan, Reference Vaughan, Hurt and Rakita2001; Lycett and Gowlett, Reference Lycett and Gowlett2008; Emery, Reference Emery2010; Key, Reference Key2019; Méndez-Quintas et al., Reference Méndez-Quintas, Santonja, Arnold, Cunha-Ribeiro, da Silva, Demuro, Duval, Gomes, Meireles, Monteiri-Rodrigues and Pérez-González2020; McNabb, Reference McNabb2022; García-Medrano et al., Reference García-Medrano, Moncel, Maldonado-Garrido, Ollé and Ashton2023).
Acheulean cultural information is not ubiquitous in Europe after c. 900 ka. Its absence from eastern and central Europe has long been known (Klein, Reference Klein1966; Dennell and Roebroeks, Reference Dennell and Roebroeks1996; Rocca et al., Reference Rocca, Abruzzese and Aureli2016; Sharon and Barsky, Reference Sharon and Barsky2016), potentially due to the influence of climatic factors, including low temperatures and low precipitation, on ecology (Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review). Climatically linked cycles of Acheulean presence and absence have also long been proposed in northern Europe (Roe, Reference Roe1981; Wymer, Reference Wymer1999; White and Schreve, Reference White and Schreve2000). During warmer interglacial periods, populations with Acheulean culture have been suggested to occupy northwestern Europe, only to be driven out during colder glacial periods (Ashton and Lewis, Reference Ashton and Lewis2012; Moncel et al., Reference Moncel, Ashton, Lamotte, Tuffreau, Cliquet and Despriee2015; Shipton and White, Reference Shipton and White2020; Ashton and Davis, Reference Ashton and Davis2021). As evidenced through the warm marine isotope stages (MIS) associated with nearly all biface sites in the region (Table 1; Supplementary material S1) (although see Moncel et al., Reference Moncel, Antoine, Herisson, Locht, Hurel and Bahain2022).
Other episodes of Acheulean cultural loss – or extirpation (localised, regional extinction) – have been proposed in southern Europe (MacDonald et al., Reference MacDonald, Martinon-Torres, Dennell and de Castro2012). The recent discovery of early, temporally outlying, biface sites in Iberia have created a substantial gap in the region’s Acheulean record between c. 850 and 500 ka (Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; Walker et al., Reference Walker, Uriarte, Jimenez, Martinez, Lerma, Van der Made, Duval and Grun2020; Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023). The Italian peninsula has evidence of sparsity in its Acheulean record too, with early sites such as Notarchirico (Moncel et al., Reference Moncel, Santagata, Pereira, Nomade, Bahain, Voinchet and Piperno2019) and Valle Giumentina (Villa et al., Reference Villa, Nicoud, Guibert-Cardin, Tomasso, Chausse, Boschian, Degeai, Fusco and Limondin-Lozouet2024) evidencing an 80,000-year gap to later occurrences such as Fontana Ranuccio (Muttoni et al., Reference Muttoni, Scardia, Kent, Swisher and Manzi2009). Elsewhere in Europe, temporal breaks of 30,000 years or more are evidenced in the Acheulean record (Table 1; Supplementary Tables S1 and S2). These breaks do not necessitate an absence of hominin populations – as in the case of the northern ‘Clactonian’ (Ashton and Davis, Reference Ashton and Davis2021) or the diverse Middle Pleistocene flake and core sites in southern Europe (Martínez and Garriga, Reference Martínez and GArriga2016) – but instead a loss of populations retaining the cultural information required for the production of bifaces (Lycett and von Cramon-Taubadel, Reference Lycett and von Cramon-Taubadel2008).
Our ability to understand why there may have been regional or continental-level breaks in the archaeological record is, however, dependent on gaining an accurate picture of where these gaps occur. We may infer an Acheulean absence based on a 50,000-year gap in the archaeological record, but search intensity biases, taphonomic processes and past demographic variation, among other factors, could all have plausibly created the perception of a gap, when in reality the cultural information was present (Surovell et al., Reference Surovell, Finley, Smith, Brantingham and Kelly2009; Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a; Pope et al., Reference Pope, Blundell, Scott and Cutler2016; Key and Ashton, Reference Key and Ashton2023). Even when these processes were equal, archaeologists may infer a cultural absence simply because a temporal break is subjectively perceived to be large. In 2005, Solow and Smith (Reference Solow and Smith2005) introduced the ‘surprise test’ to Palaeolithic archaeology, a statistical method capable of assessing the temporal exceptionality of an outlying occurrence (dated site) relative to a sample or earlier or later occurrences. The technique assesses the scale of a break in the known archaeological record relative to the sites preceding or following it, and objectively records how likely it is to represent an absence of the cultural phenomenon under investigation.
Following Solow and Smith (Reference Solow and Smith2005) and Roberts et al. (Reference Roberts, Jarić, Lycett, Flicker and Key2023), Key (Reference Key2022) used the surprise test to demonstrate the early and late Acheulean records of Africa and Eurasia to be temporally cohesive. This included Europe, where no significant breaks in the Acheulean record were identified between 300 and 160 ka (Key, Reference Key2022). In turn, it became possible to infer that during this time a continuous lineage of Acheulean cultural information was likely present on the continent. Here, the temporal cohesiveness of the entire European Acheulean archaeological record is statistically assessed. Using a comprehensive sample of reliably dated biface-retaining sites, gathered from an exhaustive review of published literature, the relative scale of regional and continental-level breaks in the continent’s entire Acheulean record is examined. Significant breaks and long periods of continuity are observed, providing new insight into the loss (extirpation) and persistence of Acheulean cultural information in Europe.
Methods
Following Solow and Smith (Reference Solow and Smith2005), the ‘surprise test’ was used to identify whether temporal gaps in the European Acheulean archaeological record should be considered representative of cultural absence. The surprise test asks whether a new, potentially outlying record was generated by the same process that created previous or later consecutive records (Solow and Smith, Reference Solow and Smith2005). In the present context, it asks whether a dated Acheulean occurrence (site) can be considered part of the same lineage of cultural information that preceded or followed it, or alternatively, whether it represents a culturally distinct, temporally ‘surprising’, occurrence (Key, Reference Key2022; Roberts et al., Reference Roberts, Jarić, Lycett, Flicker and Key2023). Rejection of the null hypothesis – cultural continuity between the preceding or following occurrences and the occurrence of interest – indicates a relative temporal gap sufficient to infer cultural absence.
The surprise test uses range and spacing data across a series of k consecutive temporal occurrences – here, dated Acheulean sites. The record of interest could be an outlier site or a site at the start of a new series of consecutive occurrences, preceded by a temporal gap of any scale. The occurrences (sample) against which the record of interest is tested are assumed to represent the k largest or smallest records of a larger collection of records generated from a distribution from the Gumbel domain of attraction (Solow and Smith, Reference Solow and Smith2005). The Gumbel distribution can be used to fit diverse scenarios, including those characterised by symmetrical, skewed, unimodal and bimodal data (Al-Aqtash et al., Reference Al-Aqtash, Lee and Famoye2014).
As a generalised extreme value distribution, the Gumbel distribution can be used to model the range limits of scaled linear data. In this case, years before present, represented by the age of dated Acheulean sites, with the youngest or oldest records in this sample feasibly representing the start or end of a lineage of cultural information. If the record of interest is identified as being statistically surprising relative to the larger sample, it can be considered part of a separate lineage of cultural information. Thus, the temporal gap evidencing this cultural distinction could represent Acheulean absence and subsequent re-emergence. As outlined by Key (Reference Key2022), in Europe this would not necessarily represent an episode of cultural convergence, but most likely implies an extirpation event followed by the Acheulean’s later reintroduction from elsewhere in Eurasia or Africa (Figure 3). If the record of interest is not statistically exceptional or surprising relative to the main site sample, then persistence of the Acheulean across the investigated temporal gap can be supported.
Surprise test
As described by Solow and Smith (Reference Solow and Smith2005), and more recently by Roberts et al. (Reference Roberts, Jarić, Lycett, Flicker and Key2023), for tests in the forward temporal direction, let t 1 > t 2 > … > t k be the k most recent Acheulean records ordered from the most recent to the earliest. With the record of interest being dated at time y, the test assesses the exceptionality of this more recent occurrence. Following the null hypothesis that the later record of interest was generated by the same process as the earlier occurrences (i.e. the main sample), Solow and Smith (Reference Solow and Smith2005) demonstrated the quantity,
to have a β distribution with parameters 1 and k–1 so that the P-value corresponding to an observed value Sk is
Solow and Smith (Reference Solow and Smith2005) demonstrate that the power of the surprise test does not heavily depend on k, with k of 5 and 10 performing adequately; both are applied here. Due to the finite record of Acheulean occurrences, both forward and reverse versions of the model are also used when possible (Table 2). Thus, both ‘origination’ and ‘extinction’ Gumbel distribution tails are modelled, dependent on whether the test is run in the reverse or forward temporal direction (respectively). Alter the above instruction as appropriate in the case of tests in the reverse temporal direction (see Supplementary material S2). In all instances α = 0.05.
Note: See Supplementary Tables S3–S7 for the results pertaining to each investigated temporal gap. A clear inference of Acheulean absence (i.e. ‘yes’) is based on a significant result being returned and a majority of the sites in the main sample not displaying identical (or near-identical) ages assigned principally through MIS-stage associations only. When the latter criteria cannot be met, the absence of cultural information is determined to be ‘unclear’. Note that only the Iberia scenario returned a significant result via the resampling procedure.
As all dated European Acheulean sites are represented by age ranges (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a), the resampling procedure applied by Roberts et al. (Reference Roberts, Jarić, Lycett, Flicker and Key2023) was followed here to account for this uncertainty. Dates were drawn randomly from a normal distribution bounded by the defined age range for a given occurrence. The central age (see below) of each Acheulean occurrence was used as the mean value, while the standard deviation was half the difference between the central age and the relevant range boundary. These randomly generated datasets were investigated using the surprise test as outlined above, with the process being repeated 5,000 times. The mean across all iterations was used as the resampling result. The resampling procedure was used in addition to running the surprise test using each occurrence’s central age value. R version 4.3.2 was used throughout (R Core Team, 2013). Associated code is available in Supplementary material S2 and Roberts et al. (Reference Roberts, Jarić, Lycett, Flicker and Key2023).
European Acheulean site sample and data scenarios
An exhaustive review of Acheulean sites in Europe was undertaken (Figure 1). European Acheulean sites were spatially defined as belonging to the European continent (or outlying islands) up to the western borders of modern Russia and Turkey. Technologically, sites were only included if they displayed the presence of bifacially flaked core-tools (handaxes or cleavers) and were associated with the Acheulean tradition by the individuals who excavated and/or dated the site (i.e. those who know the site best; Supplementary material S1). Sites or archaeological layers described as Acheulean but also displaying prepared core technologies were excluded due to the presence of Middle Palaeolithic-defining cultural information.
The surprise test procedure, as outlined above, required three pieces of temporal data for each occurrence: a central age, as well as upper and lower range boundaries. Sites without these data were excluded from the sample, as were all sites with unreliable age determinations. Individual sites could return multiple occurrences for inclusion, so long as no date-range overlap was observed between each Acheulean layer (e.g. la Noira, France [Moncel et al., Reference Moncel, Despise, Voinchet, Tissoux, Moreno, Bahain, Courcimault and Falgueres2013]).
Age determination reliability was graded between three and zero for each site. Three represented a securely dated site, while zero represented a site with age associations that could potentially negatively impact model accuracy. Sites graded ‘zero’ are not listed in Table 1. Central age values often represented the author’s ‘preferred’ site age (Supplementary material S1); sometimes reported as a date range’s central value, or a centralised date determined based on sediment accumulation rates or other evidence. Date ranges were almost exclusively determined using radiometric methods. Occasionally, dating procedures resulted in marine isotope stage (MIS) age associations. In these instances, MIS boundaries were used for date ranges, following Lisiecki and Raymo (Reference Lisiecki and Raymo2005). Data reflect current understanding in September 2023, but it is important to recognise that some age determinations are subject to ongoing debate and/or research. Additional justification for the inclusion or exclusion of specific sites can be seen in Supplementary material S1. The site review was intended to be exhaustive, but it is acknowledged that a small number of sites could have been missed. Middle Pleistocene archaeological sites often display poor chronological resolution and unclear technological comparability between sites (MacDonald and Roebroeks, Reference MacDonald and Roebroeks2012), but the present investigation represents an analysis of the field’s current understanding. As new data and refined understanding come to light the analyses should be repeated.
Longitudinal and latitudinal data were recorded for each site to facilitate investigation of five spatially defined data scenarios. These scenarios reflect current understanding concerning the presence and absence of Acheulean culture in Europe.
Scenario 1 (S1): Europe S1. Every securely dated Acheulean site in Europe (n = 67). This represents all sites assigned a date reliability value of three. Sites with reliability values from two to zero were excluded. S1 examines whether the Acheulean was ever absent from Europe after its earliest known presence at Barranc de la Boella (Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023).
Scenario 2 (S2): Europe S2. Every dated Acheulean site in Europe (n = 82; Table 1). This represents all sites assigned a date reliability value from three to one. Sites with reliability values of zero were excluded. S2 similarly examines whether the Acheulean was ever absent from Europe after its earliest known presence.
Scenario 3 (S3): Northern Europe. Defined as all sites above 49° latitude with a date reliability of three to one (n = 36). This scenario investigates the widely held view that Acheulean hominins were repeatedly forced from northern latitudes during glacial periods due to inhospitable climatic conditions. The most southerly located site in this sample is St Pierre-les-Elbeuf (Cliquet et al., Reference Cliquet, Lautridou, Antoine, Lamothe, Leroyer, Limondin-Lozouet and Mercier2009).
Scenario 4 (S4): Southern Europe. Defined as all sites below 45° latitude with a date reliability of three to one (n = 38). In this scenario, Barbas 1 (Boëda et al., Reference Boëda, Kervazo, Mercier and Valladas1996) represents the most northerly site included in the sample. S4 examines Acheulean cultural continuity across southern Europe, on the basis that it is widely thought of as habitable by hominins across glacial and interglacial periods.
Scenario 5 (S5): Iberia. Defined as all sites south and west of La Cansaladeta (Ollé et al., Reference Ollé, Verges, Rodriguez-Alvarez, Caceres, Angelucci, Vallverdu, Demuro, Arnold, Falgueres, Bennasar, Lopez-Garcia, Blain, Banuls-Cardona, Burjachs, Exposito, Lopez-Polin and Lopez-Ortega2016b) and Atapuerca (García-Medrano et al., Reference García-Medrano, Ollé, Mosquera, Caceres, Diez and Carbonell2014) on the Iberian Peninsula with a date reliability of three to one (n = 21). This scenario investigated Acheulean cultural continuity at a localised, regional level due to the peninsula’s geographic isolation and its potential role as a refugium for Acheulean populations.
In each spatial scenario, the model tests the null hypothesis that temporal gaps in the archaeological record are not the result of the loss of Acheulean cultural information, be it at a localised (S3, S4, S5) or continental (S1, S2) level. While these regional categorisations are artificial and come with inherent inferential limitations (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a), they reflect a latitudinal and geographic reality that would have impacted demographic processes, and aid current understanding given the field often focuses on northern, southern and Iberian spatial scenarios.
Results
Table 2 displays all significant results across all spatial scenarios and versions of the model. From a total of 239 investigated temporal gaps (S1, S2, S3, S4 and S5 had 66, 81, 35, 37 and 20 gaps, respectively), 18 potential instances of Acheulean cultural absence were identified (Supplementary Tables S3 to S7). Some, such as the gap between Abbeville (France) and Old Park (UK), were significant in all relevant spatial scenarios (Supplementary Tables S3 to S5). Others, such as the temporal gap succeeding the site of Cueva Negra (Spain), were only significant in one scenario; in this case the Iberian Peninsula scenario (Supplementary Table S7). The European Acheulean archaeological record is, therefore, relatively cohesive with few periods when Acheulean cultural information may be lacking at a regional level.
Of the 18 significant results, only five can be considered reliable indicators of cultural absence (Table 2). This is due to the dating approaches used in northern Europe, where ESR, OSL, IRFR and other radiometric techniques are used to determine sediment/site ages, and often these are subsequently used to associate artefacts with an interglacial MIS stage (e.g. Bridgland, Reference Bridgland1994; Antoine et al., Reference Antoine, Moncel, Locht, Limondin-Lozouet, Auguste, Stoetzel, Dabkowski, Voinchet, Bahain and Falgueres2015; Davis et al., Reference Davis, Ashton, Hoare and LEwis2021). This occasionally results in several sites with radiometrically determined, but identical, MIS stage age ranges and central age estimates clustering together. When three or more of these sites cluster, a significant result can be returned even when the investigated temporal gap is relatively small. A phenomenon more likely to occur when k = 5. By combining the results in Table 2 with the dating techniques and age associations in Supplementary material S1, these instances can be identified. When this is considered, 13 significant temporal gaps are revealed to be the product of this phenomenon, meaning only five can be considered a reliable (clear) indicator of Acheulean absence (Table 2). A ‘clear’ inference of Acheulean absence is based on (1) a significant result and (2) a majority of the relevant main sample sites not displaying identical (or near-identical) ages assigned through MIS-stage associations. Only one clear instance of cultural absence included a significant resampling result, emphasising how our understanding of the phenomenon’s temporal character is limited by the date ranges associated with many sites.
Of these five instances, the 423,000-year gap between Cueva Negra and Sima de los Huesos on the Iberian Peninsula (S5) is the largest and clearest period of Acheulean absence. All four model versions were significant (Table 2), strongly supporting the inference that, based on current evidence, Acheulean culture was not present in Iberia during this time. The northern European scenario (S3) returned the four other significant and seemingly reliable periods of Acheulean absence. The earliest, between Moulin Quignon (660 ka) and Rampart Field (592 ka), broadly aligns with the cold glacial MIS 16. This is followed by the period between Maidscross Hill (580 ka) and Abbeville (525 ka), which aligns with the cold MIS 14. High Lodge (492 ka) and Beeches Pit (414 ka) bound the next significant period of Acheulean absence in northern Europe, which can be associated with MIS 12. Finally, a period of biface absence is inferred between St Pierre-les-Elbeuf (385 ka) and Stoke Newington (318 ka), another glacial stage; in this case MIS 10. The current middle Pleistocene archaeological record of northern Europe therefore supports repeated cycles of Acheulean presence and absence in-line with interglacial and glacial marine isotope stages.
Discussion
The European Acheulean was overwhelmingly a story of cultural persistence. Only five regionally-defined breaks in the European archaeological record are great enough, on a relative basis, to reliably infer a period of Acheulean cultural absence. The majority of these appear to have been driven by glacial cycles in northern Europe, where current site-dating suggests the region became too cold for populations with Acheulean culture to survive during MIS 16, 14, 12 and 10 (although see counter argument below). A substantial break in the Acheulean record is evidenced in Iberia between c. 850 and 500 ka. At a continental level, neither investigated scenario (S1, S2) returned a reliable significant result, suggesting Acheulean cultural information to have been permanently present on the continent after its first introduction. The Acheulean cultural phenomenon appears to have only ended at a continental level once Middle Palaeolithic (e.g. Levallois) technologies start to emerge in Europe (Moncel et al., Reference Moncel, Despriee, Courcimaut, Voinchet and Bahain2020b; Key et al., Reference Key, Roberts and Jarić2021), potentially due to functional and economic (raw material) advantages (Brantingham and Kuhn, Reference Brantingham and Kuhn2001; Lycett and Eren, Reference Lycett, Eren, Vonk and Shackelford2022).
Discussion of trends in the European Lower Palaeolithic must acknowledge the diverse natural (e.g. geological, taphonomic) and human-led (e.g. search intensity biases, variation in funding) factors impacting where and when we see evidence of hominin populations (Surovell et al., Reference Surovell, Finley, Smith, Brantingham and Kelly2009; MacDonald and Roebroeks, Reference MacDonald and Roebroeks2012; Pope et al., Reference Pope, Blundell, Scott and Cutler2016; Key and Ashton, Reference Key and Ashton2023). Undoubtedly, present trends will vary in some ways relative to those realised in the Middle Pleistocene. Nonetheless, more than 160 years of archaeological discovery has informed these analyses and many trends will be correct, particularly for the denser portions of the Acheulean record, and as with the present study, archaeologists have a wealth of statistical means at their disposal to help navigate such challenges (e.g. Surovell et al., Reference Surovell, Finley, Smith, Brantingham and Kelly2009; Faith et al., Reference Faith, Du, Behrensmeyer, Davies, Patterson, Rowan and Wood2021; Key et al., Reference Key, Roberts and Jarić2021; Vidal-Cordasco et al., Reference Vidal-Cordasco, Ocio, Hickler and Marin-Arroyo2022). New archaeological sites, new dating methods and revised dating efforts are also constantly updating Acheulean temporal records. Future reanalysis is therefore encouraged, and while results may vary from the present data, these differences should be slight (assuming all test assumptions are met equally).
Regional Acheulean extinction(s)
Iberia
The present data provide empirical support for longstanding inferences concerning the loss of Acheulean cultural information in Iberia and northern Europe. In the former, a notable break in the Acheulean record has been evident since the discovery of Barranc de la Boella and Cueva Negra, two temporally outlying Acheulean sites in Spain (Walker et al., Reference Walker, Uriarte, Jimenez, Martinez, Lerma, Van der Made, Duval and Grun2020; Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023). Indeed, a c. 423,000-year break in the archaeological record is substantial; greater even than the remainder of the Iberian Acheulean after its reappearance at Sima de los Huesos (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a). Potential explanation for this absence includes an early but short-lived Acheulean dispersal event from North Africa (Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023), the extinction of H. antecessor and any associated cultural information c. 800 ka (Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013), and geological and/or research biases impacting the known archaeological record (Vallverdu et al., Reference Vallverdu, Saladie, Rosas, Huguet, Caceres, Mosquera, Garci-Tabernero, Estalrrich, Lozano-Fernandez, Pineda-Alcala, Carrancho, Villalain, Bourles, Braucher, Lebatard, Vilalta, Esteban-Nadal, Bennasar, Bastir, Lopez-Polin, Olle, Verges, Ros-Montoya, Martinez-Navarro, Grcia, MArtinell, Exposito, Burjachs, Agusti and Carbonell2014). A lack of ecological (utilitarian) selective pressures could feasibly have restricted the production of bifaces (Gowlett, Reference Gowlett2011; Key and Lycett, Reference Key and Lycett2017), but habitats suitable for Acheulean-retaining populations were present in Iberia during this period (Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review) and hominin presence is evidenced through flake-and-core technologies and palaeoclimatic data (Martínez and Garriga, Reference Martínez and GArriga2016; Rodríguez et al., Reference Rodríguez, Willmes, Sommer and Mateos2022). Diverse other demographic and climatic factors could also have played a role, with Acheulean cultural information either going regionally extinct or populations dispersing elsewhere (Lycett and von Cramon-Taubadel, Reference Lycett and von Cramon-Taubadel2008; MacDonald et al., Reference MacDonald, Martinon-Torres, Dennell and de Castro2012; Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; French, Reference French2021; Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023).
Of course, substantial date ranges are attached to the sites that bound this early Iberian temporal gap, but these reinforce the inference of Acheulean absence; there is no overlap in their ranges and the resampling technique identified a significant break. An inference of Acheulean absence therefore appears robust given the known archaeological record. It is however unclear what portions of the Iberian Acheulean record remain unknown (Vallverdu et al., Reference Vallverdu, Saladie, Rosas, Huguet, Caceres, Mosquera, Garci-Tabernero, Estalrrich, Lozano-Fernandez, Pineda-Alcala, Carrancho, Villalain, Bourles, Braucher, Lebatard, Vilalta, Esteban-Nadal, Bennasar, Bastir, Lopez-Polin, Olle, Verges, Ros-Montoya, Martinez-Navarro, Grcia, MArtinell, Exposito, Burjachs, Agusti and Carbonell2014; Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023), with happenstance and/or sites excluded from the present analyses potentially playing a role. It is feasible that a site bridging the Cueva Negra to Sima de los Huesos temporal gap may be discovered in the future, even if the high density of Iberian sites currently known at c. 400 ka suggests this to be unlikely (following an assumption of even likelihood of discovery through time).
Northern Europe
Northern Europe is widely thought to have become uninhabitable for hominin populations during the MIS 16, 14, 12 and 10 glacial periods (White and Schreve, Reference White and Schreve2000; Ashton and Lewis, Reference Ashton and Lewis2012; Moncel et al., Reference Moncel, Ashton, Lamotte, Tuffreau, Cliquet and Despriee2015; Ashton and Davis, Reference Ashton and Davis2021). Such has been the strength of narrative concerning an absence of hominins during glacial periods, radiometric and terrace stratigraphy-based dating of artefacts are often applied in combination with, or superseded by, interglacial MIS stage associations (e.g. Bridgland, Reference Bridgland1994; Keen et al., Reference Keen, Hardaker and Lang2006; Parfitt et al., Reference Parfitt, Ashton, Lewis, Abel, Coope, Field, Gale, Hoare, Larkin, Lewis, Karloukovski, Maher, Pegler, Preece, Whittaker and Stringer2010; Antoine et al., Reference Antoine, Moncel, Locht, Limondin-Lozouet, Auguste, Stoetzel, Dabkowski, Voinchet, Bahain and Falgueres2015; Davis et al., Reference Davis, Ashton, Hoare and LEwis2021; Key et al., Reference Key, Lauer, Skinner, Pope, Bridgland, Nobel and Proffitt2022). Paleoclimatic data supports a contraction of the northern range for Acheulean populations in glacial periods (Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review), but a notable north-western persistence was present with areas as far north as southern Britain (~52° latitude) appearing habitable (Rodríguez et al., Reference Rodríguez, Willmes, Sommer and Mateos2022; Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review). The site of Moulin Quignon in northern France further supports an Acheulean presence above 49° latitude during the MIS 16 glacial period (Antoine et al., Reference Antoine, Moncel, Voinchet, Locht, Amselem, Herisson, Hurel and Bahain2019).
There is, therefore, discrepancy between the present results, which suggest Acheulean absence in northern Europe during MIS 14, 12 and 10, along with most of MIS 16 (Table 2), and the aforementioned palaeoclimatic modelling (Rodríguez et al., Reference Rodríguez, Willmes, Sommer and Mateos2022; Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review). Given the present episodes of cultural absence have been identified using sites where some ages were determined via interglacial MIS-stage associations (Supplementary material S1), it is not surprising that significant gaps in the archaeological record were identified during glacial periods. The present northern European results can therefore be interpreted in two ways. If the inference of Acheulean absence in northern Europe during glacial stages is upheld, then the present results provide theoretically grounded empirical support in favour of these absences.
Alternatively, palaeoclimatic data and Moulin-Quignon provide evidence of a need to revise our understanding of Acheulean presence in northern Europe during glacial periods. In short, hominins could have been present in northern Europe during glacial periods – potentially with reduced population levels or as part of seasonal migratory patterns (Figure 3B) (Hosfield, Reference Hosfield2016; Rodríguez et al., Reference Rodríguez, Willmes, Sommer and Mateos2022; Moncel et al., Reference Moncel, Antoine, Herisson, Locht, Hurel and Bahain2022; Leonardi et al., Reference Leonardi, Lycett, Manica and Keyin review) – and Acheulean sites dated using interglacial associations may not necessarily be reliable. In this scenario, the surprise test would only return accurate northern European results after the age of many sites has been re-evaluated. Caution is therefore essential when interpreting the present northern European data. They accurately represent current understanding concerning the presence and absence of Acheulean cultural information in the region, but there is a growing need to reassess the theoretical foundation on which this understanding is based.
Regional Acheulean continuity
Southern Europe, Iberia and northern Europe all identified long periods of temporal cohesion, and therefore Acheulean presence, in the Middle-to-Late Pleistocene. Subsequent to the Acheulean’s emergence in Europe prior to 885 ka (Ollé et al., Reference Ollé, Lombao, Asryan, García-Medrano, Arroyo, Fernandez-Marchena, Yesilova, Caceres, Huguet, Lopez-Polin, Pineda, Garcia-Tabernero, Fidalgo, Rosas, Saladie and Vallverdu2023) it appears to be continually present in southern Europe until soon after the 160.5 ka dated occurrence of Arbo (Spain) (Méndez-Quintas et al., Reference Méndez-Quintas, Demuro, Arnold, Duval, Perez-Gonzalez and Santonja2019; Key et al., Reference Key, Roberts and Jarić2021). The Italian peninsula may have facilitated this continuity by acting as source population for Acheulean information during Iberia’s early period of absence, as evidenced by the sites of Notarchirico and Loreto (Lefèvre et al., Reference Lefèvre, Raynal, Vernet, Kieffer and Piperno2010; Moncel et al., Reference Moncel, Santagata, Pereira, Nomade, Bahain, Voinchet and Piperno2019) (Figure 3). Potentially hinting at a barrier to the flow of Acheulean cultural information to southwest Europe between c. 800 and 500 ka. The re-emergence of bifaces in Iberia after a >400,000-year break could represent an episode of cultural convergence as opposed to a dispersal from elsewhere in southern Europe, but it is potentially more likely that Acheulean populations dispersed from the Italian peninsula given the close proximity (Lycett and von Cramon-Taubadel, Reference Lycett and von Cramon-Taubadel2008; Shipton, Reference Shipton and Groucutt2020; Key, Reference Key2022). North Africa could have also contributed additional cultural information at this point too (Sharon, Reference Sharon2011; Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; Méndez-Quintas et al., Reference Méndez-Quintas, Santonja, Arnold, Cunha-Ribeiro, da Silva, Demuro, Duval, Gomes, Meireles, Monteiri-Rodrigues and Pérez-González2020), leading to an Iberian Acheulean ‘melting pot’. Alternatively, the Acheulean may have been present in Iberia during this c. 300,000-year period, but site discovery rates in Iberia and the Italian Peninsula may vary. Acheulean cultural persistence in southern Europe during MIS 16, 14, 12 and 10 would have allowed the region to act as a source population for northern dispersals (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a; Figure 3).
In Iberia, Acheulean bifaces are constantly present after the tradition’s re-emergence c. 427 ka, while in northern Europe the tradition is present after MIS 10 (c. 337 ka [Lisiecki and Raymo, Reference Lisiecki and Raymo2005]). Suggesting Acheulean hominins to have occupied Iberia through the MIS 10 and 8 glacial periods, and northern Europe through MIS 8. During this period there are notable biological and technological changes in populations, with an early Neanderthal phenotype and the earliest Levallois tools appearing (Arsuaga et al., Reference Arsuaga, Martinez, Arnold, Aranburu, Gracia-Tellez, Sharp, Quam, Falgueres, Pantoja-Perez, Bischoff, Poza-Rey, Pares, Carretero, Demuro, Lorenzo, Sala, Martinon-Torres, Garcia, Alcazar, Cuenca-Bescos, Gomez-Olivencia, Moreno, Pablos, Shen, Rodriguez, Ortega, Garcia, Bonmati, Bermudez de Castro and Carbonell2014; Moncel et al., Reference Moncel, Ashton, Arzarello, Fontana, Lamotte, Scott, Muttillo, Berruti, Nenzioni, Tuffreau and Peretto2020a). Potentially these cultural and biological changes helped populations maintain durable (larger, more genetically diverse) and permanent populations through glacial periods. The precise nature of how, when, where and why the Acheulean in Europe ceased to be a distinct cultural entity free from prepared core technologies remains debated (Moncel et al., Reference Moncel, Moigne and Combier2012; Lycett et al., Reference Lycett, von Cramon-Taubadel and Eren2016; Malinsky-Buller, Reference Malinsky-Buller2016; de Lombera-Hermida et al., Reference de Lombera-Hermida, Rodriguez-Alvarez, Mosquera, Olle, Garcia-Medrano, Pedergnana, Terradillos-Bernal, Lopez-Ortego, Beargallo, Rodriguez-Hidalgo, Saladie, Bermudez de Castro and Carbonell2020; Moncel et al., Reference Moncel, Ashton, Arzarello, Fontana, Lamotte, Scott, Muttillo, Berruti, Nenzioni, Tuffreau and Peretto2020a; Kuhn et al., Reference Kuhn, Moncel, Weinstein-Evron and Zaidner2021), but based on current site definition, considerable overlap is present in Europe between these technologically distinct phenomena (Key et al., Reference Key, Roberts and Jarić2021).
Acheulean persistence in Europe
Both continental-scale scenarios (Europe S1, Europe S2) identified the European Acheulean archaeological record to be temporally cohesive. Once Acheulean cultural information was introduced to the continent, it appears to have never left. This reading of the European record differs from some previous studies, where Acheulean reintroduction events have been proposed from c. 700 to 500 ka (MacDonald et al., Reference MacDonald, Martinon-Torres, Dennell and de Castro2012; Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; Vallverdu et al., Reference Vallverdu, Saladie, Rosas, Huguet, Caceres, Mosquera, Garci-Tabernero, Estalrrich, Lozano-Fernandez, Pineda-Alcala, Carrancho, Villalain, Bourles, Braucher, Lebatard, Vilalta, Esteban-Nadal, Bennasar, Bastir, Lopez-Polin, Olle, Verges, Ros-Montoya, Martinez-Navarro, Grcia, MArtinell, Exposito, Burjachs, Agusti and Carbonell2014; Moncel et al., Reference Moncel, Santagata, Pereira, Nomade, Voinchet, Bahain, Daujeard, Curci, Lemorini, Hardy, Eramo, Berto, Raynal, Arzarello, Mecozzi, Iannucci, Sardella, Allegretta, Delluniversita, Terzano, Dugas, Jouanic, Queffelec, d’Andrea, Valentini, Minucci, Carpentiero and Piperno2020c; French, Reference French2021). Relevant population levels and hominin spatial presence would have varied through this c. 720,000-year period, and potentially we see demographic signals through reduced discovery rates in the archaeological record; for example, between 700 and 800 ka (Figure 3B). Indeed, the temporal gap between Cueva Negra and la Noira returned lower P values compared to most (~0.200 (Supplementary Tables S3 and S4)). The regional-scale losses identified above provide additional evidence of demographic fluctuation within the continent. Diverse factors, including climatic and ecological change, pressures from non-Acheulean hominin populations and disease, could have created pressures leading to lower presence at times (Bar-Yosef and Belfer-Cohen, Reference Bar-Yosef and Belfer-Cohen2001; Lycett and Norton, Reference Lycett and Norton2010; French, Reference French2016). Fundamentally, however, temporal evidence suggests that after c. 885 ka, Acheulean populations were continuously present in Europe and the relevant cultural information was remarkably durable.
Cultural evolutionary mechanisms, including the introduction of new Acheulean information from central Eurasia, and potentially Africa, would have impacted the way the phenomenon was expressed in Europe between 885 and 160 ka (Lycett et al., Reference Lycett, von Cramon-Taubadel and Eren2016; Kuhn, Reference Kuhn2020). This was likely to have been in regionally dependent ways (Sharon and Barsky, Reference Sharon and Barsky2016; Key, Reference Key2019; Shipton and White, Reference Shipton and White2020; García-Medrano et al., Reference García-Medrano, Moncel, Maldonado-Garrido, Ollé and Ashton2023). The fundamental bauplan that allows bifaces at both Barranc de la Boella and Arbo to be defined as part of the same cultural tradition did not, however, change (Lycett and Gowlett, Reference Lycett and Gowlett2008). The present results therefore support the presence of a single, but variable, Acheulean tradition in Europe by supporting the presence of a single branching lineage of Acheulean cultural information (Isaac, Reference Isaac1977; Crompton and Gowlett, Reference Crompton and Gowlett1993; Lycett and Gowlett, Reference Lycett and Gowlett2008; Shipton, Reference Shipton and Groucutt2020; Key, Reference Key2022). At times it may have been regionally absent, but any later reintroduction would have been from another channel in the larger ‘braided stream’ of European Acheulean cultural information (Figure 3).
It is important to re-emphasise that any change to the age of European Acheulean sites, particularly in the north (see above), or the discovery of new sites, could adjust some of the results reported here. Use of the surprise test is also dependent on the assumption that all sites in the main sample represent a continuous lineage of cultural information (i.e. they themselves do not contain an episode of cultural absence). Finally, it is again worth restressing that any Lower Palaeolithic temporal analyses are inevitably limited by the large date ranges associated with most sites (Ollé et al., Reference Ollé, Mosquera, Rodriguez-Alvarez, Garcia-Medrano, Barsky, de Lombera-Hermida and Carbonell2016a; Pope et al., Reference Pope, Blundell, Scott and Cutler2016). Our understanding of prehistoric material culture is, however, bounded by the sum of all information currently known, and as it stands, the most accurate interpretation for European Acheulean continuity is one of cultural persistence that is only rarely punctuated.
Conclusion
Demographic trends in Lower Palaeolithic Europe have recently been argued to represent ‘discontinuous, fragmented European populations who, like those of the Early Pleistocene, visited rather than occupied the continent’ (French, Reference French2021, 128). Other studies have returned similar conclusions for the European Middle Pleistocene based on the fragmented evidence we currently possess (e.g. MacDonald et al., Reference MacDonald, Martinon-Torres, Dennell and de Castro2012; Mosquera et al., Reference Mosquera, Ollé and Rodriguez2013; Moncel et al., Reference Moncel, Santagata, Pereira, Nomade, Voinchet, Bahain, Daujeard, Curci, Lemorini, Hardy, Eramo, Berto, Raynal, Arzarello, Mecozzi, Iannucci, Sardella, Allegretta, Delluniversita, Terzano, Dugas, Jouanic, Queffelec, d’Andrea, Valentini, Minucci, Carpentiero and Piperno2020c; Ashton and Davis, Reference Ashton and Davis2021; Margari et al., Reference Margari, Hodell, Parfitt, Ashton, Grimalt, Kim, Yun, Gibbard, Sytringer, Timmermann and Tzedakis2023). A sparse archaeological record is not, however, tantamount to an absence of hominin populations. What appears to be a substantial temporal gap may in fact represent cultural continuity when contextualised against the rest of the archaeological record.
Here, the temporal cohesion of the European Acheulean archaeological record has been statistically assessed. Five regionally-defined breaks in the Acheulean record were identified: between c. 850 and 500 ka in Iberia, and during MIS 16, 14, 12 and 10 in northern Europe. For each, hominins retaining Acheulean cultural information are inferred to have been absent; be it due to cultural extirpation or populations dispersing to alternative regions. The northern results should, however, be used with caution given increasing evidence that hominins may have been present in northern latitudes during glacial periods. At a continental level, the Acheulean was identified as being continuously present. No breaks were substantial enough, on a relative basis, to infer an absence of hominin populations retaining Acheulean cultural information. The European Acheulean is therefore overwhelmingly characterised as a period of cultural persistence; it was likely a single, braided lineage of cultural information that appears to have always been present in Europe after its first introduction. Regional extinctions occurred and variable technological and morphological trajectories developed, but cultural information would have flowed between populations and dispersal events would have reintroduced the overarching ‘tradition’ (Lycett and Gowlett, Reference Lycett and Gowlett2008) back into unoccupied regions.
Open peer review
To view the open peer review materials for this article, please visit http://doi.org/10.1017/ext.2024.13.
Supplementary material
The supplementary material for this article can be found at http://doi.org/10.1017/ext.2024.13.
Data availability statement
All data and code are included in the Supplementary material.
Acknowledgements
Thanks are extended to Michela Leonardi for assistance during the creation of the site database and to Andreu Ollé for interesting discussions that helped formulate ideas for the present research.
Author contribution
A.K. undertook all aspects of the work.
Financial support
No financial support was received.
Competing interest
The author declares none.